Assessing the Influence of Roads on Fire Ignition : Does Land Cover Matter ?

In human-affected fire environments, assessing the influence of human activities on the spatial distribution of wildfire ignitions is of paramount importance for fire management planning. Previous studies have shown that roads have significant effects on fire ignition. However, since different land cover classes are subject to different levels of ignition risk, roads in different land cover classes may differently affect fire ignition. The aim of this paper is thus to assess the influence of roads on fire ignition in selected land cover classes subjected to different levels of anthropogenic pressure in Sardinia (Italy). Our results show that fires are preferentially ignited close to roads in all land cover classes. However, the influence of roads is much stronger in less impacted land uses, where the availability of human-induced ignitions highly depends on the accessibility networks. Our approach represents a first step towards the systematic integration of interacting fire ignition drivers such as roads and land cover into fire risk analysis.


Introduction
Fire characteristics such as type, intensity, spread, and temporal distribution depend on several factors such as climate, fuel load and continuity and ignition sources [1].In this respect, humans play a manifold role, acting indirectly as fuel shapers through land use and related vegetation, as well as directly by setting and controlling fire [2].As a consequence, human influence on fire regimes highly varies from region to region worldwide [3].
Several studies demonstrated that in the Mediterranean areas intentional or accidental human-induced fire ignitions account for more than 90% of forest fires [4], thus far exceeding natural fires [5,6].In such human-affected fire environments, a detailed understanding of the influence of anthropogenic activities on the temporal and spatial distribution of wildfires is a crucial factor for fire management planning [7][8][9].Given the complexity of the human component of wildfires, quantitative analyses of fire risk have generally paid less attention to the socioeconomic aspects of fire occurrence compared to the biological and physical factors, such as elevation, slope, temperature, rainfall, or vegetation [10][11][12][13][14][15].In recent years, however, fire risk studies have started to incorporate human-related variables that directly or indirectly influence spatial patterns of fire occurrence.At the landscape scale, a number of authors showed that fire ignitions exhibit strong preferences for given land use/land cover (LULC) types [16][17][18][19][20]. LULCs are the result of human activities, such as agriculture, farming, silviculture, recreation, etc. and are subjected to a different fire ignition risk depending on fuel load and human-induced ignitions.The intrinsic fire selectivity of each LULC is further modulated by the presence of infrastructure such as housing or roads, which are also well known to influence fire ignition [21][22][23][24][25][26][27][28][29].However, while LULC data and distance from roads have been already used for modeling fire ignition risk e.g., [30][31][32][33], to the best of our knowledge, a detailed analysis of the combined effect of road network and land cover is provided here for the first time.
The objective of this paper is thus to assess the influence of roads on fire ignition in selected land cover classes subjected to different levels of anthropogenic pressure in Sardinia (Italy).To achieve this aim, we referred to a selectivity index according to which certain resources (i.e., land cover classes) are more fire-prone than others, being differentially selected by fire, depending on fuel load, structure, spatial continuity and human activities [5,[16][17][18].Fire is considered selective towards a resource (i.e., land cover type) when such resource is used disproportionately to their availability [16,17].By integrating the selectivity behaviour of fire towards LULC and the road network, this paper can provide insights into the interaction between these two fire drivers.

Study Area
Sardinia is the second largest island of the Mediterranean basin, with a total surface area of approximately 24,000 km 2 (Figure 1).The island is characterized by a prevalently hilly topography and extreme heterogeneity in geological and geomorphological features.The highest peak is Punta La Marmora (1834 m a.s.l.), located in central Sardinia.The climate of Sardinia is typically Mediterranean with a mild and rainy period from October to May and a warm and dry period from June to September when temperatures may reach peaks above 30 • C. The average annual amount of precipitation ranges from 500 mm along the coasts to more than 1200 mm in the higher regions [34].
Sardinia hosts a very high ecosystem diversity, with habitats ranging from coastal sand dunes to extensive maquis and forests in the mountains.Land use is shaped by a long human history.Urban areas cover 3% of Sardinia, and the two main cities of Cagliari and Sassari are where most of the 1.7 million inhabitants are concentrated.Agricultural areas represent 25% of the island, whereas approximately 15% is covered by broadleaf forests, mainly composed of Mediterranean oaks, such as Quercus ilex, Q. suber and Q. pubescens.Mediterranean maquis and garrigue cover an additional 23%.
Sardinia is also considered to be a wildfire hotspot, with more than 2500 events per year on average [35].Like in many other regions of southern Europe, the majority of fires directly originate from human activities.However, in recent decades the pronounced summer aridity combined with the abandonment of the traditional agriculture and the resulting progressive built up of fuel increased the risk of large and severe fires.Mean yearly burnt area is approximately 14,500 ha.More than 86% of fires occur from June to September with a peak in July.Fire size ranges over many orders of magnitude from 0.01 ha to >1000 ha, although less than 1% of total events get larger than 100 ha.However, fires > 100 ha account for roughly 51% of the total area burned in Sardinia between 2000 and 2015.Detailed statistics on fire causes are not available.Like in most Mediterranean regions, we can however assume that almost all fire events are of anthropogenic origin.

Data Preprocessing
A dataset of 45,976 wildfire events that occurred in Sardinia between 2000 and 2015 was extracted from the digital database of the Regional Forest Service (Figure 1).For each record, the dataset includes the date and the geographic coordinates of fire ignition and a field estimate of the burned area.The ignition points of the most recent fires (from 2006 onwards) are usually accurately geolocated with GPS, while the positional accuracy of the older ignition points is more uncertain.
For analyzing the relationship between fires, land cover and road network, we used a road network map containing all state, regional and county roads of Sardinia, together with a LULC map of the island at scale 1:25,000.The LULC map was built by the Regional Administration of Sardinia according to the standard classification methodology of the CORINE Land Cover legend [36] based on 2003-2004 digital orthophotos and high-resolution Ikonos satellite images (2005-2006) combined with field surveys.Both maps were downloaded from the Regional geographic information portal of Sardinia (www.sardegnageoportale.it).
We next selected eight land cover classes with a sufficient number of fire records for studying fire incidence patterns at the landscape scale: (i) discontinuous urban areas, (ii) arable land, (iii) mixed agriculture, (iv) olive groves, (v) deciduous forests, (vi) maquis, (vii) garrigue, (viii) pastures.The selected classes were grouped into two main categories with different levels of anthropogenic pressure according to the CORINE Land Cover classification [36]: Artificial and agricultural classes (i-iv) and Natural and semi-natural classes (v-viii).

Data Preprocessing
A dataset of 45,976 wildfire events that occurred in Sardinia between 2000 and 2015 was extracted from the digital database of the Regional Forest Service (Figure 1).For each record, the dataset includes the date and the geographic coordinates of fire ignition and a field estimate of the burned area.The ignition points of the most recent fires (from 2006 onwards) are usually accurately geolocated with GPS, while the positional accuracy of the older ignition points is more uncertain.
For analyzing the relationship between fires, land cover and road network, we used a road network map containing all state, regional and county roads of Sardinia, together with a LULC map of the island at scale 1:25,000.The LULC map was built by the Regional Administration of Sardinia according to the standard classification methodology of the CORINE Land Cover legend [36] based on 2003-2004 digital orthophotos and high-resolution Ikonos satellite images (2005-2006) combined with field surveys.Both maps were downloaded from the Regional geographic information portal of Sardinia (www.sardegnageoportale.it).
We next selected eight land cover classes with a sufficient number of fire records for studying fire incidence patterns at the landscape scale: (i) discontinuous urban areas, (ii) arable land, (iii) mixed agriculture, (iv) olive groves, (v) deciduous forests, (vi) maquis, (vii) garrigue, (viii) pastures.The selected classes were grouped into two main categories with different levels of anthropogenic pressure according to the CORINE Land Cover classification [36]: Artificial and agricultural classes (i-iv) and Natural and semi-natural classes (v-viii).

Data Analysis
We first created a sequence of overlapping buffers around the road network, ranging from 0 to 200 m in four steps of 50 m each (i.e., 0-50 m, 0-100 m, 0-150 m and 0-200 m).The buffers were then intersected with the selected LULCs to calculate the area A kb of each land cover class k within each buffer b.Next, the number of fires N kb within each land cover class in each buffer were counted (Table 1).The actual number of fires N kb was then compared with the expected number of fires assuming complete spatial randomness in fire distribution within each land cover class (CSR k , implying that distance to roads does not influence the spatial distribution of fires within single LULC classes) and within each buffer (CSR b , implying that distance to roads does not influence the relative proportion of fires among different LULC classes).For land cover class k, the number of fires in a given buffer b expected under CSR k is obtained by multiplying the total number of fires in that land cover class N k+ by the proportional area of that class within the buffer relative to the total area of k over the entire study area (p k b ).The expected number of fires E k b in class k within buffer b is then Likewise, for a given land cover class k, the number of fires in a given buffer b expected under CSR b is obtained by multiplying the total number of fires in that buffer N +b by the proportional area of that class within the buffer relative to the total area of buffer b (p kb ).The expected number of fires E kb in class k under CSR b is then To compare the actual number of fires N kb in class k within buffer b with the corresponding expected values E k b and E kb , we finally calculated the normalized selectivity indices And σ k b thus summarizes the degree of fire clustering of a given land cover class k in a given buffer b compared to the fire distribution within the whole land cover class, while σ kb summarizes the degree of fire clustering of class k in buffer b compared to the fire distribution within the whole buffer.Both indices range from −1 to +1.Values of σ > 0 represent areas where fire occurs more frequently than expected under complete spatial randomness in fire distribution within each land cover class (σ k b ) or within each buffer (σ kb ), meaning that these areas are positively selected by fire.By contrast, values of σ < 0 represent areas where ignitions are less frequent than expected under complete spatial randomness, meaning that these areas are negatively selected by fires [18].

Results
The selectivity values of σ k b and σ kb for all buffers and land cover classes are shown in Table 2, while Figure 2 shows the plot of σ k b vs. σ kb for the 100-m buffer only.Note that, as shown in Table 2, the results obtained for the 100-m buffer are qualitatively and quantitatively consistent with the results obtained for the other buffers.
According to Figure 2, the values of σ k b show that relative fire occurrence close to roads is higher than expected by chance alone in all LULC classes, meaning that for every land cover class, the proportion of fires close to roads is significantly higher compared to with respect to more interior areas of the same class.Nonetheless, this 'clustering effect' tends to be much lower for classes with high human pressure (i.e., urban areas, arable land, mixed agricultural land, and olive groves) than for natural and semi-natural classes (i.e., pastures, deciduous forests, garrigue, maquis) where the occurrence of human-related ignitions highly depends on the proximity to accessibility networks.
At the same time, the values of σ kb show that land uses with high anthropogenic pressure experience a disproportionately higher number of ignitions than expected under complete spatial randomness in fire distribution within each buffer (σ kb > 0).By contrast, in forests and semi-natural LULC classes, the number of ignitions is generally lower than expected by chance alone (σ kb < 0).As a result, all natural and semi-natural land uses are located in the upper-left portion of the plot of σ k b vs. σ kb , which is characterized by positive values of σ k b and negative values of σ kb , whereas all anthropogenic land uses are located in the upper-right quadrant in which both variables assume positive values.

Discussion
There has been a long research tradition focusing on the socio-economic aspects of wildfire [17,30,33,[37][38][39][40].In this study, we analyzed existing interactions between roads and land cover in driving fire ignition.Despite existing uncertainties in terms of the geolocation of a subset of ignition points, the very consistent patterns of the influence of roads and LULC classes on fire ignitions obtained throughout the buffer sizes allow us to assume a strong consistency of the overall results.According to our results, fires are preferentially ignited close to roads in all LULC classes.This effect is consistent with a number of previous studies in various ecosystems, which showed a high spatial clustering of human-caused fire ignitions around accessibility networks [8,9,22,25,[41][42][43].However, the observed influence of roads on the spatial pattern of fire ignitions is not the same in all LULC classes.It is particularly noticeable for the natural and semi-natural land uses, which usually experience a lower number of fires compared to more human-impacted land uses and for which roads represent the major vector of anthropogenic ignition energy [29].By contrast, in LULC classes with high anthropogenic pressure, the ubiquity of the ignition energy related to the constant and widespread presence of humans reduces the clustering effect of the road network on fire ignitions.
In LULC classes with coarser fuels, such as deciduous forests or maquis, road corridors may directly influence ignition risk by reducing canopy closure, thus affecting the environmental conditions along road edges through increased solar radiation, temperature and wind speed, and reduced fuel moisture content [25].These edge effects are not limited to the immediate neighborhoods of roads, but extend over a wider area in which fuel conditions are altered compared to continuous forest sites.The peculiar fuel characteristics of road edges have led [25]) to consider them as a separate and specific landscape component in terms of fuel proprieties that differ from the adjacent interior areas.
These results are of relevance for land and forest managers who have to be aware that planned roads and extensions of the road network may have different impacts on the potential of fire ignitions according to the land uses in the concerned area.In this sense, the increased ignition risk related to

Discussion
There has been a long research tradition focusing on the socio-economic aspects of wildfire [17,30,33,[37][38][39][40].In this study, we analyzed existing interactions between roads and land cover in driving fire ignition.Despite existing uncertainties in terms of the geolocation of a subset of ignition points, the very consistent patterns of the influence of roads and LULC classes on fire ignitions obtained throughout the buffer sizes allow us to assume a strong consistency of the overall results.According to our results, fires are preferentially ignited close to roads in all LULC classes.This effect is consistent with a number of previous studies in various ecosystems, which showed a high spatial clustering of human-caused fire ignitions around accessibility networks [8,9,22,25,[41][42][43].However, the observed influence of roads on the spatial pattern of fire ignitions is not the same in all LULC classes.It is particularly noticeable for the natural and semi-natural land uses, which usually experience a lower number of fires compared to more human-impacted land uses and for which roads represent the major vector of anthropogenic ignition energy [29].By contrast, in LULC classes with high anthropogenic pressure, the ubiquity of the ignition energy related to the constant and widespread presence of humans reduces the clustering effect of the road network on fire ignitions.
In LULC classes with coarser fuels, such as deciduous forests or maquis, road corridors may directly influence ignition risk by reducing canopy closure, thus affecting the environmental conditions along road edges through increased solar radiation, temperature and wind speed, and reduced fuel moisture content [25].These edge effects are not limited to the immediate neighborhoods of roads, but extend over a wider area in which fuel conditions are altered compared to continuous forest sites.The peculiar fuel characteristics of road edges have led [25]) to consider them as a separate and specific landscape component in terms of fuel proprieties that differ from the adjacent interior areas.
These results are of relevance for land and forest managers who have to be aware that planned roads and extensions of the road network may have different impacts on the potential of fire ignitions according to the land uses in the concerned area.In this sense, the increased ignition risk related to increasing housing and road density as postulated by [44] may be additionally modulated by the vegetation types crossed by roads.
To conclude, the proposed approach represents a first step towards the integration and the detailed assessment of two major fire ignition drivers, such as land cover and road network, which represent a basic input to fire risk analysis and modern fire management planning [45,46].A thorough knowledge of the joint effect of LULC and roads on fire ignition patterns is a key point, yet future studies should also take into account fire size and severity in order to fully understand the role of anthropogenic presence in fire regimes.

Figure 1 .
Figure 1.The study area of Sardinia with the road network used in this paper and the fire ignitions recorded during 2000-2015.

Figure 1 .
Figure 1.The study area of Sardinia with the road network used in this paper and the fire ignitions recorded during 2000-2015.

Figure 2 .
Figure 2. Plot of   ′  (fire selectivity related to land cover) vs. (  ′ fire selectivity related to buffer size) for the 100-m buffer.Gray symbols represent artificial and agricultural Classes, open symbols represent natural and semi-natural classes.

Figure 2 .
Figure 2. Plot of σ k b (fire selectivity related to land cover) vs. (σ kb fire selectivity related to buffer size) for the 100-m buffer.Gray symbols represent artificial and agricultural Classes, open symbols represent natural and semi-natural classes.

Table 1 .
Number of fires and areal extent of the selected land cover classes within each buffer.

Table 2 .
Values of the selectivity indices σ k b and σ kb for the selected land cover classes within each buffer.