Search

Vegetation distribution and composition are expected to change in northern high latitudes under rapid warming, which regulates ecosystem functions but remains challenging to predict. Vegetation change arises from the interplay of chronic climate trends such as warming and transient demographic processes of recruitment, growth, competition, and mortality. Most predictive models overlooked the role of demographic dynamics controlled by plant traits. Here, we simulate vegetation dynamics at the Kougarok Hillslope site in Alaska under historical and future climates using the E3SM Land Model coupled to the Functionally Assembled Terrestrial Simulator (ELM-FATES). To evaluate the roles of plant traits, we parameterize the model with 5,265 trait configurations representing diverse physiological and demographic strategies. Results show current modeled biomass, composition, and productivity are most sensitive to traits controlling photosynthetic capacity, carbon allocation, allometry, and phenology. Among all trait configurations, ∼5% reproduce in situ biomass and plant functional type (PFT) composition measured in 2016, that are indistinguishable from these two observed ecosystem states. Notably, these same trait configurations produce diverging biomass, composition, and productivity under future climate, where the uncertainty attributable to traits is twice the change attributable to climate change. The variation of projected productivity arises from emerging PFT composition under novel climate regimes, primarily explained by traits controlling cold-induced mortality, recruitment, and allometry. Our findings highlight the importance and uncertainty of demographic dynamics and its interaction with climate change in shaping Arctic vegetation change. Improved model predictions will likely benefit from explicit consideration of vegetation demography and better constraints of critical traits.

Understanding forest water limitation during droughts within a warming climate is essential for accurate predictions of forest-climate interactions. In hyperdiverse ecosystems like the Amazon forest, the mechanisms shaping hysteresis patterns in transpiration relative to environmental factors are not well understood. From this perspective, we investigated these dynamics by conducting in situ leaf-level measurements throughout and after the 2015 El Niño-Southern Oscillation (ENSO) drought. Our findings indicate a substantial increase in the hysteresis area (Harea) among transpiration (E), vapor pressure deficit (VPD), and stomatal conductance (gs) at canopy level during the ENSO peak, attributed to both temporal lag and differences in magnitude between gs and VPD peaks. Specifically, the canopy species Pouteria anomala exhibited an increased Harea, due to earlier maximum gs rates leading to a greater temporal lag with VPD compared to the post-drought period. Additionally, leaf water potential (ΨL) and canopy temperature (Tcanopy) showed larger Harea during the ENSO peak compared to post-drought conditions across all studied species, suggesting that stomatal closure, particularly during the afternoon, acts to minimize water loss and may explain the counterclockwise hysteresis observed between ΨL and Tcanopy. The pronounced Harea during the drought points to a potential imbalance between water supply and demand, underlining the role of stomatal behavior of isohydric species in response to drought.

Biomass allocation in plants is the foundation for understanding dynamics in ecosystem carbon balance, species competition, and plant-environment interactions. However, existing work on plant allometry has mainly focused on trees, with fewer studies having developed allometric equations for grasses. Grasses with different life histories can vary in their carbon investment by prioritizing the growth of specific organs to survive, outcompete co-occurring plants, and ensure population persistence. Further, because grasses are important fuels for wildfire, the lack of grass allocation data adds uncertainty to process-based models that relate plant physiology to wildfire dynamics. To fill this gap, we conducted a greenhouse experiment with 11 common California grasses varying in photosynthetic pathway and growth form. We measured plant sizes and harvested above- and belowground biomass throughout the life cycle of annual species, while for the establishment stage of perennial grasses to quantify allometric relationships for leaf, stem, and root biomass, as well as plant height and canopy area. We used basal diameter as a reference measure of plant size. Overall, basal diameter is the best predictor for leaf and stem biomass, height, and canopy area. Including height as another predictor can improve model accuracy in predicting leaf and stem biomass and canopy area. Fine root biomass is a function of leaf biomass alone. Species vary in their allometric relationships, with most variation occurring for plant height, canopy area, and stem biomass. We further explored potential trade-offs in biomass allocation across species between leaf and fine root, leaf and stem, and allocation to reproduction. Consistent with our expectation, we found that fast-growing plants allocated a greater fraction to reproduction. Additionally, plant height and specific leaf area negatively influenced the leaf-to-stem ratio. However, contrary to our hypothesis, there were no differences in root-to-leaf ratio between perennial and annual or C4 and C3 plants. Our study provides species-specific and functional-type-specific allometry equations for both above- and belowground organs of 11 common California grass species, enabling nondestructive biomass assessment in California grasslands. These allometric relationships and trade-offs in carbon allocation across species can improve ecosystem model predictions of grassland species interactions and environmental responses through differences in morphology.

We present a representation of nitrogen and phosphorus cycling in the Functionally Assembled Terrestrial Ecosystem Simulator, a demographic vegetation model within the Energy Exascale Earth System land model. This representation is modular, and designed to allow testing of multiple hypothetical approaches for carbon-nutrient coupling in plants. Novel model hypotheses introduced in this work include, (a) the controls on plant acquisition of aqueous mineralized nutrients in the soil and (b) fairly straight forward methods of allocating nutrients to specific plant organs and their losses through live plant turnover as well as litter fluxes generated through plant mortality. This combines the new with pre-existing hypotheses (such as nitrogen fixation and soil decomposition) into a system that can accommodate plant-soil dynamics for a large number of size- and functional-type-resolved plant cohorts within a time-since-disturbance-resolved ecosystem. Root uptake of nutrients is governed by fine root biomass, and plants vary in their fine root biomass allocation in order to balance carbon and nutrient limitations to growth. We test the sensitivity of the model to a wide range of parameter variations and structural representations, and in the context of observations at Barro Colorado Island, Panama. A key model prediction is that plants in the high-light-availability canopy positions allocate more carbon to fine roots than plants in low-light understory environments, given the widely different carbon versus nutrient constraints of these two niches within a given ecosystem. This model provides a basis for exploring carbon-nutrient coupling with vegetation demography within Earth system models.

Tropical forest dynamics play a crucial role in the global carbon, water, and energy cycles. However, realistically simulating the dynamics of competition and coexistence between different plant functional types (PFTs) in tropical forests remains a significant challenge. This study aims to improve the modeling of PFT coexistence in the Functionally Assembled Terrestrial Ecosystem Simulator (FATES), a vegetation demography model implemented in the Energy Exascale Earth System Model (E3SM) land model (ELM), ELM-FATES. Specifically, we explore (1) whether plant trait relationships established from field measurements can constrain ELM-FATES simulations and (2) whether machine learning (ML)-based surrogate models can emulate the complex ELM-FATES model and optimize parameter selections to improve PFT coexistence modeling. We conducted three ensembles of ELM-FATES experiments at a tropical forest site near Manaus, Brazil. By comparing the ensemble experiments without (Exp-CTR) and with (Exp-OBS) consideration of observed trait relationships, we found that accounting for these relationships slightly improves the simulations of water, energy, and carbon variables when compared to observations but degrades the simulation of PFT coexistence. Using ML-based surrogate models trained on Exp-CTR, we optimized the trait parameters in ELM-FATES and conducted another ensemble of experiments (Exp-ML) with these optimized parameters. The proportion of PFT coexistence experiments significantly increased from 21 % in Exp-CTR to 73 % in Exp-ML. After filtering the experiments that allow for PFT coexistence to agree with observations (within 15 % tolerance), 33 % of the Exp-ML experiments were retained, which is a significant improvement compared to the 1.4 % in Exp-CTR. Exp-ML also accurately reproduces the annual means and seasonal variations in water, energy, and carbon fluxes and the field inventory of aboveground biomass. This study represents a reproducible method that utilizes machine learning to identify parameter values that improve model fidelity against observations and PFT coexistence in vegetation demography models for diverse ecosystems. Our study also suggests the need for new mechanisms to enhance the robust simulation of coexisting plants in ELM-FATES and has significant implications for modeling the response and feedbacks of ecosystem dynamics to climate change.

We explore the response of the Earth’s coupled climate and carbon system to an idealized sequential addition and removal of CO2 to the atmosphere, following a symmetric and continuous emissions pathway, in contrast to the discontinuous emissions pathways that have largely informed our understanding of the climate response to net zero and net negative emissions to date. We find, using both an Earth system model and an ensemble of simple climate model realizations, that warming during the emissions reduction and negative emissions phases is defined by a combination of a proportionality of warming to cumulative emissions characterized by the transient climate response to emissions (TCRE), and a deviation from that proportionality that is governed by the zero emissions commitment (ZEC). About half of the ZEC is realized before reaching zero emissions, and the ZEC thus also controls the timing between peak cumulative CO2 emissions and peak temperature, such that peak temperature may occur before peak cumulative emissions if ZEC is negative, underscoring the importance of ZEC in climate policies aimed to limit peak warming. Thus we argue that ZEC is better defined as the committed warming relative to the expected TCRE proportionality, rather than as the additional committed warming that will occur after reaching net zero CO2 emissions. Once established, the combined TCRE and ZEC relationship holds almost to complete removal of prior cumulative CO2 emissions. As cumulative CO2 emissions approach zero through negative CO2 emissions, CO2 concentrations drop below preindustrial values, while residual long-term climate change continues, governed by multicentennial dynamical processes.

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 model 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.

The growth and survival of individual trees determine the physical structure of a forest with important consequences for forest function. However, given the diversity of tree species and forest biomes, quantifying the multitude of demographic strategies within and across forests and the way that they translate into forest structure and function remains a significant challenge. Here, we quantify the demographic rates of 1961 tree species from temperate and tropical forests and evaluate how demographic diversity (DD) and demographic composition (DC) differ across forests, and how these differences in demography relate to species richness, aboveground biomass (AGB), and carbon residence time. We find wide variation in DD and DC across forest plots, patterns that are not explained by species richness or climate variables alone. There is no evidence that DD has an effect on either AGB or carbon residence time. Rather, the DC of forests, specifically the relative abundance of large statured species, predicted both biomass and carbon residence time. Our results demonstrate the distinct DCs of globally distributed forests, reflecting biogeography, recent history, and current plot conditions. Linking the DC of forests to resilience or vulnerability to climate change, will improve the precision and accuracy of predictions of future forest composition, structure, and function.

Scientific communities are increasingly publishing data to evaluate, accredit, and build on published research. However, guidelines for curating data for publication are sparse for model-related research, limiting the usability of archived simulation data. In particular, there are no established guidelines for archiving data related to terrestrial models that simulate land processes and their coupled interactions with climate. Terrestrial modelers have a unique set of challenges when publishing data due to the diversity of scientific domains, research questions, and the types and scales of simulations. Researchers in the U.S. Department of Energy’s (DOE) projects use a variety of multiscale models to advance robust predictions of terrestrial and subsurface ecosystem processes. Here, we synthesize archiving needs for data associated with different DOE models, and provide guidelines for publishing terrestrial model data components following FAIR (Findable, Accessible, Interoperable, Reusable) principles. The guidelines recommend archiving model inputs and testing data used in final simulation runs along with associated codes, workflow scripts, and metadata in public repositories. Researchers should consider archiving model outputs if they are within the storage limits of the repository. We also provide considerations for how to bundle files into different data publications with citable digital object identifiers. Finally, we identify repository features and tools that would enable storage and reuse of model data. Given the diversity of DOE terrestrial models, these guidelines are transferable to other model types and will enable efficient reuse of simulation data for purposes such as model intercomparisons, initialization, benchmarking, synthesis, and comparisons with field observations.

Future climate projections from Earth system models (ESMs) typically focus on the timescale of this century. We use a set of five ESMs and one Earth system model of intermediate complexity (EMIC) to explore the dynamics of the Earth's climate and carbon cycles under contrasting emissions trajectories beyond this century to the year 2300. The trajectories include a very-high-emissions, unmitigated fossil-fuel-driven scenario, as well as a mitigation scenario that diverges from the first scenario after 2040 and features an "overshoot", followed by a decrease in atmospheric CO2 concentrations by means of large net negative CO2 emissions. In both scenarios and for all models considered here, the terrestrial system switches from being a net sink to either a neutral state or a net source of carbon, though for different reasons and centered in different geographic regions, depending on both the model and the scenario. The ocean carbon system remains a sink, albeit weakened by carbon cycle feedbacks, in all models under the high-emissions scenario and switches from sink to source in the overshoot scenario. The global mean temperature anomaly is generally proportional to cumulative carbon emissions, with a deviation from proportionality in the overshoot scenario that is governed by the zero emissions commitment. Additionally, 23rd century warming continues after the cessation of carbon emissions in several models in the high-emissions scenario and in one model in the overshoot scenario. While ocean carbon cycle responses qualitatively agree in both globally integrated and zonal mean dynamics in both scenarios, the land models qualitatively disagree in zonal mean dynamics, in the relative roles of vegetation and soil in driving C fluxes, in the response of the sink to CO2, and in the timing of the sink-source transition, particularly in the high-emissions scenario. The lack of agreement among land models on the mechanisms and geographic patterns of carbon cycle feedbacks, alongside the potential for lagged physical climate dynamics to cause warming long after CO2 concentrations have stabilized, points to the possibility of surprises in the climate system beyond the 21st century time horizon, even under relatively mitigated global warming scenarios, which should be taken into consideration when setting global climate policy. Copyright:

Arctic warming is affecting snow cover and soil hydrology, with consequences for carbon sequestration in tundra ecosystems. The scarcity of observations in the Arctic has limited our understanding of the impact of covarying environmental drivers on the carbon balance of tundra ecosystems. In this study, we address some of these uncertainties through a novel record of 119 site-years of summer data from eddy covariance towers representing dominant tundra vegetation types located on continuous permafrost in the Arctic. Here we found that earlier snowmelt was associated with more tundra net CO2 sequestration and higher gross primary productivity (GPP) only in June and July, but with lower net carbon sequestration and lower GPP in August. Although higher evapotranspiration (ET) can result in soil drying with the progression of the summer, we did not find significantly lower soil moisture with earlier snowmelt, nor evidence that water stress affected GPP in the late growing season. Our results suggest that the expected increased CO2 sequestration arising from Arctic warming and the associated increase in growing season length may not materialize if tundra ecosystems are not able to continue sequestering CO2 later in the season.

Tropical forest photosynthesis can decline at high temperatures due to (1) biochemical responses to increasing temperature and (2) stomatal responses to increasing vapor pressure deficit (VPD), which is associated with increasing temperature. It is challenging to disentangle the influence of these two mechanisms on photosynthesis in observations, because temperature and VPD are tightly correlated in tropical forests. Nonetheless, quantifying the relative strength of these two mechanisms is essential for understanding how tropical gross primary production (GPP) will respond to climate change, because increasing atmospheric CO2 concentration may partially offset VPD‐driven stomatal responses, but is not expected to mitigate the effects of temperature‐driven biochemical responses. We used two terrestrial biosphere models to quantify how physiological process assumptions (photosynthetic temperature acclimation and plant hydraulic stress) and functional traits (e.g., maximum xylem conductivity) influence the relative strength of modeled temperature versus VPD effects on light‐saturated GPP at an Amazonian forest site, a seasonally dry tropical forest site, and an experimental tropical forest mesocosm. By simulating idealized climate change scenarios, we quantified the divergence in GPP predictions under model configurations with stronger VPD effects compared with stronger direct temperature effects. Assumptions consistent with stronger direct temperature effects resulted in larger GPP declines under warming, while assumptions consistent with stronger VPD effects resulted in more resilient GPP under warming. Our findings underscore the importance of quantifying the role of direct temperature and indirect VPD effects for projecting the resilience of tropical forests in the future, and demonstrate that the relative strength of temperature versus VPD effects in models is highly sensitive to plant functional parameters and structural assumptions about photosynthetic temperature acclimation and plant hydraulics. [ABSTRACT FROM AUTHOR]

Summary: Terrestrial biosphere models (TBMs) include the representation of vertical gradients in leaf traits associated with modeling photosynthesis, respiration, and stomatal conductance. However, model assumptions associated with these gradients have not been tested in complex tropical forest canopies.We compared TBM representation of the vertical gradients of key leaf traits with measurements made in a tropical forest in Panama and then quantified the impact of the observed gradients on simulated canopy‐scale CO2 and water fluxes.Comparison between observed and TBM trait gradients showed divergence that impacted canopy‐scale simulations of water vapor and CO2 exchange. Notably, the ratio between the dark respiration rate and the maximum carboxylation rate was lower near the ground than at the top‐of‐canopy, leaf‐level water‐use efficiency was markedly higher at the top‐of‐canopy, and the decrease in maximum carboxylation rate from the top‐of‐canopy to the ground was less than TBM assumptions.The representation of the gradients of leaf traits in TBMs is typically derived from measurements made within‐individual plants, or, for some traits, assumed constant due to a lack of experimental data. Our work shows that these assumptions are not representative of the trait gradients observed in species‐rich, complex tropical forests. [ABSTRACT FROM AUTHOR]

Extreme droughts are a major determinant of ecosystem disturbance, which impact plant communities and feed back to climate change through changes in plant functioning. However, the complex relationships between above- and belowground plant hydraulic traits, and their role in governing plant responses to drought, are not fully understood. In this study, we use a plant hydraulics model, FATES-Hydro, to investigate ecosystem responses to the 2012-2015 California drought, in comparison with observations, for a site in the southern Sierra Nevada that experienced widespread tree mortality during this drought. We conduct a sensitivity analysis to explore how different plant water sourcing and hydraulic strategies lead to differential responses during normal and drought conditions. The analysis shows that: 1) deep roots that sustain productivity through the dry season are needed for the model to capture observed seasonal cycles of ET and GPP in normal years, and that deep- rooted strategies are nonetheless subject to large reductions in ET and GPP when the deep soil reservoir is depleted during extreme droughts, in agreement with observations. 2) risky stomatal strategies lead to greater productivity during normal years as compared to safer stomatal control, but lead to high risk of xylem embolism during the 2012-2015 drought. 3) for a given stand density, the stomatal and xylem traits have a stronger impact on plant water status than on ecosystem level fluxes. Our study reveals the importance of resolving plant water sourcing strategies in order to represent drought impacts on plants, and consequent feedbacks, in models. [ABSTRACT FROM AUTHOR]

Quantifying the responses of forest disturbances to climate warming is critical to our understanding of carbon cycles and energy balances of the Earth system. The impact of warming on bark beetle outbreaks is complex as multiple drivers of these events may respond differently to warming. Using a novel model of bark beetle biology and host tree interactions, we assessed how contemporary warming affected western pine beetle (Dendroctonus brevicomis) populations and mortality of its host, ponderosa pine (Pinus ponderosa), during an extreme drought in the Sierra Nevada, California, United States. When compared with the field data, our model captured the western pine beetle flight timing and rates of ponderosa pine mortality observed during the drought. In assessing the influence of temperature on western pine beetles, we found that contemporary warming increased the development rate of the western pine beetle and decreased the overwinter mortality rate of western pine beetle larvae leading to increased population growth during periods of lowered tree defense. We attribute a 29.9% (95% CI: 29.4%–30.2%) increase in ponderosa pine mortality during drought directly to increases in western pine beetle voltinism (i.e., associated with increased development rates of western pine beetle) and, to a much lesser extent, reductions in overwintering mortality. These findings, along with other studies, suggest each degree (°C) increase in temperature may have increased the number of ponderosa pine killed by upwards of 35%–40% °C−1 if the effects of compromised tree defenses (15%–20%) and increased western pine beetle populations (20%) are additive. Due to the warming ability to considerably increase mortality through the mechanism of bark beetle populations, models need to consider climate's influence on both host tree stress and the bark beetle population dynamics when determining future levels of tree mortality. [ABSTRACT FROM AUTHOR]

Terrestrial biosphere models (TBMs) include the representation of vertical gradients in leaf traits associated with modeling photosynthesis, respiration, and stomatal conductance. However, model assumptions associated with these gradients have not been tested in complex tropical forest canopies. We compared TBM representation of the vertical gradients of key leaf traits with measurements made in a tropical forest in Panama and then quantified the impact of the observed gradients on simulated canopy-scale CO 2 and water fluxes. Comparison between observed and TBM trait gradients showed divergence that impacted canopy-scale simulations of water vapor and CO 2 exchange. Notably, the ratio between the dark respiration rate and the maximum carboxylation rate was lower near the ground than at the top-of-canopy, leaf-level water-use efficiency was markedly higher at the top-of-canopy, and the decrease in maximum carboxylation rate from the top-of-canopy to the ground was less than TBM assumptions. The representation of the gradients of leaf traits in TBMs is typically derived from measurements made within-individual plants, or, for some traits, assumed constant due to a lack of experimental data. Our work shows that these assumptions are not representative of the trait gradients observed in species-rich, complex tropical forests., (© 2023 The Authors New Phytologist © 2023 New Phytologist Foundation This article has been contributed to by U.S. Government employees and their work is in the public domain in the USA.)