Workflow used for species distribution modelling. Presence data was transformed to a raster grid (5 × 5 km), spatiotemporally analysed, accordingly grouped into five levels of persistence and then combined with environmental data and generated absences (using four different sampling strategies) to fit logistic regression models. The absence sampling strategy that derived the best evaluation measures was used for predicting suitability under current (1979–2013) and future climatic conditions after going through a model selection and averaging procedure.