Our work aims to assess how butterfly communities in the Italian Maritime Alps changed over the past 40 years, in parallel with altitudinal shifts occurring in plant communities. In 2019, we sampled butterflies at 7 grassland sites, between 1300–1900 m, previously investigated in 2009 and 1978, by semi-quantitative linear transects. Fine-scale temperature and precipitation data elaborated by optimal interpolation techniques were used to quantify climate changes. The changes in the vegetation cover and main habitat alterations were assessed by inspection of aerial photographs (1978–2018/1978–2006–2015). The vegetation structure showed a marked decrease of grassland habitats and an increase of woods (1978–2009). Plant physiognomy has remained stable in recent years (2009–2019) with some local exceptions due to geomorphic disturbance. We observed butterfly ‘species substitution’ indicating a general loss in the more specialised and a general gain in more tolerant elements. We did not observe any decrease in species richness, but rather a change in guild compositions, with (i) an overall increased abundance in some widespread and common lowland species and (ii) the disappearance (or strong decrease) of some alpine (high elevation) species, so that ‘resilience’ could be just delusive. Changes in butterfly community composition were consistent with predicted impacts of local warming.

Mountain biodiversity is associated with rare and fragile biota that are highly sensitive to climate change. To estimate the vulnerability of biodiversity to temperature rise, long-term field data are crucial. Species distribution models are an essential tool, in particular for invertebrates, for which detailed information on spatial and temporal distributions is largely missing. We applied presence-only distribution models to field data obtained from a systematic survey of 5 taxa (birds, butterflies, carabids, spiders, staphylinids), monitored in the northwestern Italian Alps. We estimated the effects of a moderate temperature increase on the multi-taxa distributions. Only small changes in the overall biodiversity patterns emerged, but we observed significant differences between groups of species and along the altitudinal gradient. The effects of temperature increase could be more pronounced for spiders and butterflies, and particularly detrimental for high-altitude species. We observed significant changes in community composition and species richness, especially in the alpine belt, but a clear separation between vegetation levels was retained also in the warming scenarios. Our conservative approach suggests that even a moderate temperature increase (about 1 °C) could influence animal biodiversity in mountain ecosystems: only long-term field data can provide the information to improve quantitative predictions, allowing us to readily identify the most informative signals of forthcoming changes.

Mountain ecosystems are biodiversity hotspots that are increasingly threatened by climate and land use/land cover changes. Long-term biodiversity monitoring programs provide unique insights into resulting adverse impacts on plant and animal species distribution. Species distribution models (SDMs) in combination with satellite remote sensing (SRS) data offer the opportunity to analyze shifts of species distributions in response to these changes in a spatially explicit way. Here, we predicted the presence probability of three different rove beetles in a mountainous protected area (Gran Paradiso National Park, GPNP) using environmental variables derived from Landsat and Aster Global Digital Elevation Model data and an ensemble modelling approach based on five different model algorithms (maximum entropy, random forest, generalized boosting models, generalized additive models, and generalized linear models). The objectives of the study were (1) to evaluate the potential of SRS data for predicting the presence of species dependent on local-scale environmental parameters at two different time periods, (2) to analyze shifts in species distributions between the years, and (3) to identify the most important species-specific SRS predictor variables. All ensemble models showed area under curve (AUC) of the receiver operating characteristics values above 0.7 and true skills statistics (TSS) values above 0.4, highlighting the great potential of SRS data. While only a small proportion of the total area was predicted as highly suitable for each species, our results suggest an increase of suitable habitat over time for the species Platydracus stercorarius and Ocypus ophthalmicus, and an opposite trend for Dinothenarus fossor. Vegetation cover was the most important predictor variable in the majority of the SDMs across all three study species. To better account for intra- and inter-annual variability of population dynamics as well as environmental conditions, a continuation of the monitoring program in GPNP as well as the employment of SRS with higher spatial and temporal resolution is recommended.