Climate versus weather extremes: Temporal predictor resolution matters for future rather than current regional species distribution models

Aim

Climate is considered a major driver of species distributions. Long‐term climatic means are commonly used as predictors in correlative species distribution models (SDMs). However, this coarse temporal resolution does not reflect local conditions that populations experience, such as short‐term weather extremes, which may have a strong impact on population dynamics and local distributions. We here compare the performance of climate‐ and weather‐based predictors in regional SDMs and their influence on future predictions, which are increasingly used in conservation planning.

Location

South‐western Germany.

Methods

We built different SDMs for 20 Orthoptera species based on three predictor sets at a regional scale for current and future climate scenarios. We calculated standard bioclimatic variables and yearly and seasonal sets of climate change indicating variables of weather extremes. As the impact of extreme events may be stronger for habitat specialists than for generalists, we distinguished species’ degrees of specialization. We computed linear mixed‐effects models to identify significant effects of algorithm, predictor set and specialization on model performance and calculated correlations and geographical niche overlap between spatial predictions.

Results

Current predictions were rather similar among all predictor sets, but highly variable for future climate scenarios. Bioclimatic and seasonal weather predictors performed slightly better than yearly weather predictors, though performance differences were minor. We found no evidence that specialists are more sensitive to weather extremes than generalists.

Main conclusions

For future projections of species distributions, SDM predictor selection should not solely be based on current performances and predictions. As long‐term climate and short‐term weather predictors represent different environmental drivers of a species’ distribution, we argue to interpret diverging future projections as complements. Even if similar current performances and predictions might imply their equivalency, favouring one predictor set neglects important aspects of future distributions and might mislead conservation decisions based on them.

     

Incorporating spatial autocorrelation into species distribution models alters forecasts of climate‐mediated range shifts

Species distribution models (SDMs) are widely used to forecast changes in the spatial distributions of species and communities in response to climate change. However, spatial autocorrelation (SA) is rarely accounted for in these models, despite its ubiquity in broad‐scale ecological data. While spatial autocorrelation in model residuals is known to result in biased parameter estimates and the inflation of type I errors, the influence of unmodeled SA on species' range forecasts is poorly understood. Here we quantify how accounting for SA in SDMs influences the magnitude of range shift forecasts produced by SDMs for multiple climate change scenarios. SDMs were fitted to simulated data with a known autocorrelation structure, and to field observations of three mangrove communities from northern Australia displaying strong spatial autocorrelation. Three modeling approaches were implemented: environment‐only models (most frequently applied in species' range forecasts), and two approaches that incorporate SA; autologistic models and residuals autocovariate (RAC) models. Differences in forecasts among modeling approaches and climate scenarios were quantified. While all model predictions at the current time closely matched that of the actual current distribution of the mangrove communities, under the climate change scenarios environment‐only models forecast substantially greater range shifts than models incorporating SA. Furthermore, the magnitude of these differences intensified with increasing increments of climate change across the scenarios. When models do not account for SA, forecasts of species' range shifts indicate more extreme impacts of climate change, compared to models that explicitly account for SA. Therefore, where biological or population processes induce substantial autocorrelation in the distribution of organisms, and this is not modeled, model predictions will be inaccurate. These results have global importance for conservation efforts as inaccurate forecasts lead to ineffective prioritization of conservation activities and potentially to avoidable species extinctions.     

The fate of European breeding birds under climate,land‐use and dispersal scenarios

Many species have already shifted their distributions in response to recent climate change. Here, we aimed at predicting the future breeding distributions of European birds under climate, land‐use, and dispersal scenarios. We predicted current and future distributions of 409 species within an ensemble forecast framework using seven species distribution models (SDMs), five climate scenarios and three emission and land‐use scenarios. We then compared results from SDMs using climate‐only variables, habitat‐only variables or both climate and habitat variables. In order to account for a species’ dispersal abilities, we used natal dispersal estimates and developed a probabilistic method that produced a dispersal scenario intermediate between the null and full dispersal scenarios generally considered in such studies. We then compared results from all scenarios in terms of future predicted range changes, range shifts, and variations in species richness. Modeling accuracy was better with climate‐only variables than with habitat‐only variables, and better with both climate and habitat variables. Habitat models predicted smaller range shifts and smaller variations in range size and species richness than climate models. Using both climate and habitat variables, it was predicted that the range of 71% of the species would decrease by 2050, with a 335 km median shift. Predicted variations in species richness showed large decreases in the southern regions of Europe, as well as increases, mainly in Scandinavia and northern Russia. The partial dispersal scenario was significantly different from the full dispersal scenario for 25% of the species, resulting in the local reduction of the future predicted species richness of up to 10%. We concluded that the breeding range of most European birds will decrease in spite of dispersal abilities close to a full dispersal hypothesis, and that given the contrasted predictions obtained when modeling climate change only and land‐use change only, both scenarios must be taken into consideration.     

Benchmarking novel approaches for modelling species range dynamics

Increasing biodiversity loss due to climate change is one of the most vital challenges of the 21st century. To anticipate and mitigate biodiversity loss, models are needed that reliably project species’ range dynamics and extinction risks. Recently, several new approaches to model range dynamics have been developed to supplement correlative species distribution models (SDMs), but applications clearly lag behind model development. Indeed, no comparative analysis has been performed to evaluate their performance. Here, we build on process‐based, simulated data for benchmarking five range (dynamic) models of varying complexity including classical SDMs, SDMs coupled with simple dispersal or more complex population dynamic models (SDM hybrids), and a hierarchical Bayesian process‐based dynamic range model (DRM). We specifically test the effects of demographic and community processes on model predictive performance. Under current climate, DRMs performed best, although only marginally. Under climate change, predictive performance varied considerably, with no clear winners. Yet, all range dynamic models improved predictions under climate change substantially compared to purely correlative SDMs, and the population dynamic models also predicted reasonable extinction risks for most scenarios. When benchmarking data were simulated with more complex demographic and community processes, simple SDM hybrids including only dispersal often proved most reliable. Finally, we found that structural decisions during model building can have great impact on model accuracy, but prior system knowledge on important processes can reduce these uncertainties considerably. Our results reassure the clear merit in using dynamic approaches for modelling species’ response to climate change but also emphasize several needs for further model and data improvement. We propose and discuss perspectives for improving range projections through combination of multiple models and for making these approaches operational for large numbers of species.     

Biotic interactions influence the projected distribution of a specialist mammal under climate change

Aim

To measure the effects of including biotic interactions on climate‐based species distribution models (SDMs) used to predict distribution shifts under climate change. We evaluated the performance of distribution models for an endangered marsupial, the northern bettong (Bettongia tropica), comparing models that used only climate variables with models that also took into account biotic interactions.

Location

North‐east Queensland, Australia.

Methods

We developed separate climate‐based distribution models for the northern bettong, its two main resources and a competitor species. We then constructed models for the northern bettong by including climate suitability estimates for the resources and competitor as additional predictor variables to make climate + resource and climate + resource + competition models. We projected these models onto seven future climate scenarios and compared predictions of northern bettong distribution made by these differently structured models, using a ‘global’ metric, the I similarity statistic, to measure overlap in distribution and a ‘local’ metric to identify where predictions differed significantly.

Results

Inclusion of food resource biotic interactions improved model performance. Over moderate climate changes, up to 3.0 °C of warming, the climate‐only model for the northern bettong gave similar predictions of distribution to the more complex models including interactions, with differences only at the margins of predicted distributions. For climate changes beyond 3.0 °C, model predictions diverged significantly. The interactive model predicted less contraction of distribution than the simpler climate‐only model.

Main conclusions

Distribution models that account for interactions with other species, in particular direct resources, improve model predictions in the present‐day climate. For larger climate changes, shifts in distribution of interacting species cause predictions of interactive models to diverge from climate‐only models. Incorporating interactions with other species in SDMs may be needed for long‐term prediction of changes in distribution of species under climate change, particularly for specialized species strongly dependent on a small number of biotic interactions.     

Habitat availability and gene flow influence diverging local population trajectories under scenarios of climate change: a place‐based approach

Ecological niche theory holds that species distributions are shaped by a large and complex suite of interacting factors. Species distribution models (SDMs) are increasingly used to describe species’ niches and predict the effects of future environmental change, including climate change. Currently, SDMs often fail to capture the complexity of species’ niches, resulting in predictions that are generally limited to climate‐occupancy interactions. Here, we explore the potential impact of climate change on the American pika using a replicated place‐based approach that incorporates climate, gene flow, habitat configuration, and microhabitat complexity into SDMs. Using contemporary presence–absence data from occupancy surveys, genetic data to infer connectivity between habitat patches, and 21 environmental niche variables, we built separate SDMs for pika populations inhabiting eight US National Park Service units representing the habitat and climatic breadth of the species across the western United States. We then predicted occurrence probability under current (1981–2010) and three future time periods (out to 2100). Occurrence probabilities and the relative importance of predictor variables varied widely among study areas, revealing important local‐scale differences in the realized niche of the American pika. This variation resulted in diverse and – in some cases – highly divergent future potential occupancy patterns for pikas, ranging from complete extirpation in some study areas to stable occupancy patterns in others. Habitat composition and connectivity, which are rarely incorporated in SDM projections, were influential in predicting pika occupancy in all study areas and frequently outranked climate variables. Our findings illustrate the importance of a place‐based approach to species distribution modeling that includes fine‐scale factors when assessing current and future climate impacts on species’ distributions, especially when predictions are intended to manage and conserve species of concern within individual protected areas.     

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.     

Modeling plant species distributions under future climates: how fine scale do climate projections need to be?

Recent studies suggest that species distribution models (SDMs) based on fine‐scale climate data may provide markedly different estimates of climate‐change impacts than coarse‐scale models. However, these studies disagree in their conclusions of how scale influences projected species distributions. In rugged terrain, coarse‐scale climate grids may not capture topographically controlled climate variation at the scale that constitutes microhabitat or refugia for some species. Although finer scale data are therefore considered to better reflect climatic conditions experienced by species, there have been few formal analyses of how modeled distributions differ with scale. We modeled distributions for 52 plant species endemic to the California Floristic Province of different life forms and range sizes under recent and future climate across a 2000‐fold range of spatial scales (0.008–16 km²). We produced unique current and future climate datasets by separately downscaling 4 km climate models to three finer resolutions based on 800, 270, and 90 m digital elevation models and deriving bioclimatic predictors from them. As climate‐data resolution became coarser, SDMs predicted larger habitat area with diminishing spatial congruence between fine‐ and coarse‐scale predictions. These trends were most pronounced at the coarsest resolutions and depended on climate scenario and species' range size. On average, SDMs projected onto 4 km climate data predicted 42% more stable habitat (the amount of spatial overlap between predicted current and future climatically suitable habitat) compared with 800 m data. We found only modest agreement between areas predicted to be stable by 90 m models generalized to 4 km grids compared with areas classified as stable based on 4 km models, suggesting that some climate refugia captured at finer scales may be missed using coarser scale data. These differences in projected locations of habitat change may have more serious implications than net habitat area when predictive maps form the basis of conservation decision making.     

Climate‐driven range shifts explain the distribution of extant gene pools and predict future loss of unique lineages in a marine brown alga

The climate‐driven dynamics of species ranges is a critical research question in evolutionary ecology. We ask whether present intraspecific diversity is determined by the imprint of past climate. This is an ongoing debate requiring interdisciplinary examination of population genetic pools and persistence patterns across global ranges. Previously, contrasting inferences and predictions have resulted from distinct genomic coverage and/or geographical information. We aim to describe and explain the causes of geographical contrasts in genetic diversity and their consequences for the future baseline of the global genetic pool, by comparing present geographical distribution of genetic diversity and differentiation with predictive species distribution modelling (SDM) during past extremes, present time and future climate scenarios for a brown alga, Fucus vesiculosus. SDM showed that both atmospheric and oceanic variables shape the global distribution of intertidal species, revealing regions of persistence, extinction and expansion during glacial and postglacial periods. These explained the distribution and structure of present genetic diversity, consisting of differentiated genetic pools with maximal diversity in areas of long‐term persistence. Most of the present species range comprises postglacial expansion zones and, in contrast to highly dispersive marine organisms, expansions involved only local fronts, leaving distinct genetic pools at rear edges. Besides unravelling a complex phylogeographical history and showing congruence between genetic diversity and persistent distribution zones, supporting the hypothesis of niche conservatism, range shifts and loss of unique genetic diversity at the rear edge were predicted for future climate scenarios, impoverishing the global gene pool.     

Combining physiological threshold knowledge to species distribution models is key to improving forecasts of the future niche for macroalgae

Species distribution models (SDM) are a useful tool for predicting species range shifts in response to global warming. However, they do not explore the mechanisms underlying biological processes, making it difficult to predict shifts outside the environmental gradient where the model was trained. In this study, we combine correlative SDMs and knowledge on physiological limits to provide more robust predictions. The thermal thresholds obtained in growth and survival experiments were used as proxies of the fundamental niches of two foundational marine macrophytes. The geographic projections of these species’ distributions obtained using these thresholds and existing SDMs were similar in areas where the species are either absent‐rare or frequent and where their potential and realized niches match, reaching consensus predictions. The cold‐temperate foundational seaweed Himanthalia elongata was predicted to become extinct at its southern limit in northern Spain in response to global warming, whereas the occupancy of southern‐lusitanic Bifurcaria bifurcata was expected to increase. Combined approaches such as this one may also highlight geographic areas where models disagree potentially due to biotic factors. Physiological thresholds alone tended to over‐predict species prevalence, as they cannot identify absences in climatic conditions within the species’ range of physiological tolerance or at the optima. Although SDMs tended to have higher sensitivity than threshold models, they may include regressions that do not reflect causal mechanisms, constraining their predictive power. We present a simple example of how combining correlative and mechanistic knowledge provides a rapid way to gain insight into a species’ niche resulting in consistent predictions and highlighting potential sources of uncertainty in forecasted responses to climate change.     

Uncertainty in predictions of range dynamics: black grouse climbing the Swiss Alps

Empirical species distribution models (SDMs) constitute often the tool of choice for the assessment of rapid climate change effects on species’ vulnerability. Conclusions regarding extinction risks might be misleading, however, because SDMs do not explicitly incorporate dispersal or other demographic processes. Here, we supplement SDMs with a dynamic population model 1) to predict climate‐induced range dynamics for black grouse in Switzerland, 2) to compare direct and indirect measures of extinction risks, and 3) to quantify uncertainty in predictions as well as the sources of that uncertainty. To this end, we linked models of habitat suitability to a spatially explicit, individual‐based model. In an extensive sensitivity analysis, we quantified uncertainty in various model outputs introduced by different SDM algorithms, by different climate scenarios and by demographic model parameters. Potentially suitable habitats were predicted to shift uphill and eastwards. By the end of the 21st century, abrupt habitat losses were predicted in the western Prealps for some climate scenarios. In contrast, population size and occupied area were primarily controlled by currently negative population growth and gradually declined from the beginning of the century across all climate scenarios and SDM algorithms. However, predictions of population dynamic features were highly variable across simulations. Results indicate that inferring extinction probabilities simply from the quantity of suitable habitat may underestimate extinction risks because this may ignore important interactions between life history traits and available habitat. Also, in dynamic range predictions uncertainty in SDM algorithms and climate scenarios can become secondary to uncertainty in dynamic model components. Our study emphasises the need for principal evaluation tools like sensitivity analysis in order to assess uncertainty and robustness in dynamic range predictions. A more direct benefit of such robustness analysis is an improved mechanistic understanding of dynamic species’ responses to climate change.     

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.     

To Be Or Not to Be? Variable selection can change the projected fate of a threatened species under future climate

Species distribution models (SDMs) are commonly used to project future changes in the geographic ranges of species, to estimate extinction rates and to plan biodiversity conservation. However, these models can produce a range of results depending on how they are parameterized, and over‐reliance on a single model may lead to overconfidence in maps of future distributions. The choice of predictor variable can have a greater influence on projected future habitat than the range of climate models used. We demonstrate this in the case of the Ptunarra Brown Butterfly, a species listed as vulnerable in Tasmania, Australia. We use the Maxent model to develop future projections for this species based on three variable sets; all 35 commonly used so‐called ‘bioclimatic’ variables, a subset of these based on expert knowledge, and a set of monthly climate variables relevant to the species’ primary activity period. We used a dynamically downscaled regional climate model based on three global climate models. Depending on the choice of variable set, the species is projected either to experience very little contraction of habitat or to come close to extinction by the end of the century due to lack of suitable climate. The different conclusions could have important consequences for conservation planning and management, including the perceived viability of habitat restoration. The output of SDMs should therefore be used to define the range of possible trajectories a species may be on, and ongoing monitoring used to inform management as changes occur.     

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.     

Disregarding the edaphic dimension in species distribution models leads to the omission of crucial spatial information under climate change: the case of Quercus pubescens in France

Species distribution modelling is an easy, persuasive and useful tool for anticipating species distribution shifts under global change. Numerous studies have used only climate variables to predict future potential species range shifts and have omitted environmental factors important for determining species distribution. Here, we assessed the importance of the edaphic dimension in the niche‐space definition of Quercus pubescens and in future spatial projections under global change over the metropolitan French forest territory. We fitted two species distribution models (SDM) based on presence/absence data (111 013 plots), one calibrated from climate variables only (mean temperature of January and climatic water balance of July) and the other one from both climate and edaphic (soil pH inferred from plants) variables. Future predictions were conducted under two climate scenarios (PCM B2 and HadCM3 A2) and based on 100 simulations using a cellular automaton that accounted for seed dispersal distance, landscape barriers preventing migration and unsuitable land cover. Adding the edaphic dimension to the climate‐only SDM substantially improved the niche‐space definition of Q. pubescens, highlighting an increase in species tolerance in confronting climate constraints as the soil pH increased. Future predictions over the 21st century showed that disregarding the edaphic dimension in SDM led to an overestimation of the potential distribution area, an underestimation of the spatial fragmentation of this area, and prevented the identification of local refugia, leading to an underestimation of the northward shift capacity of Q. pubescens and its persistence in its current distribution area. Spatial discrepancies between climate‐only and climate‐plus‐edaphic models are strengthened when seed dispersal and forest fragmentation are accounted for in predicting a future species distribution area. These discrepancies highlight some imprecision in spatial predictions of potential distribution area of species under climate change scenarios and possibly wrong conclusions for conservation and management perspectives when climate‐only models are used.     

Palaeoclimatic models help to understand current distribution of Caucasian forest species

Spatial and temporal constraints on dispersal explain the absence of species from areas with potentially suitable conditions. Previous studies have shown that post‐glacial recolonization has shaped the current ranges of many species, yet it is not completely clear to what extent interspecific differences in range size depend on different dispersal rates. The inferred boundaries of glacial refugia are difficult to validate, and may bias spatial distribution models (SDMs) that consider post‐glacial dispersal constraints. We predicted the current distribution of 12 Caucasian forest plants and animals, factoring in the effective geographical distance from inferred glacial refugia as an additional predictor. To infer glacial refugia, we tested the transferability of the current SDMs based on the distribution of climatic variables, and projected the most transferable ones onto two climate scenarios simulated for the Last Glacial Maximum (LGM). We then calculated least‐cost distances from the inferred refugia, using elevation as a friction surface, and recalculated the current SDMs incorporating the distances as an additional variable. We compared the predictive powers of the initial with the final SDMs. The palaeoclimatic simulation that best matched the distribution of species was assumed to represent the closest fit to the true palaeoclimate. SDMs incorporating refugial distance performed significantly better for all but one studied species, and the Model for Interdisciplinary Research on Climate (MIROC) climatic simulation provided a more convincing pattern of the LGM climate than the Community Climate System Model (CCSM) simulation. Our results suggest that the projection of suitable habitat models onto past climatic conditions may yield realistic boundaries of glacial refugia, and that the current distribution of forest species in the study region is strongly associated with locations of former refugia. We inferred six major forest refugia throughout western Asia: (1) Colchis; (2) western Anatolia; (3) western Taurus; (4) the upper reaches of the Tigris River; (5) the Levant; and (6) the southern Caspian basin. The boundaries of the modelled refugia were substantially broader than the refugia boundaries inferred solely from pollen records. Thus, our method could be used to: (1) improve models of current species distributions by considering the dispersal histories of the species; and (2) validate alternative reconstructions of palaeoclimate with current distribution data. © 2011 The Linnean Society of London, Biological Journal of the Linnean Society, 2012, 105 , 231–248.     

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.     

The Combined Use of Correlative and Mechanistic Species Distribution Models Benefits Low Conservation Status Species

Species can respond to climate change by tracking appropriate environmental conditions in space, resulting in a range shift. Species Distribution Models (SDMs) can help forecast such range shift responses. For few species, both correlative and mechanistic SDMs were built, but allis shad (Alosa alosa), an endangered anadromous fish species, is one of them. The main purpose of this study was to provide a framework for joint analyses of correlative and mechanistic SDMs projections in order to strengthen conservation measures for species of conservation concern. Guidelines for joint representation and subsequent interpretation of models outputs were defined and applied. The present joint analysis was based on the novel mechanistic model GR3D (Global Repositioning Dynamics of Diadromous fish Distribution) which was parameterized on allis shad and then used to predict its future distribution along the European Atlantic coast under different climate change scenarios (RCP 4.5 and RCP 8.5). We then used a correlative SDM for this species to forecast its distribution across the same geographic area and under the same climate change scenarios. First, projections from correlative and mechanistic models provided congruent trends in probability of habitat suitability and population dynamics. This agreement was preferentially interpreted as referring to the species vulnerability to climate change. Climate change could not be accordingly listed as a major threat for allis shad. The congruence in predicted range limits between SDMs projections was the next point of interest. The difference, when noticed, required to deepen our understanding of the niche modelled by each approach. In this respect, the relative position of the northern range limit between the two methods strongly suggested here that a key biological process related to intraspecific variability was potentially lacking in the mechanistic SDM. Based on our knowledge, we hypothesized that local adaptations to cold temperatures deserved more attention in terms of modelling, but further in conservation planning as well.     

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.     

Selecting from correlated climate variables: a major source of uncertainty for predicting species distributions under climate change

Correlative species distribution models are frequently used to predict species’ range shifts under climate change. However, climate variables often show high collinearity and most statistical approaches require the selection of one among strongly correlated variables. When causal relationships between species presence and climate parameters are unknown, variable selection is often arbitrary, or based on predictive performance under current conditions. While this should only marginally affect current range predictions, future distributions may vary considerably when climate parameters do not change in concert. We investigated this source of uncertainty using four highly correlated climate variables together with a constant set of landscape variables in order to predict current (2010) and future (2050) distributions of four mountain bird species in central Europe. Simulating different parameterization decisions, we generated a) four models including each of the climate variables singly, b) a model taking advantage of all variables simultaneously and c) an un‐weighted average of the predictions of a). We compared model accuracy under current conditions, predicted distributions under four scenarios of climate change, and – for one species – evaluated back‐projections using historical occurrence data. Although current and future variable‐correlations remained constant, and the models’ accuracy under contemporary conditions did not differ, future range predictions varied considerably in all climate change scenarios. Averaged models and models containing all climate variables simultaneously produced intermediate predictions; the latter, however, performed best in back‐projections. This pattern, consistent across different modelling methods, indicates a benefit from including multiple climate predictors in ambiguous situations. Variable selection proved to be an important source of uncertainty for future range predictions, difficult to control using contemporary information. Small, but diverging changes of climate variables, masked by constant overall correlation patterns, can cause substantial differences between future range predictions which need to be accounted for, particularly when outcomes are intended for conservation decisions.