Introduction

Invasive species threaten biodiversity through habitat alteration, increased predation pressure, and competition with native species (Wilcove et al 1998; Sala et al 2000; Gaertner et al 2009; Vilà et al 2011; Early et al 2016; Dueñas et al 2018). Range expansions of invasive wetland plants, such as Lythrum salicaria, Typha × glauca, Myriophyllum spicatum, Phalaris arundinacea, and Phragmites australis australis have dramatically modified many Nearctic wetlands (Galatowitsch et al 1999). By outcompeting native flora, invasive wetland macrophytes can influence plant species richness and composition, but they may also affect trophic interactions and habitat quality for wetland faunas (Chambers et al 1999; Keller 2000; Levine et al 2003; Minchinton and Bertness 2003; Minchinton et al 2006). Understanding whether endangered, native faunas can persist in wetland habitats that are modified by invasive plants is essential for accurately predicting the future impacts of invasive plants, and for prioritizing control actions. This requires a mechanistic understanding of how invasive species affect the behaviours of the native species that encounter them (Stewart et al 2021).

Palearctic Phragmites australis australis is a robust graminoid that has been spreading through Nearctic wetlands since its introductions to eastern North America in the late 1800s (Van Der Putten 1997; Chambers et al 1999; Orson 1999; Tulbure et al 2007; Vilà et al 2010; Hirtreiter and Potts 2012; Kiviat 2013). A native, non-invasive subspecies, Phragmites australis americanus, also occurs in Nearctic wetlands (Orson 1999; Saltonstall et al 2004), but the invasive P. a. australis (hereafter “Phragmites”) is now the dominant lineage in many areas of North America including regions surrounding the Laurentian Great Lakes (LGL) (Saltonstall 2002; Tulbure and Johnston 2010; Kirk et al 2011b; Meyerson and Cronin 2013). In wetlands where it is well established, Phragmites affects habitat structure and connectivity by altering vegetation composition, hydrology, canopy structure, and microclimates (Weinstein and Balletro 1999). Phragmites stands can be extremely dense, with up to 100 live culms per square meter and heights of up to four meters; these dense stands can maintain relatively cool temperatures close to soil and water surfaces in wetlands (Windham and Meyerson 2003; Rogalski and Skelly 2012; Mifsud 2014).

Across much of its introduced range, invasive Phragmites is broadly sympatric with Typha, which in Canada and the northern USA comprises native Typha latifolia L., introduced T. angustifolia L., and their hybrid T. × glauca (Godr.) (Galatowitsch et al 1999; Ciotir et al 2013; Ciotir and Freeland 2016). In regions surrounding the LGL and St. Lawrence Seaway, T. × glauca is considered invasive because it can outcompete its progenitor species (Frieswyk and Zedler 2007; Larkin et al 2012; Freeland et al 2013), and often forms dense stands that can significantly modify abiotic and biotic habitat characteristics such as water depth, nutrient concentration, species abundance and species richness (e.g. Tulbure et al. 2007; Tuchman et al. 2009; Bellavance and Brisson 2010; Lishawa et al. 2010; Mitchell et al. 2011; Hirtreiter and Potts 2012). In recent decades T. × glauca has invaded wetlands in the LGL region (Farrer and Goldberg 2009; Tuchman et al 2009; Lishawa et al 2010; Mitchell et al 2011; Geddes et al 2014; Bansal et al 2019), and is now the dominant Typha in those areas even though it often co-exists with one or both of its progenitor species (Travis et al 2010; Kirk et al 2011a; Freeland et al 2013; Pieper et al 2020). As a result, Typha stands around the LGL have been substantially modified over the past few decades, and we hereafter refer to stands that typically comprise mixed Typha species and hybrids as ‘Typha’.

Evidence for the impact of invasive Phragmites and Typha on Nearctic wetland faunas is equivocal (Bernstein and Mclean 1980; Kantrud 1986; Weis and Weis 2003; Hershner and Havens 2008; Martin and Blossey 2013; Gagnon Lupien et al 2015; Robichaud and Rooney 2017). Fowler’s toad (Anaxyrus fowleri) abundance declined following Phragmites invasion in a coastal wetland (Greenberg and Green 2013), but that invasion co-occurred with other habitat alterations that may also have affected toad abundance. Abundance and diversity of invertebrates and waterfowl is lower in monodominant stands of Typha compared to wetland habitats with higher heterogeneity (Kantrud 1986; Linz et al 1996b, 1999; Lishawa et al 2010), and abundance and diversity of birds are lower in Phragmites stands compared to native marsh habitats (Kantrud 1986; Benoit and Askins 1999; Robichaud and Rooney 2017). In contrast, dense Typha stands provide suitable roosting habitats for some bird species including red-winged blackbirds (Agelaius phoeniceus), yellow-headed blackbirds (Xanthocephalus xanthocephalus) and marsh wrens (Cistothorus palustris) (Linz et al 1996a; Linz and Homan 2011). Typha also provides food and cover for some mammals (e.g. muskrats, Ondatra zibethicus) (Boyce 1978; Grace and Harrison 1986; Campbell and MacArthur 1994), although it is unclear whether this is true of both native and invasive Typha. Invasive Phragmites stands provide habitat for native species such as mummichog (Fundulus heteroclitus; Fell et al. 1998), other freshwater fish (Wynia 2019), and American bullfrog (Lithobates catesbeianus; Rogalski and Skelly 2012).

As the ranges of invasive Phragmites and T. × glauca expand, understanding their effects on endangered turtles is crucial for mitigating the inadvertent effects of control measures for invasive macrophytes. Habitat restoration efforts in Canada and the United States currently control invasive Phragmites with measures that include herbicide applications (typically glyphosate-based products), mechanical rolling of stands, cutting and removal, and/or burning (Kettenring and Adams 2011; Wagner et al 2017; Quirion et al 2018). Typha × glauca is controlled in the United States with some of the same methods (Linz and Homan 2011; Lawrence et al 2016). If habitats dominated by invasive Phragmites or Typha are not suitable for endangered turtles, then controlling these macrophytes may be a priority for recovery of endangered turtle populations. For example, the range of invasive Phragmites is expected to completely overlap with the Canadian ranges of Blanding’s turtle (Emydoidea blandingii) and spotted turtle (Clemmys guttata) by 2030 (Catling and Mitrow 2011; Environment and Climate Change Canada 2018a; Environment and Climate Change Canada 2018b), and both species are listed as globally endangered by the International Union for the Conservation of Nature (van Dijk 2011; van Dijk and Rhodin 2011). Conversely, if turtles do not avoid these invasive-dominated habitats, then careful planning of macrophyte control is essential to minimize risk to resident turtles, while slowing invasion by non-native plant species.

There has been relatively little investigation into whether turtles actively avoid stands of large, invasive macrophytes, although tracking of Blanding’s turtles in two coastal wetlands revealed dynamic home ranges (i.e. home ranges can shift among years), and suggested that turtles did not avoid Phragmites when selecting their home ranges (i.e. second-order habitat selection, sensu Johnson 1980; Markle and Chow-Fraser 2018). However, at the individual level (third-order selection; Johnson 1980), males and non-gravid females avoided untreated Phragmites stands within their home range (Markle and Chow-Fraser 2018). Turtles did not avoid Typha, despite structural similarities between Typha and Phragmites stands (Farnsworth and Meyerson 2003; Paradis et al 2014).

The high shoot densities and substantial litter accumulation in Phragmites and T. × glauca stands (Tulbure et al 2007; Vaccaro et al 2009; Holdredge and Bertness 2011; Larkin et al 2012) could restrict movement of turtles. These dense stands could also reduce access to essential resources such as basking sites, food, and protective cover, and even affect nest success by invading and shading previously unshaded nesting sites (Bolton and Brooks 2010; COSEWIC 2016a; Cook et al 2018; Environment and Climate Change Canada 2018b). Markle and Chow-Fraser (2018) hypothesized that the individual-level avoidance of Phragmites by E. blandingii may have arisen because Phragmites either obstructs the movements of turtles, or provides inadequate thermal conditions.

The objective of our study was to directly test these two hypotheses: that stands of invasive macrophytes (Phragmites and Typha spp.) are barriers to movement and provide inadequate thermal conditions for turtles. Assumed support for these hypotheses is implicit in recent conservation status assessments for E. blandingii and the sympatric, highly endangered C. guttata (COSEWIC 2014, 2016b). Clemmys guttata is much smaller than E. blandingii, and we speculated that it might therefore be better able to navigate dense macrophyte stands. We used GPS tracking to collect hourly or bi-hourly relocations from E. blandingii and C. guttata across ~ 70 days in the peak active season of the two species. We used these high-resolution datasets to quantify fine-scale movement and habitat use of each species in a wetland mosaic that included stands of invasive Phragmites and Typha. We predicted that if Phragmites and Typha stands were barriers to movement, turtle movements through Phragmites or Typha would be slower than movements through other habitat types, and turtle movements through invasive macrophyte stands would have higher relative turning angles compared to movements in other habitat types (i.e. if macrophytes are barriers then they should impede navigation). We predicted that if the stands act as barriers or provide inadequate thermal conditions, then turtles moving through their home ranges would avoid habitats dominated by these macrophytes when selecting locations within a home-range, and would select microhabitats with lower densities of Phragmites and Typha shoots than they would encounter if moving randomly. Finally, we compared shell temperatures of turtles that were using Phragmites and Typha stands to those of turtles using uninvaded habitats in order to test the hypothesis that Phragmites and Typha stands interfere with a turtle’s ability to thermoregulate due to increased shading. Our study fills critical knowledge gaps to inform adaptive management of endangered turtle populations during biological invasions of wetlands.

Materials and methods

Characterization of available habitat

We conducted this study in a coastal wetland (~ 8.31 km2) on the north shore of Lake Erie, Ontario, Canada, one of the Laurentian Great Lakes (“Site A” in Markle and Chow-Fraser 2018). The wetland is a mosaic of habitat dominated by Typha and Phragmites, interspersed among open meadow marsh, forested areas, and a coastal dune system. To confirm that Typha stands at this site include the invasive hybrid T. × glauca, we collected 45 pollen samples from dehiscent staminate spikes at six locations across the study site on 28 June and 3 July 2019. We examined pollen grains under a compound microscope at 40 × magnification to identify plants as either T. latifolia (tetrad pollen grains; Online Resource 1; Fig. S1), T. angustifolia (monad pollen grains; Online Resource 1; Fig. S1) or T. × glauca (mixture of pollen grain types, including tetrads and monads; Online Resource 1; Fig. S1) (Krattinger 1975).

We used WorldView 2 multispectral satellite imagery (1.8 m spatial resolution, 8 multispectral bands) acquired on 5th September 2018 to map the habitat types using ENVI 5.5 (Exelis Visual Information Solutions, Boulder, CO). We used maximum likelihood classification to map seven habitat types: meadow marsh, Typha-dominated marsh, forest, open land, open water, Phragmites-dominated marsh, and shrub (see Table 1 for habitat descriptions). There have been efforts to control Phragmites at the study site, and the classification of Phragmites-dominated habitats included areas where Phragmites had been sprayed in previous years (but not rolled or burned). For the image classification and accuracy assessment, we used ground reference locations from vegetation surveys conducted in July and August 2018 (further details provided below), and from visual interpretations of high-resolution satellite imagery (pansharpened WorldView 2; 50 cm). We used only 192 out of 810 sampling locations collected in the field for this purpose, because these were the locations that had one land cover class covering over 80% of a 4 m2 vegetation survey quadrat (details of vegetation surveys are provided below). For the image classification, we used 100 locations from the field data and 98 locations from image interpretation. We used 92 of the ground-reference locations plus 233 locations selected from high-resolution image interpretation as test locations to assess the accuracy of the image classification. The test locations gave an overall accuracy of 84.27% for the image classification which allowed us to proceed with further analyses.

Table 1 Habitat classification within the minimum convex polygon of tracked Emydoidea blandingii (Blanding’s turtles; n = 14) and Clemmys guttata (spotted turtles; n = 12) in a coastal wetland in the Laurentian Great Lakes

Telemetry

Turtles were captured between 2 May and 21 June 2018 either using baited hoop and minnow traps (Memphis Net and Twine Company, Memphis, Tennessee, USA), or hand-capturing along line transects. Trapping sites and transects targeted all available habitat types that contained standing water, including stands of Phragmites and Typha, marsh meadow, and open water areas.

We affixed dual, very high frequency, global positioning system (VHF–GPS) transmitters (Pinpoint VHF 120 and 240; Lotek Inc., Newmarket, Ontario, Canada) to the rear costal scutes of turtles using epoxy (J-B Weld WaterWeld, J-B Weld Company, Texas, USA; Fig. 1). Turtles weighing < 180 g were outfitted with a 5 g Pinpoint120 (GPS location taken and archived every 60 min when the transmitter antenna was exposed). Heavier turtles (≥ 180 g) were outfitted with a 7 g Pinpoint 240 (GPS locations archived every 30 min when the transmitter antenna was exposed). Transmitters archived a temperature measurement alongside each GPS location, and at scheduled times (every 30 or 60 min) even where the GPS location could not be taken (i.e. the turtle was underwater or under cover at the scheduled location time). Total attachment weight for E. blandingii (transmitter plus attachment and epoxy) weighed 0.6 – 5% of the turtles’ body mass (mean: 2% of body mass). Total attachment weight for C. guttata ranged from 4—6% of the turtles’ body weight (mean: 5% of body weight).

We located turtles every 1–3 weeks, using a 3-element Yagi antenna (Lotek Inc.) and a Communication Specialists receiver (Communication Specialists Inc., Orange, California, USA) to download archived GPS data. We also replaced the transmitter on located turtles with a fully charged unit. We tracked each turtle from its capture date until its transmitter stopped functioning.

We fitted 27 E. blandingii and 17 C. guttata with transmitters. Unfortunately, the transmitters stopped functioning earlier than predicted by the manufacturer, precluding us from continuing the study into hibernation as initially planned, and from using data from all 44 turtles that were fitted with transmitters. We collected movement data between 2 May and 12 July 2018. Searches of wetland habitat within turtles’ home ranges were used to recapture turtles after their tags ceased functioning, and to remove non-functional transmitters.

We estimated home range kernels for 14 E. blandingii (7 males and 7 females; 15–416 relocations/individual; mean ± SE: 95 ± 35 relocations/turtle), and 12 C. guttata (7 males and 5 females; 38–516 relocations/individual; mean ± SE: 191 ± 43 relocations/turtle). Data spanned a range of 11–59 days/individual for E. blandingii (1328 total relocations), and 9–58 days/individual for C. guttata (2295 total relocations). The mean time difference between relocations was 4.5 h (range 0.5–336 h).

Habitat selection

We analyzed habitat selection using the package ‘adehabitatHR’ (Calenge 2006) in R 3.3.1 (R Core Team 2018) and ArcMap 10.6.1 (ESRI 2018). We explored habitat selection independently for E. blandingii and C. guttata and considered two spatial scales. Landscape scale habitat selection (2nd order selection, Johnson 1980) described selection of home ranges within the available habitat (Aebisher et al 1993; Rasmussen and Litzgus 2010). At this scale, we quantified home ranges (Aebisher et al 1993) by creating home range kernels for each turtle that had ≥ 15 relocations (GPS locations; n = 14 E. blandingii and 12 C. guttata). We used the method of Row and Blouin-Demers (2006) to select an appropriate smoothing function for the kernels. We delineated minimum available habitat for each species at the study site by creating a minimum convex polygon around all the turtle relocations obtained from GPS-tracking locations. Next, we simulated comparable home ranges that reflected available habitat. We simulated 20 “available habitat kernels” per turtle in ArcMap 10.6.1 (ESRI 2018) by creating polygons of the same size and shape as the actual home range kernels previously calculated for each turtle. We used the ‘Create random points’ tool in ArcMap 10.6.1 to locate centroids of available habitat kernels randomly in the minimum available habitat. We calculated the proportions of different habitat types in the actual home range kernels and the simulated, available habitat kernels using the dissolve tool and field calculator in ArcMap 10.6.1 (ESRI 2018).

We also examined home-range-scale habitat selection (3rd order selection; Johnson 1980). At this scale, we defined used habitats as those habitats in which turtles were located within individual home range kernels (Aebisher et al 1993; Rasmussen and Litzgus 2010). We defined available habitats as the habitats within home range kernels of each turtle (Aebisher et al 1993). We compared the proportion of relocations of each individual of each species in each used habitat (number of relocations in each habitat / total number of relocations in that turtle’s home range kernel) with the availability of each habitat (proportion of the turtle’s home range kernel represented by each habitat type).

We calculated Manly’s selection ratios (Manly et al 2002; Calenge and Dufour 2006) for each individual to describe the use of each habitat type in the home range relative to its availability in the minimum available habitat (landscape scale selection). Selection ratios (± confidence intervals, CI) between 0 and 1 indicate a habitat type used less than its availability (avoidance), while selection ratios (± CI) > 1 indicate habitat preference. Habitat types with a selection ratio (± CI) including 1 are neither avoided nor preferred (Calenge and Dufour 2006). We used Manly’s selection ratios to quantify home-range-scale habitat selection by comparing used locations for each individual to habitats available within its home range.

To ensure our conclusions were not a consequence of our choice of analyses, we also compared used and available habitats using compositional analyses (Aebisher et al 1993) and Eigen analyses of selection ratios (‘adehabitatHS’; Calenge and Dufour 2006), summarized in the Online Resource 1.

Microhabitat selection

To investigate the impact of macrophyte shoot density on microhabitat selection, we compared Phragmites and Typha shoot densities at locations used by each turtle to shoot densities at randomly selected locations within its home range. We used the R packages ‘sp’ (Pebesma and Bivand 2005) and ‘tripack’ (Renka et al 2020) to randomly select 10–15 GPS locations for each turtle, with a horizontal dilution of precision < 2. For turtles with more than 15 GPS points, we randomly sampled 15 GPS locations and 15 random locations evenly distributed across an individual home range. For turtles with fewer than 15 GPS points, we used every location and randomly sampled an equal number of random locations evenly distributed across an individual home range. We counted the number of Phragmites and Typha stems and recorded % cover of other vegetation types within a 4 m2 quadrat centered at each used and random location.

We compared the densities of Phragmites and Typha shoots in used and random locations with binomial generalized linear mixed models (GLMM, ‘lme4’ package, Bates et al. 2015). We ran an independent model for each species. Each model included a binary response variable (used or random location), a random factor (turtle ID, to account for the underlying variance in vegetation between different turtle home ranges), and fixed main effects of Phragmites density, Typha density, and their interaction.

Movement

Turtle movement speed could differ among marsh habitat types if it is physically more challenging for turtles to move through dense vegetation, or if turtle behaviour differs among habitats. If turtles avoid particular habitats, then we might also expect different trajectories in some habitats compared to others. We built trajectories for each turtle by connecting consecutive locations to quantify “steps”. We used GLMMs to test whether turtles’ movement speeds (m/hr) and turning angles (in radians, relative to previous step) differed in Phragmites or Typha patches compared to other habitats. We used separate models for each species and response (step speed and turning angle). We used only locations collected between 06:00 and 21:00, because nocturnal step lengths were biased towards zero. In our analyses for each turtle, we included only those steps for which the difference in time between locations was < 4 h. This approach excluded steps separated by long periods with no transmitter GPS reception, because we could not differentiate between scenarios in which these periods indicated time spent stationary under cover, or movement under water during which GPS points could not be recorded. We used random intercepts for turtle ID to model inter-individual behavioural variation, and to account for multiple observations per turtle. We compared turning angles and speeds between habitats using Tukey corrected pair-wise comparisons.

Temperature

We tested the hypothesis that Phragmites and Typha limit turtles’ diurnal access to thermal resources (i.e. exposure to sunshine) due to increased shading, which predicts that turtles will have lower diurnal shell temperatures on average when using Typha or Phragmites, compared to other habitat types. We compared shell temperatures of turtles in different habitats using temperatures recorded by transmitters at each daytime GPS point (hourly or bi-hourly locations collected from 07:00 to 20:00). We used a generalized additive mixed-effects model to compare turtle shell temperature between habitat types (fixed effect) at different times (smoothed effect of time). We included random intercepts for turtle ID and date to account for repeated measures of each turtle and the changing availability of temperature between days. We used separate models for E. blandingii and C. guttata.

Results

Characterization of available habitat

Examination of 45 pollen samples from Typha-dominated habitats used by turtles identified all three Typha taxa (Online Resource 1; Table S1; Fig. S1). Although the numbers of Typha precluded an exhaustive investigation into the proportions of each taxon, these data allow us to conclude that, as anticipated (e.g. Kirk et al. 2011a; Freeland et al. 2013; Pieper et al. 2020), the Typha stands in this wetland comprise multiple taxa including the invasive hybrid T. × glauca. The mean habitat patch size was 322 ± 34 m2 for habitat available to E. blandingii and 290 ± 30 m2 for habitat available to C. guttata. Dense Phragmites patches (i.e. > 70% cover) represented ~ 12% and 11% of the minimum available habitat of E. blandingii and C. guttata, respectively (Table 1). Dense Typha stands (> 70% cover) represented ~ 21% and 30% of the minimum available habitat of E. blandingii and C. guttata, respectively (Table 1).

Habitat selection

At the landscape scale, E. blandingii did not exhibit preference for, or avoidance of, any of the seven habitat types available to them when selecting home ranges from the minimum available habitat. (Fig. 2a). At the home range scale, E. blandingii avoided open land, forest and open water (CI for selection ratios < 1.0; Fig. 2b). They neither avoided nor preferred Typha, Phragmites, shrub and meadow marshes.

Fig. 1
figure 1

An Emydoidea blandingii (Blanding’s turtle) (a) and a Clemmys guttata (spotted turtle) (b) with dual VHF-GPS transmitter attached to rear costal scutes

Clemmys guttata selected home ranges containing more meadow marsh and less open water than available in simulated, randomly placed home ranges (Fig. 2c). We found no evidence of preference for, or avoidance of, the other habitat types by C. guttata at the landscape scale. At the home range scale, C. guttata preferred meadow marsh, and avoided forest, open land and open water (Fig. 2d). C. guttata neither avoided nor preferred Phragmites, Typha, or shrub within their home ranges. Compositional and Eigenvector analyses for both species produced similar results (Online Resource 1; Table S2; Table S3; Fig. S2).

Fig. 2
figure 2

Selection ratios (± 95 % confidence intervals) of habitat selection for Emydoidea blandingii (Blanding’s turtles; n = 14) at the landscape scale (a) and the home range scale (b). Selection ratios (± 95 % confidence intervals) of habitat selection for Clemmys guttata (spotted turtles; n = 12) at the landscape scale (c) and at the home range scale (d). Overlap of selection ratios with the horizontal dashed line indicates no selection or avoidance.

Microhabitat selection

Phragmites density was similar in quadrats used by E. blandingii (mean ± SE; 6.75 ± 1.17 stems/m2) and random quadrats (6.02 ± 0.85 stems/m2; GLMM: χ2 = 0.18, p = 0.67; Fig. 3a). Typha density was similar in E. blandingii quadrats (9.68 ± 1.33 stems/m2) and random quadrats (11.30 ± 2.02 stems/m2; GLMM: χ2 = 0.38, p = 0.54; Fig. 3a). There was no evidence that an interaction between Phragmites and Typha density (χ2 = 0.67, p = 0.41) influenced plot choice by E. blandingii.

Fig. 3
figure 3

Shoot density (stems/m2) of Phragmites australis and Typha measured in 4m2 quadrats was similar in quadrats used by turtles, compared to randomly selected quadrats within the turtle’s home ranges, for (a) Emydoidea blandingii (Blanding’s turtles; n = 200 used, 200 random quadrats) and (b) Clemmys guttata (spotted turtles; n = 180 used, 180 random quadrats)

Phragmites density was similar in C. guttata quadrats (7.08 ± 1.40 stems/m2) and random quadrats (4.80 ± 0.92 stems/m2; GLMM: χ2 = 1.74, p = 0.19; Fig. 3b). Typha density was similar in C. guttata quadrats (15.7 ± 1.57 stems/m2) and random quadrats (16.1 ± 2.02 stems/m2; GLMM: χ2 = 0.001, p = 0.97; Fig. 3b). There was no evidence that an interaction between Phragmites and Typha density (χ2 = 0.75, p = 0.39) influenced plot choice by C. guttata.

Movement

Movement speeds and turning angles for E. blandingii and C. guttata were similar in Phragmites and Typha stands compared to other habitat types. Emydoidea blandingii in open land moved at higher speeds than in other habitats (F = 2.71, df = 6, 794, p = 0.01; Fig. 4a; all other pairwise comparisons p > 0.05). The turning angles of E. blandingii were similar between steps within Phragmites, Typha, and all other habitat types (F = 1.22, df = 6, 451, p = 0.30; Fig. 4b). Movement speed of C. guttata was similar in all habitat types (F = 1.02, df = 6, 1302, p = 0.41). Turning angles of C. guttata differed among habitat types (F = 2.65, df = 6, 611, p = 0.02; Fig. 4b). Pairwise comparisons of each pair of habitat types suggested that the turning angles of C. guttata were slightly lower in meadow marshes (mean = 1.92 radians, 95% CI: 1.84 – 2.00) than in Typha patches (mean = 2.12 radians, 95% CI: 1.99 – 2.24, pairwise-adjusted p = 0.08), but was not statistically significant following Tukey adjustment of p-values to correct for multiple pairwise comparisons. Clemmys guttata turning angles were similar between all other habitat types (all other pairwise comparisons: p > 0.20).

Fig. 4
figure 4

The mean and 95% confidence intervals for (a) movement speed and (b) turning angles for Emydoidea blandingii (Blanding’s turtles; n = 14) and Clemmys guttata (spotted turtles; n = 12) in seven habitat types in a southern Ontario wetland

Temperature

Turtles in Phragmites and Typha were not cooler, on average, than turtles in other habitats. There was an effect of habitat type on the shell temperatures of E. blandingii (F = 2.62, df = 6, 1073, p = 0.02, n = 1084 temperature readings; Fig. 5a, b) where turtles in the forest habitat were an estimated 2.9 °C cooler than turtles in Typha habitat. No other pairwise comparisons between habitats were significant (all other pairwise-adjusted p > 0.17). There was no effect of habitat type on shell temperatures of C. guttata (F = 1.19, df = 6, 1527, p = 0.31, n = 1541 temperature readings; Fig. 5c, d).

Fig. 5
figure 5

Turtle shell temperatures of Emydoidea blandingii (Blanding’s turtles; n = 1084 temperature readings) in seven habitat types (a) and comparing changing temperatures through the day (b). Turtle shell temperatures of Clemmys guttata (spotted turtles; n = 1541 temperature readings) in seven habitat types (c) and comparing changing temperatures through the day (d). Solid lines are model predictions for each habitat. Dashed lines are ± 1 SE

Discussion

Our results do not support the hypotheses that invasive Phragmites or mixed-species stands of Typha impede the movements or limit thermoregulation of threatened turtles. Previous year-round tracking of turtles allowed estimates of the turtles’ complete home ranges and seasonal shifts in habitat use, and inferred individual-level avoidance of Phragmites by E. blandingii within their home ranges (Markle and Chow-Fraser 2018). We tested two potential mechanisms driving that avoidance, using datasets with high temporal resolution over a short study period spanning the time of greatest activity for temperate freshwater turtles. Stands of Phragmites and Typha did not restrict the movements of mature, active E. blandingii or C. guttata, and use of these stands did not result in cooler shell temperatures than use of other wetland habitat types. We also found no evidence that density of Phragmites or Typha influenced turtle’s use of these habitats. Our results suggest that turtle avoidance of dense macrophytes is not the result of physical barriers or inadequate thermal conditions, and further research should test other potential explanations. From the perspective of managing biological invasions, our study demonstrates regular use of invasive Phragmites and Typha stands by endangered E. blandingii and C. guttata. This finding has implications for the management of invasive macrophytes in wetlands where endangered turtles are present, as control efforts should mitigate potential harm to turtles.

Our results provide a data-rich snapshot of active-season habitat use by C. guttata and E. blandingii. However, our study period captured only a part of the turtles’ yearly cycle, and habitat selection varies among seasons. We initially intended to track turtles to over-wintering sites and continue our study for an entire year, to capture seasonal variation in habitat selection. Unfortunately, the transmitters we used stopped working much earlier than we expected (i.e. stopped recharging and holding their charge). The active season of freshwater turtles is often divided into ‘behavioural’ seasons (i.e., emergence from hibernation, nesting, post-nesting, overwintering;Rasmussen and Litzgus 2010; Markle and Chow-Fraser 2014). Different biological and behavioural needs related to mating and thermoregulation are directly linked to specific habitat types used during different ‘behavioural’ seasons (Joyal et al 2001; Litzgus and Mousseau 2004; Beaudry et al 2009), and previous studies have stressed the importance of the temporal scale in assessing habitat selection by turtles. We have not attempted to generalize our results to other seasons, as we did not capture the full annual cycle of E. blandingii or C. guttata. Instead, we used the high frequency of relocations in our dataset to test two potential mechanisms driving active-season habitat selection in invaded habitats.

Autocorrelation between relocations is widespread in tracking and habitat selection studies (Fieberg et al 2010), but can be accounted for with appropriate analytical approaches. Autocorrelation between successive relocations is unlikely to have affected the conclusions from our habitat selection analyses because we recorded turtles moving up to 208 m within 0.5 h, which would allow turtles to move between habitat types based on the mean habitat patch size in our study area for each species. Our results demonstrate that invasive Phragmites and mixed-species stands of Typha can provide suitable summer habitat for E. blandingii and C. guttata when these habitats are embedded in a larger matrix of available habitat types.

We also acknowledge that turtle use of macrophyte stands may be influenced by differences among wetlands, including the distribution of macrophyte stands within the available turtle habitat. Habitat selection documented for a species at one or a few sites is often generalized across sites or regions during conservation planning (e.g. COSEWIC 2014, 2016b), yet our results differ slightly from those obtained by Markle and Chow-Fraser (2018) in the same wetland. Variation exists not only among wetlands, but also among years (for example, as shifting water levels alter the depth of coastal wetlands). Water levels at our study site were relatively low in 2011 and 2013, when turtles were tracked by Markle and Chow-Fraser (2018). In 2018 when we tracked turtles at the same site, the water level was at least 0.5 m higher, which may have made macrophyte stands more accessible, and ongoing erosion, invasion by macrophytes, and treatment of invasive macrophytes had changed the shape and habitat composition of the wetland. The effects of shifting water levels and erosion on habitat use by turtles in coastal wetlands deserves further study. Finally, while our data suggest that turtles at our study site use invasive macrophyte stands, at least during their peak active season, our study site was invaded over a decade ago. The responses of turtles to invasion (i.e. shifts in habitat use in response to rapidly changing habitat composition) deserves robust analysis, as also suggested by Markle and Chow-Fraser (2018).

Emydoidea blandingii are relatively large turtles, (adult carapace length between 150 and 240 mm; Tyning 1990; Congdon et al. 1993). In our study, they showed no evidence of avoiding Phragmites or Typha, nor a clear preference for any habitat types when they chose a home range or moved between locations within their home range. The smaller C. guttata (max. carapace length approx. 130 mm;Tyning 1990; Ernst et al. 1994; Litzgus 2006) also did not avoid Phragmites or Typha stands. However, C. guttata showed a preference for meadow marshes within their home ranges, as previously reported (Rasmussen and Litzgus 2010). Home range size varies widely among individual E. blandingii (between 12 and 60 ha in Canada; (COSEWIC 2016b), but E. blandingii home ranges are generally much larger than those of C. guttata (between 0.7 and 8.8 ha; COSEWIC 2014). Thus, E. blandingii may tolerate more habitat loss within their larger home ranges before losing access to a particular habitat type, while the smaller home ranges of C. guttata may leave them more vulnerable to loss of particular, key habitats (e.g. meadow marshes).

Clemmys guttata may prefer meadow marshes because these are structurally complex habitats which provide access to shallow water, a range of basking sites, and high densities of aquatic and emergent vegetation (Rasmussen and Litzgus 2010). Meadow marshes may also harbour relatively high concentrations of prey, such as juvenile amphibians (Fell et al 2003; Meyer 2003). Invasive Phragmites and T. × glauca both tend to replace meadow marsh habitats (Wilcox et al 2003; Frieswyk and Zedler 2007; Catling and Mitrow 2011; Markle and Chow-Fraser 2018), which could reduce overall habitat quality for C. guttata even if the turtles do not avoid the Phragmites and Typha stands themselves. Additional research is needed to determine the differences in resource availability and quality in invaded versus non-invaded habitats used by freshwater turtles with different spatial and biological requirements.

Given that invasive Phragmites and Typha do not provide a barrier to movement, or limit access to thermal resources, further research should test alternate hypotheses to explain the avoidance of these stands observed in previous research (Markle and Chow-Fraser 2018). It is possible that these stands provide reduced foraging opportunities compared to other habitat types (Mifsud 2014; Markle and Chow-Fraser 2018). Access to food resources can influence habitat selection by turtles (Rasmussen and Litzgus 2010), as turtles require sufficient food during the spring and summer to fuel their seasonal metabolic activities (e.g., growth, reproduction and storage), and prepare for brumation (Congdon 1989). However, Phragmites patches can support aquatic invertebrates, fish, and amphibians (Fell et al 1998, 2003; Angradi et al 2001; Meyer 2003), all of which are prey for E. blandingii and C. guttata (Lagler 1943; Kofron and Schreiber 1985; Rowe 1992; Rasmussen et al 2009).

Although Phragmites stands are often home to a myriad of prey suitable for turtles, additional investigations have concluded that these invasive plants are detrimental to turtles. For example, the replacement of shallow aquatic habitats and mixed meadow marshes by invasive Phragmites was considered a potential driver of extirpation of C. guttata at Point Pelee National Park, Ontario (Browne and Hecnar 2007; Markle et al 2018). The extent of Phragmites- and Typha-dominated habitats at our study site might not have reached a threshold that critically limits turtles’ access to meadow marshes or other key habitats. Our study was not designed to test whether there is a threshold proportion of meadow marsh habitat that C. guttata require within their home ranges, below which a population may not persist. Monitoring the trends and distributions of turtle populations in relation to the expansion of Phragmites and T. × glauca is necessary before we can determine whether the spread of invasive Phragmites and Typha correlates with the declines of freshwater turtles, and if so, whether there is a threshold beyond which Phragmites and T. × glauca expansion cause turtle populations to decline.

Although we found no evidence for negative impacts on turtle habitat in this study, invasive Phragmites and Typha do impact some taxa, and can modify wetland structure and function (Tulbure and Johnston 2010; Mitchell et al 2011; Lishawa et al 2014; Lawrence et al 2016; Robichaud and Rooney 2017; Yuckin and Rooney 2019). Coastal wetlands in the Laurentian Great Lakes host high biodiversity and provide a range of ecosystem services (Sierszen et al 2012), and control of invasive macrophytes may be required to protect these. Macrophytes such as Phragmites and T. × glauca can be controlled mechanically by rolling, mowing and cutting, and are often also treated with herbicides such as glyphosate or imazapyr (Mozdzer et al 2008; Hazelton et al 2014; Wagner et al 2017; Judd and Francoeur 2019). Wetlands invaded by Phragmites and Typha are often repeatedly colonized because populations of these taxa can expand through both seed and rhizome dispersal (Kirk et al 2011b; Pieper et al 2020). Management of invasive macrophytes is therefore an ongoing process.

One objective of controlling invasive Phragmites and Typha is to protect coastal wetland habitats for the native species that rely on them. Understanding how control measures might instead inadvertently harm native species can inform control efforts and thus reduce risk to endangered turtles and other species using invaded habitats. Mechanical control and some herbicide application methods require the use of machinery that is heavy enough to injure turtles or other non-target species, as has been observed at some control sites (Cross 2009; Zarnetske et al. 2010; Howey and Roosenburg 2013; S. Gillingwater, pers. comm.). Limiting the use of heavy machinery in macrophyte stands to the winter months may reduce the risk of encountering turtles. Over-wintering turtles require access to water with specific environmental characteristics that will not freeze solid (Litzgus et al 1999; Edge et al 2009), which macrophyte stands in shallow water may not provide, but the risk of encountering and potentially harming overwintering turtles cannot be completely eliminated. Control measures that are undertaken during the winter months should plan for potential encounters with overwintering turtles, and ideally, identify local overwintering sites ahead of time (i.e. through telemetry), so that these can be avoided. Finally, from the perspective of turtle habitat restoration in invaded areas, restoration of meadow marsh habitats is likely to benefit endangered species such as C. guttata.