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.

Abstract: Numerous groundwater analyses from the crystalline bedrock in the Forsmark area have been performed between 2002 and 2019, together with thorough geological, geophysical, and hydrogeological studies, within the site investigations carried out by the Swedish Nuclear Fuel and Waste Management Company. The groundwater samples have been taken from boreholes down to ≈ 1000 m and the analysis include major- and trace-elements, stable and radiogenic isotopes, gases and microbes. The chemical and isotopic composition of these groundwaters evidences the presence of non-marine brackish to saline groundwaters with very long residence times (many hundreds of thousands of years) and a series of complex mixing events resulting from the recharge of different waters over time: glacial meltwaters, probably from different glaciations of which the latest culminated some 20,000 years ago, and marine waters from the Baltic starting some 7000 years ago. Later, meteoric water and present Baltic Sea water have recharged in different parts of the upper 100 m. These mixing events have also triggered chemical and microbial reactions that have conditioned some of the important groundwater parameters and, together with the structural complexity of the area, they have promoted a heterogeneous distribution of groundwater compositions in the bedrock. Due to these evident differences in chemistry, residence time and origin of the groundwater, several groundwater types were defined in order to facilitate the visualisation and communication. The differentiation (linked to the paleohydrological history of the area) was based on Cl concentration, Cl/Mg ratio (marine component), and δ18O value (glacial component). The work presented in this paper increases the understanding of the groundwater evolution in fractured and compartmentalised aquifers where mixing processes are the most important mechanisms. The model proposed to characterise the present groundwater system of the Forsmark area will also help to predict the future hydrogeochemical behaviour of the groundwater system after the construction of the repositories for the nuclear wastes.

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.