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.     

CMBEAR: Python-Based Recharge Estimator Using the Chloride Mass Balance Method in Australia

The chloride mass balance (CMB) method is widely used to estimate long-term rates of groundwater recharge. In regions where surface water runoff is negligible, recharge can be estimated using measurements of chloride concentrations of groundwater and precipitation, and an estimate of long-term average rainfall. This paper presents the Chloride Mass Balance Estimator of Australian Recharge (CMBEAR), a Jupyter (Python) Notebook that is set up to rapidly apply the CMB method using gridded maps of chloride deposition rates across the Australian continent. For an Australian context, the chloride deposition rate and rainfall maps have been provided. Thus, CMBEAR requires only a spreadsheet with the groundwater chloride concentration, the latitude and longitude of the sample location, and some simple user inputs. CMBEAR may be easily applied in other regions, providing that a gridded chloride deposition map is available. Recharge estimates from CMBEAR are compared against published applications of the CMB method. CMBEAR is also applied to a large dataset from the Northern Territory and is used to produce a gridded map of recharge for western Victoria. CMBEAR provides a reproducible and straightforward approach to apply the CMB method to estimate groundwater recharge.     

Regional estimation of total recharge to ground water in Nebraska

Naturally occurring long-term mean annual recharge to ground water in Nebraska was estimated by a novel water-balance approach. This approach uses geographic information systems (GIS) layers of land cover, elevation of land and ground water surfaces, base recharge, and the recharge potential in combination with monthly climatic data. Long-term mean recharge > 140 mm per year was estimated in eastern Nebraska, having the highest annual precipitation rates within the state, along the Elkhorn, Platte, Missouri, and Big Nemaha River valleys where ground water is very close to the surface. Similarly high recharge values were obtained for the Sand Hills sections of the North and Middle Loup, as well as Cedar River and Beaver Creek valleys due to high infiltration rates of the sandy soil in the area. The westernmost and southwesternmost parts of the state were estimated to typically receive < 30 mm of recharge a year.     

Comparative performance of soil water balance models in computing semi-arid aquifer recharge

Abstract

Estimating groundwater recharge is essential to ensure the sustainable use of groundwater resources, particularly in arid and semi-arid regions. Soil water balances have been frequently advocated as valuable tools to estimate groundwater recharge. This article compares the performance of three soil water balance models (Hydrobal, Visual Balan v2.0 and Thornthwaite) in the Ventós-Castellar aquifer, Spain. The models were used to simulate wet and dry years. Recharge estimates were transformed into water table fluctuations by means of a lumped groundwater model. These, in turn, were calibrated against piezometric data. Overall, the Hydrobal model shows the best fit between observed and calculated levels (r² = 0.84), highlighting the role of soil moisture and vegetation in recharge processes.

Editor D. Koutsoyiannis; Associate editor X. Chen

Citation Touhami, I., et al., 2014. Comparative performance of soil water balance models in computing semi-arid aquifer recharge. Hydrological Sciences Journal, 59 (1), 193–203.     

Influence of Seasonal Pumping on Groundwater Sources and Flow System,Nagaoka Plain,Japan

Intensive groundwater development in the urban area of the Nagaoka Plain, Japan, has induced changes in the pH and saturation index of calcite in groundwater. To account for these chemical changes, it is important to determine seasonal variations of recharge and the groundwater flow system in the aquifer. This study identified the sources and flow system of groundwater in this urban area by a comprehensive method using stable isotope data and a numerical groundwater model of the Nagaoka Plain. Stable isotope evidence shows that the groundwater is recharged by meteoric water originating from low‐elevation areas rather than the mountains surrounding the plain. The water table in the study area is drawn down during the winter and recovers in the other seasons. Numerical modeling shows that discharge occurs primarily along the Shinano River during the recovery period, whereas discharge is centered in urbanized areas during the drawdown period, when a conical depression of the water table stimulates recharge from the immediate area. These results are indications of a local groundwater flow system, with its recharge area between the Shinano River and the urban areas, which is governed by intensive seasonal groundwater extraction.     

Spatio‐temporal analysis of groundwater recharge and mound dynamics in an unconfined aquifer: a GIS‐based approach

Groundwater recharge and mounding of water‐table is a complex phenomenon involving time‐ and space‐dependent hydrologic processes. The effect of long‐term groundwater mounding in the aquifer depends on soil, aquifer geometry and the area contributing to recharge. In this paper, a GIS‐based spatio‐temporal algorithm has been developed for the groundwater mound dynamics to estimate the potential rise in the water‐table and groundwater volume balance residual in an unconfined aquifer. The recharge and mound dynamics as predicted using the methodology recommended here were compared with those using the Hantush equation, and the differences were quite significant. The significance of the study is to assess the effectiveness of the basin in terms of its hydrologic and hydraulic properties for sustainable management of groundwater recharge. Copyright © 2007 John Wiley & Sons, Ltd.     

Hydrodynamic and isotopic investigations for evaluating the mechanisms and amount of groundwater seepage through a rockslide dam

This study addresses the influence of landslide dams on surface water drainage and groundwater flow. In the study area of Scanno Lake and Sagittario River (Central Italy), a limestone rockslide‐avalanche formed a lake, which has an outlet that is occasionally active, showing infiltration into the rockslide dam. Several springs are present at the lake's base and are partly fed by seepage through the rockslide debris. Piezometric surveys, discharge measurements, pumping tests and chemical analyses are tools used to build a conceptual model of the groundwater flow and to evaluate the flow through the rockslide debris. Seasonal water isotopic signatures validate the assumed model, showing a mixing of infiltration recharge and groundwater seepage throughout the rockslide debris. Various recharge areas have been found for springs, pointing out those directly fed by the rockslide debris aquifer. Hypotheses about seasonal groundwater mixing between the regional carbonate aquifer and the rockslide debris aquifer are supported by isotope results. Seasonal changes in groundwater table level due to recharge and surface losses from seasonal outlet have been correlated with isotopic groundwater composition from the rockslide debris aquifer and the downstream springs; this relationship highlights the role of the rockslide dam body on the hydrodynamics of the studied area. Relationships between surface waters and groundwater in the area have been completely understood on the basis of water isotopic fingerprinting, finally obtaining a complete evaluation of groundwater renewable resources and its regimen. Copyright © 2010 John Wiley & Sons, Ltd.     

MODIS‐Aided Statewide Net Groundwater‐Recharge Estimation in Nebraska

Monthly evapotranspiration (ET) rates (2000 to 2009) across Nebraska at about 1‐km resolution were obtained by linear transformations of the MODIS (MODerate resolution Imaging Spectroradiometer) daytime surface temperature values with the help of the Priestley–Taylor equation and the complementary relationship of evaporation. For positive values of the mean annual precipitation and ET differences, the mean annual net recharge was found by an additional multiplication of the power‐function‐transformed groundwater vulnerability DRASTIC‐code values. Statewide mean annual net recharge became about 29 mm (i.e., 5% of mean annual precipitation) with the largest recharge rates (in excess of 100 mm/year) found in the eastern Sand Hills and eastern Nebraska. Areas with the largest negative net recharge rates caused by declining groundwater levels due to large‐scale irrigation are found in the south‐western region of the state. Error bounds of the estimated values are within 10% to 15% of the corresponding precipitation rates and the estimated net recharge rates are sensitive to errors in the precipitation and ET values. This study largely confirms earlier base‐flow analysis‐based statewide groundwater recharge estimates when considerations are made for differences in the recharge definitions. The current approach not only provides better spatial resolution than available earlier studies for the region but also quantifies negative net recharge rates that become especially important in numerical modeling of shallow groundwater systems.     

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.     

Distinguishing Mountain Front and Mountain Block Recharge in an Intermontane Basin of the Himalayan Region

Intermontane basin aquifers worldwide, particularly in the Himalayan region, are recharged largely by the adjoining mountains. Recharge in these basins can occur either by water infiltrating from streams near mountain fronts (MFs) as mountain front recharge (MFR) or by sub-surface mountain block infiltration as mountain block recharge (MBR). MFR and MBR recharge are challenging to distinguish and are least quantified, considering the lack of extensive understanding of the hydrological processes in the mountains. This study used oxygen and hydrogen isotopes (δ¹⁸O and δ²H), electrical conductivity (EC) data, hydraulic head, and groundwater level data to differentiate MFR and MBR. Groundwater level data provide information about the groundwater-surface water interactions and groundwater flow directions, whereas isotopes and EC data are used to distinguish and quantify different recharge sources. The present methodology is tested in an intermontane basin of the Himalayan region. The results suggest that karst springs (KS) and deep groundwater (DGW) recharge are dominated by snowmelt (47% ± 10% and 46% ± 9%) as MBR from adjacent mountains, insignificantly affected by evaporation. The hydraulic head data and isotopes indicate Quaternary shallow groundwater (SGW) aquifer system recharge as MFR of local meteoric water with significant evaporation. The results indicate several flow paths in the aquifer system, a local flow for KS, intermediate flow for SGW, and regional flow for DGW. The findings will significantly impact water resource management in the area and provide vital baseline knowledge for sustainable groundwater management in other Himalayan intermontane basins.     

Evidence of short-term groundwater recharge signal propagation from the Andes to the central Atacama Desert: a singular spectrum analysis approach

The effect of pressure-driven groundwater recharge signal propagation in the Andean-Atacamenian environment is investigated by assessing a record of 15 years of water table fluctuations of an unconfined–confined aquifer system. Based on a singular spectrum analysis of water table time series, it is shown that, in the given case, groundwater levels in the central Atacama Desert are hydraulically controlled by two distant recharge areas associated with the Andes. The maximum observed range of the pressure signal propagation is ~50 km over an elevation difference of more than 3000 m at a lag of ~25 months. Several findings indicate that an often-cited study misinterpreted a water level rise at the same site as an in-situ alluvial fan recharge. Thus, the effect’s impact on groundwater dynamics in complex aquifer systems can easily be overlooked. Singular spectrum analysis could be of use to investigate pressure effects at hydrologically comparable sites.     

Estimating recharge using relations between precipitation and yield in a mountainous area with large variability in precipitation

Estimates of recharge to bedrock aquifers from infiltration of precipitation can be difficult to obtain, especially in areas with large spatial and temporal variability in precipitation. In the Black Hills area of western South Dakota and eastern Wyoming, streamflow yield is highly influenced by annual precipitation, with yield efficiency (annual yield divided by annual precipitation) increasing with increasing annual precipitation. Spatial variability in annual yield characteristics for Black Hills streams is predictably influenced by precipitation patterns. Relations between precipitation and yield efficiency were used to estimate annual recharge from long-term records of annual precipitation. A series of geographic information system algorithms was used to derive annual estimates for 1000- by 1000-m grid cells. These algorithms were composited to derive estimates of annual recharge rates to the Madison and Minnelusa aquifers in the Black Hills area of western South Dakota and eastern Wyoming during water years 1931–1998 and an estimate of average recharge for water years 1950–1998. This approach provides a systematic method of obtaining consistent and reproducible estimates of recharge from infiltration of precipitation. Resulting estimates of average annual recharge (water years 1950–1998) ranged from 1 cm in the southern Black Hills to 22 cm in the northwestern Black Hills. Recharge rates to these aquifers from infiltration of precipitation on outcrops was estimated to range from 0.9 m³/s in 1936 to 18.8 m³/s in 1995.     

Application of grey correlation method to evaluate potential groundwater recharge sites

Artificial recharge is a practical tool available for increasing the groundwater storage capacity. The efficiency of artificial recharge is related to various hydrogeological factors of the target area. In this study, a variable saturated groundwater flow model, FEMWATER, was used to evaluate the arrival times of recharged water that infiltrates from an artificial recharge pond to the groundwater table under various hydrogeological conditions. Forty-five arrival times were generated by FEMWATER. The relationships between the arrival times and hydrogeological factors used in the simulation of FEMWATER were analyzed by the grey correlation method. The results show the order of importance of the factors as they influence the arrival time. In order from high to low importance, they are α, D _g, θ _e, D _p, K _S and β. D _g and D _p are interpreted as the potential for movement of the recharge water; θ_e is the water storage capacity of soil, and K _S represents the ability of soil to transport water. α and β describe the characteristic curve of the unsaturated soil. The method was applied to evaluate a suitable site for artificial recharge in the Yun-Lin area. Grey correlation analysis was performed to obtain the grey correlation grade using the minimum arrival time as a reference sequence. An index is proposed herein to determine the recharge efficiency of 20 sampling sites. A contour mapping of index values at the 20 sampling sites identified three areas for artificial aquifer recharge in Yun-Lin. Area A in the upper plain is considered more appropriate for groundwater recharge than areas B and C in the coast.     

Dynamics of a headwater system and peatland under current conditions and with climate change

Interactions between headwater aquifers and peatlands have received limited scientific attention. Hydrological stresses, including those related to climate change, may adversely impact these interactions. In this study, the dynamics of a southern Québec headwater system where a peatland is present is simulated under current conditions and with climate change. The model is calibrated in steady state on field‐measured data and provides satisfactory results for transient‐state conditions. Under current conditions, simulations confirm that the peatland is fed by the fractured bedrock aquifer year‐round and provides continuous baseflow to its outlets. Climate change is simulated through its impact on groundwater recharge. Predicted precipitation and temperature data from a suite of regional climate model scenarios provide a net precipitation variation range from +10% to ?30% for the 2041–2070 horizon. Calibrated recharge is modified within this range to perform a sensitivity analysis of the headwater model to recharge variations (+10%, ?15% and ?30%). Total contribution from the aquifer to rivers and streams varies from +14% to ?44% of the baseline for +10% to ?30% recharge changes from spring 2010 data, for example. With higher recharge, the peatland receives more groundwater, which could significantly change its vegetation pattern and eventually ecosystem functions. For a ?30% recharge, the peatland becomes perched above the aquifer during the summer, fall and winter. Recharge reductions also induce sharp declines in groundwater levels and drying streams. Copyright © 2013 John Wiley & Sons, Ltd.     

Noble gas tracers for characterisation of flow dynamics and origin of groundwater: A case study in Switzerland

The Grenchen aquifer system in the Swiss Plateau was extensively investigated in order to determine the extent of groundwater contamination and to assess the natural attenuation capacity. Environmental tracer data were applied to estimate groundwater travel times, mixing ratios, and evaluate groundwater origin. Recharge is basically possible in two distinct topographical areas, the immediate vicinity of the town of Grenchen and the elevated plateau of the first Jura Mountain ridge. Groundwater dating was performed with the ³H/³He dating method and supplemented by ⁸⁵Kr measurements. Stable isotope data (δ¹⁸O, δ²H) and dissolved noble gas concentrations allow the determination of the recharge temperature, which is correlated to the recharge elevation. Noble gas temperatures (NGT) decrease in the direction of groundwater flow and range from 10 to 13 °C in the upstream area of the town to 7–9 °C in the downstream river plain. This trend could suggest the admixture of water from the underlying limestone aquifer recharged under cooler infiltration conditions, e.g. at higher recharge elevations. However, it is shown in this study that the difference in NGT does not require such a recharge. Rather, increasing air temperatures over the last 40 years and the urban heat island effect could possibly explain most of the observed temperature shift. Furthermore, it is concluded that the downstream river plain is hydrologically disconnected from the upstream town area. Consequently most water from the town area is drained by the creek Witibach and recharge in the river plain is higher than previously assumed.     

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.     

Soil profile evolution following land‐use change: implications for groundwater quantity and quality

Soil and vadose zone profiles are used as an archive of changes in groundwater recharge and water quality following changes in land use in an area of the Loess Plateau of China. A typical rain‐fed loess‐terrace agriculture region in Hequan, Guyuan, is taken as an example, and multiple tracers (chloride mass balance, stable isotopes, tritium and water chemistry) are used to examine groundwater recharge mechanisms and to evaluate soil water chloride as an archive for recharge rate and water quality. Results show that groundwater recharge beneath natural uncultivated grassland, used as a baseline, is about 94–100 mm year^?1 and that the time it takes for annual precipitation to reach water table through the thick unsaturated zone is from decades to hundreds of years (tritium free). This recharge rate is 2–3 orders of magnitude more than in the other semiarid areas with similar annual rainfall but with deep‐rooted vegetation and relatively high temperature. Most of the water that eventually becomes recharge originally infiltrated in the summer months. The conversion from native grassland to winter wheat has reduced groundwater recharge by 42–50% (50–55 mm year^?1 for recharge), and the conversion from winter wheat to alfalfa resulted in a significant chloride accumulation in the upper soil zone, which terminated deep drainage. The paper also evaluates the time lag between potential recharge and actual recharge to aquifer and between increase in solute concentration in soil moisture and that in the aquifer following land‐use change due to the deep unsaturated zone. Copyright © 2012 John Wiley & Sons, Ltd.     

Interactions between perched and saprolite aquifers on a small,salt-affected and deeply weathered hillslope

A small hillslope was chosen to investigate the role of throughflow as a mechanism responsible for the movement of soil water and solutes towards a saline seep and as a source of recharge to a permanent, regional aquifer at depth. The hydraulic properties, chemical characteristics and physical responses of both systems were studied on a deeply weathered, salt-affected hillslope. Additional data were also obtained from other sites in south-western Australia. Regional groundwater flow occurred in a variably textured, deeply weathered material in which the hydraulic conductivity varied from < 0·001 to 0·14m day^?1. Perched groundwater flow (throughflow) occurred in the higher permeability (? 1·5 m day^?1), near-surface soil materials. Throughflow occurred throughout winter, contributing approximately 530 m³ of fresh (? 160 mg l^?1 Cl) water to a saline seep. By contrast, the deep aquifer discharged approximately 1100 m³ of waters with salt concentrations of 2000–6000 mg l^?1 Cl. Recharge and discharge rates to and from the deep aquifer, were estimated to be of the order of 5–20 mm a^?1 and 50–300 mm a^?1 respectively. Saturated conditions existed throughout winter within the seep and the immediately adjacent non-saline area, with up to 60 per cent of the hillslope soils becoming saturated after major rainfall events ( > 20 mm day^?1). In the mid-slopes, in particular along a central depression, saturation of the shallow soils caused macropore channel recharge to take place through the clay-textured subsoils. Water-level responses suggest that approximately 25–30 per cent of annual recharge occurred from one storm studied in September 1984. Recharge through macropore channels is a significant mechanism in the concave slope areas on the hillslope. Throughflow was found to be a major source of water, but not salt, contributing to the saline seep. In general, the contribution of throughflow was found to decrease further inland at other sites studied. However, at inland sites where perennial, perched aquifers have developed in deep sands, saline areas have been caused by throughflow and not by deep aquifer discharge.     

Pervasive,Preferential Flow through Mega-Thick Unsaturated Zones in the Southern Great Basin

Recharge from preferential flow through mega-thick (100–1000 m) unsaturated zones is a pervasive phenomenon, as demonstrated with a case study of volcanic highland recharge areas in the Great Basin province in southern Nevada, USA. Statistically significant rising water-level trends occur for most study-area wells and resulted from a relatively wet period (1969–2005) in south-central Nevada. Wet and dry winters control water-level trends, with water levels rising within a few months to a year following a wet-winter recharge event and declining during sustained dry periods. Even though a megadrought has persisted since 2000, this drought condition did not preclude major recharge events. Modern groundwater reaching the water table is consistent with previous geochemical studies of the study area that indicate mixing of modern and late Pleistocene recharge water. First-order approximations and simple mixing models of modern and late Pleistocene water indicate that 10% to 40% of recharge is preferential flow and that modern recharge may play a larger role in the water budget than previously thought.