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.

; ; ; ; ; ; ;
Award ID(s):
1926423 1637630 1926591
Publication Date:
Journal Name:
Nature Communications
Nature Publishing Group
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, suggestingmore »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.« less
  2. 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,more »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.« less
  3. Hypoxia and associated acidification are growing concerns for ecosystems and biogeochemical cycles in the coastal zone. The northern Gulf of Mexico (nGoM) has experienced large seasonal hypoxia for decades linked to the eutrophication of the continental shelf fueled by the Mississippi River nutrient discharge. Sediments play a key role in maintaining hypoxic and acidified bottom waters, but this role is still not completely understood. In the summer 2017, when the surface area of the hypoxic zone in the nGoM was the largest ever recorded, we investigated four stations on the continental shelf differentially influenced by river inputs of the Mississippi-Atchafalaya River System and seasonal hypoxia. We investigated diagenetic processes under normoxic, hypoxic, and nearly anoxic bottom waters by coupling amperometric, potentiometric, and voltammetric microprofiling with high-resolution diffusive equilibrium in thin-films (DET) profiles and porewater analyses. In addition, we used a time-series of bottom-water dissolved oxygen from May to November 2017, which indicated intense O 2 consumption in bottom waters related to organic carbon recycling. At the sediment-water interface (SWI), we found that oxygen consumption linked to organic matter recycling was large with diffusive oxygen uptake (DOU) of 8 and 14 mmol m –2 d –1 , except when the oxygenmore »concentration was near anoxia (5 mmol m –2 d –1 ). Except at the station located near the Mississippi river outlet, the downcore pore water sulfate concentration decrease was limited, with little increase in alkalinity, dissolved inorganic carbon (DIC), ammonium, and phosphate suggesting that low oxygen conditions did not promote anoxic diagenesis as anticipated. We attributed the low anoxic diagenesis intensity to a limitation in organic substrate supply, possibly linked to the reduction of bioturbation during the hypoxic spring and summer.« less
  4. While hydroelectric dams play a significant role in meeting the increasing energy demand worldwide, they pose a significant risk to riverine biodiversity and food security for millions of people that mainly depend upon floodplain fisheries. Dam structures could affect fish populations both directly and indirectly through loss of accessible spawning and rearing habitat, degradation of habitat quality (e.g., changes in temperature and discharge), and/or turbine injuries. However, our understandings of the impacts of dam life span and the initial fishery conditions on restoration time and hence the dynamic hydropower (energy)-fish (food) nexus remain limited. In this study, we explored the temporal energy-food tradeoffs associated with a hydroelectric dam located in the Penobscot River basin of the United States. We investigated the influence of dam life span, upstream passage rate, and downstream habitat area on the energy-food tradeoffs using a system dynamics model. Our results show that around 90% of fish biomass loss happen within 5 years of dam construction. Thereafter, fish decline slowly stabilizes and approaches the lowest value at around the 20th year after dam construction. Fish restoration period is highly sensitive even to a short period of blockage. The biomass of alewife spawners need 18 years to recovermore »with only 1-year of blockage to the upstream critical habitats. Hydropower generation and loss of fish biomass present a two-segment linear relationship under changes in dam life span. When the dam life span is less than 5 years, generating 1 GWh energy cause around 0.04 million kg loss of fish biomass; otherwise, the loss of fish biomass is 0.02 million kg. The loss of fish biomass could be significantly decreased with minimal energy loss through increasing upstream passage rate and/or the size of downstream habitat area.« less
  5. Abstract. Watersheds are the fundamental Earth surface functioning units that connect the land to aquatic systems. Many watershed-scale models represent hydrological processes but not biogeochemical reactive transport processes. This has limited our capability to understand and predict solute export, water chemistry and quality, and Earth system response to changing climate and anthropogenic conditions. Here we present a recently developed BioRT-Flux-PIHM (BioRT hereafter) v1.0, a watershed-scale biogeochemical reactive transport model. The model augments the previously developed RT-Flux-PIHM that integrates land-surface interactions, surface hydrology, and abiotic geochemical reactions. It enables the simulation of (1) shallow and deep-water partitioning to represent surface runoff, shallow soil water, and deeper groundwater and of (2) biotic processes including plant uptake, soil respiration, and nutrient transformation. The reactive transport part of the code has been verified against the widely used reactive transport code CrunchTope. BioRT-Flux-PIHM v1.0 has recently been applied in multiple watersheds under diverse climate, vegetation, and geological conditions. This paper briefly introduces the governing equations and model structure with a focus on new aspects of the model. It also showcases one hydrology example that simulates shallow and deep-water interactions and two biogeochemical examples relevant to nitrate and dissolved organic carbon (DOC). These examples are illustrated in twomore »simulation modes of complexity. One is the spatially lumped mode (i.e., two land cells connected by one river segment) that focuses on processes and average behavior of a watershed. Another is the spatially distributed mode (i.e., hundreds of cells) that includes details of topography, land cover, and soil properties. Whereas the spatially lumped mode represents averaged properties and processes and temporal variations, the spatially distributed mode can be used to understand the impacts of spatial structure and identify hot spots of biogeochemical reactions. The model can be used to mechanistically understand coupled hydrological and biogeochemical processes under gradients of climate, vegetation, geology, and land use conditions.« less