Flexible reaction norms to environmental variables along the migration route and the significance of stopover duration for total speed of migration in a songbird migrant
Frontiers in Zoology volume 14, Article number: 17 (2017)
Predicting the consequences of continuing anthropogenic changes in the environment for migratory behaviours such as phenology remains a major challenge. Predictions remain particularly difficult, because our knowledge is based on studies from single-snapshot observations at specific stopover sites along birds’ migration routes. However, a general understanding on how birds react to prevailing environmental conditions, e.g. their ‘phenotypic reaction norm’, throughout the annual cycle and along their entire migration routes is required to fully understand how migratory birds respond to rapid environmental change.
Here, we provide direct evidence that northern wheatears (Oenanthe oenanthe) from a breeding population in Alaska adjusted their probability to resume migration as well as the distance covered per night, i.e. travel speed, to large-scale environmental conditions experienced along their 15,000 km migratory route on both northwards and southwards migrations. These adjustments were found to be flexible in space and time. At the beginning of autumn migration, northern wheatears showed high departure probabilities and high travel speeds at low surface air temperatures, while far away from Alaska both traits decreased with increasing air temperatures. In spring, northern wheatears increasingly exploited flow assistance with season, which is likely a behavioural adjustment to speed up migration by increasing the distance travelled per night. Furthermore, the variation in total stopover duration but not in travel speed had a significant effect on the total speed of migration, indicating the prime importance of total stopover duration in the overall phenology of bird migration.
Northern wheatears from Alaska provide evidence that the phenotypic reaction norm to a set of environmental conditions cannot be generalized to universal and persistent behavioural reaction pattern across entire migratory pathways. This highlights the importance of full annual-cycle studies on migratory birds to better understand their response to the environment. Understanding the mechanisms behind phenotypic plasticity during migration is particularly important in the assessment of whether birds can keep pace with the potentially increasing phenological mismatches observed on the breeding grounds.
Human-induced climate change  has altered temporal patterns in the environment , resulting in a general advance in the arrival times of migratory birds to their breeding grounds [3–6]. However, the responses of species to rapidly changing environments are not synchronized which often leads to serious phenological mismatches between trophic levels [7–9]. And although many long-distance migratory bird species have significantly advanced their spring phenology [4, 10, 11], their current worldwide decline  indicates remaining or even increasing phenological mismatches with their primary food sources [13–17]; but see [18, 19]. Three non-mutually exclusive proximate causes can explain variation in arrival dates: start of migration, total migration distance, and total speed of migration, i.e. the average migratory distance covered per day including periods of stopover during the entire migratory period [km/day], [20, 21]. Furthermore, ontogenetic factors  as well as exogenous phenological events at the breeding area may influence the timing of the various life-history stages in migratory birds . Conditions at the wintering grounds may also be used as a cue to initiate spring migration [23, 24], which is a good predictor for breeding area arrival timing [25–29]. However, high within-individual repeatability in the timing of spring migration in some species [27, 30] suggest that conditions at the wintering grounds may not necessarily influence arrival dates [31, 32]. If migrants experience favourable environmental conditions along the entire migratory route, the arrival at breeding grounds can potentially be adjusted by increasing the total speed of migration [20, 33, 34]. Yet, this adjustment might be limited since the rate of accumulating energy at a stopover is much slower than the rate of spending energy in flight, indicating that stopovers are likely the crucial periods of the migratory journey affecting total speed of migration .
To date, our knowledge on the effects of the environment, e.g. prevailing weather conditions, on stopover duration and travel speed, i.e. the daily migratory distance covered during periods of flight, is primarily inferred from single-snapshot studies detailing the behaviour of individuals at a specific site [36, 37] or along a single migratory leg [38–40]. Yet, individuals may react differently to a similar set of environmental conditions  depending on e.g. the distance remaining to the destination . Information from one radar-tracking study from various species and two locations suggested that individuals tracked closed to the breeding area had a higher migration speed than those tracked further away . This suggests that the biological significance of environmental parameters may change during migration. In contrast to large species [44–51], it remains difficult to follow individual songbirds over large distances and entire annual cycles [52, 53]. Due to this constraint, we are still restricted in our ability to identify environmental parameters that determine the entire migratory schedule and whether stopover duration significantly affects the travel speed in long-distance migratory songbirds. In this study, we used light-level geolocators to track northern wheatears (Oenanthe oenanthe, wheatear hereafter), a small nocturnal migratory songbird, from a population breeding in Alaska. Wheatears are exclusive nocturnal migrants  and individuals breeding in Alaska are known to travel approximately 15,000 km from Alaska to their African wintering grounds . Our general aim was to investigate the extent to which birds adjust their departure timing and travel speed in response to local environmental conditions, and whether such decisions are influenced by the proximity of the breeding areas and/wintering grounds.
Our first objective was the quantification of the direct effects of environmental variation on daily departure probability and daily travel speed along the migration route. Using remotely sensed environmental parameters to describe weather changes and conditions along the migration route and within the regional stopover sites, we aim to model how each of these parameters influenced the daily probability of an individual to resuming migration, as well as the travel speed of each migration bout. We hypothesize that the daily departure probability increases with decreasing temperature in autumn [54, 56], as a mechanism to avoid high energy costs at low surface temperature during the night . Furthermore, in spring, we expect birds to slow down their migration or even reverse migration, if the temperature drops significantly due to high energy costs at cold conditions . The likelihood of this process may increase towards the breeding area in the Arctic and with even colder conditions. In addition, we predict that departure probability from stopover sites will increase with decreasing precipitation [36, 58, 59] and high air pressure [60–62]. Lower wind speeds [58, 59, 63] and favourable wind assistance [63, 64] are also known to increase the probability to resume migration and to also affect the individuals’ travel speed . Given the restrictions of nocturnal migrants to gain distance ‘solely’ during the night, early departures would allow birds to extent flight duration and potentially achieve larger distances per night. It is therefore assumed that favourable environmental conditions on the ground would results in early night departures, whereas less favourable conditions may lead to a higher variation in the nocturnal departure time . Finally, during spring migration we expect that prevailing climatic conditions would become less important to decision making the closer, the individual is to the breeding area: timely arrival at the breeding grounds is crucial , and individuals might be forced to resume migration, even if conditions are unfavourable once they get closer to the breeding area [42, 43]. Here we will provide evidence whether the phenotypic reaction norm to a set of environmental conditions is fixed or flexible across entire migratory pathways of northern wheatears in both seasons and in accordance with the above described predictions from single snap-shot studies. Our second objective was to quantify the direct effects of stopover duration and travel speed on total speed of migration and estimate the importance for arrival timing at the migratory destination. According to previous findings , we expect that seasonal variation in stopover duration of individual weathers and not mean travel speed determines their total speed of migration.
Study site and light-level geolocators
Wheatears were ringed and tagged at two sites in Alaska (Eagle Summit: 65.6° N, 145.4° W and Toolik: 68.6° N, 149.6° W, U.S.A.; Additional file 1). In 2009, 15 adult male and female wheatears were tagged with Mk10S geolocators from British Antarctic Survey (Cambridge, UK; 1.4 g representing max. 6% of bird’s body mass) at Eagle Summit . In 2013, 28 adult male, 32 adult female and 60 juvenile wheatears were tagged at Eagle Summit and 23 adult male and female and 74 juvenile wheatears were tagged at Toolik. In 2013, we used ML6190 geolocators from Biotrack (Wareham, UK; 0.9 g representing max. 4% of bird’s body mass) and Intigeo P65A9-11 from Migrate Technology (Cambridge, UK; 0.9 g). Adult body mass of tagged individuals was 26.0 ± 2.2 g (mean ± s.d., n = 135) compared to juveniles with a body mass of 25.9 ± 1.8 g (n = 134). Further information on the light-level geolocator devices as well as explanations on recapture probabilities can be found in Additional files 1 and 2.
Analyses of light-level geolocation data
Light-level geolocators archive light intensity over time. For any given site, light intensity changes specifically over the year with respect to a standard time, allowing estimation of positions [69, 70]. The uncertainty of those estimates is frequently described to be within 100–200 km, but is highly dependent on mainly the bird’s behaviour which affects the shading of the light sensor, and the time of the year . In addition, the accuracy and precision of the estimates depends on the analytical approach, see [72–74].
In this study, light-level geolocator data, available upon request from “Movebank” (www.movebank.org, project: Oenanthe oenanthe Alaska), were analysed using the software R 3.1.2 (Vienna, Austria)  with the freely available R packages “SGAT”  and “BAStag” , cf. [74, 78]. First, daily sunrise and sunset times were identified using the arbitrary light intensity of “1.35” for the BAS and Biotrack geolocators and 1.35 lux for the Intigeos, similar as in other studies [74, 78, 79] (Additional file 1). Ninety-three (0.02%) out of 4,606 sunrise and sunset times were discarded and 99 (0.02%) of the remaining ones were the difference was > 30 min from the otherwise stable surrounding twilight events were adjusted (Additional files 1 and 3). Initial locations were estimated using the simple threshold method with a zenith angle of 96.8 ± 1.3° (mean ± s.d.; n = 8). Next, a Bayesian approach was used to refine the simple threshold estimates and describe the posterior distribution of bird locations . Within a Markov Chain Monte Carlo simulation we incorporated a priori knowledge in terms of a species-specific land mask, a twilight model, and a species-specific movement model to obtain a likelihood surface for location estimates ; see  for further details. The land mask biased potential locations towards land, given that the wheatear is a typical land bird ; thus, the relative probability to be over land was higher than that to be over water. This mask was applied for daytime and night-time location estimates. The twilight model described the expected error in the detection of sunrise and sunset times. Since natural light cannot be detected before sunrise or after sunset, we used a log-normal density distribution with the parameters meanlog = 2.2 and sdlog = 1.0 . Parameters for the log-normal distribution were estimated using light intensity recordings from the breeding site (known location) pooled across all individuals. The zenith angle was then estimated for each individual separately. Songbirds are known to spend most of their time during migration at stopover sites [56, 57, 81]. The movement model has been parameterised to account for high probabilities of low speeds and lower probability of high speeds by assuming a gamma distribution with shape = 0.7 and scale = 0.05. The expected mean airspeed of wheatear’s is approximately 50 km/h . For further information about the movement model see Additional file 3. Parameters of these three models were used for each individual to describe the posterior distribution of estimated locations. Additional file 3 provides the corresponding R code.
Defining migration schedule
To discriminate between stopovers and periods of migratory flights, we used the “changeLight” function of the R package “GeoLight” . This method is based on a “changepoint” analysis that quantifies the probability of each sunrise and sunset to be different than the surrounding sunrise and sunset times and provided evidence for shifts/movement in the underlying locations. The sunrise and sunset times associated with a “changepoint” probability greater than the 0.75 quantile of all probabilities were used to define stationary periods separated by periods of movement. Using the 0.75 quantile of all “changepoint” probabilities and only considering stopovers of two days or longer is a conservative approach . If the mean of the location estimates of consecutive stationary periods fall within a radius of 200 km, we merged those considered sites to be one. However, merging of stationary periods occurred only once in two birds at the beginning of autumn migration. By definition, birds resumed migration on the last evening of a stopover. Because the selected zenith influenced latitude estimates, the onset of autumn migration was defined as the start of the period of consecutive migratory flights preceding the last stationary site at > 2° of longitude westward of the breeding area. For stopover sites closer than 2° longitude, the last day of the subsequent period of consecutive migratory flights was defined as the onset of autumn migration. Termination of autumn migration was defined by the beginning of the first stationary site in eastern Africa. The onset of spring migration was defined as the first date of the period of consecutive migratory flights preceding the first stationary period 200 km away from the wintering ground. In spring, the last day of the last period of consecutive migratory flights before reaching the breeding area was considered as the termination of spring migration. However, this date could be determined only for one wheatear because for the others, the light intensity was above the threshold of 1.35 for detecting twilight events (see above) when Alaska was reached. Total speed of migration [km/day] was estimated as the total migration distance [km], as tracked by light-level geolocation data, divided by the duration [days] of the entire migration. In spring, total speed of migration was simplified as the distance from bird’s wintering ground to its last location estimate divided by the duration for this subset of the journey. The travel speed [km/days] was defined as the great circle distance [km] covered between daily locations on travel days [day], excluding days of stopover.
National Centres for Environmental Prediction Reanalysis II data from the National Oceanic and Atmospheric Administration (NOAA, Boulder, CO, USA; ) were obtained via the R package ”RNCEP” . This dataset consists of various climatic measurements with a 2.5x2.5° spatial resolution corresponding to ground surface area cells of 120 × 120 km2 and 290 × 290 km2 for the study area. This area size reflects the general uncertainty of location estimates derived from light-level geolocation; resolution of the weather data and the accuracy and precision of the light-level geolocator data were therefore in the same order of magnitude (see above). Several studies have shown that RNCEP wind data correlated highly significantly with real wind measurements in terms of speed and direction [86–89]. To describe the general environmental conditions that has most likely been experienced by the birds, we downloaded daily surface air temperatures [°C], surface air pressure levels [hPa], and surface wind speed [m/s] for the time of local sunsets for each day of the migration period [84, 85]. The amount of precipitation [kg/m2/s] was summed from three hours before until three hours after local sunset. Environmental values were interpolated in space and time . Wheatears perform nocturnal exploratory flights  possibly to assess wind conditions aloft . The best flow assistance of four different pressure levels was considered (1000, 925, 850, and 700 mbar) [65, 87, 90, 91]. Flow assistance was quantified using equationAirspeed:
with θ is the angular difference between the direction into which the wind is blowing and the bird’s preferred direction of migratory movement, z is bird’s speed [m/s] relative to the wind, and y is wind speed [m/s], following , by applying the “NCEP.flight” function of the R package ”RNCEP” , with “flow.assist” set to “NCEP.Airspeed”, “air speed” set to 13 m/s following Bruderer and Boldt , and “path” to “great circle”. Thus, the best flow assistance of the four pressure levels to the next migratory goal was considered for each day, cf. . On 48 of the 1038 “bird-days” during both autumn and spring migration, equationAirspeed did not provide an estimate because the wind speed was higher than the assumed airspeed. We assumed that birds perceived these modelled environmental conditions en route and considered these for their daily departure decision. For modelling travel speed the same environment conditions were incorporated. However, flow assistance was modelled for each hour and corresponding wind conditions were interpolated in space and time and separately for the different pressure levels. Thus, birds did not experience an abrupt change in wind conditions when crossing two grid cells.
Statistics were calculated using the R 3.1.2 statistical software package . All dates and times are given in Greenwich Mean Time.
Numeric explanatory variables were z-transformed to achieve homogeneity of variance. All numeric explanatory variables considered in a model were tested against one another for collinearity with the “vif” function of the R package “usdm” . If collinearity exceeded the conservative threshold of 2, one of the corresponding covariates was removed. For values that were lower than 2, the variables were treated as not collinear . Autumn and spring migrations were modelled separately.
The variation in the individual departure probability along the migration route was modelled using a mixed effect logistic regression model (MELR) . Variation in the individual departure probability was related to the environmental conditions at sunset (surface air temperature, precipitation, surface air pressure, surface wind speed, flow assistance, and remaining migratory distance, for details see “Environmental data”). Year of deployment was included as a random factor. Scale parameters ranged between 0.75 and 1.4 (1.06 in the autumn model and 0.98 in the spring model) indicating that overdispersion was not a confounding factor . Furthermore, corresponding goodness of fit plots showed no violation of model assumptions (Additional file 4).
The variation in the travel speed was modelled using a linear mixed effect model (LMM) assuming normally distributed errors [95, 97]. To avoid violation of model assumptions stopover duration was log10-transformed . Surface air temperature, precipitation, surface air pressure, surface wind speed (all sunset values), flow assistance, and remaining migration distance were used as explanatory variables. In this model, the flow assistance was taken as the mean of the hourly flow assistance experienced by the bird towards the next estimated location. Year was included as a random effect. To account for the inherent uncertainty of single location estimates derived from light-level geolocation , the data values were weighted by the inverse of the standard deviation of the latitude estimate.
The surface air temperature was collinear with the remaining migration distance in all spring models. Hence, the potential effect of the remaining migration distance on the variation in departure probability and travel speed could not be separated from the potential effect of surface air temperature in the spring models. As a result, surface air temperature was removed from all spring models.
The variation in the total speed of migration was modelled by total stopover duration and season using a LMM and assuming normally distributed errors. The variation in the total speed of migration was also separately modelled by the mean travel speed and season again using a LMM [95, 97]. For both models, individual was included as a random factor and all numeric parameters, i.e. the response and the predictor variables, were log10-transformed to avoid violation of model assumptions  which we assessed graphically (Additional file 4).
We used improper prior distributions, namely p(β) ~ 1 for the coefficients, and p(β) ~ 1/σ for the variance parameters following , which were applied to all models described above. To obtain the posterior distribution, we directly simulated 2,000 values from the joint posterior distribution of the model parameters using the function sim of the package “arm” . The median of the simulated values from the joint posterior distributions of the model parameters was used as an estimate, with the 2.5 and 97.5% quantiles as the lower and upper limits of the 95% credible intervals (CrI), respectively. We declared an effect significant when the corresponding 95% CrI did not include zero. Statistics were detailed for significant interactions only but not for their corresponding primary effects. The mean value and corresponding 95% confidence intervals for the interactions of mixed models were estimated using the “effect” function of the R package “effects” . The values of other explanatory variables were fixed at their means. For further information see . Details on statistical models, including residual analyses that show no violation of model assumptions, are provided in Additional file 4.
Of the 13 tagged birds that return to their breeding sites in the year after deployment (5 in 2010; 8 in 2014), only eight tracks were obtained because two birds lost their geolocator device, and three individuals could not be recaptured (Additional file 2). Wheatears initiated autumn migration between the 4th and 17th of Aug., crossed the Bering Strait and migrated westwards across Russia and Kazakhstan before heading south and westwards towards Africa in the Middle East (Fig. 1, Table 1, and Additional file 5). After travelling 14,800 ± 1,800 km (n = 8), birds arrived at their wintering grounds between the 3rd and 27th of Nov. In spring, birds started their migration between the 24th of Mar. and the 14th of Apr. (Table 1). Geolocation by light failed to estimate locations for birds experiencing 24-h daylight, which was the case for wheatears approaching their breeding area in spring (Fig. 1). Therefore, spring migration routes were not complete and individuals were tracked for 8,000 ± 1,900 km. Most stopover periods occurred along the western part of the route (Fig. 1b, c). The experienced environmental conditions along the route were similar between years, e.g. during autumn migration 2009 and 2013 and during spring 2010 and 2014 (Additional file 6).
In autumn, departure probability generally decreased with increasing surface air temperature, improving flow assistance, and decreasing remaining migration distance (Table 2). We found no significant effect of the main effects precipitation, surface air pressure, and surface wind speed on the individual departure probability (Table 2). The interaction between the surface air temperature and the remaining migration distance on the birds departure probability was found significant, indicating that the reaction to temperature changed with the remaining migration distance (Table 2; Fig. 2). During their first 8,000 km wheatears commonly experienced temperatures ranging from -5 °C to +20 °C (Additional file 7). Surface air temperature significantly affected departure probability with low temperatures being associated with high departure probabilities during the first two thirds of migration distance (Fig. 2). In contrast, relatively low temperatures were associated with low departure probabilities and high temperatures with high probabilities along the remaining one third of the migration distance (Fig. 2). Note that wheatears did not encounter temperatures above 20 °C on their first two thirds of migration and below 5 °C on their final third of migration. Corresponding model predictions should therefore be treated with caution. The interaction between the surface air pressure and the remaining migration distance was significant, suggesting that low surface air pressures were generally associated with high departure probabilities at the beginning of autumn migration, while far from Alaska departure probability increased with increasing air pressures (Table 2, Fig. 2). In spring, departure probability increased with rising surface air pressure, improving flow assistance, and decreasing remaining migration distance (Table 2, Fig. 3). No significant effects were found for precipitation, surface wind speed, or any other interactions (Table 2).
In autumn, travel speed increased with air pressure and was generally higher, if surface air temperature was low (Table 3). The effect of the surface air temperature on travel speed diminished, once the birds approached the wintering grounds (Fig. 4). In spring, precipitation had an apparent positive effect on the travel speed at the beginning of spring migration, whereas the reverse effect was observed during the remaining half of spring migration. However, occurrence of precipitation was strongly biased towards the second half of spring migration, which may have potentially influenced the interaction (Additional file 8). The significant interaction between the flow assistance and the remaining migration distance indicated that birds took advantage of flow assistance as they approached the breeding areas (Fig. 5). Additional file 8 shows the flow assistance experienced by wheatears.
Neither the timing of autumn departure nor the total migration distance were related to arrival timing at the wintering grounds (Wilcoxon signed rank sum test: S > 61, p > 0.5, n = 8). As wheatears experienced 24 h of daylight at locations close to the breeding areas in spring, arrival timing could not be related either to start of migration or to total migration distance. The total stopover duration had a significant negative effect on the total speed of migration (LMM: log10(total stopover duration), 95% CrI: -0.4 – -0.20 d, n = 16; in a more complex model, the additional fixed effect ‘season’ was not significant, 95% CrI: -0.17–0.01, and therefore omitted in the final model; Table 4, Fig. 6a, Additional file 4). Due to the truncated tracks in spring, we would not compare the total stopover duration between seasons (Fig. 1, Additional file 5). In autumn, individuals showed a median stopover duration of 2.7 days per 1,000 km (25% quantile: 2.5 days, 75% quantile: 3.7 days), compared to 2.3 days (1.2 days, 2.8 days) in spring (Wilcoxon signed-rank test comparing autumn with spring values: V = 28, p = 0.20, n = 8). We found no effect of mean travel speed on total speed of migration (LMM: log10(mean travel speed), 95% CrI: -0.50–0.99 km/day, n = 16; Table 4, Fig. 6b). However, in this model the effect of the season was significant (95% CrI: 0.02–0.22) with higher values during spring migration (Table 4).
In this study, wheatears were tracked with light-level geolocators from their breeding areas in Alaska to their wintering grounds in eastern Africa and back (Fig. 1, Additional file 5). Despite the relatively low spatial resolution of both the tracking [70, 71, 100] and the environmental data , as well as the limited sample size of eight individuals, our results show clear dependencies of migratory behaviour and perceived environmental conditions. Most importantly, we found that individual birds adjusted their departure decisions and travel speeds to the encountered environmental conditions along their 15,000 km migration routes (Figs. 2, 3, 4 and 5). These adjustments were flexible in space and time indicating that the individual phenotypic reactions to the same set of environmental conditions may change along the migration route. Hence, we cannot assume that reactions to changes in the environment at specific sites somewhere along the pathway are general behavioural patterns that can be extrapolated to the entire migration route. In addition, we showed that variation in total stopover duration is the main factor in determining total speed of migration in the tracked wheatears. Hence, total stopover duration crucially affected the individual’s phenology (Fig. 6). Detailed understanding on how migratory birds react to changes in the environment along their entire migration route and how this affects total speed of migration is essential to improve our ability to predict the consequences of continuing and rapid environmental changes on bird migration.
Wheatears in this study responded differently to prevailing environmental conditions in autumn and spring (Figs. 2, 3 and 4). In autumn, birds did not time their migratory flights in accordance with increasing flow assistance, which would maximize their travel speed . No association was found with wind speed and precipitation (Table 2), contrasting findings of former studies [58, 59, 64, 65, 101, 102]. One reason could be, that the weather conditions at a local scale were different to the modelled regional weather conditions used in this study, or that other extrinsic conditions, e.g. surface air temperature, or intrinsic factors, e.g. body condition, had a stronger effect on birds’ decision to resume migration than precipitation and/or wind [37, 103]. The latter has often been observed in cases where flow assistance was not significantly adverse [54, 58, 104–106]. In our study flow assistance was indeed mostly favourable (> -5 m/s; Additional file 9). At the beginning of autumn migration, the tracked wheatears resumed migration and travelled long distances during conditions characterised by low surface air temperatures and low surface air pressure (Figs. 2 and 4). We know that energy costs associated with thermoregulation during stopover are significantly minimized by avoiding cold temperature on the ground , flying long stretches towards likely warmer regions and possibly eluding fronts potentially carrying rain. During stopover diurnal movements, e.g. in search for food, were shown to be inversely associated with rain . Due to increased mortality and energy expenditure, is expected to be costly to migrating during rain events . Hence, wheatears seem to minimize energy expenditure at stopovers during autumn. Models on departure probability and travel speed include predictions for environmental conditions not experienced by wheatears en route, e.g. high temperatures at the beginning and low temperatures at the final part of autumn migration (Figs. 2 and 4). This change in the distribution of environmental predictors (e.g. air temperature) with remaining migration distance may affect the significant interactions. Nevertheless, neglecting these extreme conditions still showed that wheatears adjusted their departure decision and travel speed to a similar set of environmental conditions depending on the remaining migration distance (Figs. 2 and 4).
In spring, surface air pressure positively influenced departure probability (Table 2). This affect is most likely attributed to the fact that high surface air pressures generally represents favourable conditions for migration . Flow assistance influenced the departure probability in spring (Fig. 3), supporting the findings of multiple other studies [58, 59, 64, 65, 101, 102]. Both, minimization of energy or time can lead to such behavioural pattern . Interestingly, reducing stopovers (Fig. 1) and increasing departure probabilities along the migration route (Fig. 3) weakened the importance of stopover duration on the total speed of migration but simultaneously strengthened the significance of travel speed. Therefore, the increase in exploiting flow assistance with season (Fig. 5) was likely a behavioural adjustment to speed up migration, achieved by an increase in the distance travelled per night. Why flow assistance had the opposite effect on travel speed at the beginning of spring (Fig. 5) remains unknown. In addition to weather conditions, migrants are able to use cues and experience associated to the landscape of stopover sites to adjust migration phenology [108, 109]. The facts that we found a concentration of stopover sites just before the Arabian and Sahara deserts in autumn as well as in spring (Fig. 1) indicate that these sites may have been used to prepare and/or recover from crossings those ecological barriers. It remains to be investigated whether the few stopovers in the far-east of Russia during autumn migration indicate infrequent patches of suitable habitats for wheatears to refuel, as demonstrated for bluethroats (Lusinia sveccia) , or highly favourable feeding conditions allowing quick replenishment of fuel resources used during the previous migratory flight.
The correct timing of arrival at breeding areas is predicted to have generally stronger, more immediate fitness consequences than the timing of arrival at the wintering grounds . Hence, minimizing time rather than energy use is likely more beneficial in spring than during autumn migration [35, 43]. This trade-off between minimizing either time or energy and the resulting seasonally different selection pressures should lead to seasonally different reactions to the environment, as observed in the wheatears of this study (Figs. 2, 3 and 4).
Total stopover duration was a significant factor explaining total speed of migration (Fig. 6), which is in line with a recent tracking study of greater white-fronted geese (Anser albifrons) . In contrast to other migration studies in songbirds, in which the seasonal start of migration was correlated with the arrival date at the final destination [20, 25–29, 78], total stopover duration in the tracked wheatears was found more important in determining arrival dates. It seems possibly that the extremely long distance of the migration route and the resulting high number of stopovers (Fig. 1) moderates the effect of the onset of migration on arrival timing (Table 1) in comparison to shorter distance migrants that rest at fewer stopover sites [29, 111, 112]. The arrival timing at migratory destination, i.e. breeding area or wintering ground, is also explained by the environmental conditions experienced somewhere en route [10, 113]. Thus, environmental conditions affecting departure probability and travel speed (Figs. 2, 4 and 5) might have a stronger cumulative effect on the phenology in wheatears breeding in Alaska than in shorter distance migrants. In the former, the cumulative effect of the environment on the total stopover duration is, therefore, likely an important driver of the variation in arrival timing via total speed of migration, cf. Fig. 6.
In addition to the variation in the environment affecting departure probability and travel speed (Tables 2 and 3), the endogenously controlled changes in bird physiology throughout an annual cycle influences the overall rate of energy accumulation [114–116]. However, the within-season variation in this rate cannot be attributed to an endogenously controlled change in the motivation to refuel over the season ; therefore, variation in the environment is the most likely source for the variation in the rate of accumulating energy within a season. The innate migration programme also contributes in the regulation of departure probability, which at least in spring increases with the progression of the season  and with the birds getting closer to their breeding area (Fig. 3). Further support for the importance of the innate migration program on birds’ phenology comes from a high within-individual repeatability in the timing of onset of migration and arrival at the migratory destination across years [27, 30]. Therefore, evolutionary processes acting on onset of spring migration, total speed of migration, and/or total migration distance are likely required for birds to adapt to future environmental conditions to arrive on time [5, 6, 22, 118]. However, an adaptive evolutionary response to an environmental change is not immediate [6, 119], and the extent of delay is a function of the rate of contemporary microevolution and generation times, i.e., at least one year in birds . Because climate change is characterized by an increase in temperature combined with greater between-year variability , among other variables, microevolutionary responses are unlikely to explain alone the twice a year seasonal variation that occurs in bird phenology. Therefore, phenotypic flexibility is generally accepted as an important proximate cause underlying recent changes in the timing of migration [20, 118, 122].
To ultimately assess the predictions, that phenotypic flexibility is one of the major drivers for the current shifts in migratory phenology , we need to find strong spatial autocorrelations along the migratory routes allowing migrants en-route to predict future conditions on the breeding area and/or wintering ground . Furthermore, we need a deeper understanding on how birds respond to these changes on site. Our results suggests that wheatears potentially adjusted total speed of migration to changes in the environment via a plastic phenotypic response in the departure probability, affecting stopover duration, and the travel speed. However, endogenous constraints and a spatial-specific reaction norm (e.g. not taking advantage of tail wind conditions during the first half of spring migration; Fig. 5) may limit the magnitude of this reaction (e.g. departure probability; Table 2, Fig. 2). Thus, changes in a certain environmental parameter may potentially evoke a strong response at one site, but a weaker or even a different reaction at another site. To improve understanding how phenological change  is proximately caused [15, 20], an elaborate effort to study migrants’ reaction to environmental conditions along the entire migration route as well as at the breeding areas and wintering grounds is indispensable.
Our results are based on location estimates derived from light-level geolocation data recorded during autumn and spring migration of eight wheatears breeding in Alaska. The presented results needs to be treated cautiously given the inherent problems of accuracy and precision in light-level geolocation as well as the coarse resolution of the environmental data . However, and despite these issues we found that wheatears phenotypically adjusted their departure probability and travel speed according to environmental conditions experienced along the entire migration route. More specifically, individuals showed flexible adjustments over time and space and between the seasons. The consideration of this complex phenotypic flexibility in bird reaction norms to the environment and the importance of total stopover duration on total speed of migration are essential to understand the influence of climate change on migration in general and timing in particular. Such knowledge might be crucial to understand why there is an apparent phenological mismatch in some [13–17] but not all species [18, 19] and their food resources, e.g. if the importance of total stopover duration on total speed of migration varies between species. In some species, total speed of migration might be unrelated or only weakly associated to variation in total stopover duration . In this case, and because environmental variation experienced during migration is unlikely to alter phenology via a change in total speed of migration, changes at the wintering grounds have likely the highest significant impact on their arrival times at the breeding area. These changes concern location of wintering grounds and onset of spring migration. In other species, total speed of migration is a function of total stopover duration (Fig. 6). In those, environmental conditions can have a significant effect on departure probability from stopover sites, which in turn shapes arrival timing at the migratory destination via total stopover duration. There are two non-mutually exclusive, potential solutions of how birds could overcome the current phenological mismatch with the primary food source at the breeding area. Either the phenotypic responses to ongoing changes in the environment are sufficient or evolutionary processes are required for migration phenology to keep pace. To predict the natural solutions to this challenge, future models need to incorporate the spatiotemporal-specific reaction norm with the environment and the species-specific significance of total stopover duration on total speed of migration.
Karl TR, Trenberth KE. Modern global climate change. Science. 2003;302:1719–23.
Thackeray SJ, Henrys PA, Hemming D, Bell JR, Botham MS, Burthe S, Helaouet P, Johns DG, Jones ID, Leech DI, et al. Phenological sensitivity to climate across taxa and trophic levels. Nature. 2016;535:241–5.
Walther G-R, Post E, Menzel A, Parmesan C, Beebee TJC, Fromentin J-M, Hoegh-Guldberg O, Bairlein F. Ecological responses to recent climate change. Nature. 2002;416:389–95.
Jenni L, Kéry M. Timing of autumn bird migration under climate change: advances in long-distance migrants, delays in short-distance migrants. Proc R Soc Lond B. 2003;270:1467–71.
Gill JA, Alves JA, Sutherland WJ, Appelton GF, Potts PM, Gunnarsson TG. Why is timing of bird migration advancing when individuals are not? Proc R Soc Lond B. 2014;281:20132161.
Visser ME, Gienapp P, Husby A, Morrisey M, de la Hera I, Pulido F, Both C. Effects of spring temperatures on the strength of selection on timing of reproduction in a long-distance migratory bird. PLoS Biol. 2015;13:e1002120.
Edwards M, Richardson AJ. Impact of climate change on marine pelagic phenology and trophic mismatch. Nature. 2004;430:881–4.
Both C, van Asch M, Bijlsma RG, van den Burg AB, Visser ME. Climate change and unequal phenological changes across four trophic levels: constraints or adaptations? J Anim Ecol. 2009;78:73–83.
van Gils JA, Lisovski S, Lok T, Meissner W, Ożarowska A, de Fouw J, Rakhimberdiev E, Soloviev MY, Piersma T, Klaassen M. Body shrinkage due to Arctic warming reduces red knot fitness in tropical wintering range. Science. 2016;352:819–21.
Hüppop O, Hüppop K. North Atlantic Oscillation and timing of spring migration in birds. Proc R Soc Lond B. 2003;270:233–40.
Jonzén N, Lindén A, Ergon T, Knudsen E, Vik JO, Rubolini D, Piacentini D, Brinch C, Spina F, Karlsson L, et al. Rapid advance of spring arrival dates in long-distance migratory birds. Science. 2006;312:1959–61.
Bairlein F. Migratory birds under threat. Science. 2016;354:547–8.
Møller AP, Rubolini D, Lehikoinen E. Populations of migratory bird species that did not show a phenological response to climate change are declining. PNAS. 2008;105:16195–200.
Both C, Van Turnhout CAM, Bijlsma RG, Siepel H, Van Strien AJ, Foppen RPB. Avian population consequences of climate change are most severe for long-distance migrants in seasonal habitats. Proc R Soc Lond B. 2010;277:1259–66.
Vickery JA, Ewing SR, Smith KW, Pain DJ, Bairlein F, Škorpiloá J, Gregory RD. The decline of Afro-Palaearctic migrants and an assessment of potential causes. Ibis. 2014;156:1–22.
Sanderson FJ, Donald PF, Pain DJ, Burfield IJ, van Bommel FPJ. Long-term population declines in Afro-Palearctic migrant birds. Biol Cons. 2006;131:93–105.
Rushing CS, Ryder TB, Marra PP. Quantifying drivers of population dynamics for a migratory bird throughout the annual cycle. Proc R Soc Lond B. 2016; 283.10.1098/rspb.2015.2846.
Dunn PO, Møller AP. Changes in breeding phenology and population size of birds. J Anim Ecol. 2014;83:729–39.
Senner NR, Stager M, Sandercock BK. Ecological mismatches are moderated by local conditions for two populations of a long-distance migratory bird. Oikos. 2017;126:61–72.
Knudsen E, Lindén A, Both C, Jonzén N, Pulido F, Saino N, Sutherland WJ, Bach LA, Coppack T, Ergon T, et al. Challenging claims in the study of migratory birds and climate change. Biol Rev. 2011;86:928–46.
Coppack T, Both C. Predicting life-cycle adaptations of migratory birds to global climate change. Ardea. 2002;90:369–78.
Durant JM, Hjermann DO, Otterseon G, Stenseth NC. Climate and the match or mismatch between predator requirements and resource availability. Clim Res. 2007;33:271–83.
Studds CE, Marra PP. Rainfall-induced changes in food availability modify the spring departure programme of a migratory bird. Proc R Soc Lond B. 2011;278:3437–43.
Cooper NW, Sherry TW, Marra PP. Experimental reduction of winter food decreases body condition and delays migration in a long-distance migratory bird. Ecol. 2015;96:1933–42.
Schmaljohann H, Meier C, Arlt D, Bairlein F, van Oosten HH, Morbey YE, Åkesson S, Buchmann M, Chernetsov N, Desaever R, et al. Proximate causes of avian protandry differ between subspecies with contrasting migration challenges. Behav Ecol. 2016;27:321–31.
Tøttrup AP, Klaassen RHG, Strandberg R, Thorup K, Willemoes Kristensen M, Søgaard Jørgensen P, Fox JW, 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 R Soc Lond B. 2012;279:1008–16.
Stanley CQ, MacPherson M, Fraser KC, McKinnon EA, Stutchbury BJM. Repeat tracking of individual songbirds reveals consistent migration timing but flexibility in route. PLoSONE. 2012;7:e40688.
Jahn AE, Cuteo VR, Fox JW, Husak MS, Kim DH, Landoll DV, Ledezma JP, LePage HK, Levey DJ, Murphy MT, Renfrew RB. Migration timing and wintering areas of three species of flycatchers (Tyrannus) breeding in the Great Plains of North America. Auk. 2013;130:247–57.
Ouwehand J, Both C. African departure rather than migration speed determines variation in spring arrival in pied flycatchers. J Anim Ecol. 2017;86:88–97.
Conklin JR, Battley PF, Potter MA. Absolute consistency: Individual versus population variation in annual-cycle schedules of a long-distance migrant bird. PLoSONE. 2013;8:e54535.
McKinnon EA, Stanley CQ, Stutchbury BJM. Carry-over effects of nonbreeding habitat on start-to-finish spring migration performance of a songbird. PLoSONE. 2015;10:e0141580.
Senner NR, Hochachka WM, Fox JW, Afanasyev V. An exception to the rule: carry-over effects do not accumulate in a long-distance migratory bird. PLoSONE. 2014;9:e86588.
Both C. Flexibility of timing of avian migration to climate change masked by environmental constraints en route. Curr Biol. 2010;20:1–6.
Marra PP, Francis CM, Mulvihill RS, Moore FR. The influence of climate on the timing and rate of spring bird migration. Oecol. 2005;142:307–15.
Alerstam T, Lindström Å. Optimal bird migration: the relative importance of time, energy, and safety. In: Gwinner E, editor. Bird migration: physiology and ecophysiology. Berlin Heidelberg: Springer; 1990. p. 331–51.
Smith AD, McWilliams SR. What to do when stopping over: behavioral decisions of a migrating songbird during stopover are dictated by initial change in their body condition and mediated by key environmental conditions. Behav Ecol. 2014;25:1423–35.
Jenni L, Schaub M. Behavioural and physiological reactions to environmental variation in bird migration: a review. In: Berthold P, Gwinner E, Sonnenschein E, editors. Avian migration. Berlin Heidelberg: Springer; 2003. p. 155–71.
Deppe JL, Ward MP, Bolus RT, Diehl RH, Celis-Murillo A, Zenzal TJ, Moore FR, Benson TJ, Smolinsky JA, Schofield LN, et al. Fat, weather, and date affect migratory songbirds’ departure decisions, routes, and time it takes to cross the Gulf of Mexico. PNAS. 2015;112:E6331–8.
Crysler ZJ, Ronconi RA, Taylor PD. Differential fall migratory routes of adult and juvenile Ipswich Sparrows (Passerculus sandwichensis princeps). Mov Ecol. 2016;4:3.
Mitchell GW, Woodworth BK, Taylor PD, Norris DR. Automated telemetry reveals age specific differences in flight duration and speed are driven by wind conditions in a migratory songbird. Mov Ecol. 2015;3:19.
Liechti F, Bruderer B. The relevance of wind for optimal migration theory. J Avian Biol. 1998;29:561–8.
Alerstam T. Detours in bird migration. J Theor Biol. 2001;209:319–31.
Karlsson H, Nilsson C, Bäckman J, Alerstam T. Nocturnal passerine migrants fly faster in spring than in autumn: a test of the time minimization hypothesis. Anim Behav. 2012;83:87–93.
Sergio F, Tanferna A, De Stephanis R, Jiménez LJ, Blas J, Tavecchia G, Preatoni D, Hiraldo F. Individual improvements and selective mortality shape lifelong migratory performance. Nature. 2014;515:410–3.
Chevallier D, Handrich Y, Georges J-Y, Baillon F, Brossault P, Aurouet A, Le Maho Y, Massemin S. Influence of weather conditions on the flight of migrating black storks. Proc R Soc Lond B. 2010;277:2755–64.
Klaassen RHG, Hake M, Strandberg R, Alerstam T. Geographical and temporal flexibility in the response to crosswinds by migrating raptors. Proc R Soc Lond B. 2010;278:1339–46.
Flack A, Fiedler W, Blas J, Pokrovsky I, Kaatz M, Mitropolsky M, Aghababyan K, Fakriadis I, Makrigianni E, Jerzak L, et al. Costs of migratory decisions: a comparison across eight white stork populations. Sci Adv. 2016;2:e1500931.
Mandel JT, Bildstein KL, Bohrer G, Winkler DW. Movement ecology of migration in turkey vultures. PNAS. 2008;105:19102–7.
Vansteelant WMG, Bouten W, Klaassen RHG, Koks BJ, Schlaich AE, van Diermen J, Van Loon EE, Shamoun-Baranes J. Regional and seasonal fl ight speeds of soaring migrants and the role of weather conditions at hourly and daily scales. J Avian Biol. 2015;46:25–39.
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. 2016;86:179–91.
Åkesson S, Bianco G, Hedenström A. Negotiating an ecological barrier: crossing the Sahara in relation to winds by common swifts. Philos T Roy Soc B. 2016;371:20150393.
Wikelski M, Kays R, Kasdin NJ, Thorup K, Smith JA, Swenson Jr GW. Going wild: what a global small-animal tracking system could do for experimental biologists. J Exp Biol. 2007;210:181–6.
McKinnon EA, Fraser KC, Stutchbury BJM. New discoveries in landbird migration using geolocators, and a flight plan for the future. Auk. 2013;130:211–22.
Schmaljohann H, Korner-Nievergelt F, Naef-Daenzer B, Nagel R, Maggini I, Bulte M, Bairlein F. Stopover optimization in a long-distance migrant: the role of fuel load and nocturnal take-off time in Alaskan northern wheatears (Oenanthe oenanthe). Front Zool. 2013;10:26.
Bairlein F, Norris DR, Nagel R, Bulte M, Voigt CC, Fox JW, Hussell DJT, Schmaljohann H. Cross-hemisphere migration of a 25-gram songbird. Biol Lett. 2012;8:505–7.
Schmaljohann H, Fox JW, Bairlein F. Phenotypic response to environmental cues, orientation and migration costs in songbirds flying halfway around the world. Anim Behav. 2012;84:623–40.
Wikelski M, Tarlow EM, Raim A, Diehl RH, Larkin RP, Visser GH. Costs of migration in free-flying songbirds. Nature. 2003;423:704.
Erni B, Liechti F, Underhill LG, Bruderer B. Wind and rain govern the intensity of nocturnal bird migration in central Europe - a log-linear regression analysis. Ardea. 2002;90:155–66.
Schaub M, Liechti F, Jenni L. Departure of migrating European robins, Erithacus rubecula, from a stopover site in relation to wind and rain. Anim Behav. 2004;229–237.
Dau CP. The fall migration of pacific flyway Brent Branta bernicla in relation to climatic conditions. Wildfowl. 1992;43:80–95.
Richardson WJ. Timing of bird migration in relation to weather: updated review. In: Gwinner E, editor. Bird migration. Berlin Heidelberg: Springer; 1990. p. 78–101.
Zehnder S, Åkesson S, Liechti F, Bruderer B. Nocturnal autumn bird migration at Falsterbo, south Sweden. J Avian Biol. 2001;32:239–48.
Delingat J, Bairlein F, Hedenström A. Obligatory barrier crossing and adaptive fuel management in migratory birds: the case of the Atlantic crossing in Northern Wheatears (Oenanthe oenanthe). Behav Ecol Sociobiol. 2008;62:1069–78.
Åkesson S, Hedenström A. Wind selectivity of migratory flight departures in birds. Behav Ecol Sociobiol. 2000;47:140–4.
Liechti F. Birds: blowin’ by the wind? J Ornithol. 2006;147:202–11.
Müller F, Taylor PD, Sjöberg S, Muheim R, Tsvey A, Schmaljohann H. Towards a conceptual framework for explaining variation in the nocturnal departure time of songbird migrants. Mov Ecol. 2016;4:24.
Kokko H. Competition for early arrival in birds. J Anim Ecol. 1999;68:940–50.
Kölzsch A, Müskens GJDM, Kruckenberg H, Glazov P, Weinzierl R, Nolet BA, Wikelski M. Towards a new understanding of migration timing: slower spring than autumn migration in geese reflects different decision rules for stopover use and departure. Oikos. 2016;125:1496–507.
Wilson RP, Ducamp JJ, Rees G, Culik BM, Niekamp K. Estimation of location: global coverage using light intensity. In: Priede IM, Swift SM, editors. Wildlife telemetry: remote monitoring and tracking of animals. Chichester: Ellis Horward; 1992. p. 131–4.
Hill RD. Theory of geolocation by light levels. In: Boeuf L, Burney J, Laws RM, editors. Elephant seals: population ecology, behaviour, and physiology. Berkeley: University of California Press; 1994. p. 228–37.
Lisovski S, Hewson CM, Klaassen RHG, Korner-Nievergelt F, Kristensen MW, Hahn S. Geolocation by light: accuracy and precision affected by environmental factors. Meth Ecol Evol. 2012;3:603–12.
Rakhimberdiev E, Senner NR, Verhoeven MA, Winkler DW, Bouten W, Piersma T. Comparing inferences of solar geolocation data against high-precision GPS data: annual movements of a double-tagged black-tailed godwit. J Avian Biol. 2016;47:589–96.
Sumner MD, Wotherspoon SJ, Hindell MA. Bayesian estimation of animal movement from archival and satellite tags. PLoSONE. 2009;4:e7324.
Lisovski S, Gosbell K, Christie M, Hoye BJ, Klaassen M, Steward ID, Taysom AJ, Minton C. Movement patterns of Sanderling (Calidris alba) along the East Asian Australasian flyway and a comparison of methods to identify crucial areas for conservation. Emu. 2016;116:168–77.
R Core Team. R: a language and environment for statistical computing. 2015.
Wotherspoon SJ, Sumner MD, Lisovski S. R package SGAT: Solar/Satellite geolocation for animal tracking. GitHub Repository - http://github.com/swotherspoon/sgat 2013. Accessed 29 May 2015.
Wotherspoon SJ, Sumner M, Lisovski S. BAStag: basic data processing for light based geolocation archival tags, R package version 01–3. 2013.
Yamaura Y, Schmaljohann H, Lisovski S, Senzaki M, Kawamura K, Fujimaki Y, Nakamura F. Tracking the Stejneger’s stonechat Saxicola stejnegeri along the East Asian–Australian flyway from Japan via China to Southeast Asia. J Avian Biol. 2016;48:197–202.
Salewski V, Flade M, Poluda A, Kiljan G, Liechti F, Lisovski S, Hahn S. An unknown migration route of the ‘globally threatened’ aquatic warbler revealed by geolocators. J Ornithol. 2013;154:549–52.
Cramp S. Handbook of the birds of Europe, the Middle East and North Africa. The birds of the Western Palearctic, vol. 5. Oxford: Oxford University Press; 1988.
Green M, Alerstam T, Clausen P, Drent R, Ebbinge BS. Dark-bellied Brent Geese Branta bernicla bernicla, as recorded by satellite telemetry, do not minimize flight distance during spring migration. Ibis. 2002;144:106–21.
Bruderer B, Boldt A. Flight characteristics of birds: I. radar measurements of speeds. Ibis. 2001;143:178–204.
Lisovski S, Hahn S. GeoLight - processing and analysing light-based geolocator data in R. Meth Ecol Evol. 2012;3:1055–9.
Kanamitsu M, Ebisuzaki W, Woollen J, Yang S-K, Hnilo JJ, Fiorino M, Potter GL. NCEP-DOE AMIP-II reanalysis (R-2). Bull Am Meteorol Soc. 2002;77:437–70.
Kemp MU, van Loon E, Shamoun-Baranes J, Bouten W. RNCEP: global weather and climate data at your fingertips. Meth Ecol Evol. 2012;3:65–70.
Schmaljohann H, Becker PJJ, Karaardic H, Liechti F, Naef-Daenzer B, Grande C. Nocturnal exploratory flights, departure time, and direction in a migratory songbird. J Ornithol. 2011;152:439–52.
Schmaljohann H, Liechti F, Bruderer B. Trans-Sahara migrants select flight altitudes to minimize energy costs rather than water loss. Behav Ecol Sociobiol. 2009;63:1609–19.
Chelton DB, GFreeilich MH. Scatterometer-based assessment of 10-m wind analyses from the operational ECMWF and NCEP numerical weather prediction models. Mon Weather Rev. 2005;133:409–29.
Bromwich DH, Wang S-H. Evaluation of the NCEP–NCAR and ECMWF 15- and 40-Yr reanalyses using rawinsonde data from two independent Arctic field experiments. Mon Weather Rev. 2005;133:3562–78.
Schmaljohann H, Bruderer B, Liechti F. Sustained bird flights occur at temperatures beyond expected limits of water loss rates. Anim Behav. 2008;76:1133–8.
Mateos-Rodríguez M, Liechti F. How do diurnal long-distance migrants select flight altitude in relation to wind? Behav Ecol. 2012;23:403–9.
Kemp MU, Shamoun-Baranes J, van Loon EE, McLaren JD, Dokter AM, Bouten W. Quantifying flow-assistance and implications for movement research. J Theor Biol. 2012;308:56–67.
Babak N. usdm: Uncertainty analysis for species distribution models. R package version 11–12. 2013.http://CRAN.R-project.org/package=usdm. Accessed 29 May 2015.
Zuur AE, Irwin DE, Elphick CS. A protocol for data exploration to avoid common statistical probelms. Meth Ecol Evol. 2010;1:3–14.
Bates D, Mächler M, Bolker B, Walker S. lme4: Linear mixed-effects models using Eigen and S4. R package version 11–7. 2014. https://cran.r-project.org/web/packages/lme4/index.html. Accessed 29 May 2015.
Korner-Nievergelt F, Roth T, von Felten S, Guélat J, Almasi B, Korner-Nievergelt P. Bayesian data analysis in ecology using linear models with R, BUGS, and Stan. London: Elsevier; 2015.
Forstmeier W, Wagenmakers E-J, Parker TH. Detecting and avoiding likely false-positive findings – a practical guide. Biol Rev. 2016:10.1111/brv.12315.
Gelman A, Hill J. Data analysis using regression and multilevel/hierarchical models. Cambridge: Cambridge University Press; 2007.
Fox J. Effect displays in R for generalised linear models. J Stat Softw. 2003;8:1–27.
Fudickar AM, Wikelski M, Partecke J. Tracking migratory songbirds: accuracy of light-level loggers (geolocators) in forest habitats. Meth Ecol Evol. 2012;3:47–52.
Åkesson S, Walinder G, Karlsson L, Ehnbom S. Nocturnal migratory flight initiation in reed warblers Acrocephalus scirpaceus: effect of wind on orientation and timing of migration. J Avian Biol. 2002;33:349–57.
Dierschke V, Delingat J. Stopover behaviour and departure decision of northern wheatears, Oenanthe oenanthe, facing different onward non-stop flight distances. Behav Ecol Sociobiol. 2001;50:535–45.
Bairlein F, Eikenaar C, Schmaljohann H. Routes to genes: unravelling the control of avian migration—an integrated approach using Northern Wheatear Oenanthe oenanthe as model organism. J Ornithol. 2015;156:S3–S14.
Sapir N, Horvitz N, Wikelski M, Avissar R, Mahrer Y, Nathan R. Migration by soaring or flapping: numerical atmospheric simulations reveal that turbulence kinetic energy dictates bee-eater flight mode. Proc R Soc Lond B. 2011;278:3380–6.
Fransson T. Patterns of migratory fuelling in Whitethroats Sylvia communis in relation to departure. J Avian Biol. 1998;29:569–73.
Tsvey A, Bulyuk VN, Kosarev V. Influence of body condition and weather on departures of first-year European robins, Erithacus rubecula, from an autumn migratory stopover site. Behav Ecol Sociobiol. 2007;61:1665–74.
Richardson WJ. Timing and amount of bird migration in relation to weather: a review. Oikos. 1978;30:224–72.
Cohen EB, Pearson SM, Moore FR. Effects of landscape composition and configuration on migrating songbirds: inference from an individual-based model. Ecol Appl. 2014;24:169–80.
Cohen EB, Moore FR, Fischer RA. Experimental evidence for the interplay of exogenous and endogenous factors on the movement ecolgy of a migrating songbird. PLoSONE. 2012;7:e41818.
Panov IN, Chernetsov N. Migratory strategy of bluethroats, Luscinia svecica, in eastern fennoscandia. Part 2: response to acoustic markers and habitat selection at stopover. Proc Zool Inst RAS. 2010;314:173–83.
Bibby CJ, Green RE. Autumn migration strategies of reed and sedge warblers. Ornis Scand. 1981;12:1–12.
DeLuca WV, Woodworth BK, Rimmer CC, Marra PP, Taylor PD, McFarland KP, Mackenzie SA, Norris DR. Transoceanic migration by a 12 g songbird. Biol Lett. 2015;11:20141045.
Tøttrup AP, Klaassen RHG, Kristensen MW, Strandberg R, Vardanis Y, Lindström Å, Rahbek C, Alerstam T, Thorup K. Drought in Africa caused delayed arrival of European songbirds. Science. 2012;338:1307.
Bairlein F. How to get fat: nutritional mechanisms of seasonal fat accumulation in migratory songbirds. Naturwissenschaften. 2002;89:1–10.
Bairlein F, Gwinner E. Nutritional mechanisms and temporal control of migratory energy accumulation in birds. Ann Rev Nutr. 1994;14:187–215.
Maggini I, Bairlein F. Endogenous rhythms of seasonal migratory body mass changes and nocturnal restlessness in different populations of northern wheatear Oenanthe oenanthe. J Biol Rhyth. 2010;25:268–76.
Eikenaar C, Tsvey A, Schmaljohann H. Faster spring migration in northern wheatears is not explained by an endogenous seasonal difference in refueling rates. J Avian Biol. 2015;46:616–21.
Charmantier A, Gienapp P. Climate change and timing of avian breeding and migration: evolutionary versus plastic changes. Evol Appl. 2014;7:15–28.
Charmantier A, McCleery RH, Cole LR, Perrins CM, Kruuk LEB, Sheldon BC. Adaptive phenotypic plasticity in response to climate change in a wild bird population. Science. 2008;320:800–3.
Newton I. The migration ecology of birds. London: Academic; 2008.
Easterling DR, Meehl GA, Parmesan C, Changnon SA, Karl TR, Mearns LO. Climate extremes: observations, modeling, and impacts. Science. 2000;289:2068–74.
Gienapp P, Teplitsky C, Alho JS, Mills JA, Merilä J. Climate change and evolution: disentangling environmental and genetic responses. Mol Ecol. 2008;17:167–78.
Newson SE, Moran NJ, Musgrove AJ, Pearce-Higgins JW, Gillings S, Atkinson PW, Miller R, Grantham MJ, Baillie SR. Long-term changes in the migration phenology of UK breeding birds detected by large-scale citizen science recording schemes. Ibis. 2016;158:481–95.
NCEP Reanalysis data were provided by the NOAA/OAR/ESRL PSD, Boulder, CO, USA, from the website at http://www.esrl.noaa.gov/psd/. We are deeply indebted to Beth Pattinson, USFWS, and Robert Gill, USGS, for advice with permits and Susan Sharbaugh, Alaska Bird Observatory, and Kevin Winker, University of Alaska Museum, for logistic support. M. Bulte, T. Klinner, F. Müller, R. Nagel, T. Obracay and J. L. Schläfke provided excellent support in the field. We thank Axios Review for feedback and comments of four anonymous reviewers and Alice Risely helped to improve the manuscript.
The Deutsche Forschungsgemeinschaft (SCHM 2647, BA 816) provided financial support for this work, but it had no influence on the design of the study, analyses and interpretation of data and in writing the manuscript.
Availability of data and materials
The datasets generated and/or analysed during the current study are available at the GitHub page of Simeon Lisovski [https://github.com/slisovski]. All raw light-level geolocation data are available upon request from “Movebank” (www.movebank.org, project: Oenanthe oenanthe Alaska).
HS conceived the study. HS and FB collected the data. HS and SL analysed the light-level geolocation data. HS analysed output data. HS wrote the first draft of the manuscript, and all authors contributed substantially to revisions. All authors read and approved the final manuscript.
The authors declare that they have no competing interests.
Consent for publication
Ethics approval and consent to participate
Northern Wheatears were caught, ringed, and tagged with light-level geolocators under a license of the U.S. Fish and Wildlife Service (Federal Fish and Wildlife Permits: MB207892-0, MB97904A-0) and the State of Alaska Department of Fish and Game (Permits: 13–103, 14–009).
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Individual data of birds, table. (DOCX 39 kb)
Recapture probability of tagged Northern Wheatears, documentation. (DOCX 159 kb)
R code for analyzing light-level geolocation data, documentation. (DOCX 186 kb)
R code for analyzing phenotypic response to environmental cues, documentation. (DOCX 147 kb)
Migration routes of the Northern Wheatears, figures. (DOCX 239 kb)
Environmental conditions during autumn and spring migration, figures. (DOCX 108 kb)
Air temperature over remaining migration distance in autumn, figure. (DOCX 67 kb)
Precipitation and flow assistance before performing a migratory flight in spring, figures. (DOCX 63 kb)
Flow assistance during autumn migration, figure. (DOCX 47 kb)
About this article
Cite this article
Schmaljohann, H., Lisovski, S. & Bairlein, F. Flexible reaction norms to environmental variables along the migration route and the significance of stopover duration for total speed of migration in a songbird migrant. Front Zool 14, 17 (2017). https://doi.org/10.1186/s12983-017-0203-3