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

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, 2006), 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; Inger et al., 2008; 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 non-breeding 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 non-breeding 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 non-breeding 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.

Materials and Methods

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).

Figure 1. Location of the four northern gannet study colonies across the UK, Ireland and France.

Table 1. Summary of samples collected from 187 northern gannets across four breeding colonies between 2008 and 2012 including geolocation loggers and feather stable isotope analysis (SIA).

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.

Figure 2. Centroid December location of northern gannets (n = 77) tracked from four colonies in the North East Atlantic. Dark lines connect locations of individuals (n = 43) from Bass Rock and Rouzic tracked over two consecutive non-breeding periods.

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 ¹³C/¹²C and ¹⁵N/¹⁴N isotope analysis. Stable isotope ratios are reported in δ notation, expressed as parts per thousand (%0) deviation according to the equation δX = [(R_sample/R_standard)-1], where X is ¹³C or ¹⁵N, R is the corresponding ratio ¹³C/¹²C or ¹⁵N/¹⁴N, and R_standard is the ratio of the international references VPDB for carbon and AIR for nitrogen. At set intervals, standards of GEL, ¹⁴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 δ¹³C and ± 0.2 ‰ for δ¹⁵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. To estimate the consistency in stable isotope ratios during the non-breeding season we calculated the repeatability of δ¹³C and δ¹⁵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 δ¹³C and δ¹⁵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).

Results

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 (R = 0.91; CI = 0.83, 0.95; P = 0.001) and longitude (R = 0.92; CI = 0.87, 0.96; P = 0.001; Figure 2).

Stable isotope values in primary feathers from individuals sampled in two consecutive years at Grassholm (n = 13) were repeatable with respect to both δ¹³C (R = 0.73; CI = 0.35, 0.90; P = 0.004) and δ¹⁵N (R = 0.57; CI = 0.06, 0.82; P = 0.028). Individuals from Bass Rock (n = 27) also showed significant repeatability in both δ¹³C (R = 0.59; CI = 0.27, 0.79; P = 0.002), and δ¹⁵N (R = 0.52; CI = 0.20, 0.74; P = 0.002).

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 δ¹³C and 13.2 δ¹⁵N, and a second cluster centered on −16.1 δ¹³C and 15.8 δ¹⁵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 δ¹⁵N to higher δ¹⁵N) and one male from Grassholm switched from cluster 2 to cluster 1 (higher δ¹⁵N to lower δ¹⁵N).

Table 2. Comparison of the log-likelihoods of multivariate normal mixture models fitted with k distributions.

Figure 3. Primary feather δ¹³C and δ¹⁵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).

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 (χ² = 0.51, P = 0.78) or non-breeding region (χ² = 0.13, P = 0.94). However, individuals that wintered in the northern region were all assigned to cluster 2 (higher δ¹⁵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.

Figure 4. Centroid December locations of 56 northern gannets colored by assignment to one of two isotopic clusters based on a multivariate mixture model of primary feather δ¹³C and δ¹⁵N. Dark lines connect individuals (n = 22) sampled in two consecutive years.

Consequences of Non-breeding Strategy

Female gannets were estimated to be on average 89.4 g (95% CI; 43.9, 134.9) heavier than males (${\chi }_1^2$ = 12.9, P < 0.001). Scaled mass differed slightly between the four colonies (${\chi }_3^2$ = 9.7, P = 0.02) and post-hoc comparisons indicated that, when averaging over sex differences, individuals sampled at Grassholm were 179.6 g ± 58.0 g lighter than individuals sampled at Bass Rock (z = 3.1, P = 0.01) with no other significant differences between colonies.

There was no difference in scaled mass between individuals using either the northern or southern non-breeding region (${\chi }_1^2$ = 0.58, P = 0.44), nor were there differences in scaled mass between birds in the two isotopic clusters (${\chi }_1^2$ = 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 δ¹⁵N cluster (z = 2.53, P = 0.03). Of these five switching individuals, four were female and switched from the low to high δ¹⁵N cluster and one was male and switched from the high to low δ¹⁵N cluster. Scaled mass was unavailable for two other switching individuals.

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 δ¹³C and δ¹⁵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.

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 non-breeding 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 non-breeding site fidelity. The range of ~4 %0 in δ¹³C and ~6 %0 in δ¹⁵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 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 (δ¹³C, δ¹⁵N) and trophic (δ¹⁵N) characteristics. One cluster was described by higher δ¹⁵N and depleted δ¹³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 δ¹³C) prey at a lower trophic level (depleted δ¹⁵N). Although δ¹⁵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 δ¹⁵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 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 (Garthe et al., 2012). These individuals were all assigned to the higher δ¹⁵N cluster which suggests a greater consumption of fisheries discards or a lack lower δ¹⁵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 δ¹⁵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 δ¹⁵N cluster. The only individual to switch from the higher to lower δ¹⁵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 δ¹⁵N cluster, so this may indicate a short-term benefit of switching to a diet based on fisheries discards or other alternative higher δ¹⁵N resources within non-breeding region. The higher δ¹⁵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).

Data Availability

Telemetry data are available through the BirdLife International Seabird Tracking Database: http://www.seabirdtracking.org. Biometric data are available through the University of St Andrews Research Portal: https://doi.org/10.17630/b3c6dc92-13eb-447d-82d3-acf01d029bc9 (Grecian et al., 2019).

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).

Conflict of Interest Statement

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Acknowledgments

We thank Sir Hew Hamilton-Dalrymple (Bass Rock), the RSPB (Grassholm), the Neale family (Great Saltee) and the Ligue pour la Protection des Oiseaux and Réserve Naturelle Nationale des Sept- Îles (Rouzic) for permitting and facilitating fieldwork. We are grateful to Fraser Bell, Nadja Christen, Armel Deniau, Richard Inger, Mark Jessopp, Greg & Lisa Morgan, Claudia Stauss, Alyn Walsh, Emma Wood, the Scottish Seabird Center, Maggie Sheddan, and Venture Jet for field assistance.

References