skip to main content


Title: Increasing Functional Diversity in a Global Land Surface Model Illustrates Uncertainties Related to Parameter Simplification
Abstract

Simulations of the land surface carbon cycle typically compress functional diversity into a small set of plant functional types (PFT), with parameters defined by the average value of measurements of functional traits. In most earth system models, all wild plant life is represented by between five and 14 PFTs and a typical grid cell (≈100 × 100 km) may contain a single PFT. Model logic applied to this coarse representation of ecological functional diversity provides a reasonable proxy for the carbon cycle, but does not capture the non‐linear influence of functional traits on productivity. Here we show through simulations using the Energy Exascale Land Surface Model in 15 diverse terrestrial landscapes, that better accounting for functional diversity markedly alters predicted total carbon uptake. The shift in carbon uptake is as great as 30% and 10% in boreal and tropical regions, respectively, when compared to a single PFT parameterized with the trait means. The traits that best predict gross primary production vary based on vegetation phenology, which broadly determines where traits fall within the global distribution. Carbon uptake is more closely associated with specific leaf area for evergreen PFTs and the leaf carbon to nitrogen ratio in deciduous PFTs.

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

    Tropical forests are a key determinant of the functioning of the Earth system, but remain a major source of uncertainty in carbon cycle models and climate change projections. In this study, we present an updated land model (LM3PPA‐TV) to improve the representation of tropical forest structure and dynamics in Earth system models (ESMs). The development and parameterization of LM3PPA‐TV drew on extensive datasets on tropical tree traits and long‐term field censuses from Barro Colorado Island (BCI), Panama. The model defines a new plant functional type (PFT) based on the characteristics of shade‐tolerant, tropical tree species, implements a new growth allocation scheme based on realistic tree allometries, incorporates hydraulic constraints on biomass accumulation, and features a new compartment for tree branches and branch fall dynamics. Simulation experiments reproduced observed diurnal and seasonal patterns in stand‐level carbon and water fluxes, as well as mean canopy and understory tree growth rates, tree size distributions, and stand‐level biomass on BCI. Simulations at multiple sites captured considerable variation in biomass and size structure across the tropical forest biome, including observed responses to precipitation and temperature. Model experiments suggested a major role of water limitation in controlling geographic variation forest biomass and structure. However, the failure to simulate tropical forests under extreme conditions and the systematic underestimation of forest biomass in Paleotropical locations highlighted the need to incorporate variation in hydraulic traits and multiple PFTs that capture the distinct floristic composition across tropical domains. The continued pressure on tropical forests from global change demands models which are able to simulate alternative successional pathways and their pace to recovery. LM3PPA‐TV provides a tool to investigate geographic variation in tropical forests and a benchmark to continue improving the representation of tropical forests dynamics and their carbon storage potential in ESMs.

     
    more » « less
  2. 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
  3. Abstract

    The spatio‐temporal variation of stomatal conductance directly regulates photosynthesis, water partitioning, and biosphere‐atmosphere interactions. While many studies have focused on stomatal response to stresses, the spatial variation of unstressed stomatal conductance remains poorly determined, and is usually characterized in land surface models (LSMs) simply based on plant functional type (PFT). Here, we derived unstressed stomatal conductance at the ecosystem‐scale using observations from 115 global FLUXNET sites. When aggregated by PFTs, the across‐PFT pattern was highly consistent with the parameterizations of LSMs. However, PFTs alone captured only 17% of the variation in unstressed stomatal conductance across sites. Within the same PFT, unstressed stomatal conductance was negatively related to climate dryness and canopy height, which explained 45% of the total spatial variation. Our results highlight the importance of plant‐environment interactions in shaping stomatal traits. The trait‐environment relationship established here provides an empirical approach for improved parameterizations of stomatal conductance in LSMs.

     
    more » « less
  4. Summary

    Process‐based vegetation models attempt to represent the wide range of trait variation in biomes by grouping ecologically similar species into plant functional types (PFTs). This approach has been successful in representing many aspects of plant physiology and biophysics but struggles to capture biogeographic history and ecological dynamics that determine biome boundaries and plant distributions. Grass‐dominated ecosystems are broadly distributed across all vegetated continents and harbour large functional diversity, yet most Land Surface Models (LSMs) summarise grasses into two generic PFTs based primarily on differences between temperate C3grasses and (sub)tropical C4grasses. Incorporation of species‐level trait variation is an active area of research to enhance the ecological realism of PFTs, which form the basis for vegetation processes and dynamics in LSMs. Using reported measurements, we developed grass functional trait values (physiological, structural, biochemical, anatomical, phenological, and disturbance‐related) of dominant lineages to improve LSM representations. Our method is fundamentally different from previous efforts, as it uses phylogenetic relatedness to create lineage‐based functional types (LFTs), situated between species‐level trait data and PFT‐level abstractions, thus providing a realistic representation of functional diversity and opening the door to the development of new vegetation models.

     
    more » « less
  5. Abstract Widespread changes in the distribution and abundance of plant functional types (PFTs) are occurring in Arctic and boreal ecosystems due to the intensification of disturbances, such as fire, and climate-driven vegetation dynamics, such as tundra shrub expansion. To understand how these changes affect boreal and tundra ecosystems, we need to first quantify change for multiple PFTs across recent years. While landscape patches are generally composed of a mixture of PFTs, most previous moderate resolution (30 m) remote sensing analyses have mapped vegetation distribution and change within land cover categories that are based on the dominant PFT; or else the continuous distribution of one or a few PFTs, but for a single point in time. Here we map a 35 year time-series (1985–2020) of top cover (TC) for seven PFTs across a 1.77 × 10 6 km 2 study area in northern and central Alaska and northwestern Canada. We improve on previous methods of detecting vegetation change by modeling TC, a continuous measure of plant abundance. The PFTs collectively include all vascular plants within the study area as well as light macrolichens, a nonvascular class of high importance to caribou management. We identified net increases in deciduous shrubs (66 × 10 3 km 2 ), evergreen shrubs (20 × 10 3 km 2 ), broadleaf trees (17 × 10 3 km 2 ), and conifer trees (16 × 10 3 km 2 ), and net decreases in graminoids (−40 × 10 3 km 2 ) and light macrolichens (−13 × 10 3 km 2 ) over the full map area, with similar patterns across Arctic, oroarctic, and boreal bioclimatic zones. Model performance was assessed using spatially blocked, nested five-fold cross-validation with overall root mean square errors ranging from 8.3% to 19.0%. Most net change occurred as succession or plant expansion within areas undisturbed by recent fire, though PFT TC change also clearly resulted from fire disturbance. These maps have important applications for assessment of surface energy budgets, permafrost changes, nutrient cycling, and wildlife management and movement analysis. 
    more » « less