Digital elevation model with the ground-based SAR IBIS-L as basis for volcanic deformation monitoring

Within the project Exupéry, a hybrid deformation observation system is developed, which is part of a Volcano Fast Response System (VFRS). The VFRS shall be a mobile, robust, real-time, easy to use early warning system, which, in case of a volcanic crisis, will support the locale authorities in their decisions about hazard mitigation provisions, especially about the evacuation of people.The hybrid deformation observation system combines a ground-based Synthetic Aperture Radar named IBIS-L with a network of GPS receivers and will allow the continuous, weather independent determination of areal 3D displacements.In this paper we present the results of first tests with IBIS-L in an active quarry near Dieburg, Germany. A Digital Elevation Model (DEM) was determined with IBIS-L and compared with a DEM derived by terrestrial laser scanning and photogrammetry. The ability to determine accurate DEMs with IBIS-L is the basis for the combination of different observation techniques, as without a DEM the displacements observed with IBIS-L cannot be georeferenced.The standard deviation of the differences between the DEM by IBIS-L and the DEM by laser scanning and photogrammetry depends on the surface characteristics. We found a standard deviation of 0.8 m at a slope of rocks and debris, 2.0 m at mining terraces with low vegetation and 3.0 m in a vegetation covered area.     

Slope angle and aspect as influencing factors on the accuracy of the SRTM and the ASTER GDEM databases

Presently, the application of digital elevation or surface models have increasing relevance in all areas of scientific research and in practical engineering applications. The ASTER GDEM and SRTM databases are the most widely used digital surface models, due to their free accessibility and global coverage. The SRTM model was produced using a radar-based technique and the ASTER GDEM was developed using optical stereo image-pairs. Therefore, as all models contain errors (i.e. differences stemming from real surface or vertical biases), errors in these models will also differ. Our aim was to examine these vertical biases and to calculate the rate of error variance. A TIN (Triangulated Irregular Network) model was used as a reference surface, derived from the contour lines of a large scale topographic map. Errors were evaluated with statistical and geoinformation techniques. We discovered significant differences between the surfaces. The mean difference between topographic elevations minus the SRTM-V2 is +2.6 ± 4 m, while the mean difference between topographic elevations minus the SRTM-V3 is +2.7 ± 2.5 m. With the GDEM, the mean difference was 2.7 ± 9.1 m. Furthermore, we found that in the case of SRTM, the differences were significant considering the aspects and the steepness of the slopes: southern and eastern directions and larger slope angles showed greater differences compared to the reference data. The GDEM V2 DEM had a larger error variance, but the error did not vary significantly with slope angle.     

Evaluation of ASTER GDEM and SRTM and their suitability in hydraulic modelling of a glacial lake outburst flood in southeast Tibet

The resolution and accuracy of digital elevation models (DEMs) can affect the hydraulic simulation results for predicting the effects of glacial lake outburst floods (GLOFs). However, for the Tibetan Plateau, high‐quality DEM data are often not available, leaving researchers with near‐global, freely available DEMs, such as the Advanced Spaceborne Thermal Emission and Reflection Radiometer Global Digital Elevation Model (ASTER GDEM) and the Shuttle Radar Topography Mission data (SRTM) for hydraulic modelling. This study explores the suitability of these two freely available DEMs for hydraulic modelling of GLOFs. Our study focused on the flood plain of a potentially dangerous glacial lake in southeastern Tibet, to evaluate the elevation accuracy of ASTER GDEM and SRTM, and their suitability for hydraulic modelling of GLOFs. The elevation accuracies of ASTER GDEM and SRTM were first validated against field global position system (GPS) survey points, and then evaluated with reference to the relatively high precision of 1:50 000 scale DEM (DEM5) constructed from aerial photography. Moreover, the DEM5, ASTER GDEM and SRTM were used as basic topographic data to simulate peak discharge propagation, as well as flood inundation extent and depth in the Hydrologic Engineering Center's River Analysis System one‐dimensional hydraulic model. Results of the three DEM predictions were compared to evaluate the suitability of ASTER GDEM and SRTM for GLOF hydraulic modelling. Comparisons of ASTER GDEM and SRTM each with DEM5 in the flood plain area show root‐mean‐square errors between the former two as ± 15·4 m and between the latter two as ± 13·5 m. Although SRTM overestimates and ASTER GDEM underestimates valley floor elevations, both DEMs can be used to extract the elevations of required geometric data, i.e. stream centre lines, bank lines and cross sections, for flood modelling. However, small errors still exist in the cross sections that may influence the propagation of peak discharge. The flood inundation extent and mean water depths derived from ASTER GDEM predictions are only 2·2% larger and 2·3‐m deeper than that of the DEM5 predictions, whereas the SRTM yields a flood zone extent 6·8% larger than the DEM5 prediction and a mean water depth 2·4‐m shallower than the DEM5 prediction. The modelling shows that, in the absence of high‐precision DEM data, ASTER GDEM or SRTM DEM can be relied on for simulating extreme GLOFs in southeast Tibet. Copyright © 2011 John Wiley & Sons, Ltd.     

Raised Holocene paleo-shorelines along the Capo Vaticano coast (western Calabria,Italy): Evidence of co-seismic and steady-state deformation

Detailed mapping of geomorphological and biological sea-level markers around the Capo Vaticano promontory (western Calabria, Italy), has documented the occurrence of four Holocene paleo-shorelines raised at different altitudes. The uppermost shoreline (PS1) is represented by a deeply eroded fossiliferous beach deposit, reaching an elevation of ∼2.2 m above the present sea-level, and by a notch whose roof is at ∼2.3 m. The subjacent shoreline PS2 is found at an elevation of ∼1.8 m and is represented by a Dendropoma rim, a barnacle band and by a wave-cut platform. Shoreline PS3 includes remnants of vermetid concretions, a barnacle band, a notch and a marine deposit, and reaches an elevation of ∼1.4 m. The lowermost paleo-shoreline (PS4) includes a wave-cut platform and a notch and reaches an elevation of ∼0.8 m. Radiocarbon dating of material from individual paleo-shorelines points to an average uplift rate of 1.2–1.4 mm/yr in the last ∼6 ka at Capo Vaticano. Our data suggest that Holocene uplift was asymmetric, with a greater magnitude in the south-west sector of the promontory, in a manner similar to the long-term deformation attested by Pleistocene terraces. The larger uplift in the south-western sector is possibly related to the additional contribution, onto a large-wavelength regional signal, of co-seismic deformation events, which are not registered to the north-east. We have recognized four co-seismic uplift events at 5.7–5.4 ka, 3.9–3.5 ka, ∼1.9 ka and <1.8 ka ago, superposed on a regional uplift that in the area, is occurring at a rate of ∼1 mm/yr. Our findings places new constrains on the recent activity of border faults south of the peninsula and on the location of the seismogenic source the 1905 destructive earthquake.     

Optical properties and light penetration in a deep,naturally acidic,iron rich lake: Lago Caviahue (Patagonia,Argentina)

The optical properties and light climate in the deep and extremely acid Lake Caviahue have been studied in order to better understand its characteristics and the possible influence upon the phytoplankton community. The absorption coefficients for the dissolved fraction were maximal in the ultraviolet (UV) region and the water absorption spectra showed a shoulder around 300 nm, which was attributed to the concentration of Fe(III). No radiation was detected in the water column below 360 nm. The depth of the 1% incident radiation was dependent of wavelength, showing its maximum of 13.3 m at 565 nm, compared to 1.7 m and 4.8 m at 400 nm and 700 nm, respectively. Phytoplankton biomass was low and showed an almost constant profile with depth despite the relative darkness of the water column. Optical climate of Lake Caviahue is not typical of high elevation lakes but is more similar to low elevation shallow lakes of the Andean region. The chemical composition of the water, mainly Fe oxidation state and concentration, is the responsible for the high attenuation of the UV radiation (UVR). Living organisms are protected of UVR because Lake Caviahue waters are a shield against UV-B.     

Laurentia crustal motion observed using TOPEX/POSEIDON radar altimetry over land

A new method to estimate the vertical crustal motion from satellite altimetry over land was developed. The method was tested around Hudson Bay, where the observed vertical motion is largely caused by the incomplete glacial isostatic adjustment (GIA) as a result of the Laurentide ice sheet deglaciation since the last glacial maximum (LGM). Decadal (1992–2003) TOPEX/POSEIDON radar altimetry data over land surfaces were used. The results presented here are improved compared to a previous study (Lee, H., Shum, C.K., Kuo, C.Y., Yi, Y., Braun, A., 2008. Application of TOPEX altimetry for solid Earth deformation studies. Terr. Atmos. Ocean. Sci. 19, 37–46. doi:10.3319/TAO.2008.19.1-2.37(SA).) which estimated vertical motion only over relatively flat land surfaces (standard deviation of the height variation <40 cm). In this study, we extended the concept of traditional 1-Hz (one-per-frame) radar altimeter ocean stackfiles to build 10-Hz (10-per-frame) land stackfiles over Hudson Bay land regions, and succeeded in obtaining vertical motion estimates over much rougher surfaces (standard deviation of the height variation <2 m). 90-m C-band Shuttle Radar Topography Mission (SRTM) Digital Elevation Model (DEM) is used as a reference surface to select an optimal waveform retracker, to correct surface gradient errors, and to calculate land surface anomalies. Here, we developed an alternative retracker, called the modified threshold retracker, resulting in decadal vertical motion time series over a 1500 km by 1000 km region covering northern Ontario, northeastern Manitoba, and the Great Lakes region which is at the margin of the former Laurentide ice sheet. The average of the estimated uncertainties for the vertical motion is 2.9 mm/year which is comparable to 2.1 mm/year of recent GPS solutions. The estimated vertical motion is compared with other geodetic observations from GPS, tide gauge/altimetry, GRACE, and several GIA models. The data agree best with the laterally varying 3D GIA model, RF3S20 (β = 0.4) whereas the combination of land altimetry solution with other measurements match best with the models RF3S20 (β = 0.0) or RF3S20 (β = 0.2) in terms of mean and standard deviation of the differences. It is anticipated that this innovative technique could potentially be used to provide additional constraints for GIA model improvement, and be applied to other geodynamics studies.     

Improvement of streams hydro‐geomorphological assessment using LiDAR DEMs

Hydro‐geomorphological assessments are an essential component for riverine management plans. They usually require costly and time‐consuming field surveys to characterize the spatial variability of key variables such as flow depth, width, discharge, water surface slope, grain size and unit stream power throughout the river corridor. The objective of this research is to develop automated tools for hydro‐geomorphological assessments using high‐resolution LiDAR digital elevation models (DEMs). More specifically, this paper aims at developing geographic information system (GIS) tools to extract channel slope, width and discharge from 1 m‐resolution LiDAR DEMs to estimate the spatial distribution of unit stream power in two contrasted watersheds in Quebec: a small agricultural stream (Des Fèves River) and a large gravel‐bed river (Matane River). For slope, the centreline extracted from the raw LiDAR DEM was resampled at a coarser resolution using the minimum elevation value. The channel width extraction algorithm progressively increased the centerline from the raw DEM until thresholds of elevation differences and slopes were reached. Based on the comparison with over 4000 differential global positioning system (GPS) measurements of the water surface collected in a 50 km reach of the Matane River, the longitudinal profile and slope estimates extracted from the raw and resampled LiDAR DEMs were in very good agreement with the field measurements (correlation coefficients ranging from 0 · 83 to 0 · 87) and can thus be used to compute stream power. The extracted width also corresponded very well to the channel as seen from ortho‐photos, although the presence of bars in the Matane River increased the level of error in width estimates. The estimated maximum unit stream power spatial patterns corresponded well with field evidence of bank erosion, indicating that LiDAR DEMs can be used with confidence for initial hydro‐geomorphological assessments. Copyright © 2013 John Wiley & Sons, Ltd.     

Extracting drainage networks and their connectivity using LiDAR data

Policies, measures, and models geared towards flood prevention and managing surface waters benefit from high quality data on the presence and characteristics of drainage ditches. As a cost and labour effective alternative for acquiring such data through field surveys, we propose a method (a) to extract vector data representing ditch drainage networks based on local morphologic features derived from high resolution digital elevation models (DEM) and (b) to identify possible connections in the ditch network by calculating a probability of the connectivity using a logistic regression where the predictor variables are characteristics of the ditch centre lines or derived from the DEM. Using Light Detection and Ranging (LiDAR) derived DEMs with a 1 m resolution, the method was developed and tested for a mixed agricultural residential area in north‐eastern Belgium. The derived ditch segments had an error of omission of 8% and an error of commission of 5%. The original positional accuracy of the centre lines of the extracted ditches was 0.6 m and could be improved to 0.4 m by shifting each vertex to the position of the lowest LiDAR point located within a radius equal to the spatial resolution of the used DEM. About 69% of the false disconnections in the network were identified and corrected leading to a reduction of the unconnected parts of the ditch network by 71%. The extracted and connected network approximated the reference ditch network fairly well.     

Timing for high lake levels of Qinghai Lake in the Qinghai-Tibetan Plateau since the Last Interglaciation based on quartz OSL dating

Qinghai Lake is situated in the northeast of the Qinghai-Tibetan Plateau (QTP). Its size and proximity to the junction of three major climate systems make it sensitive to climate changes. Some investigations on shorelines of Qinghai Lake suggested highstands during MIS 3, but to what extent the lake level was higher than today is yet undetermined. Others proposed that the maximum highstands probably dated to MIS 5. It has also been shown that the lake level 120 m higher than today occurred at around 12 ka. Most of these previous ages were obtained using ¹⁴C dating or multiple-aliquot IRSL/OSL dating. For ¹⁴C dating, because of the dating limit (<40 ka) and the lack of suitable dating materials in this arid area, it is difficult to establish reliable chronological control. In the present study, seven samples collected from lacustrine deposits (five samples) and sand wedges (two samples) were dated using quartz optically stimulated luminescence (OSL) with the single aliquot regenerative-dose (SAR) protocol. OSL dating results showed that (1) the lake had experienced two high lake levels, one was in MIS 5 and another in early to middle MIS 3; (2) no evidence of high lake levels in MIS 4 has been found; (3) the alluvial gravels, whose surface is at an elevation of ～3246 m, were formed at least 28.8 ± 2.3 ka ago, and the widespread sand wedges within the alluvial gravels were formed during the period of 15.1–28.8 ka, which implied that the lake level had not reached an elevation of ～3240 m after 28.8 ± 2.3 ka.     

Active tectonics of the Northern Rif (Morocco) from geomorphic and geochronological data

We present results of a geomorphological and morphotectonic analysis of the northeastern part of the Rif. We show that the present day kinematics of the Rif is characterized by active deformation along the Trougout and Nekor faults in the North-East. Digital Elevation Models of offset drainage features (streams, fluvial terraces) allow determining a normal-left-lateral motion along the Trougout fault and a left-lateral strike-slip motion along the Nekor fault. Preliminary ³He cosmogenic dates of tectonic markers yield vertical and horizontal slip rates of ∼0.9 mm/yr and ∼0.5 mm/yr, respectively along the Trougout fault. The present-day localized transtension seen in the north-eastern Rif morphology (Ras Tarf) is coeval with uplifted marine terraces near the Al Hoceima Bay. U/Th dating of shells yield an average uplift rate of ∼0.2 mm/yr during the past 500 ka. These data show that active transtension in the northeastern Rif is also associated with uplift. These new morphotectonic constraints are consistent with the GPS measurements showing southwestward overall motion of most of the Rif belt with respect to stable Africa.     

TEC variations over the Mediterranean before and during the strong earthquake (M = 6.5) of 12th October 2013 in Crete,Greece

In this paper, the total electron content (TEC) data from eight global positioning system (GPS) stations of the EUREF network, provided by IONOLAB (Turkey), were analyzed using discrete Fourier analysis to investigate the TEC variations over the Mediterranean before and during the strong earthquake of 12th October 2013, which occurred west of Crete, Greece. In accordance with the results of similar analyses in the area, the main conclusions of this study are the following: (a) TEC oscillations in a broad range of frequencies occur randomly over an area of several hundred km from the earthquake and (b) high frequency oscillations (f ⩾ 0.0003 Hz, periods T ⩽ 60 m) may point to the location of the earthquake with questionable accuracy. The fractal characteristics of the frequency distribution may point to the locus of the earthquake with higher accuracy. We conclude that the lithosphere–atmosphere–ionosphere coupling (LAIC) mechanism through acoustic or gravity waves could explain this phenomenology.     

Effects of hypoxia on benthic organisms in Tokyo Bay,Japan: A review

Bottom hypoxia (dissolved oxygen concentration ⩽2 ml l⁻¹) from anthropogenic eutrophication is a growing global concern. Here, we summarized characteristics of hypoxia and its effects on benthic organisms in Tokyo Bay. Despite recent decreases in nutrient inputs, hypoxia has been increasing in duration and spatial extent, suggesting that the substantial loss of tidal flats from reclamation is contributing to a decrease in the ability of Tokyo Bay to recycle nutrients. Hypoxia develops in the central to northern part of the bay and persists from spring to autumn, causing defaunation of benthic organisms. After the abatement of hypoxia in autumn, the defaunated area is recolonized, either through migration or larval settlement. Some megabenthic species with a spawning peak in spring and summer experience failure of larval settlement, which is probably due to hypoxia. The adverse effects of hypoxia are an impediment to recovery of benthic organisms in Tokyo Bay.     

Spatial–temporal changes in Andean plateau climate and elevation from stable isotopes of mammal teeth

Paleoelevation constraints from fossil leaf physiognomy and stable isotopes of sedimentary carbonate suggest that significant surface uplift of the northern Andean plateau, on the order of 2.5 ± 1 km, occurred between ～ 10.3 and 6.4 Ma. Independent spatial and temporal constraints on paleoelevation and paleoclimate of both the northern and southern plateau are important for understanding the distribution of rapid surface uplift and its relation to climate evolution across the plateau. This study focuses on teeth from modern and extinct mammal taxa (including notoungulates, pyrotheres, and litopterns) spanning ～ 29 Ma to present, collected from the Altiplano and Eastern Cordillera of Bolivia (16.2°S to 21.4°S), and lowland Brazil. Tooth enamel of large, water-dependent mammals preserves a record of surface water isotopes and the type of plants that animals ingested while their teeth were mineralizing. Previous studies have shown that the δ¹⁸O of modern precipitation and surface waters decrease systematically with increasing elevations across the central Andes. Our results from high elevation sites between 3600 and 4100 m show substantially more positive δ¹⁸O values for late Oligocene tooth samples compared to < 10 Ma tooth δ¹⁸O values. Late Oligocene teeth collected from low elevation sites in southeast Brazil show δ¹⁸O values similar (within 2‰) to contemporaneous teeth collected at high elevation in the Eastern Cordillera. This affirms that the Andean plateau was at a very low elevation during the late Oligocene. Late Oligocene teeth from the northern Eastern Cordillera also yield consistent δ¹³C values of about ? 9‰, indicating that the environment was semi-arid at that time. Latitudinal gradients in δ¹⁸O values of late Miocene to Pliocene fossil teeth are similar to modern values for large mammals, suggesting that by ～ 8 Ma in the northern Altiplano and by ～ 3.6 Ma in the southern Altiplano, both regions had reached high elevation and established a latitudinal rainfall gradient similar to modern.     

Evaluation of on-line DEMs for flood inundation modeling

Recent and highly accurate topographic data should be used for flood inundation modeling, but this is not always feasible given time and budget constraints so the utility of several on-line digital elevation models (DEMs) is examined with a set of steady and unsteady test problems. DEMs are used to parameterize a 2D hydrodynamic flood simulation algorithm and predictions are compared with published flood maps and observed flood conditions. DEMs based on airborne light detection and ranging (LiDAR) are preferred because of horizontal resolution, vertical accuracy (∼0.1 m) and the ability to separate bare-earth from built structures and vegetation. DEMs based on airborne interferometric synthetic aperture radar (IfSAR) have good horizontal resolution but gridded elevations reflect built structures and vegetation and therefore further processing may be required to permit flood modeling. IfSAR and shuttle radar topography mission (SRTM) DEMs suffer from radar speckle, or noise, so flood plains may appear with non-physical relief and predicted flood zones may include non-physical pools. DEMs based on national elevation data (NED) are remarkably smooth in comparison to IfSAR and SRTM but using NED, flood predictions overestimate flood extent in comparison to all other DEMs including LiDAR, the most accurate. This study highlights utility in SRTM as a global source of terrain data for flood modeling.     

New insights into the origin and evolution of Lake Vida,McMurdo Dry Valleys,Antarctica — A noble gas study in ice and brines

Unlike other lakes in the McMurdo Dry Valleys, Antarctica, Lake Vida has a thick (~ 19 m) ice cover sealing a liquid brine body of unusually high salinity (~ 245 g/L) from the atmosphere. To constrain the conditions under which the atypical Lake Vida ice cover formed and evolved, 19 ice samples were collected down to a depth of ~ 14 m, together with three brine samples trapped in the ice at ~ 16 m for analysis of helium, neon, argon, krypton, and xenon concentrations. The broad pattern of noble gas concentrations for Lake Vida samples is fundamentally different from that of air saturated water (ASW) at 0 °C and an elevation of 340 m for salinities of 0 (ice) and 245 g/L (brine). Overall, ice samples are enriched in He and depleted in Ne with saturation relative to ASW averages of 1.38 and 0.82, respectively, and strongly depleted in Ar, Kr, and Xe with relative saturations of 0.10, 0.06, and 0.05, respectively. By contrast, brine samples are generally depleted in He and Ne (relative saturation averages of 0.33 and 0.27, respectively) but enriched in Ar, Kr, and Xe, with relative saturation averages of 1.45, 3.15, and 8.86, respectively. A three-phase freezing partitioning model generating brine, ice and bubble concentrations for all stable noble gases was tested and compared with our data. Measured brine values are best reproduced for a salinity value of 175 g/L, a pressure of 1.1 atm, and a bubble volume of 20 cm³ kg^?1. Sensitivity tests for ice + bubble samples show an ideal fit for bubble volumes of ~ 1–2 cm³ kg^?1. Our results show that the conditions under which ice and brine formed and evolved at Lake Vida are significantly different from other ice-covered lakes in the area. Our brine data suggest that Lake Vida may be transitioning from a wet to a dry-based lake, while the ice + bubble data suggest at least partial re-equilibration of residual liquid with the atmosphere as ice forms at the top of Lake Vida ice cover.     

APPLICATION OF DEM GENERATION TECHNOLOGY FROM HIGH RESOLUTION SATELLITE IMAGE IN QUANTITATIVE ACTIVE TECTONICS STUDY: A CASE STUDY OF FAULT SCARPS IN THE SOUTHERN MARGIN OF KUMISHI BASIN

Traditional method to generate Digital Elevation Model (DEM)through topographic map and topographic measurement has weak points such as low efficiency, long operating time and small range. The emergence of DEM-generation technology from high resolution satellite image provides a new method for rapid acquisition of large terrain and geomorphic data, which greatly improves the efficiency of data acquisition. This method costs lower compared with LiDAR (Light Detection and Ranging), has large coverage compared with SfM (Structure from Motion). However, there is still lack of report on whether the accuracy of DEM generated from stereo-imagery satisfies the quantitative research of active tectonics. This research is based on LPS (Leica Photogrammetry Suit)software platform, using Worldview-2 panchromatic stereo-imagery as data source, selecting Kumishi Basin in eastern Tianshan Mountains with little vegetation as study area. We generated 0.5m resolution DEM of 5-km swath along the newly discovered rupture zone at the south of Kumishi Basin, measured the height of fault scarps on different levels of alluvial fans based on the DEM, then compared with the scarp height measured by differential GPS survey in the field to analyze the accuracy of the extracted DEM. The results show that the elevation difference between the topographic profiles derived from the extracted DEM and surveyed by differential GPS ranges from -2.82 to 4.87m. The shape of the fault scarp can be finely depicted and the deviation is 0.30m after elevation correction. The accuracy of measuring the height of fault scarps can reach 0.22m, which meets the need of high-precision quantitative research of active tectonics. It provides great convenience for rapidly obtaining fine geometry, profiles morphology, vertical dislocations of fault and important reference for sites selection for trench excavation, slip rate, and samples. This method has broad prospects in the study of active tectonics.     

The equivalent dose of different grain size quartz fractions from lakeshore sediments in the arid region of north China

Equivalent dose (D_e) values were measured by using medium aliquots of different grain size quartz fractions of five lakeshore sediments from the arid region of north China. There are two different relationships between D_e values and grain sizes of these five samples. The first relationship is that the D_e values obtained from various grain sizes are in agreement within 1 delta errors. The second relationship is that D_e values are similar to each other for fractions between 125 and 300 μm, while the D_e value of the 63–90 μm fraction is 40～55% smaller than others. For example, the D_e values obtained for sample #3 are 20.15 ± 1.19 Gy, 19.80 ± 0.83 Gy and 20.93 ± 1.06 Gy for fractions of 90–125, 125–150 and 250–300 μm respectively, but are 10.79 ± 0.84 Gy for the 63–90 μm fraction. The second relationship can't be interpreted by previous studies of both dosimetry and heterogeneous bleaching. It is deduced for sample #2, #3 and #6 that fine particles (<90 μm) intruded after the dominant sedimentation. Comparison of OSL ages from different grain size fractions of sample #2 with a radiocarbon age from the same lithologic layer supports that fractions coarser than 125 μm yield more reliable burial ages, while the fraction finer than 90 μm yields underestimated ages for some lakeshore sediments from this arid region.     

Volume and mass changes of the Greenland ice sheet inferred from ICESat and GRACE

This study examines the recent evolution of the Greenland ice sheet and its six major drainage basins. Based on laser altimetry data acquired by the Ice, Cloud and Land Elevation Satellite (ICESat), covering the period September–November 2003 to February–March 2008, ice surface height changes and their temporal variations were inferred. Our refined repeat track analysis is solely based on ICESat data and is independent of external elevation models, since it accounts for both ice height changes and the local topography. From the high resolution ice height change pattern we infer an overall mean surface height trend of −0.12 ± 0.006 m yr⁻¹. Furthermore, the largest changes could be identified at coastal margins of the ice sheet, exhibiting rates of more than −2 m yr⁻¹. The total ice volume change of the entire ice sheet amounts to −205.4 ± 10.6 km³ yr⁻¹. In addition, we assessed mass changes from 78 monthly Gravity Recovery and Climate Experiment (GRACE) solutions. The Release-04 gravity field solutions of GeoForschungsZentrum Potsdam cover the period between August 2002 and June 2009. We applied an adjusted regional integration approach in order to minimize the leakage effects. Attention was paid to an optimized filtering which reduces error effects from different sources. The overall error assessment accounts for GRACE errors as well as for errors due to imperfect model reductions. In particular, errors caused by uncertainties in the glacial isostatic adjustment models could be identified as the largest source of errors. Finally, we determined both seasonal and long-term mass change rates. The latter amounts to an overall ice mass change of −191.2 ± 20.9 Gt yr⁻¹ corresponding to 0.53 ± 0.06 mm yr⁻¹ equivalent eustatic sea level rise. From the combination of the volume and mass change estimates we determined a mean density of the lost mass to be 930 ± 11 kg m⁻³. This value supports our applied density assumption 900 ± 30 kg m⁻³ which was used to perform the volume–mass-conversion of our ICESat results. Hence, mass change estimates from two independent observation techniques were inferred and are generally in good agreement.     

Comparison of OSL ages from young dune sediments with a high-resolution independent age model

In order to test the accuracy of quartz optically stimulated luminescence (OSL) dating for young dune sediments (<100 a), a series of aerial images of a migrating sand dune is used to cross validate OSL ages. The investigated dune is situated on the northern part of the island of Sylt (southern North Sea). Based on aerial images and a map from 1925 to 2009 and the internal architecture of the dune obtained by ground-penetrating radar (GPR), an independent age model has been developed to attribute sedimentary-architectural elements of the dune to time. The annual rate of dune migration is calculated to be around 4.1 m/a. Along a 245 m transect oriented parallel to the direction of dune movement, 13 samples for OSL dating were collected at equidistant locations. Sand-sized quartz (150–250 μm) was used for determining the equivalent dose (D_e) applying a single-aliquot regenerative-dose (SAR) protocol. Results show that the oldest OSL age from the investigated recent dune appeared to be 110 ± 10 a, whereas the modern analogue was dated to 34 ± 3 a. In comparison with the aerial images, the OSL ages show a systematic overestimation of 10–40 a for six out of seven younger samples, which are expected to be younger than ∼60 a. This offset is negligible for older samples, but a substantial error in these younger ages. The overestimation is originated from a combination of small thermal transfer of 4–12 mGy during preheat and incomplete bleaching in medium OSL component causing a residual dose of about 15 mGy. The contribution of the incompletely bleached medium component cannot be removed totally by an early background subtraction approach. Despite the observed offset for youngest samples, this study corroborates the suitability of the OSL technique to date young dune sediments (<100 a).     

Helicopter-borne measurements of sea ice thickness,using a small and lightweight,digital EM system

Sea ice is an important climate variable and is also an obstacle for marine operations in polar regions. We have developed a small and lightweight, digitally operated frequency-domain electromagnetic-induction (EM) system, a so-called EM bird, dedicated for measurements of sea ice thickness. It is 3.5 m long and weighs only 105 kg, and can therefore easily be shipped to remote places and operated from icebreakers and small helicopters. Here, we describe the technical design of the bird operating at two frequencies of f1 = 3.68 kHz and f2 = 112 kHz, and study its technical performance. On average, noise amounts to ± 8.5 ppm and ± 17.5 ppm for f1 and f2, respectively. Electrical drift amounts to 200 ppm/h and 2000 ppm/h for f1 and f2, during the first 0.5 h of operation. It is reduced by 75% after 2 h. Calibration of the Inphase and Quadrature ppm signals varies by 2 to 3%. A sensitivity study shows that all these signal variations do affect the accuracy of the ice thickness retrieval, but that it remains better than ± 0.1 m over level ice in most cases. This accuracy is also confirmed by means of comparisons of the helicopter EM data with other thickness measurements. The paper also presents the ice thickness retrieval from single-component Inphase data of f1.