Long-Term Phenological Shifts in Raptor Migration and Climate

Abstract

Climate change is having a discernible effect on many biological and ecological processes. Among observed changes, modifications in bird phenology have been widely documented. However, most studies have interpreted phenological shifts as gradual biological adjustments in response to the alteration of the thermal regime. Here we analysed a long-term dataset (1980-2010) of short-distance migratory raptors in five European regions. We revealed that the responses of these birds to climate-induced changes in autumn temperatures are abrupt and synchronous at a continental scale. We found that when the temperatures increased, birds delayed their mean passage date of autumn migration. Such delay, in addition to an earlier spring migration, suggests that a significant warming may induce an extension of the breeding-area residence time of migratory raptors, which may eventually lead to residency.

About the Authors:

Mikaël Jaffré

* E-mail: [email protected]

Affiliation: Centre National de la Recherche Scientifique, Laboratoire d’Océanologie et de Géosciences UMR LOG CNRS 8187, Université Lille 1, Station Marine de Wimereux, Wimereux, France

Grégory Beaugrand

Affiliation: Centre National de la Recherche Scientifique, Laboratoire d’Océanologie et de Géosciences UMR LOG CNRS 8187, Université Lille 1, Station Marine de Wimereux, Wimereux, France

Éric Goberville

Affiliation: Centre National de la Recherche Scientifique, Laboratoire d’Océanologie et de Géosciences UMR LOG CNRS 8187, Université Lille 1, Station Marine de Wimereux, Wimereux, France

Frédéric Jiguet

Affiliation: Muséum National d’Histoire Naturelle, Centre de Recherches sur la Biologie des Populations d’Oiseaux, UMR MNHN-CNRS-UPMC 7204, Paris, France

Nils Kjellén

Affiliation: Department of Biology, University of Lund, Lund, Sweden

Gerard Troost

Affiliation: Trektellen.org, SOVON Dutch Centre for Field Ornithology, Nijmegen, The Netherlands

Philippe J. Dubois

Affiliation: Ligue pour la Protection des Oiseaux, Rochefort, France

Alain Leprêtre

Affiliation: Université Lille 1, Écologie Numérique et Écotoxicologie, UPRES EA 4515, Villeneuve d’Ascq, France

Christophe Luczak

Affiliations Centre National de la Recherche Scientifique, Laboratoire d’Océanologie et de Géosciences UMR LOG CNRS 8187, Université Lille 1, Station Marine de Wimereux, Wimereux, France, Université d’Artois, IUFM, Gravelines, France

Introduction

Current global climate change has already caused consistent patterns of change in the phenology and biogeography of numerous species, ranging from plants to vertebrates [1,2]. Such changes include pronounced shifts in the timing of the annual cycle events of birds, in particular breeding and migration [3-6]. Most studies have focused on spring migration, showing that migration occurs earlier [4] when temperatures increases [7]. The autumn migration has been less investigated and the responses in the timing of birds seem to be more complex, depending on the ecology and life history traits of the species [8,9]. The fluctuations in the timing of bird migration in response to climate change differs however between long-distance and short-distance migrants [3]. The later birds, spending the winter close to their breeding area, can alter the timing of their migration in response to climate change more rapidly than long-distance migrants [3], which overwinter in the tropics and whose migration is under endogenous control [10]. Several hypotheses have been proposed to explain the variability in post-breeding migration. A first hypothesis assumes an unchanged Breeding-Area Residence Time (so-called BART, Figure 1A) [11] and stipulates that the timing of post-breeding migration takes place earlier [9,12] as a consequence of a premature spring arrival (Figure 1B). A second hypothesis propounds that a more climatically suitable autumn may lead to a later departure because adverse environmental conditions are delayed (Figure 1C) [13,14].

[Figure omitted. See PDF.]

Figure 1. Schematic diagram of the length of the Breeding Area Residence Time (BART) depending on the initial condition and both hypotheses for changes in migration dates.

A, Initial state: migration dates are unchanged; the length of the BART is stable. B, Hypothesis 1: both spring and autumn migration dates occur earlier; the length of the BART is stable but lagged. C, Hypothesis 2: spring migration date is earlier, autumn migration date is later; the length of the BART is extended.

https://doi.org/10.1371/journal.pone.0079112.g001

Materials and Methods

Biological data

To investigate long-term (1980-2010) changes in raptor migrations in Western Europe, we used data collected on a daily basis at thirteen watchsites located along the most important western european migration flyways. To reduce the number of missing data and optimise seasonal coverage [15], the thirteen watchsites were rearranged into five sub-regions (Figure 2). The bird abundances for each monitoring day and for each sub-region were calculated using the mean of the abundances of the corresponding day in the corresponding sites (Table 1).

[Figure omitted. See PDF.]

Figure 2. Location of the five sub-regions used in the study:

A, Sweden. B, Netherlands. C, French Alps. D, Massif Central. E, Pyrenees.

Illustrations of the seven species of raptor from left to right: red kite, marsh Harrier, hen harrier, merlin, sparrowhawk, common buzzard, and kestrel.

https://doi.org/10.1371/journal.pone.0079112.g002

Sub-region Count season Mean days season^-1 (SD) Number of years Names of watch-sites Latitude Longitude

Sweden 01Aug.-10Nov. 110 (4) 29 Falsterbo 55.38°N 12.82°E

Netherlands 01Jul.-31Dec. 137 (11) 31 Camperduin 52.73°N 4.64°E

Kennemerduinen 52.42°N 4.55°E

Scheveningen 52.12°N 4.28°E

French Alps 15Jul.-30Nov. 96 (36) 26 Les Conches 46.18°N 5.32°E

Défilé de l’Écluse 46.12°N 5.90°E

Massif Central 01Aug.-31Oct. 92 (26) 24 Saint-Gervais-d’Auvergne 46.03°N 2.82°E

Montagnes de la Serre 45.66°N 3.09°E

Col de Baracuchet 45.58°N 3.92°E

Creste 45.55°N 3.04°E

Pyrenees 15Jul.-15Nov. 113 (8) 30 Col de Lizarrieta 43.26°N 1.62°W

Redoute de Lindux 43.03°N 1.36°W

Col d’Organbidexka 43.03°N 1.02°W

Table 1. Main characteristics of the sub-regions and watch-sites used in this study.

CSV

Download CSV

We chose seven partial migrant raptor species [16] with a sufficient abundance on the monitored watchsites: eurasian sparrowhawk Accipiter nisus, common buzzard Buteo buteo (subspecies buteo), marsh harrier Circus aeruginosus (subspecies aeruginosus), hen harrier Circus cyaneus, merlin Falco columbarius, common kestrel Falco tinnunculus and red kite Milvus milvus. Only years with at least 75% of days with counts during the migration window (2.5% to 97.5% of the passage) were retained.

Air temperatures data

Mean monthly air surface temperatures data (1980-2010) originated from the National Oceanic & Atmospheric Administration (NOAA, http://www.esrl.noaa.gov/psd/). These data were retrieved on a grid 2.5°longitude x 2.5°latitude for an area corresponding to Western Europe (12°W - 30°E and 40°N - 70°N).

Estimation of long-term phenological changes (1980-2010)

The study of long-term phenological changes in raptor post-breeding migration was undertaken using the Mean Passage Date (MPD) [17]. This technique provides accurate estimations of phenological shifts, being robust to variation in sample sizes and imperfect detectability [17].

Only a few statistical techniques exist to analyse such complex tables, e.g. co-inertia analysis [18] and three-mode principal component analysis (PCA; [19]). The later method was applied here. The three-mode PCA first calculates 3 classical PCAs on 2-dimensional tables after having transformed one table to ensure that the total inertia is identical in each mode [20]. The analysis then relates the different modes by assessing a core matrix calculated from the eigenvectors of each mode (see Beaugrand et al. 2000 [20] for a detailed description of the algorithm). In the present study, the application of this statistical technique allowed in a single analysis (1) the characterisation of temporal changes in the timing of autumn migration by the examination of the first two principal components (PCs; Figure 3A and Figure S1A) and (2) both the identification of species that contribute to the changes and the recognition of sites mainly influenced by the temporal patterns (associated normalised eigenvectors – Figure 3B and Figure S1B). To assess long-term spatial and temporal changes in the MPD of the seven species in western Europe, we only calculated one standardised PCA (i.e. subtracting the mean and dividing by the standard deviation (SD) to give all parameters the same variation) on the deployed 3-way matrix 31 years x (7 species x 5 sub-regions), which represents the first stage performed in a three-mode PCA [20-22].

[Figure omitted. See PDF.]

Figure 3. Long-term changes in the timing of autumn migration in relation to autumn temperatures in Europe from 1980 to 2010.

A, B, PCA of the year-to-year changes in the mean passage date of seven species of raptors from 1980 to 2010. (A) Long-term changes in the first principal component (PC1), which represents 20.24% of the total variability. Abrupt shifts, identified by the Rodionov’s test, are superimposed (grey line). (B) Bar diagram of the first eigenvector normalised between -1 and 1. Negative values are in blue and positive values are in red. C, D, PCA of the year-to-year changes in the August-October temperatures from 1980 to 2010. (C) Long-term changes in the first principal component (PC1), which represents 44.25% of the total variability. Abrupt shifts, identified by the Rodionov’s test, are superimposed (grey line). (D) Mapping of the first normalised eigenvector, which shows the correlation between autumn temperatures and the first component. Negative values are in blue and positive values are in red.

https://doi.org/10.1371/journal.pone.0079112.g003

Long-term changes in air temperatures

To evaluate the influence of climate on the phenology of post-breeding migration, we examined long-term changes in autumn and spring temperatures in Western Europe by standardized PCA (Joliffe 1986). The analyses were performed in the spatial domain ranging from 10°W to 30°E and from 40°N to 70°N and data were averaged on a grid of 2.5° longitude x 2.5° latitude for each year of the period 1980-2010. We selected (1) August to October temperatures (autumn) because this period corresponds to the time between the emancipation of young and the migration period and (2) April temperatures (spring) as the timing of breeding during this month could be linked with the timing of autumn migration [23]. We then performed two standardised PCAs with the double objective of identifying major long-term changes in both autumn and spring temperatures (examination of the first PCs, Figure 3C) and locating their geographical pattern (mapping of the first normalised eigenvectors, Figure 3D).

Identification and quantification of temporal discontinuities

We applied the Rodionov's sequential algorithm to detect the temporal discontinuities (1980-2010) occurring in the first PCs originating from the PCA performed on post-breeding migration (Figure 3A and Figure S1A), autumn temperatures (Figure 3C) and spring temperatures (Figure S1C). We used a first-order autoregressive model to consider temporal autocorrelation [24], a probability threshold of p=0.1 and cut-off length of l=10 to control the magnitude and the scale of the shift, respectively. To consider outliers above 1 standard deviation, the Huber's weight parameter was fixed to h=1. As guidance for choosing parameters p and l, setting parameters were adjusted following the recommendations of previous studies to detect shifts with meaningful environmental and biological implications [24-26]. The choice of a significance level p=0.1 is a good compromise between robustness of the detected shift and the flexibility of the detection [24,25]. A cut-off length of l=10 discounts regimes of shorter length to detect only statistically significant long regime shifts depending on the studied period (i.e. 30 years - [26]). A Huber parameter of h=1 is such that the number of detected shifts remains stable for higher values of h. More details about the methods and parameters it takes can be found at www.beringclimate.noaa.gov/regimes/. The magnitude of the shifts was subsequently quantified by selecting Sweden and Pyrenees (Table 1), sub-regions for which complete datasets exist for all considered raptors. For each regime, MPD anomalies were assessed for every year and species with respect to the period 1980-2010. A histogram was then calculated to reveal the statistical distribution of all anomalies for each regime (Figure 4). The average anomalies of autumn temperatures and post-breeding migration were also assessed for each regime (Figure 4A).

[Figure omitted. See PDF.]

Figure 4. Effect of the abrupt shifts in autumn temperatures on the timing of migration.

A, mean anomalies of autumn temperatures and Mean Passage Date (MPD; in day) for all periods revealed by the Rodionov’s test (see Figure 2). Both autumn temperatures and MPD are expressed in anomalies with respect to the period 1980-2010. b, c, d, Repartition of the anomalies of MPD during (B) the first period 1980-1994 (T1=12.29°C), (C) the second period 1995-2006 (T2=12.91°C) and (D) the third period 2007-2010 (T3=12.66°C).

https://doi.org/10.1371/journal.pone.0079112.g004

Relationships between long-term changes in post-breeding migration and temperatures

Linear correlations were calculated to assess the relationships between long-term changes in the first two PCs of post-breeding migration and changes in both April and autumn temperatures. The probabilities of the correlations were corrected to account for temporal autocorrelation [21].

Results and Discussion

Year-to-year changes in the first Principal Component performed on MPD (PC1; 20.24% of the total variability) revealed a period of relative early passage until 1994, followed by a period of later bird migration until 2006 and finally a return to an earlier MPD (Figure 3A). The examination of the first normalised eigenvector indicates that the long-term phenological shifts exhibited by PC1 concerned virtually all species and sub-regions (Figure 3B). The Rodionov's test applied on PC1 shows that short-distance migrant raptors have abruptly delayed their MPD in Western Europe in 1995 and returned to an earlier autumn migration in 2006. The second Principal Component (PC2; 13.03% of the total variability) revealed a tendency for these birds (except red kite Milvus milvus) to an earlier MPD interrupted in 2006 (Figure S1A). The examination of the second normalised eigenvector indicates an opposition between the long-term phenological shifts in northern sites (the Netherlands, the French Alps and to a lesser extent in Sweden) on one side, and in southern sites (the Massif Central and Pyrenees) on the other side (Figure S1B).

Mapping of the first eigenvectors of the PCA performed on temperatures in Western Europe indicated that the trends detected by the PCAs reflect well changes in autumn temperatures in the north-western part of Europe (Figure 3D) and mostly changes in spring temperatures in the northern part of Europe (Figure S1D).

We subsequently applied the Rodionov's test on the PC1 of autumn temperatures (44.25% of the total variability) and on the PC1 of spring temperatures (34.32% of the total variability) to detect temporal discontinuities. For autumn temperatures, the analysis detected the same temporal discontinuities than those identified previously by the first PC of the MPD (see Figure 3, A and C): a relative cold regime until 1994, followed by a warmer regime until 2006 and a return to cooler conditions from 2007 to 2010 (Figure 3C). For spring temperatures, the analysis only found one single temporal discontinuity: in 1999. A relative cold regime occurred until 1998 and was subsequently followed by a warmer regime (Figure S1C). No synchronous discontinuities were detected between the second PC of the MPD and the first PC of spring and autumn temperatures. We also tested for linear relationships between both spring and autumn temperatures and the timing of post-breeding migration. For PC1 of the MPD, no linear correlation was detected between either spring (r=-0.23, PACF=0.28) or autumn (r=0.31, PACF=0.18) temperatures. For PC2 of the MPD, while no linear correlation was detected with autumn temperatures (r=0.44, PACF=0.12), a significant relationship was found with spring temperatures (r=-0.55, PACF=0.03). These results provide evidence that, although part of phenological changes were correlated with the spring temperatures (i.e. linear relationship with the CP2 of the MPD), most short-distance migrants responded abruptly to pronounced shifts in autumn temperatures during the period 1980-2010 (i.e. synchronous discontinuities with the CP1 of the MPD; Figure 3A and Figure 3C), suggesting that they adjusted their postnuptial migration timing to the thermal regime.

Three periods were detected: 1980-1994, 1995-2006 and 2007-2010 (Figure 4A). We quantified the magnitude of the shifts (1994-1995 and 2006-2007) by calculating the average of MPD and autumn air temperature for each period (Methods). While between 1980-1994 and 1995-2006 the autumn temperatures increased by 0.61°C, post-breeding migration expressed as MPD was delayed for all species but one (marsh harrier) from 0.2 to 8.9 days (mean = 3.23 days, or 5.28 days/°C; Figure 4). Marsh harrier exhibited a negative response probably because northern populations of this species are long-distance migrants [27]. This result is in agreement with a study, which showed that post-breeding migration became earlier during the last three decades for this species [12]. Between 1995-2006 and 2007-2010, the autumn temperatures diminished by 0.25°C and paralleled an earlier MPD ranging from 1.6 to 8.9 days (mean = -4.52 days, or -18.28 days/°C; Figure 4). Sparks and Menzel reported that a 1°C increase in temperature could advance arrival time by days [28]. Therefore, our observed changes of 5.28 days/°C (first shift) and -18.28 days/°C (second shift) are less conservative and suggest a high sensitivity of post-breeding migration of short-distance migrant raptors to climate-induced changes in autumn temperatures over their breeding area. It has been shown that genetic factors play an important role in the control of migratory behaviour of birds [29] and that this behaviour is partly hereditary [30]. Such a result implies that rapid changes in migratory behaviour can occur, especially for species with short generation time. However, the birds studied here are long-lived species and have a long generation time (e.g. five years for common buzzard Buteo buteo, i.e. up to six generations during the three decades [23]), which is inconsistent with the hypothesis of a positive selection of late migrants. Hence, abrupt shifts observed in this study are probably not a micro-evolutionary response of species. This statement is in accordance with the weak evolutionary potential found in raptors [31]. Instead, the rapid phenological response of short-distance migrant raptors observed here in autumn is probably a phenotypic response of these birds to improved foraging and environmental conditions at middle-high latitudes.

Other studies found a relationship between earlier autumn migration and spring temperatures [8,9,14]. Here we also found a significant correlation between the second PC, which exhibit an earlier MPD, and the increase of spring temperatures. This earlier passage of raptors is mostly prominent over northern sub-regions, where birds are known to be more strictly migratory. This means that the first hypothesis (unchanged BART) is also valid with our data, but is of lesser importance compared to the second hypothesis of increased BART (13.03 < 20.24% of the total variability). Most studies, including ours, failed to find a link between the autumn MPD of birds and the autumn temperatures on the breeding grounds [14,32]. However, our analysis on temporal discontinuity revealed an abrupt response of birds to changes in autumn temperatures. Our results suggest that not only BART may increase as a result of earlier spring arrival due to warmer conditions in spring [23,33] but also because European short-distance migratory raptors tend to delay their departure as a result of warmer conditions in autumn. In this context, it is likely that global warming will enable several species to reside on their breeding grounds some favourable years. According to the Intergovernmental Panel on Climate Change, a moderate scenario (A1B scenario) of global warming will increase the temperatures between 2.2°C and 5.3°C in Europe by 2100 [34]. Assuming a linear relationship between BART and autumn temperatures, it is plausible that such a rise in temperatures leads to a mean extension ranging from 11.6 (+2.2°C) to 28.0 (+5.3°C) days of BART in autumn for the seven studied species and by [32.1 - 77.4] days for the extreme case of hen harrier. Improved foraging in winter could enable birds overwintering closer to their breeding range which therefore reduce their migration distance [35,36]. The present manuscript adds a phenological aspect to the distance reduction, and both phenological changes and reduction of the migration distance suggest that climate change, associated with a decrease in seasonality, could have potential implications for some birds to forgo migration altogether, leading to residency [37]. This hypothesis is supported by the ‘Shifting Home’ model [36], a model of the origin of the migration system which suggests that the periods of global cooling during the Cenozoic drove shifts of the wintering area of some Palearctic birds and led to the emergence of migration [36]. In an opposite climatic context (a global warming instead of a global cooling), we observe an opposite feature. While the ‘Shifting Home’ model proposes that decreasing winter temperatures led birds to migrate towards more favourable areas, the current increase in temperature could result in a reduction of the migratory behaviour, resulting in the residency of species.

Supporting Information

[Figure omitted. See PDF.]

Figures S1.

Long-term changes in the timing of Autumn migration in relation to Spring temperatures in Europe from 1980 to 2010. A, B, PCA of the year-to-year changes in the mean passage date of seven species of raptors from 1980 to 2010. (A) Long-term changes in the second principal component (PC2), which represents 13.03% of the total variability. Abrupt shifts, identified by the Rodionov’s test, are superimposed (grey line). (B) Bar diagram of the second eigenvector normalised between -1 and 1. Negative values are in blue and positive values are in red. C, D, PCA of the year-to-year changes in the April temperatures from 1980 to 2010. (C) Long-term changes in the first principal component (PC1), which represents 34.32% of the total variability. Abrupt shifts, identified by the Rodionov’s test, are superimposed (grey line). (D) Mapping of the first normalised eigenvector, which shows the correlation between spring temperatures and the first component. Negative values are in blue and positive values are in red.

https://doi.org/10.1371/journal.pone.0079112.s001

(TIF)

Acknowledgments

We would like to thank the many birdwatchers and volunteers counting birds at the European watchsites for many years, and among them David Erterius (Falsterbo), Nick van der Ham (Camperduin ), Hans Groot, Eef Kieft and Pim de Nobel (Kennemerduinen), Michel de Lange, Rinse van der Vliet, Rob Westerduijn (Scheveningen), Pierre Crouzier and Marie-Paule de Thiersant (Les Conches), Xavier Birot-Colomb and Dr Pierre Charvoz (Défilé de l’Écluse), Madeleine Chaffraix (Saint-Gervais-d’Auvergne), Christophe Eymard (Montagnes de la Serre), Nicolas Lorenzini (Col de Baracuchet), Thomas Bernard (Creste), Jean-Paul Urcun, Laurent Couzi and the Observatoire Régional de la Migration des Oiseaux en Aquitaine (ORMO) (3 Pyrenean sites). The authors also thank the photographers who provided the photographs of Figure 1: Stephen Peten, Thierry Tancrez, Guy Flohart and Christophe Luczak. We are grateful to SOVON, Falsterbo Birdwatch, and Ligue pour la Protection des Oiseaux (LPO – www.migraction.net) for providing the biological data.

Author Contributions

Conceived and designed the experiments: MJ CL GB. Performed the experiments: MJ CL GB. Analyzed the data: MJ CL GB EG. Contributed reagents/materials/analysis tools: MJ GB EG NK GT PJD CL. Wrote the manuscript: MJ GB EG FJ NK GT PJD AL CL.

Citation: Jaffré M, Beaugrand G, Goberville É, Jiguet F, Kjellén N, Troost G, et al. (2013) Long-Term Phenological Shifts in Raptor Migration and Climate. PLoS ONE 8(11): e79112. https://doi.org/10.1371/journal.pone.0079112

References

1. Devictor V, van Swaay C, Brereton T, Chamberlain D, Heliölä J et al. (2012) Differences in the climatic debts of birds and butterflies at a continental scale. Nature. Clim Change 2: 121-124. doi:10.1038/nclimate1347.

2. Parmesan C (2006) Ecological and evolutionary responses to recent climate change. Annu Rev Ecol Syst 37: 637-669. doi:10.1146/annurev.ecolsys.37.091305.110100.

3. Lehikoinen E, Sparks TIM, Zalakevicius M (2004) Arrival and departure dates. Adv Ecol Res 35: 1-31. doi:10.1016/S0065-2504(04)35001-4.

4. Jonzén N, Lindén A, Ergon T, Knudsen E, Vik JO et al. (2006) Rapid advance of spring arrival dates in long-distance migratory birds. Science 312: 1959-1961. doi:10.1126/science.1126119. PubMed: 16809542.

5. Both C, Visser ME (2001) Adjustment to climate change is constrained by arrival date in a long-distance migrant bird. Nature 411: 296-298. doi:10.1038/35077063. PubMed: 11357129.

6. Coppack T, Both C (2002) Predicting life-cycle adaptation of migratory birds to global climate change. Ardea 90: 369-378.

7. Hüppop O, Hüppop K (2003) North Atlantic Oscillation and timing of spring migration in birds. Proc R Soc Lond B 270: 233-240. doi:10.1098/rspb.2002.2236. PubMed: 12614571.

8. Cotton PA (2003) Avian migration phenology and global climate change. Proc Natl Acad Sci U_S_A 100: 12219-12222. doi:10.1073/pnas.1930548100. PubMed: 14519854.

9. Jenni L, Kéry M (2003) Timing of autumn bird migration under climate change: advances in long–distance migrants, delays in short–distance migrants. Proc R Soc Lond B 270: 1467-1471. doi:10.1098/rspb.2003.2394.

10. Gwinner E (1996) Circannual clocks in avian reproduction and migration. Ibis 138: 47-63.

11. Thorup K, Tøttrup AP, Rahbek C (2007) Patterns of phenological changes in migratory birds. Oecologia 151: 697-703. doi:10.1007/s00442-006-0608-8. PubMed: 17160398.

12. Filippi-Codaccioni O, Moussus J-P, Urcun J-P, Jiguet F (2010) Advanced departure dates in long-distance migratory raptors. J Ornithol 151: 687-694. doi:10.1007/s10336-010-0500-5.

13. Mills AM (2005) Changes in the timing of spring and autumn migration in North American migrant passerines during a period of global warming. Ibis 147: 259-269. doi:10.1111/j.1474-919X.2005.00380.x.

14. Sokolov L (2006) Effect of global warming on the timing of migration and breeding of passerine birds in the 20th century. Entomol Rev 86: S59-S81. doi:10.1134/S0013873806100058.

15. Hussell DJT, Ralph СJ (2005) Recommended methods for monitoring change in landbird populations by counting and capturing migrants. North Birds Bander 30: 6-20.

16. Porter RF, Beaman MAS (1985) A résumé of raptor migration in Europe and the Middle East. In: I. NewtonRD Chancellor. Conservation Studies on Raptors: ICBP Tech. Pub. pp. 237-242.

17. Moussus JP, Julliard R, Jiguet F (2010) Featuring 10 phenological estimators using simulated data. Methods Ecol Evol 1: 140-150. doi:10.1111/j.2041-210X.2010.00020.x.

18. Dolédec S, Chessel D (1994) Co-inertia analysis: an alternative method for studying species–environment relationships. Freshw Biol 31: 277-294. doi:10.1111/j.1365-2427.1994.tb01741.x.

19. Hohn ME (1993) Principal component analysis of three-way data. In: JC DavisUC Herzfeld. Computers in Geology - 25 years of progress. New York: Oxford University Press. pp. 181-194.

20. Beaugrand G, Ibañez F, Reid PC (2000) Spatial, seasonal and long-term fluctuations of plankton in relation to hydroclimatic features in the English Channel, Celtic Sea and Bay of Biscay. Mar Ecol Prog S 200: 93-102. doi:10.3354/meps200093.

21. Beaugrand G, Brander KM, Lindley JA, Souissi S, Reid PC (2003) Plankton effect on cod recruitment in the North Sea. Nature 426: 661-664. doi:10.1038/nature02164. PubMed: 14668864.

22. Goberville E, Beaugrand G, Sautour B, Treguer P, Team S (2010) Climate-driven changes in coastal marine systems of western Europe. Marine Ecology Progress Series 408: 129-148.

23. Lehikoinen A, Byholm P, Ranta E, Saurola P, Valkama J et al. (2009) Reproduction of the common buzzard at its northern range margin under climatic change. Oikos 118: 829-836. doi:10.1111/j.1600-0706.2008.17440.x.

24. Rodionov SN (2006) Use of prewhitening in climate regime shift detection. Geophys Res Lett 33: L12707. doi:10.1029/2006GL025904.

25. Rodionov SN (2004) A sequential algorithm for testing climate regime shifts. Geophys Res Lett 31: L09204. doi:10.1029/2004GL019448.

26. Luczak C, Beaugrand G, Jaffré M, Lenoir S (2011) Climate change impact on Balearic shearwater through a trophic cascade. Biol Lett 7: 702-705. doi:10.1098/rsbl.2011.0225. PubMed: 21508020.

27. Strandberg R, Klaassen RHG, Hake M, Olofsson P, Thorup K et al. (2008) Complex timing of Marsh Harrier Circus aeruginosus migration due to pre-and post-migratory movements. Ardea 96: 159-171. doi:10.5253/078.096.0202.

28. Sparks TH, Menzel A (2002) Observed changes in seasons: an overview. Int J Climatol 22: 1715-1725. doi:10.1002/joc.821.

29. Pulido F, Berthold P, Van Noordwijk AJ (1996) Frequency of migrants and migratory activity are genetically correlated in a bird population: evolutionary implications. Proc Natl Acad Sci U_S_A 93: 14642–14647. doi:10.1073/pnas.93.25.14642. PubMed: 8962107.

30. Berthold P, Pulido F (1994) Heritability of migratory activity in a natural bird population. Proc R Soc Lond B 257: 311-315. doi:10.1098/rspb.1994.0131.

31. Mihoub JB, Mouawad NG, Pilard P, Jiguet F, Low M et al. (2012) Impact of temperature on the breeding performance and selection patterns in lesser kestrels Falco naumanni. J Avian Biol 43: 472-480. doi:10.1111/j.1600-048X.2012.05529.x.

32. Péron G, Henry PY, Provost P, Dehorter O, Julliard R (2007) Climate changes and post-nuptial migration strategy by two reedbed passerines. Clim Res 35: 147-157. doi:10.3354/cr00721.

33. Sergio F (2003) Relationship between laying dates of black kites Milvus migrans and spring temperatures in Italy: rapid response to climate change? J Avian Biol 34: 144-149. doi:10.1034/j.1600-048X.2003.03104.x.

34. Christensen JH, Hewitson B, Busuioc A, Chen A, Gao X et al. (2007) Regional Climate Projections. In: S. SolomonD. QinM. ManningZ. ChenM. Marquis. Climate; Change (2007) The Physical Science Basis. Cambridge, UK: Cambridge University Press. pp. 847-940.

35. La Sorte FA, Thompson FR (2007) Poleward shifts in winter ranges of North American birds. Ecology 88: 1803-1812. doi:10.1890/06-1072.1. PubMed: 17645026.

36. Louchart A (2008) Emergence of long distance bird migrations: a new model integrating global climate changes. Naturwissenschaften 95: 1109-1119. doi:10.1007/s00114-008-0435-3. PubMed: 18712337.

37. Pulido F, Berthold P (2010) Current selection for lower migratory activity will drive the evolution of residency in a migratory bird population. Proc Natl Acad Sci U_S_A 107: 7341-7346. doi:10.1073/pnas.0910361107. PubMed: 20368446.