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

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

Using hydrochemical,stable isotope,and river water recharge data to identify groundwater flow paths in a deeply buried karst system

The deeply buried river‐connected Xishan karst aquifer (XKA) in western Beijing, China, has been suffering from diminishing recharge for several decades, which in turn leads to the disappearing of spring water outflows and continuously lowering of groundwater level in the area. Thus, it is important to correctly recognize the groundwater recharge and flow paths for the sustainable development of the XKA. To investigate these issues, the hydrochemical and isotopic compositions are analysed for both surface water and groundwater samples collected over an area of about 280 km². Results show that (a) the river water is characterized by high Na contents; (b) the δ²H and δ¹⁸O values in the river water are distinctively higher than those of groundwater samples, after experiencing the long‐time evaporative enrichment in the upstream reservoir; (c) the Sr concentrations and ⁸⁷Sr/⁸⁶Sr ratios of groundwater clearly indicated the interaction between water and carbonate minerals but excluded the water–silicate interaction; and (d) the groundwater samples in the direct recharge area of the XKA have the lowest Na concentrations and the δ²H and δ¹⁸O values. Based on the large differences in the Na contents and ¹⁸O values of groundwater and surface water, a simple two‐component mixing model is developed for the study area and the fractions of the river water are estimated for groundwater samples. We find that the distribution pattern of the river water fractions in the XKA clearly shows a change of directions in the preferential flow path of the groundwater from its source zone to the discharge area. Overall, our results suggest that the recharged surface water can be a useful evidence for delineating the groundwater flow path in river‐connected karst aquifer. This study improves our understanding of the heterogeneity in karst groundwater systems.     

Estimating groundwater recharge for an arid karst system using a combined approach of time‐lapse camera monitoring and water balance modelling

Groundwater is the principal water resource in semi‐arid and arid environments. Therefore, quantitative estimates of its replenishment rate are important for managing groundwater systems. In dry regions, karst outcrops often show enhanced recharge rates compared with other surface and sub‐surface conditions. Areas with exposed karst features like sinkholes or open shafts allow point recharge, even from single rainfall events. Using the example of the As Sulb plateau in Saudi Arabia, this study introduces a cost‐effective and robust method for recharge monitoring and modelling in karst outcrops. The measurement of discharge of a representative small catchment (4.0 · 10⁴ m²) into a sinkhole, and hence the direct recharge into the aquifer, was carried out with a time‐lapse camera. During the monitoring period of two rainy seasons (autumn 2012 to spring 2014), four recharge events were recorded. Afterwards, recharge data as well as proxy data about the drying of the sediment cover are used to set up a conceptual water balance model. The model was run for 17 years (1971 to 1986 and 2012 to 2014). Simulation results show highly variable seasonal recharge–precipitation ratios between 0 and 0.27. In addition to the amount of seasonal precipitation, this ratio is influenced by the interannual distribution of rainfall events. Overall, an average annual groundwater recharge for the doline (sinkhole) catchment on As Sulb plateau of 5.1 mm has estimated for the simulation period. Copyright © 2015 John Wiley & Sons, Ltd.     

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

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

A model to simulate recharge processes of karst massifs

The recharge processes have been evaluated for two karst massifs of southern Italy, the Mt Terminio and Mt Cervialto, characterized by wide endorheic areas. The annual mean recharge has been estimated by Geographic Information System (GIS) tools, from regression of annual mean values of different ground‐elevated rain gauges and thermometers. The recharge has been distinguished for endorheic areas and the other areas of spring catchment, and the ratio between the output spring and input rainfall has been also estimated (recharge coefficient). The annual recharge has been used to calibrate a daily scale model, which allows to estimate the amount of effective rainfall, which is retained as soil moisture; the amount reaching the water table (recharge s.s.); and the amount of rainfall, which develops the runoff and leaves the catchment. All these amounts vary through the hydrological year, in function of soil moisture deficit and daily rainfall intensity. The model allows estimating the recharge conditions through the hydrological year, and it is a useful tool for water management. Copyright © 2014 John Wiley & Sons, Ltd.     

Estimating groundwater recharge using the chloride mass-balance method in the West Bank,Palestine

Abstract

The quantification of natural recharge rate is a prerequisite for efficient and sustainable groundwater resources management. Since groundwater is the only source of water supply in the West Bank, it is of utmost importance to estimate the rate of replenishment of the aquifers. The chloride mass-balance method was used to estimate recharge rates at different sites representing the three groundwater basins of the Mountain Aquifer in the West Bank. The recharge rate for the Eastern Basin was calculated as between 130.8 and 269.7 mm/year, with a total average replenishment volume of 290.3 × 10⁶ m³/year. For the Northeastern Basin, the calculated recharge rate ranged between 95.2 and 269.7 mm/year, with a total average recharge volume of 138.5 × 10⁶ m³/year. Finally, the recharge rate for the Western Basin was between 122.6 and 323.6 mm/year, with a total average recharge volume of 324.9 × 10⁶ m³/year. The data reveal a replenishment potential within the estimated replenishment volumes of previous studies for the same area. Also, the range was between 15 and 50% of total rainfall, which is still within the range of previous studies. The geological structure and the climate conditions of the western slope were clearly play an important role in the increment of total volume. In some cases, such as the geological formations in the Northeastern Basin, the interaction between Eocene and Senonian chalk formations result in minimum recharge rates.

Citation Marei, A., Khayat, S., Weise, S., Ghannam, S., Sbaih, M. & Geyer, S. (2010) Estimating groundwater recharge using the chloride mass-balance method in the West Bank, Palestine. Hydrol. Sci. J. 55(5), 780–791.     

Estimating groundwater recharge in lowland watersheds

Little attention has been given to the role of groundwater in the hydrological cycle of lowland watersheds. Our objective in this study was to estimate total recharge to groundwater by analysing water table response to storm events and the rate at which water was transferred into the shallow aquifer. This was conducted at three sites in a rural watershed in the lower Atlantic coastal plain near Charleston, South Carolina, USA. A novel version of the water table fluctuation method was used to estimate total recharge to the shallow aquifer by comparing hourly data of water table position following storm events and measuring water table recession behavior, rather than subjective graphical analysis methods. Also, shallow aquifer recharge rates (vertical fluxes) were estimated using Darcy's Law by comparing static water levels in a water table well and in a shallow piezometer during dry periods. The total annual recharge estimated ranged from 107 ± 39 mm·yr^–1 (5–10% of annual precipitation) at a poorly drained topographic low area to 1140 ± 230 mm·yr^–1 (62–94% of annual precipitation) for a moderately well‐drained upland site. The average aquifer recharge rate was 114 ± 60 mm·yr^–1, which is similar to previous estimations of base flow for the ephemeral third‐order streams in this watershed. The difference in the two methods may have been caused by processes not accounted for in the Darcy flux method, soil moisture deficits, and average evapotranspiration demand, which is about 1100 mm·yr^–1 for this region. Although other factors also can affect partitioning of recharge, an integrated approach to inspecting easily gathered groundwater data can provide information on an often neglected aspect of water budget estimation. We also discuss the effects of land use change on recharge reduction, given a typical development scenario for the region. Copyright © 2011 John Wiley & Sons, Ltd.     

Estimating parameters of groundwater recharge model in frequency domain: Karst springs Jadro and Žrnovnica

The spectral method is based on the assumption of exponential forms of autocorrelation functions of the input and output time series. The results of practical applications of the spectral method on several karst springs in Croatia show that this assumption cannot always be reliably applied to the time series of rainfall and spring discharge, especially if short time periods are analysed. For this reason, the method is modified. The scale factor between the transfer function of total rainfall and the transfer function of groundwater recharge is considered as an additional model parameter which is determined in the optimization procedure. In addition to the theoretical background, this paper also presents the results of application of the modified spectral method on the springs Jadro and ?rnovnica located in the Dinaric karst area in Croatia. The parameters of groundwater recharge model are determined separately for both the springs and the obtained values are discussed and compared with the results of previous investigations. The values of the scale factors determined in the optimization procedure are compared to the values calculated using the theoretical equation. The comparison shows that the theoretical equation underestimates the optimal values of the scale factor. The underestimations are especially evident for periods of one year. Finally, the characteristics of rainfall–runoff relations of two springs are compared using the identified transfer functions and the results of simulations of the periodogram of observed discharge by the parametric periodogram. Copyright © 2008 John Wiley & Sons, Ltd.     

Geochemical modeling of groundwater in karst area and its application at Pingdingshan coalfield

Several approaches including hydrogeochemistry and isotope hydrogeology have been used to investigate the karst groundwater systems at Pingdingshan coalfield in recent years. The results of the modeling and evaluation of groundwater chemistry, as parts of recent research progress at the area, are presented. The characteristics of Cambrian karst groundwater flow was analyzed in terms of tritium distribution of groundwater based on the fact that the Guodishan fault, the largest one within the coalfield, is divided into permeable (southern and northern) and impermeable (middle) sections. The evolution of groundwater chemistry, the suitability of geothermometers and the feature of karst development were deduced and discussed using the speciation modeling and mass balance approach. Project partially supported by China Postdoctoral Science Foundation.     

Reliability of recharge rates estimated from groundwater age with a simplified analytical approach

Groundwater age is often used to estimate groundwater recharge through a simplified analytical approach. This estimated recharge is thought to be representative of the mean recharge between the point of entry and the sampling point. However, given the complexity in actual recharge, whether the mean recharge is reasonable is still unclear. This study examined the validity of the method to estimate long-term average groundwater recharge and the possibility of obtaining reasonable spatial recharge pattern. We first validated our model in producing reasonable age distributions using a constant flux boundary condition. We then generated different flow fields and age patterns by using various spatially varying flux boundary conditions with different magnitudes and wavelengths. Groundwater recharge was estimated and analysed afterwards using the method at the spatial scale. We illustrated the main findings with a field example in the end. Our results suggest that we can estimate long-term average groundwater recharge with 10% error in many parts of an aquifer. The size of these areas decreases with the increase in both the amplitude and the wavelength. The chance of obtaining a reasonable groundwater recharge is higher if an age sample is collected from the middle of an aquifer and at downstream areas. Our study also indicates that the method can also be used to estimate local groundwater recharge if age samples are collected close to the water table. However, care must be taken to determine groundwater age regardless of conditions.     

A numerical investigation of bedrock groundwater recharge and exfiltration on soil mantled hillslopes

Here we use Richards Equation models of variably saturated soil and bedrock groundwater flow to investigate first-order patterns of the coupling between soil and bedrock flow systems. We utilize a Monte Carlo sensitivity analysis to identify important hillslope parameters controlling bedrock recharge and then model the transient response of bedrock and soil flow to seasonal precipitation. Our results suggest that hillslopes can be divided into three conceptual zones of groundwater interaction, (a) the zone of lateral unsaturated soil moisture accumulation (upper portion of hillslope), (b) the zone of soil saturation and bedrock recharge (middle of hillslope) and (c) the zone of saturated-soil lateral flow and bedrock groundwater exfiltration (bottom of hillslope). Zones of groundwater interaction expand upslope during periods of precipitation and drain downslope during dry periods. The amount of water partitioned to the bedrock groundwater system a can be predicted by the ratio of bedrock to soil saturated hydraulic conductivity across a variety of hillslope configurations. Our modelled processes are qualitatively consistent with observations of shallow subsurface saturation and groundwater fluctuation on hillslopes studied in our two experimental watersheds and support a conceptual model of tightly coupled shallow and deep subsurface circulation where groundwater recharge and discharge continuously stores and releases water from longer residence time storage.     

A combined methodology for estimating the potential natural aquifer recharge in an arid environment

ABSTRACT

An innovative methodology that combines an indirect physiography-based method for determining the runoff coefficient at a sub-basin scale and a water balance model applied on a daily time scale was developed to calculate the natural groundwater recharge in three watersheds within the Oum Zessar arid area, Tunisia. The effective infiltration was calculated as part of the water surplus by considering the average available water content (AWC) of soil and an average runoff coefficient for each sub-basin. The model indicates that the sub-basins covered mainly by the “artificial” soils of tabias and jessour, characterized by average AWC values greater than 150 mm, did not contribute to natural groundwater recharge over the 10-year period (2003–2012) considered. The estimated volume for the Triassic aquifer amounted to about 4.5 hm³ year^?1, which is consistent with previous studies. For the Jurassic and Cretaceous aquifers, the estimated volumes amounted to about 200 dm³ year^?1.     

A groundwater recharge field study: site characterization and initial results

A field study site was installed in east‐central Pennsylvania to examine processes controlling groundwater recharge. It was instrumented to monitor climatic inputs, soil water dynamics and groundwater response. Characterization of the layered fractured bedrock underlying the site by rock coring, seismic surveys and interval packer testing showed consistencies between layer depths, fracture frequencies, seismic velocities and hydraulic conductivities. Monthly summaries of rainfall and percolate over two years showed that percolate rates were generally high and closely related to precipitation during the dormant season. During the growing season, however, the relationship became erratic with large variabilities occurring between individual lysimeter measurements. Eight dormant season rainfall events were examined in detail. Smaller events produced similar responses from 1 m deep percolate lysimeters. Approximately 10–15 mm of rain was required to initiate percolate, with the time delay in response dependent on how long it took this depth to accumulate; 5 to 6 mm of the rain was retained in storage, with the remainder becoming percolate. Larger rains, from 30–110 mm, caused correspondingly larger depths of percolate and larger water table responses, but generally similar patterns of site response. Groundwater at the site was typically about 6 m below the land surface during the dormant season. It responded 1–2 hours after the onset of percolate, and reached its maximum elevation anywhere from 4 to 16 hours after that, even though percolate was still occurring. Based on causative depth of recharge and amount of water level rise in wells, the specific yield of the aquifer was found to be of the order of 0·01. This value is characteristic of fracture geometry rather than matrix properties of the bedrock. Copyright © 1999 John Wiley & Sons, Ltd.     

Assessment of CHADFDM satellite-based input dataset for the groundwater recharge estimation in arid and data scarce regions

Aquifer natural recharge estimations are a prerequisite for understanding hydrologic systems and sustainable water resources management. As meteorological data series collection is difficult in arid and semiarid areas, satellite products have recently become an alternative for water resources studies. A daily groundwater recharge estimation in the NW part of the Lake Chad Basin, using a soil–plant-atmosphere model (VisualBALAN), from ground- and satellite-based meteorological input dataset for non-irrigated and irrigated land and for the 2005–2014 period is presented. Average annual values were 284 mm and 30°C for precipitation and temperature in ground-based gauge stations. For the satellite-model-based Lake Chad Basin Flood and Drought Monitor System platform (CHADFDM), average annual precipitation and temperature were 417 mm and 29°C, respectively. Uncertainties derived from satellite data measurement could account for the rainfall difference. The estimated mean annual aquifer recharge was always higher from satellite- than ground-based data, with differences up to 46% for dryland and 23% in irrigated areas. Recharge response to rainfall events was very variable and results were very sensitive to: wilting point, field capacity and curve number for runoff estimation. Obtained results provide plausible recharge values beyond the uncertainty related to data input and modelling approach. This work prevents on the important deviations in recharge estimation from weighted-ensemble satellite-based data, informing in decision making to both stakeholders and policy makers.     

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

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

A numerical model for estimating groundwater flux from subsurface temperature profiles

A simple numerical model is presented for estimating vertical groundwater flux from transient subsurface temperature profiles obtained from field measurements. The model developed utilizes the MacCormack scheme, which is based on the Finite Difference Method (FDM), for solving the governing partial differential equation of convection–diffusion heat transport with appropriate initial and boundary conditions within the subsurface. In order to validate the model, numerical solutions obtained for the study area located in the Nagoka plain, Japan are compared with the published measured data and results obtained by others. Results obtained show good agreement and fit the observed data with a correlation coefficient, R², of 0·88. The estimated groundwater flux is 1·85 × 10⁻⁷ m s⁻¹. Sensitivity analyses were also carried out to investigate the effect of variations in groundwater fluxes, thermal properties and the annual thermal variability due to climatic changes on the transient subsurface temperature profiles and to have a better understanding of the subsurface thermal dynamics. A substantial effect of annual climatic variability is observed on the temporal distributions of temperature depth profiles, and a better estimate of thermal parameters is required to estimate vertical groundwater flux. The largest change in subsurface temperature depth profiles due to groundwater flux over a year is within ± 4 °C. The influence of groundwater flux on subsurface temperature distributions in space and time may be more pronounced in areas where the top of the saturated layer fluctuates considerably. Variation in thermal diffusivity results in temperature change up to ± 1·5% and may cause change in groundwater flux estimate by ± 18%. The model presented has merits over analytical solutions (type curve matching techniques) in terms of suitability and applicability to real field problems, and can be a good asset to hydrological models as quantifying groundwater recharge or deducing it from other quantities, such as rainfall, evapotranspiration and runoff, is often complicated. Copyright © 2007 John Wiley & Sons, Ltd.     

Interpretation of groundwater level monitoring results in karst aquifers: examples from the Dinaric karst

The paper presents an attempt to determine the characteristics of karst aquifers using information on groundwater level (GWL) in natural holes and boreholes with different data quantity and time resolution of GWL measurements. In this paper the particulars of karst aquifers were analysed for four examples from the Dinaric karst. In all four study areas, aquifers are formed in bare, deep and well‐developed Dinaric karst consisting of Cretaceous limestones. The first example represents a wide area of Imotsko polje in the karst. The aquifer was analysed on the basis of infrequent water level monitoring in natural karst water features (jamas, lakes, wells) and discharges of springs and rivers. The karst aquifer in this example is complex, non‐homogenous and variable in space and time, which is frequent in the Dinaric karst. Regardless of the aforementioned it was possible to determine its elementary characteristics. The second example represents 10 wells used for the water supply for the city of Pula. The GWL and salinity were measured once a week in the period between 1981 and 1996. Even though these measurements were relatively infrequent in space and time, they served as bases for assessment of average and maximum aquifer conditions as well as boundaries of saltwater intrusion. In the third example only a portion of aquifer of the karst spring Blaz, which is in the contact with the Adriatic Seas, has been analyzed. It is a spring with an intrusion of salt water. For purposes of study of saltwater intrusion, 26 piezometers were drilled in its vicinity in which GWL, salinity and temperature were measured once a day during 168 days, a period comprising one complete cycle of seawater intrusion and retreat. These measurements proved the existence of dispersed discharge from the aquifer into the sea and its non‐homogeneity in space. In the fourth example GWL was measured continuously in 10 deep (up to 300 m) piezometers in the hinterland of the Ombla Spring catchment. The measurement period lasted 2 years (January 1988 to December 1989). The analyses are made with hourly data. The results made it possible to determine numerous characteristics of the karst aquifer and a significant non‐homogeneity of groundwater distribution in karst aquifers, depending more on the underground karst phenomena than the surface karst forms. Copyright © 2000 John Wiley & Sons, Ltd.     

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

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

Modelling regional effects of artificial groundwater recharge in a multilayer aquifer characterized by perched water tables

This paper presents an approach to estimate the effects of a managed recharge experiment in a multilayer aquifer characterized by the presence of perched water tables in the Medina del Campo groundwater body, Douro basin, central Spain. A numerical model was developed to evaluate the effect of artificial recharge on the shallow sector of a regional-scale aquifer and on formerly active wetlands. The model was developed in the Visual MODFLOW Pro v.2011.1 environment in order to represent and analyse the regional impact of this artificial recharge event. Results suggest that the assumption of a single perched system may prove useful in regional contexts where data is limited. From a study site perspective, managed recharge is observed to increase shallow storage along the riverbanks, which is considered valuable for environmental purposes. However, downstream wetlands are unlikely to experience a significant recovery. Furthermore, only a small percentage of artificial recharge is expected to reach the deep regional aquifer. This method can be exported to settings characterized by the presence of perched aquifers and associated groundwater dependent ecosystems.     

Unique episodic groundwater recharge event in a South American sedimentary aquifer and its long-term impact on baseflow

The anomalous entrance of water into groundwater systems can affect storage throughout long periods and normally relies on infrequent and irregular pulses of groundwater recharge defined by the term episodic recharge. Recently there was a groundwater recharge of large magnitude with unknown circumstances in the Caiuá aquifer. This unique event was explored in detail here and allowed to better understand the occurrence of such events in humid subtropical climates in South America. For this study, groundwater monitoring daily data from the Integrated Groundwater Monitoring Network was used combined with a specific yield obtained from geophysical wireline logging to obtain groundwater recharge rates. To improve the investigation, we also used a baseflow separation method to obtain the groundwater contribution into local rivers. The groundwater storage variations were also assessed by remote sensing with the GRACE data. Results showed the importance of high soil moisture storage on the occurrence of large episodic recharge events. We estimated that the groundwater recharger volumes derived from 1 year that included the unique episodic recharge observed (total of 866 mm for April 2015–March 2016) were comparable with the sum of 7 years of groundwater recharge (total of 867 mm). Atypical rainfall in winter periods were responsible for the increase in soil moisture that explained that unique event. GRACE-based GWS showed concordance detecting the occurrence of the unique episodic recharge. However, the variation in terms of volumes obtained by GRACE does not represent the behaviour observed in the aquifer by the WTF method. The results also indicated that changes in aquifer storage caused by episodic recharge events directly affect low flows in rivers over long periods. The main knowledge gap addressed here relates to exploring a unique episodic recharge event quite rare to observe with its long-term impacts on hydroclimatic variability over a humid subtropical portion of the Caiuá aquifer.