Modeling Benthic Versus Hyporheic Nutrient Uptake in Unshaded Streams With Varying Substrates

Assessments of riverine ecosystem health and water quality require knowledge of how headwater streams transport and transform nutrients. Estimates of nutrient demand at the watershed scale are commonly inferred from reach-scale solute injections, which are typically reported as uptake velocities (v(f)). Multiple interacting processes control v(f), making it challenging to predict how v(f) responds to physical changes in the stream. In this study, we link v(f) to a continuous time random walk model to quantify how v(f) is controlled by in-stream (velocity, dispersion, and benthic reaction) and hyporheic processes (exchange rate, residence times, and hyporheic reaction). We fit the model to conservative (NaCl) and nitrate (NO3--N) pulse tracer injections in unshaded replicate streams at the Notre Dame Linked Experimental Ecosystem Facility, which differed only in substrate size and distribution. Experiments were conducted over the first 25days of biofilm colonization to examine how the interaction between substrate type and biofilm growth influenced modeled processes and v(f). Model fits of benthic reaction rates were approximate to 8x greater than hyporheic reaction rates for all experiments and did not vary with substrate type or over time. High benthic reactivity was associated with filamentous green algae coverage on the streambed, which dominated total algal biomass. Finally, v(f) was most sensitive to benthic reaction rate and stream velocity, and sensitivity varied with stream conditions due to its nonlinear dependence on all modeled processes. Together, these results demonstrate how reach-scale nutrient demand reflects the relative contributions of biotic and abiotic processes in the benthic layer and the hyporheic zone.