Considering heterogeneity in porous media pore size and connectivity is essential to predicting reactive solute transport across interfaces. However, exchange with less‐mobile porosity is rarely considered in surface water/groundwater recharge studies. Previous research indicates that a combination of pore‐fluid sampling and geoelectrical measurements can be used to quantify less‐mobile porosity exchange dynamics using the time‐varying relation between fluid and bulk electrical conductivity. For this study, we use macro‐scale (10 s of cm) advection–dispersion solute transport models linked with electrical conduction in COMSOL Multiphysics to explore less‐mobile porosity dynamics in two different types of observed sediment water interface porous media. Modeled sediment textures contrast from strongly layered streambed deposits to poorly sorted lakebed sands and cobbles. During simulated ionic tracer perturbations, a lag between fluid and bulk electrical conductivity, and the resultant hysteresis, is observed for all simulations indicating differential loading of pore spaces with tracer. Less‐mobile exchange parameters are determined graphically from these tracer time series data without the need for inverse numerical model simulation. In both sediment types, effective less‐mobile porosity exchange parameters are variable in response to changes in flow direction and fluid flux. These observed flow‐dependent effects directly impact local less‐mobile residence times and associated contact time for biogeochemical reaction. The simulations indicate that for the sediment textures explored here, less‐mobile porosity exchange is dominated by variable rates of advection through the domain, rather than diffusion of solute, for typical low‐to‐moderate rate (approximately 3–40 cm/day) hyporheic fluid fluxes. Overall, our model‐based results show that less‐mobile porosity may be expected in a range of natural hyporheic sediments and that changes in flowpath orientation and magnitude will impact less‐mobile exchange parameters. These temporal dynamics can be assessed with the geoelectrical experimental tracer method applied at laboratory and field scales.
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- Award ID(s):
- 2142686
- NSF-PAR ID:
- 10371263
- Publisher / Repository:
- DOI PREFIX: 10.1029
- Date Published:
- Journal Name:
- Water Resources Research
- Volume:
- 58
- Issue:
- 7
- ISSN:
- 0043-1397
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
more » « less
Abstract -
more » « less
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
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.
-
null (Ed.)Minimum free energy prediction of RNA secondary structures is based on the Nearest Neighbor Thermodynamics Model. While such predictions are typically good, the accuracy can vary widely even for short sequences, and the branching thermodynamics are an important factor in this variance. Recently, the simplest model for multiloop energetics—a linear function of the number of branches and unpaired nucleotides—was found to be the best. Subsequently, a parametric analysis demonstrated that per family accuracy can be improved by changing the weightings in this linear function. However, the extent of improvement was not known due to the ad hoc method used to find the new parameters. Here we develop a branch-and-bound algorithm that finds the set of optimal parameters with the highest average accuracy for a given set of sequences. Our analysis shows that the previous ad hoc parameters are nearly optimal for tRNA and 5S rRNA sequences on both training and testing sets. Moreover, cross-family improvement is possible but more difficult because competing parameter regions favor different families. The results also indicate that restricting the unpaired nucleotide penalty to small values is warranted. This reduction makes analyzing longer sequences using the present techniques more feasible.more » « less
-
more » « less
Abstract Hyporheic exchange has the potential to significantly influence river temperatures in regions of continuous permafrost under low‐flow conditions given the strong thermal gradients that exist in river bed sediments. However, there is limited understanding of the impacts of hyporheic exchange on Arctic river temperatures. To address this knowledge gap, heat fluxes associated with hyporheic exchange were estimated in a fourth‐order Arctic river using field observations coupled with a river temperature model that accounts for hyporheic exchange influences. Temperature time series and tracer study solute breakthrough curves were measured in the main channel and river bed at multiple locations and depths to characterize hyporheic exchange and provide parameter bounds for model calibration. Model results for low‐flow periods from 3 years indicated that hyporheic exchange contributed up to 27% of the total river energy balance, reduced the main channel diel temperature range by up to 1.7 °C, and reduced mean daily temperatures by up to 0.21 °C over a 13.1‐km study reach. These influences are due to main channel heat loss during the day and gain at night via hyporheic exchange and heat loss from the hyporheic zone to the ground below via conduction. Main channel temperatures were found to be sensitive to simulated changes in ground temperatures due to changes in hyporheic exchange heat flux and deeper ground conduction. These results suggest that the moderating influence of hyporheic exchange could be reduced if ground temperatures warm in response to projected increases in permafrost thaw below rivers.
