Bahamas connection: residence areas selected by breeding female loggerheads tagged in Dry Tortugas National Park, USA

Abstract

Background

Delineation of home ranges, residence and foraging areas, and migration corridors is important for understanding the habitat needs for a given species. Recently, many population segments of Northwest Atlantic loggerhead sea turtles (Caretta caretta) were designated as endangered or threatened; the smallest subpopulation is in the Dry Tortugas. Foraging and residence areas for this subpopulation have not been defined outside the Gulf of Mexico. Here, for Dry Tortugas loggerheads that traveled to the Bahamas, we use a combination of switching state-space modeling (SSM) and home-range estimators to determine migration period, spatially delineate and describe residence areas, and examine inter-annual home-range repeatability.

Results

In 5,973 tracking days, migration dates for Dry Tortugas loggerheads traveling to the Bahamas occurred during July–September, with turtles tracked twice showing remarkably similar migration paths and timing of departure from nesting sites. Core-use residence areas for 19 loggerheads ranged from 3.7 to 179.5 km² (mean ± 1 SD = 56.2 ± 49.5 km²). For three turtles, we found inter-annual home-range repeatability, with centroids of core areas only 0.7–2.9 km apart and significant overlap of inter-annual 50% kernel contours.

Conclusions

We demonstrate a previously unknown link between Dry Tortugas nesting beaches and Bahamas residence areas; 17/39 (43.6%) of nesting loggerheads tagged in and tracked from the Dry Tortugas take up residence at sites in the Bahamas. Residence area estimates for these turtles were similar in size to previous foraging area estimates for two turtles tracked to the Bahamas in other studies. We show inter-annual residence area repeatability, and that residence areas of different individuals generally did not overlap. We suggest that these loggerheads possibly establish territories.

Background

Delineation of home ranges, residence and foraging areas, and migration corridors is essential to understand the spatial extent of habitats necessary for a given species. Defining these areas can be challenging in the marine environment given their often separate geographic locations, but electronic tags (i.e., satellite, GPS, and geolocation tags) have allowed delineation of high-use at-sea residence sites for marine megafauna, including sea turtles [1-6]. Hamann et al. [5] ranked efforts to define important in-water habitats as one of the top priorities worldwide for imperiled marine turtles.

As satellite-tracking studies have grown in number and scope over the past several decades [7-9], evidence is emerging that marine turtles occupy discrete habitats (e.g., residence and foraging areas) between nesting seasons, where individual turtles have been documented to reside for one year or more [9,10]. Connections between rookeries and foraging grounds are only now being quantified [2-4,11,12]. Loggerhead sea turtles nest every 2–5 years, depositing 2–6 clutches of 75–120 eggs approximately every 2 weeks during the nesting season [13-18]. After nesting, they migrate back to their residence or foraging site, to which they show a high level of fidelity [3,9,16,19-21]. For loggerhead sea turtles specifically, locations of oceanic ‘hotspots’ representing key residence and foraging habitats have recently been delineated in the southeast US [11,22] and the Gulf of Mexico (GoM) [3,12], yet the mechanisms for selection of these areas remain poorly understood. Additionally, while foraging likely takes place within areas designated as ‘foraging’ throughout sea turtle literature due to long-term site fidelity, the quantification of this behavior is lacking and the extent to which these areas are foraged remains largely unknown.

Loggerhead numbers are considerably reduced from historic estimates [23]. Imperiled loggerhead sea turtles in the Northwest Atlantic are listed as threatened under the U.S. Endangered Species Act. Consistent interactions with fisheries and nest number declines at major rookeries throughout Florida instigated a proposal to raise the level of protection for loggerheads in the USA [24,25]. Recently, nine distinct population segments of loggerheads were designated as endangered or threatened [26], with proposals following for critical habitat [27,28].

In the Northwest Atlantic, loggerheads exist as five subpopulations [29] and ten management units [30,31] based on analyses of mitochondrial DNA. The five subpopulations are: 1. Northern (southern Virginia to Florida/Georgia border); 2. Peninsular Florida (Florida/Georgia border through Pinellas County, Florida); 3. Dry Tortugas (DRTO; islands west of Key West, Florida); 4. Northern GoM (Franklin County, Florida through Texas); and 5. Greater Caribbean (all other nesting beaches throughout the Caribbean and Mexico) [30,32,33]. The smallest subpopulation is in DRTO, with an estimated nesting subpopulation of 258–496 females (50 percentile distribution = 331; [34]). The current status of the DRTO subpopulation is unknown, yet recent tracking studies are beginning to define the characteristics of this subpopulation. Through satellite tracking, Hart et al. [18] characterized inter-nesting areas for DRTO loggerheads (seven adult females) and Hart et al. [3] defined foraging areas in the GoM for four DRTO-nesting loggerheads.

Residence areas for the DRTO subpopulation have not previously been defined outside the GoM, although loggerhead flipper tag returns from other populations have indicated connections between the Bahamas and both Cuba [35] and Florida [36]. More recent satellite telemetry studies have corroborated these connections, demonstrating that some adult female loggerheads from three subpopulations (Northern, Peninsular Florida, Northern GoM) traveled to the Bahamas after nesting on various Florida beaches [37 (2 turtles), 23 (6 turtles), 14 and 12 (same 2 turtles), 38 (1 turtle), 13 (5 turtles)], and one adult male traveled there from a breeding area near Florida [37,38]. Despite 15 tracks/turtles tracked to the Bahamas, estimates of residence area size and characterizations of occupancy patterns at specific sites are limited (but see Additional file 1 showing previous estimates of loggerhead foraging area size in the Bahamas). Quantifying residence area size and verifying foraging behavior for this threatened species are especially important as the Bahamas are within the Caribbean Islands Biodiversity Hotspot which is further listed as one of the top eight ‘hottest hotspots’ due to the high number of endemic species and habitat loss [39].

Here, we use a combination of switching state-space modeling (SSM) and home-range estimators (i.e., kernel density estimation [KDE] and minimum convex polygons [MCP]) to determine periods of migration and residency and delineate the spatial extent of residence areas for post-nesting DRTO loggerheads that traveled to the Bahamas. SSM has been used to identify locations when turtles show directed movements versus area-restricted search (ARS) patterns—deemed previously as migration and inter-nesting or foraging ‘modes’ [1,3,4,40-50]. When combined with MCP (simple polygon created with home-range locations [51,52]) or KDE (a non-parametric method used to identify one or more areas of disproportionately heavy use [i.e. core areas] within a home-range boundary, see [53-55]), SSM has been suggested as a way to enhance analysis of Argos tracking data sets [56]. SSM was recently used to delineate foraging zones for post-nesting Kemp’s ridleys in the GoM [4], define the migrations and foraging areas for leatherback turtles nesting in the Pacific Ocean [47,57], determine foraging areas for juvenile green turtles in the SW Atlantic [58], and describe residence or foraging areas in the GoM for female loggerheads from three separate subpopulations [3].

Our specific objectives were to determine for DRTO loggerheads: (1) timing of post-nesting migrations; (2) location and spatial extent of residence areas; (3) characteristics of residence areas (i.e., bathymetry, distance from shore); and (4) inter-annual home-range repeatability at residence areas for a subset of turtles tracked in successive nesting seasons.

Results

Turtles

Upon first capture, turtles (n = 19, four of which were tagged twice) that we intercepted and tagged after nesting in DRTO on either East Key or Loggerhead Key (Figure 1) ranged in size from 86.0 to 111.6 cm straight carapace length (SCL; mean ± SD = 93.8 ± 5.7 cm; Table 1). Recaptures had remigration intervals of 3 years (turtle B) and 2 years (turtles E, G, L). In a total of 5,973 tracking days, individual turtle tracking durations ranged from 68 to 704 days (mean ± SD = 259.7 ± 194.0 days; Table 1).

Figure 1

Study site. Female loggerhead sea turtles (Caretta caretta) were tagged at Loggerhead Key and East Key in Dry Tortugas National Park (DRTO, shown as zoomed-in area within red box; thick blue line represents outer boundary of the park, thin blue line represents different management zones within the park). Turtles migrated from DRTO through Florida Straits to the Bahamas.

Table 1 Summary of satellite-tracking details for nesting loggerheads tracked to residence areas in the Bahamas 2008–2013

State-space modeling and migration time

The SSM predicted 408 total migration days across all 23 tracks (range 9–27 days; mean ± SD = 17.7 ± 5.6 days; Table 2) from DRTO to the Bahamas (approximately 400–700+ km traveled). Migration dates ranged from 7 July to 2 September with the exception of one turtle (turtle D) for which the SSM predicted additional migration dates in late September (16–22 September, see Figure 2). For all years combined, by the 2nd week of July, approximately 50% of the tracked loggerheads were in migration mode on their way to the Bahamas (Figure 3). Turtles traveled through the Florida Straits to their residence areas; most followed the Straits around the south/southeast of Florida before veering south from there (Figure 4). The SSM estimated a more southern course for turtles D and F, with turtle F’s track following the coast of Cuba; however, it should be noted that for both of these turtles, we received no locations between the Florida Keys and the Bahamas. For turtles that were tracked twice (n = 4), their migration dates overlapped in all cases (Figure 2). Turtle B had the same dates of migration for 8 of 15 total migration days (total = first month/day of earliest migration to last month/day of latest migration); turtle E had the same dates of migration for 13 of 37 total migration days, turtle G had the same dates of migration for 3 of 10 total migration days, and turtle L had the same dates of migration for 13 of 27 total migration days (see Figure 2 for migration dates). These four turtles tracked twice also had very similar migration paths across years (Figure 5).

Table 2 Migration dates for Dry Tortugas loggerheads, predicted from switching state-space modeling

Figure 2

Migration timing for post-nesting loggerheads. Migration days for female loggerhead sea turtles (Caretta caretta) that migrated to the Bahamas after nesting at Dry Tortugas National Park (DRTO), Florida. Darker shading indicates turtles that were tracked twice from DRTO in different years.

Figure 3

Proportion of post-nesting loggerheads in migration mode. Migration dates from 23 available tracks of female loggerhead sea turtles (Caretta caretta). Turtles were tagged while nesting at Dry Tortugas National Park and then migrated to residence areas in the Bahamas.

Figure 4

Migration paths for post-nesting loggerheads. Migration paths for 22 female loggerhead sea turtles (Caretta caretta) from nesting grounds in Dry Tortugas National Park (DRTO) to residence areas in the Bahamas. Dates of migration were determined through SSM; locations received on those dates shown here, filtered for swim speed, and then connected with a line. Turtle D had only one received location during SSM-determined migration dates and is not shown here.

Figure 5

Inter-annual repeatability of migration paths. State-space model (SSM) predicted tracks for four female loggerhead sea turtles (Caretta caretta) tagged while nesting in Dry Tortugas National Park (DRTO) during their migration to the Bahamas after two separate nesting seasons. The first predicted migration is shown as a blue line and the second, 2–3 years later, is shown as orange dots. For these turtles, the paths were very similar across years.

State-space modeling and residence areas

We obtained SSM results for 23 tracks, representing the paths of 19 individual turtles, four turtles tracked twice in successive nesting seasons (see Additional file 2 for an example SSM prediction figure and Additional file 3 for example model parameters). Of these, 21 tracks had periods classified as ARS mode after a migration (two turtles remained in migration mode after inter-nesting for the period in which they were tracked); however, we deemed these ‘residence areas’ because we have not yet directly observed foraging behavior for these individuals in the Bahamas. Eighteen tracks settled into only one residence area whereas three tracks settled into more than one residence area throughout the study period (i.e., periods at residence areas were interspersed with SSM-determined migration points). Nineteen tracks had enough mean daily locations during residence periods for KDE analysis; all 19 displayed site fidelity to these residence area locations (proportion of tracks that were more constrained than random movement paths >99.0099, Table 1).

The overall size of 50% core-use residence areas ranged from 3.7 to 179.5 km² (mean ± 1 SD = 56.2 ± 49.5 km²; see Additional file 4 for all KDEs and Table 1 for 95% KDE sizes). In sum, we obtained 2,160 total mean daily locations for analyses; however, because not every tracking day provided a turtle location, the time period during which turtles were at residence areas differed slightly from the number of mean daily locations. Overall, turtles occupied Bahamian residence area locations for a total period of 4,958 days (range 2–667; mean ± 1 SD = 206.6 ± 205.3 days); several tags were still transmitting daily location data at the time of data synthesis and write-up (1 December 2013) and no turtles showed movement away from the residence areas before transmission ceased or at the time of write-up (with the exception of turtle G that was tracked intermittently back to DRTO nesting grounds and recaptured in 2012; return track not shown). We also calculated three MCPs (see Additional file 4 for figure showing the MCP areas); the time periods at MCP residence areas ranged from 3 to 48 days (mean ± 1 SD = 22.7 ± 23.0 d) and the size of the areas ranged from 112.1 to 899.5 km² (mean ± 1 SD = 583.97 ± 416.3 km²; Table 3).

Table 3 Summary of minimum convex polygon (MCP) residence areas

Spatial configuration of residence areas

Distances to the nearest land from centroids of 50% KDEs at residence areas ranged from 1.2 to 72.9 km (mean ± SD = 46.2 ± 24.5 km). Bathymetry values (i.e., a proxy for water depths) at these centroid locations were relatively shallow and ranged from −15.0 to −2.0 m (with the exception of turtle R at −88.0 m; mean ± SD = −10.5 ± 19.2 for all locations: Table 1). Similarly, distances to the nearest land from centroids of the three MCPs at residence areas ranged from 29.4 to 46.1 km (mean ± SD = 39.4 ± 8.8). Bathymetry values at these locations ranged from −9.0 to −3.0 m (mean ± SD = −5.6 ± 3.1 m; Table 3). For all centroids combined (both KDE and MCP), the distance to the nearest centroid of another turtle ranged from 2.8 to 192.3 km (mean ± SD = 52.2 ± 50.3 km; Table 4), with minimal overlap of residence areas (see Additional file 5 for figure showing overlap of residence areas); this result suggests at least some separation in residence areas among individuals. In addition, turtle D had two periods of residency within the same year (2010) and we successfully created an MCP for one period and a KDE for the other; the distance between these two centroids was 6.1 km (Table 5).

Table 4 Distances to closest centroid of another turtle for Dry Tortugas loggerheads foraging in the Bahamas

Table 5 Centroid distances for turtles with more than one centroid location

Repeat tracks

Four turtles were recaptured 2–3 years after initial capture and tracking periods, and we re-tagged each one after they returned to nest in the Dry Tortugas. One of these turtles (turtle E/E*, tagged in 2010 and 2012) did not have enough locations in 2010 at Bahamian residence areas for KDE or MCP analysis but did have enough locations for KDE analysis in 2012 (see Table 1 and Figure 6). For the other three turtles, we successfully created KDE areas and calculated centroids for both tracking years (Figure 6). Successive residence areas selected were remarkably close together; the distance between KDE centroid locations from 2010 and 2012 for turtle G/G* was 1.5 km, the distance between KDE centroid locations for turtle L/L* from 2011 and 2013 was 2.9 km, and the distance between KDE centroid locations from 2009 and 2012 for turtle B/B* was 0.7 km. These turtles essentially selected the same residence area after two separate nesting seasons (Table 5; Figure 6). Further, for the three turtles tagged in different years with KDEs (B/B*, G/G* and L/L*), the average annual pairwise overlap of 50% kernel contours at residence areas was 0.28 (SD = 0.09).

Figure 6

Inter-annual repeatability of residence areas. Residence areas for four loggerhead sea turtles (Caretta caretta) recaptured 2–3 years after initial capture in Dry Tortugas National Park (DRTO). Turtle E’s points from 2010 did not produce a KDE or MCP, so the last filtered location received from that tracking period is shown as a small square. 50% Kernel Density Estimates (KDE) outlines are shown for the other three turtles for both tracking years with centroids (matching colored points). The distance between KDE centroid locations is 1.5 km for turtle G/G*, 2.9 km for L/L*, and 0.7 km for turtle B/B*.

Reserves and residence areas

We obtained coordinates for seven turtles for plotting from other investigators, and 19 to plot from visual estimation (Figure 7). These additional coordinates were located proximal to DRTO loggerhead residence area centroids. Overall, only one core-use residence area intersected a current Marine protected area (MPA): turtle O overlapped with a Wild Bird Reserve named Mammy Rhoda Cay. However, this reserve is small and only covers 0.01 km² of turtle O’s 107.3 km² 50% KDE area (Figure 8). Our 22 centroids (19 KDE and 3 MCP) ranged from 2.8 to 141.5 km from the nearest MPA (mean ± SD = 84.8 ± 36.3 km). These MPAs included Central Andros (National Park, 3 centroids), Channel Cays and Flat Cay (Wild Bird Reserve, 5 centroids), Little Inagua (National Park, 1 centroid), Mammy Rhoda Cay (Wild Bird Reserve, 1 centroid), Washerwomans Cut Cays (Wild Bird Reserve, 11 centroids), and Water Cay (Wild Bird Reserve, 1 centroid; Figure 8).

Figure 7

Loggerhead residence areas in the Bahamas. Locations of residence area centroids for loggerhead sea turtles (Caretta caretta) in the Bahamas including 19 tracks [50% KDEs, purple circles, this study], three MCPs (yellow triangles, this study), and locations from other studies (green squares). Residence area locations were provided by other authors (light green squares) or visually estimated from publications (dark green squares and underlined number). Letters indicate turtle designation associated with this study. Numbers indicate source and are as follows: 1) M. Arendt [pers. comm], 2) A. Foley [pers. comm], 3) D. Griffin [pers. comm], 4) [21], 5) [59], 6) [60].

Figure 8

Loggerhead residence areas and MPAs. Marine protected areas in the Bahamas marine ecoregion [61] with 50% KDEs (variously colored polygons) for 19 female loggerhead sea turtle (Caretta caretta) tracks and MCPs for three (yellow stars).

Discussion

Our results provide a previously unknown link for loggerheads between Dry Tortugas nesting beaches and Bahamas residence areas; out of 39 loggerhead nesting females satellite-tagged during 2008–2013 in DRTO, 17 (43.6%) migrated to the Bahamas and took up residence at distinct sites [see www.seaturtle.org/tracking/?project_id=402]. Using SSM, we determined that these DRTO loggerheads traveled along similar migration paths through the Florida Straits to reach the Bahamas. Additionally, turtles tagged in more than one nesting season showed remarkably similar migration paths and timing (both throughout and beyond the Florida Straits) as compared to their own previous migrations.

Residence areas

Once at their residence areas in the Bahamas, turtles selected sites that were relatively shallow (mean of −9.9 m for KDE centroids) and close to land (on average 46.7 km from land for KDE centroids). Hart et al. [3] described foraging areas for loggerheads in the GoM, four of which nested in DRTO. Those ten GoM-foraging loggerheads that nested on three different beaches throughout the GoM had a larger mean core-use area (50% KDE; 91.8 km²) than those in the Bahamas (56.2 km²). Similarly, GoM-foraging loggerheads had larger mean home ranges (95% KDE; 498.27 km²) than loggerheads taking up residence in the Bahamas (321.2 km²). Despite the variation in area, Bahamas residence area characteristics are generally similar to foraging sites defined by Hart et al. [3] in the GoM; they are also located in relatively shallow waters close to shore (GOM: water less than −50 m and within a mean distance of 58.5 m of the coastline).

Previous researchers that have tracked nesting loggerheads to the Bahamas did not all assess residence area core-use or home-range size [10,21,60]. However, residence area estimates from one adult female tracked to the Bahamas (56.6 km², 50% KDE; [62]) and one adult male (174.7 km², MCP; [10]), were both within the range of areas calculated for the 50% KDEs for our Bahamas loggerheads (3.7–179.5 km²). Hawkes et al. [22] provided home-range estimates for 13 foraging loggerheads (‘year-round’ turtles), two of which traveled to the Bahamas. The median home-range estimates for these 13 turtles that selected areas in South Carolina, Georgia, Florida, and the Bahamas ranged from 1,889.9 km² (α-hull) to 4,371.9 km² (MCP) depending on estimator method; Bahamas-only estimates were not published. Compared to these previously published estimates, our loggerhead home-range estimates (95% KDE) for residence areas in the Bahamas were smaller and ranged from 16.5 to 1,053.3 km² (mean 321.2 km²).

Possible residence area territories

Repeat tracking of several individuals showed remarkable re-selection of a particular residence area patch after completion of a migratory/nest/migratory cycle between DRTO and the Bahamas. The centroid distances between a turtle’s first residence area and the same turtle’s residence area selected after a nesting trip to DRTO ranged from only 0.7 to 2.9 km apart, and their KDEs overlapped substantially. This result demonstrates that loggerheads utilized the same discrete residence area upon returning to the Bahamas after nesting.

In addition to residence area repeatability, each turtle seemed to select a discrete patch of residence habitat; most turtles did not have substantial overlap of their residence areas with a different turtle. The one exception occurred with turtles K and Q; however, while turtle K’s KDE overlapped substantially with turtle Q, turtle Q’s KDE was much larger and these turtles only shared approximately 22 km² of turtle Q’s 160 km² residence area (14%); their residence area centroids were also 2.8 km apart, suggesting that these two turtles were not consistently sharing the same location. Centroid distances across all turtles (different individuals) ranged from 2.8 to 192.3 km apart. Site fidelity was observed at all residence areas, which covered almost 5,000 turtle days. Satellite tags are criticized for their higher locational error as compared to GPS tags [63]. However, we were able to discern these patterns in residence area separation and repeatability and to remotely observe (i.e., track) some turtles for long periods at distinct residence sites (e.g., turtle A, 744/817 tracking days (91.1%) at residence area). Further refinement of residence area separation could entail use of GPS tags deployed on individuals at putative foraging grounds.

Based on delimited residence areas that do not substantially overlap, we suggest that loggerheads at this putative foraging ground may establish territories in which they remain for long periods. Repeated tracking of individuals revealing consistent migration corridors to these sites and selection of the same discrete residence areas in successive tracking periods also supports this hypothesis. Recent studies have found long-term site fidelity of loggerheads at foraging grounds (adult and sub-adult loggerheads captured <5 km from original capture sites over time periods >1 year in Australia [64]) and have also observed little overlap in residence areas (49.0 km mean distance between individual core-use areas in the GoM [3]).

A territory, defined as an area within an animal’s home range over which it may have exclusive or priority use [65], may offer optimal quality food resources, rest sites, and access to mates. Fretwell and Lucas’ (1970) theory of ‘ideal free distribution’ [66] predicts that foraging animals will act to maximize foraging efficiency, due to a level of knowledge about forage resource profitability [67], and therefore possibly establish territories. While loggerheads may establish foraging territories and return to these sites after nesting, it remains unknown whether they ‘defend’ these discrete territories against conspecifics. Female-female aggression including threat displays, biting, and chasing, has been documented in Mediterranean loggerheads during breeding periods [68], but observations of loggerheads interacting at Bahamian residence areas are lacking. Determining available food resources and rest sites through benthic characterization as well as ascertaining fine-scale loggerhead behavior at these residence areas would therefore be useful for understanding how distribution of turtles may match that of ‘optimal’ resources and for verifying not only foraging behavior but also existence of potential loggerhead territories.

Conservation implications

There is little human development near these residence areas to threaten loggerheads; however, an Environmental Impact Assessment was completed in 2012 for the Bahamas Petroleum Company (BPC) to begin exploratory drilling on the continental shelf southwest of Andros, adjacent to the Cuba/Bahamas marine border and within BPC license blocks covering about 3 million acres [69]. These license blocks overlap with many of the loggerhead residence areas identified in this study. Although plans are focused on waters >300 m [69], if drilling begins, any spill could have far-reaching effects for this biodiversity hotspot.

Only one residence area intersected a MPA—a Wild Bird Reserve—and that overlap was minimal (<1% of the turtle’s core-use area). Thus, current MPAs do not offer protection to resident loggerheads tagged in DRTO. Additionally, during migration, the loggerheads passed through the Florida Straits, a major shipping fairway, and possible ship strikes could threaten the turtles on their migratory pathway. These turtles are therefore migrating through non-protected, potentially dangerous waters of the Florida Straits and later primarily settling at residence areas in non-protected waters once they reach the Bahamas. To inform better management of the DRTO turtles within international exclusive economic zones (EEZs) and on high-seas habitats, dynamic ocean management, in which policies can shift in space and time depending on resource conditions and human use [6], may be necessary.

According to the Food and Agriculture Organization (FAO) of the United Nations, commercial fishing takes place on the continental shelves of the Bahamas (both the Great Bahama Bank and the Little Bahama Bank). The most recent FAO summary [70] estimates that 95% of fishing efforts target spiny lobsters and primarily use casitas (sheets of zinc-coated metal placed over wood or concrete blocks), but also lobster hooks, compressors, spears, lobster traps, and fish traps. The live-weight estimate in 2007 for recorded landings of spiny lobster totaled 6,976 tons. In addition to lobster, the Bahamas fisheries also targets conch, crabs, and various species of fish. As loggerheads frequently forage on benthic invertebrates [71], the Bahamas fishery has the potential to impact loggerhead food resources, although further study is warranted to assess the extent of this impact. The fishery also poses a direct threat to loggerheads as they (and other turtles) can become entangled in lines attached to gear. Currently, direct turtle harvest is illegal in the Bahamas and has been since 2009 (K. Bjorndal, pers. comm.). The Convention on International Trade in Endangered Species of Wild Fauna and Flora [72], a global treaty regulating the trade in wildlife and wildlife products, entered into force in 1975 and prohibited commercial international trade in all Caribbean marine turtles and their products by 1977. However, Buchan [73] clearly indicated that turtle harvest continued into at least 1998 in his compilation of total landings of marine products and the FAO listed 2007 loggerhead catches in the Bahamas equivalent to about 1-ton live-weight [70]. Turtles are still harvested in the Turks and Caicos Islands, adjacent to the southeastern Bahamas [74].

Within the Bahamas EEZ, the government only allows Bahamian-owned fishing vessels; however, the FAO summary [70] claims the Bahamian government faces challenges in terms of enforcement and assessing stocks over the vast area of the continental shelf. These challenges cause some conservation concern as loggerheads have also been harvested from the Bahamas by Cuban fishermen [60].

Conclusions

Female loggerheads traveling from DRTO to the Bahamas show a high level of repeatability in both migration corridors and residence area use, and a low amount of residence-area overlap with conspecifics, suggesting they may establish territories. These residence areas primarily occur outside protected areas in the Bahamas. Therefore, international collaborative conservation efforts between the USA and Bahamas that focused on protection of loggerhead residence areas could offer a significant level of protection for DRTO loggerheads, a genetically distinct subpopulation of a species threatened with extinction. The Bahamas is also an important residence area and putative foraging site for two other subpopulations of NW Atlantic loggerheads (i.e., those from Northern and Peninsular FL subpopulations) that would benefit from these conservation efforts. Whereas the U.S. is a signatory on the Inter-American Convention (IAC) for the Protection and Conservation of Sea Turtles, the Bahamas has not yet committed to being a contracting party. Thus, designation of protected areas for critical turtle habitat and restriction of human activities that could harm turtles has not yet occurred. Our residence area location data would be valuable for consideration of critical turtle habitat should measures be implemented to reduce threats to loggerheads in the Bahamas. In addition, future studies to characterize the environment (i.e., oceanography) and environmental influences on the movements and behaviors of loggerheads will help to inform potential real-time efforts to further management and conservation efforts. Direct observations of foraging behavior for loggerheads resident at these sites are also needed to refine our understanding of foraging resources and extent. As well, information on movement and habitat use of male loggerheads in the Bahamas is lacking and in-water captures of males in the Bahamas is warranted to determine any sex-based differences in residence area selection.

Methods

Ethics statement

Experimental research on turtles followed internationally recognized guidelines. All turtle handling and sampling were performed according to NMFS-SESC 2008 and the USGS Institutional Animal Care Protocol USGS-SESC-IACUC # 2011–05.

Turtle capture and transmitter deployment

We tagged and outfitted loggerhead females with satellite transmitters after they nested in Dry Tortugas National Park on either East Key or Loggerhead Key between 2008 and 2013 (Table 1, Figure 1); as in Hart et al. [18], all tagging followed established protocols [75]. We tagged each animal with a passive integrated transponder (PIT) tag in the right shoulder region and affixed an individually numbered flipper tag to each of the front flippers. Immediately after tagging each animal, we took standard carapace measurements, including curved and straight carapace lengths (CCL and SCL, respectively).

We fitted a Wildlife Computers SPOT5 platform terminal transmitter (PTT) to each turtle; each tag (2.5xAA model, 71 mm long × 54 mm wide × 24 mm high) had a saltwater switch and output of 0.5 W and a mass of 115 g in air. Prior to transmitter application with either Power-Fast™ or SuperBond™ 2-part cool-setting marine epoxies, we removed epibionts (e.g., barnacles, algae) from the carapace of each turtle and sanded and cleaned the carapace with isopropyl alcohol. We streamlined attachment materials and minimized the epoxy footprint. The anticipated battery life of each tag was 1 year, and each tag in 2008, 2009, and 2010 was programmed to operate continuously; in 2011, 2012, and 2013, we adjusted winter tracking (November through April) to every 3rd day to prolong battery life. All tagged turtles were released within 2 h at their capture location.

Sea turtle tracking

Location data were filtered using Satellite Tracking and Analysis Tool [76] available on www.seaturtle.org/stat. Location classes (LC) 3, 2, 1, 0, A, and B were used to reconstruct routes and calculate straight-line and total distances that the turtles traveled. Locations were rejected if they were LC Z (for which no error estimation was available). Argos assigns accuracy estimates of <250 m for LC 3, 250 to <500 m for LC 2, 500 to <1,500 m for LC 1, and >1,500 m for LC 0 [77]. The estimated accuracy is unknown for LCs A and B, and locations failing the Argos plausibility tests are tagged as class LC Z. Argos performed both traditional least-squares location processing (1998–2010) as well as Kalman filtering (initiated in 2011, [78]) on location data. This newly implemented Kalman-filtering algorithm provides more estimated positions and significantly improves position accuracy, most significantly for locations obtained in LCs A and B [79].

Switching state-space modeling

We used switching state-space modeling [40,44] to estimate location and differentiate between migration and area-restricted search (ARS; to represent residence areas and inter-nesting) movements of loggerhead turtles by modeling and accounting for uncertainty due to data gaps and positional errors. The model was described in 2005 [41] and has previously been applied to the movement of marine animals including turtles [1,3,4,42,43,46-50,56,80]. Location data obtained through satellite transmitters are often received at irregular time intervals and sometimes involve large gaps and positional errors. Ad hoc filtering of location data based on location quality is not sufficient to remove erroneous locations and also results in loss of information [42]. Switching SSM estimates location and behavioral mode at regular time intervals, accounting for satellite positional errors and dynamics of the animal movement pattern [41] and is recommended as the best analytical technique for Argos tracking data once post processed by removing land points and adding back in good Argos locations [56].

We used the switching SSM approach to determine time periods in migration and the arrival date of each turtle at its respective residence area following Hart et al. [3] and Shaver et al. [4]. The behavioral mode was binary, defined as ‘ARS’ or ‘migration’ as in earlier applications [42,44,45]; however, since we tagged animals during nesting seasons, we defined the behavioral mode as ‘residence areas and/or nesting’ or ‘migration’. We summarized the Argos tracking data for the periods after migration away from nesting beaches, and during time periods with ‘residence area’ locations. Extended temporal gaps between locations can cause erroneous model predictions; therefore, we used satellite locations before gaps if there were long data gaps. We applied a model used by Breed et al. [45], which is a modified version of a model described by Jonsen et al. [41] that estimates model parameters by Markov Chain Monte Carlo (MCMC) using WinBUGS via the software program R. We used all tracking data except for LC Z. We fit the model to tracks of each individual turtle to estimate location and behavioral mode every 8 h from two independent and parallel chains of MCMC. Our samples from the posterior distribution were based on 10,000 iterations after a burn-in of 7,000 and thinned by five.

Delineation of migration periods

After fitting the switching SSM to individual loggerhead tracks, we identified locations where turtles were in migration or ARS mode. For migration periods, we summarized migration both temporally (number of days) and spatially.

Delineation of in-water residence areas

During ARS periods of residence in the Bahamas identified by SSM, we considered raw data locations to be biologically plausible if they were retained after passing through a speed and topography filter (points on land or with an implausible pathway across land, very distant, or with speeds >5 kph were removed). If an individual residence area period was at least 20 days in length, we generated mean daily locations to minimize autocorrelation using the filtered locations; the resulting coordinates provided raw data for KDE analysis. We used the Home-Range Tools for ArcGIS extension [81] and fixed-kernel least-squares cross-validation smoothing factor (h _cv ) for each KDE [59,82]. We calculated KDEs for turtles that had at least n = 20 mean daily locations at a residence area. When we observed unequal variance of the x and y coordinates, we followed procedures in previous studies and rescaled the data to select the best bandwidth [59,61]. We used ArcGIS 9.3 [83] to calculate the in-water area (km²; UTM or North America Albers Equal Area Conic projections were used to minimize error in area measurements) within each kernel density contour (50% and 95%) and to plot the areas; we used 95% KDEs to represent the overall home residence area and the 50% KDEs to represent core zone of activity at residence areas [84]. We summarized data for residence periods until the transmitters stopped sending information or at the time of data synthesis (n = 23, 23 October, 2013).

Some turtles continued transmitting after the completion of the SSM analysis. Once turtles were in ARS mode, and if they continued to transmit from a residence area, we assumed they remained in ARS mode. We combined these additional points with SSM-identified ARS locations and used the additional data in KDE and MCP estimates of residence areas and in site-fidelity tests (see below); this added 2,353 location days to the analysis.

Because the data set includes several turtles that were re-tracked in successive nesting seasons, following [85], we quantified the inter-annual variation in core residence areas as (2a_overlap/(a1 + a2), where a1 and a2 are the areas of the two respective 50% kernel contours and a_overlap is the area of overlap.

We tested for and quantified site fidelity within residence areas using the Animal Movement Analysis Extension for ArcView 3.2. Using Monte Carlo Random Walk simulations (100 replicates), we tested tracks during a turtle’s time at the residence area for spatial randomness against randomly generated walks [84]. We bounded the range for random walks from 2,500 m to 0 m bathymetry to include only the realistic extent of the in-water habitat for our animals during the study period. The Bahamas consists of deep oceanic canyons interspersed by atolls with more shallow neritic areas and so this boundary included depths deeper than the generally accepted ‘foraging’ depths of loggerheads (within the neritic zone, ≤ 200 m, see [11]). Tracks exhibiting site fidelity show movements that are more spatially constrained rather than randomly dispersed [84].

To further characterize at-sea residence areas selected by individual loggerheads, we calculated the centroid of each turtle’s 50% KDE; if a 50% KDE included multiple activity centers, we calculated the centroid of the largest activity center. We also calculated the depth at each centroid and the distance from each centroid to the nearest land; we used the North America Equidistant Conic projection for distance measurements in order to minimize error. For bathymetry, we used the GEBCO_08 Grid [86], a 30 arc-second continuous terrain model of both ocean and land.

For residence areas (as defined by SSM) without 20 mean daily locations, and more than three satellite locations, we performed MCP analysis (100% of points [51,53]) using ArcMap 9.3 [83]. Points for MCP analysis were filtered in the same way as those for KDE analysis (by speed and topography). We then calculated the centroid of these MCPs, as well as the depth and distance from the land for each centroid. If an individual turtle had multiple centroids (MCP and/or KDE), we summarized the spatial separation between them.

To further depict the distribution of loggerhead residence locations in the Bahamas, we added residence locations from other studies onto our map. This included one location from M. Arendt [pers. comm], four from A. Foley [pers. comm.], and two from D. Griffin and L. Hawkes [pers. comm]. We also visually estimated five locations from Girard et al. [21], two from Dodd and Byles [60], and 12 from Tucker et al. [87]. We plotted MPAs in the Bahamas to determine the level of habitat protection available to DRTO loggerheads putatively foraging in and taking up residence in this area. We downloaded the 33 available MPA layers from Protected Planet [88] for the Bahamas marine ecoregion.