Accounting for dispersal and biotic interactions to disentangle the drivers of species distributions and their abundances

Although abiotic factors, together with dispersal and biotic interactions, are often suggested to explain the distribution of species and their abundances, species distribution models usually focus on abiotic factors only. We propose an integrative framework linking ecological theory, empirical data and statistical models to understand the distribution of species and their abundances together with the underlying community assembly dynamics. We illustrate our approach with 21 plant species in the French Alps. We show that a spatially nested modelling framework significantly improves the model's performance and that the spatial variations of species presence-absence and abundances are predominantly explained by different factors. We also show that incorporating abiotic, dispersal and biotic factors into the same model bring new insights to our understanding of community assembly. This approach, at the crossroads between community ecology and biogeography, is a promising avenue for a better understanding of species co-existence and biodiversity distribution.     

Context‐dependent biotic interactions control plant abundance across altitudinal environmental gradients

Many biotic interactions influence community structure, yet most distribution models for plants have focused on plant competition or used only abiotic variables to predict plant abundance. Furthermore, biotic interactions are commonly context‐dependent across abiotic gradients. For example, plant–plant interactions can grade from competition to facilitation over temperature gradients. We used a hierarchical Bayesian framework to predict the abundances of 12 plant species across a mountain landscape and test hypotheses on the context‐dependency of biotic interactions over abiotic gradients. We combined field‐based estimates of six biotic interactions (foliar herbivory and pathogen damage, fungal root colonization, fossorial mammal disturbance, plant cover and plant diversity) with abiotic data on climate and soil depth, nutrients and moisture. All biotic interactions were significantly context‐dependent along temperature gradients. Results supported the stress gradient hypothesis: as abiotic stress increased, the strength or direction of the relationship between biotic variables and plant abundance generally switched from negative (suggesting suppressed plant abundance) to positive (suggesting facilitation/mutualism). For half of the species, plant cover was the best predictor of abundance, suggesting that the prior focus on plant–plant interactions is well‐justified. Explicitly incorporating the context‐dependency of biotic interactions generated novel hypotheses about drivers of plant abundance across abiotic gradients and may improve the accuracy of niche models.     

Revisiting the Holy Grail: using plant functional traits to understand ecological processes

One of ecology's grand challenges is developing general rules to explain and predict highly complex systems. Understanding and predicting ecological processes from species' traits has been considered a ‘Holy Grail’ in ecology. Plant functional traits are increasingly being used to develop mechanistic models that can predict how ecological communities will respond to abiotic and biotic perturbations and how species will affect ecosystem function and services in a rapidly changing world; however, significant challenges remain. In this review, we highlight recent work and outstanding questions in three areas: (i) selecting relevant traits; (ii) describing intraspecific trait variation and incorporating this variation into models; and (iii) scaling trait data to community‐ and ecosystem‐level processes. Over the past decade, there have been significant advances in the characterization of plant strategies based on traits and trait relationships, and the integration of traits into multivariate indices and models of community and ecosystem function. However, the utility of trait‐based approaches in ecology will benefit from efforts that demonstrate how these traits and indices influence organismal, community, and ecosystem processes across vegetation types, which may be achieved through meta‐analysis and enhancement of trait databases. Additionally, intraspecific trait variation and species interactions need to be incorporated into predictive models using tools such as Bayesian hierarchical modelling. Finally, existing models linking traits to community and ecosystem processes need to be empirically tested for their applicability to be realized.     

Incorporating biotic factors in species distribution modeling: are interactions with soil microbes important?

It is increasingly recognized that species distributions are driven by both abiotic factors and biotic interactions. Despite much recent work incorporating competition, predation, and mutualism into species distribution models (SDMs), the focus has been confined to aboveground macroscopic interactions. Biotic interactions between plants and soil microbial communities are understudied as potentially important drivers of plant distributions. Some soil bacteria promote plant growth by cycling nutrients, while others are pathogenic; thus they have a high potential for influencing plant occurrence. We investigated the influence of soil bacterial clades on the distributions of bryophytes and 12 vascular plant species in a high elevation talus‐field ecosystem in the Rocky Mountain Front Range, Colorado, USA. We used an information‐theoretic criterion (AICc) modeling approach to compare SDMs with the following different sets of predictors: abiotic variables, abiotic variables and other plant abundances, abiotic variables and soil bacteria clade relative abundances, and a full model with abiotic factors, plant abundances, and bacteria relative abundances. We predicted that bacteria would influence plant distributions both positively and negatively, and that these interactions would improve prediction of plant species distributions. We found that inclusion of either plant or bacteria biotic predictors generally improved the fit, deviance explained, and predictive power of the SDMs, and for the majority of the species, adding information on both other plants and bacteria yielded the best model. Interactions between the modeled species and biotic predictors were both positive and negative, suggesting the presence of competition, parasitism, and facilitation. While our results indicate that plant–plant co‐occurrences are a stronger driver of plant distributions than plant–bacteria co‐occurrences, they also show that bacteria can explain parts of plant distributions that remain unexplained by abiotic and plant predictors. Our results provide further support for including biotic factors in SDMs, and suggest that belowground factors be considered as well.     

Enhanced effects of biotic interactions on predicting multispecies spatial distribution of submerged macrophytes after eutrophication

Water eutrophication creates unfavorable environmental conditions for submerged macrophytes. In these situations, biotic interactions may be particularly important for explaining and predicting the submerged macrophytes occurrence. Here, we evaluate the roles of biotic interactions in predicting spatial occurrence of submerged macrophytes in 1959 and 2009 for Dianshan Lake in eastern China, which became eutrophic since the 1980s. For the four common species occurred in 1959 and 2009, null species distribution models based on abiotic variables and full models based on both abiotic and biotic variables were developed using generalized linear model (GLM) and boosted regression trees (BRT) to determine whether the biotic variables improved the model performance. Hierarchical Bayesian‐based joint species distribution models capable of detecting paired biotic interactions were established for each species in both periods to evaluate the changes in the biotic interactions. In most of the GLM and BRT models, the full models showed better performance than the null models in predicting the species presence/absence, and the relative importance of the biotic variables in the full models increased from less than 50% in 1959 to more than 50% in 2009 for each species. Moreover, co‐occurrence correlation of each paired species interaction was higher in 2009 than that in 1959. The findings suggest biotic interactions that tend to be positive play more important roles in the spatial distribution of multispecies assemblages of macrophytes and should be included in prediction models to improve prediction accuracy when forecasting macrophytes’ distribution under eutrophication stress.     

Linking macroecology and community ecology: refining predictions of species distributions using biotic interaction networks

Macroecological models for predicting species distributions usually only include abiotic environmental conditions as explanatory variables, despite knowledge from community ecology that all species are linked to other species through biotic interactions. This disconnect is largely due to the different spatial scales considered by the two sub‐disciplines: macroecologists study patterns at large extents and coarse resolutions, while community ecologists focus on small extents and fine resolutions. A general framework for including biotic interactions in macroecological models would help bridge this divide, as it would allow for rigorous testing of the role that biotic interactions play in determining species ranges. Here, we present an approach that combines species distribution models with Bayesian networks, which enables the direct and indirect effects of biotic interactions to be modelled as propagating conditional dependencies among species’ presences. We show that including biotic interactions in distribution models for species from a California grassland community results in better range predictions across the western USA. This new approach will be important for improving estimates of species distributions and their dynamics under environmental change.     

Can functional groups be used to indicate estuarine ecological status?

International legislation demands that statutory bodies report on the health of aquatic ecosystems. Traditionally, ecosystem components have been characterised according to species assemblages but with limited success in predicting health. On the other hand, many studies based upon functional groupings that include trophic relationships and bioturbation potential have shown response to pollution. However, these and other functional group responses have not yet been linked to broad scale physical variables. To date this has hindered the development of a predictive model of function based on abiotic factors. In addition, most functional studies ignore any potential role of body size when assessing the importance of each species to overall functional group measures. By weighting all species that belong to the same guild equally, the investigator risks overestimating the true importance of any one guild to the environment. This study compared the ability of different functional group approaches to discriminate between separate estuarine sites, whilst linking biotic data with abiotic factors. Using data for the Tamar Estuary, we show that no two methods of classifying the biotic data, according to function, produce the same groupings of sites; nor did any method produce groupings that matched clusters based on abiotic factors alone. Instead, results show that not only can choice of functional method alter our perception of site associations but also, can influence the strength of similarity relationships between abiotic and biotic datasets. Both the use of bioturbation measures and weighting species abundance data by body size provided better relationships between biotic and abiotic data than the use of trophic groups. Thus both methods merit further research to produce algorithms for modelling studies.     

Biotic interactions boost spatial models of species richness

Biotic interactions are known to affect the composition of species assemblages via several mechanisms, such as competition and facilitation. However, most spatial models of species richness do not explicitly consider inter‐specific interactions. Here, we test whether incorporating biotic interactions into high‐resolution models alters predictions of species richness as hypothesised. We included key biotic variables (cover of three dominant arctic‐alpine plant species) into two methodologically divergent species richness modelling frameworks – stacked species distribution models (SSDM) and macroecological models (MEM) – for three ecologically and evolutionary distinct taxonomic groups (vascular plants, bryophytes and lichens). Predictions from models including biotic interactions were compared to the predictions of models based on climatic and abiotic data only. Including plant–plant interactions consistently and significantly lowered bias in species richness predictions and increased predictive power for independent evaluation data when compared to the conventional climatic and abiotic data based models. Improvements in predictions were constant irrespective of the modelling framework or taxonomic group used. The global biodiversity crisis necessitates accurate predictions of how changes in biotic and abiotic conditions will potentially affect species richness patterns. Here, we demonstrate that models of the spatial distribution of species richness can be improved by incorporating biotic interactions, and thus that these key predictor factors must be accounted for in biodiversity forecasts.     

Habitat shifts of endangered species under altered climate conditions: importance of biotic interactions

Predicting changes in potential habitat for endangered species as a result of global warming requires considering more than future climate conditions; it is also necessary to evaluate biotic associations. Most distribution models predicting species responses to climate change include climate variables and occasionally topographic and edaphic parameters, rarely are biotic interactions included. Here, we incorporate biotic interactions into niche models to predict suitable habitat for species under altered climates. We constructed and evaluated niche models for an endangered butterfly and a threatened bird species, both are habitat specialists restricted to semiarid shrublands of southern California. To incorporate their dependency on shrubs, we first developed climate‐based niche models for shrubland vegetation and individual shrub species. We also developed models for the butterfly's larval host plants. Outputs from these models were included in the environmental variable dataset used to create butterfly and bird niche models. For both animal species, abiotic–biotic models outperformed the climate‐only model, with climate‐only models over‐predicting suitable habitat under current climate conditions. We used the climate‐only and abiotic–biotic models to calculate amounts of suitable habitat under altered climates and to evaluate species' sensitivities to climate change. We varied temperature (+0.6, +1.7, and +2.8 °C) and precipitation (50%, 90%, 100%, 110%, and 150%) relative to current climate averages and within ranges predicted by global climate change models. Suitable habitat for each species was reduced at all levels of temperature increase. Both species were sensitive to precipitation changes, particularly increases. Under altered climates, including biotic variables reduced habitat by 68–100% relative to the climate‐only model. To design reserve systems conserving sensitive species under global warming, it is important to consider biotic interactions, particularly for habitat specialists and species with strong dependencies on other species.     

Assessing congruence of opportunistic records and systematic surveys for predicting Hispaniolan mammal species distributions

Comparative assessment of the relative information content of different independent spatial data types is necessary to evaluate whether they provide congruent biogeographic signals for predicting species ranges. Opportunistic occurrence records and systematically collected survey data are available from the Dominican Republic for Hispaniola’s surviving endemic non‐volant mammals, the Hispaniolan solenodon (Solenodon paradoxus) and Hispaniolan hutia (Plagiodontia aedium); opportunistic records (archaeological, historical and recent) exist from across the entire country, and systematic survey data have been collected from seven protected areas. Species distribution models were developed in maxent for solenodons and hutias using both data types, with species habitat suitability and potential country‐level distribution predicted using seven biotic and abiotic environmental variables. Three different models were produced and compared for each species: (a) opportunistic model, with starting model incorporating abiotic‐only predictors; (b) total survey model, with starting model incorporating biotic and abiotic predictors; and (c) reduced survey model, with starting model incorporating abiotic‐only predictors to allow further comparison with the opportunistic model. All models predict suitable environmental conditions for both solenodons and hutias across a broadly congruent, relatively large area of the Dominican Republic, providing a spatial baseline of conservation‐priority landscapes that might support native mammals. Correlation between total and reduced survey models is high for both species, indicating the substantial explanatory power of abiotic variables for predicting Hispaniolan mammal distributions. However, correlation between survey models and opportunistic models is only moderately positive. Species distribution models derived from different data types can provide different predictions about habitat suitability and conservation‐priority landscapes for threatened species, likely reflecting incompleteness and bias in spatial sampling associated with both data types. Models derived using both opportunistic and systematic data must therefore be applied critically and cautiously.     

How traits determine species responses to environmental gradients

The role traits play in species distributions is central to developing a mechanistic and predictive plant ecology. Developments in generalized linear mixed models described by Jamil et al. in this issue provide a significant step in this endeavour. Trait‐based GLMMs are intuitive parametric models for exploring and predicting the distribution of species along environmental gradients, as affected by species traits.     

Mechanistic forecasts of species responses to climate change: The promise of biophysical ecology

A core challenge in global change biology is to predict how species will respond to future environmental change and to manage these responses. To make such predictions and management actions robust to novel futures, we need to accurately characterize how organisms experience their environments and the biological mechanisms by which they respond. All organisms are thermodynamically connected to their environments through the exchange of heat and water at fine spatial and temporal scales and this exchange can be captured with biophysical models. Although mechanistic models based on biophysical ecology have a long history of development and application, their use in global change biology remains limited despite their enormous promise and increasingly accessible software. We contend that greater understanding and training in the theory and methods of biophysical ecology is vital to expand their application. Our review shows how biophysical models can be implemented to understand and predict climate change impacts on species' behavior, phenology, survival, distribution, and abundance. It also illustrates the types of outputs that can be generated, and the data inputs required for different implementations. Examples range from simple calculations of body temperature at a particular site and time, to more complex analyses of species' distribution limits based on projected energy and water balances, accounting for behavior and phenology. We outline challenges that currently limit the widespread application of biophysical models relating to data availability, training, and the lack of common software ecosystems. We also discuss progress and future developments that could allow these models to be applied to many species across large spatial extents and timeframes. Finally, we highlight how biophysical models are uniquely suited to solve global change biology problems that involve predicting and interpreting responses to environmental variability and extremes, multiple or shifting constraints, and novel abiotic or biotic environments.     

The Importance of Biotic vs. Abiotic Drivers of Local Plant Community Composition Along Regional Bioclimatic Gradients

We assessed if the relative importance of biotic and abiotic factors for plant community composition differs along environmental gradients and between functional groups, and asked which implications this may have in a warmer and wetter future. The study location is a unique grid of sites spanning regional-scale temperature and precipitation gradients in boreal and alpine grasslands in southern Norway. Within each site we sampled vegetation and associated biotic and abiotic factors, and combined broad- and fine-scale ordination analyses to assess the relative explanatory power of these factors for species composition. Although the community responses to biotic and abiotic factors did not consistently change as predicted along the bioclimatic gradients, abiotic variables tended to explain a larger proportion of the variation in species composition towards colder sites, whereas biotic variables explained more towards warmer sites, supporting the stress gradient hypothesis. Significant interactions with precipitation suggest that biotic variables explained more towards wetter climates in the sub alpine and boreal sites, but more towards drier climates in the colder alpine. Thus, we predict that biotic interactions may become more important in alpine and boreal grasslands in a warmer future, although more winter precipitation may counteract this trend in oceanic alpine climates. Our results show that both local and regional scales analyses are needed to disentangle the local vegetation-environment relationships and their regional-scale drivers, and biotic interactions and precipitation must be included when predicting future species assemblages.     

Exploring macroinvertebrate species distributions at regional and local scales across a sandy beach geographic continuum

Exposed sandy beaches are highly dynamic ecosystems where macroinvertebrate species cope with extremely variable environmental conditions. The majority of the beach ecology studies present exposed beaches as physically dominated ecosystems where abiotic factors largely determine the structure and distribution of macrobenthic communities. However, beach species patterns at different scales can be modified by the interaction between different environmental variables, including biotic interactions. In this study, we examined the role of different environmental variables for describing the regional and local scale distributions of common macrobenthic species across 39 beaches along the North coast of Spain. The analyses were carried out using boosted regression trees, a relatively new technique from the field of machine learning. Our study showed that the macroinvertebrate community on exposed beaches is not structured by a single physical factor, but instead by a complex set of drivers including the biotic compound. Thus, at a regional scale the macrobenthic community, in terms of number of species and abundance, was mainly explained by surrogates of food availability, such as chlorophyll a. The results also revealed that the local scale is a feasible way to construct general predictive species-environmental models, since relationships derived from different beaches showed similar responses for most of the species. However, additional information on aspects of beach species distribution can be obtained with large scale models. This study showed that species-environmental models should be validated against changes in spatial extent, and also illustrates the utility of BRTs as a powerful analysis tool for ecology data insight.     

Predicting future coexistence in a North American ant community

Global climate change will remodel ecological communities worldwide. However, as a consequence of biotic interactions, communities may respond to climate change in idiosyncratic ways. This makes predictive models that incorporate biotic interactions necessary. We show how such models can be constructed based on empirical studies in combination with predictions or assumptions regarding the abiotic consequences of climate change. Specifically, we consider a well‐studied ant community in North America. First, we use historical data to parameterize a basic model for species coexistence. Using this model, we determine the importance of various factors, including thermal niches, food discovery rates, and food removal rates, to historical species coexistence. We then extend the model to predict how the community will restructure in response to several climate‐related changes, such as increased temperature, shifts in species phenology, and altered resource availability. Interestingly, our mechanistic model suggests that increased temperature and shifts in species phenology can have contrasting effects. Nevertheless, for almost all scenarios considered, we find that the most subordinate ant species suffers most as a result of climate change. More generally, our analysis shows that community composition can respond to climate warming in nonintuitive ways. For example, in the context of a community, it is not necessarily the most heat‐sensitive species that are most at risk. Our results demonstrate how models that account for niche partitioning and interspecific trade‐offs among species can be used to predict the likely idiosyncratic responses of local communities to climate change.     

A General and Dynamic Species Abundance Model, Embracing the Lognormal and the Gamma Models

One aspect of community ecology that has been given particular attention is the pattern of species abundances in a community. The species may have a wide range of abundances; some are very common and others rare. When species abundance models are fitted to observations, the lognormal model and one of the gamma models (e.g., the log-series model) are usually applied. The model that gives the best fit according to some goodness-of-fit test is then chosen. By applying a diffusion approximation for each species' dynamics with density regulation of the straight theta-logistic type, we here present a general species abundance model that embraces the two most widely applied species abundance models, the lognormal and the gamma. Our general model will, therefore, provide a better fit than the two special cases, except when it corresponds to one of them. In contrast to the classical models, ours is also dynamic, making it possible to evaluate the fluctuations in species abundance over time through both biotic and abiotic factors. The model is fitted to several species abundance data sets and our results compared to previous attempts to fit a model, usually either the lognormal or the log-series.     

Plant Phenology Supports the Multi-emergence Hypothesis for Ebola Spillover Events

Ebola virus disease outbreaks in animals (including humans and great apes) start with sporadic host switches from unknown reservoir species. The factors leading to such spillover events are little explored. Filoviridae viruses have a wide range of natural hosts and are unstable once outside hosts. Spillover events, which involve the physical transfer of viral particles across species, could therefore be directly promoted by conditions of host ecology and environment. In this report, we outline a proof of concept that temporal fluctuations of a set of ecological and environmental variables describing the dynamics of the host ecosystem are able to predict such events of Ebola virus spillover to humans and animals. We compiled a data set of climate and plant phenology variables and Ebola virus disease spillovers in humans and animals. We identified critical biotic and abiotic conditions for spillovers via multiple regression and neural network-based time series regression. Phenology variables proved to be overall better predictors than climate variables. African phenology variables are not yet available as a comprehensive online resource. Given the likely importance of phenology for forecasting the likelihood of future Ebola spillover events, our results highlight the need for cost-effective transect surveys to supply phenology data for predictive modelling efforts.     

种间作用对物种分布模拟的影响及建模方法

物种分布模型(SDMs)通过量化物种分布和环境变量之间的关系,并将其外推到未知的景观单元,模拟、预测地理空间中生物的潜在分布,是生态学、生物地理学、保护生物学等研究领域的重要工具.然而,目前物种分布模型主要采用非生物因素作为预测变量,由于数据量化和建模表达困难,生物因素特别是种间作用在物种分布模型中常被忽略,将种间作用...     

No species is an island: testing the effects of biotic interactions on models of avian niche occupation

Traditionally, the niche of a species is described as a hypothetical 3D space, constituted by well‐known biotic interactions (e.g. predation, competition, trophic relationships, resource–consumer interactions, etc.) and various abiotic environmental factors. Species distribution models (SDMs), also called “niche models” and often used to predict wildlife distribution at landscape scale, are typically constructed using abiotic factors with biotic interactions generally been ignored. Here, we compared the goodness of fit of SDMs for red‐backed shrike Lanius collurio in farmlands of Western Poland, using both the classical approach (modeled only on environmental variables) and the approach which included also other potentially associated bird species. The potential associations among species were derived from the relevant ecological literature and by a correlation matrix of occurrences. Our findings highlight the importance of including heterospecific interactions in improving our understanding of niche occupation for bird species. We suggest that suite of measures currently used to quantify realized species niches could be improved by also considering the occurrence of certain associated species. Then, an hypothetical “species 1” can use the occurrence of a successfully established individual of “species 2” as indicator or “trace” of the location of available suitable habitat to breed. We hypothesize this kind of biotic interaction as the “heterospecific trace effect” (HTE): an interaction based on the availability and use of “public information” provided by individuals from different species. Finally, we discuss about the incomes of biotic interactions for enhancing the predictive capacities on species distribution models.