ZooKeys 341: 49–76, doi: 10.3897/zookeys.341.5502
Diversity and distribution of reptiles in Romania
Dan Cogălniceanu 1, Laurentiu Rozylowicz 2, Paul Székely 1, Ciprian Samoilă 1, Florina Stănescu 1, Marian Tudor 1, Diana Székely 1, Ruben Iosif 1
1 University Ovidius Constanţa, Faculty of Natural Sciences and Agricultural Sciences, Al. Universităţii nr. 1, corp B, 900470, Constanţa, Romania
2 University of Bucharest, Center for Environmental Research and Impact Studies, Bd. N. Bălcescu nr. 1, 010041, Bucharest, Romania

Corresponding author: Laurenţiu Rozylowicz (laurentiu.rozylowicz@g.unibuc.ro)

Academic editor: J. Penner

received 11 May 2013 | accepted 29 September 2013 | Published 8 October 2013


(C) 2013 Dan Cogălniceanu. This is an open access article distributed under the terms of the Creative Commons Attribution License 3.0 (CC-BY), which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.


For reference, use of the paginated PDF or printed version of this article is recommended.

Citation: Cogălniceanu D, Rozylowicz L, Székely P, Samoilă C, Stănescu F, Tudor M, Székely D, Iosif R (2013) Diversity and distribution of reptiles in Romania. ZooKeys 341: 49–76. doi: 10.3897/zookeys.341.5502

Abstract

The reptile fauna of Romania comprises 23 species, out of which 12 species reach here the limit of their geographic range. We compiled and updated a national database of the reptile species occurrences from a variety of sources including our own field surveys, personal communication from specialists, museum collections and the scientific literature. The occurrence records were georeferenced and stored in a geodatabase for additional analysis of their spatial patterns. The spatial analysis revealed a biased sampling effort concentrated in various protected areas, and deficient in the vast agricultural areas of the southern part of Romania. The patterns of species richness showed a higher number of species in the warmer and drier regions, and a relatively low number of species in the rest of the country. Our database provides a starting point for further analyses, and represents a reliable tool for drafting conservation plans.

Keywords

Reptilia, species distribution, species range, biodiversity data, species richness, rarity

Introduction

Reptiles are declining worldwide at an alarming rate (Gibbons et al. 2000, Böhm et al. 2013) and, along with amphibians, are considered among the most threatened vertebrate groups (Stuart et al. 2008, Hof et al. 2011). The decline of reptiles has been induced by a variety of threats such as habitat loss, degradation and fragmentation, pet trade, invasive species, pollution and diseases (Gibbons et al. 2000, Cox and Temple 2009). Currently, it is considered that climate changes severely affect reptiles (Böhm et al. 2013). Therefore, we need to assess the potential impacts of the forecasted global changes by assessing anticipated changes in the ecological niche of species. Although complex species distribution models have been developed, with some managing to predict the species ecological niches from just a few known localities (Phillips et al. 2006, Peterson et al. 2011), we still face difficulties to accurately assess the impact of these future threats. Among the limitations encountered when modeling ecological niches is the quality of the biological data used to train the models, which is the first condition to meet certain standards (Pineda and Lobo 2009). Data quality issues are triggered by use of improper spatial resolution, misidentification, missing data, and the lack of a proper sampling design (Araújo and Guisan 2006, Peterson et al. 2011). In this manner, it is mandatory to check for any bias in the sampling effort when mapping species distribution (e.g., Proess 2003, Loureiro et al. 2010).

Romania has five biogeographic regions (i.e., Alpine, Continental, Pannonian, Black Sea, and Steppic) out of the nine regions recognized by the European Union (Cogălniceanu and Cogălniceanu 2010, Iojă et al. 2010, Evans 2012). Romania still maintains a significant proportion of habitats with high conservation value, such as boreal coniferous forests, mesophilous, hygrophilous and xerothermic broadleaved forests, grassland and shrubby ecosystems. There is also a rich diversity of aquatic ecosystems including mountain springs and rivers, river floodplains, glacial lakes, coastal wetlands and bogs (Rey et al. 2007).

The diversity of reptiles is surprisingly high for a country with a mostly temperate and continental climate. There are 23 reptile species in Romania, out of which 12 species reach the limit of their geographic range (i.e., Testudo hermanni, Testudo graeca, Ablepharus kitaibelii, Lacerta trilineata, Podarcis tauricus, Podarcis muralis, Eremias arguta, Darevskia praticola, Eryx jaculus, Dolichophis caspius, Elaphe sauromates, Vipera ammodytes), while other two species are near their range edge in Romania (i.e., Zamenis longissimus, Vipera ursinii) (Gasc et al. 1997).

No updated distribution maps have been published since the publication of the milestone volume on Reptiles in the series Fauna of Romania fifty-two years ago (Fuhn and Vancea 1961), despite a substantial increase in the inventory effort over the following years. However, there are published papers compiling species distributions of single species such as Darevskia praticola (Sós et al. 2012), Vipera ursinii (Krecsák and Zamfirescu 2008), Natrix tessellata (Strugariu et al. 2011), Emys orbicularis (Sós 2011), Testudo hermanni (Rozylowicz and Dobre 2010), Eryx jaculus (Krecsák and Iftime 2006), or restricted to particular geographical areas (e.g., Ghira et al. 2002, Covaciu-Marcov et al. 2006).

We here present an updated overview of the distribution and diversity of reptiles in Romania. The aims of our study are to (1) map the distribution of reptile species in Romania and (2) analyze the spatial pattern of distribution data.

Materials and methods
Mapping species occurrences

We extracted the occurrence records from four major sources: published data, museum collections, personal communications from specialists, and our own unpublished field data. The records were primarily stored and managed in a Microsoft Access database, and later imported it in an ESRI file geodatabase using ArcGIS Desktop 10.1 (ESRI, CA). We checked for data quality by (1) filtering the database for doubtful and erroneous records, (2) aggregating the known localities to a finer resolution, and (3) assessing the bias in sampling effort. Our own data were collected over a period of almost 25 years and it involved a large variety of methods. Since the majority of studies carried out were of ecology, the detailed distribution data was not made available in the resulting publications. Our own few publications presenting species distribution data were included in Appendix I. No voucher specimens were collected during our studies.

The distribution records that could not be georeferenced to an actual locality or toponym (e.g., occurrences assigned to mountain ranges, geographical provinces or hydrographic basins) or records with unspecified taxa within genera were not included in the geodatabase. Other doubtful or erroneous records such as species out of their known range or vagrant individuals sensu IUCN (2001) were also discarded.

The species taxonomy considered in the present paper is based on Speybroeck et al. (2010). Due to rapid changes in taxonomy, we did not analyze the subspecies in our study, except for Vipera (Acridophaga) ursinii, to which detailed studies have confirmed the relevance of taxonomic unit (Ferchaud et al. 2012). While the taxonomic status of Anguis fragilis is still under debate, we considered the species complex as a single species (Gvoždík et al. 2010). We encountered a similar problem for Vipera (berus) nikolskii; therefore we did not examine it separately from Vipera berus (Zinenko et al. 2010).

We aggregated the occurrence records to the Universal Transverse Mercator (UTM) grid system at a spatial resolution of 25 km2 (UTM 5 × 5 km). The records with a spatial resolution of ≤ 25 km2 were assigned the corresponding UTM 5 × 5 km grid cell code using primarily the UTM index of localities (Lehrer and Lehrer 1990). The species occurrences with a spatial resolution of > 25 km2 were assigned only one grid cell code based on expert knowledge of the species’ habitat requirements (Scott et al. 2002, Franklin 2009) and visual help from the available satellite imagery and an overlaid KMZ file with the UTM 5 × 5 km grid in Google Earth v. 7.0.2 (Google Inc., CA). In order to georeference all records in the geodatabase in ArcGIS Desktop, we created a relationship between the table with species occurrence records and the UTM 5 × 5 km polygon feature class based on the grid cell code as a common attribute.

The occurrences were characterized based on the year of observation into old (if recorded before 1990) and recent (if after 1990, i.e. post Cold-War period) records. The distinction is based on the type of studies published and the style of work. Most papers before 1990 did not indicate the year of observation/collection and some did not even differentiate between new and old records. After 1990 the trend was to publish new locations in pure faunistic studies done within administrative units. Moreover, after 1990 access was not restricted anymore in areas close to the border and it involved also major changes in landscape use and an overall reduction of threats, mostly caused by a reduction in the use of pesticides and fertilizers in agriculture, land abandonment, changes in local hydrology caused by the decline in irrigation intensity and the decline of mining and industrial activities. If the year of observation was not mentioned in the publications, mostly in our national and local journals, then we subtracted 3 to 5 years from the date of publication, based on our estimate of the delay between actual fieldwork and time of publication.

Spatial patterns in species occurrences

Based on the number of reptile records per UTM 5 × 5 km grid cell, we tested for spatial autocorrelation using Global Moran’s I statistic under the null hypothesis that the occurrence records were evenly distributed. If the null hypothesis is rejected, the occurrence records are more spatially clustered (Z > 0) or dispersed (Z < 0) than expected (Fortin and Dale 2005). We assessed the local patterns of sampling bias using Getis Ord Gi* spatial statistic (Ord and Getis 1995). The metric allowed us to identify the clusters of UTM 5 × 5 cells where the sampling effort was significantly lower (GiZScore < -1.96; i.e., coldspots of occurrences) or higher (GiZScore > 1.96; i.e., hotspots of occurrences) than expected by chance. We set up the distance threshold to 7100 m in order to include the surrounding eight UTM grid cells (Getis and Ord 1992).

In order to assess the altitudinal distribution of each species, we extracted the mean altitude per grid cell from the SRTM digital topographic database (Jarvis et al. 2008) using ArcGIS Desktop Zonal Statistics geoprocessing tool. The grid cells intersecting the Romanian border were excluded from the analysis.

We calculated the species richness at a spatial resolution of 50 × 50 km and 10 × 10 km without considering subspecies separately. Mapping the species richness at a coarser resolution reduced the potential bias in sampling effort and allowed us a better understanding and visualization of regional patterns (Graham and Hijmans 2006).

We then calculated the Extent of Occurrence (EOO) as the minimum convex polygon, and the Area of Occupancy (AOO) as the total area of UTM 5 × 5 km grid cells where a species was reported (IUCN 2001). The exact area of UTM grid cells situated along the border and only partially located in Romania was taken into account in order to avoid overestimation.

Finally, we calculated a rarity index at a 50 × 50 km resolution (RI: Romano et al. 2012) as a measure of species range size in relation to the country area. The rarity index took values between 0 for widespread species and close to 100 for species with a restricted distribution.

All spatial statistics analyses were performed using geoprocessing tools under ArcGIS for Desktop 10.1 (ESRI, CA). Non-spatial statistical analyses were performed in Minitab 16 (Minitab Inc., PA).

Results
Species occurrences and ranges

We extracted and georeferenced 18036 reptile records from published papers (66%), museum collections (2%), personal communication from specialists, and our own field surveys (32%). The published papers with the reptile distribution records used in this study are listed in Appendix I, while the dynamics in time of their publication is presented in Figure 1. Among the records, 87.1% were dated after 1990, and 12.9% before 1990 (Appendix II and Table 1). Compared to the national density of reptile records per 100 km2 of 0.029 in the GBIF dataset (http://data.gbif.org, accessed 30.03.2013), our database increased it to 7.5 reptile records per 100 km2.

The species’ Extent of Occurrence ranged from 6439 km2 (Testudo hermanni) to 233497 km2 (Natrix natrix), while species’ Area of Occupancy ranged from 125 km2 (Eryx jaculus) to 37346 km2 (Lacerta agilis, Table 1).

Figure 1.

Number of publications with distribution data of reptiles in Romania (1823–2013). The reference line indicates the last country wide assessment published in 1961 (Fuhn and Vancea 1961).

Table 1.

The occurrence records for reptile species in Romania. EOO was estimated as 100% minimum convex polygon, and AOO was estimated as the total area of all UTM grid cell containing species records. Since not all UTM cells matched the 25 km2 area, the computed AOO is not a multiple of this value. Rarity index was calculated relative to a 50 × 50 km cell in order to check for regional patterns.

Species Total number of records Old records (before 1990) New records (after 1990) Number of UTM5 squares EOO (km2) AOO (km2) Rarity index
Emys orbicularis 753 131 622 561 232748 13491 22.7
Testudo graeca 1159 117 1042 154 14441 3748 91.1
Testudo hermanni 891 83 808 96 6439 2120 93.4
Anguis fragilis 1006 111 895 728 211520 17994 38.2
Eremias arguta 85 23 62 30 10851 625 95.1
Lacerta agilis 2554 268 2286 1521 231768 37346 13.8
Darevskia praticola 93 34 59 60 80161 1449 84.5
Lacerta trilineata 133 26 107 64 14028 1589 91.1
Lacerta viridis 2737 214 2523 1101 228111 26886 19.5
Zootoca vivipara 805 105 700 463 114513 11488 60.9
Podarcis muralis 672 111 561 374 167436 9093 56.9
Podarcis tauricus 1399 188 1211 358 163668 8671 70.7
Ablepharus kitaibelii 250 56 194 100 63728 2413 80.4
Eryx jaculus 5 4 1 5 10692 125 96.7
Coronella austriaca 527 98 429 389 213819 9581 39.8
Zamenis longissimus 500 73 427 363 201245 8868 50.4
Elaphe sauromates 47 15 32 27 17099 667 91.8
Dolichophis caspius 351 50 301 202 83585 4810 79.6
Natrix natrix 2180 163 2017 1361 233497 33202 15.4
Natrix tessellata 843 139 704 406 193909 9758 47.9
Vipera ammodytes 305 95 210 142 94177 3423 82.9
Vipera berus 663 138 525 464 144079 11515 47.1
Vipera ursinii 78 35 43 33 67560 729 91.1
Total 18036 2277 15759 - - -
Spatial patterns in species occurrences

Only 27.7% of the total UTM 5 × 5 km grid cells in Romania contained the sighting of at least one reptile species (Figure 2). The number of cumulative reptile records per cell had a clustered pattern, supporting the hypothesis of an overall bias in sampling effort (Z = 19.98, p < 0.001). Getis Ord Gi* spatial statistic identified several isolated sites, mainly protected areas, as hotspots of sampling effort (Figure 3).

Figure 2.

The presence records of reptiles per UTM 5 × 5 km grid cell in Romania. Records reported before 1990 were plotted as old records whereas those reported after 1990 were considered new records.

Figure 3.

The hotspots of sampling efforts within Romania. The p value was < 0.05 when Z scores took values between 1.96 and 61.60, suggesting a highly clustered pattern in the number of reptile occurrences per UTM 5 × 5 km grid cell. The numbered hotspots were: 1 Măcin Mountains 2 Jiului Gorges 3 Iron Gates and Mehedinţi Plateau 4 Sweet chestnut Arboretum of Baia Mare and the surroundings of Baia Mare town 5 Goşman Mountains and the surroundings of Piatra Neamţ town 6 the Danube Delta and Babadag Forest 7 the Danube Delta and Histria Archaeological Complex 8 Canaraua Fetii – Iortmac 9 Hagieni – Cotu Văii Forest 10 Nera river mouth and Baziaş.

The species richness map highlighted a lower sampling effort in the vast agricultural areas in the southern part of Romania (Figure 4 and Appendix III). By contrast, the southeastern and southwestern part of Romania, with Mediterranean influences (Rey et al. 2007), presented the highest diversity of reptile species with a maximum richness of 17 species per 50 × 50 km grid cell.

Figure 4.

The reptile species richness at a 50 × 50 km grid resolution within Romania.

The altitudinal range of reptiles varied between 0 and 2075 m (Figure 5), with only Zootoca vivipara and Vipera berus occurring above 2000 m, while Eremias arguta (occurring at maximum 34 m) and Eryx jaculus (occurring at maximum 88 m) are clearly restricted to lowlands. The updated distribution maps of reptile species in Romania are presented in Figures 628.

Figure 5.

The box and whisker plot of altitudinal distribution of reptile species in Romania. The boxes represent 25th–75th percentiles, upper and lower whiskers extends minimum and maximum data point within 1.5 box heights from the bottom and from the top of the boxes. Asterisks indicate outliers, lines and dots inside the boxes denote medians and means, respectively.

Figure 6.

Emys orbicularis.

Figure 7.

Testudo graeca.

Figure 8.

Testudo hermanni.

Figure 9.

Anguis fragilis.

Figure 10.

Eremias arguta.

Figure 11.

Lacerta agilis.

Figure 12.

Darevskia praticola.

Figure 13.

Lacerta trilineata.

Figure 14.

Lacerta viridis.

Figure 15.

Zootoca vivipara.

Figure 16.

Podarcis muralis.

Figure 17.

Podarcis tauricus.

Figure 18.

Ablepharus kitaibelii.

Figure 19.

Eryx jaculus.

Figure 20.

Coronella austriaca.

Figure 21.

Zamenis longissimus.

Figure 22.

Elaphe sauromates.

Figure 23.

Dolichophis caspius.

Figure 24.

Natrix natrix.

Figure 25.

Natrix tessellata.

Figure 26.

Vipera ammodytes.

Figure 27.

Vipera berus.

Figure 28.

Vipera ursinii moldavica and Vipera ursinii rakosiensis.

Discussion

We here present the first comprehensive distribution database of the reptile species occurring in Romania. We fulfilled the major requirements of data quality by (1) filtering the large amount of data for doubtful and erroneous records, (2) aggregating the known localities to a fine resolution of 25 km2, and (3) assessing the bias in sampling effort and thus providing useful information for further analyses.

We addressed a series of issues related to the quality and relevance of occurrence records. The most conspicuous issue was the overall biased density of records, which might appear due to differences in species detectability. For instance, our database showed that the snake species with low detectability (Hartel et al. 2009, Durso et al. 2011) had a low number of occurrences despite their wide range. Notable exceptions are the grass snake and the dice snake (genus Natrix), both of them being active during the day and easy to identify. The higher detectability of most lizards and both tortoise species resulted in a large number of occurrences and thus, more detailed distribution maps.

Another issue of concern was misidentification (e.g., Török 2012), due to difficulties in distinguishing the actual species based only on sightings. For example, it is rather difficult to distinguish among green lizards, especially between Lacerta viridis and Lacerta trilineata and we assumed that some of the occurrence records for those species were misallocated.

Although individuals of Testudo hermanni were reported from south-east Romania (Iftime 2002, Sós et al. 2008), we chose not to map the respective occurrences, considering them a case of vagrancy (Rozylowicz and Dobre 2010). The presence report of Testudo hermanni from Hunedoara county (Togănel 1993) was also not considered since it was the case of a recent human introduction into the wild (Rozylowicz 2008).

The analysis of sampling effort revealed a significantly higher number of sightings per grid cell within certain areas, mostly protected areas, such as the Jiului Gorges National Park, the Iron Gates Natural Park or the Măcin Mountains National Park. That was a common pattern previously reported from other countries and for different taxa (e.g., Loureiro et al. 2010, Botts et al. 2011). Despite those hotspots, the sampling effort was balanced across the country; therefore the dataset might be useful for additional analyses with only a simple trimming procedure required (Peterson et al. 2011). In the case of amphibians (Cogălniceanu et al. 2013), the sampling effort was biased towards the same hotspots but their number was significantly less numerous for reptiles. Nevertheless, the patterns in species richness support a similar statement revealing a higher richness in the warmer and drier climate in the south-west and the south-east parts of Romania (e.g., Rodríguez et al. 2005) and a constant richness in the rest of the country. Two gaps that require further investigations were revealed in the southern part of Romania (Oltenia and Bărăgan plains), probably being determined by the lack of research interest in those regions, due to their dominant agricultural landscape (Rey et al. 2007, Iojă et al. 2011). The gaps in reptile richness were similar to those of amphibian richness (Cogălniceanu et al. 2013).

The number of reptile species observed in Romania might increase in the near future with two species, namely Pseudopus apodus (Lepşi 1926, Cogălniceanu et al. 2008) and Mediodactylus kotschyi, bothpresent in Bulgaria and near the Romanian border (Gasc et al. 1997).

Several species were widespread across the country (e.g., Emys orbicularis, Lacerta viridis, Lacerta agilis, Natrix natrix, Natrix tessellata, Zamenis longissimus, Vipera berus), while others occurred only in the south of Romania (e.g., Testudo hermanni, Testudo graeca, Eremias arguta, Lacerta trilineata, Eryx jaculus, Elaphe sauromates, Dolichophis caspius). Lacerta viridis had the largest number of sightings mainly because it is a highly detectable species. Its AOO was surpassed only by Lacerta agilis and Natrix natrix, two other highly detectable species. Eryx jaculus is the rarest reptile in the country (see Table 1) with a single new report after 1990, based on a road-kill specimen (Covaciu-Marcov et al. 2012).

The different biogeographic regions overlapping in Romania will face different types and levels of these changes (Popescu et al. in press). Climate change may cause a shift in reptile ranges, although the Carpathian Mountains and the Steppic and Black Sea province of Dobrogea will act as a refuge, being considered critical areas for conservation (Araújo et al. 2006, Popescu et al. in press). To alleviate these threats by conservation activities (e.g., by establishing new protected areas), a key step is predicting the range responses to different climate change scenarios (Franklin 2009, Araújo et al. 2011). Our dataset provides a robust starting point for such analyses. Climate change alone does not result in a rapid range shift; instead, additionally, habitat loss due to the abandonment of traditional agricultural activities such as manual mowing, low-intensity grazing, small-size orchards and vineyards (Plieninger et al. 2006, Angelstam et al. 2013) might lead to a rapid contraction of the reptiles ranges (e.g., Testudo hermanni, Rozylowicz and Popescu 2013). Thus, our database also allows for further investigations resulting in concrete conservation activities.

Acknowledgements

We are grateful to the following persons for sharing their distribution data with us: Dr Arntzen Jan, Dr Bănărescu Petru (†), Dr Bereş Iosif, Buhaciuc Elena, Dr Gâldean Nicolae, Dr Hartel Tibor, Dr Kyek Martin, Dr Pârvulescu Lucian, Dr Oţel Vasile, Dr Skolka Marius, Sós Tibor, Talbot Neil, and Tallowin Oliver. We thank three anonymous reviewers for their constructive comments on this paper. This work was supported by two grants of the Romanian National Authority for Scientific Research, CNCS-UEFISCDI, project number PN-II-RU-TE-2011-3-0183 (principal investigator Laurenţiu Rozylowicz), and CNCS-UEFISCDI, project number PN-II-ID-PCE-2011-3-0173 (principal investigator Dan Cogălniceanu).

References
Angelstam P, Elbakidze M, Axelsson R, Čupa P, Halada L, Molnar Z, Pătru-Stupariu I, Perzanowski K, Rozylowicz L, Standovar T (2013) Maintaining cultural and natural biodiversity in the Carpathian Mountain ecoregion: need for an integrated landscape approach. In: Kozak J, Ostapowicz K, Bytnerowicz A, Wyżga B (Eds) The Carpathians: Integrating Nature and Society Towards Sustainability. Springer, Berlin Heidelberg, 393-424. doi: 10.1007/978-3-642-12725-0_28
Araújo MB, Guisan A (2006) Five (or so) challenges for species distribution modelling. Journal of Biogeography 33: 1677-1688. doi: 10.1111/j.1365-2699.2006.01584.x
Araújo MB, Thuiller W, Pearson RG (2006) Climate warming and the decline of amphibians and reptiles in Europe. Journal of Biogeography 33: 1712-1728. doi: 10.1111/j.1365-2699.2006.01482.x
Araújo MB, del Dedo-Garcimartín M, Pozo I (2011) Biodiversidade e Alterações Climáticas / Biodiversidad y Alteraciones Climát. Ministério do Ambiente e Ordenamento do Território & Ministerio de Medio Ambiente y Medio Rural y Marino, Lisboa, Madrid, 329 pp.
Böhm M, Collen B, Baillie JEM, Bowles P, Chanson J, Cox N, Hammerson G, Hoffmann M, Livingstone SR, Ram M, Rhodin AGJ, Stuart SN, van Dijk PP, Young BE, Afuang LE, Aghasyan A, García A, Aguilar C, Ajtic R, Akarsu F, Alencar LRV, Allison A, Ananjeva N, Anderson S, Andrén C, Ariano-Sánchez D, Arredondo JC, Auliya M, Austin CC, Avci A, Baker PJ, Barreto-Lima AF, Barrio-Amorós CL, Basu D, Bates MF, Batistella A, Bauer A, Bennett D, Böhme W, Broadley D, Brown R, Burgess J, Captain A, Carreira S, Castañeda MDR, Castro F, Catenazzi A, Cedeño-Vázquez JR, Chapple DG, Cheylan M, Cisneros-Heredia DF, Cogalniceanu D, Cogger H, Corti C, Costa GC, Couper PJ, Courtney T, Crnobrnja-Isailovic J, Crochet PA, Crother B, Cruz F, Daltry JC, Daniels RJR, Das I, de Silva A, Diesmos AC, Dirksen L, Doan TM, Dodd CK, Doody JS, Dorcas ME, Duarte de Barros Filho J, Egan VT, El Mouden EH, Embert D, Espinoza RE, Fallabrino A, Feng X, Feng ZJ, Fitzgerald L, Flores-Villela O, França FGR, Frost D, Gadsden H, Gamble T, Ganesh SR, Garcia MA, García-Pérez JE, Gatus J, Gaulke M, Geniez P, Georges A, Gerlach J, Goldberg S, Gonzalez JCT, Gower DJ, Grant T, Greenbaum E, Grieco C, Guo P, Hamilton AM, Hare K, Hedges SB, Heideman N, Hilton-Taylor C, Hitchmough R, Hollingsworth B, Hutchinson M, Ineich I, Iverson J, Jaksic FM, Jenkins R, Joger U, Jose R, Kaska Y, Kaya U, Keogh JS, Köhler G, Kuchling G, Kumlutaş Y, Kwet A, La Marca E, Lamar W, Lane A, Lardner B, Latta C, Latta G, Lau M, Lavin P, Lawson D, LeBreton M, Lehr E, Limpus D, Lipczynski N, Lobo AS, López-Luna MA, Luiselli L, Lukoschek V, Lundberg M, Lymberakis P, Macey R, Magnusson WE, Mahler DL, Malhotra A, Mariaux J, Maritz B, Marques OAV, Márquez R, Martins M, Masterson G, Mateo JA, Mathew R, Mathews N, Mayer G, McCranie JR, Measey GJ, Mendoza-Quijano F, Menegon M, Métrailler S, Milton DA, Montgomery C, Morato SAA, Mott T, Muñoz-Alonso A, Murphy J, Nguyen TQ, Nilson G, Nogueira C, Núñez H, Orlov N, Ota H, Ottenwalder J, Papenfuss T, Pasachnik S, Passos P, Pauwels OSG, Pérez-Buitrago N, Pérez-Mellado V, Pianka ER, Pleguezuelos J, Pollock C, Ponce-Campos P, Powell R, Pupin F, Quintero Díaz GE, Radder R, Ramer J, Rasmussen AR, Raxworthy C, Reynolds R, Richman N, Rico EL, Riservato E, Rivas G, da Rocha PLB, Rödel MO, Rodríguez Schettino L, Roosenburg WM, Ross JP, Sadek R, Sanders K, Santos-Barrera G, Schleich HH, Schmidt BR, Schmitz A, Sharifi M, Shea G, Shi HT, Shine R, Sindaco R, Slimani T, Somaweera R, Spawls S, Stafford P, Stuebing R, Sweet S, Sy E, Temple HJ, Tognelli MF, Tolley K, Tolson PJ, Tuniyev B, Tuniyev S, üzüm N, van Buurt G, Van Sluys M, Velasco A, Vences M, Veselý M, Vinke S, Vinke T, Vogel G, Vogrin M, Vogt RC, Wearn OR, Werner YL, Whiting MJ, Wiewandt T, Wilkinson J, Wilson B, Wren S, Zamin T, Zhou K, Zug G (2013) The conservation status of the world’s reptiles. Biological Conservation 157: 372-385. doi: 10.1016/j.biocon.2012.07.015
Botts EA, Erasmus BFN, Alexander GJ (2011) Geographic sampling bias in the South African Frog Atlas Project: implications for conservation planning. Biodiversity and Conservation 20: 119-139. doi: 10.1007/s10531-010-9950-6
Cogălniceanu D, Samoilă C, Tudor M, Skolka M (2008) Herpetofauna zonei costiere Cap Midia-Cap Kaliakra. In: Făgăraş M (Ed) Biodiversitatea zonei costiere a Dobrogei dintre Capul Midia şi Capul Kaliakra. Editura Ex Ponto, Constanţa, 323-370.
Cogălniceanu D, Cogălniceanu GC (2010) An enlarged European Union challenges priority settings in conservation. Biodiversity and Conservation 19: 1471-1483. doi: 10.1007/s10531-010-9777-1
Cogălniceanu D, Székely P, Samoilă C, Iosif R, Tudor M, Plăiaşu R, Stănescu F, Rozylowicz L (2013) Diversity and distribution of amphibians in Romania. ZooKeys 296: 35-57. doi: 10.3897/zookeys.296.4872
Covaciu-Marcov SD, Ghira I, Cicort-Lucaciu AS, Sas I, Strugariu A, Bogdan HV (2006) Contributions to knowledge regarding the geographical distribution of the herpetofauna of Dobrudja, Romania. North-Western Journal of Zoology 2: 88-25.
Covaciu-Marcov SD, Ferenţi S, Cicort-Lucaciu AS, Sas I (2012) Eryx jaculus (Reptilia, Boidae) north of Danube: a road-killed specimen from Romania. Acta Herpetologica 7: 41-47.
Cox NA, Temple HJ (2009) European red list of reptiles. Office for Official Publications of the European Communities, Luxembourg, 32 pp.
Durso AM, Willson JD, Winne CT (2011) Needles in haystacks: Estimating detection probability and occupancy of rare and cryptic snakes. Biological Conservation 144: 1508-1515. doi: 10.1016/j.biocon.2011.01.020
Evans D (2012) Building the European Union’s Natura 2000 network. Nature Conservation 1: 11-26. doi: 10.3897/natureconservation.1.1808
Ferchaud A-L, Ursenbacher S, Cheylan M, Luiselli L, Jelić D, Halpern B, Major Á, Kotenko T, Keyan N, Behrooz R (2012) Phylogeography of the Vipera ursinii complex (Viperidae): mitochondrial markers reveal an east–west disjunction in the Palaearctic region. Journal of Biogeography 39: 1836-1847. doi: 10.1111/j.1365-2699.2012.02753.x
Fortin MJ, Dale M (2005) Spatial Analysis. A Guide for Ecologists. Cambridge University Press, Cambridge, 365 pp.
Franklin J (2009) Mapping species distributions: spatial inference and prediction. University Press Cambridge, Cambridge, 336 pp.
Fuhn IE, Vancea Ş (1961) Reptilia (Ţestoase, Şopîrle, Şerpi). Fauna RPR, 14(2). Editura Academiei RPR, Bucureşti, 352 pp.
Gasc J-P, Cabela A, Crnobrnja-Isailovic J, Dolmen D, Grossenbacher K, Haffner P, Lescure J, Martens H, Martinez Rica JP, Maurin H, Oliveira ME, Sofianidou TS, Veith M, Zuiderwijk A (Eds) (1997) Atlas of amphibians and reptiles in Europe. Paris, 496 pp.
Getis A, Ord JK (1992) The analysis of spatial association by use of distance statistics. Geographical analysis 24: 189-206. doi: 10.1111/j.1538-4632.1992.tb00261.x
Ghira I, Venczel M, Covaciu-Marcov SD, Mara G, Ghile P, Hartel T, Török Z, Farkas L, Rácz T, Farkas Z, Brad T (2002) Mapping of Transylvanian herpetofauna. Nymphaea Folia naturae Bihariae 29: 145-201.
Gibbons J, Scott DE, Ryan TJ, Buhlmann KA, Tuberville TD, Metts BS, Greene JL, Mills T, Leiden Y, Poppy S, Wiine C (2000) The Global Decline of Reptiles, Déjà Vu Amphibians. BioScience 50: 653-666. doi: 10.1641/0006-3568(2000)050[0653:TGDORD]2.0.CO;2
Graham CH, Hijmans RJ (2006) A comparison of methods for mapping species ranges and species richness. Global Ecology and Biogeography 15: 578-587. doi: 10.1111/j.1466-8238.2006.00257.x
Gvoždík V, Jandzik D, Lymberakis P, Jablonski D, Moravec J (2010) Slow worm, Anguis fragilis (Reptilia: Anguidae) as a species complex: Genetic structure reveals deep divergences. Molecular Phylogenetics and Evolution 55: 460-472. doi: 10.1016/j.ympev.2010.01.007
Hartel T, Öllerer K, Farczády L, Monga C, Băncila R (2009) Using species detectability to infer distribution, habitat use and absence of a cryptic species: the smooth snake (Coronella austriaca) in Saxon Transylvania. Acta Scientiarum Transylvanica 17: 61-76.
Hof C, Araújo MB, Jetz W, Rahbek C (2011) Additive threats from pathogens, climate and land-use change for global amphibian diversity. Nature 480: 516-519. doi: 10.1038/nature10650
Iftime A (2002) Testudo hermanni Gmelin, 1789 in Dobroudja (SE Romania), with comments on conservation. Herpetozoa 15: 183-186.
Iojă CI, Pătroescu M, Rozylowicz L, Popescu VD, Vergheleţ M, Zotta MI, Felciuc M (2010) The efficacy of Romania’s protected areas network in conserving biodiversity. Biological Conservation 143: 2468-2476. doi: 10.1016/j.biocon.2010.06.013
Iojă CI, Rozylowicz L, Pătroescu M, Niţă MR, Onose DA (2011) Agriculture and Conservation in the Natura 2000 Network: A Sustainable Development. In: Andreopoulou Z, Manos BM, Polman N, Viaggi D (Eds) Agricultural and Environmental Informatics, Governance and Management: Emerging Research Applications. IGI Global, Hershey, Pennsylvania, 339–358. doi: 10.4018/978-1-60960-621-3.ch018
IUCN (2001) IUCN Red list categories and criteria: version 3.1. http://www.iucnredlist.org/technical-documents/categories-and-criteria/2001-categories-criteria [accessed 5.02.2013].
Jarvis A, Reuter H, Nelson A, Guevara E (2008) Hole-filled SRTM for the globe Version 4. CGIAR-SXI SRTM 90m database: http://srtm.csi.cgiar.org [accessed 5.02.2013]
Krecsák L, Iftime A (2006) A review of the records of the Sand boa (Eryx jaculus) in Romania. Herpetological Bulletin 98: 31–34.
Krecsák L, Zamfirescu Ş (2008) Vipera (Acridophaga) ursinii in Romania: historical and present distribution. North-Western Journal of Zoology 4: 339-359.
Lehrer AZ, Lehrer MM (1990) Cartografierea faunei şi florei României (Coordonate arealografice). Ed. Ceres, Bucureşti, 295 pp.
Lepşi I (1926) Ophisaurus apus Pall. neu für Rumänien. Verhandlungen und Mitteilungen des Siebenbürgischen Vereines für Naturwissenschaften zu Hermannstadt 75/76: 45–48
Loureiro A, Ferrand N, Carretero MA, Paulo O (2010) Atlas dos Anfíbios e Répteis de Portugal. Sfera do Caos, Lisboa, 256 pp.
Ord JK, Getis A (1995) Local spatial autocorrelation statistics: distributional issues and an application. Geographical analysis 27: 286–306. doi: 10.1111/j.1538-4632.1995.tb00912.x
Peterson AT, Soberon J, Pearson RG, Martinez-Meyer E, Nakamura M, Araújo MB (2011) Ecological niches and geographic distributions (MPB-49). Princeton University Press, 328 pp.
Phillips SJ, Anderson RP, Schapire RE (2006) Maximum entropy modeling of species geographic distributions. Ecological modelling 190: 231-259. doi: 10.1016/j.ecolmodel.2005.03.026
Pineda E, Lobo JM (2009) Assessing the accuracy of species distribution models to predict amphibian species richness patterns. Journal of Animal Ecology 78: 182-190. doi: 10.1111/j.1365-2656.2008.01471.x
Plieninger T, Höchtl F, Spek T (2006) Traditional land-use and nature conservation in European rural landscapes. Environmental science & policy 9: 317-321. doi: 10.1016/j.envsci.2006.03.001
Popescu DV, Rozylowicz L, Cogălniceanu D, Niculae IM, Cucu AL (in press) Moving into protected areas? Setting conservation priorities for Romanian reptiles and amphibians at risk from climate change. PLoS ONE.
Proess R (Ed) (2003) Verbreitungsatlas der Amphibien des Großherzogtums Luxemburg. Musée national d’histoire naturelle, Luxembourg, 92 pp.
Rey V, Groza O, Ianoş I, Pătroescu M (2007) Atlas de la Roumanie. Reclus, Montpelier and Paris, France, 211 pp.
Rodríguez MÁ, Belmontes JA, Hawkins BA (2005) Energy, water and large-scale patterns of reptile and amphibian species richness in Europe. Acta Oecologica 28: 65–70. doi: 10.1016/j.actao.2005.02.006
Romano A, Bartolomei R, Conte AL, Fulco E (2012) Amphibians in Southern Apennine: distribution, ecology and conservation notes in the “Appennino Lucano, Val d’Agri e Lagonegrese” National Park (Southern Italy). Acta Herpetologica 7: 203-219.
Rozylowicz L (2008) Metode de analiză a distribuţiei areal-geografice a ţestoasei lui Hermann (Testudo hermanni Gmelin, 1789) în România. Studiu de caz: Parcul Natural Porţile de Fier. Editura Universităţii din Bucureşti, Bucureşti, 169 pp.
Rozylowicz L, Dobre M (2010) Assessing the threatened status of Testudo hermanni boettgeri Mojsisovics, 1889 (Reptilia: Testudines: Testudinidae) population from Romania. North-Western Journal of Zoology 6: 190-202.
Rozylowicz L, Popescu VD (2013) Habitat selection and movement ecology of eastern Hermann’s tortoises in a rural Romanian landscape. European Journal of Wildlife Research 59: 47-55. doi: 10.1007/s10344-012-0646-y
Scott JM, Heglund PJ, Morrison ML, Raven PH (2002) Predicting species occurrences: issues of scale and accuracy. Island Press, Washington D.C., 840 pp.
Sós T, Daróczi S, Zeitz R, Pârâu L (2008) Notes on morphological anomalies observed in specimens of Testudo hermanni boettgeri Gmelin, 1789 (Reptilia: Chelonia: Testudinidae) from Southern Dobrudja, Romania. North-Western Journal of Zoology 4: 154-160.
Sós T (2011) În obiectiv: Ţestoasa de apă europeană, Emys orbicularis. Asociaţia Ecouri Verzi, Cluj Napoca, 110 pp.
Sós T, Kecskés A, Hegyeli Z, Marosi B (2012) New data on the distribution of Darevskia pontica (Lantz and Cyrén, 1919) (Reptilia: Lacertidae) in Romania: filling a significant gap. Acta Herpetologica 7: 175-180.
Speybroeck J, Beukema W, Crochet PA (2010) A tentative species list of the European herpetofauna (Amphibia and Reptilia) - an update. Zootaxa 2492: 1-27.
Stuart SN, Hoffmann M, Chanson JS, Cox NA, Berridge RJ, Ramani P, Young BE (2008) Threatened Amphibians of the World. Lynx Edicions/IUCN/Conservation International, Barcelona/Gland/Arlington, 758 pp.
Strugariu A, Gherghel I, Ghira I, Covaciu-Marcov SD, Mebert K (2011) Distribution, Habitat Preferences and Conservation of the Dice Snake (Natrix tessellata) in Romania. Mertensiella 18: 272-287.
Togănel F (1993) Contribuţii la cunoaşterea răspândirii speciei Testudo hermanni Gmelin (Fam. Testudinae) în Transilvania. In: Sesiunea Ştiinţifică Secţiunea Biologie. Facultatea Biologie Geografie Geologie, Secţia Biologie Universitatea „Babeş-Bolyai”, Cluj-Napoca, 100–101 pp.
Török Z (2012) Doubtful records of reptile species in some areas of the Danube Delta Biosphere Reserve (Romania). Analele Ştiinţifice ale Institutului Delta Dunării 18: 223-232.
Zinenko O, Turcanu V, Strugariu A (2010) Distribution and morphological variation of Vipera berus nikolskii Vedmederja, Grubant et Rudaeva, 1986 in Western Ukraine, The Republic of Moldova and Romania. Amphibia-Reptilia 31: 51-67. doi: 10.1163/156853810790457885
Appendix I

The publications used to compile the distribution of reptile species native to Romania. (doi: 10.3897/zookeys.341.5502.app1) File format: Microsoft Word file (doc).

Appendix II

The reptiles occurrence records within Romania per time intervals (1961 represented the publishing year of the last country wide assessment). (doi: 10.3897/zookeys.341.5502.app2) File format:Microsoft Word file (doc).

Appendix III

The reptile species richness at a 10 × 10 km grid resolution within Romania (SCI = Natura 2000 Sites of Community Importance). (doi: 10.3897/zookeys.341.5502.app3) File format: Microsoft Word file (doc).