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: It is proposed that a geothermal artesian karstic system at the central part of the Yarkon–Taninim aquifer creates the ‘Ayalon Saline Anomaly’ (ASA), whose mechanism has been under debate for several decades. A 4-year-long detailed groundwater monitoring was carried out at 68 new shallow boreholes in the Ayalon region, accompanied by a comprehensive survey of karstic voids. Results indicate the rising of warm-brackish groundwater through highly permeable swarms of karstic shafts, serving as an outflow of the artesian geothermal system. The ASA area contains ‘hot spots’, where groundwater contrasts with ‘normal’ water hundreds of meters away. The ASA temperature reaches 30°C (∼5°C warmer than its surroundings), chloride concentration reaches 528mg/l (50–100mg/l in the surrounding), H2S concentration reaches 5.6mg/l (zero all around) and pH value is 7.0 (compared with 7.8 around). Subsequently, the hydrothermal water flows laterally of at the watertable horizon through horizontal conduits, mixing with ‘normal’ fresh water which had circulated at shallow depth. Following rainy seasons, maximal watertable rise is observed in the ASA compared to its surroundings. Regional hydrogeology considerations suggest that the replenishment area for the ASA water is at the Samaria Mountains, east of the ASA. The water circulates to a great depth while flowing westward, and a cross-formational upward flow is then favored close the upper sub-aquifer’s confinement border.

Abstract: Study region Wadi Tayyibah is located in south Sinai, Egypt, in a region called Abou Zenima, and it is used to develop this study. Study focus Flash floods tremendously impact many facets of human life due to their destructive consequences and the costs associated with mitigating efforts. This study aims to evaluate the harvesting of Runoff by delineating the watersheds using the Hydrologic Engineering Center-1 (HEC-1) model and ArcGIS software in trying to benefit from it in different ways. All morphometric parameters of the basin were considered, and the risk degree of the different sub-basins was determined. The suitable locations of dams were identified using a Geographical Information System (GIS) using the basin’s morphometric characteristics. New hydrological insights for the region The study proposed a total number of eight dams, including five dams that were recommended for sub-basin (1) and three dams in sub-basin (4), while sub-basins (2) and (3) are not suitable locations to build dams according to the contour map of Wadi Tayyibah. Results indicate that, based on the constructed flash flood hazard maps and the basin’s detailed morphometric characteristics, the best locations of dams are Dam (3) in sub-basin (1) and Dam (7) in sub-basin (4), where the runoff volume reached 3.13 million cubic meters (Mm3) and 5.56 Mm3 for return period 100, respectively. This study is useful for decision-makers and designers for using morphometric parameters and flash flood hazard degree maps to select dam locations. Also, the cost-benefit analysis for using the morphometric parameters is required to be investigated.

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.

Abstract: The migration process of leached uranium in the in-situ leaching of uranium is considered a typical reactive transport problem. During this process, the lixiviant concentration and injection rate are important in-situ leaching process parameters (ISLPP) to efficiently recover uranium. However, several uncertain factors affect the outcomes of the ISLPP design. In addition, the repeated use of the reactive transport model (RTM) for investigating the acid in-situ leaching of uranium with the application of the Monte Carlo method leads to a substantial computational load. For this reason, a machine learning (ML)–based surrogate model was developed with the backpropagation neural network (BPNN) method to replace the RTM under the condition of uncertain parameters. Moreover, the simulated annealing optimisation model for ISLPP was created based on the proposed surrogate model. The optimal ISLPP was achieved that generated maximum profits from uranium recovery under different lixiviant prices, uranium prices and exploitation times. The optimal design framework of ISLPP based on the proposed ML algorithm was then applied in the Bayan-Uul sandstone-type uranium deposit in Inner Mongolia, China. From our analysis, it was demonstrated that the ML-based surrogate model exhibited great fitness with the RTM. The optimal results of the ISLPP indicated that the lixiviant concentration and injection rate could be adjusted based on the fluctuations in lixiviant price, uranium price and exploitation time. If the prices of sulphuric acid were high, a specific concentration of hydrogen peroxide could be injected into the injection well to promote the oxidation and dissolution of the uranium ore to increase the income from the uranium recovery. The optimisation model can also use the ISLPP scheme to boost the revenues from different lixiviant prices, uranium prices and exploitation times under the uncertainty of porosity, illustrating the applicability of the ML-based optimal design method of ISLPP in ISL mining. A general framework for developing surrogate models, as well as for conducting uncertainty analyses for a wide range of groundwater models was proposed here yielding valuable insights.