Abstract: In this paper, a new numerical model has been developed for simulation of Qanat-aquifer system. This model employs quasi-3D mixed-form of Richards’ equation and 1D fully-hydrodynamic form of Saint-Venant equations to simulate subsurface and overland flow, respectively. In order to handle non-orthogonal grids, subsurface flow module benefits from coordinate transformation technique. Using the above-mentioned governing equations, the presented model is able to simulate water flow inside both collection and conveyance sections of the gallery as well as dynamics of groundwater and vadose zone from impermeable bed rock to the soil-air interface. Since measured data corresponding to the hydraulics of Qanats is scarce, the overland and subsurface modules have been validated with analytical, numerical and experimental benchmarks in the literature. Subsequently, the model was employed to simulate ten different hypothetical aquifer-Qanat systems with different properties including the depth of groundwater aquifer, roughness of the gallery and saturated hydraulic conductivity of the gallery-aquifer boundary and the influence of each the parameters was monitored on the outflow rate at the appearance point of each Qanat. Furthermore, the advance of water inside two initially dry galleries were simulated at different time levels up to steady state. Eventually, the streamlines have been shown at the steady state for two Qanat-aquifer systems. Although, the presented study sheds light on some aspects of Qanat-aquifer hydraulics, the validation of the presented model with in-lab or on-field data remains ongoing for the future researches.

Abstract: Krypton-85 and other radioactive noble gases are widely used for groundwater dating purposes. 85Kr analysis require large volumes of water to reach the analytical requirements. Conventionally, this water is pumped to the surface to be degassed with a gas extraction system. The large pumping rate may disturb the natural flow field and requires substantial field logistics. Hence, we propose a new in-situ degassing method, in which membrane contactors are used to degas the groundwater directly in the well and gas is collected at the surface. This way, field work is facilitated, groundwater system disturbance is minimized, and the gas sample is collected at a specific depth. We demonstrate the tightness of the system regarding atmospheric air contamination for a collection times of 24 h, which is sufficient for both low-level counting and laser-based counting methods for 85Kr. The minimal borehole diameter is 7.5 cm for the prototype presented in this research but can easily be reduced to smaller diameters. In a case study, we compare the results obtained with the new passive method with those from a conventional packer setup sampling. Additionally, 3H/3He samples were collected for both sampling regimes and the dating results were compared with those from 85Kr. A good agreement between tracer ages is demonstrated and the age stratigraphy is consistent with the expected age distribution for a porous unconfined aquifer. In addition, our study emphasizes the differences between the age information sampled with various methods. In conclusion, we demonstrate that the new in situ quasi-passive method provides a more representative age stratigraphy with depth in most cases.

Abstract: Many studies highlight the decrease in precipitation due to climate change in the Mediterranean region, making it a prominent hotspot. This study examines the combined impacts of climate change and three groundwater demand scenarios on the water resources of the Western Mountain Aquifer (WMA) in Israel and the West Bank. While commonly used methods for quantifying groundwater recharge and water resources rely on regression models, it is important to acknowledge their limitations when assessing climate change impacts. Regression models and other data-driven approaches are effective within observed variability but may lack predictive power when extrapolated to conditions beyond historical fluctuations. A comprehensive assessment requires distributed process-based numerical models incorporating a broader range of relevant physical flow processes and, ideally, ensemble model projections. In this study, we simulate the dynamics of dual-domain infiltration and precipitation partitioning using a HydroGeoSphere (HGS) model for variably saturated water flow coupled to a soil-epikarst water balance model in the WMA. The model input includes downscaled high-resolution climate projections until 2070 based on the IPCC RCP4.5 scenario. The results reveal a 5% to 10% decrease in long-term average groundwater recharge compared to a 30% reduction in average precipitation. The heterogeneity of karstic flow and increased intensity of individual rainfall events contribute to this mitigated impact on groundwater recharge, underscoring the importance of spatiotemporally resolved climate models with daily precipitation data. However, despite the moderate decrease in recharge, the study highlights the increasing length and severity of consecutive drought years with low recharge values. It emphasizes the need to adjust current management practices to climate change, as freshwater demand is expected to rise during these periods. Additionally, the study examines the emergence of hydrogeological droughts and their propagation from the surface to the groundwater. The results suggest that the 48-month standardized precipitation index (SPI-48) is a suitable indicator for hydrogeological drought emergence due to reduced groundwater recharge.