How to kill a tree: empirical mortality models for 18 species and their performance in a dynamic forest model

Dynamic Vegetation Models (DVMs) are designed to be suitable for simulating forest succession and species range dynamics under current and future conditions based on mathematical representations of the three key processes regeneration, growth, and mortality. However, mortality formulations in DVMs are typically coarse and often lack an empirical basis, which increases the uncertainty of projections of future forest dynamics and hinders their use for developing adaptation strategies to climate change. Thus, sound tree mortality models are highly needed. We developed parsimonious, species-specific mortality models for 18 European tree species using >90,000 records from inventories in Swiss and German strict forest reserves along a considerable environmental gradient. We comprehensively evaluated model performance and incorporated the new mortality functions in the dynamic forest model ForClim. Tree mortality was successfully predicted by tree size and growth. Only a few species required additional covariates in their final model to consider aspects of stand structure or climate. The relationships between mortality and its predictors reflect the indirect influences of resource availability and tree vitality, which are further shaped by species-specific attributes such as maximum longevity and shade tolerance. Considering that the behavior of the models was biologically meaningful, and that their performance was reasonably high and not impacted by changes in the sampling design, we suggest that the mortality algorithms developed here are suitable for implementation and evaluation in DVMs. In the DVM ForClim, the new mortality functions resulted in simulations of stand basal area and species composition that were generally close to historical observations. However, ForClim performance was poorer than when using the original, coarse mortality formulation. The difficulties of simulating stand structure and species composition, which were most evident for Fagus sylvatica L. and in long-term simulations, resulted from feedbacks between simulated growth and mortality as well as from extrapolation to very small and very large trees. Growth and mortality processes and their species-specific differences should thus be revisited jointly, with a particular focus on small and very large trees in relation to their shade tolerance.