Distribution models representing areas of climate suitability and predicted diversity indices over time. (a) Values of the response curves shown are averages over 60 runs; grey margins show standard deviation over all replicates. (b) Relative importance that evaluated using AUC, with standard deviation calculated for each variable. (c) Relative suitability under current climate and the two emission models (i.e., RCP 2.6 and RCP 8.5) in 2070. The color spectrum represents suitability on the maps, from areas with low (0) to high (1) values. (d) Comparison of diversity indices over time. The predictions are grouped by different data subsets using sample locations overlapping with moderately suitable (P > 0.5) and highly suitable (P > 0.75) habitats. The error bars denote the standard deviation for each index. AR, allelic richness; FIS, inbreeding coefficient; He, expected heterozygosity.