Spatial Autocorrelation Can Generate Stronger Correlations between Range Size and Climatic Niches Than the Biological Signal — A Demonstration Using Bird and Mammal Range Maps

Abstract

Species’ geographic ranges could primarily be physiological tolerances drawn in space. Alternatively, geographic ranges could be only broadly constrained by physiological climatic tolerances: there could generally be much more proximate constraints on species’ ranges (dispersal limitation, biotic interactions, etc.) such that species often occupy a small and unpredictable subset of tolerable climates. In the literature, species’ climatic tolerances are typically estimated from the set of conditions observed within their geographic range. Using this method, studies have concluded that broader climatic niches permit larger ranges. Similarly, other studies have investigated the biological causes of incomplete range filling. But, when climatic constraints are measured directly from species’ ranges, are correlations between species’ range size and climate necessarily consistent with a causal link? We evaluated the extent to which variation in range size among 3277 bird and 1659 mammal species occurring in the Americas is statistically related to characteristics of species’ realized climatic niches. We then compared how these relationships differed from the ones expected in the absence of a causal link. We used a null model that randomizes the predictor variables (climate), while retaining their broad spatial autocorrelation structure, thereby removing any causal relationship between range size and climate. We found that, although range size is strongly positively related to climatic niche breadth, range filling and, to a lesser extent, niche position in nature, the observed relationships are not always stronger than expected from spatial autocorrelation alone. Thus, we conclude that equally strong relationships between range size and climate would result from any processes causing ranges to be highly spatially autocorrelated.

About the Authors:

Véronique Boucher-Lalonde

* E-mail: [email protected]

Current address: Department of Biology, Université Sherbrooke, Sherbrooke, Québec, Canada

Affiliation: Department of Biology, University of Ottawa, Ottawa, Ontario, Canada

David J. Currie

Affiliation: Department of Biology, University of Ottawa, Ottawa, Ontario, Canada

Introduction

Species’ geographic ranges have long been thought to be strongly limited by their tolerances of climate [1]. Temperature and water availability, in particular, are generally thought to be important axes of the fundamental niche [2], which encompasses all environmental conditions in which a species can survive [3]. The realized niche, i.e. the environmental conditions experienced within the species’ range, is set by the fundamental niche and reduced, or expanded, by dispersal, biotic interactions and the available set of conditions [4]. At macro-scales and coarse grains, climate is generally assumed to be an important range determinant, with range limits close to the species’ climatic tolerances [5–7].

Strong climatic constraints on individual species’ ranges are also thought to lead to higher-order patterns in nature such as the richness–climate relationship [8, 9] and inter-specific variations in range size [10, 11]. In particular, species that are able to tolerate broader environmental conditions are hypothesized to be more widely distributed: niche breadth positively affects range size [11, 12] through an increase in species’ potential range size. Biotic interactions (e.g., greater parasitism in warm climates) and dispersal limitation (e.g., colonization lags following glaciations) are hypothesized to exclude species from occurring in all suitable climatic conditions, such that species whose biology allows them to fill a greater proportion of their fundamental niche should have larger ranges [13–16]. Last, niche position, i.e. where species’ niche occurs along a given climatic gradient, has also been hypothesized to affect range size through various mechanisms [17–19].

Tests of these hypotheses relating range size to species’ biological characteristics typically rely on measuring climatic tolerances from observed geographic ranges. Early and Sax [20] noted that “the use of species’ distributions to measure species’ environmental tolerances is one of the most fundamental techniques in biogeography”. Physiological climatic tolerances are costly and time-consuming to measure (but see [9, 21, 22]) and in-lab measurements are not necessarily transferable to natural settings because of behaviour, habitat complexity or other confounding factors [23, 24]. It is therefore appealing to measure climatic tolerances from the observed species’ ranges themselves [13, 20, 25].

In particular, to test the hypothesis that wider tolerances increase species’ range size, species’ environmental or habitat niche breadths have generally been measured from the same geographic data as range size was calculated [10, 11]. As predicted, range size consistently increases with species' realized niche breadth [10, 11]. Slayter et al.’s [11] meta-analysis “Niche breadth predicts geographical range size: a general ecological pattern” called for a deeper investigation of the underlying ecological and evolutionary causes. But, in many of the studies compiled by Slayter et al. [11], niche breadth was measured from the same macro-scale data as range size; not independently. While four out of the five studies in Slatyer et al.’s [11] meta-analysis that measured range size and “environmental tolerance breadth” from the same species’ distribution data found the predicted positive relationship, only two out of six studies that used independent data found significant positive relationships (i.e. calcium requirement in Briers [26] and germination temperature in Calosi [21]). The results for “habitat breadth” in Slatyer et al. [11] were similar: 24/26 studies measuring range size and habitat breadth from the same species’ distribution data detected the predicted positive relationship, while only 7/18 studies that used independent data (e.g. microhabitat use, field guide classification, independent census data) did so. It is therefore possible that correlations between climatic niche breadth and range size are inflated when measured from the same distributional data.

Similarly, range-filling, the proportion of a species’ potential range that is occupied, has also been measured directly from the observed range of occupied climates, assuming that species can tolerate the macroclimatic conditions anywhere within this potential range, but not beyond it [27]; but see [24]. Hypotheses have then been proposed regarding the biological factors that could affect the extent of range-filling, e.g. post-glacial migration [27–29], competition (Pearson & Dawson 2003) or colonization ability [30]. Cold- and drought-tolerance have also been measured directly from species’ ranges in order to reach conclusions about potential causal mechanisms explaining variations in range size [17, 31].

Measuring climatic tolerance directly from species geographic ranges to test hypotheses about the biological determinants of species’ ranges rests on the critical assumption that species’ occupied climatic niche are generally closely related to their tolerances [25]. Although physiological climatic limits beyond which any given species cannot survive undoubtedly exist [32], the extent to which they actually impose range boundaries is less clear [33–38]. Most species ranges do not encompass all apparently suitable climates in the surrounding regions [13, 38, 39]. And, terrestrial species’ geographic ranges often fall short of (or extend beyond) the species’ physiological climatic limits [22, 35] such that there are not in climatic equilibrium [40].

If a circumscribed area such as a species’ range is placed randomly on a continent, it will encompass a range of climatic conditions—a realized climatic niche—even if there is no deterministic relationship between that range and the climate on which it was placed. It is already well recognized that the observed niche breadths of widespread species are taken from more locations than those of small-ranged species, which can generate an artefactual positive relationship between range size and niche breadth [41]. Studies have found that the relationship between range size and niche breadth holds even after correcting for this sampling bias [10, 11, 42]. But, what has been underappreciated is that species’ ranges, as well as extant climatic conditions, are highly spatially autocorrelated. Superimposing any spatially autocorrelated surfaces can lead to a correlation between them, even if there is no causal relationship [43–45].

Here, we evaluate the extent to which inter-specific variations in range size would be statistically related to species’ occupied climatic niches (i.e., climatic niche breadth, range filling and niche position) in the absence of any mechanistic connection between these variables. We then ask whether real (i.e. observed) range size–climatic niche relationships are stronger than expected from a randomization of the predictor variable (i.e., climatic niche breadth, range filling or niche position). To do this, we studied the breeding ranges of 3277 non-migratory birds and 1659 mammals occurring in the Americas. Note that, although animals may adjust their behaviour to occupy microclimatic conditions that are very different from the regional average, much of the literature on variation in range size and species’ climatic niche largely rely on coarse-grained environmental and animal distribution data, e.g. [46]. Birds and mammals sometimes show amongst the strongest effect sizes, e.g. Table S3 in [11]. We used macro-scale data on species’ geographic ranges, and we measured (macro)climatic niche characteristics from these. We used a null model that randomizes the spatial variation in climate, while maintaining its spatial autocorrelation [43–45]. This allowed us to examine the extent to which inter-specific variations in range size are related to characteristics of the climatic niche, versus spatial autocorrelation alone. Relationships stronger than the null expectations would be consistent with deterministic effects of climatic tolerances on species ranges.

Methods

The breeding ranges of 3277 non-migratory birds and 1659 mammals of the Americas were obtained from NatureServe [47, 48] and converted into raster format. We divided the Americas into 4141 equal-area cells of 10,000 km² following Boucher-Lalonde et al. [38], which is roughly the effective resolution of expert-drawn range maps [49]. A cell was considered occupied by a species when its breeding range overlapped the cell at least partly. Mean annual temperature and total precipitation were obtained from the WorldClim database at a 30 arc-seconds resolution [50]. Their means were calculated for each cell.

We measured the properties of species’ occupied temperature and precipitation niche as follows. We measured niche breadth as the range of occupied climate [51–53]. We measured climatic range filling as the realized/potential range size ratio following Svenning & Skov [13] where the potential range is assumed to include all cells with climatic conditions that fall within of conditions that the species’ occupies somewhere within its range. We restricted the potential range to be within the zoogeographical region(s) in which the species occurs (i.e. North American, Mexican, Panamanian, Amazonian, South American and/or Artico-Siberian; Holt et al. 2013), such that climatically suitable areas outside the zoogeographic regions in which the species occurs do not affect the measured potential range (see Fig F in S1 File for results where the potential range is not restricted). We measured niche position as the environmental centroid [54], which is here simply the mean climate of occupied cells. We also calculated the minimum and maximum occupied temperature and precipitation.

We generated null expectations for the empirical relationships between range size and measures of the climatic niche (described below). This was done by randomizing the climatic variables, while maintaining their autocorrelation structure among cells following Chapman [44]. The method is a simple extension of a resampling test, where the climate values are randomly rearranged among cells but where the resampled gradient is constrained to have a similar spatial correlation structure (estimated by empirical variograms) as the real gradients. This procedure removes any causal link between species’ ranges and climate. This could have also been done by placing species’ ranges randomly on the observed climatic gradients in North America[45, 55], but that procedure creates an artificial mid-domain effect [43] and is very sensitive to slight changes in how the null distribution is modelled [56], whereas we found that our conclusions were robust to the details of how the empirical variograms were modelled.

We then assessed how inter-specific variation in range sizes is related to species' climatic niche properties, using general linear models to explore these relationships. We then compared each range size–niche property relationship to the null expectation based on randomized climates (i.e. in which there is no causal link between range size and climatic niches). We also assessed how species’ climatic niche properties (niche breadth, range filling and niche position) covary among species, and we tested whether these correlations were different than those expected from spatial autocorrelation alone. In the literature, niche properties are sometimes correlated, and potential causal mechanisms are then inferred from these correlations [6]. Our aim here was therefore to test whether these correlations are stronger than what we can expect from spatial autocorrelation alone.

To implement Chapman’s [44] procedure to generate randomized gradients in temperature and precipitation, we first quantified spatial autocorrelation in temperature and precipitation in the Americas, at the spatial grain of 10⁴ km². For both temperature and precipitation, we modelled the spatial correlation structure of the data through an empirical variogram, using a powered-exponential function in the 'geoR' package [57] in R [58]. Then, we created a new variable that reproduced this spatial structure using stationary Gaussian Random Fields generated in the 'RandomFields' package [59]. This is a stochastic process by which random numbers are generated and mapped in space. The probability that a randomly generated value is assigned to a particular cell depends on the neighboring values through the variogram model. We then adjusted for the observed mean difference between coastal and inland cells following Chapman [44]. The values for each cell generated by the Gaussian Random Fields where then ranked. In each cell, the rank was replaced with the exact temperature or precipitation that has the same rank order in the real data. This reproduced the general, broad-scale spatial structure of the climatic variables (Fig A in S1 File), but obviously could not account for all its complexity. A detailed description of the algorithm and the R code to generate the simulated null gradients is available in Chapman’s [44] Appendix S3. We created one thousand iterations of these randomized temperature and precipitation surfaces within the Americas to serve as a null model in our analyses. Correlations between the real and randomized variables varied greatly, but they were centered on zero, as were the correlations between the randomized temperature and precipitation variables (Fig M in S1 File).

We then superimposed the observed breeding ranges of North American birds and mammals on each set of randomized climatic variables, and we calculated the niche properties for each species on the randomized climatic surfaces. Since the potential range was restricted to the zoogeographic region(s) in which each species occurs, our null model for the distribution of range filling was calculated by randomizing the climatic variables only within the zoogeographic region(s) in which the species occurs (but see Fig F in S1 File in which we relax this constraint). We then compared the slope and R² of the observed relationships to those obtained from our set of randomizations. Because some randomizations produced climate surfaces that were strongly correlated with observed climate, we repeated all our analyses with a set of 1000 randomizations in which we excluded all randomized climatic surfaces that were correlated with |r| > 0.2, with the real climatic variables. By restricting the randomizations to be uncorrelated with the real variables, we increase our ability to distinguish the real relationships from the null model. These results are presented in Figs I-L in S1 File and were qualitatively similar to those obtained from the full set of randomizations.

To test whether our approach has sufficient power to detect a signal of climate on species’ ranges and whether this is affected by intrinsic (endogenous) spatial autocorrelation (i.e. that arises from population processes such as dispersal), we simulated species’ ranges that are causally related to climate. Specifically, we performed two different simulation models in which we simulated species’ geographic ranges. In these, each species’ range size and climatic niche breadth were constrained by those measured in the real data. Consequently, the simulated niche characteristics were highly correlated to those observed (0.71≤ r ≤ 0.99). In the first simulation model, we probabilistically sampled from all cells that fell within the species’ climatic niche until the species’ range size was reached In this model, any spatial autocorrelation in the range is solely induced from the spatially autocorrelated climatic variables (i.e. exogenous spatial autocorrelation only). In the second simulation model, we additionally modelled strong intrinsic spatial autocorrelation in the species’ range through a simple spreading-dye model. Amongst the cells that fell within the species’ climatic niche, one was randomly seeded, and neighboring cells that fell within the species’ climatic niche were then sequentially selected to be occupied until the species’ range size was attained (if all neighboring cells within the species’ climatic niche were occupied and total range size was not yet attained, a new cell was randomly seeded). For both simulation models, we then compared the range size–climatic niche properties relationships and tested whether they were stronger than those that would have been observed in the 1000 randomizations of the climatic surfaces. This allowed us to test whether, given the variance in the data, a causal relationship between species’ ranges and climate could have been detected in the presence of, versus in the absence of, intrinsic spatial autocorrelation.

Results

Correlations among niche properties

Temperature and precipitation niche position are highly positively correlated for mammals (r = 0.68) and for birds (r = 0.64) in the Americas [6], which presumably at least partly reflects the fact that temperature and precipitation in the Americas are highly correlated (r = 0.62). Additionally, temperature and precipitation niche breadth are highly correlated (r = 0.77 for mammals and 0.65 for birds), as well as range filling of the temperature and precipitation niches (r = 0.91) (Table 1). Climatic niche breadth is weakly negatively related to temperature and precipitation niche position (Table 1). The strongest observed correlations were between niche breadth and the minimum temperature or precipitation occupied by the species (Fig N in S1 File). However, all the observed correlations fall either within the range expected from the randomization procedure, or are weaker than expected (Table 1; Figs H-G in S1 File).

[Figure omitted. See PDF.]

Table 1. Correlation coefficients from pairwise correlations between the different occupied niche properties for n = 3277 bird and 1659 mammal species within their American breeding ranges.

The measured niche properties are the temperature (T) and precipitation (P) breadth, mean position of occupied climates and range filling (i.e., the proportion of cells with suitable climates, in the zoogeographical region(s) in which the species occurs, that fall within the species’ range). The percentile of the simulated coefficients (n = 1000) in which falls the observed coefficient is presented in parentheses. All observed correlations either fell within the range of simulated coefficients or closer to a correlation of zero; none were stronger.

https://doi.org/10.1371/journal.pone.0166243.t001

Correlations with range size

Since the conclusions were similar for birds and mammals (see S1 File for detailed results), we present the pooled results below. Variations in range sizes are, by definition, entirely explained by the interaction of potential range size (i.e. the spatial extant of climates that are suitable based on occupancy) and range filling. Thus, it is unsurprising that observed range size was positively related to both species’ range filling (r² = 0.41 for temperature and 0.38 for precipitation) and climatic niche breath (r² = 0.34 for temperature and 0.22 for precipitation). Inter-specific variations in range size were only very weakly positively related to niche position (r² = 0.01 for both temperature and precipitation).

However, comparing the above results to null expectations, we found that these observed relationships either fell within, or were weaker (closer to zero) than, the relationships expected from a randomization of the predictor variables. Specifically, the observed slopes (Fig 1) and coefficients of determination (Fig 2) were within or below the range expected under a null model in which there is no causal link between range size and the characteristics of the climatic niche, but which conserves the spatial structure in both species’ ranges and climatic conditions. These conclusions were maintained when we analyzed the bird and mammal data separately (Figs B-E in S1 File). When we restricted our set of randomizations to be uncorrelated with the real climatic variables, we also found very similar results (Figs I-L in S1 File). In sum, we cannot reject the null hypothesis that the observed correlation between range size and niche breadth arose from spatial autocorrelation alone.

[Figure omitted. See PDF.]

Fig 1. Histograms of expected distribution of slopes between range size and a set of climatic niche characteristics under the null model.

The observed slope of the relationship between range size, measured as number of 10⁴-km² cells occupied, and (a, b) niche breadth, (c, d) niche position and (e, f) range filling of the occupied (a, c, e) temperature (°C) and (c, d, f) precipitation (mm) is shown by a vertical solid line. The data are from the American geographic ranges of 3277 non-migratory bird and 1659 mammal species, pooled together. The expected distribution of slopes was generated by randomizing the predictor variable, while maintaining spatial autocorrelation in the climatic variables. The vertical dashed lines represent the 2.5% and 97.5% thresholds of the randomizations.

https://doi.org/10.1371/journal.pone.0166243.g001

[Figure omitted. See PDF.]

Fig 2. Histograms of expected distribution of coefficients of determination (R²) between range size and a set of climatic niche characteristics under the null model.

The observed R² of the relationship between range size and (a, b) niche breadth, (c, d) niche position and (e, f) range filling of the occupied (a, c, e) temperature (°C) and (c, d, f) precipitation (mm) is shown by a vertical solid line. The data are from the American geographic range of 3277 non-migratory bird and 1659 mammal species, pooled together. The expected distribution of R² has been generated by randomizing the predictor variable, while maintaining spatial autocorrelation in the climatic variables. The vertical dashed lines represent the 2.5% and 97.5% thresholds of the randomizations.

https://doi.org/10.1371/journal.pone.0166243.g002

Our randomization procedures has the statistical power to detect the relationships between range size and climatic niche characteristics as being stronger than expected under our null model only in some cases (Figs 3 and 4). Specifically, when we simulated species’ ranges that are probabilistically determined by climate, while maintaining both range size and the range of climates occupied, the slopes and R² of relationships between range size and temperature niche breadth, temperature niche position as well as temperature and precipitation range filling were higher than those obtained by randomizing these simulated ranges (Figs 3 and 4). However, when we simulated species’ ranges through a spreading-dye model that introduces strong intrinsic spatial autocorrelation, the expected relationships between range size and climatic niche breadth became much stronger. Therefore, the observed slopes and R² of range size as a function of temperature niche breadth becomes undistinguishable from those obtained from the randomization procedure (Figs 3A and 4A). In other words, spatial autocorrelation, and intrinsic spatial autocorrelation in particular, makes causal relationships between range size and species’ climatic niches especially difficult to distinguish from those expected in the absence of causality.

[Figure omitted. See PDF.]

Fig 3. Histograms of expected slopes from power analyses that either incorporate intrinsic spatial autocorrelation (black) or not (gray).

Species geographic ranges were simulated to be probabilistically related to temperature and precipitation while maintaining the range size and climatic niche of real species. The observed slope of the relationship between range size, measured in number of cells occupied, and (a, b) niche breadth, (c, d) niche position and (e, f) range filling of the occupied (a, c, e) temperature (°C) and (c, d, f) precipitation (mm) are shown by vertical solid lines. Cells within a species’ climatic niche either had an equal probability of being occupied (i.e., ranges are not necessarily cohesive, shown in gray), or range cohesion, and thereby intrinsic spatial autocorrelation, was modelled through a spreading-dye algorithm (black). For each of the two sets of simulations, we then generated a histogram of the expected distribution of slope values by randomizing the predictor variable, while maintaining its spatial autocorrelation. The vertical dashed lines represent the 2.5% and 97.5% thresholds of the randomizations. This figure is equivalent to Fig 1, but for geographic ranges that have been simulated to be causally linked to macroclimatic variables.

https://doi.org/10.1371/journal.pone.0166243.g003

[Figure omitted. See PDF.]

Fig 4. Histograms of expected R² from power analyses that either incorporate intrinsic spatial autocorrelation (black) or not (gray).

Species geographic ranges that are probabilistically related to temperature and precipitation and that maintain the range size and climatic niche of real species were simulated. The observed R² of the relationship between range size, measured in number of cells occupied, and (a, b) niche breadth, (c, d) niche position and (e, f) range filling of the occupied (a, c, e) temperature (°C) and (c, d, f) precipitation (mm) are shown by vertical solid lines. Cells within a species’ climatic niche either had an equal probability of being occupied (gray) or intrinsic spatial autocorrelation was modelled through a spreading-dye algorithm (black). For each of the two sets of simulations, we then generated a histogram of the expected distribution of slope values by randomizing the predictor variable, while maintaining its spatial autocorrelation function. The vertical dashed lines represent the 2.5% and 97.5% thresholds of the randomizations. This figure is equivalent to Fig 2, but for geographic ranges that have been simulated to be causally linked to macroclimatic variables.

https://doi.org/10.1371/journal.pone.0166243.g004

Discussion

Temperature and precipitation are considered to be among the most important factors determining species’ range boundaries, especially at broad spatial scales [60]. They statistically explain much of the variation in the probability of occupancy of a region by a species [38, 39] as well as inter-specific variations in range size [12, 17]. Consistent with earlier literature, we find that the range size of bird and mammals species in the Americas is positively correlated with their climatic niche breadth, their range filling and their niche position. But, importantly, we found that similarly strong correlations between range size and species’ climatic niche properties (climatic niche breadth, niche position, and range filling) at macro-scales could arise from spatial autocorrelation alone. In the presence of strong spatial autocorrelation, correlations between a species’ range size and the climate niche characteristics that one can extract from this range are expected, regardless of any causal link, therefore providing an inconclusive test of causality.

Of course, species’ ranges could be highly spatially autocorrelated because they are structured by climate in the first place [41, 61]. However, geographic ranges could also be highly spatially autocorrelated because of processes that are independent of climate such as dispersal limitation. Here, we simply argue that the observed correlations between range sizes and niche properties cannot be used to support hypothesized causal links with climatic tolerance. This is because, although the relationships have slopes significantly steeper than zero, they are not stronger than those obtained from a null model that retains spatial autocorrelation. We note that we have only considered a single form of non-independence (spatial autocorrelation), but that other forms likely exacerbate the issue (e.g. phylogenetic relatedness). We argue that, to perform a strong test of the hypothesis that inter-specific differences in climatic niches characteristics cause variations in range size, it is critical that one measures range size and climatic tolerances independently (e.g. in lab measurements). Such tests are less common and the tested correlations have often been weaker or even absent [28, 62–64].

The conflation of environment and geography is a common problem in ecology that canlead to spurious inferences about the importance of deterministic processes in generating spatial patterns in species distributions [25]. It had already been noted that, in the presence of strong spatial autocorrelation, strong relationships between species’ geographic distributions and climate are expected [43–45], making evidence of a strong statistical relationship between the two (e.g., [38, 39]) particularly weak evidence of a causal relationship. Similarly, the conclusion that range size depends on climate largely rests on evidence of a statistical relationship where a species’ range size and climate niche are measured from the same coarse scale data, which also provides particularly weak evidence of causality. When relating range size to potential predictor variables, it is already common practice to account for the fact that wide-ranging species are sampled from more locations than small-ranging species, such that positive relationships are expected in the absence of any causal relationship [41]. We propose that future studies that rely on non-independent data should also explicitly account for spatial autocorrelation in the independent variable (when present).

Our results do raise the possibility that other factors than tolerance of the macro-climate generally explain variations in range size. Slatyer et al. [11] did find that diet breadth positively relate to range size (in addition to tolerance and habitat breadth). It is possible that geographic ranges are determined by a large set of factors. If the importance of these factors varies a lot among species, and through space and time, variations in range size could be at best only weakly predictable by a single predictor. Another (not mutually exclusive) hypothesis is that inter-specific variations in range size are partly driven by the ability of species to maintain high abundances in local communities. Local abundance has a strong negative effect on extinction rates [65] and positive effect on colonization rates of surrounding regions [66] such that range size generally increases with local abundance [67]. But, to rigorously test hypotheses about determinants of range size on equal grounds, and to avoid the risk of inflating the apparent effect of climate and other strongly spatially autocorrelated variables, range size and potential predictor variables need to be measured independently.

We have here shown that the relationship between range size and climate in the range could be explained by the spatial autocorrelation of the climate variables themselves and of the species’ ranges. Spatial autocorrelation greatly inflates the observed relationships between geographic ranges and climate [43, 45]. In sum, while climate obviously excludes many species from occurring in some areas [22], other factors could be more proximately determining species’ ranges, and we would nonetheless observe strong correlations between geographic ranges and climate. Climate could definitely be a major determinant of species’ range and range size, but we caution that some of the correlations commonly presented in the literature do not provide strong evidence of this.

Supporting Information

S1 File. Additional results Figures presenting examples of simulated climatic gradients, results for birds and mammals separately, results where the potential range is unconstrained, randomization distributions for the correlations between niche properties, results where randomizations are constrained to be poorly correlated with real climate and distributions of correlations between real and randomized gradients.

https://doi.org/10.1371/journal.pone.0166243.s001

(DOCX)

Author Contributions

1. Conceptualization: VBL DJC.

2. Formal analysis: VBL.

3. Methodology: VBL.

4. Writing – original draft: VBL.

5. Writing – review & editing: VBL DJC.

Citation: Boucher-Lalonde V, Currie DJ (2016) Spatial Autocorrelation Can Generate Stronger Correlations between Range Size and Climatic Niches Than the Biological Signal — A Demonstration Using Bird and Mammal Range Maps. PLoS ONE 11(11): e0166243. https://doi.org/10.1371/journal.pone.0166243

References

1. Andrewartha HG, Birch LC. The Distribution and Abundance of Animals. Chicago, Illinois: The University of Chicago Press; 1954. 782 p.

2. Kearney M, Porter WP. Mapping the fundamental niche: physiology, climate, and the distribution of a nocturnal lizard. Ecology. 2004;85(11):3119–31.

3. Hutchinson G. Concluding remarks. Cold Spring Harbour Symposia on Quantitative Biology1957. p. 415–27.

4. Peterson ATP, Soberón J, Pearson RG, Anderson RP, Martínez-Meyer E, Nakamura M, et al. Ecological Niches and Geographical Distributions: A Modeling Perspective: Princeton University Press; 2011.

5. Bozinovic F, Calosi P, Spicer JI. Physiological Correlates of Geographic Range in Animals. Annual Review of Ecology, Evolution, and Systematics, Vol 42. 2011;42:155–79.

6. Chejanovski ZA, Wiens JJ. Climatic niche breadth and species richness in temperate treefrogs. Journal of Biogeography. 2014;41(10):1936–46.

7. Strubbe D, Beauchard O, Matthysen E. Niche conservatism among non-native vertebrates in Europe and North America. Ecography. 2014:n/a-n/a.

8. Wiens JJ, Graham CH. Niche conservatism: Integrating evolution, ecology, and conservation biology. Annual Review of Ecology Evolution and Systematics. 2005;36:519–39.

9. Romdal TS, Araujo MB, Rahbek C. Life on a tropical planet: niche conservatism and the global diversity gradient. Global Ecology and Biogeography. 2013;22(3):344–50.

10. Morin X, Lechowicz MJ. Niche breadth and range area in North American trees. Ecography. 2013;36(3):300–12.

11. Slatyer RA, Hirst M, Sexton JP. Niche breadth predicts geographical range size: a general ecological pattern. Ecology Letters. 2013;16(8):1104–14. pmid:23773417

12. Stevens GC. The latitudinal gradient in geographical range—How so many species coexist in the Tropics. American Naturalist. 1989;133(2):240–56.

13. Svenning JC, Skov F. Limited filling of the potential range in European tree species. Ecology Letters. 2004;7(7):565–73.

14. Svenning JC, Skov F. Could the tree diversity pattern in Europe be generated by postglacial dispersal limitation? Ecology Letters. 2007;10(6):453–60. pmid:17498144

15. Soberon J, Peterson TA. Interpretation of models of fundamental ecological niches and species' distributional areas. Biodiversity Informatics. 2005;2:1–10.

16. Nogués-Bravo D, Pulido F, Araújo MB, Diniz-Filho JAF, García-Valdés R, Kollmann J, et al. Phenotypic correlates of potential range size and range filling in European trees. Perspectives in Plant Ecology, Evolution and Systematics. 2014.

17. Pither J. Climate tolerance and interspecific variation in geographic range size. Proceedings of the Royal Society of London Series B-Biological Sciences. 2003;270(1514):475–81.

18. Barnagaud J-Y, Devictor V, Jiguet F, Barbet-Massin M, Le Viol I, Archaux F. Relating Habitat and Climatic Niches in Birds. Plos One. 2012;7(3).

19. Laube I, Graham CH, Boehning-Gaese K. Intra-generic species richness and dispersal ability interact to determine geographic ranges of birds. Global Ecology and Biogeography. 2013;22(2):223–32.

20. Early R, Sax DF. Climatic niche shifts between species' native and naturalized ranges raise concern for ecological forecasts during invasions and climate change. Global Ecology and Biogeography. 2014;23(12):1356–65.

21. Calosi P, Bilton DT, Spicer JI, Votier SC, Atfield A. What determines a species’ geographical range? Thermal biology and latitudinal range size relationships in European diving beetles (Coleoptera: Dytiscidae). Journal of Animal Ecology. 2010;79(1):194–204. pmid:19761459

22. Sunday JM, Bates AE, Dulvy NK. Global analysis of thermal tolerance and latitude in ectotherms. Proceedings of the Royal Society B: Biological Sciences. 2011;278(1713):1823–30. pmid:21106582

23. Hoffmann AA, Chown SL, Clusella-Trullas S. Upper thermal limits in terrestrial ectotherms: how constrained are they? Functional Ecology. 2013;27(4):934–49.

24. Sunday JM, Bates AE, Kearney MR, Colwell RK, Dulvy NK, Longino JT, et al. Thermal-safety margins and the necessity of thermoregulatory behavior across latitude and elevation. Proceedings of the National Academy of Sciences. 2014;111(15):5610–5.

25. Warren DL, Cardillo M, Rosauer DF, Bolnick DI. Mistaking geography for biology: inferring processes from species distributions. Trends in Ecology & Evolution. 2014;29(10):572–80.

26. Briers RA. Range size and environmental calcium requirements of British freshwater gastropods. Global Ecology and Biogeography. 2003;12(1):47–51.

27. Skov F, Svenning JC. Potential impact of climatic change on the distribution of forest herbs in Europe. Ecography. 2004;27(3):366–80.

28. Sunday JM, Bates AE, Dulvy NK. Thermal tolerance and the global redistribution of animals. Nature Climate Change. 2012;2(9):686–90.

29. Svenning J-C, Normand S, Skov F. Postglacial dispersal limitation of widespread forest plant species in nemoral Europe. Ecography. 2008;31(3):316–26.

30. Schurr FM, Midgley GF, Rebelo AG, Reeves G, Poschlod P, Higgins SI. Colonization and persistence ability explain the extent to which plant species fill their potential range. Global Ecology and Biogeography. 2007;16(4):449–59.

31. Stephenson NL. Climatic control of vegetation distribution—The role of the water-balance. American Naturalist. 1990;135(5):649–70.

32. Araujo MB, Ferri-Yanez F, Bozinovic F, Marquet PA, Valladares F, Chown SL. Heat freezes niche evolution. Ecology Letters. 2013;16(9):1206–19. pmid:23869696

33. Bahn V, McGill BJ. Testing the predictive performance of distribution models. Oikos. 2013;122(3):321–31.

34. Ordonez A. Realized climatic niche of North American plant taxa lagged behind climate during the end of the Pleistocene. American Journal of Botany. 2013;100(7):1255–65. pmid:23825136

35. Hargreaves AL, Samis KE, Eckert CG. Are Species' Range Limits Simply Niche Limits Writ Large? A Review of Transplant Experiments beyond the Range. American Naturalist. 2014;183(2):157–73. pmid:24464192

36. Gouveia SF, Hortal J, Tejedo M, Duarte H, Cassemiro FAS, Navas CA, et al. Climatic niche at physiological and macroecological scales: the thermal tolerance–geographical range interface and niche dimensionality. Global Ecology and Biogeography. 2014;23:446–56.

37. Boucher-Lalonde V, Kerr JT, Currie DJ. Does climate limit species richness by limiting individual species' ranges? Proceedings of the Royal Society B-Biological Sciences. 2014;281(1776).

38. Boucher-Lalonde V, Morin A, Currie D. A consistent occupancy–climate relationship across birds and mammals of the Americas Oikos. 2014;123:1029–36.

39. Boucher-Lalonde V, Morin A, Currie D. How are tree species distributed in climatic space? A simple and general pattern. Global Ecology and Biogeography. 2012;21(12):1157–66.

40. Araujo MB, Pearson RG. Equilibrium of species' distributions with climate. Ecography. 2005;28(5):693–5.

41. Brown JH. On the relationship between abundance and distribution of species. American Naturalist. 1984;124(2):255–79.

42. Boulangeat I, Gravel D, Thuiller W. Accounting for dispersal and biotic interactions to disentangle the drivers of species distributions and their abundances. Ecology Letters. 2012;15(6):584–93. pmid:22462813

43. Beale CM, Lennon JJ, Gimona A. Opening the climate envelope reveals no macroscale associations with climate in European birds. Proceedings of the National Academy of Sciences of the United States of America. 2008;105(39):14908–12. pmid:18815364

44. Chapman DS. Weak climatic associations among British plant distributions. Global Ecology and Biogeography. 2010;19(6):831–41.

45. Bahn V, McGill BJ. Can niche-based distribution models outperform spatial interpolation? Global Ecology and Biogeography. 2007;16(6):733–42.

46. Munguia M, Peterson AT, Sanchez-Cordero V. Dispersal limitation and geographical distributions of mammal species. Journal of Biogeography. 2008;35(10):1879–87.

47. Patterson B, Ceballos G, Sechrest W, Tognelli M, Brooks T, Luna L, et al. Digital Distribution Maps of the Mammals of the Western Hemisphere, version 3.0. NatureServe, Arlington, Virginia, USA2007.

48. Ridgely RS, Allnutt TF, Brooks T, McNicol DK, Mehlman DW, Young BE, et al. Digital Distribution Maps of the Birds of the Western Hemisphere, version 3.0. Nature Serve, Arlington, Virginia, USA2007.

49. Hurlbert AH, Jetz W. Species richness, hotspots, and the scale dependence of range maps in ecology and conservation. Proceedings of the National Academy of Sciences of the United States of America. 2007;104(33):13384–9. pmid:17686977

50. Hijmans RJ, Cameron SE, Parra JL, Jones PG, Jarvis A. Very high resolution interpolated climate surfaces for global land areas. International Journal of Climatology. 2005;25(15):1965–78.

51. Pannek A, Ewald J, Diekmann M. Resource-based determinants of range sizes of forest vascular plants in Germany. Global Ecology and Biogeography. 2013;22(8):1019–28.

52. Quintero I, Wiens JJ. What determines the climatic niche width of species? The role of spatial and temporal climatic variation in three vertebrate clades. Global Ecology and Biogeography. 2013;22(4):422–32.

53. Nakazawa Y. Niche breadth, environmental landscape, and physical barriers: their importance as determinants of species distributions. Biological Journal of the Linnean Society. 2013;108(2):241–50.

54. Soininen J, Heino J, Lappalainen J, Virtanen R. Expanding the ecological niche approach: Relationships between variability in niche position and species richness. Ecological Complexity. 2011;8(1):130–7.

55. McGill B, Collins C. A unified theory for macroecology based on spatial patterns of abundance. Evolutionary Ecology Research. 2003;5(4):469–92.

56. Beale CM, Lennon JJ, Gimona A. European bird distributions still show few climate associations. Proceedings of the National Academy of Sciences. 2009;106(16):E41–E3.

57. Ribeiro PJ Jr. geoR: a package for geostatistical analysis R-NEWS. 2001;1(2):15–8.

58. Team RC. R: A language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria, URL http://www.R-project.org/. 2013.

59. Schlather M, Malinowski A, Oesting M, Boecker D, Strokorb K, Engelke S, et al. RandomFields: Simulation and Analysis of Random Fields. R package version 3.0.10. http://CRAN.R-project.org/package=RandomFields 2014.

60. Wiens JJ. The causes of species richness patterns across space, time and clades and the role of "ecological limits". Quarterly Review of Biology. 2011;86(2):75–96. pmid:21800635

61. Diniz‐Filho JAF, Bini LM, Hawkins BA. Spatial autocorrelation and red herrings in geographical ecology. Global Ecology and Biogeography. 2003;12(1):53–64.

62. Calosi P, Bilton DT, Spicer JI. Thermal tolerance, acclimatory capacity and vulnerability to global climate change2008 2008-02-23 00:00:00. 99–102 p.

63. Luna B, Moreno J. Range-size, local abundance and germination niche-breadth in Mediterranean plants of two life-forms. Plant Ecology. 2010;210(1):85–95.

64. Bates AE, McKelvie CM, Sorte CJB, Morley SA, Jones NAR, Mondon JA, et al. Geographical range, heat tolerance and invasion success in aquatic species. Proceedings of the Royal Society of London B: Biological Sciences. 2013;280(1772).

65. Boucher-Lalonde V, Theriault FL, Currie DJ. Can climate explain interannual local extinctions among bird species? Journal of Biogeography. 2014;41(3):443–51.

66. Hubbell SP. The Unified Neutral Theory of Biodiversity and Biogeography. Princeton, N.J.: Princeton University Press; 2001.

67. Gaston KJ, Blackburn TM, Lawton JH. Interspecific abundance range size relationships: An appraisal of mechanisms. Journal of Animal Ecology. 1997;66(4):579–601.