Explicit simulation of environmental gas tracers with integrated surface and subsurface hydrological models

Abstract

Environmental gas tracers allow inferring groundwater travel times and mixing ratios. Their concentrations are commonly interpreted with simplified and indirect approaches that are conceptually at odds with the high degree of complexity found in natural systems. However, the information content of the tracers can potentially be fully explored through the explicit simulation of an advection-dispersion transport equation, for example using integrated surface-subsurface hydrological models (ISSHMs). These integrated models can be used to explicitly simulate environmental tracers in complex environments. ISSHMs are usually variably saturated flow models. However, these models do not explicitly simulate gas partitioning with the aqueous phase, restricting explicit simulation of gas tracers to fully saturated conditions or to tracers with very low solubilities. We propose a mathematical formulation for the production of environmental gas tracers that are emanated in the subsurface. The production is scaled according to gas/water partitioning and water saturation, which is already computed by the model. Therefore, ISSHMs can now be used to their full potential to explicitly simulate tracer concentrations under variably saturated and dynamic conditions. The new formulation has been successfully verified against reference simulations provided with a multi-phase flow and transport model. In addition, explicit simulation of 222Rn and 37Ar groundwater concentrations in a synthetic alluvial river-groundwater system was demonstrated, for the first time, with an ISSHM.