# High individual repeatability of the migratory behaviour of a long-distance migratory seabird

Nathalie Kürten, Heiko Schmaljohann, Coraline Bichet, Birgen Haest, Oscar

Vedder, Jacob González-Solís, Sandra Bouwhuis

► To cite this version:

Nathalie Kürten, Heiko Schmaljohann, Coraline Bichet, Birgen Haest, Oscar Vedder, et al. High individual repeatability of the migratory behaviour of a long-distance migratory seabird. Movement Ecology, 2022, 10 (1), pp.5. 10.1186/s40462-022-00303-y. hal-03576912

HAL Id: hal-03576912

https://hal.science/hal-03576912v1

Submitted on 16 Jul 2024

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

# High individual repeatability of the migratory behaviour of a long-distance migratory seabird

Nathalie Kürten<fnref n="1" />,²* , Heiko Schmaljohann<fnref n="1" />,², Coraline Bichet<fnref n="1" />,³, Birgen Haest⁴, Oscar Vedder<fnref n="1" />, Jacob González-Solís⁵ and Sandra Bouwhuis<fnref n="1" />

**Abstract**

**Background:** Understanding the evolution of migration requires knowledge of the patterns, sources, and consequences of variation in migratory behaviour, a need exacerbated by the fact that many migratory species show rapid population declines and require knowledge-based conservation measures. We therefore need detailed knowledge on the spatial and temporal distribution of individuals across their annual cycle, and quantify how the spatial and temporal components of migratory behaviour vary within and among individuals.

**Methods:** We tracked 138 migratory journeys undertaken by 64 adult common terns (*Sterna hirundo*) from a breeding colony in northwest Germany to identify the annual spatiotemporal distribution of these birds and to evaluate the individual repeatability of eleven traits describing their migratory behaviour.

**Results:** Birds left the breeding colony early September, then moved south along the East Atlantic Flyway. Wintering areas were reached mid-September and located at the west and south coasts of West Africa as well as the coasts of Namibia and South Africa. Birds left their wintering areas late March and reached the breeding colony mid-April. The timing, total duration and total distance of migration, as well as the location of individual wintering areas, were moderately to highly repeatable within individuals (repeatability indexes: 0.36–0.75, 0.65–0.66, 0.93–0.94, and 0.98–1.00, respectively), and repeatability estimates were not strongly affected by population-level inter-annual variation in migratory behaviour.

**Conclusions:** We found large between-individual variation in common tern annual spatiotemporal distribution and strong individual repeatability of several aspects of their migratory behaviour.

**Keywords:** Animal movement, Bird migration, Migratory behaviour, Individual consistency, Repeatability, Geolocation, Spatial ecology, Phenology

**Background**

Migration constitutes an important part of the annual cycle of many species across taxa [1, 2]. Spatiotemporal patterns, however, vary remarkably among species and populations, as well as among and within individuals

within populations [1]. Such variation should generally reflect differences in the selective ecological conditions experienced by these species, populations, or individuals [3], but detailed knowledge of the patterns, sources, and consequences of variation in migratory behaviour is needed to understand the evolution of migration [4, 5]. This need is exacerbated by the fact that many migratory species show rapid population declines (e.g. [6]), such that we are in urgent need of knowledge-based

<footnote>*Correspondence: nathalie.kuerten@ifv-vogelwarte.de</footnote>

<footnote n="1">¹ Institute of Avian Research, An der Vogelwarte 21, 26386 Wilhelmshaven,
Germany</footnote>

<footnote>Full list of author information is available at the end of the article</footnote>
conservation measures to protect decreasing populations (e.g. [7, 8]).

To understand the causes and consequences of variation in migratory behaviour, we need detailed information on the spatial and temporal distribution of individuals across their annual cycle. Ideally, this information should comprise a description of both the between-individual variation in, and within-individual consistency of, migratory behaviour. This knowledge can then be used to (i) learn whether migratory behaviour is underpinned by (epi)genetic heritability, developmental plasticity or phenotypic flexibility [4, 5, 9], (ii) assess the consequences of migratory behaviour at different timescales by quantifying seasonal, multi-seasonal or multi-year carry-over effects [10, 11] on fitness parameters (e.g. [12]), and (iii) inform or improve conservation management (see [13, 14]).

Between-individual variation in behaviour is widespread in natural populations [15, 16]. With respect to migratory behaviour, most studies pertain to birds. Studies on raptors (e.g. [17]), songbirds (e.g. [18]), geese (e.g. [19]), shorebirds (e.g. [20]), and seabirds (e.g. [21]) have reported varying levels of between-individual variation in migratory behaviour, linking it to traits such as age (e.g. [22, 23]) or sex (e.g. [24–26]). Juvenile honey buzzards (*Pernis apivorus*), for example, were found to migrate slower than adults [22] and female black-browed albatrosses (*Thalassarche melanophris*) were found to winter further north than males [25].

Variation in migratory behaviour can also occur within individuals [1]. The within-individual consistency of behaviour is often calculated using the repeatability index "R" [27], which is defined as the proportion of the total variation in a trait that can be attributed to differences between (compared to within) individuals [28]. As for the between-individual variation, the repeatability of migratory behaviour has been quantified across taxa (e.g. reptiles [29]), fish [30], and insects [31]), albeit most frequently in birds (reviewed by [32, 33]). Here, studies on raptors (e.g. [34, 35]), shorebirds (e.g. [36–38]), geese (e.g. [19, 39]), songbirds (e.g. [40, 41]) and seabirds (e.g. [42–44]) have reported repeatabilities as low as 0.03 for the duration of migration in Scopoli's shearwater (*Calonec tris diomedea*) [42], or as high as 0.99 for stopover sites (longitude and latitude) in oriental honey buzzards [35]. The observed variability in repeatabilities between traits, seasons, populations and species shows that general conclusions are hard to draw and that any study aiming at understanding the extent, causes, and consequences of variation in migratory behaviour needs to first quantify variation within and between individuals.

Here, we report on a 5-year study in which we tracked 138 migratory journeys undertaken by 64 common terns

(Sterna hirundo) from a breeding colony in northwest Germany. Although internationally the common tern is listed as being of least concern [45], in Germany it is locally endangered [46], and understanding its migratory behaviour may help to learn why. Using one to four tracks per individual, we therefore (i) identify the general spatiotemporal distribution during the annual cycle, (ii) evaluate the individual repeatability of eleven traits describing their migratory behaviour, and (iii) assess whether our repeatability estimates are affected by correcting for annual variation in migratory behaviour at the population level.

# Methods

## Study species and site

Common terns are Polarctic colonially breeding and long-distance migratory seabirds [47]. They display high breeding site fidelity (adult local return rate is ca. 90%; [48, 49]), are relatively easily caught during incubation, and are large enough to carry small tracking devices without detectable detrimental effects on their reproductive performance or survival [50]. This facilitates efficient use of such devices to monitor the (repeated) migratory behaviour of many individuals.

We studied the migratory behaviour of common terns breeding at a monospecific colony located at the Banter See in Wilhelmshaven, at the German North Sea coast (53° 30' 40" N, 08° 06' 20" E). This colony is the focus of a long-term individual-based study, in which all local fledglings have been marked with a transponder since 1992, allowing for automatic and remote detection of recruits using an antenna system (for more details see [51]). The sex of tracked birds has been molecularly determined following Becker & Wink [52].

## Deployment and recovery of light-level geolocators

Between mid-May and early July 2016–2019, we used the antenna system to identify 24, 36, 50, and 54 focal common terns of both sexes as they incubated their clutches (see Additional file 1: Table S1). We caught them using an electronically released drop trap, on average 17 days ± 3.4 SD after the first egg was laid. Before initial capture (and for birds carrying a tracking device from a previous year, on the day of laying), eggs were replaced by model eggs to prevent potential catching- (or tag-) induced damage to the real eggs [53]. Real eggs were incubated using digital incubators (Rcom max 50 and Rcom pro 20; Autoelex Co., Ltd., South Korea; [54]). The captured birds were weighed (average body mass: 129.2 g ± 7.8 SD) using a digital balance (± 1.0 g accuracy; MAULalpha, Jakob Maul GmbH, Germany) and tagged with a light-level geolocator (Intigeo-C65, Migrate Technology, UK). The geolocator was attached to the leg of the bird using a
10 mm aluminium ring. The total mass of the ring, glue and geolocator was 1.6 g, i.e. 1.2% ± 0.1 SD of the body mass of the birds at tagging and below the recommended threshold of 3% [55]. Tags did not have a detectable effect on the behaviour, reproductive performance, or survival of the birds [50]. Total handling time was 5.8 min ± 2.7 SD and all birds resumed normal incubation after being handled (i.e. no clutch was abandoned).

In the breeding seasons of 2017–2020, we recaptured 22, 29, 42, and 49 of the ‘geolocator birds’ that returned to remove the geolocator and extract the data. The other birds either did not survive or return to the colony (n=16), returned but did not attempt to reproduce, i.e. were impossible to catch (n=4), or had lost their geolocator (n=2). Out of the 142 geolocators that we retrieved from 30 males and 35 females, 105 (74%) were still recording at recapture, 2 (1%) stopped working during spring migration, 31 (22%) stopped working in the wintering area and 4 (3%) did not record any data (see Additional file 1: Table S1 for more details).

## Analysis of light-level geolocation data

We set our geolocators to sample ambient light intensity every minute and to archive the maximum light intensity every 5 min (mode 10, Migrate Technology). The retrieved light intensity data were analysed using the software R (version 4.0.3, [56]) and the R package “FLightR” (version 0.5.0, [57]), following the workflow detailed in the supplementary material from Rakhimberdiev et al. [58].

First, daily sunrise and sunset (i.e. twilight events) were identified with the function “`preprocessLight`” of the R package “`BAStag`” [59] using a light-level threshold of 1.5. Extreme outliers (e.g. > 30 min difference with the previous and subsequent twilight) were adjusted or excluded by means of a visual inspection of each individual sunrise and sunset. Next, calibration periods were individually determined by visual inspection of the output plotted using the “`plot_slopes_by_location`” function in `FLightR` [60], using the weeks the birds were known to be at the breeding colony (i.e. “on-bird” calibration). The first calibration period started after incubation (a period of high shading) and ended prior to migration. For geolocators that still recorded data at recapture, we defined a second calibration period, which started after the return to the colony and ended at the onset of incubation.

For the subsequent movement analyses, location estimates were not restricted to either land or sea (as terns may migrate over-sea and across-land [61]), but positions were spatially constrained using the function “`make.grid`” to the area between 40° W, 45° S, 45° E and 65° N to adhere to ring recovery data [53]. We also limited the maximum flight distance between twilights to

1200 ± 300 km, i.e. 24 h times 50 km/h [62]. Finally, the “`run.particle.filter`” function was set to 1 million particles to optimize the track of each individual and its uncertainty [60] (the R code of the analysis is provided in the supplementary information – see Additional file 2).

## Defining migration traits

We defined each individual’s year-specific migratory and stationary periods (stopover site(s) and wintering area(s)) using the “`stationary.migration summary`” function of “`FLightR`”. To detect migratory periods, we set the minimum probability that defines movement (“`prob.cutoff`”) to 0.4 (see Additional file 1: Tables S2 and S3). Stationary periods are assigned by the function if an individual stays in a given area for a minimum number of days, set using the argument “`min.stay`”. For the terns, we set “`min.stay`” to ten twilights, i.e. five days, as common terns moving along the East Atlantic Flyway are thought to use stopover sites for more than five consecutive days [53, 63]. Detected stopover sites within a range of 250 km around the breeding colony were discarded (n=10 stopover sites of 9 individuals) to account for potential movements within the breeding period, including local movement to nearby locations used after breeding and loafing (i.e. after breeding failure, but before migration).

From the output provided by the “`stationary.migration summary`” function, we extracted the longitude and latitude of each stopover site and each wintering area. Due to an uncertainty of dates of arrival to, and departure from, the stopover site(s) (i.e. no length of stay at stopover site(s) calculable; see Additional file 1: Tables S2 and S3), we only extracted the dates of arrival to, and departure from, the breeding colony and each wintering area. The length of stay at the wintering area(s) was calculated by subtracting the arrival date at the wintering area(s) from the departure date from the wintering area(s). Wintering area(s) were defined as sites south of 27° N (the Canary Islands, located at 28° N, are a well-known stopover area for common terns [53]) at which individuals stayed more than two months. Since birds are not totally stationary within their wintering area(s), the function detected movements in 19 tracks from 16 individuals and split their wintering area(s) into two or more small wintering areas, characterised by a small distance between them (average: 252 km; with a maximum cut-off of 500 km) and overlapping departure and arrival dates. In these cases, we manually merged the wintering areas and calculated the location of the resulting overall wintering area by using all available positions between arrival at the first, and departure from the last, wintering area. In an additional 5 tracks from 4 individuals, however, the distance between wintering areas was > 500 km and the
departure and arrival dates did not overlap, such that we retained both wintering areas (see Results).

The total migration distance (km) covered by an individual in a given year was calculated as the great circle distance (i.e. the shortest distance between two points; orthodrome) between the breeding area, any potential stopover site(s) and the wintering area(s) or vice versa using the function "distVincentyEllipsoid" of the package "geosphere" [64]. We used the great circle distance instead of the rhumb line (loxodrome), because (i) the inherent error associated with the measurement is the same for every bird, and (ii) there was only a marginal difference between these two calculations (0.47% and 0.40% for total migration distance in autumn and spring, respectively). For the 5 cases in which 4 individuals had more than one wintering area, the total migration distance was calculated to the first wintering area for autumn migration and from the last wintering area for spring migration.

We considered two different estimates of migration duration and speed. The total duration of migration (days) was estimated by subtracting the arrival date from the departure date during autumn and spring migration. The actual duration of migration (days) was estimated similarly, but using only tracks of birds for which no stopover sites were detected. The total speed of migration (km/day) was calculated by dividing the total migration distance by the total duration of migration. The (minimum) travel speed (km/day) was calculated by dividing the total migration distance by the actual duration of migration.

Migration schedules were summarized using mean values of the longitude and latitude for location estimates, mean values for estimates of the length of stay at the wintering area(s) (days), the dates of departure and arrival, total migration distance (km), total duration of migration (days), actual duration of migration (days), total speed of migration (km/day) and travel speed (km/day), and the range for the number of stopover sites.

## Statistical analyses

For autumn migration, our dataset comprised 138 tracks from 64 birds. For the wintering area, where 31 geolocators stopped working (see above), we could still analyse the longitude and latitude using all 138 tracks, as our estimates of the repeatability of longitude and latitude did not differ when we excluded tracks based on data for less than two months (see Additional file 1: Table S4). For the length of stay and departure date from the wintering area, however, our dataset comprised 107 tracks from 60 birds. For the remaining traits during spring migration, our dataset comprised 105 tracks from 60 birds.

## Spatiotemporal distribution

We first tested whether migratory traits differed between autumn and spring migration and/or between the sexes by running a set of (generalised) linear mixed models that included season and sex as two 2-level categorical fixed effects. We also included the interaction between sex and season, but removed it if not significant to allow for a more straightforward interpretation of the main effects. All models included year as a categorical fixed effect and individual identity as a random intercept to account for between-year variation and the non-independence of repeated tracks from the same birds, respectively. Variation in total migration distance, total duration of migration, actual duration of migration, total speed of migration and travel speed were analysed using five linear mixed models, since log transformation of the data facilitated the fit of a normal error distribution. Variation in stopover probability and the number of stopover sites were analysed using two generalized linear mixed models with a binomial error distribution and a "logit" link function, and a Poisson distribution and a "log" link function, respectively.

Second, we tested for sex differences in the departure date from the colony, arrival date at the (first) wintering area, departure date from the (last) wintering area, arrival date at the colony, longitude and latitude of the (first) wintering area (all not season-dependent) using six linear mixed models with a normal error distribution. These models also included year as a categorical fixed effect and individual identity as a random intercept.

The models were run with the functions "lmer" and "glmer" of the R package "lme4" [65]. P-values for linear mixed models were obtained using the "lmerTest" package [66], those for generalised linear mixed models were extracted from the "glmer" summary output. The level of significance was set to *p* < 0.05 and parameter estimates are given as mean ± SE.

## Repeatability of the spatiotemporal distribution

The intra-individual repeatability "R" was calculated for spring and autumn total migration distance, spring and autumn total duration of migration, autumn departure date from the breeding colony, autumn arrival date at the (first) wintering area, longitude of the (first) wintering area, latitude of the (first) wintering area, length of stay at the (first) wintering area, spring departure date from the (last) wintering area and spring arrival date at the breeding colony using eleven models, each fitted with a normal error distribution, year and sex as categorical fixed effects and individual identity as a random intercept. In addition, we calculated the repeatability without accounting for year as a fixed effect. While analyses "with
year" provide information on the repeatability of the relative trait expression of an individual in comparison to its conspecifics (e.g. whether an individual generally departs earlier than its conspecifics), analyses "without year" provide information on the repeatability of its absolute trait expression (e.g. whether an individual generally departs on a similar day of the year, while its conspecifics depart on another day of the year).

We also tried to calculate the intra-individual repeatability of stopover probability during autumn and spring migration using models fitted with a binominal error distribution and "logit" link function. These models again included (year and) sex as (a) categorical fixed effect(s) and individual identity as a random intercept. Since these models, however, did not produce reliable output (see Additional file 1: Table S5), we instead investigated the output provided by the "stationary.migration summar" function of individuals for which we had full data for two or more years, and report the number of individuals that showed consistent or inconsistent stopover behaviour.

All repeatability models were run using the "rpt" and the "rptBinary" functions of the R package "rptR" [67], which provided mean (±SE) R estimates, 95% confidence intervals obtained from 1000 bootstrap iterations and p-values. We visualised our data using the R package "ggplot2" [68] and the program QGIS [69] making use of the "heatmap" function (quartic kernel with a radius of 2 degrees, and an output grid size of 0.1 degrees; Figs. 1, 3).

# Results

## Spatiotemporal distribution

With respect to the *average* spatiotemporal distribution, common terns left the breeding colony in the northwest of Germany on 5 September (range 24 July–1 October), with females leaving earlier than males (Table 1, pt. A). Birds of both sexes moved south along the East Atlantic Flyway, flying both over-sea and across-land (Fig. 1). In 49 tracks of 26 individuals, 1–2 autumn stopover sites at the coasts of France, Portugal, Morocco, Western Sahara or Mauretania, or at the Canary Islands were used (see Additional file 1: Table S6).

Birds arrived at their wintering areas on 18 September (range 5 September–15 November), with females arriving earlier than males (Table 1, pt. B). These wintering areas were mainly located at (i) the west coast of West Africa, (ii) the south coast of West Africa and (iii) the coast of Namibia and South Africa (Fig. 1; see also Additional file 1: Table S7).

The total autumn migration distance was 5281 km (range 3974–11,027). This distance was covered in 12 days (range 2–67) with a total speed of 602 km/day (range 81–2231). The actual duration of the period

in which birds covered the distance was 7 days (range 2–16), the travel speed 742 km/day (range 323–2231).

Most birds used a single wintering area in which they stayed for 192 days (range 116–247; n=103 tracks of 58 individuals). Four birds showed an exception to this pattern in some years by using two wintering areas (distance: 877 km; range 616–1328; n=5 tracks of 4 individuals). For the subset of birds for which we could also quantify the length of their stay at their multiple wintering areas, this was 101 days (range 63–130; n=4 tracks of 3 individuals) for the first and 105 days (range 66–146; n=3 tracks of 2 individuals) for the second wintering area. With respect to the wintering distribution, we found no significant difference between the sexes (Table 1, pt. C and D). Furthermore, while some pair members overwintered in the same wintering area (n=3 couples), others went to different wintering areas (n=4 couples).

<figure><img src="image_3.png" /><figcaption>**Fig. 1** Migration routes and wintering areas of 64 common terns tracked with a light-level geolocator. Wintering areas are based on the individuals' locations (estimated using FLightR) at daily twilights from one week after the estimated arrival up to one week before the estimated departure from the wintering area. For the birds of which the geolocator stopped working at the wintering area, data were used until the last estimated position. The heatmaps of the twilight positions were produced in QGIS using quartic kernel density with a 2° radius</figcaption></figure>
<table><caption>**Table 1** Summary of models testing for effects of sex (female as a reference) on (A) departure from colony, (B) arrival at wintering area, (C) wintering longitude, (D) wintering latitude, (E) wintering latitude and (F) arrival at colony of common terns deployed with geolocators during 2016–2020</caption><thead><tr><th>Model</th><th>estimate ± SE</th><th>df</th><th>t-value</th><th>p-value</th></tr></thead><tbody><tr><td colspan="5"><strong>A. departure from colony (days)<sup>a</sup></strong></td></tr><tr><td>Intercept</td><td>233.73 ± 2.70</td><td>120.49</td><td>86.64</td><td><strong>&lt;0.001</strong></td></tr><tr><td>year_2017/18</td><td>4.16 ± 2.15</td><td>76.01</td><td>1.93</td><td>0.057</td></tr><tr><td>year_2018/19</td><td>12.82 ± 2.10</td><td>79.59</td><td>6.11</td><td><strong>&lt;0.001</strong></td></tr><tr><td>year_2019/20</td><td>11.36 ± 2.17</td><td>86.55</td><td>5.24</td><td><strong>&lt;0.001</strong></td></tr><tr><td>sex</td><td>13.19 ± 3.09</td><td>73.58</td><td>4.27</td><td><strong>&lt;0.001</strong></td></tr><tr><td colspan="5"><strong>B. arrival at wintering area (days)<sup>a</sup></strong></td></tr><tr><td>Intercept</td><td>243.86 ± 3.64</td><td>125.43</td><td>67.09</td><td><strong>&lt;0.001</strong></td></tr><tr><td>year_2017/18</td><td>7.13 ± 3.18</td><td>79.23</td><td>2.25</td><td><strong>0.028</strong></td></tr><tr><td>year_2018/19</td><td>12.60 ± 3.08</td><td>83.93</td><td>4.09</td><td><strong>&lt;0.001</strong></td></tr><tr><td>year_2019/20</td><td>14.75 ± 3.16</td><td>92.54</td><td>4.67</td><td><strong>&lt;0.001</strong></td></tr><tr><td>sex</td><td>13.60 ± 3.95</td><td>70.80</td><td>3.44</td><td><strong>&lt;0.001</strong></td></tr><tr><td colspan="5"><strong>C. wintering area longitude<sup>a</sup></strong></td></tr><tr><td>Intercept</td><td>-11.97 ± 1.17</td><td>1.17</td><td>69.83</td><td><strong>&lt;0.001</strong></td></tr><tr><td>year_2017/18</td><td>0.07 ± 0.13</td><td>0.13</td><td>70.01</td><td>0.555</td></tr><tr><td>year_2018/19</td><td>0.02 ± 0.12</td><td>0.12</td><td>70.06</td><td>0.875</td></tr><tr><td>year_2019/20</td><td>0.02 ± 0.13</td><td>0.13</td><td>70.12</td><td>0.131</td></tr><tr><td>sex</td><td>0.17 ± 0.57</td><td>0.57</td><td>78.94</td><td>0.769</td></tr><tr><td colspan="5"><strong>D. wintering area latitude<sup>a</sup></strong></td></tr><tr><td>Intercept</td><td>13.94 ± 1.55</td><td>92.06</td><td>9.02</td><td><strong>&lt;0.001</strong></td></tr><tr><td>year_2017/18</td><td>-1.36 ± 0.46</td><td>70.11</td><td>-2.96</td><td><strong>0.004</strong></td></tr><tr><td>year_2018/19</td><td>-0.41 ± 0.45</td><td>70.49</td><td>-0.90</td><td>0.370</td></tr><tr><td>year_2019/20</td><td>-1.19 ± 0.48</td><td>71.06</td><td>-3.99</td><td><strong>&lt;0.001</strong></td></tr><tr><td>sex</td><td>-1.18 ± 1.68</td><td>130.31</td><td>-0.70</td><td>0.480</td></tr><tr><td colspan="5"><strong>E. departure from wintering area (days)<sup>b</sup></strong></td></tr><tr><td>Intercept</td><td>85.41 ± 3.70</td><td>101.82</td><td>23.11</td><td><strong>&lt;0.001</strong></td></tr><tr><td>year_2017/18</td><td>3.93 ± 3.61</td><td>57.57</td><td>1.09</td><td>0.281</td></tr><tr><td>year_2018/19</td><td>3.01 ± 4.29</td><td>70.68</td><td>0.70</td><td>0.485</td></tr><tr><td>year_2019/20</td><td>-1.92 ± 3.52</td><td>74.29</td><td>-0.55</td><td>0.587</td></tr><tr><td>sex</td><td>5.59 ± 3.66</td><td>63.53</td><td>1.53</td><td>0.131</td></tr><tr><td colspan="5"><strong>F. arrival at colony (days)<sup>c</sup></strong></td></tr><tr><td>Intercept</td><td>105.06 ± 1.82</td><td>99.03</td><td>57.86</td><td><strong>&lt;0.001</strong></td></tr><tr><td>year_2017/18</td><td>4.44 ± 1.45</td><td>45.87</td><td>3.07</td><td><strong>0.004</strong></td></tr><tr><td>year_2018/19</td><td>5.81 ± 1.85</td><td>54.09</td><td>3.14</td><td><strong>0.003</strong></td></tr><tr><td>year_2019/20</td><td>6.08 ± 1.48</td><td>55.07</td><td>4.11</td><td><strong>&lt;0.001</strong></td></tr><tr><td>sex</td><td>1.38 ± 2.04</td><td>69.22</td><td>0.68</td><td>0.501</td></tr></tbody></table>

Between-year variation in the dependent variables was accounted for by adding year as a fixed effect to the models (2016/2017 as a reference). P-values ≤ 0.05 are presented in bold

<table><caption>**Table 2** Summary of models testing for effects of season (spring as a reference) and sex (female as a reference) on (A) total migration distance, (B) total duration of migration, (C) actual duration of migration, (D) total speed of migration and (E) travel speed of common terns deployed with geolocators during 2016–2020</caption><thead><tr><th>Model</th><th>estimate ± SE*</th><th>df</th><th>t-value</th><th>p-value</th></tr></thead><tbody><tr><td colspan="5"><strong>A. total migration distance (km)<sup>a</sup></strong></td></tr><tr><td>Intercept</td><td>8.51 ± 0.04</td><td>102.20</td><td>237.88</td><td><strong>&lt;0.001</strong></td></tr><tr><td>season</td><td>0.03 ± 0.01</td><td>174.30</td><td>3.19</td><td><strong>0.002</strong></td></tr><tr><td>year_2017/18</td><td>0.03 ± 0.01</td><td>175.60</td><td>2.08</td><td><strong>0.039</strong></td></tr><tr><td>year_2018/19</td><td>0.01 ± 0.01</td><td>178.10</td><td>0.55</td><td>0.583</td></tr><tr><td>year_2019/20</td><td>0.05 ± 0.01</td><td>180.00</td><td>3.64</td><td><strong>&lt;0.001</strong></td></tr><tr><td>sex</td><td>0.01 ± 0.04</td><td>167.10</td><td>0.18</td><td>0.859</td></tr><tr><td colspan="5"><strong>B. total duration of migration (days)<sup>a</sup></strong></td></tr><tr><td>Intercept</td><td>2.16 ± 0.13</td><td>169.59</td><td>17.16</td><td><strong>&lt;0.001</strong></td></tr><tr><td>season</td><td>0.55 ± 0.07</td><td>182.36</td><td>8.31</td><td><strong>&lt;0.001</strong></td></tr><tr><td>year_2017/18</td><td>0.08 ± 0.11</td><td>195.17</td><td>0.74</td><td>0.460</td></tr><tr><td>year_2018/19</td><td>0.05 ± 0.11</td><td>210.31</td><td>0.47</td><td>0.643</td></tr><tr><td>year_2019/20</td><td>0.27 ± 0.11</td><td>223.66</td><td>2.47</td><td><strong>0.014</strong></td></tr><tr><td>sex</td><td>-0.04 ± 0.13</td><td>69.97</td><td>-0.33</td><td>0.741</td></tr><tr><td colspan="5"><strong>C. actual duration of migration (days)<sup>b</sup></strong></td></tr><tr><td>Intercept</td><td>1.82 ± 0.12</td><td>96.30</td><td>17.14</td><td><strong>&lt;0.001</strong></td></tr><tr><td>season</td><td>0.28 ± 0.09</td><td>106.84</td><td>3.21</td><td><strong>0.001</strong></td></tr><tr><td>year_2017/18</td><td>-0.09 ± 0.11</td><td>103.66</td><td>-0.75</td><td>0.454</td></tr><tr><td>year_2018/19</td><td>-0.05 ± 0.12</td><td>110.90</td><td>-0.42</td><td>0.673</td></tr><tr><td>year_2019/20</td><td>0.09 ± 0.11</td><td>113.67</td><td>0.85</td><td>0.396</td></tr><tr><td>sex</td><td>0.16 ± 0.08</td><td>38.13</td><td>1.92</td><td>0.063</td></tr><tr><td colspan="5"><strong>D. total speed of migration (km/day)<sup>a</sup></strong></td></tr><tr><td>Intercept</td><td>6.38 ± 0.10</td><td>183.03</td><td>61.82</td><td><strong>&lt;0.001</strong></td></tr><tr><td>season</td><td>-0.52 ± 0.06</td><td>187.88</td><td>-8.38</td><td><strong>&lt;0.001</strong></td></tr><tr><td>year_2017/18</td><td>-0.05 ± 0.10</td><td>206.74</td><td>-0.47</td><td>0.637</td></tr><tr><td>year_2018/19</td><td>-0.06 ± 0.11</td><td>225.32</td><td>-0.60</td><td>0.547</td></tr><tr><td>year_2019/20</td><td>-0.22 ± 0.10</td><td>235.23</td><td>-2.27</td><td><strong>0.024</strong></td></tr><tr><td>sex</td><td>0.01 ± 0.10</td><td>63.16</td><td>0.12</td><td>0.902</td></tr><tr><td colspan="5"><strong>E. travel speed (km/day)<sup>b</sup></strong></td></tr><tr><td>Intercept</td><td>6.59 ± 0.10</td><td>89.08</td><td>68.91</td><td><strong>&lt;0.001</strong></td></tr><tr><td>season</td><td>-0.29 ± 0.08</td><td>112.05</td><td>-3.51</td><td><strong>&lt;0.001</strong></td></tr><tr><td>year_2017/18</td><td>0.10 ± 0.11</td><td>110.38</td><td>0.94</td><td>0.351</td></tr><tr><td>year_2018/19</td><td>0.05 ± 0.11</td><td>113.81</td><td>0.41</td><td>0.680</td></tr><tr><td>year_2019/20</td><td>-0.06 ± 0.10</td><td>110.34</td><td>-0.62</td><td>0.536</td></tr><tr><td>sex</td><td>-0.12 ± 0.07</td><td>34.55</td><td>-1.80</td><td>0.080</td></tr></tbody></table>

Between-year variation in the dependent variables was accounted for by adding year as a fixed effect to the models (2016/2017 as a reference). P-values ≤ 0.05 are presented in bold

Birds of both sexes left their wintering areas on 29 March (range 10 February–25 May) (Table 1, pt. E). They travelled a total distance of 5463 km (range 3497–11,447)

to their breeding colony, i.e. further than during autumn migration (Table 2, pt. A). This distance was covered in 22 days (range 4–99) with a total speed of 358 km/day

<footnote>a n = 138 tracks of 64 individuals (29 ♂ + 35 ♀)</footnote>

<footnote>b n = 107 tracks of 60 individuals (27 ♂ + 33 ♀)</footnote>

<footnote>c n = 105 tracks of 60 individuals (27 ♂ + 33 ♀)</footnote>

<footnote>a n = 243 (138 + 105) tracks of 64 individuals</footnote>

<footnote>b n = 120 (89 + 31) tracks of 47 individuals</footnote>

<footnote>c Data were log10-transformed</footnote>
(range 61–1277), such that birds travelled longer and more slowly in spring than in autumn (Table 2, pt. B and D). Similarly, the actual duration of spring migration, at 9 days (range 4–23), was longer than that in autumn, and covered with a lower travel speed of 560 km/day (range 209–1277) (Table 2, pt. C and E).

In 74 tracks of 48 individuals, 1–4 spring stopover sites were used at the coasts of Mauretania, Western Sahara, Morocco or Portugal, at the Canary Islands or in the Bay of Biscay (see Additional file 1: Table S6). The number of stopover sites did not differ between spring and autumn migration ($0.25 \pm 0.16$ SE, $z = 1.56, p = 0.119$; $n = 123$ tracks of 53 individuals), but the probability that birds used stopover sites was higher in spring than in autumn ($2.24 \pm 0.42$ SE, $z = 5.35, p < 0.001$; $n = 243$ tracks of 64 individuals).

Spring migration ended on 20 April (range 1 April–8 May), and arrival at the breeding colony did not differ between the sexes (Table 1, pt. F).

## Repeatability of the spatiotemporal distribution

Both the annual departure date from the breeding colony and arrival date at the wintering area were highly repeatable ($R = 0.72, p < 0.001$ and $R = 0.63, p < 0.001$; Fig. 2a, d; Table 3, pt. A), such that the total duration of autumn migration was highly repeatable too ($R = 0.66, p < 0.001$; Fig. 2c; Table 3, pt. A). During autumn migration, most birds (78%) showed consistent stopover behaviour, with 22 birds stopping in each of their 60 tracks and 9 birds never stopping in any of their 27 tracks. Only 9 birds were inconsistent in whether or not they stopped in their 27 tracks.

The total migration distance ($R = 0.94, p < 0.001$; Fig. 2b; Table 3, pt. A), and longitude and latitude of the wintering areas were highly repeatable ($R = 1.00, p < 0.001$ and $R = 0.98, p < 0.001$; Fig. 2e, f; Fig. 3a–f; Table 3, pt. A), although the individual migration routes showed variability between years (see Additional file 1: Fig. S1). The individual length of stay at the wintering area was highly repeatable too ($R = 0.64, p < 0.001$; Fig. 2g; Table 3, pt. A).

The timing of spring migration was moderately repeatable with respect to the departure from the wintering area ($R = 0.36, p < 0.001$; Fig. 2h; Table 3, pt. A), although arrival date at the breeding colony and total duration of spring migration were again highly repeatable ($R = 0.75, p < 0.001$ and $R = 0.65, p < 0.001$; Fig. 2k, j; Table 3, pt. A).

As in autumn, the individual total spring migration distance was highly repeatable ($R = 0.93, p < 0.001$; Fig. 2i; Table 3, pt. A). Spring stopover behaviour was consistent in 63% of birds, since 6 birds stopped in each of their 13 tracks, 13 birds never stopped in any of their 33 tracks, and 11 birds were inconsistent in whether they stopped or not (29 tracks).

Removing year as a fixed effect from the models led to minor or no reductions in the repeatability estimates for migration distance, the spatial distribution of the wintering areas or the phenology of spring migration, but to somewhat stronger reductions of the repeatability estimates for the phenology of autumn migration (Table 3, pt. B).

## Discussion

Understanding the evolution of migration requires knowledge of the patterns, sources, and consequences of variation in migratory behaviour. Moreover, since many populations of migratory species show rapid population declines (e.g. [6]), with populations of long-distance migrants declining faster than those of short-distance migrants, at least among birds (e.g. [70]), this knowledge also is crucial for facilitating implementation of knowledge-based conservation measures (e.g. [7, 8]). As a first step towards understanding the causes and consequences of variation in migratory behaviour of a locally endangered species, the common tern, we tracked 138 migratory journeys undertaken by 64 individual birds. We used these data to (i) identify the annual spatiotemporal distribution of these birds, (ii) show the individual repeatability of their migratory behaviour to be at the high end of the range reported for birds, and (iii) demonstrate that only the repeatability of the phenology of autumn migration is substantially affected by inter-annual variation in migratory behaviour across the population.

## Spatiotemporal distribution

The common terns from our study population in northwest Germany migrated along the East Atlantic Flyway. In autumn, females departed from the breeding colony, and arrived at the wintering areas ca. two weeks before males did – a pattern referred to as protogyny [71]. While the opposite pattern, protandry, seems to be the norm

(See figure on next page.)

**Fig. 2** Repeatability of the a departure date from the colony, b total migration distance in autumn, c total duration of autumn migration, d arrival date at the wintering area, e longitude of the wintering area, f latitude of the wintering area, g length of stay at the wintering area, h departure date from the wintering area, i total migration distance in spring, j total duration of spring migration and k arrival date at the colony. For the purpose of visualisation, we randomly selected 2 years of data per individual and plotted them against each other (a–f: $n = 38$; g–h: $n = 31$; i–k: $n = 30$). All models were, however, run using all data available for each individual. Dotted lines represent lines of equality. JD, Julian Day
<figure><img src="image_4.png" /><figcaption>**Fig. 2** (See legend on previous page.)</figcaption></figure>
<table><caption>**Table 3** Repeatability (R) of the distance, phenology and spatial distribution of common tern migration, with and without including year as a fixed effect in the models. Note that all models included sex as a fixed effect</caption><thead><tr><th rowspan="2">Trait</th><th colspan="3">A. sex- & year-specific repeatability</th><th colspan="3">B. sex-specific repeatability</th></tr><tr><th>R ± SE</th><th>95% CI</th><th>p-value</th><th>R ± SE</th><th>95% CI</th><th>p-value</th></tr></thead><tbody><tr><td><strong>migration distance</strong></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>total migration distance in autumn<sup>a</sup></td><td>0.937 ± 0.015</td><td>0.906–0.963</td><td><strong>&lt;0.001</strong></td><td>0.934 ± 0.017</td><td>0.895–0.959</td><td><strong>&lt;0.001</strong></td></tr><tr><td>total migration distance in spring<sup>c</sup></td><td>0.925 ± 0.021</td><td>0.881–0.961</td><td><strong>&lt;0.001</strong></td><td>0.907 ± 0.025</td><td>0.851–0.947</td><td><strong>&lt;0.001</strong></td></tr><tr><td><strong>phenology</strong></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>departure date from colony<sup>a</sup></td><td>0.719 ± 0.063</td><td>0.579–0.827</td><td><strong>&lt;0.001</strong></td><td>0.571 ± 0.081</td><td>0.391–0.717</td><td><strong>&lt;0.001</strong></td></tr><tr><td>arrival date at wintering area<sup>a</sup></td><td>0.630 ± 0.080</td><td>0.473–0.764</td><td><strong>&lt;0.001</strong></td><td>0.533 ± 0.090</td><td>0.318–0.692</td><td><strong>&lt;0.001</strong></td></tr><tr><td>total duration of autumn migration<sup>a</sup></td><td>0.661 ± 0.072</td><td>0.503–0.795</td><td><strong>&lt;0.001</strong></td><td>0.630 ± 0.080</td><td>0.448–0.759</td><td><strong>&lt;0.001</strong></td></tr><tr><td>length of stay at wintering area<sup>b</sup></td><td>0.638 ± 0.089</td><td>0.451–0.796</td><td><strong>&lt;0.001</strong></td><td>0.598 ± 0.095</td><td>0.392–0.751</td><td><strong>&lt;0.001</strong></td></tr><tr><td>departure date from wintering area<sup>b</sup></td><td>0.358 ± 0.130</td><td>0.100–0.618</td><td><strong>&lt;0.001</strong></td><td>0.368 ± 0.130</td><td>0.084–0.611</td><td><strong>&lt;0.001</strong></td></tr><tr><td>arrival date at colony<sup>c</sup></td><td>0.745 ± 0.064</td><td>0.617–0.872</td><td><strong>&lt;0.001</strong></td><td>0.713 ± 0.075</td><td>0.550–0.837</td><td><strong>&lt;0.001</strong></td></tr><tr><td>total duration of spring migration<sup>c</sup></td><td>0.646 ± 0.085</td><td>0.492–0.804</td><td><strong>&lt;0.001</strong></td><td>0.662 ± 0.085</td><td>0.475–0.806</td><td><strong>&lt;0.001</strong></td></tr><tr><td><strong>spatial distribution</strong></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>wintering area longitude<sup>a</sup></td><td>0.998 ± 0.001</td><td>0.997–0.999</td><td><strong>&lt;0.001</strong></td><td>0.998 ± 0.001</td><td>0.997–0.999</td><td><strong>&lt;0.001</strong></td></tr><tr><td>wintering area latitude<sup>a</sup></td><td>0.979 ± 0.005</td><td>0.970–0.988</td><td><strong>&lt;0.001</strong></td><td>0.973 ± 0.007</td><td>0.958–0.983</td><td><strong>&lt;0.001</strong></td></tr></tbody></table>

P-values ≤ 0.05 are presented in bold

for autumn migration in birds [26], protogyny has previously been observed in common terns ([72–75], but see 53]), as well as in other migratory seabirds (e.g. [25, 76]), raptors (e.g. [77]), and songbirds (e.g. [71]). In common terns, it may be explained by males performing longer post-fledging parental care than females [72], a pattern that is interesting in itself, as it may lead to offspring traits developing during this phase to be more strongly influenced by the father. An example may have recently been found in oystercatchers (*Haematopus ostralegus*), in which fledglings socially inherit their migratory behaviour from their fathers [78].

We found common terns to migrate both over-sea and across-land (Fig. 1). Similar across-land migration was recently reported in tracking studies on various sea- and shorebird species, among which arctic terns (*Sterna paradisaea*) [61], lesser black-backed gulls (*Larus fus cus*) [79], reddish egrets (*Egretta rufescens*) [80], brown pelicans (*Pelecanus occidentalis*) [80], red knots (*Calidris canutus*) [80–82], as well as common terns [75]. As such, our results add to the growing evidence for the hypothesis that across-land migration is a much more widespread migration strategy for waterbirds than previously thought, which may have important implications for the protection of these birds [80].

Some terns were found to use stopover sites at the coasts and seas of Morocco, Western Sahara, Mauritania, and the Canary Islands (see Additional file 1: Table S6),

hence showing similar stopover behaviour to many other seabird species [83]. Most birds (76%) then wintered along the coast of West Africa, but some (19%) wintered further south or even as far south as Namibia and South Africa (5%) (Fig. 1). These regions are located in three different Large Marine Ecosystems: the Canary Current Large Marine Ecosystem in the west of West Africa, the Guinea Current Large Marine Ecosystem in the south of West Africa, and the Benguela Current Large Marine Ecosystem in Namibia and South Africa. These ecosystems are characterised by strong upwelling and show high rates of primary productivity (e.g. [84, 85]), thereby providing abundant food (e.g. [83, 86–88]) and being of extreme importance for a wide range of marine animals [89]. These ecosystems, however, are known to be affected by current climate change (for a review see [90]), such that consequences for primary and fish production are expected [91]. Such consequences, in turn, are likely to affect many marine animals (for a review see [92]), including seabirds (e.g. [93, 94]).

On the return journey to the breeding colony in spring, the terns showed a higher probability to use stopover sites, although the number of stopover sites used did not differ from that in autumn. We also found the total duration of migration to be longer in spring than in autumn. This pattern corresponds with that found in roseate terns (*Sterna dougallii*) wintering in West-Africa [88], but contrasts with the more general pattern of avian migration

<footnote><sup>a</sup> n = 138 tracks of 64 individuals</footnote>

<footnote><sup>b</sup> n = 107 tracks of 60 individuals</footnote>

<footnote><sup>c</sup> n = 105 tracks of 60 individuals</footnote>
<figure><img src="image_5.png" /><figcaption>**Fig. 3** Wintering areas of common terns tracked with light-level geolocators for four (a–c, n = 12) or three years (d–f, n = 11). Dots represent an individual's locations (estimated using FLightR) at daily twilights from one week after the estimated arrival up to one week before the estimated departure from the wintering area. For the birds of which the geolocator stopped working at the wintering area, data for that year were used until the last estimated position. The heatmaps of the twilight positions were produced in QGIS using quartic kernel density with a 2º radius, and visually scaled from 1 (transparent) to 300 (most intense) twilights</figcaption></figure>

being faster in spring than in autumn [95, 96]. Perhaps the mating process [97] and benefit of a short interval between arrival and egg laying to facilitate early breeding [98] exert a stronger selection pressure on the terns' condition at arrival in the breeding area than on their arrival date per se, leading them to spend more time foraging during spring than autumn migration. Alternatively, feeding conditions may be better at the breeding than wintering area, such that there is less need to forage during autumn than spring migration. Or perhaps the pattern can be explained by the prevailing wind

conditions during migration. Winds rotate clockwise in the North Atlantic, thereby offering tailwind in autumn, but headwind in spring [99], potentially leading terns, who are very susceptible to wind [100], to fly more slowly in spring than in autumn. Effects of other environmental conditions that differ consistently between the seasons, such as temperature or rainfall, could pose an additional explanation [101, 102] and future work should investigate whether and how such environmental factors influence the various aspects of migratory behaviour.
## Repeatability of the spatiotemporal distribution

We found an extremely high repeatability in common tern wintering area locations. These birds thus are highly site-faithful to both their breeding colony and individual wintering areas. Similar observations have been made in various other migratory bird species (e.g. [33, 103]), suggesting many birds to generally use a genetically inherited (see [104]) or familiar (previously encountered or socially learned) wintering area. The latter may be facilitated by spatial memory formation [105], develop through ontogeny as birds change from exploratory to more consistent migratory behaviour (e.g. [106, 107]), and come with the benefit of increased local knowledge (e.g. [108]). On the other hand, 'simply' returning to a familiar location could also pose a risk if conditions at a location were to become unsuitable but birds remained site-faithful (see [13]).

In our case, more than 75% of the terns we tracked wintered along the west coast of West Africa. A large proportion of our study population thus appears to be dependent on a relatively small geographic (wintering) region. Given this region is predicted to become less suitable with further change of the environmental conditions (for climate predictions see [109]; for fish stock predictions see [91]), if these birds do not change wintering area, we may expect them to start performing increasingly poorly, i.e. to survive less well, or to show reduced breeding performance in the case of carry-over effects to the breeding season (e.g. [110, 111]). This, in turn, could lead to (i) a local population decline if birds with another migratory phenotype (i.e. birds wintering in the other regions) are not able to compensate for such a decrease in performance, or even to (ii) a general population decline if terns from other populations use the same wintering areas [112, 113], if no new suitable wintering areas emerge, or other existing wintering areas decline in suitability too, or if other existing wintering areas cannot sustain an increasing number of birds (i.e. if there is a lack of alternatives [114]).

Contrary to the high individual repeatability in wintering area, we found migration routes to be more inconsistent (see Additional file 1: Fig. S1). Route inconsistency has been demonstrated for various migratory seabirds (e.g. [115]), shorebirds (e.g. [116]) and raptors (e.g. [117]), as well as songbirds (e.g. [118]), and suggests migration to be influenced by variable or unpredictable environmental conditions experienced *en route* (e.g. [119, 120]). Interestingly, in our study, the observed inconsistency in migration route did not result in a low repeatability of the total migration time. This suggests that the birds either managed to find what they needed by adjusting their routes to the environmental conditions, or that they prioritised their time keeping over finding what they needed. In the latter case, we would expect environmental conditions

experienced *en route* to have carry-over effects on body condition at arrival [121], which first (population-level) analyses of the study colony did not support [122], but which warrants further investigation.

In terms of phenology, departure date from the colony, arrival date at the wintering area, and arrival date at the breeding colony were highly repeatable. While consistency in the timing of migration has been demonstrated for several migratory bird species [32], the degree of repeatability in our study is at the high end of the reported spectrum. Such strong individual consistency may partly result from a genetic basis to migration phenology [123]. For one of the traits, arrival date at the breeding colony, quantitative genetic analyses of our long-term dataset (7436 first annual registrations of 1630 individuals across 1994–2019) indeed indicate an additive genetic component of c. 20% (Moiron et al. unpublished). Any non-genetic part of the repeatability in phenology should reflect what quantitative geneticists refer to as ‘permanent environment or permanent individual effects’, i.e. non-genetic stable differences between individuals for example resulting from (i) life-long carry-over effects of natal conditions (e.g. hatching date [5, 124]) or (ii) stochasticity (e.g. initial social and environmental effects) having turned into routine via reinforcement through memory models (see [105] for a spatial context, and [125] for the role of stochasticity in determining breeding phenology).

Interestingly, the repeatability of the departure date from the wintering area was only moderate repeatable and half that of the arrival date at the breeding colony (with the 95% CI overlapping only 0.1%). As such, birds relatively inconsistently timed their departure for spring migration, but then sped up or slowed down to arrive at the colony at their individual-specific date. The additional fact that the repeatability of spring arrival date itself was largely independent of general inter-annual variation (since the repeatability with or without accounting for year as a fixed effect differed by only 3.2%; Table 2) suggests there may be adjustment to keep an internal clock (synchronisation, [123]) or that a late departure, for example, is associated with more favourable conditions *en route*.

When assessing whether our estimates of the repeatability of migratory traits other than arrival date at the breeding colony were affected by general inter-annual variation, we found only the repeatability of autumn migration phenology (i.e. departure from the breeding colony and arrival at the wintering area) to be increased when correcting for year (by 15 and 10%, respectively). This increase suggests that the (tagged) population as a whole departs and arrives substantially earlier or later
in some years, while the phenological ranking of individual birds remains largely the same. A possible mechanism underlying this inter-annual population variation in autumn migration could be annual variation in breeding success, which indeed is substantial in the study population [126]. In years with poor herring abundance, common terns show strong brood reduction [127, 128]. As a result, most birds can perhaps depart (and therefore arrive) earlier in years of relatively unsuccessful reproduction, such as observed in other species (e.g. [129, 130]). Alternatively, birds may all display largely similar reaction norms to drivers of autumn migration, or there may exist a strong social influence on departure [131].

## Conclusions

Using a multi-year individual-based dataset of migrating common terns, we showed these birds to display remarkable within-individual consistency in their migratory behaviour, with repeatability values ranging from 0.36 for the departure date from the wintering area to 1.00 (longitude) and 0.98 (latitude) for the wintering area. This knowledge can be used to (i) learn whether migratory behaviour is underpinned by (epi)genetic heritability, developmental plasticity or phenotypic flexibility, (ii) assess the consequences of migratory behaviour at different timescales, and (iii) use spatiotemporal knowledge of risks individuals or populations may face to inform or improve conservation management. Hence, our study provides essential knowledge (also for future work), and exemplifies how intensive individual-based tracking studies can improve our understanding of migratory behaviour.

## Supplementary Information

The online version contains supplementary material available at https://doi.org/10.1186/s40462-022-00303-y.

**Additional file 1.** Tables S1–S7 and Figure S1.

**Additional file 2.** R code for the analysis of our light-level geolocator data.

### Acknowledgements

We would like to thank Götz Wagenknecht and Kristin Piening for their help in the field. We would also like to thank Theunis Piersma for his constructive comments, which helped us improve our manuscript.

### Authors' contributions

SB, HS and JGS initiated the study. NK and SB carried out the fieldwork. NK analysed the data with support of HS, CB, OV and SB. NK and BH visualised the data. NK drafted the manuscript with support of SB. All authors read and approved the final manuscript.

### Funding

Open Access funding was enabled and organized by Projekt DEAL. NK was supported by the German Federal Environmental Foundation (DBU) and the German Ornithologists' Society (DO-G) and JGS by the Ministerio de Economia, Industria y Competitividad of Spain (CGL2016-78530 R) and fondos FEDER.

## Availability of data and materials

The datasets underlying the results and conclusions reported in this article are provided in Additional file 1; the R code for the analysis of light-level data in Additional file 2. Light-level data are available upon request to the corresponding author.

## Declarations

### Ethics approval and consent to participate

The study was performed under licenses of the city of Wilhelmshaven and the Lower Saxony State Office for Consumer Protection and Food Safety, Germany.

### Consent for publication

Not applicable.

### Competing interests

The authors declare to have no competing interests.

### Author details

¹Institute of Avian Research, An der Vogelwarte 21, 26386 Wilhelmshaven, Germany. ²Institute of Biology and Environmental Sciences, University of Oldenburg, Carl-von-Ossietzky-Str. 9–11, 26129 Oldenburg, Germany.

³Centre d'Etudes Biologiques de Chizé, UMR 7372, CNRS-Université de La Rochelle, 79360 Villiers-en-Bois, France. ⁴Department of Bird Migration, Swiss Ornithological Institute, 6204 Sempach, Switzerland. ⁵Institut de Recerca de la Biodiversitat and Departament de Biologia Evolútiva, Ecologia i Ciències Ambientals, Universitat de Barcelona, Av. Diagonal 643, 08028 Barcelona, Spain.

Received: 2 January 2022 Accepted: 23 January 2022

Published online: 05 February 2022

## References

1. Newton I. The migration ecology of birds. London: Academic Press; 2008.

2. Fudickar AM, Jahn AE, Ketterson ED. Animal migration: an overview of one of nature's great spectacles. Annu Rev Ecol Evol. 2021;52:479–97. https://doi.org/10.1146/annurev-ecolsys-012021-031035.

3. Alerstam T, Hedenström A, Åkesson S. Long-distance migration: evolution and determinants. Oikos. 2003;103:247–60. https://doi.org/10.1034/j.1600-0706.2003.12559.x.

4. Piersma T. Flyway evolution is too fast to be explained by the modern synthesis: Proposals for an 'extended' evolutionary research agenda. J Ornithol. 2011;152:151–9. https://doi.org/10.1007/s10336-011-0716-z.

5. Verhoeven MA, Loonstra AHJ, McBride AD, Kaspersma W, Hooijmeijer JCEW, Both C, Senner NR, Piersma T. Age-dependent timing and routes demonstrate developmental plasticity in a long-distance migratory bird. J Anim Ecol. 2021. https://doi.org/10.1111/1365-2656.13641.

6. Vickery JA, Ewing SR, Smith KW, Pain DJ, Bairlein F, Škorpilová J, Gregory RD. The decline of Afro-Palaearctic migrants and an assessment of potential causes. Ibis. 2014;156:1–22. https://doi.org/10.1111/ibi.12118.

7. Verhoeven MA, Loonstra AJ, McBride AD, Both C, Senner NR, Piersma T. Migration route, stopping sites, and non-breeding destinations of adult Black-tailed Godwits breeding in southwest Fryslân, the Netherlands. J Ornithol. 2021;162:61–76. https://doi.org/10.1007/s10336-020-01807-3.

8. Strøm H, Descamps S, Ekker M, Fauchald P, Moe B. Tracking the movements of North Atlantic seabirds: steps towards a better understanding of population dynamics and marine ecosystem conservation. Mar Ecol Prog Ser. 2021;676:97–116. https://doi.org/10.3354/meps13801.

9. Charmantier A, Gienapp P. Climate change and timing of avian breeding and migration: evolutionary versus plastic changes. Evol Appl. 2014;7:15–28. https://doi.org/10.1111/eva.12126.

10. Harrison XA, Blount JD, Inger R, Norris DR, Bearhop S. Carry-over effects as drivers of fitness differences in animals. J Anim Ecol. 2011;80:4–18. https://doi.org/10.1111/j.1365-2656.2010.01740.x.

11. O'Connor CM, Norris DR, Crossin GT, Cooke SJ. Biological carryover effects: linking common concepts and mechanisms in ecology and evolution. Ecosphere. 2014;5:1–11. https://doi.org/10.1890/ES13-00388.1.
12. Rakhimberdiev E, Duijs N, Karagicheva J, Camphuysen CJ, Dekinga A, Dekker R, Gavrilov A, ten Horn J, Jukema J, Saveliev A, Soloviev M, Tibbitts TL, van Gils JA, Piersma T. Fuelling conditions at staging sites can mitigate Arctic warming effects in a migratory bird. Nat Commun. 2018;9:1–10. https://doi.org/10.1038/s41467-018-0.

13. Rakhimberdiev E, van den Hout PJ, Brugge M, Spaans B, Piersma T. Seasonal mortality and sequential density dependence in a migratory bird. J Avian Biol. 2015;46:332–41. https://doi.org/10.1111/j.00701.

14. Piersma T. Ornithology from the flatlands. Seasonal selective forces and innovations in bird migration. Ardea. 2019;107:115–7. https://doi.org/10.5253/arde.v107i2.a9.

15. Bolnick DI, Svanbäck R, Fordyce JA, Yang LH, Davis JM, Hulsey CD, Forister ML. The ecology of individuals: incidence and implications of individual specialization. Am Nat. 2003;161:1–28. https://doi.org/10.1086/343878.

16. Dall SRX, Bell AM, Bolnick DI, Ratnieks FLW. An evolutionary ecology of individual differences. Ecol Lett. 2012;15:1189–98. https://doi.org/10.1111/j.1461-0248.2012.01846.x.

17. Catry I, Dias MP, Catry T, Afanasyev V, Fox J, Frankco AMA, Sutherland WJ. Individual variation in migratory movements and winter behaviour of Iberian Lesser Kestrels *Falco naumanni* revealed by geolocators. Ibis. 2011;153:154–64. https://doi.org/10.1111/j.1474-919X.2010.01073.x.

18. Delmore KE, Fox JW, Darren EI. Dramatic intraspecific differences in migratory routes, stopover sites and wintering areas, revealed using light-level geolocators. Proc R Soc B. 2012;279:4582–9. https://doi.org/10.1098/rspb.2012.1229.

19. Bêty J, Giroux JF, Gauthier G. Individual variation in timing of migration: causes and reproductive consequences in greater snow geese (*Anser caerulescens atlanticus*). Behav Ecol Sociobiol. 2004;57:1–8. https://doi.org/10.1007/s00265-004-0840-3.

20. Tedeschi A, Sorrenti M, Bottazzo M, Spagnesi M, Telletxea I, Ibàñez R, Tormen N, De Pascalis F, Guidolin L, Rubolini D. Interindividual variation and consistency of migratory behavior in the Eurasian woodcock. Curr Zool. 2020;66:155–63. https://doi.org/10.1093/cz/zo038.

21. Brown JM, van Loon EE, Bouten W, Camphuysen KC, Lens L, Müller W, Thaxter CB, Shamoun-Baranes J. Long-distance migrants vary migratory behaviour as much as short-distance migrants: an individual-level comparison from a seabird species with diverse migration strategies. J Anim Ecol. 2021;90:1058–70. https://doi.org/10.1111/1365-2656.13431.

22. Hake M, Kjellén N, Alerstam T. Age-dependent migration strategy in honey buzzards Pernis apivorus tracked by satellite. Oikos. 2003;103:385–96. https://doi.org/10.1034/j.1600-0706.2003.12145.x.

23. Wobker J, Heim W, Schmaljohann H. Sex, age, molt strategy, and migration distance explain the phenology of songbirds at a stopover along the East Asian flyway. Behav Ecol Sociobiol. 2021;75:1–14. https://doi.org/10.1007/s00265-020-02957-3.

24. Rubolini D, Spina F, Saino N. Protandry and sexual dimorphism in trans-Saharan migratory birds. Behav Ecol. 2004;15:92–601. https://doi.org/10.1093/beheco/arh048.

25. Phillips RA, Silk JR, Croxall JP, Afanasyev V, Bennett VJ. Summer distribution and migration of nonbreeding albatrosses: individual consistencies and implications for conservation. Ecology. 2005;86:2386–96. https://doi.org/10.1890/04-1885.

26. Briedis M, Bauer S, Adamik P, Alves JA, Costa JS, Emmenegger T, Gustafsson L, Koleček J, Liechti F, Meier CM, Procházka P, Hahn S. A full annual perspective on sex-biased migration timing in long-distance migratory birds. Proc RR Soc B. 2019;286:20182821. https://doi.org/10.1098/rspb.2018.2821.

27. Bell AM, Hankison SJ, Laskowski KL. The repeatability of behaviour: a meta-analysis. Anim Behav. 2009;77:771–83. https://doi.org/10.1016/j.anbehav.2008.12.022.

28. Nakagawa S, Schielzeth H. Repeatability for Gaussian and non-Gaussian data: a practical guide for biologists. Biol Rev. 2010;85:935–56. https://doi.org/10.1111/j.1469-185X.2010.00141.x.

29. Semlitsch RD, Scott DE, Pechmann JH, Gibbons JW. Phenotypic variation in the arrival time of breeding salamanders: individual repeatability and environmental influences. J Anim Ecol. 1993;62:334–40. https://doi.org/10.2307/5364.

30. Birnie-Gauvin K, Koed A, Aarestrup K. Repeatability of migratory behaviour suggests trade-off between size and survival in a wild iteroparous salmonid. Funct Ecol. 2021;35:2717–27. https://doi.org/10.1111/1365-2435.13917.

31. Kent JW, Rankin MA. Heritability and physiological correlates of migratory tendency in the grasshopper *Melanoplus sanguinipes*. Physiol Entomol. 2001;26:371–80. https://doi.org/10.1046/j.0307-6962.2001.00257.x.

32. Both C, Bijlsma RG, Ouwehand J. Repeatability in spring arrival dates in pied flycatchers varies among years and sexes. Ardea. 2016;104:3–21. https://doi.org/10.5253/arde.v104i1.a1.

33. Phillips RA, Lewis S, González-Solis J, Daunt F. Causes and consequences of individual variability and specialization in foraging and migration strategies of seabirds. Mar Ecol Prog Ser. 2017;578:117–50. https://doi.org/10.3354/meps12217.

34. Vardanis Y, Nilsson JÄ, Klaassen RH, Strandberg R, Alerstam T. Consistency in long-distance bird migration: contrasting patterns in time and space for two raptors. Anim Behav. 2016;113:177–87. https://doi.org/10.1016/j.anbehav.2015.12.014.

35. Sugasawa S, Higuchi H. Seasonal contrasts in individual consistency of oriental honey buzzards' migration. Biol Lett. 2019;15:20190131. https://doi.org/10.1098/rsbl.2019.0131.

36. Carneiro C, Gunnarsson TG, Alves JA. Why are whimbrels not advancing their arrival dates into Iceland? Exploring seasonal and sex-specific variation in consistency of individual timing during the annual cycle. Front Ecol Evol. 2019;7:248. https://doi.org/10.3389/fevo.2019.00248.

37. Senner NR, Verhoeven MA, Abad-Gómez JM, Alves JA, Hooijmeijer JC, Howison RA, Kentie R, Loonstra AHJ, Masero JA, Rocha A, Stager M, Piersma T. High migratory survival and highly variable migratory behavior in black-tailed godwits. Front Ecol Evol. 2019;7:96. https://doi.org/10.3389/fevo.2019.00096.

38. Verhoeven MA, Loonstra AH, Senner NR, McBride AD, Both C, Piersma T. Variation from an unknown source: large inter-individual differences in migrating black-tailed godwits. Front Ecol Evol. 2019;7:31. https://doi.org/10.3389/fevo.2019.00031.

39. Petersen MR. Reproductive ecology of emperor geese: annual and individual variation in nesting. Condor. 1992;94:383–97. https://doi.org/10.2307/1369211.

40. Hasselquist D, Montràs-Janer T, Tarka M, Hansson B. Individual consistency of long-distance migration in a songbird: significant repeatability of autumn route, stopovers and wintering sites but not in timing of migration. J Avian Biol. 2017;48:91–102. https://doi.org/10.1111/jav.01292.

41. Fraser KC, Shave A, de Greef E, Siegrist J, Garroway CJ. Individual variability in migration timing can explain long-term, population-level advances in a songbird. Front Ecol Evol. 2019;7:324. https://doi.org/10.3389/fevo.2019.00324.

42. Müller MS, Massa B, Phillips RA, Dell'omo G. Individual consistency and sex differences in migration strategies of Scopoli's shearwaters *Calonectris diomedea* despite year differences. Curr Zool. 2014;60:631–41. https://doi.org/10.1093/czoolo/60.5.631.

43. Delord K, Barbraud C, Pinaud D, Ruault S, Patrick SC, Weimerskirch H. Individual consistency in the non-breeding behavior of a long-distance migrant seabird, the Grey Petrel *Procellaria cinerea*. Mar Ornithol. 2019;47:93–103.

44. Léandri-Breton D-J, Tarroux A, Elliott KH, Legagneux P, Angelier F, Blévin P, Bråthen VS, Fauchald P, Goutte A, Jouanneau W, Tartu S, Moe B, Chaste O. Long-term tracking of an Arctic-breeding seabird indicates high fidelity to pelagic wintering areas. Mar Ecol Prog Ser. 2021;676:205–18. https://doi.org/10.3354/meps13798.

45. BirdLife International. *Sterna hirundo* (amended version of 2018 assessment). The IUCN Red List of Threatened Species 2019: e.T22694623A155537726. https://doi.org/10.2305/IUCN.UK.2019-3.RLT. T22694623A155537726.en

46. Ryslavy T, Bauer H-G, Gerlach B, Hüppop O, Stahmer J, Südbeck P, Sudfeldt C. The Red List of breeding birds of Germany. Ber Vogelschutz. 2020;57:13–112.

47. Becker PH, Ludwigs J-D. *Sterna hirundo* common tern. In: Parkin D, editors. BWP Update Vol. 6 Nos 1 and 2. New York: Oxford University Press; 2004. p. 93–139.
48. Szostek KL, Becker PH. Terns in trouble: demographic consequences of low breeding success and recruitment on a common tern population in the German Wadden Sea. J Ornithol. 2012;153:313–26. https://doi.org/10.1007/s10336-011-0745-7.

49. Zhang H, Rebke M, Becker PH, Bouwhuis S. Fitness prospects: effects of age, sex and recruitment age on reproductive value in a long-lived seabird. J Anim Ecol. 2015;84:199–207. https://doi.org/10.1111/1365-2656.12259.

50. Kürten N, Vedder O, González-Solís J, Schmaljohann H, Bouwhuis S. No detectable effect of light-level geolocators on the behaviour and fitness of a long-distance migratory seabird. J Ornithol. 2019;160:1087–95. https://doi.org/10.1007/s10336-019-01686-3.

51. Becker PH, Wendeln H. A new application for transponders in population ecology of the common tern. Condor. 1997;99:534–8. https://doi.org/10.2307/1369963.

52. Becker PH, Wink M. Influences of sex, sex composition of brood and hatching order on mass growth in common terns *Sterna hirundo*. Behav Ecol Sociobiol. 2003;54:136–46. https://doi.org/10.1007/s00265-003-0605-4.

53. Becker PH, Schmaljohann H, Riechert J, Wagenknecht G, Zajková Z, González-Solís J. Common terns on the East Atlantic flyway: temporal-spatial distribution during the non-breeding period. J Ornithol. 2016;157:927–40. https://doi.org/10.1007/s10336-016-1346-2.

54. Vedder O, Kürten N, Bouwhuis S. Intraspecific variation in and environment-dependent resource allocation to embryonic development time in common terns. Physiol Biochem Zool. 2017;90:453–60. https://doi.org/10.1086/691690.

55. Kenward RE. A manual for wildlife radio tagging. London: Academic Press; 2001.

56. R Core Team. R: A language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria. 2020. Available from: https://www.r-project.org/.

57. Rakhimberdiev E, Saveliev A. FLightR: SSM for solar geolocation. 2020. R package version 0.5.0. Available from: http://github.com/eldarrak/FLightR.

58. Rakhimberdiev E, Saveliev A, Piersma T, Karagicheva J. FLightR: an R package for reconstructing animal paths from solar geolocation loggers. Methods Ecol Evol. 2017;8:1482–7. https://doi.org/10.1111/2041-210X.12765.

59. Wotherspoon S, Sumner M, Lisovski S. BAStag: basic data processing for light based geolocation archival tags. 2016. R package version 0.1.3. Available from https://github.com/SWotherspoon/BAStag.

60. Rakhimberdiev E, Winkler DW, Bridge E, Seavy NE, Sheldon D, Piersma T, Saveliev A. A hidden Markov model for reconstructing animal paths from solar geolocation loggers using templates for light intensity. Mov Ecol. 2015;3:25. https://doi.org/10.1186/s40462-015-0062-5.

61. Redfern CP, Bevan RM. Overland movement and migration phenology in relation to breeding of Arctic Terns *Sterna paradisaea*. Ibis. 2020;162:373–80. https://doi.org/10.1111/ibl.12723.

62. Bruderer B, Boldt A. Flight characteristics of birds: I. Radar measurements of speeds. Ibis. 2001;143:178–204. https://doi.org/10.1111/j.1474-919X.2001.tb04475.x.

63. Kralj J, Martinović M, Jurinović L, Szinai P, Sütó S, Preiszner B. Geolocator study reveals east African migration route of Central European Common Terns. Avian Res. 2020;11:1–11. https://doi.org/10.1186/s40657-020-00191-z.

64. Hijmans RJ. Geosphere: spherical trigonometry. 2019. R package version 1.5–10. Available from: https://CRAN.R-project.org/package=geosphere

65. Bates D, Maechler M, Bolker B, Walker S. Fitting linear mixed-effects models using lme4. J Stat Softw. 2015;67:1–48. https://doi.org/10.18637/jss.v067.i01.

66. Kuznetsova A, Brockhoff PB, Christensen RHB. ImerTest package: tests in linear mixed effects models. J Stat Softw. 2017;82:1–26. https://doi.org/10.18637/jss.v082.i13.

67. Stoffel MA, Nakagawa S, Schielzeth H. rptR: repeatability estimation and variance decomposition by generalized linear mixed-effects models. Methods Ecol Evol. 2017;8:1639–44. https://doi.org/10.1111/2041-210X.12797.

68. Wickham H. ggplot2: elegant graphics for data analysis. New York: Springer; 2016.

69. QGIS Development Team. QGIS Geographic Information System. QGIS Association. 2021. Available from: https://www.qgis.org/.

70. Sanderson FJ, Donald PF, Pain DJ, Burfield IJ, van Bommel FP. Long-term population declines in Afro-Palearctic migrant birds. Biol Conserv. 2006;131:93–105. https://doi.org/10.1016/j.biocon.2006.02.008.

71. Mills AM. Protogyny in autumn migration: do male birds "play chicken"? Auk. 2005;122:71–81. https://doi.org/10.1093/auk/122.1.71.

72. Nisbet ICT, Szczys P, Mostello CS, Fox JW. Female common terns *Sterna hirundo* start autumn migration earlier than males. Seabird. 2011;24:103–6.

73. Nisbet ICT, Mostello CS, Veit RR, Fox JW, Afanasyev V. Migrations and winter quarters of five common terns tracked using geolocators. Water-birds. 2011;34:32–9. https://doi.org/10.1675/063.034.0104.

74. Neves VC, Nava CP, Cormons M, Bremer E, Castresana G, Lima P, Azevedo Junior SM, Phillips RA, Magalhães MC, Santos RS. Migration routes and non-breeding areas of common terns (*Sterna hirundo*) from the Azores. Emu. 2015;115:158–67. https://doi.org/10.1071/MU13112.

75. Bracey A, Lisovski S, Moore D, McKellar A, Craig E, Matteson S, Strand F, Costa J, Pekarik C, Curtis P, Niemi G, Cuthbert F. Migratory routes and wintering locations of declining inland North American Common Terns. Auk. 2018;135:385–99. https://doi.org/10.1642/AUK-17-210.1.

76. Fifield DA, Montevecchi WA, Garthe S, Robertson GJ, Kubetzki U, Rail JF. Migratory tactics and wintering areas of northern gannets (*Morus bassanus*) breeding in North America. Ornithol Monogr. 2014;79:1–63.

77. Mueller HC, Mueller NS, Berger DD, Allez G, Robichaud W, Kaspar JL. Age and sex differences in the timing of fall migration of hawks and falcons. Wilson Bull. 2000;112:214–24. https://doi.org/10.1676/0043-5643(2000)112[0214:AASDIT]2.0.CO;2.

78. Méndez V, Gill JA, Pórisson B, Vignisson SR, Gunnarsson TG, Alves JA. Paternal effects in the initiation of migratory behaviour in birds. Sci Rep. 2021;11:2782. https://doi.org/10.1038/s41598-021-81274-9.

79. Schmaljohann H, Liechti F, Bruderer B. First records of Lesser Black-backed Gulls (*Larus fuscus*) crossing the Sahara non-stop. J Avian Biol. 2008;39:233–7. https://doi.org/10.1111/j.2007.0908-8857.04174.x.

80. Lamb JS, Newstead DJ, Kozkur LM, Ballard BM, Green MC, Jodice PG. A bridge between oceans: overland migration of marine birds in a wind energy corridor. J Avian Biol. 2018;49:e01474. https://doi.org/10.1111/jav.01474.

81. Kok EM, Tibbitts TL, Douglas DC, Howey PW, Dekinga A, Gnep B, Piersma T. A red knot as a black swan: how a single bird shows navigational abilities during repeat crossings of the Greenland Icecap. J Avian Biol. 2020;51:e02464. https://doi.org/10.1111/jav.02464.

82. Piersma T, Kok EM, Hassell CJ, Peng HB, Verkuil YI, Lei G, Karagicheva J, Rakhimberdiev E, Howey PW, Tibbitts L, Chan YC. When a typical jumper skips: itineraries and staging habitats used by Red Knots (*Calidris canutus piersmai*) migrating between northwest Australia and the New Siberian Islands. Ibis. 2021;163:1235–51. https://doi.org/10.1111/ibi.12964.

83. Grecian WJ, Witt MJ, Attrill MJ, Bearhop S, Becker PH, Egevang C, Furness RW, Godley BJ, González-Solís J, Grémillet D, Kopp M, Lescroël A, Matthiopoulos J, Patrick SC, Peter H-P, Phillips RA, Stenhouse IJ, Votier SC. Seabird diversity hotspot linked to ocean productivity in the Canary Current Large Marine Ecosystem. Biol Lett. 2016;12:20160024. https://doi.org/10.1098/rsbl.2016.0024.

84. Crawford RJM. Food and population variability in five regions supporting large stocks of anchovy, sardine and horse mackerel. S Afr J Mar Sci. 1987;5:735–57. https://doi.org/10.2989/025776187784522243.

85. Binet D. Climate and pelagic fisheries in the Canary and Guinea currents 1964–1993: the role of trade winds and the southern oscillation. Oceanol Acta. 1997;20:177–90.

86. Wynn RB, Knefelkamp B. Seabird distribution and oceanic upwelling off northwest Africa. Br Birds. 2004;97:323–35.

87. Camphuysen CJ, van der Meer J. Wintering seabirds in West Africa: foraging hotspots off Western Sahara and Mauritania driven by upwelling and fisheries. Afr J Mar Sci. 2005;27:427–37. https://doi.org/10.2989/18142320590540101.

88. Redfern CP, Kinchin-Smith D, Newton S, Morrison P, Bolton M, Piec D. Upwelling systems in the migration ecology of Roseate Terns (*Sterna hirundo*).
dougallii) breeding in northwest Europe. Ibis. 2021;163:549–65. https://doi.org/10.1111/ibi.12915.

89. Block BA, Jonsen ID, Jorgensen SJ, Winship AJ, Shaffer SA, Bograd SJ, Hazen EL, Foley DG, Breed GA, Harrison A-L, Ganong JE, Swithenbank A, Castleton M, Dewar H, Mate BR, Shillinger GL, Schaefer KM, Benson SR, Weise MJ, Henry RW, Costa DP. Tracking apex marine predator movements in a dynamic ocean. Nature. 2011;475:86–90. https://doi.org/10.1038/nature10082.

90. Hoegh-Guldberg O, Bruno JF. The impact of climate change on the world's marine ecosystems. Science. 2010;328:1523–8. https://doi.org/10.1126/science.1189930.

91. Blanchard JL, Jennings S, Holmes R, Harle J, Merino G, Allen JI, Holt J, Dulvy NK, Barange M. Potential consequences of climate change for primary production and fish production in large marine ecosystems. Philos Trans R Soc B. 2012;367:2979–89. https://doi.org/10.1098/rstb.2012.0231.

92. Sydeman WJ, Poloczanska E, Reed TE, Thompson SA. Climate change and marine vertebrates. Science. 2015;350:772–7. https://doi.org/10.1126/science.aac9874.

93. Crawford RJ, Makhado AB, Whittington PA, Randall RM, Oosthuizen WH, Waller LJ. A changing distribution of seabirds in South Africa—the possible impact of climate and its consequences. Front Ecol Evol. 2015;3:10. https://doi.org/10.3389/fevo.2015.00010.

94. Sydeman WJ, Schoeman DS, Thompson SA, Hoover BA, García-Reyes M, Daunt F, Agnew P, Anker-Nilssen T, Barbraud C, Barrett R, Becker PH, Bell E, Boersma PD, Bouwhuis S, Cannell B, Crawford RJM, Dann P, Delord K, Elliott G, Erikstad KE, Flint E, Furness RW, Harris MP, Hatch S, Hilwig K, Hinke JT, Jahncke J, Mills JA, Reiertsen TK, Renner H, Sherley RB, Surman C, Taylor G, Thayer JA, Trathan PN, Velarde E, Walker K, Wanless S, Warzybok P, Watanuki Y. Hemispheric asymmetry in ocean change and the productivity of marine ecosystem sentinels. Science. 2021;372:980–3. https://doi.org/10.1126/science.abf1772.

95. Nilsson C, Klaassen RH, Alerstam T. Differences in speed and duration of bird migration between spring and autumn. Am Nat. 2013;181:837–45. https://doi.org/10.1086/670335.

96. Schmaljohann H. Proximate mechanisms affecting seasonal differences in migration speed of avian species. Sci Rep. 2018;8:4106. https://doi.org/10.1038/s41598-018-22421-7.

97. Moiron M, Araya-Ajoy YG, Teplitsky C, Bouwhuis S, Charmantier A. Understanding the social dynamics of breeding phenology: indirect genetic effects and assortative mating in a long distance migrant. Am Nat. 2020;196:566–76. https://doi.org/10.1086/711045.

98. Dobson FS, Becker PH, Arnaud CM, Bouwhuis S, Charmantier A. Plasticity results in delayed breeding in a long-distance migrant seabird. Ecol Evol. 2017;7:3100–9. https://doi.org/10.1002/ece3.2777.

99. Liechti F. Birds: Blowin' by the wind? J Ornithol. 2006;147:202–11. https://doi.org/10.1007/s10336-006-0614-7.

100. Barrett RT. Upwind or downwind: the spring arrival of Arctic Terns *Sterna paradisaea* at Troms, north Norway. Ring Mig. 2016;31:23–9. https://doi.org/10.1080/03078698.2016.1190610.

101. Haest B, Hüppop O, Bairlein F. The influence of weather on avian spring migration phenology: what, where and when? Glob Chang Biol. 2018;24:5769–88. https://doi.org/10.1111/gcb.14450.

102. Haest B, Hüppop O, Pol M, Bairlein F. Autumn bird migration phenology: a potpourri of wind, precipitation and temperature effects. Glob Chang Biol. 2019;25:4064–80. https://doi.org/10.1111/gcb.14746.

103. Cresswell W. Migratory connectivity of Palaearctic-African migratory birds and their responses to environmental change: the serial residency hypothesis. Ibis. 2014;156:493–510. https://doi.org/10.1111/ibi.12168.

104. Wynn J, Guilford T, Padget O, Perrins CM, Mckee N, Gillies N, Tyson C, Dean B, Kirk H, Fayet AL. Early-life development of contrasting outbound and return migration routes in a long-lived seabird. Ibis. 2021. https://doi.org/10.1111/ibi.13030.

105. Smouse PE, Focardi S, Moorcroft PR, Kie JG, Forester JD, Morales JM. Stochastic modelling of animal movement. Philos Trans R Soc B. 2010;365:2201–11. https://doi.org/10.1098/rstb.2010.0078.

106. Sergio F, Tanferna A, De Stephanis R, Jiménez LL, Blas J, Tavecchia G, Preatoni D, Hiraldo F. Individual improvements and selective mortality shape lifelong migratory performance. Nature. 2014;515:410–3. https://doi.org/10.1038/nature13696.

107. Campioni L, Dias MP, Granadeiro JP, Catry P. An ontogenetic perspective on migratory strategy of a long-lived pelagic seabird: timings and destinations change progressively during maturation. J Anim Ecol. 2020;89:29–43. https://doi.org/10.1111/1365-2656.13044.

108. Robertson GJ, Cooke F. Winter philopatry in migratory waterfowl. Auk. 1999;116:20–34. https://doi.org/10.2307/4089450.

109. BirdLife International and Durham University. Species climate change impacts factsheet: hirundo. 2021. http://www.birdlife.org. Accessed 10 Apr 2021.

110. Saino N, Szép T, Ambrosini R, Romano M, Møller AP. Ecological conditions during winter affect sexual selection and breeding in a migratory bird. Proc R Soc B. 2004;271:681–6. https://doi.org/10.1098/rspb.2003.2656.

111. Rockwell SM, Boцetti CI, Marra PP. Carry-over effects of winter climate on spring arrival date and reproductive success in an endangered migratory bird, Kirtland's Warbler (Setophaga kirtlandii). Auk. 2012;129:744–52. https://doi.org/10.1525/auk.2012.12003.

112. Wernham CV, Toms MP, Marchant JH, Clark JA, Siriwardena GM, Baillie SR. The migration atlas: movements of the birds of Britain and Ireland. London: T & AD Poyser; 2002.

113. Bairlein F, Dierschke J, Dierschke V, Salewski V, Geiter O, Hüppop K, Köppen U, Fielder W. Atlas des Vogelzugs. Wiebelsheim: AULA; 2014.

114. Zhang SD, Ma Z, Choi CY, Peng HB, Bai QQ, Liu WL, Tan K, Melville DS, He P, Chan Y-C, van Gils JA, Piersma T. Persistent use of a shore-bird staging site in the Yellow Sea despite severe declines in food resources implies a lack of alternatives. Bird Conserv Int. 2018;28:534–48. https://doi.org/10.1017/S0959270917000430.

115. Dias MP, Granadeiro JP, Catry P. Individual variability in the migratory path and stopovers of a long-distance pelagic migrant. Anim Behav. 2013;86:359–64. https://doi.org/10.1016/j.anbehav.2013.05.026.

116. Loonstra AJ, Verhoeven MA, Zbyryt A, Schaaf E, Both C, Piersma T. Individual Black-tailed Godwits do not stick to single routes: a hypothesis on how low population densities might decrease social conformity. Ardea. 2020;107:251–61. https://doi.org/10.5253/arde.v107i3.a11.

117. López-López P, García-Ripollés C, Urios V. Individual repeatability in timing and spatial flexibility of migration routes of trans-Saharan migratory raptors. Curr Zool. 2014;60:642–52. https://doi.org/10.1093/czoolo/60.5.642.

118. Stanley CQ, MacPherson M, Fraser KC, McKinnon EA, Stutchbury BJ. Repeat tracking of individual songbirds reveals consistent migration timing but flexibility in route. PLoS ONE. 2012;7:e40688. https://doi.org/10.1371/journal.pone.0040688.

119. Mellone U, López-López P, Limiñana R, Urios V. Weather conditions promote route flexibility during open ocean crossing in a long-distance migratory raptor. Int J Biometeorol. 2011;55:463–8. https://doi.org/10.1007/s00484-010-0368-3.

120. González-Solís J, Felicísimo A, Fox JW, Afanasyev V, Kolbeinsson Y. Muñoz J. Influence of sea surface winds on shearwater migration detours. Mar Ecol Prog Ser. 2009;391:221–30. https://doi.org/10.3354/meps08128.

121. Robson D, Barriocanal C. The influence of environmental conditions on the body mass of Barn Swallows (*Hirundo rustica*) during spring migration. J Ornithol. 2008;149:473–8. https://doi.org/10.1007/s10336-008-0297-7.

122. Szostek KL, Bouwhuis S, Becker PH. Are arrival date and body mass after spring migration influenced by large-scale environmental factors in a migratory seabird. Front Ecol Evol. 2015;3:42. https://doi.org/10.3389/fevo.2015.00042.

123. Åkesson S, Ilieva M, Karagicheva J, Rakhimberdiev E, Tomotani B, Helm B. Timing avian long-distance migration: from internal clock mechanisms to global flights. Phil Trans R Soc B. 2017;372:20160252. https://doi.org/10.1098/rstb.2016.0252.

124. Gill JA, Alves JA, Sutherland WJ, Appleton GF, Potts PM, Gunnarsson TG. Why is timing of bird migration advancing when individuals are not? Proc R Soc B. 2014;281:20132161. https://doi.org/10.1098/rspb.2013.2161.

125. Youngflesh C, Jenouvrier S, Hinke JT, DuBois L, St. Leger J, Trivelpiece WZ, Trivelpiece SG, Lynch HJ. Rethinking "normal": the role of stochasticity in the phenology of a synchronously breeding seabird. J Anim Ecol. 2018;87:682–90. https://doi.org/10.1111/1365-2656.12790.
126. Vedder O, Pen I, Bouwhuis S. How fitness consequences of early-life conditions vary with age in a long-lived seabird: a Bayesian multivariate analysis of age-specific reproductive values. J Anim Ecol. 2021;90:1505–14. https://doi.org/10.1111/1365-2656.13471.

127. Vedder O, Zhang H, Bouwhuis S. Early mortality saves energy: estimating the energetic cost of excess offspring in a seabird. Proc R Soc B. 2017;284:20162724. https://doi.org/10.1098/rspb.2016.2724.

128. Vedder O, Zhang H, Dänhardt A, Bouwhuis S. Age-specific offspring mortality economically tracks food abundance in a piscivorous seabird. Am Nat. 2019;193:88–597. https://doi.org/10.1086/702304.

129. Bogdanova MI, Daunt F, Newell M, Phillips RA, Harris MP, Wanless S. Seasonal interactions in the black-legged kittiwake, *Rissa tridactyla*: links between breeding performance and winter distribution. Proc R Soc B. 2011;278:2412–8. https://doi.org/10.1098/rspb.2010.2601.

130. Bogdanova MI, Butler A, Wanless S, Moe B, Anker-Nilssen T, Frederiksen M, Boulinier T, Chivers LS, Christensen-Dalsgaard S, Descamps S, Harris MP, Newell M, Olsen B, Phillips RA, Shaw D, Steen H, Strøm H, Thórarinsson TL, Daunt F. Multi-colony tracking reveals spatio-temporal variation in carry-over effects between breeding success and winter movements in a pelagic seabird. Mar Ecol Prog Ser. 2017;578:167–81. https://doi.org/10.3354/meps12096.

131. Helm B, Piersma T, van der Jeugd H. Sociable schedules: interplay between avian seasonal and social behaviour. Anim Behav. 2006;72:245–62. https://doi.org/10.1016/j.anbehav.2005.12.007.

**Publisher's Note**

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Ready to submit your research? Choose BMC and benefit from:

* fast, convenient online submission

* thorough peer review by experienced researchers in your field

* rapid publication on acceptance

* support for research data, including large and complex data types

* gold Open Access which fosters wider collaboration and increased citations

* maximum visibility for your research: over 100M website views per year

At BMC, research is always in progress.

Learn more biomedcentral.com/submissions

<figure><img src="image_6.png" /></figure>