Vadose zone-attenuated artificial recharge for input to a ground water model

Accurate representation of artificial recharge is requisite to calibration of a ground water model of an unconfined aquifer for a semiarid or arid site with a vadose zone that imparts significant attenuation of liquid transmission and substantial anthropogenic liquid discharges. Under such circumstances, artificial recharge occurs in response to liquid disposal to the vadose zone in areas that are small relative to the ground water model domain. Natural recharge, in contrast, is spatially variable and occurs over the entire upper boundary of a typical unconfined ground water model. An improved technique for partitioning artificial recharge from simulated total recharge for inclusion in a ground water model is presented. The improved technique is applied using data from the semiarid Hanford Site. From 1944 until the late 1980s, when Hanford's mission was the production of nuclear materials, the quantities of liquid discharged from production facilities to the ground vastly exceeded natural recharge. Nearly all hydraulic head data available for use in calibrating a ground water model at this site were collected during this period or later, when the aquifer was under the diminishing influence of the massive water disposals. The vadose zone is typically 80 to 90 m thick at the Central Plateau where most production facilities were located at this semiarid site, and its attenuation of liquid transmission to the aquifer can be significant. The new technique is shown to improve the representation of artificial recharge and thereby contribute to improvement in the calibration of a site-wide ground water model.     

Recharge estimation for transient ground water modeling

Reliable ground water models require both an accurate physical representation of the system and appropriate boundary conditions. While physical attributes are generally considered static, boundary conditions, such as ground water recharge rates, can be highly variable in both space and time. A practical methodology incorporating the hydrologic model HELP3 in conjunction with a geographic information system was developed to generate a physically based and highly detailed recharge boundary condition for ground water modeling. The approach uses daily precipitation and temperature records in addition to land use/land cover and soils data. The importance of the method in transient ground water modeling is demonstrated by applying it to a MODFLOW modeling study in New Jersey. In addition to improved model calibration, the results from the study clearly indicate the importance of using a physically based and highly detailed recharge boundary condition in ground water quality modeling, where the detailed knowledge of the evolution of the ground water flowpaths is imperative. The simulated water table is within 0.5 m of the observed values using the method, while the water levels can differ by as much as 2 m using uniform recharge conditions. The results also show that the combination of temperature and precipitation plays an important role in the amount and timing of recharge in cooler climates. A sensitivity analysis further reveals that increasing the leaf area index, the evaporative zone depth, or the curve number in the model will result in decreased recharge rates over time, with the curve number having the greatest impact.     

Effects of Linking a Soil‐Water‐Balance Model with a Groundwater‐Flow Model

A previously published regional groundwater‐flow model in north‐central Nebraska was sequentially linked with the recently developed soil‐water‐balance (SWB) model to analyze effects to groundwater‐flow model parameters and calibration results. The linked models provided a more detailed spatial and temporal distribution of simulated recharge based on hydrologic processes, improvement of simulated groundwater‐level changes and base flows at specific sites in agricultural areas, and a physically based assessment of the relative magnitude of recharge for grassland, nonirrigated cropland, and irrigated cropland areas. Root‐mean‐squared (RMS) differences between the simulated and estimated or measured target values for the previously published model and linked models were relatively similar and did not improve for all types of calibration targets. However, without any adjustment to the SWB‐generated recharge, the RMS difference between simulated and estimated base‐flow target values for the groundwater‐flow model was slightly smaller than for the previously published model, possibly indicating that the volume of recharge simulated by the SWB code was closer to actual hydrogeologic conditions than the previously published model provided. Groundwater‐level and base‐flow hydrographs showed that temporal patterns of simulated groundwater levels and base flows were more accurate for the linked models than for the previously published model at several sites, particularly in agricultural areas.     

Distinguishing sources of ground water recharge by using delta2H and delta18O

Stable isotope values of hydrogen and oxygen from precipitation and ground water samples were compared by using a volumetrically based mixing equation and stable isotope gradient to estimate the season and location of recharge in four basins. Stable isotopes were sampled at 11 precipitation sites of differing elevation during a 2-year period to quantify seasonal stable isotope contributions as a function of elevation. Supplemental stable isotope data collected by the International Atomic Energy Association during a 14-year period were used to reduce annual variability of the mean seasonal stable isotope data. The stable isotope elevation relationships and local precipitation elevation relationships were combined by using a digital elevation model to calculate the total volumetric contribution of water and stable isotope values as a function of elevation within the basins. The results of these precipitation calculations were compared to measured ground water stable isotope values at the major discharge points near the terminus of the basins. Volumetric precipitation contributions to recharge were adjusted to isolate contributing elevations. This procedure provides an improved representation of recharge contributions within the basins over conventional stable isotope methods. Stable isotope values from wells and springs at the terminus of each basin were used to infer the elevations of precipitation important for recharge of the regional ground water flow system. Ancillary climatic, geologic, and stable isotope values were used to further constrain the location where precipitation is entering the ground water flow system.     

Ground water budget analysis and cross-formational leakage in an arid basin

Ground water budget analysis in arid basins is substantially aided by integrated use of numerical models and environmental isotopes. Spatial variability of recharge, storage of water of both modern and pluvial age, and complex three-dimensional flow processes in these basins provide challenges to the development of a good conceptual model. Ground water age dating and mixing analysis with isotopic tracers complement standard hydrogeologic data that are collected and processed as an initial step in the development and calibration of a numerical model. Environmental isotopes can confirm or refute a priori assumptions of ground water flow, such as the general assumption that natural recharge occurs primarily along mountains and mountain fronts. Isotopes also serve as powerful tools during postaudits of numerical models. Ground water models provide a means of developing ground water budgets for entire model domains or for smaller regions within the model domain. These ground water budgets can be used to evaluate the impacts of pumping and estimate the magnitude of capture in the form of induced recharge from streams, as well as quantify storage changes within the system. The coupled analyses of ground water budget analysis and isotope sampling and analysis provide a means to confirm, refute, or modify conceptual models of ground water flow.     

Transient analysis of the source of water to wells: Cape Cod, Massachusetts

A transient flow modeling analysis for potential public-supply wells on western Cape Cod, Massachusetts, demonstrates the difference between transient and steady-state recharge areas can have important implications for wellhead protection. An example of a single pumping well illustrates that commonly used steady-state time-related capture areas do not represent the recharge area and travel times of water being pumped from the well until sufficient time has elapsed for steady-state flow conditions to be established. Until that time, transient recharge areas are needed to account for the portion of water discharging from the well that entered the aquifer before pumping started. An example of two pumping wells demonstrates the same area at the water table cannot supply water to more than one well under steady-state flow conditions. Transient recharge areas to multiple wells can overlap, however, until steady-state flow conditions are established. The same area can, therefore, be a source of water to more than one well during early pumping times, and the water pumped from a given well may derive from source areas, including contaminated areas, that do not lie within the well's steady-state recharge area.     

Recharge process and aquifer models of a small watershed

Abstract

Kanchanapally watershed covering an area of about 11 km² in Nalgonda district, Andhra Pradesh, India is located in granitic terrain. Groundwater recharge has been estimated from a water balance model using hydrometeorological data from 1978–1994. The monthly recharge estimates obtained from the water balance model formed input for the groundwater flow model during transient model testing. The groundwater flow model has been prepared to simulate steady state groundwater conditions of 1977 using the nested squares finite difference method. The transient groundwater flow model has been tested during 1977–1994 using the estimated recharge values. The present study helped verify the usefulness of monthly recharge estimates for accounting dynamic variations in recharge as reflected in water level fluctuations in hydrographs.     

Development and application of an analytical model of stream/aquifer interaction

A mathematical model to simulate stream/aquifer interactions in an unconfined aquifer subjected to time varying river stage was developed from the linearized Boussinesq equation using the principle of superposition and the concept of semigroups. The mathematical model requires an estimate of three parameters to simulate ground-water elevations; transmissivity, specific yield, and recharge. The solution has physical significance and includes terms for the steady-state water level, the steady-state water level as influenced by a change in river stage, a transient redistribution of water levels in the aquifer from the previous day, and a transient change in water level caused by a change in river stage. The mathematical model was tested using observed water table elevations at three locations across a 2-km-wide alluvial valley aquifer. The average absolute deviation between observed and simulated daily water levels was 0.09 m. The difference in river stage over the test year was 4.9 m.     

Ground water flow analysis of a mid-Atlantic outer coastal plain watershed, Virginia, U.S.A

Models for ground water flow (MODFLOW) and particle tracking (MODPATH) were used to determine ground water flow patterns, principal ground water discharge and recharge zones, and estimates of ground water travel times in an unconfined ground water system of an outer coastal plain watershed on the Delmarva Peninsula, Virginia. By coupling recharge and discharge zones within the watershed, flowpath analysis can provide a method to locate and implement specific management strategies within a watershed to reduce ground water nitrogen loading to surface water. A monitoring well network was installed in Eyreville Creek watershed, a first-order creek, to determine hydraulic conductivities and spatial and temporal variations in hydraulic heads for use in model calibration. Ground water flow patterns indicated the convergence of flow along the four surface water features of the watershed; primary discharge areas were in the nontidal portions of the watershed. Ground water recharge zones corresponded to the surface water features with minimal development of a regional ground water system. Predicted ground water velocities varied between < 0.01 to 0.24 m/day, with elevated values associated with discharge areas and areas of convergence along surface water features. Some ground water residence times exceeded 100 years, although average residence times ranged between 16 and 21 years; approximately 95% of the ground water resource would reflect land use activities within the last 50 years.     

Important observations and parameters for a salt water intrusion model

Sensitivity analysis with a density-dependent ground water flow simulator can provide insight and understanding of salt water intrusion calibration problems far beyond what is possible through intuitive analysis alone. Five simple experimental simulations presented here demonstrate this point. Results show that dispersivity is a very important parameter for reproducing a steady-state distribution of hydraulic head, salinity, and flow in the transition zone between fresh water and salt water in a coastal aquifer system. When estimating dispersivity, the following conclusions can be drawn about the data types and locations considered. (1) The "toe" of the transition zone is the most effective location for hydraulic head and salinity observations. (2) Areas near the coastline where submarine ground water discharge occurs are the most effective locations for flow observations. (3) Salinity observations are more effective than hydraulic head observations. (4) The importance of flow observations aligned perpendicular to the shoreline varies dramatically depending on distance seaward from the shoreline. Extreme parameter correlation can prohibit unique estimation of permeability parameters such as hydraulic conductivity and flow parameters such as recharge in a density-dependent ground water flow model when using hydraulic head and salinity observations. Adding flow observations perpendicular to the shoreline in areas where ground water is exchanged with the ocean body can reduce the correlation, potentially resulting in unique estimates of these parameter values. Results are expected to be directly applicable to many complex situations, and have implications for model development whether or not formal optimization methods are used in model calibration.     

Transient Recharge Estimability Through Field‐Scale Groundwater Model Calibration

The estimation of recharge through groundwater model calibration is hampered by the nonuniqueness of recharge and aquifer parameter values. It has been shown recently that the estimability of spatially distributed recharge through calibration of steady‐state models for practical situations (i.e., real‐world, field‐scale aquifer settings) is limited by the need for excessive amounts of hydraulic‐parameter and groundwater‐level data. However, the extent to which temporal recharge variability can be informed through transient model calibration, which involves larger water‐level datasets, but requires the additional consideration of storage parameters, is presently unknown for practical situations. In this study, time‐varying recharge estimates, inferred through calibration of a field‐scale highly parameterized groundwater model, are systematically investigated subject to changes in (1) the degree to which hydraulic parameters including hydraulic conductivity (K) and specific yield (S_y) are constrained, (2) the number of water‐level calibration targets, and (3) the temporal resolution (up to monthly time steps) at which recharge is estimated. The analysis involves the use of a synthetic reality (a reference model) based on a groundwater model of Uley South Basin, South Australia. Identifiability statistics are used to evaluate the ability of recharge and hydraulic parameters to be estimated uniquely. Results show that reasonable estimates of monthly recharge (<30% recharge root‐mean‐squared error) require a considerable amount of transient water‐level data, and that the spatial distribution of K is known. Joint estimation of recharge, S_y and K, however, precludes reasonable inference of recharge and hydraulic parameter values. We conclude that the estimation of temporal recharge variability through calibration may be impractical for real‐world settings.     

Delineating ground water recharge from leaking irrigation canals using water chemistry and isotopes

Across the Great Plains irrigation canals are used to transport water to cropland. Many of these canals are unlined, and leakage from them has been the focus of an ongoing legal, economic, and philosophical debate as to whether this lost water should be considered waste or be viewed as a beneficial and reasonable use since it contributes to regional ground water recharge. While historically there has been much speculation about the impact of canal leakage on local ground water, actual data are scarce. This study was launched to investigate the impact of leakage from the Interstate Canal, in the western panhandle of Nebraska, on the hydrology and water quality of the local aquifer using water chemistry and environmental isotopes. Numerous monitoring wells were installed in and around a small wetland area adjacent to the canal, and ground water levels were monitored from June 1992 until January 1995. Using the water level data, the seepage loss from the canal was estimated. In addition, the canal, the monitoring wells, and several nearby stock and irrigation wells were sampled for inorganic and environmental isotope analysis to assess water quality changes, and to determine the extent of recharge resulting from canal leakage. The results of water level monitoring within study wells indicates a rise in local ground water levels occurs seasonally as a result of leakage during periods when the canal is filled. This rise redirects local ground water flow and provides water to nearby wetland ecosystems during the summer months. Chemical and isotopic results were used to delineate canal, surface, and ground water and indicate that leaking canal water recharges both the surface alluvial aquifer and upper portions of the underlying Brule Aquifer. The results of this study indicate that lining the Interstate Canal could lower ground water levels adjacent to the canal, and could adversely impact the local aquifer.     

Importance of unsaturated zone flow for simulating recharge in a humid climate

Transient recharge to the water table is often not well understood or quantified. Two approaches for simulating transient recharge in a ground water flow model were investigated using the Trout Lake watershed in north-central Wisconsin: (1) a traditional approach of adding recharge directly to the water table and (2) routing the same volume of water through an unsaturated zone column to the water table. Areas with thin (less than 1 m) unsaturated zones showed little difference in timing of recharge between the two approaches; when water was routed through the unsaturated zone, however, less recharge was delivered to the water table and more discharge occurred to the surface because recharge direction and magnitude changed when the water table rose to the land surface. Areas with a thick (15 to 26 m) unsaturated zone were characterized by multimonth lags between infiltration and recharge, and, in some cases, wetting fronts from precipitation events during the fall overtook and mixed with infiltration from the previous spring snowmelt. Thus, in thicker unsaturated zones, the volume of water infiltrated was properly simulated using the traditional approach, but the timing was different from simulations that included unsaturated zone flow. Routing of rejected recharge and ground water discharge at land surface to surface water features also provided a better simulation of the observed flow regime in a stream at the basin outlet. These results demonstrate that consideration of flow through the unsaturated zone may be important when simulating transient ground water flow in humid climates with shallow water tables.     

Modeling a Ground Water Circulation Well Alternative

Ground water circulation wells (GCWs) provide an appealing alternative to typical pump-and-treat ground water remediation systems because of the inherent resource-conservative nature of the GCW systems. GCW performance prediction is challenging because the consideration of extraction and recharge in a single well is unusual for most practitioners, the technology is relatively new, and a meaningful body of literature has not been published. A three-part evaluation process using state-of-the-practice numerical ground water flow and mass transport models was developed for application during GCW pilot studies at the former Nebraska Ordnance Plant site. A small-scale ground water flow model was developed during the pilot study planning process to predict the system performance and to locate performance-measuring monitoring wells. Key predictions included the capture zone predicted to develop upgradient of the GCW, the downgradierit GCW recharge zone, and the circulation zone centered on the GCW. The flow model was subsequently verified using ground water elevation data and contaminant concentration data collected during pilot study operation. Aquifer parameters were reestimated as a result of the verification process. Those parameter values were used as input to a larger scale model, which was used to develop a remedial alternative consisting of multiple GCW systems.     

Influence of calibration methodology on ground water flow predictions

We constructed a numerical model of transient ground water flow and solute transport for a portion of the Biscayne Aquifer in Florida, and calibrated the model with three different combinations of data from a 193-day period: head (h) data alone, data on h and ground water discharge to a canal (q), and data on h, q, and ground water chloride concentration (C). We used each of the three calibrated models to predict h and q during a 182-day test period separate from the calibration period. All three calibrated models predicted h equally well during the test period (r = 0.95, where r = 1 indicates perfect agreement between measured and simulated values), though the model calibrated on h alone had significantly different parameter values than the other two models. Predictions of q during the test period depended on calibration methodology; models calibrated with multiple targets simulated q more accurately than the model calibrated on h alone (r = 0.79 compared to r = 0.49). Based on the results of these simulations, we conclude: (1) Post-calibration prediction is important in assessing the value of different data types in automated calibration; (2) inverse-solution uniqueness is not a requirement for accurate h predictions; (3) relatively simple models can predict with reasonable accuracy transient ground water flow in a complex aquifer, and parameters governing this prediction can be estimated by nonlinear regression methods that incorporate both h and q data; (4) addition of C data to the calibration did not improve model predictive capacity because the information in the C data was similar to that in the q data, from the perspective of model calibration (the subsurface chemical signal in question was controlled mainly by seepage of high-chloride canal water into the low-chloride ground water system).     

Estimating recharge rates with analytic element models and parameter estimation

Quantifying the spatial and temporal distribution of recharge is usually a prerequisite for effective ground water flow modeling. In this study, an analytic element (AE) code (GFLOW) was used with a nonlinear parameter estimation code (UCODE) to quantify the spatial and temporal distribution of recharge using measured base flows as calibration targets. The ease and flexibility of AE model construction and evaluation make this approach well suited for recharge estimation. An AE flow model of an undeveloped watershed in northern Wisconsin was optimized to match median annual base flows at four stream gages for 1996 to 2000 to demonstrate the approach. Initial optimizations that assumed a constant distributed recharge rate provided good matches (within 5%) to most of the annual base flow estimates, but discrepancies of >12% at certain gages suggested that a single value of recharge for the entire watershed is inappropriate. Subsequent optimizations that allowed for spatially distributed recharge zones based on the distribution of vegetation types improved the fit and confirmed that vegetation can influence spatial recharge variability in this watershed. Temporally, the annual recharge values varied >2.5-fold between 1996 and 2000 during which there was an observed 1.7-fold difference in annual precipitation, underscoring the influence of nonclimatic factors on interannual recharge variability for regional flow modeling. The final recharge values compared favorably with more labor-intensive field measurements of recharge and results from studies, supporting the utility of using linked AE-parameter estimation codes for recharge estimation.     

Analysis of recharge conceptualization in inverse groundwater modelling

ABSTRACT

This study presents a systematic illustration quantifying how misleading the calibration results of a groundwater simulation model can be when recharge rates are considered as the model parameters to be estimated by inverse modelling. Three approaches to recharge estimation are compared: autocalibration (Model 1), the empirical return coefficient method (Model 2), and distributed hydrological modelling using the Soil and Water Assessment Tool, SWAT (Model 3). The methodology was applied in the Dehloran Plain, western Iran, using the MODFLOW modular flow simulator and the PEST method for autocalibration. The results indicate that, although Model 1 performed the best in simulating water levels at observation wells in the calibration stage, it did not perform satisfactorily in real future scenarios. Model 3, with SWAT-based recharge rates, performed better than the other models in the validation stage. By not evaluating the model performance solely on calibration results, we demonstrate the relative significance of using more accurate recharge estimates when calibrating groundwater simulation models.

EDITOR D. Koutsoyiannis; ASSOCIATE EDITOR M. Besbes     

Estimation of recharge from floods in disconnected stream-aquifer systems

Stream-aquifer interaction has been the subject of much research for cases of good hydraulic connection (continuous saturated zone) between a river and an aquifer. Under these conditions, floods do not represent a very large net input to the aquifer because most of the water that enters the aquifer during the flood returns to the river when its stage recedes. The situation is different in disconnected stream-aquifer systems, where the streambed lies above the water level in the aquifer, thus preventing return flow from the aquifer. Under these conditions, floods may represent large, but hard to quantify, water inputs. Here, we present a methodology to estimate recharge from floods for disconnected stream-aquifer systems. Recharge is estimated as the product of a flood time function (dependent on the streamflow) and an unknown factor, which is obtained from calibrating a ground water flow model to aquifer heads. The approach can also benefit from concentration data, which can be very informative when river water concentrations vary over time. This methodology is applied to a field situation where recharge from river flooding is found to amount to nearly 15 million m(3)/year on the average, which represents 40% of the total aquifer inputs. Recharge from flooding helps explain major head recoveries, suggesting that basin water management programs should allow some floods to occur.     

The role of hand calculations in ground water flow modeling

Most ground water modeling courses focus on the use of computer models and pay little or no attention to traditional analytic solutions to ground water flow problems. This shift in education seems logical. Why waste time to learn about the method of images, or why study analytic solutions to one-dimensional or radial flow problems? Computer models solve much more realistic problems and offer sophisticated graphical output, such as contour plots of potentiometric levels and ground water path lines. However, analytic solutions to elementary ground water flow problems do have something to offer over computer models: insight. For instance, an analytic one-dimensional or radial flow solution, in terms of a mathematical expression, may reveal which parameters affect the success of calibrating a computer model and what to expect when changing parameter values. Similarly, solutions for periodic forcing of one-dimensional or radial flow systems have resulted in a simple decision criterion to assess whether or not transient flow modeling is needed. Basic water balance calculations may offer a useful check on computer-generated capture zones for wellhead protection or aquifer remediation. An easily calculated "characteristic leakage length" provides critical insight into surface water and ground water interactions and flow in multi-aquifer systems. The list goes on. Familiarity with elementary analytic solutions and the capability of performing some simple hand calculations can promote appropriate (computer) modeling techniques, avoids unnecessary complexity, improves reliability, and is likely to save time and money. Training in basic hand calculations should be an important part of the curriculum of ground water modeling courses.     

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.