Sandy riverbeds are very rarely flat. They are typically covered by ripples and dunes. Because of their topography, these ripples and dunes drive variations in water pressure across their surfaces ...due to deflection, acceleration, and deceleration of the river flow. These pressure variations drive river water to infiltrate into the porous and permeable sediment where pressure is high and exit from the sediment where it is low. This pressure-driven flow, called hyporheic exchange, is critical to the water quality of rivers since it allows river water to undergo biogeochemical reactions that take place within the sediment. Ripples are highly dynamic however and respond readily to changes in river flow. How the migration and variable shape of ripples affect hyporheic exchange and the biogeochemical reactions it dictates is poorly understood and seldom studied. Here we bring concepts from ripple dynamics, river and groundwater hydraulics, and biogeochemistry into a unified modeling framework. The modeling was used to assess the effects of ripple migration on hyporheic zone biogeochemistry. We found that migrating ripples generally process less nitrate, a widespread pollutant, compared to their stationary counterparts. Thus, investigations and applications of hyporheic zone biogeochemical processes should pay attention to the dynamics of ripples.
River corridor metabolomes reflect organic matter (OM) processing that drives aquatic biogeochemical cycles. Recent work highlights the power of ultrahigh-resolution mass spectrometry for ...understanding metabolome composition and river corridor metabolism. However, there have been no studies on the global chemogeography of surface water and sediment metabolomes using ultrahigh-resolution techniques. Here, we describe a community science effort from the Worldwide Hydrobiogeochemistry Observation Network for Dynamic River Systems (WHONDRS) consortium to characterize global metabolomes in surface water and sediment that span multiple stream orders and biomes. We describe the distribution of key aspects of metabolomes including elemental groups, chemical classes, indices, and inferred biochemical transformations. We show that metabolomes significantly differ across surface water and sediment and that surface water metabolomes are more rich and variable. We also use inferred biochemical transformations to identify core metabolic processes shared among surface water and sediment. Finally, we observe significant spatial variation in sediment metabolites between rivers in the eastern and western portions of the contiguous United States. Our work not only provides a basis for understanding global patterns in river corridor biogeochemical cycles but also demonstrates that community science endeavors can enable global research projects that are unfeasible with traditional research models.
Hyporheic exchange in riverbeds is driven by current‐bed topography interactions. Because riverbeds exhibit topographic roughness across scales, from individual grains to bedforms and bars, they can ...exhibit fractal patterns. This study analyzed the influence of fractal properties of riverbed topography on hyporheic exchange. A set of synthetic fractal riverbeds with different scaling statistics was used as inputs to sequentially coupled numerical simulations of turbulent channel flow and hyporheic flow. In the analysis, the maximum power spectrum (dune size) and the fractal dimension (topographic complexity) were considered as independent variables and we then investigated how interfacial fluxes and hyporheic travel times are functionally related to these variables. As the maximum power spectrum increases (i.e., dune height to flow depth ratio), the average interfacial flux increases logarithmically whereas it increases exponentially with an increase in fractal dimension. Hyporheic exchange is more sensitive to additional roughness (larger fractal dimensions) than to bedform size (larger maximum power). Our results imply that fractal properties of riverbeds are crucial to predicting hyporheic exchange. The predictive relationships we propose could be integrated with reduced complexity, large‐scale models. They can also be used to design artificial topographies that target hyporheic ecosystem services.
Key Points
A series of numerical simulations were conducted to explore how the fractal properties of bedforms are related to hyporheic exchange
The results show that hyporheic exchange flux increases exponentially with respect to the bed topography fractal dimension
Fractal properties of riverbeds are crucial to predicting hyporheic exchange
Hydrologic exchange flows (HEFs) across the river‐aquifer interface have important implications for biogeochemical processes and contaminant plume migration in the river corridor, yet little is known ...about the hydrogeomorphic factors that control HEFs dynamics under dynamic flow conditions. Here, we developed a 3‐D numerical model for a large regulated river corridor along the Columbia River to study how HEFs are controlled by the interplays between dam‐regulated flow conditions and hydrogeomorphic features of such river corridor system. Our results revealed highly variable intra‐annual spatiotemporal patterns in HEFs along the 75‐km river reach, as well as strong interannual variability with larger exchange volumes in wet years than dry years. In general, the river was losing during late spring to early summer when the river stage was high, and river was gaining in fall and winter when river stage was low. The magnitude and timing of river stage fluctuations controlled the timing of high exchange rates. Both river channel geomorphology and the thickness of a highly permeable river bank geologic layer controlled the locations of exchange hot spots, while the latter played a dominant role. Dam‐induced, subdaily to daily river stage fluctuations drove high‐frequency variations in HEFs across the river‐aquifer interfaces, resulting in greater overall exchange volumes as compared to the case without high‐frequency flows. Our results demonstrated that upstream dam operations enhanced the exchange between river water and groundwater with strong potential influence on the associated biogeochemical processes and on the fate and transport of groundwater contaminant plumes in such river corridors.
Plain Language Summary
River stage fluctuations due to upstream dam operations enhance water exchange between rivers and aquifers. These exchange flows can stimulate biogeochemical processes in the subsurface. Here, we developed a 3‐D numerical model to simulate the flow and transport in a flood plain aquifer within the Hanford Reach, a free‐flowing section of Columbia River (WA). Our simulations revealed large spatial and temporal variability in exchange flow magnitude and direction in response to different river flow conditions. Furthermore, river channel morphology and subsurface geology together controlled the locations of high exchange flow rates. Our study provides important insights on the behavior of groundwater contaminant plumes in the river corridors of large regulated river systems.
Key Points
Hydrologic exchange flows across the riverbed show strong intra‐annual spatiotemporal patterns as well as interannual variabilities
Hot moments of HEFs are controlled by magnitude and timing of river stage fluctuations
Hot spots of HEFs are primarily controlled by the thickness of the underlying highly permeable sediment layer
Fine particles are critical to stream ecosystem functioning, influencing in‐stream processes from pathogen transmission to carbon cycling, all of which depend on particle immobilization. However, our ...ability to predict particle immobilization is limited by (1) availability of combined solute and particle tracer data and (2) identifying parameters that appropriately represent fine particle immobilization, due to the myriad of objective functions and model formulations. We found that improved predictions of the full distribution of possible fine particle residence times requires using an objective function that assesses both the peak and tailing of breakthrough curves together with solute tracers to constrain in‐stream transport processes. The representation of immobilization processes was significantly improved when solute tracer data were combined with a particle model, starkly contrasting the common assumption that fine particles transport as washload. We develop a clear strategy for improving fine particle transport predictions, reshaping the potential role of fine particles in water quality management.
Plain Language Summary
Fine particles, a general term that can be used to describe inorganic material like clays, particulate organic carbon, and harmful bacteria (i.e., pathogens), are important to stream functioning and water quality. The time it takes fine particles to move through a stream is difficult to predict primarily because there is no general guidance regarding the required data types and modeling approaches. Through comprehensive computational experiments and data analysis, we found that fine particles remain in streams much longer than commonly assumed, reshaping understanding of the role of fine particles in stream functioning and how waterborne pathogens are transmitted.
Key Points
Applying a balanced objective function and two tracers improves accuracy of particle retention up to 57% versus applying a washload assumption
Improved parameter identifiability found with an objective function that considers both peak and tailing of breakthrough curves
Particle immobilization can be predicted directly from solute data by adding particle transport terms to a model with hyporheic exchange
Braided reaches were common along near‐natural Alpine rivers, and the associated habitat dynamics supported plant and animal species specialized on early‐successional stages. The extensive riparian ...zones could mitigate climate change by absorbing floods and by retaining water during droughts. Human impacts largely reduced active river corridors through altered discharge and construction of dykes, while recent restoration projects aim at increasing river dynamics. The causes and consequences of Alpine river degradation are well understood, but there are only few quantitative studies on floodplain degradation and restoration. Thus, we have reconstructed historical changes of gravel bars along five Alpine rivers (Iller, Inn, Isar, Lech, and Wertach) in Southern Germany in the period 1808–2009, based on historical maps and aerial images. We found losses of >90% in gravel bar area along these rivers since the mid‐19th century. The decline was caused by a reduction of the active river corridor and by ongoing succession of the remaining open habitats. Within the past 30 years, at the Isar River, restoration measures were realized with the aim to widen the active river corridor and to recreate gravel bars. In four restored reaches, we found that 5% of the historical gravel bar area recovered, and that the proportion of restored gravel bar area was highest after intermediate flooding. We conclude that the active river corridors of German Alpine rivers are almost completely lost, and that more extensive restoration needs to be done to preserve the habitat dynamics and biodiversity of these systems, and to adapt Alpine rivers to climate change.
River stage fluctuations drive surface water‐groundwater exchanges within river corridors. This study evaluates how repeated daily stage fluctuations, representative of hydropeaking conditions, ...influence aerobic respiration of river‐sourced dissolved organic carbon (DOC) in the riparian exchange zone using reactive flow and transport simulations. Over 50 hypothetical scenarios were modelled to evaluate how the duration of the daily flood signal, river DOC concentration, aquifer hydraulic conductivity and ambient groundwater flow condition affect the fate and transport of DOC and DO in the riparian aquifer. Time series subsurface snapshots highlight how the various factors influence the subsurface distribution of DOC and DO. The total mass of DOC respired per meter of river had a wide range depending on the parameters, spanning from 1.4 to 71 g over 24‐h, with high hydraulic conductivity and losing ambient groundwater flow conditions favouring the largest amount of DOC respired. The ratio of DOC mass entering the riparian zone with the mass returning to the river showed that as little as 5% to as much as 76% of the DOC that enters the bank during stage fluctuations returns to the river. This return ratio is dependent on river DOC concentration, hydraulic conductivity and ambient groundwater flow. The results illustrate that stage variations due to river regulation can be a significant control on aerobic respiration in riparian exchange zones.
River stage variations due to river regulation is a significant control on aerobic respiration in riparian exchange zones.
The total mass of DOC respired per meter of river had a wide range with high hydraulic conductivity and losing ambient groundwater flow conditions favouring the largest amount of DOC respired.
5%–76% of the DOC mass entering the riparian zone returns to the river depending on river DOC concentration, hydraulic conductivity and ambient groundwater flow.
The gravel-bed Biała River, Polish Carpathians, was heavily affected by channelization and channel incision in the twentieth century. Not only were these impacts detrimental to the ecological state ...of the river, but they also adversely modified the conditions of floodwater retention and flood wave passage. Therefore, a few years ago an erodible corridor was delimited in two sections of the Biała to enable restoration of the river. In these sections, short, channelized reaches located in the vicinity of bridges alternate with longer, unmanaged channel reaches, which either avoided channelization or in which the channel has widened after the channelization scheme ceased to be maintained. Effects of these alternating channel morphologies on the conditions for flood flows were investigated in a study of 10 pairs of neighbouring river cross sections with constrained and freely developed morphology. Discharges of particular recurrence intervals were determined for each cross section using an empirical formula. The morphology of the cross sections together with data about channel slope and roughness of particular parts of the cross sections were used as input data to the hydraulic modelling performed with the one-dimensional steady-flow HEC-RAS software. The results indicated that freely developed cross sections, usually with multithread morphology, are typified by significantly lower water depth but larger width and cross-sectional flow area at particular discharges than single-thread, channelized cross sections. They also exhibit significantly lower average flow velocity, unit stream power, and bed shear stress. The pattern of differences in the hydraulic parameters of flood flows apparent between the two types of river cross sections varies with the discharges of different frequency, and the contrasts in hydraulic parameters between unmanaged and channelized cross sections are most pronounced at low-frequency, high-magnitude floods. However, because of the deep incision of the river, both cross section types are typified by a similar, low potential for the retention of floodwater in floodplain areas. The study indicated that even though river restoration has only begun here, it already brings beneficial effects for flood risk management, reducing flow energy and shear forces exerted on the bed and banks of the channel in unmanaged river reaches. Only within wide, unmanaged channel reaches can the flows of low-frequency, high-magnitude floods be conveyed with relatively low shear forces exerted on the channel boundary. In contrast, in channelized reaches, flow velocity and shear forces are substantially higher, inevitably causing bank erosion and channel incision.
•1D steady-flow model was used to examine the effect of mountain river restoration on hydraulic conditions of flood flows.•Unmanaged and channelized cross sections differ in hydraulic parameters of flood flows.•Hydraulic contrasts between both cross section types increase with flood magnitude.•River widening in freely developed sections reduced shear forces in the channel.•Floodwater retention in floodplain areas is low in unmanaged and channelized cross sections.
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.
While the hyporheic zone (HZ) accounts for a significant portion of whole stream CO2 concentrations, HZ respiration modeling studies are lacking in quantifying their contributions to the total CO2 at ...large watershed/basin scales. Quantifying the contribution of anaerobic respiration is also underappreciated. This study used a carbon‐nitrogen‐coupled river corridor model to quantify HZ aerobic and anaerobic respiration and determined the key factors controlling their spatial variability within the Columbia River Basin (CRB). The modeled respiration patterns showed high spatial variability. Among the nine sub‐basins composing the CRB, the Lower Columbia and the Willamette, which receive higher precipitation, had higher respiration. Medium‐sized rivers (fourth to sixth orders) produced the highest aerobic and anaerobic respiration among reaches of different sizes. At the basin scale, aerobic respiration is dominant, representing approximately 98.7% of the total respiration across the CRB. While most of the reaches were dominant with aerobic respiration, reaches in agricultural land showed a relatively higher anaerobic respiration (18%) ratio. A variable importance analysis showed that hyporheic exchange flux controlled most of the spatial variability of HZ respiration, dominating over other physical variables such as residence time, stream dissolved organic carbon (DOC), nitrate, and dissolved oxygen (DO). The influence of substrate concentration (DOC and DO) is larger in modeling anaerobic respiration than aerobic respiration. Future efforts will focus on improving the estimation of the HZ exchange flux and the implementation of spatially explicit parameterizations for the reactions of interest to reduce model uncertainty.
Plain Language Summary
Riverbeds generate high amounts of CO2, but their contribution to the total CO2 budgets in rivers is not well quantified via numerical simulation models for large regions. This study used a numerical simulation model to estimate the CO2 emissions from riverbeds into water columns in the presence and absence of oxygen and identified important variables explaining the spatial variation of riverbed CO2 emissions within the Columbia River Basin (CRB). Our modeling study found that CO2 emissions from riverbeds showed high spatial variability. Within the CRB, wetter sub‐basins showed higher CO2 emissions than drier sub‐basins. Medium‐sized rivers generated the highest CO2 emissions. Most channel CO2 emissions occurred in the presence of oxygen. However, reaches in agricultural lands generated relatively high CO2 emissions in the absence of oxygen. The water exchange rate between channel and riverbed can explain the spatial variation of CO2 emissions over other physical variables.
Key Points
Hyporheic exchange flux controls the spatial variation of respiration across reaches with different sizes and land uses
At the basin scale, hyporheic zone aerobic respiration accounts for about 98.7% of total simulated respiration
18% of total respiration is from anaerobic respiration in agricultural streams