Socio

Scientific Reports volume 13, Article number: 8168 (2023) Cite this article

1046 Accesses

231 Altmetric

Metrics details

Volcanism can cause major impacts, including climate change and mass extinctions. However, the impact of monogenetic volcanism is often considered as limited in volcanological research. This work provides for the first time an interdisciplinary approach to the socio-ecological impact of monogenetic volcanism in a key region, the La Garrotxa Volcanic Field (GVF, Girona, NE Iberia), where intense monogenetic volcanic activity occurred in the past. The analyses of a sedimentary sequence from the GVF enabled identifying previously unknown volcanic eruptions in the time interval 14–8.4 ka cal BP, constrain their volcanic stratigraphy and age, and unfold the effects of environmental change on geomorphology, vegetation, aquatic organisms and humans. Moreover, we reconstruct the major palaeoenvironmental changes caused by the eruptions in terms of fire episodes and subsequent disturbance on vegetation, hydrology and limnological conditions. When put in context with the archaeological record, it appears that the last hunter–gatherer communities were resilient at an extra-local scale, facing episodes of vulnerability due to volcanic activity, suggesting that their flexible nomadic patterns and foraging economies were an efficient source of risk management against the volcanic eruptions and their ecological impacts.

There is ample evidence to support the view that volcanic eruptions have caused in the past major impacts on environment and societies, both directly (e.g. lava flows, tephra deposition and earthquakes1,2,3), as well as indirectly (e.g. contributing to climate change4,5,6, mass extinctions7,8 and disturbance on human communities9,10,11). However, volcanological and related palaeoecological research has been mainly focused on explosive eruptions from polygenetic volcanoes (with repeated eruptions), such as Laacher See (Germany1), Santorini (Greece12), Somma-Vesuvius and Campanian Ignimbrite (Italy13,14), Mount Pinatubo (Philipines15) or Toba and Krakatoa (Indonesia16), but the impact of monogenetic (single small basaltic eruptions) volcanism has been generally neglected or considered as limited.

The La Garrotxa monogenetic Volcanic Field (GVF) is the youngest (Middle Pleistocene (ca. 350 ka BP) to early Holocene) volcanic area of the Iberian Peninsula. Recent studies in this area have addressed the characterization of the volcanic activity and petrology of volcanic products17,18,19,20,21 and the geological and structural controls of the magmatism and related volcanism22,23,24. Still, the impact of these eruptions on the environment and past human societies has not yet been addressed. Past human communities in volcanic areas have been repeatedly threatened by eruptive activity and subjected to short-term catastrophic events leading to major landscape changes13. It becomes crucial to consider the socio-ecological characteristics of the human societies, since their vulnerability or resilience relies on different aspects such as their settlement patterns, demography and socio-political organization and economic activities25. In fact, volcanism does not necessarily drive to disaster and societal collapses deterministically, as there are examples of resilient societies that coped with volcanic events26. In this context, the mobility of hunter–gatherers may have been an effective strategy to deal with the consequences of environmental constraints27,28,29. To understand past hunter–gatherer interaction with volcanic eruptions is important to consider post-eruption resilience measures but also the geological parameters of the eruption30.

Successive lava flows in the La Garrotxa Volcanic Field dammed the Fluvià river and led to the formation of a lacustrine basin in the deepest part of the Vall d’en Bas valley (La Garrotxa, NE Iberia) (Fig. 1). The youngest lava flows were likely associated with the most recent eruptions from the Puig Jordà (17 ka BP22) and the Croscat (15.7–13.2 ka cal BP31) volcanoes. However, the ages of the eruptions are not well constrained as they have high uncertainty and low reliability. Recently, the Pla de les Preses sediment succession, at the Vall d’en Bas valley (Fig. 1) enabled precise radiocarbon dating of macrofossils and bulk sediment associated with tephra (Suppl. Files 1), identifying several eruptions in the period 14.0–8.4 ka cal BP, showing the potential of barrier-lake deposits to date volcanic eruptions. These lacustrine sediments also offered the possibility of conducting detailed palaeoenvironmental reconstructions involving regional landscape (vegetation) and local lacustrine environments (aquatic organisms). In that sense, the objective of this paper is to show the potential socio-ecological impact of monogenetic volcanism, integrating geochemical (XRF) and palaeobiological (pollen and non-pollen palynomorphs, sedimentary charcoal, ostracods, charophyte gyrogonites and diatoms) proxies and the archaeological record in the region, to assess how volcanism affected last hunter–gatherer communities and their surrounding environment.

Location of La Garrotxa Volcanic Field in NE Iberia. Top left, archaeological sites and palaeoecological cores mentioned in the text are indicated. In the main map (right), lava flows and volcanos for the last 20,000 years are indicated. Map created with QGIS (https://www.qgis.org/es/site/), version 3.22, and modified in Adobe Illustrator version 23.0.5.

The monogenetic Quaternary La Garrotxa Volcanic Field18,22 forms part of the Catalan Volcanic Zone (north-east Iberian Peninsula)32,33, one of the alkaline volcanic provinces of the European rift system. It harbours more than 50 basaltic monogenetic cones that range in age from the Middle Pleistocene (ca. 350 ka BP) to the early Holocene and include cinder and scoria cones, lava flows, tuff rings and maars.

The 15-m long sediment core from Pla de les Preses (Fig. 1) provides palaeoenvironmental data for the last 14,000 years (Suppl. Files 1), but we focus on the time interval 14–8 ka cal BP, since palaeobiological studies had richer results in lacustrine and wetland sediments in that period, and volcanic activity was not detected in fluvial sediments from 8.4 ka cal BP onwards (Suppl. Files 2). At the beginning of the sedimentary record, fluvial sedimentation was disrupted by coarse sands and volcanic tephra, linked with volcanic activity at ca. 14.0 ka cal BP. Lava flows from this eruption (by the Puig Jordà volcano, Fig. 1) would have built a barrier on the Fluvià river enabling the formation of a lake. The lacustrine environment lasted from 13.6 to 9.3 ka cal BP, thus spanning the time interval from the Late-Glacial (Bølling-Allerød and Younger Dryas) to the early Holocene. Several tephra layers (Suppl. Files 3) were deposited during a phase of more intense volcanic activity at 13.0–12.0 ka cal BP and another eruptive episode at 10.4 ka cal BP (Fig. 2). A shallowing process was accentuated by a cool and dry episode (Bond event 6)34, leading to the transition from lacustrine to palustrine conditions at 9.3 ka cal BP (Fig. 2). Tephra layers are also identified intercalated between these palustrine peaty layers, which correspond to the most recent eruptions at 9.4–8.4 ka cal BP (Figs. 2, 4). Around 8.2 ka cal BP, wetland areas converted to a fluvial floodplain again, suggesting a shallower water body and the incision of the volcanic dam by the river (Suppl. Files 2).

Paleoenvironmental evolution in Vall d’en Bas valley based on the main geochemical PCA and aquatic organisms (pollen, non-pollen palynomorphs, ostracods, gyrogonites and diatoms). The main environmental phases are indicated in colour frames (Green: palustrine before lake formation; Light blue: shallow lake; Dark blue: deep lake; Grey: lake margin; Grey to pale brown: palustrine; Orange: fluvial floodplain). The red bars indicate volcanic tephra layers.

Pollen analysis provided data about vegetation history from 14.0 to 8.0 ka cal BP. During the Late-Glacial, the prevailing landscape ecosystems were steppes and grasslands (40–60%) and forests were dominated by Pinus and woodlands of Betula, Acer and Juniperus. This characteristic landscape is consistent with cool climate conditions during the Late-Glacial. However, the transition from the warmer Bølling-Allerød to the colder Younger Dryas period did not cause substantial change of vegetation (Fig. 3), confirming a milder Younger Dryas in NE Iberia35, as attested previously in the Central Pyrenean region36,37,38. After the first evidence of links between eruption and fire episodes in 14.0 and 13.5 ka cal BP, the major Late-Glacial phase of intense fire and volcanic activity occurred during 13.0–12.0 ka cal BP (Fig. 3), confirming the significant role of volcanic activity in initiating fire episodes39, among other factors such as dry Late-Glacial climatic conditions, as attested in other Iberian records37,40. Volcanic eruptions and subsequent forest fires caused short-term episodes of expansion of pine forest and decline of Late-Glacial woodlands (Betula-Acer-Juniperus, Fig. 3) and steppes. However, the Late-Glacial vegetation was not disrupted dramatically, showing recovery during the following 50–100 years after the eruptions (Fig. 4).

Percentage pollen diagram of selected taxa and categories. Data is plotted on an age scale (cal BP) and tephra layers were excised. Categories: broadleaf deciduous (Quercus deciduous, Corylus, Ulmus, Fagus), steppe (Helianthemum-t, Sanguisorba, Artemisia, Amaranthaceae, Plantago, Apiaceae, Rumex, Galium-t, Filipendula). In the broadleaf deciduous curve, the discontinuous line shows values in Les Palanques core (2 km distance from PdP) and the grey line shows SB2 Banyoles core (25 km away). Sedimentary charcoal data is plotted as accumulation rate of charcoal particles, n of particles > 1 mm width and 0.5–0.99 mm width, and charcoal peaks identified by numerical analyses. Red shadings indicate periods of intense eruptive and fire activity, grey lines show events of tephra deposition. Greenland isotopic curves (GRIP and GISP2)42,43 are plotted and cooling episodes/cold phases (blue) and warmer periods (orange) are indicated.

Multi-proxy diagram showing environmental shifts in relation to selected tephra. Data is plotted on a depth scale (cm) and grey shadings indicate tephra deposition events. Occupation layers from Late Mesolithic sites in the region (15–40 km from the La Garrotxa Volcanic Field) are shown in the top-right corner: Sota Palou (10,200–9100 cal BP), Bauma del Serrat del Pont IV.1 (9400–9100 cal BP), IV.2 (9100–8800 cal BP), IV.3 (8600–8400 cal BP), Bauma dels Fadrins (8700 cal BP).

The Pinus and Betula-Acer-Juniperus dominated Late-Glacial woodlands, as well as the grasslands and steppes, were replaced by broadleaf deciduous forests (mainly Quercus and Corylus) due to a warmer and wetter early Holocene climate between 11.7 and 10.5 ka cal BP41,42 (Fig. 3). This rapid expansion of broadleaf deciduous forest, which reached maximum values around 10.3–9.2 ka cal BP (Fig. 3) was interrupted by volcanic activity, combined with the cooling Bond Event 6 (9.3 ka cal BP). The first evidence of early Holocene volcanism is observed at 10.4–10.3 ka cal BP, provoking fire episodes and a limited impact on deciduous woodlands (see tephra 13 in Fig. 4). Later, at 9.4–8.3 ka cal BP the broadleaf deciduous forests collapsed in conjunction with more frequent fire episodes and volcanic eruptions. This suggests that fire episodes induced by volcanic eruptions dramatically impacted on broadleaf deciduous forests, the predominant vegetation in this region during the early and middle Holocene43, triggering the expansion of grasslands and of secondary Pinus-dominated woodlands. Although palaeoecological evidence attested landscape transformation by the use of fire by last hunter–gatherer communities in Western and Central Europe44,45,46, the characteristic sparse Mesolithic populations in NE Iberia43 would have had a limited impact on the landscape. The detailed records for eruptions at 9.15–9.1 ka cal BP and 8.84–8.77 ka cal BP (tephra 17 and 19 in Fig. 4) show major impacts of eruptions and wildfires in a two-step process: firstly, dramatic impact on broadleaf deciduous forests, followed by pine woodlands expansion; secondly, the burning of pine forest and the expansion of clearings dominated by grasslands (Poaceae) (Fig. 4) and a moderate expansion of disturbance-sensitive trees (Abies and Tilia47) when the frequency of volcanic eruptions decreased and climatic conditions became warmer (Fig. 3). Recurrent eruptive activity and more intense wildfires prevented deciduous forests to rapidly recover. Nevertheless, a previous study in the Vall d’en Bas valley (Les Palanques, 2 km to the west, Fig. 1) showed that broadleaf deciduous forests declined in conjunction with peaks in Pinus and Poaceae at 9.5–9.0 and 8.5–8.3 ka cal BP and recovered from 8.1 ka cal BP onwards43,48 (Fig. 2), pointing to a significant impact of volcanic activity on vegetation at the local scale and slower recovery processes during the early Holocene than during the Late-Glacial. Likewise, an expansion of Pinus and Poaceae and decline in broadleaf deciduous forests is documented around 8.8–8.6 ka cal BP (Fig. 2) at extra-local scale, 15 km away, in Bauma del Serrat del Pont (highest values of Pinus and Poaceae in Mesolithic layers49,50), and 25 km away, in Lake Banyoles51 (Fig. 1). However, the pollen record from Lake Banyoles clearly shows maxima of broadleaf deciduous vegetation in conjunction with low fire activity during the period 9.0–7.5 ka cal BP51, when climatic conditions were wetter and warmer (Holocene Climate Optimum). Overall, the evidence of dramatic environmental changes from the Pla de les Preses records support the view that local factors (e.g. volcanic activity) rather than large-scale factors such as climate drove vegetation and fire dynamics in the Vall d’en Bas valley. Cooling episodes during the Younger Dryas and the Holocene (10.3 and 9.3 kyr cal BP) acted as amplifiers, configuring drier landscapes that are more prone to wildfire occurrence and spread.

The multi-proxy approach developed in this work clearly evidenced the impact of volcanic eruptions on limnological conditions, as shown by the abrupt declines in aquatic plants (genus Myriophyllum in the Late-Glacial, genus Nymphaea in the early Holocene), algae (Botryococcus, Tetraedron, charophytes, diatoms), cyanobacteria and ostracods (Candona candida, Ilyocypris gibba, Neglecandona neglecta) after eruptive episodes (Fig. 2). During the Late-Glacial, the deposition of volcanic tephra, charcoal and ash from wildfires altered water conditions, resulting in lower alkalinity, as inferred by the predominance of the acid-resistant algae Botryococcus52 in deep water phases and the predominance of the aquatic plant Myriophyllum in shallower phases (tephras 5 & 7 and 4 & 6, respectively, Fig. 4). Limnological changes imposed by volcanism limited life on aquatic organisms in the short term, but Gloeotrichia played a significant role as pioneer in phases of poor nutrient availability53, spreading fast perhaps due to fixing nitrogen and the high availability of phosphorus. These conditions were also favourable for Myriophyllum, a macrophyte growing in waters with high nitrogen content and phosphorus availability54 and that prefers acidic waters55 (shallow lake phase in Unit 2a, Fig. 2). Peaks in cyanobacteria and Chlorophyta show how the lake trophism changed from oligotrophic, immediately after volcanic eruptions, to eutrophic in periods of 50 to 100 years (Tephra 5, 6, 7 in Fig. 4). Similarly, the planktonic diatom Lindavia radiosa was only observed after the volcanic eruption at 10.24 ka cal BP, which indicates an increase in the lake trophic state and the presence of a water column. However, the most abundant benthic diatoms in the pre-tephra sample (10.42 ka cal BP) were still abundant after the volcanic eruption sample (10.24 ka cal BP), despite the scarcity of diatoms observed in the tephra layer (Suppl. Mat. 2). These results indicate a rapid recovery of diatom assemblages and other aquatic organisms after the 10.35 ka cal BP eruption (Fig. 2).

Therefore, volcanism enhanced water acidity and aquatic organisms adapted to these new conditions. However, these disturbance events were followed by fast recovery processes, especially during the Late-Glacial, when the deeper lake showed high resilience. Nevertheless, during the early Holocene, the combination of a shallowing process enhanced by climate change, volcanic activity and high-intensity local wildfires affected local aquatic communities, which did not recover after the eruptive activity at 9.3–8.8 ka cal BP (Figs. 2, 4).

The archaeological record shows a gap of human settlement in the GVF area since the Upper Palaeolithic (La Rodona, Olot, ca. 33–24 ka cal BP56) and human communities did not settle again until the Neolithic (Codella and La Dou, 6.7 ka cal BP57,58), suggesting that this area was hostile for the last hunter–gatherers during the Late-Glacial and early Holocene. In this context, it is worth mentioning that the low Mesolithic population density is not exclusive to the GVF area. This period has not attracted especial attention in archaeological research in the NE Iberian Peninsula given the scarcity of archaeological records with evidence of occupations during the early Holocene, as opposed to neighbouring areas (Ebro basin, eastern coast of Iberia and Mediterranean coastal area of France)59. The scarcity of archaeological evidence in this region of NE Iberia during the Late-Glacial and the early Holocene is probably not driven by climate, as the archaeologically inferred gap of occupations is not restricted to colder phases (e.g. the Younger Dryas). Instead, the occupation gap could have been driven by settlement patterns or post-depositional processes affecting the preservation of archaeological sites. Mesolithic occupations are documented 15 km away from the GVF, in the Llierca valley (Fig. 1). There, the rock shelter site Bauma del Serrat del Pont (BSP) (Fig. 1) shows 4 layers of Mesolithic occupations, from 9.4 to 8.0 ka cal BP (9.45–9.1, 9.1–8.8, 8.6–8.4, 8.3–8.0 ka cal BP)60, contemporary to the most frequent early Holocene volcanic eruptions in the GVF area. These communities were located in ecotones providing diverse resources for their terrestrial foraging economies, including many wild hunted species such as red deer, roe deer, wild boar and ibex61 and plants and fruits such as acorns, sorb apples, hazelnuts and strawberry tree fruits62. Overall, the coincidence of eruptive episodes and changes in archaeological layers (Fig. 4) suggests that hunter–gatherer communities at BSP abandoned temporarily the site at 9.1 (transition from layers IV.1 to IV.2) and 8.8 ka cal BP (200 years of gap between layers IV.2 and IV.3) in response to enhanced volcanic activity. Different hazard agents likely acted in the proximal impact zone (up to at least 50 km), including lava flows, tephra deposition, ash storms, gases, aerosols, pyroclastic flows and earthquakes30, affecting flora and fauna (bioresources) as well as the quality of air and water. In that context, other Mesolithic sites were occupied during this period of environmental disturbance in the region, including Sota Palou (10.2–9.1 ka cal BP63) and Bauma dels Fadrins (8.7 ka cal BP64) (28 and 38 km from the GVF, respectively), suggesting a spatially restricted impact of the last volcanic eruptions in GVF. BSP was re-settled at 8.6–8.4 and 8.3–8.0 ka cal BP, suggesting that the overall social system was not challenged, and Mesolithic communities did not collapse. Hunter–gatherer societies dwelling in the proximal impact zone may have abandoned the area temporarily during high volcanic activity periods, but returned afterwards, proving high reorganization capacity. Thus, the low magnitude monogenetic volcanic eruptions in GVF did not have such a dramatic effect as larger magnitude explosive volcanic eruptions causing the collapse on past hunter–gatherer communities30. What emerges from the analysis of the archaeological and palaeoecological records in the La Garrotxa region is that the last hunter–gatherer communities were resilient at an extra-local scale (15 to 40 km from the GVF) against the local-scale impact of monogenetic volcanism, which indeed affected the settlement in the Vall d’en Bas valley since the Upper Palaeolithic to the early Neolithic (gap of human occupations in 24–6.7 ka cal BP at local scale). Their flexible nomadic strategies and foraging economies were likely an efficient source of risk management to face episodes of vulnerability caused by volcanic eruptions. Mobility was the risk mitigation strategy enabling the resilience of past hunter–gatherers in NE Iberia against volcanic eruptions, as there is no clear evidence in the archaeological record of other possible crisis management strategies such as storage, exchange, diversification or intensification30.

This work provides insights on how environmental disturbance by monogenetic volcanism influenced geomorphology, vegetation, aquatic organisms and past human societies. New Late-Glacial-early Holocene volcanic eruptions, not previously reported in GVF, are presented, and their volcanic stratigraphy and age are constrained, reconstructing their major palaeoenvironmental impacts (Fig. 5). Volcanic activity had a significant environmental impact producing intense fires and causing dramatic changes in landscape vegetation on a local scale. Volcanism initiated intense wildfire episodes affecting woodlands (Acer-Betula-Juniperus in the Late-Glacial; Corylus-Quercus deciduous in the early Holocene), except for pine forests, which expanded colonising the disturbed areas. This study shows that Late-Glacial cold steppes recovered quicker from volcanic disturbance than the Holocene broadleaf deciduous forests. In addition, the records indicate that tephra deposits also altered the lake ecosystem, enhancing water acidity and boosting the aquatic organisms adapted to these conditions. These disturbance events were followed by fast recovery at decadal to centennial scales, especially during the Late-Glacial, when the deeper lake showed higher resilience.

Synthetic diagram including climate change events, sedimentary evolution, volcanic eruptions (tephra), fire episodes, charcoal accumulation rate, vegetation dynamics and human settlement during Late Mesolithic in the region.

This study enables a better understanding of palaeoenvironmental dynamics and ecological changes that occurred during the youngest volcanic eruptions in GVF. Volcanism had a large socio-ecological impact at local scale (radius < 15 km) causing intense fire activity and remarkable changes in vegetation (forests), aquatic communities (lake) and a gap of human settlements until the Neolithic (6.7 ka cal BP). Most importantly, the new records add important evidence to explain the Late-Mesolithic population resilience against volcanism at an extra-local scale, using mobility as a source to face vulnerability episodes, since hunter–gatherer groups abandoned (during high volcanic activity periods) and re-occupied (during volcanic quiescence) the BSP site four times during the period 9.4–8.0 ka cal BP. Thereby, this study proves the socio-ecological impact and the interest in developing interdisciplinary palaeoecological research to address local-scale landscape transformations caused by the most frequent volcanism on Earth, monogenetic volcanism.

A 15-m long core was obtained using a mechanical rotary drilling machine (TP-50/D) from the lowest part of the valley of La Vall d’en Bas (Girona, Spain) (UTM 455189.0 X/4667356.0 Y/458.1 m asl), in the area known as Pla de les Preses (PdP).

The depth-age model was established with RBacon65 and uses 15 control points. Thirteen of these points are based on radiocarbon (14C) dates measured on bulk sediments and two points are based on 14C dates measured on terrestrial plant remains (Pinus sp. needles and seeds) (Suppl. Files 1). As the latter 14C dates very likely provide more accurate sediment-deposition ages66, narrower student‐t error distributions were applied were applied for plant-macrofossil dates. Abrupt sedimentation events (e.g. tephra deposits) were excised, and a sedimentation hiatus was set at 945 cm depth with a maximum hiatus length of 1500 years. The hiatus, which was probably caused by poor sedimentation recovery during the coring, led to the absence of sediments for the period between 11.7 and 10.5 ka cal BP.

The lithostratigraphic study of the core was defined taking into account the different sedimentary facies (Suppl. Files 2). Sedimentary facies were defined by visual macroscopic description and microscopic observation of smear slides following LRC procedures67 (and by mineralogical, organic and geochemical compositions). Different stratigraphic units were defined along the core and their depositional environments and processes inferred based on their sedimentological characteristics68.

A high-resolution geochemical analysis (1 cm step size) of the core was performed using an Avaatech XRF Core-scanner at the Corelab Laboratory (University of Barcelona). The analysis was performed using a Rhodium source under two different working conditions: 1) with an X-ray current of 800 μA, at 10 s count time and 10 kV x-ray voltage for the measurement of Al, Si, P, S, Cl, Ar, K, Ca, Ti, V, Rh, Cr, Mn and Fe; 2) with an X-ray current of 2000 μA, at 25 s count time, 30 kV x-ray voltage and using a Pd filter, for the measurement of Ni, Cu, Zn, Ga, Ge, As, Br, Rb, Sr, Y, Zr, Nb and Pb. This method allowed a semi-quantitative analysis of the elemental chemical composition from Al to U, based on the proportion of counts per second (cps) for each element compared to the rest. The most abundant and significant elements (Al, Si, Cl, K, Ca, Ti, V, Mn, Fe, Ni, Cu, Zn, Ga, Br, Rb, Sr and Pb) were selected for multivariant statistical analysis (PCA) to reduce the number of variables and define the main phases and processes involved in the formation of the core record according to their chemostratigraphy. Before the PCA analysis all unreliable measurements were removed, so as not to obscure the statistical treatment of the data. XRF geochemical data was normalized using centred log-ratio transformation69,70 using CoDaPack software71 and processed with multivariate statistics. A Principal Component Analysis was performed using SPSS 23.0 software at correlation mode, factor scores were calculated, and rotated (Varimax) and not rotated solutions were evaluated and the most suitable to geochemical data variance selected (Suppl. Files 2).

Pollen samples were obtained each 3–5 cm in organic clayish and peaty facies and each 10 cm in inorganic silts in the fluvial layers. Samples were processed following standard methods72 including treatment with HCl and NaOH, flotation in Thoulet heavy liquid, treatment in HF, and finally mounting in glycerine. At least 300 pollen grains of terrestrial taxa were counted using an Olympus Bx43 microscope fitted with 10× oculars and 40/60× objectives. Hygrophytic plants (Cyperaceae, Ranunculaceae Typha latifolia and Typha/Sparganium) and aquatic plants (Myriophyllum, Nuphar, Nymphaea, Potamogeton) were excluded from the pollen sum. Pollen grains were identified using a pollen atlas73. The identification of non-pollen palynomorphs (NPPs) followed van Geel74, van Geel et al.75, Revelles et al.76 and Revelles and van Geel77. Pollen percentages were calculated with respect to the pollen sum and diagrams plotted using Tilia software78. Despite pollen analysis was applied to the whole succession, some samples were poor in pollen for the top 650 cm, and data is only provided for the bottom part of the core (14.0–8.0 ka cal BP).

Quantification of charcoal particles was performed with the sieving method79 with a 150 μm mesh size80 in order to identify fire episodes. Samples of 1 cm3 were retrieved each 1 cm from the whole sediment sequence. Samples were first soaked in 10% H2O2 for 12 h to deflocculate and bleach the sediment and then sieved on a 150 μm sieve under a soft-water jet. Very organic samples were additionally soaked in 10% NaOCl for 4 h to further bleach the organic material. The bleached sieving residue was analysed under a stereoscopic microscope (Leica M80 at 60x) equipped with a camera CMEX DC 5000 connected to a computer with an image-analysis software (WinSeedle, Regent Instruments Canada, Inc.) that allowed the measurement of charcoal concentration, charcoal areas of individual particles and the cumulative sum of charcoal-particle areas81. The CharAnalysis software82 was used to calculate charcoal accumulation rates (charcoal cm−2 yr−1) and to detect fire episodes. The analysis was run into two time-windows, one for Holocene samples (8.2–10.5 ka cal BP), one for Late Glacial ones (11.7–14.0 ka cal BP). Charcoal counts, sample volume and sample depths were interpolated to a constant temporal resolution of 10/6 (Holocene/Late Glacial) years before calculating charcoal accumulation rate to account for unequal sampling intervals resulting from variable sediment accumulation rates. No data transformation was performed before applying numerical analysis. The slowly varying mean or background component (Cback) was modelled through a LOWESS function with two different smoothing windows (400 yrs in Holocene, 300 yrs in Late Glacial) based on the highest values of goodness of fit test. Finally, two categories of charcoal particles width were plotted (Fig. 3) to explore the local signal of fires.: > 1 mm width and 0.5–0.99 mm width.

Some macrofossils were recovered within sedimentary charcoal samples. Despite organic macrofossils were affected by H2O2 processing, most seeds could be identified using a stereoscopic microscope (Leica M80 at 60x). Identifications were made based on literature83,84 and the reference collection of seeds at the University of Montpellier).

Samples of around 20 g of sediment were retrieved for ostracods and charophyte analyses every 10 cm in the peaty and lacustrine facies (600–1470 cm). The samples were rinsed in water (with use of H2O2 to disaggregate clayish samples) and sieved through 250 μm. Finally, the samples were dried and all ostracod remains and charophyte gyrogonites were picked up with a fine brush.

All ostracod remains (shells and disarticulated valves) were identified to species level whenever possible, following mainly Meisch85 and Fuhrmann86. Densities were estimated as the number of valves per gram of dry sediment. The qualitative taxonomic characteristics used to identify the charophyte gyrogonites were: apical zone, basal structures (presence or absence of a basal column, shape of the basal plate) as well as other features such as overall outline and number of spiral turns (or ridges) visible in lateral view87. The observation and measurements were made with a stereomicroscope at 400x. Length was measured as the longest polar axis (LPA = vertical axis); width as the largest equatorial diameter (LED = horizontal axis at the largest diameter). The length/width ratio was also calculated and expressed as the isopolarity index (ISI = LPA/LED × 100).

Twelve samples were selected for diatoms analysis before and after 4 tephra layers to assess the impact and the recovery of the diatoms community to tephra deposition in the lake ecosystem. All these samples were treated with 33% hydrogen peroxide (H2O2) and HCL (1 M). Subsequently samples were mounted in Naphrax (R.I. = 1.7) following the method described in Battarbee et al.88. Diatom identification was followed mainly Krammer and Lange-Bertalot89,90,91,92, but the diatoms nomenclature (basionym) were updated to accepted names following currently accepted nomenclature93. Unfortunately, diatom preservation was sufficient only in two Holocene samples (before and after T6) to identify and count enough, at least 10 diatom valves per slide. Diatom dissolution, and therefore their absence in sediment records, could result from an increase in lake salinity or alkalinity94.

The datasets generated during the current study are available in the Neotoma Palaeoecological Database repository (https://data.neotomadb.org/56691).

Birks, H. J. B. & Lotter, A. F. The impact of the Laacher See Volcano (11000 yr BP) on terrestrial vegetation and diatoms. J. Paleolimnol. 11(3), 313–322 (1994).

Article ADS Google Scholar

Allen, J. R. & Huntley, B. Effects of tephra falls on vegetation: A Late-Quaternary record from southern Italy. J. Ecol. 106(6), 2456–2472 (2018).

Article Google Scholar

Antoniades, D. et al. The timing and widespread effects of the largest Holocene volcanic eruption in Antarctica. Sci. Rep. 8, 17279 (2018).

Article ADS PubMed PubMed Central Google Scholar

Robock, A. Volcanic eruptions and climate. Rev. Geophys. 38(2), 191–219 (2000).

Article ADS CAS Google Scholar

Baldini, J. U., Brown, R. J. & McElwaine, J. N. Was millennial scale climate change during the Last Glacial triggered by explosive volcanism?. Sci. Rep. 5(1), 1–9 (2015).

Article Google Scholar

Cooper, C. et al. Evaluating the relationship between climate change and volcanism. Earth Sci. Rev. 177, 238–247 (2018).

Article ADS Google Scholar

Keller, G. Impacts, volcanism and mass extinction: Random coincidence or cause and effect?. Aust. J. Earth Sci. 52(4–5), 725–757 (2005).

Article ADS Google Scholar

Percival, L. M. et al. Mercury evidence for pulsed volcanism during the end-Triassic mass extinction. PNAS 114(30), 7929–7934 (2017).

Article ADS CAS PubMed PubMed Central Google Scholar

Sigurdsson, H., Cashdollar, S. & Sparks, S. R. J. The eruption of vesuvius in A.D. 79: Reconstruction from Historical and Volcanological evidence. Am. J. Archaeol. 86(1), 39–51 (1982).

Article Google Scholar

Stothers, R. B. Volcanic dry fogs, climate cooling, and plague pandemics in Europe and the Middle East. Clim. Chang. 42(4), 713–723 (1999).

Article Google Scholar

Riede, F. Changes in mid-and far-field human landscape use following the Laacher See eruption (c. 13,000 BP). Quat. Int. 394, 37–50 (2016).

Article Google Scholar

Druitt, T. H., Mellors, R. A., Pyle, D. M. & Sparks, R. S. J. Explosive volcanism on Santorini. Greece. Geol. Mag. 126(2), 95–126 (1989).

Article ADS CAS Google Scholar

Vogel, S. et al. From a stratigraphic sequence to a landscape evolution model: Late Pleistocene and Holocene volcanism, soil formation and land use in the shade of Mount Vesuvius (Italy). Quat. Int. 394, 155–179 (2016).

Article Google Scholar

Silleni, A., Giordano, G., Isaia, R. & Ort, M. H. The magnitude of the 39.8 ka Campanian Ignimbrite eruption, Italy: Method, uncertainties and errors. Front. Earth Sci. 8, 66 (2020).

Article Google Scholar

McCormick, M. P., Thomason, L. W. & Trepte, C. R. Atmospheric effects of the Mt. Pinatubo eruption. Nature 75, 399–404 (1995).

Article ADS Google Scholar

Robock, A. et al. Did the Toba volcanic eruption of ∼74 ka BP produce widespread glaciation?. J. Geophys. Res. Atmos. https://doi.org/10.1029/2008JD011652 (2009).

Article Google Scholar

Gisbert, G., Gimeno, D. & Fernandez-Turiel, J. L. Eruptive mechanisms of the Puig De La Garrinada volcano (Olot, Garrotxa volcanic field, Northeastern Spain): A methodological study based on proximal pyroclastic deposits. J. Volcanol. Geotherm. Res. 180(2–4), 259–276 (2009).

Article ADS CAS Google Scholar

Martí, J., Planagumà, L., Geyer, A., Canal, E. & Pedrazzi, D. Complex interaction between Strombolian and phreatomagmatic eruptions in the Quaternary monogenetic volcanism of the Catalan Volcanic Zone (NE of Spain). J. Volcanol. Geotherm. Res. 201(1–4), 178–193 (2011).

Article ADS Google Scholar

Cimarelli, C., Di Traglia, F., De Rita, D. & Torrente, D. G. Space–time evolution of monogenetic volcanism in the mafic Garrotxa Volcanic Field (NE Iberian Peninsula). Bull. Volcanol. 75(11), 1–18 (2013).

Article Google Scholar

Pedrazzi, D., Bolós, X., Barde-Cabusson, S. & Martí, J. Reconstructing the eruptive history of a monogenetic volcano through a combination of fieldwork and geophysical surveys: the example of Puig d’Àdri (Garrotxa Volcanic Field). J. Geol. Soc. 173(6), 875–888 (2016).

Article ADS CAS Google Scholar

Di Traglia, F. Hydrogeomorphic and sedimentary response to the Late Pleistocene violent Strombolian eruption of the Croscat volcano (Garrotxa Volcanic Field, Spain). Med. Geosci. Rev. 2(2), 217–231 (2020).

Article Google Scholar

Bolos, X., Planagumà, L. & Marti, J. Volcanic stratigraphy of the Quaternary La Garrotxa Volcanic Field (north-east Iberian Peninsula). J. Quat. Sci. 29(6), 547–560 (2014).

Article Google Scholar

Bolós, X. et al. Geophysical exploration on the subsurface geology of La Garrotxa monogenetic volcanic field (NE Iberian Peninsula). Int. J. Earth Sci. 103(8), 2255–2269 (2014).

Article Google Scholar

Bolós, X. et al. Volcano-structural analysis of La Garrotxa Volcanic Field (NE Iberia): Implications for the plumbing system. Tectonophysics 642, 58–70 (2015).

Article ADS Google Scholar

Riede, F. Doing palaeo-social volcanology: Developing a framework for systematically investigating the impacts of past volcanic eruptions on human societies using archaeological datasets. Quat. Int. 499, 266–277 (2019).

Article Google Scholar

Torrence, R. Social resilience and long-term adaptation to volcanic disasters: The archaeology of continuity and innovation in the Willaumez Peninsula, Papua New Guinea. Quat. Int. 394, 6–16 (2016).

Article Google Scholar

Halstead, P. Bad Year Economics: Cultural Responses to Risk and Uncertainty (Cambridge University Press, 1989).

Book Google Scholar

Bradtmöller, M., Pastoors, A., Weninger, B. & Weniger, G.-C. The repeated replacement model e rapid climate change and population dynamics in Upper Pleistocene Europe. Quat. Int. 247, 38–49 (2012).

Article Google Scholar

Solich, M. & Bradtmöller, M. Socioeconomic complexity and the resilience of hunter–gatherer societies. Quat. Int. 446, 109–127 (2017).

Article Google Scholar

Riede, F., Oetelaar, G.A. & Van der Hoek, R. From crisis to collapse in hunter–gatherer societies. A comparative investigation of the cultural impacts of three large volcanic eruptions on past hunter–gatherers. In Crisis to Collapse. The Archaeology of Social Breakdown, edited by T. Cunningham and J. Driessen, 23–39. AEGIS 11 (UCL Presses Universitaires de Louvain, 2017).

Puiguriguer, M. et al. 14C dating of the last Croscat volcano eruption (Garrotxa Region, NE Iberian Peninsula). Geol. Acta 10(1), 1–5 (2012).

Google Scholar

Martí, J., Mitjavila, J., Roca, E. & Aparicio, A. Cenozoic magmatism of the Valencia trough (western Mediterranean): Relationship between structural evolution and volcanism. Tectonophysics 203(1–4), 145–165 (1992).

Article ADS Google Scholar

Martí, J. & Bolós, X. The Neogene-Quaternary Alkaline Volcanism of Iberia. In The Geology of Iberia: A Geodynamic Approach 167–182 (Springer, 2019).

Bond, G. et al. A pervasive millennial-scale cycle in North Atlantic Holocene and Glacial Climates. Science 278, 1257–1266 (1997).

Article ADS CAS Google Scholar

Morales, J. I. et al. Paleogeografía humana durante el Tardiglaciar y Holoceno inicial en el ámbito mediterráneo del NE Ibérico. Cuaternario y Geomorfología 26(3–4), 11–28 (2012).

Google Scholar

González-Sampériz, P. et al. Glacial and Lateglacial vegetation in northeastern Spain: New data and a review. Quat. Int. 140, 4–20 (2005).

Article Google Scholar

González-Sampériz, P. et al. Environmental and climate change in the southern Central Pyrenees since the Last Glacial Maximum: A view from the lake records. CATENA 149, 668–688 (2017).

Article Google Scholar

Gil-Romera, G. et al. Biomass-modulated fire dynamics during the last glacial–interglacial transition at the Central Pyrenees (Spain). Palaeogeogr. Palaeoclimatol. Palaeoecol. 402, 113–124 (2014).

Article Google Scholar

Torres-Rodríguez, E. et al. Last Glacial droughts and fire regimes in the central Mexican highlands. J. Quat. Sci. 30(1), 88–99 (2015).

Article Google Scholar

López-Sáez, J. A. et al. Late Glacial-early holocene vegetation and environmental changes in the western Iberian Central System inferred from a key site: The Navamuño record, Béjar range (Spain). Quat. Sci. Rev. 230, 106167 (2020).

Article Google Scholar

Rasmussen, S. O. et al. A stratigraphic framework for abrupt climàtic changes during the Last Glacial period based on three synchronized Greenland ice-core records: refining and extending the INTIMATE event stratigraphy. Quat. Sci. Rev. 106, 14–28 (2022).

Article ADS Google Scholar

Seierstad, I. K. et al. Consistently dated records from the Greenland GRIP, GISP2 and NGRIP ice cores for the past 104 ka reveal regional millenial-scale d18O gradients with possible Heinrich event imprint. Quat. Sci. Rev. 106, 29–46 (2022).

Article ADS Google Scholar

Revelles, J. et al. Human–environment interaction during the Mesolithic-Neolithic transition in the NE Iberian Peninsula. Vegetation history, climate change and human impact during the Early-Middle Holocene in the Eastern Pre-Pyrenees. Quat. Sci. Rev. 184, 183–200 (2018).

Article ADS Google Scholar

Innes, J. B. & Blackford, J. J. The ecology of late mesolithic woodland disturbances: Model testing with fungal spore assemblage data. J. Arch. Sci. 30, 185–194 (2002).

Article Google Scholar

Dietze, E. et al. Holocene fire activity during low-natural flammability periods reveals scale-dependent cultural human-fire relationships in Europe. Quat. Sci. Rev. 201, 44–56 (2018).

Article ADS Google Scholar

Heidgen, S. et al. Palaeoecological signals for Mesolithic land use in a Central European landscape. J. Quat. Sci. 37(6), 1164–1179 (2022).

Article Google Scholar

Morales-Molino, C. et al. Long-term responses of mediterranean mountain forests to climate change, fire and human activities in the Northern Apennines (Italy). Ecosystems 24, 1361–1377. https://doi.org/10.1007/s10021-020-00587-4 (2021).

Article CAS PubMed Google Scholar

Pérez-Obiol, R. Histoire Tardiglaciaire et Holocène de la végétation de la région volcanique d’Olot (NE PéninsulaIbèrique). Pollen Spores 30(2), 189–202 (1988).

Google Scholar

Burjachs, F. Bauma del Serrat del Pont. Girona, in CARRIÓN, J. (coord), Paleoflora y Paleovegetación de la Península Ibérica e Islas Baleares: Plioceno-Cuaternario. Madrid, MINECO, 306–307 (2012).

Piqué, R., Revelles, J., Burjachs, F., Caruso Fermé, L. & Pérez-Obiol, R. Interdisciplinary approach to the landscape and firewood exploitation during the Holocene at La Garrotxa (Girona, NE Iberia). Quat. Int., 463, 401–413 (2018).

Revelles, J. et al. Mid-holocene vegetation history and Neolithic land-use in the Lake Banyoles area (Girona, Spain). Palaeogeogr. Palaeoclimatol Palaeoecol. 435, 70–85 (2015).

Article Google Scholar

Mehringer, P. J., Blinman, E. & Petersen, K. L. Pollen influx and volcanic ash. Science 198(4314), 257–261 (1977).

Article ADS PubMed Google Scholar

van Geel, B., de Lange, L. & Wiegers, J. Reconstruction and interpretation of the local vegetational succession of a Lateglacial deposit from Usselo (The Netherlands), based on the analysis of micro-and macrofossils. Acta Bot. Neerl. 33(4), 535–546 (1984).

Article Google Scholar

Gross, E. M., Groffier, H., Pestelard, C. & Hussner, A. Ecology and environmental impact of Myriophyllum heterophyllum, an aggressive invader in European Waterways. Diversity 12(4), 127 (2020).

Article Google Scholar

Brunel, S., Schrader, G., Brundu, G. & Fried, G. Emerging invasive alien plants for the Mediterranean Basin. EPPO Bull. 40(2), 219–238 (2010).

Article Google Scholar

Alcalde, G., Saña, M., Mateu, J., Palomo, A. & Buch, M. L. intervenció arqueològica d’urgència en el jaciment a l’aire lliure de la Rodona (Olot, Garrotxa). Tribuna d’Arqueologia 1992–1993, 25–30 (1994).

Google Scholar

Alcalde, G., Saña, M. & Tornero, C. Memòria de les intervencions arqueològiques realitzades al jaciment de la Codella (Les Preses, la Garrotxa). Anys 2003–2006. Servei d’Arqueologia, Direcció General del Patrimoni Cultural, Departament de Cultura, Barcelona (Unpublished report) (2009).

Alcalde, G. et al. Vuit anys de recerca al jaciment arqueològic de la Dou (Sant Esteve d’en Bas, Garrotxa) (2006–2013): del neolític antic al bronze final. Tribuna d’Arqueologia 2013–2014, 196–209 (2017).

Google Scholar

Oms, F. X., Terradas, X., Morell, B. & Gibaja, J. F. Mesolithic-Neolithic transition in the northeast of Iberia: Chronology and socioeconomic dynamics. Quat. Int. 470, 383–397 (2018).

Article Google Scholar

Alcalde, G. & Saña, M. Procés d'ocupació de la Bauma del Serrat del Pont (la Garrotxa) entre 7400 i 5480 cal aC (Publicacions Eventuals d'Arqueologia de la Garrotxa Museu Comarcal de la Garrotxa, 2008).

Saña, M. La gestió dels recursos animals. In: Alcalde, G., Saña, M. (Eds.), Procés d’ocupació de la Bauma del Serrat del Pont (La Garrotxa) entre 7400 i 5480 cal AC 53–59 (Museu Comarcal de la Garrotxa, 2008).

Antolín, F. et al. Changes in the interaction between society and the environment from the mesolithic (10300–8500 cal Bc) to the early neolithic (c. 5400 cal Bc) in can Sadurní cave (Begues, Barcelona province, Spain) a view from the archaeobotanical data. In: Dambon, F. (Ed.), Proceedings of the Fourth International Meeting of Anthracology-Brussels, 8–13 September 2008-. Archaeopress (BAR eInt. Series-2486) 19–29 (Oxford, 2013).

Carbonell, E. & Mora, R. Anàlisi espacial d’un campament prehistòric del Post-Glacial. Sota Palou (Campdevànol). Cypsela revista de prehistòria i protohistòria 66, 35–45 (1985).

Google Scholar

Díez-Canseco, C., Ramírez-Pedraza, I., Carbonell, E. & Tornero, C. Nous treballs de recerca al Torrent d’Estremera: prospecció de l’àrea i excavació de la Bauma dels Fadrins (Queralbs, Ripollès). Actes de les Quinzenes Jornades d’Arqueologia de les Comarques de Girona 66, 43–46 (2020).

Google Scholar

Blaauw, M. & Christen, J. A. Flexible paleoclimate age-depth models using an autoregressive gamma process. Bayesian Anal. 6(3), 457–474 (2011).

Article MathSciNet MATH Google Scholar

Finsinger, W. et al. Fire on ice and frozen trees? Inappropriate radiocarbon dating leads to unrealistic reconstructions. New Phytol. 222(2), 657–662. https://doi.org/10.1111/nph.15354 (2019).

Article PubMed Google Scholar

Schnurrenberger, D., Russell, J. & Kelts, K. Classification of lacustrine sediments based on sedimentary components. J. Paleolimnol. 29(2), 141–154 (2003).

Article Google Scholar

Reading, H. G. Clastic coast. In Sedimentary Environments: Processes, Facies and Stratigraphy, 154–231 (1996)

Weltje, G. J. & Tjallingii, R. Calibration of XRF core scanners for quantitative geochemical logging of sediment cores: Theory and application. Earth Planet. Sci. Lett. 274(3–4), 423–438 (2008).

Article ADS CAS Google Scholar

Weltje, G. J., Bloemsma, M. R., Tjallingii, R., Heslop, D., Röhl, U. & Croudace, I. W. Prediction of geochemical composition from XRF core scanner data: a new multivariate approach including automatic selection of calibration samples and quantification of uncertainties. In Micro-XRF Studies of Sediment Cores 507–534 (Springer, 2015).

Comas Cufí, M. & Thió I Fernández de Henestrosa, S. CoDaPack 2.0: A Stand-Alone, Multi-platform Compositional Software (2011).

Girard, M. & Renault-Miskovsky, J. Nouvelles techniques de préparation en Palynologie appliqués à trois sédiments du Quaternaire final de l’Abri Cornille (Istres-Bouches du Rhône). Bull. AFEQ 4, 275–284 (1969).

Google Scholar

Reille, M. Pollen et spores d'Europe et d'Afrique du nord (URA, CNRS, Laboratoire de Botanique Historique et Palynologie, 1992).

Van Geel, B. A palaeoecological study of Holocene peat bog sections in Germany and The Netherlands. Rev. Palaeobot. Palynol. 25, 1–120 (1978).

Article Google Scholar

Van Geel, B. et al. Environmental reconstruction of a Roman Period settlement site in Uitgeest (The Netherlands), with special reference to coprophilous fungi. J. Archaeol. Sci. 30, 873–883 (2003).

Article Google Scholar

Revelles, J., Burjachs, F. & Van Geel, B. Pollen and non-pollen palynomorphs from the Early Neolithic settlement of La Draga (Girona, Spain). Rev. Palaeobot. Palynol. 225, 1–20 (2016).

Article Google Scholar

Revelles, J. & Van Geel, B. Human impact and ecological changes in lakeshore environments. The contribution of non-pollen palynomorphs in Lake Banyoles (NE Iberia). Rev. Palaeobot. Palynol., 232, 81–97 (2016).

Grimm, E.C. Tilia, Tilia∙Graph and TGView (Illinois State Museum, 1991–2011). http://museum.state.il.us/pub/grimm/tilia/.

Carcaillet, C. et al. Change of fire frequency in the eastern Canadian boreal forests during the Holocene: Does vegetation composition or climate trigger the fire regime?. J. Ecol. 89, 930–946 (2001).

Article Google Scholar

Clark, J. S. Particle motion and theory of charcoal analysis: Source area, transport, deposition, and sampling. Quat. Res. 30, 67–80 (1988).

Article ADS Google Scholar

Finsinger, W., Kelly, R., Fevre, J. & Magyari, E. K. A guide to screening charcoal peaks in macrocharcoal-area records for fire-episode reconstructions. The Holocene 24(8), 1002–1008 (2014).

Article ADS Google Scholar

Higuera, P. E., Gavin, D. G., Bartlein, P. J. & Hallet, D. J. Peak detection in sediment-charcoal records: Impacts of alternative data analysis methods on fire-history interpretations. Int. J. Wildl. Fire 19, 996–1014 (2010).

Article Google Scholar

Cappers, R. T. J., Bekker, R. M. & Jans, J. E. A. Digitale Zaden atlas van Nederland (Barkhuis Publishing & Groningen University Library, 2006).

Google Scholar

Mauquoy, D. & van Geel, B., Mire and peatmacros. In: Elias, S.A. (Ed.), Encyclopedia of Quaternary Science vol. 3 2315–2336 (Elsevier, 2007).

Meisch, C. Freshwater Ostracoda of Western and Central Europe (Spektrum Akademischer Verlag, 2000)

Fuhrmann, R. Atlas quartärer und rezenter Ostrakoden Mitteldeutschlands (Altenburger Naturwissenschaftliche Forschungen, Naturk und liches Museum Mauritianium, 2012).

Soulié-Märsche, I. Etudecomparée de gyrogonites de Charophytes actuelleset fossiles et phylogénie des genresactuels 237 (Imprimerie des Tilleuls, 1989).

Battarbee, R.W. et al. Diatoms. In: Smol, J.P., Birks, H.J.B., Last, W.M. (eds) Tracking Environmental Change Using Lake Sediments, vol 3. Terrestrial, Algal, and Siliceous Indicators pp 155–202 (Kluwer, 2001).

Krammer, K. & Lange-Bertalot, H. Bacillariophyceae. Naviculaceae 2/1 876 (Gustav Fischer Verlag, 1986).

Krammer, K. & Lange-Bertalot, H. Bacillariophyceae. Bacillariaceae, EpithemiaceaeSurirellaceae 2/2 596 (Gustav Fischer Verlag, 1988).

Krammer, K. & Lange-Bertalot, H. Bacillariophyceae. Achnanthaceae 2/4 473. (Gustav Fischer Verlag, 1991).

Krammer K. & Lange-Bertalot, H. Bacillariophyceae. Centrales, Fragilariaceae, Eunotiaceae 2/3 576 (Gustav Fischer Verlag, 1991).

Rivera-Rondón, C. A. & Catalan, J. Diatom diversity in the lakes of the Pyrenees: An iconographic reference. Limnetica 36, 127–396 (2017).

Google Scholar

Ryves, D. B. et al. Physical and chemical predictors of diatom dissolution in freshwater and saline lake sediments in North America and West Greenland. Limnol. Oceanogr. 51(3), 1355–1368 (2006).

Article ADS CAS Google Scholar

Download references

In memoriam of Dr. Gabriel Alcalde, who contributed significantly to the conception of this research based of his broad expertise in archaeological research in the region of La Garrotxa. J.R. developed this research with a contract Juan de la Cierva Incorporación (IJC2020) (MCINN, Spain), in the research group GAPS (2017 SGR 836). The Institut Català de Paleoecologia Humana i Evolució Social (IPHES-CERCA) has been funded by MCINN in the programme of Units of Excellence ‘Maria de Maeztu’ (CEX2019-000945-M). This research has been undertaken in the frame of the projects ‘2014/100638-Evolució del poblament i ús del territori al Prepirineu Oriental durant la prehistòria recent (8000–900 cal ANE): Anàlisi arqueoecològica de les dinàmiques de canvi social i de la gestió dels recursos naturals (2014-2017)’and‘CLT009/18/00023-Evolució del poblament i ús del territori al Prepirineu oriental durant la prehistòria recent (10000-900 CAL ANE): Anàlisi arqueoecològica de les dinàmiques de canvi social (2018-2021)’.

Gabriel Alcalde is deceased.

Institut Català de Paleoecologia Humana i Evolució Social (IPHES-CERCA), Zona Educacional 4, Campus Sescelades URV (Edifici W3), 43007, Tarragona, Spain

Jordi Revelles & Francesc Burjachs

Universitat Rovira i Virgili (URV), Àrea de Prehistòria, Avinguda de Catalunya 35, 43002, Tarragona, Spain

Jordi Revelles & Francesc Burjachs

Department of Geosciences, Institute of Environmental Assessment and Water Research (IDAEA-CSIC), Jordi Girona 18-26, 08034, Barcelona, Spain

Joan Martí Molist

ICREA, Pg. Lluís Companys 23, 08010, Barcelona, Spain

Francesc Burjachs

ISEM, University of Montpellier, CNRS, IRD, EPHE, Montpellier, France

Walter Finsinger

Laboratorio de Evolución Humana/IsoTOPIK, Departamento de Historia, Geografía y Comunicación, Universidad de Burgos, Plaza Misael Bañuelos s/n, Edificio de I+d+i, 09001, Burgos, Spain

Eneko Iriarte

"Cavanilles" Institute of Biodiversity and Evolutionary Biology, University of Valencia, Catedrático José Beltrán Martínez, 2, 46980, Paterna, Spain

Francesc Mesquita-Joanes & Maria A. Rodrigo

Unitat d’Ecologia, Departament de Biologia Animal, de Biologia vegetal i Ecologia, Universitat Autònoma de Barcelona, 08193, Bellaterra, Catalonia, Spain

Sergi Pla-Rabés

CREAF, Center for Ecological and Forestry Applications, 08193, Cerdanyola del Vallès, Catalonia, Spain

Sergi Pla-Rabés

Tosca, Environment Services of Education, Casal dels Volcans, Av. Santa Coloma, 17800, Olot, Spain

Llorenç Planagumà

Departament Història i Històriadel’Art, Universitat de Girona, 17071, Girona, Spain

Gabriel Alcalde

Departament de Prehistòria Edifici B, Facultat de Filosofia i Lletres, Universitat Autònoma de Barcelona, 08193, Barcelona, Spain

Maria Saña

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

J.R., M.S. and G.A. conceived the research; J.R. analysed pollen and sedimentary charcoal; E.I. analysed sedimentology and geochemistry; J.M.M. characterised tephra; F.M.J. analysed ostracods; M.A.R. analysed charophyte gyrogonites; S.P. analysed diatoms; W.F. contributed in numerical analyses for fire peak detection and age-depth modelling; all authors contributed to the final interpretation and writing of the manuscript with major contributions by J.R., J.M.M and E.I.

Correspondence to Jordi Revelles.

The authors declare no competing interests.

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

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

Reprints and Permissions

Revelles, J., Martí Molist, J., Burjachs, F. et al. Socio-ecological impact of monogenetic volcanism in the La Garrotxa Volcanic Field (NE Iberia). Sci Rep 13, 8168 (2023). https://doi.org/10.1038/s41598-023-35072-0

Download citation

Received: 07 December 2022

Accepted: 12 May 2023

Published: 20 May 2023

DOI: https://doi.org/10.1038/s41598-023-35072-0

Anyone you share the following link with will be able to read this content:

Sorry, a shareable link is not currently available for this article.

Provided by the Springer Nature SharedIt content-sharing initiative

By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.