Diversity Promotes Temporal Stability across Levels of Ecosystem Organization in Experimental Grasslands

Abstract

The diversity–stability hypothesis states that current losses of biodiversity can impair the ability of an ecosystem to dampen the effect of environmental perturbations on its functioning. Using data from a long-term and comprehensive biodiversity experiment, we quantified the temporal stability of 42 variables characterizing twelve ecological functions in managed grassland plots varying in plant species richness. We demonstrate that diversity increases stability i) across trophic levels (producer, consumer), ii) at both the system (community, ecosystem) and the component levels (population, functional group, phylogenetic clade), and iii) primarily for aboveground rather than belowground processes. Temporal synchronization across studied variables was mostly unaffected with increasing species richness. This study provides the strongest empirical support so far that diversity promotes stability across different ecological functions and levels of ecosystem organization in grasslands.

About the Authors:

Raphaël Proulx

* E-mail: [email protected]

Affiliation: Max Planck Institute for Biogeochemistry, Hans-Knöll Strasse, Jena, Germany

Christian Wirth

Affiliations Max Planck Institute for Biogeochemistry, Hans-Knöll Strasse, Jena, Germany, Department of Special Botany and Functional Biodiversity Research, University of Leipzig, Leipzig, Germany

Winfried Voigt

Affiliation: Hydrology, Institute of Ecology, University of Jena, Jena, Germany

Alexandra Weigelt

Affiliations Department of Special Botany and Functional Biodiversity Research, University of Leipzig, Leipzig, Germany, Hydrology, Institute of Ecology, University of Jena, Jena, Germany

Christiane Roscher

Affiliation: Max Planck Institute for Biogeochemistry, Hans-Knöll Strasse, Jena, Germany

Sabine Attinger

Affiliation: UFZ Helmholtz Center for Environmental Research, Leipzig, Germany

Jussi Baade

Affiliation: Institute of Geography, Friedrich Schiller University Jena, Jena, Germany

Romain L. Barnard

Affiliation: Institute for Plant, Animal and Agroecosystem Sciences, ETH Zurich, Zürich, Switzerland

Nina Buchmann

Affiliation: Institute for Plant, Animal and Agroecosystem Sciences, ETH Zurich, Zürich, Switzerland

François Buscot

Affiliation: Soil Ecology, UFZ Helmholtz Centre for Environmental Research, Halle, Germany

Nico Eisenhauer

Affiliations Institute of Zoology, Darmstadt University of Technology, Darmstadt, Germany, J.F. Blumenbach Institute of Zoology and Anthropology, Georg-August-University Göttingen, Göttingen, Germany

Markus Fischer

Affiliations Institute of Plant Sciences, University of Berne, Berne, Switzerland, Institute of Biochemistry and Biology, University of Potsdam, Potsdam, Germany

Gerd Gleixner

Affiliation: Max Planck Institute for Biogeochemistry, Hans-Knöll Strasse, Jena, Germany

Stefan Halle

Affiliation: Hydrology, Institute of Ecology, University of Jena, Jena, Germany

Anke Hildebrandt

Affiliation: UFZ Helmholtz Center for Environmental Research, Leipzig, Germany

Esther Kowalski

Affiliation: Hydrology, Institute of Ecology, University of Jena, Jena, Germany

Annely Kuu

Affiliations Hydrology, Institute of Ecology, University of Jena, Jena, Germany, Tartu College, Tallinn University of Technology, Tartu, Estonia

Markus Lange

Affiliation: Hydrology, Institute of Ecology, University of Jena, Jena, Germany

Alex Milcu

Affiliation: Center for Population Biology, Imperial College London, Silwood Park, United Kingdom

Pascal A. Niklaus

Affiliations Institute for Plant, Animal and Agroecosystem Sciences, ETH Zurich, Zürich, Switzerland, Institute of Evolutionary Biology and Environmental Studies, University of Zürich, Zürich, Switzerland

Yvonne Oelmann

Affiliation: Geographic Institute, Johannes Gutenberg University, Mainz, Germany

Stephan Rosenkranz

Affiliation: Geographic Institute, Johannes Gutenberg University, Mainz, Germany

Alexander Sabais

Affiliation: Institute of Zoology, Darmstadt University of Technology, Darmstadt, Germany

Christoph Scherber

Affiliation: Department of Crop Sciences, University of Göttingen, Göttingen, Germany

Michael Scherer-Lorenzen

Affiliations Institute for Plant, Animal and Agroecosystem Sciences, ETH Zurich, Zürich, Switzerland, Faculty for Biology, University of Freiburg, Freiburg, Germany

Stefan Scheu

Affiliation: J.F. Blumenbach Institute of Zoology and Anthropology, Georg-August-University Göttingen, Göttingen, Germany

Ernst-Detlef Schulze

Affiliation: Max Planck Institute for Biogeochemistry, Hans-Knöll Strasse, Jena, Germany

Jens Schumacher

Affiliation: Institute for Stochastics, University of Jena, Jena, Germany

Guido Schwichtenberg

Affiliation: UFZ Helmholtz Center for Environmental Research, Leipzig, Germany

Jean-François Soussana

Affiliation: Grassland Ecosystem Research, Institut National de la Recherche Agronomique, Clermont-Ferrand, France

Vicky M. Temperton

Affiliations Max Planck Institute for Biogeochemistry, Hans-Knöll Strasse, Jena, Germany, Phytosphere Institute ICG-3, Forschungszentrum Jülich GmbH, Jülich, Germany

Wolfgang W. Weisser

Affiliation: Hydrology, Institute of Ecology, University of Jena, Jena, Germany

Wolfgang Wilcke

Affiliations Geographic Institute, Johannes Gutenberg University, Mainz, Germany, Geographic Institute, University of Berne, Berne, Switzerland

Bernhard Schmid

Affiliation: Institute of Evolutionary Biology and Environmental Studies, University of Zürich, Zürich, Switzerland

Introduction

Ecosystems are subjected to natural environmental perturbations ranging from small- to large-scale processes. Species-rich communities host a variety of life strategies that can respond differently to environmental perturbations and contribute to ecological functioning in various ways, thus increasing ecosystem stability [1]–[4]. However, the mechanisms by which species-specific variations in response to perturbations translate into ecosystem stability are still debated [5], [6] and it remains largely unknown whether the diversity–stability hypothesis holds for ecological functions other than plant community biomass. In this context, temporal stability is defined as the capacity of an ecosystem to dampen environmental perturbations while retaining the ecological function of interest [1]–[4].

In a recent meta-analysis, Jiang and Pu [7] showed that biodiversity stabilizes community variables and may additionally stabilize species populations in multi-trophic systems. Their result led to the proposition that in the presence of temporal synchronization —when different variables respond similarly to environmental perturbations— stabilization of ecological functions at the component level (e.g., population) of organization may promote stabilization at the system level (e.g., community). This proposition contrasts with previous theories that predicted populations to be destabilized and communities to be stabilized [4], [8], but is in agreement with other conceptual models [9]–[11]. Acknowledging a possible bottom-up effect of species richness on the temporal stability of ecological functions at the community level of organization, two corollaries follow. First, species richness decreases the temporal variance of ecological functions at both the component level of species (populations of single species or groups of functionally similar species) and the system level of the community (plant community or ecosystem property). Second, the temporal co-variance between variables characterizing an ecological function (e.g., between the biomass production of plant functional groups, such as herbs, grasses and legumes, contributing to the aboveground plant biomass production) is positive and not affected by species richness. In other words, ecological functions at lower (population) levels of organization may be stabilized by species richness and, if variables co-vary positively in time, this stabilizing effect would propagate through the whole ecosystem. However, a critical evaluation of the above two corollaries can only be accomplished with an experimental approach. Furthermore, in spite of recent meta-analyses indicating that ecological functioning is enhanced by higher biodiversity [12]–[14], it is not clear whether biodiversity concurrently increases the temporal stability of these functions across many levels of organization. We analyzed this question for multiple ecological functions measured in a single experiment.

Using 7 years of data collected on 82 large plots in a grassland biodiversity experiment (the Jena Experiment) [15], [16], we tested to what degree higher species richness of plant communities translates into temporal stability when functions are characterized at the organizational level of the population, the functional and phylogenetic group, the community, or the ecosystem (Table 1). Naturally-occurring or management-caused environmental perturbations in the Jena Experiment (50°57′3.09″N, 11°37′23.49″E) include summer droughts, winter soil freezing, mowing of the vegetation twice a year, and spatial heterogeneity in vegetation cover through weeding. In this experiment, the lower resistance of less diverse communities against spontaneously invading species was shown to occur independently of the spatial heterogeneity in vegetation cover [17]. Functions at different organizational levels refer to: the species level (populations of species), the functional and phylogenetic level (groups of similar species), the community level of primary producers and secondary consumers, and the ecosystem level (mainly biogeochemical soil properties).

[Figure omitted. See PDF.]

Table 1. Summary of 42 variables used to characterize the stability of twelve ecological functions.

https://doi.org/10.1371/journal.pone.0013382.t001

Kolasa and Li [9] demonstrated that habitat specialization usually increases with species richness and, because the abundance of habitat-specialist species is generally more variable over time, this may mask the stabilizing effect of other species. Here we equated specialization with low abundance and thus excluded sub-dominant species, as well as sub-dominant functional or phylogenetic groups of species; a procedure compatible with Kolasa and Li's formalism. Our approach has the main advantage that the same set of target variables is used to evaluate the effect that varying plant species richness has on temporal stability for different ecological functions and levels of organization (see Materials and Methods).

As a measure of temporal stability [3], [4], we calculated the multivariate coefficient of variation (CV2) by summing the variances and co-variances of scaled variables characterizing the same ecological function over time as follows (see Materials and Methods): CV2 =  [Σ Variances +Σ Co-variances] / [Σ Means]2. The reciprocal of the numerator reflects the first part of the temporal stability definition of the introduction paragraph: “the capacity of an ecosystem to dampen environmental perturbations…”. The denominator reflects the second part of that definition: … “ while retaining the ecological function of interest.” A comparatively lower sum of variances scaled by the square of the summed means (hereafter the variance CV) among variables characterizing an ecological function indicates higher stability, since it implies a dampening of response to a common set of environmental perturbations. Whereas a comparatively lower sum of co-variances scaled by the square of the summed means (hereafter the co-variance CV) among these variables also indicates higher stability since it implies that different facets of an ecological function synchronize less over time, whereby minimizing the risk of functional breakdown. From the above approach, the population is characterized by variables A, B, C measured at the species level, while the community is characterized by another set of variables D, E, F measured at the community level. This notably contrasts with previous approaches [3], [4], [9] where the community level variable D is the sum of A, B, and C.

Our methodology aims at determining how plant species richness influences the temporal variances and co-variances of target variables grouped according to their level of organization and ecological function. However, the CV2 measure of variability rests on the mathematical postulate that the variance and the squared mean of a measured variable scale positively and linearly [18]. When more than one variable is used to characterize an ecological function, both the summed variance and co-variance components of the CV2 measure are expected to scale linearly with the squared sum of means. More specifically, if the variables are either positively or negatively correlated with one another, the summed co-variances (expressed in absolute terms) reflect the summed variances and both components scale with the squared sum of means. On the other hand, if variables are de-correlated (independent) of one another, the summed co-variances tend towards zero and only the summed variances scale with the squared sum of means. Thus, to be a valid comparative measure of temporal stability, the CV ratio must account for the fact that summed variances and co-variances are not mathematically independent of the summed means. In our study, the use of variance and co-variance CV ratios was justified on two points. i) For all twelve ecological functions, at the alpha rejection rate of 0.05, we found a significantly positive Pearson's r correlation across experimental plots between the squared sum of means and both the summed variances and the summed co-variances expressed in absolute terms. The mean ± 1SD correlation across the twelve ecological functions was 0.743±0.136 (min 0.484, max 0.923) for the variances and 0.525±0.235 (min 0.230, max 0.849) for the co-variances. ii) The slope of the log-log relationships between summed variances (or summed co-variances) and the squared sum of means approached one, indicating that a vast majority of correlations were linear. The mean ± 1SD linear slope across the twelve ecological functions was 0.767±0.183 (min 0.398, max 1.081) for the variances and 0.873±0.236 (min 0.4284, max 1.128) for the co-variances.

Results

Our synthesis of the diversity–stability hypothesis in the Jena Experiment emphasizes two main results. Firstly, ecological functions at various organizational levels were stabilized with increasing plant species richness, as indicated by a decrease of the variance CV for the abundance of parasitic hymenoptera (food-web complexity), the suppression of non-resident plant species (bioregulation), vegetation structure and biomass production (primary producers), the abundance and diversity of invertebrates (secondary consumers), as well as trace gas fluxes (ecosystem properties). We found that the variance CV of ecological functions processes (Production of Earthworm Biomass and Below-Ground Invertebrates, Soil Nutrient and Water Content) associated to belowground processes was not significantly decreased with increasing species richness.

Secondly, the co-variance CV of a majority (9 out of 12) of ecological functions did not show a relationship with species richness. However, the variables characterizing ecological functions at the community level were more synchronized (more dependent on one another) at low than high species richness, as indicated by the increased co-variance CV in monocultures (Figure 1). Furthermore, the proportion of positive co-variance CV across all ecological functions and experimental plots attained 85% (Figure 2), suggesting that temporal synchronization among variables is common.

[Figure omitted. See PDF.]

Figure 1. Diversity–stability relationships for the twelve ecological functions grouped by levels of organization.

Each dot represents the mean of temporal stability values (vertical axes) obtained across experimental plots sowed with the same number of plant species (1, 2, 4, 8, 16 or 60 species). Error bars represent the 95% confidence interval around the mean. Each temporal stability measure was standardized with a mean of zero and variance of one to ease the comparison of diversity–stability relationships within and between organizational levels. The left panels (A–D) show the relationships between plant species richness and the variance CV, while the right panels (E–H) show the relationships to the co-variance CV. No measurements were available for species richness treatments 2 and 8 in Earthworm Biomass, Parasitic Hymenoptera, Below- and Aboveground Invertebrates, and Arthropod Diversity production functions. For each component of temporal stability we report Pearson's correlation coefficient r estimated by a linear fit between the logarithm of plant species richness and one component of temporal stability. The number of stars next to Pearson's r values gives the probability of accepting the null hypothesis following a distribution-free randomization test [37] (df are in Table 1): Blank (p>0.05), *(p<0.05), **(p<0.01), ***(p<0.001).

https://doi.org/10.1371/journal.pone.0013382.g001

[Figure omitted. See PDF.]

Figure 2. Proportion of positive co-variance CV across experimental plots in each ecological function.

EB: Earthworm Biomass; PH: Parasitic Hymenoptera; IPBP: Invasive Plant Bioregulation of Populations; BGI: Below-Ground Invertebrates; AGI: Above-Ground Invertebrates; PFGB: Plant Functional Group Biomass; PSS: Plant Stand Structure; IPBC: Invasive Plant Bioregulation of Communities; AD: Arthropod Diversity; SWC: Soil Water Content; SNC: Soil Nutrient Concentration; TGF: Trace Gas Fluxes. See Table 1 for additional details on the number of experimental plots in each ecological function.

https://doi.org/10.1371/journal.pone.0013382.g002

Discussion

Our finding that the variance CV of ecological functions decreased with increasing species richness corroborates previous finding about the stabilization of aboveground productivity by biodiversity [3], [19]–[21] and extends it to further functions. We furthermore show that the diversity–stability relationship holds across multiple ecological functions in a single experiment, which is remarkable considering the spectrum of organizational (from population to ecosystem) and trophic levels (from primary producers to higher consumers) investigated. Our results lend support to the proposition that temporal stabilization at the lower levels of organization can promote stabilization at the higher levels (communities, ecosystems) [7], [21]. Some examples from our experiment illustrate this point in more detail. The temporal stability of parasitic wasps could depend on the continuous presence of a diversity of hosts, which in turn would depend on the vegetation structure and the biomass of plant functional groups. Temporal stability against plant invasion by non-resident species may be linked to resource pre-emption; i.e., the environmental resources are exploited by the resident plant species. At the sub-community level, the temporal stability of aboveground invertebrate abundances could be associated with specific plant–insect interactions. Temporal stability in arthropod biodiversity may represent a circumstance where the diversity in grassland plant species promotes diversity at other trophic levels [22].

Unfortunately, our data did not allow disentangling among these alternative hypotheses because time series between ecological functions were either too short or not temporally aligned. Nevertheless, the general picture suggests that species groups that are dominant in the plant community can exert a bottom-up control on the stability of aboveground ecological functions at higher organizational levels. A recent meta-analysis [12] reported that out of sixteen candidate hypotheses, the organizational level had the strongest influence in explaining the magnitude of the biodiversity effect on ecological functions. Our finding that multiple ecological functions are concurrently stabilized by plant species richness is a strong indication that identifying diversity–stability mechanisms requires a complete overview of the identities of at least all dominant species and ecological interactions within and between organizational levels.

Previous studies on grassland systems suggested that diversity promotes temporal stability at the community level, but that the stability of individual species within communities may show the opposite trend [3], [19], [20], [23]. Our study is particular in the sense that, at organizational levels below the community, we considered only variables measured on the dominant populations of species and the dominant groups of functionally similar species. Our results suggest that the increasing response of the variance CV for plant species biomass reported in earlier studies may have been due to the accumulation of habitat specialists at low abundance with increasing species richness [9]. In plant communities, the accumulation of specialists is paralleled by a reduction in the individual species biomass with increasing richness, which follows from the partitioning of total biomass among species per unit area. The variance CV of a few dominant plant species may in fact decrease with increasing species richness, but this stabilizing effect would be confounded by the opposite response of a larger number of sub-dominant species. This interpretation is also in agreement with the “weak interaction effect”, which states that the dynamics of dominant, strongly interacting, species is stabilized by the presence of less dominant species [2], [24]. Experiments on both natural [25] and theoretical systems [24] suggest that distributions of interaction strengths in species-rich communities are typically skewed towards many weak and few strong species–species interactions.

We could detect only a subtle stabilizing effect of plant species richness on the variance CV of the trace-gas-flux function at the ecosystem level, while the observed stabilizing effect was consistently stronger on ecological functions at the community level (Figure 1). One possibility is that weaker diversity–stability relationships at the ecosystem level are due to the presence of soil compensatory mechanisms, such as inherited carbon pools from previous land management and time-lag responses of the belowground communities. This explanation is supported by at least two recent studies from the Jena Experiment [26], [27], revealing that microbial activity and carbon sequestration in soils were significantly increased only 4 years after the beginning of the experiment. However, we also found that the variance CV of other belowground ecological functions (earthworm biomass production, belowground invertebrate abundances) did not decrease significantly with increasing plant species richness. This suggests a possible link between the nature of the experimental manipulation and temporal stability. In managed grasslands, the belowground rooting system and soil texture remain comparatively unaltered by periodically harvesting the aboveground plant biomass, therefore creating different perturbation regimes.

Following their review of the literature on the diversity–stability hypothesis, Jiang and Pu [7] proposed that temporal stabilization of variables at lower organizational levels could promote stabilization at the higher levels (community, ecosystem); i.e. they proposed what we here called bottom-up effects of species richness on the temporal stability of ecological functions at the community and ecosystem levels of organization. Their proposition is valid even when temporal synchronization remains unaffected with increasing species richness; an assertion largely supported by our results on the co-variance CV (Figures 1 and 2). When different variables characterizing an ecological function are strongly correlated with one another (either positively or negatively), all co-variances reflect the variances. This statistical effect alone could explain the decrease of the co-variance CV of ecological functions with increasing species richness that we observed at the community level (Figure 1), where variables are strongly correlated. Conversely, we observed no effect of increasing plant species richness on the co-variance CV of ecological functions below the community level of organization. This indicates that although a majority of variables characterizing an ecological function at the species and phylogenetic (functional) group levels were in general positively correlated (Figure 2), these variables are not as strongly correlated as those measured at the community and ecosystem levels.

The strong positive correlations among variables measured at the community and ecosystem levels may also result from allometric relationships between plant height, cover and biomass (Plant Stand Structure and Invasive Plant Bioregulation functions), from mathematical relations linking species abundance and richness in ecological communities (Arthropod Diversity function), and from spatial auto-correlation among observations (Soil Water and Nutrient Content functions). Although other factors may cause variables to be more strongly correlated at the higher levels of organization, this should only strengthen the possibility that increasing plant species richness can stabilize ecosystem functions in a bottom-up fashion.

What would cause populations of dominant group of species in a community to co-vary positively (synchronize) in time? Through an extensive review of 41 natural and experimental communities, Houlahan and colleagues [28] concluded that species–environment interactions largely dominate community dynamics, hence driving species populations to co-vary positively in their abundance and biomass. However, it may be premature to conclude whether a positive co-variance indeed reflects the synchronized response of species populations sharing a common set of environmental perturbations [6], [29]. More work will be needed to disentangle whether co-variances are the result of: i) species–species interactions (e.g., synchronous response of several sub-dominant species to a few dominant ones), ii) species–environment interactions (e.g., synchronous response of species to a common set of environmental perturbations), or iii) both interaction types (e.g., synchronous response of consumers to temporal variations in the producers).

Our study advocates a bottom-up effect of plant species populations on the temporal stability of community- and ecosystem-level functions in managed grasslands. Novel conceptual frameworks for the diversity–stability hypothesis should now attempt to go beyond the grouping of ecological functions into species and community levels of organization to fully account for the notion of spatial and temporal observation scale. At small observation scales the population dynamics of a few individuals is captured, while ecosystem dynamics need larger observation scales to consider feedbacks on the population dynamics. Additionally, a more practical quantification of biological interactions is needed that embraces the vast array of species–species competitive and multi-trophic interactions. This will require a better theoretical and empirical foundation of how the pattern of interaction strength between species affects the functioning of whole ecosystems. While managing ecosystem stability in a stochastic and changing world is beyond current abilities [30], preserving or creating diverse ecosystems so far remains our best option to achieve the stabilization of ecosystems.

Materials and Methods

Experimental Design

In the Jena Experiment, regularly mown and weeded grassland plots of 20×20 m were established with plant species richness of 1, 2, 4, 8, 16, and 60 species. At each species richness level, 16 experimental plots were created randomly from a pool of 60 species, with the exception of richness levels 16 and 60 for which 14 and 4 plots were available, respectively. Of the measurements conducted in the course of the study, 42 variables were measured on six or more dates and fulfilled the following criteria: the variable contained less than 50% of zeros across all experimental plots and temporal observations, the variable has a temporal variance sufficiently far from zero to prevent the inclusion of processes that do not respond to environmental perturbations or have ceased to operate, and the variable could be temporally aligned along with other variables characterizing the same ecological function. Additionally, for the quantification of belowground invertebrates from Kempson cores (Table 1) two temporal observations were made within a few days in experimental subplots subjected to different treatments: with and without the application of an equal amount of chlorpyrifos.

The 42 variables used in this study were grouped into twelve ecological functions according to their similarity and level of organization (a short description of each ecological function is given in parentheses): 1) Earthworm Biomass production (soil engineering; population level), 2) Parasitic Hymenoptera Abundance (food-web complexity; population level), 3) Invasive Plant Bioregulation (suppression of spontaneously invading plant species; population level), 4) Belowground Invertebrate Abundance (density of dominant soil macrofauna functional groups; phylogenetic clade level), 5) Aboveground Invertebrate Abundance (density of dominant aboveground, mainly saprophagous, phytophagous and parasitic, insect functional groups; phylogenetic clade level), 6) Plant Functional Group Biomass production (aboveground plant productivity; phylogenetic clade level), 7) Plant Stand Structure (vegetation development in cover, height, and biomass; community level), 8) Invasive Plant Bioregulation (suppression of spontaneously invading vegetation in cover, species diversity, and biomass; community level), 9) Arthropod Diversity (arthropod abundance and richness; community level), 10) Soil Water Content (use of water by organisms; ecosystem level), 11) Soil Nutrient Concentration (use of ammonium and nitrate by organisms; ecosystem level), and 12) Trace Gas Fluxes (net ecosystem–atmosphere exchange of CO2, N2O and CH4; ecosystem; ecosystem level).

The main references for the field sampling protocols listed in Table 1 are: Trace Gas Fluxes [31]; Soil Nutrient Concentration [32]; Soil Water Content [33]; Arthropod Diversity, Above-Ground Invertebrates, and Parasitic Hymenoptera [34]; Plant Stand Structure and Plant Functional Group Biomass [16]; Invasive Plant Bioregulation [35]; Below-Ground Invertebrate and Earthworm Biomass [36].

Our primary aim was to include only species, or more generally variables, observed at least once in all experimental plots. This approach has two main advantages. It focuses on a common set of variables that capture the functioning of all grassland plots across the species richness gradient. It avoids confounding the effect of habitat specialists (i.e., sub-dominant species in the present context) with measures of temporal stability at the population level [9]. Furthermore, because plant species were randomly attributed to the different communities in our experiment, each species is only found in a limited number of experimental plots. Therefore, we did not include population variables characterizing the productivity of sown (resident) plant species.

Temporal stability measures

We linearly scaled variables between zero and one to remove the effect of measurement units. The scaling of a variable x was performed across all experimental plots and temporal observations as follows: x-min(x) /range(x). Variables characterizing the Invasive Plant Bioregulation functions (population and community levels), as well as the Soil Water Content and Nutrient Concentration functions (ecosystem level), were expressed in terms of ecological services and were therefore transformed by taking the reciprocal as follows: 1/ (x+0.01). Adding 0.01 to these scaled x variable prevented division by zero and preserved a log-linear relationship between the squared mean and variance. Then, in each experimental plot, and for each of the twelve ecological functions, we constructed a matrix containing temporal observations in rows and variables in columns. We obtained from each matrix the multivariate coefficient of variation as follows: CV2 =  [Σ Variances +Σ Co-variances]/ [Σ Means]2. We partitioned the squared coefficient of variation CV2 into two additive components: an insurance effect (Variance CV =  [Σ Variances] / [Σ Means]2) and a synchronizing effect (Co-variance CV =  [Σ Co-variances] / [Σ Means]2). Although summed variances and co-variances have been recently criticized as a means of precisely identifying ecological mechanisms [6], [29], their ability to capture generic aspects of temporal stability in ecosystems is not in question.

The above measures of temporal stability can be defined in the context of previous ecological theories of the diversity–stability hypothesis. An insurance (variance dampening) effect results from functional complementarity among co-occurring species. By increasing the number of plant species, the insurance effect predicts a decrease in the variance CV of ecological functions, which may occur through a differentiation in resource use [1], [4] and a complex network of many weak interactions among species [2], [24]. With increasing species richness, the insurance effect is associated with an increase in the mean level of functionality, which is not exceeded by a similar increase in the sum of variances. On the other hand, the synchronizing (positive co-variance) effect is associated to a temporal coupling among co-occurring species. By increasing the number of plant species, the co-variance CV of ecological functions is predicted to decrease because compensatory interactions intensify among species competing for limited resources [4], [28] or because species are differently sensitive to environmental perturbations [5], [10]. With increasing species richness, the synchronizing effect is associated with a decrease in the co-variances among variables jointly characterizing the same ecological function.

Acknowledgments

We would like to thank the gardeners and numerous student helpers for weeding and mowing of the experimental plots and help with data collection. We also thank Michel Loreau for providing helpful comments on a previous version of this manuscript. We acknowledge Teja Tscharntke, Elisabeth Marquard, Holger Bessler, and Christof Engels for their assistance in designing the Jena Experiment.

Author Contributions

Conceived and designed the experiments: CW WV NB FB MF GG SH ML MSL SS EDS JS JFS WW WW BS. Performed the experiments: RP AW CR RLB NE EK AK AM PAN YO AS CS VMT. Analyzed the data: RP PAN CS JS. Contributed reagents/materials/analysis tools: CW WV AW CR SA JB RLB NE GG AH EK AK AM SR GS. Wrote the paper: RP CW WV AW CR SA JB RLB NB FB NE MF GG SH AH EK AK ML AM PAN YO SR AS CS MSL SS EDS JS GS JFS VMT WW WW BS.

Citation: Proulx R, Wirth C, Voigt W, Weigelt A, Roscher C, Attinger S, et al. (2010) Diversity Promotes Temporal Stability across Levels of Ecosystem Organization in Experimental Grasslands. PLoS ONE5(10): e13382. https://doi.org/10.1371/journal.pone.0013382

References

1. Yachi S, Loreau M (1999) Biodiversity and ecosystem productivity in a fluctuating environment: The insurance hypothesis. Proc Natl Acad Sci U S A 96: 1463–1468.S. YachiM. Loreau1999Biodiversity and ecosystem productivity in a fluctuating environment: The insurance hypothesis.Proc Natl Acad Sci U S A9614631468

2. McCann KS (2000) The diversity-stability debate. Nature 405: 228–233.KS McCann2000The diversity-stability debate.Nature405228233

3. Tilman D, Reich PB, Knops JHM (2006) Biodiversity and ecosystem stability in a decade-long grassland experiment. Nature 441: 629–632.D. TilmanPB ReichJHM Knops2006Biodiversity and ecosystem stability in a decade-long grassland experiment.Nature441629632

4. Lehman CL, Tilman D (2000) Biodiversity, stability, and productivity in competitive communities. Am Nat 156: 534–552.CL LehmanD. Tilman2000Biodiversity, stability, and productivity in competitive communities.Am Nat156534552

5. Ives AR, Carpenter SR (2007) Stability and diversity of ecosystems. Science 317: 58–62.AR IvesSR Carpenter2007Stability and diversity of ecosystems.Science3175862

6. Loreau M, de Mazancourt C (2008) Species synchrony and its drivers: neutral and nonneutral community dynamics in fluctuating environments. Am Nat 172: E48–E66.M. LoreauC. de Mazancourt2008Species synchrony and its drivers: neutral and nonneutral community dynamics in fluctuating environments.Am Nat172E48E66

7. Jiang L, Pu Z (2009) Different effects of species diversity on temporal stability in single-trophic and multitrophic communities. Am Nat 174: 651–659.L. JiangZ. Pu2009Different effects of species diversity on temporal stability in single-trophic and multitrophic communities.Am Nat174651659

8. May RM (1972) What is the chance that a large complex system will be stable? Nature 237: 413–414.RM May1972What is the chance that a large complex system will be stable?Nature237413414

9. Kolasa J, Li B-L (2003) Removing the confounding effect of habitat specialization reveals the stabilizing contribution of diversity to species variability. Proc R Soc Lond B 270: S198–S201.J. KolasaB-L Li2003Removing the confounding effect of habitat specialization reveals the stabilizing contribution of diversity to species variability.Proc R Soc Lond B270S198S201

10. Ives AR, Gross K, Klug JL (1999) Stability and variability in competitive communities. Science 286: 542–544.AR IvesK. GrossJL Klug1999Stability and variability in competitive communities.Science286542544

11. Romanuk TN, Kolasa J (2004) Population variability is lower in diverse rock pools when the obscuring effects of local processes are removed. Ecoscience 11: 455–462.TN RomanukJ. Kolasa2004Population variability is lower in diverse rock pools when the obscuring effects of local processes are removed.Ecoscience11455462

12. Balvanera P, Pfisterer AB, Buchmann N, He JS, Nakashizuka T, et al. (2006) Quantifying the evidence for biodiversity effects on ecosystem functioning and services. Ecol Lett 9: 1146–1156.P. BalvaneraAB PfistererN. BuchmannJS HeT. Nakashizuka2006Quantifying the evidence for biodiversity effects on ecosystem functioning and services.Ecol Lett911461156

13. Cardinale BJ, Srivastava DS, Duffy JE, Wright JP, Downing AL, et al. (2006) Effects of biodiversity on the functioning of trophic groups and ecosystems. Nature 443: 989–992.BJ CardinaleDS SrivastavaJE DuffyJP WrightAL Downing2006Effects of biodiversity on the functioning of trophic groups and ecosystems.Nature443989992

14. Schmid B, Balvanera P, Cardinale BJ, Godbold J, Pfisterer AB, et al. (2009) Consequences of species loss for ecosystem functioning: Meta-analyses of data from biodiversity experiments. In: Naeem S, Bunker DE, Hector A, Loreau M, Perrings C, editors. Biodiversity, Ecosystem Functioning, and Human Wellbeing: An Ecological and Economic Perspective. Oxford: Oxford University Press. pp. 14–29.B. SchmidP. BalvaneraBJ CardinaleJ. GodboldAB Pfisterer2009Consequences of species loss for ecosystem functioning: Meta-analyses of data from biodiversity experiments.S. NaeemDE BunkerA. HectorM. LoreauC. PerringsBiodiversity, Ecosystem Functioning, and Human Wellbeing: An Ecological and Economic PerspectiveOxfordOxford University Press1429

15. Roscher C, Schumacher J, Baade J, Wilcke W, Gleixner G, et al. (2004) The role of biodiversity for element cycling and trophic interactions: an experimental approach in a grassland community. Basic Appl Ecol 5: 107–121.C. RoscherJ. SchumacherJ. BaadeW. WilckeG. Gleixner2004The role of biodiversity for element cycling and trophic interactions: an experimental approach in a grassland community.Basic Appl Ecol5107121

16. Weigelt A, Marquard E, Temperton VM, Roscher C, Scherber C, et al. (2010) The Jena Experiment: six years of data from a grassland biodiversity experiment. Ecology 91: 930–931.A. WeigeltE. MarquardVM TempertonC. RoscherC. Scherber2010The Jena Experiment: six years of data from a grassland biodiversity experiment.Ecology91930931

17. Roscher C, Temperton VM, Buchmann N, Schulze E-D (2009) Community assembly and biomass production in regularly and never weeded experimental grasslands. Acta Oecologia 35: 206–217.C. RoscherVM TempertonN. BuchmannE-D Schulze2009Community assembly and biomass production in regularly and never weeded experimental grasslands.Acta Oecologia35206217

18. Tokeshi M (1995) On the Mathematical Basis of the Variance-Mean Power Relationship. Res Popul Ecol 37: 43–48.M. Tokeshi1995On the Mathematical Basis of the Variance-Mean Power Relationship.Res Popul Ecol374348

19. Caldeira MC, Hector A, Loreau M, Pereira JS (2005) Species richness, temporal variability and resistance of biomass production in Mediterranean grassland. Oikos 110: 115–123.MC CaldeiraA. HectorM. LoreauJS Pereira2005Species richness, temporal variability and resistance of biomass production in Mediterranean grassland.Oikos110115123

20. Flynn DFB, Schmid B, He J-S, Wolfe-Bellin KS, Bazzaz FA (2008) Hierarchical reliability in experimental plant assemblages. J Plant Ecol 1: 59–65.DFB FlynnB. SchmidJ-S HeKS Wolfe-BellinFA Bazzaz2008Hierarchical reliability in experimental plant assemblages.J Plant Ecol15965

21. Bai Y, Han X, Wu J, Chen Z, Li L (2004) Ecosystem stability and compensatory effects in the Inner Mongolia grassland. Nature 431: 181–184.Y. BaiX. HanJ. WuZ. ChenL. Li2004Ecosystem stability and compensatory effects in the Inner Mongolia grassland.Nature431181184

22. Siemann E, Tilman D, Haarstad J, Ritchie M (1998) Experimental tests of the dependence of arthropod diversity on plant diversity. Am Nat 152: 738–750.E. SiemannD. TilmanJ. HaarstadM. Ritchie1998Experimental tests of the dependence of arthropod diversity on plant diversity.Am Nat152738750

23. Hector A, Hautier Y, Saner P, Wacker L, Bagchi R, et al. (2010) General stabilizing effects of plant diversity on grassland productivity through population asynchrony and overyielding. Ecology 1: 2213–2220.A. HectorY. HautierP. SanerL. WackerR. Bagchi2010General stabilizing effects of plant diversity on grassland productivity through population asynchrony and overyielding.Ecology122132220

24. Christianou M, Kokkoris GD (2008) Complexity does not affect stability in feasible model communities. J Theor Biol 253: 162–169.M. ChristianouGD Kokkoris2008Complexity does not affect stability in feasible model communities.J Theor Biol253162169

25. Wootton JT, Emmerson M (2005) Measurement of interaction strength in nature. Ann Rev Ecol Evol Syst 36: 419–444.JT WoottonM. Emmerson2005Measurement of interaction strength in nature.Ann Rev Ecol Evol Syst36419444

26. Steinbeiss S, Beβler H, Engels C, Temperton VM, Buchmann N, et al. (2008) Plant diversity positively affects short-term soil carbon storage in experimental grasslands. Glob Change Biol 14: 2937–2949.S. SteinbeissH. BeβlerC. EngelsVM TempertonN. Buchmann2008Plant diversity positively affects short-term soil carbon storage in experimental grasslands.Glob Change Biol1429372949

27. Eisenhauer N, Bessler H, Engels C, Gleixner G, Habekost M, et al. (2010) Plant diversity effects on soil microorganisms support the singular hypothesis. Ecology 91: 485–496.N. EisenhauerH. BesslerC. EngelsG. GleixnerM. Habekost2010Plant diversity effects on soil microorganisms support the singular hypothesis.Ecology91485496

28. Houlahan JE, Currie DJ, Cottenie K, Cumming GS, Ernest SK, et al. (2007) Compensatory dynamics are rare in natural communities. Proc Natl Acad Sci U S A 104: 3273–3277.JE HoulahanDJ CurrieK. CottenieGS CummingSK Ernest2007Compensatory dynamics are rare in natural communities.Proc Natl Acad Sci U S A10432733277

29. Ranta E, Kaitala V, Fowler MS, Laakso J, Ruokolainen L, et al. (2008) Detecting compensatory dynamics in competitive communities under environmental forcing. Oikos 117: 1907–1911.E. RantaV. KaitalaMS FowlerJ. LaaksoL. Ruokolainen2008Detecting compensatory dynamics in competitive communities under environmental forcing.Oikos11719071911

30. Srivastava DS, Vellend M (2005) Biodiversity-ecosystem function research: Is it relevant to conservation. Ann Rev Ecol Evol Syst 36: 267–294.DS SrivastavaM. Vellend2005Biodiversity-ecosystem function research: Is it relevant to conservation.Ann Rev Ecol Evol Syst36267294

31. Niklaus PA, Wardle DA, Tate KR (2006) Effects of plant species diversity and composition on nitrogen cycling and the trace gas balance of soils. Plant Soil 282: 83–98.PA NiklausDA WardleKR Tate2006Effects of plant species diversity and composition on nitrogen cycling and the trace gas balance of soils.Plant Soil2828398

32. Oelmann Y, Wilcke W, Temperton VM, Buchmann N, Roscher C, et al. (2007) Soil and plant nitrogen pools as related to plant diversity in an experimental grassland. Soil Sci Soc Am J 71: 720–729.Y. OelmannW. WilckeVM TempertonN. BuchmannC. Roscher2007Soil and plant nitrogen pools as related to plant diversity in an experimental grassland.Soil Sci Soc Am J71720729

33. Delta-T Devices Ltd (1999) Theta Probe Soil Moisture Sensor Type ML2x User Manual. Cambridge, England: Delta-T Devices, Inc. Delta-T Devices Ltd1999Theta Probe Soil Moisture Sensor Type ML2x User Manual.Cambridge, EnglandDelta-T Devices, Inc

34. Duelli P, Obrist MK, Schmatz DR (1999) Biodiversity evaluation in agricultural landscapes: above-ground insects. Agr Ecosyst Environ 74: 33–64.P. DuelliMK ObristDR Schmatz1999Biodiversity evaluation in agricultural landscapes: above-ground insects.Agr Ecosyst Environ743364

35. Roscher C, Bessler H, Oelmann Y, Engels C, Wilcke W, et al. (2009) Resources, recruitment limitation and invader species identity determine pattern of spontaneous invasion in experimental grasslands. J Ecol 97: 32–47.C. RoscherH. BesslerY. OelmannC. EngelsW. Wilcke2009Resources, recruitment limitation and invader species identity determine pattern of spontaneous invasion in experimental grasslands.J Ecol973247

36. Eisenhauer N, Milcu A, Nitschke N, Sabais ACW, Scherber C, et al. (2009) Earthworm and below-ground competition effects on plant productivity in a plant diversity gradient. Oecologia 161: 291–301.N. EisenhauerA. MilcuN. NitschkeACW SabaisC. Scherber2009Earthworm and below-ground competition effects on plant productivity in a plant diversity gradient.Oecologia161291301

37. Legendre P, Legendre L (1998) Numerical Ecology. Amsterdam: Elsevier. 853 p.P. LegendreL. Legendre1998Numerical Ecology.AmsterdamElsevier853