skip to main content

Title: Optimal model complexity for terrestrial carbon cycle prediction
Abstract. The terrestrial carbon cycle plays a critical role in modulating the interactions of climate with the Earth system, but different models often make vastly different predictions of its behavior. Efforts to reduce model uncertainty have commonly focused on model structure, namely by introducing additional processes and increasing structural complexity. However, the extent to which increased structural complexity can directly improve predictive skill is unclear. While adding processes may improve realism, the resulting models are often encumbered by a greater number of poorly determined or over-generalized parameters. To guide efficient model development, here we map the theoretical relationship between model complexity and predictive skill. To do so, we developed 16 structurally distinct carbon cycle models spanning an axis of complexity and incorporated them into a model–data fusion system. We calibrated each model at six globally distributed eddy covariance sites with long observation time series and under 42 data scenarios that resulted in different degrees of parameter uncertainty. For each combination of site, data scenario, and model, we then predicted net ecosystem exchange (NEE) and leaf area index (LAI) for validation against independent local site data. Though the maximum model complexity we evaluated is lower than most traditional terrestrial biosphere models, the more » complexity range we explored provides universal insight into the inter-relationship between structural uncertainty, parametric uncertainty, and model forecast skill. Specifically, increased complexity only improves forecast skill if parameters are adequately informed (e.g., when NEE observations are used for calibration). Otherwise, increased complexity can degrade skill and an intermediate-complexity model is optimal. This finding remains consistent regardless of whether NEE or LAI is predicted. Our COMPLexity EXperiment (COMPLEX) highlights the importance of robust observation-based parameterization for land surface modeling and suggests that data characterizing net carbon fluxes will be key to improving decadal predictions of high-dimensional terrestrial biosphere models. « less
; ; ; ; ; ; ; ; ; ; ; ;
Award ID(s):
Publication Date:
Journal Name:
Sponsoring Org:
National Science Foundation
More Like this
  1. As the Arctic region moves into uncharted territory under a warming climate, it is important to refine the terrestrial biosphere models (TBMs) that help us understand and predict change. One fundamental uncertainty in TBMs relates to model parameters, configuration variables internal to the model whose value can be estimated from data. We incorporate a version of the Terrestrial Ecosystem Model (TEM) developed for arctic ecosystems into the Predictive Ecosystem Analyzer (PEcAn) framework. PEcAn treats model parameters as probability distributions, estimates parameters based on a synthesis of available field data, and then quantifies both model sensitivity and uncertainty to a given parameter or suite of parameters. We examined how variation in 21 parameters in the equation for gross primary production influenced model sensitivity and uncertainty in terms of two carbon fluxes (net primary productivity and heterotrophic respiration) and two carbon (C) pools (vegetation C and soil C). We set up different parameterizations of TEM across a range of tundra types (tussock tundra, heath tundra, wet sedge tundra, and shrub tundra) in northern Alaska, along a latitudinal transect extending from the coastal plain near Utqiaġvik to the southern foothills of the Brooks Range, to the Seward Peninsula. TEM was most sensitive tomore »parameters related to the temperature regulation of photosynthesis. Model uncertainty was mostly due to parameters related to leaf area, temperature regulation of photosynthesis, and the stomatal responses to ambient light conditions. Our analysis also showed that sensitivity and uncertainty to a given parameter varied spatially. At some sites, model sensitivity and uncertainty tended to be connected to a wider range of parameters, underlining the importance of assessing tundra community processes across environmental gradients or geographic locations. Generally, across sites, the flux of net primary productivity (NPP) and pool of vegetation C had about equal uncertainty, while heterotrophic respiration had higher uncertainty than the pool of soil C. Our study illustrates the complexity inherent in evaluating parameter uncertainty across highly heterogeneous arctic tundra plant communities. It also provides a framework for iteratively testing how newly collected field data related to key parameters may result in more effective forecasting of Arctic change.« less
  2. Abstract

    Global estimates of the land carbon sink are often based on simulations by terrestrial biosphere models (TBMs). The use of a large number of models that differ in their underlying hypotheses, structure and parameters is one way to assess the uncertainty in the historical land carbon sink. Here we show that the atmospheric forcing datasets used to drive these TBMs represent a significant source of uncertainty that is currently not systematically accounted for in land carbon cycle evaluations. We present results from three TBMs each forced with three different historical atmospheric forcing reconstructions over the period 1850–2015. We perform an analysis of variance to quantify the relative uncertainty in carbon fluxes arising from the models themselves, atmospheric forcing, and model-forcing interactions. We find that atmospheric forcing in this set of simulations plays a dominant role on uncertainties in global gross primary productivity (GPP) (75% of variability) and autotrophic respiration (90%), and a significant but reduced role on net primary productivity and heterotrophic respiration (30%). Atmospheric forcing is the dominant driver (52%) of variability for the net ecosystem exchange flux, defined as the difference between GPP and respiration (both autotrophic and heterotrophic respiration). In contrast, for wildfire-driven carbon emissions modelmore »uncertainties dominate and, as a result, model uncertainties dominate for net ecosystem productivity. At regional scales, the contribution of atmospheric forcing to uncertainty shows a very heterogeneous pattern and is smaller on average than at the global scale. We find that this difference in the relative importance of forcing uncertainty between global and regional scales is related to large differences in regional model flux estimates, which partially offset each other when integrated globally, while the flux differences driven by forcing are mainly consistent across the world and therefore add up to a larger fractional contribution to global uncertainty.

    « less
  3. Globally, planted forests are rapidly replacing naturally regenerated stands but the implications for canopy structure, carbon (C) storage, and the linkages between the two are unclear. We investigated the successional dynamics, interlinkages and mechanistic relationships between wood net primary production (NPPw) and canopy structure in planted and naturally regenerated red pine (Pinus resinosa Sol. ex Aiton) stands spanning ≥ 45 years of development. We focused our canopy structural analysis on leaf area index (LAI) and a spatially integrative, terrestrial LiDAR-based complexity measure, canopy rugosity, which is positively correlated with NPPw in several naturally regenerated forests, but which has not been investigated in planted stands. We estimated stand NPPw using a dendrochronological approach and examined whether canopy rugosity relates to light absorption and light–use efficiency. We found that canopy rugosity increased similarly with age in planted and naturally regenerated stands, despite differences in other structural features including LAI and stem density. However, the relationship between canopy rugosity and NPPw was negative in planted and not significant in naturally regenerated stands, indicating structural complexity is not a globally positive driver of NPPw. Underlying the negative NPPw-canopy rugosity relationship in planted stands was a corresponding decline in light-use efficiency, which peaked in themore »youngest, densely stocked stand with high LAI and low structural complexity. Even with significant differences in the developmental trajectories of canopy structure, NPPw, and light use, planted and naturally regenerated stands stored similar amounts of C in wood over a 45-year period. We conclude that widespread increases in planted forests are likely to affect age-related patterns in canopy structure and NPPw, but planted and naturally regenerated forests may function as comparable long-term C sinks via different structural and mechanistic pathways.« less
  4. Abstract. Land models are essential tools for understanding and predicting terrestrial processes and climate–carbon feedbacks in the Earth system, but uncertainties in their future projections are poorly understood. Improvements in physical process realism and the representation of human influence arguably make models more comparable to reality but also increase the degrees of freedom in model configuration, leading to increased parametric uncertainty in projections. In this work we design and implement a machine learning approach to globally calibrate a subset of the parameters of the Community Land Model, version 5 (CLM5) to observations of carbon and water fluxes. We focus on parameters controlling biophysical features such as surface energy balance, hydrology, and carbon uptake. We first use parameter sensitivity simulations and a combination of objective metrics including ranked global mean sensitivity to multiple output variables and non-overlapping spatial pattern responses between parameters to narrow the parameter space and determine a subset of important CLM5 biophysical parameters for further analysis. Using a perturbed parameter ensemble, we then train a series of artificial feed-forward neural networks to emulate CLM5 output given parameter values as input. We use annual mean globally aggregated spatial variability in carbon and water fluxes as our emulation and calibrationmore »targets. Validation and out-of-sample tests are used to assess the predictive skill of the networks, and we utilize permutation feature importance and partial dependence methods to better interpret the results. The trained networks are then used to estimate global optimal parameter values with greater computational efficiency than achieved by hand tuning efforts and increased spatial scale relative to previous studies optimizing at a single site. By developing this methodology, our framework can help quantify the contribution of parameter uncertainty to overall uncertainty in land model projections.« less
  5. Abstract. Plant transpiration downregulation in the presence of soil water stress is a critical mechanism for predicting global water, carbon, and energy cycles. Currently, many terrestrial biosphere models (TBMs) represent this mechanism with an empirical correction function (β) of soil moisture – a convenient approach that can produce large prediction uncertainties. To reduce this uncertainty, TBMs have increasingly incorporated physically based plant hydraulic models (PHMs). However, PHMs introduce additional parameter uncertainty and computational demands. Therefore, understanding why and when PHM and β predictions diverge would usefully inform model selection within TBMs. Here, we use a minimalist PHM to demonstrate that coupling the effects of soil water stress and atmospheric moisture demand leads to a spectrum of transpiration responses controlled by soil–plant hydraulic transport (conductance). Within this transport-limitation spectrum, β emerges as an end-member scenario of PHMs with infinite conductance, completely decoupling the effects of soil water stress and atmospheric moisture demand on transpiration. As a result, PHM and β transpiration predictions diverge most for soil–plant systems with low hydraulic conductance (transport-limited) that experience high variation in atmospheric moisture demand and have moderate soil moisture supply for plants. We test these minimalist model results by using a land surface model atmore »an AmeriFlux site. At this transport-limited site, a PHM downregulation scheme outperforms the β scheme due to its sensitivity to variations in atmospheric moisture demand. Based on this observation, we develop a new “dynamic β” that varies with atmospheric moisture demand – an approach that overcomes existing biases within β schemes and has potential to simplify existing PHM parameterization and implementation.« less