A Tracer Test to Characterize Treatment of TCE in a Permeable Reactive Barrier

A tracer test was conducted to characterize the flow of groundwater across a permeable reactive barrier constructed with plant mulch (a biowall) at the OU‐1 site on Altus Air Force Base, Oklahoma. This biowall is intended to intercept and treat groundwater contaminated by trichloroethylene (TCE) in a shallow aquifer. The biowall is 139‐m long, 7.3‐m deep, and 0.5‐m wide. Bromide was injected from an upgradient well into the groundwater as a conservative tracer, and was subsequently observed breaking through in monitoring wells within and downgradient of the biowall. The bromide breakthrough data demonstrate that groundwater entering the biowall migrated across it, following the slope of the local groundwater surface. The average seepage velocity of groundwater was approximately 0.06 m/d. On the basis of the Darcy velocity of groundwater and geometry of the biowall, the average residence time of groundwater in the biowall was estimated at 10 d. Assuming all TCE removal occurred in the biowall, the reduction in TCE concentrations in groundwater across the biowall corresponds to a first‐order attenuation rate constant in the range of 0.38 to 0.15 per d. As an independent estimate of the degradation rate constant, STANMOD software was used to fit curves through data on the breakthrough of bromide and TCE in selected wells downgradient of the injection wells. Best fits to the data required a first‐order degradation rate constant for TCE removal in the range of 0.13 to 0.17 per d. The approach used in this study provides an objective evaluation of the remedial performance of the biowall that can provide a basis for design of other biowalls that are intended to remediate TCE‐contaminated groundwater.     

Permeable Asphalt: A New Tool to Reduce Road Salt Contamination of Groundwater in Urban Areas

Chloride contamination of groundwater in urban areas due to deicing is a well‐documented phenomenon in northern climates. The objective of this study was to evaluate the effects of permeable pavement on degraded urban groundwater. Although low impact development practices have been shown to improve stormwater quality, no infiltration practice has been found to prevent road salt chlorides from entering groundwater. The few studies that have investigated chlorides in permeable asphalt have involved sampling directly beneath the asphalt; no research has looked more broadly at surrounding groundwater conditions. Monitoring wells were installed upgradient and downgradient of an 860 m² permeable asphalt parking lot at the University of Connecticut (Storrs, Connecticut). Water level and specific conductance were measured continuously, and biweekly samples were analyzed for chloride. Samples were also analyzed for sodium (Na), calcium (Ca), and magnesium (Mg). Analysis of variance analysis indicated a significantly (p < 0.001) lower geometric mean Cl concentration downgradient (303.7 mg/L) as compared to upgradient (1280 mg/L). Concentrations of all alkali metals increased upgradient and downgradient during the winter months as compared to nonwinter months, indicating that cation exchange likely occurred. Despite the frequent high peaks of chloride in the winter months as well as the increases in alkali metals observed, monitoring revealed lower Cl concentrations downgradient than upgradient for the majority of the year. These results suggest that the use of permeable asphalt in impacted urban environments with high ambient chloride concentrations can be beneficial to shallow groundwater quality, although these results may not be generalizable to areas with low ambient chloride concentrations.     

Oxygen-Containing Compounds Identified in Groundwater from Fuel Release Sites Using GCxGC-TOF-MS

This research continues a 7-year study of oxygen-containing organic compounds present in groundwater at gasoline and diesel fuel release sites that are quantified as diesel-range “total petroleum hydrocarbons” when measured by methods utilizing solvent extraction and gas chromatography. Two-dimensional gas chromatography with time-of-flight mass spectrometry was used to tentatively identify 1162 compounds (TICs) in 113 groundwater samples from 22 sites. Samples were collected from wells either upgradient of the release, within the source zone, or downgradient of the source but still within the plume of dissolved organics associated with release. The names and formulas of all TICs found in samples from each well type are presented and the results from upgradient and downgradient locations are compared in detail. About 60% of the most frequently detected TICs in downgradient wells were also detected in upgradient wells. A majority of these were saturated straight chain alkyl acids, commonly called fatty acids, or fatty acid esters. Of TICs frequently detected in downgradient wells but not upgradient wells, over half were branched alkyl alcohols. Hierarchical cluster analysis results suggest about 80% of the chemical composition of downgradient samples is more similar to upgradient samples than to source area samples. This similarity is due to the presence of the same types of fatty acids and esters. Principal component analysis indicates a continuum of biodegradation between the source area and downgradient samples with the latter becoming more consistent with upgradient samples. Results suggest some TICs may not be petroleum degradation intermediates but compounds synthesized by microorganisms through secondary production and carbon cycling.     

Macropores as preferential flow paths in meander bends

Macropores are subsurface connected void spaces caused by processes such as fracture of soils, micro‐erosion, and fauna burrows. They are common near streams (e.g. hyporheic and riparian zones) and may act as preferential flow paths between surface and groundwaters, affecting hydrologic and biogeochemical processes. We tested the hydrologic function of macropores by constructing an artificial macropore within the saturated zone of a meander bend (open macropore, ‘OM’) and later filling its upstream end (partially filled macropore, ‘PFM’). For each treatment, we injected saline tracer at an upgradient monitoring well within the meander and monitored downgradient hydraulics and tracer transport. Pressure transducers in monitoring wells indicated hydraulic gradients within the meander were 32% higher perpendicular to and 6% higher parallel to the macropore for the OM than for the PFM. Additionally, hydraulic conductivities measured via falling head tests were 29 to 550 times higher along the macropore than in nearby sediment. We used electrical conductivity probes in wells and electrical resistivity imaging to track solute transport. Transport velocities through the meander were on average 9 and 21% higher (per temporal moment analysis and observed tracer peak, respectively) for the OM than for the PFM. Furthermore, temporal moments of tracer breakthrough analysis indicated downgradient longitudinal dispersion and breakthrough tracer curve tailing were on average 234% and 182% higher for the OM, respectively. This suggests the OM enabled solute transport at overall shorter timescales than the matrix but also increased tailing. Our results demonstrate the importance of macropores to meander bend hydrology and solute transport. Copyright © 2012 John Wiley & Sons, Ltd.     

Determination of longitudinal dispersivity in an unconfined sandy aquifer

Study of the mobility of contaminants in an aquifer is an important issue for the proper remediation of contaminated groundwater. Determination of associated solute transport parameters therefore is essential for investigation of the extent to which groundwater can be contaminated. This study aimed at determining solute transport parameters for an unconfined sandy aquifer at a laboratory scale through various tracer tests using a conservative solute as a tracer. Tracer tests consisted of both well‐tracer tests (single and double wells) and an aquifer tracer test using a plume‐capturing device such as time domain reflectometry (TDR). The results showed that longitudinal dispersivities estimated from the single and double well‐tracer tests were 2·2 cm and 13·5 cm for a travel distance of 9·3 cm and 13·5 cm from the injection point respectively. These results agreed reasonably well with the results of the aquifer tracer test. The solute transport parameters obtained at multiple points in the aquifer through the aquifer tracer test revealed that the dispersivity length was proportional to the travel distance by a factor of 0·3, which was moderately higher than the value of 0·1 given in the literature. Copyright © 2002 John Wiley & Sons, Ltd.     

Single Well Tracer Tests in Aquifer Characterization

When the purpose of aquifer testing is to yield data for modeling aqueous mass transport, pumping tests and gradient measurement can only partially satisfy characterization requirements. Effective porosity, ground water flow velocity, and the vertical distribution of hydraulic conductivity within the aquifer are left as unknowns. Single well tracer methods, when added to the testing program, can be used to estimate these parameters. A drift, and pumpback test yields porosity and velocity, and point-dilution testing yields depth-discrete hydraulic information, A single emplacement of tracer into a test well is sufficient to conduct both tests. The tracer tests are facilitated by a simple method for injecting and evenly distributing the tracer solution into a wellbore, and by new ion-selective electrode instrumentation, specifically designed for submersible service, for monitoring the concentration of tracers such as bromide.     

Geochemical imaging of flow near an artificial recharge facility, Orange County, California

Critical for the management of artificial recharge operations is detailed knowledge of ground water dynamics near spreading areas. Geochemical tracer techniques including stable isotopes of water, tritium/helium-3 (T/3He) dating, and deliberate gas tracer experiments are ideally suited for these investigations. These tracers were used to evaluate flow near an artificial recharge site in northern Orange County, California, where approximately 2.5 x 10(8) m3 (200,000 acre-feet) of water are recharged annually. T/3He ages show that most of the relatively shallow ground water within 3 km of the recharge facilities have apparent ages < 2 years; further downgradient apparent ages increase, reaching > 20 years at approximately 6 km. Gas tracer experiments using sulfur hexafluoride and xenon isotopes were conducted from the Santa Ana River and two spreading basins. These tracers were followed in the ground water for more than two years, allowing subsurface flow patterns and flow times to be quantified. Results demonstrate that mean horizontal ground water velocities range from < 1 to > 4 km/year. The leading edges of the tracer patch moved at velocities about twice as fast as the center of mass. Leading edge velocities are important when considering the potential transport of microbes and other "time sensitive" contaminants and cannot be determined easily with other methods. T/3He apparent ages and tracer travel times agreed within the analytical uncertainty at 16 of 19 narrow screened monitoring wells. By combining these techniques, ground water flow was imaged with time scales on the order of weeks to decades.     

An analytical solution for ground water transit time through unconfined aquifers

An exact, closed-form analytical solution is developed for calculating ground water transit times within Dupuit-type flow systems. The solution applies to steady-state, saturated flow through an unconfined, horizontal aquifer recharged by surface infiltration and discharging to a downgradient fixed-head boundary. The upgradient boundary can represent, using the same equation, a no-flow boundary or a fixed head. The approach is unique for calculating travel times because it makes no a priori assumptions regarding the limit of the water table rise with respect to the minimum saturated aquifer thickness. The computed travel times are verified against a numerical model, and examples are provided, which show that the predicted travel times can be on the order of nine times longer relative to existing analytical solutions.     

Line-source multi-tracer test for assessing high groundwater velocity

Segmented line-source multi-tracer injection is suggested as an effective method for assessing groundwater velocities and flow directions in subsurfaces characterized by high water flux. Modifying the common techniques of injecting a tracer into a well became necessary after point-source natural and forced gradient tracer tests ended with no reliable information on the local groundwater flow. The tracer's line-source increases the likelihood of success of the test and could provide additional information regarding the lateral heterogeneity of the aquifer. In a field experiment conducted in the northwestern part on the Dead Sea coast, tracers were injected into an 8-m-long line injection system perpendicular to the assumed flow direction. The injection system was divided into four separate segments with four different tracers. An array of five boreholes located within a 10 × 10 m area downstream was used for monitoring the tracers' transport. Two dye tracers (uranine and Na naphthionate) were injected in a long pulse of several hours into two of the injection pipe segments. Two other tracers (Rhenium oxide and Gd-DTPA) were instantaneously injected into the other two segments. The tracers were detected 0.7 to 2.3 h after injection in four of the five observation wells, located 2.3 to 10 m away from the injection system. The groundwater velocity was determined to be ～80 to 170 m/d, based on the recoveries of the tracers. The groundwater flow direction was derived based on the arrival of the tracers and was found to be quite consistent with the apparent direction of the hydraulic gradient.     

Backward probabilistic model of groundwater contamination in non-uniform and transient flow

Backward location and travel time probabilities, which provide information about the former location of contamination in an aquifer, can be used to identify unknown contamination sources. Backward location probability describes the possible upgradient positions of contamination at a known time in the past, and backward travel time probability describes the time required for contamination to travel from a known upgradient location to an observation point. These probabilities are related to adjoint states of resident concentration, and their governing equation is the adjoint of a forward contaminant transport model. Using adjoint theory to obtain the appropriate governing equation, we extend the backward probability model for conservative solutes to more general non-uniform and transient flow fields. In particular, we address three important extensions, spatially-varying porosity, transient flow and temporally-varying porosity, and internal distributed sources and sinks of solute and water. For the first time we learn that forward and backward location and travel time probabilities are not necessarily equivalent to adjoint states, but are related to them. The extensions are illustrated using a vertically-integrated groundwater model, creating transient flow by a step change in pumping and using areal recharge as an internal distributed source. Both the movement and spread of probabilities are affected. With internal sources of water, there are two interpretations of backward probability, depending on whether or not the source of water is also a source of solute. The results demonstrate how the backward probability model can be applied to other, perhaps more important, non-uniform and transient flow conditions, with time- and space-varying water storage, such as time-varying pumping or unsaturated (or saturated–unsaturated) flow and transport with spatially- and temporally-varying moisture content.     

Relationship Between Subsurface Landfill Gas and Arsenic Mobilization into Groundwater

The maximum contaminant level for arsenic was reduced by the U.S. Environmental Protection Agency (U.S. E.P.A.) for Drinking Water Standards from 50 micrograms per liter (µg/L) to 10 µg/L, effective January 23, 2006. The subject site is a double-lined sanitary landfill facility located in the mid-Atlantic region of the United States. Arsenic was reported above the maximum contaminant level in a downgradient monitoring well (MW-18) in July 2005. Since July 2005, arsenic levels in MW-18 fluctuated above and below 10 µg/L. This research focuses on determining whether reducing conditions in groundwater, enhanced by subsurface landfill gas emissions, were causing naturally-occurring arsenic to mobilize from the native variably-saturated vadose zone soils into groundwater. The groundwater data collected from the impacted well (MW-18) were compared to an upgradient well (MW-8) to determine whether significant differences existed during the time period of April 2004 to April 2007. Linear regression analysis was also used to determine whether other parameters had a significant relationship with the arsenic concentrations detected in MW-18. The groundwater located in MW-18 was consistently more reduced than groundwater located in the upgradient/background well MW-8, and this was most likely attributed to the presence of subsurface landfill gas in the area. According to the U.S. E.P.A., oxidation-reduction potential (ORP) values in groundwater less than 50 millivolts (mv) suggest that a reducing environment may be present. The data presented in this study indicate that arsenic can mobilize into groundwater under moderately reducing conditions, with ORP measurements averaging 53 mv.     

Simultaneous analyses and applications of multiple fluorobenzoate and halide tracers in hydrologic studies

An analytical method that employs ion chromatography has been developed to exploit the use of fluorobenzoic acids (FBAs) and halides more fully as hydrologic tracers. In a single run, this reliable, sensitive, and robust method can simultaneously separate and quantify halides (fluoride, chloride, bromide, and iodide) and up to seven FBAs from other common groundwater constituents (e.g. nitrate and sulphate). The usefulness of this analytical method is demonstrated in both field and laboratory tracer experiments. The field study examines the hydrologic response of fractures and the matrix to different flow rates and the contribution of matrix diffusion in chemical transport. Laboratory tracer experiments with eight geologic media from across the USA—mostly from Department of Energy facilities where groundwater contamination is prevalent and where subsurface characterization employing tracers has been ongoing or is in need—reveal several insights about tracer transport behaviour: (1) bromide and FBAs are not always transported conservatively; (2) the delayed transport of these anionic tracers is likely related to geologic media characteristics, such as organic matter, pH, iron oxide content, and clay mineralogy; (3) use of iodine as a hydrologic tracer should take into account the different sorption behaviours of iodide and iodate and the possible conversion of iodine's initial chemical form; (4) the transport behaviour of potential FBA and halide tracers under relevant geochemical conditions should be evaluated before beginning ambitious, large‐scale field tracer experiments. Copyright © 2005 John Wiley & Sons, Ltd.     

Identification of Multiple Sources of Groundwater Contamination by Dual Isotopes

Chlorinated solvents are one of the most commonly detected groundwater contaminants in industrial areas. Identification of polluters and allocation of contaminant sources are important concerns in the evaluation of complex subsurface contamination with multiple sources. In recent years, compound‐specific isotope analyses (CSIA) have been employed to discriminate among different contaminant sources and to better understand the fate of contaminants in field‐site studies. In this study, the usefulness of dual isotopes (carbon and chlorine) was shown in assessments of groundwater contamination at an industrial complex in Wonju, Korea, where groundwater contamination with chlorinated solvents such as trichloroethene (TCE) and carbon tetrachloride (CT) was observed. In November 2009, the detected TCE concentrations at the study site ranged between nondetected and 10,066 µg/L, and the CT concentrations ranged between nondetected and 985 µg/L. In the upgradient area, TCE and CT metabolites were detected, whereas only TCE metabolites were detected in the downgradient area. The study revealed the presence of separate small but concentrated TCE pockets in the downgradient area, suggesting the possibility of multiple contaminant sources that created multiple comingling plumes. Furthermore, the variation of the isotopic (δ¹³C and δ³⁷Cl) TCE values between the upgradient and downgradient areas lends support to the idea of multiple contamination sources even in the presence of detectable biodegradation. This case study found it useful to apply a spatial distribution of contaminants coupled with their dual isotopic values for evaluation of the contaminated sites and identification of the presence of multiple sources in the study area.     

Environmental tracer transport (H and SF6) in the saturated and unsaturated zones and its use in nitrate pollution management

An important quantity in groundwater protection is the residence time of water in an aquifer. It relates to both the travel time of a pollutant to arrive at a well and the time span required for self-purification of a polluted aquifer after removal of pollutant inputs. Time scales for aquifers can be gained from artificial tracer experiments or from environmental tracer data, the latter offering the only realistic alternative if time scales of years or decades have to be taken into account.

Different tracers show different time scales due to their different transport mechanisms especially in the unsaturated zone. While solute tracers are moved advectively with the seepage water, gas tracers pass the unsaturated zone diffusively through the air phase. Depending on the properties of the unsaturated zone (hydraulic properties, thickness) this difference in behavior can be used to separate the subsurface transport process into the unsaturated and the saturated parts.

In a field study in Germany, SF₆ and ³H were used as environmental tracers. Both have a relatively well-known input function. Interpretation of data from observation wells by a box model approach led to spatially and temporally varying residence times. This was an indication that the influence of the unsaturated zone could not be neglected. While the gas tracer SF₆ shows only residence times in the saturated zone, the tracer ³H reflects the whole travel time of water including both the unsaturated and saturated zones. Using a one-dimensional plug-flow model for the unsaturated zone combined with a detailed two-dimensional flow and transport model for the saturated zone leads to a holistic and consistent interpretation of the measured tracer concentrations. The observed pattern of old water under thick loess cover and younger water under areas where the fractured basalt aquifer crops out is reproduced after adjusting only two parameters: the effective porosity of the saturated aquifer and the product of field capacity and thickness of the unsaturated zone. While the effective porosity of the saturated zone is adjusted by means of the SF₆ data, the field capacity of the loess layer is adjusted by means of the ³H observations. The thickness of the unsaturated zone is deduced from geological and pedological maps. All flow data are obtained from a calibrated flow model, which is based on geological data, observed heads and pumping tests only.

The transport model for the saturated zone was calibrated by fitting the porosity by means of gaseous tracer concentrations (SF₆). The combined saturated–unsaturated zone model was then calibrated by fitting the field capacity of the unsaturated zone by means of ³H concentrations. With this model it was possible to verify the observed NO₃ concentrations at the drinking water wells and to develop predictions for their future development under various scenarios of fertilizer input reduction in specific areas.     

Viruses as Groundwater Tracers: Using Ecohydrology to Characterize Short Travel Times in Aquifers

Viruses are attractive tracers of short (<3 year) travel times in aquifers because they have unique genetic signatures, are detectable in trace quantities, and are mobile in groundwater. Virus “snaphots” result from infection and disappearance in a population over time; therefore, the virus snapshot shed in the fecal wastes of an infected population at a specific point in time can serve as a marker for tracking virus and groundwater movement. The virus tracing approach and an example application are described to illustrate their ability to characterize travel times in high‐groundwater velocity settings, and provide insight unavailable from standard hydrogeologic approaches. Although characterization of preferential flowpaths does not usually characterize the majority of other travel times occurring in the groundwater system (e.g., center of plume mass; tail of the breakthrough curve), virus approaches can trace very short times of transport, and thus can fill an important gap in our current hydrogeology toolbox.     

Point dilution tests to calculate groundwater velocity: an example in a porous aquifer in northeast Italy

ABSTRACT

The point dilution test is a single-well technique for estimating horizontal flow velocity in the aquifer surrounding a well. The test is conducted by introducing a tracer into a well section and monitoring its decreasing concentration over time. When using a salt tracer, the method is easy and inexpensive. Traditionally, the horizontal Darcy velocity is calculated as a function of the rate of dilution and is based on the simple assumption that the decreasing tracer concentration is proportional both to the apparent velocity into the test section and to the Darcy velocity in the aquifer. In this article, an alternative approach to analyse the results of point dilution tests is proposed and verified using data acquired at a test site in the middle Venetian plain, northeast Italy. In this approach, the one-dimensional equilibrium advection–dispersion equation is inverted using the CXTFIT model to estimate the apparent velocity inside the test section. Analysis of the field data obtained by the two approaches shows good agreement between the methods and suggests that it is possible to use the equilibrium advection–dispersion equation to estimate apparent velocity over a wide range of velocities.

Editor D. Koutsoyiannis; Associate editor K. Heal     

Estimating Reaction Rate Coefficients Within a Travel‐Time Modeling Framework

A generalized, efficient, and practical approach based on the travel‐time modeling framework is developed to estimate in situ reaction rate coefficients for groundwater remediation in heterogeneous aquifers. The required information for this approach can be obtained by conducting tracer tests with injection of a mixture of conservative and reactive tracers and measurements of both breakthrough curves (BTCs). The conservative BTC is used to infer the travel‐time distribution from the injection point to the observation point. For advection‐dominant reactive transport with well‐mixed reactive species and a constant travel‐time distribution, the reactive BTC is obtained by integrating the solutions to advective‐reactive transport over the entire travel‐time distribution, and then is used in optimization to determine the in situ reaction rate coefficients. By directly working on the conservative and reactive BTCs, this approach avoids costly aquifer characterization and improves the estimation for transport in heterogeneous aquifers which may not be sufficiently described by traditional mechanistic transport models with constant transport parameters. Simplified schemes are proposed for reactive transport with zero‐, first‐, nth‐order, and Michaelis‐Menten reactions. The proposed approach is validated by a reactive transport case in a two‐dimensional synthetic heterogeneous aquifer and a field‐scale bioremediation experiment conducted at Oak Ridge, Tennessee. The field application indicates that ethanol degradation for U(VI)‐bioremediation is better approximated by zero‐order reaction kinetics than first‐order reaction kinetics.     

The Horizontal Reactive Media Treatment Well (HRX Well®) for Passive In‐Situ Remediation

A new in‐situ remediation concept termed a Horizontal Reactive Media Treatment Well (HRX Well®) is presented that utilizes horizontal wells filled with reactive media to passively treat contaminated groundwater in‐situ. The approach involves the use of large‐diameter directionally drilled horizontal wells filled with granular reactive media generally installed parallel to the direction of groundwater flow. The design leverages natural “flow‐focusing” behavior induced by the high in‐well hydraulic conductivity of the reactive media relative to the aquifer hydraulic conductivity to passively capture and treat proportionally large volumes of groundwater within the well. Clean groundwater then exits the horizontal well along its downgradient sections. Many different types of solid granular reactive media are already available (e.g., zero valent iron, activated carbon, ion exchange resins, zeolite, apatite, chitin); therefore, this concept could be used to address a wide range of contaminants. Three‐dimensional flow and transport simulations were completed to assess the general hydraulic performance, capture zones, residence times, effects of aquifer heterogeneity, and treatment effectiveness of the concept. The results demonstrate that capture and treatment widths of up to tens of feet can be achieved for many aquifer settings, and that reductions in downgradient concentrations and contaminant mass flux are nearly immediate. For a representative example, the predicted treatment zone width for the HRX Well is approximately 27 to 44 feet, and contaminant concentrations immediately downgradient of the HRX Well decreased an order of magnitude within 10 days. A series of laboratory‐scale physical tests (i.e., tank tests) were completed that further demonstrate the concept and confirm model prediction performance. For example, the breakthrough time, peak concentration and total mass recovery of methylene blue (reactive tracer) was about 2, 35, and 20 times (respectively) less than chloride (conservative tracer) at the outlet of the tank‐scale HRX Well.     

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.     

On Solving Groundwater Flow and Transport Models with Algebraic Multigrid Preconditioning

Iterative solvers preconditioned with algebraic multigrid have been devised as an optimal technology to speed up the response of large sparse linear systems. In this work, this technique was implemented in the framework of the dual delineation approach. This involves a single groundwater flow linear solution and a pure advective transport solution with different right-hand sides. The new solver was compared with other preconditioned iterative methods, the MODFLOW's GMG solver, and direct sparse solvers. Test problems include two- and three-dimensional benchmarks spanning homogeneous and highly heterogeneous and anisotropic formations. For the groundwater flow problems, using the algebraic multigrid preconditioning speeds up the numerical solution by one to two orders of magnitude. The algebraic multigrid preconditioner efficiency was preserved for the three dimensional heterogeneous and anisotropic problem unlike for the MODFLOW's GMG solver. Contrarily, a sparse direct solver was the most efficient for the pure advective transport processes such as the forward travel time simulations. Hence, the best sparse solver for the more general advection-dispersion transport equation is likely to be Péclet number dependent. When equipped with the best solvers, processing multimillion grid blocks by the dual delineation approach is a matter of seconds. This paves the way for its routine application to large geological models. The paper gives practical hints on the strategies and conditions under which algebraic multigrid preconditioning would remain competitive for the class of nonlinear and/or transient problems.