Reintroducing bison results in long-running and resilient increases in grassland diversity.

The widespread extirpation of megafauna may have destabilized ecosystems and altered biodiversity globally. Most megafauna extinctions occurred before the modern record, leaving it unclear how their loss impacts current biodiversity. We report the long-term effects of reintroducing plains bison (Bison bison) in a tallgrass prairie versus two land uses that commonly occur in many North American grasslands: 1) no grazing and 2) intensive growing-season grazing by domesticated cattle (Bos taurus). Compared to ungrazed areas, reintroducing bison increased native plant species richness by 103% at local scales (10 m2) and 86% at the catchment scale. Gains in richness continued for 29 y and were resilient to the most extreme drought in four decades. These gains are now among the largest recorded increases in species richness due to grazing in grasslands globally. Grazing by domestic cattle also increased native plant species richness, but by less than half as much as bison. This study indicates that some ecosystems maintain a latent potential for increased native plant species richness following the reintroduction of native herbivores, which was unmatched by domesticated grazers. Native-grazer gains in richness were resilient to an extreme drought, a pressure likely to become more common under future global environmental change.

Grazing herbivores reduce herbaceous biomass and fire activity across African savannas.

Fire and herbivory interact to alter ecosystems and carbon cycling. In savannas, herbivores can reduce fire activity by removing grass biomass, but the size of these effects and what regulates them remain uncertain. To examine grazing effects on fuels and fire regimes across African savannas, we combined data from herbivore exclosure experiments with remotely sensed data on fire activity and herbivore density. We show that, broadly across African savannas, grazing herbivores substantially reduce both herbaceous biomass and fire activity. The size of these effects was strongly associated with grazing herbivore densities, and surprisingly, was mostly consistent across different environments. A one-zebra increase in herbivore biomass density (~100 kg/km2 of metabolic biomass) resulted in a ~53 kg/ha reduction in standing herbaceous biomass and a ~0.43 percentage point reduction in burned area. Our results indicate that fire models can be improved by incorporating grazing effects on grass biomass.

Accounting for herbaceous communities in process-based models will advance our understanding of "grassy" ecosystems.

Grassland and other herbaceous communities cover significant portions of Earth's terrestrial surface and provide many critical services, such as carbon sequestration, wildlife habitat, and food production. Forecasts of global change impacts on these services will require predictive tools, such as process-based dynamic vegetation models. Yet, model representation of herbaceous communities and ecosystems lags substantially behind that of tree communities and forests. The limited representation of herbaceous communities within models arises from two important knowledge gaps: first, our empirical understanding of the principles governing herbaceous vegetation dynamics is either incomplete or does not provide mechanistic information necessary to drive herbaceous community processes with models; second, current model structure and parameterization of grass and other herbaceous plant functional types limits the ability of models to predict outcomes of competition and growth for herbaceous vegetation. In this review, we provide direction for addressing these gaps by: (1) presenting a brief history of how vegetation dynamics have been developed and incorporated into earth system models, (2) reporting on a model simulation activity to evaluate current model capability to represent herbaceous vegetation dynamics and ecosystem function, and (3) detailing several ecological properties and phenomena that should be a focus for both empiricists and modelers to improve representation of herbaceous vegetation in models. Together, empiricists and modelers can improve representation of herbaceous ecosystem processes within models. In so doing, we will greatly enhance our ability to forecast future states of the earth system, which is of high importance given the rapid rate of environmental change on our planet.

Multiple global change drivers show independent, not interactive effects: a long-term case study in tallgrass prairie.

Ecosystems are faced with an onslaught of co-occurring global change drivers. While frequently studied independently, the effects of multiple global change drivers have the potential to be additive, antagonistic, or synergistic. Global warming, for example, may intensify the effects of more variable precipitation regimes with warmer temperatures increasing evapotranspiration and thereby amplifying the effect of already dry soils. Here, we present the long-term effects (11 years) of altered precipitation patterns (increased intra-annual variability in the growing season) and warming (1 °C year-round) on plant community composition and aboveground net primary productivity (ANPP), a key measure of ecosystem functioning in mesic tallgrass prairie. Based on past results, we expected that increased precipitation variability and warming would have additive effects on both community composition and ANPP. Increased precipitation variability altered plant community composition and increased richness, with no effect on ANPP. In contrast, warming decreased ANPP via reduction in grass stems and biomass but had no effect on the plant community. Contrary to expectations, across all measured variables, precipitation and warming treatments had no interactive effects. While treatment interactions did not occur, each treatment did individually impact a different component of the ecosystem (i.e., community vs. function). Thus, different aspects of the ecosystem may be sensitive to different global change drivers in mesic grassland ecosystems.

Richness, not evenness, varies across water availability gradients in grassy biomes on five continents.

We sought to understand the role that water availability (expressed as an aridity index) plays in determining regional and global patterns of richness and evenness, and in turn how these water availability-diversity relationships may result in different richness-evenness relationships at regional and global scales. We examined relationships between water availability, richness and evenness for eight grassy biomes spanning broad water availability gradients on five continents. Our study found that relationships between richness and water availability switched from positive for drier (South Africa, Tibet and USA) vs. negative for wetter (India) biomes, though were not significant for the remaining biomes. In contrast, only the India biome showed a significant relationship between water availability and evenness, which was negative. Globally, the richness-water availability relationship was hump-shaped, however, not significant for evenness. At the regional scale, a positive richness-evenness relationship was found for grassy biomes in India and Inner Mongolia, China. In contrast, this relationship was weakly concave-up globally. These results suggest that different, independent factors are determining patterns of species richness and evenness in grassy biomes, resulting in differing richness-evenness relationships at regional and global scales. As a consequence, richness and evenness may respond very differently across spatial gradients to anthropogenic changes, such as climate change.

Global change effects on plant communities are magnified by time and the number of global change factors imposed.

Global change drivers (GCDs) are expected to alter community structure and consequently, the services that ecosystems provide. Yet, few experimental investigations have examined effects of GCDs on plant community structure across multiple ecosystem types, and those that do exist present conflicting patterns. In an unprecedented global synthesis of over 100 experiments that manipulated factors linked to GCDs, we show that herbaceous plant community responses depend on experimental manipulation length and number of factors manipulated. We found that plant communities are fairly resistant to experimentally manipulated GCDs in the short term (<10 y). In contrast, long-term (≥10 y) experiments show increasing community divergence of treatments from control conditions. Surprisingly, these community responses occurred with similar frequency across the GCD types manipulated in our database. However, community responses were more common when 3 or more GCDs were simultaneously manipulated, suggesting the emergence of additive or synergistic effects of multiple drivers, particularly over long time periods. In half of the cases, GCD manipulations caused a difference in community composition without a corresponding species richness difference, indicating that species reordering or replacement is an important mechanism of community responses to GCDs and should be given greater consideration when examining consequences of GCDs for the biodiversity-ecosystem function relationship. Human activities are currently driving unparalleled global changes worldwide. Our analyses provide the most comprehensive evidence to date that these human activities may have widespread impacts on plant community composition globally, which will increase in frequency over time and be greater in areas where communities face multiple GCDs simultaneously.

Determinants of community compositional change are equally affected by global change.

Global change is impacting plant community composition, but the mechanisms underlying these changes are unclear. Using a dataset of 58 global change experiments, we tested the five fundamental mechanisms of community change: changes in evenness and richness, reordering, species gains and losses. We found 71% of communities were impacted by global change treatments, and 88% of communities that were exposed to two or more global change drivers were impacted. Further, all mechanisms of change were equally likely to be affected by global change treatments-species losses and changes in richness were just as common as species gains and reordering. We also found no evidence of a progression of community changes, for example, reordering and changes in evenness did not precede species gains and losses. We demonstrate that all processes underlying plant community composition changes are equally affected by treatments and often occur simultaneously, necessitating a wholistic approach to quantifying community changes.

Temporal variability in production is not consistently affected by global change drivers across herbaceous-dominated ecosystems.

Understanding how global change drivers (GCDs) affect aboveground net primary production (ANPP) through time is essential to predicting the reliability and maintenance of ecosystem function and services in the future. While GCDs, such as drought, warming and elevated nutrients, are known to affect mean ANPP, less is known about how they affect inter-annual variability in ANPP. We examined 27 global change experiments located in 11 different herbaceous ecosystems that varied in both abiotic and biotic conditions, to investigate changes in the mean and temporal variability of ANPP (measured as the coefficient of variation) in response to different GCD manipulations, including resource additions, warming, and irrigation. From this comprehensive data synthesis, we found that GCD treatments increased mean ANPP. However, GCD manipulations both increased and decreased temporal variability of ANPP (24% of comparisons), with no net effect overall. These inconsistent effects on temporal variation in ANPP can, in part, be attributed to site characteristics, such as mean annual precipitation and temperature as well as plant community evenness. For example, decreases in temporal variability in ANPP with the GCD treatments occurred in wetter and warmer sites with lower plant community evenness. Further, the addition of several nutrients simultaneously increased the sensitivity of ANPP to interannual variation in precipitation. Based on this analysis, we expect that GCDs will likely affect the magnitude more than the reliability over time of ecosystem production in the future.

Temporal heterogeneity increases with spatial heterogeneity in ecological communities.

Heterogeneity is increasingly recognized as a foundational characteristic of ecological systems. Under global change, understanding temporal community heterogeneity is necessary for predicting the stability of ecosystem functions and services. Indeed, spatial heterogeneity is commonly used in alternative stable state theory as a predictor of temporal heterogeneity and therefore an early indicator of regime shifts. To evaluate whether spatial heterogeneity in species composition is predictive of temporal heterogeneity in ecological communities, we analyzed 68 community data sets spanning freshwater and terrestrial systems where measures of species abundance were replicated over space and time. Of the 68 data sets, 55 (81%) had a weak to strongly positive relationship between spatial and temporal heterogeneity, while in the remaining communities the relationship was weak to strongly negative (19%). Based on a mixed model analysis, we found a significant but weak overall positive relationship between spatial and temporal heterogeneity across all data sets combined, and within aquatic and terrestrial data sets separately. In addition, lifespan and successional stage were negatively and positively related to temporal heterogeneity, respectively. We conclude that spatial heterogeneity may be a predictor of temporal heterogeneity in ecological communities, and that this relationship may be a general property of many terrestrial and aquatic communities.

Ambient changes exceed treatment effects on plant species abundance in global change experiments.

The responses of species to environmental changes will determine future community composition and ecosystem function. Many syntheses of global change experiments examine the magnitude of treatment effect sizes, but we lack an understanding of how plant responses to treatments compare to ongoing changes in the unmanipulated (ambient or background) system. We used a database of long-term global change studies manipulating CO2 , nutrients, water, and temperature to answer three questions: (a) How do changes in plant species abundance in ambient plots relate to those in treated plots? (b) How does the magnitude of ambient change in species-level abundance over time relate to responsiveness to global change treatments? (c) Does the direction of species-level responses to global change treatments differ from the direction of ambient change? We estimated temporal trends in plant abundance for 791 plant species in ambient and treated plots across 16 long-term global change experiments yielding 2,116 experiment-species-treatment combinations. Surprisingly, for most species (57%) the magnitude of ambient change was greater than the magnitude of treatment effects. However, the direction of ambient change, whether a species was increasing or decreasing in abundance under ambient conditions, had no bearing on the direction of treatment effects. Although ambient communities are inherently dynamic, there is now widespread evidence that anthropogenic drivers are directionally altering plant communities in many ecosystems. Thus, global change treatment effects must be interpreted in the context of plant species trajectories that are likely driven by ongoing environmental changes.

Ecological consequences of forest elephant declines for Afrotropical forests.

Poaching is rapidly extirpating African forest elephants (Loxodonta cyclotis) from most of their historical range, leaving vast areas of elephant-free tropical forest. Elephants are ecological engineers that create and maintain forest habitat; thus, their loss will have large consequences for the composition and structure of Afrotropical forests. Through a comprehensive literature review, we evaluated the roles of forest elephants in seed dispersal, nutrient recycling, and herbivory and physical damage to predict the cascading ecological effects of their population declines. Loss of seed dispersal by elephants will favor tree species dispersed abiotically and by smaller dispersal agents, and tree species composition will depend on the downstream effects of changes in elephant nutrient cycling and browsing. Loss of trampling and herbivory of seedlings and saplings will result in high tree density with release from browsing pressures. Diminished seed dispersal by elephants and high stem density are likely to reduce the recruitment of large trees and thus increase homogeneity of forest structure and decrease carbon stocks. The loss of ecological services by forest elephants likely means Central African forests will be more like Neotropical forests, from which megafauna were extirpated thousands of years ago. Without intervention, as much as 96% of Central African forests will have modified species composition and structure as elephants are compressed into remaining protected areas. Stopping elephant poaching is an urgent first step to mitigating these effects, but long-term conservation will require land-use planning that incorporates elephant habitat into forested landscapes that are being rapidly transformed by industrial agriculture and logging.

Asynchrony among local communities stabilises ecosystem function of metacommunities.

Temporal stability of ecosystem functioning increases the predictability and reliability of ecosystem services, and understanding the drivers of stability across spatial scales is important for land management and policy decisions. We used species-level abundance data from 62 plant communities across five continents to assess mechanisms of temporal stability across spatial scales. We assessed how asynchrony (i.e. different units responding dissimilarly through time) of species and local communities stabilised metacommunity ecosystem function. Asynchrony of species increased stability of local communities, and asynchrony among local communities enhanced metacommunity stability by a wide range of magnitudes (1-315%); this range was positively correlated with the size of the metacommunity. Additionally, asynchronous responses among local communities were linked with species' populations fluctuating asynchronously across space, perhaps stemming from physical and/or competitive differences among local communities. Accordingly, we suggest spatial heterogeneity should be a major focus for maintaining the stability of ecosystem services at larger spatial scales.

Deadwood stocks increase with selective logging and large tree frequency in Gabon.

Deadwood is a major component of aboveground biomass (AGB) in tropical forests and is important as habitat and for nutrient cycling and carbon storage. With deforestation and degradation taking place throughout the tropics, improved understanding of the magnitude and spatial variation in deadwood is vital for the development of regional and global carbon budgets. However, this potentially important carbon pool is poorly quantified in Afrotropical forests and the regional drivers of deadwood stocks are unknown. In the first large-scale study of deadwood in Central Africa, we quantified stocks in 47 forest sites across Gabon and evaluated the effects of disturbance (logging), forest structure variables (live AGB, wood density, abundance of large trees), and abiotic variables (temperature, precipitation, seasonality). Average deadwood stocks (measured as necromass, the biomass of deadwood) were 65 Mg ha-1 or 23% of live AGB. Deadwood stocks varied spatially with disturbance and forest structure, but not abiotic variables. Deadwood stocks increased significantly with logging (+38 Mg ha-1 ) and the abundance of large trees (+2.4 Mg ha-1 for every tree >60 cm dbh). Gabon holds 0.74 Pg C, or 21% of total aboveground carbon in deadwood, a threefold increase over previous estimates. Importantly, deadwood densities in Gabon are comparable to those in the Neotropics and respond similarly to logging, but represent a lower proportion of live AGB (median of 18% in Gabon compared to 26% in the Neotropics). In forest carbon accounting, necromass is often assumed to be a constant proportion (9%) of biomass, but in humid tropical forests this ratio varies from 2% in undisturbed forest to 300% in logged forest. Because logging significantly increases the deadwood carbon pool, estimates of tropical forest carbon should at a minimum use different ratios for logged (mean of 30%) and unlogged forests (mean of 18%).

Asymmetric responses of primary productivity to precipitation extremes: A synthesis of grassland precipitation manipulation experiments.

Climatic changes are altering Earth's hydrological cycle, resulting in altered precipitation amounts, increased interannual variability of precipitation, and more frequent extreme precipitation events. These trends will likely continue into the future, having substantial impacts on net primary productivity (NPP) and associated ecosystem services such as food production and carbon sequestration. Frequently, experimental manipulations of precipitation have linked altered precipitation regimes to changes in NPP. Yet, findings have been diverse and substantial uncertainty still surrounds generalities describing patterns of ecosystem sensitivity to altered precipitation. Additionally, we do not know whether previously observed correlations between NPP and precipitation remain accurate when precipitation changes become extreme. We synthesized results from 83 case studies of experimental precipitation manipulations in grasslands worldwide. We used meta-analytical techniques to search for generalities and asymmetries of aboveground NPP (ANPP) and belowground NPP (BNPP) responses to both the direction and magnitude of precipitation change. Sensitivity (i.e., productivity response standardized by the amount of precipitation change) of BNPP was similar under precipitation additions and reductions, but ANPP was more sensitive to precipitation additions than reductions; this was especially evident in drier ecosystems. Additionally, overall relationships between the magnitude of productivity responses and the magnitude of precipitation change were saturating in form. The saturating form of this relationship was likely driven by ANPP responses to very extreme precipitation increases, although there were limited studies imposing extreme precipitation change, and there was considerable variation among experiments. This highlights the importance of incorporating gradients of manipulations, ranging from extreme drought to extreme precipitation increases into future climate change experiments. Additionally, policy and land management decisions related to global change scenarios should consider how ANPP and BNPP responses may differ, and that ecosystem responses to extreme events might not be predicted from relationships found under moderate environmental changes.

Characterizing differences in precipitation regimes of extreme wet and dry years: implications for climate change experiments.

Climate change is intensifying the hydrologic cycle and is expected to increase the frequency of extreme wet and dry years. Beyond precipitation amount, extreme wet and dry years may differ in other ways, such as the number of precipitation events, event size, and the time between events. We assessed 1614 long-term (100 year) precipitation records from around the world to identify key attributes of precipitation regimes, besides amount, that distinguish statistically extreme wet from extreme dry years. In general, in regions where mean annual precipitation (MAP) exceeded 1000 mm, precipitation amounts in extreme wet and dry years differed from average years by ~40% and 30%, respectively. The magnitude of these deviations increased to >60% for dry years and to >150% for wet years in arid regions (MAP<500 mm). Extreme wet years were primarily distinguished from average and extreme dry years by the presence of multiple extreme (large) daily precipitation events (events >99th percentile of all events); these occurred twice as often in extreme wet years compared to average years. In contrast, these large precipitation events were rare in extreme dry years. Less important for distinguishing extreme wet from dry years were mean event size and frequency, or the number of dry days between events. However, extreme dry years were distinguished from average years by an increase in the number of dry days between events. These precipitation regime attributes consistently differed between extreme wet and dry years across 12 major terrestrial ecoregions from around the world, from deserts to the tropics. Thus, we recommend that climate change experiments and model simulations incorporate these differences in key precipitation regime attributes, as well as amount into treatments. This will allow experiments to more realistically simulate extreme precipitation years and more accurately assess the ecological consequences.

Interactive effects of grazing, drought, and fire on grassland plant communities in North America and South Africa.

Grazing, fire, and climate shape mesic grassland communities. With global change altering all three factors, understanding how grasslands respond to changes in these combined drivers may aid in projecting future changes in grassland ecosystems. We manipulated rainfall and simulated grazing (clipping) in two long-term fire experiments in mesic grasslands in North America (NA) and South Africa (SA). Despite their common drivers, grasslands in NA and SA differ in evolutionary history. Therefore, we expected community structure and production in NA and SA to respond differently to fire, grazing, and drought. Specifically, we hypothesized that NA plant community composition and production would be more responsive than the SA plant communities to changes in the drivers and their interactions, and that despite this expected stability of SA grasslands, drought would be the dominant factor controlling production, but grazing would play the primary role in determining community composition at both sites. Contrary to our hypothesis, NA and SA grasslands generally responded similarly to grazing, drought, and fire. Grazing increased diversity, decreased grass cover and production, and decreased belowground biomass at both sites. Drought alone minimally impacted plant community structure, and we saw similar treatment interactions at the two sites. Drought was not the primary driver of grassland productivity, but instead drought effects were similar to or less than grazing and fire. Even though these grasslands differed in evolutionary history, they responded similarly to our fire, grazing, and climate manipulations. Overall, we found community and ecosystem convergence in NA and SA grasslands. Grazing and fire are as important as climate in controlling mesic grassland ecosystems on both continents.

Plant community response to loss of large herbivores differs between North American and South African savanna grasslands.

Herbivory and fire shape plant community structure in grass-dominated ecosystems, but these disturbance regimes are being altered around the world. To assess the consequences of such alterations, we excluded large herbivores for seven years from mesic savanna grasslands sites burned at different frequencies in North America (Konza Prairie Biological Station, Kansas, USA) and South Africa (Kruger National Park). We hypothesized that the removal of a single grass-feeding herbivore from Konza would decrease plant community richness and shift community composition due to increased dominance by grasses. Similarly, we expected grass dominance to increase at Kruger when removing large herbivores, but because large herbivores are more diverse, targeting both grasses and forbs, at this study site, the changes due to herbivore removal would be muted. After seven years of large-herbivore exclusion, richness strongly decreased and community composition changed at Konza, whereas little change was evident at Kruger. We found that this divergence in response was largely due to differences in the traits and numbers of dominant grasses between the study sites rather than the predicted differences in herbivore assemblages. Thus, the diversity of large herbivores lost may be less important in determining plant community dynamics than the functional traits of the grasses that dominate mesic, disturbance-maintained savanna grasslands.

Loss of a large grazer impacts savanna grassland plant communities similarly in North America and South Africa.

Large herbivore grazing is a widespread disturbance in mesic savanna grasslands which increases herbaceous plant community richness and diversity. However, humans are modifying the impacts of grazing on these ecosystems by removing grazers. A more general understanding of how grazer loss will impact these ecosystems is hampered by differences in the diversity of large herbivore assemblages among savanna grasslands, which can affect the way that grazing influences plant communities. To avoid this we used two unique enclosures each containing a single, functionally similar large herbivore species. Specifically, we studied a bison (Bos bison) enclosure at Konza Prairie Biological Station, USA and an African buffalo (Syncerus caffer) enclosure in Kruger National Park, South Africa. Within these enclosures we erected exclosures in annually burned and unburned sites to determine how grazer loss would impact herbaceous plant communities, while controlling for potential fire-grazing interactions. At both sites, removal of the only grazer decreased grass and forb richness, evenness and diversity, over time. However, in Kruger these changes only occurred with burning. At both sites, changes in plant communities were driven by increased dominance with herbivore exclusion. At Konza, this was caused by increased abundance of one grass species, Andropogon gerardii, while at Kruger, three grasses, Themeda triandra, Panicum coloratum, and Digitaria eriantha increased in abundance.

Responses of C:N:P stoichiometric correlations among plants, soils and microorganisms to warming: A meta-analysis.

Plants, soils and microorganisms play important roles in maintaining stable terrestrial stoichiometry. Studying how nutrient balances of these biotic and abiotic players vary across temperature gradients is important when predicting ecosystem changes on a warming planet. The respective responses of plant, soil and microbial stoichiometric ratios to warming have been observed, however, whether and how the stoichiometric correlations among the three components shift under warming has not been clearly understood and identified. In the present study, we have performed a meta-analysis based on 600 case studies from 74 sites or locations to clarify whether and how warming affects plant, soil and microbial stoichiometry, respectively, and their correlations. Our results indicated that: (1) globally, plants had higher C:N and C:P values compared to soil and microbial pools, but their N:P distributions were similar; (2) warming did not significantly alter plant, soil and microbial C:N and C:P values, but had a noticeable effect on plant N:P ratios. When ecosystem types, duration and magnitude of warming were taken into account, there was an inconsistent and even inverse warming response in terms of the direction and magnitude of changes in the C:N:P ratios occurring among plants, soils and microorganisms; (3) despite various warming responses of the stoichiometric ratios detected separately for plants, soils and microorganisms, the stoichiometric correlations among all three parts remained constant even under different warming scenarios. Our study highlighted the complexity of the effect of warming on the C:N:P stoichiometry, as well as the absence and importance of simultaneous measurements of stoichiometric ratios across different components of terrestrial ecosystems, which should be urgently strengthened in future studies.

CoRRE Trait Data: A dataset of 17 categorical and continuous traits for 4079 grassland species worldwide.

In our changing world, understanding plant community responses to global change drivers is critical for predicting future ecosystem composition and function. Plant functional traits promise to be a key predictive tool for many ecosystems, including grasslands; however, their use requires both complete plant community and functional trait data. Yet, representation of these data in global databases is sparse, particularly beyond a handful of most used traits and common species. Here we present the CoRRE Trait Data, spanning 17 traits (9 categorical, 8 continuous) anticipated to predict species' responses to global change for 4,079 vascular plant species across 173 plant families present in 390 grassland experiments from around the world. The dataset contains complete categorical trait records for all 4,079 plant species obtained from a comprehensive literature search, as well as nearly complete coverage (99.97%) of imputed continuous trait values for a subset of 2,927 plant species. These data will shed light on mechanisms underlying population, community, and ecosystem responses to global change in grasslands worldwide.