Acta Herpetologica 11(2): 197-212, 2016 ISSN 1827-9635 (print) © Firenze University Press ISSN 1827-9643 (online) www.fupress.com/ah DOI: 10.13128/Acta_Herpetol-18061 Olim palus, where once upon a time the marsh: distribution, demography, ecology and threats of amphibians in the Circeo National Park (Central Italy) Antonio Romano1,*, Riccardo Novaga2, Andrea Costa1 1 Consiglio Nazionale delle Ricerche, Istituto di Biologia Agroambientale e Forestale, Via Salaria Km 29,300 I-00015 Monterotondo Scalo, Rome, Italy. * Corresponding author. E-mail: antonioromano71@gmail.com 2 Viale dello Statuto 37, Latina, Italy. Submitted on 2016, 27th February; revised on 2016, 27th June; accepted on 2016, 7th July Editor: Gentile Francesco Ficetola Abstract. The Circeo National Park lies in a territory that was deeply shaped by human activity, and represents one of the few remaining patches of plain wetland habitat in Central Italy. In this study distribution and few demographic information of the amphibians in the Park were provided. Seven species and 25 bibliographic and 84 original breeding sites were recorded, and population size estimations were carried out for a population of these three species: Pelophylax sinkl esculentus, Bufo balearicus and Rana dalmatina. For the studied populations of pool frog and green toad the oper- ational sex ratio and the demographic effective population size was also estimated. For Rana dalmatina, which is strictly associated to forest environment, a positive and significant correlation between the number of egg clutches and maxi- mum depth of the swamps was found. The State plain forest is the most important habitat for amphibians’ conservation in the park. The occurrence of dangerous alien species was investigated and they are evaluated as the major threat for amphibians in the park, especially the crayfish Procambarus clarkii in the State plain forest. Index of Calling Survey were performed for anurans and the medians did not differ among species. The potential distribution of amphibians in the Park was evaluated by building a species distribution model. Finally, the absence of three species reported in litera- ture in the 60’s of the last century (Bombina pachypus, Salamandrina perspicillata, Rana italica) is also discussed. Keywords. Alien species, Capture-Marking-Recapture, effective population size, Index of Calling Survey, land recla- mation, Species Distribution Models, swamps. INTRODUCTION Until the early twenties of the past century, marsh- lands dominated the landscapes along the coast about 45 km southeast of Rome, from Anzio to Terracina towns between the  Tyrrhenian sea and the Volsci Chain, at inland distances from  15 to 25  km. The territory, which now belongs to the province of Latina (Latium region), was an extensive marsh at about sea level originated in an alluvial plain (e.g., Linoli, 2005).  Forested swamps domi- nated above sea level while areas below the sea level were covered by mud flats and pools. These so called “Paludi Pontine” (Pontine Marshes) were the subject of land rec- lamation works, performed periodically since the pre- Roman period, initially by the Italic tribe  of Latins, but with scarce success. Land reclamation was conducted extensively with considerable success by Fascist regime, starting in the 1920s and radically changing the land- scape of the area, which was converted into an extensive agricultural plain and new towns were founded and built (Littoria, renamed Latina, Pontinia, Sabaudia, Aprila and Pomezia are the most important ones). Out of the origi- 198 Antonio Romano et alii nal 20,700 hectares of forest and swampland, about 3,200 were put under protection with the creation in 1934 of the Circeo National Park (CNP), which includes also other not forested areas. Among terrestrial vertebrates, amphibians are the class most strictly associated to wet- lands. They are an important component of biodiversity and suffer a recent worldwide decline (Stuart et al., 2004; Wake and Vredenburg, 2008). In Italy protected areas play a key role in conservation of amphibians, and act as stepping stones in the face of climate change (D’Amen et al., 2011). The area where CNP falls is classified in the highest irreplaceability naturalistic values (Maiorano et al., 2006). However data on the fauna before the land rec- lamation was available only for mammals and birds (Lep- ri, 1935). Consequently, information on which herpeto- logical species occurred in this area before land reclama- tion could be only deducted by current species distribu- tion on wider areas (Bruno, 1973; Bruno 1981; Bologna et al., 2000), by specific works on the fauna of the CNP (Carpaneto, 1986; Ravenna, 2013; Cinquegranelli et al., 2015), and by herpetological census of surrounding areas which were also subjected to land reclamation in the 30’s (Novaga et al., 2013). In particular, the recent paper of Cinquegranelli et al. (2015) provided few updated data on distribution of amphibians in the Park. However Cinque- granelli an co-authors surveyed only 15 sites, and infor- mation on habitat use, occupancy level and species detec- tion probability were the main goals of their work. The aim of our study was fourfold. First, by perform- ing an extensive survey on aquatic habitats, we provide an updated species distribution of amphibians. Second, we estimated some demographic and abundance param- eters. Third, we evaluated habitat preferences of amphib- ians and we provided potential distribution information. Finally, we detected relevant threats for species, popula- tion or habitats and we provided information on conser- vation measures. MATERIALS AND METHODS Study area The Circeo National Park covers 8,484 hectares of a coastal area of the Central Italy, and it consists of five main environ- ments: the plain forest, four coastal lakes, the coastal dune area, the limestone massif of Mount Circeo (541 m asl, a promon- tory that marks the southwestern limit of the former  Pontine Marshes) and the island of Zannone. The plain State Forest covers about 3,190 hectares and mainly consists of deciduous woods. With many areas few meters below the sea level, most of the Park does not exceed the 30 meters a.s.l. The Park ranges from latitude 41°13’N to 41°24’N, and from longitude 12°50’E to 13°07’E. CNP is covered by wooded areas and semi-natural habitats for the 56%, agricultural fields cover the 18%, water bodies the 13%, artificial territories the 11% and, finally, a small portion is covered by wetlands (2%; Giagnacovo et al., 2003). The climate falls in the Lower Mesomediterranean Thermotype, Upper Subhumid Ombrotype (Blasi, 1994). Mean temperature of this area is from 9.5 °C to 17.1 °C, and temperatures below zero are uncommon. Precipitation is mainly concentrated in autumn and early winter  (October-December); relative humid- ity is high all over the year, wind is frequent with a South-West- ern dominance (Padula , 1985). Distribution and species occurrence The data reported in literature, when they were provid- ed with enough accuracy, were georeferenced. Field surveys (March-September 2015) were preceded by a careful analy- sis of the maps produced by the Istituto Geografico Militare (I.G.M, 1:25000) and by a analysis of satellite images to detect water bodies not reported in the maps (see Romano et al., 2012). Information on methods used to detect the occurrence of amphibians are reported in detail by Romano et al. (2010; 2012). To describe species rarity and their diffusion we used a method that weights both diffusion (W: wide; M: medium; L: limited) and density (C: common; F: frequent; R: rare) and it consists of a graph of the relationship between the coverage (%) of the UTM grid (we used a mesh of 2x2 km) and the mean number of observations for a square occupied by each species, according to the method proposed by Doria and Salvidio (1994) and used in other herpetological studies (Turrisi and Vaccaro, 2004; Romano et al., 2012). In the computation of the score to build the graph of rarity and diffusion we decided to use only a subset of meshes. We included only the meshes where at least one species record occurred (i.e., 35 meshes), excluding, for example, meshes in highly urbanised areas with no data. Ecology Breeding aquatic sites were assigned to the following typol- ogies: (i) ponds and marshes; (ii) slow running waters: ditches, streams and artificial channels (iii) rheocrenic springs (which were checked s up to 50 meters from they source); (iv) forest swamps (including those whose filling is partially due to lim- nocrenic springs); (v) artificial tanks; (vi) brackish coastal lakes. Correspondence Analysis (CoA) was used to identify associa- tions among amphibian species and aquatic habitats. Consider- ing that the variance of the data was homogenous (Levene’s test for homogeneity, based on means, P = 0.758), the hypothesis that habitat categories may host different syntopic number of species was tested using one-way ANOVA. Correlation between the habitat availability and the number of species for each habi- tat type was tested using non-parametric Spearman’s rank corre- lation. Analyses concerning habitat typologies were performed both on original and bibliographic data, when the latter could be certainly assigned to a given habitat typology. Ecological analyses were performed in the statistical package PAST (Ham- mer et al., 2001). 199Amphibians of the Circeo National Park Population abundance estimation Populations estimates were performed using four different methods: Index of calling survey (ICS), capture-mark-recap- ture  (CMR), removal sampling (RES), and egg mass counts (EMC). For the ICS, which is the relative measure of calling den- sity (Mossman et al., 1998; Weir and Mossman, 2005), surveys were performed as reported in Dorcas et al. (2009) and Romano et al. (2012). ICS may vary among 0 and 3. Considering that ICS provides scores as ordinal measures and with many ties, the scores were analysed using the median test as performed in Statistica® ver. 5.0 (Statistica package, Statsoft Inc., USA). Differ- ences in calling survey scores among anuran species were com- pared using the non parametric Kruskall-Wallis test and relative post hoc comparisons. Using photo-identification, CMR analysis was applied on a population of Bufo balearicus breeding in an artificial tank in the inner of Sabaudia town; we used the software CAPTURE (Otis et al. 1978) to estimate adult population size (N) analysing data from the four sampling sessions performed in a short time range (23 September – 2 October; light rain, between 9 and 11 p.m.). RES was applied on a population of Pelophylax sinkl. escu- lentus breeding in the artificial pond of the headquarter of the Park (Sabaudia). We used the jackknife estimator of Pollock and Otto (1983) as performed in the program CAPTURE to estimate N on the basis of three removal sessions performed in about three hours. Both for Bufo balearicus and Pelophylax sinkl. esculentus populations, the operational sex ratio (calculated just on the number of males and females captured and not on population estimates), in accordance with Wilson and Hardy (2002), was expressed as the proportion of mature males, i.e. males/(males + females). The demographic effective population size (Ne) was estimated as Ne = (4*Nm*Nf)/( Nm+Nf), where N is the number of mature males (m) or females (f ) individuals (Wright, 1938), which is a widely used equation to obtain demographic estimates of Ne (e.g. Jehle et al., 2001; Schmeller and Merila, 2007). EMC was used in March to estimate breeding populations of Rana dalmatina in seven natural ponds in forest environ- ment; egg counts were performed in a unique sampling ses- sion per pond; as R. dalmatina is an explosive breeder (Sofia- nidou and Kyriakopoulou-Sklavounou, 1983; Guarino and Bellini, 1993) and each female lays a single egg mass per sea- son (Nollert and Nollert, 1992), consequently we considered the counts of egg masses as a good proxy of the minimum female population size (Griffiths and Raper, 1994; Grossenbacher et al., 2002). Surface areas of breeding ponds and swamps were esti- mated by walking the  perimeter  of each site with a GPS, auto- matically calculating the area inside the resulting shape. Where egg masses of R. dalmatina were counted, Spearman  rank  cor- relation  was used to test the association between density of egg masses and size or maximum depth of the swamps. Correla- tion between water surface and their maximum depth was also tested for all 15 swamps where R. dalmatina breed (correlations were performed in PAST; Hammer et al., 2001). Potential distribution The potential distribution of amphibians in the CNP was evaluated by building a species distribution model (SDM), using the algorithm of maximum entropy (Maxent, Phillips et al., 2006). Selection of environmental data layers, to be employed for these analysis, was based on availability and  a priori  expec- tation of influences on amphibian population. Considering that the Park area runs along the coast and it is a small area, pre- cipitations and other climatic variables were considered homog- enous in that area and were not included among variables. We used the Digital Elevation Model (DEM) with a spatial resolu- tion of 90m to obtain other topographic variables. All vari- ables were resampled in order to match the 90m resolution of the DEM. The environmental variables used were: Corine land cover data 2006; Landsat tree cover, representing the percent- age of canopy cover of trees higher than 5 m; distance from for- est swamps; distance from running waters; distance from lakes; topographic variables calculated from DEM were elevation, aspect (Northness or Eastness), valley depth, Topographic Wet- ness Index (TWI; Sörensen et al., 2006), Topographic Rugged- ness Index (TRI; Riley et al., 1999), wind exposure, direct inso- lation (kw/h per square meter). To build the SDM, both original and bibliographic data were pooled together. The whole data set was split in two subgroups by random selection: 70% of points were used for building the model (training), while the remain- ing 30% of point-data were employed to evaluate its predictive power. This procedure was repeated 10 times for each species, generating an averaged prediction of amphibians’ distribution. Finally, the predictive power of the model was evaluated by calculating the area under the receiver operating characteris- tic curve (AUC). Analysis were conducted in software MaxEnt 3.3.3k and default software settings were used, with the excep- tion of the employment of bootstrapping procedure and num- ber of iterations (1000). All GIS processing to obtain the above mentioned layers was performed using software SAGA Gis. Threats Considering that in the Park there is a high-density road network and that road mortality may be considered as an addi- tional factor contributing to the amphibian decline (Puky, 2006; Glista et al., 2008), roads were checked systematically. In particu- lar we controlled the two parallel coastal road and waterfront roads (both about 25 km) after rains in spring and summer. The causes for amphibian declines are many (e.g., Col- lins, 2010) but habitat loss and alteration, predator alien species and emerging diseases are among the leading. Habitat loss and habitat alterations in progress was recorded and we searched for exotic animal species for which is well established that they are a threats for amphibians (i.e., predator fishes, crayfishes). Alien species were searched in every site where amphibians surveys 200 Antonio Romano et alii were performed, both using visual surveys (using a binocular too) and blind dip netting (20-30 dip netting per site). Informa- tion on fish was also obtained by fishermen. RESULTS Distribution and species occurrence Distribution was the result of both bibliographic (when they could be georeferenced with a good approxi- mation, i.e. 500 m; n = 25) and original data (n =84). All available literature (Bruno, 1973; Bruno 1981; Bologna et al., 2000; Carpaneto, 1986; Ravenna, 2013; Cinquegranelli et al., 2015) is consistent in reporting the following spe- cies in the park: Triturus carnifex (Laurenti, 1768), Lisso- triton vulgaris (Linnaeus, 1758), Bufo balearicus Boettger, 1880, Bufo bufo (Linnaeus, 1758), Hyla intermedia Bou- lenger, 1882, Pelophylax sinkl. esculentus (Linnaeus, 1758) and Rana dalmatina Fitzinger in Bonaparte, 1838. The pool frog synklepton is formed by two entities: the paren- tal species, P. lessonae (Camerano, 1882) and its hemi- clonal hybrid, the klepton P. kl. easculentus (Linnaeus, 1758). Our research confirmed the occurrence of all these species. Conversely, the occurrence of Salamandrina per- spicillata (Savi, 1821), Bombina pachypus (Bonaparte, 1838) and Rana italica Dubois, 1987 reported by Bruno (1981) in many localities of the park was not confirmed by our samplings, consistently with other researches (Carpaneto, 1986; Bologna et al., 2000). Species distribu- tions are reported in Figs. 1 and 2 The graph related to diffusion and density (Fig. 3) showed that Triturus carnifex is the species having the most critical situation in the Park, with very limited dif- fusion and occurrence (i.e., low number of sites). No spe- cies are positioned in the upper right corner of the graph, however Rana dalmatina and Pelophylax sinkl. esculentus are of little concern. All the other species, considering their distribution and rarity, do not appear to be particu- larly threatened. For the particular situation and position in the Fig. 3 of Bufo balearicus, at the boundary between 4 quadrants, see Discussion. Ecology Fig. 4 shows aquatic habitat preferences of each spe- cies. Both newts occurred in a limited number of habitat typologies while the more ecologically plastic species was the tree frog that breeds in all the habitat typologies. Var- iation in amphibian species composition among habitats is highly explained (more than 80%) by the first two axes of the CoA scatter plot (Fig. 5). Associations were found between the amphibian species and the various aquatic habitat typologies. The results show that all the species are quite different from each other for habitat preference. Few species are strictly associated to only one aquatic habitat typology: Rana dalamatina, Hyla intermedia and Pelophylax sinkl. esculentus are closely linked to forest swamps, ponds and slow running waters respectively. Bufo balearicus is associated with ponds and lakes while B. bufo seems to be associated with springs and artifi- cial tanks. The newts display a moderate and comparable association with forest swamps, ponds and slow running waters. The number of syntopic species did not differ sig- nificantly among the six habitat categories (one-way ANOVA, F1,5 = 0.32, P = 0.90). Number of species in the different aquatic habitat typologies ranges between 3 (lakes) and 7 (ponds and slow running waters). Spear- man’s  rank  correlation did not detect a significant asso- ciation between the habitat availability and the number of species for each habitat type (r = 0.806; P = 0.083). Population abundance estimation Index of calling survey (ICS) was performed to sev- eral populations of B. balearicus (16 populations), H. intermedia (16), P. sinkl. esculentus (52) and B. bufo (16). For the three first species all the four ICS scores were recorded (0-3), but for the latter the highest score lacked (Fig. 6). The score 1 for both toads clearly exceeded oth- er scores proportionally. Medians of ICS did not differ among species (χ2 = 5.678, df = 3, P = 0.128). Capture-Marking-Recapture (CMR) was applied on a population of Bufo balearicus living in a urban mead- ow (5800 m2) surrounding the breeding site (a concrete tank). During the four sampling sessions, we performed a total of 56 captures in which 26 different individuals were marked. Closure test confirmed that the population was closed (z = 1.89; p = 0.97). The recapture rate was high and the estimated population was 27 ± 1.32 adult toads (estimate ± SE; CI 95% = 27-33). Considering that the operational sex ratio was extremely balanced (0.46), the estimate of effective population size (Ne=25.8) was similar to that of the adult population census size (N). The N of P. sinkl. esculentus breeding in an artificial ponds (88 m2), was 60 ± 6 individuals (estimate ± SE; CI 95% = 53-78). The operational sex ratio was strongly male biased (0.94) and, as a result, effective population size was much lower than N (Ne=11.25). The seven swamps where egg mass counting (EMC) of R. dalmatina was performed had greatly variable water surface area (mean ± SD = 4637.43 ± 5691.36 m2; range = 124-15860 m2) and maximum depth (mean 201Amphibians of the Circeo National Park ± SD =64.28 ± 26.37 cm; range = 40-110 cm). Surface area and maximum depth were not significantly corre- lated (r=0.654; p=0.128). Correlation remained above the significance threshold even considering all the swamps (N = 15) where reproductive activity was recorded (i.e., pooling sites where egg counts was performed and these where it was not performed) (r = 0.506; P = 0.053). Number and density of egg masses greatly varied among sites (Fig. S7). A total of 1419 egg masses were recorded (mean ± SD = 202.71 ± 199.09 SD; range 13-604) and their density ranged from 0.01 to 0.47 egg clutches per square meter (mean ± SD = 0.14/m2 ± 0.18). A signifi- Lissotriton vulgaris Salamandrina perspicillata Bombina pachypus Rana italica Bufo balearicus a b c d Figure 1 o Triturus carnifex Lissotriton vulgaris Bufo balearicus a b c d Salamandrina perspicillata Bombina pachypus Rana italica o m Fig. 1. Distribution of amphibians in the Circeo National Park (Central Italy). a = distribution of three species reported just once in the literature but not recorded during further researches. For figs. b, c and d, circles = original data; triangles = bibliographic data. Grid reports 10x10 km UTM squares. Dashed line: Park boundary. Dotted lines: surface hydrography. Dotted areas: lakes. The urban area of Sabaudia is showed in grey. 202 Antonio Romano et alii cant correlation among eggs parameters and swamps features was detected only between the number of egg masses and maximum depth of the swamps (r = -0.509, P = 0.249; r = -0.055, P = 0.919, for egg density vs. swamp size or depth respectively; r = 0.643, P = 0.109; r = 0.836, P = 0.025 for egg number vs. swamp size or depth respectively. See Fig. 7). Potential distribution Species distribution models (SDM) were built for all seven amphibian species occurring in the Park. The number of available data to be employed (both as train- ing and test data) for model building was dependent on the species (see Tab. 1 for details) and ranged between Hyla intermediaBufo bufo Figure 2 a b c d Pelophylax sinkl esculentus Rana dalmatina Fig. 2. Distribution of amphibians in the Circeo National Park (Central Italy). Circles = original data; triangles = bibliographic data. Grid reports 10x10 km UTM squares. Dashed line: Park boundary. Dotted lines: surface hydrography. Dotted areas: lakes. The urban area of Sabaudia is showed in grey. 203Amphibians of the Circeo National Park 11 (Triturus carnifex) and 57 (Pelophylax sinkl. escu- lentus). All models shown a high predictive power, as revealed by AUC: indeed all averaged models received an AUC value > 0.970 (ranging from 0.974 for Bufo balearicus to 0.997 for Rana dalmatina) with the excep- tion of the models regarding Triturus carnifex, which shown a slightly lower predictive power (AUC = 0.947). As a rule, among all variables included in the analyses, the most important ones for the major part of the spe- cies are the distance from water bodies (both lakes and swamps), tree cover and insolation, altogether with many Corine categories. The detailed results, regarding the list of variable effect, together with variable percentage con- tribution, for each species, is reported in Tab. 1. Poten- tial distribution maps are presented as supplementary materials (Figs. S1-S4) Threats In the Circeo National Park we found three ponds and a swamps which are breeding sites of Lissotriton vul- garis, Pelophylax sinkl. esculentus and Rana dalmatina that suffer progressive filling with soil, reduction of the depth and earlier dry up. In April 2015, remains of sev- eral R. dalmatina eggs desiccated before hatching were recorded in these ponds. The water body with highest danger of disappearing was the one in locality “Cerreto Fig. 3. Relationship between percentage of amphibian species occurrence in UTM square grids (2x2 km) and mean number of observations per UTM square. Tricar = Triturus carnifex; Lisvul = Lissotriton vulgaris; Bufbuf = Bufo bufo; Bufbal = Bufo baleari- cus; Hylint = Hyla intermedia; Randal = Rana dalmatina; Pelesc = Pelophylax synkl. esculentus. Figure 4 Fig. 4. Habitat partitioning (number of sites on the left vertical axe) of amphibians in the Circeo National Park (Central Italy). Codes of species are as reported in Fig. 3. Fig. 5. Correspondence Analysis (CoA) scatter plot illustrating vari- ations of amphibians species (black dots) distribution with aquatic habitat typologies (white squares). The percentages of variation explained by each axis are given in round brackets (codes of species are as reported in Fig.3). 204 Antonio Romano et alii Fontana”, a pond with surface area of 240 m2 and maxi- mum depth of about 20 cm. No relevant evidence of road killing on amphibians were recorded and, during about 20 surveys, just few individuals of B. balearicus were found crushed in sum- mer, after a rainy day, in a small stretch of the seafront road (loc. Bufalara, Sabaudia). Reproductive populations of five alien species, iden- tified by the literature as threat to amphibians (see Dis- cussion), were detected. In ascending order of threat to amphibians they were one reptile, three fishes and a crus- tacean: Trachmeys scripta; Carassius auratus, Lepomis gib- bosus, Gambusia sp., Procambarus clarkii. On the whole they occurred at least in 14 UTM mashes (2x2 km), that is 32% of the total meshes occupied by the Park for at least 10% of their surface (n= 44; Fig. 8). Amphibian occurred in all 2x2 km UTM meshes where alien species were recorded. Original data concerning the distributions of these alien species are shown in supplementary materi- als (Figs S5-S6). DISCUSSION Species occurrence and their actual and potential distribu- tion In the 84 sites with amphibians, we found seven species in the Circeo National Park as reported in lit- erature (Carpaneto, 1986; Cinquegranelli et al., 2015). However the occurrence ratios among species we found (L. vulgaris: 17%; T. carnifex 4%; B. balearicus: 18%; B. bufo: 25%, H. intermedia: 19%; P. synkl. esculentus: 62%; Fig. 6. Percentage of Index of Calling Survey (ICS) scores for four anuran species in the Circeo National Park (code of species as in Fig. 3). N = number of sites for each species. Fig. 7. Trend of the relationship between the number of Rana dal- matina egg masses and maximum depth of seven swamps. Alpha- numerical codes refer to swamps: VER (Piscina della Verdesca), BAG (Piscina delle Bagnature), CAR (Piscina del Carpino). Corre- lation was statistically significant (Spearman correlation, r = 0.836, P = 0.025). 0 5 10 15 20 25 Trachemys scripta Carassius auratus Lepomis gibbosus Gambusia affinis Poracambarus clarkii Alien species sites UTM (2x2 km) % UTM with amphibians (n=35) Fig. 8. Number of sites and UTM mashes (2x2 km) where alien species were recorded in the Circeo National Park. 205Amphibians of the Circeo National Park R. dalmatina: 30%) differed, in some cases, from those reported by Cinquegranelli et al., 2015 (20%, 13%, 20%, 20%, 53%, 67%, 20%, respectively). These authors, per- forming 14 visit in each site, provide an interesting con- tribution testing the species detection probabilities (p), misdetection rates (Mr) and minimum number of visits (Nm maximum=10.3) necessary to be 95 % certain that an unrecorded species is in fact absent from a given site. We performed 1-5 samplings in each site and the species presence we recorded is, obviously, affected by species Table 1. Contribution to species distribution models of amphibians in of the Circeo National Park (central Italy). AUC, that may range from 0 to 1 (null-maximum predictive power of the model) is also shown. Lissotriton vulgaris AUC = 0.995 Occurrence locations = 24 Variable % Contribution Effect / Corine categories Distance from swamps 54.1 Negative Tree Cover 13.1 Positive for values > 50% Corine 12.1 Agricultural areas with significant portions of natural vegetation, Forested areas – Scrubs and herbaceous – Sclerophyllous vegetation, Inland water-bodies – Water courses Triturus carnifex AUC = 0.947 Occurrence locations = 11 Variable % Contribution Effect / Corine categories Corine 50.3 Agricultural areas with significant portions of natural vegetation, Agro-Forestry areas, Forested areas – Scrub and herbaceous – Sclerophyllous vegetation Distance from Swamps 24.3 Negative Distance from Lakes 14.2 Negative Bufo balearicus AUC = 0.974 Occurrence locations = 20 Variable % Contribution Effect / Corine categories Corine 42.6 Agricultural areas with significant portions of natural vegetation, Forested areas, Open spaces with little vegetation Distance from Lakes 36.1 Negative Insolation 11.2 Negative Bufo bufo AUC = 0.993 Occurrence locations = 24 Variable % Contribution Effect / Corine categories Corine 25.7 Forested areas Distance from Lakes 21.9 Negative TRI 20.3 Positive Hyla intermedia AUC = 0.989 Occurrence locations = 23 Variable % Contribution Effect / Corine categories Tree Cover 19.2 Positive Distance from Lakes 18.7 Negative Corine 18 Inland water-bodies - Wetlands, Forested areas Distance from Swamps 17.5 Negative Pelophylax sinkl. esculentus AUC = 0.981 Occurrence locations = 57 Variable % Contribution Effect / Corine categories Distance from Swamps 44.4 Negative Tree Cover 15.6 Positive Distance from Lakes 14.4 Negative Rana dalmatina AUC = 0.997 Occurrence locations = 29 Variable % Contribution Effect / Corine categories Distance from Swamps 75.1 Negative Tree Cover 10.7 Positive for intermediate values 206 Antonio Romano et alii detection probabilities. However in the 40% of the sites sampled by Cinquegranelli et al. (2015) we found at least one additional species (1-4). We think that this discrep- ancy among the species’ ratios could be probably due to the high difference in the number of sampling sites, to different sampling protocols, and to different years (dif- ference in annual precipitation may affect, intuitively, amphibian species detectability) Our original data showed that urodelans have a more limited diffusion and occurrence than anurans. In particu- lar, T. carnifex has been detected only in three sites dur- ing our survey, resulting as the rarest species in the Park (Fig. 3). Its distribution appeared associated with the deep swamps and slow running waters of the plain forest (Fig. 4; Tab. 1), and this datum agrees with the habitat prefer- ences known for this species (Andreone and Marconi, 2006). However, data for the surrounding areas (Novaga et al., 2013) indicate that vernal ponds and marshes are often colonized by both T. carnifex and L. vulgaris, sug- gesting that the limited distribution of newts outside the State forest might also be affected by the occurrence of alien species (Fig. S5), especially P. clarkii (Fig. S6). Bufo balearicus was recorded in open habitat with sandy and clay soils, as typically showed by this species in other Ital- ian areas (Balletto et al., 2007), around the coastal lakes (Fig. 1) where retrodunal ponds and marshes in grazing lands are the most preferred breeding sites (Fig 4; Tab. 1). The concentration of its elective habitats along the coast- al areas explains the overall limited diffusion of this toad in the Park area. On the contrary, B. bufo appears wide- spread in different habitat typologies (Fig. 2), and uses a greater variety of breeding sites (Fig. 4). Hyla intermedia has been detected both in open and forest habitats but the ICS scores showed higher concentrations in the wooded habitats, as could be predicted for a semi-arboreal species. Rana dalmatina can be considered to be the most repre- sentative species of the hygrophilous plain forest and no breeding sites were found outside the forest boundary. Finally, P. sinkl. esculentus is the commonest species in the park (Fig. 2 and 3), as expected for a species which is ecologically plastic and tolerant to anthropic disturbance in a highly urbanised context. During our surveys we did not find S. perspicillata, R. italica and B.pachypus. The question is whether these spe- cies actually occurred in the 60s’ as reported by Bruno (1981) and became extinct in the last decades, or their records reported by this author have to be considered as erroneous. These three species are characterised by dif- ferent ecological requirements, with Salamandrina and R. italica strictly associated to clean running waters, and mainly shady, cool and damp areas (Utzeri et al., 2004, Angelini et al., 2007). Conversly, Bombina pachypus has a realtively wide ecological niche, but it is an heliophilous and thermophilus species tipically linked to open, lentic and shallow waters (Guarino et al., 2007). While for the latter species suitbale habitat were actually identified in the park and its occurrence in past decades cannot be exclud- ed, for the first two amphibians we did not find terrestrial and aquatic habitats matching their ecological require- ments. Suitable breeding sites could presumably be avail- able before land reclamation (20s’ of the XX century) in a restricted piedmont area on the northern slope of the Cir- ceo massif (loc. “sorgente Mezzomonte” and “Rio Torto”). Ecology The full range of aquatic habitats available in the Park is largely exploited by toads, tree frog and pool frog, while the two newts and Rana dalmatina exhibited a narrow habitat niche (Fig. 4). Larger pools are gener- ally deeper in environments similar to that of the CNP (Brooks, 2005) but in the forest we studied this correla- tion was not significant. Wet phase duration is generally correlated with both pond surface area and maximum depth (Schneider and Frost, 1996; Brooks and Hayashi, 2002) but, if these two features are considered indepen- dently, maximum water depth is generally the best pre- dictor of hydroperiod (Calhoun et al., 2003; Skidds and Golet, 2005). The Habitat Suitability Index for R. dalmati- na elaborated by Radiguet (2012) showed that the date of drying of the pond is one of the key components to make habitats highly suitable for this frog. The positive and significant correlation we found between the number of egg clutches of R. dalmatina and maximum depth of the swamps (Fig. 7) is probably related to wet phase duration. Furthermore, as in other Italian areas (e.g., Bernini et al., 2004), also in the CNP Rana dalmatina was strictly asso- ciated with swamps and eggs were preferentially spawned in water bodies with intermediate values of tree cover (Tab. 1; Fig. 4 and 5). Distance from swamps also nega- tively affected the occurrence of both newt species (Tab. 1; Fig. 5). Toads are associated with the environment sur- rounding the lakes (Tab. 1), but the two species differed in their canopy requirements, as the green toad was asso- ciated with open spaces with little vegetation, while the common toad was associated with forested areas. Population abundance estimation Abundance of B. bufo, B. balearicus, H. interme- dia and P. sinkl. esculentus greatly varied among sites, as emerged from the ICS (Fig. 6); the maximum score (3) was recorded for almost all species (2-19% of the surveyed 207Amphibians of the Circeo National Park sites, see Fig. 6), except the common toad. The question of whether Index of Calling Survey (ICS) can be consid- ered a good proxy of actual population size is controver- sial (see for instance Jansen, 2009 and Corn et al., 2011 for articulated discussions on this topic). The reliability of ICS depends on species and is higher for species with loud calls (i.e., higher detectability), as that we studied, and for tree frogs and pool frogs in particular (Pellet and Schmidt, 2005; Tanadini and Schimdt, 2011). For B. bufo, B. balear- icus, H. intermedia and P. sinkl. esculentus in the Circeo National Park, ICS could be easily used in extensive moni- toring programs with relatively low effort. Furthermore, monitoring data collected for ICS may be elaborated using a recent class of statistic models that provide abundance estimations and that consider detectability (N-mixture models; Royle, 2004; Royle and Link, 2005). Fort three species we elaborated demographic esti- mates. Capture-Marking-Recapture (CMR) method estimated about 27 green toads in a population, display- ing  an even sex ratio; a population of pool frogs, using removal method, was estimated to consist of about 60 adults, with strongly male biased sex ratio. In several Italian populations of B. balearicus sex ratio is typically male biased and their Ne are highly variable (see table 1 in Giacoma, 1999). In 26 Italian populations demographic parameters are (range, mean ± s.d., 25th- 75th  percentile): N= 3-292; 89.46 ± 79.66, 28-144.5; sex ratio = 0.53-1, 0.82 ± 0.13, 0.73-0.94, Ne = 0-254.94, 50.27 ± 61.85, 12.84-64.05; (data elaboration from synoptic table 1 in Giacoma, 1999). These two anurans populations showed a contrasting situation if N and Ne are considered independently. One season of data collection provides the size of effectively breeding individuals (Nb) which is directly connected and derived by Ne because Nb times the generation time approximates Ne (Waples, 1990). Thus both Nb and Ne are connected to population’s persistence probability, and may be used as indicators of a popula- tion’s viability (e.g., Frankham et al., 2002). Considering both the high accessibility of the two sites and the ease of the sampling, these two populations of B. balearicus and P. sinkl. esculentus could monitored in be long term to cor- roborate our data and to assess populations trends. Rana dalmatina may be considered the most repre- sentative species of the Park, because this frog is strictly associated with the swamps in the State Forest which rep- resents the residual environment of the pre-reclamation (Tab. 1). Although this frog has a limited distribution it is common (Fig. 3) and abundant (Figs 2d and 7) in its elective environment. The sex ratio in Rana dalmatina is male biased in 90% of the breeding ponds (Lodé et al., 2004; Lodé et al., 2005; Lodé 2009) ranging from about 0.8 to 3. Considering this range of sex ratio, in the seven swamps we studied in the state forest, the whole popula- tion size might approximately range between 2550 and 5700 adults. Threats and conservation strategies Habitat loss and alteration. The CNP is among the most urbanized protected areas at national level. As a consequence, aquatic habitat loss and alteration (which are the main causes of amphibian decline at global level; see Collins, 2010) associated with land development, present the greatest challenge to the persistence of these habitats and their animals. However we did not find sig- nificant evidence of aquatic habitat loss due to current anthropic pressure. The aquatic habitat characterizing the PNC are swamps and ponds with a typical semi- perennial or seasonal hydroperiod. They are lentic shal- low water bodies that are deep, on average, about 60 cm (see results). The main conservation problem of swamps and ponds was an habitat evolution toward dryer situa- tions in shorter time, which is largely due to the dramatic change in the hydrological regime resulting from the past land reclamation. An excavation to increase water depth and the hydroperiod duration was planned for the pond in the locality “Cerreto Fontana” where three amphibian species spawn and that is exposed to a fast drying (see results). Road mortality. Road mortality in the CNP does not seem to be a problem for amphibians and, probably, only one situation (loc. Bufalara) deserves further researches to estimate the actual impact on the green toad popula- tion. Alien species. Invasive alien species are among the key factors threatening biodiversity (EEA, 2012). We found five alien species that, considering the available lit- erature, may be considered as threats to amphibians (Fig. 8, Fig. S5-S6). The red-eared  terrapin Trachemys scripta may have a large negative impact on amphibian popula- tions (tadpoles; Polo-Cavia et al., 2010) in water bod- ies with high numbers of alien turtles. Trachemys scripta seems to be more diffused in the northern part of the Park (Fig. S5a). Its occurrence in the swamps of the State Forest was never recorded, even in the past (A. Romano pers. obs.). We did not find massive aggregations of Tra- chemys, although some sites (e.g., the surroundings of the Fogliano lake) support higher population density than southernmost areas. As a consequence of the apparently low density of Trachemys, we think that its management, for the conservation of amphibians, is less urgent than that of other aliens species. Three alien fishes were detected in small water basin too. The gold fish, Carassius auratus can strongly affect 208 Antonio Romano et alii amphibian populations either by predation at different life stages (e.g., Monello and Wright, 2001) and influenc- ing reproductive behavior (Winandy and Denoël, 2013). However when gold fish does not reach high demograph- ic density, it seems that its presence is compatible with persistence of native amphibians (Hartel, 2004). Therefore to well understand the threat level to amphibians, the occurrence of gold fish in the Park should be evaluated both in the distribution and in population size. The pres- ence data here reported (Fig. S5b), which shows a scat- tered and limited distribution, are the only one recorded during this study, but it should be considered that this fish is more widespread in the park (Zerunian, 1984; Zerunian and Leone, 1996). Gold fish did not occur in the water bodies of the State Forest, probably because these aquatic habitats have a seasonal hydroperiod. The pumpkinseed sunfish, Lepomis gibbosus, was introduced for the first time in Italy in the Varano Lake in 1900 (Central Italy) and experienced an impressive increase of distribution in these last years (Zerunian, 2002). Due to its relatively small size L. gibbosus mainly preys upon amphibian larvae and eggs, but can severely damage also the adults (Hartel et al., 2007). Information about its neg- ative impact on amphibian populations are corroborated from experiments with  controlled conditions (Adams, 2000). L. gibbosus, that we recorded in several ditches, is probably widespread in this park characterised by a high connectivity among linear water bodies. It is worth to mention that L. gibbous was recorded at high density in a large concrete artificial tank (30x5 m, about 2 m depth) on the southern slope of the Cir- ceo massif, (Fig. S5c), which is the area with lowest water habitat availability in the Park. We found a female of Bufo bufo in that tank but no breeding activity was recorded. The artificial tank is a potential habitat for at least three amphibians species, because B bufo, P. sinkl. esculentus, and H. intermedia breed in a similar tank 1.6 km away on the same slope (Fig. 2). The eradication of Lepomis gib- bosus for these sites is planned and will be carried out in 2016 by the Park. Vredenburg (2004) demonstrated that removing introduced fish can enable amphibian popula- tions to recover to pre-decline levels. The mosquitofishes of the genus Gambusia were introduced into natural or artificial water environments in many parts of the world as a  biocontrol  to mosquito populations, in particular where there are (or there were) malaria infections. The effectiveness of this fish in com- bating malaria is still debated and this discussion is out- side the scope of this paper. By the way, it was imported for the first time in Italy, in Pontine marshes on 1922 by G.B. Grassi (Sella, 1926; Ronchetti, 1968) where, as the other Italian areas, is widespread. The introduction of Gambusia sp. strongly depresses all amphibian popula- tions (Adams, 2000; Katz and Ferrer, 2003; Wells, 2007; Segev et al., 2009). The observations in the Park indi- cate widespread presence (Fig. S5d), as the mosquitofish occurs in a least in a quarter of aquatic sites that host amphibians (Fig. 8). The mosquitofish seems to be absent from the water bodies in the State Forest, probably because they have of a seasonal wet regime. Eradication of mosquitofish was planned and carried out in 2015, in an artificial ponds where only pool frogs breed. Around this aquatic site, located in a meadow garden, six rock piles were also placed to offer additional refugia to small vertebrates. After fish removal and the placing of the arti- ficial shelters (November 2015), on February 2016, the pond was colonised by B. bufo, but it is a potential breed- ing site also for tree frogs that are in the surroundings. The pond is in the area devoted to tourist and visitor reception; information panels about the performed con- servation action were also placed. The last alien species we recorded, and the most dangerous one, is the red swamp crayfish, Procambarus clarkii, which is an efficient predator of amphibian larvae of several European species (Gherardi et al, 2001; Cruz and Rebelo, 2005; Cruz et al., 2006a; Ficetola et al., 2011). The presence of this crayfish is a deterrent for the coloni- zation of potentially suitable aquatic habitats by amphib- ians (Cruz et al., 2006a, b; Ficetola et al., 2011). Further- more, P. clarkii is the cause of massive local extinction of amphibians, as happened for instance in a Portugal natu- ral reserve where crayfish caused the disappearance of more than 50% of amphibian species (Cruz et al., 2006b), or, in Italy, the extinction of Rana latastei from one part of its already small distribution range (Mazzotti et al., 2007). It is also a vector of the pathogen Batrachochytrium den- drobatidis which is capable to depressing or to extinct of amphibian populations. In the CNP it is the alien species with more records and widest distribution (Figs 8 and S6) and its occurrence is likely underestimated. Procambarus clarkii is the only alien species reported within the State Forest (Fig. S6), so all amphibians occurring in the forest are severely threatened. Eradication of P. clarkii from the site in the State Forest, is an aim of the Park. To achieve this goal, an integrate strategy is under consideration, using both intensive removal (Hein et al., 2007) and bio- logical control (Aquiloni et al., 2010). An annual monitor- ing to detect further invasions into other areas of the for- est was also planned. A pilot study to test the effectiveness of active removal (i.e., the less expensive method) on the State Forest populations is currently in progress. An updated, large, and geo-referenced database of species distribution is the essential prerequisite for any protected area to effectively manage its resources and 209Amphibians of the Circeo National Park to plan appropriate conservation strategies. In the CNP seven species were recorded, few of them strictly associ- ated with particular aquatic habitats and thus with lim- ited distribution, as Rana dalamatina is. For this species the maintenance of swamps is an essential prerequisite for conservation. The debated issue about the presumed presence of Salamandrina perspicillata, Bombina pachypus and Rana italica can reasonably be regarded as solved: these species are absent in the Park and their closest sites are at least 10 km away, on the Volsci chain. No reliable information are available to assess if they occurred before land reclamation or not. Despite the highly urbanised territory, habitat loss and alteration seem to be limited, and few practical and rather simple actions can be made to improve the current situation. The main threats to amphibians in the park, in our opinion, are the spread of alien species. Particular concern deserves the invasion of the red swamp crayfish in the State plain forest which is the area with highest level of species richness. ACKNOWLEDGEMENT This research was carried out under the project “Pro- getto di Sistema dei Parchi Nazionali Italiani; Action 6: Monitoraggio delle specie  di  ambiente umido  /  acqua- tico”, funded by the Italian Ministry of Environment (Direttiva MATTM ex cap. 1551). Capture permit and manipulation of individuals were approved by the Italian Ministry of Environment with the authorisation number PNM-2015-0016824/PNM. Ester Del Bove (PNC) and Alessandra Noel (Corpo Forestale dello Stato) greatly support this research; Luigi Loffredi contributed to field researches. Thanks to Marta Biaggini for her help in the field study on green toad. The budding naturalist Franc- esco Maria Romano has contributed with great enthusi- asm to sampling activities. We are indebt with two anon- ymous reviewers who greatly improved the ms. SUPPLEMENTARY MATERIAL Supplementary material associated with this article can be found at < http://www.unipv.it/webshi/appendix >. REFERENCES Adams, M.J. (2000): Pond permanence and the effects of exotic vertebrates on anurans. Ecol. Appl. 10: 559-568. Andreone, F., Marconi, M. (2006): Triturus carnifex (Laurenti, 1768). In: Atlante degli Anfibi e dei Rettili d’Italia, pp. 220-225. Sindaco, R., Doria, G., Razzetti, E., Bernini, F., Eds, Ed. Polistampa, Firenze, Italy. Angelini, C., Vanni, S., Vignoli, L. (2007): Salamandrina terdigitata (Bonnaterre, 1789), Salamandrina per- spicillata (Savi, 1821). In: Fauna d’Italia. Vol. XLII. Amphibia, pp. 227-236. Lanza, B., Andreone, F., Bologna, M.A.. Corti, C., Razzetti, E., Eds, Calderini, Bologna, Italy. Aquiloni, L., Brusconi, S., Cecchinelli, E., Tricarico, E., Mazza, G., Paglianti, A., Gherardi, F. (2010): Biologi- cal control of invasive populations of crayfish: the European eel (Anguilla anguilla) as a predator of Pro- cambarus clarkii. Biol. Invasions 12: 3817-3824. Balletto, E., Bologna, M.A., Giacoma, C. (2007): Bufo lineatus Ninni, 1879. In: Fauna d’Italia. Vol. XLII. Amphibia, pp. 299-304. Lanza, B., Andreone, F., Bologna, M.A.. Corti, C., Razzetti, E., Eds, Calderini, Bologna, Italy. Bernini, F., Gentilli, A., Merli, E., Razzetti, E. (2004): Rana dalmatina and R. latastei: Habitat selection, fluctuation in egg clutch deposition and response to exceptional floods in northern Italy, It. J. Zool. 71(S2): 147-149. Blasi, C. (1994): Fitoclimatologia del Lazio. Fitosociologia 27: 151-175. Bologna, M.A., Capula, M., Carpaneto, G.M., Eds, (2000): Anfibi e Rettili del Lazio. Fratelli Palombi Editori, Roma. Brooks, R.T., Hayashi, M. (2002): Depth-area-volume and hydroperiod relationships of ephemeral (vernal) forest pools in Southern New England. Wetlands 22: 247- 255. Brooks, R.T. (2005): A review of basin morphology and pool hydrology of isolated ponded wetlands: impli- cations for seasonal forest pools of the northeastern United States. Wetl. Ecol. Manag. 13: 335-348. Bruno, S. (1973): Anfibi d’Italia: Caudata. Studi sulla fau- na erpetologia italiana – XVII. Natura 64: 209-450. Bruno, S. (1981): Anfibi e Rettili di alcune stazioni del litorale tirrenico tra la foce dell’Arno e il Circeo. In. Ricerche ecologiche, flogistiche e faunistiche sulla fascia costiera mediotirrenica italiana. Quad. Accad. Naz. Lincei 254: 31-76. Calhoun, A.J.K., Walls, T.E., Stockwell, S.S., McCollough, M. (2003): Evaluating vernal pools as a basis for con- servation strategies: a Maine case study. Wetlands 23: 70-81. Carpaneto, G.M. (1986): osservazioni preliminari sugli Anfibi e sui Rettili del Parco Nazionale del Circeo (Amphibia et Reptilia). In: Aspetti faunistici e prob- lematiche zoologiche del Parco Nazionale del Circeo, pp. 145-156. Biondi, M., Ed, Ministero dell’Agricoltura 210 Antonio Romano et alii e delle Foreste, S.G.P. Communication S.r.l., Sabaudia, Italy Cinquegranelli, A., Salvi, D., Vignoli, L. (2015): The amphibians of the Circeo National Park, central Ita- ly: distribution and aquatic habitat use. Herpetozoa 28: 63-73. Collins, J.P.  (2010):  Amphibian decline and extinc- tion:  what we know and what we need to learn. Dis. Aquat. Organ. 92: 93-99. Corn, P.S., Muths, E., Kissel, A.M, Scherer, R.D. (2011): Breeding chorus indices are weakly related to estimat- ed abundance of boreal chorus frogs. Copeia 2011: 365-371. Cruz, M.J., Rebelo, R. (2005): Vulnerability of southwest Iberian amphibians to an introduced crayfish, Pro- cambarus clarkii. Amphibia-Reptilia 26: 293-303. Cruz, M.J., Pascoal, S., Tejedo, M., Rebelo, R. (2006a): Predation by an exotic crayfish, Procambarus clarkii, on natterjack toad, Bufo calamita, embryos: its role on the exclusion of this amphibian from its breeding ponds. Copeia 2006: 274-280. Cruz, M.J., Rebelo, R. and Crespo, E.G. (2006b): Effects of an introduced crayfish, Procambarus clarkii, on the distribution of south-western Iberian amphibians in their breeding habitats. Ecography 29: 329-338. D’Amen, M., Bombi, P., Pearman, P.B., Schmatz, D.R., Zimmermann, N.E. and Bologna, M.A. (2011): Will climate change reduce the efficacy of protected areas for amphibian conservation in Italy? Biol. Conserv. 144: 989-997. Dorcas, M.E., Price, S.J., Walls, S.C., Barichivich, W.J. (2009): Auditory monitoring of anuran populations. In: Amphibian ecology and conservation: a handbook of techniques, pp. 281-298. Dodd, K.Jr., Ed, Oxford University Press, Oxford, UK. Doria, G., Salvidio, S. (1994): Atlante degli anfibi e rettili della Liguria. Regione Liguria, Genova. EEA - European Environment Agency (2012): The impacts of invasive alien species in Europe. EEA Technical Report. ISBN 978-92-9213-345-0 Ficetola, G.F., Siesa, M.E., Manenti, R., Bottoni, L., De Bernardi, F., Padoa-Schioppa, E. (2011): Early assess- ment of the impact of alien species: differential con- sequences of an invasive crayfish on adult and larval amphibians. Diversity Distrib. 17: 1141-1151. Frankham, R., Ballou, J.D., Briscole, D.A. (2002): Intro- duction to conservation genetics. Cambridge Univer- sity Press, Cambridge, United Kingdom. Gherardi, F., Aquiloni, L., Diéguez-Uribeondo, J., Tricar- ico, E. (2011): Managing invasive crayfish: is there a hope? Aquat. Sci. 73: 185-200. Giacoma, C. (1999): Simposio su Bufo viridis: struttura di popolazione. Riv. Idrobiol. 38: 251-265. Giagnacovo, G., Tomassini, A. (2003): Realizzazione del- la Carta dell’Uso del Suolo del Parco Nazionale del Circeo e Aree Limitrofe (Progetto “Parchi in qualità”, “Applicazione pilota del Sistema di Gestione Ambien- tale nelle aree naturali protette”. Accordo di program- ma Ministero dell’Ambiente e del Territorio ed ENEA giugno 1999 − marzo 2003, ENEA. Glista, D.J., DeVault, T.L. De Woody, J.A. (2008): Verte- brate road mortality predominantly impacts Amphib- ians. Herpetol. Conserv. Biol. 3: 77-87. Griffiths, R.A., Raper S.J. (1994): A review of current techniques for sampling amphibian communities. JNCC report N. 210, Joint Nature Conservation Com- mittee, Peterborough. Grossenbacher, K., Lippuner, M., Zumbach, S., Bourgula, A., Luscher, B. (2002): Phenology and reproduction of the three brown frog species Rana latastei, R. dalmati- na, R. temporaria; development and status of the R. latastei populations in Mendrisiotto, Southern Ticino, Switzerland. In: Atti III Conv. Salvaguardia Anfibi, pp. 91-100. Ferri, V., Ed, Mus. Cant. Stor. Nat. Lugano, Cogecstre, Penne. Guarino, F.M., Bellini, L. (1993): Reproductive activity and plasma androgen concentrations in the male of Rana dalmatina. Boll. Zool. 60: 281-286. Guarino, F.M., Picariello, O., Venchi, A. (2007): Bombina pachypus (Bonaparte, 1838). In: Fauna d’Italia. Vol. XLII. Amphibia, pp. 277-280. Lanza B., Andreone, F., Bologna, M.A.. Corti, C., Razzetti, E., Eds, Calderini, Bologna, Italy. Hammer, Ø., Harper, D.A.T., Ryan, P.D. (2001): PAST: Paleontological statistics software package for educa- tion and data analysis. Palaeont. Electr. 4: 1-9. Hartel, T. (2004): The long term trend and the distribu- tion of amphibian populations in a semi-natural pond in the middle section of the Tarnava-Mare Valley (Romania). Biota, J. Biol. Ecol. 5: 25-36. Hartel, T., Nemes, S., Cogalniceanu, D., Ollerer, K., Sch- weiger, O., Moga, C-I., Laszlo, D. (2007): The effect of fish and aquatic habitat complexity on amphibians. Hydrobiologia 583: 173-182. Hein, C. L., Vander Zanden, M. J.andMagnuson, J. J. (2007): Intensive trapping and increased fish preda- tion cause massive population decline of an invasive crayfish. Freshw. Biol. 52: 1134-1146. Jansen, M (2009): Measuring temporal variation in call- ing intensity of a frog chorus with a data logging sound level meter: results from a pilot study in Boliv- ia. Herpetol. Notes 2: 143-149. Jehle, R., Arntzen, J.W., Burke, T., Krupa, A.P., Hödl, W. (2001): The annual number of breeding adults and the effective population size of syntopic newts (Triturus cristatus, T. marmoratus). Mol. Ecol. 10: 839-850. 211Amphibians of the Circeo National Park Kats, L., Ferrer, R. (2003): Alien predators and amphib- ian declines: review of two decades of science and the transition to conservation. Divers. Distrib. 9: 99-110. Lepri, G. (1935):. La fauna pontina prima e dopo la bon- ifica. In: La bonifica delle paludi pontine. pp. 93-102. AA.VV., Eds, Ist. Di Studi Romani. Leonardo da Vin- ci, Casa ed., Roma. Linoli, A. (2005): Twenty-six centuries of reclamation and agricultural improvement on the Pontine marshes. In: Integrated Land and Water Resources Management in History, pp.  27–56. Ohleg, C., Ed., Schriften der Deutschen Wasserhistorischen Gesellschaft (DWhG) Sonderband 2. Lodé, T. (2009): For a few males more: do changes in sex ratio predict reproductive success and offspring sur- vival? Evol. Ecol. Res. 11: 95-107.  Lodé, T., Holveck, M.J., Lesbarrères, D., Pagano, A. (2004): Sex-biased predation by polecats influences the mating system of frogs. Proc. Roy Soc B: Biol. Sci. 271: S399-S401. Lodé, T., Holveck, M.J., Lesbarrères, D. (2005): Asynchro- nous arrival pattern, operational sex ratio and occur- rence of multiple paternities in a territorial breeding anuran, Rana dalmatina. Biol. J. Linn. Soc. 86: 191- 200. Maiorano, L., Falcucci, A., Boitani, L. (2006): Gap analy- sis of terrestrial vertebrates in Italy: priorities for con- servation planning in a human dominated landscape. Biol. Cons. 133: 455-473. Mazzotti, S., Mantovani, R., Penazzi, A., Cavalieri d’Oro, V., Gentile, M., Rossini, L., Lizzio, E., Rizzati, F., Fras- son, V., Mingozzi, A., Noferini, A. (2007): Le comu- nità degli Anfibi del Parco del Delta del Po. In: HERP HELP - Status e strategie di conservazione degli Anfi- bi e dei Rettili del Parco regionale del Delta del Po, pp. 49-58. Mazzotti S., Ed., Quad. Staz. Ecol. civ. Mus. St. nat. Ferrara 17, Ferrara, Italy. Monello, R.J., Wright, R.G. (2001): Predation by goldfish (Carassius auratus) on eggs and larvae of the eastern long-toed salamander (Ambystoma macrodactylum columbianum). J. Herpetol. 35: 350-353. Mossman, M.J., Hartman, L.M., Hay, R., Sauer, J.R., Dhuey, B.J. (1998): Monitoring longterm trends in Wisconsin frog and toad populations. In: Status and conservation of midwestern amphibians, pp. 169-198. Lannoo, M.J., Ed, University of Iowa Press, Iowa City, Iowa. Nollert, A., Nollert, C. (1992): Die Amphibien Europas. Kosmos, Stuttgart. Novaga, R., Corsetti, L., Romano, A. (2013): Dove c’era la palude: batracofauna urbana della città di Latina (Lazio, Italia centrale). In: Atti IX Congresso Nazion- ale della Societas Herpetologica Italica, pp. 163-167. Scillitani, G., Liuzzi, C., Lo russo, L., Mastropasqua, F., Ventrella, P., Eds, Ed. Pineta, Conversano, Italy. Otis, D., Burnham, K.P., White, G.C., Anderson, D.R. (1978): Statistical inference from capture data on closed animal populations. Wildl. Monogr. 62: 1-135. Padula, M. (1985): Aspetti della vegetazione del Parco Nazionale del Circeo. Webbia 39: 29-110. Pellet, J., Schmidt, B.R. (2005): Monitoring distributions using call surveys: estimating site occupancy, detec- tion probabilities and inferring absence. Biol. Con- serv. 123: 27-35. Phillips, S.J., Anderson, R.P., Schapire, R.E. (2006): Maxi- mum entropy modeling of species geographic distri- butions. Ecol. Model. 190: 231-259. Polo-Cavia, N., Gonzalo, A., Lopez, P., Martin, J. (2010): Predator recognition of native but not invasive turtle predators by naïve anuran tadpoles. Anim. Behav. 80: 461-466. Pollock, K.H, Otto, M.C. (1983): Robust estimation of pop- ulation size in closed animal populations from capture- recapture experiments. Biometric 39: 1035-1049 Puky, M. (2006): Amphibian road kills: a global perspec- tive. In: Proceedings of the 2005 International Con- ference on Ecology and Transportation, pp. 325-338. Irwin, C.L., Garrett, P., McDermott, K.P., Eds., Center for Transportation and the Environment, North Caro- lina State University (Raleigh, NC, USA). Radiguet, F. (2012): The creation of a Habitat Suitabil- ity Index for the Agile frog, Rana dalmatina, in Jer- sey. BSc Thesis, Durrell Institute of Conservation and Ecology, University of Kent. Ravenna, C. (2013): Circeo - Il parco nazionale - Natura, storia, itinerari guidati. Edizioni Belvedere, Latina. Riley, S.J., De Gloria, S.D., Elliot, R. (1999): A terrain ruggedness index that quantifies topographic hetero- geneity. Intermount. J. Sci. 5: 23-27. Romano, A., Ventre, N, De Riso, L., Pignataro, C., Spilin- ga, C. (2010): Amphibians of the “Cilento e Vallo di Diano” National Park (Campania, Southern Italy): update check list, distribution and conservation notes. Acta Herpetol. 5: 233-244. Romano, A., Bartolomei, R., Conte, L.A. and Fulco, E. (2012): Amphibians in Southern Apennine: distribu- tion, ecology and conservation notes in the “Appen- nino Lucano, Val d’Agri e Lagonegrese” National Park (Southern Italy). Acta Herpetol. 7: 203-219. Ronchetti, G. (1968): L’azione antianofelica dei pesci del genere Gambusia, utilizzati per la lotta biologica con- tro la malaria. Natura, LVIII: 25-41. Royle, J.A. (2004): Modeling abundance index data from anuran calling surveys. Cons. Biol. 18: 1378-1385. 212 Antonio Romano et alii Royle, J.A., Link, W.A (2005): A general class of multino- mial mixture models for anuran calling survey data. Ecology 86: 2505-2512. Schmeller, D.S, Merila, J. (2007): Demographic and genetic estimates of effective population and breeding size in the amphibian Rana temporaria. Cons. Biol. 21: 142-151. Schneider, D.W., Frost, T.M. (1996): Habitat duration and community structure in temporary ponds. J. N. Am. Benthol. Soc. 15: 64-86. Segev, O., Mangel, M., Blaustein, L. (2009): Deleteri- ous effects by mosquitofish (Gambusia affinis) on the endangered fire salamander (Salamandra infraimmac- ulata). Anim. Conserv. 12: 29-37. Sella, M. (1926): A new means of combating Anopheles in Italy: a account of the acclimatization and progress of Gambusia. In: C. R. Cong. Int. Paludisme, Rome. Skidds, D.E, Golet, F.C. (2005): Estimating hydroperiod suitability for breeding amphibians in southern Rhode Island seasonal forest ponds. Wetl. Ecol. Manag. 13: 349-366. Sofianidou, T.S., Kyriakopoulou-Sklavounou, P. (1983): Studies on the biology of the frog Rana dalmatina Bonaparte during the breeding season in Greece (Amphibia: Anura: Ranidae). Amphibia-Reptilia 4: 125-136. Sörensen, R., Zinko, U., Seibert, J. (2006). On the calcu- lation of the topographic wetness index: evaluation of different methods based on field observations. Hydrol. Earth Syst. Sci. Discus. 10: 101-112. Stuart, S.N., Chanson, J.S., Cox, N.A., Young, B.E., Rodri- gues, A.S.L., Fischman, D.L., Waller, R.W. (2004): Sta- tus and trends of amphibian declines and extinctions worldwide. Science 306: 1783-1786. Tanadini; L.G., Schmidt, B.R. (2011): Population size influences amphibian detection probability: implica- tions for biodiversity monitorino. Programs. PLoS ONE 6: e28244. Turrisi, G., Vaccaro, A. (2004): Status and conservation of herpetofauna from the Iblean area (south eastern Sic- ily). It. J. Zool. 71 (s2): 185-189. Utzeri, C., Antonelli, D., Angelini, C. (2004): A note on terrestrial activity and feeding in the Spectacled sala- mander, Salamandrina terdigitata (Urodela, Salaman- dridae). Herpetol. Bull. 90: 27-31. Vredenburg, V.T. (2004): Reversing introduced species effects: Experimental removal of introduced fish leads to rapid recovery of declining frog. P. Natl., Acad. Sci. USA 101: 7646-7650. Wake, D.B., Vredenburg, V.T. (2008): Are we in the midst of the sixth mass extinction? A view from the world of amphibians. P. Natl., Acad. Sci. USA 105: 11466- 11473. Waples, R. S. (1990): Conservation genetics of Pacific salmon. II. Effective population size and the rate of loss of genetic variability. J. Hered. 81: 267-276. Weir, L.A., Mossman, M.J. (2005): North American Amphibian Monitoring Program (NAAMP). In: Amphibian declines: conservation status of United States species, pp. 307–313. Lanoo, M.J., Ed, Univer- sity of California Press, Berkeley. Wells, K.D. (2007): The Ecology and behavior of amphib- ians. The University of Chicago Press. Wilson, K., Hardy, I.C.W. (2002): Statistical analysis of sex ratios: an introduction. In: Sex Ratios – Concepts and Research Methods, pp. 48-92. Hardy, I.C.W., Ed., Cambridge, Cambridge University Press. Winandy, L., Denoël, M. (2013): Introduced Goldfish affect amphibians through inhibition of sexual behav- iour in risky habitats: an experimental approach. Plo- sOne 8: e82736. Wright, S. (1938): Size of population and breeding struc- ture in relation to evolution. Science 87: 430-431. Zerunian, S. (1984): I Pesci del Fiume Amaseno e dei corsi d’acqua della Pianura Pontina (Lazio). Quad. Ist. Idrobiol. Acquacolt. Brunelli 4: 26-67.  Zerunian, S. (2002): Condannati all’estinzione? Biodiver- sità, biologia, minacce e strategie di conservazione dei pesci d’acqua dolce indigeni in Italia. Bologna: Il Sole 24 Ore Ed. Agricole. Zerunian, S., Leone, M., Eds, (1996): Monitoraggio delle acque interne e Carta ittica della Provincia di Latina: i bacini campione del Fiume Amaseno e del Lago di Fondi. Amm. Prov. Latina. Acta Herpetologica Vol. 11, n. 2 - December 2016 Firenze University Press Predator-prey interactions between a recent invader, the Chinese sleeper (Perccottus glenii) and the European pond turtle (Emys orbicularis): a case study from Lithuania Vytautas Rakauskas1,*, Rūta Masiulytė1, Alma Pikūnienė2 Effective thermoregulation in a newly established population of Podarcis siculus in Greece: a possible advantage for a successful invader Grigoris Kapsalas1, Ioanna Gavriilidi1, Chloe Adamopoulou2, Johannes Foufopoulos3, Panayiotis Pafilis1,* The unexpectedly dull tadpole of Madagascar’s largest frog, Mantidactylus guttulatus Arne Schulze1,*, Roger-Daniel Randrianiaina2,3, Bina Perl3, Frank Glaw4, Miguel Vences3 Thermal ecology of Podarcis siculus (Rafinesque-Schmalz, 1810) in Menorca (Balearic Islands, Spain) Zaida Ortega*, Abraham Mencía, Valentín Pérez-Mellado Growth, longevity and age at maturity in the European whip snakes, Hierophis viridiflavus and H. carbonarius Sara Fornasiero1, Xavier Bonnet2, Federica Dendi1, Marco A.L. Zuffi1,* Redescription of Cyrtodactylus fumosus (Müller, 1895) (Reptilia: Squamata: Gekkonidae), with a revised identification key to the bent-toed geckos of Sulawesi Sven Mecke1,*,§, Lukas Hartmann1,2,§, Felix Mader3, Max Kieckbusch1, Hinrich Kaiser4 The castaway: characteristic islet features affect the ecology of the most isolated European lizard Petros Lymberakis1, Efstratios D. Valakos2, Kostas Sagonas2, Panayiotis Pafilis3,* Sources of calcium for the agamid lizard Psammophilus blanfordanus during embryonic development Jyoti Jee1, Birendra Kumar Mohapatra2, Sushil Kumar Dutta1, Gunanidhi Sahoo1,3,* Mediodactylus kotschyi in the Peloponnese peninsula, Greece: distribution and habitat Rachel Schwarz1,*, Ioanna-Aikaterini Gavriilidi2, Yuval Itescu1, Simon Jamison1, Kostas Sagonas3, Shai Meiri1, Panayiotis Pafilis2 Swimming performance and thermal resistance of juvenile and adult newts acclimated to different temperatures Hong-Liang Lu, Qiong Wu, Jun Geng, Wei Dang* Olim palus, where once upon a time the marsh: distribution, demography, ecology and threats of amphibians in the Circeo National Park (Central Italy) Antonio Romano1,*, Riccardo Novaga2, Andrea Costa1 On the feeding ecology of Pelophylax saharicus (Boulenger 1913) from Morocco Zaida Ortega1,*, Valentín Pérez-Mellado1, Pilar Navarro2, Javier Lluch2 Notes on the reproductive ecology of the rough-footed mud turtle (Kinosternon hirtipes) in Texas, USA Steven G. Platt1, Dennis J. Miller2, Thomas R. Rainwater3,*, Jennifer L. Smith4 Heavy traffic, low mortality - tram tracks as terrestrial habitat of newts Mikołaj Kaczmarski*, Jan M. Kaczmarek Book Review: Sutherland, W.J., Dicks, L.V., Ockendon, N., Smith, R.K. (Eds). What works in conservation. Open Book Publishers, Cambridge, UK. http://dx.doi.org/10.11647/OBP.0060 Sebastiano Salvidio