Articles | Volume 18, issue 5
Research article
25 May 2022
Research article |  | 25 May 2022

Reconstructing burnt area during the Holocene: an Iberian case study

Yicheng Shen, Luke Sweeney, Mengmeng Liu, Jose Antonio Lopez Saez, Sebastián Pérez-Díaz, Reyes Luelmo-Lautenschlaeger, Graciela Gil-Romera, Dana Hoefer, Gonzalo Jiménez-Moreno, Heike Schneider, I. Colin Prentice, and Sandy P. Harrison

Charcoal accumulated in lake, bog or other anoxic sediments through time has been used to document the geographical patterns in changes in fire regimes. Such reconstructions are useful to explore the impact of climate and vegetation changes on fire during periods when human influence was less prevalent than today. However, charcoal records only provide semi-quantitative estimates of change in biomass burning. Here we derive quantitative estimates of burnt area from vegetation data in two stages. First, we relate the modern charcoal abundance to burnt area using a conversion factor derived from a generalised linear model of burnt area probability based on eight environmental predictors. Then, we establish the relationship between fossil pollen assemblages and burnt area using tolerance-weighted weighted averaging partial least-squares regression with a sampling frequency correction (fxTWA-PLS). We test this approach using the Iberian Peninsula as a case study because it is a fire-prone region with abundant pollen and charcoal records covering the Holocene. We derive the vegetation–burnt area relationship using the 31 records that have both modern and fossil charcoal and pollen data and then reconstruct palaeoburnt area for the 113 records with Holocene pollen records. The pollen data predict charcoal-derived burnt area relatively well (R2= 0.44), and the changes in reconstructed burnt area are synchronous with known climate changes through the Holocene. This new method opens up the possibility of reconstructing changes in fire regimes quantitatively from pollen records, after regional calibration of the vegetation–burnt area relationship, in regions where pollen records are more abundant than charcoal records.

1 Introduction

Fire is an important element in many ecosystems and in the Earth system (Bowman et al., 2009; Resco de Dios, 2020). It impacts vegetation dynamics, ecosystem functioning and biodiversity (Harrison et al., 2010; Ward et al., 2012; Keywood et al., 2013). It also affects climate through vegetation changes and the release of trace gases and aerosols. Fire directly impacts socio-economic assets (e.g. Stephenson et al., 2013; Thomas et al., 2017) and has deleterious effects on human health though the release of smoke and particulates into the atmosphere (e.g. Johnston et al., 2012; Yu et al., 2020). These impacts make it important to understand what controls the incidence and severity of fires.

Analyses of fire regimes during the satellite era have shown that multiple factors play a role in determining the occurrence of fire, including climate and fire weather, vegetation properties and human activities (e.g. Harrison et al., 2010; Brotons et al., 2013; Bistinas et al., 2014; Knorr et al., 2014; Andela et al., 2017; Forkel et al., 2019a, b; Kuhn-Régnier et al., 2020). However, the satellite record only covers a short time period (ca. 20 years), and the impact of anthropogenic changes to land cover on suppressing fire during this interval is strong (Andela et al., 2017). Reconstructing changing fire regimes during the pre-industrial Holocene (12 000 BP to ca. CE 1850) provides an opportunity to investigate the controls on fire over timescales when human influences on the landscape, including fire regimes, were more localised and less profound than they have become during the industrial era.

Sedimentary charcoal, preserved in lakes, peatbogs and other anoxic environments, has been widely used as an indicator of past changes in fire regimes (Marlon et al., 2008, 2016; Power et al., 2008; Daniau et al., 2012; Vannière et al., 2016; Connor et al., 2019). Evaluations that combine charcoal-inferred palaeofire reconstructions with past hydrological, vegetation and archaeological data support the idea that there are strong relationships among climate, fire, vegetation and human activities (Carrión et al., 2007; Marlon et al., 2008; Gil-Romera et al., 2010; Turner et al., 2010; Vannière et al., 2011; López-Sáez et al., 2018; Morales-Molino et al., 2018). However, charcoal records only provide a semi-quantitative index of fire activity rather than quantitative estimates of burnt area or biomass loss. Attempts to calibrate the charcoal record to provide quantitative estimates of proximity or area burnt are either site-specific (Duffin et al., 2008; Hennebelle et al., 2020) or rely on modelling (Higuera et al., 2007). Furthermore, although the number of charcoal records is increasing, there are still comparatively few sites compared to other types of palaeoenvironmental data, and this can make it difficult to make regional reconstructions of changing fire regimes.

Although the occurrence of fire is influenced by multiple factors, analyses of present-day fire relationships globally using satellite-derived data have shown that vegetation properties determining fuel availability are the strongest determinants of fire occurrence (Bistinas et al., 2014; Forkel et al., 2019a, b; Kuhn-Régnier et al., 2020). This suggests that palaeovegetation data could provide a way of reconstructing burnt area in the past, particularly at times when human influences on land cover were less important. This would also allow us to capitalise on the more extensive site networks for palaeovegetation.

In this study, we present a new method to reconstruct quantitative changes in fire regimes over the Holocene. We relate the relative scale of modern charcoal abundance to absolute burnt area using a conversion factor derived from a generalised linear model (GLM) of fire probability based on burnt area data. We then derive quantitative relationships between pollen assemblages and inferred burnt area using tolerance-weighted weighted averaging partial least-squares regression with a sampling frequency correction (fxTWA-PLS; Liu et al., 2020). The vegetation–burnt area relationship is then used to reconstruct changes in burnt area through time from pollen assemblages, including at sites with no charcoal record. We use the Iberian Peninsula as a test case. The Iberian Peninsula is the most fire-affected region in southern Europe (Jesus et al., 2019; Molina-Terrén et al., 2019). Although the modern fire regime is partly driven by human activities, the patterns also reflect the strong climate gradients across the region. Although much of the Iberian Peninsula has a typical Mediterranean climate, parts of the region are influenced by proximity to the Atlantic Ocean or the Mediterranean Sea and by the mountainous topography, giving rise to complex weather and climate patterns and large gradients in vegetation diversity (Loidi, 2017). We reconstruct fire regimes across the Iberian Peninsula through the Holocene and discuss the implications of the reconstructed changes.

Figure 1Flow chart of the methodology of burnt area reconstructions.


2 Methods

The central premise of our approach is that fire frequency is one of the factors that influences vegetation assemblages (see the Supplement) and therefore that specific aspects of differences in vegetation assemblages – identified by a numerical technique that can isolate the effects of any one controlling factor on taxon composition – can be used to reconstruct fire. The vegetation–fire relationship can be derived by comparing changes in pollen assemblages and charcoal records through time. However, since the charcoal records from different sites consist of different size fractions, and the records must be normalised to facilitate comparisons, it is necessary to derive site-specific conversion factors between modern charcoal abundance and present-day burnt area fraction. This calibration is then applied to the charcoal record in order to derive an estimate of the palaeoburnt area for each pollen sample.

2.1 Iberian pollen and charcoal data

Pollen data were obtained from the European Pollen Database (EPD;, (last access: 5 April 2021) or provided by the authors (Table S2 in the Supplement). Non-pollen palynomorphs (e.g. fungi and algae), introduced species and fire-insensitive plants (e.g. obligate aquatics) were removed from the assemblages before analysis. Some pollen taxa are not identified consistently by palynologists or occur at very few sites, so some pollen types were amalgamated to higher taxonomic groups (mostly genera for trees and families for herbaceous taxa) for consistency across the records (Table S3). Charcoal data were obtained from the Global Charcoal Database (Power et al., 2010; Marlon et al., 2016) or provided by the authors (Table S4). The original age models for both the pollen and the charcoal records were constructed using different methods and different calibrations of radiometric to calendar ages. We created new age models for all the records using the IntCal20 calibration curve (Reimer et al., 2020) and the BACON Bayesian age-modelling tool in the rbacon package (2.5.0) in CRAN (Blaauw and Christeny, 2011). Charcoal concentration data were converted to charcoal accumulation rate (influx: particlescm-2yr-1) before analysis by multiplying concentration by the background sedimentation rate.

2.2 Development of the generalised linear model

We obtained modern burnt area for the Iberian Peninsula from the fourth version of the Global Fire Emissions Database (GEFD4) (Randerson et al., 2017). The GLM was initially developed using 13 environmental variables covering climate, vegetation and human activities (Table S5). Some environmental data sets were only available at 0.5× 0.5 resolution, so all data sets were aggregated to this resolution using bilinear interpolation prior to analysis. Analyses were done for the common period between the data sets (January 2001 to December 2016) using annual values of all variables. The GLM was run using the stats package in R (version R.3.6.0) and used the logit link function and assumed a quasi-binomial distribution (R Core Team, 2019). We tested combinations of environmental predictors and selected the most parsimonious model with statistically significant variables and high prediction ability as assessed using pseudo-R2 (McFadden, 1973). The GLM-fitted burnt area was disaggregated from 0.5× 0.5 to 0.0083× 0.0083 by bilinear interpolation in order to extract present-day burnt area at each of the sites with modern charcoal records.

Figure 2Map showing the location of the 31 entities (31 sites) with modern charcoal used to derive the fire–vegetation relationship (red triangles) and the 113 entities (111 sites) used for burnt area reconstructions (blue circles) in the Iberian Peninsula.

2.3 Quantitative reconstructions of burnt area

We derived the relationship between the pollen assemblage and burnt area using 31 records with modern pollen and modern charcoal (Fig. 2). Rare pollen taxa with fewer than or equal to five occurrences in the data set were removed because they have been shown to have little predictive power in WA-PLS climate reconstructions (Turner et al., 2021). The charcoal records included some sites with only macroscopic and some with only microscopic charcoal. Since this had little impact on the patterns of change through time (Fig. S4 in the Supplement), we used both types, although we used macroscopic charcoal at sites with both macroscopic and microscopic charcoal. The sampling resolution varies between the individual records. To ensure comparability across records, the charcoal and pollen data were temporally binned prior to analysis: the modern bin covers the post-industrial period (CE 1850 to the present); a 100-year bin width was used for earlier intervals. Pollen counts were summed and converted to percentage of the total count in each bin. To standardise the values for different charcoal measurement units, we used a maximum transformation to convert mean charcoal accumulation rates to a 0–1 range (Eq. 1).

(1) x i , j = x i , j x j , max ,

where xi,j is the transformed value of the ith sample (xi) in the jth entity. xj,max is the maximum value of all samples in this entity.

The maximum transformation resulted in a similar scale of variability between entities in fire-prone areas and areas with little fire. We therefore applied a conversion factor to rescale the relative charcoal abundance to absolute burnt area for each of the records:

(2) conversion factor j = present-day burnt area fraction j modern charcoal data j ,

where modern charcoal data are the core-top binned charcoal data in the jth entity, and the present-day burnt area fraction in the jth entity was obtained from the GLM.

The palaeoburnt area fraction for the ith sample in the jth entity was then derived by multiplying the conversion factor of the jth entity by the charcoal value for this sample (Eq. 3):

(3) palaeoburnt area fraction i , j = conversion  factor j × charcoal  data i , j .

We applied Box–Cox transformation (Box and Cox, 1964) with λ= 0.25 to the palaeoburnt area fraction in order to reduce skewness prior to the fxTWA-PLS analyses (see the Supplement). The fire–vegetation relationship was determined using the last significant component in fxTWA-PLS, assessed using the p value, to avoid overfitting.

We applied the fxTWA-PLS-derived relationship between pollen abundance and burnt area to the binned pollen data from the 113 pollen records available from the Iberian Peninsula (Fig. 2) to reconstruct changes in fire regimes through the Holocene. Some of these 113 entities included pollen taxa that were not present in the data used to derive the vegetation–burnt area relationship; these taxa were therefore removed prior to analysis. We used composite plots and maps of specific times to show the spatial and temporal changes of reconstructed palaeofire regimes through the Holocene. We used loess smoothing with a window half-width of 300 years to construct the composite plots, with the uncertainty of the reconstruction estimated by bootstrap resampling of the individual reconstructions 1000 times (Efron, 1979; Efron and Tibshirani, 1993). We tested the robustness of our method by comparing the reconstructed burnt area composite with the trends shown by raw charcoal data for those records with fossil charcoal, where the uncertainty is estimated again by bootstrap resampling of the individual charcoal records 1000 times. Maps were created using the reconstructed burnt area for individual sites in the bin covering the time period of interest.

Table 1Generalised linear model of the modern burnt area fraction.

Notes:  p<0.1. * p<0.05. ** p<0.01. *** p<0.001.

Download Print Version | Download XLSX

Figure 3Mean (over 16 years) of observed (a) and fitted (b) values of burnt area fraction.


3 Results

3.1 The GLM

Several of the environmental predictors of burnt area were highly correlated to one another (Fig. S1). We tested the impact of including/removing highly and moderately correlated variables before selecting the final GLM (see the Supplement). The final model was constructed using eight variables (Tables 1 and S6 in the Supplement) and has a pseudo-R2= 0.20. However, all of these variables show statistically significant relationships (P< 0.1) with burnt area, and most have p values < 0.05. Gross primary production (GPP) shows a very strong positive relationship (t= 10.10) with burnt area fraction (Table 1, Fig. S2). Dry days per month (t= 8.46) and non-tree cover (t= 7.34) also show strong positive relationships with burnt area (Table 1, Fig. S2). These relationships make sense given that much of the Iberian Peninsula is relatively arid: increasing GPP and increasing non-tree cover are indices of increased fuel availability in arid, fuel-limited regions and promote increased burnt area. The number of dry days per month determines fuel dryness, and hence there is a positive relationship between the number of dry days and the burnt area. Predictions of burnt area from the final model show reasonably good agreement with the observed average burnt area (Fig. 3). Hotelling's T-squared test shows that there is no statistically significant difference between observed and fitted values (p value = 1). Both the observations and the model show the highest burnt area in northern Portugal and moderate burnt area in southern Portugal. Both observed and simulated burnt area are low along the northern coast and in the Pyrenees where fire is limited by wet conditions and in the dry interior where fire is limited by fuel availability.

Table 2Leave-out cross-validation fitness of the fxTWA-PLS method, showing results for all the components. The last significant number of components are shown in bold. RMSEP is the root mean square error of prediction. ΔRMSEP is the percent change of RMSEP using the current number of components (ncomp) than using one component less. b0, b1, b0⋅SE, and b1⋅SE are the intercept, slope, standard error of the intercept and standard error of the slope of the linear regression using the cross-validation result and burnt area data converted from charcoal abundance.

Download Print Version | Download XLSX

Figure 4The fitted plot and residual plot of TWA-PLS method, with fx correction. Panel (a) is the reconstructed burnt area using the last significant number of components, which is four here. The x axis is the burnt area fraction derived from charcoal data, and the y axis is the burnt area fraction reconstructed from pollen data using TWA-PLS with fx correction. The 1:1 line is shown in black, and the linear regression line is shown in red to show the degree of overall compression. Panel (b) shows the residuals of reconstructed burnt area fraction using the last significant number of components. The x axis is the burnt area fraction derived from charcoal data, and the y axis is the residual of burnt area reconstruction using TWA-PLS with fx correction. The zero line is shown in black, and the locally estimated scatterplot smoothing is shown in red to show the degree of local compression. The low compression zone is shown by grey shading.


3.2 The pollen–burnt area relationship

The fxTWA-PLS-derived relationship, based on the last (fourth) significant component, has good predictive power (R2= 0.44) (Table 2). A linear regression of the cross-validation results and the burnt area data has a slope of 0.526 (Table 2, Fig. 4a), which shows that the degree of overall compression towards the centre of the sampled range is relatively low. The degree of local compression, which is assessed by whether the residuals are around zero across the burnt area range in locally estimated scatterplot smoothing, indicates that the low-compression zone where reconstructed values after Box–Cox transformation are most reliable is between 3.25 and 2.5, in other words, between 0.12 % and 1.98 % of the grid cell area (Fig. 4b). Comparison with results using WA-PLS and tolerance-weighted WA-PLS (TWA-PLS) confirms that fxTWA-PLS produces a large reduction in compression in the central part of the burnt-area range and has a higher predictive power (Table S8, Figs. S5 and S6). However, although fxTWA-PLS reduces the compression bias, it does not remove it completely: burnt area is overestimated at the low end and underestimated at the high end of the burnt area (Fig. 4).

Figure 5Composite plots comparing maximum-transformed charcoal values and the reconstructed burnt area for these entities for the 51 entities with charcoal. Maximum-transformed charcoal is shown in blue; burnt area fraction is shown in red. The loess smoothing is done with a span of 0.04.


Charcoal values are not expected to be directly comparable with the reconstructed burnt area but should show comparable temporal trends. A composite plot of reconstructed burnt area for the 51 entities that have both pollen records used to reconstruct burnt area and charcoal records, and therefore can be compared, shows similar trends to the composite plot derived from the maximum-transformed charcoal (Fig. 5). This suggests there is little distortion of the signal caused by deriving burnt area using the fxTWA-PLS relationship.

Figure 6Composite curve of reconstructed burnt area using fxTWA-PLS, using the locfit() function with a half-width of 300 and 1000 bootstrap samples. The locally estimated scatterplot smoothing is shown in blue. The upper and lower 95th percentile confidence intervals are shown in grey.


Figure 7Spatial patterns of reconstructed burnt area fraction at key times in the Holocene.

3.3 Fire history of the Iberian Peninsula through the Holocene

The composite plot based on all 113 pollen records from the Iberian Peninsula (Fig. 6) shows a moderate peak in burnt area around 13 ka followed by a marked increase at 12.5 ka and a subsequent peak in burnt area around 11.5 ka. Although this early part of the record is based on relatively few sites, and so the confidence intervals are large, the pattern corresponds to high fire activity during the Bølling–Allerød (14.6–12.9 ka) warm interval and low fire activity during much of the Younger Dryas (12.9–11.7 ka) cold phase, with an increase in burnt area associated with the rapid warming at the end of the Younger Dryas. Burnt area is relatively low at the beginning of the Holocene. Although there is a gradual increase in burnt area between 9 and 0.6 ka, the burnt area fraction is lower than present until at least 2 ka. The increase in burnt area is quite marked after around 4.5 ka and peaks at 0.6 ka. The burnt area fraction at 0.6 ka is larger than at any time in the record. Burnt area declines after 0.6 ka, although the modern reconstructed value is still higher than the values obtained for most of the Holocene.

The spatial coverage of sites (Fig. 7) for the earlier part of the record is sparse, but coverage is good from 7 ka onwards. The pattern of lower burnt area in eastern than in western Iberia, seen in the modern observations, is generally preserved both in high and low fire intervals. However, some of the records from northern Iberia (e.g. Saldropo, Puerto de Los Tornos) show extremely high burnt area which exceeds the scope of the low-compression zone during the last millennium, and particularly at 0.6 ka. This may reflect the persistent bias at the high end of the fx-TWA-PLS reconstruction range.

4 Discussion

We have shown that it is possible to derive trends in burnt area through time by applying a quantitative relationship between pollen assemblages and charcoal-derived burnt area to palaeovegetation records from the Iberian Peninsula. Our analyses exploit the multivariate nature of vegetation and hence pollen assemblages. Vegetation patterns, and the distribution of individual species, are controlled by many factors including seasonal temperature and precipitation regimes, disturbance (including wildfires) and human activities. Pollen-based palaeoclimate methods have long exploited the multivariate nature of pollen assemblages to reconstruct different aspects of climate (see, for example, the discussion in Bartlein et al., 2011). The canonical correspondence analysis (CCA) shows that there is sufficient information in the pollen assemblages to assess the independent contribution of fire to vegetation assemblages. The overall relationship between pollen and charcoal-derived burnt area is reasonably strong (R2= 0.44), reflecting the importance of vegetation properties (gross primary production and non-tree cover) in driving the occurrence of fire – as seen in the GLM analysis of satellite-derived modern burnt area patterns. The overwhelming importance of vegetation properties in influencing modern fire occurrence is consistent with results from global analyses (e.g. Moritz et al., 2012; Pausas and Ribeiro, 2013; Bistinas et al., 2014; Forkel et al., 2019b). Nevertheless, the GLM analysis shows that climate factors, in particular the occurrence of dry intervals, are important controls on modern fire patterns in Iberia. Again, this is consistent with global analyses of the modern drivers of fire occurrence.

We have used fxTWA-PLS (Liu et al., 2020) to make the burnt area reconstructions because this technique reduces the compression bias characteristic of many other reconstruction techniques by accounting for differences in the tolerance of individual taxa in an assemblage and for the frequency of the reconstructed variable in the training data set. However, although the bias is apparently reduced, there is still an overestimation at the low end and an underestimation at the high end of the burnt area range. This is reflected by the extremely high burnt area values reconstructed for some sites in northern Iberia in the recent millennium which exceed the upper limit of the low-compression zone. This remaining bias may also be explained by the comparatively small sample size (1106 binned samples) compared with the much larger data set of 6458 samples used by Liu et al. (2020) to carry out climate reconstructions for Eurasia. It would be useful to test whether the problem of compression bias in the reconstruction of burnt area could be overcome by expanding the training data set to cover a wider range of vegetation types and fire regimes.

Although palaeoburnt area reconstructions have only been obtained from a limited number of records, they nevertheless show interesting patterns over the past ca. 15 kyr. The high fire intervals at the beginning of the record, between 14–13 ka and between 12–11 ka, correspond to the Bølling–Allerød (14.6–12.9 ka) warming interval and to rapid warming at the end of the subsequent Younger Dryas (12.9–11.7 ka) cold phase (Alley et al., 1993). A similar response to these climate events has been seen in charcoal records from eastern North America (Marlon et al., 2009). Burnt area is less than today through the Early and Middle Holocene (10–5 ka), an interval when pollen, speleothem and lake records suggest the Mediterranean region was wetter than today (Prentice et al., 1996; Magny et al., 2002; Bartlein et al., 2011; Roberts et al., 2011). Reconstructions of fire activity anomalies (FAAs) for the south-eastern part of the Iberia Peninsula also indicate low-level fire activity in the Mid-Holocene between 7.5 and 6 ka (Gil-Romera et al., 2010). Burnt area continuously increases during the Medieval Warm Period (MWP; 1–0.7 ka) and peaks at 0.6 ka, consistent with the warm and dry conditions recorded during this period in the Iberian Peninsula (Moreno et al., 2012). During the Little Ice Age (LIA; 0.6–0.1 ka), the reconstructed fire indicates a sharp downturn, which may be associated with subsequent cold and wet climate (Ramos-Román et al., 2016; Abrantes et al., 2017). Thus, the broadscale patterns of trends in reconstructed burnt area are consistent with known Holocene climate changes in this region.

There is a distinct west–east gradient in burnt area across Iberia today, and this gradient of high fire in the west and less fire in the east is also present during other intervals of the Holocene. This pattern is likely related to the regional gradient in fuel availability and drought (Pausas and Fernández-Muñoz, 2012). However, the west–east gradient in burnt area is less pronounced during the Mid-Holocene, consistent with a less pronounced gradient in precipitation and moisture availability shown by other studies (e.g. González-Sampériz et al., 2017; Liu, 2019). Reconstructed patterns in fire were also more homogenous after 1 ka, and again this is consistent with the fact that temperature and humidity gradients were less pronounced at that time than they are today (Sánchez-López et al., 2016).

Our analyses show that climate, and climate-induced changes in vegetation, have influenced the fire regimes of the Iberian Peninsula during the Holocene. However, many studies have suggested that human activities could also have been important (Blanco-González et al., 2018; Connor et al., 2019; Feurdean et al., 2020). Land clearance during the Neolithic agricultural transition has been associated with increases in fire activity in some sites from the Iberian Peninsula (e.g. García-Ruiz et al., 2016; Carracedo et al., 2018). Although initiation of agriculture was not synchronous everywhere, the regional onset of agriculture is registered around 7.5 ka (Zapata et al., 2004; Fyfe et al., 2019; Harrison et al., 2020) when the burnt area reconstructions do not indicate high fire activity. However, the gradual increase in reconstructed burnt area between 5 and 0.6 ka may be an indication of increasing human activity, since the initial increase is broadly consistent with increased population shown by summed probability distributions (SPDs) of radiocarbon dates (Balsera et al., 2015; Lillios et al., 2016; Harrison et al., 2020). Human activities, such as deforestation and appropriation of land for agriculture, may have been an important driver of fire patterns from the Bronze Age onwards (Morales-Molino et al., 2013; Morales-Molino and García-Antón, 2014; González-Sampériz et al., 2017), while the competing effects of land abandonment and fire suppression may have contributed to the changes in burnt area in recent times (Turco et al., 2016; Silva et al., 2019). Nevertheless, our GLM analysis indicates that the intensity of human influence, as measured by crop or grazing land area or by population density, consistently has a negative effect on burnt area under modern conditions. It seems likely that human influence on Holocene fire regimes may have been complex, with agricultural expansion both promoting and suppressing fire occurrence. More detailed comparisons of the reconstructed burnt area and archaeological data are required to test this.

The limited availability of charcoal records has meant that the analysis of past fire regimes has tended to focus on large-scale zonal or continental-scale patterns (e.g. Marlon et al., 2008; Power et al., 2008; Daniau et al., 2010; Vannière et al., 2011). Our new methodology opens up the possibility of reconstructing changes in fire regimes from pollen data and thus of examining finer-scale patterning that might reflect climate or human influences on fire. Spatially explicit reconstructions of burnt area would also be useful to evaluate the simulated response of fire to changing environmental drivers in the past (Thonicke et al., 2005; Brücher et al., 2014; Martin Calvo et al., 2014; Marlon et al., 2016; Kraaij et al., 2020) since comparisons based on qualitative inferences from charcoal are inconclusive (e.g. Brücher et al., 2014).

5 Conclusion

We have developed a novel method to reconstruct palaeoburnt area quantitatively from vegetation records, based on fire–vegetation relationships derived using fxTWA-PLS and the calibration of modern charcoal against GLM modelling of modern burnt area. We have applied this approach to reconstruct changes in burnt area through the Holocene for the Iberian Peninsula. The good predictive power of the fxTWA-PLS-derived fire–vegetation relationship and the plausibility of the palaeofire reconstructions with respect to known climate changes in the region suggest that this calibration approach could be applied more generally to provide quantitative reconstructions of past fire regimes in other regions where there are limited charcoal data, and pollen data are more abundant.

Code and data availability

The pollen and charcoal data from the Iberian Peninsula used in this analysis are available from Harrison et al. (2022; All other data used are publicly accessible. The code used to generate the new-age models (ageR) was created by Villegas-Diaz et al. (2021; and is available from (last access: 9 May 2022). This Github repository contains code scripts created by Shen (2022;


The supplement related to this article is available online at:

Author contributions

YS, ICP and SPH designed this study. JALS, SPD, RLL, GJM, DH, HS and GGR contributed pollen and charcoal data. YS and LS developed the new pollen and charcoal age models. YS carried out the analyses. YS and SPH wrote the first draft of the manuscript, and all authors contributed to the final version.

Competing interests

The contact author has declared that neither they nor their co-authors have any competing interests.


Publisher's note: Copernicus Publications remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.


Yicheng Shen and Sandy P. Harrison acknowledge support from the ERC-funded project GC 2.0 (Global Change 2.0: Unlocking the past for a clearer future; grant no. 94481). I. Colin Prentice acknowledges support from the ERC under the European Union Horizon 2020 Research and Innovation programme (grant no. 787203 REALM). Luke Sweeney acknowledges support from the Leverhulme Centre for Wildfires, Environment and Society. Mengmeng Liu acknowledges support from Imperial College through the Lee Family Scholarship. José Antonio López-Sáez acknowledges support from the REDISCO-HAR2017-88035-P (Plan Nacional I+D+I, Spanish Ministry of Economy and Competitiveness) project. Reyes Luelmo is funded by a FPU grant. Some of the pollen data used in the analyses were extracted from the European Pollen Database (EPD;, last access: 5 April 2021), and the work of the data contributors and the EPD community is gratefully acknowledged. Some of the charcoal data were extracted from the Global Charcoal Database (, last access: 5 April 2021), and we gratefully acknowledge contributors to this effort and the curators of the database. We thank colleagues in the Leverhulme Centre for Wildfires, Environment and Society (, last access: 10 May 2022) and from the SPECIAL group at the University of Reading (, last access: 10 May 2022) for discussions during the development of this work.

Financial support

This research has been supported by the European Research Council (GC2.0 (grant no. 694481)), the European Research Council (REALM (grant no. 787203)), Imperial College through the Lee Family Scholarship, the Leverhulme Centre for Wildfires, Environment and Society (grant no. RC-2018-023), and the REDISCO (grant no. HAR2017-88035-P) project.

Review statement

This paper was edited by Keely Mills and reviewed by two anonymous referees.


Abrantes, F., Rodrigues, T., Rufino, M., Salgueiro, E., Oliveira, D., Gomes, S., Oliveira, P., Costa, A., Mil-Homens, M., Drago, T., and Naughton, F.: The climate of the Common Era off the Iberian Peninsula, Clim. Past, 13, 1901–1918,, 2017. 

Alley, R. B., Meese, D. A., Shuman, C. A., Gow, A. J., Taylor, K. C., Grootes, P. M., White, J. W. C., Ram, M., Waddington, E. D., Mayewski, P. A., and Zielinski, G. A.: Abrupt increase in Greenland snow accumulation at the end of the Younger Dryas event, Nature, 362, 527–529,, 1993. 

Andela, N., Morton, D. C., Giglio, L., Chen, Y., Van Der Werf, G. R., Kasibhatla, P. S., DeFries, R. S., Collatz, G. J., Hantson, S., Kloster, S., Bachelet, D., Forrest, M., Lasslop, G., Li, F., Mangeon, S., Melton, J. R., Yue, C., and Randerson, J. T.: A human-driven decline in global burned area, Science, 356, 1356–1362,, 2017. 

Balsera, V., Díaz-del-Río, P., Gilman, A., Uriarte, A., and Vicent, J. M.: Approaching the demography of late prehistoric Iberia through summed calibrated date probability distributions (7000–2000 cal BC), Quatern. Int., 386, 208–211,, 2015. 

Bartlein, P. J., Harrison, S. P., Brewer, S., Connor, S., Davis, B. A. S., Gajewski, K., Guiot, J., Harrison-Prentice, T. I., Henderson, A., Peyron, O., Prentice, I. C., Scholze, M., Seppä, H., Shuman, B., Sugita, S., Thompson, R. S., Viau, A. E., Williams, J., and Wu, H.: Pollen-based continental climate reconstructions at 6 and 21 ka: A global synthesis, Clim. Dynam., 37, 775–802,, 2011. 

Bistinas, I., Harrison, S. P., Prentice, I. C., and Pereira, J. M. C.: Causal relationships versus emergent patterns in the global controls of fire frequency, Biogeosciences, 11, 5087–5101,, 2014. 

Blaauw, M. and Christeny, J. A.: Flexible paleoclimate age-depth models using an autoregressive gamma process, Bayesian Anal., 6, 457–474,, 2011. 

Blanco-González, A., Lillios, K. T., López-Sáez, J. A., and Drake, B. L.: Cultural, demographic 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., 31, 1–79,, 2018. 

Bowman, D. M. J. S., Balch, J. K., Artaxo, P., Bond, W. J., Carlson, J. M., Cochrane, M. A., D'Antonio, C. M., DeFries, R. S., Doyle, J. C., Harrison, S. P., Johnston, F. H., Keeley, J. E., Krawchuk, M. A., Kull, C. A., Marston, J. B., Moritz, M. A., Prentice, I. C., Roos, C. I., Scott, A. C., Swetnam, T. W., Van Der Werf, G. R., and Pyne, S. J.: Fire in the earth system, Science, 324, 481–484,, 2009. 

Box, G. E. P. and Cox, D. R.: An analysis of transformations, J. R. Stat. Soc. B, 26, 211–234, 1964. 

Brotons, L., Aquilué, N., de Cáceres, M., Fortin, M. J., and Fall, A.: How fire history, fire suppression practices and climate change affect wildfire regimes in Mediterranean landscapes, PLoS One, 8, e62392,, 2013. 

Brücher, T., Brovkin, V., Kloster, S., Marlon, J. R., and Power, M. J.: Comparing modelled fire dynamics with charcoal records for the Holocene, Clim. Past, 10, 811–824,, 2014. 

Carracedo, V., Cunill, R., García-Codron, J. C., Pèlachs, A., Pérez-Obiol, R., and Soriano, J. M.: History of fires and vegetation since the Neolithic in the Cantabrian Mountains (Spain), Land Degrad. Dev., 29, 2060–2072,, 2018. 

Carrión, J. S., Fuentes, N., González-Sampériz, P., Sánchez Quirante, L., Finlayson, J. C., Fernández, S., and Andrade, A.: Holocene environmental change in a montane region of southern Europe with a long history of human settlement, Quaternary Sci. Rev., 26, 1455–1475,, 2007. 

Connor, S. E., Vanniè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., Revelles, J., Schneider, H., van der Knaap, W. O., van Leeuwen, J. F., and Woodbridge, J.: Humans take control of fire-driven diversity changes in Mediterranean Iberia's vegetation during the mid–late Holocene, Holocene, 29, 886–901,, 2019. 

Daniau, A. L., D'Errico, F., and Sánchez Goñi, M. F.: Testing the hypothesis of fire use for ecosystem management by Neanderthal and Upper Palaeolithic modern human populations, PLoS One, 5,, 2010. 

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., Mooney, S., Power, M. J., Stevenson, J., Tinner, W., Andrič, M., Atanassova, J., Behling, H., Black, M., Blarquez, O., Brown, K. J., Carcaillet, C., Colhoun, E. A., Colombaroli, D., Davis, B. A. S., D'Costa, D., Dodson, J., Dupont, L., Eshetu, Z., Gavin, D. G., Genries, A., Haberle, S., Hallett, D. J., Hope, G., Horn, S. P., Kassa, T. G., Katamura, F., Kennedy, L. M., Kershaw, P., Krivonogov, S., Long, C., Magri, D., Marinova, E., McKenzie, G. M., Moreno, P. I., Moss, P., Neumann, F. H., Norstrm, E., Paitre, C., Rius, D., Roberts, N., Robinson, G. S., Sasaki, N., Scott, L., Takahara, H., Terwilliger, V., Thevenon, F., Turner, R., Valsecchi, V. G., Vannière, B., Walsh, M., Williams, N., and Zhang, Y.: Predictability of biomass burning in response to climate changes, Global Biogeochem. Cy., 26, GB4007,, 2012. 

Duffin, K. I., Gillson, L., and Willis, K. J.: Testing the sensitivity of charcoal as an indicator of fire events in savanna environments: quantitative predictions of fire proximity, area and intensity, Holocene, 18, 279–291,, 2008. 

Efron, B.: Bootstrap Methods: Another Look at the Jackknife, Ann. Stat., 7, 1–26,, 1979. 

Efron, B. and Tibshirani, R. J.: An Introduction to the Bootstrap, Chapman and Hall/CRC, Boca Raton, Florida, ISBN 0-412-04231-2, 1993. 

Feurdean, A., Vannière, B., Finsinger, W., Warren, D., Connor, S. C., Forrest, M., Liakka, J., Panait, A., Werner, C., Andrič, M., Bobek, P., Carter, V. A., Davis, B., Diaconu, A.-C., Dietze, E., Feeser, I., Florescu, G., Gałka, M., Giesecke, T., Jahns, S., Jamrichová, E., Kajukało, K., Kaplan, J., Karpińska-Kołaczek, M., Kołaczek, P., Kuneš, P., Kupriyanov, D., Lamentowicz, M., Lemmen, C., Magyari, E. K., Marcisz, K., Marinova, E., Niamir, A., Novenko, E., Obremska, M., Pȩdziszewska, A., Pfeiffer, M., Poska, A., Rösch, M., Słowiński, M., Stančikaitė, M., Szal, M., Świȩta-Musznicka, J., Tanţău, I., Theuerkauf, M., Tonkov, S., Valkó, O., Vassiljev, J., Veski, S., Vincze, I., Wacnik, A., Wiethold, J., and Hickler, T.: Fire hazard modulation by long-term dynamics in land cover and dominant forest type in eastern and central Europe, Biogeosciences, 17, 1213–1230,, 2020. 

Forkel, M., Andela, N., Harrison, S. P., Lasslop, G., van Marle, M., Chuvieco, E., Dorigo, W., Forrest, M., Hantson, S., Heil, A., Li, F., Melton, J., Sitch, S., Yue, C., and Arneth, A.: Emergent relationships with respect to burned area in global satellite observations and fire-enabled vegetation models, Biogeosciences, 16, 57–76,, 2019a. 

Forkel, M., Dorigo, W., Lasslop, G., Chuvieco, E., Hantson, S., Heil, A., Teubner, I., Thonicke, K., and Harrison, S. P.: Recent global and regional trends in burned area and their compensating environmental controls, Environ. Res. Commun., 1, 051005,, 2019b. 

Fyfe, R. M., Woodbridge, J., Palmisano, A., Bevan, A., Shennan, S., Burjachs, F., Legarra Herrero, B., García Puchol, O., Carrión, J.-S., Revelles, J., and Roberts, C. N.: Prehistoric palaeodemographics and regional land cover change in eastern Iberia, Holocene, 29, 799–815,, 2019. 

García-Ruiz, J. M., Sanjuán, Y., Gil-Romera, G., González-Sampériz, P., Beguería, S., Arnáez, J., Coba-Pérez, P., Gómez-Villar, A., Álvarez-Martínez, J., Lana-Renault, N., Pérez-Cardiel, E., and López de Calle, C.: Mid and late Holocene forest fires and deforestation in the subalpine belt of the Iberian range, northern Spain, J. Mt. Sci., 13, 1760–1772,, 2016. 

Gil-Romera, G., Carrión, J. S., Pausas, J. G., Sevilla-Callejo, M., Lamb, H. F., Fernández, S., and Burjachs, F.: Holocene fire activity and vegetation response in South-Eastern Iberia, Quaternary Sci. Rev., 29, 1082–1092,, 2010. 

González-Sampériz, P., Aranbarri, J., Pérez-Sanz, A., Gil-Romera, G., Moreno, A., Leunda, M., Sevilla-Callejo, M., Corella, J. P., Morellón, M., Oliva, B., and Valero-Garcés, B.: 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. 

Harrison, S. P., Marlon, J. R., and Bartlein, P. J.: Fire in the Earth System, in: Changing Climates, Earth Systems and Society, edited by: Dodson, J., Springer Netherlands, Dordrecht, 21–48,, 2010. 

Harrison, S., Shen, Y., and Sweeney, L.: Pollen data and charcoal data of the Iberian Peninsula (version 3), University of Reading [data set],, 2022. 

Harrison, S. P., Gaillard, M.-J., Stocker, B. D., Vander Linden, M., Klein Goldewijk, K., Boles, O., Braconnot, P., Dawson, A., Fluet-Chouinard, E., Kaplan, J. O., Kastner, T., Pausata, F. S. R., Robinson, E., Whitehouse, N. J., Madella, M., and Morrison, K. D.: Development and testing scenarios for implementing land use and land cover changes during the Holocene in Earth system model experiments, Geosci. Model Dev., 13, 805–824,, 2020. 

Hennebelle, A., Aleman, J. C., Ali, A. A., Bergeron, Y., Carcaillet, C., Grondin, P., Landry, J., and Blarquez, O.: The reconstruction of burned area and fire severity using charcoal from boreal lake sediments, Holocene, 30, 1400–1409,, 2020. 

Higuera, P. E., Peters, M. E., Brubaker, L. B., and Gavin, D. G.: Understanding the origin and analysis of sediment-charcoal records with a simulation model, Quaternary Sci. Rev., 26, 1790–1809,, 2007. 

Johnston, F. H., Henderson, S. B., Chen, Y., Randerson, J. T., Marlier, M., DeFries, R. S., Kinney, P., Bowman, D. M. J. S., and Brauer, M.: Estimated global mortality attributable to smoke from landscape fires, Environ. Health Persp., 120, 695–701,, 2012. 

Keywood, M., Kanakidou, M., Stohl, A., Dentener, F., Grassi, G., Meyer, C. P., Torseth, K., Edwards, D., Thompson, A. M., Lohmann, U., and Burrows, J.: Fire in the air: biomass burning impacts in a changing climate, Crit. Rev. Env. Sci. Tech., 43, 40–83,, 2013. 

Knorr, W., Kaminski, T., Arneth, A., and Weber, U.: Impact of human population density on fire frequency at the global scale, Biogeosciences, 11, 1085–1102,, 2014. 

Kraaij, T., Engelbrecht, F., Franklin, J., and Cowling, R. M.: A fiery past: A comparison of glacial and contemporary fire regimes on the Palaeo-Agulhas Plain, Cape Floristic Region, Quaternary Sci. Rev., 235, 106059,, 2020. 

Kuhn-Régnier, A., Voulgarakis, A., Nowack, P., Forkel, M., Prentice, I. C., and Harrison, S. P.: The importance of antecedent vegetation and drought conditions as global drivers of burnt area, Biogeosciences, 18, 3861–3879,, 2021. 

Lillios, K. T., Blanco-González, A., Drake, B. L., and López-Sáez, J. A.: Mid-late Holocene climate, demography, and cultural dynamics in Iberia: A multi-proxy approach, Quaternary Sci. Rev., 135, 138–153,, 2016. 

Liu, M.: A theory of palaeoclimate reconstruction from biotic indicators: Application to Holocene pollen records from the Iberian Peninsula, thesis, Imperial College London, 2019. 

Liu, M., Prentice, I. C., Ter Braak, C. J. F., and Harrison, S. P.: An improved statistical approach for reconstructing past climates from biotic assemblages: Improved palaeoclimate reconstruction, P. Roy. Soc. A-Math., 476,, 2020. 

Loidi, J. (Ed.): The Vegetation of the Iberian Peninsula, 1st edn., Springer International Publishing, Cham, Switzerland,, 2017. 

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., and Abel-Schaad, D.: Paleofire dynamics in Central Spain during the Late Holocene: The role of climatic and anthropogenic forcing, Land Degrad. Dev., 29, 2045–2059,, 2018. 

Magny, M., Miramont, C., and Sivan, O.: Assessment of the impact of climate and anthropogenic factors on Holocene Mediterranean vegetation in Europe on the basis of palaeohydrological records, Palaeogeogr. Palaeocl., 186, 47–59,, 2002. 

Marlon, J. R., Bartlein, P. J., Carcaillet, C., Gavin, D. G., Harrison, S. P., Higuera, P. E., Joos, F., Power, M. J., and Prentice, I. C.: Climate and human influences on global biomass burning over the past two millennia, Nat. Geosci., 1, 697–702,, 2008. 

Marlon, J. R., Bartlein, P. J., Walsh, M. K., Harrison, S. P., Brown, K. J., Edwards, M. E., Higuera, P. E., Power, M. J., Anderson, R. S., Briles, C., Brunelle, A., Carcaillet, C., Daniels, M., Hu, F. S., Lavoie, M., Long, C., Minckley, T., Richard, P. J. H., Scott, A. C., Shafer, D. S., Tinner, W., Umbanhowar, C. E., and Whitlock, C.: Wildfire responses to abrupt climate change in North America, P. Natl. Acad. Sci. USA, 106, 2519–2524,, 2009. 

Marlon, J. R., Kelly, R., Daniau, A.-L., Vannière, B., Power, M. J., Bartlein, P., Higuera, P., Blarquez, O., Brewer, S., Brücher, T., Feurdean, A., Romera, G. G., Iglesias, V., Maezumi, S. Y., Magi, B., Courtney Mustaphi, C. J., and Zhihai, T.: Reconstructions of biomass burning from sediment-charcoal records to improve data–model comparisons, Biogeosciences, 13, 3225–3244,, 2016. 

Martin Calvo, M., Prentice, I. C., and Harrison, S. P.: Climate versus carbon dioxide controls on biomass burning: a model analysis of the glacial–interglacial contrast, Biogeosciences, 11, 6017–6027,, 2014. 

McFadden, D. L.: Conditional logit analysis of qualitative choice behavior.: Chapter 4, in: Frontiers in Econometrics, edited by: Zarembka, P., Academic Press Inc, New York, 105–142, 1973. 

Molina-Terrén, D. M., Xanthopoulos, G., Diakakis, M., Ribeiro, L., Caballero, D., Delogu, G. M., Viegas, D. X., Silva, C. A., and Cardil, A.: Analysis of forest fire fatalities in Southern Europe: Spain, Portugal, Greece and Sardinia (Italy), Int. J. Wildl. Fire, 28, 85–98,, 2019. 

Morales-Molino, C. and García-Antón, M.: Vegetation and fire history since the last glacial maximum in an inland area of the western Mediterranean Basin (Northern Iberian Plateau, NW Spain), Quaternary Res., 81, 63–77,, 2014. 

Morales-Molino, C., García-Antón, M., Postigo-Mijarra, J. M., and Morla, C.: Holocene vegetation, fire and climate interactions on the westernmost fringe of the Mediterranean Basin, Quaternary Sci. Rev., 59, 5–17,, 2013. 

Morales-Molino, C., Colombaroli, D., Tinner, W., Perea, R., Valbuena-Carabaña, M., Carrión, J. S., and Gil, L.: Vegetation and fire dynamics during the last 4000 years in the Cabañeros National Park (central Spain), Rev. Palaeobot. Palyno., 253, 110–122,, 2018. 

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., Belmonte, Á., Sancho, C., Cacho, I., Herrera, G., Canals, M., Grimalt, J. O., Jiménez-Espejo, F., Martínez-Ruiz, F., Vegas-Vilarrúbia, T., and Valero-Garcés, B. L.: The Medieval Climate Anomaly in the Iberian Peninsula reconstructed from marine and lake records, Quaternary Sci. Rev., 43, 16–32,, 2012. 

Moritz, M. A., Parisien, M.-A., Batllori, E., Krawchuk, M. A., Van Dorn, J., Ganz, D. J., and Hayhoe, K.: Climate change and disruptions to global fire activity, Ecosphere, 3, 49,, 2012. 

Pausas, J. G. and Fernández-Muñoz, S.: Fire regime changes in the Western Mediterranean Basin: From fuel-limited to drought-driven fire regime, Climatic Change, 110, 215–226,, 2012. 

Pausas, J. G. and Ribeiro, E.: The global fire-productivity relationship, Global. Ecol. Biogeogr., 22, 728–736,, 2013. 

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., Mooney, S., Moreno, P. I., Prentice, I. C., Thonicke, K., Tinner, W., Whitlock, C., Zhang, Y., Zhao, Y., Ali, A. A., Anderson, R. S., Beer, R., Behling, H., Briles, C., Brown, K. J., Brunelle, A., Bush, M., Camill, P., Chu, G. Q., Clark, J., Colombaroli, D., Connor, S., Daniau, A. L., Daniels, M., Dodson, J., Doughty, E., Edwards, M. E., Finsinger, W., Foster, D., Frechette, J., Gaillard, M. J., Gavin, D. G., Gobet, E., Haberle, S., Hallett, D. J., Higuera, P., Hope, G., Horn, S., Inoue, J., Kaltenrieder, P., Kennedy, L., Kong, Z. C., Larsen, C., Long, C. J., Lynch, J., Lynch, E. A., McGlone, M., Meeks, S., Mensing, S., Meyer, G., Minckley, T., Mohr, J., Nelson, D. M., New, J., Newnham, R., Noti, R., Oswald, W., Pierce, J., Richard, P. J. H., Rowe, C., Sanchez Goñi, M. F., Shuman, B. N., Takahara, H., Toney, J., Turney, C., Urrego-Sanchez, D. H., Umbanhowar, C., Vandergoes, M., Vanniere, B., Vescovi, E., Walsh, M., Wang, X., Williams, N., Wilmshurst, J., and Zhang, J. H.: Changes in fire regimes since the last glacial maximum: An assessment based on a global synthesis and analysis of charcoal data, Clim. Dynam., 30, 887–907,, 2008. 

Power, M. J., Marlon, J. R., Bartlein, P. J., and Harrison, S. P.: Fire history and the global charcoal database: A new tool for hypothesis testing and data exploration, Palaeogeogr. Palaeocl., 291, 52–59,, 2010. 

Prentice, I. C., Guiot, J., Huntley, B., Jolly, D., and Cheddadi, R.: Reconstructing biomes from palaeoecological data: A general method and its application to European pollen data at 0 and 6 ka, Clim. Dynam., 12, 185–194,, 1996. 

R Core Team: R: A language and environment for statistical computing, R Foundation for Statistical Computing, Vienna, Austria, (last access: 5 April 2021), 2019. 

Ramos-Román, M. J., Jiménez-Moreno, G., Anderson, R. S., García-Alix, A., Toney, J. L., Jiménez-Espejo, F. J., and Carrión, J. S.: Centennial-scale vegetation and North Atlantic Oscillation changes during the Late Holocene in the southern Iberia, Quaternary Sci. Rev., 143, 84–95,, 2016. 

Randerson, J. T., van der Werf, G. R., Giglio, L., Collatz, G. J., and Kasibhatla, P. S.: Global Fire Emissions Database, Version 4.1 (GFEDv4), ORNL DAAC [data set], Oak Ridge, Tennessee, USA,, 2017. 

Reimer, P. J., Austin, W. E. N., Bard, E., Bayliss, A., Blackwell, P. G., Bronk Ramsey, C., Butzin, M., Cheng, H., Edwards, R. L., Friedrich, M., Grootes, P. M., Guilderson, T. P., Hajdas, I., Heaton, T. J., Hogg, A. G., Hughen, K. A., Kromer, B., Manning, S. W., Muscheler, R., Palmer, J. G., Pearson, C., Van Der Plicht, J., Reimer, R. W., Richards, D. A., Scott, E. M., Southon, J. R., Turney, C. S. M., Wacker, L., Adolphi, F., Büntgen, U., Capano, M., Fahrni, S. M., Fogtmann-Schulz, A., Friedrich, R., Köhler, P., Kudsk, S., Miyake, F., Olsen, J., Reinig, F., Sakamoto, M., Sookdeo, A., and Talamo, S.: The IntCal20 Northern Hemisphere Radiocarbon Age Calibration Curve (0–55 cal kBP), Radiocarbon, 62, 725–757,, 2020. 

Resco de Dios, V.: Plant-Fire Interactions: Applying Ecophysiology to Wildfire Management, 1st edn., edited by: von Gadow, K., Pukkala, T., and Tomé, M., Springer Nature, Cham, Switzerland,, 2020. 

Roberts, N., Brayshaw, D., Kuzucuoğlu, C., Perez, R., and Sadori, L.: The mid-Holocene climatic transition in the Mediterranean: Causes and consequences, Holocene, 21, 3–13,, 2011. 

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., and Giralt, S.: Climate reconstruction for the last two millennia in central Iberia: The role of East Atlantic (EA), North Atlantic Oscillation (NAO) and their interplay over the Iberian Peninsula, Quaternary Sci. Rev., 149, 135–150,, 2016. 

San-Miguel-Ayanz, J., Durrant, T., Boca, R., Libertà, G., Branco, A., de Rigo, D., Ferrari, D., Maianti, P., Artés Vivancos, T., Oom, D., Pfeiffer, H., Nuijten, D., and Leray, T.: Forest fires in Europe, Middle East and North Africa 2018, Scientific and Technical Research series, 107 pp.,, 2019. 

Shen, Y.: Yicheng-Shen/Burnt-area-reconstruction: Burnt area reconstruction using R, Version v1.0.0, Zenodo [code],, 2022. 

Silva, J. M. N., Moreno, M. V., Le Page, Y., Oom, D., Bistinas, I., and Pereira, J. M. C.: Spatiotemporal trends of area burnt in the Iberian Peninsula, 1975–2013, Reg. Environ. Change, 19, 515–527,, 2019. 

Stephenson, C., Handmer, J., and Betts, R.: Estimating the economic, social and environmental impacts of wildfires in Australia, Environ. Hazards-UK, 12, 93–111,, 2013.  

Thomas, D., Butry, D., Gilbert, S., Webb, D., and Fung, J.: The costs and losses of wildfires: A literature survey, NIST Special Publication 1215, 72 pp.,, 2017. 

Thonicke, K., Prentice, I. C., and Hewitt, C.: Modeling glacial-interglacial changes in global fire regimes and trace gas emissions, Global Biogeochem. Cy., 19, 1–10,, 2005. 

Turco, M., Bedia, J., Di Liberto, F., Fiorucci, P., von Hardenberg, J., Koutsias, N., Llasat, M.-C., Xystrakis, F., and Provenzale, A.: Decreasing fires in Mediterranean Europe, edited by: Carcaillet, C., PLoS One, 11, e0150663,, 2016. 

Turner, M. G., Wei, D., Prentice, I. C., and Harrison, S. P.: The impact of methodological decisions on climate reconstructions using WA-PLS, Quaternary Res., 99, 341–356,, 2021. 

Turner, R., Roberts, N., Eastwood, W. J., Jenkins, E., and Rosen, A.: Fire, climate and the origins of agriculture: Micro-charcoal records of biomass burning during the last glacial-interglacial transition in Southwest Asia, J. Quaternary Sci., 25, 371–386,, 2010. 

Vannière, B., Power, M. J., Roberts, N., Tinner, W., Carrión, J., Magny, M., Bartlein, P., Colombaroli, D., Daniau, A. L., Finsinger, W., Gil-Romera, G., Kaltenrieder, P., Pini, R., Sadori, L., Turner, R., Valsecchi, V., and Vescovi, E.: Circum-Mediterranean fire activity and climate changes during the mid-Holocene environmental transition (8500–2500 cal. BP), Holocene, 21, 53–73,, 2011. 

Vannière, B., Blarquez, O., Rius, D., Doyen, E., Brücher, T., Colombaroli, D., Connor, S., Feurdean, A., Hickler, T., Kaltenrieder, P., Lemmen, C., Leys, B., Massa, C., and Olofsson, J.: 7000-year human legacy of elevation-dependent European fire regimes, Quaternary Sci. Rev., 132, 206–212,, 2016. 

Villegas-Diaz, R., Cruz-Silva, E., and Harrison, S. P.: ageR: Supervised Age Models, Zenodo [code],, 2021. 

Ward, D. S., Kloster, S., Mahowald, N. M., Rogers, B. M., Randerson, J. T., and Hess, P. G.: The changing radiative forcing of fires: global model estimates for past, present and future, Atmos. Chem. Phys., 12, 10857–10886,, 2012. 

Yu, P., Xu, R., Abramson, M. J., Li, S., and Guo, Y.: Bushfires in Australia: a serious health emergency under climate change, Lancet Planet. Heal., 4, e7–e8,, 2020. 

Zapata, L., Peña-Chocarro, L., Pérez-Jordá, G., and Stika, H. P.: Early neolithic agriculture in the iberian peninsula, J. World Prehist., 18, 283–325,, 2004. 

Short summary
We present a method to reconstruct burnt area using a relationship between pollen and charcoal abundances and the calibration of charcoal abundance using modern observations of burnt area. We use this method to reconstruct changes in burnt area over the past 12 000 years from sites in Iberia. We show that regional changes in burnt area reflect known changes in climate, with a high burnt area during warming intervals and low burnt area when the climate was cooler and/or wetter than today.