On the cumulant expansion up scaling of ground water contaminant transport equation with nonequilibrium sorption

The laboratory-scale ground water transport equation with nonequilibrium sorption reaction subjected to unsteady, nondivergence-free, and nonstationary velocity fields is up-scaled to the field-scale by using the ensemble-averaged equations obtained from the cumulant expansion ensemble-averaging method. It is found that existing ensemble-averaged equations obtained with the help of the cumulant expansion method for the system of linear partial differential equations are not second-order exact. Although the cumulant expansion methodology is designed for noncommuting operators, it is found that there are still commudativity requirements that need to be satisfied by the functions and constants exist in the coefficient matrix of the system of ordinary/partial differential equations. A reversibility requirement, which covers the commudativity requirements, is also proposed when applying the cumulant expansion method to a system of partial differential equations/a partial differential equation. The significance of the new velocity correction obtained in this study due to the applied second-order exact cumulant expansion is investigated on a numerical example with a linear trend in the distribution coefficient. It is found that the effect of the new velocity correction can be significant enough to affect the maximum concentration values and the plume center of mass in the case of a trending distribution coefficient in a physically heterogeneous environment.     

IMPACT: an implicit time integration scheme for chemical species and families

The implicit time integration scheme of Stott and Harwood (1993) was proposed as an efficient scheme for use in three-dimensional chemical models of the atmosphere. The scheme was designed for chemistry schemes using chemical families, in which species with short lifetimes are grouped into longer-lived families. Further study with more complex chemistry, more species and reactions showed the scheme to be non-convergent and unstable under certain conditions; particularly for the perturbed chemical scenarios of polar stratospheric winters. In this work the scheme has been improved by revising the treatment of families and the convergence properties of the scheme. The new scheme has been named IMPACT (IMPlicit Algorithm for Chemical Time-stepping). It remains easy to implement and produces simulations that compare well with integrations using more accurate higher order schemes.     

An Open,Object-Based Framework for Generating Anisotropy in Sedimentary Subsurface Models

The spatial distribution of hydraulic properties in the subsurface controls groundwater flow and solute transport. However, many approaches to modeling these distributions do not produce geologically realistic results and/or do not model the anisotropy of hydraulic conductivity caused by bedding structures in sedimentary deposits. We have developed a flexible object-based package for simulating hydraulic properties in the subsurface—the Hydrogeological Virtual Realities (HyVR) simulation package. This implements a hierarchical modeling framework that takes into account geological rules about stratigraphic bounding surfaces and the geometry of specific sedimentary structures to generate realistic aquifer models, including full hydraulic-conductivity tensors. The HyVR simulation package can create outputs suitable for standard groundwater modeling tools (e.g., MODFLOW), is written in Python, an open-source programming language, and is openly available at an online repository. This paper presents an overview of the underlying modeling principles and computational methods, as well as an example simulation based on the Macrodispersion Experiment site in Columbus, Mississippi. Our simulation package can currently simulate porous media that mimic geological conceptual models in fluvial depositional environments, and that include fine-scale heterogeneity in distributed hydraulic parameter fields. The simulation results allow qualitative geological conceptual models to be converted into digital subsurface models that can be used in quantitative numerical flow-and-transport simulations, with the aim of improving our understanding of the influence of geological realism on groundwater flow and solute transport.     

Schotcrete as permanent lining for the Furka base tunnel

Summary The lining concept of the 15.4 km long one-track Furka railway tunnel is described with particular emphasis on shotcrete application as a permanent support Altogether 70000 m³ dry aggregate/ cement mix and 115000 rock bolts have been applied. The tunnel was lined at a rate of approx. 400 tunnel metres per month using up to five shotcrete machines simultaneously. Systematic quality control throughout the construction work indicated an average uniaxial compressive strength of shotcrete of 40 N/mm². The support requirements were grouped into 16 different classes and an additional special class accounting for extreme geological conditions in some relatively short stretches.     

Albedos and size distribution of meteoroids from 0.3 to 4.8 AU

Comparison is made between the run of number density of meteoroids from penetration detectors aboard Helios A (masses below 10^?8 g) and Pioneer 10 (masses near and above 3 × 10^?9 g), the source function of the zodiacal light deduced from photometric observations aboard Helios A and Pioneer 10, counts versus brightness of objects passing by Pioneer 10 from the Sisyphus experiment and the distribution of meteoroids deduced from radar and optical meteors at the Earth. The Sisyphus experiment on Pioneer 10 observed reflecting glints on meteoroids rather than the meteoroids themselves and the counting statistics refer not to the effective radii of the meteoroids but to the effective radii of curvature of the reflecting glints on the meteoroids. The penetration detectors appear to find some increase in number density toward the Sun and a flat distribution outward to 5.2 AU. The overall behavior of the zodiacal light is that the relative distribution over direction is unchanged while the source scattering function diminishes as the inverse 1.4 power of distance from the Sun. The fit to the brightness of the zodiacal light obtained from these statistics can be combined with the mass distribution results from the optical meteors to deduce a mean geometric albedo of meteoroids of 0.006 at 1 AU from the Sun. Combination of the space distribution from radar meteors with the scattering source function of the zodiacal light yields geometric albedos for meteoroids running from 0.07 at 0.1 AU, from the Sun through 0.006 at 1 AU down to about 0.0001 at 3.3 AU which may run flat thence outward. This result is imposed by the indicated modest increase in density of meteoroids very near the Sun, a minimum between the Sun and the Earth near 0.4 AU and rising density outward to somewhere beyond 3.3 AU which is very different from the inverse 1.4 power of the distance shown for scatterers (product of number density and albedo) by the zodiacal light. A check on the distribution at very large sizes is possible if a search is made for fireballs in Jupiter's atmosphere by the Mariner Jupiter Saturn 1977 television cameras during the two encounters with Jupiter in 1979. An easy detection of such activity would put the maximum in the meteoroid distribution out near Jupiter and lend further confirmation to the indicated drop in albedo.     

Spectral reflectance of Martian meteorites: Spectral signatures as a template for locating source region on Mars

Abstract— We report the spectral reflectance of Martian meteorites from 0.3‐2.6 microns for the purpose of cataloguing spectra and the association of their spectral properties with mineralogy and petrology. We fit the spectra to a series of overlapping, modified Gaussian absorptions using least squares fitting. The results are validated against established relationships between photon interactions with mineral chemistry and the band parameters. These resultant band parameters can be used to constrain interpretations of Martian reflectance spectra in the search for the source region of meteorites from Mars. The limitations of the fitting method are discussed.     

Low-temperature thermodynamic properties of disordered zeolites of the natrolite group

 The heat capacity of paranatrolite and tetranatrolite with a disordered distribution of Al and Si atoms has been measured in the temperature range of 6–309 K using the adiabatic calorimetry technique. The composition of the samples is represented with the formula (Na_1.90K_0.22Ca_0.06)[Al_2.24Si_2.76O₁₀]·nH₂O, where n=3.10 for paranatrolite and n=2.31 for tetranatrolite. For both zeolites, thermodynamic functions (vibrational entropy, enthalpy, and free energy function) have been calculated. At T=298.15 K, the values of the heat capacity and entropy are 425.1 ± 0.8 and 419.1 ±0.8 J K⁻¹ mol⁻¹ for paranatrolite and 381.0 ± 0.7 and 383.2 ± 0.7 J K⁻¹ mol⁻¹ for tetranatrolite. Thermodynamic functions for tetranatrolite and paranatrolite with compositions corrected for the amount of extraframework cations and water molecules have also been calculated. The calculation for tetranatrolite with two water molecules and two extraframework cations per formula yields: C _p (298.15)=359.1 J K⁻¹ mol⁻¹, S(298.15) −S(0)=362.8 J K⁻¹ mol⁻¹. Comparing these values with the literature data for the (Al,Si)-ordered natrolite, we can conclude that the order in tetrahedral atoms does not affect the heat capacity. The analysis of derivatives dC/dT for natrolite, paranatrolite, and tetranatrolite has indicated that the water- cations subsystem within the highly hydrated zeolite may become unstable at temperatures above 200 K. Received: 30 July 2001 / Accepted: 15 November 2001