The Mediterranean basin is considered a hotspot of biological diversity with a long history of modification of natural ecosystems by human activities, and is one of the regions that will face extensive changes in climate. For 181 terrestrial mammals (68% of all Mediterranean mammals), we used an ensemble forecasting approach to model the future (approx. 2100) potential distribution under climate change considering five climate change model outputs for two climate scenarios. Overall, a substantial number of Mediterranean mammals will be severely threatened by future climate change, particularly endemic species. Moreover, we found important changes in potential species richness owing to climate change, with some areas (e.g. montane region in central Italy) gaining species, while most of the region will be losing species (mainly Spain and North Africa). Existing protected areas (PAs) will probably be strongly influenced by climate change, with most PAs in Africa, the Middle East and Spain losing a substantial number of species, and those PAs gaining species (e.g. central Italy and southern France) will experience a substantial shift in species composition.
The Mediterranean basin is one of the richest and most complex regions on the Earth from a geological, biological and cultural point of view . Widely considered as one of the global hotspots of biodiversity [2,3], the region hosts about 25 000 plant species (50% of which are endemic), more than 150 000 insect species (on average 15–20% of endemics; up to 90% of endemics in cave systems) and more than 1100 terrestrial vertebrates (endemism rates range from 17% for breeding birds up to 64% for amphibians; ). A number of biogeographic, geological and historical factors contribute to this exceptional biological diversity [1,4–8], but the tight integration between traditional human activities and Mediterranean ecosystems may lead to coevolution as an additional important driver .
In fact, landscapes in the Mediterranean have been (re)shaped by humans for at least the last 10 000 years [10,11], significantly longer than in any other global biodiversity hotspot. The interaction of natural ecosystems with traditional human activities is one of the reasons for the high environmental diversity that characterizes the region, but it is also the reason why the Mediterranean basin is considered one of the four most significantly altered hotspots on Earth [3,12].
Currently, the high biodiversity that characterizes the Mediterranean basin is highly endangered. Changes in land-use and human population density in the last 50 years have been dramatic [13,14], with both positive and negative consequences for biodiversity [15–17]. At the same time, the Mediterranean region is already home to some 455 million human inhabitants and it is the first vacation destination in the world with 31 per cent of all tourists worldwide visiting the region in 2006 . Moreover, future projections predict a more homogeneous and more densely populated landscape, with tourism expected to double by 2025  with associated overall negative consequences for biodiversity, especially in North Africa and in the Middle East.
The Mediterranean basin is also considered to be one of the regions that will face the largest changes in climate worldwide . The region has already experienced large climate shifts in the past . However, by the end of the twenty-first century, Gao & Giorgi  have projected a substantial northward increase of dry and arid lands, with the frequency of hot extremes increasing by 200–500% through the region , with most of the changes occurring in summer and spring , and with an observed trend of recent warming that is much higher than was simulated by the state-of-the-art global circulation models over the recent past .
All of these factors support an evaluation of the available conservation options that include projections of ongoing global change. A number of studies have clearly demonstrated the negative effects that climate change will have and is already having on biodiversity [26,27], and some species have already gone extinct owing to recent climate change . However, only limited attention has been given to the effects of future climate changes on mammals [29–33], and to our knowledge no one has considered the entire Mediterranean basin specifically.
Recently, predictive species distribution models (SDMs) have become increasingly important to address a number of issues in ecology, conservation biology and climate change research . We applied the state-of-the-art bioclimatic SDMs  to evaluate the potential effects of climate change on Mediterranean mammals. In particular, we calculated changes in climatically suitable areas for all species, classifying the species into different threat categories depending on the amount of (positive or negative) changes. Furthermore, we calculated species turnover for each pixel in the Mediterranean basin, considering also the existing protected areas (PAs). All results were analysed considering averages across climate change projections originating from five general circulation model (GCM) outputs using two emission scenarios (A2 and B1) and considering two contrasting assumptions on dispersal ability (species unable to disperse versus species dispersing with no constraint [36,37]).
2. Material and methods
(a) Study area and protected areas
We defined our study area (approx. 3 500 000 km2; figure 1) following the boundaries of the Mediterranean biogeographic area . The region, stretching over approximately 5300 km west–east and 2200 km north–south, touches three continents (Europe, Asia and Africa) and 24 different countries.
From the World Database on Protected Areas (http://protectedplanet.net/), we extracted 1390 national and international PAs occurring inside our study area. Moreover, for the European part of the Mediterranean basin, we also obtained a database on the Natura 2000 network (1260 areas inside our study area), a network of areas whose aim is to conserve an extensive range of habitat types and wildlife species throughout Europe under the framework of the European Union Birds and Habitat Directives (http://www.eea.europa.eu/data-and-maps/data/natura). Overall, PAs plus the Natura 2000 cover 8 per cent of the study area, with the large majority of the areas being in Europe (only 26 areas in Africa, corresponding to 0.25% of the total extensions of PAs). PAs and Natura 2000 were overlaid to the final maps (see below) to identify the most critical areas for mammal species conservation.
(b) Species data
To avoid including in the analyses species only marginally present in the Mediterranean basin, we initially considered only the 212 species (electronic supplementary material, table S1) that according to Schipper et al.  have at least 5 per cent of their distribution range inside the region (79.8% of all species occurring in the Mediterranean). On average, more than 47 per cent of the distribution range for these 212 species fall inside the Mediterranean basin, with 18 species being sub-endemic of the study area (i.e. more than 90% of the range is included in the study area) and 35 being endemic.
Rodentia represents the order with the highest number of species (92), and Muridae is the most represented family (43 species). Considering the most recent IUCN Red List available , three species (one of which is endemic) are classified as Critically Endangered (CE), seven as Endangered (EN) (three endemics), 19 as Vulnerable (VU) (eight endemics or sub-endemics) and one endemic species is considered Extinct in the Wild (electronic supplementary material, table S1).
For each species, we obtained the current distribution range from Schipper et al.  that was used later in the analysis to define different dispersal scenarios. Moreover, we collected all the freely and readily available points of presence from a number of published data sources (main references [40–43]; complete list available from L.M.) and online databases (http://data.gbif.org), and, when available, we also obtained the positional error associated with each point (from 300 m to 50 km). To obtain a more general representation of each species' climatic niche, presence points were collected considering an area larger than the final study area, corresponding roughly to the western Palearctic. We excluded 30 species (6 endemics or sub-endemics) from further analysis for which we found very few points of presence to cover most of the distribution range included in the study area. We also excluded Oryx dammah, which is extinct in the wild, obtaining a final set of 181 species (85.4% of the initial set of 212 species; table 1).
(c) Current climate and climate projections
We considered a set of five bioclimatic variables that are important to shape the climate of the study area, and which are at the same time not highly correlated (i.e. |Pearson r| ≤ 0.8). The Mediterranean climate is transitional between cold temperate and dry tropical, with a unique combination of hot and dry summers, and cool (or cold) and humid winters. Given the limited availability of surface water during the months of maximal solar insolation, the long and dry summer periods (two months in the west and five or six months in the east with basically no precipitation) represent the most unfavourable period of the year for both plants and animals, contrary to the regions further north where the unfavourable period is winter owing to low temperatures. Accordingly, the weather during the short spring and autumn seasons is extremely important, representing a ‘climatic buffer’ against the extreme summer climate.
In a recent review on climate change projections for the Mediterranean region, Giorgi & Lionello  considered the most recent and comprehensive ensembles of global and regional climate change simulations and obtained a robust and consistent picture for the future Mediterranean climate. The main features consist of an important decrease in precipitation, especially in summer, and of a pronounced warming, also to occur mostly during the summer season. Moreover, inter-annual variability is projected to increase with a greater occurrence of extremely high temperature events.
Given these characteristics, we considered the following variables: the average moisture index of summer months (MS); the mean temperature of the coldest quarter (MTC); the mean temperature of the driest quarter (MTD); annual precipitation (ANNP); and precipitation of the driest quarter (PDQ). We obtained all climate variables at 1 km resolution from WORLDCLIM , with the exception of the moisture index that was calculated as the difference between precipitation and potential evapotranspiration each summed over the three months June, July and August. Potential evapotranspiration was calculated following the empirical equation of Turc , which considers monthly average temperature and potential global radiation. The latter was derived in a geographic information system using the method by Kumar et al. . More information is given in Zimmermann et al. .
Future projections were derived using climate model outputs made available through the Intergovernmental Panel on Climate Change (IPCC) Data Distribution Centre (http://ipcc-ddc.cru.uea.ac.uk). In particular, we used five climate change model outputs that form part of the fourth assessment report  (CGCM31 run by the CCMA, CSIRO's MK35, CM21 run by the GFDL, PCM1 operated by NCAR and UKMO's HADCM3) for two of the IPCC's climate scenarios, the B1 scenario, describing a world with reduced use of natural resources and the use of clean and resource-efficient technologies, and the A2 scenario, where the greenhouse gas emission rate continues to increase . Climate variables were averaged over the period 2071–2100.
The original climate change models came with varying resolutions of roughly 2 × 2°, roughly corresponding in the study area to 180 × 200 km. The Mediterranean basin is characterized by very high topographic complexity, which renders such a resolution too coarse to adequately depict the orographically induced fine-scale structure of the local climate. Thus, we downscaled the available climatic layers to 1 km2 by calculating anomalies of the future monthly temperature and precipitation values against the 1950–2000 means of each GCM model output, where the latter represents the WORLDCLIM base period. Anomalies represent absolute temperature difference (Δ°C) and relative precipitation differences (% change) per coarse resolution pixel measured directly at the model output. Second, we down-scaled these anomalies in two interpolation steps to 1 km2 of spatial resolution, using bilinear resampling. Third, we added temperature anomalies to WORLDCLIM and multiplied precipitation with the respective anomalies to obtain monthly maps of future temperature and precipitation for each model and scenario at a 1 km2 spatial resolution. Finally, we re-calculated all derived climate maps mentioned above. The applied procedure has the advantage that a possible model offset under current climate is not added to the projected climate trends.
(d) Bioclimatic distribution models
We predicted mammal distributions using an ensemble forecasting approach, as implemented in BIOMOD, a bioclimatic niche modelling package , within the R environment. We used the following eight models: (i) generalized linear models, (ii) generalized additive models, (iii) classification tree analysis, (iv) artificial neural networks, (v) generalized boosted models, (vi) random forests, (vii) flexible discriminant analysis, and (viii) multivariate adaptive regression spline. See Thuiller et al.  for more details.
The location error associated with the available points of presence was often larger than the resolution adopted in our analysis, especially for points collected through the Global Biodiversity Information Facility. Thus, around each point of presence, we considered a circular buffer with a radius corresponding to the location error. We then generated a random point inside each buffer to which we associated the value of the climatic variables. We repeated the same procedure 10 times for each species, obtaining 10 sets of points representing species' presence and accounting for positional errors. We also generated 10 sets of 20 000 random background points (drawn from the western Palearctic, i.e. the same area over which we collected presences) that we used to characterize the climate of the study area and to represent pseudo-absences. We combined the 10 sets of points representing species presence with the 10 sets of background points, obtaining 100 sets of points for each species to be used for the calibration of each of the eight niche models. For each of the 100 sets of points, we calibrated the eight models using a 70 per cent random sample, and we evaluated model performance using the remaining 30 per cent of the data. In particular, we measured the area under the receiver-operating characteristic curve (AUC) . Only models with AUC values greater than 0.7  were projected onto the available future climate scenarios.
Probability of occurrence for each species and for each model was transformed into presence–absence values using a threshold maximizing the sum of sensitivity (the percentage of presence correctly predicted) and specificity (the percentage of absence correctly predicted) . Then, a final consensus model for each species was obtained defining as presences all those cells where more than 50 per cent of the models (i.e. five out of eight) predict species' presence. Using the same thresholds, we obtained binary models also for climate change projections, for a total of 160 models for each species (five climate change models, by two IPCC scenarios, by eight modelling techniques, by two dispersal scenarios). We produced four consensus forecasts for each species, one for each possible combination of two IPCC emission scenarios (A2 and B1) and of two dispersal scenarios, i.e. no-dispersal (assuming that species will persist only in suitable areas inside their current range) versus unlimited dispersal (assuming that a species could disperse to all new suitable areas inside the Mediterranean basin). These assumptions on dispersal are commonly used in the literature (e.g. ), and represent two extremes, while the realized patterns at the end of the twenty-first century will necessarily be somewhere in between depending on the species' capacity to track climate change. For each ensemble forecast we defined as presences all those cells where more than 50 per cent of the models (i.e. 21 out of 40 available for each scenario) predict species' presence.
(e) Changes in climatically suitable areas and threat categories
We calculated the percentage of cells gained or lost by each species under the B1 and A2 scenarios with respect to the initial climatically suitable range, corresponding to the current distribution range under the no-dispersal hypothesis, and to the full study area under the unlimited dispersal hypothesis. Moreover, following Thuiller et al. , we defined the following threat categories based on the IUCN Red List Criteria (, but see  for a critique of this approach): threat-level 1 (Threat-1; species with a projected range loss of 100%), threat-level 2 (Threat-2; species with a projected range loss >80%), threat-level 3 (Threat-3; species with a projected range loss between 50 and 80%) and threat-level 4 (Threat-4; species with a projected range loss between 30% and 50%). All other species were classified either as threat-level 5 (Threat-5; species losing less than 30%) or as gain (species actually gaining climatic suitability).
(f) Species richness, species turnover and protected areas
Considering all species, we calculated both the current and the future potential species richness and species turnover (measured as percentage of species gained or lost at each cell over the corresponding potential species richness). All calculations were performed following recommendations in Thuiller et al. .
To assess the sensitivity of PAs to climate change, we tested if PAs will experience a higher or lower species turnover if compared with a random set of areas sharing the same characteristics as PAs. In particular, we generated 1000 sets of random areas occupying for each country the same area as the existing PAs. For each set of random areas, we calculated the average turnover, obtaining the distribution that would be expected for our study area under the null hypothesis stated above. We calculated the same metric for existing PAs and we rejected the null hypothesis when the metric measured inside PAs was in the top 5 per cent of the null distribution. Always considering PAs, we also compared the current species richness with the richness projected under the different emission and dispersal scenarios. All calculations were performed for both the B1 and the A2 scenarios, and for both dispersal scenarios.
(a) Model evaluation
All models showed a good predictive power with a mean AUC of 0.96 (±0.042). The minimum AUC value was 0.79 for the Tarabul's Gerbil (Gerbillus tarabuli), while the maximum AUC value was 0.99 for the James's Gerbil (Gerbillus jamesi). The average AUC value for endemic and sub-endemic species was 0.98, with a maximum of 0.99 and a minimum of 0.94. The same values were obtained for species classified as Critically Endangered, Endangered or Vulnerable, with a mean AUC of 0.98, a minimum of 0.94 and a maximum of 0.99, respectively.
(b) Changes in climatically suitable areas and threat categories: unlimited dispersal
Assuming unlimited dispersal, both positive and negative changes in climatically suitable area are possible (table 1), with Rhinopomatidae being the family with highest gains in climatically suitable areas, and Talpidae and Cervidae showing the highest losses. The comparably mild B1 climate scenario projects that changes in the extent of climatically suitable habitats would vary on average between −91 and +880 per cent (mean = −11%; s.d. = 76%) among all mammals, while under the more severe A2 scenario the projected changes are more extreme, ranging from −100 to +990 per cent (mean = −25%; s.d. = 96%). With few exceptions (e.g. Lepus capensis and Mus cypriacus), all species gaining climatically suitable areas are strongly associated with desert and/or steppe environments.
Considering the overall average across all climate change models, one species only (Sorex arunchi) was predicted to lose completely its climatically suitable range in the future (Threat-1) under the A2 climate scenario. Considering the single climate change models, another 15 species were predicted to become extinct in one or more climate change models under the A2 scenario, while S. arunchi was classified as extinct under the MK35 climate model for the B1 scenario.
Under the A2 scenario (figure 2 and electronic supplementary material, table S2), 8.3 per cent of species were classified in Threat-2, 37.2 per cent in Threat-3 and 20.6 per cent in Threat-4, percentages much higher than the comparable classes in the current red listing (1.0% CE; 3.3% EN; 9.4% VU). Considering endemic and sub-endemic species, 73.9 per cent were classified as Threat-2, Threat-3 or Threat-4 as opposed to 21.7 per cent currently in the comparable IUCN categories (i.e. CE, EN or VU).
Projections for the B1 scenario (figure 2 and electronic supplementary material, table S2) produced only 0.5 per cent of species predicted in Threat-2, 8.8 per cent in Threat-3 and 23.8 per cent in Threat-4, percentages that in any case are bigger than the current red list classifications. Considering endemic and sub-endemic species, 30.4 per cent were classified as Threat-2, Threat-3 or Threat-4. For both A1 and B2, some species (respectively, 19.9 and 22.7%) are predicted to gain climatically suitable areas.
(c) Changes in climatically suitable areas and threat categories: no dispersal
Assuming no dispersal, even if some species would still gain climatically suitable areas, the results predict stronger negative trends (table 1), with Hyaenidae showing the highest gains in climate suitability and Ctenodactylidae and Talpidae showing the highest losses. Under the B1 climate scenario, changes in the extent of climatically suitable areas would vary on average between −100 and +110 per cent (mean = −22%; s.d. = 28%), while under the A2 changes would range from −100 to +90 per cent (mean = −45%; s.d. = 33%). Also in this case, species predicted to gain climatically suitable areas are mostly associated with deserts and/or steppes, while S. arunchi is the only species projected as extinct under A2 (i.e. Threat-1; figure 2 and electronic supplementary material, table S2). Moreover, 17 species were predicted to lose all their climatically suitable range (corresponding to Threat-1) under one or more climate change models for the A2 scenario, and four showed the same results under the B1 scenario.
Considering the A2 scenario (figure 2 and electronic supplementary material, table S2), the projected losses are much stronger if compared with the unlimited dispersal hypothesis, with 14.9 per cent of Threat-2 species, 37.2 per cent of Threat-3 and 22.2 per cent of Threat-4. Comparable results were obtained for the B1 scenario (figure 2 and electronic supplementary material, table S2), with 2.8 per cent of species in Threat-2, 11.1 per cent in Threat-3 and 18.2 per cent in Threat-4. The same trend towards stronger negative changes was measured for endemic and sub-endemic species, with 34.8 per cent of species in Threat-2, Threat-3 or Threat-4 under the B1 scenario, and 91.3 per cent under the A2 scenario.
(d) Species richness, species turnover and protected areas
Under the current climate, the highest values of potential species richness were located in the northern part of the Mediterranean (electronic supplementary material, figure S1a), and particularly in the Italian peninsula, in northern Spain, in the internal part of Greece and in northern Turkey. Other areas of high species richness were predicted in the main islands (Sardinia, Sicily and Corsica) and along the coasts of Morocco, Algeria and Tunisia in Africa. The pattern of species richness is similar under both future climate scenarios (electronic supplementary material, figures S1b,c), with a strong decrease in richness that is predicted for Africa and Spain, and with Italy that should be hosting most of the areas with the highest species richness, representing by far the main centre for mammal diversity in the future.
Considering the unlimited dispersal scenario, our models predict that the highest species turnover, corresponding to a net species lost, will be localized especially in Spain and in northern Africa (B1 climate scenario; figure 3a), while areas with high species turnover in the Middle East, in Turkey, in Greece and in central Italy represent mainly areas of positive changes in species richness (figure 3a). Under the A2 climate scenario (figure 3b), the turnover of species richness is predicted to be more extensive, but it should correspond to a positive change in species richness only for a few small areas in central Italy, France, Turkey and in the Middle East. Species turnover under the no-dispersal scenario (figure 3c,d) would basically show the same pattern described for the unlimited dispersal scenario, but with lower intensities as was expected given the assumption on lack of dispersal.
Considering the average species turnover measured inside PAs, we found no statistically significant difference and thus we could not reject the null hypothesis that PAs are significantly different from a set of random areas. The same general pattern can be inferred considering species richness (figure 4), even with the obvious differences existing throughout the study area. PAs located in drier parts of the Mediterranean basin were most sensitive to climate change, including mainly PAs in northern Africa and the Middle East, where a very high number of species may be lost (on average from 37% to 25% decrease in species richness depending on the scenario), but also in southern Europe, mainly Spain (with an average decrease ranging from 43% to 29%). PAs located in central and northern Italy (along the Apennines), as well as in France (excluding the PAs along the coasts) were projected to gain a substantial number of species (average increase in species richness from 10% to 20%).
In this study, we presented a comprehensive analysis of the potential effects of climate change on Mediterranean mammals. Our results should be considered with particular attention to three main limitations of the approach we adopted [34,54–56]. First, we focused on climatic variables only, even though climate is only one of the many determinants of species distribution [34,57–59]. However, considering the large spatial scale of our analysis, species distribution, even for mammals, is likely to be dominantly determined by climate [29,60,61]. Moreover, the direct effect of climate, especially rainfall, on mammal distribution and biology has been widely documented for a number of species and study areas [31,62,63], and Kearney et al.  clearly demonstrated, although for a single species limited to Australia, that physiologically based mechanistic models and correlative distribution models may provide congruent forecasts under climate change, further supporting the approach we adopted. Finally, the Mediterranean basin is going to be an important hotspot for climate change , with significant climate shifts being projected by 2100 [22–24]. In addition, human pressure is already stressing many mammal species, possibly enhancing the already negative influence that climate change may have.
Second, we did not consider any evolutionary and/or phenological responses to climate change . Réale et al.  demonstrated that the changes in breeding dates as a function of climate change for a population of Red Squirrel in the Arctic was for 62 per cent a result of phenotypic plasticity and for 13 per cent a result of genetic change. More examples are reviewed in Parmesan , but despite the existence of data supporting the idea of local adaptations to climate change evolving at specific sites, the fossil record shows only little evidence for an evolutionary response. During Pleistocene glaciations, for example, flora and fauna experienced climate changes 5–10 times stronger than the magnitude of the twentieth century global warming, and the main response was the shift of geographical distributions to track the changing climate .
A final limitation is related to the assumptions about the dispersal capacity of species through the Mediterranean basin, an assumption that will strongly influence the possibility to track climate change. The possible consequences for ecological communities are likely to be unpredictable  but the current knowledge is so poor that it is difficult to get a reliable estimate of potential dispersal distance for the species we considered. Assuming two extreme scenarios (no dispersal versus full dispersal) is debatable but at least covers the full range of possible outcomes. The available literature, however, leans more towards a limited-dispersal scenario for the Mediterranean basin. In fact, all projection of future land use and human pressure on the Mediterranean basin (e.g. [18,19]) predict a landscape with less and less natural areas and with an increased human pressure, making future habitat connectivity an important issue and questioning the effective possibility of species to track the changing climate.
We considered in our analysis only species whose distribution model gave an AUC score greater than 0.7, and we actually obtained AUC scores greater than 0.79 (average of 0.96) for our models. The lower AUC values were obtained for species only marginally present in the Mediterranean, such as Gerbillus tarabuli (Tarabul's Gerbil) with only 10 per cent of the distribution range inside the study area, or for species with particular habitat requirements, such as Lutra lutra (European Otter).
We produced our models at a fine spatial resolution (1 km2) to avoid the problem of hidden local refugia, one of the most important issues in climate change projections based on correlative distribution models . Moreover, we developed our models considering an ensemble modelling framework, and accounted for uncertainty in both statistical models and climate changes models . However, it should be recognized that the location errors associated with the available points of presence were often larger than the 1 km2 resolution that we adopted. Although we explicitly accounted for this factor, we cannot avoid possible negative effects in the analyses originating from these positional errors, which possibly result in smoothing the response to climate that we had measured for all species using SDMs.
Of the 181 species that we considered, only the Lowland Shrew (S. arunchi; but see  for the taxonomic validity of the species) was committed to extinction considering the final ensemble forecasts, whatever scenario was considered for dispersal (no-dispersal versus full dispersal). However, the variability between the different climate change models was large (see electronic supplementary material, figure S1d,e for a representation of the uncertainty in estimates of species richness), with the number of species predicted to become extinct (i.e. Threat-1) under one or more climate change model, but not under the final ensemble forecast. Overall, the uncertainty linked to the different climate change models that we considered is important for all scenarios, with the majority of species being projected into more than one threat level (90.6% under A2 no-dispersal; 85.6% under A2 unlimited dispersal; 74.4% under B1 no-dispersal; 71.1% under B1 unlimited dispersal) and many species spanning at least three different threat levels (43.3% under A2 no-dispersal; 22.8% under A2 unlimited dispersal; 5.6% under B1 no-dispersal; 2.8% under B1 unlimited dispersal). The species Mesocricetus auratus provides an extreme example of this variability: under the B1, unlimited dispersal scenario, the species was predicted to lose more than 80 per cent of its climatically suitable area according to the MK35 climate change model, more than 50 per cent under the PCM1 model, while under the CM21 model its climatically suitable range was predicted to increase. These results confirm what has already been found in previous analyses (e.g. [67,69]) and further stress the importance of adopting a full ensemble forecasting approach in studying the effects of climate change on species distribution .
Overall, a substantial number of Mediterranean mammals will be severely threatened by future climate change (electronic supplementary material, table S2), particularly endemic species. Moreover, additional drivers could enhance the risk of species extinction or at least could increase further the uncertainty linked to our results (e.g. parasites, fire frequency, diseases; ). Particularly important is the lack of consideration for the effects of land use and other types of human pressure. Currently, only 5 per cent of the original natural vegetation is remaining and, while in the northern part of the basin land use changes are going towards a more natural landscape (e.g. ), in the Middle East, as well as across North Africa, land use change and human pressure are intensifying, with potentially important and negative effects on biodiversity. This is obviously going to exacerbate the effects of climate change, at least for some species, particularly considering the almost complete absence of PAs in the African part of the Mediterranean basin.
Our models also predict that a number of species will gain climatically suitable areas as a result of climate change (electronic supplementary material, table S2). Although this is possible for many species, a number of factors may actually influence these results. For example, a species may not be able to colonize all climatically suitable areas because they may be limited by dispersal capabilities, by behavioural patterns or biotic interactions, by habitat fragmentation and natural barriers, or by the absence of environmentally suitable habitats. This may be the case of species inhabiting the Mediterranean peninsulas, who will have difficulty moving northwards owing to east–west-oriented mountain ranges such as the Pyrenees and the Alps. The same is true for species inhabiting Mediterranean islands (among which there are many endemics), which will face even greater barriers, making it basically impossible for them to benefit from gains in climatic suitability.
We found important changes in potential species richness owing to climate change (but see Bartolino et al.  for a discussion on the threshold approach that we used to identify hotspots). As a result of shifts in species ranges, mostly towards higher elevations but also northward, the areas that will gain species are predicted to concentrate around mountainous regions (e.g. central Italy), as already predicted by other studies (e.g. ; but see Thuiller et al.  for the opposite result). In the European part of the study area, our models predict species losses from peninsulas and islands, with the main hotspots of species loss being concentrated in the western Mediterranean basin (mostly Spain and North Africa). Species turnover will be higher under the unlimited dispersal assumption (and especially for the A2 climate scenario), a result that was clearly expected given the complete lack of movement constraints. The areas with the highest turnover rates (figure 3) basically correspond to the hotspots of species loss (European peninsulas and islands, and northern African mountain ranges). Obviously, these estimates should be considered carefully as it is possible that species from regions around the Mediterranean might be able to immigrate to these newly available areas.
The existing PAs would be strongly influenced by climate change according to our results. Basically all PAs in Africa, in the Middle East and in Spain would lose a substantial amount of their current mammalian species diversity, with significant losses that are not compensated by species influxes. Overall, also those PAs that would gain species or remain stable in central Italy and southern France should experience a substantial shift in species composition of a magnitude unprecedented in recent geological times , with potentially important consequences because of changes in the pattern of interspecific interactions . These conclusions are based on the assumption that species will be able to track shifts in habitats/climate with no lagged response. These assumptions are debatable in general , although rapid range adjustments are not entirely improbable for mammals . Moreover, climate warming is likely to result in changes in species phenology, which can further disrupt species associations . The outcome of such new species assemblages may be particularly difficult to predict, for example because of nonlinearities that may emerge.
Managing natural systems in the face of such changes is a challenging task, especially as decisions need to be taken using limited and highly uncertain projections of future climate impacts . However, a few indications are possible given the results of our analyses for Mediterranean mammals. In particular, it is clear that the adopted conservation strategy should follow different paths in the European and in the African/Asian parts of the Mediterranean. In the Mediterranean Europe, the number and extent of PAs make it almost impossible to propose expanding the existing networks. However, new strategies could be considered. First of all, it seems important to avoid non-climate-related threats, like invasive species, further habitat loss and increased fragmentation, overharvest, etc. . Moreover, we suggest using tools such as the Natura 2000 network to increase habitat connectivity between PAs . Finally, we also propose considering adaptive management strategies both inside and outside existing reserves as a critical step forward .
The PA system in northern Africa is completely different. Here, the PA network is limited to few small and mostly spatial disconnected areas. Establishing new PAs is the first and most important need. Particular attention should be given to the inclusion of high habitat heterogeneity and, possibly, on large reserves across biotic transition zones that would allow for the continued protection under ongoing climate change . Furthermore, a focus on the socio-economic components (i.e. increasing socio-economic conditions over a simple subsistence level, introducing eco-tourism, enforcing conservation measures, etc.) would represent a clear step forward to ensure long-term conservation success.
Our results suggest that the effects of climate change on mammal species distributions and communities they form may not be a drastic loss of species from their current distribution, but a change in community structure. The results of this study should be considered as a first approximation but the main pattern that we obtained should remain stable even considering possible improvements (e.g. better presence/absence data ; mechanistic models ). Species richness is predicted to shift and the extinction risk for many species is going to increase. In this context, the current network of PAs is not expected to effectively conserve existing species, as most PAs will potentially lose species according to our analyses. Expanding the network (e.g. more PAs are urgently needed in northern Africa) and complementing it with other strategies should be urgently considered. In particular, we suggest a stronger emphasis on off-reserve conservation, explicitly considering connectivity among PAs, as well as management and mitigation actions, among which going back to the traditional agricultural and husbandry practices may represent an important first step.
L.M., A.G., N.E.Z., A.P. and J.P. were funded under the FP6 EU-funded ECOCHANGE project (Challenges in assessing and forecasting biodiversity and ecosystem changes in Europe; number 066866 GOCE). David M. Theobald and three anonymous referees provided helpful comments on a previous version of the manuscript.
- This journal is © 2011 The Royal Society