Isotope variations in a Sierra Nevada snowpack and their relation to meltwater

Isotopic variations in melting snow are poorly understood. We made weekly measurements at the Central Sierra Snow Laboratory, California, of snow temperature, density, water equivalent and liquid water volume to examine how physical changes within the snowpack govern meltwater δ¹⁸O. Snowpack samples were extracted at 0.1 m intervals from ground level to the top of the snowpack profile between December 1991 and April 1992. Approximately 800 mm of precipitation fell during the study period with δ¹⁸O values between −21.35 and −4.25‰. Corresponding snowpack δ¹⁸O ranged from −22.25 to −6.25‰. The coefficient of variation of δ¹⁸O in snowpack levels decreased from −0.37 to −0.07 from winter to spring, indicating isotopic snowpack homogenization. Meltwater δ¹⁸O ranged from −15.30 to −8.05‰, with variations of up to 2.95‰ observed within a single snowmelt episode, highlighting the need for frequent sampling. Early snowmelt originated in the lower snowpack with higher δ¹⁸O through ground heat flux and rainfall. After the snowpack became isothermal, infiltrating snowmelt displaced the higher δ¹⁸O liquid in the lower snowpack through a piston flow process. Fractionation analysis using a two-component mixing model on the isothermal snowpack indicated that δ¹⁸O in the initial and final half of major snowmelt was 1.30‰ lower and 1.45‰ higher, respectively, than the value from simple mixing. Mean snowpack δ¹⁸O on individual profiling days showed a steady increase from −15.15 to −12.05‰ due to removal of lower δ¹⁸O snowmelt and addition of higher δ¹⁸O rainfall. Results suggest that direct sampling of snowmelt and snow cores should be undertaken to quantify tracer input compositions adequately. The snowmelt sequence also suggests that regimes of early lower δ¹⁸O and later higher δ¹⁸O melt may be modeled and used in catchment tracing studies.     

Predicting average annual groundwater levels from climatic variables: an empirical model

On the basis of one-dimensional theoretical water flow model, we demonstrate that the groundwater level variation follows a pattern similar to recharge fluctuation, with a time delay that depends on the characteristics of aquifer, recharge pattern as well as the distance between the recharge and observation locations. On the basis of a water budget model and the groundwater flow model, we propose an empirical model that links climatic variables to groundwater level. The empirical model is tested using a partial data set from historical records of water levels from more than 80 wells in a monitoring network for the carbonate rock aquifer, southern Manitoba, Canada. The testing results show that the predicted groundwater levels are very close to the observed ones in most cases. The overall average correlation coefficient between the predicted and observed water levels is 0.92. This proposed empirical statistical model could be used to predict variations in groundwater level in response to different climate scenarios in a climate change impact assessment.     

Seismic scattering at the top of the mantle Transition Zone

We document strong seismic scattering from around the top of the mantle Transition Zone in all available high resolution explosion seismic profiles from Siberia and North America. This seismic reflectivity from around the 410 km discontinuity indicates the presence of pronounced heterogeneity in the depth interval between 320 and 450 km in the Earth’s mantle. We model the seismic observations by heterogeneity in the form of random seismic scatterers with typical scale lengths of kilometre size (10-40 km by 2-10 km) in a 100-140 km thick depth interval. The observed heterogeneity may be explained by changes in the depths to the α-β-γ spinel transformations caused by an unexpectedly high iron content at the top of the mantle Transition Zone. The phase transformation of pyroxenes into the garnet mineral majorite probably also contributes to the reflectivity, mainly below a depth of 400 km, whereas we find it unlikely that the presence of water or partial melt is the main cause of the observed strong seismic reflectivity. Subducted oceanic slabs that equilibrated at the top of the Transition Zone may also contribute to the observed reflectivity. If this is the main cause of the reflectivity, a substantial amount of young oceanic lithosphere has been subducted under Siberia and North America during their geologic evolution. Subducted slabs may have initiated metamorphic reactions in the original mantle rocks.     

Flat mantle reflectors in Eastern China: possible evidence of lithospheric thinning

Recent 24 s deep seismic reflection records revealed five flat reflectors in the lithospheric mantle in Eastern China. With increasing depth, they are named M1 to M5 and can be seen on both field single-shot and stacked records. Reflector M1 corresponds to the Moho discontinuity, whereas M5 may be the reflection from the bottom of the current lithosphere, which is about 78 km deep according to geothermal measurements. The other three reflectors seem peculiar and might result from interactions between the lithosphere and deeper mantle. Based on lithological and geochemical data, it is suggested that the lithosphere has been thinned from about 150 km to about 60 km in the Late Mesozoic, and then has been thickened to about 78 km during the Cenozoic. The thinning process produced a granulite layer in the old lower crust caused by magmatic underplating, whereas an eclogite layer formed beneath owing to the subduction of the Paleo-Tethys and Yangtze Craton during the Permian and Early Mesozoic. Reflector M2 at about 12 s two-way traveltime (TWT) might result from the Paleozoic Moho, which represents the boundary between the previous granulite and eclogite facies. Reflector M3 at about 14 s might correspond to the bottom of the eclogite layer, beneath which the old lithospheric mantle remained. The old and the newly developed mantle may have different compositions, resulting in reflector M4. The multi-layered mantle reflectors demonstrate a mantle structure that possibly correlates with the lithospheric thinning process that occurred in Eastern China during the Late Mesozoic. The discovery of multi-layered mantle reflectors in the studied areas indicates a high heterogeneity of the upper mantle. Reflection seismology with improved technology, together with velocity and resistivity imaging and rock-physics measurements, can provide more details of the heterogeneity and related dynamic processes that occurred in the lithospheric mantle.     

Crossover from Nonlinearity Controlled to Heterogeneity Controlled Mixing in Two-Phase Porous Media Flows

We use high resolution Monte Carlo simulations to study the dispersive mixing in two-phase, immiscible, porous media flow that results from the interaction of the nonlinearities in the flow equations with geologic heterogeneity. Our numerical experiments show that distinct dispersive regimes occur depending on the relative strength of nonlinearity and heterogeneity. In particular, for a given degree of multiscale heterogeneity, controlled by the Hurst exponent which characterizes the underlying stochastic model for the heterogeneity, linear and nonlinear flows are essentially identical in their degree of dispersion, if the heterogeneity is strong enough. As the heterogeneity weakens, the dispersion rates cross over from those of linear heterogeneous flows to those typical of nonlinear homogeneous flows.     

Time-domain computation of scattering from a heterogeneous layer

Adopting Born and ray approximations, time-domain synthetic seismograms for P-P and P-S scattering from a plane wave incident on a thin, laterally heterogeneous layer are presented in this paper. The time-domain P coda is a convolution between a structure function and the second-order derivative of the time function of the incident P wave. Examples of synthetic seismograms are given using a time function from a computed short-period seismogram for a point explosive source in a half-space. These show that it is impossible, with realistic values of the parameters involved, to generate significant codas when only single scattering is involved.     

Genesis and distribution of mineral waters as a consequence of recent lithospheric dynamics: the Rhenish Massif,Central Europe

Three major, interdependent processes control the genesis and distribution of mineral and thermal waters in the Rhenish Massif, Central Europe: (a) Magmatic processes in the upper mantle provide most of the CO₂ to produce bicarbonate waters in shallow aquifers. (b) Extension of the brittle upper crust enables the ascent of sodium chloride waters. (c) Uplift and erosion shape the massif's relief, which determines the extent of flow systems and the distribution of thermal springs. The chemistry of mineral waters further depends on the aquifers' mineral composition. A comprehensive set of hydrological, chemical, tectonic and geophysical data on the Rhenish Massif has been compiled. It was used to classify the mineral waters and to map the spatial distribution of water properties. The composition of cuttings from several representative wells producing different water types shows that the hydrothermal alteration of the aquifer rocks consists mainly of kaolinization of chlorite and dissolution of feldspar. Numerical transport simulations favour two modes of groundwater flow: topography-driven flow and the pressure-driven ascent of basement brines along active faults. Thermal convection is less important.     

Flexible spectral methods for the generation of random fields with power-law semivariograms

Random field generators serve as a tool to model heterogeneous media for applications in hydrocarbon recovery and groundwater flow. Random fields with a power-law variogram structure, also termed fractional Brownian motion (fBm) fields, are of interest to study scale dependent heterogeneity effects on one-phase and two-phase flow. We show that such fields generated by the spectral method and the Inverse Fast Fourier Transform (IFFT) have an incorrect variogram structure and variance. To illustrate this we derive the prefactor of the fBm spectral density function, which is required to generate the fBm fields. We propose a new method to generate fBm fields that introduces weighting functions into the spectral method. It leads to a flexible and efficient algorithm. The flexibility permits an optimal choice of summation points (that is points in frequency space at which the weighting function is calculated) specific for the autocovariance structure of the field. As an illustration of the method, comparisons between estimated and expected statistics of fields with an exponential variogram and of fBm fields are presented. For power-law semivariograms, the proposed spectral method with a cylindrical distribution of the summation points gives optimal results.     

Observations of particle movement in a monitoring well using the colloidal borescope

The colloidal borescope consists of a set of lenses and miniature video cameras capable of observing natural particles in monitoring wells. Based on field observations of these particles, it appears possible to measure in situ groundwater velocity in a well bore. Field observations have shown that directional measurements using the colloidal borescope are generally in good agreement with expected flow directions. However, the magnitude of flow velocity is higher compared with values based on conventional test methods. High relative flow velocities, even after correction factors have been applied to compensate for well bore effects, are believed to be due to preferential flow zones in the surrounding aquifer. Low flow zones exhibit swirling multidirectional flow that does not allow for a linear velocity measurement. Consequently, groundwater flow velocities measured by the colloidal borescope in heterogeneous aquifers will be biased toward the maximum velocity values present in the aquifer. A series of laboratory experiments was conducted to assess the reliability of the instrument. Based on this work, a seepage velocity correction factor ( ) of 1–4 was found for quantifying groundwater seepage velocity in the adjacent aquifer from observations in a well bore. Laboratory measurements also indicate that preferential flow in the surrounding aquifer dominates flow in the well. Results of this work suggest the possibility of quantifying higher-flow velocities associated with preferential flow zones in the subsurface.