skip to main content


Title: Improving estimates and forecasts of lake carbon dynamics using data assimilation
Abstract

Lakes are biogeochemical hotspots on the landscape, contributing significantly to the global carbon cycle despite their small areal coverage. Observations and models of lake carbon pools and fluxes are rarely explicitly combined through data assimilation despite successful use of this technique in other fields. Data assimilation adds value to both observations and models by constraining models with observations of the system and by leveraging knowledge of the system formalized by the model to objectively fill observation gaps. In this article, we highlight the utility of data assimilation in lake carbon cycling research by using the ensemble Kalman filter to combine simple lake carbon models with observations of lake carbon pools and fluxes. We demonstrate that data assimilation helps reduce uncertainty in estimates of lake carbon pools and fluxes and more accurately estimate the true carbon pool size compared to estimates derived from observations alone. Data assimilation techniques should be embraced as valuable tools for lake biogeochemists interested in learning about ecosystem dynamics and forecasting ecosystem states and processes.

 
more » « less
Award ID(s):
1725386
NSF-PAR ID:
10461875
Author(s) / Creator(s):
 ;  ;  ;  ;  
Publisher / Repository:
Wiley Blackwell (John Wiley & Sons)
Date Published:
Journal Name:
Limnology and Oceanography: Methods
Volume:
17
Issue:
2
ISSN:
1541-5856
Page Range / eLocation ID:
p. 97-111
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    The Western United States is dominated by natural lands that play a critical role for carbon balance, water quality, and timber reserves. This region is also particularly vulnerable to forest mortality from drought, insect attack, and wildfires, thus requiring constant monitoring to assess ecosystem health. Carbon monitoring techniques are challenged by the complex mountainous terrain, thus there is an opportunity for data assimilation systems that combine land surface models and satellite‐derived observations to provide improved carbon monitoring. Here, we use the Data Assimilation Research Testbed to adjust the Community Land Model (CLM5.0) with remotely sensed observations of leaf area and above‐ground biomass. The adjusted simulation significantly reduced the above‐ground biomass and leaf area, leading to a reduction in both photosynthesis and respiration fluxes. The reduction in the carbon fluxes mostly offset, thus both the adjusted and free simulation projected a weak carbon sink to the land. This result differed from a separate observation‐constrained model (FLUXCOM) that projected strong carbon uptake to the land. Simulation diagnostics suggested water limitation had an important influence upon the magnitude and spatial pattern of carbon uptake through photosynthesis. We recommend that additional observations important for water cycling (e.g., snow water equivalent, land surface temperature) be included to improve the veracity of the spatial pattern in carbon uptake. Furthermore, the assimilation system should be enhanced to maximize the number of the simulated state variables that are adjusted, especially those related to the recommended observed quantities including water cycling and soil carbon.

     
    more » « less
  2. Abstract

    More frequent and severe droughts are driving increased forest mortality around the globe. We urgently need to describe and predict how drought affects forest carbon cycling and identify thresholds of environmental stress that trigger ecosystem collapse. Quantifying the effects of drought at an ecosystem level is complex because dynamic climate–plant relationships can cause rapid and/or prolonged shifts in carbon balance. We employ the CARbon DAta MOdel fraMework (CARDAMOM) to investigate legacy effects of drought on forest carbon pools and fluxes. Our Bayesian model‐data fusion approach uses tower observed meteorological forcing and carbon fluxes to determine the response and sensitivity of aboveground and belowground ecological processes associated with the 2012–2015 California drought. Our study area is a mid‐montane mixed conifer forest in the Southern Sierras. CARDAMOM constrained with gross primary productivity (GPP) estimates covering 2011–2017 show a ~75% reduction in GPP, compared to negligible GPP change when constrained with 2011 only. Precipitation across 2012–2015 was 45% (474 mm) lower than the historical average and drove a cascading depletion in soil moisture and carbon pools (foliar, labile, roots, and litter). Adding 157 mm during an especially stressful year (2014, annual rainfall = 293 mm) led to a smaller depletion of water and carbon pools, steering the ecosystem away from a state of GPP tipping‐point collapse to recovery. We present novel process‐driven insights that demonstrate the sensitivity of GPP collapse to ecosystem foliar carbon and soil moisture states—showing that the full extent of GPP response takes several years to arise. Thus, long‐term changes in soil moisture and carbon pools can provide a mechanistic link between drought and forest mortality. Our study provides an example for how key precipitation threshold ranges can influence forest productivity, making them useful for monitoring and predicting forest mortality events.

     
    more » « less
  3. Abstract

    Soil organic matter decomposition and its interactions with climate depend on whether the organic matter is associated with soil minerals. However, data limitations have hindered global-scale analyses of mineral-associated and particulate soil organic carbon pools and their benchmarking in Earth system models used to estimate carbon cycle–climate feedbacks. Here we analyse observationally derived global estimates of soil carbon pools to quantify their relative proportions and compute their climatological temperature sensitivities as the decline in carbon with increasing temperature. We find that the climatological temperature sensitivity of particulate carbon is on average 28% higher than that of mineral-associated carbon, and up to 53% higher in cool climates. Moreover, the distribution of carbon between these underlying soil carbon pools drives the emergent climatological temperature sensitivity of bulk soil carbon stocks. However, global models vary widely in their predictions of soil carbon pool distributions. We show that the global proportion of model pools that are conceptually similar to mineral-protected carbon ranges from 16 to 85% across Earth system models from the Coupled Model Intercomparison Project Phase 6 and offline land models, with implications for bulk soil carbon ages and ecosystem responsiveness. To improve projections of carbon cycle–climate feedbacks, it is imperative to assess underlying soil carbon pools to accurately predict the distribution and vulnerability of soil carbon.

     
    more » « less
  4. Abstract

    Most earth system models fail to capture the seasonality of carbon fluxes in radiation‐limited tropical evergreen forests (TEF) in the Amazon. Kim et al. (2012,https://doi.org/10.1111/j.1365-2486.2011.02629.x) first statistically incorporated a light‐controlled phenology module into an ecosystem model to improve carbon flux simulations at one TEF site. However, it is not clear how their approach can be extended to other TEF sites with different climatic conditions. Here we evaluated temporal variability in plant functional traits at three different TEF sites using a data‐conditioned stochastic parameterization method. We showed that previously studied links—between seasonal photosynthetically active radiation (PAR) and the traitsVcmax25and leaf longevity—occur across sites. We further determined that seasonal PAR could similarly drive variations in the stomatal conductance slope parameter. Differences found in temporal trait estimates among sites indicate that dynamic trait parameters cannot be applied uniformly over space, but it may be possible to extrapolate them based on climatic factors. Motivated by recent observations that physiological capacity develops as leaves mature, we built new regression models for predicting traits that not only include PAR but also an autoregressive lag term to capture observed physiological delays behind PAR‐driven phenology shifts. With our stochastic parameterization, we predicted the three sites to be carbon neutral or carbon sinks under the RCP 8.5 future climate scenario. In contrast, projections using standard static trait parameters show most of the Amazonian TEF region becoming a carbon source. We further approximated that variable traits may allow at least a third of the radiation‐limited TEF region in the Amazon to serve as a future net carbon sink.

     
    more » « less
  5. Abstract

    Wild large herbivores are declining worldwide. Despite extensive use of exclosure experiments to investigate herbivore impacts, there is little consensus on the effects of wild large herbivores on ecosystem function.

    Of the ecosystem functions likely impacted, we reviewed the five most‐studied in exclosure experiments: ecosystem resilience/resistance to disturbance, nutrient cycling, carbon cycling, plant regeneration, and primary productivity.

    Experimental data on large wild herbivores' effects on ecosystem functions were predominately derived from temperate grasslands (50% grasslands, 75% temperate zones). Additionally, data were from experiments that may not be of adequate size (median size 400 m2despite excluding all experiments below 25 m2) or duration (median duration 6 years) to capture ecosystem‐scale responses to these low‐density and wide‐ranging taxa.

    Wild herbivore removal frequently impacted ecosystem functions; for example, net carbon uptake increased by three times in some instances. However, the magnitude and direction of effects, even within a single function, were highly variable.

    A focus on carbon cycling highlighted challenges in interpreting effects on a single function. While the effect of large herbivore exclusion on carbon cycling was slightly positive when its components (e.g. pools vs. fluxes of carbon) were aggregated, effects on individual components were variable and sometimes opposed.

    Given modern declines in large wild herbivores, it is critical to understand their effects on ecosystem function. However, this synthesis highlights strong variability in direction, magnitude, and modifiers of these effects. Some variation is likely due to disparity in what components are used to describe a given function. For example, for the carbon cycle we identified eight distinctly meaningful components, which are not easily combined yet are potentially misrepresentative of the larger cycle when considered alone. However, much of the observed difference in responses likely reflects real ecological variability across complex systems.

    To move towards a general predictive framework we must identify where variation in effect is due to methodological differences and where due to ecosystem context. Two critical steps forward are (a) additional quantitative synthetic analyses of large herbivores' effects on individual functions, and (b) improved, increased systematic exclosure research focusing on effects of large herbivores' exclusion on functions.

    A freePlain Language Summarycan be found within the Supporting Information of this article.

     
    more » « less