Beyond climate envelopes: effects of weather on regional population trends in butterflies

Butterfly abundance changes in England are well associated with extreme climate events

Climate change exerts significant impacts on ecosystems through extreme climate events (ECEs), which are linked to various climate variables and can occur in different seasons. However, previous studies predominantly focus on a single type of ECE within specific seasons. We address this research gap by examining four typical types of ECEs: droughts, pluvials, heatwaves, and coldwaves, which are consistently defined and identified on a daily scale using a recently proposed statistical method. Butterflies in England serve as our study subject since their life stages occur throughout the year, and we had access to a 45-year dataset encompassing 57 butterfly species. First, we analyzed the correlation between abundance changes and the severity of ECEs, revealing varying sensitivity to ECEs across different life stages. Notably, abundance changes are negatively correlated with the severity of heatwaves in the hibernation, egg, and larval stages, as well as with the severity of pluvials in the larval, pupal, and adult stages. Second, we identified the most extreme climate events (MECEs) in England in the period from 1950 to 2020 and linked them with synchronized abundance changes for all species. MECEs in the sensitive stages of butterflies are frequently (for eight out of 11 MECEs) associated with synchronized decreases in their abundance. Our findings underscore the importance of considering diverse types of ECEs across all seasons to gain insights into their potential ecosystem effects. We draw attention to the fact that our analyses are primarily data-driven, with limited ecophysiological interpretation.

Weather explains the decline and rise of insect biomass over 34 years

Insects have a pivotal role in ecosystem function, thus the decline of more than 75% in insect biomass in protected areas over recent decades in Central Europe1 and elsewhere2,3 has alarmed the public, pushed decision-makers4 and stimulated research on insect population trends. However, the drivers of this decline are still not well understood. Here, we reanalysed 27 years of insect biomass data from Hallmann et al.1, using sample-specific information on weather conditions during sampling and weather anomalies during the insect life cycle. This model explained variation in temporal decline in insect biomass, including an observed increase in biomass in recent years, solely on the basis of these weather variables. Our finding that terrestrial insect biomass is largely driven by complex weather conditions challenges previous assumptions that climate change is more critical in the tropics5,6 or that negative consequences in the temperate zone might only occur in the future7. Despite the recent observed increase in biomass, new combinations of unfavourable multi-annual weather conditions might be expected to further threaten insect populations under continuing climate change. Our findings also highlight the need for more climate change research on physiological mechanisms affected by annual weather conditions and anomalies.

Winter Rains Support Butterfly Diversity, but Summer Monsoon Rainfall Drives Post-Monsoon Butterfly Abundance in the Arid Southwest of the US

Butterfly populations are declining worldwide, reflecting our current global biodiversity crisis. Because butterflies are a popular and accurate indicator of insect populations, these declines reflect an even more widespread threat to insects and the food webs upon which they rely. As small ectotherms, insects have a narrow range of habitable conditions; hence, extreme fluctuations and shifts caused by climate change may increase insects' risk of extinction. We evaluated trends of butterfly richness and abundance and their relationship with relevant climate variables in Arizona, U.S.A., using the past 40 years of community science data. We focused on precipitation and temperature as they are known to be influential for insect survival, particularly in arid areas like southwestern U.S.A. We found that preceding winter precipitation is a driver of both spring and summer/fall butterfly richness and spring butterfly abundance. In contrast, summer/fall butterfly abundance was driven by summer monsoon precipitations. The statistically significant declines over the 40-year period were summer/fall butterfly abundance and spring butterfly richness. When controlling for the other variables in the model, there was an average annual 1.81% decline in summer/fall season butterfly abundance and an average annual decline of 2.13 species in the spring season. As climate change continues to negatively impact winter precipitation patterns in this arid region, we anticipate the loss of butterfly species in this region and must consider individual butterfly species trends and additional management and conservation needs.

Extinction risk modeling predicts range-wide differences of climate change impact on Karner blue butterfly (Lycaeides melissa samuelis)

The Karner blue butterfly (Lycaeides melissa samuelis, or Kbb), a federally endangered species under the U.S. Endangered Species Act in decline due to habitat loss, can be further threatened by climate change. Evaluating how climate shapes the population trend of the Kbb can help in the development of adaptive management plans. Current demographic models for the Kbb incorporate in either a density-dependent or density-independent manner. We instead created mixed density-dependent and -independent (hereafter "endo-exogenous") models for Kbbs based on long-term count data of five isolated populations in the upper Midwest, United States during two flight periods (May to June and July to August) to understand how the growth rates were related to previous population densities and abiotic environmental conditions, including various macro- and micro-climatic variables. Our endo-exogenous extinction risk models showed that both density-dependent and -independent components were vital drivers of the historical population trends. However, climate change impacts were not always detrimental to Kbbs. Despite the decrease of population growth rate with higher overwinter temperatures and spring precipitations in the first generation, the growth rate increased with higher summer temperatures and precipitations in the second generation. We concluded that finer spatiotemporally scaled models could be more rewarding in guiding the decision-making process of Kbb restoration under climate change.

Precipitation buffers temperature-driven local extinctions of moths at warm range margins

Species' distributions are moving polewards in response to climate change, and although range expansions of relatively warm-adapted species are widely reported, reports of range retractions in cool-adapted species are less common. Here, we analysed species' distribution shifts for 76 cool-adapted moths in Great Britain using citizen science occurrence records from the National Moth Recording Scheme over a 40-year period. Although we find evidence for trailing edge shifts to higher latitudes, shifts in species' range centroids are oriented towards the north-west, and are more closely correlated with directional changes in total precipitation than average temperature. We also found that species' local extinction risk is higher in areas where temperature is high and precipitation is low, but this risk diminishes as precipitation increases. Adaptation efforts should therefore focus on maintaining or increasing water availability as the climate continues to change.

Measuring and modelling microclimatic air temperature in a historically degraded tropical forest

Climate change is predicted to cause widespread disruptions to global biodiversity. Most climate models are at the macroscale, operating at a ~ 1 km resolution and predicting future temperatures at 1.5-2 m above ground level, making them unable to predict microclimates at the scale that many organisms experience temperature. We studied the effects of forest structure and vertical position on microclimatic air temperature within forest canopy in a historically degraded tropical forest in Sikundur, Northern Sumatra, Indonesia. We collected temperature measurements in fifteen plots over 20 months, alongside vegetation structure data from the same fifteen 25 × 25 m plots. We also performed airborne surveys using an unmanned aerial vehicle (UAV) to record canopy structure remotely, both over the plot locations and a wider area. We hypothesised that old-growth forest structure would moderate microclimatic air temperature. Our data showed that Sikundur is a thermally dynamic environment, with simultaneously recorded temperatures at different locations within the canopy varying by up to ~ 15 °C. Our models (R² = 0.90 to 0.95) showed that temperature differences between data loggers at different sites were largely determined by variation in recording height and the amount of solar radiation reaching the topmost part of the canopy, although strong interactions between these abiotic factors and canopy structure shaped microclimate air temperature variation. The impacts of forest degradation have smaller relative influence on models of microclimatic air temperature than abiotic factors, but the loss of canopy density increases temperature. This may render areas of degraded tropical forests unsuitable for some forest-dwelling species with the advent of future climate change.

Weather and butterfly responses: a framework for understanding population dynamics in terms of species' life-cycles and extreme climatic events

Understanding population responses to environmental conditions is key in the current context of climate change and the extreme climatic events that are threatening biodiversity in an unprecedented way. In this work, we provide a framework for understanding butterfly population responses to weather and extreme climatic seasons by taking into account topographic heterogeneity, species' life-cycles and density-dependent processes. We used a citizen-science database of Mediterranean butterflies that contains long-term population data (28 years) on 78 butterfly species from 146 sites in the Mediterranean mesic and alpine climate regions. Climatic data were obtained from 93 meteorological stations operating during this period near the butterfly sites. We studied how seasonal precipitation and temperature affect population growth while taking into account the effects of density dependence. Our results reveal (i) the beneficial effects of winter and spring precipitation for butterfly populations, which are most evident in the Mediterranean region and in univoltine species, and mainly affect the larval stage; (ii) a general negative effect of summer rain in the previous year, which affects the adult stage; and (iii) a consistent negative effect of mild autumns and winters on population growth. In addition, density dependence played a major role in the population dynamics of most species, except for those with long-term negative population trends. Our analyses also provide compelling evidence that both extreme population levels in previous years and extreme climatic seasons in the current year provoke population crashes and explosions, especially in the Mediterranean mesic region.

Modulating Effects of Landscape Characteristics on Responses to Warming Differ Among Butterfly Species

Understanding and predicting biodiversity responses to climate change are vital to inform conservation strategies, but this is not straightforward as climate change responses depend on the landscape context and differ among species. Here, we quantified changes in the distribution and abundance of 30 butterfly species in the Netherlands in relation to climate change and in landscapes that vary in the amount and connectivity of (semi-)natural vegetation (SNV). We obtained yearly counts of well-monitored butterfly species from 327 time series over 27 years (1992–2018). We used these counts to build mixed effect hurdle models to relate species’ occurrence and abundance to temperature and the amount and connectivity of SNV around the sites. For 55% of the butterfly species, an increased amount or connectivity of SNV corresponded with stronger increases or reduced decreases in occurrence in response to warming, indicating that SNV may facilitate range expansion or mitigate extirpations due to climate change. However, for the occurrence of the other species we found no or a negative interaction between warming and SNV. Further, we did not find indications of a mitigating effect of SNV on abundance responses to warming. Our results thus suggest that increasing the amount and connectivity of SNV does not offer a “one-size-fits-all” solution, highlighting the need for additional measures if butterfly diversity is to be conserved.

Exploring Cold Hardiness within a Butterfly Clade: Supercooling Ability and Polyol Profiles in European Satyrinae

The cold hardiness of overwintering stages affects the distribution of temperate and cold-zone insects. Studies on Erebia, a species-rich cold-zone butterfly genus, detected unexpected diversity of cold hardiness traits. We expanded our investigation to eight Satyrinae species of seven genera. We assessed Autumn and Winter supercooling points (SCPs) and concentrations of putatively cryoprotective sugars and polyols via gas chromatography–mass spectrometry. Aphantopus hyperantus and Hipparchia semele survived freezing of body fluids; Coenonympha arcania, C. gardetta, and Melanargia galathea died prior to freezing; Maniola jurtina, Chazara briseis, and Minois dryas displayed a mixed response. SCP varied from −22 to −9 °C among species. Total sugar and polyol concentrations (TSPC) varied sixfold (2 to 12 μg × mg−1) and eightfold including the Erebia spp. results. SCP and TSPC did not correlate. Alpine Erebia spp. contained high trehalose, threitol, and erythritol; C. briseis and C. gardetta contained high ribitol and trehalose; lowland species contained high saccharose, maltose, fructose, and sorbitol. SCP, TSPC, and glycerol concentrations were affected by phylogeny. Species of mountains or steppes tend to be freeze-avoidant, overwinter as young larvae, and contain high concentrations of trehalose, while those of mesic environments tend to be freeze-tolerant, overwinter as later instars, and rely on compounds such as maltose, saccharose, and fructose.

The influence of chalk grasslands on butterfly phenology and ecology

The influence of large-scale variables such as climate change on phenology has received a great deal of research attention. However, local environmental factors also play a key role in determining the timing of species life cycles. Using the meadow brown butterfly Maniola jurtina as an example, we investigate how a specific habitat type, lowland calcareous grassland, can affect the timing of flight dates. Although protracted flight periods have previously been reported in populations on chalk grassland sites in the south of England, no attempt has yet been made to quantify this at a national level, or to assess links with population genetics and drought tolerance. Using data from 539 sites across the UK, these differences in phenology are quantified, and M. jurtina phenology is found to be strongly associated with both site geology and topography, independent of levels of abundance. Further investigation into aspects of M. jurtina ecology at a subset of sites finds no genetic structuring or drought tolerance associated with these same site conditions.

The role of climate change in pollinator decline across the Northern Hemisphere is underestimated

Pollinator biodiversity loss occurs at unprecedented rates globally, with particularly sharp declines documented in the North Temperate Zone. There is currently no consensus on the main drivers of the decline. Although climate change is expected to drive biodiversity loss in the future, current warming is often suggested to have positive impacts on pollinator assemblages in higher latitudes. Consequently, pollinator conservation initiatives in Europe and the USA tend to lack climate adaptation initiatives, an omission of which may be risky if climate change has significant negative impacts on pollinators. To gain an understanding of the impacts of climate change on pollinator biodiversity in the Northern Hemisphere, we conducted a literature review on genetic, species and community level diversity. Our findings suggest that global heating most likely causes homogenization of pollinator assemblages at all levels of pollinator biodiversity, making them less resilient to future stochasticity. Aspects of biodiversity that are rarely measured (e.g. genetic diversity, β-diversity, species evenness) tend to be most affected, while some dimensions of climate change, such as fluctuations in winter weather conditions, changes in the length of the vegetational season and increased frequency of extreme weather events, that seldom receive attention in empirical studies, tend to be particularly detrimental to pollinators. Negative effects of global heating on pollinator biodiversity are most likely exacerbated by homogenous and fragmented landscapes, widespread across Europe and the US, which limit opportunities for range-shifts and reduce micro-climatic buffering. This suggests the need for conservation initiatives to focus on increasing landscape connectivity and heterogeneity at multiple spatial scales.

Insects and recent climate change

Insects have diversified through more than 450 million y of Earth's changeable climate, yet rapidly shifting patterns of temperature and precipitation now pose novel challenges as they combine with decades of other anthropogenic stressors including the conversion and degradation of land. Here, we consider how insects are responding to recent climate change while summarizing the literature on long-term monitoring of insect populations in the context of climatic fluctuations. Results to date suggest that climate change impacts on insects have the potential to be considerable, even when compared with changes in land use. The importance of climate is illustrated with a case study from the butterflies of Northern California, where we find that population declines have been severe in high-elevation areas removed from the most immediate effects of habitat loss. These results shed light on the complexity of montane-adapted insects responding to changing abiotic conditions. We also consider methodological issues that would improve syntheses of results across long-term insect datasets and highlight directions for future empirical work.

Predicted Future Benefits for an Endemic Rodent in the Irano-Turanian Region

Climate change is expected to have an impact on the geographical distribution ranges of species. Endemic species and those with a restricted geographic range may be especially vulnerable. The Persian jird (Meriones persicus) is an endemic rodent inhabiting the mountainous areas of the Irano-Turanian region, where future desertification may form a threat to the species. In this study, the species distribution modelling algorithm MaxEnt was used to assess the impact of future climate change on the geographic distribution range of the Persian jird. Predictions were made under two Representative Concentration Pathways and five different climate models for the years 2050 and 2070. It was found that both bioclimatic variables and land use variables were important in determining potential suitability of the region for the species to occur. In most cases, the future predictions showed an expansion of the geographic range of the Persian jird which indicates that the species is not under immediate threat. There are however uncertainties with regards to its current range. Predictions may therefore be an over or underestimation of the total suitable area. Further research is thus needed to confirm the current geographic range of the Persian jird to be able to improve assessments of the impact of future climate change.

Well-managed grassland heterogeneity promotes butterfly conservation in a corridor network

Habitat degradation is a major concern in transformed landscapes, as it reduces complexity by removing species, interactions, and ultimately biodiversity. Degradation is also of concern for ecological networks (ENs) composed of an interconnected system of conservation corridors among South Africa's commercial forestry compartments. These corridors are predominantly grasslands, and used as rangeland, so managed to optimize grazing conditions. Yet, how this management approach influences biodiversity remains unknown. Here, we studied how butterfly assemblages respond to local differences in rangeland quality (low, high and reference sites), and how this effect compared to that of local environmental variables (e.g. rockiness and bare ground), meso environmental gradients (e.g. topographic position and aspect), and landscape composition (i.e. proportion of different land cover types in the surrounding matrix). We calculated species richness and composition, Shannon's diversity index (H'), and the Butterfly Conservation Index (BCIn) representing the proportion of sensitive and range-restricted butterfly species per site. Rangeland quality was considered less important for butterflies than other environmental variables, but it was also significantly confounded with other environmental variables. At the landscape scale, proportion of grassland in the landscape matrix influenced butterfly assemblage composition, while proportion of thicket had a significant positive effect on BCIn. Moreover, the effect of elevation on assemblage composition emphasizes the value of maintaining environmental gradients within these conservation corridors. At the meso spatial scale, butterfly species richness and diversity (H') declined with increased dominance by a single plant species, which usually occurs late in a normal fire cycle. This suggests a reliance by butterflies on recurring natural disturbances for long-term persistence. We recommend moderate patch burning and grazing, as well as occasional hot burns to reduce thicket in Afromontane grassland. This approach would improve local scale vegetation patterns, and increase heterogeneity across the landscape for conserving these butterflies into the future.

Inadequacy of typical physiological experimental protocols for investigating consequences of stochastic weather events emerging from global warming

Increasingly variable, extreme, and nonpredictable weather events are predicted to accompany climate change, and such weather events will especially affect temperate, terrestrial environments. Yet, typical protocols in comparative physiology that examine environmental change typically employ simple step-wise changes in the experimental stressor of interest (e.g., temperature, water availability, oxygen, nutrition). Such protocols fall short of mimicking actual natural environments and may be inadequate for fully exploring the physiological effects of stochastic, extreme weather events. Indeed, numerous studies from the field of thermal biology, especially, indicate nonlinear and sometimes counterintuitive findings associated with variable and fluctuating (but rarely truly stochastic) protocols for temperature change. This Perspective article suggests that alternative experimental protocols should be employed that go beyond step-wise protocols and even beyond variable protocols employing circadian rhythms, for example, to those that actually embrace nonpredictable elements. Such protocols, though admittedly more difficult to implement, are more likely to reveal the capabilities (and, importantly, the limitations) of animals experiencing weather, as distinct from climate. While some possible protocols involving stochasticity are described as examples to stimulate additional thought on experimental design, the overall goal of this Perspective article is to encourage comparative physiologists to entertain incorporation of nonpredictable experimental conditions as they design future experimental protocols.

Sensitivity of UK butterflies to local climatic extremes: which life stages are most at risk?

There is growing recognition as to the importance of extreme climatic events (ECEs) in determining changes in species populations. In fact, it is often the extent of climate variability that determines a population's ability to persist at a given site. This study examined the impact of ECEs on the resident UK butterfly species (n = 41) over a 37-year period. The study investigated the sensitivity of butterflies to four extremes (drought, extreme precipitation, extreme heat and extreme cold), identified at the site level, across each species' life stages. Variations in the vulnerability of butterflies at the site level were also compared based on three life-history traits (voltinism, habitat requirement and range). This is the first study to examine the effects of ECEs at the site level across all life stages of a butterfly, identifying sensitive life stages and unravelling the role life-history traits play in species sensitivity to ECEs. Butterfly population changes were found to be primarily driven by temperature extremes. Extreme heat was detrimental during overwintering periods and beneficial during adult periods and extreme cold had opposite impacts on both of these life stages. Previously undocumented detrimental effects were identified for extreme precipitation during the pupal life stage for univoltine species. Generalists were found to have significantly more negative associations with ECEs than specialists. With future projections of warmer, wetter winters and more severe weather events, UK butterflies could come under severe pressure given the findings of this study.

Increasing frequency of low summer precipitation synchronizes dynamics and compromises metapopulation stability in the Glanville fritillary butterfly

Climate change is known to shift species' geographical ranges, phenologies and abundances, but less is known about other population dynamic consequences. Here, we analyse spatio-temporal dynamics of the Glanville fritillary butterfly (Melitaea cinxia) in a network of 4000 dry meadows during 21 years. The results demonstrate two strong, related patterns: the amplitude of year-to-year fluctuations in the size of the metapopulation as a whole has increased, though there is no long-term trend in average abundance; and there is a highly significant increase in the level of spatial synchrony in population dynamics. The increased synchrony cannot be explained by increasing within-year spatial correlation in precipitation, the key environmental driver of population change, or in per capita growth rate. On the other hand, the frequency of drought during a critical life-history stage (early larval instars) has increased over the years, which is sufficient to explain the increasing amplitude and the expanding spatial synchrony in metapopulation dynamics. Increased spatial synchrony has the general effect of reducing long-term metapopulation viability even if there is no change in average metapopulation size. This study demonstrates how temporal changes in weather conditions can lead to striking changes in spatio-temporal population dynamics.

Seasonal life history trade-offs in two leafwing butterflies: Delaying reproductive development increases life expectancy

Surviving inhospitable periods or seasons may greatly affect fitness. Evidence of this exists in the prevalence of dormant stages in the life cycles of most insects. Here I focused on butterflies with distinct seasonal morphological types (not a genetic polymorphism) in which one morphological type, or form, delays reproduction until favorable conditions return, while the other form develops in an environment that favors direct reproduction. For two butterflies, Anaea aidea and A. andria, I tested the hypothesis that the development of each seasonal form involves a differential allocation of resources to survival at eclosion. I assayed differences in adult longevity among summer and winter forms in either a warm, active environment or a cool, calm environment. Winter form adults lived 40 times longer than summer form but only in calm, cool conditions. The magnitude of this difference provided compelling evidence that the winter form body plan and metabolic strategy (i.e. resource conservatism) favor long term survival. This research suggests that winter form adults maintain lowered metabolic rate, a common feature of diapause, to conserve resources and delay senescence while overwintering.

No evidence of the effect of extreme weather events on annual occurrence of four groups of ectothermic species

Weather extremes may have strong effects on biodiversity, as known from theoretical and modelling studies. Predicted negative effects of increased weather variation are found only for a few species, mostly plants and birds in empirical studies. Therefore, we investigated correlations between weather variability and patterns in occupancy, local colonisations and local extinctions (metapopulation metrics) across four groups of ectotherms: Odonata, Orthoptera, Lepidoptera, and Reptilia. We analysed data of 134 species on a 1×1 km-grid base, collected in the last 20 years from the Netherlands, combining standardised data and opportunistic data. We applied dynamic site-occupancy models and used the results as input for analyses of (i) trends in distribution patterns, (ii) the effect of temperature on colonisation and persistence probability, and (iii) the effect of years with extreme weather on all the three metapopulation metrics. All groups, except butterflies, showed more positive than negative trends in metapopulation metrics. We did not find evidence that the probability of colonisation or persistence increases with temperature nor that extreme weather events are reflected in higher extinction risks. We could not prove that weather extremes have visible and consistent negative effects on ectothermic species in temperate northern hemisphere. These findings do not confirm the general prediction that increased weather variability imperils biodiversity. We conclude that weather extremes might not be ecologically relevant for the majority of species. Populations might be buffered against weather variation (e.g. by habitat heterogeneity), or other factors might be masking the effects (e.g. availability and quality of habitat). Consequently, we postulate that weather extremes have less, or different, impact in real world metapopulations than theory and models suggest.

Spatial climate patterns explain negligible variation in strength of compensatory density feedbacks in birds and mammals

The use of long-term population data to separate the demographic role of climate from density-modified demographic processes has become a major topic of ecological investigation over the last two decades. Although the ecological and evolutionary mechanisms that determine the strength of density feedbacks are now well understood, the degree to which climate gradients shape those processes across taxa and broad spatial scales remains unclear. Intuitively, harsh or highly variable environmental conditions should weaken compensatory density feedbacks because populations are hypothetically unable to achieve or maintain densities at which social and trophic interactions (e.g., competition, parasitism, predation, disease) might systematically reduce population growth. Here we investigate variation in the strength of compensatory density feedback, from long-term time series of abundance over 146 species of birds and mammals, in response to spatial gradients of broad-scale temperature precipitation variables covering 97 localities in 28 countries. We use information-theoretic metrics to rank phylogenetic generalized least-squares regression models that control for sample size (time-series length) and phylogenetic non-independence. Climatic factors explained < 1% of the remaining variation in density-feedback strength across species, with the highest non-control, model-averaged effect sizes related to extreme precipitation variables. We could not link our results directly to other published studies, because ecologists use contrasting responses, predictors and statistical approaches to correlate density feedback and climate--at the expense of comparability in a macroecological context. Censuses of multiple populations within a given species, and a priori knowledge of the spatial scales at which density feedbacks interact with climate, seem to be necessary to determine cross-taxa variation in this phenomenon. Despite the availability of robust modelling tools, the appropriate data have not yet been gathered for most species, meaning that we cannot yet make any robust generalisations about how demographic feedbacks interact with climate.

Estimating how inflated or obscured effects of climate affect forecasted species distribution

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.