Isotopes and sustainability of ground water resources, North China Plain

Ground water in deep confined aquifers is one of the major water resources for agricultural, industrial, and domestic uses in the North China Plain. Detailed information on ground water age and recharge is vital for the proper management of these water resources, and to this end, we used carbon 14 of dissolved inorganic carbon and tritium in water to measure the age and determine the recharge areas of ground water in the North China Plain. These isotopic data suggest that most ground water in the piedmont part of the North China Plain is <40 years old and is recharged locally. In contrast, ground water in the central and littoral portions of the North China Plain is 10,000 to 25,000 years old. The delta18O (deltaD) values of this ground water are 1.7 per thousand (11 per thousand) less than that in the piedmont plain ground water and possibly reflect water recharged during a cooler climate during the last glaciation. The temperature of this recharge, based on delta18O values, ranges from 3.7 degrees C to 8.4 degrees C, compared to 12 degrees C to 13 degrees C of modern recharge water. The isotopic data set combined indicates that ground water in the central and littoral part of the North China Plain is being mined under non-steady state conditions.     

Spatial and temporal variability of ground water recharge in central Australia: a tracer approach

Two environmental tracer methods are applied to the Ti-Tree Basin in central Australia to shed light on the importance of recharge from floodouts of ephemeral rivers in this arid environment. Ground water carbon-14 concentrations from boreholes are used to estimate the average recharge rate over the interval between where the ground water sample first entered the saturated zone and the bore. Environmental chloride concentrations in ground water samples provide estimates of the recharge rate at the exact point in the landscape where the sample entered the saturated zone. The results of the two tracer approaches indicate that recharge rates around one of the rivers and an extensive floodplain are generally higher than rates of diffuse recharge that occurs in areas of lower topographic relief. Ground water 2H/1H and 18O/16O compositions are all depleted in the heavier isotopes (delta2H = -67 per thousand to -50 per thousand; delta18O = -9.2 per thousand to -5.7%o) compared with the long-term, amount-weighted mean isotopic composition of rainfall in the area (delta2H = -33.8 per thousand; delta18O = -6.3 per thousand). This indicates that recharge throughout the basin occurs only after intense rainfall events of at least 150 to 200 mm/month. Finally, a recharge map is developed to highlight the spatial extent of the two recharge mechanisms. Floodout recharge to the freshest ground water (TDS <1,000 mg/L) is approximately 1.9 mm/year compared with a mean recharge rate of approximately 0.2 mm/year to the remainder of the basin. These findings have important implications for management of the ground water resource.     

Ground water recharge and flow characterization using multiple isotopes

Stable isotopes of delta(18)O, delta(2)H, and (13)C, radiogenic isotopes of (14)C and (3)H, and ground water chemical compositions were used to distinguish ground water, recharge areas, and possible recharge processes in an arid zone, fault-bounded alluvial aquifer. Recharge mainly occurs through exposed stream channel beds as opposed to subsurface inflow along mountain fronts. This recharge distribution pattern may also occur in other fault-bounded aquifers, with important implications for conceptualization of ground water flow systems, development of ground water models, and ground water resource management. Ground water along the mountain front near the basin margins contains low delta(18)O, (14)C (percent modern carbon [pmC]), and (3)H (tritium units [TU]), suggesting older recharge. In addition, water levels lie at greater depths, and basin-bounding faults that locally act as a flow barrier may further reduce subsurface inflow into the aquifer along the mountain front. Chemical differences in ground water composition, attributed to varying aquifer mineralogy and recharge processes, further discriminate the basin-margin and the basin-center water. Direct recharge through the indurated sandstones and mudstones in the basin center is minimal. Modern recharge in the aquifer is mainly through the broad, exposed stream channel beds containing coarse sand and gravel where ground water contains higher delta(18)O, (14)C (pmC), and (3)H (TU). Spatial differences in delta(18)O, (14)C (pmC), and (3)H (TU) and occurrences of extensive mudstones in the basin center suggest sluggish ground water movement, including local compartmentalization of the flow system.     

Boron isotopes as an artificial tracer

A field study was conducted using a combination of intrinsic and artificial tracers to estimate travel times and dilution during transport of infiltrate from a reclaimed water infiltration basin to nearby monitoring wells. A major study objective was to validate boric acid enriched in (10)B as an artificial tracer. Basin 10E at the Rio Hondo Spreading Grounds in Whittier, California, was the site of the test. The basin normally receives a mixture of treated municipal waste water, purchased State Project water, and local runoff from the San Gabriel River. Approximately 3.5 kg of (10)B-enriched boric acid was dispersed among 2.05 x 10(5) m(3) of basin water to initiate the experiment. The resultant median delta(11)B in the infiltration basin was -71 per thousand. Prior to tracer addition, the basin water had an intrinsic delta(11)B of +2 per thousand. Local monitoring wells that were used to assess travel times had delta(11)B values of +5 per thousand and +8 per thousand at the time of tracer addition. Analytic results supported an assumption that boron is conserved during ground water transport and that boron enriched in (10)B is a useful artificial tracer. Several intrinsic tracers were used to reinforce the boric acid tracer findings. These included stable isotopes of oxygen (delta(18)O) and hydrogen (deltaD), sulfate concentration, and the boron to chloride ratio. Xenon isotopes, (136)Xe and (124)Xe, also supported boron isotope results. Xenon isotopes were added to the recharge basin as dissolved gases by investigators from the Lawrence Livermore National Laboratory.     

Use of NADP archive samples to determine the isotope composition of precipitation: characterizing the meteoric input function for use in ground water studies

Stable oxygen and hydrogen isotopes have been used in ground water studies to investigate recharge, mixing, ground water/surface water interaction, advective-diffusive transport, paleohydrogeologic interactions and to estimate ground water ages. Such studies require that the isotopic composition of precipitation be known, as precipitation is a major input to ground water and surface water systems. As oxygen-18 and deuterium data for precipitation are lacking across much of the United States, there is need to establish additional local meteoric water lines as isotope input functions across the region, as well as to develop better understanding of the isotopic climate linkages that control oxygen and hydrogen isotope ratios in precipitation. In the absence of long-term monitoring stations, one possible solution to this problem is to determine the delta 18O and delta 2H values of precipitation using archive samples collected at monitoring stations managed by the National Atmospheric Deposition Program (NADP). This study describes and interprets the seasonal delta 18O and delta 2H composition of archived precipitation samples collected in eastern Nebraska near the town of Mead during the years 1992-1994. Values for delta 18O range from -23.6 to -0.7@1000. Values for delta 2H range from -172 to 0@1000. Yearly arithmetic mean delta 18O and delta 2H values for the Mead station are -8.1@1000 and -53@1000, respectively. Weighted yearly means for delta 18O and delta 2H were -7.4@1000 and -48@1000, respectively. Mead values show a strong isotopic enrichment between winter and summer precipitation, and a strong delta 18O-T correlation (r2 = 0.91) for mean monthly values of about 0.5@1000 per degree Celsius. The local meteoric water line for the Mead site is delta 2H = 7.40 delta 18O + 7.32. Deuterium excess values suggest that most of the moisture across the region is derived primarily from a Gulf of Mexico source. The results of this study demonstrate that in the absence of long-term monitoring stations such as those operated globally by the International Atomic Energy Association, NADP archive samples can be used to determine the isotopic composition of precipitation, to characterize the local meteoric water line and establish the various climatic relationships, and define the meteoric input function for use in ground water studies.     

Environmental isotopes as indicators for ground water recharge to fractured granite

To assess the contribution of accumulated winter precipitation and glacial meltwater to the recharge of deep ground water flow systems in fracture crystalline rocks, measurements of environmental isotope ratios, hydrochemical composition, and in situ parameters of ground water were performed in a deep tunnel. The measurements demonstrate the significance of these ground water recharge components for deep ground water flow systems in fractured granites of a high alpine catchment in the Central Alps, Switzerland. Hydrochemical and in situ parameters, as well as delta(18)O in ground water samples collected in the tunnel, show only small temporal variations. The precipitation record of delta(18)O shows seasonal variations of approximately 14% and a decrease of 0.23% +/- 0.03% per 100 m elevation gain. delta(2)H and delta(18)O in precipitation are well correlated and plot close to the meteoric water line, as well as delta(2)H and delta(18)O in ground water samples, reflecting the meteoric origin of the latter. The depletion of 18O in ground water compared to 18O content in precipitation during the ground water recharge period indicates significant contributions from accumulated depleted winter precipitation to ground water recharge. The hydrochemical composition of the encountered ground water, Na-Ca-HCO3-SO4(-F), reflects an evolution of the ground water along the flowpath through the granite body. Observed tritium concentrations in ground water range from 2.6 to 16.6 TU, with the lowest values associated with a local negative temperature anomaly and anomalous depleted 18O in ground water. This demonstrates the effect of local ground water recharge from meltwater of submodern glacial ice. Such localized recharge from glaciated areas occurs along preferential flowpaths within the granite body that are mainly controlled by observed hydraulic active shear fractures and cataclastic faults.     

Sources of ground water salinity on islands using 18O, 2H, and 34S

Stable isotopes of 18O and 2H in water, and 34S and 18O in dissolved SO4, are used to verify the interpretation of the chemical evolution and proposed sources of salinity for two islands that have undergone postglacial rebound. Results for delta18O and delta34S in dissolved SO4 on the Gulf Islands, southwest British Columbia, Canada, suggest a three-component mixing between (1) atmospheric SO4 derived largely from recharge of meteoric origin, (2) modern marine SO4 associated with either modern-day salt water intrusion or Pleistocene age sea water, and (3) terrestrial SO4. The age of the marine SO4 is uncertain based on the geochemistry and SO4 isotopes alone. Two options for mixing of saline ground waters are proposed--either between current-day marine SO4 and atmospheric SO4, or between older (Pleistocene age) marine SO4 and atmospheric SO4, delta18O and delta2H compositions are relatively consistent between both islands, with a few samples showing evidence of mixing with water that is a hybrid mixture of Fraser River water and ocean water. The isotopic composition of this hybrid water is approximately delta18O = 10 per thousand. delta18O and delta2H values for many saline ground waters plot close to the global meteoric water line, which is distinctly different from the local meteoric water line. This suggests a meteoric origin for ground waters that is different from the current isotopic composition of meteoric waters. It is proposed these waters may be late Pleistocene in age and were recharged when the island was submerged below sea level and prior to rebound at the end of the last glaciation.     

Clay aquitards as archives of holocene paleoclimate: delta18O and thermal profiling

Paleoclimatic conditions in the Holocene were reconstructed from a detailed vertical profile of pore water delta(18)O and a series of downhole thermal profiles at a thick, hydrogeologically simple, aquitard research site in the Northern Great Plains of Saskatchewan. Reconstructions were obtained using the theory of one-dimensional diffusive transport and an empirical Bayesian inversion technique. Inversion of the delta(18)O profile shows that input signal consists of a sudden increase of +6 per thousand (corresponding to a temperature increase of about 18 degrees C) at about 12,000 years before present (BP), after which no measurable change in delta(18)O is apparent to present day. This research shows, at this location, that there is no evidence of large amplitude temperature shifts in the Holocene and supports the commonly used assumption in ground water studies of constant atmospheric delta(18)O values throughout the Holocene. Inversion of the temperature profiles suggests the ground surface temperature increased primarily in the last half of the 20th century, with a peak temperature (about 3 degrees C) circa 1990. For both profiles, the ability to resolve historical variability decays rapidly with time. For the temperature profiles, the decay in resolution precluded the development of reliable estimates of climatic conditions prior to about 1950 and, in the case of the delta(18)O profile, it prevented the precise definition of climate changes (e.g., Hypsithermal and Little Ice Age) in the Holocene.     

Forensic isotope analysis to refine a hydrologic conceptual model

Water resources in the arid southwestern United States are frequently the subject of conflict from competing private and public interests. Legal remedies may remove impasses, but the technical analysis of the problem often determines the future success of legal solutions. In Owens Valley, California, the source of water for the Los Angeles Aqueduct (LAA) is flow diverted from the Owens River and its tributaries and ground water from valley aquifers. Future management of ground water delivered to the LAA needs technical support regarding quantity available, interconnection of shallow and confined aquifers, impact on local springs, and rate of recharge. Ground water flow models and ground water composition are tools already in use, but these have large uncertainty for local interpretations. This study conducted targeted sampling of springs and wells to evaluate the hydrologic system to corroborate conceptual and numerical models. The effort included measurement of intrinsic isotopic composition at key locations in the aquifers. The stable isotopic data of boron (delta(11)B), sulfur (delta(34)S), oxygen (delta(18)O), hydrogen (delta D), and tritium ((3)H) supported by basic chemical data provided rules for characterizing the upper and the lower aquifer system, confirmed the interpretation of ground water flow near faults and flow barriers, and detected hydraulic connections between the LAA and the perennial springs at key locations along the unlined reach of the LAA. This study exemplifies the use of forensic isotopic approaches as independent checks on the consistency of interpretations of conceptual models of a ground water system and the numerical hydrologic simulations.     

Isotopic composition (δ18O and δD) of precipitation and groundwater in a semi‐arid,mountainous area (Guadiana Menor basin,Southeast Spain)

We characterize the precipitation and groundwater in a mountainous (peaks slightly above 3000 m a.s.l.), semi‐arid river basin in SE Spain in terms of the isotopes ¹⁸O and ²H. This basin, with an extension of about 7000 km², is an ideal site for such a study because fronts from the Atlantic and the Mediterranean converge here. Much of the land is farmed and irrigated both by groundwater and runoff water collected in reservoirs. A total of approximately 100 water samples from precipitation and 300 from groundwater have been analysed. To sample precipitation we set up a network of 39 stations at different altitudes (800–1700 m a.s.l.), with which we were able to collect the rain and snowfall from 29 separate events between July 2005 and April 2007 and take monthly samples during the periods of maximum recharge of the aquifers. To characterize the groundwater we set up a control network of 43 points (23 springs and 20 wells) to sample every 3 months the main aquifers and both the thermal and non‐thermal groundwater. We also sampled two shallow‐water sites (a reservoir and a river). The isotope composition of the precipitation forms a local meteoric water line (LMWL) characterized by the equation δD = 7·72δ¹⁸O + 9·90, with mean values for δ¹⁸O and δD of − 10·28‰ and − 69·33‰, respectively, and 12·9‰ for the d‐excess value. To correlate the isotope composition of the rainfall water with groundwater we calculated the weighted local meteoric water line (WLMWL), characterized by the equation δD = 7·40δ¹⁸O + 7·24, which takes into account the quantity of water precipitated during each event. These values of (dδD/dδ¹⁸O)< 8 and d‐excess (δD–8δ¹⁸O)< 10 in each curve bear witness to the ‘amount effect’, an effect which is more manifest between May and September, when the ground temperature is higher. Other effects noted in the basin were those of altitude and the continental influence. The isotopic compositions of the groundwater are represented by the equation δD = 4·79δ¹⁸O − 18·64. The groundwater is richer in heavy isotopes than the rainfall, with mean values of − 8·48‰ for δ¹⁸O and − 59·27‰ for δD. The isotope enrichment processes detected include a higher rate of evaporation from detrital aquifers than from carbonate ones, the effects of recharging aquifers from irrigation return flow and/or from reservoirs' leakage and enrichment in δ¹⁸O from thermal water. Copyright © 2010 John Wiley & Sons, Ltd.     

Advantages and challenges of using soil water isotopes to assess groundwater recharge dominated by snowmelt at a field study located in Canada

Seasonal signals of stable isotopes in precipitation, combined with measurements of isotope ratios in soil water, can be used for quantitative estimation of groundwater recharge rates. This study investigates the applicability of using the piston flow principle and the peak shift displacement method to estimate actual groundwater recharge rates in a humid Nordic region located in the province of Quebec, Canada. Two different sites with and without vegetation (C1 and C2) in an unconfined aquifer were tested by measuring soil water isotope ratios (¹⁸O/¹⁶O and ²H/¹H) and volumetric pore water content. Core samples were obtained along the vadose zone down to the groundwater table at the two sites (2.45 m for Site C1 and 4.15 m for Site C2). The peak shift method to estimate groundwater recharge rates was shown to be accurate only in certain specific conditions inherent to the soil properties and the topographical situation of the investigated sites. Indeed, at Site C2, recharge from the snowmelt could not be estimated because of heterogeneity in the lower part of the vadose zone. At this same site the later recharge after the snowmelt (in the period from late spring to early autumn) could be estimated accurately because the upper part of the vadose zone was homogeneous. Furthermore, at site C1, runoff/runon phenomena hampered calculations of actual infiltration and thus produced inaccurate results for recharge. These two different site effects (heterogeneity in the first site and runoff/runon in the other site) were identified as being limiting factors in the accurate assessment of actual recharge. This study therefore recommends the use of the peak shift method for (1) humid Nordic regions, (2) homogeneous and thick vadose zones, and (3) areas with few or limited site effects (runoff/runon).     

Extinction depth and evapotranspiration from ground water under selected land covers

In many landscapes, vegetation extracts water from both the unsaturated and the saturated zones. The partitioning of evapotranspiration (ET) into vadose zone evapotranspiration and ground water evapotranspiration (GWET) is complex because it depends on land cover and subsurface characteristics. Traditionally, the GWET fraction is assumed to decay with increasing depth to the water table (DTWT), attaining a value of 0 at what is termed the extinction depth. A simple assumption of linear decay with depth is often used but has never been rigorously examined using unsaturated-saturated flow simulations. Furthermore, it is not well understood how to relate extinction depths to characteristics of land cover and soil texture. In this work, variable saturation flow theory is used to simulate GWET for three land covers and a range of soil properties under drying soil conditions. For a water table within half a meter of the land surface, nearly all ET is extracted from ground water due to the close hydraulic connection between the unsaturated and the saturated zones. For deep-rooted vegetation, the decoupling of ground water and vadose zone was found to begin at water table depths between 30 and 100 cm, depending on the soil texture. The decline of ET with DTWT is better simulated by an exponential decay function than the commonly used linear decay. A comparison with field data is consistent with the findings of this study. Tables are provided to vary the extinction depth for heterogeneous landscapes with different vegetation cover and soil properties.     

Isotope methods for management of shared aquifers in northern Africa

Access to fresh water is one of the major issues of northern and sub-Saharan Africa. The majority of the fresh water used for drinking and irrigation is obtained from large ground water basins where there is minor contemporary recharge and the aquifers cross national borders. These aquifers include the Nubian Aquifer System shared by Chad, Egypt, Libya, and Sudan; the Iullemeden Aquifer System, extending over Niger, Nigeria, Mali, Benin, and Algeria; and the Northwest Sahara Aquifer System shared by Algeria, Libya, and Tunisia. These resources are subject to increased exploitation and may be severely stressed if not managed properly as witnessed already by declining water levels. In order to make appropriate decisions for the sustainable management of these shared water resources, planners and managers in different countries need an improved knowledge base of hydrological information. Three technical cooperation projects related to aquifer systems will be implemented by the International Atomic Energy Agency, in collaboration with the United Nations Educational, Scientific and Cultural Organization and United Nations Development Programme/Global Environmental Facility. These projects focus on isotope hydrology studies to better quantify ground water recharge and dynamics. The multiple isotope approach combining commonly used isotopes 18O and 2H together with more recently developed techniques (chlorofluorocarbons, 36Cl, noble gases) will be applied to improve the conceptual model to study stratification and ground water flows. Moreover, the isotopes will be an important indicator of changes in the aquifer due to water abstraction, and therefore they will assist in the effort to establish a sustainable ground water management.     

Constraining a Flow Model with Field Measurements to Assess Water Transit Time Through a Vadose Zone

The modeling of thick vadose zones is particularly challenging because of difficulties in collecting a variety of measured sediment properties, which are required for parameterizing the model. Some models rely on synthetic data, whereas others are simplified by running as homogeneous sediment domains and relying on a single set of sediment properties. Few studies have simulated flow processes through a thick vadose zone using real and comprehensive data sets comprising multiple measurements. Here, we develop a flow model for a 7-m-thick vadose zone. This model, combining the numerical codes CTRAN/W with SEEP/W, includes the measured sediment hydraulic properties of the investigated vadose zone and incorporates the actual climate and subsurface conditions of the study site (precipitations, water-table elevations, and stable isotope data). The model is calibrated by fitting the simulated and measured vertical profiles of water content. Our flow model calculates a transit time of 1 year for the travel of water through the 7-m vadose zone; this estimate matches stable isotope-based results obtained previously for this site. A homogeneous sediment domain flow model, which considers only a single set of sediment properties, produces a transit time that is approximately half the duration of that of the heterogeneous flow model. This difference highlights the importance of assuming heterogeneous material within models of thick vadose zones and testifies to the advantage gained when using real sediment hydraulic properties to parametrize a flow model.     

Application of VLEACH to Vadose Zone Transport of VOCs at an Arizona Superfund Site

Cleanup standards for volatile organic compounds in thick vadose zones can be based on indirect risk (transport to ground water) when contamination is below depths of significant direct risk. At one Arizona Superfund site, a one-dimensional vadose zone transport model (VLE-ACH) was used to estimate the continued transport of VOCs from the vadose zone to ground water. VLEACH is a relatively simple and readily available model that proved useful for estimating indirect risk from VOCs in the vadose zone at this site. The estimates of total soil concentrations used as initial conditions for VLF.ACH incorporated a variety of data from the site. Soil gas concentrations were found to be more useful than soil matrix data for estimating total soil concentrations at this arid-zone site. A simple mixing cell model was used with the VLEACH-derived mass loading estimates from the vadose zone over time to estimate the resulting changes in ground water concentrations. For this site, the results of the linked VLEACH/mixing cell simulations indicate it is likely that the federal MCI. for TCE will be exceeded in underlying ground water if remedial action on I he vadose zone is not pursued.     

Development of saline ground water through transpiration of sea water

As vegetation usually excludes salt during water uptake, transpiration will increase the salinity of the residual water. If the source water is sea water, then the residual water may become highly saline. In the unconfined coastal aquifer of the tropical Burdekin River delta, northeastern Australia, areas of highly saline ground water with chloride concentrations up to almost three times that of sea water occur up to 15 km from the present coastline, and are attributed to transpiration by mangrove vegetation during periods of high sea level. Radiogenic ((14)C) carbon isotope analyses indicate that ground water with chloride concentrations between 15,000 and 35,000 mg/L is mostly between 4000 and 6000 years old, at which time sea level was 2 to 3 m higher than present. Stable isotope analyses of oxygen-18 and deuterium show no evidence for evaporative enrichment of this water. Oxygen-18, deuterium, and stable (delta(13)C) carbon isotope analyses of ground water and soil water point to a recharge environment beneath the mangrove forests during this postglacial sea level high stand. During that period, transpiration of the mangrove forests would have led to high chloride concentrations in the residual ground water, without inducing isotopic fractionation. Due to the higher density, this hypersaline water moved downward through the aquifer by gravity and has formed lenses of highly saline ground water at the bottom of the unconfined aquifer.     

Influence of the Danube River inputs on C and N stable isotope ratios of the Romanian coastal waters and sediment (Black Sea)

The main aim of this study was to estimate the influence of the seasonal variations of the Danube River particulate organic matter (POM) inputs on the Black Sea surface seawater POM and upper layer of sediments along the Romanian coast. Ratios of carbon ((13)C/(12)C) and nitrogen ((15)N/(14)N) stable isotopes allowed differentiating river and marine organic matter sources. Danube River POM presented significantly lower average values of delta(13)C (-27.52+/-0.88 per thousand) and delta(15)N (4.88+/-1.45 per thousand) than seawater POM (delta(13)C=-24.70+/-2.37 per thousand and delta(15)N=6.75+/-1.96 per thousand), whereas surface sediment presented average values similar to seawater POM (delta(13)C=-24.02+/-2.39 per thousand and delta(15)N=7.29+/-2.16 per thousand). Stable isotope values showed that the Danube River influence on marine ecosystems decreased from the North to the South of the Romanian coast. Strong seasonal variations of C and N isotopic signatures were observed in all compartments studied with generally higher values in spring when the river was flooding.     

Strontium isotopic identification of water-rock interaction and ground water mixing

87Sr/86Sr ratios of ground waters in the Bighorn and Laramie basins' carbonate and carbonate-cemented aquifer systems, Wyoming, United States, reflect the distinctive strontium isotope signatures of the minerals in their respective aquifers. Well water samples from the Madison Aquifer (Bighorn Basin) have strontium isotopic ratios that match their carbonate host rocks. Casper Aquifer ground waters (Laramie Basin) have strontium isotopic ratios that differ from the bulk host rock; however, stepwise leaching of Casper Sandstone indicates that most of the strontium in Casper Aquifer ground waters is acquired from preferential dissolution of carbonate cement. Strontium isotope data from both Bighorn and Laramie basins, along with dye tracing experiments in the Bighorn Basin and tritium data from the Laramie Basin, suggest that waters in carbonate or carbonate-cemented aquifers acquire their strontium isotope composition very quickly--on the order of decades. Strontium isotopes were also used successfully to verify previously identified mixed Redbeds-Casper ground waters in the Laramie Basin. The strontium isotopic compositions of ground waters near Precambrian outcrops also suggest previously unrecognized mixing between Casper and Precambrian aquifers. These results demonstrate the utility of strontium isotopic ratio data in identifying ground water sources and aquifer interactions.     

Carbon biogeochemistry of ground water, Guiyang, southwest China

Variations in the concentrations and isotopic compositions (delta13C(DIC)) of dissolved inorganic carbon (DIC) reflect contamination and biogeochemical cycling of the carbon in ground water. In order to understand contamination and biogeochemical cycling of DIC, we carried out research on the geochemistry of ground water of Guiyang, the capital city of Guizhou Province, China. Results show that ground water is mainly characterized by SO4.HCO3-Ca.Mg and HCO3-Ca.Mg chemical compositions. The hydrochemical characteristics of these types of water are mainly controlled by lithology of the aquifers. HCO3- is the dominant species of DIC in ground water and has lower concentrations and more negative values of delta13C(DIC) in the high-flow (summer monsoon) season, as compared to the low-flow season. This indicates that DIC is relatively enriched in carbon of biological origin in the high-flow season as compared to the low-flow season and that biological activities are the predominant control on shifts of stable carbon isotope values. The evidence that the delta13C(DIC) values of ground water decrease with increasing concentrations of anthropogenic species shows that the carbon isotopic composition of DIC can be a useful tracer of contamination, in addition to biogeochemical cycling of inorganic carbon in ground water. Results from this study show that ground water is impacted by significant levels of contamination from human activities, especially in the urban areas, as well as the northeast and west suburbs, in Guiyang city, southwest China.     

An analytical solution for predicting the transient seepage from a subsurface drainage system

Subsurface drainage systems have been widely used to deal with soil salinization and waterlogging problems around the world. In this paper, a mathematical model was introduced to quantify the transient behavior of the groundwater table and the seepage from a subsurface drainage system. Based on the assumption of a hydrostatic pressure distribution, the model considered the pore-water flow in both the phreatic and vadose soil zones. An approximate analytical solution for the model was derived to quantify the drainage of soils which were initially water-saturated. The analytical solution was validated against laboratory experiments and a 2-D Richards equation-based model, and found to predict well the transient water seepage from the subsurface drainage system. A saturated flow-based model was also tested and found to over-predict the time required for drainage and the total water seepage by nearly one order of magnitude, in comparison with the experimental results and the present analytical solution. During drainage, a vadose zone with a significant water storage capacity developed above the phreatic surface. A considerable amount of water still remained in the vadose zone at the steady state with the water table situated at the drain bottom. Sensitivity analyses demonstrated that effects of the vadose zone were intensified with an increased thickness of capillary fringe, capillary rise and/or burying depth of drains, in terms of the required drainage time and total water seepage. The analytical solution provides guidance for assessing the capillary effects on the effectiveness and efficiency of subsurface drainage systems for combating soil salinization and waterlogging problems.