Inconsistent Range Shifts within Species Highlight Idiosyncratic Responses to Climate Warming

Climate in part determines species’ distributions, and species’ distributions are shifting in response to climate change. Strong correlations between the magnitude of temperature changes and the extent of range shifts point to warming temperatures as the single most influential factor causing shifts in species’ distributions species. However, other abiotic and biotic factors may alter or even reverse these patterns. The importance of temperature relative to these other factors can be evaluated by examining range shifts of the same species in different geographic areas. When the same species experience warming in different geographic areas, the extent to which they show range shifts that are similar in direction and magnitude is a measure of temperature’s importance. We analyzed published studies to identify species that have documented range shifts in separate areas. For 273 species of plants, birds, mammals, and marine invertebrates with range shifts measured in multiple geographic areas, 42-50% show inconsistency in the direction of their range shifts, despite experiencing similar warming trends. Inconsistency of within-species range shifts highlights how biotic interactions and local, non-thermal abiotic conditions may often supersede the direct physiological effects of temperature. Assemblages show consistent responses to climate change, but this predictability does not appear to extend to species considered individually.


Introduction
Climate change has altered ecological phenomena across the globe [1], [2].Meta-analyses of range shifts from climate change show a strong, generalized pattern of climate tracking with many species shifting their ranges uphill or poleward to stay within suitable climatic conditions [3]- [6].Demographically, synchronous population fluctuations in response to climatic fluctuations (e.g., the Moran effect; [7]) have been documented as a consequence of rapid climate change [8], [9].
Temperature appears to be a dominant factor causing range shifts, although many factors influence species' distributions [10].A critical aspect of the most recent meta-analysis of terrestrial range shifts emphasized that warming was sufficient to predict the extent of range shifts: mean distances shifted and temperature increases were significantly related (r 2 = 0.59 and r 2 = 0.37 for latitudinal and elevational range shifts, respectively [4]).Thus, changes in species' populations and ranges appear to be largely predictable on the basis of temperature when many species are considered together.
Prior analyses have focused on analyzing range shifts across species within the same area [4].In essence, such approaches examine the consistency of range shifts while holding the environmental context constant and varying the species.A different approach is to examine the consistency of range shifts within species.If the species are held constant, but the environmental setting varies, how often will the same species show the same response?The key question in this analysis is whether species' responses in one area are predictive of species' responses in another area.Range shifts are caused by a variety of factors, but temperature, precipitation, species' interactions, and dispersal are among the most important [11].Temperature data are often available in conjunction with studies of range shifts, whereas data for precipitation, species' interactions, and dispersal are less common.Here, we use paired populations of the same species, in different environmental contexts, where temperature is increasing in both situations.Few ecologists would suggest that temperature explains either all or none of the variation in species' range shifts.However, how often warming temperatures lead to uphill or poleward range shifts within the same species has not yet been evaluated, and doing so would help ecologists understand the relative influence of temperature changes against a background of other ecological factors.If responses are consistent within species, then species moving uphill/poleward in one area should also be moving uphill/poleward in other areas.Similarly, species moving downhill/toward the equator in one area should also be moving downhill/toward the equator in different areas.As a result, few species will be moving in opposite directions between areas (Fig 1).Alternatively, factors such as species' interactions or dispersal barriers may restrict species from tracking shifting climatic conditions.
An increasing number of studies document changes in species' ranges as a consequence of climate change, providing the basis for meta-analyses that examine the consistency of responses to climate change across species [4], [5].As the number of studies increases, some species have been measured in different geographic areas, providing the opportunity to examine responses to climate change within the same species.We identified species with documented responses to temperature increases in distinct geographic areas to evaluate the consistency of range shifts within species.We expected that the within-species consistency would meet or exceed the among-species consistency in range shifts associated with warming.

Selection of studies
We conducted a literature search to identify studies documenting climate-induced range shifts.We based our search on papers cited within two previous meta-analyses [3], [4] or any papers that have subsequently cited those studies, using Web of Science, in January 2014.Studies were excluded if they did not include information on individual species range shifts or if they did not include information on temperature trends within the study area; however, we included some studies that did not report temperature data in the original publication if temperature data for the area were included in a recent meta-analysis [4].We also excluded species if their range shifts in different areas were assessed with different metrics (e.g., the upper elevational limit was measured in one area, and the optimal elevational limit was measured in another area).All paired comparisons therefore either changes at the upper limits or changes at a measure of the range centroid (e.g., mean or optimal elevation).

Quantification of range shifts
The method of calculating range shifts in this analysis matches the methods used in a previous meta-analysis [4].The range shift was measured as the difference between the recent range limit (or range centroid) and the historic range limit (or range centroid), with uphill or poleward shifts receiving positive values and equatorial or downhill shifts receiving negative values [4].All studies included in the analysis experienced warming trends, eliminating the possibility that inconsistent range shifts could be caused by species shifting in response to different thermal trends.Both the magnitude of range shifts and the temperature change during the study period were based on the same values as a prior analysis, and the metric of distance shifted per years, standardized as m Ã year -1 for elevational shifts and km Ã year -1 for latitudinal shifts [4].Thus, the values for range shifts and temperature change were standardized to reveal acrossspecies predictability on the basis of temperature.We also standardized the range shifts by the rate of temperature change during the study period to account for different rates of warming in different areas.The resulting metric was distance shifted Ã °C-1 , which incorporates the rate of range shift and the rate of temperature change:

Consistency of the range shift
We evaluated the consistency of range shifts within species in two ways (full details on how we extracted data from each study are available in S1 File).Because the estimates of the rate of temperature change and the rate of range shifts are both subject to uncertainty, we first analyzed only the consistency of the direction in which species moved.The direction species shift is based only on the sign of the distance shifted (e.g., positive for uphill/poleward, negative for downhill/toward the equator), and this will be the same whether analyzed as distance Ã year -1 or distance Ã °C-1 .There were three possible categories: species could move in the direction predicted by warming temperatures (uphill or poleward), opposite the direction predicted by warming (downhill or toward the equator), or show no change in their distribution (given as exactly 0 in published studies).Species with populations showing responses that fell into more than one of the above categories across geographic areas were classified as inconsistent.To provide a measure of inconsistency that used only shifts in opposite directions (e.g., uphill vs. downhill), we also analyzed a subset of the data in which all species with range shift of exactly 0 were removed.
We also examined the consistency of species' range shifts (measured as distance Ã °C-1 ) by examining how well the magnitude of range shifts in one area predicted the magnitude of range shifts for the same species in another area.We did this by measuring the amount of variation in one area that could be explained by responses in another area, for which we used the coefficient of determination.It is possible that there are two classes of species: one that shows idiosyncratic responses to climate change and another that shows consistent responses to climate change.As another way of examining the predictive value of range shifts, we also analyzed the data when species that shifted in different directions were excluded.This method examined only species that shifted in the same direction to see if these species were also more consistent in the magnitude of their range shifts (measured as distance Ã °C-1 ) in different areas.We used the ratio of the larger range shift to the smaller range shift to estimate the difference in magnitude of range shifts, with a value of 1 representing identical range shifts in both areas.Some species were measured in three different geographic areas.For these species, we created three paired comparisons to measure consistency.To compensate for the greater number of comparisons possible when a response was measured in three areas, we weighted the results for each study by the total number of species.For this reason, the total number of species examined is sometimes less than the number of paired comparisons.
We compared observed patterns of paired range shifts to expected values under a null model where species would shift independently in the two areas (i.e., idiosyncratically).There are four different ways a species can shift its range in paired areas: uphill in both areas, downhill in both areas, and two ways in which it can shift uphill in one area and downhill in another (Fig 1).The null model calculates the expected number of species shifting in these four ways based only on the overall proportion of species shifting uphill or downhill in each study (e.g., calculations for Areas 1 and 2 in Fig 2).Under this scenario, species identity has no influence on how species respond, and the direction any given species shifts is unrelated to how it has shifted elsewhere.Alternatively, if range shifts accompanying warming are consistent within species, there should be more species shifting uphill in both areas or downhill in both areas but fewer species shifting uphill in one area and downhill in another.We used a χ 2 test in R (version 3.0.2[12]) to determine if the observed number of species in each of the four quadrats in .We omitted all paired comparisons in which there was a species range shift of exactly 0, as reported in the original study.We also omitted pairs of studies in which one study had all species range shifts in only one direction and the other study had species range shifts in both directions, as this produces a scenario in which the observed and expected values must be equal.We restricted the χ 2 analysis to paired comparisons in which the expected values in each cell were five or higher in most cells, as this is a common practice for χ 2 tests.This resulted in four paired comparisons.In two instances, the minimum expected values were four in one of the cells but we retained these analyses because the minimum values met or exceeded five in all other cells.

Direction of Range Shifts
We identified 13 studies [11], [13][14][15][16][17][18][19][20][21][22][23][24] encompassing 273 species of plants, birds, mammals, and marine invertebrates across North America, Europe, and Australia with documented range shifts in two or three geographic areas (Table 1; data for all studies analyzed are presented in S1 File).Across all the species included in the analysis, 61% (426 out of 700 populations) shifted in the direction expected from warming, closely matching an estimate that 65% of species shift uphill with warming [25].However, when we analyzed the species in paired comparisons, 50% shifted their distributions in categorically different directions in different geographic areas.The proportion of species shifting inconsistently was similar across plants, birds, marine invertebrates, and mammals (47%, 54%, 46%, and 60%, respectively).Omitting comparisons in which a species showed no range shift in at least one area, 42% of species moved in opposite directions (n = 257 paired comparisons).The observed number of species shifting in each of the four categories (uphill/uphill, downhill/downhill, uphill/downhill, or downhill/uphill) was not significantly different from the expected numbers under the null model in three of the four χ 2 analyses (see S2 File).
We evaluated whether the proportion of species shifting inconsistently depended on how the range shifts were evaluated.As noted previously (Table 1), 50% of species displayed inconsistent responses when range shifts were based on changes in the cold-edge limits (upper elevation or poleward latitude).However, the results were nearly identical when range shifts were based on changes in the species' range centroids (optimal or mean elevations), with 49% displaying inconsistent responses.

Magnitude of range shifts
When all 273 species were considered together, the weighted mean coefficient of determination was minimal (r 2 = 0.05), as illustrated for birds in Europe (Fig 3) and birds in the southwestern U.S. (Fig 4; data for all comparisons are available in S1-S14 Figs).The results were not appreciably different when analyzed in terms of distance shifted per year, and we report the results as distance per degree of warming because we believe this better standardizes across areas with different levels of warming.Some species were consistent in the direction they shifted when measured in paired areas.For this subset of species, which included only species with non-zero range shifts, the mean ratio of the larger range shift divided by the smaller range shift was 8.8 (SD = 17.5).

Discussion
The main result of this study is somewhat counterintuitive: climate-induced range shifts are generally predictable in their direction and magnitude, but this predictability declines substantially when the unit of analysis is an individual species, not a large group of species.This is a surprising result because a coherent pattern emerged from individual components that behaved idiosyncratically.The proportion of species showing inconsistent range shifts was fairly similar across a broad range of taxonomic groups (plants, birds, mammals, and marine invertebrates) and locations (North America, Europe, and Australia).Therefore, inconsistent range shifts seem to be a widespread response to climate change rather than a phenomenon in a single area or taxonomic group.
While previous research has highlighted the variability of species' range shifts (e.g., [4], [25]), quantifying the within-species variability had not been accomplished at such a broad geographic scale.As such, the results do not have an empirical benchmark to which they are compared.On one hand, within-species consistency of 50-58% could be viewed as high.Range shifts are governed by many factors, and for a single factor to explain so much of how species shift is a testament to the importance of temperature.On the other hand, 50-58% within-species consistency also suggests that species' responses to climate change will be modulated by other factors that do not necessarily co-vary with changes in temperature.We were unable to Table 1.Studies used to compare the direction and magnitude of paired range shifts.
The percentage of species with inconsistent range shifts is reported including all species, and the percentage of inconsistent range shifts when species with range shifts of zero are excluded is reported in italics; numbers of species in the analysis are reported in parentheses.
Temperature trends are reported as rates, calculated by the total warming divided by the study duration or as the rate reported in the original study.Precipitation data, where available, were calculated in the same way.For each paired comparison, we list the coefficient of determination to describe the amount of variation in range shifts in one area that can be explained by range shifts in the other area.For both comparisons involving bird species, there are three paired comparisons summarized in each row; we report the mean value, weighted by the number of species.The overall weighted mean 1 is listed as a footnote.Weighted mean across all studies is 50% (42% when species with range shifts of exactly 0 are excluded).doi:10.1371/journal.pone.0132103.t001 include any information from the world's largest mountain ranges, in South American and Asia, because no data were available; therefore, studies that document paired range shifts in these areas would help to evaluate the cosmopolitan nature of the patterns we observed.Consistent species-level responses (scenario A in Fig 1 ) would suggest information about species' thermal physiology or other traits could sufficient for estimating species' responses in areas without data on species' range shifts.For instance, range shifts in an endangered species in one area could be used to estimate its response throughout its range, even where data on its distribution was scarce.We used χ 2 analyses to test this idea.In three out of four χ 2 analyses, the observed within-species consistency was more supportive of the null expectation that species respond independently between areas (scenario B in Fig 1).We also asked if variation in range shifts in one area can explain range shifts in another area by examining the coefficient of determination of range shifts in paired comparisons.Our results indicate that knowing a Some species were consistent in the direction they shifted when measured in paired areas.For this subset of species, which included only species with non-zero range shifts, the mean ratio of the larger range shift divided by the smaller range shift was 8.8.Even considering uncertainties in temperature measurements and range shift measurements, the observed ratio of 8.8 indicates there is almost an order of magnitude difference in the rate at which populations of a species respond to similar climate change.Further, this disparity exists for the subset of species that are moving in consistent directions, and up to half of all species are not even shifting in consistent directions.

Ecological mechanisms to explain range shifts
Explaining the idiosyncratic nature of range shifts is an urgent question that we cannot fully assess using the data available; however, there are plausible explanations for the responses we observed.Two mechanisms, which are not mutually exclusive, may explain the inconsistency in range shifts.First, biotic interactions may exert a greater influence than climatic factors in determining distribution limits for some species.It is generally assumed that biotic interactions shape the warm-edge limits of species' distributions and physiological responses to climate shape the cold-edge limits [26] (i.e., the poleward or uphill range limits).However, this may not result in uniform poleward or uphill shifts of species as temperatures increase.Theory suggests species with warm-edge limits (i.e., their equatorial or downhill range limits) set by competition may expand downhill when climate change acts a disturbance that decreases the strength of competitive interactions [25].In addition, experimental manipulations have demonstrated that species interactions can override predictions of species' responses to climate change.Species' responses to experimental warming plots in a California grassland initially proceeded as expected based on physiological responses to temperature and water availability, but species interactions reversed the outcomes after several years [27].Similarly, the effects of competition and predation in mesocosm experiments led to limits on species ranges that differed substantially from those expected based only on physiological responses [28].
As species track climate change to new regions, they may encounter species with which they had not previously interacted.In such novel assemblages, or groups of co-occurring species with no modern analog [29], [30], species' distributions may be heavily influenced by strong biotic interactions.In the late Quaternary, North American mammals shifted in different directions during a period of climate change [31], similar to modern patterns documented here.In both cases, interactions between species whose ranges did not previously overlap may explain part of the inconsistency in species' distributions.As more rapid climate change occurs in the 21st century, the occurrence of novel assemblages is likely to be greater, as well as the impact of novel interactions on species' distributions.
A second mechanism that could explain inconsistent range shifts involves counter-gradients shifts such that species' optimum climatic conditions move downhill during warming [11].This type of response may be particularly likely to happen in ecosystems where water availability is more limiting than temperature.In such ecosystems, rising temperatures may be accompanied by increased precipitation, resulting in water availability isoclines moving downhill with warming.As a result, species are "pushed" uphill by rising temperatures but "pulled" downhill by greater water availability, as demonstrated in birds in California [11].Precipitation may be the single most important factor in arid or Mediterranean ecosystems [11].Therefore, species in arid ecosystems that display range shifts in opposite directions in different areas may be responding in a consistent fashion if the optimum climatic conditions are moving in different directions in different parts of their geographic distribution.
However, precipitation is unlikely to explain all of the variation observed in the present study.For instance, inconsistent range shifts occurred in marine ecosystems, where water availability is not an issue.Studies of marine [32] and freshwater [33] fishes have also shown counter-intuitive range shifts, suggesting range shifts opposite the direction expected by warming is a response that occurs even where water is not limiting.The lack of precipitation data in the studies analyzed here precludes any substantial analysis of this factor.Only in three instances were precipitation data available for direct comparison between studies; in all three instances, both temperature and precipitation increased in both areas, and these instances, the trend of species' response was similar to the overall mean (Table 1).Examining the relative influence of temperature and precipitation is an important part of understanding range shifts, but our analysis was restricted to temperature.
Barriers to dispersal may account for some of the instance in which a population showed no range shift.Documenting barriers to disperal in studies of climate-induced range shifts is difficult, but a study of elevational range shifts in tropical moths identified geology as a potential barrier that may have prevented species from shifting their ranges uphill [34].Moths use host plants, and their host plants may be unable to shift uphill if geological formations uphill of their current limits contain mineral formations that are unsuitable.

Examining potential confounding factors
One potential explanation for species shifting in different directions between studies is that highly variable responses to climate change occur over short periods, but more consistent responses will emerge over longer timescales [1].However, there was little correlation between the duration of studies, which ranged from a mean of 18 to 89 years between sampling intervals, and the proportion of species showing inconsistent range shifts (r 2 < 0.01; S15 Fig) .Thus, the temporal scale did not influence the results.A potential explanation for equatorial or downhill range shifts is that populations may be tracking local, counter-gradient temperature trends [6], so these directional shifts are actually consistent with the concept of climate tracking at local scales.However, the spatial scale in the analysis of within-species range shifts was the same as the spatial scale in the analysis of across-species range shifts.In both instances, the sources of data for range shifts and temperature change were the same, so the results cannot be explained by using different spatial scales for different analyses.It is possible that microclimatic refugia account for some of the species that did not shift their distributions in either direction (e.g., [35]), but further analysis of this was not possible with the data we analyzed.
In any analysis of published literature, there is concern that bias in the selection of studies can shape the trends observed.However, it is unlikely that selective publication or inclusion of studies is influencing the trends we measured.When range shifts are considered across species, the proportion of uphill or poleward range shifts closely matches previous estimates [4], [25].Therefore, we are not analyzing a subset of species that show unusual responses but rather a large (n = 700 populations) group of species that are representative of climate-induced range shifts in general.Similarly, publication bias is unlikely to be influential in this analysis.We excluded studies that published single species accounts, as these are more likely to selectively focus on exceptional responses.Moreover, our analysis could only be done for species whose distributions had been studied in earlier time periods meaning the selection of which species are included was set by ecologists working several decades ago, before climate-induced range shifts were a consideration.

Predictable assemblages, unpredictable populations
We interpret the collective results as evidence that temperature is sufficient to predict mean range shifts across an assemblage of species but insufficient to explain the results at the level of individual species.Each method of analyzing range shifts within species produced similar results.At the broadest level, the direction species moved was frequently inconsistent.The magnitude of a population's range shift in one area was a poor predictor of the magnitude in another area, and even for the subset of species shifting in the same direction, there was about an order of magnitude difference in the response.The χ 2 test indicates the populations appear to be much less predictable than assemblages.This implies that rising temperatures act as a force pushing most species to track temperatures, although which species do so in an area is strongly influenced by non-thermal factors.
Importantly, the results we observed do not contradict previous meta-analyses: most populations (61%) in our study shifted in the direction expected from warming, which agrees well with estimates from other analyses [4], [25].An intriguing aspect of climate-induced range shifts is that a coherent pattern emerges when species are considered collectively even though species appear to behave less predictably when considered individually (e.g., scenario B in Fig 1).Although it is a fundamental premise that species' physiological requirements will drive their response to climate change [36], the low within-species consistency in response to rising temperatures suggests biotic interactions and local climatic conditions may often overwhelm the influence of temperature, at least at the warming rates observed in the 20th century.While we expected that non-thermal conditions would be somewhat influential in shaping range shifts, we were surprised at how influential they appear to be for the studies we examined.

Conservation implications
If species shift their distributions one way in one region and shift their distributions in another direction in a separate region, what is the effect on their rangewide distribution?Presently, there is little evidence to suggest changes in geographic distribution in a single area could be used to extrapolate to a species' entire range.Studies documenting changes in distribution throughout a species' entire range are rare (but see [37], [38]).Consequently, many of the inferences regarding climate change impacts on species' distributions come from transect studies, partly because mountain transects contain large temperature gradients over a short distance and partly because these are often the best baseline data from past decades.Yet, the results of this analysis suggest that measuring changes in distributions in a single area may provide little information about the total change in a species' distribution as a consequence of climate change.Much of the knowledge regarding climate-related range shifts comes from syntheses of many studies [4], [39], most of which do not document range shifts across a species' entire range.The present research suggests that the most general conclusions of range shift syntheses are correct (most species in any given area shift to cooler regions), but what happens along a transect may not be representative of what happens at a much larger spatial scale, such as an entire mountain range.
This study highlights a problem for managing and conserving biodiversity at the level of an individual species, as is often mandated through laws such as the Endangered Species Act of the United States.Forecasting will probably work well for assemblages, because the signal from warming becomes apparent when averaged across a large numbers of species.But forecasting may be less accurate for single species since local, non-thermal climatic factors and species interactions may override the effects expected due to regional warming.Similarly, projections of temperature trends for the 21 st century may be less variable than projections of precipitation trends.Locally-refined estimates of climate change and incorporation of biotic interactions will produce better predictive frameworks but will require more effort and cost to collect the necessary data [40], [41].

Fig 1 .
Fig 1. Different levels of within-species consistency can produce the same level of across-species consistency.(A) Consistent within-species range shifts from climate cahnge in two geographic areas versus idiosyncratic (B) within-species range shifts.For each scenario, the responses of 20 species (A through T) are depicted.In scenario A, 100% of the species show a consistent response between areas whereas in scenario B, only 50% of the species show a consistent response between areas.However, in both scenarios, 60% of the species show the expected uphill response to climate warming.doi:10.1371/journal.pone.0132103.g001 Fig 2 exceeded the null expectation (see Fig 2 for examples of how expected values were calculated)

Fig 2 .
Fig 2. Derivation of expected values for the χ 2 test.Each of four possible scenarios for paired range shifts are depicted.doi:10.1371/journal.pone.0132103.g002

Fig 3 .
Fig 3. Idiosyncratic elevational range shifts among 17 bird species in France and Italy.Data from France (blue bars) are from [13], and data from Italy (red bars) are from [17].Positive numbers indicate uphill shifts.Numbers on x axis refer to species codes given in the Supplementary Information (see S1 File for sources of data).doi:10.1371/journal.pone.0132103.g003

Fig 4 .
Fig 4. Elevational range shifts among 58 species of birds in California.Positive numbers indicate uphill shifts.Yosemite National Park is area A in the inset map, and Southern California is area B in the inset map.Both areas were reported in [11].Scatterplots for all paired studies are available in the Supporting Information Figures (S1-S14 Figs).doi:10.1371/journal.pone.0132103.g004