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.     

Application of WetSpass model to estimate groundwater recharge variability in the Nile Delta aquifer

Long-term groundwater recharge from rainfall in the Nile Delta is needed as an input for integrated groundwater modelling in the Nile Delta aquifer for more accurate simulation. The main objective is to estimate the spatial and temporal variation of groundwater recharge from rainfall in the Nile Delta aquifer. Water and Energy Transfer between Soil, Plants and Atmosphere under quasi-Steady State (WetSpass) model parameters were identified for the Nile Delta based on the available meteorological data for the area collected in 1991 and 2000. The collected data were rainfall, temperature, wind speed and evapotranspiration. Geomorphological characteristics, such as soil type, topography, groundwater depth and slope, were also collected as input data for the WetSpass model. ENVI software was used to come up with land use classification based on available land cover images of the Nile Delta for 1972, 1984, 1990, 2000 and 2009. The WetSpass model was calibrated by comparing the simulated groundwater recharge with the calculated one by using the water balance equation model. The results indicated close agreement in groundwater recharge between the two model outputs with R ² of 0.99 and 0.94, while the root-mean-square errors (RMSEs) were 4.86 and 9.39 mm for 1991 and 2000, respectively. The WetSpass model was then applied in respect of 1970, 1980, 1990 and 2010 for the purpose of validation. The overall RMSE and R ² for the 6 years were 8.83 mm and 0.88, respectively. The results of the WetSpass calibrated model provide information to support integrated groundwater modelling. The results reveal that WetSpass works well in simulating the components of the hydrological balance in the Nile Delta.     

Holistic approach of GIS based Multi-Criteria Decision Analysis (MCDA) and WetSpass models to evaluate groundwater potential in Gelana watershed of Ethiopia

Appropriate quantification and identification of the groundwater distribution in a hydrological basin may provide necessary information for effective management, planning and development of groundwater resources. Groundwater potential assessment and delineation in a highly heterogeneous environment with limited Spatiotemporal data derived from Gelana watershed of Abaya Chamo lake basin is performed, using integrated multi-criteria decision analysis (MCDA), water and energy transfer between soil and plant and atmosphere under quasi-steady state (WetSpass) models. The outputs of the WetSpass model reveal a favorable structure of water balance in the basin studied, mainly using surface runoff. The simulated total flow and groundwater recharge are validated using river measurements and estimated baseflow at two gauging stations located in the study area, which yields a good agreement. The WetSpass model effectively integrates a water balance assessment in a geographical information system (GIS) environment. The WetSpass model is shown to be computationally reputable for such a remote complex setting as the African rift, with a correlation coefficient of 0.99 and 0.99 for total flow and baseflow at a significant level of p-value<0.05, respectively. The simulated annual water budget reveals that 77.22% of annual precipitation loses through evapotranspiration, of which 16.54% is lost via surface runoff while 6.24% is recharged to the groundwater. The calibrated groundwater recharge from the WetSpass model is then considered when determining the controlling factors of groundwater occurrence and formation, together with other multi-thematic layers such as lithology, geomorphology, lineament density and drainage density. The selected five thematic layers through MCDA are incorporated by employing the analytical hierarchy process (AHP) method to identify the relative dominance in groundwater potential zoning. The weighted factors in the AHP are procedurally aggregated, based on weighted linear combinations to provide the groundwater potential index. Based on the potential indexes, the area then is demarcated into low, moderate, and high groundwater potential zones (GWPZ). The identified GWPZs are finally examined using the existing groundwater inventory data (static water level and springs) in the region. About 70.7% of groundwater inventory points are coinciding with the delineated GWPZs. The weighting comparison shows that lithology, geomorphology, and groundwater recharge appear to be the dominant factors influence on the resources potential. The assessment of groundwater potential index values identify 45.88% as high, 39.38% moderate, and 14.73% as low groundwater potential zones. WetSpass model analysis is more preferable in the area like Gelana watershed when the topography is rugged, inaccessible and having limited gauging stations.     

Effects of climate change on the groundwater system in the Grote-Nete catchment,Belgium

The effects of climate change on the groundwater systems in the Grote-Nete catchment, Belgium, covering an area of 525 km², is modeled using wet (greenhouse), cold or NATCC (North Atlantic Thermohaline Circulation Change) and dry climate scenarios. Low, central and high estimates of temperature changes are adopted for wet scenarios. Seasonal and annual water balance components including groundwater recharge are simulated using the WetSpass model, while mean annual groundwater elevations and discharge are simulated with a steady-state MODFLOW groundwater model. WetSpass results for the wet scenarios show that wet winters and drier summers are expected relative to the present situation. MODFLOW results for wet high scenario show groundwater levels increase by as much as 79 cm, which could affect the distribution and species richness of meadows. Results obtained for cold scenarios depict drier winters and wetter summers relative to the present. The dry scenarios predict dry conditions for the whole year. There is no recharge during the summer, which is mainly attributed to high evapotranspiration rates by forests and low precipitation. Average annual groundwater levels drop by 0.5 m, with maximum of 3.1 m on the eastern part of the Campine Plateau. This could endanger aquatic ecosystem, shrubs, and crop production.     

城市扩张下的北京平原区降雨入渗补给量变化

结合WetSpass模型与地理信息系统(geographic information system,GIS)、遥感(remote sense,RS)技术分析了城市扩张引起的土地利用类型变化对北京平原区降水入渗补给量的影响.在估算出1982年和2007年降水入渗补给量的基础上,将2007年土地利用类型还原成1982年的情景重新估算,利用转移矩阵分析两年土地利用类型的相互转化关系,同时,基于GIS空间数据统计功能,计算出不同土地利用类型下的地下水补给量.结果表明,1982年至2007年,研究区内水浇地减少874 km2,其中517 km2转变为城镇建设用地.相对于1982年,2007年城镇建设用地扩张了831 km2,区内降水入渗补给量减少约3 000万m3.研究成果可以为北京平原区的地下水资源保护及土地资源配置提供较为科学的参考.      

Impact of recharge variations on water quality as indicated by excess air in groundwater of the Kalahari, Botswana

Groundwater is an important and often exclusive water resource in arid and semi-arid regions. The aim of the present paper was to gain insight into the processes and conditions that control the deterioration of groundwater quality in the semi-arid Kalahari of Botswana. Measurements of ³He, ⁴He, ²⁰Ne, ²²Ne, and of ¹⁴C of dissolved inorganic carbon (DIC) were combined with existing isotopic and hydrochemical data to investigate groundwater from the Ntane Sandstone Aquifer, which is affected by high nitrate concentrations of non-anthropogenic origin. All groundwater samples revealed neon concentrations in excess to air-saturated water, which we attributed to the addition of excess air during recharge. Neon concentrations ranged from values close to air saturation for ¹⁴C DIC rich samples (up to 80.5%MC) up to values of 90% in excess to air-saturated water for lower ¹⁴C DIC contents (2.6-61.3%MC). A strong linear correlation of excess Ne with nitrate concentrations suggests an intimate connection between groundwater quality and the processes and conditions during groundwater recharge. Low groundwater recharge rates under present-day semi-arid conditions are associated with low amounts of excess Ne and elevated nitrate concentrations. In contrast to this, higher excess Ne values in groundwater of lower ¹⁴C DIC and nitrate contents indicate that the high quality groundwater end-member presumably is related to higher groundwater table fluctuations during wetter climatic conditions in the past. We attribute the decline in groundwater quality with respect to nitrate to a decreasing rate and temporal variability of groundwater recharge, and to concurrent changes in biogeochemical activities following a transition to a drier climate during the Holocene. Under such conditions, a much stronger decrease in groundwater recharge compared to the release of nitrate from soil organic matter may result in elevated nitrate concentrations in the vadose zone and groundwater. This implies a strong impact of climate change on the transport of solutes like nitrate through the vadose zone which needs to be considered in predictions of future groundwater quality.     

Sustainable groundwater development in the coastal Tra Vinh province in Vietnam under saltwater intrusion and climate change

Three-dimensional transient groundwater flow and saltwater transport models were constructed to assess the impacts of groundwater abstraction and climate change on the coastal aquifer of Tra Vinh province (Vietnam). The groundwater flow model was calibrated with groundwater levels (2007–2016) measured in 13 observation wells. The saltwater transport model was compared with the spatial distribution of total dissolved solids. Model performance was evaluated by comparing observed and simulated groundwater levels. The projected rainfalls from two climate models (MIROC5 and CRISO Mk3.6) were subsequently used to simulate possible effects of climate changes. The simulation revealed that groundwater is currently depleted due to overabstraction. Towards the future, groundwater storage will continue to be depleted with the current abstraction regime, further worsening in the north due to saltwater intrusion from inland trapped saltwater and on the coast due to seawater intrusion. Notwithstanding, the impact from climate change may be limited, with the computed groundwater recharge from the two climate models revealing no significant change from 2017 to 2066. Three feasible mitigation scenarios were analyzed: (1) reduced groundwater abstraction by 25, 35 and 50%, (2) increased groundwater recharge by 1.5 and 2 times in the sand dunes through managed aquifer recharge (reduced abstraction will stop groundwater-level decline, while increased recharge will restore depleted storage), and (3) combining 50% abstraction reduction and 1.5 times recharge increase in sand dune areas. The results show that combined interventions of reducing abstraction and increasing recharge are necessary for sustainable groundwater resources development in Tra Vinh province.

     

Sensitivity Analysis of a Combined Groundwater Flow and Solute Transport Model Using Local-Grid Refinement: A Case Study

Combining groundwater flow models with solute transport models represents a common challenge in groundwater resources assessments and contaminant transport modeling. Groundwater flow models are usually constructed at somewhat larger scales (involving a coarser discretization) to include natural boundary conditions. They are commonly calibrated using observed groundwater levels and flows (if available). The groundwater solute transport models may be constructed at a smaller scale with finer discretization than the flow models in order to accurately delineate the solute source and the modeled target, to capture any heterogeneity that may affect contaminant migration, and to minimize numerical dispersion while still maintaining a reasonable computing time. The solution that is explored here is based on defining a finer grid subdomain within a larger coarser domain. The local-grid refinement (LGR) implemented in the Modular 3D finite-difference ground-water flow model (MODFLOW) code has such a provision to simulate groundwater flow in two nested grids: a higher-resolution sub-grid within a coarse grid. Under the premise that the interface between both models was well defined, a comprehensive sensitivity and uncertainty analysis was performed whereby the effect of a parameter perturbation in a coarser-grid model on transport predictions using a higher-resolution grid was quantified. This approach was tested for a groundwater flow and solute transport analysis in support of a safety evaluation of the future Belgian near-surface radioactive waste disposal facility. Our reference coarse-grid groundwater flow model was coupled with a smaller fine sub-grid model in two different ways. While the reference flow model was calibrated using observed groundwater levels at a scale commensurate with that of the coarse-grid model, the fine sub-grid model was used to run a solute transport simulation quantifying concentrations in a hypothetical well nearby the disposal facility. When LGR coupling was compared to a one-way coupling, LGR was found to provide a smoother flow solution resulting in a more CPU-efficient transport solution. Parameter sensitivities performed with the groundwater flow model resulted in sensitivities at the head observation locations. These sensitivities identified the recharge as the most sensitive parameter, with the hydraulic conductivity of the upper aquifer as the second most sensitive parameter in regard to calculated groundwater heads. Based on one-percent sensitivity maps, the spatial distribution of the observations with the highest sensitivities is slightly different for the upper aquifer hydraulic conductivity than for recharge. Sensitivity analyses were further performed to assess the prediction scaled sensitivities for hypothetical contaminant concentrations using the combined groundwater flow and solute transport models. Including all pertinent parameters into the sensitivity analysis identified the hydraulic conductivity of the upper aquifer as the most sensitive parameter with regard to the prediction of contaminant concentrations.     

Spatial and temporal constraints on regional-scale groundwater flow in the Pampa del Tamarugal Basin,Atacama Desert,Chile

Aquifers within the Pampa del Tamarugal Basin (Atacama Desert, northern Chile) are the sole source of water for the coastal city of Iquique and the economically important mining industry. Despite this, the regional groundwater system remains poorly understood. Although it is widely accepted that aquifer recharge originates as precipitation in the Altiplano and Andean Cordillera to the east, there remains debate on whether recharge is driven primarily by near-surface groundwater flow in response to periodic flood events or by basal groundwater flux through deep-seated basin fractures. In addressing this debate, the present study quantifies spatial and temporal variability in regional-scale groundwater flow paths at 20.5°S latitude by combining a two-dimensional model of groundwater and heat flow with field observations and δ¹⁸O isotope values in surface water and groundwater. Results suggest that both previously proposed aquifer recharge mechanisms are likely influencing aquifers within the Pampa del Tamarugal Basin; however, each mechanism is operating on different spatial and temporal scales. Storm-driven flood events in the Altiplano readily transmit groundwater to the eastern Pampa del Tamarugal Basin through near-surface groundwater flow on short time scales, e.g., 10⁰–10¹ years, but these effects are likely isolated to aquifers in the eastern third of the basin. In addition, this study illustrates a physical mechanism for groundwater originating in the eastern highlands to recharge aquifers and salars in the western Pampa del Tamarugal Basin over timescales of 10⁴–10⁵ years.     

Evaluation and numerical modeling of seawater intrusion in the Gaza aquifer (Palestine)

A numerical assessment of seawater intrusion in Gaza, Palestine, has been achieved applying a 3-D variable density groundwater flow model. A two-stage finite difference simulation algorithm was used in steady state and transient models. SEAWAT computer code was used for simulating the spatial and temporal evolution of hydraulic heads and solute concentrations of groundwater. A regular finite difference grid with a 400 m² cell in the horizontal plane, in addition to a 12-layer model were chosen. The model has been calibrated under steady state and transient conditions. Simulation results indicate that the proposed schemes successfully simulate the intrusion mechanism. Two pumpage schemes were designed to use the calibrated model for prediction of future changes in water levels and solute concentrations in the groundwater for a planning period of 17 years. The results show that seawater intrusion would worsen in the aquifer if the current rates of groundwater pumpage continue. The alternative, to eliminate pumpage in the intruded area, to moderate pumpage rates from water supply wells far from the seashore and to increase the aquifer replenishment by encouraging the implementation of suitable solutions like artificial recharge, may limit significantly seawater intrusion and reduce the current rate of decline of the water levels.     

Assessment of factors influencing groundwater-level change using groundwater flow simulation,considering vertical infiltration from rice-planted and crop-rotated paddy fields in Japan

Assessing factors that influence groundwater levels such as land use and pumping strategy, is essential to adequately manage groundwater resources. A transient numerical model for groundwater flow with infiltration was developed for the Tedori River alluvial fan (140 km²), Japan. The main water input into the groundwater body in this area is irrigation water, which is significantly influenced by land use, namely paddy and upland fields. The proposed model consists of two models, a one-dimensional (1-D) unsaturated-zone water flow model (HYDRUS-1D) for estimating groundwater recharge and a 3-D groundwater flow model (MODFLOW). Numerical simulation of groundwater flow from October 1975 to November 2009 was performed to validate the model. Simulation revealed seasonal groundwater level fluctuations, affected by paddy irrigation management. However, computational accuracy was limited by the spatiotemporal data resolution of the groundwater use. Both annual groundwater levels and recharge during the irrigation periods from 1975 to 2009 showed long-term decreasing trends. With the decline in rice-planted paddy field area, groundwater recharge cumulatively decreased to 61 % of the peak in 1977. A paddy-upland crop-rotation system could decrease groundwater recharge to 73–98 % relative to no crop rotation.     

Recharge source and hydrogeochemical evolution of shallow groundwater in a complex alluvial fan system,southwest of North China Plain

Many cities around the world are developed at alluvial fans. With economic and industrial development and increase in population, quality and quantity of groundwater are often damaged by over-exploitation in these areas. In order to realistically assess these groundwater resources and their sustainability, it is vital to understand the recharge sources and hydrogeochemical evolution of groundwater in alluvial fans. In March 2006, groundwater and surface water were sampled for major element analysis and stable isotope (oxygen-18 and deuterium) compositions in Xinxiang, which is located at a complex alluvial fan system composed of a mountainous area, Taihang Mt. alluvial fan and Yellow River alluvial fan. In the Taihang mountainous area, the groundwater was recharged by precipitation and was characterized by Ca–HCO₃ type water with depleted δ¹⁸O and δD (mean value of −8.8‰ δ¹⁸O). Along the flow path from the mountainous area to Taihang Mt. alluvial fan, the groundwater became geochemically complex (Ca–Na–Mg–HCO₃–Cl–SO₄ type), and heavier δ¹⁸O and δD were observed (around −8‰ δ¹⁸O). Before the surface water with mean δ¹⁸O of −8.7‰ recharged to groundwater, it underwent isotopic enrichment in Taihang Mt. alluvial fan. Chemical mixture and ion exchange are expected to be responsible for the chemical evolution of groundwater in Yellow River alluvial fan. Transferred water from the Yellow River is the main source of the groundwater in the Yellow River alluvial fan in the south of the study area, and stable isotopic compositions of the groundwater (mean value of −8.8‰ δ¹⁸O) were similar to those of transferred water (−8.9‰), increasing from the southern boundary of the study area to the distal end of the fan. The groundwater underwent chemical evolution from Ca–HCO₃, Na–HCO₃, to Na–SO₄. A conceptual model, integrating stiff diagrams, is used to describe the spatial variation of recharge sources, chemical evolution, and groundwater flow paths in the complex alluvial fan aquifer system.     

Determining groundwater degradation from irrigation in desert-marginal northern China

Groundwater degradation from irrigated agriculture is of concern in semi-arid northern China. Data-scarcity often means the causes and extent of problems are not fully understood. An irrigated area in Inner Mongolia was studied, where abstraction from an unconfined Quaternary aquifer has increased threefold over 20 years to 20 million m³/year; groundwater levels are falling at up to 0.5 m/year; and groundwater is increasingly mineralised (TDS increase from 400 to 700–1,900 mg/L), with nitrate concentrations up to 137 mg/L N. Residence-time (chlorofluorocarbons), stable-isotope and hydrogeochemical indicators helped develop a conceptual model of groundwater system evolution, demonstrating a direct relationship between modern water proportion and the degree of groundwater mineralisation, indicating that irrigation-water recycling is reducing groundwater quality. The investigations suggest that before irrigation development, active recharge to the aquifer from wadis significantly exceeded groundwater inflow from nearby mountains, previously held to be the main groundwater input. Away from active wadis, groundwater is older with a probable pre-Holocene component. Proof-of-concept groundwater modelling supports geochemical evidence, indicating the importance of wadi recharge and irrigation return flows. Engineering works protecting the irrigated area from flooding have reduced good quality recharge; active recharge is now dominated by irrigation returns, which are degrading the aquifer.     

Delineating volcanic aquifer recharge areas using geochemical and isotopic tools

Relative recharge areas are evaluated using geochemical and isotopic tools, and inverse modeling. Geochemistry and water quality in springs discharging from a volcanic aquifer system in Guatemala are related to relative recharge area elevations and land use. Plagioclase feldspar and olivine react with volcanically derived CO₂ to produce Ca-montmorillonite, chalcedony and goethite in the groundwater. Alkalinity, Mg, Ca, Na, and SiO_2(aq) are produced, along with minor increases in Cl and SO₄ concentrations. Variations in groundwater δD and δ¹⁸O values are attributed to recharge elevation and used in concert with geochemical evolution to distinguish local, intermediate, and regional flow systems. Springs with geochemically inferred short flow paths provided useful proxies to estimate an isotopic gradient for precipitation (??.67 δ¹⁸O/100?m). No correlation between spring discharge and relative flow-path length or interpreted recharge elevation was observed. The conceptual model was consistent with evidence of anthropogenic impacts (sewage and manure) in springs recharged in the lower watershed where livestock and humans reside. Spring sampling is a low-budget approach that can be used to develop a useful conceptual model of the relative scale of groundwater flow (and appropriate watershed protection areas), particularly in volcanic terrain where wells and boreholes are scarce.     

Environmental impact assessment of risk associated with groundwater overdraft remediation in cone of depression,Jining, China

In this paper, a simple scenario and probabilistic approach is used to assess the potential groundwater risk due to proposed overdraft remedial actions in cone of depression, Jining City, China. Focusing on the concentrations of Chloride ions (Cl⁻) and total hardness (TH), the impact of artificial recharge and reduced pumping on groundwater quality and quantity is analysed by using the three-dimensional finite difference groundwater flow and transport model, Visual MODFLOW, to simulate groundwater flow and transport within the study area based on scenarios, and utilizes SURFER software to map risk levels. Although 5, 10 or 15% reduced pumping with artificial recharge leads to more decrease in Cl⁻ and TH concentrations than the 25%, less volume increase is achieved for the remediation of land subsidence and other environmental problems in the cone of depression. The Cl⁻ concentrations in recovered groundwater are within the desired concentration of 200 mg/l; however, TH in some cases are above the maximum permissible limit of 500 mg/l, with an exceedence probability of about 0.67 for recharge and recharge with reduced pumping at 25%. The presence of fractures and hydrogeological complexity greatly determines impacts of remediation, and the 22% reduced pumping with artificial recharge offers an optimum strategy for overdraft remediation in the Jining cone of depression.     

Parallel simulation of groundwater non-point source pollution using algebraic multigrid preconditioners

The simulation of non-point source pollution in agricultural basins is a computationally demanding process due to the large number of individual sources and potential pollution receptors (e.g., drinking water wells). In this study, we present an efficient computational framework for parallel simulation of diffuse pollution in such groundwater basins. To derive a highly detailed velocity field, we employed algebraic multigrid (AMG) preconditioners to solve the groundwater flow equation. We compare two variants of AMG implementations, the multilevel preconditioning provided by Trilinos and the BoomerAMG provided by HYPRE. We also perform a sensitivity analysis on the configuration of AMG methods to evaluate the application of these libraries to groundwater flow problems. For the transport simulation of diffuse contamination, we use the streamline approach, which decomposes the 3D transport problem into a large number of 1D problems that can be executed in parallel. The proposed framework is applied to a 2,600-km² groundwater basin in California discretized into a grid with over 11 million degrees of freedom. Using a Monte Carlo approach with 200 nitrate loading realizations at the aquifer surface, we perform a stochastic analysis to quantify nitrate breakthrough prediction uncertainty at over 1,500 wells due to random, temporally distributed nitrate loading. The results show that there is a significant time lag between loading and aquifer response at production wells. Generally, typical production wells respond after 5–50 years depending on well depth and screen length, while the prediction uncertainty for nitrate in individual wells is very large—approximately twice the drinking water limit for nitrate.     

Groundwater resource sustainability in the Wadi Watir delta, Gulf of Aqaba, Sinai, Egypt

The Wadi Watir delta, in the arid Sinai Peninsula, Egypt, contains an alluvial aquifer underlain by impermeable Precambrian basement rock. The scarcity of rainfall during the last decade, combined with high pumping rates, resulted in degradation of water quality in the main supply wells along the mountain front, which has resulted in reduced groundwater pumping. Additionally, seawater intrusion along the coast has increased salinity in some wells. A three-dimensional (3D) groundwater flow model (MODFLOW) was calibrated using groundwater-level changes and pumping rates from 1982 to 2009; the groundwater recharge rate was estimated to be 1.58?×?10⁶ m³/year. A variable-density flow model (SEAWAT) was used to evaluate seawater intrusion for different pumping rates and well-field locations. Water chemistry and stable isotope data were used to calculate seawater mixing with groundwater along the coast. Geochemical modeling (NETPATH) determined the sources and mixing of different groundwaters from the mountainous recharge areas and within the delta aquifers; results showed that the groundwater salinity is controlled by dissolution of minerals and salts in the aquifers along flow paths and mixing of chemically different waters, including upwelling of saline groundwater and seawater intrusion. Future groundwater pumping must be closely monitored to limit these effects.     

Assessment of paleo-recharge under the Fennoscandian Ice Sheet and its impact on regional groundwater flow in the northern Baltic Artesian Basin using a numerical model

The study investigates the mechanism of glacial meltwater recharge under the Fennosciandian Ice Sheet during the last glacial maximum (LGM) and its impact on regional groundwater flow in the northern Baltic Artesian Basin (BAB) in Estonia and Latvia. The current hypothesis is that a flow reversal occurred in the BAB due to subglacial recharge during the LGM. This hypothesis is supported by an extensive dataset of geochemical and isotopic measurements in the groundwater of northern Estonia, exhibiting significant depletion in δ¹⁸O with respect to modern precipitation. To verify the consistency of this hypothesis and better understand groundwater flow dynamics during the LGM period, a numerical model is developed for this area. Two cross-sectional models have been created across the northern BAB, in which groundwater flow and the transport of δ¹⁸O have been simulated from the beginning of the LGM to present-day. Several simulations were performed with different subglacial boundary conditions, to investigate the uncertainty related to subglacial recharge of meltwater during the LGM and the subsequent flow reversal in the northern BAB. Several simulations provide a satisfying fit between computed and observed values of δ¹⁸O, which means that the hypothesis of subglacial recharge of meltwater is consistent with δ¹⁸O distribution. The numerical model suggests that preservation of meltwater in northern Estonia is controlled by confining layers and the proximity to the outcrop area of aquifers, located in the Gulf of Finland. The results also suggest that glacial meltwater has been preserved under the Baltic Sea in the Gulf of Riga.     

Effects of rainwater-harvesting-induced artificial recharge on the groundwater of wells in Rajasthan, India

In light of the increasing deterioration of groundwater supplies in Rajasthan, India, rainwater harvesting practices in southern Rajasthan were studied to determine the effects of artificially recharged groundwater on the supply and quality of local groundwater. A physical and geochemical investigation utilizing environmental tracers (δ¹⁸O and Cl^–), groundwater level and groundwater quality measurements, and geological surveys was conducted with two objectives: (1) to quantify the proportion of artificially recharged groundwater in wells located near rainwater harvesting structures and (2) to examine potential effects of artificial recharge on the quality of groundwater in these wells. A geochemical mixing model revealed that the proportion of artificial recharge in these wells ranged from 0 to 75%. Groundwater tracer, water table, and geological data provided evidence of complex groundwater flow and were used to explain the spatial distribution of artificial recharge. Furthermore, wells receiving artificial recharge had improved groundwater quality. Statistical analysis revealed a significant difference between the water quality in these wells and wells determined not to receive artificial recharge, for electrical conductivity and SO ₄ ^– . The findings from this study provide quantitative evidence that rainwater harvesting structures in southern Rajasthan influence the groundwater supply and quality of nearby wells by artificially recharging local groundwater.     

An evaluation of the relative importance of the effects of climate change and groundwater extraction on seawater intrusion in coastal aquifers in Atlantic Canada

To investigate the relative importance of projected sea-level rise, climate change effects on recharge, and groundwater extraction on seawater intrusion in important coastal aquifers in Atlantic Canada, a three-dimensional numerical model of density-dependent groundwater flow coupled with solute transport was developed for the Richibucto region of New Brunswick. The model was used, with an efficient 2^k factorial design approach, to perform simulations for the period 2011–2100. The results of the factorial analyses indicate that the relative importance of the three factors investigated varies depending on the model location considered. The effect of declining recharge is most significant at shallow to intermediate depths along the freshwater–seawater transition zone, while the effect of increasing pumping rates dominates at a location relatively close to the well field. The effect of sea-level rise is shown to be significant only at the much deeper inland toe of the transition zone. The spatial variation in importance is related to how different model boundary conditions influence freshwater flow at the different locations within the model domain. This investigation indicates that sea-level rise has the least significant effect (of the three factors considered) on future seawater intrusion in sandstone aquifers in the Richibucto region.