Strontium and Oxygen Isotope Analyses Reveal Late Cretaceous Shark Teeth in Iron Age Strata in the Southern Levant

Skeletal remains in archaeological strata are often assumed to be of similar ages. Here we show that combined Sr and O isotope analyses can serve as a powerful tool for assessing fish provenance and even for identifying fossil fish teeth in archaeological contexts. For this purpose, we established a reference Sr and O isotope dataset of extant fish teeth from major water bodies in the Southern Levant. Fossil shark teeth were identified within Iron Age cultural layers dating to 8–9th century BCE in the City of David, Jerusalem, although the reason for their presence remains unclear. Their enameloid 87Sr/86Sr and δ18OPO4 values [0.7075 ± 0.0001 (1 SD, n = 7) and 19.6 ± 0.9‰ (1 SD, n = 6), respectively], are both much lower than values typical for modern marine sharks from the Mediterranean Sea [0.7092 and 22.5–24.6‰ (n = 2), respectively]. The sharks’ 87Sr/86Sr are also lower than those of rain- and groundwater as well as the main soil types in central Israel (≥0.7079). This indicates that these fossil sharks incorporated Sr (87Sr/86Sr ≈ 0.7075) from a marine habitat with values typical for Late Cretaceous seawater. This scenario is in line with the low shark enameloid δ18OPO4 values reflecting tooth formation in the warm tropical seawater of the Tethys Ocean. Age estimates using 87Sr/86Sr stratigraphy place these fossil shark teeth at around 80-million-years-old. This was further supported by their taxonomy and the high dentine apatite crystallinity, low organic carbon, high U and Nd contents, characteristics that are typical for fossil specimens, and different from those of archaeological Gilthead seabream (Sparus aurata) teeth from the same cultural layers and another Chalcolithic site (Gilat). Chalcolithic and Iron Age seabream enameloid has seawater-like 87Sr/86Sr of 0.7091 ± 0.0001 (1 SD, n = 6), as expected for modern marine fish. Fossil shark and archaeological Gilthead seabream teeth both preserve original, distinct enameloid 87Sr/86Sr and δ18OPO4 signatures reflecting their different aquatic habitats. Fifty percent of the analysed Gilthead seabream teeth derive from hypersaline seawater, indicating that these seabreams were exported from the hypersaline Bardawil Lagoon in Sinai (Egypt) to the Southern Levant since the Iron Age period and possibly even earlier.


INTRODUCTION
Past fish habitats have traditionally been reconstructed based on taxonomic identification of fish remains recovered in archaeological sites (e.g., Wheeler and Jones, 1989;Zohar and Biton, 2011;Zohar, 2017). Recently, stable isotope analyses of fish bioapatite have also been used as palaeoenvironmental proxies for fish habitat and provenance (e.g., Dufour et al., 2007;Otero et al., 2011;Leuzinger et al., 2015;Sisma-Ventura et al., 2015, 2018Fetner and Iwaszczuk, 2020). In this study, we analyse the Sr and O isotope composition of teeth from modern fish derived from different marine and freshwater bodies in the Levant to provide a reference dataset for assessing the provenance of ancient fish teeth from Chalcolithic and Iron Age cultural layers in Israel.

The Use of Strontium and Oxygen Isotopes as Provenance Proxies for Fish
The strontium isotope composition ( 87 Sr/ 86 Sr) of the wellpreserved bioapatite of fish tooth enameloid has been widely used as a chemostratigraphic dating method for marine sediments (e.g., Ingram, 1995;Kocsis et al., 2009;Harrell et al., 2016) as a proxy to infer the palaeosalinity levels of oceanic basins and brackish water bodies (e.g., Schmitz et al., 1991Schmitz et al., , 1997Bryant et al., 1995;Reinhardt et al., 1998;Kocsis et al., 2009), as well as for tracing past freshwater habitat use in sharks (e.g., Fischer et al., 2013) and to track their migration from the marine realm into fresh water systems (Kocsis et al., 2007. This is possible because the 87 Sr/ 86 Sr of ambient (sea)water is incorporated into the bioapatite of the teeth of Elasmobranchii (subclass of cartilaginous fish) and Teleostei (bony fish) in vivo, without any apparent measurable isotope fractionation (Vennemann et al., 2001). Marine fish tooth enameloid thus has a Sr isotope composition that reflects seawater 87 Sr/ 86 Sr, which today is a globally constant value of 0.70918 (Hodell et al., 1990;Mokadem et al., 2015;McArthur et al., 2020), but has varied over geological time (i.e., Phanerozoic) between values of around 0.7072-0.7092 (McArthur et al., 2020). Across the last glacial-interglacial cycle the seawater 87 Sr/ 86 Sr did not deviate measurably from the modern-day value ( 87 Sr/ 86 Sr = 0.7091792 ± 0.00000212 SE, n = 17; Mokadem et al., 2015). In contrast, the 87 Sr/ 86 Sr of Sr dissolved in surface freshwater bodies (rivers and lakes) depends both on the Sr isotope composition of the local bedrock (age and lithology; e.g., Goldstein and Jacobsen, 1987) and groundwater and precipitation (Herut et al., 1993;Chadwick et al., 1999;Vengosh et al., 1999) in the drainage area of rivers and lakes. Freshwater will thus have 87 Sr/ 86 Sr below or above the contemporaneous global seawater value (Dufour et al., 2007). The 87 Sr/ 86 Sr in fish teeth can be used as a proxy to distinguish between freshwater and marine habitats (Figure 1A; e.g., Schmitz et al., 1997;Kocsis et al., 2007Kocsis et al., , 2014Tütken et al., 2011;Fischer et al., 2013), whereas for brackish water bodies the 87 Sr/ 86 Sr depends on the sea-to-freshwater mixing ratio (Bryant et al., 1995;Reinhardt et al., 1998).
The phosphate oxygen isotope composition ( 18 O/ 16 O expressed as δ 18 O PO4 value) of bony fish and shark (we refer to these collectively as 'fish' elsewhere in the paper) tooth enameloid has been used to reconstruct past seawater temperatures (e.g., Kolodny and Raab, 1988;Kolodny and Luz, 1991;Pucéat et al., 2003), seawater salinity changes (Zacke et al., 2009;Leuzinger et al., 2015) and past aquatic fish habitats such as rivers, lakes, and lagoons (Kolodny et al., 1983;Dufour et al., 2007;Klug et al., 2010;Otero et al., 2011;Kocsis et al., 2015;Leuzinger et al., 2015;Sisma-Ventura et al., 2015, 2018. Fish bones and teeth are mineralised in oxygen isotope equilibrium with their body fluids. The body fluids of fish have the same oxygen isotope composition as the ambient water (δ 18 O Water ) in which they reside and the body (i.e., tooth formation) temperature of fish is identical to the ambient water temperature ( Figure 1A). Thus, the bioapatite δ 18 O PO4 values of fish bones and teeth are controlled by the local δ 18 O Water value and water temperature according to the temperature-dependent water-phosphate oxygen isotope fractionation (Longinelli and Nuti, 1973;Kolodny et al., 1983;Pucéat et al., 2010;Lécuyer et al., 2013). The δ 18 O Water values of the modern oceans are ≈ 0.0 in tropical and subtropical regions and are more positive in evaporative basins such as the Mediterranean Sea (≈ 1.5 ), while fresh-and meltwater influx at high latitudes results in negative surface water values (between −2.0 and −4.0 ; Rohling, 2007).
The 87 Sr/ 86 Sr of rivers and lakes is controlled by dissolved Sr released from the weathering of bedrock in catchment areas (Goldstein and Jacobsen, 1987;Woodward et al., 2015). This enables us to distinguish fish from bodies of freshwater draining bedrocks of different ages and lithologies, which have distinct 87 Sr/ 86 Sr ( Figure 1A). Brackish water, which is a mixture of seawater and freshwater, will have variable 87 Sr/ 86 Sr, depending on the mixing ratio of the waters, Sr content, and the 87 Sr/ 86 Sr of the freshwater input (e.g., Bryant et al., 1995;Reinhardt et al., 1998;Hobbs et al., 2019). Brackish water is typical in estuaries of rivers and coastal wetlands.
Combined Sr and O isotope analysis of fish remains therefore represents a promising geochemical approach for assessing fish provenance in archaeological and palaeontological contexts. In this study we establish a reference dataset of Sr and O isotope values by analysing the teeth of modern fish from different water bodies in Egypt and Israel. We then use this reference dataset to reconstruct past fish habitats by analysing fish remains from Chalcolithic (Gilat) and Iron Age (City of David, Jerusalem) archaeological sites in the Southern Levant.

Fish as a Commodity
Identification of fish trade in antiquity is traditionally based on the presence of "exotic" fish remains in archaeological sites, which are distantly located from the original habitat of the fish (Davis, 1985; . In the Eastern Mediterranean region, several key species have been identified as indicators of fish trade as they occur only in one of three main aquatic habitats: the Nile (Egypt), Eastern Mediterranean and the Red Sea ( Figure 1B). For example, Nilotic species such as Lates niloticus (Latidae, Nile perch), Oreochromis niloticus (Cichlidae; Nile tilapia), Bagrus sp. (Bagridae, bagrid catfish) (Van Neer et al., 1997, 2005Van Neer and Depraetere, 2005) are considered typical indicators of fish trade from Egypt to the Southern Levant. Recently, stable isotope studies of archaeological fish teeth demonstrated that Gilthead seabream (Sparidae; Sparus aurata) were also exported from Egypt to this area of modern-day Israel as they derive from the hypersaline Bardawil Lagoon in northern Sinai (Figure 2A) and not from local fishing activity (Sisma-Ventura et al., 2015, 2018. The occurrence of imported fish in archaeological layers demonstrates that Nilotic and lagoonal fish were common trade goods from the Bronze Age onwards (Van Neer and Ervynck, 2004;, 2005Raban-Gerstel et al., 2008;Sisma-Ventura et al., 2015, 2018Zohar and Artzy, 2019).
However, most of these studies have concentrated on bony fish, while shark (Elasmobranchii) exploitation and trade routes have received little attention. In general, Elasmobranchii (Chondrichthyes) do not preserve well because, with the exception of their teeth, dermal denticles, spines and vertebrae centra, the rest of their skeleton is composed of cartilage (Wheeler and Jones, 1989;Rick et al., 2002). Teeth and vertebrae centra are typically the main shark remains found in archaeological contexts. While shark vertebrae centra are commonly found and usually well preserved, shark teeth are rarely recovered, even when the soil is carefully sieved ( Table 1 and references therein). Interestingly, previous studies (Cione and Bonomo, 2003;Betts et al., 2012;Charpentier et al., 2020) determined that fossil shark teeth recovered from archaeological contexts have been employed as ornaments and tools, while shark vertebrae centra are only rarely identified as jewelry (Boulanger et al., 2020) and are usually interpreted as food remnants.
A review of published data on ca. 30,000 fish remains recovered in Holocene archaeological sites in the Southern Levant (Table 1) shows that 511 elasmobranch remains (sharks or rays), were recovered from 37 sites dated from the Pre-Pottery Neolithic A (ca. 10,000 BC) until the Crusader period (ca. 1099 AD). Most of the remains are vertebrae centra, and the earliest appearance of elasmobranchs occur in the Pre-Pottery Neolithic A (Sultanian) level of the inland site of Hatoula , and in the Pre-Pottery Neolithic C (PPNC) coastal sites of Ashkelon (Garfinkel et al., 2005;Lernau, 2008) and Atlit-Yam (Galili et al., 2004). Shark remains from the Late Chalcolithic are reported only from Namir road (Van den Brink et al., 2016) and from the Early Bronze Age layer of coastal Ashkelon Afridar (Lernau, 2004). Elasmobranch remains occur in most sites dating to the Late Bronze Age until the Crusader period, regardless of the sieving method applied. However, in all cases shark teeth are rarely recovered (8.4% of all shark remains), and out of 43 identified shark teeth, 38 (88%) were excavated in the City of David, Jerusalem (Table 1).
From an economic and dietary perspective, elasmobranchs (sharks and rays) are regarded as attractive food source since they are rich in protein, vitamin A, and oil. Additionally, their skin may be used as sandpaper and their spines and teeth (including FIGURE 2 | (A) Map of the main water bodies and potential fish habitats of Egypt and the Southern Levant: Mediterranean Sea, Red Sea, upper Nile River, the brackish Lake Manzala on the Nile delta, the hypersaline Bardawil Lagoon in northern Sinai, and the freshwater Lake Kinneret and Jordan River, its main tributary, in northern Israel. These represent the key modern and past fish habitats from which we analysed modern fish samples to establish a Sr and O isotope reference dataset with which to assess the provenance of archaeological fish specimens. The ancient (i.e., fossil and historic) fish teeth were analysed from three different archaeological sites (black circles): Gilat (Chalcolithic), Ashkelon (Iron Age I), and Jerusalem/City of David (Iron Age I-II). Localities where fossil Cretaceous shark teeth were found in Israel (stars) after Kolodny and Raab (1988) and Retzler et al. (2013): 1 Ein Yabrud, 2 Letaot Mt., 3 Mishor Rotem, 4 HaMakhtesh HaKatan, 5 Ein Aqrabbim, 6 Hor Hahar, 7 Makhtesh Ramon. (B) Location map of the "Rock Cut Pool" adjacent to the Gihon spring, City of David, Jerusalem, and the nearby Siloam stream. (C) Photograph of the eastern wall of the "Rock Cut Pool" with the limestone boulders of the infill. Both (B,C) modified after Reich et al. (2007Reich et al. ( , 2008 and Sneh et al. (2010).
Frontiers in Ecology and Evolution | www.frontiersin.org fossil teeth) can be used as raw material for tools (Borhegyi, 1961;Noe-Nygaard, 1971;Luomala, 1984;Rick et al., 2002). Therefore, despite rare preservation in archaeological sites, when recovered in cultural contexts shark teeth have been interpreted to reflect exploitation of coastal water bodies, lagoons and estuaries and shed new light on the habitat the sharks exploited and fishing practices in the past (Serena, 2005).

The Archaeological Assemblage
Large numbers of fish remains (>10,000) were recovered in the City of David, located outside and south of the southern walls of Jerusalem Old City (present day Silwan village) (Mazar, 2015). These remains were recovered from a structure located adjacent to the Gihon spring named "The Rock Cut Pool" and dated to the Iron Age II period, 8-7th century BCE ( Figure 2B; Reich et al., 2007Reich et al., , 2008De Groot and Fadida, 2011;Reich and Shukron, 2011). This large rectangular structure (15 × 10 m) cut into the limestone bedrock has straight, almost vertical rock walls, and during the Bronze Age period, prior to the operation of Siloam water Tunnel, was probably part of a complex underground water system allowing access to water during times of siege. During the end of the 9th century BCE/beginning of the 8th century BCE, it was converted to a private house. For this purpose, the pool was filled with large limestone boulders and soil ( Figure 2C) in order to raise the floor-level by 3 m, and then sealed by the floor of the building above it. During excavation, the sediments beneath the floor were dry and wet sieved through a 5 mm mesh. More than 10,000 fish remains, together with 6.5 tons of pottery and several hundred broken clay seals, known as "bullae" were recovered and dated to Iron Age IIB. The "bullae" were used to seal letters or packages and thus would have been broken and thrown away when the packages and letters were opened. Fish taxonomic composition included 14 families of fish, most of them (ca. 90%) originating from the Mediterranean Sea (mainly Sparidae and Mugilidae) as well as some from freshwater habitats (ca. 8.5%) including the Nile (Reich et al., 2007(Reich et al., , 2008. Of these fish, 39 teeth and vertebrae centra belonged to Chondrichthyes fish (sharks and rays; Table 1). Several of these shark teeth were analysed for this study alongside Sparus aurata (Gilthead seabream) teeth from this site and from the Chalcolithic site of Gilat (5800-5500 BP; Rowan and Golden, 2009).

Shark Tooth Formation
Extant elasmobranchs (a subclass of cartilaginous fish that includes sharks, skates and rays) exhibit huge diversity in tooth shape. Their dentition is comprised of a large number of highly specialised single teeth that are organised in rows, and undergo continuous replacement throughout the life of the shark ( Figure 3A; Debiais-Thibaud et al., 2007Jernvall and Thesleff, 2012;Rasch et al., 2016;Meredith Smith et al., 2018). Tooth replacement and dentition patterns are distinctive for sharks and rays Meredith Smith et al., 2018). In sharks, only teeth in the outer row are fully upright and functional while replacement teeth lie posteriorly with their tips angled toward the rear. Tooth development and regeneration in Chondrichthyes is still poorly understood and characterised from a developmental and genetic perspective (Rasch et al., 2016). In different modern shark species, tooth replacement cycles vary and can be as rapid as 8-10 days, as in Lemon sharks (Negaprion brevirostris), 28 days as in Nurse sharks (Ginglymostoma cirratum), 38 days as in Leopard sharks (Triakis semifasciata) or as long as 150 days as in Sandy dogfish (Scyliorhinus canicular) (Moss, 1967;Reif et al., 1978;Luer et al., 1990;Botella et al., 2009 and references therein). Thus, while all sharks replace their teeth continuously, the timing of tooth formation and replacement varies depending on the species and age of the individual (younger sharks replace their teeth faster and thus more frequently). The teeth of sharks and bony fish consist of a compact outer layer of hard, highly mineralised enameloid and an interior core of softer, less mineralised dentine ( Figure 3A; Enault et al., 2015;Jambura et al., 2018Jambura et al., , 2020. The tooth root is made up of porous osteodentine, a specific kind of dentine, which contains dentinal osteons and inter-osteonal tissue, and superficially resembles osteonal bone Jambura et al., 2018Jambura et al., , 2020. Enameloid is harder and has a lower solubility than hydroxylapatite characteristic of enamel, because the bioapatite often has a high fluorine content (i.e., fluorapatite; Miake et al., 1991;Suga et al., 1993;Enax et al., 2014). Thus, enameloid is highly resistant to diagenetic alteration following post mortem burial, even over geological time periods (Sharp et al., 2000). In Chondrichthyes, enameloid is composed of two units: a dense tissue of single crystallite enameloid (SCE), and an inner tissue of bundled crystallite enameloid (BCE) (Figure 3A; Enault et al., 2015). In this study, we sampled material from SCE, wherever possible, but in some cases sampled a mixture of both of these layers. Except for specimens (a,f), seven of these fossil shark teeth were analysed for their Sr and O isotope composition (the unidenfied lamniform shark (a) and specimen (f) did not contain enough enameloid to sample). All of these teeth are presumably from the Late Cretaceous (

MATERIALS AND METHODS
The provenance of archaeological fish remains (sharks and bony fish) recovered in the Southern Levant both from inland and coastal sites (Figure 2) was examined by analysing 87 Sr/ 86 Sr and δ 18 O PO4 of the bioapatite of their teeth (i.e., enameloid). To assess past fish habitats and provenance we first establish a reference database of 87 Sr/ 86 Sr and δ 18 O PO4 values for modern sharks as well as marine and freshwater bony fish from the major water bodies in Egypt and in the Southern Levant (Nile, Lake Kinneret, Jordan River, Bardawil Lagoon, Figure 2A and Table 2). The isotope data from archaeological fish teeth are then compared with those of modern fish to assess their origin. Additional δ 18 O values and 87 Sr/ 86 Sr of rainwater and bedrock for each archaeological site are considered in order to infer the soil water Sr and O isotope composition and identify possible effects of post mortem diagenetic alteration. Additionally, the dentine apatite crystallinity, organic carbon content and trace element composition were analysed to further characterise the degree of diagenesis.

Archaeological Assemblage
Analysed archaeological material included six Sparus aurata teeth recovered from the Chalcolithic site of Gilat (Rowan and Golden, 2009;Van den Brink et al., 2016) and from the Iron Age period of the City of David (Jerusalem), as well as 10 shark teeth from the Iron Age period of the City of David, Jerusalem (Reich et al., 2007(Reich et al., , 2008 (Table 3 and Figure 3B).

Taxonomic Identification
Identification of S. aurata was performed using the first molariform tooth (Sisma-Ventura et al., 2019). Photos of these teeth are provided in the supplements (Supplementary Figure S1). Shark teeth were identified by Jürgen Kriwet from the University of Vienna. Two of these teeth were identified to the order of Lamniformes (tooth a and c in Figure 3B), the latter probably from the family Odontaspididae (SK-79, tooth c in Figure 3B), and two teeth (SK-85, tooth d in Figure 3B and S07815 Supplementary Figure S1B) were identified as Squalicorax an extinct taxon only occurring in the Late Cretaceous Cicimurri, 2005, 2006;Cappetta, 2012;Siversson et al., 2018). All identified taxa represent extinct fossil sharks belonging to the Late Cretaceous. The other shark teeth could not be identified due to poor preservation but displayed the same macroscopic appearance and preservation as the other fossil teeth (Supplementary Figure S1).

Isotope Analysis of Dental Tissue
Enameloid was carefully separated from the dentine by manual drilling. For S. aurata the cleaned tooth crown was crushed in an agate mortar after dentine removal to obtain a homogenous enameloid powder that was used for O and solution mode Sr isotope analysis. In situ LA-MC-ICP-MS Sr isotope analyses were performed directly on the tooth surface (in the SCE layer) of the shark teeth ( Figure 3A). The teeth of modern Conger eel, Tilapia, Nile Perch and Gulper shark were crushed whole (due to their very small size) and analysed for Sr and O isotopes.

Screening for Diagenetic Alteration
The preservation of fish teeth excavated from the Chalcolithic and Iron Age cultural layers was determined using various analytical methods including X-ray diffraction analysis (XRD), total organic carbon (TOC) content determination, and laser ablation inductively coupled plasma mass spectrometry (LA-ICP-MS) analysis of trace elements. These analyses were performed on the dentine of the shark and Gilthead seabream teeth, as dentine is much more sensitive to diagenetic alteration than enamel(oid) (Ayliffe et al., 1994;Kohn et al., 1999;Tütken et al., 2008). Diagenesis leads to recrystallization of the bioapatite, loss of collagen (i.e., organic matter), and post mortem uptake of fluorine, uranium and rare earth elements (REE) (Staudigel et al., 1985;Toyoda and Tokonami, 1990;Kohn et al., 1999). Therefore, million-year-old fossil fish teeth are expected to have a much higher apatite crystallinity (i.e., larger apatite crystallite size), lower TOC and up to several orders of magnitude higher REE and U content in the dentine relative to the archaeological fish teeth (which are only a few thousand years old). In modern fish teeth, REE and U occur in very low ng/g-level concentrations (Vennemann et al., 2001;Kocsis et al., 2015; this study) -often below the detection limit of LA-ICP-MS -and are only strongly enriched in cases of post mortem alteration. Such mineralogical and geochemical changes enable us to clearly distinguish whether fish teeth are contemporaneous to the cultural layer being excavated, or are fossil fish teeth. For comparison, dentine from a modern Great white shark and Gilthead seabream were analysed with the same methods as unaltered reference points.

X-Ray Diffraction (XRD) Analysis
Dentine powder samples were placed on a monocrystalline, background-free silicon sample holder and X-ray diffraction patterns were analysed using a Seifert XRD 3000 TT diffractometer using a Cu cathode at 40 kV and 30 mA and Cu K α radiation of 0.15406 nm with 0.03 • (2 ) steps in 2 s over the range from 5 • to 70 • 2 at the Institute of Geosciences, University of Mainz.

Total Organic Carbon (TOC) Analysis
The TOC content was determined in two steps. First, by analysing the dentine powders (about 10-15 mg/per sample) for total carbon (TC) content using a CN-analyser (Vario EL Cube, Elementar GmbH, Hanau, Germany) at the Institute for Geography, University of Mainz. Then, the structurally bound carbonate total inorganic carbon content (TIC) of dentine was measured using a Gasbench II coupled to a Delta V continuous flow isotope ratio gas mass spectrometer at the Max Planck Institute for Chemistry, Mainz. An internal modern elephant (Loxodonta africana) enamel standard (AG-Lox) of known TIC content (7.5 wt.%; Wacker et al., 2016) was analysed in different amounts to establish a regression between TIC and peak area of the mass-44 CO 2 peak. The dentine TOC content was then calculated as the difference between the TC and TIC content.

LA-ICP-MS Analysis of Trace Elements
Shark and Gilthead seabream teeth were embedded in epoxy resin to produce polished one-inch mounts. To determine the trace element compositions of enameloid and dentine, line scans were measured across the teeth (including both enameloid and dentine) by laser ablation (Supplementary Figure S2). The trace element measurements were performed using an ArF Excimer laser system (193 nm wavelength, NWR193 by ESI/NewWave) equipped with a TwoVol2 ablation cell and coupled to an Agilent 7500ce quadrupole ICP-MS at the Institute of Geoscience, Mainz. Laser repetition rate was set to 10 Hz with a transition speed of 10 µm/s, using a fluence of 3.5 J/cm 2 . Line scans were carried out with a rectangular spot size of 130 µm × 30 µm. Pre-ablation was performed prior to each analysis to clean the surface. Background signals were acquired for 15 s during laser warm-up prior to each scan. 43 Ca was used as internal standard and NIST SRM 612 was used for calibration, using the preferred values from the GeoReM Database (Jochum et al., 2005). NIST SRM 610, NIST SRM 1400 and NIST SRM 1486 were used as quality control materials to test accuracy and reproducibility of the measurements (Supplementary Table S1). All reference materials were measured prior to and after each sample block.

Strontium Isotope ( 87 Sr/ 86 Sr) Analysis of Fish Teeth
Sr isotope analyses of all teeth were performed either by solutionbased or LA-MC-ICP-MS (laser ablation inductively coupled plasma mass spectrometry) using a Nu Plasma MC-ICP-MS (Nu Instruments TM ) at the Max Planck Institute for Chemistry, Mainz, following the methods described by Weber et al. (2017). For solution-based analysis, the MC-ICP-MS was coupled to a CETAC Aridus 2 desolvator nebulizer system and samples were diluted to ≈ 100 ng/g Sr. A standard-sample bracketing approach was applied, using a static multi-collection mode with 100 cycles and 5 s integration time. All measurements were corrected using NIST SRM 987 ( 87 Sr/ 86 Sr of 0.710248; McArthur et al., 2020). Prior to analysis, Sr was separated by ion-exchange chromatography using Sr spec resin, following the protocol presented in Weber et al. (2018). For in situ laser ablation analysis, the MC-ICP-MS was coupled to a 213 nm Nd:YAG laser ablation system (New Wave Research TM UP-213). Signals of the following masses were monitored during the analysis: 82 Kr, 83 Kr, 84 Sr, 85 Rb, 86 Sr, 87 Sr, 88 Sr and half masses 85.5 and 86.5 (half masses were only monitored during LA-MC-ICP-MS). Mass bias correction was performed using the exponential law (Ingle et al., 2003) and assuming a constant 86 Sr/ 88 Sr of 0.1194. Laser ablation analysis was performed by linescan analysis, each linescan with a length of 1000 µm, a circular spot size of 65-100 µm (depending on the sample and daily tuning parameter), a transition rate of 5 µm/s and a laser repetition rate of 10 Hz. The resulting fluence varied between 10-30 J/cm 2 , depending on the chosen energy output (70-100%). Prior to each analysis, a pre-ablation was performed to prevent surface contamination.
To correct for Kr in the plasma, we used the "on peak zero" method by measuring a 45 s gas background prior to each analysis without the laser firing. After background subtraction, the remaining signal on mass 82 was used to check for Ca argide and dimer formation (<0.1 mV). Half masses were checked for doubly charged rare earth elements (Yb and Er) and the signals were corrected using constant isotope ratios (Berglund and Wieser, 2011). To correct for the occurrence of 87 Rb on mass 87, we used a 87 Rb/ 85 Rb of 0.3857 (Berglund and Wieser, 2011).

Phosphate Oxygen Isotope (δ 18 O PO4 ) Analysis of Fish Tooth Enameloid
The phosphate fraction of the fish tooth enameloid was separated using a method modified after Dettmann et al. (2001) and described in detail by Tütken et al. (2006). In summary, approximately 5 mg of pretreated sample powder was digested in 0.8 ml HF (2 M) on a vibrating table for ca. 12 h. After centrifugation (3 min at 10,000 rpm), the supernatant sample solution was separated from the CaF precipitates into a new vial. After neutralising the HF solution with NH 4 OH (25%) in the presence of bromothymol blue as pH indicator, Ag 3 PO 4 was precipitated by adding 0.8 ml of 2 M AgNO 3 . After settling of the Ag 3 PO 4 crystals, the samples were centrifuged and the supernatant solution containing excess AgNO 3 was removed by pipette. The Ag 3 PO 4 precipitate was then rinsed five times with Milli-Q water and dried overnight in an oven at 50 • C.
Ag 3 PO 4 aliquots of 0.5 mg were placed into silver capsules and analysed in triplicate by means of high temperature reduction using a Finnigan TC-EA coupled via a Conflo III to a Micromass 100 GC-IRMS at the University of Mainz, or to a Finnigan Delta Plus XL GC-IRMS at the Universities of Tübingen and Lausanne, according to the method of Vennemann et al. (2002). The raw δ 18 O PO4 values were normalised to an Ag 3 PO 4 reference material produced by Elemental Microanalysis with a certified value of 21.7 (silver phosphate P/N IVA33802207, batch no. 180097, distributed by IVA Analysentechnik, Germany). The analytical precision for this reference material was better than ± 0.3 (1 SD). For untreated NIST SRM 120c Florida phosphate rock, we obtained a δ 18 O PO4 value of 21.9 ± 0.3 (1 SD, n = 15). This value compares well with the values from most other laboratories as compiled in Chenery et al. (2010).

Bioapatite δ 18 O PO4 Value Reconstruction for Fish From Levantine Water Bodies
The expected isotope equilibrium range of δ 18 O PO4 values for fish tooth bioapatite was calculated for each of the water bodies in Egypt and Israel from which fish were sampled using published and a few measured δ 18 O H2O values ( Table 4). The calculations are based on the temperature-dependent isotope fractionation during biomineralisation of apatite by Lécuyer et al. (2013): T • C = 117.4 -4.5 (δ 18 O PO4 -δ 18 O Water ), where δ 18 O PO4 and δ 18 O Water correspond to the isotope compositions of bioapatite and ambient water relative to VSMOW, respectively. This relationship is valid for the temperature range of 8 • C < T < 32 • C thus covering water temperatures typically encountered in the Mediterranean region. For Mediterranean surface water, we used a δ 18 O Water of 1.6 ± 0.3 and an annual temperature range of 17-30 • C (Sisma-Ventura et al., 2014). The estimated δ 18 O PO4 of Mediterranean deep water was calculated for a temperature range of 13-15 • C and a δ 18 O Water value of 1.5 ± 0.2 . For the Nile delta water bodies, we used the post-Aswan dam δ 18 O Water range between 0.0 and 3.0 (Nile near Minia) and between 2.0 and 5.0 (Lake Manzala) and a water temperature range between 18 to 30 • C (Kolodny et al., 1983;Aly et al., 2004). For Lake Kinneret and its main tributary, the Jordan River, we used a δ 18 O Water range between −6.0 and 2.0 and between −8.0 and −5.0 and a temperature range between 16-28 • C and 10-22 • C, respectively (Stiller and Magaritz, 1974;Talbot, 1990;Zohary et al., 1994). We also estimated the δ 18 O PO4 equilibrium range of fish from past Nile delta water bodies during the Iron Age, the New Kingdom period, by using the Nile δ 18 O Water range of −2.0 to −1.0 , reconstructed from the δ 18 O PO4 values of tooth enamel from human mummies  and assuming a temperature range similar to today.

The 87 Sr/ 86 Sr and δ 18 O PO4 Expectation Fields for Levantine Fish Habitats
For each of the major water bodies in Egypt and the Southern Levant we also compiled a range of expected 87 Sr/ 86 Sr water Nile ( Kolodny and Raab (1988). This reference dataset was then used as an interpretive framework to assess the provenance of the archaeological S. aurata and sharks based on their enameloid Sr and O isotope composition using Linear Discriminant Analysis (LDA).

Linear Discriminant Analysis (LDA)
For assessing the fish provenance from different waterbodies, we applied a LDA, a commonly used supervised pattern datareduction and recognition method that assigns new observations to certain classes. This is achieved through a linear discriminant function (LDF) that maximises ratios of between-group variance and within-group variance. We used LDA (R package MASS version 7.3-51.6, Ripley, 2002) in order to assign provenance to archaeological fish tooth samples based on their enameloid Sr and O isotope composition. The 87 Sr/ 86 Sr ratios and δ 18 O PO4 values from published and measured water samples of known provenance and associated modern fish bioapatite (bone, dentine, enameloid) samples were used as predictors. Note that the expected fish bioapatite δ 18 O PO4 values for each waterbody were calculated from the δ 18 O H2O values and water temperatures using the phosphate-water oxygen isotope fractionation equation of Lécuyer et al. (2013). In order to compute the LDF, a total of 82 observations from six different water bodies in the Southern Levant (Bardawil Lagoon, Lake Manzala, Nile River, Lake Kinneret/Jordan River, and modern Mediterranean seawater) were used. These comprise 44 observations from water samples (using min, max and mean  Lécuyer et al. (2013). Note that modern Nile water is 18 O-enriched due to evaporation after building of the Aswan High Dam (Bialik and Sisma-Ventura, 2016) and thus has even slightly higher δ 18 O Water than the Mediterranean seawater. Therefore, calculated bioapatite δ 18 O PO4 values for fish from Nile and Mediterranean seawater overlap. However, Sr isotopes are clearly distinct for both settings. The enameloid δ 18 O PO4 values of fossil shark teeth from the City of David fall well into the expectation field for Late Cretaceous seawater based on Late Cretaceous shark teeth δ 18 O PO4 values of Kolodny and Raab (1988) and Pucéat et al. (2003). values taken from literature data and few measurements) and 38 observations from 14 modern fish specimens of known provenance (Table 1 and Supplementary Tables S4-S6). In order to increase the sample size of data used to compute the LDF, specimens that underwent several Sr isotope analyses per tooth but only a single O isotope analysis were used. Thus a single δ 18 O PO4 value was used for each 87 Sr/ 86 Sr measured in the same individual. Seven Bronze and Iron Age archaeological Gilthead seabream teeth were subsequently assigned a provenance to one of the six different water bodies using LDA. Provenance of the fossil, presumably Late Cretaceous sharks (n = 7;  (Kolodny and Raab, 1988;Kolodny and Luz, 1991) were added to the reference dataset used to compute the LDF to assess the robusticity of the model (Supplementary Tables S7-S9). All values were transformed by subtracting the mean of the predictor's data from the observation's value and then by dividing by the standard deviation.
Accuracy of the groups assigned by the LDF were assessed with a confusion matrix using a "leave-one-out" cross-validation method, where an observation is omitted and then assigned by the LDF. The procedure is then repeated for each observation and allows for the success-rate of classification to be calculated (Friedman et al., 2001). All statistical analyses were performed using the statistical program R (version 4.0.2; R Core Team, 2018).

Shark Taxonomy and Biostratigraphic Age
All of the shark teeth displayed a grey-beige colouration and macroscopic preservation typical for Cretaceous fossil shark teeth from the Near East (Supplementary Figure S1). Several shark teeth were identified to the order of Lamniformes and some were assigned to Squalicorax sp. (Anacoracidae; n = 2), a taxon that only occurred during the Late Cretaceous Cicimurri, 2005, 2006). Another tooth was assigned to the family Odontaspididae (Sand shark; n = 1), also from the Cretaceous (Figure 3B and Table 3). The fossil nature of the shark teeth was further tested by analysing different mineralogical and geochemical proxies for diagenesis.

X-Ray Diffraction Analysis
The differences in crystallinity of the analysed fish tooth dentine samples are quite significant, and any peak broadening contribution from the instrument and sample preparation are considered to be negligible. The X-ray diffraction pattern of the dentine from four shark teeth exhibited three well-separated diffraction peaks (211, 112, 300) typical of recrystallised bioapatite in fossil dentine or bone (Supplementary Figure S3). In contrast, the dentine of the Chalcolithic Gilthead seabream teeth still had one broad diffraction peak, similar to the unaltered dentine of a modern Great white shark (Supplementary Figure S3). Therefore, the Gilthead seabream have a lower apatite crystallinity than the shark tooth samples. Fluorapatite was identified as the primary mineral phase, indicating diagenetic fluorine uptake into the dentine that would have originally consisted of a fluor-bearing hydroxylapatite in teleost fishes (<0.3 wt.% F; Suga et al., 1993). Modern shark dentine can contain high fluorine content (0.6-1.5 wt.%; Enax et al., 2012) in its carbonated bioapatite but this content is still much lower than in enameloid that consists of fluorapatite (3.1 wt.%; Enax et al., 2012;Kocsis et al., 2015). Some dentine samples (e.g., SK-84) also contained small amounts of calcite as a secondary diagenetic mineral infilling.

Total Organic Carbon (TOC) Analysis
The dentine of the four shark teeth had low TOC contents of 0.23 ± 0.26 wt.% (n = 4) which is more than 40 times lower than the TOC content of dentine in the modern Great white shark (9.54 wt.%) and more than eight times lower than in the dentine of Chalcolithic Sparus aurata teeth (1.99 ± 0.47 wt.%; n = 3) (Supplementary Table S2 and Supplementary Figure S4).

Trace Element in situ LA-ICP-MS Analysis
Profiles of trace element content were measured across both dental tissues (enameloid and dentine) of all S. aurata (n = 5) and shark teeth (n = 7) from Chalcolithic and Iron Age cultural layers as well as one modern Great white shark and two S. aurata teeth as comparison to obtain unaltered dental tissue compositions (Supplementary Figure S5 and Supplementary Table S3). Here we focus on two trace elements that are highly sensitive to diagenetic alteration -uranium and neodymium -to identify fossil specimens. Uranium and Nd are both typically only present in very low, ng/g-level concentrations in modern mammalian tooth bioapatite (Kohn et al., 1999) and also in modern fish teeth (Shaw and Wasserburg, 1985;Vennemann et al., 2001;Kocsis et al., 2015). The analysed modern Gilthead seabream and Great white shark dentine yield low Nd and U contents of 0.001-0.006 µg/g and 0.05-0.09 µg/g, respectively ( Figure 5). This is in good agreement with the Nd content of 0.011 ± 0.011 µg/g in modern shark teeth (Shaw and Wasserburg, 1985;Vennemann et al., 2001). In contrast, the dentine of the shark teeth excavated from the City of David site had ≈ 100-1000 times higher U and Nd concentrations of ≈ 10-100 µg/g, typical for fossil dentine of shark teeth (Vennemann and Hegner, 1998;Kocsis et al., 2009). The dentine of the seabream teeth had intermediate Nd and U concentrations of ≈ 1-10 µg/g, and thus are ≈ 10-100 times enriched in U but less enriched in Nd as compared to the shark teeth (Figure 5). Enameloid is much less enriched (i.e., altered) in Nd and U compared to the dentine of the same tooth, both for Gilthead seabream and shark teeth.

Sr and O Isotope Data for Modern Fish Teeth and Bones
The 87 Sr/ 86 Sr and δ 18 O PO4 data for modern sharks and bony fish from different water bodies in Egypt and the Modern fish tooth dentine Nd and U field drawn after data from Vennemann et al. (2001) and this study. Fossil shark teeth have much higher dentine U and Nd contents than archaeological Gilthead seabream teeth. (B) Uranium (U) and strontium (Sr) concentration in fish teeth. Note the logarithmic scale of the y-axis. The grey Sr and U field for modern shark teeth is based on Sr concentrations from Vennemann et al. (2001) and Kocsis et al. (2015), U concentrations from this study. Note, despite U enrichment enameloid Sr concentrations of both archaeological Gilthead seabream and fossil shark teeth are still in the range of modern sharks. Table 2 and are presented with respect to the expected values for each water body (Figure 6). A typical modern marine 87 (Figure 6 and Table 4). The 87 Sr/ 86 Sr of extant freshwater fish from the Nile (Egypt) and Lake Kinneret (Israel), exhibited lower values between 0.707 and 0.708. The 87 Sr/ 86 Sr of extant Nilotic fish: Nile tilapia (Cichlidae; Oreochromis niloticus), upside-down catfish (Mochokidae: Synodontis schall) and bagrid catfish (Bagridae; Bagrus bajad) caught at different places in the Nile are very similar, yielding an average of 0.7071 ± 0.0002 (1 SD) and likewise for extant freshwater fish from Lake Kinneret (0.7073 ± 0.0003, 1 SD; Figure 6 and Table 4). This similarity results from the fact that the drainage area of both freshwater bodies is comprised of extended volcanic bedrocks with low 87 Sr/ 86 Sr (Krom et al., 1999;Hartman and Richards, 2014).

Southern Levant are summarised in
Relatively high δ 18 O PO4 values were observed in extant Nilotic fish, ranging between 20.1 and 24.5 , reflecting the Nile's source water which derives from monsoonal rains (δ 18 O Water range of −3.0 to 0.0 ; Schilman et al., 2001) with additional 18 Oenrichment due to evaporation effects along the Nile trajectory through the desert and in the Aswan High Dam .
Three species of extant freshwater cyprinids were analysed for their 87 Sr/ 86 Sr: Luciobarbus longiceps, Carasobarbus canis, and Capoeta damascina. The first two species have 87 Sr/ 86 Sr that vary between 0.70731 and 0.70757, agreeing well with water from different depths of Lake Kinneret which have an average 87 Sr/ 86 Sr of 0.70750 ± 0.00005 (2 SE) (Fruchter et al., 2017). Even lower values, between 0.70681 ± 0.00008 and 0.70729 ± 0.00006 (2 SE) were obtained for the C. damascina samples. Modern Lake Kinneret fish display a large δ 18 O PO4 range between 13.8 and 20.1 , reflecting the local hydrological conditions with an average rainwater δ 18 O Water value around −7.0 to −6.0 (Gat and Dansgaard, 1972) and the strong evaporation effects in the Kinneret basin that lead to an 18 O-enrichment in the  Table 4). Note that all modern fish teeth fall within the according expected isotopic ranges of the water bodies in which they were captured. The modern blue fish from the brackish Kishon River estuary falls within the calculated range of mixed Mediterranean seawater and Israeli groundwater (blue stippled lines), with Kishon River as the primary freshwater source. Note the fossil marine sharks fall within the expectation field for Late Cretaceous seawater and do not overlap with any modern marine influenced water body. lake water resulting in δ 18 O Water values between −2.5 and 0.0 (Figure 6; Gat, 1970).
Linear discriminant analysis was performed on water samples of known provenance and associated modern fish tooth samples for cross-validation, producing an error rate of 18%. This is due to the multiple partial overlaps in 87 Sr/ 86 Sr and δ 18 O PO4 values from the different bodies of water: Bardawil Lagoon with the Mediterranean Sea, and the Nile River with both Lake Manzala as well as Lake Kinneret/Jordan River (Figures 4, 6 and Table 4). Nevertheless, more than 80% of the modern fish of known provenance were correctly assigned to their respective water bodies using LDA (Supplementary Table S4).

Sr and O Isotope Data of Archaeological and Fossil Fish Teeth
The 87 Sr/ 86 Sr and δ 18 O PO4 data of sharks and S. aurata from Chalcolithic and Iron Age archaeological layers of Israel are summarised in Table 3 and shown in Figure 6. The shark teeth obtained from the Iron Age I-II periods of Jerusalem, all yielded unexpectedly low enameloid 87 Sr/ 86 Sr and δ 18 O PO4 values for contemporaneous sharks, ranging between 0.7074 and 0.7076, and between 18.3 and 20.7 , respectively. Placing the 87 Sr/ 86 Sr of the seven shark teeth on the Sr seawater curve (McArthur et al., 2020) yields Late Cretaceous Sr stratigraphic ages ranging from 86.5 to 76.3 Ma (Table 3) with an average late Santonian/early Campanian age of around 80.29 Ma +3.24, −2.64 Ma. The age uncertainties include analytical uncertainties, the 95% confidence interval of the Sr seawater curve, as well as the standard deviation of the calculated Sr stratigraphic ages. Six shark teeth cluster between 76.3 to 80.7 Ma (early Campanian). Only specimen SK-86 yields an older age of 86.5 Ma (latest Cenomanian/early Santonian) ( Table 3). Their δ 18 O PO4 are in good agreement with those of Campanian fossil fish teeth (teleosts and sharks) from Israel and Jordan ranging from 18.3 to 20.0 (Kolodny and Raab, 1988). In contrast, S. aurata teeth from similar archaeological contexts exhibited marine 87 Sr/ 86 Sr of around 0.7092, typical for modern seawater, and δ 18 O PO4 values that agree well with the range expected for fish from Mediterranean surface water or the hypersaline water of the Bardawil Lagoon (Figures 4A, 6), which formed as the result of seawater evaporation at the north Sinai coast.
To assign a provenance to archaeological fish samples, LDA was performed using the modern water and fish tooth reference dataset. All Chalcolithic/Iron Age seabream specimens were assigned to either the Bardawil Lagoon (50%, n = 3) or to the Mediterranean Sea (50%, n = 3) (Supplementary Tables S4-S6). Both Iron Age seabream teeth SP-11 and SP-13 from the City of David as well as one Chalcolithic specimen (EB-5) from Gilat (also with a high, Bardawil-like, δ 18 O PO4 value) derive from the hypersaline Bardawil Lagoon. The other three seabream specimens from Gilat were all assigned a Mediterranean provenance (Supplementary Table S5). For the fossil shark teeth, all specimens were assigned to the Late Cretaceous seawater field, except for sample SK-83 (Supplementary Tables S7-S9). The LDA gave the same result when Late Cretaceous seawater values were restricted to those of the Campanian period. In both scenarios the sample SK-83 was assigned to Lake Kinneret/Jordan River, which is obviously incorrect for an 80 million-year-old fossil shark, and thus is best explained by the overlap of the two water bodies in Sr and O isotope space. Note also the incorrect assignment of the Nile perch (Lates niloticus) to the Mediterranean (Supplementary Table S5), resulting from diagenetic alteration of its bone Sr isotope composition (Figure 7). FIGURE 7 | Strontium isotope ratios of pairs of enameloid and dentine of fossil shark and archaeological Gilthead seabream (Sparus aurata) teeth from Iron Age (Jerusalem, City of David) and Chalcolithic (Gilat) archaeological sites ( Table 3). A Nile perch vertebra bone from the Iron Age site of Ashkelon and one jawbone from S. aurata from the City of David are also plotted. Both are clearly diagenetically altered and fall within the range of values expected due to diagenetic alteration (grey field) defined by the range of 87 Sr/ 86 Sr for soil and rainwater compositions from Israel (Table 4). In contrast, none of the enameloid samples seem to be affected by diagenesis and instead preserve the original water 87 Sr/ 86 Sr of the Gilthead seabream and shark habitats. Note that all seven shark enameloid samples from different taxa have low 87

DISCUSSION
Using fish tooth enameloid 87 Sr/ 86 Sr and δ 18 O PO4 values collected from extant fish from different water bodies in Egypt and the Southern Levant enabled us to establish an isotopic reference dataset, which we then used as an interpretive framework for inferring the provenance of ancient S. aurata and different shark species recovered in coastal and inland archaeological layers from the Chalcolithic and Iron Age I-II periods (Figure 2A). Post mortem diagenetic alteration was carefully assessed to ensure that the original Sr and O isotope compositions are well preserved and still reflect the isotope values of the ambient water of the fish habitat. Fish taxonomic identification, combined with XRD, TOC and trace element (i.e., REE, U) analysis of dentine clearly demonstrated that the S. aurata teeth are sub-fossil archaeological fish teeth while the shark teeth from the same Iron Age cultural layers are all fossil teeth, most likely from the Late Cretaceous.

Variation of 87 Sr/ 86 Sr in Extant Shark and Bony Fish as Marker of Their Habitat
Marine shark and bony fish tooth enameloid both record the 87 Sr/ 86 Sr of ambient seawater. Hence, extant sharks from different oceanic basins and water depths with different hydrological conditions in temperature and salinity have rather uniform, seawater-like tooth 87 Sr/ 86 Sr values, averaging 0.709167 (range: 0.709160-0.709172, n = 8; Vennemann et al., 2001). In this study we obtained a similar value of 0.70919 ± 0.00005 (2 SE, n = 2; Table 2) for extant sharks from the Mediterranean (Centrophorus granulosus) and from the Red Sea (Carcharinus plumbeus), which matches the modern-day global seawater value of 0.70918 (McArthur et al., 2020). Similarly, tooth enameloid of modern S. aurata from the Eastern Mediterranean Sea also yields a typical marine values ( Table 2). Much lower 87 Sr/ 86 Sr, between 0.70681 ± 0.00008 and 0.70757 ± 0.00005 (2 SE), were measured in the teeth of three species of modern Cyprinidae from Lake Kinneret, which receives its water from the Golan Heights drainage area that is mainly underlain by basaltic bedrocks (Hartman and Richards, 2014). The two freshwater cyprinids Luciobarbus longiceps and Carasobarbus canis have 87 Sr/ 86 Sr that vary between 0.70731 and 0.70757, agreeing well with Lake Kinneret water with an average 87 Sr/ 86 Sr of 0.70750 ± 0.00005 (2 SE) (Fruchter et al., 2017). Even lower values, between 0.70681 ± 0.00008 and 0.70729 ± 0.00006 (2 SE), were obtained for the C. damascina samples that are in good agreement with the Sr isotope composition of the main fresh water source of Lake Kinneret and the Jordan River, which has an average 87 Sr/ 86 Sr of 0.70678 ± 0.00006 (1 SD) (Fruchter et al., 2017). This is well in line with the known migration behavior of some Cyprinidae species, including C. damascina, between Lake Kinneret and the Jordan River (Fishelson et al., 1996). All of these results further support the argument that fish, including sharks, record the 87 Sr/ 86 Sr of their ambient water. All extant marine and freshwater fish fall well within their respective expected isotopic ranges for the aquatic habitats in which they were caught (Figure 6). This was supported by Linear Discriminant (LDA) analysis, which correctly assigned more than 80% of the modern fish to the correct water bodies based on their tooth Sr and O isotope composition (Supplementary Table S4).
Nilotic fish teeth and bones have 87 Sr/ 86 Sr between 0.70681 and 0.70737 (Table 2 and Figure 6), which are in good agreement with modern Nile water values ranging from 0.70600 to 0.70718 (Brass, 1976;Gerstenberger et al., 1997; this study; Table 4). Similar values occur in freshwater fish from Lake Kinneret and the Jordan River. This similarity results from the fact that the drainage area for each of these water bodies is underlain by volcanic bedrocks. In this case, the coupling of 87 Sr/ 86 Sr and δ 18 O PO4 data from the same teeth can be used to further differentiate fish from different aquatic habitats with similar bedrock types (and hence dissolved Sr of similar 87 Sr/ 86 Sr) in their drainage area, but different hydrological regimes (i.e., water sources, salinities and/or degree of evaporation).

Assessing Post Mortem Diagenetic Alteration and Fish Provenance
Post mortem diagenetic alteration in the burial environment may have overprinted and shifted the original marine Sr isotope signature of fish teeth toward lower 87 Sr/ 86 Sr. We assess diagenetic alteration by comparing the 87 Sr/ 86 Sr of dentine, which is much more prone to alteration than enamel (oid) (Ayliffe et al., 1994;Hoppe et al., 2003;Tütken et al., 2008), to that of enameloid from the same tooth (Figure 7). Both dental tissues record in vivo the same 87 Sr/ 86 Sr of the ambient water ( Table 2; Vennemann et al., 2001;Tütken et al., 2011) but post mortem dentine is more prone to exchange of Sr with water in the soil (Tuross et al., 1989;Becker et al., 2008), similar to bone (Nelson et al., 1986;Hoppe et al., 2003). Thus, 87 Sr/ 86 Sr of dentine and bone are both expected to shift toward diagenetic soil water values more easily than enameloid. To infer Sr isotopic compositions of potential diagenetic fluids, rain and bedrock 87 Sr/ 86 Sr in Israel are considered (Figure 7 and Table 4).
The main soil types in the area where the fish teeth were excavated are Rendzina, which typically has 87 Sr/ 86 Sr between 0.7079 and 0.7084 and Terra rossa, with a mean value of 0.70857 ± 0.00026 (1 SD) (Hartman and Richards, 2014). Both soils developed by weathering over Cretaceous-Eocene marine sedimentary rocks, thus soil pore water values are expected to have similar 87 Sr/ 86 Sr (Table 4). Several studies have demonstrated that Sr from pore water in soils can exchange with bioapatite (Nelson et al., 1986;Hoppe et al., 2003). Thus, diagenetic alteration would be expected to shift 87 Sr/ 86 Sr toward soil water values. Indeed, the dentine of the Chalcolithic S. aurata teeth shows slight modification toward the soil water field, as do the fossil shark teeth from Jerusalem.
Furthermore, the 87 Sr/ 86 Sr measured in modern rainwater samples from Israel range between 0.7079 and 0.7092 (mean 0.7087 ± 0.0004, 1 SD; n = 18; Herut et al., 1993;Frumkin and Stein, 2004). The major sources contributing to the 87 Sr/ 86 Sr in rainfall are Cretaceous-Eocene marine sedimentary rocks, dissolved Mid-Cretaceous-Eocene chalk dust from the Sahara (0.7078) and small quantities of sea spray (0.7092) (Herut et al., 1993). However, none of these Sr sources has low enough 87 Sr/ 86 Sr to explain the values we obtained from the fossil shark teeth from the Iron Age cultural layers. Groundwater from Late Cretaceous aquifers of Israel have low 87 Sr/ 86 Sr of around 0.7075 but these deep aquifers are located further inland in the Judean and Samarian arches (Starinsky et al., 1980) and thus cannot have contributed Sr to cause diagenetic alteration in the archaeological sites investigated.
Thus, the low enameloid 87 Sr/ 86 Sr of the shark teeth cannot be explained by diagenesis as no potential Sr source with such low values is present in the burial setting. In contrast, diagenesis should increase the 87 Sr/ 86 Sr of the shark teeth, as occurs in the dentine of the same teeth, which are clearly altered (Figure 5) and shifted to higher 87 Sr/ 86 Sr (Figure 7). Even if some degree of diagenesis is invoked, the original shark enameloid 87 Sr/ 86 Sr are even lower, suggesting that measured data should be considered as maximum values. Regardless, any significant alteration is rather unlikely as enameloid of marine S. aurata still reflects expected seawater 87 Sr/ 86 Sr within an uncertainty of 0.00002 (2 SE) and is not shifted toward lower values of the diagenetic soil water endmember, as is the case for the dentine/bone of the same teeth/fish (Figure 7). Therefore, their original enameloid 87 Sr/ 86 Sr is very likely preserved.
Moreover, the dentine in the shark teeth has higher 87 Sr/ 86 Sr than the enameloid, approaching 87 Sr/ 86 Sr of bedrocks and rainfall in the study area (Figure 7). In contrast, dentine 87 Sr/ 86 Sr of S. aurata are shifted to lower values compared to the seawater-like enameloid values. Thus, dentine 87 Sr/ 86 Sr of both sharks and S. aurata converge on the diagenetic fluid values, which are controlled by the local bedrock substrate (Figure 7). The opposing directions of the isotopic shifts in 87 Sr/ 86 Sr between the enameloiddentine/jawbone pairs of shark and S. aurata teeth further support preservation of original enameloid 87 Sr/ 86 Sr in fish teeth. Therefore, we can conclude that original, low enameloid 87 Sr/ 86 Sr ≈ 0.7075 values in the shark teeth and high, modern seawater-like enameloid 87 Sr/ 86 Sr ≈ 0.7092 in S. aurata teeth from the same layers are still preserved. The shark enameloid 87 Sr/ 86 Sr fall well into the Late Cretaceous seawater field (Figure 6), as expected based on the taxonomic identification of some of the shark teeth (i.e., Squalicorax sp.) to Late Cretaceous taxa. Thus, these shark teeth are clearly fossils although they were found together with thousands of other archaeological fish remains in Iron Age cultural deposits ( Table 1).
Sparus aurata is a marine fish with a unique life cycle that includes exploitation of hypersaline lagoons (Tancioni et al., 2003). We observed elevated δ 18 O PO4 values in the teeth of S. aurata, implying that these fish originated from closed or semi-closed water bodies with a high degree of evaporation, such as the hypersaline Bardawil Lagoon situated on the Northern Sinai coast in Egypt (Kolodny et al., 1983;Sisma-Ventura et al., 2015, 2018. A high rate of evaporation increases the δ 18 O Water of the Mediterranean seawater from which the shallow, nearshore Bardawil Lagoon forms. It does not, however, change its water 87 Sr/ 86 Sr. This was confirmed by strontium isotope analysis of a water sample taken from the Bardawil Lagoon in January 2019 yielding a 87 Sr/ 86 Sr of 0.70916 ± 0.00002 (2 SE) which is the same (within error) as the value for hypersaline eastern Mediterranean seawater (0.709172; Reinhardt et al., 2001 and references therein; Table 4). One S. aurata tooth from Gilat (EB-5) and two from the City of David (SP-11, SP-13) are characterised by hypersaline δ 18 O PO4 (Figure 4A) and seawater-like 87 Sr/ 86 Sr (Figure 4 and Table 3). Based on their Sr and O isotope signature, the LDA assigned 50% (three out of six) of the archaeological S. aurata teeth to a hypersaline water origin, most likely the Bardawil Lagoon, and the remainder to the Mediterranean (Supplementary Table S5). This agrees well with previous δ 18 O PO4 results for S. aurata teeth from other Late Bronze Age to Iron Age sites in Israel and indicates a Bardawil Lagoon origin for these fish (Sisma-Ventura et al., 2018). Thus, seabreams were presumably exported from Northern Egypt (Bardawil Lagoon) to the Southern Levant (specifically to the City of David, Jerusalem) during the Iron Age, and perhaps even as early as the Chalcolithic. Since there are no indications for hypersaline lagoons (only for some brackish lagoons) along the Israeli coast during the Holocene (Sivan et al., 2011(Sivan et al., , 2016, the Bardawil and/or other hypersaline lagoons must have existed elsewhere in the Southern Levant. Finally, while the 87 Sr/ 86 Sr of the dentine of S. aurata is only somewhat altered by diagenesis, the 87 Sr/ 86 Sr (0.70892 ± 0.00006, 2 SE) of the Iron Age Nile perch vertebra from the coastal site of Ashkelon (Figure 7) is strongly altered and clearly shifted by >0.002 from expected Nilotic values ≈ 0.706 to 0.707 (Brass, 1976;Gerstenberger et al., 1997; Table 4) toward a much higher 87 Sr/ 86 Sr within in the soil water alteration field (Figure 7). This indicates that fish vertebrae, although abundant in archaeological sites, are prone to diagenetic alteration and thus not suitable to assess fish provenance.

Assessing the Origin and Age of Fossil Shark Teeth From Jerusalem
Some of the shark teeth recovered from strata dated to the 8-9 th century BCE (Iron Age I-II) (Reich et al., 2007(Reich et al., , 2008 in the City of David, Jerusalem, were identified (based on morphology) as extinct species from the Upper Cretaceous (e.g., Squalicorax). All shark teeth (including those which were not taxonomically identifiable), display trace element signatures of high U and Nd contents (Figure 5A), typical for fossilised fish teeth (Staudigel et al., 1985;Pucéat et al., 2003;Kocsis et al., 2009). This was further supported by the high apatite crystallinity (Supplementary Figure S3) and low organic carbon content of the dentine (Supplementary Figure S4) typical for fossil specimens but different from archaeological S. aurata teeth from the same cultural layers (thus, providing a special case for which the 87 Sr/ 86 Sr can be used as a tool to date marine fossils deposited in much younger archaeological contexts). This could improve our understanding regarding the extent and the meaning of this phenomenon in archaeological sites.
The 87 Sr/ 86 Sr values of these fossil shark teeth do not fit the modern reference dataset compiled for contemporaneous fish in the southern Levant. The mean enameloid 87 Sr/ 86 Sr value of 0.70751 ± 0.00006 (1 SD; n = 7) of the fossil shark teeth from Jerusalem is much lower than the 87 Sr/ 86 Sr measured in modern sharks (≈ 0.70917, Table 2 and Figure 6), but agrees well with Late Cretaceous seawater 87 (Kolodny and Raab, 1988). Assuming a seawater δ 18 O of −1 for a continental ice-free Earth, marine temperatures in the subtropical zone (i.e., 30-35 • N) were estimated to around 28-29 • C in the Cenomanian-Turonian period (Pucéat et al., 2003). The low shark enameloid δ 18 O PO4 values thus reflect tooth formation in the warm shallow, tropical seawater of the Tethys Ocean shelf (Kolodny and Raab, 1988).
The recovery of a relatively large number of shark teeth in the Iron Age assemblage of the City of David together with diverse fish remains is puzzling for two main reasons: (1) shark teeth are rare in archaeological sites in the Southern Levant (Table 1); (2) this study clearly demonstrates that these shark teeth are fossils and hence do not represent shark consumption. Instead the presence of these teeth may reflect collection for other cultural purposes (tools, pendants, etc.), as has been observed in different contexts worldwide (Cione and Bonomo, 2003;Betts et al., 2012;Charpentier et al., 2020). Alternatively, the fossils were simply a component of the local soil substrate and derive from weathering processes of the Late Cretaceous limestone substrate. However, this is not very likely as the Late Cretaceous marine sedimentary strata beneath the City of David are stratigraphically older (Turonian; Figure 8) than the Sr-stratigraphic age of the fossil shark teeth (late Santonian/early Campanian). However, fossil shark teeth from the Late Cretaceous are abundant in the marine sediments of the Mount Scopus group in Israel and in the Menuha Formation of the Southern Negev (Figure 2A; Gvirtzman, 2004;Gvirtzman et al., 2008;Wilson et al., 2012;Retzler et al., 2013).
Shark teeth recovered in the City of David were mainly from the "Rock Cut Pool" (Figure 2B; note that the complete area of this structure was fully excavated and sieved, Reich et al., 2007Reich et al., , 2008, and from area G, a mixed dump deposition, dated to ca. 586 BCE (Lernau, 2015). In area G, nine shark teeth were recovered, but these were not used in this study. However, except for these two locations, no fossil shark teeth have been reported from geological strata in the city of David or anywhere else in Jerusalem. Thus, the fossil shark teeth found in the City of David are probably not derived from the underlying marine limestone and dolomite of the Shivta Formation, which is of Turonian age (Rosenfeld and Hirsch, 2005; Figure 8), and is thus older than the Santonian/Campanian age inferred by Sr stratigraphy for the FIGURE 8 | Geological map of Jerusalem with the City of David modified from Rosensaft and Sneh (2011). Therefore, rock formation names are provided in both English and Hebrew (as provided in the official geological map of Jerusalem). The City of David is situated on marine sedimentary strata of the Shivta Formation being Turonian in age. Note that the fossil shark teeth found in the Iron Age strata yielded younger late Santonian/early Campanian Sr stratigraphic ages. Such fossil shark teeth are only known from the Menuha Formation of the Mount Scopus group in the Negev Desert (Kolodny and Raab, 1988;Retzler et al., 2013). These shark teeth bearing outcrops are 90 km up to 180 km south of Jerusalem in the Negev Desert ( Figure 2A). Note that sediments of the Menuha Formation also crop out at the hillslopes a few hundred meters east/southeast of the City of David, however, no fossil shark teeth were reported from these strata.
fossil shark teeth. However, Turonian (0.70732 ± 0.00002 1 SD) and Campanian (0.70758 ± 0.00008 1 SD) seawater 87 Sr/ 86 Sr differ on average by only 0.00026 (McArthur et al., 2020). Therefore, only a slight diagenetic alteration of the enameloid 87 Sr/ 86 Sr by +0.00016 would be sufficient to shift the value of Turonian shark teeth to the observed value 0.70751 and hence a late Santonian/early Campanian age. The dentine of two of the shark teeth was diagentically shifted by 0.0001-0.0002, yielding higher 87 Sr/ 86 Sr values compared to enameloid. Strontium isotopes in sub-mm sized Cenozoic fossil fish teeth from marine sediment drill cores are prone to diagenetic alteration (Martin and Scher, 2004). However, the fluorapatite of the enameloid of larger cm-sized shark teeth is less prone to diagenesis and still has low U and Nd contents (about 10 times lower compared to the dentine; Figure 5A) and no correlation was observed between Sr and U contents (Figure 5B), which could be taken as an indication against significant diagenetic alteration of the Sr isotope composition. The Sr stratigraphic age of the fossil shark teeth from the City of David fits perfectly to the Santonian to early Campanian age of the fossil shark teeth-bearing Menuha Formation in the southern Negev (Retzler et al., 2013). Therefore, we can conclude that these fossil shark teeth may have been collected in the Negev Desert up to 90 km south of Jerusalem (Figure 2A). For instance from the Letaot Mountain fossil shark teeth of late Santonian/early Campanian age are known (Kolodny and Raab, 1988). The Menuha Formation also crops out in hill slopes a few hundred meters to the east/southeast of the City of David (Figure 8), however, no finds of fossil shark teeth are reported from this or other areas in Jerusalem thus far. Therefore, it is entirely feasible that the fossil teeth were transported to the City of David although the reason for bringing these teeth to the city and their function remains elusive.

CONCLUSION
We establish a Sr and O isotope reference dataset for bioapatite of modern marine and freshwater fish from the Mediterranean, Bardawil Lagoon, Lake Kinneret, Jordan River and Nile River. Fish from these different water bodies in Egypt and the Southern Levant have distinct Sr and O isotope signatures that enable us to distinguish their provenance from Nilotic, hypersaline, marine and fresh water bodies using LDA. This will be an important asset for future provenance studies of archaeological fish remains in the Levant and will shed new light on the Levantine fish provenance and trade. Unexpectedly, multiple Late Cretaceous (80.3 ± 3.2 Ma) fossil shark teeth were encountered in the same Iron Age cultural layers of the City of David, Jerusalem together with a wide diversity of archaeological fish originating from the Mediterranean Sea and the Nile. These fossil shark teeth clearly do not represent food remnants although their use remains elusive. It may be they were brought to the city intentionally as the nearest fossil shark teeth bearing outcrops of Campanian age are situated in the Negev Desert (≈ 90 km south of Jerusalem), although similar aged strata also crop out near the City of David. Thus, in archaeological fish assemblages, fossil fish teeth can also be encountered. This should be carefully evaluated, and we here provide a toolset comprised of combined dentine XRD, TOC and trace element (i.e., U, REE) analysis to distinguish fossil from archaeological teeth based on degree of diagenetic alteration. In contrast, the enameloid of the fossil shark and the archaeological Gilthead seabream (S. aurata) teeth still preserves its original 87 Sr/ 86 Sr and δ 18 O PO4 values. The Sr and O isotope signatures indicate that some of the Gilthead seabreams were not caught locally in the Mediterranean, but instead in the hypersaline water of the Bardawil Lagoon and were imported from Egypt to coastal and inland sites in Israel during the Iron Age and probably already since the Chalcolithic. We demonstrate that this dual-isotope approach, frequently used to assess provenance of terrestrial vertebrates, is also a powerful toolset for tracing the provenance of fish and provides new insights into past human fish exploitation (even of fossil fish teeth) and trade patterns.

DATA AVAILABILITY STATEMENT
The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author/s.

ETHICS STATEMENT
Ethical review and approval was not required for the animal study because we analysed modern staple food fish obtained from fishermen as well as archaeological and fossil fish teeth.

AUTHOR CONTRIBUTIONS
TT and GS-V designed the research. TT, GS-V, IZ, and MW wrote the manuscript. GS-V and TT interpreted the isotope data. HH, IZ, and OL provided and prepared the fish samples. IZ determined the fish species and provided the all fish related information. MW performed the Sr isotope and LA-ICP-MS analysis and evaluated the Sr isotope data. TT evaluated TOC, XRD and trace element data. NB performed the statistics using Linear Discriminant Analysis (LDA). GS-V prepared the tooth samples for oxygen isotope analysis. KJ provided the analytical facilities for Sr isotope analysis. All authors contributed to the article and approved the submitted version.

FUNDING
This project received funding from the European Research Council (ERC) under the European Union's Horizon 2020 Research and Innovation Programme (grant agreement no. 681450). MW is thankful to the Max Planck Graduate Center for funding.

ACKNOWLEDGMENTS
We thank Torsten Vennemann and Michael Maus for performing phosphate oxygen isotope measurements and Hubert Vonhof for running water samples for oxygen isotope analyses and the dentine for carbonate content. We thank B. Stoll and U. Weis for assistance in the LA-MC-ICP-MS laboratory. Modern Mediterranean fish were collected as part of the Israel Oceanographic and Limnological Research deep-sea national monitoring program. We thank Sabine Fiedler and Björn Glasner for total carbon analysis, Sven Brömme for TIC analysis, and Ralf Meffert for the XRD analysis of the dentine samples. Special thanks to Jürgen Kriwet for taxonomic shark tooth identification. Or Bilak is kindly acknowledged for discussing Cretaceous geology of Israel and fossil shark tooth occurrence. Menachem Goren, from the Steinhardt Museum, Tel Aviv University, provided us information regarding the biogeographic distribution of Oreochromis niloticus. We also thank especially Jennifer Leichliter for meticulously proofreading the English.
Finally, yet importantly, we want to thank the two reviewers who helped with their constructive comments to significantly improve the manuscript. Especially the reviewer who pointed us to the possibility that the shark teeth from the archaeological layers could be fossil teeth, which was a game changer. It led us to perform additional analysis and, ultimately, to rework the narrative of the manuscript, changing the story from shark exportation from the Nile delta to fossil shark finds in cultural layers. This was a textbook example how a good review should work and we are very appreciative.