Groundwater Withdrawal Research Articles

Linked Data-Driven, Physics-Based Modeling of Pumping-Induced Subsidence with Application to Bangkok, Thailand.

Research into land subsidence caused by groundwater withdrawal is hindered by the availability of measured heads, subsidence, and forcings. In this paper, a parsimonious, linked data-driven and physics-based approach is introduced to simulate pumping-induced subsidence; the approach is intended to be applied at observation well nests. Time series analysis using response functions is applied to simulate heads in aquifers. The heads in the clay layers are simulated with a one-dimensional diffusion model, using the heads in the aquifers as boundary conditions. Finally, simulated heads in the layers are used to model land subsidence. The developed approach is applied to the city of Bangkok, Thailand, where relatively short time series of head and subsidence measurements are available at or near 23 well nests; an estimate of basin-wide pumping is available for a longer period. Despite the data scarcity, data-driven time series models at observation wells successfully simulate groundwater dynamics in aquifers with an average root mean square error (RMSE) of 2.8 m, relative to an average total range of 21 m. Simulated subsidence matches sparse (and sometimes very noisy) land subsidence measurements reasonably well with an average RMSE of 1.6 cm/year, relative to an average total range of 5.4 cm/year. Performance is not good at eight out of 23 locations, most likely because basin-wide pumping is not representative of localized pumping. Overall, this study demonstrates the potential of a parsimonious, linked data-driven, and physics-based approach to model pumping-induced subsidence in areas with limited data.

Long-term investigation of the irrigation intervals and supplementary irrigation strategies effects on winter wheat in the U.S. Central High Plains based on a combination of crop modeling and field studies

Implementing optimized irrigation strategies to achieve acceptable wheat yield while conserving groundwater resources is critical for extensively irrigated regions. Field experiments regarding winter wheat irrigation management were conducted for two growing seasons (2013–2014 and 2014–2015) to calibrate and validate the AquaCrop model. To determine proper parameterization of the AquaCrop model for various data availability conditions, the model performance was tested for calibration based on three different datasets, including a) calibration based on a full irrigation condition (CA 1), b) calibration based on a dryland condition (CA 2), and c) calibration based on diverse irrigation conditions, which included full irrigation, different deficit irrigation levels, and dryland (CA 3). The CA 3 calibration scenario resulted in the model's highest accuracy in simulating biomass, grain yield, total soil water, and seasonal evapotranspiration during the calibration and validation process. However, the model performance was also convenient when calibration based on full irrigation conditions was pursued. The calibrated and validated AquaCrop model based on the CA3 was used to analyze long-term (36-year) irrigation management scenarios for identifying the optimized water-conservative winter wheat irrigation strategies. The long-term analysis emphasized that increasing irrigation intervals from 7 to 15 or 20 days could reduce groundwater withdrawal by 52–64 % while expecting a 14–19 % reduction in grain yield. By implementing one 25 mm irrigation after planting and two 25 mm depth irrigation events at jointing or leaf growth stages, 91 % of Kansas's average winter wheat yield (2.88 tons/ha) could be achieved. During wet years, the 1.2 ton/ha reduction in biomass and 0.27–0.62 ton/ha reduction in grain yield is expected irrespective of irrigation management strategies. Considerable uncertainty was detected in grain yield production during dry years under dryland (rainfed) conditions. The maximum winter wheat production could be achieved under normal years, establishing a balance between precipitation and heat units. Our agro-hydrological analysis revealed that the highest irrigation water use efficiency could be attained by the application of two 25 mm irrigation events at jointing or flag leaf growth stages during normal years, the application of two irrigation events during jointing or heading during wet years, and two 25 mm irrigations at heading or flowering during dry years. The results of this research could be used as a baseline for producers of the U.S. Central High Plains and semi-arid regions with similar climate characteristics to cope with water scarcity in winter wheat production.

Visual MODFLOW, solute transport modeling, and remote sensing techniques for adapting aquifer potentiality under reclamation and climate change impacts in coastal aquifer

Global environmental changes, such as climate change and reclamation alterations, significantly influence hydrological processes, leading to hydrologic nonstationarity and challenges in managing water availability and distribution. This study introduces a conceptual underpinning for the rational development and sustainability of groundwater resources. As one of the areas intended for the development projects within the Egyptian national plan for the reclamation of one and a half million acres; hundreds of pumping wells were constructed in the Moghra area to fulfill the reclamation demand. This study investigates the long-term impacts of exploiting the drilled pumping wells under climate change. The approach is to monitor the groundwater levels and the salinity values in the Moghra aquifer with various operational strategies and present proposed sustainable development scenarios. The impact of global warming and climate change is estimated for a prediction period of 30 years by using satellite data, time series geographical analysis, and statistical modeling. Using MODFLOW and Solute Transport (MT3DMS) modules of Visual MODFLOW USGS 2005 software, a three-dimensional (3D) finite-difference model is created to simulate groundwater flow and salinity distribution in the Moghra aquifer with the input of forecast downscaling (2020–2050) of main climatic parameters (PPT, ET, and Temp). The optimal adaptation-integrated scenario to cope with long-term groundwater withdrawal and climate change impacts is achieved when the Ministry of Irrigation and Water Resources (MWRI) recommends that the maximum drawdown shouldn’t be more significant than 1.0 m/ year. In this scenario, 1,500 pumping wells are distributed with an equal space of 500 m, a pumping rate of 1,200 m3/day and input the forecast of the most significant climatic parameters after 30 years. The output results of this scenario revealed a drawdown level of 42 m and a groundwater salinity value of 16,000 mg/l. Climate change has an evident impact on groundwater quantity and quality, particularly in the unconfined coastal aquifer, which is vulnerable to saltwater intrusion and pollution of drinking water resources. The relationship between climate change and the hydrologic cycle is crucial for predicting future water availability and addressing water-related issues.

A method for estimating maximum safe installable power for groundwater extraction with application to Africa

New groundwater development is a likely way to meet growing global water demand but needs careful management. To help inform the sustainable development of groundwater resources, a novel method based on the maximum safe installable power for water pumping systems and the maximum safe remaining installable power (considering current abstraction) is developed. The proposed model couples energy, technology and hydrogeological parameters, and is then developed to compute the maximum power that can be safely installed per km2 without exceeding a maximum annual pumpable volume, calculated through available recharge and storage. The model is applied to estimate the maximum safe installable power across Africa with a 0.2-degree resolution, using available energy and hydrogeological data. Constrained by recharge (considering that 25 % of the annual recharge is available for utilization), the maximum safe installable power ranges between 0 and 9960 W/km2 across Africa with regions such as the Congo Basin (~340 W/km2), and western Africa between the Ivory Coast and Nigeria (~230 W/km2) identified as having high potential for sustainable pumping system development. Constrained by storage (considering that 0.1 % of storage can be withdrawn per year), it ranges between 0 and 13,425 W/km2 highlighting the ability to harness storage for pumping system development in the large aquifers of Northern Sahara (~8720 W/km2). When considering limitations posed by both groundwater recharge and aquifer storage (considering that 25 % of the annual recharge and 0.1 % of storage can be withdrawn per year), along with current groundwater withdrawal, 93 % of the maximum safe installable power remains still installable on average across Africa. Nevertheless, 2 % of the locations are estimated to be already experiencing overexploitation, particularly in Sudan, northern Africa, and northeastern South Africa. These findings provide a novel and adaptable way to examine water security, which can assist institutions in targeting investments to meet water demand sustainably.

Groundwater in the city of Pesaro (Marche, Italy): anthropic impact and interference with the urban environment

Groundwater in the urban area of Pesaro is a significant local resource, supplementing the main surfacewater supply from the external Metauro watershed. Nevertheless, since Pesaro’s groundwater is tapped from a coastal alluvial aquifer in a heavily anthropized environment, its use involves mutual impacts with the urban environment, which poses several constraints. Our work schematized the main interferences of groundwater with building works such as housing (i.e. foundations and underground works), hydraulic infrastructures, wells’ network development, and geoexchange projects. Moreover, the presence of a seawater wedge, historically intruding the aquifer, must be taken into account, alongside its potential evolution in changing climate scenarios. Contaminated sites were also analyzed, as they represent a critical issue when further groundwater withdrawals are planned. Due to these complex interactions, there is a need for renewed commitment to studying and monitoring the local water resource. This is also essential for promoting awareness among both citizens and the local government that its use will remain crucial for a long time.

Modern Drip Irrigation Technology on Tomato and Head Cabbage Production and Its Economic Feasibility at Misrak Silti District, Siltie Zone, Ethiopia

In Ethiopia various techniques were applied to improve on-farm irrigation water management under surface irrigation, especially on furrow irrigation system for last many years, however, it was very difficult to achieve threshold limit of water use efficiency, wisely use of scarce water resources in irrigated field, crop water productivity, precisely controlled application of irrigation water to plant roots, uniform water delivery to all plants, crop yield and its quality, regulate flow, deliver optimum crop water requirement, field water losses, groundwater withdrawal and save labor. It is great practical significance and series of measures using drip irrigation system to solve mentioned problems in irrigated field, because the irrigated agriculture is largest water-consuming sector in this area. This study was aimed to demonstrate drip irrigation technology, prove its economic feasibility and create skills of farmers and extension experts on implementation of the drip system on tomato (Galila 555) and head cabbage (Copen Hagen) productions within groundwater source area on Balo koriso main station in Misrak Silti woreda, Siltie zone, southern Ethiopia. The activity was done for four consecutive years (2019 to 2022 G. C) on fixed plot of 2500 square meter area in collaboration with Agricultural Research Institute and Techno serve project. During demonstration all the costs were considered to economic feasibility including initial investment cost. The study found that the use of drip irrigation saved 3690 m<sup>3</sup>/ha of water compared with use of furrow irrigation. This saved water may irrigate additional 2ha area of land by drip irrigation. This study also revealed that, a net income of 313066.04 ETB/ha in single irrigation season was obtained using modern drip irrigation technology. The farmers and local experts recognized noticeable saving water, fuel cost, irrigating time and labour, good crop performance from demonstrated drip irrigation system and shown interest to use the technology abundantly. Therefore, scaling up locally manageable drip irrigation system around the study area through comprehensive training of farmers and supportive staffs and allowing sufficient local market for drip kits and incentives or loan will be economically feasible and affordable to increase income and saves irrigation water.

GROUNDWATER CONDITION AND VULNERABILITY IN THE ATHABASCA AND COLD LAKE OIL SANDS REGIONS: GAPS, OPPORTUNTIES AND CHALLENGES

Oil sands development in the Athabasca and Cold Lake oil sands regions of Alberta has raised concerns about potential impacts to groundwater and groundwater dependent ecosystems. This review summarizes the current state of understanding as to how oil sands mining and in situ activities can affect groundwater systems using a stressor-mechanism-response framework. Specific oil sands activities and practices are reviewed, and where possible, described in terms of how they can impact hydraulic head, the hydraulic properties of aquifers, recharge and transport of constituents of concern and linked to observed or potential impacts to groundwater quantity and quality. Groundwater is an important component of the water balance in boreal ecosystems, and specific vulnerabilities related to development are reviewed, including water use, landscape disturbance, groundwater withdrawals, tailings pond seepage, deep well disposal and thermal impacts. Knowledge gaps include lack of baseline data and monitoring of the quantity and quality of groundwater discharge to rivers, lakes and wetlands. One key monitoring challenge is attribution of hydrogeologic responses to specific oil sands stressors given the range of other natural and anthropogenic factors contributing to their variability. Quantitative groundwater exchange mapping, regional-scale isotope mass balance assessment, and broader incorporation of isotopic and geochemical tracers for fingerprinting water sources and incorporation of Indigenous Knowledge appear promising for improved effectiveness of monitoring.

Are regional groundwater models suitable for simulating wetlands, rivers and intermittence? The example of the French AquiFR platform

Predicting and managing water resources at regional scale under different climate and socio-economic scenarios is crucial to support drinking water supply and other sectors. At the same time, protecting rivers and wetlands from pollutions and droughts is essential and must include groundwater given its contribution to surface water. Yet, assessing temporal and spatial variability of groundwater contributions to surface water is constrained due to limited observations. This study aims to quantify the spatio-temporal distribution of groundwater discharge (hereafter called “GW discharge”) zones, i.e. the groundwater flow to rivers and wetlands, estimated by calibrated regional groundwater flow models (French AquiFR platform). We compare simulation results with two types of surface observations: (i) the spatial distribution of surface water (BD TOPO, RAMSAR and Natura 2000 database) and (ii) an innovative datasets, the river intermittence observed in headwaters (ONDE network). Results show that simulated GW discharge zones are consistent with the observed location of rivers and wetlands. Time variations in GW discharge are well correlated with the intermittence observed at 396 of the 515 selected stations. Of these, groundwater model continues to feed surface water upstream of the station for ∼75 % of observed river drying up events, which may be consistent with a small alluvial flow. The groundwater withdrawals are shown to have a strong impact on the GW discharge and thus on the river intermittence.

Response to Pumping of Wells in Carbonate and Karst Aquifers and Effect on the Assessment of Sustainable Well Yield: Some Examples from Southern Italy

Carbonate and karst aquifers are of great importance for human water supplies, for supporting aquatic habitats and providing ecosystem services. Optimizing the groundwater withdrawals is therefore essential for obtaining the maximum flow rate for human purposes while minimizing the negative effects on the environment. In particular, when the abstraction of groundwater occurs through wells, the problem of defining the sustainable yield arises. This study analyzes pumping tests conducted in carbonate and karst aquifers in southern Italy to derive indications for defining the sustainable yield of yields. The four examined cases concern the Mesozoic–Cenozoic platform and transition pelagic carbonate rocks characterized by different degree of fracturing and karstification and hosting a carbonate aquifer with variable average groundwater yields. The analysis compared drawdown–time trends and their derivatives for 35 pumping tests with theoretical curves to identify the flow dimension. Parameters useful for examining the well yields were then determined. The results show that the response to the pumping of the investigated aquifers is very variable, both among the different sites and within the same site. Well yields are very different due to aquifer heterogeneity, local hydrostratigraphy and structural setting, and position of the pumping center within the groundwater flow system. To determine the operational pumping rate for a well in this environment, this study emphasizes the importance of analyzing drawdown trends over time to correctly predict the well’s long-term response to pumping. Specifically, when pumping induces a steady-state drawdown response, the focus for defining the sustainable abstraction shifts to the basin or aquifer scale. Conversely, when a transient drawdown response to pumping results, the well’s capacity to capture groundwater becomes the primary factor for well yield and its sustainability.

Enhancing smallholder agricultural production through sustainable use of shallow groundwater in the Borkena catchment, Awash River Basin, Ethiopia

ABSTRACT With increasing shallow groundwater use for agricultural purposes, understanding the spatiotemporal variability in recharge rates, storage capacity, and its interaction with surface waters becomes crucial for its sustainable management. An integrated SWAT–MODFLOW model is developed to assess shallow groundwater availability in the Borkena catchment. The model is calibrated using streamflow and static groundwater level data. Results show that groundwater recharge in the catchment is 85 mm/a, representing 11% of the mean annual rainfall. Shallow groundwater resources exist across nearly 42% of the Borkena catchment. The percentage of shallow groundwater withdrawal to groundwater recharge is very low (0.1%), signifying the potential for increased shallow groundwater development. However, caution must be taken as its uncontrolled expansion may result in a high risk of depletion. This integrated modeling is one of the few efforts conducted to provide important information regarding shallow groundwater potential in the Borkena catchment, which is essential for the resilience of small-scale producers in the continued growing water demand and climate change.

Estimating irrigation water use from remotely sensed evapotranspiration data: Accuracy and uncertainties at field, water right, and regional scales

Irrigated agriculture is the dominant user of water globally, but most water withdrawals are not monitored or reported. As a result, it is largely unknown when, where, and how much water is used for irrigation. Here, we evaluated the ability of remotely sensed evapotranspiration (ET) data, integrated with other datasets, to calculate irrigation water withdrawals and applications in an intensively irrigated portion of the United States. We compared irrigation calculations based on an ensemble of satellite-driven ET models from OpenET with reported groundwater withdrawals from hundreds of farmer irrigation application records and a statewide flowmeter database at three spatial scales (field, water right group, and management area). At the field scale, we found that ET-based calculations of irrigation agreed best with reported irrigation when the OpenET ensemble mean was aggregated to the growing season timescale (bias = 1.6–4.9 %, R2 = 0.53–0.74), and agreement between calculated and reported irrigation was better for multi-year averages than for individual years. At the water right group scale, linking pumping wells to specific irrigated fields was the primary source of uncertainty. At the management area scale, calculated irrigation exhibited similar temporal patterns as flowmeter data but tended to be positively biased with more interannual variability. Disagreement between calculated and reported irrigation was strongly correlated with annual precipitation, and calculated and reported irrigation agreed more closely after statistically adjusting for annual precipitation. The selection of an ET model was also an important consideration, as variability across ET models was larger than the potential impacts of conservation measures employed in the region. From these results, we suggest key practices for working with ET-based irrigation data that include accurately accounting for changes in soil moisture, deep percolation, and runoff; careful verification of irrigated area and well-field linkages; and conducting application-specific evaluations of uncertainty.

Semi-analytical solutions for land subsidence due to groundwater withdrawal

This paper introduces novel semi-analytical models tailored for estimating land subsidence resulting from groundwater extraction in confined aquifers. These models offer high scalability, allowing them to be applied to various well configurations and pumping schedules. Their development involves the numerical integration of two key analytical solutions: the “nucleus of strain” (NoS) (Mindlin and Chen, 1950), which represents a localised zone within the aquifer where a unit change in pore pressure leads to deformation and subsequent surface displacement, and the classic Theis equation (Theis, 1935) for the pore pressure changes induced by a constant-rate well pumping from a laterally unbounded aquifer. These integrations yield surface displacement components, both horizontal and vertical, expressed as functions of two dimensionless spatial–temporal variables, which encompass aquifer depth, thickness, well placement, pumping schedules, and critical hydro-geomechanical parameters like hydraulic conductivity, porosity, vertical compressibility, and water compressibility. Proposed are two distinct modelling approaches: one employing a lookup table (LT) derived from numerical integration results, and the other providing direct closed-form surface displacement solutions by fitting LT data with “hinge models”, which use piecewise-linear functions linked by sigmoidal curves for computational efficiency. In both cases, surface displacement components are estimated by plugging in the dimensionless variables. Conditions of variable pumping from multiple wells can be addressed by applying superposition of solutions. In essence, these semi-analytical models offer swift computational capabilities for understanding and forecasting land subsidence dynamics. Their scalability makes them adaptable to a wide array of well configurations and scheduling scenarios, rendering them valuable for numerous applications. They are particularly significant for providing preliminary estimates of the impacts of groundwater development, conducting “what-if” tests, and performing sensitivity analyses to identify key factors affecting land subsidence risk. This underscores the importance of these models in sustainable groundwater resource management and in mitigating land subsidence and its associated consequences.

REUSE OF TREATED WATER FROM BOJONGSOANG DOMESTIC WASTEWATER TREATMENT INSTALLATION AS WATER FOR INDUSTRIAL PURPOSES

Subdistrict Bojongsoang and its surroundings are based on the 2016-2036 Bandung Regency Spatial Planning, the designated area where the industry is. Currently, there are 125 existing industries operating. The area is also included in groundwater withdrawal through well artesian on a massive scale, including into the national protected area (CAT: Groundwater Basin) Bandung - Soreang Area. Excessive groundwater extraction poses challenges coupled with significant water demand from society and industry. This will result in deteriorated groundwater conditions that require prompt restoration, having already suffered damage. On the other hand, there is Bojongsoang WWTP, which can treat wastewater at 2800 l/s; meanwhile, at the moment, this capacity processing is utilized about 37-40% or about 100m3/day. This research uses a descriptive method. The research results prove that reusing water from Bojongsoang WWTP could produce raw water for the necessity industry of 300 l/s.

The influence of withdrawal-recharging pattern on the deformation characteristics of sand in confined aquifer

Land subsidence caused by excessive groundwater withdrawal has become a global hazard, which demands further researches and the potential measures to control. Using the FlowTrac Ⅱ consolidation test system, six compression tests were designed to investigate the stress state and stress paths of sand within confined aquifers under conditions of withdrawal and recharging. The deformation characteristics of aquifer sand were studied under different withdrawal-recharging patterns. During pumping and recharge processes, sand deformation responses were observed to lag behind changes in applied stress. The characteristics of this hysteresis effect on deformation were summarized. The alternating phenomenon of rebound and compression of sand deformation under the recharging process is analyzed. When the recharging effect was relatively small than withdrawing effect under the stable withdrawal-recharging pattern, the compression deformation was observed in the recharging process. The research results provide a rational explanation for the continuous compression deformation of the aquifer during groundwater level recovery and offer experimental evidence for the rational design of artificial groundwater recharge in engineering construction.

Using neural network modeling to improve the detection accuracy of land subsidence due to groundwater withdrawal

Using neural network modeling to improve the detection accuracy of land subsidence due to groundwater withdrawal

Cover collapse sinkhole formation is delayed in time and uncorrelated to distance from quarry in a long-term study of a karst basin

Collapse sinkholes are problematic for the hazards they bring to property, infrastructure, and human life but the prediction of their formation in time and space remains elusive. We make use of long-term data gathered in a small, highly monitored, hydrologically constrained basin; Primrose Creek, Bucks County, Pennsylvania, USA, to examine the relation between groundwater extraction, precipitation, terrain position, and cover collapse sinkhole development. The selected site is unique as a natural laboratory to explore induced sinkhole processes in that it is a clearly demarcated topographic basin with geological structural and lithologic boundaries; the stressor is extreme, has been plainly identified, and effectively captures all upstream runoff and groundwater; the groundwater withdrawals and water-level elevations have been closely monitored over a period of 20 years; and the general record of land conditions goes back more than a century. We evaluated the history of a carbonate-rock quarry, its development and expansion over many decades, as well as associated precipitation and ground-water level data, and cover collapse sinkhole occurrence within the basin. Interception of discrete high conductivity zones in the Paleozoic carbonate rocks by quarry expansion was an important factor in propagating collapse occurrence. Collapses occurred mainly, though not exclusively, along drainageways and at lithologic boundaries. They formed in time-based clusters with delay of up to several years. Surprisingly, distance of collapse formation from the pumping area was not directly correlated with time elapsed; i.e., collapses did not begin to form close to the quarry, and then later form farther away. Monitoring wells were unsuccessful in delineating the pumping zone of influence although some documented water table decline and response to weather and anthropogenic events. This comprehensive study shows that in order to successfully monitor and control impacts from quarry dewatering in Paleozoic carbonate rocks it is critical to 1) develop a detailed conceptual model of groundwater pathways, 2) try to avoid excavating into high permeability zones, 3) clearly document karst features when exposed, and 4) carefully site a suitable number of monitoring wells and collect data at short intervals.

Climate-Informed Management of Irrigated Cotton in Western Kansas to Reduce Groundwater Withdrawals

The Ogallala aquifer, underlying eight states from South Dakota to Texas, is practically non-recharging south of Nebraska, and groundwater withdrawals for irrigation have lowered the aquifer in western Kansas. Subsequent well-yield declines encourage deficit irrigation, greater reliance on precipitation, and producing profitable drought-tolerant crops like upland cotton (Gossypium hirsutum (L.)). Our objective was to evaluate deficit irrigated cotton growth, yield, and water productivity (CWP) in northwest, west-central, and southwest Kansas in relation to El Niño southern oscillation (ENSO) phase effects on precipitation and growing season cumulative thermal energy (CGDD). Using the GOSSYM crop growth simulator with actual 1961–2000 location weather records partitioned by the ENSO phase, we modeled crop growth, yield, and evapotranspiration (ET) for irrigation capacities of 2.5, 3.75, and 5.0 mmd−1 and periods of 4, 6, and 8 weeks. Regardless of location, the ENSO phase did not influence CGDD, but precipitation and lint yield decreased significantly in southwest Kansas during La Niña compared with the Neutral and El Niño phases. Simulated lint yields, ET, CWP, and leaf area index (LAI) increased with increasing irrigation capacity despite application duration. Southwestern Kansas producers may use ENSO phase information with deficit irrigation to reduce groundwater withdrawals while preserving desirable cotton yields.

Assessment to locate potential deep aquifer systems using lithological logs, pumping tests and electromagnetic surveys in hard rock terrain of Dodoma urban area, Tanzania

The Dodoma region located in Central Tanzania is a semi-arid area characterised by scarce surface water resources. Because of climate change and dependence on groundwater resources in the Dodoma urban area, there is currently an increased demand for water due to low groundwater recharge and high groundwater withdrawals due to the growing population. There is therefore a need to explore potential deep aquifer systems around the Dodoma urban area to meet the increasing water demand. This study uses lithological logs, pumping tests, and electromagnetic survey data to delineate the subsurface lithologies to locate potential deep aquifer system in the area. Results show that the shallow aquifer system is unconfined in some areas with resistivity values ranging from 11 to 28 Ωm and semi-confined in other areas with resistivity values ranging from 19 to 27 Ωm. The unconfined aquifer extends up to 60 m, while the semi-confined aquifer extends between 50–55 m and 65–120 m. The study found a potential deep aquifer system at a depth of between 200 and 290 m, with resistivity values ranging from 11 to 20 Ωm. Currently, this deep aquifer system has not been exploited as most boreholes in the area are only up to 150 m deep. The main lithology in the study area is weathered and fractured granite, with different degrees of weathering and fracturing, indicating a hard-rock aquifer system. This study adds valuable knowledge on the location of potential deep aquifer systems in the area for proper groundwater utilisation and management.

Application of VIC-WUR model for assessing the spatiotemporal distribution of water availability in anthropogenically-impacted basins

Climate change, rapid increase in population, and urbanization have given rise to over-exploitation and shortage of freshwater resources. Acquiring accurate representations of fresh water availability is therefore crucial to sustainable development. Based on the underlying hydrological processes and storage types, fresh water can be divided into blue and green water components. Hence, quantification of the spatio-temporal distribution of blue and green water requires estimating the components of terrestrial water cycle. In this study, a framework is proposed to evaluate the water scarcity using the blue and green water concept. In this framework, the Wageningen Institute and Research (WUR) version of the Variable Infiltration Capacity Model (VIC) is employed to simulate terrestrial water cycle. Compared to VIC model, VIC-WUR model employs additional modules to account for sectoral water demands, return flows, and groundwater withdrawals, however, very limited number of its applications have been reported so far. The Divandareh-Bijar (DBJ) basin in west of Iran is selected because it has experienced severe anthropogenic alterations in the past three decades. To address precipitation and parameter uncertainty, the VIC-WUR model is separately calibrated for three precipitation datasets. As a part of this framework, the calibrated model is used to simulate blue and green water components and their corresponding scarcity metrics for the first time. According to the results, blue and green water resources in the DBJ are relatively scarce. The results also highlight the importance of considering anthropogenic impacts on the spatiotemporal dynamics of green and blue water resources. Additionally, quantification of the uncertainty in blue and green water components can help produce more reliable estimations for intelligent management of water resources. The ANalysis Of VAriance (ANOVA) method is therefore used to quantify the contribution of three major factors, including precipitation inputs, model parameters, and the interaction of these two, on the uncertainty of calculated water scarcity metrics. Structural uncertainty is implicitly addressed through the use of the modified version of the VIC model. The results indicate seasonality of these contributions for the two metrics under investigation and the differences in relative contribution of different sources to the overall uncertainty. The proposed framework could be employed to identify water stress in various spatiotemporal scales across the globe.

Mechanism controlling groundwater chemistry in the hyper-arid basin with intermittent river flow: insights from long-term observations (2001–2023) in the lower Heihe River, Northwest China

The geochemical processes of groundwater in arid regions are generally influenced by both natural hydrological processes and human activities. However, impacts of water-rock interactions on groundwater recharge via hydrological processes, controlled by both intermittent river water flow and groundwater withdrawals, is still poorly understood. In this study, 327 groundwater chemistry datasets collected from the upper, middle (including Gobi and riparian zones), and lower regions of the Ejina Delta in Northwest China from 2001 to 2023 were analyzed. Our results revealed that the total dissolved solids (TDS) concentration of groundwater in Ejina Delta ranged from approximately 881.5 ± 331.6 mg/L in the upper regions to 1,953.6 ± 1,208.5 mg/L in the lower regions, with an increasing trend observed. Ecological water conveyance (EWC), recharging aquifer through intermittent river water flow, resulted in a decrease in TDS concentrations from 2001 to 2023 mainly in the upper region. While irrigation notably affected groundwater chemistry in the lower region, resulting in a substantial increase in groundwater salinity. Groundwater chemistry in the Middle Gobi region remained relatively stable over the study period. Generally, the hydrochemical composition shifted from the Na-Mg-SO4-HCO3 and Na-Mg-Ca-SO4-HCO3 types in the upper region to Na-Mg-SO4-HCO3 and Na-Mg-SO4-Cl types in the lower region, with Na-SO4-Cl predominant in the Middle Gobi. These shifts were likely be attributed to the interplay of water-rock interactions, coupled with evaporation-crystallization processes. Inverse modeling using PHREEQC revealed that in the upper-middle region, primary water-rock interactions involved calcite dissolution and the precipitation of dolomite, gypsum, halite, and sylvite salts, as well as cation exchange reactions (2NaX+Ca2+→CaX2+2Na+). In contrast, the hydrogeological system in the middle-lower region exhibited an opposite pattern of water-rock interactions. Overall, ecological water conveyance partially facilitated water-rock interactions during lateral groundwater flow, while irrigation disrupted the natural hydrogeochemical equilibrium, involving halite dissolution and opposite cation exchange reactions compared to other regions.