Articles | Volume 12, issue 4
Earth Syst. Sci. Data, 12, 3641–3652, 2020
Earth Syst. Sci. Data, 12, 3641–3652, 2020

Data description paper 23 Dec 2020

Data description paper | 23 Dec 2020

Bioavailable soil and rock strontium isotope data from Israel

Bioavailable soil and rock strontium isotope data from Israel
Ian Moffat1,2, Rachel Rudd1, Malte Willmes3,4, Graham Mortimer2, Les Kinsley2, Linda McMorrow2, Richard Armstrong2, Maxime Aubert5,6,2, and Rainer Grün6,2 Ian Moffat et al.
  • 1Archaeology, College of Humanities, Arts and Social Sciences, Flinders University,Bedford Park, 5042, Australia
  • 2Research School of Earth Sciences, The Australian National University, Canberra, 2600, Australia
  • 3University of California Santa Cruz, Institute of Marine Science, Santa Cruz, CA, 95064, USA
  • 4Southwest Fisheries Science Center, National Marine Fisheries Service, Santa Cruz, USA
  • 5Griffith Centre for Social and Cultural Research, Griffith University, Southport, 4222, Australia
  • 6Australian Research Centre for Human Evolution, Griffith University, Nathan, 4111, Australia

Correspondence: Ian Moffat (


Strontium isotope ratios (87Sr 86Sr) of biogenic material such as bones and teeth reflect the local sources of strontium ingested as food and drink during their formation. This has led to the use of strontium isotope ratios as a geochemical tracer in a wide range of fields including archaeology, ecology, food studies and forensic sciences. In order to utilise strontium as a geochemical tracer, baseline data of bioavailable 87Sr 86Sr in the region of interest are required, and a growing number of studies have developed reference maps for this purpose in various geographic regions, and over varying scales. This study presents a new data set of bioavailable strontium isotope ratios from rock and soil samples across Israel, as well as from sediment layers from seven key archaeological sites. This data set may be viewed and accessed both in an Open Science Framework repository (, Moffat et al., 2020) or via the IRHUM (Isotopic Reconstruction of Human Migration) database.

1 Introduction

Strontium (Sr) isotope geochemistry has applications in many fields of research including archaeology (Bentley, 2006; Slovak and Paytan, 2012), ecology (Barnett-Johnson et al., 2008; Hobson et al., 2010), food traceability (Voerkelius et al., 2010) and forensic sciences (Beard and Johnson, 2000). Strontium is widely distributed within geological and biological materials, and the strontium isotope ratios (87Sr 86Sr) of these materials reflect the sources of strontium in the environment during their formation (Dasch, 1969). The use of traditional isotope systems such as hydrogen and oxygen systems in provenance studies is limited by their broad gradients across the Earth's surface, whereas 87Sr 86Sr ratios show predictable geographic variability determined by lithology, with some limited temporal variability (Aggarwal et al., 2008; Willmes et al., 2018). To utilise the potential of strontium isotopes as a geochemical tracer, prior knowledge of baseline data in the region of interest is required. Bedrock strontium isotope ratios are a product of the age, mineralogy and origin of source material (Capo et al., 1998; Faure and Mensing, 2005). The amount of 87Sr within a geological unit typically increases over time through the radioactive decay of 87Rb to 87Sr, while the amounts of 84Sr, 86Sr and 88Sr remain constant (Capo et al., 1998). In general, older rocks such as granite in continental crust have higher 87Sr 86Sr ratios than younger rocks with a lower Rb  Sr ratio, such as oceanic basalt (Capo et al., 1998). The 87Sr 86Sr ratios of carbonates reflect the water from which they precipitated, either seawater for marine carbonates or a local water body for non-marine carbonates (Faure, 1986; Neat et al., 1979). Clastic sedimentary rocks typically have 87Sr 86Sr ratios which reflect the lithology of their source material, tempered by the relative ease with which minerals containing strontium are removed by weathering (McDermott and Hawkesworth, 1990). The 87Sr 86Sr ratios observed for geological units of varying compositions and ages are discussed further by Capo et al. (1998), Faure and Mensing (2005), Faure and Powell (1972), and others.

The strontium isotope composition of regolith is principally derived from local weathering of bedrock (Capo et al., 1998) but may be augmented by windblown marine, mineral and anthropogenic aerosols; different mineral weathering rates; transportation of regolith; and pore water chemistry (Frumkin and Stein, 2004; Goede et al., 1998). Bioavailable strontium is the component of strontium present in an environment which is available for incorporation into biological systems. Regolith, in combination with water, is the principal source of bioavailable strontium, although processes including precipitation, sea spray and fertiliser use may also affect the isotopic composition of bioavailable strontium (Bentley, 2006; Frei and Frei, 2011; Price et al., 2002; Slovak and Paytan, 2012).

1.1 Strontium isotope mapping

There are a range of approaches to mapping bioavailable strontium, as summarised below and in greater detail in Bentley (2006), Maurer et al. (2012), Slovak and Paytan (2012) and Bataille et al. (2020). Archaeological human or faunal samples are considered to be one of the best indicators of the local bioavailable range, although their geographic origin may be uncertain (Maurer et al., 2012; Price et al., 2002). Modern faunal samples from known localities may also be used, although non-local food sources and fertiliser use need to be considered (Bentley, 2006; Maurer et al., 2012). Soils, plants and water may also be measured for local bioavailable strontium isotope ratios (Evans et al., 2010; Maurer et al., 2012; Price et al., 2002). While soils and plants average strontium over a smaller area, ground water and surface waters may provide an estimate of bioavailable strontium over a wider area, depending on catchment size (Evans et al., 2010; Willmes et al., 2014). As discussed previously, the strontium isotope ratios of plants, soils and water may be affected by numerous factors including precipitation, sea spray and fertiliser use, which may vary over seasonal and annual timescales (Bentley, 2006; Hoogewerff et al., 2019; Price et al., 2002). Local bioavailable strontium isotope ratios may also be modelled from bedrock taking into account age, mineralogy and weathering rates (Bataille et al., 2012; Bataille and Bowen, 2012). In the case of soil and rock samples, an ammonium acetate or ammonium chloride digestion method, which removes cations from electrostatically bound pore water in addition to exchange sites on minerals and organic matter, can be used as a means of extracting the bioavailable component of these samples (Stewart et al., 1998). This approach has been found to extract between 0.1 % and 62 % of whole soil strontium (Chadwick et al., 2009).

At a country scale, strontium isoscapes have been developed for Mesoamerica (Hodell et al., 2004), the United Kingdom (Evans et al., 2010), France (Willmes et al., 2014, 2018) and Denmark (Frei and Frei, 2011) from a combination of bedrock, soil, plant and water samples. Isoscapes of bioavailable strontium have also been modelled for the contiguous USA (Bataille and Bowen, 2012), western Europe (Bataille et al., 2018; Hoogewerff et al., 2019) and the Caribbean region (Bataille et al., 2012). At a larger scale, Voerkelius et al. (2010) sampled mineral waters across western Europe to predict bedrock 87Sr 86Sr but found local processes to result in unexpected values.

1.2 Strontium isotope studies and mapping in Israel

Previous studies in Israel and across the Levant region have utilised strontium isotope analyses in a variety of applications in addition to undertaking small-scale mapping of bioavailable strontium. Herut et al. (1993) measured the 87Sr 86Sr ratios of rainwater samples across Israel and used these results in conjunction with the chemical composition of the sample to determine the sources of Sr and soluble salts. Shewan (2004) developed a map of bioavailable strontium throughout Israel using modern faunal bones and grass samples, for comparison to bones from archaeological sites in the region. The work of Shewan (2004) and Perry et al. (2008), who measured 87Sr 86Sr from archaeological faunal dental samples in western Jordan, was collated to create a bioavailable strontium map for Israel and Jordan, combining modern flora and fauna with material from archaeological sites (Perry et al., 2009). These studies highlight several distinct provinces in the region, with samples along the Jordan Valley distinguished from those in the eastern highlands of Israel and the Western Highlands in Jordan, as well as several smaller areas with varying isotopic signatures in the north of the study area (Perry et al., 2009). Rosenthal et al. (1989) report the strontium isotope composition of gastropod shells to reconstruct water sources in the Dead Sea Rift area of the Jordan Valley. Spiro et al. (2011) measured 87Sr 86Sr ratios of water and mollusc shells in the Hula catchment in the Upper Jordan Valley, adjacent to the Golan Heights, and found distinct water sources and aquifers. Stein et al. (1997) sampled water and sediments around the Dead Sea to determine water sources and the evolution of the Dead Sea and its precursor, Lake Lisan, identifying two distinct periods of lake evolution. Hartman and Richards (2014) sampled plants and invertebrates in northern Israel and the Golan Heights to produce a map of modern bioavailable strontium isotope ratios, and to investigate potential sources of variability, including inter-site variability and the influence of precipitation. Arnold et al. (2016) sampled plants in the vicinity of Tell es-Safi/Gath to create a local bioavailable strontium map as a baseline for use in interpreting the mobility of domestic animals from archaeological sites in the region. Moffat et al. (2012) undertook spatially resolved strontium isotope analysis using laser ablation techniques to demonstrate annual migration using a bovid tooth from the archaeological site of Skhul. These studies have highlighted the value and potential of bioavailable strontium mapping in this region, and the data presented in this study aim to add to this growing knowledge base.

2 Methods

2.1 Sample collection

Soil and rock samples were collected throughout Israel in September and October 2008. Sample locations were chosen opportunistically, with reference to a digital 1 : 200 000 geological map of Israel using the Old Israel Grid co-ordinate system (Sneh et al., 1998), to provide the greatest representation of the stratigraphic units present in Israel that were accessible via road. A single rock and soil sample was taken at each site, with no replicates. Soil samples were collected from the surface, at a single point, with no attempt made to sample multiple soil horizons where present. Rock samples were taken for the principal geological unit present at each site and a brief description of the lithology made in the field, which can be viewed in the data file. Sediments from stratigraphic layers in seven archaeological sites, Amud, the Atlit Railway Bridge site, Neve David, Qafzeh, Sefunim Cave, Skhul and Tabun, were also sampled for analysis via fieldwork or from archival collections. Where possible, sediments from archaeological sites were collected for multiple stratigraphic units. Sample locations and archaeological sites are illustrated in Fig. 1.

Figure 1Sampling locations overlain on a lithological map of Israel (adapted from Sneh and Rosensaft, 2014).

2.2 Analytical methods

2.2.1 Sample treatment

Rock and soil samples were heated to 60 C for a minimum of 24 h prior to sample preparation to comply with Australian quarantine procedures. After heating, rock samples were crushed to a powder using a hand piston. Rock and soil samples were passed through a 2 mm sieve, and 1 g aliquot of the sieved sample was leached by adding 2.5 mL of 1 M ammonium nitrate (NH4NO3), following DIN (German Industrial Standard) ISO (International Organization for Standardization) 197310 (2009). Samples were shaken for 24 h and subsequently centrifuged at 3000 rpm for 5–10 min. 1 mL of supernatant was extracted and evaporated until dry, before being dissolved in 2 mL of 2 M high-purity nitric acid (HNO3).

The concentration of strontium was determined by ICP-AES (inductively coupled plasma atomic emission spectrometry). Ion exchange chromatography was then used to isolate strontium from other elements, particularly 87Rb to prevent isobaric interference with 87Sr (Dickin, 2000), using two columns filled with Eichrom strontium-specific resin (prefilter resin and strontium specific resin). The strontium concentration determined by ICP-AES was used to determine the amount of sample added to the ion exchange column. Samples were diluted prior to MC-ICP-MS (multi-collector inductively coupled plasma mass spectrometer) analysis to provide sufficient additional samples to allow for reanalysis if necessary.

2.2.2 Neptune MC-ICP-MS analysis

A Neptune MC-ICP-MS was used to measure strontium isotope ratios in the Environmental Geochemistry and Geochronology Laboratory at the Research School of Earth Sciences, Australian National University (ANU). The isotopes measured and Faraday cup configuration used for the analysis are shown in Table 1. Data reduction was performed offline and includes a blank, 87Rb isobar and an exponential mass bias correction.

Table 1Standard cup configuration and analysed masses (amu or isotope mass) employed for solution strontium isotope analysis on the Neptune MC-ICP-MS at the Research School of Earth Sciences, the Australian National University.

Download Print Version | Download XLSX

To ensure precision, accuracy and reproducibility in the data produced, the samples were run in a sequence, which included blank and standard samples. The strontium carbonate reference material SRM987 (National Institute of Standards and Technology) was measured on the Neptune MC-ICP-MS during the sample sequence to quantify instrument drift. During the analysis period, measurements of SRM987 ranged from 0.71012±0.00001 to 0.71028±0.00001, with a mean of 0.71022 ± 0.00003 (n=32). The soil and rock results were not calibrated based on these SRM987 measurements.

Figure 2Bioavailable 87Sr 86Sr of rock samples, excluding samples with high 87Sr 86Sr ratios.

Figure 3Bioavailable 87Sr 86Sr of soil samples.

3 Results

A total of 60 soil samples and 48 rock samples were analysed for 87Sr 86Sr. The 87Sr 86Sr ratios of the soil samples range from 0.70577±0.00001 to 0.71020±0.00003. The 87Sr 86Sr ratios of the rock samples range from 0.70529±0.00001 to 0.74072±0.00001. These results are illustrated on satellite images of Israel in Fig. 2 (rock samples) and Fig. 3 (soil samples). Gross lithologies, defined by observations in the field and geological maps of the region (Sneh et al., 1998), are used to group the results of Sr isotope analysis. The soil and rock samples analysed show some differences in 87Sr 86Sr ratios between the different lithologies sampled, as illustrated in Fig. 4, although there is some overlap between most lithologies. Figure 4a displays all samples analysed, while Fig. 4b displays an inset removing a high 87Sr 86Sr granite sample (site IS044) and a high 87Sr 86Sr rhyolite sample (site IS047), to display the rest of the data set more clearly. Basalt samples are generally less radiogenic in 87Sr 86Sr ratio than the other lithologies sampled for both soil and rock samples. Soil samples from carbonate (limestone, dolostone, chalk and marl), granite, kurkar (aeolian quartz sandstone with carbonate cement), rhyolite, siliciclastic lithologies and areas with no bedrock have comparable 87Sr 86Sr ratios (Fig. 4b). The median 87Sr 86Sr ratios from rock samples are slightly lower than those from soils for basalt, carbonate and kurkar samples (Fig. 4b). The 87Sr 86Sr ratios measured from granite rock samples (two samples) and the rhyolite rock sample are substantially elevated compared to other units, and compared to the soil samples from the same lithologies.

Figure 4Bioavailable 87Sr 86Sr of rock and soil samples in Israel, grouped by gross lithology observed in the field. Points indicate individual samples. (a) All data, (b) subset of data with the two highest 87Sr 86Sr ratio samples, one each from the granite (site IS044) and rhyolite (site IS047) lithologies, removed to better display variability between other samples.


At 43 of the sample locations, both soil and rock samples were collected, and the 87Sr 86Sr ratios of these samples are compared in Fig. 5. As with Fig. 4, Fig. 5a illustrates the entire data set, while Fig. 5b has a granite rock sample (site IS044) and a rhyolite sample (site IS047) removed to better display the variation in the rest of the data set. The variation between soil and rock 87Sr 86Sr ratios for samples collected from the same site ranges from 0.00001 to 0.03130. This offset between rock and soil samples is also visualised in Fig. 6, in which the reference line indicates where rock and soil samples have the same 87Sr 86Sr ratios for both rock and soil samples (y=x). Points which lie above this reference line have higher 87Sr 86Sr ratios for soil samples than for rock samples, which in this study is the case for most basalt, carbonate, siliciclastic and kurkar samples. Points which lie below the y=x reference line have higher 87Sr 86Sr ratios for rock samples compared to soil samples; in this study this is the case for granite and rhyolite samples. The mean offsets from rock to soil 87Sr 86Sr ratios of samples collected from the same sites are shown in Table 2, grouped by lithology.

Sediments from seven archaeological sites were also analysed for 87Sr 86Sr: 4 from Amud, 1 from the Atlit Railway Bridge site, 2 from Neve David, 5 from Qafzeh, 3 from Sefunim Cave, 3 from Skhul and 16 from Tabun. These results are illustrated in Fig. 7, with samples R_1024 and R_1026 from Amud highly elevated compared to the rest of the data set (Fig. 7a). Figure 7b displays a subset of the data excluding these elevated samples to better display the variability in the rest of the samples. For the sites other than Amud, the samples from Qafzeh show the largest range in 87Sr 86Sr ratios between layers, while for Neve David, Sefunim Cave, Skhul and Tabun, most samples have similar 87Sr 86Sr ratios.

Table 2Mean offset by lithology from rock to soil 87Sr 86Sr ratios collected from the same sample locations.

Download Print Version | Download XLSX

4 Discussion

4.1 Comparison between rock and soil 87Sr 86Sr ratios

The difference between the 87Sr 86Sr ratios of soil and rock sampled from the same site is in most cases greater than the analytical uncertainty (Fig. 5), which may be due to inputs to the soil from sources other than bedrock, such as sea spray, irrigation or Saharan aeolian dust, or alternately due to differences in the weathering rates of minerals. Rainwater collected across Israel has been found to contain strontium associated with sea spray, marine minerals and dust from carbonate minerals, the relative proportions of which affect the 87Sr 86Sr ratio (Herut et al., 1993).

Figure 5Bioavailable 87Sr 86Sr of rock and soil samples collected from the same sample locations displayed with 2σ error bars. Sample locations are grouped based on gross lithology observed in the field. (a) All data, (b) subset of data with high 87Sr 86Sr granite sample (site IS044) and rhyolite sample (site IS047) removed to display variability between other samples.


The input of aeolian dust is expected to be particularly important to the 87Sr 86Sr ratios of these soil samples, as it has been reported that aeolian material may make up to 50 % of soils formed on hard limestone rocks in Israel (Yaalon, 1997). The 87Sr 86Sr isotope ratio of dust in the region is known to have varied over time, which is an important consideration in the use of modern isoscapes for the interpretation of archaeological studies. 87Sr 86Sr ratios of dust in this region varied from 0.711–0.712 during the marine isotope stage 2 and 4, and 0.709–0.710 during MIS 1, 3 and 5 (Haliva-Cohen et al., 2012). The delivery of aeolian dust is also affected by climate variations (Frumkin et al., 2011), with Sr isotope ratios found to increase at major climate transitions that correspond to sapropel formation in the Mediterranean (Stein et al., 2007). Speleothem analyses have shown that colder and drier conditions are associated with higher Sr isotope ratios due to an increase in sea spray and aeolian dust (Ayalon et al., 1999; Frumkin and Stein, 2004).

Figure 6Bioavailable 87Sr 86Sr of rock and 87Sr 86Sr of soil samples collected from the same sample locations. Sample locations are grouped based on gross lithology observed in the field. Points which lie along the reference line (y=x) have the same 87Sr 86Sr ratios for both soil and rock samples from the same sample site, while points above the line indicate a higher soil 87Sr 86Sr ratio than rock, and points which lie below this line indicate a rock 87Sr 86Sr ratio higher than the soil sample from the same site.


4.2 Sediment 87Sr 86Sr ratios from archaeological sites

The sediments analysed from archaeological sites in this study are mostly in carbonate units (limestone, dolostone, marl, chalk, chert) which are Mid-Eocene to Albian–Cenomanian in age, while the Atlit Railway Bridge site is located in natural caves formed in Quaternary kurkar (aeolian quartz sandstone with carbonate cement) on the coastal plain (Porat et al., 2018). The single sample analysed from the Atlit Railway Bridge site has a 87Sr 86Sr ratio of 0.70952±0.00001, which is slightly elevated compared to the other kurkar soils analysed in this study, which range from 0.70894 ± 0.00001 to 0.70934 ± 0.00013 (n=4). Soils associated from carbonate units in this study have 87Sr 86Sr ratios which range from 0.70851±0.00011 to 0.70925±0.00001 (n=35), although not all of the sediments from archaeological sites in this study from carbonate units fall within this range. Sediment samples from Tabun and Qafzeh in particular have variable 87Sr 86Sr ratios. Samples from stratigraphic unit B1 at Amud are the most varied; samples R_1024 and R_1026 have highly elevated 87Sr 86Sr ratios compared to samples R_1025 and R_1027. Unit B at Amud is composed of interbedded calcareous silt, calcareous concretions and ash derived from anthropogenic activity (Madella et al., 2002). These elevated 87Sr 86Sr ratios are likely due to the ash component in these sediments, and the variable nature of the interbedded sediments at this site.

Figure 787Sr 86Sr of soils collected from archaeological sites in Israel. (a) All data, (b) subset excluding the highest 87Sr 86Sr ratios of samples R_1024 and R_1026 from Amud to display the range of variability in the rest of the data set.


4.3 Comparison with regional data sets

Several previous studies analysing strontium isotopes have been undertaken in Israel, particularly in northern Israel and the Golan Heights region. The results obtained in this study of soils derived from basalts in the Golan Heights region have 87Sr 86Sr ratios in the range of 0.70577 to 0.70681 for samples (sample numbers: 175, 203, 204, 205 and 216) with robust geological provenance. Basalt rock samples (sample numbers: 119, 143, 144, 147 and 278) in this study from the Golan Heights region have 87Sr 86Sr ratios ranging from 0.70529 to 0.70681. For the same region, Shewan (2004) report 87Sr 86Sr ratios ranging from 0.70529 to 0.70571 (n=4), while the results of Spiro et al. (2011) of 87Sr 86Sr measured from water samples in the Hula Valley and Golan Heights region range from 0.70467 to 0.70790 (n=37). The small difference between the results of this study and those of Shewan (2004) may reflect the limited number of samples analysed, as the larger data set of Spiro et al. (2011) encompasses the results of both this study and Shewan (2004). Weinstein (2006) measured bedrock 87Sr 86Sr from basalt units in the region, ranging from 0.7031 to 0.7034, lower than those measured in this study. Hartman and Richards (2014) measured 87Sr 86Sr ratios from plants and invertebrates, and for samples taken from areas with basalt bedrock, ligneous (woody) plants were found to have 87Sr 86Sr ratios ranging from 0.70456 to 0.70851, non-ligneous plants ranging from 0.70473 to 0.70872, and invertebrates from 0.70494 to 0.70868. Rosenthal et al. (1989) summarise 87Sr 86Sr isotope sampling of basalt-derived groundwater in the Jordan Valley, with 87Sr 86Sr ratios of 0.7045 to 0.705. Rainwater in the Golan Heights region was found to have a 87Sr 86Sr isotope composition in the range of 0.70804 to 0.70923 (Herut et al., 1993). The wide range of 87Sr 86Sr isotope results in the Golan Heights region may be explained by the rapid depletion of Sr from rocks and soil in this region due to weathering and the large contribution of aeolian dust to soil profiles (Singer, 2007, 202–206).

Carbonate units were the most sampled lithology in this study, due to their widespread distribution across the country. This lithology groups a range of carbonate units (limestones, dolostones, chalk and marl), and the 87Sr 86Sr ratios of rock samples range from 0.70733 to 0.70911, while soil 87Sr 86Sr ratios range from 0.70803 to 0.71020. Hartman and Richards (2014) report bedrock 87Sr 86Sr ratios from northern Israel for carbonate units to range from 0.7073 to 0.7078, comparable to the lower range of rock samples measured in this study. Ligneous plant samples have 87Sr 86Sr ratios ranging from 0.70789 to 0.70910, non-ligneous plants from 0.70791 to 0.70919 and invertebrates from 0.70807 to 0.70876 (Hartman and Richards, 2014), comparable to the results of this study from rock and soil samples. Rosenthal et al. (1989) summarise 87Sr 86Sr isotopes from carbonate-derived water sources in the Jordan Valley as having 87Sr 86Sr ratios ranging from 0.7070 to 0.7080.

Arnold et al. (2016) sampled plants in the vicinity of Tell es-Safi/Gath to create a local bioavailable strontium map as a baseline for use in interpreting the mobility of domestic animals from archaeological sites in the region. Of the 10 plant samples collected, 4 were from kurkar soils, with 87Sr 86Sr ratios ranging from 0.70863 to 0.70881. The kurkar soil 87Sr 86Sr ratios reported in this study are slightly more elevated, ranging from 0.70894 to 0.70934. Hartman and Richards (2014) also measured 87Sr 86Sr ratios from plants and invertebrates from kurkar units, ligneous plants ranged from 0.70899 to 0.70903, non-ligneous plants ranged from 0.70900 to 0.70903, and two mollusc shells were measured with 87Sr 86Sr ratios of 0.70897 and 0.70904, comparable to the kurkar soil samples reported in this study. Coastal plain grasses were analysed by Shewan (2004) in the region of El Wad cave, south of Haifa. These 87Sr 86Sr ratios ranged from 0.70886 to 0.70965, with an outlying 87Sr 86Sr ratio of 0.71003 attributed to influence from the road base (Shewan, 2004). The results from the kurkar samples in this study are comparable to those of coastal plain grasses as reported by Shewan (2004).

Samples from granite, rhyolite and siliciclastic lithologies, as well as samples from regions with no bedrock, are also reported in this study. Hartman and Richards measured 87Sr 86Sr ratios from plants growing in alluvium, with a range of 0.70712 to 0.70830 for combined ligneous and non-ligneous (n=4) plants. From this study, samples sites where no bedrock was found had 87Sr 86Sr ratios ranging from 0.70783 to 0.70952, comparable to the results of Hartman and Richards (2014). There are no comparable studies which sample bioavailable 87Sr 86Sr for the granite and rhyolite samples included in this study.

5 Data availability

The data set can be viewed and downloaded on the IRHUM (Isotopic Reconstruction of Human Migration) database (, last access: 23 November 2020). The IRHUM database architecture and functionality are described by Willmes et al. (2014). The data are also available at the Open Science Framework data repository “Data Associated with Bioavailable Soil and Rock Strontium isotope data from Israel” (Moffat et al., 2020).

6 Conclusion

This data set represents a substantial contribution of bioavailable soil and rock 87Sr 86Sr ratios for Israel, to complement and build on previous research. The 87Sr 86Sr results are shown to be principally controlled by lithology and are in broad agreement with previous smaller-scale studies in the region. Soil and rock 87Sr 86Sr ratios from the same site are generally offset, and determining the mechanisms behind this offset could be the focus of future research to supplement this study.

Author contributions

RG conceived the study. IM collected sample material. IM, MW, GM, LK, LM, RA and MA undertook the analysis. IM and RR wrote the first draft and produced the figures. All authors reviewed the paper.

Competing interests

The authors declare that they have no conflict of interest.


Thank you to Avraham Ronen and Liora Horwitz for their assistance with field work in Israel and Steven Savage for providing the Landsat background for Figs. 2 and 3.

Financial support

This research has been supported by the Australian Research Council (grant nos. DP0664144, DP110101417, and DE160100703) and the Flinders University (Research Investment Fund Grant).

Review statement

This paper was edited by Attila Demény and reviewed by Jurian Hoogewerff and one anonymous referee.


Aggarwal, J., Habicht-Mauche, J., and Juarez, C.: Application of heavy stable isotopes in forensic isotope geochemistry: A review, Appl. Geochemistry, 23, 2658–2666,, 2008. 

Arnold, E. R., Hartman, G., Greenfield, H. J., Shai, I., Babcock, L. E., and Maeir, A. M.: Isotopic evidence for early trade in animals between Old Kingdom Egypt and Canaan, PLoS One, 11, 1–12,, 2016. 

Ayalon, A., Bar-Matthews, M., and Kaufman, A.: Petrography, strontium, barium and uranium concentrations, and strontium and uranium isotope ratios in speleothems as palaeoclimatic proxies: Soreq Cave, Israel, Holocene, 9, 715–722, 1999. 

Barnett-Johnson, R., Pearson, T. E., Ramos, F. C., Grimes, C. B., and Bruce MacFarlane, R.: Tracking natal origins of salmon using isotopes, otoliths, and landscape geology, Limnol. Oceanogr., 53, 1633–1642,, 2008. 

Bataille, C. P. and Bowen, G. J.: Mapping 87Sr 86Sr variations in bedrock and water for large scale provenance studies, Chem. Geol., 304–305, 39–52,, 2012. 

Bataille, C. P., Laffoon, J., and Bowen, G. J.: Mapping multiple source effects on the strontium isotopic signatures of ecosystems from the circum-Caribbean region, Ecosphere, 3, 1–24,, 2012. 

Bataille, C. P., von Holstein, I. C. C., Laffoon, J. E., Willmes, M., Liu, X. M., and Davies, G. R.: A bioavailable strontium isoscape for Western Europe: A machine learning approach, PLoS One, 13, 1–27,, 2018. 

Bataille, C. P., Crowley, B. E., Wooller, M. J., and Bowen, G. J.: Advances in global bioavailable strontium isoscapes, Palaeogeogr. Palaeoclimatol. Palaeoecol., 81, 104265,, 2020. 

Beard, B. L. and Johnson, C. M.: Strontium Isotope Composition of Skeletal Material Can Determine the Birth Place and Geographic Mobility of Humans and Animals, J. Forensic Sci., 45, 1049–1061,, 2000. 

Bentley, R. A.: Strontium isotopes from the earth to the archaeological skeleton: A review, J. Archaeol. Method Theory, 13, 135–187,, 2006. 

Capo, R. C., Stewart, B. W., and Chadwick, O. A.: Strontium isotopes as tracers of earth surface processes: theory and methods, Geoderma, 82, 197–225, 1998. 

Chadwick, O. A., Derry, L. A., Bern, C. R., and Vitousek, P. M.: Changing sources of strontium to soils and ecosystems across the Hawaiian Islands, Chem. Geol., 267, 64–76,, 2009. 

Dasch, E. J.: Strontium isotopes in weathering profiles, deep-sea sediments, and sedimentary rocks, Geochim. Cosmochim. Act., 33, 1521–1552,, 1969. 

Dickin, A. P.: Radiogenic Isotope Geology, Cambridge University Press, Cambridge, 2000. 

DIN: ISO 19730: Soil quality – Extraction of trace elements from soil using ammonium nitrate solution (ISO 19730:2008), DIN, 2009. 

Evans, J. A., Montgomery, J., Wildman, G., and Boulton, N.: Spatial variations in biosphere 87Sr 86Sr in Britain, J. Geol. Soc. London., 167, 1–4,, 2010. 

Faure, G.: Principles of Isotope Geology, John Wiley and Sons, New York, 1986. 

Faure, G. and Mensing, T.: Isotopes: Principals and Applications, 3rd edn., John Wiley and Sons Inc., 2005. 

Faure, G. and Powell, J.: Strontium Isotope Geology, Springer, Berlin, Heidelberg, 1972. 

Frei, K. M. and Frei, R.: The geographic distribution of strontium isotopes in Danish surface waters – A base for provenance studies in archaeology, hydrology and agriculture, Appl. Geochemistry, 26, 326–340,, 2011. 

Frumkin, A. and Stein, M.: The Sahara-East Mediterranean dust and climate connection revealed by strontium and uranium isotopes in a Jerusalem speleothem, Earth Planet. Sc. Lett., 217, 451–464,, 2004. 

Frumkin, A., Bar-Yosef, O., and Schwarcz, H. P.: Possible paleohydrologic and paleoclimatic effects on hominin migration and occupation of the Levantine Middle Paleolithic, J. Hum. Evol., 60, 437–451,, 2011. 

Goede, A., McCulloch, M., McDermott, F., and Hawkesworth, C.: Aeolian contribution to strontium and strontium isotope variations in a tasmanian speleothem, Chem. Geol., 149, 37–50,, 1998. 

Haliva-Cohen, A., Stein, M., Goldstein, S. L., Sandler, A., and Starinsky, A.: Sources and transport routes of fine detritus material to the Late Quaternary Dead Sea basin, Quat. Sci. Rev., 50, 55–70,, 2012. 

Hartman, G. and Richards, M.: Mapping and defining sources of variability in bioavailable strontium isotope ratios in the Eastern Mediterranean, Geochim. Cosmochim. Act., 126, 250–264,, 2014. 

Herut, B., Starinsky, A., and Katz, A.: Strontium in rainwater from Israel: Sources, isotopes and chemistry, Earth Planet. Sc. Lett., 120, 77–84,, 1993. 

Hobson, K. A., Barnett-Johnson, R., and Cerling, T.: Using Isoscapes to Track Animal Migration, in: Isoscapes: Understanding Movement, Pattern, and Process on Earth Through Isotope Mapping, edited by: West, J. B., Bowen, G. J., Dawson, T. E., and Tu, K. P., Springer, 273–298, 2010. 

Hodell, D. A., Quinn, R. L., Brenner, M., and Kamenov, G.: Spatial variation of strontium isotopes (87Sr 86Sr) in the Maya region: A tool for tracking ancient human migration, J. Archaeol. Sci., 31, 585–601,, 2004. 

Hoogewerff, J. A., Reimann, C., Ueckermann, H., Frei, R., Frei, K. M., van Aswegen, T., Stirling, C., Reid, M., Clayton, A., and Ladenberger, A.: Bioavailable 87Sr 86Sr in European soils: A baseline for provenancing studies, Sci. Total Environ., 672, 1033–1044,, 2019. 

Madella, M., Jones, M., Goldberg, P., Goren, Y., and Hovers, E.: The Exploitation of Plant Resources by Neanderthals in Amud Cave (Israel): The Evidenve from Phytolith Studies, J. Archaeol. Sci., 29, 703–719,, 2002. 

Maurer, A. F., Galer, S. J. G., Knipper, C., Beierlein, L., Nunn, E. V., Peters, D., Tütken, T., Alt, K. W., and Schöne, B. R.: Bioavailable 87Sr 86Sr in different environmental samples – Effects of anthropogenic contamination and implications for isoscapes in past migration studies, Sci. Total Environ., 433, 216–229,, 2012. 

McDermott, F. and Hawkesworth, C.: The evolution of strontium isotopes in the upper continental crust, Nature, 344, 850–853, 1990. 

Moffat, I., Stringer, C., and Grün, R.: Spatially Resolved LA-MC-ICPMS Strontium Isotope Microanalysis of Archaeological Fauna, Palaios, 27, 667–670,, 2012. 

Moffat, I., Rudd, R., Willmes, M., Mortimer, G., Kinsley, L., McMorrow, L., and Grün, R.: Data Associated with “Bioavailable Soil and Rock Strontium Isotope Data from Israel', Open Science Framework,, 2020. 

Neat, P., Faure, G., and Pegram, W.: The isotopic composition of strontium in non-marine carbonate rocks: The Flagstaff Formation of Utah, Sedimentology, 26, 271–282, 1979. 

Perry, M. A., Coleman, D., and Delhopital, N.: Mobility and exile at 2nd century A.D. Khirbet edh-Dharih: Strontium isotope analysis of human migration in Western Jordan, Geoarchaeology, 23, 528–549,, 2008. 

Perry, M. A., Coleman, D. S., Dettman, D. L., and Al-Shiyab, A. H.: An isotopic perspective on the transport of byzantine mining camp laborers into southwestern Jordan, Am. J. Phys. Anthropol., 140, 429–441,, 2009. 

Porat, N., Jain, M., Ronen, A., and Horwitz, L. K.: A contribution to late Middle Paleolithoc chronology of the Levant: New luminescence ages fo the Atlit Railway Bridge site, Coastal Plain, Israel, Quat. Int., 464, 32–42,, 2018. 

Price, T. D., Burton, J. H., and Bentley, R. A.: The characterization of biologically available strontium isotope ratios for the study of prehistoric migration, Archaeometry, 44, 117–135,, 2002. 

Rosenthal, Y., Katz, A., and Tchernov, E.: The reconstruction of quaternary freshwater lakes from the chemical and isotopic composition of gastropod shells: The dead sea rift, Israel, Palaeogeogr. Palaeoclimatol. Palaeoecol., 74, 241–253,, 1989. 

Shewan, L.: Natufian settlement systems and adaptive strategies: The issue of sedentism and the potential of strontium isotope analysis, in The Last Hunter-Gatherer Societies in the Near East, edited by: Delage, C., BAR International Series, John and Erica Hedges, Oxford, 2004. 

Singer, A.: The Soils of Israel, Springer, Berlin, 2007. 

Slovak, N. M. and Paytan, A.: Applications of Sr Isotopes in Archaeology, in Handbook of Environmental Isotope Geochemistry, edited by: Baskaran, M., 743–768, 2012. 

Sneh, A., Bartov, Y., Weissbrod, T., and Rosensaft, M.: Geological Map of Israel, Geological Survey of Israel, 4 pp., 1998. 

Sneh, A. and Rosensaft, M.: Major exposed lithologic units of Israel and environs (1 : 500 000), Geological Survey of Israel, 2014. 

Spiro, B., Ashkenazi, S., Starinsky, A., and Katz, A.: Strontium isotopes in Melanopsis sp. as indicators of variation in hydrology and climate in the Upper Jordan Valley during the Early-Middle Pleistocene, and wider implications, J. Hum. Evol., 60, 407–416,, 2011. 

Stein, M., Starinsky, I. A., Katz, I. A., Goldstein, J. S. L., Machlus, M., and Schramm, A.: Strontium isotopic, chemical, and sedimentological evidence for the evolution of Lake Lisan and the Dead Sea, Geochim. Cosmochim. Act., 61, 3975–3992, 1997. 

Stein, M., Almogi-Labin, A., Goldstein, S. L., Hemleben, C., and Starinsky, A.: Late Quaternary changes in desert dust inputs to the Red Sea and Gulf of Aden from 87Sr 86Sr ratios in deep-sea cores, Earth Planet. Sc. Lett., 261, 104–119,, 2007. 

Stewart, B. W., Capo, R. C., and Chadwick, O. A.: Quantitative strontium isotope models for weathering, pedogenesis and biogeochemical cycling, Geoderma, 82, 173–195,, 1998. 

Voerkelius, S., Lorenz, G. D., Rummel, S., Quétel, C. R., Heiss, G., Baxter, M., Brach-Papa, C., Deters-Itzelsberger, P., Hoelzl, S., Hoogewerff, J., Ponzevera, E., Van Bocxstaele, M., and Ueckermann, H.: Strontium isotopic signatures of natural mineral waters, the reference to a simple geological map and its potential for authentication of food, Food Chem., 118, 933–940,, 2010. 

Weinstein, Y., Navon, O., Altherr, R., and Stein, M.: The Role of Lithospheric Mantle Heterogeneity in the Generation of Plio-Pleistocene Alkali Basaltic Suites from NW Harrat Ash Shaam (Israel), J. Petrol., 47, 1017–1050,, 2006. 

Willmes, M., McMorrow, L., Kinsley, L., Armstrong, R., Aubert, M., Eggins, S., Falguères, C., Maureille, B., Moffat, I., and Grün, R.: The IRHUM (Isotopic Reconstruction of Human Migration) database – bioavailable strontium isotope ratios for geochemical fingerprinting in France, Earth Syst. Sci. Data, 6, 117–122,, 2014. 

Willmes, M., Bataille, C. P., James, H. F., Moffat, I., McMorrow, L., Kinsley, L., Armstrong, R. A., Eggins, S., and Grün, R.: Mapping of bioavailable strontium isotope ratios in France for archaeological provenance studies, Appl. Geochemistry, 90, 75–86,, 2018. 

Yaalon, D. H.: Soils in the Mediterranean region: What makes them different?, Catena, 28, 157–169,, 1997. 

Short summary
This research has measured the bioavailable strontium isotope composition of soil and rock samples from the major geological units in Israel. Plants, animals and people take up this value into their bones and teeth via food and water and so mapping these values has important implications for understanding where archaeological, ecological, food science and forensic samples are from.