Monitoring of drawdown pattern during pumping in an unconfined physical aquifer model

In this study, we attempted to analyse a drawdown pattern around a pumping well in an unconfined sandy gravelly aquifer constructed in a laboratory tank by means of both experimental and numerical modelling of groundwater flow. The physical model consisted of recharge, aquifer and discharge zones. Permeability and specific yield of the aquifer material were determined by Dupuit approximation under steady‐state flow and stepwise gravitational drainage of groundwater, respectively. The drawdown of water table in pumping and neighbouring observation wells was monitored to investigate the effect of no‐flow boundary on the drawdown pattern during pumping for three different boundary conditions: (i) no recharge and no discharge with four no‐flow boundaries (Case 1); (ii) no recharge and reservoir with three no‐flow boundaries (Case 2); (iii) recharge and discharge with two no‐flow boundaries (Case 3). Based on the aquifer parameters, numerical modelling was also performed to compare the simulated drawdown with that observed. Results showed that a large difference existed between the simulated drawdown and that observed in wells for all cases. The reason for the difference could be explained by the formation of a curvilinear type water table between wells rather than a linear one due to a delayed response of water table in the capillary fringe. This phenomenon was also investigated from a mass balance study on the pumping volume. The curvilinear type of water table was further evidenced by measurement of water contents at several positions in the aquifer between wells using time domain reflectometry (TDR). This indicates that the existing groundwater flow model applicable to an unconfined aquifer lacks the capacity to describe a slow response of water table in the aquifer and care should be taken in the interpretation of water table formation in the aquifer during pumping. Copyright © 2001 John Wiley & Sons, Ltd.     

USE OF TRITIUM IN HYDROLOGY

SUMMARY

The Coastal Plain aquifer of Israel, of Plio-Pleistocene Age, stretches from Binyamina in the North to the Gaza Strip in the South-a distance of about 112 km and has an average width of about 15 Km. The allowed withdrawal is estimated at about 200 MCM/year.

As a result of an average yearly withdrawal of 426 MCM/year during the last 10 years the water levels dropped to a dangerously low position (-2)-(-4) m below sea level at distances of 3–5 Km from the coast, causing sea water intrusion which, in Tel Aviv and Emek Hefer, endangered water supply wells.

As a counter-measure, artificial groundwater recharge through wells was practiced in Emek Hefer since 1959. Recharge was practiced in 7 wells at a rate of 6 MCM/year, the water coming from adjacent Cretaceous limestone aquifers.

In Tel Aviv a fresh water barrier was established in 1964 by injecting Lake Kinereth water into 17 wells during winter at a rate of 6 MCM/winter. In the rest of the Coastal Plain water was injected to the aquifer through about 40–45 wells at a total yearly rate of about 10–12 MCM.

Recharge by spreading is practiced in Yavneh at a rate of about 10–13 MCM per winter, also recharge by spreading is practiced with flood water of Nahal Shikma at a rate of up to 8 MCM/winter.     

Atrazine in a Stream-Aquifer System: Estimation of Aquifer Properties from Atrazine Concentration Profiles

Lincolns municipal wellfield consists of 44 wells developed in an alluvial aquifer adjacent to the Platte River near Ashland, Nebraska Induced recharge from the river is the primary source of water for the wellfield. Wafer samples were collected on a periodic basis from the Platte River arid two transects of monitoring wells. These samples were analyzed for the herbicide atrazine, which was used as a tracer of induced recharge in this stream-aquifer system. Atrazine concentrations in the river and aquifer were much less than 1.0 ppb during late fall and winter, but increased to as high as 18.9 ppb during spring and summer, associated with runoff from upgradient agricultural lands. There was approximately a 21-day lag time from the first detection of increasing atrazine concentration in the river to the first detection in monitoring wells immediately adjacent to the river. This lag time was relatively constant throughout the year and from one year to the next, even with major fluctuations of river stage and wellfield production. This consistency of lag time indicated that the travel times from the river to the first set of monitoring wells immediately adjacent to the river were fairly constant.
Paths of preferential flow were identified in the aquifer at a depth of 25 to 35 feet below land surface. This aquifer zone appeared to play a significant role in movement of water from beneath the river into the wellfield.
Aquifer dispersivity was calculated using a method described by Hoehn and Santschi (1987). Macrodispersivity (A_L) was shown to increase linearly over the scale of the wellfield. Calculated values of A_L were within limits of other reported values for this type of aquifer material and agreed well with values reported by Hoehn and Santschi (1987); These findings will be extremely beneficial for planning and management of the municipal wellfield.     

Spatio-temporal analysis of potential aquifer recharge: Application to the Basin of Mexico

Regional estimates of aquifer recharge are needed in data-scarce regions such as the Basin of Mexico, where nearly 20 million people are located and where the Basin’s aquifer system represents the main water source. In order to develop the spatio-temporal estimates of aquifer recharge and to analyze to what extent urban growth has affected aquifer recharge, this work presents a daily soil water balance which uses different vegetation and soil types as well as the effect of topography on climatological variables and evapotranspiration. The soil water balance was applied on a daily time step in the Basin of Mexico for the period 1975–1986, obtaining an annually-lumped potential recharge flow of 10.9–23.8 m³/s (35.9–78.1 mm) in the entire Basin, while the monthly values for the year with the largest lumped recharge value (1981 = 78.1 mm) range from 1 m³/s (0.3 mm) in December to 87.9 m³/s (23.7 mm) in June. As aquifer recharge in the Basin mainly occurs by subsurface flow from its enclosing mountains as Mountain Block Recharge, urban growth has had a minimal impact on aquifer recharge, although it has diminished recharge in the alluvial plain.     

Influence of aquifer properties and the spatial and temporal distribution of recharge and abstraction on sustainable yields in semi‐arid regions

This paper aims to contribute to understanding the importance of four factors on the determination of sustainable yields: (i) aquifer properties; (ii) temporal distribution of recharge; (iii) temporal distribution of groundwater pumping; and (iv) spatial distribution of pumping wells. It is important to comprehend how the present‐day and future vulnerability of groundwater systems to pumping activities depend on these critical factors and what the risks are of considering sustainable yield as a fixed percentage of mean annual recharge (MAR). A numerical model of the Querença–Silves aquifer in Portugal is used to develop hypothetical scenarios with which these factors are studied. Results demonstrate the aquifer properties, particularly the storage coefficient, have an important role in determining the resilience of an aquifer and therefore to which degree it is dependent on the spatial and temporal distribution of abstraction and recharge, as well as the occurrence of extreme events. Sustainable yields are determined for the developed scenarios based on specific criteria rather than a fraction of MAR. Under simplified current recharge and abstraction conditions, the sustainable yield was determined at approximately 73% of MAR or 76 million m³. When considering a concentration of rainfall in time, as predicted by climate scenarios for the region, sustainable yield could drop to ca 70% of MAR. However, a more even distribution of pumping volumes throughout the year could increase this value. The location of the pumping wells is seen to affect the distribution of hydraulic heads in the aquifer, albeit without significant changes in sustainable yield. Copyright © 2011 John Wiley & Sons, Ltd.     

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.     

Strategies for offsetting seasonal impacts of pumping on a nearby stream

Ground water pumping from aquifer systems that are hydraulically connected to streams depletes streamflow. The amplitude and timing of stream depletion depend on the stream depletion factor (SDF(i)) of the pumping wells, which is a function of aquifer hydraulic characteristics and the distance from the wells to the stream. Wells located at different locations, but having the same SDF and the same rate and schedule of pumping, will deplete streamflow equally. Wells with small SDF(i) deplete streamflow approximately synchronously with pumping. Wells with large SDF(i) deplete streamflow at approximately a constant rate throughout the year, regardless of the pumping schedule. For large values of SDF(i), artificial recharge that occurs on a different schedule from pumping can offset streamflow depletion effectively. The requirements are (1) that the pumping and recharge wells both have the same SDF(i) and (2) that the annual total quantities of recharge and pumping be equal. At larger SDF(i) values, it takes longer for pumping to impact streamflow in a wide aquifer than it does in a narrow aquifer. In basins that are closed to further withdrawals because streamflow is fully allocated, water-use changes replace new allocations as the source of water for new developments. Ground water recharge can be managed to offset the impacts of new ground water developments, allowing for changes in the timing and source of withdrawals from a basin without injuring existing users or instream flows.     

Episodic and Continuous Recharge Estimation from High-Resolution Well Records

Water table fluctuation (WTF) methods are a primary and well-established way to determine groundwater recharge based on the direct response of the water table to precipitation input. An emerging complexity of recharge is whether it occurs as an episodic and transient process, or a continuous steady-state process, however, most studies have not focused on these short-term vs. long-term timescales, in part because of a lack of data resolution. Here, high-resolution (subhourly) precipitation and water-level data are analyzed for wells in the suburbs of New York City using two contrasting WTF approaches, with a common mathematical basis, that are suited to episodic and continuous processes. The resulting hourly recharge results, like the individual water-level records from comparable wells, are sensitive indicators of subtle differences in aquifer conditions such as thickness of the unsaturated zone, position in the flow system and localized preferential flow. While the episodic, transient approach excludes diffuse recharge by design, the continuous, steady-state approach is influenced by short-term precipitation events, and therefore integrates transient processes to some extent. However, the continuous, steady-state approach is subject to its own limitations relating to position in the aquifer system, and may overestimate recharge if aquifer conditions are not well understood. More widespread use of higher resolution data as well as understanding aquifer conditions and refining aquifer parameters would improve WTF recharge estimation.     

Atrazine in a Stream-Aquifer System: Transport of Atrazine and Its Environmental Impact Near Ashland, Nebraska

Lincoln's municipal wellfield consists of 44 wells located adjacent to the Platte River near Ashland, Nebraska. The herbicide atrazine was monitored in the river and two transects of monitoring, wells. The amount of atrazine transported down the Platte River in 1989, 1990, and 1991 was shown to increase each year. Induced recharge from the Platte River results in movement of atrazine from the river into the aquifer. A 21-day lag time was determined for the movement of atrazine from the river to a transect of monitoring wells 10 feet West of the bank. The role that colloids play on the transport of atrazine was determined to be insignificant. A small percentage of atrazine found in the river was determined to come from rain water. The infiltration of agrichemical-contaminated river water was shown to significantly reduce the quality of raw water and finished water being produced by the adjacent aquifer.     

Numerical modelling of stream–aquifer interaction: Quantifying the impact of transient streambed permeability and aquifer heterogeneity

Stream–aquifer interaction plays a vital role in the water cycle, and a proper study of this interaction is needed for understanding groundwater recharge, contaminants migration, and for managing surface water and groundwater resources. A model‐based investigation of a field experiment in a riparian zone of the Schwarzbach river, a tributary of the Rhine River in Germany, was conducted to understand stream–aquifer interaction under alternative gaining and losing streamflow conditions. An equivalent streambed permeability, estimated by inverting aquifer responses to flood waves, shows that streambed permeability increased during infiltration of stream water to aquifer and decreased during exfiltration. Aquifer permeability realizations generated by multiple‐point geostatistics exhibit a high degree of heterogeneity and anisotropy. A coupled surface water groundwater flow model was developed incorporating the time‐varying streambed permeability and heterogeneous aquifer permeability realizations. The model was able to reproduce varying pressure heads at two observation wells near the stream over a period of 55 days. A Monte Carlo analysis was also carried out to simulate groundwater flow, its age distribution, and the release of a hypothetical wastewater plume into the aquifer from the stream. Results of this uncertainty analysis suggest (a) stream–aquifer exchange flux during the infiltration periods was constrained by aquifer permeability; (b) during exfiltration, this flux was constrained by the reduced streambed permeability; (c) the effect of temporally variable streambed permeability and aquifer heterogeneity were found important to improve the accurate capture of the uncertainty; and (d) probabilistic infiltration paths in the aquifer reveal that such pathways and the associated prediction of the extent of the contaminant plume are highly dependent on aquifer heterogeneity.     

A first estimate of ground water ages for the deep aquifer of the Kathmandu Basin, Nepal, using the radioisotope chlorine-36

The Kathmandu Basin in Nepal contains up to 550 m of Pliocene-Quaternary fluvio-lacustrine sediments which have formed a dual aquifer system. The unconfined sand and gravel aquifer is separated by a clay aquitard, up to 200 m thick, from the deeper, confined aquifer, comprised of Pliocene sand and gravel beds, intercalated with clay, peat, and lignite. The confined aquifer currently provides an important water supply to the central urban area but there are increasing concerns about its sustainability due to overexploitation. A limited number of determinations of the radioisotope 36Cl have been made on bore waters in the basin, allowing us to postulate on the age of ground water in the deeper, confined aquifer. Ground water evolution scenarios based on radioisotope decay, gradual dissolution of formational salts as the ground waters move downgradient, and flow velocity estimations produce comparable ground water ages for the deep waters, ranging from 200,000 to 400,000 years. From these ages, we deduce a mean ground water flow velocity of only 45 mm/year from recharge in the northeast to the main extraction region 15 km to the southwest. We thus estimate current recharge at about 5 to 15 mm/year, contributing 40,000 to 1.2 million m3/year to the ground water system. Current ground water extraction is estimated to be 20 times this amount. The low specific discharge confirms that the resource is being mined, and, based on current projections, reserves will be used up within 100 years.     

Estimate of recharge of a rising water table in semiarid Niger from 3H and 14C modeling

A hydrodynamic survey carried out in semiarid southwest Niger revealed an increase in the unconfined ground water reserves of approximately 10% over the last 50 years due to the clearing of native vegetation. Isotopic samplings (3H, 18O, 2H for water and 14C, 13C for the dissolved inorganic carbon) were performed on about 3500 km2 of this silty aquifer to characterize recharge. Stable isotope analyses confirmed the indirect recharge process that had already been shown by hydrodynamic surveys and suggested the tracers are exclusively of atmospheric origin. An analytical model that takes into account the long-term rise in the water table was used to interpret 3H and 14C contents in ground water. The natural, preclearing median annual renewal rate (i.e., recharge as a fraction of the saturated aquifer volume) lies between 0.04% and 0.06%. For representative characteristics of the aquifer (30 m of saturated thickness, porosity between 10% and 25%), this implies a recharge of between 1 and 5 mm/year, which is much lower than the estimates of 20 to 50 mm/year for recent years, obtained using hydrological and hydrodynamic methods and the same aquifer parameters. Our study, therefore, reveals that land clearing in semiarid Niger increased ground water recharge by about one order of magnitude.     

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

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

Interaction of Aquifer and River‐Canal Network near Well Field

The article presents semi‐analytical mathematical models to asses (1) enhancements of seepage from a canal and (2) induced flow from a partially penetrating river in an unconfined aquifer consequent to groundwater withdrawal in a well field in the vicinity of the river and canal. The nonlinear exponential relation between seepage from a canal reach and hydraulic head in the aquifer beneath the canal reach is used for quantifying seepage from the canal reach. Hantush's (1967) basic solution for water table rise due to recharge from a rectangular spreading basin in absence of pumping well is used for generating unit pulse response function coefficients for water table rise in the aquifer. Duhamel's convolution theory and method of superposition are applied to obtain water table position due to pumping and recharge from different canal reaches. Hunt's (1999) basic solution for river depletion due to constant pumping from a well in the vicinity of a partially penetrating river is used to generate unit pulse response function coefficients. Applying convolution technique and superposition, treating the recharge from canal reaches as recharge through conceptual injection wells, river depletion consequent to variable pumping and recharge is quantified. The integrated model is applied to a case study in Haridwar (India). The well field consists of 22 pumping wells located in the vicinity of a perennial river and a canal network. The river bank filtrate portion consequent to pumping is quantified.     

Applicability of spectral analysis to determine hydraulic diffusivity

This study is to evaluate the applicability of estimating the one-dimensional horizontal hydraulic diffusivity of an unconfined aquifer with time-dependent fluctuation of lateral head and vertical recharge boundaries using observed water level spectra. Different models of boundary condition are imposed to evaluate the statistical significance between the calculated hydraulic diffusivity (ξ^′) with the given hydraulic diffusivity (ξ). The auto-spectra of the water level in observation wells tapping the same aquifer are closely related to those at the disturbed boundaries. For an aquifer with a constant hydraulic diffusivity, the water level fluctuation in the monitoring wells is linearly related to the water level spectra observed at the boundaries. The spectral density function of aquifer hydraulic head varies inversely with specific yield (S _y) and directly with recharge. Given small variation in water level spectra at the disturbed boundaries, the water level fluctuation in the aquifer is affected by the recharge condition and the aquifer spectral density function is sensitive to S _y. Using an iterative technique to estimate ξ from 1400 sets of given parameters, 99% of the ξ^′/ξ values deviated within only one order of magnitude with the model length (L) being equal to 1 km and 10 km. For L equal to 100 m, approximately 82% of the ξ^′/ξ population falls within two orders of magnitude. Therefore, spectral analysis of aquifer hydraulic head response can be used to estimate the hydraulic diffusivity of an unconfined aquifer which is affected by periodic variations in recharge and head at boundaries.     

Modeling decadal timescale interactions between surface water and ground water in the central Everglades, Florida, USA

Surface-water and ground-water flow are coupled in the central Everglades, although the remoteness of this system has hindered many previous attempts to quantify interactions between surface water and ground water. We modeled flow through a 43,000 ha basin in the central Everglades called Water Conservation Area 2A. The purpose of the model was to quantify recharge and discharge in the basin's vast interior areas. The presence and distribution of tritium in ground water was the principal constraint on the modeling, based on measurements in 25 research wells ranging in depth from 2 to 37 m. In addition to average characteristics of surface-water flow, the model parameters included depth of the layer of ‘interactive’ ground water that is actively exchanged with surface water, average residence time of interactive ground water, and the associated recharge and discharge fluxes across the wetland ground surface. Results indicated that only a relatively thin (8 m) layer of the 60 m deep surfical aquifer actively exchanges surface water and ground water on a decadal timescale. The calculated storage depth of interactive ground water was 3.1 m after adjustment for the porosity of peat and sandy limestone. Modeling of the tritium data yielded an average residence time of 90 years in interactive ground water, with associated recharge and discharge fluxes equal to 0.01 cm d⁻¹. ³H/³He isotopic ratio measurements (which correct for effects of vertical mixing in the aquifer with deeper, tritium-dead water) were available from several wells, and these indicated an average residence time of 25 years, suggesting that residence time was overestimated using tritium measurements alone. Indeed, both residence time and storage depth would be expected to be overestimated due to vertical mixing. The estimate of recharge and discharge (0.01 cm d⁻¹) that resulted from tritium modeling therefore is still considered reliable, because the ratio of residence time and storage depth (used to calculated recharge and discharge) is much less sensitive to vertical mixing compared with residence time alone. We conclude that a small but potentially significant component of flow through the Everglades is recharged to the aquifer and stored there for years to decades before discharged back to surface water. Long-term storage of water and solutes in the ground-water system beneath the wetlands has implications for restoration of Everglades water quality.     

Steady, annual, and monthly recharge implied by deep unconfined aquifer flow

We consider the response of a deep unconfined horizontal aquifer to steady, annual, and monthly recharge. A groundwater divide and a zero head reservoir constrain the aquifer, so that sinusoidal monthly and aperiodic annual recharge fluctuations create transient specific discharge near the reservoir and an unsteady water table elevation inland. One existing and two new long-term data sets from the Plymouth-Carver Aquifer in southeastern Massachusetts calibrate and confirm hydraulic properties in a set of analytical models. [Geohydrology and simulated groundwater flow, 1992] data and a new power law for tritiugenic helium to tritium ratios calibrate the steady recharge that drives the classical parabolic model of steady hydraulics [Applied Hydrogeology, 2001]. Observed water table and gradient fluctuations calibrate the transient recharge models. In the latter regard, monitoring wells within 1 km of Buttermilk Bay exhibit appreciable specific discharge and reduced water table fluctuations. We apply [Trans Am Geophys Union 32(1951)238] periodic model to the monthly hydraulics and a recharge convolution integral [J Hydrol 126(1991)315] to annual flow. An infiltration fraction of 0.79 and a consumptive use coefficient of 1.08×10⁻⁸ m/s °C relate recharge to precipitation and daylight weighted temperature across all three time scales. Errors associated with this recharge relation decrease with increasing time scale.     

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.     

Electromagnetic‐Induction Logging to Monitor Changing Chloride Concentrations

Water from the San Joaquin Delta, having chloride concentrations up to 3590 mg/L, has intruded fresh water aquifers underlying Stockton, California. Changes in chloride concentrations at depth within these aquifers were evaluated using sequential electromagnetic (EM) induction logs collected during 2004 through 2007 at seven multiple‐well sites as deep as 268 m. Sequential EM logging is useful for identifying changes in groundwater quality through polyvinyl chloride‐cased wells in intervals not screened by wells. These unscreened intervals represent more than 90% of the aquifer at the sites studied. Sequential EM logging suggested degrading groundwater quality in numerous thin intervals, typically between 1 and 7 m in thickness, especially in the northern part of the study area. Some of these intervals were unscreened by wells, and would not have been identified by traditional groundwater sample collection. Sequential logging also identified intervals with improving water quality—possibly due to groundwater management practices that have limited pumping and promoted artificial recharge. EM resistivity was correlated with chloride concentrations in sampled wells and in water from core material. Natural gamma log data were used to account for the effect of aquifer lithology on EM resistivity. Results of this study show that a sequential EM logging is useful for identifying and monitoring the movement of high‐chloride water, having lower salinities and chloride concentrations than sea water, in aquifer intervals not screened by wells, and that increases in chloride in water from wells in the area are consistent with high‐chloride water originating from the San Joaquin Delta rather than from the underlying saline aquifer.     

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.