Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Leatherback Turtle Movements, Dive Behavior, and Habitat Characteristics in Ecoregions of the Northwest Atlantic Ocean

  • Kara L. Dodge ,

    kara.dodge@post.harvard.edu

    Affiliations Department of Biological Sciences, University of New Hampshire, Durham, New Hampshire, United States of America, Large Pelagics Research Center, University of Massachusetts Amherst, Gloucester, Massachusetts, United States of America

  • Benjamin Galuardi,

    Affiliation Large Pelagics Research Center, University of Massachusetts Amherst, Gloucester, Massachusetts, United States of America

  • Timothy J. Miller,

    Affiliation National Marine Fisheries Service, Northeast Fisheries Science Center, Woods Hole, Massachusetts, United States of America

  • Molly E. Lutcavage

    Affiliation Large Pelagics Research Center, University of Massachusetts Amherst, Gloucester, Massachusetts, United States of America

Abstract

Leatherback sea turtles, Dermochelys coriacea, are highly migratory predators that feed exclusively on gelatinous zooplankton, thus playing a unique role in coastal and pelagic food webs. From 2007 to 2010, we used satellite telemetry to monitor the movements and dive behavior of nine adult and eleven subadult leatherbacks captured on the Northeast USA shelf and tracked throughout the Northwest Atlantic. Leatherback movements and environmental associations varied by oceanographic region, with slow, sinuous, area-restricted search behavior and shorter, shallower dives occurring in cool (median sea surface temperature: 18.4°C), productive (median chlorophyll a: 0.80 mg m−3), shallow (median bathymetry: 57 m) shelf habitat with strong sea surface temperature gradients (median SST gradient: 0.23°C km−1) at temperate latitudes. Leatherbacks were highly aggregated in temperate shelf and slope waters during summer, early fall, and late spring and more widely dispersed in subtropical and tropical oceanic and neritic habitat during late fall, winter and early spring. We investigated the relationship of ecoregion, satellite-derived surface chlorophyll, satellite-derived sea surface temperature, SST gradient, chlorophyll gradient and bathymetry with leatherback search behavior using generalized linear mixed-effects models. The most well supported model showed that differences in leatherback search behavior were best explained by ecoregion and regional differences in bathymetry and SST. Within the Northwest Atlantic Shelves region, leatherbacks increased path sinuosity (i.e., looping movements) with increasing SST, but this relationship reversed within the Gulf Stream region. Leatherbacks increased path sinuosity with decreasing water depth in temperate and tropical shelf habitats. This relationship is consistent with increasing epipelagic gelatinous zooplankton biomass with decreasing water depth, and bathymetry may be a key feature in identifying leatherback foraging habitat in neritic regions. High-use habitat for leatherbacks in our study occurred in coastal waters of the North American eastern seaboard and eastern Caribbean, putting turtles at heightened risk from land- and ocean-based human activity.

Introduction

Highly migratory marine predators such as leatherback sea turtles encounter a diversity of habitats during their long-distance movements. Oceanographic processes create regional ecosystems with distinct rates of primary productivity and community structure [1]. Predators may exhibit different behaviors in response to region-specific environmental conditions, with some regions optimal for foraging and (or) breeding while others serve as migratory habitat between breeding and feeding grounds. Obtaining direct measurements of foraging behavior in migratory marine species is challenging since the animals are difficult to observe for extended periods of time. Studies often rely on measures of search behavior to distinguish foraging from transiting, with the underlying assumption that a foraging animal should increase time and search effort in resource-rich areas (i.e., area-restricted search behavior) and decrease search effort in areas with fewer resources [2]. Marine animal tracking data has been used to measure area-restricted search behavior through analyses of speed, turning angle, path straightness and first passage time [3][5], while switching state-space models have been used to statistically estimate animal behavioral modes (e.g., foraging vs. transiting) [6][8].

Leatherback sea turtles (Dermochelys coriacea) are far-ranging marine predators, capable of swimming thousands of kilometers between boreal and tropical latitudes [9][11]. In recent decades, satellite telemetry has demonstrated that leatherbacks can undertake annual migrations (defined here as the seasonal movement between regions/habitats based on favorable versus unfavorable conditions, after Dingle & Drake [12]) in the Atlantic, Pacific and Indian Oceans [13][19]. These extensive migrations take leatherbacks through a heterogeneous seascape where they experience strong differences in biological and physical oceanographic conditions. In oceanic habitat, inferred foraging behavior and foraging success of some leatherbacks have been linked to fronts, upwelling and downwelling zones, and mesoscale features [15], [19][26], while some individuals move continuously without associating with particular oceanographic features [27]. Although leatherbacks are most often associated with an oceanic lifestyle, some individuals make seasonal use of highly productive continental shelf and slope habitats [16], [19], [26], [28][31], residing in near-shore areas for several months [18], [32][33]. Continental shelf and slope waters are productive regions where spring bloom conditions can lead to increased seasonal abundance of plankton [29], [34][36]. Increased nutrient input from land, tidal-mixing, and wind-driven upwelling can trigger increases in scyphozoan populations, while physical discontinuities in shelf waters and along ocean fronts promote aggregation and retention of gelatinous organisms [37][41]. Productive water masses and fronts in oceanic and neritic regions are also targeted by commercial fishing operations with fixed and mobile gear, often resulting in incidental captures of sea turtles and other non-target species [42][46].

Broad-scale tracking studies over the past decade have given new insight into the relationship between leatherback behavior and their environment [8], [17], [19], [21], [27], [47][52], but most of these studies focused on the post-nesting migrations of adult females. Comparatively few studies have been initiated in leatherback foraging grounds where different sexes and age classes mix [18][19], [24], [47]. We need additional data on males and females at multiple life history stages to fully understand the habitat-use and environmental associations of leatherback in-water populations.

In the present study, we deployed satellite tags on adult male, female and sub-adult leatherbacks turtles captured off Massachusetts, USA. This is the first in-water tagging study conducted in leatherback foraging grounds in the US Atlantic, and one of three direct-capture studies of leatherback turtles worldwide. We collected geolocation and dive data to: 1) determine leatherback occupancy of distinct oceanographic regions in the North Atlantic; 2) characterize leatherback regional movements, dive behavior, and environmental associations; 3) identify seasonal high-use habitat and 4) determine key environmental features associated with leatherback search behavior in the NW Atlantic.

Methods

Ethics statement

This work was conducted under the authority of the National Marine Fisheries Service Endangered Species Act Section 10 Permit # 1557-03 and the University of New Hampshire IACUC # 060501 and #090402. Turtle disentanglement was conducted under the authority of NOAA 50 CFR Part 222.310.

Satellite telemetry

Leatherbacks were located off the coast of Massachusetts, USA (∼41°N, 70°W) from August 2007 to September 2009, and captured with a breakaway hoopnet (n = 11) [53] or accessed through the Massachusetts sea turtle disentanglement network (n = 9) (Fig. 1; Table 1). The breakaway hoopnet was deployed from a modified bowsprit on commercial fishing or research vessels, and the net was pursed over individual turtles at the surface. We moved netted turtles alongside the vessel to the stern and brought them on board using a custom-built stern ramp. Veterinary personnel conducted health evaluations on all turtles prior to tag attachment [54].

thumbnail
Figure 1. Leatherback turtle tagging locations off Cape Cod, Massachusetts, USA, 2007–2009.

Letters represent satellite tag deployments as listed in Table 1.

https://doi.org/10.1371/journal.pone.0091726.g001

thumbnail
Table 1. Summary data for twenty leatherback sea turtles equipped with satellite tags off Massachusetts, USA, 2007–2009.

https://doi.org/10.1371/journal.pone.0091726.t001

Twenty leatherbacks were fitted with Wildlife Computers, Inc. (Redmond, WA, USA) model MK10-A (n = 8) and MK10-AF (n = 12) ARGOS-linked satellite time depth recorders. The tags deployed in 2007 had flexible, plastic-coated metal baseplates that conformed to the leatherbacks' medial ridge. We worked with Wildlife Computers to improve this design in future seasons (2008–2009), resulting in the “ridge-mount” tag model specifically developed for leatherback turtles. We attached satellite tags directly to the leatherbacks' medial ridge following methods first developed by Lutcavage et al. [55]. Prior to attachment we cleaned and anesthetized the attachment site with Betadine (three applications), isopropyl alcohol and a topical freezing agent (ethyl chloride) [55]. We made two horizontal drill holes (4.5 mm diameter) in the ridge bone using an electric drill and orthopedic drill bit. Drill hole sites were marked with furazone ointment and spaced to match the spacing of the MK10-A and MK10-AF ridge-mount tag bases. Plastic-coated braided stainless steel tethers were passed through tygon sheathing in the two drill tracts. We created a base for the tag that conformed to the turtle's ridge using a two-part cold-curing silicone putty (Equinox Silicone Putty, Smooth-On, Inc., Easton, PA, USA), and the tag was positioned on top of the putty base. The two tethers were threaded through holes in the tag base and secured using corrodible stainless steel crimps to ensure eventual release of the tag.

Leatherbacks were measured to the nearest 0.1 cm (curved carapace length: CCL and curved carapace width: CCW) with a flexible fiberglass measuring tape, and ranged from 123.0 to 161.5 cm CCL (Table 1). We used CCL to classify turtles as adults (CCL≥145 cm) or sub-adults (CCL<145 cm), and we determined gender based on tail length of adult turtles [56]. Five sub-adult turtles were sexed based on presence of a penis, subsequent necropsy or evidence of nesting. We collected samples of blood and skin [54], [57], and all turtles were photographed, scanned for passive integrated transponder (PIT) tags, and given PIT and flipper tags if none were present.

The satellite tags transmitted Fastloc GPS locations (MK10-AF model only), ARGOS-derived locations (all tags) and dive information (depth resolution ±0.5m and temperature resolution ±0.05°C) via Service ARGOS (Toulouse, France) (Table 1). Ninety-five percent of Fastloc GPS locations are accurate to ±55 m [58], while ARGOS-derived location error varies by location class (LC) as follows: LC3<150 m, LC2 150–350 m, LC1 350–1000 m, and LC0>1000 m. ARGOS does not provide accuracy estimates for LCA and LCB locations, and LCZ are considered invalid.

We defined a dive as vertical movement below two meters for at least one minute. The number of dives within specified depth and duration ranges and the time spent within depth and temperature ranges were collected as frequency histograms based on preprogrammed bins (Table 2). Histograms were aggregated over four 6-hour periods in GMT: 0:00–5:59, 6:00–11:59, 12:00–17:59, 18:00–23:59. Tags deployed in 2007 and 2008 (n = 17) were programmed to transmit daily while tags deployed in 2009 (n = 3) were programmed to transmit daily from July to December, every other day from January to April, and every third day from May to June.

thumbnail
Table 2. Bin ranges of dive parameters from satellite tags deployed on leatherback sea turtles, 2007–2009.

https://doi.org/10.1371/journal.pone.0091726.t002

Environmental data

We selected environmental data likely to influence production and distribution of gelatinous prey [39], [59] and thus leatherback movements. Sea surface temperature (SST), surface chlorophyll a concentration (chl a), SST gradient, chl a gradient and bathymetry were used as potential predictors. We used SST and chl a gradients as a proxy for the presence of fronts [60]. SST data were obtained as a blended product available at the GHRSST website (http://www.ghrsst.org) as 9 km, daily averages. The blended SST product was derived from microwave SST data from three sources (Advanced Microwave Scanning Radiometer, Tropical Rainfall Measuring Mission Microwave Imager and WindSAT Polarimetric Radiometer) and InfraRed SST data from Aqua Moderate Resolution Imaging Spectroradiometer (MODIS). Chl a data were obtained from MODIS data as 2.5 km, 8 day averages at the NOAA ERDDAP website (http://www.coastwatch.pfeg.noaa.gov/erddap/info/erdMWchla8day/index.html). Bathymetry was determined using 1-minute gridded global relief data (ETOPO1) from the National Geophysical Data Center (http://www.ngdc.noaa.gov/mgg/global/). SST and chl a gradients were generated using the Belkin-O'Reilly oceanic front detection algorithm [60].

Data analysis

We filtered 30,173 raw ARGOS and GPS locations using Kalman filter methods outlined in Royer & Lutcavage [61]. Since our analysis included GPS data, we extended the original error covariance structure to include this information. Data were interpolated to a three-hour time step and smoothed. Environmental data were then extracted for the 28,253 filtered turtle locations. Rate of travel was determined using the distance function in Matlab (Mathworks, Natick, MA) and a daily straightness index was calculated as the ratio of straight-line distance to total distance traveled by each turtle per day (in km d−1), resulting in a dimensionless index from near 0 (sinuous) to 1 (straight) [62][63]. Travel rate and straightness were not calculated from January to June for the three turtles tagged in 2009 since tags were duty-cycled and did not transmit daily during these months. Leatherback positions and dive data were assigned to distinct biogeographic provinces or “ecoregions” defined by Longhurst [1]. We calculated the duration of leatherback occupancy in each ecoregion, and assessed variability in leatherback search effort, dive behavior and environmental associations across ecoregions. Averages are shown as mean ± standard deviation (mean ± SD) for normally distributed data, and median, interquartile range (Q1–Q3) where data are not normally distributed.

To investigate variation in seasonal habitat use, we created density utilization maps of filtered leatherback positions for pooled data across all turtles by season. Seasons were defined as: July – September (summer), October – December (autumn), January – March (winter), and April – June (spring). Daily locations were summed into hexagonal area bins, with the area of each hexagon approximately 669 km2 (or 4 hexagons per degree). These bins are larger than the error associated with our filtered ARGOS and GPS location data, but small enough to identify regional high-use areas. Density utilization maps were produced using R [64] and Generic Mapping Tools [65] (see http://www.scott.sherrillmix.com/R-GMT_HexPlot.php).

We applied generalized linear mixed-effects models to investigate the influence of ecoregion, SST, chl a, SST gradient magnitude, chl a gradient magnitude and bathymetry (fixed effects) on leatherback path straightness while accounting for the correlation of repeated observations from individual turtles [66]. Changes in path straightness have been used to identify purported search behavior associated with foraging in leatherbacks [25], as well as other marine predators [67][69]. Density plots of the environmental data showed that logarithmic transformation was required for SST gradient, chl a, chl a gradient and bathymetry. We were primarily interested in the influence of SST, SST gradient, chl a and chl a gradient in regions where these surface features are most variable (least homogenous), and we were mainly interested in the effect of bathymetry in neritic habitats where leatherbacks can access the entire water column. Therefore, we estimated region-specific regression parameters for these variables: SST (Northwest Atlantic Shelves, Gulf Stream), SST gradient (Northwest Atlantic Shelves, Gulf Stream) and bathymetry (Northwest Atlantic Shelves, Guianas Coastal). We compared models where effects of SST, SST gradient, chl a, chl a gradient and bathymetry on path straightness are the same for the respective groups of regions with less parsimonious models that allow the effects to differ for each ecoregion (Table 3). As path straightness ranges from near 0 to 1, we assumed a Gaussian error structure for the logit-transformed path straightness with continuous first-order autoregressive model correlation structure for repeated observations as a function of time between observations [70].

thumbnail
Table 3. Definition of fitted mean logit path straightness models for observation j of turtle i.

https://doi.org/10.1371/journal.pone.0091726.t003

We fit the models by maximum marginal likelihood in R [64] using the lme4 package and compared relative performance using Akaike Information Criterion (AIC) [71] of the fitted models. As a measure of evidence for relative performance of each model that we fit, we used Akaike weights [72]. We chose to use AIC over a criteria adjusted for sample size (AICc) because AICc requires a known number of independent observations, and this is not straightforward for mixed-effects models [73]. For first-order autoregressive models as used here, the effective number of independent observations depends on the correlation of the observations within each individual. Since the autocorrelation we estimated for our models is generally low, we had a large effective number of independent observations where differences between AIC and AICc are negligible. For the model that provided the best fit, we used restricted maximum likelihood to obtain parameter estimates and predict changes in steepness with various covariates and factors [74].

Results

Satellite telemetry

We received data from all tagged leatherbacks: four tags transmitted for less time than expected and 16 tags met or exceeded predicted battery life. Tags reported between 16 and 414 days, with a median tracking duration of 184 (152 to 219; Q1–Q3) days (Fig. 2; Table 1). We pooled the percent frequency distributions of the four dive parameters: dive-depth (n = 19 tags), dive-duration (n = 17 tags), time-at-depth (n = 15 tags), and time-at-temperature (n = 15 tags), shown in Fig. 3. Three tags deployed in 2007 reported spurious dive-duration data (i.e., the total number of dives recorded in the >52 min bin exceeded the 6-hour time period) and one 2007 tag recorded insufficient dive-depth data, so these tags were excluded from the dive-duration and dive-depth analysis. We only included time-at-depth and time-at-temperature data for tags with the same lower bin ranges (2008 and 2009).

thumbnail
Figure 2. Reconstructed movements of 20 satellite-tagged leatherback turtles in the North Atlantic Ocean, 2007–2010.

Tracks show turtle movements from point of release (Cape Cod) to point of last Argos transmission (red triangles). Tags were deployed on adult males (F, I, L, N, P, T), adult females (G, H, M), and sub-adults (A, B, C, D, E, J, K, O, Q, R, S).

https://doi.org/10.1371/journal.pone.0091726.g002

thumbnail
Figure 3. Frequency distributions of dive parameters from satellite tags on leatherback turtles in the North Atlantic Ocean.

a) Dive-depth (n = 19) and b) dive-duration (n = 17) from turtles tagged 2007–2009. c) Time-at-depth (n = 15) and d) time-at-temperature (n = 15) from turtles tagged from 2008–2009.

https://doi.org/10.1371/journal.pone.0091726.g003

Of the 210,556 dives reported for the dive-depth parameter, over 28% were to depths less than 5 m and 90% were shallower than 75 m (Fig. 3a). Fifteen turtles dove deeper than 500 m during the study period, with males recording the deepest dives (>1200 m, n = 3 turtles). The pooled dive-duration data showed that close to 75% of all leatherback dives were less than 12 min, and more than 90% were shorter than 32 min (Fig. 3b). Sixteen turtles recorded extended dives lasting over 52 minutes, but these represent less than 3% of the total. The pooled frequency distribution of the time-at-depth shows that turtles spent over 25% of their time within 2 m of the surface, and over 50% of their time shallower than 10 m (Fig. 3c). Over 90% of their time was spent in the top 100 m of the water column. The pooled frequency distribution of the time-at-temperature shows that turtles spent 86% of their time between 16°C and 28°C (Fig. 3d).

Habitat use and environmental associations

Leatherbacks ranged widely between 39°W and 83°W, and between 9°N and 47°N (Fig. 2), over six oceanographically distinct ecoregions defined by Longhurst [1]: the Northwest Atlantic Shelves (n = 20), the Gulf Stream (n = 16), the North Atlantic Subtropical Gyral West (hereafter referred to as the Subtropical Atlantic, n = 15), the North Atlantic Tropical Gyral (hereafter referred to as the Tropical Atlantic, n = 15), the Caribbean (n = 6) and the Guianas Coastal (n = 7) (Fig. 4). All leatherbacks were tagged in the Northwest Atlantic Shelves, and 16 turtles (Turtles D, F–T) were tracked long enough to determine an average minimum residency of 79 days (±39 days) in this region post-tagging (Table 1). Fifteen turtles (Turtles D, G–T) left the Northwest Atlantic Shelves between late September and mid-November, with the majority leaving between mid-October and mid-November (n = 11, Turtles D, G–K, N–O, R–T). Most turtles spent less than a week in the Gulf Stream region (median 6 days), but two individuals (Turtles F & D; Table 1 & Fig. 2) made more extensive use of this region (95 and 59 days, respectively). Between October and February, leatherbacks transited rapidly through the Subtropical Atlantic (median 29 days), entering the Tropical Atlantic between November and early February. Leatherbacks either remained in the Tropical Atlantic for the remainder of the tracking period (n = 7, Turtles D, G, J–K, O–P, R) or continued on to breeding and (or) foraging areas in Caribbean and Guianas Coastal regions (n = 8, Turtles H–I, L–N, Q, S–T). Three sub-adults and one small adult male were tracked long enough to observe a complete (32 and 92 days; Turtles D & R, Table 1) or partial (83 and 56 days; Turtles P & S, Table 1) overwintering period in the Tropical Atlantic before they returned to the Subtropical Atlantic between late March and mid-May. One turtle (Turtle R) returned to the Northwest Atlantic Shelves in mid-May, remaining in the region for 96 days before the tag stopped transmitting in late August.

thumbnail
Figure 4. Seasonal habitat use of adult male, female and sub-adult leatherback turtles, 2007–2010.

Summer, July – September (n = 19 turtles), autumn, October – December (n = 17 turtles), winter, January – March (n = 16 turtles), and spring, April – June (n = 5 turtles). There are four hexagons per degree; each hexagon represents approximately 669 km2. Color scale shows the number of track days per hexagon. Ecoregions from Longhurst [1]: NWCS, Northwest Atlantic Shelves; GFST, Gulf Stream; NASW, North Atlantic Subtropical Gyral West; NATR, North Atlantic Tropical Gyral; CARB, Caribbean; GUIA, Guianas Coastal.

https://doi.org/10.1371/journal.pone.0091726.g004

Turtles modified their movements and dive behavior while occupying different ecoregions. Leatherbacks in the Northwest Atlantic Shelves had the lowest travel rates and path straightness of all regions (Fig. 5a,b), and they combined slow, sinuous swimming with short, shallow dives (Fig. 6a,b). They spent most of their time in the top 10 m of the water column at temperatures between 16°C and 20°C (Fig. 6c,d). Outside of the Northwest Atlantic Shelves, turtles increased their travel rate and path straightness (Fig. 5a,b). Leatherbacks continued making shallow dives (<5 m), but increasingly made deeper, longer dives as they traveled south (Fig. 6a,b). As turtles moved into subtropical and tropical ecoregions, they began spending more time at temperatures over 22°C and at depths over 50 m, experiencing the warmest temperatures and making the deepest, longest dives in the Tropical Atlantic and Caribbean (Fig. 6c,d). Diving patterns in the primarily shelf waters of the Guianas Coastal region were distinct from other tropical regions: leatherbacks spread their diving effort throughout the top 200 m and most dives were less than 24 min (Fig. 6a,b). Leatherback surface times (0–2 m) were similar across most regions, with the greatest surface times recorded in the Subtropical Atlantic (mean 30%), Northwest Atlantic Shelves (mean 27%), Caribbean (mean 26%) and Gulf Stream (mean 25%). Leatherbacks spent the least amount of time at the surface in the Guianas Coastal (mean 14%) and Tropical Atlantic (mean 21%).

thumbnail
Figure 5. Leatherback behavioral indices in Longhurst regions of the North Atlantic Ocean.

a) Rate of travel, b) straightness, c) SST, d) SST gradient magnitude, e) chl a, f) chl a gradient magnitude, and g) bathymetric depth. NWCS, Northwest Atlantic Shelves; GFST, Gulf Stream; NASW, North Atlantic Subtropical Gyral West; NATR, North Atlantic Tropical Gyral; CARB, Caribbean; GUIA, Guianas Coastal. Boxplots: centerline, median; edges of box, 1st and 3rd quartiles; whiskers, data points within the range Q1−1.5(Q3−Q1) to Q3+1.5(Q3−Q1).

https://doi.org/10.1371/journal.pone.0091726.g005

thumbnail
Figure 6. Cumulative frequency plots of leatherback dives in Longhurst regions of the North Atlantic Ocean.

Depth and duration bins (top panels) and leatherback hours in depth and temperature bins (bottom panels). NWCS, Northwest Atlantic Shelves; GFST, Gulf Stream; NASW, North Atlantic Subtropical Gyral West; NATR, North Atlantic Tropical Gyral; CARB, Caribbean; GUIA, Guianas Coastal. The second y-axis (red) corresponds to NWCS (red line) while all other regions reference the first y-axis (black). This highlights the increased dive activity of tagged turtles in the NWCS region.

https://doi.org/10.1371/journal.pone.0091726.g006

Turtles experienced highly variable environmental conditions across ecoregions, where bathymetry ranged from shallow bays and sounds on the continental shelf to deep oceanic waters, SST from 9.6°C to 28.9°C and chl a from near zero to 64.36 mg m−3 (Fig. 5). Leatherbacks occupied areas with SST and chl a gradients of varying magnitudes, with SST gradients of 0 to 1.85°C km−1 and chl a gradients of 1.052 to 1.144 (ratio) (Fig. 5). Turtles in the Northwest Atlantic Shelves used relatively shallow habitat, staying mostly within the 80 m isobath and associating with the highest chl a of all regions (Fig. 5e,g). Turtles experienced the coolest, most variable sea surface temperatures and strongest SST gradients in the Northwest Atlantic Shelves and the Gulf Stream (Fig. 5c,d). While the chl a gradients were similar across ecoregions, leatherbacks in the Guianas Coastal region used habitat with the most variable and strongest chl a gradients (Fig 5f).

Based on AIC values, the most well supported model showed that differences in leatherback search behavior (represented by logit-transformed path straightness) were best explained by ecoregion and effects of bathymetry and SST, with effects of SST depending on the ecoregion (Table 4). In the mainly neritic ecoregions (Northwest Atlantic Shelves and Guianas Coastal), there was a positive relationship between path straightness and bathymetry, with leatherback movements becoming more sinuous as water depth decreased (Table 5; Fig. 7a). In the regions where SST was most variable (Northwest Atlantic Shelves and Gulf Stream), the relationship between path straightness and SST differed. In the Gulf Stream, the relationship between path straightness and SST was positive, with leatherback path sinuosity increasing with decreasing SST. In the Northwest Atlantic Shelves, the relationship between path straightness and SST was negative, with leatherbacks increasing path sinuosity with increasing SST (Table 5; Fig. 7b). The relationship between path straightness and SST in the Gulf Stream region was slightly positive, reflecting the fact that this slope parameter was not significant in the model.

thumbnail
Figure 7. Predicted straightness (solid lines) and 95% confidence intervals (dashed lines) from the best performing model in Table 4.

Leatherback path straightness is shown in relation to: a) observed log-transformed bathymetry (m) in each of two distinct ecoregions of the northwest Atlantic and b) observed sea surface temperature (SST) in each of two distinct ecoregions of the northwest Atlantic; Northwest Atlantic Shelves (NWCS); Gulf Stream (GFST); Guianas Coastal (GUIA). Mean SST value in each region was used to create the bathymetry plot and mean bathymetry value in each region was used to create the SST plot.

https://doi.org/10.1371/journal.pone.0091726.g007

thumbnail
Table 4. AIC results for linear mixed-effects models where mean logit-transformed path straightness is a function of ecoregion, sea surface temperature (SST), logarithm of surface chlorophyll a concentration (chla), logarithm of SST gradient (SSTg), logarithm of surface chla gradient (chlag), and logarithm of bathymetry (bathy).

https://doi.org/10.1371/journal.pone.0091726.t004

thumbnail
Table 5. Fixed effects parameter estimates of final model.

https://doi.org/10.1371/journal.pone.0091726.t005

Seasonal density utilization maps showed leatherback movements were the least extensive during summer, with turtles tagged off Massachusetts showing a strong preference for the Northeast US continental shelf, concentrating movements off southern New England and Virginia/North Carolina (Fig. 4 ‘Summer’). Leatherbacks expanded their range in autumn, increasing their use of the Mid-Atlantic Bight and Gulf Stream before initiating a rapid, directed southward migration, following widely dispersed pathways through the Subtropical and Tropical Atlantic (Fig. 4 ‘Autumn’). Winter movements were restricted to tropical oceanic habitat and neritic waters of the Antilles, South America, and Central America, with the exception of two sub-adults (Turtles D & R) and one small adult male (Turtle F) that occupied subtropical oceanic waters for part of the season. Densely occupied winter habitat occurred off breeding beaches in the Windward Islands, particularly the north/northeast coast of Trinidad and the western half of the Tobago Basin (Fig. 4 ‘Winter’). Our limited tracking data for the spring showed that leatherbacks either remained in tropical breeding areas in the Windward Islands or began northward migrations, with one individual (Turtle R) occupying the North Carolina shelf (Fig. 4 ‘Spring’).

Overall, dispersal patterns differed between adult and sub-adult leatherbacks. Most adults followed widely spaced but highly oriented south/southeast headings during their southward migration until they reached latitudes between 10°N and 13°N (Fig. 2). Sub-adult leatherbacks had more variable headings and did not disperse as far south, with most turtles remaining north of 15°N (Fig. 2). Two sub-adults (Turtles K & S) were tracked near land (the Bahamas and the Lesser Antilles) during a portion of their migrations, but most sub-adults occupied offshore tropical and subtropical habitat during winter and early spring. The four largest adult males (Turtles I, L, N, T) and two adult females (Turtles H & M) traveled to areas off nesting beaches in the Lesser Antilles and Central America where they remained until tag transmissions ceased. Two smaller adult males and one adult female did not travel to known breeding areas; one male (Turtle F) remained in the Gulf Stream into early January and the other two turtles (Turtles G & P) overwintered in a tropical region near the convergence of the North Equatorial Current and North Equatorial Counter-Current.

Discussion

We deployed GPS-linked and conventional ARGOS STDRs to simultaneously collect data on movements and dive behavior of adult and sub-adult leatherbacks in the North Atlantic. This study is one of the first to obtain highly accurate GPS locations from leatherback turtles, allowing us to identify high use habitat, movement patterns and environmental associations with less observation error. We also used novel design (e.g. “ridge-mount tag”) and direct attachment techniques to deploy low profile, hydrodynamic tags. In an analysis of transmitter drag and tag attachment procedures, Jones et al. [75] found that Wildlife Computers ridge-mount tags directly attached to the leatherback medial ridge resulted in the least amount of drag (0.6–1.8% increase in drag coefficient) whereas telemetry studies employing backpack harnesses resulted in the greatest amount of drag (92–112% increase in drag co-efficient). Field studies comparing the behavior of harnessed leatherbacks and leatherbacks with directly attached tags found significant differences in leatherback swim speed [51], [76] and dive duration [76]. We believe our tag model and tag attachment methods minimized impacts to the turtles' natural behavior, resulting in a more accurate portrayal of leatherback behavior than studies using harness attachments.

Dispersal patterns and seasonal habitat use

Leatherbacks tagged off Massachusetts showed a strong affinity to the Northeast US continental shelf before dispersing widely throughout the northwest Atlantic. One individual tracked for >1 year exhibited site fidelity to the US shelf, returning in late spring and remaining through late summer. Surprisingly, only one Massachusetts-tagged leatherback moved onto the eastern Canada shelf, an important and well-documented leatherback foraging ground [16], [32], [77][78]. In contrast, leatherbacks tagged off eastern Canada spent extended periods in both Canadian and Northeast US shelf waters within a single foraging season [16]. Over half of the turtles tagged off Massachusetts were classified as sub-adults (55%) whereas sub-adults were the minority in the Canada tagging study (16%, n = 38). The sub-adults in our study may have been at the northern part of their foraging range. Migration distance depends partly on an animal's body size and their capacity to store energy (e.g., adipose tissue) [79]. Larger animals are able to swim longer distances at lower energetic cost [80], as has been predicted for cod [81] and bluefin tuna [82]. Smaller body size and lower lipid stores may limit the migratory range of some sub-adult leatherbacks. Larger leatherbacks also have an enhanced capacity for retaining metabolically produced heat and heat they absorb from their environment as a result of a small surface area to volume ratio and greater peripheral insulation [83][84]. Their superior thermoregulatory abilities would allow larger adult leatherbacks to more effectively exploit the cooler waters of the Canadian shelf than sub-adults. Alternatively, resources may have been sufficient on the Northeast US shelf during the years of our study, precluding a longer migration to more northerly foraging areas. Adult female leatherbacks tagged on nesting beaches in the northwest Atlantic showed different habitat utilization patterns than those tagged in temperate foraging grounds: high-use habitat tended to occur close to the nesting beaches where turtles were tagged or turtles were more widely distributed in oceanic and neritic regions throughout the North Atlantic [85]. Although there is inherent bias in determining habitat utilization from tracking data from one deployment location, bias is reduced with increasing deployment time, and alternative methods may be used to address this in the future [86][87].

There was a strong seasonal component to habitat selection, with most leatherbacks remaining in temperate latitudes in the summer and early autumn and moving into subtropical and tropical habitat in the late autumn, winter and spring. This latitudinal shift is consistent with previous studies of leatherbacks tracked from foraging grounds in the North Atlantic [24], [32] and similar to seasonal migration patterns of other sea turtle species [88][89] and large pelagic species such as ocean sunfish Mola mola [90], basking sharks Cetorhinus maximus [91], bluefin tuna Thunnus thynnus [92][93], and swordfish Xiphias gladius [94]. Latitudinal shifts in habitat-use likely reflect seasonal changes in temperature, productivity and prey availability at high latitudes, and, for some species, the necessity to spawn or nest in tropical regions. At temperate latitudes, gelatinous zooplankton populations are highest during summer months following the spring phytoplankton bloom, and large scyphozoan species may be more abundant in these cooler, temperate waters due to enhanced trophic transfer efficiencies [26]. These factors likely drive leatherbacks' seasonal exploitation of temperate coastal habitat, while their thermoregulatory and breeding constraints compel them to return to subtropical and tropical habitat from late autumn to spring [26].

During the over-wintering period, sub-adults, small adult males, and a single inter-nesting-year female primarily remained in oceanic habitat, while large adult males and two reproductive females moved into coastal breeding areas. Little is known about the demographics of male leatherbacks, but there may be a size constraint whereby smaller males are unable to compete for females, and are effectively displaced from breeding areas by larger, more dominant individuals [79]. Smaller males may direct energy toward growth rather than reproduction, and select over-wintering habitat to maximize limited foraging opportunities in the tropics. Adult females in an inter-nesting year are likely to avoid breeding areas where they would be subject to aggressive mating attempts by males [95]; by overwintering offshore, they can save energy and accumulate fat stores for return migration and future reproductive effort. Sub-adults also largely avoided breeding areas and most did not travel as far south as adults. The highly oriented paths taken by adult leatherbacks suggests movement toward a goal (e.g., specific breeding and/or foraging areas), while the more variable headings taken by sub-adults may indicate an opportunistic overwintering strategy, or lack of experience locating consistent resource patches in oceanic habitat.

Regional movements, dive behavior, and habitat characteristics

Temperate neritic habitat.

Our density utilization maps demonstrate that the Northeast US shelf, particularly southern New England, provides important seasonal habitat for leatherback turtles tagged of Massachusetts. The Northeast US shelf is one of the most well-studied and productive large marine ecosystems in the world [1], [96]. The region is characterized by a nutrient-limited, spring production peak; increasing irradiance and stratification in spring lead to a phytoplankton bloom, with peak surface chl a biomass in April [1]. The spring bloom in this productive coastal ecosystem is dominated by large cell phytoplankton (diatoms), likely leading to increased trophic transfer efficiencies (reviewed by Saba [26]). Food resources resulting from the predictable bloom feed juvenile gelatinous zooplankton that will be become fully mature by summer when leatherbacks are present [26], [97]. In addition to the high productivity of the system, the Northeast US shelf provides the hard bottom substrate needed by the benthic polyp stage of most schyphozoans [26], [97], the favored prey of leatherbacks in coastal waters of the Northeast US [57] and eastern Canada [78].

Previous studies have identified the Northwest Atlantic Shelves as important foraging habitat for leatherbacks [16], [26], [33], [57], [77], [98], and we found their behavior (slow, sinuous swimming) in this region to be consistent with area-restricted search. Leatherback locations in the Northwest Atlantic Shelves coincided with higher surface chl a and stronger SST gradients (with the exception of the Gulf Stream) relative to other regions, and most dives were in the euphotic zone within the average seasonal mixed layer depth on the Northeast US shelf (10–20 m) [1]. Highly productive water masses and frontal zones influence the spatial distribution and movements of some top predators, aggregating them in relatively small areas or “hotspots” on the shelf, shelf break, slope, offshore and at depths where prey is concentrated [36], [38], [99][100]. Other sea turtle species associate with enhanced frontal activity [101][102], though some cheloniids may face thermal constraints, limiting their access to cooler temperate frontal zones [102]. Weak winds and strong stratification are typical summer conditions on the Northeast US shelf that could lead to increased retention of gelatinous organisms [26]. Leatherback movements coincided spatially and temporally with the persistent Shelf-Slope Front and tidal-mixing fronts north of Nantucket Shoals, in the Gulf of Maine, and around Georges Bank [103]. The tidal-mixing fronts occur during peak leatherback presence in summer and early autumn, and may play an important role in consolidating seasonally abundant patches of the leatherback's gelatinous prey. Our model results showed that leatherbacks increased path sinuosity at shallower depths and warmer surface temperatures within the Northwest Atlantic Shelves. Decreasing water depth has been linked to increases in epipelagic gelatinous zooplankton biomass on a global scale, with greatest biomass found in shallow locations (<10 m average depth) [59]. Shallow shelf habitat such as shoals, banks, and ledges may be important for leatherback prey searching and (or) foraging in this region, and is consistent with our field observations of leatherbacks feeding in shoal habitat off Massachusetts [57] (Fig. 1).

The average mixed layer depth on the Northeast US shelf is also shallowest (10–20 m) during the summer and early autumn [1], potentially aggregating gelatinous prey at or above the pycnocline [104]-[105]. This would reduce ascent and descent times for foraging leatherbacks, and minimize time spent in cool waters below seasonal thermoclines. The percentage of time that leatherbacks spent at the surface in the Northwest Atlantic Shelves (27%) was much lower than that observed by James et al. [106] for turtles tagged off eastern Canada (mean 43% night and 50% day, n = 12 (0–2 m) and n = 3 (0–3 m)). The Canada-tagged turtles used habitat off the Northeast US as well as eastern Canada, so the disparity in percent surface times is somewhat surprising, but the colder water temperatures experienced by leatherbacks off eastern Canada may contribute to increased surface times in that region. James et al. [32], [106] described regular observations of leatherbacks basking at the surface off Nova Scotia, a behavior that we rarely observed off Massachusetts. It is possible that leatherbacks off eastern Canada spend greater time at the surface for thermoregulation [32] but this behavior is less important in the comparatively warm waters of the Northeast US shelf.

Subtropical oceanic habitat.

We observed marked behavioral changes as leatherbacks left continental shelf habitat and began their southward migrations through subtropical oceanic habitat. As turtles moved through the Gulf Stream and Subtropical Atlantic, they showed rapid, directed travel and began spending more time at depths >50 m. Most turtles spent minimal time in these ecoregions, suggesting that these are less important feeding areas for Massachusetts-tagged leatherbacks and are primarily used for transiting between temperate (i.e., foraging) and tropical (i.e., breeding) habitat. However, two individuals did make more extensive use of the Gulf Stream region during summer, fall and early winter. The Gulf Stream's strong horizontal SST gradient, particularly in fall and winter, is evident in the strong SST fronts encountered by leatherbacks there. Leatherback movements in the Gulf Stream became slightly more sinuous at lower SST, possibly associated with upwelling along the Gulf Stream front, but this relationship was weak in our model. The Gulf Stream has been previously identified as probable foraging habitat for leatherback turtles [20], [25], [49], and turtles exploiting this area may take advantage of the enhanced productivity of energetic, mesoscale eddies characteristic of this region [1].

Leatherback movement patterns and dive behavior in the Subtropical Atlantic were consistent with other studies of this species in the North Atlantic [22], [25], [32], [49], [106]. Observed changes in dive behavior in the Subtropical Atlantic may be partially explained by cycles of seasonal stratification and depth of the mixed layer and (or) deep chlorophyll maximum. In the Subtropical Atlantic, the mixed layer deepens in the fall and winter when leatherbacks are present, and average mixed layer depth is over 50 m [1]. If leatherbacks target the pycnocline to search for prey, they would increase their diving activity to depth strata >50 m, with longer durations associated with increasing ascent and descent times. Leatherbacks in Subtropical Atlantic may also be exploiting the intense trophic activity associated with the deep chlorophyll maximum, which occurs at about 100 m [1]. Consumption at the deep chlorophyll maximum follows a strong diel cycle [1], and leatherbacks may be capitalizing on enhanced nighttime prey availability by making deeper, longer dives [106]. However, the rapid transit rate and limited time spent in the Subtropical Atlantic suggests that this region is not as important to over-wintering leatherbacks as tropical regions. Reproductive adults would have additional incentive to move quickly through the Subtropical Atlantic to reach breeding and nesting areas in the tropics.

The percentage of time that leatherbacks spent at the surface in the Gulf Stream and the Subtropical Atlantic was similar to the percent surface times recorded by James et al. [106] in the morning (06:00–12:00 GMT; mean 29%) and evening (18:00–0:00 GMT; mean 29%) periods of the southern migration, but much lower than their percent surface time during the day (12:00–18:00 GMT; mean 77%). Considering only the day period in our data set (12:00–18:00 GMT), the average percent surface time is still much lower (44%). Turtles tagged off eastern Canada and Massachusetts had similar dispersal and migratory patterns in the subtropical gyre, so the observed difference in surface times may be due to other factors such as different demographics (e.g., sex and body size) of our turtle sample or tagging technique (harness vs. direct attachment). Comparisons of these techniques showed that leatherbacks had lower travel rates [51], [76] and shorter dive durations [76] when wearing a harness. While surface time was not directly addressed in these comparisons, it's possible that leatherbacks wearing harnesses increase their surface time to recover from the energetic costs of increased drag caused by the harness (91–112% increase in drag coefficient) [75].

Tropical oceanic and neritic habitats.

Leatherbacks overwintered in tropical ecoregions, with reproductively active adults primarily occupying the Guianas Coastal and Caribbean regions while non-reproductively active adults and sub-adults mainly used oceanic habitat in the Tropical Atlantic. Turtles slowed down in the tropics compared to the subtropical gyre but travel was still directed compared to the sinuous movements we observed in the Northwest Atlantic Shelves, suggesting a mix of behaviors that may include foraging, transiting and breeding. The Guianas Coastal and Caribbean regions encompass important breeding and nesting habitat for leatherbacks [85]. Two adult females in our study nested in the Guianas Coastal region (Trinidad and Costa Rica/Panama) [57], and two adult males remained in coastal waters off Trinidad. In the primarily shallow, shelf Guianas Coastal region, leatherback movements became more sinuous in response to decreasing water depth, probably linked to breeding activity [16], [22], [107] rather than feeding [108], although some leatherbacks do forage during the nesting season [109][111].

The average mixed layer depth varies throughout the tropics, with deepest depths occurring in winter when leatherbacks are present [1]. Leatherbacks occupied the western side of the Tropical Atlantic where average winter mixed layer depths are 70–80 m [1]. This could explain the deeper, longer dives that leatherbacks made there if prey accumulates near the pycnocline and (or) nutricline [112][113]. In the Tropical Atlantic, the night-time depths of diel migrants is consistently within the upper 50–75 m [1], and leatherbacks may target prey aggregated in this layer [106], though nocturnal foraging may be light-limited based on studies of leatherback ocular morphology [114] and feeding behavior [110].

The convergence of the North Equatorial Current, North Equatorial Counter-Current and the North Brazil Current appears to play an important role for overwintering leatherbacks in the southern part of their range. From June to January, the upper North Brazil Current joins the meandering North Equatorial Counter-Current at a retroflection zone near 5–10°N, where large, anti-cyclonic eddies are formed [115]. The eddies, known as North Brazil Current rings, have loop diameters of hundreds of kilometers, with a lifespan on the order of months, and they propagate northwest along the Brazil coast towards the Lesser Antilles [115][116]. North Brazil Current rings, and the convergence of the westward-flowing North Equatorial Current and eastward-flowing North Equatorial Counter-Current at 10–12°N, are associated with enhanced wintertime productivity [1]. Three leatherbacks resided in this region during winter months, associating with meanders from the North Equatorial Current and Counter-Current convergence and mesoscale eddies. Fossette et al. [25] inferred high foraging success year-round for leatherbacks at the southern boundary of the Tropical Atlantic, while a single tagged leatherback resided in this area for several months [117]. Since overall production in open tropical oceans is low compared to temperate and boreal latitudes [1], leatherbacks may rely on enhanced productivity there to maximize foraging opportunities during overwintering periods in the tropics.

Both adult and sub-adult leatherbacks in our study adjusted their movements and dive behavior in response to regional differences in environmental features. Leatherbacks increased their path sinuosity with decreasing water depth in temperate and tropical shelf habitats. This relationship is consistent with increases in gelatinous zooplankton biomass with decreasing water depth [59], and bathymetry may be a key feature in identifying leatherback foraging habitat in neritic regions. Coastal ecosystems are under intense pressure worldwide, with some of the highest predicted cumulative impact in the North American eastern seaboard and the eastern Caribbean [118]. Parts of these regions constituted high-use habitat for leatherbacks in our study, putting turtles at heightened risk from both land- and ocean-based human activity.

Acknowledgments

We thank A. Myers, C. Merigo, C. Innis, M. Dodge, G. Purmont, M. Leach, B. Sharp, S. Landry, M. Murphy, G. Tomasian, N. Fragoso, K. Sampson, R. Smolowitz, K. Hirokawa, J. Casey, S. Leach, J. Wilson, and E. Eldredge for invaluable assistance in the field. We thank A. Rhodin for surgical expertise and guidance in development of humane and effective direct tag attachment procedures for leatherbacks. S. Benson, P. Dutton and M. James gave helpful suggestions and advice regarding leatherback captures at sea. We thank F. Royer for insightful discussions and expert assistance with track reconstruction. M. Coyne provided valuable assistance with oceanographic data extraction. We are grateful to B. Prescott and K. Dourdeville for timely leatherback sightings information in southern New England (www.seaturtlesightings.org). Massachusetts Environmental Police, US Coast Guard, and staff from Provincetown Center for Coastal Studies and New England Aquarium provided logistical support. The authors acknowledge the use of SEATURTLE.ORG's Maptool program (http://www.seaturtle.org/maptool/) and Scott Sherrill-Mix's R-GMT HexMap code (http://scott.sherrillmix.com/R-GMT_HexPlot.php). We thank M. Baumgartner, H. Haas, H. Howell, L. Madin, P. Rago, V. Saba, and two anonymous reviewers for helpful discussions and comments that improved earlier drafts of this manuscript.

Author Contributions

Conceived and designed the experiments: KD ML. Performed the experiments: KD. Analyzed the data: KD BG TM. Wrote the paper: KD BG TM ML.

References

  1. 1. Longhurst AR (2007) Ecological geography of the sea. Second edition. Burlington: Academic Press. 542 p.
  2. 2. Kareiva P, Odell G (1987) Swarms of predators exhibit preytaxis if individual predators use area-restricted search. Am Nat 130(2): 233–270.
  3. 3. Pinaud D, Weimerskirch H (2007) At-sea distribution and scale-dependent foraging behaviour of petrels and albatrosses: a comparative study. J Anim Ecol 76: 9–19.
  4. 4. Robinson PW, Tremblay Y, Crocker DE, Kappes MA, Kuhn CE, et al. (2007) A comparison of indirect measures of feeding behaviour based on ARGOS tracking data. Deep-Sea Res II 54: 356–368.
  5. 5. Weng KC, Foley DG, Ganong JE, Perle C, Shillinger GL, et al. (2008) Migration of an upper trophic level predator, the salmon shark Lamna ditropis, betweet distant ecoregions. Mar Ecol Prog Ser 372: 253–264.
  6. 6. Jonsen ID, Myers RA, James MC (2007) Identifying leatherback turtle foraging behaviour from satellite telemetry using a switching state-space model. Mar Ecol Prog Ser 337: 255–264.
  7. 7. Bailey H, Shillinger G, Palacios D, Bograd S, Spotila J, et al. (2008) Identifying and comparing phases of movement by leatherback turtles using state-space models. J Exp Mar Biol Ecol 356: 128–135.
  8. 8. Bailey H, Benson SR, Shillinger GL, Bograd SJ, Dutton PH, et al. (2012) Identification of distinct movement patterns in Pacific leatherback turtle populations influenced by ocean conditions. Ecol App 22: 735–747.
  9. 9. Carr AF, Ogren LH (1959) The ecology and migrations of sea turtles. 3. Dermochelys in Costa Rica, 1955. Am Mus Novit 1835.
  10. 10. Musick JA, Limpus CJ (1997) Habitat utilization and migration in juvenile sea turtles. In: Lutz PL, Musick JA, editors. The biology of sea turtles. Boca Raton: CRC Press. pp. 137–163.
  11. 11. Plotkin P (2002) Adult migrations and habitat use. In: Lutz PL, Musick JA, Wyneken J, editors. The biology of sea turtles volume II. Boca Raton: CRC Press. pp. 225–233.
  12. 12. Dingle H, Drake VA (2007) What is migration? BioScience 57: 113–121.
  13. 13. Morreale SJ, Standora EA, Spotila JR, Paladino FV (1996) Migration corridor for sea turtles. Nature 384: 319–320.
  14. 14. Hughes GR, Luschi P, Mencacci R, Papi F (1998) The 7000-km journey of a leatherback turtle tracked by satellite. J Exp Mar Biol Ecol 229: 209–217.
  15. 15. Ferraroli S, Georges JY, Gaspar P, Maho YL (2004) Endangered species: where leatherback turtles meet fisheries. Nature 429: 521–522.
  16. 16. James MC, Ottensmeyer CA, Myers RA (2005) Identification of high-use habitat and threats to leatherback sea turtles in northern waters: new directions for conservation. Ecol Lett 8: 195–201.
  17. 17. Shillinger GL, Palacios DM, Bailey H, Bograd SJ, Swithenbank AM, et al. (2008) Persistent leatherback turtle migrations present opportunities for conservation. PLoS Biol 6(7): e171
  18. 18. López-Mendilaharsu M, Rocha CFD, Miller P, Domingo A, Prosdocimi L (2009) Insights on leatherback turtle movements and high use areas in the Southwest Atlantic Ocean. J Exp Mar Biol Ecol 378: 31–39.
  19. 19. Benson SR, Eguchi T, Foley DG, Forney KA, Bailey H, et al. (2011) Large-scale movements and high-use areas of western Pacific leatherback turtles, Dermochelys coriacea. Ecosphere 2(7): art84
  20. 20. Lutcavage ME (1996) Planning your next meal: leatherback travel routes and ocean fronts. In: Keinath JA, Barnard DE, Musick JA, Bell BA, editors. Proc 15th Ann Symp Sea Turtle Biology and Conservation. NOAA Tech Memo NMFS-SEFSC-387, p 355.
  21. 21. Luschi P, Sale A, Mencacci R, Hughes GR, Lutjeharms JRE, et al. (2003) Current transport of leatherback sea turtles (Dermochelys coriacea) in the ocean. Proc R Soc Lond B (Suppl) 270: S129–S132.
  22. 22. Eckert SA (2006) High-use oceanic areas for Atlantic leatherback sea turtles (Dermochelys coriacea) as identified using satellite telemetered location and dive information. Mar Biol 149: 1257–1267.
  23. 23. Gaspar P, Georges JY, Fossette S, Lenoble A, Ferraroli S, et al. (2006) Marine animal behaviour: neglecting ocean currents can lead us up the wrong track. Proc R Soc B 273: 2697–2702.
  24. 24. Doyle TK, Houghton JDR, O'Suilleabhain PF, Hobson VJ, Marnell F, et al. (2008) Leatherback turtles satellite-tagged in European waters. Endang Species Res 4: 23–31.
  25. 25. Fossette S, Hobson VJ, Girard C, Calmettes B, Gaspar P, et al. (2010) Spatio-temporal foraging patterns of a giant zooplanktivore, the leatherback turtle. J Mar Sys 81: 225–234.
  26. 26. Saba VS (2013) Oceanic habits and habitats: Dermochelys coriacea. In: Wyneken J, Musick JA, Lohmann KJ, editors. The biology of sea turtles volume III. Boca Raton: CRC Press. pp. 163–188.
  27. 27. Hays GC, Hobson VJ, Metcalf JD, Righton D, Sims DW (2006) Flexible foraging movements of leatherback turtles across the North Atlantic Ocean. Ecology 87: 2647–2656.
  28. 28. Shoop CR, Kenney RD (1992) Seasonal distributions and abundances of loggerhead and leatherback sea turtles in the waters of the northeastern United States. Herp Monogr 6: 43–67.
  29. 29. Houghton JDR, Doyle TK, Wilson MW, Davenport J, Hays GC (2006) Jellyfish aggregations and leatherback turtle foraging patterns in a temperate coastal environment. Ecology 87: 1967–1972.
  30. 30. Benson SR, Forney KA, Harvey JT, Carretta JV, Dutton PH (2007) Abundance, distribution, and habitat of leatherback turtles (Dermochelys coriacea) off California, 1990-2003. Fish Bull 105: 337–347.
  31. 31. Witt MJ, Broderick AC, Johns DJ, Martin C, Penrose R, et al. (2007) Prey landscapes help identify potential foraging habitats for leatherback turtles in the NE Atlantic. Mar Ecol Prog Ser 337: 231–244.
  32. 32. James MC, Myers RA, Ottensmeyer CA (2005) Behaviour of leatherback sea turtles, Dermochelys coriacea, during the migratory cycle. Proc R Soc B 272: 1547–1555.
  33. 33. Eckert S, Bagley D, Kubis S, Ehrhart L, Johnson C, et al. (2006) Internesting and postnesting movements and foraging habitats of leatherback sea turtles (Dermochelys coriacea) nesting in Florida. Chelonian Conserv Biol 5: 239–248.
  34. 34. Larson RJ (1976) Marine flora and fauna of the northeastern United States. Cnidaria:Scyphozoa. NOAA Tech Rep NMFS Circ 397..
  35. 35. Madin LP, Kremer P, Wiebe PH, Purcell JE, Horgan EH, et al. (2006) Periodic swarms of the salp Salpa aspera in the slope waters of the NE United States: biovolume, vertical migration, grazing, and vertical fluz. Deep-Sea Res 53: 804–819.
  36. 36. Mann KH, Lazier JRN (2006) Dynamics of marine ecosystems: biological-physical interactions in the oceans (3rd ed). Cambridge: Blackwell Science Inc.
  37. 37. Shenker JM (1984) Scyphomedusae in surface waters near the Oregon coast, May-August, 1981. Estuar Coast Shelf Sci 19: 619–632.
  38. 38. Olson DB, Hitchcock GL, Mariano AJ, Ashjian CJ, Peng G, et al. (1994) Life on the edge: marine life and fronts. Oceanography 7: 52–59.
  39. 39. Graham WM, Pages F, Hammer WM (2001) A physical context for gelatinous zooplankton aggregations: a review. Hydrobiologia 451: 199–212.
  40. 40. Purcell JE, Graham WM, Dumont HJ (2001) Jellyfish blooms: ecological and societal importance. Kluwer Academic Publishers.
  41. 41. Deibel D, Paffenhöfer GA (2009) Predictability of patches of neritic salps and doliolids (Tunicata, Thaliacea). J Plankton Res 31: 1571–1579.
  42. 42. Lewison RL, Freeman SA, Crowder LB (2004) Quantifying the effects of fisheries on threatened species: the impact of pelagic longlines on loggerhead and leatherback sea turtles. Ecol Lett 7: 221–231.
  43. 43. Lum LL (2006) Assessment of incidental sea turtle catch in the artisanal gillnet fishery in Trinidad and Tobago, West Indies. Appl Herp 3: 357–368.
  44. 44. Phillips RA, Silk JRD, Croxall JP, Afanasyev V (2006) Year-round distribution of white-chinned petrels from South Georgia: relationships with oceanography and fisheries. Biol Conser 129: 336–347.
  45. 45. Alfaro-Shigueto J, Dutton PH, Van Bressem M, Mangel J (2007) Interactions between leatherback turtles and Peruvian artisanal fisheries. Chelonian Conserv Biol 6: 129–134.
  46. 46. Alfaro-Shigueto J, Mangel JC, Pajuelo M, Dutton PH, Seminoff JA, et al. (2010) Where small can have a large impact: structure and characterization of small-scale fisheries in Peru. Fish Res 106: 8–17.
  47. 47. James MC, Sherrill-Mix SA, Martin K, Myers RA (2006) Canadian waters provide critical foraging habitat for leatherback sea turtles. Biol Conserv 133: 347–357.
  48. 48. McMahon CR, Hays GC (2006) Thermal niche, large-scale movements and implications of climate change for a critically endangered marine vertebrate. Global Change Biol 12: 1330–1338.
  49. 49. Fossette S, Girard C, López-Mendilaharsu M, Miller P, Domingo A, et al. (2010) Atlantic leatherback migratory paths and temporary residence areas. PLoS ONE 5(11): e13908
  50. 50. Shillinger GL, Swithenbank AM, Bailey H, Bograd SJ, Castelton MR, et al. (2011) Vertical and horizontal habitat preferences of post-nesting leatherback turtles in the South Pacific Ocean. Mar Ecol Prog Ser 422: 275–289.
  51. 51. Witt MJ, Bonguno EA, Broderick AC, Coyne MS, Formia A, et al. (2011) Tracking leatherback turtles from the world's largest rookery: assessing threats across the South Atlantic. Proc R Soc B 278: 2338–2347.
  52. 52. Bailey H, Fossette S, Bograd SJ, Shillinger GL, Swithenbank AM, et al. (2012) Movement patterns for a critically endangered species, the leatherback turtle (Dermochelys coriacea), linked to foraging success and population status. PLoS ONE 7(5): e36401
  53. 53. Asper ED (1975) Techniques of live capture of smaller cetacea. J Fish Res Board Can 32: 1191–1196.
  54. 54. Innis C, Merigo C, Dodge K, Tlusty M, Dodge M, et al. (2010) Health evaluation of leatherback turtles (Dermochelys coriacea) in the Northwestern Atlantic during direct capture and fisheries gear disentanglement. Chelonian Conserv Biol 9: 205–222.
  55. 55. Lutcavage M, Rhodin AGJ, Sadove SS, Conroy CR (2001) Direct carapacial attachment of satellite tags using orthopedic bioabsorbable mini-anchor screws on leatherback turtles in Culebra, Puerto Rico. Mar Turt Newsl 95: 9–12.
  56. 56. James MC, Sherrill-Mix SA, Myers RA (2007) Population characteristics and seasonal migrations of leatherback sea turtles at high latitudes. Mar Ecol Prog Ser 337: 245–254.
  57. 57. Dodge KL, Logan JM, Lutcavage ME (2011) Foraging ecology of leatherback sea turtles in the Northwest Atlantic determined through multi-tissue stable isotope analyses. Mar Biol 158: 2813–2824.
  58. 58. Bryant E (2007) 2D location accuracy statistics for Fastloc cores running firmware versions 2.2 and 2.3. Technical Report TR01. Wildtrack Telemetry Systems, Leeds, UK. Available: http://www.wildtracker.com/results_files/Technical%20Report%20TR01.pdf.
  59. 59. Lilley MKS, Beggs SE, Doyle TK, Hobson VJ, Stromberg KHP, et al. (2011) Global patterns of epipelagic gelatinous zooplankton biomass. Mar Biol 158: 2429–2436.
  60. 60. Belkin IM, O'Reilly JE (2009) An algorithm for oceanic front detection in chlorophyll and SST imagery. J Mar Sys 78: 319–326.
  61. 61. Royer F, Lutcavage M (2008) Filtering and interpreting location errors in satellite telemetry of marine animals. J Exp Mar Biol Ecol 359: 1–10.
  62. 62. Batschelet E (1981) Circular statistics in biology. NewYork: Academic Press.
  63. 63. Benhamou S (2004) How to reliably estimate the tortuosity of an animal's path: straightness, sinuosity, or fractal dimension? J Theor Biol 229: 209–220.
  64. 64. R Core Team (2013) R: A language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria. Available: http://www.R-project.org/.
  65. 65. Wessel P, Smith WHF (1991) Free software helps map and display data. EOS Trans AGU 72: 441.
  66. 66. Zuur AF, Ieno EN, Walker NJ, Saveliev AA, Smith GM (2009) Mixed effects models and extensions in ecology with R. NewYork: Springer.
  67. 67. Weimerskirch H, Bonadonna F, Bailleul F, Mabille G, Dell'Omo G, et al. (2002) GPS tracking of foraging albatrosses. Science 295: 1259.
  68. 68. Kuhn CE, Tremblay Y, Ream RR, Gelatt TS (2010) Coupling GPS tracking with dive behavior to examine the relationship between foraging strategy and fine-scale movements of northern fur seal. Endang Species Res 12: 125–139.
  69. 69. McCarthy AL, Heppell S, Royer F, Freitas C, Dellinger T (2010) Identification of likely foraging habitat of pelagic loggerhead sea turtles (Caretta caretta) in the North Atlantic through analysis of telemetry track sinuosity. Prog Oceanogr 86: 224–231.
  70. 70. Chi EM, Reinsel GC (1989) Models for longitudinal data with random effects and AR(1) errors. J Am Stat Assoc 84: 452–459.
  71. 71. Akaike H (1973) Information theory and an extension of the maximum likelihood principle. In: Petrov BN, Csaki F, editors. International Symposium on Information Theory. Budapest: Akademia Kiado, pp. 267–281.
  72. 72. Burnham KP, Anderson DR (1998) Model selection and inference: A practical information-theoretic approach. NewYork: Springer-Verlag. 353 p.
  73. 73. Faes C, Molenbergh G, Aerts M, Verbeke G, Kenward MG (2009) The effective sample size and an alternative small-sample degrees-of-freedom method. The American Statistician 63: 389–399.
  74. 74. Pinheiro JC, Bates DM (2000) Mixed-effects models in S and S-PLUS. NewYork: Springer-Verlag. 528 p.
  75. 75. Jones TT, Bostrom B, Carey M, Imlach B, Mikkelsen J, et al.. (2011) Determining transmitter drag and best-practice attachment procedures for sea turtle biotelemetry. NOAA Tech Memo NMFS-SWFSC-480. 58 p.
  76. 76. Fossette S, Corbel H, Gaspar P, LeMaho Y, Georges JY (2007) An alternative technique for the long-term satellite tracking of leatherback turtles. Endang Species Res 4: 33–41.
  77. 77. Bleakney JS (1965) Reports of marine turtles from New England and Eastern Canada. Can Field Nat 79(2): 120–128.
  78. 78. Heaslip SG, Iverson SJ, Bowen WD, James MC (2012) Jellyfish support high energy intake of leatherback sea turtles (Dermochelys coriacea): video evidence from animal-borne cameras. PLoS ONE 7(3): e33259
  79. 79. Alerstam T, Hedenström A, Åkesson S (2003) Long-distance migration: evolution and determinants. Oikos 103: 247–260.
  80. 80. Schmidt-Nielsen K (1972) Locomotion: energy cost of swimming, flying, and running. Science 177: 222–228.
  81. 81. Jørgensen C, Ernande B, Fiksen Ø, Dieckmann U (2006) The logic of skipped spawning in fish. Can J Fish Aquat Sci 63: 200–211.
  82. 82. Chapman EW, Jørgensen C, Lutcavage ME (2011) Atlantic bluefin tuna (Thunnus thynnus): a state-dependent energy allocation model for growth, maturation, and reproductive investment. Can J Fish Aquat Sci 68: 1934–1951.
  83. 83. Frair W, Ackman RG, Mrosovsky N (1972) Body temperature of Dermochelys coriacea: warm turtle from cold water. Science 177: 791–793.
  84. 84. Paladino FV, O'Connor MP, Spotila JR (1990) Metabolism of leatherback turtles, gigantothermy, and thermoregulation of dinosaurs. Nature 344: 858–860.
  85. 85. TEWG (Turtle Expert Working Group) (2007) An assessment of the leatherback turtle population in the Atlantic Ocean. NOAA Tech Memo NMFS-SEFSC-555.
  86. 86. Walli A, Teo SLH, Boustany A, Farwell CJ, Williams T, et al. (2009) Seasonal movements, aggregations and diving behavior of Atlantic bluefin tuna (Thunnus thynnus) revealed with archival tags. PLoS ONE 4(7): e6151
  87. 87. Whitehead H, Jonsen ID (2013) Inferring animal densities from tracking data using Markov Chains. PLoS ONE 8(4): e60901
  88. 88. Polovina JJ, Balaz GH, Howell EA, Parker DM, Seki MP, et al. (2004) Forage and migration habitat of loggerhead (Caretta caretta) and olive ridley (Lepidochelys olivacea) sea turtles in the central North Pacific. Fish Oceanogr 13: 36–51.
  89. 89. McClellan CM, Read AJ (2007) Complexity and variation in loggerhead sea turtle life history. Biol Lett 3: 592–594.
  90. 90. Potter IF, Galuard B, Howell WH (2011) Horizontal movement of ocean sunfish, Mola mola, in the northwest Atlantic. Mar Biol 158: 531–540.
  91. 91. Skomal GB, Zeeman SI, Chisholm JH, Summers EL, Walsh HJ, et al. (2009) Transequatorial migrations by basking sharks in the western Atlantic Ocean. Curr Biol 19: 1019–1022.
  92. 92. Sibert JR, Lutcavage ME, Nielsen A, Brill RW, Wilson SG (2006) Inter-annual variation in large-scale movement of Atlantic bluefin tuna (Thunnus thynnus) determined from pop-up satellite archival tags. Can J Fish Aquat Sci 63: 2154–2166
  93. 93. Galuardi B, Royer F, Golet W, Logan J, Neilson J, et al. (2010) Complex migration routes of Atlantic bluefin tuna (Thunnus thynnus) question current population structure paradigm. Can J Fish Aquat Sci 67: 966–976.
  94. 94. Neilson JD, Smith S, Royer F, Paul SD, Porter JM, et al.. (2009) Investigations of horizontal movements of Atlantic swordfish using pop-up satellite archival tags. Tagging and Tracking of Marine Animals with Electronic Devices. Available: http//dx.doi.org/10.1007/978-1-4020-9640-2_18
  95. 95. Reina RD, Abernathy KJ, Marshall GJ, Spotila JR (2005) Respiratory frequency, dive behaviour and social interactions of leatherback turtles, Dermochelys coriacea during the inter-nesting interval. J Exp Mar Biol Ecol 316: 1–16.
  96. 96. Sherman K, Skjoldal HR (2002) Large marine ecosystems of the North Atlantic: changing states and sustainability. Amsterdam: Elsevier Science.
  97. 97. Saba VS, Spotila JR, Chavez FP, Musick JA (2008) Bottom-up and climatic forcing on the worldwide population of leatherback turtles. Ecology 89: 1414–1427.
  98. 98. Lazell JD (1980) New England waters: critical habitat for marine turtles. Copeia 1980: 290–295.
  99. 99. Schick RS, Goldstein J, Lutcavage ME (2004) Bluefin tuna (Thunnus thynnus) distribution in relation to sea surface temperature fronts in the Gulf of Maine (1994-96). Fish Oceanogr 13: 225–238.
  100. 100. Bost CA, Cotté C, Bailleul F, Cherel Y, Charrassin JB, et al. (2009) The importance of oceanographic fronts to marine birds and mammals of the southern oceans. J Mar Sys 78: 363–376.
  101. 101. Polovina JJ, Kobayashi DR, Parker DM, Seki JP, Balaz GH (2000) Turtles on the edge: movements of loggerhead turtles (Caretta caretta) along oceanic fronts, spanning longline fishing grounds in the central North Pacific, 1997-1998. Fish Oceanogr 9: 71–82.
  102. 102. Seminoff JA, Zárate P, Coyne M, Foley DG, Parker D, et al. (2008) Post-nesting migrations of Galápagos green turtles Chelonia mydas in relation to oceanographic conditions: integrating satellite telemetry with remotely sensed ocean data. Endang Species Res 4: 57–72.
  103. 103. Belkin IM, Cornillon PC, Sherman K (2009) Fronts in large marine ecosystems. Prog Oceanogr 81: 223–236.
  104. 104. Graham WM, Martin DL, Martin JC (2003) In situ quantification and analysis of large jellyfish using a novel video profiler. Mar Ecol Prog Ser 254: 129–140.
  105. 105. Rakow KC, Graham WM (2006) Orientation and swimming mechanics by the scyphomedusa Aurelia sp. in shear flow. Limnol Oceanogr 51: 1097–1106.
  106. 106. James MC, Ottensmeyer CA, Eckert SA, Myers RA (2006) Changes in diel diving patterns accompany shifts between northern foraging and southward migration in leatherback turtles. Can J Zool 84: 754–765.
  107. 107. James MC, Eckert SA, Myers RA (2005) Migratory and reproductive movements of male leatherback turtles (Dermochelys coriacea). Mar Biol 147: 845–853.
  108. 108. Plot V, Jenkins T, Robin JP, Fossette S, Georges JY (2013) Leatherback turtles are capital breeders: morphometric and physiological evidence from longitudinal monitoring. Physiol Biochem Zool 86: 385–397.
  109. 109. Myers AE, Hays GC (2006) Do leatherback turtles Dermochelys coriacea forage during the breeding season? A combination of data-logging devices provide new insights. Mar Ecol Prog Ser 322: 259–267.
  110. 110. Casey J, Garner J, Garner S, Williard Southwood A (2010) Diel foraging behavior of gravid leatherback sea turtles in deep waters of the Caribbean Sea. J Exp Biol 213: 3961–3971.
  111. 111. Fossette S, Gleiss AC, Casey JP, Lewis AR, Hays GC (2012) Does prey size matter? Novel observations of feeding in the leatherback turtle (Dermochelys coriacea) allow a test of predator-prey size relationships. Biol Lett 8: 351–354.
  112. 112. Eckert SA, Eckert KL, Ponganis P, Kooyman GL (1989) Diving and foraging behavior of leatherback sea turtles (Dermochelys coriacea). Can J Zool 67: 2834–2840.
  113. 113. Hays GC, Houghton DJ, Issacs C, King RS, Lloyd C, et al. (2004) First records of oceanic dive profiles of leatherback turtles, Dermochelys coriacea, indicate behavioral plasticity associated with long-distance migration. Anim Behav 67: 733–743.
  114. 114. Brudenall DK, Schwab IR, Fritsches KA (2008) Ocular morphology of the Leatherback sea turtle (Dermochelys coriacea). Vet Ophthalmol 11: 99–110.
  115. 115. Johns WE, Lee TN, Schott FA, Zantopp RJ, Evans RH (1990) The North Brazil current retroflection: seasonal structure and eddy variability. J Geophys Res 95: 22103–22120.
  116. 116. Fratantoni DM, Richardson PL (2006) The evolution and demise of North Brazil current rings. J Phys Oceanogr 36: 1241–1264.
  117. 117. Hays GC, Houghton JD, Myers AE (2004) Pan-Atlantic leatherback turtle movements. Nature 429: 522.
  118. 118. Halpern BS, Walbridge S, Selkoe KA, Kappel CV, Micheli F, et al. (2008) A global map of human impact on marine ecosystems. Science 319: 948–952.