Towards mechanistic integration of the causes and consequences of biodiversity

The global biodiversity crisis has stimulated decades of research on three themes: species coexistence, biodiversity-ecosystem functioning relationships (BEF), and biodiversity-ecosystem functional stability relationships (BEFS). However, studies on these themes are largely independent, creating barriers to an integrative understanding of the causes and consequences of biodiversity. Here we review recent progress towards mechanistic integration of coexistence, BEF, and BEFS. Mechanisms underlying the three themes can be linked in various ways, potentially creating either positive or negative relationships between them. That said, we generally expect positive associations between coexistence and BEF, and between BEF and BEFS. Our synthesis represents an initial step towards integrating causes and consequences of biodiversity; future developments should include more mechanistic approaches and broader ecological contexts.

Mechanistic Models of Trophic Interactions: Opportunities for Species Richness and Challenges for Modern Coexistence Theory

AbstractMany potential mechanisms promote species coexistence, but we know little about their relative importance. To compare multiple mechanisms, we modeled a two-trophic planktonic food web based on mechanistic species interactions and empirically measured species traits. We simulated thousands of possible communities under realistic and altered interaction strengths to assess the relative importance of three potential drivers of phytoplankton and zooplankton species richness: resource-mediated coexistence mechanisms, predator-prey interactions, and trait trade-offs. Next, we computed niche and fitness differences of competing zooplankton to obtain a deeper understanding of how these mechanisms determine species richness. We found that predator-prey interactions were the most important driver of phytoplankton and zooplankton species richness and that large zooplankton fitness differences were associated with low species richness, but zooplankton niche differences were not associated with species richness. However, for many communities we could not apply modern coexistence theory to compute niche and fitness differences of zooplankton because of conceptual issues with the invasion growth rates arising from trophic interactions. We therefore need to expand modern coexistence theory to fully investigate multitrophic-level communities.

Clarifying the effect of biodiversity on productivity in natural ecosystems with longitudinal data and methods for causal inference

Causal effects of biodiversity on ecosystem functions can be estimated using experimental or observational designs - designs that pose a tradeoff between drawing credible causal inferences from correlations and drawing generalizable inferences. Here, we develop a design that reduces this tradeoff and revisits the question of how plant species diversity affects productivity. Our design leverages longitudinal data from 43 grasslands in 11 countries and approaches borrowed from fields outside of ecology to draw causal inferences from observational data. Contrary to many prior studies, we estimate that increases in plot-level species richness caused productivity to decline: a 10% increase in richness decreased productivity by 2.4%, 95% CI [-4.1, -0.74]. This contradiction stems from two sources. First, prior observational studies incompletely control for confounding factors. Second, most experiments plant fewer rare and non-native species than exist in nature. Although increases in native, dominant species increased productivity, increases in rare and non-native species decreased productivity, making the average effect negative in our study. By reducing the tradeoff between experimental and observational designs, our study demonstrates how observational studies can complement prior ecological experiments and inform future ones.

Do plant-soil feedbacks promote coexistence in a sagebrush steppe?

Recent studies have shown the potential for negative plant-soil feedbacks (PSFs) to promote stable coexistence, but have not quantified the stabilizing effect relative to other coexistence mechanisms. We conducted a field experiment to test the role of PSFs in stabilizing coexistence among four dominant sagebrush steppe species that appear to coexist stably, based on previous work with observational data and models. We then integrated the effects of PSF treatments on focal species across germination, survival, and first-year growth. To contribute to stable coexistence, soil microbes should have host-specific effects that result in negative feedbacks. Over two replicated growing seasons, our experiments consistently showed that soil microbes have negative effects on plant growth, but these effects were rarely host-specific. The uncommon host-specific effects were mostly positive at the germination stage, and negative for growth. Integrated effects of PSF across early life-stage vital rates showed that PSF-mediated self-limitation occasionally had large effects on projected plant biomass, but occurred inconsistently between years. Our results suggest that while microbially-mediated PSF may not be a common mechanism of coexistence in this community, it may still affect the relative abundance of dominant plant species via changes in host fitness. Our work also serves as a blueprint for future investigations that aim to identify underlying processes and test alternative mechanisms to explain important patterns in community ecology.

Climate disequilibrium dominates uncertainty in long-term projections of primary productivity

Rapid climate change may exceed ecosystems' capacities to respond through processes including phenotypic plasticity, compositional turnover and evolutionary adaption. However, consequences of the resulting climate disequilibria for ecosystem functioning are rarely considered in projections of climate change impacts. Combining statistical models fit to historical climate data and remotely-sensed estimates of herbaceous net primary productivity with an ensemble of climate models, we demonstrate that assumptions concerning the magnitude of climate disequilibrium are a dominant source of uncertainty: models assuming maximum disequilibrium project widespread decreases in productivity in the western US by 2100, while models assuming minimal disequilibrium project productivity increases. Uncertainty related to climate disequilibrium is larger than uncertainties from variation among climate models or emissions pathways. A better understanding of processes that regulate climate disequilibria is essential for improving long-term projections of ecological responses and informing management to maintain ecosystem functioning at historical baselines.

Water availability dictates how plant traits predict demographic rates

A major goal in ecology is to make generalizable predictions of organism responses to environmental variation based on their traits. However, straightforward relationships between traits and fitness are rare and likely vary with environmental context. Characterizing how traits mediate demographic responses to the environment may enhance predictions of organism responses to global change. We synthesized 15 years of demographic data and species-level traits in a shortgrass steppe to determine whether the effects of leaf and root traits on growth and survival depend on seasonal water availability. We predicted that (1) species with drought-tolerant traits, such as lower leaf turgor loss point (TLP) and higher leaf and root dry matter content (LDMC and RDMC), would be more likely to survive and grow in drier years due to higher wilting resistance, (2) these traits would not predict fitness in wetter years, and (3) traits that more directly measure physiological mechanisms of water use such as TLP would best predict demographic responses. We found that graminoids with more negative TLP and higher LDMC and RDMC had higher survival rates in drier years. Forbs demonstrated similar yet more variable responses. Graminoids grew larger in wetter years, regardless of traits. However, in both wet and dry years, graminoids with more negative TLP and higher LDMC and RDMC grew larger than less negative TLP and low LDMC and RDMC species. Traits significantly mediated the impact of drought on survival, but not growth, suggesting survival could be a stronger driver of species' drought response in this system. TLP predicted survival in drier years, but easier-to-measure LDMC and RDMC were equal or better predictors. These results advance our understanding of the mechanisms by which drought drives population dynamics, and show that abiotic context determines how traits drive fitness.

Nutrient identity modifies the destabilising effects of eutrophication in grasslands

Nutrient enrichment can simultaneously increase and destabilise plant biomass production, with co-limitation by multiple nutrients potentially intensifying these effects. Here, we test how factorial additions of nitrogen (N), phosphorus (P) and potassium with essential nutrients (K+) affect the stability (mean/standard deviation) of aboveground biomass in 34 grasslands over 7 years. Destabilisation with fertilisation was prevalent but was driven by single nutrients, not synergistic nutrient interactions. On average, N-based treatments increased mean biomass production by 21-51% but increased its standard deviation by 40-68% and so consistently reduced stability. Adding P increased interannual variability and reduced stability without altering mean biomass, while K+ had no general effects. Declines in stability were largest in the most nutrient-limited grasslands, or where nutrients reduced species richness or intensified species synchrony. We show that nutrients can differentially impact the stability of biomass production, with N and P in particular disproportionately increasing its interannual variability.

Quadrat-based monitoring of desert grassland vegetation at the Jornada Experimental Range, New Mexico, 1915-2016

The data set covers a 101-yr period (1915-2016) of quadrat-based plant sampling at the Jornada Experimental Range in southern New Mexico. At each sampling event, a pantograph was used to record the location and perimeter of living plants within permanent quadrats. Basal area was recorded for perennial grass species, canopy cover area was recorded for shrub species, and all other perennial species were recorded as point data. The data set includes 122 1 × 1 m permanent quadrats, although not all quadrats were sampled in each year of the study and there is a gap in monitoring from 1980 to 1995. These data provide a unique opportunity to investigate changes in the plant community over 100 yr of variation in precipitation and other environmental conditions. We provide the following data and data formats: (1) the digitized maps in shapefile format; (2) a data table containing coordinates (x, y) of perennial species within quadrats, including cover area for grasses and shrubs; (3) a data table of counts of annual plant individuals per quadrat; (4) a species list indicating growth form and habit of recorded species; (5) a table of dates when each quadrat was sampled; (6) a table of the pasture each quadrat was located within (note that pasture boundaries have changed over time); (7) a table of depth to petrocalcic layer measurements taken at quadrat locations; (8) a table of particle size analysis of soil samples taken at quadrat locations; (9) a table of topographic characteristics of quadrat locations (e.g., concave or convex topography). Pantograph sampling is currently conducted at 5-yr intervals by USDA-ARS staff, and new data will be added periodically to the EDI Data Portal Repository (see section V.E.2). This information is released under the Creative Commons license-Attribution-CC BY and the consumer of these data is required to cite it appropriately in any publication that results from its use.

Biotic vs abiotic controls on temporal sensitivity of primary production to precipitation across North American drylands

Dryland net primary productivity (NPP) is sensitive to temporal variation in precipitation (PPT), but the magnitude of this 'temporal sensitivity' varies spatially. Hypotheses for spatial variation in temporal sensitivity have often emphasized abiotic factors, such as moisture limitation, while overlooking biotic factors, such as vegetation structure. We tested these hypotheses using spatiotemporal models fit to remote-sensing data sets to assess how vegetation structure and climate influence temporal sensitivity across five dryland ecoregions of the western USA. Temporal sensitivity was higher in locations and ecoregions dominated by herbaceous vegetation. By contrast, much less spatial variation in temporal sensitivity was explained by mean annual PPT. In fact, ecoregion-specific models showed inconsistent associations of sensitivity and PPT; whereas sensitivity decreased with increasing mean annual PPT in most ecoregions, it increased with mean annual PPT in the most arid ecoregion, the hot deserts. The strong, positive influence of herbaceous vegetation on temporal sensitivity indicates that herbaceous-dominated drylands will be particularly sensitive to future increases in precipitation variability and that dramatic changes in cover type caused by invasions or shrub encroachment will lead to changes in dryland NPP dynamics, perhaps independent of changes in precipitation.

Species loss due to nutrient addition increases with spatial scale in global grasslands

The effects of altered nutrient supplies and herbivore density on species diversity vary with spatial scale, because coexistence mechanisms are scale dependent. This scale dependence may alter the shape of the species-area relationship (SAR), which can be described by changes in species richness (S) as a power function of the sample area (A): S = cA^z , where c and z are constants. We analysed the effects of experimental manipulations of nutrient supply and herbivore density on species richness across a range of scales (0.01-75 m² ) at 30 grasslands in 10 countries. We found that nutrient addition reduced the number of species that could co-occur locally, indicated by the SAR intercepts (log c), but did not affect the SAR slopes (z). As a result, proportional species loss due to nutrient enrichment was largely unchanged across sampling scales, whereas total species loss increased over threefold across our range of sampling scales.

A practical guide to selecting models for exploration, inference, and prediction in ecology

Selecting among competing statistical models is a core challenge in science. However, the many possible approaches and techniques for model selection, and the conflicting recommendations for their use, can be confusing. We contend that much confusion surrounding statistical model selection results from failing to first clearly specify the purpose of the analysis. We argue that there are three distinct goals for statistical modeling in ecology: data exploration, inference, and prediction. Once the modeling goal is clearly articulated, an appropriate model selection procedure is easier to identify. We review model selection approaches and highlight their strengths and weaknesses relative to each of the three modeling goals. We then present examples of modeling for exploration, inference, and prediction using a time series of butterfly population counts. These show how a model selection approach flows naturally from the modeling goal, leading to different models selected for different purposes, even with exactly the same data set. This review illustrates best practices for ecologists and should serve as a reminder that statistical recipes cannot substitute for critical thinking or for the use of independent data to test hypotheses and validate predictions.

Increasing effects of chronic nutrient enrichment on plant diversity loss and ecosystem productivity over time

Human activities are enriching many of Earth's ecosystems with biologically limiting mineral nutrients such as nitrogen (N) and phosphorus (P). In grasslands, this enrichment generally reduces plant diversity and increases productivity. The widely demonstrated positive effect of diversity on productivity suggests a potential negative feedback, whereby nutrient-induced declines in diversity reduce the initial gains in productivity arising from nutrient enrichment. In addition, plant productivity and diversity can be inhibited by accumulations of dead biomass, which may be altered by nutrient enrichment. Over longer time frames, nutrient addition may increase soil fertility by increasing soil organic matter and nutrient pools. We examined the effects of 5-11 yr of nutrient addition at 47 grasslands in 12 countries. Nutrient enrichment increased aboveground live biomass and reduced plant diversity at nearly all sites, and these effects became stronger over time. We did not find evidence that nutrient-induced losses of diversity reduced the positive effects of nutrients on biomass; however, nutrient effects on live biomass increased more slowly at sites where litter was also increasing, regardless of plant diversity. This work suggests that short-term experiments may underestimate the long-term nutrient enrichment effects on global grassland ecosystems.

What processes must we understand to forecast regional-scale population dynamics?

An urgent challenge facing biologists is predicting the regional-scale population dynamics of species facing environmental change. Biologists suggest that we must move beyond predictions based on phenomenological models and instead base predictions on underlying processes. For example, population biologists, evolutionary biologists, community ecologists and ecophysiologists all argue that the respective processes they study are essential. Must our models include processes from all of these fields? We argue that answering this critical question is ultimately an empirical exercise requiring a substantial amount of data that have not been integrated for any system to date. To motivate and facilitate the necessary data collection and integration, we first review the potential importance of each mechanism for skilful prediction. We then develop a conceptual framework based on reaction norms, and propose a hierarchical Bayesian statistical framework to integrate processes affecting reaction norms at different scales. The ambitious research programme we advocate is rapidly becoming feasible due to novel collaborations, datasets and analytical tools.

Global impacts of fertilization and herbivore removal on soil net nitrogen mineralization are modulated by local climate and soil properties

Soil nitrogen (N) availability is critical for grassland functioning. However, human activities have increased the supply of biologically limiting nutrients, and changed the density and identity of mammalian herbivores. These anthropogenic changes may alter net soil N mineralization (soil net N_min ), that is, the net balance between N mineralization and immobilization, which could severely impact grassland structure and functioning. Yet, to date, little is known about how fertilization and herbivore removal individually, or jointly, affect soil net N_min across a wide range of grasslands that vary in soil and climatic properties. Here we collected data from 22 grasslands on five continents, all part of a globally replicated experiment, to assess how fertilization and herbivore removal affected potential (laboratory-based) and realized (field-based) soil net N_min . Herbivore removal in the absence of fertilization did not alter potential and realized soil net N_min . However, fertilization alone and in combination with herbivore removal consistently increased potential soil net N_min. Realized soil net N_min , in contrast, significantly decreased in fertilized plots where herbivores were removed. Treatment effects on potential and realized soil net N_min were contingent on site-specific soil and climatic properties. Fertilization effects on potential soil net N_min were larger at sites with higher mean annual precipitation (MAP) and temperature of the wettest quarter (T.q.wet). Reciprocally, realized soil net N_min declined most strongly with fertilization and herbivore removal at sites with lower MAP and higher T.q.wet. In summary, our findings show that anthropogenic nutrient enrichment, herbivore exclusion and alterations in future climatic conditions can negatively impact soil net N_min across global grasslands under realistic field conditions. This is an important context-dependent knowledge for grassland management worldwide.

Synchrony matters more than species richness in plant community stability at a global scale

The stability of ecological communities is critical for the stable provisioning of ecosystem services, such as food and forage production, carbon sequestration, and soil fertility. Greater biodiversity is expected to enhance stability across years by decreasing synchrony among species, but the drivers of stability in nature remain poorly resolved. Our analysis of time series from 79 datasets across the world showed that stability was associated more strongly with the degree of synchrony among dominant species than with species richness. The relatively weak influence of species richness is consistent with theory predicting that the effect of richness on stability weakens when synchrony is higher than expected under random fluctuations, which was the case in most communities. Land management, nutrient addition, and climate change treatments had relatively weak and varying effects on stability, modifying how species richness, synchrony, and stability interact. Our results demonstrate the prevalence of biotic drivers on ecosystem stability, with the potential for environmental drivers to alter the intricate relationship among richness, synchrony, and stability.

Technical Comment on Pande et al. (2020): Why invasion analysis is important for understanding coexistence

Pande et al. (2020) point out that persistence time can decrease even as invader growth rates (IGRs) increase, which potentially undermines modern coexistence theory. However, because persistence time increases rapidly with system size only when IGR > 0, to understand how any real community persists, we should first identify the mechanisms producing positive IGR.

Author Correction: Leaf nutrients, not specific leaf area, are consistent indicators of elevated nutrient inputs

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Nutrient availability controls the impact of mammalian herbivores on soil carbon and nitrogen pools in grasslands

Grasslands are subject to considerable alteration due to human activities globally, including widespread changes in populations and composition of large mammalian herbivores and elevated supply of nutrients. Grassland soils remain important reservoirs of carbon (C) and nitrogen (N). Herbivores may affect both C and N pools and these changes likely interact with increases in soil nutrient availability. Given the scale of grassland soil fluxes, such changes can have striking consequences for atmospheric C concentrations and the climate. Here, we use the Nutrient Network experiment to examine the responses of soil C and N pools to mammalian herbivore exclusion across 22 grasslands, under ambient and elevated nutrient availabilities (fertilized with NPK + micronutrients). We show that the impact of herbivore exclusion on soil C and N pools depends on fertilization. Under ambient nutrient conditions, we observed no effect of herbivore exclusion, but under elevated nutrient supply, pools are smaller upon herbivore exclusion. The highest mean soil C and N pools were found in grazed and fertilized plots. The decrease in soil C and N upon herbivore exclusion in combination with fertilization correlated with a decrease in aboveground plant biomass and microbial activity, indicating a reduced storage of organic matter and microbial residues as soil C and N. The response of soil C and N pools to herbivore exclusion was contingent on temperature - herbivores likely cause losses of C and N in colder sites and increases in warmer sites. Additionally, grasslands that contain mammalian herbivores have the potential to sequester more N under increased temperature variability and nutrient enrichment than ungrazed grasslands. Our study highlights the importance of conserving mammalian herbivore populations in grasslands worldwide. We need to incorporate local-scale herbivory, and its interaction with nutrient enrichment and climate, within global-scale models to better predict land-atmosphere interactions under future climate change.

Soil net nitrogen mineralisation across global grasslands

Soil nitrogen mineralisation (Nmin), the conversion of organic into inorganic N, is important for productivity and nutrient cycling. The balance between mineralisation and immobilisation (net Nmin) varies with soil properties and climate. However, because most global-scale assessments of net Nmin are laboratory-based, its regulation under field-conditions and implications for real-world soil functioning remain uncertain. Here, we explore the drivers of realised (field) and potential (laboratory) soil net Nmin across 30 grasslands worldwide. We find that realised Nmin is largely explained by temperature of the wettest quarter, microbial biomass, clay content and bulk density. Potential Nmin only weakly correlates with realised Nmin, but contributes to explain realised net Nmin when combined with soil and climatic variables. We provide novel insights of global realised soil net Nmin and show that potential soil net Nmin data available in the literature could be parameterised with soil and climate data to better predict realised Nmin.

Climate change, snow mold and the Bromus tectorum invasion: mixed evidence for release from cold weather pathogens

Climate change is reducing the depth and duration of winter snowpack, leading to dramatic changes in the soil environment with potentially important ecological consequences. Previous experiments in the Intermountain West of North America indicated that loss of snowpack increases survival and population growth rates of the invasive annual grass Bromus tectorum; however, the underlying mechanism is unknown. We hypothesized that reduced snowpack might promote B. tectorum population growth by decreasing damage from snow molds, a group of subnivean fungal pathogens. To test this hypothesis, we conducted greenhouse and field experiments to investigate the interaction between early snowmelt and either fungicide addition or snow mold infection of B. tectorum. The greenhouse experiment confirmed that the snow mold Microdochium nivale can cause mortality of B. tectorum seedlings. In the field experiment, early snowmelt and fungicide application both increased B. tectorum survival, but their effects did not interact, and snow mold inoculation had no effect on survival. We did find interactive effects of snowmelt and fungal treatments on B. tectorum seed production: with ambient snowpack, M. nivale inoculation reduced seed production and fungicide increased it, whereas in the early snowmelt treatment seed production was high regardless of fungal treatment. However, treatment effects on seed production did not translate directly to overall population growth, which did not respond to the snow melt by fungal treatment interaction. Based on our mixed results, the hypothesis that reduced snowpack may increase B. tectorum fitness by limiting the effects of plant pathogens deserves further investigation.

Sensitivity of global soil carbon stocks to combined nutrient enrichment

Soil stores approximately twice as much carbon as the atmosphere and fluctuations in the size of the soil carbon pool directly influence climate conditions. We used the Nutrient Network global change experiment to examine how anthropogenic nutrient enrichment might influence grassland soil carbon storage at a global scale. In isolation, enrichment of nitrogen and phosphorous had minimal impacts on soil carbon storage. However, when these nutrients were added in combination with potassium and micronutrients, soil carbon stocks changed considerably, with an average increase of 0.04 KgCm^-2  year^-1 (standard deviation 0.18 KgCm^-2  year^-1 ). These effects did not correlate with changes in primary productivity, suggesting that soil carbon decomposition may have been restricted. Although nutrient enrichment caused soil carbon gains most dry, sandy regions, considerable absolute losses of soil carbon may occur in high-latitude regions that store the majority of the world's soil carbon. These mechanistic insights into the sensitivity of grassland carbon stocks to nutrient enrichment can facilitate biochemical modelling efforts to project carbon cycling under future climate scenarios.

Environmental responses, not species interactions, determine synchrony of dominant species in semiarid grasslands

Temporal asynchrony among species helps diversity to stabilize ecosystem functioning, but identifying the mechanisms that determine synchrony remains a challenge. Here, we refine and test theory showing that synchrony depends on three factors: species responses to environmental variation, interspecific interactions, and demographic stochasticity. We then conduct simulation experiments with empirical population models to quantify the relative influence of these factors on the synchrony of dominant species in five semiarid grasslands. We found that the average synchrony of per capita growth rates, which can range from 0 (perfect asynchrony) to 1 (perfect synchrony), was higher when environmental variation was present (0.62) rather than absent (0.43). Removing interspecific interactions and demographic stochasticity had small effects on synchrony. For the dominant species in these plant communities, where species interactions and demographic stochasticity have little influence, synchrony reflects the covariance in species' responses to the environment.

Spatial heterogeneity in species composition constrains plant community responses to herbivory and fertilisation

Environmental change can result in substantial shifts in community composition. The associated immigration and extinction events are likely constrained by the spatial distribution of species. Still, studies on environmental change typically quantify biotic responses at single spatial (time series within a single plot) or temporal (spatial beta diversity at single time points) scales, ignoring their potential interdependence. Here, we use data from a global network of grassland experiments to determine how turnover responses to two major forms of environmental change - fertilisation and herbivore loss - are affected by species pool size and spatial compositional heterogeneity. Fertilisation led to higher rates of local extinction, whereas turnover in herbivore exclusion plots was driven by species replacement. Overall, sites with more spatially heterogeneous composition showed significantly higher rates of annual turnover, independent of species pool size and treatment. Taking into account spatial biodiversity aspects will therefore improve our understanding of consequences of global and anthropogenic change on community dynamics.

Competition and coexistence in plant communities: intraspecific competition is stronger than interspecific competition

Theory predicts that intraspecific competition should be stronger than interspecific competition for any pair of stably coexisting species, yet previous literature reviews found little support for this pattern. We screened over 5400 publications and identified 39 studies that quantified phenomenological intraspecific and interspecific interactions in terrestrial plant communities. Of the 67% of species pairs in which both intra- and interspecific effects were negative (competitive), intraspecific competition was, on average, four to five-fold stronger than interspecific competition. Of the remaining pairs, 93% featured intraspecific competition and interspecific facilitation, a situation that stabilises coexistence. The difference between intra- and interspecific effects tended to be larger in observational than experimental data sets, in field than greenhouse studies, and in studies that quantified population growth over the full life cycle rather than single fitness components. Our results imply that processes promoting stable coexistence at local scales are common and consequential across terrestrial plant communities.

Weak interspecific interactions in a sagebrush steppe? Conflicting evidence from observations and experiments

Stable coexistence requires intraspecific limitations to be stronger than interspecific limitations. The greater the difference between intra- and interspecific limitations, the more stable the coexistence, and the weaker the competitive release any species should experience following removal of competitors. We conducted a removal experiment to test whether a previously estimated model, showing surprisingly weak interspecific competition for four dominant species in a sagebrush steppe, accurately predicts competitive release. Our treatments were (1) removal of all perennial grasses and (2) removal of the dominant shrub, Artemisia tripartita. We regressed survival, growth, and recruitment on the locations, sizes, and species identities of neighboring plants, along with an indicator variable for removal treatment. If our "baseline" regression model, which accounts for local plant-plant interactions, accurately explains the observed responses to removals, then the removal coefficient should be non-significant. For survival, the removal coefficients were never significantly different from zero, and only A. tripartita showed a (negative) response to removals at the recruitment stage. For growth, the removal treatment effect was significant and positive for two species, Poa secunda and Pseudoroegneria spicata, indicating that the baseline model underestimated interspecific competition. For all three grass species, population models based on the vital rate regressions that included removal effects projected 1.4- to 3-fold increases in equilibrium population size relative to the baseline model (no removal effects). However, we found no evidence of higher response to removal in quadrats with higher pretreatment cover of A. tripartita, or by plants experiencing higher pre-treatment crowding by A. tripartita, raising questions about the mechanisms driving the positive response to removal. While our results show the value of combining observations with a simple removal experiment, more tightly controlled experiments focused on underlying mechanisms may be required to conclusively validate or reject predictions from phenomenological models.

The response of big sagebrush (Artemisia tridentata) to interannual climate variation changes across its range

Understanding how annual climate variation affects population growth rates across a species' range may help us anticipate the effects of climate change on species distribution and abundance. We predict that populations in warmer or wetter parts of a species' range should respond negatively to periods of above average temperature or precipitation, respectively, whereas populations in colder or drier areas should respond positively to periods of above average temperature or precipitation. To test this, we estimated the population sensitivity of a common shrub species, big sagebrush (Artemisia tridentata), to annual climate variation across its range. Our analysis includes 8,175 observations of year-to-year change in sagebrush cover or production from 131 monitoring sites in western North America. We coupled these observations with seasonal weather data for each site and analyzed the effects of spring through fall temperatures and fall through spring accumulated precipitation on annual changes in sagebrush abundance. Sensitivity to annual temperature variation supported our hypothesis: years with above average temperatures were beneficial to sagebrush in colder locations and detrimental to sagebrush in hotter locations. In contrast, sensitivity to precipitation did not change significantly across the distribution of sagebrush. This pattern of responses suggests that regional abundance of this species may be more limited by temperature than by precipitation. We also found important differences in how the ecologically distinct subspecies of sagebrush responded to the effects of precipitation and temperature. Our model predicts that a short-term temperature increase could produce an increase in sagebrush cover at the cold edge of its range and a decrease in cover at the warm edge of its range. This prediction is qualitatively consistent with predictions from species distribution models for sagebrush based on spatial occurrence data, but it provides new mechanistic insight and helps estimate how much and how fast sagebrush cover may change within its range.

Ecosystem functional response across precipitation extremes in a sagebrush steppe

Background

Precipitation is predicted to become more variable in the western United States, meaning years of above and below average precipitation will become more common. Periods of extreme precipitation are major drivers of interannual variability in ecosystem functioning in water limited communities, but how ecosystems respond to these extremes over the long-term may shift with precipitation means and variances. Long-term changes in ecosystem functional response could reflect compensatory changes in species composition or species reaching physiological thresholds at extreme precipitation levels.

Methods

We conducted a five year precipitation manipulation experiment in a sagebrush steppe ecosystem in Idaho, United States. We used drought and irrigation treatments (approximately 50% decrease/increase) to investigate whether ecosystem functional response remains consistent under sustained high or low precipitation. We recorded data on aboveground net primary productivity (ANPP), species abundance, and soil moisture. We fit a generalized linear mixed effects model to determine if the relationship between ANPP and soil moisture differed among treatments. We used nonmetric multidimensional scaling to quantify community composition over the five years.

Results

Ecosystem functional response, defined as the relationship between soil moisture and ANPP, was similar among irrigation and control treatments, but the drought treatment had a greater slope than the control treatment. However, all estimates for the effect of soil moisture on ANPP overlapped zero, indicating the relationship is weak and uncertain regardless of treatment. There was also large spatial variation in ANPP within-years, which contributes to the uncertainty of the soil moisture effect. Plant community composition was remarkably stable over the course of the experiment and did not differ among treatments.

Discussion

Despite some evidence that ecosystem functional response became more sensitive under sustained drought conditions, the response of ANPP to soil moisture was consistently weak and community composition was stable. The similarity of ecosystem functional responses across treatments was not related to compensatory shifts at the plant community level, but instead may reflect the insensitivity of the dominant species to soil moisture. These species may be successful precisely because they have evolved life history strategies that buffer them against precipitation variability.

Direct effects of warming increase woody plant abundance in a subarctic wetland

Both the direct effects of warming on a species' vital rates and indirect effects of warming caused by interactions with neighboring species can influence plant populations. Furthermore, herbivory mediates the effects of warming on plant community composition in many systems. Thus, determining the importance of direct and indirect effects of warming, while considering the role of herbivory, can help predict long-term plant community dynamics. We conducted a field experiment in the coastal wetlands of western Alaska to investigate how warming and herbivory influence the interactions and abundances of two common plant species, a sedge, Carex ramenskii, and a dwarf shrub, Salix ovalifolia. We used results from the experiment to model the equilibrium abundances of the species under different warming and grazing scenarios and to determine the contribution of direct and indirect effects to predict population changes. Consistent with the current composition of the landscape, model predictions suggest that Carex is more abundant than Salix under ambient temperatures with grazing (53% and 27% cover, respectively). However, with warming and grazing, Salix becomes more abundant than Carex (57% and 41% cover, respectively), reflecting both a negative response of Carex and a positive response of Salix to warming. While grazing reduced the cover of both species, herbivory did not prevent a shift in dominance from sedges to the dwarf shrub. Direct effects of climate change explained about 97% of the total predicted change in species cover, whereas indirect effects explained only 3% of the predicted change. Thus, indirect effects, mediated by interactions between Carex and Salix, were negligible, likely due to use of different niches and weak interspecific interactions. Results suggest that a 2°C increase could cause a shift in dominance from sedges to woody plants on the coast of western Alaska over decadal timescales, and this shift was largely a result of the direct effects of warming. Models predict this shift with or without goose herbivory. Our results are consistent with other studies showing an increase in woody plant abundance in the Arctic and suggest that shifts in plant-plant interactions are not driving this change.

Multi-model comparison highlights consistency in predicted effect of warming on a semi-arid shrub

A number of modeling approaches have been developed to predict the impacts of climate change on species distributions, performance, and abundance. The stronger the agreement from models that represent different processes and are based on distinct and independent sources of information, the greater the confidence we can have in their predictions. Evaluating the level of confidence is particularly important when predictions are used to guide conservation or restoration decisions. We used a multi-model approach to predict climate change impacts on big sagebrush (Artemisia tridentata), the dominant plant species on roughly 43 million hectares in the western United States and a key resource for many endemic wildlife species. To evaluate the climate sensitivity of A. tridentata, we developed four predictive models, two based on empirically derived spatial and temporal relationships, and two that applied mechanistic approaches to simulate sagebrush recruitment and growth. This approach enabled us to produce an aggregate index of climate change vulnerability and uncertainty based on the level of agreement between models. Despite large differences in model structure, predictions of sagebrush response to climate change were largely consistent. Performance, as measured by change in cover, growth, or recruitment, was predicted to decrease at the warmest sites, but increase throughout the cooler portions of sagebrush's range. A sensitivity analysis indicated that sagebrush performance responds more strongly to changes in temperature than precipitation. Most of the uncertainty in model predictions reflected variation among the ecological models, raising questions about the reliability of forecasts based on a single modeling approach. Our results highlight the value of a multi-model approach in forecasting climate change impacts and uncertainties and should help land managers to maximize the value of conservation investments.

The effects of intransitive competition on coexistence

Coexistence theory has been developed with an almost exclusive focus on interactions between two species, often ignoring more complex and indirect interactions, such as intransitive loops, that can emerge in competition networks. In fact, intransitive competition has typically been studied in isolation from other pairwise stabilising processes, and thus little is known about how intransitivity interacts with more traditional drivers of species coexistence such as niche partitioning. To integrate intransitivity into traditional coexistence theory, we developed a metric of growth rate when rare, Δri¯, to identify and quantify the impact of intransitive competition against a backdrop of pairwise stabilising niche differences. Using this index with simulations of community dynamics, we demonstrate that intransitive loops can both stabilise or destabilise species coexistence, but the strength and importance of intransitive interactions are significantly affected by the length and the topology of these loops. We conclude by showing how Δri¯ can be used to evaluate effects of intransitivity in empirical studies. Our results emphasise the need to integrate complex mechanisms emerging from diverse interactions into our understanding of species coexistence.