Dietary reconstruction of the Bronze Age necropolis of Cova des Pas (Minorca Island): evidence from δ13C and δ15N analyses

Current paleodietary studies about the Naviform and Talayotic groups that took place in Minorca (Balearic Islands) during the Late Bronze Age–Early Iron Age (ca.1600–850/800 BC) have suggested a mixed and variable diet, largely depending on terrestrial sources of vegetables and meat. This study explores the nutritional pattern of the individuals buried in the Cova des Pas site (Minorca Island, Spain), a cave used as a collective sepulcher and the most exceptional and major human assemblage found in the Balearic Islands during this period. Carbon (δ13C) and nitrogen (δ15N) stable isotopic signatures were measured on extracted bone collagen from 49 individuals. Further, faunal remains from the Son Mercer de Baix site, the closest contemporaneous village to the collective sepulcher, were also analyzed to provide a baseline corpus of data to interpret human isotopic data. The results indicate a human diet based mainly on C3 plants with an important consumption of animal protein. The δ15N values of infants up to 4 years were high, indicating the enriched isotopic signature of breast milk, and weaning is assumed to occur around this age. Differences between sexes and the age subcategories were not statistically significant, assuming that the different groups of society had the same access to food. The data obtained in this isotopic study provides insight into the palaeodietary pattern of the human groups dated to the Late Bronze Age–Early Iron Age ages in the island of Minorca, contributing to the present debate on the emergence and development of complex societies on the Balearic archipelago.


Introduction
Diet is one of the most important factors to understand ancient human populations, which is not only a means of survival but a cultural and social phenomenon as it reflects the socioeconomic dynamics involved in subsistence strategies and access to food consumption (Gumerman 1997). The necessity of archaeological research to better approach issues of economy and ideology in antiquity quickly led to the study of food, its origins, its production, and its consumption. The analysis of stable isotopes from bones offers a direct way to reconstitute an organism's diet and hence is a very optimal way to overcome the current limitations of traditional osteological approaches (Tykot 2004;Knudson & Stojanowski 2008).
During the Bronze Age (ca.1600-850/800 cal BC), the Balearics witnessed the emergence and development of complex societies, including the Naviform culture (documented in the whole archipelago) and, later during the Early Iron Age, the Talayotic culture (ca.1100-600 cal BC, exclusively in Majorca and Minorca) (see Lull et al. 1999;Guerrero Ayuso et al. 2002;Anglada et al. 2014;Depalmas 2014). The dietary pattern of these periods corresponds to a controversial issue. The traditional hypothesis that had been launched was that its subsistence strategy focused on animal husbandry (Cerdá, 1978;Roselló-Bordoy 1968), with plant cultivation remaining of minor significance. However, bioarchaeological studies dealing with zooarchaeological, archaeobotanical, and anthropological remains (Hernández-Gasch et al. 2002;Pérez-Jordà et al. 2018) pointed out a mixed farming and herding economy, without a systematical use of marine or freshwater food items in Majorca and Minorca. Contrarily, in the island of Formentera, C 4 plants such as millet and marine resources, including fish and mollusks, were apparently consumed by prehistoric human groups (Sureda et al. 2017a). First studies using the isotopic δ 13 C and δ 15 N measurements from human and animal bone samples measured for radiocarbon dates ( 14 C) suggested similar results (Van Strydonck et al. 2005). Palaeodietary studies using elemental analysis on human bones recovered from the necropolis of the Cova des Càrritx (1400-800 cal BC) in Minorca show that contribution on diet was proportional between vegetables and meat (Pérez-Pérez et al. 1999). Multi-element analysis on human and animal bones of the neighboring necropolis of S'Illot des Porros (ca. 1400-850 cal. BC) revealed a mixed diet with a high consumption of marine food, seeds, and cereals (Subirà & Malgosa 1992a). Dental microwear studies of a Talayotic population such as the Son Real (ca. 800-600 cal. BC) population (Majorca, Spain) found also a mixed diet, but depending on abrasive plant foods and a reduced amount of meat intake (Jarošová et al. 2006). Very interesting in those studies was that archeobiochemical data presented no nutritional difference between males and females in Cova des Càrritx (Pérez-Pérez et al. 1999) and S'Illot des Porros (Subirà & Malgosa 1992b), while an increased vegetal or cereal food intake by females in comparison to males is hypothesized in Son Real (Jarošová et al. 2006).
The Cova des Pas (CdP) is one of the most important prehistoric sites in Minorca and the Balearic archipelago ( Fig. 1). During its excavation, a unique collection of burials was recovered. Although other massive contemporaneous burial deposits have been documented up today in this area (e.g., Son Olivaret, Cova des Càrritx), this deposit allows the major representation of complete individuals . The possibility to study each individual as a whole allowed performing a variety of different diet studies, such as dental or boning analysis, and comparing between them. Furthermore, all the paleoanthropological information available from previous studies (Fullola et al. 2008;Armentano et al. 2010;Van Strydonck et al. 2010;Malgosa et al. 2011;Cañas Cortés 2012;Simón et al. 2016) permits a better approach and comparative analysis of the samples of the site and witnesses its uniqueness.
CdP is a small burial cave located on the left side of the Trebalúger ravine, near the town of Ferreries in Minorca (Balearic Islands, Spain). It is a small cavity, approximately 6.5 m wide and 4.5 m long, situated in the wall of the ravine about 15 m above the ground (Fullola et al. 2008) (again see Fig. 1). It was used as a collective sepulcher during the Late Bronze Age (1100-800 BC) by the Talayotic culture of Minorca (Bergadà et al 2015). The cave was discovered by three speleologists, Pere Arnau, Pep Riera, and Monica Zubillaga during the spring of 2005. During the excavation works in 2005-2006, the site was completely dug up, and an important set of burials in the cave was recovered (Fullola et al. 2007;Armentano et al 2010).
The environmental conditions of the cave allowed an exceptional state of preservation of the skeletal remains (Armentano et al. 2010). Anthropological studies identify a total set of 66 different individuals (Fig. 2). The burial shows a picture of age (35 adults and 31 subadults) and sex structure (25 women, 26 men, and 15 more of unknown sex) highly representative of the population of the island more than 3000 years ago, associated with the Talayotic period (Armentano et al 2010;Simón et al. 2016). Health status was also reconstructed revealing an extensive record of diseases that include other pathologies, from congenital alterations to tumors, fractures, inflammatory alterations, or injuries due to activity . On the taphonomic level, the studies carried out have been able to recognize a complex treatment of the corpse (Cañas Cortés 2012) and the conditioned burial space itself. The great relevance of this burial chamber lies in the relative preservation of human organic tissues, compatible with muscles, brain, and lungs, for instance (Armentano et al. 2012;Prats-Muñoz et al 2013). At the genetic level, the results indicate the existence of important endogamy at the level of mitochondrial DNA, which could indicate a closed and isolated population (Simon et al., 2016).
Son Mercer de Baix (SMB) is an open-air settlement located on the edge of the Fideu ravine (Ferreries, Minorca Island). It is, possibly, one of the most emblematic naviform villages in the Balearic Islands, due to the state of conservation of its Naviform I, or "Cova des Moro" (Fig. 3). The village is composed of, at least, 4 naviform structures (boat-shaped houses with cyclopean walls) and 2 other rectangular rooms connected to them (rooms 1 and 2). It was initially excavated by Mª Lluïsa Serra Belabre in 1962 (Naviform II) and later in the 1980s by Cristina Rita and Lluís Plantalamor (Naviforms III, IV, and the connected rooms) (Plantalamor & Rita 1983;Rita 1982, 1987, Rita and Topp 1988. Finally, in 2001, conservation and consolidation actions were performed at Naviform I, including the excavation of some exterior areas. Most of those materials and archaeological contexts are still unpublished, although recently, some research related to their naviform structures or metallurgical remains have been started (Garcia Amengual 2006;Sureda 2022).
Balearic Naviform villages are known to have a large occupation sequence that could last several centuries as the cases of Closos de Can Gaià (Majorca) (Servera 2005) or Cap de Barbaria II (Formentera) (Sureda et al. 2017b). SMB presents a first Bronze Age (Early Naviform) occupation phase (ca. 1600-1200/1100 cal BC), represented by the naviform buildings but also by several archaeological materials, including typical Naviform ceramic and metallic assemblages (globular bowls, truncated conical cups, quadrangular chisel, riveted knife, bronze bracelet, or various crucibles).
Furthermore, a specific Late Bronze Age/Early Iron Age (ca. 1200/1100-800 cal BC) occupation phase can be inferred in the site, re-using some of the old buildings (Naviforms II, III, IV, rooms 1 and 2) and generating new archaeological features and sedimentary deposits, especially  -Muñoz et al. (2012) represented by the presence ceramic sets of large barrels and large pots composed by globular shape and flat base, known to be "index fossils" for Balearic LBA (Guerrero, 2008). Besides, other Minorcan naviforms such as Cala Morell or Son Blanc have provided similar archaeological contexts with radiocarbon dates (ca. 1200-900 cal BC) associated with this set of materials (see SI Table 1).
The faunal remains that were found at SMB include Ovis sp., Capra sp., and a considerable number of Sus sp. No signs of agriculture were initially detected (Fernández Miranda, 1997), but this should be reconsidered as several grinding stones were recovered from the site, today hosted at the Museu de Menorca (Minorca, Balearic Islands; Spain). The SMB represents an appropriate context to be used as an isotopic referent for CdP cave. It is the closest contemporaneous site largely excavated in the area (less than 600 m distanced from the CdP and over 3 km from the sea), and it represents a domestic habitat co-existing with the CdP funerary cave.
The present research aims to study the dietary pattern of the prehistoric population buried in the CdP site, investigating the type of protein contribution to the diet and the degree of exploitation of marine resources for human consumption, as well as the potential differences in access to food resources among population groups (considering sex and age of sampled individuals). For this purpose, direct information will be obtained from the isotopic composition of the major number of individuals by measuring the isotopic signal of carbon (δ 13 C) and nitrogen (δ 15 N) on the skeletal remains. At the same time, this study will investigate a corpus of faunal remains recovered from the SMB open-air site in order to interpret the isotopic data retrieved from the anthropological collection from CdP. Due to the absence of reliable chronological data from SMB, we also present here two new radiocarbon data to assure the chronological utility of CdP and SMB selected samples.
The results of all the previous studies carried out on this site, providing information on the population groups (i.e., age, sex, funeral rite, deposit, health status), render this study exceptional and allow much deeper project research. In addition, isotopic values of faunal remains were also analyzed to provide baseline values of the local trophic food chain and to explore whether domesticated animals had access to marine foods.

Stable isotopes and palaeodietary reconstruction
Carbon (δ 13 C) and nitrogen (δ 15 N) stable isotope ratio analyses of the bone collagen are the method most used for the reconstruction of dietary protein patterns in archaeological populations (Schoeninger 1995;Katzenberg 2000).
Collagen isotopic values reflect the diet of the last 10 years, approximately, before the death of the organism being studied. This is because in every living organism, the bone tissue is continuously remodeled by the interaction of osteoblasts and osteoclasts, cells responsible for this process (Hedges et al. 2007). For the biochemical reconstitution of the human diet, carbon and nitrogen stable isotopes are the most used isotopic markers (Katzenberg 2000). The basic principle of these approaches is that "we are what we eat" (Kohn 1999); that is, the basic units that make up all body tissues of any animal, including bones, come from the foods they have ingested during their lifetime. Several experimental studies have shown that as an organism processes the amino acids, it has consumed to create the protein it needs, and a difference in the carbon isotope ratio between the dietary protein and collagen is created by ~ 4 to 5‰ from bulk diet to collagen consumer (Ayliffe et al 2004;Ambrose & Norr 1993;Schoeninger & DeNiro 1984) and ~ 1‰ from collagen diet to collagen consumer in a mono-isotopic diet (Bocherens & Drucker 2003). However, the nitrogen isotopic fractionation increases by ~ 3-5‰ for each trophic level (Minagawa & Wada 1984).
Stable carbon isotope values are able to distinguish the presence of plants with different photosynthetic pathways (i.e., C 3 and C 4 plants) in human diet since the ranges of their δ 13 C values do not overlap each other (Deines 1980;Van der Merwe & Vogel 1978): C 3 plants from temperate and cold areas with values of about -26.5‰ such as wheat, and C 4 plants from tropical, dry, or semi-dry regions with values around -13 ± 2 ‰ such as millet Winter 1981;Kohn 2010). Moreover, carbon isotopes can also detect a diet based on terrestrial or marine foods with values approximately at − 20‰ and -12 ± 1‰, respectively (Schwarcz & Schoeninger 1991;Richards & Hedges 1999). Something similar occurs with the stable nitrogen isotopes. The δ 15 N composition of human bone can determine the trophic level position on the food chain, with the values Fig. 3 Son Mercer de Baix and its landscape: a SMB naviforms and associated structures distribution (modified from Plantalamor and Rita 1983). b, d, and e Different views of the "Cova del Moro" Naviform in SMB. c General view of fhe Fideu's ravine where SMB is located Table 1 List of human specimens sampled for collagen analysis Type of bone, sex and age of sampled specimens, collagen carbon (%C) and nitrogen (%N) contents, collagen carbon (δ 13 C) and nitrogen (δ 15 N) isotope compositions, and collagen carbon:nitrogen atomic ratio (C:N). The samples that were excluded from further discussion due to their poor collagen preservation are presented with an asterisk increasing by 3-5‰ with each ascending step of the food chain (Bocherens & Drucker 2003;O'Connell et al. 2012).
The local baseline background of the trophic food chain (faunal and floral remains) is required to build a more reliable interpretation of the nitrogen isotope values in the human diet (O'Connell et al. 2012). In the same time, carbon and nitrogen isotopic values in human and animal tissues could also reflect environmental factors affecting the soil-plant system at the base of the food chain, such as temperature, precipitation, aridity, salinity, soil composition, or manuring, as well as physiologic and metabolic factors (Szpak 2014;Reitsema 2013).

Cova des Pas (CdP)
A total of 49 human specimens were measured by stable isotope analysis. The rest of the samples recovered from the site did not proceed to analysis due to the lack of enough yields during the extraction procedure. Samples were selected considering previous paleoanthropological data from Armentano et al. (2010) Table 1). The sampling procedure was carried out by careful selection of rib bone fragments for almost all individuals, except in five cases (samples 2-5, 3, 33, 36, 58, 64, and 70) where rib bones were not available and then skull bone or phalanx or fibula was selected, depending on the bone available for each specimen. Ribs have a relatively fast cortical turnover rate, depicting an individual's diet from a more recent period prior to death (Cox & Sealy 1997;Fahy et al. 2017).

Son Mercer de Baix (SMB)
Faunal samples for isotopic analyses belonged from rooms 1 (level III) and 2 (sector 1, Capa 1) and UE 1008. These archaeological contexts are stratigraphically linked to the Late Bronze Age/Early Iron Age (ca. 1200/1100-800 cal BC) occupation of the site. Before selection, faunal samples were identified and classified following zooarchaeological criteria. The osteological analysis focused on the study of the taxonomic and body part representation, age-at-death, and potential pathologies observed in order to obtain as much information as possible about these animals when they were alive. Sex information was not possible to record due to the high fragmentation of the bones that prevented observing possible diagnostic traits (such as in the pelvis) or taking measurements for comparison. Ovis sp. and Capra sp. samples were classified following osteological criteria in Boessneck (1980), Payne (1985), and Prummel and Frisch (1986). The faunal set is compound of 33 samples (Table 2) and comprised samples of Caprinae (Capra hircus/Ovis aries) (n = 11), Bos Taurus (n = 5), Capra hircus (n = 6), Ovis aries (n = 7), and Sus domesticus (n = 4). Further, two faunal samples recovered from room 1 in SMB were selected for radiocarbon dating and measured at the CNA-CSIC research center (Spain). The first of these (SMB03) corresponded to an upper Ovis sp./Capra sp. molar, while the second (SMB06) was a second phalanx of Sus sp.

Bone collagen extraction
Sample preparation and collagen extractions were performed at the biomolecular laboratory at the Catalan Institute of Human Palaeoecology and Social Evolution (IPHES) held in Tarragona . This protocol allows working with very small samples of bones while ensuring optimal extraction of fractions of collagen fibers, including those samples that by a diagenetic or fossilization process preserve very little quantity. A small bone fragment of each specimen was mechanically cleaned using a Dremmel rotating tool equipped with a circular diamond-coated blade to remove all visible contaminants. Briefly, the bone collagen was purified according to the acid-base-acid protocol proposed by Longin (1971) and modified by Bocherens et al. (1991). Bone shards (ca. 300 to 350 mg) were soaked in 0.5 M HCl for demineralization, in NaOH (0.125 M) to remove contaminants, rinsed with distilled water several times, and gelatinized with 0.01 M HCl at 100 °C for 17 h (pH 2). Once filtered after the gelatinization process and frozen, samples were freeze-dried at the Institute of Chemical Research in Catalonia (ICIQ). Gelatin-collagen samples weighing about 300 μg were analyzed using a Thermo Flash 1112 elemental analyzer (EA) coupled to a Thermo Delta V Advantage isotope ratio mass spectrometer (IRMS) with a Conflo III interface, at the Institute of Environmental Science and Technology (ICTA), Universitat Autònoma de Barcelona (Barcelona, Spain). The international standard laboratory IAEA 600 (caffeine) was used as a control. Analytical precision was checked using the international standard laboratory IAEA 600 (caffeine), and the average analytical error was < 0.15‰ (1σ) calculated for each isotopic measurement, δ 13 C and δ 15 N separately. The standard utilized for δ 13 C was Vienna PeeDee Belemnite (V-PDB), and the standard used for δ 15 N was air N 2 (AIR). The isotope ratios are expressed for carbon as δ 13 C vs. Vienna Pee Dee Belemnite (V-PDB) and for nitrogen as δ 15 N vs. atmospheric nitrogen (AIR): δ χ¼ (R sample /R standard −1 ) × 1000‰, where χ stands for 13 C or 15 N and R stands for 13 C/ 12 C or 15 N/ 14 N. The reliability of the isotopic signatures of the collagen extracts was addressed using several criteria (yield of extraction ≥ 10 mg·g −1 , percentages of C ≥ 30% and N ≥ 10%, and the atomic C/N ratio 2.9 < C/N > 3.6).
When evaluating the differences in the stable isotope results between samples, statistical significance was tested using the T test and the ANOVA and MANCOVA tests, setting the significance p value < 0.05. Nonparametric tests were used when the data did not follow a normal distribution. As for the age, the samples were grouped as two different subcategories, one including the individuals over 13 years old and the other those with age under 13 years old.

Selection of radiocarbon dates from Son Mercer de Baix
Two faunal samples from Son Mercer de Baix were selected for radiocarbon dating. These samples proceed from room 1 and level 4, collected during the 1978 excavation campaign. It is expected that these samples correspond to the oldest archaeological context documented in this area. Type of bone, species identification, age of sampled specimens, collagen carbon (%C) and nitrogen (%N) contents, collagen carbon (δ 13 C) and nitrogen (δ 15 N) isotope compositions, and collagen carbon:nitrogen atomic ratio (C:N) (Pr = probably)

Radiocarbon dates from Son Mercer de Baix
From the results of both dates, it can be proposed that room I of SMB was in use, at least, between 1643 and 1412 cal BC (Table 3). In the same way, it is interesting to note that the construction of rooms I and II as well as that of the Naviform IV, to which they have adhered, would necessarily be prior to 1504 cal BC and its abandonment, necessarily posterior to 1514 cal BC. Also, these dates define a terminus post quem for all the selected fauna remains which are necessarily younger than those levels.

Collagen preservation
The results of the stable isotopes and collagen quality indicators are reported in Tables 1 and 2. In the faunal samples, yields were ranging from 25 to 159 mg/g. In the human samples, all the selected bones yielded collagen which was ranging from 6 to 169 mg/g. Five samples 5 (24, 44, 61, 62, and 68) show slightly low C or N% values lower than these limits but showed C:N ratios ranging within the internationally accepted limits and were considered for discussion. However, nine samples (6, 26, 33, 34, 50, 55, 58, 64, and 66) presented both C and N contents, and C:N ratios out of the accepted limits were discarded for final interpretation. Excluding these nine samples, the mean carbon and nitrogen contents (%C and %N) of the bone collagen extracts were 35.7 ± 4.64% and 12.8 ± 1.68%, respectively. The mean atomic C:N ratio of samples used for interpretation is 3.2 ± 0.09. These results are indicative of collagen suitable for isotopic analysis supporting markers of good collagen preservation (DeNiro 1985;Ambrose 1990).

Human isotope results
Results from the isotope analysis of bone collagen from human specimens are presented in Table 4. The mean of carbon and nitrogen isotope values are − 19.7 ± 0.34‰ and 9.5 ± 0.91‰, respectively. The carbon isotope results range from − 20.5 to − 18.9‰, while those of nitrogen isotope range from 8.3 to 12.4‰. The δ 13 C and δ 15 N data for adult sample are plotted in Fig. 4. 95CIs and mean values are presented in Fig. 5, where sexes do not show statistically significant differences. Considering both isotopic signals separately, as well as, jointly, there were not significant differences between sexes in both cases (δ 13 C: T test: 0.228; p = 0.822; δ 15 N: T test: 1.306; p = 0.205; and MANCOVA test: Wilks' Lambda = 0.903, p = 0.342).
The δ 13 C and δ 15 N data for all samples between different age categories are plotted in Fig. 6. The differences between the three subcategories (16 to 25 years old, 25 to 45 years old, and over 45 years old) for the individuals  (Fig. 7). The δ 15 N results show that three out of six specimens (39, 49, and 65), classified in the age category under 4 years old, present high values (12.4‰, 12‰, and 11.4‰, respectively). This observation is shown in Fig. 8, where δ 15 N values are distributed only for individuals under 13 years old. In this case, the graph shows that young specimens under 4 years old have high δ 15 N signatures and that this value is decreasing on sampled specimens between 0 and 13 years old.

The local faunal baseline
The isotopic analysis of faunal samples offers a more reliable interpretation of the human isotopic values. In this study, considering the mean and ranging δ 13 C values (− 20.6 ± 0.72‰; max − 19.1, min − 21.8‰) of all faunal samples, it can be demonstrating a livestock diet based on C 3 plants. However, it is noteworthy the enrichment of the carbon isotope values of Capra sp. (samples 37.1, 43.1, 43.7, 43.8, 39.2, 43.6). This could be explained by the broad range of variation in C 3 plants, since this plant group includes most of the species known, nearly all trees, herbs, and cold climate grasses. Further, C 3 plants can also vary due to a wide range of environmental effects like radiation, water availability, soil nutrients, temperature, or salinity, among others (O'Leary, 1995;Farquhar et al. 1989;Tieszen and Boutton 1989;Ehleringer and Monson 1993;Heaton, 1999), and plants adapted to dry environments tend to show greater water use efficiency (hence, higher δ 13 C). δ 15 N values in faunal adult samples provided a mean value of 6.5 ± 0.60‰ and maximum and minimum values of 7.5‰ and 5.3‰, respectively. However, it is interesting the fact that juvenile Caprinae (samples 32.2, 42.1, 42.4, 42.5, 42.6, 42.7) show high values (max. value 9.6‰). Considering their low δ 13 C value (from − 20.3 to − 21.8), the values of these specimens could be explained by the breastfeeding effect. Moreover, some subadult animal specimens (samples 2.1, 39.1, 1.1, 41.1, 43.5, 44.1, and 44.3) show slightly enriched δ 15 N values (7.1‰, 7‰, 7.1‰, 6.9‰, 7.1‰, 7‰, 7.4‰, and 7.5‰, respectively) in relation to the adult samples. These values could indicate the weaning effect and the introduction to the adult diet. However, their difference is not statistically significant and strong enough to make a powerful conclusion (Kruskal-Wallis: x 2 = 3.51; p = 0.06).

Human dietary reconstruction
The isotopic results of the human samples from the CdP site indicate a main diet based on terrestrial resources. Considering data from faunal remains, the mean δ 13 C value of the bone collagen in adult samples (− 19.7‰ ± 0.31‰) fits well with a consumption mainly based on C 3 plants compounds in their diet, while the mean δ 15 N value (9.2‰ ± 0.57‰) will suggest a significant contribution of animal protein in the human diet (e.g., Lanting and van der Plicht 1996a, b).
These values in the case of the samples 32, 48, 51, and 69 (− 19.4‰, 19.3‰, 19.4‰, and − 19.4‰ respectively), accompanied by low δ 15 N values (8.5‰, 8.7‰, 8.9‰, and 8.6‰ respectively) not expected in marine ecosystems, could suggest a minor proportion of C 4 plants in the diet of those specimens, although current uncertainties on bulk collagen isotope must be taken into account. Because C 4 plants are naturally scarce or absent in Minorca Island today and in the past (Mateu, 1993;Pyankov et al 2010), a simple explanation will suggest consumption of domestic C 4 plants on the diet of Talayotic individuals buried in the CdP site. Millet is the early C 4 plant introduced in the Iberian Peninsula, but its introduction on the human diet during the Middle and Late Bronze Age is still an open debate (see López-Costas et al. 2015, for instance), and its introduction in the Balearic Islands is only suggested for the Naviform period in Formentera (Sureda et al. 2017a), while its presence in the other islands is still unknown.
Alternatively, some of these enriched δ 13 C values (e.g., sample 27; − 19.2‰) could be the consequence of another carbon source diet. Indeed, these adult specimens like 3, 8, 18, and 25 show also higher δ 15 N values in collagen (≥ 10‰), potentially revealing the influence of marine food products on their diet.
At this point, it must be stressed that previous δ 13 C and δ 15 N analyses (Van Strydonck et al. 2005) clearly indicated that the diet in the Balearic Islands during the prehistory was not predominantly based on marine resources. This lack of a clear evidence of consumption of marine foods has also observed in other studies from Majorca (Davis 2002;Garcia et al. 2004) and the Mediterranean region in general (Craig et al. 2006), all of them during prehistoric times, although that does not mean that human groups did not consume a small amount of it, undetectable by current bulk collagen isotope analysis. Comparing the adult samples of the CdP with other contemporaneous samples (1100-800 BC) (Fig. 9) from Minorca and Majorca Islands (Fig. 10), it can be observed that the samples from the CdP share, basically, the same mixed diet of animal products and vegetables as the other contemporaneous sites (Van Strydonck et al. 2005) (data about samples used is provided in Table C1 of SI).
However, considering both isotopic signals jointly, differences between the different sites were statistically significant (MANCOVA test: Lambda de Wilks = 0.331, p < 0.001). When evaluating the values of δ 13 C and δ 15 N separately, it seems that these differences are statistically significant for both δ 13 C and δ 15 N values (ANOVA test: δ 13 C: F = 4.57; p < 0.001; δ 15 N: F = 3.34; p = 0.002), where CdP samples tend to show higher δ 15 N but lower δ 13 C values than the majority of the other sites. It is interesting the fact that CdP shows considerably lower values than those obtained in the two human samples from the close Ses Arenes site from Minorca (Fig. 10). Unfortunately, although the number of sites is important to elaborate this comparison, the limited number of samples per site and the quality of the information provided for each sample (i.e., unknown data about sex and age of specimens sampled) make it difficult to reach more deep interpretations.
The age information about the samples of Son Blanc and Binipati Nou is unknown (Van Strydonck et al. 2005), but taking into account that they show the same high δ 15 N values as the infant samples of CdP, and it can be assumed that they share the same range of age. Otherwise, it could indicate the consumption of an unknown terrestrial dietary item or a diet containing juvenile animals not fully weaned and thus N-enriched (Jay & Richards 2006).

Other socioeconomic glints of light
Differences on diet between males and females are not statistically significant, the fact that does not allow concluding for variation between sexes in the access to food. The scatter plot in Fig. 4 confirms that males and females up to date measured share close values in both δ 13 C and δ 15 N results. These results agree with those of previous studies in different sites of Minorca such as the Cova des Cárritx (Pérez-Pérez et al. 1999) or the Binipati, Es Tudons, and southern Rafal Rubí (Van Strydonck et al. 2002).
This can be explained by the general trend that the Naviform and Talayotic societies-both in Majorca and Minorca Islands-consisted of egalitarian communities (e.g., Gasull et al. 1984;Lull et al. 1999) with no sign point the existence of a hierarchical society or constituted by units of population differentiated within a settlement by its social function or its economic resources (Fernández-Miranda 1991).
However, Anglada et al. (2014) proposed that some of the human communities that inhabited Minorca between the late 2nd and early 1st millennium cal BC were immersed in an incipient process of increasing complexity and maybe social inequality, suggesting that a kind of an elite group was connected to the control of storage and redistribution of agricultural products. This seems not to relate to the population of the CdP, since differences in the access to food do not appear to exist either between sexes or the different age categories. Paleodemographic and taphonomical studies of the CdP indicating the same life expectancy and burial treatment between sexes (Armentano et al. 2010; confirm this social equality. The nitrogen isotopic signal of the infant samples #37 (2 years ± 8 months old), #38 (1.5 years old), #39 (3 years old), #49 (3-4 years old), and #65 (1,5-2 years old), 9.8‰, 10‰, 12.4‰, 12‰, and 11.4‰, respectively, was highly indicating the enriched isotopic signature of breast milk (see again Fig. 8). Fogel et al. (1989), attempting to evaluate nursing and weaning, demonstrated that infants feeding on mothers' milk exhibit enrichment δ 15 N values and a decrease to levels similar to those of mothers shortly after breastfeeding stopped. Taking into account that samples #25 (9-10 years old), #48 (7 ± 1 years old), and #57 (8 ± 1 years old) do not follow the same pattern as #37, #38, #39, #49, and #65, and it can be assumed that breastfeeding was occurring until around the age of 4 years, at least. A major number of samples are required to investigate this point to exclude other reasons for these values, such as potentially ill or malnourished infants whose isotope values could have been influenced by metabolic processes or not enough isotopic information from infants under and above the age of 4 years old. However, these data could add some light on the duration of the breastfeeding period and weaning process of the prehistoric population of the Cova des Pas.

Conclusions
The foregoing analysis uses the stable isotope of carbon (δ 13 C) and nitrogen (δ 15 N) in order to make inferences about the diet pattern of the individuals buried in CdP site. The isotopic data from human samples show a human diet based largely on C 3 plants with an important contribution of meat revealing a mixed diet of plants and animals. The δ 15 N values of infants up to 4 years were high, indicating the enriched isotopic signature of breast milk, while it is hypothesized that weaning occurred around this age in this population. Differences between males and females were not statistically significant, assuming that both men and women had the same access to food. The same pattern seems to occur between the age subcategories. The paleodietary data obtained here is important in broadening the debate on the emergence and development of complex societies in the Balearic archipelago, providing direct information about the way how these societies managed diet habits and how insularity conditions affected the way they were socially organized.
Acknowledgements This study was supported by the Institut Menorquí d'Estudis (IME) from the Consell Insular de Menorca (Balearic Islands, Spain) through a research grants called "Estudi isotòpic del sepulcre funerari de la Cova des Pas: reconstrucció del patró d'alimentació de Menorca fa 3000 anys" led by CT, AM, and TS and "La producció metal·lúrgica menorquina en època Naviforme. Una aproximació arqueometal·lúrgica" led by PS. A preliminary version of this study was presented as Final Project of TS at the Master in Biological Anthropology (Autonomous University of Barcelona, Spain). Research at IPHES is framed in the CERCA program. This research was supported by the Spanish Ministry of Science and Innovation through the "María de Maeztu" excellence accreditation (CEX2019-000945-M). CT was supported by the Spanish Ministry of Science and Innovation through the "Ramon y Cajal" program (RYC2020-029404-I). AM and XJ are members of the GREAB supported by Generalitat de Catalunya (Ref. 2017SGR 1630. PS was supported by the Xunta de Galicia-GAIN (IN606C-2021/002) and the Spanish Ministry of Science and  (Table 1 of Supplementary Information) and comprised samples from Cova des Pas (n = 40), Binipati Nou (n = 1), Rafal Rubí (n = 6), Son Blanc (n = 1), Cala Morell (n = 1), Cova Gregoria B (n = 1), Es Tudons (n = 1), Ses Arenes (n = 2), Ses Aritges (n = 5), Ses Roques Llises (n = 4), and Alcaïdus (n = 2). The appropriate bars for standard deviation have also been included Innovation through the "Juan de la Cierva-Incorporación" programme (IJC2020-042714-I/ AEI/https:// doi. org/ 10. 13039/ 50110 00110 33). We want to thank Damià Ramis for providing some of the SMB pictures.
Funding Open Access Funding provided by Universitat Autonoma de Barcelona.

Data Availability
The authors confirm that the data supporting the findings of this study are available within the article and its supplementary materials.

Declarations
Competing interests The authors declare no competing interests.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http:// creat iveco mmons. org/ licen ses/ by/4. 0/.