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 more » 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. « less
; ; ; ; ; ; ; ; ; ; ; ;
Award ID(s):
Publication Date:
Journal Name:
Proceedings of the National Academy of Sciences
Sponsoring Org:
National Science Foundation
More Like this
  1. 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 sparselycovermore »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.« less
  2. Peatlands store substantial amounts of carbon and are vulnerable to climate change. We present a modified version of the Organising Carbon and Hydrology In Dynamic Ecosystems (ORCHIDEE) land surface model for simulating the hydrology, surface energy, and CO2 fluxes of peatlands on daily to annual timescales. The model includes a separate soil tile in each 0.5° grid cell, defined from a global peatland map and identified with peat-specific soil hydraulic properties. Runoff from non-peat vegetation within a grid cell containing a fraction of peat is routed to this peat soil tile, which maintains shallow water tables. The water table position separates oxic from anoxic decomposition. The model was evaluated against eddy-covariance (EC) observations from 30 northern peatland sites, with the maximum rate of carboxylation (Vcmax) being optimized at each site. Regarding short-term day-to-day variations, the model performance was good for gross primary production (GPP) (r2 =  0.76; Nash–Sutcliffe modeling efficiency, MEF  =  0.76) and ecosystem respiration (ER, r2 =  0.78, MEF  =  0.75), with lesser accuracy for latent heat fluxes (LE, r2 =  0.42, MEF  =  0.14) and and net ecosystem CO2 exchange (NEE, r2 =  0.38, MEF  =  0.26). Seasonal variations in GPP, ER, NEE, and energy fluxes on monthly scales showed moderate to high r2 values (0.57–0.86). For spatial across-site gradients of annual meanmore »GPP, ER, NEE, and LE, r2 values of 0.93, 0.89, 0.27, and 0.71 were achieved, respectively. Water table (WT) variation was not well predicted (r2 &lt; 0.1), likely due to the uncertain water input to the peat from surrounding areas. However, the poor performance of WT simulation did not greatly affect predictions of ER and NEE. We found a significant relationship between optimized Vcmax and latitude (temperature), which better reflects the spatial gradients of annual NEE than using an average Vcmax value.« less
  3. Accelerated warming in the Arctic has led to concern regarding the amount of carbon emission potential from Arctic water bodies. Yet, aquatic carbon dioxide (CO 2 ) and methane (CH 4 ) flux measurements remain scarce, particularly at high resolution and over long periods of time. Effluxes of methane (CH 4 ) and carbon dioxide (CO 2 ) from Toolik Lake, a deep glacial lake in northern Alaska, were measured for the first time with the direct eddy covariance (EC) flux technique during six ice-free lake periods (2010–2015). CO 2 flux estimates from the lake (daily average efflux of 16.7 ± 5.3 mmol m −2 d −1 ) were in good agreement with earlier estimates from 1975–1989 using different methods. CH 4 effluxes in 2010–2015 (averaging 0.13 ± 0.06 mmol m −2 d −1 ) showed an interannual variation that was 4.1 times greater than median diel variations, but mean fluxes were almost one order of magnitude lower than earlier estimates obtained from single water samples in 1990 and 2011–2012. The overall global warming potential (GWP) of Toolik Lake is thus governed mostly by CO 2 effluxes, contributing 86–93% of the ice-free period GWP of 26–90 g CO 2,eq mmore »−2 . Diel variation in fluxes was also important, with up to a 2-fold (CH 4 ) to 4-fold (CO 2 ) difference between the highest nighttime and lowest daytime effluxes. Within the summer ice-free period, on average, CH 4 fluxes increased 2-fold during the first half of the summer, then remained almost constant, whereas CO 2 effluxes remained almost constant over the entire summer, ending with a linear increase during the last 1–2 weeks of measurements. Due to the cold bottom temperatures of this 26 m deep lake, and the absence of ebullition and episodic flux events, Toolik Lake and other deep glacial lakes are likely not hot spots for greenhouse gas emissions, but they still contribute to the overall GWP of the Arctic.« less
  4. Streams in the southeastern United States Coastal Plains serve as an essential source of energy and nutrients for important estuarine ecosystems, and dissolved organic matter (DOM) exported from these streams can have profound impacts on the biogeochemical and ecological functions of fluvial networks. Here, we examined hydrological and temperature controls of DOM during low-flow periods from a forested stream located within the Coastal Plain physiographic region of Alabama, USA. We analyzed DOM via combining dissolved organic carbon (DOC) analysis, fluorescence excitation–emission matrix combined with parallel factor analysis (EEM-PARAFAC), and microbial degradation experiments. Four fluorescence components were identified: terrestrial humic-like DOM, microbial humic-like DOM, tyrosine-like DOM, and tryptophan-like DOM. Humic-like DOM accounted for ~70% of total fluorescence, and biodegradation experiments showed that it was less bioreactive than protein-like DOM that accounted for ~30% of total fluorescence. This observation indicates fluorescent DOM (FDOM) was controlled primarily by soil inputs and not substantially influenced by instream production and processing, suggesting that the bulk of FDOM in these streams is transported to downstream environments with limited in situ modification. Linear regression and redundancy analysis models identified that the seasonal variations in DOM were dictated primarily by hydrology and temperature. Overall, high discharge and shallowmore »flow paths led to the enrichment of less-degraded DOM with higher percentages of microbial humic-like and tyrosine-like compounds, whereas high temperatures favored the accumulation of high-aromaticity, high-molecular-weight, terrestrial, humic-like compounds in stream water. The flux of DOC and four fluorescence components was driven primarily by water discharge. Thus, the instantaneous exports of both refractory humic-like DOM and reactive protein-like DOM were higher in wetter seasons (winter and spring). As high temperatures and severe precipitation are projected to become more prominent in the southeastern U.S. due to climate change, our findings have important implications for future changes in the amount, source, and composition of DOM in Coastal Plain streams and the associated impacts on downstream carbon and nutrient supplies and water quality.« less
  5. Robust estimates for the rates and trends in terrestrial gross primary production (GPP; plant CO 2 uptake) are needed. Carbonyl sulfide (COS) is the major long-lived sulfur-bearing gas in the atmosphere and a promising proxy for GPP. Large uncertainties in estimating the relative magnitude of the COS sources and sinks limit this approach. Sulfur isotope measurements ( 34 S/ 32 S; δ 34 S) have been suggested as a useful tool to constrain COS sources. Yet such measurements are currently scarce for the atmosphere and absent for the marine source and the plant sink, which are two main fluxes. Here we present sulfur isotopes measurements of marine and atmospheric COS, and of plant-uptake fractionation experiments. These measurements resulted in a complete data-based tropospheric COS isotopic mass balance, which allows improved partition of the sources. We found an isotopic (δ 34 S ± SE) value of 13.9 ± 0.1‰ for the troposphere, with an isotopic seasonal cycle driven by plant uptake. This seasonality agrees with a fractionation of −1.9 ± 0.3‰ which we measured in plant-chamber experiments. Air samples with strong anthropogenic influence indicated an anthropogenic COS isotopic value of 8 ± 1‰. Samples of seawater-equilibrated-air indicate that the marine COSmore »source has an isotopic value of 14.7 ± 1‰. Using our data-based mass balance, we constrained the relative contribution of the two main tropospheric COS sources resulting in 40 ± 17% for the anthropogenic source and 60 ± 20% for the oceanic source. This constraint is important for a better understanding of the global COS budget and its improved use for GPP determination.« less