skip to main content


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):
1840243
NSF-PAR ID:
10380059
Author(s) / Creator(s):
; ; ; ; ; ; ; ; ; ; ; ;
Date Published:
Journal Name:
Proceedings of the National Academy of Sciences
Volume:
119
Issue:
11
ISSN:
0027-8424
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. 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
  2. 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
  3. 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
  4. 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 athttps://fluxnet.org/data/fluxnet-ch4-community-product/, 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 athttps://doi.org/10.5281/zenodo.4672601 (Delwiche et al., 2021).Additionally, raw FLUXNET-CH4 data used to extract seasonality parameterscan be downloaded from https://fluxnet.org/data/fluxnet-ch4-community-product/ (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
  5. Abstract. Despite their small spatial extent, fluvial ecosystems play a significant role in processing and transporting carbon in aquatic networks, which results in substantial emission of methane (CH4) into the atmosphere. For this reason, considerable effort has been put into identifying patterns and drivers of CH4 concentrations in streams and rivers and estimating fluxes to the atmosphere across broad spatial scales. However, progress toward these ends has been slow because of pronounced spatial and temporal variability of lotic CH4 concentrations and fluxes and by limited data availability across diverse habitats and physicochemical conditions. To address these challenges, we present a comprehensive database of CH4 concentrations and fluxes for fluvial ecosystems along with broadly relevant and concurrent physical and chemical data. The Global River Methane Database (GriMeDB; https://doi.org/10.6073/pasta/f48cdb77282598052349e969920356ef, Stanley et al., 2023) includes 24 024 records of CH4 concentration and 8205 flux measurements from 5029 unique sites derived from publications, reports, data repositories, unpublished data sets, and other outlets that became available between 1973 and 2021. Flux observations are reported as diffusive, ebullitive, and total CH4 fluxes, and GriMeDB also includes 17 655 and 8409 concurrent measurements of concentrations and 4444 and 1521 fluxes for carbon dioxide (CO2) and nitrous oxide (N2O), respectively. Most observations are date-specific (i.e., not site averages), and many are supported by data for 1 or more of 12 physicochemical variables and 6 site variables. Site variables include codes to characterize marginal channel types (e.g., springs, ditches) and/or the presence of human disturbance (e.g., point source inputs, upstream dams). Overall, observations in GRiMeDB encompass the broad range of the climatic, biological, and physical conditions that occur among world river basins, although some geographic gaps remain (arid regions, tropical regions, high-latitude and high-altitude systems). The global median CH4 concentration (0.20 µmol L−1) and diffusive flux (0.44 mmolm-2d-1) in GRiMeDB are lower than estimates from prior site-averaged compilations, although ranges (0 to 456 µmol L−1 and −136 to 4057 mmolm-2d-1) and standard deviations (10.69 and 86.4) are greater for this larger and more temporally resolved database. Available flux data are dominated by diffusive measurements despite the recognized importance of ebullitive and plant-mediated CH4 fluxes. Nonetheless, GriMeDB provides a comprehensive and cohesive resource for examining relationships between CH4 and environmental drivers, estimating the contribution of fluvial ecosystems to CH4 emissions, and contextualizing site-based investigations.

     
    more » « less