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

Nonnative tree invaders lead to declines in native tree species richness

Biological invasions are profoundly altering Earth's ecosystems, but generalities about the effects of nonnative species on the diversity and productivity of native communities have been elusive. This lack of generality may reflect the limited spatial and temporal extents of most previous studies. Using >5 million tree measurements across eastern US forests from 1995 to 2023, we quantified temporal trends in tree diversity and biomass. We then analyzed community-level changes in native tree diversity and biomass in relation to nonnative tree invasion and native species colonization. Across the entire eastern United States, native tree species richness decreased over time in plots where nonnatives occurred, whereas nonnative species richness and the biomass of both natives and nonnatives increased over time. At the community scale, native richness tended to decline following nonnative invasion, whereas native biomass and richness-independent measures of trait and phylogenetic diversity tended to remain stable. These patterns can be explained by the rarity of the displaced native species and their functional and phylogenetic similarity to native species that survived nonnative invasions. In contrast, native survivors tended to be functionally distinct from nonnative invaders, suggesting an important role for niche partitioning in community dynamics. Colonization by previously absent native species was associated with an increase in native richness (beyond the addition of native colonizers), which contrasts with declines in native richness that tended to follow nonnative invasion. These results suggest a causal role for nonnative species in the native richness decline of invaded communities.

Enhancing economic multifunctionality without compromising multidiversity and ecosystem multifunctionality via forest enrichment

Enriching tree species-poor and less productive forests by introducing economically valuable species is a strategy proposed for achieving multipurpose forest management. However, empirical evidence from managed and mature forests on the impact of this enrichment on ecological (multidiversity and ecosystem multifunctionality) and economic dimensions remains scarce, particularly when nonnative species are used. Here, we propose and test a framework that integrates economic multifunctionality, encompassing timber production-oriented goals and resistance against disturbances, with multidiversity and ecosystem multifunctionality in European beech forest stands enriched with conifers. Our results show that enriched beech forest stands (~80 years old) can provide high levels of economic multifunctionality without compromising multidiversity and ecosystem multifunctionality. In comparison to pure beech stands, enriched stands with Douglas-fir supported win-win-win situations for these three dimensions. Our findings contribute to the discussion of integrating biodiversity, ecosystem, and economic functions, providing empirical evidence for future forest management.

Testing the ectomycorrhizal-dominance hypothesis for ecosystem multifunctionality in a subtropical mountain forest

Mycorrhizal associations are key mutualisms that shape the structure of forest communities and multiple ecosystem functions. However, we lack a framework for predicting the varying dominance of distinct mycorrhizal associations in an integrated proxy of multifunctionality across ecosystems. Here, we used the datasets containing diversity of mycorrhizal associations and 18 ecosystem processes related to supporting, provisioning, and regulating services to examine how the dominance of ectomycorrhiza (EcM) associations affects ecosystem multifunctionality in subtropical mountain forests in Southwest China. Meanwhile, we synthesized the prevalence of EcM-dominant effects on ecosystem functioning in forest biomes. Our results demonstrated that elevation significantly modified the distributions of EcM trees and fungal dominance, which in turn influenced multiple functions simultaneously. Multifunctionality increased with increasing proportion of EcM associations, supporting the ectomycorrhizal-dominance hypothesis. Meanwhile, we observed that the impacts of EcM dominance on individual ecosystem functions exhibited different relationships among forest biomes. Our findings highlight the importance of ectomycorrhizal dominance in regulating multifunctionality in subtropical forests. However, this ectomycorrhizal feedback in shaping ecosystem functions cannot necessarily be generalized across forests. Therefore, we argue that the predictions for ecosystem multifunctionality in response to the shifts of mycorrhizal composition could vary across space and time.

Ecosystem services at risk from disturbance in Europe's forests

Global change impacts on disturbances can strongly compromise the capacity of forests to provide ecosystem services to society. In addition, many ecosystem services in Europe are simultaneously provided by forests, emphasizing the importance of multifunctionality in forest ecosystem assessments. To address disturbances in forest ecosystem policies and management, spatially explicit risk analyses that consider multiple disturbances and ecosystem services are needed. However, we do not yet know which ecosystem services are most at risk from disturbances in Europe, where the respective risk hotspots are, nor which of the main disturbance agents are most detrimental to the provisioning of multiple ecosystem services from Europe's forests. Here, we quantify the risk of losing important ecosystem services (timber supply, carbon storage, soil erosion control and outdoor recreation) to forest disturbances (windthrows, bark beetle outbreaks and wildfires) in Europe on a continental scale. We find that up to 12% of Europe's ecosystem service supply is at risk from current disturbances. Soil erosion control is the ecosystem service at the highest risk, and windthrow is the disturbance agent posing the highest risk. Disturbances challenge forest multifunctionality by threatening multiple ecosystem services simultaneously on 19.8 Mha (9.7%) of Europe's forests. Our results highlight priority areas for risk management aiming to safeguard the sustainable provisioning of forest ecosystem services.

Trophic interactions in soil micro-food webs drive ecosystem multifunctionality along tree species richness

Rapid biodiversity losses under global climate change threaten forest ecosystem functions. However, our understanding of the patterns and drivers of multiple ecosystem functions across biodiversity gradients remains equivocal. To address this important knowledge gap, we measured simultaneous responses of multiple ecosystem functions (nutrient cycling, soil carbon stocks, organic matter decomposition, plant productivity) to a tree species richness gradient of 1, 4, 8, 16, and 32 species in a young subtropical forest. We found that tree species richness had negligible effects on nutrient cycling, organic matter decomposition, and plant productivity, but soil carbon stocks and ecosystem multifunctionality significantly increased with tree species richness. Linear mixed-effect models showed that soil organisms, particularly arbuscular mycorrhizal fungi (AMF) and soil nematodes, elicited the greatest relative effects on ecosystem multifunctionality. Structural equation models revealed indirect effects of tree species richness on ecosystem multifunctionality mediated by trophic interactions in soil micro-food webs. Specifically, we found a significant negative effect of gram-positive bacteria on soil nematode abundance (a top-down effect), and a significant positive effect of AMF biomass on soil nematode abundance (a bottom-up effect). Overall, our study emphasizes the significance of a multitrophic perspective in elucidating biodiversity-multifunctionality relationships and highlights the conservation of functioning soil micro-food webs to maintain multiple ecosystem functions.

Rhizosphere microbial community assembly and association networks strongly differ based on vegetation type at a local environment scale

Introduction

Rhizosphere microbes perform critical functions for their hosts, and their structure is strongly influenced by vegetation type. Although studies on the effects of vegetation on rhizosphere microbial community structure have been conducted at large and global environment scales, studies at local environment scales would eliminate numerous external factors such as climate and soil type, while highlighting the potential influence of local vegetation type.

Methods

Here, we compared rhizosphere microbial communities using 54 samples under three vegetation types (herb, shrubs, and arbors, with bulk soil as the control) at the campus of Henan University. 16S rRNA and ITS amplicons were sequenced using Illumina high throughput sequencing.

Results and Discussion

Rhizosphere bacterial and fungal community structures were influenced considerably by vegetation type. Bacterial alpha diversity under herbs was significantly different from that under arbors and shrubs. The abundance of phyla such as Actinobacteria was extremely higher in bulk soil than in the rhizosphere soils. Herb rhizosphere harbored more unique species than other vegetation type soils. Furthermore, bacterial community assembly in bulk soil was more dominated by deterministic process, whereas the rhizosphere bacterial community assembly was dominated by stochasticity and the construction of fungal communities was all dominated by deterministic processes. In addition, rhizosphere microbial networks were less complex than bulk soil networks, and their keystone species differed based on vegetation type. Notably, bacterial community dissimilarities were strongly correlated with plant phylogenetic distance. Exploring rhizosphere microbial community patterns under different vegetation types could enhance our understanding of the role of rhizosphere microbes in ecosystem function and service provision, as well as basic information that could facilitate plant and microbial diversity conservation at the local environment scale.

Size matters for linking traits to ecosystem multifunctionality

A priority research field addresses how to optimize diverse ecosystem services to people, including biodiversity support, regulatory, utilitarian and cultural services. This field may benefit from linking ecosystem services to the sizes of different body parts of organisms, with functional traits as the go-between. Using woody ecosystems to explore such linkages, we hypothesize that across stem diameter classes from trunk via branches to twigs, key wood and bark functional traits (especially those defining size-shape and resource economics spectra) vary both within individual trees and shrubs and across woody species, thereby together boosting ecosystem multifunctionality. While we focus on woody plants aboveground, we discuss promising extensions to belowground organs of trees and shrubs and analogs with other organisms, for example, vertebrate animals.

Long-Term Productivity of Monospecific and Mixed Oak (Quercus petraea [Matt.] Liebl. and Quercus robur L.) Stands in Germany: Growth Dynamics and the Effect of Stand Structure

Wood production is one of the most important ecosystem service that forests provide to society. However, under changing climatic conditions, this appears to be subject to increasing uncertainties. In the present study we analyzed how long-term productivity of oak (Quercus petraea [Matt.] Liebl. and Quercus robur L.) stands has developed, how oak behaved on tree and stand level depending on the stand structure and which trade-offs can be observed. For the analyses, data from 147 long-term monospecific and mixed stands were investigated, which have been regularly recorded since 1898. Firstly, long-term stand productivity has increased up to 21% until 2020 as compared to 1960. This trend was observed for both, monospecific as well as mixed oak stands. Secondly, stand productivity was on average 19% higher in mixed compared to monospecific oak stands. This superiority can be explained by higher stand densities, a vigorous understory and the admixture of beech in particular. With increasing age, the observed positive effect of stand density was higher. Thirdly, individual oak productivity slowed down under interspecific competition, especially in young to mid-aged stands. In this context, the productivity of individual oaks depended strongly on their social position within the stand. Fourthly, in terms of growth partitioning larger trees contributed most in young oak stands, regardless of mixture. In order to preserve oak as a productive component of future mixed forests, the results suggest a silvicultural promotion of oak. Consistent management of dominant and vital oaks can achieve high productive trees while maintaining the positive characteristics of highly structured and mixed forests. A vigorous secondary stand can increase overall stand productivity at lower densities and allows silvicultural flexibility at the stand level. Creating vertical stand structure to reduce competition has only a limited positive effect on productivity of individual oaks that is highly related to its social status. Special attention should still be paid to beech as admixed tree species, which can continue to crowd oak even at higher stand ages.

Temporal Soundscape Patterns in a Panamanian Tree Diversity Experiment: Polycultures Show an Increase in High Frequency Cover

In this ecoacoustic study we used the setting of a tropical tree diversity planted forest to analyze temporal patterns in the composition of soundscapes and to test the effects of tree species richness on associated biodiversity measured as acoustic diversity. The analysis of soundscapes offers easy, rapid and sustainable methods when assessing biodiversity. During the last years the quantification of regional or global acoustic variability in sounds and the analysis of different soundscapes has been evolving into an important tool for biodiversity conservation, especially since case studies confirmed a relationship between land-use management, forest structure and acoustic diversity. Here we analyzed soundscapes from two seasons (dry and rainy season) and aurally inspected a subset of audio recordings to describe temporal patterns in soundscape composition. Several acoustic indices were calculated and we performed a correlation analysis and a non-metric multidimensional scaling analysis to identify acoustic indices that: (i) were complementary to each other and such represented different aspects of the local soundscapes and (ii) related most strongly to differences in acoustic composition among tree species richness, season and day phase. Thus, we chose “High Frequency Cover,” “Bioacoustic Index,” and “Events Per Second” to test the hypothesis that acoustic diversity increases with increasing tree species richness. Monocultures differed significantly from polycultures during night recordings, with respect to High Frequency Cover. This index covers sounds above 8 kHz and thus represents part of the orthopteran community. We conclude that increasing tree species richness in a young tropical forest plantation had positive effects on the vocalizing communities. The strongest effects were found for acoustic activity of the orthopteran community. In contrast to birds, orthopterans have smaller home ranges, and are therefore important indicator species for small scale environmental conditions.

Taxonomic, structural diversity and carbon stocks in a gradient of island forests

Assessment of forest ecosystems and their services is seen as a key action for the advancement of biodiversity objectives, and to inform the development and implementation of related policies and planning. The Azorean forest is important for timber production, the protection of soil and water resources, and for its recreational and aesthetic value. However, its role in carbon accumulation has not been fully addressed. We assessed plant diversity, forest structure and carbon stocks in a gradient of three forest types (Natural Forest-NF; Exotic Woodland-EW and Production Forest-PF) in three of the Azores islands. We used biodiversity indices and found that NF harbored the highest plant diversity levels and PF the lowest. Diversity levels were lower for structural than for taxonomic data, particularly for PF. The highest tree carbon stock was found at EW in one of the islands, while PF consistently exhibited relatively high tree carbon stocks in the three islands. The largest soil carbon stocks were found at EW, while leaf litter carbon stocks were higher at PF. We concluded that NF play a fundamental role as plant diversity hotspots but have lower relevance as carbon stocks what might be associated with montane environmental conditions. PFs provide economic assets and act as carbon sinks, while EWs play a major role as carbon sinks in soil, but also at tree level in the oldest forests.

Identifying Forest Structural Types along an Aridity Gradient in Peninsular Spain: Integrating Low-Density LiDAR, Forest Inventory, and Aridity Index

Forest structure is a key driver of forest functional processes. The characterization of forest structure across spatiotemporal scales is essential for forest monitoring and management. LiDAR data have proven particularly useful for cost-effectively estimating forest structural attributes. This paper evaluates the ability of combined forest inventory data and low-density discrete return airborne LiDAR data to discriminate main forest structural types in the Mediterranean-temperate transition ecotone. Firstly, we used six structural variables from the Spanish National Forest Inventory (SNFI) and an aridity index in a k-medoids algorithm to define the forest structural types. These variables were calculated for 2770 SNFI plots. We identified the main species for each structural type using the SNFI. Secondly, we developed a Random Forest model to predict the spatial distribution of structural types and create wall-to-wall maps from LiDAR data. The k-medoids clustering algorithm enabled the identification of four clusters of forest structures. A total of six out of forty-one potential LiDAR metrics were utilized in our Random Forest, after evaluating their importance in the Random Forest model. Selected metrics were, in decreasing order of importance, the percentage of all returns above 2 m, mean height of the canopy profile, the difference between the 90th and 50th height percentiles, the area under the canopy curve, and the 5th and the 95th percentile of the return heights. The model yielded an overall accuracy of 64.18%. The producer’s accuracy ranged between 36.11% and 88.93%. Our results confirm the potential of this approximation for the continuous monitoring of forest structures, which is key to guiding forest management in this region.

Remarkable Effects of Urbanization on Forest Landscape Multifunctionality in Urban Peripheries: Evidence from Liaoyuan City in Northeast China

Forest landscape multifunctionality (FLM) provides multiple benefits, such as climate regulation, water storage, and biodiversity maintenance. However, the external factors limiting FLM have not been fully identified, although addressing them could contribute to sustainable development. The present study aimed to identify and quantify the role of urbanization as an external factor that affects FLM. To this end, impervious area changes in Liaoyuan, China, were observed from 2000 to 2018, and 10 buffer zones at 500 m intervals were established outside the city. Within each buffer zone, we analyzed changes in forest landscape functions, including habitat maintenance, carbon sequestration, and water yield, as well as changes in the multifunctionality of their composition. The urbanization of Liaoyuan was significant in 2000–2018. The functions of the forest landscape became stronger and more stable as they were located further away from the urban edge. We refer to this pattern as the gradient effect of urbanization. Specifically, urbanization affected the investigated functions at a distance of 1000–2500 m. The FLM showed a more significant gradient effect of urbanization. The impact distance of urbanization on the FLM increased from 3000 m in 2000 to over 5000 m in 2018. This impact distance increased significantly whenever urbanization strengthened significantly (i.e., in 2005–2010 and 2015–2018). These findings are instructive for forest and urban managers working to achieve multiple Sustainable Development Goals.

Bat responses to changes in forest composition and prey abundance depend on landscape matrix and stand structure

Despite the key importance of the landscape matrix for bats, we still not fully understand how the effect of forest composition interacts at combined stand and landscape scales to shape bat communities. In addition, we lack detailed knowledge on the effects of local habitat structure on bat-prey relationships in forested landscapes. We tested the assumptions that (i) forest composition has interacting effects on bats between stand and landscape scales; and (ii) stand structure mediates prey abundance effects on bat activity. Our results indicated that in conifer-dominated landscapes (> 80% of coniferous forests) bat activity was higher in stands with a higher proportion of deciduous trees while bats were less active in stands with a higher proportion of deciduous trees in mixed forest landscapes (~ 50% of deciduous forests). Moth abundance was selected in the best models for six among nine bat species. The positive effect of moth abundance on Barbastella barbastellus was mediated by vegetation clutter, with dense understory cover likely reducing prey accessibility. Altogether, our findings deepen our understanding of the ecological processes affecting bats in forest landscapes and strengthen the need to consider both landscape context and trophic linkage when assessing the effects of stand-scale compositional and structural attributes on bats.

Application of the Conservation Planning Tool Zonation to Inform Retention Planning in the Boreal Forest of Western Canada

Retention forestry is an approach in which live trees and other components of forest structure are retained within harvested areas. A primary objective of retention forestry is to maintain biodiversity and to hasten post-harvest recovery of forest structure and function. Retention is now a key element in sustainable forest management practices in many regions of the world. However, locating where retention should be placed to best achieve management objectives is a challenging problem, and evidence-based approaches to operational applications are rare. We suggest here that harvest planners could benefit from the use of systematic conservation planning principles and methods to inform retention design. Specifically, we used a conservation planning—or prioritization—tool, Zonation, to create spatially-explicit scenarios of retention harvesting in a boreal mixedwood forest in northwestern Alberta, Canada. Scenarios were informed by several environmental variables related to site productivity; in particular, we used a metric of wetness (depth-to-water from the Wet Areas Mapping algorithm) that is based on airborne lidar-derived terrain models previously shown to correlate with patterns in post-harvest forest regeneration and biodiversity. The nine retention scenarios examined here related to the placement of retention focused to drier, mesic, or wetter sites in combination with other prioritization constraints. Results were compared with an existing harvest plan to assess differences in the spatial pattern of retention (e.g., percent overlapping area, number of patches, size of the patches). We also tested for the homogeneity of forest attributes (e.g., tree species, deciduous density) between scenarios and the existing harvest plan using multivariate dispersion analysis. Our results showed limited commonalities among scenarios compared to the existing harvest plan; they were characterized as having limited spatial overlap, and more and smaller patches with the use of a timber-cost constraint further affecting retention patterns. While modeling results significantly differed from current retention practices, the approach presented here offers flexibility in testing different scenarios and assessing trade-offs between timber production and conservation goals using a standardized conservation planning toolkit.

Resilience of Spanish forests to recent droughts and climate change

A widespread increase in forest cover is underway in northern Mediterranean forests because of land abandonment and decreased wood demand, but the resilience of these successional forests to climate change remains unresolved. Here we use 18-year time series of canopy greenness derived from satellite imagery (NDVI) to evaluate the impacts of climate change on Spain's forests. Specifically, we analyzed how NDVI was influenced by the climatic water balance (i.e. Standardized Precipitation-Evapotranspiration Index, SPEI), using monthly time-series extracted from 3,100 pixels of forest, categorized into ten forest types. The forests increased in leaf area index by 0.01 per year on average (from 1.7 in 2000 to 1.9 in 2017) but there was enormous variation among years related to climatic water balance. Forest types varied in response to drought events: those dominated by drought-avoiding species showed strong covariance between greenness and SPEI, while those dominated by drought-tolerant species showed weak covariance. Native forests usually recovered more than 80% of greenness within the 18 months and the remainder within 5 years, but plantations of Eucalyptus were less resilient. Management to increase the resilience of forests-a key goal of forestry in the Mediterranean region-appears to have had a positive effect: canopy greenness within protected forests was more resilient to drought than within non-protected forests. In conclusion, many of Spain's successional forests have been resilient to drought over the past 18 years, from the perspective of space. Future studies will need to combine remote sensing with field-based analyses of physiological tolerances and mortality processes to understand how Mediterranean forests will respond to the rapid climate change predicted for this region in the coming decades.

Neighbourhood-mediated shifts in tree biomass allocation drive overyielding in tropical species mixtures

Variations in crown forms promote canopy space-use and productivity in mixed-species forests. However, we have a limited understanding on how this response is mediated by changes in within-tree biomass allocation. Here, we explored the role of changes in tree allometry, biomass allocation and architecture in shaping diversity-productivity relationships (DPRs) in the oldest tropical tree diversity experiment. We conducted whole-tree destructive biomass measurements and terrestrial laser scanning. Spatially explicit models were built at the tree level to investigate the effects of tree size and local neighbourhood conditions. Results were then upscaled to the stand level, and mixture effects were explored using a bootstrapping procedure. Biomass allocation and architecture substantially changed in mixtures, which resulted from both tree-size effects and neighbourhood-mediated plasticity. Shifts in biomass allocation among branch orders explained substantial shares of the observed overyielding. By contrast, root-to-shoot ratios, as well as the allometric relationships between tree basal area and aboveground biomass, were little affected by the local neighbourhood. Our results suggest that generic allometric equations can be used to estimate forest aboveground biomass overyielding from diameter inventory data. Overall, we demonstrate that shifts in tree biomass allocation are mediated by the local neighbourhood and promote DPRs in tropical forests.

Blind spots in global soil biodiversity and ecosystem function research

Soils harbor a substantial fraction of the world's biodiversity, contributing to many crucial ecosystem functions. It is thus essential to identify general macroecological patterns related to the distribution and functioning of soil organisms to support their conservation and consideration by governance. These macroecological analyses need to represent the diversity of environmental conditions that can be found worldwide. Here we identify and characterize existing environmental gaps in soil taxa and ecosystem functioning data across soil macroecological studies and 17,186 sampling sites across the globe. These data gaps include important spatial, environmental, taxonomic, and functional gaps, and an almost complete absence of temporally explicit data. We also identify the limitations of soil macroecological studies to explore general patterns in soil biodiversity-ecosystem functioning relationships, with only 0.3% of all sampling sites having both information about biodiversity and function, although with different taxonomic groups and functions at each site. Based on this information, we provide clear priorities to support and expand soil macroecological research.

The reclaimed coal mine ecosystem diverges from the surrounding ecosystem and reaches a new self-sustaining state after 20-23 years of succession in the Loess Plateau area, China

Whether or not the completely destroyed ecosystem would follow a succession trajectory towards the surrounding forest ecosystem after restoration remains debatable. Here, a comprehensive dataset of thirty-five ecosystem functions were measured on five reclaimed opencast coal mine forest plots (two Robinia pseudoacacia - Pinus tabuliformis mixed forests with different technosol conditions: RP_T and RP_M; one R. pseudoacacia - Ulmus pumila - Ailanthus altissima mixed forest: RUA; one Picea meyeri - Picea wilsonii - Hippophae rhamnoides mixed forest: PPH; one R. pseudoacacia monoculture forest: RM) and one natural forest plot (Populus simonii monoculture forest: PM) in Pingshuo opencast coal mine, Shanxi Province, China. These functions were employed to examine the reclamation effects among plots in terms of four management scenarios (i.e., biomass productivity, carbon sequestration, general biodiversity conservation and nutrient accumulation) and to determine the affinities between reclaimed ecosystem and its native counterpart according to Bray Curtis Distance Algorithm. The results showed that after 20-23 years of succession, thirty-five ecosystem functions differed among plots and eight highest ecosystem functions were found in RP_T and RP_M (P < 0.05). In scenarios of biomass productivity, carbon sequestration, general biodiversity conservation and nutrient accumulation scenarios, RP_M scored 0.645, 0.470, 0.467 and 0.578, respectively. Accordingly, RP_T scored 0.458, 0.447, 0.405 and 0.515, respectively. Consequently, RP_T and RP_M had the highest scores in the four management objectives compared to other plots. With regard to the ecosystem affinities, RP_T was analogous to RP_M. Similarly, RUA and PPH resembled PM and RM, respectively. Overall, the pioneer tree species determined the reclamation quality in the dump in the long term, which could be referred as the "trigger effect" in the succession trajectory in reclaimed ecosystem, and the reclaimed coal mine ecosystem diverged from the surrounding ecosystem and reached a new self-sustaining state after 20-23 years of succession.

Plant diversity and soil stoichiometry regulates the changes in multifunctionality during pine temperate forest secondary succession

The shift in ecosystem multifunctionality during ecosystem succession (years to decades) remains largely unexplored. In this study, we used a 120-year-old pine temperate forest chronosequence (1: 1-19 years, stage 2: 20-39 years, stage 3: 40-59 years, stage 4: 60-79 years, stage 5: 80-99 years, stage 6: 100-120 years) to evaluate the role that time plays in shaping ecosystem multifunctionality (nutrient cycling, carbon stocks, water regulation, decomposition and wood production), and found that, over the first century, ecosystem functioning gradually increased every ~50 years. Such a result was maintained for individual groups of ecosystem functions and services including nutrient cycling, carbon stocks, decomposition and wood production. Plant diversity and soil stoichiometry (C:N ratio) were the major environmental predictors for the changes in ecosystem multifunctionality during forest secondary succession. Plant diversity increased during ecosystem succession and was positively related to ecosystem multifunctionality. The soil C:N ratio decreased during ecosystem succession and was negatively related to multifunctionality. Our results suggest that increases in aboveground resource heterogeneity (higher plant diversity) and organic matter quality (lower soil C:N ratios) could help explain the increases in multifunctionality over a century of forest development. Our work illustrates the importance of time in shaping multifunctionality during the first century of ecosystem succession, and further provide important insights for the management of temperate forest ecosystems.

A multitrophic perspective on biodiversity-ecosystem functioning research

Concern about the functional consequences of unprecedented loss in biodiversity has prompted biodiversity-ecosystem functioning (BEF) research to become one of the most active fields of ecological research in the past 25 years. Hundreds of experiments have manipulated biodiversity as an independent variable and found compelling support that the functioning of ecosystems increases with the diversity of their ecological communities. This research has also identified some of the mechanisms underlying BEF relationships, some context-dependencies of the strength of relationships, as well as implications for various ecosystem services that mankind depends upon. In this paper, we argue that a multitrophic perspective of biotic interactions in random and non-random biodiversity change scenarios is key to advance future BEF research and to address some of its most important remaining challenges. We discuss that the study and the quantification of multitrophic interactions in space and time facilitates scaling up from small-scale biodiversity manipulations and ecosystem function assessments to management-relevant spatial scales across ecosystem boundaries. We specifically consider multitrophic conceptual frameworks to understand and predict the context-dependency of BEF relationships. Moreover, we highlight the importance of the eco-evolutionary underpinnings of multitrophic BEF relationships. We outline that FAIR data (meeting the standards of findability, accessibility, interoperability, and reusability) and reproducible processing will be key to advance this field of research by making it more integrative. Finally, we show how these BEF insights may be implemented for ecosystem management, society, and policy. Given that human well-being critically depends on the multiple services provided by diverse, multitrophic communities, integrating the approaches of evolutionary ecology, community ecology, and ecosystem ecology in future BEF research will be key to refine conservation targets and develop sustainable management strategies.

Thorny Shrubs Limit the Browsing Pressure of Large Herbivores on Tree Regeneration in Temperate Lowland Forested Landscapes

Thorny shrubs are considered as an important driver in the natural development of temperate forest structures, particularly in European lowland forests. We assessed the current role of thorny shrubs in the regeneration of deciduous tree species under heavy browsing pressure in a central European temperate forested landscape. The study’s military training area offered a unique opportunity to investigate the processes in which deciduous tree seedlings grew under thorny shrubs and in the close vicinity of thorny shrubs in a landscape with a high density of large herbivores (red deer and sika deer). We assessed the number of seedlings, species composition, seedling height, and degree of browsing damage, and their relationship to study plots elevation, thorny shrub species, coverage, and height. The regenerated tree seedlings were mostly detected as common ash (Fraxinus excelsior) and wild cherry (Cerasus avium). The species of thorny shrubs were blackthorn (Prunus spinosa), hawthorn (Crataegus sp.), and wild rose (Rosa sp.). We found that the thorny shrubs protected the tree seedlings from browsers to a large extent. However, the effects of thorny shrubs on the tree seedlings' characteristics varied among the shrub species. While results revealed significant effects of hawthorn and wild rose on the tree seedlings' abundance and survival, blackthorn’s negative effect of shading the tree seedlings outweighed its protective role. These results indicated a possible mechanism that enabled the regeneration of deciduous tree species under large herbivore pressure. These results can be applied in the landscape planning and forest management of deciduous tree regeneration and forest restoration in temperate forested lowland landscapes, where high densities of large herbivores (without the presence of large predators) usually occur.

Functional diversity drives ecosystem multifunctionality in a Pinus yunnanensis natural secondary forest

It is essential to understand how the loss of biodiversity impacts both ecosystem function (EF) and multifunctionality (EMF). Previous studies have mostly focused on predicting how species richness (SR) impacts EMF, while the effect of functional diversity (FD) on EMF remains unclear. Specifically, we know little about the primary functional drivers impacting EMF compared with SR. Therefore, we analysed 8 ecosystem functions within 58 natural secondary forest plots to investigate the effect of FD on both individual EF and EMF. Our results suggest that SR and FD had very significant positive effects on plant phosphorus, soil available phosphorus, and soil total nitrogen. FD explained significantly more variations in these functional responses than SR for individual ecosystem functioning. We also used a multiple threshold approach to test the effect of SR and FD on EMF. We found that FD and SR were positively related to EMF regardless of whether low-level function or high-level function was desired, but FD had a larger effect than SR. Based on the averaging approach, OLS regression, multivariate linear regression model and random forest analysis, we found that SR and FD were both drivers of EMF but that FD had a stronger effect and could explain more variation. As such, we conclude that FD drives ecosystem multifunctionality more than SR.

A Multifactorial Approach to Value Supporting Ecosystem Services in Spanish Forests and Its Implications in a Warming World

Carbon storage and sequestration are key ecosystem services critical to human well-being and biodiversity conservation. In a warming context, the quantification and valuation of carbon storage and sequestration is important in ensuring that effective incentives are put in place to tackle climate change. The quantification and valuation of ES such as carbon storage and sequestration requires the calculus of actual values and prediction, however, it usually does not include key processes that can indirectly influence carbon dynamics (i.e., risk, conservation or management). Here, we define a multifactorial approach to value ecosystem services based on two stages: (1) a biophysical approximation that integrates yearly supporting ecosystem services (i.e., quantification of carbon storage and sequestration) and (2) a weighing approach including factors that indirectly influence carbon storage and sequestration or that deserve specific attention (i.e., risk, conservation or management factors). The quantification of carbon storage and sequestration indicated that Spanish forests store on average 43 Mg C ha−1 and sequestrate on average 1.02 Mg C ha−1 year−1. Forest structure was a strong determinant of carbon storage and sequestration in Iberian forests, hence there was a strong spatial variation in the carbon sink. We adapted the weighting values to a financial cap and the monetary value of carbon increased more than four times when the weighting factors were taken into account. Finally, we argue that a multifactorial approach to value supporting ecosystem services incorporating aspects related to conservation and risk prevention can facilitate ecosystem service valuation and assist policy makers and stakeholders to establish payment service policies.

Biodiversity across trophic levels drives multifunctionality in highly diverse forests

Human-induced biodiversity change impairs ecosystem functions crucial to human well-being. However, the consequences of this change for ecosystem multifunctionality are poorly understood beyond effects of plant species loss, particularly in regions with high biodiversity across trophic levels. Here we adopt a multitrophic perspective to analyze how biodiversity affects multifunctionality in biodiverse subtropical forests. We consider 22 independent measurements of nine ecosystem functions central to energy and nutrient flow across trophic levels. We find that individual functions and multifunctionality are more strongly affected by the diversity of heterotrophs promoting decomposition and nutrient cycling, and by plant functional-trait diversity and composition, than by tree species richness. Moreover, cascading effects of higher trophic-level diversity on functions originating from lower trophic-level processes highlight that multitrophic biodiversity is key to understanding drivers of multifunctionality. A broader perspective on biodiversity-multifunctionality relationships is crucial for sustainable ecosystem management in light of non-random species loss and intensified biotic disturbances under future environmental change.

Biodiversity change can impact ecosystem functioning, though this is primarily studied at lower trophic levels. Here, Schuldt et al. find that biodiversity components other than tree species richness are particularly important, and higher trophic level diversity plays a role in multifunctionality.