Appendix A Supplementary data

Transcription

1 Appendix A Supplementary data Table A1 Data sources used to document Cabrera vole occurrences ( ) 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: 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.: 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): 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), 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): 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 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 Central Central 3 South-eastern Spain North-eastern North-eastern

2 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 SNB46 RI SNB62 RI SNB65 RE SNB78 RE SNB81 RE SNB90 RI SNC55 RI SNC75 RE SNC76 RI SNC80 RE SNC93 RE SND01 FE SND24 FE SND44 FI SND79 FE SND84 FI SPC08 RI SPC27 RE SPC29 RI SPC46 RE SPD00 RI SPD05 FI SPD23 RI SPD24 FI SPD29 FE SPD32 RI SPD36 FE SPD40 RE SPE42 FI SPE50 FI TPE24 FE TPE58 FE TPF33 FE TPF61 FI TPF65 FI TPF77 FI TPG43 FE TPG70 FI TPG71 FE

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 bio bio bio bio bio bio bio bio bio bio bio bio bio bio bio bio bio bio

4 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 = 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.

5 Figure A1 Sampled 10km x10km UTM grid squares. Inside the potential niche (black) and outside the potential niche (white), as given by ENM1.

6 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).