Attention:The NSF Public Access Repository (PAR) system and access will be unavailable from 11:00 PM ET on Thursday, October 15th until 12:00 AM ET on Friday, October 16th due to maintenance. We apologize for the inconvenience.


Title: Superlinear scaling of riverine biogeochemical function with watershed size
Abstract 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
PAR ID:
10318690
Author(s) / Creator(s):
; ; ; ; ; ; ;
Date Published:
Journal Name:
Nature Communications
Volume:
13
Issue:
1
ISSN:
2041-1723
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract Hyporheic zones regulate biogeochemical processes in streams and rivers, but high spatiotemporal heterogeneity makes it difficult to predict how these processes scale from individual reaches to river basins. Recent work applying allometric scaling (i.e., power‐law relationships between size and function) to river networks provides a new paradigm for understanding cumulative hyporheic biogeochemical processes. We used previously published model predictions of reach‐scale hyporheic aerobic respiration to explore patterns in allometric scaling across two climatically divergent basins with differing characteristics in the Pacific Northwest, United States. In the model, hydrologic exchange fluxes (HEFs) regulate hyporheic respiration, so we examined how HEFs might influence allometric scaling of respiration. We found consistent scaling behaviors where HEFs were either very low or very high, but differences between basins when HEFs were moderate. Our findings provide initial model‐generated hypotheses for factors influencing allometric scaling of hyporheic respiration. These hypotheses can be used to optimize new data generation efforts aimed at developing predictive understanding of allometries that can, in turn, be used to scale biogeochemical dynamics across watersheds. 
    more » « less
  2. 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
  3. Abstract Nitrogen (N) wet deposition chemistry impacts watershed biogeochemical cycling. The timescale and magnitude of (a)synchrony between wet deposition N inputs and watershed N outputs remains unresolved. We quantify deposition‐river N (a)synchrony with transfer entropy (TE), an information theory metric enabling quantification of lag‐dependent feedbacks in a hydrologic system by calculating directional information flow between variables. Synchrony is defined as a significant amount of TE‐calculated reduction in uncertainty of river N from wet deposition N after conditioning for antecedent river N conditions. Using long‐term timeseries of wet deposition and river DON, NO3−, and NH4+concentrations from the Lamprey River watershed, New Hampshire (USA), we constrain the role of wet deposition N to watershed biogeochemistry. Wet deposition N contributed information to river N at timescales greater than quick‐flow runoff generation, indicating that river N losses are a lagged non‐linear function of hydro‐biogeochemical forcings. River DON received the most information from all three wet deposition N solutes while wet deposition DON and NH4+contributed the most information to all three river N solutes. Information theoretic algorithms facilitated data‐driven inferences on the hydro‐biogeochemical processes influencing the fate of N wet deposition. For example, signals of mineralization and assimilation at a timescale of 12 to 21‐weeks lag display greater synchrony than nitrification, and we find that N assimilation is a positive lagged function of increasing N wet deposition. Although wet deposition N is not the main driver of river N, it contributes a significant amount of information resolvable at time scales of transport and transformations. 
    more » « less
  4. null (Ed.)
    Compound flooding is a physical phenomenon that has become more destructive in recent years. Moreover, compound flooding is a broad term that envelops many different physical processes that can range from preconditioned, to multivariate, to temporally compounding, or spatially compounding. This research aims to analyze a specific case of compound flooding related to tropical cyclones where the compounding effect is on coastal flooding due to a combination of storm surge and river discharge. In recent years, such compound flood events have increased in frequency and magnitude, due to a number of factors such as sea-level rise from warming oceans. Therefore, the ability to model such events is of increasing urgency. At present, there is no holistic, integrated modeling system capable of simulating or forecasting compound flooding on a large regional or global scale, leading to the need to couple various existing models. More specifically, two more challenges in such a modeling effort are determining the primary model and accounting for the effect of adjacent watersheds that discharge to the same receiving water body in amplifying the impact of compound flooding from riverine discharge with storm surge when the scale of the model includes an entire coastal line. In this study, we investigated the possibility of using the Advanced Circulation (ADCIRC) model as the primary model to simulate the compounding effects of fluvial flooding and storm surge via loose one-way coupling with gage data through internal time-dependent flux boundary conditions. The performance of the ADCIRC model was compared with the Hydrologic Engineering Center- River Analysis System (HEC-RAS) model both at the watershed and global scales. Furthermore, the importance of including riverine discharges and the interactions among adjacent watersheds were quantified. Results showed that the ADCIRC model could reliably be used to model compound flooding on both a watershed scale and a regional scale. Moreover, accounting for the interaction of river discharge from multiple watersheds is critical in accurately predicting flood patterns when high amounts of riverine flow occur in conjunction with storm surge. Particularly, with storms such as Hurricane Harvey (2017), where river flows were near record levels, inundation patterns and water surface elevations were highly dependent on the incorporation of the discharge input from multiple watersheds. Such an effect caused extra and longer inundations in some areas during Hurricane Harvey. Comparisons with real gauge data show that adding internal flow boundary conditions into ADCIRC to account for river discharge from multiple watersheds significantly improves accuracy in predictions of water surface elevations during coastal flooding events. 
    more » « less
  5. River deltas are densely populated, ecologically vital landscapes threatened by rising sea levels. Distributary channel networks disperse sediment to build deltaic land, yet the relationship between the network organization and land building remains elusive. Inspired by Hack’s law, which shows that watershed drainage area scales with channel length in tributary networks, we analyzed a global dataset of distributary networks and found a nearly identical scaling relationship between distributary channel length and nourishment area, the land-building counterpart to drainage area. Despite this apparent global scaling, we further identified two distinct local land-building patterns: uniform delta networks consistently follow Hack’s law, whereas composite delta networks exhibit a scale break, transitioning from space-filling growth around the delta apex to quasi-linear growth near the coast. The unexpected growth patterns suggest that global simplicity and local variability coexist in how river deltas grow and organize. 
    more » « less