The study of individual movement patterns using stable hydrogen isotopes (δ 2H) as a natural marker has grown; however, recent studies have suggested that measurement of δ 2H in feathers (δ 2Hf) may prove unreliable as a means for identifying region of origin of migrating or wintering birds, particularly raptors. In this study, we examine whether differences in body condition could explain some of the variability in δ 2H in feathers. We analyzed growing feathers of 21 Swainson's Hawks breeding in northern CA for δ 2H, nitrogen (δ 15N), and carbon (δ 13C) stable isotopes in relation to body condition. We found that δ 2H was variable (range = 40‰), and that variability was significantly associated with body condition. Raptors derive most or all of their moisture from prey. Therefore, we suggest that individuals in poor condition have an enriched pool of body water relative to individuals in good condition, due to fractionation of body water stores during respiratory water loss and metabolic processes. Body condition was also negatively correlated with δ 15Nf. However, δ 2Hf, δ 15Nf, and δ 13Cf were not correlated, suggesting that the relationship between δ 2Hf and body condition is a result of physiological processes rather than differences in dietary δ 2H. We used an isotopic basemap of δ 2Hf values to assess individual origin as if they were encountered naively on the migration or wintering grounds, and all individuals fell within the 95% confidence interval of our study area. Conversely, the 95% confidence interval of δ 2Hf values obtained encompassed almost the entire breeding range of this species, indicating little ability to differentiate origins of this species.
The use of naturally occurring stable isotopes within animal tissues as a tool to answer ecological questions has helped explain questions about diet, movements, and physiology of many bird species (Thompson et al. 2005, Fox and Bearhop 2008, Inger and Bearhop 2008, Kelly et al. 2008). The isotopes of carbon (δ 13C), nitrogen (δ 15N), and hydrogen (δ 2H) have been particularly useful in ecological studies (Thompson et al. 2005). In particular, δ 2H is often used to assess movement patterns of individuals because δ 2H varies consistently across the landscape due to variation in rainfall (Hobson and Wassenaar 1997). In fact, δ 2H has been used to investigate a myriad of ecological questions associated with animal movement, such as migratory connectivity (Hobson 2005, Sarasola et al. 2008), migratory behavior (e.g., Cardador et al. 2015, Domenech and Pitz 2015) and natal dispersal (Hobson et al. 2004) in avian species. Birds provide an ideal study organism for such investigations because feathers provide researchers a source of isotopic material that is metabolically inert and retains the isotopic signature of the area in which it was grown (Hobson and Wassenaar 1997).
Despite increased interest and utilization of δ 2H analyses, recent studies have suggested that stable isotope compositions may not be reliably measured within a population or even across an individual feather in some species (Smith et al. 2008, 2009, Wunder et al. 2009). These results indicate a lack of precision when trying to assess origin of tissues using δ 2H or other isotopes. In spite of these questions, some studies advocate grouping observations into regions of origin, rather than providing precise estimates of geographic origin, to account for some of the variability between samples collected in the same location (Langin et al. 2007); however, the source of this interindividual variability is largely unknown. We hypothesize that some of this variation of feather-δ 2H is the result of body condition-dependent processes, particularly for individuals that derive some or all of their water from food.
Other isotopes, such as δ 13C and δ 15N, have been used to examine diet at the population and individual levels (reviewed in Inger and Bearhop 2008). Nitrogen is often used to reveal trophic positions of individuals because nitrogen becomes more enriched with increased trophic level due to fractionation (DeNiro and Epstein 1981, Hobson et al. 1993). However, in addition to trophic enrichment, several studies have demonstrated an enrichment of δ 15N in nutritionally stressed individuals (Ambrose and DeNiro 1986, Hobson and Clark 1992, Hobson et al. 1993, Castillo and Hatch 2007). For example, Hobson et al. (1993) demonstrated that tissues from incubating Ross's Geese (Chen rossii) were enriched after the incubation period (i.e., post-fasting during egg-laying and incubation). The authors suggest that the enrichment is the result of catabolization of tissues that have already been enriched to meet energetic demands. Thus, catabolization causes enrichment of δ 15N values above a baseline. However, not all studies have found δ 15N enrichment correlated with measure(s) of body condition (Ben-David et al. 1999). Similarly, δ 13C may vary by diet, as different photosynthetic pathways differ in their discrimination of the heavier isotope (Fry 2006) and may be depleted when food is restricted in some species (Robb et al. 2015).
For species that acquire most or all of their water from food resources, we predict that δ 2H will also be isotopically enriched. In this study, we examined feathers collected from a population of breeding Swainson's Hawks (Buteo swainsoni) nesting in northern CA U.S.A. By only sampling feathers of nesting adults that were still in sheath, we ensured that all feather material was grown within a few km of the nest site (Woodbridge 1991). Therefore, all feathers should have similar isotopic signatures if the isotopic composition of precipitation (δ 2Hp) is the primary driver of the isotopic composition of feathers (δ 2Hf), as they were all grown within 30 km of each other with only 60 m of elevation change across the study area (Woodbridge 1991). The value of δ 2Hf should be correlated with δ 2Hp, as δ 2H from rainfall is incorporated into plants and subsequently prey. Fractionation can occur at all stages of the process, wherein lighter isotopes may be used preferentially in metabolic processes (e.g., evaporative water loss, respiration, etc.; Fry 2006). Therefore, individuals can become isotopically enriched relative to precipitation, but they should reflect δ 2Hp. In contrast, if other factors (e.g., body condition) play a role in controlling δ 2Hf, then we expect wide range of δ 2Hf values, disrupting the correlation between δ 2Hp and δ 2Hf. If such wide variability exists, δ 2Hf would be a poor surrogate of molt location. Swainson's Hawks provide a good study species because they have relatively small home ranges (Woodbridge 1991) and generally do not drink standing water (Roest 1957, Cooper 1968). Like most birds of prey, they acquire most or all of their body moisture from their prey (Bartholomew and Cade 1957). Therefore, δ 2Hf signatures of feathers grown on the breeding territory should reflect local prey, and should not vary substantially between individuals because there would not be differences in δ 2Hp across territories.
Methods
Study Site and Species. We monitored a population of breeding Swainson's Hawks in Butte Valley, CA U.S.A. (41°45.7′N, 121°48.37′W) from April through August, 2008–2010. We monitored territories annually and located nest sites by watching for nest-building, mating behavior, and territoriality (April–May). In the summer months, we found nests by watching for prey deliveries to the nest site. We trapped adults between June 29 and August 15 near the nest site using dho-gaza-style net with a Great Horned Owl lure (Bubo virginianus, Bloom et al. 1992) or a bal-chatri baited with mice (Berger and Mueller 1959). We color-marked adults using unique two-digit color bands and a U.S. Geological Survey (U.S.G.S.) aluminum band, and determined sex by presence/absence of a brood patch, by size measurements, and by observations of copulatory behavior.
We measured wing chord to the nearest mm and weight to the nearest 1 g. We also opportunistically collected a growing secondary covert (i.e., the feather was still in sheath) from either wing. However, growing secondary coverts were not in the same position on each individual; therefore, we could only collect these feathers if the individual was molting a secondary feather. We recorded whether each individual had any food present in its crop, which could bias our weight measurements. All procedures were approved by the University of Nevada, Reno, IACUC (protocol no. 000115).
Isotopic Analysis. We washed each feather to remove oils following the recommendations of Paritte and Kelly (2009). We collected all feather material from the distal portion of each feather approximately 50 mm from the tip of the feather and we avoided using the central rachis to avoid potential contamination with blood. We weighed out between 500–600 μg of feather material for δ 2H analysis, and analyzed it following the technique of Hilkert et al. (1999), using δ 2H standards obtained from L.I. Wassenaar (Environment Canada) to adjust for the exchangeable portion of hydrogen in keratin (Wassenaar and Hobson 2003). Results are reported in standard δ notation in units of ‰ versus VSMOW.
We weighed out an additional 500–600 μg of feather material for δ 13C and δ 15N analysis, and analyzed it following the technique of Werner et al. (1999). Results are reported in standard δ notation in units of ‰ versus VPDB and air, respectively. All samples were analyzed at the Nevada Stable Isotope Laboratory in Reno, Nevada U.S.A.
Statistical Analysis. We estimated an index of body condition by using a standard major axis (SMA) regression to predict body mass based on wing chord (Peig and Green 2009). We regressed the natural log of weight against the natural log of wing chord for all birds of one sex and obtained an index of condition for each individual by subtracting the mass predicted by the SMA regression from the observed weight. The difference between the actual weight of the individual and the predicted weight was considered the index of condition (i.e., a positive value indicates a heavier weight than would be expected for a given size of the individual). We calculated condition indices separately for males and females. We used Pearson correlations to examine potential relationships among δ 2H, carbon (δ 13C), and nitrogen (δ 15N).
We used a linear mixed model to perform a series of three regressions examining δ 2Hf, δ 13Cf, and δ 15Nf and their relationships with (1) the residuals of the SMA regression (i.e., condition index), (2) capture date (Julian date), and (3) sex in R 3.3.1 (R Development Core Team 2009) and the nlme package (Pinheiro et al. 2016). Specifically, we regressed each of the isotopes against body condition and capture date as fixed effects and year was used as a random effect for all models to account for potential unmeasured annual differences. For all analyses, no individuals showed any signs of food in the crop. Due to low sample sizes, we did all regressions univariately to avoid overparameterization of the model. We set our significant threshold to α = 0.05 for all models.
Geographic Assignment. We followed Hobson et al. (2009) to create a basemap of δ 2H values for raptors of North America to determine where each individual would be naively located if it were encountered outside of the breeding range (e.g., Sarasola et al. 2008). Because there were too few Swainson's Hawks measured in the Lott and Smith (2006) study to create a basemap specific for Swainson's Hawks (Hobson et al. 2009), we created a basemap using all species measured in Lott and Smith (2006). Following Hobson et al. (2009) we used a reduced major axis regression to determine the average fractionation of δ 2Hf from δ 2Hp, where δ 2Hp was calculated for the growing season (Meehan et al. 2004). We used the SD of the estimate from the SMA regression to create 1000 new estimates of δ 2Hf. We created percentiles (5th through 95th) based on those simulated results and used those estimates create a map of potential origins. We used these measures to assess the accuracy of using δ 2H as a predictor in describing region of feather growth.
Results
The estimated growing-season δ 2Hp within the study area was −102‰, and δ 2Hf values ranged from −71‰ to −114‰. There were no correlations among δ 2Hf and δ 15Nf, or δ 13Cf isotopic compositions (P > 0.24 for all comparisons). There was a significant relationship between δ 2Hf and body condition (P < 0.001, n = 21) as well as δ 15N and condition (P < 0.05, n = 21; Table 1; Fig. 1). All other relationships were not significant (P > 0.1).
Table 1.
Model results (estimates ± SE) from a linear mixed model of hydrogen (δ 2H), nitrogen (δ 15N), and carbon (δ 13C) stable isotopes in breeding Swainson's Hawk feathers from Butte Valley, CA, from 2008–2010 against body condition, capture date, and sex.
We calculated that the mean δ 2Hf value for our study area should be −106‰ (n = 21). The 95% CI for our error-propagated results for Butte Valley, CA was −147 to −64‰. Despite the variability in δ 2Hf composition, all individuals sampled fell within the 95% confidence interval based on the residuals from the SMA regression of δ 2Hf, and the mean δ 2Hf observed from growing Swainson's Hawk feathers (−95 ± 3‰) were within the 50th percentile of predicted δ 2Hf values for our study area (Fig. 2).
Discussion
We found isotopic enrichment of both δ 2Hf and δ 15Nf in Swainson's Hawks negatively related to body condition of breeding individuals. Isotopic enrichment of δ 15N in tissues has been described in several species (Ambrose and DeNiro 1986, Hobson and Clark 1992, Hobson et al. 1993, but see Ben-David et al. 1999). For example, Hobson et al. (1993) found that juvenile Japanese Quail (Coturnix japonica) that were food-deprived had more enriched δ 15N values compared to individuals fed ad lib. This relationship was also observed in Swainson's Hawks, with individuals in relatively poorer body condition having enriched feathers (Bearhop et al. 2002). The underlying cause of the relationship is still unknown; however, it is reasonable to suspect that the catabolization of tissues to meet energetic demands, and the subsequent fractionation as those tissues are used to fuel the individual, result in more enriched tissues. These enriched tissues are then incorporated into newly growing feathers, obscuring the relationship between δ 2Hf and δ 2Hp.
However, a similar relationship among body condition and isotopic enrichment has not been observed previously in δ 2H tissue values. Raptors derive most or all of their body water from prey, and are rarely observed to drink water (Bartholomew and Cade 1957). Individuals in poor condition (i.e., individuals that do not eat or do not eat enough to cover energetic costs) will therefore have a net water loss due to evaporative, respiratory, and metabolic water losses. Evaporative and respiratory water is expected to be isotopically light versus body water, due to the isotopic fractionation associated with the water liquid to water vapor phase change (Horita and Wesolowski 1994). Hence, without additional water input (in the form of prey intake), isotopic enrichment of remaining body water would occur, and this isotopically heavy signature would subsequently be recorded in the δ 2Hf composition. In fact, this is a mechanism proposed to potentially cause difference in δ 2Hf between nestling and adult American Kestrels (Falco sparverius; Greenwood and Dawson 2011) and Northern Saw-whet Owls (Aegolius acadicus; De Ruyck et al. 2013). The hypothesis of reduced intake of prey may also be consistent with the observed isotopic enrichment of δ 15Nf at low body condition, due to the nitrogen isotope fractionation associated with the deamination of amino acids (Macko et al. 1986), and depending upon the ratio of dietary nitrogen lost via excretion to dietary nitrogen uptake (Fry 2006).
It is possible that the relationship between body condition and δ 2H may explain some of the variation observed in raptor feathers in other studies (e.g., Smith et al. 2008, 2009). Variation in δ 2H within raptors has been observed in a number of species, and has been documented at >40‰ (e.g., Smith et al. 2008), which is similar to the variation among individuals observed in this population (43‰). Although such variation may be typical of raptors, it confounds our ability to accurately assess origins of feather growth in wintering or migratory species, even in regions where there is significant isotopic variation in precipitation. Although all samples fell within the 95% confidence interval for the study area, the 95% confidence interval also encompassed almost the entire breeding range of this species, indicating little ability to differentiate origins of this species. The only area excluded from our 95% confidence interval was the extreme northern part of this species' breeding range. Thus, the wide confidence intervals of the model, coupled with large intrapopulation variability, suggest that studies examining migratory connectivity of this species (e.g., Sarasola et al. 2008) may need more validation before definitive conclusions can be drawn, regardless of the mechanism underlying the intrapopulation differences in δ 2H.
There were no significant correlations among any of the isotopes we measured (i.e., δ 2Hf, δ 15Nf, and δ 13Cf). A lack of correlations among δ 2Hf and δ 15Nf, and δ 13Cf indicates that the correlation between δ 2Hf and body condition may not be mediated by differences in prey base or trophic position of breeding individuals. Thus, the isotopic enrichment of δ 2Hf, and likely other tissues (Bearhop et al. 2002) in Swainson's Hawks, is more likely due to water loss without replacement if an individual does not eat, rather than isotopic differences or fractionation rates of differing prey items. Similarly, enrichment of δ 15Nf may come from fractionation as tissues are catabolized.
Almost 90% of the diet of Swainson's Hawks in this area consists of small mammals, primarily Belding's ground squirrels (Spermophilus beldingi), Mazama pocket gophers (Thomomys mazama), and montane voles (Microtus montanus; Woodbridge 1991). Therefore, prey should also reflect the local δ 2H values, as the species that are consumed in large numbers either are nonmigratory, or have not yet had time to migrate (i.e., young of the year or juveniles). However, because we did not sample prey, we could not explicitly rule out prey differences driving relationships among body condition and δ 2Hf or δ 15Nf. Similarly, we could not exclude differential use of exogenous versus endogenous reserves within the individual (Oppel et al. 2010). In fact, Swainson's Hawks may refuel on their pre-breeding migration (Bechard et al. 2006, Kochert et al. 2011), which could provide outside δ 2H or δ 15N for fuel later in the season. Therefore, some variability in δ 2Hf or δ 15Nf values could reflect catabolization of body reserves collected months prior.
In addition, differences in diet and prolonged dietary restriction may reduce metabolic rate (Cherel et al. 1988), which can influence isotope fractionation factors and the observed δ 2H values (Wassenaar and Hobson 2006). For example, individuals eating prey with higher fat content require more water per caloric unit (Kirkley and Gessaman 1990). Isotopic fractionation of body water for metabolic processes associated with processing high fat content could lead to enriched δ 2H of available water stores within an individual. For populations that solely or primarily acquire water through their diet and whose location during molt is known, δ 2H may provide a longer-term measure of body condition relative to mass adjusted for weight. This may be particularly beneficial in species whose ptilochronology cannot be assessed (Grubb 2006). However, more research is necessary on the underlying physiological mechanisms that cause increased δ 2H compositions in individuals with low body weight relative to size (e.g., metabolic versus evaporative water loss and fractionation).
Acknowledgments
We thank the private landowners of Butte Valley for access to nests; B. Smucker, J. Barnes, and C. Cheyne for field assistance; and L. Wassenaar for isotope standards for δ 2H. We thank M. Ben-Hamo, C. Downs, B. Pinshow, and three anonymous reviewers for comments on earlier drafts of this report. This work was conducted under federal Bird Banding lab permit number 21368 and California Scientific Collecting permit 007333.