Anthropogenic land-use legacies underpin climate change-related risks to forest ecosystems

Forest ecosystems with long-lasting human imprints can emerge worldwide as outcomes of land-use cessation. However, the interaction of these anthropogenic legacies with climate change impacts on forests is not well understood. Here, we set out how anthropogenic land-use legacies that persist in forest properties, following alterations in forest distribution, structure, and composition, can interact with climate change stressors. We propose a risk-based framework to identify anthropogenic legacies of land uses in forest ecosystems and quantify the impact of their interaction with climate-related stress on forest responses. Considering anthropogenic land-use legacies alongside environmental drivers of forest ecosystem dynamics will improve our predictive capacity of climate-related risks to forests and our ability to promote ecosystem resilience to climate change.

Quantifying carbon storage and sequestration by native and non-native forests under contrasting climate types

Non-native trees may have significant impacts on the carbon sink capacity of forested lands. However, large-scale patterns of the relative capacity of native and non-native forests to uptake and store carbon remain poorly described in the literature, and this information is urgently needed to support management decisions. In this study, we analyzed 17,065 plots from the Spanish Forest Inventory (covering c. 30 years) to quantify carbon storage and sequestration of natural forests and plantations of native and non-native trees under contrasting climate types, while controlling for the effects of environmental factors (forest structure, climate, soil, topography, and management). We found that forest origin (non-native vs. native) highly influenced carbon storage and sequestration, but such effect was dependent on climate. Carbon storage was greater in non-native than in native forests in both wet and dry climates. Non-native forests also had greater carbon sequestration than native ones in the wet climate, due to higher carbon gains by tree growth. However, in the dry climate, native forests had greater carbon gains by tree ingrowth and lower carbon loss by tree mortality than non-native ones. Furthermore, forest type (classified by the dominant species) and natural forests versus tree plantations were important determinants of carbon storage and sequestration. Native and non-native Pinus spp. forests had low carbon storage, whereas non-native Eucalyptus spp. forests and native Quercus spp., Fagus sylvatica, and Eurosiberian mixed forests (especially not planted ones) had high carbon storage. Carbon sequestration was greatest in Eucalyptus globulus, Quercus ilex, and Pinus pinaster forests. Overall, our findings suggest that the relative capacity of native and non-native forests to uptake and store carbon depends on climate, and that the superiority of non-native forests over native ones in terms of carbon sequestration declines as the abiotic filters become stronger (i.e., lower water availability and higher climate seasonality).

Functional traits and climate drive interspecific differences in disturbance-induced tree mortality

With climate change, natural disturbances such as storm or fire are reshuffled, inducing pervasive shifts in forest dynamics. To predict how it will impact forest structure and composition, it is crucial to understand how tree species differ in their sensitivity to disturbances. In this study, we investigated how functional traits and species mean climate affect their sensitivity to disturbances while controlling for tree size and stand structure. With data on 130,594 trees located on 7617 plots that were disturbed by storm, fire, snow, biotic or other disturbances from the French, Spanish, and Finnish National Forest Inventory, we modeled annual mortality probability for 40 European tree species as a function of tree size, dominance status, disturbance type, and intensity. We tested the correlation of our estimated species probability of disturbance mortality with their traits and their mean climate niches. We found that different trait combinations controlled species sensitivity to disturbances. Storm-sensitive species had a high height-dbh ratio, low wood density and high maximum growth, while fire-sensitive species had low bark thickness and high P50. Species from warmer and drier climates, where fires are more frequent, were more resistant to fire. The ranking in disturbance sensitivity between species was overall consistent across disturbance types. Productive conifer species were the most disturbance sensitive, while Mediterranean oaks were the least disturbance sensitive. Our study identified key relations between species functional traits and disturbance sensitivity, that allows more reliable predictions of how changing climate and disturbance regimes will impact future forest structure and species composition at large spatial scales.

Tree growth response to drought partially explains regional-scale growth and mortality patterns in Iberian forests

Tree-ring data has been widely used to inform about tree growth responses to drought at the individual scale, but less is known about how tree growth sensitivity to drought scales up driving changes in forest dynamics. Here, we related tree-ring growth chronologies and stand-level forest changes in basal area from two independent data sets to test if tree-ring responses to drought match stand forest dynamics (stand basal area growth, ingrowth, and mortality). We assessed if tree growth and changes in forest basal area covary as a function of spatial scale and tree taxa (gymnosperm or angiosperm). To this end, we compared a tree-ring network with stand data from the Spanish National Forest Inventory. We focused on the cumulative impact of drought on tree growth and demography in the period 1981-2005. Drought years were identified by the Standardized Precipitation Evapotranspiration Index, and their impacts on tree growth by quantifying tree-ring width reductions. We hypothesized that forests with greater drought impacts on tree growth will also show reduced stand basal area growth and ingrowth and enhanced mortality. This is expected to occur in forests dominated by gymnosperms on drought-prone regions. Cumulative growth reductions during dry years were higher in forests dominated by gymnosperms and presented a greater magnitude and spatial autocorrelation than for angiosperms. Cumulative drought-induced tree growth reductions and changes in forest basal area were related, but initial stand density and basal area were the main factors driving changes in basal area. In drought-prone gymnosperm forests, we observed that sites with greater growth reductions had lower stand basal area growth and greater mortality. Consequently, stand basal area, forest growth, and ingrowth in regions with large drought impacts was significantly lower than in regions less impacted by drought. Tree growth sensitivity to drought can be used as a predictor of gymnosperm demographic rates in terms of stand basal area growth and ingrowth at regional scales, but further studies may try to disentangle how initial stand density modulates such relationships. Drought-induced growth reductions and their cumulative impacts have strong potential to be used as early-warning indicators of regional forest vulnerability.

Linking functional traits with tree growth and forest productivity in Quercus ilex forests along a climatic gradient

Plant functional traits are highly plastic to changes in climatic factors and nutrient availability. However, the intraspecific plant response to abiotic factors and the overall effect on tree growth and productivity is still under debate. We studied forest productivity for 30 Quercus ilex subsp. ballota forests in Spain along a broad climatic gradient of aridity (mean annual precipitation from 321 to 858 mm). We used linear mixed models to quantify the effect of climatic and edaphic (soil nutrients, topography, and texture) factors on tree functional traits (leaf and branch traits), and subsequently, the effect of such functional traits and abiotic factors on the relative growth rate (RGR) of adult trees. We used piecewise structural equation models (SEMs) to determine the causal effect of intrinsic and extrinsic factors on forest productivity. Our results showed that tree functional traits were mainly explained by climatic and edaphic factors. Functional traits and tree biomass explained forest biomass and RGR, respectively, which ultimately explained forest productivity. In conclusion, intraspecific variability of functional traits has a significant effect on plant biomass and growth, which ultimately may explain forest productivity in Quercus ilex forests.

Climate reverses directionality in the richness-abundance relationship across the World's main forest biomes

More tree species can increase the carbon storage capacity of forests (here referred to as the more species hypothesis) through increased tree productivity and tree abundance resulting from complementarity, but they can also be the consequence of increased tree abundance through increased available energy (more individuals hypothesis). To test these two contrasting hypotheses, we analyse the most plausible pathways in the richness-abundance relationship and its stability along global climatic gradients. We show that positive effect of species richness on tree abundance only prevails in eight of the twenty-three forest regions considered in this study. In the other forest regions, any benefit from having more species is just as likely (9 regions) or even less likely (6 regions) than the effects of having more individuals. We demonstrate that diversity effects prevail in the most productive environments, and abundance effects become dominant towards the most limiting conditions. These findings can contribute to refining cost-effective mitigation strategies based on fostering carbon storage through increased tree diversity. Specifically, in less productive environments, mitigation measures should promote abundance of locally adapted and stress tolerant tree species instead of increasing species richness.

Evidence of non-stationary relationships between climate and forest responses: Increased sensitivity to climate change in Iberian forests

Climate and forest structure are considered major drivers of forest demography and productivity. However, recent evidence suggests that the relationships between climate and tree growth are generally non-stationary (i.e. non-time stable), and it remains uncertain whether the relationships between climate, forest structure, demography and productivity are stationary or are being altered by recent climatic and structural changes. Here we analysed three surveys from the Spanish Forest Inventory covering c. 30 years of information and we applied mixed and structural equation models to assess temporal trends in forest structure (stand density, basal area, tree size and tree size inequality), forest demography (ingrowth, growth and mortality) and above-ground forest productivity. We also quantified whether the interactive effects of climate and forest structure on forest demography and above-ground forest productivity were stationary over two consecutive time periods. Since the 1980s, density, basal area and tree size increased in Iberian forests, and tree size inequality decreased. In addition, we observed reductions in ingrowth and growth, and increases in mortality. Initial forest structure and water availability mainly modulated the temporal trends in forest structure and demography. The magnitude and direction of the interactive effects of climate and forest structure on forest demography changed over the two time periods analysed indicating non-stationary relationships between climate, forest structure and demography. Above-ground forest productivity increased due to a positive balance between ingrowth, growth and mortality. Despite increasing productivity over time, we observed an aggravation of the negative effects of climate change and increased competition on forest demography, reducing ingrowth and growth, and increasing mortality. Interestingly, our results suggest that the negative effects of climate change on forest demography could be ameliorated through forest management, which has profound implications for forest adaptation to climate change.

Multiple abiotic and biotic pathways shape biomass demographic processes in temperate forests

Forests play a key role in regulating the global carbon cycle, and yet the abiotic and biotic conditions that drive the demographic processes that underpin forest carbon dynamics remain poorly understood in natural ecosystems. To address this knowledge gap, we used repeat forest inventory data from 92,285 trees across four large permanent plots (4–25 ha in size) in temperate mixed forests in northeast China to ask the following questions: (1) How do soil conditions and stand age drive biomass demographic processes? (2) How do vegetation quality (i.e., functional trait diversity and composition) and quantity (i.e., initial biomass stocks) influence biomass demographic processes independently from soil conditions and stand age? (3) What is the relative contribution of growth, recruitment, and mortality to net biomass change? Using structural equation modeling, we showed that all three demographic processes were jointly constrained by multiple abiotic and biotic factors and that mortality was the strongest determinant on net biomass change over time. Growth and mortality, as well as functional trait diversity and the community‐weighted mean of specific leaf area (CWM_SLA), declined with stand age. By contrast, high soil phosphorous concentrations were associated with greater functional diversity and faster dynamics (i.e., high growth and mortality rates), but associated with lower CWM_SLA and initial biomass stock. More functionally diverse communities also had higher recruitment rates, but did not exhibit faster growth and mortality. Instead, initial biomass stocks and CWM_SLA were stronger predictors of biomass growth and mortality, respectively. By integrating the full spectrum of abiotic and biotic drivers of forest biomass dynamics, our study provides critical system‐level insights needed to predict the possible consequences of regional changes in forest diversity, composition, structure and function in the context of global change.

Evolutionary history predicts the response of tree species to forest loss: A case study in peninsular Spain

Evolutionary history can explain species resemblance to a large extent. Thus, if closely related species share combinations of traits that modulate their response to environmental changes, then phylogeny could predict species sensitivity to novel stressors such as increased levels of deforestation. To test this hypothesis, we used 66,949 plots (25-m-radius) of the Spanish National Forest Inventory and modelled the relationships between local (plot-level) stem density of 61 Holarctic tree species and forest canopy cover measured at local and landscape scales (concentric circles centred on the plots with radiuses of 1.6, 3.2 and 6.4 km, respectively). Then, we used the output model equations to estimate the probability of occurrence of the species as a function of forest canopy cover (i.e. response to forest loss), and quantified the phylogenetic signal in their responses using a molecular phylogeny. Most species showed a lower probability of occurrence when forest canopy cover in the plots (local scale) was low. However, the probability of occurrence of many species increased when forest canopy cover decreased across landscape scales. We detected a strong phylogenetic signal in species response to forest loss at local and small landscape (1.6 km) scales. However, phylogenetic signal was weak and non-significant at intermediate (3.2 km) and large (6.4 km) landscape scales. Our results suggest that phylogenetic information could be used to prioritize forested areas for conservation, since evolutionary history may largely determine species response to forest loss. As such, phylogenetically diverse forests might ensure contrasted responses to deforestation, and thus less abrupt reductions in the abundances of the constituent species.

Continental mapping of forest ecosystem functions reveals a high but unrealised potential for forest multifunctionality

Humans require multiple services from ecosystems, but it is largely unknown whether trade-offs between ecosystem functions prevent the realisation of high ecosystem multifunctionality across spatial scales. Here, we combined a comprehensive dataset (28 ecosystem functions measured on 209 forest plots) with a forest inventory dataset (105,316 plots) to extrapolate and map relationships between various ecosystem multifunctionality measures across Europe. These multifunctionality measures reflected different management objectives, related to timber production, climate regulation and biodiversity conservation/recreation. We found that trade-offs among them were rare across Europe, at both local and continental scales. This suggests a high potential for 'win-win' forest management strategies, where overall multifunctionality is maximised. However, across sites, multifunctionality was on average 45.8-49.8% below maximum levels and not necessarily highest in protected areas. Therefore, using one of the most comprehensive assessments so far, our study suggests a high but largely unrealised potential for management to promote multifunctional forests.

Climate- and successional-related changes in functional composition of European forests are strongly driven by tree mortality

Intense droughts combined with increased temperatures are one of the major threats to forest persistence in the 21st century. Despite the direct impact of climate change on forest growth and shifts in species abundance, the effect of altered demography on changes in the composition of functional traits is not well known. We sought to (1) quantify the recent changes in functional composition of European forests; (2) identify the relative importance of climate change, mean climate and forest development for changes in functional composition; and (3) analyse the roles of tree mortality and growth underlying any functional changes in different forest types. We quantified changes in functional composition from the 1980s to the 2000s across Europe by two dimensions of functional trait variation: the first dimension was mainly related to changes in leaf mass per area and wood density (partially related to the trait differences between angiosperms and gymnosperms), and the second dimension was related to changes in maximum tree height. Our results indicate that climate change and mean climatic effects strongly interacted with forest development and it was not possible to completely disentangle their effects. Where recent climate change was not too extreme, the patterns of functional change generally followed the expected patterns under secondary succession (e.g. towards late-successional short-statured hardwoods in Mediterranean forests and taller gymnosperms in boreal forests) and latitudinal gradients (e.g. larger proportion of gymnosperm-like strategies at low water availability in forests formerly dominated by broad-leaved deciduous species). Recent climate change generally favoured the dominance of angiosperm-like related traits under increased temperature and intense droughts. Our results show functional composition changes over relatively short time scales in European forests. These changes are largely determined by tree mortality, which should be further investigated and modelled to adequately predict the impacts of climate change on forest function.

Structural overshoot of tree growth with climate variability and the global spectrum of drought-induced forest dieback

Ongoing climate change poses significant threats to plant function and distribution. Increased temperatures and altered precipitation regimes amplify drought frequency and intensity, elevating plant stress and mortality. Large-scale forest mortality events will have far-reaching impacts on carbon and hydrological cycling, biodiversity, and ecosystem services. However, biogeographical theory and global vegetation models poorly represent recent forest die-off patterns. Furthermore, as trees are sessile and long-lived, their responses to climate extremes are substantially dependent on historical factors. We show that periods of favourable climatic and management conditions that facilitate abundant tree growth can lead to structural overshoot of aboveground tree biomass due to a subsequent temporal mismatch between water demand and availability. When environmental favourability declines, increases in water and temperature stress that are protracted, rapid, or both, drive a gradient of tree structural responses that can modify forest self-thinning relationships. Responses ranging from premature leaf senescence and partial canopy dieback to whole-tree mortality reduce canopy leaf area during the stress period and for a lagged recovery window thereafter. Such temporal mismatches of water requirements from availability can occur at local to regional scales throughout a species geographical range. As climate change projections predict large future fluctuations in both wet and dry conditions, we expect forests to become increasingly structurally mismatched to water availability and thus overbuilt during more stressful episodes. By accounting for the historical context of biomass development, our approach can explain previously problematic aspects of large-scale forest mortality, such as why it can occur throughout the range of a species and yet still be locally highly variable, and why some events seem readily attributable to an ongoing drought while others do not. This refined understanding can facilitate better projections of structural overshoot responses, enabling improved prediction of changes in forest distribution and function from regional to global scales.

Tree mortality across biomes is promoted by drought intensity, lower wood density and higher specific leaf area

Drought events are increasing globally, and reports of consequent forest mortality are widespread. However, due to a lack of a quantitative global synthesis, it is still not clear whether drought-induced mortality rates differ among global biomes and whether functional traits influence the risk of drought-induced mortality. To address these uncertainties, we performed a global meta-analysis of 58 studies of drought-induced forest mortality. Mortality rates were modelled as a function of drought, temperature, biomes, phylogenetic and functional groups and functional traits. We identified a consistent global-scale response, where mortality increased with drought severity [log mortality (trees trees^-1  year^-1 ) increased 0.46 (95% CI = 0.2-0.7) with one SPEI unit drought intensity]. We found no significant differences in the magnitude of the response depending on forest biomes or between angiosperms and gymnosperms or evergreen and deciduous tree species. Functional traits explained some of the variation in drought responses between species (i.e. increased from 30 to 37% when wood density and specific leaf area were included). Tree species with denser wood and lower specific leaf area showed lower mortality responses. Our results illustrate the value of functional traits for understanding patterns of drought-induced tree mortality and suggest that mortality could become increasingly widespread in the future.

Complementarity effects on tree growth are contingent on tree size and climatic conditions across Europe

Neglecting tree size and stand structure dynamics might bias the interpretation of the diversity-productivity relationship in forests. Here we show evidence that complementarity is contingent on tree size across large-scale climatic gradients in Europe. We compiled growth data of the 14 most dominant tree species in 32,628 permanent plots covering boreal, temperate and Mediterranean forest biomes. Niche complementarity is expected to result in significant growth increments of trees surrounded by a larger proportion of functionally dissimilar neighbours. Functional dissimilarity at the tree level was assessed using four functional types: i.e. broad-leaved deciduous, broad-leaved evergreen, needle-leaved deciduous and needle-leaved evergreen. Using Linear Mixed Models we show that, complementarity effects depend on tree size along an energy availability gradient across Europe. Specifically: (i) complementarity effects at low and intermediate positions of the gradient (coldest-temperate areas) were stronger for small than for large trees; (ii) in contrast, at the upper end of the gradient (warmer regions), complementarity is more widespread in larger than smaller trees, which in turn showed negative growth responses to increased functional dissimilarity. Our findings suggest that the outcome of species mixing on stand productivity might critically depend on individual size distribution structure along gradients of environmental variation.

Revealing patterns of local species richness along environmental gradients with a novel network tool

How species richness relates to environmental gradients at large extents is commonly investigated aggregating local site data to coarser grains. However, such relationships often change with the grain of analysis, potentially hiding the local signal. Here we show that a novel network technique, the “method of reflections”, could unveil the relationships between species richness and climate without such drawbacks. We introduced a new index related to potential species richness, which revealed large scale patterns by including at the local community level information about species distribution throughout the dataset (i.e., the network). The method effectively removed noise, identifying how far site richness was from potential. When applying it to study woody species richness patterns in Spain, we observed that annual precipitation and mean annual temperature explained large parts of the variance of the newly defined species richness, highlighting that, at the local scale, communities in drier and warmer areas were potentially the species richest. Our method went far beyond what geographical upscaling of the data could unfold, and the insights obtained strongly suggested that it is a powerful instrument to detect key factors underlying species richness patterns, and that it could have numerous applications in ecology and other fields.

Climatic Stress during Stand Development Alters the Sign and Magnitude of Age-Related Growth Responses in a Subtropical Mountain Pine

The modification of typical age-related growth by environmental changes is poorly understood, In part because there is a lack of consensus at individual tree level regarding age-dependent growth responses to climate warming as stands develop. To increase our current understanding about how multiple drivers of environmental change can modify growth responses as trees age we used tree ring data of a mountain subtropical pine species along an altitudinal gradient covering more than 2,200 m of altitude. We applied mixed-linear models to determine how absolute and relative age-dependent growth varies depending on stand development; and to quantify the relative importance of tree age and climate on individual tree growth responses. Tree age was the most important factor for tree growth in models parameterised using data from all forest developmental stages. Contrastingly, the relationship found between tree age and growth became non-significant in models parameterised using data corresponding to mature stages. These results suggest that although absolute tree growth can continuously increase along tree size when trees reach maturity age had no effect on growth. Tree growth was strongly reduced under increased annual temperature, leading to more constant age-related growth responses. Furthermore, young trees were the most sensitive to reductions in relative growth rates, but absolute growth was strongly reduced under increased temperature in old trees. Our results help to reconcile previous contrasting findings of age-related growth responses at the individual tree level, suggesting that the sign and magnitude of age-related growth responses vary with stand development. The different responses found to climate for absolute and relative growth rates suggest that young trees are particularly vulnerable under warming climate, but reduced absolute growth in old trees could alter the species' potential as a carbon sink in the future.

Patterns and drivers of tree mortality in iberian forests: climatic effects are modified by competition

Tree mortality is a key process underlying forest dynamics and community assembly. Understanding how tree mortality is driven by simultaneous drivers is needed to evaluate potential effects of climate change on forest composition. Using repeat-measure information from c. 400,000 trees from the Spanish Forest Inventory, we quantified the relative importance of tree size, competition, climate and edaphic conditions on tree mortality of 11 species, and explored the combined effect of climate and competition. Tree mortality was affected by all of these multiple drivers, especially tree size and asymmetric competition, and strong interactions between climate and competition were found. All species showed L-shaped mortality patterns (i.e. showed decreasing mortality with tree size), but pines were more sensitive to asymmetric competition than broadleaved species. Among climatic variables, the negative effect of temperature on tree mortality was much larger than the effect of precipitation. Moreover, the effect of climate (mean annual temperature and annual precipitation) on tree mortality was aggravated at high competition levels for all species, but especially for broadleaved species. The significant interaction between climate and competition on tree mortality indicated that global change in Mediterranean regions, causing hotter and drier conditions and denser stands, could lead to profound effects on forest structure and composition. Therefore, to evaluate the potential effects of climatic change on tree mortality, forest structure must be considered, since two systems of similar composition but different structure could radically differ in their response to climatic conditions.