Equilibrium of species’ distributions with climate

Abstract

The degree to which species’ distributions are at equilibrium with current climate is an important issue in the ecological literature (Davis 1986, Gaston 2003). Distinguishing the relative roles of present and past climates in determining species’ distributions is not only of theoretical interest, but is key to understanding responses of species to climate change. Recent work has demonstrated significant departures from equilibrium with current climate for tree distributions in Europe (Svenning and Skov 2004), yet differences in the degree of equilibrium can be expected across organisms whose dispersal abilities (and hence ability to track changing climate conditions) vary greatly. Here we assess differences in equilibrium between present-day climate and vascular plant, breeding bird, amphibian and reptile assemblages in Europe. Our analyses included all know European species of breeding birds (Hagemeijer and Blair 1997), amphibians and reptiles (Gasc et al. 1997), and ca 20% of the European vascular flora (Jalas and Suominen 1972–1996, Lahti and Lampinen 1999). Assemblages of plants and breeding birds were found to be relatively closer to equilibrium than assemblages of reptiles and amphibians. Results suggest that responses of plant and bird species to climate change are more likely to be accurately forecasted by models correlating present-day distributions with climate, and that reptile and amphibian species will be least capable of shifting distributions, making them most vulnerable to rapid environmental changes. Species are said to be at equilibrium with climate if they occur in all climatically suitable areas whilst being absent from all unsuitable ones (Hutchinson 1957). The observation that species are absent from many suitable areas is trivial. The critical question is how distant from equilibrium are current distributions. Accurate estimates of species’ limits of tolerance to climate can be obtained by means of controlled experiments (Chuine and Beaubien 2001). However, these are expensive and time consuming and attempts to measure departures from equilibrium for large numbers of species have previously used the “bioclimatic envelope” approach, whereby observed species’ distributions are correlated with environmental variables to approximate the ecological requirements of organisms (for review see Pearson and Dawson 2003). With this approach, the ratio between observed and modelled distributions is interpreted as indicating departures of species distributions from equilibrium (Svenning and Skov 2004). By relying on observed distributions (which inherently reflect multiple range determinants, both historical and ecological) to determine ecological requirements, this methodology is likely to underestimate the true range of climate variation that species are able to tolerate. An alternative approach consists of measuring patterns of covariation between species assemblages and climate (Ferrier et al. 2002). It is assumed that if species are at equilibrium with current climate, then patterns of covariation between assemblage composition and climate will be high. A weak covariation between climate and biota may indicate that species tend not to occupy all available climate spaces due to historical and/or ecological factors (notably competition and limited dispersal ability). In some cases it could also imply that the selected climate parameters do not adequately describe the climatic determinants of the distribution of assemblages. However this possibility is inherent to any correlational approach (including bioclimate-envelope modelling) and the choice of climate variables needs to be supported from existing knowledge of species ecologies. This pattern-based approach can be explored by measuring the degree of association between two dissimilarity matrices using a suitable correlation coefficient (also known as mantel tests, see Manly 1998). Here we use Spearman's rank correlation coefficients (ρ) to measure the degree of association between species (Bray-Curtis) dissimilarity matrices and (standardised Euclidean) distances in a climate matrix. The climate matrix was constructed with seven climate variables (described in Araújo et al. 2004a) identified from the literature as being key determinants of physiological processes limiting distributions of plant (Prentice et al. 1992), bird (Lennon et al. 2000) and amphibian (Carey and Alexander 2003) distributions. An exploratory analysis of the potential impacts of climate change on European herptile species distributions (Araújo et al. unpubl.), also shows that climate variables selected by models tend to be similar for reptiles and amphibians. The significance of the climate-biota relationship can be determined by comparison with the distribution of ρ found by randomly reallocating the order of the elements in one of the matrices. Under the null hypothesis of no relation between two matrices, ρ will be approximately zero. Here, its null distribution, eith