Emphasising the community: demographic composition of an exceptional tomb—the Chalcolithic burial site of Camino del Molino, Caravaca de la Cruz, Murcia

Reconstructing the biological profile of a skeletal sample is essential for defining a particular demographic group or classifying isolated remains. These results allow us to complete the population pyramid of a settlement, analyse mortality trends and relate individuals of a particular sex or age-at-death category to possible funerary rituals, lifestyles and/or states of health and disease. In this work, we carry out a paleodemographic analysis of a singular tomb: the chalcolithic burial site of Camino del Molino, Murcia, SE Spain. The tomb hosted 1348 individuals (30.7% non-adults and 69.3% adults) over two contiguous funerary phases, spanning a large part of the third millennium BC, which makes it a reference site for knowledge of the Recent Prehistoric populations. For this purpose, we estimated different paleodemographic parameters (life tables, mortality rates and sex ratios) and compared them to model life tables of preindustrial populations and data from other contemporary peninsular series to evaluate possible demographic anomalies. The results suggest that Camino del Molino was home to individuals of all ages and sex. However, there is a clear under-representation of newborns and nursing/breastfeeding infants and an over-representation of 5–15-year-old individuals. These findings could indicate potential issues related to diet/weaning, disease and early inclusion in the economic activities of the group.


Introduction
Paleodemography is a subdiscipline that attempts to identify the demographic parameters of populations that precede archaeological contexts from a data set of skeletal remains and historical documentary sources (Hoppa 2002: 9). The first studies on archaeological collections were developed in the late 1940s, although they were at their peak during the 1970s and 1980s (e.g., Angel 1947Angel , 1969Acsádi and Nemeskéri 1970;Masset 1973;Armelagos and Medina 1977;Bocquet-Appel 1977, 1986Bocquet-Appel and Masset 1977;Gage 1985). At the same time, a profound scientific debate emerged about the validity of these analyses in archaeological populations due to possible biases related to the nature of the reference samples used in the model life tables and age-at-death estimation methods (Bocquet-Appel and Masset 1982Masset , 1985Masset , 1996Van Gerven and Armelagos 1983;Buikstra and Konigsberg 1985). This discussion led to a methodological review of the discipline, with skeletal collections of which sex and age at death were already known, testing their validity on a larger number of archaeological samples. Subsequently, life tables were improved through the development of mathematical adjustments that offset the known biases (e.g., Gage 1988;Konigsberg and Frankenberg 2002;Paine and Boldsen 2002;Séguy and Buchet 2013).
Studies that attempt to reconstruct the population composition of prehistoric burials are becoming increasingly common. The Neolithic demographic transition has been one of the main topics in Recent Prehistoric population research (e.g., Piontek et al. 1996;Meiklejohn et al. 1997;Bocquet-Appel 2002, 2009Bocquet-Appel and Dubouloz 2003;Séguy and Buchet 2013;Dubouloz et al. 2017; 1 3 140 Page 2 of 16 Eshed and Gopher 2018). In contrast, recent studies from the Iberian Peninsula have mainly focused on skeletal collections from megalithic ossuaries (e.g., Silva 2003Silva , 2012Costa Caramé et al. 2010;Rivera 2011;Tomé and Silva 2013;Boaventura et al. 2014;Cunha et al. 2015;Silva and Ferreira 2007;de la Rúa 2015, 2016;Díaz-Navarro 2021), enriching the debate on selection trends by documenting the demographic anomalies of the individuals deposited there, with particular focus on age and sex.
It is important to take into account the limitations of the peninsular osteoarchaeological record of Recent Prehistory. Firstly, the fact that the skeletal remains are generally found in collective burials with disjointed remains hinders an accurate record and recovery. Furthermore, the poor bone preservation typical in this type of context and high fragmentation render it impossible to produce accurate estimates of sex and age. On the other hand, anthropologists often analyse past excavated collections, with little documentation and which are often poorly preserved as a result of old and obsolete storage protocols. Furthermore, until recently, human bones were rarely directly radiocarbon dated. The small size of prehistoric skeletal assemblages, coupled with their limited representation of the population as a whole, poses another challenge that further restricts our understanding of the communities that inhabited a specific time and place. This is directly linked to possible selection biases arising from mortuary practices. Finally, the temporal average of the skeletal groups that the prehistoric cemeteries comprise is another limitation.
On this basis, this research aims to contribute to the knowledge of the populations of the third millennium BC through the paleodemographic and chronometric analysis of a skeletal sample that is unique due to its volume, varied population composition, the excellent preservation of the skeletal remains and its recent excavation, with exhaustive documentation and recording methods. This is the collective burial of Camino del Molino (CMOL), Murcia, SE Spain, which contains the largest prehistoric funerary record known to date (Fig. 1a).

Material
This study analyses the full skull collection from the CMOL site, which was discovered by chance in December 2007 and excavated during 2008 (Lomba et al. 2009a(Lomba et al. , 2009b. The burial site was defined by a circular structure excavated from travertine with a diameter of c.6-7 m. The lower 2 m of this structure, where the archaeological deposit was located, was preserved intact (Haber et al. 2012;Lomba et al. 2009aLomba et al. , 2009b (Fig. 1b). Geomorphological analysis indicates that the structure was part of a cave, although a clear anthropic alteration can be noted both in the walls-which were regularised, adopting a bell-shaped form-and in the floor (Lomba et al. 2009b). This leads us to believe that it was originally an artificial cave or hypogeum with an entrance at the top and a perishable structure as a support, of which only the postholes are preserved .
The crania count carried out during the fieldwork indicated a minimum of 1300 individuals (Lomba et al. 2009a(Lomba et al. , 2009b)-a figure that has now increased after the study of the remains-making CMOL the largest cemetery in European Prehistory. The funerary sequence is characterised by the repeated introduction of corpses, in many cases simultaneously, which were grouped around the walls of the structure as the space filled up (Haber et al. 2012). This has led to the recovery of 167 complete and articulated individuals, mainly from the perimeter (Díaz-Navarro et al. 2023).

Radiocarbon dating and Bayesian analysis
A total of 28 14C dates from different articulated individuals (Table 1) are included. The main objective in selecting the samples was to have a well-balanced representation of the different stratigraphic units (SUU) and tomb areas (Fig. 2).
Bayesian inference was the analytical tool implemented to precisely define the phases of use and estimate their time boundaries with a high percentage of probability. The phase model, the overlapping phase model and the contiguous phase model were carried out, with an agreement of over 60%, for both the general models and for each of the data analysed (Supplementary Information S1).

Determination of the minimum number of individuals
Due to the large volume of skeletal remains, we conducted a quantitative analysis focusing on the most commonly found anatomical element in the sample, which was the skull. Only complete crania or fragmentary crania with preserved parietal, occipital, pars petrosa and frontal bones were taken into account.

Sex and age-at-death estimation
In non-adult individuals, the degree of dental eruption was analysed according to Ubelaker's (1989) chart method as the main criteria. Where teeth were not preserved, morphometric variables of different cranial regions such as the pars basilaris, pars squama and pars lateralis (Fazekas and Kosa 1978) were analysed, as well as the degree of ossification of temporal or sphenoid bones (Scheuer and Black 2000). Sex estimation in non-adults has been categorized as indeterminable, except for some individuals older than 16 years, where clear skull sex indicators were observed. The conservation of 82 articulated individuals with skull and pelvis of both sexes and all age categories has allowed us to establish the morphological and morphometric skull traits that best discriminate between sexes (Herrmann et al. 1990;Buikstra and Ubelaker 1994). For this purpose, 12 cranial and mandibular features were analysed, assigning each of them a sex estimation category (not preserved, female, male, indeterminate). The results were then compared with the most reliable estimates obtained from the pelvis, which were further validated through DNA analysis on 20 articulated skeletons. As a result, the most reliable indicators for sex estimation from isolated skulls were determined, as well as those that should be used with caution due to their wide margins of error ( Supplementary Information S2).
Likewise, the data on the degree of dental wear according to the method of Zoubov (1968) in 93 articulated skeletons whose age at death has been determined by other more reliable methods analysed in the pelvic Lovejoy et al. 1985) and costal (Işcan et al. 1985a(Işcan et al. , 1985b region were cross-checked. This has made it possible to establish a series of wear patterns that facilitate the ascription of the isolated skulls to an age-at-death category (Supplementary Information S2). The method designed by Zoubov (1968) offers a scale comprising five degrees of occlusal wear for incisors and canines and six for premolars and molars.

Paleodemography
There are a few limitations to be considered before commencing a paleodemographic analysis on an archaeological sample. These include (1) the size of the sample and how it represents the original population of the buried and the living (Bocquet-Appel and Masset 1977;Waldron 1994), (2) the stationary nature of the population (Sellier 1994) and (3) the existence of an archaic mortality pattern (Ledermann 1969).
In the CMOL case, although we accept the Waldron scheme (1994) and we are aware of the lack of information regarding the destruction of the north-east zone of the burial site by machinery (Lomba et al. 2009a(Lomba et al. , 2009b, we believe that we are dealing with a reliable and representative sample of the original living population, taking into account (1) the large volume of individuals, (2) the prolonged funerary use of the grave, (3) the complete excavation of the site with modern methods, (4) the population trends according to dates and the results of mobility analyses (Merner 2017) and (5) the presence of individuals of both sexes and all age categories, suggesting that there was apparently no selection process to access the grave.
From the estimates of the sex and age at-death of the crania, mortality tables have been produced based on stationary populations (Acsádi and Nemeskéri 1970) to determine the demographic composition of CMOL throughout the time of use of the grave. In addition, life expectancy (e 0 ) at birth was calculated for the overall sample, by sex and funerary use.
For this purpose, the age at death has been classified in 5-year intervals (Acsádi and Nemeskéri 1970) considering men and women both as a whole and separately. For the adult individuals, who were generally classified as young (20-39 years) or middle-aged adults (40-59 years), 5-year intervals have been established by cubic interpolation (Burden and Faires 1985) from the cumulative d x , while the relative d x was calculated from the representation of the interpolated cumulative d x (Valverde and Bush 1992). This is based on the premise that the distribution of a sample of archaeological adult individuals over short age intervals is currently impossible (Bocquentin 2003).
For individuals with a large standard error in age-atdeath estimation, located between two or more intervals, the distribution was made according to Sellier's (1994) principle of minimisation of demographic anomalies.
To detect demographic anomalies, several mortality coefficients have been calculated, which are useful to address problems arising from bias in adult age-at-death estimates and the frequent under-representation of non-adults in archaeological populations (Bocquet-Appel and Masset 1977;Bocquet-Appel 1979).
Finally, the sex ratio has been calculated to record irregularities related to the sex of the buried individuals (Ledermann 1969).
The demographic data has been compared to other contemporary archaeological series of the site under study, as well as to the Ledermann (1969) model life tables for preindustrial populations with a 30-year life expectancy, which closely resembles the estimated life expectancy of the CMOL population.

Chronometric analysis
The calibrated date graph suggests that the time span of this tomb covers almost the entirety of the third millennium BC (Table 1; Supplementary Information S1). However, the uncertainty intervals of each of the radiocarbon dates and the inaccuracy of the calibration curve in these chronological ranges make the temporal sequences appear much longer than they were.
Bayesian modelling shows that the contiguous phase model (Bayliss and Bronk Ramsey 2004) is the most robust as it has the highest agreement index (Amodel: 99.3) and more acceptable individual agreement indices (> 60%) (Supplementary Information S1).
The results of this analysis suggest that the first funerary phase of use started at the beginning of the third millennium (2971-2711 years cal BC, 2σ) and spanned a maximum of three centuries (0-259/82-377 years cal BC, 2σ). The second funerary phase was shorter (0-220/20-282 years cal BC, 2σ), and the end of the sequence is marked between 2451 and 2251 cal BC 2σ). The transition period between the two phases was 80 years at most (Fig. 3).

Paleodemography
In order to define the CMOL population, a mortality table has been prepared by dividing the total population into 5-year age groups, without taking into account the different levels of funerary use.
However, if we consider the mortality coefficient of the non-adult bulk ( 20 q 0 ), the CMOL figures are lower than those calculated for preindustrial populations (306.38‰ versus 557.5‰) (Fig. 5a).
Adults are the most numerous, accounting for 65% of the total population, and both subcategories (young and mature) have a similar number of individuals, with a slightly higher number of deaths in the mature age group (35% vs 31%). Old adults are the age category with the lowest mortality (3.7%) ( Table 2).
The sex ratio offers a value of 1.05:1, close to parity (1:1). Within the male sample, the highest number of deaths occurs in middle-aged adults (37%), surpassing the percentage of deaths among young adults (27%). In contrast, the higher percentage of women deaths occurs in youth (36% young adults compared to 29% of middle-aged adult females). Table 2 Life table for the CMOL population, which provides information on life expectancy at birth (e x ) n D x absolute frequency of individuals in the age range; n d x absolute frequency of individuals in the age interval among the total population per 100; n d xi relative frequency of individuals in the age interval after proportional distribution of indeterminates by cubic interpolation per 100 divided by the total population ( n d x ·100/ΣD x ); l n surviving population at an age interval. It is calculated by subtracting the sum of the percentages of the preceding age intervals from the total population (l n -1n d x -1); n q x probability of death for each age group obtained by dividing the percentage of deaths in an interval (d x ) by the number of survivors (l n ); n L x total years lived by the survivors of the age interval (n· (l n +l n +1)/2); n T x total years left to live in an age range ( n L x + n L x +1+ n L x +2+…+ n L x +n); e x life expectancy for each age interval ( n T x / n L x )  Furthermore, the sample has been divided by funerary phases (according to the archaeological results, the radiocarbon and Bayesian analysis) to determine possible differences in the demographic patterns of both phases of funerary use. The division of the sample allows us to observe a significant difference in the number of subjects deposited in both levels, documenting twice as many individuals in the second phase (Tables 3 and 4).
Looking at the life tables, we identify an increase in life expectancy of 2 years between the first and second level of use (Tables 3 and 4). Focusing on the 5-year age intervals, a similar distribution is documented in both phases, starting with a mortality peak among 5-9-year-olds that is more remarkable in the first phase (14.3%) than in the second (11.4%). From this category onwards, the figures decrease in both funerary levels until the young adult stage, with a mortality rate of 33.2% in the first phase and 29.7% in the second. The middle-aged adult's rate is 30.7% for the first level and 36.9% for the second. In the old adult's category, we did not observe any differences between the two phases (3.4% and 3.9%, respectively).
Mortality coefficients for both types of funerary use do not differ from those observed in the full sample ( Fig. 5b and Table 5). It should be noted that there is under-representation of perinatal and under-five individuals at both levels ( 1 q 0 = 4.26‰ and 9.1‰ vs. 288‰ of preindustrial populations (Ledermann 1969); 5 q 0 = 83.1‰ and 105.8‰ vs. 521‰; 4 q 1 = 78.9‰ and 96.7‰ vs. 306.5‰). However, mortality in this age range is higher in the second level of funerary use. Again, in both phases, the numbers increase considerably as the non-adults grow older. This is perceived in the 5 q 10 coefficient, documenting 142.85‰ in the first funerary level and 113.7‰ in the second, and likewise in the 5 q 15 coefficient, which yields values of 196.16‰ in the first funerary phase and 157‰ in the second, figures that are far from the parameters calculated by Ledermann (1969) ( Table 5). In this case, mortality is more pronounced in the first phase ( Fig. 5b and Table 5).
If we consider the mortality of individuals under 15 years of age at death and of all non-adults, the figures for both CMOL levels are lower than those established by Ledermann (1969) ( Table 5). A difference of 30‰ is observed between the two phases, with the first phase expressing a higher non-adult mortality rate (Fig. 5b).
In terms of sex distribution, there are no significant differences by funerary phases and, as in the general sample, parity is documented (sex ratio of the first level is 1.1:1 and that of the second is 1.03:1). However, a higher proportion of women dying in youth than in middle-aged is documented in both phases (64 mature women compared to 78 young women in the first phase and 117 mature women compared to 149 young women in the second phase, Fig. 6a, b). This contrasts with the pattern observed in male deaths, especially in the second funerary phase, with higher mortality after the age of 40 (75 young adults and 79 middle adult males in the first phase and 106 young adults vs 165 middle adult males in the second phase, Fig. 6a, b).

Discussion and conclusions
The CMOL collective burial housed the largest prehistoric ossuary known to date, with a minimum of 1348 individuals, deposited in two contiguous phases that span a large part of the third millennium (2971-2711 to 2451-2251 years cal. BC). The first burial phase housed the remains of 469 individuals, while during the somewhat shorter second phase, almost twice as many people (N 879) were buried in this grave. This seems to suggest a considerable demographic increase, supported by evidence of a decrease in infant mortality and a slight increase in life expectancy. Additionally, the results of the 87 Sr/ 86 Sr mobility analysis carried out on 93 individuals (Merner 2017) point to the presence of 12 non-local individuals at CMOL, 10 of them (83.3%) from the most recent burial phase.
To understand the magnitude of this tomb, it is necessary to investigate the funerary record of other contemporary sites. In Europe, although there are collective burials with a large number of individuals, such as the hypogeum of Crottes in Roaix with 136 individuals (Chambon 2003), or the hypogeum of Boileau in Vaucluse (Mahieu 1987) with 350, none come close to the size of CMOL. Other tombs with a large volume of human remains have also been identified in the Iberian Peninsula, such as the rock shelter of San Juan Ante Portam Latinam in Alava-MNI 388- (Etxeberría and Herrasti 2007), the tholos of Paimogo I in Lourinhã-MNI 413- (Silva 2012), the artificial cave of São Paulo II in Almada-MNI 255- (Silva 2012), the tholos of La Pijotilla in Badajoz-MNI 308-(Díaz-Zorita et al. 2017a), or the hypogeum of Can Martorell-MNI 194-(Mercadal et al. 2005) and La Sagrera-MNI 207- (Balaguer et al. 2015) in Barcelona.
The exceptional nature of the site lies not only in the size of the skeletal series, but also in the excellent conservation of the skeletal remains and the preservation of a large sample of articulated individuals. This is unusual for third millennium 1 q 0 (‰) 5 q 0 (‰) 4 q 1 (‰) 5 q 10 (‰) 5 q 15 (‰) 15 q 0 (‰) 20 q 0 (‰) CMOL first phase e° (0)  collective tombs, which generally have reduced or displaced primary deposits-in which the remains undergo continuous removal once they have been skeletonised-or secondary deposits. These aspects make CMOL a site of reference for studying the paleodemographic composition of the community deposited there. A series of demographic patterns can be traced in this work. A homogeneous population composed of individuals of both sexes and different ages was deposited at CMOL. However, anomalies have been identified in the frequencies of non-adults regarding the mortality pattern established for preindustrial populations. This is characterised by high infant mortality, most pronounced in the first 5 years of life, which decreases with increasing age (Masset 1973: 96-97).
Only 10 perinatal subjects were found at the study site, accounting for 0.7% of the total population. Therefore, their presence could be described as anecdotal. Likewise, there is an under-representation of subjects aged 1 to 4 years, three times below the parameters established in Ledermann's (1969) model life tables for preindustrial populations. This is especially significant at the first moment of use. Similar demographic patterns have been documented in other third millennium peninsular tombs (Botella 1973;Mercadal et al. 2005; Díaz-Zorita 2013; Fernández-Crespo and de la Rúa 2015; Díaz-Zorita et al. 2016, 2017a, 2017bEvangelista 2018;Herrero Corral 2019).
This demographic anomaly could be due to several factors. Firstly, we must consider a possible methodological bias in the fieldwork, related to an incorrect identification of the infant cranial remains. However, an anthropological specialist was present during the entire excavation, so methodological bias is unlikely. Another option to be considered is the differential preservation of the immature remains due to their more cartilaginous composition (Sellier 1994;Lewis 2006). The fact that a considerable number of skulls of this age category were preserved allows us to reject this possibility. We also ruled out a possible bias in age-ofdeath estimation methods, as this does not pose difficulties or large margins of error in non-adults when teeth are used.
The absence of infants in prehistoric burials has already been mentioned by other authors and has been related to cultural practices (e.g., Dedet et al. 1991Dedet et al. , 2022Gusi et al. 2008;Lancy 2014;Matney 2018), such as the possible exclusion of newborns and infants, whose social recognition would come later, or to the existence of a specific funerary ritual for this population group (Scott 1999;Delibes 2010;Fernández-Crespo and de la Rúa 2015;Burgos 2016;Blanco 2020). Traditionally, if a child died at an early age, they did not receive the same funeral treatment as other members of their community, renouncing an emotional bond with infants until they passed through critical developmental periods (Gennep 1981;Scott 1999). Lancy (2014) suggests that the infants' differential funerary practices must be related to the delay in the status recognition of these individuals and their acceptance as persons, so that they would not be integrated into the group until they reached a certain age. Similar patterns have been documented in other Iberian IV-III millennium graves, including the megalithic tombs of the Rioja Alavesa region (Fernández-Crespo y de la Rúa 2015), the necropolis of Humanejos (Madrid) (Herrero Corral 2019), the tholos I of Perdigões (Evora) (Evangelista 2018), the hypogeum of Costa de Can Martorell (Barcelona) (Mercadal et al. 2005) as well as southern sites such as La Pijotilla (Badajoz), Valencina-Castilleja (Seville), Panoría (Granada) and Los Millares and El Barranquete (Almeria) (Botella 1973;Peña 2011;Díaz-Zorita 2013;Díaz-Zorita et al. 2016, 2017b. On the other hand, there is a significant mortality peak in individuals of 5-15 years, with a higher incidence at the first funerary phase. This data, again, contrasts that estimated by Ledermann (1969) for preindustrial populations. This mortality peak is especially representative at 5-9 years, which is significant considering that mortality should decrease from the age of 5, when individuals are no longer as vulnerable (Masset 1973;Sellier 1994). This anomaly is not only found in CMOL but can also be traced in other collective burials in the Iberian Peninsula (Silva 2003;Mercadal et al. 2005;Cunha et al. 2015; Fernández-Crespo and de la Rúa 2015; Evangelista 2018; Herrero Corral 2019; Díaz-Navarro 2021).
The methodological approach applied renders it impossible for the over-representation of individuals aged 5-15 years to be the result of a compensation derived from a lack of such individuals in preceding groups, also ruling out a methodological bias in the age estimates for reasons previously explained for the preceding age group.
Therefore, this over-representation seems to allude, once again, to cultural practices. The analysis of the dentition of the skulls has allowed us to observe a large number of deaths between 3 and 6 years of age (36% of non-adults), which could be related to weaning. This could explain the higher mortality documented in individuals aged 5-6 years, which decreases considerably in the following age range (7-8 years, 4.1%) (Díaz-Navarro 2023). In most archaic societies, weaning occurs at 2-4 years of age (Nájera et al. 2010;Pearson et al. 2010;Fulminante 2015), and it has been shown to be a critical factor in an individual's life, with high mortality rates identified in the years immediately following (Schurr 1998;Pearson et al. 2010). This is because after breastfeeding, immunological benefits decrease and subjects are exposed to environmental pathogens and are therefore more likely to contract diseases (Katzenberg and Saunders 2008). Furthermore, a recent study on the Moro de Alins cave identified intermittent post-weaning stress episodes in individuals aged 4 to 7 years, suggesting that this period was as stressful or more stressful than weaning itself (Fernández-Crespo et al. 2022). Therefore, this dietary change could explain the over-representation of individuals in the 5-9-year category in both funerary phases. The paleopathological data also points to a special prevalence of metabolic indicators of non-specific origin in this age range, which seems to corroborate this hypothesis .
The mortality of older non-adults does not seem to be associated with the same factor as, although some authors have proposed that children may have had a different and/ or poorer diet than adults (Waterman 2012), analyses at certain sites show no significant differences in diet between age groups (Fernández-Crespo and Schulting 2017;Díaz del Río et al. 2017). In the absence of specific CMOL analyses, it is worth considering whether this demographic anomaly could be related to the community's rites of passage, i.e., the development of certain activities at an early age that endangered the lives of individuals, as well as their inclusion in the daily work of the community.
After weaning, non-adults would go from being protected in the family nucleus and directly linked to their mothers to acquiring freedom and exposure to the dangers of the outside world and the harsh conditions of physical labour (De Miguel 2010;Herrero Corral 2019). Different authors have suggested that children and adolescents would have been an important source of labour in multiple tasks such as sibling care, handicrafts (Lancy 2018), mining (Villalba et al. 1986), or livestock care, which could have led to increased exposure to animal pathogens (Waterman and Thomas 2011). In this regard, a high frequency of infectious diseases is documented in individuals of this age from CMOL, which could suggest close contact with livestock (Díaz-Navarro 2023).
These same patterns can be traced in present-day societies. For example, Lapp children receive their first leather belt and knife at the age of eight, when they are able to drive a reindeer to pasture, while Nigerian children take care of goats and camels (Lancy 2018: 90). Thus, the early development of physical labour or certain dangerous activities in children with still under-developed bodies could explain the higher mortality of 10-15-year-old individuals. We must also consider the impact of pubertal skeletal development. This process requires a greater demand for energy (body formation, growth spurts, hormonal processes and sexual development) which in CMOL was possibly not compensated for by a proper diet, leading to a growth retardation (Díaz-Navarro 2023) and, ultimately, a shortage of defences that exposed children, possibly already active, to a greater risk of contracting diseases (Cusick and Kuck 2012;Rogol et al. 2000).
The last 5-year category of non-adult individuals (15-19 years) exhibits the lowest mortality rates after individuals in their first year of life. These figures are similar to those observed in certain archaeological populations (González-Martín 2007;Livi-Bacci 2009;Fernández-Crespo and de la Rúa 2015;Díaz-Navarro 2021). Researchers suggest that individuals in this age range are at their optimal developmental stage and acquire maximum skills for the development of activities that in previous years could be lethal. However, these values differ from those observed by Ledermann (1969) in pre-Jennerian populations, which points to an increase in mortality from the 10-14-year age category to the subsequent group. This should be related to possible problems related to pregnancy and childbirth, as well as to the involvement of individuals in violent conflicts or risky activities that could lead to accidents. The contemporary samples to CMOL that report high juvenile mortality rates generally show numerous signs of violent deaths (Mercadal et al. 2005;Etxeberría and Herasti 2007;Rivera 2011). Even though the first possibility is plausible, we must also consider a potential methodological bias in the age-at-death estimation of skulls within 15-20-year range, which could explain the absence of this population age group. Determining an accurate age for individuals within this age becomes challenging once the process of tooth eruption and dental development is completed, especially if the measurement of open apices cannot be analysed, or if agenesis of the third molar is documented (AlQahtani et al. 2010).
The overall count of immature subjects calculated from the 20 q 0 coefficient is lower than that estimated for preindustrial populations, given the low proportion of individuals aged 0-4 years. Despite this, the total number of non-adults is 30%, which is within the limit of normal parameters for ancient populations (Bocquet-Appel and Masset 1977). This rate is higher in the first level of funerary use, suggesting a slight decrease in infant mortality in the community over the time of the grave's use, which could point to an improvement in living conditions.
The widespread presence of non-adults over 5 years of age in the grave shows that they were considered a member of the community, and therefore, it can be implied that, while the death of an infant affected only the most direct family sphere, the death of a child active in communal work had an impact on the whole community (Waterman and Thomas 2011).
The paleodemographic study also revealed that life expectancy at birth in CMOL was approximately 30 years. This age is at the limit of the values established for pre-Jennerian populations, which is estimated at 20-30 years (Ledermann 1969;Masset 1973). This is higher than that of earlier, contemporary, and later skeletal series, which we relate to the under-representation of non-adults in relation to the sample size, which clearly exceeds that estimated at other sites. In the megalithic ossuaries of the Duero Valley, the calculated life expectancy is 26 years (Díaz-Navarro 2021), while in sites from the fourth to third millennium BC transition, such as San Juan Ante Portam Latinam or Longar in Navarra, it is lower (20.3 years and 23.8 years, respectively) (Etxeberría and Herasti 2007;Rivera 2011), in accordance with the higher proportion of immature individuals compared to adults, which may be related to the violent acts previously mentioned. In Chalcolithic and Bronze Age tombs, this indicator varies from 20 to 27 years (Kunter 1990;Nájera et al. 2010;Pérez Villa 2013;Herrero Corral 2019).
The sex ratio allows us to rule out demographic anomalies related to the sex of the individuals, as it offers a value close to the parity found in natural populations (Ledermann 1969). When the sample is divided into funerary phases, a similar distribution of individuals of both sexes is observed, so there is no change in the conditions of access to the grave influenced by the sex of the individuals.
However, we documented differences in the age at death of males and females. When we cross-checked the sex and age-at-death estimations in both funerary phases, we observe that most males died in middle age, while females generally died in their youth. The higher mortality of females while young is related to obstetric risk, as childbirth is considered to be one of the most significant causes of mortality in preindustrial populations, especially the first birth (Pfeiffer et al. 2014). Complications in childbirth are generally related to anomalies between pelvic dimensions and baby size or infant positioning (Pfeiffer et al. 2014). The association of short maternal height with greater difficulties in childbirth is relevant (Guzmán et al. 2001;Toh-Adam et al. 2012;Stulp et al. 2011), given that CMOL women are significantly short (148 cm) (Díaz-Navarro 2023). Furthermore, some authors argue that plasticity in maternal skeletal growth may lead to a reduction in pelvic size in response to poor health and nutrition during this period (Wells 2012), which is also noticeable among CMOL women, as they show a much higher morbidity status than men in the community (Díaz-Navarro 2023).
There is an imbalance between the sexes in favour of males in the tombs dated to the transition from the fourth to the third millennium in the north of the Iberian Peninsula (Etxeberría and Herasti 2007;Fernández-Crespo and de la Rúa 2015). In contrast, the sex ratio values are inverted in caves from the same geographical area such as Las Yurdinas II, Peña Larga and Peña de Marañón, which leads researchers to suggest the existence of burial practices differentiated by sex (Fernández-Crespo and de la Rúa 2016). There is also a slight male superiority in Chalcolithic tombs from the south of the peninsula (Jiménez-Brobeil 1988;Nielsen et al. 2010;Díaz-Zorita 2013;Díaz-Zorita et al. 2016), while in other sites in the south and west, the values are reversed (Silva 2012;Pecero 2016;Beck 2016).
The sexual composition of the contemporary peninsular tombs presents a heterogeneous panorama, although there are no patterns of selective access related to the sex of the individuals that are visible in previous periods (Fernández-Crespo and de la Rúa 2015, 2016; Díaz-Navarro 2021). The socioeconomic and political transformation observed in the transition from the fourth to the third millennium appears to be linked to a change in the rules of access to collective tombs, which is reflected in homogeneous population compositions in Chalcolithic tombs. This is especially visible in CMOL, where individuals of both sexes and all age categories are included, with the exception of clearly underrepresented children under 5 years of age, a pattern well documented in prehistoric and protohistoric societies.

Acknowledgements
The contribution of Sonia Díaz-Navarro in this study was supported by a Predoctoral Fellowship awarded by the Castilla & León Regional Government and European Social Fund (ORDEN EDU/574/2018). The contribution of Cristina Tejedor-Rodríguez was supported by a Postdoctoral Fellowship awarded by Juan de la Cierva Programme-Incorporation Grants, of the Ministry of Science and Innovation, Government of Spain. Some of the unpublished radiocarbon dates presented in this work were funded by the University of Valladolid Research Project "METOO" (PROYEMER-2021-45).
Author contributions All authors contributed to the study conception and design. Material preparation, data collection and analysis were performed by SDN, MHU and CTR. The archaeological data were provided by JLM, director of the excavation of the tomb analysed. The first draft of the manuscript was written by SDN and all authors commented on previous versions of the manuscript. All authors read and approved the final manuscript.
Funding Open Access funding provided thanks to the CRUE-CSIC agreement with Springer Nature.

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/.