Document Detail

Estimating how inflated or obscured effects of climate affect forecasted species distribution.
Jump to Full Text
MedLine Citation:
PMID:  23349726     Owner:  NLM     Status:  MEDLINE    
Climate is one of the main drivers of species distribution. However, as different environmental factors tend to co-vary, the effect of climate cannot be taken at face value, as it may be either inflated or obscured by other correlated factors. We used the favourability models of four species (Alytes dickhilleni, Vipera latasti, Aquila fasciata and Capra pyrenaica) inhabiting Spanish mountains as case studies to evaluate the relative contribution of climate in their forecasted favourability by using variation partitioning and weighting the effect of climate in relation to non-climatic factors. By calculating the pure effect of the climatic factor, the pure effects of non-climatic factors, the shared climatic effect and the proportion of the pure effect of the climatic factor in relation to its apparent effect (ρ), we assessed the apparent effect and the pure independent effect of climate. We then projected both types of effects when modelling the future favourability for each species and combination of AOGCM-SRES (two Atmosphere-Ocean General Circulation Models: CGCM2 and ECHAM4, and two Special Reports on Emission Scenarios (SRES): A2 and B2). The results show that the apparent effect of climate can be either inflated (overrated) or obscured (underrated) by other correlated factors. These differences were species-specific; the sum of favourable areas forecasted according to the pure climatic effect differed from that forecasted according to the apparent climatic effect by about 61% on average for one of the species analyzed, and by about 20% on average for each of the other species. The pure effect of future climate on species distributions can only be estimated by combining climate with other factors. Transferring the pure climatic effect and the apparent climatic effect to the future delimits the maximum and minimum favourable areas forecasted for each species in each climate change scenario.
Raimundo Real; David Romero; Jesús Olivero; Alba Estrada; Ana L Márquez
Related Documents :
16105896 - Oligonucleotide properties determination and primer designing: a critical examination o...
24732386 - Can technology help to reduce underage drinking? evidence from the false id laws with s...
24530806 - Susceptible-infectious-recovered models revisited: from the individual level to the pop...
23024186 - Human demographic processes and genetic variation as revealed by mtdna simulations.
9358266 - Sample size required for predefined linkage decision quality.
20975986 - Recognizing independent and joint activities among multiple residents in smart environm...
Publication Detail:
Type:  Journal Article; Research Support, Non-U.S. Gov't     Date:  2013-01-11
Journal Detail:
Title:  PloS one     Volume:  8     ISSN:  1932-6203     ISO Abbreviation:  PLoS ONE     Publication Date:  2013  
Date Detail:
Created Date:  2013-01-25     Completed Date:  2013-07-09     Revised Date:  2013-07-11    
Medline Journal Info:
Nlm Unique ID:  101285081     Medline TA:  PLoS One     Country:  United States    
Other Details:
Languages:  eng     Pagination:  e53646     Citation Subset:  IM    
Biogeography, Diversity, and Conservation Research Team, Department of Animal Biology, Faculty of Sciences, University of Malaga, Malaga, Spain.
Export Citation:
APA/MLA Format     Download EndNote     Download BibTex
MeSH Terms
Climate Change / statistics & numerical data
Conservation of Natural Resources / statistics & numerical data*
Models, Statistical*

From MEDLINE®/PubMed®, a database of the U.S. National Library of Medicine

Full Text
Journal Information
Journal ID (nlm-ta): PLoS One
Journal ID (iso-abbrev): PLoS ONE
Journal ID (publisher-id): plos
Journal ID (pmc): plosone
ISSN: 1932-6203
Publisher: Public Library of Science, San Francisco, USA
Article Information
Download PDF
Copyright: 2013 Real et al
Received Day: 24 Month: 5 Year: 2012
Accepted Day: 3 Month: 12 Year: 2012
collection publication date: Year: 2013
Electronic publication date: Day: 11 Month: 1 Year: 2013
ecorrected publication date: Day: 18 Month: 1 Year: 2013
Volume: 8 Issue: 1
E-location ID: e53646
PubMed Id: 23349726
ID: 3548625
Publisher Id: PONE-D-12-15314The Pure Effect of Climate on Species Distribution
DOI: 10.1371/journal.pone.0053646

Estimating How Inflated or Obscured Effects of Climate Affect Forecasted Species Distribution
Raimundo Real1
David Romero1
Jesús Olivero1
Alba Estrada2
Ana L. Márquez1*
Kimberly Patraw Van Nieledit1 Role: Editor
1Biogeography, Diversity, and Conservation Research Team, Department of Animal Biology, Faculty of Sciences, University of Malaga, Malaga, Spain
2Instituto de Investigación en Recursos Cinegéticos IREC, (CSIC-UCLM), Ciudad Real, Spain
University of Western Australia, Australia
Correspondence: * E-mail:
[conflict] Competing Interests: The authors have declared that no competing interests exist.
Contributed by footnote: Conceived and designed the experiments: RR DR. Analyzed the data: ALM JO AE. Wrote the paper: ALM RR.


Species distribution models (SDMs) are becoming increasingly important tools for conservation biology, because determining which factors drive the distribution patterns can help to adopt more specific and appropriate strategies for the management and conservation species [1]. This knowledge is also the basis for making good forecasts on the effect of climate change on future species distributions or suitable areas, which is a new challenge for environmental managers [2, 3]. However, the estimation of impact of the climate change on future species distribution is complex and related to different kinds of uncertainties [4–9], including the inability to assess the weight of climate as a driver of species distribution.

Climate envelope models are widely used to forecast future species distributions under climate change scenarios [10–13]. Some authors argue against the validity of using SDMs based on climatic variables alone as tools to forecast future species distributions, because they consider that other factors play a role in the distributions and that these factors are not taken into account in the models [3, 5, 14–16]. Apart from climate, species distributions may be controlled by spatial trends, topography, human activity, biotic interactions, history, and population dynamics, among others [17–20]. As species may show differential responses to these factors [3, 5, 21], their combined importance should be assessed before projecting species distribution models to the future. In addition, the effect of climate can only be assessed in the context of the other influential factors, because its pure effect on species distributions could be obscured or overrated by correlated aspects, becoming evident only when all the relevant factors are considered together [18, 19, 22, 23].

Variation partitioning techniques have been used to separate the effects of different factors on species richness [24, 25], on abundance [26], on ecological communities [27], on species assemblages [28] and also on species distributions [29, 30]. These techniques have also been used to segregate the pure effect of different factors (topography, climate, human activity, spatial situation and lithology) on species distributions [31]. However, these techniques have not been used to relate the pure climatic effect to its apparent effect, being the latter in correlation with other factors. This is of fundamental importance, because the apparent climatic effect could be misrepresenting the true role of climate on species distributions due to the effect of other correlated factors. Therefore, the potential changes in species distributions related to climate change could be distorted and lead to misleading conclusions about the species vulnerability or their risk of extinction.

Mountains are areas of interest regarding the early detection and study of the signals of climate change and its impact on ecological systems [32]. Mountain species are particularly sensitive to climate change [33–36], because mountain areas have more pristine habitats than lowland landscapes, and because these species can track climate change over shorter distances [37, 38]. Predicting the possible effects of climate change on the future distributions of these kinds of species is of fundamental importance in conservation plans.

We evaluated the relative contribution of climate in mainland Spain to the forecasted favourability of four vertebrate mountain species whose distributions are related to climate and to altitude or slope, and for which published distribution models are available both for the present and for the future according to the apparent effect of climate [8, 19]. The aim of this work is to propose a method by which to analyze the relative contribution of climate in relation to non-climatic factors (spatial, topographic, and human) and to distinguish between its apparent and its pure effect in models designed to forecast how favourable areas for species could vary because of climate change. Our results underline the possibly misleading outcome of not considering the pure climatic effect in the projections of the SDMs to the future.

The study area

Mainland Spain is located in southwestern Europe and covers an area of 493,518 km2. Its latitude (40° N), geographical position and complex orography make its climate heterogeneous, with a precipitation gradient (100–2500 mm) decreasing mainly eastward and southward, and a temperature gradient (6°–18°C) decreasing mainly northward [39]. In Spain, five homogeneous climatic precipitation regions [39] can be distinguished: 1) the North Atlantic coast, which has abundant and regular precipitation due to the continuous arrival of Atlantic frontal systems; 2) the central area that receives wet and cold air intrusions from frontal Atlantic systems and presents low precipitations; 3) the eastern coast, which is characterized by irregular and scanty annual precipitation, with large variability due to severe rainfall events produced by wet and warm air intrusions from the Mediterranean Sea [40]; 4) the southeastern region, which is a dry desert-like area with very little rainfall; and 5) the southwestern region, which has more regular and abundant rainfall influenced by Atlantic winds.

This makes Spain particularly appropriate for analyzing the effect of different climate change scenarios on species distributions (e.g., [19]).

The species

We analysed the distribution in mainland Spain of four vertebrate species whose distributions are positively associated with altitude or slope. We chose an amphibian (Baetic midwife toad, Alytes dickhilleni), a reptile (Lataste’s viper, Vipera latasti), a bird (Bonelli’s eagle, Aquila fasciata) and a mammal (Iberian wild goat, Capra pyrenaica). A. dickhilleni is a small toad, between 32.8 and 56.5 mm in length, endemic to Spain and located exclusively in the mountainous systems of the southwestern part of the Iberian Peninsula. It lives on rough and steep terrains, in cracks and crevices next to streams, springs and pools. Species reproduction occurs in permanent water points. V. latasti is a venomous viper species found in southwestern Europe and northwestern Africa, which can reach 70 cm in length. Its distribution in Spain is relegated by human activity to mountainous and sparsely populated areas. A. fasciata is a small to medium-size eagle, 55–65 cm in length. It is one of the rarest raptors in Europe and is a priority target species for special conservation measures in Spain (Council Directive 2009/147/EC and National Real Decreto 439/1990). It occupies mountain ranges, small hills and plains, where it breeds mainly in cliffs. C. pyrenaica is an endemic goat, only found in the mountainous areas of Spain. It is a species with strong sexual dimorphism, males are larger than females and their horns are three times longer and thicker than those of females. It lives in both forests and grassy expanses in mountains at altitudes between 500 and 2500 meters.

Baseline models

As starting point, we used the current favourability models obtained for A. dickhilleni, V. latasti, A. fasciata and C. pyrenaica using different climate change scenarios in mainland Spain, available in Márquez et al. [19]. The distributions of the four species were modelled using variables related to climate, spatial situation, topography, and human activity (see Table S1 for more specific details of the variables) and with the favourability function as the modelling technique [41–43]. Climatic variables were provided by the Agencia Estatal de Meteorología (AEMET), which regionalized the general circulation models CGCM2 (Canadian Climate Centre for Modeling) and ECHAM4 (Max Planck Institut für Meteorologie) and the A2 and B2 emission scenarios for Spain. Mean values of the climatic variables were obtained for the periods: 1961–1990, 2011–2040, 2041–2070 and 2071–2100 (Figure S1). For each species a total of four factor models were obtained, related to climate alone, spatial situation alone, topography alone, and human activity alone, respectively. A combined model was then obtained which took into account the four environmental factors (climatic and non-climatic) together (see Márquez et al. [19] for more details). In this way, the favourability values for each species in each cell at the present time (Fp) were obtained. Future favourability values for each species according to the apparent effect of climate (FfClim) in each cell, as well as an analysis of the impact of climate model choice and scenario choice on expected favourability, were obtained by Real et al. [8] by replacing the current (1961–1990) climatic values in the combined favourability models with those expected according to each AOGCM and SRES for the following time periods (2011–2040, 2041–2070, 2071–2100).

Variation partitioning

We segregated the pure effect of climate from the effect of the other factors in the models using a variation partitioning procedure similar to that of Borcard et al. [44], Barbosa et al. [29] and Muñoz et al. [45] with some modifications. Thus, we specified how much of the variation of the combined favourability model was explained by the pure effect of climate (not affected by the covariation with other factors in the model), and which proportion of the climate effect cannot be distinguished from that of the other factors (shared effect) [45].

The portion of the variation in the model apparently explained by climate was estimated using the coefficient of determination of the linear regression of the logit function (y) of the model on the climatic variables included in it (R2Clim) for the period 1961–1990; the part apparently explained by the non-climatic factors (R2NClim) was obtained in a similar manner. The pure effect of climate (R2pClim) was then assessed by subtracting from 1 the variation of the combined model explained by the non climatic factors (R2pClim  =  1 - R2NClim). The pure effect of the non-climatic factors was obtained by subtracting from 1 the variation explained by climate (R2pNClim  =  1 - R2Clim). The effect shared by climate and non-climatic factors was obtained by subtracting from 1 both pure effects (R2ClimNClim  =  1 - R2pClim - R2pNClim). We used the adjusted R2 in all cases [46] although, given the high number of cells used (n = 5167 10x10 km2 UTM cells), the difference between R2 and adjusted R2 was very small. This partitioning procedure was applied only to the portion of variation explained by the model, not over the total variation of the species distribution [44, 46], as it is the explanatory model which is calibrated to be transferred to the future. We expressed these effects as percentages and considered them to be the percentage of model variation attributable to the pure climatic effect (PCE), the pure non-climatic effect (PNCE) and the shared effect of climate and non-climatic factors (SCE).

We estimated the proportion of the apparent climatic effect represented by the pure effect of climate as . We calculated the logit function expected for the future in each cell according to the pure effect of climate (yfClim) by applying the expression , where yp is the logit function at the present time and yf is the logit function expected for the future according to the apparent effect of climate. We obtained the future favourability according to the pure climatic effect (FfPClim) using the expression:

[Formula ID: pone.0053646.e004]

where n1 is the number of presences and n0 the number of absences (see formula 7 in Real et al. [41]).

This way of inferring the effect of climate differs from usual projections according to climate change scenarios, which are customarily based on the apparent effect of climate. The difference between the sum of areas forecasted to be favourable according to the apparent climatic effect (FAC), calculated using the expression , and those forecasted according to the pure climatic effect (FPC), calculated using the expression , was computed and expressed as the relative proportion of discrepancy R = (FAC-FPC )/FAC.


Climate had a more significant effect than non-climatic factors on A. fasciata (Table 1). However, for the other species (A. dickhilleni, V. latasti and C. pyrenaica), the non-climatic effect was more important than the climatic effect (Table 1). In addition, regarding A. dickhilleni, the percentage of variation of the model attributable to a Shared Climatic Effect (SCE) was very important, which means that the apparent effect of climate could be due in large part to other correlated factors. Regarding A. fasciata and C. pyrenaica, the values of SCE were negative in the majority of the models. These negative values measure the amount of reciprocal obscuring caused by factors that have opposite geographic effects on the explained favourability, so that in these species the apparent climatic effect under-represents the pure climatic effect.

Figure 1 shows three examples of the differences between forecasted favourabilities taking into account the apparent climatic effect (FAC) and the pure climate effect (FPC), when FPC is lower, similar or higher than FAC. The differences between the two future forecasted favourabilities (FAC-FPC) and the relative proportion of discrepancy between both types of effects (R =  (FAC-FPC) / FAC) are shown in Table 2. Figure 2 shows the spatial distribution of the difference between the forecasted favourabilities taking into account the apparent climatic effect and the pure climatic effect (|FfClim-FfPClim|) for the species and situations described in Figure 1.


The inclusion of climatic and non-climatic factors in SDMs is recommended not only because it can improve fit and increase their predictive accuracy [47, 48], but also because the effect of climate can only be assessed in the context of the other influential factors [18, 19, 22, 23]. Our results show that the correlation of influential non-climatic factors with temperature and precipitation could either inflate or obscure the apparent effect of climate, and that this modification of the apparent effect of climate would remain hidden if non-climatic variables were not included in the SDM. Even the use of the latitude and longitude of every cell alone may pinpoint certain areas of origin, dispersion, or past vicariance events driving current distributions, which results in a historically-caused spatial pattern that may coincide with specific climatic characteristics [49]. Consequently, the true effect of climate should be assessed in the context of spatial influences both on species distributions and on climate [50]. We used human, topographic and spatial variables as non-climatic predictors that, although correlated with climatic variables, can influence species distributions for reasons not directly linked to climate [18]. It was by taking into account these non-climatic factors and removing their effects statistically that we identified the underlying pure climate-distribution relationships, which could then be used in forecasting their distribution shifts under climate change [3].

However, the inclusion of climatic variables together with non-climatic, static variables entails other kinds of problems. Stanton et al. [51] considered that static variables such as elevation, latitude or longitude may hinder the accuracy of future predictions, as the relationships between them and climatic variables is likely to change in the future, and that including such variables in the SDM is likely to result in models which underestimate the effects of climate change. Our results confirm that this may be the case, although these effects may be under- or over-estimated. Our procedure is a way to gauge these relationships and assess the maximum extent to which the current correlation between these static variables and climate may affect the climatic parameters in the SDM.

On the other hand, the shared climatic effect is equally attributable to climate or to other correlated factors, so in our current state of knowledge the exact effect of climate cannot be determined with precision, although it lies somewhere between the apparent effect and the pure effect. The uncertainty related to the differences between both effects vary spatially and their intensity depends on the species (see Figure 2). This kind of uncertainty is added to other sources of uncertainty associated with forecasting future species distributions [8, 52, 53], among them those derived from assuming that the species' climate tolerances will remain constant through time, which is one serious limitation to the customary use of SDMs. However, despite these uncertainties, the estimation of species range shifts is the basis to predict where the species are likely to move under different future conditions [54]. More reliable predictions of species distribution responses to future climate conditions depend on developing more rigorous statistical analyses of the available data and on the combination of different factors, as well as on placing limits on the different uncertainties involved in the scientific forecasting of future events [19]. Transferring the pure climatic effect and the apparent climatic effect to the future allows us to delimit the maximum and minimum effect of climate on the species distributions.

Most of the favourability models of all the species considered in this study included three or four factors (spatial situation, topography, human activity and climate) (see Table 3 in Márquez et al. [19]) that we summarized into a climatic factor and a non-climatic factor. Pure climatic effect (PCE), pure non-climatic effect (PNCE) and shared climatic effect (SCE) on the future favourable areas for A. dickhilleni, V. latasti, A. fasciata and C. pyrenaica differed substantially (see Table 1). For A. dickhilleni, V. latasti and C. pyrenaica the PNCE was more important than the PCE, which suggests that their future distributions will be more related to non-climatic environmental variables, such as biotic interactions, past human activities or past contingent events, than with the climatic factor [31, 55]. For A. fasciata, and in most scenarios for C. pyrenaica, the SCE was negative, that is, the climatic effect and the non-climatic effect can be reciprocally obscured by their opposite effect on the explained favourability (Table 1) [24, 46]. In these cases the apparent climatic effect under-represents the pure climatic effect. Future favourabilitiy for these cases taking into account the pure climatic effect would represent their maximum future favourable area, which is higher than that forecasted according to the apparent climatic effect (Table 2). The future areas favourable to A. dickhilleni forecasted according to the apparent climatic effect differed from those forecasted according to the pure climatic effect by 61% on average, which was the highest difference in the four species considered (Table 2). In this case, the apparent climatic effect was highly inflated by non-climatic factors.

The SCE could be a measure of uncertainty related to the complex interactions existing between climate and non-climatic factors. In some cases it represents the obscured climatic effect (when it has a negative value) and in other cases the inflated climatic effect (when it has a positive value). In any case, the SCE represents the uncertainty associated with the possibility of misunderstanding the effect of climate due to the effect of other correlated factors. This improves the usefulness of this kind of model for understanding species’ potential responses to climate change, although possible changes in species-environment correlations through time can, nevertheless, place a limit on the predictive performance of these models [56]. According to Pearson and Dawson [10], understanding the complex interaction between the many factors affecting distributions is needed for the performance of more realistic simulations of the effect of climate change on species distributions. Models that take the effect of climate at face value yield future potential favourable areas that are, depending on the species, overestimated or underestimated. Using the method proposed in this paper, models may more realistically assess the levels of potential threat or opportunities to species of climate change.


In contrast to the tendency of not using correlated variables in spatial distribution models due to the possibility of the resulting coefficients can being unstable, we have to deal with the fact that in nature most factors are correlated; thus analyses that separate the pure and combined effect of the relevant factors should be performed. Given that the apparent effect of climate can be either inflated or obscured by other correlated factors, transferring both the pure climatic effect and the apparent climatic effect to the future allows us to delimit the maximum and minimum favourable areas forecasted for each species in each climate change scenario, thus permitting us to assess the uncertainty associated with the possibility of misrepresenting the effect of climate. This also allows us to detect and control the over- or under-estimation of the effect of climate change (either positive or negative) on future species distributions that is implicit in current climate envelope models. This may make models more complex and harder to perform, but the output would be closer to what may be scientifically forecasted taking into account this kind of uncertainty.

Supporting Information Figure S1

Precipitations and temperatures. Annual precipitations and mean annual maximum temperatures for each period and combination of circulation model and scenario used.


Click here for additional data file (pone.0053646.s001.png)

Table S1

Factors and variables. Explanatory factors and variables used to model the species distributions of Márquez et al 2011 [19].


Click here for additional data file (pone.0053646.s002.doc)

We are grateful to P. Acevedo for his comments on earlier drafts of the manuscript. We thank the Agencia Estatal de Meteorología of Spain for providing the climatic data.

1. Guisan A,, Zimmermann NE, (Year: 2000) Predictive habitat distribution models in ecology. Ecol Model135: 147–186
2. Huntley B,, Collingham YC,, Willis SG,, Green RE, (Year: 2008) Potential impacts of climatic change on European breeding birds. PLoS One3 (1) e1439 doi:10.1371/journal.pone.0001439. 18197250
3. Duncan RP,, Cassey P,, Blackburn TM, (Year: 2009) Do climate envelope models transfer? A manipulative test using dung beetle introductions. P Roy Soc B-Biol Sci276: 1449–1457
4. Beaumont LJ,, Hughes L,, Pitman AJ, (Year: 2008) Why is the choice of future climate scenarios for species distribution modelling important?Ecol Lett11: 1135–114618713269
5. Dormann CF,, Schweiger O,, Arens P,, Augenstein I,, Aviron ST,, et al. (Year: 2008) Prediction uncertainty of environmental change effects on temperate European biodiversity. Ecol Lett11: 235–24418070098
6. Baer P,, Risbey JS, (Year: 2009) Uncertainty and assessment of the issues posed by urgent climate change. An editorial comment. Climatic Change92: 31–36
7. Diniz-Filho JAF,, Bini LM,, Rangel TF,, Loyola RD,, Hof C,, et al. (Year: 2009) Partitioning and mapping uncertainties in ensembles of forecasts of species turnover under climate change. Ecography32: 897–906
8. Real R,, Márquez AL,, Olivero J,, Estrada A, (Year: 2010) Species distribution models in climate change scenarios are not useful yet for informing emission policy planning: an uncertainty assessment using fuzzy logic. Ecography33: 304–314
9. Engler R,, Randin CF,, Thuiller W,, Dullinger S,, Zimmermannk NE,, et al. (Year: 2011) 21st century climate change threatens mountain flora unequally across Europe. Global Change Biol17: 2330–2341
10. Pearson RG,, Dawson TP, (Year: 2003) Predicting the impacts of climate change on the distribution of species: are bioclimate envelope models useful?Global Ecol Biogeogr12: 361–371
11. Araujo MB,, Thuiller W,, Pearson RG, (Year: 2006) Climate warming and the decline of amphibians and reptiles in Europe. J Biogeogr33: 1712–1728
12. Thuiller W,, Broennimann O,, Hughes G,, Alkemade JRM,, Midgley GF,, et al. (Year: 2006) Vulnerability of African mammals to anthropogenic climate change under conservative land transformation assumptions. Global Change Biol12: 424–440
13. Green RE,, Collingham YC,, Willis SG,, Gregory RD,, Smith KW,, et al. (Year: 2008) Performance of climate envelope models in retrodicting recent changes in bird population size from observed climatic change. Biol Lett4: 599–60218664422
14. Stefanescu C,, Herrando S,, Páramo F, (Year: 2004) Butterfly species richness in the north-west Mediterranean Basin: the role of natural and human-induced factors. J Biogeogr31: 905–916
15. Ritchie EG,, Martin JK,, Johnson CN,, Foz BJ, (Year: 2009) Separating the influences of environment and species interactions on patterns of distribution and abundance: competition between large herbivores. J Anim Ecol78: 724–73119175695
16. WallisDeVries MF,, Baxter W,, Van Vliet AJH, (Year: 2011) Beyond climate envelopes: effect of weather on regional population trends in butterflies. Oecologia167: 559–57121590332
17. Real R,, Márquez AL,, Estrada A,, Muñoz AR,, Vargas JM, (Year: 2008) Modelling chorotypes of invasive vertebrates in Mainland Spain. Divers Distrib14: 364–373
18. Aragón P,, Lobo JM,, Olalla-Tárrega MA,, Rodríguez MA, (Year: 2010) The contribution of contemporary climate to ectothermic and endothermic vertebrate distributions in glacial refuge. Global Ecol Biogeogr19: 40–49
19. Márquez AL,, Real R,, Olivero J,, Estrada A, (Year: 2011) Combining climate with other influential factors for modelling climate change impact on species distribution. Climatic Change108: 135–157
20. Acevedo P,, Jiménez-Valverde A,, Melo-Ferreira J,, Real R,, Alves PC, (Year: 2012) Parapatric species and the implications for climate change studies: a case study on hares in Europe. Global Change Biol18: 1509–1519
21. Willis SG,, Thomas CD,, Hill JK,, Collingham YC,, Telfer MG,, et al. (Year: 2009) Dynamic distribution modelling: predicting the present from the past. Ecography32: 5–12
22. Lavergne S,, Thuiller W,, Molina J,, Debussche M, (Year: 2005) Environmental and human factors influencing rare plant local occurrence, extinction and persistence: a 115-year study in the Mediterranean region. J Biogeogr32: 799–811
23. Dubey S,, Shine R, (Year: 2011) Predicting the effects of climate change on reproductive fitness of an endangered montane lizard, Eulamprus leuraensis (Scincidae). Climatic Change107: 531–547
24. Jiménez-Valverde A,, Lobo JM, (Year: 2007) Determinants of local spider (Araneidae and Thomisidae) species richness on a regional scale: climate and altitude vs. habitat structure. Ecol Entomol32: 113–122
25. Hortal J,, Rodríguez J,, Nieto-Díaz M,, Lobo JM, (Year: 2008) Regional and environmental effects on the species richness of mammal assemblages. J Biogeogr35: 1202–1214
26. Randin CF,, Jaccard H,, Vittoz P,, Yoccoz NG,, Guisan A, (Year: 2009) Land use improves spatial predictions of mountain plant abundance but not presence-absence. J Veg Sci20: 996–1008
27. Peres-Neto PR,, Legendre P, (Year: 2010) Estimating and controlling for spatial structure in the study of ecological communities. Global Ecol Biogeogr19: 174–184
28. Godinho C,, Rabaça JE,, Segurado P, (Year: 2010) Breeding bird assemblages in riparian galleries of the Guadiana River basin (Portugal): the effect of spatial structure and habitat variables. Ecol Res25: 283–294
29. Barbosa AM,, Real R,, Márquez AL,, Rendón MA, (Year: 2001) Spatial, environmental and human influences on the distribution of otter (Lutra lutra) in the Spanish provinces. Divers Distrib7: 137–144
30. Galantinho A,, Mira A, (Year: 2009) The influence of human, livestock, and ecological features on the occurrence of genet (Genetta genetta): a case study on Mediterranean farmland. Ecol Res24: 671–685
31. Acevedo P,, Real R, (Year: 2011) Biogeographical differences between the two Capra pyrenaica subspecies, C. p. victoriae and C. p. hispanica, inhabiting the Iberian Peninsula: Implications for conservation. Ecol Model222: 814–823
32. Beniston M, (Year: 2003) Climatic change in mountain regions: a review of possible impacts. Climatic Change59: 5–31
33. Klanderud K,, Birks HJB, (Year: 2003) Recent increases in species richness and shifts in altitudinal distributions of Norwegian mountain plants. Holocene13: 1–6
34. Peñuelas J,, Boada MA, (Year: 2003) Global change-induced biome shift in the Montseny mountains (NE Spain). Global Change Biol9: 131–140
35. Wilson RJ,, Gutiérrez D,, Gutiérrez J,, Martínez D,, Agudo R,, et al. (Year: 2005) Changes to the elevational limits and extent of species ranges associated with climate change. Ecol Lett8: 1138–114621352437
36. Pauli H,, Gottfried M,, Reiter K,, Klettner C,, Grabherr G, (Year: 2007) Signals of range expansions and contractions of vascular plants in the high Alps: observations (1994–2004) at the GLORIA* master site Schrankogel, Tyrol, Austria. Global Change Biol13: 147–156
37. Wilson RJ,, Gutiérrez D,, Gutiérrez J,, Monserrat VJ, (Year: 2007) An elevational shift in butterfly species richness and composition accompanying recent climate change. Global Change Biol13: 1873–1887
38. Gasner MR,, Jankowski JE,, Ciecka AL,, Kyle KO,, Rabenold KN, (Year: 2010) Projecting the local impacts of climate change on a Central American montane avian community. Biol Conser143: 1250–1258
39. Muñoz-Díaz D,, Rodrigo F, (Year: 2004) Spatio-temporal patterns of seasonal rainfall in Spain (1912–2000) using cluster and principal component analysis: comparison. Ann. Geophys22: 1435–1448
40. García-Ortega E,, Fita L,, Romero R,, López L,, Ramis C,, et al. (Year: 2007) Numerical simulation and sensitivity study of a severe hail-storm in northeast Spain. Atmos Res83: 225–241
41. Real R,, Barbosa AM,, Vargas JM, (Year: 2006) Obtaining environmental favourability functions from logistic regression. Environ Ecol Stat13: 237–245
42. Nielsen C,, Hartvig P,, Kollmann J, (Year: 2008) Predicting the distribution of the invasive alien Heracleum mantegazzianum at two different spatial scales. Diver Distrib14: 307–317
43. Barbosa AM,, Real R,, Vargas JM, (Year: 2010) Use of coarse-resolution models of species' distributions to guide local conservation inferences. Conserv Biol24: 1378–138720455912
44. Borcard D,, Legendre P,, Drapeau P, (Year: 1992) Partialling out the spatial component of ecological variation. Ecology87: 2614–262517089669
45. Muñoz AR,, Real R,, Barbosa AM,, Vargas JM, (Year: 2005) Modelling the distribution of Bonelli's eagle in Spain: implications for conservation planning. Divers Distrib11: 477–486
46. Peres-Neto PR,, Legendre P,, Dray S,, Borcard D, (Year: 2006) Variation partitioning of species data matrices: estimation and comparison of fractions. Ecology87: 2614–262517089669
47. Heikkinen R,, Luoto M,, Araújo MB,, Virkkala R,, Thuiller W,, et al. (Year: 2006) Methods and uncertainties in bioclimatic envelope modelling under climate change. Prog Phys Geog30: 751–777
48. Triviño M,, Thuiller W,, Cabeza M,, Hickler T,, Araújo MB, (Year: 2011) The contribution of vegetation and landscape configuration for predicting environmental change impacts on Iberian birds. PLoS One6: e29373 doi:10.1371/journal.pone.0029373. 22216263
49. Pereira MC,, Itami RM, (Year: 1991) GIS-based habitat modeling using logistic multiple regression: A study of the Mt. Graham red squirrel. Photogramm Eng Rem S57: 1475–1486
50. Márquez AL,, Real R,, Vargas JM, (Year: 2004) Dependence of broad-scale geographical variation in fleshy-fruited plant species richness on disperser bird species richness. Global Ecol Biogeogr13: 295–304
51. Stanton JC,, Pearson RG,, Horning N,, Ersts P,, Akçakaya HR, (Year: 2012) Combining static and dynamic variables in species distribution models under climate change. Methods Ecol Evol3: 349–357
52. Kharouba HM,, Algar AC,, Kerr JT, (Year: 2009) Historically calibrated predictions of butterfly species' range shift using global change as a pseudo-experiment. Ecology90: 2213–222219739383
53. Buisson L,, Thuiller W,, Casajus N,, Lek S,, Grenouillet G, (Year: 2010) Uncertainty in ensemble forecasting of species distribution. Global Change Biol16: 1145–1157
54. Lawler JJ,, Shafer SL,, White D,, Kareiva P,, Maurer ED,, et al. (Year: 2009) Projected climate impacts for the amphibians of the Western Hemisphere. Conser Biol24: 38–50
55. Santos X,, Brito JC,, Sillero N,, Pleguezuelos JM,, Llorente GA,, et al. (Year: 2006) Inferring habitat-suitability areas with ecological modelling techniques and GIS: A contribution to assess the conservation status of Vipera latastei. Biol Conser130: 416–425
56. Rubidge EM,, Monahan WB,, Parra JL,, Camero SE,, Brashares JS, (Year: 2010) The role of climate, habitat, and species co-occurrence as drivers of change in small mammal distributions over the past century. Global Change Biol17: 696–708


[Figure ID: pone-0053646-g001]
doi: 10.1371/journal.pone.0053646.g001.
Figure 1  Forecasted favourability.

Distribution of the future favourability forecasted according to the apparent climatic effect (FAC) and that forecasted according to the pure climatic effect (FPC). These maps represent three situations: FAC < FPC, FAC FPC, FAC > FPC. The examples shown are those where the difference or similarity were most evident.

[Figure ID: pone-0053646-g002]
doi: 10.1371/journal.pone.0053646.g002.
Figure 2  Differences between forecasted local favourabilities.

Distribution of the uncertainty associated with differences between the favourabilities forecasted according to the apparent and the pure climatic effect (|FfClim-FfPClim|) for the three species and situations represented in Figure 1.

[TableWrap ID: pone-0053646-t001] doi: 10.1371/journal.pone.0053646.t001.
Table 1  Variation partitioning of the combined favourability models for the period 1961–1990.
AOGCM-SRES A. dickhilleni V. latasti A. fasciata C. pyrenaica
CGCM2-A2 PNCE 36.20 56.13 67.07 78.63
PCE 6.80 24.19 66.73 38.40
SCE 57.00 19.68 −33.80 −17.03
ρ 0.107 0.551 2.026 1.797
CGCM2-B2 PNCE 37.27 59.27 57.78 77.99
PCE 7.02 37.62 72.82 37.82
SCE 55.71 3.10 −30.60 −15.80
ρ 0.112 0.924 1.725 1.718
ECHAM4-A2/B2 PNCE 66.43 43.61 46.92 80.08
PCE 12.11 19.45 79.28 6.40
SCE 21.46 36.94 −26.19 13.53
ρ 0.361 0.345 1.493 0.321

Values shown are the percentages of variation explained by the Pure Non-Climatic Effect (PNCE), the Pure Climatic Effect (PCE) and the Shared Climatic Effect (SCE). ρ: Proportion of pure climatic factor in relation to the whole climatic factor.

[TableWrap ID: pone-0053646-t002] doi: 10.1371/journal.pone.0053646.t002.
Table 2  Differences between the favourability forecasted according to the apparent climatic effect (FAC) and to the pure climatic effect (FPC) for each species and period of time.
A. dickhilleni V. latasti A. fasciata C. pyrenaica
2011–2040 CGCM2-A2 21.874 0.043 −501.387 −0.518 −339.467 −0.151 −429.507 −0.209
CGCM2-B2 68.640 0.122 −30.569 −0.017 −353.400 −0.146 −413.958 −0.204
ECHAM4-A2 241.694 0.297 −385.007 −0.255 −215.186 −0.092 192.548 0.092
ECHAM4-B2 282.001 0.303 −279.382 −0.167 −336.028 −0.125 216.448 0.121
2041–2070 CGCM2-A2 −150.240 −0.489 −416.164 −0.339 −721.372 −0.277 −753.816 −0.294
CGCM2-B2 −28.377 −0.063 −62.350 −0.048 −517.740 −0.195 −531.307 −0.230
ECHAM4-A2 216.533 0.270 −406.433 −0.277 −656.160 −0.158 406.952 0.197
ECHAM4-B2 15.849 0.034 50.292 0.023 −724.439 −0.202 559.879 0.246
2071–2100 CGCM2-A2 −335.978 −4.230 −328.051 −0.225 −864.559 −0.309 724.801 0.176
CGCM2-B2 −269.066 −1.262 512.495 0.173 −613.793 −0.217 −551.656 −0.237
ECHAM4-A2 −9.056 −0.022 350.620 0.133 −285.609 −0.114 943.755 0.335
ECHAM4-B2 148.574 0.227 214.544 0.089 −578.274 −0.177 813.794 0.311
Mean absolute percentage of change 61.33 18.86 18.05 22.10

R: Relative proportion of change ((FAC-FPC) / FAC).

Article Categories:
  • Research Article
Article Categories:
  • Biology
    • Computational Biology
      • Ecosystem Modeling
    • Ecology
      • Ecosystems
        • Ecosystem Modeling
      • Biodiversity
      • Biogeography
      • Conservation Science
      • Ecological Environments
      • Global Change Ecology
      • Macroecology
      • Spatial and Landscape Ecology
    • Zoology
Article Categories:
  • Computer Science
    • Geoinformatics
      • Environmental Systems Modeling
Article Categories:
  • Earth Sciences
    • Atmospheric Science
      • Climatology
        • Climate Change
        • Climate Modeling
    • Geography
      • Geoinformatics
        • Environmental Systems Modeling
Article Categories:
  • Social and Behavioral Sciences
    • Geography
      • Geoinformatics
        • Environmental Systems Modeling

Previous Document:  Analysis of V2 antibody responses induced in vaccinees in the ALVAC/AIDSVAX HIV-1 vaccine efficacy t...
Next Document:  The different role of Notch1 and Notch2 in astrocytic gliomas.