An official website of the United States government
Abstract Robust ecological forecasting of tree growth under future climate conditions is critical to anticipate future forest carbon storage and flux. Here, we apply three ingredients of ecological forecasting that are key to improving forecast skill: data fusion, confronting model predictions with new data, and partitioning forecast uncertainty. Specifically, we present the first fusion of tree‐ring and forest inventory data within a Bayesian state‐space model at a multi‐site, regional scale, focusing onPinus ponderosavar.brachypterain the southwestern US. Leveraging the complementarity of these two data sources, we parsed the ecological complexity of tree growth into the effects of climate, tree size, stand density, site quality, and their interactions, and quantified uncertainties associated with these effects. New measurements of trees, an ongoing process in forest inventories, were used to confront forecasts of tree diameter with observations, and evaluate alternative tree growth models. We forecasted tree diameter and increment in response to an ensemble of climate change projections, and separated forecast uncertainty into four different causes: initial conditions, parameters, climate drivers, and process error. We found a strong negative effect of fall–spring maximum temperature, and a positive effect of water‐year precipitation on tree growth. Furthermore, tree vulnerability to climate stress increases with greater competition, with tree size, and at poor sites. Under future climate scenarios, we forecast increment declines of 22%–117%, while the combined effect of climate and size‐related trends results in a 56%–91% decline. Partitioning of forecast uncertainty showed that diameter forecast uncertainty is primarily caused by parameter and initial conditions uncertainty, but increment forecast uncertainty is mostly caused by process error and climate driver uncertainty. This fusion of tree‐ring and forest inventory data lays the foundation for robust ecological forecasting of aboveground biomass and carbon accounting at tree, plot, and regional scales, including iterative improvement of model skill.
more »
« less
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 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.
more »
« less
- Award ID(s):
- 1942133
- PAR ID:
- 10314190
- Date Published:
- Journal Name:
- Biogeosciences
- Volume:
- 18
- Issue:
- 8
- ISSN:
- 1726-4189
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Abstract Terrestrial ecosystems contribute most of the interannual variability (IAV) in atmospheric carbon dioxide (CO2) concentrations, but processes driving the IAV of net ecosystem CO2exchange (NEE) remain elusive. For a predictive understanding of the global C cycle, it is imperative to identify indicators associated with ecological processes that determine the IAV of NEE. Here, we decompose the annual NEE of global terrestrial ecosystems into their phenological and physiological components, namely maximum carbon uptake (MCU) and release (MCR), the carbon uptake period (CUP), and two parameters, α and β, that describe the ratio between actual versus hypothetical maximum C sink and source, respectively. Using long‐term observed NEE from 66 eddy covariance sites and global products derived from FLUXNET observations, we found that the IAV of NEE is determined predominately by MCU at the global scale, which explains 48% of the IAV of NEE on average while α, CUP, β, and MCR explain 14%, 25%, 2%, and 8%, respectively. These patterns differ in water‐limited ecosystems versus temperature‐ and radiation‐limited ecosystems; 31% of the IAV of NEE is determined by the IAV of CUP in water‐limited ecosystems, and 60% of the IAV of NEE is determined by the IAV of MCU in temperature‐ and radiation‐limited ecosystems. The Lund‐Potsdam‐Jena (LPJ) model and the Multi‐scale Synthesis and Terrestrial Model Inter‐comparison Project (MsTMIP) models underestimate the contribution of MCU to the IAV of NEE by about 18% on average, and overestimate the contribution of CUP by about 25%. This study provides a new perspective on the proximate causes of the IAV of NEE, which suggest that capturing the variability of MCU is critical for modeling the IAV of NEE across most of the global land surface.more » « less
-
Abstract A significant warming effect on arctic tundra is greening. Although this increase in predominantly woody vegetation has been linked to increases in gross primary productivity, increasing temperatures also stimulate ecosystem respiration. We present a novel analysis from small-scale plot measurements showing that the shape of the temperature- and light-dependent sink-to-source threshold (where net ecosystem exchange (NEE) equals zero) differs between two tussock tundra ecosystems differing in leaf area index (LAI). At the higher LAI site, the threshold is exceeded (i.e the ecosystem becomes a source) at relatively higher temperatures under low light but at lower temperatures under high light. At the lower LAI site, the threshold is exceeded at relatively lower temperatures under low light but at higher temperatures under high light. We confirmed this response at a single site where LAI was experimentally increased. This suggests the carbon balance of the tundra may be sensitive to small increases in temperature under low light, but that this effect may be significantly offset by increases in LAI. Importantly, we found that this LAI effect is reversed under high light, and so in a warming tundra, greater vegetation cover could have a progressively negative effect on net carbon uptake.more » « less
-
Abstract Climate change is having significant impacts on Earth’s ecosystems and carbon budgets, and in the Arctic may drive a shift from an historic carbon sink to a source. Large uncertainties in terrestrial biosphere models (TBMs) used to forecast Arctic changes demonstrate the challenges of determining the timing and extent of this possible switch. This spread in model predictions can limit the ability of TBMs to guide management and policy decisions. One of the most influential sources of model uncertainty is model parameterization. Parameter uncertainty results in part from a mismatch between available data in databases and model needs. We identify that mismatch for three TBMs, DVM-DOS-TEM, SIPNET and ED2, and four databases with information on Arctic and boreal above- and belowground traits that may be applied to model parametrization. However, focusing solely on such data gaps can introduce biases towards simple models and ignores structural model uncertainty, another main source for model uncertainty. Therefore, we develop a causal loop diagram (CLD) of the Arctic and boreal ecosystem that includes unquantified, and thus unmodeled, processes. We map model parameters to processes in the CLD and assess parameter vulnerability via the internal network structure. One important substructure, feed forward loops (FFLs), describe processes that are linked both directly and indirectly. When the model parameters are data-informed, these indirect processes might be implicitly included in the model, but if not, they have the potential to introduce significant model uncertainty. We find that the parameters describing the impact of local temperature on microbial activity are associated with a particularly high number of FFLs but are not constrained well by existing data. By employing ecological models of varying complexity, databases, and network methods, we identify the key parameters responsible for limited model accuracy. They should be prioritized for future data sampling to reduce model uncertainty.more » « less
-
Abstract Secondary forest regrowth shapes community succession and biogeochemistry for decades, including in the Upper Great Lakes region. Vegetation models encapsulate our understanding of forest function, and whether models can reproduce multi‐decadal succession patterns is an indication of our ability to predict forest responses to future change. We test the ability of a vegetation model to simulate C cycling and community composition during 100 years of forest regrowth following stand‐replacing disturbance, asking (a) Which processes and parameters are most important to accurately model Upper Midwest forest succession? (b) What is the relative importance of model structure versus parameter values to these predictions? We ran ensembles of the Ecosystem Demography model v2.2 with different representations of processes important to competition for light. We compared the magnitude of structural and parameter uncertainty and assessed which sub‐model–parameter combinations best reproduced observed C fluxes and community composition. On average, our simulations underestimated observed net primary productivity (NPP) and leaf area index (LAI) after 100 years and predicted complete dominance by a single plant functional type (PFT). Out of 4,000 simulations, only nine fell within the observed range of both NPP and LAI, but these predicted unrealistically complete dominance by either early hardwood or pine PFTs. A different set of seven simulations were ecologically plausible but under‐predicted observed NPP and LAI. Parameter uncertainty was large; NPP and LAI ranged from ~0% to >200% of their mean value, and any PFT could become dominant. The two parameters that contributed most to uncertainty in predicted NPP were plant–soil water conductance and growth respiration, both unobservable empirical coefficients. We conclude that (a) parameter uncertainty is more important than structural uncertainty, at least for ED‐2.2 in Upper Midwest forests and (b) simulating both productivity and plant community composition accurately without physically unrealistic parameters remains challenging for demographic vegetation models.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.5194/bg-18-2727-2021
