Next Article in Journal / Special Issue
Land-Cover Dependent Relationships between Fire and Soil Moisture
Previous Article in Journal
Emergence and Evolution of Anthropogenic Landscapes in the Western Mediterranean and Adjacent Atlantic Regions
Previous Article in Special Issue
Half-Century Changes in LULC and Fire in Two Iberian Inner Mountain Areas
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

The Iberian Peninsula’s Burning Heart—Long-Term Fire History in the Toledo Mountains (Central Spain)

1
Instituto de Historia, Consejo Superior de Investigaciones Científicas, C/Albasanz 26-28, 28037 Madrid, Spain
2
Departamento de Geografía, Universidad Autónoma de Madrid, C/ Francisco Tomás y Valiente 1, 28049 Madrid, Spain
3
Département de Géographie, Université de Montréal, Pavillon 520, Chemin Côte Sainte-Catherine, C. P. 6128, Succursale Centre-ville, Montréal, QC H3C 3J7, Canada
4
Departamento de Geografía, Urbanismo y Ordenación del Territorio, Universidad de Cantabria, Avenida Los Castros 44, 39005 Santander, Spain
5
Institute of Plant Sciences and Oeschger Centre for Climate Change Research, University of Bern, Altenbergrain 21, 3013 Bern, Switzerland
*
Author to whom correspondence should be addressed.
Submission received: 17 June 2019 / Revised: 3 October 2019 / Accepted: 12 October 2019 / Published: 16 October 2019
(This article belongs to the Special Issue Land-Use and Fire around the World from the Past to the Present)

Abstract

:
Long-term fire ecology can help to better understand the major role played by fire in driving vegetation composition and structure over decadal to millennial timescales, along with climate change and human agency, especially in fire-prone areas such as the Mediterranean basin. Investigating past ecosystem dynamics in response to changing fire activity, climate, and land use, and how these landscape drivers interact in the long-term is needed for efficient nature management, protection, and restoration. The Toledo Mountains of central Spain are a mid-elevation mountain complex with scarce current anthropic intervention located on the westernmost edge of the Mediterranean basin. These features provide a perfect setting to study patterns of late Holocene fire activity and landscape transformation. Here, we have combined macroscopic charcoal analysis with palynological data in three peat sequences (El Perro, Brezoso, and Viñuelas mires) to reconstruct fire regimes during recent millennia and their linkages to changes in vegetation, land use, and climatic conditions. During a first phase (5000–3000 cal. BP) characterized by mixed oak woodlands and low anthropogenic impact, climate exerted an evident influence over fire regimes. Later, the data show two phases of increasing human influence dated at 3000–500 cal. BP and 500 cal. BP–present, which translated into significant changes in fire regimes increasingly driven by human activity. These results contribute to prove how fire regimes have changed along with human societies, being more related to land use and less dependent on climatic cycles.

1. Introduction

Fire has been perceived as a dramatic hazard despite its early and deep presence in the Mediterranean basin. In fact, it has played an essential role in defining the landscapes and biota present there, so much so that fire became an intrinsic component of the Mediterranean cultures [1,2,3,4]. On top of that, fire has long been used in the Mediterranean region as a major tool for landscape transformation and management by human communities [5,6,7,8]. Indeed, a great number of human activities are linked to the use of fire in the Mediterranean basin, which has led to extremely short fire return intervals, usually below 100 years and sometimes even shorter [9,10]. Therefore, it is clear that fire is part of Mediterranean ecosystems. This makes it impossible to understand their evolution and current dynamics without looking at the role played by fire through time [3,11,12,13,14,15,16].
Fire has driven the distribution, structure, and composition of Mediterranean ecosystems, building time-since-fire mosaic patterns and increasing the biodiversity in some areas [6,10,17,18]. It also plays an important role in determining communities’ susceptibility to many other hazards, such as floods or resprouting ability [19,20].
However, fire regimes do not show a strict and uniform pattern in time and space [21,22,23]. Fire occurrence is determined by many interlinked factors, such as vegetation dynamics, climatic conditions, and land use, among others. Vegetation composition and structure regulate the availability, flammability, and amount of fuel present in every landscape [10,13,24]. At the same time, it is necessary to point out those climatic drivers such as temperature and precipitation, regulating the amount of biomass and fuel available in the ecosystem and humidity or dryness controlling fuel ignition [13,21]. Furthermore, fire can start by natural causes or mediated by humans, which also plays an important role in the fire regime. Despite the importance of this element in human management, some activities such as grazing or different land use (e.g., agricultural and urban patterns) modify the occurrence of fires [10,13,24,25]. Therefore, there is a need to understand fire history in the Mediterranean area and the relationship between all the actors involved in the fire regime. The current fire dynamic is a consequence of the past fire regimes but despite the adaptations, fire has become a threat for ecosystems and human communities because of the changes in the fire regime drivers [3,12,26]. This circumstance brings uncertainty to the future because some new changes are expected and this makes it difficult to choose the best management goals which could take into account the many involved factors and helped to preserve ecological richness and land use interests, for example [3,7,26,27,28].
Different methods to reconstruct past fire regimes have been developed, such as fire scars and historical data, but these rarely allow researchers to go back longer than several centuries [15,29]. Here is where sedimentary charcoal records emerge as the best tool for reconstructing fire history, given that multi-millennial records are relatively frequent [30,31]. Fire occurrence can be inferred from the charcoal deposited in lake sediments and peatlands, among other natural deposits. Charcoal accumulation rates are actually a proxy for the biomass burnt in an area during one fire episode, which may in turn embrace one or several fire events depending on the time resolution of sediment record and the subsequent sampling [28,32]. Both theoretical and empirical research on charcoal dispersal has shown that peaks in the deposition of macroscopic charcoal particles are commonly associated with local fire events mostly occurring in the study catchment and adjacent areas; meanwhile sieved charcoal is more accurate at showing local fires [33,34,35,36,37,38]. Therefore, the so-called “peak detection analysis” is often used to reconstruct local fire regimes [10,14,26,33].
In this paper, we use macroscopic charcoal series from three mires located in different areas of the Toledo Mountains, in the heart of the Iberian Peninsula, to reconstruct historical fire regimes in the region. The diverse locations and historical development could highlight unequal fire patterns related to the physical and socio-economical features of these lands. We have focused on the possible occurrence of common trends and their causes as well as the divergences between the reconstructions. By comparing the fire history reconstructions with vegetation and land use records from the study sites, we have assessed whether climate change or human activities took the lead in driving fire activity in these mid-elevation mountains of central Spain and the impact of fire on long-term vegetation dynamics. Furthermore, we have compared the regional fire activity of the Toledo Mountains with the general trends in Spain. Finally, we provide some management guidelines inspired by the inferences taken from the fossil charcoal records.

2. Materials and Methods

2.1. Study Area

The Toledo Mountains is a west-to-east-aligned mid-mountain range complex located in central Spain, between Toledo and Ciudad Real provinces (Figure 1). This mountain chain extends over more than 3800 km2 at the western side of the southern Iberian plateau, where it separates the River Tagus and River Guadiana’s basins. Eastwards, the Toledo Mountains border La Mancha plain [39,40]. The average height of this mid-mountain complex is 800–1100 m.a.s.l. [39,41]. One of the most typical geomorphological features of this mid-mountain area is the locally named rañas, ramp alluvial/colluvial deposits located at the foothills of the main massifs [39].
The climate of the Toledo Mountains is predominantly Mediterranean, with dry and warm summers, cold and wet winters, and the usual occurrence of extreme events such as severe droughts and floods. The western side of the range receives certain oceanic influence related to the westerlies (lower seasonality, higher rainfall). The average annual temperature is 17 °C and the mean annual rainfall oscillates between 600 and 800 mm [40,41,44]. The vegetation in the Toledo Mountains meso-Mediterranean foothills is mainly composed of holm oak (Quercus ilex subsp. ballota) and cork oak (Q. suber) woodlands. Holm oak-dominated communities usually include meso-thermophilous taxa such as Arbutus unedo, Phillyrea angustifolia, Pistacia terebinthus, or Pyrus bourgaeana. The Quercus suber woodlands are also often associated with deciduous trees such as Quercus faginea subsp. broteroi, Q. pyrenaica and Acer monspessulanum, on moister settings. Deciduous Quercus pyrenaica forests are dominant in the supra-Mediterranean belt (roughly above 1000 m a.s.l.). There also are some pine reforestations in the Toledo Mountains, carried out during the second half of the 20th century. The Toledo Mountains host some of the southernmost populations of several temperate and Atlantic species such as Myrica gale, Corylus avellana, Betula pendula, or B. pubescens [42,43,45,46,47].
Previous research has identified more than 130 peatlands in the Toledo Mountains, most of them located in the Ciudad Real province. Despite this optimistic number, these peatlands present a limited extension (the smallest mire presents a surface of 0.02 ha and the largest 113.6 ha, but the average extension is around 4.5 ha) [40]. The Toledo Mountains peatlands are mainly minerotrophic mires formed by terrestrialization and connected to groundwater discharge in springs or to impeded stream flow in the rañas [40].
For this study, three mires have been selected (Table 1, Figure 1) based on their different situations into the mountainous complex. They are in the Ciudad Real province, on the northern and southern sectors of the Toledo Mountains, helping us to discern whether fire history significantly differed in these two areas during recent millennia. In addition, Viñuelas and Brezoso are in higher and more difficult-to-access locations, meanwhile El Perro mire is close to the valley, nearer to the populated area. Some of the mires present in these provinces are protected under the Natura 2000 Network and by specific strategies created inside the Spanish and Provincial laws, such as micro-reserves, river reserves, or the most important figure: national parks. Brezoso and Viñuelas mires are under this particular protection plan [40]. Despite efforts the different authorities have made to guarantee the protection of these special ecosystems, they have to face several threats such as drainage (Brezoso and El Perro), erosion (Viñuelas, Brezoso, and El Perro), and trampling/overgrazing by wild ungulates (Viñuelas, Brezoso), among others. Just to illustrate the relevance of these degradation processes, intensive drainage to reclaim land for agricultural purposes has led peat accumulation to stop at El Perro mire [40,42,43].

2.2. Sampling and Chronology

A Russian peat sampler (5 cm diameter) was used to core El Perro mire. In the field, peat sections were placed in PVC tubes, and later, in the laboratory, they were stored under dark and cold conditions (4 °C). The core was then subsampled into contiguous 1 cm thick slices. The chronology of the El Perro peat sequence was established by means of 11 AMS (accelerator mass spectrometry) 14C dates on peat samples. AMS dating was conducted at the Poznán (Poznán, Poland) and Belfast (United Kingdom) radiocarbon laboratories. The radiocarbon dates were converted to calendar years using the software CALIB 7.1 with the IntCal13 calibration curve [48]. Post-bomb samples (i.e., younger than 1950 CE = 0 cal. BP) were treated according to Hua and Barbetti [49]. Finally, an age–depth model was produced fitting a smoothing spline (smoothing parameter=0.2) to the accepted radiocarbon dates using the R-Package “clam” [50,51] running in R 3.2.2. (Table 2, Figure 2). Confidence intervals of the calibrations and the age–depth model were calculated at 95% (2σ) with 1000 iterations (Table 2). The chronology of the Brezoso and Viñuelas paleo-ecological records is also based on AMS radiocarbon dating of terrestrial plant macrofossils and peat (the reader is referred to [42,43] for details), and the age–depth relationships were also modelled using smoothing splines fitted with clam (details in [42,43] for Brezoso and in [43] for Viñuelas).

2.3. Pollen Analysis

Pollen analysis was carried out along El Perro, Brezoso, and Viñuelas following the methods described by Moore [52] treated with hydrochloric acid (HCl), potassium hydroxide (KOH), and hydrofluoric acid (HF) for removing carbonates, silicates, and organic matter, respectively. The El Perro pollen samples were recovered using Thoulet solution [53]; meanwhile Brezoso and Viñuelas were sieved using a 250 µm mesh and decanted. Lycopodium tablets were added to the samples to estimate the pollen concentration [42,43,54].
Palynomorphs were identified using different pollen atlases and keys [52,55,56,57]. Betula, Myrica and Corylus pollen types were identified following [46,58], and Erica pollen types according to [59]. Pollen diagrams have been plotted against depth and age using TGview [60].

2.4. Charcoal Analysis

For macroscopic charcoal analysis, we treated contiguous 1 cm thick subsamples of 1 cm3 volume throughout the El Perro sequence [33,61], soaking them in a 10% KOH solution for 24 h and then adding 15% H2O2 for 24 h to bleach the uncharred organic matter. The sediment was then sieved through a 125 µm mesh, and we finally counted the number of charcoal particles in each sample under the stereomicroscope. We also included in this study the previously published macroscopic charcoal series from Brezoso and Viñuelas [43], which also consist of contiguous 1 cm thick 1 cm3 samples, treated in a different laboratory following a slightly different but comparable protocol (100 µm instead of 125 µm mesh size and sodium hexametaphosphate as deflocculant agent instead of KOH; see details in [43]). In any case, those differences do not affect either the statistical treatment of the data nor the obtained results, so they are completely comparable.
Charcoal analysis has been run with MATLAB and R, using CharAnalysis software [26,30,62] and Paleofire R-Package [63]. These analyses are based on the charcoal accumulation rates (CHAR) curve decomposition by two main elements: CHARbackground (Cback) and CHARpeak (Cpeak) [27]. The first one is the slow trend variation charcoal income, related with the regional fire event charcoal production and adds some noise to the charcoal records; meanwhile, Cpeak, superimposed data over the Cback, is related to the local fire events [22,64,65].
The analysis requires a pre-treatment step which interpolates to the median resolution CHAR data to correct the temporal variability and the bias produced by the different sedimentation rates which could shadow many of the smallest charcoal accumulation data. Once the CHAR has been interpolated, it is necessary to apply some smooth trends for detrending CHAR, obtaining Cback and Cpeak. In this case, all the 5 possible smoothing methods available within CharAnalysis have been applied for finding the Cback [31] using a smoothing window width ranging from 100 to 1500 years such as t=100, 175, 200, …, 1500. In total, for each site we obtained a 5 × 55 = 275 reconstructions ensemble member. Cpeak has been obtained by subtracting Cback to interpolated CHAR (Ci-Cback) [26].
Cpeak curve is still composed by the Cnoise, which is the statistical noise derived from charcoal accumulation, mixing, and redeposition, and Cfire, related to the occurrence of one or more fires in the site’s catchment area, hereafter referred as fire events. They have been obtained and separated by using a Gaussian mixture model, following the recommendations made by Gavin et al. [27], using a 99th-percentile threshold in each reconstruction for separate both components. The local threshold has been applied using a 500-year time window [26]. We identified fire events corresponding to one or more fires in the mire catchment area. Given the short fire return intervals in the Mediterranean area, the time resolution of paleo-ecological analysis is too coarse to identify, without uncertainty, single fire events. Then fire peak may correspond to one or multiple fire events. A signal-to-noise index (SNI) has been used to select all members (i.e., analysis reiterations) with a median SNI > 3 [66]. We used 75% agreement between the members for a fire episode that is considered to be a statistically robust fire event [10].
Fire frequency was estimated using kernel density estimation with a width of 250 years. Bootstrapped 95% confidence intervals around the fire frequency were established by a bootstrap resampling of the kernel density estimation using the Paleofire R-Package [50,67]. Those data help to better understand the changes in the fire regime trends [10,63]. Fire return interval was calculated using the lagged differences function in R using the robust fire events found in previous steps in the analysis.

2.5. Fire Synchrony

We used the Ripley K-function simplified for one dimension to test whether fire events in the three sites co-occur more than would be expected by chance. To do this, we calculated the bivariate K-function and its transform (L) using an edge correction; the function was computed for T/2 years where T is the length in years of the longest record (here T/2 was 2000 years). Then, we assessed 95% confidence envelope around the L-function by randomizing the fire events of the three sites 200 times. The analysis was run using the K1D software v1.2 [64]. The L-function lying above or below the confidence envelope indicated significant synchrony or asynchrony. When the L-function was between the envelope it indicated an independence of fire occurrence for the three sites. For details on the Ripley K-function the reader may refer to [64,68,69].

2.6. Sites Comparison

Compositing of multiple charcoal records from different settings, with different methods and trends, has been recognized for almost two decades as a way to understand trends in biomass burning and global past fire emissions [33,70]. Site comparison has been made using R and the pf.CompositeLF function included in the Paleofire R-Package [63,70]. As for the CHAR analysis it is necessary to interpolate the data before any comparison to reduce the bias due to the different sedimentation rates and the temporal resolution anomalies present in the studied cores due to their different natures. The data produced in the Toledo Mountains have been compared to those 33 available sites for the Iberian Peninsula in the Global Charcoal Database. We used the same procedure as for comparing the three selected Toledo Mountains sites (see above). At the same time, we conducted a Kruskal–Wallis test to understand the statistical differences between the compared composites (Table S1) [13,63].

3. Results and Discussion

The results show different charcoal accumulation trends for the three studied mires. Mean CHAR are variable in the studied sequences (Figure 3), with the highest mean values at Viñuelas (26.22 # cm−2 year−1), followed by Brezoso (8.64 # cm−2 year−1), and finally the lowest accumulation rate at El Perro (4.1 # cm−2 year−1). In all the records, it is possible to see a major fire activity in the oldest samples, a second period between 3000 and 500 cal. BP with variable results, and one last significant change in the last five hundred years in the three records. (Figure 3, Figure A1 and Figure A3 in Appendix A).

3.1. Fire along the Time: First Sights

The Brezoso mire charcoal record shows an overall high fire activity, with very high CHAR values at the bottom of the sequence and comparatively moderate to low values in the rest of the samples (Figure 3 and Figure A2 in Appendix A). Altogether, 17 robust fire episodes have been identified along the sequence. On the other hand, CHAR values are highest in the middle section of the Viñuelas sequence (Figure 3). Apart from a minor peak around 3900 cal. BP, CHAR are mostly low (3–4 # cm−2 year−1) until ca. 3200 cal. BP, when an increasing trend begins. The Kruskal–Wallis test corroborates this perception, separating three major phases regarding the values of CHAR (Figure A3 in Appendix A). In this case, 7 fire episodes have been detected along the record (Figure 3, Figure 4 and Figure 5). Finally, El Perro CHAR and fire frequency curves follow an overall gently ascending trend, especially evident from 1000–500 cal. BP. It is particularly remarkable the low amount of biomass burning registered at the bottom of the record, with CHAR values ranging from 0.7 to 2.12 # cm−2 year−1 between ca. 5000 to 3750 cal. BP, from when CHAR begins to increase. The second part of the record, until ca. 400 cal. BP, shows rather continuous moderate CHAR values indicative of enhanced biomass burning. Subsequently, a sharp increase in CHAR starts at ca. 400 cal. BP, which leads to the highest CHAR values of the whole sequence (28.28 # cm−2 year−1 in the peak at ca. 90 cal. BP, Figure 3).
There are noticeable differences in the mean fire return intervals shown by the mires: 244 years in Brezoso, 296 years in Viñuelas, and 548 years in El Perro (Figure A1 in Appendix A). However, despite the large differences in fire return intervals and fire frequencies (Figure A1 in Appendix A), it is possible to discern some patterns, as shown in the L-function (transform of the K-function) [64]. It has been run for the fire events represented in (a) with 95% confidence envelope based on 1000 randomizations and it shows a fire recurrence within 150–250 years, which points through a growing human influence in the fire regime in the Toledo Mountains since shorter fire intervals are related with human land management more than climatic trends, something visible in the pollen data [13] (Figure 4).
Previous studies have widely proved the importance of fire in defining the structure, composition, and distribution of Mediterranean ecosystems [24,71], and this holds particularly true for the Iberian Peninsula. During the early and middle Holocene, a drier climate favored the occurrence of particularly frequent and severe fire events modifying the entire landscape [4,9,12], whereas late Holocene fire history is mainly human-mediated. Fire has been used as a tool to manage the environment, increasing or improving the available resources, creating new spaces for crops or livestock [10,72]. This use of fire as a management tool is especially relevant in mountainous lands, although some authors have questioned the importance of human influence in shaping mountain landscapes during the middle Holocene [73]. In fact, many studies have demonstrated the essential role that human societies have had in these spaces, by using the resources, clearing the vegetation for crop cultivation and livestock, raising even to the point of developing an unquestionable co-evolution [10,74,75].
Despite the above-mentioned importance of this element in the Iberian Peninsula, fire history in the Toledo Mountains does not match the general trends reconstructed from the Iberian Peninsula fire records, as is evident in Figure 5, nor shows a regular pattern in the whole territory. As Figure 5 reflects, there is a lack of resemblance between general trends of the Iberian charcoal records and the selected sequences from the Toledo Mountains. Although through the core it is possible to find some minimum correlations as it highlights the Kruskal–Wallis tests, it is easy to see how different the trends are. Nevertheless, it is possible to find agreement between the periods centered in 1250 cal. BP, where the Kruskal–Wallis tests shows statistical significance between the Iberian Peninsula and the Toledo Mountains records. Apart from that, it is possible to find partial matches between the trends in the Toledo Mountains and the Iberian Peninsula, i.e., in the period centered at 3500 cal. BP and 2500 cal. BP.
In the Toledo Mountains records, there is certain resemblance in the timing of the reconstructed changes in fire regimes. For instance, between ca. 4000–3000 cal. BP, fire occurrence tends to increase in the Iberian Peninsula as a whole while it seems to diminish in the Toledo Mountains. This disagreement is easily visible in Figure 5 and statistically different, as the Kruskal–Wallis tests highlighted, finding only partial agreement in the period centered in 3500 cal. BP. There is more disagreement in the oldest moments of the record, as is possible to see in Figure 5. Until around 2500 cal. BP the fire activity recorded from the Iberian Peninsula and the Toledo Mountains is statistically different one from the other, except for the beginning of the core, when some significance is registered. From around 2000 cal. BP onwards, it is possible to find more agreement between the composites, although only there is a perfect match between the fire records in the above-mentioned events around 1500–1250 cal. BP. In the last part of the record, the agreement is once again scarce, although some similarities can be recognized in the fire activity in the comparison, as the Kruskal–Wallis tests show (Figure 5).

3.2. Climatic Influence: from 5000–3000 cal. BP

The first part of the El Perro charcoal record, from 5000 to 3000 cal. BP, shows an overall, although soft, rise in fire activity with increasing CHAR values (Figure 3). Additionally, the fire regime reconstructions depict increasing fire frequencies in all the sites for this period (Figure A1 in Appendix A). Dry pulses during the 4.2 ka event (ca. 4300–3800 cal. BP) [76,77,78,79,80] provide a favorable setting for fire spread as suggested by the Brezoso and, to a lesser extent, the Viñuelas charcoal records [76,77,78]. Indeed, the largest CHAR values of the last 4000 years are recorded at Brezoso from 3900 to 3800 cal. BP (53.29–230.06 # cm−2 year−1), which is corroborated by the Kruskal–Wallis test (Figure A3 in Appendix A). Peak detection suggests the occurrence of two fire episodes with a severe impact on the surrounding vegetation [42,43]. Thus, Quercus woodlands retreated while disturbance-tolerant communities such as Betula woods and Mediterranean shrublands dominated by Erica expanded [42,43,81]. Minor increases in Asteraceae, Poaceae, and anthropic herbs (see Figure 6 and [42]) point to a discrete but constant presence of human activity in the area. Similarly, discontinuous occurrences of cereal pollen grains suggest that cereal cultivation occurred not far from the mire but was not widespread (Figure 6).
El Perro mire registers the lowest CHAR values in the initial part of the sequence in comparison to Brezoso and Viñuelas mires, averaging values 1.86 # cm−2 year−1 (Figure 3). Fire frequency does not experience major changes during this phase but a final increase starting at ca. 3500 cal. BP (Figure A1 in Appendix A), given that only one fire episode has been identified during this period (at ca. 3661 cal. BP) (Figure 3 and Figure A1 in Appendix A). Following this fire episode, deciduous and evergreen Quercus communities were partly replaced with Mediterranean shrublands and grasslands (Figure 6). Minor increases in grazing pollen indicators and coprophilous fungi suggest that probably fire events in this period were human-induced, in particular by Bronze Age communities (ca. 4300–3800 cal. BP) [78,81,83].
CHAR values at Viñuelas mire (Figure 3) are considerably low at the beginning of the record except for the small peak at the very bottom of the sequence. Low fire activity at Viñuelas during this period is somehow unexpected, considering the dominant dry conditions [83,84]. This behavior is marked in the Kruskal–Wallis test as different, and it is coherent with the other results (Figure A3 in Appendix A). This likely fire episode, although not identified as robust by peak detection analysis, probably influenced the local spread of Betula at the expense of Quercus [43]. The palynological record indicates that human activities around this mire were rather limited during this period [43] (Figure 6).
This period roughly corresponds to the Bronze Age (ca. 4150–3350/3300 cal. BP), when human communities developed livestock husbandry and agriculture as their main economic activities and land use becomes more intense, producing real impact in the environment [15,42,85], despite the many changes in human communities’ behavior due to climatic variability [40,83].

3.3. Climatic and Human Influence: 3000–500 cal. BP

In this period, the tendency presented by the three studied sites is mainly different. Nevertheless, there is an evident reflection in the statistical tests which separates this part of the record from the extremes (Figure 3 and Figure A3 in Appendix A). This long time includes several climatic and cultural periods, with potential influence on fire history such as the warm and humid Iberian–Roman Humid Period (2660–1450 cal. BP) [86], the colder and drier Early Medieval Cold Episode (1450–1050 cal. BP) or the warm Medieval Climate Anomaly (1050–650 cal. BP). Finally, it also includes the Little Ice Age (650–150 cal. BP), a long time when the cold and dry conditions characterized the climate in the northern hemisphere [87,88].
The Brezoso mire shows the lowest CHAR values during this phase, especially at 2500–1500 cal. BP (Figure 3). Nevertheless, fire frequency is still higher than in Viñuelas and El Perro despite it following a decreasing trend throughout this period (Figure A1, supplementary material). In contrast to Brezoso and El Perro mires, the Viñuelas charcoal record depicts significant biomass burning throughout this period, following a rising trend starting at ca. 3200 cal. BP (Figure 3) [43]. It is interesting to highlight the deceleration in fire activity between ca. 1150 cal. BP and ca. 850 cal. BP, a noteworthy reduction which is also recorded in El Perro mire (Figure 3, and Figure A1 and Figure A2 in Appendix A). In both mires between 3000–2000 cal. BP there is still a remarkable amount of charcoal accumulation, although in El Perro mire it begins to increase a bit earlier (Figure 3). Interestingly, the ensemble member procedure for peak detection analysis detects a nearly synchronous fire episode dated at ca. 2651 cal. BP at the three study sites (Figure 3 and Figure 4), at the beginning of the Iron Age, when an intensification of human activities is usually recorded [9,42,89].
This fire episode or episodes occurred after the climatic event 2.8 ka. cal. BP. [90,91,92] and within the Iberian–Roman Humid Period, when climatic conditions were mostly warm and humid [86,87]. At Brezoso and El Perro, Betula slightly spread following this fire episode alongside Erica under low anthropogenic impact, just indicated by minor increases in Cardueae, Asteraceae, and Poaceae. Lack of Cerealia-type pollen grains and grazing indicators indicated that farming was limited around the mires [42,43] (Figure 6). In marked contrast, at Viñuelas CHAR values are rather high from 2542–662 cal. BP, with a couple of fire episodes recorded (2839, 2539 cal. BP). Enhanced biomass burning (Figure 3) and increased fire frequency (Figure A1 in Appendix A) seem to have favored certain local spread of Betula around the mire. The paucity of anthropogenic indicators during this period suggests climate–vegetation feedbacks promoting the increase in fire activity.
From 2000 to 500 cal. BP the three mire records show different trends (Figure A1 in Appendix A). Viñuelas mire shows an evident increase in the CHAR values, (mean values for this period 42.28 particles cm−2 year−1 (Figure 3)). The Kruskal–Wallis test individualizes these charcoal trends in the whole record as statistically different (Figure A3 in Appendix A). On the contrary, El Perro and Brezoso do not accumulate a great amount of CHAR (Figure 3). During this sequence, the Roman Period (2000–1500 cal. BP) promoted the intensification of cereal cultivation and the livestock. This has its reflection in the pollen diagrams (Figure 6) with the spread of Poaceae, Rumex, and Plantago among others, indicating the success of these activities, which are kept in the first part of the Early Middle Age [42,43]. The Visigothic Kingdom, which maintained strong persistence in the rural areas [77], successfully developed a slash-and-burn system, preparing the land for crops and grazing by using fire as the main tool for managing the environment [42,93,94,95,96].
Between all these fire events, the moments where there is a lack of fire activity are also relevant (Figure A1, supplementary material) although it is precisely in these kinds of situations where the differences between the three records are more visible. Viñuelas is the mire which most of these periods without fire or little fire activity registers (Figure 3 and Figure A1 in Appendix A). Nevertheless, Viñuelas’ record gathers in this period three great peaks in 1019, 1296, and 1479 cal. BP, with CHAR values between 60 and 80 particles cm−2 year−1. The fire frequency increases in this interval and tends to descend from ca. 1000 cal. BP. El Perro and Brezoso mires, by their part, share a moment of fire absence from ca.1200–900 cal. BP. This period with less fire activity shows how the lack of fire can help the forest to develop and be useful for the human communities living around the study area. Thus, it is possible to appreciate an increase in the woodlands with the spread of Quercus ilex, Quercus pyrenaica, and Quercus suber, which are species with an evident economic interest for the human communities. At the same time, the reduction trends of Corylus and Betula, opportunistic species, are visible in the diagram. In this moment, Brezoso mire gathers the highest percentages of cereal pollen, also present in El Perro mire, and there is an evident increase in livestock and human pressure [42,43] (Figure 6).
The population is scarce in these lands during the Middle Ages (1500–500 cal. BP) although some development is visible along the time. The Islamic Period (1239–850 cal. BP) brings a small demographic growth and the Muslim communities living in this mid-mountain system maintain spread cereal cultivation and livestock, since grazing is the most important economic resource for these communities. These activities are linked to the development of some roads connecting different Islamic territories crossing these mountains, used by flock movements and commercial routes [42,43,94,97,98,99]. However, the war is also present in the fire regime because fire is an element used for clearing the landscape and avoiding ambushes or forcing the enemy to move, as well as for destroying enemy resources. It is not until the ca. 650 cal. BP, when the peace is established in the territory, that the Castilian Kingdom can repopulate the taken lands and develop their economic activities. Among them, livestock is the most important practice for the economy. Especially relevant are the transhumance movements through the pacified territory, because the movement of those great flocks helped the landscape shaping [42,43,99,100].

3.4. The Fire That Creates: from 500 cal. BP to Present

In the Iberian Peninsula, the general trend shows an evident growth in fire activity in the three studied sites. This final period includes the Modern and Contemporary age (500–150 cal. BP; 150 cal. BP–present) and the climate is determined by the final stages of the Little Ice Age (600–150 cal. BP) [101]. However, human activities will become the most important driver in the fire regime during this period. This fact will be of such importance that these fires recorded in the first part of the Modern Age will be also present in the legal contemporary literature, since many of those documents try to regulate and control the use of fire [42,43,99,102].
El Perro mire shows a sharp increase in CHAR values (Figure 3) despite there being only one although large robust fire episode according to peak detection analysis. This is a consequence of the registration of the biggest fire events in this mire. In fact, this part of El Perro record is statistically different from the rest of the core, showing just a few similarities with the above-mentioned fire event in 2651 cal. BP. The Viñuelas and Brezoso also records moderate to high CHAR values (Figure 3), although fire frequencies seem to have descended (Figure A1, Supplementary material). In El Perro mire, the two registered episodes trigger different vegetation responses. There is a fire event in 333 cal. BP, which opens the woodlands meanwhile the shrubs show a descending trend. The human activity increases around this space from the Modern Age and the 19th century, with an evident growth in the anthropic herbs but especially in the grazing for livestock. This population growth forces the neighbors to include in the cultivation circuits some lands never exploited before due to their location, insufficient quality, or low fertility, putting pressure on the mountains and other spaces not exploited in the past [74]. On the other hand, the fire event in 92 cal. BP does not disentangle such a landscape change, but the woodlands closed and got denser, meanwhile the anthropic herbs show an evident descent.
Brezoso mire shows an increase in the fire activity, slight at the beginning of the period and sharper in the end. The final stage shows a further gentle increase, from ca. 300 cal. BP until the present, especially evident in the fire events recorded (Figure 3). However, despite the opening landscape it implies, the drop of fire events at the very end helped the arboreal layer take back the space. Thus, at the end of the period, the landscape is dominated by Quercus with a herb layer dominated by Poaceae and anthropic species. Cereal cultivation is also present around the mire surroundings, although not in a representative way [42,43] (Figure 6).
Finally, Viñuelas mire shows a very irregular CHAR trend (Figure 3) and the fire activity in this moment is also statistically different from the other events recorded, although the fire frequency descends (Figure A1 in Appendix A). The landscape is dominated by deciduous Quercus, but the evergreen Quercus is scarce. The shrublands have a representative presence and the herbs dominate the ensemble, although the livestock does not show a representative impact in the landscape [42,43] (Figure 6). The demographic growth and the political circumstances, and the owning changes coming from the Confiscation Laws promulgated against some ecclesiastic and public properties in the 19th century, promoted some intense changes in land use and forest exploitation. The 20th century, on the other hand, implies some population shift and rural abandonment, which has the above-mentioned consequences in the landscape. After some time with a descending fire regime and woodland spreading, the trend is increasing, favored by warm and dry climatic conditions and human hand or, in this case, desertion [21,42,43,74,103,104,105].

4. Conclusions

The reconstructed fire history for the Toledo Mountains highlights how human activities take on progressively more importance over time. Socio-economic changes were essential for defining fire regimes in this mid-elevation range, especially during recent times. Land use has determined vegetation composition, structure, and recent fires. Thus, fire history runs along the different cultural stages and human interests and actions [7,21,105,106].
This study dives into the fire regimes of the Toledo Mountains, following the trends present by El Perro, Brezoso, and Viñuelas mires’ records. This study shows the evident disagreement between the Iberian Peninsula general trend and the fire history recorded in the Toledo Mountains, where it is possible to distinguish three moments with different fire trends. In the first one (ca. 5000–3000 cal. BP), the climate has an indisputable effect on fire regimes, while in the second one (ca. 3000–500 cal. BP), human activities also regulate the fire trends, becoming the most important fire driver in the last five hundred years. In the very last one, from 500 cal. BP onwards, human activities—or the lack of them—determined the fire activity.
It is necessary to understand the fire regime of the past and the way it affected ecosystem dynamics to better understand the current fire events and their impact in the landscape [107,108]. More long-term fire ecology studies are needed to complete this essential task for environment management, especially in the Toledo Mountains, where there is still a lot of work to do.

Supplementary Materials

The following are available online at https://0-www-mdpi-com.brum.beds.ac.uk/2571-6255/2/4/54/s1.

Author Contributions

Conceptualization, R.L.-L., S.P.-D. and J.A.L.-S.; Data curation, R.L.-L. and C.M.-M.; Formal analysis, R.L.-L. and O.B.; Funding acquisition, J.A.L.-S.; Investigation, R.L.-L.; Resources, C.M.-M.; Software, O.B.; Supervision, O.B., S.P.-D. and J.A.L.-S.; Writing—original draft, R.L.-L.; Writing—review & editing, O.B., S.P.-D., C.M.-M. and J.A.L.-S.

Funding

This research was funded the project REDISCO-HAR2017-88035-P (Plan Nacional I+D+I, Spanish Ministry of Economy and Competitiveness). Reyes Luelmo-Lautenschlaeger is funded by an Formación del Profesorado Universitario (FPU) grant (Spanish Ministry of Education, Culture, and Sports).

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Figure A1. Fire Frequency in the selected mires.
Figure A1. Fire Frequency in the selected mires.
Fire 02 00054 g0a1
Figure A2. Fire densities and Robust Fire Events in the selected mires (red crosses): Brezoso, Viñuelas, El Perro.
Figure A2. Fire densities and Robust Fire Events in the selected mires (red crosses): Brezoso, Viñuelas, El Perro.
Fire 02 00054 g0a2
Figure A3. Kruskal–Wallis test from the studied sites.
Figure A3. Kruskal–Wallis test from the studied sites.
Fire 02 00054 g0a3

References

  1. Naveh, Z. The evolutionary significance of fire in the Mediterranean Region. Vegetatio 1975, 29, 199–208. [Google Scholar] [CrossRef]
  2. Caldadaro, N. Human ecological intervention and the role of forest fires in human ecology. Sci. Total Environ. 2002, 292, 141–165. [Google Scholar] [CrossRef]
  3. Pausas, J.; Llovet, J.; Rodrigo, A.; Vallejo, R. Are wildfires a disaster in the Mediterranean basin?—A review. Int. J. Wildland Fire 2008, 17, 713–723. [Google Scholar] [CrossRef]
  4. Pausas, J.G.; Keeley, J.E. A burning story: The role of fire in the history of life. BioScience 2009, 59, 593–601. [Google Scholar] [CrossRef]
  5. Pérez, B.; Cruz, A.; Fernández-González, F.; Moreno, J.M. Effects of the recent land-use history on the postfire vegetation of uplands in Central Spain. For. Ecol. Manag. 2003, 182, 273–283. [Google Scholar] [CrossRef]
  6. Rodrigo, A.; Retana, J.; Pico, F.J. Direct regeneration is not the only response of mediterranean forests to large fires. Ecology 2004, 85, 716–729. [Google Scholar] [CrossRef]
  7. Montiel-Molina, C.; Galiana-Martín, L. Fire Scenarios in Spain: A Territorial Approach to Proactive Fire Management in the Context of Global Change. Forests 2016, 7, 273. [Google Scholar] [CrossRef]
  8. Goren-Inbar, N. Evidence of Hominin Control of Fire at Gesher Benot Ya’aqov, Israel. Science 2004, 304, 725–727. [Google Scholar] [CrossRef] [PubMed]
  9. Burjachs, F.; Expósito, I. Charcoal and pollen analysis: Examples of Holocene fire dynamics in Mediterranean Iberian Peninsula. Catena 2015, 135, 340–349. [Google Scholar] [CrossRef]
  10. López-Sáez, J.A.; Vargas, G.; Ruiz-Fernández, J.; Blarquez, O.; Alba-Sánchez, F.; Oliva, M.; Pérez-Díaz, S.; Robles-López, S.; Abel-Schaad, D. Paleofire dynamics in central Spain during the Late Holocene: The role of climatic and antrhopogenic forcing. Land Degrad. Dev. 2018, 29, 2045–2059. [Google Scholar] [CrossRef]
  11. Bond, W.J.; Woodland, F.I.; Midgley, G.F. The global distribution of ecosystems in a world without fire. New Phytol. 2004, 165, 525–538. [Google Scholar] [CrossRef] [PubMed]
  12. Gil-Romera, G.; Carrión, J.S.; Pausas, J.G.; Sevilla-Callejo, M.; Lamb, H.F.; Fernández, S.; Burjarchs, F. Holocene fire activity and vegetation response in South-Eastern Iberia. Quat. Sci. Rev. 2010, 29, 1082–1092. [Google Scholar] [CrossRef] [Green Version]
  13. Vanière, B.; Blarquez, O.; Rius, D.; Doyen, E.; Brücher, T.; Colombaroli, D.; Connor, S.; Feurdean, A.; Hickler, T.; Kaltenrieder, P.; et al. 7000-year human legacy of elevation-dependent European fire regimes. Quat. Sci. Rev. 2016, 132, 206–212. [Google Scholar] [CrossRef]
  14. Morales-Molino, C.; Tinner, W.; García-Antón, M.; Colombaroli, D. The historical demise of Pinus nigra forests in the Northern Iberian Plateau (south-western Europe). J. Ecol. 2017, 105, 634–646. [Google Scholar] [CrossRef]
  15. Connor, S.E.; Vanière, B.; Colombaroli, D.; Anderson, R.S.; Carrión, J.S.; Ejarque, A.; Gil Romera, G.; González-Sampériz, P.; Hoefer, D.; Morales-Molino, C.; et al. Humans take control of fire-driven diversity changes in Mediterranean Iberia’s vegetation during the mid-late Holocene. Holocene 2019, 29, 886–901. [Google Scholar] [CrossRef]
  16. Daniau, A.L.; Sánchez-Goni, M.F.; Beaufort, D.; Loggoun-Défarge, F.; Loutre, M.F.; Duprat, J. Dansgaard-Oeschger climatic variability revealed by fire emissions in southwestern Iberia. Quat. Sci. Rev. 2007, 26, 1369–1383. [Google Scholar] [CrossRef]
  17. Delarze, R.; Caldelari, D.; Hainard, P. Effects of fire on forest dynamics in southern Switzerland. J. Veg. Sci. 1992, 3, 55–60. [Google Scholar] [CrossRef]
  18. Withlock, C.; Bartlein, P.J. Holocene fire activity as a record of past environmental change. Dev. Quat. Sci. 2004, 1, 479–490. [Google Scholar] [CrossRef]
  19. Bradtmöller, M.; Grimm, S.; Riel-Salvatore, J. Resilience Theory in Archaeological Practice. An Annotated Review. Quat. Int. 2017, 446, 3–16. [Google Scholar] [CrossRef]
  20. Pausas, J.G.; Keeley, J.E. Wildfires as an ecosystem service. Front. Ecol. Environ. 2019. [Google Scholar] [CrossRef]
  21. Ortega, M.; Saura, S.; González-Ávila, S.; Gómez-Sanz, V.; Elena-Rosselló, R. Landscape vulnerability to wildfires at the forest-agriculture interface: Half-century patterns in Spain assessed through the SISPARES monitoring framework. Agrofor. Syst. 2012, 85, 331–349. [Google Scholar] [CrossRef]
  22. Adolf, C.; Wunderle, S.; Colombaroli, D.; Weber, H.; Gobet, E.; Heiri, O.; van Leeuwen, J.F.N.; Bigler, C.; Connor, S.E.; Galka, M.; et al. The sedimentary and remote-sensing reflection of biomass burning in Europe. Glob. Ecol. Biogeogr. 2018, 27, 199–212. [Google Scholar] [CrossRef]
  23. Tianhua, H.; Lamont, B.B.; Pausas, J.G. Fire as a key driver of Earth’s biodiversity. Biol. Rev. 2009. [Google Scholar] [CrossRef]
  24. Vanière, B.; Power, M.J.; Roberts, N.; Tinner, W.; Carrión, J.S.; Magny, M.; Bartlein, P.; Colombaroli, D.; Daniau, A.L.; Fisinger, W.; et al. Circum-Mediterranean fire activity and climate changes during the mid-Holocene environmental transition (8500–2500 cal. BP). Holocene 2011, 21, 53–73. [Google Scholar] [CrossRef]
  25. Iglesias, V.; Withlock, C. Fire responses to postglacial climate change and human impact in northern Patagonia (41–43°S). Proc. Natl. Acad. Sci. USA 2014. [Google Scholar] [CrossRef] [PubMed]
  26. Higuera, P.; Brubaker, L.B.; Anderson, P.A.; Sheng Hu, F.; Brown, T.A. Vegetation mediated the impacts of postglacial climate change on fire regimes in the south-central Brooks Range, Alaska. Ecol. Monogr. 2009, 79, 201–209. [Google Scholar] [CrossRef]
  27. Gavin, D.G.; Sheng Hu, F.; Lertzman, K.; Corbert. Weak climatic control of stand-scale fire history during the Late Holocene. Ecology 2006, 87, 1722–1732. [Google Scholar] [CrossRef]
  28. Power, M.J.; Marlon, J.; Ortiz, N.; Bartlein, P.J.; Harrison, S.P.; Mayle, F.E.; Ballouche, A.; Bradshaw, R.H.W.; Carcaillet, C.; Cordova, C.; et al. Changes in fire regimes since the Last Glacial Maximum: An assessment based on a global synthesis and analysis of charcoal data. Clim. Dyn. 2008, 30, 887–907. [Google Scholar] [CrossRef]
  29. Whitlock, C.; Shafer, S.L.; Marlon, J. The role of climate and vegetation change in shaping past and future fire regimes in the northwestern US and the implications for ecosystem management. For. Ecol. Manag. 2003, 178, 5–21. [Google Scholar] [CrossRef]
  30. Higuera, P.E.; Gaving, 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. Wildland Fire 2010, 19, 996–1014. [Google Scholar] [CrossRef]
  31. Blarquez, O.; Girardin, M.P.; Leys, B.; Ali, A.A.; Aleman, J.C.; Bergeron, Y.; Carcaillet, C. Paleofire reconstruction based on an ensemble-member strategy applied to sedimentary charcoal. Geophys. Res. Lett. 2013, 40, 2667–2672. [Google Scholar] [CrossRef]
  32. Long, C.J.; Withlock, C. Fire and Vegetation History from the Coastal Rain Forest of the Western Oregon Coast Range. Quat. Res. 2002, 58, 215–225. [Google Scholar] [CrossRef] [Green Version]
  33. Carcaillet, C.; Bouvier, M.; Fréchette, B.; Laroouche, A.C.; Richard, P.J.H. Comparison of pollen-slide and sieving methods in lacustrine charcoal analyses for local and regional fire history. Holocene 2001, 11, 467–476. [Google Scholar] [CrossRef]
  34. Mooney, S.D.; Tinner, W. The analysis of charcoal in peat and organic sediments. Mir. Peat. 2011, 7, 1–18. [Google Scholar]
  35. Rémy, C.C.; Fouquemberg, C.; Asselin, H.; Andrieux, B.; Magnan, G.; Brossier, B.; Grondin, P.; Bergeron, Y.; Talon, B.; Girardin, M.P.; et al. Guidelines for the use and interpretation of palaeofire reconstructions based on various archives and proxies. Quat. Sci. Rev. 2018, 193, 312–322. [Google Scholar] [CrossRef]
  36. Whitlock, C.; Larsen, C. Charcoal as a fire proxy. In Tracking Environmental Change Using Lake Sediments, vol. 3, Terrestrial, Algal, and Siliceous Indicators; Smol, J.P., Birks, H.J.B., Last, W.M., Eds.; Kluwer: Dordrecht, The Netherlands, 2001; pp. 75–97. [Google Scholar] [CrossRef]
  37. Lynch, J.A.; Clark, J.S.; Stocks, B.J. Charcoal production, dispersal, and deposition from the Fort Providence experimental fire: Interpreting fire regimes from charcoal records in boreal forests. Can. J. For. Res. 2004, 34, 1642–1656. [Google Scholar] [CrossRef]
  38. Higuera, P.E.; Petters, M.E.; Brubaker, L.B.; Gavin, D.G. Understanding the origin and analysis of sediment-charcoal records with a simulation model. Quat. Sci. Rev. 2007, 26, 1790–1809. [Google Scholar] [CrossRef]
  39. Martín-Serrano, A.; Molina, E.; Nozal, F.; Carral, M.P. Itinerario A2. Transversal en los Montes de Toledo. In Itinerarios Geomorfológicos por Castilla-La Mancha: Libro de las Excursiones Desarrolladas Durante la VIII Reunión Nacional de Geomorfología, Celebrada en Toledo, 22–25 de Septiembre de 2004; Benito, G., Díez Herrero, A., Eds.; Sociedad Española de Geomorfología-CSIC Centro de Ciencias Medioambientales: Madrid, Spain, 2004; pp. 51–82. ISBN 849219586X. [Google Scholar]
  40. López-Sáez, J.A.; García-Río, R.; Alba-Sánchez, F.; García-Gómez, E.; Pérez-Díaz, S. Peatlands in the Toledo Mountains (central Spain): Characterisation and conservation status. Mires Peat 2014, 15, 1–23. [Google Scholar]
  41. San Miguel, A.; Rodríguez-Vigal, C.; Perea Garcí a-Calvo, R. Los Quintos de Mora. Gestión integral del monte mediterráneo. In Pastos, Paisajes Culturales Entre Tradición y Nuevos Paradigmas del Siglo XXI. Visitas de Campo; López-Carrasco, C., Rodríguez, M.P., San Miguel, A., Fernández, F., Roig, S., Eds.; Sociedad Española para el Estudio de los Pastos: Madrid, Spain, 2011; p. 704. ISBN 9788461487134. [Google Scholar]
  42. Morales-Molino, C.; Colombaroli, C.; Tinner, W.; Perea, R.; Valbuena-Carabaña, M.; Carrión, J.S.; Gil, L. Vegetation and fire dynamics during the last 4000 years in the Cabañeros National Park (central Spain). Rev. Paleobot. Paleoec. 2018, 253, 110–122. [Google Scholar] [CrossRef]
  43. Morales-Molino, C.; Tinner, W.; Perea, R.; Carrión, J.S.; Colombaroli, D.; Valbuena-Carabaña, M.; Zafra, E.; Gil, L. Unprecedented herbivory threatens rear-edge populations of Betula in southwestern Eurasia. Ecology 2019. [Google Scholar] [CrossRef] [PubMed]
  44. Ninyerola, M.; Roure, J.M.; Fernández, X.P. Atlas climático digital de la Península Ibérica: Metodología y aplicaciones en bioclimatología y geobotánica; Universitat Autònoma de Barcelona: Bellaterra, Spain, 2005; p. 45. ISBN 932860-8-7. [Google Scholar]
  45. Perea, D.F.; Perea, R. Vegetación y Flora de los Montes de Toledo; Ediciones Covarrubias: Toledo, Spain, 2008; p. 296. ISBN 9788493603519. [Google Scholar]
  46. Punt, W.; Marcks, A.; Hoen, P.P. Myricaceae. Rev. Paleaeobot. Palynol. 2002, 123, 99–105. [Google Scholar] [CrossRef]
  47. Luengo-Nicolau, E.; Sánchez-Mata, D. A hazel tree relict community (Corylus avellana L., Betulaceae) from the Guadiana River Middle Basin (Ciudad Real, Spain). Lanzaroa 2015, 36, 133–137. [Google Scholar] [CrossRef]
  48. Reimer, P.J.; Bard, E.; Bayliss, A.; Beck, J.W.; Blackwell, P.G.; Bronk Ramsey, C.; Buck, C.E.; Cheng, H.; Edwards, R.L.; Friedrich, M.; et al. Intcal13 and marine13 radiocarbon age calibration curves 0–50,000 years cal. BP. Radiocarbon 2013, 55, 1869–1887. [Google Scholar] [CrossRef]
  49. Hua, Q.; Barbetti, M. Review of tropospheric bomb 14C data for carbon cycle modelling and age calibration purposes. Radiocarbon 2004, 46, 1273–1298. [Google Scholar] [CrossRef]
  50. Blaauw, M. Methods and code for classical age-modelling of radiocarbon sequences. Quat. Geochron. 2010, 5, 512–518. [Google Scholar] [CrossRef]
  51. Blaaw, M. Classical Age-Depth Modelling of Cores from Deposits R library. Available online: https://CRAN.R-project.org/package=clam (accessed on 23 April 2019).
  52. Moore, P.D.; Webb, J.A.; Collinson, M.E. Pollen Analysis; Blackwell: London, UK, 1991; p. 216. ISBN 9780632021765. [Google Scholar]
  53. Goeury, C.; de Beaulieu, J.L. À propos de la concentration du pollen à l’aide de la liqueur de Thoulet dans les sédiments minéraux. Pollen Spores 1979, 21, 239–251. [Google Scholar]
  54. Stockmarr, J. Tablets with spores used in absolute pollen analysis. Pollen Spores 1971, 13, 614–621. [Google Scholar]
  55. Ramil-Rego, P.; Aira-Rodríguez, M.J.; Saá-Otero, P. Clave polínica de las Ericaceae gallegas. Lazaroa 1992, 13, 33–40. [Google Scholar]
  56. Beug, H. Leiftaden der Pollenbestimmung fürMitteleuropa un angrezende Gebiete; F. Pfeil Verlag: Munich, Germany, 2004. [Google Scholar]
  57. Reille, M. Pollen et spores d’Europe et d’Afrique du Nord, 2nd ed.; Laboratoire de Botanique Historique et Palynologie: Marseille, France, 1999; p. 543. ISBN 2950717500. [Google Scholar]
  58. Blackmore, S.; Steinmann, J.A.J.; Hoen, P.P.; Punt, W. Betulaceae and Corylaceae. Rev. Paleaeobot. Palynol. 2003, 123, 71–98. [Google Scholar] [CrossRef]
  59. Mateus, J.E. Pollen Morphography of Portuguese Ericales. Rev. Biol. 1989, 14, 135–208. [Google Scholar]
  60. Grimm, E.D. TGView; Illinois State Museum, Research and Collection Center: Springfield, MA, USA, 2004. [Google Scholar]
  61. Tinner, W.; van Leeuwen, J.F.N.; Colombaroli, D.; Vescovi, E.; van der Knaap, W.O.; Henne, P.D.; Pasta, S.; D’Angelo, S.; La Mantia, T. Holocene environmental and climatic changes at Gorgo Basso, a coastal lake in southern Sicily, Italy. Quat. Sci. Rev. 2009, 28, 1498–1510. [Google Scholar] [CrossRef]
  62. Higuera, P. CharAnalysis R library. Diagnostic and Analytical Tools for Sediment-Charcoal Analysis. Available online: http://Charanalysis.googlepages.com (accessed on 2 April 2019).
  63. Blarquez, O.; Vannière, B.; Marlon, J.R.; Daniau, A.-L.; Power, M.J.; Brewer, S.; Bartlein, P.J. Paleofire: An R package to analyse sedimentary charcoal records from the Global Charcoal Database to reconstruct past biomass burning. Comput. Geosci. 2014, 72, 255–261. [Google Scholar] [CrossRef]
  64. Gavin, D.G. K1D: Multivariate Ripley’s K-Function for One-Dimensional Data. 2010. Available online: http://citeseerx.ist.psu.edu/viewdoc/download?doi=10.1.1.483.6528&rep=rep1&type=pdf (accessed on 2 April 2019).
  65. Van Bellen, S.; Garneau, M.; Ali, A.A.; Bergeron, Y. Did fires drive Holocene carbon sequestration in boreal ombrotrophic peatlands of eastern Canada? Quat. Res. 2012, 78, 50–59. [Google Scholar] [CrossRef]
  66. Kelly, R.F.; Higuera, P.E.; Barrett, C.M.; Sheng Hu, F. A signal-to-noise index to quantify the potential for peak detection in sediment–charcoal records. Quat. Res. 2011, 75, 11–17. [Google Scholar] [CrossRef]
  67. Mudelsee, M.; Börbgen, M.; Tetzlaff, G.; Grünewald, U. Extreme floods in central Europe over the past 500 years. Role of cyclone pathway. “Zugstrasse Vb”. J. Geophys. Res. Atmos. 2004, 109, D23101. [Google Scholar] [CrossRef]
  68. Carcaillet, C.; Ali, A.A.; Blarquez, O.; Genries, A. Spatial variability of fire history in subalpine forests: From natural to cultural regimes. Ecoscience 2009, 16, 1–12. [Google Scholar] [CrossRef]
  69. Ali, A.A.; Carcaillet, C.; Bergeron, Y. Long-term fire-frequency variability in the eastern Canadian boreal forest: The influences of climate vs. local factors. Glob. Chang. Biol. 2009, 15, 1230–1241. [Google Scholar] [CrossRef]
  70. Daniau, A.L.; Bartlein, P.J.; Harrison, S.P.; Prentice, I.C.; Brewer, S.; Friedlingstein, P.; Harrison-Prentice, T.I.; Inoue, J.; Izumi, K.; Marlon, J.R.; et al. Predictability of biomass burning in response to climate changes. Glob. Biogeochem. Cycles 2012, 26. [Google Scholar] [CrossRef] [Green Version]
  71. Pausas, J.G. Simulating Mediterranean landscape pattern and vegetation dynamics under different fire regimes. Plant Ecology 2006, 187, 249–259. [Google Scholar] [CrossRef]
  72. Turner, R.; Roberts, N.; Jones, M.D. Climatic pacing of Mediterranean fire histories from lake sedimentary microcharcoal. Glob. Planet. Chang. 2008, 63, 317–324. [Google Scholar] [CrossRef]
  73. López-Sáez, J.A.; Abel-Schaad, D.; Pérez-Díaz, S.; Blanco-González, A.; Alba-Sánchez, F.; Dorado-Valiño, M.; Ruiz-Zapata, B.; Gil-García, M.J.; Gómez-González, C.; Franco-Múgica, F. Vegetation history, climate and human impact in the Spanish Central System over the last 9000 years. Quat. Int. 2014, 353, 98–122. [Google Scholar] [CrossRef]
  74. Lasanta-Martínez-Martínez, T. La transformación del paisaje en montaña media por la actividad agrícola en relación con las condiciones ambientales. In Acción Humana y Desertificación en Ambientes Mediterráneos; García-Ruiz, J.M., López-García, P., Eds.; Instituto Pirenaico de Ecología: Zaragoza, Spain, 1997; pp. 145–172. ISBN 9788492184224. [Google Scholar]
  75. Kohler, T.; Giger, M.; Hurni, H.; Ott, C.; Wiesmann, U.; Wyman-von Dach, S.; Maselli, D. Mountains and climate change: A global concern. Mt. Res. Dev. 2010, 30, 53–55. [Google Scholar] [CrossRef]
  76. Magny, M. Holocene climate variability as reflected by mid-European lake-level fluctuations and its probable impact on prehistoric human settlements. Quat. Int. 2004, 113, 65–79. [Google Scholar] [CrossRef]
  77. Lillios, K.T.; Blanco-González, A.; Lee, B.; López-Sáez, J.A. Mid-late Holocene climate, demography, and cultural dynamics in Iberia: A multiproxy approach. Quat. Sci. Rev. 2016, 135, 138–153. [Google Scholar] [CrossRef]
  78. López-Sáez, J.A.; Alba-Sánchez, F.; Nájera-Colino, T.; Molina-González, F.; Pérez-Díaz, S.; Sabariego-Ruiz, S. Paleoambiente y sociedad en la Edad del Bronce de la Mancha: La Motilla del Azuer. CPAG 2014, 24, 391–422. [Google Scholar]
  79. Bini, M.; Zanchetta, G.; Per+ou, A.; Cartier, R.; Català, A.; Cacho, I.; Dean, J.R.; Di Rita, F.; Drysdale, R.N.; Finnè, M.; et al. The 4.2 la BP Event in the Mediterranean Region: An overview. Clim. Past. 2019, 15, 555–577. [Google Scholar] [CrossRef]
  80. Schirrmacher, J.; Weinelt, M.; Blanz, T.; Andersen, N.; Salgueiro, E.; Schleider, R.R. Multi-decadal atmospheric and marine climate variability in southern Iberia during the mid-to late- Holocene. Clim. Past. 2019, 15, 617–634. [Google Scholar] [CrossRef]
  81. Sánchez-del Álamo, C.; Sardinero, S.; Bouso, V.; Hernández-Palacios, G.; Pérez-Badía, R.; Fernández-González, F. Los abedulares del Parque Nacional de Cabañeros: Sistemática, demografía, biología reproductiva y estrategias de conservación. In Proyectos de Investigación en Parques Nacionales: 2006–2009; Organismo Autónomo Parques Nacionale: Madrid, Spain, 2010; pp. 275–310. [Google Scholar]
  82. Behre, K.E. The interpretation of anthropogenic indicators in pollen diagrams. Pollen et Spores 1981, 23, 225–245. [Google Scholar]
  83. Blanco-González, A.; Lillios, K.T.; López-Sáez, J.A.; Drake, B.L. Cultural, demo-graphic and environmental dynamics of the Copper and Early Bronze Age in Iberia (3300–1500 BC): Towards an interregional multiproxy comparison at the time of the 4.2 ky BP event. J. World Prehist. 2018, 31, 1–79. [Google Scholar] [CrossRef]
  84. López-Sáez, J.A.; Pérez-Díaz, S.; Rodríguez-Ramírez, A.; Blanco-González, A.; Villarçias-Robles, J.R.; Luelmo-Lautenschlaeger, R.; Jiménez-Moreno, G.; Celestino-Pérez, S.; Cerrillo-Cuenca, E.; Pérez-Asensio, J.N.; et al. Mid-late Holocene environmental and cultural dynamics at the south-west tip of Europe (Doñana National Park, SW Iberia, Spain). J. Archaeol. Sci. Rep. 2018, 22, 58–78. [Google Scholar] [CrossRef]
  85. Moreno, A.; Pérez, A.; Frigola, J.; Nieto-Moreno, V.; Rodrigo-Gámiz, M.; Martrat, B.; González-Sampériz, P.; Morellón, M.; Martín-Puertas, C.; Corella, J.P.; et al. The Medieval Climate Anomaly in the Iberian Peninsula reconstructed from marine and lake records. Quat. Sci. Rev. 2012, 43, 16–32. [Google Scholar] [CrossRef]
  86. Martín-Puertas, C.; Valero-Garcés, B.L.; Brauer, A.; Mata, P.; Delgado-Huertas, A.; Dulski, P. The Iberian-Roman humid Period (2600-1600 cal yr BP) in the Zoñar Lake varve record (andalucía, southern Spain). Quat. Res. 2009, 71, 108–120. [Google Scholar] [CrossRef]
  87. Gil-García, M.J.; Ruiz Zapata, M.B.; Santisteban, J.I.; Mediavilla, R.; López-Pamo, E.; Dabrio, C.J. Late holocene environments in Las Tablas de Daimiel (south central Iberian Peninsula, Spain). Veget. Hist. Archaeobot. 2007, 16, 241–250. [Google Scholar] [CrossRef]
  88. Sánchez-López, G.; Hernández, A.; Pla-Rabes, S.; Trigo, R.M.; Toro, M.; Granados, I.; Sáez, A.; Masqué, P.; Pueyo, J.J.; Rubio-Inglés, M.J.; et al. Climate reconstruction for the las two millennia in central Iberia: The role of East Atlantic (EA), North Atlantic Oscillation (NAO) and their interplay over the Iberian Peninsula. Quat. Sci. Rev. 2016, 149, 135–150. [Google Scholar] [CrossRef]
  89. Fernández-Dávila, A. Paisaje y poblamiento en la Carpetania. Un territorio en definición. Zona Arqueol. 2014, 17, 45–70. [Google Scholar]
  90. López Sáez, J.A.; Blanco González, A.; López, L.; Ruiz, B.; Dorado, M.; Pérez, S.; Valdeolmillos, A.; Burjachs, F. Landscape and Climatic Changes during the End of the Late Prehistory in the Amblés Valley (Ávila, central Spain) from 1200 to 400 cal BC. Quat. Int. 2009, 90–101. [Google Scholar] [CrossRef]
  91. Joanin, S.; Magny, M.; Peyron, O.; Vannière, B.; Galop, D. Climate and land-use change during the late Holocene at Lake Ledro (southern Alps, Italy). Holocene 2014, 24, 591–602. [Google Scholar] [CrossRef]
  92. Van Geel, B.; Heijnis, H.; Charman, D.J.; Thompson, G.; Engels, S. Bog burst in the eastern Netherlands triggered by the 2.8 kyr BP climate event. Holocene 2014, 24, 1465–1477. [Google Scholar] [CrossRef]
  93. Jiménez-Morán, C. El poblamiento Visigodo en la zona oriental de la provincial de Toledo. In La Mancha occidental y la Mesa de Ocaña; Diputación Provincial de Toledo: Toledo, Spain, 2001; Volume 1, pp. 333–356. [Google Scholar]
  94. Blanco-González, A.; López-Sáez, J.A.; López-Merino, L. Ocupación y uso del territorio en el sector centromeridional de la cuenca del Duero entre la Antigüedad y la Alta Edad Media (siglos I-XI d.C.). Arch. Esp. Arqueol. 2009, 82, 275–300. [Google Scholar] [CrossRef]
  95. López-Sáez, J.A.; Peña-Chocarro, L.; López-Merino, L.; García-Gómez, E.; Pérez-Díaz, S.; García-Entero, V.; Castelo- Ruano, R. Paisajes culturales de las villas romanas de Toledo. Cuad. Soc. Esp. Cien. For. 2009, 30, 101–106. [Google Scholar]
  96. López-Sáez, J.A.; Serra-González, C.; Alba-Sánchez, F.; Robles-López, S.; Pérez-Díaz, S.; Abel-Schaad, D.; Glais, A. Exploring seven hundred years of transhumance, dynamic, fire and human activity through a historical mountain pass in central Spain. J. Mt. Sci. 2016, 13, 1139–1153. [Google Scholar] [CrossRef]
  97. Blanco-González, A.; López-Sáez, J.A.; Alba, F.; Abel-Schaad, D.; Pérez-Díaz, S. Medieval landscapes in the Spanish Central System (450–1350): A palaeoenvironmental and historical perspective. J. Mediev. Iber. Stud. 2015, 7, 1–17. [Google Scholar] [CrossRef]
  98. Ruibal, A. El Camino de Toledo a Córdoba por el paso de Alhover y su defensa. In Actas del II Congreso Internacional de Caminería Hispánica; Criado de Val, M.C., Ed.; AACHE Ediciones: Guadalajara, Spain, 1996; Tomo 2; pp. 37–52. ISBN 84-87743-66-8. [Google Scholar]
  99. Luelmo-Lautenschlaeger, R.; Pérez-Díaz, S.; Alba-Sánchez, F.; Abel-Schaad, D.; López-Sáez, J.A. Vegetation History in the Toledo Mountains (Central Iberia): Human Impact during the Last 1300 Years. Sustainability 2018, 10, 2575. [Google Scholar] [CrossRef]
  100. Izquierdo-Benito, R. Reconquista y Repoblación de la Tierra Toledana; Diputación provincial de Toledo-Instituto Provincial de Investigaciones y Estudios Toledanos: Toledo, Spain, 1983; 45p, ISBN 0211-4607. [Google Scholar]
  101. Oliva, M.; Ruiz-Fernández, J.; Barriendos, M.; Benito, G.; Cuadrat, J.M.; Domínguez-Castro, F.; García-Ruiz, J.M.; Giralt, S.; Gómez-Ortiz, A.; Hernández, A.; et al. The Little Ice Age in Iberian mountains. Earth-Sci. Rev. 2018, 177, 175–208. [Google Scholar] [CrossRef]
  102. Jiménez-de Gregorio, F. La Comarca Histórica Toledana de los Montes de Toledo; Instituto Provincial de Investigación y Estudios Toledanos-Diputación de Toledo: Toledo, Spain, 2008; p. 149. ISBN 8487103943. [Google Scholar]
  103. De Linares, V.G.G. Los bosques en España a lo largo de la Historia. In Historia de Los Bosques. El Significado de la Madera en el Desarrollo de la Civilización; Perlin, J., Ed.; Gaia Proyecto 2050: Madrid, Spain, 1999; pp. 429–480. ISBN 9788493023218. [Google Scholar]
  104. Bauer-Maanderscheid, E. Los Montes de España en la Historia; Ministerio de Agricultura: Madrid, Spain, 1980; 610p, ISBN 8474790840. [Google Scholar]
  105. Martínez, J.; Vega-García, C.; Chuvieco, E. Human-caused wildfire risk rating for prevention planning in Spain. J. Environ. Manag. 2009, 90, 1241–1252. [Google Scholar] [CrossRef] [PubMed]
  106. Blarquez, O.; Ali, A.A.; Girardin, M.P.; Grondin, P.; Fréchette, B.; Bergeron, Y.; Hély, C. Regional paleofire regimes affected by non-uniform climate, vegetation and human drivers. Sci. Rep. 2015, 5, 13356. [Google Scholar] [CrossRef] [PubMed]
  107. Conedera, M.; Tinner, W.; Nelf, C.; Meurer, M.; Dickens, A.F.; Krebs, P. Reconstructing past fire regimes: Methods, applications, and relevance to fire management and conservation. Quat. Sci. Rev. 2009, 28, 555–576. [Google Scholar] [CrossRef]
  108. Hennebelle, A.; Grondin, P.; Aleman, J.C.; Ali, A.A.; Bergeron, Y.; Borcard, D.; Blarquez, O. Using paleoecology to improve reference conditions for ecosystem-based management in western spruce-moss subdomain of Québec. For. Ecol. Manag. 2018, 430, 157–165. [Google Scholar] [CrossRef]
Figure 1. Map of the Toledo Mountains in central Spain with location of the selected peatlands (red dots): 1. Brezoso [42,43]; 2. Viñuelas [43]; 3. El Perro [this paper].
Figure 1. Map of the Toledo Mountains in central Spain with location of the selected peatlands (red dots): 1. Brezoso [42,43]; 2. Viñuelas [43]; 3. El Perro [this paper].
Fire 02 00054 g001
Figure 2. Age–depth model of El Perro mire. The red mark represents an outlier.
Figure 2. Age–depth model of El Perro mire. The red mark represents an outlier.
Fire 02 00054 g002
Figure 3. Charcoal accumulation rates (CHAR) and robust fire events for (top) Brezoso, (middle) Viñuelas, and (bottom) El Perro mires. “+” denotes statistically robust fire events. The grey areas correspond to raw charcoal data while the black lines correspond to CHAR interpolated to the median resolution of the records. Note the different scale of the Y axes.
Figure 3. Charcoal accumulation rates (CHAR) and robust fire events for (top) Brezoso, (middle) Viñuelas, and (bottom) El Perro mires. “+” denotes statistically robust fire events. The grey areas correspond to raw charcoal data while the black lines correspond to CHAR interpolated to the median resolution of the records. Note the different scale of the Y axes.
Fire 02 00054 g003
Figure 4. Fire synchronicity. (a) Timing of fire events in the three records. “+” denotes statistically robust fire events. (b) The L-function (transform of the K-function) for the events in (a) with 95% confidence envelope (thin grey lines) based on 1000 randomizations. The function exceeds the upper confidence envelop from 150 to 250 years, indicating strong correlation of event times within windows of that scale, but lack of long-term patterns in the three records for longer windows, i.e., independence in fire occurrence.
Figure 4. Fire synchronicity. (a) Timing of fire events in the three records. “+” denotes statistically robust fire events. (b) The L-function (transform of the K-function) for the events in (a) with 95% confidence envelope (thin grey lines) based on 1000 randomizations. The function exceeds the upper confidence envelop from 150 to 250 years, indicating strong correlation of event times within windows of that scale, but lack of long-term patterns in the three records for longer windows, i.e., independence in fire occurrence.
Fire 02 00054 g004
Figure 5. Composite curve of fire history in the three studied sites (a) and the entire Iberian Peninsula (b). Data obtained from the Global Charcoal Database [63], the complete list of sites used to construct the Iberian Peninsula composite curve are in supplementary Table S1. This study was combined with a Kruskal–Wallis tests to better perceive the similarities and differences between trends in both composites.
Figure 5. Composite curve of fire history in the three studied sites (a) and the entire Iberian Peninsula (b). Data obtained from the Global Charcoal Database [63], the complete list of sites used to construct the Iberian Peninsula composite curve are in supplementary Table S1. This study was combined with a Kruskal–Wallis tests to better perceive the similarities and differences between trends in both composites.
Fire 02 00054 g005
Figure 6. Summary percentage pollen diagrams of the (top) Brezoso, (middle) Viñuelas, and (bottom) El Perro mires. Dots represent percentages below 0.5%. Mediterranean shrubs: Cystus; Phillyrea. Anthropic herbs: Apiaceae; Asphodelus; Aster type; Boraginaceae; Cardueae; Centaurea; Cichorioideae; Convolvulus type; Dipsacus type; Echium. Anthropozoogenous herbs: Chenopodiaceae; Plantaginaceae; Urtica. Coprophilous fungi: Sordaria type (HdV-55); Sporomiella type (HdV-113). Taxonomical groups following Behre [82]. Grey lines are x3 exaggeration curves.
Figure 6. Summary percentage pollen diagrams of the (top) Brezoso, (middle) Viñuelas, and (bottom) El Perro mires. Dots represent percentages below 0.5%. Mediterranean shrubs: Cystus; Phillyrea. Anthropic herbs: Apiaceae; Asphodelus; Aster type; Boraginaceae; Cardueae; Centaurea; Cichorioideae; Convolvulus type; Dipsacus type; Echium. Anthropozoogenous herbs: Chenopodiaceae; Plantaginaceae; Urtica. Coprophilous fungi: Sordaria type (HdV-55); Sporomiella type (HdV-113). Taxonomical groups following Behre [82]. Grey lines are x3 exaggeration curves.
Fire 02 00054 g006
Table 1. Main features of the studied peatlands. Nutrients source: M = minerotrophic. Protection: RNP = reserve national park; NP = Non-protected.
Table 1. Main features of the studied peatlands. Nutrients source: M = minerotrophic. Protection: RNP = reserve national park; NP = Non-protected.
NameMunicipalityLatitudeLongitudeAltitude
(m.a.s.l.)
Area
(ha)
NutrientsProtection
BrezosoAlcoba de los Montes39°20′56.69″ N4°21′39.80″ W73310.62MRNP
ViñuelasAlcoba de los Montes39°22′28.00″ N4°29′18.08″ W7610.02MRNP
El PerroPuebla de don Rodrigo39°3′51.49″ N4°45′20.25″ W6903.74MNP
Table 2. AMS radiocarbon data with 2σ range of calibrated ages from El Perro mire.
Table 2. AMS radiocarbon data with 2σ range of calibrated ages from El Perro mire.
Laboratory CodeDepth (cm)14C Age (BP)Calibrated Age (cal. BP, 2σ)Calibrated age (cal. BP, Average Probability)
Ua-5529020122.3 ± 0.3−6.04–(−5.56)−5.8
Ua-5529140185 ± 250–294181
Poz-8425452955 ± 30796–927855
Ua-55292602345 ± 272324–24392352
Poz-84255682485 ± 302438–27242585
Ua-55293752594 ± 272719–27622743
Ua-55294843445 ± 313632–38283704
Poz-84256903830 ± 354099–44064232
Ua-55295954148 ± 314575–48244693
Poz-84257996470 ± 407293–74577376

Share and Cite

MDPI and ACS Style

Luelmo-Lautenschlaeger, R.; Blarquez, O.; Pérez-Díaz, S.; Morales-Molino, C.; López-Sáez, J.A. The Iberian Peninsula’s Burning Heart—Long-Term Fire History in the Toledo Mountains (Central Spain). Fire 2019, 2, 54. https://0-doi-org.brum.beds.ac.uk/10.3390/fire2040054

AMA Style

Luelmo-Lautenschlaeger R, Blarquez O, Pérez-Díaz S, Morales-Molino C, López-Sáez JA. The Iberian Peninsula’s Burning Heart—Long-Term Fire History in the Toledo Mountains (Central Spain). Fire. 2019; 2(4):54. https://0-doi-org.brum.beds.ac.uk/10.3390/fire2040054

Chicago/Turabian Style

Luelmo-Lautenschlaeger, Reyes, Olivier Blarquez, Sebastián Pérez-Díaz, César Morales-Molino, and José Antonio López-Sáez. 2019. "The Iberian Peninsula’s Burning Heart—Long-Term Fire History in the Toledo Mountains (Central Spain)" Fire 2, no. 4: 54. https://0-doi-org.brum.beds.ac.uk/10.3390/fire2040054

Article Metrics

Back to TopTop