skip to main content


Title: A Terrestrial‐Aquatic Model Reveals Cross‐Scale Interactions Regulate Lateral Dissolved Organic Carbon Transport From Terrestrial Ecosystems
Abstract

Lateral carbon transport (LCT), the flux of terrestrial C transported to aquatic ecosystems, displaces carbon (C) across the terrestrial‐aquatic continuum and is on the same order of magnitude as terrestrial net ecosystem production. However, few continental scale C models include LCT or the C‐hydrology linkages necessary for modeling LCT. Those that do exist, borrow processes and conceptual understanding from watershed scale models, assuming that large‐scale and small‐scale drivers of LCT are the same. We develop a conceptual framework of LCT, which focuses on lateral dissolved organic carbon (DOC) transport (LCT‐DOC), and operationalize it with a coupled terrestrial‐aquatic C and hydrology model. After comparing our model LCT‐DOC to previous estimates derived from a summation of landscape scale fluxes for the Contiguous U.S., we use model experiments to partition the importance of LCT‐DOC drivers including total annual precipitation, air temperature, and plant traits, which interact across regional and local scales. We find that climate is the strongest driver of LCT‐DOC, where LCT‐DOC is positively related to precipitation but inversely related to temperature at continental scales. However, the net effect of climate on LCT‐DOC is the product of cross‐scale interactions between climate and vegetation. Plant traits also interact strongly with climate and have a measurable influence on LCT‐DOC, with water use efficiency as the most influential plant trait because it couples terrestrial water and C cycling. We demonstrate that our conceptual framework and relatively simple linked C‐hydrology process model of LCT‐DOC can inform hypotheses and predict LCT‐DOC.

 
more » « less
Award ID(s):
1754561 1754363
NSF-PAR ID:
10445477
Author(s) / Creator(s):
 ;  ;  ;  
Publisher / Repository:
DOI PREFIX: 10.1029
Date Published:
Journal Name:
Journal of Geophysical Research: Biogeosciences
Volume:
127
Issue:
5
ISSN:
2169-8953
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    Climate change is thawing and potentially mobilizing vast quantities of organic carbon (OC) previously stored for millennia in permafrost soils of northern circumpolar landscapes. Climate‐driven increases in fire and thermokarst may play a key role in OC mobilization by thawing permafrost and promoting transport of OC. Yet, the extent of OC mobilization and mechanisms controlling terrestrial‐aquatic transfer are unclear. We demonstrate that hydrologic transport of soil dissolved OC (DOC) from the active layer and thawing permafrost to headwater streams is extremely heterogeneous and regulated by the interactions of soils, seasonal thaw, fire, and thermokarst. Repeated sampling of streams in eight headwater catchments of interior Alaska showed that the mean age of DOC for each stream ranges widely from modern to ∼2,000 years B.P. Together, an endmember mixing model and nonlinear, generalized additive models demonstrated that Δ14C‐DOC signature (and mean age) increased from spring to fall, and was proportional to hydrologic contributions from a solute‐rich water source, related to presumed deeper flow paths found predominantly in silty catchments. This relationship was correlated with and mediated by catchment properties. Mean DOC ages were older in catchments with >50% burned area, indicating that fire is also an important explanatory variable. These observations underscore the high heterogeneity in aged C export and difficulty of extrapolating estimates of permafrost‐derived DOC export from watersheds to larger scales. Our results provide the foundation for developing a conceptual model of permafrost DOC export necessary for advancing understanding and prediction of land‐water C exchange in changing boreal landscapes.

     
    more » « less
  2. BACKGROUND The availability of nitrogen (N) to plants and microbes has a major influence on the structure and function of ecosystems. Because N is an essential component of plant proteins, low N availability constrains the growth of plants and herbivores. To increase N availability, humans apply large amounts of fertilizer to agricultural systems. Losses from these systems, combined with atmospheric deposition of fossil fuel combustion products, introduce copious quantities of reactive N into ecosystems. The negative consequences of these anthropogenic N inputs—such as ecosystem eutrophication and reductions in terrestrial and aquatic biodiversity—are well documented. Yet although N availability is increasing in many locations, reactive N inputs are not evenly distributed globally. Furthermore, experiments and theory also suggest that global change factors such as elevated atmospheric CO 2 , rising temperatures, and altered precipitation and disturbance regimes can reduce the availability of N to plants and microbes in many terrestrial ecosystems. This can occur through increases in biotic demand for N or reductions in its supply to organisms. Reductions in N availability can be observed via several metrics, including lowered nitrogen concentrations ([N]) and isotope ratios (δ 15 N) in plant tissue, reduced rates of N mineralization, and reduced terrestrial N export to aquatic systems. However, a comprehensive synthesis of N availability metrics, outside of experimental settings and capable of revealing large-scale trends, has not yet been carried out. ADVANCES A growing body of observations confirms that N availability is declining in many nonagricultural ecosystems worldwide. Studies have demonstrated declining wood δ 15 N in forests across the continental US, declining foliar [N] in European forests, declining foliar [N] and δ 15 N in North American grasslands, and declining [N] in pollen from the US and southern Canada. This evidence is consistent with observed global-scale declines in foliar δ 15 N and [N] since 1980. Long-term monitoring of soil-based N availability indicators in unmanipulated systems is rare. However, forest studies in the northeast US have demonstrated decades-long decreases in soil N cycling and N exports to air and water, even in the face of elevated atmospheric N deposition. Collectively, these studies suggest a sustained decline in N availability across a range of terrestrial ecosystems, dating at least as far back as the early 20th century. Elevated atmospheric CO 2 levels are likely a main driver of declines in N availability. Terrestrial plants are now uniformly exposed to ~50% more of this essential resource than they were just 150 years ago, and experimentally exposing plants to elevated CO 2 often reduces foliar [N] as well as plant-available soil N. In addition, globally-rising temperatures may raise soil N supply in some systems but may also increase N losses and lead to lower foliar [N]. Changes in other ecosystem drivers—such as local climate patterns, N deposition rates, and disturbance regimes—individually affect smaller areas but may have important cumulative effects on global N availability. OUTLOOK Given the importance of N to ecosystem functioning, a decline in available N is likely to have far-reaching consequences. Reduced N availability likely constrains the response of plants to elevated CO 2 and the ability of ecosystems to sequester carbon. Because herbivore growth and reproduction scale with protein intake, declining foliar [N] may be contributing to widely reported declines in insect populations and may be negatively affecting the growth of grazing livestock and herbivorous wild mammals. Spatial and temporal patterns in N availability are not yet fully understood, particularly outside of Europe and North America. Developments in remote sensing, accompanied by additional historical reconstructions of N availability from tree rings, herbarium specimens, and sediments, will show how N availability trajectories vary among ecosystems. Such assessment and monitoring efforts need to be complemented by further experimental and theoretical investigations into the causes of declining N availability, its implications for global carbon sequestration, and how its effects propagate through food webs. Responses will need to involve reducing N demand via lowering atmospheric CO 2 concentrations, and/or increasing N supply. Successfully mitigating and adapting to declining N availability will require a broader understanding that this phenomenon is occurring alongside the more widely recognized issue of anthropogenic eutrophication. Intercalibration of isotopic records from leaves, tree rings, and lake sediments suggests that N availability in many terrestrial ecosystems has steadily declined since the beginning of the industrial era. Reductions in N availability may affect many aspects of ecosystem functioning, including carbon sequestration and herbivore nutrition. Shaded areas indicate 80% prediction intervals; marker size is proportional to the number of measurements in each annual mean. Isotope data: (tree ring) K. K. McLauchlan et al. , Sci. Rep. 7 , 7856 (2017); (lake sediment) G. W. Holtgrieve et al. , Science 334 , 1545–1548 (2011); (foliar) J. M. Craine et al. , Nat. Ecol. Evol. 2 , 1735–1744 (2018) 
    more » « less
  3. Abstract

    In the Arctic waterbodies are abundant and rapid thaw of permafrost is destabilizing the carbon cycle and changing hydrology. It is particularly important to quantify and accurately scale aquatic carbon emissions in arctic ecosystems. Recently available high-resolution remote sensing datasets capture the physical characteristics of arctic landscapes at unprecedented spatial resolution. We demonstrate how machine learning models can capitalize on these spatial datasets to greatly improve accuracy when scaling waterbody CO2and CH4fluxes across the YK Delta of south-west AK. We found that waterbody size and contour were strong predictors for aquatic CO2emissions, attributing greater than two-thirds of the influence to the scaling model. Small ponds (<0.001 km2) were hotspots of emissions, contributing fluxes several times their relative area, but were less than 5% of the total carbon budget. Small to medium lakes (0.001–0.1 km2) contributed the majority of carbon emissions from waterbodies. Waterbody CH4emissions were predicted by a combination of wetland landcover and related drivers, as well as watershed hydrology, and waterbody surface reflectance related to chromophoric dissolved organic matter. When compared to our machine learning approach, traditional scaling methods that did not account for relevant landscape characteristics overestimated waterbody CO2and CH4emissions by 26%–79% and 8%–53% respectively. This study demonstrates the importance of an integrated terrestrial-aquatic approach to improving estimates and uncertainty when scaling C emissions in the arctic.

     
    more » « less
  4. Terrestrial vegetation and soils hold three times more carbon than the atmosphere. Much debate concerns how anthropogenic activity will perturb these surface reservoirs, potentially exacerbating ongoing changes to the climate system. Uncertainties specifically persist in extrapolating point-source observations to ecosystem-scale budgets and fluxes, which require consideration of vertical and lateral processes on multiple temporal and spatial scales. To explore controls on organic carbon (OC) turnover at the river basin scale, we present radiocarbon ( 14 C) ages on two groups of molecular tracers of plant-derived carbon—leaf-wax lipids and lignin phenols—from a globally distributed suite of rivers. We find significant negative relationships between the 14 C age of these biomarkers and mean annual temperature and precipitation. Moreover, riverine biospheric-carbon ages scale proportionally with basin-wide soil carbon turnover times and soil 14 C ages, implicating OC cycling within soils as a primary control on exported biomarker ages and revealing a broad distribution of soil OC reactivities. The ubiquitous occurrence of a long-lived soil OC pool suggests soil OC is globally vulnerable to perturbations by future temperature and precipitation increase. Scaling of riverine biospheric-carbon ages with soil OC turnover shows the former can constrain the sensitivity of carbon dynamics to environmental controls on broad spatial scales. Extracting this information from fluvially dominated sedimentary sequences may inform past variations in soil OC turnover in response to anthropogenic and/or climate perturbations. In turn, monitoring riverine OC composition may help detect future climate-change–induced perturbations of soil OC turnover and stocks. 
    more » « 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 two 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. 
    more » « less