Abstract: Study Region: Bandung groundwater basin, Indonesia. Study focus: Groundwater abstraction of various magnitudes, pumped out from numerous depths in a multitude of layers of aquifers, stimulates different changes in hydraulic head distribution, including ones under vertical cross-sections. This generates groundwater flow in the vertical direction, where groundwater flows within its storage from the shallow to the underlying confined aquifers. In the Bandung groundwater basin, previous studies have identified such processes, but quantitative evaluations have never been conducted, with data scarcity mainly standing as one of the major challenges. In this study, we utilize the collated (1) environmental water tracer data, including major ion elements (Na+/K+, Ca2＋, Mg2+, Cl−, SO42−,HCO3−), stable isotope data (2H and δ18O), and groundwater age determination (14C), in conjunction with (2) groundwater flow modeling to quantify the aquifer interaction, driven mainly by the multi-layer groundwater abstraction in the Bandung groundwater basin, and demonstrate their correspondence. In addition, we also use the model to quantify the impact of multi-layer groundwater abstraction on the spatial distribution of the groundwater level changes. New hydrological insights for the region: In response to the limited calibration data availability, we expand the typical model calibration that makes use of the groundwater level observations, with in-situ measurement and a novel qualitative approach using the collated environmental water tracers (EWT) data for the model evaluation. The analysis in the study area using EWT data and quantitative methods of numerical groundwater flow modeling is found to collaborate with each other. Both methods show agreement in their assessment of (1) the groundwater recharge spatial distribution, (2) the regional groundwater flow direction, (3) the groundwater age estimates, and (4) the identification of aquifer interaction. On average, the downwelling to the deeper aquifer is quantified at 0.110 m/year, which stands out as a significant component compared to other groundwater fluxes in the system. We also determine the unconfined aquifer storage volume decrease, calculated from the change in the groundwater table, resulting in an average declining rate of 51 Mm3/year. This number shows that the upper aquifer storage is dwindling at a rate disproportionate to its groundwater abstraction, hugely influenced by losses to the deeper aquifer. The outflow to the deeper aquifer contributes to 60.3% of the total groundwater storage lost, despite representing only 32.3% of the total groundwater abstraction. This study shows the possibility of quantification of aquifer interaction and groundwater level change dynamics driven by multi-layer groundwater abstraction in a multi-layer hydrogeological setting, even in a data-scarce environment. Applying such methods can assist in deriving basin-scale groundwater policies and management strategies under the changing anthropogenic and climatic factors, thereby ensuring sustainable groundwater management.

Abstract: This study combines historical records of 14C and 3H in the atmosphere and soil with renewal rate and groundwater lumped parameter models to predict the abundance of 14C and 3H in groundwater over time. 624 groundwater samples from numerous studies, over four decades (1978–2019), in South Africa, Namibia, Botswana and Mozambique were collated to compare with predicted groundwater activities of 14C and 3H within the South African Development Community (SADC) region. Spatial datasets of carbonate bearing lithology, C3/C4 vegetation, summer/winter rainfall and coastal proximity were used to apply corrections to 14C and 3H data. Corrected values of 14C and 3H were compared with the theoretical abundance of these tracers, derived from the lumped parameter models, to estimate the general mean residence times and presence of groundwater mixing between modern recharge and older groundwaters. This study found that corrected values produced varying mean residence times derived from 14C ages (∼500–28500 years) and a wide range of potentially mixed waters within each aquifer system (0–100 % of tested wells) across the study area. The largest proportions of mixed groundwater, as well as the youngest mean residence times, were found in alluvial and primary fractured rock aquifers (e.g., western coast of South Africa and southern Mozambique). The smallest proportions of mixed groundwater were predicted in deep confined clay-rich aquifers as well as layered coal bearing carbonate sequences (e.g., Orapa, Malwewe and Serowe, Botswana). Insights into the proportions of mixed groundwater and mean residence times can help assess hydrological resilience on a regional scale. Such information is pertinent in promoting socio-economic development and increased water/food security in the SADC region. By understanding the resilience of groundwater resources, robust and informed strategies for water equality and GDP growth in the SADC region can be envisioned and implemented.

Abstract: Study region this study focuses on Madagascar. This island is characterized by a great diversity of climate, due to trade winds and the varying topography. This country is also undergoing extreme rainfall events such as droughts and cyclones. Study focus the rain gauge network of Madagascar is limited (about 30 stations). Consequently, we consider relevant satellite-based precipitation datasets to fill gaps in ground-based datasets. We assessed the reliability of 21 satellite-based and reanalysis precipitation products (P-datasets) through a direct comparison with 24 rain gauge station measurements at the monthly time step, using four statistical indicators: Kling-Gupta Efficiency (KGE), Correlation Coefficient (CC), Root Mean Square Error (RMSE), and Bias. Based on this first analysis, we produced a merged dataset based on a weighted average of the 21 products. New hydrological insights for the region based on the KGE and the CC scores, WFDEI (WATCH Forcing Data methodology applied to ERA-Interim), CMORPH-BLD (Climate Prediction Center MORPHing satellite-gauge merged) and MSWEP (Multi-Source Weighted Ensemble Precipitation) are the most accurate for estimating rainfall at the national scale. Additionally, the results reveal a high discrepancy between bio-climatic regions. The merged dataset reveals higher performance than the other products in all situations. These results demonstrate the usefulness of a merging approach in an area with a deficit of rainfall data and a climatic and topographic diversity.