Assessing the relative potential ecological impacts and invasion risks of emerging and future invasive alien species
expand article infoJames W.E. Dickey, Ross N. Cuthbert, Michael Rea, Ciaran Laverty, Kate Crane, Josie South§, Elizabeta Briski|, Xuexiu Chang, Neil E. Coughlan, Hugh J. MacIsaac#, Anthony Ricciardi¤, Gillian E. Riddell, Meng Xu«, Jaimie T.A. Dick
‡ Queen’s University Belfast, Belfast, United Kingdom
§ South African Institute for Aquatic Biodiversity, Grahamstown, South Africa
| GEOMAR, Helmholtz-Zentrum für Ozeanforschung Kiel, Kiel, Germany
¶ Yunnan University, Kunming, China
# University of Windsor, Windsor, Canada
¤ McGill University, Montreal, Canada
« Pearl River Fisheries Research Institute, Chinese Academy of Fishery Sciences, Guangzhou, China
Open Access


Invasive alien species (IAS) cause myriad negative impacts, such as ecosystem disruption, human, animal and plant health issues, economic damage and species extinctions. There are many sources of emerging and future IAS, such as the poorly regulated international pet trade. However, we lack methodologies to predict the likely ecological impacts and invasion risks of such IAS which have little or no informative invasion history. This study develops the Relative Impact Potential (RIP) metric, a new measure of ecological impact that incorporates per capita functional responses (FRs) and proxies for numerical responses (NRs) associated with emerging invaders. Further, as propagule pressure is a determinant of invasion risk, we combine the new measure of Pet Propagule Pressure (PPP) with RIP to arrive at a second novel metric, Relative Invasion Risk (RIR). We present methods to calculate these metrics and to display the outputs on intuitive bi- and triplots. We apply RIP/RIR to assess the potential ecological impacts and invasion risks of four commonly traded pet turtles that represent emerging IAS: Trachemys scripta scripta, the yellow-bellied slider; T. s. troostii, the Cumberland slider; Sternotherus odoratus, the common musk turtle; and Kinosternon subrubrum, the Eastern mud turtle. The high maximum feeding rate and high attack rate of T. s. scripta, combined with its numerical response proxies of lifespan and fecundity, gave it the highest impact potential. It was also the second most readily available according to our UK surveys, indicating a high invasion risk. Despite having the lowest maximum feeding rate and attack rate, S. odoratus has a high invasion risk due to high availability and we highlight this species as requiring monitoring. The RIP/RIR metrics offer two universally applicable methods to assess potential impacts and risks associated with emerging and future invaders in the pet trade and other sources of future IAS. These metrics highlight T. s. scripta as having high impact and invasion risk, corroborating its position on the EU list of 49 IAS of Union Concern. This suggests our methodology and metrics have great potential to direct future IAS policy decisions and management. This, however, relies on collation and generation of new data on alien species functional responses, numerical responses and their proxies, and imaginative measures of propagule pressure.


ecological impacts, functional response, invasive alien species, numerical response, pet propagule pressure, relative impact potential, relative invasion risk, risk assessment


Invasive alien species (IAS), i.e. those introduced to areas outside their native range (which may or may not then have impact; see Ricciardi and Cohen 2007), can be major drivers of global biodiversity loss and cause a range of other negative impacts (Tilman et al. 2017). While many species fail to establish, i.e. are not capable of reproducing and becoming self-sustaining populations, some do establish, spread and exert strong, negative ecological effects (Colautti et al. 2014; Dick et al. 2017a, b). Propagule pressure, the number, frequency and viability of individuals of each species released, is a key factor determining whether a species establishes, with the aquarium and ornamental trades shown to be a major determinant of propagule pressure (Gertzen et al. 2008; García-Díaz et al. 2015). This is deemed responsible for a third of aquatic IAS (Padilla and Williams 2004), due to poor industry regulation (Raghavan et al. 2013) and fuelled by misconceptions surrounding the ethics and consequences of “mercy” release (Liu et al. 2013). The global trade of reptiles and amphibians in particular has led to the spread of some high profile IAS with serious ecological impacts (Kraus 2015). The numbers of alien species introductions are unlikely to saturate in the near future (Seebens et al. 2018), but understanding, assessing and predicting the associated ecological impacts has proven difficult, especially for species without invasion histories. Some high profile failed attempts at controlling IAS (Courchamp et al. 1999; Rayner et al. 2007) have led to a consensus being reached that prevention is the best option for IAS management (Piria et al. 2017). However, prevention strategies require new methods to quantify likely ecological impact, here defined as negatively affecting the abundance of one or more native species, and invasion risk, which combines impact with likelihood of establishment, for emerging and potential future invaders (Dick et al. 2017b).

Comparative functional responses (CFRs) have been successful in characterising damaging IAS and have proven predictive for those without invasion impact history (Dick et al. 2014, 2017b). CFR can also take into account a wealth of context-dependencies and interactive effects, such as dissolved oxygen levels (Laverty et al. 2015), habitat complexity (Wasserman et al. 2016), temperature (Xu et al. 2016; Cuthbert et al. 2018), higher order predators (Barrios-O’Neill et al. 2014) and presence of parasites (Laverty et al. 2017b). However, despite the many advantages of the CFR methodology, the impacts of a consumer on resources (e.g. predator on prey) will clearly be the product of these per capita effects and the population response of the consumer (i.e. the numerical response (NR); Solomon, 1949). A simpler proxy for NR is the abundance of a species, which has recently been combined with CFR into the Relative Impact Potential (RIP) metric to successfully predict invader ecological impact (Dick et al. 2017b).

RIP might prove particularly valuable for the study of IAS emerging from the pet trade, with global trade of freshwater turtles in particular a pressing problem (Nori et al. 2017). For example, Trachemys scripta elegans (the red-eared slider) appears both on the EU List of IAS of Union concern and the IUCN’s 100 Worst Invasive Species list, having invaded 73 countries (García-Díaz et al. 2015; Capinha et al. 2017) and living on every continent except Antarctica (Rödder et al. 2009). This species was once widely traded due to its small size as juveniles, ease of maintenance, relative affordability (Teillac-Deschamps et al. 2009) and, remarkably, the global craze of Teenage Mutant Ninja Turtles (Somma et al. 2009a). Upon release, T. s. elegans has been shown to compete with native turtles for food and basking habitat (Cadi and Joly 2003; Pearson et al. 2015). We thus urgently require methods to predict the impacts of these and similar potential IAS that currently lack both invasion history and invasion range population abundance data.

Here, we utilise the RIP metric and biplots of Dick et al. (2017b) and Laverty et al. (2017b) and develop these with alternative available proxies for the numerical response, namely “lifespan”, “fecundity” and their product, “lifetime fecundity”. These proxies may not necessarily be accurate reflections of the numerical response per se, but rather are useful for comparative purposes, as were abundance and biomass in the original RIP metric (Dick et al. 2017b). Unlike the true numerical response and even abundance/biomass (see Dick et al. 2017b), the three measures above are generally available in literature for most species and are so for the current study species. Further, as fecundity is clearly linked to the numerical response and lifespan may determine ecological impact over time, these proxies, in a comparative sense, are useful multipliers of per capita effects (see Dick et al. 2017b). Also included in the present study is the functional response “attack rate”, which offers insights into predatory impact at low prey densities (i.e. destabilising Type II FRs), as well as the “maximum feeding rate”, calculated as the reciprocal of the “handling time” (Dick et al. 2017b). Since propagule pressure is a major determinant of invasion risk (Briski et al. 2012) and possibly impact (Ricciardi and Cohen 2007), we further modify the RIP biplots of Laverty et al. (2017b) to include propagule pressure on a third axis (i.e. triplots) to give a combined measure of Relative Invasion Risk (RIR). Using these metrics, we investigate the relative potential ecological impacts and invasion risks of four commonly traded freshwater turtles: two Trachemys scripta subspecies that have recently been added to the EU list of 49 IAS of Union Concern, namely T. s. scripta, the yellow-bellied slider and T. s. troostii, the Cumberland slider; as well as Sternotherus odoratus, the common musk turtle and Kinosternon subrubrum, the Eastern mud turtle, both of which are also widely available in the pet trade.


Study species

Trachemys scripta scripta, T. s. troostii, Sternotherus odoratus and Kinosternon subrubrum originate from North America and are sold around the world (Polo-Cavia et al. 2011), with T. scripta and S. odoratus being amongst the most commonly imported turtle species into the UK and Ireland (García-Díaz et al. 2014). Their longevity and lack of predators and competitors mean that, even if current environmental conditions prevent reproduction in some of their introduced ranges, there is the potential for ecological impact and invasion risk. This may be exacerbated if climate change and/or adaptation leads to these populations becoming viable (Bugter et al. 2011).

T. s. scripta and T. s. troostii can live for 36 years (Frazer et al. 1990) and reproduction is expected to be similar to T. s. elegans, with females laying up to five clutches per year, with between 2 and 23 eggs in each clutch (Somma et al. 2009a). Sternotherus odoratus can live for 30 years (Bugter et al. 2011), with females laying 2 to 4 clutches per year, with between 1 and 9 eggs per clutch (Somma and Fuller 2009). Kinosternon subrubrum can live for 46 years (Frazer et al. 1990) and females tend to lay between 1 and 3 clutches per year, with 1 to 6 eggs per clutch (Somma et al. 2009b).

Animal collection and maintenance

The four turtle species were provided by Maidenhead Aquatics, Northern Ireland (carapace lengths 35–50 mm; mixed sexes) and maintained in holding tanks containing a water heater (150W Eheim thermocontrol, Germany) and water cooler to ensure water temperature was maintained at 16 °C. Two basking platforms and basking lights created a hot spot of 23 °C, controlled by an automatic temperature controller (Habistat Classic, England). Substrate (0.8 mm grain size) was added to the bottom 30 mm of each holding tank (JBL, Germany). Nine individuals of each species were acquired for each experimental batch and each species received its own holding tank and was quarantined for one week prior to experiments. During this period, no illness or deaths were recorded and the animals were fed daily with commercial floating turtle food (JBL, Germany). Before experimental FR trials, all turtles were starved for 24 hours to standardise hunger levels. Focal prey, the amphipod crustacean Gammarus pulex (15–17mm body length; unparasitised), upon which all turtle species were observed to feed readily and represents a general prey item, were collected from the Minnowburn River, N. Ireland (N54.546, W5.594) two days before the experiments and acclimatised to the experimental temperature.

FR procedure

Experiments were performed 22 February–27 April 2016. Experimental tanks (250 mm × 120 mm × 90 mm) with 30 mm of substrate at 16 °C were supplied with prey 15 minutes prior to the turtles being introduced. Prey densities were 2, 4, 8, 16, 32, 64, 128, 256 (n = 6 per experimental group). For each turtle species, individuals were randomly selected and assigned to a random prey density and allowed to feed for thirty minutes. Controls were performed for each prey density (n = 3 each) with the same experimental conditions but in the absence of turtle predators, to quantify prey mortality for any other reasons.

Statistical methods

Data were analysed using R version 3.2.3. (R Core Team 2015). Logistic regression of the proportion of prey killed as a function of prey density was used to discern functional response types (see Juliano 2001). Where a significant negative first order linear coefficient was detected, a Type II response was ascribed; conversely, a Type III form was considered when a significant positive first order linear coefficient was followed by a significant negative second order coefficient (Juliano 2001). Rogers’ random predator equation for non-replacement of prey was applied to estimate values of ‘h’ (handling time) and ‘a’ (attack rate):

Ne=N0 (1–exp (a (Neh–T))) (1)

where Ne is the number of prey eaten, N0 is the initial density of prey, h is the handling time, a is the attack constant and T is the total experimental period. Model fitting used the Lambert W function (Bolker 2008) in R due to the implicit nature of the random predator equation. This relates to the fact that the random predator equation is not solvable on its own and requires the Lambert W function for this. See the opening paragraph of Data were bootstrapped (n = 30) to calculate multiple estimates of the handling time h, maximum feeding rate (1/h) and attack rate a with standard error (SE).

Metrics and measures

Relative Impact Potential (RIP) was originally developed using population abundance/density/biomass as a proxy for the consumer numerical response (NR: Dick et al. 2017b), but in certain situations such data are lacking, as, for example, with potential future invasions by turtles and other potential IAS. Since our metrics are comparative rather than absolute, we substituted the NR with other life history measures of the invaders that may influence impact. Here, we thus propose three other alternative proxies for the numerical response:

(1) Lifespan (L). With temperatures in many temperate regions high enough for turtles to survive in the wild for many years, but as yet too low to facilitate reproduction (Teillac-Deschamps et al. 2009; but see Standfuss et al. 2016), lifespan offers a proxy multiplier of per capita effects (when other proxies such as abundance are not available), giving a metric quantifying the ecological impact that a species may exert over space and time. Thus, Impact Potential utilising lifespan data (IPL) is:

IPL = FR × L (2)

where FR is the functional response (estimated maximum feeding rate, 1/h, from equation 1) and L is maximum lifespan (Table 1);

Numerical response proxies of lifespan, fecundity and lifetime fecundity, plus Pet Propagule Pressure (PPP; see Text and Table 3) for our study turtle species. Pet Propagule Pressure (Northern Ireland, NI) is a function of the proportional availability of each species across pet shops and the total number of pet shops surveyed and PPP (Great Britain, GB) is the proportional availability of the four species based on online classified advertisements and the total number of online advertisements surveyed.

Turtle Lifespan (L) (maximum years) Ref. Fecundity (F) (eggs per year) Ref. Lifetime fecundity (LxF) Pet Prop. Press. (PPPNI) Pet Prop. Press. (PPPGB)
T. s. scripta 36 Frazer et al. 1990 115 Somma et al. (2009a) 4140 0.05 0.1
T. s. troostii 36 Frazer et al. 1990 115 Somma et al. (2009a) 4140 0.00 0.05
S. odoratus 30 Bugter et al. 2011 36 Somma and Fuller (2009) 1080 0.30 0.39
K. subrubrum 46 Frazer et al. 1991 18 Somma et al. (2009b) 828 0.05 0.02

(2) Fecundity (F). Where reproduction of the turtles occurs or may occur in future, fecundity offers another proxy multiplier of per capita effects, since reproductive output is clearly an element of the true numerical response. Thus, Impact Potential utilising fecundity data (IPF) is:

IPF = FR × F (3)

where FR is as above and F is the product of clutch size and number of clutches per annum (Table 1);

(3) Lifespan fecundity (LF). Where suitable data are available, a third proxy for the numerical response may be constructed as the product of maximum lifespan and fecundity, that is lifetime fecundity (LF), as this captures both reproductive output per bout and over time and thus Impact Potential is:

IPLF = FR × LF (4)

The RIP calculations of Dick et al. (2017b) and RIP biplots of Laverty et al. (2017b) use maximum feeding rate (1/h) as the FR measure (see above), based on the curve parameter “h” (handling time). Here, we propose the additional use of the other classic FR parameter, “a” (attack rate), since this quantifies the initial gradient of the functional response curve, giving insights into the critical population level impact that a consumer can exert at low resource densities (Dick et al. 2014). Thus, Impact Potential utilising attack rate and lifespan data (IPL) is:

IPL = a × L (5)

and with fecundity data is:

IPF = a × F (6)

and with lifetime fecundity:

IPLF = a × LF (7)

In addition, to enter the propagule pressure argument to measure overall invasion risk, we qualify each IP equation with Pet Propagule Pressure (PPP). We propose two PPP methods. First, we quantified the availability of the four species on a local (i.e. Northern Ireland, NI) level via a survey of twenty pet shops between the 31 January and 1 March 2017 (Suppl. material 1: Table S1; PPP values as per Table 1). PPP in this context is calculated as follows:

PPP (NI) = Np/Tp (8)

where Pet Propagule Pressure (Northern Ireland, NI) is a function of the proportional availability of each species across pet shops (Np) and the total number of pet shops surveyed (Tp).

The second version of PPP involved a survey of online classified advertisements (Suppl. material 2: Table S2; PPP values as per Table 1). A major advantage of this method is that each advertisement represents an unwanted pet, a key feature of pet releases into the wild. This second survey was conducted on 21 November 2017 and used two websites, and Search terms of ‘terrapin’ and ‘turtle’ were used for both websites and screenshots were taken of each advertisement claiming to be selling our study species. Using the location of the seller and the screenshot images, we prevented the double-counting of advertisements. The PPP values were calculated as follows:

PPP (GB) = Na/Ta (9)

where PPP (Great Britain, GB) is the proportional availability of the four species based on online classified advertisements (Na) and the total number of online advertisements surveyed (Ta).

By incorporating these two measures of propagule pressure, the Impact Potential (IP) equations (equations 2–7) can incorporate both risk of introduction and its ecological consequences to become Invasion Risk (IR):

IRL = FR × L × PPP (10)

IRF = FR × F × PPP (11)

IRLF = FR × LF × PPP (12)

IRL = a × L × PPP (13)

IRF = a × F × PPP (14)

IRLF = a × LF × PPP (15)

We present biplots to illustrate Relative Impact Potential (equations 2–7) and triplots for Relative Invasion Risk (equations 10–15) of the four turtle species to give visual representations of relative ecological impact and invasion risk (see Suppl. material 3, 4 for R scripts, and Suppl. material 5 for associated .csv file).


Prey survival in control treatments was 98–100%, therefore mortality during FR experiments was attributed to predation, which was also directly observed. Type II functional responses were observed for all turtle species (Table 2; Figure 1). Handling times were lowest and, hence, maximum feeding rates were highest, for the two T. scripta subspecies, with the order T. s. scripta > T. s. troostii > K. subrubrum > S. odoratus (Fig. 2a, b). T. s. scripta had the highest attack rate, S. odoratus the lowest and T. s. troostii and K. subrubrum were intermediate (Fig. 2c), with K. subrubrum having a higher attack rate, but also higher handling time (i.e. lower maximum feeding rate) than T. s. troostii (Fig. 2a–c).

First order terms calculated from logistic regression to denote functional response type across all predator treatments. The significant negative first order term values across all four turtles indicate Type II functional responses for each predator. Handling time (h), maximum feeding rate (1/h) and attack rate (a) parameter mean estimates (bootstrapped, n = 30), derived using Rogers’ random predator equation (eqn 1).

Predator First term, P Handlingtime, h Maximum feeding rate, 1/h (G. pulex consumed per 30 mins) Attack rate, a
T. s. scripta -0.011, <0.001 0.027 37.036 2.678
T. s. troostii -0.011, <0.001 0.028 35.405 2.038
S. odoratus -0.011, <0.001 0.039 25.468 1.847
K. subrubrum -0.012, <0.001 0.037 27.142 2.314
Figure 1.

Functional responses of T. s. scripta, T. s. troostii, S. odoratus and K. subrubrum towards G. pulex prey. Values are mean ±SE.

Figure 2.

Parameter estimates (±SE) of: a handling time h b maximum feeding rate 1/h c attack rate a, for bootstrapped (n = 30) Type II functional response curves of T. s. scripta (Tss), T. s. troostii (Tst), S. odoratus (So) and K. subrubrum (Ks) towards G. pulex prey.

The numerical response proxy values are given in Table 1 and the IP and IR values in Table 3. Impact Potential (IP) scores using maximum feeding rate (see Tables 2, 3) with all three numerical response proxies were higher for T. s. scripta and T. s. troostii relative to K. subrubrum and S. odoratus (equations 2–4; Table 3; Fig. 3a–c). However, when attack rate was used with lifespan, the IP score was highest for K. subrubrum, closely followed by T. s. scripta and T. s. troostii and then S. odoratus (equation 5; Table 3; Fig. 3d). If we consider fecundity and lifetime fecundity, IP was again highest for T. s. scripta and T. s. troostii versus S. odoratus and K. subrubrum, which were similar (equation 6, 7; Table 3; Fig. 3e, f).

Impact Potential (IP) and Invasion Risk (IR) calculations, whereby: IP(FR) = Maximum feeding rate (FR) × NRproxy i.e. lifespan (L), fecundity (F) or lifetime fecundity (LF); IP(a) = Attack rate (a) × NRproxy; IR(FR) = IP(FR) × Pet Propagule Pressure (PPP); IR(a) = IP(a) × PPP. PPPNI (Northern Ireland) is a function of the proportional availability of each species across pet shops and the total number of pet shops surveyed and PPPGB (Great Britain) is the proportional availability of the four species based on online advertisements and the total number of advertisements surveyed.

T. s. scripta 1,333.30 4,259.14 153,329.04 96.41 307.97 11,086.92
T. s. troostii 1,274.76 4,072.15 146,597.40 73.37 234.37 8,437.32
S. odoratus 764.04 916.85 27,505.44 55.41 66.49 1,994.76
K. subrubrum 1,248.53 488.56 22,473.58 106.44 41.65 1,915.99
T. s. scripta 66.66 212.96 7,666.45 4.82 15.40 554.35
T. s. troostii 0 0 0 0 0 0
S. odoratus 229.21 275.05 8,251.63 16.62 19.95 598.43
K. subrubrum 62.43 24.43 1,123.68 5.32 2.08 95.80
T. s. scripta 133.33 425.91 15,332.90 9.64 30.80 1,108.69
T. s. troostii 63.74 203.61 7,329.87 3.67 11.72 421.87
S. odoratus 297.98 357.57 10,727.12 21.61 25.93 777.96
K. subrubrum 24.97 9.77 449.47 2.13 0.83 38.32
Figure 3.

Biplots showing Relative Impact Potential of T. s. scripta, T. s. troostii, S. odoratus and K. subrubrum towards G. pulex prey. Impact potential calculated as a product of maximum feeding rate and lifespan (a), fecundity (b) and lifetime fecundity (c); then attack rate and lifespan (d), fecundity (e) and lifetime fecundity (f). Impact increases from bottom left to top right.

Pet Propagule Pressure (PPP) of each species was similar in both the Northern Ireland (NI) and Great Britain (GB) surveys, with respective orders of S. odoratus > T. s. scripta = K. subrubrum > T. s. troostii and S. odoratus > T. s. scripta > T. s. troostii > K. subrubrum (see Table 1 and details of the two surveys are outlined in Suppl. material 1: Table S1 and Suppl. material 2: Table S2).

For all six of the NI Relative Invasion Risk (RIR) triplots, values were highest for S. odoratus, with the order S. odoratus > T. s. scripta > K. subrubrum > T. s. troostii across all derivations of RIR (equations 10–15; Table 3; Fig. 4a–f). For the GB RIR triplots, the order using both lifespan calculations was S. odoratus > T. s. scripta > T. s. troostii > K. subrubrum (Fig. 5a, d), but for all other RIR metrics, the order was T. s. scripta > S. odoratus > T. s. troostii > K. subrubrum (equations 10–15; Table 3; Fig. 5b, c, e, f).

Figure 4.

Triplots showing Relative Invasion Risk of T. s. scripta, T. s. troostii, S. odoratus and K. subrubrum in a Northern Irish context. Invasion Risk calculated as a product of maximum feeding rate and Pet Propagule Pressure (PPP) with lifespan (a), with fecundity (b) and with lifetime fecundity (c); then attack rate and PPP with lifespan (d), with fecundity (e) and with lifetime fecundity (f). PPP for each species calculated by surveying 20 local pet shops and determining proportions of each species sold. Invasion Risk increases from bottom left to top right of each plot, with species ranked 1–4.

Figure 5.

Triplots showing Relative Invasion Risk of T. s. scripta, T. s. troostii, S. odoratus and K. subrubrum in a Great British context. Invasion Risk is calculated as a product of maximum feeding rate, lifespan and Pet Propagule Pressure (PPP) (a), with fecundity (b) and with lifetime fecundity (c); then attack rate, lifespan and PPP (d), with fecundity (e) and with lifetime fecundity (f). PPP for each species calculated by surveying classified advertisements online and finding what proportion were selling the species in question. Invasion Risk increases from bottom left to top right of each plot, with species ranked 1–4.


Invasion ecology has long lacked a unifying methodology that predicts ecological impacts and overall invasion risks of invasive species (Dick et al. 2017a). Propagule pressure alone is often unable to incorporate context-dependencies (Dick et al. 2017a) and invasion history is not always available for emerging and potential future invaders (Kulhanek et al. 2011). Furthermore, no single species trait or combination of traits offers consistent predictive efficacy (Ricciardi et al. 2013), although comparative functional responses (CFR) have, up until now, given the most effective predictions of invader ecological impact across a wide range of abiotic and biotic contexts (Dick et al. 2014, 2017b). These relationships between resource availability and resource uptake rate have proven robust in comparing invaders with trophically analogous natives or other invaders, even in laboratory conditions that do not necessarily mimic natural conditions (Laverty et al. 2017b). Here, we found that each turtle species followed classic Type II functional responses, whereby, at low prey densities, high proportions of prey are consumed and hence prey populations may be destabilised. Thus, the introduction, establishment and spread of alien turtles could have marked top-down impacts on aquatic fauna in freshwater ecosystems. We predict that T. s. scripta and T. s. troostii will have the greatest ecological impacts. The two T. scripta subspecies were shown to have the highest maximum feeding rates, with T. s. scripta and K. subrubrum shown to have the highest attack rates, an indication of impact at low prey densities. T. s. scripta displayed a “perfect storm” of highest maximum feeding rates and highest attack rates, suggesting potential for significant impacts upon freshwater invertebrate communities.

While FR analyses have been used mainly to determine the impacts of alien predators (see Dick et al. 2014; Alexander et al. 2014; Dick et al. 2017b), they may also be used to assess impacts of herbivores (e.g. Xu et al. 2016) or any other trophic/taxonomic group (Dick et al. 2014). Additionally, FRs can inform about the competitive abilities of species and, indeed, competition theory was heavily based on the “functional resource-utilisation responses” of plants (see Tilman 1977; Dick et al. 2017a). While the use of comparative FRs in assessing inter-specific competition amongst animal taxa has been slow in its adoption, we foresee development of this method and our metrics to assess ecological impact through competition, as occurs between alien and native turtles (Cadi and Joli 2003; Polo-Cavia et al. 2011).

The Impact Potential (IP) and Invasion Risk (IR) metrics, plus our illustrative bi- and triplots giving Relative IP and Relative IR, retain the benefits of CFR, but bolster these per capita measures with proxies for the numerical response (NR), that is, the consumer population response. Emerging invaders may lack NR data and have no data for their abundances/densities/biomass in potential invasion regions. With these latter NR proxies not available, we hence require alternatives and use lifespan, fecundity and lifetime fecundity as comparative multipliers of per capita effects. The resulting impact potentials were subsequently combined with values for our two versions of Pet Propagule Pressure (PPP) to give Invasion Risk (IR), which assesses which species are currently the most likely candidates for introduction, combined with potential impact. Using these approaches, we determined that T. s. scripta and T. s. troostii have the highest RIP, but the more commonly traded S. odoratus has the greatest RIR and thus should be of great concern. Our approach illustrates the potential use of combinatorial metrics to guide policy and intervention and exploits inherent life-history traits of invaders with their feeding impacts and their likelihood of introduction.

The use of maximum lifespan as a numerical response proxy offers a readily available multiplier of per capita effects. With reproduction by these turtles not yet possible in most temperate regions, the longer the species survive in the wild, the greater the ecological impact that will accrue. One caveat is the combination of unfamiliar climate, flora and fauna, combined with invader naiveté, which could alter the estimations of lifespan should these species be released into the wild. A caveat exists for lifetime fecundity too, as fecundity does not remain constant over the course of a lifetime. However, as our metric is comparative and all species should be equally affected, such reductions may not affect predictions of relative ecological impact. Here, using lifespan, we find that T. s. scripta and T. s. troostii had the highest impact potentials on the maximum feeding rate biplot, while K. subrubrum and T. s. scripta had the highest and second highest impacts on the attack rate biplot. However, as illustrated by Dick et al. (2017b), it is maximum feeding rate combined with NR proxies that give highest success in ecological impact prediction and hence we conclude here that T. s. scripta and T. s. troostii will have the greatest ecological impacts as invaders.

Our second NR proxy, fecundity, defined as the number of offspring born over a given period of time (Lamb et al. 2009), can be a key determinant of whether or not an introduced population establishes and subsequently persists (Pöckl 2007). The inability to reproduce in northern European climates has dampened the impact of T. scripta elegans, but where this barrier to reproduction does not apply, this and other introduced reptiles have strong impacts (Kraus 2015). Our biplots illustrate the two T. scripta subspecies having the highest of such impact values, with the higher maximum feeding rate and attack rate of T. s. scripta indicative of higher impact at both high and low prey densities. Our third NR proxy, lifetime fecundity, clearly distinguished the two T. scripta subspecies from the other turtles, with T. s. scripta having the greatest potential impact.

Our novel Relative Invasion Risk (RIR) triplots used Pet Propagule Pressure (PPP) to give a third dimension for invasion risk assessment. The first PPP calculations are based on a survey of 20 pet shops, ranging from small independent traders to UK-wide chains across Northern Ireland (NI). This offered vital data on which species are currently being sold in the NI pet trade and, by proxy, which species are likely candidates for future release and escape (Bugter et al. 2011). Pet shop availability has previously been used as an indicator of propagule pressure (Rixon et al. 2005; Gertzen et al. 2008), but the combination of this measure with impact potential allows users to counteract the shortcomings of propagule pressure for invasion risk prediction on its own (Dick et al. 2017a). Our survey found that S. odoratus was the most readily available of the four, with K. subrubrum and T. s. scripta featuring rarely. T. s. troostii, acquired for the study in the prior year, was not found in the survey. From all six of the triplots, the prevalence of S. odoratus in NI gave it the highest RIR of the four species, with all of the NR proxies and both maximum feeding rate and attack rate parameters. T. s. scripta came second to S. odoratus in all contexts with the exception of attack rate-lifespan, for which K. subrubrum had the second highest RIR.

The second measure of PPP was derived from online classified advertisements for unwanted pet turtles in GB and showed a similar result to that of the NI pet shop survey, with S. odoratus again found to be much more available than the other three turtle species. S. odoratus had the highest RIR for the lifespan calculations, but T. s. scripta, found to be twice as common in GB as it was in NI, had the highest RIR when fecundity and lifetime fecundity were taken into account. Trachemys scripta troostii was also more available than in the prior NI survey and, as a result, poses a greater risk in GB. Monitoring needs to occur in the future as changes in supply and demand will lead to the study species shifting their relative availabilities (Kitowski and Pachol 2009) and, as a result, RIR. Indeed, this novel metric illustrates the need to reduce propagule pressure driven by the pet trade and private holdings of such species.

Gertzen et al. (2008) cited size and aggression as major reasons for fish releases and similarly ownership difficulty was found to be a key determinant for establishment success of exotic reptiles (Fujisaki et al. 2010; García-Díaz et al. 2015). T. scripta spp. have been shown to be aggressively territorial (Polo-Cavia et al. 2011), with males reaching sizes of up to 24 cm and females 29 cm (van Dijk et al. 2011). Though smaller, S. odoratus still reach 14 cm (van Dijk 2015) and K. subrubrum 12 cm (van Dijk 2011). The likelihood of outgrowing their tanks, combined with longevity and their use in religious ceremonial release (Liu et al. 2013), suggests that all four are likely candidates for release.

Assessing potential for long-term impact requires information on which species will likely establish. Temperature is crucial for embryonic development and offspring phenotype (Booth 2006) and there are two types of temperature-dependent sex determination (TSD). Species with TSD II, female-male-female, might have an invasion advantage, producing females at both cool and warm temperatures with males at intermediate temperatures (Ewert et al. 1994). S. odoratus has TSD II, with the lowest pivotal temperature at 20–21 °C (Bugter et al. 2011), whereas the T. scripta subspecies exhibit TSD Ia, male-female, with the sole pivotal temperature of 29–29.5 °C (Wibbels et al. 1998). However, potential to adapt to local thermal conditions (Zhao et al. 2015) and maternal behaviour buffering environmental conditions (Jackson 1988) highlight the need for effective monitoring across Europe.

The Relative Impact Potential metric, in its original form (Dick et al. 2017b), was based on the total response equation, i.e. the product of functional and numerical responses. With the complexity of calculating numerical responses well documented, using proxies of population abundance, density or biomass has been suggested (Dick et al. 2017b). In situations where such data are lacking, or inappropriate, such as when a species has no invasion history and using native population abundance could be misleading, there is a need to find other numerical response proxies. While we addressed longevity and fecundity, alternatives exist, for example, climatic similarity between native and potential invaded ranges (Filz et al. 2018), age of sexual maturity, embryological development time and number of generations per year (Grabowski et al. 2007). Further, embryonic development and hence key life history traits can be affected by toxicants (Zi et al. 2018), with the potential for such effects to vary amongst invaders and natives. We also highlight the potential for alternative proxies for propagule pressure, with the two methods applied here appropriate for species already available in the pet trade of the recipient area, but of limited use for those yet to arrive. We open the metric to alternative sources such as the US Law Enforcement Management Information System (LEMIS) database which records live wildlife imports and exports (Fujisaki et al. 2010; García-Díaz et al. 2015; Tingley et al. 2016) or studies from the native range that assess survivability during transit in hostile conditions, e.g. ship ballast water (Gollasch et al. 2000). With no likely single predictor of invasion risk, we propose that combining functional responses with proxies for numerical responses, plus propagule pressure, offers an effective three-pronged assessment that spans the invasion process.


Using our impact potential metric, the turtle warranting management priority is Trachemys scripta scripta. While uncertainty surrounds the ability and timeframe of all four turtles to adapt to more temperate climates, the potential for high relative clutch size, high feeding rate on a locally abundant prey, large body size and aggression to ward off potential competitors and predators, mean the likelihood of establishment and ecological impact of T. s. scripta is high. However, the widespread availability of S. odoratus, combined with a wide range of habitats, TSDII reproduction and the same population destabilising Type II functional response as T. s. scripta, mean the Relative Invasion Risk (RIR) triplots highlight this as a species that would otherwise have been overlooked solely on the basis of Comparative Functional Response (CFR) and Relative Impact Potential (RIP) studies.

With the pet trade likely to continue to be the main driver of any turtle species arriving, either by release or escape, knowing the species being imported into local pet shops and being sold by owners, is of vital importance. For that reason, the RIR triplots, which combine IP data with a measure of propagule pressure (i.e. PPP), offer an informative way of prioritising potential invasive species for management interventions. Going forward, there is vital need for regular surveys and for assessing the potential impact and risks of newly arrived specimens. With some turtle species encountered in both PPP surveys that were not investigated in this study, what we have provided is an avenue for further research and a starting point for the compilation of a user-friendly database of potential pet shop invaders to help decision-makers worldwide to assess IAS impact and invasion risk. More broadly, in terms of future research and management directions, there is a need for compilation of data on alien species functional responses, that is, existing and new data similar to those collated for biocontrol agents. We also need better estimates of numerical responses, as this latter measure is the gold standard for providing the total response and hence impact of invaders (see Dick et al. 2017b). Further, where our numerical response proxies are used in metrics, more ground-truthing of their predictive capacities is required, although this may only come with new invasions and monitoring of their actual impacts. Propagule measures also need more careful and imaginative derivations and examined for actual performance against real invasions.


JTAD, MR, CL, KC, JS and GR conceived the study, performed the experiments and produced initial results; JWED and RNC conducted statistical analyses; JWED, RNC and JTAD prepared the initial manuscript, led by JWED. AR, HJM, NEC, MX, XC and EB contributed vital input to the development of concepts. All contributed to production of final manuscript and gave approval for publication. JWED supported by Inland Fisheries Ireland (IFI), RNC by Department for the Economy Northern Ireland, NEC by Environmental Protection Agency (EPA), EB by Alexander von Humboldt Sofja Kovalevskaja Award and HJM and TR by NSERC Discovery grants. Thanks also to the Natural Environment Research Council (NERC). The authors would like to give credit to Pablo García-Díaz for insightful, constructive comments that improved the quality of the manuscript.


  • Barrios-O’Neill D, Dick JTA, Emmerson MC, Ricciardi A, MacIsaac HJ, Alexander ME, Bovy HC (2014) Fortune favours the bold: A higher predator reduces the impact of a native but not an invasive intermediate predator. Journal of Animal Ecology 83: 693–701.
  • Bolker BM (2008) Emdbook: ecological models and data in R. Princeton University Press, Princeton.
  • Booth DT (2006) Influence of incubation temperature on hatchling phenotype in reptiles. Physiological and Biochemical Zoology 79: 274–281.
  • Briski E, Bailey SA, Casas-Monroy O, DiBacco C, Kaczmarska I, Levings C, MacGillivary ML, McKindsey CW, Nasmith LE, Parenteau M, Piercey GE, Rochon A, Roy S, Simard N, Villac MC, Weise AM, MacIsaac HJ (2012) Relationship between propagule pressure and colonization pressure in invasion ecology: a test with ships’ ballast. Proceedings of the Royal Society B: Biological Sciences 279: 2990–2997.
  • Cadi A, Joly P (2003) Competition for basking places between the endangered European pond turtle (Emys orbicularis galloitalica) and the introduced red-eared slider (Trachemys scripta elegans). Canadian Journal of Zoology-Revue Canadienne De Zoologie 81: 1392–1398.
  • Capinha C, Seebens H, Cassey P, García-Díaz P, Lenzner B, Mang T, Moser D, Pyšek P, Rödder D, Scalera R, Winter M, Dullinger S, Essl F (2017) Diversity, biogeography and the global flows of alien amphibians and reptiles. Diversity and Distributions 23: 1313–1322.
  • Cuthbert RN, Dick JTA, Callaghan A, Dickey JWE (2018) Biological control agent selection under environmental change using functional responses, abundances and fecundities; the Relative Control Potential (RCP) metric. Biological Control 121: 50–57.
  • Dick JTA, Alexander ME, Jeschke JM, Ricciardi A, MacIsaac HJ, Robinson TB, Kumschick S, Weyl OLF, Dunn AM, Hatcher MJ, Paterson RA, Farnsworth KD, Weyl OLF, Richardson DM (2014) Advancing impact prediction and hypothesis testing in invasion ecology using a comparative functional response approach. Biological Invasions 16: 735–753.
  • Dick JTA, Alexander ME, Ricciardi A, Laverty C, Downey PO, Xu M, Jeschke JM, Saul W-C, Hill MP, Wasserman R, Barrios-O’Neill D, Weyl OLF, Shaw RH (2017a) Functional responses can unify invasion ecology. Biological Invasions: 19: 1667–1672.
  • Dick JTA, Laverty C, Lennon JJ, Barrios-O’Neill D, Mensink PJ, Britton JR, Médoc V, Boets P, Alexander ME, Taylor NG, Dunn AM, Hatcher MJ, Rosewarne PJ, Crookes S, MacIsaac HJ, Xu M, Ricciardi A, Wasserman RJ, Ellender BR, Weyl OLF, Lucy FE, Banks PB, Dodd JA, MacNeil C, Penk MR, Aldridge DC, Caffrey JM (2017b) Invader Relative Impact Potential: a new metric to understand and predict the ecological impacts of existing, emerging and future invasive alien species. Journal of Applied Ecology 54: 1259–1267.
  • Ewert M, Jackson DR, Nelson CE (1994) Patterns of Temperature-Dependent Sex Determination in Turtles. The Journal of Experimental Zoology 15: 3–15.
  • Filz KJ, Bohr A, Lötters S (2017) Abandoned Foreigners: is the stage set for exotic pet reptiles to invade Central Europe? Biodiversity and Conservation 27: 1–19.
  • Frazer NB, Whitfield Gibbons J, Greene JL (1990) Life Tables of a Slider Turtle Population. In Whitfield Gibbons, J (Ed.) Life History and Ecology of the Slider Turtle. Smithsonian Institution Press, Washington, 183–200.
  • Fujisaki I, Hart KM, Mazzotti FJ, Rice KG, Snow S, Rochford M (2010) Risk assessment of potential invasiveness of exotic reptiles imported to south Florida. Biological Invasions 12: 2585–2596.
  • García-Díaz P, Ross JV, Ayres C, Cassey P (2015) Understanding the biological invasion risk posed by the global wildlife trade: Propagule pressure drives the introduction and establishment of Nearctic turtles. Global Change Biology 21: 1078–1091.
  • García-Díaz P, Ross JV, Ayres C, Cassey P (2014) Data from: Understanding the biological invasion risk posed by the global wildlife trade: propagule pressure drives the introduction and establishment of Nearctic turtles. Dryad Digital Repository.
  • Gertzen E, Familiar O, Leung B (2008) Quantifying invasion pathways: fish introductions from the aquarium trade. Canadian Journal of Fisheries and Aquatic Sciences 65: 1265–1273.
  • Gollasch S, Lenz J, Dammer M, Andres H-G (2000) Survival of tropical ballast water organisms during a cruise from the Indian Ocean to the North Sea. Journal of Plankton Research 22: 923–937.
  • Grabowski M, Bacela K, Konopacka A (2007) How to be an invasive gammarid (Amphipoda: Gammaroidea) – Comparison of life history traits. Hydrobiologia 590: 75–84.
  • Juliano SA (2001) Nonlinear curve fitting: predation and functional response curves. In: Scheiner SM, Gurevitch J (Eds) Design and analysis of ecological experiments. Oxford University Press, Oxford, 178–196.
  • Kitowski I, Pachol D (2009) Monitoring the trade turnover of Red-eared terrapins (Trachemys scripta elegans) in pet shops of the Lublin region, East Poland. North-Western Journal of Zoology 5: 34–39.
  • Kulhanek SA, Ricciardi A, Leung B (2011) Is invasion history a useful tool for predicting the impacts of the world’s worst aquatic invasive species? Ecological Applications 21: 189–202.
  • Lamb RJ, MacKay PA, Migui SM (2009) Measuring the performance of aphids: fecundity versus biomass. The Canadian Entomologist 141: 401–405.
  • Laverty C, Dick JTA, Alexander ME, Lucy FE (2015) Differential ecological impacts of invader and native predatory freshwater amphipods under environmental change are revealed by comparative functional responses. Biological Invasions 17: 1761–1770.
  • Laverty C, Brenner D, McIlwaine C, Lennon JJ, Dick JTA, Lucy FE, Christian K A (2017a) Temperature rise and parasitic infection interact to increase the impact of an invasive species. International Journal for Parasitology 47: 291–296.
  • Laverty C, Green KD, Dick JTA, Barrios-O’Neill D, Mensink PJ, Médoc V, Spataro T, Caffrey JM, Lucy FE, Boets P, Britton JR, Pegg J, Gallagher C (2017b) Assessing the ecological impacts of invasive species based on their functional responses and abundances. Biological Invasions 19: 165–1665.
  • Nori J, Tessarolo G, Ficetola GF, Loyola R, Di Cola V, Leynaud G (2017) Buying environmental problems: The invasive potential of imported freshwater turtles in Argentina. Aquatic Conservation: Marine and Freshwater Ecosystems 27: 685–691.
  • Pearson SH, Avery HW, Spotila JR (2015) Juvenile invasive red-eared slider turtles negatively impact the growth of native turtles: Implications for global freshwater turtle populations. Biological Conservation 186: 115–121.
  • Piria M, Copp G, Dick JTA, Duplić A, Groom Q, Jelić D, Lucy FE, Roy HE, Sarat E, Simonović P, Tomljanović T, Tricarico E, Weinlander M, Adámek Z, Bedolfe S, Coughlan NE, Davis E, Dobrzycka-Krahel A, Grgić Z, Kırankaya SG, Ekmekçi FG, Lajtner J, Lukas JAY, Koutsikos N, Mennen GJ, Mitić B, Pastorino P, Ruokonen TJ, Skóra ME, Smith ERC, Šprem N, Tarkan AS, Treer T, Vardakas L, Vehanen T, Vilizzi L, Zanella D, Groom Q (2017) Tackling invasive alien species in Europe II : threats and opportunities until 2020. Management of Biological Invasions 8: 273–286.
  • Pöckl M (2007) Strategies of a successful new invader in European fresh waters: Fecundity and reproductive potential of the Ponto-Caspian amphipod Dikerogammarus villosus in the Austrian Danube, compared with the indigenous Gammarus fossarum and G. roeseli. Freshwater Biology 52: 50–63.
  • Polo-Cavia N, López P, Martín J (2011) Aggressive interactions during feeding between native and invasive freshwater turtles. Biological Invasions 13: 1387–1396.
  • R Core Team (2016) R: A Language and Environment for Statistical Computing. R Foundation for Statistical Computing, Vienna.
  • Raghavan R, Dahanukar N, Tlusty MF, Rhyne AL, Krishna Kumar K, Molur S, Rosser AM (2013) Uncovering an obscure trade: Threatened freshwater fishes and the aquarium pet markets. Biological Conservation 164: 158–169.
  • Rayner MJ, Hauber ME, Imber MJ, Stamp RK, Clout MN (2007) Spatial heterogeneity of mesopredator release within an ocean island system. Proceedings of the National Academy of Sciences of the United States of America 14: 20862–20865.
  • Ricciardi A, Hoopes MF, Marchetti MP, Lockwood JL (2013) Progress toward understanding the ecological impacts of non-native species. Ecological Applications 83: 263–282.
  • Rixon CAM, Duggan IC, Bergeron NMN, Ricciardi A, MacIsaac HJ (2005) Invasion risks posed by the aquarium trade and live fish markets on the Laurentian Great Lakes. Biodiversity and Conservation 14: 1365–1381.
  • Rödder D, Schmidtlein S, Veith M, Lötters S (2009) Alien invasive slider turtle in unpredicted habitat: A matter of niche shift or of predictors studied? PLoS ONE 4: E7843.
  • Seebens H, Blackburn TM, Dyer E, Genovesi P, Hulme PE, Jeschke JM, Pagad S, Pyšek P, van Kleunen M, Winter M, Ansong M, Arianoutsou M, Bacher S, Blasius B, Brockerhoff EG, Brundu G, Capinha C, Causton CE, Celesti-Grapow L, Dawson W, Dullinger S, Economo EP, Fuentes N, Guénard B, Jäger H, Kartesz J, Kenis M, Kühn I, Lenzner B, Liebhold AM, Mosena A, Moser D, Nentwig W, Nishino M, Pearman D, Pergl J, Rabitsch W, Rojas-Sandoval J, Roques A, Rorke S, Rossinelli S, Roy HE, Scalera R, Schindler S, Štajerová K, Tokarska-Guzik B, Walker K, Ward DF, Yamanaka T, Essl F (2018) Global rise in emerging alien species results from increased accessibility of new source pools. Proceedings of the National Academy of Sciences 115: E2264–E2273.
  • Standfuss B, Lipovšek G, Fritz U, Vamberger M (2016) Threat or fiction: is the pond slider (Trachemys scripta) really invasive in Central Europe? A case study from Slovenia. Conservation Genetics 17: 557–563.
  • Teillac-Deschamps P, Lorrillière R, Servais V, Delmas V, Cadi A, Prévot-Julliard A (2009) Management strategies in urban green spaces: Models based on an introduced exotic pet turtle. Biological Conservation 142: 2258–2269.
  • Tilman D, Clark M, Williams DR, Kimmel K, Polasky S, Packer C (2017) Future threats to biodiversity and pathways to their prevention. Nature 546: 73–81.
  • Tingley R, Mahoney PJ, Durso AM, Tallian AG, Morán-Ordóñez A, Beard KH (2016) Threatened and invasive reptiles are not two sides of the same coin. Global Ecology and Biogeography 25: 1050–1060.
  • Wasserman RJ, Alexander ME, Weyl OLF, Barrios-O’Neill D, Froneman PW, Dalu T (2016) Emergent effects of structural complexity and temperature on predator-prey interactions. Ecosphere 7: E01239.
  • Xu M, Dick JTA, Ricciardi A, Fang M, Zhang C, Gu D, Mu X, Luo D, Wei H, Hu Y (2016) Warming mediates the relationship between plant nutritional properties and herbivore functional responses. Ecology and Evolution 6: 8777–8784.
  • Zhao B, Chen Y, Lu HL, Zeng ZG, Du WG (2015) Latitudinal differences in temperature effects on the embryonic development and hatchling phenotypes of the Asian yellow pond turtle, Mauremys mutica. Biological Journal of the Linnean Society 114: 35–43.
  • Zi J, Pan X, MacIsaac HJ, Yang J, Xu R, Chen S, Chang X (2018) Cyanobacteria blooms induce embryonic heart failure in an endangered fish species. Aquatic Toxicology 194: 78–85.