Appendix A Supplementary data Table A1 Data sources used to document Cabrera vole occurrences (1970-2011) in UTM 10 km x 10 km grid squares across the Iberian Peninsula. Species occurrence in some squares were recorded by more than one study, and so the number of grid squares summed across studies is larger than the total number of grid squares with confirmed occurrences. Source Published Inventario Nacional de Biodiversidad. Vertebrados. url: http://www.marm.es/es/biodiversidad/temas/inventarios-nacionales/inventario-nacionaldebiodiversidad/mamiferos_m.aspx. Conservation Biology Unit Database (Universidade de Évora, ) Mira, A.; Marques, C.C.; Santos, S.; Rosário, I.T.; & Mathias, M.L. (2008). Environmental determinants of the distribution of the Cabrera vole (Microtus cabrerae) in : Implications for conservation. Mammalian Biology - Zeitschrift für Säugetierkunde, 73(2), pp.: 102-110. Santos, S. M., Mathias, M., Mira, A., & Simões, M. P. (2007). Vegetation structure and composition of road verge and meadow sites colonized by Cabrera vole (Microtus cabrerae Thomas). Polish Journal of Ecology, 55(3), 481. Santos, S. M., Simões, M. P., Mathias, M. D. L., & Mira, A. (2006). Vegetation analysis in colonies of an endangered rodent. Ecological Research, 21(2): 197-207. Rosário, I. T., Cardoso, P. E., & Mathias, M. D. L. (2008). Is habitat selection by the Cabrera vole (Microtus cabrerae) related to food preferences? Mammalian Biology-Zeitschrift für Säugetierkunde, 73(6), 423-429. Rosário. I.T. Towards a conservation strategy for na endangered rodent, the Cabrera vole (Microtus cabrerae Thomas) insights from ecological data. PhD Thesis. Faculdade de Ciências, Departamento de Biologia Animal, Universidade de Lisboa. Lisboa,. Pita, R. (2010). Persistence and coexistence of spatially structured populations in heterogeneous environments: The case of Cabrera and water voles in Mediterranean farmland. PhD Thesis, Universidade de Évora. Évora,. Mira, A.; Ascensão, F. & Alcobia, S. (2003). Distribuição das espécies de roedores e insectívoros Relatório final, Unidade de Biologia da Conservação, Universidade de Évora. Gonçalves P., Alcobia S. & Santos-Reis M., Eds. (2013). Atlas dos Mamíferos na Charneca do Infantado. Companhia das Lezírias S.A. / Centro de Biologia Ambiental (FCUL), Benavente e Lisboa, 92 pp. Ortuño, A. (2009). Nuevos datos sobre la distribución del Topillo de Cabrera Microtus cabrerae Thomas, 1906 en Murcia. Galemys, 21(2): 71-72. Profico Ambiente (2009). RECAPE do Aproveitamento Hidroeléctrico de Foz Tua; Anexo Y - Estudo Complementar Sobre o Rato-de-Cabrera; Relatório de Conformidade Ambiental do Projecto de Execução do Aproveitamento Hidroeléctrico de Foz Tua. Mathias, M.L.; Mira, A.; Pereira, M.; Pereira, P.; Nunes, A.C.; Marques, C.C.; Figueiredo, C.; Carvalho, F.N.; Sousa, I.; Perestrello, M.; Santos, M.J.; Santos, S. & Oliveira, V. (2004). Programa de monitorização do património natural (Área do regolfo de Alqueva e Pedrógão) - Projecto Pmo 6.2 - Monitorização de roedores - Relatório Final, Centro de Biologia Ambiental (Faculdade de Ciências da Universidade de Lisboa) e Centro de Ecologia Aplicada (Universidade de Évora). Anonym (2008).Estudo de Impacte Ambiental do Estudo Prévio da Ocupação Turística da UNOP 4 DE Tróia, Resumo Não Técnico. Unpublished Frederico Mestre (presonal observations). Célia Gomes unpublished data. Geographical range covered Continental Spain Continental Continental Nr. grid squares 276 51 35 Central 25 15 10 6 Central 3 South-eastern Spain North-eastern North-eastern 3 2 1 1 3 3
Table A2 - Summary results of non-invasive genetic analysis. Grid square location: RE Rear edge, out of the predicted potential niche; RI - Rear edge, inside the predicted potential niche; FE Front edge, out of the predicted potential niche; FI - Front edge, inside the predicted potential niche. See the main text for details. UTM code Grid square location Sample collection Genetic analysis Cabrera vole Other species Genetic ID failed Sites Samples Sites Samples Sites Samples Sites Samples 29SNB11 RE 0 0 - - - - - - - 29SNB46 RI 1 1 1 1 1 1 0 0 0 29SNB62 RI 1 2 1 1 1 1 0 0 0 29SNB65 RE 0 0 - - - - - - - 29SNB78 RE 0 0 - - - - - - - 29SNB81 RE 0 0 - - - - - - - 29SNB90 RI 0 0 - - - - - - - 29SNC55 RI 1 1 1 1 0 0 0 0 1 29SNC75 RE 1 1 1 1 0 0 1 1 0 29SNC76 RI 1 1 1 1 0 0 0 0 1 29SNC80 RE 0 0 - - - - - - - 29SNC93 RE 1 1 1 1 0 0 1 1 0 29SND01 FE 2 6 2 6 0 0 2 6 0 29SND24 FE 5 12 5 10 0 0 5 10 0 29SND44 FI 2 10 2 6 1 3 1 2 1 29SND79 FE 9 27 6 8 3 3 4 5 0 29SND84 FI 1 2 1 2 1 2 0 0 0 29SPC08 RI 1 1 1 1 1 1 0 0 0 29SPC27 RE 0 0 - - - - - - - 29SPC29 RI 1 2 1 2 1 2 0 0 0 29SPC46 RE 0 0 - - - - - - - 29SPD00 RI 1 2 1 2 1 2 0 0 0 29SPD05 FI 1 3 1 2 1 2 0 0 0 29SPD23 RI 0 0 - - - - - - - 29SPD24 FI 1 3 1 3 0 0 1 3 0 29SPD29 FE 8 22 2 2 2 2 0 0 0 29SPD32 RI 0 0 - - - - - - - 29SPD36 FE 2 5 1 2 1 2 0 0 0 29SPD40 RE 1 1 1 1 0 0 1 1 0 29SPE42 FI 2 11 1 3 1 3 0 0 0 29SPE50 FI 4 12 1 2 1 2 0 0 0 29TPE24 FE 6 18 5 8 2 2 3 5 1 29TPE58 FE 10 46 10 12 0 0 9 10 2 29TPF33 FE 9 44 2 3 1 1 1 1 1 29TPF61 FI 5 20 1 3 1 3 0 0 0 29TPF65 FI 3 12 1 2 1 2 0 0 0 29TPF77 FI 7 13 2 2 2 2 0 0 0 29TPG43 FE 8 33 7 12 0 0 7 11 1 29TPG70 FI 7 20 1 2 1 1 1 1 0 29TPG71 FE 8 19 7 13 0 0 6 10 3
Table A3 Summary results of pairwise Spearman correlation analysis between climatic variables. Large ( r > 0.7) pairwise correlations are highlighted in dark grey, and variables retained for subsequent distribution modelling are highlighted in light grey. bio1 - Annual Mean Temperature, bio2 - Mean Diurnal Range (Mean of monthly (max temp - min temp)), bio3 - Isothermality (bio2/bio7) (* 100), bio4 - Temperature Seasonality (standard deviation *100), bio5 - Max Temperature of Warmest Month, bio6 - Min Temperature of Coldest Month, bio7 - Temperature Annual Range (bio5-bio6), bio8 - Mean Temperature of Wettest Quarter, bio9 - Mean Temperature of Driest Quarter, bio10 - Mean Temperature of Warmest Quarter, bio11 - Mean Temperature of Coldest Quarter, bio12 - Annual Precipitation, bio13 - Precipitation of Wettest Month, bio14 - Precipitation of Driest Month, bio15 - Precipitation Seasonality (Coefficient of Variation), bio16 - Precipitation of Wettest Quarter, bio17 - Precipitation of Driest Quarter, bio18 - Precipitation of Warmest Quarter, bio19 - Precipitation of Coldest Quarter. bio2 bio3 bio4 bio5 bio6 bio7 bio8 bio9 bio10 bio11 bio12 bio13 bio14 bio15 bio16 bio17 bio18 bio19 bio1 0.068 0.162-0.050 0.695 0.881 0.008 0.426 0.814 0.894 0.938-0.222-0.037-0.833 0.785-0.004-0.798-0.756 0.066 bio2-0.042 0.848 0.715-0.334 0.954-0.043 0.478 0.431-0.207-0.594-0.631-0.377 0.066-0.588-0.427-0.467-0.427 bio3-0.515-0.101 0.307-0.316-0.100 0.048-0.084 0.323 0.208 0.284-0.123 0.276 0.306-0.094-0.097 0.353 bio4 0.627-0.447 0.962 0.026 0.343 0.375-0.355-0.606-0.674-0.233-0.100-0.650-0.291-0.324-0.544 bio5 0.323 0.699 0.211 0.886 0.923 0.434-0.571-0.488-0.822 0.557-0.429-0.834-0.839-0.264 bio6-0.407 0.304 0.552 0.601 0.982 0.130 0.308-0.600 0.743 0.333-0.538-0.497 0.361 bio7-0.030 0.432 0.419-0.295-0.617-0.676-0.317-0.015-0.638-0.375-0.414-0.492 bio8 0.165 0.382 0.382-0.499-0.398-0.206 0.008-0.458-0.216-0.137-0.571 bio9 0.914 0.637-0.375-0.250-0.867 0.699-0.193-0.854-0.852-0.021 bio10 0.694-0.476-0.336-0.869 0.671-0.295-0.861-0.835-0.187 bio11 0.004 0.194-0.679 0.765 0.219-0.626-0.579 0.254 bio12 0.953 0.370 0.045 0.952 0.448 0.418 0.875 bio13 0.210 0.241 0.986 0.291 0.281 0.892 bio14-0.862 0.143 0.987 0.978-0.018 bio15 0.301-0.816-0.806 0.412 bio16 0.227 0.208 0.939 bio17 0.988 0.059 bio18 0.011
Table A4 Summary of BIOMOD modelling options used for running each of the eight algorithms considered in the ensemble to predict current and future distribution range of the Cabrera vole Model algorithm Code Description of modelling options adopted in BIOMOD Generalised linear model Generalised boosting model Generalised additive model Classification tree analysis Artificial neural networks Flexible discriminant analysis Multivariate adaptive regression spline GLM GBM GAM CTA ANN FDA MARS TypeGLM = "poly", Test = "AIC" - ordinary polynomial terms (second and third order) were used, with AIC as the stepwise selection procedure for the selection of the most significant variables to be included in the model (Akaike, 1974). No.trees = 4000 - GBM was run with 4000 trees, which can be considered a good compromise value (Thuiller et al. 2009a), and with 5 cross validations to select the optimal number of trees (by default). Spline = 4 - GAM in BIOMOD uses a cubic spline smoother, which is set of polynomials of degree less than or equal to 3 (Thuiller et al. 2009a). It was run with 4 degrees of freedom (by default). CV.tree = 30 - CTA was run with a 30-fold cross-validation to select the optimal length of the tree (which represents a trade-off between number of leaves and the explained deviance). CV.ann = 3 - ANN was run with the default option of 3 cross validations to select the best amount of weigh decay and the number of units in the hidden layer. BIOMOD s default options. BIOMOD s default options. Random forest RF BIOMOD s default options.
Figure A1 Sampled 10km x10km UTM grid squares. Inside the potential niche (black) and outside the potential niche (white), as given by ENM1.
Figure A2 Maps of predicted distribution of the Cabrera vole in 2020, 2050, and 2080 (IPCC Scenarios B2 and A1b), based on the projecti on of models developed considering the original dataset (ENM1) and the dataset including additional field data (ENM2).