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.     

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.     

Effects of variations of river stage and hydraulic conductivity on temporal scaling of groundwater levels: numerical simulations

It was found in previous studies that groundwater levels may fluctuate as a temporal fractal. In this study numerical simulations of groundwater level fluctuations in an unconfined aquifer near a river were conducted to investigate the effects of aquifer heterogeneity and river stage variations on the fractal behavior of the water levels, h(t). Groundwater recharge was taken to be a white-noise process. The aquifer heterogeneity was simulated with a second-order stationary field of hydraulic conductivity (K) with an exponential variogram model. The results showed that groundwater levels fluctuate as a temporal fractal in both homogeneous and heterogeneous aquifers as long as K is less than 10 m/d. Most aquifers may indeed act as a fractal filter which takes a random non-fractal recharge inputs and produces a fractal responses of groundwater level fluctuations. A crossover in temporal scaling of h(t) may appear in more permeable aquifers. Fluctuations of the groundwater level in a homogeneous aquifer are dominated by the recharge process when the river stage is constant or by the river stage variations when the river stage varies in highly permeable aquifers. Heterogeneity plays an important role in the temporal scaling of h(t) in more permeable aquifers: the stronger the heterogeneity, the stronger the temporal scaling of h(t).     

Proglacial groundwater storage dynamics under climate change and glacier retreat

Proglacial aquifers are an important water store in glacierised mountain catchments that supplement meltwater-fed river flows and support freshwater ecosystems. Climate change and glacier retreat will perturb water storage in these aquifers, yet the climate-glacier-groundwater response cascade has rarely been studied and remains poorly understood. This study implements an integrated modelling approach that combines distributed glacio-hydrological and groundwater models with climate change projections to evaluate the evolution of groundwater storage dynamics and surface-groundwater exchanges in a temperate, glacierised catchment in Iceland. Focused infiltration along the meltwater-fed Virkisá River channel is found to be an important source of groundwater recharge and is projected to provide 14%–20% of total groundwater recharge by the 2080s. The simulations highlight a mechanism by which glacier retreat could inhibit river recharge in the future due to the loss of diurnal melt cycling in the runoff hydrograph. However, the evolution of proglacial groundwater level dynamics show considerable resilience to changes in river recharge and, instead, are driven by changes in the magnitude and seasonal timing of diffuse recharge from year-round rainfall. The majority of scenarios simulate an overall reduction in groundwater levels with a maximum 30-day average groundwater level reduction of 1 m. The simulations replicate observational studies of baseflow to the river, where up to 15% of the 30-day average river flow comes from groundwater outside of the melt season. This is forecast to reduce to 3%–8% by the 2080s due to increased contributions from rainfall and meltwater runoff. During the melt season, groundwater will continue to contribute 1%–3% of river flow despite significant reductions in meltwater runoff inputs. Therefore it is concluded that, in the proglacial region, groundwater will continue to provide only limited buffering of river flows as the glacier retreats.     

Periodic recharge to finite aquifiers from rectangular areas

Analytical solutions for groundwater flow in unconfined rectangular aquifers are presented for the case of recharge from rectangular areas. The linearised differential equation of the flow is solved using Laplace and finite Fourier transforms. The problem concerns the response of finite aquifers to periodic (seasonal) recharge schemes of variable duration. Results can be obtained for various combinations of Dirichlet and Neumann boundary conditions and can be easily used in a preliminary study of water resources management.     

Numerical modeling groundwater recharge and its implication in water cycles of two interdunal valleys in the Sand Hills of Nebraska

The topography and geomorphology of the sand dunes and interdunal valleys in the Nebraska Sand Hills play important roles in regional water cycle by influencing groundwater recharge and evapotranspiration (ET). In this study, groundwater recharge, associated with precipitation and ET as well as soil hydraulics, and its spatial variations owing to the topography of dunes and valleys are examined. A method is developed to describe the recharge as a function of the storage capacity of dunes of various heights. After the method is tested using observations from a network of wells in the Sand Hills, it is used in the MODFLOW model to simulate and describe recharge effects on groundwater table depth at two different dune-valley sites. Analysis of modeled groundwater budget shows that the groundwater table depth in the interdunal valleys is critically influenced by vertical groundwater flows from surrounding dunes. At the site of higher dunes there are steadier and larger vertical groundwater flows in the dunes from their previous storage of precipitation. These vertical flows change to be horizontal converging groundwater flows and create upwelling in the interdunal valleys, where larger ET loss at the surface further enhances groundwater upwelling. Such interdunal valley is the major concentration area of the surface water and groundwater flow in the Sand Hills. At the site of shallow dunes and a broad interdunal valley the supply of groundwater from the dunes is trivial and inadequate to support upwelling of groundwater in the valley. The groundwater flows downward in the valley, and the valley surface is dry. Weak ET loss at the surface has a smaller effect on the groundwater storage than the precipitation recharge, making such area a source for groundwater.     

Using Domestic Well Records to Determine Fractured Bedrock Watersheds and Recharge Rates

This study presents an approach for delineating groundwater basins and estimating rates of recharge to fractured crystalline bedrock. It entailed the use of completion report data (boring logs) from 2500 domestic wells in bedrock from the Coventry Quadrangle, which is located in northeastern Connecticut and characterized by metamorphic gneiss and schist. Completion report data were digitized and imported into ArcGIS^® for data analysis. The data were processed to delineate groundwater drainage basins for the fractured rock based on flow conditions and to estimate groundwater recharge to the bedrock. Results indicate that drainage basins derived from surface topography, in general, may not correspond with bedrock drainage basins due to scale. Estimates of recharge to the bedrock for the study area indicate that only a small fraction of the precipitation or the amount of water that enters the overburden recharges the rock. The approach presented here can be a useful method for water resource‐related assessments that involve fractured rock aquifers.     

Closed‐Form Approximations for Two‐Dimensional Groundwater Age Patterns in a Fresh Water Lens

Simple closed‐form approximations are presented for calculating the steady‐state groundwater age distribution in two‐dimensional vertical cross sections of idealized fresh water lenses overlying salt water, for aquifers that are vertically semi‐infinite and of finite thickness. The approximations are developed on the basis of existing one‐dimensional analytical solutions for travel‐time calculation in fresh water lenses and approximate streamline formulations. The two‐dimensional age distributions based on the closed‐form solutions match convincingly with numerical simulations. As expected, notable deviations from the numerical solution are encountered at the groundwater flow divide and when submarine groundwater discharge occurs. Ratios of recharge over hydraulic conductivities are varied to explore how the magnitude of the deviations changes, and it is found that the approximate closed‐form solutions perform well over a range of conditions found in natural systems.     

Stream depletion in alluvial valleys using the SDF semianalytical model

A semianalytical method commonly used for quantifying stream depletion caused by ground water pumping was reviewed for applicability in narrow alluvial aquifers. This stream depletion factor (SDF) method is based on the analytic Glover model, but uses a numerical model-derived input parameter, called the SDF, to partly account for mathematically nonideal conditions such as variable transmissivity and nearby aquifer boundaries. Using the SDF can improve and simplify depletion estimates. However, the method's approximations introduce error that increases with proximity to the impermeable aquifer boundary. This article reviews the history of the method and its assumptions. New stream depletion response curves are presented as functions of well position within bounded aquifers. A simple modification to modeled SDF values is proposed that allows the impermeable boundary to be accounted for with image wells, but without overaccounting for boundary effects that are already reflected in modeled SDFs. It is shown that SDFs for locations closer to the river than to the aquifer boundary do not reflect impermeable-boundary effects, and thus need no modification, and boundary effects in the other portion of the aquifer follow a predictable removable pattern. This method is verified by comparing response curves using modified SDFs with response curves from an extensively calibrated numerical model of a managed ground water recharge site. The modification improves SDF-based stream depletion estimates in bounded aquifers while still benefiting from the additional information contained in SDF maps and retaining their value as standardized references for water rights administration.     

Predicting Aquifer Response Time for Application in Catchment Modeling

It is well established that changes in catchment land use can lead to significant impacts on water resources. Where land‐use changes increase evapotranspiration there is a resultant decrease in groundwater recharge, which in turn decreases groundwater discharge to streams. The response time of changes in groundwater discharge to a change in recharge is a key aspect of predicting impacts of land‐use change on catchment water yield. Predicting these impacts across the large catchments relevant to water resource planning can require the estimation of groundwater response times from hundreds of aquifers. At this scale, detailed site‐specific measured data are often absent, and available spatial data are limited. While numerical models can be applied, there is little advantage if there are no detailed data to parameterize them. Simple analytical methods are useful in this situation, as they allow the variability in groundwater response to be incorporated into catchment hydrological models, with minimal modeling overhead. This paper describes an analytical model which has been developed to capture some of the features of real, sloping aquifer systems. The derived groundwater response timescale can be used to parameterize a groundwater discharge function, allowing groundwater response to be predicted in relation to different broad catchment characteristics at a level of complexity which matches the available data. The results from the analytical model are compared to published field data and numerical model results, and provide an approach with broad application to inform water resource planning in other large, data‐scarce catchments.     

Influence of Seasonal Variations in Sea Level on the Salinity Regime of a Coastal Groundwater–Fed Wetland

Seasonal variations in sea level are often neglected in studies of coastal aquifers; however, they may have important controls on processes such as submarine groundwater discharge, sea water intrusion, and groundwater discharge to coastal springs and wetlands. We investigated seasonal variations in salinity in a groundwater‐fed coastal wetland (the RAMSAR listed Piccaninnie Ponds in South Australia) and found that salinity peaked during winter, coincident with seasonal sea level peaks. Closer examination of salinity variations revealed a relationship between changes in sea level and changes in salinity, indicating that sea level–driven movement of the fresh water‐sea water interface influences the salinity of discharging groundwater in the wetland. Moreover, the seasonal control of sea level on wetland salinity seems to override the influence of seasonal recharge. A two‐dimensional variable density model helped validate this conceptual model of coastal groundwater discharge by showing that fluctuations in groundwater salinity in a coastal aquifer can be driven by a seasonal coastal boundary condition in spite of seasonal recharge/discharge dynamics. Because seasonal variations in sea level and coastal wetlands are ubiquitous throughout the world, these findings have important implications for monitoring and management of coastal groundwater–dependent ecosystems.     

Temporal variation of stable isotopes in a precipitation–groundwater system: implications for determining the mechanism of groundwater recharge in high mountain–hills of the Loess Plateau,China

The groundwater in shallow loess aquifers in high mountain–hills in the western Loess Plateau in China is almost the sole water resource for local residents. However, the question about how the loess groundwater naturally circulates in these high mountain–hills, characterized by low precipitation and high potential evaporation, remains unclear. The objectives of this study are to evaluate the application of hydrogen and oxygen isotopes to (1) examine temporal variations of the isotopic composition of precipitation and shallow groundwater and (2) uncover the mechanism of groundwater recharge in high mountain–hills. Results from 2 years of monitoring data show a difference in the stable isotopes for groundwater and local precipitation between the winter and summer periods. Similar to precipitation, stable isotopes in groundwater are observed to be depleted in winter and enriched in summer, particularly in oxygen isotope. A prominent characteristic is that H and O isotopes of groundwater show a very clear response to strong precipitation in the rainy season in 2013. The results highlight that local precipitation is the likely recharge source for groundwater in shallow loess aquifers. Annual recharge from local precipitation maintains the groundwater resource in the shallower loess aquifer. The mechanisms governing shallow loess groundwater recharge in high mountain–hills were evaluated. In addition to possible vertical slow percolation of soil water through the unsaturated zone, rapid groundwater recharge mechanisms have been identified as temporal preferential infiltration through sinkholes, slip surface or landslide surface and through the interface of loess layer and palaeo‐soils. Most groundwater can be recharged after a heavy rainy season. Copyright © 2015 John Wiley & Sons, Ltd.     

Assessment of groundwater recharge processes through karst vadose zone by cave percolation monitoring

Recharge processes of karst aquifers are difficult to assess given their strong heterogeneity and the poorly known effect of vadose zone on infiltration. However, recharge assessment is crucial for the evaluation of groundwater resources. Moreover, the vulnerability of karst aquifers depends on vadose zone behaviour because it is the place where most contamination takes place. In this work, an in situ experimental approach was performed to identify and quantify flow and storage processes occurring in karst vadose zone. Cave percolation monitoring and dye tracing were used to investigate unsaturated zone hydrological processes. Two flow components (diffuse and quick) were identified and, respectively, account for 66% and 34% of the recharge. Quickflow was found to be the result of bypass phenomenon in vadose zone related to water saturation. We identify the role of epikarst as a shunting area, most of the storage in the vadose zone occurring via the diffuse flow component in low permeability zones. Relationship between rainfall intensity and transit velocity was demonstrated, with 5 times higher velocities for the quick recharge mode than the diffuse mode. Modelling approach with KarstMod software allowed to simulate the hybrid recharge through vadose zone and shows promising chances to properly assess the recharge processes in karst aquifer based on simple physical models.     

A low‐dimensional hillslope‐based catchment model for layered groundwater flow

Despite the strong interaction between surface and subsurface waters, groundwater flow representation is often oversimplified in hydrological models. For instance, the interplay between local or shallow aquifers and deeper regional‐scale aquifers is typically neglected. In this work, a novel hillslope‐based catchment model for the simulation of combined shallow and deep groundwater flow is presented. The model consists of the hillslope‐storage Boussinesq (hsB) model representing shallow groundwater flow and an analytic element (AE) model representing deep regional groundwater flow. The component models are iteratively coupled via a leakage term based on Darcy's law, representing delayed recharge to the regional aquifer through a low conductivity layer. Simulations on synthetic single hillslopes and on a two‐hillslope open‐book catchment are presented, and the results are compared against a benchmark three‐dimensional Richards equation model. The impact of hydraulic conductivity, hillslope plan geometry (uniform, convergent, divergent), and hillslope inclination (0.2%, 5%, and 30%) under drainage and recharge conditions are examined. On the single hillslopes, good matches for heads, hydrographs, and exchange fluxes are generally obtained, with the most significant differences in outflows and heads observed for the 30% slope and for hillslopes with convergent geometry. On the open‐book catchment, cumulative outflows are overestimated by 1–4%. Heads in the confined and unconfined aquifers are adequately reproduced throughout the catchment, whereas exchange fluxes are found to be very sensitive to the hillslope drainable porosity. The new model is highly efficient computationally compared to the benchmark model. The coupled hsB/AE model represents an alternative to commonly used groundwater flow representations in hydrological models, of particular appeal when surface–subsurface exchanges, local aquifer–regional aquifer interactions, and low flows play a key role in a watershed's dynamics. Copyright © 2011 John Wiley & Sons, Ltd.     

Estimating groundwater recharge from water isotope (δ2H, δ18O) depth profiles in the Densu River basin,Ghana

Abstract

Accurate estimation of groundwater recharge is essential for the proper management of aquifers. A study of water isotope (δ²H, δ¹⁸O) depth profiles was carried out to estimate groundwater recharge in the Densu River basin in Ghana, at three chosen observation sites that differ in their altitude, geology, climate and vegetation. Water isotopes and water contents were analysed with depth to determine water flow in the unsaturated zone. The measured data showed isotope enrichment in the pore water near the soil surface due to evaporation. Seasonal variations in the isotope signal of the pore water were also observed to a depth of 2.75 m. Below that depth, the seasonal variation of the isotope signal was attenuated due to diffusion/dispersion and low water flow velocities. Groundwater recharge rates were determined by numerical modelling of the unsaturated water flow and water isotope transport. Different groundwater recharge rates were computed at the three observation sites and were found to vary between 94 and 182 mm/year (± max. 7%). Further, the approximate peak-shift method was applied to give information about groundwater recharge rates. Although this simple method neglects variations in flow conditions and only considers advective transport, it yielded mean groundwater recharge rates of 110–250 mm/year (± max. 30%), which were in the same order of magnitude as computed numerical modelling values. Integrating these site-specific groundwater recharge rates to the whole catchment indicates that more water is potentially renewed than consumed nowadays. With increases in population and irrigation, more clean water is required, and knowledge about groundwater recharge rates – essential for improving the groundwater management in the Densu River basin – can be easily obtained by measuring water isotope depth profiles and applying a simple peak-shift approach.

Citation Adomako, D., Maloszewski, P., Stumpp, C., Osae, S. & Akiti, T. T. (2010) Estimating groundwater recharge from water isotope (δ²H, δ¹⁸O) depth profiles in the Densu River basin, Ghana. Hydrol. Sci. J. 55(8), 1405–1416.     

Application of an integrated surface water‐groundwater model to multi‐aquifers modeling in Choushui River alluvial fan,Taiwan

This research proposes a combination of SWAT and MODFLOW, MD‐SWAT‐MODFLOW, to address the multi‐aquifers condition in Choushui River alluvial fan, Taiwan. The natural recharge and unidentified pumping/recharge are separately estimated. The model identifies the monthly pumping/recharge rates in multi‐aquifers so that the daily streamflow can be simulated correctly. A multi‐aquifers condition means a subsurface formation composed of at least the unconfined aquifer, the confined aquifer, and an in‐between aquitard. In such a case, the variation of groundwater level is related to pumping/recharge activities in vertically adjacent aquifer and the river‐aquifer interaction. Both factors in turn affect the streamflow performance. Results show that MD‐SWAT‐MODFLOW performs better than SWAT alone in terms of simulated streamflow, especially during low flow period, when pumping/recharge rates are properly estimated. A sensitivity analysis of individual parameter suggests that the vertical leakance may be the most sensitive among all investigated MODFLOW parameters in terms of the estimated pumping/recharge among aquifers, and the Latin‐Hypercube‐One‐factor‐At‐a‐Time sensitivity analysis indicates that the hydraulic conductivity of channel is the most sensitive to the model performance. It also points out the necessity to simultaneously estimate pumping/recharge rates in multi‐aquifers. The estimated net pumping rate can be treated as a lower bound of the actual local pumping rate. As a whole, the model provides the spatio‐temporal groundwater use, which gives the authorities insights to manage groundwater resources. Copyright © 2012 John Wiley & Sons, Ltd.     

Is there a superior conceptual groundwater model structure for baseflow simulation?

Previous work has shown that streamflow response during baseflow conditions is a function of storage, but also that this functional relationship varies among seasons and catchments. Traditionally, hydrological models incorporate conceptual groundwater models consisting of linear or non‐linear storage–outflow functions. Identification of the right model structure and model parameterization however is challenging. The aim of this paper is to systematically test different model structures in a set of catchments where different aquifer types govern baseflow generation processes. Nine different two‐parameter conceptual groundwater models are applied with multi‐objective calibration to transform two different groundwater recharge series derived from a soil‐atmosphere‐vegetation transfer model into baseflow separated from streamflow data. The relative performance differences of the model structures allow to systematically improve the understanding of baseflow generation processes and to identify most appropriate model structures for different aquifer types. We found more versatile and more aquifer‐specific optimal model structures and elucidate the role of interflow, flow paths, recharge regimes and partially contributing storages. Aquifer‐specific recommendations of storage models were found for fractured and karstic aquifers, whereas large storage capacities blur the identification of superior model structures for complex and porous aquifers. A model performance matrix is presented, which highlights the joint effects of different recharge inputs, calibration criteria, model structures and aquifer types. The matrix is a guidance to improve groundwater model structures towards their representation of the dominant baseflow generation processes of specific aquifer types. Copyright © 2014 John Wiley & Sons, Ltd.     

Analytic element modelling for strip aquifers

The Analytic Element Method (AEM) provides a convenient tool for groundwater flow analysis in unbounded continuous domains. The AEM is based on the superposition of analytic functions, known as elements, useful at both regional and local scales. In this study, analytic elements for strip aquifers are presented. Such aquifers occur in riverine or coastal deposits and in outcrop zones of confined aquifers. Local flow field is modelled indirectly, using a reference plane related to the aquifer domain through the Schwarz‐Christoffel transform. The regional flow is obtained as a solution of the one‐dimensional flow equation. The proposed methodology was tested by modelling two hypothetical situations, which were compared to exact solutions. It is shown that regional boundaries can be reproduced exactly while local fields are adequately reproduced with analytic elements. The developed elements are applied to simulate a real flow field in northeastern Brazil showing good agreement with measured water levels. Copyright © 2011 John Wiley & Sons, Ltd.     

Modification of the SWAT model to simulate regional groundwater flow using a multicell aquifer

The soil and water assessment tool (SWAT) has been widely used and thoroughly tested in many places in the world. The application of the SWAT model has pointed out that 2 of the major weaknesses of SWAT are related to the nonspatial reference of the hydrologic response unit concept and to the simplified groundwater concept, which contribute to its low performance in baseflow simulation and its inability to simulate regional groundwater flow. This study modified the groundwater module of SWAT to overcome the above limitations. The modified groundwater module has 2 aquifers. The local aquifer, which is the shallow aquifer in the original SWAT, represents a local groundwater flow system. The regional aquifer, which replaces the deep aquifer of the original SWAT, represents intermediate and regional groundwater flow systems. Groundwater recharge is partitioned into local and regional aquifer recharges. The regional aquifer is represented by a multicell aquifer (MCA) model. The regional aquifer is discretized into cells using the Thiessen polygon method, where centres of the cells are locations of groundwater observation wells. Groundwater flow between cells is modelled using Darcy's law. Return flow from cell to stream is conceptualized using a non‐linear storage–discharge relationship. The SWAT model with the modified aquifer module, the so‐called SWAT‐MCA, was tested in 2 basins (Wipperau and Neetze) with porous aquifers in a lowland area in Lower Saxony, Germany. Results from the Wipperau basin show that the SWAT‐MCA model is able (a) to simulate baseflow in a lowland area (where baseflow is a dominant source of streamflow) better than the original model and (b) to simulate regional groundwater flow, shown by the simulated groundwater levels in cells, quite well.     

Using stable water isotopes to identify spatio-temporal controls on groundwater recharge in two contrasting East African aquifer systems

Understanding the spatio-temporal variability in groundwater recharge is a prerequisite to sustainable management of aquifers. Spatial analysis of groundwater stable isotopes uncovered predominant controls on groundwater recharge in the Nairobi aquifer system (NAS) and the South Coast aquifer (SC), two exemplar East African aquifers relied upon by 7 million people. We analysed 368 samples for stable isotopes and basic physico-chemical parameters. The NAS groundwater isotopes are controlled by precipitation orographic effects and enriched recharge from impounded lakes/wetlands; the SC isotopes are correlated with water-table depth influencing evapotranspiration. Global Network of Isotopes in Precipitation (GNIP) data revealed groundwater recharge during months of heavy rains in the NAS, whilst the SC experiences spatio-temporally diffuse recharge. Inferred “isoscapes” show: in the NAS, (1) direct, rapid recharge favoured by faults, well-drained soils and ample rainfall in uplands, (2) delayed recharge from impounded lakes and wetlands in mid-lands, and (3) focused, event-based recharge in floodplains; and in the SC, diffuse recharge complicated by significant water-table evapotranspiration processes.