Research Article |
Corresponding author: Jelka Crnobrnja-Isailović ( jelka@pmf.ni.ac.rs ) Academic editor: Günter Gollmann
© 2019 Bogdan Jovanović, Jelka Crnobrnja-Isailović.
This is an open access article distributed under the terms of the Creative Commons Attribution License (CC BY 4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Citation:
Jovanović B, Crnobrnja-Isailović J (2019) Fluctuations in population abundance in two anurans from Central Serbia. Herpetozoa 32: 65-71. https://doi.org/10.3897/herpetozoa.32.e35660
|
We monitored the population size of the agile frog (Rana dalmatina) and the common toad (Bufo bufo) from 2011 or 2012, respectively, to the year 2017 at a syntopic breeding site in the vicinity of Belgrade. Adult R. dalmatina population size had minor fluctuations during the years of study (from 351 to 108 frogs). On the contrary, the adult B. bufo population was widely fluctuating towards decline (from 1158 to 141 toads). In both species, population fluctuations were not significantly related to variation of meteorological parameters (air temperature, humidity, precipitation). Density dependence effects on population size were not detected in either species. Apart from possible effects of climate change, the indicated trend towards decline of the monitored B. bufo population could also be the outcome of common population fluctuations or of increasing anthropogenic impact (vicinity of settlement and agricultural land). More years of monitoring more than one population are required to obtain precise information. Nevertheless, our results seem to be coherent with other studies that recommend conservation action for this species.
Bufo bufo, climate variation, population growth rate, population size variation, Rana dalmatina, south-eastern Europe
Many researchers all over the world report evidence on amphibian decline (
Monitoring is a key activity for establishing amphibian population trends and planning conservation actions. However, monitoring actions are more typically carried out on rare, endangered species within protected areas rather than on common species from human-altered habitats and which can lead to the extinction of local populations before any action is taken (see in
Rana dalmatina is widespread in Southern Europe. It is classified as “Least Concern” by the IUCN (
Both species are widespread in Serbia and are strictly protected by national legislation (regulation on the proclamation and protection of strictly protected and protected wild species of plants, animals and fungi – appendix I,
The study site is a permanent pond on the outskirts of Belgrade (Zuce village, 44°40.93’N, 20°33.12’E, altitude 240 m) in Central Serbia (
We analysed data starting from 2012 to 2017 and from 2011 to 2017 for the R. dalmatina and B. bufo, respectively. We made several visits from 2011–2015 to the site from early March onwards. From 2015, visits became more frequent e.g. on a weekly level, starting in January or February. For R. dalmatina, spawning begins one or two weeks before the mating period of B. bufo and ends within this period. We counted egg clutches of R. dalmatina and adults of both species regularly throughout the breeding season and the number of counting days varied from 2 to 20 (10.4 on average). On each visit, we inspected the pond and recorded anuran presence by walking slowly along the bank, always in a clockwise direction from a fixed starting point. From 2015, for R. dalmatina, we counted the number of egg clutches per day until all had hatched, while for the B. bufo, we performed regular daily scanning of reproductive activity when the number of amplexuses seen in the pond reached 100. Visits became less frequent when the number of toads started to decrease significantly but were not stopped until the last adult toad left the pond. We searched for adults and egg clutches across the whole surface of the pond, but found the vast majority of both near the pond edge.
We collected meteorological data used in this study from the Republic Hydrometeorological Institute of Serbia (http://www.hidmet.gov.rs) with an included period from 1 December to 28 or 29 February, except data on mean daily precipitation for the period 01.01.2015 – 19.01.2015 which are missing. In the raw database, the mean daily humidity was sometimes slightly higher than 100%, but we treated all those records in our analyses as 100%.
For R. dalmatina, we estimated an annual adult population size (N) as the maximum number of egg clutches counted during the reproductive season of a particular year following
The set of meteorological data included annual average (AVG), minimum (MIN) and maximum (MAX) values of mean daily temperature (T, °C), humidity (H, %) and precipitation (P, mm) for three winter months (December, January and February) for each year (set: AVGT, AVGH, AVGP, MINT, MINH, MINP, MAXT, MAXH, MAXP). We chose meteorological data from that part of the year for comparison with the fluctuation of population size in the two analysed species because winter time is very important for anuran survival in the temperate climate zone: frogs and toads spend those months in hibernation and the quality of hibernation is reflected in the body condition of post-hibernating individuals (see in
We also made an estimate of population growth rate for both species following
The number of R. dalmatina egg clutches varied from 351 in year 2012 to 108 in year 2015 (Table
Population parameters, meteorological parameters, and descriptive statistics of variables for Rana dalmatina population.
Year | Population parameters | Meteorological parameters | |||||||||
N | Δ N | Temperature °C | Humidity % | Precipitation mm | |||||||
AVG | MIN | MAX | AVG | MIN | MAX | AVG | MIN | MAX | |||
2012 | 351 | -0.32 | 1.66 | -12.10 | 14.90 | 76.01 | 52.10 | 98.00 | 2.16 | 0.00 | 36.80 |
2013 | 254 | -0.71 | 3.25 | -4.80 | 13.30 | 78.31 | 47.80 | 94.90 | 2.06 | 0.00 | 18.50 |
2014 | 124 | -0.14 | 5.36 | -5.90 | 16.90 | 75.62 | 45.00 | 98.20 | 0.58 | 0.00 | 8.90 |
2015 | 108 | 0.91 | 4.17 | -9.20 | 11.20 | 76.94 | 49.60 | 97.60 | 2.01 | 0.00 | 19.00 |
2016 | 271 | -0.76 | 4.95 | -9.20 | 15.00 | 79.44 | 38.00 | 100.00 | 1.03 | 0.00 | 12.80 |
2017 | 126 | 0.42 | -13.10 | 16.40 | 76.50 | 31.00 | 100.00 | 0.66 | 0.00 | 17.60 | |
Descriptive statistics | |||||||||||
Valid Number of years | 6 | 5 | 6 | 6 | 6 | 6 | 6 | 6 | 6 | 6 | 6 |
Mean | 206 | -0.20 | 3.30 | -9.05 | 14.62 | 77.14 | 43.92 | 98.12 | 1.42 | 0.00 | 18.93 |
MIN | 108 | -0.76 | 0.42 | -13.10 | 11.20 | 75.62 | 31.00 | 94.90 | 0.58 | 0.00 | 8.90 |
MAX | 351 | 0.91 | 5.36 | -4.80 | 16.90 | 79.44 | 52.10 | 100.00 | 2.16 | 0.00 | 36.80 |
SD | 100.28 | 0.68 | 1.94 | 3.28 | 2.10 | 1.46 | 7.97 | 1.88 | 0.74 | 0.00 | 9.59 |
In the analysed B. bufo population, the total number of adult individuals, as well as the number of adult males, seemed to be declining over the monitoring years (Table
Population parameters, meteorological parameters, and descriptive statistics of variables for Bufo bufo population.
Year | Population parameters | Meteorological parameters | ||||||||||||
N | Nmmax | Nfmax | SR | ΔN | Temperature °C | Humidity % | Precipitation mm | |||||||
AVG | MIN | MAX | AVG | MIN | MAX | AVG | MIN | MAX | ||||||
2011 | 950 | 885 | 65 | 14:1 | 0.20 | 1.69 | -8.10 | 16.60 | 80.01 | 56.30 | 95.90 | 1.83 | 0.00 | 27.20 |
2012 | 1158 | 1034 | 124 | 8:1 | -0.40 | 1.66 | -12.10 | 14.90 | 76.01 | 52.10 | 98.00 | 2.16 | 0.00 | 36.80 |
2013 | 777 | 706 | 170 | 4:1 | -0.32 | 3.25 | -4.80 | 13.30 | 78.31 | 47.80 | 94.90 | 2.06 | 0.00 | 18.50 |
2014 | 566 | 551 | 59 | 9:1 | 0.50 | 5.36 | -5.90 | 16.90 | 75.62 | 45.00 | 98.20 | 0.58 | 0.00 | 8.90 |
2015 | 933 | 871 | 62 | 14:1 | -1.88 | 4.17 | -9.20 | 11.20 | 76.94 | 49.60 | 97.60 | 2.01 | 0.00 | 19.00 |
2016 | 141 | 139 | 8 | 17:1 | 0.65 | 4.95 | -9.20 | 15.00 | 79.44 | 38.00 | 100.00 | 1.03 | 0.00 | 12.80 |
2017 | 270 | 264 | 6 | 44:1 | 0.42 | -13.10 | 16.40 | 76.50 | 31.00 | 100.00 | 0.66 | 0.00 | 17.60 | |
Descriptive statistics | ||||||||||||||
Valid number of years | 7 | 7 | 7 | 7 | 6 | 7 | 7 | 7 | 7 | 7 | 7 | 7 | 7 | 7 |
Mean | 685 | 636 | 71 | 16:1 | -0.21 | 3.07 | -8.91 | 14.90 | 77.55 | 45.69 | 97.80 | 1.48 | 0.00 | 20.11 |
MIN | 141 | 139 | 6 | 4:1 | -1.88 | 0.42 | -13.10 | 11.20 | 75.62 | 31 | 94.90 | 0.58 | 0.00 | 8.90 |
MAX | 1158 | 1034 | 170 | 44:1 | 0.65 | 5.36 | -4.80 | 16.90 | 80.01 | 56.30 | 100 | 2.16 | 0.00 | 36.80 |
SD | 375.55 | 334.95 | 59.25 | 13.23 | 0.92 | 1.87 | 3.01 | 2.06 | 1.72 | 8.65 | 1.91 | 0.69 | 0.00 | 9.29 |
The population of R. dalmatina, analysed in this study, can be considered as stable during the monitoring period of six years, despite minor fluctuations in size. On the contrary, the size of the syntopic B. bufo population apparently varied between consecutive years. Fluctuation patterns of this population size suggested its decline. This B. bufo population might be threatened by extinction as the number of adult members (especially females) was very low in the last two years of the study.
Literature data revealed a variety of information on numbers of egg clutches in R. dalmatina: they varied from 6 to 82 (average of 26) and from 7 to 68 (average of 37) at two breeding sites in Austria (
Literature data show that B. bufo population size estimates varied from 238 to 614 females and from 532 to 1045 males in Sweden (
The two syntopic anuran populations from this study are situated in the vicinity of Belgrade, the capital of Serbia. Continuous urbanisation of the city has also been affecting villages close to this amphibian breeding place by causing inevitable fragmentation and pristine habitat loss over time (see in Vegetation Continuous Fields – VCF,
Another possible cause of the B. bufo population decline could be climate change but our data include too short a period in time for testing relationships between changes in population size and meteorological data. Literature data provide a number of examples on how autumn and winter temperature affect amphibians, from positive effects of higher temperatures on the agile frog (
We are indebted to Günter Gollmann and two anonymous reviewers for comments and suggestions that significantly improved the quality of this manuscript. We are also grateful to the Hydrometeorological Institute of Serbia for sharing meteorological data and to our dear colleague Mario Schweiger for help with the literature. Mrs Esther Helaizen improved the quality of English language. This study was funded by the Ministry of Education, Science and Technological Development of the Republic of Serbia, Grant 173025. Permits for fieldwork were issued by the Ministry of Environment and Spatial Planning of Republic of Serbia No.353–01–29/2011–03, Ministry of Environment, Mining and Spatial Planning of Republic of Serbia No. 353–01–505/2012–03, Ministry of Energetics, Development and Nature Protection No.353–01–54/2013–08, No. 353–01–312/2014–08, No. 353–01–42/2014 Ministry of Agriculture and Nature Protection No. 353–01–170/2016–17 and No.353–01–2666/2016–17.