Estimating ground water recharge from topography, hydrogeology, and land cover

Proper management of ground water resources requires knowledge of the rates and spatial distribution of recharge to aquifers. This information is needed at scales ranging from that of individual communities to regional. This paper presents a methodology to calculate recharge from readily available ground surface information without long-term monitoring. The method is viewed as providing a reasonable, but conservative, first approximation of recharge, which can then be fine-tuned with other methods as time permits. Stream baseflow was measured as a surrogate for recharge in small watersheds in southeastern Wisconsin. It is equated to recharge (R) and then normalized to observed annual precipitation (P). Regression analysis was constrained by requiring that the independent and dependent variables be dimensionally consistent. It shows that R/P is controlled by three dimensionless ratios: (1) infiltrating to overland water flux, (2) vertical to lateral distance water must travel, and (3) percentage of land cover in the natural state. The individual watershed properties that comprise these ratios are now commonly available in GIS data bases. The empirical relationship for predicting R/P developed for the study watersheds is shown to be statistically viable and is then tested outside the study area and against other methods of calculating recharge. The method produces values that agree with baseflow separation from streamflow hydrographs (to within 15% to 20%), ground water budget analysis (4%), well hydrograph analysis (12%), and a distributed-parameter watershed model calibrated to total streamflow (18%). It has also reproduced the temporal variation over 5 yr observed at a well site with an average error < 12%.     

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.     

Groundwater recharge and palaeoclimate in the Sirte and Kufra basins, Libya

Stable and radio-isotope results (C, H, O) for groundwaters from the Sirte and northern Kufra basins are used to determine the recharge history during the Holocene and late Pleistocene. Radiocarbon ages have been corrected on the basis of their stable carbon isotope ratios and on environmental samples from the areas, and two groups may be recognised: (1) low ¹⁴C activity groundwaters (13000–34000 yr. BP) with δ ¹³C-5.6 to −11.7‰; and (2) higher ¹⁴C activity groundwaters (5000–7800 yr. BP) enriched in ¹³C up to δ ¹³C = −3.2‰. There is a general correlation of age with depth.

A well defined freshwater (< 50 mg/l Cl⁻) channel can be traced within the aquifer for some 130 km through the region, which is considered to represent recharge from a former wadi. This water with an age of ± 7800 yr. BP is chemically and isotopically distinct from the regional groundwaters and provides direct evidence of a significant recharge event during the Holocene.

The stable isotope (O and H) composition of groundwater from the Kufra and Sirte basins are all related by an evaporative line with slope δ D = 4.5δ ¹³O − 35 with an intercept on the meteoric line of -11‰. This suggests a recharge source continuing into the Holocene from air masses, analogous to current heavy monsoon rain derived from south of the Sahara. The spatial and temporal distribution of groundwaters in relation to the evaporative line suggests a progressive change in character of the recharge which is controlled by a shift towards strongly convective rainfall during the Holocene.

The direct hydrogeological and geochemical evidence supports climatic models proposed by several workers in which discrete humid episodes during the Holocene are inferred.     

Artificial recharge through a thick, heterogeneous unsaturated zone

Thick, heterogeneous unsaturated zones away from large streams in desert areas have not previously been considered suitable for artificial recharge from ponds. To test the potential for recharge in these settings, 1.3 x 10(6) m(3) of water was infiltrated through a 0.36-ha pond along Oro Grande Wash near Victorville, California, between October 2002 and January 2006. The pond overlies a regional pumping depression 117 m below land surface and is located where thickness and permeability of unsaturated deposits allowed infiltration and saturated alluvial deposits were sufficiently permeable to allow recovery of water. Because large changes in water levels caused by nearby pumping would obscure arrival of water at the water table, downward movement of water was measured using sensors in the unsaturated zone. The downward rate of water movement was initially as high as 6 m/d and decreased with depth to 0.07 m/d; the initial time to reach the water table was 3 years. After the unsaturated zone was wetted, water reached the water table in 1 year. Soluble salts and nitrate moved readily with the infiltrated water, whereas arsenic and chromium were less mobile. Numerical simulations done using the computer program TOUGH2 duplicated the downward rate of water movement, accumulation of water on perched zones, and its arrival at the water table. Assuming 10 x 10(6) m(3) of recharge annually for 20 years, a regional ground water flow model predicted water level rises of 30 m beneath the ponds, and rises exceeding 3 m in most wells serving the nearby urban area.     

A recharge model for high altitude,arid, Andean aquifers

Evidence for groundwater recharge in arid zones is mounting, despite early ideas that recharge was unlikely where evaporation greatly exceeded precipitation. The mechanisms and magnitude of groundwater recharge in the Andes and Atacama Desert are not well known but the subject of current research. Diffuse recharge is expected to be limited to high altitude areas with coarse‐grained soils devoid of vegetation. A recharge model for this environment is developed based on a simple soil moisture budgeting technique and the calculation of actual evaporation based on empirical studies. The model is run with data for the Linzor basins, over 4000 m elevation at 22·2°S on the west slope of the Andes. It is checked against independent estimates based on the chloride mass balance (CMB) method and flood events measured downstream in the Río Salado and found to provide robust and reliable results. The results indicate that irregular and volumetrically limited amounts of diffuse recharge occur at high elevations in half of all years, with a tendency to cluster during La Niña episodes. For the Linzor Basins, mean annual recharge is found to be equivalent to 28 mm a^?1, although no recharge occurs in years with precipitation less than 120 mm, and increases proportionately with annual rainfall amounts above this limit. Copyright © 2009 John Wiley & Sons, Ltd.     

Assessment of recharge to aquifers / Evaluation de recharge d'aquifères

Abstract

This paper describes a research project based in Cyprus to investigate a number of different methods for determining quantitative recharge to aquifers in semiarid areas. Large lysimeters and geochemical profile techniques have been used in the southeastern Mesaoria area of Cyprus where there has already been a quantitative water resources investigation and where a groundwater model is in preparation. Further joint geochemical and lysimeter studies are also being undertaken in the Akrotiri peninsula (average annual rainfall 430 mm). Excavation and piling techniques have been used to install the lysimeters which are up to 100 m² in area by 4 m deep. The lysimeters have been made as large as possible to minimize edge effects whilst maintaining the soil and aquifer material within them in an ‘undisturbed’ state. A dry drilling method has been developed to obtain sand samples from the unsaturated zone together with a method for processing the core material to obtain data on moisture and solutes. Chloride profiles have been interpreted to provide estimates of the direct recharge component using a steady state, mass balance approach and initial results compare favourably with tritium estimates. The preliminary results for recharge obtained by the different techniques are compared with each other and with estimates obtained from conventional methods. Initial results from the lysimeters indicate a current annual recharge of 5 mm in southeastern Mesaoria (1978), whilst a mean annual recharge of 50 mm at Akrotiri has been obtained by the chloride profile method.     

Groundwater recharge by channel infiltration in El Barbon basin,Baja California,Mexico

The amount of groundwater recharge by channel infiltration is estimated for El Barbon basin, in Baja California, Mexico. The basin’s lower portion includes the valleys of Ojos Negros and Real del Castillo Viejo, which are crossed by several ephemeral washes, including the mainstem El Barbon Wash. A distributed catchment model with the capability for nonlinear channel routing and channel abstraction is used to calculate groundwater recharge by channel infiltration for storm events of 2-, 5-, 10-, 25-, 50-, and 100-yr return period. The results confirm that event channel infiltration can be a substantial component of the vertical recharge.     

Groundwater recharge amidst focused stormwater infiltration

Distributed, infiltration‐based approaches to stormwater management are being implemented to mitigate effects of urban development on water resources. One of the goals of this type of storm water management, sometimes called low impact development or green infrastructure, is to maintain groundwater recharge and stream base flow at predevelopment levels. However, the connection between infiltration‐based stormwater management and groundwater recharge is not straightforward. Water infiltrated through stormwater facilities may be stored in soil moisture, taken up by evapotranspiration or contribute to recharge and eventually base flow. This study focused on a 1.1 km² suburban, low impact development watershed in Clarksburg, Maryland, USA, that was urbanized and contained 73 infiltration‐based stormwater facilities. Continuous water table measurements were used to quantify the movement of infiltrated stormwater. Time series analyses were performed on hydrographs of 7 wells, and the episodic master recession method was used. Persistence in water levels, as measured by autocorrelation function, was found to be positively related to depth to water. Storm properties (precipitation rate and duration) and well location (proximity to the nearest stream) were significant in driving episodic recharge to precipitation ratios. The well that had the highest recharge to precipitation ratios and water table rises of up to 1.5 m in response to storm events was located furthest from the stream and down gradient of stormwater infiltration locations. This work may be considered in evaluating the effects of planned watershed‐scale infiltration‐based stormwater management on groundwater flow systems.     

Groundwater recharge: A hydrogeologic thought

The ensuing paper imparts vital information on an important component of the hydrologic cycle: recharge. Water flows through the porous media and forms a definite flow pattern that can be depicted in an elaborate manner through the micro level studies at the small watershed level. The estimation of recharge is indispensable for the groundwater budgeting studies. The advantages and disadvantages of some of the techniques have been reviewed. In the present paper an attempt is made to develop a fundamental understanding of the spatial and temporal distribution of the recharge component and to attract excellent research in the same field.     

Geochemical quantification of semiarid mountain recharge

Analysis of a typical semiarid mountain system recharge (MSR) setting demonstrates that geochemical tracers help resolve the location, rate, and seasonality of recharge as well as ground water flowpaths and residence times. MSR is defined as the recharge at the mountain front that dominates many semiarid basins plus the often-overlooked recharge through the mountain block that may be a significant ground water resource; thus, geochemical measurements that integrate signals from all flowpaths are advantageous. Ground water fluxes determined from carbon-14 ((14)C) age gradients imply MSR rates between 2 x 10(6) and 9 x 10(6) m(3)/year in the Upper San Pedro Basin, Arizona, USA. This estimated range is within an order of magnitude of, but lower than, prior independent estimates. Stable isotopic signatures indicate that MSR has a 65% +/- 25% contribution from winter precipitation and a 35% +/- 25% contribution from summer precipitation. Chloride and stable isotope results confirm that transpiration is the dominant component of evapotranspiration (ET) in the basin with typical loss of more than 90% of precipitation-less runoff to ET. Such geochemical constraints can be used to further refine hydrogeologic models in similar high-elevation relief basins and can provide practical first estimates of MSR rates for basins lacking extensive prior hydrogeologic measurements.     

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

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

An artificial recharge experiment in the San Jacinto basin, Riverside, southern California

During a 3 month recharge experiment related to conjunctive use of water resources, 1.5 × 10⁶ T of imported water were percolated through a pond of 128 m by 128 m in the San Jacinto basin. The infiltration rates, which declined with time, averaged 1.9 m day⁻¹, equivalent to four times the lowest laboratory-measured hydraulic conductivity of the fluvial deposits. Ponding altered the unimodal grain-size distribution at the ground surface to types without a dominant mode, but this redistribution did not always lead to reduction in conductivities, which varied over at least three orders of magnitude. The water table 80 m downstream from the ponding edge began to rise slowly 1 month after the start of ponding; it leveled off at 8 m above the pre-recharge water table depth of 75 m and did not recede 2 months after termination of ponding. Water levels in wells bottomed in the original vadose zone suggested that an inverted water table migrated downward to meet the rising water table. Minor, local perching occurred at 14 m depth, as indicated by the presence of moist ground near one monitoring well and by hydraulic responses during a 20 day intermission in percolation. As it percolated through the sediments, the imported northern California water gained Ca but lost Mg, so that the Mg/Ca ratio resembled that of local ground water. Such cation exchange has also been demonstrated by leaching experiments in the laboratory. However, the characteristics of the original source waters appear to be retained by D/H isotope ratios and Cl concentrations, as well as cross-plots of SO₄ vs. Cl and B vs. Cl. Such unreactive tracers could serve to monitor transport and mixing of the chemically diverse water used in future recharge programs in the San Jacinto basin.     

Groundwater recharge assessment in a rural,arid, mid-mountain basin in North-Central Chile

ABSTRACT

A basic component of any hydrogeological study is the magnitude and temporal variation of groundwater recharge. This can be difficult to assess accurately, particularly in arid and semi-arid rainfed mid-mountain zones, as is the situation in the rural, low population density zones of North-Central Chile. In this study, recharge in the Punitaqui Basin, North-Central Chile, was characterized, contrasting the results of two methods: a modified Thornthwaite-Mather (MTM) and discharge recession analysis (DRA). We found a recharge rate of between 1 and 4% of average annual precipitation. Average recharge estimated by the MTM method is consistently higher than that estimated by DRA. Also, DRA tends to smooth the recharge values, resulting in a lower inter-annual variation coefficient. Both methods identified a threshold value of total annual precipitation, above which recharge can be expected to occur, of the order of 180 mm year^?1, consistent with values reported in similar areas.     

Estimating groundwater recharge in lowland watersheds

Little attention has been given to the role of groundwater in the hydrological cycle of lowland watersheds. Our objective in this study was to estimate total recharge to groundwater by analysing water table response to storm events and the rate at which water was transferred into the shallow aquifer. This was conducted at three sites in a rural watershed in the lower Atlantic coastal plain near Charleston, South Carolina, USA. A novel version of the water table fluctuation method was used to estimate total recharge to the shallow aquifer by comparing hourly data of water table position following storm events and measuring water table recession behavior, rather than subjective graphical analysis methods. Also, shallow aquifer recharge rates (vertical fluxes) were estimated using Darcy's Law by comparing static water levels in a water table well and in a shallow piezometer during dry periods. The total annual recharge estimated ranged from 107 ± 39 mm·yr^–1 (5–10% of annual precipitation) at a poorly drained topographic low area to 1140 ± 230 mm·yr^–1 (62–94% of annual precipitation) for a moderately well‐drained upland site. The average aquifer recharge rate was 114 ± 60 mm·yr^–1, which is similar to previous estimations of base flow for the ephemeral third‐order streams in this watershed. The difference in the two methods may have been caused by processes not accounted for in the Darcy flux method, soil moisture deficits, and average evapotranspiration demand, which is about 1100 mm·yr^–1 for this region. Although other factors also can affect partitioning of recharge, an integrated approach to inspecting easily gathered groundwater data can provide information on an often neglected aspect of water budget estimation. We also discuss the effects of land use change on recharge reduction, given a typical development scenario for the region. Copyright © 2011 John Wiley & Sons, Ltd.     

Groundwater recharge and discharge in a hyperarid alluvial plain (Akesu,Taklimakan Desert,China)

Groundwater recharge and discharge in the Akesu alluvial plain were estimated using a water balance method. The Akesu alluvial plain (4842 km²) is an oasis located in the hyperarid Tarim River basin of central Asia. The land along the Akesu River has a long history of agricultural development and the irrigation area is highly dependent on water withdrawals from the river. We present a water balance methodology to describe (a) surface water and groundwater interaction and (b) groundwater interaction between irrigated and non‐irrigated areas. Groundwater is recharged from the irrigation system and discharged in the non‐irrigated area. Uncultivated vegetation and wetlands are supplied from groundwater in the hyperarid environment. Results show that about 90% of groundwater recharge came from canal loss and field infiltration. The groundwater flow from irrigated to non‐irrigated areas was about 70% of non‐irrigated area recharge and acted as subsurface drainage for the irrigation area. This desalinated the irrigation area and supplied water to the non‐irrigated area. Salt moved to the non‐irrigation area following subsurface drainage. We conclude that the flooding of the Akesu River is a supplemental groundwater replenishment mechanism: the river desalinates the alluvial plain by recharging fresh water in summer and draining saline regeneration water in winter. Copyright © 2007 John Wiley & Sons, Ltd.     

Well recharge in idealized rectangular aquifers

Analytical solutions for groundwater flow in rectangular aquifers are presented in the case of a single-well recharge. The problem concerns the impact of a seasonal recharge scheme of variable duration on aquifers with various boundary conditions. The results obtained from these idealized aquifers can be used in a preliminary assessment of the groundwater response to artificial recharge schemes.     

Eruption episodes and magma recharge events in andesitic systems: Mt Taranaki,New Zealand

Eruption episodes, where a series of eruption events are generically related, can include the eruption of a wide spectrum of volcanic activity over decadal periods. This paper concentrates on the opening phases of an eruption episode which occurred approximately 1800 yrs BP from Mt Taranaki, New Zealand. These events spanned the eruption of differing bulk compositions and styles from two distinct vent locations; an andesitic sub-plinian eruption from the summit vent and a scoria cone-building eruption of basaltic magma from a satellite vent. Compositional profiles and zoning textures of plagioclase, amphibole and clinopyroxene phenocrysts from the opening andesitic event show evidence of magma mixing and subsequent crystallisation just prior to the initiation of the eruption episode. Titanomagnetite grain morphology and Ti variation suggest that the magma mixing event occurred within a few days to weeks before the eruption acting as a trigger for it. We present a magmatic model which is constrained by the petrological observations and eruptions of the episode. In this model magma differentiation at depth causes its rise and recharging of a mid-crustal magma storage area at 5–7 km. Although the recharging magma differed slightly in oxygen fugacity and temperature, it was compositionally and physically similar enough to the residing andesitic magma to allow efficient mixing. The petrological characteristics described here can be readily observed and enable identification of mixing events in other recent eruption episodes.     

Spatial and temporal infiltration dynamics during managed aquifer recharge

Natural groundwater recharge is inherently difficult to quantify and predict, largely because it comprises a series of processes that are spatially distributed and temporally variable. Infiltration ponds used for managed aquifer recharge (MAR) provide an opportunity to quantify recharge processes across multiple scales under semi-controlled conditions. We instrumented a 3-ha MAR infiltration pond to measure and compare infiltration patterns determined using whole-pond and point-specific methods. Whole-pond infiltration was determined by closing a transient water budget (accounting for inputs, outputs, and changes in storage), whereas point-specific infiltration rates were determined using heat as a tracer and time series analysis at eight locations in the base of the pond. Whole-pond infiltration, normalized for wetted area, rose rapidly to more than 1.0 m/d at the start of MAR operations (increasing as pond stage rose), was sustained at high rates for the next 40 d, and then decreased to less than 0.1 m/d by the end of the recharge season. Point-specific infiltration rates indicated high spatial and temporal variability, with the mean of measured values generally being lower than rates indicated by whole-pond calculations. Colocated measurements of head gradients within saturated soils below the pond were combined with infiltration rates to calculate soil hydraulic conductivity. Observations indicate a brief period of increasing saturated hydraulic conductivity, followed by a decrease of one to two orders of magnitude during the next 50 to 75 d. Locations indicating the most rapid infiltration shifted laterally during MAR operation, and we suggest that infiltration may function as a "variable source area" processes, conceptually similar to catchment runoff.     

Geochemical evidence of hydrothermal recharge in Lake Baringo,central Kenya Rift Valley

Lake Baringo, a freshwater lake in the central Kenya Rift Valley, is fed by perennial and ephemeral rivers, direct rainfall, and hot springs on Ol Kokwe Island near the centre of the lake. The lake has no surface outlet, but despite high evaporation rates it maintains dilute waters by subsurface seepage through permeable sediments and faulted lavas. New geochemical analyses (major ions, trace elements) of the river, lake, and hot spring waters and the suspended sediments have been made to determine the main controls of lake water quality. The results show that evaporative concentration and the binary mixing between two end members (rivers and thermal waters) can explain the hydrochemistry of the lake waters. Two zones are recognized from water composition. The southern part of the lake near sites of perennial river inflow is weakly influenced by evaporation, has low total dissolved species (TDS), and has a seasonally variable load of mainly detrital suspended sediments. In contrast, waters of the northern part of the lake show evidence for strong evaporation (TDS of up to eight times inflow). Authigenic clay minerals and calcite may be precipitating from those more concentrated fluids. The subaerial hot‐spring waters have a distinctive chemistry and are enriched in some elements that are also present in the lake water. Comparison of the chemical composition of the inflowing surface waters and lake water shows (1) an enrichment of some species (HCO₃^?, Cl, SO₄^2?, F, Na, B, V, Cr, As, Mo, Ba and U) in the lake, (2) a depletion in SiO₂ in the lake, and (3) a possible hydrothermal origin for most F. The rare earth element distribution and the F/Cl and Na/Cl ratios give valuable information on the rate of mixing of the river and hydrothermal fluids in the lake water. Calculations imply that thermal fluids may be seeping upward locally into the lake through grid‐faulted lavas, particularly south of Ol Kokwe Island. Copyright © 2006 John Wiley & Sons, Ltd.