Unveiling the factors shaping the distribution of widely distributed alpine vertebrates, using multi-scale ecological niche modelling of the bat Plecotus macrobullaris
© Alberdi et al.; licensee BioMed Central Ltd. 2014
Received: 5 March 2014
Accepted: 7 October 2014
Published: 29 October 2014
Several alpine vertebrates share a distribution pattern that extends across the South-western Palearctic but is limited to the main mountain massifs. Although they are usually regarded as cold-adapted species, the range of many alpine vertebrates also includes relatively warm areas, suggesting that factors beyond climatic conditions may be driving their distribution. In this work we first recognize the species belonging to the mentioned biogeographic group and, based on the environmental niche analysis of Plecotus macrobullaris, we identify and characterize the environmental factors constraining their ranges. Distribution overlap analysis of 504 European vertebrates was done using the Sorensen Similarity Index, and we identified four birds and one mammal that share the distribution with P. macrobullaris. We generated 135 environmental niche models including different variable combinations and regularization values for P. macrobullaris at two different scales and resolutions. After selecting the best models, we observed that topographic variables outperformed climatic predictors, and the abruptness of the landscape showed better predictive ability than elevation. The best explanatory climatic variable was mean summer temperature, which showed that P. macrobullaris is able to cope with mean temperature ranges spanning up to 16C. The models showed that the distribution of P. macrobullaris is mainly shaped by topographic factors that provide rock-abundant and open-space habitats rather than climatic determinants, and that the species is not a cold-adapted, but rather a cold-tolerant eurithermic organism. P. macrobullaris shares its distribution pattern as well as several ecological features with five other alpine vertebrates, suggesting that the conclusions obtained from this study might be extensible to them. We concluded that rock-dwelling and open-space foraging vertebrates with broad temperature tolerance are the best candidates to show wide alpine distribution in the Western Palearctic.
KeywordsAlpine distribution Alpine long-eared bat Biogeography Chiroptera Distribution Modelling Mountain long-eared bat Zoogeography
Studying the causes shaping distribution patterns is one of the major tasks in biogeography . Many organisms exhibit similar distributions that despite their current ranges could be the result of different historical, ecological and physiological processes. Identifying the common factors constraining their range is key to understanding these geographic analogies. The joint characteristic attribute of alpine species is the geographic restriction to high mountain environments, which were covered with ice-sheets during long periods in the Pleistocene . The study of the biogeographic patterns of Eurasian alpine species has mainly targeted plants and invertebrates - that due to either low or no mobility , are largely conditioned by edaphic and climatic factors . Most vertebrates, however, are unique in their higher mobility, and due to their distinct ecological and evolutionary traits, the biogeographic patterns of vertebrates are likely to differ from those of plants and invertebrates.
The distribution patterns of high mountain vertebrates in the Western Palearctic vary considerably. Endemics restricted to one or few mountain systems are common mainly among vertebrates with low mobility ,. Most volant species however show wider distributions. Many of them are ubiquitous and can be commonly found in other habitats in addition to alpine environments , and few arctic-alpine species show broad boreal distributions with isolated relict populations in alpine areas at lower latitudes . Finally, there is a group of species sharing a broad geographic distribution restricted to Southern Palearctic mountain systems, hereafter referred to as palealpine species, whose distribution pattern has so far received little attention. The distribution of a single bat species, the alpine long-eared bat (Plecotus macrobullaris) is also known to fit this biogeographic pattern . P. macrobullaris is a recently described species - known to forage in alpine meadows  and roost mainly in rock crevices and talus slopes . Unlike other bat species in the Western Palearctic , this species has breeding populations in the subalpine and alpine belts of mountain systems in the South-western Palearctic .
The factors explaining why palealpine vertebrates like P. macrobullaris are restricted to mountain environments, yet absent from the rest of the region, have so far not been addressed. Although species inhabiting alpine environments are commonly regarded as cold-adapted species , this relationship is not always accurate. Whilst species with arctic-alpine distributions are physiologically restricted to cold conditions and remain active in the alpine belt during the winter , this is not necessarily the case for all alpine species . In fact, the current distribution of many palealpine species also includes relatively warm areas in South-eastern Europe and the Middle East , suggesting that factors beyond climatic conditions likely dictate the distribution of these species.
Identifying the factors constraining geographic ranges requires using the correct resolution and extent  because the distribution of species is shaped by processes acting at different scales -. Climatic factors mainly drive the continental scale distribution due to the physiological limits of species . Ecological features, instead, usually act more locally and require finer resolution ,. Additionally, the effect of topographic factors such as slope or aspect also depends on scale . Hence, a multiple resolution and extent approach allows the input of variables in the correct scale in which they act, providing more representative outputs of the environmental factors limiting the distribution of species .
In this study we test whether climatic variables or habitat and topographic factors are the main drivers of the distribution of palealpine vertebrates, so that the main factors constraining their distribution are identified. First, we identify the vertebrate species that exhibit palealpine distribution analysing range similarities. And secondly, in order to understand why palealpine species are closely linked to mountain environments, we apply our modelling approach at two spatial scales using P. macrobullaris as model species. At the broad-scale, across the Western Palearctic, we use coarse-grain presence records of the species for continental-level modelling, while at the fine-scale, around the Pyrenees, we use precise roosting records and finer spatial resolution to understand the factors shaping the local distribution.
Broad-scale distribution modelling
Composition, evaluation scores (AUCtest, AICc and MPA) and variable contribution ranks of the best models
Variable contribution rank
B8, B10, B12, B17
B10 (45.3) > B12 > B8 > B17
B4, B10, B12, B15, B17
B10 (43.7) > B12 > B17 > B15 > B4
ABR (84.4) > ELEV
ABR (86.1) > ELEV
Climatic + Topographic
B4, B8, B10, B12, B15, B17, ELEV, ABR
ABR (71.7) > ELEV > B10 > B17 > B12 > B15 > B4 > B8
B4, B8, B10, B12, B15, B17, ABR
ABR (79.9) > B10 > B12 > B15 > B17 > B4 > B8
Climatic + Habitat
B4, B8, B10, B12, B15, B17, LAND
B10 (32.7) > B12 > LAND > B8 > B17 > B15 > B4
B4, B10, B12, B15, B17, LAND
B10 (37.0) > B12 > LAND > B17 > B15 > B4
Topographic + Habitat
ELEV, ABR, LAND
ABR (86.7) > ELEV > LAND
ABR (96.6) > LAND
Climatic + Topographic + Habitat
B4, B8, B10, B12, B15, B17, ABR, LAND
ABR (70.6) > B10 > B17 > LAND > B15 > B12 > B8 > B4
B4, B10, B12, B15, B17, ABR, LAND
ABR (75.1) > B10 > BIO12 > B17 > B15 > LAND > B4
Variable contribution rank
ELEV, SLOP, ORI
ELEV (56.3) > SLO > ORI
ELEV (54.9) > SLO
DIS-URBAN, DIS-FOREST, DIS-ROCK
dis-rock (74.3) > dis-urban > dis-forest
dis-rock (89.3) > dis-urban
Topographic + Habitat
ELEV, SLOP, ORI, DIS-URBAN, DIS-FOREST, DIS-ROCK
dis-rock (32.0) > SLO > ELEV > dis-urban > dis-forest > ORI
ELEV, SLOP, DIS-URBAN, DIS-ROCK
dis-rock (34.4) > SLO > ELEV > dis-urban
Different metrics on the contribution of variables to the best models
Jackknife training gain with only variable
Jackknife training gain without this variable
Fine-scale roosting habitat suitability modelling
The best model was built using two habitat (distance to rocky and urban areas) and two topographic variables (elevation and slope) (Table1, Additional file 1: Table S2). The ENM showed good predictive ability (AUC =0.956). The three main factors shaping roosting habitat suitability for P. macrobullaris were the distance to rock areas (DIS-ROCK), slope and elevation. DIS-ROCK was the variable with the highest relative contribution to the best model, and Jackknife tests stressed the importance of slope and elevation (Table2, Additional file 2: Figure S3). The closer the rocks and higher the slope the greater the probability of finding suitable roosting sites (Figure3C and E). The response curve for elevation showed suitable areas for roosting in an elevation range from 1300 to 2400m, with maximum values between 1500 and 2000m (Figure3D).
The generated broad-scale model predicted that suitable areas for P. macrobullaris are found in almost all mountain ranges in the Western Palearctic, including some ranges where the species has not being hitherto reported. The species might be present in the Atlas Mts. in Morocco as well as in several mountain chains in the Iberian and Italian Peninsulas, as occurs with some other palealpine species (Figure4). Two of the palealpine birds inhabit the Atlas while the Cantabrian Mts. and the Apennines host all the mentioned species. On the other hand, the presence of the species in the suitable mountain ranges of Scotland and Scandinavia seems implausible. Both ranges were covered by ice caps during Pleistocene glaciations , and the large unsuitable flatland separating them from the southern mountain ranges could have prevented their postglacial colonisation by palealpine species. Additionally, harder climatic conditions in winter than in Southern Europe might be another reason that hampers colonisation. This view is supported by the distribution of other palealpine species, which are absent from both Scotland and Scandinavia despite the relatively high mobility of most of them (Figure4).
Our results differ from the previously published ecological niche model of P. macrobullaris in Switzerland , which predicted a distribution restricted to lower elevations than what our model did for the Alps. Nevertheless, both research are hardly comparable as they only used localities of colonies detected in buildings at valley bottoms for generating their ENM, which might led to a biased model of the distribution of P. macrobullaris, since the species also uses natural roosting resources , and it has been also reported above the treeline in the Alps .
The role of topographic variables
Abruptness likely acts as a surrogate of several ecological variables. One of them is high rock availability, an important resource for palealpine species. In fact, the fine-scale modelling disclosed the importance of the proximity of rocky environments and the steepness of the terrain, surfacing the significance of rock walls and their debris. The five vertebrates identified as palealpine species depend on rocks for sheltering -. Cliffs and talus slopes are also a key roosting resource for P. macrobullaris, even though it also uses buildings for sheltering ,,. The relative area of rock outcrops is extensive at high mountain environments, favouring roost availability for saxicolous species. This behaviour probably allowed palealpine species to expand their ranges throughout the narrow and steep mountain massifs in the Western Palearctic. For instance, most snowfinches (genera Montifringilla, Pyrgilauda and Onychostruthus) roost on the ground and are behaviourally adapted to living in steppes , but only the saxicolous M. nivalis has expanded its distribution towards the European mountain massifs.
Another characteristic shared by all palealpine species is their preference for foraging in open areas, mainly in habitats with short grassy vegetation and rocks ,-. Climatic factors hamper the development of wood patches in alpine areas, providing suitable grounds for open-space foragers. Steep environments at lower elevation can also provide similar open habitats aided by the outcrop of bedrock and the decrease in soil thickness.
Finally, the high elevation range of mountain areas allows species to perform short-distance altitudinal migrations. Boreal species cannot avoid adverse atmospheric conditions by performing short migrations, and instead need to be physiologically adapted to extreme cold temperatures. In contrast, palealpine species tend to commute to lower elevations when conditions become harsher ,,, and therefore physiological adaptations to the extreme cold may not be a compulsory requirement for them. Despite altitudinal migration have been documented in bats , such a behaviour has not been studied in P. macrobullaris so far. Nevertheless, it seems implausible bats could remain the whole year in alpine environments, and ecological features as well as comparative biogeography suggest that P. macrobullaris may actually perform seasonal altitudinal migrations .
The role of climatic variables
Despite the lower contribution of climatic variables compared to topographic factors, the former also influenced the broad-scale model. Temperature mainly encompasses the altitudinal domains of the species distribution. Different climatic patterns along the elevation gradients of each mountain range (with differences of up to 13C in mean summer temperature at the same elevation; Additional file 2: Figure S4) produced contrasting suitability-elevation relationships between mountain chains across the predicted range of P. macrobullaris (Figure5). For instance, suitable areas in the Zagros Mountains (Iran) were restricted to areas above 1100m, whereas near-sea-level areas of the Alps and the Dinaric Alps were considered suitable. This fact explains the poor predictive ability of elevation in the continental scale models. The fine-scale model also reflected the limiting effect of temperature in the Pyrenees, since no areas above 2600m were considered suitable despite containing abundant steep areas.
Broad-scale modelling showed that P. macrobullaris, and by extension palealpine species, are able to cope with mean summer temperatures that span from 6 to 22C (Figure3B), revealing their eurithermic nature. Even though palealpine vertebrates were previously considered as cold-adapted species , their distribution is not restricted to cold environments, suggesting they are cold-tolerant rather than cold-adapted. Climatic models have predicted that alpine species like P. graculus and M. nivalis might extend their distribution to high latitude areas due to climate change . Even though those areas might become climatically suitable for alpine species, our results suggest that such a scenario is highly improvable. We show that topography and its correlating ecological factors are shaping mostly the distribution of such species, and thus, the large unsuitable flatland between boreal areas and Southern European mountain ranges would hamper putative colonisation events. The wide thermal tolerance likely facilitated the colonisation of most of the mountain massifs in the Central and Southern Palearctic. Conversely, the southern isolated populations of an arctic-alpine vertebrate such as the rock ptarmigan (Lagopus muta) are limited to mountain ranges that hold extensive areas with summer mean temperatures below 10C, i.e. the Pyrenees and the Alps in Europe .
The human factor
In addition to topographic and climatic factors, human activities probably affected the distribution of palealpine vertebrates. Similar to other palealpine species ,,, P. macrobullaris is also an open-space forager . The natural treeline in the Pyrenees is expected to be around 20002200m; however, historical landscape management for pasture has lowered the upper forest limit to around 1500m . It is noteworthy that the most suitable elevation range predicted by both models for the Pyrenees coincides with this modified belt (1500 2000m), which suggests that it may offer the best ecological conditions for P. macrobullaris and other palealpine vertebrates in terms of greater food abundance and diversity  and milder climatic conditions than at higher elevations . Hence, future changes in landscape management due to pastoral abandonment will likely affect the foraging ecology of P. macrobullaris and other alpine vertebrates .
The importance of biotic interactions in the distribution of species is a largely recognised phenomenon . Processes like competitive exclusion  or predatorprey relationships  can affect the distribution of species, in addition to climatic and habitat-related factors. Accordingly, Rutishauser et al.  reported that the mostly allopatric distribution of P. macrobullaris and P. austriacus in Switzerland might be the result of competitive exclusion. In fact, both species share ecological traits ,,. However, under this scenario sympatric populations would be the exception rather than the rule and the niche of P. macrobullaris would expand where it is released from interspecific competition. Neither one nor the other seems to occur. There are broad areas where both species are found in sympatry, such as the Pyrenees (own data), the French Prealps  and the Dinaric Alps , and the ecological niche of P. macrobullaris as well as its link to mountain environments is kept beyond the range of P. austriacus, e.g. the eastern part of the distribution of P. macrobullaris. Competitive exclusion might occur in some specific areas like Switzerland, or across the altitudinal gradient in precise contact zones . We do not discard that competitive exclusion shows up where resources are limiting, but in our opinion there are no signs indicating a general effect on the ecological niche and continental-level distribution of P. macrobullaris.
The environmental niche of palealpine species
Our results show that climatic factors are not the main drivers of the distribution of palealpine vertebrates. Generally, the broad-scale geographic distribution of species tend to be limited by climatic variables , while non-climatic factors usually shape their distribution at the regional or local level . However, our models show that in the case of palealpine species topographic factors act at a broader scale. Unlike arctic-alpine organisms, palealpine species are able to thrive in relatively warm as well as cold temperatures, and this eurithermicity is responsible for the limited ability of climatic variables to predict their distribution at the broad-scale.
Although palealpine species share the same biogeographic pattern driven by topographic factors, they are expected to show differential responses to climatic conditions. Variations in behaviour as well as geographical and altitudinal distribution suggest that some species cope better with cold than others, e.g. the white-winged snowfinch compared to the wallcreeper . Differences in thermal tolerance between groups of organisms are also expected to occur due to differences in their ecology, reproduction and physiology. Major differences between birds and bats probably stem from their distinct circadian rhythm, thermoregulation abilities and reproductive systems. Birds are diurnal, homeothermic and oviparous animals , while bats are highly heterothermic nocturnal organisms with viviparous reproduction . Even though facultative heterothermy may offer an advantage over homeothermic organisms when exploiting alpine environments, this advantage is only applicable to non-reproductive bats. Breeding females cannot decrease their body temperature as much as non-breeding bats due to the metabolic demands of foetal development and lactation , so they may be hindered from exploiting high mountain habitats. Consequently, elevational segregation might occur in the coldest regions of their distribution, which would explain the high number of P. macrobullaris maternity colonies found below the treeline in the Alps ,. Nevertheless, such a pattern was not so noticeable in the slightly milder range of the Pyrenees . Finally, models generated in this study reflect the environmental niche during the favourable season. During winter, palealpine vertebrates behave differently depending on species and taxa: bats hibernate while birds migrate to lower elevations.
The use of multi-scale modelling allowed us to identify features shaping the biogeographic pattern shared by alpine vertebrates with wide geographic distribution, which could pave the way for future studies and help fine-tune conservation measures for this group of animals. The wide temperature tolerance identified in this study suggests that rather than physiological limitations relating to warmer conditions, ecological factors, such as changes in treeline elevation and the incursion of lowland species, may be the main challenges palealpine species will have to cope with under future conditions resulting from climate change . Palealpine species are found in high-mountain environments not because they are constrained by cold temperatures, but because these cold-tolerant species find in alpine environments their preferred foraging and sheltering habitats, characterized by abundant open-space and rocky areas.
We calculated the distribution area overlap of P. macrobullaris and 503 European terrestrial vertebrates using the Srensen Similarity Index (SSI - ) in order to identify the species sharing the distribution pattern with P. macrobullaris. We generated a reference GIS layer based on updated information about the distribution of P. macrobullaris and compared it with breading season distribution layers of 203 mammals available at IUCN Red List  and 298 birds obtained from BirdLife International . We also compared the distribution area of P. macrobullaris we generated with the area available at IUCN Red List to assess the uncertainty of the data and create a corrected SSI (cSSI) that reflects better the distribution similarities. Visual inspection of the geographical distributions of the species with the highest cSSI values showed that cSSI = 0.60 marks a resemblance threshold between widely distributed species limited to mountain ranges and species with a range that extends beyond the main mountain chains.
Environmental Niche Modelling (ENM)
The extent of the broad-scale models was set as the Western Palearctic region (20–70° N, −20−60° E), while fine-scale models were limited to the Pyrenees and the surrounding area (40–42.1° N, −1.1−2.2° E). All predictive models were generated with the presence-only species distribution modelling software Maxent (Maximum Entropy Algorithm) version 3.3.3 ,. When running the software, we set maximum iterations to 5000, used 10000 background points, and ran 50 replicates of each model, following a resampling method randomly selecting a subsample of 25% of location records for model validation in each round. The output of all models was set in raw format in order to analyse them using the software ENMTools. Correlation between variables was analysed using ENMTools v1.3 , and the variable with a priori the highest ecological meaning was selected when the Spearmans correlation coefficient (ρ) between two variables was greater than 0.65. The remaining spatial and statistical analyses were performed using ArcView GIS 3.2 and R 2.9.2 software (packages raster and dismo) available from CRAN .
We used presence location records from all of the known distribution area of P. macrobullaris. In order to minimise errors and ensure the best model performance, we took several measures following Merow et al. . We applied 5 filters to all the available records with the following criteria: (1) Molecular verification: species identification was checked by molecular means to avoid misidentifications with morphologically similar congeners; (2) Resolution: maximum location uncertainty of 1km; (3) Pseudoreplication: a minimum distance of 10km between presence records was set based on the home-range information of the species ,; (4) Record source-type homogeneity: similar amount of different record types (e.g. roost observations, net captures); and (5) Geographical homogeneity: similar amount of records in each geographic area. Sampling effort East of the Bosphorus has been considerably lower than in Europe, and the number of European records were four-fold the number of eastern locations after applying the filtering. In order to correct the sampling bias , a bias grid was introduced to the algorithm where the relative sampling probability West of the Bosphorus was considered four times higher than the eastern area. This helped outweigh the relative value of Anatolia, the Near East and the Caucasus compared to Europe. Spatial autocorrelation analyses were discarded because the disjunct distribution of alpine species is actually grouped in clusters limited to high mountain massifs.
Characteristics of the variables used in the two modelling scales
Generated from elevation data
Worldclim Bioclim 4
Mean temperature of the wettest quarter
Worldclim Bioclim 8
Mean temperature of the warmest quarter
Worldclim Bioclim 10
Worldclim Bioclim 12
Worldclim Bioclim 15
Precipitation of the Driest Quarter
Worldclim Bioclim 17
SRTM 90m DEM (CGIAR-CSI)
Generated from elevation data
Generated from elevation data
Distance to rock areas
Obtained from Corine LandCover 2006
Distance to forest areas
Obtained from Corine LandCover 2006
Distance to urban areas
Obtained from Corine LandCover 2006
We used roost location data for modelling the fine-scale environmental niche of P. macrobullaris (Additional file 2: Figure S1). Roosts were identified using the homing-in radio-tracking method . Only roosts located in natural structures were used, as they represent the original as well as the predominant roost type of P. macrobullaris. In order to minimise pseudoreplication, one roost for each grid cell was used for the analysis, resulting in a total of 43 records. We used 6 variables (3 topographic and 3 habitat-related Table3) to generate 45 different models. Climatic variables were not used in the fine-scale modelling for being highly correlated with elevation. All variables were edited and homogenised to 100m resolution. Elevation data was obtained from CGIAR-CSI , and the rest of topographic variables were derived from it. Habitat variables were obtained from Corine Land Cover 2006  and categories were reclassified in ArcView GIS 3.2.
Abruptness was used in the broad-scale modelling for reflecting the elevation range available for the bats in each geographic location, and it was based on the home-range information of P. macrobullaris,. Conversely, we used slope in the fine-scale modelling because we considered it might play a key role providing suitable rock roosts to P. macrobullaris.
Model evaluation and selection
Models were evaluated using both threshold-independent and dependent means that deal with several aspects of model performance, including model accuracy, model complexity and prediction success -. Accuracy was evaluated using the Area Under the Curve (AUC) of the Receiver Operator Characteristics (ROC) ,-. Corrected Akaikes Information Criteria (AICc)  was used for evaluating model complexity in ENMTools. For that task a single model out of the 50 replicates was selected using a Principal Components Analysis (PCA), in order to identify the replicate with the closest values to the average. The PCA was performed using five matrices: regularised training gain, test gain, test AUC, 10 percentile training presence logistic threshold and 10 percentile training presence area. We also used the threshold-dependent minimal predicted area (MPA) evaluator . Finally, we validated the best broad-scale model using all published distribution records (349) of P. macrobullaris.
Model output and variable importance
The best models were repeated and plotted in a logistic format to provide the estimates of the probability of occurrence as predicted by the variables in a map ,. We used the threshold selection method max SSS to discern suitable and unsuitable areas . The relative importance of each variable was checked using heuristic (percent contribution), permutation (permutation importance) and jackknife approaches (Table2) ,. Response curves for the best explanatory variables were plotted in order to determine the response of suitability values to changes in specific predictors, and thus identify the most suitable conditions for the species (Figure3).
AA, JA and IG designed the study. All authors participated in the fieldwork. AA performed the analyses. All authors drafted the manuscript. All authors have read and approved the final manuscript.
We wish to thank Orly Razgour for her insightful comments on the manuscript and for proofreading the text. The Basque Government (research projects IT385-07, IT301-10 and IT754-13, and predoctoral grants BFI-2010-190 and BFI-2009-252) and the University of the Basque Country UPV/EHU supported this work.
- Lomolino MV, Riddle BR, Whittaker RJ, Brown JH: Biogeography. 2010, Sinauer Associates, SunderlandGoogle Scholar
- Hughes PD, Woodward JC, van Calsteren PC, Thomas LE, Adamson KR: Pleistocene ice caps on the coastal mountains of the Adriatic Sea. Quat Sci Rev. 2010, 29: 3690-3708. 10.1016/j.quascirev.2010.06.032.View ArticleGoogle Scholar
- Schmitt T, Muster C, Schnswetter P: Are Disjunct Alpine and Arctic-Alpine Animal and Plant Species in the Western Palearctic Really Relics of a Cold Past?. Relict Species Phylogeography and Conservation Biology. Edited by: Habel JC, Assmann T. 2009, Springer, Berlin, 239-252.Google Scholar
- Ehrich D, Gaudeul M, Assefa A, Koch MA, Mummenhoff K, Nemomissa S, Consortium I, Brochmann C: Genetic consequences of Pleistocene range shifts: contrast between the Arctic, the Alps and the East African mountains. Mol Ecol. 2007, 16: 2542-2559. 10.1111/j.1365-294X.2007.03299.x.PubMedView ArticleGoogle Scholar
- Varga ZS, Schmitt T: Types of oreal and oreotundral disjunctions in the western Palearctic. Biol J Linn Soc. 2008, 93: 415-430. 10.1111/j.1095-8312.2007.00934.x.View ArticleGoogle Scholar
- Ozenda P: Die Vegetation Der Alpen Im Europischen Gebirgsraum. 1988, Fischer, StuttgartGoogle Scholar
- Roland J, Keyghobadi N, Fownes S: Alpine parnassius butterfly dispersal: effects of landscape and population size. Ecology. 2000, 81: 1642-1653. 10.1890/0012-9658(2000)081[1642:APBDEO]2.0.CO;2.View ArticleGoogle Scholar
- Finn DS, Theobald DM, Black WC, Poff NL: Spatial population genetic structure and limited dispersal in a Rocky Mountain alpine stream insect. Mol Ecol. 2006, 15: 3553-3566. 10.1111/j.1365-294X.2006.03034.x.PubMedView ArticleGoogle Scholar
- Dubuis A, Giovanettina S, Pellissier L, Pottier J, Vittoz P, Guisan A: Improving the prediction of plant species distribution and community composition by adding edaphic to topo-climatic variables. J Veg Sci. 2012, 24: 593-606. 10.1111/jvs.12002.View ArticleGoogle Scholar
- Carranza S, Amat F: Taxonomy, biogeography and evolution of Euproctus (Amphibia: Salamandridae), with the resurrection of the genus Calotriton and the description of a new endemic species from the Iberian Peninsula. Zool J Linnean Soc. 2005, 145: 555-582. 10.1111/j.1096-3642.2005.00197.x.View ArticleGoogle Scholar
- Krytufek B, BUAN EV: Rarity and decline in palaeoendemic Martinos vole Dinaromys bogdanovi. Mammal Rev. 2008, 38: 267-284. 10.1111/j.1365-2907.2008.00127.x.View ArticleGoogle Scholar
- Meffert PJ, Marzluff JM, Dziock F: Unintentional habitats: value of a city for the wheatear (Oenanthe oenanthe). Landscape Urban Plan. 2012, 108: 49-56. 10.1016/j.landurbplan.2012.07.013.View ArticleGoogle Scholar
- Favaron M, Scherini GC, Preatoni D, Tosi G, Wauters LA: Spacing behaviour and habitat use of rock ptarmigan (Lagopus mutus) at low density in the Italian Alps. J Ornithol. 2006, 147: 618-628. 10.1007/s10336-006-0087-z.View ArticleGoogle Scholar
- Alberdi A, Garin I, Aizpurua O, Aihartza J: Review on the geographic and elevational distribution of the mountain long-eared bat Plecotus macrobullaris, completed by utilising a specific mist-netting technique. Acta Chiropterol. 2013, 15: 451-461. 10.3161/150811013X679071.View ArticleGoogle Scholar
- Kiefer A, Veith M: A new species of long-eared bat from Europe (Chiroptera: Vespertilionidae). Myotis. 2002, 39: 5-16.Google Scholar
- Spitzenberger F, Haring E, Tvrtkovic N: Plecotus microdontus (Mammalia, Vespertilionidae), a new bat species from Austria. Natura Croatica. 2002, 11: 1-18.Google Scholar
- Spitzenberger F, Strelkov P, Haring E: Morphology and mitochondrial DNA sequences show that Plecotus alpinus Kiefer & Veith, 2002 and Plecotus microdontus Spitzenberger, 2002 are synonyms of Plecotus macrobullaris Kuzjakin, 1965. Natura Croatica. 2003, 12: 39-53.Google Scholar
- Alberdi A, Garin I, Aizpurua O, Aihartza J: The foraging ecology of the mountain long-eared Bat Plecotus macrobullaris revealed with DNA mini-barcodes. PLoS One. 2012, 7: e35692-10.1371/journal.pone.0035692.PubMedPubMed CentralView ArticleGoogle Scholar
- Alberdi A, Aihartza J, Aizpurua O, Salsamendi E, Brigham M, Garin I: Living above the treeline: roosting ecology of the alpine batPlecotus macrobullaris.Eur J Wildl Res 2014. doi:10.1007/s10344-014-0862-8.,Google Scholar
- Dietz C, Dietmar N, von Helversen O: Bats of Britain, Europe and Northwest Africa. 2009, A&C Black, LondonGoogle Scholar
- Smith SE, Gregory RD, Anderson BJ, Thomas CD: The past, present and potential future distributions of cold-adapted bird species. Divers Distrib. 2013, 19: 352-362. 10.1111/ddi.12025.View ArticleGoogle Scholar
- Bienkowski P, Marszalek U: Metabolism and energy budget in the snow vole. Acta Theriol. 1974, 19: 55-67. 10.4098/AT.arch.74-4.View ArticleGoogle Scholar
- Arlott N: Collins Field Guide: Birds of the Palearctic: Passerines. 2007, Collins, LondonGoogle Scholar
- Boyce M: Scale for resource selection functions. Divers Distrib. 2006, 12: 269-276. 10.1111/j.1366-9516.2006.00243.x.View ArticleGoogle Scholar
- Levin SA: The problem of pattern and scale in ecology. Ecology. 1992, 73: 1943-1967. 10.2307/1941447.View ArticleGoogle Scholar
- Hortal J, Roura-Pascual N, Sanders NJ, Rahbek C: Understanding (insect) species distributions across spatial scales. Ecography. 2010, 33: 51-53. 10.1111/j.1600-0587.2009.06428.x.View ArticleGoogle Scholar
- Lomba A, Pellissier L, Randin C, Vicente J, Moreira F, Honrado J, Guisan A: Overcoming the rare species modelling paradox: a novel hierarchical framework applied to an Iberian endemic plant. Biol Cons. 2010, 143: 2647-2657. 10.1016/j.biocon.2010.07.007.View ArticleGoogle Scholar
- Gaston KJ: The Structure and Dynamics of Geographic Ranges. 2003, Oxford University Press on Demand, OxfordGoogle Scholar
- Pearson RG, Dawson TP, Liu C: Modelling species distributions in Britain: a hierarchical integration of climate and land-cover data. Ecography. 2004, 27: 285-298. 10.1111/j.0906-7590.2004.03740.x.View ArticleGoogle Scholar
- Graf RF, Bollmann K, Suter W, Bugmann H: The importance of spatial scale in habitat models: capercaillie in the Swiss Alps. Landscape Ecol. 2005, 20: 703-717. 10.1007/s10980-005-0063-7.View ArticleGoogle Scholar
- Piedallu C, Ggout J-C: Efficient assessment of topographic solar radiation to improve plant distribution models. Agric For Meteorol. 2008, 148: 1696-1706. 10.1016/j.agrformet.2008.06.001.View ArticleGoogle Scholar
- Razgour O, Hanmer J, Jones G: Using multi-scale modelling to predict habitat suitability for species of conservation concern: The grey long-eared bat as a case study. Biol Cons. 2011, 144: 2922-2930. 10.1016/j.biocon.2011.08.010.View ArticleGoogle Scholar
- Mangerud J, Jakobsson M, Alexanderson H, Astakhov V, Clarke GKC, Henriksen M, Hjort C, Krinner G, Lunkka J-P, Mller P, Murray A, Nikoloskaya O, Saarnisto M, Svendsen JI: Ice-dammed lakes and rerouting of the drainage of northern Eurasia during the last glaciation. Quat Sci Rev. 2004, 23: 1313-1332. 10.1016/j.quascirev.2003.12.009.View ArticleGoogle Scholar
- Rutishauser MD, Bontadina F, Braunisch V, Ashrafi S, Arlettaz R: The challenge posed by newly discovered cryptic species: disentangling the environmental niches of long-eared bats. Divers Distrib. 2012, 18: 1107-1119. 10.1111/j.1472-4642.2012.00904.x.View ArticleGoogle Scholar
- Hof AR, Jansson R, Nilsson C: The usefulness of elevation as a predictor variable in species distribution modelling. Ecol Model. 2012, 246: 86-90. 10.1016/j.ecolmodel.2012.07.028.View ArticleGoogle Scholar
- Hrsch B: Modelling the spatial distribution of montane and subalpine forests in the central Alps using digital elevation models. Ecol Model. 2003, 168: 267-282. 10.1016/S0304-3800(03)00141-8.View ArticleGoogle Scholar
- Newbold T, Gilbert F, Zalat S, El-Gabbas A, Reader T: Climate-based models of spatial patterns of species richness in Egypts butterfly and mammal fauna. J Biogeogr. 2009, 36: 2085-2095. 10.1111/j.1365-2699.2009.02140.x.View ArticleGoogle Scholar
- Yan-Hua Q, Zuo-Hua Y, De Ritis SF: Distribution patterns of snow finches (genus Montifringilla) in the Tibetan Plateau of China. Avocetta. 2002, 26: 11-18.Google Scholar
- Saniga M: Seasonal distribution, altitudinal and horizontal migration of the wallcreeper (Tichodroma muraria) in the Mal Fatra Mountais, Slovak Carpathians. Folia Zool. 1995, 44: 237-246.Google Scholar
- Luque-Larena JJ, Lpez P, Goslbez J: Microhabitat use by the snow vole Chionomys nivalis in alpine environments reflects rock-dwelling preferences. Can J Zool. 2002, 80: 36-41. 10.1139/z01-197.View ArticleGoogle Scholar
- Mattei-Roesli M: Situazione del genere Plecotus (Chiroptera) nel Cantone Ticino (Svizzera). Boll Soc Ticin Sci Nat. 2010, 98: 31-34.Google Scholar
- Presetnik P, Koselj K, Zagmajster M: Atlas Netopirjev (Chiroptera) Slovenije. Atlas of Bats (Chiroptera) of Slovenia. 2009, Miklav na Dravskem polju: Center za kartografijo favne in flore, LjubianaGoogle Scholar
- Rolando A, Laiolo P, Formica M: A comparative analysis of the foraging behaviour of the Chough Pyrrhocorax pyrrhocorax and the Alpine Chough Pyrrhocorax graculus coexisting in the Alps. Ibis. 1997, 139: 461-467. 10.1111/j.1474-919X.1997.tb04659.x.View ArticleGoogle Scholar
- Laiolo P, Dondero F, Ciliento E, Rolando A: Consequences of pastoral abandonment for the structure and diversity of the alpine avifauna. J Appl Ecol. 2004, 41: 294-304. 10.1111/j.0021-8901.2004.00893.x.View ArticleGoogle Scholar
- Prez-Aranda D: Biologa, ecologa, gentica y conservacin del topillo nival (Chionomys nivalis) en Pealara y Sierra Nevada. In PhD Thesis. Madrid: Complutense University, Department of Zoology and Physical Anthropology; 2009.,Google Scholar
- Henry PY: Differential migration in the polygynandrous Alpine Accentor Prunella collaris. Bird Study. 2011, 58: 160-170. 10.1080/00063657.2010.543948.View ArticleGoogle Scholar
- Laiolo P, Rolando A, Carisio L: Winter movements of the alpine chough: implications for management in the Alps. J Mount Ecol. 2001, 6: 21-30.Google Scholar
- McGuire LP, Boyle WA: Altitudinal migration in bats: evidence, patterns, and drivers. Biol Rev. 2013, 88: 767-786.PubMedView ArticleGoogle Scholar
- Bech N, Boissier J, Drovetski S, Novoa C: Population genetic structure of rock ptarmigan in the sky islands of French Pyrenees: implications for conservation. Anim Conserv. 2009, 12: 138-146. 10.1111/j.1469-1795.2008.00233.x.View ArticleGoogle Scholar
- Janeau G, Aulagnier S: Snow vole - Chionomys nivalis. Ibex JME. 1997, 4: 1-11.Google Scholar
- Janiga M, Novotn M: Feeding preferences and foraging behaviour in the Alpine Accentor Prunella collaris. Ornis Fennica. 2006, 83: 170-180.Google Scholar
- Ninot JM, Carrillo E, Font X, Carreras J, Ferr A, Masalles RM, Soriano I, Vigo J: Altitude zonation in the Pyrenees. A geobotanic interpretation. Phytocoenologia. 2007, 37: 371-398. 10.1127/0340-269X/2007/0037-0371.View ArticleGoogle Scholar
- Beck J, Altermatt F, Hagmann R, Lang S: Seasonality in the altitudediversity pattern of Alpine moths. Basic Appl Ecol. 2010, 11: 714-722. 10.1016/j.baae.2010.08.009.View ArticleGoogle Scholar
- Barry RG: Mountain Weather and Climate. 2008, Cambridge University Press, New YorkView ArticleGoogle Scholar
- Arajo MB, Luoto M: The importance of biotic interactions for modelling species distributions under climate change. Global Ecol Biogeogr. 2007, 16: 743-753. 10.1111/j.1466-8238.2007.00359.x.View ArticleGoogle Scholar
- Wisz MS, Pottier J, Kissling WD, Pellissier L, Lenoir J, Damgaard CF, Dormann CF, Forchhammer MC, Grytnes J-A, Guisan A, Heikkinen RK, Hye TT, Khn I, Luoto M, Maiorano L, Nilsson M-C, Normand S, ckinger E, Schmidt NM, Termansen M, Timmermann A, Wardle DA, Aastrup P, Svenning J-C: The role of biotic interactions in shaping distributions and realised assemblages of species: implications for species distribution modelling. Biol Rev. 2012, 88: 15-30. 10.1111/j.1469-185X.2012.00235.x.PubMedPubMed CentralView ArticleGoogle Scholar
- de Arajo CB, Marcondes-Machado LO, Costa GC: The importance of biotic interactions in species distribution models: a test of the Eltonian noise hypothesis using parrots. J Biogeogr. 2013, 41: 513-523. 10.1111/jbi.12234.View ArticleGoogle Scholar
- Dietrich S, Szameitat DP, Kiefer A, Schnitzler HU, Denzinger A: Echolocation signals of the plecotine bat, Plecotus macrobullaris Kuzyakin, 1965. Acta Chiropterol. 2006, 8: 465-475. 10.3161/1733-5329(2006)8[465:ESOTPB]2.0.CO;2.View ArticleGoogle Scholar
- Razgour O, Clare EL, Zeale M, Hanmer J, Schnell I, Rasmussen M, GIlbert T, Jones H: High-throughput sequencing offers insight into mechanisms of resource partitioning in cryptic bat species. Ecol Evol. 2011, 1 (4): 556-570. 10.1002/ece3.49.PubMedPubMed CentralView ArticleGoogle Scholar
- Arthur L, Lemaire M: Les Chauves-Souris De France, Belgique, Luxembourg Et Suisse. 2009, Biotope Editions, MezeGoogle Scholar
- Tvrtkovic N, Pavlinic I, Haring E: Four species of long-eared bats (Plecotus, Geoffroy, 1818; Mammalia, Vespertilionidae) in Croatia: field identification and distribution. Folia Zool. 2005, 54: 75-88.Google Scholar
- Pavlinic I, Tvrtkovic N: Altitudinal distribution of four Plecotus species (Mammalia, Vespertilionidae) occurring in Croatia. Croatian Natural History Museum. 2004, 13: 395-401.Google Scholar
- Gaston KJ: Rarity (Population and Community Biology Series). 1994, Springer, Falmouth, UKGoogle Scholar
- Bicudo JE, Buttemer WA, Chappell MA, Pearson JT, Bech C: Ecological and Environmental Physiology of Birds. 2010, Oxford University Press, New YorkView ArticleGoogle Scholar
- Kunz TH, Fenton MB: Bat Ecology. 2006, University of Chicago Press, ChicagoGoogle Scholar
- Dzal YA, Brigham RM: The tradeoff between torpor use and reproduction in little brown bats (Myotis lucifugus). J Comp Physiol B. 2012, 183: 279-288. 10.1007/s00360-012-0705-4.PubMedView ArticleGoogle Scholar
- Reif J, Flousek J: The role of species ecological traits in climatically driven altitudinal range shifts of central European birds. Oikos. 2012, 121: 1053-1060. 10.1111/j.1600-0706.2011.20008.x.View ArticleGoogle Scholar
- Srensen T: A method of establishing groups of equal amplitude in plant sociology based on similarity of species and its application to analyses of the vegetation on Danish commons. Biol Skr. 1948, 5: 1-34.Google Scholar
- IUCN Red List of Threatened Species. , [http://www.iucnredlist.org]
- BirdLife International. , [http://www.birdlife.org]
- Phillips SJ, Anderson RP, Schapire RE: Maximum entropy modeling of species geographic distributions. Ecol Model. 2006, 190: 231-259. 10.1016/j.ecolmodel.2005.03.026.View ArticleGoogle Scholar
- Phillips S, Dudk M: Modeling of species distributions with Maxent: new extensions and a comprehensive evaluation. Ecography. 2008, 31: 161-175. 10.1111/j.0906-7590.2008.5203.x.View ArticleGoogle Scholar
- Warren D, Glor R, Turelli M: ENMTools: a toolbox for comparative studies of environmental niche models. Ecography. 2010, 33: 607-611. 10.1111/j.1600-0587.2009.06041.x.View ArticleGoogle Scholar
- The Comprehensive R Archive Network. , [http://cran.r-project.org/]
- Merow C, Smith MJ, Silander JA: A practical guide to MaxEnt for modeling species distributions: what it does, and why inputs and settings matter. Ecography. 2013, 36: 1058-1069. 10.1111/j.1600-0587.2013.07872.x.View ArticleGoogle Scholar
- Preatoni D, Spada M, Wauters L, Tosi G, Martinoli A: Habitat use in the female Alpine long-eared bat (Plecotus macrobullaris): does breeding make the difference?. Acta Chiropterol. 2011, 13: 355-364. 10.3161/150811011X624820.View ArticleGoogle Scholar
- Kramer-Schadt S, Niedballa J, Pilgrim JD, Schrder B, Lindenborn J, Reinfelder V, Stillfried M, Heckmann I, Scharf AK, Augeri DM, Cheyne SM, Hearn AJ, Ross J, Macdonald DW, Mathai J, Eaton J, Marshall AJ, Semiadi G, Rustam R, Bernard H, Alfred R, Samejima H, Duckworth JW, Breitenmoser-Wuersten C, Belant JL, Hofer H, Wilting A: The importance of correcting for sampling bias in MaxEnt species distribution models. Divers Distrib. 2013, 19: 1366-1379. 10.1111/ddi.12096.View ArticleGoogle Scholar
- Lobo JM, Jimnez-Valverde A, Real R: AUC: a misleading measure of the performance of predictive distributingon models. Global Ecol Biogeogr. 2008, 17: 145-151. 10.1111/j.1466-8238.2007.00358.x.View ArticleGoogle Scholar
- VanDerWal J, Shoo L, Graham CH, Williams S: Selecting pseudo-absence data for presence-only distribution modeling: How far should you stray from what you know?. Ecol Model. 2009, 220: 589-594. 10.1016/j.ecolmodel.2008.11.010.View ArticleGoogle Scholar
- Hijmans RJ, Graham CH: The ability of climate envelope models to predict the effect of climate change on species distributions. Glob Chang Biol. 2006, 12: 1-10. 10.1111/j.1365-2486.2006.01256.x.View ArticleGoogle Scholar
- Global Land Cover by National Mapping Organizations (GLCNMO). , [http://www.iscgm.org/]
- Consortium for spatial information (CGIAR-CSI). , [http://www.cgiar-csi.org]
- Corine Land Cover 2006. , [http://sia.eionet.europa.eu/CLC2006]
- Murphy AH, Winkler RL: A general framework for forecast verification. Mon Weather Rev. 1987, 115: 1330-1338. 10.1175/1520-0493(1987)115<1330:AGFFFV>2.0.CO;2.View ArticleGoogle Scholar
- Pearce J, Ferrier S: Evaluating the predictive performance of habitat models developed using logistic regression. Ecol Model. 2000, 133: 225-245. 10.1016/S0304-3800(00)00322-7.View ArticleGoogle Scholar
- Warren D, Seifert S: Ecological niche modeling in Maxent: the importance of model complexity and the performance of model selection criteria. Ecol Appl. 2011, 21: 335-342. 10.1890/10-1171.1.PubMedView ArticleGoogle Scholar
- Rebelo H, Jones G: Ground validation of presence-only modelling with rare species: a case study on barbastelles Barbastella barbastellus (Chiroptera: Vespertilionidae). J Appl Ecol. 2010, 47: 410-420. 10.1111/j.1365-2664.2009.01765.x.View ArticleGoogle Scholar
- Boubli J, de Lima M: Modeling the geographical distribution and fundamental niches of Cacajao spp. And Chiropotes israelita in Northwestern Amazonia via a maximum entropy algorithm. Int J Primatol. 2009, 30: 217-228. 10.1007/s10764-009-9335-4.View ArticleGoogle Scholar
- Wisz M, Hijmans R, Li J, Peterson A, Graham C, Guisan A: Effects of sample size on the performance of species distribution models. Divers Distrib. 2008, 14: 763-773. 10.1111/j.1472-4642.2008.00482.x.View ArticleGoogle Scholar
- Giovanelli J, De SM, Haddad C, Alexandrino J: Modeling a spatially restricted distribution in the Neotropics: How the size of calibration area affects the performance of five presence-only methods. Ecol Model. 2010, 221: 215-224. 10.1016/j.ecolmodel.2009.10.009.View ArticleGoogle Scholar
- Burnham KP, Anderson DR: Model Selection and Multi-Model Inference. 2002, Springer, BerlinGoogle Scholar
- Engler R, Guisan A, Rechsteiner L: An improved approach for predicting the distribution of rare and endangered species from occurrence and pseudo-absence data. J Appl Ecol. 2004, 41: 263-274. 10.1111/j.0021-8901.2004.00881.x.View ArticleGoogle Scholar
- Elith J, Phillips SJ, Hastie T, Dudk M, Chee YE, Yates CJ: A statistical explanation of MaxEnt for ecologists. Divers Distrib. 2010, 17: 43-57. 10.1111/j.1472-4642.2010.00725.x.View ArticleGoogle Scholar
- Liu C, White M, Newell G: Selecting thresholds for the prediction of species occurrence with presence-only data. J Biogeogr. 2013, 40: 778-789. 10.1111/jbi.12058.View ArticleGoogle Scholar
- Yost AC, Petersen SL, Gregg M, Miller R: Predictive modeling and mapping sage grouse (Centrocercus urophasianus) nesting habitat using Maximum Entropy and a long-term dataset from Southern Oregon. Ecol Inform. 2008, 3: 375-386. 10.1016/j.ecoinf.2008.08.004.View ArticleGoogle Scholar
- Moratelli R, Andreazzi CS D, de Oliveira JA, Cordeiro J: Current and potential distribution of Myotis simus (Chiroptera, Vespertilionidae). Mammalia. 2011, 75: 227-234. 10.1515/mamm.2011.028.View ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.