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: The hydrogeological functioning of four different areas in a complex evaporite-karst unit of predominantly aquitard behavior in S Spain was investigated. Environmental dating tracers (3H, 3He, 4He, CFC-12, SF6) and hydrochemical data were determined from spring samples to identify and characterize groundwater flow components of different residence times in the media. Results show a general geochemical evolution pattern, from higher (recharge areas) to lower positions (discharge areas), in which mineralization rises as well as the value of the rCl−/SO42−, evidencing longer water-rock interaction. Ne values show degassing of most of the samples, favored by the high salinity of groundwater and the development of karstification so that the concentration of all the considered gases were corrected according to the difference between the theoretical and the measured Ne. The presence of modern groundwater in every sample was proved by the detection of 3H and CFC-12. At the opposite, the higher amount of radiogenic 4He in most samples also indicates that they have an old component. The 3H/3He dating method does not give reliable ages as a consequence of degassing and the large uncertainty of the 3He/4He ratios of the sources for the radiogenic Helium. The large SF6 concentrations suggest terrigenic production related to halite and dolomite. Binary Mixing and Free Shape Models were created based on 3H and CFC-12 data to interpret the age distribution of the samples. Two parameters (GA50 and >70%) were proposed as an indicator of that distribution, as they provide further information than the mean age. Particularly, GA50 is derived from the median groundwater age and is presented as a new way of interpreting mixed groundwater age data. A greater fraction of old groundwater (3H and CFC-12 free) was identified in discharge areas, while the proportion and estimated infiltration date of the younger fractions in recharge areas were higher and more recent, respectively. The application of different approaches has been useful to corroborate previous theoretical conceptual model proposed for the study area and to test the applicability of the used environmental tracer in dating brine groundwater and karst springs.