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.     

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.     

Water mass variation in the Mediterranean and Black Seas

The mass-induced sea level variability and the net mass transport between Mediterranean Sea and Black Sea are derived for the interval between August 2002 and July 2008 from satellite-based observations and from model data. We construct in each basin two time series representing the basin mean mass signal in terms of equivalent water height. The first series is obtained from steric-corrected altimetry while the other is deduced from GRACE data corrected for the contamination by continental hydrology. The series show a good agreement in terms of annual and inter-annual signals, which is in line with earlier works, although different model corrections influence the consistency in terms of seasonal signal and trend.In the Mediterranean Sea, we obtain the best agreement using a steric correction from the regional oceanographic model MFSTEP and a continental hydrological leakage correction derived from the global continental hydrological model WaterGAP2. The inter-annual time series show a correlation of 0.85 and a root mean square (RMS) difference of 15 mm. The two estimates have similar accuracy and their annual amplitude and phase agree within 3 mm and 23 days respectively. The GRACE-derived mass-induced sea level variability yields an annual amplitude of 27 ± 5 mm peaking in December and a trend of 5.3 ± 1.9 mm/yr, which deviates within 3 mm/yr from the altimetry-derived estimate.In the Black Sea, the series are less consistent, with lower accuracy of the GRACE-derived estimate, but still show a promising agreement considering the smaller size of the basin. The best agreement is realized choosing the corrections from WaterGAP2 and from the regional oceanographic model NEMO. The inter-annual time series have a correlation and RMS differences of 0.68 and 55 mm, their annual amplitude and phase agree within 4 mm and 6 days respectively. The GRACE-derived seawater mass signal has an annual amplitude of 32 ± 4 mm peaking in April. On inter-annual time scales, the mass-induced sea level variability is stronger than in the Mediterranean Sea, with an increase from 2003 to 2005 followed by a decrease from 2006 to 2008.Based on mass conservation, the mass-induced sea level variations, river runoff and precipitation minus evaporation are combined to derive the strait flows between the basins and with the Atlantic Ocean. At the Gibraltar strait, the net inflow varies annually with an amplitude of 52 ± 10 × 10⁻³ Sv peaking end of September (1 Sv = 10⁶ m³ s⁻¹). The inflow through the Bosphorus strait displays an annual amplitude of 13 ± 3 ×10⁻³ Sv peaking in the middle of March. Additionally, an increase of the Gibraltar net inflow (3.4 ± 0.8 × 10⁻³ Sv/yr) is detected.     

Characteristics of narrow bipolar pulses observed in Malaysia

Narrow bipolar pulses (NBPs) are considered as isolated intracloud events with higher peak amplitude and strong high frequency emission compared to the first return strokes and other intracloud discharges. From 182 NBPs recorded in Malaysia in the tropic, 75 were narrow negative bipolar pulses (NNBPs) while 107 were narrow positive bipolar pulses (NPBPs). The mean duration of NNBPs was 24.6±17.1 μs, while 30.2±12.3 μs was observed for NPBPs. The mean full-width at half-maximum (FWHM) was 2.2±0.7 and 2.4±1.4 μs for NNBPs and NPBPs, respectively. The mean peak amplitude of NPBPs normalized to 100 km was 22.7 V/m, a factor of 1.3 higher than that of NNBPs which is 17.6 V/m. In contrast to the previous studies, it was observed that the electric field change was characterized by a bipolar pulse with a significant amount of fine structures separated by a few tens of nanoseconds intervals, embedded on it.     

Evaluation of DEM generation based on Interferometric SAR using TanDEM-X data in Tokyo

This study is focused on the evaluation of a Digital Elevation Model (DEM) for Tokyo, Japan from data collected by the recently launched TerraSAR add-on for Digital Elevation Measurements (TanDEM-X), satellite of the German Aerospace Center (DLR). The aim of the TanDEM-X mission is to use Interferometric SAR techniques to generate a consistent high resolution global DEM dataset. In order to generate an accurate global DEM using TanDEM-X data, it is important to evaluate the accuracy at different sites around the world. Here, we report our efforts to generate a high-resolution DEM of the Tokyo metropolitan region using TanDEM-X data. We also compare the TanDEM-X DEM with other existing DEMs for the Tokyo region. Statistical techniques were used to calculate the elevation differences between the TanDEM-X DEM and the reference data. Two high-resolution LiDAR DEMs are used as independent reference data. The vertical accuracy of the TanDEM-X DEM evaluated using the Root Mean Square Error (RMSE) is considerably higher than the existing global digital elevation models. However, the local area DEM generated by Geospatial Information Authority of Japan (GSI DEM) showed the highest accuracy among all non-LiDAR DEM’s. The vertical accuracy in terms of RMSE estimated using the 2 m LiDAR as reference is 3.20 m for TanDEM-X, 2.44 m for the GSI, 7.00 m for SRTM DEM and 10.24 m for ASTER-GDEM. We also compared the accuracy of TanDEM-X with the other DEMs for different types of land cover classes. The results show that the absolute elevation error of TanDEM-X is higher for urban and vegetated areas, likewise to those observed for other global DEM’s. This is probably because the radar signals used by TanDEM-X tend to measure the first reflective surface that is encountered, which is often the top of the buildings or canopy. Hence, the TanDEM-X based DEM is more akin to a Digital Surface Model (DSM).     

40Ar / 39Ar geochronology constraints on late miocene weathering rates in Minas Gerais,Brazil

⁴⁰Ar / ³⁹Ar incremental heating ages for twenty one grains of cryptomelane, collected at 0, 42, 45, and 60 m depths in the Cachoeira Mine weathering profile, Minas Gerais, permit calculating long-term (10 Ma time scale) weathering rate (saprolitization rate) in SE Brazil. Pure well-crystallized cryptomelane grains with high K contents (3–5 wt.%) yield reliable geochronological results. The ⁴⁰Ar / ³⁹Ar plateau ages obtained decrease from the top to the bottom of the profile (12.7 ± 0.1 to 7.6 ± 0.1 Ma at surface; 7.6  ± 0.2 to 6.1 ± 0.2 Ma at 42 m; and 7.1 ± 0.2 to 5.9 ± 0.1 Ma at 45 m; 6.6 ± 0.1 to 5.2 ± 0.1 Ma at 60 m), yielding a weathering front propagation rate of 8.9 ± 1.1 m/m.y. From the geochronological results and the mineral transformations implicit by the current mineralogy in the weathering profiles, it is possible to calculate the saprolitization rate for the Cachoeira Mine lithologies and for adjacent weathering profiles developed on granodiorites and schists. The measured weathering front propagation rate yields a saprolitization rate of 24.9 ± 3.1 t/km²/yr. This average long-term (> 10 Ma) saprolitization rate is consistent with mass balance calculations results for present saprolitization rates in weathering watersheds. These results are also consistent with long-term saprolitization rates estimated by combining cosmogenic isotope denudation rates with mass balance calculations.     

Miocene magmatism and tectonics within the Peri-Alboran orogen (western Mediterranean)

The aim of this paper concerns Miocene igneous activity in the Alboran Sea and Peri-Alboran area (northern Morocco, western Algeria and Betic Cordilleras in Spain), considering its age and its location with regard to major tectonics structures.We have compiled previous K-Ar isotopic ages of lavas and plutonic boulders and intrusives with an error of ±1σ and completed this set by a new K-Ar isotopic age for andesitic tuffites from Alboran Island. Geochemistry of most of these samples has been considered after previous analyses completed with new data for Spain magmatism. These two sets of data allow us to place the magmatic activity within the regional stratigraphy and tectonics and their chronological framework of the three major tectonic phases of the Maghrebian orogen, at 17 Ma (Burdigalian), 15 Ma (Langhian) and 9 Ma (Tortonian). Petro-geochemical characteristics are compared through time and geographical locations. A major goal of this coupled approach is to help the elaboration of possible geodynamical processes.As an application, we present the case study of the Dellys, Djinet and Thenia region (east of Algiers) where the successive magmatic events between 19.4 ± 1 and 11.6 ± 0.5 Ma are closely related to the local tectonics and sedimentation.The Peri-Alboran igneous activity is placed in a multidisciplinary framework. Timing of activity is defined according to the ages of the neighbouring sedimentary units and the K-Ar ages of igneous rocks.In Spain, the Cabo de Gata-Carboneras magmatic province displays late Oligocene and early Miocene leucogranitic dikes, dated from 24.8 ± 1.3 to 18.1 ± 1.2 Ma; three following andesitic to rhyolitic events took place around 15.1 ± 0.8 to 14.0 ± 0.7 Ma, 11.8 ± 0.6 to 9.4 ± 0.4 Ma, 8.8 ± 0.4 to 7.9 ± 0.4 Ma; this last event displays also granitic rocks. Lamproitic magmas dated between 8.4 ± 0.4 and 6.76 ± 0.04 Ma were emplaced after the Tortonian phase.In Morocco, after the complex building of the Ras Tarf volcanic edifice, major calc-alkaline to shoshonitic volcanoes were built between 9.0 ± 0.5 and 4.8 ± 0.5 Ma, in particular the large Gourougou volcanic complex. Near Oujda, volcanic activity of alkaline affinity leads to multiple emissions of basalts throughout Pliocene times until the beginning of Pleistocene, between 6.2 ± 0.3 and 1.5 ± 0.1 Ma.In the Alboran domain, an age of 19.7 ± 0.8 Ma is reported (this study) for the andesitic tuffites that form the emergent part of the Alboran Island. This age is comparable to that of the Algerian tuffites and cherts “silexites” and the Malaga ones in Spain. Younger activity, completely separated from the previous one, forms the low-K basaltic andesitic dikes from Alboran Island, dated between 9.1 ± 0.5 and 7.5 ± 0.3 Ma. Along the Alboran Ridge both low-K and high-K andesites to dacites were emitted in the estimated range of 10.7–8.7 Ma. Low-K and high-K andesites to dacites sampled at ODP sites 977 and 978 into the East Alboran Basin, are dated between 12.1 ± 0.2 and 9.3 ± 0.1 Ma.We propose to relate with the Trans-Alboran lineament only the post-Tortonian igneous activity younger than 9 Ma.     

Thermal history and origin of the Tanzanian Craton from Pb isotope thermochronology of feldspars from lower crustal xenoliths

Common and radiogenic Pb isotopic compositions of plagioclase and anti-perthitic feldspars from granulite-facies lower crustal xenoliths from the Labait Volcano on the eastern margin of the Tanzanian Craton have been measured via laser ablation MC-ICP-MS. Common Pb in plagioclase and a single stage Pb evolution model indicate that the lower crust of the Tanzanian Craton was extracted from mantle having a ²³⁸U/²⁰⁴Pb of 8.1 ± 0.3 and a ²³²Th/²³⁸U of 4.3 ± 0.1 at 2.71 ± 0.09 Ga (all uncertainties are 2σ). Since 2.4 Ga, some orthoclase domains within anti-perthites have evolved with a maximum ²³⁸U/²⁰⁴Pb of 6 and ²³²Th/²³⁸U of 4.3. The spread in Pb isotopic composition in the anti-perthitic feldspars yields single crystal Pb–Pb isochrons of ～ 2.4 Ga, within uncertainty of U–Pb zircon ages from the same sample suite. The Pb isotopic heterogeneities imply that these granulites resided at temperatures < 600 °C in the lower crust of the Tanzanian Craton from ca. 2.4 Ga to the present. In concert with the chemistry of surface samples, mantle xenoliths, and lower crustal xenoliths, our data imply that the cratonic lithosphere in Tanzania formed ca. ～ 2.7 Ga, in a convergent margin setting, and has remained undisturbed since 2.7 Ga.     

Geochemical and geo-statistical assessment of selected heavy metals in the surface sediments of the Gorgan Bay,Iran

We investigated heavy metal concentrations of zinc (Zn), copper (Cu), chromium (Cr), and lead (Pb), their spatial distribution and enrichment factor index in surface sediments of the Gorgan Bay. Sediment Quality Guidelines were also applied to assess adverse biological effects of these metals. Heavy metals were determined by inductively coupled plasma-mass spectroscopy (ICP-MS). The results indicated mean concentrations (ppm) of heavy metals were (mean ± S.D.) Pb: 11.5 ± 4.88, Cu: 18 ± 8.83, Zn: 42 ± 22.15 and Cr: 32 ± 15.19. Based on Enrichment index, the Gorgan Bay is a low-enriched to non-enriched bay. Heavy metal contents were lower than the standard limits of PEL, ERL, and ERM that reveal no threatening influence of the metals in the Bay.     

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.     

The penultimate deglaciation: High-resolution paleoceanographic evidence from a north–south transect along the eastern Nordic Seas

The penultimate termination has been studied with focus on oceanographic changes in the eastern Nordic Seas and the influence of these changes on the surrounding ice sheets and vice versa. Repeatedly, major changes in the strength of the Atlantic Meridional Overturning Circulation (AMOC) occurred during the studied interval. Times of strong overturning and increased heat transport northwards were of importance in triggering one minor and two major disintegration events. The two major disintegration events were separated by a deglacial pause, characterized by a strong AMOC. The same variability is seen throughout the eastern Nordic Seas, from the Faeroe–Shetland Channel in the south to the Fram Strait in the north. Some of the oceanographic changes occurring during the penultimate termination are comparable with changes seen through the last deglaciation. Reduced winter moisture flux and increased summer melting due to the present insolation forcing further amplified the rate of ice sheet disintegration. Calculated sea-level change through TII shows a mean change of 121 ± 4 m, 41 ± 16 m in the first step and 80 ± 13 m in the last step.     

Uncertainty in seasonal snow reconstruction: Relative impacts of model forcing and image availability

There are many areas of uncertainty when solving the inverse problems of snow water equivalent (SWE) reconstruction. These include (i) the ability to infer the Final Date of the Seasonal Snow (FDSS) cover, particularly from remote sensing; (ii) errors in model forcing data (such as air temperature or radiation fluxes); and (iii) weaknesses in the snow model used for the reconstruction, associated with both the fidelity of the equations used to simulate snow processes (structural uncertainty) and the parameter values selected for use in the model equations. We investigate the trade-offs among these sources of uncertainty using 10,000 station-years worth of data from the western US SNOTEL network. Model structural and parameter uncertainty are eliminated by using a perfect model scenario i.e. comparing results to modelled control runs. The model was calibrated for each station-year to ensure that the model simulations reflect reality. Results indicate that for a temperature index model, a ±5 days error in FDSS gives a median −25%/+32% error in maximum SWE. A 1 °C air temperature bias produces a SWE error larger than a 5 days error in the FDSS for 50% of the 10,000 cases. Similarly, a 5 days error in FDSS could be accounted for by a net radiation error of 13 W m⁻² or less during the melt period, in 50% of cases. Mean absolute errors of 1 °C or more are typically reported in the literature for air temperature interpolations at high elevations. Observed solar radiation during the melt season can differ by 30 W m⁻² over relatively short distances, while estimates from reanalysis (NARR, ERA-Interim, MERRA, CFSRR) and GOES satellites typically span more than 40 W m⁻². Using data from both MODIS sensors (Terra & Aqua) at all snow covered points in the western US, a consecutive 5 days gap in imagery at time of FDSS is likely to occur only 5–10% of the time. This work shows that errors in model forcing data are at least as important, if not more, than image availability when reconstructing SWE.     

Vertical distribution of HTO and 125I− transport parameters in Boom Clay in the Essen-1 borehole (Mol,Belgium)

In the frame of the R&D activities performed on the Boom Clay for assessing the suitability of deep clayey formations for radioactive waste disposal, the transferability of the scientific results obtained on the Boom Clay in Mol to the whole Campine Basin is investigated. Boreholes were drilled at different locations (e.g. Mol, Doel, Essen) and cores were sampled over the entire thickness of the Boom Clay formation on which the migration parameters for iodide and tritiated water (HTO) are determined.At Essen, the transport parameters in the Boom Clay can be considered as homogeneous in the range from 159 m to 241 m Below Drilling Table. The average hydraulic conductivity is (5.4 ± 1.7) × 10⁻¹² m/s. The average ηR value for iodide is 0.25 ± 0.03 and 0.42 ± 0.05 for HTO. For HTO, this high value is mainly due to a higher value in the Putte Member (0.46 ± 0.03) compared to the other members (0.39 ± 0.02). The apparent diffusion coefficient is (1.3 ± 0.1) × 10⁻¹⁰ m²/s for HTO and (1.1 ± 0.2) × 10⁻¹⁰ m²/s for iodide. The expected effect of ionic strength (increasing with depth) on the ηR value of iodide is of the same size as the measurement error, what might explain why it was not observed.On a lateral (horizontal) level, the hydraulic conductivity at the Essen-1 borehole (5.4 × 10⁻¹² m/s) lies between that of Boom Clay in Mol-1 (2.5 × 10⁻¹² m/s) and that of Boom Clay in Doel-2b (1.4 × 10⁻¹¹ m/s). For iodide, the higher ηR value in Essen-1 and Doel-2b (ηR ≈ 0.25) than in Mol-1 (ηR ≈ 0.16) can partly be explained by a higher ionic strength of the pore water. Apart from the Putte Member at Essen-1, the HTO porosities of the Terhagen Member and the Transition zone in Essen are in the range of the average porosities observed in Mol and Doel (ηR ≈ 0.37–0.39). For both iodide and HTO, the value of the apparent diffusion coefficient D_app is similar in Mol-1 and in Doel-2b, with a clearly higher value for HTO than for iodide. In Essen-1, the apparent diffusion coefficients for iodide and HTO are nearly equal, and slightly smaller than the iodide value in Mol-1/Doel-2b. Accordingly, the HTO apparent diffusion coefficient is considerably smaller in Essen-1 than in Mol-1/Doel-2b.     

Global-scale validation of model-based load deformation of the Earth's crust from continental watermass and atmospheric pressure variations using GPS

Temporal mass variations in the continental hydrosphere and in the atmosphere lead to changes in the gravitational potential field that are associated with load-induced deformation of the Earth’s crust. Therefore, models that compute continental water storage and atmospheric pressure can be validated by measured load deformation time series. In this study, water mass variations as computed by the WaterGAP Global Hydrology Model (WGHM) and surface pressure as provided by the reanalysis product of the National Centers for Environmental Prediction describe the hydrological and atmospheric pressure loading, respectively. GPS observations from 14 years at 208 stations world-wide were reprocessed to estimate admittance factors for the associated load deformation time series in order to determine how well the model-based deformation fits to real data. We found that such site-specific scaling factors can be identified separately for water mass and air pressure loading. Regarding water storage variation as computed by WGHM, weighted global mean admittances are 0.74 ± 0.09, 0.66 ± 0.10, 0.90 ± 0.06 for the north, east and vertical component, respectively. For the dominant vertical component, there is a rather good fit to the observed displacements, and, averaged over all sites, WGHM is found to slightly overestimate temporal variations of water storage. For Europe and North America, with a dense GPS network, site-specific admittances show a good spatial coherence. Regarding regional over- or underestimation of WGHM water storage variations, they agree well with GRACE gravity field data. Globally averaged admittance estimates of pre-computed atmospheric loading displacements provided by the Goddard Geodetic VLBI Group were determined to be 0.88 ± 0.04, 0.97 ± 0.08, 1.13 ± 0.01 for the north, east and vertical, respectively. Here, a relatively large discrepancy for the dominant vertical component indicates an underestimation of corresponding loading predictions.     

Selected trace metals (As,Cd and Hg) distribution and contamination in the coastal wetland sediment of the northern Beibu Gulf,South China Sea

Contamination with As, Cd and Hg, their spatial and temporal distribution are reported from the coastal wetland sediments of the northern Beibu Gulf, South China Sea. The content of As, Cd, Hg and TOC in surface sediments is 8.1 ± 5.8 μg g^?1, 0.08 ± 0.14 μg g^?1, 0.034 ± 0.028 μg g^?1 and 0.45 ± 0.39%, respectively. The mean sedimentation rates are 0.93–1.37 cm year^?1 during 1920s to 2008 determined by ²¹⁰Pb and ¹³⁷Cs dating in three cores. The vertical profiles of As, Cd and Hg content in the cores retrieved from Qin and Nanliu River estuaries show increasing trends during 1985–2008 due to anthropogenic impact caused by local economic development. Locally the surface sediments have potential ecological risk of As to benthos according to the NOAA sediment quality guidelines.     

Oxygen isotope constraints on the origin and differentiation of the Moon

We report new high-precision laser fluorination three-isotope oxygen data for lunar materials. Terrestrial silicates with a range of δ¹⁸O values (− 0.5 to 22.9‰) were analyzed to independently determine the slope of the terrestrial fractionation line (TFL; λ = 0.5259 ± 0.0008; 95% confidence level). This new TFL determination allows direct comparison of lunar oxygen isotope systematics with those of Earth. Values of Δ¹⁷O for Apollo 12, 15, and 17 basalts and Luna 24 soil samples average 0.01‰ and are indistinguishable from the TFL. The δ¹⁸O values of high- and low-Ti lunar basalts are distinct. Average whole-rock δ¹⁸O values for low-Ti lunar basalts from the Apollo 12 (5.72 ± 0.06‰) and Apollo 15 landing sites (5.65 ± 0.12‰) are identical within error and are markedly higher than Apollo 17 high-Ti basalts (5.46 ± 0.11‰). Evolved low-Ti LaPaz mare-basalt meteorite δ¹⁸O values (5.67 ± 0.05‰) are in close agreement with more primitive low-Ti Apollo 12 and 15 mare basalts. Modeling of lunar mare-basalt source composition indicates that the high- and low-Ti mare-basalt mantle reservoirs were in oxygen isotope equilibrium and that variations in δ¹⁸O do not result from fractional crystallization. Instead, these differences are consistent with mineralogically heterogeneous mantle sources for mare basalts, and with lunar magma ocean differentiation models that result in a thick feldspathic crust, an olivine–pyroxene-rich mantle, and late-stage ilmenite-rich zones that were convectively mixed into deeper portions of the lunar mantle. Higher average δ¹⁸O (WR) values of low-Ti basalts compared to terrestrial mid ocean ridge basalts (Δ=0.18‰) suggest a possible oxygen isotopic difference between the terrestrial and lunar mantles. However, calculations of the δ¹⁸O of lunar mantle olivine in this study are only 0.05‰ higher than terrestrial mantle olivine. These observations may have important implications for understanding the formation of the Earth–Moon system.     

In situ cosmogenic 10Be production-rate calibration from the Southern Alps,New Zealand

We present a ¹⁰Be production-rate calibration derived from an early Holocene debris-flow deposit at about 1000 m above sea level in the central Southern Alps, New Zealand, in the mid-latitude Southern Hemisphere. Ten radiocarbon ages on macrofossils from a soil horizon buried by the deposit date the deposit to 9690 ± 50 calendar years before AD2008. Surface ¹⁰Be concentrations of seven large boulders partially embedded in the stable surface of the deposit are tightly distributed, yielding a standard deviation of ～2%. Conversion of the ¹⁰Be measurements to sea level/high-latitude values using each of five standard scaling methods indicates ¹⁰Be production rates of 3.84 ± 0.08, 3.87 ± 0.08, 3.83 ± 0.08, 4.15 ± 0.09, and 3.74 ± 0.08 atoms g^?1 a^?1, relative to the ‘07KNSTD’ ¹⁰Be AMS standard, and including only the local time-integrated production-rate uncertainties. When including a sea level high-latitude scaling uncertainty the overall error is ～2.5% (1σ) for each rate. To test the regional applicability of this production-rate calibration, we measured ¹⁰Be concentrations in a set of nearby moraines deposited before 18 060 ± 200 years before AD2008. The ¹⁰Be ages are only consistent with minimum-limiting ¹⁴C age data when calculated using the new production rates. This also suggests that terrestrial in situ cosmogenic-nuclide production did not change significantly from Last Glacial Maximum to Holocene time in New Zealand. Our production rates agree well with those of a recent calibration study from northeastern North America, but are 12–14% lower than other commonly adopted values. The production-rate values presented here can be used elsewhere in New Zealand for rock surfaces exposed during or since the last glacial period.     

Radiogenic heat production in the lithosphere of Sulu ultrahigh-pressure metamorphic belt

The Chinese Continental Scientific Drilling (CCSD) project is located at the Sulu ultrahigh-pressure metamorphic (UHPM) belt. It offers a unique opportunity for studying the radiogenic heat production of both shallower and deeper rocks. Based on the concentrations of radiogenic elements U, Th and K on 349 samples from main hole of CCSD (CCSD MH), pilot holes and exposures, we determined radiogenic heat productions of all major rock types in the Sulu UHPM belt. Results show the mean values of orthogneiss and paragneiss are respectively 1.65 ± 0.81 and 1.24 ± 0.61 µW m^? 3. Due to different composition and grade of retrogressive metamorphism, the eclogites display significant scatter in radiogenic heat production, ranging from 0.01 to 2.85 µW m^? 3, with a mean of 0.44 ± 0.55 µW m^? 3. The radiogenic heat production in ultramafic rocks also varies within a large range of 0.02 to 1.76 µW m^? 3, and the average turns out to be 0.18 ± 0.31 µW m^? 3. Based on the measurements and crustal petrologic model, the vertical distribution model of heat production in Sulu crust is established. The resulting mean heat production (0.76 µW m^? 3) contributes 24 mW m^? 2 to the surface heat flow. 1-D thermal model indicates that the temperature at the Moho reaches above 750 °C, and the thermal thickness of the lithosphere is ~ 75 km, in good agreement with the geophysical results. The high teat flow (~ 75 mW m^? 2) together with thin lithosphere presents strong support for the extension events during the late Cretaceous and Cenozoic.     

Sea level rise in the north-western part of the Arabian Gulf

Relative sea level variations in the north-western part of the Arabian Gulf have been estimated in the past using no more than 10 to 15 years of observations. In this study, we have almost doubled the period to 28.7 years by examining all available tide gauge data in the area and constructing a mean gauge time-series from seven coastal tide gauges. We found for the period 1979–2007 a relative sea level rise of 2.2 ± 0.5 mm/year. Using the subsidence observed at 6 GPS stations within a radius of 100 km of the tide gauges as an indication of the vertical land motion, the corresponding absolute sea level rise is 1.5 ± 0.8 mm/year that is in agreement with the global estimate of 1.9 ± 0.1 mm/year (Church and White, 2011) for the same studied period. By taking into account the temporal correlations we conclude that previous published results underestimate the true sea level rate uncertainty in this area by a factor of 5–10.