skip to main content


Title: Utilizing Environmental Tracers to Reduce Groundwater Flow and Transport Model Parameter Uncertainties
Abstract

Non‐uniqueness in groundwater model calibration is a primary source of uncertainty in groundwater flow and transport predictions. In this study, we investigate the ability of environmental tracer information to constrain groundwater model parameters. We utilize a pilot point calibration procedure conditioned to subsets of observed data including: liquid pressures, tritium (3H), chlorofluorocarbon‐12 (CFC‐12), and sulfur hexafluoride (SF6) concentrations; and groundwater apparent ages inferred from these environmental tracers, to quantify uncertainties in the heterogeneous permeability fields and infiltration rates of a steady‐state 2‐D synthetic aquifer and a transient 3‐D model of a field site located near Riverton, Wyoming (USA). To identify the relative data worth of each observation data type, the post‐calibration uncertainties of the optimal parameters for a given observation subset are compared to that from the full observation data set. Our results suggest that the calibration‐constrained permeability field uncertainties are largest when liquid pressures are used as the sole calibration data set. We find significant reduction in permeability uncertainty and increased predictive accuracy when the environmental tracer concentrations, rather than apparent groundwater ages, are used as calibration targets in the synthetic model. Calibration of the Riverton field site model using environmental tracer concentrations directly produces infiltration rate estimates with the lowest uncertainties, however; permeability field uncertainties remain similar between the environmental tracer concentration and apparent groundwater age calibration scenarios. This work provides insight on the data worth of environmental tracer information to calibrate groundwater models and highlights potential benefits of directly assimilating environmental tracer concentrations into model parameter estimation procedures.

 
more » « less
Award ID(s):
1633831
NSF-PAR ID:
10446509
Author(s) / Creator(s):
 ;  ;  
Publisher / Repository:
DOI PREFIX: 10.1029
Date Published:
Journal Name:
Water Resources Research
Volume:
57
Issue:
7
ISSN:
0043-1397
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    Submarine groundwater discharge (SGD) is an important driver of coastal biogeochemical budgets worldwide. Radon (222Rn) has been widely used as a natural geochemical tracer to quantify SGD, but field measurements are time consuming and costly. Here, we use deep learning to predict coastal seawater radon in SGD‐impacted regions. We hypothesize that deep learning could resolve radon trends and enable preliminary insights with limited field observations of groundwater tracers. Two deep learning models were trained on global coastal seawater radon observations (n = 39,238) with widely available inputs (e.g., salinity, temperature, water depth). The first model used a one‐dimensional convolutional neural network (1D‐CNN‐RNN) framework for site‐specific gap filling and producing short‐term future predictions. A second model applied a fully connected neural network (FCNN) framework to predict radon across geographically and hydrologically diverse settings. Both models can predict observed radon concentrations withr2 > 0.76. Specifically, the FCNN model offers a compelling development because synthetic radon tracer data sets can be obtained using only basic water quality and meteorological parameters. This opens opportunities to attain radon data from regions with large data gaps, such as the Global South and other remote locations, allowing for insights that can be used to predict SGD and plan field experiments. Overall, we demonstrate how field‐based measurements combined with big‐data approaches such as deep learning can be utilized to assess radon and potentially SGD beyond local scales.

     
    more » « less
  2. Abstract

    Time‐lapse electrical resistivity tomography (ERT) data are increasingly used to inform the hydrologic dynamics of mountainous environments at the hillslope scale. Despite their popularity and recent advancements in hydrogeophysical inversion methods, few studies have shown how time‐lapse ERT data can be used to determine hydraulic parameters of subsurface water flow models. This study uses synthetic and field‐collected, hillslope‐scale, time‐lapse ERT data to determine subsurface hydraulic properties of a two‐layer, physics‐based, 2‐D vertical flow model with predefined layer and boundary locations. Uncoupled and coupled hydrogeophysical inversion methods are combined with a fine‐earth fraction optimization scheme to reduce the number of parameters needing calibration and interpret the influence of the hydraulic parameters on the hydrologic model predictions. Inversions of synthetic ERT data recover the prescribed fine‐earth fraction bulk density to within 0.1 g cm−3. Field‐collected ERT data from a mountain hillslope result in hydrologic model dynamics that are consistent with previous studies and measured water content data but struggle to capture measured groundwater levels. The uncoupled hydrogeophysical inversion method is more sensitive to changes in hydraulic parameter values of the lower hydrologic model layer than the coupled hydrogeophysical inversion method. Time series of minimum objective function value simulations indicate that periodically collected ERT data may recover hydraulic parameters to a similar level of uncertainty as daily ERT data. Using simple hydrologic model domains within hydrogeophysical inversions shows promise for providing reasonable hydrologic predictions while maintaining relatively simple calibration schemes and should be explored further in future studies.

     
    more » « less
  3. null (Ed.)
    Abstract. The acceleration of urbanization requires sustainable, adaptive management strategies for land and water use in cities. Although the effects of buildings and sealed surfaces on urban runoff generation and local climate are well known, much less is known about the role of water partitioning in urban green spaces. In particular, little is quantitatively known about how different vegetation types of urban green spaces (lawns, parks, woodland, etc.) regulate partitioning of precipitation into evaporation, transpiration and groundwater recharge and how this partitioning is affected by sealed surfaces. Here, we integrated field observations with advanced, isotope-based ecohydrological modelling at a plot-scale site in Berlin, Germany. Soil moisture and sap flow, together with stable isotopes in precipitation, soil water and groundwater recharge, were measured over the course of one growing season under three generic types of urban green space: trees, shrub and grass. Additionally, an eddy flux tower at the site continuously collected hydroclimate data. These data have been used as input and for calibration of the process-based ecohydrological model EcH2O-iso. The model tracks stable isotope ratios and water ages in various stores (e.g. soils and groundwater) and fluxes (evaporation, transpiration and recharge). Green water fluxes in evapotranspiration increased in the order shrub (381±1mm) more » « less
  4. We introduce EcH2O-iso, a new development of the physically-based, fully-distributed ecohydrological model EcH2O where the tracking of water isotopic tracers (2H and 18O) and age has been incorporated. EcH2O-iso is evaluated at a montane, low-energy experimental catchment in eastern Scotland using 16 independent isotope time series from various landscape positions and compartments; encompassing soil water, groundwater, stream water, and plant xylem. We find a good model-observation match in most cases, despite having only calibrated the model using hydrometric data and energy fluxes. These results provide further validation of the physical basis of the model for successfully capturing catchment hydrological functioning, both in terms of the celerity in energy propagation (e.g. runoff generation under prevailing hydraulic gradients) and flow velocities of water molecules (e.g., in consistent tracer concentrations at given locations and times). We also show that the spatially-distributed formulation of EcH2O-iso provides a powerful tool for quantitatively linking water stores and fluxes with spatio-temporal patterns of isotopes ratios and water ages. Finally, our study highlights some model development and benchmarking needs, refined using isotope-based calibration, for hypothesis testing and improved simulations of catchment dynamics that is transferable beyond the catchment landscape studied here. 
    more » « less
  5. Abstract

    Uncertainty quantification of groundwater (GW) aquifer parameters is critical for efficient management and sustainable extraction of GW resources. These uncertainties are introduced by the data, model, and prior information on the parameters. Here, we develop a Bayesian inversion framework that uses Interferometric Synthetic Aperture Radar (InSAR) surface deformation data to infer the laterally heterogeneous permeability of a transient linear poroelastic model of a confined GW aquifer. The Bayesian solution of this inverse problem takes the form of a posterior probability density of the permeability. Exploring this posterior using classical Markov chain Monte Carlo (MCMC) methods is computationally prohibitive due to the large dimension of the discretized permeability field and the expense of solving the poroelastic forward problem. However, in many partial differential equation (PDE)‐based Bayesian inversion problems, the data are only informative in a few directions in parameter space. For the poroelasticity problem, we prove this property theoretically for a one‐dimensional problem and demonstrate it numerically for a three‐dimensional aquifer model. We design a generalized preconditioned Crank‐Nicolson (gpCN) MCMC method that exploits this intrinsic low dimensionality by using a low‐rank‐based Laplace approximation of the posterior as a proposal, which we build scalably. The feasibility of our approach is demonstrated through a real GW aquifer test in Nevada. The inherently two‐dimensional nature of InSAR surface deformation data informs a sufficient number of modes of the permeability field to allow detection of major structures within the aquifer, significantly reducing the uncertainty in the pressure and the displacement quantities of interest.

     
    more » « less