The influence of climate on the distribution of taxa has been extensively investigated in the last two decades through Habitat Suitability Models (HSMs). In this context, the Worldclim database represents an invaluable data source as it provides worldwide climate surfaces for both historical and future time horizons. Thousands of HSMs‐based papers have been published taking advantage of Worldclim 1.4, the first online version of this repository. In 2017, Worldclim 2.1 was released. Here, we evaluated spatially explicit prediction mismatch at continental scale, focusing on Europe, between HSMs fitted using climate surfaces from the two Worldclim versions (between‐version differences). To this aim, we simulated occurrence probability and presence‐absence across Europe of four virtual species (VS) with differing climate‐occurrence relationships. For each VS, we fitted HSMs upon uncorrelated bioclimatic variables derived from each Worldclim version at three grid resolutions. For each factor combination, HSMs attaining sufficient discrimination performance on spatially independent test data were projected across Europe under current conditions and various future scenarios, and importance scores of the single variables were computed. HSMs failed in accurately retrieving the simulated climate‐occurrence relationships for the climate‐tolerant VS and the one occurring under a narrow combination of climatic conditions. Under current climate, noticeable between‐version prediction mismatch emerged across most of Europe for these two VSs, whose simulated suitability mainly depended upon diurnal or yearly variability in temperature; differently, between‐version differences were more clustered toward areas showing extreme values, like mountainous massifs or southern regions, for VSs responding to average temperature and precipitation trends. Under future climate, the chosen emission scenarios and Global Climate Models did not evidently influence between‐version prediction discrepancies, while grid resolution synergistically interacted with VSs' niche characteristics in determining extent of such differences. Our findings could help in re‐evaluating previous biodiversity‐related works relying on geographical predictions from Worldclim‐based HSMs.

Boosted Regression Trees (BRT) is one of the modelling techniques most recently applied to biodiversity conservation and it can be implemented with presence-only data through the generation of artificial absences (pseudo-absences). In this paper, three pseudo-absences generation techniques are compared, namely the generation of pseudo-absences within target-group background (TGB), testing both the weighted (WTGB) and unweighted (UTGB) scheme, and the generation at random (RDM), evaluating their performance and applicability in distribution modelling and species conservation. The choice of the target group fell on amphibians, because of their rapid decline worldwide and the frequent lack of guidelines for conservation strategies and regional-scale planning, which instead could be provided through an appropriate implementation of SDMs. Bufo bufo, Salamandrina perspicillata and Triturus carnifex were considered as target species, in order to perform our analysis with species having different ecological and distributional characteristics. The study area is the “Gran Sasso—Monti della Laga” National Park, which hosts 15 Natura 2000 sites and represents one of the most important biodiversity hotspots in Europe. Our results show that the model calibration ameliorates when using the target-group based pseudo-absences compared to the random ones, especially when applying the WTGB. Contrarily, model discrimination did not significantly vary in a consistent way among the three approaches with respect to the tree target species. Both WTGB and RDM clearly isolate the highly contributing variables, supplying many relevant indications for species conservation actions. Moreover, the assessment of pairwise variable interactions and their three-dimensional visualization further increase the amount of useful information for protected areas’ managers. Finally, we suggest the use of RDM as an admissible alternative when it is not possible to individuate a suitable set of species as a representative target-group from which the pseudo-absences can be generated.

Climate is often considered as a key ecological factor limiting the capability of expansion of most species and the extent of suitable habitats. In this contribution, we implement Species Distribution Models (SDMs) to study two parapatric amphibians, Lissotriton vulgaris meridionalis and L. italicus, investigating if and how climate has influenced their present and past (Last Glacial Maximum and Holocene) distributions. A database of 901 GPS presence records was generated for the two newts. SDMs were built through Boosted Regression Trees and Maxent, using the Worldclim bioclimatic variables as predictors. Precipitation-linked variables and the temperature annual range strongly influence the current occurrence patterns of the two Lissotriton species analyzed. The two newts show opposite responses to the most contributing variables, such as BIO7 (temperature annual range), BIO12 (annual precipitation), BIO17 (precipitation of the driest quarter) and BIO19 (precipitation of the coldest quarter). The hypothesis of climate influencing the distributions of these species is also supported by the fact that the co-occurrences within the sympatric area fall in localities characterized by intermediate values of these predictors. Projections to the Last Glacial Maximum and Holocene scenarios provided a coherent representation of climate influences on the past distributions of the target species. Computation of pairwise variables interactions and the discriminant analysis allowed a deeper interpretation of SDMs’ outputs. Further, we propose a multivariate environmental dissimilarity index (MEDI), derived through a transformation of the multivariate environmental similarity surface (MESS), to deal with extrapolation-linked uncertainties in model projections to past climate. Finally, the niche equivalency and niche similarity tests confirmed the link between SDMs outputs and actual differences in the ecological niches of the two species. The different responses of the two species to climatic factors have significantly contributed to shape their current distribution, through contractions, expansions and shifts over time, allowing to maintain two wide allopatric areas with an area of sympatry in Central Italy. Moreover, our SDMs hindcasting shows many concordances with previous phylogeographic studies carried out on the same species, thus corroborating the scenarios of potential distribution during the Last Glacial Maximum and the Holocene emerging from the models obtained.

The pond turtle Emys trinacris is an endangered endemic species of Sicily showing a fragmented distribution throughout the main island. In this study, we applied “Ensemble Niche Modelling”, combining more classical statistical techniques as Generalized Linear Models and Multivariate Adaptive Regression Splines with machine-learning approaches as Boosted Regression Trees and Maxent, to model the potential distribution of the species under current and future climatic conditions. Moreover, a “gap analysis” performed on both the species’ presence sites and the predictions from the Ensemble Models is proposed to integrate outputs from these models, in order to assess the conservation status of this threatened species in the context of biodiversity management. For this aim, four “Representative Concentration Pathways”, corresponding to different greenhouse gases emissions trajectories were considered to project the obtained models to both 2050 and 2070. Areas lost, gained or remaining stable for the target species in the projected models were calculated. E. trinacris’ potential distribution resulted to be significantly dependent upon precipitation-linked variables, mainly precipitation of wettest and coldest quarter. Future negative effects for the conservation of this species, because of more unstable precipitation patterns and extreme meteorological events, emerged from our analyses. Further, the sites currently inhabited by E. trinacris are, for more than a half, out of the Protected Areas network, highlighting an inadequate management of the species by the authorities responsible for its protection. Our results, therefore, suggest that in the next future the Sicilian pond turtle will need the utmost attention by the scientific community to avoid the imminent risk of extinction. Finally, the gap analysis performed in GIS environment resulted to be a very informative post-modeling technique, potentially applicable to the management of species at risk and to Protected Areas’ planning in many contexts.