Predicting species distributions based on incomplete survey data: the trade‐off between precision and scale

Systematic species surveys over large areas are mostly not affordable, constraining conservation planners to make best use of incomplete data. Spatially explicit species distribution models (SDM) may be useful to detect and compensate for incomplete information. SDMs can either be based on standardized, systematic sampling in a restricted subarea, or – as a cost‐effective alternative – on data haphazardly collated by “volunteer‐based monitoring schemes” (VMS), area‐wide but inherently biased and of heterogeneous spatial precision. Using data on capercaillie Tetrao urogallus, we evaluated the capacity of SDMs generated from incomplete survey data to localise unknown areas inhabited by the species and to predict relative local observation density. Addressing the trade‐off between data precision, sample size and spatial extent of the sampling area, we compared three different sampling strategies: VMS‐data collected throughout the whole study area (7000 km²) using either 1) exact locations or 2) locations aggregated to grid cells of the size of an average individual home range, and 3) systematic transect counts conducted within a small subarea (23.8 km²). For each strategy, we compared two sample sizes and two modelling methods (ENFA and Maxent), which were evaluated using cross‐validation and independent data. Models based on VMS‐data (strategies 1 and 2) performed equally well in predicting relative observation density and in localizing “unknown” occurrences. They always outperformed strategy 3‐models, irrespective of sample size and modelling method, partly because the VMS‐data provided the more comprehensive clues for setting the discrimination‐threshold for predicting presence or absence. Accounting for potential errors due to extrapolation (e.g. projections outside the environmental domain or potentially biasing variables) reduced, but did not fully compensate for the observed discrepancies. As they cover a broader range of species‐habitat relations, the area‐wide data achieved a better model quality with less a‐priori knowledge. Furthermore, in a highly mobile species like capercaillie a sampling resolution corresponding to an individuals' home range can lead to equally good predictions as the use of exact locations. Consequently, when a trade‐off between the sampling effort and the spatial extent of the sampling area is necessary, less precise data unsystematically collected over a large representative region are preferable to systematically sampled data from a restricted region.     

Intra‐specific variability and plasticity influence potential tree species distributions under climate change

Aim To assess the effect of local adaptation and phenotypic plasticity on the potential distribution of species under future climate changes. Trees may be adapted to specific climatic conditions; however, species range predictions have classically been assessed by species distribution models (SDMs) that do not account for intra‐specific genetic variability and phenotypic plasticity, because SDMs rely on the assumption that species respond homogeneously to climate change across their range, i.e. a species is equally adapted throughout its range, and all species are equally plastic. These assumptions could cause SDMs to exaggerate or underestimate species at risk under future climate change. Location The Iberian Peninsula. Methods Species distributions are predicted by integrating experimental data and modelling techniques. We incorporate plasticity and local adaptation into a SDM by calibrating models of tree survivorship with adaptive traits in provenance trials. Phenotypic plasticity was incorporated by calibrating our model with a climatic index that provides a measure of the differences between sites and provenances. Results We present a new modelling approach that is easy to implement and makes use of existing tree provenance trials to predict species distribution models under global warming. Our results indicate that the incorporation of intra‐population genetic diversity and phenotypic plasticity in SDMs significantly altered their outcome. In comparing species range predictions, the decrease in area occupancy under global warming conditions is smaller when considering our survival–adaptation model than that predicted by a ‘classical SDM’ calibrated with presence–absence data. These differences in survivorship are due to both local adaptation and plasticity. Differences due to the use of experimental data in the model calibration are also expressed in our results: we incorporate a null model that uses survival data from all provenances together. This model always predicts less reduction in area occupancy for both species than the SDM calibrated with presence–absence. Main conclusions We reaffirm the importance of considering adaptive traits when predicting species distributions and avoiding the use of occurrence data as a predictive variable. In light of these recommendations, we advise that existing predictions of future species distributions and their component populations must be reconsidered.     

Effects of ecohydrological variables on current and future ranges,local suitability patterns,and model accuracy in big sagebrush

Forecasting of species and ecosystem responses to novel conditions, including climate change, is one of the major challenges facing ecologists at the start of the 21st century. Climate change studies based on species distribution models (SDMs) have been criticized because they extend correlational relationships beyond the observed data. Here, we compared conventional climate‐based SDMs against ecohydrological SDMs that include information from process‐based simulations of water balance. We examined the current and future distribution of Artemisia tridentata (big sagebrush) representing sagebrush ecosystems, which are widespread in semiarid western North America. For each approach, we calculated ensemble models from nine SDM methods and tested accuracy of each SDM with a null distribution. Climatic conditions included current conditions for 1970–1999 and two IPCC projections B1 and A2 for 2070–2099. Ecohydrological conditions were assessed by simulating soil water balance with SOILWAT, a daily time‐step, multiple layer, mechanistic, soil water model. Under current conditions, both climatic and ecohydrological SDM approaches produced comparable sagebrush distributions. Overall, sagebrush distribution is forecasted to decrease, with larger decreases under the A2 than under the B1 scenario and strong decreases in the southern part of the range. Increases were forecasted in the northern parts and at higher elevations. Both SDM approaches produced accurate predictions. However, the ecohydrological SDM approach was slightly less accurate than climatic SDMs (?1% in AUC, ?4% in Kappa and TSS) and predicted a higher number of habitat patches than observed in the input data. Future predictions of ecohydrological SDMs included an increased number of habitat patches whereas climatic SDMs predicted a decrease. This difference is important for understanding landscape‐scale patterns of sagebrush ecosystems and management of sagebrush obligate species for future conditions. Several mechanisms can explain the diverging forecasts; however, we need better insights into the consequences of different datasets for SDMs and how these affect our understanding of future trajectories.     

Minimum required number of specimen records to develop accurate species distribution models

Species distribution models (SDMs) are widely used to predict the occurrence of species. Because SDMs generally use presence‐only data, validation of the predicted distribution and assessing model accuracy is challenging. Model performance depends on both sample size and species’ prevalence, being the fraction of the study area occupied by the species. Here, we present a novel method using simulated species to identify the minimum number of records required to generate accurate SDMs for taxa of different pre‐defined prevalence classes. We quantified model performance as a function of sample size and prevalence and found model performance to increase with increasing sample size under constant prevalence, and to decrease with increasing prevalence under constant sample size. The area under the curve (AUC) is commonly used as a measure of model performance. However, when applied to presence‐only data it is prevalence‐dependent and hence not an accurate performance index. Testing the AUC of an SDM for significant deviation from random performance provides a good alternative. We assessed the minimum number of records required to obtain good model performance for species of different prevalence classes in a virtual study area and in a real African study area. The lower limit depends on the species’ prevalence with absolute minimum sample sizes as low as 3 for narrow‐ranged and 13 for widespread species for our virtual study area which represents an ideal, balanced, orthogonal world. The lower limit of 3, however, is flawed by statistical artefacts related to modelling species with a prevalence below 0.1. In our African study area lower limits are higher, ranging from 14 for narrow‐ranged to 25 for widespread species. We advocate identifying the minimum sample size for any species distribution modelling by applying the novel method presented here, which is applicable to any taxonomic clade or group, study area or climate scenario.     

Species distribution models as a tool to estimate reproductive parameters: a case study with a passerine bird species

1.?Correlative species distribution models (SDMs) assess relationships between species distribution data and environmental features, to evaluate the environmental suitability (ES) of a given area for a species, by providing a measure of the probability of presence. If the output of SDMs represents the relationships between habitat features and species performance well, SDM results can be related also to other key parameters of populations, including reproductive parameters. To test this hypothesis, we evaluated whether SDM results can be used as a proxy of reproductive parameters (breeding output, territory size) in red-backed shrikes (Lanius collurio). 2.?The distribution of 726 shrike territories in Northern Italy was obtained through multiple focused surveys; for a subset of pairs, we also measured territory area and number of fledged juveniles. We used Maximum Entropy modelling to build a SDM on the basis of territory distribution. We used generalized least squares and spatial generalized mixed models to relate territory size and number of fledged juveniles to SDM suitability, while controlling for spatial autocorrelation. 3.?Species distribution models predicted shrike distribution very well. Territory size was negatively related to suitability estimated through SDM, while the number of fledglings significantly increased with the suitability of the territory. This was true also when SDM was built using only spatially and temporally independent data. 4.?Results show a clear relationship between ES estimated through presence-only SDMs and two key parameters related to species' reproduction, suggesting that suitability estimated by SDM, and habitat quality determining reproduction parameters in our model system, are correlated. Our study shows the potential use of SDMs to infer important fitness parameters; this information can have great importance in management and conservation.     

Not seeing the forest for the trees: Generalised linear model out-performs random forest in species distribution modelling for Southeast Asian felids

Species Distribution Models (SDMs) are a powerful tool to derive habitat suitability predictions relating species occurrence data with habitat features. Two of the most frequently applied algorithms to model species-habitat relationships are Generalised Linear Models (GLM) and Random Forest (RF). The former is a parametric regression model providing functional models with direct interpretability. The latter is a machine learning non-parametric algorithm, more tolerant than other approaches in its assumptions, which has often been shown to outperform parametric algorithms. Other approaches have been developed to produce robust SDMs, like training data bootstrapping and spatial scale optimisation. Using felid presence-absence data from three study regions in Southeast Asia (mainland, Borneo and Sumatra), we tested the performances of SDMs by implementing four modelling frameworks: GLM and RF with bootstrapped and non-bootstrapped training data. With Mantel and ANOVA tests we explored how the four combinations of algorithms and bootstrapping influenced SDMs and their predictive performances. Additionally, we tested how scale-optimisation responded to species' size, taxonomic associations (species and genus), study area and algorithm. We found that choice of algorithm had strong effect in determining the differences between SDMs' spatial predictions, while bootstrapping had no effect. Additionally, algorithm followed by study area and species, were the main factors driving differences in the spatial scales identified. SDMs trained with GLM showed higher predictive performance, however, ANOVA tests revealed that algorithm had significant effect only in explaining the variance observed in sensitivity and specificity and, when interacting with bootstrapping, in Percent Correctly Classified (PCC). Bootstrapping significantly explained the variance in specificity, PCC and True Skills Statistics (TSS). Our results suggest that there are systematic differences in the scales identified and in the predictions produced by GLM vs. RF, but that neither approach was consistently better than the other. The divergent predictions and inconsistent predictive abilities suggest that analysts should not assume machine learning is inherently superior and should test multiple methods. Our results have strong implications for SDM development, revealing the inconsistencies introduced by the choice of algorithm on scale optimisation, with GLM selecting broader scales than RF.     

Sensitivity of predictive species distribution models to change in grain size

Predictive species distribution modelling (SDM) has become an essential tool in biodiversity conservation and management. The choice of grain size (resolution) of environmental layers used in modelling is one important factor that may affect predictions. We applied 10 distinct modelling techniques to presence-only data for 50 species in five different regions, to test whether: (1) a 10-fold coarsening of resolution affects predictive performance of SDMs, and (2) any observed effects are dependent on the type of region, modelling technique, or species considered. Results show that a 10 times change in grain size does not severely affect predictions from species distribution models. The overall trend is towards degradation of model performance, but improvement can also be observed. Changing grain size does not equally affect models across regions, techniques, and species types. The strongest effect is on regions and species types, with tree species in the data sets (regions) with highest locational accuracy being most affected. Changing grain size had little influence on the ranking of techniques: boosted regression trees remain best at both resolutions. The number of occurrences used for model training had an important effect, with larger sample sizes resulting in better models, which tended to be more sensitive to grain. Effect of grain change was only noticeable for models reaching sufficient performance and/or with initial data that have an intrinsic error smaller than the coarser grain size.     

Why georeferencing matters: Introducing a practical protocol to prepare species occurrence records for spatial analysis

Species Distribution Models (SDMs) are widely used to understand environmental controls on species’ ranges and to forecast species range shifts in response to climatic changes. The quality of input data is crucial determinant of the model's accuracy. While museum records can be useful sources of presence data for many species, they do not always include accurate geographic coordinates. Therefore, actual locations must be verified through the process of georeferencing. We present a practical, standardized manual georeferencing method (the Spatial Analysis Georeferencing Accuracy (SAGA) protocol) to classify the spatial resolution of museum records specifically for building improved SDMs. We used the high‐elevation plant Saxifraga austromontana Wiegand (Saxifragaceae) as a case study to test the effect of using this protocol when developing an SDM. In MAXENT, we generated and compared SDMs using a comprehensive occurrence dataset that had undergone three different levels of georeferencing: (1) trained using all publicly available herbarium records of the species, minus outliers (2) trained using herbarium records claimed to be previously georeferenced, and (3) trained using herbarium records that we have manually georeferenced to a ≤ 1‐km resolution using the SAGA protocol. Model predictions of suitable habitat for S. austromontana differed greatly depending on georeferencing level. The SDMs fitted with presence locations georeferenced using SAGA outperformed all others. Differences among models were exacerbated for future distribution predictions. Under rapid climate change, accurately forecasting the response of species becomes increasingly important. Failure to georeference location data and cull inaccurate samples leads to erroneous model output, limiting the utility of spatial analyses. We present a simple, standardized georeferencing method to be adopted by curators, ecologists, and modelers to improve the geographic accuracy of museum records and SDM predictions.     

Niche variability and its consequences for species distribution modeling

When species distribution models (SDMs) are used to predict how a species will respond to environmental change, an important assumption is that the environmental niche of the species is conserved over evolutionary time-scales. Empirical studies conducted at ecological time-scales, however, demonstrate that the niche of some species can vary in response to environmental change. We use habitat and locality data of five species of stream fishes collected across seasons to examine the effects of niche variability on the accuracy of projections from Maxent, a popular SDM. We then compare these predictions to those from an alternate method of creating SDM projections in which a transformation of the environmental data to similar scales is applied. The niche of each species varied to some degree in response to seasonal variation in environmental variables, with most species shifting habitat use in response to changes in canopy cover or flow rate. SDMs constructed from the original environmental data accurately predicted the occurrences of one species across all seasons and a subset of seasons for two other species. A similar result was found for SDMs constructed from the transformed environmental data. However, the transformed SDMs produced better models in ten of the 14 total SDMs, as judged by ratios of mean probability values at known presences to mean probability values at all other locations. Niche variability should be an important consideration when using SDMs to predict future distributions of species because of its prevalence among natural populations. The framework we present here may potentially improve these predictions by accounting for such variability.     

Incorporating physiology into species distribution models moderates the projected impact of warming on selected Mediterranean marine species

Species distribution models (SDMs) correlate species occurrences with environmental predictors, and can be used to forecast distributions under future climates. SDMs have been criticized for not explicitly including the physiological processes underlying the species response to the environment. Recently, new methods have been suggested to combine SDMs with physiological estimates of performance (physiology-SDMs). In this study, we compare SDM and physiology-SDM predictions for select marine species in the Mediterranean Sea, a region subjected to exceptionally rapid climate change. We focused on six species and created physiology-SDMs that incorporate physiological thermal performance curves from experimental data with species occurrence records. We then contrasted projections of SDMs and physiology-SDMs under future climate (year 2100) for the entire Mediterranean Sea, and particularly the ‘warm’ trailing edge in the Levant region. Across the Mediterranean, we found cross-validation model performance to be similar for regular SDMs and physiology-SDMs. However, we also show that for around half the species the physiology-SDMs substantially outperform regular SDM in the warm Levant. Moreover, for all species the uncertainty associated with the coefficients estimated from the physiology-SDMs were much lower than in the regular SDMs. Under future climate, we find that both SDMs and physiology-SDMs showed similar patterns, with species predicted to shift their distribution north-west in accordance with warming sea temperatures. However, for the physiology-SDMs predicted distributional changes are more moderate than those predicted by regular SDMs. We conclude, that while physiology-SDM predictions generally agree with the regular SDMs, incorporation of the physiological data led to less extreme range shift forecasts. The results suggest that climate-induced range shifts may be less drastic than previously predicted, and thus most species are unlikely to completely disappear with warming climate. Taken together, the findings emphasize that physiological experimental data can provide valuable supplemental information to predict range shifts of marine species.     

Assessing multitemporal calibration for species distribution models

Species distribution models (SDMs), the most prominent tool in modern biogeography, rely on the assumptions that (i) species distribution is in equilibrium with the environment and (ii) that climatic niche has been conserved throughout recent geological time. These issues affect the spatial and temporal transferability of SDMs, limiting their reliability for applications such as when studying effects of past climate change on species distribution and extinctions. The integration of paleontological and neontological data for a multitemporal calibration and validation of SDMs has been suggested for improving SDMs flexibility. Here, we provide an empirical test for a multitemporal calibration, employing virtual species (i.e., with perfectly-known distributions) and comparing them directly with monotemporal SDMs (i.e., SDM calibrated in a single time layer). We used 1kyr-interval scenarios throughout the last 22 kyr BP for two ecologically different species in South America (a “hot and wet” species and a “cold and dry” species). Models with multitemporal calibration performed similarly to models with monotemporal calibration, regardless of species, sample sizes, and time frame. However, multitemporal calibration performed better when dealing with non-analogous climates among time layers. By improving the temporal SDMs transferability, multitemporal calibration opens new avenues for integrating fossil and recent occurrence data, which may substantially benefit biogeography and paleoecology.     

Keep collecting: accurate species distribution modelling requires more collections than previously thought

Aim Species distribution models (SDMs) use the locations of collection records to map the distributions of species, making them a powerful tool in conservation biology, ecology and biogeography. However, the accuracy of range predictions may be reduced by temporally autocorrelated biases in the data. We assess the accuracy of SDMs in predicting the ranges of tropical plant species on the basis of different sample sizes while incorporating real‐world collection patterns and biases. Location Tropical South American moist forests. Methods We use dated herbarium records to model the distributions of 65 Amazonian and Andean plant species. For each species, we use the first 25, 50, 100, 125 and 150 records collected and available for each species to analyse changes in spatial aggregation and climatic representativeness through time. We compare the accuracy of SDM range estimates produced using the time‐ordered data subsets to the accuracy of range estimates generated using the same number of collections but randomly subsampled from all available records. Results We find that collections become increasingly aggregated through time but that additional collecting sites are added resulting in progressively better representations of the species’ full climatic niches. The range predictions produced using time‐ordered data subsets are less accurate than predictions from random subsets of equal sample sizes. Range predictions produced using time‐ordered data subsets consistently underestimate the extent of ranges while no such tendency exists for range predictions produced using random data subsets. Main conclusions These results suggest that larger sample sizes are required to accurately map species ranges. Additional attention should be given to increasing the number of records available per species through continued collecting, better distributed collecting, and/or increasing access to existing collections. The fact that SDMs generally under‐predict the extent of species ranges means that extinction risks of species because of future habitat loss may be lower than previously estimated.     

Convolutional neural networks improve species distribution modelling by capturing the spatial structure of the environment

Convolutional Neural Networks (CNNs) are statistical models suited for learning complex visual patterns. In the context of Species Distribution Models (SDM) and in line with predictions of landscape ecology and island biogeography, CNN could grasp how local landscape structure affects prediction of species occurrence in SDMs. The prediction can thus reflect the signatures of entangled ecological processes. Although previous machine-learning based SDMs can learn complex influences of environmental predictors, they cannot acknowledge the influence of environmental structure in local landscapes (hence denoted “punctual models”). In this study, we applied CNNs to a large dataset of plant occurrences in France (GBIF), on a large taxonomical scale, to predict ranked relative probability of species (by joint learning) to any geographical position. We examined the way local environmental landscapes improve prediction by performing alternative CNN models deprived of information on landscape heterogeneity and structure (“ablation experiments”). We found that the landscape structure around location crucially contributed to improve predictive performance of CNN-SDMs. CNN models can classify the predicted distributions of many species, as other joint modelling approaches, but they further prove efficient in identifying the influence of local environmental landscapes. CNN can then represent signatures of spatially structured environmental drivers. The prediction gain is noticeable for rare species, which open promising perspectives for biodiversity monitoring and conservation strategies. Therefore, the approach is of both theoretical and practical interest. We discuss the way to test hypotheses on the patterns learnt by CNN, which should be essential for further interpretation of the ecological processes at play.     

Incorporating intraspecific variation into species distribution models improves distribution predictions,but cannot predict species traits for a wide-spread plant species

The most common approach to predicting how species ranges and ecological functions will shift with climate change is to construct correlative species distribution models (SDMs). These models use a species’ climatic distribution to determine currently suitable areas for the species and project its potential distribution under future climate scenarios. A core, rarely tested, assumption of SDMs is that all populations will respond equivalently to climate. Few studies have examined this assumption, and those that have rarely dissect the reasons for intraspecific differences. Focusing on the arctic-alpine cushion plant Silene acaulis, we compared predictive accuracy from SDMs constructed using the species’ full global distribution with composite predictions from separate SDMs constructed using subpopulations defined either by genetic or habitat differences. This is one of the first studies to compare multiple ways of constructing intraspecific-level SDMs with a species-level SDM. We also examine the contested relationship between relative probability of occurrence and species performance or ecological function, testing if SDM output can predict individual performance (plant size) and biotic interactions (facilitation). We found that both genetic- and habitat-informed SDMs are considerably more accurate than a species-level SDM, and that the genetic model substantially differs from and outperforms the habitat model. While SDMs have been used to infer population performance and possibly even biotic interactions, in our system these relationships were extremely weak. Our results indicate that individual subpopulations may respond differently to climate, although we discuss and explore several alternative explanations for the superior performance of intraspecific-level SDMs. We emphasize the need to carefully examine how to best define intraspecific-level SDMs as well as how potential genetic, environmental, or sampling variation within species ranges can critically affect SDM predictions. We urge caution in inferring population performance or biotic interactions from SDM predictions, as these often-assumed relationships are not supported in our study.     

Where is positional uncertainty a problem for species distribution modelling?

Species data held in museum and herbaria, survey data and opportunistically observed data are a substantial information resource. A key challenge in using these data is the uncertainty about where an observation is located. This is important when the data are used for species distribution modelling (SDM), because the coordinates are used to extract the environmental variables and thus, positional error may lead to inaccurate estimation of the species–environment relationship. The magnitude of this effect is related to the level of spatial autocorrelation in the environmental variables. Using local spatial association can be relevant because it can lead to the identification of the specific occurrence records that cause the largest drop in SDM accuracy. Therefore, in this study, we tested whether the SDM predictions are more affected by positional uncertainty originating from locations that have lower local spatial association in their predictors. We performed this experiment for Spain and the Netherlands, using simulated datasets derived from well known species distribution models (SDMs). We used the K statistic to quantify the local spatial association in the predictors at each species occurrence location. A probabilistic approach using Monte Carlo simulations was employed to introduce the error in the species locations. The results revealed that positional uncertainty in species occurrence data at locations with low local spatial association in predictors reduced the prediction accuracy of the SDMs. We propose that local spatial association is a way to identify the species occurrence records that require treatment for positional uncertainty. We also developed and present a tool in the R environment to target observations that are likely to create error in the output from SDMs as a result of positional uncertainty.     

Transferable species distribution modelling: Comparative performance of Generalised Functional Response models

Predictive species distribution models (SDMs) are becoming increasingly important in ecology, in the light of rapid environmental change. However, the predictions of most current SDMs are specific to the habitat composition of the environments in which they were fitted. This may limit SDM predictive power because species may respond differently to a given habitat depending on the availability of all habitats in their environment, a phenomenon known as a functional response in resource selection. The Generalised Functional Response (GFR) framework captures this dependence by formulating the SDM coefficients as functions of habitat availability. The original GFR implementation used global polynomial functions of habitat availability to describe the functional responses. In this study, we develop several refinements of this approach and compare their predictive performance using two simulated and two real datasets. We first use local radial basis functions (RBF), a more flexible approach than global polynomials, to represent the habitat selection coefficients, and balance bias with precision via regularization to prevent overfitting. Second, we use the RBF-GFR and GFR models in combination with the classification and regression tree CART, which has more flexibility and better predictive powers for non-linear modelling. As further extensions, we use random forests (RFs) and extreme gradient boosting (XGBoost), ensemble approaches that consistently lead to variance reduction in generalization error. We find that the different methods are ranked consistently across the datasets for out-of-data prediction. The traditional stationary approach to SDMs and the GFR model consistently perform at the bottom of the ranking (simple SDMs underfit, and polynomial GFRs overfit the data). The best methods in our list provide non-negligible improvements in predictive performance, in some cases taking the out-of-sample R2 from 0.3 up to 0.7 across datasets. At times of rapid environmental change and spatial non-stationarity ignoring the effects of functional responses on SDMs, results in two different types of prediction bias (under-prediction or mis-positioning of distribution hotspots). However, not all functional response models perform equally well. The more volatile polynomial GFR models can generate biases through over-prediction. Our results indicate that there are consistently robust GFR approaches that achieve impressive gains in transferability across very different datasets.     

Terrestrial or marine species distribution model: Why not both? A case study with seabirds

Species reliant on both the terrestrial and marine realms present a challenge for conventional species distribution models (SDMs). For such species, standard single‐realm SDMs may omit key information that could result in decreased model accuracy and performance. Existing approaches to habitat suitability modeling typically do not effectively combine information from multiple realms; this methodological gap can ultimately hamper management efforts for groups such as seabirds, seals, and turtles. This study, for the first time, jointly incorporates both terrestrial information and marine information into a single species distribution model framework. We do this by sampling nearby marine conditions for a given terrestrial point and vice versa using parameters set by each species’ mean maximum foraging distance and then use standard SDM methods to generate habitat suitability predictions; therefore, our method does not rely on post hoc combination of several different models. Using three seabird species with very different ecologies, we investigate whether this new multi‐realm approach can improve our ability to identify suitable habitats for these species. Results show that incorporating terrestrial information into marine SDMs, or vice versa, generally improves model performance, sometimes drastically. However, there is considerable variability between species in the level of improvement as well as in the particular method that produces the most improvement. Our approach provides a repeatable and transparent method to combine information from multiple ecological realms in a single SDM framework. Important advantages over existing solutions include the opportunity to, firstly, easily combine terrestrial and marine information for species that forage large distances inland or out to sea and, secondly, consider interactions between terrestrial and marine variables.     

Climatic associations of British species distributions show good transferability in time but low predictive accuracy for range change

Conservation planners often wish to predict how species distributions will change in response to environmental changes. Species distribution models (SDMs) are the primary tool for making such predictions. Many methods are widely used; however, they all make simplifying assumptions, and predictions can therefore be subject to high uncertainty. With global change well underway, field records of observed range shifts are increasingly being used for testing SDM transferability. We used an unprecedented distribution dataset documenting recent range changes of British vascular plants, birds, and butterflies to test whether correlative SDMs based on climate change provide useful approximations of potential distribution shifts. We modelled past species distributions from climate using nine single techniques and a consensus approach, and projected the geographical extent of these models to a more recent time period based on climate change; we then compared model predictions with recent observed distributions in order to estimate the temporal transferability and prediction accuracy of our models. We also evaluated the relative effect of methodological and taxonomic variation on the performance of SDMs. Models showed good transferability in time when assessed using widespread metrics of accuracy. However, models had low accuracy to predict where occupancy status changed between time periods, especially for declining species. Model performance varied greatly among species within major taxa, but there was also considerable variation among modelling frameworks. Past climatic associations of British species distributions retain a high explanatory power when transferred to recent time--due to their accuracy to predict large areas retained by species--but fail to capture relevant predictors of change. We strongly emphasize the need for caution when using SDMs to predict shifts in species distributions: high explanatory power on temporally-independent records--as assessed using widespread metrics--need not indicate a model's ability to predict the future.