Geostatistical mapping of leakance in a regional aquitard, Oak Ridges Moraine area, Ontario, Canada

The Newmarket Till forms a regionally extensive aquitard separating two major aquifer systems in the Greater Toronto area, Canada. The till is incised, and sometimes eroded entirely, by a network of sand- and gravel-filled channels forming productive aquifers and, locally, high-conductivity windows between aquifer systems. Leakage through the till may also be substantial in places. This study investigates the spatial variability of aquitard leakance in order to assess the relative importance of recharge processes to the lower aquifers. With a large database derived from water-well records and containing both hard and soft information, the Sequential Indicator Simulation method is used to generate maps of aquitard thickness and window probability. These can be used for targeting channel aquifers and for identifying potential areas of recharge to the lower aquifers. Conductivities are modeled from sparse data assuming that their correlation range is much smaller than the grid spacing. Block-scale leakances are obtained by upscaling nodal values based on simulated conductivity and thickness fields. Under the "aquifer-flow' assumption, upscaling is performed by arithmetic spatial averaging. Histograms and maps of upscaled leakances show that heterogeneities associated with aquitard windows have the largest effect on regional groundwater flow patterns. Electronic Publication     

Estimation of groundwater recharge using a GIS-based distributed water balance model in Dire Dawa, Ethiopia

Sustainable groundwater management requires knowledge of recharge. Recharge is also an important parameter in groundwater flow and transport models. Spatial variation in recharge due to distributed land-us.e, soil texture, topography, groundwater level, and hydrometeorological conditions should be accounted for in recharge estimation. However, conventional point-estimates of recharge are not easily extrapolated or regionalized. In this study, a spatially distributed water balance model WetSpass was used to simulate long-term average recharge using land-use, soil texture, topography, and hydrometeorological parameters in Dire Dawa, a semiarid region of Ethiopia. WetSpass is a physically based methodology for estimation of the long-term average spatial distribution of surface runoff, actual evapotranspiration, and groundwater recharge. The long-term temporal and spatial average annual rainfall of 626 mm was distributed as: surface runoff of 126 mm (20%), evapotranspiration of 468 mm (75%), and recharge of 28 mm (5%). This recharge corresponds to 817 l/s for the 920.12 km² study area, which is less than the often-assumed 1,000 l/s recharge for the Dire Dawa groundwater catchment.     

Hydrogeochemical and Flow Modeling of Aquitard Percolation in the Cariri Valley-Northeast Brazil

Flow modeling and hydrogeochemical inversion calculations were performed on the aquifer systems of the Crato-Juazeiro Graben in the Araripe Sedimentary Basin of the Cariri Valley, Northeast Brazil, in order to assess aquifer interconnections through an aquitard of the Araripe Plateau. Steady state hydraulic head numeric calculations with the finite difference method used Modflow software with hydrologically based boundary conditions and a piezometric surface as initial condition. Hydrogeochemical inversion with PHREEQC was performed for two water samples representing typical initial and final waters along the flow path. Flow simulation shows that approximately 20% of all recharge water from the superior aquifer system passes through the Santana aquitard reaching the middle aquifer system. Quantification of mineral dissolution/precipitation processes, ion exchange and microbiological redox reactions taking place in the aquifer systems also indicates hydraulic interconnection of the superior and middle aquifer systems through the Santana aquitard.     

Provenance of buried esker groundwater: the case of Vars-Winchester esker aquifer,Eastern Ontario,Canada

An innovative mode of groundwater recharge to a buried esker aquifer is considered. The current conceptual model affords a natural safeguard to underlying aquifers from the overlying muds. A hypothesis of groundwater recharge to a buried esker aquifer via preferential pathways across its overlying muds is tested here by heuristic numerical one-dimensional and two-dimensional modeling simulations. The hypothesis has been tested against two other conventionally accepted scenarios involving: (1) distal esker outcrop areas and (2) remote shallow-bedrock recharge areas. The main evidence comes from documented recharge pressure pulses in the overlying mud aquitard and in the underlying esker hydraulic-head time series for the Vars-Winchester esker aquifer in Eastern Ontario, Canada. These perturbations to the potentiometric surface are believed to be the aquifer response to recharge events. The migration rate of these pressure pulses is directly related to the hydraulic diffusivity of the formation. The measured response time and response amplitude between singular radar precipitation events and well hydrographs constituted the heuristic model calibration targets. The main evidence also includes mud-layering deformation (water escape features) which was observed in seismic surveys of the over-esker muds. These disturbed stratigraphic elements provide a realistic mechanism for migrating water to transit through the muds. The effective hydraulic conductivities of these preferential pathways in the muds were estimated to be between 2?×?10^?6 and 7?×?10^?6 m/s. The implications of these findings relate to the alleged natural safeguard of these overlying muds.     

Estimation of groundwater recharge in a major sand and gravel aquifer in Ireland using multiple approaches

Groundwater recharge was investigated in the most extensive sand and gravel aquifer (area of approximately 200 km²) in the Republic of Ireland as part of a wider study seeking to derive recharge estimates using aquifer vulnerability mapping. The proportion of effective rainfall (total rainfall minus actual evapotranspiration) that leads to recharge is known as the recharge coefficient. The recharge investigation involved a variety of approaches, including soil moisture budgeting, well hydrograph analysis, numerical modelling and a catchment water balance. The adoption of multiple techniques provided insights on recharge and also on aquifer properties. Comparison of two soil moisture budgeting approaches (FAO Penman-Monteith with Penman-Grindley) showed how variations in the effective rainfall values from these methods influence groundwater levels simulated in a numerical groundwater model. The catchment water balance estimated the recharge coefficient to be between 81 and 85%, which is considered a reasonable range for this aquifer, where overland flow is rarely observed. The well hydrograph analysis, using a previous estimate of specific yield (0.13), gave recharge coefficients in the range of 40–80%, considered low for this aquifer: a revised specific yield of 0.19 resulted in a more reasonable range of recharge coefficients of between 70 and 100%.     

Aquifer system response to intensive pumping in urban areas of the Gangetic plains,India: the case study of Patna

A significant component of domestic demand for water of urban areas located in the Gangetic plains is met by heavy pumping of groundwater. The present study is focused on the Patna municipal area, inhabited by 17 million people and spanning over 134 km², where entire urban water demand is catered from pumping by wells of various capacities and designs. The present study examines the nature of the aquifer system within the urban area, the temporal changes in the water/piezometric level and the recharge mechanism of the deeper aquifers. The aquifer system is made up of medium-to-coarse unconsolidated sand, lying under a ~40-m-thick predominantly argillaceous unit holding 8- to 13-m-thick localised sand layers and continues up to 220 m below ground. Groundwater occurs under semi-confined condition, with transmissivity of aquifers in 5,500–9,200 m² day^?1 range. Hydraulic head of the deeper aquifer remains in 9–19 m range below ground, in contrast to 1–9 m range of that of the upper aquitard zone. The estimated annual groundwater extraction from the deeper aquifer is ~212.0 million m³, which has created a decline of 3.9 m in the piezometric level of the deeper aquifer during the past 30 years. Unregulated construction of deep tube wells with mushrooming of apartment culture may further exacerbate the problem. The sand layers within the aquitard zone are experiencing an annual extraction of 14.5 million m³ and have exhibited stable water level trend for past one and half decades. This unit is recharged from monsoon rainfall, besides contribution from water supply pipe line leakage and seepage from unlined storm water drains.     

Numerical groundwater-flow modeling to evaluate potential effects of pumping and recharge: implications for sustainable groundwater management in the Mahanadi delta region,India

Process-based groundwater models are useful to understand complex aquifer systems and make predictions about their response to hydrological changes. A conceptual model for evaluating responses to environmental changes is presented, considering the hydrogeologic framework, flow processes, aquifer hydraulic properties, boundary conditions, and sources and sinks of the groundwater system. Based on this conceptual model, a quasi-three-dimensional transient groundwater flow model was designed using MODFLOW to simulate the groundwater system of Mahanadi River delta, eastern India. The model was constructed in the context of an upper unconfined aquifer and lower confined aquifer, separated by an aquitard. Hydraulic heads of 13 shallow wells and 11 deep wells were used to calibrate transient groundwater conditions during 1997–2006, followed by validation (2007–2011). The aquifer and aquitard hydraulic properties were obtained by pumping tests and were calibrated along with the rainfall recharge. The statistical and graphical performance indicators suggested a reasonably good simulation of groundwater flow over the study area. Sensitivity analysis revealed that groundwater level is most sensitive to the hydraulic conductivities of both the aquifers, followed by vertical hydraulic conductivity of the confining layer. The calibrated model was then employed to explore groundwater-flow dynamics in response to changes in pumping and recharge conditions. The simulation results indicate that pumping has a substantial effect on the confined aquifer flow regime as compared to the unconfined aquifer. The results and insights from this study have important implications for other regional groundwater modeling studies, especially in multi-layered aquifer systems.     

Groundwater level fluctuation in the Waimea Plains, New Zealand: changes in a coastal aquifer within the last 30 years

Data for the Waimea Plains, New Zealand indicate that the lower confined groundwater aquifer is hydraulically homogeneous and that shallow groundwater levels inland are affected mostly by anthropogenic processes, while those near the coast are affected more by sea level variation. Analysis of long-term data for New Zealand indicates that sea level has increased continuously, but trends are not spatially uniform. Results from non-parametric trend analysis show that rising trends for groundwater levels are predominant in the shallow aquifer both inland on the Waimea Plains and, for recent years, near the coast, while decreasing trends are evident in the underlying confined aquifer near the coast. Groundwater level change in the shallow aquifer appears to be more affected by climate change than the lower confined aquifer. Correlation analysis indicated that groundwater levels are more affected by rainfall during the rainy season than the dry season and more influenced by rainfall inland than near the coast. Groundwater level declines in the lower confined aquifer near the coast, which has its major recharge area inland in the catchment, may be substantially affected by groundwater abstraction in inland areas as well as sea level variation, but there are little evidences of seawater intrusion. Meanwhile, groundwater recharge over the catchment area has great influence on rising groundwater levels in the shallow aquifer and its recharge is estimated to be 417.8 mm/year using chloride concentrations of precipitation and groundwater.     

Hydrogeology and hydrochemistry of the Parmelia aquifer,northern Perth Basin,Western Australia

This study aims to characterise the hydrogeology and hydrochemistry of the Parmelia aquifer and to understand controls on recent water-level changes as these are needed to underpin a quantitative analysis of recharge. The Parmelia aquifer, a layered sequence of sand, silt and discontinuous lenses of clay, receives diffuse rainfall recharge on its outcrop and groundwater recharge occurs across the Dandaragan Plateau at different rates. Water levels have risen steadily over the last three decades between 10 and 55 cm/y in response to the replacement of native vegetation with pasture and annual crops. The mean aquifer properties from sediment analyses indicate a very wide range of porosity (8.9 – 49.5 %) with an arithmetic mean of 26% and consequently a very broad range of specific yield (0.0004 – 0.4) with an arithmetic mean of 0.14. Groundwater in the Parmelia aquifer has an underlying meteoric origin with compositional changes due to reactions with silicate minerals and leaching of chloride that has concentrated in the soil by evapotranspiration. The hydrochemistry sampled at different depths and locations in the aquifer indicates that the groundwater is not well mixed, and variations arise due to relatively recent recharge that has undergone evaporation in some areas.     

Delineating spring recharge areas in a fractured sandstone aquifer (Luxembourg) based on pesticide mass balance

A simple method to delineate the recharge areas of a series of springs draining a fractured aquifer is presented. Instead of solving the flow and transport equations, the delineation is reformulated as a mass balance problem assigning arable land in proportion to the pesticide mass discharged annually in a spring at minimum total transport cost. The approach was applied to the Luxembourg Sandstone, a fractured-rock aquifer supplying half of the drinking water for Luxembourg, using the herbicide atrazine. Predictions of the recharge areas were most robust in situations of strong competition by neighbouring springs while the catchment boundaries for isolated springs were extremely sensitive to the parameter controlling flow direction. Validation using a different pesticide showed the best agreement with the simplest model used, whereas using historical crop-rotation data and spatially distributed soil-leaching data did not improve predictions. The whole approach presents the advantage of integrating objectively information on land use and pesticide concentration in spring water into the delineation of groundwater recharge zones in a fractured-rock aquifer.     

Mapping an aquitard breach using shear-wave seismic reflection

In multi-layered hydrostratigraphic systems, aquitard breaches caused by faulting or paleo-erosion can allow substantial quantities of water of differing quality to be exchanged between aquifers. Seismic reflection technology was used to map the extent and orientation of an aquitard breach connecting a shallow alluvial aquifer to the deeper semi-confined Memphis aquifer in southwestern Tennessee, USA. Geophysical well logs indicate the presence of the aquitard at borehole locations that define the beginning and end points on two seismic survey lines, which intersect at a borehole where the aquitard is absent. A SE–NW-oriented paleochannel, 350 m wide and approximately 35–40 m deep, is interpreted from the seismic reflection surveys. The paleochannel cuts through the aquitard and into the upper part of the Memphis aquifer, thus creating a hydraulic connection between the shallow unconfined and deeper, semi-confined aquifers. The results indicate the potential of the shear-wave seismic reflection methods to resolve shallow breaches through fine-grained aquitards given availability of sufficient well control.     

Use of the subsurface thermal regime as a groundwater-flow tracer in the semi-arid western Nile Delta,Egypt

Temperature profiles from 25 boreholes were used to understand the spatial and vertical groundwater flow systems in the Western Nile Delta region of Egypt, as a case study of a semi-arid region. The study area is located between the Nile River and Wadi El Natrun. The recharge areas, which are located in the northeastern and the northwestern parts of the study area, have low subsurface temperatures. The discharge areas, which are located in the western (Wadi El Natrun) and southern (Moghra aquifer) parts of the study area, have higher subsurface temperatures. In the deeper zones, the effects of faults and the recharge area in the northeastern direction disappear at 80 m below sea level. For that depth, one main recharge and one main discharge area are recognized. The recharge area is located to the north in the Quaternary aquifer, and the discharge area is located to the south in the Miocene aquifer. Two-dimensional groundwater-flow and heat-transport models reveal that the sealing faults are the major factor disturbing the regional subsurface thermal regime in the study area. Besides the main recharge and discharge areas, the low permeability of the faults creates local discharge areas in its up-throw side and local recharge areas in its down-throw side. The estimated average linear groundwater velocity in the recharge area is 0.9 mm/day to the eastern direction and 14 mm/day to the northwest. The average linear groundwater discharge velocities range from 0.4 to 0.9 mm/day in the southern part.     

A poly morphological landform approach for hydrogeological applications in heterogeneous glacial sediments

The ‘poly morphological (PM) concept’ is used to assess geological heterogeneity in glacial sediments and to evaluate groundwater recharge to aquifers, especially those underlying clay till. The novelty of the PM concept is that it enhances the knowledge of, especially, heterogeneity of till units in geological models by using superimposed geomorphological units with typical sediment types related to the individual landforms. The PM concept is demonstrated through the construction of a poly morphological map for Zealand in Denmark and an aggregated map comparing clay-till-thickness distribution within individual PM types. The hydrogeological applicability of the PM concept was tested with a physical-based distributed 3D hydrological model. The aggregated PM map was compared with fracture and redox conditions at 21 field sites to evaluate the PM-type correlation to: (1) till thickness, (2) thickness of the reduced-till zone, and (3) depth to the zone with fracture spacing >1 m The results show that the till thickness is a critical parameter for the formation of fracture networks in clay till, and that fracture density decreases in till thicker than 8–10 m. However, the amount of data is still inadequate for statistically proving the PM concept.     

Consideration of the use of direct recharge in analytical models of flood-wave response: a case study of the Merguellil alluvial aquifer,central Tunisia

In a stream–aquifer system, adequate knowledge of the aquifer hydraulic parameters is paramount for the analysis of groundwater/surface-water interaction. In such systems, analytical flood-wave response models are commonly used to assess these parameters. However, when an important means of direct recharge (precipitation) occurs simultaneously with the stream flood waves, the classical flood-wave response model (i.e. when only stream flood waves are considered as an input signal) leads to overestimated and variable values of the aquifer hydraulic parameters. In this study, an analytical flood-wave response model was implemented to determine hydraulic parameters of the Merguellil alluvial aquifer (central Tunisia) taking into account the recharge rate as a second input signal. An array of three piezometers and a river-stage gauge equipped with an automated Diver data logger were installed in the floodplain of the Merguellil Wadi. The results showed a good agreement between the observed and modelled hydrographs when the recharge rate was considered, which allows the model to provide the same value of diffusivity (50–70 m² h^?1) in the three piezometers and for all flood events being tested. However, when the recharge rate was not taken into account, the model gave an overestimated diffusivity with a very poor hydrograph match. Overall, the methodology described in this paper can be applicable to any area in the world because of its simplicity and its better estimation of aquifer hydraulic parameters.     

Modeling approaches and strategies for data-scarce aquifers: example of the Dar es Salaam aquifer in Tanzania

Management of groundwater resources can be improved by using groundwater models to perform risk analyses and to improve development strategies, but a lack of extensive basic data often limits the implementation of sophisticated models. Dar es Salaam in Tanzania is an example of a city where increasing groundwater use in a Pleistocene aquifer is causing groundwater-related problems such as saline intrusion along the coastline, lowering of water-table levels, and contamination of pumping wells. The lack of a water-level monitoring network introduces a problem for basic data collection and model calibration and validation. As a replacement, local water-supply wells were used for measuring groundwater depth, and well-top heights were estimated from a regional digital elevation model to recalculate water depths to hydraulic heads. These were used to draw a regional piezometric map. Hydraulic parameters were estimated from short-time pumping tests in the local wells, but variation in hydraulic conductivity was attributed to uncertainty in well characteristics (information often unavailable) and not to aquifer heterogeneity. A MODFLOW model was calibrated with a homogeneous hydraulic conductivity field and a sensitivity analysis between the conductivity and aquifer recharge showed that average annual recharge will likely be in the range 80–100 mm/year.     

Holocene estuarine sediments as a source of arsenic in Pleistocene groundwater in suburbs of Hanoi,Vietnam

Groundwater pollution by arsenic is a major health threat in suburban areas of Hanoi, Vietnam. The present study evaluates the effect of the sedimentary environments of the Pleistocene and Holocene deposits, and the recharge systems, on the groundwater arsenic pollution in Hanoi suburbs distant from the Red River. At two study sites (Linh Dam and Tai Mo communes), undisturbed soil cores identified a Pleistocene confined aquifer (PCA) and Holocene unconfined aquifer (HUA) as major aquifers, and Holocene estuarine and deltaic sediments as an aquitard layer between the two aquifers. The Holocene estuarine sediments (approximately 25–40 m depth, 9.6–4.8 cal ka BP) contained notably high concentrations of arsenic and organic matter, both likely to have been accumulated by mangroves during the Holocene sea-level highstand. The pore waters in these particular sediments exhibited elevated levels of arsenic and dissolved organic carbon. Arsenic in groundwater was higher in the PCA (25–94 μg/L) than in the HUA (5.2–42 μg/L), in both the monitoring wells and neighboring household tubewells. Elevated arsenic concentration in the PCA groundwater was likely due to vertical infiltration through the arsenic-rich and organic-matter-rich overlying Holocene estuarine sediments, caused by massive groundwater abstraction from the PCA. Countermeasures to prevent arsenic pollution of the PCA groundwater may include seeking alternative water resources, reducing water consumption, and/or appropriate choice of aquifers for groundwater supply.     

Assessment of groundwater inflows into Kuteshwar Limestone Mines through flow modeling study, Madhya Pradesh, India

The development of limestone mining activities in Katni, Madhya Pradesh becomes necessary to increase the depth of exploration to produce ore. Increase in the exploration depth means that mining pits were subjected to water inrush. A hydrological and a hydrogeological model for the Katni area have been developed using USGS flow code, MODFLOW 2000. Collected GIS-based information was synthesized in a finite difference numerical model. The regional steady flow was calibrated under pre-development conditions assuming an equivalent porous medium approach. Water budget calculations show that the total groundwater flow into the aquifer system due to interaction with river amounts to 14,783 m³/day. Infiltration from precipitation provides 1,600 m³/day of the groundwater supply, while 1,446 m³/day comes from lateral inflow and the remaining. The inflows into mine pit area amounts to 15,725 m³/day. Although the karstic nature of the limestone aquifer the equivalent porous medium flow model is appropriate to represent hydraulic heads and recharge/discharge relationships on a regional scale. The results of this study can be used to predict the required amounts of pumping and the possible locations to dewater the groundwater in the mining pits.     

Multiple-methods investigation of recharge at a humid-region fractured rock site,Pennsylvania, USA

Lysimeter-percolate and well-hydrograph analyses were combined to evaluate recharge for the Masser Recharge Site (central Pennsylvania, USA). In humid regions, aquifer recharge through an unconfined low-porosity fractured-rock aquifer can cause large magnitude water-table fluctuations over short time scales. The unsaturated hydraulic characteristics of the subsurface porous media control the magnitude and timing of these fluctuations. Data from multiple sets of lysimeters at the site show a highly seasonal pattern of percolate and exhibit variability due to both installation factors and hydraulic property heterogeneity. Individual event analysis of well hydrograph data reveals the primary influences on water-table response, namely rainfall depth, rainfall intensity, and initial water-table depth. Spatial and seasonal variability in well response is also evident. A new approach for calculating recharge from continuous water-table elevation records using a master recession curve (MRC) is demonstrated. The recharge estimated by the MRC approach when assuming a constant specific yield is seasonal to a lesser degree than the recharge estimate resulting from the lysimeter analysis. Partial reconciliation of the two recharge estimates is achieved by considering a conceptual model of flow processes in the highly-heterogeneous underlying fractured porous medium.     

Recharge characteristics of a phreatic aquifer as determined by storage accumulation

The cumulative storage accumulation curve (CSAC) is a tool for saturated-volume fluctuation (SVF) analysis of transient recharge to shallow phreatic aquifers discharging only to springs. The method assumes that little underflow or phreatic evapotranspiration occurs. The CSAC is a modified water-table hydrograph that distinguishes storage increase caused by recharge from loss due to springflow-induced recession. Required for the analysis are water-table fluctuations at a single representative location within the catchment of a single spring and either direct measurements or robust interpolations of springflows at different aquifer stages. The method employs empirical manipulation of head observations, varying spring catchment area to minimize CSAC water-level changes in late portions of long recessions. Results include volumetric estimates of recharge integrated over individual events and instantaneous rates of recharge to the water table, at the temporal resolution of the water-level sampling interval. The analysis may also yield physically realistic information on spring catchment and recharge focusing. In a test case in West Virginia, USA, recharge estimates by this technique were consistent with integrated springflow time series but greater than estimates based on potential evapotranspiration. Results give insight into dynamic recharge behavior over time as well as an indication of recharge catchment size. Electronic Publication     

Unusual catchment runoff in a high alpine karst environment influenced by a complex geological setting (Northern Calcareous Alps,Tyrol, Austria)

Garber Schlag (Q-GS) is one of the major springs of the Karwendel Mountains, Tyrol, Austria. This spring has a unique runoff pattern that is mainly controlled by the tectonic setting. The main aquifer is a moderately karstified and jointed limestone of the Wetterstein Formation that is underlain by nonkarstified limestone of the Reifling Formation, which acts as an aquitard. The aquifer and aquitard of the catchment of spring Q-GS form a large anticline that is bound by a major fault (aquitard) to the north. Discharge of this spring shows strong seasonal variations with three recharge origins, based on δ¹⁸O and electrical conductivity values. A clear seasonal trend is observed, caused by the continuously changing portions of water derived from snowmelt, rainfall and groundwater. At the onset of the snowmelt period in May, the discharge is composed mainly of groundwater. During the maximum snowmelt period, the water is dominantly composed of water derived from snowmelt and subordinately from rainfall. During July and August, water derived from snowmelt continuously decreases and water derived from rainfall increases. During September and October, the water released at the spring is mainly derived from groundwater and subordinately from rainfall. The distinct discharge plateau from August to December and the following recession until March is likely related to the large regional groundwater body in the fissured and moderately karstified aquifer of the Wetterstein Formation and the tectonic structures (anticline, major fault). Only a small portion of the water released at spring Q-GS is derived from permafrost.