skip to main content


Title: Exploring Tracer Information and Model Framework Trade‐Offs to Improve Estimation of Stream Transient Storage Processes
Abstract

Novel observation techniques (e.g., smart tracers) for characterizing coupled hydrological and biogeochemical processes are improving understanding of stream network transport and transformation dynamics. In turn, these observations are thought to enable increasingly sophisticated representations within transient storage models (TSMs). However, TSM parameter estimation is prone to issues with insensitivity and equifinality, which grow as parameters are added to model formulations. Currently, it is unclear whether (or not) observations from different tracers may lead to greater process inference and reduced parameter uncertainty in the context of TSM. Herein, we aim to unravel the role of in‐stream processes alongside metabolically active (MATS) and inactive storage zones (MITS) using variable TSM formulations. Models with one (1SZ) and two storage zones (2SZ) and with and without reactivity were applied to simulate conservative and smart tracer observations obtained experimentally for two reaches with differing morphologies. As we show, smart tracers are unsurprisingly superior to conservative tracers when it comes to partitioning MITS and MATS. However, when transient storage is lumped within a 1SZ formulation, little improvement in parameter uncertainty is gained by using a smart tracer, suggesting the addition of observations should scale with model complexity. Importantly, our work identifies several inconsistencies and open questions related to reconciling time scales of tracer observation with conceptual processes (parameters) estimated within TSM. Approaching TSM with multiple models and tracer observations may be key to gaining improved insight into transient storage simulation as well as advancing feedback loops between models and observations within hydrologic science.

 
more » « less
Award ID(s):
1652293
NSF-PAR ID:
10453258
Author(s) / Creator(s):
 ;  ;  ;  ;  ;  ;  ;  ;  ;  ;  ;  ;  ;  
Publisher / Repository:
DOI PREFIX: 10.1029
Date Published:
Journal Name:
Water Resources Research
Volume:
55
Issue:
4
ISSN:
0043-1397
Page Range / eLocation ID:
p. 3481-3501
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    Hyporheic exchange in streams is critical to ecosystem functions such as nutrient cycling along river corridors, especially for slowly moving or small stream systems. The transient storage model (TSM) has been widely used for modeling of hyporheic exchange. TSM calibration, for hyporheic exchange, is typically used to estimate four parameters, including the mass exchange rate coefficient, the dispersion coefficient, stream cross‐sectional area, and hyporheic zone cross‐sectional area. Prior studies have raised concerns regarding the non‐uniqueness of the inverse problem for the TSM, that is, the occurrence of different parameter vectors resulting in TSM solution that reproduces the observed in‐stream tracer break through curve (BTC) with the same error. This leads to practical non‐identifiability in determining the unknown parameter vector values even when global‐optimal values exist, and the parameter optimization becomes practically non‐unique. To address this problem, we applied the simulated annealing method to calibrate the TSM to BTCs, because it is less susceptible to local minima‐induced non‐identifiability. A hypothetical (or synthetic) tracer test data set with known parameters was developed to demonstrate the capability of the simulated annealing method to find the global minimum parameter vector, and it identified the “hypothetically‐true” global minimum parameter vector even with input data that were modified with up to 10% noise without increasing the number of iterations required for convergence. The simulated annealing TSM was then calibrated using two in‐stream tracer tests conducted in East Fork Poplar Creek, Tennessee. Simulated annealing was determined to be appropriate for quantifying the TSM parameter vector because of its search capability for the global minimum parameter vector.

     
    more » « less
  2. Abstract. Resazurin (Raz) and its reaction product resorufin (Rru) have increasingly been used as reactive tracers to quantify metabolic activity and hyporheic exchange in streams. Previous work has indicated that these compounds undergo sorption in stream sediments. We present laboratory experiments on Raz and Rru transport, sorption, and transformation, consisting of 4 column and 72 batch tests using 2 sediments with different physicochemical properties under neutral (pH = 7) and alkaline (pH = 9) conditions. The study aimed at identifying the key processes of reactive transport of Raz and Rru in streambed sediments and the experimental setup best suited for their determination. Data from column experiments were simulated by a travel-time-based model accounting for physical transport, equilibrium and kinetic sorption, and three first-order reactions. We derived the travel-time distributions directly from the breakthrough curve (BTC) of the conservative tracer, fluorescein, rather than from fitting an advective-dispersive transport model, and inferred from those distributions the transfer functions of Raz and Rru, which provided conclusive approximations of the measured BTCs. The most likely reactive transport parameters and their uncertainty were determined by a Markov chain–Monte Carlo approach. Sorption isotherms of both compounds were obtained from batch experiments. We found that kinetic sorption dominates sorption of both Raz and Rru, with characteristic timescales of sorption in the order of 12 to 298 min. Linear sorption models for both Raz and Rru appeared adequate for concentrations that are typically applied in field tracer tests. The proposed two-site sorption model helps to interpret transient tracer tests using the Raz–Rru system.

     
    more » « less
  3. 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
  4. Abstract

    Land use within a watershed impacts stream channel morphology and hydrology and, therefore, in‐stream solute transport processes and associated transient storage mechanisms. This study evaluated transport processes in two contrasting stream sites where channel morphology was influenced by the surrounding land use, land cover, climate and geologic controls: Como Creek, CO, a relatively undisturbed, high gradient, forested stream with a gravel bed and complex channel morphology, and Clear Creek, IA, an incised, low‐gradient stream with low‐permeability substrate draining an agricultural landscape. We performed conservative stream tracer injections at these sites to address the following questions: (1) How does solute transport vary between streams with differing morphologies? and (2) How does solute transport at each stream site change as a function of discharge? We analysed in‐stream tracer time series data and compared results quantifying solute attenuation in surface and subsurface transient storage zones. Significant trends were observed in these metrics with varying discharge conditions at the forested site but not at the agricultural site. There was a broad range of transport mechanisms and evidence of substantial exchange with both surface and hyporheic transient storage in the relatively undisturbed, forested stream. Changing discharge conditions activated or deactivated different solute transport mechanisms in the forested site and greatly impacted advective travel time. Conversely in the simplified agricultural stream, there was a narrow range of solute transport behaviour across flows and predominantly surface transient storage at all measured discharge conditions. These results demonstrate how channel simplification inhibits available solute transport mechanisms across varying discharge conditions.

     
    more » « less
  5. null (Ed.)
    High concentrations of trace metal(loid)s exported from abandoned mine wastes and acid rock drainage pose a risk to the health of aquatic ecosystems. To determine if and when the hyporheic zone mediates metal(loid) export, we investigated the relationship between streamflow, groundwater–stream connectivity, and subsurface metal(loid) concentrations in two ~1-km stream reaches within the Bonita Peak Mining District, a US Environmental Protection Agency Superfund site located near Silverton, Colorado, USA. The hyporheic zones of reaches in two streams—Mineral Creek and Cement Creek—were characterized using a combination of salt-tracer injection tests, transient-storage modeling, and geochemical sampling of the shallow streambed (<0.7 m). Based on these data, we present two conceptual models for subsurface metal(loid) behavior in the hyporheic zones, including (1) well-connected systems characterized by strong hyporheic mixing of infiltrating stream water and upwelling groundwater and (2) poorly connected systems delineated by physical barriers that limit hyporheic mixing. The comparatively large hyporheic zone and high hydraulic conductivities of Mineral Creek created a connected stream–groundwater system, where mixing of oxygen-rich stream water and metal-rich groundwater facilitated the precipitation of metal colloids in the shallow subsurface. In Cement Creek, the precipitation of iron oxides at depth (~0.4 m) created a low-hydraulic-conductivity barrier between surface water and groundwater. Cemented iron oxides were an important regulator of metal(loid) concentrations in this poorly connected stream–groundwater system due to the formation of strong redox gradients induced by a relatively small hyporheic zone and high fluid residence times. A comparison of conceptual models to stream concentration–discharge relationships exhibited a clear link between geochemical processes occurring within the hyporheic zone of the well-connected system and export of particulate Al, Cu, Fe, and Mn, while the poorly connected system did not have a notable influence on metal concentration–discharge trends. Mineral Creek is an example of a hyporheic system that serves as a natural dissolved metal(loid) sink, whereas poorly connected systems such as Cement Creek may require a combination of subsurface remediation of sediments and mitigation of upstream, iron-rich mine drainages to reduce metal export. 
    more » « less