Climate warming drives divergence of montane butterfly communities in Southern Italy

Because of climate change, many insect species are shifting their altitudinal and latitudinal ranges, including Mediterranean butterflies, particularly in mountainous regions. In this study, we evaluated changes in butterfly communities over time, sampled in 1975, 2004, and 2012, in relation to their altitude and two indices representing the climate envelopes of species within a given community: CTI (Community Temperature Index) and CPI (Community Precipitation Index). The study took place in a protected area where we found strong changes in community compositions over the 37-year study period. There was no vertical stratification of communities in 1975, but became significantly so in 2004 and 2012. Likewise, CTI and CPI were correlated with altitude only in 2004 and 2012. Over time, CTI increased at lower altitudes, indicating an increase in species associated with higher temperatures, and was stable or decreased at higher altitudes. CPI showed opposing trends, decreasing at lower altitude of communities and increasing in higher altitude communities. This resulted in asymmetric changes along the altitudinal gradient. The highest elevations (>1900 m) shifted towards butterfly species that are more associated with colder, wetter habitats, and lower elevations shifted towards species more associated with hotter, drier habitats. In conclusion, changes in butterfly communities were consistent with expectations from observed changes of temperatures and precipitations at low altitudes and mid-altitudes, but not at the highest altitudes. This counter-intuitive result may be due to land-use changes following creation of a national park in 1993 that encompassed the sample sites, but we lack the data to test this hypothesis.


Introduction
This climatic trend of the past decades and the increase of frequent extreme weather events (EWEs) have significantly modified the distribution of biota around the World (Thomas 2010). Climatic projections for the next future depict an increase of warm and dry conditions, at a much higher rate than those observed in the recent past (Siegmund et al. 2019). The causal effects of EWEs can be detected more easily than those of global climate change as the former can be explored by planning ad hoc studies also after their occurrence (Altwegg et al. 2017), whilst the latter can only be inferred, thanks to the presence of longterm monitoring programs or, at least, when past standardized data are available. The Mediterranean Basin is regarded as one of the hotspots of global warming on Earth (Giorgi 2006) because it is subject to an increase of temperature and a decrease of precipitation, and to an increasing frequency of EWEs (Diffenbaugh et al. 2007;Gao and Giorgi 2008). Thanks to the presence of several glacial refugia, the Mediterranean Basin also represents an important biodiversity hotspot (Myers et al. 2000) hosting several endemic species (Médail and Quézel 1999;Médail and Diadema 2009;Blondel et al. 2010). Despite the importance of this geographic area for the conservation of biodiversity, longterm monitoring schemes involving the effects of climatic changes on animal species are only available for few Mediterranean countries, even in the case of appealing taxa such as butterflies (van Swaay et al. 2008 Lepidoptera monitoring has proven to be a useful tool for evaluating the effects of climate change around the World (Parmesan et al. 1999;Roy and Sparks 2000;Visser and Holleman 2001;Wilson et al. 2005;Schweiger et al. 2008;Bell et al. 2019). A study carried out in Spain demonstrated that over the past decades, butterfly communities shifted upward by approximately 293 m in mountainous regions, resulting in a net loss of diversity (Wilson et al. 2007). At lower altitudes, a decline of high-elevation species and a general increase in semi-natural habitats of thermophilous ones have been observed, suggesting a causal effect of climate warming on the composition of butterfly communities (Zografou et al. 2014). Short-term changes are less pronounced and cannot be univocally attributed to the recent increases in temperature, as individual species population dynamics can conceal such low level of changes (Zografou et al. 2014;Cerrato et al. 2019). Studies concerning individual species confirm the results obtained for the entire species assemblages, showing an elevational upward shift of about 200 m of their range (Wilson et al. 2005;Scalercio et al. 2014). Nevertheless, the upward shift of species distribution is not always coupled with a reduction of population size. Interestingly, Erebia cassioides showed an increased population size coupled with a range contraction to the top of the Pollino Massif in southern Italy, suggesting asymmetric effects of climate changes along the altitudinal gradient (Scalercio et al. 2014).
An increasing number of papers are devoted to describing the effects of climate change on the latitudinal and altitudinal range of butterflies. Studies across southern Europe (Wilson et al. 2005(Wilson et al. , 2007Zografou et al. 2014;Cerrato et al. 2019), and other geographic areas (Pounds et al. 1999;Shoo et al. 2006;Despland et al. 2012;Molina-Martínez et al. 2013;Dewan et al. 2019), tend to lump high-elevation sites together into one "high" altitudinal band which fails to consider potential differences in species' responses to climate change within montane habitats. Specifically, mountaintop alpine environments often differ considerably from those just at or below treeline, although they are all typically considered "high elevation" habitats.
The Calabria region, the southernmost tip of the Italian Peninsula, is a perfect area to study the consequences of global warming as in southern Italy precipitations decreased by 10% relative to the mean of the standard period 1961-1990 and the mean minimum temperatures increased by 1.3°C (Brunetti et al. 2006). Climate projections for the 21st century also show a significant warming and a general reduction of precipitation (Bucchignani et al. 2016) for this geographic area. Moreover, it represents one of the most important glacial refugia for continental biota (Konnert and Bergmann 1995;Médail and Diadema 2009;Grill et al. 2009;Andersen et al. 2019), with a strong altitudinal gradient, hosting several populations of Lepidoptera at their southern range limit (Scalercio 2009;Infusino et al. 2016;Greco et al. 2018). According to available climatic scenarios, in this area, we can expect a very high rate of local extinctions along with the expected loss of habitat suitable for cold-adapted species (Scalercio 2009;Scalercio et al. 2014). However, the key point that allows us to carry out this study is the availability of standardized quantitative butterfly data collected on the Pollino Massif, the highest southern mountain with relict alpine prairies, in the year 1975 by Balletto et al. (1977). The aim of the current study is to describe ongoing altitudinal changes in butterfly communities in a protected area of a Mediterranean massif, at a finer altitudinal resolution than most other studies. As the establishment of protected areas could be expected to improve butterfly habitat conditions, we hypothesize the mitigation of global warming effects on butterfly communities in the Pollino National Park. After the establishment of the national park in 1993, there was a marked decrease in grazing mainly at high altitudes where the core area of the park is located (Scalercio S., personal observation), but there are no official data for this. Unfortunately, we also lack data concerning changes in vegetation cover, but it is known that decreased cow pressure on vegetation can led to a taller grass cover and to an increased forested surface (Carmel and Kadmon 1999;Sternberg et al. 2000).
To achieve the aims of this study, in 2004 and 2012, we resampled most of the sites monitored in 1975, comparing abundance and presence/absence data along the altitudinal gradient. Then, in order to assess the association between main climatic parameters and butterfly communities, we computed the Community Temperature Index (CTI) and the Community Precipitation Index (CPI), functioning as biothermometers that make it easier to summarize and communicate changes in community compositions.

Study area
The study area was included entirely within the Pollino National Park, one of the largest National Parks in South Italy with a surface of about 190,000 ha ( Fig. 1). It is a mountainous territory with a large area above 1800 m of altitude and five large peaks over 2000 m (Serra Dolcedorme, Monte Pollino, Serra del Prete, Serra delle Ciavole, and Serra Crispo). Alpine prairies are present above the tree line. High-altitude natural grasslands are characterized by Festuca bosniaca and Carex kitaibeliana. Pasturelands inside the beech belt are composed of large clearings covered by Festuca violacea-Meum athamanticum grasslands. The treeline (1800-1900 m) is largely influenced by grazing or other anthropic factors that have led to the development of a wide variety of xeric or mesic herbaceous biotopes in its surroundings. However, although cow grazing was strongly reduced by the establishment of the National Park, the tree line was quite stable as demonstrated by old pictures showing that sub-vertical rocky slopes slow down the expected upward shift of European beech forest (see Tamanini 1961). Karstic depressions around 2000-2100 m host scattered patches of acidophilic Nardus grass mats (Nardo-Luzuletum pindicae, known in the Corine classification as 36.381: Subalpine southern Italian mat-grass swards) (Brandmayr et al. 2002;Pizzolotto et al. 2010).
The climate of the Pollino National Park varies according to the altitude from summer-arid true Mediterranean to mountain Mediterranean. The heaviest rains are concentrated from November to February. Climate data for the area are fragmentary and available for sites at lower altitude than that of the study area, but only few kilometers apart. Temperatures are available for the Castrovillari meteorological station (400 m a.s.l.), while rainfall data are available for Castrovillari and Campo Tenese (1400 m a.s.l.) meteorological stations. Temperature and precipitation trends of the past decades were analyzed by Scalercio et al. (2014), and they seem to be consistent with those found by Brunetti et al. (2006) in Southern Italy as a whole. At low altitudes, temperatures show an evident increase of maxima with a value of −1.5°C below the mean of the 1925-2012 period in 1985 and a value of +1.0°C in 2012. The same occurred for the minima with a value of −1.5°C in 1986 and of +2.1°C in 2012. The rainfall dropped rapidly starting from the 1970s, with about 1700mm (10-year running mean) of precipitation registered in the year 1970 at the Campo Tenese meteorological station, dropping to about 950mm in 2012.

Butterfly monitoring
Communities were sampled for the first time by Balletto et al. (1977) during the year 1975 within 23 sites. We repeated the sampling in most of these sites in 2004 and 2012, after 29 and 37 years respectively. As a result of excluding sites below 1400m and three sites whose exact locations we could not find, we sampled a total of 16 sites ranging from 1470 to 2220 m in altitude (Table 1). Butterfly monitoring was carried out from early-July to late-August as in Balletto et al. (1977), corresponding to the period of maximum abundance and richness of butterflies in Mediterranean mountainous areas. During this period, four sampling sessions were carried out, separated by no less than 14, and no more than 20 days during optimal weather conditions for butterfly activity (sunny days, no or low wind speed). In previous studies on the effects of climate change on Mediterranean butterfly communities along altitudinal gradients, not enough details were provided on changes occurring in the highest altitude communities as sites were grouped in wide altitudinal band (Wilson et al. 2005(Wilson et al. , 2007Cerrato et al. 2019). In this study, we focused our attention only on the highest sites sampled by Balletto and collaborators, grouping sites into altitudinal band: (i) low altitudinal band (1470-1575m a.s.l.) includes 7 sites mainly composed of beech forest margins and pastures, (ii) medium altitudinal band (1610-1825m a.s.l.) includes 5 sites mainly composed of beech forest margins and pastures, and (iii) high altitudinal band (1900-2220m a.s.l.) includes 4 sites all dominated by Seslerio-Brometum grasslands located above the tree-line.
We utilized a different sampling method than that used by Balletto et al. (1977). In fact, whilst Balletto et al. (1977) Fig. 1 Study area is indicated by a white arrow. Numeric codes represent monitored sites and indicate their location quantified species abundances as the number of observed individuals per hectare producing density values, we applied a time-constrained method by sampling for 20 min along a fixed transect within homogeneous habitats according to the Pollard walk protocol (Pollard 1977), producing the number of individuals as abundance data. As our data were not comparable per se with those collected in 1975, we used relative abundances of species to homogenize the data.
Species are named according to Wiemers et al. (2018).

Data analysis
The unavailability of a species/individual's matrix for 1975 allowed us to use only presence/absence non-parametric estimators of species richness to compute sampling completeness. We selected the Chao2 richness estimator (Chao2 Mean) and its standard deviation computed after 100 randomizations of samples running EstimateS 9.1.0 (Colwell 2013). Sampling completeness was computed in terms of proportion of species diversity sampled versus the species diversity estimated. Sampling efficiency was higher than the 94% in the year 1975, being between 84 and 98% in 2004 and between 71 and 94% in 2012. In order to explore the dissimilarities among butterfly communities, different non-metric multidimensional scaling (NMDS) using the Bray-Curtis index was performed. In particular, we analyzed butterfly species composition (quantitative and presence-absence) among the different altitudinal bands in each year. Differences in butterfly communities over the years and altitudinal bands were tested by permutational multivariate analysis of variance using distance matrices (ADONIS function).
To assess which species significantly affected similarity along the altitudinal gradient, we searched for the changes across the years of altitudinal band indicator species. Indicator species are defined as the most characteristic species  Caceres et al. 2010). This is an extension of the original Individual Indicator Value method (IndVal) proposed by Dufrêne and Legendre (1997), which can be used to measure the degree of association between the target species for both individual site groups and each of the possible combinations of site groups. To run statistical analyses, we used the software R 3.6.1 (R Core Team 2019), using the "vegan" package (Oksanen et al. 2014) for NMDS and Permanova analyses and using the "multipatt" function of the package "indicspecies" (De Cáceres and Jansen 2016), with "IndVal.g" as association index on presence/absence data to identify indicator species for each altitudinal band in a multilevel pattern analysis.
We computed the Community Temperature Index (CTI) and the Community Precipitation index (CPI) increasingly in use to measure the responses to climate change of animal communities. They were computed starting from the climatic niche of individual species, determined by abiotic parameters registered in their occurrence sites. In this paper, we used data mined from CLIMBER (Schweiger et al. 2014), where several climatic variables for European butterfly species were reported, at a resolution of 50 km. The scale at which climatic requirements are estimated can affect the possibility of finding a link between population trends and climate change, especially on a regional scale such as in this study, where in the same grid cell of 50 km, we can find alpine prairies and coastal vegetation. Among available climatic variables, we selected the maximum temperature (T.max, mean of monthly maximum temperature) and the minimum precipitation (Prec.min, mean of monthly minimum precipitation) measured across the entire European range of butterflies. The maximum temperature is largely known to determine a rapid range expansion of thermophilous species (Sparks et al. 2005;Battisti et al. 2006), whereas the minimum precipitation can contribute to the decrease of butterfly diversity in a water deficient area such as the Mediterranean Basin (Stefanescu et al. 2011). CTI and CPI are the means of these variables computed according to the relative abundance of any species composing a given community and their formulas can be written as: Prec:min Â Relative:abundance where n is the number of species in the community and the relative abundance is the species' abundance divided by the abundance of all species.
We carried out Pearson correlations between these community climatic indices and altitude for each sampling year, also searching for significant differences of the slopes of the regression lines among years using a t-test with the standard error of the regression models as the error term (R Core Team 2019).
Recently some authors have stated that these indices should be used with caution as they can respond not only to changes in the climate spaces of species but also to some covariates linked to significant changes in land use (Barnagaud et al. 2012). A model-based approach has been proposed for improving community weighted indices (Bowler and Böhning-Gaese 2017). In our case, the bias due to land use change can be considered weaker than those observed in cases of land use changes like those due to the conversion of forest to agricultural lands. This study was conducted in a natural protected territory where we only observed the reduction of cow grazing pressure on vegetation during the study period.

Results
During the 3 years of surveys, we observed a total of 72 species (Table 2). According to the Chao2 richness estimator, the overall richness estimated during 2004 and 2012 was 69.9 ± 5.4 (95% confidence intervals) and 62.2 ± 8.7, respectively, higher than that estimated during 1975 (47.4 ± 1.8).
A total of 46 species was observed in 1975, with a richness in individual sites comprised between 4 (site 21) and 25 (site 3), and a median value of 13.5 species. During 2004, we registered 63 species and 3491 individuals. The lowest richness values increased to 15 species in sites 11, 15, and 19, whilst the highest richness values remained stable (26 species in sites 3, 4, and 5). This resulted in an increase of median richness to 21.5 species. During 2012, the sample was composed of 50 species and 1172 individuals, with the number of species comprised between 9 (site 20) and 26 (site 4), and a median value of 13 species (Fig. 2). The richness registered during 2004 was significantly higher than that registered during 1975 (Kolmogorov-Smirnov test, p = 0.0019) and 2012 (p = 0.0019), whilst no differences were detected between 1975 and 2012 (p = 0.6325).
Seventeen species were registered for only one year, 4 in 1975, 11 in 2004, and 2 in 2012 ( Table 2). The list of species grew during the last research years, i.e., those registered only in both 2004 and 2012, and was longer than those lost since 1975 (Table 2). In fact, 13 new entries appeared to be stable in recent communities and, with the sole exception of Polyommatus bellargus, all of them are more common at lower altitudes than those of the study area.
The composition of butterfly communities was in general quite stable among years, always showing Erebia cassioides and Argynnis niobe among the five most abundant species, whilst Lysandra coridon, Satyrus ferula, and Hipparchia hermione appeared twice among them (Table 2). On the other Table 2 Relative abundance, as percentage, and occupancy, as number of sites where species were recorded, during study years. Trends in relative abundance were recognized as negative or positive (-, +) when the same sign was registered across all years; they were considered strong (--, + +) when their magnitude is higher than the 100% at least between 2 study years, and higher than the 1% in absolute terms between the lowest and the highest value. Trends in occupancy were recognized as negative or positive (-, +) when the same sign was registered across all years. A positive trend was also recognized when a species was absent in the first year and constantly recognized during last two study years. Their changes were considered strong (--, + +) when at least the half of sites were lost or gained compared to the first sampled year and when a given species was present in at least three sites during one of the study years     As observed for relative abundance values, occupancy of species also changed over the years, sometimes corresponding to changes in abundance. In fact, 4 out of 7 species with a strong decrease in occupancy also had a strongly decreased abundance, and only in Erebia cassioides, Hipparchia semele, and Melitaea didyma, the occupancy decreased more than the relative abundance (Table 2). On the other hand, all the 5 species which showed a considerable increase in relative abundance also showed a strong increase in occupancy (Colias crocea, Issoria lathonia, Lycaena hippothoe, Aricia agestis, and Pieris napi), whilst a further 10 species greatly increased only their occupancy (Table 2).
Changes in abundance and distribution of individual species led to changes of indicator species for altitudinal bands. In 1975, the low altitudinal band was the only one with indicator species and they were two mesophilous ones (Table 3). In recent years, marked changes occurred as indicator species of the low altitudinal band were clearly much more thermophilous (Table 3), while new indicator species appeared for the high altitudinal band which previously had been absent.
The communities observed in the different altitudinal bands showed major changes in recent years at community level using both quantitative and presence/absence similarity matrices (Fig. 3). In 1975, no significant differences among communities of the different altitudinal bands were observed, whilst in 2004 and 2012, we observed highly significant differences among them, mainly due to the clear differentiation of the high altitudinal band from the others. In particular, according to quantitative and presence/absence analysis of 2004 (ADONIS, quantitative: p < 0.001; presence/absence: p < 0.01) and 2012 (ADONIS, quantitative: p < 0.001; presence/absence: p < 0.01), species composition of low and medium altitude seemed to be more similar to each other than to high altitude, as they were closer in the Cartesian plane (Fig. 3).
Analyzing data by altitudinal bands, we found that both CTI and CPI were quite stable at low altitude, whilst significant changes of opposite signs occurred in the remaining altitudinal bands. Communities became hotter and dryer at medium altitude, colder and wetter at high altitude (Table 4)

Discussion
In this paper, we found that over the past 37 years, butterfly communities of the Pollino Massif, sampled above 1400 m of altitude, showed substantial changes in species abundance and composition. This has led to a vertical stratification of communities only observed in recent years. This pattern stemmed from opposite trends in community composition over time below and above the tree line (near 1900 m a.s.l.). This was seen as an increasing correlation between Community Temperature Index and Precipitation Community Index with altitude. Despite having only 3 years of censuses, the changes noted since 1975 appear to be qualitatively similar in the two latter censuses (in both 2004 and 2012), indicating that the changes were not driven by a single nor temporary event, but represent true long-term changes in community Fig. 3 Graphic displaying the first two axes of the non-metric multidimensional scaling (NMDS) on butterfly species composition found at different mountain altitudes (low, medium, high) in each year of sampling (1975,2004,2012). Statistical differences in species compositions among the mountain altitude were reported as p-value in the graph (permutational multivariate analysis of variance -ADONIS) composition. This is supported by the fact that both CTI and CPI were significantly different between 1975 and both of the more recent censuses, and not significantly different between the 2004 and 2012 censuses. Across the European latitudinal gradient, it has been observed that limiting factors for butterfly species richness differ, with water availability as the prevailing factor in southern countries (Stefanescu et al. 2011;Herrando et al. 2019) and sunshine hours and temperature in northern ones (Turner et al. 1987). It has been observed that in Mediterranean countries species richness of butterfly communities varied according to observational scale and local characteristics. Wilson et al. (2007) estimated a regional decline of species richness of 90% in the Sierra de Guadarrama, Spain, because the uphill shift of the few low-elevation species presents in that region failed to compensate for the retreat of mountain species at higher altitude. In the Dadia National Park, Greece, the weak altitudinal gradient of this area (20-650m) favored a decrease of high-elevation species and an increase of low-elevation species (Zografou et al. 2014). In Alpine sites, Cerrato et al. (2019) observed an increase in species richness over a short time period, similar to the one observed in northern Scandinavia over a long time period (Franzén and Öckinger 2012). Scarce evidence of a similar trend also arose from our study as species richness was higher in 2004than in 1975. However, 1975 and 2012 seemed to be more similar than 1975 and 2004, as during 2012, median and maximum values of diversity were the same of those registered in 1975, showing an increase only in minimum values. This could be due to our study design, which is affected by a phenological sampling bias. Although it is known that phenological changes occur as a response to climate warming, shifting the presence in the field of the winged stages (Roy and Sparks 2000;Bell et al. 2019), our sampling only covered 2 months of the butterfly activity season in order to cover the same time span as Balletto et al. (1977). A truncated monitoring season can produce data affected by differences in the phenology of species between years due to different climatic conditions. In fact, low and median values of species richness were significantly higher during 2004, when it was likely to have increased due to the presence of species usually flying earlier in the study area such as Boloria euphrosyne, Cupido minimus, and Melitaea cinxia. It is also interesting to note that in 2003, a major heat wave occurred, likely affecting the butterfly communities observed during 2004, e.g., favoring the observed outbreak of Lysandra coridon. In any case, although 2004 appears to be anomalous from quantitative (about half of the sample composed by L. coridon specimens) and qualitative (highest species richness) viewpoints, we found very similar results to those in 2012 In this study, we found an evident change in community compositions across years, which led to a significant vertical stratification absent in 1975, observed firstly in 2004, and successively confirmed in 2012. This was the most interesting finding of this study. Results of abundance-and presence/ absence-based analyses are consistent, allowing us to exclude the possibility that observed patterns could be significantly affected by individual species outbreaks and or by the phenological shift of the flying period of few and not abundant species. In fact, presence/absence-based analyses do not take into account abundance of species, making these results unbiased by anomalous abundances of one or few species, whilst abundance based analyses are only slightly modified by the presence of few rare species in flight later than usual in a given period.
The vertical stratification of communities seems to be due to two concurrent changes in medium and low altitudinal bands: the increase of generalist and widespread species (Warren et al. 2001;Menéndez et al. 2007) and the reduction of cold-adapted ones (Zografou et al. 2014). In fact, most of the species having a positive trend of abundance and occupancy in the Pollino Massif are generalist ones, such as Colias crocea, Issoria lathonia, Aricia agestis, Pieris napi, and Papilio machaon. However, their increase is not homogeneous in all altitudinal bands. For example, the last three species have appeared in recent years only in low and medium altitude communities, with A. agestis and P. napi representing indicator species for the low altitudinal band. This likely happens because of the persistence of extreme environmental conditions above the treeline that confers resistance to the colonization of these species. On the other hand, coldadapted species suffer from the increase of temperature only at low altitude as demonstrated by the retreat of Erebia cassioides towards the top of the Pollino Massif where it increased its populations (Scalercio et al. 2014) and became the indicator species for the high altitudinal band.
Following the creation of the Pollino National Park in 1993, grazing pressure by cattle has been reduced. In other systems, this has been shown to lead to higher meadow vegetation and woody encroachment into meadow areas (Carmel and Kadmon 1999;Sternberg et al. 2000). In the Swiss Alps, the abandonment of hay-making and grazing on the mountaintops led to an upward shift in treeline, reclaiming historically forested lands, and dominating any effect of climate change (Gehrig-Fasel et al 2007). Higher meadow vegetation and woody encroachment can lead to local cooling and an increase in available moisture. For example, Körner (2003), in a European-wide study of alpine habitats, found that soil temperatures at 5-cm depth was an average of 2°C warmer 200-300 m above treeline than at treeline itself. If a similar process were occurring in the Pollino Massif, then the mountaintop alpine habitats would have undergone local cooling and become wetter. This would be consistent with the lowing of CTI and increase in CPI at the mountaintop sites in our study. However, there is a need of further studies to test this hypothesis. Unfortunately, historical data is lacking, but using, indirect evidence, such as the variation of the number of grazing cows across years, or comparisons of old and recent photographs of the same place may be possible in future studies. This possibility highlights the potential importance of expanding natural protected areas to help buffer wild mountaintop species against the detrimental effects of temperature increase in the Mediterranean Basin.

Conclusions
In the Mediterranean Basin, our study, coupled with others, demonstrate that butterfly communities have already been strongly modified by recent climate change and that these changes are especially prominent in mountainous areas. In the Pollino Massif, we found asymmetrical trends in butterfly communities along altitudinal gradients, with lower and mid-slope communities shifting towards hotter/drier associated species, and mountaintop communities shifting towards cooler/wetter associated species. Additionally, butterfly communities were more stable above treeline than below, in spite of general regional warming, possibly due to buffering effects of land-use changes, but this hypothesis could not be tested.
Author contribution SS conceived and designed the experiments. SS performed the experiments. LS, AI, and SS analyzed the data. SS, AI, and LS wrote the manuscript; all authors provided editorial advice.
Funding Open access funding provided by Università degli Studi di Napoli Federico II within the CRUI-CARE Agreement. Research supported by the Ministry of University and Education of Italy (PRIN D.M. 19 marzo 2010, n. 51, prot. 200947YRB9 Impact of global change on ecosystems, animal communities, and species of alpine and Mediterranean areas of Italy: models, scenarios, and evaluation from macro-scale to microscale, based on ecology and philogeography of vertebrates and invertebrates).
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.