Identifying patterns in foraging-area origins in breeding aggregations of migratory species: Loggerhead turtles in the Northwest Atlantic

Population assessments conducted at reproductive sites of migratory species necessitate understanding the foraging-area origins of breeding individuals. Without this information, efforts to contextualize changes in breeding populations and develop effective management strategies are compromised. We used stable isotope analysis of tissue samples collected from loggerhead sea turtles (Caretta caretta) nesting at seven sites in the Northern Recovery Unit (NRU) of the eastern United States (North Carolina, South Carolina and Georgia) to assign females to three separate foraging areas in the Northwest Atlantic Ocean (NWA). We found that the majority of the females at NRU nesting sites (84.4%) use more northern foraging areas in the Mid-Atlantic Bight, while fewer females use more proximate foraging areas in the South Atlantic Bight (13.4%) and more southerly foraging areas in the Subtropical Northwest Atlantic (2.2%). We did not find significant latitudinal or temporal trends in the proportions of NRU females originating from different foraging areas. Combining these findings with previous data from stable isotope and satellite tracking studies across NWA nesting sites showed that variation in the proportion of adult loggerheads originating from different foraging areas is primarily related differences between recovery units: individuals in the NRU primarily use the Mid-Atlantic Bight foraging area, while individuals from the three Florida recovery units primarily use the Subtropical Northwest Atlantic and Eastern Gulf of Mexico foraging areas. Because each foraging area is associated with its own distinct ecological characteristics, environmental fluctuations and anthropogenic threats that affect the abundance and productivity of individuals at nesting sites, this information is critical for accurately evaluating population trends and developing effective region-specific management strategies.

affect the abundance and productivity of individuals at nesting sites, this information is critical for accurately evaluating population trends and developing effective region-specific management strategies.
in lower trophic levels, which are then transferred up the food web and incorporated into the tissues of turtles using each area [20]. Because isotopic turnover is relatively slow in sea turtles (at least 4 months in epidermal tissue) [21], skin samples collected at nesting sites reflect the diet of turtles at foraging areas prior to their breeding migration [14][15]18,[22][23]. Similar to skin samples, the isotopic composition of egg yolks collected at nesting beaches can be used to assign nesting females to distinct foraging areas [24][25][26][27][28].
Previous studies on NRU loggerheads show that the majority of turtles (70-80%) nesting on Bald Head Island, NC and Wassaw Island in Georgia (GA) use the MAB foraging area, while fewer use the SAB foraging area (15-25%) and even fewer use the SNWA foraging area (~5%) [14][15]29]. The first goal of this study was to evaluate the consistency of this pattern across the nesting area used by NRU loggerheads, including nesting sites located within the MAB foraging area (north of Cape Hatteras, NC). The second goal of this study was to combine these new data with previous data from stable isotope and satellite tracking studies across NRU and NWA nesting and breeding sites to gain a deeper understanding of the spatial patterns in habitat use and migration of this important and threatened population.

Ethics statement
The animal use protocol for this research was reviewed and approved by the University of Florida Institutional Animal Care and Use Committee (IACUC protocol #20081985 and #201101985). Procedures were permitted by the individual state sea turtle management agencies under the authority of the United States Fish and Wildlife Service (North Carolina Wildlife Resources Commission, South Carolina Department of Natural Resources, and Georgia Department of Natural Resources).

Study sites and sample collection
Egg yolk samples were collected from 596 clutches laid by loggerhead turtles during the 2012 and 2013 nesting seasons (May-August) at six nesting sites (  Fig 1). Egg yolk samples were collected within 10-12h of clutch deposition and placed in 95% ethanol, while skin samples (6mm biopsy punches) were collected from the shoulder region of each turtle following oviposition or failed nesting and placed in 70% ethanol. All samples were stored at room temperature prior to sample preparation and analysis in January 2014 (6-19 months after collection).
At the two southernmost sampling sites, Wassaw and Jekyll Island, GA (Fig 1), individuals encountered during nocturnal patrols were identified using either physical tags (two Inconel metal tags and one PIT tag) or unique multilocus genetic tags determined from DNA extracted from skin samples collected during each turtle's first encounter in the study [30][31]. At these two sites, tissue samples-yolk from WAS or epidermis from JEK-were collected from each female one time per season (i.e., samples were not collected during within-season recaptures of previously sampled females identified by their tags). At the five sites in SC and NC where nocturnal tagging patrols are not conducted (ORE, PEA, STH, KWH and HHI), yolk samples were identified to individual females using the aforementioned multilocus genetic tags extracted from maternal DNA found in the shell of freshly laid eggs [32]. See Shamblin et al. [32][33] for details on genetic tagging across this subpopulation. All samples were included in subsequent stable isotope analysis (SIA) to crosscheck foraging-ground assignments when individual females were sampled at more than one site within the same season or within more (B) Previously published foraging-area assignments and proportions for 273 skin samples collected from NRU loggerheads nesting at one additional site (BHI) [14] and six additional years from WAS [15,29]. (C) Previously published foraging-area assignments and proportions for 310 skin samples collected from loggerhead turtles nesting at four sites in the Peninsular Florida Recovery Unit (PFRU; listed north to south) [14].
Abbreviations: MAB, Mid-Atlantic Bight; SAB, South-Atlantic Bight; SNWA, Subtropical Northwest Atlantic (Fig 1) a Samples collected from the same individual, either at different sites within the same season or in different seasons, and assigned to one of the three foraging areas with posterior probabilities �0.80 were excluded from analyses. All duplicate samples from the same individual were assigned to the same foraging area.
https://doi.org/10.1371/journal.pone.0231325.t001 than one season, but each individual was included only once when comparing the contribution of different foraging grounds to each nesting site.

Sample preparation and analysis
Ethanol-preserved yolk samples were mixed into a homogeneous solution, then portions of each sample were pipetted into weigh boats and dried at 60˚C for 24-48h. Yolk samples were not lipid-extracted prior to SIA. Instead, previous work demonstrates that untreated yolk beaches) that were included because foraging areas were assigned using the same training data and statistical constraints [14]. (B) There were no significant trends between the latitude of NRU nesting sites and the proportion of females originating from each of the three different foraging areas (simple linear regression; one analysis for each of the three foraging areas). preserved in ethanol can be converted to maternal epidermis values using the following equations: δ 13 C skin = 0.58 x δ 13 C yolk -4.27 and δ 15 N skin = 1.16 x δ 15 N yolk -1.90 [27]. These equations account for the lipid normalization step, as there are separate equations for yolk that was frozen or lipid-extracted prior to SIA [27].
Skin samples were rinsed with distilled water and cleaned with isopropyl alcohol swabs, then the epidermis was separated from the dermis, homogenized using a scalpel blade and dried at 60˚C for 24-48h. Skin samples were not lipid extracted prior to SIA because lipid extraction does not affect the stable isotope values of loggerhead epidermis [15].
For SIA, 0.5-0.6 mg of each sample was analyzed for δ 13 C and δ 15 N composition by combustion in an ECS 4010 elemental analyzer (Costech) interfaced via a ConFlo III to a DeltaPlus XL isotope ratio mass spectrometer (ThermoFisher Scientific) at the University of Florida Light Stable Isotope Mass Spectroscopy Laboratory. Delta notation is used to express stable isotope abundances, as defined as parts per thousand (‰) relative to the standard: δX = [(R sample /R standard )-1] × 1000, where R sample and R standard are the corresponding ratios of rare to common isotopes ( 13 C/ 12 C and 15 N/ 14 N) in the sample and international standard, respectively. Vienna Pee Dee Belemnite was used as the standard for 13 C, and atmospheric N 2 for 15 N. Working standard, L-glutamic acid USGS40 (δ 13 C = −26.39‰ and δ 15 N = −4.52‰), was used to normalize results. In addition, a reference laboratory standard, homogenized loggerhead scute (δ 13 C = −18.36‰ and δ 15 N = 7.68‰), was used to examine consistency in isotopic values in a sample similar to the samples used in this study. The analytical precision of measurements-calculated as the standard deviation of replicates from laboratory standards-was 0.14‰ for δ 13 C and 0.24‰ for δ 15 N of Lglutamic acid USGS40 (N = 77) and 0.20‰ for δ 13 C and 0.26‰ for δ 15 N of loggerhead scute (N = 29).  [16,22], Pajuelo et al. [14], Price et al. [29], Vander Zanden et al. [15], and this study (S3 Table). Base map generated using the SEATURTLE.ORG Maptool [17]. https://doi.org/10.1371/journal.pone.0231325.g004

Foraging area determination
Quadratic discriminant function analysis of stable isotope values was used to assign females to one of the three foraging areas known for NRU loggerheads: (1) MAB, (2) SAB, or (3) SNWA (Fig 1). The training data for the quadratic discriminant function analysis consisted of skin stable isotope values of 60 adult loggerheads from known foraging areas (either through satellite tracking or from sample collection at a foraging site) collected between 2004 and 2011 [14; S1 Table]. A uniform prior probability distribution was assumed, meaning it was equally likely the individuals could have originated from each foraging area, rather than weighting the prior probability by sample size from each area. Only females assigned to one of the three foraging areas with posterior probabilities �0.80 were considered. This constraint translates to an eightfold improvement in odds over random assignment [34]. Using stable isotope values to determine foraging origin of animals depends on the use of isotopic measurements at temporal and spatial scales relevant to the movement of the animals of interest [35]. Several factors may affect the isotopic values of sea turtles over time (e.g., baseline variation, prey composition, food stress, etc.), which could confound foraging area determination. The effect of these factors on the isotopic values turtle tissues are not yet well understood or have not been described. However, despite these potential difficulties, the highly structured and predictable isotopic values encountered in the turtles with known foraging areas sampled over various years allowed for foraging area determination of sea turtles in the NWA in this study.

Statistical analyses
Pearson's Chi-squared tests with p-values computed through Monte Carlo simulation (based on 2,000 replicates) were used to compare the proportion of turtles originating from different foraging areas. First, we compared the inter-annual variation within each applicable nesting site in this study: 2012-2013 for PEA, STH, KWH, HHI and WAS, and 2011-2013 for JEK (Bonferroni-corrected alpha = 0.008 for six comparisons; ORE sampled only in 2012). Second, we pooled all years within each site and tested for differences between sites by making all pairwise comparisons between the seven sites in this study, Bald Head Island (BHI), NC (2004)(2005), and an eighth "site" called WAS � that included additional data for WAS collected between 2004-2014 (Table 1B) (Bonferroni-corrected alpha = 0.0014 for 34 comparisons). Data for BHI and WAS published in Pajuelo et al. [14], Vander Zanden et al. [15], and Price et al. [29] provide direct comparisons because these studies use the same training data and statistical constraints to assign females to different foraging areas. From these data, we also tested for trends between the latitude of each nesting site and the proportion of turtles originating from each foraging area using simple linear regressions-one analysis for each of the three foraging areas. Last, we compared the inter-annual variation across years on WAS and tested for differences between years by making all pairwise comparisons between sampled years from 2004-2014 (Bonferroni-corrected alpha = 0.002 for 28 comparisons). From these data, we also estimated trends in the proportion of turtles originating from each foraging area over time using generalized least squares with first order autocorrelated error structure [36]-one analysis for each of the three foraging areas.

Review of stable isotope and satellite tracking studies
We conducted a two-tiered literature search to collect published articles that use either stable isotopes or satellite tracking to identify the foraging-area origins/destinations of nesting (females) or breeding (males) loggerhead turtles in the NWA. First, a structured search was conducted in Google Scholar and Web of Science using the following search terms: stable isotope, satellite, loggerhead, Caretta and Atlantic. Then, an unstructured search was conducted by reviewing the reference lists of all relevant publications from the structured search. Theses and dissertations were included, but conference presentations and reports were not. We reviewed all potential references and retained only those involving loggerheads nesting or breeding in the NWA that were either assigned to foraging areas using stable isotope analysis or satellite tracked to foraging areas during post-breeding migrations.
For each nesting or breeding site, we collated the number of individuals that were assigned or tracked to each foraging area (Fig 1), keeping stable isotope and satellite tracking studies separate. We combined the new dataset in this study with other studies that use stable isotopes. For satellite-tracking studies, we counted each satellite-tracked individual only once (to the best of our knowledge). When satellite tracking data were included in multiple publications, we used the foraging-area destinations from the most recent study because they frequently involve different and/or more robust statistical analyses, as well as additional unpublished tracks collected in subsequent years.

New foraging-area assignments in the NRU
Of the 596 yolk samples (ORE, PEA, STH, KWH, HHI and WAS) and 150 skin samples (JEK) collected from individual females and analyzed for stable isotopes, 493 yolk samples (83%) and 120 skin samples (80%) were assigned to one of the three foraging areas with posterior probabilities �0.80 (Table 1A; S1 Fig). A total of 133 samples (18%) could not be assigned to a foraging area-posterior probability <0.80 for all foraging areas. Twenty-five individuals were sampled twice either at more than one site within the same season or within more than one season: 15 individuals had both samples assigned to a foraging area (posterior probability �0.80 for one area), six had only one sample assigned, and four had neither sample assigned (posterior probability <0.80 for all areas). All 15 individuals with assignable duplicate samples were assigned to the same foraging area in both samples (14 MAB and 1 SNWA). After removing unassignable duplicates and assignable duplicates collected during a second sampling occasion, 482 yolk samples (81%) and 116 skin samples (77%) were attributed to 598 individuals, each assigned to one of the three foraging areas with posterior probabilities �0.80 and included only once in the dataset. Stable isotope ratios (following appropriate egg yolk to skin conversions) and posterior probabilities for foraging-area assignment of each sample are provided in S2 Table (metadata for S2 Table provided in S1 File).
Within each nesting site sampled in this study, there were no annual differences in the pro- Overall, the majority of the assigned females nesting at the seven study sites was assigned to the northern MAB foraging area (505 out of 598 females; 84.4%) ( Table 1). Eighty females (13.4%) were assigned to the central SAB foraging area, and 13 females (2.2%) were assigned to the southern SNWA foraging area. These proportions were qualitatively similar to previous studies conducted on NRU loggerheads in different years, but very different than loggerheads nesting in Peninsular Florida (Table 1; Fig 2). Florida data from Pajuelo et al. [14] are shown in Table 1C and Fig 2 because the same training data and statistical constraint (posterior probability �0.80) were used to assign females to different foraging areas. Ceriani et al. [16,22] used different training data and a less stringent statistical constraint (posterior probability �0.60) for foraging-area assignments and are therefore not included in Table 1 or Fig 2 but are included in the stable isotope summary below.
Among NRU nesting sites, the proportion of females using different foraging areas was significantly different between (1) JEK and KWH and (2) JEK and WAS (dark gray cells, Table 2; Fig 2A). While the relative proportions remained the same, these differences are primarily attributable to a slightly higher proportion of females nesting on JEK (the southernmost sampling site) that use the central SAB foraging area versus the more northerly MAB foraging area. Despite some significant differences between sites, as well as smaller samples sizes at the two northern sites (ORE and PEA), there were no significant trends between the latitude of the nesting site and the proportion of females originating from each of the three different foraging areas (simple linear regression; MAB: y = -0.034x + 1.94, F 6 = 3.32, P = 0.12; SAB: y = 0.026x -0.69, F 6 = 2.44, P = 0.17; SNWA: y = 0.009x -0.27, F 6 = 6.30, P = 0.05; Fig 2B).
Among the eight sampled nesting seasons between 2004-2014 on WAS, there was a significant difference in the proportion of females using different foraging areas (Pearson's chisquare tests; χ 2 = 36.25, df = 14, P = 0.0035). The proportion of females using different foraging areas was significantly different between (1) Table 3; Fig 3A). Despite some differences among years, as well as small samples sizes in two of the early years (2004 and 2006), there were no significant trends over time in the proportion of WAS females originating from the three foraging areas (generalized least squares: MAB, t 8 = 1.39, p-value = 0.22; SAB, t 8 = -1.16, p-value = 0.29; SNWA, t 8 = -1.33, p-value = 0.23; Fig 3B).

Patterns within the NRU and across the NWA
Among 24 published articles identified in our literature search (S3 Table), we found 1,398 loggerheads from 13 nesting sites that were assigned to an NWA foraging area using stable isotope analysis (Fig 4) and 408 loggerheads from 16 nesting/breeding sites that were tracked to an NWA foraging area using satellite telemetry (Fig 5). Spatial patterns in the proportion of individuals that use different foraging areas were mostly congruent between stable isotope and satellite tracking studies. Stable isotopes have not been used to infer foraging-area origins for female loggerheads nesting in Dry Tortugas or Gulf of Mexico (Fig 4), making satellite tracking data the only means to compare patterns across all recovery units in the US. Foraging-area proportions from nesting sites with small samples sizes (smaller pie charts on Figs 4 and 5) are less robust because the foraging-area origin of each individual carries proportionately more weight for describing the overall pattern. Because satellite-tracking studies often suffer from small sample sizes, these patterns should be interpreted with more caution.

Patterns within the NRU
Our data highlight the importance of the MAB foraging area for the long-term productivity and stability of the NRU loggerhead subpopulation. Variation in the proportion of females assigned to different foraging areas was subtle based on isotopic data with the majority of females at all eight nesting sites originating from the northern MAB foraging area (Fig 4). The importance of the MAB foraging area to NRU loggerheads is also reflected in the post-nesting migrations of satellite-tracked females (Fig 5). Whereas we detected some differences among NRU nesting sites and among nesting seasons at one long-term NRU site (Wassaw Island, GA) in our new dataset, these differences were not associated with significant latitudinal or temporal trends in the proportions of females originating from different foraging areas (Figs 3B and 4B). However, in the satellite-tracking data there was a significant latitudinal trend in  [37][38], Ceriani et al. [9,16,22], Dodd and Byles [39], Evans et al. [40], Foley et al. [7], Girard et al. [41], Griffin et al. [42], Hardy et al. [43], Hart et al. [12,[44][45], Hawkes et al. [6], Lamont et al. [46], Mansfield [47], Phillips [48], Plotkin and Spotila [49], Tucker et al. [8], Vander Zanden et al. [50] (S3 Table). Base map generated using the SEATURTLE.ORG Maptool [17].
https://doi.org/10.1371/journal.pone.0231325.g005 which the proportion of post-nesting females that migrate to the MAB vs. SAB and SNWA decreased from south to north. These differences might be attributable to a subtle, long-term increase in the proportion of NRU females using the MAB foraging area that has swamped any latitudinal pattern in more recent years.
Satellite tracking data that were used to determine the post-nesting foraging-area destinations in Fig 5 were collected between 1992-2008 (79% before 2005), whereas studies that used isotopic data to assign nesting females to foraging areas (Fig 4) were conducted between 2004-2014 (72% after 2012) (S3 Table). Although the trend was not statistically significant and early years may be biased by small sample sizes, the proportion of females nesting on Wassaw Island that were assigned to the MAB shows a similar increasing pattern over time (Fig 3B). Therefore, the temporal increase in the proportion of NRU females using the MAB foraging area might be a true biological shift, not an artifact of using different methodologies in different years (satellite tracking from 1992-2008 vs. stable isotopes from 2004-2014). Because this shift coincides with a substantial increase in the total number of nests laid at NRU nesting sites since the early-to mid-2000s [51][52], the health and productivity of the MAB has clearly supported the ongoing recovery of the NRU loggerhead subpopulation. Patterns among isotopic data over this 10-year period assume that there has not been a shift in the isotopic baseline, a hypothesis that has not been tested. Nevertheless, even in the absence of an increasing trend, management actions and conservation efforts aimed to maintain the viability of foraging habitats in the MAB for NRU loggerheads and minimize threats to sea turtle survival in the MAB are therefore critical for this subpopulation's long-term stability and continued recovery [3,[53][54].
Post-nesting loggerheads that use the MAB foraging area migrate north within the continental shelf through a relatively narrow corridor along the North Carolina coast before occupying localized 'summertime' foraging sites between North Carolina and New Jersey, including productive habitats within the Chesapeake and Delaware bays [6,42]. In the winter months (November-March), these turtles migrate south along the same corridor to overwintering sites within the northern end of the SAB, where they can enter warmer waters adjacent to the Gulf Stream while minimizing the migratory distance, time and energy required to return to their northern foraging sites when water temperatures rise in the spring [6,42,55]. Similar seasonal MAB-SAB shuttling along this "migratory bottleneck" is exhibited by immature loggerheads [54].
Stable isotope data suggest that although these females occupy habitats in the SAB for part of the year, the majority of nutrients consumed and assimilated in the time period represented by the tissue reflect habitats in the MAB [15]. The frequency of this 'seasonal' MAB foraging strategy in NRU loggerheads suggests that the quality of foraging habitats in the MAB offsets the negative consequences associated with seasonal migrations [42,54], as well as possible bouts of winter fasting [55]. Indeed, the MAB is one of the most productive marine regions in the world [13], and turtles nesting in Georgia that forage in the MAB are larger and have larger clutch sizes and the same average remigration interval than turtles that remain in the SAB year-round [15,29]. Therefore, despite their longer northward migrations and seasonal movements, NRU females that forage in the MAB tend to have higher fecundity and likely contribute disproportionately more offspring to future generations.
NRU loggerheads that use the SAB foraging area move relatively short distances after nesting-both to the north and south of their nesting beaches-to foraging sites distributed within the continental shelf from North Carolina to central Florida. These turtles either remain in localized foraging sites year-round or undergo seasonal shuttling between summertime foraging sites along the inner shelf (<20m) and adjacent overwintering sites in the mid-or outer shelf (>20m) [42]. Similar to the seasonal north-south movements by turtles in the MAB, seasonal west-east movements by turtles in the SAB are associated with changes in water temperature and winter movements towards the warmer waters adjacent to the Gulf Stream for thermoregulation [42]. While benthic community characteristics and productivity in these areas are poorly understood, turtles that nest and forage in the SAB ('residents') do not appear to reinvest the energy needed for longer migrations into increases in growth and reproduction. Turtles using the SAB are in fact smaller and less fecund than turtles using the MAB [15]. Persistent differences in regional productivity might help explain why the proportion of females using the SAB is lower, and possibly decreasing, relative to the MAB.
The Subtropical Northwest Atlantic (SNWA) and Eastern Gulf of Mexico (EGoM) foraging areas support the majority of adult loggerheads in the Peninsular FL Recovery Unit (PFRU), making these areas globally important to the survival of the species (see below). However, very few NRU loggerheads use these foraging areas based on data from stable isotopes (Fig 4), satellite telemetry (Fig 5), as well as tag returns [56][57][58]. Those NRU turtles that do use the SNWA foraging area undergo lengthy post-nesting migrations southward along the continental shelf, often through a narrow corridor between south Florida and the Gulf Stream, before taking up year-round residence in localized home ranges in either the FL Keys or The Bahamas [6,42]. The latter of these females, along with a small number that make temporary oceanic forays adjacent to the MAB and SAB, represent rare instances when adult NRU loggerheads leave US waters.
Similar to NRU females that use the MAB foraging area, NRU females that use the SNWA are larger and more fecund than females that use the SAB [15]. Despite substantial differences in migratory direction and foraging strategy, NRU females that use MAB and SNWA foraging areas do not show differences in body size or clutch size. While NRU females from the MAB tend to have shorter remigration intervals than those from the SNWA [15], these differences may be due to small sample size (too few SNWA females) or to differences in site fidelity that artificially increase remigration intervals when nesting events/seasons are missed at a focal nesting site (e.g., Wassaw Island, GA). Additional work is needed to more thoroughly evaluate not only the reproductive carry-over effects of different foraging areas used by NRU loggerheads (e.g., clutch size, clutch frequency and the length of remigration intervals), but also how the use of different foraging areas affects nest-site fidelity, recruitment and adult survival.

Patterns across the NWA
The foraging-area origins of adult loggerheads shifts quite dramatically between NRU and PFRU nesting sites (Figs 2 and 4). While male PFRU turtles at a mating site in eastern FL predominately migrate to the MAB and SAB foraging areas after breeding [38], the proportion of PFRU female turtles that use the MAB foraging area decreases moving south along the east coast of FL as the predominant foraging area shifts from the SAB to the SNWA (Figs 4 and 5). Genetic subdivision between PFRU turtles nesting in central eastern vs. southeastern FL [3,59] is accompanied by this latitudinal shift in the predominant foraging area used by nesting females-MAB/SAB to SNWA. Nevertheless, data from both stable isotopes and satellite telemetry indicate that the MAB represents an important, although secondary, foraging area for PFRU loggerheads breeding along the east coast of Florida, a globally important nesting area for this species [60].
Like NRU females, these turtles use the same continental corridor adjacent to North Carolina during pre-and post-breeding migrations, as well as seasonal movements between summertime foraging grounds north of Cape Hatteras and overwintering sites south of Cape Hatteras [22,38]. Relative to NRU females, PFRU adults that use the MAB foraging area migrate greater distances to reach their foraging grounds. Conversely, PFRU females that use the SNWA foraging area migrate shorter distances than NRU females using the same foraging area (SNWA). While these differences in migratory distance are not associated with apparent differences in body size, clutch size, or clutch success (hatching or emergence), there is evidence for a tradeoff associated with remigration interval. Females that migrate to more distant foraging areas-PFRU females using the MAB and NRU females using the SNWA-exhibit longer remigration intervals (~5yr between nesting seasons) than females that migrate to more proximate foraging areas (~3yr between nesting seasons)-PFRU females using the SNWA and NRU females using the MAB [15][16]. Assuming that differences in nest-site fidelity do not drive apparent differences in remigration interval, this tradeoff between migratory distance and remigration interval might play an important role in regulating the frequency of nesting turtles that use different foraging areas.
The northern end of the SAB provides both an important migratory corridor and overwintering refuge for NRU and PFRU females that use productive foraging sites in the MAB. However, the southern end of the SAB, specifically the intercoastal waterways and adjacent continental shelf habitats around Cape Canaveral, FL, represents a foraging hotspot for PFRU turtles (both male and female) that use the adjacent breeding/nesting sites in east central Florida [9,22,37,40]. While a non-migratory or 'resident' foraging strategy might allow turtles to reallocate resources otherwise used in migration towards a greater investment in growth and reproduction, this behavior might also come at a cost if local foraging sites exhibit lower habitat quality and productivity or greater competition. Indeed, turtles (males and females) that remain resident in east central FL are smaller than migratory males [38,40] and 'resident' females show lower emergence success-a proxy for egg quality-than females that migrate south to the SNWA [9]. These differences, again, suggest a tradeoff between the cost of migration and the benefits of seeking out more productive foraging areas.
For nesting sites on the west coast of Florida, the prevalence of PFRU females that use the SNWA foraging area shifts to a greater contribution of females that use the EGoM foraging area, adjacent to their nesting beaches (Fig 5). Genetic subdivision between PFRU turtles nesting in southwestern vs. central western FL [59] is accompanied by subtle differences in the proportions of females originating from different foraging areas. Additionally, females from the distinct subpopulation that nest in the Dry Tortugas, FL predominantly use the SNWA foraging area, especially The Bahamas [12], while those in the Northern Gulf of Mexico Recovery Unit predominantly use the EGoM and Northern Gulf of Mexico (NGoM) foraging areas ( Fig  5). Because these patterns are based satellite tracking data alone, sample sizes are relatively low and population-wide inferences are potentially less robust. As exemplified by studies in the NRU and eastern PFRU, stable isotope analysis might provide a robust and cost-effective alternative for management agencies seeking to identify and monitor the foraging-area composition of loggerheads at nesting sites in the western PFRU, DTRU, and NGMRU.
Satellite-tracking studies suggest that like the aforementioned 'resident' PFRU turtles that forage and breed in east central FL (southern SAB), but unlike most NRU females, many turtles nesting in southern Florida and the Gulf of Mexico tend to migrate relatively short distances and use nesting sites located within their foraging area. Whether a more 'resident' strategy of using nesting sites within the same foraging area or the alternative strategy of using nesting sites in a different foraging area ('non-resident') entails a reproductive advantage (or disadvantage) for turtles nesting in southern FL and the Gulf of Mexico remains unknown. Reproductive carry-over effects (e.g., clutch size, clutch frequency, and remigration interval) between 'resident' and 'non-resident' females have yet to be analyzed. However, the quality and productivity of more proximate foraging habitats might be sufficiently high to support the energy demands of reproduction such that fewer females exhibit a 'non-resident' strategy.
Although never the most predominant foraging area used by turtles in this study, foraging sites along and adjacent to the Yucatan Peninsula in the Southern Gulf of Mexico (SGoM) support females from three of the four loggerhead Recovery Units in the United States (Fig 5). The SGoM may very well be the predominate foraging area used by the loggerheads nesting in Mexico, including beaches in Quintana Roo, Yucatan, an important rookery for the Greater Caribbean Recovery Unit [3]. However, neither stable isotope nor satellite tracking data have been used to identify the post-nesting foraging destinations of females in this rookery or other rookeries across the Greater Caribbean. Whereas tag returns indicate that females nesting in Mexico do use foraging areas outside the SGoM [61], the "major" foraging area designated for loggerheads in the SGoM [60][61][62][63] likely supports a mixed aggregation of adult turtles that predominately nest/breed in Mexico and northwest Caribbean, as well as turtles migrating from the southern United States and northern Caribbean (i.e., Cuba and The Bahamas, including Cay Sal).

Conservation implications
Patterns of foraging-area origins of nesting loggerheads in the NWA vary among recovery units: females in the NRU primarily use more northerly foraging areas in the MAB, while females from the three Florida recovery units (PFRU, DTRU, and NGMRU) primarily use more southerly foraging areas in the SNWA and EGoM (Figs 4 and 5). Recovery of NRU loggerheads appears to be associated with an increasing contribution of females foraging in the MAB, therefore managing anthropogenic threats to turtle survival and habitat quality in the MAB, as well as migratory corridors and overwintering sites in the SAB, remain conservation priorities for this subpopulation [3,[64][65]. Similarly, the recent rebound in PFRU nesting following a period of decline in the early 2000's [66][67] is likely linked to conditions in the SNWA and EGoM foraging areas. Managing threats and productivity in these foraging areas will continue to support the health and stability of the PFRU, a globally important subpopulation of loggerheads [3,[64][65].
Identifying the primary foraging area for each recovery unit is clearly important. However, the variation among foraging-area contributions displayed across nesting sites in Figs 4 and 5 confirm two important patterns that will continue inform future interpretations and actions. First, the patterns indicate that each recovery unit is composed of females that originate from multiple foraging areas, not just one. For this reason, any attempt to evaluate region-specific threats or management strategies by assessing changes in abundance and productivity at nesting sites would be confounded without monitoring changes in the foraging-area composition of the breeding population [9,50]. Second, the patterns indicate that each foraging area is used by females from multiple recovery units. Such 'mixed stock' aggregations characterize the foraging areas used by both adult and juvenile loggerheads [68][69][70]. Conservation actions aimed to mitigate prominent threats in one foraging area would therefore enhance the viability of multiple genetic stocks and life stages. For this reason, because PFRU loggerheads use all six foraging areas in the NWA, efforts to protect this globally important subpopulation would likely include management actions and policies that encompass the foraging areas also used by loggerheads from the less prominent recovery units: NRU, DTRU, and NGMRU. Protection of such 'umbrella' stocks in this and other marine turtle populations will help preserve other smaller, more vulnerable stocks, including stocks that receive less protection in other areas of their geographic range (e.g., loggerheads in the Greater Caribbean Recovery Unit).
Marine turtles, like other migratory species, are difficult to access during most of their lives. Breeding sites where migratory species congregate for reproduction (e.g., nesting beaches, spawning sites, rookeries, and calving grounds) offer invaluable opportunities to conduct population assessments of these threatened species [1]. As the patterns summarized in this study show, breeding aggregations are frequently composed of individuals that use different foraging areas, each associated with a different migratory distance and route, and distinguished by its own ecological characteristics, environmental conditions, anthropogenic threats, and management actions. Understanding the foraging-area composition of breeding aggregations is therefore critical for contextualizing changes in population abundance, productivity, and stability over time, as well as evaluating population-level responses to management strategies.  Table. Stable isotope data. Stable isotope ratios (following appropriate conversions) and posterior probabilities for foraging-area assignment of the 596 yolk samples and 150 skin samples collected from individual loggerhead turtles in this study. (CSV) S3 Table. Literature review. Summary of published articles that use either stable isotopes or satellite tracking to identify the foraging-area origins/destinations of nesting (females) or breeding (males) loggerhead turtles in the NWA. (CSV) S1 File. Metadata associated with S1 Table, S2 Table, and S3 Table. (PDF)