Study Species
Black spiny-tailed iguanas
Black spiny-tailed iguanas (Ctenosaura similis) are native from southern Mexico to Panama but are considered invasive in Florida (Avery et al. 2014). The species was first introduced to Florida on Gasparilla Island in 1979 through the exotic pet trade and they have since spread across the island. They have become established throughout southwestern Florida and there are also now populations along the southeastern coast and the Florida Keys in residential, commercial, and natural areas (Krysko et al. 2003).
Juveniles mainly feed on invertebrates and their diet shifts towards herbivory with age. However, adults will also prey on invertebrates, other lizards, small birds, small mammals, and turtle hatchlings, including the gopher tortoise (Fitch and Hackforth-Jones 1983; Avery et al. 2009; Krysko et al. 2009; García-Rosales et al. 2020; Furness 2021). Within their native range black spiny-tailed iguanas are known to burrow into earthen banks, and under rocks and logs (Burger and Gochfield 1991). In Florida they burrow around sea walls and rock piles, under houses, and will also inhabit gopher tortoise burrows. They are diurnal and seek these shelters before nightfall, emerging again in the morning. Within tortoise burrows they dig side tunnels into the wall of the main burrow excavated by the tortoise (personal observation, SLM). They are large, highly fecund lizards exceeding 30 cm SVL with clutch sizes up to 62 eggs recorded on Gasparilla Island, FL (Avery et al. 2014).
Gopher Tortoises
Gopher tortoises (Gopherus polyphemus) occur throughout the Southeastern Coastal Plain of the U.S. from Louisiana east to South Carolina and south to Florida. They occur in dry, upland habitats and require well-drained, sandy soils for burrowing and nesting, abundant herbaceous plants for forage, and open canopies that provide adequate sunlight for basking and nesting (Hipes et al. 2000). They are terrestrial turtles, reaching carapace lengths of 31 cm (Ernst et al. 1994). They mature slowly, with females becoming reproductive at 9–21 years depending on geographic location and environmental conditions (Landers et al. 1982; Godley 1989; Mushinsky et al. 1994). They have shovel-like forelimbs adapted for digging long burrows in which they live. Adult burrows average 4.5 m long and 2 m deep. Burrows typically consist of one slightly curving tunnel with a single entrance (Hansen 1963). Soil from the burrow is deposited near the burrow mouth, forming a conspicuous mound (i.e., apron) that is visible even long after a tortoise has abandoned a burrow.
Gopher tortoises spend as much as 80% of their time in their burrows, which they excavate for shelter from extreme temperatures, desiccation, fire, and predators (Hipes et al. 2000). They are diurnal, entering their burrow before nightfall and emerge to forage, thermoregulate, and carry out other activities during the day. Gopher tortoises are active year-round in southern Florida but are least active, spending the most time in their burrows, from December-February. In southern Florida, surface activity occurs between 0700–2000 h during summer months, but only between 1100–1800 h on warmer days during December-February. Activity peaks occur during the hottest hours of afternoon (1300–1600 h) throughout the year (Douglass and Layne 1978).
The number of burrows used per individual tortoise over the course of a year ranges from 1–7 and varies by sex, age class, and population. Older tortoises use more burrows than younger ones and males use more than females (McRae et al. 1981; Diemer 1992; Wilson et al. 1994). Reasons for this variation may include any combination of habitat quality, soil composition, temperature extremes at differing latitudes, breeding behavior differences, and number of disturbances to burrows, such as that of potential predators (Wilson et al. 1994).
Gopher tortoises are ecosystem engineers because their burrows provide refugia, foraging, or reproduction sites to at least 360 other species, some of which are obligate commensals (Jackson and Milstrey 1989). They are considered a keystone species, as tortoise burrow density has been shown to be a strong predictor of diversity and relative abundance of herpetofauna and small mammals (Catano and Stout 2015). Gopher tortoises are federally listed as Threatened in the western portion of their range (U.S. Fish and Wildlife Service 1987), designated as a Threatened species in Florida, and are state protected on some level in every other part of their range (Florida Fish and Wildlife Conservation Commission 2020).
Study Sites
Florida boasts more established nonnative herpetofauna species than any other region in the world. As many as 180 introduced amphibian and reptile taxa have been documented within the state, and of those, at least 63 (dominated by reptiles) are considered established (Krysko et al. 2016). In fact, Florida has three times as many established species of introduced lizards as native lizards (Hardin 2007). The magnitude of herpetofauna invasions in Florida can be attributed to several factors. Florida is a shipping hub with major ports of entry for arrival of wildlife species (both legal and illegal), thus resulting in high propagule pressure, a process known to contribute to invasion success (Lockwood et al 2005). The state has a popular captive wildlife industry, and facilities housing nonnative species (e.g., breeders, retail vendors, individual pet owners) are susceptible to escapes and sometimes intentional releases of unwanted or “seeding” animals (Fujisaki et al. 2010; Episcopio-Sturgeon and Pienaar 2019). Additionally, Florida’s mild climate and abundant rainfall results in a strong climate match for many introduced amphibians and reptiles and climate is recognized an important predictor of invasion success (Bomford et al. 2008; Engeman et al. 2011).
We conducted our study on Gasparilla Island, an 11-km barrier island off the southwest coast of Florida (Fig. 1). The island includes six distinct natural communities: mangrove swamp, beach dune, coastal strand, maritime hammock, estuarine unconsolidated substrate, and marine unconsolidated substrate (Florida Natural Areas Inventory 2010). Most of the island is private property and highly developed, containing many homes, resorts, businesses, and a golf course. We sampled tortoise burrows at two study sites on Gasparilla Island within the beach dune habitat, which we designated as our Private Property (PP) and Gasparilla Island State Park (SP) sites (Fig. 1). At PP, USDA Wildlife Services has conducted ongoing iguana removal efforts, with over 15,000 iguanas removed from an area of about 106-ha since 2008 (Avery et al. 2014; Samantha Watson, personal communication, Feb. 3, 2021; Glenn Mitchell, personal communication, July 18, 2021). On the other hand, iguana management at the SP site has been limited due to a lack of resources and high visibility by park guests to iguana removal efforts. Only ~ 600 iguanas were removed from the 51-ha area since 2012, and there are no records of removals prior to then (Karen Rogers, personal communication Dec. 2, 2021).
Experimental Design
We delineated beach dune habitat within our study sites visually using the World Imagery layer in ArcGIS® Online by Esri (Source: Esri, Maxar, Earthstar Geographics, and the GIS User Community). Although gopher tortoises can occur in coastal strand, we excluded this habitat because of low tortoise burrow density and the difficultly of confidently locating all burrows due to dense vegetation (Lau and Dodd 2015). We generated 200 random points spaced at least 10 m apart within the beach dune at each study site and selected the burrow closest to each random point. We selected the closest burrow by walking to each random point using an ArcGIS® Survey123 map (Esri), staking a long PVC tube in the ground at the point, and walking outwards from there in a gradually expanding spiral until the first burrow was found within ~ 10 m of the point, marking each burrow with a metal ID tag and stake. We discarded points farther than ~ 10 m from a burrow, or where the closest burrow was already marked from another PVC point. We measured burrow widths by inserting calipers 50 cm into the burrow or as deep as possible to compare burrow size-class distributions between sites. Burrow width is strongly correlated with tortoise carapace length (CL) of the resident tortoise (Alford 1980; Martin and Layne 1987) and CL is associated with age in gopher tortoises (Landers et al. 1982). During pilot surveys, we noted iguanas and tortoises in burrows that appeared inactive (e.g., having an entrance filled with debris, no tracks on the burrow apron). Therefore, we considered tortoise burrows with an entrance open enough such that a tortoise could enter without modification as potentially occupied and included them in our surveys regardless of apparent activity status.
We determined the occupancy of tortoise burrows by tortoises and iguanas by inspecting them with a burrow camera (Environmental Management Systems, Canton, Georgia, USA). The camera kit consisted of a 7.6 m scope with a 5 cm diameter camera-head and a “juvenile” camera with a 3.6 m scope with a 2.5 cm camera-head. The camera feed was displayed on a color monitor housed in a watertight Pelican™ case. We scoped burrows from January to March 2022, between sunrise and noon to maximize the likelihood that tortoises and iguanas were inside the burrows. In an occupancy modeling framework, site occupancy is assumed to be constant between surveys, so repeat surveys were conducted back-to-back such that animals would not change locations between burrow surveys on a given day (i.e., we inserted the scope as far as we could down the burrow and any side-tunnels, pulled out the scope, and scoped it again). We treated each pair of surveys at a burrow as a “season” where occupancy in a burrow is assumed to be constant between the two consecutive surveys. Two surveys per season is considered appropriate when detection probability of a species is high (ρ > 0.6) (MacKenzie et al. 2018).
Iguanas frequently dug side-tunnels in the tortoise burrow at angles that were impossible to scope with a standard burrow scope. To help address this, we attached the string of a dive reel to the camera to make it more maneuverable. We positioned the head of the scope in front of a side tunnel and by pulling the string we could lift the head of the scope or turn it to the side to enter a side tunnel that otherwise could not be scoped due to the angle. We recorded the occupancy state of the tortoise burrow as neither (0), only iguana (1), only tortoise (2), or both (3) for each survey, and we recorded the presence of iguana side-tunnels. Iguanas frequently excavated side-tunnels off the main tortoise burrow (Fig. 2) and we fully scoped as many as possible. Most side-tunnels had obvious crest-shaped marks along their roof, and many had obvious iguana tail drags and/or claw marks on their floor.
We surveyed burrows at the PP site for three seasons with two surveys per season and at the SP site for one season, also with two surveys per season. The third season at PP and the single season at SP overlapped in time. Surveys at the PP site in Season 1 occurred from Jan. 3–18, Season 2 from Jan. 19 to Feb. 9, and Season 3 from Mar. 7–22. Burrow surveys for the single season at the SP site occurred Mar. 10–20.
Analysis
We used MSOMs to examine species co-occurrence patterns between black spiny-tailed iguanas and gopher tortoises to test our predictions relating to gopher tortoise burrow use given knowledge of their natural histories. We modelled co-occurrence patterns between the two species while accounting for imperfect detection (MacKenzie et al. 2018; Richmond et al. 2010; Waddle et al. 2010). We interpreted the results of our models using 95% confidence intervals. We included the detection data input and PRESENCE output for each model in Online Resource 2.
To address prediction 1, that tortoise occupancy would be lower given the presence of iguanas, we conducted a single-season two-species occupancy analysis in Program Presence (version 2.13.18) using data from the third season at the PP site and the one season at SP that overlapped temporally, estimating detection and occupancy probabilities of iguanas and detection and marginal occupancy probabilities of tortoises (Rota et al. 2016). This model allowed detection probability to be different for each species. We assumed that this detection probability was fixed and the same for both surveys of a burrow. We did not allow detection probability of tortoises to differ conditional on the presence or absence of iguanas because the presence of an iguana would not affect our ability to detect a tortoise with the burrow scope. We also conducted a model with the same parameterizations except with tortoise occupancy conditional upon the presence of iguanas. We compared the two models using ΔAIC to determine if the model including iguana occupancy as a predictor of tortoise occupancy was more parsimonious (see: “Conditional and marginal occupancy models” in Results.
To address prediction 2, that presence of iguana side-tunnels within tortoise burrows is positively associated with iguana occupancy and negatively associated with tortoise occupancy, we used these same data to conduct a single-season two-species analyses estimating the effects of binary covariate of iguana side-tunnel presence within burrows (present or absent, see “Iguana side-tunnel model” in Results. We did the same to address prediction 3, that iguana occupancy is lower and tortoise occupancy is higher in sites with ongoing iguana removal, but instead including the binary site covariate of site management (PP or SP, see “Private property management model” in Results). We also conducted a multi-season analysis using the three seasons of data from the PP, but the data were insufficient to provide clear results. We included this analysis in Online Resource 3.