Isotope analyses to explore diet and mobility in a medieval Muslim population at Tauste (NE Spain)

Abstract

The Islamic necropolis discovered in Tauste (Zaragoza, Spain) is the only evidence that a large Muslim community lived in the area between the 8th and 10th centuries. A multi-isotope approach has been used to investigate the mobility and diet of this medieval Muslim population living in a shifting frontier region. Thirty-one individuals were analyzed to determine δ¹⁵N, δ¹³C, δ¹⁸O and ⁸⁷Sr/⁸⁶Sr composition. A combination of strontium and oxygen isotope analysis indicated that most individuals were of local origin although three females and two males were non-local. The non-local males would be from a warmer zone whereas two of the females would be from a more mountainous geographical region and the third from a geologically-different area. The extremely high δ¹⁵N baseline at Tauste was due to bedrock composition (gypsum and salt). High individual δ¹⁵N values were related to the manuring effect and consumption of fish. Adult males were the most privileged members of society in the medieval Muslim world and, as isotope data reflected, consumed more animal proteins than females and young males.

About the Authors:

Iranzu Guede

Affiliation: Department of Mineralogy and Petrology, Faculty of Science and Technology, University of the Basque Country-UPV/EHU, Vizcaya, Spain

Luis Angel Ortega

* E-mail: [email protected]

Affiliation: Department of Mineralogy and Petrology, Faculty of Science and Technology, University of the Basque Country-UPV/EHU, Vizcaya, Spain

ORCID http://orcid.org/0000-0002-8929-0328

Maria Cruz Zuluaga

Affiliation: Department of Mineralogy and Petrology, Faculty of Science and Technology, University of the Basque Country-UPV/EHU, Vizcaya, Spain

Ainhoa Alonso-Olazabal

Affiliation: Department of Mineralogy and Petrology, Faculty of Science and Technology, University of the Basque Country-UPV/EHU, Vizcaya, Spain

Xabier Murelaga

Affiliation: Department of Stratigraphy and Palaeontology, Faculty of Science and Technology, University of the Basque Country-UPV/EHU, Vizcaya, Spain

Miriam Pina

Affiliation: “El Patiaz” Cultural Association, Cuesta de la Cámara 12, Tauste, Zaragoza, Spain

Francisco Javier Gutierrez

Affiliation: “El Patiaz” Cultural Association, Cuesta de la Cámara 12, Tauste, Zaragoza, Spain

Paola Iacumin

Affiliation: Department of Physics and Earth Sciences, University of Parma, Parma, Italy

Introduction

Muslims invaded most of the Iberian Peninsula in the Early Middle Ages (AD 711) and remained for the next seven centuries, until 1492 when the Christian Kingdoms totally reconquered the peninsula. The northern frontier of the country captured by the Muslims, known as al-Andalus, extended eastward on the southern slopes of the Cantabrian range from the present Galicia to Catalonia. Following the Muslim conquest, al-Andalus was at first (711–750) a province of the Umayyad Caliphate centered on Damascus. From 740 a series of civil wars between various Muslim groups resulted in the breakdown of the Arab empire and the Emirate of Cordova (c. 750–929) emerged. In 929 the emir of Cordova proclaimed himself Caliph and the period of the Caliphate of Cordova was established (929–1031). The Cordova Caliphate collapsed during a civil war and Al-Andalus broke up into a number of mostly independent states called taifa kingdoms. The independent taifas were too weakened to defend themselves against the Christian Kingdoms in the north and west, allowing the Reconquest. The Christian reconquest of Iberia ended with the final assault on the Emirate of Granada in 1492. From 711 to 1492, as political dominions changed, the boundaries between the Christian north and the Islamic south shifted constantly.

In the Ebro Valley, the first Muslims arrived in the early 8th century conquering the main towns without any relevant or attested resistance and the Upper March (or northern frontier) was established along the Ebro basin. Thus began the Muslim period in the Ebro Valley that, for four centuries, was centered on the metropolis of Saragossa. Shortly after the Muslim conquest, the nobleman Count Cassius converted to Islam to preserve his lands and political power and founded the Banu Qasi dynasty. In the 9th century the Upper March was under the dominion of the Banu Qasi dynasty as a semi-autonomous territory within the Cordova caliphate [1,2]. During the 9th century, the Banu Qasi lineage was successively loyal and rebellious toward the Cordova emir. In the second half of the century, the Banu Qasi domains increased considerably, extending north to the Pyrenees and east nearly to the Mediterranean coastline. However, in the later 9th century the Cordova emir recovered most of the Upper March territories and in the early 10th century, harassed by its Christian neighbors and without the support of Cordova, the Banu Qasi dynasty lost all its territories [3].

The society of Al-Andalus was made up of three main religious groups: Christians, Muslims and Jews, who inhabited distinct neighborhoods in the cities. Islamic society stratification was mainly by ethnic division. The kinship system ascribed importance only to relationships through males and endogamous marriages were viewed as the ideal system [4]. The more powerful a tribal group was, the more women it would attract from outside and the fewer it would lose, and the more endogamous it would become. Under Islamic law, the most privileged members of society were devout Muslim men, and women were treated as second-class citizens [4]. In particular, women’s rights were contingent on their place within society on several levels, including their religious, economic and marital status. Under Islamic law, other groups in society such as Jews and Christians had fewer rights and privileges, to varying degrees.

Within this framework, Tauste was placed midway between the two most significant cities: Saragossa, metropolis of the Upper March, and Tudela, the center of the Banu Qasi territory.

The Muslim occupation of Tauste (Zaragoza, Spain) has been considered incidental and even non-existent, according to traditional and written sources. However, recent excavations suggest a large stable Muslim population lived in the town from the early Islamic period in the Iberian Peninsula. In 2010, a cemetery with several human skeletons aligned perpendicular to Mecca was discovered. The bodies were placed on their right side, facing towards Mecca, as is characteristic of a Muslim cemetery [5]. In contrast, no remains of the Muslim village associated with this necropolis have yet been found. Multi-isotopic studies, including radiogenic strontium, stable oxygen, carbon, and nitrogen, have been used to reconstruct the geographic origin, mobility and dietary practices of the Tauste individuals during the Islamic period in the Iberian Peninsula. Stable isotope composition of bone collagen reveals information about nutrition, life history, and mobility in past populations [6–10].

Isotope analyses in bioarchaeology

The analysis of carbon and nitrogen isotope composition in bone collagen constitutes an approach to palaeodietary reconstruction. It can provide information about the protein portion of the diet averaged over roughly the last 10 years prior to death and also about different protein sources [11,12].

Carbon isotope analysis provides information about the ecosystem that foodstuffs come from, distinguishing between terrestrial and marine ecosystems. In the case of a terrestrial diet, it informs about the plants that were consumed. Two classes of plants are distinguished according to their photosynthetic pathways: C₃ plants and C₄ plants. C₃ plants are most vegetables, wheat (Tritium) and barley (Hordeum vulgare), while C₄ plants include millet (Pennisetum), maize (Zea mays) and sugar cane (Saccharum officinarum). C₄ plants exhibit more enriched carbon values than C₃ plants, so that the mean δ¹³C values are -13‰ and -27‰ respectively [13,14]. Marine plants are all C₃ plants and their average values are about 7.5‰ higher than terrestrial C₃ plants. Carbon isotope composition can be used to distinguish marine protein consumption in terrestrial C₃-based diets, but when C₄ plants are involved marine and terrestrial values can overlap [15,16]. Carbon fractionates in δ¹³C by only about 1‰ throughout the food chain [6–8,17]. In freshwater ecosystems the δ¹³C composition of plants is variable and consequently freshwater fish exhibit a broad range of δ¹³C values that are largely depleted [18,19]. Therefore, δ¹³C ratios more negative than -22‰, the value corresponding to the low-end of a diet based only on C₃ terrestrial plants, suggest freshwater fish consumption.

Nitrogen isotope values reflect the intake of animal proteins and inform about the trophic level of an individual [10,20,21]. Thus, nitrogen isotopes in terrestrial ecosystems are enriched in δ¹⁵N by 2–5‰ (on average, 3‰) from food to body tissue as trophic levels increase [8,22,23] Terrestrial protein sources have δ¹⁵N values ranging from 5‰ to 12‰, while aquatic food sources range from about 12‰ to 22‰ for marine fish and 7.2‰ to 16.7‰ for freshwater fish [24–28]. When C₃ plants are consumed, nitrogen isotope analysis is combined with carbon isotope analysis to distinguish between proteins derived from terrestrial, freshwater and marine resources. Other reasons for variability in δ¹⁵N ratios of plants and animals include natural environmental conditions such as salinity and aridity or anthropogenic factors like manuring [29,30]. In general, human diet corresponds to a mixture of food with different isotope signatures. Plots of collagen δ¹³C vs δ¹⁵N values can be interpreted as mixtures of multiple components [31,32] that do not yield unique solutions, but may outline the dominant components in the diet of the studied individuals. Besides, stable nitrogen isotope analysis can also be used to investigate breastfeeding and weaning practices. In fact, during breastfeeding, children exhibit δ¹⁵N values enriched about 2–3‰ over that of their mothers [33]. Strontium and oxygen isotopes are two independent isotopic systems in which strontium reflects local geology and oxygen reflects geography and can be used to reconstruct movements of past populations. The combination of these two isotopic systems is able to constrain possible areas and provide information about an individual’s area of origin and thus determine mobility patterns [34–36].

Oxygen and strontium are fixed in phosphate in teeth and bones through ingested food and water. Strontium isotopes appear by substituting calcium in biogenic phosphate [37–39]. After formation during infancy, tooth enamel does not incorporate other elements and thus will reflect the geological composition of the place of residence during childhood, assuming that childhood residence and food production area coincide, at least for the majority of the food intake [37,40]. However, these patterns are not perfectly predictable at any level, because of vagaries in available food over time, and because the strontium ratio reflects an average value that synthesizes the geological composition of the different food provenances ingested during childhood. The average expected patterns are used to predict the most likely geographic links between tissue and location. In contrast, bones are continuously remodelled throughout an individual’s lifetime.

The radiogenic strontium isotopes are related to geology and vary according to the composition and age of bedrock. The strontium concentration in organisms varies according to the trophic level but the ⁸⁷Sr/⁸⁶Sr isotope signature of humans and fauna has negligible metabolic fractionation and will reflect the isotope signature of the underlying bedrock [39,41–48]. ⁸⁷Sr/⁸⁶Sr ratios in bedrock, soils, water and plants will be reflected in humans and animals that consume food and water from those sources [21,38,49]. Since the ⁸⁷Sr/⁸⁶Sr isotope ratio is inherited from the local environment, it is necessary to define the local bioavailable strontium isotope signature to evaluate residential mobility of individuals. There are several methods to establish the local baseline of the isotope signature by analyzing environmental samples including freshwater, soil leachates, ancient fauna and present-day small wild animals [38,50–52]. However anthropogenic activities such as the use of fertilizers could modify the strontium isotope ratios of modern ecosystems [53–56].

In contrast, the oxygen isotope reflects the isotopic composition of ingested water that is derived from meteoric water. The δ¹⁸O in precipitation varies regionally according to temperature and other climatic parameters, such as distance from the coastline, altitude and latitude [57–60]. Oxygen isotopes in the body are subject to several steps of metabolic fractionation. The fractionation mechanisms are relatively well known, allowing the calculation of approximate drinking water (δ¹⁸O_w) values from the δ¹⁸O_p of biogenic phosphate by means of conversion equations [57,60–66]. Despite complexities in the calculation of meteoric water isotope composition in the past, the oxygen isotope composition of human remains allows the identification of palaeomobility patterns.

The aim of this study was to reconstruct palaeomobility and palaeodiet patterns in the medieval Muslim population at Tauste. Tauste Muslim necropolis constitutes a suitable site to examine human mobility since it was located on the northern frontier during a very convulsive period of time. In addition, the palaoedietary pattern can illustrate the basic dynamics of medieval Muslim social life. For these purposes, stable isotopes (δ¹³C, δ¹⁵N, δ¹⁸O) and radiogenic strontium (⁸⁷Sr/⁸⁶Sr) were investigated to obtain information about nutrition and social stratification.

Archaeological setting

Tauste archaeological site is located in the town with the same name in the province of Zaragoza (northern Spain) (Fig 1). Tauste is in the Ebro basin, on the River Arba, a tributary of the River Ebro. The Muslim archaeological site of Tauste is formed only by the cemetery, with a total absence of other vestiges of Islamic population. All the graves were aligned SW-NE and the human bodies were carefully placed on their right side, facing Mecca, indicative of a Muslim necropolis [67] (Fig 2). All individuals were found in anatomical connection. Graves were dug in clay soil without any structure, or only a minimum structure formed by rammed earth on the sides according to Muslim burial rituals. More complex tombs corresponding to double grave burials (shaq or ladj) were found. A similar burial system has been documented in other Muslim necropolises in the Iberian Peninsula, such as Marroquíes Bajos (Jaen) [68], Tossal de Manises (Alicante) [69] or the recent find at Valdeherrera (Calatayud). The human remains extended over an area of two hectares and the density of graves (0.25–0.30 individuals/m²) indicates a minimum of 4.500 burials, excluding children [70,71]. Only a simple bronze hoop earring was found in a female’s grave, and the lack of grave goods is also indicative of Islamic funeral rituals. The excavations have found at least two levels of burials, indicating this cemetery was in use during an extended period of time.

[Figure omitted. See PDF.]

Fig 1. Location of Tauste archaeological site and the excavated areas.

Reprinted from under CC by license, with permission from [Instituto Geográfico Nacional (IGN)], original copyright [2015].

https://doi.org/10.1371/journal.pone.0176572.g001

[Figure omitted. See PDF.]

Fig 2. Aerial view of some burials showing individuals placed in the graves following the Muslim burial rituals (facing east).

https://doi.org/10.1371/journal.pone.0176572.g002

Radiocarbon dating of human bones dates the graveyard in the 8th to 10th centuries and it could be one of the oldest Muslim necropolises in the Iberian Peninsula (Table 1, Fig 3). Calendar ages were determined using the Oxcal v 4.2.4 program [72] with the latest IntCal13 calibration curve for atmospheric data [73]. Calibrated age ranges correspond to 95.4% probability (2σ) and are expressed in years cal AD. The age and extent of the necropolis suggest Tauste was a thriving village in the times of the Banu Qasi dynasty, when the northernmost limit of Al-Andalus was established [74,75].

[Figure omitted. See PDF.]

Fig 3. Radiocarbon dating of human bone samples from Tauste calibrated with OxCal v4.2.4 [72] and IntCal13 atmospheric data [73].

https://doi.org/10.1371/journal.pone.0176572.g003

[Figure omitted. See PDF.]

Table 1. Calibrated radiocarbon dating of the Tauste site.

https://doi.org/10.1371/journal.pone.0176572.t001

Geologically, the Muslim necropolis of Tauste is located in the Ebro Basin, composed by Tertiary (Miocene) and Quaternary sedimentary rocks of continental origin (Fig 4) [76]. Miocene materials around Tauste are composed by claystones with interbedded gypsum layers. These materials correspond to evaporite lacustrine facies, i.e. sediments deposited in the centre of a continental sedimentary basin. Miocene deposits are overlaid by Quaternary materials consisting mainly of river terrace deposits and fluvial sediments. The evaporitic nature of the bedrock causes a large increment in salt contents in the environment. In fact, high levels of sodium chloride and sulphate ions have been found in the freshwater River Arba [77,78].

[Figure omitted. See PDF.]

Fig 4. Geological map of the Tauste region showing evaporitic nature of the bedrock.

Reprinted from [76] under CC by license, with permission from [Instituto Geológico y Minero de España (IGME)], original copyright [2015].

https://doi.org/10.1371/journal.pone.0176572.g004

Materials and methods

This study deals with archaeological skeletal material and all necessary permits were obtained for the described study, which complied with all relevant regulations. The excavation licenses were issued by the General Director of Cultural Heritage of the Gobierno de Aragon (Spain) and are stored in its archives. Following excavation campaigns 2010/2013, the bones and teeth samples were transferred to the Heritage and Cultural Landscape Research Group (GIPyPAC) at the University of Basque Country-UPV/EHU, Spain for investigation. At present all archaeological remains, including the human bones, recovered at the site of Tauste are stored in the Museo de Zaragoza.

Carbon and nitrogen isotope measurements have been performed in bone collagen extracted from 31 human individuals corresponding to the sectors excavated in 2012 and 2013 (Fig 5) and nine faunal bone samples. Additionally, 23 teeth and 8 bone samples were analyzed for strontium and oxygen isotope studies. In order to define the strontium isotope baseline, four soil samples and one freshwater sample were analyzed. The soil samples were collected in different parts of the cemetery, while the freshwater sample was collected from the River Arba near Tauste. Surface waters were collected from the banks of the river. Before analysis, water samples were filtered to remove suspended particles.

[Figure omitted. See PDF.]

Fig 5. Detailed map of two excavated areas (sectors 2012 and 2013) showing studied individuals.

https://doi.org/10.1371/journal.pone.0176572.g005

Water aliquot was filtered with disposable syringe filters (0.45 um) and 10 mL aliquot was transferred to a 15 mL Teflon (Savillex™) vial, evaporated down on a hot plate at 80°C overnight and then dissolved in 1.5 mL of 2N HNO₃.

The measured individuals were corresponded to 12 males, 10 females and 9 of indeterminate sex [71]. The individuals were categorized by age into infants (0–3 years), children (6–12 years), juveniles (12–17 years), young adults (18–34 years), middle-aged adults (35–50 years), and older adults (older than 50 years). According to these categories, 2 individuals are older adults, 11 are middle-aged adults, 6 are juveniles, 7 are young adults, 4 are infants and 1 is of indeterminate age. Sex determination was carried out according to the classical patterns of dimorphism and age was defined by the most reliable markers: changes in auricular surface and pubic symphysis, epiphyseal closure, cranial sutures and dental eruption [79]. The faunal samples corresponded to three wood mice (Apodemus sylvaticus), three shrews (Crocidura russula) and two common barbels (Barbus barbus) and a madrilla (Parachondrostoma miegii).

For carbon and nitrogen isotope analyses, bone collagen was extracted following the procedure in Bocherens et al. [80]. 300 mg of bone sample powder were demineralised in 1M HCl for 20 min at room temperature until the sample dissolved. To remove humic acid the samples were rinsed with distilled water and treated with 0.125 M NaOH. The resulting insoluble fraction after being rinsed again with distilled water was gelatinized in HCl solution at pH3 for 17 h at 90°C. Then, samples were filtered with disposable syringe filters (5 um), freeze-dried and finally lyophilized. Lyophilized collagens (2.5–3.5 mg) were enclosed in tin capsules for isotopic analysis. Carbon and nitrogen isotope analyses were performed using an elemental analyzer on line with a continuous-flow isotope ratio mass spectrometer (EA-IRMS) at Iso-Analytical (Cheshire, UK). Replicate measurements of the liver standard NBS-1577B and ammonium sulphate IA-R045 working standard were run to confirm instrument accuracy. Replicate analysis of the NBS-1577B δ¹³C standard during runs gave a ¹³C/¹²C of −21.62 ± 0.02 (1σ, n = 7) and ¹⁵N/¹⁴N of 7.62 ± 0.13 (1σ, n = 7), and the IA-R045 working standard during runs gave a ¹⁵N/¹⁴N of −4.56 ± 0.17 (1σ, n = 4) and ¹³C/¹²C of −26.1 ± 0.03 (1σ, n = 4). Isotopic values are reported as δ values in per thousand (‰) relative to international defined standards for carbon (VPDB: Vienna Pee Dee Belemnite) and nitrogen (AIR: Ambient Inhalable Reservoir). The instrumental precision for δ¹³C was ± 0.06‰ or better and for δ¹⁵N was between ± 0.06‰ and ± 0.08‰, determined by replicated analyses of internal standards.

Tooth enamel and bones were used to determine strontium and oxygen isotope composition. The samples were washed in an ultrasonic bath to remove impurities and further cleaned by mechanical abrasion to remove the outer surface and avoid potential contamination.

For strontium isotope analysis a fraction of dental enamel was collected mechanically with a diamond-coated trepanation drill (MF-perfect, W & H Dentalwork, Bürmoos, Austria). The enamel sample was taken transversally. Enamel and bone samples (~10 mg) were dissolved in 7 mL Savillex® vials (Minnetonka, MN, USA) with 1.5 mL of 2N HNO₃ (analytical grade purified by subboiling distillation). In order to establish the local isotopic composition, two water samples and four soil samples were also analyzed. 15 mL of freshwater was evaporated to dryness and then dissolved in 2 mL HNO₃. A 1g aliquot soil sample was leached by adding 2.5 mL 1 M ammonium nitrate (NH₄NO₃) and shaking for 8 h to obtained the bioavailable Sr. After samples were centrifuged at 3000 rpm for 15 min, the supernatant was extracted (~1–2 mL) and evaporated to dryness and then redissolved in 2 mL HNO₃. The solutions were loaded into cation exchange columns filled with Sr.spec® (ElChroM industries, Dariel, IL, USA), a strontium selective resin. The resin was used once to elute the sample and then discarded. Strontium procedural blanks were less than 100 pg and hence provided a negligible contribution.

The radiogenic strontium isotope samples were analyzed on a Neptune multi-collector inductively coupled plasma mass spectrometer (MC-ICP-MS) at the Advanced Research Facilities (SGIker) of the University of the Basque Country (UPV/EHU). ⁸⁷Sr/⁸⁶Sr measurements were corrected for krypton (Kr) and rubidium (Rb) interferences and normalized for instrumental mass bias using ⁸⁷Sr/⁸⁶Sr = 8.375209. Repeated analyses of NIST SRM-987 international standard yielded a value of ⁸⁷Sr/⁸⁶Sr = 0.710262 ± 0.000026 (2σ, n = 3). Long-term ⁸⁷Sr/⁸⁶Sr value, determined over a twenty-two month period, was 0.710266 ± 0.000021 (2σ, n = 47).

Tooth enamel was also prepared for oxygen isotope analysis following the procedure described in Stephan [81]. 60 mg of dental enamel powder was processed. The organic matter was removed with a solution of 2.5% NaOCl for 24 h at room temperature followed by a 48 h treatment in 0.125M NaOH at room temperature. The hydroxyapatite powder free of organic matter was dissolved in 2 mL of HF for 24 h. The phosphate solution and the residue composed of CaF₂ were separated by centrifugation, pipetted into a 100 mL glass tube and neutralized with 3 mL 2M KOH. Silver phosphate (Ag₃PO₄) was precipitated by adding 30 mL of a buffered silver amide solution (0.2 M AgNO₃; 1.16 M NH₄NO₃; 0.75 M NH₄OH) gradually warmed up to 70°C, holding the temperature for 5–6 h and cooling down slowly. Silver phosphate crystals were filtered on a weighed 0.2μm filter and washed several times with double distilled water, then dried at 50°C for 1–2 h. Subsequently, 0.3 mg of Ag₃PO₄ was mixed with 0.5–1 mg of AgCl and 0.3 mg of graphite in silver capsules. The capsules were transferred into the autosampler carousel of the Temperature Conversion Elemental Analyser (TCEA) and degassed for 30 minutes at 80°C in a vacuum. The oxygen isotope analyses were performed on a Thermo Finnigan TCEA coupled to a Delta Plus XP Spectrometer at the University of Parma. Isotopic compositions were given in the conventional δ-notation relative to V-SMOW (Vienna-Standard Mean Ocean Water). Normalization to the V-SMOW scale was based on four replicated international reference materials provided by the International Atomic Energy Agency (IAEA): IAEA-601, IAEA-602, IAEA-CH6, and IAEA-SO-6. The analytical precision of a single determination was better than ±0.4‰.

To identify outliers in δ¹⁸O_PO4 and in ⁸⁷Sr/⁸⁶Sr within the Tauste population two statistical techniques were used. Boundaries of intra-sample variation based on two measurements of scales were defined: ± 2 standard deviation (2SD) from the mean and Tuke’s inter-quartile range method (IQR) considering 1.5xIQR and 3xIQR [82]. Parametric statistics were used to describe isotope distribution and compare isotope values between groups. Differences between sample groups were analyzed by applying an unpaired Student’s t-test. Statistical significance was accepted as p < 0.05. Statistical analysis was performed with SPSS v.20 (Statistical package for Social Sciences).

Results and discussion

Residential mobility

Strontium and oxygen isotope data for 23 tooth enamel samples, 8 rib bones and local geological materials to establish the strontium baseline signature at Tauste are shown in Table 2.

[Figure omitted. See PDF.]

Table 2. Strontium, carbon, nitrogen and oxygen isotope results for human bones and teeth, and freshwater and soil samples from Tauste.

https://doi.org/10.1371/journal.pone.0176572.t002

To establish local bioavailable strontium, bedrock, fauna, soils and surface water are used but archaeological microfauna or snail shells are considered the most appropriate material [32,43,45,38,83,84]. Since the site is a Muslim necropolis, no fauna remains are associated with the burials. Therefore, to define the Tauste bioavailable strontium isotope baseline, surface water and soils were considered. Local ⁸⁷Sr/⁸⁶Sr isotope composition determined by the local soils varies between 0.70867 and 0.70869, while the freshwater composition is 0.70843. The ⁸⁷Sr/⁸⁶Sr ratios of enamel vary between 0.70837 and 0.70867 and human bone samples range from 0.70867 to 0.70883 (Fig 6).

[Figure omitted. See PDF.]

Fig 6. Strontium isotope variation of studied samples.

Grey area corresponds to local strontium background defined by local freshwater and soils. Abbreviations (e): enamel sample, (b): bone sample.

https://doi.org/10.1371/journal.pone.0176572.g006

According to the defined local baseline, most individuals buried in Tauste have a ⁸⁷Sr/⁸⁶Sr signature consistent with local origin. Only two individuals plot outside the estimated local compositional range (1σ) according to outliers identifying by 1.5xIQR and 2SD method. Enamel of young adult female T-24 displays a lower radiogenic strontium isotope value (⁸⁷Sr/⁸⁶Sr = 0.70837) and the rib of adult female T-32(b) presents a higher strontium value (⁸⁷Sr/⁸⁶Sr = 0.70883) (Fig 6). These compositions indicate two different mobility patterns for these females. The female T-24’s strontium value suggests she was born in another place and died in Tauste. Whereas female T-32(b)’s isotope value suggests that she spent her childhood in Tauste, during adulthood moved to another location and came back to Tauste a few years before she died.

The enamel phosphate oxygen ratios (δ¹⁸O_p) cover a broad range of values from 14.25‰ to 19.30‰. Based on kernel density estimations (Fig 7), data can be split into three groups: a larger group (n = 18) with δ¹⁸O_p ratios between 16.4‰ and 18‰ and two smaller groups, one formed by two males with δ¹⁸O_p signature higher than 19.2‰ and the other formed by two females with δ¹⁸O_p signature lower than 14.4‰. Local meteoric water δ¹⁸O_dw is -5.6‰ (average values from 2000 to 2006) considering data from Zaragoza airport station [85]. δ¹⁸O_dw values for humans from Tauste were calculated using the available phosphate/drinking water equations [58,60,61,63,66], and comparing them with expected local water values derived from IAEA/WISER data set [85]. The equation by Iacumin and Venturelli [66] was used to estimate the drinking water isotope value. The larger group of individuals show calculated δ¹⁸O_dw ranging from -3.7‰ to -6.9‰, consistent with local meteoric water which ranges between -4.1‰ and -6.3‰ as annual average [85]. The number of outliers identified was determined using 1.5xIQR and 2SD statistical methods. The two males and two females whose isotope values fall outside the larger group may come from a warmer, more coastal or possible more arid climate, and from a colder or higher altitude region, respectively. When oxygen isotope data are compared with ⁸⁷Sr/⁸⁶Sr ratios of the same teeth the individuals in the three groups display isotopic values compatible with bio-accessible strontium measured to establish the local signature. Only the female T-24 falls strictly outside the expected strontium range for local origin although she falls into the expected range of the calculated drinking water values for Tauste (Fig 7). Similar values of δ¹⁸O for Tauste meteoric water showed a broad geographic distribution (Fig 8) overlapping different geological areas. The lack of a strontium isotope composition database in the Iberian Peninsula prevents a determination of the regional provenance of the non-local female T-24.

[Figure omitted. See PDF.]

Fig 7. δ¹⁸O versus ⁸⁷Sr/⁸⁶Sr in the tooth enamel of Tauste individuals.

Kernel density contour lines represent 10%.

https://doi.org/10.1371/journal.pone.0176572.g007

[Figure omitted. See PDF.]

Fig 8.

A. Large-scale map of oxygen isotope signatures in the Iberian Peninsula [86,87]. B. Detailed map of oxygen isotope signature of the Ebro Valley [86,87].

https://doi.org/10.1371/journal.pone.0176572.g008

During Muslim period the Ebro valley was a trade route between Mediterranean coast towards north of the Iberian Peninsula and trans-Pyrenean Europe. Taking on account both oxygen and strontium isotope data the origin areas for these two outliers groups would be confined to the Ebro Valley. Though, during the Early Middle Ages few people regularly travelled because it was simply too difficult and too dangerous. Muslim females were subject to inbreeding marriages so that Tauste outlier females would move by patrilocal marriages. Males would move to get better economic opportunities and possibility of improving status. Tauste non-local males would move from farms to a large village to find new kinds of works. The political instability of the Upper March frontier region also favored displacement of people towards larger and safer urban centers [88].

Dietary reconstruction

The bone collagen obtained was very well preserved with an average yield of 6.27±4.19%wt (1s.d.). The content of carbon and nitrogen in bone collagen was about 13.1–16.8% and 25.6%–36.2%, respectively, so well-preserved bone collagen should display a carbon/nitrogen molar ratio, based on the content (in %) of these elements in the sample between 2.9 and 3.6 [89,90]. The individual data are given in Table 2.

The δ¹³C ratio for human bone samples ranges between -17.0‰ and -19.9‰ with a mean value of -19.1±0.50‰, and δ¹⁵N ratios range between 9.6‰ and 17.5‰ with a mean value of 14.9±1.74‰. The rather strong δ¹⁵N signal of nearly all individuals is noteworthy, with values more than 5‰ over the terrestrial ecosystem baseline (Fig 7), an offset unusually large for a single trophic level effect [91,92,21].

Individuals from Tauste were compared with broadly coetaneous Muslim and Christian populations at several locations in the Iberian Peninsula (Fig 9). There are no significant differences in δ¹³C ratios between Tauste (-19.1±0.5‰) and most neighboring Muslim populations (δ¹³C -19.0±0.3‰ in Zaragoza and -19±0.2‰ in Albarracín), or between Tauste and Christian populations (δ¹³C -18.4±1.1‰ in Jaca, -18.4 ±0.6‰ in Valencia, -19.0±1‰ in Aistra, -19.8±0.7‰ in Zaballa, -18.1±1.1‰ in Zornoztegi and -19.6±0.7‰ in Treviño [93–95]) (Table 3, Fig 10). However, the Benipeixcar Muslim population showed enrichment in δ¹³C (-16.36±0.9‰) attributed to marine resource consumption [95]. In contrast, the δ¹⁵N ratios at Tauste (average 15.0±1.7‰) are unusually high compared with contemporaneous Christian and Muslim populations in the Iberian Middle Ages, whose average values are lower 11‰ (Table 3).

[Figure omitted. See PDF.]

Fig 9. Location of Muslim and Christian archaeological sites and the Upper March or Muslim northern frontier during the 9th century.

https://doi.org/10.1371/journal.pone.0176572.g009

[Figure omitted. See PDF.]

Fig 10. Plot of δ¹³C and δ¹⁵N values for the individuals at Tauste.

Ecosystem boxes are based on faunal data reported by [90] and [92]. Field values are corrected for trophic level effect. 1–9 corresponds to coetaneous Muslim and Christian populations in the Iberian Peninsula: (1) Zornostegui 12th-14th centuries, (2) Aistra 8th-9th centuries, (3) Zaballa 10th-15th centuries, (4) Treviño 12th-14th centuries, (5) Jaca 13th-15th centuries (6) Valencia 14th-15th centuries (7) Albarracín 10th-12th centuries (8) Zaragoza 10th-12th centuries; (9) Benipeixcar 15th-16th centuries (1)-(4) [93], (5)-(8) [94], (9) [95].

https://doi.org/10.1371/journal.pone.0176572.g010

[Figure omitted. See PDF.]

Table 3. Mean collagen δ¹³C and δ¹⁵N of individuals at Tauste and coeval Christian and Islamic sites.

https://doi.org/10.1371/journal.pone.0176572.t003

Faunal samples are required to strengthen conclusions about human diet. Establishing the local isotope composition baseline was problematic since the Islamic burial ritual forbade any objects being buried with the body. Additionally, there is no evidence of Muslim settlement in Tauste to provide coeval fauna. Thus it was not possible to obtain the local baseline with archaeozoological data. Furthermore, present-day local faunal must be discarded because livestock are fed with non-local resources and the isotopic signal will not correspond to local plant resources. Pasture-fed livestock will also exhibit nitrogen isotope depletion due to the widespread use of mineral fertilizers [96,97]. In addition to the fodder and fertilizer effect, livestock trade is another factor affecting isotope composition. All these factors prevented the use of present-day macro-mammals to establish the carbon and nitrogen isotope local baseline.

For these reasons, two local species of small mammals and freshwater fish were analyzed to establish the dietary baseline for Tauste medieval population. Small mammals were selected since present low mobility with restricted home ranges more accurately reflects the local isotopic composition baseline. The analyzed species were wood mice (Apodemus sylvaticus), as they are herbivorous (seed eaters), and shrews (Crocidura russula), which are one level higher than wood mice in the trophic chain.

To address this issue, the isotope enrichment of one trophic level has to be established. For this purpose, both faunal and human values from the nearest coeval Christian and Muslim settlements were compared (Table 4, Fig 11). Livestock from Alava archaeological sites and Benipeixcar Muslim site displayed different average δ¹³C and δ¹⁵N values. The nitrogen values of plants and fauna are strongly influenced by local climate conditions [98]. Fauna at the Alava sites, located in the northern Iberian Peninsula, displayed lower nitrogen isotope values due to a more humid climate than Benipeixcar, on the warmer and drier Mediterranean coastline.

[Figure omitted. See PDF.]

Fig 11. Plot of δ¹³C and δ¹⁵N values of present-day fauna and medieval human at Tauste and the comparative Alava and Benipeixcar archaeofauna and human dataset [91,93].

Colored rectangles indicate local carbon and nitrogen baseline. Doted arrows show the offset between local baseline and human δ¹⁵N due to one trophic level enrichment. Grey rectangle indicates estimated local carbon and nitrogen baseline at Tauste.

https://doi.org/10.1371/journal.pone.0176572.g011

[Figure omitted. See PDF.]

Table 4. Present-day fauna carbon and nitrogen isotope data from Tauste region and archaeofauna data from their coeval sites of Alava and Beneipeixcar.

https://doi.org/10.1371/journal.pone.0176572.t004

Comparison of human values with fauna revealed the average offset between human and livestock of c. 1.5‰ in δ¹³C and 4–4.5‰ in δ¹⁵N (Fig 11) typical of one trophic level [99].

In an attempt to establish the trophic level offset in Tauste, small mammal isotope values were considered. The average of wood mice values (-18.4±0.9‰ in δ¹³C and 8.7±0.7‰ in δ¹⁵N) and shrew values (-19.1±2.8‰ in δ¹³C and 12.4±2.2‰ in δ¹⁵N) exhibited an offset of c. -1‰ in δ¹³C and 4‰ in δ¹⁵N. Considering the predictable trophic level offset, the isotope values of the Tauste population should fall within the average of c. -21±1‰ in δ¹³C and 11±0.5‰ in δ¹⁵N, which is the enrichment typical of one trophic level. Such an enriched δ¹⁵N baseline can be explained by regional environmental conditions. Regional aridity and the local high relative salinity led to this δ¹⁵N enrichment [100]. In fact, chemical analysis performed in teeth from Tauste individuals showed the log (Ba/Sr) = -2.35; similar high concentrations of strontium are found in both arid/semiarid region soils and in marine food sources [101,102]. Tauste site is located far from the coastline and the possibility large-scale consumption of marine food can be discarded. Thus Ba/Sr values and the δ¹⁵N enrichment can be explained by the local bedrock (gypsum and salt) and environmental factors in the Tauste region.

The offset between the isotope ratios of the mice collagen (as equivalent to livestock) and the average of the expected range for human diet suggested a mixture of protein from terrestrial and freshwater sources in their foodstuffs. The relative proportion of terrestrial and freshwater resources in diet can be estimated c. 50–50% using a simple linear mixing model.

However, a particularly high consumption of freshwater fish does not appear justified from either historical or anthropological points of view. Islamic texts about the daily diet in Muslim medieval Spain indicate it was based mainly on cereals: wheat (Tritium), barley (Hordeum vulgare) and rye (Secale cereale), together with such C₄ grains as millet (Pennisetum) and sorghum (Sorghum), vegetables and pulses, such as chickpeas (Cicer arietinum), lentils (Lens culinaris), and peas (Pisum sativum), with some regional differences [103–105]. The main sources of proteins were meat and pulses, and the type and quantity of protein consumed varied according to social status and gender. The most highly regarded meats were lamb and poultry. Pork and any animal not slaughtered in a way considered halal were excluded from the diet because Islamic law forbade it. People from lowly backgrounds consumed little meat and often made do with offal because it was cheap. Pulses such as broad beans, chickpeas and lentils constituted another source of protein and were classified as medical food. Although fish was not considered food of great dietetic value, it was part of the diet of the people from more humble backgrounds, particularly in river or coastal areas [106].

Another possible interpretation for the high δ¹⁵N enrichment in the Tauste population is to consider the manuring effect on plants. Since the use of manure was an advanced agricultural method introduced by Muslims in the Iberian Peninsula [107] all plants from manured soils showed a δ¹⁵N enrichment that can be about 5‰ in cereal δ¹⁵N [30]. Consequently, domestic animals foddered with manured chaff and grains will present higher δ¹⁵N values. Assuming the consumption of plants enriched in δ¹⁵N by the manure effect, the contribution of freshwater resources may decrease and become less than 10% of dietary protein intake.

Another interesting aspect of the reconstruction of the Tauste inhabitants’ palaeodiet is the variation in isotope composition by sex and age. For the comparison, female T-32 was excluded because her strontium isotope composition indicates a return to Tauste in the last years of her life preventing δ¹⁵N and δ¹³C remodeling to the local isotope signature. The mean values of the male samples are 15.5±1.00‰ for δ¹⁵N and -19.2±0.31‰ for δ¹³C. In female samples, the mean values are 14.6±1.15‰ for δ¹⁵N and -19.0±0.32‰ for δ¹³C. Comparison of the results with the Student's t-test for males and females confirms the significant differences in δ¹⁵N (t₁₈, p = 0.026). These differences suggest that females and males had different access to animal protein, probably due to the sexual division of labor [108]. In fact, written evidence suggests the diet of Muslim females and males differed. Medical treatises make recommendations about diet, and special recommendations were made for females during pregnancy or lactation [109]. Two meals a day were recommended but in practice a smaller meal as breakfast was taken by working males. Lower dietetic needs were expected of females since female labor was restricted to domestic tasks and other activities in the household [110]. However, the differences in nitrogen composition related to sex are determined by the middle-aged adult male (ages 35–50 years) composition. In fact, excluding middle-aged adult males differences between sexes are not significant. Middle-aged adult males show higher δ¹⁵N values than young males and females of any age. Variations between adult males and the other age segments of population indicate a different dietary intake which may be related, for instance, to a greater meat protein consumption.

Finally, the highest δ¹⁵N values were found in children younger than 3 years, whereas children older than 4 years old have similar δ¹⁵N values to those of adults. In younger infants, higher δ¹⁵N values are due to a “nursing effect” indicating a diet based mainly on maternal milk.

However considering the significantly wide chronological span and the small number of radiocarbon dates, it is not possible to rule out temporal dietary differences, which is to say that there may have been variations in the diet in the course of time.

Conclusions

Isotope signatures in archaeological human remains have been used to investigate palaeomobility and the palaeodiet of the medieval Muslim population at Tauste, on the changing frontier between Muslim and Christian kingdoms in the 8th to 10th centuries.

The combination of strontium and oxygen isotope analyses was able to discriminate non-local and local individuals. Although Tauste was located on the northern al-Andalus frontier most individuals were of local origin and only three females and two males were non-locals. Establishing the provenance of incoming individuals is difficult as strontium isotope ratios indicate a similar geological region. According to the oxygen isotope composition, two males would come from a warmer region while two females would come from a more mountainous geographical area. Also T-24 was of non-local origin since the different strontium values indicate a different geological provenance. Within the local population, the female T-32 stand out because she was born in Tauste, some years later she moved and lived in another geological region and then returned to Tauste a few years before she died.

As regards the medieval Muslim diet, the δ¹³C and δ¹⁵N results illustrate not only differences in diet according to sex and age but also the environmental conditions. The extremely high δ¹⁵N values in Tauste population (δ¹⁵N = 15.0‰ on average) indicate an anomalously high δ¹⁵N baseline that can be explained by the concurrence of (1) geological and environmental conditions, (2) the manuring effect on vegetables, cereals and livestock and (3) the consumption of freshwater fish. The amount of fish in the diet varies from 50% to <20% as the manuring effect increases. Significant differences were observed in δ¹³C and δ¹⁵N by sex, indicating different diets that may be related to the sexual division of labor since Muslim female work was restricted to the household. The main dietary differences between males and females were amongst adult individuals, suggesting adult males were differentially valued in medieval Muslim society and consumed more animal protein than females and young males. The lower δ¹⁵N values of the elder females indicated lower protein consumption due to lesser dietary needs. In contrast, the significant higher δ¹⁵N values in the children younger than 4 years could be related to the “nursing effect”.

Supporting information

[Figure omitted. See PDF.]

S1 Table. Radiocarbon data with the 1 sigma and 2 sigma calibrations.

Dating of human bone samples from Tauste Muslim necropolis was carried out in Roca Solano Laboratory (CSIC Madrid).

https://doi.org/10.1371/journal.pone.0176572.s001

(XLSX)

S2 Table. Strontium, oxygen, carbon and nitrogen isotopes analyses on human samples.

https://doi.org/10.1371/journal.pone.0176572.s002

(XLSX)

S3 Table. Strontium isotope analyses on water and soils samples used as comparative data for the local range at Tauste.

https://doi.org/10.1371/journal.pone.0176572.s003

(XLSX)

S4 Table. Results of the carbon and nitrogen isotope analysis of the present-days faunal remains from Tauste used as comparative data for the local isotope baseline.

https://doi.org/10.1371/journal.pone.0176572.s004

(XLSX)

Acknowledgments

The authors would like to thank the two anonymous referees for their comments and suggestions on the manuscript. They also would like to thank Peter Smith for language assistance.

Author Contributions

1. Conceptualization: IG LO MCZ.

2. Funding acquisition: LO MCZ AAO XM.

3. Investigation: IG LO MCZ AAO.

4. Resources: XM MP FJG.

5. Writing – original draft: IG LO MCZ.

6. Writing – review & editing: IG LO MCZ MP AAO PI.

Citation: Guede I, Ortega LA, Zuluaga MC, Alonso-Olazabal A, Murelaga X, Pina M, et al. (2017) Isotope analyses to explore diet and mobility in a medieval Muslim population at Tauste (NE Spain). PLoS ONE 12(5): e0176572. https://doi.org/10.1371/journal.pone.0176572

References

1. García de Cortázar JA, Sesma JA. Manual de Historia Medieval. Madrid: Alianza Universidad; 2008.

2. Hitchcock R. Muslim Spain reconsidered from 711 to 1502. Edinburgh: Edinburgh University Press; 2014.

3. Cañada Juste A. Los Banu Qasi (714–924). Principe de Viana. 1980; 158–159: 5–96.

4. Coope JA. An Etiquette for Women: Women’s Experience of Islam in Muslim Spain. Essays in Medieval Studies. 2014; 29: 76.

5. Petersen A. The Archaeology of Death and Burial in the Islamic World. In: Nilsson Stutz L, Tarlow S, editors. The Oxford Handbook of the Archaeology of Death and Burial; 2013. pp. 241–258.

6. Ambrose SH. Isotopic analysis of paleodiets: methodological and interpretive considerations. In: Sandford MK, editor. Investigations of Ancient Human Tissue: Chemical Analysis in Anthropology. Langthorne: Gordon and Breach; 1993. pp. 59–30.

7. DeNiro MJ, Epstein S. Influence of diet on the distribution of carbon isotopes in animals. Geochim Cosmochin Ac. 1978; 42: 495–506.

8. Hedges REM, Reynard LM. Nitrogen isotopes and the trophic level of humans in archaeology. J Archaeol Sci. 2007; 34 (8):1240–1251.

9. Katzenberg A. Stable Isotope Analysis: A Tool for Studying Past Diet, Demography, and Life History. In: Katzenberg MA, Saunders SR, editors. Biological Anthropology of the Human Skeleton. 2nd ed. Hoboken, NJ: John Wiley and Sons; 2008. pp. 413–441.

10. Lee-Thorp JA. On Isotopes and Old Bones. Archaeometry. 2008; 50 (6): 925–950.

11. Schwarcz HP, Schoeninger MJ. Stable isotope analyses in human nutritional ecology. Yearb Phys Anthropol. 1991; 34: 283–321.

12. Meier-Augenstein W. Stable isotope forensics: an introduction to the forensic application of stable isotope analysis. Oxford: Wiley-Blackwell; 2010.

13. Smith BN, Epstein S. Two categories of 13C/12C ratios for higher plants. Plant Physiol. 1971; 47: 380–384. pmid:16657626

14. O’Leary MH. Carbon isotope fractionation in plants. Phytochemistry. 1981; 20: 552–567.

15. Hoefs J. Stable Isotope Geochemistry. 1st ed. Berlin: Springer-Verlag; 2009.

16. White WM. Isotope Geochemistry. Chichester: Wiley-Blackwell; 2015.

17. Malainey ME. A consumer's guide to archaeological science: analytical techniques. New York:Springer; 2011.

18. Dufour E, Bocherens A, Mariotti A. Palaeodietary implications of isotopic variability in Eurasian lacustrine fish. J Archaeol Sci. 1999; 26: 617–627.

19. Pazdur A, Goslarl T, Pawlyt M, Hercman H, Gradzinski M. Variations of isotopic composition of carbon in the karst environment from southern Poland, present and past. Radiocarbon. 1999; 41 (1): 81–97.

20. Sandford MK. Investigations of Ancient Human Tissue: Chemical Analyses. Langhorne: Gordon and Breach; 1993.

21. Bocherens H, Drucker D. Trophic level isotopic enrichment of carbon and nitrogen in bone collagen: Case studies from recent and ancient terrestrial ecosystems. Int J Osteoarchaeol. 2003; 13: 46–53.

22. Schoeninger M.J, Moore K, Bone stable isotope studies in archaeology. J. World Prehist. 1992; 6: 247–296.

23. Schoeninger MJ, Deniro MJ, Tauber H. Stable nitrogen isotope ratios of bone collagen reflect marine and terrestrial components of prehistoric human diet. Science. 1983; 220: 1381–1383.24. pmid:6344217

24. Schoeninger MJ, DeNiro MJ. Nitrogen and carbon isotopic composition of bone collagen from marine and terrestrial animals. Geochim Cosmochim Acta. 1984; 48: 625–639.

25. Walker PL, DeNiro MJ. Stable nitrogen and carbon isotope ratios in bone collagen as indices of prehistoric dietary dependence on marine and terrestrial resources in southern California. Am J Phys Anthropol. 1986; 71: 51–61. pmid:3777147

26. Katzenberg AM, Weber A. Stable isotope ecology and paleodiet in the Lake Baikal region of Siberia. J Archaeol Sci. 1999; 26: 651–659.

27. Fuller BT, De Cupere B, Marinova E, Van Neer W, Waelkens M, Richards MP. Isotopic reconstruction of human diet and animal husbandry practices during the Classical-Hellenistic, Imperial and Byzantine periods at Sagalassos, Turkey. Am J Phys Anthropol. 2012; 149: 157–171. pmid:22729657

28. Robson HK, Andersen SH, Clarke L, Craig OE, Gron KJ, Jones AKG, et al. Carbon and nitrogen stable isotope values in freshwater, brackish and marine fish bone collagen from Mesolithic and Neolithic sites in central and northern Europe. Environmental Archaeology; 2015.

29. Bogaard A, Heaton THE, Poulton P, Merbach I. The impact of manuring on nitrogen isotope ratios in cereals: archaeological implications for reconstruction of diet and crop management practices. J Archaeol Sci. 2007; 34: 335–343.

30. Fraser RA, Bogaard A, Heaton T, Charles M, Jones G, Christensen BT, et al. Manuring and stable nitrogen isotope ratios in cereals and pulses: towards a new archaeobotanical approach to the inference of land use and dietary practices. J Archaeo Sci. 2011; 38: 2790–2804.

31. Schwarcz HP, Schoeninger MJ. Stable isotope analyses in human nutritional ecology. Yearb Phys Anthropol. 1991; 34: 283–321.

32. Phillips DL, Gregg JW. Source partitioning using stable isotopes: coping with too many sources. Oecol. 2003; 136 (2): 261–269.

33. Fogel ML, Tuross N, Owsley DW. Nitrogen Isotope Tracers of Human Lactation in Modern and Archaeological Populations. In: Yearbook of the Carnegie Institution Geophysical Laboratory. Washington DC: Carnegie Institution of Washington; 1989. pp. 111–117.

34. Bentley RA, Knipper C. Geographical patterns in biologically available strontium, carbon and oxygen isotope signatures in prehistoric SW Germany. Archaeometry. 2005; 47: 629–644.

35. Evans J, Chenery C, Fitzpatrick AP. Bronze age childhood migration of individuals near Stonehenge, revealed by strontium and oxygen isotope tooth enamel analysis. Archaeometry. 2006a; 48: 309–321.

36. Evans J, Stoodly N, Chenery N. A strontium and oxygen isotope assessment of a possible fourth century immigrant population in a Hampshire cemetery, southern England. J Archaeol Sci. 2006b; 33: 265–272.

37. Hillson S. Dental Anthropology. Cambridge: Cambridge University Press; 1996.

38. Price TD, Burton JH, Bentley RA. The characterization of biological available strontium isotope ratios for the study of prehistoric migration. Archaeometry 2002; 44: 117–136.

39. Hoppe KA, Koch PL, Furutani TT. Assessing the preservation of biogenic strontium in fossil bones and tooth enamel. Int J Osteoarchaeol. 2003; 13 (1–2): 20–28.

40. Hillson S. Teeth. Cambridge: Cambridge University Press; 1986.

41. Burton JH. Trace elements in bone as paleodietary indicators. In: Orna MV, editors. Archaeological Chemistry: Organic, Inorganic, and Biochemical Analysis. Washington, DC: American Chemical Society; 1996. pp. 327–333.

42. Hoppe KA, Koch PL, Carlson RW, Webb SD. Tracking mammoths and mastodons: Reconstruction of migratory behaviour using strontium isotope ratios. Geology. 1999; 27 (5): 439–442.

43. Blum JD, Taliaferro EH, Weisse MT, Holmes RT. Changes in Sr/Ca, Ba/Ca and Sr-87/Sr-86 ratios between trophic levels in two forest ecosystems in the northeastern USA. Biogeochemistry. 2000; 49 (1): 87–101.

44. Balter V, Person A, Labourdette N, Drucker D, Renard M, Vandermeersch B, Were Neandertalians essentially carnivores? Sr and Ba preliminary results of the mammalian palaeobiocoenosis of Saint-Cesaire. Comptes Rendus De L Academie Des Sciences Serie Ii Fascicule a-Sciences De La Terre Et Des Planetes. 2001; 332 (1): 59–65.

45. Balter V, Bocherens H, Person A, Labourdette N, Renard M, Vandermeersch B, Ecological and physiological variability of Sr/Ca and Ba/Ca in mammals of West European mid-Wurmian food webs. Palaeogeography, Palaeoclimatology, Palaeoecology. 2002; 186 (1–2): 127–143.

46. Balter V. Allometric constraints on Sr/Ca and Ba/Ca partitioning in terrestrial mammalian trophic chains. Oecol. 2004; 139 (1): 83–88.

47. Faure G, Mensing TM. Isotopes: principles and applications. Hoboken, NJ: Wiley; 2005.

48. Hoppe KA, Koch PL. Reconstructing the migration patterns of late Pleistocene mammals from northern Florida, USA. Quaternary Res. 2007; 68 (3): 347–352.

49. Bentley RA. Strontium isotopes from the earth to the archaeological skeleton: a review. J Archaeol Method Th. 2006; 13 (3): 135–187.

50. Bentley RA, Price TD, Stephan E. Determining the ‘local’ 87Sr/86Sr range for archaeological skeletons: A case study from Neolithic Europe. J Arch Sci. 2004; 31: 365–375.

51. Evans JA, Montgomery J, Wildman G, Boulton N. Spatial variations in biosphere 87Sr/86Sr in Britain. J Geol Soc London. 2010; 167 (1): 1–4.

52. Tütken T, Vennemann TW, Pfretzschner HU. Nd and Sr isotope compositions in modern and fossil bones-proxies for vertebrate provenance and taphonomy. Geochim Cosmochim Acta. 2011; 75: 5951–5970.

53. Böhlke JK, Horan M. Strontium isotope geochemistry of ground waters and streams affected by agriculture, Locust Grove, Maryland. Appl Geochem. 2000; 15: 599–609.

54. West JB, Hurley JM, Dudas FO, Ehleringer JR. The stable isotope ratios of Marijuana. II. Strontium isotopes relate to geographic origin. J Forensic Sci. 2009; 54: 1261–1269. pmid:19804525

55. Tichomirowa M, Heidel C, Junghans M, Haubrich F, Matschullat J. Sulfate and strontium water source identification by O, S and Sr isotopes and their temporal changes (1997–2008) in the region of Freiberg, central-eastern Germany. Chem Geol. 2010; 276: 104–118.

56. Christian LN, Banner JL, Mack LE. Sr isotopes as tracers of anthropogenic influences on stream water in the Austin, Texas, area. Chem Geol. 2011; 282: 84–97.

57. Longinelli A. Oxygen isotopes in mammal bone phosphate: A new tool for paleohydrological and paleoclimatological search? Geochim Cosmochim Acta. 1984; 48: 385–390.

58. White CD, Spence MW, Stuart-Williams HLQ, Schwarcz HP. Oxygen isotopes and the identification of geographical origins: the Valley of Oaxaca versus the Valley of Mexico. J Archaeol Sci. 1998; 25: 643–655.

59. Darling WG, Bath AH, Gibson JJ, Rozanski K. Isotopes in water. In: Leng MJ, editor. Isotopes in Palaeoenvironmental Research. London: Springer; 2006. pp. 1–52.

60. Daux V, Lécuyer C, Heran MA, Amiot R, Simon L, Fourel F, et al. Oxygen isotope fractionation between human phosphate and water revisited. J Hum Evol. 2008; 55: 1138–1147. pmid:18721999

61. Luz B, Kolodny Y, Horowitz M. Fractionation of oxygen isotopes between mammalian bone-phosphate and environmental drinking water. Geochim Cosmochim Acta. 1984; 48: 1689–1693.

62. Luz B, Kolodny Y. Oxygen isotope variations in phosphate of biogenic apatite, IV. Mammal teeth and bones. Earth Planet Sci Letters. 1985; 75:29–36.

63. Levinson AA, Luz B, Kolodny Y. Variations in oxygen isotopic compositions of human teeth and urinary stones. Appl Geochem. 1987; 2: 367–371.

64. Bryant JD, Froelich PN. A model of oxygen isotope fractionation in body water of large mammals. Geochim Cosmochim Acta. 1995; 59: 4523–4537.

65. Kohn MJ. Predicting animal δ18O: Accounting for diet and physiological adaption. Geochim Cosmochim Acta. 1996; 60 (23): 4811–4829.

66. Iacumin P, Venturelli G. The δ18O of phosphate of ancient human biogenic apatite can really be used for quantitative palaeoclimate reconstruction? European Scientific Journal. 2015; 11 (9): 221–235.

67. Stutz LN, Tarlow S. The Oxford Handbook of the Archaeology of Death and Burial. Oxford university Press; 2013.68.

68. Serraro Peña JL, Jimenez Morillas Y, Alcala Lirio F, Cano Carrillo J. Intervención arqueológica de urgencia en la Calle A del SUNP1 de Marroquíes Bajos (Jaén). Archivo de la Delegación Provincial de Cultura de Jaén, 2000.

69. Olcina Doménech MH, Tendero Porras E, Guilabert Mas AP. Lucentum: Anales de la universidad de Alicante. Prehistoria, arqueología e historia antigua. 2008; 27:213–228.

70. Gutiérrez FJ, Pina M. El cementerio andalusí de Tauste. Tauste en su historia. Actas de las XII Jornadas sobre la Historia de Tauste. Asociación Cultural “El Patiaz”. 2011; 67–113.

71. Gutiérrez FJ, Pina M, Laliena C. La maqbara de Tauste. Primeras investigaciones. Acta I Congreso Aragonés de Arqueología y Patrimonio. 2016; 415–424.

72. Bronk Ramsey C, Lee S. Recent and Planned developments of the program Oxcal. Radiocarbon. 2013, 55 (2–3): 720–730.

73. Reimer PJ, Bard E, Bayliss A, Beck J W, Blackwell PG, Ramsey CB, et al. Intcal13 and marine13 radiocarbon age calibration curves 0–50,000 years cal bp. Radiocarbon. 2013; 55 (4): 1869–1887.

74. James D. Early Islamic Spain: The History of Ibn Al-Qutiyah, New York: Routledge; 2009.

75. Safran JM. Defining Boundaries in al-Andalus: Muslims, Christians, and Jews in Islamic. New York: Cornel University Spress; 2013.

76. Rodríguez Fernández LR, López Olmedo F, Oliveira JT, Medialdea T, Terrinha P, Matas J, et al. Mapa Geológico de la Península Ibérica, Baleares y Canarias a escala 1:1.000.000. IGME, Madrid; 2015.

77. Navas A. la participación de los yesos en la salinización de las aguas superficiales de la cuenca de Ebro. Metodología de cuantificación. An Aula Dei. 1988; 19 (3–4): 345–359.

78. Causapé Valenzuela J. Calidad de los ríos Riguel y Arba (Zaragoza). Influencia del regadío de Bardenas y modelización geoquímica del sistema. Rev. Real Academia de Ciencias. 2003; 58: 7–36.

79. White T, Black M, Folkens P. Human osteology. San Diego: Academic Press; 1991.

80. Bocherens H., Fizet M., Mariotti A., Lange-Badre B., Vandermeersch B., Borel J.P., Bellon G., Isotopic biogeochemistry (13C, 15N) of fossil vertebrate collagen: application to the study of a past food web including Neanderthal man, J. Hum. Evol. 1991; 20 (6): 481–492.

81. Stephan E. Oxygen isotope analysis of animal bone phosphate: method refinement, influence of consolidants, and reconstruction of palaeotemperatures for Holocene sites. J Archaeol Sci. 2000; 25: 523–535.

82. Lightfood E, O’Connell TC. On the use of biomineral oxygen isotope data to identify human migrants in the archaeological record: intra-sample variation, statistical, methods and geographical considerations. Plos One. 2016; 11(4): e0153850. pmid:27124001

83. Maurer AF, Galer SJG, Knipper C, Beierlein L, Nunn EV, Peters D et al. Bioavailable 87Sr/86Sr in different environmental samples-Effects of anthropogenic contamination and implications for isoscapes in past migration studies. Sci Total Environ. 2012; 433: 216–29. pmid:22796412

84. Voerkelius S, Lorenz GD, Rummel S, Quétel CR, Heiss G, Baxter M. Strontium isotopic signatures of natural mineral waters, the reference to a simple geological map and its potential for authentication of food. Food Chem. 2010; 118: 933–40.

85. Global Network of Isotopes in Precipitation [Internet]. The GNIP Database. Available from: http://www.iaea.org/water. 2015.

86. Bittel LM. Women in Early Medieval Europe, 400–1100. Cambrige: Cambrige University Press; 2002.

87. Gridded maps of the isotopic composition of meteoric waters [Internet]. Available from: http://www.waterisotopes.org. 2015.

88. Capilla JE, Rodriguez Arevalo J, Castaño Castaño S, Díaz Teijeiro MF, Sanchez del Moral R, Heredia Diaz J. Mapping oxygen-18 in meteoric precipitation over Peninsular Spain using geostatistical tools. In: Fall Meeting, AGU (abstract No. H34C-04). San Francisco, CA, USA: American Geophysical Union; 2011.

89. DeNiro MJ. Post-mortem preservation and alteration of in vivo bone collagen isotope ratios in relation to paleodietary reconstruction. Nature. 1985; 31:806–809.

90. Ambrose SJ. Preparation and characterization of bone and tooth collagen from isotopic analysis. J Archaeol Sci. 1990; 17: 431–451.

91. Barrett ELB, Moore AJ, Moore PJ. Diet and social conditions during sexual maturation have unpredictable influences on female life history trade-offs. J Evolution Biol. 2002; (3): 571–581.

92. Müldner G, Richards MP. Diet and diversity at later medieval Fishergate: the isotopic evidence. Am J Phys Anthropol. 2007; 134: 162–74. pmid:17568445

93. Quirós Castillo JA. Los comportamientos alimentarios del campesinado medieval en el País Vasco y su entorno (siglos VIII-XIV). Historia agraria: Revista de agricultura e historia rural. 2013; 59: 13–41.

94. Mundee M. Exploring Diet and Society in Medieval Spain: New Approaches Using Stable Isotope Analysis. PhD Thesis, Durham University. 2010. Available from: http://etheses.dur.ac.uk/434/.

95. Alexander MM, Gerard CM, Gutiérrez A, Millard AR. Diet, society, and economy in late medieval Spain: Stable isotope evidence from Muslims and Christians from Gandía, Valencia. Am J Phys Anthropol. 2014; 156 (2): 263–273. pmid:25351146

96. Bol R, Jorgen Eriksen J, Smith P, Garnett MH, Coleman K, Christensen BT. The natural abundance of 13C, 15N, 34S and 14C in archived (1923–2000) plant and soil samples from the Askov long-term experiments on animal manure and mineral fertilizer Rapid Commun Mass Spectrom. 2005; 19: 3216–3226. pmid:16220463

97. White WM. Geochemistry. New York: Willey Backwell; 2013. https://doi.org/10.1016/j.chemer.2012.09.001

98. Ambrose SH. Effects of diet, climate and physiology on nitrogen isotope abundances in terrestrial foowebs. J Archaeol Sci. 1991; 18: 293–317.

99. Druker D, Bocherens H. Carbon and nitrogen isotopes as tracers of change in diet breadth during Middle and Upper Paleolithic in Europe. Int J Osteoarchaeol. 2004; 14: 162–177.

100. Heaton TH. The 15N/14N ratio of plants in South Africa and Namibia: relationship to climate and coastal/saline environments. Oecol. 1987; 74 (2): 236–246.

101. Guede I, Zuluaga MC, Ortega LA, Alonso-Olazabal A, Murelaga X, Piña M, et al. Analyses of human dentine and tooth enamel by laser ablation-inductively coupled plasma-mass spectrometry (LA-ICP-MS) to study the diet of medieval Muslim individuals from Tauste (Spain). Microchem J. 2017; 130: 287–294.

102. Burton JH, Price TD. The ratio of barium to strontium as a paleodietary indicator of consumption of marine resources. J Archaeol Sci. 1990; 17 (5): 547–557.

103. García Sánchez E. La gastronomía andalusí. In: Salvatierra Cuenca V, editor. El zoco: vida económica y artes tradicionales en el al-Andalus y Marruecos. Madrid: Lunwerg; 1995. pp. 49–57.

104. Salas-Salvadó J, Huetos-solano Maria D.;Garcia-Lorda Pilar, Bullo Monica. Diet and dietetics in al-Andalus. British Journal of Nutrition 2006; 96 (1): S100–S114.

105. Garcia Sanchez E. Dietetic aspect of food in al–Andalus. In: Waines D, editors. Patterns of Everyday Life. Aldeshot: Ashgate; 2002. pp. 275–288.

106. Anderson G.D. Food and diet. In: Meri JW, editor. Medieval Islamic Civilization: an Encyclopedia. New York: Routlegde; 2006. pp. 264–266.

107. Glick T. Islamic and Christian Spain in the Early Middle Ages. 2nd ed. Leiden: Brill; 2005.

108. Coope JA. Marriage, Kinship, and Islamic Law in Al-Andalus: Reflections on Pierre Guichard's Al-Andalus. Al-Masāq: Journal of the Medieval Mediterranean 2010; 20 (2): 171.

109. García-Sánchez E. Comidas de mujeres en la sociedad andalusí. In: Calero-Secall I, editor. Mujeres y sociedad islámica: una visión plural. Málaga: Universidad de Málaga; 2006. pp. 201–222.

110. Shatzmiller M. Labour in the Medieval Islamic World. Netherlands: Brill; 1994.