As temperatures continue rising, the direction, magnitude, and tempo of change in disturbance‐prone forests remain unresolved. Even forests long resilient to stand‐replacing fire face uncertain futures, and efforts to project changes in forest structure and composition are sorely needed to anticipate future forest trajectories. We simulated fire (incorporating fuels feedbacks) and forest dynamics on five landscapes spanning the Greater Yellowstone Ecosystem (GYE) to ask the following questions: (1) How and where are forest landscapes likely to change with 21st‐century warming and fire activity? (2) Are future forest changes gradual or abrupt, and do forest attributes change synchronously or sequentially? (3) Can forest declines be averted by mid‐21st‐century stabilization of atmospheric greenhouse gas (GHG) concentrations? We used the spatially explicit individual‐based forest model iLand to track multiple attributes (forest extent, stand age, tree density, basal area, aboveground carbon stocks, dominant forest types, species occupancy) through 2100 for six climate scenarios. Hot‐dry climate scenarios led to more fire, but stand‐replacing fire peaked in mid‐century and then declined even as annual area burned continued to rise. Where forest cover persisted, previously dense forests were converted to sparse young woodlands. Increased aridity and fire drove a ratchet of successive abrupt declines (i.e., multiple annual landscape‐level changes ≥20%) in tree density, basal area, and extent of older (>150 yr) forests, whereas declines in carbon stocks and mean stand age were always gradual. Forest changes were asynchronous across landscapes, but declines in stand structure always preceded reductions in forest extent and carbon stocks. Forest decline was most likely in less topographically complex landscapes dominated by fire‐sensitive tree species (Picea engelmannii, Abies lasiocarpa, Pinus contorta var. latifolia) and where fire resisters (Pseudotsuga menziesii var. glauca) were not already prevalent. If current GHG emissions continue unabated (RCP 8.5) and aridity increases, a suite of forest changes would transform the GYE, with cascading effects on biodiversity and myriad ecosystem services. However, stabilizing GHG concentrations by mid‐century (RCP 4.5) would slow the ratchet, moderating fire activity and dampening the magnitude and rate of forest change. Monitoring changes in forest structure may serve as an operational early warning indicator of impending forest decline.

Tree regeneration is a key process in forest dynamics, particularly in the context of forest resilience and climate change. Models are pivotal for assessing long‐term forest dynamics, and they have been in use for more than 50 years. However, there is a need to evaluate their capacity to accurately represent tree regeneration. We assess how well current models capture the overall abundance, species composition, and mortality of tree regeneration. Using 15 models built to capture long‐term forest dynamics at the stand, landscape, and global levels, we simulate tree regeneration at 200 sites representing large environmental gradients across Central Europe. The results are evaluated against extensive data from unmanaged forests. Most of the models overestimate recruitment levels, which is compensated only in some models by high simulated mortality rates in the early stages of individual‐tree dynamics. Simulated species diversity of recruitment generally matches observed ranges. Models simulating higher stand‐level species diversity do not feature higher species diversity in the recruitment layer. The effect of light availability on recruitment levels is captured better than the effects of temperature and soil moisture, but patterns are not consistent across models. Increasing complexity in the tree regeneration modules is not related to higher accuracy of simulated tree recruitment. Furthermore, individual model design is more important than scale (stand, landscape, and global) and approach (empirical and process‐based) for accurately capturing tree regeneration. Despite the mismatches between simulation results and data, it is remarkable that most models capture the essential features of the highly complex process of tree regeneration, while not having been parameterized with such data. We conclude that much can be gained by evaluating and refining the modeling of tree regeneration processes. This has the potential to render long‐term projections of forest dynamics under changing environmental conditions much more robust.

Mountain forests are plant diversity hotspots, but changing climate and increasing forest disturbances will likely lead to far‐reaching plant community change. Projecting future change, however, is challenging for forest understory plants, which respond to forest structure and composition as well as climate. Here, we jointly assessed the effects of both climate and forest change, including wind and bark beetle disturbances, using the process‐based simulation model iLand in a protected landscape in the northern Alps (Berchtesgaden National Park, Germany), asking: (1) How do understory plant communities respond to 21st‐century change in a topographically complex mountain landscape, representing a hotspot of plant species richness? (2) How important are climatic changes (i.e., direct climate effects) versus forest structure and composition changes (i.e., indirect climate effects and recovery from past land use) in driving understory responses at landscape scales? Stacked individual species distribution models fit with climate, forest, and soil predictors (248 species currently present in the landscape, derived from 150 field plots stratified by elevation and forest development, overall area under the receiving operator characteristic curve = 0.86) were driven with projected climate (RCP4.5 and RCP8.5) and modeled forest variables to predict plant community change. Nearly all species persisted in the landscape in 2050, but on average 8% of the species pool was lost by the end of the century. By 2100, landscape mean species richness and understory cover declined (−13% and −8%, respectively), warm‐adapted species increasingly dominated plant communities (i.e., thermophilization, +12%), and plot‐level turnover was high (62%). Subalpine forests experienced the greatest richness declines (−16%), most thermophilization (+17%), and highest turnover (67%), resulting in plant community homogenization across elevation zones. Climate rather than forest change was the dominant driver of understory responses. The magnitude of unabated 21st‐century change is likely to erode plant diversity in a species richness hotspot, calling for stronger conservation and climate mitigation efforts.

As 21st‐century climate and disturbance dynamics depart from historic baselines, ecosystem resilience is uncertain. Multiple drivers are changing simultaneously, and interactions among drivers could amplify ecosystem vulnerability to change. Subalpine forests in Greater Yellowstone (Northern Rocky Mountains, USA) were historically resilient to infrequent (100–300 year), severe fire. We sampled paired short‐interval (125‐year) post‐fire plots most recently burned between 1988 and 2018 to address two questions: (1) How do short‐interval fire, climate, topography, and distance to unburned live forest edge interact to affect post‐fire forest regeneration? (2) How do forest biomass and fuels vary following short‐interval versus long‐interval severe fires? Mean post‐fire live tree stem density was an order of magnitude lower following short‐interval versus long‐interval fires (3240 vs. 28,741 stems ha⁻¹, respectively). Differences between paired plots were amplified at longer distances to live forest edge. Surprisingly, warmer–drier climate was associated with higher seedling densities even after short‐interval fire, likely relating to regional variation in serotiny of lodgepole pine (Pinus contorta var. latifolia). Unlike conifers, density of aspen (Populus tremuloides), a deciduous resprouter, increased with short‐interval versus long‐interval fires (mean 384 vs. 62 stems ha⁻¹, respectively). Live biomass and canopy fuels remained low nearly 30 years after short‐interval fire, in contrast with rapid recovery after long‐interval fire, suggesting that future burn severity may be reduced for several decades following reburns. Short‐interval plots also had half as much dead woody biomass compared with long‐interval plots (60 vs. 121 Mg ha⁻¹), primarily due to the absence of large snags. Our results suggest differences in tree regeneration following short‐interval versus long‐interval fires will be especially pronounced where serotiny was high historically. Propagule limitation will also interact with short‐interval fires to diminish tree regeneration but lessen subsequent burn severity. Amplifying driver interactions are likely to threaten forest resilience under expected trajectories of a future fire.