Individual Spatial Consistency and Dietary Flexibility in the Migratory Behavior of Northern Gannets Wintering in the Northeast Atlantic

to the rest of the population, suggesting there may be beneﬁts to ﬂexibility during the non-breeding period. Our results indicate that a consistent migratory strategy can be effective regardless of wintering region or diet, but that there may be beneﬁts to those individuals able to display ﬂexibility. This appears to be an important behavioral strategy that may enhance individual condition.

Migration is a fundamental behavioral process prevalent among a wide variety of animal taxa. As individuals are increasingly shown to present consistent responses to environmental cues for breeding or foraging, it may be expected that approaches to migration would present similar among-individual consistencies. Seabirds frequently show consistent individual differences in a range of traits related to foraging and space-use during both the breeding and non-breeding seasons, but the causes and consequences of this consistency are poorly understood. In this study, we combined analysis of geolocation and stable isotope data across multiple years to investigate individual variation in the non-breeding movements and diets of northern gannets Morus bassanus, and the consequences for changes in body condition. We found that individuals were highly repeatable in their non-breeding destination over consecutive years even though the population-level non-breeding distribution spanned >35 • of latitude. Isotopic signatures were also strongly repeatable, with individuals assigned to one of two dietary clusters defined by their distinct trophic (δ 15 N) and spatial (δ 13 C) position. The only non-breeding destination in which the two dietary clusters co-occurred was off the coast of northwest Africa. The majority of individuals adopted a consistent foraging strategy, as they remained within the same dietary cluster across years, with little variation in body mass corrected for size among these consistent individuals. In contrast, the few individuals that switched clusters between years were in better condition relative to the rest of the population, suggesting there may be benefits to flexibility during the non-breeding period. Our results indicate that a consistent migratory strategy can be effective regardless of wintering region or diet, but that there may be benefits to those individuals able to display flexibility. This appears to be an important behavioral strategy that may enhance individual condition.

INTRODUCTION
Animal migration is a fundamental behavioral process that involves seasonal movements between habitats in response to resource heterogeneity. Although prevalent among a wide variety of animal taxa, there is enormous variation in migration strategy in terms of the distance traveled and the degree of seasonal site fidelity, which ranges from strong philopatry to the loose tracking of seasonal resources (Webster et al., 2002;Newton, 2008). Variation among individuals is often attributed to age, sex or morphology (Marra, 2000;Alerstam et al., 2003;Bailleul et al., 2010) but may also be the result of differences in foraging behavior, breeding success or endogenous control (Bradshaw et al., 2004;Phillips et al., 2005;Broderick et al., 2007;Dias et al., 2011). Despite increasing evidence for individual differences in migratory behavior, the degree of consistency or plasticity and their causes and consequences remain incompletely understood (Chapman et al., 2011;Phillips et al., 2017).
A high degree of consistency in migration strategy with strong individual non-breeding site fidelity can be advantageous if this allows access to predictable foraging resources (Bradshaw et al., 2004;Weimerskirch, 2007) and reduces risks associated with exploring novel habitat (McNamara and Dall, 2010). Relatively inflexible strategies are seen in a number of groups, including passerines (Cuadrado et al., 1995), waterfowl (Hestbeck et al., 1991), cetaceans (Calambokidis et al., 2001), pinnipeds (Bradshaw et al., 2004), seabirds (Phillips et al., 2005, turtles (Broderick et al., 2007), and sharks (Jorgensen et al., 2010). Conversely, if food availability is unpredictable, or environmental conditions are prone to deteriorate in particular regions during the non-breeding period, selection should favor migratory flexibility or nomadism (Andersson, 1980) and facilitate plastic responses within individuals (Switzer, 1993;Sutherland, 1998). Such strategies are seen in groups including seabirds (Dias et al., 2011), waterbirds (Pedler et al., 2014), ungulates (Morrison and Bolger, 2012), and fish (Tibblin et al., 2016). Thus, the extent to which individuals respond to biotic and abiotic variation across time and space can select for clear individual differences in both movement and foraging strategies.
Recent studies of migrant birds show that individual differences in habitat selection and foraging behavior can influence diet quality during the non-breeding season and impact subsequent breeding traits such as body condition, timing of breeding, egg volume, and breeding success (Bearhop et al., 2004;Sorensen et al., 2009;Hoye et al., 2012), with important fitness consequences (Marra et al., 1998;Crossin et al., 2010;Inger et al., 2010;Harrison et al., 2011).
Thus, individuals that pursue a non-breeding strategy that produces strong negative carry-over effects might be expected to preferentially switch strategies in subsequent years, reducing within-individual consistency (Switzer, 1993;Dias et al., 2011;Morrison and Bolger, 2012). Understanding the incidence and implications of individual consistency or flexibility in nonbreeding behavior is therefore a key issue in animal ecology, yet there are few long-term studies that quantify these individual differences over multiple seasons or migration periods (Araújo et al., 2011;Phillips et al., 2017).
Marine predators such as seabirds provide an ideal model for examining such questions as they exhibit a broad spectrum of individual differences in behavior (Votier et al., 2010;Patrick et al., 2014). Recent work suggests these differences are likely to develop through ontogeny (Votier et al., 2017) as individuals learn to target profitable habitat (Grecian et al., 2018). In addition, many species display site fidelity to broadly productive regions during the non-breeding period (Grecian et al., 2016;Phillips et al., 2017). Disentangling individual differences in non-breeding foraging behavior and site fidelity may provide insights into how carry-over effects shape the annual cycle of an individual (Furness et al., 2006). For example, when local conditions are poor, individuals may switch non-breeding region while targeting the same preferred prey or, alternatively, may remain within the same preferred non-breeding region and instead switch prey types (Orben et al., 2015).
In this study, we combine multi-year deployments of geolocation loggers with stable carbon and nitrogen isotope analysis of winter-grown feathers to investigate the degree of individual consistency in the non-breeding destination and foraging behavior of a generalist marine predator, the northern gannet (Morus bassanus), tracked from four breeding colonies in the NE Atlantic. Gannets exhibit a southward-oriented chain migration following a flyway running along the coast of western Europe and Africa (Fort et al., 2012). Variation in migratory behavior, the migration path, final non-breeding destination and foraging behavior during these periods, occurs both among and within populations (Kubetzki et al., 2009;Fort et al., 2012;Deakin et al., 2019), and one recent study has shown that individuals in the NW Atlantic exhibit consistent behavioral strategies in successive years (Fifield et al., 2014). Additionally, the nonbreeding distributions of gannets may have changed in recent decades (Kubetzki et al., 2009), suggesting a degree of plasticity in migratory behavior. Such shifts could be linked to changes in human fishing activity as many seabirds are attracted to the foraging opportunities afforded by fisheries (Pichegru et al., 2007;Votier et al., 2010;Bodey et al., 2014a;Patrick et al., 2015). This behavior may come at a cost; as well as increasing the risk of bycatch (Bicknell et al., 2013), diets high in discards can have reduced lipid content compared to pelagic fishes, with the potential for adverse effects on body condition and breeding success (Grémillet et al., 2008;Votier et al., 2010). Should dependency on this resource also be evident in the nonbreeding season, there may be further fitness consequences via carry-over effects. We therefore examine whether differences in non-breeding destination and diet affect individual body condition (as a short-term fitness proxy) during the subsequent breeding season.

Study System and Data Collection
We collected data between 2008 and 2012 from gannets at four colonies in the northeast Atlantic: Bass Rock, Scotland; Grassholm, Wales; Great Saltee, Ireland; and Rouzic, France (Figure 1). In total, 187 breeding adults with chicks aged between 2 and 7 weeks (egg laying is poorly synchronized across the breeding colony) were caught at the nest during changeover of brood-guard duties using a brass noose or crook attached to the end of a carbon fiber pole ( Table 1). On capture, the mass (to the nearest 50 g) and bill length (to the nearest 0.1 mm) of each individual was measured, and sex was subsequently assigned from DNA using 2550F, 2718R, or 2757R primers (Griffiths et al., 1998;Fridolfsson and Ellegren, 1999) following Stauss et al. (2012).

Non-breeding Destination
Combined geolocation-immersion loggers (Mk 19, 15, and 5, British Antarctic Survey, Cambridge UK) were deployed on 77 of these individuals across the four colonies. Loggers were attached with two cable ties to a plastic ring, which was then fitted to the tarsus and remained in place for up to 2 years before the bird was recaptured at the breeding colony. The total mass of the attachment did not exceed 10 g, representing <0.35% of average adult body mass, and so unlikely to have any adverse effects (Bodey et al., 2018a). The loggers sampled ambient light every minute and recorded the maximum value every 2, 5, or 10 min (Mk 19, 15, and 5 loggers, respectively). Positions were calculated from logger data following established methods (Wilson et al., 1992;Phillips et al., 2004). Briefly, the timings of sunset and sunrise were estimated using TransEdit2 (British Antarctic Survey, Cambridge UK) using a light-intensity threshold of 16. A minimum dark period of 4 h was set to remove any light-dark transitions created by shading or cloud cover. Latitude was derived from day length, and longitude from the timing of local midday and midnight, with respect to Greenwich Mean Time and Julian day, providing two positions per day with an accuracy of ∼200 km (Phillips et al., 2004). Examination of individual migration tracks revealed latitude to be the major axis of movement, with birds tending to migrate southward from the breeding colonies toward northwest Africa (Fort et al., 2012; Figure 2). Plots of displacement from the colony indicated that all individuals reached their final non-breeding destinations by December and remained in this region for a minimum of 1 month before commencing their return migration. The mean latitude and longitude for December was therefore used as the non-breeding destination of each bird.

Non-breeding Stable Isotopes
Small samples from the 8th primary feather were taken from 148 individuals for stable isotope analysis, with 43 of these individuals sampled a second time when loggers were removed the following year ( Table 1). Gannets perform a complete annual molt after the breeding season (from September; Ginn and Melville, 1983), suspending molt by the time the return to the breeding colony (January to March) to invest in nest attendance and foraging trips (see Nelson, 2006). Thus, as feathers are metabolically inert after formation and larger feathers grow over a protracted period, the stable isotope ratios of primary feathers were assumed to largely represent prey consumed at the non-breeding grounds (between October and December).
Feather samples were thoroughly washed with distilled water and placed in a drying oven at ∼40 • C until dry. The barbules were cut into fine pieces and subsamples of 0.7 ± 0.1 mg were weighed into tin cups. Stable isotope analysis of these subsamples was then conducted at the East Kilbride Node of the Natural Environment Research Council Life Sciences Mass Spectrometry Facility via continuous flow isotope ratio mass spectrometry, using a Thermo Fisher Scientific Delta V Plus with a Costech ECS 4010 elemental analyser configured for simultaneous 13 C/ 12 C and 15 N/ 14 N isotope analysis. Stable isotope ratios are reported in δ notation, expressed as parts per thousand (‰) deviation according to the equation δX = [(R sample /R standard )-1], where X is 13 C or 15 N, R is the corresponding ratio 13 C/ 12 C or 15 N/ 14 N, and R standard is the ratio of the international references VPDB for carbon and AIR for nitrogen. At set intervals, standards of GEL, 14 N ALA, glycine and tryptophan were analyzed between feather samples in the IRMS. The measurement precision, calculated as the standard deviation of multiple analyses of these standards, was ± 0.1 ‰ for δ 13 C and ± 0.2 ‰ for δ 15 N.

Consistency in Non-breeding Strategies and Isotopic Clustering
To examine the consistency of non-breeding destination and stable isotope ratios we calculated the repeatability of these traits based on the intra-class correlation coefficient from linear mixed-effect models fitted with bird ID as a random intercept, using the package "rptR" v. 0.9.21 (Stoffel et al., 2017). We used repeatability as a proxy for behavioral consistency, testing the hypothesis that between-individual variance in a particular trait was greater than within-individual variance (Patrick et al., 2014). To estimate the consistency of non-breeding destinations we calculated the repeatability of mean December latitude and longitude for those individuals from Bass Rock (n = 22) and Rouzic (n = 21), that were tracked over two consecutive years.
Frontiers in Ecology and Evolution | www.frontiersin.org To estimate the consistency in stable isotope ratios during the non-breeding season we calculated the repeatability of δ 13 C and δ 15 N in primary feathers of individuals from Bass Rock (n = 27) and Grassholm (n = 13) sampled in two consecutive years. The three individuals from Great Saltee were excluded from the estimate of isotopic repeatability due to the small, multi-year sample size (Table 1).
To test for the occurrence of distinct dietary clusters in stable isotope ratios we fitted a multivariate normal mixture model to δ 13 C and δ 15 N values using the package "mixtools" v. 1.1.0 (Benaglia et al., 2009). The best-fitting model was selected by comparing the log-likelihood of candidate models with differing numbers of clusters. Feather samples were assigned to a dietary cluster based on a probability of assignment >0.5.

Consequences of Non-breeding Strategy
We estimated body condition using a scaled mass index (SMI, Peig and Green, 2009) with bill length as a linear measurement of body size in relation to body mass. However, given that 53 of the individuals were measured in multiple years, we extended this approach to a mixed-effects model with an individual level random intercept, fitted using the package "lme4" v. 1.1-18-1 (Bates et al., 2015). The correlation between body mass and bill length accounting for repeated measures was estimated using the package "rmcorr" v. 0.3.0 (Bakdash and Marusich, 2017). SMI allows for the comparison of the relative size of energy reserves of individuals within a population, avoiding the assumption that larger animals have better body condition due to a higher absolute mass (Peig and Green, 2009). While reproductive success would be a more robust measure of fitness, chick survival from hatching to fledging is over 90% and the majority of offspring mortality occurs during the post-fledging and juvenile period at sea (Nelson, 1966).
The implications of alternative non-breeding strategies (destination and dietary cluster) at the individual level were explored by examining the effects of sex, breeding colony, dietary cluster and non-breeding destination on scaled mass using linear regressions. In cases where there were two observations of an individual's scaled mass in consecutive years, these were fitted as mixed-effects models with individual as a random intercept term. Model selection of linear regressions was based on the F statistic, model selection of mixed-effects linear regression was based on the Chi-squared statistic using likelihood ratio tests. Post-hoc comparisons were made using the package "lsmeans" v. 2.30-0 (Lenth, 2016). All analyses were carried out in R v. 3.4.3 (R Core Team, 2018).

Consistency in Non-breeding Destination and Stable Isotope Ratios
Gannets spent the month of December in one of three regions: a northern region (>36 • N), around the British Isles and the Bay of Biscay (n = 26); a southern region (<36 • N) from Gibraltar to Mauritania (n = 47); and the Mediterranean Sea (n = 4, all from Rouzic; Figure 2). Individuals from all four colonies were present in both the northern and southern regions during the non-breeding period. The 43 birds from Bass Rock and Rouzic that were tracked over two consecutive years exhibited a high degree of consistency in non-breeding destination and were highly repeatable in both mean December latitude (

Isotopic Clustering
Stable isotope ratios in primary feathers sampled from 148 individuals were best described by a mixture of k = 2 multivariate normal distributions ( Table 2). One cluster centered on −13.9 δ 13 C and 13.2 δ 15 N, and a second cluster centered on −16.1 δ 13 C and 15.8 δ 15 N. The 95% ellipses of the two multivariate normal distributions did not overlap (Figure 3). Seventy-three individuals were assigned to cluster 1 and 75 individuals to cluster 2. Of the 43 individuals (Bass Rock n = 27; Grassholm n = 13; Great Saltee n = 3) that were sampled in consecutive years, most were consistent in their cluster assignment, with 16 assigned to cluster 1 and 20 assigned to cluster 2 in both years. Nevertheless, seven individuals switched between clusters from 1 year to the next (Bass Rock n = 3; Grassholm n = 2; Great Saltee n = 2). Six of these were female and switched from cluster 1 to cluster 2 (lower δ 15 N to higher δ 15 N) and one male from Grassholm switched from cluster 2 to cluster 1 (higher δ 15 N to lower δ 15 N).

Isotopic Clustering Controlling for Winter Destination
Both non-breeding destination and stable isotope data were available for 56 individuals (Bass Rock n = 35, Grassholm n = 13, Great Saltee n = 8). Colony of origin was unrelated to cluster assignment (χ 2 = 0.51, P = 0.78) or non-breeding region (χ 2 = 0.13, P = 0.94). However, individuals that wintered in the northern region were all assigned to cluster 2 (higher δ 15 N) whereas individuals that wintered in the southern region were assigned to either isotopic cluster (Figure 4). No isotope data were available for the four individuals that wintered in the Mediterranean.
Frontiers in Ecology and Evolution | www.frontiersin.org FIGURE 3 | Primary feather δ 13 C and δ 15 N values of 148 individual gannets sampled at four colonies in the northeast Atlantic. Crosses represent the mean of the two isotopic clusters identified from a multivariate mixture model and dotted lines represent the 95% ellipses for each distribution. Points are colored based on a probability of group assignment >0.5. Gray lines connect individuals sampled over two consecutive non-breeding periods (n = 43).
There was no difference in scaled mass between individuals using either the northern or southern non-breeding region (χ 2 1 = 0.58, P = 0.44), nor were there differences in scaled mass between birds in the two isotopic clusters (χ 2 1 = 0.0, P = 0.96). Data on scaled mass and feather stable isotope ratios in two consecutive years were available for 34 individuals sampled at Bass Rock (n = 21) and Grassholm (n = 13). In this sample, five individuals switched between the two isotopic clusters and had higher scaled mass (∼200 g heavier) compared to those that did not switch isotopic cluster, after accounting for both colony and sex differences (Figure 5). Post-hoc comparisons of marginal means indicated individuals that switched were significantly heavier than individuals in the high δ 15 N cluster (z = 2.53, P = 0.03). Of these five switching individuals, four were female and switched from the low to high δ 15 N cluster and one was male and switched from the high to

DISCUSSION
In this study, we reveal that the non-breeding behavior of individual northern gannets is highly repeatable over consecutive years, with a high degree of site fidelity and consistency in stable isotope ratios during successive non-breeding seasons. Despite substantial differences in destination and variation in foraging strategy among individuals, consistent behaviors during the nonbreeding period had no apparent carry-over effect on scaled mass in the subsequent breeding season.

Consistency in Non-breeding Destination and Stable Isotope Ratios
Gannets tracked in this study tended to migrate uniformly southward on a known flyway (Kubetzki et al., 2009;Fort et al., 2012), and spent the non-breeding period in a wide variety of marine habitats including the North Sea, Bay of Biscay, Mediterranean Sea, and Canary Current Upwelling region (Grecian et al., 2016 ; Figure 2). These three regions differ in their environmental conditions, yet individuals tracked over two consecutive years displayed a high degree of nonbreeding site fidelity. The range of ∼4 ‰ in δ 13 C and ∼6 ‰ in δ 15 N in stable isotope data from the broader sample of the population are larger than the estimates of baseline isotopic variation across the differing non-breeding Frontiers in Ecology and Evolution | www.frontiersin.org 6 June 2019 | Volume 7 | Article 214 FIGURE 5 | Partial effect of dietary strategy on scaled mass (g) for northern gannets sampled in two consecutive years at Bass Rock and Grassholm (n = 34) after accounting for effect of Colony and Sex. Individuals were assigned to isotopic clusters in each year based on a multivariate normal mixture model of primary feather δ 13 C and δ 15 N. Individuals (n = 5) that switched between cluster one and cluster two were heavier relative to those individuals that remained within either cluster one (n = 15) or cluster two (n = 14). Post-hoc comparisons (*) P < 0.05.
destinations (McMahon et al., 2013;Magozzi et al., 2017). This suggests that while the population winters across a range of locations, prey are targeted across trophic levels within locations (Inger and Bearhop, 2008). Adult gannets display consistency in foraging movements and diet within breeding seasons (Patrick et al., 2014;Wakefield et al., 2015;Votier et al., 2017;Bodey et al., 2018b), and the high isotopic consistency observed in individuals in our study that were sampled in consecutive years suggests a similar degree of consistency in both wintering region and the trophic level of prey consumed. Non-breeding site fidelity has been documented in other migratory marine vertebrates (Bradshaw et al., 2004;Broderick et al., 2007;Jorgensen et al., 2010;Phillips et al., 2017) and could allow individuals to increase knowledge of a specific area and thus improve foraging efficiency (Dall et al., 2012).

Isotopic Clustering
Pooling the stable isotope data from all colonies indicated two clusters, indicative of alternative foraging strategies that differed in both spatial (δ 13 C, δ 15 N) and trophic (δ 15 N) characteristics. One cluster was described by higher δ 15 N and depleted δ 13 C, consistent with a higher trophic level diet and offshore prey, respectively (Hobson et al., 1994;Post, 2002;Inger and Bearhop, 2008). In contrast, the second cluster was more representative of a diet of inshore (higher δ 13 C) prey at a lower trophic level (depleted δ 15 N). Although δ 15 N can also vary with geographic location (Seminoff et al., 2012;McMahon et al., 2013), the observed difference between these two clusters is greater than the baseline variation between non-breeding destinations (McMahon et al., 2013;Magozzi et al., 2017). In addition, the co-occurrence of individuals from both dietary clusters in the southern wintering area indicates that cluster assignment is not purely driven by the local environment. However, there may be other drivers of the observed isotopic differences, for example individual variation in molt location or feather growth rate could result in a shift in feather isotope signature.
While we lack conventional samples of gannet diet during the winter, the higher trophic level cluster may represent prey obtained primarily as fisheries discards, as δ 15 N values are elevated in demersal relative to pelagic fish (Votier et al., 2010;Bicknell et al., 2013). In contrast, the second cluster is suggestive of a more inshore diet in pursuit of small forage fish (Garthe et al., 2000;Nelson, 2002). The majority of individuals that were sampled in two consecutive years remained in the same isotopic cluster from one year to the next. Therefore, these clusters may represent foraging strategies that reduce competition among individuals though niche differentiation (Phillips et al., 2009;Young et al., 2010;Bodey et al., 2014b). Foraging specializations have been documented during the breeding season for many seabird species (Annett and Pierotti, 1999;Bearhop et al., 2006;Woo et al., 2008;Phillips et al., 2017) including northern gannets where individuals can vary in the extent of their reliance on high trophic level prey such as fisheries discards (Votier et al., 2010;Patrick et al., 2014;Wakefield et al., 2015;Bodey et al., 2018b).

Consequences of Non-breeding Strategy
Based on our metric of scaled mass, we did not detect any consequences for individuals consistently pursuing different non-breeding strategies. Neither non-breeding destination nor isotopic cluster was significantly related to scaled mass; instead, sex and colony effects drove the observed differences. This is in contrast to patterns seen in thick-billed murres Uria lomvia, where over-wintering foraging strategies are strongly dependant on body size (Orben et al., 2015). Differences in energetic demands over the breeding season may also lead to variation in body condition (Moe et al., 2002) but all individuals in this study were sampled during the chick provisioning period. Sex-linked differences in scaled mass have been observed previously in Northern gannets and may reflect the differing physiological demands of reproduction, and breeding role specialization, as well as more subtle differences between the sexes in prey-capture techniques, nutritional requirements and fine-scale habitat and prey selection (Stauss et al., 2012;Cleasby et al., 2015;Machovsky-Capuska et al., 2016;Bodey et al., 2018b). The difference in scaled mass between individuals at Grassholm and Bass Rock may reflect variation in the prey resources and environmental conditions Frontiers in Ecology and Evolution | www.frontiersin.org 7 June 2019 | Volume 7 | Article 214 accessible to individuals from their respective colonies. For example, the North Sea differs in oceanography to the Celtic Sea and supports fewer competing gannet colonies, though the colony at Bass Rock is much larger than at Grassholm and so within-colony competition will be more severe (Nelson, 2002;Wakefield et al., 2013). Some individuals remained close to the breeding colonies during the non-breeding period and were consistent in this behavior over the two years (Figure 2). Remaining in these areas may decrease energy expenditure by reducing migration costs (Flack et al., 2016), however, this may be offset by the increased energetic requirement for thermoregulation in these colder more northerly waters . These individuals were all assigned to the higher δ 15 N cluster which suggests a greater consumption of fisheries discards or a lack lower δ 15 N prey available during the winter period (e.g., shoaling forage fish).
The small number of individuals that switched between the two dietary clusters were in better body condition, after accounting for colony and sex effects, than those that were consistent in their cluster assignment. This difference was largest compared to individuals in the higher δ 15 N cluster, which had relatively low scaled mass. The switching strategy was observed in seven of the 43 individuals (ca. 16%) that were sampled in two consecutive years. Six of these were female and all switched from the lower to higher δ 15 N cluster. The only individual to switch from the higher to lower δ 15 N cluster was male. Although foraging on discards brings an additional risk of mortality via incidental bycatch (Bicknell et al., 2013), previous work suggests such a diet may not be detrimental to adult body condition in Cape gannets (Grémillet et al., 2008). Our findings suggest that individuals capable of switching between higher and lower trophic level diets between non-breeding seasons may benefit when compared to individuals specializing in a diet likely to consist of a high proportion of fisheries discards. Alternatively, individual in better condition may be the only ones capable of investing in more risky behaviors (Geary et al., 2019). The majority of individuals switched to the higher δ 15 N cluster, so this may indicate a short-term benefit of switching to a diet based on fisheries discards or other alternative higher δ 15 N resources within non-breeding region. The higher δ 15 N cluster represents one third of those individuals wintering off the coast of northwest Africa, a region known to have experienced a recent intensification of fishing activity (Worm et al., 2009;Grecian et al., 2016).

CONCLUSIONS
Our results reveal strong individual consistencies in movement and diet during the non-breeding season, and it is this consistency rather than the strategy itself, that may be important for long-lived species (Ceia et al., 2012;Gilmour et al., 2015). Indeed, such consistency has been demonstrated to result in similar life-time reproductive success among Brünnich's guillemots (Uria lomvia) pursuing different foraging strategies (Woo et al., 2008). Individual repeatability is frequently seen in marine vertebrates despite strong between-year variation in environmental variables and prey fields (Cherel et al., 2007). Importantly however, such consistency could come at a price for highly specialized individuals; for example, changes to anthropogenic subsidies disproportionately affect sub-sections of populations that specialize on such resources (Whitehead and Reeves, 2005;Bicknell et al., 2013). The findings here further highlight the importance of research that links different aspects of behavior between seasons or across annual cycles to understand ecological differentiation among individuals, populations and species (Friesen et al., 2007;Bodey et al., 2014b;Wakefield et al., 2015); and the need to consider the degree of flexibility of individuals and populations to changes in resource availability (Grémillet and Boulinier, 2009

ETHICS STATEMENT
Birds were ringed and loggers deployed with permits and approval from the British Trust for Ornithology and Scottish Natural Heritage. Tissue samples were collected under license from the UK Home Office (PPL 30/3065 and 40/3408).

AUTHOR CONTRIBUTIONS
SV, SB, and KH conceived the study. WG, HW, and TB wrote the first draft of the manuscript. WG, HW, SV, SB, IC, DG, KH, ML, AL, SP, EW, and TB collected data. WG, HW, IC, JN, RP, EW, and TB conducted analyses. All authors contributed to manuscript revision, read and approved the submitted version.

FUNDING
This study was funded by the UK Natural Environment Research Council (NERC Standard Grant NE/H007466/1). SB is currently funded by an ERC consolidator's grant (310820: STATEMIG). EW is funded by the NERC (Independent Research Fellowship NE/M017990/1). Tissue samples were collected under license from the UK Home Office (PPL 40/3408).