skip to main content


Title: Improving Predictions of Stream CO 2 Concentrations and Fluxes Using a Stream Network Model: A Case Study in the East River Watershed, CO, USA
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
Award ID(s):
2103520
NSF-PAR ID:
10362044
Author(s) / Creator(s):
 ;  
Publisher / Repository:
DOI PREFIX: 10.1029
Date Published:
Journal Name:
Global Biogeochemical Cycles
Volume:
35
Issue:
12
ISSN:
0886-6236
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    Streams and rivers are major sources of greenhouse gases (GHGs) to the atmosphere, as carbon and nitrogen are converted and outgassed during transport. Although our understanding of drivers of individual GHG fluxes has improved with numerous site‐specific studies and global‐scale compilations, our ability to parse out interrelated physical and biogeochemical drivers of gas concentrations is limited by a lack of consistently collected, temporally continuous samples of GHGs and their associated drivers. We present a first analysis of such a dataset collected by the National Ecological Observatory Network across 27 streams and rivers across ecoclimatic domains of the United States. Average concentrations of CO2ranged from 36.9 ± 0.88 to 404 ± 33 μmol L−1, CH4from 0.003 ± 0.0003 to 4.99 ± 0.72 μmol L−1, and N2O from 0.015 to 0.04 μmol L−1and spanned ranges of previous global compilations. Both CO2and CH4were strongly affected by physical drivers including mean air temperature and stream slope, as well as by dissolved oxygen and total nitrogen concentrations. N2O was exclusively correlated with total nitrogen concentrations. Results suggested that potential for gas exchange dominated patterns in gas concentrations at the site level, but contributions of in‐stream aerobic and anaerobic metabolism, and groundwater also likely varied across sites. The highest gas concentrations as well as highest variability occurred in low‐gradient, warmer, and nonperennial systems. These results are a first step in providing unprecedented, continuous estimates of GHG flux constrained by temporally variable physical and biogeochemical drivers of GHG production.

     
    more » « less
  2. Abstract

    Rivers and streams play an important role within the global carbon cycle, in part through emissions of carbon dioxide (CO2) to the atmosphere. However, the sources of this CO2and their spatiotemporal variability are difficult to constrain. Recent work has highlighted the role of carbonate buffering reactions that may serve as a source of CO2in high alkalinity systems. In this study, we seek to develop a quantitative framework for the role of carbonate buffering in the fluxes and spatiotemporal patterns of CO2and the stable and radio‐ isotope composition of dissolved inorganic carbon (DIC). We incorporate DIC speciation calculations of carbon isotopologues into a stream network CO2model and perform a series of simulations, ranging from the degassing of a groundwater seep to a hydrologically‐coupled 5th‐order stream network. We find that carbonate buffering reactions contribute >60% of emissions in high‐alkalinity, moderate groundwater‐CO2environments. However, atmosphere equilibration timescales of CO2are minimally affected, which contradicts hypotheses that carbonate buffering maintains high CO2across Strahler orders in high alkalinity systems. In contrast, alkalinity dramatically increases isotope equilibration timescales, which acts to decouple CO2and DIC variations from the isotopic composition even under low alkalinity. This significantly complicates a common method for carbon source identification. Based on similar impacts on atmospheric equilibration for stable and radio‐ carbon isotopologues, we develop a quantitative method for partitioning groundwater and stream corridor carbon sources in carbonate‐dominated watersheds. Together, these results provide a framework to guide fieldwork and interpretations of stream network CO2patterns across variable alkalinities.

     
    more » « less
  3. Abstract

    Streams in high‐elevation tropical ecosystems known as páramos may be significant sources of carbon dioxide (CO2) to the atmosphere by transforming terrestrial carbon to gaseous CO2. Studies of these environments are scarce, and estimates of CO2fluxes are poorly constrained. In this study, we use two independent methods for measuring gas transfer velocity (k), a critical variable in the estimation of CO2evasion and other biogeochemical processes. The first method, kinematick600(k600‐K), is derived from an empirical relationship between temperature‐adjustedk(k600) and the physical characteristics of the stream. The second method, measuredk600(k600‐M), estimates gas transfer velocity in the stream by in situ measurements of dissolved CO2(pCO2) and CO2evasion to the atmosphere, adjusting for temperature. Measurements were collected throughout a 5‐week period during the wet season of a peatland‐stream transition within a páramo ecosystem located above 4000 m in elevation in northeastern Ecuador. We characterized the spatial heterogeneity of the 250‐m reach on five occasions, and both methods showed a wide range of variability ink600at small spatial scales. Values ofk600‐Kranged from 7.42 to 330 m/d (mean = 116 ± 95.1 m/d), whereas values ofk600‐Mranged from 23.5 to 444 m/d (mean = 121 ± 127 m/d). Temporal variability ink600was driven by increases in stream discharge caused by rain events, whereas spatial variability was driven by channel morphology, including stream width and slope. The two methods were in good agreement (less than 16% difference) at high and medium stream discharge (above 7.0 L/s). However, the two methods considerably differed from one another (up to 73% difference) at low stream discharge (below 7.0 L/s, which represents 60% of the observations collected). Our study provides the first estimates ofk600values in a high‐elevation tropical catchment across steep environmental gradients and highlights the combined effects of hydrology and stream morphology in co‐regulating gas transfer velocities in páramo streams.

     
    more » « less
  4. Abstract

    Rivers and streams are control points for CO2emission to the air (fCO2), with emission rates often exceeding internal metabolism (net ecosystem production, NEP). The difference is usually attributed to CO2‐supersaturated groundwater inputs from upland soil respiration and rock weathering, but this implies a terrestrial‐to‐aquatic C transfer greater than estimated by terrestrial mass balance. One explanation is that riparian zones—rich in organic and inorganic C but mostly neglected in terrestrial mass balances—contribute disproportionately tofCO2. To test this hypothesis, we measuredfCO2, NEP, and the lateral CO2contributions from both terrestrial uplands (TER) and riparian wetlands (RIP) for seven reaches in a lowland river network in Florida, USA. NEP contributed about half offCO2, but the remaining CO2emission was generally much larger than measured TER. The relative importance of RIP versus TER varied markedly between contrasting hydrogeologic settings: RIP contributed 49% offCO2where geologic confinement forced lateral drainage through riparian soils, but only 12% where unconfined karst allowed deeper groundwater flowpaths that bypassed riparian zones. On a land area basis, the narrow riparian corridor yielded far more CO2than the terrestrial uplands (33.1 vs. 1.4 g‐C m−2 yr−1), resulting in river corridors (i.e., stream channel plus adjacent wetlands, NEP + RIP) sourcing 87% offCO2to streams. Our findings imply that true terrestrial CO2subsidies to streams may be smaller than previously estimated by aquatic mass balance and highlight the importance of explicitly integrating riparian zones into the conceptual model for terrestrial‐to‐aquatic C transfer.

     
    more » « less
  5. Abstract

    Riparian zones and the streams they border provide vital habitat for organisms, water quality protection, and other important ecosystem services. These areas are under pressure from local (land use/land cover change) to global (climate change) processes. Woody vegetation is expanding in grassland riparian zones worldwide. Here we report on a decade‐long watershed‐scale mechanical removal of woody riparian vegetation along 4.5 km of stream channel in a before–after control impact experiment. Prior to this removal, woody plants had expanded into grassy riparian areas, associated with a decline in streamflow, loss of grassy plant species, and other ecosystem‐scale impacts. We confirmed some expected responses, including rapid increases in stream nutrients and sediments, disappearance of stream mosses, and decreased organic inputs to streams via riparian leaves. We were surprised that nutrient and sediment increases were transient for 3 years, that there was no recovery of stream discharge, and that areas with woody removal did not shift back to a grassland state, even when reseeded with grassland species. Rapid expansion of shrubs (Cornus drummondii,Prunus americana) in the areas where trees were removed allowed woody vegetation to remain dominant despite repeating the cutting every 2 years. Our results suggest woody expansion can fundamentally alter terrestrial and aquatic habitat connections in grasslands, resulting in inexorable movement toward a new ecosystem state. Human pressures, such as climate change, atmospheric CO2increases, and elevated atmospheric nitrogen deposition, could continue to push the ecosystem along a trajectory that is difficult to change. Our results suggest that predicting relationships between riparian zones and the streams they border could be difficult in the face of global change in all biomes, even in well‐studied sites.

     
    more » « less