Uncertainty of bioclimate envelope models based on the geographical distribution of species

Aim We explored the effects of prevalence, latitudinal range and spatial autocorrelation of species distribution patterns on the accuracy of bioclimate envelope models of butterflies. Location Finland, northern Europe. Methods The data of a national butterfly atlas survey (NAFI) carried out in 1991–2003 with a resolution of 10 × 10 km were used in the analyses. Generalized additive models (GAM) were constructed, for each of 98 species, to estimate the probability of occurrence as a function of climate variables. Model performance was measured using the area under the curve (AUC) of a receiver operating characteristic (ROC) plot. Observed differences in modelling accuracy among species were related to the species’ geographical attributes using multivariate GAM. Results Accuracies of the climate–butterfly models varied from low to very high (AUC values 0.59–0.99), with a mean of 0.79. The modelling performance was related negatively to the latitudinal range and prevalence, and positively to the spatial autocorrelation of the species distribution. These three factors accounted for 75.2% of the variation in the modelling accuracy. Species at the margin of their range or with low prevalence were better predicted than widespread species, and species with clumped distributions better than scattered dispersed species. Main conclusions The results from this study indicate that species’ geographical attributes highly influence the behaviour and uncertainty of species–climate models, which should be taken into account in biogeographical modelling studies and assessments of climate change impacts.     

Does seasonal fine‐tuning of climatic variables improve the performance of bioclimatic envelope models for migratory birds?

We examined the influence of 'seasonal fine-tuning' of climatic variables on the performance of bioclimatic envelope models of migrating birds. Using climate data and national bird atlas data from a 10 × 10 km uniform grid system in Finland, we tested whether the replacement of one 'baseline' set of variables including summer (June–August) temperature and precipitation variables with climate variables tailored ('fine-tuned') for each species individually improved the bird-climate models. The fine-tuning was conducted on the basis of time of arrival and early breeding of the species. Two generalized additive models (GAMs) were constructed for each of the 63 bird species studied, employing (1) the baseline climate variables and (2) the fine-tuned climate variables. Model performance was measured as explanatory power (deviance change) and predictive power (area under the curve; AUC) statistics derived from cross-validation. Fine-tuned climate variables provided, in many cases, statistically significantly improved model performance compared to using the same baseline set of variables for all the species. Model improvements mainly concerned bird species arriving and starting their breeding in May–June. We conclude that the use of the fine-tuned climate variables tailored for each species individually on the basis of their arrival and critical breeding periods can provide important benefits for bioclimatic modelling.     

Determinants of the biogeographical distribution of butterflies in boreal regions

Aim  We explored the relative contributions of climatic and land-cover factors in explaining the distribution patterns of butterflies in a boreal region.
Location  Finland, northern Europe.
Methods  Data from a national butterfly atlas survey carried out during 1991–2003, with a 10-km grain grid system, were used in these analyses. We used generalized additive models (GAM) and hierarchical partitioning (HP) to explore the main environmental correlates (climate and land-cover) of the realized niches of 98 butterfly species. The accuracy of the distribution models (GAMs) was validated by resubstitution and cross-validation approaches, using the area under the curve (AUC) derived from the receiver operating characteristic (ROC) plots.
Results  Predictive accuracies of the 98 individual environment–butterfly models varied from low to very high (cross-validated AUC values 0.48–0.99), with a mean of 0.79. The results of both the GAM and HP analyses were broadly concordant. Most of the variation in butterfly distributions is associated with growing degree-days, mean temperature of the coldest month and cover of built-up area in all six phylogenetic groups (butterfly families). There were no statistically significant differences in predictive accuracy among the different butterfly families.
Main conclusions  About three-quarters of the distributions of butterfly species in Finland appear to be governed principally by climatic, predominantly temperature-related, factors. This indicates that many butterfly species may respond rapidly to the projected climate change in boreal regions. By determining the ecological niches of multiple species, we can project their range shifts in response to changes in climate and land-cover, and identify species that are particularly sensitive to forecasted global changes.     

Disregarding topographical heterogeneity biases species turnover assessments based on bioclimatic models

We investigated whether the inclusion of topographical heterogeneity in bioclimatic envelope models would significantly alter predictions of climate change – induced broad-scale butterfly species range size changes in Europe. Using generalized additive models, and data on current climate and species distributions and two different climate scenarios (HadCM3A2 and HadCM3B2) for the period 2051–2080, we developed predictions of the currently suitable area and potential range size changes of 100 European butterfly species. The inclusion of elevation range increased the predictive accuracy of climate-only models for 86 of the 100 species. The differences in projected future distributions were most notable in mountainous areas, where the climate–topography models projected only ca. half of the species losses than the climate-only models. By contrast, climate–topography models estimated double the losses of species than climate-only models in the flatlands regions. Our findings suggest that disregarding topographical heterogeneity may cause a significant source of error in broad-scale bioclimatic modelling. Mountainous regions are likely to be even more important for future conservation of species than had until now been predicted, based on bioclimatic envelope models that did not take an explicit account of elevational range of grid squares.     

Heritability of flight and resting metabolic rates in the Glanville fritillary butterfly

Dispersal capacity is a key life‐history trait especially in species inhabiting fragmented landscapes. Evolutionary models predict that, given sufficient heritable variation, dispersal rate responds to natural selection imposed by habitat loss and fragmentation. Here, we estimate phenotypic variance components and heritability of flight and resting metabolic rates (RMRs) in an ecological model species, the Glanville fritillary butterfly, in which flight metabolic rate (FMR) is known to correlate strongly with dispersal rate. We modelled a two‐generation pedigree with the animal model to distinguish additive genetic variance from maternal and common environmental effects. The results show that FMR is significantly heritable, with additive genetic variance accounting for about 40% of total phenotypic variance; thus, FMR has the potential to respond to selection on dispersal capacity. Maternal influences on flight metabolism were negligible. Heritability of flight metabolism was context dependent, as in stressful thermal conditions, environmentally induced variation dominated over additive genetic effects. There was no heritability in RMR, which was instead strongly influenced by maternal effects. This study contributes to a mechanistic understanding of the evolution of dispersal‐related traits, a pressing question in view of the challenges posed to many species by changing climate and fragmentation of natural habitats.     

Differences in spatial predictions among species distribution modeling methods vary with species traits and environmental predictors

Prediction maps produced by species distribution models (SDMs) influence decision‐making in resource management or designation of land in conservation planning. Many studies have compared the prediction accuracy of different SDM modeling methods, but few have quantified the similarity among prediction maps. There has also been little systematic exploration of how the relative importance of different predictor variables varies among model types and affects map similarity. Our objective was to expand the evaluation of SDM performance for 45 plant species in southern California to better understand how map predictions vary among model types, and to explain what factors may affect spatial correspondence, including the selection and relative importance of different environmental variables. Four types of models were tested. Correlation among maps was highest between generalized linear models (GLMs) and generalized additive models (GAMs) and lowest between classification trees and GAMs or GLMs. Correlation between Random Forests (RFs) and GAMs was the same as between RFs and classification trees. Spatial correspondence among maps was influenced the most by model prediction accuracy (AUC) and species prevalence; map correspondence was highest when accuracy was high and prevalence was intermediate (average prevalence for all species was 0.124). Species functional type and the selection of climate variables also influenced map correspondence. For most (but not all) species, climate variables were more important than terrain or soil in predicting their distributions. Environmental variable selection varied according to modeling method, but the largest differences were between RFs and GLMs or GAMs. Although prediction accuracy was equal for GLMs, GAMs, and RFs, the differences in spatial predictions suggest that it may be important to evaluate the results of more than one model to estimate the range of spatial uncertainty before making planning decisions based on map outputs. This may be particularly important if models have low accuracy or if species prevalence is not intermediate.     

Co-occurrence of tree species at fine spatial scale in a woodland cerrado,southeastern Brazil

Species co-occurrence at fine spatial scales is expected to be nonrandom in relation to species phylogenetic relatedness and functional similarity. On the one hand, closely related species that occur together and experience similar environmental conditions are likely to share phenotypic traits due to the process of environmental filtering. On the other hand, species that are too similar are unlikely to co-occur due to competitive exclusion. We surveyed a woodland cerrado, southeastern Brazil, to test whether co-occurrence in tree species shows functional or phylogenetic structuring at fine spatial scale. Searching for correlations between an index of species co-occurrence and both functional trait differences and phylogenetic distances, we provided evidence for a predominant role of environment filters in determining the co-occurrence of functionally similar tree species in cerrado. However, we did not find any effect of phylogenetic relatedness on tree species co-occurrence. We suggest that the phylogenetic relatedness of co-occurring cerrado tree species did not present a pattern, because the species functional traits were randomly distributed on the phylogeny. Thus, phylogenetic relatedness and functional similarity do not seem to limit the co-occurrence at fine spatial scale of cerrado tree species.     

The role of land cover in bioclimatic models depends on spatial resolution

Aim We explored the importance of climate and land cover in bird species distribution models on multiple spatial scales. In particular, we tested whether the integration of land cover data improves the performance of pure bioclimatic models. Location Finland, northern Europe. Methods The data of the bird atlas survey carried out in 1986–89 using a 10 × 10 km uniform grid system in Finland were employed in the analyses. Land cover and climatic variables were compiled using the same grid system. The dependent and explanatory variables were resampled to 20‐km, 40‐km and 80‐km resolutions. Generalized additive models (GAM) were constructed for each of the 88 land bird species studied in order to estimate the probability of occurrence as a function of (1) climate and (2) climate and land cover variables. Model accuracy was measured by a cross‐validation approach using the area under the curve (AUC) of a receiver operating characteristic (ROC) plot. Results In general, the accuracies of the 88 bird–climate models were good at all studied resolutions. However, the inclusion of land cover increased the performance of 79 and 78 of the 88 bioclimatic models at 10‐km and 20‐km resolutions, respectively. There was no significant improvement at the 40‐km resolution. In contrast to the finer resolutions, the inclusion of land cover variables decreased the modelling accuracy at 80km resolution. Main conclusions Our results suggest that the determinants of bird species distributions are hierarchically structured: climatic variables are large‐scale determinants, followed by land cover at finer resolutions. The majority of the land bird species in Finland are rather clearly correlated with climate, and bioclimate envelope models can provide useful tools for identifying the relationships between these species and the environment at resolutions ranging from 10 km to 80 km. However, the notable contribution of land cover to the accuracy of bioclimatic models at 10–20‐km resolutions indicates that the integration of climate and land cover information can improve our understanding and model predictions of biogeographical patterns under global change.     

Predicting the distribution of shrub species in southern California from climate and terrain-derived variables

Abstract. Generalized additive, generalized linear, and classification tree models were developed to predict the distribution of 20 species of chaparral and coastal sage shrubs within the southwest ecoregion of California. Mapped explanatory variables included bioclimatic attributes related to primary environmental regimes: averages of annual precipitation, minimum temperature of the coldest month, maximum temperature of the warmest month, and topographically-distributed potential solar insolation of the wettest quarter (winter) and of the growing season (spring). Also tested for significance were slope angle (related to soil depth) and the geographic coordinates of each observation. Models were parameterized and evaluated based on species presence/absence data from 906 plots surveyed on National Forest lands. Although all variables were significant in at least one of the species’ models, those models based only on the bioclimatic variables predicted species presence with 3–26% error. While error would undoubtedly be greater if the models were evaluated using independent data, results indicate that these models are useful for predictive mapping – for interpolating species distribution data within the ecoregion. All three methods produced models with similar accuracy for a given species; GAMs were useful for exploring the shape of the response functions, GLMs allowed those response functions to be parameterized and their significance tested, and classification trees, while some-times difficult to interpret, yielded the lowest prediction errors (lower by 3–5%).     

Predicting range expansion of the map butterfly in Northern Europe using bioclimatic models

The two main goals of this study are: (i) to examine the range shifts of a currently northwards expanding species, the map butterfly (Araschnia levana), in relation to annual variation in weather, and (ii) to test the capability of a bioclimatic envelope model, based on broad-scale European distribution data, to predict recent distributional changes (2000–2004) of the species in Finland. A significant relationship between annual maximum dispersal distance of the species and late summer temperature was detected. This suggests that the map butterfly has dispersed more actively in warmer rather than cooler summers, the most notable dispersal events being promoted by periods of exceptionally warm weather and southerly winds. The accuracy of the broad-scale bioclimatic model built for the species with European data using Generalized Additive Models (GAM) was good based on split-sample evaluation for a single period. However, the model’s performance was poor when applied to predict range shifts in Finland. Among the many potential explanations for the poor success of the transferred bioclimatic model, is the fact that bioclimatic envelope models do not generally account for species dispersal. This and other uncertainties support the view that bioclimatic models should be applied with caution when they are used to project future range shifts of species.     

Modelling the occurrence of threatened plant species in taiga landscapes: methodological and ecological perspectives

Aim Understanding the spatial patterns of species distribution and predicting the occurrence of high biological diversity and rare species are central themes in biogeography and environmental conservation. The aim of this study was to model and scrutinize the relative contributions of climate, topography, geology and land‐cover factors to the distributions of threatened vascular plant species in taiga landscapes in northern Finland. Location North‐east Finland, northern Europe. Methods The study was performed using a data set of 28 plant species and environmental variables at a 25‐ha resolution. Four different stepwise selection algorithms [Akaike information criterion (AIC), Bayesian information criterion (BIC), adaptive backfitting, cross selection] with generalized additive models (GAMs) were fitted to identify the main environmental correlates for species occurrences. The accuracies of the distribution models were evaluated using fourfold cross‐validation based on the area under the curve (AUC) derived from receiver operating characteristic plots. The GAMs were tentatively extrapolated to the whole study area and species occurrence probability maps were produced using GIS techniques. The effect of spatial autocorrelation on the modelling results was also tested by including autocovariate terms in the GAMs. Results According to the AUC values, the model performance varied from fair to excellent. The AIC algorithm provided the highest mean performance (mean AUC = 0.889), whereas the lowest mean AUC (0.851) was obtained from BIC. Most of the variation in the distribution of threatened plant species was related to growing degree days, temperature of the coldest month, water balance, cover of mire and mean elevation. In general, climate was the most powerful explanatory variable group, followed by land cover, topography and geology. Inclusion of the autocovariate only slightly improved the performance of the models and had a minor effect on the importance of the environmental variables. Main conclusions The results confirm that the landscape‐scale distribution patterns of plant species can be modelled well on the basis of environmental parameters. A spatial grid system with several environmental variables derived from remote sensing and GIS data was found to produce useful data sets, which can be employed when predicting species distribution patterns over extensive areas. Landscape‐scale maps showing the predicted occurrences of individual or multiple threatened plant species may provide a useful basis for focusing field surveys and allocating conservation efforts.     

Refining predictions of climate change impacts on plant species distribution through the use of local statistics

Bioclimate envelope models are often used to predict changes in species distribution arising from changes in climate. These models are typically based on observed correlations between current species distribution and climate data. One limitation of this basic approach is that the relationship modelled is assumed to be constant in space; the analysis is global with the relationship assumed to be spatially stationary. Here, it is shown that by using a local regression analysis, which allows the relationship under study to vary in space, rather than conventional global regression analysis it is possible to increase the accuracy of bioclimate envelope modelling. This is demonstrated for the distribution of Spotted Meddick in Great Britain using data relating to three time periods, including predictions for the 2080s based on two climate change scenarios. Species distribution and climate data were available for two of the time periods studied and this allowed comparison of bioclimate envelope model outputs derived using the local and global regression analyses. For both time periods, the area under the receiver operating characteristics curve derived from the analysis based on local statistics was significantly higher than that from the conventional global analysis; the curve comparisons were also undertaken with an approach that recognised the dependent nature of the data sets compared. Marked differences in the future distribution of the species predicted from the local and global based analyses were evident and highlight a need for further consideration of local issues in modelling ecological variables.     

Inclusion of soil data improves the performance of bioclimatic envelope models for insect species distributions in temperate Europe

Aim To analyse the effect of the inclusion of soil and land‐cover data on the performance of bioclimatic envelope models for the regional‐scale prediction of butterfly (Rhopalocera) and grasshopper (Orthoptera) distributions. Location Temperate Europe (Belgium). Methods Distributional data were extracted from butterfly and grasshopper atlases at a resolution of 5 km for the period 1991–2006 in Belgium. For each group separately, the well‐surveyed squares (n = 366 for butterflies and n = 322 for grasshoppers) were identified using an environmental stratification design and were randomly divided into calibration (70%) and evaluation (30%) datasets. Generalized additive models were applied to the calibration dataset to estimate occurrence probabilities for 63 butterfly and 33 grasshopper species, as a function of: (1) climate, (2) climate and land‐cover, (3) climate and soil, and (4) climate, land‐cover and soil variables. Models were evaluated as: (1) the amount of explained deviance in the calibration dataset, (2) Akaike’s information criterion, and (3) the number of omission and commission errors in the evaluation dataset. Results Information on broad land‐cover classes or predominant soil types led to similar improvements in the performance relative to the climate‐only models for both taxonomic groups. In addition, the joint inclusion of land‐cover and soil variables in the models provided predictions that fitted more closely to the species distributions than the predictions obtained from bioclimatic models incorporating only land‐cover or only soil variables. The combined models exhibited higher discrimination ability between the presence and absence of species in the evaluation dataset. Main conclusions These results draw attention to the importance of soil data for species distribution models at regional scales of analysis. The combined inclusion of land‐cover and soil data in the models makes it possible to identify areas with suitable climatic conditions but unsuitable combinations of vegetation and soil types. While contingent on the species, the results indicate the need to consider soil information in regional‐scale species–climate impact models, particularly when predicting future range shifts of species under climate change.     

Phenology responses of temperate butterflies to latitude depend on ecological traits

Global change influences species’ seasonal occurrence, or phenology. In cold‐adapted insects, the activity is expected to start earlier with a warming climate, but contradictory evidence exists, and the reactions may be linked to species‐specific traits. Using data from the GBIF database, we selected 105 single‐brooded Holarctic butterflies inhabiting broad latitudinal ranges. We regressed patterns of an adult flight against latitudes of the records, controlling for altitude and year effects. Species with delayed flight periods towards the high latitudes, or stable flight periods across latitudes, prevailed over those that advanced their flight towards the high latitudes. The responses corresponded with the species’ seasonality (flight of early season species was delayed and flight of summer species was advanced at high latitudes) and oceanic vs. continental climatic niches (delays in oceanic, stability in continental species). Future restructuring of butterfly seasonal patterns in high latitudes will reflect climatic niches, and hence the evolutionary history of participating species.     

Do traits of plant species predict the efficacy of species distribution models for finding new occurrences?

Species distribution models (SDMs) are used to test ecological theory and to direct targeted surveys for species of conservation concern. Several studies have tested for an influence of species traits on the predictive accuracy of SDMs. However, most used the same set of environmental predictors for all species and/or did not use truly independent data to test SDM accuracy. We built eight SDMs for each of 24 plant species of conservation concern, varying the environmental predictors included in each SDM version. We then measured the accuracy of each SDM using independent presence and absence data to calculate area under the receiver operating characteristic curve (AUC) and true positive rate (TPR). We used generalized linear mixed models to test for a relationship between species traits and SDM accuracy, while accounting for variation in SDM performance that might be introduced by different predictor sets. All traits affected one or both SDM accuracy measures. Species with lighter seeds, animal‐dispersed seeds, and a higher density of occurrences had higher AUC and TPR than other species, all else being equal. Long‐lived woody species had higher AUC than herbaceous species, but lower TPR. These results support the hypothesis that the strength of species–environment correlations is affected by characteristics of species or their geographic distributions. However, because each species has multiple traits, and because AUC and TPR can be affected differently, there is no straightforward way to determine a priori which species will yield useful SDMs based on their traits. Most species yielded at least one useful SDM. Therefore, it is worthwhile to build and test SDMs for the purpose of finding new populations of plant species of conservation concern, regardless of these species’ traits.     

A meta-analysis of the traits affecting dispersal ability in butterflies: can wingspan be used as a proxy?

1. Dispersal ability of a species is a key ecological characteristic, affecting a range of processes from adaptation, community dynamics and genetic structure, to distribution and range size. It is determined by both intrinsic species traits and extrinsic landscape-related properties. 2. Using butterflies as a model system, the following questions were addressed: (i) given similar extrinsic factors, which intrinsic species trait(s) explain dispersal ability? (ii) can one of these traits be used as a proxy for dispersal ability? (iii) the effect of interactions between the traits, and phylogenetic relatedness, on dispersal ability. 3. Four data sets, using different measures of dispersal, were compiled from published literature. The first data set uses mean dispersal distances from capture-mark-recapture studies, and the other three use mobility indices. Data for six traits that can potentially affect dispersal ability were collected: wingspan, larval host plant specificity, adult habitat specificity, mate location strategy, voltinism and flight period duration. Each data set was subjected to both unifactorial, and multifactorial, phylogenetically controlled analyses. 4. Among the factors considered, wingspan was the most important determinant of dispersal ability, although the predictive powers of regression models were low. Voltinism and flight period duration also affect dispersal ability, especially in case of temperate species. Interactions between the factors did not affect dispersal ability, and phylogenetic relatedness was significant in one data set. 5. While using wingspan as the only proxy for dispersal ability maybe problematic, it is usually the only easily accessible species-specific trait for a large number of species. It can thus be a satisfactory proxy when carefully interpreted, especially for analyses involving many species from all across the world.     

Equilibrium of global amphibian species distributions with climate

A common assumption in bioclimatic envelope modeling is that species distributions are in equilibrium with contemporary climate. A number of studies have measured departures from equilibrium in species distributions in particular regions, but such investigations were never carried out for a complete lineage across its entire distribution. We measure departures of equilibrium with contemporary climate for the distributions of the world amphibian species. Specifically, we fitted bioclimatic envelopes for 5544 species using three presence-only models. We then measured the proportion of the modeled envelope that is currently occupied by the species, as a metric of equilibrium of species distributions with climate. The assumption was that the greater the difference between modeled bioclimatic envelope and the occupied distribution, the greater the likelihood that species distribution would not be at equilibrium with contemporary climate. On average, amphibians occupied 30% to 57% of their potential distributions. Although patterns differed across regions, there were no significant differences among lineages. Species in the Neotropic, Afrotropics, Indo-Malay, and Palaearctic occupied a smaller proportion of their potential distributions than species in the Nearctic, Madagascar, and Australasia. We acknowledge that our models underestimate non equilibrium, and discuss potential reasons for the observed patterns. From a modeling perspective our results support the view that at global scale bioclimatic envelope models might perform similarly across lineages but differently across regions.     

The role of niche breadth, resource availability and range position on the life history of butterflies

We analysed the relationship between three life history characteristics (mobility, length of flight period and body size) and niche breadth (larval host plant specificity and adult habitat breadth), resource availability (distribution and abundance of host plants) and range position (distance between the northernmost distribution record and southernmost point of Finland) of the butterfly fauna of Finland. The data is based on literature and questionnaires. Often in across species studies phylogeny may create spurious relationships between life-history and ecological variables. We took the phylogenetic relatedness of butterfly species into account by analysing the data with phylogenetically independent contrasts (CAIC method). Butterfly mobility was positively related to the niche breadth, resource availability and range position. The length of the flight period was negatively related to the range position, indicating that the species at the northern edge of their distribution range have shorter flight period than species which are further way from the range edge. After controlling for the phylogenetic relatedness we found no significant correlations between body size and niche breadth, resource availability or range position. We suggest that the relationship between the length of the flight period and range position may arise as a consequence of lower hatching asynchrony in edge species as a result of lower environmental variance in larval growth conditions. Our results on the mobility suggest that there is selection pressure towards lower migration rate in species that have restricted niche breadth, low resource availability and in species that are on the northern edge of their geographical distribution range. In such species, selection against mobile individuals is likely to result from the decreased probability of finding another habitat patch suitable for egg laying.     

Thermoregulation and microhabitat use in mountain butterflies of the genus Erebia: Importance of fine-scale habitat heterogeneity

Mountain butterflies have evolved efficient thermoregulation strategies enabling their survival in marginal conditions with short flight season and unstable weather. Understanding the importance of their behavioural thermoregulation by habitat use can provide novel information for predicting the fate of alpine Lepidoptera and other insects under ongoing climate change. We studied the link between microhabitat use and thermoregulation in adults of seven species of a butterfly genus Erebia co-occurring in the Austrian Alps. We captured individuals in the field and measured their body temperature in relation to microhabitat and air temperature. We asked whether closely related species regulate their body temperature differently, and if so, what is the effect of behaviour, species traits and individual traits on body to air and body to microhabitat temperature differences. Co-occurring species differed in mean body temperature. These differences were driven by active microhabitat selection by individuals and also by species–specific habitat preferences. Species inhabiting grasslands and rocks utilised warmer microclimates to maintain higher body temperature than woodland species. Under low air temperatures, species of rocky habitats heated up more effectively than species of grasslands and woodlands which allowed them to stay active in colder weather. Species morphology and individual traits play rather minor roles in the thermoregulatory differences; although large species and young individuals maintained higher body temperature. We conclude that diverse microhabitat conditions at small spatial scales probably contribute to sympatric occurrence of closely related species with different thermal demands and that preserving heterogeneous conditions in alpine landscapes might mitigate detrimental consequences of predicted climate change.     

Modelling the distribution of a threatened habitat: the California sage scrub

Aim Using predictive species distribution and ecological niche modelling our objectives are: (1) to identify important climatic drivers of distribution at regional scales of a locally complex and dynamic system – California sage scrub; (2) to map suitable sage scrub habitat in California; and (3) to distinguish between bioclimatic niches of floristic groups within sage scrub to assess the conservation significance of analysing such species groups. Location Coastal mediterranean‐type shrublands of southern and central California. Methods Using point localities from georeferenced herbarium records, we modelled the potential distribution and bioclimatic envelopes of 14 characteristic sage scrub species and three floristic groups (south‐coastal, coastal–interior disjunct and broadly distributed species) based upon current climate conditions. Maxent was used to map climatically suitable habitat, while principal components analysis followed by canonical discriminant analysis were used to distinguish between floristic groups and visualize species and group distributions in multivariate ecological space. Results Geographical distribution patterns of individual species were mirrored in the habitat suitability maps of floristic groups, notably the disjunct distribution of the coastal–interior species. Overlap in the distributions of floristic groups was evident in both geographical and multivariate niche space; however, discriminant analysis confirmed the separability of floristic groups based on bioclimatic variables. Higher performance of floristic group models compared with sage scrub as a whole suggests that groups have differing climate requirements for habitat suitability at regional scales and that breaking sage scrub into floristic groups improves the discrimination between climatically suitable and unsuitable habitat. Main conclusions The finding that presence‐only data and climatic variables can produce useful information on habitat suitability of California sage scrub species and floristic groups at a regional scale has important implications for ongoing efforts of habitat restoration for sage scrub. In addition, modelling at a group level provides important information about the differences in climatic niches within California sage scrub. Finally, the high performance of our floristic group models highlights the potential a community‐level modelling approach holds for investigating plant distribution patterns.