Individual consistency in migration strategies of a tropical seabird, the Round Island petrel
Movement Ecology volume 10, Article number: 13 (2022)
In migratory species, the extent of within- and between-individual variation in migratory strategies can influence potential rates and directions of responses to environmental changes. Quantifying this variation requires tracking of many individuals on repeated migratory journeys. At temperate and higher latitudes, low levels of within-individual variation in migratory behaviours are common and may reflect repeated use of predictable resources in these seasonally-structured environments. However, variation in migratory behaviours in the tropics, where seasonal predictability of food resources can be weaker, remains largely unknown.
Round Island petrels (Pterodroma sp.) are tropical, pelagic seabirds that breed all year round and perform long-distance migrations. Using multi-year geolocator tracking data from 62 individuals between 2009 and 2018, we quantify levels of within- and between-individual variation in non-breeding distributions and timings.
We found striking levels of between-individual variation in at-sea movements and timings, with non-breeding migrations to different areas occurring across much of the Indian Ocean and throughout the whole year. Despite this, repeat-tracking of individual petrels revealed remarkably high levels of spatial and temporal consistency in within-individual migratory behaviour, particularly for petrels that departed at similar times in different years and for those departing in the austral summer. However, while the same areas were used by individuals in different years, they were not necessarily used at the same times during the non-breeding period.
Even in tropical systems with huge ranges of migratory routes and timings, our results suggest benefits of consistency in individual migratory behaviours. Identifying the factors that drive and maintain between-individual variation in migratory behaviour, and the consequences for breeding success and survival, will be key to understanding the consequences of environmental change across migratory ranges.
Rapid changes in environmental conditions are affecting ecosystems, communities and species worldwide . For migratory species that are dependent upon the availability of habitats, resources and conditions in multiple locations across migratory ranges, differing rates of change in these areas can greatly increase the potential for deleterious impacts at some point in the annual cycle . Consequently, species lacking variability and flexibility in their migratory traits which might facilitate responses to changing environmental conditions could be at a disadvantage . For example, among European breeding birds, species which have shown little or no phenological change in recent decades tend to also be those that are currently in population decline . In order to identify potential constraints on migratory species’ responses to environmental change, we therefore need to understand variability and flexibility in migratory behaviour.
Recent advances in remote tracking technology have facilitated numerous studies following the movements of birds across multiple annual cycles . These repeated measures of individuals allow investigation of spatiotemporal consistency (or, conversely, flexibility) in migration strategies within and among individuals (e.g., [8, 35]). As a result, there are now many avian studies which have investigated individual consistency in migratory timings (reviewed by Franklin et al., in review, ), and non-breeding locations [13, 42] and, more recently, year-to-year fidelity in migratory routes  and stopovers [8, 22]. Most studies have revealed high individual time-  and site-fidelity across years , with variation existing between individuals. However, some studies have also recorded individuals changing migratory behaviours between years [14, 35]. The occurrence and persistence of individual consistency in migratory behaviour may have emerged as a result of spatial and temporal predictability of resource availability , with familiarity of conditions at known locations and times being more beneficial than trying to locate optimal conditions at any given time [8, 60]. Predictable resource distributions might therefore be expected to result in low within-individual variation in space and time across years, and resource landscapes that are heterogeneous and predictable might be expected to result in between-individual variation and within-individual consistency . By contrast, migratory species in environments in which resources are less predictable in space and time might be expected to show higher levels of within-individual flexibility. As most tracking studies have been conducted on species breeding at temperate and higher latitudes, levels of migratory consistency in less seasonal and unpredictable environments, such as tropical systems, remains unclear  (Franklin et al., in review).
Seabirds wintering in temperate and polar regions often associate with physical oceanographic features, such as oceanic fronts, shelf and ice edges or upwellings. These features, along with seasonal temperature and salinity gradients, tend to lead to temporally and spatially predictable prey aggregations . Individual birds in these systems tend to have predictable migrations to one or more of these high-productivity ocean areas and show high levels of migration fidelity between years (e.g., ). By contrast, large areas of tropical oceans are often considered low in productivity and prey abundance, and have less marked seasonal variation in temperature, making prey aggregations unpredictable , but see . Seabirds foraging at lower latitudes also often rely on subsurface marine predators, such as dolphin and tuna, that drive prey to the surface [3, 24, 49], and these events are likely to be less predictable than static oceanographic variables related to marine productivity. Non-breeding migrations of tropical species may therefore be less predictable both within- and between- individuals. However, only a very restricted number of studies have tracked individual tropical seabirds on multiple migrations [25, 44], and typically in such small sample sizes that quantifying within-individual variability has not yet been possible. So far for tropical seabirds, studies have mainly focused on variation in behaviours at breeding grounds, when individuals are constrained to a central location, and have revealed no consistent patterns in foraging site-fidelity [11, 40]. Consequently, it remains unclear whether consistent individual migratory behaviours occur in less predictable, tropical environments.
We addressed these questions using a tropical, pelagic seabird, the Round Island petrel, as a model species. This population of gadfly petrels (genus Pterodroma) breeds all-year round on Round Island, Mauritius, in the western Indian Ocean, and has been the focus of a long-term geolocator tracking project . The mid-ocean location of Round Island means that petrel migrations could potentially occur in any compass direction and for a huge range of distances. Here we estimate how repeatable individual Round Island petrels are in their migratory timings (arrival to and departure from Round Island), and migration duration (time away from Round Island). We then use the earth mover’s distance (EMD), an algorithm originally developed for image comparison  and subsequently adapted to quantify similarity between spatial distributions , to quantify a) spatial consistency of petrel migrations across the entire non-breeding distribution, and b) whether individuals consistently occur in the same locations at the same stages (~ monthly) of their migration schedules.
Study site and species
This study was carried out at Round Island Nature Reserve (19.85° S, 57.78° E), a 219-ha island situated 22.5 km off the North coast of mainland Mauritius. The climate of Mauritius and the surrounding ocean is strongly seasonal and can be divided into two broad seasons: the austral winter (hereafter ‘winter’), and the austral summer (hereafter ‘summer’). The former typically runs from May to September and the latter from October to April, which reflects the influence of the monsoon circulation of the Indian Ocean . Round Island is the only confirmed colony in the Indian Ocean of an unusual population of Pterodroma petrels, known as the Round Island petrel. Genetic evidence has shown this population to comprise at least three species of Pterodroma petrel (Trindade petrel, P. arminjoniana, Kermadec petrel, P. neglecta and Herald petrel, P. heraldica), which breed and extensively hybridise on the island [4, 6, 7, 52, 53]. Round Island petrels nest on the ground, typically under rock ledges or among piles of boulders. Eggs and chicks can be present in any month of the year, although there is a peak in egg-laying in August-October . As petrel breeding activity is typically lowest on the island in May each year, we selected 1 June as the start of the petrel breeding calendar and numbered days sequentially from this origin .
Petrel surveys have been undertaken monthly since 2001, and involve regular visits to known nesting areas, ringing of adults and chicks (with South African Bird Ringing Unit numbered rings), and their subsequent recapture. Between 2009 and 2016, 421 light-level geolocators (GLS) were deployed on adult petrels. GLS were attached to the tarsus via a 1 or 0.75 mm thick Salbex (an industrial grade PVC; Sallu Plastics, Redditch, UK) colour-ring. Between 2009 and 2012, MK15 British Antarctic Survey geolocators (Cambridge, UK) were deployed, and during 2014 and 2016, Intigeo C250 and Intigeo C330 (Migrate Technology, Cambridge, UK) were deployed. The total device weight (including plastic ring) across all three tag types amounted to 3.6—3.9 g, which represents approximately 1.0% of the mean body mass of adult petrels (374 g). For the first three deployment periods, the tagged petrels were caught during targeted searches, whereas later GLS deployments occurred during the standard monthly petrel breeding surveys; these birds were predominantly resting on the island and not directly observed in a breeding attempt. Loggers were opportunistically recovered a minimum of one year later, during breeding surveys or during occasional specific searches. All loggers underwent a 3- to 5-day calibration period at a known location (Round Island, or mainland Mauritius (20.25° S, 57.44° E)) pre- and post- deployment. Details of the numbers of tags deployed and recovered until the end of December 2019 are provided in Additional File 1: Table S1.
Geolocation data processing
At-sea locations for each individual were estimated from raw light-level data using the threshold method of estimating positions based on twilight events (i.e., sunrise and sunset transitions; ). Twilight events were defined using the preprocessLight function in the R package TwGeos  using a light intensity threshold of 4 and 1, for MK15 and Intigeo tags, respectively. The corresponding zenith angle was defined separately for each tag from the sunrise and sunset times recorded during the pre-deployment calibration period (range of 94.0° to 96.9° for MK15 tags and 96.8° to 99.0° for Intigeo tags). A Bayesian framework was used to refine the initial positions estimated from the threshold method and to derive uncertainty estimates. The R package SGAT  uses Markov Chain Monte Carlo (MCMC) simulations allowing the incorporation of: (1) a spatial probability mask, (2) sea-surface temperatures (SST) recorded by the tag in relation to global remote sensed SST maps, (3) prior definition of the error distribution of twilight events (twilight model) and (4) a flight speed distribution (behavioural model), to refine location estimates. The twilight model should reflect the expected error in detecting the real time of sunrise and sunset. Since the petrels spend a substantial amount of time sitting when at the breeding site (obscuring the light sensor) and can travel many hundreds of kilometres when still associated with the island, we could not use twilight times from a known location (i.e., Round Island) to parameterise the twilight model. We therefore used a rather conservative prior (log-normal distribution: meanlog = 2.2, sdlog = 1.0) describing a large variation in the discrepancy between the real and recorded twilight events. The spatial probability mask was constructed on the premise that birds only use marine environments when away from the colony (probability of 0 for positions on land). The probability at sea was further refined for each individual using remotely sensed sea-surface temperatures (weekly means on a 1.0° × 1.0° resolution downloaded from NOAA’s Optimum Interpolation Sea Surface Temperature V2 dataset: https://www.esrl.noaa.gov/psd/data/gridded/data.noaa.oisst.v2.html) and the SST values recorded by the loggers. The potential flight speeds were modelled following a gamma distribution (shape = 2.0, rate = 0.1). For each individual, we used these parameters and started by drawing an initial 2000 samples for burn-in and tuning of the proposal distribution using a modified model with relaxed assumptions. Then a further 300 samples were drawn three times to evaluate chain convergence before drawing another 4 MCMC chains of 3000 samples each to describe the posterior distribution. Tracks were summarised to produce median tracks and 95% credible intervals.
To identify the start and end dates of migration from GLS data, each day of GLS tracking was classified as one of two behavioural states (ashore on Round Island or at-sea) using a hidden Markov model (HMM). HMMs are a type of state‐space model, which decompose observed time‐series data into an observed sequence of discrete behavioural states. Behaviours were classified as either ashore or at-sea from the proportion of light interference during core daylight hours, the sum of daily wet records, and the distance from Round Island, using the R package depmixS4 . These variables were chosen as when petrels are frequenting the colony, we expect; (1) a high degree of shading of the tag as individuals are sitting, (2) longer periods without immersion, and (3) petrels to remain close to the colony. Distance from colony (km) was derived from the SGAT processed median location estimates and calculated using the distGeo function in the R package geosphere , whereas light and immersion data were from the respective raw GLS files. All GLS models sampled light every minute and logged the maximum light recorded at 10-min intervals. To calculate light interference, we used the twilight times defined from the preprocessLight function to select out light recordings during daylight hours only. In contrast to the MK15 GLS, which record only low light levels, Intigeo GLS record the entire light range. Therefore, light samples with values of less than 64 or less than 100 were classed as interference for MK15 and Intigeo, respectively. Due to the different number of light samples that fall within each daylight period, this variable was calculated as the proportion of light samples with interference for each calendar day. GLS tested for saltwater immersion every 3 or every 30 s and stored the sum of positive tests at 10-min intervals, resulting in values between 0 (entirely dry), and 200 or 20 (entirely wet), for the MK15 and Intigeo GLS, respectively. The sum of these 144 values (number of 10-min periods in 24-h) were calculated for each calendar day. A Gaussian distribution was used to describe both distance from Round Island and immersion data, and a binomial distribution to describe light interference. Migration periods were defined as a sequence of consecutive days that were assigned to the same behaviour (at-sea) by the HMM for a period of at least three months, before switching to the other behavioural state (ashore). The duration of the migratory period was calculated using these dates.
To investigate consistency in migratory locations within- and between- individuals, we used the median locations identified as the migration period from the HMM in combination with the earth mover’s distance (EMD; ) to create a matrix of space-use similarity estimates. In contrast to spatial overlap indices (where distributions with no spatial overlap have the same similarity value irrespective of their distance from one another), the EMD integrates a measure of spatial proximity in the similarity between different space-use distributions, by calculating the effort it takes to transform, for example, one migration track into another. EMD ‘effort’ values are therefore on a continuous scale, starting at zero for two identical migrations and increase with increasing dissimilarity. EMD was calculated for all petrel and year combinations, using Haversine distance with the ‘emd’ function in the R package move , using geographical coordinates directly as inputs.
In order to show how EMD relates to a widely used method for quantifying home range overlap, we also calculated relative overlap of petrel migrations using Bhattacharyya's affinity (BA; ). BA estimates range between 0 (no overlap) to 1 (identical distributions) and therefore do not quantify how dissimilar migratory distributions are if they do not overlap. For each individual and year, migratory locations were used to generate kernel utilisation distributions (UDs) using the R package adehabitatHR . A fixed smoothing parameter (h) of 200 km was used to account for precision error around location estimates . BA was calculated between all possible paired combinations of petrel UDs using the ‘kerneloverlap’ function in the R package adehabitatHR. This analysis was carried out in a Lambert azimuthal equal-area projection centred on the centroid of all positions. By comparing with BA, we were also able to investigate the impact of any scale-dependence on our analyses of EMD. As EMD is an absolute metric, individuals covering large migratory distances have greater capacity to differ and are thus more likely to generate large EMD values.
To investigate the spatiotemporal similarity of the same individual’s migrations in different years, two different approaches were used to define time: each petrel’s migrations was split into (1) 30-day periods irrespective of the migration start date (Additional File 1: Fig. S1a), meaning that a short period of time may be excluded for within-individual migrations of different durations; and into (2) six stages of equal duration (average of 29.4 days ± 4.7 SD; Additional File 1: Fig. S1b), meaning that the corresponding within-individual stage between migrations may be of a different length but in total spans the whole migratory period. EMD values were calculated for within-individual migrations for the same 30-day period/stage only (i.e., comparing the first 30-day period/stage of the first migration to the first 30-day period/stage of the second migration etc.), using the same method as described above.
Repeatability (R) of (a) arrival to, (b) departure from Round Island and (c) migration duration were each estimated in a mixed-effects model framework, using 1000 bootstrap iterations to estimate confidence intervals (CIs), with the R package rptR . Arrival and departure dates were recorded as days from 1 June; individuals with consecutive migrations spanning 1 June, (e.g., arriving on 30 May (day 364) in year one, and 2 June (day 2) in year two) had the first date converted to a negative day value to overcome the circular nature of this variable (the inclusion or exclusion of negative values had no significant impact on the estimated repeatability values). All three models were fitted with a Gaussian error family and included ‘individual identity’ as a random effect.
To investigate whether petrels arriving/departing in the different seasons differ in their levels of variability in migratory timings, the number of days between the earliest and latest date for each individual (i.e., within-individual variation) was calculated separately for arrival and departure dates, for petrels which consistently arrived or departed in either the summer or winter only. These values were included as the response variable in generalised linear models (GLMs) with gamma error family and log link function for arrival and departure dates separately, with season as a fixed effect.
To quantify spatial repeatability in individual petrel migrations and whether this varied with differences in departure timing, the EMD values comparing whole petrel migrations were included as the response variable in a GLM with the binary factor ‘same individual’ (0 as ‘no’ and 1 as ‘yes’), and the difference in days between the compared individuals’ departure dates from Round Island as fixed factors. An interaction between ‘same individual’ and difference in departure dates was also included. A significantly negative coefficient for the ‘same individual’ classification would indicate higher overlap (i.e., greater spatial similarity) within-individuals than between-individuals. The EMD values were continuous, non-negative and right-skewed; therefore, a gamma error family and identity link function were used. To examine if scale-dependence in EMD influenced these findings, we ran a second GLM with binomial error family and log link function with the same fixed factors, but this time with the BA values as the response variable.
To investigate within-individual consistency in distribution at different stages during the migration period, GLMs with gamma error family and identity link function were run with EMD values for each (a) stage or (b) period as the response variable, and difference in days between departure dates, and stage/period of migration (categorical variable including a level for whole migration) as fixed factors. An interaction between the two fixed factors was also included, to explore whether similar departure timing between years results in more consistent use of the same locations at the same time in each year. For all models, non-significant (p < 0.05) terms were sequentially removed using backwards stepwise deletion and significance of terms was determined using ANOVA. Pairwise comparisons among levels were calculated based on estimated marginal means and adjusted using post-hoc Tukey correction using the R package emmeans . All analyses were conducted using R version 3.6.2 .
A total of 337 retrieved GLS loggers provided data on 267 complete migration tracks of 198 Round Island petrel individuals between 2009 and 2019 (Additional File 1: Table S1). This includes repeated tracking of the same individuals providing complete migration tracks during two (n = 57), three (n = 3) or four (n = 2) migratory periods. The five individuals with three or four migratory periods were a result of two separate non-consecutive GLS deployments, rather than birds evading capture for multiple years following deployment. These tracks provided arrival dates over a period of 2 to 4 years, and migration durations for up to 4 years, for 62 petrels. Due to logger failure during deployment, partial data recovery was possible for a number of GLS devices, providing a total of 76 petrels with between 2 and 5 departure dates from Round Island (Table 1).
Timing of migration
GLS tracking of sampled individuals covered departure and arrival dates of petrel migrations spanning all months. Repeated tracking of individual petrels indicated a much higher degree of consistency within- than between-individuals in all three migratory timings (arrival to, and departure from Round Island, and migration duration), with individual arrival dates to Round Island being the most repeatable (Fig. 1; Table 1). This repeatability was consistent across all years of the study period (Fig. 1). Although significantly repeatable, there is some variation within-individuals (median difference between latest and earliest departure date for each individual = 56 days (range 1–240 days) and for arrival = 47 days (range 2–220 days)), and individuals with a larger departure date range were also more likely to have a larger arrival date range (Pearson's correlation: r = 0.64, t = 6.38, df = 60, p < 0.0001). This within-individual repeatability in arrival and departure also differed between the two seasons which petrels arrived and departed, with the 17 birds departing in the winter (i.e. off-peak breeding period) having significantly higher within-individual variation than the 39 departing in the summer (i.e. peak breeding period, GLM: adjusted R2 = 0.25, β = 0.79, SE = 0.23, t = 6.55, p < 0.001), and the 10 birds arriving in the summer, having significantly higher within-individual variation than the 37 arriving in the winter (GLM: adjusted R2 = 0.48, β = 56.21, SE = 8.58, t = 6.55, p < 0.0001).
Migration durations of the 131 tracks (from the 62 repeatedly tracked petrels) lasted on average 175 days ± 28 SE (range 104–256 days). Despite a lower repeatability value for migration duration (Table 1), individuals showed relatively small mean individual ranges (difference between the largest and smallest record for each individual, in days; 22.8 ± 19.6 SD), highlighting that this low R value may reflect relatively low variance in migration duration among individuals. The variance components for each of the repeatability estimates can be found in Additional File 1: Table S8.
Spatial consistency in migratory journeys
Round Island petrels showed striking levels of between-individual variation in migratory journeys, with individuals undertaking non-breeding migrations to different areas across much of the Indian Ocean north of ~ 35°S (Fig. 2; Additional File 1: Fig. S2). Most birds migrated north to the Somali Basin and further into the Arabian Sea, or east across the central Indian Ocean between 10 and 20°S. Comparatively few travelled into the Bay of Bengal or Western Australian Basin, but all petrels largely avoided nearshore/shelf waters (Fig. 2; Additional File 1: Fig. S2). Despite this large between-individual spatial variation, repeated tracking of individual petrels for up to 4 years indicates remarkably high levels of spatial consistency (Figs. 2, 3), with within-individual petrel migrations being significantly more similar than between-individual migrations (Figs. 3, 4a; Table 2; Additional File 1: Table S2). EMD values for within-individual petrel migration comparisons ranged from 156.0 to 1618.7 compared to a range of 226.4 to 5419.0 for between-individual migrations (Figs. 3, 4a). The highest EMD value of 5419.0 indicates, in this case, the difference between individuals travelling north to the Arabian Sea versus east to the Western Australian Basin (Fig. 3f). The relative overlap of this pair of migrations when calculated using BA revealed a very low overlap of 0.06 (Additional File 1: Fig. S3a). However, as BA does not take into account the spatial distance between migrations, the same BA value can be observed across a range of EMD values, which represent paired migrations that are relatively close together in space (EMD = 2090.8; Additional File 1: Fig. S3c) or far apart (EMD = 5419.0; Additional File 1: Fig. S3a).
Within-individual spatial consistency was also significantly higher for individuals that departed Round Island at more similar times in each year (Fig. 4b; Table 2; Additional File 1: Table S2). However, this effect is very small in comparison to the large amount of variation in EMD values across the range of differences in departure timing (R2 = 0.020; Fig. 4). There was no significant interaction between ‘same individual’ and difference in departure dates meaning that the greater similarity of migratory ranges for birds departing at similar times applies both within- and between- individuals (Additional File 1: Table S2). Overall, comparison with the BA analysis showed both methods yielded broadly consistent results suggesting any scale-dependence in EMD values does not vary systematically in relation to difference in departure date, and thus does not affect our conclusions (Additional File 1: Tables S2, S3).
Despite high within-individual consistency in entire migratory journeys, individual Round Island petrels do not always use the same at-sea locations at the same stage or period within their migratory journeys. EMD values calculated for each stage/period of each petrel’s migrations were slightly, but significantly, higher than the EMD values for the whole migration (Fig. 5; Table 2; Additional File 1: Tables S2, S4–S6). These differences in the EMD values across stages/periods are not a result of individuals departing later or earlier in one migration compared to the other, as there was no significant interaction between stage/period and the difference in departure dates between years. Nonetheless, the overall pattern of petrels that departed Round Island at more similar times of year having more similar migrations still exists (Table 2).
Tropical, migratory Round Island petrels have striking levels of between-individual variation in at-sea movement patterns, with individuals undertaking migrations across much of the Indian Ocean and throughout the whole year. However, repeat tracking of individuals across the non-breeding period revealed very little within-individual variation, with migratory journeys being remarkably consistent in both space and time.
The areas frequented by Round Island petrels during migration cover much of the Indian Ocean north of ~ 35°S. Compared with most other studies tracking seabirds breeding in the western Indian Ocean, albeit in smaller numbers, the range of areas used by Round Island petrels during the non-breeding period is particularly large [30, 44]. GLS tracking of sooty terns (Onychoprion fuscatus) from Bird Island, Seychelles, also revealed use of a range of different non-breeding areas, three of which (the Bay of Bengal, northeast to an area straddling the Chagos-Laccadive plateau, southeast to an area on each side of the 90 East Ridge; ) are also used by Round Island petrels. However, the majority of Indian Ocean seabirds for which migratory journeys have been tracked, including the closely related Barau’s petrel (Pterodroma baraui), which nests on nearby Réunion Island, tend to show consistent eastward migrations to specific areas of the central and eastern Indian Ocean . Systems with large non-breeding distributions, as seen in the Round Island petrel system, therefore provide ideal opportunities to explore the degree of within- and between-individual variation in migratory behaviour.
Although migration strategies were highly variable across the population, individuals used distinct areas within the overall range in a consistent manner across years. Individual consistency in space use has been shown for many migratory bird species, from great reed warblers (Acrocephalus arundinaceus) to Atlantic puffins (Fratercula arctica) [20, 22], however, this topic has received very little attention for migratory birds breeding in the tropics (e.g., ). Our results suggest that familiarity with locations may be more beneficial than tracking current environmental conditions even when resources may be patchy and unpredictable, such as those in tropical systems [2, 59]. Further, the fact that individuals which depart closer together in time are likely to have more similar non-breeding distributions may suggest that the meteorological or oceanographic conditions at departure can influence individual non-breeding distributions (but see ). However, we also observed temporal variation in space use, with EMD values being higher (i.e., migrations were less similar) when split into stages/periods, compared to the migration as a whole. This suggests a degree of temporal flexibility, such that individuals use the same areas in different years, but not necessarily at the same time during the non-breeding period, which has also been shown for other seabirds [37, 54]. This temporal flexibility therefore seems to only occur within the range of known areas for a particular individual, suggesting that relying on familiar areas is more beneficial than switching to a new location , and implies that temporal variation in resource availability may not be very large at these scales. Additionally, as gadfly petrel flying behaviour is strongly affected by wind conditions [50, 55] and thus how fast petrels move through their environment may vary with the environmental conditions that they experience in each year. The EMD metric provided an effective method for measuring spatial dissimilarity in non-overlapping distributions. This is particularly valuable for systems with large, non-overlapping variation in the possible range of individual migratory distributions, such as Round Island petrels. This framework could easily be applied to other species to compare space-use patterns within and across taxa (although effects of scale-dependence may need to be considered).
Repeated tracking of individual petrels also indicated a high degree of consistency in migratory timings. This has been reported for many other species (mean ICC = 0.408 (95% confidence interval: 0.3–0.5; reviewed by Franklin et al., in review) suggesting consistent individual differences in migration phenology to be a common feature of migratory systems. Recently, Trindade petrels breeding on Trindade Island in the South Atlantic Ocean have also shown to maintain their breeding schedules year-to-year . For tropical seabirds, breeding phenology can range from seasonal and synchronised breeders (e.g., Barau’s petrel; ), to aseasonal breeders, albeit in varying numbers across the year (e.g., Round Island petrels; ). Consequently, it is important to note that repeatability of individual phenologies may be naturally inflated when a large number of viable phenologies exist in a population. Despite this, Round Island petrels were still remarkably repeatable in their migratory timings with much lower within- than between-individual variation. While the exact breeding status and/or outcome of individual petrels on Round Island is rarely known, other studies on seabirds have shown that failed breeders and non-breeders may depart earlier from their colony in comparison to successful breeders [42, 62]. Regardless of this, we still found high repeatability in migratory departure (although not as high as arrival) without accounting for breeding outcome (i.e., success, failure, or if breeding was attempted). Calculating repeatability in migratory timings requires consistent methods of classifying phenological events. In our study, we used HMMs in order to assign dates of arrival and departure from Round Island in an objective and reproducible manner, which was particularly important given that the low spatial resolution available from GLS tracks can make identification of departure times from location data alone problematic.
While most individuals seem to follow a consistent migratory schedule, the differing levels of within-individual variation between the seasons suggests an influence of prevailing environmental conditions on timings of departure and arrival. Petrels arriving at Round Island in the winter (which show more consistent timings of arrival) are likely to breed during the peak breeding period, and could thus experience greater competition for resources, including nest sites, than birds arriving during the summer period (which are less consistent in arrival times). Arriving at a consistent time each year may facilitate synchronous mate arrival , which may be particularly important if pairs are to compete for nest sites. Although the petrels nest in a range of conditions across the island, most nests occur within a relatively small number of colonies. Nesting sites within these colonies are likely to be in high demand, particularly during the peak breeding period, when broken eggs with peck marks and young chicks with head wounds are often observed . Observations from camera traps have also shown intra- and inter- specific fights at petrel nest sites (Franklin, pers. obs.). As the Round Island petrel is a hybrid species complex [4, 6, 7], it means that the population comprises individuals with a great deal of genetic variation dictated by evolutionary histories. This, together with the fact that petrels have an asynchronous breeding period, means that interactions with different environmental conditions (such as the semi-annual wind reversals as a result of the two monsoon periods in the Indian Ocean ), may have given rise to the diverse range of migration patterns. However, it is not yet clear whether there is any temporal structuring in the genotypes of petrels on Round Island (i.e., are certain petrel hybrids on the island at certain times of year), or if different genotypes have different migratory distributions, which may contribute to individual phenological and spatial variation.
The small amount of within-individual variation suggests that consistency in migratory behaviours is favoured even in comparatively patchy and unpredictable tropical systems . This consistency, together with the fact that birds can be found breeding on Round Island all year round, means that different individuals are potentially exposed to different environmental conditions and human-associated impacts, with potentially important consequences for breeding success (e.g., ), survival and, ultimately, the status of this population. The Round Island petrel population appears to have arisen relatively recently in time through range expansions of different Pterodroma taxa [4, 6, 7]. The high level of individual migratory consistency means that future changes in non-breeding distributions and timings will most likely reflect changes in the numbers of individuals undertaking different journeys. Determining what is driving the large levels of between-individual variation in this system will be key in revealing the implications of individual consistency for population demography, and the potential consequences of future environmental changes across the migratory range.
Availability of data and materials
Tracking data are deposited on the Seabird Tracking database (www.seabirdtracking.org) handled by BirdLife International (ID: 1810). R scripts to reproduce analysis can be found on Github (https://github.com/kirstyfranklin/RIpetrel-consistency).
Earth mover’s distance
Generalised linear model
Hidden Markov model
Markov Chain Monte Carlo
Sea surface temperature
Abrahms B, Hazen EL, Aikens EO, Savoca MS, Goldbogen JA, Bograd SJ, et al. Memory and resource tracking drive blue whale migrations. Proc Natl Acad Sci USA. 2019;116:5582–7. https://doi.org/10.1073/pnas.1819031116.
Ashmole N. Seabird ecology and the marine environment. New York: Academic Press; 1971.
Au DWK, Pitman RL. Seabird interactions with dolphins and tuna in the eastern Tropical Pacific. Condor. 1986;88:304–17. https://doi.org/10.2307/1368877.
Booth Jones KA, Nicoll MAC, Raisin C, Dawson DA, Hipperson H, Horsburgh GJ, et al. Widespread gene flow between oceans in a pelagic seabird species complex. Mol Ecol. 2017;26:5716–28. https://doi.org/10.1111/mec.14330.
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.
Brown RM, Nichols RA, Faulkes CG, Jones CG, Bugoni L, Tatayah V, et al. Range expansion and hybridization in Round Island petrels (Pterodroma spp.): evidence from microsatellite genotypes. Mol Ecol. 2010;19:3157–70. https://doi.org/10.1111/j.1365-294X.2010.04719.x.
Brown RM, Jordan WC, Faulkes CG, Jones CG, Bugoni L, Tatayah V, et al. Phylogenetic relationships in Pterodroma petrels are obscured by recent secondary contact and hybridization. PLoS ONE. 2011;6: e20350. https://doi.org/10.1371/journal.pone.0020350.
Brown JM, van Loon EE, Bouten W, Camphuysen KCJ, Lens L, Müller W, et al. 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.
Calenge C. The package adehabitat for the R software: a tool for the analysis of space and habitat use by animals. Ecol Model. 2006;197:516–9.
Ceia FR, Ramos JA. Individual specialization in the foraging and feeding strategies of seabirds: a review. Mar Biol. 2015;162:1923–38. https://doi.org/10.1007/s00227-015-2735-4.
Cerveira LR, Ramos JA, Rodrigues I, Almeida N, Araújo PM, dos Santos I, et al. Inter-annual changes in oceanic conditions drives spatial and trophic consistency of a tropical marine predator. Mar Environ Res. 2020;162: 105165. https://doi.org/10.1016/j.marenvres.2020.105165.
Dell’Ariccia G, Benhamou S, Dias MP, Granadeiro JP, Sudre J, Catry P, et al. Flexible migratory choices of Cory’s shearwaters are not driven by shifts in prevailing air currents. Sci Rep. 2018;8:3376. https://doi.org/10.1038/s41598-018-21608-2.
Delord K, Barbraud C, Pinaud D, Ruault S, Patrick SC, Weimerskirch H. Individual consistency in the non-breeding behaviour of a long-distance migrant seabird, the Grey petrel Procellaria cinerea. Mar Ornithol. 2019;47:93–103.
Dias MP, Granadeiro JP, Phillips RA, Alonso H, Catry P. Breaking the routine: individual Cory’s shearwaters shift winter destinations between hemispheres and across ocean basins. Proc R Soc B. 2011;278:1786–93. https://doi.org/10.1098/rspb.2010.2114.
Fayet AL, Freeman R, Shoji A, Boyle D, Kirk HL, Dean BJ, et al. Drivers and fitness consequences of dispersive migration in a pelagic seabird. Behav Ecol. 2017;2017(27):1061–72. https://doi.org/10.1093/beheco/arw013.
Fieberg J, Kochanny CO. Quantifying home- range overlap: the importance of the utilization distribution. J Wildl Manag. 2015;69:1346–59. https://doi.org/10.2193/0022-541X(2005)69[1346:QHOTIO]2.0.CO;2.
Gill JA, Alves JS, 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.
Gilroy JJ, Gill JA, Butchart SH, Jones VR, Franco AM. Migratory diversity predicts population declines in birds. Ecol Lett. 2016;19:308–17. https://doi.org/10.1111/ele.12569.
Grist H, Daunt F, Wanless S, Nelson EJ, Harris MP, Newell M, et al. Site fidelity and individual variation in winter location in partially migratory European shags. PLoS ONE. 2014;9: e98562. https://doi.org/10.1371/journal.pone.0098562.
Guilford T, Freeman R, Boyle D, Dean B, Kirk H, Phillips R, et al. A dispersive migration in the Atlantic puffin and its implications for migratory navigation. PLoS ONE. 2011;6: e21336. https://doi.org/10.1371/journal.pone.0021336.
Gunnarsson TG, Gill JA, Sigurbjörnsson T, Sutherland WJ. Pair bonds: arrival synchrony in migratory birds. Nature. 2004;431:646. https://doi.org/10.1038/431646a.
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 timing of migration. J Avian Biol. 2017;48:91–102. https://doi.org/10.1111/jav.01292.
Hijmans RJ, Williams E, Vennes C. Geosphere: spherical trigonometry. R package version 1.5-10; 2019
Jaquemet S, Le Corre M, Weimerskirch H. Seabird community structure in a coastal tropical environment: importance of natural factors and fish aggregating devices (FADs). Mar Ecol Prog Ser. 2004;268:281–92. https://doi.org/10.3354/meps268281.
Jaeger A, Feare CJ, Summers RW, Lebarbenchon C, Larose CS, Le Corre M. Geolocation reveals year-round at-sea distribution and activity of a superabundant tropical seabird, the sooty tern Onychoprion fuscatus. Front Mar Sci. 2017;4:394. https://doi.org/10.3389/fmars.2017.00394.
Kumar SP, Narvekar J, Nuncio M, Gauns M, Sardesai S. What drives the biological productivity of the Northern Indian Ocean? In: Wiggert JD, Hood RR, Naqvi SWA, Brink KH, Smith SL, editors. Indian Ocean biogeochemical processes and ecological variability. San Francisco: American Geophysical Union; 2009. p. 33–56.
Kranstauber B, Smolla M, Safi K. Similarity in spatial utilization distributions measured by the earth mover’s distance. Methods Ecol Evol. 2017;8:155–60. https://doi.org/10.1111/2041-210X.12649.
Kranstauber B, Smolla M, Scharf AK. move: visualising and analysing animal track data. R package version 4.0.6; 2020. https://CRAN.R-project.org/package=move
Leal GR, Bugoni L. Individual variability in habitat, migration routes and niche used by Trindade petrels. Pterodroma arminjoniana. Mar Biol. 2021;168:134.
Le Corre M, Jaeger A, Pinet P, Kappes MA, Weimerskirch H, Catry T, et al. Tracking seabirds to identify potential Marine Protected Areas in the tropical Western Indian Ocean. Biol Conserv. 2012;156:83–93. https://doi.org/10.1016/j.biocon.2011.11.015.
Lenth RV. emmeans: Estimated Marginal Means, aka least-squares means. R package version 1.6.1; 2021. https://CRAN.R-project.org/package=emmeans
Lisovski S, Hewson CM, Klaassen RH, Korner-Nievergelt F, Kristensen MW, Hahn S. Geolocation by light: accuracy and precision affected by environmental factors. Methods Ecol Evol. 2012;3:603–12. https://doi.org/10.1111/j.2041-210X.2012.00185.x.
Lisovski S, Wotherspoon S, Sumner M. TwGeos: Basic data processing for light-level geolocation archival tags. R package version 0.1.2; 2016. https://github.com/slisovski/TwGeos/
López-López P. Individual-based tracking systems in ornithology: welcome to the era of big data. Ardeola. 2016;63:103–136. https://doi.org/10.13157/arla.63.1.2016.rp5
McFarlane Tranquilla LA, Montevecchi WA, Fifield DA, Hedd A, Gaston AJ, Roberston GJ, et al. Individual winter movement strategies in two species of murre (Uria spp.) in the Northwest Atlantic. PLoS ONE. 2014;9:90583. https://doi.org/10.1371/journal.pone.0090583.
McNamara JM, Dall SRX. Information is a fitness enhancing resource. Oikos. 2010;119:231–6. https://doi.org/10.1111/j.1600-0706.2009.17509.x.
Merkel B, Descamps S, Yoccoz NG, Grémillet D, Daunt F, Erikstad KE, et al. Individual migration strategy fidelity but no habitat specialization in two congeneric seabirds. J Biogeogr. 2020;48:263–75. https://doi.org/10.1111/jbi.13883.
Møller AP, Rubolini D, Lehikoinen E. Populations of migratory bird species that did not show a phenological response to climate change are declining. Proc Natl Acad Sci USA. 2008;105:16195–200. https://doi.org/10.1073/pnas.0803825105.
Nicoll MAC, Nevoux M, Jones CG, Ratcliffe N, Ruhomaun K, Tatayah V, et al. Contrasting effects of tropical cyclones on the annual survival of a pelagic seabird in the Indian Ocean. Glob Chang Biol. 2017;23:550–65. https://doi.org/10.1111/gcb.13324.
Oppel S, Weber S, Weber N, Fox D, Leat E, Sim J, et al. Seasonal shifts in foraging distribution due to individual flexibility in a tropical pelagic forager, the Ascension frigatebird. Mar Ecol Prog Ser. 2017;585:199–212. https://doi.org/10.3354/meps12377.
Phillips RA, Silk JRD, Croxall JP, Afanasyev V, Briggs DR. Accuracy of geolocation estimates for flying seabirds. Mar Ecol Prog Ser. 2004;266:265–72. https://doi.org/10.3354/meps266265.
Phillips RA, Silk JRD, 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.
Pinet P, Salamolard M, Probst J, Russell J, Jaquemet S, Le Corre M. Barau’s Petrel Pterodroma baraui: history, biology and conservation of an endangered endemic petrel. Mar Ornithol. 2009;37:107–13.
Pinet P, Jaquemet S, Pinaud D, Weimerskirch H, Phillips RA, Le Corre M. Migration, wintering distribution and habitat use of an endangered tropical seabird, Barau’s petrel Pterodroma baraui. Mar Ecol Prog Ser. 2011;423:291–302. https://doi.org/10.3354/meps08971.
Robinson RA, Crick HQ, Learmonth JA, Maclean I, Thomas CD, Bairlein F, et al. Travelling through a warming world: climate change and migratory species. Endanger Species Res. 2009;7:87–99. https://doi.org/10.3354/esr00095.
Rubner Y, Tomasi C, Guibas LJ. The earth mover’s distance as a metric for image retrieval. Int J Comput Vis. 2000;40:99–121. https://doi.org/10.1023/A:1026543900054.
R Core Team. R: a language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria; 2019.
Schott FA, McCreary JP. The monsoon circulation of the Indian Ocean. Prog Oceanogr. 2001;51:1–123. https://doi.org/10.1016/S0079-6611(01)00083-0.
Spear LB, Balance LT, Ainley DG. Response of seabirds to thermal boundaries in the tropical Pacific: the thermocline versus the Equatorial Front. Mar Ecol Prog Ser. 2001;219:275–89. https://doi.org/10.3354/meps219275.
Spear LB, Ainley DG. Flight behaviour of seabirds in relation to wind direction and wing morphology. Ibis. 1997;139:221–33. https://doi.org/10.1111/j.1474-919X.1997.tb04620.x.
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.
Tatayah V. The breeding biology of the Round Island Petrel (Pterodroma arminjoniana) and the factors determining breeding success. MPhil Thesis, University of Mauritius; 2007.
Tatayah V. The breeding biology of the Round Island petrel. PhD Thesis, University of Mauritius; 2010.
van Bemmelen R, Moe B, Hansen SA, Schmidt NM, Hansen J, Lang J, et al. Flexibility in otherwise consistent non-breeding movements of a long-distance migratory seabird, the long-tailed skua. Mar Ecol Prog Ser. 2017;578:197–211. https://doi.org/10.3354/meps12010.
Ventura F, Granadeiro JP, Padget O, Catry P. Gadfly petrels use knowledge of the windscape, not memorized foraging patches, to optimize foraging trips on ocean-wide scales. Proc R Soc B. 2020;287:20191775. https://doi.org/10.1098/rspb.2019.1775.
Visser I, Speekenbrink M. depmixS4: an R package for hidden markov models. J Stat Softw. 2010;36:1–21.
Wakefield ED, Phillips RA, Matthiopoulos J. Quantifying habitat use and preferences of pelagic seabirds using individual movement data: a review. Mar Ecol Prog Ser. 2009;391:165–82. https://doi.org/10.3354/meps08203.
Walther GR, Post E, Convey P, Menzel A, Parmesan C, Beebee TJC, et al. Ecological responses to recent climate change. Nature. 2002;416:389–95. https://doi.org/10.1038/416389a.
Weimerskirch H. Are seabird foraging for unpredictable resources? Deep-Sea Res Part II Top Stud Oceanogr. 2007;54:211–23. https://doi.org/10.1016/j.dsr2.2006.11.013.
Winger BM, Auteri GG, Pegan TM, Weeks BC. A long winter for the Red Queen: rethinking the evolution of seasonal migration. Biol Rev. 2018;94:737–52. https://doi.org/10.1111/brv.12476.
Wotherspoon SJ, Sumner MD, Lisovski S. R package SGAT: solar/satellite geolocation for animal tracking. GitHub Repository; 2013. http://github.com/swotherspoon/sgat
Yamamoto T, Takahashi A, Sato K, Oka N, Yamamoto M, Trathan PN. Individual consistency in migratory behaviour of a pelagic seabird. Behaviour. 2014;151:683–701. https://doi.org/10.1163/1568539X-00003163.
This work would not have been possible without the many staff and volunteers from the Mauritian Wildlife Foundation (MWF) and National Parks and Conservation Service (NPCS) who have helped with the Round Island petrel research programme. In particular, we would like to thank Vimul Nundlaul, Richard Baxter, Pat Banville, Katherine Booth Jones, Lucy Rouse, Helen Gath, Luke A’Bear and Johannes Chambon in facilitating the deployment and recovery of GLS on the island. We would also like to thank Stephen Vickers and Martin Stoffel for statistical advice, and the associate editor, Matthieu Le Corre and an anonymous reviewer who provided valuable feedback that greatly improved the manuscript. The main EMD analysis was conducted on the High Performance Computing Cluster supported by the Research and Specialist Computing Support service at the University of East Anglia.
KAF was funded by a John and Pat Warham PhD studentship from the British Ornithologists’ Union (BOU). The petrel tracking programme was supported by the Natural Environmental Research Council (NERC; UK) (Grant NE/H5081500) and Research England. JAG was also supported by NERC (Grant NE/M012549/1).
Ethics approval and consent to participate
This work involved the tagging (rings and geolocators) and recapture of petrels. All work followed the ethical standards set out by The Mauritian Wildlife Foundation and its partner and consulting organizations, the North of England Zoological Society, the Durrell Wildlife Conservation Trust and the International Zoo Vet Group. In addition, geolocator deployment and recovery were approved by the Zoological Society of London Animal Ethics Committee.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
The original version of this article was revised: the incorrect version of Additional file 1 was replaced
: Details of all geolocator deployments and recoveries on adult Round Island petrels from 2009 to 2019, and number of complete migrations which took place in each petrel year for the 62 petrels with repeat migrations. Table S2: Results of ANOVA tests for generalised linear model selection for the similarity of petrel migrations within- and between- individuals, using a) the earth mover’s distance (EMD) ‘effort’ values, and b) Bhattacharyya’s affinity (BA), and the spatiotemporal similarity of within-individual petrel migrations when split into c) six equal size stages, and d) 30-day periods, both using EMD. Significant effects (p < 0.05) are highlighted in bold. Table S3: Results of generalised linear model to investigate the similarity of petrel whole migrations within- and between- individuals, using Bhattacharyya's affinity (BA) values. Minimum adequate model is shown. Note, the binary categorical variable ‘same individual’ is being compared to the reference level of that variable, which is 0 (different individuals). Significant effects (p < 0.05) are highlighted in bold. Table S4: Pairwise comparisons between each level of the categorical variable ‘stage’ from the generalised linear model examining the spatiotemporal similarity of petrel migrations when split into six equal size stages, using earth mover’s distance (EMD) ‘effort’ values. Significant effects (p < 0.05) are highlighted in bold. Table S5: Results of generalised linear model to investigate the spatiotemporal similarity of within-individual petrel migrations when split into 30-day periods, using earth mover’s distance (EMD) ‘effort’ values. Minimum adequate model is shown and categorical variable ‘period’ is being compared to reference level of whole migration. Significant effects (p < 0.05) are highlighted in bold. Table S6: Pairwise comparisons between levels of the categorical variable ‘30-day period’ from the generalised linear model examining the spatiotemporal similarity of petrel migrations when split into 30-day periods, using earth mover’s distance (EMD) ‘effort’ values. Significant effects (p < 0.05) are highlighted in bold. Table S7: Colour of geographic coordinates in figures (main text and Additional File 1) and the corresponding year/s of each migration. Table S8. : The variance components and repeatability estimates (R) from adult Round Island petrels with repeated tracks (2-5 years) for departure date from the breeding colony, arrival at the colony and duration of the migratory period. Figure S1: Earth mover’s distance (EMD) values are calculated for both A) consecutive 30-day periods, irrespective of start date, and B) by splitting the migration into six equal stages. For A), this typically means that a short period of longer within-individual migrations is not included, however, there was evidence of individual consistency in migration duration and so this often only comprises a small proportion of the whole migration. Whereas for B), if one migration is longer than the other then each stage will comprise of more days than the other corresponding stage. Figure S2: Geographic coordinates (a) and tracks (b) from 62 adult Round Island petrels that have been tracked over multiple complete migrations (n = 131) with geolocators. Colours of lines and points represent the year of tracking (2009-17). Black diamond indicates the location of Round Island, Mauritius. Figure S3. Example tracks of between-individual Round Island petrel migration comparisons with the same relative overlap value (Bhattacharyya's affinity (BA) = 0.06), but with comparatively high (A), moderate (B), and low (C) earth mover’s distance ‘effort’ values. Positions denote twice-daily median locations with the two different years illustrated in different colours. Black diamond indicates the location of Round Island, Mauritius. The tracking year that each set of colours represents can be found in Additional File 1: Table S7.
About this article
Cite this article
Franklin, K.A., Norris, K., Gill, J.A. et al. Individual consistency in migration strategies of a tropical seabird, the Round Island petrel. Mov Ecol 10, 13 (2022). https://doi.org/10.1186/s40462-022-00311-y