Diversity and Productivity in a Long-Term Grassland Experiment

See allHide authors and affiliations

Science  26 Oct 2001:
Vol. 294, Issue 5543, pp. 843-845
DOI: 10.1126/science.1060391


Plant diversity and niche complementarity had progressively stronger effects on ecosystem functioning during a 7-year experiment, with 16-species plots attaining 2.7 times greater biomass than monocultures. Diversity effects were neither transients nor explained solely by a few productive or unviable species. Rather, many higher-diversity plots outperformed the best monoculture. These results help resolve debate over biodiversity and ecosystem functioning, show effects at higher than expected diversity levels, and demonstrate, for these ecosystems, that even the best-chosen monocultures cannot achieve greater productivity or carbon stores than higher-diversity sites.

Recent demonstrations that greater plant diversity can lead to greater productivity (1–5) have generated considerable debate (6–14). This has been fueled by uncertainty about which of many alternative hypotheses (6, 9, 15–20) is operating in nature and about the number of species required to maintain ecosystem functioning (6, 7, 12,14). We report results of a long-term experiment that allows tests of these alternative hypotheses.

It has been hypothesized that productivity may be greater at higher diversity because of “niche complementarity” among particular combinations of species and the greater chance of occurrence of such combinations at higher diversity (15–20). Niche complementarity, which results from interspecific differences in resource requirements and in spatial and temporal resource and habitat use, or from positive interactions (21), is predicted to allow stable multispecies coexistence and sustainably greater productivity at higher diversity (17, 18). Alternatively, it has been hypothesized that reported diversity effects might be short-lived transients caused solely by the presence of some species with high growth rates (6); be experimental artifacts resulting solely from species pools containing some low-viability species (6); or result from the most productive species being the best competitor (6,9, 17). These “sampling effects” all result from the greater chance of any given species being present at higher diversity and from dynamics that cause a single species to dominate and determine ecosystem functioning (6, 9,17, 19).

Sampling and complementarity have different signatures (6, 9, 12,17–20). Sampling effects limit the maximal productivity of higher-diversity plots to that of the best monoculture, giving an upper bound of variation in community performance that is independent of diversity. With niche complementarity, the upper bound increases with diversity because no monoculture is as productive as some combinations of two species and no combination of N species is as productive as some combinations of N + 1 species.

In a 7-year experiment [(4), supplement A (22)], we controlled one component of diversity, the number of plant species, in 168 plots, each 9 m by 9 m. We seeded the plots, in May 1994, to have 1, 2, 4, 8, or 16 species, with 39, 35, 29, 30, and 35 replicates, respectively. The species composition of each plot was chosen by random draw from a pool of 18 grassland perennials that included four C4 (warm-season) grasses, four C3 (cool-season) grasses, four legumes, four nonlegume forbs, and two woody species. All species occurred in monoculture, and all but three were in at least two monoculture plots, allowing comparison of responses of each species in monoculture to higher-diversity combinations of these same species. We do not use 76 additional plots that had functional group compositions drawn from an augmented species pool or 46 plots planted to 32 species (4) because the additional species were not grown in monoculture, and combining results from different species pools could introduce bias. We focus analyses on species number, because it was directly controlled, and functional group composition (23), because of its hypothesized importance. Other measures of diversity, including number of species per functional group and the presence or absence of species or functional groups, are highly correlated with species number and show similar responses.

In our grasslands, plant aboveground living biomass, because it is all produced within a growing season, is an index of primary productivity. In contrast, total biomass (aboveground plus belowground plant biomass) measures carbon accumulated in living tissues. Both aboveground and total biomass increased highly significantly with species number each year (Fig. 1, A and B, and Table 1), and functional group composition explained a highly significant amount of the residual variation (Table 1). Moreover, when the effect of each variable was determined after controlling for effects of the other (type III regressions), effects of functional group composition predominated in the early years [as for (24, 25)], but species number had highly significant positive effects on both aboveground and total biomass by 1999 and 2000, showing the simultaneous importance of species number and functional group composition in the long-term (Table 1).

Figure 1

The dependence of (A) plant aboveground biomass and (B) total biomass (aboveground plus belowground living plant mass) on the number of planted species. Data are shown as the mean ± SE. (C) The relation between the number of species planted in a plot and the actual number (mean ± SE) of planted species visually observed in a 2 m−2 area of each plot. (D) The percentage of all plots of a given planted diversity level, on average for 1999 and 2000 combined, or on average for 1997, that had greater biomass than the single monoculture plot with the greatest biomass.

Table 1

Analyses, using general linear models, of the effects of number of planted species (continuous variable; entered first using type I SS) and of functional group composition (categorical variable; entered second) on total biomass and on aboveground biomass, showing results for each year when measured. N = 168. Overall model df = 28 and error df = 139. Species number df = 1, composition df = 27. The last columns show type III effects of species number (entered second, after functional group composition).

View this table:

The initial saturating dependence of aboveground and total biomass on species number (Fig. 1, A and B) became, by 2000, a linear increase for species number ≥ 2. In 2000, 16-species plots had 22% greater aboveground biomass total and 27% greater total biomass than 8-species plots (both differences significant; t tests:P = 0.018, P = 0.002, respectively). The dependence of biomass on species number and functional group composition became progressively stronger, explaining about one-third of variance in 1997 and two-thirds in 2000 (Table 1). This strengthening of the effect of diversity and the increasingly steep and linear trends (Fig. 1, A and B) fail to support the hypothesis (6) that diversity effects were short-lived transients. Comparable and significant (P < 0.01) dependences of total and aboveground biomass on diversity and composition were observed when analyses used the actual number of planted species observed in each plot (Fig. 1C) or the Shannon diversity index [supplement B (22)].

We tested the low-viability sampling hypothesis by identifying the five species that attained least total biomass in monoculture in 2000 and excluding from analysis plots containing any combinations of just these species. Total biomass was still significantly dependent on species number and functional group composition in the remaining 131 plots [general linear model (GLM) type III regression:F OVERALL = 8.39, P < 0.001; F DIV = 11.6,F COMP = 4.36, P < 0.001 for each]. Similar results occurred when we excluded from analysis of aboveground biomass plots containing any combinations of the five species with least aboveground biomass in monoculture (GLM type III regression:F OVERALL = 7.08, P < 0.001; F DIV = 10.6, P = 0.0014; F COMP = 3.84,P < 0.001). In another analysis, we excluded the 30 plots with the lowest total biomass in 2000 (total biomass < 400 g m−2). Species number and composition still had highly significant effects on total biomass in the remaining plots (GLM type III regression: F OVERALL = 5.05, P < 0.001; F DIV = 12.3, P < 0.001;F COMP = 2.82, P < 0.001). Similarly, when we excluded the 31 plots with lowest aboveground biomass (<100 g m−2) in 2000, effects of composition and species number remained highly significant (GLM type III regression: F OVERALL = 3.77,P < 0.001; F DIV = 10.5,P < 0.0016; F COMP = 2.28,P = 0.0015). Similar results occurred with lower and higher cutoffs (including 50% higher) for aboveground and total biomass. In total, the dependence of biomass on species number and composition was not explained solely by sampling effects for a species pool containing some poorly performing species.

We tested the sampling hypothesis that the most productive species determined the effects of diversity (6,9, 17) by retaining in analyses of year 2000 results only plots containing at least one of the nine species with the highest monoculture total biomass in 2000. Total biomass remained significantly dependent on species number and functional group composition in these 145 plots, and in the subset of 95 plots that contained at least two of these nine species [type III regressions; supplement C (22)]. Similar results occurred for aboveground biomass in 2000 [supplement C (22)]. These analyses fail to support the sampling hypothesis. Another test comes from examining performance of higher-diversity plots relative to the best monoculture (Fig. 2). In 1999 and 2000, many higher-diversity plots had greater aboveground and total biomass than the single best-performing monoculture (Fig. 2). The percentage of such plots was an increasing function of diversity, on average for 1999 and 2000, with about half of the 16-species plots having greater aboveground or total biomass than the best monocultures (Fig. 1D). The strength and repeatability of this increasing upper bound in both aboveground and total biomass support the importance of niche effects and refute the hypothesis that sampling effects were the sole explanation for the long-term effects of diversity. The coexistence of most species, with about 12 planted species per 2 m2persisting in each 16-species plot (Fig. 1C), further supports niche complementarity. However, in earlier years, such as 1997, few high-diversity plots had greater biomass than top monocultures, and the percentage was independent of species number (Fig. 1D), which is consistent with the sampling hypothesis (6,7, 12).

Figure 2

The dependence of aboveground (A andB) and of total (C and D) biomass of each plot on planted species number for 1999 and 2000. The broken line shows the biomass of the top monoculture for a given year. The solid line is a regression of biomass on the logarithm of species number. Logarithm of species number was used in the figure because it gave slightly better fits, but was not used in Table 1 because it often gave slightly lower R 2 values than species number.

The increasing importance of complementarity and the increasingly linear effects of species number raise another question. Did complementarity occur among most species—i.e., did most species contribute to increasing community biomass—or is there a smaller set of species with complementary interactions, with this set being increasingly likely to co-occur at higher diversity (19)?

We used analysis of variance (ANOVA) to determine the simultaneous effects of the presence or absence of each species (entered as main effects) on aboveground or total biomass (one test per year, with 4 years for total biomass, 5 years for aboveground biomass). Three or four species had significant (P < 0.05) positive effects on aboveground or total biomass in most years. Among legumes, Lupinus perennis had significant effects in all nine tests, Lespedeza capitata in six tests, andPetalostemum purpureum in two tests. Schizachyrium scoparium and Sorghastrum nutans, both C4 grasses, were significant in five tests each. These are five of the six most abundant species in mixtures. A rarer forb also had a significant effect. Similarly, when plots were characterized by the presence or absence of functional groups in ANOVAs, in 2000 there were significant positive effects of legumes (P < 0.001), forbs (P < 0.05), and C4 grasses (P < 0.01) on aboveground biomass, and significant positive effects of legumes (P < 0.001) and C4 grasses (P < 0.001) on total biomass. For aboveground biomass, only the legume × C4 grass interaction was significant (P < 0.05), and it was positive. For total biomass, the legume × C4 grass interaction was marginally significant (P = 0.068) and biased toward positive (26), suggesting complementarity or facilitation among legumes and C4 grasses (4). However, even after controlling for the presence or absence of all functional groups, there were positive (P < 0.02) effects of species number on both aboveground and total biomass in 2000, indicating that biomass also depended on species number rather than on just the presence of functional groups.

Although these analyses suggest that the presence of about five dominant species might explain much of the effects of diversity, there may be small but additive effects of rarer species. To test this possibility, we ranked all species on the basis of their abundance (percentage cover) in 16-species plots in 2000, and created 17 new diversity indices. Each index states how many of the Nmost abundant species (N = 2, 3,…18) had been planted in each plot. We then determined which diversity index (log transformed) explained the most variance in aboveground or total biomass in 2000. Aboveground biomass was most dependent on how many of the 9 to 13 most abundant species were planted in each plot [supplement D (22)], showing that many rarer species contributed detectably to aboveground biomass. However, total biomass in 2000 was most dependent on how many of the four most abundant species were planted [supplement D (22)], likely because three of these four are C4 grasses, the species that accumulate the greatest root mass.

In summary, diversity effects were neither transients, nor explained in the long-term solely by other sampling-effect hypotheses, nor solely by the presence of legumes on a low-Nsoil. Rather, niche complementarity contributed significantly. Plant species number [as in (1–5,20, 21)] and functional group composition [as in (4, 5, 24, 25)] became simultaneously and approximately equally important in our long-term experiment. Compared with the average of the single best species in monoculture, our 16-species plots had 39% greater aboveground biomass and 42% greater total biomass on average for 1999 and 2000. Moreover, 16-species plots in 1999 and 2000 had 2.7 to 2.9 times greater aboveground and total biomass than the average for all species in monoculture (Fig. 1A). The nonsaturating increase in aboveground biomass with diversity likely reflects niche effects among about 9 to 13 species and their greater chance of co-occurrence at higher diversity (19), whereas such effects among about four species seem to account for total biomass responses.

The demonstration that diversity effects strengthened through time and were not the result solely of sampling effects or functional group composition should resolve aspects of the biodiversity debate (6–14). Moreover, our results suggest, for ecosystems in which niche complementarity occurs, that even with the wisest choices, monocultures may be less productive and accumulate less living carbon than many higher-diversity species combinations. Our results show that ecosystem processes are simultaneously influenced by diversity and composition, but long-term work in additional ecosystems is needed to determine the generality of our results, to better understand the effects of nonrandom community assembly and disassembly, and to better determine the implications of biodiversity for ecosystem management.

  • * To whom correspondence should be addressed. E-mail: tilman{at}


View Abstract

Navigate This Article