Author Correction: Environment and host as large-scale controls of ectomycorrhizal fungi.

Change history: In the HTML version of this Article, author 'Filipa Cox' had no affiliation in the author list, although she was correctly associated with affiliation 3 in the PDF. In addition, the blue circles for 'oak' were missing from Extended Data Fig. 1. These errors have been corrected online.

Environment and host as large-scale controls of ectomycorrhizal fungi.

Explaining the large-scale diversity of soil organisms that drive biogeochemical processes-and their responses to environmental change-is critical. However, identifying consistent drivers of belowground diversity and abundance for some soil organisms at large spatial scales remains problematic. Here we investigate a major guild, the ectomycorrhizal fungi, across European forests at a spatial scale and resolution that is-to our knowledge-unprecedented, to explore key biotic and abiotic predictors of ectomycorrhizal diversity and to identify dominant responses and thresholds for change across complex environmental gradients. We show the effect of 38 host, environment, climate and geographical variables on ectomycorrhizal diversity, and define thresholds of community change for key variables. We quantify host specificity and reveal plasticity in functional traits involved in soil foraging across gradients. We conclude that environmental and host factors explain most of the variation in ectomycorrhizal diversity, that the environmental thresholds used as major ecosystem assessment tools need adjustment and that the importance of belowground specificity and plasticity has previously been underappreciated.

Responses of soil fauna communities to the individual and combined effects of multiple global change factors.

Soil fauna plays a key role in regulating biogeochemical cycles, but how multiple global change factors (GCFs) may affect faunal communities remains poorly studied. We conducted a meta-analysis using 1154 observations to evaluate the individual and combined effects of elevated CO2 , nitrogen (N) addition, warming, increased rainfall and drought on soil fauna density and diversity. Here we show that, overall, individual and combined effects of GCFs had negligible effects on soil fauna density and diversity, except that density was negatively affected by drought (-27.4%) and positively affected by increased rainfall individually (+24.9%) and in combination with N addition (+67.3%) or warming (+70.4%). GCF effects varied among taxonomic groups both in magnitude and direction. Variables such as latitude, elevation and experimental setting significantly impacted both individual and combined effects. Our results suggest that soil fauna density is affected by changed rainfall regimes, while diversity is resistant against individual and combined effects of multiple GCFs.

Biotic homogenization can decrease landscape-scale forest multifunctionality.

Many experiments have shown that local biodiversity loss impairs the ability of ecosystems to maintain multiple ecosystem functions at high levels (multifunctionality). In contrast, the role of biodiversity in driving ecosystem multifunctionality at landscape scales remains unresolved. We used a comprehensive pan-European dataset, including 16 ecosystem functions measured in 209 forest plots across six European countries, and performed simulations to investigate how local plot-scale richness of tree species (α-diversity) and their turnover between plots (ß-diversity) are related to landscape-scale multifunctionality. After accounting for variation in environmental conditions, we found that relationships between α-diversity and landscape-scale multifunctionality varied from positive to negative depending on the multifunctionality metric used. In contrast, when significant, relationships between ß-diversity and landscape-scale multifunctionality were always positive, because a high spatial turnover in species composition was closely related to a high spatial turnover in functions that were supported at high levels. Our findings have major implications for forest management and indicate that biotic homogenization can have previously unrecognized and negative consequences for large-scale ecosystem multifunctionality.

The response of soil solution chemistry in European forests to decreasing acid deposition.

Acid deposition arising from sulphur (S) and nitrogen (N) emissions from fossil fuel combustion and agriculture has contributed to the acidification of terrestrial ecosystems in many regions globally. However, in Europe and North America, S deposition has greatly decreased in recent decades due to emissions controls. In this study, we assessed the response of soil solution chemistry in mineral horizons of European forests to these changes. Trends in pH, acid neutralizing capacity (ANC), major ions, total aluminium (Altot ) and dissolved organic carbon were determined for the period 1995-2012. Plots with at least 10 years of observations from the ICP Forests monitoring network were used. Trends were assessed for the upper mineral soil (10-20 cm, 104 plots) and subsoil (40-80 cm, 162 plots). There was a large decrease in the concentration of sulphate (SO42-) in soil solution; over a 10-year period (2000-2010), SO42- decreased by 52% at 10-20 cm and 40% at 40-80 cm. Nitrate was unchanged at 10-20 cm but decreased at 40-80 cm. The decrease in acid anions was accompanied by a large and significant decrease in the concentration of the nutrient base cations: calcium, magnesium and potassium (Bc = Ca2+  + Mg2+  + K+ ) and Altot over the entire dataset. The response of soil solution acidity was nonuniform. At 10-20 cm, ANC increased in acid-sensitive soils (base saturation ≤10%) indicating a recovery, but ANC decreased in soils with base saturation >10%. At 40-80 cm, ANC remained unchanged in acid-sensitive soils (base saturation ≤20%, pHCaCl2 ≤ 4.5) and decreased in better-buffered soils (base saturation >20%, pHCaCl2 > 4.5). In addition, the molar ratio of Bc to Altot either did not change or decreased. The results suggest a long-time lag between emission abatement and changes in soil solution acidity and underline the importance of long-term monitoring in evaluating ecosystem response to decreases in deposition.

Tree species diversity affects decomposition through modified micro-environmental conditions across European forests.

Different tree species influence litter decomposition directly through species-specific litter traits, and indirectly through distinct modifications of the local decomposition environment. Whether these indirect effects on decomposition are influenced by tree species diversity is presently not clear. We addressed this question by studying the decomposition of two common substrates, cellulose paper and wood sticks, in a total of 209 forest stands of varying tree species diversity across six major forest types at the scale of Europe. Tree species richness showed a weak but positive correlation with the decomposition of cellulose but not with that of wood. Surprisingly, macroclimate had only a minor effect on cellulose decomposition and no effect on wood decomposition despite the wide range in climatic conditions among sites from Mediterranean to boreal forests. Instead, forest canopy density and stand-specific litter traits affected the decomposition of both substrates, with a particularly clear negative effect of the proportion of evergreen tree litter. Our study suggests that species richness and composition of tree canopies modify decomposition indirectly through changes in microenvironmental conditions. These canopy-induced differences in the local decomposition environment control decomposition to a greater extent than continental-scale differences in macroclimatic conditions.

Does species richness affect fine root biomass and production in young forest plantations?

Tree species diversity has been reported to increase forest ecosystem above-ground biomass and productivity, but little is known about below-ground biomass and production in diverse mixed forests compared to single-species forests. For testing whether species richness increases below-ground biomass and production and thus complementarity between forest tree species in young stands, we determined fine root biomass and production of trees and ground vegetation in two experimental plantations representing gradients in tree species richness. Additionally, we measured tree fine root length and determined species composition from fine root biomass samples with the near-infrared reflectance spectroscopy method. We did not observe higher biomass or production in mixed stands compared to monocultures. Neither did we observe any differences in tree root length or fine root turnover. One reason for this could be that these stands were still young, and canopy closure had not always taken place, i.e. a situation where above- or below-ground competition did not yet exist. Another reason could be that the rooting traits of the tree species did not differ sufficiently to support niche differentiation. Our results suggested that functional group identity (i.e. conifers vs. broadleaved species) can be more important for below-ground biomass and production than the species richness itself, as conifers seemed to be more competitive in colonising the soil volume, compared to broadleaved species.

Plant movements and climate warming: intraspecific variation in growth responses to nonlocal soils.

Most range shift predictions focus on the dispersal phase of the colonization process. Because moving populations experience increasingly dissimilar nonclimatic environmental conditions as they track climate warming, it is also critical to test how individuals originating from contrasting thermal environments can establish in nonlocal sites. We assess the intraspecific variation in growth responses to nonlocal soils by planting a widespread grass of deciduous forests (Milium effusum) into an experimental common garden using combinations of seeds and soil sampled in 22 sites across its distributional range, and reflecting movement scenarios of up to 1600 km. Furthermore, to determine temperature and forest-structural effects, the plants and soils were experimentally warmed and shaded. We found significantly positive effects of the difference between the temperature of the sites of seed and soil collection on growth and seedling emergence rates. Migrant plants might thus encounter increasingly favourable soil conditions while tracking the isotherms towards currently 'colder' soils. These effects persisted under experimental warming. Rising temperatures and light availability generally enhanced plant performance. Our results suggest that abiotic and biotic soil characteristics can shape climate change-driven plant movements by affecting growth of nonlocal migrants, a mechanism which should be integrated into predictions of future range shifts.

Afforestation effects on SOC in former cropland: oak and spruce chronosequences resampled after 13 years.

Chronosequences are commonly used to assess soil organic carbon (SOC) sequestration after land-use change, but SOC dynamics predicted by this space-for-time substitution approach have rarely been validated by resampling. We conducted a combined chronosequence/resampling study in a former cropland area (Vestskoven) afforested with oak (Quercus robur) and Norway spruce (Picea abies) over the past 40 years. The aims of this study were (i) to compare present and previous chronosequence trends in forest floor and top mineral soil (0-25 cm) C stocks; (ii) to compare chronosequence estimates with current rates of C stock change based on resampling at the stand level; (iii) to estimate SOC changes in the subsoil (25-50 cm); and (iv) to assess the influence of two tree species on SOC dynamics. The two chronosequence trajectories for forest floor C stocks revealed consistently higher rates of C sequestration in spruce than oak. The chronosequence trajectory was validated by resampling and current rates of forest floor C sequestration decreased with stand age. Chronosequence trends in topsoil SOC in 2011 did not differ significantly from those reported in 1998, however, there was a shift from a negative rate (1998: -0.3 Mg C ha(-1) yr(-1) ) to no change in 2011. In contrast SOC stocks in the subsoil increased with stand age, however, not significantly (P = 0.1), suggesting different C dynamics in and below the former plough layer. Current rates of C change estimated by repeated sampling decreased with stand age in forest floors but increased in the topsoil. The contrasting temporal change in forest floor and mineral soil C sequestration rates indicate a shift in C source-sink strength after approximately 40 years. We conclude that afforestation of former cropland within the temperate region may induce soil C loss during the first decades followed by a recovery phase of yet unknown duration.

Fungal community composition predicts forest carbon storage at a continental scale.

Forest soils harbor hyper-diverse microbial communities which fundamentally regulate carbon and nutrient cycling across the globe. Directly testing hypotheses on how microbiome diversity is linked to forest carbon storage has been difficult, due to a lack of paired data on microbiome diversity and in situ observations of forest carbon accumulation and storage. Here, we investigated the relationship between soil microbiomes and forest carbon across 238 forest inventory plots spanning 15 European countries. We show that the composition and diversity of fungal, but not bacterial, species is tightly coupled to both forest biotic conditions and a seven-fold variation in tree growth rates and biomass carbon stocks when controlling for the effects of dominant tree type, climate, and other environmental factors. This linkage is particularly strong for symbiotic endophytic and ectomycorrhizal fungi known to directly facilitate tree growth. Since tree growth rates in this system are closely and positively correlated with belowground soil carbon stocks, we conclude that fungal composition is a strong predictor of overall forest carbon storage across the European continent.

Dramatic changes in ectomycorrhizal community composition, root tip abundance and mycelial production along a stand-scale nitrogen deposition gradient.

• Nitrogen (N) availability is known to influence ectomycorrhizal fungal components, such as fungal community composition, biomass of root tips and production of mycelia, but effects have never been demonstrated within the same forest. • We measured concurrently the abundance of ectomycorrhizal root tips and the production of external mycelia, and explored the changes in the ectomycorrhizal community composition, across a stand-scale N deposition gradient (from 27 to 43 kg N ha⁻¹ yr⁻¹) at the edge of a spruce forest. The N status was affected along the gradient as shown by a range of N availability indices. • Ectomycorrhizal root tip abundance and mycelial production decreased five and 10-fold, respectively, with increasing N deposition. In addition, the ectomycorrhizal fungal community changed and the species richness decreased. The changes were correlated with the measured indices of N status, in particular N deposition and N leaching. • The relationship between the altered ectomycorrhizal community, root tip abundance and mycelial production is discussed in the context of the N parameters. We suggest that increased N deposition to forests will cause large changes in ectomycorrhizal fungal community structure and functioning, which, in turn, may result in reduced N uptake by roots and fungi, and increased losses of N by leaching.

Four decades of post-agricultural forest development have caused major redistributions of soil phosphorus fractions.

Fertilisation of agricultural land causes an accumulation of nutrients in the top soil layer, among which phosphorus (P) is particularly persistent. Changing land use from farmland to forest affects soil properties, but changes in P pools have rarely been studied despite their importance to forest ecosystem development. Here, we describe the redistributions of the P pools in a four-decadal chronosequence of post-agricultural common oak (Quercus robur L.) forests in Belgium and Denmark. The aim was to assess whether forest age causes a repartitioning of P throughout the various soil P pools (labile P, slowly cycling P and occluded P); in particular, we addressed the time-related alterations in the inorganic versus organic P fractions. In less than 40 years of oak forest development, significant redistributions have occurred between different P fractions. While both the labile and the slowly cycling inorganic P fractions significantly decreased with forest age, the organic fractions significantly increased. The labile P pool (inorganic + organic), which is considered to be the pool of P most likely to contribute to plant-available P, significantly decreased with forest age (from >20 to <10% of total P), except in the 0-5 cm of topsoil, where labile P remained persistently high. The shift from inorganic to organic P and the shifts between the different inorganic P fractions are driven by biological processes and also by physicochemical changes related to forest development. It is concluded that the organic labile P fraction, which is readily mineralisable, should be taken into account when studying the bioavailable P pool in forest ecosystems.

Climatic conditions, not above- and belowground resource availability and uptake capacity, mediate tree diversity effects on productivity and stability.

Tree species diversity promotes multiple ecosystem functions and services. However, little is known about how above- and belowground resource availability (light, nutrients, and water) and resource uptake capacity mediate tree species diversity effects on aboveground wood productivity and temporal stability of productivity in European forests and whether the effects differ between humid and arid regions. We used the data from six major European forest types along a latitudinal gradient to address those two questions. We found that neither leaf area index (a proxy for light uptake capacity), nor fine root biomass (a proxy for soil nutrient and water uptake capacity) was related to tree species richness. Leaf area index did, however, enhance productivity, but negatively affected stability. Productivity was further promoted by soil nutrient availability, while stability was enhanced by fine root biomass. We only found a positive effect of tree species richness on productivity in arid regions and a positive effect on stability in humid regions. This indicates a possible disconnection between productivity and stability regarding tree species richness effects. In other words, the mechanisms that drive the positive effects of tree species richness on productivity do not per se benefit stability simultaneously. Our findings therefore suggest that tree species richness effects are largely mediated by differences in climatic conditions rather than by differences in above- and belowground resource availability and uptake capacity at the regional scales.

Simulating tree growth response to climate change in structurally diverse oak and beech forests.

This study aimed to simulate oak and beech forest growth under various scenarios of climate change and to evaluate how the forest response depends on site properties and particularly on stand characteristics using the individual process-based model HETEROFOR. First, this model was evaluated on a wide range of site conditions. We used data from 36 long-term forest monitoring plots to initialize, calibrate, and evaluate HETEROFOR. This evaluation showed that HETEROFOR predicts individual tree radial growth and height increment reasonably well under different growing conditions when evaluated on independent sites. In our simulations under constant CO2 concentration ([CO2]cst) for the 2071-2100 period, climate change induced a moderate net primary production (NPP) gain in continental and mountainous zones and no change in the oceanic zone. The NPP changes were negatively affected by air temperature during the vegetation period and by the annual rainfall decrease. To a lower extent, they were influenced by soil extractable water reserve and stand characteristics. These NPP changes were positively affected by longer vegetation periods and negatively by drought for beech and larger autotrophic respiration costs for oak. For both species, the NPP gain was much larger with rising CO2 concentration ([CO2]var) mainly due to the CO2 fertilisation effect. Even if the species composition and structure had a limited influence on the forest response to climate change, they explained a large part of the NPP variability (44% and 34% for [CO2]cst and [CO2]var, respectively) compared to the climate change scenario (5% and 29%) and the inter-annual climate variability (20% and 16%). This gives the forester the possibility to act on the productivity of broadleaved forests and prepare them for possible adverse effects of climate change by reinforcing their resilience.

Forest tree growth is linked to mycorrhizal fungal composition and function across Europe.

Most trees form symbioses with ectomycorrhizal fungi (EMF) which influence access to growth-limiting soil resources. Mesocosm experiments repeatedly show that EMF species differentially affect plant development, yet whether these effects ripple up to influence the growth of entire forests remains unknown. Here we tested the effects of EMF composition and functional genes relative to variation in well-known drivers of tree growth by combining paired molecular EMF surveys with high-resolution forest inventory data across 15 European countries. We show that EMF composition was linked to a three-fold difference in tree growth rate even when controlling for the primary abiotic drivers of tree growth. Fast tree growth was associated with EMF communities harboring high inorganic but low organic nitrogen acquisition gene proportions and EMF which form contact versus medium-distance fringe exploration types. These findings suggest that EMF composition is a strong bio-indicator of underlying drivers of tree growth and/or that variation of forest EMF communities causes differences in tree growth. While it may be too early to assign causality or directionality, our study is one of the first to link fine-scale variation within a key component of the forest microbiome to ecosystem functioning at a continental scale.

Highly comparable metabarcoding results from MGI-Tech and Illumina sequencing platforms.

With the developments in DNA nanoball sequencing technologies and the emergence of new platforms, there is an increasing interest in their performance in comparison with the widely used sequencing-by-synthesis methods. Here, we test the consistency of metabarcoding results from DNBSEQ-G400RS (DNA nanoball sequencing platform by MGI-Tech) and NovaSeq 6000 (sequencing-by-synthesis platform by Illumina) platforms using technical replicates of DNA libraries that consist of COI gene amplicons from 120 soil DNA samples. By subjecting raw sequencing data from both platforms to a uniform bioinformatics processing, we found that the proportion of high-quality reads passing through the filtering steps was similar in both datasets. Per-sample operational taxonomic unit (OTU) and amplicon sequence variant (ASV) richness patterns were highly correlated, but sequencing data from DNBSEQ-G400RS harbored a higher number of OTUs. This may be related to the lower dominance of most common OTUs in DNBSEQ data set (thus revealing higher richness by detecting rare taxa) and/or to a lower effective read quality leading to generation of spurious OTUs. However, there was no statistical difference in the ASV and post-clustered ASV richness between platforms, suggesting that additional denoising step in the ASV workflow had effectively removed the 'noisy' reads. Both OTU-based and ASV-based composition were strongly correlated between the sequencing platforms, with essentially interchangeable results. Therefore, we conclude that DNBSEQ-G400RS and NovaSeq 6000 are both equally efficient high-throughput sequencing platforms to be utilized in studies aiming to apply the metabarcoding approach, but the main benefit of the former is related to lower sequencing cost.

Stability of buried carbon in deep-ploughed forest and cropland soils - implications for carbon stocks.

Accumulation of soil organic carbon (SOC) may play a key role in climate change mitigation and adaptation. In particular, subsoil provides a great potential for additional SOC storage due to the assumed higher stability of subsoil SOC. The fastest way in which SOC reaches the subsoil is via burial, e.g. via erosion or deep ploughing. We assessed the effect of active SOC burial through deep ploughing on long-term SOC stocks and stability in forest and cropland subsoil. After 25-48 years, deep-ploughed subsoil contained significantly more SOC than reference subsoils, in both forest soil (+48%) and cropland (+67%). However, total SOC stocks down to 100 cm in deep-ploughed soil were greater than in reference soil only in cropland, and not in forests. This was explained by slower SOC accumulation in topsoil of deep-ploughed forest soils. Buried SOC was on average 32% more stable than reference SOC, as revealed by long-term incubation. Moreover, buried subsoil SOC had higher apparent radiocarbon ages indicating that it is largely isolated from exchange with atmospheric CO2. We concluded that deep ploughing increased subsoil SOC storage and that the higher subsoil SOC stability is not only a result of selective preservation of more stable SOC fractions.

Jack-of-all-trades effects drive biodiversity-ecosystem multifunctionality relationships in European forests.

There is considerable evidence that biodiversity promotes multiple ecosystem functions (multifunctionality), thus ensuring the delivery of ecosystem services important for human well-being. However, the mechanisms underlying this relationship are poorly understood, especially in natural ecosystems. We develop a novel approach to partition biodiversity effects on multifunctionality into three mechanisms and apply this to European forest data. We show that throughout Europe, tree diversity is positively related with multifunctionality when moderate levels of functioning are required, but negatively when very high function levels are desired. For two well-known mechanisms, 'complementarity' and 'selection', we detect only minor effects on multifunctionality. Instead a third, so far overlooked mechanism, the 'jack-of-all-trades' effect, caused by the averaging of individual species effects on function, drives observed patterns. Simulations demonstrate that jack-of-all-trades effects occur whenever species effects on different functions are not perfectly correlated, meaning they may contribute to diversity-multifunctionality relationships in many of the world's ecosystems.

Influences of evergreen gymnosperm and deciduous angiosperm tree species on the functioning of temperate and boreal forests.

It has been recognized for a long time that the overstorey composition of a forest partly determines its biological and physical-chemical functioning. Here, we review evidence of the influence of evergreen gymnosperm (EG) tree species and deciduous angiosperm (DA) tree species on the water balance, physical-chemical soil properties and biogeochemical cycling of carbon and nutrients. We used scientific publications based on experimental designs where all species grew on the same parent material and initial soil, and were similar in stage of stand development, former land use and current management. We present the current state of the art, define knowledge gaps, and briefly discuss how selection of tree species can be used to mitigate pollution or enhance accumulation of stable organic carbon in the soil. The presence of EGs generally induces a lower rate of precipitation input into the soil than DAs, resulting in drier soil conditions and lower water discharge. Soil temperature is generally not different, or slightly lower, under an EG canopy compared to a DA canopy. Chemical properties, such as soil pH, can also be significantly modified by taxonomic groups of tree species. Biomass production is usually similar or lower in DA stands than in stands of EGs. Aboveground production of dead organic matter appears to be of the same order of magnitude between tree species groups growing on the same site. Some DAs induce more rapid decomposition of litter than EGs because of the chemical properties of their tissues, higher soil moisture and favourable conditions for earthworms. Forest floors consequently tend to be thicker in EG forests compared to DA forests. Many factors, such as litter lignin content, influence litter decomposition and it is difficult to identify specific litter-quality parameters that distinguish litter decomposition rates of EGs from DAs. Although it has been suggested that DAs can result in higher accumulation of soil carbon stocks, evidence from field studies does not show any obvious trend. Further research is required to clarify if accumulation of carbon in soils (i.e. forest floor + mineral soil) is different between the two types of trees. Production of belowground dead organic matter appears to be of similar magnitude in DA and EG forests, and root decomposition rate lower under EGs than DAs. However there are some discrepancies and still are insufficient data about belowground pools and processes that require further research. Relatively larger amounts of nutrients enter the soil-plant biogeochemical cycle under the influence of EGs than DAs, but recycling of nutrients appears to be slightly enhanced by DAs. Understanding the mechanisms underlying forest ecosystem functioning is essential to predicting the consequences of the expected tree species migration under global change. This knowledge can also be used as a mitigation tool regarding carbon sequestration or management of surface waters because the type of tree species affects forest growth, carbon, water and nutrient cycling.