Introduction

The stable oxygen (18O/16O) isotope values (δ18O) of human tissues can be utilized to identify and reconstruct an individual’s region-of-residence or origin as δ18O of human tissue relates largely to the δ18O of an individual’s drinking water, which varies with geography1. Numerous studies have applied the δ18O of human tissues to determine the likely origin of modern, historic, and prehistoric humans2,3,4,5,6,7,8,9. However, predictions of origin using δ18O can be geographically broad10,11. Similar to δ18O, variations in the strontium (Sr) isotope ratios (87Sr/86Sr) relate to geography and the analysis of 87Sr/86Sr in human tissues has been utilized to address many historical6,12,13,14 and prehistorical human provenance questions14,15,16,17,18. With 87Sr/86Sr, the age and geology of the underlying bedrock control the 87Sr/86Sr in soils, waters, plants, and animals19,20,21, which are incorporated into human tissues22. Thus, combining the independent δ18O and 87Sr/86Sr systems may allow for greater resolution in human provenance predictions.

Measurements of both δ18O and 87Sr/86Sr in tooth enamel have been successfully applied in several modern23,24 and prehistorical applications25 to identify the probable locations where an individual spent their early life. The measurement of δ18O and 87Sr/86Sr in hydroxylapatite tissues from humans is straightforward with tooth enamel having an average of 544 μg g−1 Sr26. Where tooth enamel represents early periods in an individual’s life, the isotope values of hair keratin represent much more recent intervals of time. Thus, the combined measurement of δ18O and 87Sr/86Sr in hair is emerging as a very attractive approach to understand an individual’s very recent life history27. This combination of isotopes has many potential applications, especially recognizing known variations in tap waters that exist within and among metropolitan regions28,29,30. However, while the δ18O systematics in keratin is relatively well understood and commonly applied1, the Sr isotopic system in human hair is not22.

The various means of Sr incorporation into internal human tissues are ingestion of food and beverages in addition to inhalation of aerosols and particles31,32; however, Sr in hair – an external tissue – represents a complex mixture of both these endogenous and exogenous sources of Sr33. Endogenous sources of Sr to hair originate from the body’s Sr pools within blood and bones31,32, while exogenous sources represent external environmental influences from aerosols, particulates, and environmental waters27,33,34. The overwhelming importance of exogenous contributions is evidenced by elevated Sr concentrations in the cortex of the hair35,36 as well as continual increases in Sr concentration from the proximal to distal terminus of the hair27,37. Chemical separation and isolation of the exterior and interior Sr pools demonstrated that these inputs may impart unique 87Sr/86Sr ratios and bias the overall 87Sr/86Sr value34. Distinguishing the impact of these two potentially conflicting Sr sources on the overall 87Sr/86Sr value is required to fully utilize 87Sr/86Sr of hair for the reconstruction of region-of-residence and remains a fundamental question in the field.

To begin to distinguish the importance of endogenous and exogenous Sr to hair and create a better understanding of how these sources of Sr relate to the spatial signal encoded within 87Sr/86Sr of hair, we collected and analyzed the 87Sr/86Sr and δ18O of human hair from three locations within a single city (Salt Lake City, Utah, USA). In addition to 87Sr/86Sr and δ18O, the sex of the hair donors were also gathered, as well as, the 87Sr/86Sr and δ18O of tap waters collected from six sites near hair collection locations. Our aim was to explore the relationships of 87Sr/86Sr and δ18O of hair at the city-scale and to begin to develop a better understanding of the geospatial signal encoded in the 87Sr/86Sr of hair.

Results

Oxygen isotope values of hair

Sixty-seven hair samples were obtained from students from three public schools (Fig. 1). We analyzed thirty-one hair samples for δ18O and found hair samples ranged in δ18O from 6.7 to 12.8‰ with a mean δ18O of 9.9 ± 0.8‰ (1σ, n = 31). The mean δ18O of hair samples from females was 9.8 ± 1.4‰ (1σ, n = 28), while the mean δ18O of hair samples from males was 10.3 ± 0.4‰ (1σ, n = 3) (Fig. 2). While our dataset was weighed toward female participants, the δ18O of hair from females was not significantly different than δ18O of hair from males (Student’s t-test, p = 0.5370). When grouped by collection location, we found the mean δ18O of hair samples from individuals at Schools X, Y, and Z were 9.7 ± 1.8‰ (1σ, n = 10), 10.0 ± 1.5‰ (1σ, n = 7), and 10.0 ± 0.9‰ (1σ, n = 14), respectively (Fig. 2). Here, the δ18O of hair from individual schools were not significantly different from one another (one-way ANOVA, p = 0.8438).

Figure 1
figure 1

Location map showing the collection locations of hair and tap water. Salt Lake City is highlighted. School locations are shown with flagged symbols and tap water collection locations are shown with points. Point shading indicates groupings with group A sites shown in grey, group B sites shown in black, and the group C site shown in white. Figure produced using ArcGIS 10.4 (ESRI; Redlands, CA, USA).

Figure 2
figure 2

Box and whisker plot showing the δ18O values of hair from female and male donors and hair data isolated by collection location.

Oxygen isotope values of water

Water samples were collected from taps in six public buildings four times a year for four years (Fig. 1). The mean δ18O of collected tap water was −16.3 ± 0.2‰ (1σ, n = 66) with a range from −15.3 to −16.8‰ (Table 1). Following Ehleringer et al.1, the δ18O of drinking water predicted from measured δ18O of hair ranged between −12.5‰ to −17.0‰ with a mean δ18O of −14.8 ± 3.7‰ (1σ, n = 31) (Fig. 3). Here, we found the measured δ18O of tap water was not significantly different than δ18O of drinking water predicted from hair (Welch’s t-test, p = 0.0707).

Table 1 Oxygen and strontium isotope ratios of tap waters from sites within Salt Lake City, Utah.
Figure 3
figure 3

Histogram showing the measured δ18O values of tap water at collection locations and the predicted δ18O values of drinking water from the measured δ18O values of hair.

Strontium abundance of hair

Sixty-seven hair samples were analyzed for Sr elemental abundance [Sr]. We found [Sr] ranged from 0.08–27.5 μg g−1, while the mean [Sr] of hair samples from female and male participants was 7.8 ± 5.0 μg g−1 (1σ, n = 60) and 3.7 ± 3.5 μg g−1 (1σ, n = 7), respectively. [Sr] of hair samples from female participants was significantly higher than [Sr] of hair from males (Welch’s t-test, p = 0.0003). The mean [Sr] of hair samples from students at Schools X, Y, and Z were 7.4 ± 4.7 μg g−1 (1σ, n = 11), 8.1 ± 4.8 μg g−1 (1σ, n = 29), and 6.7 ± 4.0 μg g−1 (1σ, n = 27), respectively. We did not find any significant differences in [Sr] of hair samples from collection locations (one-way ANOVA, p = 0.5510).

Strontium isotope ratios of hair

Forty-one hair samples were measured for 87Sr/86Sr and ranged from 0.70910 to 0.71509 (Fig. 4). The mean 87Sr/86Sr of hair samples from female and male participants was 0.71215 ± 0.00193 (1σ, n = 37) and 0.71150 ± 0.00139 (1σ, n = 4), respectively. Here, we observed the 87Sr/86Sr of hair from female and male students were not significantly different (Student’s t-test, p = 0.5210). The mean 87Sr/86Sr of hair samples from students at Schools X, Y, and Z were 0.70965 ± 0.00028 (1σ, n = 3), 0.71342 ± 0.00152 (1σ, n = 18), and 0.71125 ± 0.00142 (1σ, n = 20), respectively. These differences in 87Sr/86Sr of hair samples at individual schools were significantly different from one another (one-way ANOVA, p < 0.0001).

Figure 4
figure 4

Dot plot showing the 87Sr/86Sr of hair from female and male donors, 87Sr/86Sr of hair from each collection location, and 87Sr/86Sr of tap water from the three tap water Site Groupings. The line represents the mean 87Sr/86Sr ratio.

Strontium isotope ratios of water

Sixty-one tap water samples from six collection sites were analyzed for 87Sr/86Sr; they ranged from 0.70906 to 0.71399 (Table 1). Tap water samples were grouped based on proximity to hair sampling locations into Site Groups. The mean 87Sr/86Sr of tap water samples from Sites 1, 2, and 3 (Site Group A), Sites 4 and 5 (Site Group B), and Site 6 (Site Group C) were 0.71028 ± 0.00028 (1σ, n = 29), 0.71192 ± 0.00181 (1σ, n = 20), and 0.71014 ± 0.00060 (1σ, n = 12), respectively (Fig. 4). We found the mean 87Sr/86Sr of Site Groups were significantly different from one another (one-way ANOVA, p < 0.0001).

Discussion

Oxygen isotope analysis of human hair is becoming a relatively common tool in determining the region-of-origin of modern humans3,5,38. Here, we observed the δ18O of hair from individuals to be indistinguishable at the sites sampled in Salt Lake City, Utah, likely because the origins of urban tap waters share a common montane origin29,39. These data are consistent with previous reports of the δ18O of hair from humans and animals living in or near this region1,27. Previously, Ehleringer, et al.1 found that the average δ18O of human hair from individuals living in Utah was 9.7 ± 1.5‰ (1σ, n = 12). When compared to our dataset, we found no significant difference between our data and these previous data (Student’s t-test, p = 0.6478). The δ18O of hair has been linked to the δ18O of drinking water in humans40,41,42 and to assess the relationship between individuals and their drinking water, we translated the δ18O of hair to δ18O of drinking water following Ehleringer, et al.1. Consistent with expectations, the δ18O of drinking water predicted from measured δ18O of hair was not significantly different than measured δ18O of tap water. This is in agreement with previous studies that indicate the δ18O of an individual’s hair largely relates to the δ18O of their drinking water, and thus the environment and community in that they reside3,43,44,45. Further, it is not unexpected that there were no differences in the δ18O of hair from these individuals (Figs 2 and 3), given the relatively small range of δ18O of tap water within the regions where these individuals were residing.

Studies of the δ18O of tap water have described spatial and temporal dynamics within single metropolitan areas, including in those around Salt Lake City28,29,30. These dynamics may confound the application of δ18O of human hair as a geographic predictor, particularly within regions that utilize transported water46,47. Within the Salt Lake metropolitan area, Jameel et al.29 showed that different regions within the study area utilized different management strategies to supply culinary water, resulting in subtle, but unique δ18O values of tap water across the area. During the three studied seasons, Jameel et al.29 found the δ18O of tap water delivered to Salt Lake City had a relatively small range of isotope values (i.e., −16.5 to −14.0‰), consistent with our findings. This was not unexpected as all waters delivered to Salt Lake City ultimately originated from cold season, high elevation precipitation29,39. However, the utility supplying Salt Lake City utilized both surficial and groundwater sources depending on season and supply level48. While surficial and groundwater sources would have similar δ18O values given their common source, other geochemical tracers likely vary due to water-rock interactions and other processing. If there are underlying geographic controls of the supply of groundwater to an urban region (e.g., proximity to groundwater wells, etc.), then additional geochemical or isotopic tracers may be useful to uncover intra-city scale variations in water supply and possibly chemical signatures within the individuals residing within these regions.

Strontium in human tissue is generally considered a conservative tracer of geographic origin22, however, in human hair the [Sr] and 87Sr/86Sr are controlled by both diet and environmental contamination33,34. The endogenous and exogenous Sr pools in human hair may have unique or quite similar 87Sr/86Sr depending on the Sr source and either of these pools may or may not preferentially influence the final Sr isotope value34,49. Multiple methods have been developed to potentially isolate endogenous and exogenous Sr for dietary and provenancing studies in humans, yet no consensus has been reached on the preferred methodology34,36,50,51,52,53,54,55,56. Regardless of these known uncertainties, the 87Sr/86Sr of human and animal hair have been shown to trace movement histories27,37,57. However, it remains unclear if the geographic signal within 87Sr/86Sr in hair relates to the endogenous or exogenous Sr signature, or some combination of both.

Here, we found the range of [Sr] in hair was in good agreement with previous studies of hair from students58,59. We observed [Sr] of hair from female participants was significantly higher than [Sr] of hair from males. This too was consistent with other studies that have found that [Sr] in hair a key element to differentiate females and males58. In previous studies, [Sr] has been linked to sex, with females most often having greater [Sr] than males58,59,60,61,62,63. While the reason(s) behind this pattern remains unclear, it has been suggested that differences in physiology, biochemistry, diet, and/or activity levels may be the causes of these patterns31,33,64. [Sr] is known to increase along the length of hair from proximal to distal portions due to the incorporation of exogenous strontium on and/or into the hair cuticle27,35,65. Given that our dataset is weighed towards female participants and females are more likely to wear their hair longer than males, these differences may be related to female participants’ hair incorporating exogenous strontium for extended periods of time, and thereby increasing the [Sr].

While we observed a difference between [Sr] of hair from females and males, we did not find any significant differences in [Sr] of hair samples from collection locations within Salt Lake City. Given the small geographic area from that the individuals resided and our dataset bias towards female participants, it may not be unexpected that there were no differences in [Sr] of hair samples between the three collection locations.

In contrast to [Sr] and δ18O values of hair, we observed an unexpectedly wide range in the 87Sr/86Sr of hair. Here, we found the 87Sr/86Sr of hair from individuals residing in Salt Lake City varied by more than 0.005. This variation was larger than any previously reported 87Sr/86Sr from hair specimens from individuals residing in the same city or geographic region34,37. In a previous study using hair specimens collected from a single salon in Taylorsville, Utah located less than 10 km south of Salt Lake City (Fig. 1), Tipple and colleagues34 found the 87Sr/86Sr of hair ranged from 0.70909 to 0.71469 with an average of 0.71203 ± 0.00140 (1σ, n = 22)34. Our data were not statistically different than these previous data from Taylorsville, Utah (Student’s t-test; p = 0.9050). Given the constrained geographic area with identical climate and underlying geology, these variations point towards an additional control for the wide range of 87Sr/86Sr of hair observed within this single metropolitan area.

There remains no consensus regarding the dominant geospatial control on the 87Sr/86Sr of hair27,34,37,57. Several studies have assessed the relationship between the 87Sr/86Sr of modern hair and residence with some studies suggesting that the endogenous signal to hair is most important for geospatial relationships57, while others found compelling relationships with exogenous sources, such as water, atmospheric dust, etc.27,37. Unraveling the significant contributions of Sr to hair for geospatial reconstructions is becoming more important in both modern cold case66 and archeological67,68 contexts.

We found the overall range of 87Sr/86Sr of hair and tap water was similar at the locations sampled in Salt Lake City (Fig. 4). Yet, the 87Sr/86Sr of hair from participants varied systematically by collection location within the city (Fig. 4). We noted that participants from School Y had higher 87Sr/86Sr ratios and a much larger range than the other two collection locations. Previously, Vautour and colleagues37 noted that four individuals that traveled from Paris or St. Benoit sur Loire, France to Montreal, Canada showed a transition in 87Sr/86Sr along a length of hair corresponding to the change in location. They found that the 87Sr/86Sr of hair from these individuals converged towards an average value with some slight differences between the four individuals and used this finding as evidence that the Sr in hair is not only controlled by endogenous inputs, but exogenous contributions as well. Furthermore, they noted that the 87Sr/86Sr of hair was similar, but distinct from the 87Sr/86Sr of local tap water and argued that this was due to reworking and buffering of the Sr isotopic signal from Sr incorporated within bone37. Similarly, Font and associates57 argued that the 87Sr/86Sr of hair from two individuals that traveled from Kanpur, India to Amsterdam, The Netherlands slowly reached isotopic equilibrium with the new location due to recycling and incorporation of Sr from bone32. While both dietary Sr and remobilization of bone Sr likely contribute to overall 87Sr/86Sr of hair, the impact of the 87Sr/86Sr of water on hair may not have been fully appreciated in these previous studies. Here, we analyzed sixty-one water samples in a single municipality and found an astonishingly wide range of 87Sr/86Sr ratios. These previous studies analyzed five or fewer tap waters from each of the locations that the individuals resided, and thus, may have not captured all the possible 87Sr/86Sr variability that could exist. Here, we found spatial relationships between the 87Sr/86Sr of hair and tap water from the neighborhoods study participants resided in, suggesting interactions with local tap waters may contribute to the 87Sr/86Sr of hair. While the specific mechanism Sr uptake remains unclear, it may take place through imbibing local tap waters or when the hair is wetted during bathing or showering. While this study cannot specifically distinguish between the incorporation of endogenous or exogenous Sr through interactions with water, previous research has shown that significant amounts of Sr can be added to hair after it exits the scalp and is exposed to the environment33. Thus, we hypothesize that interactions with bathing water may be an important contributor of Sr to hair and that small spatial variations in the 87Sr/86Sr of water are possibly recorded in the 87Sr/86Sr of hair. Geospatial modeling of bioavailable 87Sr/86Sr for human provenancing remains nascent19,20,21,69,70,71 and our results suggest that these modeling exercises should consider additional factors or layers for applications towards modern humans72, including the 87Sr/86Sr of tap waters and the processes by which communities acquire and transport their water resources. Nonetheless, additional controlled dietary studies that consider the 87Sr/86Sr of food, beverages, and external Sr sources will be required to test this hypothesis and establish a mechanism for Sr incorporation from bathing water.

Methods

Ethics statement

The Institutional Review Board (IRB) of the University of Utah approved this research program under protocol number [00035524]. Specifically, all sampling and analytical methods used were in accordance with these IRB regulations. Informed consent was obtained from all subjects or from their legal guardians in accordance with and maintained under IRB regulations.

Hair and water samples

Sixty-seven hair samples were obtained from students from three public schools following Valenzuela, et al.73 (Fig. 1). Schools X and Y are public schools that pull students from defined areas near the school, while School Z is a public magnet school that has students from throughout the public school district. All study participants lived in Salt Lake City, Utah. Participants were assigned a randomized sample identifier and all personal and identifying information was stored according to IRB. Hair samples were placed in paper coin envelopes and labeled with sample identifier. Hair samples were returned to laboratory and stored. In addition to hair samples, participants provided their sex, along with any specific dietary or health information.

Water samples were collected from taps in six public buildings four times a year for four years (summer, 2012 to summer, 2015). Water collection sites were located within a few blocks of hair collection locations (Fig. 1). The plumbing was flushed prior to sample collection by running the cold-water tap for 10 s. Water samples were collected in 4-ml glass vials and the cap seal was with Parafilm®. Water samples were returned to laboratory and stored in 4 °C refrigerator prior to analysis.

Hair cleaning

Hair samples were cleaned with 2:1 chloroform:methanol (v/v). Samples were enclosed in filter paper and completely submerged in the solvents. Samples were agitated using a shaker plate and after 5 min the supernatant was decanted and discarded. This process was repeated a total of three times. Chloroform (OmniSolv®, EMD; Darmstadt, Germany) and methanol (OmniSolv®, EMD; Darmstadt, Germany) used in hair cleaning were HPLC grade. Cleaned hair samples were allowed to dry at room temperature for 72 hrs within a laminar flow hood. After dried, cleaned hair samples were ground in a ball mill and stored in ashed glass vials until further analysis-specific preparation procedures.

Oxygen isotope analysis of hair

Hair samples equilibrated with ambient laboratory atmospheric water vapor for 48 hrs alongside in-house keratin reference materials (RMs) following Chau, et al.27. Hair samples (~150 mg), as well as RMs, were weighed in 3.5 × 5 mm silver capsules (Costech Analytical Technologies, Inc.; Valencia, California, USA). Weighed samples and RMs were stored under vacuum for 7 days before being analyzed. The stable oxygen isotope values of keratin were determined with a continuous flow isotope ratio mass spectrometer (IRMS) (MAT 253, Thermo Finnigan; Bremen, Germany), housed at IsoForensics, Inc. in Salt Lake City, Utah. Samples were introduced to the IRMS via a zero-blank autosampler attached to a high temperature conversion elemental analyzer (TC/EA, ThermoFinnigan; Breman, Germany). Samples were analyzed alongside sets of natural keratin RM for primary quality assurance and secondary quality control (QA/QC). Keratin RMs were previously calibrated to Vienna Standard Mean Ocean Water (VSMOW)-Standard Light Antarctic Precipitation (SLAP) international isotope scale. Sets included two primary QA RMs for slope/intercept normalization and a secondary QC RMs to insure suitable calibration. Primary keratin QA RMs were DS and ORX (δ18O = 6.02‰ and 25.09‰, respectively), while POW was used for QC (long-term mean δ18O = 12.44‰, 1σ = 0.54‰, n = 335). A minimum of four sets of RMs was analyzed along side samples. All δ18O values of hair are expressed relative to VSMOW. All light stable isotope values are reported in parts per thousand (‰) and in delta notation:

$$\delta =[({R}_{samp}/{R}_{std})-1],$$

where R represents the 18O/16O abundance ratio, and R samp and R std are the ratios in the sample and standard, respectively.

Oxygen isotope analysis of water

Water samples were analyzed using cavity ring-down water isotope spectroscopy (L1102-I, Picarro; Sunnyvale, California). Each sample was analyzed four times (four consecutive replicate injections) alongside a set of three liquid water laboratory RMs, (ZE = −0.2‰, EV = −10.2‰, DI = −16.5‰ for δ18O values) that had previously been calibrated to the VSMOW-SLAP international isotope scale. Two QA RMs were used for data normalization (ZE and DI) and a QC RM for quality control (EV). Analytical precision of the QA RM was ± 0.3‰ for δ18O values. All δ18O values of water are expressed relative to VSMOW.

Digestion for Sr abundance and isotope analysis

Hair samples were digested using an Ethos EZ® microwave digestion system (Milestone, Inc.; Shelton, Connecticut, USA). Approximately 50 mg of hair was weighed into a Teflon® digestion microvessel. Two milliliters of concentrated ultrapure concentrated HNO3 (Aristar® ULTRA, BDH Chemical; Darmstadt, Germany) were added to the microvessel containing the hair and the microvessel was then sealed, submerged in 10 ml of milli-Q water and 50 ml H2O2 (30% v/v), and then placed within an outer vessel. The outer vessel was then placed in the digester carousel. Two certified reference materials (TORT-2, National Research Council, Ottawa, Canada; Human Hair No. 13, National Institute for Environmental Studies, Tsukuba, Japan) and a method blank of reagents were digested along with the hair samples. The microwave program used for hair digestion was 13.3 °C/min ramp to 200 °C, followed by an isothermal at 200 °C for 15 min with a 60 min cool down to room temperature. The microwave was operated at full power (1500 W) for all heating cycles. Once cooled to room temperature, the hair digests were transferred to acid-leached 2-ml snap-cap centrifuge tubes. A 100-ml aliquot of the primary hair digest was transferred to a 15-ml tube and the volume was brought to 10 ml with ultrapure water. The ultrapure water used for sample cleaning and acid dilutions was from a Milli-Q Academic A10® system (EMD Millipore; Billerica, Massachusetts, USA) with a resistivity >18 MΩ. A standard solution containing 10 ppb In was added to each sample as an internal concentration standard.

Strontium abundance analysis

All strontium elemental abundances were measured via inductively coupled plasma quadrupole-mass spectrometry (ICP-MS) (Agilent 7500ce, Agilent Technologies; Santa Clara, California, USA) at the Department of Geology & Geophysics at the University of Utah, Salt Lake City, Utah. A double-pass spray chamber with perfluoroalcoxy fluorocarbon (PFA) nebulizer (0.1 mL/min), a quartz torch, and nickel cones were used. A calibration solution containing Sr was prepared gravimetrically using a single-element standard (Inorganic Ventures, Inc.; Christiansburg, Virginia, USA). Standard reference solution T-205 (USGS; Reston, Virginia, USA) was measured as an external calibration standard at least five times within each analytical run. The long-term reproducibility for T-205 and differences relative to the accepted values indicated that the Sr concentrations were accurate within 10%. TORT-2 has a certified Sr concentration of 45.2 ± 1.9 mg g−1 (1σ) and the measured Sr concentration of TORT-2 was 41.2 ± 4.4 mg g−1 (1σ, n = 6).

Strontium isotope analysis

All strontium isotope measurements were made using a Neptune Plus multi-collector ICP-MS (Thermo Fisher Scientific; Bremen, Germany) housed in the Department of Geology & Geophysics at the University of Utah, Salt Lake City, Utah, USA. Digests were introduced using an online Sr purification method following Tipple, et al.34 for hair samples and Chesson, et al.74 for water samples. This online system automates the purification of Sr by utilizing a peristaltic pump, a pair of 6-way valves, an in-line separation column, and a SC-2 DX autosampler with a FAST2 valve block (Elemental Scientific; Omaha, Nebraska, USA). The in-line separation column was packed with crown ether Sr resin (Eichrom Technologies; Lisle, Illinois, USA). Variable speed settings on the peristaltic pump allowed samples to be rapidly loaded into the purification column where Sr was trapped while all other elements were rinsed away; the column flow was then reversed and purified Sr was eluted into the spray chamber. A timing solution containing 66 ppb Sr was analyzed daily to insure proper chromatography and to assess column chemistry and efficiency. The instrument was operated at an RF power of 1200 W with nickel sampling and skimmer cones (1.1 mm and 0.8 mm apertures, respectively) and was optimized daily for signal intensity and stability. Cool, auxiliary and sample gas flow rates were 16 L/min, 0.85 L/min, and 0.91 L/min, respectively. The instrument was tuned for sensitivity daily with a solution containing 20 ppb Sr. For 87Sr/86Sr analysis, a static multi-collector routine was used that consisted of 1 block of 170 cycles with an integration time of 1.032 sec per cycle for an individual analysis. Each analysis was followed by a blank to monitor the efficiency of the crown ether Sr resin column. Sr isotope ratios of samples and references were blank- and interference-corrected and then normalized for instrumental mass discrimination using a defined 86Sr/88Sr of 0.1194.

Samples were analyzed for 87Sr/86Sr alongside sets of reference materials and blanks. Reference sets included one primary QA reference for normalization and a keratin QC references to insure measurement reproducibility. The primary reference material was SRM® 987 (0.71034 ± 0.00026 [95% CI]; National Institute of Standards and Technology; Gaithersburg, Maryland, USA), while Human Hair No. 13 was used for QC. Samples and SRM® 987 were analyzed at 5:1. The measured 87Sr/86Sr of SRM 987 and Human Hair No. 13 was 0.71030 ± 0.00004 (1σ, n = 92) and 0.70827 ± 0.00004 (1σ, n = 6), respectively.

Statistical analysis

Statistical analysis was completed using JMP Pro 13® (SAS Institute Inc.; Cary, NC, USA). Normality of the distributions were tested with the Shapiro-Wilk test. If the distributions were normal, then the Student’s t-test was used to compare means at α = 0.005. If the distributions were not normally distributed, then the Welch’s t-test was used at differences of α = 0.005. One-way ANOVA test was used to assess differences between groups at α = 0.005.

Mapping

Mapping and original figure creation was conducted using ArcGIS 10.4 (ESRI; Redlands, CA, USA).

Tap water sources

The culinary water delivered to Salt Lake City originates from mountain streams, surface water reservoirs, and groundwater wells and springs48 with each source having a distinct 87Sr/86Sr value (Supplemental Data Table S1). Groundwater in the Salt Lake Valley has more a radiogenic 87Sr/86Sr than surface water sources (Supplemental Data Table S1). School Y is located near several groundwater wells with groundwater being most likely utilized in the autumn and winter (personal communication, M. Hubbard-Rice, Salt Lake City Public Utilities).