Influence of HCl/HF treatment on organic matter in aquifer sediments: A Rock-Eval pyrolysis study

Rock–Eval pyrolysis is increasingly used for the routine characterization of natural organic matter in soils and sediments. In this work the bulk composition of sedimentary organic matter (SOM) in sandy aquifer sediments is studied, as well as purified samples (isolation of SOM) by HCl/HF treatment. This treatment is necessary to avoid detection limit problems for samples with low SOM contents, but the results presented here indicate that this treatment influences the organic geochemistry of the aquifer sediment samples. The FID and CO₂/CO pyrograms show a shift of 10–40 °C of the major peak to a lower temperature. Organic matter alteration or removal of components containing O-bearing groups may explain this. It is also suggested that destruction of the mineral matrix may lead to the reduced retention of the material. For the change of the CO₂/CO pyrograms of the RC fraction only organic matter alteration seems to be likely. Concentrated organic matter samples may also accelerate the release of exothermic energy and influence the pyrograms. Results indicate that the organic matter concentration in the sample influences the measured total organic matter (TOM) content and the T_max of the FID pyrogram, while the sample loading (absolute organic matter amount) up to 80 mg in the Rock–Eval apparatus does not. The FID pyrograms can be deconvoluted into four subpeaks, which allows comparison of samples at various depths. Rock–Eval pyrolysis may only be routinely applied to characterize SOM in aquifer sediments when such systematic and analytical phenomena are taken into account.     

Three-dimensional geostatistical inversion of flowmeter and pumping test data

We jointly invert field data of flowmeter and multiple pumping tests in fully screened wells to estimate hydraulic conductivity using a geostatistical method. We use the steady-state drawdowns of pumping tests and the discharge profiles of flowmeter tests as our data in the inference. The discharge profiles need not be converted to absolute hydraulic conductivities. Consequently, we do not need measurements of depth-averaged hydraulic conductivity at well locations. The flowmeter profiles contain information about relative vertical distributions of hydraulic conductivity, while drawdown measurements of pumping tests provide information about horizontal fluctuation of the depth-averaged hydraulic conductivity. We apply the method to data obtained at the Krauthausen test site of the Forschungszentrum Jülich, Germany. The resulting estimate of our joint three-dimensional (3D) geostatistical inversion shows an improved 3D structure in comparison to the inversion of pumping test data only.     

Modelling the water balance of a mesoscale catchment basin using remotely sensed land cover data

Hydrological modelling of mesoscale catchments is often adversely affected by a lack of adequate information about specific site conditions. In particular, digital land cover data are available from data sets which were acquired on a European or a national scale. These data sets do not only exhibit a restricted spatial resolution but also a differentiation of crops and impervious areas which is not appropriate to the needs of mesoscale hydrological models. In this paper, the impact of remote sensing data on the reliability of a water balance model is investigated and compared to model results determined on the basis of CORINE (Coordination of Information on the Environment) Land Cover as a reference. The aim is to quantify the improved model performance achieved by an enhanced land cover representation and corresponding model modifications. Making use of medium resolution satellite imagery from SPOT, LANDSAT ETM+ and ASTER, detailed information on land cover, especially agricultural crops and impervious surfaces, was extracted over a 5-year period (2000–2004). Crop-specific evapotranspiration coefficients were derived by using remote sensing data to replace grass reference evapotranspiration necessitated by the use of CORINE land cover for rural areas. For regions classified as settlement or industrial areas, degrees of imperviousness were derived. The data were incorporated into the hydrological model GROWA (large-scale water balance model), which uses an empirical approach combining distributed meteorological data with distributed site parameters to calculate the annual runoff components. Using satellite imagery in combination with runoff data from gauging stations for the years 2000–2004, the actual evapotranspiration calculation in GROWA was methodologically extended by including empirical crop coefficients for actual evapotranspiration calculations. While GROWA originally treated agricultural areas as homogeneous, now a consideration and differentiation of the main crops is possible. The accuracy was determined by runoff measurements from gauging stations. Differences in water balances resulting from the use of remote sensing data as opposed to CORINE were analysed in this study using a representative subcatchment. Resulting Nash–Sutcliff model efficiencies improved from 0.372 to 0.775 and indicate that the enhanced model can produce thematically more accurate and spatially more detailed local water balances. However, the proposed model enhancements by satellite imagery have not exhausted the full potential of water balance modelling, for which a higher temporal resolution is required.     

Ion microprobe uranium‐thorium‐lead dating of Shergotty phosphates

Abstract— We report ion microprobe U‐Th‐Pb dating of Shergotty phosphates by means of the sensitive high‐resolution ion microprobe (SHRIMP) recently installed at Hiroshima University, Japan. ten analyses of whitlockite (merrillite) and three analyses of apatite indicate a ²³⁸u/²⁰⁶pb isochron age of 225 ± 200 ma and a tera‐wasserburg concordia‐constrained linear three‐dimensional isochron age of 217 ± 110 ma in the ²³⁸u/²⁰⁶pb‐²⁰⁷pb/²⁰⁶pb²⁰⁴pb/²⁰⁶pb diagram. These ages agree well with the 232Th‐208pb age of 189 ± 83 Ma, which suggests that primary crystallization or a shock metamorphic event defined the formation age of the phosphate minerals. The average of the later two ages, 204 ± 68 Ma, is consistent with the previously published Rb‐Sr age of 165 ± 11 Ma and U‐Th‐Pb age of ～200 Ma. These show marginal agreement with the ⁴⁰Ar‐³⁹Ar age of 254 ± 10 Ma but are significantly different from the Sm‐Nd age of 360 ± 16 Ma. Taking into account the closure temperature of the U‐Pb system in apatite, we suggest the time that Shergotty last experienced a temperature of ～900 °C was 204 ± 68 Ma.     

Crystallization of the basaltic shergottites: Insights from crystal size distribution (CSD) analysis of pyroxenes

Abstract— Quantitative petrographic analysis, using the crystal size distribution (CSD) method, provides a novel approach for examining the crystallization histories of basaltic shergottites. Grain number densities at different sizes are plotted against grain size, and the resulting curve relates to the geologic processes involved with the crystallization of the grain population. Most basaltic shergottites are dominated by pigeonite and augite; and because plagioclase is primarily interstitial, and therefore constrained in its growth by the surrounding pyroxenes, we limited our size measurements to the pyroxene phases. The groundmasses of Elephant Moraine (EET) A79001 lithology A and Dar al Gani (DaG) 476 are fine grained with cumulus pyroxene and interstitial plagioclase glass. Their simple linear CSD plots record a single stage of pyroxene crystallization under steady‐state conditions of continuous nucleation and growth. The textures of Queen Alexandra Range (QUE) 94201 and EETA79001 lithology B are quite different from the other shergottites, with intergrown pyroxene and plagioclase. Likewise, their CSD plots are also distinct, with curved trends that suggest a lack of large grains, most likely because of interference between simultaneously growing silicate phases. However, the CSD plot shapes are smooth, also implying a single stage of growth. Shergotty and Zagami, with coarser cumulus textures, display CSD plots that are generally linear over most grain sizes. This implies that conditions of nucleation and growth were dominant during formation of the pyroxene populations. Both plots, however, also display kinks, implying multiple stages of growth. A similar kink is also visible in a CSD plot of only the Mg‐rich cores of Shergotty pyroxenes, which suggests the feature represents changes in conditions during core crystallization, rather than an event coincident with the change in composition to the Fe‐rich rims. The plot may be interpreted as representing two stages of core growth with an intervening short hiatus of nucleation, with continued crystallization associated with ascent of the magma. Eruption onto the surface probably triggered the compositional change to Fe‐rich rims. The CSD analysis of products from a controlled crystallization study agree with experimental and petrologic estimates that cooling rates for Zagami were on the order of a few tenths of a degree per hour. Growth rates derived from these cooling rates suggest crystallization of Shergotty and Zagami pyroxenes occurred over a period of a few weeks to months.     

On the origin of rim textures surrounding anhydrous silicate grains in CM carbonaceous chondrites

Abstract— Fine‐grained, optically opaque rims coat individual olivine and pyroxene grains in CM matrices and chondrules. Bulk chemical analyses and observations of these rims indicate the presence of phyllosilicates and disseminated opaques. Because phyllosilicates could not have survived the chondrule formation process, chondrule silicate rims must have formed entirely by late‐state aqueous reactions. As such, these textures provide a useful benchmark for isolating alteration features from more complex CM matrix materials. Both chondrule silicate and matrix silicate rims exhibit morphological features commonly associated with advancing stages of replacement reactions in terrestrial serpentinites. Contacts between many matrix silicate rims and the adjacent matrix materials suggest that these rims formed entirely by aqueous reactions in a parent‐body setting. This contrasts with previous assertions that rim textures can only form by the accretion of nebular dust but does not imply an origin for the rims surrounding other types of CM core components, such as chondrules.     

Rapid changes of sediment dynamic processes in Yalu River Estuary under anthropogenic impacts

The response of the Yalu River Estuary to human activities was investigated.Changes of sediment dynamics during the past 10 years were explored through hydrodynamic calculation,as well as heavy mineral...     

Lava channel roofing,overflows, breaches and switching: insights from the 2008–2009 eruption of Mt. Etna

During long-lived basaltic eruptions, overflows from lava channels and breaching of channel levées are important processes in the development of extensive 'a'ā lava flow-fields. Short-lived breaches result in inundation of areas adjacent to the main channel. However, if a breach remains open, lava supply to the original flow front is significantly reduced, and flow-field widening is favoured over lengthening. The development of channel breaches and overflows can therefore exert strong control over the overall flow-field development, but the processes that determine their location and frequency are currently poorly understood. During the final month of the 2008–2009 eruption of Mt. Etna, Sicily, a remote time-lapse camera was deployed to monitor events in a proximal region of a small ephemeral lava flow. For over a period of ~10 h, the flow underwent changes in surface elevation and velocity, repeated overflows of varying vigour and the construction of a channel roof (a required prelude to lava tube formation). Quantitative interpretation of the image sequence was facilitated by a 3D model of the scene constructed using structure-from-motion computer vision techniques. As surface activity waned during the roofing process, overflow sites retreated up the flow towards the vent, and eventually, a new flow was initiated. Our observations and measurements indicate that flow surface stagnation and flow inflation propagated up-flow at an effective rate of ~6 m h⁻¹, and that these processes, rather than effusion rate variations, were ultimately responsible for the most vigorous overflow events. We discuss evidence for similar controls during levée breaching and channel switching events on much larger flows on Etna, such as during the 2001 eruption.     

Simulations of spatial variability in particle‐size emissions during wind erosion events

This study simulates how spatial variations in particle‐size emissions from a playa affect bulk and size‐resolved dust concentration profiles during two contrasting wind erosion events (a small local and a large regional event) in the Channel Country, Lake Eyre Basin, Australia. The regional event had higher dust concentration as a result of stronger frontal winds and higher erodibility across the playa. For each event, two emission scenarios are simulated to determine if measured size‐resolved dust concentration profiles can be explained by spatial variability in source area emissions. The first scenario assumes that particle‐size emissions from source areas occur at a uniform rate, while the second scenario assumes that particle‐size emissions vary between and within source areas. The uniform emission scenario, reproduced measured bulk dust concentration profiles (R² = 0·93 regional and R² = 0·81 local), however simulated size‐resolved dust concentration profiles had poor statistical fits to measured size‐resolved profiles for each size class (the highest were R² = 0·5 regional and R² = 0·3 local). For the differential particle‐size emission scenario, the fit to the measured bulk dust concentration profiles is improved (R² = 0·97 regional and R² = 0·83 local). However, the fit to the size‐resolved profiles improved dramatically, with the lowest being R² = 0·89 (regional) and R² = 0·80 (local). Particle‐size emission models should therefore be tested against both bulk and size‐resolved dust concentration profiles, since if only bulk dust concentration profiles are used model performance may be over‐stated. As the source areas in the first 90 m upwind of the tower were similar for both events, the percentage contributions of each particle‐size class to total emissions can be compared. The contribution of each particle‐size class was similar even though the wind speed, turbulence and dust concentrations were significantly different; suggesting that the contribution of each particle‐size to the total emitted dusts is not related to wind speed and turbulence. Copyright © 2012 John Wiley & Sons, Ltd.