•Implementation of a mechanistic, dynamic model of network expansion and contraction.•Importance of geologic setting and hydrologic forcing change through a water year.•Network expansion is ...insensitive to hydrologic forcing under wet conditions.•Geologic setting matters most under low and moderate discharge conditions.•Prediction of channel network dynamics may inform management of river corridors.
Headwater stream networks expand and contract in response to changes in stream discharge. The changes in the extent of the stream network are also controlled by geologic or geomorphic setting – some reaches go dry even under relatively wet conditions, other reaches remain flowing under relatively dry conditions. While such patterns are well recognized, we currently lack tools to predict the extent of the stream network and the times and locations where the network is dry within large river networks. Here, we develop a perceptual model of the river corridor in a headwater mountainous catchment, translate this into a reduced-complexity mechanistic model, and implement the model to examine connectivity and network extent over an entire water year. Our model agreed reasonably well with our observations, showing that the extent and connectivity of the river network was most sensitive to hydrologic forcing under the lowest discharges (Qgauge < 1 L s−1), that at intermediate discharges (1 L s−1 < Qgauge < 10 L s−1) the extent of the network changed dramatically with changes in discharge, and that under wet conditions (Qgauge > 10 L s−1) the extent of the network was relatively insensitive to hydrologic forcing and was instead determined by the network topology. We do not expect that the specific thresholds observed in this study would be transferable to other catchments with different geology, topology, or hydrologic forcing. However, we expect that the general pattern should be robust: the dominant controls will shift from hydrologic forcing to geologic setting as discharge increases. Furthermore, our method is readily transferable as the model can be applied with minimal data requirements (a single stream gauge, a digital terrain model, and estimates of hydrogeologic properties) to estimate flow duration or connectivity along the river corridor in unstudied catchments. As the available information increases, the model could be better calibrated to match site-specific observations of network extent, locations of dry reaches, or solute break through curves as demonstrated in this study. Based on the low initial data requirements and ability to later tune the model to a specific site, we suggest example applications of this parsimonious model that may prove useful to both researchers and managers.
Lakes, reservoirs, and other ponded waters are ubiquitous features of the aquatic landscape, yet their cumulative role in nitrogen removal in large river basins is often unclear. Here we use ...predictive modeling, together with comprehensive river water quality, land use, and hydrography datasets, to examine and explain the influences of more than 18,000 ponded waters on nitrogen removal through river networks of the Northeastern United States. Thresholds in pond density where ponded waters become important features to regional nitrogen removal are identified and shown to vary according to a ponded waters' relative size, network position, and degree of connectivity to the river network, which suggests worldwide importance of these new metrics. Consideration of the interacting physical and biological factors, along with thresholds in connectivity, reveal where, why, and how much ponded waters function differently than streams in removing nitrogen, what regional water quality outcomes may result, and in what capacity management strategies could most effectively achieve desired nitrogen loading reduction.
Measurement of planktonic chlorophyll‐a—a proxy for algal biomass—in rivers may represent local production or algae transported from upstream, confounding understanding of algal bloom development in ...flowing waters. We modeled 3 years of chlorophyll‐a transport through a 394‐km portion of the Illinois River and found that although algal biomass is longitudinally widespread, most net production occurs at river control points in the upper reaches (up to 3.7 Mg chlorophyll‐a y−1 km−1). Up to 69% of the algal biomass in the upper river was a result of within‐reach production, with the remainder recruited from headwaters and tributaries. High chlorophyll‐a measured farther downstream was largely because of transport from source‐area control points, with substantial net losses of algal biomass occurring in the lower river. Modeling the often‐overlooked river transport component is necessary to characterize where, when, and why planktonic algae grow and predict how far and fast they move downstream.
Plain Language Summary
Planktonic algae in rivers may accumulate during periods of high productivity stimulated by favorable light, temperature, nutrient, and flow conditions, which can disrupt ecological processes and affect human uses including recreation and drinking water supply. Planktonic algae observed in rivers may occur because of local growth or transport from upstream source areas. Therefore, considering both local and upstream conditions may improve early warnings of potentially harmful blooms. Along a 394‐km stretch of the Illinois River, we found that most of the algae grew in the upper reaches and was then transported to downstream reaches, contributing to potential downstream harms such as excessive turbidity, organic carbon, biological oxygen demand, and algal toxins. We demonstrate how the often‐overlooked river transport component can be quantified to better identify where, when, and why algae grow in river networks.
Key Points
Planktonic algal biomass is pervasive in the Illinois River, yet production is favored at certain locations and times
Most planktonic algal biomass was produced in upper‐reach control points that supplied downstream areas
Transport analysis using local and upstream data improves understanding of river algal blooms
Abstract
Nutrients that have gradually accumulated in soils, groundwaters, and river sediments in the United States over the past century can remobilize and increase current downstream loading, ...obscuring effects of conservation practices aimed at protecting water resources. Drivers of storage accumulation and release of nutrients are poorly understood at the spatial scale of basins to watersheds. Predicting water quality outcomes in large river basins demands modeling storage lags and time varying reactivity that models of mean conditions typically cannot elucidate. We developed a seasonally dynamic approach to large-scale nutrient modeling based on a multiscale framework and nutrient storage lags were quantified for the nearly 190 000 small catchments that feed the rivers across the northeastern United States where catchment mean transit times were found to be around 4.7 (2–10) years for nitrogen and 1.3 (0.7–2) years for phosphorus. Nutrient loads carried in river flow in the current season contained a significant—and sometimes dominant—portion of mass lagged in its release from catchment storage repositories. Our approach of integrating storage releases with seasonally dynamic hydroclimatic drivers sets the stage to assess the accumulated effects of nutrient storage and lagged releases to the river interacting with seasonally varying nutrient reactivity and societal management actions throughout large river basins.
The spatial and temporal scales of hyporheic exchange within the stream corridor are controlled by stream discharge and groundwater inflow interacting with streambed morphology. While decades of ...study have resulted in a clear understanding of how morphologic form controls hyporheic exchange at the feature scale, we lack comparable predictive power related to stream discharge and the spatial structure of groundwater inflows at the reach scale, where spatial heterogeneity in both geomorphic setting and hydrologic forcing are present. In this study, we simulated vertical hyporheic exchange along a 600 m mountain stream reach under high, medium, and low stream discharge while considering groundwater inflow as negligible, spatially uniform, or proportional to upslope accumulated area. Most changes to hyporheic flow path residence time or length in response to stream discharge were small (<5%), suggesting that discharge is a secondary control relative to morphologically driven hyporheic exchange. Groundwater inflow was a primary control and mostly caused decreases in hyporheic flow path residence time and length. This finding generally agrees with expectations from the literature; however, flow path response was not consistent across the study reach. Instead, we found that flow paths driven by large hydraulic gradients coinciding with large morphologic features were less sensitive to changes in groundwater inflow than those driven by hydraulic gradients similar to the valley gradient. Our results indicate that consideration of heterogeneous arrangement of morphologic features is necessary to differentiate between hyporheic flow paths that persist in time and those that are sensitive to changing hydrologic conditions.
Key Points
Hyporheic flow paths do not respond uniformly to different stream discharge and groundwater inflow conditions
Hyporheic flow paths with hydraulic gradients similar to the valley gradient are the most sensitive to changing hydrologic conditions
Hyporheic flow paths with hydraulic gradients set by large morphologic features are the least sensitive to changing hydrologic conditions
The relative roles of dynamic hydrologic forcing and geomorphology as controls on the timescales and magnitudes of stream‐aquifer exchange and hyporheic flow paths are unknown but required for ...management of stream corridors. We developed a comprehensive framework relating diel hydrologic fluctuations to hyporheic exchange in the absence of geomorphic complexity. We simulated groundwater flow through an aquifer bounded by a straight stream and hillslope and under time‐varying boundary conditions. We found that diel fluctuations can produce hyporheic flow path lengths and residence times that span orders of magnitude. With these results, hyporheic flow path residence times and lengths can be predicted from the timing and magnitude of diel fluctuations and valley slope. Finally, we demonstrated that dynamic hydrologic boundary conditions can produce spatial and temporal scales of hyporheic flow paths equivalent to those driven by many well‐studied geomorphic features, indicating that these controls must be considered together in future efforts of upscaling to stream networks.
Key Points
Diel hydrologic fluctuations create hyporheic flow path lengths and residence times that span orders of magnitude
Diel hydrologic fluctuations can produce hyporheic flow path residence times and lengths equivalent to those driven by geomorphic features
Hyporheic flow path residence time and geometry are not directly coupled
Exchange of water between streams and their hyporheic zones is known to be dynamic in response to hydrologic forcing, variable in space, and to exist in a framework with nested flow cells. The ...expected result of heterogeneous geomorphic setting, hydrologic forcing, and between‐feature interaction is hyporheic transit times that are highly variable in both space and time. Transit time distributions (TTDs) are important as they reflect the potential for hyporheic processes to impact biogeochemical transformations and ecosystems. In this study we simulate time‐variable transit time distributions based on dynamic vertical exchange in a headwater mountain stream with observed, heterogeneous step‐pool morphology. Our simulations include hyporheic exchange over a 600 m river corridor reach driven by continuously observed, time‐variable hydrologic conditions for more than 1 year. We found that spatial variability at an instance in time is typically larger than temporal variation for the reach. Furthermore, we found reach‐scale TTDs were marginally variable under all but the most extreme hydrologic conditions, indicating that TTDs are highly transferable in time. Finally, we found that aggregation of annual variation in space and time into a “master TTD” reasonably represents most of the hydrologic dynamics simulated, suggesting that this aggregation approach may provide a relevant basis for scaling from features or short reaches to entire networks.
Plain Language Summary
The exchange of water between streams and the shallow subsurface of their valleys is important to ecological functions, many of which depend upon the time it takes water to travel down the valley. For decades, researchers have focused on how river form affects this travel time, leaving a limited understanding of how travel times may change throughout the year due to different wetness conditions. Here we directly address this limitation by estimating travel times through a stream and its entire valley for a period of more than 1 year. We found that change in travel times caused by different wetness conditions throughout the year was smaller than variation caused by river form. Exceptions to this rule are during very dry times of year when streamflow is low or during very wet times of the year when streamflow is high such as those resulting from spring snowmelt flooding.
Key Points
Hyporheic transit time distributions are more variable in space than in time
Transit time distributions are not transferrable during periods of rapid changes in stream discharge
Relatively short study reaches may not capture representative spatial variation along a study reach
Small ponds—farm ponds, detention ponds, or impoundments below 0.01 km2—serve important human needs throughout most large river basins. Yet the role of small ponds in regional nutrient and sediment ...budgets is essentially unknown, currently making it impossible to evaluate their management potential to achieve water quality objectives. Here we used new hydrography data sets and found that small ponds, depending on their spatial position within both their local catchments and the larger river network, can dominate the retention of nitrogen, phosphorus, and sediment compared to rivers, lakes, and reservoirs. Over 300,000 small ponds are collectively responsible for 34%, 69%, and 12% of the mean annual retention of nitrogen, phosphorus, and sediment in the Northeastern United States, respectively, with a dominant influence in headwater catchments (54%, 85%, and 50%, respectively). Small ponds play a critical role among the many aquatic features in long‐term nutrient and sediment loading to downstream waters.
Plain Language Summary
Reservoirs created by river damming have extensive impacts on downstream water quality but are not necessarily the most important elements of a diverse aquatic landscape. Many more small ponds have been constructed to serve important human needs ranging from farm irrigation in agricultural areas to flood control and trapping of nutrients and fine sediment in urban areas. The number of human‐influenced small ponds is projected to rise worldwide, yet their role in the delivery of nutrients and sediment from headwaters to oceans is currently unresolved. Here we used new data sets and found that small ponds are collectively responsible for trapping a substantial amount of the nutrients and sediment that are exported annually from headwaters. These findings support the need to jointly consider features such as urban detention ponds, farm ponds, and beaver ponds in managing headwaters to decrease long‐term nutrient and sediment loading to downstream waters and sensitive coastal areas.
Key Points
Small ponds located in headwater catchments dominate nutrient and sediment retention compared to streams, rivers, lakes, and reservoirs
Small ponds located directly adjacent to streams or away in upland positions have distinct effects on nitrogen, phosphorus, and sediment
Finer‐scale small streams are minor net sources of phosphorus and major net sources of sediment where soil erodibility is high
The hydrologic connectivity between streams and their valley bottoms (stream corridor) is a critical determinant of their ecological function. Ecological functions are known to be spatially and ...temporally variable, but spatial dimensions of the problem are not easily quantified and thus they are usually overlooked. To estimate the spatial patterns of connectivity, and how connectivity varies with changes in discharge, we developed the hyporheic potential model. We used the model to interpret a series of solute tracer injections in two headwater mountain streams with contrasting valley bottom morphologies to estimate connectivity in the stream corridor. The distributions of flow path origination locations and the lengths of hyporheic flow paths appear to vary with base flow recession, even in cases where transport timescales are apparently unchanged. The modeled distribution of origination locations further allowed us to define a spatial analog to the temporal window of detection associated with solute tracer studies, and enables assessment of connectivity dynamics between streams and their corridors. Altogether, the reduced complexity hyporheic potential model provides an easy way to anticipate the spatial distribution and origination locations of hyporheic flow paths from a basic understanding of the valley bottom characteristics and solute transport timescales.
Plain Language Summary
The manuscript details a simple method to assess the spatial connectivity of streams and their riparian zones. While the timescales of exchange in the river corridor have been broadly studied, the complimentary spatial dimension (i.e., the geometry of exchange flowpaths) remains largely unknown. The major challenge in assessing the spatial dimensions of exchange is the limited information available in the subsurface. Here, we develop a reduced complexity model of valley bottom transport to overcome these information limitations. With this model, relatively simple field site characterization and solute tracer data are combined to assess the spatial distribution of downwelling along a headwater mountain stream. We validate the model with a numerical experiment, and demonstrate its application in two watersheds of contrasting geology, repeated through baseflow recession.
Key Points
Hyporheic flow path geometry varies with discharge, even in cases where transport times remain unchanged
In‐stream discharge and along‐stream morphology cannot be used to identify flow path origination locations
Observations of solute tracers have a spatial window of detection in addition to the more broadly recognized temporal window of detection
Solute transport along riparian and hyporheic flow paths is broadly expected to respond to dynamic hydrologic forcing by streams, aquifers, and hillslopes. However, direct observation of these ...dynamic responses is lacking, as is the relative control of geologic setting as a control on responses to dynamic hydrologic forcing. We conducted a series of four stream solute tracer injections through base flow recession in each of two watersheds with contrasting valley morphology in the H.J. Andrews Experimental Forest, monitoring tracer concentrations in the stream and in a network of shallow riparian wells in each watershed. We found hyporheic mean arrival time, temporal variance, and fraction of stream water in the bedrock‐constrained valley bottom and near large roughness elements in the wider valley bottom were not variable with discharge, suggesting minimal control by hydrologic forcing. Conversely, we observed increases in mean arrival time and temporal variance and decreasing fraction stream water with decreasing discharge near the hillslopes in the wider valley bottom. This may indicate changes in stream discharge and valley bottom hydrology control transport in less constrained locations. We detail five hydrogeomorphic responses to base flow recession to explain observed spatial and temporal patterns in the interactions between streams and their valley bottoms. Models able to account for the transition from geologically dominated processes in the near‐stream subsurface to hydrologically dominated processes near the hillslope will be required to predict solute transport and fate in valley bottoms of headwater mountain streams.
Key Points:
Hyporheic transport in constrained valleys is controlled by geologic setting
Hyporheic transport near large in‐stream features is not variable with discharge
Transport in broad riparian zones is controlled by dynamic hydrologic forcing