Next Article in Journal
Spatial Simulation of Land-Use Development of Feixi County, China, Based on Optimized Productive–Living–Ecological Functions
Previous Article in Journal
Influence of Miscanthus Rhizome Pyrolysis Operating Conditions on Products Properties
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

GIS-Based Modeling for Vegetated Land Fire Prediction in Qaradagh Area, Kurdistan Region, Iraq

by
Sarkawt G. Salar
1,*,
Arsalan Ahmed Othman
2,3,*,
Sabri Rasooli
4,
Salahalddin S. Ali
5,
Zaid T. Al-Attar
6 and
Veraldo Liesenberg
7
1
Department of Geography, College of Education, University of Garmian, Sulaymaniyah 46021, Iraq
2
Iraq Geological Survey, Al-Andalus Square, Baghdad 10068, Iraq
3
Department of Petroleum, College of Engineering, Komar University of Science and Technology, Sulaimaniyah 46013, Iraq
4
Department of Forestry, Faculty of Natural Resources, University of Guilan, Someh Sara 41996-13776, Iran
5
Civil Engineering Department, College of Engineering, Komar University of Science and Technology, Sulaimaniyah 46013, Iraq
6
Department of Geology, University of Baghdad, Al-Jadiryah Street, Baghdad 10071, Iraq
7
Department of Forest Engineering, Santa Catarina State University (UDESC), Lages 88520-000, SC, Brazil
*
Authors to whom correspondence should be addressed.
Sustainability 2022, 14(10), 6194; https://0-doi-org.brum.beds.ac.uk/10.3390/su14106194
Submission received: 2 April 2022 / Revised: 10 May 2022 / Accepted: 17 May 2022 / Published: 19 May 2022

Abstract

:
This study aims to estimate the susceptibility of fire occurrence in the Qaradagh area of the Iraqi Kurdistan Region, by examining 16 predictive factors. We selected these predictive factors, dependent on analyzing and performing a comprehensive review of about 57 papers related to fire susceptibility. These papers investigate areas with similar environmental conditions to the arid environments as our study area. The 16 factors affecting the fire occurrence are Normalized Difference Vegetation Index (NDVI), slope gradient, slope aspect, elevation, Topographic Wetness Index (TWI), Topographic Position Index (TPI), distance to roads, distance to rivers, distance to villages, distance to farmland, geology, wind speed, relative humidity, annual temperature, annual precipitation, and Land Use and Land Cover (LULC). To extract fires that occurred between 2015 and 2020, 121 scenes of satellite images (most of them are scenes of Sentinel-2) were used, with the aid of a field survey. In total, 80% of the data (185,394 pixels) were used for the training dataset in the model, and 20% of the data (46,348 pixels) were used for the validation dataset. Conversely, 20% of these data were used for the training dataset in the model, and 80% of the data were used for the validation dataset to check the model’s overfitting. We used the logistic regression model to analyze the multi-data sites obtained from the 16 predictive factors, to predict the forest and vegetated lands that suffer from fire. The receiver operating characteristic (ROC) curve and the area under the curve (AUC) were used to evaluate the accuracy of the proposed models. The AUC value is more than 84.85% in all groups, which shows very high accuracy for both the model and the factors selected for preparing fire zoning maps in the studied area. According to the factor weight results, classes of LULC and wind speed gained the highest weight among all groups. This paper emphasizes that the used approach is useful for monitoring shrubland, grassland, and cropland fires in other similar areas, which are located in the Mediterranean climate zone. Besides, the model can be applied in other regions, taking the local influencing factors into consideration, which contribute to forest fire mitigation and prevention planning. Hence, the mentioned results can be applied to primary warning, fire suppression resource planning, and allocation work. The mentioned results can be used as prior warnings of the outbreak of fires, taking the necessary measures and methods to prevent and extinguish fires.

1. Introduction

Vegetated lands are important geomorphological phenomena. They play a significant role in surface processes, such as erosion and landslides [1,2,3]. Vegetation is a vital natural resource, which aids in the preservation of environmental balance [4]. Wildfires are the greatest serious hazard to the world’s grasslands and forests, causing severe ecological, economic, and social consequences [5]. Wildfires may occur as a result of long dry seasons; however, the majority of these events are caused by anthropogenic actions, such as accumulated fuel load in unmanaged vegetated environments as well as vegetation removal and clearing for increased pasture and agricultural lands [5,6]. They have not, yet, been described adequately because of the complexity of the processes determining fire behavior, the great amount of data required, and the difficulty of the extraction and collection of the data [7]. Every year, wildfires cause a great deal of damage to natural vegetation, property, the biotic elements of a given ecosystem, and the environment itself, thereby threatening people’s lives and significant human resources [8]. Fire leads to environmental degradation, by reducing vegetation coverage and, consequently, food resources for fauna, while modifying soil properties leading to degradation, mass movement, and diversity exhaustion [9].
It is vital to identify the factors affecting the spread and occurrence of fires in vegetated areas, to reduce the fire risk in a given environment [10]. Generally, the causes of fires in the vegetation include natural (environmental) and artificial (anthropogenic) factors. Environmental factors affecting fire include biological, physiographic, and climatic factors [11]. Vegetation is one of the most common biological factors affecting fire occurrence, which exerts its effects depending on the species, density, humus, and moisture [12,13,14,15,16]. The factor of topography is composed of the sub-factors of slope, aspect, altitude, plan curvature, and Topographic Position Index (TPI). TPI is another vital factor affecting fire occurrences [17,18,19,20,21].
Among climatic factors, temperature increases, relative humidity decreases, precipitation amounts, and wind speed and direction affect fire occurrences [19,22,23]. Aside from natural factors, anthropogenic factors as fire ignition sources have caused many fires worldwide [24,25], such as roads, farmlands, and residential areas located within or around forests [14,26,27]. In addition, many fires are produced as a means of clearing lands for agriculture or afforestation purposes [28,29,30]. Therefore, we carried out an inventory of previous research that has a similar topographic environment to the study area, to select the effective factors in fire occurrence. These factors can be divided into two categories: environmental and anthropogenic factors. Environmental factors include slope, slope aspect, elevation, topographic Wetness Index (TWI), TPI, annual temperature, precipitation, wind speed, Normalized Difference Vegetation Index (NDVI), distance from streams, relative humidity, and geology. On the other hand, anthropogenic factors are the distance from roads, the distance from settlements, the distance from farmlands, and LULC.
Various studies have been done on the factors affecting fire intensity and occurrence [2,5,31,32,33,34,35,36]. The geographical distribution of fires in central Spain has been studied, in terms of environmental factors. This study showed a significant relationship between the area of the burned regions and the topography of those regions. In the same vein, Lozano et al. [37] showed a high correlation between fire occurrence and environmental factors on different geographical scales. Zumbrunnen et al. [38] showed that climate, roads, and livestock play a significant role in fire occurrence in Switzerland. In addition, a non-linear relationship was established between fire occurrence and anthropogenic factors, with the fire occurrence rate being higher in warmer climates. In another study led by Justus [39], the vital factor affecting forest fire distribution is topography, followed by anthropogenic factors.
Humans cannot control fire completely, but they can reduce its intensity and the damage to natural resources, to a great extent. To this end, each area can be zoned, given the points available in terms of fire likelihood. In other words, areas with high fire likelihood must be identified; then, necessary measures should be adopted to prevent or fight fires [40,41]. Several studies have been conducted, so far, to model and assess fire risks in the world, using different methods [42,43,44,45]. These studies have used the analytical hierarchical process (AHP) [2,5,11,13,43,44,45], fuzzy methods [21,46,47,48,49], Analytical Network Process (ANP) [50], logistic regression [51,52,53,54], artificial neural networks [20,55], support vector machine c, and random forests [15,56] to model fire risks. These methods require accurate preparation and implementation of the produced models, for the factors affecting fire occurrence. Given the complexity of the fire processes in various regions and the effects of various factors on its occurrence, each of the mentioned methods has potential and limitations. Moreover, in many cases and places, researchers may not have the necessary efficiency for the accurate prediction and geographical zoning of fire risks [57]. Therefore, examining the capabilities of these methods and employing them in different regions can help with choosing the best and most effective method. Moreover, the use of the capabilities of the Geographic Information System (GIS) and remote sensing has facilitated the zoning and modeling of fire occurrence [2,12,14,47,58].
In recent years, frequent fires have occurred in the vegetated lands of the Iraqi Kurdistan Region. In fact, fires are among the major causes that destroy these vegetated areas, and few studies have been reported so far. Thus, given the limited natural resources, including shrublands, in this region, the importance of trying to protect them against fires is maximized. However, despite the occurrence of frequent fires, no research has been conducted on the factors affecting the occurrence, intensity, and distribution of fires in this region. Accordingly, in this study, the shrublands of the Qaradagh area were studied, as a sample of the Kurdistan Region forests. Our study investigates the effects of natural and anthropogenic factors on fire risks, by integrating satellite images, logistic regression, and GIS techniques. Therefore, fire risk potentials were mapped and modeled. Preparing the map of fire potentials using a suitable method plays an important role in understanding and predicting fire occurrence in vegetated areas of this region, thereby helping environmental agencies, farmers, and land managers to prevent and deal with this problem through better supervision of their natural resources.

2. Materials and Methods

2.1. Study Area

The Qaradagh area encompasses the southern part of Sulaymaniyah governate, in the Kurdistan Region of Iraq (Figure 1). It is an important area in terms of ecology, biodiversity, and recreation, extending along with the two-mountain series from the northeast to the southeast. The Qaradagh, Sagrma, and Gulan mountain series (anticlines) forms the southern and western parts of the region, while Baranan Mountain (anticline) forms the northern and eastern parts. The region’s altitude ranges between 337 m and 1853 m above sea level. This region has a Mediterranean climate [59,60]. According to data from the Tropical Precipitation Measuring Mission (TRMM), for the period of 1998 to 2017 [61], the average annual precipitation of the study area varies between 591 and 703 mm; whereas, the average annual temperature varies between 17 and 19 °C, for the period of 1979 to 2014 [62].
Fires constitute one of the most important challenges that threaten the ecosystem and biodiversity in the Qaradagh area. Through the available satellite images of fires (between May 2015 and October 2020), we surveyed the area and recorded (660) vegetation fires (Appendix B). As a result, 12.35% of the whole study area was exposed to fire during the mentioned period. Sagrma Mountain, at the western part of the study area, was exposed to fire five times, whereas Baranan Mountain, at the eastern part of the study area, was subjected to fire four times.

2.2. Material and Factors Affecting Fire Susceptibility

Many natural and socioeconomic factors affect fire susceptibility in vegetated and pastures of a region, with the impact of each of these factors being specific and different in each region [49]. However, due to the large number and high variety of the factors affecting fire occurrences and data availability, and, also, to facilitate the better modeling and understanding of the calculation process, part of these factors—the main factors affecting fire occurrences in shrublands and pastures—are selected. In the present study, the factors were selected after the comprehensive review of 57 papers (Appendix A) selected from the Scopus (Elsevier) database and related to fire susceptibility in forests and pastures in different countries, published between 2008 and 2020. Accordingly, 16 predictive factors, confirmed in most sources, were employed to produce a fire susceptibility map.
In general, the factors effective in fire occurrence in this study were divided into two categories: environmental and anthropogenic factors. Environmental factors include slope, slope aspect, elevation, Topographic Wetness Index (TWI), TPI, annual temperature, annual precipitation, wind speed, NDVI, distance from streams, relative humidity, and geology. On the other hand, anthropogenic factors are the distance from roads, the distance from settlements, the distance from farmlands, and LULC. From among the aforementioned factors, the three discrete factors, i.e., slope aspect, LULC, and geology, are classified into 9, 5, and 10 classes, respectively. Thus, each class in these three factors is considered a variable, and the total number of layers was 37 factors.
We selected the factors used in the previous papers at the rate of ≥10% (16 factors for this study). In other words, we eliminated factors decreasing the accuracy of vegetation fire susceptibility mapping and retained those increasing it.
The required and available information were employed, in light of the factors impacting vegetation fires that were discovered and mentioned (Figure 2). In total, 121 pieces of Sentinel-2 images were downloaded from the EO-Browser website, to map burned areas at their burning time, during the period 2015–2020. The Sentinel-2 imagery was utilized to determine the NDVI, using Sentinel Application Platform (SNAP) software (Table 1). NDVI is defined based on Equation (1) [63]:
NDVI = NIR     R NIR   +   R
where NIR and R values are the infrared and red portion of the electromagnetic spectrum, respectively.
In addition, the Digital Elevation Model (DEM) of the Shuttle Radar Topography Mission (STRM), with 1 arc-second resolution, was employed to derive slope, slope aspect, elevation, TWI, and TPI maps. The TWI and TPI are defined, based on Equation (2) [64] and Equation (3) [65], respectively:
TWI = ln ( α tan β )
TPI = E c ( 1 nM i m E i )
where α is the area of the drained zone, and tan β is the slope angle in degrees. Where Ec is the elevation of the central pixel, Ei is the elevation of the grid within the local window, and n is the total number of surrounding points employed in the evaluation.
Moreover, GIS world imagery was used to create the road, settlements, and farmlands maps, to extract the maps of distance from roads, distance from settlements, and distance from farmlands. More details of the datasets used and their specifications are fully described and presented in Table 2. Figure 3 and Figure 4 show the 16 predictive factors, which were used in this study.
Due to the low number of climatological stations in the study area, we used the TRMM monthly global precipitation data, from 1998 to 2017 [61], and air temperatures acquired from Climate Forecast System Reanalysis (CFSR), from 1979 to July 2014 [62]. The TRMM data within Iraq and the Kurdistan Region have been used and verified by several researchers, such as [75,76], where their accuracy shows a strong direct relationship with the climatological stations’ data (coefficient of determination (R2) > 0.8). Although the CFSR was not verified before, it was compared to the data obtained from Darbandikhan station, located southeast of the study area (Figure 4A). For this purpose, we used the nearest pixel that was located to the climatological stations. In Figure 4C,D, the R2 indicates a strong direct relationship (>0.8) between the CFSR and the climatological stations’ data. The spatial resolution of the CFSR is ~38 km [77]. The TRMM data have a resolution of 0.1 degrees (~9 km). The TRMM and the CFSR data have been converted to vectors and interpolated with a 30 m spatial resolution.
The geological map was obtained from the Iraqi Geological Survey [78]. Whereas, the LULC and geological maps were derived from USGS EarthExplorer website and GEOSURV-Iraq [79].

2.3. Burn Scar Inventory

A map of past fires is required to investigate the relationship between fire locations and the factors influencing fire occurrence. Thus, the manual delineation of the burn scar inventory was created for 121 satellite images, for the period from 2015 to 2020, in addition to the field studies. These data have been delivered from Sentinel Hub services, which, efficiently, provide long-term analysis. Next, the locations of all burn scars, which occurred during this time interval, were observed on these images, using RGB:432 color composite. Finally, the burn scars were drawn as polygons, and were, thereafter, captured and used for modeling and validation purposes.
When all required data and information were collected from different sources, and after necessary corrections and processing were applied, all factors were converted to the raster format. Moreover, given the fact that the input layers have varied units, which are not suitable as direct inputs for the logistic regression, the input parameters were normalized within the range of 0 to 1 [80]. In the end, making use of the models and coefficients achieved from the logistic regression method, the maps were overlaid with fire risk maps of the study area produced.

2.4. Data Analysis by Logistic Regression

After the parameters were prepared, fire susceptibility maps were prepared, using the logistic regression method. Logistic regression is a statistical method devoted to the group of generalized linear statistical models, which predicts the probability of an event (fires) using independent variables. The key point in logistic regression is that the dependent variable is Boolean (a two-state variable that can be only number 0, meaning non-occurrence, or number 1, indicating the occurrence of an event). Concerning fire probability mapping, logistic regression is used to find the best model for describing the relationships between the presence or absence of a dependent variable (vegetation fires) and a set of independent variables. Logistic regression utilizes the maximum likelihood estimation (MLE) method to find the set of parameters fitting the model best. The model’s output yields coefficients between 0 and 1, which through fuzzy theory assigns value 1 to probabilities more than 0.5 (fire occurrence) and value 0 to probabilities less than 0.5 (non-fire occurrence). Finally, it produces a Boolean map for fire susceptibility.
We used the fitting generalized linear models (glm), with a binomial family and logit link in R. Logistic regression is used, by assuming that the probability of the dependent variable being 1 follows the logarithmic curve. Besides, its value is estimated by Function 4 [81], as follows:
P ( y = 1 X ) = exp ( SBX ) 1 + exp ( SBX )
where, P is the possibility of being 1 for the dependent variables, X represents the independent variables, B is the estimated parameters, and y represents the dependent variables. To linearize the function above, the following transformation is applied:
l o g e ( P 1 + P ) = b 0 + b 1 x 1 + b 2 x 2 + b 3 x 3 + b k x k + e r r o r   t e r m
where, P is the possibility of being 1 for the dependent variables, x1, x2 …, xk represents the independent variables, b0 is the coefficient of the regression equation, and b1, b2, …, bk are the coefficients of each independent variable.
This logarithmic change makes the predicted probability become continuous within the range 0 to 1 and makes the model’s output be presented as a geographical prediction map for risk probability [82].

2.5. Preparation of Training and Validation Dataset

In total, 103 scenes of satellite images were taken from Sentinel-2, 6 scenes from Landsat-8, and 3 scenes from Landsat-7 (Appendix B). The Landsat images were used to avoid cloudy days and cover the time before the Sentinel-2 launch in 2015. All burned locations, during the mentioned period, were merged to show all areas exposed to the fire on one map. We separated the burned area randomly into training and validation data subsets, as recommended in the literature [83]. A percentage of 80% (185,394 points) was used for the training dataset in the model and 20% (46,348 points) for the validation dataset. A similar number to the training points (185,394 points) from the burned areas were picked at random from the unburned areas. Hence, 185,394 points from the burned areas and 185,394 points from the unburned areas were entered into the model as a training dataset.

2.6. Omission of the Factors Less Effective in Fire Occurrence to Map Fire-Susceptible Areas

To map the fire-susceptible areas, all the training data collected from the 16 main factors were entered as a table into RStudio. The table includes a column that shows the case of the burn scar that exists. An existed burn was coded as 1 and a non-existed burn was coded as 0. The logistic regression model was employed to analyze the data and find the estimation constants (α and β) for probability calculation.
Accordingly, it was observed that in most of the studies, factors, such as slope, slope aspect, and altitude, were used. However, factors such as relative humidity, distance from agricultural lands, and TPI (less than 20% of the reviewed articles) were utilized in fewer studies. Apart from a map, in which all of the mentioned factors were used, we prepared other maps, from which we removed one of the less important factors from the total factors every time. Based on their number of frequencies in the studies, factors of less importance included geology, TWI, relative humidity, distance from farmlands, and TPI. Besides, another map was prepared, from which all six factors were removed. Less important factors were omitted, to stress the importance of the impact of the factors examined on preparing the fire susceptibility map. This is because the careful selection of factors is more important than the methods used for the purpose of zoning [75]. In the end, we will compare the aforementioned maps.

2.7. Validation

Several performance metrics, such as the Roc Curve, Mean Absolute Error (MAE), Relative Error (RE), and Percentage Error (PE) have been used to assess the predictive accuracy of the implemented model in the current study.
Validation or accurate evaluation is required in modeling, to determine the scientific significance or evaluation of a study [3,84]. The efficiency of the model presented in the present study was evaluated, making use of the receiver operating characteristic (ROC) curve and the area under the curve (AUC). This is the most common method used for evaluating modeling efficiency, being used in various research fields [85]. The ROC method is the graphical representation of the balance between the False Positive Rate and the True Positive Rate, for every value of probability. The area beneath the AUC-ROC curve indicates the value predicted by the model, by describing its capability of accurately estimating events or their non-occurrence. In general, an AUC value of 1–0.8 indicates very high efficiency, 0.7–0.8 indicates high efficiency, 0.6–0.7 denotes medium efficiency, and 0.5–0.6 indicates low efficiency, of the model.
In this study, 80% of the fires that have occurred in the study period were used as the training dataset and assigned to model production (Appendix C), while the remaining 20% of the dataset was assigned as a validation. In contrast, we used 20% of the dataset to predict the models (Appendix D), while the remaining 80% of the dataset was utilized to validate the results, to check the model’s overfitting.
We used MAE for evaluating the performance of a regression model (Equation (6)). It is a measurable indicator to validate the closeness of predicted values to actual values, as the name suggests [86]. The RE was calculated (Equation (7)), where the ratio of the absolute value to the approximate value is relative error. If the approximate error is unknown, the actual value is used [87]. Equation (8) used the reference value to calculate the PE [88].
MAE = 1 n i = 1 n | P O |
RE = | P O | O
PE = 100 * | P O | O
where n represents the represents the number of observations; O represents the actual value; and P represents the predicted value from the susceptibility model.

3. Results

3.1. Vegetation Fire Susceptibility Mapping

After preparing all the factors required and determining the geographical locations of the fires, which took place during the six years within the study area, logistic regression was utilized to determine the relationship between factors affecting fire occurrence and the related spatial modeling. The logistic regression relationship was established between fire occurrences, as dependent variables, and the mentioned parameters, as independent variables. Since the studied factors were selected from the literature review, the regression model was repeated seven times, once to validate the dataset and once, again, to train the dataset to examine the importance of selecting these factors (Appendix C and Appendix D). The p-values for all predictor factors used in all models are less than 0.05. Therefore, they have a statistically significant relationship with the fire occurrence in all models. In group 1, we selected those factors that were used in less than 20% of the studies in various regions, which were entered into the regression. In group 2, geology was omitted. Besides, TPI, TWI, and the distance from agricultural lands were removed in groups 3, 4, and 5, respectively. In group 6, relative humidity was omitted, and, finally, all respective factors were entered into the logistic regression in group 7 (Appendix C and Appendix D). As a result, seven zoning maps of fire susceptible areas were prepared. The importance of the order of the factors, based on the coefficients achieved in each of the seven cases, is according to the order of the following diagrams, with the maps indicating seven zoning times (Figure 5 and Figure 6).

3.2. Validation and Comparison of the Maps

Figure 7 shows the validation results for the fire potential map. This diagram shows that the success rate in all groups is quite apt, ranging from 84.85 to 97.68%. Among all groups, group 3, i.e., the group with TPI excluded from the model calculations, has the highest stability and success rates (AUC), which are 97.30% and 97.31% for both 20% and 80%, as training dataset evaluations, respectively (Figure 7). In addition, when assessing the 20% as a training dataset, group 7 had the highest accuracy (97.68%), compared to the findings of the other groups (Figure 7A), and it came in second place (96.72%) in terms of accuracy, when evaluating the 80% as a validation dataset (Figure 7B), which means less stability than that of group 3.
The results of the MAE, RE, and the PE of the seven groups show that the best model is group 7. It has a lesser error, while MAE, RE, and PE are 0.62, 2.32, and 0.64, respectively. The worst model (highest errors) come in group 4, where the MAE, RE, and PE are 2.54, 15.15, and 2.99, respectively (Table 3).

3.3. Key Factor

The predictive ranges of all seven models show that LULC classes, which are related to vegetated areas, are the most common factors that affect fire occurrence (Appendix C and Appendix D). This factor includes six classes, which, in order of weight importance, include: cropland, grassland, urban and built-up land, open shrub, and barren area. The wind speed factor has the second-most impact on vegetation-fire occurrence, followed by the geological formations factor, which is composed of clastic sediments (i.e., Tanjero, Aqra, and Bekhme formations).

4. Discussion

4.1. Model Comparison

This study aimed to examine the factors affecting fire occurrence by using the logistic regression method to zone fire occurrence in the Qaradagh area of the Iraqi Kurdistan Region. Based on the AUC of the seven models of the logistic regression (Figure 7), group 3, i.e., the group in which the factor of TPI was not considered, had the highest success stability rate (97.30% and 97.31%) for both evaluations. This reveals that it is not necessary to include the TPI factor in the applied models, to find the susceptible areas to fires. This group ranks second, behind group 7, based on the accuracy results from the MAE, RE, and PE methods, which indicates the effectiveness of all the selected criteria in determining the sensitive areas to fires, as confirmed by the outputs of Figure 7A (97.68%). After groups 3 and 7, the accuracy of all used methods revealed that the values of groups 6 and 5 (excluding relative humidity and distance from agricultural land) were stable, indicating that the effect of the factors on vegetation fires was maintained at the same level and sequence. In terms of stability and sequence, groups 2 and 1 were similar to groups 6 and 5, although with less accuracy, as shown in Figure 7 and Table 3. As for groups 2 and 1 (group 2 with the geology factor removed, and group 1 with all factors that were used in less than 20% of past case studies removed), they were similar to groups 6 and 5 in terms of stability and sequence, but with less accuracy. The low level of accuracy is due to the importance of the geological factor, and its impact on the spatial distribution of plants that are exposed to fires. The accuracy of all the applied statistical methods shows that the least accurate model was represented by group 4, from which the TWI factor was removed. This highlights the effectiveness and strength of the impact of this factor on the accuracy level of the model, which must take into account the extent of its importance. As demonstrated, the values of the statistical accuracy methods show the very strong performance of the models, as well as the selection of apt factors, for preparing a fire zoning map in the studied area. However, it is crucial to be careful enough, when selecting factors affecting fire occurrence. The low success rate in group 4 implies that the factors being more important in fire zoning should not be removed to increase the model’s success rate. This methodology is also applicable to other areas with similar geographical and environmental conditions, and decision-makers could benefit from such an approach.

4.2. Factors Impacting Vegetation Fire

The spatial relation strength of affecting and impacting factors on the event was shown by the logistic regression model. It determined the relative importance of each factor influencing the vegetated fire in the study area. Hence, it revealed that the increase in wind speed, NDVI, and geological formation, which includes vegetation, as well as LULC class, which is related to vegetated areas, are the most common factors affecting fire occurrence (Appendix C and Appendix D). According to the factor weight results (Figure 5, Figure 6 and Figure 7), classes of LULC and, then, wind speed gained the highest weight among all groups, with their weight being significantly different from that of other factors. Among LULC classes, the highest weight has been assigned to the cropland class, which can be considered the key factor that enables us to better anticipate and, even, prevent land fires, requiring, therefore, some proper land management practices such as reducing fuel loads periodically and, also, installing firebreaks. Similarly, the results of Pourtaghi et al. [34] show that LULC is the most effective factor in fire occurrence. The presence of adequate fuel is the main cause of most fires occurring in pastures, forests, and agricultural lands. Apart from natural fires in these areas, the role of anthropogenic factors trying to eliminate agricultural residues, repel destructive pests from agricultural products, and speed up the revival of grass species in pastures by creating deliberate fires, is noteworthy. Besides, some unofficial reports relate many wildfires to ethnic disputes, which need more in-depth investigation. Given these reasons, and considering that the highest weight in all seven groups was related to the agricultural land class, one can say that humans are the major cause of fires in this particular region.
After fire occurrence, natural factors can positively or negatively affect fire spread. For example, wind speed is one of the most effective factors in its spread. Besides, wind affects both the fire triangle of fuel, heat, and oxygen as well as the fire behavior triangle of fuel, topography, and climate [89]. Wind, in the fire triangle, transfers more oxygen to the fuel, thereby causing the initial fire spark. On the other hand, the climate factor is seen in the fire behavior triangle, with wind being one of the main components of climate. Wind spreads fire further to surrounding areas and makes its extinguishment more difficult.
The data about wind direction in the study area are poor, as they are missing detailed information. Local winds dominate topographic fires more than prevailing winds, due to slope and differences in the solar heating of the Earth’s surface. Forest fires are more sensitive to small changes, so slight differences in topographic winds, slopes, or sides have a greater impact on fire behavior [90]. From the general information, the effect of wind direction on vegetation fires in each pixel is close to equal, in areas characterized by an almost equal geographical distribution of vegetation density; therefore, we ignored the wind direction factor. In the future, we will consider the wind direction factor once local measurements become available for better spatial analysis.
High temperatures, as well as low precipitation and relative humidity, dry the environment, thereby providing suitable environmental conditions for fire occurrence and spread. In regions with these attributes, in case there are combustible materials with, usually, a small amount of humidity, fire occurrence potentials are very high [91]. Due to the fact that two climatic factors, temperature and precipitation, are highly effective in causing fires in terms of time, and given the investigation that was made in the studied area, they cannot cause significant climatic changes in a habitat on such a limited scale. Research on the correlation between fires and soil types indicates that the soil moisture content and its effects on the diversity of plant species are the most effective factors in controlling fire occurrence or non-occurrence [34].
The factors of the distance from roads, residential areas, and rivers reflect humans’ role in fire occurrence [92]. These factors can play a dual role in fire occurrence and spread. In the vicinity of these factors, due to more human traffic and activity, fire occurrence is more likely. However, the presence of these factors can be effective in preventing fire spread, after its occurrence. Roads can act as firebreaks, by preventing the fire from spreading to surrounding areas. As fire can be contained by humans, in case it occurs in the vicinity of residential areas, it will be extinguished sooner due to easy human access. In the present study, locations of residential areas were considered as a polygon, not a point, and, then, a map of the distance from these areas was prepared. In other words, one could say that population density has been somehow considered. Accordingly, the higher the population density is, the higher the fire probability will be. In the vicinity of rivers and waterways, due to the higher humidity of combustibles, the fire spread rate decreases [45].
The TPI implies that flat areas are not as susceptible to fire as sloping areas. Moreover, fire occurrence is higher in convex slopes than in concave slopes. This case can be justified, both at the beginning of the fire and at the time of its spread. Humidity is lower in sloping areas, so the fire spread speed increases with increasing the slope angle after fire occurrence. The TWI determines the effects of topography on the rate of saturated surfaces for producing runoff, which plays a significant role in the soil moisture content and slope stability. Hence, there is an inverse correlation between an increase in this index and fire occurrence and spread [93].
The high density of the Gramineae species results in increased vegetation density and the supply of fuels effective in surface fires. Moreover, fire statistics indicate that the highest frequency of fires occurs almost from early summer to early autumn. Given the increasing temperatures, decreased relative humidity, dryness of vegetation in summer, and concomitant effects of human activities, fire is intensified in the region, being consistent with the studies of Gholami et al. [94] and Rafiee et al. [95].
In the Qaradagh area, further work is necessary to estimate the local biomass and total carbon levels, which, in turn, can give an idea of the health of the available crops that could feed random fires. Therefore, we suggested using the method applied by Oliveira et al. [96] to integrate the field-observed biomass data with lidar data from various sensors, to estimate the biomass and carbon levels obtained in the dry vegetated land within our study area. Moreover, lidar metrics can be used to predict wildfires behaviors and risk assessment measures, by combining lidar data with Landsat-derived indices [96,97] to map canopy cover, canopy height, canopy base height, and canopy bulk density as well as to evaluate the aboveground carbon density through shrubland allometric relationships.
Furthermore, the addition of hyperspectral data can be beneficial in terms of distinguishing the moisture content of crops more accurately [98], which, in turn, can act as a predictor of fires rate. The implementation of such data, in addition to our model, will provide holistic yet easy-to-apply methods for agriculture managers, which will aid them in their production rates while minimizing losses through raising awareness on how to combat wildfires.

4.3. Advantages and Disadvantages

One of the major merits of this study is providing information about the statistical weight of every individual variable on the overall result. Moreover, the high accuracy achieved in the present research could have been due to the use of the majority of the factors affecting fire occurrence, which, in turn, allowed the preparation of a more accurate fire susceptibility map. The model’s limitation is that the factors impacting the vegetation fire differ in space and time; in other words, the model cannot be applied everywhere, at any time, using the same factors influencing the event. Therefore, this model can be applied in arid and semi-arid mountainous regions. For other regions that have different environments and topography, the priority should be on careful choice of the predictive factors, which is far more significant than the methods utilized. Interestingly, multitemporal measurements of Synthetic Apertura Radar (SAR) and Light Detection and Ranging (LIDAR) are suggested for further studies, to estimate carbon release and better assess soil erosion and vegetation recovery.

5. Conclusions

Fire events, caused by natural or anthropogenic factors, are a major threat to landscape resources, so they must be controlled and prevented. To minimize this threat to forests, managers should be fully knowledgeable about past fire processes. Accordingly, spatial analysis of fires must be performed in every region, and conservation strategies should be formulated for every region, after a thorough understanding of temporal–spatial processes. The prerequisite for such endeavor is the ability to make predictions, perform spatial analyses, and prepare a fire susceptibility map. The use of the logistic regression method is one of the risk-modeling methods, which helps to determine the relationship between independent and dependent variables. Additionally, combining remote sensing data with geographic information data and statistical models permits fire managers and personnel to predict ‘where and when’ vegetation fires will occur.
In the current study, a fire susceptibility map was prepared for the Qaradagh district, using the logistic regression method. To examine the importance of the factors selected, the regression model was repeated seven times. Each time, one of the factors used less frequently in previous studies was omitted. Moreover, all of the factors were entered once into the regression for modeling. The AUC value was shown to be higher than 96%, in all groups. This indicates the remarkable performance of the model and the selection of suitable factors, for preparing a fire zoning map in the study area. According to the results, to increase the model’s success rate, factors less effective in fire zoning should be removed. Based on the coefficients obtained, it was found out that the land-use factor exerted the greatest impact on fire occurrence, among all seven cases, in the current study area.
This approach can be applied to other areas having the same geographical conditions. However, this requires selecting the variables and their weights that work for the tested region, based on its environmental conditions.
More detailed land cover classification (particularly classifying other vegetation types and trees) can increase the understanding and accuracy of vegetated land fire prediction and is encouraged for future studies. Finally, we suggest considering the Qaradagh district as a natural protection area, since forest resources are rare and vegetation fires are high in the Kurdistan region.

Author Contributions

Conceptualization, S.G.S., A.A.O. and S.R.; methodology, S.G.S., A.A.O. and S.R.; software, S.G.S., A.A.O. and S.R.; validation, S.G.S., A.A.O., Z.T.A.-A. and S.R.; formal analysis, S.G.S., A.A.O. and S.R.; investigation, S.G.S., A.A.O. and S.R.; resources, S.G.S. and S.R.; data curation, S.G.S., A.A.O. and S.R.; writing—original draft preparation, S.G.S., A.A.O., S.R. and Z.T.A.-A.; writing—review and editing, S.G.S., A.A.O., S.R., S.S.A. and V.L.; visualization, S.G.S., A.A.O., S.R., S.S.A. and V.L. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

We would like to thank the USGS and ESA for their free data, which were used in this study. Moreover, we would like to thank ESA and SEOM (Scientific Exploitation of Operational Missions) for providing SNAP software. V.L. thanks the Brazilian Council for Scientific and Technological Development (CNPq) and the Research and Innovation Support Foundation of Santa Catarina State (FAPESC), for individual support. We are thankful to Ahmed K. Obaid for helping in the discussion.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A. Dependable Factors

ReferencesSlopeSlope AspectElevation or AltitudeDistance from the
Residential or Urban Area
PrecipitationDistance from RoadLand Use and Land CoverVegetation CoverAnnual Air TemperatureWind SpeedNDVIBurned PixelsDistance from StreamsSoil TypePlan CurvatureTopographic Wetness Index (TWI)Relative HumidityDistance from Agriculture LandTPI
[2]****** * *
[4]******* *** **** *
[5]****** ** *
[9] *
[12]**** *
[13]** * * *
[14] * ** * * * *
[15]******* *** ****
[19]******* *** ***** *
[21]*********** ****
[22]*** *** *** ****
[23]***** * *** *
[26]******* * * *
[34]******* *** **** *
[35]*** * ** **
[42]*** * *
[43] *
[45]**** * * *
[48]******* ***
[50]******* ***
[53]* * * * * *
[54]** **** ** *
[56]*** * * *
[58]** * * * *
[99]***
[100] *
[101] *
[102] *
[103] *
[104] *
[105]******* * * * *
[106] *
[107] *
[108]*** ***
[109] ** *
[110] * *
[111] * *
[112] *
[113] *
[114] * *
[115]*** ****** * * * *
[116]**** ***
[117] *
[118]******* ****
[119] * *
[120] *
[121] * * *
[122] *
[123]****** ** * *
[124] * ** *
[125] * * *** *
[126]*** * ***
[127]****** *** ** **
[128] * *
[129] *
[130] *
[131]***** * *** *
Sum31302924232323232218161613101010876
Average
rate
54.2452.5450.8540.6840.6840.6838.9838.9838.9832.2028.8127.1223.7316.9516.9516.9513.5611.8610.17
* The factor existed in the paper.

Appendix B. Recorded Vegetation Fires during the Period (May 2015 to October 2020)

No.YearMonthDayDateImageryArea (m2)No.YearMonthDayDateImageryArea (m2)
1201552121052015Landsat-727,270331201792626092017Sentinel-2762,500
2201552121052015Landsat-72816332201792626092017Sentinel-250,670
3201552121052015Landsat-752,620333201792626092017Sentinel-232,170
4201552121052015Landsat-722,3203342017101111102017Sentinel-224,900
5201552121052015Landsat-763,6803352017111111102017Sentinel-2387,900
6201552121052015Landsat-717,6503362017111111102017Sentinel-21,534,000
7201552121052015Landsat-718,4903372017111111102017Sentinel-23992
8201552121052015Landsat-771653382017111111102017Sentinel-24350
9201552121052015Landsat-770073392017111111102017Sentinel-2989.6
10201552121052015Landsat-713,970340201861313062018Sentinel-2118,000
11201552121052015Landsat-732,250341201861313062018Sentinel-21,723,000
12201552121052015Landsat-710,980342201861313062018Sentinel-25774
13201552121052015Landsat-75936343201861313062018Sentinel-29821
142015666062015Landsat-742,470344201861616062018Sentinel-22742
152015666062015Landsat-73,587,000345201861616062018Sentinel-24328
16201561414062015Landsat-85524346201861818062018Sentinel-217,240
17201561414062015Landsat-8383,400347201861818062018Sentinel-211,090
18201561414062015Landsat-894,130348201861818062018Sentinel-213,350
19201561414062015Landsat-8128,600349201861818062018Sentinel-213,620
20201561414062015Landsat-85757350201861818062018Sentinel-225,740
21201561414062015Landsat-8129,000351201862121062018Sentinel-211,070
22201561414062015Landsat-87679352201862323062018Sentinel-285,820
23201561414062015Landsat-81,021,000353201862323062018Sentinel-220,690
24201561414062015Landsat-8758,00354201862323062018Sentinel-242,380
25201561414062015Landsat-8148,600355201862323062018Sentinel-217,320
26201563030062015Landsat-85,338,000356201862828062018Sentinel-2122,200
27201563030062015Landsat-8169,200357201862828062018Sentinel-26324
28201563030062015Landsat-815,2303582018711072018Sentinel-26329
29201563030062015Landsat-822,0403592018733072018Sentinel-2878,000
30201563030062015Landsat-896,9803602018733072018Sentinel-265,790
31201563030062015Landsat-81,430,0003612018766072018Sentinel-22,348,000
32201563030062015Landsat-8106,9003622018788072018Sentinel-2305,300
33201563030062015Landsat-831,6803632018788072018Sentinel-2534,800
34201563030062015Landsat-85,050,0003642018788072018Sentinel-215,860
35201563030062015Landsat-837,950,000365201871313072018Sentinel-2150,700
36201563030062015Landsat-8774,700366201871313072018Sentinel-237,080
37201571616072015Landsat-820,260,000367201871313072018Sentinel-280,570
38201571616072015Landsat-831,120368201871313072018Sentinel-27072
39201571616072015Landsat-8201,900369201871313072018Sentinel-218,240
40201571616072015Landsat-817,850370201871616072018Sentinel-214,810
41201571616072015Landsat-85852371201871818072018Sentinel-230,930
42201571616072015Landsat-85,139,000372201871818072018Sentinel-214,110
43201581717082015Landsat-85,744373201871818072018Sentinel-224,300
44201581717082015Landsat-8570,000374201871818072018Sentinel-232,030
45201581717082015Landsat-829,520375201872121072018Sentinel-22570
46201581717082015Landsat-8425,100376201872121072018Sentinel-218,520
47201581717082015Landsat-814,340377201872323072018Sentinel-212,020
48201581717082015Landsat-830,870378201872828072018Sentinel-222,140
49201581717082015Landsat-822,3403792018822082018Sentinel-218,320
50201581717082015Landsat-885,6803802018822082018Sentinel-220,130
51201581717082015Landsat-881,8503812018822082018Sentinel-21439
52201581717082015Landsat-810,690,0003822018822082018Sentinel-216,770
53201581717082015Landsat-8163,6003832018822082018Sentinel-230,940
54201581717082015Landsat-81,174,0003842018877082018Sentinel-2329,800
55201581717082015Landsat-832,3403852018877082018Sentinel-273,950
56201581717082015Landsat-8504,100386201881212082018Sentinel-231,160
57201581717082015Landsat-818,910387201881515082018Sentinel-251,780
58201581717082015Landsat-815,570388201881515082018Sentinel-2186,900
59201581717082015Landsat-86939389201881717082018Sentinel-270,270
60201581717082015Landsat-81,029,000390201881717082018Sentinel-21312
61201581717082015Landsat-840,340391201881717082018Sentinel-213,390
62201582121082015Sentinel-21,231,000392201882222082018Sentinel-243,870
63201582121082015Sentinel-257,210393201882222082018Sentinel-247,230
642015922092015Landsat-8473,800394201882525082018Sentinel-2710,500
652015922092015Landsat-866,060395201882727082018Sentinel-26,754,000
662015922092015Landsat-871,79039620189919092018Sentinel-259,990
672015922092015Landsat-8373,2003972018999092018Sentinel-2222,100
682015922092015Landsat-862303982018999092018Sentinel-215,700
692015922092015Landsat-85183399201892121092018Sentinel-227,820
702015922092015Landsat-88226400201892929092018Sentinel-2266,400
712015922092015Landsat-8789440120181066102018Sentinel-24,619,000
72201591818092015Landsat-8176,300402201942929042019Sentinel-218,330
73201591818092015Landsat-870,4604032019577052019Sentinel-223,560
742015101212102015Landsat-7303,800404201951717052019Sentinel-28492
752015101212102015Landsat-7329,500405201951919052019Sentinel-25064
76201661313062016Sentinel-2554,200406201952222052019Sentinel-24721
77201661313062016Sentinel-2204,500407201952222052019Sentinel-27206
78201661313062016Sentinel-265,400408201952424052019Sentinel-25133
79201661313062016Sentinel-292,880409201952424052019Sentinel-21989
80201661313062016Sentinel-294,020410201952424052019Sentinel-229,080
81201661313062016Sentinel-293,330411201952727052019Sentinel-29199
82201661313062016Sentinel-251,230412201952727052019Sentinel-22695
83201661313062016Sentinel-263,120413201952929052019Sentinel-237,140
84201661313062016Sentinel-219,170414201952929052019Sentinel-256,290
85201661313062016Sentinel-251,950415201952929052019Sentinel-25268
86201661313062016Sentinel-252,380416201952929052019Sentinel-26547
87201661313062016Sentinel-25987417201952929052019Sentinel-212,360
88201661313062016Sentinel-24187418201952929052019Sentinel-212,600
89201661313062016Sentinel-210,980419201952929052019Sentinel-29172
90201661313062016Sentinel-29773420201952929052019Sentinel-219,130
91201661313062016Sentinel-218,890421201952929052019Sentinel-223,230
92201661313062016Sentinel-214,860422201952929052019Sentinel-210,570
93201661313062016Sentinel-213,040423201952929052019Sentinel-212,020
94201661313062016Sentinel-235,750424201952929052019Sentinel-29245
95201661313062016Sentinel-222,730425201952929052019Sentinel-25655
96201661313062016Sentinel-210,2304262019611062019Sentinel-25789
97201661313062016Sentinel-220734272019611062019Sentinel-29754
98201661313062016Sentinel-264314282019611062019Sentinel-29351
99201661313062016Sentinel-262604292019611062019Sentinel-23766
100201661313062016Sentinel-293794302019611062019Sentinel-25845
101201661313062016Sentinel-255584312019611062019Sentinel-25484
102201661313062016Sentinel-220,5804322019611062019Sentinel-22018
103201661313062016Sentinel-261654332019633062019Sentinel-2283,100
104201661313062016Sentinel-222,4104342019633062019Sentinel-25181
105201661313062016Sentinel-217,1704352019633062019Sentinel-278,140
106201661313062016Sentinel-228724362019633062019Sentinel-221,100
107201661313062016Sentinel-217,0404372019633062019Sentinel-28376
108201661313062016Sentinel-234334382019666062019Sentinel-26335
109201661313062016Sentinel-228334392019688062019Sentinel-216,510
110201662626062016Sentinel-27,536,0004402019688062019Sentinel-212,400
111201662626062016Sentinel-21,304,0004412019688062019Sentinel-237,900
112201662626062016Sentinel-2280,8004422019688062019Sentinel-226,980
113201662626062016Sentinel-2488,3004432019688062019Sentinel-211,840
114201662626062016Sentinel-263,250444201961111062019Sentinel-252,760
115201662626062016Sentinel-236,640445201961111072019Sentinel-2157,400
116201662626062016Sentinel-291,530446201961111072019Sentinel-2236,700
117201662626062016Sentinel-266,310447201961313062019Sentinel-2303,000
118201662626062016Sentinel-25005448201961313062019Sentinel-258,870
119201662626062016Sentinel-225,710449201961313062019Sentinel-215,500
120201662626062016Sentinel-2252,400450201961313062019Sentinel-247,520
121201662626062016Sentinel-229,010451201961313062019Sentinel-23789
1222016733072016Sentinel-2168,200452201961313062019Sentinel-24943
1232016733072016Sentinel-250,910453201961313062019Sentinel-217,630
1242016733072016Sentinel-222,700454201961313062019Sentinel-21962
1252016733072016Sentinel-236,420455201961313062019Sentinel-24811
1262016733072016Sentinel-221,480456201961818062019Sentinel-220,380
1272016733072016Sentinel-253,180457201961818062019Sentinel-245,740
1282016733072016Sentinel-248,910458201961818062019Sentinel-215,340
1292016733072016Sentinel-245,660459201961818062019Sentinel-29706
1302016733072016Sentinel-221,640460201961818062019Sentinel-27033
1312016733072016Sentinel-2328,700461201961818062019Sentinel-27621
1322016733072016Sentinel-2315,300462201961818062019Sentinel-22,261,000
1332016733072016Sentinel-2114,600463201961818062019Sentinel-2188,600
1342016733072016Sentinel-28185464201962323062019Sentinel-2169,300
1352016733072016Sentinel-216,380465201962323062019Sentinel-2193,200
1362016733072016Sentinel-213,850466201962323062019Sentinel-2827,800
1372016733072016Sentinel-217,550467201962323062019Sentinel-233,610
1382016733072016Sentinel-230,480468201962323062019Sentinel-29397
1392016733072016Sentinel-220,500469201962828062019Sentinel-220,640
1402016733072016Sentinel-27337470201962828062019Sentinel-211,670
1412016733072016Sentinel-212,400471201962828062019Sentinel-2400,700
1422016733072016Sentinel-27921472201962828062019Sentinel-231,680
1432016733072016Sentinel-24858473201962828062019Sentinel-22599
1442016733072016Sentinel-213,790474201962828062019Sentinel-239,200
1452016733072016Sentinel-256914752019711072019Sentinel-25961
1462016733072016Sentinel-222344762019711072019Sentinel-25371
1472016733072016Sentinel-295174772019733072019Sentinel-268,200
1482016733072016Sentinel-227,4204782019733072019Sentinel-212,410
149201671616072016Sentinel-210,580,0004792019733072019Sentinel-2127,800
150201671616072016Sentinel-2294,7004802019733072019Sentinel-283,930
151201671616072016Sentinel-21,651,0004812019733072019Sentinel-233,420
152201671616072016Sentinel-282154822019733072019Sentinel-24297
153201671616072016Sentinel-22,646,0004832019733072019Sentinel-229,800
154201671616072016Sentinel-260524842019733072019Sentinel-210,410
155201671616072016Sentinel-236,4404852019766072019Sentinel-23221
156201671616072016Sentinel-217,3404862019788072019Sentinel-231,720
157201671616072016Sentinel-275914872019788072019Sentinel-27330
158201671616072016Sentinel-2113,5004882019788072019Sentinel-2605,900
159201671616072016Sentinel-210,4804892019788072019Sentinel-25203
160201671616072016Sentinel-257884902019788072019Sentinel-215,210
161201671616072016Sentinel-261074912019788072019Sentinel-217,080
162201671616072016Sentinel-22623492201971313072019Sentinel-211,430
163201671616072016Sentinel-211,340493201971313072019Sentinel-226,720
164201671616072016Sentinel-224,420494201971414052019Sentinel-210,360
165201671616072016Sentinel-21178495201971616072019Sentinel-284,870
166201671616072016Sentinel-214,990496201971616072019Sentinel-256,510
167201671616072016Sentinel-215,330497201971818072019Sentinel-26161
168201671616072016Sentinel-2348,700498201972323072019Sentinel-213,120
169201671616072016Sentinel-216,630499201972323072019Sentinel-215,820
170201671616072016Sentinel-26034500201972323072019Sentinel-232,250
171201671616072016Sentinel-23501501201972828072019Sentinel-265,860
172201671616072016Sentinel-27033502201972828072019Sentinel-260,370
173201671616072016Sentinel-26112503201973131072019Sentinel-210,880,000
174201671616072016Sentinel-21992504201973131072019Sentinel-2150,400
175201671616072016Sentinel-2854.9505201973131072019Sentinel-260,240
176201671616072016Sentinel-2243,100506201973131072019Sentinel-22749
177201672323072016Sentinel-279,7105072019855082019Sentinel-29158
178201672323072016Sentinel-256225082019877082019Sentinel-2530,800
179201672323072016Sentinel-240535092019877082019Sentinel-2266,100
180201672323072016Sentinel-215,1205102019877082019Sentinel-213,530,000
181201672323072016Sentinel-253745112019877082019Sentinel-2105,500
182201672323072016Sentinel-2386.8512201981717082019Sentinel-21,737,000
183201672323072016Sentinel-2479.5513201981717082019Sentinel-295,000
184201672323072016Sentinel-23878514201981717082019Sentinel-2289,000
185201672323072016Sentinel-23388515201981717082019Sentinel-234,110
186201672323072016Sentinel-211,770516201982020082019Sentinel-257,440
187201672323072016Sentinel-21466517201982222082019Sentinel-261,490
188201672323072016Sentinel-21515518201982222082019Sentinel-23327
189201672323072016Sentinel-23239519201983030082019Sentinel-28,957,000
190201672323072016Sentinel-2388,900520201983030082019Sentinel-223,200
191201672323072016Sentinel-299265212019911092019Sentinel-22,280,000
192201672323072016Sentinel-233905222019911092019Sentinel-21,438,000
193201672323072016Sentinel-2219,5005232019911092019Sentinel-29084
194201672323072016Sentinel-289155242019911092019Sentinel-2150,700
195201672323072016Sentinel-252,4105252019911092019Sentinel-22,810,000
196201672323072016Sentinel-239,2105262019911092019Sentinel-244,940
197201672323072016Sentinel-219,7105272019911092019Sentinel-226,510
198201672323072016Sentinel-225,6205282019911092019Sentinel-22835
199201672323072016Sentinel-218,4505292019911092019Sentinel-2106,400
200201672323072016Sentinel-267,4605302019911092019Sentinel-27332
201201672323072016Sentinel-248435312019966092019Sentinel-212,650
202201672323072016Sentinel-297,4205322019966092019Sentinel-2179,500
203201672323072016Sentinel-223,4505332019966092019Sentinel-212,040
204201672323072016Sentinel-220,2305342019966092019Sentinel-267,160
205201672323072016Sentinel-25118535201991111092019Sentinel-21,058,000
206201672323072016Sentinel-22605536201992121092019Sentinel-264,710
207201672626072016Sentinel-288,380537201992424092019Sentinel-227,290
208201672626072016Sentinel-222,470538201992424092019Sentinel-244,440
209201672626072016Sentinel-2672053920191044102019Sentinel-213,200
210201672626072016Sentinel-211,98054020191066102019Sentinel-262,150
211201672626072016Sentinel-2441754120191066102019Sentinel-270,370
212201672626072016Sentinel-2743754220191066102019Sentinel-248,690
213201672626072016Sentinel-211,37054320191066102019Sentinel-224,760
21420168266072016Sentinel-2887654420191066102019Sentinel-210,610
2152016855082016Sentinel-213,67054520191066102019Sentinel-218,430
2162016855082016Sentinel-231,27054620191066102019Sentinel-293,260
2172016855082016Sentinel-210,0005472019101111102019Sentinel-251,390
2182016855082016Sentinel-218,9305482019101111102019Sentinel-269,020
2192016855082016Sentinel-226,1905492019101616102019Sentinel-2317,100
220201681212082016Sentinel-264055502019101616102019Sentinel-222,580
221201681212082016Sentinel-24811551202042121042020Sentinel-221,670
222201681212082016Sentinel-222,9305522020511052020Sentinel-2121,700
223201681212082016Sentinel-224,8305532020511052020Sentinel-219,810
224201681212082016Sentinel-263,220554202051313052020Sentinel-219,930
225201681212082016Sentinel-210,610555202052323052020Sentinel-2887,700
226201681212082016Sentinel-224,650556202052323052020Sentinel-223,840
227201681212082016Sentinel-215,390557202053131052020Sentinel-21,312,000
228201681515082016Sentinel-275,710558202053131052020Sentinel-22,885,000
229201681515082016Sentinel-218,650559202053131052020Sentinel-291,550
230201682525082016Sentinel-275,080560202053131052020Sentinel-215,290
231201682525082016Sentinel-230,6905612020622062020Sentinel-226,690
232201682525082016Sentinel-218,4705622020622062020Sentinel-2115,600
233201682525082016Sentinel-225,3205632020622062020Sentinel-261,700
234201682525082016Sentinel-243545642020622062020Sentinel-226,580
235201682525082016Sentinel-249005652020622062020Sentinel-28253
236201682525082016Sentinel-242,3805662020622062020Sentinel-2101,300
2372016911092016Sentinel-22,469,0005672020655062020Sentinel-2161,800
2382016911092016Sentinel-212,790,0005682020655062020Sentinel-210,270
2392016911092016Sentinel-232,7405692020655062020Sentinel-246,650
2402016911092016Sentinel-2227,1005702020655062020Sentinel-278,100
2412016911092016Sentinel-240,9105712020655062020Sentinel-213,910
2422016911092016Sentinel-222,5905722020677062020Sentinel-212,770
2432016911092016Sentinel-211,4405732020677062020Sentinel-27744
2442017599052017Sentinel-218,1405742020677062020Sentinel-22224
2452017599052017Sentinel-217645752020677062020Sentinel-239,130
2462017599052017Sentinel-266715762020677062020Sentinel-2101,400
2472017599052017Sentinel-223,6405772020677062020Sentinel-22534
2482017599052017Sentinel-290065782020677062020Sentinel-238,610
2492017599052017Sentinel-278515792020677062020Sentinel-213,470
250201751212052017Sentinel-227,3705802020677062020Sentinel-238,620
251201752222052017Sentinel-260,6705812020677062020Sentinel-23642
252201752222052017Sentinel-250,0505822020677062020Sentinel-25200
253201752222052017Sentinel-283,4905832020677062020Sentinel-25954
254201752222052017Sentinel-218,550584202061010062020Sentinel-2132,800
255201752222052017Sentinel-216,510585202061212062020Sentinel-235,480
256201752222052017Sentinel-21706586202061212062020Sentinel-2216,000
257201752222052017Sentinel-232,730587202061212062020Sentinel-236,830
258201752222052017Sentinel-2205,600588202061212062020Sentinel-213,670
259201752222052017Sentinel-25801589202061212062020Sentinel-2116,800
260201752222052017Sentinel-234,290590202061212062020Sentinel-220,480
261201752222052017Sentinel-226,930591202061717062020Sentinel-2262,300
262201752929052017Sentinel-212,960592202061717062020Sentinel-227,250
263201752929052017Sentinel-251,440593202061717062020Sentinel-258,620
264201752929052017Sentinel-242,960594202062222062020Sentinel-238,590
265201752929052017Sentinel-273,730595202062222062020Sentinel-212,970
266201752929052017Sentinel-24667596202062525062020Sentinel-2642,900
267201752929052017Sentinel-22779597202062525062020Sentinel-27354
268201752929052017Sentinel-25866598202062525062020Sentinel-2569,400
269201752929052017Sentinel-220,300599202062727062020Sentinel-2381,700
270201752929052017Sentinel-25799600202062727062020Sentinel-2122,400
271201752929052017Sentinel-216,910601202062727062020Sentinel-29,103,000
272201752929052017Sentinel-25432602202062727062020Sentinel-27893
273201752929052017Sentinel-24607603202062727062020Sentinel-24949
274201752929052017Sentinel-249,690604202062727062020Sentinel-211,370
2752017611062017Sentinel-21,105,000605202062727062020Sentinel-23730
276201761111062017Sentinel-276,350606202063030062020Sentinel-217,380
277201761111062017Sentinel-212,840607202063030062020Sentinel-28559
2782017711072017Sentinel-2169,0006082020722072020Sentinel-21,105,000
2792017711072017Sentinel-2293,006092020722072020Sentinel-26,691,000
2802017711072017Sentinel-2440,6006102020722072020Sentinel-233,270
2812017711072017Sentinel-22,732,0006112020722072020Sentinel-21,209,000
2822017711072017Sentinel-211,5506122020722072020Sentinel-255,370
2832017711072017Sentinel-285726132020722072020Sentinel-249,600
2842017711072017Sentinel-268756142020722072020Sentinel-224,310
2852017711072017Sentinel-217,3006152020722072020Sentinel-289,960
2862017733072017Sentinel-254,7906162020722072020Sentinel-219,590
2872017733072017Sentinel-216,8906172020722072020Sentinel-26726
2882017766072017Sentinel-2497,6006182020755072020Sentinel-261,330
289201771111072017Sentinel-2134,1006192020755072020Sentinel-2328,500
290201771111072017Sentinel-2918,300620202073131072018Sentinel-27608
291201771111072017Sentinel-218,5706212020777072020Sentinel-24874
292201771111072017Sentinel-265396222020777072020Sentinel-218,690
293201771111072017Sentinel-214,8306232020777072020Sentinel-219,190
294201772323072017Sentinel-2172,6006242020777072020Sentinel-23615
295201772323072017Sentinel-231,6606252020777072020Sentinel-2108,600
296201772323072017Sentinel-2169,0006262020777072020Sentinel-210,350
297201772323072017Sentinel-22053627202071212072020Sentinel-220,180
298201772323072017Sentinel-210,480628202071515072020Sentinel-2230,700
299201772323072017Sentinel-210,980,000629202071515072020Sentinel-213,790
300201772323072017Sentinel-25,090,000630202071515072020Sentinel-2897,900
301201772323072017Sentinel-2180,800631202071515072020Sentinel-21,517,000
302201772828072017Sentinel-227,240632202071717072020Sentinel-21,579,000
303201772828072017Sentinel-2400,900633202071717072020Sentinel-2929,500
304201772828072017Sentinel-2187,700634202071717072020Sentinel-2126,200
305201781212082017Sentinel-2213,800635202071717072020Sentinel-22,455,000
306201781212082017Sentinel-2451,400636202072222072020Sentinel-21,013,000
307201781212082017Sentinel-217,280637202072222072020Sentinel-2453,700
308201781212082017Sentinel-24320638202072222072020Sentinel-22,366,000
309201781212082017Sentinel-250,260639202072727072020Sentinel-273,270
310201781212082017Sentinel-237,790640202072727072020Sentinel-21,127,000
311201781212082017Sentinel-2187,400641202073030072020Sentinel-2416,600
312201781212082017Sentinel-2786,3006422020811082020Sentinel-24773
313201781212082017Sentinel-230,9506432020866082020Sentinel-23468
314201781212082017Sentinel-249526442020866082020Sentinel-259,090
315201782020082017Sentinel-274,950645202081111082020Sentinel-23476
316201782020082017Sentinel-223,410646202081111082020Sentinel-2178,900
317201782020082017Sentinel-210,850647202081111082020Sentinel-2568,300
3182017911092017Sentinel-2688,400648202081111082020Sentinel-29427
3192017911092017Sentinel-2110,500649202081111082020Sentinel-217,940
3202017911092017Sentinel-223,360650202081111082020Sentinel-291,750
3212017966092017Sentinel-22,563,000651202081111082020Sentinel-270,120
322201792424092017Sentinel-225,980652202081414082020Sentinel-222,120
323201792424092017Sentinel-211,490653202081616082020Sentinel-2342,500
324201792626072017Sentinel-226,630654202081616082020Sentinel-2196,800
325201792626092017Sentinel-226,070655202081919082020Sentinel-245,260
326201792626092017Sentinel-257,440656202082121082020Sentinel-259,790
327201792626092017Sentinel-2151,800657202082626082020Sentinel-2533,100
328201792626092017Sentinel-247,720658202083131082020Sentinel-295,980
329201792626092017Sentinel-21555659202091515092020Sentinel-22,568,000
330201792626092017Sentinel-220,590660202091515092020Sentinel-276,720

Appendix C. The Range of Estimation for Each Factor by Logistic Regression Results of Training Dataset for Each Group

FactorGroup1Group2Group3Group4Group5Group6Group7
Intercept−86.1−118.0−125.0−160.000−162.0−171.0−160
Slope0.0310.0340.0250.0250.0260.0270.0264
Aspect -N0.1060.0940.1710.1490.1430.1490.149
Aspect -NE0.0800.0920.1250.1190.1020.1170.119
Aspect -E−0.106−0.093−0.050−0.050−0.065−0.053−0.0502
Aspect -SE−0.146−0.138−0.127−0.118−0.115−0.116−0.116
Aspect -S−0.017−0.009−0.084−0.065−0.056−0.061−0.0658
Aspect -SW0.0260.029−0.066−0.054−0.050−0.052−0.0557
Aspect -W0.1240.1100.0690.0660.0760.0640.0659
Aspect -NW−0.073−0.089−0.042−0.051−0.038−0.051−0.0487
Elevation0.0020.0020.0020.0020.0020.0020.00193
Precipitation0.0060.0100.0260.0270.0160.0210.0268
Distance to Road0.0000.0000.0000.0000.0000.000−0.00036
Distance to Settlements0.0010.0010.0010.0010.0010.0010.000581
LULC1-Open Shrub57.30073.400115.000122.00095.100110.000122
LULC2-Grassland58.10074.100116.000122.00095.800111.000122
LULC4-Cropland58.20074.300116.000122.00096.000111.000122
LULC5-Urban Area57.50073.500116.000122.00095.600110.000122
LULC6-Barren47.20063.200105.000111.00084.70099.400111
Annual Air Temperature−0.3160.000−1.700−0.9220.156−0.166−0.909
Wind Speed9.50010.80019.40018.90013.80016.10018.8
NDVI−2.110−1.800−2.440−2.200−2.360−2.230−2.22
Distance to streams−0.001−0.0020.000−0.002−0.002−0.002−0.00182
Geology-Flood Plain//1.8501.8101.7001.8001.81
Geology-Slope//−0.295−0.313−0.290−0.326−0.315
Geology-Polygenetic//−0.595−0.606−0.586−0.596−0.607
Geology-Mukdadiyah//−0.164−0.232−0.295−0.257−0.231
Geology-Injana//0.8060.7470.7160.7420.748
Geology-Fatha//0.9610.9911.0100.9850.992
Geology-Pilaspi//0.5710.5950.5660.5820.594
Geology-Gercus//0.0420.0700.0900.0840.069
Geology-Khurmala and Sinjar//−0.138−0.109−0.052−0.074−0.11
Geology-Kolosh//−0.434−0.365−0.314−0.339−0.363
Geology-Tanjero, Aqra, and Bekhme//−2.610−2.590−2.550−2.610−2.59
TWI/0.0100.012/0.0120.0120.012
Relative Humidity/0.102−0.890−0.5170.248/−0.51
Distance to Farmland/0.0000.00020.004/0.000130.000149
TPI/0.004/0.00430.0040.0040.00426

Appendix D. The Range of Estimation for Each Factor by Logistic Regression Results of Validation Dataset for Each Group

FactorGroup1Group2Group3Group4Group5Group6Group7
Intercept−305−22.31104.511.36−142−371.28.71
Slope0.028250.033290.028620.028370.029040.035510.03116
Aspect-N0.31240.24890.24230.21050.2220.23040.2112
Aspect-NE0.062580.001159−0.05362−0.05304−0.051120.01329−0.05428
Aspect-E−0.4681−0.4903−0.5028−0.5033−0.5134−0.481−0.506
Aspect-SE−0.1643−0.1617−0.1732−0.1629−0.1421−0.1555−0.1593
Aspect-S0.15070.22010.22190.25960.24130.23040.2552
Aspect-SW0.10.18350.1870.21970.19280.18020.2158
Aspect-W0.01248−0.007560.030130.011710.02498−0.001350.01194
Aspect-NW−0.011850.000390.043040.012280.02091−0.024140.01988
Elevation0.000230.0002240.0017380.0005620.000136−0.000990.000541
Precipitation0.13670.23640.26150.25390.19560.16150.2537
Distance to Road−0.00028−0.00048−0.0005856−0.00062−0.0006−0.00044−0.00062
Distance to Settlements0.000117−4.3 × 10−54.081 × 10−54.77 × 10−50.0001879.28 × 10−54.66 × 10−5
LULC1-Open Shrub10.659.8399.4519.4549.96610.419.5
LULC2-Grassland11.5810.7310.3310.3610.811.3810.41
LULC4-Cropland11.3410.469.9669.99510.3911.1210.04
LULC5-Urban Area12.0311.9311.7111.7311.5612.3311.78
LULC6-Barren1.4051.331.4351.3531.2371.8371.387
Annual Air Temperature−0.3361−11.91−16.66−13.64−7.4920.6016−13.58
Wind Speed73.28137.2157.9149.9117.184.91149.9
NDVI−3.157−2.026−2.708−2.052−2.248−2.108−2.069
Distance to streams0.001483−0.001430.000881−0.00194−0.00169−0.0019−0.00157
Geology-Flood Plain//−0.2976−0.4112−0.7447−0.4617−0.4108
Geology-Slope//−0.8035−0.9317−1.242−1.355−0.9299
Geology-Polygenetic//1.5211.5132.0381.6381.501
Geology-Mukdadiyah//1.6871.4930.90221.1171.498
Geology-Injana//−0.3107−0.4829−0.6287−0.4022−0.479
Geology-Fatha//0.4430.47520.47890.47020.4813
Geology-Pilaspi//−0.7157−0.6053−0.6636−0.4154−0.6063
Geology-Gercus//−0.9819−0.8063−0.8455−0.3796−0.8063
Geology-Khurmala and Sinjar//2.4992.6123.0262.5422.583
Geology-Kolosh//−0.7934−0.5782−0.7149−0.2305−0.5707
Geology-Tanjero, Aqra, and Bekhme//−2.256−2.287−1.617−2.532−2.27
TWI/0.032290.03427/0.038740.043880.03755
Relative Humidity/−7.941−10.84−9.169−4.852/−9.132
Distance to Farmland/0.0004450.00045330.000418/0.000210.000418
TPI/0.006574/0.00720.0086140.010520.007242

References

  1. Ziadé, R.; Abdallah, C.; Baghdadi, N. The effect of forest fire on mass movement in Lebanese mountainous areas. Int. J. Wildl. Fire 2014, 23, 845–859. [Google Scholar] [CrossRef]
  2. Salar, S.G. Identification of Forest Fire Risk Zones in the Sartak Watershed, Kurdistan Region, Iraq: Gis and Remote Sensing Applications. Iraqi Bull. Geol. Min. 2019, 15, 105–122. [Google Scholar]
  3. Pourghasemi, H.R.; Gayen, A.; Panahi, M.; Rezaie, F.; Blaschke, T. Multi-hazard probability assessment and mapping in Iran. Sci. Total Environ. 2019, 692, 556–571. [Google Scholar] [CrossRef] [PubMed]
  4. Pourtaghi, Z.S.; Pourghasemi, H.R.; Rossi, M. Forest fire susceptibility mapping in the Minudasht forests, Golestan province, Iran. Environ. Earth Sci. 2015, 73, 1515–1533. [Google Scholar]
  5. Kayet, N.; Chakrabarty, A.; Pathak, K.; Sahoo, S.; Dutta, T.; Hatai, B.K. Comparative analysis of multi-criteria probabilistic FR and AHP models for forest fire risk (FFR) mapping in Melghat Tiger Reserve (MTR) forest. J. For. Res. 2020, 31, 565–579. [Google Scholar] [CrossRef]
  6. Zanini, G.D.; Sbrissia, A.F. Prescribed burning of pastures: A management tool? Rev. Ciências Agroveterinárias 2014, 12, 94–103. [Google Scholar]
  7. Jahdi, R.; Salis, M.; Darvishsefat, A.A.; Mostafavi, M.A.; Alcasena, F.; Etemad, V.; Lozano, O.; Spano, D. Calibration of FARSITE simulator in northern Iranian forests. Nat. Hazards Earth Syst. Sci. 2015, 15, 443–459. [Google Scholar] [CrossRef] [Green Version]
  8. Glasa, J.; Halada, L. On elliptical model for forest fire spread modeling and simulation. Math. Comput. Simul. 2008, 78, 76–88. [Google Scholar] [CrossRef]
  9. Ashiagbor, G.; Oduro, W.; Gyiele, L.; Siaw, D.; Barnes, V.R.; Agbenyega, O.; Twum-Ampofo, K.; Partey, S.; Thevathasan, N.; Gordon, A.; et al. Toward sustainable land resources management with agroforestry: Empirical evidence from the Sunyani west district of Ghana. Agrofor. Syst. 2020, 94, 527–537. [Google Scholar] [CrossRef]
  10. Esmaeili Sharif, M.; Jalilvand, H.; Amoozad, M.; Jafari, A.; Sayed Mahale, S.M. The effect of ecological factors on fire in Hyrcanian forests (Case study: Forest areas of Neka, Mazandaran, Iran). J. For. Res. Dev. 2018, 4, 113–129. [Google Scholar]
  11. Sowmya, S.V.; Somashekar, R.K. Application of remote sensing and geographical information system in mapping forest fire risk zone at Bhadra wildlife sanctuary, India. J. Environ. Biol. 2010, 31, 969. [Google Scholar] [PubMed]
  12. Adab, H.; Kanniah, K.D.; Solaimani, K. Modeling forest fire risk in the northeast of Iran using remote sensing and GIS techniques. Nat. Hazards 2013, 65, 1723–1743. [Google Scholar] [CrossRef]
  13. Jung, J.; Kim, C.; Jayakumar, S.; Kim, S.; Han, S.; Kim, D.H.; Heo, J. Forest fire risk mapping of Kolli Hills, India, considering subjectivity and inconsistency issues. Nat. Hazards 2013, 65, 2129–2146. [Google Scholar] [CrossRef]
  14. Amalina, P.; Prasetyo, L.B.; Rushayati, S.B. Forest Fire Vulnerability Mapping in Way Kambas National Park. Procedia Environ. Sci. 2016, 33, 239–252. [Google Scholar] [CrossRef] [Green Version]
  15. Gigović, L.; Pourghasemi, H.R.; Drobnjak, S.; Bai, S. Testing a new ensemble model based on SVM and random forest in forest fire susceptibility assessment and its mapping in Serbia’s Tara National Park. Forests 2019, 10, 408. [Google Scholar] [CrossRef] [Green Version]
  16. Rasooli, S.B.; Rashidi, S.; Shabanian, N. Fire sensitivity of broadleaf tree species in plantations of Kurdistan, Iran. J. For. Res. 2021, 32, 1167–1176. [Google Scholar] [CrossRef]
  17. Vadrevu, K.P.; Eaturu, A.; Badarinath, K.V.S. Fire risk evaluation using multicriteria analysis-a case study. Environ. Monit. Assess. 2010, 166, 223–239. [Google Scholar] [CrossRef]
  18. Santiago, I.T.; Kheladze, N. GIS Wildland fire Hazard Modeling in Georgia. Matra Project Report, Caucasus Environmental NGO Network. 2011. Available online: http://drm.cenn.org/Local_Case_studies/GIS%20Wildland%20Fire%20Hazard%20Modeling%20in%20Georgia1.pdf (accessed on 17 July 2021).
  19. Sachdeva, S.; Bhatia, T.; Verma, A.K. GIS-based evolutionary optimized Gradient Boosted Decision Trees for forest fire susceptibility mapping. Nat. Hazards 2018, 92, 1399–1418. [Google Scholar] [CrossRef]
  20. Zhang, G.; Wang, M.; Liu, K. Forest Fire Susceptibility Modeling Using a Convolutional Neural Network for Yunnan Province of China. Int. J. Disaster Risk Sci. 2019, 10, 386–403. [Google Scholar] [CrossRef] [Green Version]
  21. Moayedi, H.; Mehrabi, M.; Bui, D.T.; Pradhan, B.; Foong, L.K. Fuzzy-metaheuristic ensembles for spatial assessment of forest fire susceptibility. J. Environ. Manag. 2020, 260, 109867. [Google Scholar] [CrossRef]
  22. Hong, H.; Naghibi, S.A.; Dashtpagerdi, M.M.; Pourghasemi, H.R.; Chen, W. A comparative assessment between linear and quadratic discriminant analyses (LDA-QDA) with frequency ratio and weights-of-evidence models for forest fire susceptibility mapping in China. Arab. J. Geosci. 2017, 10, 167. [Google Scholar] [CrossRef]
  23. Bui, D.T.; Hoang, N.-D.; Samui, P. Spatial pattern analysis and prediction of forest fire using new machine learning approach of Multivariate Adaptive Regression Splines and Differential Flower Pollination optimization: A case study at Lao Cai province (Viet Nam). J. Environ. Manag. 2019, 237, 476–487. [Google Scholar]
  24. Stolzenburg, W. Fire in the Rain Forest; National Emergency Training Center: Emmitsburg, MD, USA, 2001.
  25. Food and Agriculture Organization. Fire Management: Global Assessment 2006. FAO For. Pap. 2007, 151, 135. [Google Scholar]
  26. Pourghasemi, H.R.; Gayen, A.; Lasaponara, R.; Tiefenbacher, J.P. Application of learning vector quantization and different machine learning techniques to assessing forest fire influence factors and spatial modelling. Environ. Res. 2020, 184, 109321. [Google Scholar] [CrossRef] [PubMed]
  27. Razavi-Termeh, S.V.; Sadeghi-Niaraki, A.; Choi, S.-M. Ubiquitous GIS-based forest fire susceptibility mapping using artificial intelligence methods. Remote Sens. 2020, 12, 1689. [Google Scholar] [CrossRef]
  28. Stolle, F.; Chomitz, K.M.; Lambin, E.F.; Tomich, T.P. Human ecological intervention and the role of forest fires in human ecology. For. Ecol. Manag. 2003, 179, 277–292. [Google Scholar] [CrossRef] [Green Version]
  29. de Alvino-Rayol, F.O.; Rayol, B.P. Effect of fire in weeds of agroforestry system, Pará, Brazil. Rev. Ciências Agroveterinárias 2020, 19, 1–7. [Google Scholar] [CrossRef] [Green Version]
  30. de Alvino-Rayol, F.O.; Rayol, B.P. Effect of fire in seed banks of agroforestry systems, Pará, Brazil. Rev. Ciências Agroveterinárias 2019, 18, 489–498. [Google Scholar] [CrossRef] [Green Version]
  31. Minnich, R.A. An integrated model of two fire regimes. Conserv. Biol. 2001, 15, 1549–1553. [Google Scholar] [CrossRef]
  32. Broncano, M.J.; Retana, J. Topography and pre-fire vegetation affecting the spatial heterogeneity generated after a large forest wildfire in the Mediterranean Basin. Int. J. Wildl. Fire 2004, 13, 209–216. [Google Scholar] [CrossRef]
  33. Dlamini, W.M. A Bayesian belief network analysis of factors influencing wildfire occurrence in Swaziland. Environ. Model. Softw. 2010, 25, 199–208. [Google Scholar] [CrossRef]
  34. Pourtaghi, Z.S.; Pourghasemi, H.R.; Aretano, R.; Semeraro, T. Investigation of general indicators influencing on forest fire and its susceptibility modeling using different data mining techniques. Ecol. Indic. 2016, 64, 72–84. [Google Scholar] [CrossRef]
  35. Ghorbanzadeh, O.; Blaschke, T.; Gholamnia, K.; Aryal, J. Forest fire susceptibility and risk mapping using social/infrastructural vulnerability and environmental variables. Fire 2019, 2, 50. [Google Scholar] [CrossRef] [Green Version]
  36. Vázquez, A.; Moreno, J.M. Spatial distribution of forest fires in Sierra de Gredos (Central Spain). For. Ecol. Manag. 2001, 147, 55–65. [Google Scholar] [CrossRef]
  37. Lozano, F.J.; Suárez-Seoane, S.; Kelly, M.; Luis, E. A multi-scale approach for modeling fire occurrence probability using satellite data and classification trees: A case study in a mountainous Mediterranean region. Remote Sens. Environ. 2008, 112, 708–719. [Google Scholar] [CrossRef]
  38. Zumbrunnen, T.; Pezzatti, G.B.; Menéndez, P.; Bugmann, H.; Bürgi, M.; Conedera, M. Weather and human impacts on forest fires: 100 years of fire history in two climatic regions of Switzerland. For. Ecol. Manag. 2011, 261, 2188–2199. [Google Scholar] [CrossRef]
  39. Justus, J. A case study in concept determination: Ecological diversity. In Philosophy of Ecology; Elsevier: Amsterdam, The Netherlands, 2011; pp. 147–168. [Google Scholar]
  40. Jaiswal, R.K.; Mukherjee, S.; Raju, K.D.; Saxena, R. Forest fire risk zone mapping from satellite imagery and GIS. Int. J. Appl. Earth Obs. Geoinf. 2002, 4, 1–10. [Google Scholar] [CrossRef]
  41. Dong, X.U.; Li-min, D.A.I.; Guo-fan, S.; Lei, T.; Hui, W. Forest fire risk zone mapping from satellite images and GIS for Baihe Forestry Bureau, Jilin, China. J. For. Res. 2005, 16, 169–174. [Google Scholar] [CrossRef]
  42. Keane, R.E.; Drury, S.A.; Karau, E.C.; Hessburg, P.F.; Reynolds, K.M. A method for mapping fire hazard and risk across multiple scales and its application in fire management. Ecol. Modell. 2010, 221, 2–18. [Google Scholar] [CrossRef]
  43. Paz, S.; Carmel, Y.; Jahshan, F.; Shoshany, M. Post-fire analysis of pre-fire mapping of fire-risk: A recent case study from Mt. Carmel (Israel). For. Ecol. Manag. 2011, 262, 1184–1188. [Google Scholar] [CrossRef]
  44. Eskandari, S.; Ghadikolaei, J.O.; Jalilvand, H.; Saradjian, M.R. Detection of fire high-risk areas in northern forests of Iran using Dong model. World Appl. Sci. J. 2013, 27, 770–773. [Google Scholar]
  45. Rasooli, S.B.; Bonyad, A.E. Evaluating the efficiency of the Dong model in determining fire vulnerability in Iran’s Zagros forests. J. For. Res. 2019, 30, 1447–1458. [Google Scholar] [CrossRef]
  46. Mahdavi, A. Forests and rangelands? wildfire risk zoning using GIS and AHP techniques. Casp. J. Environ. Sci. 2012, 10, 43–52. [Google Scholar]
  47. Rasooli, S.B.; Bonyad, A.E.; Pir Bavaghar, M. Forest fire vulnerability map using remote sensing data, GIS and AHP analysis (Case study: Zarivar Lake surrounding area). Casp. J. Environ. Sci. 2018, 16, 369–377. [Google Scholar]
  48. Zarekar, A.; Zamani, B.K.; Ghorbani, S.; Moalla, M.A.; Jafari, H. Mapping spatial distribution of forest fire using MCDM and GIS (case study: Three forest zones in Guilan Province). Iran. J. For. Poplar Res. 2013, 21, 218–230. [Google Scholar]
  49. Bui, D.T.; Bui, Q.-T.; Nguyen, Q.-P.; Pradhan, B.; Nampak, H.; Trinh, P.T. A hybrid artificial intelligence approach using GIS-based neural-fuzzy inference system and particle swarm optimization for forest fire susceptibility modeling at a tropical area. Agric. For. Meteorol. 2017, 233, 32–44. [Google Scholar]
  50. Gheshlaghi, H.A. Using GIS to develop a model for forest fire risk mapping. J. Indian Soc. Remote Sens. 2019, 47, 1173–1185. [Google Scholar] [CrossRef]
  51. Martínez, J.; Vega-Garcia, C.; Chuvieco, E. Human-caused wildfire risk rating for prevention planning in Spain. J. Environ. Manag. 2009, 90, 1241–1252. [Google Scholar] [CrossRef]
  52. Jurdao, S.; Chuvieco, E.; Arevalillo, J.M. Modelling fire ignition probability from satellite estimates of live fuel moisture content. Fire Ecol. 2012, 8, 77–97. [Google Scholar] [CrossRef]
  53. Mohammadi, F.; Bavaghar, M.P.; Shabanian, N. Forest Fire Risk Zone Modeling Using Logistic Regression and GIS: An Iranian Case Study. Small-Scale For. 2014, 13, 117–125. [Google Scholar] [CrossRef]
  54. Tien Bui, D.; Le, K.-T.T.; Nguyen, V.C.; Le, H.D.; Revhaug, I. Tropical forest fire susceptibility mapping at the Cat Ba National Park Area, Hai Phong City, Vietnam, using GIS-based kernel logistic regression. Remote Sens. 2016, 8, 347. [Google Scholar] [CrossRef] [Green Version]
  55. Vasilakos, C.; Kalabokidis, K.; Hatzopoulos, J.; Matsinos, I. Identifying wildland fire ignition factors through sensitivity analysis of a neural network. Nat. Hazards 2009, 50, 125–143. [Google Scholar] [CrossRef]
  56. Pierce, A.D.; Farris, C.A.; Taylor, A.H. Use of random forests for modeling and mapping forest canopy fuels for fire behavior analysis in Lassen Volcanic National Park, California, USA. For. Ecol. Manag. 2012, 279, 77–89. [Google Scholar] [CrossRef]
  57. Mirdeilami, T.; Shataee, S.H.; Kavoosi, M.R. Forest fire risk zone mapping in the Golestan national park using regression logistic method. J. Wood For. Sci. Technol. 2015, 22, 1–16. [Google Scholar]
  58. Srivastava, P.K.; Petropoulos, G.P.; Gupta, M.; Singh, S.K.; Islam, T.; Loka, D. Deriving forest fire probability maps from the fusion of visible/infrared satellite data and geospatial data mining. Model. Earth Syst. Environ. 2019, 5, 627–643. [Google Scholar] [CrossRef]
  59. Salar, S.G. Geomorphic Analysis for Water Harvesting Using GIS Technique in Selected Basins/Garmiyan Iraqi Kurdistan Region; University of Sulaimani: Sulaymaniyah, Iraq, 2013. [Google Scholar]
  60. Salar, S.G.; Othman, A.A.; Hasan, S.E. Identification of suitable sites for groundwater recharge in Awaspi watershed using GIS and remote sensing techniques. Environ. Earth Sci. 2018, 77, 701. [Google Scholar] [CrossRef]
  61. Kummerow, C.; Barnes, W.; Kozu, T.; Shiue, J.; Simpson, J. The Tropical Rainfall Measuring Mission (TRMM) sensor package. J. Atmos. Ocean. Technol. 1998, 15, 809–817. [Google Scholar] [CrossRef]
  62. NCEP Climate Forecast System Reanalysis (CFSR) Global Weather Data. Available online: https://swat.tamu.edu/data/cfsr (accessed on 13 August 2021).
  63. Rouse, J.W.; Haas, R.H.; Schelle, J.A.; Deering, D.W.; NASA/GSFC. Monitoring the Vernal Advancement or Retrogradation of Natural Vegetation; NASA: Greenbelt, MD, USA, 1974.
  64. Moore, I.D.; Grayson, R.B.; Ladson, A.R. Digital terrain modelling: A review of hydrological, geomorphological, and biological applications. Hydrol. Process. 1991, 5, 3–30. [Google Scholar] [CrossRef]
  65. Weiss, A. Topographic position and landforms analysis, Poster presentation. In Proceedings of the ESRI user conference, San Diego, CA, USA, 9–13 July 2001; pp. 9–13. [Google Scholar]
  66. Zeng, T.; Hudson, J.; Kay, S.; Laginestra, E.; Authority, S.O.P. A fuzzy GIS approach to fire risk assessment: A case study of Sydney Olympic Park, Australia. In Proceedings of the Spatial Sciences Conferences, Sydney, Australia, January 2003; Volume 18. [Google Scholar]
  67. Sayad, Y.O.; Mousannif, H.; Al Moatassime, H. Predictive modeling of wildfires: A new dataset and machine learning approach. Fire Saf. J. 2019, 104, 130–146. [Google Scholar] [CrossRef]
  68. Koutsias, N.; Xanthopoulos, G.; Founda, D.; Xystrakis, F.; Nioti, F.; Pleniou, M.; Mallinis, G.; Arianoutsou, M. On the relationships between forest fires and weather conditions in Greece from long-term national observations (1894–2010). Int. J. Wildl. Fire 2012, 22, 493–507. [Google Scholar] [CrossRef] [Green Version]
  69. Faramarzi, H.; Hosseini, S.M.; Gholamalifard, M. Fire hazard zoning in national golestan park using logistic regression and GIS. J. Geogr. Environ. Hazards 2014, 3, 73–90. [Google Scholar]
  70. Sakr, G.E.; Elhajj, I.H.; Mitri, G. Efficient forest fire occurrence prediction for developing countries using two weather parameters. Eng. Appl. Artif. Intell. 2011, 24, 888–894. [Google Scholar] [CrossRef]
  71. Rajabi, M.; Alesheikh, A.; Chehreghan, A.; Gazmeh, H. An innovative method for forest fire risk zoning map using fuzzy inference system and GIS. Int. J. Sci. Technol. Res. 2013, 2, 57–64. [Google Scholar]
  72. Saklani, P. Forest Fire Risk Zonation: A Case Study Pauri Garhwal; ITC: Uttrakhand, India, 2008.
  73. Behzadi, H.; Mohtashamnia, S.; Gharedaghi, H. Rangelands and forests fire risk zoning using GIS and AHP model (case study: Bamu National Park, Fars province). Iran. J. Range Desert Res. 2019, 25, fa817–fa827. [Google Scholar]
  74. Jenness, J.S. The Effects of Fire on Mexican Spotted Owls in Arizona and New Mexico. 2000. Available online: http://www.jennessent.com/Literature/Thesis/literature_review.htm#Forest%20Service%20Management%20Plans (accessed on 3 February 2022).
  75. Othman, A.A.; Al-Maamar, A.F.; Al-Manmi, D.A.M.; Veraldo, L.; Hasan, S.E.; Obaid, A.K.; Al-Quraishi, A.M.F. GIS-based modeling for selection of dam sites in the Kurdistan Region, Iraq. ISPRS Int. J. Geo-Inf. 2020, 9, 244. [Google Scholar] [CrossRef] [Green Version]
  76. Othman, A.A.; Obaid, A.K.; Amin Al-Manmi, D.A.M.; Al-Maamar, A.F.; Hasan, S.E.; Liesenberg, V.; Shihab, A.T.; Al-Saady, Y.I. New insight on soil loss estimation in the northwestern region of the Zagros fold and Thrust Belt. ISPRS Int. J. Geo-Inf. 2021, 10, 59. [Google Scholar] [CrossRef]
  77. Saha, S.; Moorthi, S.; Pan, H.-L.; Wu, X.; Wang, J.; Nadiga, S.; Tripp, P.; Kistler, R.; Woollen, J.; Behringer, D. The NCEP climate forecast system reanalysis. Bull. Am. Meteorol. Soc. 2010, 91, 1015–1058. [Google Scholar] [CrossRef]
  78. Sissakian, V.K.; Fouad, S.F.A. The Geology of Sulaimaniya Quadrangle Sheet No. NI-38-3, Scale 1:250,000; GEOSURV: Baghdad, Iraq, 2008. [Google Scholar]
  79. Al-Rubaiay, A.T.; Al-Dulaimi, T.Y. Series of Land Use Land Cover Maps of Iraq Scale 1:250,000, Sulaimaniya Quadrangle Sheet NI–38–3 (LULCM 10); GEOSURV: Baghdad, Iraq, 2012. [Google Scholar]
  80. Heydari, M.; Poorbabaei, H.; Esmaaelzade, O. Indicator plant species in monitoring forest soil conditions using logistic regression model in Zagros Oak (Quercus brantii var. persica) forest ecosystems, Ilam city. J. Plant Res. 2015, 27, 811–828. [Google Scholar]
  81. Eastman, J.R. IDRISI Taiga Guide to GIS and Image Processing; Clark Labs Clark University: Worcester, MA, USA, 2009. [Google Scholar]
  82. Pontius, R.G., Jr.; Schneider, L.C. Land-cover change model validation by an ROC method for the Ipswich watershed, Massachusetts, USA. Agric. Ecosyst. Environ. 2001, 85, 239–248. [Google Scholar] [CrossRef]
  83. Othman, A.A.; Gloaguen, R.; Andreani, L.; Rahnama, M. Improving landslide susceptibility mapping using morphometric features in the Mawat area, Kurdistan Region, NE Iraq: Comparison of different statistical models. Geomorphology 2018, 319, 147–160. [Google Scholar] [CrossRef]
  84. Youssef, A.M.; Pourghasemi, H.R.; Pourtaghi, Z.S.; Al-Katheeri, M.M. Landslide susceptibility mapping using random forest, boosted regression tree, classification and regression tree, and general linear models and comparison of their performance at Wadi Tayyah Basin, Asir Region, Saudi Arabia. Landslides 2016, 13, 839–856. [Google Scholar] [CrossRef]
  85. Jaafari, A.; Najafi, A.; Rezaeian, J.; Sattarian, A.; Ghajar, I. Planning road networks in landslide-prone areas: A case study from the northern forests of Iran. Land Use Policy 2015, 47, 198–208. [Google Scholar] [CrossRef]
  86. Shrestha, R.; Di, L.; Yu, E.G.; Kang, L.; SHAO, Y.; BAI, Y. Regression model to estimate flood impact on corn yield using MODIS NDVI and USDA cropland data layer. J. Integr. Agric. 2017, 16, 398–407. [Google Scholar] [CrossRef] [Green Version]
  87. Loveridge, M.; Rahman, A. Effects of Probability-Distributed Losses on Flood Estimates Using Event-Based Rainfall-Runoff Models. Water 2021, 13, 2049. [Google Scholar] [CrossRef]
  88. Bergmeir, C.; Benítez, J.M. On the use of cross-validation for time series predictor evaluation. Inf. Sci. 2012, 191, 192–213. [Google Scholar] [CrossRef]
  89. Bennett, M.; Fitzgerald, S.A.; Parker, B.; Main, M.L.; Perleberg, A.; Schnepf, C.; Mahoney, R.L. Reducing Fire Risk on Your Forest Property; Oregon State University: Corvallis, OR, USA, 2010. [Google Scholar]
  90. Lecina-Diaz, J.; Alvarez, A.; Retana, J. Extreme fire severity patterns in topographic, convective and wind-driven historical wildfires of Mediterranean pine forests. PLoS ONE 2014, 9, e85127. [Google Scholar] [CrossRef] [Green Version]
  91. Oliveira, S.; Pereira, J.M.C.; San-Miguel-Ayanz, J.; Lourenço, L. Exploring the spatial patterns of fire density in Southern Europe using Geographically Weighted Regression. Appl. Geogr. 2014, 51, 143–157. [Google Scholar] [CrossRef]
  92. Adab, H.; Kanniah, K.D.; Solaimani, K. Modelling static fire hazard in a semi-arid region using frequency analysis. Int. J. Wildl. Fire 2015, 24, 763–777. [Google Scholar] [CrossRef]
  93. Li, R.; Wang, N. Landslide susceptibility mapping for the Muchuan county (China): A comparison between bivariate statistical models (woe, ebf, and ioe) and their ensembles with logistic regression. Symmetry 2019, 11, 762. [Google Scholar] [CrossRef] [Green Version]
  94. Gholami, P.; Ghorbani, J.; Abbasi, H. Effect of fire vegetation on some properties of soil in rangelands of Bamo national park in Shiraz. Nat. Ecosyst. Iran 2014, 5, 41–50. [Google Scholar]
  95. Rafiee, F.; Jankju, M.; Ejtehadi, H. Investigation on tolerant, adapted and sensitive plant traits to chronological wildfires in a semiarid rangeland. Iran. J. Range Desert Res. 2015, 22, Pe73–Pe84. [Google Scholar]
  96. Oliveira, C.; Ferreira, R.; da Silva, J.; Lima, R.; Silva, E.; Silva, A.; Lucena, J.; dos Santos, N.; Lopes, I.; Pessoa, M.; et al. Modeling and Spatialization of Biomass and Carbon Stock Using LiDAR Metrics in Tropical Dry Forest, Brazil. Forests 2021, 12, 473. [Google Scholar] [CrossRef]
  97. Fernández-Guisuraga, J.M.; Calvo, L.; Fernandes, P.M.; Suárez-Seoane, S. Short-Term Recovery of the Aboveground Carbon Stock in Iberian Shrublands at the Extremes of an Environmental Gradient and as a Function of Burn Severity. Forests 2022, 13, 145. [Google Scholar] [CrossRef]
  98. Zhang, F.; Zhou, G. Estimation of vegetation water content using hyperspectral vegetation indices: A comparison of crop water indicators in response to water stress treatments for summer maize. BMC Ecol. 2019, 19, 18. [Google Scholar] [CrossRef] [Green Version]
  99. González-Olabarria, J.-R.; Rodriguez, F.; Fernández-Landa, A.; Mola-Yudego, B. Mapping fire risk in the Model Forest of Urbión (Spain) based on airborne LiDAR measurements. For. Ecol. Manag. 2012, 282, 149–156. [Google Scholar] [CrossRef]
  100. Reddy, C.S.; Bird, N.G.; Sreelakshmi, S.; Manikandan, T.M.; Asra, M.; Krishna, P.H.; Jha, C.S.; Rao, P.V.N.; Diwakar, P.G. Identification and characterization of spatio-temporal hotspots of forest fires in South Asia. Environ. Monit. Assess. 2019, 191, 791. [Google Scholar] [CrossRef]
  101. Avitabile, S.C.; Callister, K.E.; Kelly, L.T.; Haslem, A.; Fraser, L.; Nimmo, D.G.; Watson, S.J.; Kenny, S.A.; Taylor, R.S.; Spence-Bailey, L.M.; et al. Systematic fire mapping is critical for fire ecology, planning and management: A case study in the semi-arid Murray Mallee, south-eastern Australia. Landsc. Urban Plan. 2013, 117, 81–91. [Google Scholar] [CrossRef]
  102. Bernier, P.Y.; Gauthier, S.; Jean, P.-O.; Manka, F.; Boulanger, Y.; Beaudoin, A.; Guindon, L. Mapping Local Effects of Forest Properties on Fire Risk across Canada. Forests 2016, 7, 157. [Google Scholar] [CrossRef]
  103. Carvalheiro, L.C.; Bernardo, S.O.; Orgaz, M.D.M.; Yamazaki, Y. Forest Fires Mapping and Monitoring of current and past forest fire activity from Meteosat Second Generation Data. Environ. Model. Softw. 2010, 25, 1909–1914. [Google Scholar] [CrossRef]
  104. Chen, D.; Pereira, J.M.C.; Masiero, A.; Pirotti, F. Mapping fire regimes in China using MODIS active fire and burned area data. Appl. Geogr. 2017, 85, 14–26. [Google Scholar] [CrossRef]
  105. Dlamini, W.M. Application of Bayesian networks for fire risk mapping using GIS and remote sensing data. GeoJournal 2011, 76, 283–296. [Google Scholar] [CrossRef]
  106. Duarte, L.; Teododo, A.C. An easy, accurate and efficient procedure to create forest fire risk maps using the SEXTANTE plugin Modeler. J. For. Res. 2016, 27, 1361–1372. [Google Scholar] [CrossRef]
  107. Hudak, A.T.; Bright, B.C.; Pokswinski, S.M.; Loudermilk, E.L.; O’Brien, J.J.; Hornsby, B.S.; Klauberg, C.; Silva, C.A. Mappingforest structure and composition from low-density LiDAR for informed forest, fuel, and fire management at Eglin Air Force Base, Florida, USA. Can. J. Remote Sens. 2016, 42, 411–427. [Google Scholar] [CrossRef]
  108. Kumar, S.; Das Bairagi, G.; Kumar, A. Identifying triggers for forest fire and assessing fire susceptibility of forests in Indian western Himalaya using geospatial techniques. Nat. Hazards 2015, 78, 203–217. [Google Scholar] [CrossRef]
  109. Lampin-Maillet, C.; Jappiot, M.; Long, M.; Bouillon, C.; Morge, D.; Ferrier, J.-P. Mapping wildland-urban interfaces at large scales integrating housing density and vegetation aggregation for fire prevention in the South of France. J. Environ. Manag. 2010, 91, 732–741. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  110. Merino-de-Miguel, S.; Huesca, M.; González-Alonso, F. Modis reflectance and active fire data for burn mapping and assessment at regional level. Ecol. Modell. 2010, 221, 67–74. [Google Scholar] [CrossRef]
  111. Mitri, G.H.; Gitas, I.Z. Mapping post-fire forest regeneration and vegetation recovery using a combination of very high spatial resolution and hyperspectral satellite imagery. Int. J. Appl. Earth Obs. Geoinf. 2013, 20, 60–66. [Google Scholar] [CrossRef]
  112. Molina, J.R.; González-Cabán, A.; y Silva, F.R. Wildfires impact on the economic susceptibility of recreation activities: Application in a Mediterranean protected area. J. Environ. Manag. 2019, 245, 454–463. [Google Scholar] [CrossRef]
  113. Ruiz, J.A.M.; Riaño, D.; Arbelo, M.; French, N.H.F.; Ustin, S.L.; Whiting, M.L. Burned area mapping time series in Canada (1984--1999) from NOAA-AVHRR LTDR: A comparison with other remote sensing products and fire perimeters. Remote Sens. Environ. 2012, 117, 407–414. [Google Scholar] [CrossRef]
  114. Morton, D.C.; DeFries, R.S.; Nagol, J.; Souza, C.M., Jr.; Kasischke, E.S.; Hurtt, G.C.; Dubayah, R. Mapping canopy damage from understory fires in Amazon forests using annual time series of Landsat and MODIS data. Remote Sens. Environ. 2011, 115, 1706–1720. [Google Scholar] [CrossRef] [Green Version]
  115. Naderpour, M.; Rizeei, H.M.; Khakzad, N.; Pradhan, B. Forest fire induced Natech risk assessment: A survey of geospatial technologies. Reliab. Eng. Syst. Saf. 2019, 191, 106558. [Google Scholar] [CrossRef]
  116. Nandy, S.; Singh, C.; Das, K.K.; Kingma, N.C.; Kushwaha, S.P.S. Environmental vulnerability assessment of eco-development zone of Great himalayan National Park, himachal Pradesh, India. Ecol. Indic. 2015, 57, 182–195. [Google Scholar] [CrossRef]
  117. Potapov, P.; Turubanova, S.; Hansen, M.C. Regional-scale boreal forest cover and change mapping using Landsat data composites for European Russia. Remote Sens. Environ. 2011, 115, 548–561. [Google Scholar] [CrossRef]
  118. Sahana, M.; Ganaie, T.A. GIS-based landscape vulnerability assessment to forest fire susceptibility of Rudraprayag district, Uttarakhand, India. Environ. Earth Sci. 2017, 76, 676. [Google Scholar] [CrossRef]
  119. Saranya, K.R.L.; Reddy, C.S.; Rao, P.V.V.P.; Jha, C.S. Decadal time-scale monitoring of forest fires in Similipal Biosphere Reserve, India using remote sensing and GIS. Environ. Monit. Assess 2014, 186, 3283–3296. [Google Scholar] [CrossRef]
  120. Weisberg, P.J.; Ko, D.; Py, C.; Bauer, J.M. Modeling fire and landform influences on the distribution of old-growth pinyon-juniper woodland. Landsc. Ecol. 2008, 23, 931–943. [Google Scholar] [CrossRef]
  121. Whitman, E.; Rapaport, E.; Sherren, K. Modeling fire susceptibility to delineate wildland--urban interface for municipal-scale fire risk management. Environ. Manag. 2013, 52, 1427–1439. [Google Scholar] [CrossRef]
  122. Yang, W.; Zhang, S.; Tang, J.; Bu, K.; Yang, J.; Chang, L. A MODIS time series data based algorithm for mapping forest fire burned area. Chin. Geogr. Sci. 2013, 23, 344–352. [Google Scholar] [CrossRef]
  123. You, W.; Lin, L.; Wu, L.; Ji, Z.; Zhu, J.; Fan, Y.; He, D. Geographical information system-based forest fire risk assessment integrating national forest inventory data and analysis of its spatiotemporal variability. Ecol. Indic. 2017, 77, 176–184. [Google Scholar] [CrossRef]
  124. Ziccardi, L.G.; Thiersch, C.R.; Yanai, A.M.; Fearnside, P.M.; Ferreira-Filho, P.J. Forest fire risk indices and zoning of hazardous areas in Sorocaba, São Paulo state, Brazil. J. For. Res. 2020, 31, 581–590. [Google Scholar] [CrossRef]
  125. Cao, Y.; Wang, M.; Liu, K. Wildfire susceptibility assessment in Southern China: A comparison of multiple methods. Int. J. Disaster Risk Sci. 2017, 8, 164–181. [Google Scholar] [CrossRef] [Green Version]
  126. Pourghasemi, H.R.; Beheshtirad, M.; Pradhan, B. A comparative assessment of prediction capabilities of modified analytical hierarchy process (M-AHP) and Mamdani fuzzy logic models using Netcad-GIS for forest fire susceptibility mapping. Geomat. Nat. Hazards Risk 2016, 7, 861–885. [Google Scholar] [CrossRef] [Green Version]
  127. Eskandari, S. A new approach for forest fire risk modeling using fuzzy AHP and GIS in Hyrcanian forests of Iran. Arab. J. Geosci. 2017, 10, 190. [Google Scholar] [CrossRef]
  128. Li, X.; Lanorte, A.; Telesca, L.; Song, W.; Lasaponara, R. Assessment of MODIS-Based NDVI-Derived Index for Fire Susceptibility Estimation in Northern China. In Proceedings of the International Conference on Computational Science and Its Applications, Banff, AB, Canada, 22–25 June 2015; pp. 193–203. [Google Scholar]
  129. Liu, W.; Wang, S.; Zhou, Y.; Wang, L.; Zhu, J.; Wang, F. Lightning-caused forest fire risk rating assessment based on case-based reasoning: A case study in DaXingAn Mountains of China. Nat. Hazards 2016, 81, 347–363. [Google Scholar] [CrossRef]
  130. Liu, Q.; Shan, Y.; Shu, L.; Sun, P.; Du, S. Spatial and temporal distribution of forest fire frequency and forest area burnt in Jilin Province, Northeast China. J. For. Res. 2018, 29, 1233–1239. [Google Scholar] [CrossRef]
  131. Tehrany, M.S.; Jones, S.; Shabani, F.; Martínez-Álvarez, F.; Bui, D.T. A novel ensemble modeling approach for the spatial prediction of tropical forest fire susceptibility using LogitBoost machine learning classifier and multi-source geospatial data. Arch. Meteorol. Geophys. Bioclimatol. Ser. B 2018, 137, 637–653. [Google Scholar] [CrossRef]
Figure 1. Location of the study area in the southern part of the Sulaymaniyah governate in the Kurdistan Region of Iraq. The total burned area mapped is shown for the entire region.
Figure 1. Location of the study area in the southern part of the Sulaymaniyah governate in the Kurdistan Region of Iraq. The total burned area mapped is shown for the entire region.
Sustainability 14 06194 g001
Figure 2. Flow chart of used methodology in vegetation fire mapping.
Figure 2. Flow chart of used methodology in vegetation fire mapping.
Sustainability 14 06194 g002
Figure 3. Input vegetation fire factors used in the logistic regression method: (A) slope gradient; (B) slope aspect; (C) elevation; (D) distance to the road; (E) distance to the settlements; (F) land use and land cover; (G) wind speed (from 1979 to 2014); (H) NDVI; (I) distance to the streams; (J) geology; (K) TWI; (L) relative humidity (from 1979 to 2014); (M) distance to the farmland; (N) TPI.
Figure 3. Input vegetation fire factors used in the logistic regression method: (A) slope gradient; (B) slope aspect; (C) elevation; (D) distance to the road; (E) distance to the settlements; (F) land use and land cover; (G) wind speed (from 1979 to 2014); (H) NDVI; (I) distance to the streams; (J) geology; (K) TWI; (L) relative humidity (from 1979 to 2014); (M) distance to the farmland; (N) TPI.
Sustainability 14 06194 g003aSustainability 14 06194 g003bSustainability 14 06194 g003c
Figure 4. Meteorological parameter maps: (A) annual precipitation (from 1998 to 2017); (B) annual temperature (from 1979 to 2014); (C,D) coefficient of determination maximum and minimum temperature, respectively.
Figure 4. Meteorological parameter maps: (A) annual precipitation (from 1998 to 2017); (B) annual temperature (from 1979 to 2014); (C,D) coefficient of determination maximum and minimum temperature, respectively.
Sustainability 14 06194 g004
Figure 5. (A) fire susceptibility map for group 1; (B) fire susceptibility map for group 2; (C) fire susceptibility map for group 3; (D) fire susceptibility map for group 4; (E) fire susceptibility map for group 5; (F) fire susceptibility map for group 6.
Figure 5. (A) fire susceptibility map for group 1; (B) fire susceptibility map for group 2; (C) fire susceptibility map for group 3; (D) fire susceptibility map for group 4; (E) fire susceptibility map for group 5; (F) fire susceptibility map for group 6.
Sustainability 14 06194 g005
Figure 6. Best fire susceptibility map.
Figure 6. Best fire susceptibility map.
Sustainability 14 06194 g006
Figure 7. Vegetation fire prediction rate curve, for both (A) 20% and (B) 80%, as a validation dataset.
Figure 7. Vegetation fire prediction rate curve, for both (A) 20% and (B) 80%, as a validation dataset.
Sustainability 14 06194 g007
Table 1. Sentinel-2 spectral bands used to determine the NDVI.
Table 1. Sentinel-2 spectral bands used to determine the NDVI.
BandsCentral Wavelength (µm)Resolution (m)
Band 4—Red0.66510
Band 8—NIR0.84210
Table 2. Factors and sub-factors effective in fires as well as their source.
Table 2. Factors and sub-factors effective in fires as well as their source.
FactorsThe Method of Affecting FiresFactors Used by Experts (%)Source
1.Slope gradientFire spread is higher on steep slopes, so that with a 10-degree increase in the slope angle, the propagation speed doubles [66].54.24DEM-30 m
https://earthexplorer.usgs.gov/ (accessed on 13 August 2021)
2.Slope aspectSince southward slopes are warmer and drier than other areas, fire risks are higher on these slopes [67]. 52.54DEM-30 m
https://earthexplorer.usgs.gov/ (accessed on 13 August 2021)
3.ElevationUnder normal conditions, at lower and middle altitudes of an area, fires are much more probable to occur than at higher altitudes because of relatively higher temperatures, lower humidity, and ease of human access [12]. Besides, by an increase in the altitude, damage caused by fires is reduced due to the slower growth of trees compared to lands of lower altitudes and the lower accumulation of resin under the bark [68].50.85DEM-30 m
https://earthexplorer.usgs.gov/ (accessed on 13 August 2021)
4.PrecipitationPrecipitation is an important factor contributing to the high humidity of fuels, so it is considered a negative indicator of fire spread. In fact, the higher the moisture content in a species’ tissue, the more heat and the longer time it will need to evaporate the moisture content, so the higher fire resistance will be. Therefore, the moisture content of the burnable matter is one of the major factors effective in fire occurrence [16]. Besides, rising precipitation in spring increases vegetation growth. This vegetation dries in summer and makes fire occurrence possible [53]. Thus, it is required to pay attention to the amount of precipitation, the precipitation season, and the system creating it. This is because precipitation accompanied by strong lightning can cause fires.40.68TRMM-NASA
5.Distance from roadsWhen the distance from roads increases, fire risks decrease because of less traffic and human activity [47].40.68World Imagery-GIS base map
6.Distance from settlementsThere is a higher fire occurrence risk in the vicinity of settlements involving human activity [58].40.68World Imagery-GIS base map
7.LULCThe effects of LULC are exerted by human activity. Thus, it uses in which there is more human activity with suitable conditions for fire, it is more probable to occur [69].38.98LULC
https://earthexplorer.usgs.gov/ (accessed on 13 August 2021)
8.Annual temperatureTemperature rises increase evaporation and transpiration, thereby drying combustible materials, which is considered one of the factors effective in fire occurrence [70].38.98Climate Data
https://globalweather.tamu.edu/ (accessed on 13 August 2021)
9.Wind speedThe higher the wind speed is, the higher the fire intensity will be. This is because when the wind moves the air, more oxygen is transferred to the burning environment [22]. In this respect, if the wind blows from the land, it will have a greater effect on fire intensity.32.2Climate Data
https://globalweather.tamu.edu/ (accessed on 13 August 2021)
10.NDVISince vegetation is the fuel itself, and NDVI shows the condition of vegetation, this factor has a strong effect on fire ignition and the spread of fire [23].28.81Sentinel-2 A
https://earthexplorer.usgs.gov/ (accessed on 13 August 2021)
11.Distance from streamsHuman activity (particularly, camping along rivers and springs in summer, monitoring cropland, and, sometimes, visiting by tourists to this area, due to the river and the good weather), could affect vegetation fires. Thus, the increase in people’s presence along rivers is one of the serious threats [71].23.73World Imagery-GIS base map
12.GeologyThe parent rock determines the soil type. Soil is considered one of the major parameters in describing vegetation fires. Besides, it, indirectly, affects the entire environment of a given region [72]. The soil type demonstrates the effects of the texture and composition of soil substances on fire occurrence [15].16.95Iraqi Geological Survey (Scale 1: 250,000) GEOSURV-IRAQ
13.TWIThe TWI shows the size of saturated areas of runoff generation and the effects of topography on the location. Thus, upon an increase in the value of this index, fire risks decrease [15].16.95DEM-30 m
https://earthexplorer.usgs.gov/ (accessed on 13 August 2021)
14.Relative humidityIn dry seasons, with an increase in the temperature as well as a decrease in relative humidity and dryness of vegetation, the probability of fire occurrence and spread increases [73].13.56Climate Data
https://globalweather.tamu.edu/ (accessed on 13 August 2021)
15.Distance from farmlandsSince many fires are created for clearing agricultural lands off crop residues or developing them, fire risks are greater near agricultural lands [28]. Traditionally, farmers used the above-mentioned method, which is not seen as a harmful method. On the contrary, it is viewed as a method that increases soil fertility. Nevertheless, the fire sometimes gets out of control, unintentionally.11.86World Imagery-GIS base map
16.TPIIt displays the difference in height between a focal cell and all cells in the neighborhood [74]. In addition, TPI shows that flat areas are not as favorable for fire occurrence as ridges and gentle slopes [5]. Besides, the more positive or negative the curvature is, the more likely the vegetation-fire occurrence will be [4]. 10.17DEM-30 mhttps://earthexplorer.usgs.gov/ (accessed on 13 August 2021)
Table 3. Additional statistical validation methods.
Table 3. Additional statistical validation methods.
Validation MethodGroup 1Group 2Group 3Group 4Group 5Group 6Group 7
Mean Absolute Error (MAE)%1.191.010.682.540.860.760.62
Relative Error (RE)%5.735.412.6915.153.803.282.32
Percentage Error (PE)%1.261.070.702.990.890.790.64
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Salar, S.G.; Othman, A.A.; Rasooli, S.; Ali, S.S.; Al-Attar, Z.T.; Liesenberg, V. GIS-Based Modeling for Vegetated Land Fire Prediction in Qaradagh Area, Kurdistan Region, Iraq. Sustainability 2022, 14, 6194. https://0-doi-org.brum.beds.ac.uk/10.3390/su14106194

AMA Style

Salar SG, Othman AA, Rasooli S, Ali SS, Al-Attar ZT, Liesenberg V. GIS-Based Modeling for Vegetated Land Fire Prediction in Qaradagh Area, Kurdistan Region, Iraq. Sustainability. 2022; 14(10):6194. https://0-doi-org.brum.beds.ac.uk/10.3390/su14106194

Chicago/Turabian Style

Salar, Sarkawt G., Arsalan Ahmed Othman, Sabri Rasooli, Salahalddin S. Ali, Zaid T. Al-Attar, and Veraldo Liesenberg. 2022. "GIS-Based Modeling for Vegetated Land Fire Prediction in Qaradagh Area, Kurdistan Region, Iraq" Sustainability 14, no. 10: 6194. https://0-doi-org.brum.beds.ac.uk/10.3390/su14106194

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop