The research procedure was approved by the Office of Research Ethics and Integrity of the University of Ottawa (Ethics File number: H10-17-10). Specifically, all sampling and analytical methods used were in accordance with relevant guidelines from the Office of Research Ethics and Integrity of the University of Ottawa. Written consents were obtained from all participants in accordance with and maintained under regulations from the Office of Research Ethics and Integrity of the University of Ottawa.

The hair sample was collected from a volunteer in on November 20, 2010. The travel history of the volunteer is shown in Table 1. The volunteer primarily resided in Ottawa, Canada and traveled to several distant locations during the timeframe of this study (May 2009 to November 2010). The four major trips longer than a week were (1) a trip to the Prairie provinces of Canada between June 30 and July 27, 2010; (2) a trip to Europe between May 2 and May 17, 2010; (3) a trip to Cuba between January 29 and February 5, 2010; and (4) a trip to the Eastern provinces of Canada between June 15 and July 24, 2009. Figure 1 shows the regions of four major trips of the individual outside Ontario on a world map.

One bundle of approximately 80 hairs was provided by the volunteer to the lab. It was cut with metal scissor near the root of hair on November 20, 2010.

∼50 individual hairs were aligned by the root and were taped on weighing paper. The aligned hairs were cut every half centimeter using an in-house built cutting apparatus. The 0–0.5 cm segment was kept with the adhesive tape that was used to hold the orientation of the hair bundle. But the adhesive tape residue and the trapped dust may contaminate isotope analysis. Thus the 0–0.5 cm segment was not analyzed. 38 samples were obtained from the hair of ∼19.5 cm long. Each hair sample was washed for 15 min with gentle agitation in a 2:1 chloroform/methanol (v/v) solution to remove surface contamination (Font et al., 2012; Tipple et al., 2018). After drying, each sample was analyzed for δ¹³C, δ¹⁵N, and δ²H. An aliquot of each sample was saved for ⁸⁷Sr/⁸⁶Sr analysis. Sr abundance was too low for ⁸⁷Sr/⁸⁶Sr analysis at half-centimeter resolution. Consequently, we combined at least two 0.5 cm segments into one sample.

For δ¹³C and δ¹⁵N analysis, the hair segment samples and isotope standards (approximately 400 μg each) were weighed into tin capsules and loaded onto an Elemental Analyser (Vario EL III, Elementar, Germany) interfaced with a ConFlo III (Thermo, Germany) to an isotope ratio mass spectrometer (IRMS, Delta^Plus Advantage, Thermo, Germany). One tin capsule was analyzed for each sample. Internal standards used for calibration were a mix of ammonium sulfate and sucrose (δ¹³C_VPDB, −11.94‰; δ¹⁵N_AIR, 16.58‰), nicotinamide (δ¹³C_VPDB, −22.95‰; δ¹⁵N_AIR, 0.07‰), and caffeine (δ¹³C_VPDB, −28.53‰; δ¹⁵N_AIR, −3.98‰). All δ¹⁵N and δ¹³C values are reported versus AIR and V-PDB, respectively. Analytical precision is based on an internal QC sample (glutamic acid, which is not used for normalization) and is usually better than ±0.2‰ for both δ¹³C and δ¹⁵N.

For δ²H analysis, the hair segment samples and isotope standards (approximately 300 μg each) were weighed into silver capsules and loaded onto a Pyrolysis Elemental Analyser with a zero-blank autosampler (TC/EA, Thermo, Germany) interfaced with a ConFlo IV (Thermo, Germany) to an IRMS (Delta^Plus XP, Thermo, Germany). Due to exchangeable hydrogen in the cut hair samples, an exchange experiment was performed on the cut hair segments to determine the “true” δ²H value, following the procedure in Meier-Augenstein et al. (2011). This procedure requires two aliquots of the same hair sample, one equilibrated in isotopically “light” water (δ²H_VSMOW, -398‰), and the other in isotopically “heavy” water (δ²H_VSMOW, + 15.6 ‰). The segmented hair samples were normalized to two materials with non-exchangeable H: IAEA-CH-7 (δ²H_VSMOW, −100.3‰), and an in-house kaolinite (δ²H_VSMOW, −58.0‰). Four ground human hair samples were used as QC (Meier-Augenstein et al., 2011): AND (δ²H_VSMOW, −71.6‰), COL (δ²H_VSMOW, −87.9‰), CAL/CAN (δ²H_VSMOW, −106.4‰) and CAL/SAL (δ²H_VSMOW, −101.2‰). All reported δ²H values are reported to the international scale VSMOW. Analytical precision for δ²H was ±2.0‰. For all δ²H hair analyses, the production of HCN was not accounted for Gehre et al. (2015).

All samples (approximately 2 mg each) for ⁸⁷Sr/⁸⁶Sr analysis were digested in 7 mL cleaned Xpress Teflon vials using a mixture of 200 μL hydrogen peroxide ≥ 30% (for trace analysis Grade; Sigma-Aldrich, Italy) and 3 mL concentrated nitric acid (HNO₃) (TraceMetal^TM Grade; Fisher Chemical, Canada) for 15 min at 150°C using the CEM MARS 5 Microwave Accelerated Reaction System. After digestion, the solutions were dried on a hot plate at 100°C in a laminar flow fume hood. When there is only a drop left, each vial was removed from the hotplate and the bottom of the vial was observed to check if the liquid was limpid. If there were still white particles, suggesting uncompleted digestion, further digestion was done by adding hydrogen peroxide to digest all the organic residue. If there was no more particle in the final drop of solution, 1 mL concentrated HNO₃ was added to each vial and then transferred to a 7 mL Savillex PTFE vial. An aliquot of 50 μL of the solution from each Savillex vial was pipetted to a Labcon MetalFree^TM centrifuge tubes and diluted with 2 mL of 2% v/v HNO₃. The diluted solution was then centrifuged for 10 min before Sr concentration analysis. Sr concentration analysis was performed by Inductively Coupled Plasma Mass Spectrometry (ICP-MS) (Agilent 8800 triple quadrupole mass spectrometer) at the Department of Earth and Environmental Sciences, University of Ottawa. Calibration standards were prepared using single element certified standards purchased from SCP Science (Montreal, Canada).

The remaining ∼1 mL aliquot of the sample in the 7 mL Savillex PTFE vial was dried down and re-dissolved in 1 mL 6 M HNO₃. The separation of Sr was processed in 100 μL microcolumn loaded with Sr-spec Resin^TM (100–150 μm; Eichrom Technologies, LLC). The matrix was rinsed out using 6 M HNO₃. The Sr was collected with 0.05 M HNO₃. After separation, the eluates were dried and re-dissolved in 200 μL 2% v/v HNO₃ for ⁸⁷Sr/⁸⁶Sr analysis. The ⁸⁷Sr/⁸⁶Sr analysis was performed at the Isotope Geochemistry and Geochronology Research Centre, Carleton University using a Thermo Scientific^TM Neptune^TM high-resolution multi-collector inductively coupled plasma mass spectrometer (MC-ICP-MS; Thermo Fisher Scientific, Bremen, Germany). Sample solutions in 2% v/v HNO₃ were introduced using a microFAST MC single-loop system (Elemental Scientific). The loading volume was 200 μL and the injecting rate was 30 μL/min. The solution was aspirated using a PFA nebulizer, a double-pass quartz spray chamber, quartz torch, and nickel sample and skimmer cones. Isotopes ⁸²Kr, ⁸³Kr, ⁸⁴Sr, ⁸⁵Rb, ⁸⁶Sr, ⁸⁷Sr, and ⁸⁸Sr were simultaneously measured in L4, L3, L2, L1, C, H1, and H2 Faraday cups, respectively. Measurements of samples were made using a static multi-collector routine that consisted of one block of 75 cycles with an integration time of 4.194 s/cycle. ⁸⁴Sr and ⁸⁶Sr have isobaric interferences from ⁸⁴Kr and ⁸⁶Kr, respectively. ⁸⁷Sr has an isobaric interference from ⁸⁷Rb. The interferences of ⁸⁴Sr and ⁸⁶Sr were corrected by subtracting the amount of ⁸⁴Kr and ⁸⁶Kr corresponding to the ⁸³Kr signal. The interference of ⁸⁷Sr was corrected by subtracting the amount of ⁸⁷Rb corresponding to the ⁸⁵Rb signal. Instrumental mass fractionation was corrected by normalizing ⁸⁶Sr/⁸⁸Sr to 0.1194 using the exponential law. Strontium isotope compositions are reported as ⁸⁷Sr/⁸⁶Sr ratios. The long-term reproducibility of the ⁸⁷Sr/⁸⁶Sr measurement for NIST SRM987 is 0.71025 ± 0.00004 (1 SD, n = 16). A 100 ng/g pure Sr standard was measured along with the sample as the in-house standard (0.70823 ± 0.00005, 1 SD, n = 4). The long-term reproducibility of the in-house standard is (0.70822 ± 0.00004, 1 SD, n = 106).

We compiled the ranges of δ²H values, δ¹³C values, δ¹⁵N values, and ⁸⁷Sr/⁸⁶Sr ratios of each location from our database, the Forensic Attribution Database (FAD; Chartrand and St-Jean, 2015; Bataille et al., 2020a), completed by literature data (Table 1). The FAD gathers isotope data in hair from 600 Canadian hair donors and tap water samples collected across Canada (Chartrand and St-Jean, 2015; Bataille et al., 2020a). In this study, we compiled those isotopic ranges from hair samples of the locations (or the nearest locations) visited by our volunteer from May 1, 2009 to Nov 20, 2010.

We estimated the isotopic ranges of H, C, and N in Canadian locations from the hair samples in the FAD. To estimate the isotopic ranges of H, C, and N in the locations visited in Europe, Cuba, and the United States, we used existing compilations of isotope data for hair (Ehleringer et al., 2008; Hülsemann et al., 2015; Lehn et al., 2015b) (Table 1).

For Sr, we have not yet measured the ⁸⁷Sr/⁸⁶Sr of the hair samples in the FAD. Data from the literature are also too sparse to characterize the ⁸⁷Sr/⁸⁶Sr of the traveled locations either. Therefore, we cannot directly estimate the range of hair ⁸⁷Sr/⁸⁶Sr for the studied individual. Theoretically, we could use the ⁸⁷Sr/⁸⁶Sr range of local food to estimate the possible ⁸⁷Sr/⁸⁶Sr profile, considering more than 90% of Sr in mammal’s body is from food (Lewis et al., 2017). However, because of the complexity of the modern supermarket diet, it is almost impossible to compile the ⁸⁷Sr/⁸⁶Sr ranges of local foods from all traveled locations. As previous studies have shown that the hair samples usually have the same ⁸⁷Sr/⁸⁶Sr as the local tap water (Tipple et al., 2019; Ammer et al., 2020), we chose the ⁸⁷Sr/⁸⁶Sr ratios of the tap water samples in the FAD (Chartrand and St-Jean, 2015) and the river or tap water samples in the literature (Voerkelius et al., 2010; Chesson et al., 2012; Bataille et al., 2020b) to estimate the range of hair ⁸⁷Sr/⁸⁶Sr for the studied individual. The assumption behind such estimation, that the ⁸⁷Sr/⁸⁶Sr profile in hair followed the ⁸⁷Sr/⁸⁶Sr of the local water, is not ideal, but it is the best we can do in this study.

The isotopic ranges compiled in Table 1 are unlikely to cover the entire range of isotopic variability at each of the visited locations but they provided a strong empirical basis to approximate these values.

In this study, we know both the travel history of the studied individual and the ranges of the isotopic values of hair samples from local residents at each traveling location (Table 1). We aim to use those data to reconstruct the expected isotopic variations in the hair of the traveling individual. One issue, however, is that the stable isotope signal of the newly grown hair segment does not reflect 100% of the local foods and drinks, because there are different metabolic pools with different turnover rates. For example, Ayliffe et al. (2004) found that, following a dietary change from C3 to C4 food, ∼41% of carbon isotope signal of the horse tail hair was from a pool (probably from the new diet) with a very fast turnover rate (half-life of ∼0.5 days), ∼15% from a pool with an intermediate turnover rate (half-life of ∼4 days), and ∼44% from a pool with very slow turnover rate (half-life of ∼140 days). In our case study, as an individual moves on the landscape, the isotopic composition of the diet is constantly changing, which complicates the reconstruction of the expected isotopic variations in the hair profile.

Several methods can be used to account for isotope turnover in hair. Traditionally, researchers used a single exponential best-fit model to describe isotope turnover in biological tissues (Podlesak et al., 2005). Since fitting data to an exponential function cannot ascertain the possibility of multiple turnover pools, Cerling et al. (2007) introduced the reaction progress model to determine the turnover rate. It has been shown that some elements have multiple turnover pools in the body, and their isotopic profiles in hair can be better reconstructed with a reaction progress model (Cerling et al., 2007; Podlesak et al., 2008). However, constructing a multiple-pool model requires a highly controlled feeding experiment and high-resolution isotopic analysis of hair to calibrate the parameters. As far as we know, there has not been such a feeding experiment on human beings. In this study, we do not have sufficiently accurate dietary and metabolic information to calibrate this type of multiple pool model. In addition, the temporal resolution in this study (∼2 weeks for hydrogen, carbon, nitrogen and ∼4 weeks for strontium) was much less precise than the dietary studies with short sampling interval (∼1 day or less) (e.g., Cerling et al., 2007). Lastly, because we were looking at a bundle of hairs growing at different rates and different phases (LeBeau et al., 2011), it can lead to further noise in the turnover rate modeling, which we cannot address in our model. Therefore, we assume that a single exponential best-fit model is a good-enough first-order approximation to predict isotopic variations in the hair profile. We represented isotope turnover time using the following equation of an exponential decay function (Cerling et al., 2007):

where δAi⁢n⁢i⁢t, δAe⁢q, and δAt are the δ-values of element A at the initial time (t = 0), at equilibrium (t = infinity), and at time t, respectively. We calculated the λ value for this equation using existing data from the literature for hydrogen isotopes (Ehleringer et al., 2008) and we obtain a λ value of 0.04 day ^–1 (half-life of 17.3 days), which was comparable to the first-stage 19 days half-life measured on woodrats (Neotoma cinerea) hair (Podlesak et al., 2008). We use this calibrated λ value for all the isotopic systems used in this study. While it is likely that different elements have different λ values, changing the λ value within a reasonable range has little influence on the general pattern and the interpretation of the isotopic profiles. Please see the sensitivity analysis in the Supplementary Material for more information.

Once calibrated, we used this λ value and the isotope forensic attribution database to reconstruct the isotopic variations in the hair through time. We use the following equation to calculate the expected isotope variations in the hair of the traveling individual at daily resolution:

where δ is the isotopic value of hydrogen (δ²H), carbon (δ¹³C), or nitrogen (δ¹⁵N).

For ⁸⁷Sr/⁸⁶Sr ratios, this equation becomes:

where R is the isotopic ratio of strontium, ⁸⁷Sr/⁸⁶Sr. As discussed in the last section, we did not have the ranges of hair ⁸⁷Sr/⁸⁶Sr in the traveled locations. We predicted the hair ⁸⁷Sr/⁸⁶Sr profile using the ranges of ⁸⁷Sr/⁸⁶Sr of tap water or river water samples in the traveled locations. The impact of such approach is evaluated in the discussion.

The minimum and maximum of isotopic values for every day is calculated from May 1, 2009 to Nov 20, 2010.

The δ²H, δ¹³C, and δ¹⁵N values of the hair segments are shown in Table 2. The δ²H values in the hair ranged from −85 to −75‰. Except for four hair segments between A8 (Sample ID) and A11, the majority of the segments fall between −80 to −75‰. The δ¹³C values in the hair ranged from −18.4 to −17.0‰. While most of the δ¹³C values of the segments from A10 to A38 were relatively stable, the δ¹³C value of the segment A16 was much lower (−18.4‰) than the rest of the hair. The δ¹³C values of the segments from A1 to A9 were higher than the rest of the hair. The δ¹⁵N values in the hair ranged from 8.8 to 11.9 ‰. The first 13 cm of the hair (from A1 to A26) had a relatively narrow range of δ¹⁵N values (9.6–10.1‰) except for the segment A1 (10.5‰), A22 (10.5‰), and A25 (10.4‰), which contained slightly higher δ¹⁵N values. The remaining hair (A27–A38) had a much larger variation of δ¹⁵N value (Table 2).

The Sr concentrations ([Sr]) and isotopic ratios (⁸⁷Sr/⁸⁶Sr) of the hair segments are shown in Table 3. The [Sr] increases from the root to the tip of the hair from 1.6 mg/kg to 12.4 mg/kg. The ⁸⁷Sr/⁸⁶Sr ratios in the segments range from 0.7110 to 0.7119. The segment B10 was lost during the isotope analysis and is thus not reported. The first two segments (B1 and B2) had lower ⁸⁷Sr/⁸⁶Sr ratios than all other segments.

A prerequisite prior to comparing the observed and modeled hair isotopic profiles is to convert hair length into time. Hair growth rate can vary a lot among individuals. The average growth rate of human scalp hair is about 1 cm/month but can range from 0.6 to 2.2 cm/month (Pötsch, 1996). As it is difficult to determine the growth rate of an individual after hair collection, we matched the most prominent δ²H peaks between the observed and predicted hair profiles. We found a growth rate of 1.5 cm per 30 days provided the best fit between observed and predicted peaks for δ²H values but also for δ¹³C values. However, the growth rate chosen in this study may not be suitable for other studies. We also added 0.8 cm to the distance of the hair segment to the hair root (Tables 2, 3) to make up the ∼0.5 cm of hair remaining at the scalp during cutting (LeBeau et al., 2011; Lehn et al., 2019) and ∼0.3 cm of hair residue left in the hair follicle (Jimenez et al., 2011). The “distance from hair root” is then used for the comparison between the measured and the predicted isotope profiles.

In this study, we used a bundle of hair to measure isotopic variations in the hair profile of the volunteer. As mentioned above, we assumed that all hair in the bundle had a constant growth rate of 1.5 cm per 30 days. In reality, scalp hairs of the same individual have different growth rates as well as different growth phases (Hayashi et al., 1991; LeBeau et al., 2011). The isotopic variability measured on multiple hairs will become progressively attenuated due to the mixing of hair segments from slightly different life periods. This issue was underlined by Remien et al. (2014) who compared the isotopic signal measured in a single hair with that of multiple hairs and found that the signal was progressively blurred and difficult to use for provenance applications. Consequently, in our study, the isotopic variations identified close to the base of the hair are likely more representative of the true isotopic variability of the water or diet. Isotopic variations identified in the distal portion of the hair have larger uncertainty in terms of timing and amplitude.

Interestingly, Lehn et al. (2019) have found that the change of isotope signal on bundled hair profile is usually provoked by the fastest-growing hairs in the hair bundle because those fast-growing single hairs tend to be thicker and heavier thus their isotopic signals outweigh the others. Considering a commonly used mean hair growth rate of 1.1 cm/month and a proposed ± 30% intra-individual variability in single hair growth rate, they found that the growth rate of 1.4 cm per 30 days may be the best approach to calculate the exact dates of the isotopic shifts along the hair profile. Therefore, although our chosen value of 1.5 cm per 30 days is the best estimate for the hair growth rate that was used to interpret the isotope profiles, it is likely to represent the fastest growth rate, instead of the average hair growth rate, of the studied individual. Lehn et al. (2019) included the uncertainty of intra-individual hair growth rate in their calculation for the expected δ¹³C values along the hair strands. In this study, we calculate the predicted isotope values using a fixed growth rate and discuss the impact of the uncertainty of growth rate on our results in the following section.

The majority of the measured δ²H values in the hair profile of the individual falls within a tight range between −75 and −78‰ (Figure 2). Those values are within the range of observed δ²H values for human hair in the Ottawa region (Table 1). This relatively constant isotopic composition is found during periods when the volunteer was residing in the Ottawa region. This limited variability demonstrates that in the absence of large movements the δ²H values remain relatively stable and are within the expected range from the city.

We observed a major shift in δ²H values between July and September 2010 (Figure 2). This shift matches well with the predicted isotopic variations and corresponds to the 4 weeks long road trip to the Prairie provinces (Manitoba, Saskatchewan, and Alberta). The Prairie provinces are located inland in comparison with Ottawa, and tap waters in those regions have much lower δ²H values (Bowen, 2010). As the δ²H value in hair mostly (∼80%) reflects the food (Topalov et al., 2019), the remarkably low δ²H value (as low as −114‰) in the hair of Prairie residents (Table 1) indicates a distinctly low δ²H signal in the local diet. As the volunteer ate local food and drank local water this isotopic difference is transmitted to the hair. We notice that the amplitude of the observed δ²H shift is slightly attenuated in comparison with the predicted values (Figure 2) which likely reflects signal attenuation associated with combining multiple strands of hair.

The volunteer also traveled to central Europe in May 2010 for 17 days and this trip is reflected in the measured isotopic profile. The volunteer traveled to Washington, DC for 6 days during mid-April 2010. Because of the distinct low δ²H values in Washington, DC (−94 to −90‰) compared to Ottawa (Table 1), a sharp negative shift of δ²H value is shown in the predicted profile. A similar negative shift of δ²H is shown in the measured profile near 11–12 cm from the hair root (Figure 2), which is about 2 weeks earlier than the April shift in the predicted profile. As mentioned above, signals of the hair bundle distant from the hair root will be attenuated, and their timing will be more uncertain due to the differential growth rate of hair (Hayashi et al., 1991; LeBeau et al., 2011; Remien et al., 2014).

The volunteer traveled to Cuba for 7 days in January and February 2010. This trip is associated with a prominent positive shift in δ²H values in the predicted profile but only a small increase in the observed profile (Figure 2). The amplitude of the measured isotopic signal is attenuated in comparison with the predicted values. However, the increase of the lower limit of the predicted δ²H profile during the Cuba trip was not as prominent as its upper limit, allowing the possibility of a small increase of the observed δ²H values in hair. We note that signals found in bundle hair > 10 cm can still provide some indication of travel.

During summer 2009, the volunteer traveled to the Eastern provinces (New Brunswick, Newfoundland and Labrador, Nova Scotia, Prince Edward Island, and Quebec) for about 40 days. The Eastern provinces have lower δ²H values than Ottawa as reflected by the predicted δ²H values. The measured δ²H values progressively decrease toward the Eastern isotopic values. However, such decrease is longer and does not match well with the timing of travel as shown in the predicted profile (Figure 2). This signal in distal part of hair has a high uncertainty in terms of timing (Lehn et al., 2019) and thus is probably highly attenuated by the isotopic signal from Ottawa.

As the turnover rate was calculated based on the measured δ²H values in hair from Ehleringer et al. (2008), the pattern matching of δ²H profiles seems to be the best among all four isotopic systems in this study. The trip to the Prairie provinces significantly decreased the δ²H values by about 10 ‰, which is shown in both measured and predicted profiles (Figure 2). The increase around January and February 2010 during Cuba trip and the decrease around March and April in 2010 during the Washington DC trip are also shown in both profiles. This suggests that hydrogen isotope analysis of hair may be helpful to reconstruct the travel history of highly mobile individuals for forensic applications of modern fauna and human, when the traveled locations have distinct δ²H values (Ehleringer et al., 2008; Mant et al., 2016).

When using a hair bundle, quantitative geolocation is possible for an isotopic signal located close to the base of the hair (<10 cm). For isotopic signal located in a more distal portion of the hair, we suggest that only qualitative interpretation of travel can be provided as the timing and attenuation become significant issues. As the forward hair bundle model in Remien et al. (2014) shows and the isotope data along human hair strands in Lehn et al. (2019) show, the uncertainty of timing in the distal portion of hair can be ±2 months or more. Keratinous tissues can be preserved for thousands of years, however, using δ²H values in archeological samples requires additional caution as δ²H values may be altered by exchanging hydrogen with the depositional environment (e.g., Saul, 2017; Kootker et al., 2020).

The majority of the measured δ¹³C values in the hair profile of the individual falls within a tight range between −18 to −17.5‰, except for the segment A16 and the segments A2–A9 (Table 2). Interestingly, none of the measured δ¹³C values falls inside the predicted profile, although the measured and predicted profiles share very similar patterns (Figure 3). There is a relatively constant 1.5‰ offset between the measured and predicted δ¹³C profiles. For example, during the time of October to December 2009 when the volunteer stayed at Ottawa, the measured δ¹³C was around −17.9‰ compared with the Ottawa local range of −18.9 to −18.6‰ (Table 1). Such offset demonstrates the existence of individual variability of hair δ¹³C values from the local ranges, which may be caused by personal dietary preference (Hobson and Clark, 1992; Mekota et al., 2006; Petzke et al., 2010). The high δ¹³C values in the hair of this study may indicate a higher intake of animal or marine-based protein in the diet. Such discrepancy between measured and predicted individual value suggests that the absolute δ¹³C values in hair are unreliable to identify local versus non-local samples because δ¹³C values in hair can be largely influenced by individual preferences of certain foods or special diets. However, the strong correspondence between the measured and predicted profiles also suggest that δ¹³C values are useful to assess mobility. For example, the large, negative excursion in δ¹³C values in the predicted profile associated with the trip to Europe in May 2010 matches well with the negative shift of the measured profile in the hair segment A16 (Table 2 and Figure 3). The difference in δ¹³C values in food consumed in Europe and North America is well-established (Valenzuela et al., 2012; Hülsemann et al., 2015). Foods in North America, including livestock, tend to have much higher δ¹³C values than those of Europe because of the higher proportion of C₄ plant such as corn or cane sugar (with higher δ¹³C value) within food systems. In Europe, livestock are primarily fed with C₃ plants (with lower δ¹³C value) and the lower δ¹³C baseline of food system propagate within the population (Valenzuela et al., 2012; Hülsemann et al., 2015). Movement from the US to Europe usually leads to a significant decrease of δ¹³C value in hair (McCullagh et al., 2005). Therefore, the negative shifts of δ¹³C value in both the measured and the predicted profile of the volunteer correspond to a change in the isotopic baseline values in food consumed during the trip to Europe in May 2010.

While the European trip is easily identifiable in the measured and predicted δ¹³C profiles, other multi-week trips show some significant variations. The trip to the Prairie provinces between June and July 2010 is associated with a negative shift of measured δ¹³C values around 6–7 cm in the hair (A10–A11) (Figure 3) as some locations in the Prairie provinces have as low as −19.6‰ of the δ¹³C value (Supplementary Table S1), lower than the Ottawa range (−18.9 to −18.6‰) (Table 1). Again, this difference is associated with different isotopic baseline in food systems with western provinces growing and incorporating more C₃ crops. For this trip to western provinces, the amplitude of the shift in δ¹³C values is less marked than the shift of δ²H values (Figure 2). However, the combined use of δ¹³C and δ²H values can strongly help constrain the location of travel within Canada. The trip to Cuba in early 2010 did not cause significant fluctuation of δ¹³C values in the measured or predicted profile because the δ¹³C values in Cuba, although wider range, have similar average value as Ottawa (Table 1). There is an increase in δ¹³C values observed in the last two segments toward the tip of the volunteer’s hair which likely correspond to the trip to the Eastern provinces as shown in the predicted profile (Figure 3). As observed for δ²H value, there is a mismatch between the timing of the observed and predicted shifts for the trip to Eastern provinces because those isotopic signals are located in the distal portion of the hair.

There is also an increase of measured δ¹³C value in the hair segments A2–A9 (Table 2), which is not predicted by the model (Figure 3). This increase of measured δ¹³C values might be related to a change in dietary habit of the volunteer which was not recorded in our study. Another possible reason for the increase of measured δ¹³C value is that the volunteer took a series of short trips to southern Ontario during that period (Niagara Falls, Toronto, and Killbear Provincial Park in August, and Sudbury in October 2010) where the δ¹³C values can be much higher (−17.1‰) than Ottawa (−18.9 to −18.6‰) (Table 1). The predicted profile did not show an increasing δ¹³C value associated with those trips because they were too short to cause a significant shift in the predicted profile. Adoption of a more accurate hair growth model like the reaction progress model (Cerling et al., 2007) might improve the prediction and reduce such discrepancy between the measured and the predicted profiles.

We suggest that while δ¹³C values in hair likely respond to dietary preferences, they can provide further geolocation constraints for modern individuals. For example, the δ¹³C and δ¹⁵N values of human hair were used to identify the recent migrants in a rural community in SW England (Bol et al., 2007), and could also be used, along with collagen samples, to identify the migrants from Ireland to 19th century London during the Great Irish Famine (1845–1852) (Beaumont et al., 2013). As the isotopic δ¹³C baseline from different countries and regions become well established (Hülsemann et al., 2015), we suggest that δ¹³C value could be added to provide a qualitative assessment of mobility and constrain δ²H values. Similar applications could be envisioned for archeological applications as far as the δ¹³C variations in food systems of different regions at the studied time are established. Early studies have shown that the ancient mummies from different locations (Egypt, Chile, and Oetztaler Alps) have different δ¹³C and δ¹⁵N values in hair likely due to dietary differences of C₃ and C₄ food as a carbon source, and marine food as a nitrogen source (Macko et al., 1999a,b). A later study demonstrated the use of δ¹³C and δ¹⁵N values in archeological human hairs from Northern Chile to reconstruct the travel history of the individual moving between the coast and highlands (Knudson et al., 2012). Paleomobility studies would thus be further benefited from the δ¹³C hair analysis if the paleo-isoscape of food sources (e.g., δ¹³C variation due to C₃/C₄ proportion) could be established like the modern isoscapes (West et al., 2010).

The predicted δ¹⁵N profile has little fluctuation (Figure 4) despite the long-distance travel of the volunteer. The majority of the measured δ¹⁵N values in the hair profile of the individual falls within a tight range of 9.5–10.25‰ (Figure 4). These values are in the upper range of other Canadians in the FAD (Chartrand and St-Jean, 2015; Bataille et al., 2020a), suggesting, as δ¹³C values do, that our volunteer had a specific diet, possibly high intake of animal or marine-based protein (O’Connell and Hedges, 1999; Petzke et al., 2005b, 2010). Most of the values are within the predicted range except for a few outliers. The segment A1 has a δ¹⁵N value of 10.5‰ (Table 2), slightly higher than the upper limit of the predicted range (Figure 4), which might reflect a change of nutritional or metabolic status of the individual (Mekota et al., 2006; Petzke et al., 2010). We can also see several outliers in the hair segments A22–A29, with the highest δ¹⁵N value of 11.9‰ in the segment A27 (Table 2 and Figure 4).

Studies have shown that δ¹⁵N value in hair can be affected by both diet and nutritional stress (O’Connell and Hedges, 1999; Fuller et al., 2005; Mekota et al., 2006; Petzke et al., 2005a, 2010). Therefore, the drastic change of δ¹⁵N value measured in the segment A27 and A29 of the volunteer’s hair might be related to a severe illness reported by the volunteer during and post Cuba trip. This disease could have led to metabolic stress causing a change of nutritional stress or diet, as reported by the volunteer. However, the timing of the δ¹⁵N excursions in the measured profile does not match exactly the timing of the Cuba trip: one of the δ¹⁵N excursion (A29) occurs a month before the Cuba trip (Figure 4). As discussed before, it is possible that because these isotopic signals are located in the distal portion of the hair, their timing are more uncertain. We also notice some rapid δ¹⁵N fluctuations (8.8–9.7‰) in the hair segments A33–A37 (Table 2 and Figure 4). These variations do not match any travel and are likely influence by physiological or dietary factors.

In general, δ¹⁵N value in hair may not be a good tool to identify geographical movements because of its geographical invariability. δ¹⁵N values mainly depend on the protein source in the individual’s diet. As evidenced by previous studies, δ¹⁵N value in hair can be potentially a powerful tool to identify the change of diet, nutritional stress, or metabolic status of the subject (Fuller et al., 2004, 2005; Mekota et al., 2006; Huelsemann et al., 2009; Petzke et al., 2010).

Previous studies demonstrated that the ⁸⁷Sr/⁸⁶Sr ratios in hair are strongly related to the ⁸⁷Sr/⁸⁶Sr ratio of local tap water (Tipple et al., 2018, 2019; Ammer et al., 2020). If the hair of the volunteer recorded the ⁸⁷Sr/⁸⁶Sr ratios of tap water at traveled locations by drinking the tap water and incorporating its ⁸⁷Sr/⁸⁶Sr ratio into hair through the bloodstream like the light stable isotope systems, we would expect large ⁸⁷Sr/⁸⁶Sr variations throughout the predicted profile shown in Figure 5. However, the measured ⁸⁷Sr/⁸⁶Sr ratios display a low range of ⁸⁷Sr/⁸⁶Sr variations remaining around 0.7115 (0.7110–0.7119) which is much smaller than the expected range from the predicted profile (0.7100–0.7143) (Figure 5). A possible explanation for this discrepancy is that Sr in hair has a much long turnover time than hydrogen isotopes used to calibrate the hair age model (Cohn et al., 1962; Kitchings et al., 1976). A more likely reason is that hair incorporates Sr from tap water through surface contamination during bathing or showering instead of drinking (Tipple et al., 2018; Hu et al., 2020). As shown in previous studies, exogenous Sr can be rapidly incorporated into hair keratin after the eruption from the scalp (Chittleborough, 1980; Kempson and Lombi, 2011), and the ⁸⁷Sr/⁸⁶Sr ratios in human keratinous tissues can rapidly inherit the ⁸⁷Sr/⁸⁶Sr of bathing or showering water (Mancuso and Ehleringer, 2018; Tipple et al., 2018; Hu et al., 2020). Therefore, we would expect the ⁸⁷Sr/⁸⁶Sr ratios in the hair of the volunteer to incorporate exogenous Sr from the last location where the volunteer resided and/or from the location where the volunteer resided the longest (Hu et al., 2018, 2020; Tipple et al., 2018, 2019). Indeed, the ⁸⁷Sr/⁸⁶Sr ratios in the hair of the volunteer (0.7110–0.7119) fall close to the ⁸⁷Sr/⁸⁶Sr ratios of tap water in Ottawa, the last city where the volunteer resided (0.7126–0.7128) (Table 1). The slight difference between measured ⁸⁷Sr/⁸⁶Sr ratios in hair and that of Ottawa tap water may reflect seasonal variations of the ⁸⁷Sr/⁸⁶Sr ratio in tap water originated from the Ottawa river, with ⁸⁷Sr/⁸⁶Sr ratios of ∼0.711 in Spring and ∼0.712 in Fall from Telmer (1998). This hypothesis of exogenous incorporation is consistent with the longitudinal increase of Sr concentration ([Sr]) in the hair profile of the volunteer (Figure 6A) and the convergence of the hair ⁸⁷Sr/⁸⁶Sr ratios toward an end-member of ∼0.7112 (Figure 6B), equal to Ottawa river ⁸⁷Sr/⁸⁶Sr ratio in Spring (Telmer, 1998).

However, we also notice that the measured ⁸⁷Sr/⁸⁶Sr ratios in the two first hair segments (B1 and B2, Table 3) were quite lower (0.7110 and 0.7112) than the rest of the hair profile and sit on the bottom-right corner of Figure 6B, suggesting another source of Sr for segment B1 and B2. Both hair segments are the most recently grown, and as evidenced by their lower Sr concentrations (Figure 6A) they have incorporated less exogeneous Sr through bathing or showering. Young hair tends to be less damaged, thus limiting water diffusion during bathing. Consequently, the most proximal portion of the hair may preserve, in part, endogenous ⁸⁷Sr/⁸⁶Sr ratio. This is supported by the distinct isotopic ratios of these segments relative to the rest of the hair where Sr is overprinted or exchanged with exogenous Sr from Ottawa tap water. This endogenous Sr must come primarily from dietary sources and not from tap water and is incorporated into hair through blood (Lewis et al., 2017). The volunteer had a typical supermarket diet and most of the food items consumed while residing in Ottawa came from a variety of regional, national, and international sources with variable ⁸⁷Sr/⁸⁶Sr ratios. This complex food sourcing makes it challenging to characterize the integrated dietary ⁸⁷Sr/⁸⁶Sr ratio. However, it is likely that the ⁸⁷Sr/⁸⁶Sr ratios of the food consumed by the volunteer in Ottawa differ from that of the local tap water explaining the isotopic difference between the young and older hair segments.

Although studies have been trying to use ⁸⁷Sr/⁸⁶Sr ratios in hair for forensic and archeologic applications (Font et al., 2012; Frei et al., 2015; Vautour et al., 2015), it has been shown that hair can exchange Sr with the ambient water, either in modern environment during bathing or showering (Tipple et al., 2018; Hu et al., 2020), or in post-mortem environment (Gordon et al., 2019). Our result generally agrees that the Sr in modern human hair above 3 cm from scalp is mostly replaced by the exogenous signal. However, the initial 3 cm of hair (segments B1 and B2) in this study contained different ⁸⁷Sr/⁸⁶Sr ratios from the rest of the hair, indicating the preservation of the most recent endogenous ⁸⁷Sr/⁸⁶Sr ratios. Therefore, the proximal end (root) of hair could preserve endogenous Sr signal which can be useful for provenance applications and archeological applications allowing in particular to verify if an individual recently moved or if human remains were displaced.