A novel experimental procedure to investigate the biodegradation of NAPL under unsaturated conditions

Soils need to be thoroughly investigated regarding their potential for the natural attenuation of non-aqueous phase liquids (NAPL). Laboratory investigations truly representative of degradation processes in field conditions are difficult to implement for porous media partially saturated with water, NAPL and air. We propose an innovative protocol to investigate degradation processes under steady-state vadose zone conditions. Experiments are carried out in glass columns filled with a sand and, as bacteria source, a soil from a diesel-fuel-polluted site. Water and NAPL (n-hexadecane diluted in heptamethylnonane (HMN)) are added to the porous medium in a two-step procedure using ceramic membranes placed at the bottom of the column. This procedure results, for appropriate experimental conditions, in a uniform distribution of the two fluids (water and NAPL) throughout the column. In a biodegradation experiment non-biodegradable HMN is used to provide NAPL mass, while keeping biodegradable n-hexadecane small enough to monitor its rapid degradation. Biodegradation is followed as a function of time by measuring oxygen consumption, using a respirometer. Degradative activity is controlled by diffusive transfers in the porous network, of oxygen from the gas phase to the water phase and of n-hexadecane from the NAPL phase to the water phase.     

Lattice Boltzmann method with two relaxation times for advection–diffusion equation: Third order analysis and stability analysis

The objectives of this study are to investigate the third order accuracy and linear stability of the lattice Boltzmann method (LBM) with the two-relaxation-time collision operator (LTRT) for the advection–diffusion equation (ADE) and compare the LTRT model with the single-relaxation-time (LBGK) model. While the LBGK has been used extensively, the LTRT appears to be a more flexible model because it uses two relaxation times. The extra relaxation time can be used to improve solution accuracy and/or stability. This study conducts a third order Chapman–Enskog expansion on the LTRT to recover the macroscopic differential equations up to the third order. The dependency of third order terms on the relaxation times is obtained for different types of equilibrium distribution functions (EDFs) and lattices. By selecting proper relaxation times, the numerical dispersion can be significantly reduced. Furthermore, to improve solution accuracy, this study introduces pseudo-velocities to develop new EDFs to reduce the second order numerical diffusion. This study also derives stability domains based on the lattice Peclet number and Courant number for different types of lattices, EDFs and different values of relaxation times, while conducting linear stability analysis on the LTRT. Numerical examples demonstrate the improvement of the LTRT solution accuracy and stability by selecting proper relaxation times, lattice Peclet number and Courant number.     

Assessment of Water Availability in a Central-European River Basin (Elbe) Under Climate Change

 The Elbe region is representative of humid to semi-humid landscapes in Central Europe, where water availability during the summer season is the limiting factor for plant growth and crop yields, especially in the loess areas with high crop productivity having annual precipitation lower than 500 mm. This paper summarizes the results of the first phase of the GLOWA (GLObal WAter)-Elbe project and tries to assess the reliability of water supply in the German part of the Elbe river basin for the next 50 years, a time scale relevant for the implementation of water and land use management plans. One focus of the study was developing scenarios which are consistent with climate and land use changes considering possible uncertainties. The concluding result of the study is that nature and communities in parts of Central Europe will have to deal with considerably lower water resources under scenario conditions.     

Death valley regional ground-water flow model calibration using optimal parameter estimation methods and geoscientific information systems

A regional-scale, steady-state, saturated-zone ground-water flow model was constructed to evaluate potential regional ground-water flow in the vicinity of Yucca Mountain, Nevada. The model was limited to three layers in an effort to evaluate the characteristics governing large-scale subsurface flow. Geoscientific information systems (GSIS) were used to characterize the complex surface and subsurface hydrogeologic conditions of the area, and this characterization was used to construct likely conceptual models of the flow system. Subsurface properties in this system vary dramatically, producing high contrasts and abrupt contacts. This characteristic, combined with the large scale of the model, make zonation the logical choice for representing the hydraulic-conductivity distribution. Different conceptual models were evaluated using sensitivity analysis and were tested by using nonlinear regression to determine parameter values that are optimal, in that they provide the best match between the measured and simulated heads and flows. The different conceptual models were judged based both on the fit achieved to measured heads and spring flows, and the plausibility of the optimal parameter values. One of the conceptual models considered appears to represent the system most realistically. Any apparent model error is probably caused by the coarse vertical and horizontal discretization.     

A closure study of sub-micrometer aerosol particle hygroscopic behaviour

The hygroscopic properties of sub-micrometer aerosol particles were studied in connection with a ground-based cloud experiment at Great Dun Fell, in northern England in 1995. Hygroscopic diameter growth factors were measured with a Tandem Differential Mobility Analyser (TDMA) for dry particle diameters between 35 and 265 nm at one of the sites upwind of the orographic cloud. An external mixture consisting of three groups of particles, each with different hygroscopic properties, was observed. These particle groups were denoted less-hygroscopic, more-hygroscopic and sea spray particles and had average diameter growth factors of 1.11–1.15, 1.38–1.69 and 2.08–2.21 respectively when taken from a dry state to a relative humidity of 90%. Average growth factors increased with dry particle size. A bimodal hygroscopic behaviour was observed for 74–87% of the cases depending on particle size. Parallel measurements of dry sub-micrometer particle number size distributions were performed with a Differential Mobility Particle Sizer (DMPS). The inorganic ion aerosol composition was determined by means of ion chromatography analysis of samples collected with Berner-type low pressure cascade impactors at ambient conditions. The number of ions collected on each impactor stage was predicted from the size distribution and hygroscopic growth data by means of a model of hygroscopic behaviour assuming that only the inorganic substances interacted with the ambient water vapour. The predicted ion number concentration was compared with the actual number of all positive and negative ions collected on the various impactor stages. For the impactor stage which collected particles with aerodynamic diameters between 0.17–0.53 μm at ambient relative humidity, and for which all pertinent data was available for the hygroscopic closure study, the predicted ion concentrations agreed with the measured values within the combined measurement and model uncertainties for all cases but one. For this impactor sampling occasion, the predicted ion concentration was significantly higher than the measured. The air mass in which this sample was taken had undergone extensive photochemical activity which had probably produced hygroscopically active material other than inorganic ions, such as organic oxygenated substances.     

Collateral consequences of the inhomogeneous distribution of short-lived radionuclides in the solar nebula

Abstract— The presence of several short-lived (now extinct) radionuclides in the early solar system demands that they were synthesized and added to preexisting solar system materials shortly (on a time scale on the order of the relevant radionuclide lifetime) before formation of solar system solids. For diverse reasons, it is often suggested that the solar system distributions of these radionuclides were radically heterogeneous, perhaps because of the late addition. Much attention has been given to the astrophysical circumstances that might govern the synthesis and distribution of these short-lived radionuclides, but comparatively little attention has been devoted to the distribution of cosynthesized isotopes. The focus of this paper is a systematic, quantitative evaluation of the collateral consequences in stable and long-lived isotopes that might be expected if short-lived radionuclides, in particular ²⁶Al or ⁵³Mn, were injected at their canonical levels and inhomogeneously distributed in the early solar system. We mix model massive star yields of Meyer et al. (1995) and Woosley and Weaver (1995) into a reservoir of cosmic composition, as tabulated by Anders and Grevesse (1989). To mitigate the effects of systematic deviations that may be present in these mixtures due to uncertainties in model stellar yields, we follow Timmes and Clayton (1996) and also mix into a “renormalized” proxy solar system composition computed from a galactic chemical evolution model based primarily on the stellar yields of Woosley and Weaver (1995). The results are very unfavorable to the likelihood of heterogeneously distributed ²⁶Al derived from supernova ejecta. If a massive star is invoked to account for ²⁶Al, its ejecta must have been rather uniformly distributed, as inferred from the lack of measured collateral anomalies in several elements, notably Ca, Cr, and Ni. Conversely, if ²⁶Al were indeed radically heterogeneously distributed, some other nucleosynthetic source more efficient at producing ²⁶Al is required. In principle, a similar statement applies to ⁵³Mn, but the situation is more complicated. The inferred anomalies at ⁵³Cr will depend not only on how much ⁵³Mn is added by a heterogeneous component, but also more sensitively on the contributions to the associated stable nuclides, ⁵³Cr, ⁵²Cr, and ⁵⁰Cr. Consideration of predicted collateral anomalies provides no direct support for heterogeneously distributed supernova-derived ⁵³Mn, but the required quantity of supernova contribution, and thus also the collateral anomalies, are much less for ⁵³Mn than for ²⁶Al. With allowance for model calculation uncertainties, it could be argued that anomalies collateral to heterogeneous ⁵³Mn might be small enough to have evaded detection.     

Long-term records of erosional change from marine ferromanganese crusts

Ferromanganese crusts from the Atlantic, Indian and Pacific Oceans record the Nd and Pb isotope compositions of the water masses from which they form as hydrogenous precipitates. The¹⁰Be/⁹Be-calibrated time series for crusts are compared to estimates based on Co-contents, from which the equatorial Pacific crusts studied are inferred to have recorded ca. 60 Ma of Pacific deep water history. Time series of ɛNd show that the oceans have maintained a strong provinciality in Nd isotopic composition, determined by terrigenous inputs, over periods of up to 60 Ma. Superimposed on the distinct basin-specific signatures are variations in Nd and Pb isotope time series which have been particularly marked over the last 5 Ma. It is shown that changes in erosional inputs, particularly associated with Himalayan uplift and the northern hemisphere glaciation have influenced Indian and Atlantic Ocean deep water isotopic compositions respectively. There is no evidence so far for an imprint of the final closure of the Panama Isthmus on the Pb and Nd isotopic composition in either Atlantic or Pacific deep water masses.     

Variations of continental discharge pattern in space and time: implications from the Laptev Sea continental margin, Arctic Siberia

Variations in sediment input and distribution to the Laptev Sea continental margin during the Holocene and Termination I could be identified based on radiocarbon dated magnetic susceptibility logs and sediment thickness in high-resolution seismic profiles. Magnetic susceptibility of surface samples reveals an increased input of magnetic grains to the Laptev Sea deriving from the Anabar and Khatanga river catchments. Exposed magnetite schists and volcanic rocks of the Anabar shield and Putoran Plateau, respectively, function as major source of magnetic material. The distribution of magnetic susceptibility in association with the thickness of the Holocene sediments indicates bottom-current induced sediment transport guided by major submarine valleys on the Laptev Sea shelf. The sites of filled paleoriver channels identified in the seismic profiles suggest that during the Late Weichselian sea-level lowstand river runoff continued through four of the major valleys on the exposed Laptev Sea shelf. The sediments at the top of the lowstand deposits in front of the Anabar-Khatanga valley, represented in the seismic profiles by prograding deltas, are characterized by outstandingly high magnetic susceptibility values. Radiocarbon datings approximate the deposition of these high magnetic sediments between 10 and 13.4 ka. It is suggested that this increased input of magnetic material is related to the deglaciation of the Anabar shield and the Putoran Plateau and thus support their glaciation during marine isotope stage (MIS) 2.