Monthly wildfire burned area frequency is here modeled with a power law distribution and scaling exponent across different European biomes are estimated. Data sets, spanning from 2000 to 2009, comprehend the inventory of monthly burned areas from the European Forest Fire Information System (EFFIS) and simulated monthly burned areas from a recent parameterization of a Land Surface Model (LSM), that is the Community Land Model (CLM). Power law exponents are estimated with a Maximum Likelihood Estimation (MLE) for different European biomes. The characteristic fire size (CFS), i.e. the area that most contributes to the total burned area, was also calculated both from EFFIS and CLM data set. We used the power law fitting and the CFS analysis to benchmark CLM model against the EFFIS observational wildfires data set available for Europe.
Results for the EFFIS data showed that power law fittings holds for 2–3 orders of magnitude in the Boreal and Continental ecoregions, whereas the distribution of the Alpine, Atlantic are fitted only in the upper tail. Power law instead is not a suitable model for fitting CLM simulations.
CLM benchmarking analysis showed that the model strongly overestimates burned areas and fails in reproducing size-frequency distribution of observed EFFIS wildfires. This benchmarking analysis showed that some refinements in CLM structure (in particular regarding the anthropogenic influence) are needed for predicting future wildfires scenarios, since the low spatial resolution of the model and differences in relative frequency of small and large fires can affect the reliability of the predictions.
Wildfire size frequency has been studied since fire data have been regularly collected and then analyzed to gain a mechanistic understanding of the underlying process and to develop predictive models of this complex natural phenomenon (Cumming, 2001; McKenzie et al., 2011; Moritz et al., 2005; Zinck and Grimm, 2009). The spatio-temporal distribution of wildfires in a certain area can be used to characterize the fire regime. Although widely used, this term still lacks a clear and widely accepted definition (Krebs et al., 2010). Fire regimes are mainly a function of land cover (e.g. type of fuel) and climate (temperature, precipitation, duration of the dry-wet season, etc.), but they also depend on socio-economic factors such as population density, agricultural practices and so on (Bowman et al., 2011; Pausas and Keeley, 2009). Recent climate change has been found to impact wildfire regimes and future increases in global temperature coupled with more frequent droughts is expected to increase fire activity (Flannigan et al., 2000), in particular in Mediterranean landscapes (Pausas, 2004). Despite some attempts based both on satellite (Chuvieco et al., 2008) and inventory data (Malamud et al., 2005a), to date there is no well-established/quantitative method for classifying fire regimes across different climates. A proper classification method is needed for the quantification of climate change impact on ecosystems and wildfire scenarios (Migliavacca et al., 2013b; Yang et al., 2014).
Many difficulties are encountered in the modeling of wildfires because of the intrinsic complexity of this system, in which many environmental and social variables are involved in the ignition, propagation and suppression of fires. Some models were formulated to describe and predict wildfire spatial and temporal frequencies. According to Zinck and Grimm (2009) these models can be classified in two categories: the fire ecology approach and the statistical physics approach. The former class of models aims to reproduce environmental interactions between the physical status of vegetation and forcing agents that lead to wildfires (Peterson, 2002; Thonicke et al., 2008, 2010); the latter class of models analyzes emergent properties of wildfire using information from fire catalogs and models them as cellular automata (Bak et al., 1987, 1990; Drossel and Schwabl, 1992).
Regarding the latter category of models, Bak et al. (1990) proposed
an interpretation of the spatial distribution of wildfires involving the
theory of Self-Organized Criticality (SOC). This theory was first formalized
by Bak et al. (1987) for the explanation of
Emergent properties, such as the frequency-size scaling, can be revealed searching for power law distributions in empirical data (Burroughs and Tebbens, 2001; Clauset et al., 2009; Deluca and Corral, 2013). Those distributions have been widely found in nature and these occurrences have often been linked to scale invariance and fractals. Power law theoretical background has been applied to different natural hazards (Turcotte and Malamud, 2004), such as earthquakes (Gutenberg and Richter, 1956), fracture systems (Bonnet et al., 2001), landslides (Frattini and Crosta, 2013; Guzzetti et al., 2002; Malamud et al., 2004), rock and snow avalanches (Crosta et al., 2007; Birkeland, 2002), tropical cyclones (Corral et al., 2010) and wildfires (Malamud, 1998). For the latter, the scaling exponent of the power law has been estimated from different data sets in various geographic contexts (Corral et al., 2008; Fiorucci et al., 2008; Lin and Rinaldi, 2009; Malamud, 1998; Ricotta et al., 1999, 2001; Song et al., 2001, 2006). Scaling exponents have been also interpreted as a proxy of fire regimes (Hantson et al., 2015) in North America (Malamud et al., 2005a; Millington et al., 2006) and in the Iberian Peninsula (Moreno et al., 2011). Nevertheless, the existence of a power law scaling in wildfires has been questioned by several authors (Reed and McKelvey, 2002; Benavent-Corai et al., 2007; McKenzie and Kennedy, 2012; Lehsten et al., 2014) and some uncertainties still remain about their frequency distribution. For example, Reed and McKelvey (2002) suggested that it is necessary to have a constant extinguishment-growth rate ratio (i.e. fire suppression and vegetation recovery are constant in space and time) to hold a power law at all sizes of burned areas. Since it is difficult for both extinguishment and growth to be spatially and temporally constant, a power law is not expected to hold for the entire range of fire size. Recently, Lehsten et al. (2014) proposed an alternative statistics, the Characteristic Fire Size (CFS), to measure the fire size that contributes the most to the total burned area. CFS in the pan-Boreal area showed to explain the distribution of wildfires better than a power law distribution, using two independent inventory datasets.
Regarding the fire models implemented in Land Surface Models (LSM), such as the Community Land Model (CLM), different comparison between observed and simulated data has been performed by many authors (Kloster et al., 2010; Migliavacca et al., 2013a; Thonicke et al., 2010). The ability of LSM to mimic realistic wildfires dynamic has been questioned, and further comparison are needed to assess the reliability of LSM prediction for future global climate scenarios.
Recently Luo et al. (2012) proposed a framework methodology to benchmark LSM. Those benchmark analyses aim to: (i) target aspects of model performance to be evaluated, (ii) define a set of benchmarks as references to test model performance, (iii) compare performance skills among models to identify model strengths and deficiencies, and (iv) propose model improvements. As scaling behavior of natural hazard represents an important insight in the complexity of this natural phenomenon, those scaling properties can be used as a benchmark for LSM. At date, only a work by Fletcher et al. (2014) incorporates power law scaling in dynamic global vegetation models (DGVMs), observing promising results for the Amazonian region.
In this paper, we benchmark a one version of the CLM model with power law scaling exponent and Characteristic Fire Size (CFS) calculations using a well-established observational data set of forest and non-forest wildfires (the European Forest Fire Information System, EFFIS) (Camia et al., 2014; McInerney et al., 2013). In particular, one criterion for model benchmarking is to evaluate how model structures are able to reproduce the relationships between a set of particular variables against a set of benchmarks, the latter derived from observed data. For this reason, we evaluated whether a LSM is able to mimic the probability of the distributions of wildfires in Europe that emerges from the analysis of a well-established observational dataset. The LSM here used is a recent parameterization of the CLM for its specific application in Europe (Migliavacca et al., 2013a, b). The comparison is performed because, although the model has been successfully applied to simulate the mean burned area and its relationship with environmental and climatic conditions (Kloster et al., 2010; Migliavacca et al., 2013a), poor performances have been reported in the description of interannual variability and severe fire seasons (Migliavacca et al., 2013a). In particular, in Europe severe fire season characterized by large events (San-Miguel-Ayanz et al., 2013) are source of a great concern in the scientific community because of the direct link with climate extremes (e.g. 2003 and 2010 heat waves).
The European Fire Database contains information on forest and grassland fires compiled by European Union Member States and by other European countries within the framework of the European Forest Fire Information System (EFFIS) (Camia et al., 2014).
EFFIS, the European Commission reference system for pan-European fire
information, is a modular decision support system that monitors fires at
a continental scale supporting fire management across Europe with the
contributions of national forest fire services in the single countries
(McInerney et al., 2013). In the EFFIS dataset,
forest fires are defined as uncontrolled vegetation fires spreading wholly
or in part on forest and/or other wooded land. Grasslands and shrublands are
also included, while agricultural land is excluded from the fire statistics.
We imposed a minimum fire size for data analysis equal to 1
Every year, contributing countries submit fire data derived from the collection of individual fire records, which are routinely compiled by local fire fighters after fires have been extinguished. Data are checked, stored and managed by the Joint Research Centre (JRC) within the European Fire Database of EFFIS.
Among the available data, we used monthly totals of burned area from 2000 to 2009 in order to have a temporal resolution equal to the LSM. The EFFIS data were aggregated at NUTS (Nomenclature of Territorial Units for Statistics, EUROSTAT, version 2006) level 3, which correspond to local administrative units of variable extent.
The model simulations conducted in this study were performed with a modified
version of CLM version 3.5 (e.g. Stöckli et al., 2008) extended with a carbon-nitrogen
biogeochemical model
(Randerson
et al., 2009; Thornton et al., 2007, 2009). The prognostic treatment of
fires is based on the algorithm developed by
Kloster et al. (2010) and modified by
Migliavacca et al. (2013b) for its application in Europe, where
the model has been successfully applied to simulate fires for the present
climate. Briefly, the model simulates burned area as a function of
population density, soil moisture, biomass available to burn and wind-speed.
The runs used in this study refer to the period 2000–2009, and were
conducted at a spatial resolution of
Both EFFIS and CLM data sets were classified using a simplified version of the climatic stratification proposed by Metzger et al. (2005), also known as the “Environmental Stratification of Europe”. This classification is based on a principal component analysis and a cluster analysis of available environmental variables. The original 13 environmental zones were merged into 5 classes (Fig. 1): Continental, Boreal, Atlantic, Alpine and Mediterranean, in order to resemble principal biomes in Europe and to avoid highly fragmented classes that are incompatible with LSM resolution.
Power law scaling exponents are expected to reveal emerging properties of complex systems. In this framework, different eco-climatic classes are expected to show different power law coefficients. Regarding wildfires: higher (absolute) values of the exponents could be interpreted as an indicator of a greater incidence of small monthly burned area with respect to the total burned area; conversely, lower exponents will reveal that the contribution of large extreme events to the total burned area is greater.
Power law exponents are here estimated using a Maximum Likelihood Estimation (MLE) (Clauset et al., 2009), in particular we used recent methodology based on MLE applied to binned data (Virkar and Clauset, 2014).
Power law distributions are characteristic probability density distribution
functions that can be expressed as:
Clauset et al. (2009) proposed a unified method to fit power law
distributions in empirical data, then modified by Virkar and Clauset (2014) for the
application to binned data. Those methods have been largely applied by
researchers interested in power law distribution and it can be summarized as
follows. The estimated power law
Different probability distributions were tested against the power law. In particular, we here tested the lognormal, the exponential, the power law with exponential cut-off and the stretched exponential (Weibull) as possible competitors to the power law. Those functions are listed in Table 1.
The comparison is made by calculating the logarithm of the ratio of the two likelihoods (Loglikelihood Ratio, LR), since it provides a simple test statistic for discriminating power law from other distributions. LR assumes positive values if the first distribution (the power law) better fits the data with respect to another distribution, it is indistinguishable from zero in the event of a tie. The significance of this comparison is then assessed using the method proposed by Vuong (1989).
Data binning represents a delicate topic in power law research. Here we applied a standardized and objective binning scheme based on 20 bins distributed as a logarithm between the minimum and the maximum value of burned area for each class. This method is consistent with previous works that analyzed power law distributions in wildfires (Malamud et al., 2005a; Moreno et al., 2011; Ricotta et al., 1999).
It is important to underline that the method of fitting here applied aims to detect the power law above a minimum bin value. The method assumes that only a minimum threshold exists and only the upper tail of the distribution is fitted by a power law or an alternative distribution.
The code for fitting power laws and testing alternative distribution was iteratively applied to each ecological class both for EFFIS and for CLM data sets. Then, the inverse cumulative distribution functions (CDF) were represented and results were compared.
Recently, Lehsten et al. (2014) proposed a novel statistic to analyze
wildfires size distribution called the Characteristic Fire Size (CFS). The
authors claim that whether the CFS is normally distributed, this feature can
be used as support for rejection the power law distribution of fire sizes.
Here we calculate the CFS with the formula:
We calculated CFS for both CLM and EFFIS data, using the same binning scheme described in Sect. 3. Results are then are compared and CFS are spatially represented at European scale.
First basic comparisons between EFFIS and CLM were made representing the two time series of sum of monthly burned area for the different ecoregions here analyzed (Fig. 2). As already reported using an administrative stratification (Migliavacca et al., 2013a), CLM strongly overestimates burned areas. We found this feature to hold true also for the climatic stratification applied here. In ecoregions characterized by a little fire seasonality (i.e. Alpine, Atlantic and Continental), CLM badly reproduce EFFIS time series. For the Boreal domain, where fires spread mostly in the summer season and are almost absent during the winter and autumn, CLM qualitatively reproduces this trend observed in EFFIS data. Regarding the Mediterranean ecoregion, the strong seasonality of monthly burned area is reproduced by CLM. Further detailed and quantitative comparisons between EFFIS and CLM can be found in Migliavacca et al. (2013b).
The inverse cumulative distribution functions (CDF) of monthly burned areas
are shown in Fig. 3. These results were obtained with the MLE method
(Virkar and Clauset, 2014) and applied respectively to EFFIS and
CLM data aggregated with the simplified Metzger classification
(Metzger et al., 2005). Results for the EFFIS
data set (red lines) showed statistical significance (
The ranges of validity for the PL reach about three orders of magnitude for the Continental biome and about two orders for the Atlantic, Alpine and Boreal ecoregions. For the Mediterranean biome, although statistically significant, only one order of magnitude of burned area is explained by the PL distribution. For the Boreal and Continental biomes we found a deviation from the power law fit in the tail of the distribution, often referred as “cut-off”.
Regarding CLM data (green lines), power law fits are limited to the upper
two bins of the whole distribution, signifying that the power law is not
a suitable model for the description of those data and only the fires larger
than 10
In Table 3 are presented the Log-likelihood Ratio (LR) for the comparison of
PL with other statistical distributions (summarized in Table 1). The
positive sign of LR values means that, in the PL validity range, this
distribution is favored against the others. The significance of the
comparison was assessed by a statistical test (Vuong,
1989) that resulted significant for the majority of the comparisons. In some
cases, the LR comparison was not statistically significant (
Figure 4 shows the CFS plots for each ecoregion in Europe. Regarding EFFIS,
the Alpine, Atlantic and Mediterranean biomes feature a similar behavior,
where the CFS increases for small burned areas, reaches a maximum value and
then decreases again for larger burned areas. For the Mediterranean
ecoregion, data show a Gaussian shape, since it is clearly discernible both
an increasing and a decreasing trend, with a peak at 2069
In Fig. 5, we show a spatial representation of power law exponents for the EFFIS data set at European scale. A latitudinal gradient is displayed in the map, where higher exponents are associated with lower latitudes and lower exponent are generally associated with higher latitudes. The map of power law exponent for the CLM is not showed here since those fitting resulted non statistically sound.
In Fig. 6, a map of CFS for the CLM simulated data is showed. Also in this case, a spatial pattern appears: higher CFS are displayed in the Mediterranean area and lower CFS in the Continental and Boreal areas. The map of CFS is here showed only for CLM since EFFIS data showed that not all the ecoregion allowed the estimation of a CFS.
The main objective of this work was to test whether a power law function approximates wildfires distribution across European biomes using two large data set (EFFIS and CLM). In particular, specific object was to use scaling coefficients to benchmark a land surface model (CLM) against a well-established observational data set (EFFIS). These objects address important questions regarding wildfire theoretical background, since whether fire size distribution follows or not a power law distribution in Europe represents a fundamental insight in the complexity of fire dynamics. We also remark that data aggregation plays a fundamental role in power law scaling analysis (Fiorucci et al., 2008), since the more the single classes of data are coherent, the more emergent behavior are likely to appear (Cristelli et al., 2012).
In order to compare EFFIS and CLM dataset, we applied a temporal (monthly) and spatial (ecoregions) aggregation of fires. This aggregation was necessary since CLM does not simulate single fires events but only temporal aggregated sums of burned area in each pixel, and for this specific dataset the aggregation is monthly. We also assumed that different European ecoregions features different fire dynamics and we aimed to characterize them with power laws, as previously applied by Malamud et al. (2005) in North America.
Results from the EFFIS data set showed that wildfires distribution in the Boreal and Continental ecoregions feature a statistically sound power law for many orders of magnitude. The Alpine, Atlantic and Mediterranean ecoregion show a power law only in the upper tail of the distribution and only for a few orders of magnitude, we considered those fit not significant for the purpose of the paper.
We interpret these results (Figs. 3 and 5) using the Self Organized Criticality (SOC) theoretical framework. In the Boreal and Continental biomes, wildfires feature a general lower incidence with respect the Mediterranean domain (De Groot et al., 2013). In the Boreal and Continental biomes, where mainly evergreen needle leaf and broad leaf forests are present, wildfires are a rare event that can spread also for large areas since the landscape is homogeneous and non-fragmented. Furthermore, those land covers are not very resilient to fire (Goetz et al., 2006; Lentile et al., 2006) and this could act like a feedback mechanism, where fire is prevented by the absence of fuel available to burn due to previous fire. The SOC framework can be applied here since Boreal and Continental ecosystems accumulates (bio)mass (e.g. litter and deadwood) and energy through the seasonal cycles of vegetation by converting the constant input of solar radiation and carbon dioxide. Mass and energy are then dissipated by fire as a fractal, since small wildfires are naturally occurring in large number (Randerson et al., 2012), and when fire spreads it is difficultly suppressed before it can reach large extension because of the small population density and little landscape fragmentation. In any case, very big fire events are prevented through fire suppression by fire-fighters, this may results in a deviation from the power law fit in the tail of the distribution (Fig. 3), often referred as a finite-size effect (Corral et al., 2010; Hergarten, 2002). These results support that Northern Europe ecosystem show a more natural fire regime, since the influence of humans is located mainly in the upper tail of the distribution, and fire ignition and spreading is more likely be dictated by climate and fuel availability.
CFS analysis on EFFIS data shows that Boreal and Continental biomes do not feature a characteristic size of monthly wildfires. In particular, the Boreal biome shows a double peak in CFS plot (Fig. 4), this could mean that small and large fires follow a different behavior, which is instead masked in the cumulative distributions representation (Fig. 3). The Continental biome features no characteristic size, this can be used as a further support to the power law scaling.
Southern and Central Europe (SCE) biomes (i.e. Alpine, Atlantic and Mediterranean) do not feature a statistically sound power law, and the CFS revealed that they show a Gaussian-like distribution featuring one characteristic dimension. SCE biomes are highly fragmented and wildfires are usually triggered and suppressed by human beings (Ganteaume et al., 2013). We claim that the high population density and landscape fragmentation play an important role in the lack of power law fit in SCE biomes. Wildfires are often triggered by arson or accidental causes (roughly in 90 % of the total fires) that involve human being. Furthermore, active fire-fighting methods and efficient fire suppression, strongly limit the spreading of wildfires in different land covers. This results in a strong deviation from a “natural fire regime”, where fire spreading is mainly dictated by climate variability and vegetation spatial distribution. Furthermore, these ecosystems are generally highly resilient to fire (Díaz-Delgado et al., 2003; Di Mauro et al., 2014) and repeated fire can happen in the same burn scar from one year to the subsequent, preventing the feedback mechanism described for the Boreal and Continental biomes. Also in these areas, large fires can occur; those events are possibly escaped from human suppression. For example, events like these can be triggered by extreme heat waves (Barriopedro et al., 2011; Trigo et al., 2006) that create environmental conditions that promote large fires spreading. The CFS analysis showed that for the Mediterranean area a clear characteristic size can be calculated, instead for the Alpine and Atlantic biomes we found some noise in the upper tail of the distribution, that can be related to the incidence of large fires, that is higher than expected in a lognormal distribution. In fact, for those two areas the power law analysis revealed that the scaling holds only in the upper tail with statistical significance, even if only a few bins are fitted by the power law.
Regarding the benchmark analysis of CLM, we used two methodologies: power law fitting and the CFS analysis. From both results we found that in its current formulation, CLM badly reproduce the properties found in the EFFIS data set. In particular, the power law function is not a suitable description of wildfires distribution for all ecoregions in Europe. The fitting is limited to the upper two bins of the distribution, and we accounted it as not enough for a statistically sound fitting. The CFS analysis showed that, using CLM, in all the European ecoregion the mean values of burned area is well behaved, since a characteristic event size can be calculated from the peak of the Gaussian-like distributions represented in Fig. 4.
The model benchmarking strategy applied in this study showed that improvements are needed in using this model for prediction of the distribution of wildfires in Europe. This can be done by developing the model structure, for instance to reduce the overestimation in Central Europe and the lack of description of interannual variability (Migliavacca et al., 2013a). At the moment, fire suppression is parameterized just as a function of population density; we can suggest that, for example, including the scaling parameter in the model structure can help in limiting fire spreading and in representing in a more realistic way the incidence of wildfires in Europe. Another option is to use a model optimization based on multiple constrains, for instance by including the observed power laws by ecoregion to describe better the spatial variability of fire regimes. A pilot work by Fletcher et al. (2014) showed that including the scaling factor in a vegetation model increased the estimation accuracy of burned areas in the Amazonian region.
We also found that spatial representation of alpha exponents calculated from EFFIS data and CFS calculated from CLM clearly depicts differences in fire regime in Europe, as previously found in North America (Malamud et al., 2005a; Millington et al., 2006), since they display a latitudinal gradient related to the incidence of climate on fire regimes. In particular, lower scaling coefficient calculated from EFFIS are related to higher latitudes, instead higher CFS values calculated from CLM are related to lower latitudes.
Power law generating mechanisms play a fundamental role in the interpretation of results (Carlson and Doyle, 1999, 2002; Moritz et al., 2005); the fact that heavy-tailed distributions occur in complex systems implies that extreme events occur more frequently than they would in other distribution (Stumpf and Porter, 2012). This is an important issue in understanding the impact and relative frequency of large wildfires in Europe (San-Miguel-Ayanz et al., 2013).
In this paper, we analyzed wildfire distribution across Europe using EFFIS and CLM data. Regarding EFFIS data, we found that cumulative distribution of Boreal and Continental biomes can be approximated by a power law function with statistical significance. Instead, the Mediterranean, Atlantic and Alpine biomes did not show significant fits. We interpreted these results as a result of a feedback mechanism between human management of wildfires and natural fire regime. We benchmarked CLM simulation with observed data, and we found that the model does not reproduce the scaling found in the Boreal and Continental biomes, and it strongly overestimates burned areas in all biomes. Possible constraints regarding the inclusion of scaling exponents in the model formulation are proposed in order to reduce the spreading of fire using EFFIS observed variability.
Based on theoretical considerations, we can also speculate that, if the power law properly describes wildfire distribution, it would be very unlikely observed at all orders of magnitude in inventory data or simulated monthly burned area that includes anthropogenic influence. Furthermore, since wildfire is a natural hazard with strong human influences both in ignition and suppression, we conclude that a power law scaling behavior is likely observed only for natural environments, and deviation from the power law can be ascribed by the influence of human being on natural fire regimes.
We acknowledge the European Commission, Joint Research Centre (JRC), for providing the dataset extracted from the European Fire Database of EFFIS.
Alternative probability distributions tested against the power law.
Summary of power law fits for the EFFIS data set. Significant
Kolmogorov–Smirnov (KS)
Logarithms of the likelihood ratio (LR) between a power law
distribution and four alternative distribution for the EFFIS data set.
Positive values of the ratio means that the power law is favored over the
alternatives. Significances are expressed for
Climatic classification of European environment (Metzger et al., 2005). The 5 class stratification used here include Continental, Boreal, Atlantic, Alpine and Mediterranean. Shaded areas represent the available EFFIS data set in the considered period (2000–2009).
Monthly burned area time series for different ecoregions across Europe for the EFFIS (orange line) and CLM (blue line) data set. Note that ordinates are represented in a logarithmic scale.
Results of power law fitting for the EFFIS (red lines) and CLM (green lines) data sets for each ecoregion across Europe. In these graphs, the abscissa represents the monthly burned area (in hectares, ha) and the ordinate represents the inverse cumulative distribution function (CDF). Note that both axes are represented in a logarithmic scale.
Characteristic Fire Size (CFS) distribution for the EFFIS (blue lines) and for the CLM (green lines) data. Each plot represents a different ecological region in Europe.
Spatial representation of power law
Spatial representation of the Characteristic Fire size (CFS) for the Community Land Model (CLM) burned area simulations