Hyporheic exchange induced by periodic river fluctuations leads to important biogeochemical processes, particularly nitrogen cycling, in riparian zones (RZs) where chemically distinct surface water and groundwater mix. We developed a two‐dimensional coupled flow, reactive transport model to study the role of bank storage induced by river fluctuations on removing river‐borne nitrate. Sensitivity analyses were conducted to quantify the effects of river amplitude, sediment hydraulic conductivity and dispersivity, and ambient groundwater flow on nitrate removal rate. The simulations showed that nitrification occurred in the shallower zone adjacent to the bank where oxic river water and groundwater interacted while denitrification occurred deeper into the aquifer and in the riverbed sediments where oxygen was depleted. River fluctuations greatly increased the amount of nitrate being removed; the nitrate removal rate increased as river amplitude increased. Similarly, increasing hydraulic conductivity increased overall nitrate removal since it expanded the denitrifying zone but decreased efficiency. In contrast, increasing sediment dispersivity increased the removal efficiency of nitrate because it promoted mixing between electron acceptors and donors. The presence and direction of ambient groundwater flow had a significant impact on nitrate removal rate when compared to neutral conditions. A losing river showed a larger nitrate removal rate, whereas a gaining river showed a smaller nitrate removal rate. Our results demonstrated that daily river fluctuations created denitrification hot spots within the RZ that would not otherwise exist under naturally neutral or gaining conditions.

[1] This paper presents the formulation and calibration of the temperature portion of a two-zone temperature and solute (TZTS) model which separates transient storage into surface (STS) and subsurface transient storage (HTS) zones. The inclusion of temperature required the TZTS model formulation to differ somewhat from past transient storage models in order to accommodate terms associated with heat transfer. These include surface heat fluxes in the main channel (MC) and STS, heat and mass exchange between the STS and MC, heat and mass exchange between the HTS and MC, and heat exchange due to bed and deeper ground conduction. To estimate the additional parameters associated with a two-zone model, a data collection effort was conducted to provide temperature time series within each zone. Both single-objective and multiobjective calibration algorithms were then linked to the TZTS model to assist in parameter estimation. Single-objective calibrations based on MC temperatures at two different locations along the study reach provided reasonable predictions in the MC and STS. The HTS temperatures, however, were typically poorly estimated. The two-objective calibration using MC temperatures simultaneously at two locations illustrated that the TZTS model accurately predicts temperatures observed in MC, STS, and HTS zones, including those not used in the calibration. These results suggest that multiple data sets representing different characteristics of the system should be used when calibrating complex in-stream models.

Abstract. We investigated how absorption of sunlight by chromophoric dissolved organic matter (CDOM) controls the degradation and export of DOM from Imnavait Creek, a beaded stream in the Alaskan Arctic. We measured concentrations of dissolved organic carbon (DOC), as well as concentrations and characteristics of CDOM and fluorescent dissolved organic matter (FDOM), during ice-free periods of 2011-2012 in the pools of Imnavait Creek and in soil waters draining to the creek. Spatial and temporal patterns in CDOM and FDOM in Imnavait Creek were analyzed in conjunction with measures of DOM degradation by sunlight and bacteria and assessments of hydrologic residence times and in situ UV exposure. CDOM was the dominant light attenuating constituent in the UV and visible portion of the solar spectrum, with high attenuation coefficients ranging from 86 ± 12 m −1 at 305 nm to 3 ± 1 m −1 in the photosynthetically active region (PAR). High rates of light absorption and thus light attenuation by CDOM contributed to thermal stratification in the majority of pools in Imnavait Creek under low-flow conditions. In turn, thermal stratification increased the residence time of water and DOM, and resulted in a separation of water masses distinguished by contrasting UV exposure (i.e., UV attenuation by CDOM with depth resulted in bottom waters receiving less UV than surface waters). When the pools in Imnavait Creek were stratified, DOM in the pool bottom water closely resembled soil water DOM in character, while the concentration and character of DOM in surface water was reproduced by experimental photo-degradation of bottom water. These results, in combination with water column rates of DOM degradation by sunlight and bacteria, suggest that photo-degradation is the dominant process controlling DOM fate and export in Imnavait Creek. A conceptual model is presented showing how CDOM amount and lability interact with incident UV light and water residence time to determine whether photo-degradation is "light-limited" or "substrate-limited". We suggest that degradation of DOM in CDOM-rich streams or ponds similar to Imnavait is typically light-limited under most flow conditions. Thus, export of DOM from this stream will be less under conditions that increase the light available for DOM photo-degradation (i.e., low flows, sunny days).

Abstract. Beaver dams affect hydrologic processes, channel complexity, and stream temperature in part by inundating riparian areas, influencing groundwater–surface water interactions, and changing fluvial processes within stream systems. We explored the impacts of beaver dams on hydrologic and temperature regimes at different spatial and temporal scales within a mountain stream in northern Utah over a 3-year period spanning pre- and post-beaver colonization. Using continuous stream discharge, stream temperature, synoptic tracer experiments, and groundwater elevation measurements, we documented pre-beaver conditions in the first year of the study. In the second year, we captured the initial effects of three beaver dams, while the third year included the effects of ten dams. After beaver colonization, reach-scale (~ 750 m in length) discharge observations showed a shift from slightly losing to gaining. However, at the smaller sub-reach scale (ranging from 56 to 185 m in length), the discharge gains and losses increased in variability due to more complex flow pathways with beaver dams forcing overland flow, increasing surface and subsurface storage, and increasing groundwater elevations. At the reach scale, temperatures were found to increase by 0.38 °C (3.8 %), which in part is explained by a 230 % increase in mean reach residence time. At the smallest, beaver dam scale (including upstream ponded area, beaver dam structure, and immediate downstream section), there were notable increases in the thermal heterogeneity where warmer and cooler niches were created. Through the quantification of hydrologic and thermal changes at different spatial and temporal scales, we document increased variability during post-beaver colonization and highlight the need to understand the impacts of beaver dams on stream ecosystems and their potential role in stream restoration.

Abstract. We investigated how absorption of sunlight by chromophoric dissolved organic matter (CDOM) controls the degradation and export of DOM from Imnavait Creek, a beaded stream in the Alaskan Arctic. We measured concentrations of dissolved organic carbon (DOC), as well as concentrations and characteristics of CDOM and fluorescent dissolved organic matter (FDOM), during ice-free periods of 2011–2012 in the pools of Imnavait Creek and in soil waters draining to the creek. Spatial and temporal patterns in CDOM and FDOM in Imnavait Creek were analyzed in conjunction with measures of DOM degradation by sunlight and bacteria and assessments of hydrologic residence times and in situ UV exposure. CDOM was the dominant light attenuating constituent in the UV and visible portion of the solar spectrum, with high attenuation coefficients ranging from 86 ± 12 m−1 at 305 nm to 3 ± 1 m−1 in the photosynthetically active region (PAR). High rates of light absorption and thus light attenuation by CDOM contributed to thermal stratification in the majority of pools in Imnavait Creek under low-flow conditions. In turn, thermal stratification increased the residence time of water and DOM, and resulted in a separation of water masses distinguished by contrasting UV exposure (i.e., UV attenuation by CDOM with depth resulted in bottom waters receiving less UV than surface waters). When the pools in Imnavait Creek were stratified, DOM in the pool bottom water closely resembled soil water DOM in character, while the concentration and character of DOM in surface water was reproduced by experimental photo-degradation of bottom water. These results, in combination with water column rates of DOM degradation by sunlight and bacteria, suggest that photo-degradation is the dominant process controlling DOM fate and export in Imnavait Creek. A conceptual model is presented showing how CDOM amount and lability interact with incident UV light and water residence time to determine whether photo-degradation is "light-limited" or "substrate-limited". We suggest that degradation, and thus export, of DOM in CDOM-rich streams or ponds similar to Imnavait is typically light-limited under most flow conditions.

Groundwater flow regimes in the seasonally thawed soils in areas of continuous permafrost are relatively unknown despite their potential role in delivering water, carbon, and nutrients to streams. Using numerical groundwater flow models informed by observations from a headwater catchment in arctic Alaska, United States, we identify several mechanisms that result in substantial surface‐subsurface water exchanges across the land surface during downslope transport and create a primary control on dissolved organic carbon loading to streams and rivers. The models indicate that surface water flowing downslope has a substantial groundwater component due to rapid surface‐subsurface exchanges across a range of hydrologic states, from unsaturated to flooded. Field‐based measurements corroborate the high groundwater contributions, and river dissolved organic carbon concentrations are similar to that of groundwater across large discharge ranges. The persistence of these groundwater contributions in arctic watersheds will influence carbon export to rivers as thaw depth increases in a warmer climate.

[1] This paper presents the multiobjective calibration results for temperature and solute from a two-zone temperature and solute (TZTS) model which separates transient storage into surface (STS) and subsurface (HTS) transient storage components. This model contains terms associated with surface heat fluxes in the main channel (MC) and STS, heat and mass exchange between the STS and MC, heat and mass exchange between the HTS and MC, and heat exchange due to bed and deeper ground conduction. To estimate the additional parameters associated with a multiple-zone model, a data collection effort was conducted to provide temperature time series and solute tracer curves representing the movement of heat and/or solute through each zone. A multiobjective calibration algorithm was linked to the TZTS model to assist in parameter estimation and to provide information about parameter uncertainty and tradeoffs associated with matching different combinations of observations (e.g., solute and/or temperature data gathered in various zones). Results generated from three different combinations of calibration data illustrated that the two-zone model accurately reproduces temperatures and tracer concentrations observed in different zones when considered independently. However, there were many parameter sets that resulted in objectively indistinguishable results. When tracer and temperature observations were considered simultaneously in model calibration, the simplistic representation of the surface and subsurface zones did not adequately reproduce both observation types in each zone. If the uncertainty in model parameters and the data are taken into account, however, the results of the study suggest that it is plausible to use temperature and tracer information simultaneously to better inform transient storage modeling approaches.

Abstract. Beaver dams affect hydrologic processes, channel complexity, and stream temperature by increasing inundated areas and influencing groundwater-surface water interactions. We explored the impacts of beaver dams on hydrologic and temperature regimes at different spatial and temporal scales within a mountain stream in northern Utah over a three-year period spanning pre- and post-beaver colonization. Using continuous stream discharge, stream temperature, synoptic tracer experiments, and groundwater elevation measurements we documented pre-beaver conditions in the first year of the study. In the second year, we captured the initial effects of three beaver dams, while the third year included the effects of ten dams. After beaver colonization, reach scale discharge observations showed a shift from slightly losing to gaining. However, at the smaller sub-reach scale, the discharge gains and losses increased in variability due to more complex flow pathways with beaver dams forcing overland flow and increasing surface and subsurface storage. At the reach scale, temperatures were found to increase by 0.38 °C (3.8%), which in part is explained by a 230% increase in mean reach residence time. At the smallest, beaver dam scale, there were notable increases in the thermal heterogeneity where warmer and cooler niches were created. Through the quantification of hydrologic and thermal changes at different spatial and temporal scales, we document increased variability during post-beaver colonization and highlight the need to understand the impacts of beaver dams on stream ecosystems and their potential role in stream restoration.