“…Following Long et al ( 2021b ), we used six bioclimatic variables at a 10 km resolution averaged for the period 1979–2013 from the CHELSA database (Karger et al, 2017 ) for model training and projection of the current potential distribution of D. involucrate : annual mean temperature (BIO1), isothermality (BIO3), temperature annual range (BIO7), precipitation of the driest month (BIO14), precipitation seasonality (BIO15), and precipitation of the warmest quarter (BIO18), for these variables have low multicollinearity (all variables with Pearson correlation coefficients | r | < .7 and variance inflation factor VIF < 5) and have the greatest ecological relevance to D. involucrate (Long et al, 2021b ). Similarly, the same six bioclimatic variables at a 2.5 arc‐min resolution for two future time periods, 2050s (averaged for 2041–2060) and 2070s (averaged for 2061–2080), under two representative concentration pathways (RCPs) scenarios, RCP2.6 and RCP8.5, from six global circulation models (GCMs): CNRM‐CM6‐1, CNRM‐ESM2‐1, CanESM5, IPSL‐CM6A‐LR, MIROC‐ES2L, and MIROC6 were obtained from the WorldClim database (Fick & Hijmans, 2017 ).…”