Understanding the impacts of changing climate and disturbance regimes on forest ecosystems is greatly aided by the use of process-based models. Such models simulate processes based on first principles of ecology, which requires parameterization. Parameterization is an important step in model development and application, defining the characteristics of trees and their responses to the environment, i.e., their traits. For species-specific models, parameterization is usually done at the level of individual species. Parameterization is indispensable for accurately modeling demographic processes, including growth, mortality, and regeneration of trees, along with their intra- and inter-specific interactions. As it is time-demanding to compile the parameters required to simulate forest ecosystems in complex models, simulations are often restricted to the most common tree species, genera, or plant-functional types. Yet, as tree species composition might change in the future, it is important to account for a broad range of species and their individual responses to drivers of change explicitly in simulations. Thus, species-specific parameterization is a critical task for making accurate projections about future forest trajectories, yet species parameters often remain poorly documented in simulation studies.We compiled and harmonized all existing tree species parameters available for the individual-based forest landscape and disturbance model (iLand). Since its first publication in 2012, iLand has been applied in 50 peer-reviewed publications across three continents throughout the Northern Hemisphere (i.e., Europe, North America, and Asia). The model operates at individual-tree level and simulates ecosystem processes at multiple spatial scales, making it a capable process-based model for studying forest change. However, the extensive number of processes and their interactions as well as the wide range of spatio-temporal scales considered in iLand require intensive parameterization, with tree species characterized by 66 unique parameters in the model. The database presented here includes parameters for 150 temperate and boreal tree species and provenances (i.e., regional variations). Excluding missing values, the database includes a total of 9,249 individual parameter entries. In addition, we provide parameters for the individual susceptibility of tree species to wind disturbance (five parameters) for a subset of 104 tree species and provenances (498 parameter entries). To guide further model parameterization efforts, we provide an estimate of uncertainty for each species based on how thoroughly simulations with the respective parameters were evaluated against independent data.Our dataset aids the future parameterization and application of iLand, and sets a new standard in documenting parameters used in process-based forest simulations. This dataset will support model application in previously unstudied areas and can facilitate the investigation of new tree species being introduced to well-studied systems (e.g., simulating assisted migration in the context of rapid climate change). Given that many process-based models rely on similar underlying processes our harmonized parameter set will be of relevance beyond the iLand community. Our work could catalyze further research into improving the parameterization of process-based forest models, increasing the robustness of projections of climate change impacts and adaptation strategies.

Anticipating consequences of disturbance interactions on ecosystem structure and function is a critical management priority as disturbance activity increases with warming climate. Across the Northern Hemisphere, extensive tree mortality from recent bark beetle outbreaks raises concerns about potential fire behavior and post-fire forest function. Silvicultural treatments (that is, partial or complete cutting of forest stands) may reduce outbreak severity and subsequent fuel loads, but longevity of pre-outbreak treatment effects on outbreak severity and post-outbreak fuel profiles remains underexplored. Further, treatments may present tradeoffs for other management objectives focused on ecosystem services (for example, carbon storage). We measured structure in old-growth subalpine forests following a recent (early 2000s) severe mountain pine beetle (MPB; Dendroctonus ponderosae ) outbreak to examine effects of historical (1940s) cutting intensity on gray stage (~10 years after peak of outbreak) post-outbreak (1) fuel profiles and (2) aboveground biomass carbon. Compared to control (uncut) stands, historically cut stands subjected to the same MPB outbreak had approximately half the post-outbreak surface fuel loads, about 2–3 x greater live canopy fuel loads, and greater within-stand spatial heterogeneity of dead canopy cover and available canopy fuel load. Post-outbreak total aboveground biomass carbon was similar across all stands, though historically cut stands had about 2 x greater carbon in live biomass compared to uncut stands. These findings suggest tradeoffs with altered post-outbreak potential fire behavior and carbon storage in cut stands. Additional implications of historical silvicultural treatments for wildlife habitat, firefighting operations, and long-term carbon trajectories highlight temporal legacies of management on directing forest response to interacting disturbances.

Growth releases of individuals that survive disturbances are important compensatory response mechanisms that contribute to ecological resilience. However, the role of fine-scale spatial heterogeneity in shaping compensatory growth responses is poorly understood for many broad-scale disturbances. We quantified how fine-scale spatial structure affects individual and aggregate tree growth leading up to and following a severe mountain pine beetle (MPB; Dendroctonus ponderosae) outbreak. We asked: (1) How does individual tree growth vary with tree- and neighborhood-scale characteristics? (2) How do within-stand aggregate growth and overstory recruitment vary with neighborhood-scale characteristics? We used a spatially explicit long-term monitoring dataset of a subalpine lodgepole pine (Pinus contorta var. latifolia) forest (in Colorado, USA) in which every tree ≥ 5 cm diameter was measured and mapped prior to (1989, 2004) and following (2018) a severe MPB outbreak (2003–2011). We used spatial regression to characterize drivers of growth. Overall, we found strong evidence for post-outbreak compensatory responses across spatial scales. Neighborhood characteristics shaped both individual and aggregate growth, with the magnitude of growth strongly mediated by pre-outbreak neighborhood structure and neighborhood mortality. Variation in tree-scale growth, combined with the spatial arrangement of surviving trees, resulted in highly variable emergent patterns of aggregate growth and recruitment. Our findings highlight the importance of fine-scale landscape configuration in shaping forest resilience. Quantifying compensatory responses in a spatially explicit framework at different scales is critical for modeling post-disturbance forest dynamics, which is increasingly important as climate warms and forest disturbance regimes change.

Tree regeneration is a critical mechanism of forest resilience to stand-replacing wildfire (i.e., where fire results in >90 % tree mortality), and post-fire regeneration is a concern worldwide as the climate becomes warmer. Although post-fire tree regeneration has been relatively well-studied in fire-prone forests across western North America, it is less understood in fire regimes characterized by large patches of stand-replacing fire at long intervals, such as the nominally infrequent, high-severity fire regimes of the western Cascades of Washington and northern Oregon, USA (northwestern Cascadia) where some of world's highest-biomass forests reside. Recent wildfire activity (2015-2020) in northwestern Cascadia provides a unique opportunity to build understanding of the mechanisms underpinning forest resilience to fire, with implications for forest management and the myriad socio-ecological values associated with these landscapes. Here, we asked: how is post-fire conifer regeneration across northwestern Cascadia affected by pre-fire stand age, burn-patch size, and topo-climatic conditions? We established 39 1-ha long-term monitoring plots across four recent wildfires in northwestern Cascadia to measure post-fire conifer abundance and species composition. We used generalized linear models to test how rates of conifer regeneration (seedlings ha(-1) yr(-1) of all heights) and establishment (seedlings ha(-1) yr(-1) with heights >10 cm) respond to differences in pre-fire stand age/seral stage (young, mid-seral, late-seral), distance to the nearest live seed source, and post-fire topo-climate setting. In general, initial post-fire seedling densities suggest tree regeneration is abundant following stand-replacing fire in most areas across the region. For example, at 3-5 years post-fire, established seedlings were present in 100 % of plots, and regeneration and establishment levels exceeded state forest practices minimum density thresholds (470 seedlings ha(-1)) in 82.1 % and 64.1 % of plots, respectively. Post-fire conifer regeneration and tree species richness were greater when pre-fire stand age was older, and was more dominated by shade tolerant species than when pre-fire stand age was younger. Regeneration decreased with farther distances to the nearest live seed source, though seedlings were still moderately abundant at distances as far as 400 m from edges of high-severity patches. Conifer establishment was greater in areas characterized by cooler and wetter macrosite (e.g., greater post-fire precipitation) and microsite (e.g., greater bryophyte ground cover) conditions. By providing information on the drivers and conditions in which land managers can expect natural regeneration to occur or not, our findings are directly relevant to pre- and post-fire forest management objectives under future climate and increased fire activity in systems structured in part by infrequent and severe wildfires.

Amplified by warming temperatures and drought, recent outbreaks of native bark beetles (Curculionidae: Scolytinae) have caused extensive tree mortality throughout Europe and North America. Despite their ubiquitous nature and important effects on ecosystems, forest recovery following such disturbances is poorly understood, particularly across regions with varying abiotic conditions and outbreak effects. To better understand post‐outbreak recovery across a topographically complex region, we synthesized data from 16 field studies spanning subalpine forests in the Southern Rocky Mountains, USA. From 1997 to 2019, these forests were heavily affected by outbreaks of three native bark beetle species ( Dendroctonus ponderosae , Dendroctonus rufipennis and Dryocoetes confusus ). We compared pre‐ and post‐outbreak forest conditions and developed region‐wide predictive maps of post‐outbreak (1) live basal areas, (2) juvenile densities and (3) height growth rates for the most abundant tree species – aspen ( Populus tremuloides ), Engelmann spruce ( Picea engelmannii ), lodgepole pine ( Pinus contorta ) and subalpine fir ( Abies lasiocarpa ). Beetle‐caused tree mortality reduced the average diameter of live trees by 28.4% (5.6 cm), and species dominance was altered on 27.8% of field plots with shifts away from pine and spruce. However, most plots (82.1%) were likely to recover towards pre‐outbreak tree densities without additional regeneration. Region‐wide maps indicated that fir and aspen, non‐host species for bark beetle species with the most severe effects (i.e. Dendroctonus spp.), will benefit from outbreaks through increased compositional dominance. After accounting for individual size, height growth for all conifer species was more rapid in sites with low winter precipitation, high winter temperatures and severe outbreaks. Synthesis . In subalpine forests of the US Rocky Mountains, recent bark beetle outbreaks have reduced tree size and altered species composition. While eventual recovery of the pre‐outbreak forest structure is likely in most places, changes in species composition may persist for decades. Still, forest communities following bark beetle outbreaks are widely variable due to differences in pre‐outbreak conditions, outbreak severity and abiotic gradients. This regional variability has critical implications for ecosystem services and susceptibility to future disturbances.

Promoting ecological resilience to increasing disturbance activity is a key management priority under warming climate. Across the Northern Hemisphere, tree mortality from widespread bark beetle outbreaks raises concerns for how forest management can foster resilience to future outbreaks. Density reduction (i.e., thinning) treatments can increase vigor of remaining trees, but the longevity of treatment efficacy for reducing susceptibility to future disturbance remains a key knowledge gap. Using one of the longest-running replicated experiments in old-growth subalpine forests, we measured stand structure following a recent (early 2000s) severe mountain pine beetle (MPB; Dendroctonus ponderosae) outbreak to examine the legacy of historical (1940s) thinning treatments on two components of resilience. We asked: 'How did historical thinning intensity affect (1) tree-scale survival probability and stand-scale survival proportion (collectively "resistance" to outbreak) for susceptible trees (lodgepole pine [Pinus contorta] >= 12 cm diameter) and (2) post-outbreak stand successional trajectories?' Overall outbreak severity was high (MPB killed 59% of susceptible individuals and 78% of susceptible basal area), and historical thinning had little effect on tree-scale and stand-scale resistance. Tree-scale survival probability decreased sharply with increasing tree diameter and did not differ from the control (uncut stands) in the historical thinning treatments. Stand-scale proportion of surviving susceptible trees and basal area did not differ from the control in historically thinned stands, except for treatments that removed nearly all susceptible trees, in which survival proportion approximately doubled. Despite limited effects on resistance to MPB outbreak, the legacy of historical treatments shifted dominance from large-diameter to small-diameter lodgepole pine by the time of outbreak, resulting in historically thinned stands with similar to 2x greater post-outbreak live basal area than control stands. MPB-driven mortality of large-diameter lodgepole pine in control stands and density-dependent mortality of small-diameter trees in historically thinned stands led to convergence in post-outbreak live tree stand structure. One exception was the heaviest historical thinning treatments (59-77% basal area removed), for which sapling dominance of shade-tolerant, unsusceptible conifers was lower than control stands. After six decades, thinning treatments have had minimal effect on resistance to bark beetle outbreaks, but leave persistent legacies in shaping post-outbreak successional trajectories.

Abstract Understanding drivers of disturbances across scales is critical as environmental constraints change in a warming climate. Outbreaks of native bark beetles (Curculionidae: Scolytinae) are key natural disturbances that shape the structure and function of conifer forests across the northern hemisphere. While drivers of bark beetle outbreaks have been studied extensively at spatial scales ranging from stands to continents, within‐stand processes governing individual tree mortality in an outbreak are less well understood. Here, we use a spatially explicit long‐term monitoring dataset of a lodgepole pine (Pinus contorta var. latifolia) forest (>9000 individually mapped trees in three 2‐ha plots) impacted by a severe mountain pine beetle (Dendroctonus ponderosae) outbreak to explore interactions among fine scale drivers of beetle‐caused tree mortality. Using a Bayesian spatial modeling approach, we evaluated how tree scale and tree neighborhood scale characteristics interact with tree size to mediate host tree susceptibility to mountain pine beetle outbreak in the Southern Rocky Mountains (USA). We found evidence that both tree growth rate preceding the outbreak and neighborhood structure (within a 10 meter radius of the host tree) mediate the effect of tree size, and that the direction and magnitude of these mediating effects vary with tree size. Tree scale mortality probability increased with pre‐outbreak growth rate for small to medium sized host trees (~10–25 cm diameter), but that same effect was not detected for large trees. Conversely, tree scale mortality probability increased with greater neighborhood density, with the most pronounced effects for medium to large sized host trees (~15–30 cm diameter). Within‐stand topographic variability was not an important predictor of mortality probability; among stands, however, the stand in the driest topographic position experienced the greatest overall mortality. By explicitly considering how within‐stand heterogeneity mediates individual tree scale susceptibility to bark beetle outbreak, our findings bridge an important gap in understanding multi‐scale drivers of disturbance dynamics.