Ground Water Monitoring at Wellhead Protection Areas

Wellhead protection programs implemented to this date have not, in general, included a comprehensive ground water monitoring strategy. This is due both to the relatively recent introduction of the wellhead protection concept into the ground water field and also to the complexity of the monitoring problem. A well-planned monitoring program, however, can contribute significantly to the success of a wellhead protection program. The elements of wellhead protection, important from a monitoring perspective, are the management objectives, the types of potential contaminant sources, and the locations of those sources with respect to the wellhead protection area boundary. General monitoring strategies applicable to wellhead protection are derived from the relationships between these elements.     

The capture efficiency map: the capture zone under time-varying flow

The capture zone or contributing area of a ground water extraction well can be defined as that portion of the aquifer from which the well draws its water. Accurate delineation of capture zones is important in many ground water remediation applications and in the definition of wellhead protection areas. Their mathematical delineation is often simplified by using quasi-steady-state models based on time-weighted average pumping rates and background hydraulic gradients. We present a new semianalytic approach for the definition of capture zones under transient-flow conditions. We then use this approach to evaluate the effects of time variations in the direction of the background hydraulic gradient on capture. Results are presented in the form of capture efficiency maps (CEMs). Although the area contributing to a given well is found to generally expand relative to the steady-state average capture zone when the gradient direction varies, the zone of 100% capture may expand or contract depending on site-specific conditions. We illustrate our CEM approach by applying it to the design of a plume containment system.     

Reducing capture zone uncertainty with a systematic sensitivity analysis

The U.S. Environmental Protection Agency has established several methods to delineate wellhead protection areas (WHPAs) around community wells in order to protect them from surface contamination sources. Delineating a WHPA often requires defining the capture zone for a well. Generally, analytical models or arbitrary setback zones have been used to define the capture zone in areas where little is known about the distribution of hydraulic head, hydraulic conductivity, or recharge. Numerical modeling, however, even in areas of sparse data, offers distinct advantages over the more simplified analytical models or arbitrary setback zones. The systematic approach discussed here calibrates a numerical flow model to regional topography and then applies a matrix of plausible recharge to hydraulic conductivity ratios ( R / K ) to investigate the impact on the size and shape of the capture zone. This approach does not attempt to determine the uncertainty of the model but instead yields several possible capture zones, the composite of which is likely to contain the actual capture zone. A WHPA based on this composite capture zone will protect ground water resources better than one based on any individual capture zone. An application of the method to three communities illustrates development of the R / K matrix and demonstrates that the method is particularly well suited for determining capture zones in alluvial aquifers.     

Probabilistic sensitivity and modeling of two-dimensional transport in porous media

A reliability approach is used to develop a probabilistic model of two-dimensional non-reactive and reactive contaminant transport in porous media. The reliability approach provides two important quantitative results: an estimate of the probability that contaminant concentration is exceeded at some location and time, and measures of the sensitivity of the probabilistic outcome to likely changes in the uncertain variables. The method requires that each uncertain variable be assigned at least a mean and variance; in this work we also incorporate and investigate the influence of marginal probability distributions. Uncertain variables includex andy components of average groundwater flow velocity,x andy components of dispersivity, diffusion coefficient, distribution coefficient, porosity and bulk density. The objective is to examine the relative importance of each uncertain variable, the marginal distribution assigned to each variable, and possible correlation between the variables. Results utilizing a two-dimensional analytical solution indicate that the probabilistic outcome is generally very sensitive to likely changes in the uncertain flow velocity. Uncertainty associated with dispersivity and diffusion coefficient is often not a significant issue with respect to the probabilistic analysis; therefore, dispersivity and diffusion coefficient can often be treated for practical analysis as deterministic constants. The probabilistic outcome is sensitive to the uncertainty of the reaction terms for early times in the flow event. At later times, when source contaminants are released at constant rate throughout the study period, the probabilistic outcome may not be sensitive to changes in the reaction terms. These results, although limited at present by assumptions and conceptual restrictions inherent to the closed-form analytical solution, provide insight into the critical issues to consider in a probabilistic analysis of contaminant transport. Such information concerning the most important uncertain parameters can be used to guide field and laboratory investigations.     

Probabilistic sensitivity and modeling of two-dimensional transport in porous media

A reliability approach is used to develop a probabilistic model of two-dimensional non-reactive and reactive contaminant transport in porous media. The reliability approach provides two important quantitative results: an estimate of the probability that contaminant concentration is exceeded at some location and time, and measures of the sensitivity of the probabilistic outcome to likely changes in the uncertain variables. The method requires that each uncertain variable be assigned at least a mean and variance; in this work we also incorporate and investigate the influence of marginal probability distributions. Uncertain variables includex andy components of average groundwater flow velocity,x andy components of dispersivity, diffusion coefficient, distribution coefficient, porosity and bulk density. The objective is to examine the relative importance of each uncertain variable, the marginal distribution assigned to each variable, and possible correlation between the variables. Results utilizing a two-dimensional analytical solution indicate that the probabilistic outcome is generally very sensitive to likely changes in the uncertain flow velocity. Uncertainty associated with dispersivity and diffusion coefficient is often not a significant issue with respect to the probabilistic analysis; therefore, dispersivity and diffusion coefficient can often be treated for practical analysis as deterministic constants. The probabilistic outcome is sensitive to the uncertainty of the reaction terms for early times in the flow event. At later times, when source contaminants are released at constant rate throughout the study period, the probabilistic outcome may not be sensitive to changes in the reaction terms. These results, although limited at present by assumptions and conceptual restrictions inherent to the closed-form analytical solution, provide insight into the critical issues to consider in a probabilistic analysis of contaminant transport. Such information concerning the most important uncertain parameters can be used to guide field and laboratory investigations.     

Radiocarbon and δ13C Values Related to Ground-Water Recharge and Mixing

Ground-water levels in the Upper Floridan aquifer beneath the southeastern coast of South Carolina have undergone pumpage-induced declines approaching 20 ft below sea level at the southern end of Hilton Head Island. This scenario suggests the potential exists for the inducement of recharge to the Upper Floridan aquifer across the island, which could affect the quality of water being pumped by wells. However, low radiocarbon concentrations in ground-water samples (0.5 to 1.4 ± 0.1 PMC) indicate that most of the water is relict ground water reflecting prepumpage ground-water flow conditions in the Upper Floridan aquifer. The isotopic data indicate long residence times and water-chemistry evolution more characteristic of ground-water recharge occurring farther inland prior to the commencement of pumpage in the late 1800s. Radiocarbon concentrations (as Percent Modern Carbon) and stable carbon isotope ratios (as δ¹³C in dissolved inorganic carbon) determined during this study and reported in other studies on and around Hilton Head Island varied in a systematic manner. Heavier δ¹³C values (–2.8 to –1.6 per mil) in ground water beneath southern Hilton Head Island reflect ground-water discharge from prepumpage flowpaths originating over 100 miles away, hence a depletion in radiocarbon concentration with corrected ground-water ages no younger than 16,000 yrs BP. In contrast, lighter δ¹³C values (–13.9 to –8.67 per mil) beneath the northern part of the island indicate recent recharge as a result of water-level declines, and recharge in areas off the island that have not changed as a result of pumpage (evidenced by enrichment in radiocarbon with corrected ground-water ages no older than 4,000 yrs BP). This suggests that the δ¹³C composition of ground water in the Upper Floridan aquifer is a useful indicator of mixing between ground waters from different sources, and can be used to delineate recharge-discharge patterns. This approach may be applicable to other aquifers of highly evolved ground-water chemistry in regional carbonate aquifer systems that may be receiving recent recharge. Moreover, this approach could prove useful in delineating the contribution of recent water being captured by pumped wells as part of wellhead protection programs designed to assess aquifer vulnerability from surficial contaminant sources.     

Estimating Persistent Mass Flux of Volatile Contaminants from the Vadose Zone to Ground Water

Contaminants may persist for long time periods within low permeability portions of the vadose zone where they cannot be effectively treated and are a potential continuing source of contamination to ground water. Setting appropriate vadose zone remediation goals typically requires evaluating these persistent sources in terms of their impact on meeting ground water remediation goals. Estimating the impact on ground water can be challenging at sites with low aqueous recharge rates where vapor-phase movement is the dominant transport process in the vadose zone. Existing one-dimensional approaches for simulating transport of volatile contaminants in the vadose zone are considered and compared to a new flux-continuity-based assessment of vapor-phase contaminant movement from the vadose zone to the ground water. The flux-continuity-based assessment demonstrates that the ability of the ground water to move contaminant away from the water table controls the vapor-phase mass flux from the vadose zone across the water table. Limitations of these approaches are then discussed with respect to the required assumptions and the need to incorporate three-dimensional processes when evaluating vapor-phase transport from the vadose zone to the ground water. The carbon tetrachloride plume at the U.S. Department of Energy Hanford Site is used as the example site where persistent vadose zone contamination needs to be considered in the context of ground water remediation.     

A Coupled Groundwater–Surface Water Modeling Framework for Simulating Transition Zone Processes

There is an identified need for fully representing groundwater–surface water transition zone (i.e., the sediment zone that connects groundwater and surface water) processes in modeling fate and transport of contaminants to assist with management of contaminated sediments. Most existing groundwater and surface water fate and transport models are not dynamically linked and do not consider transition zone processes such as bioturbation and deposition and erosion of sediments. An interface module is developed herein to holistically simulate the fate and transport by coupling two commonly used models, Environmental Fluid Dynamics Code (EFDC) and SEAWAT, to simulate surface water and groundwater hydrodynamics, while providing an enhanced representation of the processes in the transition zone. Transition zone and surface water contaminant processes were represented through an enhanced version of the EFDC model, AQFATE. AQFATE also includes SEDZLJ, a state‐of‐the‐science surface water sediment transport model. The modeling framework was tested on a published test problem and applied to evaluate field‐scale two‐ and three‐dimensional contaminant transport. The model accurately simulated concentrations of salinity from a published test case. For the field‐scale applications, the model showed excellent mass balance closure for the transition zone and provided accurate simulations of all transition zone processes represented in the modeling framework. The model predictions for the two‐dimensional field case were consistent with site‐specific observations of contaminant migration. This modeling framework represents advancement in the simulation of transition zone processes and can help inform risk assessment at sites where contaminant sources from upland areas have the potential to impact sediments and surface water.     

Monitoring strategies at phreatic wellfields: a 3D travel time approach

Ground water quality networks for monitoring phreatic drinking water wellfields are generally established for two main purposes: (1) the short-term safeguarding of public water supply and (2) signaling and predicting future quality changes in the extracted ground water. Six monitoring configurations with different well locations and different screen depths and lengths were evaluated using a numerical model of the 3D ground water flow toward a partially penetrating pumping well in a phreatic aquifer. Travel times and breakthrough curves for observation and pumping wells were used to judge the effectiveness of different design configurations for three monitoring objectives: (1) early warning; (2) prediction of future quality changes; and (3) evaluation of protection measures inside a protection zone. Effectiveness was tested for scenarios with advective transport, first-order degradation, and linear sorption. It is shown that the location and especially the depth of the observation wells should be carefully chosen, taking into account the residence time from the surface to the observation well, the residual transit times to the extraction well, and the transformation and retardation rates. Shallow monitoring was most functional for a variety of objectives and conditions. The larger the degradation rates or retardation, the shallower should the monitoring be for effective early warning and prediction of future ground water quality. The general approach followed in the current study is applicable for many geohydrological situations, tuning specific monitoring objectives with residence times and residual transit times obtained from a site-specific ground water flow model.     

A Graphical Approach for Determining Dilution-Attenuation Factors: Basic Theory and Approach for Submerged Sources

The dilution attenuation factor (DAF) is a quantity used to relate the concentration of leachate leaving a source zone (e.g., landfill, impoundment, or contaminated soils) to its impact on down-gradient ground water quality. The DAF is of importance because it plays a key role in U.S. Environmental Protection Agency (EPA)'s methodologies for developing soil cleanup goals and for managing hazardous wastes. In this work, a simplistic graphically-based approach for determining site-specific and generic DAFs was developed. In this case the DAF is based on time-and vertically-averaged concentrations along the plume centerline, and the mathematical framework employs well-known analytical and semianalytical solutions for dissolved contaminant transport. Finite sources with a range of decay characteristics are allowed for. One unique feature of this work is that the graphical approach allows for varying levels of sitespecificity, and thus can be used when one has a little, or a lot, of site-specific information. The graphs visually indicate the sensitivity to various parameters, which is valuable information not easily gleaned from most numerical software simulators. This approach is, however, not applicable to very complex hydrogeologic settings (e.g., fractured geology), or to ground water flows that cannot be reasonably approximated as one-dimensional.     

Correlation between nitrate contamination and ground water pollution potential

AQUIPRO, a PC-based method, was used to assess aquifer vulnerability using digital water well logs. The AQUIPRO model is a parameter/factor weighting system for rating the pollution potential of an aquifer. This method uses the well depth, as well as the clay and partial clay thickness in a well, to generate pollution potential scores. In this model, aquifer protection increases as the AQUIPRO vulnerability scores increase and ground water pollution potential decreases. Computerized water well records of 2435 domestic wells with partial chemistry data were used to determine the ground water pollution potential of Kalamazoo County, Michigan. Theoretically, low AQUIPRO pollution potential scores should have more frequent occurrences of ground water contamination events than areas with high AQUIPRO scores with similar land-use, well construction, and well densities. The relative AQUIPRO scores were compared with the frequency of occurrences of nitrate-N in ground water wells. The average nitrate-N concentrations within each relative AQUIPRO vulnerability scores category were also compared. The results indicate that domestic wells containing 5 mg/L or more nitrate-N showed a positive correlation between the frequency of occurrences of nitrate-N and relative decrease of AQUIPRO (r2 = 0.99) vulnerability scores. In other words, as the ground water pollution potential increases, the occurrence frequency of nitrate-N also increases. Furthermore, the results show that as the relative AQUIPRO (r2 = 0.96) vulnerability scores decrease, the mean nitrate-N concentrations also increase.     

Transport and fate of nonaqueous phase liquid (NAPL) in variably saturated porous media with evolving scales of heterogeneity

 A stochastic simulation is performed to study multiphase flow and contaminant transport in fractal porous media with evolving scales of heterogeneity. Numerical simulations of residual NAPL mass transfer and subsequent transport of dissolved and/or volatilized NAPL mass in variably saturated media are carried out in conjunction with Monte Carlo techniques. The impact of fractal dimension, plume scale and anisotropy (stratification) of fractal media on relative dispersivities is investigated and discussed. The results indicate the significance of evolving scale of porous media heterogeneity to the NAPL transport in the subsurface. In general, the fractal porous media enhance the dispersivities of NAPL mass plume transport in both the water phase and the gas phase while the influence on the water phase is more significant. The porous media with larger fractal dimension have larger relative dispersivities. The aqueous horizontal dispersivity exhibits a most significant increase against the plume scale.     

Migration of induced-infiltrated stream water into nearby aquifers due to seasonal ground water withdrawal

Analysis of stream-aquifer interaction due to ground water extraction has traditionally focused on the determination of the amount of water depleted in the stream. Less attention has been paid to the movement of infiltrated stream water inside aquifer, particularly for agricultural areas. This paper presents a method of using particle-tracking techniques to evaluate the transport of the leaked stream water in the nearby aquifers. Simple stream-aquifer conditions are used to demonstrate the usefulness of the analysis. Travel times, pathlines, and influence zones of stream water were determined between a stream and nearby pumping wells for seasonal ground water extraction areas. When water quantity is a concern, the analyses provide additional information about stream depletion; when water quality is an issue, they offer information for wellhead protection. Analyses were conducted for transient conditions, and both pumping and nonpumping periods were considered. According to the results from the simulation examples, migration of infiltrated stream water into the nearby aquifers is generally slow and most infiltrated stream water does not arrive at the pumping well at the end of a 90-day irrigation season. Infiltrated stream water may remain in the aquifer for several years before arriving at the pumping well. For aquifers with a regional hydraulic gradient toward streams, part of the infiltrated stream water may discharge back to streams during a recovery period.     

Comparison of Line-Drive and Push-Pull Flushing Schemes

The performance of cyclodextrin (CD)‐enhanced push‐pull (PP) and line‐drive (LD) approaches to remediation of a site contaminated with a multicomponent dense nonaqueous phase liquid (DNAPL) present in a surficial sandy aquifer was evaluated in this field study. The treatment techniques were compared to each other and to the projected performance of a conventional water‐flushing system. Performance was assessed based on contaminant mass removed per unit volume of extraction solution and per unit time of operation. As expected, the CD‐enhanced LD and PP approaches to remediation were more efficient than conventional flushing with water. Between the two techniques, the PP approach performed 1.5 to 2 times better than the LD approach, particularly for higher DNAPL saturation of the source zone. This result suggests that forcing the flushing solution directly into and through the DNAPL source zone minimized flow bypassing and consequently resulted in a more efficient transfer of contaminant mass between the DNAPL phase and the flushing solution. Nonuniform treatment zone contaminant concentrations and changes in contaminant composition influenced the treatment performances, but these effects were small and still permitted the comparison of successive tests. Although CD was used as the solubility‐enhancing flushing agent in this study, it is likely that the results can be transferred to other chemically enhanced flushing technologies that use, for example, surfactants or alcohols.     

Dual states estimation of a subsurface flow-transport coupled model using ensemble Kalman filtering

Modeling the spread of subsurface contaminants requires coupling a groundwater flow model with a contaminant transport model. Such coupling may provide accurate estimates of future subsurface hydrologic states if essential flow and contaminant data are assimilated in the model. Assuming perfect flow, an ensemble Kalman filter (EnKF) can be used for direct data assimilation into the transport model. This is, however, a crude assumption as flow models can be subject to many sources of uncertainty. If the flow is not accurately simulated, contaminant predictions will likely be inaccurate even after successive Kalman updates of the contaminant model with the data. The problem is better handled when both flow and contaminant states are concurrently estimated using the traditional joint state augmentation approach. In this paper, we introduce a dual estimation strategy for data assimilation into a one-way coupled system by treating the flow and the contaminant models separately while intertwining a pair of distinct EnKFs, one for each model. The presented strategy only deals with the estimation of state variables but it can also be used for state and parameter estimation problems. This EnKF-based dual state-state estimation procedure presents a number of novel features: (i) it allows for simultaneous estimation of both flow and contaminant states in parallel; (ii) it provides a time consistent sequential updating scheme between the two models (first flow, then transport); (iii) it simplifies the implementation of the filtering system; and (iv) it yields more stable and accurate solutions than does the standard joint approach. We conducted synthetic numerical experiments based on various time stepping and observation strategies to evaluate the dual EnKF approach and compare its performance with the joint state augmentation approach. Experimental results show that on average, the dual strategy could reduce the estimation error of the coupled states by 15% compared with the joint approach. Furthermore, the dual estimation is proven to be very effective computationally, recovering accurate estimates at a reasonable cost.     

Evidence of Spatio-Temporal Variations in Contaminants Discharging to a Peri-Urban Stream

Xenobiotic organic compounds can be discharged from contaminated groundwater inflow and may seep into streams from multiple pathways with very different dynamics, some not fully understood. In this study, we investigated the spatio-temporal variation of chlorinated ethenes discharging from a former industrial site (with two main contaminant sources, A and B) into a stream system in a heterogeneous clay till setting in eastern Denmark. The investigated reach and near-stream surroundings are representative of peri-urban settings, with a mix of high channel alteration and more natural stream environment. We therefore propose an approach for risk assessing impacts arising from such complex contamination patterns, accounting for potential spatio-temporal fluctuations and presence of multiple pathways. Our study revealed substantial variations in pathway contributions and overall contaminant mass discharge to the stream. Variable contaminant contributions arising from both groundwater seepage and urban drains were identified in the channelized part of the north stream, primarily from source A. Furthermore, variations in the hyporheic and shallow groundwater flows were found to enhance contaminant transport from source B. These processes result in an increase of the overall mass of contaminant discharged, correlating with the channels' flow. Thus, an in-stream control plane approach was found to be an effective method for integrating multiple and variable discharge contributions quantitatively, although information on specific contaminant sources is lost. This study highlights the complexity and variability of contaminant fluxes occurring at the interface between groundwater and peri-urban streams, and calls for the consideration of these variations when designing monitoring programs and remedial actions for contaminated sites with the potential to impact streams.     

Three-Dimensional Simulation of Volatile Organic Compound Mass Flux from the Vadose Zone to Groundwater

Low-permeability layers of the vadose zone containing volatile organic compounds (VOCs) may persist as source zones for long time periods and may provide contamination to groundwater. At sites with low recharge rates, where vapor migration is the dominant transport process, the impact of vadose zone sources on groundwater may be difficult to assess. Typical assessment methods include one-dimensional numerical and analytical techniques. The one-dimensional approaches only consider groundwater coupling options through boundary conditions at the water table and may yield artificially high mass flux results when transport is assumed to occur by gas-phase diffusion between a source and an interface with a zero concentration boundary condition. Improvements in mass flux assessments for VOCs originating from vadose zone sources may be obtained by coupling vadose zone gas transport and dissolved contaminant transport in the saturated zone and by incorporating the inherent three-dimensional nature of gas-phase transport, including the potential of density-driven advection. This paper describes a series of three-dimensional simulations using data from the U.S. Department of Energy's Hanford site, where carbon tetrachloride is present in a low-permeability zone about 30 m above the groundwater. Results show that, for most cases, only a relatively small amount of the contaminant emanating from the source zone partitions into the groundwater and that density-driven advection is only important when relatively high source concentrations are considered.     

Partitioned multiobjective risk modeling of carcinogenic compounds in groundwater

Quantifying human cancer risk arising from exposure to contaminated groundwater is complicated by the many hydrogeological, environmental, and toxicological uncertainties involved. In this study, we used Monte Carlo simulation to estimate cancer risk associated with tetrachloroethene (PCE) dissolved in groundwater by linking three separate models for: (1) reactive contaminant transport; (2) human exposure pathways; and (3) the PCE cancer potency factor. The hydrogeologic model incorporates an analytical solution for a one-dimensional advective–dispersive–reactive transport equation to determine the PCE concentration in a water supply well located at a fixed distance from a continuous source. The pathway model incorporates PCE exposure through ingestion, inhalation, and dermal contact. The toxicological model combines epidemiological data from eight rodent bioassays of PCE exposure in the form of a composite cumulative distribution frequency curve for the human PCE cancer potency factor. We assessed the relative importance of individual model variables through their correlation with expected cancer risk calculated in an ensemble of Monte Carlo simulations with 20,000 trials. For the scenarios evaluated, three factors were most highly correlated with cancer risk: (1) the microbiological decay constant for PCE in groundwater, (2) the linear groundwater pore velocity, and (3) the cancer potency factor. We then extended our analysis beyond conventional expected value risk assessment using the partitioned multiobjective risk method (PMRM) to generate expected-value functions conditional to a 1 in 100,000 increased cancer risk threshold. This approach accounts for low probability/high impact outcomes separately from the conventional unconditional expected values. Thus, information on potential worst-case outcomes can be quantified for decision makers. Using PMRM, we evaluated the cost-benefit relationship of implementing several postulated risk management alternatives intended to mitigate the expected and conditional cancer risk. Our results emphasize the importance of hydrogeologic models in risk assessment, but also illustrate the importance of integrating environmental and toxicological uncertainty. When coupled with the PMRM, models integrating uncertainty in transport, exposure, and potency constitute an effective risk assessment tool for use within a risk-based corrective action (RBCA) framework.     

Multiobjective analysis of a public wellfield using artificial neural networks

As competition for increasingly scarce ground water resources grows, many decision makers may come to rely upon rigorous multiobjective techniques to help identify appropriate and defensible policies, particularly when disparate stakeholder groups are involved. In this study, decision analysis was conducted on a public water supply wellfield to balance water supply needs with well vulnerability to contamination from a nearby ground water contaminant plume. With few alternative water sources, decision makers must balance the conflicting objectives of maximizing water supply volume from noncontaminated wells while minimizing their vulnerability to contamination from the plume. Artificial neural networks (ANNs) were developed with simulation data from a numerical ground water flow model developed for the study area. The ANN-derived state transition equations were embedded into a multiobjective optimization model, from which the Pareto frontier or trade-off curve between water supply and wellfield vulnerability was identified. Relative preference values and power factors were assigned to the three stakeholders, namely the company whose waste contaminated the aquifer, the community supplied by the wells, and the water utility company that owns and operates the wells. A compromise pumping policy that effectively balances the two conflicting objectives in accordance with the preferences of the three stakeholder groups was then identified using various distance-based methods.