Ocean-scale variation in migration schedules of a long-distance migratory seabird is fully compensated upon return to the breeding site

Background Migratory birds generally have tightly scheduled annual cycles, in which delays can have carry-over effects on the timing of later events, ultimately impacting reproductive output. Whether temporal carry-over effects are more pronounced among migrations over larger distances, with tighter schedules, is a largely unexplored question. Methods We tracked individual Arctic Skuas Stercorarius parasiticus, a long-distance migratory seabird, from eight breeding populations between Greenland and Siberia using light-level geolocators. We tested whether migration schedules among breeding populations differ as a function of their use of seven widely divergent wintering areas across the Atlantic Ocean, Mediterranean Sea and Indian Ocean. Results Breeding at higher latitudes led not only to later reproduction and migration, but also faster spring migration and shorter time between return to the breeding area and clutch initiation. Wintering area was consistent within individuals among years; and more distant areas were associated with more time spent on migration and less time in the wintering areas. Skuas adjusted the period spent in the wintering area, regardless of migration distance, which buffered the variation in timing of autumn migration. Choice of wintering area had only minor effects on timing of return at the breeding area and timing of breeding and these effects were not consistent between breeding populations. Conclusion The lack of a consistent effect of wintering area on timing of return between breeding areas indicates that individuals synchronize their arrival with others in their population despite extensive individual differences in migration strategies.


Introduction
Stages in the annual cycle of animals are intricately linked and the timing of each event can have cascad ing effects on the timing of subsequent events, usually referred to as carryover effects [1,2].Temporal carryover effects may be particularly pronounced in animals that undertake long and challenging migra tions between their breeding and wintering areas, where delays in one phase of the annual cycle can lead to a temporal mismatch with resource availability in subsequent phases, ultimately affecting survival and reproductive success [3,4].The degree to which delays can be buffered, preventing amplification of carry over effects across time, may vary as a function of the tightness of scheduling within the annual cycle [5].
Quantifying the geographical variation in timing of annual cycles and carryover effects is therefore important for our understanding of how phenology drives variation in lifehistory traits [6].In addition, considering that daily survival rates may differ between stages of the annual cycle [7,8], variation in annual schedules can drive metapopulation dynamics [2,9,10].Insight in the drivers of metapopulation dynamics is important to inform locationspecific conservation efforts.
Populationlevel differences in the scheduling of the annual cycle may arise as a function of the location of the breeding and wintering areas [11,12] and the conditions experienced there [13].Annual schedules are generally shifted to later in the season at higher breeding latitudes, where the optimal time for reproduction is associated with a later time of return to the breeding site [18].In turn, this affects the duration and timing of subsequent phases of the annual cycle, such as migration, given that the degree of seasonality in resources and photoperiod leads to differences in optimal migration speed [19,20]; however, this is little studied.In addition to an effect of the breeding area, the timing of annual cycles may also depend on the choice of wintering area and migration routes, as more distant nonbreeding areas require longer travel time.As such, unless departure is earlier or travel speeds higher, arrival will be later [21,22].When fueling conditions in the wintering area itself or along the migration route are suboptimal, delays can arise as a result of choice of wintering area.These delays can subsequently carryover to timing of return to the breeding area and egg laying [23,24] and can reduce survival [25].Thus, temporal carryover effects are expected to be more pronounced when distances between breeding and nonbreeding areas are greater, and there is less leeway in the annual cycle to compensate for delays.
We define the net carryover effect as the timing of an event for an individual or group relative to the pop ulation mean for a sequence of events (figure 1a).The strength of temporal carryover effects can be conceptualized as the degree to which timing of an event affects the scheduling of a following event (figure 1ac).If there is full compensation for a delay, each day lost in one phase translates to one less day spent in the next phase, resulting in a similar timing of the following event.In the case of complete carryover of delays, the timing of the next event is delayed by a day, i.e. there is no compensation.Delays may also be amplified at later stages, i.e. each day of delay translates to more than a day delay in following events.
Although longdistance migration has evolved in many taxonomic groups [26], geographical variation in the timing of annual cycles and carryover effects have mainly been studied in single breeding populations of migratory birds [13,27,28], that use a single wintering area or multiple wintering areas at similar latitudes [11,12,15] or in multiple breeding populations that breed at similar latitudes [29].Ideally, quantifying and testing the relative effects of breeding and wintering location on the timing of annual schedules and on carryover effects requires substantial geographic variation in both breeding and wintering areas and mixing of individuals from different breeding populations within wintering areas and vice versa (weak migratory connectivity [30]).Such studies are rare, because of the challenges of performing largescale tracking studies and because few species show sufficient variation in both breeding and wintering latitudes.Hence, the extent to which the timing of events and carryover effects in the annual cycle of migratory birds is driven by both breeding and wintering locations is largely unresolved [31].
We studied carryover effects on the timing of different phases in the annual cycle of a longdistance migra tory seabird, the Arctic Skua Stercorarius parasiticus, from multiple breeding populations.Its circumpolar breeding distribution ranges from northern temperate to high arctic zones [[32]; figure 2].Based on ring recoveries and field sightings, the main wintering areas in the Atlantic are thought to span productive conti nental shelf and pelagic areas across a large (latitudinal) range, from Iberia to Patagonia and South Africa [32,33].By tracking the full annual cycle of many individuals over several years from multiple breeding sites ranging from East Greenland to West Siberia (59°N to 79°N in latitude and 24°W to 69°E in longitude), we present the full extent of their wintering distribution.We quantify the migratory connectivity [30] and the consistency of individuals from year to year in their migration strategy.We quantify the relative effects of breeding and wintering areas on the timing of key events and duration of intervening periods in the annual cycle.Based on the northern hemisphere (breeding area) viewpoint on annual schedules, the nonbreeding area is referred to as wintering area throughout.Using the variation in annual schedules, we study com pensation versus carryover of variation in timing between subsequent events, and how these depend on breeding and wintering areas.Later timing of events and tighter schedules are predicted to be associated with more northerly breeding areas [15,16]; later timing is also expected to lead to different allocation of time to phases of the annual cycle.Time spent in more distant wintering areas is expected to be shorter due to longer time spent on migration.Among different annual schedules, we expect carryover effects to arise particularly among schedules with less leeway to compensate for delays and that the wintering period gen erally buffers the variation in timing of autumn migration, but less so when the wintering period is shorter.As such, migrations to more distant wintering areas are expected to be associated with higher risks of delayed spring return in the breeding area a net carryover effect [24,34].
Birds were captured either on the nest using bow nets, tent spring traps or walkin traps, remote release nooses, or away from the nest using net guns.Geolocators were attached to a Darvic ring fitted around the tarsus.Various models of geolocators were used: C65, C250 (Migrate Technology, Cambridge, UK), Mk9, Mk13, Mk15 and Mk18H (British Antarctic Survey, Cambridge, UK) and Mk3006 (Biotrack, Wareham, UK).
Geolocators measured ambient light in lux (Migrate Technology) or arbitrary units (BAS/Biotrack) every 1 minute and saved the maximum value every 5 minutes.In subsequent breeding seasons, geolocators were removed and data downloaded.
In total, data for 276 interbreeding periods (i.e. one migration cycle, referred to as a 'track' below) from 155 individuals were obtained in 20092021.Sample sizes per breeding area ranged from 103 tracks of 57 individuals in mainland Norway to 8 tracks of 3 individuals in East Greenland (figure 2).In total, 34, 23 and 12 individuals were tracked for 2, 3 and 46 years, respectively.Some of the tracks (n=48) were incomplete due to battery or electronics failure.

Geolocator data processing
Positions were estimated for noon and midnight from the light measurements by the geolocators (see elec tronic supplementary material for details).
We describe the phenology from the perspective of a breeding bird, i.e. 'autumn migration' is postbreeding, and 'spring migration' is prebreeding.We define the following key events in the annual cycle: 1) departure from the breeding area; 2) arrival in the wintering area; 3) departure from the wintering area; 4) arrival in the North Atlantic north of 35°N (mostly around the subpolar frontal zone where many Arctic Skuas made a stopover before returning to the breeding grounds [35,36]); 5) return to the breeding area; 6) clutch initiation.
The duration of the intervening periods (in days) were calculated for the: 1) autumn migration; 2) wintering period; 3) spring migration, including the time in the North Atlantic; 4) time in the North Atlantic; 5) pre laying period; 6) breeding period.As loggers were deployed in the colony and retrieved before departure, a composite index of the length of the breeding period for each track was calculated as the sum of the periods from 1 July (all individuals present in the breeding area) to departure in the first year, and from clutch initiation in the second year to 30 June.Timing of movements were identified by inspecting raw position estimates.
For birds breeding in areas with continuous daylight, first position estimates after leaving the zone with continuous daylight were close to the breeding site, but last position estimates, just before returning to the zone with continuous daylight, were around 60°N and distant from the breeding area.To estimate time of departure from or arrival to the breeding area under continuous daylight, we assumed birds traveled at a speed of ca.850 km d 1 , which has been recorded for the closely related Longtailed Skua Stercorarius longicaudus [37], between the breeding site and the first or last position estimate.Rapid movement to or away from the wintering area, roughly entering or leaving the boxes shown in figure 2, was taken as the timing of arrival at and departure from the wintering area.The date of arrival in the North Atlantic during spring migration was taken as the day birds crossed 35°N.For birds wintering in the Mediterranean Sea, the day after wintering area departure was taken.Clutch initiation dates were inferred from the light data, as incubation stints are apparent as alternating periods of light and darkness each lasting more than an hour [38].
Based on positions in the wintering areas, tracks were assigned to one of seven discrete areas: 1) Caribbean region, Guiyana, Gulf of Mexico and north of Brazil; 2) Patagonian Shelf; 3) Benguela region west to Tris tan da Cuñha; 4) Gulf of Guinea and west towards eastern Brazil; 5) Canary Current and off northwest Africa/Iberia; 6) Mediterranean Sea; 7) Indian Ocean and adjacent seas (figure 2).Although skuas moved during the wintering period, sometimes over substantial distances, none had arrival and departure positions in different wintering areas.

Statistical analysis
Wintering centroids were calculated as the geographical mean of positions between arrival and departure.
Minimum migration distance was calculated as the greatcircle distance between the wintering centroid and the breeding site.As a measure of migratory connectivity, the degree to which individuals from the same breeding site stay together at wintering areas [30], we calculated the Mantel correlation   using the 'MigConnectivity' package in R [39], using the first track of each individual.Significance of   was assessed by 1000 random permutations.
Timing of key events and duration of periods were modeled as a function of breeding latitude, with ran dom intercepts for each combination of breeding and wintering area as well as for individuals, in Bayesian Generalized Linear Mixedeffects Models (GLMMs).
We then assessed whether wintering in certain areas led to earlier or later arrival compared to birds from the same breeding area, by constructing GLMMs with breeding and wintering area as fixed effects and random intercepts for individuals.Data with at least three birds for each combination of breeding and wintering areas were selected.We performed two analyses to show the effect of wintering area.First, all pairwise comparisons of wintering area effects within breeding areas were made based on the posterior distributions of the contrast between each set of two wintering areas.Posterior contrast probabilities where 0 was outside the 95% high density interval (HDI) were considered important.Second, to test whether wintering areas led to advances or delays relative to birds from the same breeding area but wintering elsewhere, the 95% HDI of the posterior distributions of each breedingwintering area combination were compared to the region of practical equivalence (ROPE), defined as the mean timing of birds from the same breeding area but wintering elsewhere, ± 0.1 * sd [40].If the 95% HDI of the posterior distribution fell completely outside the ROPE, the wintering area was taken to have an important effect.Variance components of the fixed and random effects were used to show partitioning of variation between and within individuals (see supplementary material).
Using only data from breeding areas with the largest data sets (FAR, NOR and SVA), we tested whether the duration of periods between events changed according to timing of the onset of each period, by fitting GLMMs with duration as the response, the onset of the period as a fixed effect, random intercepts and slopes for breedingwintering area combinations and random intercepts for individuals.To test whether the duration of periods also changed as a function of timing within each breedingwintering area combination, GLMMs were constructed with duration as the response, slopes with timing per breedingwintering area as fixed effects and random intercepts for individuals.GLMMs were fitted using the brms package in R [41], which provides an interface to stan [42], with 6000 iterations, a burnin of 2000, a thinning interval of five and default uninformative priors.Model results were checked for mixing of and autocorrelation within the chains, and for convergence ( R ~1).

Migratory connectivity and spatial consistency of individuals
For any given colony, there was a range of wintering areas used by individual Arctic Skuas (figure 2).The largest variation was among individuals from Slettnes, mainland Norway (largest sample size, n=46 individ uals), which wintered in all seven areas identified in this study.Variation among individuals was smallest in the small samples from Russian sites that all wintered in the Indian Ocean (n=6), and from Greenland, that all wintered in the Caribbean region (n=3).Arctic Skuas breeding in Svalbard (n=42) wintered primarily in the Caribbean region (51%), and none migrated to the southern wintering areas, the Patagonian Shelf and Benguela region.In contrast to the Svalbard population, few individuals from mainland Norway wintered in the Caribbean region (9%, n=57).A considerable proportion of individuals breeding in mainland Norway and Finland wintered on the Patagonian Shelf (43%, n=63).Birds breeding in Iceland, the Faroe Islands and Scotland wintered mainly in the Benguela region (31%), Canary Current (36%) and the Patagonian shelf (23%, n = 39), but none in the Caribbean or Mediterranean Sea.Individual winteringsite fidelity was very high (supplementary material).
The considerable degree of mixing of individuals from different breeding sites in wintering areas indicated weak overall migratory connectivity, as is expressed in the   of 0.19 (95% HDI: 0.11 -0.26); exclud ing the birds from the Russian sites resulted in a slightly lower   of 0.14 (95% HDI: 0.07 -0.22).
During both autumn and spring migration, most birds wintering in the Atlantic staged in a broad area across the central North Atlantic, roughly between 35°N and 55°N (figure 2).Individuals breeding at Svalbard and wintering in the Canary Current, however, skipped this staging area in spring and followed a more eastern route along Western Europe [36].

Annual schedules in relation to migration strategies
Birds at more northerly colonies started breeding (i.e.laid the first egg) later (figure 3, supplementary mate rial, figure S1) and also started and finished migration later.The duration of spring migration, the time spent in the North Atlantic and the time between return to the breeding area and clutch initiation all decreased with breeding latitude (figure 3).However, the duration of autumn migration, wintering and breeding periods did not change with breeding latitude (figure 3).
Birds that used different wintering areas also showed differences in timing and duration of migration (figure 4).
There was a clear trend of longer migration periods and shorter wintering period when wintering further south for all breeding areas.Birds wintering on the Patagonian Shelf and the Benguela region (the two furthest wintering areas, at 1264714664 km and 980311800 km, respectively, from colonies) generally stayed for shorter periods in the wintering areas and had longer spring migrations than individuals wintering in the Canary Current (the closest wintering area, at 40217140 km from colonies), with the Guinean and Caribbean regions showing intermediate duration of migration and wintering periods.Birds from the same breeding population returned to the North Atlantic at the same time regardless of wintering area, with the exception of later timing among birds wintering on the Patagonian Shelf, and earlier timing among Norwegian birds wintering in the Canary Current.In all breeding areas except Norway, the timing of return to the breeding area was similar between birds from different wintering areas.Among Norwegian birds, we found weak evidence that return to the breeding grounds was slightly delayed among birds that wintered on the Patagonian Shelf and was advanced among those that wintered in the Canary Current.Although timing of clutch initiation varied among Norwegian birds -birds that wintered in the Caribbean region bred substantially earlier -there was no relation between wintering area and clutch initiation date across breeding areas.From wintering area arrival to clutch initiation in the next breeding season, variation between and within individuals declined (electronic supplementary material, figure S2).

Strength of temporal carryover effects
Both between and within breeding populations, the duration of each migration phase was shorter when started later (figure 5).There was no consistent pattern among breeding areas in the steepness of this relationship, and slope estimates (a measure of the strength of the carryover effect) were not consistently associated with breeding or wintering areas (electronic supplementary material, figure S3).Overall, the decline of duration with time was steepest for the wintering period (electronic supplementary material, figure S3b), when slope estimates around 1 indicate complete compensation of the timing of arrival in the wintering area (figure 1).There were two exceptions to this pattern: slope estimates around 0.5 for the two farthest wintering areas suggested only partial compensation for late arrival in a) Norwegian birds wintering at the Patagonian Shelf ( = 0.5, 95% HDI = 0.83 -0.15) and b) Faroese birds wintering in the Benguela region ( = 0.48, 95% HDI = 1.2 -0.24).Note, however, that the estimate for the Faroese birds wintering in the Benguela region was associated with a very wide HDI overlapping with both 1 and 0, and that the slope estimate for birds from mainland Norway wintering in the Benguela region indicates complete compensation ( = 0.97, 95% HDI = 1.53 -0.43).The duration of spring migration and the time spent in the North Atlantic were tightly related to their timing, with overlapping distributions across breeding and wintering areas (figure 5cd).The slope estimates between 0.5 and 0 indicated a partial to complete carryover effect of the arrival in the North Atlantic on return to the breeding area (figure 5de).The timing of return partially carriedover to the timing of clutch initiation, with most slope estimates between 1 and 0.5 (electronic supplementary material, figure S3e).

Discussion
Our study demonstrates substantial variation within and between breeding populations in both wintering areas and annual schedules of the Arctic Skua, a longdistance migratory seabird with a huge wintering range.Arctic Skuas showed weak overall migratory connectivity in the Atlantic, with very high variation among individuals, including individuals within breeding sites, in their use of wintering areas separated by thousands of kilometers.Individuals from the colony with the most tracked birds (Slettnes, mainland Norway) wintered in all productive areas of the Atlantic between 50°S and 30°N, and also in the Mediterranean Sea and the Persian Gulf (figure 2), together representing roughly a tenth of the world's ocean surface area.The lower variability in wintering destinations of birds from other breeding sites may partly reflect smaller sample sizes.The wintering range of Arctic Skuas is as large or larger than most other wellstudied longdistance migratory seabirds [43][44][45][46] and also far exceeds the variation shown by a congeneric species, the Long tailed Skua Stercorarius longicaudus [47].Such oceanscale variation in wintering area selection makes the Arctic Skua exceptionally wellsuited for directly testing the relationships between the timing of annual events and wintering site selection within and between breeding colonies.Migration timing of the tracked Arctic Skuas varied principally with breeding latitude, with latertimed sched ules of more northern breeders associated with faster spring migrations and shorter periods between spring arrival at the breeding grounds and clutch initiation.Variation in timing of autumn migration was completely compensated during the wintering period, even in the most distant wintering areas, where skuas spent less time.Despite the large variation in migration distances and schedules and the strong carryover effect of timing during spring migration, the timing of return to the breeding area was highly synchronized among individuals from the same site, leaving no substantial net carryover effect.
In Arctic Skuas, breeding latitude not only shifted the timing of annual schedules, as expected [12,15,16,48], but later timing also affected the duration of subsequent key periods within the annual cycle.Most notably, Arctic Skuas breeding at high latitudes had faster spring migrations and shorter prelaying periods compared to birds breeding at lower latitudes.The faster migration of latertimed spring migration of Arctic Skuas contrasts with little or no geographical variation in the duration of spring migration of two insectivorous songbirds [16,48] and a longdistance migratory seabird, the Arctic Tern Sterna paradisaea [49].Arctic Skuas may face faster seasonal changes in conditions encountered en route or may be more flexible in where they fuel for migration/reproduction, by regulating the amount of premigratory fueling in the wintering area, foraging en route and the time spent in staging areas, in particular the North Atlantic staging area where Arctic Terns stage for no or few days in spring [49].Arctic Skuas from Svalbard spent fewer days during spring migration at stopovers compared to birds from more southern breeding areas (figure 5d; [36]), which may be explained by deteriorate foraging conditions in the North Atlantic for Arctic Skuas during April to May, when the potential host species to kleptoparasitize [50], such as Arctic Terns and Sabine's Gull Xema sabini, return to their breeding areas [49,51,52].Shorter time spent at stopovers suggests that birds from more northern breeding areas need to fuel up before migration rather than depend on foraging en route to fuel the migration flight, which is common among migrating seabirds [53], and to return with large body reserves in the breeding area [54].These body stores are also used for egg production: in a population breeding at 61°N, a latitude similar to our most southern study sites and thus presumably with a relatively long prelaying period, 16% of the protein in eggs originated from distantlyacquired resources [55].Carrying large nutrient stores on migration is however costly and hence less beneficial as migratory distance increases.This tradeoff potentially explains why Arctic Skuas breeding in the high Arctic do not migrate to the southernmost wintering areas.However, carrying large nutrient stores on migration can also be advantageous as it increases flight speed when adopting the 'dynamic soaring' flight technique [56], which skuas utilize at high wind speeds.Furthermore, carrying nutrient stores is relatively less costly in largerbodied birds [19].Interestingly, Arctic Skuas from Svalbard were heavier and larger than those from mainland Norway, the Faroe Islands and Scotland (B.Moe, R.S.A. van Bemmelen, K.R.S. Snell, R.A.

Phillips unpublished data). Selection for larger body size to reduce the relative costs of carrying nutrient stores during migration provides an alternative hypothesis to explain latitudinal clines in body size, for which
Bergman's rule is usually invoked [57].
Arctic Skuas were able to compensate for the timing of the autumn migration during the wintering period.
Thereby, our study adds to a growing number of studies of ducks, waders and passerines reporting the wintering period to be the main period in the annual cycle where delays are compensated [17,27,58,59], although the migration periods may additionally act as time buffer in some species [60,61].The ability to shorten the wintering period to compensate for later timing of autumn migration has also been shown by experimentally increasing the duration of parental care [62].We investigated whether birds are less able to compensate for timing of autumn migration if the wintering period is shorter, as observed in more distant wintering areas, and evidence was inconsistent.Indeed, reduced compensation in the wintering period was found among skuas breeding in northern Norway and wintering at the furthest wintering area, the Patagonian Shelf (figure 4), where the wintering period is up to 50 d shorter than in the northernmost wintering area (but still about 50 d longer than the 122 d needed for the moult of the primary feathers; [63]).However, no conclusive evidence for reduced compensation was found among individuals wintering in the Benguela region, at a similar latitude and migration distance.Therefore, sitespecific effects, rather than migration distance alone, may explain the difference in the degree to which birds compensate earlier timing variation during the wintering period, which corresponds to findings for the Sanderling Calidris alba [13].
Our study shows largely synchronized timing of return to the breeding grounds as well as clutch initiation within breeding sites, despite differences between wintering areas in migration distance of thousands of kilometers and in average spring migration duration of up to 19 d (in Scotland).Despite the longer duration of spring migration from more distant wintering areas and the increasing strength of the carryover effect of timing during spring migration, the similar departure dates across wintering areas did not result in substantial net carryover effects on the breeding grounds.This indicates that, rather than departing earlier from the wintering grounds, Arctic Skuas with longer migrations increased travel speeds to ensure a timely return to the breeding grounds, which agrees with empirical data from songbirds [21].The similarity in timing of departure across wintering areas observed in Arctic Skuas suggests control by endogenous biological 'clocks', which are thought to be finetuned by selection [31].
Birds from widely separated wintering areas arrived back in the breeding area, on average, within less than a week.This synchrony underpins the strong selection on timing of spring arrival and breeding when competi tion for territories is high [64].Later arrival at the breeding grounds of individuals from more distant wintering areas has been reported for several short to longdistance migrants [28,[65][66][67][68][69], but only in European Shags Gulosus aristotelis and Eurasian Spoonbills Platalea leucorodia did this also lead to later reproduction (6 d in shags, 12 d in male spoonbills) and lower reproductive success.In Arctic Skuas, however, we found partial compensation for later return at the breeding grounds: on average, a 1 d later spring arrival led to 0.6 d later clutch initiation, but with substantial individual variation.We did not have data on reproductive success, but as wintering area had a minimal influence on timing of return to the breeding area (on average up to 6 d later) and laying dates (up to 3 d), a substantial effect of wintering area on reproductive success via timing of breeding seems unlikely.Furthermore, unless an effect on productivity is substantial, demonstrating a carryover effect on reproduction would be difficult considering the multitude of local effects at the breeding grounds -effects that have been shown to impact reproduction and thus drive population declines in Arctic Skuas, at least in Scotland and mainland Norway [70,71].
By describing the migratory connectivity and oceanscale geographic variation in annual schedules and carryover effects of a longdistance migratory seabird, our study contributed to our understanding of poten tial drivers of metapopulation dynamics in migratory species.Our results indicate that choice of wintering area is unlikely to affect population trends through a carryover effect on timing of breeding, irrespective of breeding latitude.Of course, this does not rule out effects on survival or reproduction via other mecha nisms.For example, as daily survival rate may differ between stages of the annual cycle [7,8], wintering area selection may affect annual survival rates and therefore population trends, given their concomitant annual schedules.While the weak migratory connectivity means that conditions in each wintering area can affect multiple breeding populations, each breeding population may be buffered from the effect of conditions at specific wintering areas due the large spread of individuals across wintering areas [30].Interestingly, the population decline in Svalbard (from where no individuals wintered in the two most distant wintering areas) is less than in the Faroe Islands, Scotland and mainland Norway [70][71][72] and the longevity record holder of the species, a bird from Finland, wintered in one of the closer wintering areas, the Canary Current [73].If and how breeding area, wintering area, and the corresponding annual schedules affect seasonal survival rates in Arctic Skuas, and whether these can explain differences in population trends, remains to be investigated.

Declarations Ethics
This study was carried out in accordance with guidelines in the use of wild birds for research [74].Catch ing and tagging Arctic Skuas was approved for the following countries by the appropriate licensing author Figure 1: Conceptual framework of the two interrelated measures of carryover effects between events in the annual cycle.The framework assumes the population mean represents the overall optimal timing, while throughout the range of timings different advantages and disadvantages can be attributed.(A) Here an example delay ( 1 ) relative to the population mean timing (  1 ) of event 1 affects the timing of a subsequent event  2 relative to the population mean   2 to different degrees.The resulting degree of advancement or delay in event ( 2 ) is the net of carryover effects and is described by its strength as partial, complete, or amplified carryover effect, or a complete compensation.The strength of the carryover effect can be quantified as the slope between  1 and  2 (B), or, between  1 and the duration of the period between  1 and  2 (C).S1).Filled squares on the right yaxes indicate that the 95% HDI was outside the ROPE of birds from the same colony wintering elsewhere.Colours correspond to the wintering areas in figure 2.
Scheme); Faroe Islands (Copenhagen Bird Ringing Centre and The National Museum of the Faroe Islands); Scotland (Special Methods Technical Panel of the British Trust for Ornithology Ringing Committee); Finland (Centre for Economic Development, Transport and the Environment, VARELY/418/07.01/2012) and Russia (Federal Law No. 52FZ of April 24, 1995, Article 44).

Figures
Figures

Figure 2 :
Figure 2: Wintering areas of Arctic Skuas tracked from breeding areas between East Greenland and West Siberia.Dots represent centroids of wintering positions for each track, coloured per wintering area (boxes).Pie charts at breeding sites represent the proportion of individuals wintering in each of the seven wintering areas, with the size corresponding to the total number of individuals tracked.Sample sizes are shown as the number of tracks / number of individuals.Green shaded areas reflect the breeding range.The main stopover in the North Atlantic is outlined in dark grey.Abbreviations for breeding sites are SVA = Svalbard, KVP = Karupelv, HOC = Hochstetter Forland, ICE = Iceland, FAR = Faroe Islands, FAI = Fair Isle, ROU = Rousay, SVA = Svalbard, SLE = Slettnes, BRH = Brensholmen, TOB = Tobseda, ERK = Erkuta and FIN = Finland, and for breeding areas are GRE = Greenland, SCO = Scotland, NOR = mainland Norway and RUS = Russia.Abbreviations for wintering areas are PAT = Patagonian Shelf, BEN = Benguela region, GUL = Gulf of Guinea, CAR = Caribbean region, CAN = Canary Current, MED = Mediterranean Sea and IND = Indian Ocean.Sites within countries were merged for statistical analyses.

Figure 3 :Figure 4 :
Figure 3: Relation between breeding latitude and timing (upper panels) or duration (lower panels) in Arctic Skuas tracked from breeding areas between East Greenland and West Siberia to wintering areas in the Atlantic.Shaded areas are 95% HDIs.

Figure 5 :
Figure 5:  The relationship between duration of each period and the timing of its start in Arctic Skuas tracked from breeding areas in the North Atlantic to wintering areas in the Atlantic.Thick lines show overall slope estimates, thin lines show regression lines per combination of breeding and wintering area.Note that the aspect ratio of all plots is 1, but the axis scales differ.