# Top canopy nitrogen allocation linked to increased grassland carbon uptake in stands of varying species richness

Alexandru Milcu, Arthur Gessler, Christiane Roscher, Laura Rose, Zachary

Kayler, Dörte Bachmann, Karin Pirhofer-Walzl, Saša Zavadlav, Lucía

Galiano, Tina Buchmann, et al.

► To cite this version:

Alexandru Milcu, Arthur Gessler, Christiane Roscher, Laura Rose, Zachary Kayler, et al. Top canopy nitrogen allocation linked to increased grassland carbon uptake in stands of varying species richness. Scientific Reports, 2017, 7 (1), 10.1038/s41598-017-08819-9. hal-02480457

HAL Id: hal-02480457

https://hal.science/hal-02480457

Submitted on 3 Jun 2021

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
SCIENTIFIC REPORTS

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

OPEN

# Top canopy nitrogen allocation linked to increased grassland carbon uptake in stands of varying species richness

Received: 7 March 2017

Accepted: 19 July 2017

Published online: 16 August 2017

Alexandru Milcu<sup>1,2</sup>, Arthur Gessler<sup>3,4,5</sup>, Christiane Roscher<sup>6,7</sup>, Laura Rose<sup>8</sup>, Zachary Kayler<sup>3,5,12</sup>, Dörte Bachmann<sup>9</sup>, Karin Pirhofer-Walz<sup>1,3,5</sup>, Saša Zavadlav<sup>10</sup>, Lucia Galiano<sup>4</sup>, Tina Buchmann<sup>11</sup>, Michael Scherer-Lorenzen<sup>8</sup> & Jacques Roy<sup>1</sup>

Models predict that vertical gradients of foliar nitrogen (N) allocation, increasing from bottom to top of plant canopies, emerge as a plastic response to optimise N utilisation for carbon assimilation. While this mechanism has been well documented in monocultures, its relevance for mixed stands of varying species richness remains poorly understood. We used 21 naturally assembled grassland communities to analyse the gradients of N in the canopy using N allocation coefficients ($K_N$) estimated from the distribution of N per foliar surface area ($K_{N-F}$) and ground surface area ($K_{N-G}$). We tested whether: 1) increasing plant species richness leads to more pronounced N gradients as indicated by higher $K_N$-values, 2) $K_N$ is a good predictor of instantaneous net ecosystem CO$_2$ exchange and 3) functional diversity of leaf N concentration as estimated by Rao's Q quadratic diversity metric is a good proxy of $K_N$. Our results show a negative (for $K_{N-G}$) or no relationship (for $K_{N-F}$) between species richness and canopy N distribution, but emphasize a link (positive relationship) between more foliar N per ground surface area in the upper layers of the canopy (i.e. under higher $K_{N-G}$) and ecosystem CO$_2$ uptake. Rao's Q was not a good proxy for either $K_N$.

During the last two decades, substantial progress has been made in understanding the role of biodiversity for ecosystem functioning<sup>1-3</sup>. Numerous meta-analyses<sup>4, 5</sup> and biodiversity experiments such as the Cedar Creek experiment<sup>6</sup>, the pan-European Biodepth experiment<sup>7</sup> and the Jena Experiment<sup>8, 9</sup> have attested to a positive relationship between plant diversity and ecosystem functioning. However, while there is now consensus that biodiversity begets productivity due to complementarity<sup>10, 11</sup>, insurance<sup>12, 13</sup> and selection effects<sup>14</sup>, the underlying physiological mechanisms are not fully understood. Furthermore, there is increasing evidence that functional diversity is of greater importance for ecosystem functioning than the number of taxonomic species<sup>15, 16</sup>. Several functional diversity metrics currently quantify the variety, range and evenness of community traits and assume that a greater dissimilarity of traits indicates less niche overlap, and hence, more efficient capture of resources in

<footnote n="1"><sup>1</sup>CNRS, Ecotron (UPS-3248), Campus Baillarguet, F-34980, Montferrier-sur-Lez, France. <sup>2</sup>Centre d'Ecologie Fonctionnelle et Evolutive, CEFE-CNRS, UMR 5175, Université de Montpellier - Université Paul Valéry - EPHE, 1919 route de Mende, F-34293, Montpellier Cedex 5, France. <sup>3</sup>Institute for Landscape Biogeochemistry, Leibniz-Centre for Agricultural Landscape Research (ZALF), 15374, Müncheberg, Germany. <sup>4</sup>Swiss Federal Research Institute WSL, Zuercherstrasse 111, Birmensdorf, 8903, Switzerland. <sup>5</sup>Berlin-Brandenburg Institute of Advanced Biodiversity Research (BBIB), Altensteinstrasse 6, Berlin, 14195, Germany. <sup>6</sup>UFZ, Helmholtz Centre for Environmental Research, Department Physiological Diversity, Permoserstrasse 15, 04103, Leipzig, Germany. <sup>7</sup>German Centre for Integrative Biodiversity Research (iDiv) Halle-Jena-Leipzig, Deutscher Platz 5e, 04103, Leipzig, Germany. <sup>8</sup>University of Freiburg, Faculty of Biology, Geobotany, Schaenzlestr. 1, D-79104, Freiburg, Germany. <sup>9</sup>Institute of Agricultural Sciences, ETH Zurich, Universitaetsstrasse 2, 8092, Zurich, Switzerland. <sup>10</sup>Slovenian Forestry Institute, Department of Forest Physiology and Genetics and Department of Forest Yield and Silviculture, Večna pot 2, SI-1000, Ljubljana, Slovenia. <sup>11</sup>UFZ, Helmholtz Centre for Environmental Research, Department of Community Ecology, Theodor-Lieser-Strasse 4, 06120, Halle, Germany. <sup>12</sup>Present address: Department of Soil and Water Systems, University of Idaho, 875 Perimeter Dr., Moscow, 83844, ID, USA. Correspondence and requests for materials should be addressed to A.M. (email: alex.milcu@cnrs.fr)</footnote>
time and space<sup>17-19</sup>. Two recent studies<sup>20,21</sup> found the diversity of leaf nitrogen (N) concentration in the canopy, as measured by Rao's quadratic entropy (FD<sub>Q</sub>-N)<sup>18</sup>, to be a reliable predictor of net ecosystem CO<sub>2</sub> exchange (NEE), gross primary productivity and water use efficiency in grassland communities. The authors suggested that FD<sub>Q</sub>-N might be related to the vertical distribution of N in the canopy, with higher FD<sub>Q</sub>-N values indicating canopies with more pronounced N gradients in the canopy increasing from bottom to the top of the canopy. This conjecture is in line with the optimal N allocation hypothesis<sup>22-24</sup>, which states that a link between canopy N distribution and CO<sub>2</sub> uptake should be expected because canopies with increasing leaf N concentrations from bottom to top of the canopy, following the corresponding increase in light availability, should be more optimal. Such more optimal canopies should have higher nitrogen use efficiencies (NUE) and photosynthetic carbon (C) uptake because the gain in C uptake per unit N investment is greater at higher light intensities for a given N content<sup>23</sup>. However, to date, the evidence for a link between plant diversity, canopy N allocation and canopy-level CO<sub>2</sub> uptake remains circumstantial and has not been directly tested. Furthermore, the majority of studies that investigated the links between canopy N distributions and ecosystem C uptake used computer simulations, simplified communities (i.e., monocultures)<sup>22, 24, 25</sup> or used controlled environmental conditions with small containers, hydroponics or artificial light<sup>26-28</sup>, which cautions against the generality of their results.

Canopies of monocultures have often been found to have more uniform than optimal N distributions, but a closer to optimal N allocation has been found to lead to 1–42% higher C uptake due to increased NUE<sup>22–24, 29</sup>. In mixed stands of varying species richness, where dominant and subordinate species emerge, it is less clear how the interspecific competition for light affects the N allocation profiles and, consequently, its impact on the ecosystem C uptake. Wacker *et al.* (2009)<sup>30</sup> suggested that mixed stands may be more efficient in C uptake than monoculture stands because they form “integrated” canopies with complementary allocation of leaf mass and nitrogen along the vertical light profile, and that canopies including species with a larger range of leaf mass and N values may have a higher likelihood of assembling “integrated” canopies. However, several theoretical studies suggested the contrary, namely that mixed stands may be less efficient in C uptake due to a competition driven “tragedy of the commons” owing to the fact that the most competitive strategy at the individual level is to maximise its own C gain. This can result in suboptimal light utilisation across the entire canopy in multispecies stands<sup>31, 32</sup>. Yet, multispecies stands have often been found to be more productive in biodiversity experiments<sup>33, 34</sup>. Therefore, it is currently unclear whether species-rich stands are more efficient at matching the foliar N content with the emerging light attenuation profiles, and thus, achieve more pronounced vertical gradients of foliar N allocation leading to increased C gain.

To fill the aforementioned knowledge gaps, we took advantage of 21 communities varying in plant species richness that assembled naturally after cessation of weeding in subplots of the Jena Experiment<sup>35</sup>, i.e. where the control of species richness and composition was abandoned for at least seven years. Specifically, we test the following hypotheses: 1) increasing plant species richness leads to more pronounced vertical N gradients estimated with the N allocation coefficients (K<sub>N</sub>) according to the method of Hirose and Werger (1987)<sup>23, 2</sup>, K<sub>N</sub> is a good predictor of instantaneous NEE and 3) functional diversity of leaf N concentration as estimated by Rao's Q quadratic diversity metric is a good proxy of K<sub>N</sub> (see Table 1 for a list of abbreviations). To this end we tested two versions of K<sub>N</sub>, one based on the distribution of foliar N content per leaf surface area (K<sub>N-F</sub>), the other based on the distribution of foliar N content per ground surface area (K<sub>N-G</sub>).

## Results

Since the NEE results indicated a statistically significant difference between monospecific-dominated and mixed stands, but no plant diversity effects, we focus our presentation of the results on emphasizing differences between monospecific-dominated and mixed stands.

**Vertical distribution of foliar nitrogen and light in the canopies.** The relationship between the canopy height and the percentage of foliar N was best described by a 3-parameter exponential asymptotic, showing an increase of foliar N concentrations with height, which levels off at the top of the canopy. Based on the intercepts resulting from our exponential function fit to the data, foliar N concentrations in the lowest layer were predicted to be higher in the monospecific-dominated stands relative to mixed stands (1.78 vs. 1.16% N intercept; Fig. 1a). The coefficients indicating the saturating N concentrations towards the top of the canopies were 4.31% and 2.85% N for monospecific-dominated and mixed stands, respectively (Fig. 1a). The percentage of foliar N also increased linearly with the percentage of available light (Fig. 1b), with a slightly higher slope for monospecific-dominated stands (N% = 2.13 + 0.012*%light, R<sup>2</sup> = 0.19, P < 0.0001) relative to mixed stands (N% = 2.17 + 0.007*%light, R<sup>2</sup> = 0.20, P < 0.0001). The light attenuation coefficients (K<sub>L</sub>) varied from 0.56 to 0.90, with a mean of 0.67, and a median of 0.64 (see Supplementary Fig. S1). K<sub>L</sub> values were not significantly different between monospecific and mixed stands (F<sub>1, 19</sub> = 0.02, P = 0.875) and no significant effect of RSR (F<sub>1, 19</sub> = 0.83, P = 0.371) or SR15 (F<sub>1, 19</sub> = 0.13, P = 0.724) on K<sub>L</sub> was found. Similarly, no significant effect of the type of stand on LAI was found (5.99 vs. 5.57 for monospecific and mixed stands respectively; F<sub>1, 19</sub> = 0.30, P = 0.587). Furthermore, no significant effect of RSR (F<sub>1, 19</sub> = 0.44, P = 0.512) or SR15 (F<sub>1, 19</sub> < 0.01, P = 0.979) were observed.

The relationship between the foliar N content per ground surface area (g N m<sup>-2</sup> ground surface; N<sub>G</sub>) and the height in the canopy as well as the available light exhibited a humped shape (Fig. 1c). The distribution of N<sub>G</sub> in the canopy differed between the monospecific-dominated and mixed stands. Monospecific-dominated stands exhibited lower N<sub>G</sub> at 10 (−39%, t = −4.05, P = 0.001) and 20 cm (−39%, t = −2.24, P = 0.037) height and three- to four-fold higher N<sub>G</sub> in the higher strata at 60 (t = 2.74, P = 0.013) and 70 cm (t = 2.22, P = 0.039) height.

**Nitrogen allocation coefficients (K<sub>N</sub>).** To help visualising how the K<sub>N-F</sub> and K<sub>N-G</sub> values were obtained, in Fig. 2a and b we show examples of two plots with contrasting K<sub>N-F</sub> and K<sub>N-G</sub> values (plots no. 1 and 21 from Table S1) derived from the relationship between foliar N content (per ground or foliar surface) and the LAI at
<table><caption>**Table 1.** Table explaining the most important abbreviations.</caption><thead><tr><td>Abbreviation</td><td>Description</td><td>Unit</td></tr></thead><tbody><tr><td>F<sub>BM</sub></td><td>Foliar biomass (DW) per ground surface area</td><td>g DW m<sup>-2</sup></td></tr><tr><td>FD<sub>Q-N</sub></td><td>Functional diversity of leaf N concentrations calculated based on species-level averaged foliar N</td><td>unitless</td></tr><tr><td>K<sub>L</sub></td><td>Canopy light attenuation coefficient (eq. 2)</td><td>unitless</td></tr><tr><td>K<sub>N-F</sub></td><td>Nitrogen allocation coefficient in the canopy (eq. 1) calculated based on canopy N amount per foliar surface area</td><td>unitless</td></tr><tr><td>K<sub>N-G</sub></td><td>Nitrogen allocation coefficient in the canopy (eq. 1) calculated based on the N amount per ground surface area</td><td>unitless</td></tr><tr><td>LAI</td><td>Total leaf area index</td><td>m<sup>2</sup> leaf m<sup>-2</sup> ground</td></tr><tr><td>LAI<sub>D</sub></td><td>Cumulative LAI from the top of the canopy to the depth D</td><td>m<sup>2</sup> leaf m<sup>-2</sup> ground</td></tr><tr><td>LUE</td><td>Light use efficiency</td><td>µmol CO<sub>2</sub> m<sup>-2</sup> s<sup>-1</sup> / µmol PAR m<sup>-2</sup> s<sup>-1</sup></td></tr><tr><td>Mixed</td><td>Mixed stands containing several (2 to 5) dominant plant species</td><td>type of stand</td></tr><tr><td>Mono</td><td>Stands dominated by a single plant species</td><td>type of stand</td></tr><tr><td>NEE</td><td>Instantaneous net ecosystem CO<sub>2</sub> exchange at canopy level</td><td>µmol CO<sub>2</sub> m<sup>-2</sup> s<sup>-1</sup></td></tr><tr><td>N<sub>F</sub></td><td>Foliar N per leaf surface area</td><td>g N m<sup>-2</sup> leaf</td></tr><tr><td>N<sub>G</sub></td><td>Foliar N per ground surface area</td><td>g N m<sup>-2</sup></td></tr><tr><td>N<sub>R</sub></td><td>N content present in reproductive organs (mainly inflorescences)</td><td>g N m<sup>-2</sup></td></tr><tr><td>NUE</td><td>Nitrogen use efficiency</td><td>µmol CO<sub>2</sub> m<sup>-2</sup> s<sup>-1</sup> / g N m<sup>-2</sup></td></tr><tr><td>PAR</td><td>Photosynthetic active radiation measured at the top of the canopy</td><td>µmol m<sup>-2</sup> s<sup>-1</sup></td></tr><tr><td>RSR</td><td>Realised species richness including all species present (See Table S1).</td><td>count</td></tr><tr><td>SR15</td><td>Species richness including species with a ground surface cover higher than 15% (See Table S1).</td><td>count</td></tr><tr><td>T</td><td>Air temperature in the cuvette used for NEE measurements</td><td>°C</td></tr></tbody></table>

different canopy heights using eq. (2) (see Supplementary Figs S2 and S3 to see the fitted values for $K_{N-F}$ and $K_{N-G}$ in all plots). Across the 21 communities, the values of $K_{N-F}$ ranged between 0.25 and 3.12, with an overall mean of 1.01 and a median of 0.79. No significant relationship between $K_{N-F}$ and the tested predictors (species richness, type of stand, LAI, $N_G$, and $K_L$) was found.

In contrast to $K_{N-F}$, $K_{N-G}$ showed lower values, ranging from -0.11 to 0.33, with an overall mean of 0.08 and a median of 0.07 (see Supplementary Fig. S3). $K_{N-G}$ was best predicted by models incorporating leaf N content per ground surface area ($N_G$, with a positive coefficient, Fig. 2d) and realised species richness (RSR, with a negative coefficient, Fig. 2e) (Table 2). No effect of Simpson, Shannon diversity metrics nor the evenness metric was found on any of the two N allocation coefficients (Supplementary Fig. S6).

**Relationship between stand carbon uptake and leaf nitrogen allocation.** Multiple regression analyses (Table 3) indicate that NEE was higher in mixed stands relative to monospecific ones (Fig. 3a). No significant effect of RSR or SR15 on NEE was found (Supplementary Fig. S5). Similarly, no significant effect of Shannon and Simpson diversity metrics or community evenness affected NEE (Supplementary Fig. S6). As we found a significant increase of NEE in mixed stands concomitant with a tendency of lower $K_{N-G}$ values, we further analysed the relationship between the N allocation coefficients and NEE in all communities as well as in mixed stands only.

When the monoculture and mixed stands were analysed together, the model explaining best the NEE ($R^2$ = 0.59, P = 0.001) included an interaction between the percentage of legumes in the canopy and type of stand as well as an interaction between the percentage of legumes and canopy leaf biomass ($F_{BM}$) (Table 3). Overall, $F_{BM}$ was a consistent predictor of NEE. $K_{N-G}$ was also retained among the five best models. However, when we focused on mixed stands only, the best model predicting NEE incorporated $K_{N-G}$ and an interaction between $K_{N-G}$ and $N_G$ ($R^2$ = 0.64, P = 0.002; Table 3 and Fig. 3b).

**Effects of N allocation on the light and nitrogen use efficiencies of NEE.** To explore the potential relationships between canopy N allocation and the nitrogen use efficiency (NUE) and light use efficiency (LUE) of instantaneous NEE while avoiding the intrinsic species-specific differences occurring in monocultures, we focused on mixed stand only. LUE was best predicted by a model that included a positive coefficient for $K_{N-G}$ as well as a negative temperature (T) effect (LUE = 0.018 - 0.0003*T + 0.0015*K_{N-G}*N_G; P = 0.001, $R^2$ = 0.68). In contrast, the model for nitrogen use efficiency (NUE) was best predicted by $K_{N-G}$, with a negative regression coefficient and the interaction between temperature and the percentage of legumes in the canopy, with a negative coefficient (NUE = 8.80 - 0.88*N_G - 0.008*LAI*T; P < 0.001, $R^2$ = 0.84). See Supplementary Fig. S7 for a scatter-plot and correlation matrix depicting the relationships between NUE, LUE and their most important predictors.

**Discussion**

To our knowledge, this is the first study assessing the potential links between canopy N allocation coefficients ($K_N$) and net ecosystem CO$_2$ exchange (NEE) in stands of varying species richness under field conditions. Concerning our first hypothesis, we found that increasing plant species richness did not lead to canopies with higher N
<figure><img src="image_3.png" /><figcaption>**Figure 1.** Relationship between the foliar nitrogen (N) concentration and (a) height in the canopy and (b) percentage of light transmission in monospecific-dominated and mixed stands. Relationship between leaf N content per ground surface area and (c) height in the canopy and (d) percentage of available light.</figcaption></figure>

distribution coefficients compared to monospecific-dominated stands. Instead, we found that increasing the number of species present in the canopy leads to lower $K_{N-G}$ values, meaning that relatively less N was allocated at the top of the canopy where light availability is higher at higher diversity levels. These results are in line with the studies of Anten<sup>29</sup> and Wacker *et al.*<sup>30</sup> documenting that species in mixed stands exhibit non-optimal N allocation for increasing C gain at the whole stand level. The latter study suggested that optimal N profiles in mixed stands emerge only if the coexisting species had already complementary N profiles based on the among-species variation in monocultures. This can be explained by the fact that N cannot be allocated between species (e.g. a small size species with high N content cannot offer N to a tall species), while N can be allocated to different leaves and heights within an individual. Alternatively, Hikosaka (2014) proposed that lower $K_N$ values could result from a lower predictability of light availability in dense canopies. This could be a potential explanation for the decreasing $K_{N-G}$ values with increasing species richness as the competition for light over time and during stand development might be less predictable in species rich communities compared to monospecific stands. Therefore, individual plants in species rich stands might be less able to project which leaves will receive more direct light than others. In other words, if the position and duration of direct light is less predictable in species rich stands, plants will benefit less from allocating more N towards the upper layers of the canopies<sup>36</sup>. Taken together, our results support the tragedy of the commons phenomena<sup>29,31</sup>, in which light acquisition strategies of species optimise for the individual rather than for the entire canopy.

Although vertical gradients of leaf N allocation are a common feature of plant canopies, the factors controlling the formation of non-uniform N distributions are still debated. Theory predicts that canopy N profiles are the result of the interaction between N availability and the light attenuation profiles<sup>22, 23, 37</sup>. Several studies suggested that more pronounced gradients of N allocation in the canopy (i.e. higher $K_N$ values) will be achieved if more N is available for allocation and redistribution, and a positive correlation between total canopy N content and $K_N$ has been reported<sup>23, 38</sup>. However, two other studies<sup>26, 28</sup> showed an inverse trend, namely a more uniform N allocation with increased N supply, but these two studies where performed in small pots/containers grown in greenhouses and growth chambers with artificial lighting, conditions arguably different from natural field conditions. In our study, with canopies naturally assembled in the field, we found that $K_{N-G}$ was consistently predicted (with positive regression coefficients, Table 2) by the total amount of foliar N per ground surface area ($N_G$). However, the attenuation coefficients ($K_L$) also proved to be important for $K_{N-G}$ and, as predicted by theoretical work<sup>22, 24, 37</sup>, the
<figure><img src="image_4.png" /><figcaption>**Figure 2.** (a) Examples of two plots (no. 1 and 21 from Table S1) with contrasting K<sub>N-F</sub> values estimated with eq. (2) from the relationship between foliar N per leaf surface area and leaf area index at different depths (LAI<sub>D</sub>) (b) Same as for Fig. 2a but for the N allocation coefficient based on foliar N per ground surface area (K<sub>N-G</sub>), (c) Boxplots depicting the median and distributions of K<sub>N-F</sub>, K<sub>N-G</sub> and light attenuation coefficients (K<sub>L</sub>) in monospecific and mixed stands. (d) Relationship between the total foliar N content per ground surface area (N<sub>G</sub>) and K<sub>N-G</sub>. (e) Relationship between realised species richness (RSR) and K<sub>N-G</sub>. (f) Relationship realised species richness (RSR) and K<sub>N-F</sub>.</figcaption></figure>

interaction between K<sub>L</sub> and N<sub>G</sub> was also a significant predictor of K<sub>N-G</sub> (Table 2). Overall, our results support the
theoretical findings of Hirose and Werger (1987)<sup>23</sup> indicating that the allocation of N in the canopy is influenced
by the interaction between N availability and light extinction profiles.

Two recent papers put forth the hypothesis that a more optimal N allocation in the canopy could occur in mix-
tures with higher species richness, and that this might be part of the mechanistic basis to explain the higher C and
water fluxes observed in higher diversity mixtures<sup>20, 21</sup>. The supporting evidence for this conjecture relied on the
<table><caption>**Table 2.** The five candidate models explaining the nitrogen (N) distribution coefficient per ground surface area (K<sub>N-G</sub>) as a function of light attenuation profiles (K<sub>L</sub>), realised species richness (RSR), total leaf N (N<sub>G</sub>) and the N content allocated to reproductive organs (N<sub>R</sub>). No significant predictors were found for K<sub>N-F</sub>. AICc represents the Akaike Information Criterion adjusted for sample size whereas Akaike weights (AICw) represents the probability that a particular model is the best fit to the data. See Table 1 for an abbreviation list.</caption><thead><tr><th>K<sub>N-G</sub></th><th></th><th>AICc</th><th>AICc</th><th>R<sup>2</sup></th><th>P-value</th></tr></thead><tbody><tr><td>1.</td><td>K<sub>N-G</sub> = 0.16 + 0.05*K<sub>L</sub>*N<sub>G</sub> - 0.03*K<sub>L</sub>*RSR</td><td>-36.61</td><td>0.13</td><td>0.56</td><td>&gt;0.001</td></tr><tr><td>2.</td><td>K<sub>N-G</sub> = 0.14 + 0.03*N<sub>R</sub>*N<sub>G</sub> - 0.02*N<sub>R</sub>*RSR</td><td>-36.30</td><td>0.11</td><td>0.56</td><td>&gt;0.001</td></tr><tr><td>3.</td><td>K<sub>N-G</sub> = 0.14 + 0.03*N<sub>G</sub> - 0.02*RSR</td><td>-36.15</td><td>0.10</td><td>0.56</td><td>&gt;0.001</td></tr><tr><td>4.</td><td>K<sub>N-G</sub> = 0.13 + 0.03*N<sub>G</sub> - 0.02*RSR*K<sub>L</sub></td><td>-36.11</td><td>0.10</td><td>0.55</td><td>&gt;0.001</td></tr><tr><td>5.</td><td>K<sub>N-G</sub> = 0.27-0.03*RSR + 0.003*RSR*N<sub>G</sub></td><td>-35.90</td><td>0.09</td><td>0.55</td><td>&gt;0.001</td></tr></tbody></table>

<table><caption>**Table 3.** The five best models predicting the instantaneous CO₂ net ecosystem exchange (NEE) as a function of nitrogen (N) allocation confidents (K<sub>N-F</sub> and K<sub>N-G</sub>), foliar biomass (F<sub>BM</sub>), percentage of legume cover (Leg) and type of stand (Mono = monospecific-dominated stands, Mix = mixed stands). AICc represents the Akaike Information Criterion adjusted for sample size whereas Akaike weights (AICw) represents the probability that a particular model is the best fit to the data. See Table 1 for a list of all abbreviations.</caption><thead><tr><th>No</th><th>Model</th><th>AICc</th><th>AICw</th><th>R<sup>2</sup></th><th>P-value</th></tr></thead><tbody><tr><td colspan="6">All stands (n = 21)</td></tr><tr><td>1</td><td>NEE = 13.81 + 2.88*Leg*Mono + 1.99*Leg*Mix - 0.16*Leg*F<sub>BM</sub></td><td>63.92</td><td>0.41</td><td>0.59</td><td>0.001</td></tr><tr><td>2</td><td>NEE = 14.17 + 2.57*Leg - 0.02*Leg*F<sub>BM</sub></td><td>65.75</td><td>0.17</td><td>0.48</td><td>0.003</td></tr><tr><td>3</td><td>NEE = 6.87 + 0.059*F<sub>BM</sub> - 0.003*N<sub>G</sub>*F<sub>BM</sub></td><td>66.34</td><td>0.02</td><td>0.37</td><td>0.015</td></tr><tr><td>4</td><td>NEE = 8.56 + 0.03*F<sub>BM</sub> - 0.17*Leg*N<sub>G</sub></td><td>66.70</td><td>0.02</td><td>0.37</td><td>0.015</td></tr><tr><td>5</td><td>NEE = 8.05 + 0.04*F<sub>BM</sub> - 1.78*F<sub>BM</sub>*K<sub>N-G</sub></td><td>67.79</td><td>0.02</td><td>0.37</td><td>0.016</td></tr><tr><td colspan="6">Mixed stands only (n = 15)</td></tr><tr><td>1</td><td>NEE = 12.91-73.98*K<sub>N-G</sub> + 20.78*K<sub>N-G</sub>*N<sub>G</sub></td><td>61.32</td><td>0.18</td><td>0.64</td><td>0.002</td></tr><tr><td>2</td><td>NEE = 13.43-48.89*K<sub>N-G</sub> + 0.34*K<sub>N-G</sub>*F<sub>BM</sub></td><td>61.61</td><td>0.16</td><td>0.47</td><td>0.020</td></tr><tr><td>3</td><td>NEE = 13.19-0.005*Leg*F<sub>BM</sub></td><td>63.61</td><td>0.06</td><td>0.21</td><td>0.084</td></tr><tr><td>4</td><td>NEE = 9.71 + 0.02*F<sub>BM</sub> - 0.006*Leg*F<sub>BM</sub></td><td>63.67</td><td>0.06</td><td>0.38</td><td>0.053</td></tr><tr><td>5</td><td>NEE = 8.23 + 0.03*F<sub>BM</sub> - 0.94*Leg</td><td>63.80</td><td>0.05</td><td>0.37</td><td>0.059</td></tr></tbody></table>

detected positive correlations between carbon and water fluxes and a functional diversity metric (Rao's quadratic entropy of leaf N concentrations, FD<sub>Q</sub>-N), which captures the diversity/dissimilarity of foliar N concentrations between the species of a mixture. Here, we directly tested the hypothesis that that FD<sub>Q</sub>-N could be a good proxy of K<sub>N</sub> and found no supporting evidence that FD<sub>Q</sub>-N is positively correlated with K<sub>N-F</sub> and K<sub>N-G</sub>. However, one caveat needs to be noted relative to the previous studies proposing FD<sub>Q</sub>-N as a good proxy of K<sub>N</sub>. In this study, the calculation of FD<sub>Q</sub>-N was based on values of foliar N measured in plants collected *in situ* but not directly from the plots included in our experiment. Although the literature is abundant with functional trait-based metrics calculated with literature-derived or non *in-situ* measured trait values<sup>39-41</sup>, the fact that FD<sub>Q</sub>-N was calculated with foliar N values not measured in our communities remains a notable caveat in this study because there is evidence that plant species change their trait values in response to the growth conditions experienced in mixtures varying in diversity<sup>42, 43</sup>.

Up to now, the impact of N allocation on canopy-level C gain relied exclusively on canopy models to upscale the relationship between leaf photosynthetic capacity (Amax) and leaf N content (Supplementary Fig. S8) to canopy values<sup>23, 37</sup>. Here, we opted to test the relationship between canopy C uptake and N allocation coefficients using instantaneous NEE because midday values of instantaneous NEE measured in the same experimental system were vastly dominated by photosynthetic activity<sup>20</sup>. Although no significant impact of species richness was found on the NEE, our results provide supporting evidence for our third hypothesis postulating that the allocation of N in the canopy is important for whole canopy C gain. Whilst the control of NEE is complex because is affected by multiple factors, including (but not limited to) phenological status, temperature, soil moisture and light intensity, here we found that when more pronounced gradients of canopy N allocation (increasing from bottom to top of the canopy) were found in mixed stands as indicated by higher K<sub>N-G</sub> values, it led to higher NEE and LUE. This relationship was weaker when monocultures were included, but this can be explained by the intrinsic species-specific photosynthetic efficiencies of monocultures which presumably override the effect of the N allocation coefficients.

One important finding of this study is that of the two estimated N distribution coefficients, the one expressing the canopy N distribution per ground surface area (K<sub>N-G</sub>) was actually more relevant for canopy level C uptake than the distribution of canopy N based on foliar surface area (K<sub>N-F</sub>) (Table 3). Although there is substantial theoretical and empirical evidence that K<sub>N-F</sub> (nitrogen allocation coefficient for leaf area-based leaf N) is a reliable predictor of the photosynthetic uptake of individual plant<sup>23, 24, 44</sup> and monospecific stands, its predictive power for stand-level C fluxes in mixed stands was not directly demonstrated by *in-situ* measurements. However, K<sub>N-F</sub> does
<figure><img src="image_5.png" /><figcaption>**Figure 3.** (a) Boxplots showing the median and variation in instantaneous CO₂ net ecosystem exchange (NEE) measured in monocultures and mixed stands. Multiple regression results from Table 3 show significantly different (P < 0.05, *) regression coefficients for monospecific and mixed stands. (b) Relationship between predicted NEE values by the best model from Table 3 and measured NEE in mixed stands.</figcaption></figure>

not take into account the leaf biomass distribution by height, and consequently, a high N content per leaf area at the top of the canopy leading to a high K<sub>N-F</sub> value only indicates that the leaves at the top of the canopy have high N content, but it does not indicate the amount of foliar biomass, an important determinant of stand-level C uptake. In contrast, K<sub>N-G</sub> takes into account the leaf biomass distribution by height as a high N content per ground surface area at the top of the canopy leading to a high K<sub>N-G</sub> value indicates that there is a significant amount of foliar biomass and foliar N at the top of the canopy where light is less limiting. Furthermore, we argue that K<sub>N-F</sub> is prone to a strong bias introduced by the N concentration of the leaves situated at the very top of the canopy which often represent a minor proportion of the total LAI, and therefore, will have a relatively low contribution to the overall canopy level photosynthetic activity. The K<sub>N-F</sub> values can be overinflated when the leaves of the top of the canopy have disproportionally higher N concentrations relative the rest of the canopy. An underestimation of the N gradient in the canopy can also occur when the top leaves have disproportionally lower N due to reallocation to reproductive organs. K<sub>N-G</sub>, on the other hand, captures the distribution of total foliar N per ground surface area across the height of the whole canopy, and our results indicate that is a much better predictor of stand-level NEE than K<sub>N-F</sub>.

In conclusion, using naturally assembled canopies varying in species richness in field conditions, this study shows that increasing species richness leads to less pronounced vertical gradients (increasing from top to bottom of the canopy) of N allocation in the canopy. However, when pronounced vertical gradients of N allocation (i.e. with higher K<sub>N-G</sub> values) occurred in the canopies of mixed stands, the positive relationship between K<sub>N-G</sub> values and canopy carbon uptake hold true, leading to increased canopy-level NEE and LUE. This knowledge can be important in the selection of high-performance intercropping systems with more optimal N allocations in the canopy, which will help to increase productivity.

## Methods

**Field site and sampling strategy.** The site of the Jena Experiment (50° 57.1′ N, 11° 37.5′ E, 130 m above sea level; mean air annual temperature 9.3 °C, mean annual precipitation 587 mm<sup>45</sup>, is located on the floodplain of the Saale River (Jena, Germany), and was a former arable field until 2000. After two years of fallow, in May
2002, 82 plots (20 × 20 m) varying in sown plant species richness (1, 2, 4, 8, 16 and 60 species), plant functional groups (1 to 4, grasses, small herbs, tall herbs and legumes), and plant species composition were established<sup>46</sup>. The 82 plots were randomly allocated to four blocks that were identified to vary in soil texture<sup>46</sup>. Sown plant species composition was sustained by regular weeding of all unwanted species. This study was performed in 21 subplots (4.5 × 6.5 m) with plant communities designated initially for an invasion/colonisation experiment<sup>47</sup> and which, at the time of our study, were not weeded for seven years. These subplots were preferred to the maintained experimental gradients of species richness of the Jena Experiment because we sought communities with self-assembled intact canopies that were not disturbed by weeding. Furthermore, these communities assembled under similar environmental conditions after the initial sowing, had all developed closed canopies, but were still different in species richness, composition and functional diversity<sup>35</sup>, which was an important prerequisite to test our hypotheses. A representative area of 1 m<sup>2</sup> was selected in each plot and used for the measurements of leaf area index (LAI) and light extinction profiles, while a central circular area of 0.196 m<sup>2</sup> in each plot was used for NEE, aboveground biomass measurements and species abundances (estimated based on surface cover using a modified decimal scale; &lt;1%, 1–5%, 6–15%, 16–25%, etc.). As the community composition in the invasion plots was not perfectly homogeneous, the selected representative areas (of 1 m<sup>2</sup>) were selected following two a priori defined criteria: 1) to avoid areas that had a strongly different community composition than the majority of the selected plot area, and 2) communities developed closed canopies, i.e. with minimum surface of bare ground. These communities contained between 4 and 19 realised species richness (RSR) near equally distributed in 3 of the 4 blocks of the Jena Experiment (8 plots from Block I, 6 plots from Block II and 7 plots from Block III). In addition to RSR, we derived an additional species richness metric by including species with a surface cover higher than 15% (henceforth SR15). We argue that this threshold is relevant to canopy performance measurements such as the instantaneous CO<sub>2</sub> net ecosystem exchange (NEE) which is controlled by the dominant species of a canopy. SR15 was also used to define monospecific (n = 5) and mixed stands (n = 15); see also Supplementary Table S1. In addition to RSR and SR15 we also present the results for three classical diversity metrics (Simpson, Shannon and Pielou's evenness) as supplementary materials.

**Measurements of net ecosystem CO<sub>2</sub> exchange.** The relationship between canopy C uptake and N allocation coefficients was analysed using instantaneous (NEE), which represents the difference between gross primary production (GPP) and ecosystem respiration (Reco). Whilst GPP (where GPP = NEE – Reco) represents the actual canopy C uptake, we used NEE as a proxy for canopy C uptake because midday values of NEE measured in the same experimental system were vastly dominated by GPP (due to relatively low values of Reco)<sup>20</sup>. Two simultaneously calibrated (with the same CO<sub>2</sub> calibration gas) infra-red gas analysers (IRGAs), one LI-6400 (LI-COR Environmental, USA), and one Walz GFS-3000 (Heinz Walz GmbH, Germany) were used in connection with two large cylindrical (0.5 m dia. and 1.3 m height) UV-VIS transparent cuvettes made of polycarbonate to measure NEE. The measurements of CO<sub>2</sub> drawdown were used to calculate the instantaneous NEE as µmol CO<sub>2</sub> per m<sup>-2</sup> s<sup>-1</sup> using the ideal gas law. For these measurements, the chambers were closed and the IRGAs were operated in a closed loop setup for 5 minutes (cf. Volkmann et al.<sup>48</sup>). The measurements were performed during two consecutive days (20<sup>th</sup> and 21<sup>st</sup> of May 2014) with clear sky (i.e. sunny conditions) between 11am and 15 pm local time. During the measurements of NEE, the temperature and pressure in the cuvettes were recorded with the ancillary temperature and pressure sensors of the IRGAs. Incoming radiation values (PAR) where provided by the weather station of the Jena Experiment. Furthermore, the NEE measurements together with radiation and foliar N measurements allowed us to also estimate the light use efficiency (LUE; NEE per available PAR) and nitrogen use efficiency (NUE; NEE per total foliar N) of the measured NEE.

**Measurements of leaf area index (LAI) and light extinction coefficients (K<sub>L</sub>).* An SS1-SunScan canopy analyser coupled with a BF5 sunshine sensor (Delta-T Devices, Cambridge, UK) was used to estimate the light extinction profiles and the leaf area index (LAI) of the canopies. The SS1 canopy analyser was chosen due to its advantageous dimensions (1 m long by 13mm wide probe with PAR sensors each 15.6mm) which allowed for the estimation of LAI and light extinction profile of the selected 1 m<sup>2</sup> plots at seven heights (at 3, 10, 20, 30, 40, 50 and 60 cm) by averaging 5 readings (20 cm horizontally apart) for each height. This allowed for the estimation of the light attenuation coefficients through the canopy as estimated using the Beer's law<sup>49</sup>:

$$I = Io \exp(-K_L LAI_D) \quad (1)$$

where Io and I are the photon flux density on a horizontal plane above the canopy and within the canopy, respectively, at a given leaf area index value cumulated from the top of the canopy (LAI<sub>D</sub>); K<sub>L</sub>, is the extinction coefficient.

**Aboveground plant biomass and N measurements.** Following the NEE measurements, the aboveground biomass included in the cuvette was harvested by layer, from the top towards the bottom of the canopy at 10 cm intervals, and separated into leaves, stems and inflorescences. The different components of the biomass were weighed after drying at 70 °C for two days until constant weight. The leaf biomass by layer as well as the pooled stem and inflorescence biomass were milled into a homogeneous powder and analysed for N content using an elemental analyser (Truspec CNS; LECO Instrumente GmbH, Mönchengladbach, Germany).

**Estimation of the nitrogen allocation coefficients (K<sub>N</sub>).** The measurements of leaf N content by layer were used to estimate the coefficient of N allocation in the canopy (K<sub>N</sub>) following the method of Hirose and Werger (1987)<sup>23</sup>. The following exponential function was used to estimate K:
$$N_D=N \exp(-K_N LAI_D) \quad (2)$$

where N<sub>D</sub> and N represent leaf nitrogen content per unit area within the canopy at depth D and at the top uppermost layer of the canopy, respectively and LAI<sub>D</sub> represents the cumulative leaf area from the top of the canopy to the depth D. Larger positive K<sub>N</sub> values signify a less uniform N allocation in which upper leaves in the canopy have higher N concentrations than lower leaves to match the increased light availability. K<sub>N</sub> = 0 indicates a uniform allocation of leaf N in the canopy, in which every leaf has a N content equal to the mean. K<sub>N</sub> < 0 signifies that leaves that are lower in the canopy have higher N content. Following this approach two versions of K<sub>N</sub> have been estimated. One, where canopy N was expressed per foliar surface area (K<sub>N-F</sub>), which should capture the vertical distribution of the N per leaf area, and a second in which N was expressed per ground surface area (K<sub>N-G</sub>). We expect K<sub>N-G</sub> to be a better predictor of canopy NEE because it captures the vertical distribution of the total N amount in the canopy. To minimize biasing the K<sub>N</sub> values with data points from the highest canopy layers, which have a minimal LAI and hence a minimal contribution to the carbon fluxes but widely varying N concentrations, we only included data when the cumulative LAI from the top of the canopy was higher than 5% of the total LAI.

**Functional diversity metric based on leaf nitrogen concentrations.** We used Rao's quadratic entropy¹⁸ to calculate an index of diversity of leaf N concentrations present in the canopy which was previously found to be a good predictor of carbon and water fluxes²⁰,²¹. FD<sub>Q</sub>-N incorporates information about functional distance as well as functional evenness of a community because it is abundance weighted using species-specific plant cover data. In addition to FD<sub>Q</sub>-N, we also estimated functional richness (FRic) representing the amount of functional space filled by the community and functional evenness (FEve) which describes the evenness of abundance distribution of foliar N concentrations⁵⁰. These metrics where calculated based on the average leaf N concentration per species (estimated from pooling leaves from all canopy heights) as measured in situ in 2013 (with the exception of three species measured in 2015). They were computed using the “FD” package¹⁷ available through the R⁵¹ statistical package version 3.1.3.

**Statistical analyses.** All statistical analyses were performed in R statistical package version 3.1.3. To describe the shape of the relationship between the height of the canopy and the leaf N% (or N content per surface area) on the graphs we tested linear (y = a + bx), 2-parameter (y = a (-e<sup>-bx</sup>)) and 3-parameter (y = a - be<sup>-cx</sup>) exponential functions as well as locally weighted scatter-plot smoothers (loess). The exponential functions were fitted within the “nls” package and the best fit based on Akaike Information Criterion (AIC) was plotted in the graphs. T-tests were performed to identify differences between N content at different heights between monospecific-dominated and mixed stands.

To understand what drives the N allocation in the canopy we assembled multi-predictor linear models to test whether the N allocation coefficients (K<sub>N-F</sub> and K<sub>N-G</sub>) were affected by the type of stand (monospecific-dominated vs. mixed stands), species richness (SR), FD<sub>Q</sub>-N, FRic, FEve, total N content per ground surface area (N<sub>G</sub>), and total N content per leaf area (N<sub>F</sub>), with the soil sand content as covariable (see Table 1 for a list of all abbreviations). To account for potential effects of N reallocation in vegetative and reproductive organs on the K<sub>N</sub>-values, we also included N<sub>R</sub> (the N found in reproductive organs such as inflorescences) and stem biomass as covariables.

As predictors for NEE we tested the impact of the N allocation coefficients alongside N<sub>G</sub> and N<sub>F</sub> total leaf biomass (F<sub>BM</sub>) and the LAI, with soil sand content, cuvette measured radiation and temperature as covariables. Since in the Jena Experiment the soil texture was used to define blocks⁴⁶, here we took into account any block-related effects by introducing the sand content (decreasing from 45% in block 1 to 5% in block 4) as a covariable in all models. All initial/maximal models were simplified to achieve most parsimonious models using the “glmulti” statistical package⁵² for automated model selection and model-averaging. This package offers the possibility to identify the most parsimonious models (including parameter interactions) based on Akaike Information Criterion adjusted for sample size (AIC<sub>0</sub>)⁵³.

## References

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

2. Allan, E. et al. A comparison of the strength of biodiversity effects across multiple functions. *Oecologia* **173**, 223–37 (2013).

3. Duffy, J. E. *Why biodiversity is important to the functioning of real-world ecosystems*. **7**, 437–444 (2008).

4. Reich, P. B. et al. Impacts of biodiversity loss escalate through time as redundancy fades. *Science* (80–.). **336**, 589–592 (2012).

5. Balvanera, P. et al. Quantifying the evidence for biodiversity effects on ecosystem functioning and services. *Ecol. Lett.* **9**, 1146–1156 (2006).

6. Tilman, D. et al. Diversity and productivity in a long-term grassland experiment. *Science* (80–.). **294**, 843–845 (2001).

7. Spehn, E. M. et al. Ecosystem effects of biodiversity manipulations in European grasslands. *Ecol. Monogr.* **75**, 37–63 (2005).

8. Scherber, C. et al. Bottom-up effects of plant diversity on multitrophic interactions in a biodiversity experiment. *Nature* **468**, 553–556 (2010).

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

10. Tilman, D. et al. Plant diversity and ecosystem productivity: theoretical considerations. *Proc. Natl. Acad. Sci. USA* **94**, 1857–61 (1997).

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

12. Yachi, S. & Loreau, M. Biodiversity and ecosystem productivity in a fluctuating environment: the insurance hypothesis. *Proc. Natl. Acad. Sci. USA.* **96**, 1463–1468 (1999).

13. Allan, E. et al. More diverse plant communities have higher functioning over time due to turnover in complementary dominant species. *Proc. Natl. Acad. Sci.* **108**, 17034–17039 (2011).

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

15. Cadotte, M., Cavender-Bares, J., Tilman, D. & Oakley, T. Using phylogenetic, functional and trait diversity to understand patterns of plant community productivity. *PLoS One* **4**, e5695 (2009).
16. Roscher, C. et al. Using plant functional traits to explain diversity-productivity relationships. *PLoS One* 7, e36760 (2012).

17. Laliberté, E. & Shipley, B. F. D.: Measuring functional diversity (FD) from multiple traits, and other tools for functional ecology. *R Packag. version 1.0 9*, (2010).

18. Botta Dukát, Z. Rao's quadratic entropy as a measure of functional diversity based on multiple traits. *J. Veg. Sci.* **16**, 533–540 (2005).

19. Díaz, S. et al. Incorporating plant functional diversity effects in ecosystem service assessments. *Proc. Natl. Acad. Sci. USA.* **104**, 20684–20689 (2007).

20. Milcu, A. et al. Functional diversity of leaf nitrogen concentrations drives grassland carbon fluxes. *Ecol. Lett.* **17**, 435–444 (2014).

21. Milcu, A. et al. Plant species and functional diversity increase grassland productivity-related water vapour fluxes: a combined Ecotron and modeling approach. *Ecology* **97**, 2044–2054 (2016).

22. Field, C. Allocating leaf nitrogen for the maximization of carbon gain: leaf age as a control on the allocation program. *Oecologia* **56**, 341–347 (1983).

23. Hirose, T. & Werger, M. J. A. Maximizing daily canopy photosynthesis with respect to the leaf nitrogen allocation pattern in the canopy. *Oecologia* **72**, 520–526 (1987).

24. Anten, N., Schieving, F. & Werger, M. Patterns of light and nitrogen distribution in relation to whole canopy carbon gain in C3 and C4 mono-and dicotyledonous species. *Oecologia* **101**, 504–513 (1995).

25. Hirose, T., Werger, M. J. A., Pons, T. L. & van Rheenen, J. W. A. Canopy structure and leaf nitrogen distribution in a stand of Lysimachia vulgaris L. as influenced by stand density. *Oecologia* **77**, 145–150 (1988).

26. Dreccer, M. F., van Oijen, M., Schapendonk, A. H. C. M., Potvin, C. S. & Rabbinge, R. Dynamics of vertical leaf nitrogen distribution in a vegetative wheat canopy. *Impact on canopy photosynthesis*. *Ann. Bot.* **86**, 821–831 (2000).

27. Hirose, T. & Bazzaz, F. A. Trade-off between light- and nitrogen-use efficiency in canopy photosynthesis. *Ann. Bot.* **82**, 195–202 (1998).

28. Lötscher, M., Stroh, K. & Schnyder, H. Vertical leaf nitrogen distribution in relation to nitrogen status in grassland plants. *Ann. Bot.* **92**, 679–688 (2003).

29. Anten, N. P. R. Optimal photosynthetic characteristics of individual plants in vegetation stands and implications for species coexistence. *Ann. Bot.* **95**, 495–506 (2005).

30. Wacker, L., Baudois, O., Eichenberger-Glinz, S. & Schmid, B. Effects of plant species richness on stand structure and productivity. *J. Plant Ecol.* **2**, 95–106 (2009).

31. Schieving, F. & Poorter, H. Carbon gain in a multispecies canopy: the role of specific leaf area and photosynthetic nitrogen-use efficiency in the tragedy of the commons. *New Phytol.* **143**, 201–211 (1999).

32. Anten, N. P. & During, H. J. Is analysing the nitrogen use at the plant canopy level a matter of choosing the right optimization criterion? *Oecologia* **167**, 293–303 (2011).

33. Roscher, C. et al. Overyielding in experimental grassland communities - irrespective of species pool or spatial scale. *Ecol. Lett.* **8**, 419–429 (2005).

34. Hector, A. et al. Overyielding in grasslands communities: testing the sampling effect hypothesis with replicated biodiversity experiments. *Ecol. Lett.* **5**, 502–511 (2002).

35. Roscher, C., Schumacher, J., Gerighausen, U. & Schmid, B. Different assembly processes drive shifts in species and functional composition in experimental grasslands varying in sown diversity and community history. *PLoS One* **9**, e101928 (2014).

36. Hikosaka, K. Optimal nitrogen distribution within a leaf canopy under direct and diffuse light. *Plant, Cell Environ.* **37**, 2077–2085 (2014).

37. Werger, M. & Hirose, T. Leaf nitrogen distribution and whole canopy photosynthetic carbon gain in herbaceous stands. *Vegetatio* **97**, 11–20 (1991).

38. Hirose, T., Werger, M. J. A., Pons, T. L. & van Rheenen, J. W. A. Canopy structure and leaf nitrogen distribution in a stand of Lysimachia vulgaris L. as influenced by stand density. *T. Oecologia* **77**, 145–150 (1988).

39. Kattge, J. et al. TRY - a global database of plant traits. *Glob. Chang. Biol.* **17**, 2905–2935 (2011).

40. de Vries, F. T. et al. Abiotic drivers and plant traits explain landscape-scale patterns in soil microbial communities. *Ecol. Lett.* **15**, 1230–1239 (2012).

41. Manning, P. et al. Simple measures of climate, soil properties and plant traits predict national-scale grassland soil carbon stocks. *J. Appl. Ecol.* **52**, 1188–1196 (2015).

42. Gubsch, M. et al. Differential effects of plant diversity on functional trait variation of grass species. *Ann. Bot.* **107**, 157–169 (2011).

43. Lipowsky, A. et al. Plasticity of functional traits of forb species in response to biodiversity. *Perspect. Plant Ecol. Evol. Syst.* **17**, 66–77 (2015).

44. Hirose, T. & Werger, M. J. A. Photosynthetic capacity and nitrogen partitioning among species in the canopy of a herbaceous plant community. *Oecologia* **100**, 203–212 (1994).

45. Hoffmann, K. & Bivour, W. Klimauntersuchungen in Jena für die Anpassung an den Klimawandel und seine erwarteten Folgen. *Berichte des Deutschen Wetterdienstes* **243** (2014).

46. Roscher, C. et al. The role of biodiversity for element cycling and trophic interactions: an experimental approach in a grassland community. *Basic Appl. Ecol.* **5**, 107–121 (2004).

47. Petermann, J. S. et al. Biology, chance, or history? The predictable reassembly of temperate grassland communities. *Ecology* **91**, 408–421 (2010).

48. Volkmann, T. H. M., Haberer, K., Gessler, A. & Weiler, M. High-resolution isotope measurements resolve rapid ecohydrological dynamics at the soil-plant interface. *New Phytol.* doi:10.1111/nph.13868 (2016).

49. Monsi, M. & Saeki, T. The light factor in plant communities and its significance for dry matter production. *Japanese J. Bot.* **14**, 22–52 (1953).

50. Villéger, S., Mason, N. W. H. & Mouillot, D. New multidimensionale functional diversity indices for a multifaceted framework in functional ecology. *Ecology* **89**, 2290–2301 (2008).

51. Team, R. C. R: A Language and Environment for Statistical Computing. (R Foundation for Statistical Computing, 2016).

52. Calcagno, V. The glmulti package. Model selection and multimodel inference made easy. R project for statistical computing. 1–20 (2015).

53. Burnham, K. & Anderson, D. R. *Model Selection and Multimodel Inference*. **40** (2001).

## Acknowledgements

For equipment and field measurements this study benefited from the CNRS human (Damien Landais, Sebastien Devidal, Clement Piel and Olivier Ravel) and technical resources allocated to the ECOTRONS Research Infrastructure as well as from the state allocation 'Investissement d'Avenir' ANR-11-INBS-0001. We also thank Anne Ebeling for organizational support and Annette Gockele for help in the field.
**Author Contributions**

A.M. and J.R. designed the study with substantial input from A.G., C.R. and M.S.-L. A.M. analyzed the data and wrote the manuscript. All authors (A.M., A.G., C.R., L.R., Z.K., D.B., K.P.-W., S.Z., L.G., T.B., M.S.-L. and J.R.) contributed to field work, data collection and provided input on the manuscript writing.

**Additional Information**

**Supplementary information** accompanies this paper at doi:10.1038/s41598-017-08819-9

**Competing Interests:** The authors declare that they have no competing interests.

**Publisher's note:** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_6.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) 2017