Abstract: Study region Zohreh River Basin, Southwest Iran Study focus The salinity of Zohreh River sharply increases in three salinity zones (SZs) along the river named SZ1, SZ2 (the focus of this study), and SZ3. Determining the salinity sources and salinization mechanism using an integrated approach including geological, hydrochemical, isotopic, geophysical, river sinuosity and hydrocarbon analysis are the main objectives of this study. The study focuses on the combination of evidence of regional-scale (i.e., river sinuosity and seismic data) and small-scale (i.e., drilling core analysis). New hydrologic insights for the region Among several known sources of river salinity, it was found that the water quality of the Zohreh River is mainly threatened by the salt-bearing Gachsaran Formation and oil-field brine. It is concluded that halite brine and oil-field brine simultaneously cause the salinization in SZ2, and their contributions were delineated to be 95% and 5%, respectively. The lack of reliable geological evidence to support halite dissolution in surficial layers by circulating waters suggests the possibility of a deep source of halite brine in SZ2. The results revealed that deep halite brine of the salt layers of Gachsaran Formation is mainly responsible for the salinization of SZ2. The mechanism of deep brine penetration to the river through the hidden fault failures detected by the combination of river sinuosity analysis and geophysical data for the first time.

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: 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.

Abstract: Drawing on advances in nonstationary frequency analysis and the science of causation and attribution, this study employs a newly developed nonstationary stochastic paired watershed approach to determine the effect of forest harvesting on snowmelt-generated floods. Moreover, this study furthers the application of stochastic physics to evaluate the environmental controls and drivers of flood response. Physically-based climate and time-varying harvesting data are used as covariates to drive the nonstationary flood frequency distribution parameters to detect, attribute, and quantify the effect of harvesting on floods in the snow-dominated Deadman River (878 km2) and nested Joe Ross Creek (99 km2) watersheds. Harvesting only 21% of the watershed caused a 38% and 84% increase in the mean but no increase in variability around the mean of the frequency distribution in the Deadman River and Joe Ross Creek, respectively. Consequently, the 7-year, 20-year, 50-year, and 100-year flood events became approximately two, four, six, and ten times more frequent in both watersheds. An increase in the mean is posited to occur from an increase in moisture availability following harvest from suppressed snow interception and increased net radiation reaching the snowpack. Variability was not increased because snowmelt synchronization was inhibited by the buffering capacity of abundant lakes, evenly distributed aspects, and widespread spatial distribution of cutblocks in the watersheds, preventing any potential for harvesting to increase the efficiency of runoff delivery to the outlet. Consistent with similar recent studies, the effect of logging on floods is controlled not only by the harvest rate but most importantly the physiographic characteristics of the watershed and the spatial distribution of the cutblocks. Imposed by the probabilistic framework to understanding and predicting the relation between extremes and their environmental controls, commonly used in the general sciences but not forest hydrology, it is the inherent nature of snowmelt-driven flood regimes which cause even modest increases in magnitude, especially in the upper tail of the distribution, to translate into surprisingly large changes in frequency. Contrary to conventional wisdom, harvesting influenced small, medium, and very large flood events, and the sensitivity to harvest increased with increasing flood event size and watershed area.

Abstract: In coastal salinized groundwater systems, contamination from various nitrate (NO3) inputs combined with complex hydrogeochemical processes make it difficult to distinguish NO3 sources and identify potential NO3 transformtation processes. Effective field-based NO3 studies in coastal areas are needed to improve the understanding of NO3 contamination dynamics in groundwater of such complex coastal systems. This study focuses on a typical Mediterranean coastal agricultural area, located in Tunisia, experiencing substantial NO3 contamination from multiple anthropogenic sources. Here, multiple isotopic tracers (δ18OH2O, δ2HH2O, δ15NNO3, δ18ONO3, and δ11B) combined with a Bayesian isotope MixSIAR model are used (i) to identify the major NO3 sources and their contributions, and (ii) to describe the potential NO3 transformation processes. The measured NO3 concentrations in groundwater are above the natural baseline threshold, suggesting anthropogenic influence. The measured isotopic composition of NO3 indicates that manure, soil organic matter, and sewage are the potential sources of NO3, while δ11B values constrain the NO3 contamination to manure; a finding that is supported by the results of MixSIAR model revealing that manure-derived NO3 dominates over other likely sources. Nitrate derived from manure in the study area is attributed to organic fertilizers used to promote crop growth, and livestock that deposit manure directly on the ground surface. Evidence for ongoing denitrification in groundwaters of the study area is supported by an enrichment in both 15N and 18O in the remaining NO3, although isotopic mass balances between the measured and the theoretical δ18ONO3 values also suggest the occurrence of nitrification. The simultaneous occurrence of these biogeochemical processes with heterogeneous distribution across the study area reflect the complexity of interactions within the investigated coastal aquifer. The multiple isotopic tracer approach used here can identify the effect of multiple NO3 anthropogenic activities in coastal environments, which is fundamental for sustainable groundwater resources management.