Thank you for visiting nature.com. You are using a browser version with limited support for CSS. To obtain the best experience, we recommend you use a more up to date browser (or turn off compatibility mode in Internet Explorer). In the meantime, to ensure continued support, we are displaying the site without styles and JavaScript.
Enrichment of nutrients and loss of herbivores are assumed to cause a loss of plant diversity in grassland ecosystems because they increase plant cover, which leads to a decrease of light in the understory1,2,3. Empirical tests of the role of competition for light in natural systems are based on indirect evidence, and have been a topic of debate for the last 40 years. Here we show that experimentally restoring light to understory plants in a natural grassland mitigates the loss of plant diversity that is caused by either nutrient enrichment or the absence of mammalian herbivores. The initial effect of light addition on restoring diversity under fertilization was transitory and outweighed by the greater effect of herbivory on light levels, indicating that herbivory is a major factor that controls diversity, partly through light. Our results provide direct experimental evidence, in a natural system, that competition for light is a key mechanism that contributes to the loss of biodiversity after cessation of mammalian herbivory. Our findings also show that the effects of herbivores can outpace the effects of fertilization on competition for light. Management practices that target maintaining grazing by native or domestic herbivores could therefore have applications in protecting biodiversity in grassland ecosystems, because they alleviate competition for light in the understory.
Anthropogenic nutrient enrichment and human-induced shifts in herbivore pressure are two major drivers of global change that markedly alter the composition, diversity and functioning of terrestrial plant communities3,4,5. Increased nutrient supply due to eutrophication has been connected to a strong and persistent loss of diversity6,7, and herbivory by grazing mammals is known to be one of the key factors that maintains diversity2,8,9, as described already by Charles Darwin in his On the Origin of Species in 1859: “No wonder that, as soon as the land was enclosed, it became thickly clothed…”10. The assumed central mechanism that drives these responses and links the effects of nutrients and herbivory on diversity is competition for light1,2,3,11. Fertilization disproportionally promotes the growth of taller plants with more canopy cover and better access to light, at the expense of shorter plants and seedlings in the understory, and this leads to reduced diversity12,13. The asymmetry in competition for light stems from the fact that light as a resource is unidirectional13,14. By contrast, by consuming vegetation and selectively targeting taller species that are superior in the competition for light, herbivores can directly reduce canopy cover and increase the availability of light for shorter plants1,15. Herbivory therefore has the potential to alleviate competition for light and maintain diversity2,3.
Testing the role of light in maintaining diversity requires direct experimental tests in which light—the limiting resource—is added into the understory of plant communities where competition for light is strongest. However, most studies use different indirect ways to address the role of competition for light, such as unmanipulated light measurements, tiebacks or neighbour removals3,9,13,16, which, for methodological reasons, may produce misleading results. For example, neighbour removal could release nutrients from roots and alter temperature and humidity, and might not reflect light as a causal factor. The strongest direct evidence so far comes from a controlled greenhouse experiment, which showed that adding light by lamps to the understory of plant communities prevented the loss of species that otherwise resulted from fertilization12. However, these results have not been confirmed in natural field conditions, with more complex communities and with herbivores that are predicted to interact with nutrient effects on competition for light17,18.
Here, we experimentally manipulated light, herbivory and nutrient supply, in a full-factorial design, to test the direct causal role of light limitation in driving plant diversity loss from eutrophication and loss of herbivory in a natural, species-rich grassland. We installed modern light-emitting diode (LED) lamps below the plant canopy to provide light to plants in the understory, where light should be needed most. We combined the light-addition treatment with fertilization and exclosures that prevented grazing by sheep to test the responses of plant community richness and diversity (Fig. 1a,c and Methods). Our LED lamps, with a spectrum mimicking natural sunlight, increased the availability of light in the understory compared to ambient levels (a 57% increase, on average, in the quantity of light in fertilized and fenced plots by light addition; Fig. 1b). Furthermore, our light-addition treatment had no detectable effects on humidity and air temperature near the soil surface (Fig. 1b and Extended Data Table 1), and caused minimal disturbance (Methods). Our experiment, which we term ‘eDiValo’ (ecological effects of light (‘valo’ in Finnish) on diversity) was conducted in grazed pastures at the Global Change Experimental Facility (GCEF) in Bad Lauchstädt, Germany, where species-rich natural grassland vegetation (approximately 23 vascular plant species in a 0.25-m2 area) was exposed to short-time high-intensity sheep-grazing events two to three times each growing season. We measured species richness (because richness should be sensitive to rare species becoming extinct) and Shannon diversity (because diversity should be sensitive to species becoming less abundant before they become extinct), total live plant and litter cover and plant functional traits in 0.5 × 0.5-m plots.
a, Experimental design and predictions of the effects of treatments. See Methods for more details on the NPK fertilizer, herbivore exclosures, and lamps used for light addition. b, Effects of treatments on light quantity, relative humidity and air temperature. Values are mean ± s.e.m. derived from linear mixed-effects (LME) models, in which parameter significance was assessed by F-tests (two-tailed). For sample sizes and statistics, see Extended Data Table 1. Faded dots show individual data points. C, control; F, fertilization; E, exclosure; EF, exclosure and fertilization. c, Photos of eDiValo experimental design (left), lamps bringing light to small understory plants (middle) and a flock of sheep grazing in an experimental block (right). Photo credit: A.E. Images in a were created by G. Rada (iDiv, Media and Communications).
First, in 2017, we tested whether adding light could offset the negative effect of fertilization on richness and diversity. During the 2017 growing season, we prevented grazing in the whole experimental area to allow fertilization effects to develop. Fertilization rapidly decreased species richness and species diversity, measured as Shannon diversity, in 2017 (by 9.3% and 7.3%, respectively; Fig. 2a,b and Extended Data Tables 2 and 3). Notably, restoring light to plants in the understory of the fertilization treatment offset the loss of diversity (fertilization × light addition interaction on Shannon diversity; Fig. 2b and Extended Data Table 2). Although light addition reduced the loss of species richness (Fig. 2a and Extended Data Table 2), the interaction between light addition and fertilization was weaker on species richness than on Shannon diversity, suggesting that light addition promoted more equal abundances, whereas richness—through either extinctions or gains—was slower to respond. Our findings confirm the mechanism that was reported in a previous controlled greenhouse experiment, in which added light similarly rapidly reversed the loss of species resulting from nutrient addition12. This result is consistent with the hypothesis of competition and diversity in herbaceous vegetation1, the theory of resource competition15 and quantitative models that address the asymmetric nature of competition for light13,14. Our study provides the first—to our knowledge—direct experimental demonstration that competition for light is a central mechanism that leads to the loss of plant diversity in conditions of nutrient enrichment, in species-rich real-world grassland communities that experience varying environmental conditions and complexities of trophic interactions.
a–d, Effects of fertilization and light addition on species richness (a) and Shannon diversity (b) in 2017 and on species richness (c) and Shannon diversity (d) in 2019. Values are mean ± s.e.m. derived from LME models, in which parameter significance was assessed by F-tests (two-tailed). In 2017, n = 20 for quadrats from which species richness and Shannon diversity were estimated; in 2019: fertilized, unlighted and lighted, n = 22; unfertilized, unlighted and lighted, n = 18. For statistics, see Methods and Extended Data Tables 2 and 5. In 2017, herbivore exclusion treatment was not yet applied. In 2019, data were pooled across herbivore exclusion treatment to show fertilization effects more clearly; for a full presentation of treatments, see Extended Data Fig. 2. These results are from our eDiValo field experiment; see Fig. 1 and Methods.
Second, after the first experimental year, we extended our experimental design to include herbivory by sheep, and erected fenced herbivore exclosures randomly around half of the plots that were established in 2017 (Fig. 1a). We simultaneously removed the temporary fence around the whole experimental area and allowed the grazing of sheep in unfenced plots. In 2019, we tested whether light addition could offset the negative effect of herbivore exclusion on richness and diversity and whether light, nutrients and herbivory interact. Herbivore exclusion decreased species richness by 12.5% and Shannon diversity by 11.7%, independent of fertilization (Fig. 3a,b and Extended Data Tables 2 and 4) and consistent with previous studies3,8. At the same time, herbivore exclusion increased total cover (Fig. 3c) and decreased light (Fig. 1b). These results indicate that herbivory is a dominant factor controlling light availability and plant diversity. Our key finding was that experimental addition of light in the understory mitigated the loss of richness and diversity due to the herbivore exclusion (exclosure × light interactions on richness and Shannon diversity; Fig. 3a,b, Extended Data Tables 2 and 4 and Extended Data Fig. 2b,c). This provides experimental evidence that herbivores maintain diversity by alleviating competition for light. Extirpation of large mammalian herbivores5,19 may therefore contribute to diversity loss in plant communities through increased competition for light.
a–d, The effects of fencing and light addition on species richness (a), Shannon diversity (b), total cover (c) and litter cover (d). Both total and litter cover were visually estimated from 0.5 × 0.5-m plots (see Methods). Values are mean ± s.e.m. derived from LME models, in which parameter significance was assessed by F-tests (two-tailed); n = 20 for quadrats from which species richness, Shannon diversity, total cover and litter cover were estimated. For statistics, see Methods and Extended Data Tables 2 and 5. For a full presentation of the treatments, see Extended Data Figs. 2 and 3. These results are from our eDiValo field experiment; see Fig. 1 and Methods.
Fertilization had a negative main effect on both richness and diversity in 2019 (leading to decreases of 12.7% and 8.2%, respectively, similar magnitudes to 2017; Extended Data Tables 2 and 4 and Fig. 2c,d); however, light-addition effects in fertilized plots compared to 2017 were transitory and did not restore diversity in 2019 (no significant fertilization × light interaction; Extended Data Tables 2 and 4), in contrast to the effects of light addition in herbivore exclosures (Fig. 3a,b and Extended Data Fig. 2a,b). The negative main effect of fertilization on plant diversity was independent of changes in total cover and light availability, which were not affected by fertilization (Figs. 1b and 3c and Extended Data Fig. 2c). As a result, addition of light in fertilized plots did not restore diversity. By contrast, in the continued absence of herbivory, control of light competition shifted from bottom-up effects of fertilization to top-down control by consumers that regulated total cover, light availability, competition for light and diversity. Our results suggest that herbivore-mediated processes linked to succession can outpace the effects of fertilization on competition for light. However, mechanisms other than competition for light could explain the negative main effect of fertilization on diversity. For example, litter cover was higher in fertilized than in unfertilized plots (64.6% versus 44.9%, respectively), which could have contributed to this negative effect (main effect of fertilization; Extended Data Tables 2 and 4 and Extended Data Fig. 3b). As predicted by theory15, supplying multiple nutrients could also have reduced diversity without increasing total cover and light limitation by reducing the dimensionality of belowground nutrient trade-offs20. Fertilization could also exert direct toxic effects on seedling richness, reducing diversity through decreased recruitment opportunities21. Our results are consistent with these previous findings and suggest that herbivores can outpace the adverse effects of fertilization on diversity when they occur through intensifying competition for light and affecting litter, but not when they occur through filling nutrient niches in the soil or by direct toxic effects on seedlings.
Adding light, however, did not completely offset the negative effect of herbivore exclusion on Shannon diversity. One explanation is that additional light did not fully alleviate light limitation in the understory inside herbivore exclosures (Fig. 1b). Alternatively, besides competition for light, other factors could also have contributed to the loss of diversity inside herbivore exclosures. In our experimental plots, herbivore exclusion more than doubled the amount of undecomposed plant litter (32.8% cover outside exclosures versus 78.7% cover inside exclosures, Fig. 3d), which can decrease diversity through reduced opportunities for plant recruitment from seed22. Several other factors can also change as a result of herbivore exclusion and affect diversity23. We showed that herbivores exhibited the strongest control on litter and total cover (Fig. 3c,d and Extended Data Fig. 3b), which modified the abiotic environment, leading to a lower air temperature (2.5 °C and 12.2% lower on average) and higher levels of humidity (35% higher) inside exclosures than in grazed plots (Fig. 1b). In addition, increasing the levels of light slightly reduced the total cover inside unfertilized exclosures (Fig. 3c and Extended Data Fig. 2c). This result is counterintuitive, but could be due to changes in species composition, more-even distribution of species in the canopy and/or reduced cover by dominant species in this treatment combination. Overall, our results show that changes in understory plant life in ungrazed conditions include a suite of factors, all of which may be important for plant performance23,24.
Finally, we assessed which plant functional traits made species more or less responsive to light addition under fertilization and herbivore exclusion. Short species (size-related trait) and species with low light interception (trait related to conservative resource-use; a small specific leaf area, SLA) should be at a competitive disadvantage under low light conditions25,26,27, and should therefore benefit the most from light addition. Consistent with this prediction, species with a lower SLA and lower leaf water content (LWC) had a higher probability to increase from unlighted to lighted subplots; that is, were more likely to benefit from light addition (Fig. 4 and Extended Data Fig. 4). By contrast, short species did not benefit from light addition (see Methods). It is possible that SLA and LWC better integrate plant responses to light during the whole life cycle, including seedling and juvenile stages, even if height is more important for the competitive ability of adult plants24,28.
a,b, Probability that individual species increase in response to light addition as a function of specific leaf area (SLA; a) and leaf water content (LWC; b) in 2019. Fitted lines are from generalized linear mixed-effects (GLME) models with a binomial error structure, in which parameter significance was assessed by χ2-tests (two-tailed; see Methods): a, exclosure χ2 = 0.45, P = 0.50; SLA χ2 = 52.24, P < 0.0001; exclosure × SLA χ2 = 8.08, P = 0.0045, b, LWC χ2 = 87.18, P < 0.0001. Shaded areas represent 95% confidence intervals; n = 1,580 species observations for fenced and unfenced SLA and n = 3,240 species observations for LWC. In the models, species were nested within plots (40) that were nested within blocks (10). Fenced, no exclosure; Unfenced, exclosure. Species with a lower SLA benefit more from light addition in unfenced plots compared to fenced plots (exclosure × SLA interaction). These results are from our eDiValo field experiment; see Fig. 1 and Methods.
Anthropogenic nutrient enrichment and changes in the grazing practices of livestock and native herbivores are among the most influential global-change factors that can decrease or rescue the biodiversity of global grasslands3,6,9,29. Understanding the ecological mechanisms by which nutrients and consumers operate is fundamental for our ability to maintain and manage biodiversity, and to develop better conservation actions in the Anthropocene. By supplying light to understory plants, our experimental results here provide direct evidence of the role of competition for light as a mechanism that decreases diversity under fertilization and the absence of mammalian herbivory in natural grasslands. Our findings especially highlight the importance of grazing by native and domestic herbivores as factors that foster higher biodiversity. Carefully planned livestock management practices, protection of native herbivores and re-establishing large mammalian herbivores to areas from where they have been extirpated (that is, rewilding) may therefore be key conservation strategies in the Anthropocene30,31, because they could promote plant diversity by alleviating competition for light.
Our study site is located at the Bad Lauchstädt Field Research Station, Bad Lauchstädt, Germany (51° 22060 N, 11° 50060 E), which belongs to the Helmholtz Centre for Environmental Research–UFZ. Long-term mean annual precipitation in the area is 489 mm and the mean annual temperature is 8.9 °C (ref. 32). During 2018 and 2019, Europe experienced a record-setting drought that was especially severe in 2018 (refs. 33,34); the mean annual precipitation at our study site in 2018 and 2019 was 254 mm and 353 mm, respectively, whereas 2017 was a more normal year, with a mean annual precipitation of 403 mm. Mean annual temperatures were above average: 2017, 10.5 °C; 2018, 10.8 °C; 2019, 11.2 °C (data from the weather station at the Bad Lauchstädt field station). The soils in the study area are fertile Haplic Chernozem type32,35.
Our eDiValo experiment was conducted in the GCEF, which was designed to investigate climate change effects under different land-use scenarios32. We used 10 ‘extensively’ used pastures of the GCEF in our experiment; that is, 384-m2 (16 × 24 m) areas of grassland (hereafter called ‘pastures’) that were grazed by a flock of 20 sheep 2–3 times each year. Grazing was implemented as short-time high-intensity grazing events, each lasting 24 h (ref. 32). This type of high-intensity but short-term grazing is considered better in maintaining species richness as it gives plants more time to recover between grazing events36. It is also a recommended management type for nature conservation areas in Germany37. Vegetation in the pastures was species-rich grassland vegetation that is typical of drier regions of central Germany32,38. The whole GCEF was fenced to exclude native large mammalian herbivores (for example, deer); however, European hare (Lepus europaeus), wood mice (Apodemus sylvaticus) and voles (Microtus arvalis) are common at the site.
Our experimental design was originally intended to test the dependence of light competition on nutrient and herbivory under current and future climatic scenarios. Although we included both climate treatments in our data, climate was never significant for richness and Shannon diversity, either alone or in interaction with other factors, and our focus was therefore on the other treatments. Five of the above random pastures received future climatic treatment which was based on different dynamic regional climate models for Germany, all predicting an increased mean temperature by approximately 2 °C year-round, strongly decreased summer precipitation and slightly increased spring and autumn precipitation (https://www.regionaler-klimaatlas.de/) (ref. 32). Passive night-time (after sunset and before sunrise) warming through the use of roller blinds attached to the GCEF roof and eastern and western wall structures was used to increase the air temperature. In each spring (1 March–31 May) and autumn (1 September–30 November), future climate plots received 110% of the ambient rainfall and in the summer (1 June–31 August), they received 80% of the ambient rainfall. The precipitation treatment was adjusted weekly and compensated for a possible night-time reduction in rainfall due to temperature treatment. A detailed description of the future climate treatment is provided in a previous report32.
We first tested whether adding light can offset the negative effect of fertilization on plant diversity. In May 2017, we established a full-factorial experiment of fertilization and light addition. Within each 10 pastures (5 in ambient climatic conditions, 5 in future climatic conditions), we established 4 plots of 1.4 × 1.4 m, separated by a 1-m buffer zone (hereafter called ‘blocks’), in total 40 plots and 10 blocks. At the time the experiment was established, vegetation in the whole experimental area (that is, in a block of 4 plots and the surrounding 1-m area) was trimmed to a height of 5 cm to make conditions uniform and the whole area was temporarily fenced to let the experiment establish and fertilization effects develop. The temporary fence was removed in August when the herbivore exclusion treatment was started. Therefore, there was no grazing by sheep in the experimental plots in the summer of 2017. Two randomly chosen plots received fertilizer treatment and two were controls. For the former (fertilizer-treatment plots), slow-release granular NPK fertilizer (a mixture of Haifa Multicote 2 M 40-0-0 40% N; Triple Super Phosphate (TSP) 45% P205; and potassium sulfate fertilizer 50% K2O, 45% SO3) was added twice per growing season, in a total of 10 g N, 10 g P and 10 g K per m² (see ref. 3 for a similar protocol that is used in grasslands worldwide). In 2017, the first fertilization was done at the beginning of June right after establishing the experiment and the second fertilization was done at the beginning of July. In the subsequent years, the first fertilization was done at the beginning of the growing season (late March–April) and the second fertilization was done in June. In 2019, two previously unfertilized plots were accidentally fertilized and were thereafter treated as fertilized plots.
To manipulate light, 1.4 × 1.4-m plots were further divided into two subplots, 0.7 m × 1.4 m each, and one of these was randomly assigned to the light-addition treatment, resulting in 80 subplots (Fig. 1). We installed two 120-cm-long and 3.5-cm-wide recently developed LED lamps (C65, Valoya) parallel to each other and at a 28-cm distance from each other to each light-addition subplot. To increase light for the small understory plants that are the most likely to suffer from competition for light, we installed the lamps 10 cm above the smallest plants. The lamps were gradually uplifted over the course of the growing season to follow the growth of the smallest plants. As our light-addition treatment was intended to mimic natural sunlight (that is, making a gap in a dense vegetation and allowing the sunshine in), we chose the spectrum of the lamps to include all wavelengths of sunlight, including small amounts of ultraviolet and infrared. Each lamp added roughly 350–400 µmol and did not alter the air or aboveground soil surface temperature (Fig. 1b), which is an improvement on previous studies12. Each year, we added light during the active growing season: the lamps were switched on early in the spring (March–April), when temperatures were clearly above zero, and switched off and removed when temperatures dropped close to zero in November–December and aboveground plant parts had died and formed litter. Each day, the lamps were set to switch on two hours after sunrise, and to switch off two hours before sunset, and when the temperature exceeded 28 °C to prevent overheating. We did not install unpowered lamps to unlighted plots because our modern, narrow LED lamps caused minimal disturbance (see below) and no heating (Fig. 1b), and because unpowered lamps would have added an artefact in that they create shade that does not occur when the lamps are on in lighted plots.
At the end of August 2017, after running the fertilization–light-addition experiment for one growing season, we expanded the experiment by implementing the herbivore exclusion treatment in a full-factorial combination with the other treatments. Two of the previously established 1.4 m × 1.4-m plots, one with and one without the fertilization treatment, were randomly allotted to the herbivore (sheep) exclusion treatment and fenced with rectangular metal fences of 1.8 m × 1.8 m, 82 cm height and 10 cm mesh size. At the same time, the temporary fence established in May 2017 was removed from around the whole experimental area, allowing the grazing of sheep in unfenced plots. The fences did not exclude mice, voles and hares. For the time of each grazing event, lamps in grazed subplots were removed and switched off in the ungrazed subplots. Uplifting the lamps from grazed plots did not cause disturbance because vegetation in grazed plots was always short and did not reach above the lamps. Inside exclosures, lamps were always kept in place during the growing season, and plants could freely grow around and above them.
In July 2017, we established 50 cm × 50-cm permanent quadrats in every subplot for plant community sampling. We visually estimated the per cent areal cover for all species occurring in the quadrats, and litter cover, from the beginning of June to mid-June 2019, when the vegetation was at its peak biomass. The 2017 sampling happened later, in mid-July, because vegetation in all plots and surrounding areas was trimmed to a height of 5 cm at the time of the establishment of the experiment at the end of May, and it took later for vegetation to reach its peak biomass. In 2018, the effects of drought were devastating, and most plants had senesced or died before the planned sampling date; we therefore omitted the year 2018. At the beginning of each growing season—that is, when the lamps were installed and switched on—there was very little live biomass in the plots, and the maximum height of existing plants was approximately 5 cm (in all plots). During the peak biomass the maximum plant height was up to approximately 1 m; however, it varied greatly between the treatments and was especially low in grazed plots. All vegetation surveys were done by the same trained and experienced person with a minimum estimate threshold of 0.1%. We used plant cover data to calculate species richness and Shannon diversity.
In May–June 2020, we measured plant height (centimetres), SLA (leaf area in square millimetres per milligram of dry mass), foliar C:N (based on the per cent C and N in plant leaves) and LWC (leaf water content as 1,000 − LDMC (the ratio of leaf dry mass to saturated fresh mass), expressed as milligrams per gram39) for most species occurring in the experimental plots, and complemented the trait data from the TRY Plant Trait Database40,41,42 (v.5.0; https://www.try-db.org/TryWeb/Home.php) and for one species one trait value from another source9. The trait data were collected from seven to ten individuals per species from the study site or close areas; the collection and handling followed standard protocols39. We chose these traits because they are widely documented to be associated with responsiveness to soil nutrients, herbivory and light9,26,27,43,44,45,46. We used all traits as, although they partially reflect similar ecological adaptations (for example, leaf economics spectrum43), they could also potentially reflect independent and distinctive processes, and differently mediate the responses of species to our treatments. For example, SLA and LWC in our dataset correlated weakly (r2 = 0.16), but were to a greater extent uncorrelated (Extended Data Table 6), and could function differently, for example, in light capture and drought tolerance26,39. In 2017, our trait data covered on average 97.7–98.6% of the total cover in the plots, the value slightly differing depending on the trait as we did not have all traits for all species. Our own trait collections covered on average 96.6–97.6% and TRY data covered on average 0.9–2% of the total cover. In 2019, the whole trait data covered on average 99.5% of the total cover in the plots, again slightly depending on the trait. Our own trait collections covered on average 94.2–96.5% and TRY data covered on average 2.7–5.3% of the total cover.
We measured several soil and other environmental properties from the experimental plots. Light availability (photosynthetically active radiation; PAR) in unlighted and lighted (under lamps) subplots was measured using LI-190R and LI-250A meters (LI-COR), approximately 7–10 cm under the lamps and 15–20 cm above ground level. We measured light availability from the same distance to the ground in unlighted plots. Measurements of light availability were done in mid-July 2020 on three consecutive cloudless days around noon. Note that in grazed plots, light levels between lighted and unlighted plots are more similar than inside exclosures (Fig. 1), because herbivores keep the vegetation short, and natural sunlight can therefore reach under the lamps where the light measurements were taken. Air temperature and humidity were recorded from unlighted and lighted (under lamps) subplots using loggers (HOBO MX2301A, Onset Computer Cooperation) that were installed approximately 7 cm under the lamps and to the same height from the ground in unlighted plots, and were replicated under different combinations of fertilization, herbivore exclusion and light addition in ambient climatic conditions three times (n = 3). The logger data were collected in May 2019 before the effects of drought were visible.
We analysed our data in two steps. First, to test whether competition for light mediates the effect of fertilization on diversity, we analysed the effects of fertilization and light and their interaction on species richness and Shannon diversity using data from 2017, when the herbivore exclusion treatment had not yet been implemented. We also analysed the effects of treatment on total vegetation cover and litter cover. We fit LME models in which diversity (species richness and Shannon diversity), total cover and litter cover, each in their own model, were explained by fertilization, light addition and their interaction (fixed variables). All treatments were categorical variables with two levels (treated and untreated). In each model, subplot was nested within plot, which was nested within block (nested random variable). We simplified the models using the anova() function for model comparison in the nlme and lme4 packages in R (ref. 47) (on the basis of log likelihood ratio tests; P ≥ 0.05; Extended Data Table 2). This was done to uncover the significance of the main effects and interaction terms, to avoid overparametrization47,48 and to provide model-derived parameter estimates for the figures (Extended Data Table 5). However, we also provide full model results that are qualitatively similar to the results of simplified models (Extended Data Tables 3 and 4); therefore, model choice did not affect our conclusions. Climate treatment was included in all original models but was never significant for richness and diversity, and was not considered further. Total cover and litter results for 2017 are reported in Extended Data Figs. 1a,b and 3a). As there was heterogeneity in the variance structure between treatments, we used the varIdent() function in the nlme package in R to allow each treatment combination to have a different variance. Model fit was inspected using model diagnostic plots in the package nlme. In the full design with climate included, the number of replicates per treatment combination was ten.
Second, to include herbivore exclusion to the experimental design and to test whether competition for light mediates the effect of herbivore exclusion on diversity, and whether competition for light, herbivory and fertilization interact, we analysed the effects of herbivore exclusion, fertilization, light and their interactions on species richness and Shannon diversity using data from 2019. All treatments were categorical variables with two levels (treated and untreated). We also analysed the effects of treatment on total vegetation cover and litter cover. We fit similar models to those described above, except that herbivore exclusion was an additional fixed factor in the models. We simplified the models, used the varIdent() function to account for heteroscedasticity and checked the model fit using model diagnostic plots, as above. Climate treatment was included in all original models but was significant for litter cover only, and was not considered further. In the full design with climate included, the number of replicates per treatment combination was five.
To further assess which plant traits increased the probability of species benefiting from the addition of light, we first created a binary response variable: those species that increased from unlighted to lighted plots (that is, had a higher value in a lighted than an unlighted plot) were given a value of 1 and those that did not were given a value of 0. This response variable takes into account rare species that emerged or persisted in the lighted plots but were absent in the unlighted plots (that is, species gains and losses) and changes in small, subordinate species (those that are likely to benefit from light addition) with small but consistently trait-dependent changes in response to light. It is also in line with our species richness analyses, as species gains and losses ultimately determine richness responses. We did not use different indexes (for example, lnRR or RII) because these could not handle multiple zero values and species losses or gains (that is, species having zero cover in either unlighted or lighted subplots). Second, we fit GLME models with a binomial error structure (family = “binomial”, link = “logit”) in which a probability of a species increasing from unlighted to lighted plots was explained by categorical experimental treatments (fertilization, herbivore exclusion and their interactions), traits (SLA, height, LWC, foliar C:N), and interactions between the treatments and traits. Each trait was analysed in its own model as some of the traits were correlated (Extended Data Table 6), and to avoid overly complex models and overparametrization47,48. We included all species for which we had traits in the models. As we calculated the increase in cover from unlighted to lighted plots, our smallest experimental unit in trait analyses was a plot (not a subplot, unlike in other analyses). As there were several species in the same plots, we nested species within plots, and plots within blocks. We similarly simplified the models to include only significant variables (on the basis of χ2 tests; P ≥ 0.05). We did not include a crossed random effect for species in the models because the full models with a more complex random structure did not converge; however, when we refitted the simplified models with a crossed random effect for species, we found that the models converged (with scaled data) and that the significance of the effects remained qualitatively the same. Climate was included in all original models but was never significant. In addition, C:N and height did not predict the responsiveness of species to light in either year (P ≥ 0.13 for both); results are therefore not shown. In the full design with climate included, the number of replicates per treatment combination was five; however, the number of observations was greater (see Fig. 4 and Extended Data Fig. 4). To make sure that our results for SLA and LWC were not influenced by whether they were analysed in separate models or in the same model, or by the order in which they were in the models, we also performed analyses in which both SLA and LWC were included (in both orders). Results remained qualitatively similar and are not discussed further.
Furthermore, to check whether our trait results were driven primarily by species gains and losses or changes in abundance, we ran additional trait analyses for which we calculated the change in cover between lighted and unlighted subplots (cover in lighted subplot − cover in unlighted subplot), and analysed the ‘change’ with otherwise similar trait models to those described above, except that we used Gaussian error structure. With this index, which gives a disproportionate importance to the abundant species, we found that traits were poor predictors of changes in cover between lighted and unlighted plots (all interactions were non-significant, P > 0.05, except for a marginally significant C:N × fertilization interaction in 2017 that was no longer visible in 2019; results not shown; codes and data available in the Dryad repository). We also analysed presence–absence-based species losses and gains. In these models, each species was given a value of 1 when it was present in the lighted subplot but absent from the unlighted subplot; otherwise, these models were similar to the binomial trait models described above. These models produced, to a large extent, similar results to our models using the probability of increase in response to light as a response variable (results not shown; codes and data available in the Dryad repository). These additional analyses and results support using the probability of increase in response to light as our response variable, rather than abundance-based metrics, as it includes both gains and losses and abundance aspects, and is therefore a general test that is well suited to assessing species gains and extinctions and changes in subordinate species.
All statistical analyses were performed using R v. 4.0.0 (ref. 49). We used the nlme package (v.3.1.147) for LME models50, the lme4 package (v.1.1.23) for GLME models51, and the car package52 for P values (v.3.07).
Further information on research design is available in the Nature Research Reporting Summary linked to this article.
The datasets generated and analysed during this study are available in the Dryad repository: https://doi.org/10.5061/dryad.rjdfn2zdm.
The R code and scripts used in this study are available in the Dryad repository: https://doi.org/10.5061/dryad.rjdfn2zdm.
Newman, E. I. Competition and diversity in herbaceous vegetation. Nature 244, 310 (1973).
Olff, H. & Ritchie, M. E. Effects of herbivores on grassland plant diversity. Trends Ecol. Evol. 13, 261–265 (1998).
Article CAS PubMed Google Scholar
Borer, E. T. et al. Herbivores and nutrients control grassland plant diversity via light limitation. Nature 508, 517–520 (2014).
Article ADS CAS PubMed Google Scholar
Clark, C. M. & Tilman, D. Loss of plant species after chronic low-level nitrogen deposition to prairie grasslands. Nature 451, 712–715 (2008).
Article ADS CAS PubMed Google Scholar
Estes, J. A. et al. Trophic downgrading of planet Earth. Science 333, 301–306 (2011).
Article ADS CAS PubMed Google Scholar
Bobbink, R. et al. Global assessment of nitrogen deposition effects on terrestrial plant diversity: a synthesis. Ecol. Appl. 20, 30–59 (2010).
Article CAS PubMed Google Scholar
Midolo, G. et al. Impacts of nitrogen addition on plant species richness and abundance: a global meta‐analysis. Global Ecol. Biogeogr. 28, 398–413 (2019).
Collins, S. L., Knapp, A. K., Briggs, J. M., Blair, J. M. & Steinauer, E. M. Modulation of diversity by grazing and mowing in native tallgrass prairie. Science 280, 745–747 (1998).
Article ADS CAS PubMed Google Scholar
Kaarlejärvi, E., Eskelinen, A. & Olofsson, J. Herbivores rescue diversity in warming tundra by modulating trait-dependent species losses and gains. Nat. Comm. 8, 419 (2017).
Darwin, C. R. On the Origin of Species by Means of Natural Selection (John Murray, 1859).
Holt, R. D., Grover, J. & Tilman, D. Simple rules for interspecific dominance in systems with exploitative and apparent competition. Am. Nat. 144, 741–771 (2004).
Hautier, Y., Niklaus, P. A. & Hector, A. Competition for light causes plant biodiversity loss after eutrophication. Science 324, 636–638 (2009).
Article ADS CAS PubMed Google Scholar
DeMalach, N., Zaady, E. & Kadmon, R. Light asymmetry explains the effect of nutrient enrichment on grassland diversity. Ecol. Lett. 20, 60–69 (2017).
Weiner, J. Asymmetric competition in plant populations. Trends Ecol. Evol. 5, 360–364 (1990).
Article CAS PubMed Google Scholar
Tilman, D. Resource Competition and Community Structure (Princeton University Press, 1982).
Dickson, T. L. & Foster, B. L. Fertilization decreases plant biodiversity even when light is not limiting. Ecol. Lett. 14, 380–388 (2011).
Huisman, J., Grover, J. P., Van der Wal, R. & Van Andel, J. in Herbivores: Between Plants and Predators (eds Olff, H., Brown, V. K. & Drent, R. H.) 239–269 (Blackwell Science, 1999).
Bakker, E. S., Ritchie, M. E., Olff, H., Milchunas, D. G. & Knops, J. M. H. Herbivore impact on grassland plant diversity depends on habitat productivity and herbivore size. Ecol. Lett. 9, 780–788 (2006).
Atwood, T. B. et al. Herbivores at the highest risk of extinction among mammals, birds, and reptiles. Sci. Adv. 6, eabb8458 (2020).
Article ADS PubMed PubMed Central Google Scholar
Harpole, W. S. et al. Addition of multiple limiting resources reduces grassland diversity. Nature 537, 93–96 (2016).
Article ADS CAS PubMed Google Scholar
Britto, D. T. & Kronzucker, H. J. NH4+ toxicity in higher plants: a critical review. J. Plant Physiol. 159, 567–584 (2002).
Foster, B. L. & Gross, K. L. Species richness in a successional grassland: effects of nitrogen enrichment and plant litter. Ecology 79, 2593–2602 (1998).
Harpole, W. S. et al. Out of the shadows: multiple nutrient limitations drive relationships among biomass, light and plant diversity. Funct. Ecol. 31, 1839–1846 (2017).
Valladares, F., Laanisto, L., Niinemets, U. & Zavala, M. A. Shedding light on shade: ecological perspectives of understorey plant life. Plant Ecol. Div. 9, 237–251 (2016).
Westoby, M., Falster, D. S., Moles, A. T., Vesk, P. A. & Wright, I. J. Plant ecological strategies: some leading dimensions of variation between species. Annu. Rev. Ecol. Syst. 33, 125–159 (2002).
Valladares, F. & Niinemets, U. Shade tolerance, a key plant feature of complex nature and consequences. Annu. Rev. Ecol. Evol. Syst. 29, 237–257 (2008).
Reich, P. B. The world-wide ‘fast-slow’ plant economics spectrum: a traits manifesto. J. Ecol. 102, 275–301 (2014).
Craine, J. M. & Dybzinski, R. Mechanisms of plant competition for nutrients, water and light. Funct. Ecol. 27, 833–840 (2013).
Filazzola, A. et al. The effects of livestock grazing on biodiversity are multitrophic: a meta-analysis. Ecol. Lett. 23, 1298–1309 (2020).
Guyton, J. A. et al. Trophic rewilding revives biotic resistance to shrub invasion. Nat. Ecol. Evol. 4, 712–724 (2020).
Svenning, J.-C. Rewilding should be central to global restoration efforts. One Earth 3, 657–660 (2020).
Schädler, M. et al. Investigating the consequences of climate change under different land-use regimes: a novel experimental infrastructure. Ecosphere 10, e02635 (2019).
Hari, V., Racovec, O., Markonis, M. & Kumar, R. Increased future occurrences of the exceptional 2018–2019 Central European drought under global warming. Sci Rep. 10, 12207 (2020).
Article ADS CAS PubMed PubMed Central Google Scholar
Buras, A., Rammig, A. & Zang, C. S. Quantifying impacts of the 2018 drought on European ecosystems in comparison to 2003. Biogeosciences 17, 1655–1672 (2020).
Altermann, M. et al. Chernozem—soil of the year 2005. J. Plant Nut. Soil Sci. 168, 725–740 (2005).
Launchbaugh, K. (ed.) Targeted Grazing: A Natural Approach to Vegetation Management and Landscape Enhancement (American Sheep Industry Association, 2006).
Landesamt für Landwirtschaft, Umwelt und Ländliche Räume des Landes Schleswig-Holstein (LLUR). Beweidung von Offen- und Halboffenbiotopen. Eine Adäquate Pflegemethode unter Besonderer Berücksichtigung der FFH-Lebensraumtypen und Arten (LLUR, 2010).
Schubert, R., Hilbig, R. & Klotz, S. Bestimmungsbuch der Pflanzengesellschaften Deutschlands (Spektrum Akademischer Verlag, 2009).
Pérez-Harguindeguy, N. et al. New handbook for standardized measurement of plant functional traits worldwide. Aust. J. Bot. 61, 167–234 (2013).
Kattge, J. et al. TRY plant trait database—enhanced coverage and open access. Glob. Change Biol. 26, 119–188 (2020).
Kleyer, M. et al. The LEDA Traitbase: a database of life-history traits of the Northwest European flora. J. Ecol. 96, 1266–1274 (2008).
Schroeder-Georgi, T. et al. From pots to plots: hierarchical trait-based prediction of plant performance in a mesic grassland. J. Ecol. 104, 206–218 (2016).
Wright, I. J. et al. The worldwide leaf economics spectrum. Nature 428, 821–827 (2004).
Article ADS CAS PubMed Google Scholar
Dickson, T. L., Mittelbach, G. G., Reynolds, H. L. & Gross, K. L. Height and clonality traits determine plant community responses to fertilization. Ecology 95, 2443–2452 (2014).
Eskelinen, A. & Harrison, S. Resource co-limitation governs plant community responses to altered precipitation. Proc. Natl Acad. Sci. USA 112, 13009–13014 (2015).
Article ADS CAS PubMed PubMed Central Google Scholar
Firn, J. et al. Leaf nutrients, not specific leaf area, are consistent indicators of elevated nutrient inputs. Nat. Ecol. Evol. 3, 400–406 (2019).
Crawley, M. The R Book (John Wiley & Sons, 2007).
Bates, D., Kliegl, R., Vasishth, S. & Baayen, H. Parsimonious mixed models. Preprint at arXiv https://doi.org/10.48550/arXiv.1506.04967 (2018).
R Core Team. R: A Language and Environment for Statistical Computing. https://www.R-project.org/ (R Foundation for Statistical Computing, 2020).
Pinheiro, J. & Bates, D. nlme: linear and nonlinear mixed effects models. R version 3.1-147 https://CRAN.R-project.org/package=nlme (2020).
Bates, D., Maechler, M., Bolker, B. & Walker, S. Fitting linear mixed-effects models using lme4. J. Stat. Softw. 67, 1–48 (2015).
Fox, J. & Weisberg, S. An R Companion to Applied Regression 3rd edn (Sage Publications, 2019).
We acknowledge the Helmholtz Centre for Environmental Research–UFZ, the Federal Ministry of Education and Research, the State Ministry of Science and Economy of Saxony-Anhalt and the State Ministry for Higher Education, Research and the Arts Saxony for funding the GCEF project. We thank K. Kirsch for help in maintaining the grazing treatment and taking care of the lamps at eDiValo; I. Merbach for help during the project; the staff at the Bad Lauchstädt Experimental Research Station and M. Schädler for their work in maintaining the GCEF plots and infrastructure; H. Auge, F. Buscot, S. Klotz, T. Reitz and M. Schädler for their role in setting up the GCEF; P. Portius and the UFZ workshop staff for planning and manufacturing the lamp mounts; R. Feldman for help in establishing the project; M. Krieger, E. Ramirez, J. Young, U. Rachor and A. Leyendecker for field and laboratory assistance; S. Dunker for help with choosing the lamps; H. Auge for discussions, Valoya for support in planning the lighting; G. Rada (iDiv, Media and Communications) for support with the graphics in Fig. 1 and J. S. Lefcheck for statistical advice. A.E. was supported by the Academy of Finland project no. 297191. M.-T.J .was supported by the iDiv Flexpool program grant no. 34600565-11.
Open access funding provided by University of Oulu including Oulu University Hospital.
Present address: Department of Ecology and Genetics, University of Oulu, Oulu, Finland
Department of Physiological Diversity, Helmholtz Centre for Environmental Research–UFZ, Leipzig, Germany
Anu Eskelinen, W. Stanley Harpole & Maria-Theresa Jessen
German Centre for Integrative Biodiversity Research (iDiv) Halle-Jena-Leipzig, Leipzig, Germany
Anu Eskelinen, W. Stanley Harpole & Maria-Theresa Jessen
Department of Ecology and Genetics, University of Oulu, Oulu, Finland
Anu Eskelinen & Risto Virtanen
Institute of Biology, Martin Luther University Halle-Wittenberg, Halle, Germany
Department of Community Ecology, Helmholtz Centre for Environmental Research–UFZ, Halle, Germany
Ecology and Biodiversity Group, Department of Biology, Utrecht University, Utrecht, The Netherlands
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
A.E., W.S.H. and Y.H. were responsible for the conception and design of the project. Fieldwork was performed by A.E. with help from M.-T.J. and R.V. A.E. and Y.H. were responsible for data analysis and visualization. A.E. wrote the first draft and all authors contributed to reviewing and editing.
The authors declare no competing interests.
Nature thanks Eric Allan and Niv DeMalach for their contribution to the peer review of this work.
Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Effects of fertilization and light addition on total vegetation cover % (a) and litter cover % (b) in 2017. In 2017, fencing (herbivore exclusion) treatment was not yet applied. Values are means ± s.e.m. derived from LME models where parameter significance was assessed by F-tests (2-tailed), n = 20 for quadrats from which species richness and Shannon diversity were recorded; for statistics see Methods and Extended Data Tables 2 and 5. In 2017, herbivore exclusion treatment was not yet applied. These results are from ‘eDiValo’ field experiment; see Fig. 1 and Methods.
(Species richness (a), Shannon diversity (b), and total vegetation cover (%) (c). Values are means ± s.e.m. derived from LME models where parameter significance was assessed by F-tests (2-tailed); for statistics see Methods and Extended Data Tables 2 and 5. unfenced, unlighted and lighted n = 10; fenced, unlighted and lighted n = 8; fenced and fertilized, lighted and unlighted n = 12; fertilized, unlighted and lighted n = 10. Individual data points are shown with faded colour. Sampling unit is a quadrat from which litter was estimated. The different sample sizes result from two previously unfertilized plots being accidentally fertilized in 2019 after which these plots were considered fertilized plots. These results are from ‘eDiValo’ field experiment; see Fig. 1 and Methods.
The effects of fertilization, light, and climate on total cover (%) in 2017 (a). n = 10 for quadrats from which total cover was estimated. The effects of fencing, fertilization, and climate treatment on litter cover (%) in 2019 (b). unfenced, ambient and future climate n = 10; fenced, ambient and future climate n = 8; fenced and fertilized, ambient and future climate n = 12; fertilized, ambient and future climate n = 10. Sampling unit is a quadrat from which litter was estimated. The slightly unequal samples sizes in 2019 result from two previously unfertilized plots being accidentally fertilized in 2019 after which these plots were considered fertilized plots. In both 2017 and 2019, values are means ± s.e.m. derived from LME models where parameter significance was assessed by F-tests (2-tailed); for statistics see Methods and Extended Data Tables 2 and 5. These results are from ‘eDiValo’ field experiment; see Fig. 1 and Methods.
Probability of individual species to increase in response to light addition as a function of SLA (a), and LWC (b) in 2017. Fitted lines are from GLME models with a binomial error structure where parameter significance was assessed by χ2-tests (2-tailed, see Methods): (a) SLA χ2 = 76.21, P < 0.0001, (b) LWC χ2 = 72.02, P < 0.0001. Shaded areas represent 95 % confidence intervals, n = 2440 species observations for SLA and LWC. SLA, specific leaf area; LWC, leaf water content. These results are from ‘eDiValo’ field experiment; see Fig. 1 and Methods.
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/.
Eskelinen, A., Harpole, W.S., Jessen, MT. et al. Light competition drives herbivore and nutrient effects on plant diversity. Nature (2022). https://doi.org/10.1038/s41586-022-05383-9
DOI: https://doi.org/10.1038/s41586-022-05383-9
Anyone you share the following link with will be able to read this content:
Sorry, a shareable link is not currently available for this article.
Provided by the Springer Nature SharedIt content-sharing initiative
By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.
Nature (Nature) ISSN 1476-4687 (online) ISSN 0028-0836 (print)
Sign up for the Nature Briefing newsletter — what matters in science, free to your inbox daily.