skip to main content


The NSF Public Access Repository (NSF-PAR) system and access will be unavailable from 11:00 PM ET on Thursday, June 13 until 2:00 AM ET on Friday, June 14 due to maintenance. We apologize for the inconvenience.

Title: The importance of hydrology in routing terrestrial carbon to the atmosphere via global streams and rivers
The magnitude of stream and river carbon dioxide (CO 2 ) emission is affected by seasonal changes in watershed biogeochemistry and hydrology. Global estimates of this flux are, however, uncertain, relying on calculated values for CO 2 and lacking spatial accuracy or seasonal variations critical for understanding macroecosystem controls of the flux. Here, we compiled 5,910 direct measurements of fluvial CO 2 partial pressure and modeled them against watershed properties to resolve reach-scale monthly variations of the flux. The direct measurements were then combined with seasonally resolved gas transfer velocity and river surface area estimates from a recent global hydrography dataset to constrain the flux at the monthly scale. Globally, fluvial CO 2 emission varies between 112 and 209 Tg of carbon per month. The monthly flux varies much more in Arctic and northern temperate rivers than in tropical and southern temperate rivers (coefficient of variation: 46 to 95 vs. 6 to 12%). Annual fluvial CO 2 emission to terrestrial gross primary production (GPP) ratio is highly variable across regions, ranging from negligible (<0.2%) to 18%. Nonlinear regressions suggest a saturating increase in GPP and a nonsaturating, steeper increase in fluvial CO 2 emission with discharge across regions, which leads to higher percentages of GPP being shunted into rivers for evasion in wetter regions. This highlights the importance of hydrology, in particular water throughput, in routing terrestrial carbon to the atmosphere via the global drainage networks. Our results suggest the need to account for the differential hydrological responses of terrestrial–atmospheric vs. fluvial–atmospheric carbon exchanges in plumbing the terrestrial carbon budget.  more » « less
Award ID(s):
Author(s) / Creator(s):
; ; ; ; ; ; ; ; ; ; ; ;
Date Published:
Journal Name:
Proceedings of the National Academy of Sciences
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract. Mapping in-situ eddy covariance measurements of terrestrial land-atmosphere fluxes to the globe is a key method for diagnosing the Earth system from a data-driven perspective. We describe the first global products (called X-BASE) from a newly implemented up-scaling framework, FLUXCOM-X. The X-BASE products comprise of estimates of CO2 net ecosys- tem exchange (NEE), gross primary productivity (GPP) as well as evapotranspiration (ET) and, for the first time, a novel fully data-driven global transpiration product (ETT ), at high spatial (0.05°) and temporal (hourly) resolution. X-BASE es- timates the global NEE at -5.75 ± 0.33 PgC·yr−1 for the period 2001-2020, showing a much higher consistency with independent atmospheric carbon cycle constraints compared to the previous versions of FLUXCOM. The improvement of global NEE was likely only possible thanks to the international effort to increase the precision and consistency of eddy co- variance collection and processing pipelines, as well as to the extension of the measurements to more site-years resulting in a wider coverage of bio-climatic conditions. However, X-BASE global net ecosystem exchange shows a very low inter-annual variability, which is common to state-of-the-art data-driven flux products and remains a scientific challenge. With 125 ± 2.1 P gC · yr−1 for the same period, X-BASE GP P is slightly higher than previous FLUXCOM estimates, mostly in temperate and boreal areas. X-BASE evapotranspiration amounts to 74.7x103 ± 0.9x103 km3 globally for the years 2001-2020, but ex- ceeds precipitation in many dry areas likely indicating overestimation in these regions. On average 57% of evapotranspiration are estimated to be transpiration, in good agreement with isotope-based approaches, but higher than estimates from many land surface models. Despite considerable improvements to the previous up-scaling products, many further opportunities for devel- opment exist. Pathways of exploration include methodological choices in the selection and processing of eddy-covariance and satellite observations, their ingestion into the framework, and the configuration of machine learning methods. For this the new FLUXCOM-X framework was specifically designed to have the necessary flexibility to experiment, diagnose, and converge to more accurate global flux estimates. 
    more » « less
  2. Abstract

    Dissolved organic carbon (DOC) flux from rivers in the pan‐Arctic watershed represents an important connection between major terrestrial carbon stocks and the Arctic Ocean. Previous estimates of Arctic carbon flux and dissolved organic matter (DOM) seasonal dynamics have relied predominantly on measurements from the six major Arctic rivers, yet these may not be representative of northern high‐latitude constrained smaller watersheds. Here, we evaluate DOC concentration and DOM composition in the Onega River, a small Arctic watershed, using optical measurements and ultrahigh resolution mass spectrometry. Compared to the six largest Arctic rivers, DOC, absorbance ata350, and indicators of terrestrial DOM (e.g., specific UV absorbance at 254 nm, modified aromaticity index, relative abundance of condensed aromatics and polyphenolics) were elevated in the Onega throughout the year. Seasonality was also generally muted in comparison to the major Arctic rivers with relatively elevated DOC and terrestrial markers in both spring and fall seasons. The Onega exhibits a strong relationship betweena350and DOC, and its organic‐rich nature is apparent in its high DOC yield (4.85 g m2yr−1), and higher chromophoric DOM per unit DOC than the six largest Arctic rivers. As DOC yield from the Onega may be more representative of smaller northern high‐latitude rivers, we derived a new pan‐Arctic DOC flux scaling estimate which is over 50% higher than previous estimates scaled solely from the six major Arctic rivers. These observations suggest that smaller northern high‐latitude rivers may be underrepresented in Arctic carbon flux models and highlights uncertainty around constraining the export of DOC to the Arctic Ocean.

    more » « less
  3. Abstract

    Inland waters are an important component of the global carbon budget. However, our ability to predict carbon fluxes from stream systems remains uncertain, aspCO2varies within streams at scales of 1–100 m. This makes direct monitoring of large‐scale CO2fluxes impractical. We incorporate CO2input and output fluxes into a stream network advection‐reaction model, representing the first process‐based representation of stream CO2dynamics at watershed scales. This model includes groundwater (GW) CO2inputs, water column (WC), benthic hyporheic zone (BHZ) respiration, downstream advection, and atmospheric exchange. We evaluate this model against existing statistical methods including upscaling and multiple linear regressions through comparisons to high‐resolution streampCO2data collected across the East River Watershed in the Colorado Rocky Mountains (USA). The stream network model accurately captures GW, evasion, and respiration‐drivenpCO2variability and significantly outperforms multiple linear regressions for predictingpCO2. Further, the model provides estimates of CO2contributions from internal versus external sources suggesting that streams transition from GW‐ to BHZ‐dominated sources between 3rd and 4th Strahler orders, with GW, BHZ, and WC accounting for 49.3%, 50.6%, and 0.1% of CO2fluxes from the watershed, respectively. Lastly, stream network model atmospheric CO2fluxes are 4‐12x times smaller than upscaling technique predictions, largely due to relationships between streampCO2and gas exchange velocities. Taken together, this stream network model improves our ability to predict stream CO2dynamics and efflux. Furthermore, future applications to regional and global scales may result in a significant downward revision of global flux estimates.

    more » « less
  4. Arctic rivers provide an integrated signature of the changing landscape and transmit signals of change to the ocean. Here, we use a decade of particulate organic matter (POM) compositional data to deconvolute multiple allochthonous and autochthonous pan-Arctic and watershed-specific sources. Constraints from carbon-to-nitrogen ratios (C:N), δ 13 C, and Δ 14 C signatures reveal a large, hitherto overlooked contribution from aquatic biomass. Separation in Δ 14 C age is enhanced by splitting soil sources into shallow and deep pools (mean ± SD: −228 ± 211 vs. −492 ± 173‰) rather than traditional active layer and permafrost pools (−300 ± 236 vs. −441 ± 215‰) that do not represent permafrost-free Arctic regions. We estimate that 39 to 60% (5 to 95% credible interval) of the annual pan-Arctic POM flux (averaging 4,391 Gg/y particulate organic carbon from 2012 to 2019) comes from aquatic biomass. The remainder is sourced from yedoma, deep soils, shallow soils, petrogenic inputs, and fresh terrestrial production. Climate change-induced warming and increasing CO 2 concentrations may enhance both soil destabilization and Arctic river aquatic biomass production, increasing fluxes of POM to the ocean. Younger, autochthonous, and older soil-derived POM likely have different destinies (preferential microbial uptake and processing vs. significant sediment burial, respectively). A small (~7%) increase in aquatic biomass POM flux with warming would be equivalent to a ~30% increase in deep soil POM flux. There is a clear need to better quantify how the balance of endmember fluxes may shift with different ramifications for different endmembers and how this will impact the Arctic system. 
    more » « less
  5. null (Ed.)
    Abstract. Methane (CH4) emissions from natural landscapes constituteroughly half of global CH4 contributions to the atmosphere, yet largeuncertainties remain in the absolute magnitude and the seasonality ofemission quantities and drivers. Eddy covariance (EC) measurements ofCH4 flux are ideal for constraining ecosystem-scale CH4emissions due to quasi-continuous and high-temporal-resolution CH4flux measurements, coincident carbon dioxide, water, and energy fluxmeasurements, lack of ecosystem disturbance, and increased availability ofdatasets over the last decade. Here, we (1) describe the newly publisheddataset, FLUXNET-CH4 Version 1.0, the first open-source global dataset ofCH4 EC measurements (available at, last access: 7 April 2021). FLUXNET-CH4includes half-hourly and daily gap-filled and non-gap-filled aggregatedCH4 fluxes and meteorological data from 79 sites globally: 42freshwater wetlands, 6 brackish and saline wetlands, 7 formerly drainedecosystems, 7 rice paddy sites, 2 lakes, and 15 uplands. Then, we (2) evaluate FLUXNET-CH4 representativeness for freshwater wetland coverageglobally because the majority of sites in FLUXNET-CH4 Version 1.0 arefreshwater wetlands which are a substantial source of total atmosphericCH4 emissions; and (3) we provide the first global estimates of theseasonal variability and seasonality predictors of freshwater wetlandCH4 fluxes. Our representativeness analysis suggests that thefreshwater wetland sites in the dataset cover global wetland bioclimaticattributes (encompassing energy, moisture, and vegetation-relatedparameters) in arctic, boreal, and temperate regions but only sparselycover humid tropical regions. Seasonality metrics of wetland CH4emissions vary considerably across latitudinal bands. In freshwater wetlands(except those between 20∘ S to 20∘ N) the spring onsetof elevated CH4 emissions starts 3 d earlier, and the CH4emission season lasts 4 d longer, for each degree Celsius increase in meanannual air temperature. On average, the spring onset of increasing CH4emissions lags behind soil warming by 1 month, with very few sites experiencingincreased CH4 emissions prior to the onset of soil warming. Incontrast, roughly half of these sites experience the spring onset of risingCH4 emissions prior to the spring increase in gross primaryproductivity (GPP). The timing of peak summer CH4 emissions does notcorrelate with the timing for either peak summer temperature or peak GPP.Our results provide seasonality parameters for CH4 modeling andhighlight seasonality metrics that cannot be predicted by temperature or GPP(i.e., seasonality of CH4 peak). FLUXNET-CH4 is a powerful new resourcefor diagnosing and understanding the role of terrestrial ecosystems andclimate drivers in the global CH4 cycle, and future additions of sitesin tropical ecosystems and site years of data collection will provide addedvalue to this database. All seasonality parameters are available at (Delwiche et al., 2021).Additionally, raw FLUXNET-CH4 data used to extract seasonality parameterscan be downloaded from (last access: 7 April 2021), and a completelist of the 79 individual site data DOIs is provided in Table 2 of this paper. 
    more » « less