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.     

Global patterns in the shape of species geographical ranges reveal range determinants

Aim Do species range shapes follow general patterns? If so, what mechanisms underlie those patterns? We show for 11,582 species from a variety of taxa across the world that most species have similar latitudinal and longitudinal ranges. We then seek to disentangle the roles of climate, extrinsic dispersal limitation (e.g. barriers) and intrinsic dispersal limitation (reflecting a species’ ability to disperse) as constraints of species range shape. We also assess the relationship between range size and shape. Location Global. Methods Range shape patterns were measured as the slope of the regression of latitudinal species ranges against longitudinal ranges for each taxon and continent, and as the coefficient of determination measuring the degree of scattering of species ranges from the 1:1 line (i.e. latitudinal range = longitudinal range). Two major competing hypotheses explaining species distributions (i.e. dispersal or climatic determinism) were explored. To this end, we compared the observed slopes and coefficients of determination with those predicted by a climatic null model that estimates the potential range shapes in the absence of dispersal limitation. The predictions compared were that species distribution shapes are determined purely by (1) intrinsic dispersal limitation, (2) extrinsic dispersal limitations such as topographic barriers, and (3) climate. Results  Using this methodology, we show for a wide variety of taxa across the globe that species generally have very similar latitudinal and longitudinal ranges. However, neither neutral models assuming random but spatially constrained dispersal, nor models assuming climatic control of species distributions describe range shapes adequately. The empirical relationship between the latitudinal and longitudinal ranges of species falls between the predictions of these competing models. Main conclusions We propose that this pattern arises from the combined effect of macroclimate and intrinsic dispersal limitation, the latter being the major determinant among restricted‐range species. Hence, accurately projecting the impact of climate change onto species ranges will require a solid understanding of how climate and dispersal jointly control species ranges.     

Predicting species distributions across the Amazonian and Andean regions using remote sensing data

Aim We explore the utility of newly available optical and microwave remote sensing data from the Moderate Resolution Imaging Spectroradiometer (MODIS) and QuikSCAT (QSCAT) instruments for species distribution modelling at regional to continental scales. Using eight Neotropical species from three taxonomic groups, we assess the extent to which remote sensing data can improve predictions of their geographic distributions. For two bird species, we investigate the specific contributions of different types of remote sensing variables to the predictions and model accuracy at the regional scale, where the benefits of the MODIS and QSCAT satellite data are expected to be most significant. Location South America, with a focus on the tropical and subtropical Andes and the Amazon Basin. Methods Potential geographic distributions of eight species, namely two birds, two mammals and four trees, were modelled with the maxent algorithm at 1‐km resolution over the South American continent using climatic and remote sensing data separately and combined. For each species and model scenario, we assess model performance by testing the agreement between observed and simulated distributions across all thresholds and, in the case of the two focal bird species, at selected thresholds. Results Quantitative performance tests showed that models built with remote sensing and climatic layers in isolation performed well in predicting species distributions, suggesting that each of these data sets contains useful information. However, predictions created with a combination of remote sensing and climatic layers generally resulted in the best model performance across the three taxonomic groups. In Ecuador, the inclusion of remote sensing data was critical in resolving the known geographically isolated populations of the two focal bird species along the steep Amazonian–Andean elevational gradients. Within remote sensing subsets, microwave‐based data were more important than optical data in the predictions of the two bird species. Main conclusions Our results suggest that the newly available remote sensing data (MODIS and QSCAT) have considerable utility in modelling the contemporary geographical distributions of species at both regional and continental scales and in predicting range shifts as a result of large‐scale land‐use change.     

Predicting species distributions from small numbers of occurrence records: a test case using cryptic geckos in Madagascar

Aim Techniques that predict species potential distributions by combining observed occurrence records with environmental variables show much potential for application across a range of biogeographical analyses. Some of the most promising applications relate to species for which occurrence records are scarce, due to cryptic habits, locally restricted distributions or low sampling effort. However, the minimum sample sizes required to yield useful predictions remain difficult to determine. Here we developed and tested a novel jackknife validation approach to assess the ability to predict species occurrence when fewer than 25 occurrence records are available. Location Madagascar. Methods Models were developed and evaluated for 13 species of secretive leaf‐tailed geckos (Uroplatus spp.) that are endemic to Madagascar, for which available sample sizes range from 4 to 23 occurrence localities (at 1 km² grid resolution). Predictions were based on 20 environmental data layers and were generated using two modelling approaches: a method based on the principle of maximum entropy (Maxent) and a genetic algorithm (GARP). Results We found high success rates and statistical significance in jackknife tests with sample sizes as low as five when the Maxent model was applied. Results for GARP at very low sample sizes (less than c. 10) were less good. When sample sizes were experimentally reduced for those species with the most records, variability among predictions using different combinations of localities demonstrated that models were greatly influenced by exactly which observations were included. Main conclusions We emphasize that models developed using this approach with small sample sizes should be interpreted as identifying regions that have similar environmental conditions to where the species is known to occur, and not as predicting actual limits to the range of a species. The jackknife validation approach proposed here enables assessment of the predictive ability of models built using very small sample sizes, although use of this test with larger sample sizes may lead to overoptimistic estimates of predictive power. Our analyses demonstrate that geographical predictions developed from small numbers of occurrence records may be of great value, for example in targeting field surveys to accelerate the discovery of unknown populations and species.     

Projected changes in distributions of Australian tropical savanna birds under climate change using three dispersal scenarios

Identifying the species most vulnerable to extinction as a result of climate change is a necessary first step in mitigating biodiversity decline. Species distribution modeling (SDM) is a commonly used tool to assess potential climate change impacts on distributions of species. We use SDMs to predict geographic ranges for 243 birds of Australian tropical savannas, and to project changes in species richness and ranges under a future climate scenario between 1990 and 2080. Realistic predictions require recognition of the variability in species capacity to track climatically suitable environments. Here we assess the effect of dispersal on model results by using three approaches: full dispersal, no dispersal and a partial-dispersal scenario permitting species to track climate change at a rate of 30 km per decade. As expected, the projected distributions and richness patterns are highly sensitive to the dispersal scenario. Projected future range sizes decreased for 66% of species if full dispersal was assumed, but for 89% of species when no dispersal was assumed. However, realistic future predictions should not assume a single dispersal scenario for all species and as such, we assigned each species to the most appropriate dispersal category based on individual mobility and habitat specificity; this permitted the best estimates of where species will be in the future. Under this "realistic" dispersal scenario, projected ranges sizes decreased for 67% of species but showed that migratory and tropical-endemic birds are predicted to benefit from climate change with increasing distributional area. Richness hotspots of tropical savanna birds are expected to move, increasing in southern savannas and southward along the east coast of Australia, but decreasing in the arid zone. Understanding the complexity of effects of climate change on species' range sizes by incorporating dispersal capacities is a crucial step toward developing adaptation policies for the conservation of vulnerable species.     

Incorporating biotic interactions in the distribution models of African wild silk moths (Gonometa species,Lasiocampidae) using different representations of modelled host tree distributions

Biotic interactions influence species niches and may thus shape distributions. Nevertheless, species distribution modelling has traditionally relied exclusively on environmental factors to predict species distributions, while biotic interactions have only seldom been incorporated into models. This study tested the ability of incorporating biotic interactions, in the form of host plant distributions, to increase model performance for two host‐dependent lepidopterans of economic interest, namely the African silk moth species, Gonometa postica and Gonometa rufobrunnea (Lasiocampidae). Both species are dependent on a small number of host tree species for the completion of their life cycle. We thus expected the host plant distribution to be an important predictor of Gonometa distributions. Model performance of a species distribution model trained only on abiotic predictors was compared to four species distribution models that additionally incorporated biotic interactions in the form of four different representations of host plant distributions as predictors. We found that incorporating the moth–host plant interactions improved G. rufobrunnea model performance for all representations of host plant distribution, while for G. postica model performance only improved for one representation of host plant distribution. The best performing representation of host plant distribution differed for the two Gonometa species. While these results suggest that incorporating biotic interactions into species distribution models can improve model performance, there is inconsistency in which representation of the host tree distribution best improves predictions. Therefore, the ability of biotic interactions to improve species distribution models may be context‐specific, even for species which have obligatory interactions with other organisms.     

Disentangling scale dependencies in species environmental niches and distributions

Understanding species’ responses to environmental conditions, and how these ­species–environment associations shape spatial distributions, are longstanding goals in ecology and biogeography. However, an essential component of species–environment relationships – the spatial unit, or grain, at which they operate – remains unresolved. We identify three components of scale‐dependence in analyses of species–environment associations: 1) response grain, the grain at which species respond most strongly to their environment; 2) environment spatial structure, the pattern of spatial autocorrelation intrinsic to an environmental factor; and 3) analysis grain, the grain at which analyses are conducted and ecological inferences are made. We introduce a novel conceptual framework that defines these scale components in the context of analyzing species–environment relationships, and provide theoretical examples of their interactions for species with various ecological attributes. We then use a virtual species approach to investigate the impacts of each component on common methods of measuring and predicting species–environment relationships. We find that environment spatial structure has a substantial impact on the ability of even simple, univariate species distribution models (SDMs) to recover known species–­environment associations at coarse analysis grains. For simulated environments with ‘fine’ and ‘intermediate’ spatial structure, model explanatory power, and the frequency with which simple SDMs correctly estimated a virtual species’ response to the simulated environment, dramatically declined as analysis grain increased. Informed by these results, we use a scaling analysis to identify maximum analysis grains for individual environmental factors, and a scale optimization procedure to determine the grain of maximum predictive accuracy. Implementing these analysis grain thresholds and model performance standards in an example east African study system yields more accurate distribution predictions, compared to SDMs independently constructed at arbitrary analysis grains. Finally, we integrate our conceptual framework with virtual and empirical results to provide practical recommendations for researchers asking common questions about species–environment relationships.     

Spatial autocorrelation in predictors reduces the impact of positional uncertainty in occurrence data on species distribution modelling

Aim To investigate the impact of positional uncertainty in species occurrences on the predictions of seven commonly used species distribution models (SDMs), and explore its interaction with spatial autocorrelation in predictors. Methods A series of artificial datasets covering 155 scenarios including different combinations of five positional uncertainty scenarios and 31 spatial autocorrelation scenarios were simulated. The level of positional uncertainty was defined by the standard deviation of a normally distributed zero‐mean random variable. Each dataset included two environmental gradients (predictor variables) and one set of species occurrence sample points (response variable). Seven commonly used models were selected to develop SDMs: generalized linear models, generalized additive models, boosted regression trees, multivariate adaptive regression spline, random forests, genetic algorithm for rule‐set production and maximum entropy. A probabilistic approach was employed to model and simulate five levels of error in the species locations. To analyse the propagation of positional uncertainty, Monte Carlo simulation was applied to each scenario for each SDM. The models were evaluated for performance using simulated independent test data with Cohen’s Kappa and the area under the receiver operating characteristic curve. Results Positional uncertainty in species location led to a reduction in prediction accuracy for all SDMs, although the magnitude of the reduction varied between SDMs. In all cases the magnitude of this impact varied according to the degree of spatial autocorrelation in predictors and the levels of positional uncertainty. It was shown that when the range of spatial autocorrelation in the predictors was less than or equal to three times the standard deviation of the positional error, the models were less affected by error and, consequently, had smaller decreases in prediction accuracy. When the range of spatial autocorrelation in predictors was larger than three times the standard deviation of positional error, the prediction accuracy was low for all scenarios. Main conclusions The potential impact of positional uncertainty in species occurrences on the predictions of SDMs can be understood by comparing it with the spatial autocorrelation range in predictor variables.     

Modelling invasive alien species distributions from digital biodiversity atlases. Model upscaling as a means of reconciling data at different scales

Aim

There is a wealth of information on species occurrences in biodiversity data banks, albeit presence‐only, biased and scarce at fine resolutions. Moreover, fine‐resolution species maps are required in biodiversity conservation. New techniques for dealing with this kind of data have been reported to perform well. These fine‐resolution maps would be more robust if they could explain data at coarser resolutions at which species distributions are well represented. We present a new methodology for testing this hypothesis and apply it to invasive alien species (IAS).

Location

Catalonia, Spain.

Methods

We used species presence records from the Biodiversity data bank of Catalonia to model the distribution of ten IAS which, according to some recent studies, achieve their maximum distribution in the study area. To overcome problems inherent with the data, we prepared different correction treatments: three for dealing with bias and five for autocorrelation. We used the MaxEnt algorithm to generate models at 1‐km resolution for each species and treatment. Acceptable models were upscaled to 10 km and validated against independent 10 km occurrence data.

Results

Of a total of 150 models, 20 gave acceptable results at 1‐km resolution and 12 passed the cross‐scale validation test. No apparent pattern emerged, which could serve as a guide on modelling. Only four species gave models that also explained the distribution at the coarser scale.

Main conclusions

Although some techniques may apparently deliver good distribution maps for species with scarce and biased data, they need to be taken with caution. When good independent data at a coarser scale are available, cross‐scale validation can help to produce more reliable and robust maps. When no independent data are available for validation, however, new data gathering field surveys may be the only option if reliable fine‐scale resolution maps are needed.     

Impact of landscape predictors on climate change modelling of species distributions: a case study with Eucalyptus fastigata in southern New South Wales,Australia

Aim We consider three questions. (1) How different are the predicted distribution maps when climate‐only and climate‐plus‐terrain models are developed from high‐resolution data? (2) What are the implications of differences between the models when predicting future distributions under climate change scenarios, particularly for climate‐only models at coarse resolution? (3) Does the use of high‐resolution data and climate‐plus‐terrain models predict an increase in the number of local refugia? Location South‐eastern New South Wales, Australia. Methods We developed two species distribution models for Eucalyptus fastigata under current climate conditions using generalized additive modelling. One used only climate variables as predictors (mean annual temperature, mean annual rainfall, mean summer rainfall); the other used both climate and landscape (June daily radiation, topographic position, lithology, nutrients) variables as predictors. Predictions of the distribution under current climate and climate change were then made for both models at a pixel resolution of 100 m. Results The model using climate and landscape variables as predictors explained a significantly greater proportion of the deviance than the climate‐only model. Inclusion of landscape variables resulted in the prediction of much larger areas of existing optimal habitat. An overlay of predicted future climate on the current climate space indicated that extrapolation of the statistical models was not occurring and models were therefore more robust. Under climate change, landscape‐defined refugia persisted in areas where the climate‐only model predicted major declines. In areas where expansion was predicted, the increase in optimal habitat was always greater with landscape predictors. Recognition of extensive optimal habitat conditions and potential refugia was dependent on the use of high‐resolution landscape data. Main conclusions Using only climate variables as predictors for assessing species responses to climate change ignores the accepted conceptual model of plant species distribution. Explicit statements justifying the selection of predictors based on ecological principles are needed. Models using only climate variables overestimate range reduction under climate change and fail to predict potential refugia. Fine‐scale‐resolution data are required to capture important climate/landscape interactions. Extrapolation of statistical models to regions in climate space outside the region where they were fitted is risky.     

Uncertainty in ensemble forecasting of species distribution

Species distribution modelling has been widely applied in order to assess the potential impacts of climate change on biodiversity. Many methodological decisions, taken during the modelling process and forecasts, may, however, lead to a large variability in the assessment of future impacts. Using measures of species range change and turnover, the potential impacts of climate change on French stream fish species and assemblages were evaluated. Our main focus was to quantify the uncertainty in the projections of these impacts arising from four sources of uncertainty: initial datasets (Data), statistical methods [species distribution models (SDM)], general circulation models (GCM), and gas emission scenarios (GES). Several modalities of the aforementioned uncertainty sources were combined in an ensemble forecasting framework resulting in 8400 different projections. The variance explained by each source was then extracted from this whole ensemble of projections. Overall, SDM contributed to the largest variation in projections, followed by GCM, whose contribution increased over time equalling almost the proportion of variance explained by SDM in 2080. Data and GES had little influence on the variability in projections. Future projections of range change were more consistent for species with a large geographical extent (i.e., distribution along latitudinal or stream gradients) or with restricted environmental requirements (i.e., small thermal or elevation ranges). Variability in projections of turnover was spatially structured at the scale of France, indicating that certain particular geographical areas should be considered with care when projecting the potential impacts of climate change. The results of this study, therefore, emphasized that particular attention should be paid to the use of predictions ensembles resulting from the application of several statistical methods and climate models. Moreover, forecasted impacts of climate change should always be provided with an assessment of their uncertainty, so that management and conservation decisions can be taken in the full knowledge of their reliability.     

Influences of ecology and biogeography on shaping the distributions of cryptic species: three bat tales in Iberia

To determine what shapes the distributions of cryptic species, we aimed to unravel ecological niches and geographical distributions of three cryptic bat species complexes in Iberia, Plecotus auritus/begognae, Myotis mystacinus/alcathoe and Eptesicus serotinus/isabellinus (with 44, 69, 66, 27, 121 and 216 records, respectively), considering ecological interactions and biogeographical patterns. Species distribution models (SDMs) were built using a presence‐only technique (Maxent), incorporating genetically identified species records with environmental variables (climate, habitat, topography). The most relevant variables for each species’ distribution and respective response curves were then determined. SDMs for each species were overlapped to assess the contact zones within each complex. Niche analyses were performed using niche metrics and spatial principal component analyses to study niche overlap and breadth. The Plecotus complex showed a parapatric distribution, although having similar biogeographical affinities (Eurosiberian), possibly explained by competitive exclusion. The Myotis complex also showed Eurosiberian affinities, with high overlap between niches and distribution, suggesting resource partitioning between species. Finally, E. serotinus was associated with Eurosiberian areas, while E. isabellinus occurred in Mediterranean areas, suggesting possible competition in their restricted contact zone. This study highlights the relevance of considering potential ecological interactions between similarly ecological species when assessing species distributions. © 2014 The Linnean Society of London, Biological Journal of the Linnean Society, 2014, 112 ,150–162.     

Deconstructing the mammal species richness pattern in Europe – towards an understanding of the relative importance of climate,biogeographic history,habitat heterogeneity and humans

Aim We deconstructed the mammal species richness pattern in Europe to assess the importance of large‐scale gradients in current macroclimate relative to biogeographic history, habitat heterogeneity and human influence (HHH variables) as richness determinants for total species, and for widespread and endemic species separately. Location Europe, west of 30° E. Methods We deconstructed total species richness (50‐km resolution) into its widespread and endemic species richness components. We used simultaneous autoregressive modelling (SAR) with information‐theoretic model selection and variation partitioning to assess the importance of macroclimate and HHH variables. The HHH variables included two historical factors, estimated by novel methodologies: (1) ice‐age‐driven dynamics, represented by accessibility to recolonization from hindcasting‐estimated glacial refugia, and (2) biogeographic peninsular dynamics, represented by distance to the entry region for the main European faunal source in western Asia. Results A large fraction of explained variation was shared between macroclimate and HHH in the SAR models. For total species richness, more variation could be uniquely attributed to macroclimate than to HHH, whereas for the deconstructed patterns (widespread and endemic species) the opposite was the case. Considering the individual factors, there was a strong peninsula effect on both widespread and endemic species richness but not on total richness. Main conclusions Both macroclimate and HHH variables (history, habitat heterogeneity and human influence) proved important predictors of species richness, but also difficult to disentangle. Notably, biogeographic history, in particular peninsular dynamics, is an important determinant of widespread and endemic species richness.