The effects of species pool,dispersal and competition on the diversity–productivity relationship

Aim The diversity–productivity relationship is a controversial issue in ecology. Diversity is sometimes seen to increase with productivity but a unimodal relationship has often been reported. Competitive exclusion was cited initially to account for the decrease of diversity at high productivity. Subsequently, the roles of evolutionary history (species pool size) and dispersal rate have been acknowledged. We explore how the effects of species pool, dispersal and competition combine to produce different diversity–productivity relationships. Methods We use a series of simulations with a spatially explicit, individual‐based model. Following empirical expectations, we used four scenarios to characterize species pool size along the productivity gradient (uniformly low and high, linear increase and unimodal). Similarly, the dispersal rate varied along the productivity gradient (uniformly low and high, and unimodal). We considered both neutral communities and communities with competitive exclusion. Results and main conclusions Our model predicts that competitive interactions will result in unimodal diversity–productivity relationships. The model often predicts unimodal patterns in neutral communities as well, although the decline in richness at high productivity is less than in competing communities. A positive diversity–productivity relationship is simulated for neutral communities when the species pool size increases with productivity and the dispersal rate is high. This scenario is probably more widespread in nature than the others since positive diversity–productivity relationships have been observed more frequently than previously expected, especially in the tropics and for woody species. Our simulated effects of species pool, dispersal and competition on diversity patterns can be linked to empirical observations to uncover mechanisms behind the diversity–productivity relationship.     

Effects of productivity and disturbance on species richness: a neutral model

Productivity and disturbance are major determinants of species diversity, and results from theoretical models predict that species richness should peak at intermediate levels of both factors. Such "unimodal" responses have been documented in many field and laboratory studies and have usually been attributed to differences among species in competitive ability and/or trade-offs between competitive ability and tolerance to disturbance. Here we show that most documented patterns of disturbance-richness and productivity-richness relationships, as well as the observed interactions between the two factors, can be explained by a simple neutral model where all species are ecologically identical and lack trade-offs in species characteristics. This finding suggests that current neutral theories can be extended to explain patterns of species responses to productivity and disturbance.     

Productivity–diversity patterns in arctic tundra vegetation

Productivity has long been argued to be a major driver of species richness patterns. In the present study we test alternative productivity–diversity hypotheses using vegetation data from the vast Eurasian tundra. The productivity–species pool hypothesis predicts positive relationships at both fine and coarse grain sizes, whereas the productivity–interaction hypothesis predicts unimodal patterns at fine grain size, and monotonic positive patterns at coarse grain size. We furthermore expect to find flatter positive (productivity–species pool hypothesis) or more strongly negative (productivity–interaction hypothesis) relationships for lichens and bryophytes than for vascular plants, because as a group, lichens and bryophytes are better adapted to extreme arctic conditions and more vulnerable to competition for light than the taller‐growing vascular plants. The normalised difference vegetation index (NDVI) was used as a proxy of productivity. The generally unimodal productivity–diversity patterns were most consistent with the productivity–interaction hypothesis. There was a general trend of decreasing species richness from moderately to maximally productive tundra, in agreement with an increasing importance of competitive interactions. High richness of vascular plants and lichens occurred in moderately low productive tundra areas, whereas that of bryophytes occurred in the least productive tundra habitats covered by this study. The fine and coarse grain richness trends were surprisingly uniform and no variation in beta diversity along the productivity gradient was seen for vascular plants or bryophytes. However, lichen beta diversity varied along the productivity gradient, probably reflecting their sensitivity to habitat conditions and biotic interactions. Overall, the results show evidence that productivity–diversity gradients exist in tundra and that these appear to be largely driven by competitive interactions. Our results also imply that climate warming‐driven increases in productivity will strongly affect arctic plant diversity patterns.     

Interactive effects of productivity and predation on zooplankton diversity

Recent studies suggest the necessity of understanding the interactive effects of predation and productivity on species coexistence and prey diversity. Models predict that coexistence of prey species with different competitive abilities can be achieved if inferior resource competitors are less susceptible to predation and if productivity and/or predation pressure are at intermediate levels. Hence, predator effects on prey diversity are predicted to be highly context dependent: enhancing diversity from low to intermediate levels of productivity or predation and reducing diversity of prey at high levels of productivity or predation. While several studies have examined the interactive effects of herbivory and productivity on primary producer diversity, experimental studies of such effects in predator‐prey systems are rare. We tested these predictions using an aquatic field mesocosm experiment in which initial density of the zooplankton predator Notonecta undulata and productivity were manipulated to test their interactive effects on diversity of seven zooplankton, cladoceran species that were common in surrounding ponds. Two productivity levels were imposed via phosphorus enrichment at levels comparable to low and intermediate levels found within neighboring natural ponds. We used open systems to allow for natural dispersal and behaviorally‐mediated numerical responses by the flight‐capable predator. Effects of predators on zooplankton diversity depended on productivity level. At low and high productivity, prey species richness declined while at high productivity it showed a unimodal relationship with increasing the predator density. Effects of treatments were weaker when using Pielou's evenness index or the inverse Simpson index as measures of prey diversity. Our findings are generally consistent with model predictions in which predators can facilitate prey coexistence and diversity at intermediate levels of productivity and predation intensity. Our work also shows that the functional form of the relationship between prey diversity and predation intensity can be complex and highly dependent on environmental context.     

Species' traits predict the effects of disturbance and productivity on diversity

Disturbance is an important factor influencing diversity patterns. Ecological theory predicts that diversity peaks at intermediate levels of disturbance, but this pattern is not present in a majority of empirical tests and can be influenced by the level of ecosystem productivity. We experimentally tested the effects of disturbance on diversity and show that species' autecological traits and community relations predicted species loss. We found that – alone or in concert – increasing disturbance intensity or frequency, or decreasing productivity, reduced diversity. Our species did not exhibit a clear competition-colonization trade-off, and intrinsic growth rate was a more important predictor of response to disturbance and productivity than measures of competitive ability. Furthermore, competitive ability was more important in predicting responses when, in addition to killing individuals, disturbance returned nutrients to the ecosystem. Our results demonstrate that species' traits can help resolve conflicting patterns in the response of diversity to disturbance and productivity.     

Comparing Neutral and Trade-off Community Models in Shaping the Community Biomass-Diversity Relationship Under Different Disturbance Levels

Among numerous mechanisms shaping the unimodal relationship between diversity and community biomass, the trade-off model of “CRS” theory is the most famous one. However, recent researches indicate that this relationship may also emerge under the neutral model where all species are identical with each other. By using an individual-based spatially-explicit model, we evaluated the underlying mechanisms shaping this curve for both models under different disturbance levels. We found unimodal relationships emerged for both models at low and medium disturbance levels; the richness for the trade-off community was lower than the neutral community for most of the environment severity levels, especially at the benign environment due to the strong competitive exclusions among species. Whereas under high disturbance level, the positive relationships emerged for both models; both communities had similar richness with their curves nearly overlapped with each other, that is, because the high disturbance intensity strongly decreased the competitive exclusions within the trade-off community. Our results indicate that although the underlying mechanisms are totally different, both models will produce the similar relationship between diversity and community biomass under different disturbance levels.     

Species pools and the "hump-back" model of plant species diversity: an empirical analysis at a relevant spatial scale

We investigated the relative roles of productivity, the species pool, and spatial habitat structure in determining local species richness (alpha diversity) of plant communities within a single, well-defined landscape unit, at spatial and ecological scales where the relationship between community productivity and species diversity often assumes a unimodal or "hump-back" form. At high levels of productivity, the decrease-phase of the unimodal model of the diversity-productivity relationship is typically explained as the dynamic outcome of increased competitive exclusion, but it may also be the passive consequence of a small pool of species possessing attributes necessary to competitively survive in high-fertility environments. We conducted statistical analyses of previously collected data to determine whether variations in local richness in the herbaceous vegetation of a Slovakian mountain valley were best explained by habitat productivity itself (which presumably leads to more intense competition) or by the sizes of the relevant community species pools. We also used measures of spatial habitat structure to investigate the extent to which habitat patchiness influenced patterns of species diversity. In the study system, both community biomass and size of the species pools contributed significantly to local species richness, but the positive effect of the species pools was about twice as important as the negative effect of biomass. The combined area of related associations (alliance area), association perimeter, and habitat patch geometry were all closely related to species pool size.     

Climate-induced temporal variation in the productivity–diversity relationship

Spatial scales are known to influence the form of the productivity–diversity relationship, but less attention has been given to the influence of temporal scales. Interannual climatic variation in semi-arid ponderosa pine–bunchgrass ecosystems causes significant year-to-year differences in herbaceous production. We hypothesized that unimodal (or 'hump-backed') relationships would be detected between herbaceous production and species richness in wet years, whereas positive logarithmic relationships would be detected in dry years. We analyzed nine years of herbaceous production and species richness data and used Akaike's information criterion (AICc) to determine the weight of evidence for each model in each year. As predicted, species richness exhibited a unimodal relationship to herbaceous production in wet years; however, richness exhibited a logarithmic relationship with herbaceous production in dry years. These results suggest that competitive exclusion occurred within this semi-arid plant community in years of high production when enough moisture was available to drive abundant plant growth. Thus, just as it is important to sample broad spatial variation in production to detect the full unimodal productivity–diversity relationship, it is also important to recognize that the full unimodal curve may be undetectable in less productive dry years in semi-arid ecosystems.     

The interactive effects of parasites, disturbance, and productivity on experimental adaptive radiations

Disturbance, productivity, and natural enemies are significant determinants of the evolution of diversity, but their interactive effect remains unresolved. We develop a simple, qualitative model assuming trade-offs between growth rate, competitive ability and parasite resistance, to address the interactive effects of these variables on the evolution of host diversity. Consistent with previous studies our model predicts maximum diversity at intermediate levels of disturbance and productivity in the absence of parasitism. However, parasites break down these unimodal diversity relationships with productivity and disturbance, as selection for parasite resistance reduces the importance of growth rate-competitive ability trade-offs. We tested these predictions using the bacterium Pseudomonas fluorescens, which undergoes an adaptive radiation into spatial niche specialists under laboratory conditions. This is the first study of adaptive radiation in response to experimental manipulation of the three-way interaction between productivity, disturbance, and natural enemies. As hypothesized, unimodal diversity relationships with disturbance and productivity were weakened or disappeared in the presence of parasitic phages. This was the result of phages increasing diversity at environmental extremes, by imposing selection for phage-resistant variants, but decreasing diversity in less stressful environments, probably through reductions in resource competition. Phages had a net effect of increasing host diversity. Parasites and other natural enemies are therefore likely to have a large effect in mitigating the influence of other environmental variables on the evolution and maintenance of diversity.     

What determines disturbance‐productivity‐diversity relationships? The effect of scale,species and environment on richness patterns in an Australian woodland

Much of the observed variation in relationships between diversity and disturbance or productivity may be attributed to scale, species characteristics, or environment. We used exclusion fences to create gradients of grazing (by native and introduced herbivores), cover, and standing crop in temperate Eucalypt woodlands. We investigated patterns of native, exotic, and total plant species richness at two scales (1 m² and 625 m²). Richness patterns were similar at both scales, though species richness at 1m² was more strongly affected by our grazing treatments. Season and rainfall explained more variation in richness than did surrogate measures of productivity or disturbance by herbivores. The richness-herbivory relationship depended strongly on rainfall, season, and species origin, and altering these factors produced the entire range of observed diversity-disturbance relationships. Richness-biomass and richness-cover relationships were consistently hump-shaped, and related to species origin with native richness negatively related and exotic richness positively related. The ability of weedy annuals to pre-empt space after death may have contributed to the observed unimodal responses.     

Plant species richness–productivity relationships in a low-productive boreal region

Local species richness–productivity (SR–P) relationship is usually reported as unimodal if long productivity gradients are sampled. However, it tends to be monotonically increasing in low-productive environments due to the decreasing part of the SR–P curve being truncated. Previous work indicated that this can hold true for forest herb layers, because of an upper bound on productivity caused mainly by canopy shading. Here, we ask whether the same pattern exists in a region with an upper bound on productivity caused by a harsh climate. We sampled herbaceous vegetation of boreal forests and grasslands in a low-productive region of central Yakutia (NE Siberia) with dry and winter-cool continental climate. We collected data on species composition, herb-layer productivity (aboveground herbaceous biomass), soil chemistry and light availability. We applied regression models to discriminate between monotonically increasing, decreasing and unimodal responses of herb-layer species richness to measured variables and analysed trends in the species-pool size and beta diversity along the productivity gradient. Our expectation of the monotonically increasing SR–P relationship was confirmed for neither forest herb layers nor grasslands. In the forest herb layers, no relationship was detected. In grasslands, the relationship was unimodal with species richness decline starting at much lower productivity levels than in more productive temperate grasslands. Potential causes for this decline are either limitation of local species richness by the species pool, which contains few species adapted to more productive habitats, or competitive exclusion, which can become an important control of species richness under lower levels of productivity than is the case in temperate grasslands.     

Coral species richness: ecological versus biogeographical influences

Species richness in communities varies with habitat area, productivity, disturbance level, intensity of species interactions, and regional/historical effects. All of these factors influence coral richness but their effects vary with spatial scale, position on the reef, and regional location. Species richness of corals along depth gradients shows a unimodal, hump-shaped curve that peaks at intermediate depths. Moreover, the peak of the curve is higher in regions with larger species pools. This “regional enrichment” of the local community appears in line transect samples as small as 10 m in length. The pattern suggests that ecological factors operating over scales of tens of meters and regional/historical factors operating over thousands of kilometers can both affect local richness. Regional factors probably include differences in speciation relative to extinction rates among regions and proximity of local sites to richness hotspots. Plausible factors operating at the local scale are species interactions, disturbance, and productivity which combine in different ways to produce the unimodal pattern. Shallow areas support few species because extinction rates are high due to frequent disturbance or because of environmental extremes. In addition, high productivity encourages rapid growth and thus the potential for intense interspecific competition. In areas where branching acroporids are abundant, exclusion by these dominant competitors is possible. Deep areas may be depauperate because few species can tolerate the low light levels found there. Areas of intermediate depth have the richest communities because they are open for colonization by many species and because extinction rates are low. Several theories may explain this “openness” and species persistence: 1. Occasional disturbance coupled with low growth rates results in glacially slow exclusion by the dominant competitor. 2. Aggregation of corals creates spatial variation in the intensity of competition and thus refuges from competition within a spatial landscape. Inferior competitors persist because they are superior at dispersal and refuge colonization. 3. Specialist predators focus on high-density juvenile populations near the parent, creating ecological space for colonization by non-prey. 4. Coral competitive abilities are roughly equal and recruitment into the community is a probabilistic event. The community thus exhibits random drift and exclusion is an extremely lengthy process. Based upon empirical evidence, these theories are listed in order of plausibility, but still need to be rigorously tested. Accepted: 9 September 1999     

What determines the relationship between plant diversity and habitat productivity?

The relationship between biodiversity and habitat productivity has been a fundamental topic in ecology. Although the relationship between these parameters may exhibit different shapes, the unimodal shape has been frequently encountered. The decrease in diversity at high productivity has usually been attributed to competitive exclusion. We suggest that evolutionary history and dispersal limitation may be even more important in shaping the diversity–productivity relationship. On a global scale, unimodal diversity–productivity relationships dominate in temperate regions, whereas positive relationships are more common in the tropics. This difference can be accounted for by contrasting evolutionary history. Temperate regions have smaller species pools for productive habitats since these habitats have been scarce historically for speciation, while the opposite is true for the tropics. In addition, dispersal within a region may limit diversity either due to the lack of dispersal syndromes at low productivity or the low number of diaspores at high productivity. Thereafter, biotic interactions (competition and facilitation) can shape the relationship. All these processes can act independently or concurrently. We recommend that the common approach to examining empirical diversity–environmental relationships should start with the role of large‐scale processes such as evolutionary history and dispersal limitation, followed by influences associated with ecological interactions.     

Grain size affects the relationship between species richness and above-ground biomass in semi-arid rangelands

ABSTRACT

Background: Discrepancies in the shape of the productivity–diversity relationship may arise from differences in spatial scale. We hypothesised that there is a grain size effect on the productivity–diversity relationship.

Aims: To determine the effect of three sampling grain sizes on the productivity–diversity relationship.

Methods: We applied generalised linear mixed effect models on community data from 735 vegetation plots in the Taleghan rangelands, Iran, sampled at three grain sizes (0.25, 1 and 2 m²) to ascertain plant productivity-diversity patterns, while accounting for the effects of site, plant community type, disturbance, and life form.

Results: Overall, relationships between biomass and plant species richness were unimodal at grain sizes of 0.25 and 1 m², and asymptotical at 2 m². The spurious occurrence of a single large shrub may overwhelm a small-sized sampling unit, resulting in a high estimate of the sample’s biomass relative to species richness. However, the relationship between biomass and species richness at larger grain sizes is more likely to reach an asymptote.

Conclusions: Shrubs are partly responsible for driving the relationship between plant biomass and species richness. Given that the frequency of shrubs is highly variable between small plots but not so in large plots, their presence may result in unimodal productivity–diversity relationships at small but not at large grain sizes.     

Dark diversity reveals importance of biotic resources and competition for plant diversity across habitats

Species richness is the most commonly used metric to quantify biodiversity. However, examining dark diversity, the group of missing species which can potentially inhabit a site, can provide a more thorough understanding of the processes influencing observed biodiversity and help evaluate the restoration potential of local habitats. So far, dark diversity has mainly been studied for specific habitats or large‐scale landscapes, while less attention has been given to variation across broad environmental gradients or as a result of local conditions and biotic interactions. In this study, we investigate the importance of local environmental conditions in determining dark diversity and observed richness in plant communities across broad environmental gradients. Using the ecospace concept, we investigate how these biodiversity measures relate to abiotic gradients (defined as position), availability of biotic resources (defined as expansion), spatiotemporal extent of habitats (defined as continuity), and species interactions through competition. Position variables were important for both observed diversity and dark diversity, some with quadratic relationships, for example, plant richness showing a unimodal response to soil fertility corresponding to the intermediate productivity hypothesis. Interspecific competition represented by community mean Grime C had a negative effect on plant species richness. Besides position‐related variables, organic carbon was the most important variable for dark diversity, indicating that in late‐succession habitats such as forests and shrubs, dark diversity is generally low. The importance of highly competitive species indicates that intermediate disturbance, such as grazing, may facilitate higher species richness and lower dark diversity.     

Dominance–diversity relationships in ant communities differ with invasion

The relationship between levels of dominance and species richness is highly contentious, especially in ant communities. The dominance‐impoverishment rule states that high levels of dominance only occur in species‐poor communities, but there appear to be many cases of high levels of dominance in highly diverse communities. The extent to which dominant species limit local richness through competitive exclusion remains unclear, but such exclusion appears more apparent for non‐native rather than native dominant species. Here we perform the first global analysis of the relationship between behavioral dominance and species richness. We used data from 1,293 local assemblages of ground‐dwelling ants distributed across five continents to document the generality of the dominance‐impoverishment rule, and to identify the biotic and abiotic conditions under which it does and does not apply. We found that the behavioral dominance–diversity relationship varies greatly, and depends on whether dominant species are native or non‐native, whether dominance is considered as occurrence or relative abundance, and on variation in mean annual temperature. There were declines in diversity with increasing dominance in invaded communities, but diversity increased with increasing dominance in native communities. These patterns occur along the global temperature gradient. However, positive and negative relationships are strongest in the hottest sites. We also found that climate regulates the degree of behavioral dominance, but differently from how it shapes species richness. Our findings imply that, despite strong competitive interactions among ants, competitive exclusion is not a major driver of local richness in native ant communities. Although the dominance‐impoverishment rule applies to invaded communities, we propose an alternative dominance‐diversification rule for native communities.     

Productivity-diversity relationships in lake plankton communities

One of the most intriguing environmental gradients connected with variation in diversity is ecosystem productivity. The role of diversity in ecosystems is pivotal, because species richness can be both a cause and a consequence of primary production. However, the mechanisms behind the varying productivity-diversity relationships (PDR) remain poorly understood. Moreover, large-scale studies on PDR across taxa are urgently needed. Here, we examined the relationships between resource supply and phyto-, bacterio-, and zooplankton richness in 100 small boreal lakes. We studied the PDR locally within the drainage systems and regionally across the systems. Second, we studied the relationships between resource availability, species richness, biomass and resource ratio (N:P) in phytoplankton communities using Structural Equation Modeling (SEM) for testing the multivariate hypothesis of PDR. At the local scale, the PDR showed variable patterns ranging from positive linear and unimodal to negative linear relationships for all planktonic groups. At the regional scale, PDRs were significantly linear and positive for phyto- and zooplankton. Phytoplankton richness and the amount of chlorophyll a showed a positive linear relationship indicating that communities consisting of higher number of species were able to produce higher levels of biomass. According to the SEM, phytoplankton biomass was largely related to resource availability, yet there was a pathway via community richness. Finally, we found that species richness at all trophic levels was correlated with several environmental factors, and was also related to richness at the other trophic levels. This study showed that the PDRs in freshwaters show scale-dependency. We also documented that the PDR complies with the multivariate model showing that plant biomass is not mirroring merely the resource availability, but is also influenced by richness. This highlights the need for conserving diversity in order to maintain ecosystem processes in freshwaters.     

Biodiversity–productivity relationship in ponds: Community and metacommunity patterns along time and environmental gradients

Primary production correlates with diversity in various ways. These patterns may result from the interaction of various mechanisms related to the environmental context and the spatial and temporal scale of analysis. However, empirical evidence on diversity‐productivity patterns typically considers single temporal and spatial scales, and does not include the effect of environmental variables. In a metacommunity of macrophytes in ephemeral ponds, we analysed the diversity‐productivity relationship patterns in the field, the importance of the environmental variables of pond size and heterogeneity on such relationship, and the variation of these patterns at local (community level) and landscape scales (metacommunity level) across 52 ponds on twelve occasions, over five years (2005–2009). Combining all sampling dates, there were 377 ponds and 1954 sample‐unit observations. Vegetation biomass was used as a proxy for productivity, and biodiversity was represented by species richness, evenness, and their interaction. Environmental variables comprised pond area, depth and internal heterogeneity. Productivity and species richness were not directly related at the metacommunity level, and were positively related at the community level. Taking environmental variables into account revealed positive species richness‐productivity relationships at the metacommunity level and positive quadratic relationships at the community level. Productivity showed both positive and negative linear and nonlinear relationships with the size and heterogeneity of ponds. We found a weak relationship between productivity and evenness. The identity of variables associated with productivity changed between spatial scales and through time. The pattern of relationships between productivity and diversity depends on spatial scale and environmental context, and changes idiosyncratically through time within the same ecosystem. Thus, the diversity‐productivity relationship is not only a property of the study system, but also a consequence of environmental variations and the temporal and spatial scale of analysis.     

Burial disturbance leads to facilitation among coastal dune plants

There is growing evidence that interactions among plants can be facilitative as well as competitive, but knowledge of how disturbances influence these interactions and how they vary with species diversity is lacking. We manipulated plant density, species diversity (richness), and a burial disturbance in a controlled, complete factorial experiment to test theories about the relationships among species interactions, disturbance, and richness. The hypotheses tested were 1) burial disturbance reduces plant performance at all levels of density and richness, 2) burial disturbance can cause net plant interactions to become more facilitative, and 3) facilitation increases with species richness. Burial decreased plant survival by 60% and biomass by 50%, supporting the hypothesis that burial reduces plant performance. In the control (unburied) treatment, there was no difference in proportion survival or per plant biomass between low and high density plots, meaning that neither competition nor facilitation was detected. In the buried treatment, however, high density plots had significantly greater survival and greater per plant biomass than the low density plots, indicating net facilitative interactions. Thus facilitation occurred in the buried treatment and not in the unburied control plots, supporting the hypothesis that facilitation increases with increasing disturbance severity. The hypothesis that facilitation increases with increasing species richness was not supported. Richness did not affect survival or biomass, and there was no richness by burial treatment interaction, indicating that richness did not influence the response of the community to burial. The influence of the disturbance on plant interactions was thus consistent across levels of richness, increasing the generality of the relationship between disturbance and facilitation.     

A long-term experimental test of the dynamic equilibrium model of species diversity

The dynamic equilibrium model of species diversity predicts that ecosystem productivity interacts with disturbance to determine how many species coexist. However, a robust test of this model requires manipulations of productivity and disturbance over a sufficient timescale to allow competitive exclusion, and such long-term experimental tests of this hypothesis are rare. Here we use long-term (27 years), large-scale (8 × 50-m plots), factorial manipulations of soil resource availability and sheep grazing intensity (disturbance) in grasslands to test the dynamic equilibrium model. As predicted by the model, increased productivity not only reduced plant species richness, but also moderated the effects of grazing intensity, shifting them from negative to neutral with increasing productivity. Reductions in species richness with productivity were associated with dominance by faster growing (i.e. high specific leaf area) and taller plants. Conversely, grazing favoured shorter plants and this effect became stronger with greater productivity, consistent with the view that grazing can lead to weaker asymmetric competition for light. Our study shows that the dynamic equilibrium model can help to explain changes in plant species richness following long-term increases in soil resource availability and grazing pressure, two fundamental drivers of change in grasslands worldwide.