- Research article
- Open Access
Seasonal rainfall at long-term migratory staging sites is associated with altered carry-over effects in a Palearctic-African migratory bird
BMC Ecology volume 16, Article number: 41 (2016)
An understanding of year-round habitat use is essential for determining how carry-over effects shape population dynamics in long-distance migratory songbirds. The recent discovery of long-term migratory staging sites in many species, prior to arrival at final wintering sites, adds complexity to efforts to decipher non-breeding habitat use and connections between sites. We investigated whether habitat conditions during migratory staging carry over to influence great reed warbler (Acrocephalus arundinaceus) body condition at final wintering sites in Zambia. We asked whether the presence/absence and strength of such carry-over effects were modified by contrasting rainfall conditions during 2 years.
First, we found that individuals staging in a dry year had higher corticosterone (CORTf) and stable nitrogen isotope values (suggesting higher aridity) than birds staging in a wet year, indicating that regional weather affected staging conditions. Second, we found that carry-over effects from staging habitat conditions (measured via carbon and nitrogen isotopes) to final winter site body condition (measured via scaled mass index and β-hydroxybutyrate) were only present in a dry year, suggesting that environmental factors have consequences for the strength of carry-over effects. Our results also suggest that wet conditions at final winter sites may buffer the effects of poor staging conditions, at least in the short term, since individuals that staged in a dry year had higher scaled mass indices in Zambia than individuals that staged in a wet year.
This study provides a first insight into the connections between long-term migratory staging sites and final wintering sites, and suggests that local environmental factors can modify the strength of carry-over effects for long-distance migratory birds.
The population dynamics of migratory animals may be affected by factors operating at multiple discrete sites throughout the annual cycle. Long-distance migrants must contend with extreme variation in weather conditions, habitat types and habitat quality, imposed by sites thousands of kilometers apart. This variation can be exacerbated by anthropogenic impacts that disrupt links between specific habitats and crucial periods of the annual cycle . Breeding conditions are well known to have important consequences for population dynamics, but evidence is growing that conditions experienced during non-breeding seasons may be at least as important , such that conditions during migratory staging may carry over to affect future breeding success . However, identifying such carry-over effects (defined as processes in one season influencing the success of an individual in the following season) between migratory staging sites has been challenging owing to difficulties in tracking small-bodied migrants across multiple sites throughout the annual cycle .
There is strong evidence that winter weather conditions, specifically rainfall effects on food availability, play an important role in phenology and reproductive success of migratory birds, such that low precipitation during the winter has long-term effects on population dynamics [4–7]. However, understanding individual carry-over effects from wintering ground conditions requires knowledge of movements within the non-breeding season and the specific habitats thus encountered. Recently, new light has been shed on the largely unknown winter migratory patterns of long-distance migratory songbirds, using miniature light-level geolocators . Surprisingly, geolocator data show that rather than single wintering sites or “itinerant” winter movements , the use of long-term staging sites en route to core wintering grounds may be widespread among migratory songbirds [10–14]. Long-term staging sites complicate our understanding of year-round habitat use for migratory birds because variation in habitat quality among staging sites may have important consequences for population dynamics via carry-over effects on individual survival and condition.
Studying how conditions at multiple sites cumulatively affect year-round individual success is particularly important for Palearctic-African species conservation, management, and effective reserve design, as sub-Saharan wintering migrants are in a state of severe population decline [15, 16]. To begin addressing these shortcomings we studied the great reed warbler (Acrocephalus arundinaceus), a 30 g long-distance Palearctic-African migratory songbird that winters widely across sub-Saharan Africa, in southern Zambia. During autumn migration, great reed warblers pause for 2–3 months to moult their flight feathers before continuing on to their final winter sites further south where moult is completed . We focused on rainfall as a potential driver of carry-over effects because its effect on food availability makes it a critical environmental factor for insectivorous birds . We studied great reed warblers during 2 years with contrasting rainfall conditions at staging areas (estimated from isotopic data; see “Results”). Carry-over effects are unlikely to be uniform across years and our understanding of the processes that are capable of modifying their strength remains limited. Natural variation between our study years allowed us to investigate whether the presence and strength of carry-over effects vary depending on environmental conditions.
We first established that differences in rainfall between years affected the conditions birds experienced while staging. To do so, we utilised the fact that great reed warblers replace their flight feathers during staging , such that individuals arriving in Zambia carry an isotopic and hormonal signature of the conditions they experienced previously during staging. Specifically, we analysed stable isotope ratios of nitrogen (δ 15N) and carbon (δ 13C), and corticosterone (CORTf) levels in flight feathers grown on staging grounds (Fig. 1). Both δ 15N and δ 13C isotope ratios have been used previously to assess habitat conditions experienced by migratory birds [20–22], trophic position in food webs, or marine versus freshwater ecosystems . δ 15N increases with greater aridity, especially in Africa [24–26], and δ 13C decreases as C3-dominated plant communities, which are associated with moister and cooler environments, become less common . Corticosterone is a hormone involved in energy management and high levels can indicate compromised individual condition, reflecting a response to stressors such as reduced food availability and poor habitat quality [28–30]. CORTf is an integrated measure of the hormone and reflects conditions experienced over the feather growth period ([31, 32]; for a review, see ). We expected that if birds indeed experienced drier conditions during staging in 2011, this would be recorded in their flight feathers as higher δ 15N values (reflecting higher aridity) and higher CORTf (reflecting a physiological response to poor conditions); however, we expected no difference in δ 13C values since carbon isotopes are relatively stable between years .
Second, we asked whether variation in environmental conditions experienced during staging affected the presence and strength of carry-over effects. We used several physiological metrics to estimate condition on wintering grounds. Plasma metabolites are indicators of energetic state: high circulating levels of triglycerides in plasma (hereafter ‘TRIG’) indicate fat deposition and foraging rates over short time-scales [35, 36], whereas high circulating levels of β-hydroxybutyrate (hereafter ‘BUTY’) in plasma indicate dietary fasting, as BUTY is synthesized from fatty acids and largely replaces glucose in fueling metabolism during fasting [37, 38]. We also used scaled mass index  to estimate the amount of stored fat. Finally, to gain insight into stress levels on the final wintering site, we measured CORTf levels in body feathers which, unlike flight feathers, are replaced in Zambia after staging (Fig. 1). Since poorer staging conditions should exacerbate differences in quality between individuals, we predicted that in 2011, a year with below-average rainfall, carry-over effects would be stronger than in 2012, a year with above-average rainfall (Fig. 2).
Estimated location of great reed warbler staging sites
We captured 26 great reed warblers in 2011 and 38 in 2012. The triple isotope cluster analysis assigned the majority of great reed warbler staging locations to cluster 2 (79 %), corresponding to Angola, Zambia, Zimbabwe, and southern portions of the Democratic Republic of the Congo. Some birds were assigned to cluster 1 (8 %) and cluster 3 (13 %), but none were assigned to cluster 4 which is associated with the Horn of Africa and the north-eastern Sahel region (Table 1; see  for a cluster map of Africa). A recent geolocator study of Swedish breeding great reed warblers staging in West Africa found no further migration to southern Africa after staging and moult , which concurs with ringing data  and suggests that birds wintering in southern Africa do not stage in West Africa prior to reaching their final wintering grounds, and originate from different European breeding grounds. Lemke et al.  also found that the mean distance moved after staging was 678 km. Taken together, these results suggest that the most common staging sites for our study population are within 678–2000 km north of the Zambian study site (extent of most of cluster 2 excluding West African component of the cluster), and do not incorporate cluster 4 (East Africa) or West Africa.
Staging site rainfall in 2011 and 2012
Rainfall within estimated staging sites during autumn migration differed between the 2 years of this study. Birds wintering in Zambia in 2011 experienced lower than average rainfall during autumn migration staging, whereas birds wintering in 2012 experienced higher than average rainfall (Fig. 2).
Does staging rainfall affect conditions experienced by great reed warblers?
In support of the NOAA rainfall data (reported in Fig. 2; ), isotopic results also suggest that birds indeed experienced drier conditions in 2011 than in 2012: δ 15N values during staging were higher in 2011, suggesting that more arid conditions were experienced by birds in 2011 than in 2012 (t test: t′54.4 = 2.9, p = 0.006; 2011: mean = 9.9 ‰, SE ± 0.26; 2012: mean = 8.9 ‰, SE ± 0.22). In addition, as predicted, staging CORTf (feather corticosterone) was higher for 2011 than for 2012 (t test on ranked data: t′52.9 = 6.9, p < 0.0001; 2011: mean = 83.1 pg/mm, SE ± 2.3; 2012: mean = 65.0 pg/mm, SE ± 3.3; Fig. 3) while CORTf from feathers grown on the final wintering grounds in Zambia was similar between years (t test on ranked data: t′48.5 = −0.9, p = 0.35; 2011: mean = 21.7 pg/mm, SE ± 1.86; 2012: mean = 21.7, SE ± 0.68; Fig. 3). As expected there was also no difference between years in δ 13C values from feathers grown during staging (δ 13C = t test: t′47.5 = −0.3, p = 0.75; 2011: mean = −15.7 ‰, SE ± 0.76; 2012: mean = −15.4 ‰, SE ± 0.53).
Does yearly variation in rainfall affect the strength of carry-over effects?
With scaled mass index (on the Zambian final winter site) as the response variable, interaction terms included in the top model set were δ 13C * year, and, δ 15N *year. Other models within 2 AICc scores included staging CORTf, with combinations of year or δ 15N (Additional file 1: Table S1). The interaction term between δ 13C and year was significant in the multivariate model selected using AICc, whereas the interaction between δ 15N and year was not (Table 2). To better understand the relationship between δ 13C and year, we regressed scaled mass index against δ 13C separately for 2011 and 2012. The relationship between scaled mass index and δ 13C was negative in 2011 but not in 2012 (2011: r2 = 0.17, p = 0.03; 2012: r2 = 0.03, p = 0.37; Fig. 4a), indicating that in 2011 birds staging in areas with fewer C3 plants (associated with moister, cooler habitats) had lower scaled mass indices in Zambia. δ 15N was not related to scaled mass index in either 2011 (r2 = 0.04, p = 0.32) or 2012 (r2 = 0.009, p = 0.60; Fig. 4b).
With BUTY as the response variable, every model in the top model set included the interaction term δ 15N * year (Additional file 1: Table S1). The interaction between δ 15N and year was significant in both of the selected multivariate models (Table 2). When years were analysed separately and scaled mass index and time at capture were included as covariates (see “Methods” section), the relationship between δ 15N and BUTY was weakly positive in 2011 but not in 2012 (2011: t = 1.83, p = 0.08; 2012: t = −1.34, p = 0.2; Fig. 4c). This may indicate that in 2011, individuals staging in habitats with higher δ 15N values (corresponding to higher aridity) had higher BUTY levels (corresponding to increased fat catabolism and fasting) in Zambia.
With TRIG as the response variable, the interaction between δ 13C * year was included in the top model set. Other models within 2 AICc scores included year, with combinations of δ 13C, δ 15N or staging CORTf (Additional file 1: Table S1). The interaction between δ 13C and year was not significant in the multivariate general linear model included in the top model set (Table 2). When years were analysed separately with scaled mass index and time at capture included as covariates (see “Methods” section), δ 13C was negatively related to TRIG in 2011 but not in 2012; however, these relationships were not significant (2011: t = −1.7, p = 0.1; 2012: t = 0.34, p = 0.74; Fig. 4d).
With winter CORTf as the response variable, no interaction terms were included in the top model set (Table 2). Models within 2 AICc scores included combinations of carbon, nitrogen, and staging CORTf (Additional file 1: Table S1).
The relationships between staging conditions and final winter site conditions in Zambia did not differ between males and females in either 2011 or 2012 (Additional file 1: Table S3).
Zambian final winter site rainfall and body condition
We tested whether body condition at the final wintering site was better predicted by environmental conditions experienced during staging, or experienced at the final wintering site. Local rainfall conditions at our study site in Zambia were higher in 2011 (December–February total = 736 mm) than in 2012 (595 mm), which was the reverse of rainfall conditions on the staging grounds (Fig. 2). We found that winter CORTf, BUTY, and TRIG did not differ between years (winter CORTf: t test on ranked data t′48.5 = −0.92, p = 0.36; BUTY: ANCOVA F1,51 = 0.08, p = 0.79; TRIG: ANCOVA F1,48 = 1.06, p = 0.31); however, the great reed warbler population in Zambia in 2012 had a lower scaled mass index (an estimate of stored fat) than the population in 2011 (t test t′56.3 = 2.4, p = 0.02). This may suggest that in these 2 years body condition was more strongly related to the most recently experienced rainfall conditions (i.e. at final wintering grounds) rather than those experienced during staging. See Additional file 2 for all supporting data.
Recent research has found that migratory songbirds wintering in sub-Saharan Africa frequently use long-term staging sites before reaching their final destination winter sites. Our results support the hypothesis that local conditions during staging carry over to influence individual condition on the final wintering grounds, but only in years when rainfall conditions are poor. We found that in a drier-than-average year on staging sites (2011), birds had higher CORTf levels and δ 15N values, respectively reflecting exposure to stronger (or more pervasive) stressors and higher aridity, than birds staging in a year with above-average rainfall (2012; Fig. 3). This indicated that rainfall during staging did indeed affect the local conditions experienced by the population. We then found that individual variation in staging conditions experienced by birds (estimated by δ 13C and δ 15N) in 2011 carried over to affect their winter condition in Zambia, as estimated by scaled mass index and BUTY, both of which are likely to influence their overwinter survival and subsequent success on northwards migration. However, in a wetter-than-average year (2012), such individual variation in staging conditions was unrelated to over-winter condition in Zambia, suggesting that when conditions were good, all individuals had access to high quality habitats. Our results demonstrate that differences in local long-term staging conditions may be important for species demography . Moreover, our findings emphasise the challenges in conserving multi-site migration routes and understanding how annual variation in local conditions determine the strength of carry-over effects.
The data support the predictions of our hypothesis; however, they are based on a comparison between 2 years and additional years of study are needed to confirm the generality of the pattern they suggest. It is possible that other differences between years, in addition to rainfall, may have contributed to our observed results. However, rainfall is the most plausible driver given that its connection to food supply is known to be an important environmental driver for the success of wintering migratory birds [18, 43].
Can good conditions on final wintering sites buffer the effects of poor staging conditions?
Despite having experienced poorer habitat conditions during staging, birds wintering in Zambia in 2011 had higher scaled mass indices than in 2012. One explanation may be that local rainfall conditions at the winter study site were wetter (+141 mm) in 2011 than 2012, and that this difference was enough to counteract drier-than-average staging conditions. Alternatively, poor staging conditions in 2011 might have resulted in differential survival . If only the highest quality birds were able to complete the second stage of migration to final winter sites, this would cause scaled mass index to be higher across the winter population despite the presence of carry-over effects from poor staging conditions. Either way, if rainfall underlies these patterns then good rainfall conditions in Zambia appear to be sufficient to counteract drier-than-average staging conditions, at least in the short term . Whether and how carry-over effects from staging during dry years might affect long-term indicators of success (such as spring migration departure dates and success, survival, and fitness) is an important question for determining the consequences of staging site conditions, and should be a focus of future research exploiting continual improvements in tracking technology.
Migratory staging and connectivity
Determining the location of long-term staging sites is essential for understanding how staging habitat conditions may contribute to population dynamics. In this study, we used a triple isotope (δ 15N, δ 13C, δ 2H) cluster analysis to assign birds to staging locations. While this method can assign migrants to large clusters in Africa, it lacks the specificity that can be obtained from geolocator studies. Future work that combines geolocator and isotope data will help to improve the applicability of isotopic assignment in Africa. For species with strong migratory connectivity, determining specific staging sites may be essential for determining population responses to conditions during staging; for example, extreme local staging conditions may have specific and dramatic consequences later in the annual cycle. By contrast, for species with weak migratory connectivity, those effects will be diffuse across the species’ range . Great reed warblers seem likely to experience diffuse effects, since breeding and wintering populations appear to be loosely connected . The discovery of multiple long-term staging sites in many species makes studies of migratory connectivity more difficult, since the degree of connectivity with breeding areas could change over the non-breeding season as birds shift locations. Determining how staging locations are connected to both breeding and wintering grounds is important for predicting the consequences of poor rainfall or drought conditions.
This study supports the hypothesis that habitat conditions on distant long-term staging sites can influence individual condition at final wintering sites, and suggests that these effects are likely to be strongest when conditions are relatively poor. Longitudinal studies are needed to confirm whether the latter pattern is general. Our results suggest that when conditions on final wintering sites are good, individuals may be able to buffer the effects of poor staging conditions in the short term; however, differential survival due to poor staging conditions may also account for this pattern. Research determining variation in the strength of carry-over effects between years has been limited (but see [18, 45, 46]), and to our knowledge this study is the first to assess carry-over effects between recently-discovered long-term staging sites and final wintering sites. Long-term staging sites add complexity to our understanding of population dynamics for declining migratory species especially in the context of global change, in which habitat loss (e.g. owing to agriculture) and unpredictable events outside the reproductive period such as food shortages (e.g. owing to drought) are likely to increase. Our results provide further incentive for logistically challenging endeavours that aim to quantify year-round habitat use and disentangle interactions between multiple environmental factors and sites.
Study site and field data
Fieldwork was conducted in southern Zambia on and around Muckleneuk Farm, near Choma, in ca. 280 ha of mesic grass and reed-bed habitat around each of two watercourses and associated dams (site 1: 16°39′S, 27°00′E; site 2: 16°37′S, 26°59′E), from January to April 2011 and 2012. The surrounding habitat is thornbush savannah and broad-leaved woodland with areas of tobacco, maize, and wheat cultivation. Great reed warblers are widely distributed across their European breeding grounds and sub-Saharan wintering grounds, and absent only from the southwestern half of South Africa and parts of Botswana and Namibia .
Great reed warblers were captured using 8–10 mist nets erected daily from 0530 to 1100 h. Birds were captured without the use of playbacks in order to avoid biases towards the most aggressive individuals. At the time of capture, all individuals were weighed with an electronic scale (± 0.01 g), and tarsus length, bill length, and wing length were measured with digital calipers (± 0.1 mm). To estimate body condition, we calculated the scaled mass index from our morphometric data. This approach scales the mass of all individuals to that expected if they were all of identical body size, and is preferable to the more traditionally used residuals of mass on length . With this approach, a single measure of size is recommended rather than a principal component analysis (PCA) of multiple measures , so we used wing length as our measure of structural size since it was correlated most strongly with weight (wing: r = 0.46, p < 0.001; bill: r = 0.35, p = 0.004; tarsus: r = −0.03, p = 0.85, n = 61). Within 10 min of capture, a 100 μl blood sample was taken for molecular sexing and metabolite analyses (see “Plasma metabolites” section below); thereafter, 5–8 contour feathers (collected from the upper breast) and the 5th right rectrix feather were sampled for isotope and CORT analyses (see below). Because great reed warblers are sexually monomorphic, sex was determined using a molecular technique. Primers for molecular sexing were modified from the primers 2550F and 2718R first described by ; see  for details. Juveniles (first-winter birds) may initiate staging moult on average 15 days later than adults (51 day total moult duration; ), which could cause differences in isotope signatures of feathers; however, this effect is likely to be minimal given the total length of moult. Juveniles were separated from adults by using tongue spot distinctiveness, eye, and tarsus colouration (for details see ). Indeed, we found no difference in isotope signatures between adults and juveniles in our data (δ 15N: t test t′24.1 = 0.69, p = 0.5; δ 13C: t test t′23.8 = 0.13, p = 0.9). Therefore, all individuals (2011 = 7 juveniles, 19 adults; 2012 = 9 juveniles, 28 adults) were pooled for analysis.
Figure 2 shows precipitation during migratory staging (September–November averages) in the northern-hemisphere autumn of 2010 (hereafter ‘2011’ as this is when sampling took place in Zambia) and 2011 (hereafter ‘2012’). To quantify differences between years, we used 3 month precipitation anomalies from the NOAA Climate Prediction Center’s CAMS_OPI dataset (http://www.ncdc.noaa.gov/cag/mapping/global; Accessed May 15 2014; ), which merges observations from rain gauges with precipitation estimates from a satellite algorithm. Anomalies are expressed as the deviation from the 1979–2000 mean value. September–November averages are presented in order to include rainfall just prior to autumn arrival on staging sites (end of September to early October; , since rainfall during this time is likely to be important given the time lag observed between precipitation and primary productivity . We used farm records (Muckleneuk Farm, 16°38′S, 27°00′E) to assess local rainfall at the final winter site in Zambia; rainfall was recorded daily in mm using a rain gauge (from December to February 2011 and 2012).
Carbon and nitrogen
Prior to all isotope analyses, feathers were washed in 2:1 chloroform:methanol solution for 24 h then rinsed with distilled water and left to air dry for 24 h. Sub-samples of 0.8–0.9 mg were weighed into tin capsules. Samples were combusted in a Pyrocube (Elementar, Hanau) elemental analyser. The resulting CO2 and N2 were separated by gas chromatography and admitted into the inlet of a Micromass (Manchester, UK) Isoprime isotope ratio mass spectrometer (IRMS) for determination of 13C⁄ 12C and 15N⁄ 14N ratios. Measurements are reported in δ-notation (δ 13C and δ 15N, respectively) relative to the Pee Dee Belemnite (PDB) for carbon and atmospheric N2 for nitrogen in parts per thousand deviations (‰) using the formula δ (‰) = 1000 × [Rsample/Rstandard−1]. Two sulfanilamides (Isoprime internal standards) and two Casein standards (in house standard) were used for every seven unknowns in sequence. Internal laboratory standards indicated measurement errors (SD) of ±0.05 ‰ for δ 13C, 0.12 ‰ for δ 15N.
Cleaned samples of approximately 0.2 mg were weighed into silver capsules and pyrolysed in an Elementar (Hanau, Germany) Pyrocube elemental analyser over glassy carbon. The resulting H2 was admitted into the source inlet of a Thermo (Bremen, Germany) XP Plus mass spectrometer. Measurements are reported in δ-notation relative to the international standard SMOW (Standard Mean Ocean Water). Organic materials involving H not bonded to carbon will readily exchange H with water vapour . Subtraction of this exchangeable hydrogen is attained using standards with known non-exchangeable (i.e. indigenous) hydrogen isotope compositions. These are CFS (chicken feathers), BWB-II (bowhead whale baleen), and ISB (black kittiwake feathers; for CFS and BWB-II see ; for ISB see . Replicate measurements of these standards imply measurement errors (SD) of around 2 ‰ for δ2H.
Triple-isotope cluster analysis
We used the isotope signature of great reed warbler tail feathers, grown during staging, to determine the most likely staging locations in Africa. Specifically, we used a three-isotope (δ 2H, δ 13C, δ 15N) cluster analysis developed by  to assign birds to one of four clusters in Africa. This method uses isoscapes (i.e. maps of isotopic variation produced by iteratively applying models across regions of space) to estimate the origin of moulted feathers . The clusters were developed by examining previously published plant δ 13C and δ 15N [56, 57] isoscapes and the long-term hydrologic δ 2H  isoscapes of Africa, to obtain information on geographic structure in order to find “natural” groupings in multivariate space (Fig. 2; see  for details). Linear discriminant function analysis was then used to derive algorithms that predict the posterior probability that a sample with a given multi-isotope composition could have originated from any given cluster within Africa, given the predicted ranges for feather δ 2H, δ 13C, and δ 15N (see  for details). Error estimates were derived from . All analyses were conducted within R version 2.15.0 (R Development Core Team 2012) applying algorithms derived from  by means of script provided by Steven Van Wilgenburg.
Blood samples were collected into heparinized capillary tubes from the brachial vein and were centrifuged for 10 min (maximum speed of 8000 rpm) within 4 h of collection. The separated plasma was frozen at −20 °C until analysis. Samples were transported on ice to the Schweizerische Vogelwarte, Sempach, Switzerland, for analyses in May 2012. All metabolites were determined in the plasma using standard test-combinations, the Wako Auto-Kit for 3-hydroxybutyrate (cyclic enzymatic method; BUTY), and the enzymatic colorimetric test for triglycerides (TRIG) including free glycerol (Invicon HIT 917, PAP-method). Dilution curves of the great reed warbler plasma were run for both metabolites to test for linearity. The triglyceride assay was adapted to small volumes (5–10 µl). All samples were measured in duplicates. TRIG and BUTY values were linearly negatively related (r2 = 0.06, p = 0.05); this suggests that the inclusion of free glycerol in our test of triglycerides did not obscure true triglyceride values .
Corticosterone (CORT) was extracted from feathers using a methanol-based procedure  that has been replicated successfully in passerines [60, 61]. A single tail feather per bird was collected and CORT thus extracted. Because multiple contour feathers were collected from all individuals, two contour feathers per bird were extracted to ensure that CORTf measurements were well within limits of detectablility . Following removal of the calamus, the length of each feather sample (proximal cut to distal tip of vane) was measured flat against the edge of a ruler. Feather samples were then placed in glass vials and cut into tiny pieces (c. 5 mm2). Ten mL of methanol (HPLC grade; Fisher Scientific, Fair Lawn, NJ, USA) was added to each vial and samples were sonicated at room temperature for 30 min, followed by overnight incubation at 50 °C. Using a glass funnel fitted with polyester fibre, vacuum filtration was employed to separate the methanol extract from the feather pieces. Collected methanol extracts were evaporated and residues were reconstituted in 600 µL of phosphate-buffered saline (0.05 M, pH 7.6) and frozen at −20 °C until analysed by radioimmunoassay (RIA). Samples were extracted in two batches and the recovery efficiency of each extraction was assessed by spiking three sample extracts with approximately 5000 c.p.m. of 3H-labelled CORT. The average recoveries were 98.6 and 93.2 %, and final values were adjusted for recoveries. Serial dilutions of feather extracts were parallel to the CORT standard curve. Samples were analysed in duplicate in two RIAs using a commercial antiserum (Sigma-Aldrich, St. Louis, MO, USA; product #C8784). The variability of each assay was assessed using six replicates of the same internal standard; average intra-assay coefficient of variation (CV) was 7.4 (SD = 1.1) % and inter-assay CV was 5.5 %. Average detection limit (%B/B0) was 11.0 (SD = 1.5) pg CORT/100 µl of extract and all samples were above detection limits. CORTf measurements were corrected for total feather sample length, reflecting the time-dependent deposition of CORT into feathers, and were thus expressed as pg CORT/mm [31, 59, 63]. CORTf values are likely determined by the interaction between plasma CORT values and the dynamics of feather growth . Although we cannot be sure that the higher CORTf values during staging in 2011 vs. 2012 were due to decreased feather growth rates, higher plasma CORT levels, or both, responses to poor environmental conditions would likely result in higher CORTf values in all of these scenarios , though future work is needed to fully understand the mechanism in our study species.
We analysed each response variable separately (scaled mass index, TRIG, BUTY, final winter site CORTf). We used Akaike’s Information Criterion corrected for small sample sizes (AICc; ) to select from the following predictors δ 13C, δ 15N, staging CORTf, and interaction terms between each parameter and year. Because a difference between years was the a priori prediction for our hypothesis, if an interaction term was included in the top model set (within 2 AICc scores) we constructed multivariate general linear models of all top models including interaction terms (Table 2), and performed exploratory analysis by regressing the response variable and parameter of interest separately for each year (see “Results” section; Fig. 4). For TRIG and BUTY models, time at capture was included as a covariate since this has been shown to be related to plasma metabolite levels [65, 66]; in addition, scaled mass index was included as a covariate in order to assess plasma metabolites as a dynamic measure of change in body condition rather than a static measure of current state, such as scaled mass index . All covariates were standardized prior to analyses such that effect sizes were comparable. All models residuals conformed to normality assumptions for parametric statistics; therefore, response variables were not transformed prior to analysis. Model averaging was not conducted, as sample sizes were not large enough for robust analysis .
Three TRIG data points had Cook’s distances substantially higher than the mean (mean = 0.02 ± 0.008 SE; sample 1 = 0.21, sample 2 = 0.18, sample 3 = 0.34; ). All three outlying TRIG concentrations were measured within 2 weeks of the start of spring migration, such that these high TRIG concentrations were likely from individuals that had begun pre-migratory fattening and do not represent baseline winter feeding conditions. These three values were therefore removed from the analyses, but conclusions were unchanged if they were left in the dataset.
Sex was not included in the main models due to low residual degrees of freedom. Therefore, to determine whether carry-over effects differed between females and males, separate general linear models were instead run for each year, with an interaction term between sex and all three predictor variables (δ 15N, δ 13C, and staging CORTf) included.
To test for differences between groups we used unequal variances (Welch’s) t tests; when data did not meet assumptions of normality we ranked the data, following . To test for differences in plasma metabolite values between years we used ANCOVA and included scaled mass index and time at capture as covariates.
- CORTf :
- δ 15N:
- δ 13C:
Newton I. Can conditions experienced during migration limit the population levels of birds? J Ornithol. 2006;147:146–66.
Newton I. Population limitation in migrants. Ibis. 2004;146:197–226.
Harrison XA, Blount JD, Inger R, Norris DR, Bearhop S. Carry-over effects as drivers of fitness differences in animals. J Anim Ecol. 2010;80:4–18.
Gordo O, Brotons L, Ferrer X, Comas P. Do changes in climate patterns in wintering areas affect the timing of the spring arrival of trans-Saharan migrant birds? Glob Change Biol. 2005;11:12–21.
Saino N, Szep T, Romano M, Rubolini D, Spina F, Moller AP. Ecological conditions during winter predict arrival date at the breeding quarters in a trans- Saharan migratory bird. Ecol Lett. 2004;7:21–5.
Angelier F, Tonra CM, Holberton RL, Marra PP. Short-term changes in body condition in relation to habitat and rainfall abundance in American redstarts Setophaga ruticilla during the non-breeding season. J Avian Biol. 2011;42:335–41.
Studds CE, Marra PP. Rainfall-induced changes in food availability modify the spring departure programme of a migratory bird. Proc R Soc B. 2011;278:3437–43.
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.
Newton I. The migration ecology of birds. London: Elsevier; 2010.
Stach R, Jakobsson S, Kullberg C, Fransson T. Geolocators reveal three consecutive wintering areas in the thrush nightingale. Anim Migr. 2012;1:1–7.
Delmore KE, Fox JW, Irwin DE. Dramatic intraspecific differences in migratory routes, stopover sites and wintering areas, revealed using light-level geolocators. Proc R Soc B. 2012;279:4582–9.
Jahn AE, Cueto VR, Fox JW, Husak MS, Kim DH, Landoll DV, et al. Migration timing and wintering areas of three species of flycatchers (Tyrannus) breeding in the Great Plains of North America. Auk. 2013;130:247–57.
Tøttrup AP, Klaassen RHG, Strandberg R, Thorup K, Kristensen MW, Jørgensen PS, et al. The annual cycle of a trans-equatorial Eurasian-African passerine migrant: different spatio-temporal strategies for autumn and spring migration. Proc R Soc B. 2011;279:1008–16.
Stutchbury BJM, Siddiqui R, Applegate K, Hvenegaard GT, Mammenga P, Mickle N, et al. Ecological causes and consequences of intratropical migration in temperate-breeding migratory birds. Am Nat. 2016;188:S28–40.
Vickery JA, Ewing SR, Smith KW, Pain DJ, Bairlein F, Škorpilová J, et al. The decline of Afro-Palaearctic migrants and an assessment of potential causes. Ibis. 2013;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.
Lemke HW, Tarka M, Klaassen R, Akesson M. Annual cycle and migration strategies of a trans-Saharan migratory songbird: a geolocator study in the great reed warbler. PLoS One. 2013;8:e79209.
Studds CE, Marra PP. Linking fluctuations in rainfall to nonbreeding season performance in a long-distance migratory bird, Setophaga ruticilla. Clim Res. 2007;35:115–22.
Hedenström A, Bensch S, Hasselquist D, Lockwood M, Ottosson U. Migration, stopover and moult of the Great Reed Warbler Acrocephalus arundinaceus in Ghana, West Africa. Ibis. 1993;135:177–80.
Yohannes E, Hansson B, Lee RW, Waldenstrom J, Westerdahl H, Akesson M, et al. Isotope signatures in winter moulted feathers predict malaria prevalence in a breeding avian host. Oecologia. 2008;158:299–306.
Norris DR, Marra PP, Kyser TK, Sherry TW, Ratcliffe LM. Tropical winter habitat limits reproductive success on the temperate breeding grounds in a migratory bird. Proc R Soc B. 2004;271:59–64.
Evans KL, Newton J, Mallord JW, Markman S. Stable isotope analysis provides new information on winter habitat use of declining avian migrants that is relevant to their conservation. PLoS One. 2012;7:e34542.
Minagawa M, Wada E. Stepwise enrichment of 15N along food chains: further evidence and the relation between δ15N and animal age. Geochim Cosmochim Acta. 1984;48:1135–40.
Pardo LH, Nadelhoffer KJ. Using nitrogen isotope ratios to assess terrestrial ecosystems at regional and global scales. In: West JB, Bowen GJ, Dawson TE, Tu KP, editors. Isoscapes: understanding movements, pattern and process on Earth through isotope mapping. New York: Springer; 2009. p. 221–49.
Yohannes E, Hobson KA, Pearson DJ. Feather stable-isotope profiles reveal stopover habitat selection and site fidelity in nine migratory species moving through sub-Saharan Africa. J Avian Biol. 2007;38:347–55.
Heaton THE, Vogel JC, von La Chevallerie G, Collett G. Climatic influence on the isotopic composition of bone nitrogen. Nature. 1986;322:822–3.
Farquhar GD, Ehleringer JR, Hubick KT. Carbon isotope discrimination and photosynthesis. Annu Rev Plant Physiol Plant Mol Biol. 1989;40:503–37.
Wingfield JC, Ramenofsky M. Corticosterone and facultative dispersal in response to unpredictable events. Ardea. 1997;85:155–66.
Marra PP, Holberton RL. Corticosterone levels as indicators of habitat quality: effects of habitat segregation in a migratory bird during the non-breeding season. Oecologia. 1998;116:284–92.
Wingfield JC, Ramenofsky M. Hormones and the behavioral ecology of stress. In: Balm PHM, editor. Stress physiology in animals. Sheffield: Sheffield Academic Press; 1999. p. 1–51.
Jenni-Eiermann S, Helfenstein F, Vallat A, Glauser G, Jenni L. Corticosterone: effects on feather quality and deposition into feathers. Methods Ecol Evol. 2015;6:237–46.
Fairhurst GD, Marchant TA, Soos C, Machin KL, Clark RG. Experimental relationships between levels of corticosterone in plasma and feathers in a free- living bird. J Exp Biol. 2013;216:4071–81.
Romero LM, Fairhurst GD. Measuring corticosterone in feathers: strengths, limitations, and suggestions for the future. Comp Biochem Phys A. 2016. doi:10.1016/j.cbpa.2016.05.002.
Still CJ, Berry JA, Collatz GJ, DeFries RS. Global distribution of C3 and C4 vegetation: carbon cycle implications. Global Biogeochem Cycles. 2003;17:1–14.
Guglielmo CG, Cerasale DJ, Eldermire C. A field validation of plasma metabolite profiling to assess refueling performance of migratory birds. Physiol Biochem Zool. 2005;78:116–25.
Jenni-Eiermann S, Jenni L. Plasma metabolite levels predict individual body mass- changes in a small long-distance migrant, the garden warbler. Auk. 1994;111:888–99.
Féry F, Plat L, Melot C, Balasse EO. Role of fat-derived substrates in the regulation of gluconeogenesis during fasting. Am J Physiol. 1996;270:E822–30.
Jenni-Eiermann S, Jenni L. Metabolic responses to flight and fasting in night- migrating passerines. J Comp Physiol B. 1991;161:465–74.
Peig J, Green AJ. New perspectives for estimating body condition from mass/length data: the scaled mass index as an alternative method. Oikos. 2009;118:1883–91.
Hobson KA, Van Wilgenburg SL, Wassenaar LI, Powell RL, Still CJ, Craine JM. A multi-isotope (δ13C, δ15N, δ2H) feather isoscape to assign Afrotropical migrant birds to origins. Ecosphere. 2012;3:art44.
Yohannes E, Lee RW, Jochimsen MC. Stable isotope ratios in winter-grown feathers of great reed warblers Acrocephalus arundinaceus, clamorous reed warblers A. stentoreus and their hybrids in a sympatric breeding population in Kazakhstan. Ibis. 2011;153:502–8.
Janowiak JE, Xie P. CAMS–OPI: a global satellite–rain gauge merged product for real-time precipitation monitoring applications. J Climate. 1999;12:3335–42.
Szép T. Relationship between west African rainfall and the survival of central European sand martins Riparia riparia. Ibis. 2008;137:162–8.
Taylor CM, Norris DR. Population dynamics in migratory networks. Theor Ecol. 2009;3:65–73.
Legagneux P, Fast PLF, Gauthier G, Bety J. Manipulating individual state during migration provides evidence for carry-over effects modulated by environmental conditions. Proc R Soc B. 2012;279:876–83.
Harrison XA, Hodgson D, Inger R, Colhoun K, Gudmundsson G, McElwaine G, et al. Environmental conditions during breeding modify the strength of mass- dependent carry-over effects in a migratory bird. PLoS One. 2013;8:e77783.
Cramp S, Perrins CM. The birds of the Western Palearctic. Oxford: Oxford University Press; 1994.
Fridolfsson A, Ellegren H. A simple and universal method for molecular sexing of non-ratite birds. J Avian Biol. 1999;30:116.
Round PD, Hansson B, Pearson DJ, Kennerley PR, Bensch S. Lost and found: the enigmatic large-billed reed warbler Acrocephalus orinus rediscovered after 139 years. J Avian Biol. 2007;38:133–8.
Bensch S, Hasselquist D, Nielsen B, Hansson B. Higher fitness for philopatric than for immigrant males in a semi-isolated population of great reed warblers. Evolution. 1998.
Sinclair ARE. Factors affecting the food supply and breeding season of resident birds and movements of palaearctic migrants in a tropical african savannah. Ibis. 2008;120:480–97.
Schimmelmann A, Miller RF, Leavitt SW. Hydrogen isotopic exchange and stable isotope ratios in cellulose, wood, chitin, and amino compounds. In: Swart PK, Lohmann KC, McKenzie J, Savin S, editors. Climate change in continental isotopic records. Washington: Geophysical Monograph 68, American Geophysical Union; 1993. p. 367–74.
Hobson KA, Wassenaar LI. Tracking animal migration with stable isotopes. Amsterdam: Elsevier; 2008.
Fox TAD, Christensen TK, Bearhop S, Newton J. Using stable isotope analysis of multiple feather tracts to identify moulting provenance of vagrant birds: a case study of baikal teal Anas formosa in Denmark. Ibis. 2007;142:622–5.
Bowen GJ, Wassenaar LI, Hobson KA. Global application of stable hydrogen and oxygen isotopes to wildlife forensics. Oecologia. 2005;143:337–48.
Still CJ, Powell RL. Continental-scale distributions of vegetation stable carbon isotope ratios. Isoscapes: understanding movements, pattern and process on Earth through isotope mapping. New York: Springer; 2009. p. 179–93.
Craine JM, Elmore AJ, Aidar MPM, Bustamante M, Dawson TE, Hobbie EA, et al. Global patterns of foliar nitrogen isotopes and their relationships with climate, mycorrhizal fungi, foliar nutrient concentrations, and nitrogen availability. New Phytol. 2009;183:980–92.
Hobson KA, Bairlein F. Isotopic fractionation and turnover in captive garden warblers Sylvia borin: implications for delineating dietary and migratory associations in wild passerines. Can J Zool. 2003;81:1630–5.
Bortolotti GR, Marchant TA, Blas J. Corticosterone in feathers is a long-term, integrated measure of avian stress physiology. Funct Ecol. 2008;22:494–500.
Fairhurst GD, Vögeli M, Serrano D, Delgado A, Tella JL, Bortolotti GR. Can synchronizing feather-based measures of corticosterone and stable isotopes help us better understand habitat–physiology relationships? Oecologia. 2013;173:731–43.
Fairhurst GD, Dawson RD, van Oort H, Bortolotti GR. Synchronizing feather- based measures of corticosterone and carotenoid-dependent signals: what relationships do we expect? Oecologia. 2014;174:689–98.
Fairhurst GD, Bond AL, Hobson KA, Ronconi RA. Feather-based measures of stable isotopes and corticosterone reveal a relationship between trophic position and physiology in a pelagic seabird over a 153-year period. Ibis. 2015;157:273–83.
Bortolotti GR. Flaws and pitfalls in the chemical analysis of feathers: bad news- good news for avian chemoecology and toxicology. Ecol Appl. 2010;20:1766–74.
Burnham KP, Anderson DR. Model selection and multimodel inference: a practical information-theoretic approach. New York: Springer; 2002.
Guglielmo CG, O’Hara PD, Williams TD. Extrinsic and intrinsic sources of variation in plasma lipid metabolites of free-living western sandpipers Calidris mauri. Auk. 2002;119:437–45.
Jenni L, Jenni-Eiermann S. Metabolic responses to diurnal feeding patterns during the postbreeding, moulting and migratory periods in passerine birds. Funct Ecol. 1996;10:73–80.
Grueber CE, Nakagawa S, Laws RJ, Jamieson IG. Multimodel inference in ecology and evolution: challenges and solutions. J Evol Biol. 2011;24:699–711.
Fox J. Regression diagnostics: an introduction. Newbury Park: Sage; 1991.
Ruxton GD. The unequal variance t test is an underused alternative to student’s t test and the Mann–Whitney U test. Behav Ecol. 2006;17:688–90.
Rollin BE, Kessel ML. Guidelines for the treatment of animals in behavioural research and teaching. Anim Behav. 2012;83:301–9.
MCS and CNS designed the study, MCS conducted fieldwork, MCS, GDF, SJ-E, JN, and EY conducted lab assays, and MCS wrote the manuscript with contributions from CNS, GDF, SJ-E, JN and EY. All authors read and approved the final manuscript.
We thank Kenty Mudenda, Collins Moya, Jason Boyce, and Jane Jönsson for assistance in the field and lab, Steve Van Wilgenburg and Keith Hobson for providing isotopic cluster analysis code and guidance during analysis, Ian and Emma Bruce-Miller for their generous hospitality and for supplying rainfall data in Zambia, Hamish Ross for access to part of the study area, Staffan Bensch for generous access to lab space, and the Zambia Wildlife Authority for permits. We gratefully acknowledge Tracy Marchant for allowing us to use her endocrine lab at the University of Saskatchewan.
The authors declare that they have no competing interests.
Availability of supporting data
The dataset supporting the conclusions of this article is included as Additional file 2.
Our work complies with the international animal care guidelines of the Association for the Study of Animal Behaviour . The capture and handling of animals were approved by and conducted under licence from ZAWA (Zambia Wildlife Authority).
MCS was funded by the Gates Cambridge Trust and the Natural Sciences and Engineering Research Council of Canada; CNS was supported by a Royal Society Dorothy Hodgkin Fellowship, a BBSRC David Phillips Fellowship (BB/J014109/1), and the DST-NRF Centre of Excellence at the FitzPatrick Institute. Isotope analyses at the LSMSF were funded by a NERC LSMSF Grant (EK206-16/12).