Lessons Learned from Arson Wildfire Incidence in Reforestations and Natural Stands in Spain

Wildfires are currently considered the major threat to forests in Mediterranean countries. It has been implied that a large percentage of arson-caused fires in Spain are connected with the extensive reforestation programs implemented between 1940 and 1970. However, no consistent studies have been conducted to study the relationships between arson-caused fires and stand origin. Therefore, the goal of this study was to analyze occurrences and model the influence of forest stand origin (artificial or not) on the development of wildfires in peninsular Spain. Twenty-one neural network models were trained to estimate fire incidence through fire type (surface or crown fire), burned area and total treed burned area, based on stand age (years), canopy cover (%), natural age class (from seedling to mature stages) and fuel type classification. Models were built for reforested stands and natural stands of Pinus pinaster Ait., the Mediterranean pines Pinus sylvestris L., Pinus nigra Arn., Pinus halepensis Mill. and Eucalyptus sp. L’Hér., or groups of these species, and the resulting models were compared. Reforested stands presented higher fire incidence than natural stands mainly for productive species like Pinus pinaster Ait. According to the fire type models, thickets had a large influence in the development of crown fires in reforested stands in a general model for all species, the model with the Mediterranean group of pines, and the Pinus pinaster Ait. model. Vertical continuity influenced crown fire propagation in natural Mediterranean pines and in Eucalyptus stands. Presence of shrubs, grasslands and wood slash was related to surface fires in models for both reforested and natural stands. The results suggested that stand origin was influential on fire incidence, at least with regard to fire type and commercial species in the northwestern region of Spain.


Introduction
During the last decades, more than 95% of the wildfires occurred in Spain and in other Southern European Mediterranean countries were of anthropogenic origin, and over 55% of them were classified as arson-caused fires [1,2].Arson is not random, and many reasons have been put forward to explain why arsonists burn forests [3].Among them, Simon et al. [4], Alvarez et al. [5] and Fuentes-Santos et al. [6] have cited resentment against reforestations as a motivation for arson-caused fires.A relationship between arson fires and reforestations is assumed in the Spanish Wildland Fire National Statistics of the Ministry of Agriculture, Fishing, Food and Environment, which includes in the fire reports a code for "Arson fires motivated by resentment against reforestations" [1].In the past, these arson-caused fires would have been triggered by the discontent of many rural communities for the shifting of their productive and rangelands into forest [7] and towards the wood rights' owner (and thus the beneficiary of the wood's sale) when land and wood differed in ownership [8].However, other motivations have risen over time; forcing land use changes is a potentially relevant, but difficult to prove, motivation in Southern European countries that depend on the tertiary sector [9].Regardless of underlying motivation, the relations between arson wildfires and reforestations have received little attention in the literature.Previous studies on general human-caused fires, not specific for arson, have sometimes included forest-related variables like commercial value [10] or species composition [11][12][13][14][15]. Forest stand characteristics are usually accounted for and converted to fuel type categories [16].Yet, very little scientific evidence exists on how stand origin may affect other drivers and influence arson fire incidence, despite decades-long controversies [8,17,18].
Reforestation has been frequent in many European countries in which clear-cutting [19] and natural disturbances such as windstorms [20] or wildfires [21] have affected large tracks of forest.In addition, socioeconomic drivers from the early 1940s onwards have prompted European Mediterranean countries to implement extensive governmental reforestation programs (i.e., 3 million ha in Spain during 1940-1970, 0.8 million ha after 1952 in Italy, 340,000 ha in Portugal in 1987-1999 [22]).The natural or artificial origin of stands, and later, the effects of natural and anthropogenic disturbances determine forest characteristics, namely, composition and structure.Stand composition and structure are supposed to differ widely between naturally expanding, colonizing stands, and artificially reforested stands that follow different ecological processes [23][24][25].Forestry agencies have promoted monospecific and even-aged stands, selecting species for site restoration based on their fast growing and ground cover characteristics [26].Consequently, reforested forest stands generally have higher even-age tree density and lesser variation of diameter at breast height (DBH), quantity of seedlings, and tree and genetic diversity [18,27].
In Southern European Mediterranean environments, stands arising from secondary succession often share characteristics with reforestations.Land abandonment and wildfires have contributed in Spain to a pervasive and fast expansion of shade-intolerant conifer species [28,29] like Aleppo pine (Pinus halepensis Mill.), also used in reforestation for their pioneer capacity.Several authors have stated that wildfires in South Europe have affected mainly conifers, with Pinus pinaster Ait., Pinus halepensis Mill., and Eucalyptus sp.L´Hér.stands reforested in the 1940's and 1950's being comparatively more affected than other species [11][12][13][14][15].However, in terms of forest structure, wildfires tend to favor high-density even-aged forest stands [27,29], areas of low average tree diameter [12], or undermanaged or abandoned stands, which tend to develop vertical continuity of fuels [30]; it may be argued that both types of stands, either reforested or naturally arising from secondary-succession after abandonment since the 1950s, currently share these features and are both of them potentially highly susceptible to fire [31,32].
Our work aimed to compare the past incidence (occurrence and behavior) of arson-caused fires in reforested and nonreforested fire-prone stands in Spain.We explored fire incidence patterns to discredit or support the notion that reforestations are comparatively more affected by arson, due to resentment or other motivations.Similarity of fire incidence patterns for any stand origin would indicate arson is not specific to reforested stands, and that forest characteristics are alike across natural and artificial stands of fire-prone species determining fire behavior.Dissimilarity of patterns by origin could indicate differences in stand characteristics (behavior) and/or a differential pressure from arson (occurrence), requiring further analysis.
The hypothesis of our study was that stand origin influences ignitions, forest composition and stand structure to produce different patterns in fire occurrence, fire type and burned area.We analyzed fire occurrence, fire type and burned area in reforested and nonreforested stands, based on species composition, natural class age, canopy cover and fuel type.

Materials and Methods
The study area encompassed the Spanish mainland, covering around 492,174 km 2 (Figure 1).The predominant Mediterranean climate in Spain, characterized by long summers of high temperatures and low rainfall, contributes to the high risk of fires in the country [33].The fuel environment has undergone critical changes since the 1950's, with a general increase in continuity and load due to land abandonment and socioeconomic changes, including reforestation [34].Changes in the futures of fire regimes have been detected in past decades, with noteworthy shifts after the late 1980's and the first half of the 1990's [35].After the extreme fire events in 1998 in Spain, wildfire policies and management (Law 43/2003, of November 21th) moved to reducing the gap between suppression and prevention budgets in the country, with a temporary increase in forest management interventions that was short-lived.The financial crisis that hit Spain in 2007-2008 [36] undermined the forestry sector, leading to the closing of many industries and, again, to pervasive forest abandonment.However, the crisis also produced a reduced pressure on the land resource for development and construction [37].This factor may have influenced the occurrence arson wildfires, which decreased from an average of 10,587 arson-fires/year in 1998-2008 to 7336 fires/year in 2009-2015 [1].The economic recovery of the country is currently leading to fire environmental conditions very similar to the conditions in the pre-financial-crisis period 1998-2008 [1].Analysis of past patterns may provide clues for exploring future relations between arson fires and forest stand characteristics, determined by their origin, which justifies the selection of 1998-2008 as our study period.

Reforestation Historical Data
About 54% of Spain is covered by woodlands and shrublands, but only 35% can be considered forest [38].Reforestations (Figure 1) represent almost 8% of the total land area of Spain [39].

Materials and Methods
The study area encompassed the Spanish mainland, covering around 492,174 km² (Figure 1).The predominant Mediterranean climate in Spain, characterized by long summers of high temperatures and low rainfall, contributes to the high risk of fires in the country [33].The fuel environment has undergone critical changes since the 1950's, with a general increase in continuity and load due to land abandonment and socioeconomic changes, including reforestation [34].Changes in the futures of fire regimes have been detected in past decades, with noteworthy shifts after the late 1980's and the first half of the 1990's [35].After the extreme fire events in 1998 in Spain, wildfire policies and management (Law 43/2003, of November 21th) moved to reducing the gap between suppression and prevention budgets in the country, with a temporary increase in forest management interventions that was short-lived.The financial crisis that hit Spain in 2007-2008 [36] undermined the forestry sector, leading to the closing of many industries and, again, to pervasive forest abandonment.However, the crisis also produced a reduced pressure on the land resource for development and construction [37].This factor may have influenced the occurrence arson wildfires, which decreased from an average of 10,587 arson-fires/year in 1998-2008 to 7336 fires/year in 2009-2015 [1].The economic recovery of the country is currently leading to fire environmental conditions very similar to the conditions in the pre-financial-crisis period 1998-2008 [1].Analysis of past patterns may provide clues for exploring future relations between arson fires and forest stand characteristics, determined by their origin, which justifies the selection of 1998-2008 as our study period.

Fire Historical Data
The fire registry was obtained from the General Fire Statistics of Spain (EGIF) of the Coordination Centre for National Information on Wildfires (CCINIF) of the Ministry of Agriculture and Fishing, Food and Environment [1].The total amount of fires recorded during 1998-2008 was 417,406 but only arson fires were considered for the study (130,405 fires).Each fire record provided information on fire type, burned area, total treed-area burned, cause, motivation, species affected, stand mean age, natural age class, percentage canopy cover of the affected vegetation, fuel type, and many other variables related to fire suppression (detection date, arrival time, suppression resources, etc.) and losses caused (i.e.estimated monetary value of the forest products lost).
We selected several dependent (Fire Type-FT: surface/crown fire; Burned Area-BA in ha; and Treed Burned Area-TBA in ha) and independent variables (Table 1) from the EGIF database in agreement with the stated hypothesis.Mean stand age (AGE) was available in the fire reports, but not density, so we selected the percentage of canopy cover (CC) as independent variable.Stand structure was defined through five mutually exclusive natural age classes, from seedlings to mature forest (i.e., seedling-SD, sapling-SS, thicket-TH, pole-PS and high-forest-HFS).The fuel type (FTP1-15) classification used by the EGIF database was also included to compare the relative relevance of the fuel type class versus other forest stand characteristics on the incidence outcome.The EGIF fuel type classification is an adaptation of Rothermel's [42] conducted by the Ministry of Agriculture and Fishing, Food and Environment.

Database for Model Building
Numerical and spatial analysis techniques applied with a dedicated GIS allowed matching ignition points to forest stands by origin.The number of fires located within reforestation polygons (RF) was 21,823 (Figure 2), while 108,582 fires started in nonreforested areas (NRF) in the study period.We removed missing data records and selected only those in which the Forest Map of Spain data was in agreement with the fire report in terms of forest characteristics, which considerably reduced the available data.The available data from 130,405 to 8022 usable records.We sampled arson fires in reforested stands by natural class age category to cover the whole range of characteristics, leaving 4810 RF observations and 3212 NRF observations in the database for analysis.The number of arson fire ignitions in RF and NRF stands were compared by species and weighted by area.
characteristics, leaving 4810 RF observations and 3212 NRF observations in the database for analysis.The number of arson fire ignitions in RF and NRF stands were compared by species and weighted by area.Additionally, RF and NRF databases were split by species before fire behavior modeling was carried out.Fire behavior features (Fire Type-FT, Burned Area-BA, and Treed Burned Area-TBA) were modelled separately for RF and NRF stands of the species listed in Section 2.1.RF and NRF P. pinaster and RF stands of Eucalyptus sp. were also analyzed in separate models.No arson fires in NRF stands of Eucalyptus sp. did exist, as this non-native species in Spain is always planted.The pine species Pi. halepensis, P. nigra and P. sylvestris were grouped for model building, based on their geographical distribution allowing mixing in the fringes of contiguous altitudinal levels in continental Mediterranean mountains, and because of a lower number of records available than for P. pinaster, the species most used in reforestation in Spain.Because the number of fire records of the other species used in reforestation (P.pinea, P. radiata, Populus, P. canariensis and P. uncinata) was likewise not large enough to build suitable models, they were also discarded for individual modeling.These species have been used in reforestation less frequently (<7%), but are affected by fire, so they were included in the global models in RF and NRF stands.
The number of fire observations varied for each modeling database (Table 2).The largest data set was for modeling Fire Type (surface/crown) in Reforested Stands (1524 cases) and the smallest was the dataset for modeling Treed Burned Area in Reforested pine stands (40).Databases were as much balanced as possible for numbers of observations in each category for classification.

The Modeling Approach: Artificial Neural Network (ANN) Models
Models for FT, BA and TBA were built based on any combination of the independent variables in Table 1, for the different species or groups of species (all species, continental Mediterranean mountain pine species, P. pinaster, Eucalyptus sp.), and forest stand origin (reforested/nonreforested).A total of 21 models (Table 2) were built.Additionally, RF and NRF databases were split by species before fire behavior modeling was carried out.Fire behavior features (Fire Type-FT, Burned Area-BA, and Treed Burned Area-TBA) were modelled separately for RF and NRF stands of the species listed in Section 2.1.RF and NRF P. pinaster and RF stands of Eucalyptus sp. were also analyzed in separate models.No arson fires in NRF stands of Eucalyptus sp. did exist, as this non-native species in Spain is always planted.The pine species Pi. halepensis, P. nigra and P. sylvestris were grouped for model building, based on their geographical distribution allowing mixing in the fringes of contiguous altitudinal levels in continental Mediterranean mountains, and because of a lower number of records available than for P. pinaster, the species most used in reforestation in Spain.Because the number of fire records of the other species used in reforestation (P.pinea, P. radiata, Populus, P. canariensis and P. uncinata) was likewise not large enough to build suitable models, they were also discarded for individual modeling.These species have been used in reforestation less frequently (<7%), but are affected by fire, so they were included in the global models in RF and NRF stands.
The number of fire observations varied for each modeling database (Table 2).The largest data set was for modeling Fire Type (surface/crown) in Reforested Stands (1524 cases) and the smallest was the dataset for modeling Treed Burned Area in Reforested pine stands (40).Databases were as much balanced as possible for numbers of observations in each category for classification.

The Modeling Approach: Artificial Neural Network (ANN) Models
Models for FT, BA and TBA were built based on any combination of the independent variables in Table 1, for the different species or groups of species (all species, continental Mediterranean mountain pine species, P. pinaster, Eucalyptus sp.), and forest stand origin (reforested/nonreforested).A total of 21 models (Table 2) were built.Artificial neural networks (ANN) were chosen to explore the relationship between stand characteristics and fire behavior because they had proven useful in other wildfire studies, providing good accuracies and generalization capability (i.e., [29,43,44]).When dealing with non-normal distributions of the data, or nonlinear, complex relationships, these models are as reliable as any traditional statistical method, and often more robust [45].An ANN is an information processing algorithm capable of learning patterns in data by iterative adjustment of the connections between processing elements (nodes) organized in layers.Our models were feed-forward multilayered non-linear fully connected cascade-correlation networks [46], built using NeuralWorks Predict ® v.3.24software [47].The models were computed as in Alcázar et al. [48] and Debouk et al. [49].The cascade-correlation algorithm starts model building with the simplest possible architecture, adding nodes until an optimal multi-layer structure is reached [46].Training was based on an adaptative gradient learning rule, a variant of the general algorithm of back-propagation [50].
The FT models were based on a binary classification (with one-of-N encoding output either Surface Fire (1,0) or Crown Fire (0,1)).The BA and TBA models were predictive providing one output in hectares.Categorical independent variables were input to the network by using the same one-of-N encoding [29].Interval variables (Age, CC) underwent several transformations that were tested with a genetic algorithm before the actual model training.As usual in ANN learning, each database was divided in three subsets: 10% of the observations were separated to validate the models with data not used in model building.From the rest, 70% were used for training and 30% for testing the accuracy of the model performance and to stop training [51].Several replicates were simulated for each model.Convergence to a same or similar architecture between the replicates was checked to corroborate robustness of the models.Parsimonious models were favored because of the relatively low number of observations, which generally was below 1000 [48].
The best models were selected according to their capability to explain the dependent variable and performance in terms of total percentage correctly classified and confusion matrix (FT classification networks) and R correlation values between observed and predicted values (BA and TBA prediction networks).
A sensitivity analysis based on the partial derivatives of the outputs with respect to the inputs was run to establish which independent variable had the highest impact on the predicted variable for every model.A higher frequency of selection in the variable selection genetic algorithm was also indicative of the greater relevance of an independent variable.

Results
In terms of fire occurrence, the analysis of ignitions started in RF or NRF stands in the study period is presented in Table 3.It is remarkable that the ignition rate of P. pinaster is much higher in RF stands than in NRF stands, in a contrary pattern to the other Pinus species, with higher rates in natural stands.Planted species Eucaliptus and Populus also exhibited high rates of ignition, over one per 100 ha.[52] and the II National Forest Inventory [53]).

Species
Total RF Surface (ha) As for fire behavior dependent variables, unfortunately none of the prediction models developed for Burned Area and Treed Burned Area complied with the correlation requirements set (R > 0.50) or gave satisfactory enough results to reach significant conclusions.Consequently, only the best classification networks built for Fire Type classification are presented here.

Models RF_FT and NRF_FT: Fire Type in Reforested and Natural Stands (All Species)
According to Table 4, RF_FT and NFR_FT models present a simple architecture, with no nodes in the hidden layer in the RF_FT model (6-0-2) and only one in the NRF_FT model (8-1-2).The RF_FT model had a higher accuracy (around 0.85 in the training, testing and validation samples), whereas the NRF_FT accuracy model was around 0.7.In both models, the accuracy for crown fires was higher than that for surface fires, which was below 0.5 in the training, testing and validation samples for the nonreforested stands model NRF_FT.The sensitivity analysis (Table 5) shows that FTP2 (shrubs) and AGE (Stand mean age) variables were selected in both models, being positively correlated to surface fires.FTP2 is the variable with higher weight.Its frequency of selection is very high in both models (with a percentage of 0.94 and 0.75, respectively).The best RF_FT model also shows that a higher proportion of TH (thickets) and PS (pole-stage forests) increases crown fire hazard, whereas AGE, FTP2 (shrubs) and FTP6 (forest with grasses) favor surface fires in reforested stands.
In the NRF_FT model for natural stands, only CC (canopy cover) is related with the development of crown fires.FTP1 (grasses), FTP2 (shrubs) and FTP4 (wood slash) were the most important variables for the development of surface fires in natural stands.RFPINE_FT and NFRPINE_FT models present a simple architecture (Table 6), with no nodes in the hidden layers.The RFPINE_FT model accuracy was 0.82 and the accuracy of the NFRPINE_FT model 0.67.In both models, crown fires accuracy was higher than surface fires accuracy, especially in the RFPINE_FT model.Table 7 shows that FTP2 (Shrubs) was selected in the best RFPINE_FT and NRFPINE_FT as in the general (all species) models, being positively associated to surface fires.CC (canopy cover), FTP11 (forests with shrubs and grasses) and FTP15 (forests with shrubs, grasses and slash) were also selected in both models, but were positively correlated with crown fires.
According to Table 7, the best RFPINE_FT model shows that the development of crown fires in reforested pine stands is higher in those forests with a higher AGE, CC (canopy cover), TH (thicket stage), FTP11 (forests with shrubs and grasses) and FTP15 (forests with shrubs, grasses and slash).On the contrary, a higher probability of surface fires in RF pine stands is associated to a higher percentage of FTP2 (shrubs) and FTP4 (wood slash).The variables TH (1), AGE (0.92), and CC (0.90) had the higher frequency of selection by the genetic algorithm.
The highest probability of crown fires in the best NRFPINE_FT model is mainly related to a higher percentage of CC, FTP8 (bushes and forest), FTP11 (forests with shrubs and grasses) and FTP15 (forests with shrubs, grasses and slash).Surface fires are mainly related to FTP2 (shrubs).FTP2, FTP8, TH and CC are the variables most frequently selected for these models by the genetic algorithm applied to variable selection.RFPP_FT and NFRPP_FT models present a simple architecture (Table 8), with no nodes in the hidden layers (9-0-2 and 7-0-2, respectively).The accuracy in both models followed the same pattern as in previous results (but better than the accuracies of previous models, Tables 4 and 6), and were better for artificial stands (accuracy around 0.88 in RF and 0.80 in NRF stands) and crown fires (0.98-0.90 in crown fires and 0.80-0.70 in surface fires).
Table 9 shows that FTP2 was selected in the best RF and NRF P. pinaster models as in the previous models, also being associated to surface fires.The RFPP_FT model shows that the development of crown fires in P. pinaster stands is higher in those forests with higher percentage of TH (thicket), SS (sapling) and PS (pole) stages and the development of surface fires is more likely with a higher proportion of CC, SD (seedlings) and FTP2 (Shrubs).SD (0.86) and CC (0.64) are the variables with the highest frequency of selection in the RFPP_TF model.The NFRPP_FT model identified high probability of surface fires (or low of crown fires) with higher presence of HFS (high-forest, mature stage), FTP1 (grasses), FTP2 (shrubs), FTP5 (grasses and shrubs) and FTP9 (shrubs and wood slash), with HFS (1) and FTP2 (0.92) being the most often selected variables.

Models RFE_FT: Eucalyptus Artificial Stands
The best RFE_FT model presents a simple architecture, with six inputs, no hidden nodes, and two outputs (Table 10).The best model RFE accuracies were over 0.84 for the training, testing and validation subsets, with higher crown fires accuracy (over 0.92) than surface fires accuracy (over 0.80).The RFE_FT model (Table 11) shows that crown fires in Eucalyptus stands are directly related to FTP8 (shrubs and forest), FTP3 (forest) and FTP6 (grasses and forest).Surface fires are positively linked mainly to CC and FTP9 (shrubs and slash).In the RFE_FT variable selection, FTP3 (0.91), FTP8 (0.81) and AGE (0.8) were most often selected.

Discussion
This study focused on the past incidence of arson wildfires in artificial and natural stands, within a time frame providing a fire environment that is very likely to reproduce in the next few years in Spain.In order to learn our lessons, ignition points obtained from fire history records between 1998 and 2008 were overlapped with the National Forest Map of Spain at that time, to better understand occurrence and fire behaviour patterns in this study period.

Wildfire Incidence: Occurrence in Reforested and Nonreforested Stands
When considering the number of ignition points in relation to the origin of forest stands (reforested or natural) and species composition, we found that P. pinaster indeed presented relatively more arson ignitions in reforested than in nonreforested stands than expected based on its respective area of distribution.On the contrary, proportionally more arson-caused forest fires occurred in natural stands of P. halepensis, P. nigra.and P. pinea than in their reforested stands.P. sylvestris and P. uncinata presented a similar proportion of arson fire ignitions in stands of any origin.
Though some authors stated in previous work that wildfires in South Europe have affected mainly conifers reforested in the 1940's and 1950's (P.pinaster, P. halepensis, and Eucalyptus sp.stands) [11][12][13][14][15], and others [6][7][8] indicated reforestations as a major factor of occurrence of arson wildfires in Spain, we found this only true for a few of the most productive and fast-growing species.Reforestations per se were not comparatively more affected by arson unless we considered regional location.The spatial distribution of arson fires (Figure 2 and reference [1]) reveals that the highest percentage of these was located in the Northwestern (NW) region of Spain, where there was a differential pressure of arson on stands of different origin.In NW Spain, P. pinaster and Eucalyptus sp. are the most common species for reforestation, based on commercial value, and human-caused fire occurrence is high [6].According to Cabana Iglesia [8], NW Spain has seen many occurrences of arson as a form of protest by rural communities in favor of the return of communal forests to the traditional silvopastoral use.Consequently, if this was the underlying motivation, arson occurrence may not have been species-specific, but location-specific, linked to wildlands ownership and land use issues.
On the other hand, the Mediterranean region, especially near the Eastern coastal mountains, has the highest percentage of monospecific and mixed Aleppo pine forest in the country [18], but a large part of the area can be considered as natural forests (from secondary succession), and arson seems to prefer them to reforested stands.Aleppo pine stands in the coastal Mediterranean region are more affected by human-caused fires [1] than Aleppo pine stands in the hinterland, so again, location-specific drivers could explain occurrence patterns.In the Mediterranean, drivers were probably linked to development demands for land, but the general difficulty in the certain identification of motivations and causes for arson makes risky to make assumptions, and precludes analysis beyond the scope of this work.
Thus, other risk factors related to human-caused fire occurrence, like low drought vegetation moisture [16], weather indices [54], traditional agricultural practices [15] or presence of long-term serial fire-setters as suggested by Fuentes-Santos et al. [6] could be better variables for explaining occurrence patterns than stand´s origin or composition.

Wildfire Incidence: Neural Network Models for Fire Behaviour
Modelling the relationship between fire type and a set of different independent variables was successfully achieved by means of cascade-correlation artificial neural networks.The neural networks were parsimonious, robust, and converged in all cases to same/similar solutions when the observations were randomly interchanged in the training, test and validation subsets, and when the initial random weights were shifted [45].However, no robust and reliable models were found for total burned area and treed burned area.This was probably the consequence of having one single set of forest structure parameters for all the burned area in one single fire.There may be a high diversity in forest structures typically affected by any single fire in Spain, but only one description is input in the fire reports, which prevented the successful development of these models.
The fire type models did allow us to explore stand-origin influences; models for artificial and natural stands were different in accuracy, structure, and variables selected.Tables 4, 6, 8 and 10 showed that crown fires exhibited better classification accuracy than surface fires in all the databases.
Out of all datasets and corresponding models, the best accuracy was obtained for the P. pinaster models, the species most extensively used for reforestation (24% of the reforested area), especially in the NW Spain, and the most affected by human-caused fires.Artificial seedlings of this species experimented surface fires, but thickets, young forest and pole-forest were linked to crown fires.Closed canopy cover seemed to favor surface fires [29], a trend contrary to the Mediterranean pines group [25], but understandable by the higher foliage density and thickness of this pine needles and their shadowing effects.Surface fires were related to the presence of shrubs, grasses, and slash in the under-story of high-forest in natural maritime pine stands.
The analysis performed on the continental Mediterranean pines indicated an opposite trend with regards to stand age.Crown fires were more likely to occur in pine reforestations when the stand age increased with closed canopy cover and vertical fuel continuity.Natural stands showed similar trends for canopy cover and verticality in fuels, but thickets were not selected as a crown fire factor, neither was age an input to the model.This may be explained by the fact that more natural conditions in the establishment of these stands mean that age is not coetaneous, competition among individuals is lower and thickets are not so homogeneous spatially [55].
Eucalyptus planted stands appeared to be more prone to crown fires when vertical continuity was present (forest, or forest with grasses and shrubs) but also a somewhat lower canopy cover.Surface fires seemed more likely with denser canopy cover, increasing age of the stands, and presence of shrubs and slash.In any case, presence of shrubs and slash is quite common in Eucalyptus stands, as the combination of abundant understory vegetation due to the Atlantic geographic location of the Eucalyptus plantations often combined with woody debris coming from silvicultural activities, frequently applied due to importance of the pulp market in the area (Galicia).When a flammable young tree layer is added over this understory, crown fires result [56].
Since different variables were selected in each model, the formulated hypothesis stated that the origin of a forest stand influences species, age, density and structure, and in turn fire type was upheld.According to the obtained results, fire type was influenced by different variables in stands of different origin and species.Thickets were influential on crown fires in reforestations within the models for all species.Crown fire probability on a forest stand also increased in previous work with smaller tree diameters, larger basal areas and an increasing variation in tree diameter, which suggest thicket structures with high densities [12,14,57].According to Moreno et al. [58], vertical continuity between understory vegetation (with or without grasses) and trees in pine and eucalypt stands was also linked to crown fires.The variables related to stand age, canopy cover, stand structure, and other fuel types showed variable patterns, but shrubs, grasses and wood slash proved influential for surface fires in all the models.Similar results were found by [59,60], in which fire type in mature stands was linked to the quantity and availability of fuel in their understories.The association of shrubs and surface fires, which was detected in most of the models, is a well-known fact in Spain, and is the basis for the fuel treatments usually applied for hazard reduction (i.e., [61]).General trends for all forest types and origins indicate that the presence of grasses and slash, without brushes and thickets, also favors surface fires (Tables 4, 6 and 8).
Management implications useful for silvicultural treatments and for supporting restoration strategies and policies could be derived from the historical analysis of arson incidence, reinforcing some recommendations by other authors (e.g., [62]) and introducing other actions that may be quite relevant if future conditions evolve as expected in the socioeconomic and ecological environment of Spain.Fuel management is technically well developed in European Mediterranean countries and in Spain, and abundant literature exists on the application of prescribed burning, mechanical thinning or mastication [63,64], but treatments and effects are usually analyzed for individual species (i.e., P. nigra, Piqué and Domenech [65]; P. halepensis, Palmero-Iniesta [66]; P. pinaster, Molina et al. [67] without considering stand origin.In agreement with these findings, aged, more mature forest structures and closed canopy cover generally lead to less severe fires in our models, even in Eucalyptus or P. pinaster stands, but not in P. halepensis, P. nigra and P. sylvestris forests if vertical fuel connectivity was present.However, our results indicate that controlling understory vegetation influences the development of surface fires especially for shrubs in artificial stands, so treatments should be prioritized in them.Previous work has often grouped several pine species [63], but we found that beyond species, also origin makes a difference in fire behavior.Thickets in artificial stands are the most hazardous structures, especially for our grouped pine species, but this does not seem to be the case for natural stands.Consequently, at the landscape planning scale, thinning in thicket structures should be prioritized for reforested over natural stands. Dissimilarity of models for stands of different origin (fire type explained by different explanatory variables) seemed to indicate differences in stand characteristics ruled fire behavior.However, and mainly in terms of occurrence, a differential pressure from arson can also explain incidence, so further analysis should be conducted in this regard.

Conclusions
This study aimed to learn from the historical (1998)(1999)(2000)(2001)(2002)(2003)(2004)(2005)(2006)(2007)(2008) relationship between arson fire incidence (occurrence and behavior) and reforestations in Spain, to anticipate similar future conditions in the fire environment.We intended to discredit or support the notion that reforestations were comparatively more affected by arson.The hypothesis of our study was that stand origin influenced ignitions, forest composition and stand structure to produce different patterns in fire occurrence, fire type and burned area.We analyzed fire occurrence, fire type and burned area in reforested and nonreforested stands, based on species composition, natural class age, canopy cover, and fuel type.
Arson-caused ignitions had a higher rate of occurrence in reforestations of P. pinaster, than in its natural stands, but this preference applied only to these commercial species located NW Spain, not to all considered, everywhere.
According to ANN models built for fire type, thickets had a large influence in the development of crown fires in reforested stands in a general model for all species, the model with the Mediterranean group of pines, and the P. pinaster model.Vertical continuity influenced crown fire propagation in natural Mediterranean pines and Eucalyptus stands.Presence of shrubs, grasses, and wood slash was related to surface fires in models for both reforested and natural stands.The results from this study could not confirm a direct relationship between arson fires and reforestations.However, overall outcomes suggested that stand origin and composition were influential on fire incidence, at least in regards to number of ignitions and fire type.
Future studies should consider additional variables related to location, such as ownership, land uses, topography or weather to improve the accuracy and reliability of models.Understanding arson fire incidence can provide valuable information to support restoration strategies and policies, and also

Figure 1 .
Figure 1.Map of reforested areas in Spain.(Source: adapted from Ruiz de la Torre, 1990) Figure 1.Map of reforested areas in Spain.(Source: adapted from Ruiz de la Torre, 1990)

Figure 2 .
Figure 2. Frequency of ignition points within reforestations in Spain.

Figure 2 .
Figure 2. Frequency of ignition points within reforestations in Spain.

Table 1 .
Description of the 22 independent variables.The description of the fuel types is according to the EGIF database.

Table 2 .
Database definition and number of fire observations used for model building.

Table 3 .
Comparison between the number of ignition points located in reforested and natural stand by species (1998-2008) (source: adapted from Castroviejo et al.

Table 4 .
Best model architecture and confusion matrix for the classification of Fire Type within reforested and not reforested stands (all species).

Table 5 .
Sensitivity analysis of the best models for the classification of Fire Type within reforested and not reforested stands (all species).Models RFPINE_FT and NRFPINE_FT: Fire Type in Continental Mediterranean Mountain Pines in Reforested and in Natural Stands (P.halepensis, P. nigra, P. sylvestris)

Table 6 .
Best model architecture and confusion matrix for the classification of Fire Type for continental Mediterranean mountain pines in reforested and natural stands.

Table 7 .
Sensitivity analysis of the best models for the classification of Fire Type within reforested and not reforested stands (Pine Group).

variables in the best RFPINE_FT model
3.3.Models RFPP_FT and NRFPP_FT: Fire Type in P. pinaster in Reforested and in Natural Stands

Table 8 .
Best model architecture and confusion matrix for the classification of Fire Type within reforested and not reforested stands of P. pinaster.

Table 9 .
Sensitivity analysis of the best models for the classification of Fire Type within reforested and not reforested stands (P.pinaster).

Table 10 .
Best model architecture and confusion matrix for the classification of Fire Type in Eucalyptus sp.

Table 11 .
Sensitivity analysis of the best model for the classification of Fire Type within reforested stands of Eucalyptus sp.