Research Article
Modelling spatial patterns and temporal trends of wildfires in Galicia (NW Spain)
Departamento de Fundamentos da Análise Económica, Facultad de C. Económicas de U. Santiago de Compostela, Santiago de Compostela, Spain
Departamento de Fundamentos da Análise Económica, Facultad de C. Económicas de U. Santiago de Compostela, Santiago de Compostela, Spain
Abstract Aim of study: The goal of this paper is to analyse the importance of the main contributing factors to the occurrence of wildfires. Area of study: We employ data from the region of Galicia during 20012010; although the similarities shared between this area and other rural areas may allow extrapolation of the present results. Material and Methods: The spatial dependence is analysed by using the Moran’s I and LISA statistics. We also conduct an econometric analysis modelling both, the number of fires and the relative size of afflicted woodland area as dependent variables, which depend on the climatic, land cover variables, and socioeconomic characteristics of the affected areas. Fixed effects and random effect models are estimated in order to control for the heterogeneity between the Forest Districts in Galicia. Main results: Moran’s I and LISA statistics show that there is spatial dependence in the occurrence of Galician wildfires. Econometrics models show that climatology, socioeconomic variables, and temporal trends are also important to study both, the number of wildfires and the burnedforest ratio. Research highlights: We conclude that in addition to direct forest actions, other agricultural or social public plans, can help to reduce wildfires in rural areas or wildlandurban areas. Based on these conclusions, a number of guidelines are provided that may foster the development of better forest management policies in order to reduce the occurrence of wildfires. Keywords: Causeeffect relationship; climatology; spatial and temporal indicators; fixed effects; random effects; socioeconomic factors. Citation: Barreal, J., Loureiro, M.L., (2015). Modelling Spatial Patterns and Temporal Trends of Wildfires in Galicia (NW Spain). Forest Systems, Volume 24, Issue 2, e022, 18 pages. http://dx.doi.org/10.5424/fs/201524205713. Received: 06 Feb 14. Accepted: 09 Apr 2015 Copyright © 2015 INIA. This is an open access article distributed under the Creative Commons Attribution License (CC by 3.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. Funding: The author(s) received funding for this work. Finantial support from Project EUI200803685 “Gestión de incendios para el mantenimiento de la biodiversidad y mitigación (Biodiversa CAll1: FIREMAN) (ERANET) “ is acknowledged. Competing interests: The authors have declared that no competing interests exist. Correspondence should be addressed to Maria L. Loureiro: maria.loureiro@usc.es 
CONTENTS 
IntroductionTop
It is estimated that more than 1.3 million hectares of forest are destroyed by wildfires in Europe each year (FOREST EUROPE, UNECE & FAO, 2011). Spain is one of the five southern European countries with the highest level of damage caused by wildfires, with a yearly average of 19,705 wildfires from 1998 to 2007, affecting a total of 130,714 hectares (SECF, 2010). Within Spain, the case of Galicia is particularly relevant. While only representing 6% of national surface area, between 1991 and 2010 Galicia registered an approximate average of 46% of Spanish wildfires and 21% of the total burned surface area (Figure 1), according to MARM (2012) and the regional government (Xunta de Galicia, 2011). Given the geographical concentration of this problem, we limit our analysis to the wildfires occurring in this region, also due to the lack of comparable data for other Spanish regions. We believe, however, that the current paper may provide insights which are closely applicable to other European rural and wild landurban areas.
Figure 1. Galician wildfires with respect to the total number of Spanish wildfires (19912010).

According to data provided by the Galician Institute of Statistics (IGE, 2012a), since 2001 the number of wildfires reveals an upward trend until 2005, decreasing then in number. With regards to the affected surface area, this increased gradually from 2002 to 2006, but since then it has decreased considerably as well. It should be noted that the number of affected areas reach catastrophic levels during 2006. Furthermore, the evolution of wildfires throughout Galicia varies considerably in spatial and temporal terms. Geographically, and based on the data published by the IGE (2012a), wildfires affect more severely southern districts than northern districts of Galicia, both in terms of the number of fires and forestry area affected (Figure 2). Also, western districts are the area in which forestry lands are the most affected in relation to their surface area, while the southern districts record the highest numbers of wildfires.
Figure 2. Geographic distribution of the occurrence of wildfires during 20012010. a) Representation of the number of wildfires per district with respect to the total number of Galician wildfires. b) Representation of the affected forestry area per district with respect to total affected area in Galicia. 
The wildfire risk depends on several climatological, social or environmental factors, which could be modified by public policies at short or medium term. Nevertheless, decisionmaking is characterised by the presence of dynamic risk factors (Rogalski, 1999) [1].
In order to study wildfires in depth, it is necessary to be aware of the currentsituation of the agroforestry areas in Galicia. Certain areas have scattered populations with a constant rural depopulation and continuing migration of young people to more highly populated areas (Marey et al., 2007). In addition, the thinning out of the agroforestry sector within the economy has been clear for some time now, as well as the reduction in employment in this sector. This has contributed to make much more difficult youth employment in the countryside, while at the same time, fewer farms use woodland areas to obtain productive resources (Vega, 2007; Sineiro, 2006). This trend has also caused forestry land to become increasingly neglected, allowing for an increase in the severity and spread of wildfires. Moreover, the structure of forest property is often very divided, with a high level of private property increasingly belonging to elderly people, making it more difficult to manage the lands correctly (Sineiro, 2006).
Galicia contains different climatic areas, resulting in an uneven availability of biomass that can be burnt (Martínez et al., 1999). This makes it difficult to organise the prevention and extinguishing of wildfires. Taking these circumstances into account, together with the social and environmental impacts caused by wildfires, over the last few years the government has prioritized the design of preventative policies, although most of its budget goes toward extinction activities. For these policies to work well, it is important to identify the factors that affect the occurrence of wildfires.
To this end, the proposed model must be simple, structured and easy to standardise so that it can be easily updated (King & MacGregor, 2000).
Until now, several methods have been used to identify wildfire risk factors. Some studies have used various explanatory variables in order to explain the reasons why some areas are more heavily affected than others, although they do not quantify the described relationships and/or support their arguments in a quantitative way (Lavorel et al., 2007). However, other papers use techniques based on Geographic Information Systems (GIS), using probability risk models and linking variables to the forest environment. (Cabrera, 1989; Vilar et al., 2008; Chuvieco et al., 2009; Martínez et al., 2009; RomeroCalcerrada et al., 2010). GIS techniques are used in several models in which geographic and other statistical variables are included (Pew & Larsen, 2001; VegaGarcía & Chuvieco, 2006). Relevant geographic variables include the location of roads, and industrial or recreational areas, amongst other factors (RomeroCalcerrada et al., 2010). Therefore, geographical implications in the occurrence of wildfires have also been widely studied. As a result, this study focuses on assessing the geographical differences in the occurrence of wildfires.
In some earlier work, researchers have studied the error term to identify geographical and temporal trends (Disdier & Head, 2008; Prestemon et al., 2002; Jones, 1991; Moulton, 1986). Testing the error term allows the researcher to control for the unobservable factors across the different entities, implying that this research can be used to determine whether differences across entities are significant. Therefore, an econometric model with random effects (RE) or fixed effects (FE) can be developed in order to account for specific local effects.
Other papers explore the possible relationship between wildfires and a specific group of variables (Finney et al., 2009), including forest management (Prestemon et al., 2002; Butry, 2009; Wimberly et al., 2009), meteorological variables (Aguado et al., 2007), and socioeconomic factors (Mercer & Prestemon, 2005). This research, as well as Prestemon et al. (2002), uses time series models to analyse temporal trends in wildfire risk. Other relevant research also includes socioeconomic variables such as income, machinery used, and/or number of livestock (Vilar et al., 2008). Wildfire risk is also analysed from the perspective of the different phases in the duration of a fire, and the ignition, intensity or area affected are used as dependent variables (Genton et al., 2006).
Several international studies analyze the problem of fires from a spatial context. For example, Prestemon et al. (2002) developed a model with fixed effects as to assess whether there is a spatial behavior in the occurrence of fires between administrative units from North Florida. Preisler et al. (2004) used temporal and spatial effects through a logistic regression to study the probability of fires in Oregon (USA) since 1970. Meanwhile, Brillinger et al. (2006) developed an empirical model for analyzing the evolution of fire risk. Their model contains both, FE and RE to analyze the fires occurrence in California (USA) during the years 20002003. Finally, Chen et al. (2014) also analyze the risks and causes of fires using spatial econometrics.
The aim of this study is to extend previous analyses using current data and taking into account the impact of socioeconomic factors, land cover, and climatology using spatial analysis. Thus, econometric models have been developed to analyse the possible influence of socioeconomic factors on the risk of wildfires. The Ordinary Least Square (OLS) and econometric models for counted data are used to identify these socioeconomic factors. Random effects (RE) and Fixed Effects (FE) are also estimated to assess the presence of spatial patterns. Other methods, such as the Moran’s I and LISA statistics, are included to determine whether wildfire occurrence shows spatial patterns. We expect the present results can help to improve public policy focussing on exploring spatial and temporal impacts on fire occurrence.
This research starts by explaining the data and methods used. In the next section, the results are described and discussed, and then it follows a section in which the main conclusions of the research are summarized and policy implications are provided.
Materials and methodsTop
Data
Data have been gathered from 2001 to 2010. The most uptodate data available from the 19 forest districts established by the Galician regional government were collected (Xunta de Galicia, 2011). Variability over time and between districts will be one of the desirable data properties (Cameron & Trivedi, 2009). Data have been grouped by forest districts in order to have a common geographical reference. Therefore, some variables had to be transformed prior to be included into the model by aggregating municipal data up to the district level. In the following analysis, the hectare is the unit used to measure the surface area.
The explanatory variables are shown in Table 1, and these can be grouped into seven main categories, including: the population structure, weather variables, territorial features, economic information, agroforestry situation, wildfire characteristics and time dummy variables. To avoid perfect multicollinearity in the econometric models, the dummy year for 2001 has been used as a baseline, and time effects are interpreted by using this year as a reference point. As several variables for different groups showed high correlations with each other, a limit of 70% was set for the value of the linear correlation coefficient. Furthermore, the variance inflation factor (VIF) was used to analyse the level of multicollinearity among the chosen variables (Neter et al., 1983). The VIF had values lower than 2.16 for each variable and 3.70 for the set. These values indicate that multicollinearity is not a problem in the selected variables.
Wildfires data were recorded from the Galician Forest Districts. On the other hand, meteorological data were recorded directly by the weather stations, and such data had to be linked and extrapolated to the District level. Finally, the agroforestry data are mainly recorded by Geographic Information Systems (GIS). Thus, a shapefile with the Galician Forest Districts was designed adding the municipality limits obtained from the National Geographic Institute (IGN, 2011). To conclude, agroforestry data were obtained cropping GIS information with the previous defined shapefile.
The data for the climatic variables were collected from MeteoGalicia (2012). The climate stations belonging to each district were geographically located. The average maximum temperature and rainfall recorded per month during the summer were collected[2]. The proportion of the protected areas in each district was also included to describe relevant territorial features. The protected areas were obtained from the MAGRAMA (2010). These data were provided by two maps containing the Community Interest Sites (CIS) and Special Protection Areas for Birds (SPAB). Thus, the GvSig software was used to compute the size of both protected areas by Forest District (GvSig, 2014). In this way, the ratio of protected areas is computed using the total protected area divided by total district area.
The density per hectare is used to describe the population structure. Therefore, the total population divided by the total Forest District area is used to calculate this variable. Both data were recorded from IGE (2012a) and municipality statistics. This density variable presents high correlations with the personal income, level of studies or employment rate. In order to avoid such multicollinearity problems, variables referring to personal income, education and employment rates had to be dropped from the final specification due to their high correlations among each other.
The Third Spanish National Forest Inventory (NFI3), the Corine Land Cover and the Livestock Census were the main sources to gather information about the agroforestry situation (IGE, 2012b). Tree dominant tree species were recorded in order to describe the forest plantations. The forestry areas, in which the Pinus pinaster is the main specie, were calculated from the NFI3 (MAGRAMA, 2008) while the districtforested areas were recorded from IGE (2012b). Hence, the ratio of Pinus pinaster was included in order to describe the forestry structure. The natural pasture was obtained by accounting for the lands where this activity is recorded according to the Corine Land Cover database (European Environment Agency, 2010). The GvSig software was employed to calculate both the agricultural and forestry area (GvSig, 2014). The ratio of equines was also taken into account in order to describe the livestock structure, as in Barreal et al. (2011). This variable was considered given that previous literature related the presence of equines with land management and fuel treatments (Rigueiro et al., 2002). The percentage of equines represents 1% up to 11% of total livestock according to each Forest District data (IGE, 2012b). Although cattle are the main livestock in Galicia; equines usually graze in pasturelands or forest areas. The number of agricultural cooperatives is also included in the model. This variable can be a proxy for the dynamics of the rural areas. These data were collected from IGE (2012b). However, we should note that there are some years in which yearly data are missing, so that the series had to be completed with the closest data points available.
The wildfire variables were also obtained from the IGE (2012b). For the first six years, municipalities provided the data, then the burned area and the number of wildfires had to be aggregated by forest districts. The ratio of burned area was calculated using the total forest area provided by IGE (2012a). The GeoDa software was used to obtain the spatial statistics of the dependent variables (Anselin et al., 2006; GeoDa, 2014). In order to create the final database, and conduct the estimation process, the Stata 10.1 software was used (Stata, 2010).
Methodology
Descriptive Spatial Analysis
Graphs and statistics are useful in order to identify the spatial patterns in Galician wildfires. The first one involves the representation of the data to identify the temporal trends and the heterogeneity between the Galician forest districts. Then, in case of existing temporal trends, these could be identified showing differences of each entity´s mean value. Another alternative is to represent the data for each year by a graph. The independent years could register more or less spatial differences.
The Moran’s I statistic (Moran, 1948) was used for statistical analysis, as well as the Local Indicators of Spatial Association LISA (Anselin, 1995). With both statistics, the spatial dependence can be analysed using the autocorrelation coefficients between the Galician forest districts. The first statistics analyses the spatial heterogeneity of the sample. Meanwhile, the second focuses on the relationship between each geographic unit, identifying the clusters of study.
The Moran’s I statistic takes into account the number of geographical areas (N); the analyzed areas (j and i); the study variable for each location (y); the mean of the variables of interest in all areas (); and finally the weight matrix that describes the relationship between both locations (W_{j,i}). Then, the Moran’s I statistic could be expressed by the Eq. (1).
According to the definition of the weight matrix, the relations between close forest districts are included. Therefore, the closest neighbors to each polygon are identified with this matrix. Mathematically, the weight matrix could be expressed as Eq. (2), in which ŵ_{ij} represents the spatial matrix of adjacent polygon (j), respect to the polygon that we are studying (i).
The spatial relationships in this matrix can be used with different contiguity interpretations. In other words, if a regular grid is designed, the weight matrix could be constructed according to four spatial relationships: linear (Fig.3a), Rook (Fig.3b), Bishop (Fig.3c) and Queen (Fig.3d). These relations depend on the number and directions of spatial dependences that the researchers may find. In this research, the polygons are irregular, so the criterion with more spatial directions is selected (Moreno & Vayá, 2000), and as a consequence, the queen contiguity is chosen to analyze the spatial patterns. This contiguity can be used at several levels (Lesage & Kelley, 2009). This research analyzes the direct relationships between the closed forest districts in terms of fire occurrence.
Figure 3. Types of contiguity for direct relations. 
The LISA statistics can be developed from the Moran’s I statistics (Anselin, 1995). This is described in the Eq. (3) where z_{i} represents the normalized value of the selected variable in respect to the mean and J_{i} is all polygons (districts) next to i. Therefore, the LISA statistics analyzes the spatial patterns between each entity to the selected data. In other words, the spatial autocorrelation is individually analyzed.
Econometric Analysis
In order to analyse the relationship between the previous variables and wildfires in Galicia, a baseline lineal regression estimated by OLS was used. In this baseline estimation the coefficients are controlled by the heterogeneity of each district through the HuberWhite correction of standard errors (Cameron & Trivedi, 2009). Thus, the econometric model is presented by Eq. (4), in which the variables are arranged into a panel according to each district and their respective periods of time. In this equation, the subscripts “j, k, h” represent the type of variable, “i” is the forestry district, and “t” is the period.
With this common specification, two independent equations were estimated. The first model used the ratio of forestburned area in each forest district as the dependent variable, and the second specification modelled the number of wildfires. The independent variables in both models include socioeconomic factors represented by X_{jit}, (mainly population structure, territorial features, economic information and agroforestry data for each forest district), climatology represented by X_{kit} (including the variables of average maximum temperature and average monthly precipitation); and finally, the vector X_{hit} represents the dummy yearly indicators.
Using the BoxCox test, the functional form of the Eq. (4) was selected. The BoxCox test develops a transformed dependent variable represented by the Eq. (5), in which the residual (μ_{it}) assumes a normal distribution in order to estimate the parameters β and θ.
As such, if the estimation of θ is close to zero, then the best specification to be used would be the loglineal model. However, if the respective statistics are significant and close to one, a lineal model should be used. Eq. (6) is then formulated according to the following specification.
Since, the number of wildfires is a counted data variable, the Poisson Regression Model (PRM) shown on Eq. (7) is employed, with the specification earlier presented in Eq. (4):
Given that count data can exhibit overdispersion (Cameron & Trivedi, 2005), we need to assess whether this is present by estimating Eq. (8). Overdispersion implies that the variance depends on the mean plus square parameter (α^{2}). In this case if α=0, then the variance is equal to the mean and there is no overdispersion; and thus, the PRM can be a suitable model.
On the other hand, if the coefficient α is different from zero, then the number of wildfires should be estimated by a Negative Binomial Regression model (NBRM). This model is more general than the PRM and should prove to have a better goodness of fit in case of overdispersion (Cameron & Trivedi, 2009).
In order to interpret the coefficients of the previous model, the use of the Incidence Rate Ratio (IRR) is recommended as its results are easier to interpret (Long & Freese, 2001). As such, the IRR coefficients are estimated to directly quantify the values of the respective parameter estimates. This ratio is calculated by Eq. (9), in which the results can be analysed as a change in the probability of wildfire occurrence, when there is a change in the analysed independent variable, whenever the others parameters are constant.
In this setting, two different models could be used to analyse the error term: FE and RE. The FE represented by the Eq. (10) in which the error term (ε_{it}) of the Eq. (4) is broken into two parts: one fixed term (v_{i}) and another error term (τ_{it}).
In the following RE models, the previously fixed term (v_{i}) is now random. The specification of this model is equal to Eq. (10), but in this case the random term will have a mean of v_{i} and different variance from zero (Var(v_{it})≠0). These unobservable factors are used for the OLS model but also for the MRP (Eq. 11), or NBMR in case of overdispersion.
The Hausman test (H) is used to select between RE and FE models. The specification of this test is shown in Eq. (12) and analyzes the consistency of estimators for both models. The null hypothesis states that there is no correlation between the unique errors and the independent variables. This hypothesis is tested at the 5% significance level. If the null hypothesis is not rejected, then FE are selected over RE. Otherwise, RE should be used.
ResultsTop
Spatial patterns and temporal trends analysis
The spatial patterns of the number of wildfires and burnedforest area ratio can be observed with graphical displays. The variation of the burnedforest area ratio by year is represented in the Figure 4a. In this graph the xaxis represents the Galician forest districts according to each number (Xunta de Galicia, 2011). Different values between the districts are recorded in all graphs; however its difference depends on the year. Another way to identify the existence of spatial patterns is by using the average of the burnedforest area ratio for the sample. This is included in the Figure 4b where the difference between districts can be observed. Also, the temporal trends are observed per year.
Figure 4. Graphical representation for variation of burnedforest area ratio in Galicia from 2001 to 2010. a) Data represented by year. b) Data recorded by year and the mean of each district. 
Figure 5a and Figure 5b describe the evolution of wildfires per year from 2001 to 2010. The spatial patterns can be identified in this graph. Figure 5b shows also spatial patterns in the mean of wildfires according to each district. Data show significant differences across years, therefore the number of wildfires contains also temporal effects.
Figure 5. a) The number of wildfires represented by year in Galicia from 2001 to 2010. b) The number of wildfires recorded each year according to Galician Forest Districts from 2001 to 2010. 
A weight matrix should be constructed to develop the spatial statistics. As stated earlier, the direct Queen contiguity is selected to analyse the relationship between districts (Fig. 3d) and its histogram is represented in Figure 6. In this graph, it can be observed the lowest and highest contiguity between forest districts. Thus, Galician forest districts have at the minimum two influential neighbours and six as a maximum. This histogram also highlights a big number of districts with six entities around them. However, the biggest contiguity group has only three neighbours.
Figure 6. Histogram of contiguity according to Level 1 for Galician Forest District. 
The Moran’s I and LISA statistics are used to analyse the spatial patterns of wildfires and the yearly burnedforest area ratio in Galicia. Figure 7 reports the Moran’s I statistic by year, and its average for the entire sample. It can be observed that Moran’s I registers higher values for the number of wildfires than for the ratio of burnedforest area. Thus, more spatial autocorrelation is detected for the number of wildfires than for the ratio of burnedforest area. Also, in 2007 no spatial correlation is found for burnedforest area. This may be explained because in 2006 wildfires affected many areas (Molano et al., 2007). This caused social alarm, therefore over the next year, the wildfires occurrence has been drastically reduced. Even the recorded data for the burnedforest area diminished a 97% in 2007 with respect to 2006.
Figure 7. The Moran’s I statistic for the burnedforest area ratio and the number of wildfires recorded in Galicia during 20012010. 
The LISA statistic represents the various significant spatial patterns as follows (Anselin, 1995; Lesage & Kelley, 2009):
 HighHigh (HH): a particular forest district and their neighborhoods have high values. This type of relationship is represented by the red color.
 HighLow (HL): a particular forest district has high values and their neighborhoods have lower values. This type of relationship is represented by the pink color.
 LowHigh (LH): is similar to the previous category, but in this case the forest district has high values and their neighborhoods have lower values. This type of relationship is represented by the skyblue color.
 LowLow (LL): the forest district and their neighborhoods have low values. This type of relationship is represented by the blue color.
The remaining values are represented by a grey color because these entities have a random relationship (Moreno & Vayá, 2000). Figure 8 represents the LISA statistics related to each dependent variable. The colored results are significant at the 5% level. Taking into account the burnedforest area ratio, the LISA statistics is represented in Figure 8. In each map the LISA statistic for each district is represented according to the relation with its neighborhoods. With this result, the LowLow (LL) relation is mainly recorded in the North of Galicia, although, this relationship was also recorded in the South for some particular years. However, the HighHigh (HH) relationship occurs primarily in the South. Thus, Southern forest districts and their neighbourhoods record high values of burnedforest area ratio.
Figure 8. LISA statistic for the burnedforest area ratio in Galicia during 20012010. 
In order to analyze the number of wildfires, the LISA statistics are shown in Figure 9. In these maps the LowLow (LL) relationship could be observed in the Northeast districts. On the other hand, HigherHigher (HH) relations are recorded in the South and Southeast. Therefore, there is evidence of spatial patterns in the number of wildfires.
Figure 9. LISA statistic for the number of wildfires in Galicia during 20012010. 
All previous graphs and statistics show the existence of relevant spatial patterns and temporal trends. Therefore, these should be included in the econometric model for both dependent variables. The temporal trends are included in the empirical models by using dummy variables for each year, considering 2001 as the baseline year. On the other hand, in order to correct for spatial patterns in the research, data are set according to a panel of forest districts and controlling the heterogeneity by district through standard errors correction. The spatial patterns are also analyzed using FE and RE models.
Econometric models results
Results for the burnedforest area
In order to specify the most suitable econometric model to analyse the evolution of the burnedforest area ratio, a BoxCox test was estimated (Cameron & Trivedi, 2009), being its results reported in Table 2. A logarithmic model is used in accordance with the results obtained in the BoxCox test. In other words, the statistics could not reject the logarithmic specification both for the dependent and independent variables.
Following the results displayed in Table 3, the estimation by OLS captures 69.04% of the variation of the burned forest area rate. Taking into account the Fstatistic, we find that all parameters are jointly significant. As regards the choice between the use of FE and RE in the previous models, the Hausman test recommends the use of RE to estimate the burnedforest ratio model (Prob>chi2= 0.98).
The dummy variables determine significant effects over several years. A positive trend is identified from 2002 to 2006. The majority of dummy variables are significant and positive with respect to the 2001 year. However, after 2006, the trend is clearly negative and significant for all years. These results are robust across the econometric selected models.
The climatological variables show in particular the importance of rainfall in order to reduce the burnedforest area ratio. This variable is significant carrying a negative effect in the causes of wildfire occurrence. Thus, the average effect of rainfall on burnedforest areas ratio is 0.643, when the precipitation changes by one unit over time and between districts. However, the maximum temperature has a positive effect, although this variable is not significant in order to predict the burned area. The small variability in this variable may be responsible for this finding.
In terms of socioeconomic variables, the ratio of the equines and the number of agricultural cooperatives have both a negative and significant effect on the burnedforest ratio area for the OLS and RE models. Their effects show that if the value changes over time and between districts by one unit, then the average effect of equine radio stock and the number of agricultural cooperatives over the burnedforest area will respectively decrease by a factor of 0.385 and 0.555.
On the other hand, the density of Pinus pinaster and the ratio of protected areas are positively related with this dependent variable. The coefficients are significant for the OLS and RE results. We also find that the ratio of natural pasture has no statistical impact on any of the econometric models.
Results of the number of wildfires
For the purpose of determining the functional form for the regression of wildfires, the BoxCox test does not provide conclusive evidence of the superiority of any functional form. However, in order to compare the number of wildfires with the regression of the burnedforest area ratio, the logarithmic model is selected. Also, this functional form is estimated to allow for comparability between all regressions.
Following the results displayed in Table 5, the estimation by OLS explains 78.60% of the variation of the number of wildfires. In addition, the parameters are all jointly statistically significant. In the OLS results, temporal trends are also identified. Until the year 2005, the coefficients are not significant; however, from this year onwards, all yearly dummies are significant and negative. Therefore, from 2005 onwards, the wildfire occurrence diminishes with respect to 2001. Taking into account the climatological variables, the rainfall carries a significant and negative effect on the number of wildfires (0.430). On the opposite, the maximum temperature is significant and positively related to wildfire occurrence (1.696).
Some variables, such as the ratio of equines and the agricultural cooperatives do not have a significant relationship with the number of wildfires during 20012010, are not significant in the assessment of the wildfires using the OLS models. However, the rest of the socioeconomic variables, are significant and have positive effects over the wildfires occurrence.
The number of wildfires is modelled by count data models. Therefore, overdispersion should be studied in order to select the best econometric model. Taking into account the results of Table 3, the data show overdispersion, and hence, the NBMR is selected to estimate the number of wildfires. The Hausman test recommends the use of RE to estimate the number of fires by the NBMR (Prob.>chi2=1.00).
Analysing the effects of the yearly variables, temporal trends are found according to the NBMR results, both with FE or RE. In this way, since 2006, it is observable that the wildfire occurrence diminishes with respect 2001. However, the NBRM also detects a significant growth in wildfires in 2005 with respect to the baseline year. The OLS and NBMR models, with or without RE, demonstrate the presence of temporal trends in Galician wildfires.
Furthermore, the estimator of summer rainfall is significant and carries a negative effect on the number of wildfires (0.994). According to the estimation with RE, if this independent variable changes over time and between districts by one unit, then the average effect of the average summer rainfall over the number of wildfires is significant (0.944). Otherwise, the average of the maximum temperature during the summer is not significant to explain the wildfires according with the NBMR models.
In the NBMR models, the ratio of natural pasture, Pinus pinaster and protected areas are statically significant. The effects of these variables on wildfire occurrence are positive. By analysing the IRR, if the ratio of natural pasture, the ratio of Pinus pinaster and the ratio of protected areas per landowner show an increase by one unit, then the number of wildfires increases by a factor of 149.018, 5.437 and 3.276, respectively.
In addition, socioeconomic variables are significant in the NBRM with RE. Nevertheless, the remaining variables have different impacts on wildfire occurrence. The agricultural cooperatives and population density have a positive relationship with the occurrence of Galician wildfires. If these previous variables increase by one point, the rate of the number of wildfires would be expected to increase by a factor of 1.015 and 1.180, respectively, while holding all other variables constant. Furthermore, the ratio of equines has a negative relationship with wildfires occurrence.
The summer average rainfall is significant in order to predict the wildfires occurrence. This is explained by the absence of raining, given that this increases the wildfire risk. Nevertheless, the summer maximum temperature is only significant in the OLS results.
DiscussionTop
Spatial patterns and temporal trends can be observed with graphical data representation. Furthermore, the spatial dependence of wildfires can also be determined by spatial statistics. Various econometric models are employed to assess the impact of socioeconomic, climatic and geographical variables, as well temporal and spatial effects. Following the econometric models employed, and in particular those from RE models, the number of wildfires and the affected area ratio are estimated for each Galician forest districts in 2010. The estimations portrayed in Figure 10 show the actual data for both dependent variables and predictions. The data of these variables are distributed in quantiles and represent each district. In doing so, the geographical patterns of wildfire occurrence can be clearly differentiated in these maps. It is shown that the wildfire risk depends on the forest district; and as such, regulators should focus their forestry efforts on the areas in which the prediction of wildfires is higher. In other words, and for the area of study, public policy efforts should focus more closely on the southern rather than the northern districts.
Figure 10. Estimation of wildfire occurrence in 2010. a) Actual ratio of burned area (%). b) Estimated rate of burned area (%). c) Actual number of wildfires. d) Estimated number of wildfires. 
In terms of the econometric results, it was found that the agroforestry features are important factors given that the land cover is conditioned by this activity. The type of forest plantation, the livestock used in the farms or the land assigned to agricultural activity influences the wildfire occurrence. The ratio of equines is slo important in order to reduce the wildfire occurrence (Rigueiro et al., 2002; Pasalodos et al., 2009). This species grazes freely in the surrounding farm; fed mainly with grass, bushes or seeds; keeping the land cover cleaner. Thus, the wildfire risk diminishes where there are more equines than other livestock species. The presence of agricultural cooperatives also affects the wildfires occurrence. This happens because of the traditional agricultural management practices using fires. Nevertheless, the effect is the opposite for the burned area, because in general terms, the lands are better managed when the agricultural sector is more powerful in rural areas.
We also find that the Pinus pinaster ratio has a positive influence on the occurrence of wildfires, because this species is pyrophyte, with wildfires spreading more where Pinus pinaster are being planted. Protected areas could also be expected to have a negative relationship with wildfire occurrence; however the social rejection or ineffective protection measures could cause a positive influence. Furthermore, climatology variables condition the occurrence of wildfires, how they spread or the suppression efforts. Finally, the evolution of wildfires over time demonstrates high variability. This is the justifying reason why yearly dummy variables were included.
In general terms, the population density is important in order to predict the wildfire occurrence. However, in some of the empirical models, results are not conclusive. For example, in the OLS with RE, this variable is not significant when explaining the burnedforest area ratio. The same happens in the NBRM when analyzingthe wildfire number. In the remaining models, the population density is positively related to the occurrence of wildfires. This result is explained by the progressive migratory flow from the rural to urban areas. Then, the wild landurban density around to main areas is increasing. In addition, the new residents are not involved in the agricultural sector and they are not involved with forest production (Barreiro & Hermosilla, 2013). This generates worst environmental conditions, causing an increase of wildfires (HerreroCorral et al., 2012). Public policies may supervise the surrounding environment of these areas and aware society to avoid wildfire occurrence.
The types of forest covers are represented by the ratio of Pinus pinaster. The results show a positive influence on both, the ratio of burnedforest area and wildfires number. These results are related with the species characteristics because these are more inflammable and the wildfires, when occurring, move faster than with other species. The preventive measures and supervision should also be incremented in these areas in order to avoid wildfires.
Unexpectedly, the protected areas influence positively the occurrence of wildfires. This result may show the general rejection towards having protected lands in rural areas. This result could also imply an inadequate public policy to manage these areas against wildfires (Carroll et al., 2006). Therefore, the zooning of protection areas may be revised in order to identify the possible social and environmental factors that can be improved when reducing management conflicts. These improvements will imply lower wildfire occurrence if these factors are corrected.
ConclusionTop
This research provides evidence characterizing the wildfire occurrence in the agricultural sector in relation to the climatic conditions, the forest cover, the social context, and time and spatial patterns. A relevant finding is that the forest species and the farming systems condition the wildfire risk. Hence, public policies may mitigate the factors that affect the wildfire risk. In this way, the presence of equines and extensive agricultural practices should be promoted in order to reduce the wildfire risk.
According to the main results, some guidelines could be developed as a reference for regional and local governments to help in the fight of wildfires. In particular, public policies could regulate the quality and quantity of woodland made available, as well as the plantation of different species. These regulatory agencies should also consider the geographical and spatial differences in the occurrence of wildfires in order to formulate better forest policies, and deal with possible “contagion” effects across districts.
Finally, we should remark that the current research has some limitations. In particular, additional variables would be desirable by employing more geographical disaggregated data, such as roads and infrastructures. Unfortunately, such data are not currently available, although, they are expected to be in the near future. In spite of that, many of the obtained conclusions could be applicable to other similar European areas, especially in depopulated rural areas. In particular, these main results could be implemented in the French Mediterranean basin (INSEE, 2015; PROMETHEE, 2015) and Portugal (INE, 2015), among others.
NotesTop
Following Molano et al. (2007) and Martínez et al. (2009), the causing factors of wildfires can be divided into two main categories: avoidable and unavoidable. Unavoidable causes are considered those that cannot be foreseen or dissuaded, whereas avoidable causes are those that can be prevented through individual actions or forestry policies. This implies that there are exogenous factors, which are uncontrollable, to which other endogenous factors must be added. In general, the unavoidable category contains natural phenomena, whilst the avoidable can be divided into three possible subcategories: intentional, negligent and unknown causes. Avoidable causes represent almost all of the causes, although the majority of these are classified as unknown, showing that the causality of wildfires is not recorded reliably and depends heavily on the criteria of investigators (Pérez & Delgado, 1995; Molano et al., 2007). 

In some cases, climatological data were not available for all of the time periods and/or forest districts. The unavailable data had to be supplemented with those from other forest districts according to the climatic areas established by Martínez et al. (1999). 
ReferencesTop
○  Aguado I, Chuvieco E, Borén R, Nieto H, 2007. Estimation of dead fuel moisture content from meteorological data in Mediterranean areas. Applications in fire danger assessment. International Journal of Wildland Fire 16(4): 390397. http://dx.doi.org/10.1071/WF06136 
○  Anselin L, 1995. Local Indicators of Spatial AssociationLISA. Geographical Analysis 27(2): 93115. http://dx.doi.org/10.1111/j.15384632.1995.tb00338.x 
○  Anselin L, Syabri I, Kho Y, 2006. GeoDa: An Introduction to Spatial Data Analysis. Geographical Analysis 38 (1): 522. http://dx.doi.org/10.1111/j.00167363.2005.00671.x 
○  Barreal J, Loureiro M, Picos J, 2011. Estudio de la incidencia de los incendios en Galicia: una perspectiva socioeconómica. Revista Galega de Economía 20 (special issue): 227246. 
○  Barreiro JB, Hermosilla T, 2013. Sociogeographic analysis of the causes of the 2006’s wildfires in Galicia (Spain).Forest Systems 22(3): 497509. http://dx.doi.org/10.5424/fs/201322304165 
○  Brillinger DR, Preisler HK, Benoit JW, 2006. Probabilistic risk assessment for wildfires. Environmetric 17: 623–633. http://dx.doi.org/10.1002/env.768 
○  Butry DT, 2009. Fighting fire with fire: estimating the efficacy of wildfire mitigation programs using propensity scores. Environmental and Ecological Statistics 16(2): 291319. http://dx.doi.org/10.1007/s1065100700833 
○  Cabrera AD, 1989. Metodología de estudio de los incendios forestales: el caso de la provincia de Córdoba. ERIA: 139146. 
○  Cameron AC, Trivedi PK, 2005. Microeconometrics. Methods and applications. Cambridge University Press, Cambridge, England. 1034 pp. http://dx.doi.org/10.1017/CBO9780511811241 
○  Cameron AC, Trivedi PK, 2009. Microeconometrics using Stata. Stata Press, Texas, USA. 692 pp. 
○  Carroll MS, Higgins LL, Cohn PJ, Burchfield J, 2006. Community Wildfire Events as a Source of Social Conflict. Rural Sociology 71(2): 261280. http://dx.doi.org/10.1526/003601106777789701 
○  Chen X, Goodwin BK, Prestemon JP, 2014. Is Timber Insurable? A Study of Wildfire Risks in the US Forest Sector Using SpatioTemporal Models.American Journal of Agricultural Economics 96(1): 213231. http://dx.doi.org/10.1093/ajae/aat087 
○  Chuvieco E, Aguado I, Yebra M, Nieto H, Salas J, Martína MP, Vilar L, Martínez J, Martín S, Ibarra P, De La Riva J, Baeza J, Rodríguez F, Molina JR, Herrera MA, Zamora R, 2009. Development of a framework for fire risk assessment using remote sensing and geographic information system technologies. Ecological Modelling 221(1): 4658. http://dx.doi.org/10.1016/j.ecolmodel.2008.11.017 
○  Disdier AC, Head K, 2008. The puzzling persistence of the distance effect on bilateral trade. The Review of Economics and Statistics 90(1): 37–48. http://dx.doi.org/10.1162/rest.90.1.37 
○  European Environment Agency, 2010. CORINE Land Cover. Available at: http://www.eea.europa.eu/. 
○  Finney M, Grenfell IC, McHugh CW, 2009. Modeling Containment of Large Wildfires Using Generalized Linear MixedModel Analysis. Forest Science 55(3): 249255. 
○  Forest Europe, UNECE, FAO, 2011. State of Europe’s Forests 2011. Status and Trends in Sustainable Forest. Ministerial Conference on the Protection of Forests in Europe, Aas, Norway. 342 pp. 
○  Genton MG, Butry DT, Gumpertz ML, Prestemon, JP, 2006. Spatiotemporal analysis of wildfire ignitions in the St Johns River Water Management District, Florida. International Journal of Wildland Fire 15: 8797. http://dx.doi.org/10.1071/WF04034 
○  GeoDa, 2014. 1.6.6. GeoDa Center for Geospatial Analysis and Computation. Available at: https://geodacenter.asu.edu/ 
○  GvSig, 2014. 2.0. GvSig Desktop. Asociación GvSig. Available at: http://www.gvsig.org 
○  HerreroCorral G, Jappiot M, Bouillon C, LongFournel M, 2012. Application of a geographical assessment method for the characterization of wildland–urban interfaces in the context of wildfire prevention: A case study in western Madrid. Applied Geography 35(1): 6070. http://dx.doi.org/10.1016/j.apgeog.2012.05.005 
○  IGE, 2012a. Información Territorial de Galicia (Spain) [online]. Instituto Galego de Estadística, Available at: http://www.ige.eu/web/mostrar_actividade_estatistica.jsp?idioma=gl&codigo=0101001 [17 July 2012]. 
○  IGE, 2012b. Banco de datos municipal [online]. Instituto Galego de Estadística, available at: http://www.ige.eu/web/mostrar_paxina.jsp?paxina=002001&idioma=gl [17 July 2012]. 
○  IGN, 2011. Base Cartográfica Nacional. Instituto Geográfico Nacional. Available at: http://centrodedescargas.cnig.es/CentroDescargas/catalogo.do#selectedSerie 
○  IGN, 2011. Líneas Límite de la Base de datos de Límites Jurisdiccionales de España inscritas en el Área del Registro Central de Cartografía del IGN. Instituto Geográfico Nacional. Available at: http://centrodedescargas.cnig.es/CentroDescargas [17 July 2012]. 
○  INE, 2015. Estatisticas territoriais (online). Instituto Nacional de Estatísitca (Spain), Available at: http://www.ine.pt/xportal/xmain?xpid=INE&xpgid=ine_unid_territorial&menuBOUI=13707095&contexto=ut&selTab=tab3 [02 March 2012]. 
○  INSEE, 2015. Données locales (online). Institut National de la Statistique et des Études Économiques (France). Available at: http://www.insee.fr/fr/basesdedonnees/default.asp?page=statistiqueslocales.htm [02 March 2012]. 
○  Jones K, 1991. Specifying and Estimating MultiLevel Models for Geographical Research. Transactions of the Institute of British Geographers 16(2): 148159. http://dx.doi.org/10.2307/622610 
○  King D, MacGregor C, 2000. Using Social indicators to measure community vulnerability to natural hazards. The Australian Journal of Emergency Management 15 (3): 5257. 
○  Lavorel S, Flannigan M, Lambin E, Scholes M, 2007. Vulnerability of land system to fire: Interactions among humans, climate, the atmosphere, and ecosystems. Mitigation and Adaptation Strategies for Global Change 12(1): 3353. http://dx.doi.org/10.1007/s1102700690465 
○  Lesage J, Kelley R, 2009. Introduction to Spatial Econometrics. CRC Press. Florida (EUA). http://dx.doi.org/10.1201/9781420064254 
○  Long JS, Freese J, 2001. Regression models for categorical dependent variables using Stata. Stata Press Publication, Texas, USA. 288 pp. 
○  Marey MF, Rodríguez V, Crecente R, 2007. Perfil del propietario forestal individual en Galicia: objetivos y prácticas de gestión en el noroeste de la comunidad. Revista Galega de Economía 16 (1): 124. 
○  MARM, 2012. Avance Anuario de estadística 2011 [online]. Ministerio de Medio Ambiente y Medio Rural y Marino de España, available at: http://www.magrama.gob.es/estadistica/pags/anuario/2011/AE_2011_Avance.pdf [17 July 2012]. 
○  Martínez A, Castillo F, Pérez A, Valcárcel M, Blanco R, 1999. Atlas climático de Galicia. Xunta de Galicia, Santiago de Compostela, Spain. 277 pp. 
○  Martínez J, VegaGarcía C, Chuvieco E, 2009. Humancaused wildfire risk rating for prevention planning in Spain. Journal of Environmental Management 90(2): 12411252. http://dx.doi.org/10.1016/j.jenvman.2008.07.005 
○  Mercer DE, Prestemon JE, 2005. Comparing production function models for wildfire risk analysis in the wildlandurban interface. Forest Policy and Economics 7: 782795. http://dx.doi.org/10.1016/j.forpol.2005.03.003 
○  Meteogalicia, 2012. Informes Meteorolóxicos [online]. Meteogalicia, available in: http://www.meteogalicia.es/observacion/informesclima/informesIndex.action [17 July 2012]. 
○  MAGRAMA, 2008. Tercer Inventario Forestal Nacional 19972007. Galicia. Ministerio de Medio Ambiente. Madrid, España. 
○  MAGRAMA, 2010. Red Natura 2000. Ministerio de Agricultura, Alimentación y Medio Ambiente. Dirección General de Calidad y Evaluación Ambiental y Medio Natural. Subdirección General de Medio Natural. Available at: http://www.magrama.gob.es/es/biodiversidad/servicios/bancodatosnaturaleza/informaciondisponible/rednatura2000_descargas.aspx [17 July 2012]. 
○  Molano F, Rodríguez C, Ponte JM, 2007. Informe sobre investigación de incendios en Galicia. Verano 2006. Diputación de Coruña, A Coruña, Spain. 270 pp. 
○  Moran P, 1948. The interpretation of statistical maps. Biometrika 35: 255260. 
○  Moreno R, Vayá E, 2000. Técnicas econométricas para el tratamiento de datos espaciales: La econometría espacial. Edicións de la Universitat de Barcelona. Barcelona (Spain). 156 pp. 
○  Moulton BR, 1986. Random group effects and the precision of regression estimates. Journal of Econometrics 32: 385397. http://dx.doi.org/10.1016/03044076(86)900217 
○  Neter J, Nachtsheim C, Kutner M, 1983. Applied linear regression models. Richard D. Irwin, Chicago, USA. 547 pp. 
○  PasalodosTato M, Pukkala T, RigueiroRodríguez A, FernándezNúñez E, MosqueraLosada MR, 2009. Optimal management of Pinus radiata silvopastoral systems established on abandoned agricultural land in Galicia (northwestern Spain). Silva Fenn 43(5): 831845. http://dx.doi.org/10.14214/sf.176 
○  Pérez J, Delgado JL, 1995. Análisis del riesgo de incendio forestal en Galicia (Spain). Revista Agricultura y Sociedad 77: 109124. 
○  Pew KL, Larsen CPS, 2001. GIS analysis of spatial and temporal patterns of humancaused wildfires in the temperate rain forest of Vancouver Island, Canada. Forest Ecology and Management 140(1): 1–18. http://dx.doi.org/10.1016/S03781127(00)002711 
○  Prestemon JP, Pye JM, Butry D.T, Holmes TP, Mercer DE, 2002. Understanding broadscale wildfire risks in a humandominated landscape. Forest Science 48(4): 685693. 
○  Preisler HK, Brillinger D.R, Burgan RE, Benoit JW, 2004. Probability based models for estimation of wildfire risk. International Journal of Wildland Fire 13(2): 133142. http://dx.doi.org/10.1071/WF02061 
○  PROMETHEE, 2015. Forest fires database for Mediterranean area in France. Promethee. Available at: http://www.promethee.com/ [02 March 2012]. 
○  Rigueiro A, Mosquera MR, Villari JJ, 2002. Reducción del riesgo de incendios forestales mediante el pastoreo del caballo gallego de monte. Cuadernos de la Sociedad Española de Ciencias Forestales 14. 
○  Rogalski J, 1999. Decisionmaking and management of dynamic risk. Cognition, Technology & Work 1(4): 247256. http://dx.doi.org/10.1007/s101110050021 
○  RomeroCalcerrada R, BarrioParra F, Millington JDA, Novillo CJ, 2010. Spatial modelling of socioeconomic data to understand patterns of humancaused wildfire ignition risk in the SW of Madrid (central Spain). Ecological Modelling 221(1): 3445. http://dx.doi.org/10.1016/j.ecolmodel.2009.08.008 
○  SECF (Sociedad Española de Ciencias Forestales), 2010. Situación de los bosques y del sector forestal en España [online], available in: http://www.secforestales.org/web/images/inforestal2010.pdf [17 July 2012]. 
○  Sineiro F, 2006. As causas estructurais dos incendios forestais en Galica. In: Os incendios forestais en Galicia (DíazFierros F, Baamonde P, eds). Consello da Cultura Galega, Santiago de Compostela, Spain. 7792 pp. 
○  Stata, 2010. 10. 1. Data analysis and statistical software. StataCorp LP. Available at: www.stata.com 
○  VegaGarcía C, Chuvieco E, 2006. Applying local measures of spatial heterogeneity to LandsatTM images for predicting wildﬁre occurrence in Mediterranean landscapes. Landscape Ecology 21: 595–605. http://dx.doi.org/10.1007/s1098000541195 
○  Vega JA, 2007. Impacto de los incendios sobre suelo y vegetación forestales en Galicia y desarrollo de una selvicultura preventiva. In: Por unha nova cultura forestal fronte aos incendios (DíazFierros F, Balboa X, Barreiro, XL, eds.). Consello da Cultura Galega y Fundación Caixa Galicia, Santiago de Compostela, Spain. 87126 pp. 
○  Vilar L, Martín MP, Martínez J, 2008. Empleo de técnicas de regresión logística para la obtención de modelos de riesgo humano de incendio forestal a escala regional. Boletín de la A.G.E. (Spain). 47: 529. 
○  Wimberly MC, Cochrane MA, Baer AD, Pabst, K, 2009. Assessing fuel treatment effectiveness using satellite imagery and spatial statistics. Ecological Applications 19: 1377–1384. http://dx.doi.org/10.1890/081685.1 
○  Xunta de Galicia, 2011. Pladiga 2011. Consellería de Medio Rural. Dirección Xeral de Montes, Santiago de Compostela, Spain. 123 pp. 
Webpage: www.inia.es/Forestsystems