# Above- and belowground biodiversity jointly tighten the P cycle in agricultural grasslands

Yvonne Oelmann, Markus Lange, Sophia Leimer, Christiane Roscher, Felipe Aburto, Fabian Alt, Nina Bange, Doreen Berner, Steffen Boch, Runa S Boeddinghaus, et al.

► To cite this version:

Yvonne Oelmann, Markus Lange, Sophia Leimer, Christiane Roscher, Felipe Aburto, et al. Above- and belowground biodiversity jointly tighten the P cycle in agricultural grasslands. Nature Communications, 2021, 12 (1), pp.4431. 10.1038/s41467-021-24714-4. hal-04177084

HAL Id: hal-04177084

https://hal.science/hal-04177084

Submitted on 3 Aug 2023

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

ARTICLE

https://doi.org/10.1038/s41467-021-24714-4

OPEN

<figure><img src="image_3.png" /></figure>

Check for updates

# Above- and belowground biodiversity jointly tighten the P cycle in agricultural grasslands

Yvonne Oelmann<sup>1,2,3</sup>, Markus Lange<sup>2</sup>, Sophia Leimer<sup>1,3</sup>, Christiane Roscher<sup>4,5</sup>, Felipe Aburto<sup>1,6</sup>, Fabian Alt<sup>1</sup>, Nina Bange<sup>1</sup>, Doreen Berner<sup>7</sup>, Steffen Boch<sup>1,8</sup>, Runa S. Boeddinghaus<sup>1,7</sup>, François Buscot<sup>1,9</sup>, Sigrid Dassen<sup>10</sup>, Gerlinde De Deyn<sup>1,11</sup>, Nico Eisenhauer<sup>1,5,12</sup>, Gerd Gleixner<sup>2</sup>, Kezia Goldmann<sup>1,9</sup>, Norbert Hölzel<sup>13</sup>, Malte Jochum<sup>1,5,12</sup>, Ellen Kandeler<sup>7</sup>, Valentin H. Klaus<sup>1,14</sup>, Till Kleinebecker<sup>1,15</sup>, Gaëtane Le Provost<sup>1,16</sup>, Peter Manning<sup>1,16</sup>, Sven Marhan<sup>7</sup>, Daniel Prati<sup>17</sup>, Deborah Schäfer<sup>17</sup>, Ingo Schöning<sup>2</sup>, Marion Schrumpf<sup>2</sup>, Elisabeth Schurig<sup>1</sup>, Cameron Wagg<sup>1,18,19</sup>, Tesfaye Wubet<sup>1,5,20</sup> & Wolfgang Wilcke<sup>1,3</sup>

Experiments showed that biodiversity increases grassland productivity and nutrient exploitation, potentially reducing fertiliser needs. Enhancing biodiversity could improve P-use efficiency of grasslands, which is beneficial given that rock-derived P fertilisers are expected to become scarce in the future. Here, we show in a biodiversity experiment that more diverse plant communities were able to exploit P resources more completely than less diverse ones. In the agricultural grasslands that we studied, management effects either overruled or modified the driving role of plant diversity observed in the biodiversity experiment. Nevertheless, we show that greater above- (plants) and belowground (mycorrhizal fungi) biodiversity contributed to tightening the P cycle in agricultural grasslands, as reduced management intensity and the associated increased biodiversity fostered the exploitation of P resources. Our results demonstrate that promoting a high above- and belowground biodiversity has ecological (biodiversity protection) and economical (fertiliser savings) benefits. Such win-win situations for farmers and biodiversity are crucial to convince farmers of the benefits of biodiversity and thus counteract global biodiversity loss.

<footnote n="1"><sup>1</sup>Geoecology, University of Tübingen, Tübingen, Germany. <sup>2</sup>Max Planck Institute for Biogeochemistry, Jena, Germany. <sup>3</sup>Institute of Geography and Geoecology, Karlsruhe Institute of Technology (KIT), Karlsruhe, Germany. <sup>4</sup>UFZ - Helmholtz Centre for Environmental Research, Physiological Diversity, Leipzig, Germany. <sup>5</sup>German Centre for Integrative Biodiversity Research (iDiv) Halle-Jena-Leipzig, Leipzig, Germany. <sup>6</sup>LISAB, Dep. Silvicultura, Universidad de Concepción, Concepción, Chile. <sup>7</sup>Institute of Soil Science and Land Evaluation, Soil Biology Department, University of Hohenheim, Stuttgart, Germany.</footnote>

<footnote n="8"><sup>8</sup>WSL Swiss Federal Research Institute, Birmensdorf, Switzerland. <sup>9</sup>UFZ - Helmholtz Centre for Environmental Research, Soil Ecology Department, Halle, Germany. <sup>10</sup>Department of Terrestrial Ecology, Netherlands Institute of Ecology, Wageningen, The Netherlands. <sup>11</sup>Department of Environmental Sciences, Soil Biology, University of Wageningen, Wageningen, The Netherlands. <sup>12</sup>Leipzig University, Institute of Biology, Leipzig, Germany. <sup>13</sup>Institute of Landscape Ecology, University of Münster, Münster, Germany. <sup>14</sup>Institute of Agricultural Sciences, ETH Zürich, Zürich, Switzerland. <sup>15</sup>Institute of Landscape Ecology and Resource Management, Justus-Liebig-University Gießen, Gießen, Germany. <sup>16</sup>Senckenberg Biodiversity and Climate Research Centre (SBiK-F), Frankfurt, Germany. <sup>17</sup>Institute of Plant Sciences, University of Bern, Bern, Switzerland. <sup>18</sup>Department of Evolutionary Ecology and Environmental Studies, University of Zürich, Zürich, Switzerland. <sup>19</sup>Fredericton Research and Development Center, Agriculture and Agri-Food Canada, Fredericton, NB, Canada.</footnote>

<footnote n="20"><sup>20</sup>UFZ - Helmholtz Centre for Environmental Research, Community Ecology Department, Halle, Germany. <sup>21</sup>email: yvonne.oelmann@uni-tuebingen.de</footnote>
Experiments have demonstrated that high biodiversity can increase ecosystem functioning and service provision¹⁻³. Increased productivity of plant species mixtures, in comparison to their component species in monocultures, is generally attributed to facilitation and niche complementarity of species associated with a more complete exploitation of resources⁴. Nutrients such as nitrogen (N) and phosphorus (P) are exploited more exhaustively with increasing diversity of plant mixtures⁵⁻⁹. Consequently, less N is leached to the groundwater¹⁰,¹¹. Therefore, high biodiversity contributes to a tight N cycle (contrary to a broken cycle associated with element loss¹²). However, P leaching was not reduced by increased plant species diversity⁸, partly because of the lower water solubility of mineral species of P compared to N, which forced plants to develop alternative P acquisition strategies beyond the direct uptake from soil solution by roots¹³. The low P concentrations in soil solution challenges the study of the P cycle. Accordingly, we lack a mechanistic understanding of biodiversity effects on the P cycle despite the importance of P as a limiting nutrient in terrestrial ecosystems in response to high N deposition¹⁴ and decreasing mineral P fertiliser availability, which has been identified as a threat to global food security¹⁵.

In contrast to N, P resource acquisition by nearly all plant species in the temperate climate zone is linked to a symbiosis with mycorrhizal fungi (arbuscular mycorrhizal fungi [AMF] in grasslands and ectomycorrhizal fungi [EMF] in forests)¹⁶. The plant-AMF symbiosis involves the direct transfer of P from the fungus to the plant. Species identity and the richness of plant and mycorrhizal fungal species influence productivity¹⁶,¹⁷ and thus, potentially the extent of P transfer. Furthermore, P acquisition of plants can also originate from an indirect transfer via the microbial community such as the release of P during turnover of the microbial P pool¹⁸ or from enzymatically catalyzed release from soil organic matter¹⁹. Since all three processes (P acquisition by plants via fungi, microbial and enzyme-mediated P release) were previously shown to increase with increasing plant species richness if studied separately⁸,¹⁹,²⁰, we hypothesise that effects of biodiversity—including plant and mycorrhizal species richness—on the P cycle are mediated by plant-microbe interactions.

While the mechanisms that underlie biodiversity effects on the P cycle in controlled grassland experiments are partially understood, it remains unknown whether the same mechanisms operate in permanent, managed agroecosystems such as agricultural grasslands. Approaching this knowledge gap from an experimental perspective, biodiversity and management can be manipulated simultaneously. Biodiversity experiments comprising management measures showed that positive biodiversity effects on productivity persisted on top of management effects²¹⁻²³. Furthermore, in such crossed biodiversity-fertilisation trials the exploitation of N resources by the plant community increased with increasing diversity of plant mixtures irrespective of fertilisation²⁴,²⁵. However, already during the short-term duration of these experiments, the authors mentioned that the seeded species composition (and richness) changed depending on the management treatments²³,²⁴. In long-established agricultural grasslands, biodiversity and management effects interact even more strongly²⁶. Management drives biodiversity because it can be considered as an environmental filter of the community composition²⁷. Grassland management directly or indirectly filters plant species composition because fertilisation, mowing, and grazing change abiotic site conditions. Accordingly, plant and AMF species richness are reduced by high management intensity associated with high resource availability in soil²⁸⁻³⁰. At the same time, the few plant species growing under high resource availability are known for their large productivity and thus, exploitative nutrient acquisition²⁸,³¹. Whether

resource exploitation in agricultural grasslands is dominated by biodiversity effects, i.e. the more diverse the more exhaustive is the resource exploitation, or by direct management effects, i.e. the more intensive and the more productive the more exhaustive is the resource exploitation, awaits to be shown. However, biodiversity and management are closely intertwingled, with high nutrient concentrations resulting in lower plant biodiversity in the long run, so that a fully crossed experiment cannot be maintained for very long in reality. We therefore suggest that the combination of insights gained from both, biodiversity experiments and observational data of agricultural grasslands represents an alternative promising avenue.

Here, we made use of the opportunity to combine results from 76 grassland mixtures of a biodiversity experiment with those from 100 permanent agricultural grasslands along a management gradient. All grasslands were located on soils developed from calcareous bedrock, which can be considered representative of large parts of Central Europe. We were particularly interested in the extent to which P resources are exploited by organisms (P exploitation, i.e. the part of the total bioavailable P pool in the ecosystem that is stored in the biomass) and whether P exploitation is related to biodiversity. We related organismic P stocks (in the plant and microbial community) to the total bioavailable P pool in the system (sum of both organismic P stocks and bioavailable P stocks in the soil [HPO₄²⁻ and H₂PO₄⁻ and organically bound P in labile and moderately labile P fractions accessible with moderate chemical extractants including NaHCO₃ and dilute NaOH])⁸. We first tested mechanisms underlying the biodiversity-P cycle relationship in the biodiversity experiment. Second, we explored whether biodiversity effects on the P cycle and the underlying mechanisms detected in the biodiversity experiment persisted along a management gradient (in terms of fertilisation, mowing, and grazing intensity). To do so, we used structural equation modelling (SEM) to set up an a priori defined mechanistic model of biodiversity effects on P exploitation in the biodiversity experiment and applied this model to data from the agricultural grasslands.

In this work, we show that in agricultural grasslands, management effects either overrule or modify the driving role of plant biodiversity observed in the biodiversity experiment. Nevertheless, greater above- (plants) and belowground (mycorrhizal fungi) biodiversity contributes to tightening the P cycle in agricultural grasslands. Our results demonstrate that promoting a high above- and belowground biodiversity has ecological (biodiversity protection) and economical (fertiliser savings) benefits.

## Results and discussion

**Explanation of the P exploitation.** Based on the knowledge that has been gathered so far for biodiversity effects on single processes of the P cycle (Table 1), we set up an SEM describing the relationships between biodiversity and P exploitation for the biodiversity experiment (Supplementary Fig. 1). We considered variables that are known to influence the extent to which P resources can be exploited including (i) the soil organic carbon (C<sub>org</sub>) stocks linked to microbial biomass³²,³³ and (ii) AMF species richness and relative abundance governing the symbiotic transfer of P from microbes to plants¹⁶,¹⁷. Because of the close correlation between AMF species richness and relative abundance, we introduced these two variables with correlated errors in the SEM.

We found that P exploitation was positively related to plant species richness and this relationship was mediated by C<sub>org</sub> stocks and microbial P stocks in soil (Fig. 1A). Under diverse plant mixtures, an increase of C<sub>org</sub> stocks in soil likely resulted from an increased C input from above- and belowground plant litter and
<table><caption>Table 1 Expected relationships underlying the structural equation models of the biodiversity experiment and the agricultural grasslands.</caption><thead><tr><th>Regression</th><th>Expected slope</th><th>Reason</th><th>Reference</th></tr></thead><tbody><tr><td>Plant species richness - AMF species richness</td><td>↑</td><td>Host specificity, AMF functional trait complementarity, mycelial networks</td><td>van der Heijden et al.<sup>17</sup>; Hiesalu et al.<sup>45</sup>; Dassen et al.<sup>73</sup></td></tr><tr><td>AMF species richness - Microbial P stocks</td><td>↑</td><td>More complete resource use achieved by diverse AMF community; inferred from AMF abundance - microbial P stocks relationship (see below) because abundance and species richness were closely related (Supplementary Table 1)</td><td>van der Heijden et al.<sup>79</sup></td></tr><tr><td>AMF species richness - Plant P stocks</td><td>↑</td><td>More complete resource use achieved by diverse AMF community and P subsequently provided to plant host</td><td>van der Heijden et al.<sup>17</sup>; Frew<sup>80</sup></td></tr><tr><td>Plant species richness - AMF abundance</td><td>↑</td><td>More niches for mycorrhizal colonization in diverse plant communities</td><td>van der Heijden et al.<sup>17</sup>; Hedlund et al.<sup>81</sup></td></tr><tr><td>AMF abundance - Microbial P stocks</td><td>↑</td><td>More complete resource use achieved by more abundant AMF community</td><td>van der Heijden et al.<sup>79</sup></td></tr><tr><td>AMF abundance - Plant P stocks</td><td>↑</td><td>More complete resource use achieved by more abundant AMF community and P subsequently provided to plant host</td><td>van der Heijden et al.<sup>17</sup>; Treseder et al.<sup>82</sup>; Köhl et al.<sup>83</sup></td></tr><tr><td>Plant species richness - C<sub>org</sub> stocks</td><td>↑</td><td>Increased above- and belowground plant litter input</td><td>Fornara &amp; Tilman<sup>84</sup>; Steinbeiss et al.<sup>85</sup>; Lange et al.<sup>34</sup></td></tr><tr><td>C<sub>org</sub> stocks - Microbial P stocks</td><td>↑</td><td>C<sub>org</sub> serves as an energy source and microbial homoeostasis requires a concomitant uptake of P</td><td>Hacker et al.<sup>86</sup>; Sorkau et al.<sup>32</sup></td></tr><tr><td>C<sub>org</sub> stocks - Plant P stocks</td><td>↑</td><td>C<sub>org</sub> stocks in soil are positively linked to plant biomass which forms part of plant P stock calculations</td><td>Fornara &amp; Tilman<sup>84</sup>; Steinbeiss et al.<sup>85</sup>; Lange et al.<sup>34</sup></td></tr><tr><td>Microbial P stocks - P exploitation</td><td>↑</td><td>Increased P exploitation if microbes take up and store more P</td><td></td></tr><tr><td>Plant P stocks - P exploitation</td><td>↑</td><td>Increased P exploitation if plants take up and store more P</td><td>Oelmann et al.<sup>81</sup></td></tr><tr><td>Management intensity - Plant species richness</td><td>↓</td><td>Environmental filtering: under high resource availability, dominant outcompete other species</td><td>Allan et al.<sup>28</sup></td></tr><tr><td>Management intensity - P exploitation</td><td>↑</td><td>Environmental filtering selects exploitative species</td><td>Allan et al.<sup>28</sup></td></tr></tbody></table>

↑ = positive slope of regression, ↓ = negative slope of regression, C<sub>org</sub> = organic carbon in soil.

root exudates³⁴. Given that C<sub>org</sub> serves as the energy source for the heterotrophic soil microbial community, C<sub>org</sub> concentrations were reported to be related to microbial C concentrations³⁴,³⁵. Consequently, plant species richness was also shown to be positively correlated with microbial C concentrations³⁴,³⁵. In turn, microbial C stocks were related to microbial P stocks, because microorganisms typically take up P proportional to their biomass to maintain nutrient homoeostasis³⁶. Our results may suggest nutrient homoeostasis of the microbial community since the ratios of microbial C and P (Cmic:Pmic) were not significantly related to plant species richness (Supplementary Table 1). Accordingly, our model indicates that a positive effect of plant species richness on P exploitation was related to the link between the P and the C cycle that allowed soil microorganisms to exploit P resources in soil more exhaustively.

Similarly, plant P stocks increased with increasing plant species richness (Fig. 1A). Plant P concentrations were not related to plant species richness while plant biomass was (Supplementary Table 1). In line with our mere bivariate correlation, increased biomass production of diverse grassland mixtures is well documented³⁷,³⁸. Biomass production is linked to the C cycle and thus, the simultaneous plant species richness effects on biomass production and plant P stocks reflect the coupling of plant species richness effects on the C and the P cycle. However, the biomass-induced increased P demand of the more diverse

plant communities was apparently not met by the AMF symbiosis. Although plant species richness and AMF species richness and relative abundance were related (Fig. 1A), this relationship did not directly feed back to the P cycle. There was neither a significant path connecting AMF species richness or relative abundance and plant P stocks (Fig. 1A) nor a significant relationship between plant P concentrations or plant biomass and AMF species richness or relative abundance (Supplementary Table 1). Instead, P release during microbial turnover might replenish the P pool in soil that is available for uptake by plants¹⁸. Our results from a recent study in the same biodiversity experiment corroborate the role of microbial turnover for the available P pool in soil²⁰.

In conclusion, our data suggest that plant species richness influenced the C cycle both above- and belowground which then translated into plant species richness effects on the P cycle. Ultimately, the more diverse the plant community was, the more P resources were exploited by plants and by soil microorganisms in the biodiversity experiment.

**Interactions of biodiversity and management.** In a first approach, we applied the SEM describing biodiversity effects on the P cycle in the biodiversity experiment to the observational data collected in the agricultural grasslands with the additional
<figure><img src="image_4.png" /><figcaption>**Fig. 1 Role of biodiversity for the P cycle.** Panel (A) refers to the biodiversity experiment and Panel (B) shows the agricultural grasslands. Biodiversity comprises plant and AMF (arbuscular mycorrhizal fungi) species richness. In the structural equation model (SEM), blue and red arrows display relationships with positive and negative slopes, respectively. Grey arrows indicate non-significant relationships. Solid lines of arrows refer to relationships expected according to the hypotheses (Table 1), whereas parallel lines of arrows indicate paths that were included in addition (selection procedure described in Methods). Round-shaped paths refer to the (driving) role of biodiversity while square paths indicate the role of management. Numbers on arrows indicate standardised path coefficients. Percentages in boxes give the explained variance. \*p < 0.05, \*\*p < 0.01, \*\*\*p < 0.001. The SEMs matched with the data: (A) Fisher's C = 7.12, p = 0.93, df = 14; (B) Fisher's C = 12.17, p = 0.43, df = 12. C<sub>org</sub> = organic carbon, P = phosphorus. Photo credits: The Jena Experiment (A), Jörg Hailer (B).</figcaption></figure>

consideration of management intensity (Supplementary Fig. 2). This model contained management intensity as an influential variable for plant species richness which then exerted the same influence as in the biodiversity experiment model. In addition, we inserted a direct relationship between management intensity and P exploitation. However, this initial SEM did not adequately represent the data from the agricultural grasslands (Supplementary Fig. 2).

We therefore included additional paths to optimise the model (see Methods and additional paths in Fig. 1). In the optimised model for the agricultural grasslands, four out of seven variables were directly related to management (Fig. 1B). The direct and positive influence of management intensity on plant P stocks was likely caused by P from fertilisers³⁹. Furthermore, management intensity contributed to the explanation of Corg stocks of agricultural grasslands (Fig. 1B). On the one hand, high management intensity commonly goes along with high livestock numbers and heavy machinery use and thus, increases soil bulk density⁴⁰. Because soil bulk density is used to calculate stocks (see Methods), increased soil bulk density translates into increased stocks. On the other hand, high management intensity can be associated with the application of organic fertilisers e.g., manure and slurry containing Corg and nutrients⁴¹. Compared to the biodiversity experiment, the effect of plant species richness was replaced by the influence of management on Corg stocks in the agricultural grasslands. In line with the considerations on the link between Corg stocks and microbial metabolism outlined above for the biodiversity experiment, a greater management intensity was also associated with increased microbial P stocks and thus, with an enlarged P exploitation (Fig. 1B). Similar to plants²⁸, P-exploitative microbial species seem to gain dominance under management-induced

increased P availability⁴². In general, future studies should focus on the relation of microbial taxa and their functioning in the P cycle⁴³ which could encompass the study of functional genes.

Low management intensity was related to high plant species richness (Fig. 1B). In addition, low management intensity was associated with low plant P stocks (Fig. 1B) via a decrease in plant biomass (Supplementary Table 2) because of decreased fertilizer input. These relationships are in line with well-documented beneficial effects of reduced management intensity on plant species richness²⁹. In addition, high plant species richness was linked to low plant P stocks (Fig. 1B) which likely is a result of a higher portion of plant species with low biomass production and low nutrient concentrations in these plant communities because of a conservative resource use²⁸. Accordingly, low management intensity was linked to plant P stocks directly via management measures and indirectly via filtering effects in agricultural grasslands that specifically remove plant species with a conservative resource use if land-use intensity increases. As a consequence, the relationship between plant species richness and plant P stocks in agricultural grasslands (Fig. 1B; 14 to 56 plant species per grassland plot) was reversed as compared to the positive, biomass-driven coefficient in the biodiversity experiment (Fig. 1A) covering species richness levels from monoculture to complete native communities (1-, 2-, 4-, 8-, 16- and 60-species mixtures). Furthermore, the significant relationship between plant P stocks and P exploitation in the biodiversity experiment disappeared in the agricultural grasslands (Fig. 1). Instead, there was a negative relationship between plant P stocks and microbial P stocks (Fig. 1). We suspect that microbes partly outcompeted plants in terms of P uptake⁴⁴ which might have been reinforced by less cooperative AMF species that preferentially retain P for
themselves³⁰,⁴⁵ and thus, increased microbial P stocks at the expense of plant P stocks. Yet, the positive relationship between plant species richness and P exploitation observed in the biodiversity experiment stayed the same in the agricultural grasslands: Increasing plant species richness increased P exploitation (Fig. 1). However, the mechanism underlying the plant species richness effect on P exploitation in the biodiversity experiment (i.e., the coupling of the C and P cycles above- and belowground) did not apply to the agricultural grasslands. Instead, the direct link between plant species richness and P exploitation implies that variables representing the underlying mechanism were not included in our model of the agricultural grasslands.

The less intense the management was, the more abundant and diverse were the AMF communities in the agricultural grasslands (Fig. 1B). Decreased P availability in soil due to reduced fertilisation is known to increase the species richness and abundance of the fungal community and in particular of the AMF community²⁹,³⁰. We were able to show that this increase in AMF abundance was associated with an increased P exploitation (Fig. 1B). The management effect on total P exploitation was mainly mediated via AMF and plant species richness (combined effect of all paths leading to P exploitation via AMF and plant species richness = -0.22), while the combined effect via C<sub>org</sub> stocks was negligible (0.002). Therefore, management effects in agricultural grasslands either overruled and modified the role of plant species richness observed in biodiversity experiments. Nonetheless, reduced management intensity was beneficial for above- and belowground biodiversity and for the exploitation of P resources in soil. Here we showed that both, above- and belowground biodiversity contribute to tightening the P cycle in agricultural grasslands.

**Biodiversity theory put into practice.** From an applied perspective, even small positive effects of biodiversity on productivity and P exploitation could yield benefits in terms of savings of fertilisation including labour and energy costs for application of the fertiliser. Such economic arguments should help to convince farmers to promote and conserve biodiversity. The balance between benefits of reduced management for biodiversity and potentially detrimental effects on yields remains a challenging issue and requires open-minded solutions. The identification of win-win situations is crucial to merge opposing views of representatives of agriculture – arguably one of the leading causes of biodiversity loss⁴⁶ – and nature conservation. Only if these two too frequent opponents act in concert, biodiversity theory will be put into practice and could help counteract global biodiversity loss and develop sustainable management strategies.

## Methods

**Study sites.** Grassland biodiversity experiment: Our data originate from the Jena Experiment (http://www.the-jena-experiment.de)⁴⁷. The field site is located near the German city of Jena (50°55' N, 11°35' E; 130 m above sea level). Mean annual air temperature is 9.9 °C, and mean annual precipitation amounts to 610 mm⁴⁸. The soil is an Eutric Fluvisol⁴⁹ developed from up to 2 m-thick fluvial sediments that are almost free of stones. Sediments largely originated from the same geological series (Pleistocene loess on limestone [Anisian/Ladinian] of the Middle Triassic) as described for one of the agricultural grassland regions (Hainich-Dün, see below). The systematic variation in soil texture as a consequence of fluvial dynamics was considered in the experimental design by arranging the experimental plots in four blocks at different distances to the river. The study site was converted from grassland to arable land in the early 1960s and used for cropping until the establishment of the experiment in 2002.

The main experiment comprised 82 plots that were split in different subplots with a core area of approximately 43.5 m². Each plot contains a specific combination of plant species (1, 2, 4, 8, 16, 60 species) that belong to different numbers (1, 2, 3, 4) of plant functional groups (grasses, [non-leguminous] small herbs, [non-leguminous] tall herbs, legumes). The species were chosen from a pool of 60 species typically found in mesic Molinio-Arrhenatheretea meadows. Each

plant species richness level had 16 replicates except for 14 mixtures with 16 species and 4 replicates of the 60-species mixture. Monocultures of *Bellis perennis* L. and *Cynosurus cristatus* L. included in the original design had to be abandoned because of poor establishment. We excluded the 60-species mixtures from statistical analyses (final *n* = 76), because they contain the complete species pool and thus, would potentially strengthen the sampling effect. To maintain the sown species richness levels, plots were weeded three times a year by carefully removing weeds and their roots. This minimises potential bias from weed root contributions. The management of all plots was adapted to extensive meadows used for hay production i.e., mown but not fertilised. All plots were mown twice a year in June and September with the harvest being removed from the plots.

**Agricultural grasslands.** We studied long-established, agriculturally managed (i.e. not experimentally assembled) grasslands in an interdisciplinary, large-scale and long-term programme, the Biodiversity Exploratories (http://www.biodiversity-exploratories.de)⁵⁰. The Biodiversity Exploratories comprise three study regions in Germany (Schwäbische Alb, Hainich-Dün, and Schorfheide-Chorin). We aimed to disentangle biodiversity and management effects and therefore, had to remove additional confounding factors such as geologic parent material (calcareous versus non-calcareous). Because of substantial differences in the geologic parent material⁵⁰ and of P fractions in soil⁵¹, we excluded the Schorfheide-Chorin from our analyses. Standardised field plots were located in the Schwäbische Alb and in the Hainich-Dün, both middle mountain ranges in Germany. The Schwäbische Alb is located between 460–860 m above sea level, has an annual mean temperature of 6–7 °C and an annual precipitation of 700–1,000 mm⁵⁰. The geologic parent material is calcareous bedrock of the Oxfordian Age (Epoch: Upper Jurassic) from which Leptosols and Cambisols have developed⁴⁹. The Hainich-Dün is located between 258–550 m above sea level, has an annual mean temperature of 6.5–8 °C and an annual precipitation of 500–800 mm⁵⁰. The geologic parent material in the Hainich-Dün is calcareous bedrock of the Middle Triassic Epoch partly covered by Pleistocene loess from which Cambisols, Stagnosols, Vertisols and Luvisols have developed⁴⁹. The design of the Biodiversity Exploratories was described in detail by Fischer et al.⁵⁰ We used 100 grassland plots (50 in each of Schwäbische Alb and Hainich-Dün) which can be classified as meadows (mown one to four times per year but not grazed), pastures (grazed but not mown), and mown pastures (both mown and grazed). Plots were selected to represent a gradient of management intensity with different fertilisation, frequency of mowing, and livestock units.

**Management intensity.** Management intensity of the biodiversity experiment has been the same since the establishment (no fertilisation, no grazing, mown twice a year). In the agricultural grasslands, management intensity (number of livestock units per hectare per year, frequency of mowing events per year, and amount of N-fertiliser applied per hectare per year) varied among plots and was assessed yearly via questionnaires answered by the farmers⁵². Notably, the plots were not amended with mineral P fertiliser except for two plots which received 25 kg ha⁻¹ P in 2011. However, manure and slurry were applied and also contain labile P fractions which might transform into sparingly soluble P minerals in the long run⁵³. To account for the variation in management of the agricultural grasslands, we used data of the sampling campaigns in 2011 and 2014 as described below. We first measured and calculated variables for each year individually. There were robust correlations between years for the different variables (correlation coefficient of up to *r* = 0.86) and therefore, we used means of the two years for further evaluation. We considered the mean land-use intensity of the years 2011 and 2014 as more robust than the single measurements alone, because the P availability in soil is determined by the site conditions and the long-term land use, particularly in the absence of mineral P fertilization as was the case at all our study sites, except two.

**Plant species richness and biomass harvest.** In the biodiversity experiment, sown species richness was determined on a 9 m² (3 m × 3 m) area in May 2014. In agricultural grasslands, vascular plant species richness was assessed on a 16 m² (4 m × 4 m) area by sampling all species from mid May to mid June in 2011 and 2014. In order to account for the difference in survey area between experimental and agricultural grasslands, a separate survey of 9 m² subplots nested within the 16 m² subplots was conducted in April and May 2018 in a selected number of plots in the agricultural grasslands (*n* = 18). Species richness differed significantly between 9 m² and 16 m² subplots (mean number of species ± standard error; 9 m²; 22 ± 2; 16 m²; 24 ± 2; paired t-test: T = -5.15; p < 0.001). Accordingly, we used a scaling factor of 0.91 i.e., the slope of the regression of species richness on 9 m² on species richness on 16 m² (Supplementary Fig. 3), to adjust species richness of all agricultural grasslands to the area of species richness measurements in the biodiversity experiment.

In the biodiversity experiment, aboveground plant biomass of target species was harvested in May 2014, at the estimated peak biomass before mowing. Aboveground plant biomass was sampled in all plots within a frame (0.2 × 0.5 m, height 0.03 m) at two random locations per plot and sorted by target species, weeds and detached dead plant material. Biomass was extrapolated to 1 m². In agricultural grasslands, biomass of a 2 m² large area was harvested as mixed samples of eight quadrats of 0.25 m² in close proximity to the area in which the vegetation was recorded from mid May to mid June in 2011 and 2014. Temporary fences ensured that biomass was sampled at peak standing crop but without any effect of possible
mowing or grazing events. Detached dead material was excluded from biomass sampling. Biomass was assessed on a dry-weight basis (drying at 70–80 °C for 48 h) in the experimental and agricultural grasslands.

**Soil sampling.** Soil in the biodiversity experiment was sampled in April 2014 for analyses of C<sub>org</sub> concentrations and bulk density. In each plot, three soil cores were taken to a depth of 0.3 m using a split-tube sampler (4.8 cm diameter). Soil cores were segmented into 5 cm-depth sections and pooled per depth sections and plot<sup>34</sup>. We used the mean C<sub>org</sub> concentrations of the upper three intervals (0 to 0.15 m). For the measurement of P in soil of the biodiversity experiment, sampling took place in September 2013. Nine soil cores per plot with a diameter of 0.02 m were taken at a depth of 0 to 0.15 m and combined to a composite sample considered representative for the plot. Although these samples were not taken during the identical period of time as for the agricultural grasslands (see below), we tested whether the results are applicable to the growing season of 2014. First, bioavailable Pi concentrations in soil were closely correlated between years (September 2013 and October 2014; *r* = 0.88, *p* < 0.001, *n* = 79). Second, resin-extractable P concentrations on samples of the 2013 campaign measured in the laboratory matched with P released under field conditions in May 2014 (*r* = 0.40, *p* < 0.001, *n* = 78).

In agricultural grasslands, soil samples were collected simultaneously within two weeks in May 2011 and May 2014. Samples were taken along two orthogonal transects of 20 m. Sampling points were shifted by 0.5 m in 2014 compared with 2011 to avoid an overlap of sampling positions. In each plot, 14 samples from 0 to 0.1 m soil depth were taken using core augers (diameter ~52 mm). Samples were mixed, cooled and transported to a field lab where they were sieved (<2 mm), all within 8 h of sampling. All measurements described below for agricultural grasslands were done on aliquots of samples of these joint sampling campaigns.

**Laboratory analyses and calculations.** Air-dry plant material sampled in the experimental and agricultural grasslands was ground with a mill using a 0.5-mm screen for chemical analyses. Nitrogen concentration of plant material was measured in ground samples with an elemental analyser for the biodiversity experiment and near-infrared spectroscopy (NIRS) for the agricultural grasslands. Plant samples were digested in a microwave with concentrated nitric acid and hydrogen peroxide<sup>54</sup> and P concentrations were determined by inductively-coupled plasma optical emission spectrometry (ICP-OES). In 2011, P concentrations were measured by NIRS in biomass of the agricultural grasslands. We multiplied the latter data by 0.9844 to match with the digestion method based on data of the year 2009 for which both the digestion and NIRS method had been applied (P<sub>digestion</sub> = 0.9844 × P<sub>NIRS</sub>, *r* = 0.80, *p* < 0.001, *n* = 98). Plant samples collected in 2014 in the agricultural grasslands were analysed for P concentrations by means of an X-ray fluorescence spectrometer. We can exclude a methodological shift in plant P concentrations between years, because digestion followed by ICP-OES analyses and XRF analyses were shown to match well<sup>55</sup>. Nitrogen:P ratios in plant material were calculated on an elemental mass basis. We calculated P stocks in aboveground biomass by multiplying biomass [g m<sup>-2</sup>] harvested in May with their P concentrations [mg g<sup>-1</sup>] and with the number of mowing events. In this way, we likely overestimated the absolute values of above-ground plant P stocks. But at the same time, we ensured comparability between experimental and agricultural grasslands because, for the latter, plant material of mowing events later than May was not available.

In the experimental and agricultural grasslands, soil pH values were determined with a glass electrode in a 1:2.5 soil:0.01 M CaCl<sub>2</sub> water suspension. Total C concentrations were determined on ground air-dry samples by an elemental analyser. Organic C concentrations were calculated by subtracting inorganic C concentrations (determined after removal of C<sub>org</sub> at 450 °C in a muffle furnace) from total C concentration. Soil bulk density was calculated by weighing a 100-cm³ core after drying the soil at 40 °C.

In the experimental and agricultural grasslands, various P fractions in soil were measured according to the method of Hedley et al.<sup>56</sup> modified by Kuo<sup>57</sup>. The sequential extraction scheme had three steps (bioavailable P (NaHCO<sub>3</sub> extractable), moderately labile P (NaOH extractable), and mineral P (HCl extractable)). The bioavailable P fraction comprises H<sub>2</sub>PO<sub>4</sub><sup>-</sup>/HPO<sub>4</sub><sup>2-</sup> ions in soil solution and those weakly adsorbed to mineral surfaces. More strongly adsorbed P to iron and aluminium oxides and (oxy)hydroxides is recovered in the moderately labile fraction. The mineral P fraction contains P bound in apatite and other Ca phosphates<sup>56,58,59</sup>. In all extraction solutions, Pi concentrations were analysed using the ammonium molybdate-ascorbic acid blue method described by Murphy and Riley<sup>60</sup> and measured with a continuous flow analyser. Total dissolved P concentrations in NaHCO<sub>3</sub>- and NaOH-extracts were measured with an ICP-OES. For the labile and moderately labile fractions (NaHCO<sub>3</sub>-P, NaOH-P), organic P concentrations were calculated by subtracting Pi from total dissolved P concentrations in the extracts.

In the experimental and agricultural grasslands, microbial P was measured by hexanol fumigation using a combination of the methods of McLaughlin et al.<sup>61</sup> and Kouno et al.<sup>62</sup>. Three subsamples of each soil sample were prepared by adding deionised water and one anion-exchange membrane. One of the three subsamples was additionally mixed with hexanol as fumigation reagent (P<sub>Hex</sub>) and one with a P spike (P<sub>Spike</sub>), while nothing was added to the last subsample (P<sub>H2O</sub>). Nitric acid was used to exchange the P adsorbed onto the membranes. Phosphate concentrations in solutions were measured with a continuous flow

analyser. We calculated microbial P concentrations as the difference between hexanol-fumigated and non-fumigated samples and accounted for P retention during extraction by including the P spike<sup>63</sup>. Because the calculated microbial P concentrations underestimate the amount of P stored in microbial biomass because of fumigation/extraction efficiency constraints e.g., with respect to gram-positive bacteria, we divided calculated microbial P concentrations by 0.462,<sup>64</sup>. We calculated stocks of P fractions and microbial P in soil based on bulk density [g m<sup>-2</sup> (0.15 m soil depth)<sup>-1</sup>] and multiplication with the respective P concentrations [mg g<sup>-1</sup>]. Phosphorus exploitation was calculated as the contribution [%] of organismic P stocks (either in the aboveground part of plants or microbes) to the sum of bioavailable P stocks (both organismic P stocks + labile P stocks + moderately labile P stocks) according to Oelmann et al.<sup>8</sup>. The aboveground plant P stocks represent the plant demand of P that is removed with the harvest(s) each year. It is reasonable to assume that microbial P stocks in soil can also be regarded as the annual microbial P demand since it has been shown that the microbial P stock in soil turns over once every growing season irrespective of management<sup>39</sup>.

In the biodiversity experiment, phosphomonoesterase (PASE) activity was measured according to the modified assay of Eivazi and Tabatabai<sup>65</sup>. For each soil sample, one replicate and one blank value were included in the laboratory study. We incubated soil samples with *p*-nitrophenylphosphate as an organic substrate for enzyme activity (pH = 11). *p*-nitrophenylphosphate was added to blanks only after incubation. Directly after filtration, *p*-nitrophenol concentrations were measured with a spectrophotometer. In agricultural grasslands, PASE activity was determined by fluorescence measures in a buffered solution of pH 6.1 after Marx et al.<sup>66</sup>, as described in Berner et al.<sup>67</sup>.

In the biodiversity experiment, the fungal-to-bacterial ratio was assessed by applying the phospholipid fatty acids (PLFA) method<sup>68</sup>. Within 48 h after sampling, the soil was kept at 4 °C, sieved to 2 mm, remains of roots were manually removed and the samples were stored at -20 °C until further sample processing. PLFA were extracted according to the method of Bligh and Dyer<sup>69</sup> as modified by Kramer and Gleixner<sup>70</sup>. As an indicator for fungal PLFA 18:2ω6,9 was used<sup>71,72</sup>. The bacterial PLFA was calculated as sum of the PLFA markers 14:0i, 15:0i, 15:0a, 16:0i, c16:1ω7c, 17:1, 17:0i, 17:0a and 18:1ω7 (Frostegard & Baath 1996). Furthermore, AMF species richness and relative abundance in the biodiversity experiment was investigated in soil samples (0-15 cm depth) collected in 2010 and analysed using DNA extraction and amplicon sequencing by 454-pyrosequencing as described in Dassen et al.<sup>73</sup>. We amplified 18 S rRNA gene fragments from fungi and protists with primer FR 1 and the modified version of FF390 designed to also include the Glomeromycota comprising AMF<sup>74</sup> (Supplementary Table 3). In agricultural grasslands, AMF were identified based on DNA extraction. DNA was extracted from soil of each plot using the MO BIO Power Soil DNA isolation kit (MO BIO Laboratories, Carlsbad, CA, USA) following the manufacturer's protocol. Afterwards we used a nested PCR approach to amplify fungal 18S-rDNA by using the primer pairs GlomerWT0/ Glomer1536<sup>75</sup> and NS31/AML2<sup>76,77</sup> (Supplementary Table 3), containing the Illumina adapter sequences. PCR products were then purified, cleaned and sequenced using Illumina MiSeq. The AMF sequences were processed using a customized bioinformatic pipeline following MOTHUR SOP (using Version 1.39.5) as implemented in DeltaMP (https://github.com/lentendu/DeltaMP). The taxonomical assignment was done against the Maarjam database (https://maarjam.botany.uee/). AMF OTUs were merged according to VT assignment. OTUs assigned only to genus level were kept as OTUs. Only those AMF appearing on more than five plots were considered. AMF species richness was calculated as the number of species including OTUs. The relative abundance of AMF was calculated by relating the reads per species/OTU to the total sum of reads across all plots. AMF appearing on less than five plots had low relative abundances (<1%) and thus, were considered to play a negligible role. Accordingly, only those AMF appearing on more than five plots were included. AMF species richness was calculated as the number of species including OTUs. The relative abundance of AMF was calculated by relating the reads per species/OTU to the total sum of reads across all plots.

**Statistical analyses.** Normal distribution of the variables and the homogeneity of the models were visually inspected and if necessary, variables were log-transformed to meet the prerequisite for statistical analyses (Supplementary Tables 1, 2). Pearson correlations were calculated applying the corr.test function in the R library psych<sup>66</sup>. The Pearson correlations were corrected for multiple inference using the Benjamini & Hochberg correction. Using the R library piecewiseSEM<sup>78</sup>, confirmatory path analyses were applied to test the causal relationships between plant diversity and P exploitation in the biodiversity experiment and in the agricultural grasslands and how land-use intensity impacts this relation in the latter. To identify the processes underlying the plant diversity and the land use intensity effects on P exploitation, we set up conceptual path models for each, the biodiversity experiment and the agricultural grasslands (Supplementary Figs. 1 and 2; Table 1). The models include AMF species richness, relative AMF abundance, C<sub>org</sub> stocks, plant P stocks and microbial P stocks as possible mediators. To account for the block design of the biodiversity experiment and the two regions considered in the agricultural grasslands, block was fitted as a random factor in the biodiversity experiment models and region as a random effect in the agricultural grassland models. Furthermore, we checked for spatial autocorrelation within each region. To this end, we related the residuals associated to each variable in the SEM
to the geographic coordinates of each plot. Since we found no indications of relationships or clustering of the residuals, we assumed independent data within the regions of the agricultural grasslands. In order to test whether the conceptual models adequately represent the measured data Fisher's C test statistic was used. To obtain adequate models the conceptual models were updated stepwise by including missing paths as indicated by the tests of directed separation⁷. Details of the outcome of the SEMs in addition to that provided in Fig. 1 can be found in Supplementary Tables 4, 5 and 6.

**Reporting summary.** Further information on research design is available in the Nature Research Reporting Summary linked to this article.

## Data availability

This work is based on data from several projects of the Biodiversity Exploratories programme (DFG Priority Program 1374) and The Jena Experiment (DFG FOR 456 & 1451). The data used for analyses are publicly available from the Biodiversity Exploratories Information System (https://doi.org/10.17616/R32P9Q) and the jexis database (https://jexis.idiv.de/), respectively, or will become publicly available after an embargo period of five years from the end of data assembly to give the owners and collectors of the data time to perform their analysis. The raw AMF Illumina sequences for the Biodiversity Exploratories have been deposited in the National Center for Biotechnology Information (NCBI) Sequence Read Archive (SRA) under BioProject accession number PRJNA706003. Microbial community data of the Jena Experiment have been archived in the Pangaea database (https://doi.org/10.1594/pangaea.874990). Any other relevant data are available from the corresponding author upon reasonable request.

Received: 29 June 2020; Accepted: 25 June 2021;
Published online: 21 July 2021

## References

1. Cardinale, B. J. et al. Biodiversity loss and its impact on humanity. *Nature* **486**, 59–67 (2012).

2. Hooper, D. U. et al. Effects of biodiversity on ecosystem functioning: a consensus of current knowledge. *Ecol. Monogr.* **75**, 3–35 (2005).

3. Isbell, F. et al. Linking the influence and dependence of people on biodiversity across scales. *Nature* **546**, 65–72 (2017).

4. Loreau, M. & Hector, A. Partitioning selection and complementarity in biodiversity experiments. *Nature* **412**, 72–76 (2001).

5. Bessler, H. et al. Nitrogen uptake by grassland communities: contribution of N₂ fixation, facilitation, complementarity, and species dominance. *Plant Soil* **358**, 301–322 (2012).

6. Karanika, E. D., Alifragis, D. A., Mamolos, A. P. & Veresoglou, D. S. Differentiation between responses of primary productivity and phosphorus exploitation to species richness. *Plant Soil* **297**, 69–81 (2007).

7. Lange, M. et al. How plant diversity impacts the coupled water, nutrient and carbon cycles. *Adv. Ecol. Res.* **61**, 185–219 (2019).

8. Oelmann, Y. et al. Does plant diversity influence phosphorus cycling in experimental grasslands? *Geoderma* **167-68**, 178–187 (2011).

9. Tilman, D., Wedin, D. & Knops, J. Productivity and sustainability influenced by biodiversity in grassland ecosystems. *Nature* **379**, 718–720 (1996).

10. Leimer, S., Oelmann, Y., Wirth, C. & Wilcke, W. Time matters for plant diversity effects on nitrate leaching from temperate grassland. *Agric Ecosyst. Environ.* **211**, 155–163 (2015).

11. Scherer-Lorenzen, M., Palmborg, C., Prinz, A. & Schulze, E.-D. The role of plant diversity and composition for nitrate leaching in grasslands. *Ecology* **84**, 1539–1552 (2003).

12. Elser, J. & Bennett, E. A broken biogeochemical cycle. *Nature* **478**, 29–31 (2011).

13. Lambers, H., Mougel, C., Jaillard, B. & Hinsinger, P. Plant-microbe-soil interactions in the rhizosphere: an evolutionary perspective. *Plant Soil* **321**, 83–115 (2009).

14. Wassen, M. J., Olde Venterink, H., Lapshina, E. D. & Tanneberger, F. Endangered plants persist under phosphorus limitation. *Nature* **437**, 547–550 (2005).

15. Cordell, D., Drangert, J.-O. & White, S. The story of phosphorus: Global food security and food for thought. *Glob. Environ. Change-Hum. Policy Dimens.* **19**, 292–305 (2009).

16. van der Heijden, M. G. A., Martin, F. M., Selosse, M.-A. & Sanders, I. R. Mycorrhizal ecology and evolution: the past, the present, and the future. *N. Phytol.* **205**, 1406–1423 (2015).

17. van der Heijden, M. G. A. et al. Mycorrhizal fungal diversity determines plant biodiversity, ecosystem variability and productivity. *Nature* **396**, 69–72 (1998).

18. Richardson, A. E. & Simpson, R. J. Soil microorganisms mediating phosphorus availability. *Plant Physiol.* **156**, 989–996 (2011).

19. Hacker, N. et al. Plant diversity shapes microbe-rhizosphere effects on P mobilisation from organic matter in soil. *Ecol. Lett.* **18**, 1356–1365 (2015).

20. Hacker, N., Wilcke, W. & Oelmann, Y. The oxygen isotope composition of bioavailable phosphate in soil reflects the oxygen isotope composition in soil water driven by plant diversity effects on evaporation. *Geochim. Cosmochim. Acta* **248**, 387–399 (2019).

21. Craven, D. et al. Plant diversity effects on grassland productivity are robust to both nutrient enrichment and drought. *Philos. Trans. R. Soc. B* **371**, 8 (2016).

22. Fridley, J. D. Resource availability dominates and alters the relationship between species diversity and ecosystem productivity in experimental plant communities. *Oecologia* **132**, 271–277 (2002).

23. Weigelt, A., Weisser, W. W., Buchmann, N. & Scherer-Lorenzen, M. Biodiversity for multifunctional grasslands: equal productivity in high-diversity low-input and low-diversity high-input systems. *Biogeosciences* **6**, 1695–1706 (2009).

24. Nyfeler, D. et al. Strong mixture effects among four species in fertilized agricultural grassland led to persistent and consistent transgressive overyielding. *J. Appl Ecol.* **46**, 683–691 (2009).

25. Oelmann, Y., Vogel, A., Wegener, F., Weigelt, A. & Scherer-Lorenzen, M. Management intensity modifies plant diversity effects on N yield and mineral N in soil. *Soil Sci. Soc. Am. J.* **79**, 559–568 (2015).

26. Manning P., et al. Transferring biodiversity-ecosystem function research to the management of 'real-world' ecosystems. In: *Mechanisms Underlying the Relationship between Biodiversity and Ecosystem Function* (eds Eisenhauer N., Bohan D. A., Dumbrell A. J.). Academic Press Ltd-Elsevier Science Ltd (2019).

27. Kraft, N. J. B. et al. Community assembly, coexistence and the environmental filtering metaphor. *Funct. Ecol.* **29**, 592–599 (2015).

28. Allan, E. et al. Land use intensification alters ecosystem multifunctionality via loss of biodiversity and changes to functional composition. *Ecol. Lett.* **18**, 834–843 (2015).

29. Collins, C. D. & Foster, B. L. Community-level consequences of mycorrhizae depend on phosphorus availability. *Ecology* **90**, 2567–2576 (2009).

30. Klironomos, J. N., McCune, J., Hart, M. & Neville, J. The influence of arbuscular mycorrhizae on the relationship between plant diversity and productivity. *Ecol. Lett.* **3**, 137–141 (2000).

31. Busch, V. et al. Will I stay or will I go? Plant species-specific response and tolerance to high land-use intensity in temperate grassland ecosystems. *J. Veg. Sci.* **30**, 674–686 (2019).

32. Sorkau, E. et al. The role of soil chemical properties, land use and plant diversity for microbial phosphorus in forest and grassland soils. *J. Plant Nutr. Soil Sci.* **181**, 185–197 (2018).

33. Wardele, D. A. A comparative assessment of factors which influence microbial biomass carbon and nitrogen levels in soil. *Biol. Rev. Camb. Philos. Soc.* **67**, 321–358 (1992).

34. Lange, M. et al. Plant diversity increases soil microbial activity and soil carbon storage. *Nat. Commun.* **6**, 6707 (2015).

35. Eisenhauer, N. et al. Plant diversity effects on soil microorganisms support the singular hypothesis. *Ecology* **91**, 485–496 (2010).

36. Cleveland, C. C. & Liptzin, D. C. N: P stoichiometry in soil: is there a "Redfield ratio" for the microbial biomass? *Biogeochemistry* **85**, 235–252 (2007).

37. Cardinale, B. J. et al. Impacts of plant diversity on biomass production increase through time because of species complementarity. *Proc. Natl Acad. Sci. USA* **104**, 18123–18128 (2007).

38. Marquard, E. et al. Plant species richness and functional composition drive overyielding in a 6-year grassland experiment. *Ecology* **90**, 3290–3302 (2009).

39. Liebisch, F. et al. Seasonal dynamics and turnover of microbial phosphorus in a permanent grassland. *Biol. Fertil. Soils* **50**, 465–475 (2014).

40. Boeddinghaus, R. S. et al. Plant functional trait shifts explain concurrent changes in the structure and function of grassland soil microbial communities. *J. Ecol.* **107**, 2197–2210 (2019).

41. Soussana, J. F. et al. Carbon cycling and sequestration opportunities in temperate grasslands. *Soil Use Manag.* **20**, 219–230 (2004).

42. Waldrop, M. P., Zak, D. R., Blackwood, C. B., Curtis, C. D. & Tilman, D. Resource availability controls fungal diversity across a plant diversity gradient. *Ecol. Lett.* **9**, 1127–1135 (2006).

43. Kour, D. et al. Biodiversity, current developments and potential biotechnological applications of phosphorus-solubilizing and -mobilizing microbes: a review. *Pedosphere* **31**, 43–75 (2021).

44. Dijkstra, F. A., He, M. Z., Johansen, M. P., Harrison, J. J. & Keitel, C. Plant and microbial uptake of nitrogen and phosphorus affected by drought using N-15 and P-32 tracers. *Soil Biol. Biochem.* **82**, 135–142 (2015).

45. Hiiesalu, I. et al. Species richness of arbuscular mycorrhizal fungi: associations with grassland plant richness and biomass. *N. Phytol.* **203**, 233–244 (2014).

46. Tilman, D., Cassman, K. G., Matson, P. A., Naylor, R. & Polasky, S. Agricultural sustainability and intensive production practices. *Nature* **418**, 671–677 (2002).

47. Roscher, C. et al. The role of biodiversity for element cycling and trophic interactions: an experimental approach in a grassland community. *Bas Appl. Ecol.* **5**, 107–121 (2004).
48. Hoffmann K., Bivour W., Früh B., Koßmann M., Voß P.-H. *Climate studies in Jena for adaptation to climate change and is expected consequences.* (In *German*). Selbstverlag des Deutschen Wetterdienstes (2014).

49. IUSS Working Group WRB. *World Reference Base for Soil Resources 2014, update 2015: International soil classification system for naming soils and creating legends for soil maps*. FAO (2015).

50. Fischer, M. et al. Implementing large-scale and long-term functional biodiversity research: the biodiversity exploratories. *Bas Appl Ecol.* **11**, 473–485 (2010).

51. Alt, F., Oelmann, Y., Herold, N., Schrumpf, M. & Wilcke, W. Phosphorus partitioning in grassland and forest soils of Germany as related to land-use type, management intensity, and land use-related pH. *J. Plant Nutr. Soil Sci.* **174**, 195–209 (2011).

52. Vogt, J. et al. Eleven years' data of grassland management in Germany. *Biodiver Data J.* **7**, 38 (2019).

53. Alt, F., Oelmann, Y., Schöning, I. & Wilcke, W. Phosphate release kinetics at stable pH in calcareous grassland and forest soils. *Soil Sci. Soc. Am. J.* **77**, 2060–2070 (2013).

54. Jones J. B., Wolf B., Mills H. A. *Plant analysis handbook*. Micro Macro Publishing (1991).

55. Marina, M. A. & Lopez, M. C. B. Determination of phosphorus in raw materials for ceramics: comparison between X-ray fluorescence spectrometry and inductively coupled plasma-atomic emission spectrometry. *Anal. Chim. Acta* **432**, 157–163 (2001).

56. Hedley, M. J., Stewart, J. W. B. & Chauhan, B. S. Changes in inorganic and organic soil-phosphorus fractions induced by cultivation practices and by laboratory incubations. *Soil Sci. Soc. Am. J.* **46**, 970–976 (1982).

57. Kuo S. Phosphorus. In: *Methods of Soil Analysis - Part 3 Chemical Methods* (eds Sparks D. L., et al.). SSSA (1996).

58. Cross, A. F. & Schlesinger, W. H. A literature review and evaluation of the Hedley fractionation - applications to the biogeochemical cycle of soil phosphorus in natural ecosystems. *Geoderma* **64**, 197–214 (1995).

59. Negassa, W. & Leinweber, P. How does the Hedley sequential phosphorus fractionation reflect impacts of land use and management on soil phosphorus: a review. *J. Plant Nutr. Soil Sci.* **172**, 305–325 (2009).

60. Murphy, J. & Riley, J. P. A modified single solution method for determination of phosphate in natural waters. *Anal. Chim. Acta* **26**, 31–36 (1962).

61. McLaughlin, M. J., Alston, A. M. & Martin, J. K. Measurement of phosphorus in the soil microbial biomass - a modified procedure for field soils. *Soil Biol. Biochem.* **18**, 437–443 (1986).

62. Kouno, K., Tuchiya, Y. & Ando, T. Measurement of soil microbial biomass phosphorus by an anion exchange membrane method. *Soil Biol. Biochem.* **27**, 1353–1357 (1995).

63. Bünemann, E. K., Marschner, P., Smernik, R. J., Conyers, M. & McNeill, A. M. Soil organic phosphorus and microbial community composition as affected by 26 years of different management strategies. *Biol. Fertil. Soils* **44**, 717–726 (2008).

64. Brookes, P. C., Powlson, D. S. & Jenkinson, D. S. Measurement of microbial biomass phosphorus in soil. *Soil Biol. Biochem* **14**, 319–329 (1982).

65. Eivazi, F. & Tabatabai, M. A. Phosphatases in soils. *Soil Biol. Biochem.* **9**, 167–172 (1977).

66. Marx, M. C., Wood, M. & Jarvis, S. C. A microplate fluorimetric assay for the study of enzyme diversity in soils. *Soil Biol. Biochem.* **33**, 1633–1640 (2001).

67. Berner, D. et al. Land-use intensity modifies spatial distribution and function of soil microorganisms in grasslands. *Pedobiologia* **54**, 341–351 (2011).

68. White, D. C., Davis, W. M., Nickels, J. S., King, J. D. & Bobbie, R. J. Determination of the sedimentary microbial biomass by extractable lipid phosphate. *Oecologia* **40**, 51–62 (1979).

69. Bligh, E. G. & Dyer, W. J. A rapid method of total lipid extraction and purification. *Can. J. Biochem. Physiol.* **37**, 911–917 (1959).

70. Kramer, C. & Gleixner, G. Variable use of plant- and soil-derived carbon by microorganisms in agricultural soils. *Soil Biol. Biochem.* **38**, 3267–3278 (2006).

71. Frostegard, A. & Baath, E. The use of phospholipid fatty acid analysis to estimate bacterial and fungal biomass in soil. *Biol. Fertil. Soils* **22**, 59–65 (1996).

72. Zelles, L. Identification of single cultured micro-organisms based on their whole-community fatty acid profiles, using an extended extraction procedure. *Chemosphere* **39**, 665–682 (1999).

73. Dassen, S. et al. Differential responses of soil bacteria, fungi, archaea and protists to plant species richness and plant functional group identity. *Mol. Ecol.* **26**, 4085–4098 (2017).

74. Kuramae, E. E. et al. Tracking fungal community responses to maize plants by DNA- and RNA-based pyrosequencing. *PLoS ONE* **8**, 8 (2013).

75. Wubet, T., Weiss, M., Kottke, I. & Oberwinkler, F. Two threatened coexisting indigenous conifer species in the dry Afromontane forests of Ethiopia are associated with distinct arbuscular mycorrhizal fungal communities. *Can. J. Bot.-Rev. Canadienne De Botanique* **84**, 1617–1627 (2006).

76. Lee, J., Lee, S. & Young, J. P. W. Improved PCR primers for the detection and identification of arbuscular mycorrhizal fungi. *FEMS Microbiol. Ecol.* **65**, 339–349 (2008).

77. Simon, L., Lalonde, M. & Bruns, T. D. Specific amplification of 18S fungal ribosomal genes from vesicular-arbuscular endomycorrhizal fungi colonizing roots. *Appl. Environ. Microbiol.* **58**, 291–295 (1992).

78. Lefcheck, J. S. PIECEWISESEM: Piecewise structural equation modelling in R for ecology, evolution, and systematics. *Methods Ecol. Evol.* **7**, 573–579 (2016).

79. van der Heijden, M. G. A. et al. The mycorrhizal contribution to plant productivity, plant nutrition and soil structure in experimental grassland. *N. Phytol.* **172**, 739–752 (2006).

80. Frew, A. Arbuscular mycorrhizal fungal diversity increases growth and phosphorus uptake in C-3 and C-4 crop plants. *Soil Biol. Biochem.* **135**, 248–250 (2019).

81. Hedlund, K. et al. Plant species diversity, plant biomass and responses of the soil community on abandoned land across Europe: idiosyncracy or above- belowground time lags. *Oikos* **103**, 45–58 (2003).

82. Treseder, K. K. The extent of mycorrhizal colonization of roots and its influence on plant growth and phosphorus content. *Plant Soil* **371**, 1–13 (2013).

83. Köhl, L., Oehl, F. & van der Heijden, M. G. A. Agricultural practices indirectly influence plant productivity and ecosystem services through effects on soil biota. *Ecol. Appl.* **24**, 1842–1853 (2014).

84. Fornara, D. A. & Tilman, D. Plant functional composition influences rates of soil carbon and nitrogen accumulation. *J. Ecol.* **96**, 314–322 (2008).

85. Steinbeiss, S. et al. Plant diversity positively affects short-term soil carbon storage in experimental grasslands. *Glob. Change Biol.* **14**, 2937–2949 (2008).

86. Hacker N. Phosphorus Release Mechanisms in an Experimental Grassland of Varying Biodiversity. Doctoral thesis, University of Tübingen, Germany (2017).

## Acknowledgements

We thank all people who helped with the setup and management of the Jena Experiment (biodiversity experiment) and in particular the initiators, E.-D. Schulze, B. Schmid, and W.W. Weisser. Thanks also to all the helpers who assisted during the weeding campaigns. The Jena Experiment is funded by the Deutsche Forschungsgemeinschaft (DFG, FOR 456 & 1451) with additional support from the Swiss National Science Foundation (SNF), the Friedrich Schiller University Jena and the Max Planck Society. NE acknowledges support by the German Centre for Integrative Biodiversity Research Halle-Jena-Leipzig, funded by the German Research Foundation (FZT 118). M.L. gratefully acknowledges the support of the Zwillenberg-Tietz Foundation. We further thank A. Hemp, S. Gockel, M. Gorke, K. Lorenzen, K. Reichel-Jung, S. Renner, I. Steitz, F. Straub, M. Teuscher, J. Vogt, S. Weithmann, K. Wells and K. Wiesner (local management teams), C. Fischer, M. Gleisberg, J. Mangels and S. Pfeiffer (central office), B. König-Ries, J. Nieschulze, A. Ostrowski and M. Owonibi (central database management), M. Fischer, D. Hessenmöller, E. Linsenmair, J. Nieschulze, E-D. Schulze, and the late E. Kalko for their roles in setting up the Biodiversity Exploratories program (agricultural grasslands), which was funded by the Deutsche Forschungsgemeinschaft Priority Program 1374 Infrastructure-Biodiversity Exploratories. Fieldwork permits were given by the responsible state environmental offices of Baden-Württemberg, Thüringen, and Brandenburg (according to §72 BbgNatSchG). We are grateful to Barbara Schmitt, Stefan Blaser, Verena Busch, Helmut Hillebrand and Anne Ebeling for sharing their data.

## Author contributions

Y.O. and W.W. conceived the idea of this study. Y.O. and M.L. contributed equally to the data analysis and writing of the manuscript. C.R., F.A. (Alt), N.B., D.B., S.B., R.S.B., F.B., S.D., G.D.D., N.E., G.G., K.G., N.H., M.J., E.K., V.H.K., T.K., G.L.P., P.M., S.M., D.P., D.S., I.S., M.S., E.S., C.W., T.W., namely all authors but F.A. (Aburto) and S.L. contributed data. Y.O., S.L. and M.L. performed the analyses. Y.O. performed the literature search and wrote the first draft of the manuscript and all the authors contributed substantially to the revisions.

## Competing interests

The authors declare no competing interests.

## Additional information

**Supplementary information** The online version contains supplementary material available at https://doi.org/10.1038/s41467-021-24714-4.

**Correspondence** and requests for materials should be addressed to Y.O.

**Peer review information** *Nature Communications* thanks Ning Ling, Maarja Öpik, and the other, anonymous, reviewer(s) for their contribution to the peer review of this work. Peer reviewer reports are available.

**Reprints and permission information** is available at http://www.nature.com/reprints

**Publisher's note** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
<figure><img src="image_5.png" /></figure>

**Open Access** This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.

© The Author(s) 2021