Abstract

Advances in the application of stable isotopes have allowed the quantitative evaluation of previously cryptic ecological processes. In particular, researchers have utilized the predictable spatial patterning in natural abundance of isotopes to better understand animal dispersal and migration. However, quantifying dispersal via natural abundance alone has proven to be of limited utility for species exhibiting short- or mid-distance dispersal events, including most mammalian species. In previous experimental work, we demonstrated that consumption of 1 dose of isotopically enriched baits elicited a distinct “mark” in hair of captive martens (Martes spp.). Herein, we report findings from a field test of our isotopic enrichment approach to mark free-ranging animals and quantify dispersal of martens across forest stands at sites in southeastern Alaska and northern British Columbia. In the field, we supplemented bait used in single-capture hair traps with the amino acid glycine artificially enriched in 2H, 13C, or 15N. By applying unique combinations of artificially enriched isotopic markers within discrete forest stands, the isotopic signature of collected hair reflected the forest patch where the individual originated. From our isotopic marks, we were able to infer dispersal events between forested stands and, thus, estimate rates and approximate dispersal distances. Our findings demonstrate that isotopic enrichment can be a cost-effective method to mark the hair of midsized mammals for the quantification of dispersal.

The use of stable isotopes, particularly those of nitrogen, carbon, and hydrogen, has enabled biologists to quantify a myriad of cryptic ecological and evolutionary processes: trophic interactions and dietary overlap, organismal physiol- ogy and nutrient allocation, and animal movement and behavior (Ben-David and Flaherty 2012; Newsome et al. 2012; Wunder 2012). In particular, the natural abundance of stable isotopes has provided important insights into animal migration and dispersal (Rubenstein and Hobson 2004; Wunder 2012). For example, researchers have previously employed the δ15N and δ13C signatures of animal tissue to reconstruct the timing of individual movement between foraging patches that feature different food resources (e.g., Burton and Koch 1999; Schell et al. 1988). Reliable assignment of individuals to such patches based on the natural abundance of 15N and 13C requires certain conditions that may not be attainable for many field studies; notably, habitats need to contain diet items that are isotopically distinct to yield detectable changes in animal tissue and use of the different habitats needs to be temporally discrete (Flaherty and Ben-David 2010). Consequently, naturally occurring hydrogen isotopes (2H), or deuterium (D), has gained attention as an effective alternative (Rubenstein and Hobson 2004; Wunder 2012). Because 8D values of precipitation (Bowen et al. 2005; Hobson 1999) exhibit predictable latitudinal gradients, the 2H signature of animal tissue can be assigned to a geographical location, allowing the identification of wintering grounds or” the quantification of long-distance migration events. Measuring the δD values of animal tissue and relating it back to global isoclines of δD values has been successfully used for volant animals such as butterflies (Dockx et al. 2004), passerines (Hobson and Wassenaar 1997), and microchirop- terans (Cryan et al. 2004).

Quantifying dispersal with deuterium is rare in mammalogical studies, particularly those featuring nonvolant and terrestrial species, which generally do not disperse sufficient distances to acquire detectable changes in natural abundance of δD (e.g., <100 km—Wunder and Norris 2008). Physiol- ogists (Boutton 1991; Hirons et al. 2001) and ecosystem ecologists (Hall and Tank 2003) have, however, developed procedures to artificially enrich organic and inorganic compounds with stable isotopes (such as 13C and 15N) and track the flux of resources and nutrients. More recently, isotopic enrichment has been applied to mark and subsequent- ly recapture individuals to quantify dispersal of aquatic invertebrates (Macneale et al. 2004; Wanner et al. 2006) and seeds of plants (Carlo et al. 2009). The application of enriched isotopes to mark individuals and infer dispersal distances has not been extended to the study of wild vertebrates.

Through captive feeding trials, we previously demonstrated that the hair of midsized mammals can be cost-effectively marked via artificially enriched diets (Pauli et al. 2009). For wild-caught Pacific martens (Martes caurina—MacDonald and Cook 2009) we found that during active hair growth, consumption of 1 bait containing the amino acid glycine artificially enriched with 2H, 13C, and 15N elicited a distinct “mark” in their hair. However, this approach had not been tested in the field or applied to track individual movement patterns and document long-distance dispersal events. Herein, we report findings from a test of this novel, noninvasive method to quantify dispersal rates of American (Martes americana) and Pacific (M. caurina) martens in the forested island archipelagos of southeastern Alaska and British Columbia, Canada. These 2 obligate forest carnivores diverged during the late Quaternary period in separate refugia exhibiting different levels of glacial fragmentation (Graham and Graham 1994; Stone and Cook 2002). Because the Pacific marten evolved under naturally fragmented refugia in the southwestern United States (in and around present-day eastern Colorado, Utah, and Nevada), its populations likely existed as metapopulations (den Boer 1970; Small et al. 2003; Stone et al. 2002). In contrast, American martens arose within a southeastern refugium, ranging as far south as present-day Alabama, of relatively contiguous forest (Small et al. 2003; Stone et al. 2002). Given this evolutionary context, we hypothesized that the Pacific marten would possess greater dispersal power and exhibit higher dispersal rates and distances than the American marten. By applying unique combinations of artificially enriched isotopic markers to different forest stands, and subsequently collecting hair samples via noninvasive hair snares and carcasses killed during the commercial trapping season, we evaluated the potential to mark and track the movement and dispersal of free-ranging individuals and explored differences in dispersal rates and distances between the 2 species.

Materials and Methods

Our study sites were located within temperate rain forests of southeastern Alaska, and northern British Columbia (Fig. 1). Specifically, fieldwork was conducted on Prince of Wales Island (POW [55° 38′N, 132° 48′W]; 2006), Admiralty Island (ADM [57° 43′N, 134° 22′W]; 2007), and on Misty Fjords National Monument (MFNM [55° 38′N, 130° 37′W]) in southeastern Alaska, and Haida Gwaii (Queen Charlotte Islands; QCI [53° 38′N, 132° 23′W]; 2008) in British Colum- bia. From August to September, we deployed single-capture hair traps across 4 sampling areas within each island. Sampling areas were selected from geographic information system layers based on habitat characteristics, size, and distance from one another. Within each sampling area, 3 or 4 traplines were set. To test the effectiveness of isotope labeling for quantifying dispersal rates and distances of martens between forest stands, we supplemented bait used in hair traps with the amino acid glycine enriched with 99% 2H, 13C, or 15N (ISOTEC, Miamisburg, Ohio). We applied a unique label to each forest stand so the isotopic signature of the hair reflected the location where the individual originated (Fig. 1). We estimated the quantity of glycine needed to enrich the guard hair of martens for each isotope following the equations provided in Pauli et al. (2009).

Fig. 1

Location of study sites within southeastern Alaska (Admiralty Island [ADM], Misty Fjords National Monument [MFNM], and Prince of Wales Island [POW]) and northern British Columbia, Canada (Queen Charlotte Islands [QCI]) where we used artificially enriched isotopic baits to mark American (Martes americana) and Pacific (M. caurind) martens for the quantification of dispersal, 2006–2008. Symbols denote the type of enrichment applied within each area: ▲+ = 15N, ■ = 2H, ● = 13C, ♦ = a combination of 15N and 2H.

Fig. 1

Location of study sites within southeastern Alaska (Admiralty Island [ADM], Misty Fjords National Monument [MFNM], and Prince of Wales Island [POW]) and northern British Columbia, Canada (Queen Charlotte Islands [QCI]) where we used artificially enriched isotopic baits to mark American (Martes americana) and Pacific (M. caurind) martens for the quantification of dispersal, 2006–2008. Symbols denote the type of enrichment applied within each area: ▲+ = 15N, ■ = 2H, ● = 13C, ♦ = a combination of 15N and 2H.

We simultaneously collected marten fur and isotopically labeled martens with a single-capture hair trap that allowed visiting animals to escape, but prevented subsequent access by additional animals (Pauli et al. 2008). Hairs of the departing animal were collected by wire brushes. We baited each trap with deer scraps and jam (the latter containing the unique isotopic label), and applied fish oil and essence of skunk as olfactory lures. Traps were checked every 5th day for 5 weeks. Hair samples were collected from the wire brushes using forceps, placed in 1.5-ml microcentrifuge vials containing silica bead desiccant, identified to species based on dorsal guard-hair characteristics, and stored frozen. We also opportunistically collected fresh marten scat samples found near hair traps on ADM for analyses of 13C and 15N. Finally, we purchased carcasses from commercial trappers that caught martens within or in close proximity to our study sites during the winter trapping season (December-January) that followed our autumn labeling effort. Sampling procedures adhered to the guidelines for use of mammals in research set forth by the American Society of Mammalogists (Sikes et al. 2011) and methods were approved by the Institutional Animal Care and Use Committee at the University of Wyoming. Sampling permits were obtained through the Alaska Department of Fish and Game.

Hair samples from both carcasses and traps were rinsed 3 times with 2:1 chloroform : methanol solution to remove surface oils (Cryan et al. 2004), dried for 72 h at 60° C, and homogenized with surgical scissors. Fecal samples were dried and ground to a fine powder. All samples were weighed, placed in tin (13C and 15N) or silver (2H) capsules, and submitted to the Stable Isotope Facility at the University of Wyoming. Analysis of δ13C and δ15N was conducted with a Costech 4010 elemental analyzer (Costech Analytical, Valencia, California) and δD with a Finnigan TC/EA (High Temperature Conversion Elemental Analyzer; Thermo Finnigan, Bremen, Germany) attached to a Thermo Finnigan DeltaPLUS XP Continuous Flow Isotope Ratio Mass Spectrometer (Thermo Finnigan). Results are provided as per mil (parts per thousand [‰]) ratios relative to the international standards of Peedee Belemnite (PDB; δ13C), atmospheric nitrogen (AIR; δ15N), and Standard Mean Ocean Water (SMOW; δD) with calibrated internal laboratory standards. Because all samples were processed and analyzed in the same geographic location we did not correct the raw δD data for atmospheric exchange (Wassenaar and Hobson 2000).

Natural abundance values of δ13C and δ15N are highly correlated for martens of southeastern Alaska because these terrestrial carnivores also consume salmon (Oncorhynchus spp.) and other marine organisms that are naturally enriched in both 13C and 15N (Ben-David et al. 1997). Therefore, we constructed a baseline regression model from 69 marten hair samples collected from 3 locations in southeastern Alaska (POW, ADM, and MFNM) that were either outside of our sampling areas or sampling period, and therefore could not have come from animals that ingested an isotopically labeled bait. Hair samples collected within our study area with isotopic values exceeding the 99% prediction interval of the regression equation were classified as originating from individuals that ingested a 15N-enriched bait; those with values less than the lower 99% prediction interval were classified as having been marked with 13C. We explored whether natural abundance of 13C, 15N, and deuterium differed between locations with 1-way analyses of variance. Because natural abundance values of δD are influenced by latitude and distance from the coastline (Bowen et al. 2005; Wunder 2012), we pooled samples collected along the coast of MFNM with those from coastal POW (because they were collected from the same latitude). We did not attempt to identify animals marked with deuterium from QCI because we lacked isotopically unmarked (i.e., reference) individuals. We identified marked individuals as those for which the signature of hair samples exceeded the upper 99% prediction limit of means for each location.

We identified fecal samples as isotopically marked using dual criteria of elemental abundance of nitrogen, and isotopic enrichment of 15N and 13C. First, samples were considered marked when they possessed <5.5% nitrogen, which is indicative of feces free of animal tissue (Ben-David et al. 1998), but were enriched in 15N and 13C above available vegetation in the region (i.e., >5‰ for 15N or >—25‰ for 13C—Ben-David et al. 1997). Second, scat samples with >5.5% nitrogen, indicating that animal tissues were present in the excreta, were deemed isotopically marked when 15N and 13C exceeded the signature of the most enriched animal diet items (i.e., 15N > 15‰ or 13C > -17‰).

In southeastern Alaska, the average home-range size of male martens is 7.8 km2 and females have a mean home-range size of 5.3 km2 (Flynn and Schumacher 1999). Thus, based on isotopic labeling and location of sample collection we identified potential dispersers as those animals that moved a distance of >2 home-range diameters, assuming a circular home range, away from their original location of capture (i.e., males moving >6.3 km; females moving >5.2 km). We estimated minimum dispersal distances by measuring the distance between nearest trap location supplemented with the enriched isotope to the final trapped location in ArcMap (ESRI, Redlands, California). We compared dispersal rates (the frequency of observed dispersal events relative to the number of marked animals) between American (for combined samples from POW and MFNM) and Pacific (ADM and QCI) martens with a G-test. Similarly, we compared dispersal distances for the 2 species with a 2-sample t-test (Zar 1999).

Results

In total, we collected 372 hair samples with our hair traps and obtained 366 hair samples from trapper-killed marten carcasses (Table 1). As predicted, δ13C and δ15N were strongly correlated (r = 0.930) for unmarked martens in our external samples collected across southeastern Alaska (Fig. 2A). We detected significant differences among sites in the δ13C signature of unmarked martens (F2,68 = 4.26, P = 0.018; Fig. 3A), and marginally significant differences in δ15N (F2,68 = 2.76, P = 0.071; Fig. 3B). Site-based differences were driven by samples from ADM being slightly enriched in both 15N and 13C. Nonetheless, omitting ADM samples from the regression model resulted in overlapping parameter estimates (β0 = 34.5 ± 1.4 SE and β1 = 1.2 =± 0.06) compared to the model that included ADM samples (Fig. 2A). Using the 99% prediction interval from the regression model that included ADM unmarked samples, we identified 21 carcasses marked by 13C and 9 marked by 15N (Fig. 2B). For hair obtained from the 372 samples collected in our noninvasive traps, we identified 37 instances in which the captured animal was marked by 15N and 26 instances in which the marten was marked by 13C (Fig. 2C). Of the 112 marten scat samples we collected along traplines on ADM, we deemed 35% (n = 41) as enriched in 15N (Fig. 4A) and 5% (n = 6) enriched in 13C (Fig. 4B).

Fig. 2

A) Linear regression model (β0 = 33.7 ±1.1 and β1 = 1.1 ± 0.05) and 99% prediction interval of δ15N to δ13C for 69 unmarked martens (Martes americana and M. caurina) from southeastern Alaska. δ15N and δ13C signatures of hair samples from martens in southeastern Alaska (Prince of Wales Island [POW], Misty Fjords National Monument [MFNM], and Admiralty Island [ADM]) and northern British Columbia, Canada (Queen Charlotte Islands [QCI]) collected from B) commercially trapped marten carcasses and C) noninvasively sampled free-ranging martens. We identified samples above the 99% prediction interval as marked with baits enriched in 15N; those below were considered artificially enriched by 13C. Isotopically marked samples are shown in black; unmarked samples are unfilled.

Fig. 2

A) Linear regression model (β0 = 33.7 ±1.1 and β1 = 1.1 ± 0.05) and 99% prediction interval of δ15N to δ13C for 69 unmarked martens (Martes americana and M. caurina) from southeastern Alaska. δ15N and δ13C signatures of hair samples from martens in southeastern Alaska (Prince of Wales Island [POW], Misty Fjords National Monument [MFNM], and Admiralty Island [ADM]) and northern British Columbia, Canada (Queen Charlotte Islands [QCI]) collected from B) commercially trapped marten carcasses and C) noninvasively sampled free-ranging martens. We identified samples above the 99% prediction interval as marked with baits enriched in 15N; those below were considered artificially enriched by 13C. Isotopically marked samples are shown in black; unmarked samples are unfilled.

Fig. 3

Boxplot depicting mean, standard error, and 99% prediction interval for A) 15N (‰) and B) 13C values obtained from 69 unmarked martens (Martes americana and M. caurina) from 3 sites in southeastern Alaska (Misty Fjords National Monument [MFNM], Prince of Wales Island [POW], and Admiralty Island [ADM]).

Fig. 3

Boxplot depicting mean, standard error, and 99% prediction interval for A) 15N (‰) and B) 13C values obtained from 69 unmarked martens (Martes americana and M. caurina) from 3 sites in southeastern Alaska (Misty Fjords National Monument [MFNM], Prince of Wales Island [POW], and Admiralty Island [ADM]).

Fig. 4

Isotopic signatures and percent elemental nitrogen of Pacific marten (Martes caurina) feces opportunistically collected along trapping lines in Admiralty Island, Alaska. We identified samples as isotopically marked when elemental abundance of nitrogen was <5.5% (free of animal tissue), but had >5‰ or >-25‰ signatures in 15N or 13C, respectively. Samples with >5.5% nitrogen (presence of animal tissue) were deemed isotopically marked when 15N was >15‰ or 13C was > —17‰. Thresholds are denoted with dashed lines; enriched samples are shown in black.

Fig. 4

Isotopic signatures and percent elemental nitrogen of Pacific marten (Martes caurina) feces opportunistically collected along trapping lines in Admiralty Island, Alaska. We identified samples as isotopically marked when elemental abundance of nitrogen was <5.5% (free of animal tissue), but had >5‰ or >-25‰ signatures in 15N or 13C, respectively. Samples with >5.5% nitrogen (presence of animal tissue) were deemed isotopically marked when 15N was >15‰ or 13C was > —17‰. Thresholds are denoted with dashed lines; enriched samples are shown in black.

Table 1

Number of samples, average capture success (number of hair captures/trap night), and number of carcasses collected from commercial trappers for martens (Martes americanaand M. caurina) inhabiting sampling sites during autumn for Admiralty Island (ADM), Misty Fjords National Monument (MFNM), and Prince of Wales Island (POW), Alaska, and Queen Charlotte Islands (QCI), British Columbia, Canada.

Location Species No. samples No. trap occasions No. samples/no. trap occasions No. carcasses 
ADM M. caurina 231 760 0.304 49 
QCI M. caurina 246 0.008 71 
MFNM M. americana 45 425 0.106 78 
POW M. americana 94 1,650 0.057 168 
Location Species No. samples No. trap occasions No. samples/no. trap occasions No. carcasses 
ADM M. caurina 231 760 0.304 49 
QCI M. caurina 246 0.008 71 
MFNM M. americana 45 425 0.106 78 
POW M. americana 94 1,650 0.057 168 

Values of δD were correlated with both δ13C (r2 = 0.198, F1,63 = 15.3, P = 0.001) and δ15N (r2 = 0.203, F1,63 = 15.8, P < 0.001) among the unmarked samples, indicating that martens obtained substantial quantities of structural hydrogen for hair growth from their diet. Because both δ13C and δ15N were only weakly related to δD, however, we did not include estimates of carbon or nitrogen enrichment in identifying martens that had been marked with enriched deuterium. We detected strong differences among sites in the δD signature of unmarked martens (F2,63 = 43.4, P < 0.001; Fig. 5). Samples collected from the mainland (MFNM) were most depleted, followed by samples from the northernmost (ADM) and southernmost (POW and MFNM) coastal locations. We found that deuterium for marked animals (23 noninvasive samples and 48 carcasses) ranged from — 55.6‰ to — 15.1‰ in the continental MFNM sample, -37.5‰ to +2.2‰ in the ADM sample, and — 16.7‰ to +11.5‰ in the POW and coastal MFNM samples.

Fig. 5

Distribution of deuterium (‰) for marked and unmarked martens (Martes americana and M. caurina). Boxplots depict the mean, standard error, and 99% prediction interval for 64 unmarked martens from 3 sites in southeastern Alaska (Misty Fjords National Monument [MFNM], Prince of Wales Island [POW], and Admiralty Island [ADM]) that differed in natural abundance. Hair samples exceeding the 99% prediction interval were deemed marked with 2H, and are shown as circles (black = commercially trapped marten carcasses; white = noninvasively sampled martens). Hair samples deemed as unmarked are shown with hash marks.

Fig. 5

Distribution of deuterium (‰) for marked and unmarked martens (Martes americana and M. caurina). Boxplots depict the mean, standard error, and 99% prediction interval for 64 unmarked martens from 3 sites in southeastern Alaska (Misty Fjords National Monument [MFNM], Prince of Wales Island [POW], and Admiralty Island [ADM]) that differed in natural abundance. Hair samples exceeding the 99% prediction interval were deemed marked with 2H, and are shown as circles (black = commercially trapped marten carcasses; white = noninvasively sampled martens). Hair samples deemed as unmarked are shown with hash marks.

The majority of male and female carcass samples (19 Pacific martens and 59 American martens, 56 or 72% of 78) were obtained from the forest stand in which they were originally labeled and were within our threshold value for dispersal distance. However, we identified 22 individuals (3 Pacific martens and 19 American martens) as potential dispersers. Of these, 10 were males and 12 were females, but all were either young of the year (0 years old) or yearlings. Three dispersing individuals (2 Pacific martens and 1 American marten), all 0 years old, moved >25 km, presum- ably dispersing from their natal territories.

Although the exact distances moved by martens are incalculable without individual identity for the noninvasively collected samples, we were able to detect dispersal events and estimate the minimum Euclidean distance moved. Of the 86 marked hair samples (68 Pacific martens and 18 American martens) we identified, 62 were collected within the forest stand where the isotopic label was applied and were presumed to be resident animals. The remaining 24 samples (19 Pacific martens and 5 American martens) appeared to have been collected from long-distance dispersers, moving 15-40 km. The combination of carcasses and noninvasive samples suggested that Pacific and American martens exhibited similar rates of dispersal (28% of marked individuals dispersing; G = 0.70, P = 0.40), and similar dispersal distances (X̄ = 15.5 km ± 1.01 SE; t45 = 1.01, P = 0.16).

Discussion

Our results demonstrate that enriched isotopic baits combined with hair collection can be an effective method to quantify dispersal for midsized mammals. Using this ap- proach, and without additional information, we were able to quantify dispersal rates and distances, and explore the dispersal power of 2 species of martens with putatively different evolutionary histories. Both the rates and distances of dispersal we documented using enriched isotopes generally corresponded with those recently found using radiotelemetry among juvenile martens from Ontario (Johnson et al. 2009). Our ability to identify dispersers and quantify attributes of dispersal for American and Pacific martens using isotopic labeling and noninvasive sampling should be of use to future researchers, and particularly relevant for those studying Pacific marten populations on ADM. To date, identification of individuals and, therefore, estimation of dispersal were unattainable because of low diversity in microsatellite DNA loci (Small et al. 2003). Our findings that Pacific and American martens may possess similar dispersal propensity and dispersal power should be considered with caution; site- specific differences in habitat quality, marten density, and population dynamics are currently unknown for our study sites, and likely exert a strong influence on dispersal patterns.

Our field test of isotopic labeling revealed that some elements were more effective in marking the hair of martens. In particular, we found that 15N-labeled bait produced an unambiguous mark, especially among martens on ADM, whereas 13C and 2H marks were less obvious. This finding contrasts with our previous results from captive experiments, where consumption of a single dose yielded unequivocal marking in all 3 isotopes (Pauli et al. 2009). Our results for deuterium likely reflect both high spatial variation in values of this element isotope and high individual variation in isotopic incorporation. Wolf (2011) recently documented large indi- vidual variation in the incorporation of H into several types- of tissues, including feathers, in house sparrows (Passer domesticus) and Japanese quails (Corturnix japonica). Be- cause the birds were fed controlled diets, Wolf (2011) attributed some of the variation to individual differences in the amount of water consumed and metabolic activity. Our sampling locations on POW and MFNM, which share the same latitude (55° 38′N), exhibited significant differences in deuterium between inland and coastal sites in unmarked animals. Also, the observed individual variation in δD values among animals that we considered unmarked within each study site was large, ranging from — 99‰ to - 18‰. Although these data highlight some potential limitations of using natural abundance in dispersal studies, these factors also may influence our ability to use 2H in marking studies.

It is possible that the high δ15N values we observed among martens on ADM resulted from a few individuals consuming multiple doses of enriched bait. However, repeated consump- tion of bait cannot fully explain the differences between δ13C and δ15N among individuals on this island. During our sampling on ADM, small mammal abundance was low, but capture rates of martens were high (Table 1), suggesting that these animals could have been nutritionally stressed. When undernourished, predators will allocate some consumed amino acids to respiration (Ben-David et al. 2012; Martinez del Rio and Wolf 2005; Whiteman et al. 2012). This is particularly relevant for nonessential amino acids, including glycine. Thus, martens consuming enriched glycine may have allocated a portion of it to respiration. The deamination of glycine and subsequent nitrogen recycling could have led to an elevation in δ15N (Hobson et al. 1993; but see Ben-David et al. 1999) and a pronounced signature in ADM martens. To reduce the effects of elemental routing on isotopic labeling, we recommend the use of essential amino acids in field labeling studies. Although currently cost prohibitive, cysteine, a sulfur- bearing amino acid and a main component of keratin (Richards et al. 2003; Thomas et al. 2007) would likely be a superior vehicle for isotopic enrichment. Researchers explor- ing isotopic labels for free-ranging mammals, especially with 13C, also should consider increasing the amount of enriched amino acid applied to bait in the field to ensure that unambiguous marks are achieved.

Artificially enriching keratin will be particularly useful when other emerging approaches to quantifying dispersal rates, such as satellite tracking or spatial variation in the natural abundance of deuterium- or DNA-based analyses, are intractable. Animal movement data from global positioning system or ARGOS technology are rich in detail, but can be prohibitively expensive ($2,000–8,000 per collar [Hebble-white and Haydon 2010] versus $0.45–6.60 per isotopically enriched bait [Pauli et al. 2009]) and currently are limited to the study of large or midsized mammals. Natural deuterium levels are effective at reconstructing dispersal at continental scales (Hobson 1999; Wunder 2012), but provide insufficient resolution for species that disperse short to middistances (Wunder and Norris 2008). Dispersal inferred from genetic divergence among populations (e.g., FSTWright 1951) or via the assignment of individuals to subpopulations based on genotypes (Pritchard et al. 2000) fail for systems featuring a high degree of connectivity because allele frequencies are nearly identical among subpopulations (PalsbØll et al. 2010). Genetic kinship methods, characterizing dispersal via spatial distribution of close relatives (Peery et al. 2008), are effective for functionally connected systems, but require that the majority of the population is sampled (Palsboll et al. 2010). Thus, artificial isotopic enrichment seems particularly appli- cable for small-bodied species within systems that exhibit a high degree of functional connectivity and where the majority of the population cannot be sampled. Many of the collection techniques already developed for noninvasive, DNA-based approaches (e.g., hair traps [Pauli et al. 2008] or hair snares [DePue and Ben-David 2007]) can be used with the isotopic method we describe. Isotopic labeling offers the added benefit that, unlike DNA (Roon et al. 2003), the keratin is resilient to degradation (Macko et al. 1999).

For enriched isotopic baits to be effective, researchers will need to calculate dose requirements (Pauli et al. 2009), and conduct feeding trials to ensure that isotopic marks are unambiguous. In the field, enriched baits should be deployed during a period of keratin growth, ideally at the onset of a molt. Because most amino acids incorporated in hair are from recently consumed food items (Ayliffe et al. 2004), isotopic enrichment of ingested amino acids will provide a rapid mark within those tissues. Surprisingly, we found that a reasonable proportion of scat samples collected from martens on ADM were enriched with nitrogen or carbon (deuterium was not analyzed for feces). By using elemental abundance of nitrogen in fecal samples, we 1st separated animal-free scat samples from those with animal tissue. Based on the extreme isotopic signatures of plant and animal tissues, we then discriminated scat samples isotopically enriched in 15N and 13C. Although having reference unmarked fecal samples for comparison would have been ideal, our general approach in identifying marked scat seems promising, especially for detecting scats marked with 15N. Examination of these preliminary data suggests that collection of fresh scat may provide comple- mentary information on dispersal when collected concurrently with hair sampling. Finally, site- and tissue-specific baseline values for δD, δ13C, and δ15N are needed. Where these elements in the natural diet are not correlated, individuals that have consumed an enriched bait can be identified using the upper 99% prediction intervals constructed from baseline values. In cases where the correlation between elements in the diet is high, such as the relationship we observed for salmon- consuming martens, more involved analytical methods will be necessary to identify animals marked with enriched isotopes. Of course, without quantifying the individual identity for noninvasively sampled martens, our estimates of dispersal are prone to bias, particularly if dispersing animals possess different likelihoods of visiting baited traps or exhibit strong responses to sampling (e.g., becoming trap-happy or trap- shy—Pollock 1981).

Dispersal, the main process connecting animal subpopula- tions, is a critical aspect of a species' ecology, life history, and evolution. Knowledge of dispersal characteristics is particu- larly relevant given the current level of global habitat fragmentation (Foley et al. 2005) and forecasted consequences of climate change (Thomas et al. 2004). Although we developed this method for 2 species of mesocarnivores, providing artificially enriched isotopic bait and subsequently collecting hair or other tissue represents a safe (Koletzko et al. 1997) and cost-effective (Pauli et al. 2009) approach to mark and track small and medium-sized animals, thereby, enabling the quantification of this previously cryptic process.

Acknowledgments

We thank R. J. Artery, C. Reynolds, J. Pokallus, G. Campbell, B. Neiles, J. Thomas, and E. Flaherty for assisting with fieldwork and sample processing. R. O. Hall and C. Martínez del Rio provided valuable discussions on calculation of appropriate dosages. Financial support was provided by Pacific Northwest Research Station (United States Department of Agriculture Forest Service), a National Science Foundation Doctoral Dissertation Improvement Grant (0808177), a scholarship through the Program in Ecology at the University of Wyoming, and a fellowship from the American Society of Mammalogists.

Literature Cited

Ayliffe
L. K.
et al
.
2004
.
Turnover of carbon isotopes in tail hair and breath CO2 of horses fed an isotopically varied diet
.
Oecologia
 
139
:
11
22
.
Ben-David
M.
Flaherty
E. A.
.
2012
.
Stable isotopes in mammalian research: a beginner's guide
.
Journal of Mammalogy
 
93
:
312
328
.
Ben-David
M.
Flynn
R. W.
Schell
D. M.
.
1997
.
Annual and seasonal changes in diets of martens: evidence from stable isotope analysis
.
Oecologia
 
111
:
280
291
.
Ben-David
M.
Hanley
T. A.
Schell
D. M.
.
1998
.
Fertilization of terrestrial vegetation by spawning Pacific salmon: the role of flooding and predator activity
.
Oikos
 
83
:
47
55
.
Ben-David
M.
McColl
C. J.
Boonstra
R.
Karels
T. J.
.
1999
.
15N signatures do not reflect body condition in arctic ground squirrels
.
Canadian Journal of Zoology (London)
 
77
:
1373
1378
.
Ben-David
M.
Newsome
S. D.
Whiteman
J. P.
.
2012
.
Lipid and amino acid composition influence incorporation and discrimination of 13C and 15N in mink
.
Journal of Mammalogy
 
93
:
399
112
.
Boutton
T. W.
1991
.
Tracer studies with 13C-enriched substrates: humans and large animals
. Pp.
219
242
in
Carbon isotope techniques
  (
Coleman
D. C.
, ed.).
Academic Press Inc.
,
San Diego, California
.
Bowen
G. J.
Wassenaar
L. I.
Hobson
K. A.
.
2005
.
Global application of stable hydrogen and oxygen isotopes to wildlife forensics
.
Oecologia
 
143
:
337
348
.
Burton
R. K.
Koch
P. L.
.
1999
.
Isotopic tracing of foraging and long-distance migration in northeastern Pacific pinnipeds
.
Oecologia
 
119
:
578
585
.
Carlo
T. A.
Tewksbury
J. J.
Martínez del Rio
C.
.
2009
.
A new method to track seed dispersal and recruitment using 15N isotope enrichment
.
Ecology
 
90
:
3516
3525
.
Cryan
P. M.
Bogan
M. A.
Rye
R. O.
Landis
G. P.
Kestler
C L.
.
2004
.
Stable hydrogen isotope analysis of bat hair as evidence for seasonal molt and long-distance migration
.
Journal of Mammalogy
 
85
:
995
1001
.
den Boer
P. J.
1970
.
On the significance of dispersal power for populations of carabid-beetles (Coleoptera, Carabidae)
.
Oecologia
 
4
:
1
28
.
DePue
J. E.
Ben-David
M.
.
2007
.
Hair sampling techniques for river otters
.
Journal of Wildlife Management
 
71
:
671
674
.
Dockx
C.
Brower
L. P.
Wassenaar
L. I.
Hobson
K. A.
.
2004
.
Do North American monarch butterflies travel to Cuba? Stable isotope and chemical tracer techniques
.
Ecological Applications
 
14
:
1106
1114
.
Flaherty
E. A.
Ben-David
M.
.
2010
.
Overlap and partitioning of the ecological and isotopic niches
.
Oikos
 
119
:
1409
1416
.
Flynn
R. W.
Schumacher
T. V.
.
1999
.
Ecology of martens in Southeast Alaska
 .
Alaska Department Fish and Game
,
Juneau, Alaska
,
Federal Aid in Wildlife Restoration Progress Report, Grant W-27-2
,
Study 7
.
16
:
1
XX
.
Foley
J. A.
et al
.
2005
.
Global consequences of land use
.
Science
 
309
:
570
574
.
Graham
R. W.
Graham
M. A.
.
1994
.
Late Quaternary distribution of Martes in North America
. Pp.
26
58
in
Martens, sables, and fishers: biology and conservation
  (
Buskirk
S. W.
Harestad
A. S.
Raphael
M. G
Powell
R. A.
, eds.).
Cornell University Press
,
Ithaca, New York
.
Hall
R. O.
Tank
J. L.
.
2003
.
Ecosystem metabolism controls nitrogen uptake in streams in Grand Teton National Park, Wyoming
.
Limnology and Oceanography
 
48
:
1120
1128
.
Hebblewhite
M.
Haydon
D. T.
.
2010
.
Distinguishing technology from biology: a critical review of the use of GPS telemetry data in ecology. Philosophical Transactions of the Royal Society of London, B
.
Biological Sciences
 
365
:
2303
2312
.
Hirons
A. C
Schell
D. M.
St Aubin
D. J.
.
2001
.
Growth rates of vibrissae of harbor seals (Phoca vitulina) and Steiler sea lions (Eumetopias jubatus)
.
Canadian Journal of Zoology
 
79
:
1053
1061
.
Hobson
K. A.
1999
.
Tracing origins and migration of wildlife using stable isotopes: a review
.
Oecologia
 
120
:
314
326
.
Hobson
K. A.
Alisauskas
R. T.
Clark
R. G
.
1993
.
Stable-nitrogen isotope enrichment in avian tissues due to fasting and nutritional stress: implications for isotopic analyses of diet
.
Condor
 
95
:
388
394
.
Hobson
K. A.
Wassenaar
L. I.
.
1997
.
Linking breeding and wintering grounds of Neotropical migrant songbirds using stable hydrogen isotopic analysis of feathers
.
Oecologia
 
109
:
142
148
.
Johnson
C. A.
Fryxell
J. M.
Thompson
I. D.
Baker
J. A.
.
2009
.
Mortality risk increases with natal dispersal distance in American martens. Proceedings of the Royal Society of London, B
.
Biological Sciences
 
276
:
3361
3367
.
Koletzko
B.
Sauerwald
T.
Demmelmair
H.
.
1997
.
Safety of stable isotope use
.
European Journal of Pediatrics
 
156
:
S12
S17
.
MacDonald
S. O.
Cook
J. A.
.
2009
.
Recent mammals of Alaska
 .
University of Alaska Press
,
Fairbanks, Alaska
.
Macko
S. A.
Engel
M. H.
Andrusevich
V.
Lubec
G.
O'Connell
T. C.
Hedges
R. E. M.
.
1999
.
Documenting the diet in ancient human populations through stable isotope analysis of hair. Philosophical Transactions of the Royal Society of London, B
.
Biological Sciences
 
354
:
65
76
.
Macneale
K. H.
Peckarsky
B. L.
Likens
G. E.
.
2004
.
Contradictory results from different methods for measuring direction of insect flight
.
Freshwater Biology
 
49
:
1260
1268
.
Martínez del Rio
C.
Wolf
B. O.
.
2005
.
Mass-balance models for animal isotopic ecology
. Pp.
141
174
in
Physiological and ecological adaptations to feeding in vertebrates
  (
Starck
J. M.
Wang
T.
, eds.).
Science Publishers, Enfield
,
New Hampshire
.
Newsome
S. D.
Yeakel
J. D.
Wheatley
P. V.
Tinker
M. T.
.
2012
.
Tools for quantifying isotopic niche space and dietary variation at the individual and population level
.
Journal of Mammalogy
 
93
:
329
341
.
PalsbØll
P. J.
Peery
M. Z.
Bérubé
M.
.
2010
.
Detecting populations in the ‘ambiguous’ zone: kinship-based estimation of population structure at low genetic divergence
.
Molecular Ecology Resources
 
10
:
797
805
.
Pauli
J. N.
Ben-David
M.
Buskirk
S. W.
DePue
J. E.
Smith
W. P.
.
2009
.
An isotopic technique to mark mid-sized vertebrates non-invasively
.
Journal of Zoology (London)
 
278
:
141
148
.
Pauli
J. N.
Hamilton
M. B.
Crain
E. B.
Buskirk
S. W.
.
2008
.
A single-sampling hair trap for mesocarnivores
.
Journal of Wildlife Management
 
72
:
1650
1652
.
Peery
M. Z.
et al
.
2008
.
Characterizing source-sink dynamics with genetic parentage assignments
.
Ecology
 
89
:
2746
2759
.
Pollock
K. H.
1981
.
Capture-recapture models: a review of current methods, assumptions and experimental design
. Pp.
426
435
in
Estimating numbers of terrestrial birds
  (
Ralph
C. J.
Scott
M. J.
, eds.).
Studies in Avian Biology. Cooper Ornithological Society
,
Lawrence, Kansas
.
Pritchard
J. K.
Stephens
M.
Donnelly
P.
.
2000
.
Inference of population structure using multilocus genotype data
.
Genetics
 
155
:
945
959
.
Richards
M. P.
Fuller
B. T.
Sponheimer
M.
Robinson
T.
Ayliffe
L.
.
2003
.
Sulphur isotopes in paleodietary studies: a review and results from a controlled feeding experiment
.
International Journal of Osteoarchaeology
 
13
:
37
5
.
Roon
D. A.
Waits
L. P.
Kendall
K. C.
.
2003
.
A quantitative evaluation of two methods for preserving hair samples
.
Molecular Ecology Notes
 
3
:
163
166
.
Rubenstein
D. R.
Hobson
K. A.
.
2004
.
From birds to butterflies: animal movement patterns and stable isotopes
.
Trends in Ecology & Evolution
 
19
:
256
263
.
Schell
D. M.
Saupe
S. M.
Haubenstock
N.
.
1988
.
Natural isotope abundance in bowhead whale (Balaena mysticetus) baleen: markers of aging and habitat usage
. Pp.
260
269
in
Stable isotopes in ecological research
  (
Rundel
P. W.
Ehleringer
J. R.
Nagy
K. A.
, eds.).
Ecological Studies 68
.
Springer-Verlag
,
Berlin, Germany
.
Sikes
R. S.
Gannon
W. L.
the Animal Care,Use Committee of the American Society of Mammalogists
.
2011
.
Guidelines of the American Society of Mammalogists for the use of wild mammals in research
.
Journal of Mammalogy
 
92
:
235
253
.
Small
M. P.
Stone
K. D.
Cook
J. A.
.
2003
.
American marten (Martes americana) in the Pacific Northwest: population differentiation across a landscape fragmented in time and space
.
Molecular Ecology
 
12
:
89
103
.
Stone
K. D.
Cook
J. A.
.
2002
.
Molecular evolution of Holarctic martens (genus Martes, Mammalia: Carnivora: Mustelidae)
.
Molecular Phylogenetics and Evolution
 
24
:
169
179
.
Stone
K. D
Flynn
R. W.
Cook
J. A.
.
2002
.
Post-glacial colonization of northwestern North America by the forest-associated American marten (Martes americana, Mammalia: Carnivora: Mustelidae)
.
Molecular Ecology
 
11
:
2049
2063
.
Thomas
C. D.
et al
.
2004
.
Extinction risk from climate change
.
Nature
 
427
:
145
148
.
Thomas
N.
Tivey
D. R.
Penno
N. M.
Nattrass
G.
Hynd
P. I.
.
2007
.
Characterization of transport systems for cysteine, lysine, alanine, and leucine in wool follicles of sheep
.
Journal of Animal Science
 
85
:
2205
2213
.
Wanner
H.
Gu
H.
Hattendorf
B.
Günther
D.
Dorn
S.
.
2006
.
Using the stable isotope marker 44Ca to study dispersal and host-foraging activity in parasitoids
.
Journal of Applied Ecology
 
43
:
1031
1039
.
Wassenaar
L. I.
Hobson
K. A.
.
2000
.
Improved method for determining the stable-hydrogen isotopic composition (δD) of complex organic materials of environmental interest
.
Environmental Science and Technology
 
34
:
2354
2360
.
Whiteman
J. P.
Greller
K. A.
Harlow
H. J.
Felicetti
L. A.
Rode
K. D.
Ben-David
M.
.
2012
.
Carbon isotopes in exhaled breath track metabolic substrates in brown bears (Ursus arctos)
.
Journal of Mammalogy
 
93
:
413
421
.
Wolf
N.
2011
.
An experimental exploration of the use of hydrogen and oxygen stable isotopes in animal ecology
.
Ph.D. dissertation
 ,
University of Wyoming
,
Laramie, Wyoming
.
Wright
S.
1951
.
The genetical structure of populations
.
Annals of Eugenics
 
15
:
323
354
.
Wunder
M. B.
2012
.
Determining geographic patterns of migration and dispersal using stable isotopes in keratins
.
Journal of Mammalogy
 
93
:
360
367
.
Wunder
M. B.
Norris
D. R.
.
2008
.
Improved estimates of certainty in stable-isotope-based methods for tracking migratory animals
.
Ecological Applications
 
8
:
549
559
.
Zar
J. H.
1999
.
Biostatistical analysis
 .
4
th ed.
Prentice Hall
,
Upper Saddle River, New Jersey
.

Author notes

Present address of JNP, Department of Forest and Wildlife Ecology, University of Wisconsin, 1630 Linden Drive, Madison, WI 53706, USA
Special Feature Editor was Samantha M. Wisely.