Efficient injection of gas tracers into rivers: A tool to study Surface water-Groundwater interactions.

Surface water (SW) - groundwater (GW) interactions exhibit complex spatial and temporal patterns often studied using tracers. However, most natural and artificial tracers have limitations in studying SW-GW interactions, particularly if no significant contrasts in concentrations between SW and GW exist or can be maintained for long durations. In such context, (noble) gases have emerged as promising alternatives to add to the available tracer methods, especially with the recent development of portable mass spectrometers, which enable continuous monitoring of dissolved gas concentrations directly in the field. However, long-duration gas injection into river water presents logistical challenges. To overcome this limitation, we present an efficient and robust diffusion-injection apparatus for labeling large amounts of river water. Our setup allows fine, real-time control of the gas injection rate, and is suitable for extended injection durations and different gas species. To illustrate the effectiveness of our approach, we present a case study where helium (He) is used as an artificial tracer to study river water infiltration into an alluvial aquifer. Our injection of He as a tracer increased the dissolved He concentration of the river water by one order of magnitude compared to air-saturated water concentration for 35 days. This experiment yields valuable information on travel times from the river to a pumping well and on the mixing ratios between freshly infiltrated river water and regional groundwater.

Evaluating the impact of muon-induced cosmogenic 39Ar and 37Ar underground production on groundwater dating with field observations and numerical modeling.

Groundwater dating by radioactive cosmogenic tracers such as 39Ar relies on the decay rate from a known initial atmospheric activity (100%modern). Thereby, it is assumed that cosmogenic 39Ar production in the subsurface is negligible at depths below the water table and that contributions from natural rock radioactivity are minor or missing. Here we present 39Ar data from aquifers located in quaternary glacial sediments and tertiary limestones in Denmark, which unequivocally demonstrate that cosmogenic production can induce considerable age biases. 39Ar values larger than 100%modern are observed at relatively shallow groundwater depths in non-radiogenic rocks. These activities are compared to calculations based on previously assessed depth-dependent production rates in rocks and realistic estimates of the emanated fractions to the water phase. The water residence time distribution with depth, which was determined by numerical flow modeling and particle tracking, underpinned the significance of muon-induced 39Ar production. The short-lived isotope 37Ar is produced by similar processes as 39Ar and demonstrated its usefulness as an indicator of local underground production in an aquifer. The significance of cosmogenic underground production in other possible recharge scenarios was then assessed by explicitly simulating the radioargon accumulation and decay in a 2D synthetical numerical model. These simulations demonstrated that underground production is negligible when the water infiltrates freely in a porous aquifer. However, in the presence of a confining layer impeding the infiltration at shallow depths (<30 m), as is the case in our study site in Denmark for instance, over-modern 39Ar activities (>100%modern) may occur. The age concluded from the dissolved activities is then possibly biased towards young values. Special attention should thus be paid to the recharge rates when using 39Ar for dating groundwater. 37Ar activities provide complementary information about the strength and mechanisms of underground production.

Using machine learning and remote sensing to track land use/land cover changes due to armed conflict.

Armed conflicts have detrimental impacts on the environment, including land systems. The prevailing understanding of the relation between Land Use/Land Cover (LULC) and armed conflict fails to fully recognize the complexity of their dynamics - a shortcoming that could undermine food security and sustainable land/water resources management in conflict settings. The Syrian portion of the transboundary Orontes River Basin (ORB) has been a site of violent conflict since 2013. Correspondingly, the Lebanese and Turkish portions of the ORB have seen large influxes of refugees. A major challenge in any geoscientific investigation in this region, specifically the Syrian portion, is the unavailability of directly-measured "ground truth" data. To circumvent this problem, we develop a novel methodology that combines remote sensing products, machine learning techniques and quasi-experimental statistical analysis to better understand LULC changes in the ORB between 2004 and 2022. Through analysis of the resulting annual LULC maps, we can draw several quantitative conclusions. Cropland areas decreased by 21-24 % in Syria's conflict hotspot zones after 2013, whereas a 3.4-fold increase was detected in Lebanon. The development of refugee settlements was also tracked in Lebanon and on the Syrian/Turkish borders, revealing different LULC patterns that depend on settlement dynamics. The results highlight the importance of understanding the heterogenous spatio-temporal LULC changes in conflict-affected and refugee-hosting countries. The developed methodology is a flexible, cloud-based approach that can be applied to wide variety of LULC investigations related to conflict, policy and climate.

Impacts of Subsurface Dam Construction on Downstream Groundwater Levels and Salinity in Coastal Aquifers.

Subsurface dam is a promising engineering technology for groundwater resources development. However, the possible impacts of these dams on the groundwater environment have been a major concern. Here, we used a three-dimensional (3D), variable-density, unsaturated-saturated groundwater flow model to explore how a groundwater-storage-type subsurface dam, built in the freshwater domain of an unconfined coastal aquifer, affected groundwater levels and salinity in the downstream area. Model results suggested that, after subsurface dam construction, groundwater levels in the downstream area showed intensified fluctuations in terms of phase advances, greater amplitudes, and higher frequencies following heavy rainfall events. Numerical simulations with variable subsurface dam scenarios indicated that the groundwater level fluctuations were further intensified with a higher crest elevation or a shorter distance from the coast. Moreover, during the recharging period of the subsurface reservoir, sea water in the downstream area intruded landward from its initial location, which can at least temporarily threaten water quality near the coast. A higher dam crest elevation prolonged the duration of sea water intrusion, while a dam positioned closer to the coast induced sea water intrusion with a greater horizontal extent. General implications are discussed with respect to improving assessment methodologies and engineering designs of subsurface dams.

Exploring the reliability of 222Rn as a tracer of groundwater age in alluvial aquifers: Insights from the explicit simulation of variable 222Rn production.

Knowledge of groundwater residence times (GRT; the time elapsed since surface water infiltration) between losing rivers and pumping wells is crucial for management of water resources in alluvial aquifers. The radioactive noble gas radon-222 (222Rn) has been used for decades as a natural indicator of surface water infiltration, as it can provide quantitative information on GRT. However, models using 222Rn as a tracer of GRT are often based on a set of highly simplifying assumptions, including spatially homogenous 222Rn production and exclusively advective mass transport within the aquifer. In this paper, we use the integrated surface-subsurface hydrological model HydroGeoSphere (HGS) to simulate 222Rn transport, production, and decay in a bank filtration context. Spatially variable 222Rn production, based on experimental data, is explicitly considered. We show that variable 222Rn production rates, coupled with hydrodispersive mixing of groundwater, may lead to large biases in GRT estimates. Under certain transient conditions however, changes in tracer-derived GRTs correlate well with changes in mean groundwater age. Whereas 222Rn-derived GRTs may only be reliable under a narrow range of field conditions, 222Rn may serve as a powerful tracer of changes in mean GRT even in complex and heterogenous environments.

Salix psammophila afforestations can cause a decline of the water table, prevent groundwater recharge and reduce effective infiltration.

Afforestation can reduce desertification and soil erosion. However, the hydrologic implications of afforestation are not well investigated, especially in arid and semi-arid regions. China has the largest area of afforestation in the world, with one-third of the world's total plantation forests. How the shrubs affect evapotranspiration, soil moisture dynamics, and groundwater recharge remains unclear. We designed two pairs of lysimeters, one being 1.2 m deep and the other one 4.2 m deep. Each pair consists of one lysimeter with bare soil, while on the other one a shrub is planted. The different water table depths were implemented to understand how depth to groundwater affects soil moisture and water table dynamics under different hydrological conditions. Soil moisture, water table depth, sap flow, and rainfall were measured concurrently. Our study confirms that for the current meteorological conditions in the Ordos plateau recharge is reduced or even prohibited through the large-scale plantation Salix psammophila. Shrubs also raise the threshold of precipitation required to increase soil moisture of the surface ground. For the conditions we analyzed, a minimum of 6 mm of precipitation was required for infiltration processes to commence. In addition to the hydrological analysis, the density of root distribution is assessed outside of the lysimeters for different water table depths. The results suggest that the root-density distribution is strongly affected by water table depth. Our results have important implications for the determination of the optimal shrub-density in future plantations, as well as for the conceptualization of plant roots in upcoming numerical models.

Publisher Correction: A 3D geological model of a structurally complex Alpine region as a basis for interdisciplinary research.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

A 2D hyperspectral library of mineral reflectance, from 900 to 2500 nm.

Mineral identification using machine learning requires a significant amount of training data. We built a library of 2D hyperspectral images of minerals. The library contains reflectance images of 130 samples, of 76 distinct minerals, with more than 3.9 million data points. In order to produce this dataset, various well-characterized mineral samples were scanned, using a SPECIM Short Wave Infrared (SWIR) camera, which captures wavelengths from 900 to 2500 nm. Minerals were selected to represent all the mineral classes and the most common mineral occurrences. For each sample, the following data are provided: (a) At least one hyperspectral image of the sample, consisting of 256 wavelengths between 900 and 2500 nm. The raw data, the high dynamic range (HDR) image, and the masked HDR image are provided for each scan (each of them in HDF5 format). (b) A text file describing the sample, providing supplementary information for the subsequent analysis (c) RGB images (JPEG files) and automated 3D reconstructions (Stanford Triangle PLY files). These data help the user to visualize and understand specific sample characteristics. 2D hyperspectral images were produced for each mineral, which consist of many different spectra with high diversity. The scans feature similar spectra than the ones in other available spectral libraries. An artificial neural network was trained to demonstrate the high quality of the dataset. This spectral library is mainly aimed at training machine learning algorithms, such as neural networks, but can be also used as validation data for other types of classification algorithms.

Conceptualization and Calibration of Anisotropic Alluvial Systems: Pitfalls and Biases.

Physical properties of alluvial environments typically feature a high degree of anisotropy and are characterized by dynamic interactions between the surface and the subsurface. Hydrogeological models are often calibrated under the assumptions of isotropic hydraulic conductivity fields and steady-state conditions. We aim at understanding how these simplifications affect predictions of the water table using physically based models and advanced calibration and uncertainty analysis approaches based on singular value decomposition and Bayesian analysis. Specifically, we present an analysis of the information content provided by steady-state hydraulic data compared to transient data with respect to the estimation of aquifer and riverbed hydraulic properties. We show that assuming isotropy or fixed anisotropy may generate biases both in the estimation of aquifer and riverbed parameters as well as in the predictive uncertainty of the water table. We further demonstrate that the information content provided by steady-state hydraulic heads is insufficient to jointly estimate the aquifer anisotropy together with the aquifer and riverbed hydraulic conductivities and that transient data can help to reduce the predictive uncertainty to a greater extent. The outcomes of the synthetic analysis are applied to the calibration of a dynamic and anisotropic alluvial aquifer in Switzerland (The Rhône River). The results of the synthetic and real world modeling and calibration exercises documented herein provide insight on future data acquisition as well as modeling and calibration strategies for these environments. They also provide an incentive for evaluation and estimation of commonly made simplifying assumptions in order to prevent underestimation of the predictive uncertainty.

Exploring Geological and Topographical Controls on Low Flows with Hydrogeological Models.

This study investigates how catchment properties influence low-flow dynamics. With 496 synthetic models composed of a bedrock and an alluvial aquifer, we systematically assess the impact of the hydraulic conductivity of both lithologies, of the hillslope and of the river slope on catchment dynamics. The physically based hydrogeological simulator HydroGeoSphere is employed, which allows obtaining a range of low-flow indicators. The hydraulic conductivity of the bedrock Kbedrock , a proxy for transmissivity, is the only catchment property exerting an overall control on low flows and explains 60% of the variance of Q95/Q50. The difference in dynamics of catchments with same Kbedrock depends on hillslope gradients and the alluvial aquifer properties. The buffering capacity of the bedrock is mainly related to Kbedrock and the hillslope gradient. We thus propose the dimensionless bedrock productivity index (BPI) that combines these characteristics with the mean net precipitation. For bedrock only models, the BPI explains 82% of the variance of the ratio of Q95 to mean net precipitation. The alluvial aquifer can significantly influence low flows when the bedrock productivity is limited. Although our synthetic catchment setup is simple, it is far more complex than the available analytical approaches or conceptual hydrological models. The direct application of the results to existing catchments requires nevertheless careful consideration of the local geological topographic and climatic conditions. This study provides quantitative insight into the complex interrelations between geology, topography and low-flow dynamics and challenges previous studies which neglect or oversimplify geological characteristics in the assessment of low flows.

A 3D geological model of a structurally complex Alpine region as a basis for interdisciplinary research.

Certain applications, such as understanding the influence of bedrock geology on hydrology in complex mountainous settings, demand 3D geological models that are detailed, high-resolution, accurate, and spatially-extensive. However, developing models with these characteristics remains challenging. Here, we present a dataset corresponding to a renowned tectonic entity in the Swiss Alps - the Nappe de Morcles - that does achieve these criteria. Locations of lithological interfaces and formation orientations were first extracted from existing sources. Then, using state-of-the-art algorithms, the interfaces were interpolated. Finally, an iterative process of evaluation and re-interpolation was undertaken. The geology was satisfactorily reproduced; modelled interfaces correspond well with the input data, and the estimated volumes seem plausible. Overall, 18 formations, including their associated secondary folds and selected faults, are represented at 10 m resolution. Numerous environmental investigations in the study area could benefit from the dataset; indeed, it is already informing integrated hydrological (snow/surface-water/groundwater) simulations. Our work demonstrates the potential that now exists to develop complex, high-quality geological models in support of contemporary Alpine research, augmenting traditional geological information in the process.

Channel representation in physically based models coupling groundwater and surface water: pitfalls and how to avoid them.

Recent models that couple three-dimensional subsurface flow with two-dimensional overland flow are valuable tools for quantifying complex groundwater/stream interactions and for evaluating their influence on watershed processes. For the modeler who is used to defining streams as a boundary condition, the representation of channels in integrated models raises a number of conceptual and technical issues. These models are far more sensitive to channel topography than conventional groundwater models. On all spatial scales, both the topography of a channel and its connection with the floodplain are important. For example, the geometry of river banks influences bank storage and overbank flooding; the slope of the river is a primary control on the behavior of a catchment; and at the finer scale bedform characteristics affect hyporheic exchange. Accurate data on streambed topography, however, are seldom available, and the spatial resolution of digital elevation models is typically too coarse in river environments, resulting in unrealistic or undulating streambeds. Modelers therefore perform some kind of manual yet often cumbersome correction to the available topography. In this context, the paper identifies some common pitfalls, and provides guidance to overcome these. Both aspects of topographic representation and mesh discretization are addressed. Additionally, two tutorials are provided to illustrate: (1) the interpolation of channel cross-sectional data and (2) the refinement of a mesh along a stream in areas of high topographic variability.

When can inverted water tables occur beneath streams?

Decline in regional water tables (RWT) can cause losing streams to disconnect from underlying aquifers. When this occurs, an inverted water table (IWT) will develop beneath the stream, and an unsaturated zone will be present between the IWT and the RWT. The IWT marks the base of the saturated zone beneath the stream. Although a few prior studies have suggested the likelihood of an IWT without a clogging layer, most of them have assumed that a low-permeability streambed is required to reduce infiltration from surface water to groundwater, and that the IWT only occurs at the bottom of the low-permeability layer. In this study, we use numerical simulations to show that the development of an IWT beneath an unclogged stream is theoretically possible under steady-state conditions. For a stream width of 1 m above a homogeneous and isotropic sand aquifer with a 47 m deep RWT (measured in an observation point 20 m away from the center of the stream), an IWT will occur provided that the stream depth is less than a critical value of 4.1 m. This critical stream depth is the maximum water depth in the stream to maintain the occurrence of an IWT. The critical stream depth decreases with stream width. For a stream width of 6 m, the critical stream depth is only 1 mm. Thus while theoretically possible, an IWT is unlikely to occur at steady state without a clogging layer, unless a stream is very narrow or shallow and the RWT is very deep.

An analysis of river bank slope and unsaturated flow effects on bank storage.

Recognizing the underlying mechanisms of bank storage and return flow is important for understanding streamflow hydrographs. Analytical models have been widely used to estimate the impacts of bank storage, but are often based on assumptions of conditions that are rarely found in the field, such as vertical river banks and saturated flow. Numerical simulations of bank storage and return flow in river-aquifer cross sections with vertical and sloping banks were undertaken using a fully-coupled, surface-subsurface flow model. Sloping river banks were found to increase the bank infiltration rates by 98% and storage volume by 40% for a bank slope of 3.4° from horizontal, and for a slope of 8.5°, delay bank return flow by more than four times compared with vertical river banks and saturated flow. The results suggested that conventional analytical approximations cannot adequately be used to quantify bank storage when bank slope is less than 60° from horizontal. Additionally, in the unconfined aquifers modeled, the analytical solutions did not accurately model bank storage and return flow even in rivers with vertical banks due to a violation of the dupuit assumption. Bank storage and return flow were also modeled for more realistic cross sections and river hydrograph from the Fitzroy River, Western Australia, to indicate the importance of accurately modeling sloping river banks at a field scale. Following a single wet season flood event of 12 m, results showed that it may take over 3.5 years for 50% of the bank storage volume to return to the river.

Disconnected surface water and groundwater: from theory to practice.

When describing the hydraulic relationship between rivers and aquifers, the term disconnected is frequently misunderstood or used in an incorrect way. The problem is compounded by the fact that there is no definitive literature on the topic of disconnected surface water and groundwater. We aim at closing this gap and begin the discussion with a short introduction to the historical background of the terminology. Even though a conceptual illustration of a disconnected system was published by Meinzer (1923), it is only within the last few years that the underlying physics of the disconnection process has been described. The importance of disconnected systems, however, is not widely appreciated. Although rarely explicitly stated, many approaches for predicting the impacts of groundwater development on surface water resources assume full connection. Furthermore, management policies often suggest that surface water and groundwater should only be managed jointly if they are connected. However, although lowering the water table beneath a disconnected section of a river will not change the infiltration rate at that point, it can increase the length of stream that is disconnected. Because knowing the state of connection is of fundamental importance for sustainable water management, robust field methods that allow the identification of the state of connection are required. Currently, disconnection is identified by showing that the infiltration rate from a stream to an underlying aquifer is independent of the water table position or by identifying an unsaturated zone under the stream. More field studies are required to develop better methods for the identification of disconnection and to quantify the implications of heterogeneity and clogging processes in the streambed on disconnection.

Modeling surface water-groundwater interaction with MODFLOW: some considerations.

The accuracy with which MODFLOW simulates surface water-groundwater interaction is examined for connected and disconnected losing streams. We compare the effect of different vertical and horizontal discretization within MODFLOW and also compare MODFLOW simulations with those produced by HydroGeoSphere. HydroGeoSphere is able to simulate both saturated and unsaturated flow, as well as surface water, groundwater and the full coupling between them in a physical way, and so is used as a reference code to quantify the influence of some of the simplifying assumptions of MODFLOW. In particular, we show that (1) the inability to simulate negative pressures beneath disconnected streams in MODFLOW results in an underestimation of the infiltration flux; (2) a river in MODFLOW is either fully connected or fully disconnected, while in reality transitional stages between the two flow regimes exist; (3) limitations in the horizontal discretization of the river can cause a mismatch between river width and cell width, resulting in an error in the water table position under the river; and (4) because coarse vertical discretization of the aquifer is often used to avoid the drying out of cells, this may result in an error in simulating the height of the groundwater mound. Conditions under which these errors are significant are investigated.