Habitat availability influences migration speed, refueling patterns and seasonal flyways of a fly-and-forage migrant
Movement Ecology volume 8, Article number: 10 (2020)
Despite our understanding of the principal factors that shape bird migration strategies, there is conflicting evidence regarding the role of habitat in shaping migration routes and schedules, including day and night activity and differences between autumn and spring. For fly-and-forage migrants, we predict that habitat characteristics might guide migration speed, route selection and migrating schedules.
We use solar-powered GPS transmitters, obtaining high accuracy data, to monitor the migratory movements of Eleonora’s falcon breeding in Cyprus, which is the easternmost breeding population of the species. We tested for potential preferences in habitat characteristics along the migration routes, separately for the northern, drier part and the more vegetated southern part of the trips. We also examined the relationship between migration speed and vegetative cover during day and at night, accounting for wind support.
We found that tagged individuals repeatedly exhibited an anticlockwise loop migration pattern with spring routes being more easterly than autumn ones. We identified a preference for migration through vegetation-rich areas, where during daytime tagged individuals travel at slower migration speeds compared to vegetation-poor areas, indicating fly-and-forage activity. Birds roosted during most nights, combining refueling stopovers at selected vegetation-rich areas before or after crossing ecological barriers. Conversely, both during day and night, tagged individuals overflew unsuitable habitats more quickly.
Our results suggest that habitat is an important factor in Eleonora’s falcon migratory strategies. Active selection of vegetation rich areas in combination with reduced migration speeds there, allows the migrating falcons to combine migration during the day with fly-and-forage refueling, while roosting most nights.
Migration is part of the annual cycle of many bird species that has evolved over millennia [1,2,3]. The main driving force is the exploitation of food resource fluctuations at certain times of year  avoiding resource depletion especially at breeding grounds . Current bird migration patterns are thought to have evolved within the last 15,000 years, during postglacial species range expansions [2, 4, 6], resulting in complex movement systems [1, 3, 7, 8]. Despite their complexity, migration systems are characterized by high precision and accuracy, requiring high spatiotemporal understanding by the migrating birds [3, 9, 10]. Like any form of movement, migrations emerge through the interplay of internal and external factors , including genetic instruction, physiological processes and behavioral adaptations to external conditions such as ecological barriers [1, 2, 6, 7, 11, 12]. Yet, despite our understanding of the factors that generally dictate migration, less is known regarding how landscape characteristics, such as elevation, tree cover and habitat type, influence migratory movements . Habitat might influence route selection and migration speed, resulting in differences in behavioural patterns between day and night, and autumn and spring . In addition, individual repeatability between successive migratory trips might also be attributed to vegetation characteristics .
Avian migration strategies are guided by pre-migratory fuel deposition, but refueling en route is equally critical . Integrating movement and environmental data can greatly enhance our understanding of how migrant birds balance travel with refueling opportunities [16,17,18]. An optimal strategy does not entail stopovers at all potentially fruitful refueling sites, which are instead selected according to fuel load, and perhaps using prior knowledge on food availability at specific areas over successive trips [15, 19, 20]. This combination can result in different migrating patterns, involving longer or shorter travel/refueling stints . Furthermore, fly-and-forage activity, i.e. hunting while keeping to the general direction of travel, is specifically beneficial to aerial insectivores, which catch and consume prey on the wing . Such species must therefore select flyways with insect-abundant habitats where they can keep on feeding while moving . Ecological obstacles en route that do not offer feeding opportunities are either detoured or rapidly overflown . While detours lengthen total distance travelled , they do not necessarily increase duration, as they could result from route optimization involving energy refueling and/or tailwind support [15, 26, 27]. Food availability and wind patterns are also thought to play a role in seasonal differences found in migration speed and trip duration [19, 28], as well as to observed loop patterns across autumn and spring migration [15, 29]. In addition, overall journey duration may be influenced by whether or not birds migrate at night . Nocturnal migration typically provides more favorable climatic conditions, such as lighter winds and less turbulence, making flight more efficient [11, 22, 31]. For diurnally feeding species, nocturnal flight does not involve foraging activity , thus faster migration speed at night is expected . As a strategy, nocturnal continuous flight with foraging during the day may result in shorter migration times . Raptors, however, are not in general expected to fly at night  unless they are crossing large water bodies and need to extent their flight overnight [31, 33]. Nevertheless, high resolution tracking data might reveal more about nocturnal flying activity .
Eleonora’s falcon (Falco eleonorae) is a complete, long-distance, trans-equatorial migrant species [4, 11, 23], breeding primarily on islands and islets in the Aegean Sea, which is presumed to be its original range and the centre of its distribution . The breeding population spreads from Cyprus in the east, westwards along the Mediterranean Sea, the Atlantic coast of Morocco and the Canary Islands . Almost the entire population overwinters on the island of Madagascar , with a small proportion in the Mascarene Islands  and East Africa . The species has evolved to take advantage of unique breeding and wintering niches in the Mediterranean and Madagascar respectively, utilizing abundant resources and avoiding interspecific competition . Eleonora’s falcon feeds predominantly on insects most of the year, though during the breeding season it feeds primarily on migrating birds . The delayed breeding period of Eleonora’s falcon coincides with the peak of autumn bird migration passage, with falcons catching migrant birds on the wing in order to feed themselves and their nestling during the offspring rearing period . Previous telemetry studies on Eleonora’s falcon migration have shown that tagged individuals exhibited increased migration speeds over ecological barriers, compared to the slower migration speeds before and after crossing those barriers to refuel [39, 40]. Previous work also found that Eleonora’s falcons migrate during the day and at night throughout their trips . In addition, an anti-clockwise migration pattern has been observed in Eleonora’s falcon between seasons, with spring routes being more easterly than autumn routes .
In this study, we utilized transmitters with GPS for accurate localization, providing also more consistent temporal data resolution compared to the transmitters previously used on Eleonora’s falcon migration studies [24, 39,40,41]. With those higher accuracy transmitters, we aimed to identify the extent to which migratory routes of Eleonora’s falcon target favored habitats with rich vegetation for stopovers and for refueling using a fly-and-forage strategy, adjusting their migration speeds accordingly. In order to examine specifically whether Eleonora’s falcon actively selects particular habitat characteristics, we compared use of those characteristics along the falcons’ migratory routes with the availability of those habitat features within their seasonal migration corridors. We examined if these habitat preferences could possibly vary between the northern (mostly desert) and southern (more vegetated) portions of their migration trips, as well as between seasons. In addition, if habitat influences flight and foraging decisions, we expect to find differences in migration speeds between habitats with different flying modes employed during the day and at night, such as soaring vs. flapping flight and fly-and-forage. If individuals move through different habitat types without changing their behaviour, especially their migration speed, this would suggest there is no habitat preference. Another possibility is that weather conditions may affect migration times, even resulting in stopovers when weather conditions are particularly harsh. We expected individuals to actively migrate at a faster migration speed, both during the day and at night over habitats unsuitable for foraging, such as extensive water bodies and deserts. We expect them though to migrate more slowly over insect-rich habitats and roost during the night therein. To that end we tested for the effect of vegetative cover and wind conditions on migration speeds, and examined how these effects could vary between the day and night and between seasons. In addition, obtaining data for more than one consecutive migration cycle for two individuals provided an insight into the spatial and temporal route repeatability of consecutive trips.
Cyprus hosts around 130 pairs of Eleonora’s falcon nesting on the south coast of the island, from Cape Gata on the Akrotiri Peninsula in the east, to Cape Aspro in the west . Between 2013 and 2017 we monitored the year-round movements of Eleonora’s falcons originating from Cyprus using for the first time GPS telemetry technology providing high spatial and temporal detail. In particular, we attached 12 transmitters on individuals breeding in the Akrotiri colony in Cyprus, of which eight, comprising five female adults, two male adults and one juvenile, provided us with migration related data (Table 1).
Birds were caught using mist nets at nesting and foraging locations . Mist nets were placed on a cliff top opposite nests at the southern edge of the Akrotiri peninsula. A stuffed Eagle owl (Bubo bubo) was used as a decoy . In addition, four nets were set within an olive grove, where the falcons had been observed feeding on large beetles at dusk. Birds were handled in accordance with standard procedures  and marked with metal rings , while biometric measurements including body mass were taken . Sex and age was determined using colouration of bare parts and plumage characteristics . Transmitters were attached with a teflon harness as a backpack . Two different types of solar transmitters were used: ten GPS-GSM loggers (SKUA model; Ecotone Telemetry, Poland), and two Platform Transmitter Terminals (PTTs) (5 g solar PTT model; Microwave Telemetry Inc., Columbia USA). The weight of the transmitter ranged between 1.5 and 5% of the bird’s mass at capture [41, 50, 51]. No abnormal behavior was observed after transmitter attachment , with individuals moving as expected between breeding and foraging areas over subsequent days. GPS loggers collected locations on a user defined temporal pattern (2–6 h), and sent them through the GSM network, while PTT transmitters were set on a pre-programmed standard duty cycle of 10 h on followed by 48 h off .
The positional error of GPS locations was less than 20 m for 80% of the retrieved points , we thus retained all available points for further analyses. The Kalman filtering algorithm was used to enhance the varying accuracy of PTT positions, which is divided in seven location classes (i.e., LC3 < 250 m, LC2 250–500 m, LC1 500–1500 m, LC0 > 1500 m, LCA and LCB = Unbounded accuracy, LCZ = Invalid location) . Although LC3, LC2 and LC1 are typically considered as high quality (accuracy) locations , it was shown  that LCA had similar accuracy with LC1 while LC0 were far less accurate than LCA. Thus, in cases where we had to choose between which points to exclude (see subsequent filtering) we used the following order of priority: B, 0, A, 1, 2, 3. Consequently, points of all location classes apart from LCZ were used as long as they followed the general travel direction, assessed through visual inspection of the points in Google Earth [14, 44, 57, 58]. We excluded locations for both GPS and PTT transmitters that were obtained less than an hour apart to standardize sampling interval and for excessive autocorrelation avoidance [36, 59]. Furthermore, because of lower accuracy of PTT transmitters, we additionally excluded locations less than a kilometer apart [36, 60].
Estimation of path metrics
Because of the lower accuracy of PTT obtained locations, we decided to rely only on GPS positions for path metric estimates using all available segments. Segment length between two successive telemetry fixes was calculated based on the geodesic distance , i.e. the shortest line between any two points on the earth’s surface considering the curvature of the earth ). In order to assign a position to day or night, sunset and sunrise information was obtained for each position based on local times , with half an hour before sunrise and after sunset included as daytime, because the species is also known to hunt at dawn and dusk .
Bird activity and flying mode between two successive telemetry fixes (e.g. resting, soaring, flapping and fly-and-forage) and related tortuosity of the tracks affect path metrics . Therefore in this study, we define migration speed as the geodesic distance between the endpoints (successive telemetry fixes) of a segment, divided by the time between those points. Accordingly, calculated migration speed refers to the rate at which the migrant progresses along the geodesically defined path between two successive telemetry fixes. We considered active migratory movements those path segments where migration speed was >5 km/h [24, 31, 40], though such migration speeds might also include fly-and-forage activity. Because of potential overlap between fly-and-forage activity and forage-free migration speeds, we did not set a firm travel speed threshold to distinguish between these behavioural categories. Instead we simply considered that the higher the migration speed (accounting for wind support), the less likely the bird was foraging in flight.
Habitat use at stopovers and roosts
As stopovers we considered areas with a 25 km radius where an individual stayed for at least 24 h without exhibiting directional migratory movement, allowing for one outlier point per stopover [44, 65]. This radius is sufficient to delineate stopovers as if birds exhibited directional movement towards the target destination, even with the lowest active migration speed of 5 km/h [24, 31, 40], flying only between sunrise and sunset (approximately 12 h), would cross this 25 km radius area during a 24 h period. Habitat use during night roosts, based on consecutive positions indicating that individuals were stationary, was identified by overlaying obtained migration positions onto a Moderate-resolution Imaging Spectroradiometer (MODIS) land cover layer . Two resolutions of land cover were used; a finer available one with a resolution of 500 × 500 m (MCD12Q1) to reflect high GPS accuracy locations and a coarser one with about 5 × 5 km resolution (MCD12C1) to match the lower accuracy of PTT transmitters. To further utilize the high spatial accuracy of the GPS locations, Google Earth  was used to identify habitat details at roost and stopover areas, such as single trees within croplands, which could not be identified within remotely sensed habitat datasets .
Route selection in relation to habitat characteristics
In order to determine whether Eleonora’s falcon selects specific habitat characteristics within its seasonal migration corridors, we compared actual bird positions with random locations generated within a 50 km buffer either side of each trip . This distance is expected to be well within the range of vision of tracked individuals . For each trip we therefore generated 150 random points in Geospatial Modelling Environment (GME) 7.2.1  to provide a total number similar to the total number of recorded positions for use in statistical analyses. Because of the lack of specific information on food availability especially along tracks of long-distance migrants , we used the Normalized Difference Vegetation Index (NDVI) obtained from the MODIS sensor on the Terra satellite  as a proxy for the greenness along the migratory routes [68, 73] and as an indication of food availability for the migrating falcons [40, 74]. Sixteen-day NDVI granules (nominal resolution 250 × 250 m) were downloaded covering the extent of migration period and the granule that overlapped most with the duration of each trip event was used to extract the relationships separately. Data (actual and random) were also overlaid over a percent tree cover layer (hereafter ‘tree cover’) i.e. the vegetation continuous field product (nominal resolution 250 × 250 m) also from the MODIS sensor . This more stable feature of the landscape, as well as elevation, have been found to influence movements and foraging activity of wintering falcons in Madagascar . Elevation data were obtained from the United States Geological Survey website and had a resolution of 7.5 arc-sec, i.e. about 150 × 150 m . All environment data were extracted using ArcGIS 10.1 . For all relationships of locations with environmental data (NDVI, percent tree cover and elevation), we extracted bilinear interpolated data, which incorporate adjacent cell values in the value calculation, to allow for a uniform approach for both the higher resolution GPS data and the lower resolution PTT obtained data. For habitat selection analyses we excluded actual locations and random points over the sea and other extensive water bodies for which NDVI and tree cover were not available. We investigated which environmental characteristics affect route selection of Eleonora’s falcon during migration with the use of Generalized Linear Mixed Models (GLMM) with a binomial distribution and logit link function in lmer in R 3.5.1 . Specifically, we tested for the effects of NDVI, tree cover and elevation and their interactions with season as fixed effects, individual and year as crossed random factors, with a binary value for whether a point represented an actual position or a random point as the dependent variable. We ran these models separately for the northern, drier (e.g. Sahara desert), and southern (e.g. savannah and rain forest) parts of the migration journeys. To set latitude separating these north and south sections of the trips, we used the MODIS land cover layer  as an indicator to identify the desert and non-desert areas, setting the line at latitude 16oN. We used the Corrected Akaike’s Information Criterion (AICc) to select the best models for each test, considering all models with ΔAICc < 4 compared to the best supported model.
Migration speed differences in relation to season, habitat characteristics and daily activity
Because migration speed is related to the activity of tracked individuals, for example slower migration speeds are expected when birds are foraging, we investigated which factors affect the migration speed of Eleonora’s falcon during migration with the use of GLMMs with a Gaussian distribution with identity link function using lme4 in R 3.5.1 . Specifically, we tested for the effects on migration speed of tailwind, NDVI, percent tree cover (all centred and standardized, ) and season, as well as the interactions between the environmental variables and season. Individual and year were included as crossed random factors to account for spatial and temporal autocorrelation. Because of the expected differences in flight modes between day and night, with night flights not affected by thermals, and expectations of reduced foraging at night, we ran separate models without the transitioning segments between day and night. Specifically, we investigated how NDVI, tree cover, tailwind and season affect migration speed during the day and at night, including the interactions of each factor with day and night. We ran a separate model including all recorded flights using the dataset excluding transitioning segments, and another including only flights over 5 km/h to focus specifically on active migratory and fly-and-forage migration speeds, when landscape characteristics can influence movements differently compared to stationery periods, especially in relation to food availability. Because of the smaller sample size of segments over 5 km/h, we ran Bayesian GLMMs in blme. We used AICc to select the best models for each test, while also considering all models with ΔAICc < 4 compared to the best supported model.
Accounting for confounding wind effects
Knowing that wind will have an effect on migration speeds, and to make sure that the wind effect is accounted for when estimating the independent effect of habitat characteristics, we included wind as a fixed effect in our models. The loggers that we used did not provide flying altitude data, so to calculate tailwind at each GPS location we used winds at 850 hPa pressure level, corresponding to a mean altitude of 1500 m a.s.l.. In previous studies on raptors migrating through Africa a 925 hPa pressure level (mean altitude of 750 m a.s.l.) has been used, but most of those studies were conducted on the western African flyway (e.g. [29, 68, 79]), where ground elevations are considerably lower than for the eastern flyway studied here. Indeed, mean ground elevation of telemetry locations in the present study was 750 m a.s.l.. Data downloaded from ERA5, the latest atmospheric reanalysis product of the European Centre for Medium-Range Weather Forecasts (ECMWF) . U (east-west) and V (north-south) components for each telemetry location were downloaded in R using the NCEP.interp function in RNCEP . These two wind components were combined to provide tailwind speed in relation to the overall migration direction  of 170o degrees for autumn and 350o in spring, based on the shortest direct line potential migratory route from Akrotiri colony in Cyprus to the point of the narrowest crossing over the Mozambique Channel to Madagascar and vice versa.
Description of migratory routes
From the eight Eleonora’s falcons tagged that provided migration data (Table 1), we obtained data for seven migration seasons (four in autumn and three in spring) between 2013 and 2016, from a total of 17 individual migratory events (12 autumn and five spring). Following data filtering of low-quality positions and eliminating potential sources of autocorrelation, we retained 1716 positions (1578 GPS and 138 PTT).
Autumn migration duration was on average longer (N = 8, mean = 24.75 days, SD = 11.52) than spring migration (N = 3, mean = 19.33 days, SD = 4.16), though we note that the confirmed durations in spring came from only one individual (Farofylakas) for which the comparative results indicate the same average duration (autumn: N = 4, mean = 19 days, SD = 10.03; spring: N = 3, mean = 19.33 days, SD = 4.16). In autumn, departure from Cyprus occurred between 4th October and 6th November, with the peak around 26th October (eight out of 12 confirmed departure dates falling between 24th and 27th October). Arrival in Madagascar occurred between 8 and 28th November without any obvious peak. In spring, departure from Madagascar was between 9 and 15th April, and arrival in Cyprus between 25th April and 9th May (Fig. 1).
In autumn, upon initiating migration, tagged falcons appeared to fly non-stop south. Tagged individuals crossed the Mediterranean from Cyprus and from there most routes were over the Sinai Peninsula along the western coast of the Red Sea and over the Sahara Desert at its narrowest crossing in the east, aided by predominant northerly winds (See Fig. 2a for north-south winds and Additional file 1: Figure S1a for east-west winds) [3, 80]. After crossing the Sahara, individuals roosted for most nights from there onwards, including one or two stopovers for most routes within savannah in Ethiopia, Sudan or Kenya. All routes continued southwards towards Mozambique, from where on occasion individuals headed westwards for refuelling stopovers prior to the crossing of the Mozambique Channel to Madagascar (Fig. 3a). In spring, the tagged falcons crossed the Mozambique Channel, at times straight over Zanzibar to Tanzania. Then, aided also by predominant winds (See Fig. 2b for north-south winds and Additional file 1: Figure S1a for east-west winds), moved through eastern Kenya to Ethiopia and Somalia where tagged individuals had a short stopover before following a more easterly route than in autumn, primarily along the east coast of the Red Sea, avoiding Sahara predominant headwinds (Fig. 2b) [3, 80], and then across the Mediterranean to the breeding grounds in Cyprus (Fig. 3b). The migration speed over latitude relationship shows that there were sites at specific latitudes with concentrations of stopovers for roosting and foraging. The autumn concentration of points around 10–13° N, suggests that there was a major stopover upon reaching savannah after rapidly overflying the Sahara desert. After this refuelling stop there were further autumn concentrations representing stopovers, especially around −2° S and − 15° to − 17° S. The main stopover in spring appeared to occur around 8° to 9° N. In addition, very few night roosts were observed from 13° to 32° N, indicating continuous flight day and night (Fig. 4b). In autumn, tagged individuals roosted predominantly in savannah (62% of roosting nights), while in spring they roosted mostly in shrubland (65% of roosting nights) (Fig. 5).
Route selection and migration speed in relation to habitat characteristics
The average distance covered in autumn was 7139 km (SD = 1320), compared to 7245 km (SD = 401) in spring. The shortest path distance was 5749 km (SD = 209) in autumn and 5845 (SD = 626) in spring, so actual distances traveled were about 24% further than the minimum possible distance. The longest apparently continuous flight, recorded for Ifigeneia, was 3530 km, covered in 85 h with an average migration speed of 40 km/h (maximum single segment migration speed 73 km/h). We found that during the crossing of the drier (northern) part of their migration routes, the falcons used areas with higher NDVI in lower elevations compared to generated random points within the 50 km buffer zone either side of their path (Table 2a). In the more vegetated (southern) part of their migration routes, tagged individuals did not specifically move over areas with higher NDVI, though they did select areas with denser tree cover, especially in spring, and at higher elevations particularly in autumn, compared to what was available to them along their migratory paths (Table 2b).
Birds mostly roosted at night, with an average diurnal and nocturnal migration speed of 20.02 km/h (N = 588, SD = 15.33) and 6.41 km/h (N = 474, SD = 15.94) respectively. The highest migration speeds though were recorded at night, both in autumn and spring, over the Mediterranean, Sinai and the Red Sea with a maximum migration speed of 89 km/h recorded. Overall, migration speeds were significantly slower over areas with higher NDVI and higher tree cover and also significantly slower in spring compared to autumn because of the relative effect of tailwinds (Table 3a). After excluding the transitioning segments between day and night, slower migration speeds were recorded at night and over areas with higher NDVI (Fig. 4a), with differences in migration speed between day and night being greater in autumn than in spring (Table 3b). After filtering out those segments with migration speeds < 5 km/h, we found that migration speeds were higher with increasing tailwind (Fig. 6). Migration speeds were still slower over higher NDVI, but conversely, higher during the night (N = 84, mean 35.61 km/h, SD = 19.98) compared with the day (N = 454, mean 25.42 km/h. SD = 13.25; Fig. 4a). There was a seasonal effect in this active migratory movement dataset, with the difference in migration speed between day and night being again higher in autumn compared to spring, while there was a negative interaction effect of day and night and tree cover, indicating that falcons travelled more quickly at night compared to the day over less vegetated areas (Table 3c).
In spite of the observed loop migration pattern, both average distances covered and migration duration were found to be similar between seasons, in contrast with findings for birds tagged in the western Mediterranean , though our results are based on a small sample size. On the other hand, departure dates from breeding grounds are in agreement with results from other colonies [24, 39, 41], while in spring the narrow migration window we found also coincides with that found for birds breeding in the western Mediterranean .
Despite the small sample size, we believe that the relatively small latitudinal temporal differences between consecutive trips and the relatively narrow migration corridors between consecutive routes, suggest some individual repeatability. This contradicts results for birds tagged in the western Mediterranean for which no individual repeatability was identified, with larger longitudinal distances between consecutive trips . Eleonora’s falcon overwinters in a comparatively small area, primarily in Madagascar, indicating that in addition to internal mechanisms of orientation , the species also uses external navigation cues . As soon as migrants establish an optimal route and survive, they remember and copy it for consecutive years , enhancing migration success .
The preference for vegetation-rich areas with active selection of higher NDVI over drier areas in the north and higher tree cover further south, in conjunction with the observed slower migration speeds over those vegetation-rich areas, indicate that birds spend time refueling in vegetation-rich areas, possibly with fly-and-forage strategy, thus reducing the necessity for regular stopovers . In addition, because of the high spatial accuracy and more consistent temporal position collection of this study with the use of GPS transmitters, we were able to show that in contrast to findings for birds breeding in western Mediterranean , birds roosted during most nights in vegetation-rich areas primarily in savannah in autumn and shrublands in spring, resulting in significantly slower migration speed at night. Despite that overall trend, when considering only active movements (migration speed > 5 km/h), migration speeds were significantly faster at night, when the birds would not be foraging when flying, with possibly better flying conditions as well , further supporting the daily fly-and-forage assumption. Even when we account for migration speed differences due to variation in wind support, we found that falcons travel at significantly higher migration speeds over vegetation-poor areas, day and night. This is in line with recent findings that upon facing the challenge of crossing an ecological barrier, birds are able to adapt their diurnal and nocturnal flight patterns to overcome it as quickly as possible . This is in accordance with expectations  and supports the hypothesis that vegetation-poor areas are ecological barriers for Eleonora’s falcon , which they overcome by rapidly traversing rather than seeking a detour [24, 39, 40], with higher migration speeds over sea and desert . This indicates that landscape characteristics are guiding the routes and the migration speed , with individuals behaving as ‘sprinting migrants’ when feeding opportunities are lacking, such as over the Sahara desert . Characteristically, the remarkable non-stop flight by one individual for 3530 km in 85 h, a first such record for Eleonora’s falcon and one of the longest trips recorded for any species , indicates the capacity of Eleonora’s falcon to travel extensive distances, over sea and other ecological barriers, such as deserts, without refueling. Lesser kestrels and Eurasian hobbies have been recorded flying up to 1626 km and 740 km respectively, probably non-stop, during migration [44, 84], but Eleonora’s falcon appears to be capable of much longer flights.
In autumn, the Sahel region south of the Sahara is quite green with an abundance of flood pools and insects , and birds benefit from refueling there after the Sahara crossing. On the contrary, in spring the Sahel is at its driest, with potentially less food available , possibly contributing to the observed loop pattern. The short rainy season in eastern Ethiopia in spring might facilitate movements through there because of a timely abundance of potential insect prey , as suggested similarly for Red-backed shrikes . This agrees with field records during spring migration, where Eleonora’s falcons have been observed hawking on large insects and perching on trees in Somalia following the rains . The occasional delay in crossing the Mozambique Channel in autumn could be attributed to unfavorable weather conditions, or alternatively an expected need to refuel before crossing over an ecological barrier , as found in other studies [39, 41]. The suggestion though that this stopover area towards Malawi can be considered as part of the wintering grounds of Eleonora’s falcon [36, 41] is supported by our results here and deserves further investigation.
We suggest that the loop migration pattern is most likely attributed to habitat availability and the need for refueling, in agreement with previous studies [39, 40], yet, although long distance migrants are constraint by distance and time, and are thus expected to be less selective to wind conditions compared to short distance migrants , the effect of wind cannot be overlooked. Northerly tailwinds appear to aid the southbound Sahara crossing [3, 80]. On the other hand in spring, prevailing winds might contribute to the eastwards shift towards Ethiopia and Somalia . From there tagged individuals cross the Red Sea at its narrowest point at Djibouti and remain east of the Red Sea, avoiding stronger headwinds west of the Red Sea (Fig. 2) [3, 80, 83]. Interestingly, this anticlockwise pattern in Eleonora’s falcon is contrary to the overall trend in many species, where autumn migratory routes are more eastern than spring ones in the Palearctic . Indeed, a clockwise migratory loop pattern has been found in Lesser kestrel , Marsh harrier , and Egyptian vulture , but this was attributed to dominant wind patterns between the two migration seasons in West Africa where wind patterns differ from our study area.
Multi-year monitoring provides us with a greater understanding of migration cycles  and our 4-year study is the longest telemetry study on Eleonora’s falcon thus far. Apart from duration, this study is the first to use GPS-GSM loggers on Eleonora’s falcon, providing opportunities for more detailed information on migratory routes , and allowing for more accurate explanations of bird migratory strategy in relation to wind, vegetative cover and feeding opportunities [13, 90]. More long-term migration studies with high accuracy transmitters providing flying altitude, will further contribute to our understanding of migratory strategies that in turn can better inform conservation efforts for long distance migratory birds during what is often the least well studied portion of their annual cycle [13, 91].
Availability of data and materials
All datasets used and/or analysed during the current study are available from the corresponding author on reasonable request.
Above sea level
Corrected Akaike’s Information Criterion
European Centre for Medium-Range Weather Forecasts
Generalized Linear Mixed Models
Geospatial Modelling Environment
Global Positioning System
Global System for Mobile communication
Moderate-resolution Imaging Spectroradiometer
Normalized Difference Vegetation Index
Platform Transmitter Terminals
Alerstam T. Bird migration. Cambridge: Cambridge University Press; 1993.
Sutherland WJ. Evidence for flexibility and constraint in migration systems. J Avian Biol. 1998;29:441–6.
Newton I. The migration ecology of birds. London: Academic Press; 2008.
Pulido F. The genetics and evolution of avian migration. AIBS Bull. 2007;57:165–74.
Alerstam T, Hedenström A, Åkesson S. Long-distance migration: evolution and determinants. Oikos. 2003;103:247–60.
Berthold P. Bird migration: a general survey. 2nd ed. Oxford: Oxford University Press; 2001.
Gwinner E, Wiltschko W. Endogenously controlled changes in migratory direction of the garden warbler, Sylvia borin. J Comp Physiol A. 1978;125:267–73.
Tøttrup AP, Pedersen L, Onrubia A, Klaassen HG, Thorup K. Migration of red-backed shrikes from the Iberial peninsula: optimal or sub-optimal detour? J Avian Biol. 2017;48:149–54.
Thorup K, Rabøl J. The orientation system and migration pattern of long-distance migrants: conflict between model predictions and observed patterns. J Avian Biol. 2001;32:111–9.
Nathan R, Getz WM, Revilla E, Holyoak M, Kadmon R, Saltz D, Smouse PE. A movement ecology paradigm for unifying organismal movement research. Proc Natl Acad Sci U S A. 2008;105:19052–9.
Kerlinger P. How birds migrate. Mechanicsburg: Stackpole Books; 1995.
Henshaw I, Fransson T, Jakobsson S, Kullberg C. Geomagnetic field affects spring migratory direction in a long distance migrant. Behav Ecol Sociobiol. 2010;64:1317–23.
Klaassen RH, Strandberg R, Hake M, Alerstam T. Flexibility in daily travel routines causes regional variation in bird migration speed. Behav Ecol Sociobiol. 2008;62:1427–32.
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.
Alerstam T. Optimal bird migration revisited. J Ornithol. 2011;152:5–23.
Åkesson S, Hedenström A. How migrants get there: migratory performance and orientation. AIBS Bull. 2007;57:123–33.
Garshelis DL. Delusions in habitat evaluation: measuring use, selection, and importance. In: Boitani L, Fuller TK, editors. Research techniques in animal ecology: controversies and consequences. New York: Columbia University Press; 2000. p. 111–64.
López-López P. Individual-based tracking systems in ornithology: welcome to the era of big data. Ardeola. 2016;63:103–36.
Nilsson C, Klaassen RH, Alerstam T. Differences in speed and duration of bird migration between spring and autumn. Am Nat. 2013;181:837–45.
Warnock N. Stopping vs. staging: the difference between a hop and a jump. J Avian Biol. 2019;41:621–6.
Piersma T. Hop, skip, or jump? Constraints on migration of arctic waders by feeding, fattening, and flight speed. Limosa. 1987;60:185–94.
Alerstam T. Flight by night or day? Optimal daily timing of bird migration. J Theor Biol. 2009;258:530–6.
Bildstein KL. Migrating raptors of the world. Their ecology and conservation. Ithaca: Cornell University Press; 2006.
López-López P, Limiñana R, Mellone U, Urios V. From the Mediterranean Sea to Madagascar: are there ecological barriers for the long-distance migrant Eleonora’s falcon? Landsc Ecol. 2010;25:803–13.
Hahn S, Emmenegger T, Lisovski S, Amrhein V, Zehtindjiev P, Liechti F. Variable detours in long-distance migration across ecological barriers and their relation to habitat availability at ground. Ecol Evol. 2014;4:4150–60.
Liechti F. Bird: blowin’ by the wind? J Ornithol. 2006;147:202–11.
Shamoun-Baranes J, Liechti F, Vansteelant WMG. Atmospheric conditions create freeways, detours and tailbacks for migrating birds. J Comp Physiol A. 2017;203:509–29.
Norevik G, Åkesson A, Hedenström A. Migration strategies and annual spece-use in an afro-Palearctic aerial insectivore – the European nightjar Caprimulgus europaeus. J Avian Biol. 2016;48:738–47.
Vansteelant WMG, Shamoun-Baranes J, van Manen W, van Diermen J, Bouten W. Seasonal detours by soaring migrants shaped by wind regimes along the East Atlantic flyway. J Anim Ecol. 2017;86:179–91.
Hedenström A, Åkesson A. Adaptive airspeed adjustment and compensation for wind drift in the common swift: differences between day and night. Anim Behav. 2017;127:117–23.
Strandberg R, Klaassen RH, Olofsson P, Alerstam T. Daily travel schedules of adult Eurasian hobbies Falco subbuteo – variability in flight hours and migration speed along the route. Ardea. 2009;97:287–95.
Meyburg B, Howey PW, Meyburg C, Dietrich FK. Two complete migration cycles of an adult hobby tracked by satellite. Br Birds. 2011;104:2–15.
Monti F, Grémillet D, Sforzi A, Dominici JM, Triay Bagur R, Muñoz Navarro A, Fusani L, Klaassen RHG, Alerstam T, Duriez O. Migration distance affects stopover use but not travel speed: contrasting patterns between long- and short-distance migrating ospreys. J Avian Biol. 2018. https://doi.org/10.1111/jav.01839.
Duriez O, Péron G, Grémillet D, Sforzi A, Monti F. Migrating ospreys use thermal uplifts over the open sea. Biol Lett. 2018. https://doi.org/10.1098/rsbl.2018.0687.
Walter H. Eleonora's falcon. Adaptations to prey and habitat in a social raptor. Chicago: The University of Chicago Press; 1979.
Kassara C, Gangoso L, Mellone U, Piasevoli G, Hadjikyriakou TG, Tsiopelas N, Giokas S, López-López P, Urios V, Figuerola J, Silva R, Bouten W, Kirschel ANG, Virani MZ, Fieedler W, Berthold P, Gschweng M. Current and future suitability of wintering grounds for a long-distance migratory raptor. Sci Rep. 2017;7:8798.
Orta J. Eleonora's falcon Falco eleonorae. In: del Hoyo J, Elliott A, Sargatal J, editors. Handbook of the birds of the world, vol. 2, new World vultures to Guinea fowl. Barcelona: Lynx Edicions; 1994. p. 266.
Cramp S, Simmons K, Gillmor R, Hollom P, Hudson R, Nicholson E, Ogilvie M, Olney P, Roselaar C, Voous K, Wallace DIM, Wattel J. Handbook of the birds of Europe the Middle East and North Africa. The birds of the western palearctic, vol. 2 hawks to bustards. Oxford: Oxford University Press; 1980.
Kassara C, Fric J, Gschweng M, Sfenthourakis S. Complementing the puzzle of Eleonora’s falcon (Falco eleonorae) migration: new evidence from an eastern colony in the Aegean Sea. J Ornithol. 2012;153:839–48.
Mellone U, López-López P, Limiñana R, Piasevoli G, Urios V. The trans-equatorial loop migration system of Eleonora’s falcon: differences in migration patterns between age classes, regions and seasons. J Avian Biol. 2013;44:417–26.
Gschweng M, Kalko EK, Querner U, Fiedler W, Berthold P. All across Africa: highly individual migration routes of Eleonora's falcon. Proc Biol Sci. 2008;275:2887–96.
Hadjikyriakou TG, Kirschel ANG. Video evidence confirms cannibalism in Eleonora’s falcon. J Raptor Res. 2016;50:220–3.
Bird DM, Bildstein KL, editors. Raptor research and management techniques. Surrey: Hancock House; 2007.
Limiñana R, Romero M, Mellone U, Urios V. Mapping the migratory routes and wintering areas of lesser kestrels Falco naumanni: new insights from satellite telemetry. Ibis. 2012;154:389–99.
Gosler A. Birds in the hand. In: Sutherland WJ, Newton I, Green R, editors. Bird ecology and conservation: a handbook of techniques. Oxford: Oxford University Press; 2004. p. 85–118.
Redfern C, Clark J. Ringers’ manual. Norfolk: British Trust for Ornithology; 2001.
Dawson A. Techniques in physiology and genetics. Oxford: Oxford University Press; 2004.
Forsman D. The raptors of Europe and the Middle East. London: Helm; 1999.
Kenward RE. A manual for wildlife radio tagging. London: Academic Press; 2001.
Fair JM, Paul E, Jones J, editors. Guidelines to the use of wild birds in research, 3rd ed. Washington: The Ornithological Council; 2010.
Bán M, Moskát C, Fülöp A, Hauber ME. Return migration of common cuckoos (Cuculus canorus) between breeding grounds in Hungary and wintering grounds in Africa as documented by non-PTT GPS technology. J Ornithol. 2017;159:337–44.
Barron DG, Brawn JD, Weatherhead PJ. Meta-analysis of transmitter effects on avian behaviour and ecology. Methods Ecol Evol. 2010;1:180–7.
Telemetry M. 5g solar PTT 100–5 field manual. Columbia: Microwave Telemetry Inc; 2013.
Ecotone. GPS GSM loggers. Gdynia: Ecotone Telemetry; 2013.
CLS. Argos User’s manual. Toulouse: Collecte Localisation Satellites; 2013.
Thorup K, Alerstam T, Hake M, Kjellén N. Can vector summation describe the orientation system of juvenile ospreys and honey buzzards? An analysis of ring recoveries and satellite tracking. Oikos. 2003;103:350–9.
Hays G, Åkesson S, Godley B, Luschi P, Santidrian P. The implications of location accuracy for the interpretation of satellite-tracking data. Anim Behav. 2001;61:1035–40.
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.
Limiñana R, Soutullo A, López-López P, Urios V. Pre-migratory movements of adult Montagu’s harriers Circus pygargus. Ardea. 2008;96:81–90.
Kassara C, Fric J, Sfenthourakis S. Distribution modeling of Eleonora’s falcon Falco eleonorae Géné, 1839 occurrence in its wintering grounds: a niche-based approach with satellite telemetry data. Bird Conserv Int. 2014;24:100–13.
Jorge PE, Sowter D, Marques PA. Differential annual movement patterns in a migratory species: effects of experience and sexual maturation. PLoS One. 2011;6:e22433.
ESRI. ArcGIS, computer program, ver. 10.1. Redlands: ESRI Inc; 2012.
Morrissay D. Sunsetsunrisemap. 2017. https://sunrisesunsetmap.com/ . Accessed 16 Sept 2017.
Desholm M. How much do small-scale changes in flight direction increase overall migration distance? J Avian Biol. 2003;34:155–8.
van Wijk RE, Kölzsch A, Kruckenberg H, Ebbinge BS, Müskens GJ, Nolet BA. Individually tracked geese follow peaks of temperature acceleration during spring migration. Oikos. 2012;121:655–64.
NASA. Reverb / ECHO Earth data. 2018. https://reverb.echo.nasa.gov/reverb. Accessed 8 Feb 2018.
Google Inc. Google Earth Pro, computer program. Mountain View: Google Inc; 2013.
Klaassen RH, Strandberg R, Hake M, Olofsson P, Tøttrup AP, Alerstam T. Loop migration in adult marsh harriers Circus aeruginosus, as revealed by satellite telemetry. J Avian Biol. 2010;41:200–7.
Sutherland WJ, Newton I, Green R. Bird ecology and conservation: a handbook of techniques. Oxford: Oxford University Press; 2004.
Alerstam T, Hake M, Kjellén N. Temporal and spatial patterns of repeated migratory journeys by ospreys. Anim Behav. 2006;71:555–66.
Beyer HL. Geospatial Modeling Environment, computer program, ver. 7.2.1: Spatial Ecology LLC; 2012.
Thorup K, Tøttrup AP, Willemoes M, Klaassen RHG, Strandberg R, Lomas Vega M, Dasari HP, Araújo MB, Wikelski M, Rahbek C. Resource tracking within and across continents in long-distance migrants. Sci Adv. 2017;3:e1601360.
Anyamba A, Tucker CJ. Analysis of sahelian vegetation dynamics using NOAA-AVHRR NDVI data from 1981–2003. J Arid Environ. 2005;63:596–614.
Trierweiler C, Mullie WC, Drent RH, Exo K, Komdeur J, Bairlein F, Harouna A, Bakker M, Koks BJ. A Palaearctic migratory raptor species tracks shifting prey availability within its wintering range in the Sahel. J Anim Ecol. 2013;82:107–20.
Hadjikyriakou TG, Kassara C, René de Roland LA, Giokas S, Tsiopelas N, Evangelidis A, Thorstrom R, Kirschel ANG. Phenology, variation in habitat use, and daily activity patterns of Eleonora's falcon overwintering in Madagascar. Landsc Ecol. 2019. https://doi.org/10.1007/s10980-019-00940-6.
USGS. Global Multi-resolution Terrain Data. 2016. https://topotools.cr.usgs.gov/gmted_viewer/viewer.htm. Accessed 2 Feb 2018.
R Core Team. R: a language and environment for statistical computing, computer program, ver. 3.5.1. Austria: R Foundation for Statistical Computing; 2018.
Schielzeth H. Simple means to improve the interpretability of regression coefficients. Methods Ecol Evol. 2010;1:103–13.
Mellone U, Klaassen RHG, García-Ripollés C, Limiñana R, López-López P, Pavón D, Strandberg R, Urios V, Vardakis M, Alerstam T. Interspecific comparison of the performance of soaring migrants in relation to morphology, meteorological conditions and migration strategies. PLoS One. 2012;7:e39833.
Copernicus Climate Change Service. ERA5: Fifth generation of ECMWF atmospheric reanalyses of the global climate. In: Copernicus Climate Change Service Climate Data Store; 2017. https://cds.climate.copernicus.eu/cdsapp#!/home. Accessed 12 May 2019.
Kemp MU, Emiel van Loon E, Shamoun-Baranes J, Bouten W. RNCEP: global weather and climate data at your fingertips. Methods Ecol Evol. 2012;3:65–70.
Adamík P, Emmenegger T, Briedis M, Gustafsson L, Henshaw I, Krist M, Laaksonen T, Liechti F, Procházka P, Salewski V. Barrier crossing in small avian migrants: individual tracking reveals prolonged nocturnal flights into the day as a common migratory strategy. Sci Rep. 2016;6:21560.
Tøttrup AP, Klaassen RH, Strandberg R, Thorup K, Kristensen MW, Jorgensen PS, Fox J, Afanasyev V, Rahbek C, Alerstam T. The annual cycle of a trans-equatorial Eurasian-African passerine migrant: different spatio-temporal strategies for autumn and spring migration. Proc Biol Sci. 2012;279:1008–16.
Strandberg R, Klaassen RH, Hake M, Olofsson P, Alerstam T. Converging migration routes of Eurasian hobbies Falco subbuteo crossing the African equatorial rain forest. Proc Biol Sci. 2009;276:727–33.
Alerstam T. Reoriented bird migration in coastal areas: dispersal to suitable resting grounds? Oikos. 1978;20:405–8.
Kranstauber B, Weinzierl R, Wikelski M, Safi K. Global aerial flyways allow efficient travelling. Ecol Lett. 2015;18:1338–45.
Strandberg R, Klaassen RH, Hake M, Olofsson P, Thorup K, Alerstam T. Complex timing of marsh harrier Circus aeruginosus migration due to pre-and post-migratory movements. Ardea. 2008;96:159–71.
García-Ripollés C, López-López P, Urios V. First description of migration and wintering of adult Egyptian vultures Neophron percnopterus tracked by GPS satellite telemetry. Bird Study. 2010;57:261–5.
Tomkiewicz SM, Fuller MR, Kie JG, Bates KK. Global positioning system and associated technologies in animal behaviour and ecological research. Philos trans R. Soc Lond B Biol Sci. 2010;365:2163–76.
Guilford T, Åkesson S, Gagliardo A, Holland RA, Mouritsen H, Muheim R, Wiltschko R, Wiltschko W, Bingman VP. Migratory navigation in birds: new opportunities in an era of fast-developing tracking technology. J Exp Biol. 2011;214:3705–12.
Hallworth MT, Marra PP. Miniaturized GPS tags identify non-breeding territories of a small breeding migratory songbird. Sci Rep. 2015;5:11069.
We are very grateful to C. Kassara, for her valuable help throughout this study, and for her valuable comments on draft versions of this manuscript. We would like to thank W. Cresswell. W. Vansteelant and four anonymous reviewers for helpful comments on previous versions of the manuscript. We also would like to thank P. Charilaou, SBAA Environment and Conservation Officer for his support to the project. We are also especially grateful to the late A. Crabtree for his invaluable guidance in the field.
This work was supported by the A. G. Leventis Foundation; The Peregrine Fund; the Joint Nature Conservation Committee and the Sovereign Base Areas Administration (SBAA).
Ethics approval and consent to participate
A bird ringing licence and a research licence to attach bird transmitters were obtained from the Sovereign Bases Areas Administration (SBAA) Environment Department.
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.
GLMM results for comparison between actual positions and random points generated within a 50 km buffer zone either side of the flying path, in relation to NDVI, tree cover and elevation for the north part of the migratory routes (latitude 16 N northwards) (N = 901): a) full list of candidate models; b) full model; c) candidate model (ΔAICc <4); d) candidate model (ΔAICc <4); e) candidate model (ΔAICc <4); f) candidate model (ΔAICc <4); g) best supported model. In the drier north (e.g. Sahara desert), the falcons used areas with higher NDVI compared to randomly generated points, avoiding higher elevations. Table S2. GLMM results for comparison between actual positions and random points generated within a 50 km buffer zone either side of the flying path, in relation to NDVI, tree cover and elevation for the south part of the migratory routes (latitude 16 N southwards) (N = 2811): a) full list of candidate models; b) full model, c) candidate model (Δ AICc < 4); d) best supported model. In the more vegetated south they occurred (especially in spring), in areas with higher tree cover, and subsequently higher elevations (particularly in autumn) than what was available to them alongside their routes. Table S3. GLMM results for the effect of all GPS telemetry locations (N = 1435) on migration speed, of NDVI, tailwind, tree cover and season, with interactions of season with the other fixed effects, and individual and year as crossed random factors: a) full list of candidate models; b) full model; c) Candidate model (ΔAICc < 4); d) Candidate model (ΔAICc < 4); e) Candidate model (ΔAICc < 4); f) best supported model. Migration speeds were significantly slower over areas with higher NDVI and tree cover and also significantly slower in spring compared to autumn because of wind effect. Table S4. GLMM results for the effect on migration speed, of NDVI, tailwind, tree cover, season and day/night (excluding transitioning segments between day and night) with the interactions of day/night with the other fixed effects and individual and year as crossed random factors (N = 984): a) full model; b) best supported model (Delta AICc < 4). In this dataset, slower migration speeds were recorded at night, and over areas with higher NDVI. Differences in migration speed between day and night were greater in autumn compared with spring, with significantly slower migration speeds at night during spring. Table S5. GLMM results for the effect on active migratory movements’ migration speed (> 5 km/h), of NDVI, tailwind, tree cover, season and day/night (excluding transitioning segments between day and night; N = 465), with interactions of day/night with the other fixed effects and individual and year as crossed random factors: a) full list of candidate models; b) full model; c) Candidate model (ΔAICc < 4); d) Candidate model (ΔAICc < 4); e) best supported model. In this dataset, migration speeds were higher with increasing tailwind. Migration speeds were still slower over higher NDVI, but conversely, higher during the night compared with the day. Also, the effect of increasing tree cover percentage on migration speed was higher during the day compared to the night. A seasonal effect was also identified here, with the difference in migration speed between day and night being again greater in autumn compared to spring with significantly slower migration speeds at night during spring. Figure S1. Migratory routes during (a) autumn 2013 and (b) spring 2014 migration seasons overlaid on to the respective average monthly east – west winds (U wind component) at 850 hPa (c1500 m a.s.l.). Line colors represent the different individuals. East-west wind component is weaker compared to north-south component, facilitating north-south primary movement direction to and from wintering grounds. Data downloaded from ERA5 atmospheric reanalysis product of the European Centre for Medium-Range Weather Forecasts (ECMWF).
About this article
Cite this article
Hadjikyriakou, T.G., Nwankwo, E.C., Virani, M.Z. et al. Habitat availability influences migration speed, refueling patterns and seasonal flyways of a fly-and-forage migrant. Mov Ecol 8, 10 (2020). https://doi.org/10.1186/s40462-020-0190-4