Estimating Groundwater Mounding in Sloping Aquifers for Managed Aquifer Recharge

Design of managed aquifer recharge (MAR) for augmentation of groundwater resources often lacks detailed data, and simple diagnostic tools for evaluation of the water table in a broad range of parameters are needed. In many large‐scale MAR projects, the effect of a regional aquifer base dip cannot be ignored due to the scale of recharge sources (e.g., wadis, streams, reservoirs). However, Hantush's (1967) solution for a horizontal aquifer base is commonly used. To address sloping aquifers, a new closed‐form analytical solution for water table mound accounts for the geometry and orientation of recharge sources at the land surface with respect to the aquifer base dip. The solution, based on the Dupiuit‐Forchheimer approximation, Green's function method, and coordinate transformations is convenient for computing. This solution reveals important MAR traits in variance with Hantush's solution: mounding is limited in time and space; elevation of the mound is strongly affected by the dip angle; and the peak of the mound moves over time. These findings have important practical implications for assessment of various MAR scenarios, including waterlogging potential and determining proper rates of recharge. Computations are illustrated for several characteristic MAR settings.     

Ecohydrologic Process Modeling of Mountain Block Groundwater Recharge

Regional mountain block recharge (MBR) is a key component of alluvial basin aquifer systems typical of the western United States. Yet neither water scientists nor resource managers have a commonly available and reasonably invoked quantitative method to constrain MBR rates. Recent advances in landscape-scale ecohydrologic process modeling offer the possibility that meteorological data and land surface physical and vegetative conditions can be used to generate estimates of MBR. A water balance was generated for a temperate 24,600-ha mountain watershed, elevation 1565 to 3207 m, using the ecosystem process model Biome-BGC (BioGeochemical Cycles) ( Running and Hunt 1993 ). Input data included remotely sensed landscape information and climate data generated with the Mountain Climate Simulator (MT-CLIM) ( Running et al. 1987 ). Estimated mean annual MBR flux into the crystalline bedrock terrain is 99,000 m³/d, or approximately 19% of annual precipitation for the 2003 water year. Controls on MBR predictions include evapotranspiration (radiation limited in wet years and moisture limited in dry years), soil properties, vegetative ecotones (significant at lower elevations), and snowmelt (dominant recharge process). The ecohydrologic model is also used to investigate how climatic and vegetative controls influence recharge dynamics within three elevation zones. The ecohydrologic model proves useful for investigating controls on recharge to mountain blocks as a function of climate and vegetation. Future efforts will need to investigate the uncertainty in the modeled water balance by incorporating an advanced understanding of mountain recharge processes, an ability to simulate those processes at varying scales, and independent approaches to calibrating MBR estimates.     

Modeling Stream-Aquifer Interactions Under Seasonal Groundwater Pumping and Managed Aquifer Recharge

In South Korea, a significant amount of groundwater is used for the heating of water-curtain insulated greenhouses during the winter dry season, which had led to problems of groundwater depletion. A managed aquifer recharge (MAR) project is currently underway with the goal of preventing such groundwater depletion in a typical cultivation area, located on an alluvial aquifer near the Nam River. In the present study, FEFLOW, a three-dimensional finite element model, was used to evaluate different strategies for MAR of the cultivation areas. A conceptual model was developed to simulate the stream-aquifer dynamics under the influence of seasonal groundwater pumping and MAR. The optimal rates and duration of MAR were assessed by analyzing the recovery of the groundwater levels and the change in the groundwater temperature. The simulation results indicate that a MAR rate of 8000 m³/d effectively restores the groundwater level when the injection wells are located inside the groundwater depletion area. It is also demonstrated that starting the MAR before the beginning of the seasonal pumping is more effective. Riverbank filtration is preferable for securing the injection water owing to plentiful source of induced recharge from the river. Locating the pumping wells adjacent to the river where there are thick permeable layers could be a good strategy for minimizing decreases in the groundwater level and temperature.     

Benefits and Cautions in Data Assimilation Strategies: An Example of Modeling Groundwater Recharge

Assimilating recent observations improves model outcomes for real-time assessments of groundwater processes. This is demonstrated in estimating time-varying recharge to a shallow fractured-rock aquifer in response to precipitation. Results from estimating the time-varying water-table altitude (h) and recharge, and their error covariances, are compared for forecasting, filtering, and fixed-lag smoothing (FLS), which are implemented using the Kalman Filter as applied to a data-driven, mechanistic model of recharge. Forecasting uses past observations to predict future states and is the current paradigm in most groundwater modeling investigations; filtering assimilates observations up to the current time to estimate current states; and FLS estimates states following a time lag over which additional observations are collected. Results for forecasting yield a large error covariance relative to the magnitude of the expected recharge. With assimilating recent observations of h, filtering and FLS produce estimates of recharge that better represent time-varying observations of h and reduce uncertainty in comparison to forecasting. Although model outcomes from applying data assimilation through filtering or FLS reduce model uncertainty, they are not necessarily mass conservative, whereas forecasting outcomes are mass conservative. Mass conservative outcomes from forecasting are not necessarily more accurate, because process errors are inherent in any model. Improvements in estimating real-time groundwater conditions that better represent observations need to be weighed for the model application against outcomes with inherent process deficiencies. Results from data assimilation strategies discussed in this investigation are anticipated to be relevant to other groundwater processes models where system states are sensitive to system inputs.     

A Routing Technique for Estimating Ground–Water Recharge

A technique is presented for the estimation of ground–water recharge in a region where the formations overlying an aquifer cannot transmit all of the available water to the aquifer. The method is based on a routing model which is derived from the physical nature of formations overlying the aquifer. The model takes into account precipitation, evapotranspiration, soil moisture balance and the passage of water through strata such as boulder clay. The validity of the model is tested by calculated and field stream flow hydrographs.     

Sulfur Hexafluoride and Potassium Bromide as Groundwater Tracers for Managed Aquifer Recharge

Sulfur hexafluoride (SF₆) is an established tracer for use in managed aquifer recharge projects. SF₆ exsolves from groundwater when it encounters trapped air according to Henry's law. This results in its retardation relative to groundwater flow, which can help determine porous media saturation and flow dynamics. SF₆ and the conservative, nonpartitioning tracer, bromide (Br⁻ added as KBr), were introduced to recharge water infiltrated into stacked glacial aquifers in Thurston County, Washington, providing the opportunity to observe SF₆ partitioning. Br⁻, which is assumed to travel at the same velocity as the groundwater, precedes SF₆ at most monitoring wells (MWs). Average groundwater velocity in the unconfined aquifer in the study area ranges from 3.9 to 40 m/d, except in the southwestern corner where it is slower. SF₆ in the shallow aquifer exhibits an average retardation factor of 2.5 ± 3.8, suggesting an air-to-water ratio on the order of 10⁻³ to 10⁻² in the pore space. Notable differences in tracer arrival times at adjacent wells indicate very heterogeneous conductivity. One MW exhibits double peaks in concentrations of both tracers with different degrees of retardation for the first and second peaks. This suggests multiple flowpaths to the well with variable saturation. The confining layer between the upper two aquifers appears to allow intermittent connection between aquifers but serves as an aquitard in most areas. This study demonstrates the utility of SF₆ partitioning for evaluating hydrologic conditions at prospective recharge sites.     

Groundwater Recharge Rate and Zone Structure Estimation Using PSOLVER Algorithm

The quantification of groundwater recharge is an important but challenging task in groundwater flow modeling because recharge varies spatially and temporally. The goal of this study is to present an innovative methodology to estimate groundwater recharge rates and zone structures for regional groundwater flow models. Here, the unknown recharge field is partitioned into a number of zones using Voronoi Tessellation (VT). The identified zone structure with the recharge rates is associated through a simulation‐optimization model that couples MODFLOW‐2000 and the hybrid PSOLVER optimization algorithm. Applicability of this procedure is tested on a previously developed groundwater flow model of the Tahtal? Watershed. Successive zone structure solutions are obtained in an additive manner and penalty functions are used in the procedure to obtain realistic and plausible solutions. One of these functions constrains the optimization by forcing the sum of recharge rates for the grid cells that coincide with the Tahtal? Watershed area to be equal to the areal recharge rate determined in the previous modeling by a separate precipitation‐runoff model. As a result, a six‐zone structure is selected as the best zone structure that represents the areal recharge distribution. Comparison to results of a previous model for the same study area reveals that the proposed procedure significantly improves model performance with respect to calibration statistics. The proposed identification procedure can be thought of as an effective way to determine the recharge zone structure for groundwater flow models, in particular for situations where tangible information about groundwater recharge distribution does not exist.     

The Safe Yield and Climatic Variability: Implications for Groundwater Management

Methods for calculating the safe yield are evaluated in this paper using a high‐quality and long historical data set of groundwater recharge, discharge, extraction, and precipitation in a karst aquifer. Consideration is given to the role that climatic variability has on the determination of a climatically representative period with which to evaluate the safe yield. The methods employed to estimate the safe yield are consistent with its definition as a long‐term average extraction rate that avoids adverse impacts on groundwater. The safe yield is a useful baseline for groundwater planning; yet, it is herein shown that it is not an operational rule that works well under all climatic conditions. This paper shows that due to the nature of dynamic groundwater processes it may be most appropriate to use an adaptive groundwater management strategy that links groundwater extraction rates to groundwater discharge rates, thus achieving a safe yield that represents an estimated long‐term sustainable yield. An example of the calculation of the safe yield of the Edwards Aquifer (Texas) demonstrates that it is about one‐half of the average annual recharge.     

Characteristics of Groundwater Recharge on the North China Plain

Groundwater recharge is an important component of the groundwater system. On the North China Plain (NCP), groundwater is the main water supply. Because of large‐scale overexploitation, the water table has declined, which has produced severe adverse effects on the environment and ecosystem. In this article, tracer experiment and watershed model were used to calculate and analyze NCP groundwater recharge. In the tracer experiment, average recharge was 108 mm/year and recharge coefficient 0.16. With its improved irrigation, vegetation coverage and evapotranspiration modules, the INFIL3.0 model was used for calculation of groundwater recharge. Regional modeling results showed an average recharge of 102 mm/year and recharge coefficient 0.14, for 2001–2009. These values are very similar to those from the field tracer experiment. Influences in the two methods were analyzed. The results can provide an important reference for NCP groundwater recharge.     

Simple Methodology for Estimating the Groundwater Recharge Potential of Rural Ponds and Lakes Using Remote Sensing and GIS Techniques: a Spatiotemporal Case Study of Roorkee Tehsil,India

Water Resources - Groundwater resources are depleting at an alarming rate all over the world. Therefore, ensuring water security is the foremost priority of every nation. This paper presents a...     

A New Method for Estimating Recharge to Unconfined Aquifers Using Differential River Gauging

In semiarid and arid environments, leakage from rivers is a major source of recharge to underlying unconfined aquifers. Differential river gauging is widely used to estimate the recharge. However, the methods commonly applied are limited in that the temporal resolution is event‐scale or longer. In this paper, a novel method is presented for quantifying both the total recharge volume for an event, and variation in recharge rate during an event from hydrographs recorded at the upstream and downstream ends of a river reach. The proposed method is applied to river hydrographs to illustrate the method steps and investigate recharge processes occurring in a sub‐catchment of the Murray Darling Basin (Australia). Interestingly, although it is the large flood events which are commonly assumed to be the main source of recharge to an aquifer, our analysis revealed that the smaller flow events were more important in providing recharge.     

Recognition of Regional Water Table Patterns for Estimating Recharge Rates in Shallow Aquifers

We propose a new method for groundwater recharge rate estimation in regions with stream-aquifer interactions, at a linear scale on the order of 10 km and more. The method is based on visual identification and quantification of classically recognized water table contour patterns. Simple quantitative analysis of these patterns can be done manually from measurements on a map, or from more complex GIS data extraction and curve fitting. Recharge rate is then estimated from the groundwater table contour parameters, streambed gradients, and aquifer transmissivity using an analytical model for groundwater flow between parallel perennial streams. Recharge estimates were obtained in three regions (areas of 1500, 2200, and 3300 km²) using available water table maps produced by different methods at different times in the area of High Plains Aquifer in Nebraska. One region is located in the largely undeveloped Nebraska Sand Hills area, while the other two regions are located at a transition zone from Sand Hills to loess-covered area and include areas where groundwater is used for irrigation. Obtained recharge rates are consistent with other independent estimates. The approach is useful and robust diagnostic tool for preliminary estimates of recharge rates, evaluation of the quality of groundwater table maps, identification of priority areas for further aquifer characterization and expansion of groundwater monitoring networks prior to using more detailed methods.     

The Effect of Modeled Recharge Distribution on Simulated Groundwater Availability and Capture

Simulating groundwater flow in basin‐fill aquifers of the semiarid southwestern United States commonly requires decisions about how to distribute aquifer recharge. Precipitation can recharge basin‐fill aquifers by direct infiltration and transport through faults and fractures in the high‐elevation areas, by flowing overland through high‐elevation areas to infiltrate at basin‐fill margins along mountain fronts, by flowing overland to infiltrate along ephemeral channels that often traverse basins in the area, or by some combination of these processes. The importance of accurately simulating recharge distributions is a current topic of discussion among hydrologists and water managers in the region, but no comparative study has been performed to analyze the effects of different recharge distributions on groundwater simulations. This study investigates the importance of the distribution of aquifer recharge in simulating regional groundwater flow in basin‐fill aquifers by calibrating a groundwater‐flow model to four different recharge distributions, all with the same total amount of recharge. Similarities are seen in results from steady‐state models for optimized hydraulic conductivity values, fit of simulated to observed hydraulic heads, and composite scaled sensitivities of conductivity parameter zones. Transient simulations with hypothetical storage properties and pumping rates produce similar capture rates and storage change results, but differences are noted in the rate of drawdown at some well locations owing to the differences in optimized hydraulic conductivity. Depending on whether the purpose of the groundwater model is to simulate changes in groundwater levels or changes in storage and capture, the distribution of aquifer recharge may or may not be of primary importance.     

Recharge and Groundwater Flow Within an Intracratonic Basin,Midwestern United States

The conservative nature of chloride (Cl^?) in groundwater and the abundance of geochemical data from various sources (both published and unpublished) provided a means of developing, for the first time, a representation of the hydrogeology of the Illinois Basin on a basin‐wide scale. The creation of Cl^? isocons superimposed on plan view maps of selected formations and on cross sections across the Illinois Basin yielded a conceptual model on a basin‐wide scale of recharge into, groundwater flow within and through the Illinois Basin. The maps and cross sections reveal the infiltration and movement of freshwater into the basin and dilution of brines within various geologic strata occurring at basin margins and along geologic structures. Cross‐formational movement of brines is also seen in the northern part of the basin. The maps and cross sections also show barriers to groundwater movement created by aquitards resulting in areas of apparent isolation/stagnation of concentrated brines within the basin. The distribution of Cl^? within the Illinois Basin suggests that the current chemical composition of groundwater and distribution of brines within the basin is dependent on five parameters: (1) presence of bedrock exposures along basin margins; (2) permeability of geologic strata and their distribution relative to one another; (3) presence or absence of major geologic structures; (4) intersection of major waterways with geologic structures, basin margins, and permeable bedrock exposures; and (5) isolation of brines within the basin due to aquitards, inhomogeneous permeability, and, in the case of the deepest part of the basin, brine density effects.