“…Nineteen variables were considered for the habitat suitability (4500 m) models: 14 bioclimatic variables (WorldClim database version 2.1; [56]), 1 anthropogenic (population density), and 4 environmental (20-year mean normalized difference vegetation index (NDVI), soil bulk density, soil texture, and slope). We excluded 5 of the 19 available WorldClim bioclimatic variables from the analysis based on preliminary tests of collinearity, consideration of recent variables included in ecological niche models for the European ground squirrel and its congener, namely, the Anatolian ground squirrel (Spermophilus xanthoprymnus) [43,44], and an emphasis on seasonal mean or range rather than min-max values. Details of the source, initial resolution, and model for which a variable was considered are provided in Table 1.…”