An official website of the United States government
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
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
- PAR ID:
- 10453258
- 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
-
-
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
-
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
-
Abstract Pollutant transport in discrete fracture networks (DFNs) exhibits complex dynamics that challenge reliable model predictions, even with detailed fracture data. To address this issue, this study derives an upscaled integral‐differential equation to predict transient anomalous diffusion in two‐dimensional (2D) DFNs. The model includes both transmissive and dead‐end fractures (DEFs), where stagnant water zones in DEFs cause non‐uniform flow and transient sub‐diffusive transport, as shown by both literature and DFN flow and transport simulations using COMSOL. The upscaled model's main parameters are quantitatively linked to fracture properties, especially the probability density function of DEF lengths. Numerical experiments show the model's accuracy in predicting the full‐term evolution of conservative tracers in 2D DFNs with power‐law distributed fracture lengths and two orientation sets. Field applications indicate that while model parameters for transient sub‐diffusion can be predicted from observed DFN distributions, predicting parameters controlling solute displacement in transmissive fractures requires additional field work, such as tracer tests. Parameter sensitivity analysis further correlates late‐time solute transport dynamics with fracture properties, such as fracture density and average length. Potential extensions of the upscaled model are also discussed. This study, therefore, proves that transient anomalous transport in 2D DFNs with DEFs can be at least partially predicted, offering an initial step toward improving model predictions for pollutant transport in real‐world fractured aquifer systems.more » « less
-
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
