# data[].attributes.body
Geographically, Greece belongs to Southeastern Europe and to the Eastern Mediterranean region. It is the crossroad of three continents, namely, Europe, Asia, and Africa.

The geomorphology is diverse. Greece is primarily a mountainous country, with seventy per cent of its territory covered by mountains. It has a very long coastline, with a plethora of peninsulas and islands.

The ecosystems range from the semi-desert like ones of Southeastern Crete to the cold climate mountain forests of birch, Scotch pine, and spruce. This variety often appears within small areas. For example, in a 150 km transect from the coasts of the North Aegean to the top of the Rodopi mountain one encounters ecosystems of Mediterranean, Central European and Scandinavian type.
Greece has an area of 131,940 km².
So far there are 69 species of reptiles (of which 10 are endemic) and 25 species of amphibians (of which the three are endemic).
It is important to say that the total number of reptiles hosted by Greece, per km², is the largest in Europe!
https://www.geogreece.gr/refuge_en.php

# data[].attributes.headline
Some herpetological trips (in Greece) with lists of observed species

# data[].attributes.subject
Some herpetological trips (in Greece) with lists of observed species

# data[].attributes.body
-https://www.hylawerkgroep.be/jeroen/index.php?id=8
-https://www.hylawerkgroep.be/jeroen/index.php?id=66
-https://www.hylawerkgroep.be/jeroen/index.php?id=49
-https://www.hylawerkgroep.be/jeroen/index.php?id=28
-https://www.hylawerkgroep.be/jeroen/index.php?id=7

# data[].attributes.attachment_link.description
Project Noah is a tool that nature lovers can use to explore and document local wildlife and a common technology platform that research groups can use to harness the power of citizen scientists everywhere.

# data[].attributes.body
Lake Kaiafa (Fig. A) and waterfalls of Polylimnio (Fig. B) are typical wetlands of Greece.
The former, is a lake in Peloponnese, in Ilia. The lake is located by the sea, it has thermal waters and water sports activities such as water skiing. It is administratively owned by Zacharo, which is five (5) kilometers away.
The latter, is a combination of waterfalls and small lakes in Peloponnese, in Messinia. It is 14 kilometers away from the nearest destination, Petalidi town.
In these kind of habitats, you can find many species of herpetofauna.
Like Emys orbicularis (1), Podarcis peloponnesiacus (2), Mauremys rivulata (3), Hellenolacerta graeca (4), Rana graeca (5) and Natrix natrix (6).
https://www.herpsafari.nl/tripreports-eu/2017-7-greece-peloponnese-pindos-kythira-pori/

# data[].attributes.attachment_link.title
Καταγραφή.PNG

# data[].attributes.body
Phryganic ecosystems are one of the main plant formations formed especially in southern Greece and the islands, up to the altitude of 300 meters. Phrygana covers a large part of Greece and, in general, the Mediterranean has the appropriate climatic conditions for their development. This type of vegetation is typical of the Mediterranean ecosystems and is usually considered to be the result of degradation of pre-existing vegetation, maquis and forest. It grows mainly on poor and rocky limestone and siliceous soils or on areas repeatedly burned by fires.
Οn the island of Milos, there is a large part of the phryganic ecosystems (Fig. A) where they host several species of herpetofauna.
Such as Lacerta trilineata (1), Podarcis milensis (2), Macrovipera schweizeri (3), Telescopus fallax (4) and Elaphe situla (5).
https://www.milosminingmuseum.com/en/the-wealth-of-milos/the-natural-wealth-of-milos/
https://www.youtube.com/watch?v=vcpuFbg5n_0

# data[].attributes.attachment_link.title
Καταγραφή.PNG

# data[].attributes.body
Parks in the cities are something like islands. Only instead of the sea, they are surrounded by cement and roads. So, the species found in the city parks are ensnared. Athens, is a city with several parks where there are several species of herpetofauna.
With its intense presence of stones and vegetation, Athens parks are capable of hosting both reptiles and amphibians.
The presence of people, of course, exacerbates the state of the parks, but some of the species seem to manage it.
So, in Athens parks (Fig. A) there are Testudo marginata (1), Bufotes viridis (2), Pelophylax kurtmuelleri (3) and Testudo hermanni (4).
https://www.inaturalist.org/observations?place_id=any&verifiable=any&iconic_taxa=Amphibia,Reptilia

# data[].attributes.attachment_link.title
Καταγραφή.PNG

# data[].attributes.attachment_link.description
Uetz, P. & Jirí Hošek (eds.), The Reptile Database, http://www.reptile-database.org, accessed Dec 8, 2013

# data[].attributes.body
http://www.herpetofauna.gr/

# data[].attributes.body
Remains of a varanid lizard from the middle Pleistocene of the Tourkobounia 5 locality near Athens, Greece are described. The new material comprises cranial elements only (one maxilla, one dentary, and one tooth) and is attributed to Varanus (Varanus marathonensis), the genus to which all European Neogene varanid occurrences have been assigned. Previously, the youngest undisputed varanid from Europe had been recovered from upper Pliocene sediments. The new Greek fossils therefore constitute the youngest records of this clade from the continent.
These fossils could represent a survivor of the monitor lizards of Asian origin that inhabited Europe during the Neogene.
https://www.tandfonline.com/doi/abs/10.1080/02724634.2017.1301946
-http://reptileevolution.com/index.htm
-https://journals.plos.org/plosone/article?id=10.1371%2Fjournal.pone.0207719&fbclid=IwAR2hpevPgpGybzSmEb9BD6DStdX5v3XkjtJhNT-pZIgwTtHZwkZ1ANOZ0Zg

# data[].attributes.attachment_link.title
Καταγραφή.PNG

# data[].attributes.body
The presence of the frog Latonia cf. gigantea in Greece is documented for the first time. Additionally, the presence of viperids in Aliveri implies a much wider distribution for these snakes during the early Miocene of Europe. Of special interest is the presence of a peculiar colubrid that seems to possess a hitherto unknown vertebral structure, which is herein defined as the ‘paracentral ridge’. Although incomplete, the new material has important taxonomic and biogeographic implications, as it enhances our understanding of southeastern European herpetofaunas from the early Miocene, a time period that was characterised by major dispersal and extinction events and climatic change that affected the whole continent.
Unfortunately, there is not enough information about the exact location of the species, Latonia cf. gigantea, on the phylogenetic tree.
https://www.tandfonline.com/doi/abs/10.1080/08912963.2017.1417404
-https://www.sciencedirect.com/science/article/pii/S0016699516300973
-https://www.sciencedirect.com/science/article/pii/S0016699594800588

# data[].attributes.attachment_link.title
Καταγραφή.PNG

# data[].attributes.body
During the glacial period all species migrated due to the low temperature prevailing in northern Europe. The southern Europe was not ice-covered and thus created favorable conditions for the species. So, many species settled in the southern part of Europe and created populations.
In these areas they found plenty of water, food and managed to survive.
Over the years and due to the geographical relief of southern Europe, but also a number of other factors (such as microclimate, the natural selection etc.) new species were created.
Non-migrating species were driven to extinction.
To conclusion, I do believe that these two phenomena are connected.
https://onlinelibrary.wiley.com/doi/full/10.1111/ele.12853
-https://blogs.egu.eu/divisions/cr/2016/03/
-http://ec.europa.eu/environment/nature/conservation/species/redlist/reptiles/geographic_patterns.htm
-
https://www.youtube.com/watch?v=_SuRFwrJZ9s

# data[].attributes.attachment_link.title
Καταγραφή.PNG

# data[].attributes.headline
The effect of climate change on organisms in Europe

# data[].attributes.subject
The effect of climate change on organisms in Europe

# data[].attributes.body
Over the past 40 years, species have been extending their ranges toward the poles and populations have been migrating, developing, or reproducing earlier in the spring than previously. These range expansions and changes in the timing of seasonal events have generally been attributed to "phenotypic plasticity".
Exept the "phenotypic plasticity" recent studies show that over the recent decades, climate change has led to heritable, genetic changes in populations of animals.
Blackcaps (birds) in central Europe have been increasingly overwintering in Britain rather than Iberia; the genetically distinct British subpopulation arrives earlier at the nesting grounds and thus obtains superior territories or mates. European great tits (birds) depend on caterpillars to feed their young. With earlier springs, the caterpillars have been maturing earlier, before the tit chicks hatch, leading to a decline in lifetime reproductive success of the birds.
Global warming is proceeding fastest at the most northern latitudes, where the gradient in winter cold is steepest, thereby expanding the growing season while alleviating winter cold stress without imposing summer heat stress. Hence, northern climates are becoming more like those in the south, resulting in many species in the future to emigrate to the northern part of Europe.
This would lead to disturbance of ecosystems because ecosystems are completely dependent on the diversity of species that may change in the future.
https://www.nature.com/articles/nature01286
-https://royalsocietypublishing.org/doi/10.1098/rspb.2000.1065
-http://www.cambridge.org/us/catalogue/catalogue.asp?isbn=0521864712
-https://environmentalmigration.iom.int/maps
-
https://www.eea.europa.eu/soer-2015/europe/climate-change-impacts-and-adaptation

# data[].attributes.attachment_link.title
Καταγραφή.PNG

# data[].attributes.headline
The amphibian and reptile species of Greece

# data[].attributes.subject
The amphibian and reptile species of Greece

# data[].attributes.body
As I said in "Welcome to my HerpetoArea" section, Greece has the largest number of reptiles in Europe because of the intense geographic relief and climatic conditions. There are 69 species of reptiles (of which 10 are endemic) and 25 species of amphibians (of which the three are endemic).
Below, I mention a site displaying all reptiles and amphibians in Greece.
In the icon there is a map formed by myself where all the endemic species of Greece appear.
http://www.herpetofauna.gr/
-https://www.cambridge.org/core/journals/journal-of-zoology/article/phylogeny-of-the-european-lizard-genus-algyroides-reptilia-lacertidae-based-on-dna-sequences-with-comments-on-the-evolution-of-the-group/9F5487CD8541A03F9E2E08F156E2F609
-https://www.tandfonline.com/doi/abs/10.1017/S1477200008002727
-https://www.researchgate.net/profile/Panayiotis_Pafilis/publication/328307210_Diet_composition_of_the_Karpathos_marsh_frog_Pelophylax_cerigensis_what_does_the_most_endangered_frog_in_Europe_eat/links/5bc5c309299bf17a1c559e26/Diet-composition-of-the-Karpathos-marsh-frog-Pelophylax-cerigensis-what-does-the-most-endangered-frog-in-Europe-eat.pdf
-https://www.sciencedirect.com/science/article/pii/S1055790317304979
-
https://www.mdpi.com/1424-2818/2/2/233

# data[].attributes.attachment_link.title
Καταγραφή.PNG

# data[].attributes.headline
New records and an updated list of the herpetofaunaof Kastellorizo and the adjacent islet Psomi(Dodecanese, SE Greece)

# data[].attributes.subject
New records and an updated list of the herpetofaunaof Kastellorizo and the adjacent islet Psomi(Dodecanese, SE Greece)

# data[].attributes.body
On Psomi islet, encountered individuals of Anatololacerta sp. and Mediodactylus cf. danilewskii. Based on the analysis of the mitochondrial gene cytb, the recently discovered lizard population from the settlement of Megisti was identified as Anatololacerta pelasgiana (Mertens, 1959), while that of Psomi islet as Anatololacerta budaki (Eiselt and Schmidtler, 1986), a new addition to the herpetofauna of Greece.

https://www.biotaxa.org/hn/article/view/38679/36565

# data[].attributes.attachment_link.title
Καταγραφή.PNG

# data[].attributes.headline
Physiology and adaptation of of reptile species

# data[].attributes.subject
Physiology and adaptation of of reptile species

# data[].attributes.body
Podarcis erhardii livadiacus (TL: 22cm; SVL: 7cm), is less of a climber than many other wall lizard species, but climbs more than some others (such as Podarcis tauricus). This lizard occurs, on islands running between low shurbs (from sea-level to medium altitudes, up to 1,400 m). Usually, an inquiditive species and less shy than some of its congeners and feeds on a wide variety of invertebrates.
On some islands in Aegean lives alongside Falco eleonarae. These lizards gravitate towards falcon nests when they breed in lat summer, feeding on discarded prey, flies attracted to the birds' excrement and possibly also on external parasites of the chicks.
Active throughout the year in S of range (because of the high temperature) and in N of range active between February and November (because of the low temperature during the winter).
Males occupy territories and survey them from a high vantage point.
Mating take place in March and females lay 1-5 eggs in July in mainland populations. Insular populations may have smaller clutch sizes. Juveniles emerge in September.
(Ref: Jeroen Speybroeck. Wouter Beukema, Bobby Bok. Jan Van Der Velikov. 2016. Field Guide to the Amphibians and Reptiles of Britain and Europe. British Wildlife Field Guides)

http://gallotia.de/AS/Bibliografie/BIB_8973.pdf
-https://deepblue.lib.umich.edu/bitstream/handle/2027.42/110598/evo12555.pdf?sequence=1&isAllowed=y
-
http://www.bpo-natura.de/en/index.php?id=163
-https://deepblue.lib.umich.edu/bitstream/handle/2027.42/110598/evo12555.pdf?sequence=1&isAllowed=y
-http://old.biol.uoa.gr/enphysg/pdf/hus.pdf

# data[].attributes.attachment_link.title
Καταγραφή.PNG

# data[].attributes.headline
Physiology and adaptation of amphibian species

# data[].attributes.subject
Physiology and adaptation of amphibian species

# data[].attributes.body
Pelophylax cerigensis (TL: 5.5 cm) has relative rough skin and it is found in Karpathos, Greece (and possibly also in Rhodes), at elevations around 300 m (Beerli et al. 1994). Pelophylax cerigensis is only known with certainty from a single river on northern Karpathos Island. The mountainous river is near Olimbos. The required habitat is still or slow-moving waters where mating and larval development occur (Beerli et al. 2009).
An amorous, shy frog that often shines out of the water and dives in it as soon as it feels danger. It mates the spring (between March & June) and the females give birth to eggs in several hundreds of clusters in the water.
Active throughout the year if conditions are right but generaly much more active in spring and summer.
The Karpathos marsh frog adopted a rather simple diet, comprising fewer prey taxa than other ranid frogs in the area (Cogălniceanu et al., 2000; Mollov et al., 2006; Bisa et al., 2007; Bogdan et al., 2013). The dominant prey groups were insects, spiders and isopods. Frogs from the two examined populations prey on the same invertebrate taxa (11 common groups), with limited differentiation between them.
(Ref: Pafilis, P., Kapsalas, G., Lymberakis, P., Protopappas, D., Sotiropoulos, K., 2019. Diet composition of the Karpathos marsh frog (Pelophylax cerigensis): what does the most endangered frog in Europe eat? Animal Biodiversity and Conservation, 42.1: 1–8).

https://www.researchgate.net/profile/Panayiotis_Pafilis/publication/328307210_Diet_composition_of_the_Karpathos_marsh_frog_Pelophylax_cerigensis_what_does_the_most_endangered_frog_in_Europe_eat/links/5bc5c309299bf17a1c559e26/Diet-composition-of-the-Karpathos-marsh-frog-Pelophylax-cerigensis-what-does-the-most-endangered-frog-in-Europe-eat.pdf
-https://amphibiaweb.org/species/5000
-https://www.iucnredlist.org/species/58567/11787309
-https://www.researchgate.net/publication/305113519_The_frog_of_Karpathos_Pelophylax_cerigensis_and_the_Management_Agency_of_Karpathos_Saria_in_Greek_with_English_abstract
-https://www.inaturalist.org/taxa/66317-Pelophylax-cerigensis

# data[].attributes.attachment_link.title
Καταγραφή.PNG

# data[].attributes.body
(As a scientist)
Αccording to Richard Jenkins
et al., an understanding of the conservation status of Madagascar’s endemic reptile species is needed to underpinconservation planning and priority setting in this global biodiversity hotspot, and to complement existing information onthe island’s mammals, birds and amphibians. So it is very important that we do not let the species of reptiles in Madagascar extinct.
there are several solutions to avoid this.
Continued exploration of taxonomy and diversity of reptilespecies, to gain a realistic view of Madagascar’sreptile diversity and, consequently, its threat status. To assess the vulnerability of Malagasy reptiles to climatechange, and increase scientific studies On the impact of forest destruction and its degradation on reptiles. Also, the protected areas could be reviewed as there are populations of reptiles endangered in areas outside of them.
http://www.herpconbio.org/Volume_10/Issue_1/Theisinger_Ratianarivo_2015.pdf
-https://www.ajol.info/index.php/mcd/article/view/57338
-https://www.researchgate.net/profile/Derek_Schuurman/publication/235346782_Amphibians_in_the_Mist_-_Madagascar%27s_endemic_frogs_increasingly_under_threat/links/02bfe5112401cc068a000000
-https://zslpublications.onlinelibrary.wiley.com/doi/full/10.1111/j.1469-1795.2011.00459.x
-https://www.researchgate.net/publication/264897010_Extinction_Risks_and_the_Conservation_of_Madagascar%27s_Reptiles

# data[].attributes.attachment_link.title
Καταγραφή.PNG

# data[].attributes.body
Amphibians and reptiles, like much of the world's fauna, are subject to anthropogenic pressures. Habitat lose represents not the only threat of greatest concern to the herpetofauna. Pathogens comprise a seperate category of threat. Infection with the fungus Batrachchytrium dendrobatitis has caused the collapse of amphibian populations inm any different parts of the world (Jeroen Speybroeck et al., Field Guide to the Amphibians & Reptiles of Britain and Europe. 2016. British Wildlife Field Guides). Also, there are some other threats listed in the list below:
1. Habitat loss, degradation and fragmentation
2. Invasive alien species
3. Diseases
4. Overexploitation
5. Climate change

# data[].attributes.attachment_link.description
The Amphibians & Reptiles SSC directory provides the full list and contact details for all Amphibians & Reptiles SSC Groups.

# data[].attributes.body
Recent climate change has caused the distributions of many species to shift poleward in islands.
After a scientific study it was found that populations extinction rates are higher for those species that today have more northern distributions. Northern species requiring cool, mesic habitats have less available suitable habitat among islands, implicating loss of suitable habitat in their elevated extinction rates.
These extinctions occur in the context of increasing fragmentation, with islands shrinking and separating as sea levels rose. Thus, the circumstances faced by reptiles on the islands are similar to challenges for numerous species today that must cope with a changing climate while living in an increasingly human-fragmented landscape.

# data[].attributes.attachment_link.readable_text
Abstract

Recent climate change has caused the distributions of many species to shift poleward, yet few empirical studies have addressed which species will likely be vulnerable to longer-term climate changes. To investigate past consequences of climate change, we calculated the population extinction rates of 35 reptile species from 87 Greek land-bridge islands in the Mediterranean that occurred over the last 16,000 years. Population extinction rates were higher for those species that today have more northern distributions. We further found that northern species requiring cool, mesic habitats had less available suitable habitat among islands, implicating loss of suitable habitat in their elevated extinction rates. These extinctions occurred in the context of increasing fragmentation, with islands shrinking and separating as sea levels rose. Thus, the circumstances faced by reptiles on the islands are similar to challenges for numerous species today that must cope with a changing climate while living in an increasingly human-fragmented landscape. Our island-biogeographical approach to investigating historical population extinctions gives insight into the long-term patterns of species responses to climate changes.

Keywords: habitat fragmentation, island biogeography, conservation, faunal relaxation, herpetology, Aegean Sea

Introduction

There is now evidence that many species from a wide variety of taxonomic groups have experienced poleward range shifts over the last several decades of climate change ( Walther et al. 2002, Parmesan and Yohe 2003, Root et al. 2003, Parmesan 2006). Despite this fingerprint of global climate change on biological systems, inference about climate impacts on populations over the next century is limited on several fronts. First, studies documenting recent range shifts are generally restricted to species for which extensive time-series data are available ( Parmesan et al. 1999) or which have experienced rapid shifts ( Shoo et al. 2006). Therefore, documented cases of range shifts are dominated by species that are generally well-studied and vagile, such as birds and butterflies ( Thomas and Lennon 1999, Wilson et al. 2005, Chamaille-Jammes et al. 2006, but see Hickling et al. 2006, LaSorte and Thompson 2007). Second, the majority of studies on range shifts has focused on the leading (poleward) edge of the range and hence range expansion ( Hampe and Petit 2005, Thomas et al. 2006). Movement of the trailing edge of the range, however, may be more important for species conservation, because this drives the contraction of a species' range; rapid range retraction from the trailing edge coupled with a slow expansion at the leading edge could place species in jeopardy of substantial population reduction or extinction ( Huntley 1999, Franco et al. 2006). Finally, current climate change is occurring simultaneously with habitat fragmentation as humans are changing landscapes ( Warren et al. 2001, Travis 2003). Fragmentation will inhibit range shifts not only because it limits poleward range expansion ( Honnay et al. 2002), but also because it can increase population extinction rates at the trailing edge ( Koprowski et al. 2005). At present it is unknown which species in habitat fragments (e.g., national parks, relict habitats) will be most susceptible to climate change.

Although studies on the effects of recent climate change on species have understandably focused on recent changes in the distributions and abundances of species, additional insight can be gained by investigating longer-term patterns of population extinctions and species distributional changes ( Davis et al. 2000, Barnosky et al. 2003). Since the height of the last glaciation, the world's climate has changed dramatically; mean temperatures have increased and precipitation patterns have changed regionally ( Issar 2003). By examining patterns of species distributional changes during these past climatic changes, we can identify characteristics of species that may make them vulnerable to current and future climate changes.

We investigated the natural population extinctions of 35 reptile species from 87 land-bridge islands in the northeast Mediterranean Sea that separated from the mainland starting 16,000 years ago at the end of the last ice age. These islands today harbor only impoverished species communities compared to mainland areas of equal sizes, indicating that population extinctions have been common (see Results). The environmental changes faced by these 35 species were similar to those that will be experienced by large numbers of species over the next century of climate change, including both habitat change and fragmentation. In the eastern Mediterranean the early Holocene (11,000-7,000 BP) was a period of wet, mild conditions during which there were extensive forested areas ( Caner and Algan 2002, Magny et al. 2002, Mudie et al. 2002, Eastwood et al. 2007), while the period after 6,000 BP saw increasing aridity ( Eastwood et al. 2007) and reduction in forest cover over most ( Rossignol-Strick 1999b, Roberts et al. 2001, Geraga et al. 2005), but not all ( Bottema and Sarpaki 2003), of the region. Thus, the 35 study species experienced climate change and a concomitant vegetation change. Second, increasing global temperatures led to a 120-m rise in sea level that resulted in the sequential fragmentation of a continuous coastal landscape into land-bridge islands. Hence, the study species experienced a progressive fracturing of continuous populations into ever-smaller island groups, mirroring the situation faced by many species today that live in an increasingly human-fragmented landscape.

Here, we first calculate the natural extinction rates of island populations of 35 reptile species using the current presence/absence distribution of each species. Second, we analyze the relationship between these extinction rates and the distributional patterns of the species on the mainland. In particular, we ask whether species that have more northerly ranges on the mainland are more likely to have higher population extinction rates on the study islands. We would expect this to be the case if more northern species have lower tolerance for the hotter, drier conditions and habitats at the equatorial boundary of their ranges. Thus, we are using the latitudinal range of species on the mainland as a surrogate for a number of traits that might be associated with climatic preferences or tolerances of species. Finally, we ask how the amount of suitable habitat available for each species influenced extinction rates. This question addresses the importance of climatic effects on the biotic environment for the vulnerability of species to extinction.

Methods

Species selection

We analyzed data from 87 Greek islands. Present-day reptile distributions were collected from the literature ( Chondropoulos 1986, 1989) and were supplemented by field surveys by JF. The present reptile fauna of all Aegean islands consist of 49 species. Of these, we included the 35 species that occurred on at least one of our 87 study islands and for which we could obtain data on population densities and habitat specialization. The 35 species analyzed in the present study include 28 species from an earlier study ( Foufopoulos and Ives 1999b) and 7 additional species with newly collected distributional information. A listing of the distribution and ecological characteristics of all species is available in the Online Appendix.

Calculation of extinction rates

We calculated population extinction rates for each species from the current presence/absence of species on islands that were formerly connected to a larger landmass during the last ice age ( Perissoratis and Conispoliatis 2003). To infer the past distribution of reptile species among islands, we divided islands into 5 archipelagos that were either connected to the mainland (Eastern Aegean, Argosaronic, Sporades and Ionian clusters) or formed a large continuous landmass (Cyclades) during the last glacial maximum. We inferred the sequence and timing of island formation from bathymetric charts and the known rate of sea-level rise since the height of the Wisconsin-Würm glaciation ( Pirazzoli 1991, Perissoratis and Conispoliatis 2003). This resulted in 5 island cladograms whose roots represent the ancestral landmasses and tips the present-day islands (., see also ( Foufopoulos and Ives 1999a).

We then mapped the patterns of species presence/absence onto these island cladograms under the parsimony assumption that a species occurred on an ancestral island if it occurred on any of its daughter islands that formed as the island fragmented (). The assumption that current presence indicates past occurrence and not immigration is based on a wealth of evidence demonstrating that the chances of over-water colonizations - whether natural or as the result of anthropogenic activities - are very low ( Foufopoulos and Ives 1999b). First, in contrast to tropical island regions like the Caribbean or Melanesia that have warm waters and abundant vegetation forming rafts that enable dispersal ( Calsbeek and Smith 2003), in the Aegean the relatively cold waters, general absence of floating vegetation, and substantial inter-island distances make over-water travel very difficult. Second, if substantial overwater colonization had occurred, one would expect islands closer to either the mainland or other larger islands to harbor more species. However, analyses revealed that there is no such relationship; instead, the number of species on islands was negatively associated with the time at which islands were isolated during the island fragmentation process ( Foufopoulos and Ives 1999b). Third, lack of over-water colonization is reflected in some seemingly idiosyncratic distributions where species have failed to disperse across narrow but deep straits separating islands. Taxonomic analyses based on morphological and molecular criteria resulted in the description of 5 allopatric Podarcis lizard species, more than 25 subspecies of Podarcis erhardii, and at least 13 subspecies of Cyrtopodion kotschyi in the Aegean Sea alone ( Mayer and Tiedemann 1980, Beutler 1981, Mayer and Tiedemann 1981, Gruber 1986, Kasapidis et al. 2005, Poulakakis et al. 2005a, Poulakakis et al. 2005b, Lymberakis et al. 2008). The distribution patterns of these groups mirror the time of separation of islands rather than geographic proximity, showing lack of substantial gene flow across water straits. Finally, detailed population genetic analyses of the lizard P. erhardii demonstrated lack of an isolation-by-distance effect in the genetic structure of the examined populations and a progressive loss of genetic diversity over time. Both of these relationships indicate the absence of significant gene flow even between populations separated by narrow water straights ( Poulakakis et al. 2005c, Hurston et al. in review).

To calculate extinction rates, we assume that the instantaneous rate of extinction of a species is constant between the time of island formation and fragmentation (ancestral islands) or since the time of formation (extant islands). We also assume that the extinction rate depends on the size of the island. Specifically, for island i we assume the instantaneous extinction rate is r − a( − ) where r is the overall (average) extinction rate and a measures the sensitivity of the extinction rate to island area, . For extant islands is the current area of the island, and for ancestral islands is the size of the island immediately before it separated into daughter islands, giving the minimum size of the ancestral island; this assumes that the minimum island area sets the extinction risk. The average area of islands (both past and present) from the archipelagos on which the species occurred, , is subtracted from so that r gives the average extinction rate over all islands, since the average of ( − A) over all islands i is zero. With these assumptions, the probability that a species goes extinct on island i is given by the logistic regression model

where is the fate (0 = extinction, 1 = persistence) of the species on the present or ancestral island i, and is the time between when the island was formed and either when it divided or today. This model produces estimates of extinction rates, rather than extinction probabilities, because time is incorporated; preliminary analyses (not presented) showed that time was a significant predictor of population extinctions. In summary, the logistic regression given by equation 1 estimates for each species the mean extinction rate r and the parameter a that represents the sensitivity of extinction rates to island area, with larger values of a corresponding to increased susceptibility to extinction on smaller islands. illustrates the information used in the calculations.

In fitting equation 1, rather than use the logit link function normally used in logistic regression ( McCullagh and Nelder 1989), instead we used the log link function, because this easily accommodates the exponential form of equation 1. Logistic regression is known to give biased estimates of the coefficients, and we therefore modified the method of Firth (1993) that incorporates a penalty into the likelihood proportional to the determinant of the corresponding information matrix. To obtain confidence intervals for the estimates of extinction rates r and area sensitivity coefficient a, we performed bootstrapping. This not only revealed the precision of the estimates but also any residual bias that was not corrected by the modified Firth method. Extinction rates and area sensitivities for all species are presented in the Online Appendix along with bootstrapped confidence intervals.

Our approach involving the reconstruction of the sequence of island separation and species loss gives a more-refined picture of extinction rates than would be obtained from the current presence/absence patterns of species among islands alone. For example, if a species is absent from an entire archipelago despite being present on a nearby mainland, we scored these multiple absences as only a single, ancient population extinction event. This guards, at least in part, against incorrectly inflated estimates of extinction rates for species that currently have sparse distributions among islands.

Our estimation of extinction rates assumes that species were originally uniformly distributed across the ancestral landmass, since we assume that absences of species from islands indicate extinctions. The (sub)fossil record from the islands, although far from complete, supports the existence of rich pre-fragmentation species communities with several documented cases of persistence (fossils from islands where the species still survive at present) and also demonstrates extinctions (fossils on islands of species that went subsequently extinct) ( Bachmayer et al. 1975, Schneider 1975, Szyndlar 1991). Nonetheless, to assess the sensitivity of our results to the two underlying assumptions of the estimation of extinction rates - uniform ancestral distributions of species and no colonization between islands - we performed two simulation studies ( Appendix A). In the first, we violated the assumption that species were distributed uniformly throughout the ancestral landmass from which extinctions were calculated by instead supposing that the ancestral species distributions were patchy. Second, we violated the assumption that there was no immigration between islands by incorporating colonization. Although both of these model assumptions are supported by considerable evidence, these simulation studies nonetheless make our results more secure. Even large violations of both assumptions do not change our overall conclusions ( Appendix A).

Extinction rates and species characteristics

We hypothesized that if the current latitudinal distribution of species on the mainland reflect their sensitivity to climatic conditions, then there should be a relationship between the latitudinal range of the 35 species and their extinction rates on Aegean islands. To determine the latitudinal ranges of species, we identified the northern and southern-most points of each species' present range from published distributional information ( Böhme 1981) and museum records. The northern and southern-most range points were averaged to give the latitudinal midpoint ( LM) of the species' range. Although there are multiple ways of measuring the latitudinal range of species, our approach is simple and incorporates the full range of latitudes at which species are found. We explore alternative measures of species ranges, specifically the northern and southern range limits, and the latitudinal span (= northern − southern limit) in Appendix B.

We regressed extinction rates on the latitudinal midpoints of ranges. Because the species are phylogenetically related, such regressions run the risk of type I errors, identifying a statistically significant pattern when in fact there is none (e.g., Martins and Garland 1991). Therefore, we analyzed the data with a regression method that incorporates the species phylogeny. Standard phylogenetic comparative methods assume that the residual variation from the regression can be described by a Brownian motion evolutionary process; this is the assumption of both independent contrasts ( Felsenstein 1985, Garland et al. 1992) and Generalized Least Squares approaches ( Martins and Hansen 1997, Garland and Ives 2000). Rather than make this strict assumption, instead we use the program "Regressionv2.m" from Lavin et al. (2008) in which the strength of phylogenetic signal exhibited by the residuals of the analysis is estimated simultaneously with the regression coefficients. The strength of the phylogenetic signal is given by a parameter d that varies from 0 when residuals are independent, to 1 when the phylogenetic correlations among residuals are consistent with Brownian motion evolution. Thus, the parameter d statistically adjudicates between non-phylogenetic regression and regression under a Brownian motion model of evolution. To correct for heteroscedasticity in the residuals, we square-root transformed the extinction rates prior to analyses.

In addition to the latitudinal midpoint of species ranges, we considered two other species traits that were previously shown to be related to extinction rates of reptile populations on Aegean islands, maximum population density and the breadth of habitats used by the species ( Foufopoulos and Ives 1999b). To measure population density, species were categorized into four groups according to the maximum density (number per unit area) they achieve on the study islands. We use these maximal estimates as species-specific traits, rather than as an attribute of any one population. Because we are using these estimates as traits that affect extinction rates, we are assuming that maximum densities have not changed appreciably over the course of the Holocene. We used 4 broad categories of maximum densities because uncertainty in these species traits prevents any finer-scale designations. To obtain a metric of the breadth of habitats used by species, we identified 10 habitat categories that summarize the types of habitats on the islands ( Watson 1964). Our metric of breadth of habitat use is the fraction of these 10 categories used by each species. We performed multiple regression as above, using Regressionv2.m to account for possible phylogenetic signal in the residuals ( Lavin et al. 2008).

Extinction rates and suitable habitats

We also performed analyses to determine whether availability of suitable habitat could explain differences in extinction rates. To obtain measures of the amount of suitable habitat remaining on each island for each species, we used ArcGIS to analyze a CORINE 1:100,000 land cover database ( EEA 2000) that included 71 of the 87 islands in our study. The 44 land-use/land-cover classes in CORINE were collapsed into the 10 habitat categories used previously to estimate breadth of habitat use. Various reptile groups differ broadly in the minimum area of suitable habitat they require to persist. To scale the amount of suitable habitat by the area requirements of a species, we determined the minimum area of suitable habitat on any island that was supporting a species in each order (turtles: 6.25 ha) or sub-order (snakes: 2.88 ha, lizards: 0.5 ha). The total area of suitable habitat for each species among all islands was then divided by the minimum area each species required according to their order or suborder. Finally, we scaled these values to range between 0 and 1 by dividing this value for each species by the maximum value obtained among all species, thereby giving a metric of the availability of suitable habitat for each species, SH. Because the resulting values were not normally distributed, they were arcsine-transformed SH for subsequent analyses. Note that SH is a measure of island characteristics but scaled according to the habitat requirements for the different groups of species. In contrast, the maximum population density of species and the breadth of habitat use are treated as species traits in the analyses of the last section ( Extinction rates and species characteristics).

To determine whether the loss of suitable habitat could explain the higher extinction rates of more-northern species, we regressed SH against the latitudinal midpoint of species ranges, LM, and we then regressed the extinction rates r against SH. Thus, we addressed whether more northern species have less suitable habitat currently available on the study islands, and then whether a lack of availability of suitable habitat could explain higher extinction rates. These analyses were performed using Regressionv2.m ( Lavin et al. 2008) to incorporate the possibility of phylogenetic signal.

Potentially confounding factors

We analyzed data to test two general assumptions we used to calculate extinction rates: (i) humans have played an at most minor role in the extinction or redistribution of species among islands, and (ii) species were broadly distributed among ancestral islands, so that the current absence of a species from an island implies population extinction.

(i) Human impact on the distribution of species

It is possible that humans have played a role in the colonization and extinction of species from islands, and therefore that the present distributions reflect human activity rather than natural long-term patterns of climate change and population extinctions. To assess the effect of human habitation, we regressed the number of species per island on whether or not islands were inhabited by humans, including island area in the regression because island area is known to affect the number of species per island. Because we are particularly interested in the extinction rates of northern vs. southern reptile species, we also divided the species into a northern and a southern group, with the northern group having latitudinal midpoints above the mid-latitude of the islands, 37.6°N ( N = 22 species), and the southern group below ( N = 13). If humans caused the positive relationship between extinction rates and species latitudinal midpoints, then we would expect these regressions to show a larger decrease in the occurrence of species caused by humans for northern species than for southern species. Note that for these analyses the dependent variable is the number of species on islands. Therefore, to determine whether the effect of human habitation on the distribution of species depends on the latitudinal range of species, we have to categorize species as northern or southern and perform separate analyses for both groups.

(ii) Homogeneous pre-fragmentation distribution of species

Our calculations of extinction rates ( equation 1) assume that species were distributed uniformly across the ancestral landmasses of the archipelagos on which they currently occur. If species were initially patchily distributed, then the sequential fragmentation of islands during the Holocene might not have led to extinctions. Instead, species that do not occur on an island today might not have existed in the same location prior to fragmentation. This concern is addressed in Appendix A using a simulation study. Here we address this issue empirically.

We constructed a species-area curve for the islands and compared this to a species-area curve for the mainland areas adjacent to the islands with species area data from Chondropoulos (1986, 1989) and field data from one of the authors (JF). The areas on the mainland occupy low-lying regions near the coast, and are topographically and climatically similar to the islands. If this analysis shows that islands have fewer species than mainland areas of the same size, then this would imply a strong role of extinctions in determining the pattern of species occurrences on islands ( MacArthur and Wilson 1967). If instead the distribution of species on current islands were determined solely by their former, patchy distribution across the ancestral landmass, then there is no reason to expect that the number of species on islands would be lower than comparably sized areas on the mainland.

Results

Extinction rates, r, which give the instantaneous probability that a population goes extinct from an island of average size, ranged from 0.0095 to 0.47 per 1,000 years (Online Appendix). For the 35 species, there were only 3 whose southern boundary on the mainland was further north than the southern-most island on which it was scored as being absent in the calculation of extinction rates (Online Appendix). Therefore, the extinction rates do not simply reflect the northward movement of the southern range limits of species. The area sensitivity parameter a ranged from −0.034 to 0.070 per 1,000 years per km ². All but two species had positive estimates of a, indicating that extinction rates increased on smaller islands (Online Appendix); of the two species with negative estimates, one was non-significant and the other weakly significant (0.05 > p > 0.02).

Extinction rates and species characteristics

Extinction rates r were positively correlated with the latitudinal midpoint of species ranges on the mainland, LM (, ). In the regression the estimate of the strength of phylogenetic signal (parameter d) was zero. This indicates that closely related species were not more likely to have similar extinction rates after the effect of LM was removed.

Table 1

A. Phylogenetic regression of square-root transformed extinction rates r (per
1,000 years) on the latitudinal midpoint of species ranges. The estimate of phylogenetic
signal in the residuals is d = 0, with N = 35 and R ² = 0.27. Analyses were performed
using Regressionv2.m from Lavin et al. (2008).

B. Phylogenetic regression of square-root transformed extinction rates r (per 1,000
years) on maximal species densities, the breadth habitats used by species, and LM. The
estimate of phylogenetic signal in the residuals is d = 0, with N = 35 and R ² = 0.62. Analyses
were performed using Regressionv2.m from Lavin et al. (2008).

We also performed multiple regression of extinction rates r against LM and two additional variables known to affect extinction rates, the maximum density achieved by species and their breadth of habitat use (). Both of these additional variables were statistically significant, as was LM, and there was no evidence of phylogenetic signal in the residual variation ( d = 0). The regression coefficient for LM, 0.013, was only slightly lower than in the regression with LM alone (0.016, ), indicating that the additional two variables do not explain the observed relationship between extinction rates and LM.

Extinction rates and suitable habitat

For each species, we calculated a metric SH for the relative amount of suitable habitat available across the study islands. SH was negatively related to the latitudinal midpoint of species ranges, LM, indicating that more northern species have less suitable habitat available (, ). Furthermore, extinction rates r were negatively related to SH, indicating that species with more available habitat had lower extinction rates (, ). This provides one potential explanation for the relationship between extinction rates r and LM (): northern species had less available suitable habitat, and as a consequence had higher extinction rates.

Table 2

A. Phylogenetic regression of the arcsine-transformed metric of suitable habitat
available on islands, SH, and the latitudinal midpoint of species ranges. The estimate of
phylogenetic signal in the residuals is d = 0.13, with N = 35 and R ² = 0.32. At this small
value of d, the estimates for the coefficients and t-values are within 1% of those obtained
under standard (non-phylogenetic) regression. Analyses were performed using
Regressionv2.m from Lavin et al. (2008).

B. Phylogenetic regression of square-root transformed extinction rates r (per 1,000
years) on arcsine-transformed SH. The estimate of phylogenetic signal in the residuals is d =
0, with N = 35 and R ² = 0.30. Analyses were performed using Regressionv2.m from Lavin et al. (2008).

C. Phylogenetic regression of square-root transformed extinction rates r (per 1,000
years) on arcsine-transformed SH and LM. The estimate of phylogenetic signal in the
residuals is d = 0, with N = 35 and R ² = 0.36. Analyses were performed using
Regressionv2.m from Lavin et al. (2008).

A multiple regression of r against both LM and SH provides support for this explanation. Even with the inclusion of LM, the coefficient for SH on extinction rates remains statistically significant. However, the coefficient for LM is reduced from 0.016 () to 0.0089 (), and statistical significance of this relationship is lost. This indicates that, while both LM and SH are related to extinction rates, at least part of the explanation for northern species having higher extinction rates is that they have less available habitat.

Human impact on the distribution of species

The number of reptile species on an island was not statistically related to the presence of humans in a regression analysis that accounted for the effects of island area (humans, t = −0.175, p = 0.861; area, t = 10.427, p < 0.001; total model R ² = 0.809, N = 87). These results suggest it is unlikely that the presence of humans has impacted extinction or colonization by reptile species. We performed a similar analysis on the subsets of 22 northern ( LM above the mid-latitude of the islands) and 13 southern species (northern species: t = −0.317, p = 0.752; southern species, t = −0.106, p = 0.916). Because neither the subsets of northern nor southern species showed an effect of human presence, it is unlikely that human presence was responsible for the higher extinction rates of northern species ().

Species-area relationships on islands and mainland

Species-area curves for the islands and the mainland showed lower numbers of species on islands than on comparably sized areas on the mainland (). This suggests that, similar to other land-bridge islands or habitat fragments ( Wilcox 1978, Newmark 1987, 1991, Burkey 1995), the reptile communities in the Aegean underwent relaxation following isolation from the mainland. As the steepness of the island slope indicates, the population extinctions underlying this relaxation process were most pronounced on smaller islands. This is consistent with our calculations of area sensitivity, a, in the extinction rates, for which the preponderance of negative species-specific values indicates higher extinction rates on smaller islands. The apparently strong effect of extinction in the island species-area curve relative to that from the mainland indicates that the pattern of species occurrences on islands that we used to calculate extinction rates cannot be attributed to an initially patchy distribution of species across the ancestral landmass.

Discussion

We found that island populations of more-northern reptile species had higher extinction rates r (, ). This is consistent with the expectation that species from northern climates were impacted more severely by climate change on the study islands during the Holocene ( Geraga et al. 2005, Eastwood et al. 2007, Fuchs 2007). While investigating possible explanations for this pattern, we also accounted for two species characteristics known to affect extinction rates: the maximum population density of a species ( Pimm et al. 1988, Purvis et al. 2000) and the breadth of habitat types used by the species ( Rossignol-Strick 1999a, Owens and Bennett 2000). Our goal was to determine whether these two species characteristics underlie the observed relationship between extinction rates and the midpoint of species' latitudinal range. For example, more northern species might be species that naturally occur at lower population densities, and the resulting smaller populations might be more prone to extinction. However, including these two species traits did not reduce the strength of the relationship between latitudinal range midpoints and mean extinction rates ().

Higher population extinction rates of northern species could be due to direct effects of climate (e.g., temperature) or indirect effects through changes in the biotic environment such as shifts in plant communities. Although it is not possible to completely disentangle direct and indirect effects, circumstantial evidence suggests that higher extinction rates of northern species resulted at least in part from the disproportionate loss of more cool, mesic habitats (e.g., deciduous and coniferous forest). The amount of suitable habitat for each species, SH, was negatively related to the latitudinal midpoint of species ranges, LM. Furthermore, extinction rates r were negatively related to SH. Reduced area of suitable habitat could lead to extinction by making some islands unsuitable, or by restricting the suitable area of islands that thereby reduces population sizes and makes them more prone to extinction. Whereas conditions were cool and moist at the beginning of the Holocene when most of the islands became separated, the region has since experienced increasing aridity ( Caner and Algan 2002, Magny et al. 2002, Mudie et al. 2002, Eastwood et al. 2007) and a general decrease in forest cover ( Roberts et al. 2001, Rossignol-Strick, 1999, but see Bottema and Sarpaki 2003, Geraga et al. 2005). Consequently, the paleoclimatic record, though incomplete, suggests that in the northeast Mediterranean Basin conditions have become increasingly unsuitable for taxa with cool or mesic habitat preferences.

Human land use could also influence reptile persistence on the islands, and humans caused significant changes in vegetation cover after 6,000 BP ( Fuchs 2007) on 40 (46%) of the study islands that have been inhabited. Nevertheless, in field observations (JF, unpublished data) we found that human presence had surprisingly little effect on reptile communities, with all resident species widely distributed throughout anthropogenically modified habitats and no clear differences between areas with and without human habitation. Indeed, in our statistical analyses, after correcting for island area, the presence of humans had no effect on the number of reptile species on islands, either for all species combined or for subsets of northern and southern species. Thus, the absence of a difference in human impacts on the northern and southern subsets of species implies that humans are unlikely to have been responsible for the higher extinction rates of northern species from islands. Our data indicate that while human activities have had effects on the vegetation cover of the islands in the region ( Grove and Rackham 2003), this did not translate into observable changes in the local reptile communities, possibly because the gradual expansion of non-intensive agricultural practices allowed local taxa to adapt to an anthropogenically modified landscape.

Our calculations of population extinction rates presuppose that all species were distributed throughout the proto-archipelagos of the Aegean Sea; we assumed that absence from an island indicated extinction either from that island or from its island progenitors. If instead species were patchily distributed across the pre-fragmentation landscape, the current patchwork distribution among islands could reflect the former patchy distributions rather than post-fragmentation population extinctions. However, simulation studies of our procedure for calculating extinction rates showed that our conclusions are robust to even highly patchy initial distribution of species across the ancestral landmass ( Appendix A). Furthermore, comparing the species-area curve of the islands to that of the mainland provides empirical evidence that supports the simulation studies. For a given area, islands have fewer species than the mainland; for example, small islands had 1-3 species, whereas areas of the same size, elevation, and habitat heterogeneity had 6-10 species on the mainland (). Furthermore, the slope of the species-area curve for the islands was steeper than for the mainland, which is consistent with the well known pattern of high extinction rates among small populations ( MacArthur and Wilson 1967, Groom et al. 2006). If the present distribution of species among islands were caused solely by patchy ancestral distributions, we would expect no difference between island and mainland species-area curves. Indirect evidence for the widespread distribution of reptiles across the pre-fragmentation landscapes also comes from the palynological record. It indicates that during the earliest Holocene while conditions were cool and relatively dry, climatic conditions were similar to what is typical of present-day northern latitudes were these species are widespread today ( Digerfeldt et al. 2007). Therefore, in the past northern species likely had broad geographical distributions within this island region. Finally, all of the northern species currently occur on islands spread throughout the broader geographic area ( Valakos 2008), often close to sea level, again pointing towards a widespread distribution on the pre-fragmentation landscapes.

In summary, our results suggest that species from cooler, mesic regions were more susceptible to extinction from our study islands, and their elevated susceptibility was caused at least in part by the loss of suitable habitats. This loss was likely driven in large part by climate change that resulted in higher temperatures and increasing aridity during the Holocene. These patterns of extinction occurred within the context of a landscape that was becoming increasingly fragmented during the Holocene as sea levels rose. Not only are islands today depauperate compared to similar areas on the mainland (), but 33 of 35 species had higher extinction rates on smaller islands (area sensitivity parameter a > 0, OA). Thus, fragmentation that led to smaller island areas exacerbated extinctions.

What lessons do these results hold for current climate change? Because the southern range limits of almost all our study species on the mainland lay south of the southern-most island from which they went extinct, the extinctions from islands do not simply represent a retraction of the southern range limit that matches the nearby mainland. Instead, island extinctions occurred north of the present southern range limit on the mainland. Thus, range retraction in the fragmented, island landscape is a sporadic, non-uniform process, with extinctions occurring especially on small islands as a mosaic within the equatorial side of a species range. This represents an interaction of the climate change and fragmentation processes, with more extensive fragmentation (smaller islands) leading to higher extinction rates that are further exacerbated for species that are sensitive to climatic changes. From this example in the Aegean, we would expect to see the effects of climate change on species in a fragmented landscape as inflated extinction rates on the equatorial side of species' ranges before there is range retraction. These extinctions could give a warning sign of full-blown range contraction.

Extinctions of reptiles on islands in the Aegean over the last 10,000 years share many similarities with anthropogenic climate change and habitat fragmentation as it occurs today. Human activities tend to first destroy flat, lowland habitats that are amenable to agriculture and development. As a result, many natural habitats survive today only as small habitat islands confined to mountainous regions that are surrounded by a 'sea' of low domesticated land. Such upland areas tend to be topographically diverse, as are the islands in the Aegean. Similarly, the matrix surrounding habitat islands in many human-dominated landscapes is hostile to slowly dispersing terrestrial organisms like reptiles, presenting obstacles that might be almost as impenetrable as the cold waters of the Aegean Sea. This fragmentation presents a threat in its own right, and also a major challenge for species whose distributions are shifting poleward as the climate warms.

The scenario of population extinctions suggested by our analyses allows integration of several independently observed patterns. As the climate warms, species ranges are expected to contract from their equatorial boundary ( Huntley 1999, Franco et al. 2006). Loss of suitable habitat is the primary cause of extinction risk for species throughout the world (Groombridge 1992) and for Mediterranean herptiles in particular (Cox et al. 2006). And fragmentation is well-known to increase population extinction rates, as fragmentation creates small populations and isolates them from the possibility of re-colonization ( Newmark 1987). Together these components represent a complex of processes that will likely play out over centuries, rather than decades. While there have been numerous documented cases of rapid, contemporary changes in species ranges ( Walther et al. 2002, Parmesan and Yohe 2003, Root et al. 2003, Parmesan 2006), these may represent only that subset of species susceptible to rapid and direct effects of climate change. In contrast, the extinction process of reptiles from Aegean islands has progressed through millennia, and the indirect effect of climate change through habitat loss appears to play a major role. Our results suggest that in the coming centuries, fragmentation will exacerbate the effects of a warming climate leading to differential extinctions of those species that are habitat specialists and/or that naturally occur at low densities. Unfortunately, documenting slow or delayed effects of contemporary climate change on species ranges and population extinction will be difficult. Therefore, patterns of past population extinctions, while not mimicking precisely the patterns we expect over the coming centuries, may nonetheless give insights into the types of processes that will drive changes in species populations in the future.

Acknowledgments

Funding was provided through the University of Wisconsin Department of Zoology, University of Michigan, Princeton Environmental Institute, the Cleveland Dodge Foundation, and National Science Foundation grant DEB-0816613 (ARI) and EF-0914866 (AMK) and National Institutes of Health grant 1R01AI090159-01 (A.M.K.).

Appendix A

Robustness of extinction rate estimates

In our estimation of extinction rates, we assume that species are initially uniformly distributed within the 5 landmasses that eventually divided into the 5 archipelagos we used in our study, and that colonizations did not occur. There is good circumstantial evidence for these two assumptions. Nonetheless, here we investigate the consequences of violations in these two assumptions.

Initial uniform distribution among landmasses

To estimate extinction rates, we assumed that absence of a species indicates extinction either from that island or from an ancestral island. If a species were not uniformly distributed within the initial landmass before island separations started, then we would incorrectly score the absence of that species from islands that had never previously occupied that location as an extinction. This is illustrated in that shows the initial distribution of two species, sp 1 and sp 2, and island cladograms with their current presence/absence. Sp 1 (blue) occupied the land area encompassing islands b, c, and d, but not a. The absence of sp 1 on island a today does not represent an extinction. However, absence from island c does represent an extinction. Therefore, for calculating extinction rates, only information from islands b, c, and d (shaded) should be used.

Sp 2 was initially confined to the area that later formed islands c and d. Therefore, absence from islands a and b do not represent extinctions, and only information from islands c and d should be used to calculate extinction rates. Note that in this case the method we use to calculate extinction rates would count the absence of sp 2 from islands a and b as arising from only a single extinction, on the ancestral island to islands a and b, that is, proto-island a-b. Therefore, even though there are two absences, they are scored as a single extinction. This property of our estimation procedure should reduce its sensitivity to initially non-uniform distributions of species, because recently separated islands are likely to be in close spatial proximity and therefore initial absence of species from that region will be scored as only a single extinction.

Initially non-uniform species distributions will lead to the mis-scoring of absences as extinctions. To simulate the consequences of these pseudo-extinctions, for each species we computed extinction rates after randomly removing extinction events with probability 0.5. For example, for sp 1 we calculate extinction rates assuming an extinction from island c but persistence on islands b and d, and proto-island b-c. For sp 2, we assume an extinction from island d and not from island c. This random procedure does not explicitly account for the biologically expected pattern that geographically neighboring islands are likely to have had the same initial presence/absence of species; for example, islands a and b are nearby and therefore are more likely to be within the initial range of a species (although this is not the case for either sp 1 or sp 2 in our hypothetical example). Nonetheless, because the extinction events can occur on ancestral islands as well as extant islands, removal of these ancestral extinction events is equivalent to the absence of a species from an entire region of the initial landmass. Another simplifying assumption is that we do not consider the size of an island (either extant or ancestral) in randomly removing extinction events, even though larger islands are more likely initially to have had a species. Nonetheless, our simulation to examine violation of the uniform distribution assumption is severe, removing on average 50% of the extinction events from the data set. Therefore, our procedure provides a reasonable test of the robustness of our extinction calculations.

For each species we simulated 1000 data sets by randomly removing extinction events with probability 0.5. Comparing the simulated extinction rate estimates with the base estimates ( Eq. 1) obtained assuming a uniform initial distribution of species () shows that the random removal of extinction events leads to lower estimates of extinction rates. Nonetheless, there is a strong correlation between the extinction rates simulating non-uniform initial distributions and our base estimates of extinction rates. Furthermore, replacing the base estimates of extinction rates with the simulated extinction rate estimates shows a statistically significant effect of the latitudinal midpoint of species ranges ( LM) on extinction rates ( Table A1A), only slightly weaker than obtained with the base extinction estimates (). The graph of extinction rates vs. LM remains similar (compare with in the main text).

Colonization

To estimate extinction rates, we assume that the presence of a species on an island implies that the species has persisted since the landmass started to divide. If species could colonize islands, then presence on islands could be due to colonization rather than long-term persistence. This is illustrated in that shows the colonization of islands by two hypothetical species and island cladograms with their current presence/absence. Sp 1 (blue) colonized island c and sp 2 colonized island b. For sp 1, this obscured an extinction event on island c; even though sp 1 went extinct after the separation of islands c and proto-island a-b, it is currently present on island c. This shows that colonization can hide extinction events.

The effects of colonization on the inferred pattern of extinctions can be more complicated. Because sp 2 colonized island b, it covered the extinction on proto-island a-b. In addition, applying the assumptions that we used to estimate extinction rates, colonization also made it appear that sp 2 persisted on island b. Furthermore, because sp 2 is present on island b and not island a, our method would infer that an extinction had occurred on island a. Thus, there are three changes that must be simulated to estimate extinction rates when colonization is possible: (i) persistence events must be converted to extinction events (e.g., the presence of sp 1 on island c), (ii) persistence events must be removed from the data set (e.g., the persistence of sp 2 on island b, because extinction actually occurred on proto-island a-b), and (iii) extinction events must be removed from the data set (e.g., the extinction of sp 2 from island a, because extinction actually occurred on proto-island a-b). These patterns become more complex when considering colonizations that take place on ancestral proto-islands.

As a simple procedure to mimic the effects of colonization, we simulated data sets by randomly converting persistence events into extinction events with probability 0.25, and randomly removing events (either persistence or extinction events) with probability 0.5. This procedure does not correspond to a specific set of detailed assumptions about how the colonization process occurs, such as whether colonization depends on the distance to nearby islands or island area. Nonetheless, by making the magnitude of the change large (with the assumption that 25% of extinctions are "recovered" through colonization), it provides a reasonable test for the robustness of our conclusions to our base assumption of no colonization.

The estimates of extinction rates were higher when we simulated the effects of colonization than for the base assumptions (), although the estimates were highly correlated. As a result, the relationship between extinction rates and the latitudinal midpoint of species ranges ( LM) remained the same ( Table A1B, ).

These two simulation studies of the two main assumptions under which we estimated extinction rates (the uniform initial distribution of species and no colonization) demonstrate that the relative extinction rates among the different species are broadly insensitive to even large violations of the assumptions (), and the conclusion that more northerly species have higher extinction rates remains intact ( Table A1, ).

Appendix B

Alternative measures of species ranges

In the analyses we use the midpoint of species latitudinal ranges to measure the "northerliness" of species. Here we consider three other measures of species ranges: latitudinal span (northern range limit - southern range limit), southern range boundary, and northern range boundary ( Table B1). Extinction rates were negatively and marginally significantly related to latitudinal span, positively related to the southern range boundary, and positively but non-significantly related to the northern boundary. These tests are not independent, because the latitudinal midpoints and span depend on the southern and northern boundaries. Nonetheless, they are instructive in showing the overall pattern of extinction rates and the importance of species' southern boundaries.

Footnotes

Expanded online edition: Online Appendix: Distributions, ecological data and extinction rates of species in this study.

References

Beutler A. Cyrtodactylus kotschyi (Steindachner, 1870)-Aegaeischer Bogenfingergecko. In: Boehme W, editor. Handbuch der Reptilien und Amphibien Europas. Wiesbaden, Germany: AULA Verlag; 1981. pp. 53-74.

Böhme W, editor. Handbuch der Amphibien und Reptilien Europas. Vol. 1-6. Wiesbaden, Germany: Akademische Verlagsgesellschaft; 1981.

EEA. CORINE land cover database (Version 6/1999) European Environmental Agency NATLAN; 2000.

Foufopoulos J, Ives AR. Reptile distributions and island dendrograms for the islands of the Aegean and Ionion Seas. 1999a

Groom MJ, Meffe GK, Carroll CR. Principles of Conservation Biology. 3d Ed. Sunderland, MA, USA: Sinauer Associates; 2006.

Grove AT, Rackham O. An Ecological History. New Haven, USA: Yale University Press; 2003. The Nature of Mediterranean Europe.

Gruber U. Podarcis erhardii (Bedriaga, 1876)-Aegaeische Mauereidechse. In: Boehme W, editor. Handbuch der Reptilien und Amphibien Europas. Wiesbaden, Germany: AULA Verlag; 1986. pp. 25-49.

Huntley B. The dynamic response of plants to environmental change and the resulting risks of extinction. In: Mace GM, Balmford A, Ginsberg JR, editors. Conservation in a changing world. Cambridge: Cambridge University Press; 1999. pp. 69-85.

Hurston H, Foufopoulos J, Voith L, Bonanno J, Pafilis P, Valakos E, Anthony N. Effects of fragmentation on genetic diversity in island populations of the Aegean wall lizard Podarcis erhardii (Lacertidae, Reptilia) Molecular Ecology. in review. [PubMed]

Issar AS. Climate changes during Holocene and their impact on hydrological systems. Cambridge, U.K.: Cambridge University Press; 2003.

MacArthur RH, Wilson EO. The Theory of Island Biogeography. Princeton: Princeton University Press; 1967.

McCullagh P, Nelder JA. Generalized linear models. 2 edition. London: Chapman and Hall; 1989.

Owens IPF, Bennett PM. Ecological basis of extinction risk in birds: Habitat loss versus human persecution and introduced predators. Proceedings of the National Academy of Sciences of the United States of America. 2000;97:12144-12148. [PMC free article] [PubMed]

Pirazzoli PA. World atlas of Holocene sea-level changes. Amsterdam: Elsevier Science; 1991.

Purvis A, Gittleman JL, Cowlishaw G, Mace GM. Predicting extinction risk in declining species. Proceedings of the Royal Society of London Series B-Biological Sciences. 2000;267:1947-1952. [PMC free article] [PubMed]

Travis JMJ. Climate change and habitat destruction: a deadly anthropogenic cocktail. Proceedings of the Royal Society of London Series B-Biological Sciences. 2003;270:467-473. [PMC free article] [PubMed]

Valakos E, Pafilis P, Maragou M, Sotiropoulos K, J F. The amphibians and reptiles of Greece. Frankfurt, Germany: Chimaera Verlag; 2008.

Valakos PE, Maragou M, Sotiropoulos K, Foufopoulos J. The amphibians and reptiles of Greece. Frankfurt, Germany: Chimaera Verlag; 2008.

Watson GE. PhD. New Haven, CT, USA: Yale University; 1964. Ecology and evolution of passerine birds in the islands of the Aegean Sea.

# data[].attributes.attachment_link.title
Climate Change and Elevated Extinction Rates of Reptiles from Mediterranean Islands

# data[].attributes.attachment_link.description
Recent climate change has caused the distributions of many species to shift poleward, yet few empirical studies have addressed which species will likely be vulnerable to longer-term climate changes. To investigate past consequences of climate change, we calculated the population extinction rates of 35 reptile species from 87 Greek land-bridge islands in the Mediterranean that occurred over the last 16,000 years.

# data[].attributes.body
The number of conservation efforts to save amphibian and reptile species from local extinction has increased singnificantly during the last decades in many European countries, so in Greece, with the protection of sea turtle nesting sites in Greece,

http://dro.deakin.edu.au/eserv/DU:30058241/hays-estimatingthe-2002.pdf

or the saving of Ithaca's green toad, with the creation of lakes on the island.

http://www.easystarter.com/projects/swste-ton-prasinofryno-ths-ithakhs

Even moρε extreme methods are effective As shown in the paper below, which shows a model to reduce the Lithobates catesbeianus.
Μodeling suggests that culling of metamorphs in fall is the most effective method of decreasing  Lithobates catesbeianus population growth rate.

http://www.humanities.uct.ac.za/sites/default/files/image_tool/images/75/files/res19.pdf

Conservation efforts are, however, not limited to official organizations. Instead of going out to look for amphibian and reptiles in the field, one can also encourage them to colonise gardens or other semi-natural places by making improvements that will attract them. Creating basking spots compsed of rocks, pecies of wood in a sunny spot may attract varius reptile species native to the area. Creating a small ecosystem in your own garden is an easy way to support local herpetophauna populations and to observe and admire them at close range (Jeroen Speybroeck et al., Field Guide to the Amphibians & Reptiles of Britain and Europe. 2016. British Wildlife Field Guides).

# data[].attributes.attachment_link.title
Καταγραφή.PNG

# data[].attributes.headline
Business plan of a local voluntary organisation which aims to protect the local herpetofauna.

# data[].attributes.subject
Business plan of a local voluntary organisation which aims to protect the local herpetofauna.

# data[].attributes.body
Amphibians and reptiles are among the most misunderstood groups in the animal kingdom. It is logical to this happened, because most of the people do not perceive reptiles as well as mammals. Βoth reptiles and amphibians are species that are not as expressive as mammals.
Τo be able to say that the herpetofauna of a local area will be protected, we will have to gradually focus on changing this mentality.
Information programs should first be set up in schools and later in the entire local population, as far as familiarizing people with these species. Reporting on what the reptiles and amphibians offer to ecosystems are very important and people should be known.
So, this organization can share various leaflets, shows documentaries, take apart an exhibition of live reptiles and amphibians and sharing T-shirts with logos.
A voluntary group can also carry out cleansing on habitats which hosting the herpetofauna of local area. Garbage collection, removing dry shrubs and friends to prevent fire are some tasks that can be carried out.
Also, can providing training on how to identify some of the species found on site, usually at the start of the survey season (spring). The data that will be produced is vital for local records and on a national scale. It all feds into national recording schemes to build a wider picture of what's happening to the herpetofauna of the local area.
Last but not least, they can become voluntary work, thus upgrading the habitat, such as the construction of a lake (for amphibians), or the gathering of stones and wood to encourage the colonization of reptiles in the local area.

https://www2.oxfordshire.gov.uk/cms/sites/default/files/folders/documents/environmentandplanning/countryside/protectedspecies/reptilesguidelinesdevelopers.pdf

# data[].attributes.attachment_link.provider_name
The Amphibian and Reptile Conservation Trust

# data[].attributes.attachment_link.readable_text
If you would like to find out more about volunteering with ARC check out our Frequently Asked Questions , volunteer testimonials or contact our volunteer co-ordinators .

What sort of activities can volunteers help with?

There are a number of tasks that volunteers for Amphibian and Reptile Conservation can get involved in...

A huge part of our volunteer work focuses on the management of our reserves, many of which are heathlands. Volunteers can get involved with the conservation and restoration of this special and endangered habitat on sites in both Dorset and Surrey. The regular practical tasks vary from week to week depending on the season and needs of a particular reserve. The majority of winter work is the removal of scrub and keeping the habitat in good condition, usually accompanied by a warming fire and lots of tea and biscuits - great on a chilly day! There are also a variety of summer activities including guided walks, construction, habitat work and surveying.

Take a look at this video filmed and edited by ARC volunteer Dave Meyers showing what it's like to do practical conservation volunteering with Amphibian and Reptile Conservation.

In the spring and summer months there are some opportunities to carry out surveying and monitoring of reptiles and amphibians on our reserves. Monitoring of species is a great way to get to know the sites you may have done practical work on, and find out if the management is actually working. We can provide training on how to identify some of the species found on site, usually at the start of the survey season. The data that gets produced is vital for local records and on a national scale. It all gets fed into national recording schemes to build up a wider picture of what is happening to our herpetofauna.

If you are not local to a site where surveying and monitoring is being carried out, take a look at our National Amphibian and Reptile Recording Scheme (NARRS), our targeted Natterjack toad survey and other ways to record your sightings .

Awareness raising and educating the public about amphibians, reptiles and their disappearing habitats is a vital part of our work.

Each year ARC attends a variety of public events; ranging from local fetes to county shows. Volunteering with us at events is a great way to gain experience of working with the public and a fun day out!

Amphibian and Reptile Groups of the UK (ARG UK)

As well as there being plenty of opportunity to volunteer with us, Amphibian and Reptile Conservation also works closely with the Amphibian and Reptile Groups of the UK (ARG UK) - a network of local volunteer groups, usually organised on a county basis.

The structure and activity of groups varies according to local interests and priorities but all ARGs are concerned with the conservation of native amphibians and reptiles in their local area. The primary activities of ARGs are survey and recording, but groups also carry out habitat management work, educational projects, site visits and toad patrols. ARGs are the sources of local expertise for amphibian and reptile conservation, but new members are always welcome and most groups run training events.

The ARC Weald Volunteers at Thursley NNR

Latest

Related

Most read

Tag cloud

Let us know which amphibians and reptiles you've seen Read more

Sign up and receive our e-Newsletter Read more

# data[].attributes.attachment_link.description
If you would like to find out more about volunteering with ARC check out our Frequently Asked Questions , volunteer testimonials or contact our volunteer co-ordinators . There are a number of tasks that volunteers for Amphibian and Reptile Conservation can get involved in...

# data[].attributes.body
Greece lies in the south of the Balkan peninsula and is part of the European East Mediterranean Region.
Its territory expands in approximately 132.000 km2, its coast line is 15.000 km long, and has a great variety of geomorphological formations and rocks. In general, the climate can be regarded as Mediterranean, with mild and wet winters, and hot and dry summers.

# data[].attributes.attachment_link.description
Υπουργείο Περιβάλλοντος και Ενέργειας

# data[].attributes.body
https://www.researchgate.net/figure/Map-of-forest-vegetation-zones-extant-or-potential-of-Greece_fig1_324034455

# data[].attributes.attachment_link.title
Καταγραφή.PNG

# data[].attributes.body
The climate in Greece is predominantly Mediterranean. However, due to the country's unique geography, Greece has a remarkable range of micro-climates and local variations. To the west of the Pindus mountain range, the climate is generally wetter and has some maritime features. The east of the Pindus mountain range is generally drier and windier in summer.
Mediterranean climate occurs in the Aegean Islands, especially the Cyclades and the Dodecanese, southern and Evia, low-lying areas of Attica, the western, eastern and southern low-lying Peloponnese areas, and the low-lying areas of Crete.
During the summer, the weather is most frequently sunny and dry, and any precipitation falls in the form of showers or thunderstorms from cumuliform clouds. The air is usually hot during the day and pleasantly warm at night, but there are some very windy days, especially in the Cyclades islands and around them. Heatwaves may occur, but they are usually quite mild at the coastal areas, where temperatures are moderated by the relatively cooler sea and the sea breeze. Winters are wet and any snow that falls does not last long, especially in the south-facing slopes. Rain in winter is often persistent: The west areas of this climate zone receive a relatively higher amount of precipitation.
In Alpine Mediterranean, the winter is harsh with abundant snowfalls, while the summers are cool with frequent thunderstorms. This climate is to be found on high mountains, like Pindus and Rhodope.

# data[].attributes.attachment_link.description
A specialized agency of the United Nations whose mandate covers weather, climate and water resources. The UN's scientific voice on the state and behaviour of our atmosphere and climate.

# data[].attributes.body
Species of herpetofauna across the world differ from country to country. The same thing happens with the ecosystems that host these species.
In Alicia's HerpetoArea (BE), I found very intresting the "Atlas of Amphibians and Reptiles of the Brussels Region" from 2005, cheking the herpeto biodiversity. Something that in my country does not exist!
In Andrea's HerpetoArea (DE), there are 10 reptiles and 18 amphibians.
On the other side in Niki's HerpetoArea (ZA), there is no Urodela Order at all, with the black mamba (occured in this area) be the most dangerous reptile. There is also and the desert ecosystem.
In my HerpetoArea, there are 69 reptiles and 25 amphibias, with Viperidae family be the most dangerous.
With regard to ecosystems, the biggest difference is that there are more Mediterranean ecosystems in my HerpetoArea, while in the northern countries forests. Lakes, rivers ect, is the same.

# data[].attributes.body
I compare my plan to that of Ruth, Olympia, Dimitra's Herpetofauna.
To begin with, the identification that suggest Ruth, is very important to encounter the problem.
On the other hand, Olympia suggest to give the children tools to recognize the dangers and making them feel that they have the locus of control of their own environment. This is very interesting.
Last but not least, Sergia suggest to create posters and speeches for information purposes about Herpetofauna.
Summarizing, I can see to all of us a common purpose. The children must learn about the local Herpetofauna of them. They must stop seeing the reptiles as dangerous and fed animals.

# data[].attributes.body
https://www.radford.edu/~jtso/GeologyofVirginia/Fossils/GeologyOfVAFossils3-6c.html

# data[].attributes.body
Stellagama stellio

# data[].attributes.body
Podarcis erhardii

# data[].attributes.body
Podarcis peloponesiacus

# data[].attributes.body
Podarcis muralis

# data[].attributes.body
Podarcis milensis

# data[].attributes.body
Podarcis cretensis

# data[].attributes.body
Podarcis gaigaeae

# data[].attributes.body
Podarcis ionicus

# data[].attributes.body
Podarcis levendis

# data[].attributes.body
Podarcis tauricus

# data[].attributes.body
Lacerta agilis

# data[].attributes.body
Lacerta trilineata

# data[].attributes.body
Lacerta viridis

# data[].attributes.body
Vipera ammodytes

# data[].attributes.body
Vipera berus

# data[].attributes.body
Vipera graeca

# data[].attributes.body
Macrovipera schweizeri

# data[].attributes.body
Elaphe quatuorlineata

# data[].attributes.body
Zamenis situla

# data[].attributes.body
Natrix natrix

# data[].attributes.body
Erix jaculus

# data[].attributes.body
Platyceps najadum