skip to main content

Title: Superlinear scaling of riverine biogeochemical function with watershed size

River networks regulate carbon and nutrient exchange between continents, atmosphere, and oceans. However, contributions of riverine processing are poorly constrained at continental scales. Scaling relationships of cumulative biogeochemical function with watershed size (allometric scaling) provide an approach for quantifying the contributions of fluvial networks in the Earth system. Here we show that allometric scaling of cumulative riverine function with watershed area ranges from linear to superlinear, with scaling exponents constrained by network shape, hydrological conditions, and biogeochemical process rates. Allometric scaling is superlinear for processes that are largely independent of substrate concentration (e.g., gross primary production) due to superlinear scaling of river network surface area with watershed area. Allometric scaling for typically substrate-limited processes (e.g., denitrification) is linear in river networks with high biogeochemical activity or low river discharge but becomes increasingly superlinear under lower biogeochemical activity or high discharge, conditions that are widely prevalent in river networks. The frequent occurrence of superlinear scaling indicates that biogeochemical activity in large rivers contributes disproportionately to the function of river networks in the Earth system.

more » « less
Award ID(s):
1926423 1637630 1926591
Author(s) / Creator(s):
; ; ; ; ; ; ;
Publisher / Repository:
Nature Publishing Group
Date Published:
Journal Name:
Nature Communications
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Fecal contamination is a significant source of water quality impairment globally. Aquatic ecosystems can provide an important ecosystem service of fecal contamination removal. Understanding the processes that regulate the removal of fecal contamination among river networks across flow conditions is critical. We applied a river network model, the Framework for Aquatic Modeling in the Earth System (FrAMES-Ecoli), to quantify removal of fecal indicator bacteria by river networks across flow conditions during summers in a series of New England watersheds of different characteristics. FrAMES-Ecoli simulates sources, transport, and riverine removal of Escherichia coli (E. coli). Aquatic E. coli removal was simulated in both the water column and the hyporheic zone, and is a function of hydraulic conditions, flow exchange rates with the hyporheic zone, and die-off in each compartment. We found that, at the river network scale during summers, removal by river networks can be high (19–99%) with variability controlled by hydrologic conditions, watershed size, and distribution of sources in the watershed. Hydrology controls much of the variability, with 68–99% of network scale inputs removed under base flow conditions and 19–85% removed during storm events. Removal by the water column alone could not explain the observed pattern in E. coli, suggesting that processes such as hyporheic removal must be considered. These results suggest that river network removal of fecal indicator bacteria should be taken into consideration in managing fecal contamination at critical downstream receiving waters. 
    more » « less
  2. We present a multimodel analysis for mechanistic hypothesis testing in landscape evolution theory. The study site is a watershed with well‐constrained initial and boundary conditions in which a river network locally incised 50 m over the last 13 ka. We calibrate and validate a set of 37 landscape evolution models designed to hierarchically test elements of complexity from four categories: hillslope processes, channel processes, surface hydrology, and representation of geologic materials. Comparison of each model to a base model, which uses stream power channel incision, uniform lithology, hillslope transport by linear diffusion, and surface water discharge proportional to drainage area, serves as a formal test of which elements of complexity improve model performance. Model fit is assessed using an objective function based on a direct difference between observed and simulated modern topography. A hybrid optimization scheme identifies optimal parameters and uncertainty. Multimodel analysis determines which elements of complexity improve simulation performance. Validation tests which model improvements persist when models are applied to an independent watershed. The three most important model elements are (1) spatial variation in lithology (differentiation between shale and glacial till), (2) a fluvial erosion threshold, and (3) a nonlinear relationship between slope and hillslope sediment flux. Due to nonlinear interactions between model elements, some process representations (e.g., nonlinear hillslopes) only become important when paired with the inclusion of other processes (e.g., erosion thresholds). This emphasizes the need for caution in identifying the minimally sufficient process set. Our approach provides a general framework for hypothesis testing in landscape evolution.

    more » « less
  3. Abstract

    Most terrestrial allochthonous organic matter enters river networks through headwater streams during high flow events. In headwaters, allochthonous inputs are substantial and variable, but become less important in streams and rivers with larger watersheds. As allochthonous dissolved organic matter (DOM) moves downstream, the proportion of less aromatic organic matter with autochthonous characteristics increases. How environmental factors converge to control this transformation of DOM at a continental scale is less certain. We hypothesized that the amount of time water has spent travelling through surface waters of inland systems (streams, rivers, lakes, and reservoirs) is correlated to DOM composition. To test this hypothesis, we used established river network scaling relationships to predict relative river network flow‐weighted travel time (FWTT) of water for 60 stream and river sites across the contiguous United States (3090 discrete samples over 10 water years). We estimated lentic contribution to travel times with upstream in‐network lake and reservoir volume. DOM composition was quantified using ultraviolet and visible absorption and fluorescence spectroscopy. A combination of FWTT and lake and reservoir volume was the best overall predictor of DOM composition among models that also incorporated discharge, specific discharge, watershed area, and upstream channel length. DOM spectral slope ratio (R2 = 0.77) and Freshness Index (R2 = 0.78) increased and specific ultraviolet absorbance at 254 nm (R2 = 0.68) and Humification Index (R2 = 0.44) decreased across sites as a function of FWTT and upstream lake volume. This indicates autochthonous‐like DOM becomes continually more dominant in waters with greater FWTT. We assert that river FWTT can be used as a metric of the continuum of DOM composition from headwaters to rivers. The nature of the changes to DOM composition detected suggest this continuum is driven by a combination of photo‐oxidation, biological processes, hydrologically varying terrestrial subsidies, and aged groundwater inputs.

    more » « less
  4. Abstract

    Fluvial networks integrate, transform, and transport constituents from terrestrial and aquatic ecosystems. To date, most research on water quality dynamics has focused on process understanding at individual streams, and, as a result, there is a lack of studies analyzing how physical and biogeochemical drivers scale across fluvial networks. We performed tracer tests in five stream orders of the Jemez River continuum in New Mexico, USA, to quantify reach‐scale hyporheic exchange during two different seasonal periods to address the following: How do hyporheic zone contributions to overall riverine processing change with space and time? And does the spatiotemporal variability of hyporheic exchange scale across fluvial networks? Combining conservative (i.e., bromide) and reactive (i.e., resazurin) tracer analyses with solute transport modeling, we found a dominance of reaction‐limited transport conditions and a decrease of the contributions of hyporheic processing across stream orders and flow regimes. Our field‐based findings suggest that achieving knowledge transferability of hyporheic processing within fluvial networks may be possible, especially when process variability is sampled across multiple stream orders and flow regimes. Therefore, we propose a shift in our traditional approach to investigating scaling patterns in transport processes, which currently relies on the interpretation of studies conducted in multiple sites (mainly in headwater streams) that are located in different fluvial networks, to a more cohesive, network‐centered investigation of processes using the same or readily comparable methods.

    more » « less
  5. Abstract

    River networks transport dissolved organic carbon (DOC) from terrestrial uplands to the coastal ocean. The extent to which a reach or lake within a river network uptakes DOC depends on the stream order, the seasonal conditions, and the flow. At the watershed scale, it remains unclear whether DOC uptake is dominated by biological processes such as respiration, or abiotic processes like photomineralization. The partitioning of DOC uptake in lakes vs. rivers is also unclear. In this study, we present a new model that unifies year‐round controls on DOC cycling for an entire river network, including river–lake connectivity, to elucidate the importance of biotic vs. abiotic controls on DOC uptake. We present the Catchment UPtake and Sinks by Season, Order, and Flow for DOC (CUPS‐OF‐DOC) model, which quantifies terrestrial DOC loading, gross primary productivity, and uptake via microbes and photomineralization. The model is applied to the Connecticut River Watershed, and accounts for cascading reach‐ and lake‐scale DOC cycling across 98 scenarios spanning combinations of flows, seasons, and stream orders. We show that riverine DOC uptake is nearly constant with stream order, but the proportion of DOC uptake from photomineralization varies. Photomineralization dominates in rivers in most flow conditions and stream orders, especially in winter, accounting for at least half of whole‐watershed DOC uptake in February across all flows. Whole‐watershed summer DOC uptake occurs mostly via biomineralization in lakes, accounting for 80% of DOC uptake during the growing season, despite accounting for less than 6% of watershed open water surface area.

    more » « less