Evidence of distinct contaminant transport patterns in rivers using tracer tests and a multiple domain retention model

Solute transport in rivers is controlled by surface hydrodynamics and by mass exchanges with distinct retention zones. Surface and hyporheic retention processes can be accounted for separately in solute transport models with multiple storage compartments. In the simplest two component model, short term storage can be associated to in-channel transient retention, e.g. produced by riparian vegetation or surface dead zones, and the long-term storage can be associated to hyporheic exchange. The STIR (Solute Transport In Rivers) multiple domain transport model is applied here to tracer test data from three very different Mediterranean streams with distinctive characteristics in terms of flow discharge, vegetation and substrate material. The model is used with an exponential residence time distribution (RTD) to represent surface storage processes and two distinct modeling closures are tested to simulate hyporheic retention: a second exponential RTD and a power-law distribution approximating a known solution for bedform-induced hyporheic exchange. Each stream shows distinct retention patterns characterized by different timescales of the storage time distribution. Both modeling closures lead to very good approximations of the observed breakthrough curves in the two rivers with permeable bed exposed to the flow, where hyporheic flows are expected to occur. In the one case where the occurrence of hyporheic flows is inhibited by bottom vegetation, only the two exponential RTD model is acceptable and the time scales of the two components are of the same magnitude. The significant finding of this work is the recognition of a strong signature of the river properties on tracer data and the evidence of the ability of multiple-component models to describe individual stream responses. This evidence may open a new perspective in river contamination studies, where rivers could possibly be classified based on their ability to trap and release pollutants.     

How Does Subsurface Characterization Affect Simulations of Hyporheic Exchange?

We investigated the role of increasingly well‐constrained geologic structures in the subsurface (i.e., subsurface architecture) in predicting streambed flux and hyporheic residence time distribution (RTD) for a headwater stream. Five subsurface realizations with increasingly resolved lithological boundaries were simulated in which model geometries were based on increasing information about flow and transport using soil and geologic maps, surface observations, probing to depth to refusal, seismic refraction, electrical resistivity (ER) imaging of subsurface architecture, and time‐lapse ER imaging during a solute tracer study. Particle tracking was used to generate RTDs for each model run. We demonstrate how improved characterization of complex lithological boundaries and calibration of porosity and hydraulic conductivity affect model prediction of hyporheic flow and transport. Models using hydraulic conductivity calibrated using transient ER data yield estimates of streambed flux that are three orders of magnitude larger than uncalibrated models using estimated values for hydraulic conductivity based on values published for nearby hillslopes (10^?4 vs. 10^?7 m²/s, respectively). Median residence times for uncalibrated and calibrated models are 10³ and 10⁰ h, respectively. Increasingly well‐resolved subsurface architectures yield wider hyporheic RTDs, indicative of more complex hyporheic flowpath networks and potentially important to biogeochemical cycling. The use of ER imaging to monitor solute tracers informs subsurface structure not apparent from other techniques, and helps to define transport properties of the subsurface (i.e., hydraulic conductivity). Results of this study demonstrate the value of geophysical measurements to more realistically simulate flow and transport along hyporheic flowpaths.     

Hydraulic analysis of cell-network treatment wetlands

When individual cells of a multiple-cell treatment wetland are hydraulically connected, the wetland has a cell-network structure. The hydraulic performance of treatment wetlands is often characterized using tracer residence time distributions (RTDs) measured between the wetland inlet and outlet, such that the wetland is considered as a single hydraulic unit, regardless of the extent of networking between individual internal cells. This work extends the single hydraulic unit approach to enable the specification of moments and RTD parameters for individual cells, or clusters of cells, within the cell-network based on inert tracer tests with injection only at the network inlet. Hydraulic performance is quantified in terms of hydraulic efficiency and travel time dimensionless variance using both the method of moments and RTD modeling. Cell-network analysis was applied to a case study from the Orlando Easterly Wetland (OEW), demonstrating the improvement in hydraulic performance of individual wetland cells following wetland restoration activities. Furthermore, cell-network analysis indicated that the location of water quality sampling station locations within the cell network can significantly affect the accuracy of pollutant removal effectiveness estimation when the individual sample station RTD does not represent the hydraulic unit RTD. At the OEW, it was determined that historical nutrient removal effectiveness estimation may be underestimated for one area and overestimated for another, and recommendations were provided for sample station locations to minimize future performance estimation errors.     

Seasonal variation in cascade‐driven hyporheic exchange,northern Honduras

A characterization of hyporheic exchange for dry and wet season baseflow, as well as partially dewatered discharge, was done in Prieta Creek, a first‐order cascade in northern Honduras. The cascade had discharges from 1 to 15 l s^?1, had average slopes of 12%, pool spacing of 3 m, and shallow substrate of sand and gravel. Tracer tests were conducted in a 15‐m sub‐reach, a length considered to be adequate for the experiment based on the DaI test, a ratio of exchange and transport processes. In the three tests, between 9 and 18% of tracer was not recovered, possibly due to entrainment in flowpaths passing beneath the downstream monitoring location. Tracer data were analysed by the one‐dimensional transport with inflow and storage (OTIS) transient storage model (TSM) to derive standard exchange parameters, and by the solute transport in rivers (STIR) model to examine hyporheic residence time distributions (RTDs). The best fit of the observed tracer breakthrough curves was obtained by using the STIR model with a combination of two exponential RTDs to represent hyporheic retention. With increasing discharge, the OTIS model predicted increasing storage exchange fluxes and exchange coefficients and decreasing storage zone areas and transient storage times, which are trends supported by riparian and streambed piezometric head data. Riparian water levels rose during the transition from the dry to wet season, which could constrict the hyporheic storage zone. Thirteen of the 19 streambed piezometers recorded seasonal changes in hydraulic gradients and flux direction, with fewer yet stronger upwelling zones during higher discharges. The MODFLOW model missed the observed seasonal changes, possibly due to subtle changes in the seasonal change in water surface profiles. We conclude that partially dewatered dry season exchange, compared to wet season exchange, was initiated and terminated with smaller pressure gradients and, in different streambed locations, was smaller in volume, had longer residence times, and may connect with deeper and longer flow paths. Copyright © 2010 John Wiley & Sons, Ltd.     

Continuous time random walks for non-local radial solute transport

This study formulates and analyzes continuous time random walk (CTRW) models in radial flow geometries for the quantification of non-local solute transport induced by heterogeneous flow distributions and by mobile–immobile mass transfer processes. To this end we derive a general CTRW framework in radial coordinates starting from the random walk equations for radial particle positions and times. The particle density, or solute concentration is governed by a non-local radial advection–dispersion equation (ADE). Unlike in CTRWs for uniform flow scenarios, particle transition times here depend on the radial particle position, which renders the CTRW non-stationary. As a consequence, the memory kernel characterizing the non-local ADE, is radially dependent. Based on this general formulation, we derive radial CTRW implementations that (i) emulate non-local radial transport due to heterogeneous advection, (ii) model multirate mass transfer (MRMT) between mobile and immobile continua, and (iii) quantify both heterogeneous advection in a mobile region and mass transfer between mobile and immobile regions. The expected solute breakthrough behavior is studied using numerical random walk particle tracking simulations. This behavior is analyzed by explicit analytical expressions for the asymptotic solute breakthrough curves. We observe clear power-law tails of the solute breakthrough for broad (power-law) distributions of particle transit times (heterogeneous advection) and particle trapping times (MRMT model). The combined model displays two distinct time regimes. An intermediate regime, in which the solute breakthrough is dominated by the particle transit times in the mobile zones, and a late time regime that is governed by the distribution of particle trapping times in immobile zones. These radial CTRW formulations allow for the identification of heterogeneous advection and mobile-immobile processes as drivers of anomalous transport, under conditions relevant for field tracer tests.     

Solute transport routing in a small stream

Abstract

The impact of pollution incidents on rivers and streams may be predicted using mathematical models of solute transport. Practical applications require an analytical or numerical solution to a governing solute mass balance equation together with appropriate values of relevant transport coefficients under the flow conditions of interest. This paper considers two such models, namely those proposed by Fischer and by Singh and Beck, and compares their performances using tracer data from a small stream in Edinburgh, UK. In calibrating the models, information on the magnitudes and the flow rate dependencies of the velocity and the dispersion coefficients was generated. The dispersion coefficient in the stream ranged between 0.1 and 0.9 m²/s for a flow rate range of 13–437 L/s. During calibration it was found that the Singh and Beck model fitted the tracer data a little better than the Fischer model in the majority of cases. In a validation exercise, however, both models gave similarly good predictions of solute transport at three different flow rates.     

An efficient space-fractional dispersion approximation for stream solute transport modeling

The success of transient storage (TS) modeling for natural streams depends, in part, on the ability to describe the dispersion process accurately. Evidence based on stream tracer data shows that solute transport processes often do not follow the classical second-order dispersion model (e.g., early breakthrough and faster than Fickian travel times were observed). While models based on space-fractional dispersion are a promising alternative, different definitions of fractional derivatives exist in the literature. Unlike integer-order derivatives, fractional derivatives represent convolutions of concentration with long-range spatial correlation and numerical approximations can produce dense matrices. Therefore issues of both accuracy and computational efficiency need to be examined to successfully identify model parameters for natural streams. In this paper, we first compare the performance of several numerical approaches for solving the space-fractional dispersion equation. We examine three different numerical approaches to approximate the space-fractional derivatives including: (a) a fully-implicit scheme based on the shifted Grünwald–Letnikov (GL) approximation (b) a three-point implicit representation based on the GL formula and (c) a three-point implicit scheme based on mass conservation and the Caputo definition of the fractional derivative. We then use an operator-splitting technique to evaluate a TS model based on space-fractional dispersion (the FSTS model) and test the model against analytical solutions and stream tracer data. A sequence acceleration method (Richardson extrapolation) significantly improves the performance of all schemes examined. Results indicate that the fully-implicit GL method with Richardson extrapolation produces the most accurate solutions while the three-point implicit GL scheme has a stringent time-step restriction to produce acceptable solutions. The three-point implicit scheme based on the Caputo derivative produces accurate solutions in a fraction of the time taken by the fully-implicit GL method and represents the best trade-off between accuracy and computational efficiency for practical applications. The scheme is suitable for parameter estimation and is used to successfully describe tracer data in a natural stream.     

Estimating longitudinal dispersion in rivers using Acoustic Doppler Current Profilers

The longitudinal dispersion coefficient (D) is an important parameter needed to describe the transport of solutes in rivers and streams. The dispersion coefficient is generally estimated from tracer studies but the method can be expensive and time consuming, especially for large rivers. A number of empirical relations are available to estimate the dispersion coefficient; however, these relations are known to produce estimates within an order of magnitude of the tracer value. The focus of this paper is on using the shear-flow dispersion theory to directly estimate the dispersion coefficient from velocity measurements obtained using an Acoustic Doppler Current Profiler (ADCP). Using tracer and hydrodynamic data collected within the same river reaches, we examined conditions under which the ADCP and tracer methods produced similar results. Since dead zones / transient storage (TS) are known to influence the dispersion coefficient, we assessed the relative importance of dead zones in different stream reaches using two tracer-based approaches: (1) TS modeling which explicitly accounts for dead zones and (2) the advection–dispersion equation (ADE) which does not have separate terms for dead zones. Dispersion coefficients based on the ADE tend to be relatively high as they describe some of the effects of dead zones as well. Results based on the ADCP method were found to be in good agreement with the ADE estimates indicating that storage zones play an important role in the estimated dispersion coefficients, especially at high flows. For the river sites examined in this paper, the tracer estimates of dispersion were close to the median values of the ADCP estimates obtained from multiple datasets within a reach. The ADCP method appears to be an excellent alternative to the traditional tracer-based method if care is taken to avoid spurious data and multiple datasets are used to compute a distance-weighted average or other appropriate measure that represents reach-averaged conditions.     

Characterizing solute transport with transient storage across a range of flow rates: The evidence of repeated tracer experiments in Austrian and Italian streams

Solute transport in rivers and streams with hyporheic zone exchange and/or in-stream storage is typically affected by the prevailing flow rate. The research reported here focuses on stream tracer experiments repeated many times along the same Austrian (Mödlingbach) and Italian (Torrente Lura) channel reaches to characterize parameter dependency on flow rate. Both groups of data sets showed an increase of storage zone area and main stream area with discharge. In either case, a strong negative correlation was obtained between storage zone residence time and flow rate. From the Mödlingbach data, no clear relationship with Q emerged for the dispersion coefficient and the dead zone ratio, whereas Torrente Lura showed a clear positive correlation of the dispersion coefficient with the flow rate and a slightly negative Q-dependency for the dead zone ratio. Mödlingbach and Torrente Lura results are presented against the background of other repeat experiments reported in literature.     

Bias of Apparent Tracer Ages in Heterogeneous Environments

The interpretation of apparent ages often assumes that a water sample is composed of a single age. In heterogeneous aquifers, apparent ages estimated with environmental tracer methods do not reflect mean water ages because of the mixing of waters from many flow paths with different ages. This is due to nonlinear variations in atmospheric concentrations of the tracer with time resulting in biases of mixed concentrations used to determine apparent ages. The bias of these methods is rarely reported and has not been systematically evaluated in heterogeneous settings. We simulate residence time distributions (RTDs) and environmental tracers CFCs, SF₆, ⁸⁵Kr, and ³⁹Ar in synthetic heterogeneous confined aquifers and compare apparent ages to mean ages. Heterogeneity was simulated as both K‐field variance (σ²) and structure. We demonstrate that an increase in heterogeneity (increase in σ² or structure) results in an increase in the width of the RTD. In low heterogeneity cases, widths were generally on the order of 10 years and biases generally less than 10%. In high heterogeneity cases, widths can reach 100 s of years and biases can reach up to 100%. In cases where the temporal variations of atmospheric concentration of individual tracers vary, different patterns of bias are observed for the same mean age. We show that CFC‐12 and CFC‐113 ages may be used to correct for the mean age if analytical errors are small.     

Analytical solutions for transport of decaying solutes in rivers with transient storage

Analytical solutions are presented for solute transport in rivers including the effects of transient storage and first order decay. The solute transport model considers an advection–dispersion equation for transport in the main channel linked to a first order mass exchange between the main channel and the transient storage zones. In case of a conservative tracer, it is shown that different analytical solutions presented in the literature are mathematically identical. For non-conservative solutes, first order decay reactions are considered with different reaction rate coefficients in the main river channel and in the dead zones. New analytical solutions are presented for different boundary conditions, i.e. instantaneous injection in an infinite river reach, and variable concentration time series input in a semi-infinite river reach. The correctness and accuracy of the analytical solutions is verified by comparison with the OTIS numerical model. The results of analytical and numerical approaches compare favourably and small differences can be attributed to the influence of boundary conditions. It is concluded that the presented analytical solutions for solute transport in rivers with transient storage and solute decay are accurate and correct, and can be usefully applied for analyses of tracer experiments and transport characteristics in rivers with mass exchange in dead zones.     

The impact of boundary on the fractional advection–dispersion equation for solute transport in soil: Defining the fractional dispersive flux with the Caputo derivatives

The inherent heterogeneity of geological media often results in anomalous dispersion for solute transport through them, and how to model it has been an interest over the past few decades. One promising approach that has been increasingly used to simulate the anomalous transport in surface and subsurface water is the fractional advection–dispersion equation (FADE), derived as a special case of the more general continuous time random walk or the stochastic continuum model. In FADE, the dispersion is not local and the solutes have appreciable probability to move long distances, and thus reach the boundary faster than predicted by the classical advection–dispersion equation (ADE). How to deal with different boundaries associated with FADE and their consequent impact is an issue that has not been thoroughly explored. In this paper we address this by taking one-dimensional solute movement in soil columns as an example. We show that the commonly used FADE with its fractional derivatives defined by the Riemann–Liouville definition is problematic and could result in unphysical results for solute transport in bounded domains; a modified method with the fractional dispersive flux defined by the Caputo derivatives is presented to overcome this problem. A finite volume approach is given to numerically solve the modified FADE and its associated boundaries. With the numerical model, we analyse the inlet-boundary treatment in displacement experiments in soil columns, and find that, as in ADE, treating the inlet as a prescribed concentration boundary gives rise to mass-balance errors and such errors could be more significant in FADE because of its non-local dispersion. We also discuss a less-documented but important issue in hydrology: how to treat the upstream boundary in analysing the lateral movement of tracer in an aquifer when the tracer is injected as a pulse. It is shown that the use of an infinite domain, as commonly assumed in literature, leads to unphysical backward dispersion, which has a significant impact on data interpretation. To avoid this, the upstream boundary should be flux-prescribed and located at the upstream edge of the injecting point. We apply the model to simulate the movement of Cl⁻ in a tracer experiment conducted in a saturated hillslope, and analyse in details the significance of upstream-boundary treatments in parameter estimation.     

Two-dimensional power series solution for non-axisymmetrical transport in a radially convergent tracer test with scale-dependent dispersion

It has been known for many years that dispersivity increases with solute travel distance in a subsurface environment. The increase of dispersivity with solute travel distance results from the significant variation of hydraulic properties of heterogeneous media and was identified in the literature as scale-dependent dispersion. This study presents an analytical solution for describing two-dimensional non-axisymmetrical solute transport in a radially convergent flow tracer test with scale-dependent dispersion. The power series technique coupling with the Laplace and finite Fourier cosine transform has been applied to yield the analytical solution to the two-dimensional, scale-dependent advection–dispersion equation in cylindrical coordinates with variable-dependent coefficients. Comparison between the breakthrough curves of the power series solution and the numerical solutions shows excellent agreement at different observation points and for various ranges of scale-related transport parameters of interest. The developed power series solution facilitates fast prediction of the breakthrough curves at any observation point.     

Heat Tracing in a Fractured Aquifer with Injection of Hot and Cold Water

Heat as a tracer in fractured porous aquifers is more sensitive to fracture-matrix processes than a solute tracer. Temperature evolution as a function of time can be used to differentiate fracture and matrix characteristics. Experimental hot (50 °C) and cold (10 °C) water injections were performed in a weathered and fractured granite aquifer where the natural background temperature is 30 °C. The tailing of the hot and cold breakthrough curves, observed under different hydraulic conditions, was characterized in a log–log plot of time vs. normalized temperature difference, also converted to a residence time distribution (normalized). Dimensionless tail slopes close to 1.5 were observed for hot and cold breakthrough curves, compared to solute tracer tests showing slopes between 2 and 3. This stronger thermal diffusive behavior is explained by heat conduction. Using a process-based numerical model, the impact of heat conduction toward and from the porous rock matrix on groundwater heat transport was explored. Fracture aperture was adjusted depending on the actual hydraulic conditions. Water density and viscosity were considered temperature dependent. The model simulated the increase or reduction of the energy level in the fracture-matrix system and satisfactorily reproduced breakthrough curves tail slopes. This study shows the feasibility and utility of cold water tracer tests in hot fractured aquifers to boost and characterize the thermal matrix diffusion from the matrix toward the flowing groundwater in the fractures. This can be used as complementary information to solute tracer tests that are largely influenced by strong advection in the fractures.     

Analysis of tracer tomography using temporal moments of tracer breakthrough curves

Hydraulic/partitioning tracer tomography (HPTT) was recently developed by Yeh and Zhu [Yeh T-CJ, Zhu J. Hydraulic/partitioning tracer tomography for characterization of dense nonaqueous phase liquid source zones, Water Resour Res 2007;43:W06435. doi:10.1029/2006WR004877.] for estimating spatial distribution of dense nonaqueous phase liquids (DNAPLs) in the subsurface. Since discrete tracer concentration data are directly utilized for the estimation of DNAPLs, this approach solves the hyperbolic convection–dispersion equation. Solution to the convection–dispersion equation however demands fine temporal and spatial discretization, resulting in high computational cost for an HPTT analysis. In this work, we use temporal moments of tracer breakthrough curves instead of discrete concentration data to estimate DNAPL distribution. This approach solves time independent partial differential equations of the temporal moments, and therefore avoids solving the convection–dispersion equation using a time marching scheme, resulting in a dramatic reduction of computational cost. To reduce numerical oscillations associated with convection dominated transport problems such as in inter-well tracer tests, the approach uses a finite element solver adopting the streamline upwind Petrov–Galerkin method to calculate moments and sensitivities. We test the temporal moment approach through numerical simulations. Comparing the computational costs between utilizing moments and discrete concentrations, we find that temporal moments significantly reduce the computation time. We also find that tracer moment data collected through a tomographic survey alone are able to yield reasonable estimates of hydraulic conductivity, as indicated by a correlation of 0.588 between estimated and true hydraulic conductivity fields in the synthetic case study.     

Macropores as preferential flow paths in meander bends

Macropores are subsurface connected void spaces caused by processes such as fracture of soils, micro‐erosion, and fauna burrows. They are common near streams (e.g. hyporheic and riparian zones) and may act as preferential flow paths between surface and groundwaters, affecting hydrologic and biogeochemical processes. We tested the hydrologic function of macropores by constructing an artificial macropore within the saturated zone of a meander bend (open macropore, ‘OM’) and later filling its upstream end (partially filled macropore, ‘PFM’). For each treatment, we injected saline tracer at an upgradient monitoring well within the meander and monitored downgradient hydraulics and tracer transport. Pressure transducers in monitoring wells indicated hydraulic gradients within the meander were 32% higher perpendicular to and 6% higher parallel to the macropore for the OM than for the PFM. Additionally, hydraulic conductivities measured via falling head tests were 29 to 550 times higher along the macropore than in nearby sediment. We used electrical conductivity probes in wells and electrical resistivity imaging to track solute transport. Transport velocities through the meander were on average 9 and 21% higher (per temporal moment analysis and observed tracer peak, respectively) for the OM than for the PFM. Furthermore, temporal moments of tracer breakthrough analysis indicated downgradient longitudinal dispersion and breakthrough tracer curve tailing were on average 234% and 182% higher for the OM, respectively. This suggests the OM enabled solute transport at overall shorter timescales than the matrix but also increased tailing. Our results demonstrate the importance of macropores to meander bend hydrology and solute transport. Copyright © 2012 John Wiley & Sons, Ltd.     

Simulation of less‐mobile porosity dynamics in contrasting sediment water interface porous media

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.     

Analysis of solute transport with a hyperbolic scale-dependent dispersion model

An empirical hyperbolic scale-dependent dispersion model, which predicts a linear growth of dispersivity close to the origin and the attainment of an asymptotic dispersivity at large distances, is presented for deterministic modelling of field-scale solute transport and the analysis of solute transport experiments. A simple relationship is derived between local dispersivity, which is used in numerical simulations of solute transport, and effective dispersivity, which is estimated from the analysis of tracer breakthrough curves. The scale-dependent dispersion model is used to interpret a field tracer experiment by nonlinear least-squares inversion of a numerical solution for unsaturated transport. Simultaneous inversion of concentration-time data from several sampling locations indicates a linear growth of the dispersion process over the scale of the experiment. These findings are consistent with the results of an earlier analysis based on the use of a constant dispersion coefficient model at each of the sampling depths.     

Modelling spatial variations in dissolved oxygen in fine‐grained streams under uncertainty

This paper presents a novel triple‐layer model, called VART DO‐3L, for simulation of spatial variations in dissolved oxygen (DO) in fine‐grained streams, characterized by a fluid mud (fluff or flocculent) layer (an advection‐dominated storage zone) as the interface between overlying stream water and relatively consolidated streambed sediment (a diffusion‐dominated storage zone). A global sensitivity analysis is conducted to investigate the sensitivity of VART DO‐3L model input parameters. Results of the sensitivity analysis indicate that the most sensitive parameter is the relative size of the advection‐dominated storage zones (A_s/A), followed by a lumped reaction term (R) for the flocculent layer, biological reaction rate (μ_o) in diffusive layer and biochemical oxygen demand concentration (L) in water column. In order to address uncertainty in model input parameters, Monte Carlo simulations are performed to sample parameter values and to produce various parameter combinations or cases. The VART DO‐3L model is applied to the Lower Amite River in Louisiana, USA, to simulate vertical and longitudinal variations in DO under the cases. In terms of longitudinal variation, the DO level decreases from 7.9 mg l at the Denham Springs station to about 2.89 mg l^?1 at the Port Vincent station. In terms of vertical variation, the DO level drops rapidly from the overlying water column to the advection‐dominated storage zone and further to the diffusive layer. The DO level (C_F) in the advective layer (flocculent layer) can reach as high as 40% of DO concentration (C) in the water column. The VART DO‐3L model may be applied to similar rivers for simulation of spatial variations in DO level. Copyright © 2014 John Wiley & Sons, Ltd.     

Residence time of bedform-driven hyporheic exchange

Biogeochemical and ecological transformations in hyporheic zones are dependent on the timing of hyporheic exchange. We show through linked modeling of open channel turbulent flow, groundwater flow, and solute transport that the residence time distributions of solutes advected by hyporheic flow induced by current–bedform interaction follow power-laws. This tailing behavior of solutes exiting the sediments is explained by the presence of multiple path lengths coupled with very large variability in Darcy flow velocity, both occurring without heterogeneity in sediment permeability. Hyporheic exchange through bedforms will result in short-time fractal scaling of stream water chemistry.