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: Qanat is a type of drain that extract water from aquifers by gravity. Significant amount of fresh water used in Middle East and other parts of the world are supplied by qanats. Despite their importance, discharge variation of these type of wells received almost no attention. The aim of this research is to obtain a Laplace domain solution of discharge variation of a qanat installed in an anisotropic unconfined aquifer subjected to arbitrary areal recharge and nearby pumping well(s) discharge. A new semi-analytical solution of drawdown is obtained first to implement the effects of arbitrary areal recharge and nearby pumping well(s) using the principle of superposition. Then, the discharge variation solution of the qanat is obtained from the drawdown solution. To establish a constant-head boundary condition at the qanat periphery, the qanat is discretized into several segments. The results of this study are presented in dimensionless discharge-dimensionless time curves. The effects of hydraulic as well as geometric parameters on the discharge variation of the qanat due to arbitrary areal recharge, falling of water table from its initial position and discharge of nearby wells are explored. We also investigate the influences of distance and screen depth and location of the nearby well on the discharge variation of the qanat. The results of this study can be utilized for multiple purposes: 1) to predict discharge of qanat in response to rainfall and nearby pumping well(s); 2) to estimate the aquifer parameters using hydrograph of the qanat; 3) to determine optimal location and pumping pattern of the nearby wells to minimize their influences on the discharge of the qanat; 4) to calculate water budget of aquifers drained by a qanat. The equation presented in this work can also be used to estimate discharge of a horizontal drain installed in cropland subjected to arbitrary irrigation pattern.