Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Anthracological study of a Chalcolithic funerary deposit from Perdigões (Alentejo, Portugal): A new analytical methodology to establish the wood burning temperature

  • Ginevra Coradeschi ,

    Roles Conceptualization, Data curation, Formal analysis, Investigation, Methodology, Resources, Validation, Writing – original draft, Writing – review & editing

    ginevrac@uevora.pt

    Affiliations HERCULES Laboratory, University of Évora, Évora, Portugal, CityUMacau Chair in Sustainable Heritage, Évora University, Évora, Portugal

  • Nicasio T. Jiménez Morillo,

    Roles Data curation, Methodology, Software, Validation, Writing – review & editing

    Affiliations MED – Mediterranean Institute for Agriculture, Environment and Development & Department of Biology, School of Science and Technology, Mitra (UE Campus), Évora, Portugal, Institute of Natural Resources and Agrobiology of Seville (IRNAS-CSIC), Seville, Spain

  • Cristina Barrocas Dias,

    Roles Supervision, Writing – review & editing

    Affiliations HERCULES Laboratory, University of Évora, Évora, Portugal, Chemistry and Biochemistry Department, School of Sciences and Technology, University of Évora, Évora, Portugal

  • Massimo Beltrame,

    Roles Writing – review & editing

    Affiliations HERCULES Laboratory, University of Évora, Évora, Portugal, CityUMacau Chair in Sustainable Heritage, Évora University, Évora, Portugal

  • Anabela D. F. Belo,

    Roles Investigation, Writing – review & editing

    Affiliation MED – Mediterranean Institute for Agriculture, Environment and Development & Department of Biology, School of Science and Technology, Mitra (UE Campus), Évora, Portugal

  • Arturo J. P. Granged,

    Roles Methodology

    Affiliations Department of Applied Physics I, University of Sevilla, Seville, Spain, Med_Soil Research Group, Faculty of Chemistry University of Sevilla, Seville, Spain

  • Laura Sadori,

    Roles Supervision

    Affiliation Department of Environmental Biology, Sapienza University of Rome, Rome, Italy

  • António Valera

    Roles Writing – review & editing

    Affiliations Archaeological Research Unit of Era Arqueologia S.A., Lisboa, Portugal, ICArEHB Algarve University, Gambelas Campus, Algarve University, Algarve, Portugal

Abstract

Anthracological analyses of charcoal samples retrieved from Pit 16 of Perdigões (Reguengos de Monsaraz, Portugal), a secondary deposition of cremated human remains dated back to the middle of the 3rd millennium BC, enabled the identification of 7 different taxa: Olea europaea, Quercus spp. (evergreen), Pinus pinaster, Fraxinus cf. angustifolia, Arbutus unedo, Cistus sp. and Fabaceae. All taxa are characteristic of both deciduous and evergreen Mediterranean vegetation, and this data might indicate that the gathering of woods employed for the human cremation/s occurred either on site, or in its vicinity. However, considering both the large distribution of the identified taxa and data about human mobility, it is not possible to conclusively determine the origin of the wood used in the cremation(s). Chemometric analysis were carried out to estimate the absolute burning temperature of woods employed for the human cremation/s. An in-lab charcoal reference collection was created by burning sound wood samples of the three main taxa identified from Pit 16, Olea europaea var. sylvestris, Quercus suber (evergreen type) and Pinus pinaster, at temperatures between 350 and 600 °C. The archaeological charcoal samples and the charcoal reference collection were chemically characterized by using mid-infrared (MIR) spectroscopy in the 1800–400 cm-1 range, and Partial Least Squares (PLS) regression method was used to build calibration models to predict the absolute combustion temperature of the archaeological woods. Results showed successful PLS forecasting of burn temperature for each taxon (significant (P <0.05) cross validation coefficients). The anthracological and chemometric analysis evidenced differences between the taxa coming from the two stratigraphic units within the Pit, SUs 72 and 74, suggesting that they may come from two different pyres or two different depositional moments.

1. Introduction and research aims

Perdigões is a large set of ditched enclosures, nowadays classified as a National Monument, located in Reguengos de Monsaraz municipality, Évora district, in the inner Alentejo region, centre-southern Portugal (Fig 1A and 1B). It holds at least 16 ditches and a megalithic cromlech, occupying an area of ca. 18 ha. The archaeological exploration of the area started back in 1997 and the site had a long and complex diachrony, extending from 3400 BCE (Middle—Late Neolithic transition) to 2000 BCE (Early Bronze Age).

thumbnail
Fig 1.

A. Location map of Perdigões archaeological site in the Iberian Peninsula, B) Aerial photo of Perdigões ditch enclosures.

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

Interpreted as a ceremonial and aggregation centre, these ditched enclosures have several funerary structures, revealing diversified funerary practices in terms of architecture, location within the site, body treatment and votive materials [13]. The archaeological material under study was recovered in Pit 16 (middle 3rd BC), located in the central area of the enclosures [2, 3]. Several stratigraphic units (SU) were identified during the archaeological excavation, and those that contain human remains were designated SU 72 and 74 (Fig 2).

thumbnail
Fig 2. Section of Pit 16.

In dark grey: the two stratigraphic units (SUs 72 and 74) containing human remains associated with charcoals and the grave goods.

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

In southern Portugal and more generally throughout the Iberian Peninsula, there are no other Chalcolithic funerary contexts with the characteristics of Pit 16. Burial rituals in this area during this period were quite heterogenous, including dolmens, tholoi, pits, ditches and hypogea [4, 5]. Cremation seems to have been used only for a few individuals per site and wasn’t the most common practice [6]. The different ways of human bodies processing, and the type of funerary structures and grave goods differ between cremation and non-cremation funerary contexts. Some authors [7] argued that these differences may have been related to different social and cultural backgrounds and to social emulation processes.

Perdigões, along with other recent archaeological discoveries, led to a revision of the Megalithic literature [1, 2, 411] and, at present, a more dynamic picture of the Iberian Chalcolithic funerary world is available. The reinterpretation of the world of death within this particular chronology has fostered the production of a number of studies on anthropological and archaeozoological remains recovered in these contexts, and the cultural materials associated with them [1, 4, 1223].

Overall there is a general lack of studies of funerary charcoals across Europe, with only a few studies published about Prehistoric sites [2427] with examples for Great Britain and Ireland [2830]; starting from the Roman Period, studies become more frequent and systematically conducted [3139].

In the literature a charcoal assemblage deriving from human cremation deposits is considered a single or a short time span episode [30]. The main goal of the study of these types of archaeological charcoals is to get paleo-ethnographic information, wherein the relationship between humans and their surrounding nature is investigated [25]. The analysis of charcoals from these types of funerary contexts may reflect rational and less rational human choices, linked with a ritual character of the site, where the taxonomic selection of woods employed for the pyre may have a cultural meaning [25, 30, 34, 35, 40, 41]. Indirectly, the interpretation of human cremation contexts is tied to the study and to the interpretation of the charcoal assemblage. On the other hand, as some authors have already mentioned [30], charcoal is the remainder of the funeral pyre, being the most direct way to get a glimpse on the ritual of body cremation.

The aims of this study are:

  • To characterize the different wood taxa from Pit 16.
  • To evaluate the absolute burn temperature of the tree main wood taxa identified from the Pit.
  • To obtain information regarding this unusual Chalcolithic funerary context, exploring the existence of one or different moments of secondary deposition of human bodies.

Two different approaches have been used to achieve these goals. The first one includes the anthracological analysis of the charcoal assemblage, and the second one comprises the FT-MIR/chemometric analyses of archaeological charcoals and sound wood samples (of the same species identified) combusted at different temperature (i.e., formed under controlled temperature and atmosphere using a muffle furnace in the laboratory).

We consider this approach a good option to successfully evaluate the charcoal assemblage from Pit 16. Results were useful for a better understanding of the archaeological context and the related ritual/s.

1.1. Checking out different techniques for the study of the archaeological charcoal absolute burning temperature

The partial combustion process of wood either in a forest fire, in a bonfire or in a funerary pyre is well known. It normally produces charred material which contains a significant amount of refractory carbon (C) [42, 43]. This is commonly known as “black carbon” or charcoal [44]. This material is defined as a continuum of thermally altered products, ranging from slightly charred still degradable biomass, to highly condensed, refractory soot [45, 46]. The charcoal particles are chemically heterogeneous, but typically they are made up of a conspicuous amount of unspecific chemical aromatic structures [46].

The concentration of aromatic compounds within the chemical structure of the charcoal particles is significantly associated with direct and indirect factors. Direct factors are linked to combustion characteristics, including temperature, fire permanence and heating rate (°C/min). On the other hand, indirect factors are connected to the fuel (wood) properties, like the type of wood, size, wood resilience and moisture [4749]. Thus, chemical changes in the continuum charcoal structure can be identified, and are influenced by the combustion temperature and the molecular composition of the pyrogenic material [5052].

Currently, different techniques are available for the charred organic matter evaluation of complex matrices such as soils and sediments, as well as to assess their burning temperatures. One of these techniques is the reflectance, which is considered a valuable tool to estimate the burn temperatures of modern and archaeological charcoal using a Reflected Light Microscope [49, 53]. Another analytical technique is the Fourier Transform Mid-Infrared (FT-MIR) spectroscopy, which is micro-destructive and employed for the chemical characterization of complex matrices in archaeological samples, like wood charcoal particles and bone collagen [5456].

The combination of FT-MIR spectroscopy and multivariate analytical techniques, such as Partial Least Squared (PLS) regression [57], has demonstrated to be a powerful tool for the evaluation and prediction of combustion temperature of organic matter and the composition of pyrogenic carbon in soils [44, 58]. Recently, some researchers have studied burned wood samples in terms of the Variable Importance Projection (VIP), which has allowed them to assess MIR bands significantly linked with dependent variables [44, 58].

2. Materials

2.1. Archaeological context of the charcoal samples

Pit 16 is interpreted as a secondary funerary collective deposition, a singular funerary context amongst the Iberian Chalcolithic. The pit contained human remains of at least 9 individuals (6 adults and 3 sub-adults), mixed with other burnt artefacts (faunal remains, fragments of pottery, ivory idols, and arrowheads) as well as a large number of charcoals. The pit is circular, with a profile shrinking at the top, presenting a diameter of 0.92 m at the top and 1.50 m at the base, being 0.90 m deep (Fig 2). In total, nine different stratigraphic units (SU) were identified, corresponding to three phases of deposition. The first phase corresponded to the initial filling of the pit, where two thin layers (SU126 and SU125) were formed by some residual materials, namely pottery sherds. The second phase corresponded to the secondary deposition of the cremated remains. A conic shaped deposit was formed, indicating the dumping of the remains into the pit (probably from containers). This deposition was subdivided into the two aforementioned stratigraphic units, SU72 and SU74, the later one being less compact and dustier. Both sediments were dark grey with a lot of charcoal and ashes, involving an abundant amount of human remains (5000 fragments), some cremated faunal remains and burned archaeological materials (pottery sherds, ivory schematic figurines, a copper awl, stone arrowheads) [2, 3, 5961]. Nevertheless, whether the two stratigraphic units represent a single deposition moment or not, is not clear. The third phase corresponds to the final filling of the pit, with the depositions of several layers with no human remains, incorporating mainly faunal bones and pottery sherds, most of them burned (SU76, SU71, SU90, SU68, SU58).

2.2. Current vegetation of the area

Nowadays, in natural areas, the woody vegetation cover is made up mostly of evergreen oaks (Quercus suber, Quercus rotundifolia, Quercus coccifera), Arbutus unedo, some Pinus spp. and Olea europaea var. sylvestris, along with many sclerophyll shrubs, like Cistus spp., Rosmarinus officinalis, Lavandula spp., Myrtus communis, Pistacia lentiscus or Phillyrea sp.; some deciduous trees are also present near streams and other watercourses, like Fraxinus angustifolia, Populus sp., Salix sp..

3. Methods

3.1. Sampling

Pit 16 sampling was performed based on the excavation pursued, which divided the pit into several stratigraphic units. For this work, sediments from the SU72 and the SU74 (Fig 2), corresponding to the secondary deposition of cremated human remains, were entirely sampled.

In the field, archaeologists dry sieved the sediments from SU72 and SU74 using a 2 mm mesh size to isolate human remains [60]. After the bone recovery, sediments were preserved for the anthracological study.

3.2. Anthracological analysis

Sediments from SU72 and SU74 (which had been previously dry sieved in the field to recover human remains) were additionally treated and observed in the laboratory using a stereo-zoom microscope (model LEICA M205C equipped with a camera) to remove sand grains, small ceramic, and animal bone remains and to isolate charcoal fragments (i.e., larger than 2 mm) to proceed with the anthracological analysis. A total of 754 charcoals were recovered.

The charcoal fragments were manually fractured to obtain the three anatomical sections (transversal, tangential and radial) under a stereo-zoom microscope (LEICA M205C equipped with a camera). The charcoal’s anatomical features were determined using a reflected light optical microscope (LEICA DM2500 equipped with a camera) at different magnifications (50–1000x). A selection of charcoal fragments was analysed with a Scanning Electron Microscope (SEM) to obtain high-resolution images. Wood atlases [6265] were used as comparative tools for the charcoal identification, together with an in-lab reference collection of wood specimens. Some other types of observation were also made. When present, fractures in radial directions (radial cracks) [6671], vitrification, consisting in a variable fusion of some anatomical parts of the wood (which make the wood structure homogeneous) [6974], and radial groves in the inner part of tracheids cells (known as reaction wood) [69] were recorded.

With respect to the evergreen oak’s wood samples identification, it is significant to reiterate that it was not possible to make any distinction between the evergreen oak tree species. In Portugal, nowadays, as well as in the Mediterranean region, this group includes holm oak (Quercus ilex L.), kermes oak (Quercus coccifera L.) and cork oak (Quercus suber L.) [64, 75].

Regarding the identification of the olive wood samples, the cultivated species (Olea europaea var. europaea) is not differentiable, at the anatomical level, from the wild one (Olea europaea var. sylvestris). However, studies concerning the beginning of the cultivation of this plant in the south of the Iberian Peninsula suggest that it started only in the Roman period [76]. Considering the investigation period, charcoal fragments of Olea from Pit 16 could probably be ascribed to the wild variety, O. europaea var. sylvestris.

Concerning the identification of the maritime pine wood, archaeobotanical remains showed the presence of this species since the Holocene in the northern Iberian Peninsula, in the mountainous regions of eastern Iberia as well as in the southwestern part of the peninsula [7779]. This proves the naturalness of this species in these areas [77, 80].

In this work, the identification of charcoals belonging Pinus pinaster was based on the presence and distribution of teeth in the radial tracheids, taking into consideration three levels of degree of dentation [64, 8183].

The identification of Fabaceae wood is considered particularly difficult due the large intraspecific structural variability present in this group [8486]. Few species of this group were anatomically described [65, 84]. In the present study charcoal samples belonging to this taxon were identified only at the family level.

3.3. Data quantification

For data interpretation, taxa values are interpreted using the relative frequency (%) and the mean values (%U) obtained using the “ubiquity correction” method [87]. The mean value (%U) of a specific taxon (A, for example) on n different stratigraphic units (SU) is calculated as: %U = %A (SU1) + %A (SU2) + %A (SUn) / n. This method performs a correction for ubiquity of the relative frequencies of each taxon. Charcoal size classes (0.3–0.99 cm, 1–1.99 am, and 2–2.99 cm) for each wood taxon identified were also collected.

3.4. Estimation of the absolute burn temperature of archaeological charcoals

3.4.1. Experimental combustion tests on modern reference material.

To create a statistical model able to estimate the absolute combustion temperature of archaeological charcoals an in-lab charcoal reference collection was created, and different experimental combustion tests were carried out. Sound wood samples of the three main taxa identified from Pit 16 with the anthracological analysis, namely olive, evergreen oak, and maritime pine woods have been sampled in the proximity of the archaeological site (Fig 3).

thumbnail
Fig 3. Sound woods transversal anatomical section of Olea europaea var. sylvestris (A), Quercus suber (B), and Pinus pinaster (C) heated at 350°C.

Samples were photographed with SEM.

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

Regarding the olive and considering the investigation period (as already reported), wooden samples belonging to the wild variety, Olea europaea var. sylvestris has been chosen as reference material.

For the evergreen oaks, samples belonging to Quercus suber have been chosen as reference material, considering the abundance of this species in the Alentejo region [88]. After sampling, sound wood samples were cut (dimension: 1 cm3, ≈ 3 g) using an automatic saw, then washed with distilled water and dried for 24 h. Samples were then positioned into ceramic crucibles, heated in the absence of oxygen, between 350 and 600°C, the typical burning temperature of wood in natural fires [89]. In total, 6 different heating steps were used, increasing the muffle furnace temperature of 50°C each time.

Samples remained within furnace during 90 min at maximum temperature. For each heating step, two replicates of each taxon were burnt. At the end of the heating period, a sample of each laboratory-based combusted charcoal was photographed by SEM after the gold coating procedure.

3.4.2. FT-MIR spectroscopy analysis of modern and archaeological charcoals.

An aliquot from the two charred replicas of each modern wood sample was selected for Fourier Transform Mid-Infrared analysis. The two charred aliquots of each taxon were powdered using an agata mortar and combined to create one composite sample for each different wood taxon.

Two aliquots of each archaeological sample of Olea europaea, Quercus spp. (evergreen) and Pinus pinaster were detached for FT-MIR analysis. As for modern samples, the two aliquots of each archaeological taxon were powdered and then combined to create one composite sample to be analysed.

FT-MIR spectra were obtained using an Alpha -R Spectrometer (Bruker Optics®, Germany), with Attenuated Total Reflectance (ATR) module (Platinum-ATR-sampling module, Bruker, Germany), at a wavelength range of 4000–400 cm-1 and a resolution of 2 cm-1. To improve the signal to noise ratio, 60 spectra were co-added and averaged for each recorded spectrum. Spectral data were background corrected to a reference spectrum obtained prior to every measurement, and some spurious absorptions, such as peaks from atmospheric CO2, could be removed.

3.4.3. Statistical analysis.

Partial Least Squares (PLS) regression models were generated using the ParLeS software [90]. This method has been used to acquire predictive cremation temperature models from the FT-MIR spectral intensities in the range 1800–400 cm-1 (independent variables, 250 data points). The spectral pre-processing treatments consisted of a light scatter and baseline correction by Standard Noise Variate (SNV), de-noising through a median filter, and mean centering [58]. This pre-processing has been widely used [44, 58, 91, 92] to obtain forecasting models of different environmental and physical variables, using different analytical techniques (FT-IR, Py-CSIA, Py-GC/MS and FTICR/MS) on different natural matrices (e.g., soil, sediment, foodstuff). The protocol of this pre-processing is well exposed in [58]. The Root Mean Squared Error (RMSE) and the Akaike’s Information Criterion (AIC) were used to prevent overfitting and determine the best number of factors (latent variables) for each model. Lastly, the diagnostic spectral regions of the FT-MIR spectra were studied by the combination plot of the Variable Importance for Projection (VIP) values in the spectral range under study. The VIP traces can be useful to identify the independent variables (spectral peaks), that may be linked to the temperature of the cremation. Spurious forecast models due to overfitting were discarded after comparing PLS models calculated with the randomised combustion temperature (cross-validation). The cross-validation model was previously used by [44, 58].

4. Results

4.1. Anthracological analysis

After the analysis of 754 charcoal fragments from Pit 16, 7 taxa were determined at the species level, Olea europaea, Pinus pinaster, Fraxinus cf. angustifolia, Arbutus unedo, at the section level, Quercus spp. (evergreen), at the genus level, Cistus sp., and at the family level, Fabaceae. Few charcoal fragments were identified only at the division level, Angiosperm (Table 1, Figs 4 and 5).

thumbnail
Fig 4. O. europaea, Quercus spp. (evergreen) and P. pinaster archaeological charcoal samples observed and photographed with SEM.

(A, B) Transversal and radial sections of an O. europaea sample, (C, D) Transversal and tangential sections of a Quercus spp. (evergreen) sample, (E, F) Transversal and radial sections of a P. pinaster sample.

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

thumbnail
Fig 5. Cistus sp., A. unedo, F. cf. angustifolia and Fabaceae archaeological charcoal samples observed and photographed with SEM.

(A, B) Transversal and tangential sections of a Cistus sp. sample, (C, D) Transversal and tangential sections of an A. unedo sample, (E) Transversal section of a F. cf. angustifolia sample, (F) Transversal section of a Fabaceae sample.

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

thumbnail
Table 1. Absolute number (N), relative frequency (%) and relative frequency corrected for ubiquity (%U) of taxa found in Pit 16.

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

Three taxa are dominant in the anthracological record of Pit 16, featuring the 83.8% and the 87,3%U of total of identified samples: olive tree (O. europaea), evergreen oak tree (Quercus spp. (evergreen)), and maritime pine tree (P. pinaster) (Fig 4). The most abundant taxa are olive (32.9%, 27.6%U) and evergreen oak (30.2%, 30.7%U). Finally, the maritime pine is the third taxon in abundance, accounting the 20.7% and the 29%U of the totality (Table 1, Fig 6A). From the SU72, 208 charcoal fragments were recovered, and 5 different taxa were identified, while from the SU74, a total of 546 charcoal fragments were recovered and 7 different taxa were identified. Moreover, the amount of the 3 main taxa is not similar between the two stratigraphic units, with the maritime pine tree being the main taxon of the SU72, while in SU74, olive tree is predominant (Table 1, Fig 6B).

thumbnail
Fig 6. A) Relative frequency of taxa found in Pit 16, B) Relative frequency of taxa found in the SU72 and in the SU74.

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

These differences become evident by the application of the “ubiquity correction” method, which accounts not only for the absolute number of charcoal fragments, but also for the taxa distribution within the samples (Table 1).

Regarding charcoal fragments of Quercus spp. (evergreen), samples from both SUs 74 and 72 showed the presence of vitrification. This feature was observed in 85 (37.2%) fragments of this taxon. For the P. pinaster samples from SU 72, fractures in radial directions (radial cracks) were recorded for a total of 90 (90%) charcoal fragments. Moreover, for 27 charcoal fragments of this species, radial groves in the inner part of tracheids were recorded.

With reference to charcoal size, the olive tree is the species with the highest quantity of the largest charcoal, with 71 fragments with values between 1 and 1.99 cm and 9 fragments between 2 and 2.99 cm., followed by the maritime pine and by the evergreen oak, for which the largest charcoals collected are represented by only 14 fragments (Fig 7).

The remainder of arboreal taxa it is represented by few fragments of narrow-leafed ash wood (F. cf. angustifolia). This taxon was recorded only in the SU74, and its presence does not exceed the 1.3% and the 0,9%U (Table 1, Fig 5E). Shrubs were also identified in the anthracological record of Pit 16: rock rose wood (Cistus sp.) legume wood (Fabaceae) and strawberry tree wood (A. unedo), accounting for the 8.4% and the 7,4%U of the total (Fig 5A–5D and 5F). The most present taxon is the rock rose (5.4%, 4,9%U), identified on the genus level as Cistus sp. (Table 1, Fig 5A–5D). Between shrubs, two fragments of A. unedo and three of Fabaceae have been identified as twigs, having retained their shape. The values of charcoal size showed that all charcoal fragments belonging to these species were basically smaller, with values between 0.3 and 0.99 cm if compared with the three main taxa (Fig 7). Finally, a fava bean seed (Vicia faba) was identified in SU74. It was not possible to identify the Angiosperm charcoal fragments beyond the division level due to their state of conservation.

4.2. Estimation of the absolute burn temperature

4.2.1. FT-MIR spectroscopy of archaeological samples.

Fig 8 displays the infrared spectra, in the region of 1800–400 cm-1, of the archaeological charcoal samples from SU74 (Fig 8A) and SU72 (Fig 8B) belonging to each of the three main taxa identified after anthracological analysis of archaeological charcoals. In all spectra, different bands were observed, namely the 1570 cm-1 band related to the amide groups, the 1380 cm-1 band linked to C-O stretching of phenolic OH, the 1730 cm-1 band related to stretch of groups C = O, and the bands 1270 and 1040 cm-1 belonging to lignin. These results were in line with those observed by other authors [93, 94] for burnt organic matter samples. The band at 1570 cm-1 does not have an unequivocal assignation [95], nevertheless it may be referred to the C = C stretching in aromatic and quinone structures, that is highly conjugated to carbonyl group [96]. This band is typical of carbonized, fossilized or highly humified organic matter [97]. The band at 1380 cm-1 may be related to unspecific aliphatic structures, as well as hydroaromatic compounds.

thumbnail
Fig 8. This figure represents the overlay of the infrared spectra to show possible differences among replicates.

Medium infrared spectrum (MIR) of archaeological charcoal of O. europaea (A), Quercus spp. (evergreen) (B) and P. pinaster (C) retrieved from a) SU74 and b) SU72.

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

4.2.2. FT-MIR spectroscopy and PLS analysis of modern charcoals.

The cross-validation plots (predicted vs observed values) for the burning temperature of each set of modern wood charred taxa are displayed in Fig 9. The PLS models, using exclusively the information contained in a discrete infrared spectra range (1800–400 cm-1, 250 data points), successfully predicted the crematory temperature of each taxon. Comparison with an alternative model derived from the fully randomised dependent variable that had no connection (p > 0.05) with infrared data validated the model’s accuracy. The PLS model of Olea showed the highest significant prediction (r2 = 0.998, p = 0.0000), following by Quercus and Pinus (r2 = 0.987, p = 0.001 and r2 = 0.957, p = 0.007, respectively). In all cases, the RMSE and AIC values suggested the use of four factors for all taxa. Fig 10 shows the VIP values for the different taxa used to obtain the different PLS models. The superposition of the average infrared FT-MIR spectrum of each taxon with the VIP trace indicates the best diagnostic bands to predict the combustion temperature of each taxon. The most significant FT-MIR bands to predict the burning temperature in the case of Olea model were those positioned at 1730, 1580, 1410, 800 and 740 cm-1, while for the Quercus model, the predominant bands were located at 1580, 1410, 1370, 800 and 740 cm-1. In the case of Pinus model, the main bands were identified at 1730, 1410, 1370, 1245 and 1040 cm-1.

thumbnail
Fig 9. Cross-validation plots for crematory temperature of modern charred woods taxa: A) O. europaea var. sylvestris, B) Q. suber (evergreen type) and C) P. pinaster, generated by partial least squares (PLS) regression using infrared spectral as predictors.

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

thumbnail
Fig 10. Superposition of variable importance for projection (VIP) of spectral points and average infrared spectra of modern char particles of O. europaea var. sylvestris (A), Q. suber (B) and P. pinaster (C).

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

Each forecasting model was made up of characteristic infrared bands [44]. In the case of the Olea samples, the statistic model used the bands at 1730 and 1580 cm-1 corresponding to vibrating C = C and C = O bonds, respectively, in the aromatic rings [96]. In addition, the bands at 1410 cm-1, related to C–C stretch (aromatic ring) of lignin compounds [98] and 1370 cm-1, linked to methyl C–H rock in alkene compounds [99] were also used by the forecasting model of Olea. This taxon also used the bands at 800 and 740 cm-1 which are typical of polysaccharide compounds [100]. The Quercus PLS model used the same FT-MIR bands as the model of Olea; however, in this case, the band at 1730 cm-1 was not determinant. On the other hand, the model of Pinus employed the band at 1245 cm-1 linked to acidic compounds of the resin [101]. This attribution is further supported by two characteristics bands at 800 and 740 cm-1, which indicate the presence of resinous compounds [102]. In addition, in this model, the band at 1040 cm-1, that is related to the C–O stretch, displayed a noticeable high VIP value.

4.2.3. Predicted absolute burn temperature of archaeological charcoals.

The predicted absolute burn temperature for the archaeological charcoal samples identified in SU74 and SU72 are displayed in Table 2, respectively. The validation of the forecast model by cross-validation process supported its use to determine the burn temperature experimented by archaeological charcoal samples. The temperature value is obtained from the models developed using charred modern wood taxa. In SU74, the Pinus pinaster showed the highest combustion temperature (576 ± 24°C) in comparison to Olea europaea (544 ± 20°C) and Quercus spp. (evergreen) (509 ± 6°C). The one-way ANOVA analysis showed that temperatures were statistically significant (95% confidence interval) among different wood taxa (Table 2). However, in SU72 (Table 2), the significant (p< 0.05) highest crematory temperature is registered by O. europaea (506 ± 12°C), while Quercus spp. (evergreen) and P. pinaster (492 ± 3°C and 487 ± 3°C, respectively) showed no significant (p>0.05) differences among them.

thumbnail
Table 2. Mean, minimum and maximum temperature (°C) values of the three main taxa from SU74 and SU72.

One-way ANOVA, the different letters indicate significant (p < 0.05) mean temperature (°C) differences between taxa according to the Tukey test.

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

5. Discussion

5.1. Different wood taxa

All taxa identified are characteristic of the Mediterranean pine forest and sub-Mediterranean sclerophyllous vegetation [103].

It was possible to identify three main taxa in the anthracological record of Pit 16: O. europaea, Quercus spp. (evergreen) and P. pinaster (Fig 3), with olive and evergreen oak woods prevailing (Table 1, Fig 5A). The presence of the olive wood in the studied context is remarkable, as there are only a few cases where its usage in archaeological human cremations have been reported [33, 37]. Nevertheless, macroremains of this taxa were already identified in other types of funerary archaeological contexts from the Middle Neolithic to Bronze period in Iberia [104107].

On the contrary, use of oak and pine woods is ubiquitous in different human cremations and burials of both prehistoric [24, 27, 29, 30] and historic sites [34, 35, 37, 38, 108, 109]. In all these studies the combustion properties of these two types of wood are mentioned. The simultaneous presence of oak and pine woods in a single funerary site has been observed since prehistory, suggesting a specific human choice [24]. Regarding the identification of F. cf. angustifolia, it must be emphasized that this species is numerically limited to a small number of charcoal fragments, and it is only present in the SU74 (Table 1, Fig 6B). The presence of ash wood within the Pit, despite being in small amount, is an interesting datum, because it is the only arboreal taxon of the anthracological record in addition to the three main taxa O. europaea, Quercus spp. (evergreen) and P. pinaster. There are plenty of recorded instances of use of ash wood for human archaeological cremations, although its use seems to become consistent rather than sporadic only during the historical time [34, 35, 37, 38, 108]. Regarding the presence of the evergreen taxa of Cistus sp., A. unedo and Fabaceae, numerous studies have already widely attested the use of shrubs and small trees in human cremation rituals. The Cistus sp. is a small perennial shrub that grows in dry and rocky soils of the Mediterranean area and its presence has been observed in other archaeological human cremations, especially during the Roman period [37, 108]. There are also records about the use of Fabaceae woods and of A. unedo tree wood in human cremations [38, 108]. A. unedo, whose presence Pit 16 at Perdigões is limited only to SU 74 (Table 1, Fig 6B), has also been used as a fuel in the Portuguese chalcolithic funerary site of Cabeço da Amoreira (Muge, Portugal) [24].

5.1.1. Different parts from trees.

Some anatomical alterations of charcoal from Pit 16, together with charcoal size classes were recorded. The most significant observations have been reported, although they should not be considered as conclusive, since they are not supported by trigonometric and statistical analysis.

Regarding charcoal size classes, results showed that larger charcoals belong to the three main arboreal taxa identified (Fig 7). This result seems predictable, considering that the rest of the anthracological record is mainly formed by shrubs, usually with thinner branches than a tree. Nevertheless, we cannot rule out the taphonomic processes undergone by charcoals, and the different degree of fragmentation proper to each species [110112].

For the charcoal fragments corresponding to the tree taxa identified at the anthracological record of Pit 16, i.e. O. europaea, Quercus spp. (evergreen), P. pinaster and F. cf. angustifolia (Fig 6A), a strong curvature of the ring was never. This data suggests that large branches or parts of the trunk, together with other minor parts of these trees were possibly employed for the cremation/s. In contrast, the observation of a strong ring curvature for all charcoal fragments belonging to A. unedo and Fabaceae seem to indicate the use of twigs or small branches of these taxa for the ritual/s. The same observation is attributable to all charcoal fragments of Cistus sp., for which a strong ring curvature was also observed.

Overall, the data suggest the possible use of different parts of the tree (larger and small branches, twigs, parts of the trunk), for the made up of the pyre/s formed by charcoal of the pit.

The radial cracks observed in some charcoals of P. pinaster could be the result of a conspicuous loss of humidity on burned fresh wood [113]; although some authors have argued that these features are not always directly related to green wood [66]. The absence of other signs of degradation, caused by fungi and/or insects [114] on charcoals under study, could be due to the use of P. pinaster wood shortly after the cutting of the tree.

The reaction wood feature observed on some charcoal fragments of P. pinaster, namely grooves inside tracheids visible on the longitudinal sections of charcoals, suggests a concentric growth in this part of the tree, developed during its life to prevent fall under its own weight [69].

5.2. Woods properties and their management

The characteristics and properties of the identified woods provide information regarding the selection, intentional or not, and use of wood by the Chalcolithic community which prepared the pyre/s deposited in Pit 16.

Let us first discuss the three main taxa. Nowadays the burning properties of olive wood are well appreciated, in particular due to its long-lasting fire and the pleasant aroma which derives from its burning. Furthermore, its sacredness is emphasized in the ethnographic literature [109]. It should also be emphasized that, as it is one of the main taxa of the anthracological record, this wood might have been intentionally selected by the community who made the ritual, and the selection could have been made for cultural and/or social reasons. As for the oak wood, the literature underlines how it also provides an excellent fuel, long lasting and with high calorific power [115117]. On the contrary, pine wood burns very quickly, reaching high temperatures rapidly, being usually used for fire lighting [115, 118]. The use of pine wood within funerary rituals for symbolic reasons cannot be excluded, especially given its high resin content, which would produce a pleasant smell during the incineration [37, 108].

Ash wood produces an excellent fuel, and its property of burning, both when is seasoned than when is fresh, is well known. The shrubs identified in the anthracological record of Pit 16, due to small-calibre woods, would seem particularly useful in lighting fires. Moreover, they emanate a pleasant smell during their combustion, and their use in these type of contexts may represent an anthropic symbolic choice [119]. The presence in funerary contexts of rock rose for symbolic reasons has also been contemplated in the literature [120].

Assemblages of concentrated charcoal, especially those coming from a single context, cannot be considered reliable for the reconstruction of the vegetation surrounding the archaeological site [25, 30, 35]. Despite that, the anthracological record of this type of assemblages reflects at least part of the vegetation of the site, representing a very important source of information. In the case of Pit 16 the identified taxa, if appropriately compared with the results of pollen analysis at Perdigões [121, 122] and with the current vegetation of the region, offer insights on the possible supply area of woods. These data provide indirect information on the possible area where the primary ritual could have taken place. Pollen analysis carried out on Late Neolithic layers of Perdigões identified the main three genera observed in the anthracological record of Pit 16, i.e., Olea, Pinus and Quercus. Inside the pollen records, the presence of various species belonging to Cistus and the presence of Fabaceae species has also been verified [121]. In the Chalcolithic layers of Perdigões, Quercus and Oleaceae, the olive family to which belong both Olea and Fraxinus genera, were identified [122]. Comparing the anthracological record with current vegetation data, we observe that all species present in Pit 16 are characteristic of the typically Mediterranean native vegetation, which surrounds the Perdigões archaeological complex [80, 123, 124]. According to these data, the origin of identified woods could have either been the surroundings of Perdigões or a related area of the same region. The discussion of the results of this study cannot, however, disregard the comparison with other archaeological data, which highlights how the action area of the communities of Perdigões was, actually greater than previously supposed [59]. These data are corroborated by some studies on the mobility of the communities [125, 126], which suggest that the presence of exogenous individuals in the archaeological site may have been large. Reflecting on these evidences, it is therefore necessary to report that wood species identified in the anthracological record of Pit 16 are typical of most of the Mediterranean basin [80, 127, 128].

Based on this data, it is challenging to formulate a single hypothesis concerning the origin of the wood taxa collected in the secondary cremation of Pit 16. At this stage it is not possible to exclude that the primary cremation occurred neither in the vicinity of the archaeological site nor in slightly more distant areas.

5.3. One context and two rituals?

Starting from the relative frequencies of each species identified in the anthracological record of Pit 16 and their corresponding mean value, obtained by applying the "ubiquity correction" method, results shed light on the original context of archaeological charcoals and on the ritual associated therewith. Differences in the number of charcoals and in the number of taxa within the two stratigraphic units (SU72 and SU74) (Table 1, Fig 6B), seem to suggest that woods were intentionally selected to build the pyre/s during the primary ritual of cremation/s. This consideration would strengthen the archaeological data. Moreover, these data allow to hypothesize that the two different stratigraphic units might be two different and subsequent depositional moments, likely as deriving from two different cremations. This hypothesis seems even more reasonable by looking at the differences in the relative frequencies and in the mean values of each identified taxon (Table 1, Fig 6A and 6B). Most significant differences in values are related to the species of O. europaea and of P. pinaster (Table 1, Fig 6). Regarding the rest of the taxa, F. cf. angustifolia and A. unedo are only present in SU74 (Table 1, Fig 6B).

These values confirm that the distribution of taxa through the two stratigraphic units is not homogeneous (Fig 6B).

5.4. Forecasted combustion temperature of charcoals from Pit 16

The clear relationship between FT-MIR data on the test samples and the combustion temperature, in the range tested here, has proved the potential of this technique to correctly recreate burnt temperatures for wood charcoal produced during the primary cremation ceremony.

All archaeological charcoals retrieved in SU74 showed a relative high burn temperature (> 500°C). This result was in line with data reported from other authors [53] who determined the absolute combustion temperature of archaeological charcoal samples from human cremations using reflectance. Nevertheless, in SU72 only the Olea charcoal sample particles reached the temperature of 500°C.

The burn temperature calculated by the combination of FT-MIR and chemometrics analysis, in this experiment, may suggest the existence of a single cremation for SU 74 (i.e., the human cremation was made using a mixture of different types of woods in a single pyre). On the other hand, the statistical analysis of the combustion temperature of charcoal sample particles found in SU74 determined the existence of a significant (p<0.05) difference among the temperature reached by the different types of wood. This difference may be due to different factors, such as: i) the chemical composition of each type of wood (inflammable compounds like resins or epicuticular waxes); ii) the wood moisture content [129] and/or iii) the resistance to burning (protective cover like suberin) [130]. A single deposition coming from a single pyre, employed for the primary cremation regarding the SU74 could be explained by the significant different combustion temperature of each taxon.

In the case of SU72, there were not significant differences between the burning temperature of Quercus and Pinus archaeological charcoal samples. The highest temperature registered by Olea may be due to the low wood moisture content in this kind of wood taxa. However, the burn temperature of Olea samples from SU72 was significantly different from the rest of the main taxa (Quercus and Pinus samples). In the case of Pinus samples, a significant temperature difference of 89°C (Table 2) was observed between the two stratigraphic units SUs 72 (486°C) and 74 (575°C). This fact may reflect the use of relatively green Pinus’ individuals in SU72, which usually displayed a relative low concentration of resins and waxes compounds (inflammable), burning at lower temperatures. The possible use of some green woods for the cremation would explain the low combustion temperature of Pinus samples from SU 72. In contrast, in SU 74, the highest temperature reached by Pinus may reflect the use of seasoned individuals with a relatively noticeable amount of resinous material (high degree of flammability). Therefore, in both cases, the wood seasoning of pine samples could be considered an important actor to determine the temperature of combustion.

From the combustion temperature differences of the tree main taxa of Pit 16, especially regarding charcoals fragments of P. pinaster, observed between the two stratigraphic units SUs 72 and 74, we can put forward the hypothesis of the existence of 2 different depositional moments, inside the Pit and consequently coming from two different primary rituals (cremations). This hypothesis was already raised with the anthracological analysis and it corroborates the archaeological data regarding the distinction of the two different stratigraphic units.

The absolute relative temperature reached by Quercus charcoals in SU72 is comparable to that one reached in SU74. This may be due to the use of Quercus suber wood in both crematory events. The existence of a burning protective cover, like suberin, on trunk would have been able to limit and “average” the combustion temperature of the Quercus wood in both crematory events. Moreover, Q. suber tree displays a relatively high abundance of this species in the Alentejo region [88]. Results showed in this study represent a first step, since it is not taken in account the position of the fuel (top of pyre, center of pyre, outside of pyre), as well as the existence of another type of fuel (e.g., animal fat). Some authors [131] demonstrated the significant role of these kind of factors in the combustion temperature. Therefore, the subsequent experiments will take to account the aforementioned variables.

5.5. Culturally selection vs availability selection?

The anthracological study of Pit 16 has proved to be particularly relevant for the understanding of human choices of timber for funerary purposes, demonstrating its importance in discerning the complexity of the human cremation process: do the culturally determined wood selection and the wood selection based on the territorial availability represent two different solutions, two coexistent solutions or are they intrinsic solutions?

If we consider the parameter known as the "principle of least effort" [132], we must take into consideration the possibility that this criterion was applied by the community of Perdigões during the collection of woods for cremation purposes. However, as some authors have already written [30], the taxa availability in the surroundings must be considered a priority for this type of communities. Furthermore, the identification of different types of woods can be interpreted as a decision of the community not to choose one or a few species for specific ritual reasons, but rather indicating a collection based on the easy availability of woods. On the other hand, it is necessary to consider that in the anthracological record of the Pit 16 three woods are dominant: the wild olive, the evergreen oak, and the maritime pine. We know that the first species gains a symbolic significance during historical times, while there is also numerous evidence of the use of the second and the third species in other prehistoric cremations. This is the first time that these kinds of evidence are studied in Perdigões and these results will be essential for the interpretation of the funerary practices in the archaeological site during the Chalcolithic. The data collected so far, and its comparison with other case studies [25, 133] from different prehistoric cremations and burials seems to indicate that cultural choices and choices dictated by the environment go hand in hand in all contexts investigated, emphasizing the intrinsic nature of the two variables.

Answering such a wide question by looking for a singular solution, when human choices were dictated by multiple factors in all periods of communal human existence, past and present, is an unproductive effort.

6. Conclusion

The charcoal record of Pit 16 is a mixture of different trees species. The most abundant tree taxa are O. europaea, Quercus spp. (evergreen) and P. pinaster. All the species identified in the study are characteristic of both deciduous and evergreen Mediterranean vegetation. The origin of plants selected seems to indicate the occurring of the human cremation/s either on site, in the vicinity or in an area with access to a similar vegetation.

Regarding the investigation of the charcoals burn temperature, this work has demonstrated the validity of the combination of FT-MIR and chemometric analyses to establish the absolute combustion temperature of archaeological charcoals from a human cremation. Considering both the data from chemometric and anthracological analyses, it was possible to speculate the existence of two different depositional moments inside Pit 16.

As for the understanding of human choices of wood species employed for funerary purposes, results from this study seem to indicate that both cultural choices and the ones dictated by the environment might have driven the Chalcolithic community of Pit 16 in the wood selection for the funerary ritual/s.

This study demonstrated the validity of the use of different approaches, an archaeobotanical, a chemical and a statistical approach, to shed light on the role of woods within the ritual of cremation within the funerary practices present in Perdigões during the Chalcolithic.

References

  1. 1. Evangelista LS, Valera AC. Segmenting and depositing: the manipulation of human body in ditched enclosures seen from Perdigões. NIA. In: Valera AC, editor. Fragmentation and Depositions in pre and proto-historic Portugal. NIA. Lisbon; 2019. pp. 47–70.
  2. 2. Valera AC, Silva AM, Cunha C, Evangelista LS. Funeray practices and body manipulation at Neolitich and Calcholithic Perdigões ditched enclosures (South Portugal). Recent Prehistoric Enclosures and Funerary Practices in Europe Proceedings of the International Meeting Edited by. 2014. pp. 37–58.
  3. 3. Valera AC. O Sepulcro 4 dos Perdigões. Um Tholos da segunda metade do 3° milénio AC. Perdigões Monográfica 2 Lisboa, NIA-ERA. 2020;2.
  4. 4. Valera AC, Godinho R. Ossos humanos provenientes dos fossos 3 e 4 e gestão da morte nos perdigões. Apontamentos Arqueol e Património. 2010;6: 29–40.
  5. 5. Godinho RM, Gonçalves D, Valera AC. The preburning condition of Chalcolithic cremated human remains from the Perdigões enclosures (Portugal). Int J Osteoarchaeol. 2019;29: 706–717.
  6. 6. Weiss-Krejci E. Formation processes of deposits with burned humanremains in Neolithic and Chalcolithic Portugal. ournal Iber Archaeol. 2005; 37–73.
  7. 7. Valera, AM. Bartelheim, P. Bueno Ramirez & MKC. The “exogenous” at Perdigões: Approaching interaction in the late 4th and 3rd Millennium BC in Southwest lberia. In: M. Bartelheim, P. Bueno Ramirez & MK, editor. Key resources and sociocultural developments in the Iberian chalcolithic. Tubingen: Tübingen Library Publishing; 2017. pp. 201–224.
  8. 8. Evangelista LS. Resting in Peace or in Pieces? Tomb I and Death Management in the 3rd Millennium BC at the Perdigões Enclosure (Reguengos de Monsaraz, Portugal). Resting Peace or Pieces? Tomb I Death Manag 3rd Millenn BC Perdigões Enclos (Reguengos Monsaraz, Port. 2019.
  9. 9. Valera AC. Mind the Gap: Neolithic and Chalcolithic Enclosures of South Portugal. BAR Intern. In: Gibson A, editor. Enclosing the Neolithic Recent studies in Britain and Europe. BAR Intern. Oxford: Archaeopress; 2012. pp. 165–183.
  10. 10. Valera AC. Ditches, Pits and Hypogea: new data and new problems in south Portugal late Neolitich and Chalcolitich practices. BAR Intern. In: Gibaja JF, Carvalho AF, Chambon F, editors. Funerary Practices in the Iberian Peninsula from the Mesolithic to the Chalcolithic Edited by BAR International Series 2417. BAR Intern. Oxford: Archaeopree, Publishers of British Archaeological Reports Gordon; 2012. pp. 103–112.
  11. 11. Weiss-Krejci E. Animals in mortuary contexts of Neolithic and Chalcolithic Iberia. In: Bicho N, editor. Animais na Pré-história e Arqueologia da Península Ibérica IV Congresso de Arqueologia Peninsular. Faro: Promontoria Monografica 03; 2006. pp. 35–45.
  12. 12. Duarte C, Evangelista L, Lago M, Valente MJ, Valera A. Animal remains in Chalcolithic funerary context in Portugal. The case of Pedrigōes (Reguengos de Monsaraz, Alentejo). Animais na Pré-história e Arqueol da Península Ibérica IV Congr Arqueol Penins. 2004; 47–55.
  13. 13. García Rivero D, J.O’Brien M. Among the dead. Late neolithic engraved slate plaques from the southweastern Iberian Peninsula. Estud Arqueol oeiras. 2014;21: 81–108.
  14. 14. Silva AM. Portuguese populations of late Neolithic and Chalcolithic periods exhumed from collective burials: An overview. Anthropol. 2003;41: 55–64.
  15. 15. Silva AM. Antropologia funerária e Paleobiologia das Populações Portuguesas (Litorais) do Neolítico Final/Calcolítico. Fundação Calouste Gulbenkian, editor. Lisboa; 2012.
  16. 16. Weiss-Krejci E. Changing Perspectives on Mortuary Practices in Late Neolithic / Copper Age and Early Bronze Age Iberia. Oxbow Book. In: Lillios K, editor. Comparative Archaeologies The American Southwest (AD 900–1600) and the Iberian Peninsula (3000–1500 BC). Oxbow Book. Oxford; 2011. pp. 153–174.
  17. 17. Boaventura R, Ferreira M, Neves MJ, Silva AM. Funerary practices and anthropology during middle-late Neolitich (4th and 3rd millennia BCE) in portugal: old bones, new insights. Anthropologie. 2014;LII: 183–205.
  18. 18. Cabaço N. Restos faunísticos em contexto funerário nos Perdigões, Reguengos de Monsaraz (Sepulcros 1 e 2). Arqueologica, Jia 2011 IV jornadas de jovens em Investigação. 2011. p. 267.
  19. 19. Cabrero García R, Ambrosiani Fernández J, Guijo Mauri JM, Gómez Murga E. Estudio de restos humanos procedentes del Dolmen de Cañada Real depositados en el Departamento de Anatomía y Embriología Humanas de la Facultad de Medicina de la Universidad de Sevilla. Spal Rev Prehist y Arqueol la Univ Sevilla. 2005;14: 59–74.
  20. 20. Cámara Serrano JA, Riquelme Cantal JA, Pérez Bareas C, Lizcano Prestel R, Burgos Juárez A, Torres Torres F. Sacrificio de animales y ritual en el Polideportivo de Martos-La Alberquilla (Martos, Jaén). Cuad Prehist y Arqueol la Univ Granada. 2012;20: 295–327.
  21. 21. Costa Caramé M, García Sanjuán L, Murillo-Barroso M, Parrilla Giráldez R, Wheatley D. Artefacts Produced in Rare Rocks from Funerary Contexts of the 4th-2nd Millenia Cal BCE in Southern Spain: A Review. Menga J Prehistory Andalusia. 2011;1: 253–293.
  22. 22. Cruz A, Berruti GLF, Delfino D, Gaspar F. The lithic hoard of hypogeum of colos (abrantes, middle tagus valley-Portugal): Shadows of rituals. Mediterr Archaeol Archaeom. 2016;16: 75–86.
  23. 23. Días-Zorita M, Eleazer Costa M, Garcia Sanjuan L. Funerary practices and demography from the mesolithic to the copper age in southern spain. BAR Intern. Funerary Practices in the Iberian Peninsula from the Mesolithic to the Chalcolithic. BAR Intern. Oxford: Archaeopress; 2012. pp. 51–65.
  24. 24. Diogo Monteiro PA, Cascalheira J, Marreiros JM, Pereira T, Ferreira Bicho N. Fire as a Component of Mesolithic Funerary Rituals: Charcoal analyses from a burial in Cabeço da Amoreira (Muge, Portugal). Wild things Recent Adv Palaeolithic Mesolith Res. 2015; 126–132.
  25. 25. Moskal-del Hoyo M. The use of wood in funerary pyres: Random gathering or special selection of species? Case study of three necropolises from Poland. J Archaeol Sci. 2012;39: 3386–3395.
  26. 26. Thomson GB. The analyses of wood charcoals from selected pits and funerry contexts. In: Barclay A, Halpin C, editors. The Neolithic and Bronze Age Monument Complex, Thames Valley Landscape. Oxford: Oxford Archaeological Unit; 1999. pp. 247–253.
  27. 27. Vidal Matutano P. Cueva Maturras (Ciudad Real, España): El papel del fuego en un contexto funerario del III milenio A.C. SAGVNTVM Papeles del Lab Arqueol Val. 2013;45: 39–47.
  28. 28. Gale R. Charcoal from the cemetery. In: Finn. N., editor. Bronze Age Ceremonial Enclosures and Cremation Cemetery at Eye Kettleby, Leicestershire: The Development of a Prehistoric Landscape. Leicester Archaeology Monograph; 2011.
  29. 29. O’Donnel L. People and woodlands: an investigation of charcoal remains as indicators of cultural selection and local environment in Bronze Age Ireland. In: BADAL, E., CARRIÓN, Y., GRAU, E., MACÍAS, M., NTINOU M, editor. 5th International Meeting of Charcoal Analysis: Charcoal as Cultural and Biological Heritage Sagvntvm extra 11. 2011. pp. 97–98.
  30. 30. O’Donnell L. The power of the pyre—A holistic study of cremation focusing on charcoal remains. J Archaeol Sci. 2016;65: 161–171.
  31. 31. Castiglioni E, De Carlo SM, Rottoli M. Il combustibile neue cremazioni dell’italia nordorientale. Bull la Soc Bot Fr Actual Bot. 1992;139: 311–318.
  32. 32. Cooremans B. The Roman cemeteries of Tienen and Tongeren: results from the archaeobotanical analysis of the cremation graves. Veg Hist Archaeobot. 2007;17: 3–13.
  33. 33. Coubray S. Combustibles, modes opératoires des bûchers et rituels. L ‘ analyse anthracologique. In: W. Van Andringa, H. Duday, S. Lepetz. " Mourir à Pompéi: fouille d ‘ un quartier funéraire de la n… Collection. In: Andringa V, Duday H, Lepetz S, editors. Mourir à Pompeí: fouille d’um quartier funéraire de la nécropole romaine de Porta Nocera (2003–2007). Collection. Rome; 2013. pp. 1433–1449.
  34. 34. Deforce K, Haneca K. Firewood selection in Roman cremation rituals in Northern Gaul:ritual, functional or opportunistic? Saguntum Papeles del Lab Arqueol Val. 2011;11: 193–194.
  35. 35. Deforce K, Haneca K. Ashes to ashes. Fuelwood selection in Roman cremation rituals in northern Gaul. J Archaeol Sci. 2012;39: 1338–1348.
  36. 36. Endrizzi L, Degasperi N, Castiglioni E, Rottoli M. La Sepoltura Romana rinvenuta in località Lanz (Cloz, Val di Non). AdA Archeol delle Alpi. 2014;1: 67–75.
  37. 37. Bel V, Figueiral I, Fabre L. Considerations on the nature and origin of wood-fuel from gallo-roman cremations, in the Languedoc region (southern France). Réflexions sur la nature et l’origine des bois utilisées dans la crémation gallo-romaine, en Languedoc (Sud de la France). Quaternaire. 2010;21.
  38. 38. Gutiérrez-Cuenca E, Hierro-Gárate JÁ, L. López-Dóriga I, Martín-Seijo M. Fires in the dark. Wood and charcoal analysis of the Early Medieval funerary deposits in the cave of Riocueva (Cantabria, Spain). Estud do Quaternário / Quat Stud. 2017; 73–85.
  39. 39. Aldritt D. The wood charcoal from Saltwood Tunnel (Unpublished CTRL apecialist report). Kent; 2006.
  40. 40. Kreuz A. Charcoal from ten early neolithic settlements in central europe and its interpretation in terms of woodland management and wildwood resources. Bull la Soc Bot Fr Actual Bot. 1992;139: 383–394.
  41. 41. Out , Welmoed A. Firewood collection strategies at Dutch wetland sites in the process of Neolithisation. The Holocene. 2009;20: 191–204.
  42. 42. Santín C, Doerr SH, Kane ES, Masiello CA, Ohlson M, de la Rosa JM, et al. Towards a global assessment of pyrogenic carbon from vegetation fires. Glob Chang Biol. 2016;22: 76–91. pmid:26010729
  43. 43. Santín C, Doerr SH. Fire effects on soils: The human dimension. Philos Trans R Soc B Biol Sci. 2016;371. pmid:27216528
  44. 44. De la Rosa JM, Jiménez-González MA, Jiménez-Morillo NT, Knicker H, Almendros G. Quantitative forecasting black (pyrogenic) carbon in soils by chemometric analysis of infrared spectra. J Environ Manage. 2019;251: 109567. pmid:31569023
  45. 45. Goldberg ED. Black Carbon in the Environment: Properties and Distribution. first edit. New York: Wiley; 1985.
  46. 46. Hedges JI, Eglinton G, Hatcher PG, Kirchman DL, Arnosti C, Derenne S, et al. The molecularly-uncharacterized component of nonliving organic matter in natural environments. Org Geochem. 2000;31: 945–958.
  47. 47. Asouti E. Charcoal analysis Web. Cecilia A. Wester Wood reference collection archive. http://pcwww.liv.ac.uk/~easouti/Cecilia%20A.%20Western%20Wood%20Reference%20Collection%20Notebook.html and http://pcwww.liv.ac.uk/~easouti/Cecilia%20A.%20Western.htm. 2009.
  48. 48. Braadbaart F, Wright PJ, van der Horst J, Boon JJ. A laboratory simulation of the carbonization of sunflower achenes and seeds. J Anal Appl Pyrolysis. 2007;78: 316–327.
  49. 49. Braadbaart F, Poole I. Morphological, chemical and physical changes during charcoalification of wood and its relevance to archaeological contexts. J Archaeol Sci. 2008;35: 2434–2445.
  50. 50. Antal MJ, Grønli M. The art, science, and technology of charcoal production. Ind Eng Chem Res. 2003;42: 1619–1640.
  51. 51. Hall G, Woodborne S, Scholes M. Stable carbon isotope ratios from archaeological charcoal as palaeoenvironmental indicators. Chem Geol. 2008;247: 384–400.
  52. 52. Kim NH, Hanna RB. Morphological characteristics of Quercus variabilis charcoal prepared at different temperatures. Wood Sci Technol. 2006;40: 392–401.
  53. 53. Veal R, O’Donnell L, McParland L. Reflectance—Current state of research and future directions for archaeological charcoal; Results from a pilot study on Irish Bronze Age cremation charcoals. J Archaeol Sci. 2016;75: 72–81.
  54. 54. Lebon M. Caractérisation par Spectroscopie Infrarouge à Transformée de Fourier des ossements chauffés en contexte archéologique—Comparaison entre référentiel moderne et matériel archéologique, Implication diagénétique. Muséum national d’Histoire naturelle. Paris. 2008.
  55. 55. Lebon M, Reiche I, Gallet X, Bellot-Gurlet L, Zazzo A. Rapid quantification of bone collagen content by ATR-FTIR spectroscopy. Radiocarbon. 2016;58: 131–145.
  56. 56. Pizzo B, Pecoraro E, Alves A, Macchioni N, Rodrigues JC. Quantitative evaluation by attenuated total reflectance infrared (ATR-FTIR) spectroscopy of the chemical composition of decayed wood preserved in waterlogged conditions. Talanta. 2015;131: 14–20. pmid:25281067
  57. 57. Fernández-Getino AP, Hernández Z, Piedra Buena A, Almendros G. Exploratory analysis of the structural variability of forest soil humic acids based on multivariate processing of infrared spectral data. Eur J Soil Sci. 2013;64: 66–79.
  58. 58. Jiménez-González MA, Álvarez AM, Carral P, Almendros G. Chemometric assessment of soil organic matter storage and quality from humic acid infrared spectra. Sci Total Environ. 2019;685: 1160–1168. pmid:31390706
  59. 59. Valera AC, Evangelista LS. Anthropomorphic Figurines at Perdigões Enclosure: Naturalism, Body Proportion and Canonical Posture as Forms of Ideological Language. Eur J Archaeol. 2014;17: 286–300.
  60. 60. Almeida NJ, Valera A. Faunal Remains Associated with Human Cremations: The Chalcolithic Pits 16 and 40 from the Perdigões Ditched Enclosures (Reguengos de Monsaraz, Portugal). Open Archaeol. 2022;8: 765–786.
  61. 61. Silva AM, Leandro I, Pereira D, Costa C, Valera AC. Collective secondary cremation in a pit grave: A unique funerary context in Portuguese Chalcolithic burial practices. Homo. 2015;66: 1–14. pmid:25500530
  62. 62. Crivellaro A, Schweingruber FH. Atlas of Wood, Bark and Pith Anatomy of Eastern Mediterranean Trees and Shrubs. Haupt P, editor. Berlin, Heidelberg: Springer Berlin Heidelberg; 2013.
  63. 63. Nardi Berti R. La struttura anatomica del legno ed il riconoscimento dei legnami italiani di più corrente impiego (II edizione). Ivalsa. Berti S, Fioravanti M, Macchioni N, editors. 2006.
  64. 64. Schweingruber EH. Anatomie europäischer Hölzer. Anatomy of European woods. Paul Haupt, Bern, Stuttgart; 1990.
  65. 65. Schweingruber FH. Microscopic Wood Anatomy. Swiss Federal Institute For Forest, Snow and Landscape Research; 1990.
  66. 66. Théry-Parisot I, Henry A. Seasoned or green? Radial cracks analysis as a method for identifying the use of green wood as fuel in archaeological charcoal. J Archaeol Sci. 2012;39: 381–388.
  67. 67. Vidal-Matutano P, Henry A, Théry-Parisot I. Dead wood gathering among Neanderthal groups: Charcoal evidence from Abric del Pastor and El Salt (Eastern Iberia). J Archaeol Sci. 2017;80: 109–121.
  68. 68. Schweingruber F, Börner A, Ernst Detlef S. Atlas of Woody Plant Stems: Evolution, Structure and Environmental Modifications. 2006.
  69. 69. Marguerie D, Hunot JY. Charcoal analysis and dendrology: data from archaeological sites in north-western France. J Archaeol Sci. 2007;34: 1417–1433.
  70. 70. Mas B, Allué E, Alonso ES, Vaquero M. From forest to settlement: Magdalenian hunter-gatherer interactions with the wood vegetation environment based on anthracology and intra-site spatial distribution. Archaeol Anthropol Sci. 2021;13.
  71. 71. Allué E, Mas B. The meaning of Pinus sylvestris-type charcoal taphonomic markers in Palaeolithic sites in NE Iberia. J Archaeol Sci Reports. 2020;30: 102231.
  72. 72. Vidal-Matutano P, Alberto-Barroso V, Marrero E, García JC, Pou S, de la Rosa MA. Vitrified wood charcoal and burnt bones from the pre-Hispanic site of Chasogo (Tenerife, Canary Islands, Spain). J Archaeol Sci Reports. 2019;28: 102005.
  73. 73. Henry A, Théry-Parisot I. From Evenk campfires to prehistoric hearths: Charcoal analysis as a tool for identifying the use of rotten wood as fuel. J Archaeol Sci. 2014;52: 321–336.
  74. 74. McParland LC, Collinson ME, Scott AC, Campbell G, Veal R. Is vitrification in charcoal a result of high temperature burning of wood? J Archaeol Sci. 2010;37: 2679–2687.
  75. 75. Humphries C., Press J., Sutton D. Árvores de Portugal e Europa. Fapas. 2005.
  76. 76. Rodríguez-Ariza MO, Moya EM. On the origin and domestication of Olea europaea L. (olive) in Andalucía, Spain, based on the biogeographical distribution of its finds. Veg Hist Archaeobot. 2005;14: 551–561.
  77. 77. Figueiral I. Charcoal analysis and the history of Pinus pinaster (cluster pine) in Portugal. Rev Palaeobot Palynol. 1995;89: 441–454.
  78. 78. Rubiales JM, García-Amorena I, García Álvarez S, Morla C. Anthracological evidence suggests naturalness of Pinus pinaster in inland southwestern Iberia. Plant Ecol. 2009;200: 155–160.
  79. 79. Carrió JS, Navarro C, Navarro J, Munuera M. The distribution of cluster pine (Pinuspinaster) in Spain as derived from palaeoecological data: relationships with phytosociological classification. The Holocene. 2000;10: 243–252.
  80. 80. Carrión JS, Leroy SAG. Iberian floras through time: Land of diversity and survival. Rev Palaeobot Palynol. 2010;162: 227–230.
  81. 81. Esteban LG, de Palacios P, García-Iruela A, García-Esteban L, de Vega DG, García-Fernández F. Comparative wood anatomy in pinaceae with reference to its systematic position. Forests. 2021;12.
  82. 82. Esteban LG, Palacios P, De Heinz I, Gasson P, Garc A, Garc F, et al. Softwood Anatomy: A Review genus. 2023.
  83. 83. Baas P, Blokhina N, Fujii T, Gasson PE, Grosser D, Heinz I, et al. IAWA List of microscopic features for softwood identification. IAWA J. 2004;25: 1–70.
  84. 84. Chabal L, Fabre L, Terral JF., Théry-Parisot I. La botanique. Errance (E. Collection Archéologiques. Errance (E. Paris: Revue archéologique du Centre de la France; 1999. pp. 1–66.
  85. 85. Queiroz PF, Mateus JE, Van Leeuwaarden W, Dise DP. Castro marim e o seu território imediato durante a antiguidade. paleo-etno-botânica. relatório final. Trab do Cipa. 2006;95.
  86. 86. Queiroz PF. Novos carvões dos fornos romanos do Martinhal (Sagres). Terra Scenica. 2010;17: 1–11. Available: https://www.researchgate.net/%0Apapers3://publication/uuid/DC62B856-F19A-45C1-99BE-D2E11B9F4C63
  87. 87. Moskal-del Hoyo M. The ubiquity correction as an alternative method for the quantification of charcoal. Saguntum Papeles del Lab Arqueol Val. 2011;11: 35–36.
  88. 88. Blanco E, González M, Tenorio M, Bombín R, Antón M, Fuster M, et al. Los bosques ibéricos. Una interpretación geobotánica. Barcelona: Editorial Planeta; 1997.
  89. 89. Swift LW, Elliott KJ, Ottmar RD, Vihnanek RE. Site preparation burning to improve southern Appalachian pine–hardwood stands: fire characteristics and soil erosion, moisture, and temperature. Can J For Res. 1993;23: 2242–2254.
  90. 90. Viscarra Rossel RA. ParLeS: Software for chemometric analysis of spectroscopic data. Chemom Intell Lab Syst. 2008;90: 72–83.
  91. 91. Jiménez-Morillo NT, Cabrita MJ, Barrocas Dias C, González-Vila JF, González-Pérez JA. Pyrolysis-compound-specific hydrogen isotope analysis (δ2H Py-CSIA) of Mediterranean olive oils. Food Control. 2020;110: 107023.
  92. 92. Jiménez-Morillo NT, Almendros G, Miller AZ, Hatcher PG, González-Pérez JA. Hydrophobicity of soils affected by fires: An assessment using molecular markers from ultra-high resolution mass spectrometry. Sci Total Environ. 2022;817. pmid:35016935
  93. 93. Sharma RK, Wooten JB, Baliga VL, Hajaligol MR. Characterization of chars from biomass-derived materials: Pectin chars. Fuel. 2001;80: 1825–1836.
  94. 94. Cai HY, Güell AJ, Chatzakis IN, Lim JY, Dugwell DR, Kandiyoti R. Combustion reactivity and morphological change in coal chars: Effect of pyrolysis temperature, heating rate and pressure. Fuel. 1996;75: 15–24.
  95. 95. Moreno-Castilla C, López-Ramón M V., Carrasco-Marín F. Changes in surface chemistry of activated carbons by wet oxidation. Carbon N Y. 2000;38: 1995–2001.
  96. 96. Fuente E, Menéndez JA, Díez MA, Suárez D, Montes-Morán MA. Infrared spectroscopy of carbon materials: A quantum chemical study of model compounds. J Phys Chem B. 2003;107: 6350–6359.
  97. 97. Oren A, Chefetz B. Sorptive and Desorptive Fractionation of Dissolved Organic Matter by Mineral Soil Matrices. J Environ Qual. 2012;41: 526–533. pmid:22370415
  98. 98. Michell A, Higgins H. Infrared Spectroscopy. Australian Forest Products research. Melbourn: CSIRO Forestry and Forest Products; 2002. p. 60.
  99. 99. Nakamoto K. IR Raman Spectra of Coordination Compounds. 2009.
  100. 100. Esteves B, Marques AV, Domingos I, Pereira H. Chemical changes of heat treated pine and eucalypt wood monitored by ftir. Maderas Cienc y Tecnol. 2013;15: 245–258.
  101. 101. Știrbescu NM, Ion RM, Teodorescu S, Știrbescu RM, Bucurică AI. Natural Resin Materials-Spectral and Chromatographic Investigations. Bull Transilv Univ Braşov •. 2017;10. Available: http://fsim.valahia.ro/STUCCO-MAT.html/docs/Stirbescu_NM.pdf
  102. 102. Martín-Ramos P, Fernández-Coppel IA, Ruíz-Potosme NM, Martín-Gil J. Potential of ATR-FTIR Spectroscopy for the Classification of Natural Resins. Biol Eng Med Sci Reports. 2018;4: 03–06.
  103. 103. Schultz J. Sclerophyllous vegetation. The ecozones of the world. Berlin: Springer; 2005. pp. 144–145.
  104. 104. Agustí B., Alcalde G., Burjachs F., Buxó R. J, Uan-Muns N., Oller J., et al. Dinàmica de la utilització de la Cova 120 per l’home en els darrers 6000 anys. Sèrie monogràfica 7. Girona: Centre d’Investigacions Arqueològiques.; 1987.
  105. 105. Hopf M, Catalán P. Neolithische Getreidefunde in der Höhle von Nerja (Prov. Málaga). Madrider Mitteilungen. 1970;11: 18–34.
  106. 106. Buxó IR. Seeds and Fruits. Gathering and Agriculture in the Mediterranean Area of France and Spain from Neolithic to Iron Age (en). Universitá de Montpellier 2. 1993.
  107. 107. Stika H. Botanische untersuchungen in der bronzezeitlichen höhnsiedlung Fuente Alamo. Madrieder Mitteilungen. Madrider Mitteilungen. 1988;29.
  108. 108. Martín-Seijo M, César Vila M. Oak, ash and pine: the role of firewood in funerary rituals at the Roman site of Reza Vella (Ourense, Spain). Archaeol Anthropol Sci. 2018;11: 1911–1926.
  109. 109. Hansson AM, Heiss A. Plants Used in Ritualn offerings and in Festive Contexts. In: Anderson Patricia, C., Pena-Chocarro Leonor, Andreas GH, editor. Plants and People Choise and Diversity trought Time. Oxford: Oxbow Books; 2014. pp. 610–663.
  110. 110. Chrzazvez J, Théry-Parisot I, Fiorucci G, Terral JF, Thibaut B. Impact of post-depositional processes on charcoal fragmentation and archaeobotanical implications: Experimental approach combining charcoal analysis and biomechanics. J Archaeol Sci. 2014;44: 30–42.
  111. 111. García Martínez MS, Dufraisse A. Correction factors on archaeological wood diameter estimation. Saguntum Papeles del Lab Arqueol Val Extra 13 Wood charcoal Evid Hum Nat Hist. 2012; 283–290.
  112. 112. Chabal L. Forêt et sociétés en Languedoc (Néolithique final, Antiquité final, Antiquité tardive). L’anthracologie, méthode et paléoécologie. Paris: Èditions de la Maison des Sciences de l’Homme; 1997.
  113. 113. Tereso JP. Vestígios arqueobotânicos do III milénio cal BC de Chibanes (Palmela, Setúbal). Setúbal Arqueol. 2014;15: 173–180.
  114. 114. Moskal-del Hoyo M, Wachowiak M, Blanchette RA. Preservation of fungi in archaeological charcoal. J Archaeol Sci. 2010;37: 2106–2116.
  115. 115. Gale R, Cluter D. Plants in Archaeology: Identification Manual of Vegetative Plant Materials Used in Europe and the Southern Mediterranean to c. 1500. KEW: Westbury Publications, Royal Botanic Garden; 2000.
  116. 116. Peña chocarro L, Peña Zapata L, González UJE, Ibáñez Estévez JJ. Agricultura, alimentación y uso del combustible: aplicación de modelos etnográficos en arqueobotánica botânica. Entorn—Estud Patrim Cult i Medi Ambi. 2000;3: 403–420. Available: http://www.entornsl.com/Congreso/pages/Actas/SeccionD.htm
  117. 117. Zapata Peña L, Peña Chocarro L. Uso y gestión del bosque en la Euskal Herria atlántica: aprovechamiento tradicional de los recursos forestales en Encartaciones y Gorbea. Zainak Cuad Antropol. 2003; 155–169.
  118. 118. Théry-Parisot I, Thiébault S., Delannoy JJ, Ferrier C, Feruglio V, Fritz C, et al. Illuminating the cave, drawing in black: wood charcoal analysis at Chauvet-Pont d’Arc. Antiquity. 2018;92: 320–333.
  119. 119. Segura Munguía S, Torres Ripa J. Historia de las Plantas en el mundo antiguo. CSIC PRESS, editor. 2009.
  120. 120. Figueiral I, Fabre L, Bel V. Considerations on the nature and origin of wood‑fuel from gallo‑roman cremations in the Languedoc region (Southern France). Réflexions sur la nature et l’origine des bois utilisés dans la crémation gallo-romaine, en Languedoc (Sud de la France). Quaternaire. 2010; 325–331.
  121. 121. Danielson R, Marques Mendes P. Pollen analyysis of late neolithic ditch deposits from the perdigões archaeological site. Apontamentos Arqueol e Património. 2013;9: 13–20.
  122. 122. Wheeler J. Paleoenvironmental assement of two archaeological sediment from Perdigões, Alentejo region, Portugal. Apontamentos Arqueol e Património. 2010;6: 41–46.
  123. 123. Costa JC, Neto C, Aguiar C, Capelo J, Espírito Santo MD, Honraldo J, et al. Vascular plant communities in Portugal (Continental, Azores and Madeira). Glob Geobot. 2012;2: 0–35.
  124. 124. Loidi J. The Vegetation of the Iberian Peninsula. Loidi J, editor. Cham: Springer International Publishing; 2017.
  125. 125. Žalaitė I, Maurer AF, Grimes V, Silva AM, Ribeiro S, Santos JF, et al. Diet and mobility of fauna from Late Neolithic–Chalcolithic site of Perdigões, Portugal. J Archaeol Sci Reports. 2018;19: 674–685.
  126. 126. Valera A. C; Žalaitė I.; Maurer A.F.; Grimes V.; Silva A.M.; Ribeiro S.; et al. Addressing human mobility in Iberian large prehistoric ditched enclosures: the case of Perdigões (South Portugal). J Archaeol Sci Reports. 2020;30.
  127. 127. Pignatti S. Evolutionary Trends in Mediterranean Flora and Vegetation. Plant Species and Plant Communities. Dordrecht: Springer Netherlands; 1978. pp. 157–167.
  128. 128. Thompson JD. Plant Evolution in the Mediterranean. Oxford University Press; 2005.
  129. 129. Dzurenda L, Banski A. Influence of moisture content of combusted wood on the thermal efficiency of a boiler. Arch Thermodyn. 2017;38: 63–74.
  130. 130. Jiménez-Morillo NT, de la Rosa JM, Waggoner D, Almendros G, González-Vila FJ, González-Pérez JA. Fire effects in the molecular structure of soil organic matter fractions under Quercus suber cover. Catena. 2016;145: 266–273.
  131. 131. Vanlandeghem M, Desachy B, Buonasera T, Norman L, Théry-Parisot I, Carré A, et al. Ancient arctic pyro-technologies: Experimental fires to document the impact of animal origin fuels on wood combustion. J Archaeol Sci Reports. 2020;33.
  132. 132. Kingsley Zipf G. Human Behavior and the Principle of Least Effort. Am Anthropol. 1950;52: 268–270.
  133. 133. O’Donnell L. People and woodlands: an investigation of charcoal remains as indicators of cultural selection and local environment in Bronze Age Ireland. Saguntum Papeles del Lab Arqueol Val. 2011;11: 97–98.