Predicting the distribution of oceanic-stage Kemp's ridley sea turtles

Nathan F. Putman, Katherine L. Mansfield, Ruoying He, Donna J. Shaver, Philippe Verley

Abstract

The inaccessibility of open ocean habitat and the cryptic nature of small animals are fundamental problems when assessing the distribution of oceanic-stage sea turtles and other marine animals sharing similar life-history traits. Most methods that estimate patterns of abundance cannot be applied in situations that are extremely data limited. Here, we use a movement ecology framework to generate the first predicted distributions for the oceanic stage of the Kemp's ridley sea turtle (Lepidochelys kempii). Our simulations of particle dispersal within ocean circulation models reveal substantial annual variation in distribution and survival among simulated cohorts. Such techniques can help prioritize areas for conservation, and supply inputs for more realistic demographic models attempting to characterize population trends.

1. Introduction

Knowing the distribution of a species and its variability through time is fundamental for successful species conservation and management [1]. The distribution of a species is commonly predicted by correlating occurrence records with environmental metrics to generate maps of habitat suitability or likelihood of occurrence [1,2]. Alternatively, distribution can be estimated from principles of movement ecology theory: the distribution of a species is mechanistically predicted by simulating the movement process of individuals within a realistic environmental model [3]. Compared with more frequently used tools for predicting distribution, for instance ecological niche models [1,2], the movement ecology approach may be particularly useful for species that occupy habitats which preclude sampling or that possess cryptic life-stages, for instance sea turtles [4,5]. In most sea turtle species, the young migrate from beaches into the ocean and are rarely encountered again until they return to coastal waters 2–15 years later [46]. However, predicting the distribution of oceanic-stage turtles is possible using particle-tracking software and the output of high-resolution ocean circulation models [7]. Nesting beaches, which serve as initiation points for simulations, are well defined for most species. Although swimming behaviour can affect sea turtle distributions [8], young turtles are relatively weak swimmers [6]; thus, as a first approximation, predictions of early oceanic-stage sea turtle distributions can be generated by simulating particle dispersal from nesting beaches based solely on ocean currents [412].

We use this approach to predict the distribution of the Kemp's ridley sea turtle (Lepidochelys kempii) across the Gulf of Mexico and to estimate their early survival. This critically endangered turtle nests almost exclusively in the western Gulf of Mexico, with concentrated nesting occurring in the vicinity of Ranch Nuevo in Tamaulipas, Mexico [13]. As such, the early life history of this species is probably constrained to the Gulf of Mexico [5]. The limited geographical range of these turtles combined with high-resolution ocean circulation models for this region provide a unique opportunity to examine how physical processes at the ocean surface influence the distribution of a life-stage in sea turtles that has been, thus far, impossible to estimate [4,5].

2. Material and methods

To simulate the movement of young (less than 2 years old) Kemp's ridleys, we extracted surface currents from the Gulf of Mexico Global Hybrid Coordinate Ocean Model (GOM HYCOM) [14]; GOM HYCOM output has a spatial resolution of 0.04° (approx. 4.5 km), daily snapshots of velocity and resolves mesoscale processes such as meandering currents, fronts, filaments and eddies [14]. Ichthyop (v. 2.21) particle-tracking software [15] calculated trajectories of virtual particles released within 0.04° × 0.04° cells centred 10 km offshore of nesting sites (table 1 and electronic supplementary material). Locations of hatchling releases were obtained from records for the 2009, 2010 and 2011 Kemp's ridley cohorts (Donna J. Shaver, Padre Island National Seashore, Division of Sea Turtle Science and Recovery; Patrick M. Burchfield and Luis J. Peña, Gladys Porter Zoo; Raul J. Gonzalez, Acuario de Veracruz, Mexico and Rosa C. Martinez Portugal, CONANP de Veracruz, Mexico 2012, personal communication). We used the mean relative abundance of hatchlings released at each location during the months of June, July and August to weight the proportion of particles that were released daily from each location for a particular month (table 1). Simulations were performed for the 2003–2010 cohorts, and particles were tracked for 515 days. ICHTHYOP implemented a Runge–Kutta fourth-order time-stepping method, whereby particle position was calculated every 30 min. Particles that encountered a coastline ‘bounced’ along the coast until currents changed to move them offshore [16]. To simulate the ‘frenzy-period’ of hatchling turtles, particles swam at 0.25 m s−1 in an offshore direction (±20°) during the first 48 h [5,16]; a maximum distance of 43.2 km in still water. The duration of the ‘frenzy-period’ is not described for Kemp's ridleys, this approximation is based on what is known in loggerheads (Caretta caretta) [16]. We determined distribution patterns for a cohort by summing the number of particles within each grid cell at daily intervals throughout the simulation. This highlighted locations that particles encountered, were retained, and, therefore, where Kemp's ridleys likely occur. Age structure across the Gulf was estimated by summing the ages of particles within a cell and dividing by the total number of particles counted within that cell as described above.

View this table:
Table 1.

Parameters of simulated hatchling release.

To assess the probability of eastward transport, we recorded the proportion of particles that crossed different longitudes (90° W, 85° W and 80° W). We also simulated mortality of particles based on water depths encountered during the first 185 days. Because shallower areas have greater predation risk [5], we specified increasing survival probabilities with increasing depth. At the end of each day, the likelihood of survival for particles in water less than 10 m = 98%, 10–50 m = 99%, 50–100 m = 99.5%, 100–200 m = 99.75% and greater than 200 m = 100% [16]. The actual likelihood of predation is not known; this simplification of mortality risk is meant for illustrative purposes and is therefore not represented in predicted distributions (figures 1 and 2).

Figure 1.

Predicted abundance of oceanic-stage Kemp's ridley turtles in the Gulf of Mexico by cohort. White squares indicate particle release sites. Coloration is scaled logarithmically and indicates the number of particles within each grid cell throughout the 515-day tracking period for the (a) 2003, (b) 2004, (c) 2005, (d) 2006, (e) 2007, (f) 2008, (g) 2009 and (h) 2010 cohorts.

Figure 2.

Predicted age structure of oceanic-stage Kemp's ridley turtles in the Gulf of Mexico by cohort. White squares indicate particle release sites. Coloration indicates the mean age of particles within each grid cell throughout the 515-day tracking period for the (a) 2003, (b) 2004, (c) 2005, (d) 2006, (e) 2007, (f) 2008, (g) 2009 and (h) 2010 cohorts.

Transport predictions are more accurate with finer spatio-temporal resolution of ocean velocity fields [7] and with greater numbers of particles released [17]. We examined the sensitivity of our estimates of eastward transport and survival to these parameters using hindcast output from the South Atlantic Bight–Gulf of Mexico ocean circulation model (SABGOM) [18]. SABGOM has a spatial resolution of 5 km and hourly snapshots of ocean velocity fields from 1 June 2010 to 31 December 2010. We released approximately 33 times as many particles (in the same proportion) from nesting sites (table 1) and calculated their trajectories through 31 December 2010. Finally, we also modified the behaviour of particles encountering coastlines, rather than ‘bouncing’ (which might not depict actual near shore currents) particles in SABGOM simply ‘beached’.

3. Results

Simulations indicate that oceanic-stage Kemp's ridley turtles are likely to be distributed throughout the Gulf of Mexico and into the northwestern Atlantic (figure 1). Highest abundance was predicted in the western Gulf; more than half of all particles for each cohort remained west of 90° W for the entirety of the simulation (table 2). Coastal regions from southern Texas, USA to Tabasco, Mexico consistently had high abundance. Some cohorts also contributed substantially to the coastal waters of southwest Florida, USA and the Florida Keys (figure 1). Mean particle age generally increased from west to east across the Gulf of Mexico (figure 2). At times, however, rapid transport eastward and out of the Gulf occurred (figure 2d and table 2). Variation in predicted distributions resulted in differences in cohort survival spanning an order of magnitude (table 2). Survival and transport predictions appear robust with respect to the modelling parameters that we varied (table 2).

View this table:
Table 2.

Percentage of particles transported eastward and simulated survival.

4. Discussion

Given the limited data on oceanic-stage turtles, these predicted distributions fill a substantial gap regarding sea turtle ecology. These simulations identified locations of possible conservation importance and where field-based research efforts could be focused. The waters offshore of Tamaulipas, Mexico may serve as a nursery area for numerous Kemp's ridleys less than 1 year old (figures 1 and 2). Older turtles, nearing the age when recruitment to near shore waters occurs, are more likely to be distributed in the northern Gulf, eastern Gulf and western Atlantic (figure 2). Anthropogenic stressors to these regions, such as oil spills and incidental take in fisheries, may be particularly detrimental to the Kemp's ridley's population growth [14]. By contrast, young Kemp's ridleys are less likely to be found in the Caribbean Sea and on the Campeche Bank; potential threats in these areas might be less problematic to this life stage (figure 1).

Although our modelling assumptions necessarily simplify the behaviour of young turtles, the predictions are consistent with available in situ observations. Proposed hatchling dispersal trajectories based on surface drifters [5] were observed within our models (table 2). Likewise, sporadic reports of numerous small juveniles stranding along the Texas, USA coast [19] or caught in the waters of Alabama, USA [4] imply that the number of Kemp's ridleys entering these regions varies substantially, as predicted by our simulations (figure 1). Witherington et al. captured 38 Kemp's ridleys in the pelagic Saragassum community on the West Florida Shelf. These turtles had a mean straight carapace length of 233 mm (range: 175–276 mm) corresponding to approximately 1–2 years of age [6]. Our simulations suggest that turtles could reach these waters in less than a year—with the average age of particles around 1.2 years. Particles (or turtles) entering the West Florida Shelf would probably be retained by surface currents [16,20], which might facilitate turtles' ontogenetic shift to coastal foraging grounds in this region [20].

Our analyses imply that variation in the oceanic-stage distribution among cohorts (driven by variation in ocean circulation) strongly affects predicted survival and may, therefore, profoundly influence population trends of this species. Demographic models typically assume that natural mortality is constant from year to year for oceanic-stage Kemp's ridleys; variations in population trends are often attributed to anthropogenic factors [13]. Using spatially explicit and temporally variable estimates of survival as inputs for demographic models could help quantify the role of anthropogenic and environmental factors influencing population abundance. Further development of these models in tandem with empirical data on sea turtle ecology (e.g. swimming and diving behaviour, natural and anthropogenic mortality) hold considerable promise in aiding the conservation of this and other species that share similar life-history traits.

Funding statement

K.L.M. was supported by the NOAA Oil Spill Supplemental Spend Plan. R.H. was supported by NASA grant no. 09-IDS09-0040, NOAA grant no. IOOS-2011-2002515 and GRI GISR grant no. 12-09/GoMRI-006.

Acknowledgements

Padre Island National Seashore, Division of Sea Turtle Science and Recovery; Patrick Burchfield and Jaime Peña, Gladys Porter Zoo; Raul de J. Gonzalez Diaz Miro, Acuario de Veracruz and Rosa Ciria Martinez Portugal, CONAP de Veracruz provided hatchling production information.

  • Received April 15, 2013.
  • Accepted July 19, 2013.

References

View Abstract