Present and future distribution of the European pond turtle versus seven exotic freshwater turtles, with a focus on Eastern Europe

Abstract

Freshwater turtles are often used as terrarium pets, especially juveniles of exotic species. At the adult stage they are often released by their owners into the wild despite their high invasion potential. In Europe these thermophilic potentially invasive alien species occupy the habitats of the native European pond turtle Emys orbicularis (Linnaeus, 1758), with new records from the wild being made specifically in Eastern Europe (Latvia and Ukraine) during recent decades. Assessing the potential of alien freshwater turtles to establish in new territories is of great concern for preventing invasion risks while preserving native biodiversity in the present context of climate change. We explored this issue by identifying the present and future (by 2050) suitable habitats of the European pond turtle and several potentially invasive alien species of freshwater turtle already settled in Europe, using a geographic information system (GIS) modelling approach based on datasets from CliMond for climate, Near-global environmental information (NGEI) for freshwater ecosystems (EarthEnv) and Maxent modelling using open-access databases, data from the literature and original field data. Modelling was performed for seven species of alien freshwater turtles occurring from the extreme northern to southern borders of the European range of E. orbicularis: the pond slider Trachemys scripta (Thunberg and Schoepff, 1792), the river cooter Pseudemys concinna (Le Conte, 1830), the Florida red-bellied cooter Pseudemys nelsoni (Carr, 1938), the false map turtle Graptemys pseudogeographica (Gray, 1831), the Chinese softshell turtle Pelodiscus sinensis (Wiegmann, 1835), the Caspian turtle Mauremys caspica (Gmelin, 1774) and the Balkan terrapin Mauremys rivulata (Valenciennes, 1833). In Ukraine, the most Eastern limit of E. orbicularis distribution, were previously reported northern American originated T. scripta, M. rivulata, M. caspica, whereas in Latvia, Emys' most northern limit, were additionally reported P. concinna, P. nelsoni, G. pseudogeographica and Asia originated P. sinensis. The resulting Species Distribution Models (SDM) were of excellent performance (AUC > 0.8). Of these alien species, the most potentially successful in terms of range expansion throughout Europe were T. scripta (34.3% of potential range expansion), G. pseudogeographica (24.1%), and M. caspica (8.9%) and M. rivulata (4.3%) mainly in Eastern Europe, especially in the south of Ukraine (Odesa, Kherson, Zaporizhzhia regions, and Crimean Peninsula). Correlation between the built SDMs for the native E. orbicularis and the invasive alien T. scripta was reliably high, confirming the highly likely competition between these two species in places they cooccur. Moreover, a Multiple Regression Analysis revealed that by 2050, in most of Europe (from the western countries to Ukraine), the territory overlap between E. orbicularis and potentially invasive alien species of freshwater turtles will increase by 1.2 times, confirming higher competition in the future. Importantly, by 2050, Eastern Europe and Ukraine are predicted to be the areas with most suitable habitats for the European pond turtle yet with most limited overlap with the invasive alien species. We conclude that Eastern Europe and Ukraine are the most relevant priority conservation areas for the European pond turtle where it is now necessary to take protective measures to ensure safe habitat for this native species on the long-term.

Keywords: Climate change; Competition; Freshwater turtles; Habitat; Latvia; Priority conservation areas; Stacked species distribution models; Ukraine.

Fig. 1

Potential (probabilistic) distribution model of: E. orbicularis (a—current; b—2050); T. scripta (c—current; d—2050); Graptemys pseudogeographica (Gray, 1831) (e—current; f—2050); Pseudemys concinna (Le Conte, 1830) (g—current; h—2050); the map was created in SAGA GIS software using the open-source “World Administrative Boundaries—Countries and Territories” map under Open Government Licence v3.0 (WorldClim, Fig. S2—with “relief”).

Fig. 2

Potential (probabilistic) distribution model of: P. nelsoni (a—current; b—2050); P. sinensis (c—current; d—2050); Mauremys caspica (Gmelin, 1774) (e—current; f—2050); Mauremys rivulata (Valenciennes, 1833) (g—current; h—2050); the map was created in SAGA GIS software using the open-source “World Administrative Boundaries—Countries and Territories” map under Open Government Licence v3.0 (WorldClim, Fig. S3—with “relief”).

Fig. 3

Response curves generated by Maxent – bio01 (Annual mean temperature (°C)): (a) E. orbicularis, (b) T. scripta, (c) G. pseudogeographica, (d) P. concinna, (e) P. nelsoni, (f) P. sinensis, (g) M. capsica, (n) M. rivulata.

Fig. 4

Results of cluster analysis with the use of distributional data: Mc—M. capsica, Mr—M. rivulata, Pc—P. concinna, Pn—P. nelsonii, Ps—P. sinensis, Gp—G. pseudogeographica, Eo—E. orbicularis, Ts—T. scripta.

Fig. 5

As a result of Multiple Regression Analysis, models outcomes for: (A) area of intersection of current promising habitats for the native European pond turtle E. orbicularis and 7 exotic turtle species; (B) area of current promising habitat for E. orbicularis, without seven species of invasive aquatic turtles; (C) area of intersection of further promising habitats (to 2050) of all 8 species of turtles; (D) area of further promising habitats (to 2050) for E. orbicularis, without seven species of invasive aquatic turtles.