Predicting the impacts of climate change on the distribution of species: are bioclimate envelope models useful?

Modelling strategies for predicting the potential impacts of climate change on the natural distribution of species have often focused on the characterization of a species’ bioclimate envelope. A number of recent critiques have questioned the validity of this approach by pointing to the many factors other than climate that play an important part in determining species distributions and the dynamics of distribution changes. Such factors include biotic interactions, evolutionary change and dispersal ability. This paper reviews and evaluates criticisms of bioclimate envelope models and discusses the implications of these criticisms for the different modelling strategies employed. It is proposed that, although the complexity of the natural system presents fundamental limits to predictive modelling, the bioclimate envelope approach can provide a useful first approximation as to the potentially dramatic impact of climate change on biodiversity. However, it is stressed that the spatial scale at which these models are applied is of fundamental importance, and that model results should not be interpreted without due consideration of the limitations involved. A hierarchical modelling framework is proposed through which some of these limitations can be addressed within a broader, scale‐dependent context.     

Inclusion of local environmental conditions alters high-latitude vegetation change predictions based on bioclimatic models

Current predictions of how species will respond to climate change are typically based on coarse-grained climate surfaces utilizing bioclimate envelope modelling. However, the suitability of environmental conditions for a given species might result from a variety of factors including some unrelated to climate. To address this issue, we investigated whether the inclusion of topographical and soil information in bioclimatic envelope models would significantly alter predictions of climate change—induced fine-scale tree and shrub species range size changes at the tree-limit in subarctic Europe. Using generalized additive models and data on current climate and species distributions and three different climate scenarios for the period 2040–2069, we developed predictions of the currently suitable area and potential range size changes of seven tree and shrub species in an area of 1,100 km² at a resolution of 1-ha. The inclusion of topography and soil information increased the predictive accuracy of climate-only models for all studied species. The predicted changes in species distribution volumes were contradictory, and the predicted occurrences varied greatly depending on the model used. Our results therefore support the arguments that vegetation responses to climate change can be influenced by local environmental conditions and that attention should be paid to the combined effects of these factors. We conclude that disregarding local topography and soil conditions in bioclimatic models may result in biased projections of range expansions and the associated colonization, extinction and turnover assessments.     

Uncertainty of bioclimate envelope models based on the geographical distribution of species

Aim We explored the effects of prevalence, latitudinal range and spatial autocorrelation of species distribution patterns on the accuracy of bioclimate envelope models of butterflies. Location Finland, northern Europe. Methods The data of a national butterfly atlas survey (NAFI) carried out in 1991–2003 with a resolution of 10 × 10 km were used in the analyses. Generalized additive models (GAM) were constructed, for each of 98 species, to estimate the probability of occurrence as a function of climate variables. Model performance was measured using the area under the curve (AUC) of a receiver operating characteristic (ROC) plot. Observed differences in modelling accuracy among species were related to the species’ geographical attributes using multivariate GAM. Results Accuracies of the climate–butterfly models varied from low to very high (AUC values 0.59–0.99), with a mean of 0.79. The modelling performance was related negatively to the latitudinal range and prevalence, and positively to the spatial autocorrelation of the species distribution. These three factors accounted for 75.2% of the variation in the modelling accuracy. Species at the margin of their range or with low prevalence were better predicted than widespread species, and species with clumped distributions better than scattered dispersed species. Main conclusions The results from this study indicate that species’ geographical attributes highly influence the behaviour and uncertainty of species–climate models, which should be taken into account in biogeographical modelling studies and assessments of climate change impacts.     

Model-based uncertainty in species range prediction

Aim Many attempts to predict the potential range of species rely on environmental niche (or ‘bioclimate envelope’) modelling, yet the effects of using different niche‐based methodologies require further investigation. Here we investigate the impact that the choice of model can have on predictions, identify key reasons why model output may differ and discuss the implications that model uncertainty has for policy‐guiding applications. Location The Western Cape of South Africa. Methods We applied nine of the most widely used modelling techniques to model potential distributions under current and predicted future climate for four species (including two subspecies) of Proteaceae. Each model was built using an identical set of five input variables and distribution data for 3996 sampled sites. We compare model predictions by testing agreement between observed and simulated distributions for the present day (using the area under the receiver operating characteristic curve (AUC) and kappa statistics) and by assessing consistency in predictions of range size changes under future climate (using cluster analysis). Results Our analyses show significant differences between predictions from different models, with predicted changes in range size by 2030 differing in both magnitude and direction (e.g. from 92% loss to 322% gain). We explain differences with reference to two characteristics of the modelling techniques: data input requirements (presence/absence vs. presence‐only approaches) and assumptions made by each algorithm when extrapolating beyond the range of data used to build the model. The effects of these factors should be carefully considered when using this modelling approach to predict species ranges. Main conclusions We highlight an important source of uncertainty in assessments of the impacts of climate change on biodiversity and emphasize that model predictions should be interpreted in policy‐guiding applications along with a full appreciation of uncertainty.     

Predicting species responses to climate change: demography and climate microrefugia in California valley oak (Quercus lobata)

Anticipating species movement under climate change is a major focus in conservation. Bioclimate models are one of the few predictive tools for adaptation planning, but are limited in accounting for (i) climatic tolerances in preadult life stages that are potentially more vulnerable to warming; and (ii) local‐scale movement and use of climatic refugia as an alternative or complement to large‐scale changes in distribution. To assess whether these shortfalls can be addressed with field demographic data, we used California valley oak (Quercus lobata Nee), a long‐lived species with juvenile life stages known to be sensitive to climate. We hypothesized that the valley oak bioclimate model, based on adults, would overpredict the species' ability to remain in the projected persisting area, due to higher climate vulnerability of young life stages; and underpredict the potential for the species to remain in the projected contracting area in local‐scale refugia. We assessed the bioclimate model projections against actual demographic patterns in natural populations. We found that saplings were more constricted around surface water than adults in the projected contracting area. We also found that the climate envelope for saplings is narrower than that for adults. Saplings disappeared at a summer maximum temperature 3 °C below that associated with adults. Our findings indicate that rather than a complete shift northward and upward, as predicted by the species bioclimate model, valley oaks are more likely to experience constriction around water bodies, and eventual disappearance from areas exceeding a threshold of maximum temperature. Ours is the first study we know of to examine the importance of discrete life stage climate sensitivities in determining bioclimate modeling inputs, and to identify current climate change‐related constriction of a species around microrefugia. Our findings illustrate that targeted biological fieldwork can be central to understanding climate change‐related movement for long‐lived, sessile species.     

The role of land cover in bioclimatic models depends on spatial resolution

Aim We explored the importance of climate and land cover in bird species distribution models on multiple spatial scales. In particular, we tested whether the integration of land cover data improves the performance of pure bioclimatic models. Location Finland, northern Europe. Methods The data of the bird atlas survey carried out in 1986–89 using a 10 × 10 km uniform grid system in Finland were employed in the analyses. Land cover and climatic variables were compiled using the same grid system. The dependent and explanatory variables were resampled to 20‐km, 40‐km and 80‐km resolutions. Generalized additive models (GAM) were constructed for each of the 88 land bird species studied in order to estimate the probability of occurrence as a function of (1) climate and (2) climate and land cover variables. Model accuracy was measured by a cross‐validation approach using the area under the curve (AUC) of a receiver operating characteristic (ROC) plot. Results In general, the accuracies of the 88 bird–climate models were good at all studied resolutions. However, the inclusion of land cover increased the performance of 79 and 78 of the 88 bioclimatic models at 10‐km and 20‐km resolutions, respectively. There was no significant improvement at the 40‐km resolution. In contrast to the finer resolutions, the inclusion of land cover variables decreased the modelling accuracy at 80km resolution. Main conclusions Our results suggest that the determinants of bird species distributions are hierarchically structured: climatic variables are large‐scale determinants, followed by land cover at finer resolutions. The majority of the land bird species in Finland are rather clearly correlated with climate, and bioclimate envelope models can provide useful tools for identifying the relationships between these species and the environment at resolutions ranging from 10 km to 80 km. However, the notable contribution of land cover to the accuracy of bioclimatic models at 10–20‐km resolutions indicates that the integration of climate and land cover information can improve our understanding and model predictions of biogeographical patterns under global change.     

植物种分布的模拟研究进展

从植物种水平研究植被与气候的关系一直是生态学的热点之一。该文综述了植物种与气候关系的早期研究历史和国内外近期研究进展,尤其是20世纪80年代以来,随着全球变化研究的不断发展和深入,植物种地理分布与气候因子关系研究的最新发展,汇总了最近20年来国际上模拟预测植物种潜在地理分布的模型,比较了不同模型的优缺点。统计模型主要包括以生物气候分室模型或气候分室模型为代表的相关模型、以广义线性模型和广义加性模型为代表的回归模型、以分类和回归树分析及人工神经网络为代表的基于规则的模型、以及生态位模型、气候响应面模型等。机理模型主要介绍了基于BIOME1生物地理模型和FORSKA林窗模型的STASH模型、基于过程的物候模型PHENOFIT,以及一种基于水分平衡、温度和植物物候现象的模型。总结不同模型模拟预测的不同地区植物种未来分布的格局,并介绍中国植物种潜在分布区及未来变化的模拟预测工作,从而为更加准确地模拟预测植物种在未来全球变化情景下的变化趋势提供背景知识。     

Whole range and regional‐based ecological niche models predict differing exposure to 21st century climate change in the key cool temperate rainforest tree southern beech (Nothofagus cunninghamii)

The warmer and drier climates projected for the mid‐ to late‐21st century may have particularly adverse impacts on the cool temperate rainforests of southeastern Australia. Southern beech (Nothofagus cunninghamii; Nothofagaceae), a dominant tree species in these forests, may be vulnerable to minor changes in its climate envelope, especially at the edge of the species range, with Holocene fossil evidence showing local extinction of populations in response to small climate changes. We modelled the stability of this species climate envelope using the maximum entropy algorithm implemented in Maxent and two thresholds of presence/absence by projecting the modern climate envelope onto four Global Circulation Models forecasted for two time periods (2050s and 2070s). The climate envelope, as estimated from the species present climatic range, is predicted to shrink by up to 49% by the 2050s and up to 64% by the 2070s. The greatest predicted reduction is in Victoria with 91–100% of its current range being climatically unsuitable by the 2070s. Climatically similar areas to the species present range are predicted to remain in mountainous areas of western Tasmania, the Northeast Highlands of Tasmania, and the Baw Baw Plateau in the Central Highlands of Victoria. However, region‐specific modelling approaches made very different predictions from the whole‐range based models, especially in the severity of the predicted decline for Victorian populations of N. cunninghamii which occur in much warmer climates than the rest of the species geographical range. This shows that, for widespread species that span a range of climate zones, the exposure of current populations to climate change may be better modelled using a regional based approach. How the species responds to climate change will depend on the species ability to respond to drier and warmer climates and the concomitant increase in fire intensity.     

Predicting potential climate change impacts with bioclimate envelope models: a palaeoecological perspective

Aim We assess the realism of bioclimate envelope model projections for anticipated future climates by validating ecosystem reconstructions for the late Quaternary with fossil and pollen data. Specifically, we ask: (1) do climate conditions with no modern analogue negatively affect the accuracy of ecosystem reconstructions? (2) are bioclimate envelope model projections biased towards under‐predicting forested ecosystems? (3) given a palaeoecological perspective, are potential habitat projections for the 21st century within model capabilities? Location Western North America. Methods We used an ensemble classifier modelling approach (RandomForest) to spatially project the climate space of modern ecosystem classes throughout the Holocene (at 6000, 9000, 11,000, 14,000, 16,000, and 21,000 YBP) using palaeoclimate surfaces generated by two general circulation models (GFDL and CCM1). The degree of novel arrangement of climate variables was quantified with the multivariate Mahalanobis distance to the nearest modern climatic equivalent. Model projections were validated against biome classifications inferred from 1460 palaeoecological records. Results Model accuracy assessed against independent palaeoecology data is generally low for the present day, increases for 6000 YBP, and then rapidly declines towards the last glacial maximum, primarily due to the under‐prediction of forested biomes. Misclassifications were closely correlated with the degree of climate dissimilarity from the present day. For future projections, no‐analogue climates unexpectedly emerged in the coastal Pacific Northwest but were absent throughout the rest of the study area. Main conclusions Bioclimate envelope models could approximately reconstruct ecosystem distributions for the mid‐ to late‐Holocene but proved unreliable in the Late Pleistocene. We attribute this failure to a combination of no‐analogue climates and a potential lack of niche conservatism in tree species. However, climate dissimilarities in future projections are comparatively minor (similar to those of the mid‐Holocene), and we conclude that no‐analogue climates should not compromise the accuracy of model predictions for the next century.     

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.     

Does seasonal fine‐tuning of climatic variables improve the performance of bioclimatic envelope models for migratory birds?

We examined the influence of 'seasonal fine-tuning' of climatic variables on the performance of bioclimatic envelope models of migrating birds. Using climate data and national bird atlas data from a 10 × 10 km uniform grid system in Finland, we tested whether the replacement of one 'baseline' set of variables including summer (June–August) temperature and precipitation variables with climate variables tailored ('fine-tuned') for each species individually improved the bird-climate models. The fine-tuning was conducted on the basis of time of arrival and early breeding of the species. Two generalized additive models (GAMs) were constructed for each of the 63 bird species studied, employing (1) the baseline climate variables and (2) the fine-tuned climate variables. Model performance was measured as explanatory power (deviance change) and predictive power (area under the curve; AUC) statistics derived from cross-validation. Fine-tuned climate variables provided, in many cases, statistically significantly improved model performance compared to using the same baseline set of variables for all the species. Model improvements mainly concerned bird species arriving and starting their breeding in May–June. We conclude that the use of the fine-tuned climate variables tailored for each species individually on the basis of their arrival and critical breeding periods can provide important benefits for bioclimatic modelling.     

Cetacean range and climate in the eastern North Atlantic: future predictions and implications for conservation

There is increasing evidence that the distributions of a large number of species are shifting with global climate change as they track changing surface temperatures that define their thermal niche. Modelling efforts to predict species distributions under future climates have increased with concern about the overall impact of these distribution shifts on species ecology, and especially where barriers to dispersal exist. Here we apply a bio‐climatic envelope modelling technique to investigate the impacts of climate change on the geographic range of ten cetacean species in the eastern North Atlantic and to assess how such modelling can be used to inform conservation and management. The modelling process integrates elements of a species' habitat and thermal niche, and employs “hindcasting” of historical distribution changes in order to verify the accuracy of the modelled relationship between temperature and species range. If this ability is not verified, there is a risk that inappropriate or inaccurate models will be used to make future predictions of species distributions. Of the ten species investigated, we found that while the models for nine could successfully explain current spatial distribution, only four had a good ability to predict distribution changes over time in response to changes in water temperature. Applied to future climate scenarios, the four species‐specific models with good predictive abilities indicated range expansion in one species and range contraction in three others, including the potential loss of up to 80% of suitable white‐beaked dolphin habitat. Model predictions allow identification of affected areas and the likely time‐scales over which impacts will occur. Thus, this work provides important information on both our ability to predict how individual species will respond to future climate change and the applicability of predictive distribution models as a tool to help construct viable conservation and management strategies.     

Modelling the effects of climate change on the distribution and production of marine fishes: accounting for trophic interactions in a dynamic bioclimate envelope model

Climate change has already altered the distribution of marine fishes. Future predictions of fish distributions and catches based on bioclimate envelope models are available, but to date they have not considered interspecific interactions. We address this by combining the species‐based Dynamic Bioclimate Envelope Model (DBEM) with a size‐based trophic model. The new approach provides spatially and temporally resolved predictions of changes in species' size, abundance and catch potential that account for the effects of ecological interactions. Predicted latitudinal shifts are, on average, reduced by 20% when species interactions are incorporated, compared to DBEM predictions, with pelagic species showing the greatest reductions. Goodness‐of‐fit of biomass data from fish stock assessments in the North Atlantic between 1991 and 2003 is improved slightly by including species interactions. The differences between predictions from the two models may be relatively modest because, at the North Atlantic basin scale, (i) predators and competitors may respond to climate change together; (ii) existing parameterization of the DBEM might implicitly incorporate trophic interactions; and/or (iii) trophic interactions might not be the main driver of responses to climate. Future analyses using ecologically explicit models and data will improve understanding of the effects of inter‐specific interactions on responses to climate change, and better inform managers about plausible ecological and fishery consequences of a changing environment.     

Future vulnerability mapping based on response to extreme climate events: Dieback thresholds in an endemic California oak

Aim

This study presents a bioclimate modelling approach, using responses to extreme climate events, rather than historical distributional associations, to project future species vulnerability and refugia. We aim to illustrate the compounding effects of groundwater loss and climate on species vulnerability.

Location

California, USA.

Methods

As a case study, we used the 2012–2015 California drought and resulting extensive dieback of blue oak (Quercus douglasii). We used aerial dieback surveys, downscaled climate data and subsurface water change data to develop boosted regression tree models identifying key thresholds associated with dieback throughout the blue oak distribution. We (1) combined observed dieback–climatic threshold relationships with climate futures to anticipate future areas of vulnerability and (2) used satellite‐derived measurements of subsurface water loss in drought/dieback modelling to capture the mediating effect of groundwater on species response to climatic drought.

Results

A model including climate, climate anomalies and subsurface water change explained 46% of the variability in dieback. Precipitation in 2015 and subsurface water change accounted for 62.6% of the modelled probability of dieback. We found an interaction between precipitation and subsurface water in which dieback probability increased with low precipitation and subsurface water loss. The relationship between precipitation and dieback was nonlinear, with 99% of dieback occurring in areas that received <363 mm precipitation. Based on a MIROC_rcp85 future climate scenario, relative to historical conditions, 13% of the blue oak distribution is predicted to experience more frequent years below this precipitation threshold by mid‐century and 81% by end of century.

Main conclusions

As ongoing climate change and extreme events impact ecological processes, the identification of thresholds associated with observed dieback may be combined with climate futures to help identify vulnerable populations and refugia and prioritize climate change‐related conservation efforts.     

Can we project changes in fish abundance and distribution in response to climate?

Large‐scale and long‐term changes in fish abundance and distribution in response to climate change have been simulated using both statistical and process‐based models. However, national and regional fisheries management requires also shorter term projections on smaller spatial scales, and these need to be validated against fisheries data. A 26‐year time series of fish surveys with high spatial resolution in the North‐East Atlantic provides a unique opportunity to assess the ability of models to correctly simulate the changes in fish distribution and abundance that occurred in response to climate variability and change. We use a dynamic bioclimate envelope model forced by physical–biogeochemical output from eight ocean models to simulate changes in fish abundance and distribution at scales down to a spatial resolution of 0.5°. When comparing with these simulations with annual fish survey data, we found the largest differences at the 0.5° scale. Differences between fishery model runs driven by different biogeochemical models decrease dramatically when results are aggregated to larger scales (e.g. the whole North Sea), to total catches rather than individual species or when the ensemble mean instead of individual simulations are used. Recent improvements in the fidelity of biogeochemical models translate into lower error rates in the fisheries simulations. However, predictions based on different biogeochemical models are often more similar to each other than they are to the survey data, except for some pelagic species. We conclude that model results can be used to guide fisheries management at larger spatial scales, but more caution is needed at smaller scales.     

Range position and climate sensitivity: The structure of among‐population demographic responses to climatic variation

Species’ distributions will respond to climate change based on the relationship between local demographic processes and climate and how this relationship varies based on range position. A rarely tested demographic prediction is that populations at the extremes of a species’ climate envelope (e.g., populations in areas with the highest mean annual temperature) will be most sensitive to local shifts in climate (i.e., warming). We tested this prediction using a dynamic species distribution model linking demographic rates to variation in temperature and precipitation for wood frogs (Lithobates sylvaticus) in North America. Using long‐term monitoring data from 746 populations in 27 study areas, we determined how climatic variation affected population growth rates and how these relationships varied with respect to long‐term climate. Some models supported the predicted pattern, with negative effects of extreme summer temperatures in hotter areas and positive effects on recruitment for summer water availability in drier areas. We also found evidence of interacting temperature and precipitation influencing population size, such as extreme heat having less of a negative effect in wetter areas. Other results were contrary to predictions, such as positive effects of summer water availability in wetter parts of the range and positive responses to winter warming especially in milder areas. In general, we found wood frogs were more sensitive to changes in temperature or temperature interacting with precipitation than to changes in precipitation alone. Our results suggest that sensitivity to changes in climate cannot be predicted simply by knowing locations within the species’ climate envelope. Many climate processes did not affect population growth rates in the predicted direction based on range position. Processes such as species‐interactions, local adaptation, and interactions with the physical landscape likely affect the responses we observed. Our work highlights the need to measure demographic responses to changing climate.     

Global distribution and bioclimatic characterization of alpine biomes

Although there is a general consensus on the distribution and ecological features of terrestrial biomes, the allocation of alpine ecosystems in the global biogeographic system is still unclear. Here, we delineate a global map of alpine areas above the treeline by modelling regional treeline elevation at 30 m resolution, using global forest cover data and quantile regression. We then used global datasets to 1) assess the climatic characteristics of alpine ecosystems using principal component analysis, 2) define bioclimatic groups by an optimized cluster analysis and 3) evaluate patterns of primary productivity based on the normalized difference vegetation index. As defined here, alpine biomes cover 3.56 Mkm² or 2.64% of land outside Antarctica. Despite temperature differences across latitude, these ecosystems converge below a sharp threshold of 5.9°C and towards the colder end of the global climatic space. Below that temperature threshold, alpine ecosystems are influenced by a latitudinal gradient of mean annual temperature and they are climatically differentiated by seasonality and continentality. This gradient delineates a climatic envelope of global alpine biomes around temperate, boreal and tundra biomes as defined in Whittaker's scheme. Although alpine biomes are similarly dominated by poorly vegetated areas, world ecoregions show strong differences in the productivity of their alpine belt irrespectively of major climate zones. These results suggest that vegetation structure and function of alpine ecosystems are driven by regional and local contingencies in addition to macroclimatic factors.     

The relative importance of climate and habitat in determining the distributions of species at different spatial scales: a case study with ground beetles in Great Britain

Experimental studies have shown that many species show preferences for different climatic conditions, or may die in unsuitable conditions. Climate envelope models have been used frequently in recent years to predict the presence and absence of species at large spatial scales. However, many authors have postulated that the distributions of species at smaller spatial scales are determined by factors such as habitat availability and biotic interactions. Climatic effects are often assumed by modellers to be unimportant at fine resolutions, but few studies have actually tested this. We sampled the distributions of 20 beetle species of the family Carabidae across three study sites by pitfall trapping, and at the national scale from monitoring data. Statistical models were constructed to determine which of two sets of environmental variables (temperature or broad habitat type) best accounted for the observed data at the three sites and at the national (Great Britain) scale. High‐resolution temperature variables frequently produced better models (as determined by AIC) than habitat features when modelling the distributions of species at a local scale, within the three study sites. Conversely, habitat was always a better predictor than temperature when describing species’ distributions at a coarse scale within Great Britain. Northerly species were most likely to occur in cool micro‐sites within the study sites, whereas southerly species were most likely to occur in warm micro‐sites. Effects of microclimate were not limited to species at the edges of their distribution, and fine‐resolution temperature surfaces should therefore ideally be utilised when undertaking climate‐envelope modelling.     

The tick Ixodes ricinus: distribution and climate preferences in the western Palaearctic

In this study, multivariate spatial clustering on monthly normalized difference vegetation index (NDVI) maps is used to classify ecological regions over the western Palaearctic. This classification is then used to delineate the distribution and climate preferences of populations (clades) of the tick Ixodes ricinus L. (Acari: Ixodidae) from a geographically extensive dataset of tick records and a gridded 2.5-km resolution climate dataset. Using monthly layers of the NDVI, regions of similar ecological attributes were defined and nine populations with significant differences in critical climate parameters (P< 0.005) were detected. Grouping of tick records according to other categories, such as political divisions, a 4 degrees x 4 degrees grid overlying the study area, or the CORINE) and USGS) vegetation classification schemes did not provided significantly separated populations (P = 0.094-0.304). Factor analysis and hierarchical tree clustering provided an ecological overview of these tick clades: two Mediterranean and one Scandinavian (western) clades are clearly separated from a node that includes clades of different parts of central Europe and the British Isles, with contrasting affinities between the different clades. The capture records of these ecologically separated clades produce a clear bias when bioclimate envelope modelling is applied to the mapping of habitat suitability for the tick in the western Palaearctic. The best-performing methods (Cohen's kappa = 0.834-0.912) use partial models developed with data from each ecoregion, which are then overlapped over the region of study. It is concluded that the use of ecologically derived ecoregions is an objective step in assessing the presence of ecologically different clades, and provides a guide in the development of data partitioning for habitat suitability modelling.     

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.