Caribbean tropical forests are subject to hurricane disturbances of great variability. In addition to natural storm incongruity, climate change can alter storm formation, duration, frequency, and intensity. This model-based investigation assessed the impacts of multiple storms of different intensities and occurrence frequencies on the long-term dynamics of subtropical dry forests in Puerto Rico. Using the previously validated individual-based gap model ZELIG-TROP, we developed a new hurricane damage routine and parameterized it with site- and species-specific hurricane effects. A baseline case with the reconstructed historical hurricane regime represented the control condition. Ten treatment cases, reflecting plausible shifts in hurricane regimes, manipulated both hurricane return time (i.e. frequency) and hurricane intensity. The treatment-related change in carbon storage and fluxes were reported as changes in aboveground forest biomass (AGB), net primary productivity (NPP), and in the aboveground carbon partitioning components, or annual carbon accumulation (ACA). Increasing the frequency of hurricanes decreased aboveground biomass by between 5% and 39%, and increased NPP between 32% and 50%. Decadal-scale biomass fluctuations were damped relative to the control. In contrast, increasing hurricane intensity did not create a large shift in the long-term average forest structure, NPP, or ACA from that of historical hurricane regimes, but produced large fluctuations in biomass. Decreasing both the hurricane intensity and frequency by 50% produced the highest values of biomass and NPP. For the control scenario and with increased hurricane intensity, ACA was negative, which indicated that the aboveground forest components acted as a carbon source. However, with an increase in the frequency of storms or decreased storms, the total ACA was positive due to shifts in leaf production, annual litterfall, and coarse woody debris inputs, indicating a carbon sink into the forest over the long-term. The carbon loss from each hurricane event, in all scenarios, always recovered over sufficient time. Our results suggest that subtropical dry forests will remain resilient to hurricane disturbance. However carbon stocks will decrease if future climates increase hurricane frequency by 50% or more.

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.

Frost is damaging to plants when air temperature drops below their tolerance threshold. The set of mechanisms used by cold-tolerant plants to withstand freezing is called “hardening” and typically take place in autumn to protect against winter damage. The recent incorporation of a hardening scheme in the demographic vegetation model FATES opens up the possibility to investigate frost mortality to vegetation. Previously, the hardening scheme was used to improve hydraulic processes in cold-tolerant plants. In this study, we expand upon the existing hardening scheme by implementing hardiness-dependent frost mortality into CLM5.0-FATES to study the impacts of frost on vegetation in temperate and boreal sites from 1950 to 2015. Our results show that the original freezing mortality approach of FATES, where each plant type had a fixed freezing tolerance threshold—an approach common to many other dynamic vegetation models, was restricted to predicting plant type distribution. The main results emerging from the new scheme are a high autumn and spring frost mortality, especially at colder sites, and increasing mid-winter frost mortality due to global warming, especially at warmer sites. We demonstrate that the new frost scheme is a major step forward in dynamically representing vegetation in ESMs by for the first time including a level of frost tolerance that is responding to the environment and includes some level of cost (implicitly) and benefit. By linking hardening and frost mortality in a land surface model, we open new ways to explore the impact of frost events in the context of global warming.

Arctic shrub expansion alters carbon budgets, albedo, and warming rates in high latitudes but remains challenging to predict due to unclear underlying controls. Observational studies and models typically use relationships between observed shrub presence and current environmental suitability (bioclimate and topography) to predict shrub expansion, while omitting shrub demographic processes and non-stationary response to changing climate. Here, we use high-resolution satellite imagery across Alaska and western Canada to show that observed shrub expansion has not been controlled by environmental suitability during 1984-2014, but can only be explained by considering seed dispersal and fire. These findings provide the impetus for better observations of recruitment and for incorporating currently underrepresented processes of seed dispersal and fire in land models to project shrub expansion and climate feedbacks. Integrating these dynamic processes with projected fire extent and climate, we estimate shrubs will expand into 25% of the non-shrub tundra by 2100, in contrast to 39% predicted based on increasing environmental suitability alone. Thus, using environmental suitability alone likely overestimates and misrepresents shrub expansion pattern and its associated carbon sink.

Forest disturbance and regrowth are key processes in forest dynamics, but detailed information on these processes is difficult to obtain in remote forests such as the Amazon. We used chronosequences of Landsat satellite imagery (Landsat 5 Thematic Mapper and Landsat 7 Enhanced Thematic Mapper Plus) to determine the sensitivity of surface reflectance from all spectral bands to windthrow, clear-cut, and clear-cut and burned (cut + burn) and their successional pathways of forest regrowth in the Central Amazon. We also assessed whether the forest demography model Functionally Assembled Terrestrial Ecosystem Simulator (FATES) implemented in the Energy Exascale Earth System Model (E3SM) Land Model (ELM), ELM-FATES, accurately represents the changes for windthrow and clear-cut. The results show that all spectral bands from the Landsat satellites were sensitive to the disturbances but after 3 to 6 years only the near-infrared (NIR) band had significant changes associated with the successional pathways of forest regrowth for all the disturbances considered. In general, the NIR values decreased immediately after disturbance, increased to maximum values with the establishment of pioneers and early successional tree species, and then decreased slowly and almost linearly to pre-disturbance conditions with the dynamics of forest succession. Statistical methods predict that NIR values will return to pre-disturbance values in about 39, 36, and 56 years for windthrow, clear-cut, and cut + burn disturbances, respectively. The NIR band captured the observed, and different, successional pathways of forest regrowth after windthrow, clear-cut, and cut + burn. Consistent with inferences from the NIR observations, ELM-FATES predicted higher peaks of biomass and stem density after clear-cuts than after windthrows. ELM-FATES also predicted recovery of forest structure and canopy coverage back to pre-disturbance conditions in 38 years after windthrows and 41 years after clear-cut. The similarity of ELM-FATES predictions of regrowth patterns after windthrow and clear-cut to those of the NIR results suggests the NIR band can be used to benchmark forest regrowth in ecosystem models. Our results show the potential of Landsat imagery data for mapping forest regrowth from different types of disturbances, benchmarking, and the improvement of forest regrowth 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.

Climatic extreme events are expected to occur more frequently in the future, increasing the likelihood of unprecedented climate extremes (UCEs) or record-breaking events. UCEs, such as extreme heatwaves and droughts, substantially affect ecosystem stability and carbon cycling by increasing plant mortality and delaying ecosystem recovery. Quantitative knowledge of such effects is limited due to the paucity of experiments focusing on extreme climatic events beyond the range of historical experience. Here, we present a road map of how dynamic vegetation demographic models (VDMs) can be used to investigate hypotheses surrounding ecosystem responses to one type of UCE: unprecedented droughts. As a result of nonlinear ecosystem responses to UCEs that are qualitatively different from responses to milder extremes, we consider both biomass loss and recovery rates over time by reporting a time-integrated carbon loss as a result of UCE, relative to the absence of drought. Additionally, we explore how unprecedented droughts in combination with increasing atmospheric CO2 and/or temperature may affect ecosystem stability and carbon cycling. We explored these questions using simulations of pre-drought and post-drought conditions at well-studied forest sites using well-tested models (ED2 and LPJ-GUESS). The severity and patterns of biomass losses differed substantially between models. For example, biomass loss could be sensitive to either drought duration or drought intensity depending on the model approach. This is due to the models having different, but also plausible, representations of processes and interactions, highlighting the complicated variability of UCE impacts that still need to be narrowed down in models. Elevated atmospheric CO2 concentrations (eCO2) alone did not completely buffer the ecosystems from carbon losses during UCEs in the majority of our simulations. Our findings highlight the consequences of differences in process formulations and uncertainties in models, most notably related to availability in plant carbohydrate storage and the diversity of plant hydraulic schemes, in projecting potential ecosystem responses to UCEs. We provide a summary of the current state and role of many model processes that give way to different underlying hypotheses of plant responses to UCEs, reflecting knowledge gaps which in future studies could be tested with targeted field experiments and an iterative modeling-experimental conceptual framework.

Tropical forests absorb large amounts of atmospheric CO₂ through photosynthesis but elevated temperatures suppress this absorption and promote monoterpene emissions. Using ¹³ CO₂ labeling, here we show that monoterpene emissions from tropical leaves derive from recent photosynthesis and demonstrate distinct temperature optima for five groups (Groups 1-5), potentially corresponding to different enzymatic temperature-dependent reaction mechanisms within β-ocimene synthases. As diurnal and seasonal leaf temperatures increased during the Amazonian 2015 El Niño event, leaf and landscape monoterpene emissions showed strong linear enrichments of β-ocimenes (+4.4% °C^-1 ) at the expense of other monoterpene isomers. The observed inverse temperature response of α-pinene (-0.8% °C^-1 ), typically assumed to be the dominant monoterpene with moderate reactivity, was not accurately simulated by current global emission models. Given that β-ocimenes are highly reactive with respect to both atmospheric and biological oxidants, the results suggest that highly reactive β-ocimenes may play important roles in the thermotolerance of photosynthesis by functioning as effective antioxidants within plants and as efficient atmospheric precursors of secondary organic aerosols. Thus, monoterpene composition may represent a new sensitive 'thermometer' of leaf oxidative stress and atmospheric reactivity, and therefore a new tool in future studies of warming impacts on tropical biosphere-atmosphere carbon-cycle feedbacks.