# Fire and brief human occupations in Iberia during MIS 4: Evidence from Abric del Pastor (Alcoy, Spain)

Carolina Mallol, Cristo Hernández, Norbert Mercier, Christophe Falguères, Ángel Carrancho, Dan Cabanes, Paloma Vidal-Matutano, Rory Connolly, Leopoldo Pérez, Alejandro Mayor, et al.

► To cite this version:

Carolina Mallol, Cristo Hernández, Norbert Mercier, Christophe Falguères, Ángel Carrancho, et al. Fire and brief human occupations in Iberia during MIS 4: Evidence from Abric del Pastor (Alcoy, Spain). Scientific Reports, 2019, 9 (1), 10.1038/s41598-019-54305-9. hal-02407221

HAL Id: hal-02407221

https://hal.science/hal-02407221

Submitted on 12 Dec 2019

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
OPEN

# Fire and brief human occupations in Iberia during MIS 4: Evidence from Abric del Pastor (Alcoy, Spain)

Carolina Mallol<fnref n="1" />,²*, Cristo Hernández<fnref n="1" />, Norbert Mercier D<fnref n="1" />,³, Christophe Falguères⁴, Ángel Carrancho⁵, Dan Cabanes⁶, Paloma Vidal-Matutano⁷,⁸, Rory Connolly<fnref n="1" />, Leopoldo Pérez⁹,¹⁰, Alejandro Mayor¹¹, Eslem Ben Arous D⁴ & Bertila Galván<fnref n="1" />

There is a relatively low amount of Middle Paleolithic sites in Europe dating to MIS 4. Of the few that exist, several of them lack evidence for anthropogenic fire, raising the question of how this period of global cooling may have affected the Neanderthal population. The Iberian Peninsula is a key area to explore this issue, as it has been considered as a glacial refugium during critical periods of the Neanderthal timeline and might therefore yield archaeological contexts in which we can explore possible changes in the behaviour and settlement patterns of Neanderthal groups during MIS 4. Here we report recent data from Abric del Pastor, a small rock shelter in Alcoy (Alicante, Spain) with a stratified deposit containing Middle Palaeolithic remains. We present absolute dates that frame the sequence within MIS 4 and multi-proxy geoarchaeological evidence of *in situ* anthropogenic fire, including microscopic evidence of *in situ* combustion residues and thermally altered sediment. We also present archaeostratigraphic evidence of recurrent, functionally diverse, brief human occupation of the rock shelter. Our results suggest that Neanderthals occupied the Central Mediterranean coast of the Iberian Peninsula during MIS 4, that these Neanderthals were not undergoing climatic stress and they were habitual fire users.

Archaeological evidence across the European continent from the end of MIS 5 through MIS 4 suggests that the cooling effects of the latter in the northern hemisphere influenced the Neanderthal population. There is an apparent reduction in the number of sites suggesting population shrinkage or redistribution<fnref n="1" />. At a regional scale, changes in behaviour have been proposed for Southwestern France, where Neanderthal occupations without fire during MIS 4 have been documented at different sites, motivating a hypothesis on the inability of these Neanderthals to make fire²⁻⁴. This hypothesis highlights the important role of climate in shaping human behaviour.

To better understand Neanderthal population dynamics and assess the degree of regional behavioural variability during critical conditions of environmental change it is important to seek evidence from the Iberian Peninsula, which has previously been considered as a refugium during cold periods, particularly the south of the peninsula⁵⁻⁸ and more recently the east coast as well⁹. Considerable work has been done to establish site-specific palaeoclimatic records for Iberian sites dating to MIS 5⁹⁻¹³ and MIS 3, e.g.¹⁴⁻¹⁹. In contrast, relatively few site-specific records exist for MIS 4 in the peninsula. This potentially hampers our understanding of the local environmental settings in which MIS 4 Neanderthal occupations occurred, as we depend on regional proxy data with coarse resolution.

<footnote n="1">¹UDI de Prehistoria, Arqueología e Historia Antigua, Departamento de Geografía e Historia, Facultad de Humanidades, Universidad de La Laguna, La Laguna, Spain. ²Archaeological Micromorphology and Biomarker Research Lab, University of La Laguna, La Laguna, Spain. ³Institute of Archaeomaterials Research, Université Bordeaux Montaigne, Pessac, France. ⁴UMR 7194, Département Homme et Environnement, Muséum National d'Histoire Naturelle, Paris, France. ⁵Área de Prehistoria, Departamento de Historia, Geografía y Comunicación, Universidad de Burgos, Burgos, Spain. ⁶Department of Anthropology, Rutgers University, New Brunswick, USA. ⁷Departamento de Ciencias Históricas, Universidad de Las Palmas de Gran Canaria, Las Palmas, Spain. ⁸Université Côte-d'Azur, CEPAM, CNRS, Nice, France. ⁹Institut Català de Paleoecologia Humana i Evolució Social, Tarragona, Spain. ¹⁰Área de Prehistoria, Universitat Rovira i Virgili, Tarragona, Spain. ¹¹Departament de Prehistòria, Arqueologia, Història Antiga, Filologia Grega i Filologia Llatina; Facultat de Filosofia i Lletres, Universitat d'Alacant, Sant Vicent del Raspeig, Spain. *email: cmallol@ull.es</footnote>
The Iberian Peninsula is characterised by considerable climatic diversity in which mountainous topography and Atlantic and Mediterranean effects cause local environmental conditions to vary significantly even across relatively short distances²⁰. It should be borne in mind, therefore, that local conditions may not always reflect wider regional trends. Marine records (e.g. MD95-2042) suggest that during MIS 4 the Iberian Peninsula saw the progressive development of semi-desert vegetation in line with a significant decrease in sea-surface temperature (SST)²¹,²². Palynological records from Abric Romani in the northeast indicate millennial-scale climatic variability through the course of MIS 4 into MIS 3, characterised by fluctuations between cold and warm phases driving the decline and expansion of arboreal pollen, including deciduous and evergreen *Quercus*, *Olea* and more²³,²⁴. This corresponds with terrestrial pollen data from Galicia in the northwest of the peninsula, where records for MIS 4 suggest that coniferous taxa are locally rare at this time, and the presence of *Erica*, *Calluna* and *Poaceae* suggests an environment comprised of alternating heath and temperate grasslands²⁵. Meanwhile, evidence from the Tagus Basin in central Iberia interior suggests that episodes of significant loess accumulation during MIS 4 are associated with the development of cold and dry conditions²⁶. On the other hand, palynological data from Bolomor Cave in Valencia points to the existence of a glacial refugium during MIS 6, with closed forests and high plant diversity⁹. This evidence could predict similar conditions for MIS 4 in the region and hints at a mosaic scenario for the Iberian Peninsula with different econiches shaped by local geographic factors.

Turning to the archaeological record, more Middle Palaeolithic sites in Iberia are dated to MIS 3 than to MIS 5 and MIS 4 (Supplementary Table S1). As a result, there is considerably more bioanthropological, genetic and cultural information on late Neanderthals (younger than 60 Ka). Although there is an apparent reduction of sites when compared with the preceding and subsequent periods, the significantly shorter duration of the MIS 4 period compared to MIS 5 and MIS 3 and the inaccuracy of some of the chronostratigraphic data should be considered.

The chronostratigraphic and archaeological data from Iberian MIS 4 sites is currently insufficient to understand Neanderthal population dynamics and behaviour in relation to the broader European context or to compare different proxies across MIS 5, MIS 4 and MIS 3 at a regional scale. Regarding fire, there is purported evidence of abundant charcoal, burnt sediment and burnt rocks in Navalmaçoil Layer F (central Spain), but this layer could date to late MIS 5 rather than MIS 4²⁷. The lithics from various Iberian MIS 4 sites reflect settlement patterns involving small territories, as informed by short-distance raw material procurement²⁷⁻³¹ and recurrent occupation of sites³².

The central Mediterranean region of Spain is prominent for its Neanderthal archaeological record. There are several sites in the region featuring deep stratigraphic sequences that span either the entire MIS 5 period, including Cova Negra³³ and Bolomor Cave³⁴, or the entire MIS 3 period, such as El Salt³⁵. Based on archaeological evidence from these sites, it appears that Neanderthals in this region had similar adaptations during both time intervals. Faunal assemblages are mostly composed of anthropogenic large and medium-sized ungulates, mainly deer, horse and wild goat, with minor proportions of small-sized animals (tortoise and rabbit) and reflect generalist, broad spectrum hunting traditions resulting from adaptation to biodiverse environments. Overall, the lithic technology from these MIS 5 and MIS 3 sites does not show any significant difference either³⁶⁻³⁸. Broadly, it is characterised by a predominance of recurrent centripetal Levallois (RCL) sequences and discoidal schemes, alongside a lesser representation of other Levallois modalities and non-Levallois technical strategies, such as polychedral and orthogonal procedures³⁹⁻⁴². This variability, documented throughout MIS 5 and MIS 3 contexts in the region, is possibly related to territorial aspects entailing differential raw material availability and site function determining the use of different activity-specific tools at different sites⁴³⁻⁴⁵.

Finally, current evidence shows that both MIS 5 and MIS 3 Neanderthals in the central Mediterranean region of Iberia made simple open hearths³⁵,⁴⁶,⁴⁷ and wood-gathering patterns appear to be similar irrespective of glacial, interglacial,stadial or interstadial conditions¹⁸,⁴⁸,⁴⁹.

Focusing on different aspects of behaviour possibly sensitive to changes in climate such as patterns of group mobility, site occupation duration and seasonality, data from the MIS 5 contexts in this region generally suggest that Neanderthals in the region were occupying relatively small territories, as shown by local lithic raw material procurement patterns³⁹,⁵⁰ and were recurrently occupying sites for short time periods as suggested by tool recycling and diachronic successions of discrete archaeological assemblages in Bolomor Cave³⁹,⁵¹. In Cova Negra, presence of carnivores throughout the sequence has been interpreted as evidence for short-term human occupations⁵² although no geoarchaeological or archaeostratigraphic data is available to corroborate this. For MIS 3, a lithic raw material study for El Salt shows the continued use of a river-bound territory of less than 5 km throughout the entire sequence⁵³,⁵⁴ and an archaeostratigraphic study of the faunal and lithic record from a segment of the sequence dating to early MIS 3 suggests a pattern of short-term, recurrent occupation of the site⁵⁵. Leierer et al.⁵⁶ corroborated this pattern using sedimentary evidence from the fire record and proposed relatively long periods of site abandonment, which could suggest high group mobility. Although these MIS 5 and MIS 3 data are not comparable directly due to their different degrees of temporal resolution, they generally reflect a Neanderthal population consisting of highly mobile groups occupying relatively small territories.

Considering the above-mentioned evidence, filling the gap of information for MIS 4 becomes relevant to: (1) establish the extent of environmental change in the region during MIS 4, and (2) assess the degree of Neanderthal behavioural change (including fire use) from MIS 5 to MIS 3 in a regional context. To this end, here we report archaeological evidence of human occupations with anthropogenic fire and absolute dates from Abric del Pastor, a Middle Palaeolithic rock shelter site in Alcoy, Spain (Fig. 1, Supplementary Fig. S1). The site is undergoing systematic investigation since 2006 and has yielded lithic, faunal and combustion assemblages whose study has provided valuable information regarding Neanderthal settlement dynamics⁵⁷,⁵⁸. The archaeological sequence appears to represent a succession of brief Middle Palaeolithic occupation events linked to different subsistence activities and these are all centered around single hearths, which are present throughout the sequence⁵⁷. Until now, ascribing the site to a MIS 4 chronological framework has been problematic and the *in situ* nature of the combustion features remained uncertain. A terminus *ante quem* for the sequence is given by the presence of bone remains of
<figure><img src="image_2.png" /><figcaption>**Figure 1.** Location of Abric del Pastor rock shelter.</figcaption></figure>

*Testudo Hermanni*, which diminishes in the regional record by early MIS 3⁵⁹,⁶⁰. The good states of preservation of the sedimentary deposit, together with weak frost-induced lenticular microstructures observed microscopically⁵⁷, suggest a generally cool or cold climate but is insufficient evidence to ascribe the sequence to MIS 4. A recent multiproxy study investigating the carbon and hydrogen isotopic values of sedimentary leaf waxes through the sequence (Units I-VI) point to a variable precipitation regime, while cryoturbation features revealed through soil micromorphology indicate generally cool conditions throughout⁵⁶,⁶¹. Based on the anthracological record, the Middle Palaeolithic human occupations are associated with dry/semi-arid supramediterranean conditions (mean annual temperature 8–13 °C and mean annual precipitation 200–600 mm) with juniper stands, heliophilous taxa and scattered cryophilous pine stands⁴⁸,⁴⁹.

Our ongoing investigations at Abric del Pastor are geared at understanding site formation processes and com-
bustion feature formation processes, dissecting archaeological palimpsests and providing a chronostratigraphic
context for the Middle Palaeolithic occupations. In this paper, our objectives are threefold: (1) to provide chrono-
metric dates that help establish if the archaeological deposit formed during MIS 5, MIS 4 or MIS3, and (2) to char-
acterize the combustion features and corroborate their nature as anthropogenic combustion structures, and (3)
to describe and assess the variability of the lithic, faunal and combustion record within its stratigraphic context.

**Results**

**Absolute dates.** We report three absolute dates for Stratigraphic Unit IV: (1) an ESR/U-series date (tooth sample) of 48 ± 5 for Stratigraphic Subunit IVb, (2) an OSL date (sediment sample) of 63 ± 5 for Stratigraphic Subunit IVd, and (3) an ESR/U-series date (tooth sample) of 62 + 10/−9 ka for Stratigraphic Subunit VI (Fig. 2; Table 1).

**Evidence of anthropogenic fire.** Abric del Pastor Stratigraphic Unit IV has yielded numerous burnt bone and flint fragments (Supplementary Table S2) and several sedimentary features with the appearance of combustion structures (Fig. 3). So far, 18 such features have been documented, concentrated in stratigraphic units II (1), IV (15) and VI (2) (see Fig. 2). In the field, these did not exhibit clear ash layers except for H17, which had a 1-2 cm-thick light gray, subcircular layer of silty sediment containing frequent burnt bone fragments. The others were 0.3-1 cm-thick circular or subcircular lenses of dark gray or reddish-brown sediment, with diameters between 0.5 and 1 m.

Micromorphological analysis of several of the combustion features showed the presence of microscopic com-
bustion residues including charcoal, wood ash and burnt bone (Supplementary Table S3 and Fig. S1). The gray
layer of H17 shows an aggregated, calcite-rich matrix with frequent unidentified black particles and coatings
around clasts, some of which are fissured and rubified (Fig. 4).

The calcitic fraction of the combustion features shows good preservation states. Fourier Transform Infrared Spectroscopy (FTIR) of loose sediment samples from selected combustion features shows calcite as a major mineral component, with minor amounts of clay and quartz (Supplementary Table S4). Highly disordered calcite attributable to wood ash has been identified in several samples using FTIR (Supplementary Fig. S2). The clay content in the site's sediments is too low to identify clay minerals heated above 500C unambiguously.

Magnetic susceptibility data for the area in and around feature H17 from Stratigraphic Unit IVf indicates burning (Fig. 5a and Supplementary Fig. S3). The values obtained for Unit IVf are low, as expected for a predominantly diamagnetic calcareous sedimentary deposit. However, the center of the feature and some of its adjacent zones yielded high values. The inner part of the combustion feature area and its surroundings is one order of magnitude higher. Anthracological analyses of the Unit IVf charcoal assemblage yielded a well-delimited cluster of strongly fragmented, predominantly coniferous charcoal at the center of the H17 combustion feature. The spatial distribution of this charcoal in conjunction with a zone of high magnetic susceptibility suggests presence of a hearth (Fig. 5b; Supplementary Tables S5 and S6).
<figure><img src="image_3.png" /><figcaption>**Figure 2.** Stratigraphic log of abric del Pastor showing the different lithostratigraphic units, the position of absolute dates and combustion structures.</figcaption></figure>

Lipid analysis of loose sediment samples collected from selected combustion features yielded *n*-alkane distributions with a predominance of long chain length homologues (nC27 – nC35), consistent with epicuticular leaf waxes derived from terrestrial higher plants⁶² (Fig. 6, Supplementary Table S7). Several samples yielded distributions with no evident odd-over-even predominance (OEP), which is likely associated with the thermal breakdown of long-chain homologues during combustion⁶³,⁶⁴. Other more polar lipid compounds were identified in the samples (Supplementary Table S8) but are not diagnostic of any specific biological source. It is not possible to establish if these compounds are synchronous with or related to the combustion event, or whether they entered the sediment postdepositionally. Compound specific carbon isotope analyses of long-chain *n*-alkanes nC29 and nC31 ranged between –34.3‰ to –31.7‰ and –34.3‰ to 32.8‰ respectively. These are values indicative of a vegetation source using the C3 photosynthetic pathway⁶⁵. Interestingly, when compared with control samples collected from outside the combustion features, hearth samples appear offset by ~2‰, which has been recorded elsewhere for *n*-alkane carbon signatures after exposure to high temperatures⁶⁴,⁶⁶.

**Hearth-related lithic and faunal assemblages from S.U. IVd and IVf.** In the field, the lithic, faunal and combustion remains from Unit IVd appeared to be spatially related (Fig. 7a). Zooarchaeological, technological and spatial analysis of the lithic and faunal record suggests the existence of hearth-related assemblages. The lithic assemblage is dominated by centripetal and preferential Levallois objects, with a minor representation of preferential and unidirectional strategies, discoid, and non-Levallois objects (Supplementary Table S9). Skeletal representation of the IVd faunal record is dominated by cranial remains, mainly scattered bovine and goat teeth. There are also frequent deer long bone fragments and abundant tortoise backplate and plastron plates (Supplementary Table S10). As indicated by technological and zooarchaeological features, different activities were taking place around the hearths. These analyses also show that the hearth-related assemblages were separated in time. First, a lack of bidirectional spatial relationships among different hearth-related lithic objects suggest that
<table><caption>**Table 1.** Absolute dates and associated data for teeth samples AP1601 and AP1602 (ESR/U-Series) and sediment sample OSL1 (OSL) calculated at 1-sigma confidence level.</caption><thead><tr><th rowspan="2">Sample</th><th rowspan="2">Type</th><th rowspan="2">Layer</th><th colspan="3">U Th K</th><th rowspan="2">a-efficiency</th><th rowspan="2">Internal</th><th colspan="2">external</th><th colspan="2">De</th><th colspan="2">Age</th><th colspan="4">p-value</th></tr><tr><th>(ppm)</th><th>(ppm)</th><th>(%)</th><th>+</th><th>-</th><th>dose rate (µGy/a)</th><th>+</th><th>-</th><th>(Gy)</th><th>+</th><th>-</th><th>(ka)</th><th>+</th><th>-</th><th>enamel</th><th>±</th><th>dentine</th><th>±</th><th>cement</th><th>±</th></tr></thead><tbody><tr><td>AP1601</td><td>Tooth</td><td>IV b</td><td>1,31</td><td>0,70</td><td>0,13</td><td></td><td>326</td><td>36</td><td>356</td><td>23</td><td>682</td><td>51</td><td>33</td><td>1</td><td>48</td><td>5</td><td>-0,69</td><td>0,12</td><td>-0,76</td><td>0,12</td><td>-49</td><td>0,14</td></tr><tr><td>OSL1</td><td>Sediment</td><td>IV d</td><td>1,31</td><td>0,70</td><td>0,13</td><td>5</td><td>0</td><td>0</td><td>681</td><td>21</td><td>681</td><td>21</td><td>43</td><td>1</td><td>63</td><td>5</td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>AP1602</td><td>Tooth</td><td>VI</td><td>1,31</td><td>0,70</td><td>0,13</td><td></td><td>749</td><td>155</td><td>378</td><td>23</td><td>1127</td><td>164</td><td>70</td><td>4,1</td><td>62</td><td>12/-9</td><td>-0,7</td><td>0,19</td><td>-89</td><td>0,15</td><td>-</td><td>-</td></tr></tbody></table>

<figure><img src="image_4.png" /><figcaption>**Figure 3.** Field views of two Abric del Pastor combustion structures: H2 in Unit IVb (left) and H17 in Unit IVf (right).</figcaption></figure>

they are not synchronous (Fig. 7b). Second, given the small surface area of the rock shelter, the high MNI for Unit IVd (nearly 16 individuals of different taxa) points to several human occupation events.

Unit IVf is an archaeologically poor deposit comprising a single hearth-related assemblage with frequent, strongly fragmented and variably burnt and unburnt faunal remains (n = 250), and very few (n = 6) small (<0.5 cm) flint flakes (Fig. 8). Except for an accumulation of faunal remains in the Southeast corner of the rock shelter, the bulk of the archaeological remains are concentrated around the H17 combustion feature. The low number of lithic objects and their markedly small size suggests a brief human occupation focusing in tool use rather than knapping activities, which has not been documented in the Iberian or European Middle Palaeolithic record.

Discussion

Our multidisciplinary study adds a new context to the relatively small set of Middle Palaeolithic sites dated to MIS 4. Considering the error ranges, the absolute dates obtained frame the Abric del Pastor sequence between 43 ka and 72 ka, covering the whole of MIS 4 and early MIS 3. Although the current dates are stratigraphically coherent, more dates need to be obtained to corroborate this chronological framework.

Our proxies of anthropogenic combustion in Stratigraphic Units IV and VI yielded positive results. The sampled sediment contains microscopic combustion residues indicative of burning which, based on their frequency and arrangement were slightly syn- and post-depositionally reworked. The sedimentary structure of Unit IVd is very open and loose and composed almost entirely of limestone gravel, sand and calcitic silt⁵⁷. Thus, individual small-sized particles on this type of substrate can be easily reworked by trampling. In H8 and H9 (Unit IVd), the presence of granular microstructures with sedimentary aggregates containing highly fragmented, silt-sized charcoal and ash and a mix of burnt and unburnt microscopic bone fragments suggest some reworking of the original combustion features. We did not observe any microstratigraphic features indicative of runoff or of any other water-related syn-sedimentary mechanism. Considering this evidence, the prominent absence of ash layers in the Unit IV combustion features - except for H17 in IVf -, could be a result of wind erosion in combination with trampling. This agrees with FTIR data showing random presence of wood ash in the samples, although these data should be taken with caution given the presence of microscopic secondary calcite throughout the Unit IV deposit, as observed in thin section. For now, we have not been able to isolate the signal of this (geogenic) secondary calcite and cannot be certain that it differs from our pyrogenic ash reference samples. Micromorphological samples from combustion features in Unit VI show more abundant combustion residues, including wood ash. These combustion structures, which predate a major roof collapse at the base of Unit IV, are better preserved.

Sediment samples from hearth H17 in Stratigraphic Unit IVf, which appeared in the field as a sub-circular gray zone with a few scattered bone fragments, contain microscopic char. Furthermore, the associated burnt bone remains, charcoal concentrations and magnetic susceptibility data point to *in situ* burning. High magnetic susceptibility values in the adjacent area could be explained by short-distance ash dispersal by wind, trampling or other syn-depositional agents. Charcoal analysis from H17 suggests the use of mainly juniper wood and the
<figure><img src="image_5.png" /><figcaption>**Figure 4.** Flatbed scan of a micromorphological thin section from combustion structure H17 in Unit IVf (left image) and selected photomicrographs taken from its black sediment: (a) Sedimentary matrix with abundant unidentified charred particles (black) in an ashy matrix (gray). Image taken in plane polarized light; (b) Same view in crossed polarized light; c) ashy sediment aggregate with frequent unidentified charred particles. Image taken in plane polarized light; d) Angular burnt bone fragment in plane polarized light.</figcaption></figure>

<figure><img src="image_6.png" /><figcaption>**Figure 5.** (a) Magnetic susceptibility contour map showing a high value zone at the H17 location (the H17 perimeter as recorded in the field is indicated by a dotted line); (b) Charcoal distribution contour map showing that most of the charcoal remains in Unit IVf cluster at the H17 location.</figcaption></figure>

concentration of fuel remains at the center of the hearth. Further spatial investigation is needed to assess postdepositional effects on charcoal dispersal.

All the combustion structures presented here (Units IVd, IVf and VI) date to MIS 4. Our integrated analysis of the Unit IVd-f hearth-related assemblages suggests that these represent a succession of brief, functionally diverse human occupation episodes. The Unit IVf assemblage is peculiar in that it could represent a single human occupation without knapping activity. The few lithic objects found (SI-1 Table 1; Fig. 8) are very small in size, predominantly burnt and could be from tool use rather than knapping debris (e.g. 67,68). Future use-wear analysis could clarify this issue. Our results agree with previous data from stratigraphically higher layers within Stratigraphic Unit IV, particularly Unit IVb, - dated here to early MIS 3-, which was interpreted as a succession of at least three brief human occupation episodes57. They also agree with the short-term settlement patterns observed for MIS 3 at El Salt55,56,69. El Salt is within 5 km distance from Abric del Pastor and both sites show the same lithic raw material types, all from nearby sources, pointing to exploitation of a single river-bound territory. Together, the current evidence from both sites suggests similar settlement patterns for the Neanderthal groups that occupied the this region of the Iberian Peninsula throughout MIS 4 and 3.

Regarding the palaeoclimatic context, our results also suggest that MIS 4 at Abric del Pastor was not too different from early MIS 3. Recent high-resolution geoarchaeological data point to a variable precipitation regime
<figure><img src="image_7.png" /><figcaption>**Figure 6.** Histograms displaying n-alkane concentrations in combustion features and associated control samples collected from thermally unaltered sediments immediately adjacent (x-axis = number of carbons; y-axis = concentration in µg per g of dry sediment).</figcaption></figure>

and generally cool conditions⁶¹. Previous anthracological data from stratigraphically higher layers within Unit IV (Units IVc)⁴⁸, show a rich diversity of woody taxa, pointing to supramediterranean conditions and readily available woody fuel at the site during the MIS 4 – MIS 3 transition. Instead, MIS 4 Unit IVf shows a reduction in diversity with *Juniperus* sp. as the dominant taxon. However, the unit IVf charcoal sample is associated with a single combustion structure and might reflect anthropogenic selection rather than the natural vegetation. Either way, wood was always an available resource. Further paleoenvironmental evidence from this and other mid latitude
<figure><img src="image_8.png" /><figcaption>**Figure 7.** (a) Spatial distribution of lithic and faunal remains and combustion structures from Unit IVd1; (b) Detailed spatial distribution of the refitted specimens within IVd1 related to combustion structure H9 (left) and a photograph of the refitting sets showing how they conform to an almost whole flint nodule (right).</figcaption></figure>

<figure><img src="image_9.png" /><figcaption>**Figure 8.** Spatial distribution of lithic and faunal remains associated with H17 within the context of Unit IVf.</figcaption></figure>

MIS 4 Middle Palaeolithic contexts will add to these data and help shape regional paleoenvironmental maps that will allow us to investigate Neanderthal adaptations in more detail.

Here, we have documented the presence of anthropogenic Middle Palaeolithic fire in MIS 4 in the central Mediterranean region of Iberia, which supports the prediction that southern European Neanderthals might have used fire more frequently than in higher latitudes⁷⁰. The associated local climatic context does not appear to have drastically changed during this period and Neanderthal groups living in this region were making fire and show similar settlement patterns to those from the following MIS 3 period. Comparable high-resolution geoarchaeological data is needed for MIS 5 at a regional scale, as well as a geoarchaeological focus on the MIS 5/MIS 4 and MIS4/MIS 3 stratigraphic boundaries at different sites, which might conceal valuable paleoclimatic information and contribute to our understanding of associated human dynamics.

## Methods and Materials

For this study we used two absolute dating techniques: Optically Stimulated Luminescence (OSL) on a loose sediment sample from Stratigraphic Unit IVd and ESR/U-Series on two teeth samples from stratigraphic units IVc and VI (See Supplementary Information for methodological details).

For our multiproxy investigation of the combustion features we used different techniques, which we applied to a selection of features according to their availability for sampling. Most of the features concentrate in Stratigraphic subunit IVd, a rich archaeological horizon. The techniques used were: (1) Soil micromorphology of intact sediment blocks to investigate microscopic components and their arrangement, (2) Fourier Transform Infrared Spectroscopy (FTIR) to identify calcareous wood ash and burnt clay, (3) Magnetic susceptibility to corroborate *in situ* burning, (4) Anthracology to characterize wood fuel and 4) Lipid analysis to seek combustion biomarkers (see Supplementary Information for methodological information on the different techniques).

Integrated analysis of the lithic and faunal record consisted in: (1) Technological analysis following standard analytical parameters, raw material units⁷¹ and refitting, (2) Zooarchaeological analysis using established standard
methods<sup>72,73</sup>, and (3) Spatial analysis of the lithic and faunal record using ArcGIS software (see Supplementary Information for further details).

**Site background and stratigraphic context.** Excavations and multidisciplinary investigations at Abric del Pastor have been carried out since 2006 by a research team directed by B.G., C.H. and C.M. The site comprises a small (40 m²) rock shelter containing a stratified clastic deposit dominated by coarse components (cobbles and gravel in a sandy matrix).

The Abric del Pastor stratified deposit is 1.3–1.5m-thick and has been subdivided into six lithostratigraphic units, I-VI⁵⁷, of which units II, IV, V and VI have yielded Middle Palaeolithic archaeological remains. The top of the sequence (unit I) is a recent Holocene deposit mainly composed of bioturbated, organic-rich sand with sheep dung. Units II-VI are Pleistocene calcitic deposits hosting the remains of several Middle Palaeolithic human occupation episodes. The bulk of the Middle Palaeolithic evidence retrieved and studied so far comes from Stratigraphic Unit IV, which is in turn comprised of several subunits: IVa-IVg (see Fig. 2). Of these, subunits IVa and IVc are very loose, clast-supported and archaeologically poor, while IVb and IVd are finer-grained and concentrate the bulk of archaeological remains. Excavation of the Unit IV sequence comprised a 38 m² surface area, a large portion of the total space within the small rock shelter. Units V and VI are only known from a 2 × 2 m test pit within the excavation area.

**Data availability**

The key data generated and/or analysed during this study are included in this published article (and its Supplementary Information files). The complete datasets generated during and/or analysed during the study are available from the corresponding author on reasonable request.

Received: 28 May 2019; Accepted: 7 November 2019;

Published online: 04 December 2019

**References**

1. van Andel, T. H. & Davies, W. *Neanderthals and modern humans in the European landscape during the last glaciation: archaeological results of the Stage 3 Project*. 265 (McDonald Institute for Archaeological Research monographs, 2003).

2. Sandgathe, D. M. et al. Timing of the appearance of habitual fire use. *Proceedings of the National Academy of Sciences of the United States of America* **108**, E298; author reply E299 (2011).

3. Aldeias, V. et al. Evidence for Neandertal use of fire at Roc de Marsal (France). *J. Archaeol. Sci.* **39**, 2414–2423 (2012).

4. Dibble, H. L. et al. How Did Hominins Adapt to Ice Age Europe without Fire? *Curr. Anthropol.* **58**, S278–S287 (2017).

5. Finlayson, C. et al. Late survival of Neanderthals at the southernmost extreme of Europe. *Nature* **443**, 850–853 (2006).

6. Finlayson, C. & Carrión, J. S. Rapid ecological turnover and its impact on Neanderthal and other human populations. *Trends Ecol. Evol.* **22**, 213–222 (2007).

7. Finlayson, C. et al. Gorham's Cave, Gibraltar—The persistence of a Neanderthal population. *Quat. Int.* **181**, 64–71 (2008).

8. Finlayson, C. On the importance of coastal areas in the survival of Neanderthal populations during the Late Pleistocene. *Quat. Sci. Rev.* **27**, 2246–2252 (2008).

9. Ochando, J. et al. Silvicolous Neanderthals in the far West: the mid-Pleistocene palaeoecological sequence of Bolomor Cave (Valencia, Spain). *Quaternary Science Reviews*, https://doi.org/10.1016/j.quascirev.2019.03.015 (2019).

10. Blain, H. A., et al MIS 5/4 transition in a mountain environment: herpetofaunal assemblages from Cueva del Camino, central Spain. *Boreas* (2014).

11. Torres, T., et al. The MIS 5 palaeoenvironmental record in the SE Mediterranean coast of the Iberian Peninsula (Río Antas, Almería, Spain). *Clim. Past Discuss.* **11** (2015).

12. Daura, J. et al. Cova del Rinoceront (Castelldefels, Barcelona): a terrestrial record for the Last Interglacial period (MIS 5) in the Mediterranean coast of the Iberian Peninsula. *Quat. Sci. Rev.* **114**, 203–227 (2015).

13. Zilhão, J. et al. Cueva Antón: A multi-proxy MIS 3 to MIS 5a paleoenvironmental record for SE Iberia. *Quat. Sci. Rev.* **146**, 251–273 (2016).

14. d'Errico, F. & Sánchez Goñi, M. F. Neandertal extinction and the millennial scale climatic variability of OIS 3. *Quat. Sci. Rev.* **22**, 769–788 (2003).

15. Blain, H.-A. et al. Climatic conditions for the last Neanderthals: Herpetofaunal record of Gorham's Cave, Gibraltar. *J. Hum. Evol.* **64**, 289–299 (2013).

16. López-García, J. M., Blain, H.-A., Bennàsar, M. & Fernández-García, M. Environmental and climatic context of Neanderthal occupation in southwestern Europe during MIS3 inferred from the small-vertebrate assemblages. *Quat. Int.* **326-327**, 319–328 (2014).

17. Daura, J. et al. Palaeoenvironments of the last Neanderthals in SW Europe (MIS 3): Cova del Coll Verdaguer (Barcelona, NE of Iberian Peninsula). *Quat. Sci. Rev.* **177**, 34–56 (2017).

18. Carrión, J. S. et al. Last Neanderthals in the warmest refugium of Europe: Palynological data from Vanguard Cave. *Rev. Palaeobot. Palynol.* **259**, 63–80 (2018).

19. Carrión, J. S. et al. The sequence at Carihuela Cave and its potential for research into Neanderthal ecology and the Mousterian in southern Spain. *Quaternary Science Reviews*, https://doi.org/10.1016/j.quascirev.2019.04.012 (2019).

20. López-García, J. M., Blain, H.-A., Bennàsar, M., Sanz, M. & Daura, J. Heinrich event 4 characterized by terrestrial proxies in southwestern Europe. *Clim. Past* **9**, 1053–1064 (2013).

21. Goñi, M. F. S., Sánchez Goñi, M. F., Eynaud, F., Turon, J. L. & Shackleton, N. J. High resolution palynological record off the Iberian margin: direct land-sea correlation for the Last Interglacial complex. *Earth and Planetary Science Letters* **171**, 123–137 (1999).

22. Goñi, M. F. S. & D'Errico, F. La historia de la vegetación y el clima del último ciclo climático (OIS5-OIS1, 140.000-10.000 BP) en la Península Ibérica y su posible impacto sobre los grupos paleolíticos. in *Neandertales cantábricos, estado de la cuestión: actas de la reunión científica: celebrada en el Museo de Altamira los días 20-22 de octubre de 2004* 115–129 (Subdirección General de Publicaciones, Información y Documentación, 2005).

23. Burjachs, F. & Julià, R. Abrupt Climatic Changes during the Last Glaciation Based on Pollen Analysis of the Abric Romani, Catalonia, Spain. *Quat. Res.* **42**, 308–315 (1994).

24. Fletcher, W. J. et al. Millennial-scale variability during the last glacial in vegetation records from Europe. *Quat. Sci. Rev.* **29**, 2839–2864 (2010).

25. Gómez-Orellana, L., Ramil-Rego, P. & Sobrino, C. M. The Würm in NW Iberia, A pollen record from Area Longa (Galicia). *Quaternary Research* **67**, 438–452 (2007).

26. Wolf, D. et al. Climate deteriorations and Neanderthal demise in interior Iberia. *Sci. Rep.* **8**, 7048 (2018).
27. Baquedano, E., Márquez, B., Laplana, C., Arsuaga, J. L. & Pérez-González, A. The archaeological sites at Pinilla del Valle (Madrid, Spain). *Pleistocene and Holocene hunter-gatherers in Iberia and the Gibraltar Strait: the current archaeological record*. Burgos: Universidad de Burgos-Fundación Atapuerca 577–584 (2014).

28. Terradillos-Bernal, M. et al. San Qirce (Palencia, Spain). A Neanderthal open air campsite with short term-occupation patterns. *Quat. Int.* **435**, 115–128 (2017).

29. Quam, R. M. et al. Human remains from Valdegoba Cave (Huérmeces, Burgos, Spain). *J. Hum. Evol.* **41**, 385–435 (2001).

30. Álvarez-Alonso, D., Andrés-Herrero, M., de, Diez-Herrero, A. & Rojo Hernández, J. El abrigo del Molino (Segovia, España). Ocupaciones neandertales en el norte del Sistema Central. VIII Reunión de Cuaternario Ibérico. Sevilla-La Rinconada. *El Cuaternario Ibérico. Investigación en el siglo XXI* 91–94 (2013).

31. Ramos-Muñoz, J. et al. Early use of marine resources by Middle/Upper Pleistocene human societies: The case of Benzú rockshelter (northern. Africa). *Quat. Int.* **407**, 6–15 (2016).

32. Arnold, L. J., Demuro, M., Navazo, M., Benito-Calvo, A. & Pérez-González, A. OSL dating of the Middle Palaeolithic Hotel California site, Sierra de Atapuerca, north-central Spain. *Boreas* **42**, 285–305 (2013).

33. Richard, M. et al. ESR/U-series chronology of early Neanderthal occupations at Cova Negra (Valencia, Spain). *Quat. Geochronol.* **49**, 283–290 (2019).

34. Fernández Peris, J. et al. The earliest evidence of hearths in Southern. *Europe: The case of Bolomor Cave* (Valencia, Spain). *Quat. Int.* **247**, 267–277 (2012).

35. Galván, B. et al. New evidence of early Neanderthal disappearance in the Iberian Peninsula. *Journal of Human Evolution* **75**, 16–27 (2014).

36. Peña, P., de la, de la Peña, P. & Wadley, L. Technological variability at Sibudu Cave: The end of Howiesons Poort and reduced mobility strategies after 62,000 years ago. *PLOS ONE* **12**, e0185845 (2017).

37. Faivre, J.-P., Gravina, B., Bourguignon, L., Discamps, E. & Turq, A. Late Middle Palaeolithic lithic technocomplexes (MIS 5–3) in the northeastern Aquitaine Basin: Advances and challenges. *Quat. Int.* **433**, 116–131 (2017).

38. Hilbert, Y. H., Crassard, R., Charloux, G. & Loreto, R. Nubian technology in northern Arabia: Impact on interregional variability of Middle Paleolithic industries. *Quat. Int.* **435**, 77–93 (2017).

39. Fernández-Peris, J., Barciela-González, V., Blasco, R., Cuartero, F. & Sañudo, P. El Paleolítico Medio en el territorio valenciano y la variabilidad tecno-económica de la Cova del Bolomor. *Treballs d'Arqueologia* 141–169 (2008).

40. Galván Santos, B., Hernández Gómez, C. M., Francisco Ortega, M. I., Molina Hernández, F. J. & Tarriño Vinagre, A. La producción lítica del Abric del Pastor (Alcoy, Alicante). Un ejemplo de la variabilidad musteriense. *Tabona: Revista de Prehistoria y de Arqueología* 11–62 (2008).

41. Machado J, Mayor A, Hernández CM, Galván B. Lithic refitting and the analysis of Middle Palaeolithic settlement dynamics: a high- temporal resolution example from El Pastor rock shelter (Eastern Iberia). *Archaeol. Anthropol. Sci.* accepted (2019).

42. Real, C. et al. Abrigo de la Quebrada Level IV (Valencia, Spain): Interpreting a Middle Palaeolithic Palimpsest from a Zooarchaeological and Lithic Perspective. *Journal of Paleolithic Archaeology*, https://doi.org/10.1007/s41982-018-0012-z (2018).

43. Picin, A. & Carbonell, E. Neanderthal mobility and technological change in the northeastern of the Iberian Peninsula: The patterns of chert exploitation at the Abric Romani rock-shelter. C. R. *Palevol* **15**, 581–594 (2016).

44. Romagnoli, F., Bargalló, A., Chacón, M. G., de Soler, B. G. & Vaquero, M. Testing a hypothesis about the importance of the quality of raw material on technological changes at Abric Romani (Capellades, Spain): Some considerations using a high-resolution techno-economic perspective. *Journal of Lithic Studies* **3** (2016).

45. Weiss, M., Otcherednoy, A. & Wiśniewski, A. Using multivariate techniques to assess the effects of raw material, flaking behavior and tool manufacture on assemblage variability: an example from the late Middle Paleolithic of the European Plain. *J. Archaeol. Sci.* **87**, 73–94 (2017).

46. Peris, J. F. Linking the first controlled use of fire in the late Middle Pleistocene at Bolomor cave (Valencia, Spain) to exceptional fuel resources. *Quaternary International* **279-280**, 140–141 (2012).

47. Mallol, C. et al. The black layer of Middle Palaeolithic combustion structures. Interpretation and archaeostratigraphic implications. *J. Archaeol. Sci.* **40**, 2515–2537 (2013).

48. Vidal-Matutano, P., Hernández, C. M., Galván, B. & Mallol, C. Neanderthal firewood management: evidence from Stratigraphic Unit IV of Abric del Pastor (Eastern Iberia). *Quat. Sci. Rev.* **111**, 81–93 (2015).

49. Vidal-Matutano, P. Anthracological data from Middle Palaeolithic contexts in Iberia: What do we know? *Munibe Antropologia-Arqueologia*, https://doi.org/10.21630/maa.2018.69.12 (2018).

50. Eixea, A., Giner, B., Jardón, P., Zilhão, J. & Bonilla, V. V. Elementos líticos apuntados en el yacimiento del Paleolítico Medio del abrigo de La Quebrada (Chelva, Valencia): caracterización tecno-tipológica y análisis de las macrofacturas. *Espacio Tiempo y Forma. Serie I, Prehistoria y Arqueología* 0, 153 (2015).

51. Sañudo, P., Blasco, R. & Peris, J. F. Site formation dynamics and human occupations at Bolomor Cave (Valencia, Spain): An archaeostratigraphic analysis of levels I to XII (100–200 ka). *Quaternary International* **417**, 94–104 (2016).

52. Villaverde, V., Martinez-Valle, R., Guillem, P. M. & Fumanal, M. P. Mobility and the role of small game in the Middle Paleolithic of the central region of the Spanish Mediterranean: a comparison of Cova Negra with other Paleolithic deposits. *The Last Neandertals the first anatomically modern humans. Cultural change and human evolution: The crisis at 40, 267–288 (1996).*

53. Molina, F. J., Tarriño, A., Galván, B. & Hernández, C. M. El sílex del Prebético de Alicante: tipos, variabilidad y áreas de captación y talla del Pleistoceno. *Cuadernos de Prehistoria y arqueología de la Universidad de Granada* **26**, 283–311 (2016).

54. Molina, F. J. Estudio geoarqueológico de entornos sedimentarios fluvio-lacustres y endorreicos con industrias del Paleolítico medio en el norte de la provincia de Alicante (España). *Recerques del Museu d'Alcoi* **25**, 7–29 (2016).

55. Machado, J. & Pérez, L. Temporal frameworks to approach human behavior concealed in Middle Palaeolithic palimpsests: A high-resolution example from El Salt Stratigraphic Unit X (Alicante, Spain). *Quat. Int.*, https://doi.org/10.1016/j.quaint.2015.11.050 (2015).

56. Leerer, L. et al. Insights into the timing, intensity and natural setting of Neanderthal occupation from the geoarchaeological study of combustion structures: A micromorphological and biomarker investigation of El Salt, unit Xb, Alcoy, Spain. *PLOS ONE* **14**, e0214955 (2019).

57. Machado, J., Hernández, C. M., Mallol, C. & Galván, B. Lithic production, site formation and Middle Palaeolithic palimpsest analysis: in search of human occupation episodes at Abric del Pastor Stratigraphic Unit IV (Alicante, Spain). *J. Archaeol. Sci.* **40**, 2254–2273 (2013).

58. Machado, J., Mallol, C. & Hernández, C. Insights into Eurasian Middle Palaeolithic Settlement Dynamics: The Palimpsest Problem. in *Settlement Dynamics of the Middle Paleolithic and Middle Stone Age* (ed. Nicholas J. Conard, A. D.) IV, (Kerns Verlag, 2015).

59. Morales Pérez, J. V. & Serra, A. S. The Quaternary fossil record of the genus Testudo in the Iberian Peninsula. Archaeological implications and diachronic distribution in the western Mediterranean. *J. Archaeol. Sci.* **36**, 1152–1162 (2009/5).

60. Sanchis, A., Morales, J.V., Pérez, L.J., Hernández, C.M., Galván, B. La tortuga mediterránea en yacimientos valencianos del Paleolítico medio: distribución, origen de las acumulaciones y nuevos datos procedentes del Abric del Pastor (Alcoi, Alacant). In *Preses Petites i Grups Humans En El Passat. II Jornades d'arqueozologia* (ed. Sanchis, A. Pasual, J.) 97–120 (Museu Prehistòria de València, 2015).

61. Connolly, R. et al. A multiproxy record of palaeoenvironmental conditions at the Middle Palaeolithic site of Abric del Pastor (Eastern Iberia). *Quat. Sci. Rev.* **225**, 106023 (2019).
62. Eglinton, G. & Hamilton, R. J. Leaf Epicuticular Waxes. *Science* **156**, 1322–1335 (1967).

63. Knicker, H., Hilscher, A., de la Rosa, J. M., González-Pérez, J. A. & González-Vila, F. J. Modification of biomarkers in pyrogenic organic matter during the initial phase of charcoal biodegradation in soils. *Geoderma* **197-198**, 43–50 (2013).

64. Jambrina-Enríquez, M., Herrera-Herrera, A. V. & Mallol, C. Wax lipids in fresh and charred anatomical parts of the Celtis australis tree: Insights on paleofire interpretation. *Org. Geochem.* **122**, 147–160 (2018).

65. Ficken, K. J., Barber, K. E. & Eglinton, G. Lipid biomarker, δ13C and plant macrofossil stratigraphy of a Scottish montane peat bog over the last two millennia. *Organic Geochemistry* **28**, 217–237 (1998).

66. Diefendorf, A. F., Leslie, A. B. & Wing, S. L. Leaf wax composition and carbon isotopes vary among major conifer groups. *Geochim. Cosmochim. Acta* **170**, 145–156 (2015).

67. Chabot, J., Dionne, M.-M. & Paquin, S. High magnification use-wear analysis of lithic artefacts from Northeastern America: Creation of an experimental database and integration of expedient tools. *Quaternary International* **427**, 25–34 (2017).

68. Kirgesner, S. L. *Analysis of Microwear on Flakes Used to Butcher Unfrozen and Frozen Meat*. (Kent State University, 2018).

69. Vidal-Matutano, P., Firewood and hearths: Middle Palaeolithic woody taxa distribution from El Salt, stratigraphic unit Xb (Eastern Iberia). *Quaternary International* **457**, 74–84 (2017).

70. Sandgathe, D. M. Identifying and Describing Pattern and Process in the Evolution of Hominin Use of Fire. *Curr. Anthropol.* **58**, S360–S370 (2017).

71. Conard, N. J. & Adler, D. S. Lithic Reduction and Hominid Behavior in the Middle Paleolithic of the Rhineland. *Journal of Anthropological Research* **53**, 147–175 (1997).

72. Lee Lyman, R. & Lyman, C. *Vertebrate Taphonomy*. (Cambridge University Press, 1994).

73. Reitz, E. J. & Wing, E. S. *Zooarchaeology by Elizabeth J. Reitz*, https://doi.org/10.1017/CBO9780511841354 (Cambridge University Press, 2008).

## Acknowledgements

This research was funded by a Leakey Foundation General Grant, Spanish Ministry of Science, Innovation and Universities Projects HAR2008-06117/HIST and HAR2015-68321-P, Junta de Castilla y León-FEDER Project BU235P18, the LabEx Sciences Archéologiques de Bordeaux (LaScArBx ANR-10-LABX-52) and ERC Consolidator Grant ERC-CoG-2014. Archaeological excavations at Abric del Pastor are supported by the Archaeological Museum of Alcoy and the Government of Valencia Cultural Heritage Department. The authors would like to thank the members of the Abric del Pastor excavation team for their fieldwork through the years, Caterina Rodríguez de Vera for thin section manufacture, Antonio Herrera-Herrera and Margarita Jambrina-Enríquez for laboratory guidance and assistance during biomarker analysis and Santiago Sossa Ríos for assistance in lipid extraction and injection of sediment samples.

## Author contributions

C.M., C.H. and B.G. are responsible for archaeological excavations at Abric del Pastor and the research design. N.M. performed OSL dating, C.F. and E.B.A. performed ESR dating, C.M. and R.C. performed micromorphology, R.C. performed lipid analysis, A.C. performed archaeomagnetic analysis, L.P. performed zooarchaeological analysis, P.V.-M. performed anthracological analysis, C.H., B.G. and A.M. performed lithic analysis and D.C. performed FTIR analysis. All authors wrote the main manuscript text and prepared the figures. All authors reviewed the manuscript.

## Competing interests

The authors declare no competing interests.

## Additional information

Supplementary information is available for this paper at https://doi.org/10.1038/s41598-019-54305-9.

Correspondence and requests for materials should be addressed to C.M.

Reprints and permissions information is available at www.nature.com/reprints.

Publisher's note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_10.png" /></figure>

**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 license, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons license 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 license, visit http://creativecommons.org/licenses/by/4.0/.

© The Author(s) 2019