skip to main content


Title: Structure and parameter uncertainty in centennial projections of forest community structure and carbon cycling
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
Award ID(s):
1655095
NSF-PAR ID:
10374517
Author(s) / Creator(s):
 ;  ;  ;  
Publisher / Repository:
Wiley-Blackwell
Date Published:
Journal Name:
Global Change Biology
Volume:
26
Issue:
11
ISSN:
1354-1013
Page Range / eLocation ID:
p. 6080-6096
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. null (Ed.)
    Abstract. Canopy radiative transfer is the primary mechanism by which models relate vegetation composition and state to the surface energy balance, which is important to light- and temperature-sensitive plant processes as well as understanding land–atmosphere feedbacks.In addition, certain parameters (e.g., specific leaf area, SLA) that have an outsized influence on vegetation model behavior can be constrained by observations of shortwave reflectance, thus reducing model predictive uncertainty.Importantly, calibrating against radiative transfer outputs allows models to directly use remote sensing reflectance products without relying on highly derived products (such as MODIS leaf area index) whose assumptions may be incompatible with the target vegetation model and whose uncertainties are usually not well quantified.Here, we created the EDR model by coupling the two-stream representation of canopy radiative transfer in the Ecosystem Demography model version 2 (ED2) with a leaf radiative transfer model (PROSPECT-5) and a simple soil reflectance model to predict full-range, high-spectral-resolution surface reflectance that is dependent on the underlying ED2 model state.We then calibrated this model against estimates of hemispherical reflectance (corrected for directional effects) from the NASA Airborne Visible/Infrared Imaging Spectrometer (AVIRIS) and survey data from 54 temperate forest plots in the northeastern United States.The calibration significantly reduced uncertainty in model parameters related to leaf biochemistry and morphology and canopy structure for five plant functional types.Using a single common set of parameters across all sites, the calibrated model was able to accurately reproduce surface reflectance for sites with highly varied forest composition and structure.However, the calibrated model's predictions of leaf area index (LAI) were less robust, capturing only 46 % of the variability in the observations.Comparing the ED2 radiative transfer model with another two-stream soil–leaf–canopy radiative transfer model commonly used in remote sensing studies (PRO4SAIL) illustrated structural errors in the ED2 representation of direct radiation backscatter that resulted in systematic underestimation of reflectance.In addition, we also highlight that, to directly compare with a two-stream radiative transfer model like EDR, we had to perform an additional processing step to convert the directional reflectance estimates of AVIRIS to hemispherical reflectance (also known as “albedo”).In future work, we recommend that vegetation models add the capability to predict directional reflectance, to allow them to more directly assimilate a wide range of airborne and satellite reflectance products.We ultimately conclude that despite these challenges, using dynamic vegetation models to predict surface reflectance is a promising avenue for model calibration and validation using remote sensing data. 
    more » « less
  2. 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 to 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. 
    more » « less
  3. Abstract

    Soil moisture and evapotranspiration (ET) are important components of boreal forest hydrology that affect ecological processes and land‐atmosphere feedbacks. Future trends in soil moisture in particular are uncertain. Therefore, accurate modeling of these dynamics and understanding of concomitant sources of uncertainty are critical. Here, we conduct a global sensitivity analysis, Monte Carlo parameterization, and analysis of parameter uncertainty and its contribution to future soil moisture and ET uncertainty using a physically based ecohydrologic model in multiple boreal forest types. Soil and plant hydraulic parameters and LAI have the largest effects on simulated summer soil moisture at two contrasting sites. In future scenario simulations, the selection of parameters and global climate model (GCM) choice between two GCMs influence projected changes in soil moisture and ET about as much as the projected effects of climate change in the less sensitive GCM with a late‐century, high‐emissions scenario, though the relative effects of parameters, GCM, and climate vary among hydrologic variables and study sites. Saturated volumetric water content and sensitivity of stomatal conductance to vapor pressure deficit have the most statistically significant effects on change in ET and soil moisture, though there is considerable variability between sites and GCMs. The results of this study provide estimates of: (a) parameter importance and statistical significance for soil moisture modeling, (b) parameter values for physically based soil‐vegetation‐atmosphere transfer models in multiple boreal forest types, and (c) the contributions of uncertainty in these parameters to soil moisture and ET uncertainty in future climates.

     
    more » « less
  4. Abstract

    The biogeochemical processes of carbon (C), nitrogen (N), and phosphorous (P) are fully coupled in the Earth system, which shape the structure, functioning, and dynamics of terrestrial ecosystems. However, the representation of P cycle in terrestrial biosphere models (TBMs) is still in an early stage. Here we incorporated P processes and C‐N‐P interactions into the C‐N coupled Dynamic Land Ecosystem Model (DLEM‐CNP), which had a major feature of the ability in simulating the N and P colimitation on vegetation C assimilation. DLEM‐CNP was intensively calibrated and validated against daily or annual observations from four eddy covariance flux sites, two Hawaiian sites along a chronosequence of soils, and other 13 tropical forest sites. The results indicate that DLEM‐CNP significantly improved simulations of forest gross and net primary production (R2: 0.36–0.97, RMSE:1.1–1.49 g C m−2 year−1for daily GPP at eddy covariance flux sites;R2 = 0.92, RMSE = 176.7 g C m−2 year−1for annual NPP across 13 tropical forest sites). The simulations were also consistent with field observations in terms of biomass, leaf N:P ratio and plant response to fertilizer addition. A sensitivity analysis suggests that simulated results are reasonably robust against uncertainties in model parameter estimates and the model was very sensitive to parameters of P uptake. These results suggest that incorporating P processes and N‐P interaction into terrestrial biosphere models is of critical importance for accurately estimating C dynamics in tropical forests, particularly those P‐limited ones.

     
    more » « less
  5. Abstract

    Carbon use efficiency (CUE) represents how efficient a plant is at translating carbon gains through gross primary productivity (GPP) into net primary productivity (NPP) after respiratory costs (Ra). CUE varies across space with climate and species composition, but how CUE will respond to climate change is largely unknown due to uncertainty inRaat novel high temperatures. We use a plant physiological model validated against global CUE observations and LIDAR vegetation canopy height data and find that model‐predicted decreases in CUE are diagnostic of transitions from forests to shrubland at dry range edges. Under future climate scenarios, we show mean growing season CUE increases in core forested areas, but forest extent decreases at dry range edges, with substantial uncertainty in absolute CUE due to uncertainty inRa. Our results highlight that future forest resilience is nuanced and controlled by multiple competing mechanisms.

     
    more » « less