Linking hydrogen (δ2H) isotopes in feathers and precipitation: sources of variance and consequences for assignment to isoscapes

A trans-oceanic flight of over 4,200 km by painted lady butterflies

The extent of aerial flows of insects circulating around the planet and their impact on ecosystems and biogeography remain enigmatic because of methodological challenges. Here we report a transatlantic crossing by Vanessa cardui butterflies spanning at least 4200 km, from West Africa to South America (French Guiana) and lasting between 5 and 8 days. Even more, we infer a likely natal origin for these individuals in Western Europe, and the journey Europe-Africa-South America could expand to 7000 km or more. This discovery was possible through an integrative approach, including coastal field surveys, wind trajectory modelling, genomics, pollen metabarcoding, ecological niche modelling, and multi-isotope geolocation of natal origins. The overall journey, which was energetically feasible only if assisted by winds, is among the longest documented for individual insects, and potentially the first verified transatlantic crossing. Our findings suggest that we may be underestimating transoceanic dispersal in insects and highlight the importance of aerial highways connecting continents by trade winds.

Winter-ground microhabitat use by differently coloured phenotypes affects return rate in a long-distance migratory bird

Migratory bird populations are declining globally at alarming rates. Non-breeding site conditions affect breeding populations, but generalising non-breeding habitat conditions over large spatial regions cannot address potential fine-scale differences across landscapes or local populations. Plumage characteristics can mediate the effects of environmental conditions on individual fitness. However, whether different phenotypes use distinctive non-breeding sites, and whether they respond to non-breeding site conditions differently remains largely unknown. Stable isotopes (δ13C, δ15N, δ2H) of inert tissues are useful to infer habitat characteristics and geographic origins where those tissues were grown. We collected winter-grown feathers from pied flycatchers (Ficedula hypoleuca) on their breeding grounds over several years from males whose dorsal plumage colouration ranged continuously from brown to black and assessed their stable isotope values as proxies of local habitat conditions. Based on feather δ2H profiles we found that browner males spent their non-breeding season in drier habitats than black males. Assignment to origin analysis shows potential regional non-breeding ground separation between differently coloured males. High within-individual repeatability of both δ13C and δ15N indicate the pied flycatcher males return yearly to similar areas. Blacker males were more likely to return to the breeding grounds after dry years compared with brown males. The opposite was found in wet years. Our study demonstrates that different phenotypes are exposed to different non-breeding site conditions which can differentially affect individual survivorship. This has important ramifications for population dynamics under predicted climate change scenarios where especially brown phenotype pied flycatcher males may be under a risk of decreasing.

The geographic extent of bird populations affected by renewable-energy development

Bird populations are declining globally. Wind and solar energy can reduce emissions of fossil fuels that drive anthropogenic climate change, yet renewable-energy production represents a potential threat to bird species. Surveys to assess potential effects at renewable-energy facilities are exclusively local, and the geographic extent encompassed by birds killed at these facilities is largely unknown, which creates challenges for minimizing and mitigating the population-level and cumulative effects of these fatalities. We performed geospatial analyses of stable hydrogen isotope data obtained from feathers of 871 individuals of 24 bird species found dead at solar- and wind-energy facilities in California (USA). Most species had individuals with a mix of origins, ranging from 23% to 98% nonlocal. Mean minimum distances to areas of likely origin for nonlocal individuals were as close as 97 to >1250 km, and these minimum distances were larger for species found at solar-energy facilities in deserts than at wind-energy facilities in grasslands (Cohen's d = 6.5). Fatalities were drawn from an estimated 30-100% of species' desingated ranges, and this percentage was significantly smaller for species with large ranges found at wind facilities (Pearson's r = -0.67). Temporal patterns in the geographic origin of fatalities suggested that migratory movements and nonmigratory movements, such as dispersal and nomadism, influence exposure to fatality risk for these birds. Our results illustrate the power of using stable isotope data to assess the geographic extent of renewable-energy fatalities on birds. As the buildout of renewable-energy facilities continues, accurate assessment of the geographic footprint of wildlife fatalities can be used to inform compensatory mitigation for their population-level and cumulative effects.

Determinants of spring migration departure dates in a New World sparrow: Weather variables reign supreme

Numerous factors influence the timing of spring migration in birds, yet the relative importance of intrinsic and extrinsic variables on migration initiation remains unclear. To test for interactions among weather, migration distance, parasitism, and physiology in determining spring departure date, we used the Dark-eyed Junco (Junco hyemalis) as a model migratory species known to harbor diverse and common haemosporidian parasites. Prior to spring migration departure from their wintering grounds in Indiana, USA, we quantified the intrinsic variables of fat, body condition (i.e., mass ~ tarsus residuals), physiological stress (i.e., ratio of heterophils to lymphocytes), cellular immunity (i.e., leukocyte composition and total count), migration distance (i.e., distance to the breeding grounds) using stable isotopes of hydrogen from feathers, and haemosporidian parasite intensity. We then attached nanotags to determine the timing of spring migration departure date using the Motus Wildlife Tracking System. We used additive Cox proportional hazard mixed models to test how risk of spring migratory departure was predicted by the combined intrinsic measures, along with meteorological predictors on the evening of departure (i.e., average wind speed and direction, relative humidity, and temperature). Model comparisons found that the best predictor of spring departure date was average nightly wind direction and a principal component combining relative humidity and temperature. Juncos were more likely to depart for spring migration on nights with largely southwestern winds and on warmer and drier evenings (relative to cooler and more humid evenings). Our results indicate that weather conditions at take-off are more critical to departure decisions than the measured physiological and parasitism variables.

Multi-isotope (δ2H, δ13C, δ15N) feather profiles and morphometrics inform patterns of migratory connectivity in three species of North American swallows

Aerial insectivorous birds have suffered steep population declines in North America over the last 60 years. A lack of information on migratory connectivity between breeding and non-breeding grounds for these species limits our ability to interpret factors affecting their population-specific trends. We determined likely Latin American non-breeding regions of Bank (Riparia riparia), Barn (Hirundo rustica) and Cliff (Petrochelidon pyrrhonota) swallow from populations across their breeding ranges. We used predicted feather hydrogen (δ²H_f) and carbon (δ¹³C_f) isoscapes for winter-grown feathers to indicate areas of highest probability of moult origin and incorporated these results into a cluster analysis to determine likely broad non-breeding regions. We also assessed variation in wing length among populations to determine the potential for this metric to differentiate population moult origins. We then investigated patterns of multi-isotopic (δ²H_f, δ¹³C_f, δ¹⁵N_f) and wing-length niche occupancy by quantifying niche size and overlap among populations under the assumption that broad niches were consistent with low within-species migratory connectivity and narrow and non-overlapping niches with higher connectivity. Multivariate assignment identified different non-breeding regions and potential clusters of moult origin generally corresponding to Central America and northern South America, eastern and south-central South America, and the western and southern part of that continent, with variation within and among populations and species. Separate niche space indicated different wintering habitat or areas used by species or populations whereas niche overlap indicated only potential spatial similarity. Wing length varied significantly among populations by species, being longer in the west and north for Bank and Cliff Swallow and longer in eastern Canadian Barn Swallow populations. Barn Swallow occupied consistently larger isotopic and wing length niche space than the other species. Comparisons among populations across species showed variable isotopic and wing-length niche overlap generally being greater within breeding regions and lower between western and eastern breeding populations supporting a general North American continental divide for all species with generally low migratory connectivity for all species. We present a novel approach to assessing connectivity using inexpensive and broad isotopic approaches that provides the basis for hypothesis testing using more spatially explicit expensive techniques.

Assigning harvested waterfowl to geographic origin using feather δ2H isoscapes: What is the best analytical approach?

Establishing links between breeding, stopover, and wintering sites for migratory species is important for their effective conservation and management. Isotopic assignment methods used to create these connections rely on the use of predictable, established relationships between the isotopic composition of environmental hydrogen and that of the non-exchangeable hydrogen in animal tissues, often in the form of a calibration equation relating feather (δ2Hf) values derived from known-origin individuals and amount-weighted long-term precipitation (δ2Hp) data. The efficacy of assigning waterfowl to moult origin using stable isotopes depends on the accuracy of these relationships and their statistical uncertainty. Most current calibrations for terrestrial species in North America are done using amount-weighted mean growing-season δ2Hp values, but the calibration relationship is less clear for aquatic and semi-aquatic species. Our objective was to critically evaluate current methods used to calibrate δ2Hp isoscapes to predicted δ2Hf values for waterfowl. Specifically, we evaluated the strength of the relationships between δ2Hp values from three commonly used isoscapes and known-origin δ2Hf values three published datasets and one collected as part of this study, also grouping these data into foraging guilds (dabbling vs diving ducks). We then evaluated the performance of assignments using these calibrations by applying a cross-validation procedure. It remains unclear if any of the tested δ2Hp isoscapes better predict surface water inputs into food webs for foraging waterfowl. We found only marginal differences in the performance of the tested known-origin datasets, where the combined foraging-guild-specific datasets showed lower assignment precision and model fit compared to data for individual species. We recommend the use of the more conservative combined foraging-guild-specific datasets to assign geographic origin for all dabbling duck species. Refining these relationships is important for improved waterfowl management and contributes to a better understanding of the limitations of assignment methods when using the isotope approach.

Examination of amino acid hydrogen isotope measurements of scalp hair for region-of-origin studies

RATIONALE

Hydrogen isotope (δ2 H) analysis of keratinaceous bulk tissues has been used in forensic science to reconstruct an individual's travel history or determine their region-of-origin. Here, we use a compound-specific approach to examine patterns of individual amino acid δ2 H values in relation to those of local tap water, bulk scalp hair tissues, and region-of-origin.

METHODS

We measured δ2 H values of amino acids in anonymously collected scalp hair (n = 67) and tap water from 28 locations in the United States. Samples were hydrolyzed into their constituent amino acids, derivatized alongside in-house reference materials, and analyzed in triplicate using a GC-C-IRMS system.

RESULTS

Non-essential amino acid (AANESS ) δ2 H values and their corresponding tap water samples varied systematically across continental regions. Hydrogen isotope values of alanine, glutamic acid, and glycine were significantly correlated with tap water and an estimated 42%-51% of the hydrogen atoms in these AANESS were derived from tap water. We used linear discriminate analysis (LDA) to explore regional patterns in scalp hair bulk tissue and amino acid δ2 H values. For the model that included AANESS data, 87% of the variance was explained by the first linear discriminant axis (LD1), and was driven by bulk hair tissue, alanine, and proline. This model had an overall 72% successful reclassification with samples from the south and northwest regions reclassifying correctly 92% and 78% of the time, respectively. For the model that included AAESS data, LD1 explained 81% of the variation and was driven bulk hair, threonine, valine, phenylalanine, and isoleucine. The overall reclassification rate for the model that included AAESS was 70%.

CONCLUSIONS

Our findings suggest that δ2 H analyses of AANESS and AAESS could help improve geolocation models for human and wildlife forensics by simultaneously providing information about both dietary and tap water inputs of hydrogen to tissue synthesis.

Population genetics and geographic origins of mallards harvested in northwestern Ohio

The genetic composition of mallards in eastern North America has been changed by release of domestically-raised, game-farm mallards to supplement wild populations for hunting. We sampled 296 hatch-year mallards harvested in northwestern Ohio, October-December 2019. The aim was to determine their genetic ancestry and geographic origin to understand the geographic extent of game-farm mallard introgression into wild populations in more westward regions of North America. We used molecular analysis to detect that 35% of samples were pure wild mallard, 12% were early generation hybrids between wild and game-farm mallards (i.e., F1-F3), and the remaining 53% of samples were assigned as part of a hybrid swarm. Percentage of individuals in our study with some form of hybridization with game-farm mallard (65%) was greater than previously detected farther south in the mid-continent (~4%), but less than the Atlantic coast of North America (~ 92%). Stable isotope analysis using δ2Hf suggested that pure wild mallards originated from areas farther north and west than hybrid mallards. More specifically, 17% of all Ohio samples had δ2Hf consistent with more western origins in the prairies, parkland, or boreal regions of the mid-continent of North America, with 55%, 35%, and 10% of these being genetically wild, hybrid swarm, and F3, respectively. We conclude that continued game-farm introgression into wild mallards is not isolated to the eastern population of mallards in North America, and may be increasing and more widespread than previously detected. Mallards in our study had greater incidence of game-farm hybridization than other locales in the mid-continent but less than eastern North American regions suggesting further need to understand game-farm mallard genetic variation and movement across the continent.

A feather hydrogen (δ2H) isoscape for Brazil

Spatial patterns of stable isotopes in animal tissues or “isoscapes” can be used to investigate animal origins in a range of ecological and forensic investigations. Here, we developed a feather hydrogen isotope (δ²H_f) isoscape for Brazil based on 192 samples of feathers from the family Thraupidae from scientific collections. Raw values of δ²H_f ranged from -107.3 to +5.0‰, with higher values at the Caatinga biome (northeast Brazil) and lower values at the Amazon and Pantanal. A Random Forest (RF) method was used to model the spatial surface, using a range of environmental data as auxiliary variables. The RF model indicated a negative relationship between δ²H_f and Mean Annual Precipitation, Precipitation in the Warmest Quarter, and Annual Temperature Range and positive relationships for amount-weighted February-April precipitation δ²H (δ²H_p(Feb-April)) and Mean Annual Solar Radiation. Modelled δ²H_f values ranged from -85.7 to -13.6‰. Ours is the first δ²H_f isoscape for Brazil that can greatly assist our understanding of both ecological and biogeochemical processes controlling spatial variation in δ²H for this region. This isoscape can be used with caution, due to its poor predictive power (as found in other tropical regions) and can benefit from new sample input, new GNIP data, ecological and physiological studies, and keratin standard material better encompassing the range in feather samples from Brazil. So, we encourage new attempts to build more precise feather H isoscapes, as well as isoscapes based on other elements.

Environmental and life-history factors influence inter-colony multidimensional niche metrics of a breeding Arctic marine bird

Human industrialization has resulted in rapid climate change, leading to wide-scale environmental shifts. These shifts can modify food web dynamics by altering the abundance and distribution of primary producers (ice algae and phytoplankton), as well as animals at higher trophic levels. Methylmercury (MeHg) is a neuro-endocrine disrupting compound which biomagnifies in animals as a function of prey choice, and as such bioavailability is affected by altered food web dynamics and adds an important risk-based dimension in studies of foraging ecology. Multidimensional niche dynamics (MDND; δ¹³C, δ¹⁵N, THg; total mercury) were determined among breeding common eider (Somateria mollissima) ducks sampled from 10 breeding colonies distributed across the circumpolar Arctic and subarctic. Results showed high variation in MDND among colonies as indicated by niche size and ranges in δ¹³C, δ¹⁵N and THg values in relation to spatial differences in primary production inferred from sea-ice presence and colony migratory status. Colonies with higher sea-ice cover during the pre-incubation period had higher median colony THg, δ¹⁵N, and δ¹³C. Individuals at migratory colonies had relatively higher THg and δ¹⁵N, and lower δ¹³C, suggesting a higher trophic position and a greater reliance on phytoplankton-based prey. It was concluded that variation in MDND exists among eider colonies which influenced individual blood THg concentrations. Further exploration of spatial ecotoxicology and MDND at each individual site is important to examine the relationships between anthropogenic activities, foraging behaviour, and the related risks of contaminant exposure at even low, sub-lethal concentrations that may contribute to deleterious effects on population stability over time. Overall, multidimensional niche analysis that incorporates multiple isotopic and contaminant metrics could help identify those populations at risk to rapidly altered food web dynamics.

Isotope record tracks changes in historical wintering ranges of a passerine in sub-Saharan Africa

Billions of birds migrate from the Palaearctic to sub-Saharan Africa, yet we are unaware about where exactly they stay over winter and how consistent they have been using these wintering areas over historical times. Here, we inferred the historical wintering areas of Eurasian Golden Orioles (Oriolus oriolus) from stable isotope ratios of feathers. Over the past 200 years, Golden Orioles used two major wintering grounds. Between 1895 and 1971, the relative use of these areas depended on local rainfall intensities. Golden Orioles may depend strongly on humid wintering areas in sub-Saharan Africa, which may put this species at stake when the global climate continues to change.

Migration and non-breeding ecology of the Yellow-breasted Chat Icteria virens

Detailed information spanning the full annual cycle is lacking for most songbird populations. We examined breeding, migration, and non-breeding sites for the Yellow-breasted Chat (Icteria virens, chat). We deployed archival GPS tags and light-level geolocators on breeding chats in British Columbia and light-level geolocators in California from 2013 to 2017 to determine migration routes and non-breeding sites. We examined whether chats overwintered in protected areas and characterized the percent of land cover within 1 km. We used a combination of genetics and stable hydrogen isotopes from feathers collected on non-breeding chats in Nayarit, Mexico (2017-2019) and migrating chats in Chiapas, Mexico (2018) and Veracruz, Mexico (2014-2015) to determine subspecies and infer breeding location. Endangered chats in British Columbia followed the Pacific Flyway and spent the non-breeding period in Sinaloa and Nayarit, Mexico. Two out of five chats spent the non-breeding period in protected areas, and the most common landcover type used was tropical or subtropical broadleaf deciduous forest. We found no mixing of eastern and western chats in our Mexico sites, suggesting strong migratory connectivity at the subspecies level. Western chats likely originating from multiple breeding latitudes spent the non-breeding period in Nayarit. Eastern Yellow-breasted Chats likely breeding across various latitudes migrated through Veracruz and Chiapas. Our results provide precise migration routes and non-breeding locations, and describe habitat cover types for chats, notably an endangered population in British Columbia, which may be valuable for habitat protection and conservation efforts.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s10336-021-01931-8.

Experimental Evaluation of δ2H, δ13C and δ15N Variability in Blood and Feathers of Wild and Captive Birds: Implications for Interspecific Food Web Studies

Stable-hydrogen (δ2H), nitrogen (δ15N), and carbon (δ13C) isotopes are used to decipher broad movement patterns and trophic relationships among diverse species, and an improved understanding of factors controlling natural variation in tissue-isotope measurements will enhance these applications. To evaluate the rearing environment and family-related effects on the isotopic composition of tissues, we cross-fostered nestling tree swallows (Tachycineta bicolor, Vieillot 1808) and American kestrels (Falco sparverius, Linnaeus 1758) by swapping recently hatched birds (<4 days old) among nest boxes and collecting blood and feathers prior to fledging. To assess developmental effects, we measured δ2H in blood and feathers of captive mallard (Anas platyrhynchos, Linnaeus 1758) ducklings challenged energetically during growth. Stable isotope composition was not strongly related to nest box type or natal nest (i.e., family of origin) effects in swallows and kestrels; tissue-isotope composition was related to rearing environment, indicative of differences in nest and parental quality or parental provisioning tactics. Blood and feather δ2H values in swallows were positively related to antecedent maximum ambient temperature, and unrelated to elevated energy expenditure in mallards. The average differences between δ2H in blood and feathers were similar for nestling swallows (27‰, 32‰; two sites) and mallards (26‰, 30‰; two age groups), and lower than in nestling kestrels (50‰). Strong species-specific patterns in blood-feather differences were not observed for δ15N and δ13C in swallows or kestrels; divergent δ2H results may be related to differences in nest ambient conditions, diet composition, or physiological processes affecting hydrogen assimilation during growth and feather synthesis. In swallows, tissue-isotope values reflected parental prey selection from spatially distinct food webs during nestling development with little effect(s) of family of origin, egg composition, or early growth.

Migratory connectivity then and now: a northward shift in breeding origins of a long-distance migratory bird wintering in the tropics

Temporal variation in the connectivity of populations of migratory animals has not been widely documented, despite having important repercussions for population ecology and conservation. Because the long-distance movements of migratory animals link ecologically distinct and geographically distant areas of the world, changes in the abundance and migratory patterns of species may reflect differential drivers of demographic trends acting over various spatial scales. Using stable hydrogen isotope analyses (δ²H) of feathers from historical museum specimens and contemporary samples obtained in the field, we provide evidence for an approximately 600 km northward shift over 45 years in the breeding origin of a species of songbird of major conservation concern (blackpoll warbler, Setophaga striata) wintering in the foothills of the eastern Andes of Colombia. Our finding mirrors predictions of range shifts for boreal-breeding species under warming climate scenarios and habitat loss in the temperate zone, and underscores likely drivers of widespread declines in populations of migratory birds. Our work also highlights the value of natural history collections to document the effects of global change on biodiversity.

Inferring spatial patterns of mercury exposure in migratory boreal songbirds: Combining feather mercury and stable isotope (δ2H) measurements

Migratory songbirds breeding in the Canadian Boreal forest are exposed to mercury (Hg), a potent neurotoxin that impairs avian health, however, the degree of exposure depends on many factors. As breeding grounds are geographically remote and vast, the measurement of Hg in individual birds is impractical particularly at large spatial scales. Here, we present a Canada-wide dataset of nearly 2000 migratory songbirds that were used to assess summer Hg exposure of 15 songbird species sampled during fall migration. We measured Hg concentrations in tail feathers and related those to dietary guild, geographic capture location, age, sex and probable breeding ground locations using feather δ²H. Overall mean (±SE) feather Hg concentration was 1.49 ± 0.03 μg/g (N = 1946): however, a clear geographic gradient in feather Hg concentrations emerged being highest in East and lowest in West. Dietary guild was the next strongest predictor of feather Hg with insectivorous songbirds in Eastern Canada at particular risk due to Hg exposure on summer breeding grounds. This broad-scale assessment of Hg exposure in migratory songbirds in Canada can be used to guide future studies on finer-scale determinants of Hg exposure in birds.

Empirical tests of habitat selection theory reveal that conspecific density and patch quality, but not habitat amount, drive long-distance immigration in a wild bird

Individuals that disperse long distances from their natal site must select breeding patches with no prior knowledge of patch suitability. Despite decades of theoretical studies examining which cues dispersing individuals should use to select breeding patches, few empirical studies have tested the predictions of these theories at spatial scales relevant to long‐distance dispersal in wild animal populations. Here, we use a novel assignment model based on multiple intrinsic markers to quantify natal dispersal distances of Wood Thrush (Hylocichla mustelina) breeding in forest fragments. We show that long‐distance natal dispersal in this species is more frequent than commonly assumed for songbirds and that habitat selection by these individuals is driven by density‐dependence and patch quality but not the amount of habitat surrounding breeding patches. These results represent an important contribution to understanding habitat selection by dispersing individuals, especially with regards to long‐distance dispersal.

Stable isotope (C, N, O, and H) study of a comprehensive set of feathers from two Setophaga citrina

Oxygen, hydrogen, carbon and nitrogen stable isotopes were measured on a comprehensive sampling of feathers from two spring Hooded Warblers (Setophaga citrina) in Texas to evaluate isotopic variability between feathers and during molt. Isotopic homogeneity within each bird was found across all four isotopic systems, supporting the hypothesis that molt in these neotropical migrants is fully completed on the breeding grounds. This homogeneity suggests that the isotopic composition of a single feather is may be representative of the whole songbird. However, each bird was found to have one or two outlier feathers, which could signify regrowth of lost feathers after prebasic molt.

Combining Models of Environment, Behavior, and Physiology to Predict Tissue Hydrogen and Oxygen Isotope Variance Among Individual Terrestrial Animals

Variations in stable hydrogen and oxygen isotope ratios in terrestrial animal tissues are used to reconstruct origin and movement. An underlying assumption of these applications is that tissues grown at the same site share a similar isotopic signal, representative of the location of their origin. However, large variations in tissue isotopic compositions often exist even among conspecific individuals within local populations, which complicates origin and migration inferences. Field-data and correlation analyses have provided hints about the underlying mechanisms of within-site among-individual isotopic variance, but a theory explaining the causes and magnitude of such variance has not been established. Here we develop a mechanistic modeling framework that provides explicit predictions of the magnitude, patterns, and drivers of isotopic variation among individuals living in a common but environmentally heterogeneous habitat. The model toolbox includes isoscape models of environmental isotopic variability, an agent-based model of behavior and movement, and a physiology-biochemistry model of isotopic incorporation into tissues. We compare model predictions against observed variation in hatch-year individuals of the songbird Spotted Towhee (Pipilo maculatus) in Red Butte Canyon, Utah, and evaluate the ability of the model to reproduce this variation under different sets of assumptions. Only models that account for environmental isotopic variability predict a similar magnitude of isotopic variation as observed. Within the modeling framework, behavioral rules and properties govern how animals nesting in different locations acquire resources from different habitats, and birds nesting in or near riparian habitat preferentially access isotopically lighter resources than those associated with the meadow and slope habitats, which results in more negative body water and tissue isotope values. Riparian nesters also have faster body water turnover and acquire more water from drinking (vs. from food), which exerts a secondary influence on their isotope ratios. Thus, the model predicts that local among-individual isotopic variance is linked first to isotopic heterogeneity in the local habitat, and second to how animals sample this habitat during foraging. Model predictions provide insight into the fundamental mechanisms of small-scale isotopic variance and can be used to predict the utility of isotope-based methods for specific groups or environments in ecological and forensic research.

Developing a Stopover-CORT hypothesis: Corticosterone predicts body composition and refueling rate in Gray Catbirds during migratory stopover

Migratory flight is energetically challenging, requiring alternating phases of fuel catabolism and fuel accumulation, accompanied by dramatic changes in body composition and behavior. Baseline corticosterone (CORT; the primary glucocorticoid in birds) is thought to underlie transitions between fuel catabolism during flight, fuel deposition during stopover, and the initiation of migratory flight. However, studies of CORT on stopover physiology and behavior remain disparate efforts, lacking the cohesion of a general hypothesis. Here we develop a Stopover-CORT hypothesis formalizing the relationships among CORT, body condition, and refueling rate in migratory birds. First we expect body mass to increase with triglycerides (TRIG) as birds refuel. Second, based on a synthesis of previous literature, we predict a U-shaped CORT curve over the course of stopover, postulating that elevated CORT at arrival is reactive, responding to poor body condition, while CORT elevation before departure is preparative, driving changes in behavior and body condition. We tested these predictions in Gray Catbirds (Dumetella carolinensis) following a trans-Gulf flight during spring migration. We found baseline CORT was negatively correlated with body condition and TRIG, corresponding with our predictions for arriving and refueling-but not departing-birds. It is possible catbirds undergo regional habitat translocations rather than complete the entire stopover phase at our study site. We propose the Stopover-CORT hypothesis as a useful predictive framework for future studies of the mechanistic basis of stopover physiology. By studying the regulation of stopover refueling and departure, we may better understand physiological limitations to overall migration rate and improve assessments of habitat quality for refueling birds.

Spatial variations in oxygen and hydrogen isotopes in waters and human hair across South Korea

The spatial distribution of isotopic signatures in the form of isoscape is a valuable tool to map their spatial heterogeneity in various environmental settings. However, only limited information about δ¹⁸O and δ²H in water across South Korea is available and to our knowledge no study so far has tried to examine the isotopic heterogeneity of tap water and human scalp hair in South Korea. Here, we present the first national scale analyses of stream water, groundwater, tap water, and human scalp hair isoscapes for South Korea. Stream water, groundwater, tap water, and human scalp hair samples were collected from across South Korea. These samples were analyzed for δ¹⁸O and δ²H, and the isotopic data were then used to generate interpolated δ¹⁸O and δ²H isoscapes for South Korea. The results of linear regression analyses showed strong and significant relationships between δ¹⁸O_hair and δ¹⁸O_water (R² = 0.83, P < 0.002) and between δ²H_hair and δ²H_water (R² = 0.74, P < 0.006), primarily reflecting a close co-relationship between water and hair. The slopes of linear regressions for δ¹⁸O (Δδ¹⁸O_hair/Δδ¹⁸O_water) and δ²H (Δδ²H_hair/Δδ²H_water) suggested that approximately 27% of hydrogen and 36% of oxygen in hair keratin were derived from the local drinking water. Interpolated δ¹⁸O and δ²H isotope maps of stream water, groundwater, and tap water samples collected from across South Korea showed similar spatial patterns of isotope variability. These samples showed a clear latitudinal gradient with high isotopic values in the south which progressively decrease toward the north. The same trends were observed in hair isoscapes as well, and had gradients matching the isotopic pattern of water samples. The strong relationship between water and human hair, and the consistent spatial pattern between them suggest that hair isotope signatures in South Korea can be used in provenance- and forensic-related activities.

Mediterranean precipitation isoscape preserved in bone collagen δ2H

The prehistory of the Mediterranean region has long been a subject of considerable interest, particularly the links between human groups and regions of origin. We utilize the spatial variation in the δ²H and δ¹⁸O values of precipitation (isoscapes) to develop proxies for geographic locations of fauna and humans. Bone collagen hydrogen isotope ratios (δ²H) in cattle (and to a lesser extent, ovicaprids) across the Mediterranean reflect the isotopic differences observed in rainfall (but δ¹⁸O values do not). We conclude that δ²H in herbivore bone collagen can be used as a geolocation tracer and for palaeoenvironmental studies such as tracing past isotopic variations in the global hydrological cycle. In contrast, human bone δ²H values are relatively tightly grouped and highly distinct from precipitation δ²H values, likely due to human-specific food practices and environmental modifications. Given the inter-species variability in δ²H, care should be taken in the species selected for study.

Geographical assignment of polar bears using multi-element isoscapes

Wide-ranging apex predators are among the most challenging of all fauna to conserve and manage. This is especially true of the polar bear (Ursus maritimus), an iconic predator that is hunted in Canada and threatened by global climate change. We used combinations of stable isotopes (¹³C,¹⁵N,²H,¹⁸O) in polar bear hair from > 1000 individuals, sampled from across much of the Canadian Arctic and sub-Arctic, to test the ability of stable isotopic profiles to ‘assign’ bears to (1) predefined managed subpopulations, (2) subpopulations defined by similarities in stable isotope values using quadratic discriminant analysis, and (3) spatially explicit, isotopically distinct clusters derived from interpolated (i.e. ‘kriged’) isotopic landscapes, or ‘isoscapes’, using the partitioning around medoids algorithm. A four-isotope solution provided the highest overall assignment accuracies (~80%) to pre-existing management subpopulations with accuracy rates ranging from ~30–99% (median = 64%). Assignment accuracies of bears to hierarchically clustered ecological groups based on isotopes ranged from ~64–99%. Multivariate assignment to isotopic clusters resulted in highest assignment accuracies of 68% (33–77%), 84% (47–96%) and 74% (53–85%) using two, three and four stable isotope groups, respectively. The resulting spatial structure inherent in the multiple stable isotopic compositions of polar bear tissues is a powerful forensic tool that will, in this case, contribute to the conservation and management of this species. Currently, it is unclear what is driving these robust isotopic patterns and future research is needed to evaluate the processes behind the pattern. Nonetheless, our isotopic approach can be further applied to other apex mammalian predators under threat, such as the large felids, providing that isotopic structure occurs throughout their range.

From stable isotope ecology to forensic isotope ecology - Isotopes' tales

Stable isotope ecology and forensic isotope ecology are not only linked by name. More often than not, knowledge and insights gained through the former serve as a springboard for application focused work of the latter. This review aims to offer a glimpse into the fascinating world of both though with more emphasis on forensic isotope ecology. To this end a selection of past and recent published work is presented and discussed to highlight both potential and limitations of isotopic analytical approaches to the detection of illegal trade in plants and animals.

A multi-isotope (δ13C, δ15N, δ34S, δ2H) approach to establishing migratory connectivity in lesser snow geese: Tracking an overabundant species

Expanding populations of North American midcontinent lesser snow geese (Anser caerulescens caerulescens) have potential to alter ecosystems throughout the Arctic and subarctic where they breed. Efforts to understand origins of harvested lesser snow geese to better inform management decisions have traditionally required mark-recapture approaches, while aerial photographic surveys have typically been used to identify breeding distributions. As a potential alternative, isotopic patterns that are metabolically fixed within newly grown flight feathers following summer molting could provide inferences regarding geographic breeding origin of individuals, without the need for prior capture. Our objective was to assess potential to use four stable isotopes (δ¹³C, δ¹⁵N, δ³⁴S, δ²H) from feather material to determine breeding origins. We obtained newly grown flight feathers from individuals during summer banding at three Arctic and two subarctic breeding colonies in 2014 (n = 56) and 2016 (n = 45). We used linear discriminant analyses to predict breeding origins from models using combinations of stable isotopes as predictors and evaluated model accuracy when predicting colony, subregion, or subpopulation levels. We found a strong inverse relationship between δ²H values and increasing latitude (R² = 0.83), resulting in differences (F_{4, 51} = 90.41, P < 0.0001) among sampled colonies. No differences in δ¹³C or δ¹⁵N were detected among colonies, although δ³⁴S in Akimiski Island, Baffin Island, and Karrak Lake were more enriched (F_{4, 51} = 11.25, P < 0.0001). Using δ²H values as a predictor, discriminant analyses improved accuracy in classification level as precision decreased [model accuracy = 67% (colony), 88% (subregion), 94% (subpopulation)]. Application of the isotopic methods we describe could be used to provide an alternative monitoring method of population metrics, such as overall breeding population distribution, region-specific productivity and migratory connectivity that are informative to management decision makers and provide insight into cross-seasonal effects that may influence migratory behavior.

Seasonal Migration of Pantala flavescens Across the Bohai Strait in Northern China

Pantala flavescens (Fabricius 1798) (Odonata: Libellulidae) is one of the most common species of migratory dragonflies. P. flavescens adults were captured by a searchlight trap on Beihuang Island (BH Island; 38°24'N, 120°55'E) from 2003 to 2016, where there is no freshwater. This inspired our research to analyze the pattern of seasonal migration and population dynamics. Stable hydrogen isotope measurement and the Hybrid Single Particle Lagrangian Integrated Trajectory (HYSPLIT) were used to simulate the migration pathway of P. flavescens between different breeding habitats. The results showed that there was no significant difference among population numbers of this overseas migration across years (F13, 2161 = 0.85, P = 0.604); however, the numbers were significantly different across months (F5, 2161 = 3.91, P = 0.003). Our geospatial natal assignment model suggested that P. flavescens trapped on BH were originated in different geographical regions and might have three movement strategies: wandering around northern China and north-bound (positive) and south-bound (negative) movements. Among them, the majority were engaged in wandering around northern China. Model simulations suggested that P. flavescens toured around BH. The results contribute to the knowledge of P. flavescens population ecology in a large-scale geographic region and will aid in the prediction and interpretation of insect migration patterns in response to climate change.

Stable hydrogen isotope variability within and among plumage tracts (δ2HF) of a migratory wood warbler

Hydrogen isotope analysis of feather keratin (δ²H_F) has become an essential tool for tracking movements between breeding and wintering populations of migratory birds. In particular, δ²H_F has been used to create δ²H_F isoscapes that can be used to assign the geographic origins of molt. The majority of past studies have sampled a portion of a single feather as an isotopic proxy for the entire plumage although surprisingly little is known about variation of stable isotopes within and between feather tracts of individuals in local populations. Here we examine δ²H_F variation in 24 pterylographic variables (9 primaries, 6 secondaries, 6 rectrices, and 3 patches of ventral contour feathers) in individual specimens of black-throated blue warbler (Setophaga caerulescens) breeding in the Big Santeetlah Creek watershed (5350 ha), southern Appalachian Mountains. By restricting our study to territorial ASY males (after second year) inhabiting a small watershed, we could focus on δ²H_F variation generated during the complete prebasic annual molt in a circumscribed population while factoring out age and sexual differences in foraging behavior, isotopic incorporation, and post-breeding dispersal. Summed within-individual variation (δ²H_F) across 24 pterylographic variables ranged from 12 to 60‰ (= 21.8 ± 9.4‰), with 81% of the individuals exhibiting variation ≥ 16‰ (reproducibility of analyses was ≤ 4 ‰). Adjacent feathers in feather tracts tend to have more similar δ²H_F values than feathers grown weeks apart, consistent with the stepwise replacement of flight feathers. The pooled population sample exhibited significant δ²H_F variation in primaries (-78 to -21‰), secondaries (-80 to -17‰), rectrices (-78 to -23‰), and ventral contour feathers (-92 to -32‰). Strong year effects in δ²H_F variation were observed in each of the 24 pterylographic variables. Altitudinal effects were observed only in ventral contour feathers. The current findings demonstrate that within-individual variation (δ²H_F) may be much greater than previously thought in migratory species that molt on or near breeding territories. Our study also highlights the need for greater pterylographic precision in research design of isotope-based studies of avian movement. Within-individual and within-population δ²H_F variation should be incorporated in geographic assignment models. In a broader context, the staggered Staffelmauser pattern of molt in wood warblers provides an exceptional view of the seasonal variation of hydrogen isotopes circulating in blood plasma during the six-week period of annual molt.

Isotope (δ13C, δ15N, δ2H) diet-tissue discrimination in African grey parrot Psittacus erithacus: implications for forensic studies

Diet-tissue isotopic relationships established under controlled conditions are informative for determining the dietary sources and geographic provenance of organisms. We analysed δ¹³C, δ¹⁵N, and non-exchangeable δ²H values of captive African grey parrot Psittacus erithacus feathers grown on a fixed mixed-diet and borehole water. Diet-feather Δ¹³C and Δ¹⁵N discrimination values were +3.8 ± 0.3 ‰ and +6.3 ± 0.7 ‰ respectively; significantly greater than expected. Non-exchangeable δ²H feather values (-62.4 ± 6.4 ‰) were more negative than water (-26.1 ± 2.5 ‰) offered during feather growth. There was no positive relationship between the δ¹³C and δ¹⁵N values of the samples along each feather with the associated samples of food offered, or the feather non-exchangeable hydrogen isotope values with δ²H values of water, emphasising the complex processes involved in carbohydrate, protein, and income water routing to feather growth. Understanding the isotopic relationship between diet and feathers may provide greater clarity in the use of stable isotopes in feathers as a tool in determining origins of captive and wild-caught African grey parrots, a species that is widespread in aviculture and faces significant threats to wild populations. We suggest that these isotopic results, determined even in controlled laboratory conditions, be used with caution.

Stable isotopes can be used to infer the overwintering locations of prebreeding marine birds in the Canadian Arctic

Although assessments of winter carryover effects on fitness‐related breeding parameters are vital for determining the links between environmental variation and fitness, direct methods of determining overwintering distributions (e.g., electronic tracking) can be expensive, limiting the number of individuals studied. Alternatively, stable isotope analysis in specific tissues can be used as an indirect means of determining individual overwintering areas of residency. Although increasingly used to infer the overwintering distributions of terrestrial birds, stable isotopes have been used less often to infer overwintering areas of marine birds. Using Arctic‐breeding common eiders, we test the effectiveness of an integrated stable isotope approach (13‐carbon, 15‐nitrogen, and 2‐hydrogen) to infer overwintering locations. Knowing the overwinter destinations of eiders from tracking studies at our study colony at East Bay Island, Nunavut, we sampled claw and blood tissues at two known overwintering locations, Nuuk, Greenland, and Newfoundland, Canada. These two locations yielded distinct tissue‐specific isotopic profiles. We then compared the isotope profiles of tissues collected from eiders upon their arrival at our breeding colony, and used a k‐means cluster analysis approach to match arriving eiders to an overwintering group. Samples from the claws of eiders were most effective for determining overwinter origin, due to this tissue's slow growth rate relative to the 40‐day turnover rate of blood. Despite taking an integrative approach using multiple isotopes, k‐means cluster analysis was most effective when using 13‐carbon alone to assign eiders to an overwintering group. Our research demonstrates that it is possible to use stable isotope analysis to assign an overwintering location to a marine bird. There are few examples of the effective use of this technique on a marine bird at this scale; we provide a framework for applying this technique to detect changes in the migration phenology of birds' responses to rapid changes in the Arctic.

Linking the wintering and breeding grounds of warblers along the Pacific Flyway

Long‐distance migration is a behavior that is exhibited by many animal groups. The evolution of novel migration routes can play an important role in range expansions, ecological interactions, and speciation. New migration routes may evolve in response to selection in favor of reducing distance between breeding and wintering areas, or avoiding navigational barriers. Many migratory changes are likely to evolve gradually and are therefore difficult to study. Here, we attempt to connect breeding and wintering populations of myrtle warblers (Setophaga coronata coronata) to better understand the possible evolution of distinct migration routes within this species. Myrtle warblers, unlike most other warblers with breeding ranges primarily in eastern North America, have two disjunct overwintering concentrations—one in the southeastern USA and one along the Pacific Coast—and presumably distinct routes to‐and‐from these locations. We studied both myrtle and Audubon's warblers (S. c. auduboni) captured during their spring migration along the Pacific Coast, south of the narrow region where these two taxa hybridize. Using stable hydrogen isotopes and biometric data, we show that those myrtle warblers wintering along the southern Pacific Coast of North America are likely to breed at high latitudes in Alaska and the Yukon rather than in Alberta or further east. Our interpretation is that the evolution of this wintering range and migration route along the Pacific Coast may have facilitated the breeding expansion of myrtle warblers into northwestern North America. Moreover, these data suggest that there may be a migratory divide within genetically similar populations of myrtle warblers.

Connecting the dots: Stopover strategies of an intercontinental migratory songbird in the context of the annual cycle

The phases of the annual cycle for migratory species are inextricably linked. Yet, less than five percent of ecological studies examine seasonal interactions. In this study, we utilized stable hydrogen isotopes to geographically link individual black‐and‐white warblers (Mniotilta varia) captured during spring migration with breeding destinations to understand a migrant's stopover strategy in the context of other phases of the annual cycle. We found that stopover strategy is not only a function of a bird's current energetic state, but also the distance remaining to breeding destination and a bird's time‐schedule, which has previously been linked to habitat conditions experienced in the preceding phase of the annual cycle. Birds in close proximity to their breeding destination accumulate additional energy reserves prior to arrival on the breeding grounds, as reflected by higher migratory condition upon arrival, higher refueling rates measured via blood plasma metabolites, and longer stopover durations compared to birds migrating to breeding destinations farther from the stopover site. However, late birds near their breeding destination were more likely to depart on the day of arrival (i.e., transients), and among birds that stopped over at the site, the average duration of stopover was almost half the time of early conspecifics, suggesting late birds are trying to catch‐up with the overall time‐schedule of migration for optimal arrival time on the breeding grounds. In contrast, birds with long distances remaining to breeding destinations were more likely to depart on the day of arrival and primarily used stopover to rest before quickly resuming migration, adopting similar strategies regardless of a bird's time‐schedule. Our study demonstrates that migrants adjust their en route strategies in relation to their time‐schedule and distance remaining to their breeding destination, highlighting that strategies of migration should be examined in the context of other phases of the annual cycle.

Refining Stable Oxygen and Hydrogen Isoscapes for the Identification of Human Remains in Mississippi

Isoscape refinement is an essential component for accurately predicting region-of-origin in forensic investigations involving isotope analysis of unidentified human remains. Stable oxygen (δ¹⁸ O) and hydrogen (δ² H) isotopes were measured from 57 tap water samples collected across Mississippi to model refined isoscapes for the state. A tap water conversion equation, δ¹⁸ O_tw =1.64 δ¹⁸ Op-31.35, was developed for the southeastern USA to test the prediction accuracy of the δ¹⁸ O_tw isoscape using individuals with known residential histories. A local Mississippi resident (USAFA-134) was assigned with 90% probability to the correct region-of-origin reported by the participant. Assignments for Georgia residents (USAFA-118 and USAFA-205) had variable results, predicting USAFA-118 from Mississippi and USAFA-205 as a nonlocal resident. Stable isotope values often overlap geographically and a multi-isotope approach should be used when narrowing region(s)-of-origin(s). This study demonstrates the utility of refining isoscapes and the importance of tissue calibration in prediction assignments of human remains.

Incorporating breeding abundance into spatial assignments on continuous surfaces

Determining the geographic connections between breeding and nonbreeding populations, termed migratory connectivity, is critical to advancing our understanding of the ecology and conservation of migratory species. Assignment models based on stable isotopes historically have been an important tool for studying migratory connectivity of small-bodied species, but the low resolution of these assignments has generated interest into combining isotopes with other sources in information. Abundance is one of the most appealing data sources to include in isotope-based assignments, but there are currently no statistical methods or guidelines for optimizing the contribution of stable isotopes and abundance for inferring migratory connectivity. Using known-origin stable-hydrogen isotope samples of six Neotropical migratory bird species, we rigorously assessed the performance of assignment models that differentially weight the contribution of the isotope and abundance data. For two species with adequate sample sizes, we used Pareto optimality to determine the set of models that simultaneously minimized both assignment error rate and assignment area. We then assessed the ability of the top models from these two species to improve assignments of the remaining four species compared to assignments based on isotopes alone. We show that the increased precision of models that include abundance is often offset by a large increase in assignment error. However, models that optimally weigh the abundance data relative to the isotope data can result in higher precision and, in some cases, lower error than models based on isotopes alone. The top models, however, depended on the distribution of relative breeding abundance, with patchier distributions requiring stronger downweighting of abundance, and we present general guidelines for future studies. These results confirm that breeding abundance can be an important source of information for studies investigating broad-scale movements of migratory birds and potentially other taxa.

Migration distance as a selective episode for wing morphology in a migratory insect

BACKGROUND

Selective pressures that occur during long-distance migration can influence morphological traits across a range of taxa. In flying insects, selection should favour individuals that have wing morphologies that increase energy efficiency and survival. In monarch butterflies, differences in wing morphology between migratory and resident populations suggest that migratory populations have undergone selection for larger (as measured by length and area) and more elongated (as measured by roundness and aspect ratio) forewings. However, selection on wing morphology may also occur within migratory populations, particularly if individuals or populations consistently migrate different distances.

RESULTS

Using 613 monarch butterflies that were collected on the Mexican wintering grounds between 1976 - 2014, we tested whether monarch wing traits were associated with migratory distance from their natal areas in eastern North America (migration range: 774-4430 km), as inferred by stable-hydrogen (δ²H) and -carbon (δ¹³C) isotopic measurements. Monarchs that migrated farther distances to reach their overwintering sites tended to have longer and larger wings, suggesting positive selective pressure during migration on wing length and area. There was no relationship between migration distances and either roundness or aspect ratio.

CONCLUSIONS

Our results provide correlative evidence that the migratory period may act as a selective episode on monarch butterfly wing morphology, although selection during other portions of the annual cycle, as well as extensive mixing of individuals from various natal locations on the breeding grounds, likely counteracts directional selection of migration on morphology.

Linking otolith microchemistry and dendritic isoscapes to map heterogeneous production of fish across river basins

Production patterns of highly mobile species, such as anadromous fish, often exhibit high spatial and temporal heterogeneity across landscapes. Such variability is often asynchronous in time among habitats, which stabilizes production at aggregate scales of complexity. Reconstructing production patterns explicitly in space and time across multiple scales, however, remains difficult but is important for prioritizing habitat conservation. This is especially true for fishes inhabiting river basins due to long-range dispersal, high mortality at early life stages, complex population structure and elusive life history variation. We develop a new approach for mapping production patterns of Pacific salmon across a large river basin by integrating otolith microchemistry and dendritic isoscape models. The geographically continuous Bayesian assignment framework presented here yielded high accuracies (>90%) and relatively high precisions (precisions <4%; i.e., assignment areas of <530 river km of the 13 100 km total river length) when used to determine the natal source of known-origin juvenile Chinook salmon captured throughout the study region. Integrating these methods enabled us to base estimates of provenance and habitat use of individuals on a per location basis using strontium isotopic data throughout the continuous spatial domain of a river network. Such a framework provides substantial advantages over the more common nominal approach to employing otolith microchemistry to reconstruct movement patterns of fish. In doing so, we reconstructed the spatial production patterns of adult Chinook salmon returning to a large watershed in Bristol Bay, Alaska and illustrate the power of such an approach to conservation efforts.

Multi-Isotopic (δ2H, δ13C, δ15N) Tracing of Molt Origin for Red-Winged Blackbirds Associated with Agro-Ecosystems

We analyzed stable-hydrogen (δ2H), carbon (δ13C) and nitrogen (δ 15N) isotope ratios in feathers to better understand the molt origin and food habits of Red-winged Blackbirds (Agelaius phoeniceus) near sunflower production in the Upper Midwest and rice production in the Mid-South of the United States. Outer primary feathers were used from 661 after-second-year (ASY) male blackbirds collected in Minnesota, Montana, North Dakota and South Dakota (spring collection), and Arkansas, Louisiana, Mississippi, Missouri and Texas (winter collection). The best-fit model indicated that the combination of feather δ2H, δ13C and δ15N best predicted the state of sample collections and thus supported the use of this approach for tracing molt origins in Red-winged Blackbirds. When considering only birds collected in spring, 56% of birds were classified to their collection state on the basis of δ2H and δ13C alone. We then developed feather isoscapes for δ13C based upon these data and for δ2H based upon continental patterns of δ2H in precipitation. We used 81 birds collected at the ten independent sites for model validation. The spatially-explicit assignment of these 81 birds to the δ2H isoscape resulted in relatively high rates (~77%) of accurate assignment to collection states. We also modeled the spatial extent of C3 (e.g. rice, sunflower) and C4 (corn, millet, sorghum) agricultural crops grown throughout the Upper Midwest and Mid-South United States to predict the relative use of C3- versus C4-based foodwebs among sampled blackbirds. Estimates of C3 inputs to diet ranged from 50% in Arkansas to 27% in Minnesota. As a novel contribution to blackbird conservation and management, we demonstrate how such feather isoscapes can be used to predict the molt origin and interstate movements of migratory blackbirds for subsequent investigations of breeding biology (e.g. sex-specific philopatry), agricultural depredation, feeding ecology, physiology of migration and sensitivity to environmental change.

Within-Site Variation in Feather Stable Hydrogen Isotope (δ2Hf) Values of Boreal Songbirds: Implications for Assignment to Molt Origin

Understanding bird migration and dispersal is important to inform full life-cycle conservation planning. Stable hydrogen isotope ratios from feathers (δ2Hf) can be linked to amount-weighted long-term, growing season precipitation δ2H (δ2Hp) surfaces to create δ2Hf isoscapes for assignment to molt origin. However, transfer functions linking δ2Hp with δ2Hf are influenced by physiological and environmental processes. A better understanding of the causes and consequences of variation in δ2Hf values among individuals and species will improve the predictive ability of geographic assignment tests. We tested for effects of species, land cover, forage substrate, nest substrate, diet composition, body mass, sex, and phylogenetic relatedness on δ2Hf from individuals at least two years old of 21 songbird species captured during the same breeding season at a site in northeastern Alberta, Canada. For four species, we also tested for a year × species interaction effect on δ2Hf. A model including species as single predictor received the most support (AIC weight = 0.74) in explaining variation in δ2Hf. A species-specific variance parameter was part of all best-ranked models, suggesting variation in δ2Hf was not consistent among species. The second best-ranked model included a forage substrate × diet interaction term (AIC weight = 0.16). There was a significant year × species interaction effect on δ2Hf suggesting that interspecific differences in δ2Hf can differ among years. Our results suggest that within- and among-year interspecific variation in δ2Hf is the most important source of variance typically not being explicitly quantified in geographic assignment tests using non-specific transfer functions to convert δ2Hp into δ2Hf. However, this source of variation is consistent with the range of variation from the transfer functions most commonly being propagated in assignment tests of geographic origins for passerines breeding in North America.