Next Article in Journal
Determining Characteristic Vegetation Areas by Terrestrial Laser Scanning for Floodplain Flow Modeling
Next Article in Special Issue
On the Modeling of Bank Storage in a Groundwater Model: The April, 1983, Flood Event in the Neuwieder Becken (Middle Rhine)
Previous Article in Journal
Development of a Comprehensive Fouling Model for a Rotating Membrane Bioreactor System Treating Wastewater
Previous Article in Special Issue
Effects of Streambed Conductance on Stream Depletion
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Benthic Uptake Rate due to Hyporheic Exchange: The Effects of Streambed Morphology for Constant and Sinusoidally Varying Nutrient Loads

by
Daniele Tonina
1,*,†,
Alessandra Marzadri
1,† and
Alberto Bellin
2,†
1
Center for Ecohydraulics Research, University of Idaho, 322 E Front Street, suite 322, Boise, ID 83702, USA
2
Department of Civil, Environmental and Mechanical Engineering, University of Trento, Via Mesiano 77, Trento 38123, Italy
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Water 2015, 7(2), 398-419; https://doi.org/10.3390/w7020398
Submission received: 28 September 2014 / Accepted: 16 January 2015 / Published: 23 January 2015
(This article belongs to the Special Issue Surface Water Groundwater Interactions: From Theory to Practice)

Abstract

:
Hyporheic exchange carries reactive solutes, which may include biological oxygen demand (BOD), dissolved oxygen (DO) and reactive dissolved inorganic nitrogen (Nr), into the sediment, where biochemical reactions consume DO. Here, we study the impact of streambed morphology, stream-reactive solute loads and their diel oscillations on the DO benthic uptake rate (BUR) due to hyporheic processes. Our model solves the hyporheic flow field and the solute transport equations analytically, within a Lagrangian framework, considering advection, longitudinal diffusion and reactions modeled as first order kinetics. The application of the model to DO field measurements over a gravel bar-pool sequence shows a good match with measured DO concentrations with an overall agreement of 58% and a kappa index of 0.46. We apply the model to investigate the effects of daily constant and sinusoidally time varying stream BOD, DO and Nr loads and of the morphodynamic parameters on BUR. Our modeling results show that BUR varies as a function of bedform size and of nutrient loads and that the hyporheic zone may consume up to 0.06% of the stream DO at the pool-riffle bedform scale. Daily oscillations of stream BOD and DO loads have small effects on BUR, but may have an important influence on local hyporheic processes and organisms’ distribution.

Graphical Abstract

1. Introduction

Human activities have recently increased nutrient loads in aquatic environments, which would typically be nutrient-limited [1]. This increase has led to water body eutrophication, thereby threatening the integrity of aquatic ecosystems [1,2,3,4,5]. In freshwater, dissolved inorganic reactive nitrogen, Nr, primarily under the form of nitrates, NO3 and ammonium NH4+ and dissolved organic matter have cascading effects on dissolved oxygen (DO) concentrations [6,7]. For example, nitrogenous biological oxygen demand (NBOD = LN) is associated with the transformation of NH4+ to NO3, whereas respiration and metabolic activities are characterized by the carbonaceous biological oxygen demand (CBOD = LC). High CBOD is considered the most increasing water pollution problem in the world [8].
These quantities show intertwined complex dynamics over a hierarchy of scales, from daily to seasonal [9,10]. Sinusoidal oscillations have been observed to approximate the daily variations well for DO [10], temperature [11,12] and BOD [13]. Daily oscillations of DO are due to photosynthesis/respiration plant cycles, those of temperature to the day/night cycle, while BOD variations are chiefly correlated with discharge from wastewater treatment plants.
These oscillations may affect the quality of both in-stream and pore waters, because in-stream solutes are exchanged between streams and streambed sediments through diffusive and advective processes [14,15,16], with the latter being the dominant mechanism in riverine systems due to hyporheic exchange (e.g., [17,18,19,20]). In fact, stream water downwells into the alluvium in high-head zones and upwells into the stream in low-head zones [20]. Downwelling fluxes carry reactive solutes and DO, which is consumed by the microbial communities within the hyporheic zone [21,22,23,24], resulting in upwelling stream waters with lower DO concentrations than those of stream waters [21,23,25,26,27,28,29,30].
Following the previous works of Rutherford et al. [16], we propose that this reduction in DO concentrations between downwelling and upwelling waters could be used to define the oxygen consumption due to biogeochemical activities occurring within the sediment. This oxygen consumption is also referred to as the benthic uptake rate, BUR.
Although field experiments show the importance of hyporheic residence time in riverine-hyporheic biogeochemistry [22,31,32,33], models coupling surface-subsurface exchange fluxes with biogeochemical reactions to quantify BUR are limited [16], with more publications addressing nitrate cycling (e.g., [31,33,34,35,36,37]). Rutherford et al. [16] analyzed BUR for hyporheic exchange induced by two-dimensional dune-like bedforms, while neglecting solute diffusion within the streambed and nitrification processes. Boano et al. [38] numerically analyzed biogeochemical zonation of intra-meander bars due to differences in the residence time of hyporheic water, but they did not evaluate stream oxygen consumption due to this process, nor the influence of stream bedforms on BUR. In a similar fashion, Marzadri et al. [31,35] proposed a simplified, yet effective, model of the nitrogen cycle within the hyporheic zone, but neglecting BOD loads. Similarly, other recent investigations focused on nitrate cycling and modeled dissolved oxygen transport to quantify hyporheic zone redox conditions (e.g., [34,36,37]). More recently, we investigated the effect of streambed morphology on the hyporheic thermal regime and studied the effects of daily temperature variations of stream water on dissolved oxygen distribution within the hyporheic zone [36]. In particular, we showed that mean daily temperatures suffice to characterize the hyporheic zone biogeochemistry at the bedform scale, while daily temperature oscillations may affect local biogeochemistry in sediment volumes with short hyporheic residence times. However, we did not quantify the effects of streambed size nor of BOD loads on BUR.
In the present work, we close this gap by extending the Rutherford et al. [16] model to gravel-bed rivers with a three-dimensional alternate bar morphology, which is a ubiquitous and ecologically important bedform [39,40]. We develop a semi-analytical model for BUR, coupling surface and subsurface fluxes with aerobic reactions (respiration and nitrification) within the streambeds. Our model considers solute advection, longitudinal diffusion and first-order kinetics, under the assumption of low nutrient concentrations, with temperature-dependent reaction rate coefficients.
Our first objective is to test the capability of our model to predict DO concentration patterns in a natural setting by comparing predicted and measured DO concentrations. Here, we use the field data reported by Malard et al. [41] and by Rouch [42] for a 15 m-long pool-riffle reach of the Lachein River (Pyrennees, France). Our second objective is to show how the effects of the daily time varying stream BOD, DO and Nr loads (constant or sinusoidally) change the global hyporheic response in terms of BUR and how the local oxic hyporheic pattern changes according to the DO distribution.

2. Materials and Methods

2.1. Hyporheic Hydraulic Model

We model the hyporheic flow as a Darcian flow [18,43], with the seepage velocity given by Darcy’s law: u = −Kh∇h/n, where Kh is the hydraulic conductivity, n is the sediment porosity and h is the energy head. We assume homogenous and isotropic hydraulic properties of the streambed sediment with stationary flow conditions, such that the governing equation is the following:
2 h = 0
We solve Equation (1) for a three-dimensional alternate bar morphology under the assumptions that the lateral sides (stream banks) and the bottom of the alluvium are impervious and that the hydraulic head gradient is equal to the streambed slope. The stream water surface elevation, which approximates the energy head distribution in the case of submerged bars, is decomposed through the fast Fourier transform. Notice that water surface elevation can be either measured [35,36], obtained through numerical models [40,44] or obtained from analytical hydrodynamic solutions [18,45]. Owing to the linearity of Equation (1), the main components of the Fourier decomposition are applied separately as Dirichlet boundary conditions at the upper planar surface that delimitates the subsurface domain. Previous studies showed that replacing the streambed surface with the plane coinciding with the average streambed elevation, as done here, has little effects on the flow and residence time distributions [18,45]. The solutions obtained with the single FFT components as the boundary condition are then superimposed to obtain the final head distribution within the subsurface domain. Here, we represent the head distribution from two different sources: field measured heads reported by Malard et al. [41] over a pool-riffle unit and analytically predicted by the model of Colombini et al. [46] for the case of a fully-submerged three-dimensional alternate bar in equilibrium with the stream discharge (Figure 1).
Figure 1. Sketch of (a) the planar view and (b) the cross-section of the three-dimensional alternate bar morphology; and (c) diagram of the biogeochemical model. DO, dissolved oxygen; CBOD, carbonaceous biological oxygen demand; NBOD, nitrogenous biological oxygen demand.
Figure 1. Sketch of (a) the planar view and (b) the cross-section of the three-dimensional alternate bar morphology; and (c) diagram of the biogeochemical model. DO, dissolved oxygen; CBOD, carbonaceous biological oxygen demand; NBOD, nitrogenous biological oxygen demand.
Water 07 00398 g001
Once the flow field is known, the transport equation is solved by means of the particle tracking technique [47]. The works of Marzadri et al. [45,48] on pool-riffle-induced hyporheic flow show that residence time distributions and biogeochemical transformations depend on bedform size through the ratio Y B M * = HBM/Y0 between the bedform amplitude HBM and the water hydraulic depth, Y0.

2.2. Biogeochemical Model

Following the approach of Marzadri, et al. [31], the governing equation for DO transport within the hyporheic zone is written in a Lagrangian framework by applying mass conservation along a stream tube [49]:
D O t + D O τ = α  2 D O K L L C K N L N K R D O
where α = D/u2 is the ratio between the longitudinal diffusion coefficient (D) and the module square of the mean velocity (u2) (u = |u|). In addition, KL, KN and KR are the (linear) reaction rates of carbonaceous BOD, nitrogenous BOD and respiration, respectively (Figure 1c). The linear reaction rates of the consumption of DO are used, owing to its relatively low concentration. The last term on the right-hand side of Equation (2) accounts for dissolved oxygen consumption due to other biological sources besides the downwelling carbonaceous and nitrogenous BODs. The model neglects the lateral exchange of mass between adjacent stream tubes, and the reaction rates are assumed to be temperature dependent. Because the effects of daily temperature oscillations are small on the biogeochemical response of the hyporheic zone at the bedform scale [50], reaction rates are assumed to change with the mean daily temperature, according to the Arrhenius law [51].
We assume full mixing within the stream, such that the solute concentrations are spatially constant within the downwelling surface, while they vary in time. The initial zero concentration of all solutes is assumed within the hyporheic zone. With these hypotheses, we solve Equation (2) for constant and sinusoidally varying DO, LC and LN stream water concentrations. The latter assumption derives from experimental results showing that daily variations of nutrient concentrations are well approximated by a sinusoidal function [13]. Solutions of Equation (2) under these initial and boundary conditions are reported in Appendix A.
We introduce two dimensionless numbers: the first is the ratio between the dissolved oxygen concentration in the stream water, DO0, and the total oxygen demand (L0C + L0N), RDO = DO0/(L0C + L0N), and the second is the ratio between carbonaceous and nitrogenous BOD concentrations, RBN = L0C/L0N, with the objective to facilitate the analysis of the effects of nutrient loads on BUR.

2.3. Benthic Uptake Rate

BUR is defined as the difference between mass fluxes of DO entering and leaving the streambed sediments, which are given by the first and the second term on the right-hand side of Equation (3), respectively:
B U R = 1 τ m a x [ t 0 t 0 + τ m a x 1 A b b 0 λ H ( x , y ) q ( x , y ) D O 0 ( t ) d x d y d t t 0 t 0 + τ m a x 1 A b b 0 λ H ( x , y ) q ( x , y )    D O ( τ ( x , y ) , t )    d x d y d t ]
where A = 2 b λ is the streambed area, τ(x,y) is the hyporheic residence time for a particle downwelling at x, y, z = 0 and upwelling after the (residence) time τ since injection, λ is the bedform wavelength (the distance between two pools) and 2b the width of the channel. The time τmax is the maximum hyporheic residence time for the given bedform and defines the time interval over which BUR is quantified. In addition, t0 is the initial time, which is set to a value larger than maximum residence time τmax in order to eliminate in Equation (3) the effect of the initial conditions. The step function H is unitary within the downwelling areas and is zero outside. Notice that mass balance along a single streamline has been used in the second right-hand term of Equation (3) to transform the integral over the upwelling area in an equivalent integral over the downwelling area [31,35], under the assumption that hyporheic flow is at the steady state. The steady-state flow condition allows us to rewrite Equation (3) in the following form:
B U R = D O ¯ 0    q ¯ 1 2 b  λ b b 0 λ H ( x , y ) q ( x , y )    D O ¯ ( τ ( x , y ) )    d x d y q ¯ = 1 2 b  λ b b 0 λ H ( x , y ) q ( x , y ) d x d y
where q ¯ is the mean downwelling flux averaged over the entire streambed area, and time averaged dissolved oxygen concentrations are defined as:
D O ¯ 0 = 1 τ max t 0 t 0 + τ max D O 0 ( t ) d t
and:
D O ¯ ( τ ) = 1 τ max t 0 t 0 + τ max D O ( τ , t ) d t
at the downwelling and upwelling areas, respectively.
Because the objective of the present analysis is to quantify the contribution of the hyporheic zone to the consumption of stream DO, we make the BUR dimensionless with the mean DO mass flux in the stream:
B U R = 2 b λ B U R 1 τ max t 0 t 0 + τ max D O 0 ( t ) Q s t r e a m d t = 2 b λ B U R Q s t r e a m D O 0 ¯
where Qstream is the stream discharge. Similarly, the BUR can be made dimensionless with respect to the hyporheic flux by replacing in Equation (7) Qstream with q ¯ A:
B U R H Z = B U R q ¯ D O ¯ 0
Hereafter, this dimensionless expression of BUR will be indicated with B U R H Z * .
The following two end members can be defined for BUR* when all of the upwelling fluxes are anoxic, such that all of the oxygen carried by the stream is consumed within the hyporheic zone:
B U R = q ¯ 2 b λ Q s t r e a m f o r R D O < 1 B U R = R D O 1 q ¯ 2 b λ Q s t r e a m f o r R D O > 1
which have been obtained by considering that carbonaceous and nitrogen oxygen consumption is limited by the DO of stream water, such that B U R = q ¯   D O 0 ¯ when RD0 < 1 and BUR = L0C + L0N when RD0 > 1. Equation (8) provides a first approximation analysis for BUR*, which depends only on the hyporheic ( q ¯ 2) and stream (Qstream) discharges, the hyporheic residence time distribution and the nutrient loads (RDO = DO0/(L0C + L0N)) in the stream. The hyporheic values, τmax and q ¯ , can be estimated with models presented in the literature [18,19,45,52,53,54,55].

2.4. Field Data

The works of Malard et al. [41] and Rouch [42] report both the near-bed head and the dissolved oxygen concentration distributions at the sediment-water interface, supplemented by streambed sediment hydraulic conductivity and porosity along a 15 m-long pool-riffle bar reach of the Lachein stream (Pyrennees, France). In the present work, the head distribution derived from the contour lines depicted in Figure 3c of Malard et al. [41] by superimposing a regular two-dimensional grid is first decomposed via FFT (Figure 2). Owing to the linearity of the flow Equation (1) and Darcy’s law, the most relevant FFT components are used as the upper Dirichlet boundary condition for modeling the 3D hyporheic flow field, through the solution of Laplace and Darcy equations, separately for each component, under the assumption of homogeneous and isotropic alluvium, as described in Section 2.1. The resulting velocity fields are then superimposed to obtain the final three-dimensional velocity field distribution. We set the alluvium thickness to zd = 2 m, which is deep enough to make negligible its effect on the hyporheic flow field, the hydraulic conductivity to Kh = 8.39 × 10−4 m/s, which is the mean of the hydraulic conductivity values reported in Figure 3b of Malard et al. [41], and the porosity to n = 0.2664, which is the mean of the porosity values reported by Rouch [42] for this study site. The residence time of 8192 particles injected within the downwelling area was obtained by particle tracking with the analytical velocity obtained as explained above. We do not account for particles whose trajectories do not upwell in the stream.
Figure 2. Head distribution in the 15 m-long reach of the Lachein stream measured by Malard et al.’s [41] contour lines and reconstructed by FFT (color plot).
Figure 2. Head distribution in the 15 m-long reach of the Lachein stream measured by Malard et al.’s [41] contour lines and reconstructed by FFT (color plot).
Water 07 00398 g002
For the biogeochemical model, we assumed a constant dissolved oxygen concentration in the stream water, DO0 = 10 mg/L, which is the mean of the measured values over the downwelling areas, and a saturated DO concentration of DOsat = 11.83 mg/L, which corresponds to the observed water temperature of T = 8 °C. In addition, we set the reaction rate of respiration equal to the mean value calculated from the field data reported by Beaulieu et al. [56] for streams with a water temperature near 20 °C: KR(20 °C) = 5.01 d−1. We neglect all of the other parameters, DOlim, LC, LN and α, whose values were not reported in the work of Malard et al. [41] and Rouch [42]. The coefficient KR implicitly accounts for the consumption of DO due to both carbonaceous and nitrogenous BOD, whose values were not reported, and thus, we set KL and KN equal to 0. These working assumptions allow us to use only measured values, avoiding the need to calibrate the model.
We used the error matrix technique, which is a common GIS (Geographic Information System) method for map comparison [57,58,59], to assess the accuracy of the modeled spatially-distributed DO concentrations against field observations. The error matrix provides two indexes of agreement between predicted and observed spatially-distributed quantities on a cell-by-cell basis: an overall agreement, OA, and a kappa index. The OA is calculated as the proportion of correctly-predicted cells versus the total number of cells. Differently from OA, the kappa index removes the effect of by-chance agreement. In fact, the probability that two successive random extractions of a categorical variable belong to the same category is larger than zero, by definition. Since the kappa index filters out the match of predicted and observed DO cells due to random coincidence, its value is typically lower than OA [60]. Both indexes vary from 0 to 1 with 0 being poor and 1 perfect agreement. For the K statistic, Landis and Koch [59] suggested the following classification: <0 = poor agreement, 0–0.19 slight, 0.20–0.39 fair, 0.4–0.59 moderate, 0.60–0.79 substantial and 0.80–1.00 almost perfect agreement. We calculated OA and K with the DO concentrations divided into the 3 categories reported by Malard et al. [41]. The predicted and observed DO category needs to coincide exactly to be counted as a match and not an error. This is a strict method, because it does not take into account if the match is a little bit offset spatially; for instance, a nearby cell could have the correct value [60].

2.5. Numerical Simulations

We considered 18 streams characterized by an alternate bar morphology, whose bedforms are in equilibrium with the flow discharge. The bedform dimensions, whose morphodynamic parameters are reported in Table 1, are within those typically encountered in natural streams [39,61,62,63]. For each stream, we ran 10 simulations, changing the nutrient loads (Table 2). We ran a set of simulations with constant loads and a set with sinusoidal varying CBOD, NBOD and DO concentrations with the same mean value as the constant case to explore the effects of time-dependent nutrient releases.
Table 1. Morphodynamic parameters of the pool-riffle bedform. β is the aspect ratio (β = b/Y0); θ is the Shields’ number; ds is the relative submergence of the particles; d50 is the median grain size; Y B M * is the dimensionless mean flow depth; s is the streambed slope; is the bar amplitude; λ is the bar length; and Qstream is the stream water discharge.
Table 1. Morphodynamic parameters of the pool-riffle bedform. β is the aspect ratio (β = b/Y0); θ is the Shields’ number; ds is the relative submergence of the particles; d50 is the median grain size; Y B M * is the dimensionless mean flow depth; s is the streambed slope; is the bar amplitude; λ is the bar length; and Qstream is the stream water discharge.
Testβ (-)θ (-)dS (-)d50 (m) Y B M * (-)S (%)HBM (m)λ (m)Qstream (m3·s−1)
1130.080.2000.010.332.640.158.760.057
2130.080.1500.010.391.980.1711.520.111
3130.080.1200.010.451.580.1914.230.185
4130.080.1050.010.481.390.2016.160.251
5130.080.0950.010.501.250.2117.770.314
6130.080.0850.010.531.120.2219.750.404
7130.080.0750.010.560.990.2322.260.537
8130.080.0650.010.600.860.2625.530.738
9130.080.0600.010.620.790.2727.570.883
10130.080.0550.010.650.730.2829.981.072
11130.080.0500.010.680.660.2932.871.325
12130.080.0450.010.710.590.3136.401.675
13130.080.0400.010.750.530.3340.822.174
14130.080.0350.010.810.460.3346.522.902
15130.080.0300.010.870.400.3854.134.101
16130.080.0250.010.960.330.4264.826.121
17130.080.0200.011.080.260.4680.979.977
18130.080.0100.011.740.130.56163.1445.037
Table 2. Biogeochemical parameters used in the simulations. RDO is the ratio between dissolved oxygen concentration (DO) and the total oxygen demand; RBN is the ratio between carbonaceous (LC) and nitrogenous (LN) BOD; L0C and L0N are the constant stream concentrations of carbonaceous (CBOD) and nitrogenous (NBOD) BOD, respectively. LAC and LAN are the stream daily fluctuations of CBOD and NBOD, respectively. KL, KN and KR are the decay rates of CBOD, NBOD and DO, respectively. In all of the simulations, the concentration of dissolved oxygen concentration DO0 = 10 mg/L, while the stream daily fluctuation of oxygen is DOA = 4 mg/L.
Table 2. Biogeochemical parameters used in the simulations. RDO is the ratio between dissolved oxygen concentration (DO) and the total oxygen demand; RBN is the ratio between carbonaceous (LC) and nitrogenous (LN) BOD; L0C and L0N are the constant stream concentrations of carbonaceous (CBOD) and nitrogenous (NBOD) BOD, respectively. LAC and LAN are the stream daily fluctuations of CBOD and NBOD, respectively. KL, KN and KR are the decay rates of CBOD, NBOD and DO, respectively. In all of the simulations, the concentration of dissolved oxygen concentration DO0 = 10 mg/L, while the stream daily fluctuation of oxygen is DOA = 4 mg/L.
TestRDO (-)RBN (-)L0C (mg/L)LAC (τ,0) (mg/L)L0N (mg/L)LAN (τ,0) (mg/L)KL (d1)KN (d1)KR (d1)
A1.20.52.781.115.562.220.330.260.16
B1.00.53.331.336.67.2.670.330.260.16
C0.80.54.161.678.323.330.330.260.16
D0.50.56.672.6713.335.330.330.260.16
E0.20.516.666.6733.3313.330.330.260.16
F1.22.05.562.222.781.110.330.260.16
G1.02.06.672.673.331.330.330.260.16
H0.82.08.323.334.161.670.330.260.16
I0.52.013.345.336.672.670.330.260.16
J0.22.033.3313.3316.676.670.330.260.16
These simulations are tailored to investigate the relationship between BUR*, RBN and RDO, which identify the nutrient loads, and Y B M * , an index of stream size. Large streams are typically characterized by large Y*BM and long residence times, whereas small streams by small Y*BM and short residence times [31].

3. Results and Discussion

3.1. Comparison with Field Data

Figure 3 shows the contour map of the computed hyporheic residence times (Figure 3a), the categorical map of measured dissolved oxygen concentration (red contours) (Figure 3b, which corresponds to Figure 4d in Malard et al. [41]) and the map of model-predicted dissolved oxygen concentrations (Figure 3c) at the sediment-water interface. The overall agreement value between the predicted and observed DO concentration distribution is 0.58, and the kappa index is 0.46, which correspond to a good agreement and are comparable to the values reported in the GIS literature [57,59,64]. Because, to our knowledge, this is the first comparison between observed and modeled spatially-distributed hyporheic quantities, we cannot compare these values with others reported in the hyporheic literature. This result indicates that hyporheic processes are those that chiefly control the DO concentration distribution at and below the water-sediment interface and, consequently, affect the habitat of benthic and hyporheos organisms. The predicted patterns of τ and DO concentrations also match the contours of the categorical map of DO concentrations measured in the field (Figure 4d in [41]). This further confirms that residence time distribution is a key factor in the fate of reactive solutes within the hyporheic zone. In our case, the transformations depend entirely on the reaction time, because the biogeochemical reaction rates are spatially and temporally constant.
Figure 3. Maps of: (a) the computed residence time distribution of water within the hyporheic zone, τup; (b) the measured dissolved oxygen concentrations, (DO), in % of saturation reproduced from Figure 4d in Malard et al. [41]; and (c) the model-predicted DO along the 15-m reach of the Lachein stream.
Figure 3. Maps of: (a) the computed residence time distribution of water within the hyporheic zone, τup; (b) the measured dissolved oxygen concentrations, (DO), in % of saturation reproduced from Figure 4d in Malard et al. [41]; and (c) the model-predicted DO along the 15-m reach of the Lachein stream.
Water 07 00398 g003
Malard et al. [41] observed that organisms colonized preferential areas of the streambed depending on DO concentrations and described the behavior in terms of morphological units: bar, pool and bar fringe. Patches with a high dissolved oxygen concentration (>74% of saturation) coincide with depositional areas, such as bars, which are mostly predicted as downwelling zones, and bar fringes, which are predicted as upwelling areas with a short residence time. These zones show dense communities of harpacticoid and hypogean species; whereas streambed areas with low DO concentrations, which are identified as erosional areas by Malard et al. [41], like the pool, and for which our model predicts upwelling areas with long residence times, have fewer harpacticoids and hypogean organisms. All of these results suggest that the near-bed DO concentration is chiefly controlled by stream bedforms, which therefore control the habitat within and at the streambed sediment.

3.2. Effect of Stream Morphology

Figure 4a shows BUR*, i.e., the BUR made dimensionless with respect to the stream flux, as a function of the dimensionless depth, Y*BM, for hyporheic zones induced by three-dimensional pool-riffle morphology for RDO = 1.2 and RBN = 0.5. In addition, Figure 4b shows how the BUR ( B U R H Z * ) behaves when it is normalized with respect to the hyporheic flux.
Figure 4. (a) Stream, BUR*, and (b) hyporheic, BUR*HZ, dimensionless BUR, as a function of dimensionless mean flow depth, Y B M * , for three-dimensional alternate bars with RDO = 1.2 and RBN = 0.5, with RDO being the ratio between in-stream DO (DO0) and total oxygen demand (L0C + L0N) and RBN the ratio between in-stream carbonaceous and nitrogenous BOD concentrations, RBN = L0C/L0N. BUR, benthic uptake rate.
Figure 4. (a) Stream, BUR*, and (b) hyporheic, BUR*HZ, dimensionless BUR, as a function of dimensionless mean flow depth, Y B M * , for three-dimensional alternate bars with RDO = 1.2 and RBN = 0.5, with RDO being the ratio between in-stream DO (DO0) and total oxygen demand (L0C + L0N) and RBN the ratio between in-stream carbonaceous and nitrogenous BOD concentrations, RBN = L0C/L0N. BUR, benthic uptake rate.
Water 07 00398 g004
BUR* initially increases with river size: it reaches a peak and then decreases for large streams. The peak for BUR* is caused by the initial faster increase of downwelling fluxes with long hyporheic residence times, which potentially allow more DO consumption, than stream discharge, which also increases with stream size [54]. Eventually, the increase in stream discharge is large enough to offset the increase in hyporheic exchange, such that BUR* decreases.
Differently from BUR*, BUR*HZ increases monotonically with river size, because of the longer residence time that water spends within the hyporheic zone as the size of the river increases. Whereas the BUR*HZ increase is fast around Y*BM = 1, it slows down significantly at larger values of Y*BM. This suggests that short flow paths may be always present, reducing the efficiency of the hyporheic zone to consume DO to a value smaller than 100%, even for large streams.
Figure 5 shows a snapshot of hyporheic DO concentrations along the streamlines of the 3D hyporheic flow field provided by our model. Within the alluvium, the concentration of DO decreases with depth according to the time that water spends traveling along a given streamline.
The longer the path is, the greater is the probability to find streamlines with DO < 2 mg/L (red path of the trajectory) and, thus, volumes of sediment in anaerobic conditions. The complex DO distribution depicted in Figure 5 demonstrates the importance of considering the entire three-dimensional flow patterns when studying alternate bar morphologies. As the stream dimension increases, the volume of streambed sediment saturated with stream water by the hyporheic flow increases and, consequently, so do the residence times. Thus, in large streams, the red portion of the flow lines would dominate. However, the presence of short residence times at the transitional zone between upwelling and downwelling areas would still have a blue tone. This supports the results of Figure 4, which suggest that B U R H Z * may be limited to a value below 100%. However, this upper threshold depends also on biogeochemical rates, because faster rates would require shorter times to reach anoxic conditions [32,65].
Figure 5. Snapshot of dissolved oxygen concentrations (DO) simulated within the hyporheic zone of a river with the morphological parameters of Test 5 in Table 1.
Figure 5. Snapshot of dissolved oxygen concentrations (DO) simulated within the hyporheic zone of a river with the morphological parameters of Test 5 in Table 1.
Water 07 00398 g005

3.3. Effects of Constant BOD and DO Loads

Figure 6 shows BUR* as a function of Y B M * in the case of three-dimensional alternate bars for several RDO values with RBN = 0.5 (Figure 6a) and RBN = 2 (Figure 6b). Keeping in mind the discussion in Section 3.2 concerning the mechanism leading to the peak in the BUR* (Figure 4a), here, we can observe the effect of nutrient concentrations on the location of this peak. In particular, as RDO decreases (because stream DO decreases or L0C + L0N increases), the total oxygen demand (LC + LN) is larger than the stream DO concentrations; the BUR* peak increases, and it occurs in smaller streams (smaller Y B M * ) for low rather than high RDO. For instance, the BUR* peak occurs at Y B M * = 0.75 for RDO = 1.2 and at Y*BM = 0.48 for RDO = 0.2, regardless of RBN. Both cases with a different RBN show similar patterns with slightly higher BUR* associated with larger RBN. The faster kinetics of LC than LN justifies this behavior.
Figure 6. Dimensionless stream BUR, BUR*, as a function of the dimensionless mean flow depth, Y B M * , and RDO for (a) RBN = 0.5 and (b) RBN = 2.0, with RDO the ratio between in-stream DO (DO0) and total oxygen demand (L0C + L0N) and RBN the ratio between in-stream carbonaceous and nitrogenous BOD concentrations, RBN = L0C/L0N.
Figure 6. Dimensionless stream BUR, BUR*, as a function of the dimensionless mean flow depth, Y B M * , and RDO for (a) RBN = 0.5 and (b) RBN = 2.0, with RDO the ratio between in-stream DO (DO0) and total oxygen demand (L0C + L0N) and RBN the ratio between in-stream carbonaceous and nitrogenous BOD concentrations, RBN = L0C/L0N.
Water 07 00398 g006
However, the effect of nutrient loads on BUR* is important only for small streams. As the stream size increases, both the RDO and RBN effects become negligible, and all curves tend to the same value. Our simulation suggests that the threshold is near Y B M * = 1.

3.4. Effects of Sinusoidally Varying and Constant Stream Solute Concentrations

Figure 7 shows BUR* for sinusoidally varying and constant LC and LN loads as a function of Y B M * for RBN = 0.5 with RDO = 0.2 (Figure 7a) and RDO = 1.2 (Figure 7b). The most important result shown in Figure 7a,b is that including daily variations of the stream DO provides little benefit with respect to the use of a constant stream DO for large streams. Some effects of the time-varying in-stream concentration are observed for small streams. This behavior can be explained by considering the hyporheic residence time as the main controlling factor of DO consumption for the same biological activity. Path lines with long residence times show larger portions with low DO concentrations, which are less affected by the daily DO fluctuations of stream water.
Figure 7. Dimensionless stream BUR, BUR*, as a function of the dimensionless mean flow depth, Y B M * , for periodic and constant in-stream boundary conditions and RBN = 0.5 with (a) RDO = 0.2 and (b) RDO = 1.2, with RDO the ratio between in-stream DO (DO0) and total oxygen demand (L0C + L0N) and RBN the ratio between in-stream carbonaceous and nitrogenous BOD concentrations, RBN = L0C/L0N.
Figure 7. Dimensionless stream BUR, BUR*, as a function of the dimensionless mean flow depth, Y B M * , for periodic and constant in-stream boundary conditions and RBN = 0.5 with (a) RDO = 0.2 and (b) RDO = 1.2, with RDO the ratio between in-stream DO (DO0) and total oxygen demand (L0C + L0N) and RBN the ratio between in-stream carbonaceous and nitrogenous BOD concentrations, RBN = L0C/L0N.
Water 07 00398 g007
Conversely, large portions of path lines with short residence times show DO concentrations that vary temporally, depending on stream nutrient and DO concentrations. Thus, systems with predominantly long residence times, like streams with large Y B M * , show a similar global response for both boundary conditions; whereas the BUR* response of streams with small Y*BM shows a dependence on temporal variations of stream nutrient and DO concentrations. However, as shown in Figure 7a,b, the dependence on DO fluctuations in stream water is small, yet appreciable, thereby suggesting that the mean daily concentrations of DO are sufficient information for studying the global response of the hyporheic zone in most cases of practical interest.
BUR* shows the overall system response, which is certainly important for understanding the mean behavior of the system for different types of stream nutrient loads (Figure 7a,b). However, the lifetime of many aquatic organisms is strongly related to the local behavior of the system, which could be quite different from the overall response. In this context, Figure 8 shows a three-dimensional snapshot of simulated DO concentrations within the hyporheic zone in the case of sinusoidally-periodic (Figure 8a) and constant (Figure 8b) DO, LC and LN concentrations for RDO = 1.2 and RBN = 0.5.
Figure 8. Maps of dissolved oxygen concentrations, DO, within the hyporheic zone for (a) periodic and (b) constant in-stream concentrations of DO, CBOD and NBOD. In both cases, RDO = 1.2 and RBN = 0.5 (Test 5-A), with RDO the ratio between in-stream DO (DO0) and total oxygen demand (L0C + L0N) and RBN the ratio between in-stream carbonaceous and nitrogenous BOD concentrations, RBN = L0C/L0N.
Figure 8. Maps of dissolved oxygen concentrations, DO, within the hyporheic zone for (a) periodic and (b) constant in-stream concentrations of DO, CBOD and NBOD. In both cases, RDO = 1.2 and RBN = 0.5 (Test 5-A), with RDO the ratio between in-stream DO (DO0) and total oxygen demand (L0C + L0N) and RBN the ratio between in-stream carbonaceous and nitrogenous BOD concentrations, RBN = L0C/L0N.
Water 07 00398 g008
The variation in time of the two nutrients and DO in the stream water induce differentiated spatial and temporal distributions of DO concentrations within the hyporheic zone. The temporal changes of DO concentrations are not systematic within the domain, but they are relevant for short residence times and become negligible along the streamlines as the residence time increases. Thus, the hyporheic zone may show a dynamic volume of sediment, whose DO concentrations fluctuate following those of the stream water, and a stable volume, where DO concentrations are low and stationary. We should expect a spatial specialization of organisms dwelling in the hyporheic zone following the DO profile as those observed at the sediment-water interface by Malard et al. [41].

4. Conclusions

We presented an analytical model, which accounts for advection, diffusion and both carbonaceous and nitrogenous BODs, for predicting DO concentrations within the hyporheic zone. The application of the model on a gravel-bed pool-riffle reach shows that the proposed semi-analytical tool is able to predict the spatial distribution of measured DO concentrations. The overall agreement between predicted and observed DO concentrations is 58%. Our results show that hyporheic DO consumption varies as a function of bedform morphology, size and nutrient loads. Our model uses spatially-constant reaction rates, which probably may change spatially, and homogenous hydraulic conductivities. These two limitations may be important to detect micro-biogeochemical hotspots, which may be linked to localized slow velocities [65].
Our simulations suggest that river morphology plays an important role in controlling the BOD removal capacity of the hyporheic zone. However, the system responds differently at local and global scales, depending on the temporal variations of stream DO concentrations. The bedform scale response of the hyporheic zone expressed by BUR* is only minimally influenced by daily oscillations of stream water nutrient and DO concentrations. Thus, daily mean concentrations of these solutes can be used as an index for quantifying the hyporheic response at the bedform scale. Conversely, daily oscillations of stream nutrient concentrations may influence the spatial distribution and density of aquatic organisms locally. Finally, we show that for a given nutrient concentration, there is a stream size that maximizes the removal of LC and LN and, consequently, the consumption of stream dissolved oxygen.

Acknowledgments

This research was partially supported by the Deadwood River Project, U.S. Forest Service Award 009421-01, and by the European Communities 7th Framework Programme under Grant Agreement No. 603629-ENV-2013-6.2.1-Globaqua. Any opinions, conclusions, or recommendations expressed in this work are solely those of the authors and do not necessarily reflect the views of the supporting agencies.

Appendix

DO transport and transformation within the hyporheic zone is modeled solving Equation (2) for both constant and sinusoidally varying in-stream solute concentrations. The initial and boundary conditions are:
{ DO ( 0 , t ) = DO 0 , L C ( 0 , t ) = L 0 C a n d L N ( 0 , t ) = L N 0 DO ( τ,0 ) = L C ( τ,0 ) = L N ( τ,0 ) = 0 DO τ ( , t ) = L C τ ( , t ) = L N τ ( , t ) = 0
and:
{ DO ( 0 , t ) = DO 0 + DO A sin ( ω t ) , L C ( 0 , t ) = L 0 C + L A C sin ( ω L t ) , L N ( 0 , t ) = L 0 N + L A N sin ( ω N t ) DO ( τ,0 ) = L C ( τ,0 ) = L N ( τ,0 ) = 0 DO τ ( , t ) = L C τ ( , t ) = L N τ ( , t ) = 0
for constant and sinusoidally time varying DO concentration in stream water, respectively. In Equation (A1), L0C and L0N are, respectively, the constant concentrations of LC and LN in the stream water. On the other hand, in Equation (A2), LAC and LAN and DOA are the amplitude of the in-stream oscillation signal for LC, LN and DO, respectively, and ωL, ωN and ω are their periods of oscillation.
The solution of Equation (2) with initial and boundary conditions expressed in Equation (A1), which provides the contribution of one stream tube in DO transport, assumes the following form:
D O ( τ , t ) = D O lim [ 1 1 2 exp ( τ α ) e r f c ( τ + t 2 α t ) + e r f c ( τ t 2 α t ) ] + 1 2 ( D O 0 K L L 0 C K L K R K N L 0 N K N K R ) { exp [ τ 2 α ( 1 1 + 4 K R α ) ] e r f c ( τ t 1 + 4 K R α 2 α t ) + exp [ τ 2 α ( 1 + 1 + 4 K R α ) ] e r f c ( τ + t 1 + 4 K R α 2 α t ) } + L 0 C 2 { exp [ τ 2 α ( 1 1 + 4 K L α ) ] e r f c ( τ t 1 + 4 K L α 2 α t ) + exp [ τ 2 α ( 1 + 1 + 4 K L α ) ] e r f c ( τ + t 1 + 4 K L α 2 α t ) } + L 0 N 2 { exp [ τ 2 α ( 1 1 + 4 K N α ) ] e r f c ( τ t 1 + 4 K N α 2 α t ) + exp [ τ 2 α ( 1 + 1 + 4 K N α ) ] e r f c ( τ + t 1 + 4 K N α 2 α t ) }
In addition, the solution of Equation (2) with initial and boundary conditions expressed in Equation (A2), which, similarly to the previous case, provides the contribution of one stream tube in DO transport, assumes the following form:
D O ( τ , t ) = D O lim [ 1 1 2 exp ( τ α ) e r f c ( τ + t 2 α t ) + e r f c ( τ t 2 α t ) ] + 1 2 ( D O 0 K L L 0 C K L K R K N L 0 N K N K R ) { exp [ τ 2 α ( 1 1 + 4 K R α ) ] e r f c ( τ t 1 + 4 K R α 2 α t ) + exp [ τ 2 α ( 1 + 1 + 4 K R α ) ] e r f c ( τ + t 1 + 4 K R α 2 α t ) } + 1 2 ( D O A K L L A C K L K R K N L A N K N K R ) { exp ( τ 2 α + γ D O ) [ ( 1 F 1 ) sin ( ω D O t + η R ) G 1 cos ( ω D O t + η D O ) ] + exp ( τ 2 α γ D O ) [ ( 1 F 2 ) sin ( ω D O t η R ) G 2 cos ( ω D O t η R ) ] } + K L L 0 C 2 ( K L K R ) { exp [ τ 2 α ( 1 1 + 4 K L α ) ] e r f c ( τ t 1 + 4 K L α 2 α t ) + exp [ τ 2 α ( 1 + 1 + 4 K L α ) ] e r f c ( τ + t 1 + 4 K L α 2 α t ) } + K L L A C 2 ( K L K R ) { exp ( τ 2 α + γ L ) [ ( 1 F 3 ) sin ( ω L t + η L ) G 3 cos ( ω L t + η L ) ] + exp ( τ 2 α γ L ) [ ( 1 F 4 ) sin ( ω L t η L ) G 4 cos ( ω L t η L ) ] } + K N L 0 N 2 ( K N K R ) { exp [ τ 2 α ( 1 1 + 4 K N α ) ] e r f c ( τ t 1 + 4 K N α 2 α t ) + exp [ 2 τ α ( 1 + 1 + 4 K N α ) ] e r f c ( τ + t 1 + 4 K N α 2 α t ) } + K N L A N 2 ( K N K R ) { exp ( τ 2 α + γ N ) [ ( 1 F 5 ) sin ( ω N t + η N ) G 5 cos ( ω N t + η N ) ] + exp ( τ 2 α γ N ) [ ( 1 F 6 ) sin ( ω N t η N ) G 6 cos ( ω N t η N ) ] }
where Fi (i = 1 ÷ 6) and Gi (i = 1 ÷ 6) are the functions reported in Marzadri et al. [45] and:
r j = 1 α ( 4 α K j ) 2 16 α 2 + ω j 2 ;   θ j = a r c t g ( 4 ω j α 4 K j α + 1 ) ;   w i t h   j = D O , L , N γ j = r j τ cos ( θ j 2 ) ; η j = r j τ sin ( θ j 2 ) ;   w i t h   j = D O , L , N
When the stationary conditions are reached, the solution of Equation (A4) further simplifies to:
D O ( τ , t ) = D O lim [ 1 1 2 exp ( τ α ) e r f c ( τ + t 2 α t ) + e r f c ( τ t 2 α t ) ] + + 1 2 ( D O 0 K L L 0 C K L K R K N L 0 N K N K R ) { exp [ τ 2 α ( 1 1 + 4 K R α ) ] e r f c ( τ t 1 + 4 K R α 2 α t ) + + exp [ τ 2 α ( 1 + 1 + 4 K R α ) ] e r f c ( τ + t 1 + 4 K R α 2 α t ) } + + ( D O A K L L A C K L K R K N L A N K N K R ) exp ( τ 2 α γ D O ) sin ( ω D O t η R ) + K L L 0 C 2 ( K L K R ) { exp [ τ 2 α ( 1 1 + 4 K L α ) ] e r f c ( τ t 1 + 4 K L α 2 α t ) + + exp [ τ 2 α ( 1 + 1 + 4 K L α ) ] e r f c ( τ + t 1 + 4 K L α 2 α t ) } + + K L L A C K L K R exp ( τ 2 α γ L ) sin ( ω L t η L ) + K N L 0 N 2 ( K N K R ) { exp [ τ 2 α ( 1 1 + 4 K N α ) ] e r f c ( τ t 1 + 4 K N α 2 α t ) + + exp [ τ 2 α ( 1 + 1 + 4 K N α ) ] e r f c ( τ + t 1 + 4 K N α 2 α t ) } + + K N L A N K N K R exp ( τ 2 α γ N ) sin ( ω N t η N )
Carbonaceous (LC) and nitrogenous (LN) BOD, transport and transformation within the hyporheic zone are modeled solving the transport Equations (A7) and (A8) for both constant (A1) and sinusoidally (A2) varying in-stream solute concentrations:
L C t + L C τ = α   2 L C K L L C
L N t + L N τ = α   2 L N K N L N
The solution of Equations (A7) and (A8) with initial and boundary conditions expressed in Equation (A1), which provides the contribution of one stream tube in carbonaceous and nitrogenous BOD transport, respectively, assumes the following forms:
L C ( τ , t ) = L 0 C 2 { exp [ τ 2 α ( 1 1 + 4 K L α ) ] e r f c ( τ t 1 + 4 K L α 2 α t ) + exp [ τ 2 α ( 1 + 1 + 4 K L α ) ] e r f c ( τ + t 1 + 4 K L α 2 α t ) }
L N ( τ , t ) = L 0 N 2 { exp [ τ 2 α ( 1 1 + 4 K N α ) ] e r f c ( τ t 1 + 4 K N α 2 α t ) + exp [ τ 2 α ( 1 + 1 + 4 K N α ) ] e r f c ( τ + t 1 + 4 K N α 2 α t ) }
In addition, the solution of Equations (A7) and (A8) with initial and boundary conditions expressed in Equation (A2), which, similarly to the previous case, provides the contribution of one stream tube in carbonaceous and nitrogenous BOD transport, respectively, assumes the following forms:
L C ( τ , t ) = L 0 C 2 { exp [ τ 2 α ( 1 1 + 4 K L α ) ] e r f c ( τ t 1 + 4 K L α 2 α t ) + exp [ τ 2 α ( 1 + 1 + 4 K L α ) ] e r f c ( τ + t 1 + 4 K L α 2 α t ) } + L A C 2 { exp ( τ 2 α + γ L ) [ ( 1 F 3 ) sin ( ω L t + η L ) G 3 cos ( ω L t + η L ) ] + exp ( τ 2 α γ L ) [ ( 1 F 4 ) sin ( ω L t η L ) G 4 cos ( ω L t η L ) ] }
and:
L N ( τ , t ) = L 0 N 2 { exp [ τ 2 α ( 1 1 + 4 K N α ) ] e r f c ( τ t 1 + 4 K N α 2 α t ) + exp [ 2 τ α ( 1 + 1 + 4 K N α ) ] e r f c ( τ + t 1 + 4 K N α 2 α t ) } + L A N 2 { exp ( τ 2 α + γ N ) [ ( 1 F 5 ) sin ( ω N t + η N ) G 5 cos ( ω N t + η N ) ] + exp ( τ 2 α γ N ) [ ( 1 F 6 ) sin ( ω N t η N ) G 6 cos ( ω N t η N ) ] }
where Fi (i = 1 ÷ 6) and Gi (i = 1 ÷ 6) are the functions reported in Marzadri et al. [45] and in Equation (A5). Under steady-state conditions, the solution of Equations (A11) and (A12) simplify to:
L C ( τ , t ) = L 0 C 2 { exp [ τ 2 α ( 1 1 + 4 K L α ) ] e r f c ( τ t 1 + 4 K L α 2 α t ) + + exp [ τ 2 α ( 1 + 1 + 4 K L α ) ] e r f c ( τ + t 1 + 4 K L α 2 α t ) } + + L A C exp ( τ 2 α γ L ) sin ( ω L t η L )
L N ( τ , t ) = L 0 N 2 { exp [ τ 2 α ( 1 1 + 4 K N α ) ] e r f c ( τ t 1 + 4 K N α 2 α t ) + + exp [ τ 2 α ( 1 + 1 + 4 K N α ) ] e r f c ( τ + t 1 + 4 K N α 2 α t ) } + + L A N exp ( τ 2 α γ N ) sin ( ω N t η N )

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Galloway, J.N.; Dentener, F.J.; Capone, D.G.; Boyer, E.W.; Howarth, R.W.; Seitzinger, S.P.; Asner, G.P.; Cleveland, C.C.; Green, P.A.; Holland, E.A.; et al. Nitrogen cycles: Past, present, and future. Biogeochemestry 2004, 70, 153–226. [Google Scholar] [CrossRef]
  2. Galloway, J.N.; Townsend, A.R.; Erisman, J.W.; Bekunda, M.; Cai, Z.; Freney, J.R.; Martinelli, L.A.; Seitzinger, S.P.; Sutton, M.A. Transformation of the nitrogen cycle: Recent trends, questions, and potential solutions. Science 2008, 320, 889–892. [Google Scholar] [CrossRef] [PubMed]
  3. Jelks, H.; Walsh, S.; Burkhead, N.; Contreras-Balderas, S.; Diaz-Pardo, E.; Hendrickson, D.; Lyons, J.; Mandrak, N.; McCormick, F.; Nelson, J. Conservation status of imperiled north american freshwater and diadromous fishes. Fisheries 2008, 33, 372–407. [Google Scholar] [CrossRef]
  4. McMahon, P.B.; Tindall, J.A.; Collins, J.A.; Lull, K.J.; Nuttle, J.R. Hydrologic and geochemical effect on oxygen uptake in bottom sediments of an effluent-dominated river. Water Resour. Res. 1995, 31, 2561–2569. [Google Scholar] [CrossRef]
  5. Paulsen, S.G.; Mayio, A.; Peck, D.; Stoddard, J.; Tarquinio, E.; Holdsworth, S.; Van Sickle, J.; Yuan, L.; Hawkins, C.; Herlihy, A.T. Condition of stream ecosystems in the us: An overview of the first national assessment. J. Inf. 2008, 27, 812–821. [Google Scholar]
  6. Cooper, A.B. Activities of benthic nitrifiers in streams and their role in oxygen consumption. Microb. Ecol. 1984, 10, 317–334. [Google Scholar] [CrossRef] [PubMed]
  7. Štambuk-Giljanović, N. The pollution load by nitrogen and phosphorus in the cetina river. Water Air Soil Pollut. 2010, 211, 49–60. [Google Scholar] [CrossRef]
  8. Pollock, M.S.; Clarke, L.M.J.; Dubé, M.G. The effects of hypoxia on fishes: From ecological relevance to physiological effects. Environ. Rev. 2007, 15, 1–14. [Google Scholar] [CrossRef]
  9. Yu, F.X.; Adrian, D.D.; Singh, V.P. Modeling river quality by the superposition method. J. Environ. Syst. 1991, 20, 1–16. [Google Scholar]
  10. Rusjan, S.; Mikoš, M. Seasonal variability of diurnal in-stream nitrate concentration oscillations under hydrologically stable conditions. Biogeochemestry 2010, 97, 123–140. [Google Scholar] [CrossRef]
  11. Gariglio, F.; Tonina, D.; Luce, C.H. Spatio-temporal variability of hyporheic exchange through a pool-riffle-pool sequence. Water Resour. Res. 2013, 49, 7185–7204. [Google Scholar] [CrossRef]
  12. Marzadri, A.; Tonina, D.; Bellin, A. Effects of stream morphodynamics on hyporheic zone thermal regime. Water Resour. Res. 2013, 49, 2287–2302. [Google Scholar] [CrossRef]
  13. Alshawabkeh, A.; Adrian, D.D. Analytical water quality model for a sinusoidally varying bod discharge concentration. Water Resour. Res. 1997, 31, 1207–1215. [Google Scholar]
  14. Higashino, M.; Gantzer, C.; Stefan, H.G. Unsteady diffusional mass transfer at the sediment/water interface: Theory and significance for sod measurement. Water Res. 2004, 38, 1–12. [Google Scholar] [CrossRef] [PubMed]
  15. Jorgensen, B.; Revsbech, N. Diffusive boundary layers and the oxygen uptake of sediments and detritus. Limnol. Oceanogr. 1985, 30, 111–122. [Google Scholar] [CrossRef]
  16. Rutherford, I.C.; Boyle, J.D.; Elliott, A.H.; Hatherell, T.V.J.; Chiu, T.W. Modeling benthic oxygen uptake by pumping. J. Environ. Eng. 1995, 121, 84–95. [Google Scholar] [CrossRef]
  17. Buffington, J.M.; Tonina, D. Hyporheic exchange in mountain rivers ii: Effects of channel morphology on mechanics, scales, and rates of exchange. Geogr. Compass 2009, 3, 1038–1062. [Google Scholar] [CrossRef]
  18. Elliott, A.; Brooks, N.H. Transfer of nonsorbing solutes to a streambed with bed forms: Theory. Water Resour. Res. 1997, 33, 123–136. [Google Scholar] [CrossRef]
  19. Elliott, A.; Brooks, N.H. Transfer of nonsorbing solutes to a streambed with bed forms: Laboratory experiments. Water Resour. Res. 1997, 33, 137–151. [Google Scholar] [CrossRef]
  20. Tonina, D.; Buffington, J.M. Hyporheic exchange in mountain rivers i: Mechanics and environmental effects. Geogr. Compass 2009, 3, 1063–1086. [Google Scholar] [CrossRef]
  21. Zarnetske, J.P.; Haggerty, R.; Wondzell, S.M.; Baker, M.A. Dynamics of nitrate production and removal as a function of residence time in the hyporheic zone. J. Geophys. Res. 2011, 116, G01025. [Google Scholar]
  22. Zarnetske, J.P.; Haggerty, R.; Wondzell, S.M.; Bokil, V.A.; González-Pinzón, R. Coupled transport and reaction kinetics control the nitrate source-sink function of hyporheic zones. Water Resour. Res. 2012, 48, W11508. [Google Scholar]
  23. Triska, F.J.; Duff, J.H.; Avanzino, R.J. The role of water exchange between a stream channel and its hyporheic zone in nitrogen cycling at the terrestrial-aquatic interface. Hydrobiologia 1993, 251, 167–184. [Google Scholar] [CrossRef]
  24. Triska, F.J.; Duff, J.H.; Avanzino, R.J. Patterns of hydrological exchange and nutrient transformation in the hyporheic zone of a gravel bottom stream: Examining terrestrial-aquatic linkages. Freshw. Biol. 1993, 29, 259–274. [Google Scholar] [CrossRef]
  25. Grimaldi, C.; Chaplot, V. Nitrate depletion during within-stream transport: Effects of exchange processes between streamwater, the hyporheic and riparian zones. Water Air Soil Pollut. 2000, 124, 95–112. [Google Scholar] [CrossRef]
  26. Hlaváčová, E.; Rulík, M.; Čáp, L. Anaerobic microbial metabolism in hyporheic sediment of gravel bar in a small lowland stream. River Res. Appl. 2005, 21, 1003–1011. [Google Scholar] [CrossRef]
  27. Kasahara, T.; Hill, A.R. Effects of riffle–step restoration on hyporheic zone chemistry in n-rich lowland streams. Can. J. Fish. Aquat. Sci. 2006, 63, 120–133. [Google Scholar] [CrossRef]
  28. Pinay, G.; O’Keefe, T.C.; Edwards, R.T.; Naiman, R.J. Nitrate removal in the hyporheic zone of a salmon river in alaska. River Res. Appl. 2009, 25, 367–375. [Google Scholar] [CrossRef]
  29. Krause, S.; Heathwaite, L.; Binley, A.; Keenan, P. Nitrate concentration changes at the groundwater-surface water interface of a small cumbrian river. Hydrol. Process. 2009, 23, 2195–2211. [Google Scholar] [CrossRef]
  30. Fernald, A.G.; Landers, D.H.; Wigington, P.J. Water quality changes in hyporheic fow paths between a large gravel bed river and off-channel alcoves in oregon, USA. River Res. Appl. 2006, 22, 1111–1124. [Google Scholar] [CrossRef]
  31. Marzadri, A.; Tonina, D.; Bellin, A. Morphodynamic controls on redox conditions and on nitrogen dynamics within the hyporheic zone: Application to gravel bed rivers with alternate-bar morphology. J. Geophys. Res. 2012, 117. [Google Scholar] [CrossRef]
  32. Briggs, M.A.; Lautz, L.K.; Hare, D.K. Residence time control on hot moments on net nitrate and uptake in hyporheic zone. Hydrol. Process. 2014, 28, 3741–3751. [Google Scholar] [CrossRef]
  33. Marzadri, A.; Tonina, D.; Bellin, A.; Tank, J.L. A hydrologic model demonstrates nitrous oxide emissions depend on streambed morphology. Geophys. Res. Lett. 2014, 41, 5484–5491. [Google Scholar] [CrossRef]
  34. Trauth, N.; Schmidt, J.C.; Vieweg, M.; Maier, U.; Fleckenstein, J.H. Hyporheic transport and biogeochemical reactions in pool-riffle systems under varying ambient groundwater flow conditions. J. Geophys. Res. Biogeosci. 2014, 119, 910–928. [Google Scholar] [CrossRef]
  35. Marzadri, A.; Tonina, D.; Bellin, A. A semianalytical three-dimensional process-based model for hyporheic nitrogen dynamics in gravel bed rivers. Water Resour. Res. 2011, 47, W11518. [Google Scholar]
  36. Bardini, L.; Boano, F.; Cardenas, B.M.; Revelli, R.; Ridolfi, L. Nutrient cycling in bedform induced hyporheic zones. Geochim. Cosmochim. Acta 2012, 84, 47–61. [Google Scholar] [CrossRef]
  37. Hester, E.T.; Young, K.I.; Widdowson, M.A. Controls on mixing-dependent denitrification in hyporheic zones induced by riverbed dunes: A steady state modeling study. Water Resour. Res. 2014, 50, 9048–9066. [Google Scholar] [CrossRef]
  38. Boano, F.; Demaria, A.; Revelli, R.; Ridolfi, L. Biogeochemical zonation due to intrameander hyporheic flow. Water Resour. Res. 2010, 46, W02511. [Google Scholar]
  39. Montgomery, D.R.; Buffington, J.M. Channel-reach morphology in mountain drainage basins. Geol. Soc. Am. Bull. 1997, 109, 596–611. [Google Scholar] [CrossRef]
  40. Tonina, D.; Buffington, J.M. A three-dimensional model for analyzing the effects of salmon redds on hyporheic exchange and egg pocket habitat. Can. J. Fish. Aquat. Sci. 2009, 66, 2157–2173. [Google Scholar] [CrossRef]
  41. Malard, F.; Tockner, K.; Dole-Olivier, M.-J.; Ward, J.V. A landscape perspective of surface-subsurface hydrological exchanges in river corridors. Freshw. Biol. 2002, 47, 621–640. [Google Scholar] [CrossRef]
  42. Rouch, R. Sur la répartition spatiale des crustacés dans le sous-écoulement d'un ruisseau des pyrénées. Ann. Limnol. 1988, 24, 213–234. [Google Scholar] [CrossRef]
  43. Freeze, R.A.; Cherry, J.A. Groundwater; Prentice Hall: Englewood Cliffs, NJ, USA, 1979. [Google Scholar]
  44. Cardenas, M.B.; Wilson, J.L. Hydrodynamics of coupled flow above and below a sediment–water interface with triangular bedforms. Adv. Water Resour. 2007, 30, 301–313. [Google Scholar] [CrossRef]
  45. Marzadri, A.; Tonina, D.; Bellin, A.; Vignoli, G.; Tubino, M. Effects of bar topography on hyporheic flow in gravel-bed rivers. Water Resour. Res. 2010, 46, W07531. [Google Scholar]
  46. Colombini, M.; Seminara, G.; Tubino, M. Finite-amplitude alternate bars. J. Fluid Mech. 1987, 181, 213–232. [Google Scholar] [CrossRef]
  47. Tonina, D.; Bellin, A. Effects of pore-scale dispersion, degree of heterogeneity, sampling size, and source volume on the concentration moments of conservative solutes in heterogeneous formations. Adv. Water Resour. 2007, 31, 339–354. [Google Scholar] [CrossRef]
  48. Marzadri, A.; Tonina, D.; Bellin, A. Effects of Hyporheic Fluxes on Streambed Pore Water Temperature. In Proceedings of the 10th International Conference on Hydroinformatics HIC 2012, Hamburg, Germany, 14–18 July, 2012; International Association of Hydro-Environment Engineering and Research (IAHR): Hamburg, Germany.
  49. Dagan, G.; Cvetkovic, V.; Shapiro, A.M. A solute flux approach to transport in heterogeneous formations: 1. The general framework. Water Resour. Res. 1992, 28, 1369–1376. [Google Scholar] [CrossRef]
  50. Marzadri, A.; Tonina, D.; Bellin, A. Quantifying the importance of daily stream water temperature fluctuations on the hyporheic thermal regime: Implication for dissolved oxygen dynamics. J. Hydrol. 2013, 507, 241–248. [Google Scholar] [CrossRef]
  51. Fogler, H.S. Elements of Chemical Reaction Engineering; Prentice-Hall: Englewood Cliffs, NJ, USA, 1999; p. 446. [Google Scholar]
  52. Marion, A.; Packman, A.I.; Zaramella, M.; Bottacin-Busolin, A. Hyporheic flows in stratified beds. Water Resour. Res. 2008, 44, W09433. [Google Scholar]
  53. Marion, A.; Zaramella, M.; Bottacin-Busolin, A. Solute transport in rivers with multiple storage zones: The stir model. Water Resour. Res. 2008, 44, W10406. [Google Scholar]
  54. Tonina, D. Surface water and streambed sediment interaction: The hyporheic exchange. In Fluid Mechanics of Environmental Interfaces; Gualtieri, C., Mihailović, D.T., Eds.; CRC Press, Taylor & Francis Group: London, UK, 2012; pp. 255–294. [Google Scholar]
  55. Tonina, D.; Buffington, J.M. Effects of stream discharge, alluvial depth and bar amplitude on hyporheic flow in pool-riffle channels. Water Resour. Res. 2011, 47, W08508. [Google Scholar]
  56. Beaulieu, J.J.; Tank, J.L.; Hamilton, S.K.; Wollheim, W.M.; Hall, R.O.J.; Mulholland, P.J.; Peterson, B.J.; Ashkenas, L.R.; Cooper, L.W.; Dahm, C.N.; et al. Nitrous oxide emission from denitrification in stream and river networks. Proc. Natl. Acad. Sci. USA 2011, 108, 214–219. [Google Scholar] [CrossRef] [PubMed]
  57. Benjankar, R.; Burke, M.P.; Yager, E.M.; Tonina, D.; Egger, G.; Rood, S.B.; Merz, N. Development of a spatially-distributed hydroecological model to simulate cottonwood seedling recruitment along rivers. J. Environ. Manag. 2014, 145, 277–288. [Google Scholar] [CrossRef]
  58. Congalton, R.G.; Green, K. Assessing the Accuracy of Remotely Sensed Data: Principles and Practices; Lewis Publishers: Boca Raton, FL, USA, 2008; p. 200. [Google Scholar]
  59. Landis, J.R.; Koch, G.G. The measurement of observer agreement for categorical data. Biometrics 1977, 33, 159–174. [Google Scholar] [CrossRef] [PubMed]
  60. Benjankar, R.; Glenn, N.F.; Egger, G.; Jorde, K.; Goodwin, P. Comparison of field-observed and simulated map output from a dynamic floodplain vegetation model using remote sensing and gis techniques. GIScience Remote Sens. 2010, 47, 480–497. [Google Scholar] [CrossRef]
  61. Colombini, M. Revisiting the linear theory of sand dune formation. J. Hydraul. Eng. 2004, 502, 1–16. [Google Scholar]
  62. De Silva, A.M.F.; Zhang, Y. On the steepness of dunes and determination of alluvial stream friction factor. In Proceedings of the XXVIII IAHR Congress, Graz, Austria, 22–27 August 1999; p. 7.
  63. Van Rijn, L.C. Sediment transport, part iii: Bed forms and alluvial roughness. J. Hydraul. Eng. 1984, 110, 1733–1754. [Google Scholar] [CrossRef]
  64. Benjankar, R.; Tonina, D.; McKean, J.A. One-dimensional and two-dimensional hydrodynamic modeling derived flow properties: Impacts on aquatic habitat quality predictions. Earth Surf. Process. Landf. 2014. [Google Scholar] [CrossRef]
  65. Harvey, J.W.; Böhlke, J.K.; Voytek, M.A.; Scott, D.; Tobias, C.R. Hyporheic zone denitrification: Controls on effective reaction depth and contribution to whole-stream mass balance. Water Resour. Res. 2013, 49, 6298–6316. [Google Scholar] [CrossRef]

Share and Cite

MDPI and ACS Style

Tonina, D.; Marzadri, A.; Bellin, A. Benthic Uptake Rate due to Hyporheic Exchange: The Effects of Streambed Morphology for Constant and Sinusoidally Varying Nutrient Loads. Water 2015, 7, 398-419. https://doi.org/10.3390/w7020398

AMA Style

Tonina D, Marzadri A, Bellin A. Benthic Uptake Rate due to Hyporheic Exchange: The Effects of Streambed Morphology for Constant and Sinusoidally Varying Nutrient Loads. Water. 2015; 7(2):398-419. https://doi.org/10.3390/w7020398

Chicago/Turabian Style

Tonina, Daniele, Alessandra Marzadri, and Alberto Bellin. 2015. "Benthic Uptake Rate due to Hyporheic Exchange: The Effects of Streambed Morphology for Constant and Sinusoidally Varying Nutrient Loads" Water 7, no. 2: 398-419. https://doi.org/10.3390/w7020398

APA Style

Tonina, D., Marzadri, A., & Bellin, A. (2015). Benthic Uptake Rate due to Hyporheic Exchange: The Effects of Streambed Morphology for Constant and Sinusoidally Varying Nutrient Loads. Water, 7(2), 398-419. https://doi.org/10.3390/w7020398

Article Metrics

Back to TopTop