Next Article in Journal
Application and Analysis of XCO2 Data from OCO Satellite Using a Synthetic DINEOF–BME Spatiotemporal Interpolation Framework
Next Article in Special Issue
Globally vs. Locally Trained Machine Learning Models for Landslide Detection: A Case Study of a Glacial Landscape
Previous Article in Journal
Semi-Supervised DEGAN for Optical High-Resolution Remote Sensing Image Scene Classification
Previous Article in Special Issue
Landslide Displacement Prediction Based on a Two-Stage Combined Deep Learning Model under Small Sample Condition
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Creation of Wildfire Susceptibility Maps in Plumas National Forest Using InSAR Coherence, Deep Learning, and Metaheuristic Optimization Approaches

1
Division of Science Education, Kangwon National University, Chuncheon-si 24341, Korea
2
Department of Civil and Environmental Engineering, Lamar University, 4400 MLK Blvd., Beaumont, TX 77710, USA
3
Department of Smart Regional Innovation, Kangwon National University, Chuncheon-si 24341, Korea
*
Author to whom correspondence should be addressed.
Remote Sens. 2022, 14(17), 4416; https://doi.org/10.3390/rs14174416
Submission received: 13 July 2022 / Revised: 31 August 2022 / Accepted: 1 September 2022 / Published: 5 September 2022
(This article belongs to the Special Issue Machine Learning and Remote Sensing for Geohazards)

Abstract

:
Plumas National Forest, located in the Butte and Plumas counties, has experienced devastating wildfires in recent years, resulting in substantial economic losses and threatening the safety of people. Mapping damaged areas and assessing wildfire susceptibility are necessary to prevent, mitigate, and manage wildfires. In this study, a wildfire susceptibility map was generated using a CNN and metaheuristic optimization algorithms (GWO and ICA) based on images of areas damaged by wildfires. The locations of damaged areas were identified using the damage proxy map (DPM) technique from Sentinel-1 synthetic aperture radar (SAR) data collected from 2016 to 2020. The DPMs’ depicting areas damaged by wildfires were similar to fire perimeters obtained from the California Department of Forestry and Fire Protection (CAL FIRE). Data regarding damaged areas were divided into a training set (50%) for modeling and a testing set (50%) for assessing the accuracy of the models. Sixteen conditioning factors, categorized as topographical, meteorological, environmental, and anthropological factors, were selected to construct the models. The wildfire susceptibility models were evaluated using the area under the receiver operating characteristic (ROC) curve (AUC) and root mean square error (RMSE) analysis. The evaluation results revealed that the hybrid-based CNN-GWO model (AUC = 0.974, RMSE = 0.334) exhibited better performance than the CNN (AUC = 0.934, RMSE = 0.780) and CNN-ICA (AUC = 0.950, RMSE = 0.350) models. Therefore, we conclude that optimizing a CNN with metaheuristics considerably increased the accuracy and reliability of wildfire susceptibility mapping in the study area.

Graphical Abstract

1. Introduction

A wildfire is any uncontrolled fire that spreads through vegetative fuels and threatens to destroy life, property, or resources [1]. In California, wildfires are one of the most common natural disasters and have caused significant harm to the environment, society, and economy in recent years [2]. Particularly in the last two decades, climate change and land utilization for human activities have exacerbated wildfires [3]. In addition, sustained population growth and rapid social development have led to the expansion of areas of the wildland–urban interface (WUI) and have seriously increased the number of individuals and buildings affected by wildfires, which has further exacerbated their negative impact on human life [4]. The number of acres burned, number of fires, and amount of property damaged due to wildfires have not substantially decreased in California, despite significant administrative investment in wildfire suppression and management, according to the Wildfire Redbook published by the California Department of Forestry and Fire Protection (CAL FIRE).
According to wildfire activity statistics from CAL FIRE in 2020, wildfire agencies responded to an average of 8198 fires that burned 7263 km2 per year from 2016 to 2020 [1]. Three of the deadliest, largest, and most destructive California wildfires occurred in Butte and Plumas County [5]. The Dixie Fire was caused by powerlines and burned 3898 km2, killing one person and damaging 1329 structures in July 2021. The North Complex Fire was caused by lightning and burned 1291 km2, took 15 lives, and destroyed 2352 structures in August 2020. Another devastating wildfire that occurred in November 2018 was the Camp Fire, which burned 621 km2, claimed 85 lives, and destroyed 18,804 structures. The secondary effects of wildfires include health problems, reduced air quality, and economic losses due to wildfire evacuation [6,7,8]. In terms of long-term impacts, wildfires discharge aerosols and produce carbon emissions that affect the climate and climate change-related phenomena that drive even more wildfires [9].
Taking into account the various impacts caused by the burning and spreading of wildfires, collecting data on past and current wildfires, including their start date, location, amount of area burned, and duration, serves as valuable information for developing planning programs for the prevention and response to wildfires [10]. These inventories are an essential source of information for creating a hazard, risk, and susceptibility map. The identification of wildfire-prone areas can be achieved through susceptibility mapping [11,12,13]. Susceptibility is the likelihood of a wildfire in the future in the spatial dimension, with the predisposition to burn for each spatial unit (grid cell) assessed based on terrain, meteorological, and anthropogenic features such as topography and land cover, with the risk ranked from low to high [10]. This study evaluated past and recent wildfires in particularly damaged areas and determined the anthropogenic and geo-environmental predisposing factors that increase wildfire risk.
Providing accurate data regarding the location of previous wildfires is a critical initial step in susceptibility mapping. Wildfire inventory maps can be generated by several methods, including extensive global positioning system (GPS) field surveys [12,14] or remote sensing data [11,15]. Moderate-resolution imaging spectroradiometer (MODIS) products are freely accessible and are widely employed for generating datasets for wildfire-affected areas [11,14,15,16]. However, the spatial resolution of the hotspot used for identifying wildfires is 1 km. Sulova et al. (2020) combined a MODIS dataset and a Sentinel-2 mission to obtain the locations of wildfires with a spatial resolution of 20 m. Sentinel-2 is an optical satellite with 13 spectral bands that can be utilized for fuel condition mapping (e.g., vegetation greenness, fuel moisture) and burn severity mapping. However, its function depends on the atmospheric conditions of the areas studied. Radar satellite has several advantages over optical sensors: clouds and smoke are transparent to radar signals; radar does not require sunlight; and due to its coherent character, radar signal has the ability to detect minor alterations in surface property changes. The potential of synthetic aperture radar (SAR) images for mapping burnt areas lies in the sensitivity of SAR backscattering to vegetation structure and biomass, and the changes in scattering modes caused by fire events. The InSAR technique can make a direct measurement on the decorrelation among different acquisition dates by integrating both amplitude and phase information [17]. For these reasons, this study utilized the SAR dataset and the damage proxy map (DPM) method to identify areas of Plumas National Forest that have been likely damaged by wildfires and to generate an inventory dataset. The DPM method was developed by a joint project between the California Institute of Technology and the Jet Propulsion Laboratory (JPL), known as the Advanced Rapid Imaging and Analysis (ARIA) program. Using the DPM method is beneficial for damage mapping following natural disasters and has been applied to events including the 2017 Pohang earthquake in South Korea [18], Typhoon Hagibis in Japan [19], and the 2014 eruption of the Kelud volcano in Indonesia [20].
Various probabilistic and statistical approaches have been used to assess wildfire susceptibility in recent decades, including the frequency ratio (FR) [21,22], Bayesian theory [23,24], analytic hierarchy process (AHP) [25,26], and logistic regression (LR) [27,28]. Machine learning has also gained popularity for generating wildfire susceptibility maps, including artificial neural networks (ANNs) [29], support vector machines (SVMs) [30,31], decision tree-based algorithms [12], and random forests [11,13]. These machine learning methods show great potential for evaluating non-linear and multivariate datasets. However, previous machine learning methods could not uncover more representative features from the input data to increase the precision of the forecasting process. Therefore, a new robust method has been developed to address this limitation by generating wildfire susceptibility mappings using deep learning algorithms, for instance, deep neural networks (DNNs) [32] and convolutional neural networks (CNNs) [33,34]. The CNN algorithm has a greater ability to recognize high variability patterns compared with other deep learning algorithms [35]. However, the choice of parameter settings has a substantial impact on the learning process, impacts the prediction results, and potentially leads to problems of overfitting or underfitting. Metaheuristic optimization based on iterative simulation is a widely used technique for addressing these problems; by optimizing the hyperparameters, it exhibits better performance than standalone machine learning and has been employed in the analysis of disasters such as landslides [36], subsidence [37], and floods [38].
This study aimed to produce wildfire susceptibility maps for Plumas National Forest by combining SAR data and deep learning methods based on CNN with metaheuristic optimization algorithms. Firstly, we generated DPMs after past wildfires occurring from 2016 to 2020, based on the Sentinel-1 SAR dataset as a dependent variable. The independent variables in the analysis were 16 conditioning factors from four categories: topographical, meteorological, anthropological, and environmental. Deep learning and metaheuristic optimization approaches were used in the analysis. Model performance was evaluated using root mean square error (RMSE) analysis and area under the receiver operating characteristic curve (AUC) analysis. Finally, high-risk and safe areas in Plumas National Forest were identified by creating maps based on the model.

2. Materials and Methods

Figure 1 shows the overall method that has been summarized in a graphical illustration. The first step is generating wildfire inventory utilizing the DPM and SAR datasets. The DPM-derived wildfire inventory data were then randomly separated into training (50%) and testing (50%) sets. After that, a spatial database was produced and then assessed employing spatial correlation analysis using the generated wildfire inventory and layers of the related factors.
Next, the hyperparameters of CNN were optimized using GWO and ICA metaheuristic optimization algorithms and wildfire susceptibility models were produced. Finally, the wildfire susceptibility maps were compared and evaluated using RMSE and AUC analysis.

2.1. Study Area

Plumas National Forest resides in 4638 km2 of mountain lands in the north of Sierra Nevada. Approximately 85% of the forest sits in Plumas County and the rest extends into Butte County, CA, USA, as shown in Figure 2. The forest was named after its main watershed, the Rio de las Plumas or Feather River. This area has a Mediterranean climate with hot, dry summers and cool, wet winters. Plumas National Forest has a complex topography and heterogeneous and dense vegetation (92.62% of its surface is covered by forests).
The Fire and Resources Assessment Program (FRAP) created fire perimeters and established a database to represent a comprehensive digital record of fire perimeters in California [1]. Figure 2 shows the distribution of historic fires in Plumas National Forest from 2016 to 2020. There were 77 wildfires in Plumas National Forest from 2016 to 2021, with an average of 13 wildfires burning 115,711 ha per year. Figure 3 shows the number of wildfires and burned areas from 2016 to 2021. The greatest number of wildfires occurred in 2020, with 23 fires, while the largest total burned area was 433,100 ha in 2021. Although the number of wildfires in Plumas National Forest is uncertain, the burned area is widening every year. In general, the peak season for wildfires in Plumas National Forest is late summer and early autumn. In term of frequency, from 2016–2020, the wildfire season started in April and the most frequent occurrence was observed in August.
CAL FIRE [1] also investigated and recorded the causes of wildfires and found that human causes (direct or indirect) ignited 38.57% of the wildfires in Plumas National Forest (as shown in the statistical summary in Table 1). Therefore, anthropological factors should be considered in constructing a wildfire susceptibility map. According to the report, the unknown ignition cause describes a fire that has been investigated and has insufficient information to classify further, a fire that is under investigation, or a fire that has not yet been investigated.

2.2. SAR Datasets

SAR is beneficial compared to other types of remote sensing because its radar transmission in the microwave spectrum is not affected by day and night visibility and weather conditions, including cloud cover or smoke from fires [17]. Sentinel-1 SAR C-band data (5.5 cm wavelength) from the European Space Agency (ESA) were utilized to depict damaged areas after a wildfire in Plumas National Forest. Plumas National Forest is located on path 42 and frame 462. We collected Sentinel-1 single look complex (SLC) data with interferometric wide-swath (IW) mode and vertical transmission and vertical return (VV) polarization. We collected three scenes of Sentinel-1 SLCs for each wildfire of more than 300 ha acres. Two images were obtained before the wildfire and one after the wildfire was contained. In total, we collected 33 scenes to depict the damaged areas after wildfires from 2016 to 2020 in Plumas National Forest (Table 2). The Sentinel-1 SLC level-1 data were downloaded from the Alaska Satellite Facility (ASF) where each scene was pre-processed in Gamma 2021 software, (https://www.gamma-rs.ch/uploads/media/upgrades_info_20211201.pdf (accessed on 12 July 2022)). Each scene was updated through orbit information update to reduce orbit error.

2.3. Wildfire Conditioning Factors

Selecting independent variables, also known as predictors, predisposing, or conditioning factors, is vital in predictive modeling. Based on prior studies in California [2,39,40], 16 factors related to wildfire susceptibility were selected and categorized into four categories as shown in Table 3: topographical, meteorological, environmental, and anthropological. The topographical factors included aspect, altitude, slope, and plan curvature; the meteorological factors included precipitation, maximum temperature, solar radiance, and windspeed; the environmental factors included distance to stream, drought index, soil moisture, NDVI, and topographic wetness index; and the anthropological factors included land use, distance to road, and distance to settlement. The 16 factors used for the wildfire susceptibility analysis are shown in Figure 4. All factors were organized into a raster-based spatial database with a spatial resolution of 30 m. Continuous (numeric) data were reclassified into five classes using the quantile method to identify and analyze the effect of wildfires.
The topographic conditioning factors were altitude, aspect, slope, and plan curvature, derived from the Copernicus digital elevation model (DEM) with a spatial resolution of 30 m. The Copernicus DEM is a product of radar satellite data collected during the TanDEM-X mission. Altitude is an important factor affecting the spread and severity of a wildfire and is associated with local climate variations, vegetation distribution, composition, and flammability [14,45]. Increasing the degree of inclination can increase the rate of fire spread. Fire can spread quickly into steep areas and less quickly down steep areas [11]. Aspect shows the direction the slope is facing and affects the amount of solar radiation received [15,40]
The meteorological conditioning factors were the mean values of precipitation, maximum temperature, solar radiation, and wind speed. These factors control the life cycle of flora and fauna, which contribute to producing fuels, drying fuels for ignition, or spreading wildfires [15]. Solar radiation was collected from the National Renewable Energy Laboratory (NREL) with a spatial resolution of 4 km [46]. Increasing solar radiation can reduce soil moisture and increase temperature, thus increasing wildfire risk [47]. Precipitation and maximum temperature were collected from parameter-elevation regressions on independent slopes model (PRISM) managed by Oregon State University’s PRISM Climate group. Precipitation affects the moisture level and vegetation pattern that influence the speed of fire distribution [48]. There is a direct correlation between increases in temperature and wildfires [14,41]. Windspeed was obtained from the Global Wind Atlas with a spatial resolution of 100 m. Windspeed is related to wildfires, as wind affects their speed and severity [49]. The rasters were generated from images collected from 2016 to 2020 using the mean statistic function [43].
The environmental factors were the distance to stream, soil moisture, drought index, NDVI, and topographic wetness index (TWI). Stream data were acquired from California State Geoportal. The distance to stream has a direct role in forest health, with the stream serving as a water resource [14,50]. The soil moisture and drought index were collected from the Terra Climate 2016–2020 dataset with a spatial resolution of 4 km. This monthly dataset for global terrestrial surfaces uses the interpolated time-varying anomalies from CRU Ts 4.0/JRA55, with high-spatial climatological normal obtained from the WorldClim dataset to create a high-spatial resolution dataset that covers a broader temporal record [51]. The drought index, based on the Palmer Drought Severity Index (PDSI), measures agricultural drought calculated by adding precipitation to the top two layers of the soil and employing a temperature-driven evapotranspiration algorithm to remove moisture [13]. Soil moisture directly influences the dryness of fuels and affects the dead fuels generally found in the ground; thus, it acts as a proxy for drought [52]. The normalized difference vegetation index (NDVI) was collected from MODIS/Terra with a spatial resolution of 375 m. Instead of using the mean NDVI, we calculated the median NDVI from 2016 to 2020 on a pixel-by-pixel basis to avoid bias introduced to the mean value by the greenness loss after wildfire [53]. NDVI displays the health status of vegetation. A reduction of the NDVI indicates dry grass or trees affecting the water stress and increasing the risk of fire [54]. The topographic wetness index (TWI) was generated from the Copernicus DEM and calculated with the relevant equation from Hong et al. [55]. TWI defines the aspect of steady-state soil wetness [56].
The anthropological factors were land use, distance to road, and distance to settlement. Land use and distance to settlement were acquired from the United States Geological Survey (USGS) national land cover database (NLCD-2019) with a spatial resolution of 30 m. Land use describes the landscape features of the study area. The different characteristics, such as the load and moisture content of the distinct land use types, can influence the ignition and spread of fires [11]. The distance to settlements and roads quantifies access to forest areas and, in many cases, plays an essential role in the ignition of wildfire incidents [15,41,50].

2.4. Damage Proxy Map (DPM)

To help identify damaged areas after wildfires, DPMs were produced by comparing pre- and co-event SAR images. The method depends on the decrease in the coherence of the signal scattering between SAR images taken prior and following the event to determine irregular changes in ground surface properties [57]. Coherence assesses the alteration in backscatter signal as a proxy for the ground-surface property alteration [58]. A high coherence implies little or no alteration to the ground surface reflecting the SAR radiation. Alterations can be caused by damage to the ground surface itself or damage to structures by wildfires.
Interferometric coherence is a quantification of resemblance between two radar signals that have been employed to assess the quality of the InSAR product. This statistical quantity is computed as [19]:
γ = c 1 c 2 * c 1 c 1 * c 1 c 2 * , 0 γ 1
where c1 and c2 represent the values of complex pixels of two SAR scenes and * indicates the complex conjugate. The resulting coherence ranges from 1 (coherent) to 0 (incoherent). If the two images’ observations are identical, then the coherence is equal to 1 due to stable objects in the scenes, such as buildings. However, the DPM method is not limited to building damage. As long as the predamage (reference) InSAR pair maintains reasonable coherence, this approach can be applied to any anthropogenic or natural damage that causes serious surface property changes, such as those caused by a wildfire. Major damage after a wildfire significantly intensifies the interferometric phase dissimilarity. This alteration appears as a decrease in coherence or decorrelation.
The process started with image co-registration with subpixel accuracy utilized to match scenes to one another. Three radar images were then employed to generate two pairs of interferometric coherence maps to produce a DPM: one pair prior to the damage (pre-event) and the other pair spanning the damage (co-event). The pre-event coherence pair delineated changes unrelated to the event and was expected to be the background value. In ArcGIS 10.4 software, we obtained the coherence difference (COD) by subtracting γpre-event and γco-event. The results were transformed into a georeferenced coordinate system employing the Copernicus DEM, resulting in pixels approximately 30 m across. The threshold for sufficient coherence loss was adjusted by comparing observed COD with reported damage and the area where there was no damage [58,59,60]. Here, the coherence loss threshold for DPM was chosen by considering wildfire perimeters from CAL FIRE. The fire perimeter was obtained using various collection methods, such as GPS ground, GPS air, infrared, photo interpretation, hand-drawn, and mixed collection methods.

2.5. Spatial Correlation Analysis

Variable selection is particularly important in the prediction of wildfire susceptibility. The high dimensionality of the training dataset may complicate the prediction process and decrease the prediction accuracy. In this study, the information gain ratio (IGR) and multicollinearity analysis were selected to evaluate the wildfire conditioning factors.
In addition, Pearson correlation analysis was applied to identify linear correlation relationships between pairs of variables. The Pearson correlation coefficient allows measurement of the association between variables. When a correlation is present, a change in the magnitude of one variable is associated with a change in the magnitude of the other variable, either in the same direction (positive coefficient) or in the opposite direction (negative coefficient). This coefficient is scaled and takes values between −1 and 1, where 0 is equivalent to the case in which no correlation exists [61]. A correlation coefficient greater than or equal to 0.7 is considered a correlation indicator that can lead to distortion of the modeling process and affect future predictions [62].
For the IGR method, the factor with a higher IGR value indicates a stronger predictive ability of the model. However, factors with IGR values equal to or less than 0 indicate a “null” contribution to the forest fire susceptibility model and should be excluded from further analysis [34].
The tolerance (TOLs) and variance inflation factors (VIFs) were computed to discover whether multicollinearity existed [63]. Multicollinearity occurs between variables if they have a close linear association in a regression model, which may decrease the model’s performance [64,65]. The TOL and VIF were calculated using the following equations:
T O L = 1 R j 2
V I F = 1 1 R j 2
where R j 2 denotes the regression value of j on different variables in a dataset. Thus, multicollinearity issues generally occur if the TOL score is smaller than 0.10 and the VIF score is greater than 10 [66,67].
Using the frequency ratio (FR) method, spatial relationships among the distribution of wildfire damage locations (71,371 points) and class of wildfire conditioning factors were analyzed. After classifying the numeric factors into five classes, tabulate area tools in ArcGIS 10.4 were used to handle dependent data (wildfire location) with classified factors for obtaining tabulate area cell counts indicating the number of wildfire occurrences in each class. Therefore, the FR value can be acquired by calculating the ratio between the percentage of pixels wildfire occurrence/event in each class with the percentage of pixels area of each class. If a specific class of the factor obtained a FR value greater than 1, then the class has a high probability of wildfire occurrence, high degree of damage by wildfire, and will have a substantial impact on wildfire susceptibility [68]. The FR value of each class factor was computed using the following equation:
F R = %   P i x e l   o f   w i l d f i r e   o c c u r r e n c e %   P i x e l   c l a s s   o f   w i l d f i r e   f a c t o r s
The FR value of each class of wildfire conditioning factors was utilized to produce wildfire susceptibility maps using deep learning algorithms based on CNN, CNN-GWO, and CNN-ICA.

2.6. Convolutional Neural Network (CNN)

Classification and prediction using CNNs have been increasingly applied recently in numerous disciplines as well as earth sciences [69,70,71]. CNN is one of the deep learning approaches differentiated from a traditional neural network by using multiple layers such as pooling, mutual weights, and local connections [72]. Generally, the CNN architecture is constructed of convolutional layers (CLs), activation process (AP), pooling layers (PLs), and fully connected layers (FCLs). The CLs continuously identify the connection between each feature class by extracting new features from the input images [73]. Next is AP, which chooses a rectified linear unit due to the unit having the ability to handle non-linear results from CLs and enhance the non-linear properties of the neural network [35,74]. PLs allow for stable conversion to prevent overfitting and upgrade computational performance by decreasing the number of feature class structures through the down-sampling approach [75]. Last are FCLs, which combine data in the neural network and generate the final output [73,76].

2.7. Metaheuristic Optimization Algorithms

2.7.1. Grey Wolf Optimization (GWO)

A metaheuristic algorithm was developed by mimicking the hierarchy of command and hunting strategy of grey wolves (Canis lupus) in nature [77]. The hierarchical structure of the grey wolf pack is based on four levels: alpha, beta, delta, and omega. Alpha α is the leader that controls all activities of the pack, such as resting, migration, feeding, and hunting. Beta β assists the leader in making judgements and establishing the order of the pack. Deltas ( δ ) act as hunters, sentinels, and watchers. Omegas ( ω ) are the weakest relations and care for the young; internal fights and other issues in the pack are noticed without them [78]. There are three main steps in hunting: seeking prey, surrounding the prey, and striking the prey [36]. The alpha α is considered the highest-fitted solution, followed by β ,     δ , and ω , based on Equations (5)–(8) [77]:
D α = C 1   · X α X ,   D β = C 2   · X β X ,   D δ = C 3   · X δ X ,   C = 2 × r 2
X 1 = X α A 1 × D α ,   X 2 = X β A 2 × D β ,   X 3 = X δ A 3 × D δ ,
A = 2 × a × r 1 a
X t + 1 = X 1 + X 2 + X 3 3
where for the iteration time t, D is the vector that suggests a new position and X denotes the position of the wolf; A and C are coefficient vectors and components of a , which are linearly diminished between 0 and 2 in each iteration; and r 1 and r 2 are random vectors generated for the range [0, 1]. The hunting process is completed when A takes values between −1 and 1, when an attack occurs [78].

2.7.2. Imperialist Competitive Algorithms (ICA)

ICA is one of the metaheuristic optimization algorithms imitating imperialistic competition and addresses various optimization problems [79]. ICA begins with a population of randomized solutions, where each solution is considered as a country. The population is separated into two categories: imperialist and colonies. Each imperialist initially controls a group of colonies and then competes to attract and take possession of others through its authority and influence, forming the core of the ICA algorithm [80]. The positions of imperialist countries and colonies can be exchanged if the colonies become more powerful. If an empire cannot increase its strength through interactions with other countries, it will be gradually eliminated during the competition. Thus, fragile empires lose their colonies, while substantial empires occupy more colonies, further increasing their strength. Ultimately, imperialistic competition leads to the elimination of all but one empire, which gains control of all other countries [81]. The algorithm terminates after the maximum number of iterations or when a predefined number of empires remain.

2.8. Accuracy Assessment

Evaluation is an essential step in assessing the accuracy of predictions of a model to support the scientific validity of a study [82]. This study used the area under the receiver operating characteristic (ROC) curve analysis for model assessment, which is a common way of evaluating wildfire probability models [11,15,59,83,84]. The AUC represents the performance, evaluation, utilization, and compression of model predictions [85,86,87]. This measure of accuracy ranges between 0.5 and 1 (perfect forecasting), with values near 1 indicating excellent performance and values near 0.5 denoting very poor prediction accuracy [69]. We applied this method using the testing dataset from wildfire inventory pixels (see Section 3.1) which were not used to train the applied deep learning approaches. A ROC curve is a plot of specificity (i.e., false positive on the x axis) versus sensitivity (i.e., true positive on the y axis). In wildfire modeling, sensitivity and specificity refer to the proportion of correctly predicted fire pixels and the proportion of correctly predicted non-fire pixels, respectively.
This study also used the mean square error (MSE) and root mean squared error (RMSE) approach as a second evaluation metric and as a cost function for optimizing CNN parameters as follow:
M S E = 1 n i = 1 n y i t i 2
R M S E = 1 n i = 1 n y i t i 2  
where y and t are the predicted and actual wildfire inventory values, respectively, and n denotes the number of samples. Metaheuristic algorithms were used to minimize the MSE and RMSE values through optimization of model hyperparameters. Smaller values of MSE and RMSE indicate a better performance of the models and the effectiveness of metaheuristic optimization. In addition, a wildfire susceptibility map with higher performance will be compared with recent wildfires in 2021 to evaluate the model’s predictive performance.

3. Results

3.1. Wildfire Damage Inventory Map

Most wildfires in Plumas National Forest occur during the summer or fall when most of the leaves have fallen. Thus, InSAR coherence measurements collected before a wildfire can maintain reasonable coherence, and the DPM technique can be applied to depict damaged locations after wildfires in Plumas National Forest.
DPMs were generated depicting areas in Plumas National Forest that were likely damaged by wildfires. Figure 5 reveals the DPM results after wildfires in Plumas National Forest, including the Camp Fire, Walker Fire, North Complex Fire, and Sheep Fire. Areas that experience an increase in COD are indicated by pixels with colors that stretch from yellow to red. Enlarged DPM pixels opaqueness represents significant ground or building alteration or possible damage due to wildfire. Regions where decorrelation did not alter significantly over wildfires during the period are set to be transparent, indicating no destruction. The DPMs were geocoded to the Copernicus DEM with a corresponding spatial resolution of approximately 30 m.
This study compared the distribution of coherence loss areas with the fire perimeter from CAL FIRE which showed similar damage areas after wildfires throughout Plumas National Forest from 2016 to 2020. The scattered pixels were removed using a fire perimeter to generate precise wildfire-affected areas. The combination of both sources resulted in a more accurate and dependable inventory dataset depicting the location of damage after wildfires. The total number of wildfires mapped was 11, which translated to 71,371 pixels of 30 m wildfire cell size as wildfire occurrence data. These data were then randomly divided into a 50% training and 50% testing dataset. In the context of deep learning, the creation of a wildfire susceptibility map required data of areas without wildfires. In this study, the same number of non-wildfire location data (71,371 points) were sampled through random selection by comparing the area with no prior wildfire and the very low probability class identified with the frequency ratio method. The analysis of zero data based on these results was an efficient way to aid the interpretation of the area [35]. The non-wildfire data were also divided into a training (50%) and testing (50%) dataset. The 50% training datasets for wildfire-occurrence and non-wildfire-occurrence location were then merged to generate wildfire susceptibility maps, with the remaining 50% of testing data from both datasets merged for evaluation of the performance [59,84,85].

3.2. Relationship between Damage Area and Related Factors

Figure 6 shows the results from the Pearson correlation analysis for all the conditioning factors. According to the pairwise correlation matrix, none of the correlation values among each conditioning factor were within the safe threshold. The highest correlation value was computed between solar and maximum temperature (0.63). The results suggested that there is no need to eliminate any conditioning factors and all factors would not cause distortion in the modeling process.
Based on the evaluation of the IGR analysis, all 16 conditioning factors considered in this study had predictive capabilities for wildfire modeling and contributed to the occurrence and extent of wildfires (IGR > 0), as shown in Table 4. Table 4 shows that the factor with the highest IGR value was land use (IGR = 0.39); thus, land use was the most effective wildfire conditioning factor in this study, followed by drought index (0.25) and maximum temperature (0.18). Furthermore, windspeed was found to be the least important wildfire conditioning factor with an IGR value of 0.03.
An assessment of multicollinearity was conducted to investigate the correlation between wildfire conditioning factors, as shown in Table 4. All wildfire conditioning factors had VIF scores lower than 10 and tolerance scores >0.1. The ranges of VIF and TOL were within the permissible threshold; therefore, there were no multicollinearities observed between the wildfire conditioning factors, thereby avoiding the risk of deteriorating model performance. The maximum VIF score was 5.91 for altitude and the lowest tolerance score was 0.18 for maximum temperature. The range of TOL was 0.17 to 0.99. Given that the Pearson correlation analysis, IGR, and multicollinearity values were all within their critical value, all of the conditioning factors can be used for creating a wildfire susceptibility map.
FR values can provide information about the connection between the location of damage caused by wildfires and classes of conditioning factors. The results of the FR calculation in this study are shown in Table 5. The slope aspect showed that wildfire occurrence is concentrated in the northeast, southeast, south, southwest, west, and flat portions of Plumas National Forest. The fact that south-facing aspects receive more solar radiation, which increases fuel temperature, and low fuel moisture in the North Hemisphere led to wildfire occurrence [50,88]. The areas of 13–272 and 272–1216 m in altitude had high FR values of 1.03 and 1.54, respectively, indicating that low-altitude areas in the study area are more susceptible to wildfire occurrence and implying vegetation is more burnable due to high temperature and dry weather during summer [10,50]. The high degree of slope causes wildfires to spread more quickly up the steep areas, as shown by the two classes of 7.22–12.89 and 20.37–65.75 having FR values of 1.06 and 1.44, respectively. Higher distance to road indicates a higher probability of wildfire occurrence, with the three highest classes of 152–301, 301–595, and 595–4260 m having FR values of 1.19, 1.14, and 1.27, respectively. The further away the road is, the more difficult it is for firefighters to access and extinguish the wildfire area [15]. Portions of landscape with maximum temperature > 65.67 °C, a higher drought index of −0.09, and precipitation between 25–43.35 mm were predestined places for wildfires, encompassing 52.88% (Tmax > 65.57 °C), 53.20% (drought > −0.09–0.89), and 47.44% (25–43.35 mm precipitation). High NDVI classes, namely 0.46–0.61, 0.61–0.76, and 0.76–1, had FR values of 1.04, 1.10, and 1.14, respectively, implying wildfire occurrence in high greenness and the moisture content of vegetation. The reverse trends were found for distance to settlement, with the low classes showing high FR values of 1.28 and 1.30 for distance to settlement of 287–1148 and 1148–2679 m, respectively. A closer distance to stream, namely 0–997, 997–2160, and 2160–3573 m, had high FR values of 1.17, 1.19, and 1.05, respectively. Streams cater for entertaining human interests such as tourist camps. The closer the distance to settlements and streams relating to high levels of human activity, the greater the risk of wildfire ignition [2,14]. For land use factors, herbaceous, evergreen, and mixed forest classes, with the highest FR values, experienced the highest wildfire occurrence and served as a fuel for wildfires [10,14]. Based on their FR values, plan curvature with concave and convex were associated with a high occurrence of wildfires, which is in agreement with other findings that suggest the probability of fire occurrence may be lower on flat terrain and higher on concave slopes [24]. Topographical and anthropological features also swamped the effect of soil moisture, TWI, and windspeed, all of which did not show any pattern or specific class that is strongly associated with the presence or absence of wildfires. The wildfire occurrence for TWI and windspeed features tended to be high in most classes, although our windspeed variable was intended to reflect the outcome of providing fresh oxygen and a greater potential for drying potential fuels and propelling fire across land at a faster rate.

3.3. Wildfire Susceptibility Map

The construction of the database has been performed, including: (i) the dependent variables (i.e., damage location by wildfire derived by the DPM method and Sentinel-1 product) and (ii) the independent variables (based on topographical, meteorological, environmental, and anthropological). The creation of wildfire susceptibility maps utilizing the training dataset were compiled employing the damage inventory database from the DPM of the wildfires from 2016 to 2020 and applying deep learning with metaheuristic optimization approaches, as discussed above. The wildfire susceptibility maps depict an assessment of the likelihood of experiencing a wildfire in a study area built upon the 16 conditioning factors considered. Figure 7 shows the wildfire susceptibility maps from the CNN, CNN-GWO, and CNN-ICA models. Using the quantile method, we classified the regions of the wildfire susceptibility maps into five predicted classes: very low, low, moderate, high, and very high [89,90].
The distribution of pixels in each wildfire susceptibility map was also analyzed and is shown in Figure 8. Overall, most of the study area is prone to wildfire occurrences, especially in the western part of Plumas National Forest and expanded to the center of the forest, and makes the study area a wildfire hot-spot region in California. Generally, about 40% of the study area located in areas with flat curvature and low altitude has a low to very low wildfire susceptibility, 20% of the study area has a moderate wildfire susceptibility, and 40% of the study area located in a high degree of slope has a high to very high susceptibility to future wildfire occurrence. In the CNN-ICA model, the moderate class was more distributed throughout the study area than in the other models. The distribution of pixels of standalone CNN and hybrid CNN-GWO showed similar results compared with those of CNN-ICA. CNN-ICA showed the smallest percentage of the very high class and the greatest percentage of the low wildfire susceptibility class.

3.4. Model Evaluation

A critical step in generating wildfire susceptibility models is the evaluation of the predicted models to ensure that they provide accurate and reliable wildfire susceptibility assessments. The performance of the models was evaluated and compared to assess the reliability of the wildfire susceptibility maps obtained with each algorithm. Figure 9 has two parts: errors versus number of samples and frequency versus errors. The error part of Figure 9 specifies the values of MSE and RMSE; the frequency versus errors depicts the values of the error mean and standard deviation (StD).
The results revealed that in the testing phase using the standalone CNN model, the values of MSE, RMSE, error mean, and error StD are 0.609, 0.780, −0.076, and 0.778, respectively. In the CNN-GWO model, the values of MSE, RMSE, error mean, and error StD are 0.334, 0.112, −0.014, and 0.335, respectively. Moreover, in the CNN-ICA model, the values of MSE, RMSE, error mean, and error StD are 0.129, 0.351, 0.071, and 0.344, respectively. In addition, the range of error in CNN alone was much broader (−2 < error < 2) than that in CNN-GWO and CNN-ICA (−1 < error < 1).
Figure 10 shows the ROC curves, with AUC values of 0.934, 0.950, and 0.974 for CNN (red line), CNN-GWO (blue line), and CNN-ICA (green line), respectively. The results of AUC were in agreement with the results of the model evaluation using MSE, RMSE, error mean, and error StD values in the testing phase. Therefore, the predictive ability of CNN-GWO was better than that of CNN and CNN-ICA, as indicated by lower values of MSE and RMSE and a higher value of AUC. This finding is in accordance with the results of other studies [35,91].

4. Discussion

Wildfire inventories are crucial for accurate wildfire susceptibility creation. The application of remote sensing SAR imagery is suitable for obtaining data inventories of locations damaged by wildfires due to their wide availability, independence from fieldwork, time and cost efficiency, frequent repeatability over time, cloud and sunlight independence, and high precision and sensitivity [37]. Combined with the fire perimeter generated from CAL FIRE, a more precise and reliable inventory dataset depicting the location of damage and the area of wildfires was obtained. In this study, the DPM technique was employed to determine the areas damaged by wildfires in Plumas National Forest to form inventory data utilized for training and testing the wildfire susceptibility models.
The ARIA team also created a DPM after the Camp Fire wildfire and found similar results [86]. Therefore, using SAR imagery, the DPM method can be an alternative tool to map damaged areas. However, scattered small-colored pixels were spotted outside the perimeter and may be false positives but do not necessarily mean no damage occurred. Our study revealed that the DPM technique is not limited to detecting damage to buildings after an earthquake [57,92], volcano [20], or typhoon [19]. The DPM method can be applied to any natural or anthropogenic damage that causes a substantial change in the surface property as long as a predamage (reference) InSAR pair preserves reasonable coherence. Tay et al. (2020) employed aerial imagery, obtained from the Japan Geospatial Information Authority (GSI), after Typhoon Hagibis. Nur et al. (2021) used a damage area report from the Korea Meteorological Agency (KMA) report of the Gyeongju and Pohang earthquakes.
Important conditioning factors for detecting wildfire-prone areas were identified and collected based on the literature and data availability. Spatial correlation analysis has been performed by calculating the VIF and TOL for multicollinearity analysis and the FR and IGR techniques to evaluate the relationship and influence of the factors. Similar results were acquired by [93], who indicated a higher VIF score for the altitude factor. This is due to the fluctuation of most factors, such as maximum temperature, drought index, and soil moisture, being consistent with the altitude. The output showed that no significant multicollinearity existed, with all factors used in this study having a considerable effect on wildfire susceptibility in Plumas National Forest. FR results clearly showed that the spatial relationship between each category of conditioning factors and wildfire location is not randomly distributed across Plumas National Forest and that the likelihood of wildfires is highly dependent on the characteristics of the landscapes. Although some variables were not shown to have a significant effect on wildfire occurrence when evaluated separately, as shown in other studies, we were able to identify fire-prone areas in the study area with the help of these conditioning factors. Areas with a high probability of wildfires were associated with land use, drought index, and maximum temperature. These were consistent with many studies that cite temperature, drought index, and precipitation as critical parameters affecting the relative likelihood of wildfire occurrence [14,29], presumably because fuel moisture content is largely a function of temperature, drought index, and precipitation [94]. This study also confirms previous results that found wind is the most unpredictable variable for forecasting wildfire occurrence [24].
Land cover describes the vegetation type, structure, amount, and continuity, with the fuel properties of vegetation and fuel breaks being factors that influence the characteristics of wildfires, such as their likelihood, spread, frequency, and severity, by providing different fuel abundance and settings for different time intervals. In addition, land uses such as urbanization, agriculture, and grazing and land covers such as barren land, open water, and wetland in Plumas National Forest may stop wildfires from spreading because there is insufficient continuous vegetation to carry them. The land use factors in the study area contributed to wildfires, especially in herbaceous areas with the highest FR value and evergreen forest areas with the highest wildfire occurrence. Herbaceous vegetation area is more susceptible to wildfires when the climate is drier, producing fuel, igniting fires, or spreading frequent surface fires [95]. In the forest area, tree leaf litter can also provide fuel when herbaceous groundcover fails to develop a wildfire. The combination of herbaceous vegetation and evergreen forest areas were the primary wildland classes that carried wildfires in the U.S. [96]. In recent decades, the western U.S. has warmed and the frequency and severity of heatwaves and droughts have increased [97]. Based on the National Oceanic and Atmospheric Administration (NOAA), temperatures in California have increased by approximately 3 degrees Fahrenheit (1.67 °C) since the start of the 20th century, with the 2015–2020 period being the warmest, with the highest number of extremely hot days [98]. High temperatures have caused the air to dry out. Fire seasons start earlier and end later each year. Meanwhile, snowpacks are shrinking, causing an earlier spring and an intense and longer dry season. These drier and warmer conditions are also exacerbating wildfires in the western U.S. [9].
In this study, a CNN deep learning algorithm and hybrid models, namely CNN-GWO and CNN-ICA, were used to determine the areas susceptible to wildfires. All proposed models showed good results with AUC values greater than 0.9 and yielded maps of wildfire susceptibility in the study area. The performance of CNN-GWO was the best, with AUC = 0.974 and RMSE = 0.334, followed by CNN-ICA (AUC = 0.950, RMSE = 0.351) and standalone CNN (AUC = 0.934, RMSE = 0.780). In general, the predictive ability of deep learning algorithms can vary depending on the model structure, input selection, dataset quantity and quality, and optimization of model parameters [35]. The hyperparameter adjustment of the CNN model using a metaheuristic algorithm affected the prediction performance of the model. The hybrid models obtained better results than the standalone models. Therefore, hybrid models are important techniques for improving the prediction capacity of basic classification to reduce bias and avoid the problem of overfitting. Furthermore, from the literature reviews, the GWO algorithm has a simple structure and the ability to converge quickly due to the continuous reduction of the search space and fewer tuning variables than other algorithms [77]. Therefore, GWO was both faster and had higher prediction accuracy.
Figure 11 shows the CNN-GWO models which performed better with wildfires that occurred in 2021, including the Beckwourth Complex Fire, Dixie Fire, Gunnison Fire, and Park Fire. The RMSE value for the wildfire susceptibility map from the CNN-GWO model with recent wildfires was 0.511, which was lower than the standalone CNN (0.678) and CNN-ICA (0.514).
Gunnison Fire and Park Fire burned 222 ha and 163 ha, respectively. Both wildfires were located in very high wildfire susceptibility areas. The Beckwourth Complex Fire burned 42,700 ha located outside wildfire-prone areas. This is because the model had difficulty predicting the location of wildfires where the area did not meet the required or similar conditions to high-risk areas. Additionally, the area has no history of wildfires. According to the Beckwourth Post-Fire BAER Assessment Report, the cause was lightning. The Dixie Fire burned 389,800 ha, killed one person, and damaged 1329 structures in July 2021. According to the Fire Information for Resource Management System (FIRMS) data, the initial fire was located in an area of very high wildfire susceptibility in the center of the study area (darker purple area in Figure 11). Therefore, our models managed to predict the occurrence of most wildfires in the study area.
This study has the limitation of depicting damaged areas in areas of low coherence before wildfires occurred. Although most wildfires occurred in summer and fall, some forest areas still had dense leaves showing incoherence. Other data should be combined with the DPM and fire perimeter from CAL FIRE to map more precise locations of areas damaged by wildfire. The normalized burn ratio (NBR) from Sentinel-2 imagery with 10 m resolution can reflect the burn severity level. Another limitation is the use of NDVI, solar radiance, drought index, and temperature data. NDVI was computed from MODIS images with 375 m resolution. The NDVI data were influenced by clouds and had certain spatial and temporal limitations. Future study should include the soil adjusted vegetation index (SAVI) and the normalized difference moisture index (NDMI) for more accurate analysis of the vegetation composition in Plumas National Forest [99]. The solar radiance, drought index, and soil moisture data were acquired from Terra Climate images with low resolution. Therefore, the actual climate change in the study area has not been well described, thus introducing uncertainty in the generation of wildfire susceptibility maps. Future study will integrate the capability of Sentinel-1 and Sentinel-2 for better and accurate wildfire mapping.
Along with wildfire susceptibility maps, emergency treatments and activities must be designed to decrease possible impacts on valued resources in the wildfire area, such as life and safety, property, critical natural resources, and cultural resources.

5. Conclusions

Wildfire susceptibility maps in Plumas National Forest were created for the first time by integrating SAR data using the DPM method with deep learning based on CNN with metaheuristic optimization algorithms. Identifying areas with very high and high susceptibility provides information that can be applied to risk management and the development of innovative responses to future wildfires. The Sentinel-1 SAR dataset and the DPM technique created a wildfire inventory of the study area with high accuracy and precision as well as low cost and reproducibility. CNN deep learning with metaheuristic optimization algorithms, namely GWO and ICA, were employed to model wildfire susceptibility in Plumas National Forest, California, USA. The results showed that CNN-GWO had a lower RMSE (0.334) and a higher area under the ROC curve value (0.974) than CNN-GWO and CNN alone. In conclusion, the results revealed the applicability of the CNN deep learning algorithm to wildfire susceptibility mapping and the effectiveness of the hybrid approach using metaheuristic optimization algorithms. Furthermore, the methodology is reproducible and applicable to other areas with different conditioning factors. Further study to compare the forecast accuracy of other deep learning methods and metaheuristic optimization algorithms is warranted.

Author Contributions

Conceptualization, C.-W.L.; methodology, C.-W.L. and A.S.N.; software, C.-W.L. and A.S.N.; evaluation, C.-W.L., Y.J.K. and A.S.N.; formal analysis, C.-W.L. and A.S.N.; investigation, C.-W.L., Y.J.K. and A.S.N.; resources, C.-W.L., Y.J.K. and A.S.N.; data curation, C.-W.L., Y.J.K. and A.S.N.; writing—original draft preparation, A.S.N.; writing—review and editing, Y.J.K. and C.-W.L.; visualization, C.-W.L. and A.S.N.; supervision, C.-W.L.; project administration, C.-W.L.; funding acquisition, C.-W.L. All authors have read and agreed to the published version of the manuscript.

Funding

This research was supported by a grant from the Korea Polar Research Institute (KOPRI, PE22900).

Data Availability Statement

Data supporting the reported results can be found at: https://search.asf.alaska.edu/ for Sentinel-1 data, (accessed on 8 July 2022); https://portal.opentopography.org/raster?opentopoID=OTSDEM.032021.4326.3 for Copernicus DEM data, (accessed on 8 July 2022); https://prism.oregonstate.edu/ for precipitation and temperature data, accessed on 8 July 2022; https://www.nrel.gov/gis/solar.html for solar radiance data, accessed on 8 July 2022; https://globalwindatlas.info/ for windspeed data, accessed on 8 July 2022; https://gis.data.ca.gov/ for stream data, accessed on 8 July 2022; https://www.climatologylab.org/terraclimate.html for drought and soil moisture data, accessed on 8 July 2022; https://search.earthdata.nasa.gov/ for NDVI MODIS data, accessed on 8 July 2022; and https://data.usgs.gov/datacatalog/data/USGS:60cb3da7d34e86b938a30cb9 for land use data, accessed on 8 July 2022.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study, in the collection, analyses, or interpretation of data, in the writing of the manuscript, or in the decision to publish the results.

References

  1. Porter, T.W.; Crowfoot, W.; Newsom, G. 2020 Wildfire Activity Statistics; California Department of Forestry and Fire Protection: Sacramento, CA, USA, 2019; p. 46. [Google Scholar]
  2. Li, S.; Banerjee, T. Spatial and Temporal Pattern of Wildfires in California from 2000 to 2019. Sci. Rep. 2021, 11, 8779. [Google Scholar] [CrossRef] [PubMed]
  3. Nauman, B.C. Variability in California’s Fire Activity during the Holocene, across Space and Time; University of California: Los Angeles, CA, USA, 2020. [Google Scholar]
  4. Radeloff, V.C.; Helmers, D.P.; Kramer, H.A.; Mockrin, M.H.; Alexandre, P.M.; Bar-Massada, A.; Butsic, V.; Hawbaker, T.J.; Martinuzzi, S.; Syphard, A.D.; et al. Rapid Growth of the US Wildland-Urban Interface Raises Wildfire Risk. Proc. Natl. Acad. Sci. USA 2018, 115, 3314–3319. [Google Scholar] [CrossRef] [PubMed]
  5. CAL FIRE Stats and Events. Available online: https://www.fire.ca.gov/stats-events/ (accessed on 10 January 2022).
  6. Wong, S.D.; Broader, J.C.; Shaseen, S.A. Review of California Wildfire Evacuations from 2017 to 2019; University of California: Berkeley, CA, USA, 2020. [Google Scholar]
  7. Luo, N.; Weng, W.; Xu, X.; Hong, T.; Fu, M.; Sun, K. Assessment of Occupant-Behavior-Based Indoor Air Quality and Its Impacts on Human Exposure Risk: A Case Study Based on the Wildfires in Northern California. Sci. Total Environ. 2019, 686, 1251–1261. [Google Scholar] [CrossRef] [PubMed]
  8. Aguilera, R.; Corringham, T.; Gershunov, A.; Benmarhnia, T. Wildfire Smoke Impacts Respiratory Health More than Fine Particles from Other Sources: Observational Evidence from Southern California. Nat. Commun. 2021, 12, 1493. [Google Scholar] [CrossRef]
  9. Buis, A.; The Climate Connections of a Record Fire Year in the U.S. West. Available online: https://climate.nasa.gov/ask-nasa-climate/3066/the-climate-connections-of-a-record-fire-year-in-the-us-west/ (accessed on 4 April 2022).
  10. Tonini, M.; D’andrea, M.; Biondi, G.; Esposti, S.D.; Trucchia, A.; Fiorucci, P. A Machine Learning-Based Approach for Wildfire Susceptibility Mapping. The Case Study of the Liguria Region in Italy. Geosciences 2020, 10, 105. [Google Scholar] [CrossRef]
  11. Bustillo Sánchez, M.; Tonini, M.; Mapelli, A.; Fiorucci, P. Spatial Assessment of Wildfires Susceptibility in Santa Cruz (Bolivia) Using Random Forest. Geosciences 2021, 11, 224. [Google Scholar] [CrossRef]
  12. Jaafari, A.; Zenner, E.K.; Pham, B.T. Wildfire Spatial Pattern Analysis in the Zagros Mountains, Iran: A Comparative Study of Decision Tree Based Classifiers. Ecol. Inform. 2018, 43, 200–211. [Google Scholar] [CrossRef]
  13. Saim, A.A.; Aly, M.H. Machine Learning for Modeling Wildfire Susceptibility at the State Level: An Example from Arkansas, USA. Geographies 2022, 2, 31–47. [Google Scholar] [CrossRef]
  14. Ghorbanzadeh, O.; Kamran, K.V.; Blaschke, T.; Aryal, J.; Naboureh, A.; Einali, J.; Bian, J. Spatial Prediction of Wildfire Susceptibility Using Field Survey Gps Data and Machine Learning Approaches. Fire 2019, 2, 43. [Google Scholar] [CrossRef]
  15. Piralilou, S.T.; Einali, G.; Ghorbanzadeh, O.; Nachappa, T.G.; Gholamnia, K.; Blaschke, T.; Ghamisi, P. A Google Earth Engine Approach for Wildfire Susceptibility Prediction Fusion with Remote Sensing Data of Different Spatial Resolutions. Remote Sens. 2022, 14, 672. [Google Scholar] [CrossRef]
  16. Hong, H.; Jaafari, A.; Zenner, E.K. Predicting Spatial Patterns of Wildfire Susceptibility in the Huichang County, China: An Integrated Model to Analysis of Landscape Indicators. Ecol. Indic. 2019, 101, 878–891. [Google Scholar] [CrossRef]
  17. Zhang, P.; Nascetti, A.; Ban, Y.; Gong, M. An Implicit Radar Convolutional Burn Index for Burnt Area Mapping with Sentinel-1 C-Band SAR Data. ISPRS J. Photogramm. Remote Sens. 2019, 158, 50–62. [Google Scholar] [CrossRef]
  18. Nur, A.S.; Lee, C. Damage Proxy Map (DPM) of the 2016 Gyeongju and 2017 Pohang Earthquakes Using Sentinel-1 Imagery. Korean J. Remote Sens. 2021, 37, 13–22. [Google Scholar] [CrossRef]
  19. Tay, C.W.J.; Yun, S.-H.; Chin, S.T.; Bhardwaj, A.; Jung, J.; Hill, E.M. Rapid Flood and Damage Mapping Using Synthetic Aperture Radar in Response to Typhoon Hagibis, Japan. Sci. Data 2020, 7, 100–108. [Google Scholar] [CrossRef] [PubMed]
  20. Biass, S.; Jenkins, S.; Lallemant, D.; Lim, T.N.; Williams, G.; Yun, S.-H. Remote Sensing of Volcanic Impacts. In Forecasting and Planning for Volcanic Hazards, Risks, and Disasters; Elsevier: Amsterdam, The Netherlands, 2021; pp. 473–491. [Google Scholar]
  21. 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] [CrossRef]
  22. Jaafari, A.; Gholami, D.M. Wildfire Hazard Mapping Using an Ensemble Method of Frequency Ratio with Shannon’s Entropy. Iran. J. For. Poplar Res. 2017, 25, 232–242. [Google Scholar]
  23. Nami, M.H.; Jaafari, A.; Fallah, M.; Nabiuni, S. Spatial Prediction of Wildfire Probability in the Hyrcanian Ecoregion Using Evidential Belief Function Model and GIS. Int. J. Environ. Sci. Technol. 2018, 15, 373–384. [Google Scholar] [CrossRef]
  24. Jaafari, A.; Gholami, D.M.; Zenner, E.K. A Bayesian Modeling of Wildfire Probability in the Zagros Mountains, Iran. Ecol. Inform. 2017, 39, 32–44. [Google Scholar] [CrossRef]
  25. Nikhil, S.; Danumah, J.H.; Saha, S.; Prasad, M.K.; Rajaneesh, A.; Mammen, P.C.; Ajin, R.S.; Kuriakose, S.L. Application of GIS and AHP Method in Forest Fire Risk Zone Mapping: A Study of the Parambikulam Tiger Reserve, Kerala, India. J. Geovis. Spat. Anal. 2021, 5, 14. [Google Scholar] [CrossRef]
  26. Busico, G.; Giuditta, E.; Kazakis, N.; Colombani, N. A Hybrid GIS and AHP Approach for Modelling Actual and Future Forest Fire Risk Under Climate Change Accounting Water Resources Attenuation Role. Sustainability 2019, 11, 7166. [Google Scholar] [CrossRef]
  27. Rodrigues, M.; Jiménez-Ruano, A.; Peña-Angulo, D.; de la Riva, J. A Comprehensive Spatial-Temporal Analysis of Driving Factors of Human-Caused Wildfires in Spain Using Geographically Weighted Logistic Regression. J. Environ. Manag. 2018, 225, 177–192. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  28. Jafari Goldarag, Y.; Mohammadzadeh, A.; Ardakani, A.S. Fire Risk Assessment Using Neural Network and Logistic Regression. J. Indian Soc. Remote Sens. 2016, 44, 885–894. [Google Scholar] [CrossRef]
  29. Gholamnia, K.; Nachappa, T.G.; Ghorbanzadeh, O.; Blaschke, T. Comparisons of Diverse Machine Learning Approaches for Wildfire Susceptibility Mapping. Symmetry 2020, 12, 604. [Google Scholar] [CrossRef]
  30. Al-Fugara, A.; Mabdeh, A.N.; Ahmadlou, M.; Pourghasemi, H.R.; Al-Adamat, R.; Pradhan, B.; Al-Shabeeb, A.R. Wildland Fire Susceptibility Mapping Using Support Vector Regression and Adaptive Neuro-Fuzzy Inference System-Based Whale Optimization Algorithm and Simulated Annealing. ISPRS Int. J. Geoinf. 2021, 10, 382. [Google Scholar] [CrossRef]
  31. 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]
  32. Nguyen, H.D. Hybrid Models Based on Deep Learning Neural Network and Optimization Algorithms for the Spatial Prediction of Tropical Forest Fire Susceptibility in Nghe An Province, Vietnam. Geocarto Int. 2022, 2022, 1–25. [Google Scholar] [CrossRef]
  33. Bjånes, A.; De La Fuente, R.; Mena, P. A Deep Learning Ensemble Model for Wildfire Susceptibility Mapping. Ecol. Inform. 2021, 65, 101397. [Google Scholar] [CrossRef]
  34. 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]
  35. Hakim, W.L.; Rezaie, F.; Nur, A.S.; Panahi, M.; Khosravi, K.; Lee, C.-W.; Lee, S. Convolutional Neural Network (CNN) with Metaheuristic Optimization Algorithms for Landslide Susceptibility Mapping in Icheon, South Korea. J. Environ. Manag. 2022, 305, 114367. [Google Scholar] [CrossRef]
  36. Balogun, A.L.; Rezaie, F.; Pham, Q.B.; Gigović, L.; Drobnjak, S.; Aina, Y.A.; Panahi, M.; Yekeen, S.T.; Lee, S. Spatial Prediction of Landslide Susceptibility in Western Serbia Using Hybrid Support Vector Regression (SVR) with GWO, BAT and COA Algorithms. Geosci. Front. 2021, 12, 101104. [Google Scholar] [CrossRef]
  37. Ranjgar, B.; Razavi-Termeh, S.V.; Foroughnia, F.; Sadeghi-Niaraki, A.; Perissin, D. Land Subsidence Susceptibility Mapping Using Persistent Scatterer SAR Interferometry Technique and Optimized Hybrid Machine Learning Algorithms. Remote Sens. 2021, 13, 1326. [Google Scholar] [CrossRef]
  38. Tien Bui, D.; Khosravi, K.; Li, S.; Shahabi, H.; Panahi, M.; Singh, V.; Chapi, K.; Shirzadi, A.; Panahi, S.; Chen, W.; et al. New Hybrids of ANFIS with Several Optimization Algorithms for Flood Susceptibility Modeling. Water 2018, 10, 1210. [Google Scholar] [CrossRef]
  39. Dong, L.; Leung, L.R.; Qian, Y.; Zou, Y.; Song, F.; Chen, X. Meteorological Environments Associated With California Wildfires and Their Potential Roles in Wildfire Changes During 1984–2017. J. Geophys. Res. Atmos. 2021, 126, e2020JD033180. [Google Scholar] [CrossRef]
  40. Malik, A.; Rao, M.R.; Puppala, N.; Koouri, P.; Thota, V.A.K.; Liu, Q.; Chiao, S.; Gao, J. Data-Driven Wildfire Risk Prediction in Northern California. Atmosphere 2021, 12, 109. [Google Scholar] [CrossRef]
  41. Baltar, M.; Keeley, J.E.; Schoenberg, F.P. County-Level Analysis of the Impact of Temperature and Population Increases on California Wildfire Data. Environmetrics 2014, 25, 397–405. [Google Scholar] [CrossRef]
  42. Jin, Y.; Goulden, M.L.; Faivre, N.; Veraverbeke, S.; Sun, F.; Hall, A.; Hand, M.S.; Hook, S.; Randerson, J.T. Identification of Two Distinct Fire Regimes in Southern California: Implications for Economic Impact and Future Change. Environ. Res. Lett. 2015, 10, 094005. [Google Scholar] [CrossRef]
  43. Sulova, A.; Arsanjani, J.J. Exploratory Analysis of Driving Force of Wildfires in Australia: An Application of Machine Learning within Google Earth Engine. Remote Sens. 2020, 13, 10. [Google Scholar] [CrossRef]
  44. Sazib, N.; Bolten, J.D.; Mladenova, I.E. Leveraging NASA Soil Moisture Active Passive for Assessing Fire Susceptibility and Potential Impacts over Australia and California. IEEE J. Sel. Top. Appl. Earth Obs. Remote Sens. 2022, 15, 779–787. [Google Scholar] [CrossRef]
  45. Jaafari, A.; Zenner, E.K.; Panahi, M.; Shahabi, H. Hybrid Artificial Intelligence Models Based on a Neuro-Fuzzy System and Metaheuristic Optimization Algorithms for Spatial Prediction of Wildfire Probability. Agric. For. Meteorol. 2019, 266–267, 198–207. [Google Scholar] [CrossRef]
  46. Sengupta, M.; Xie, Y.; Lopez, A.; Habte, A.; Maclaurin, G.; Shelby, J. The National Solar Radiation Data Base (NSRDB). Renew. Sustain. Energy Rev. 2018, 89, 51–60. [Google Scholar] [CrossRef]
  47. 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]
  48. Eskandari, S.; Pourghasemi, H.R.; Tiefenbacher, J.P. Fire-Susceptibility Mapping in the Natural Areas of Iran Using New and Ensemble Data-Mining Models. Environ. Sci. Pollut. Res. 2021, 28, 47395–47406. [Google Scholar] [CrossRef]
  49. Shakesby, R.A. Post-Wildfire Soil Erosion in the Mediterranean: Review and Future Research Directions. Earth Sci. Rev. 2011, 105, 71–100. [Google Scholar] [CrossRef]
  50. 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]
  51. Abatzoglou, J.T.; Dobrowski, S.Z.; Parks, S.A.; Hegewisch, K.C. TerraClimate, a High-Resolution Global Dataset of Monthly Climate and Climatic Water Balance from 1958–2015. Sci. Data 2018, 5, 170191. [Google Scholar] [CrossRef] [PubMed]
  52. Chuvieco, E.; Aguado, I.; Dimitrakopoulos, A.P. Conversion of Fuel Moisture Content Values to Ignition Potential for Integrated Fire Danger Assessment. Can. J. For. Res. 2004, 34, 2284–2293. [Google Scholar] [CrossRef]
  53. Zhang, Y.; Lim, S.; Sharples, J.J. Wildfire Occurrence Patterns in Ecoregions of New South Wales and Australian Capital Territory, Australia. Nat. Hazards 2017, 87, 415–435. [Google Scholar] [CrossRef]
  54. Tien Bui, D.; 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] [CrossRef]
  55. Hong, H.; Pourghasemi, H.R.; Pourtaghi, Z.S. Landslide Susceptibility Assessment in Lianhua County (China): A Comparison between a Random Forest Data Mining Technique and Bivariate and Multivariate Statistical Models. Geomorphology 2016, 259, 105–118. [Google Scholar] [CrossRef]
  56. Porensky, L.M.; Derner, J.D.; Pellatz, D.W. Plant Community Responses to Historical Wildfire in a Shrubland–Grassland Ecotone Reveal Hybrid Disturbance Response. Ecosphere 2018, 9, e02363. [Google Scholar] [CrossRef]
  57. Hough, S.E.; Yun, S.-H.; Jung, J.; Thompson, E.; Parker, G.A.; Stephenson, O. Near-Field Ground Motions and Shaking from the 2019 Mw 7.1 Ridgecrest, California, Mainshock: Insights from Instrumental, Macroseismic Intensity, and Remote-Sensing Data. Bull. Seismol. Soc. Am. 2020, 110, 1506–1516. [Google Scholar] [CrossRef]
  58. Agapiou, A. Damage Proxy Map of the Beirut Explosion on 4th of August 2020 as Observed from the Copernicus Sensors. Sensors 2020, 20, 6382. [Google Scholar] [CrossRef]
  59. Han, J.; Nur, A.; Syifa, M.; Ha, M.; Lee, C.-W.; Lee, K.-Y. Improvement of Earthquake Risk Awareness and Seismic Literacy of Korean Citizens through Earthquake Vulnerability Map from the 2017 Pohang Earthquake, South Korea. Remote Sens. 2021, 13, 1365. [Google Scholar] [CrossRef]
  60. Yun, S.-H.; Hudnut, K.; Owen, S.; Webb, F.; Simons, M.; Sacco, P.; Gurrola, E.; Manipon, G.; Liang, C.; Fielding, E.; et al. Rapid Damage Mapping for the 2015 Mw 7.8 r Earthquake Using Synthetic Aperture Radar Data from COSMO–SkyMed and ALOS-2 Satellites. Seismol. Res. Lett. 2015, 86, 1549–1556. [Google Scholar] [CrossRef]
  61. Schober, P.; Schwarte, L.A. Correlation Coefficients: Appropriate Use and Interpretation. Anesth. Analg. 2018, 126, 1763–1768. [Google Scholar] [CrossRef] [PubMed]
  62. Dormann, C.F.; Elith, J.; Bacher, S.; Buchmann, C.; Carl, G.; Carré, G.; Marquéz, J.R.G.; Gruber, B.; Lafourcade, B.; Leitão, P.J.; et al. Collinearity: A Review of Methods to Deal with It and a Simulation Study Evaluating Their Performance. Ecography 2013, 36, 27–46. [Google Scholar] [CrossRef]
  63. Chen, W.; Zhang, S.; Li, R.; Shahabi, H. Performance Evaluation of the GIS-Based Data Mining Techniques of Best-First Decision Tree, Random Forest, and Naïve Bayes Tree for Landslide Susceptibility Modeling. Sci. Total Environ. 2018, 644, 1006–1018. [Google Scholar] [CrossRef]
  64. Hong, H.; Liu, J.; Zhu, A.-X. Modeling Landslide Susceptibility Using LogitBoost Alternating Decision Trees and Forest by Penalizing Attributes with the Bagging Ensemble. Sci. Total Environ. 2020, 718, 137231. [Google Scholar] [CrossRef]
  65. Chen, W.; Yan, X.; Zhao, Z.; Hong, H.; Bui, D.T.; Pradhan, B. Spatial Prediction of Landslide Susceptibility Using Data Mining-Based Kernel Logistic Regression, Naive Bayes and RBFNetwork Models for the Long County Area (China). Bull. Eng. Geol. Environ. 2019, 78, 247–266. [Google Scholar] [CrossRef]
  66. Arabameri, A.; Asadi Nalivan, O.; Chandra Pal, S.; Chakrabortty, R.; Saha, A.; Lee, S.; Pradhan, B.; Tien Bui, D. Novel Machine Learning Approaches for Modelling the Gully Erosion Susceptibility. Remote Sens. 2020, 12, 2833. [Google Scholar] [CrossRef]
  67. Kadavi, P.R.; Lee, C.W.; Lee, S. Landslide-Susceptibility Mapping in Gangwon-Do, South Korea, Using Logistic Regression and Decision Tree Models. Environ. Earth Sci. 2019, 78, 116. [Google Scholar] [CrossRef]
  68. Pradhan, B.; Abokharima, M.H.; Jebur, M.N.; Tehrany, M.S. Land Subsidence Susceptibility Mapping at Kinta Valley (Malaysia) Using the Evidential Belief Function Model in GIS. Nat. Hazards 2014, 73, 1019–1042. [Google Scholar] [CrossRef]
  69. Sameen, M.I.; Pradhan, B.; Lee, S. Application of Convolutional Neural Networks Featuring Bayesian Optimization for Landslide Susceptibility Assessment. Catena 2020, 186, 104249. [Google Scholar] [CrossRef]
  70. Panahi, M.; Sadhasivam, N.; Pourghasemi, H.R.; Rezaie, F.; Lee, S. Spatial Prediction of Groundwater Potential Mapping Based on Convolutional Neural Network (CNN) and Support Vector Regression (SVR). J. Hydrol. 2020, 588, 125033. [Google Scholar] [CrossRef]
  71. Zhao, G.; Pang, B.; Xu, Z.; Peng, D.; Zuo, D. Urban Flood Susceptibility Assessment Based on Convolutional Neural Networks. J. Hydrol. 2020, 590, 125235. [Google Scholar] [CrossRef]
  72. Youssef, A.M.; Pradhan, B.; Dikshit, A.; Al-Katheri, M.M.; Matar, S.S.; Mahdi, A.M. Landslide Susceptibility Mapping Using CNN-1D and 2D Deep Learning Algorithms: Comparison of Their Performance at Asir Region, KSA. Bull. Eng. Geol. Environ. 2022, 81, 165. [Google Scholar] [CrossRef]
  73. Khosravi, K.; Panahi, M.; Golkarian, A.; Keesstra, S.D.; Saco, P.M.; Bui, D.T.; Lee, S. Convolutional Neural Network Approach for Spatial Prediction of Flood Hazard at National Scale of Iran. J. Hydrol. 2020, 591, 125552. [Google Scholar] [CrossRef]
  74. Thi Ngo, P.T.; Panahi, M.; Khosravi, K.; Ghorbanzadeh, O.; Kariminejad, N.; Cerda, A.; Lee, S. Evaluation of Deep Learning Algorithms for National Scale Landslide Susceptibility Mapping of Iran. Geosci. Front. 2021, 12, 505–519. [Google Scholar] [CrossRef]
  75. Lee, S.; Rezaie, F. Application of Statistical and Machine Learning Techniques for Habitat Potential Mapping of Siberian Roe Deer in South Korea. Proc. Natl. Inst. Ecol. Repub. Korea 2021, 2, 1–4. [Google Scholar] [CrossRef]
  76. Hakim, W.L.; Nur, A.S.; Rezaie, F.; Panahi, M.; Lee, C.W.; Lee, S. Convolutional Neural Network and Long Short-Term Memory Algorithms for Groundwater Potential Mapping in Anseong, South Korea. J. Hydrol. Reg. Stud. 2022, 39, 100990. [Google Scholar] [CrossRef]
  77. Mirjalili, S.; Mirjalili, S.M.; Lewis, A. Grey Wolf Optimizer. Adv. Eng. Softw. 2014, 69, 46–61. [Google Scholar] [CrossRef] [Green Version]
  78. Moayedi, H.; Osouli, A.; Tien Bui, D.; Foong, L.K. Spatial Landslide Susceptibility Assessment Based on Novel Neural-Metaheuristic Geographic Information System Based Ensembles. Sensors 2019, 19, 4698. [Google Scholar] [CrossRef] [PubMed]
  79. Atashpaz-Gargari, E.; Lucas, C. Imperialist Competitive Algorithm: An Algorithm for Optimization Inspired by Imperialistic Competition. In Proceedings of the 2007 IEEE Congress on Evolutionary Computation, Singapore, 25–28 September 2007; pp. 4661–4667. [Google Scholar]
  80. Wang, Y.; Hong, H.; Chen, W.; Li, S.; Panahi, M.; Khosravi, K.; Shirzadi, A.; Shahabi, H.; Panahi, S.; Costache, R. Flood Susceptibility Mapping in Dingnan County (China) Using Adaptive Neuro-Fuzzy Inference System with Biogeography Based Optimization and Imperialistic Competitive Algorithm. J. Environ. Manag. 2019, 247, 712–729. [Google Scholar] [CrossRef] [PubMed]
  81. Syifa, M.; Panahi, M.; Lee, C.-W. Mapping of Post-Wildfire Burned Area Using a Hybrid Algorithm and Satellite Data: The Case of the Camp Fire Wildfire in California, USA. Remote Sens. 2020, 12, 623. [Google Scholar] [CrossRef]
  82. Golkarian, A.; Naghibi, S.A.; Kalantar, B.; Pradhan, B. Groundwater Potential Mapping Using C5.0, Random Forest, and Multivariate Adaptive Regression Spline Models in GIS. Environ. Monit. Assess. 2018, 190, 149. [Google Scholar] [CrossRef]
  83. Fadhillah, M.F.; Lee, S.; Lee, C.-W.; Park, Y.-C. Application of Support Vector Regression and Metaheuristic Optimization Algorithms for Groundwater Potential Mapping in Gangneung-Si, South Korea. Remote Sens. 2021, 13, 1196. [Google Scholar] [CrossRef]
  84. Hakim, W.; Achmad, A.; Lee, C.-W. Land Subsidence Susceptibility Mapping in Jakarta Using Functional and Meta-Ensemble Machine Learning Algorithm Based on Time-Series InSAR Data. Remote Sens. 2020, 12, 3627. [Google Scholar] [CrossRef]
  85. Fadhillah, M.F.; Achmad, A.R.; Lee, C.-W. Integration of InSAR Time-Series Data and GIS to Assess Land Subsidence along Subway Lines in the Seoul Metropolitan Area, South Korea. Remote Sens. 2020, 12, 3505. [Google Scholar] [CrossRef]
  86. Pourghasemi, H.R.; Gayen, A.; Park, S.; Lee, C.W.; Lee, S. Assessment of Landslide-Prone Areas and Their Zonation Using Logistic Regression, LogitBoost, and Naïvebayes Machine-Learning Algorithms. Sustainability 2018, 10, 3697. [Google Scholar] [CrossRef]
  87. Pourghasemi, H.R.; Amiri, M.; Edalat, M.; Ahrari, A.H.; Panahi, M.; Sadhasivam, N.; Lee, S. Assessment of Urban Infrastructures Exposed to Flood Using Susceptibility Map and Google Earth Engine. IEEE J. Sel. Top. Appl. Earth Obs. Remote Sens. 2021, 14, 1923–1937. [Google Scholar] [CrossRef]
  88. Liu, Y.; Goodrick, S.; Heilman, W. Wildland Fire Emissions, Carbon, and Climate: Wildfire–Climate Interactions. For. Ecol. Manag. 2014, 317, 80–96. [Google Scholar] [CrossRef]
  89. Park, S.J.; Lee, C.W.; Lee, S.; Lee, M.J. Landslide Susceptibility Mapping and Comparison Using Decision Tree Models: A Case Study of Jumunjin Area, Korea. Remote Sens. 2018, 10, 1545. [Google Scholar] [CrossRef]
  90. Kadavi, P.; Lee, C.-W.; Lee, S. Application of Ensemble-Based Machine Learning Models to Landslide Susceptibility Mapping. Remote Sens. 2018, 10, 1252. [Google Scholar] [CrossRef]
  91. Panahi, M.; Gayen, A.; Pourghasemi, H.R.; Rezaie, F.; Lee, S. Spatial Prediction of Landslide Susceptibility Using Hybrid Support Vector Regression (SVR) and the Adaptive Neuro-Fuzzy Inference System (ANFIS) with Various Metaheuristic Algorithms. Sci. Total Environ. 2020, 741, 139937. [Google Scholar] [CrossRef] [PubMed]
  92. Fielding, E.J.; Liu, Z.; Stephenson, O.L.; Zhong, M.; Liang, C.; Moore, A.; Yun, S.H.; Simons, M. Surface Deformation Related to the 2019 Mw7.1 and 6.4 Ridgecrest Earthquakes in California from GPS, SAR Interferometry, and SAR Pixel Offsets. Seismol. Res. Lett. 2020, 91, 2035–2046. [Google Scholar] [CrossRef]
  93. Dodangeh, E.; Panahi, M.; Rezaie, F.; Lee, S.; Tien Bui, D.; Lee, C.W.; Pradhan, B. Novel Hybrid Intelligence Models for Flood-Susceptibility Prediction: Meta Optimization of the GMDH and SVR Models with the Genetic Algorithm and Harmony Search. J. Hydrol. 2020, 590, 125423. [Google Scholar] [CrossRef]
  94. 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]
  95. Hanberry, B.B.; Bragg, D.C.; Alexander, H.D. Open Forest Ecosystems: An Excluded State. For. Ecol. Manag. 2020, 472, 118256. [Google Scholar] [CrossRef]
  96. Hanberry, B.B. Classifying Large Wildfires in the United States by Land Cover. Remote Sens. 2020, 12, 2966. [Google Scholar] [CrossRef]
  97. Parks, S.A.; Abatzoglou, J.T. Warmer and Drier Fire Seasons Contribute to Increases in Area Burned at High Severity in Western US Forests From 1985 to 2017. Geophys. Res. Lett. 2020, 47, e2020GL089858. [Google Scholar] [CrossRef]
  98. Frankson, R.; Stevens, L.E.; Kunkel, K.E. California State Climate Summary 2022; Silver Spring: Montgomery, MD, USA, 2022. [Google Scholar]
  99. Kalantar, B.; Ueda, N.; Idrees, M.O.; Janizadeh, S.; Ahmadi, K.; Shabani, F. Forest Fire Susceptibility Prediction Based on Machine Learning Models with Resampling Algorithms on Remote Sensing Data. Remote Sens. 2020, 12, 3682. [Google Scholar] [CrossRef]
Figure 1. Flowchart of the overall methodology.
Figure 1. Flowchart of the overall methodology.
Remotesensing 14 04416 g001
Figure 2. The land use of the study area located in California, United States of America, and the distribution of historical wildfires from 2016 to 2020 in Plumas National Forest.
Figure 2. The land use of the study area located in California, United States of America, and the distribution of historical wildfires from 2016 to 2020 in Plumas National Forest.
Remotesensing 14 04416 g002
Figure 3. (a) Yearly frequency of wildfires and burned areas and (b) fire frequency by month in Plumas National Forest from 2016 to 2021.
Figure 3. (a) Yearly frequency of wildfires and burned areas and (b) fire frequency by month in Plumas National Forest from 2016 to 2021.
Remotesensing 14 04416 g003
Figure 4. Wildfire conditioning factors used in this study categorized as (ad) topographical, (eh) meteorological, (im) environmental, and (np) anthropological factors.
Figure 4. Wildfire conditioning factors used in this study categorized as (ad) topographical, (eh) meteorological, (im) environmental, and (np) anthropological factors.
Remotesensing 14 04416 g004aRemotesensing 14 04416 g004b
Figure 5. DPM depicting damage area after Camp Fire, Walker Fire, North Complex Fire, and Sheep Fire wildfires. Image from Sentinel-2 data acquired on 17 July 2021, during the Dixie Fire.
Figure 5. DPM depicting damage area after Camp Fire, Walker Fire, North Complex Fire, and Sheep Fire wildfires. Image from Sentinel-2 data acquired on 17 July 2021, during the Dixie Fire.
Remotesensing 14 04416 g005
Figure 6. Pearson correlation between wildfire conditioning factors.
Figure 6. Pearson correlation between wildfire conditioning factors.
Remotesensing 14 04416 g006
Figure 7. Wildfire susceptibility maps generated using 2016–2020 wildfire inventory and the (a) CNN, (b) CNN-GWO, and (c) CNN-ICA algorithms.
Figure 7. Wildfire susceptibility maps generated using 2016–2020 wildfire inventory and the (a) CNN, (b) CNN-GWO, and (c) CNN-ICA algorithms.
Remotesensing 14 04416 g007
Figure 8. The distribution of wildfire susceptibility levels in the CNN, CNN-GWO, and CNN-ICA models.
Figure 8. The distribution of wildfire susceptibility levels in the CNN, CNN-GWO, and CNN-ICA models.
Remotesensing 14 04416 g008
Figure 9. Analysis errors based on MSE, RMSE, error mean, and error StD using a testing dataset for: (a,d) CNN, (b,e) CNN-GWO, and (c,f) CNN-ICA models.
Figure 9. Analysis errors based on MSE, RMSE, error mean, and error StD using a testing dataset for: (a,d) CNN, (b,e) CNN-GWO, and (c,f) CNN-ICA models.
Remotesensing 14 04416 g009
Figure 10. The ROC curves and AUC values of CNN (red line), CNN-GWO (blue line), and CNN-ICA (green line) for wildfire susceptibility models using a testing dataset.
Figure 10. The ROC curves and AUC values of CNN (red line), CNN-GWO (blue line), and CNN-ICA (green line) for wildfire susceptibility models using a testing dataset.
Remotesensing 14 04416 g010
Figure 11. Comparison of the CNN-GWO model with wildfires that occurred in 2021.
Figure 11. Comparison of the CNN-GWO model with wildfires that occurred in 2021.
Remotesensing 14 04416 g011
Table 1. Statistical summary of wildfire ignition causes in Plumas National Forest from 2016 to 2020. (2021 wildfires are under investigation).
Table 1. Statistical summary of wildfire ignition causes in Plumas National Forest from 2016 to 2020. (2021 wildfires are under investigation).
CausesLightningHuman-CausedMiscellaneousUnknown
TransportationHuman ActivityConstruction
Number of wildfires253186810
Percentage35.714.2925.718.5711.4314.29
Table 2. Wildfire information and Sentinel-1 SAR data used in this study.
Table 2. Wildfire information and Sentinel-1 SAR data used in this study.
Fire NameAlarm Date(MM/DD/YYYY)Burned Area (ha)Pre-EventPost-EventFlight Direction
123
North Complex08/17/2020129,00408/01/202008/13/202009/06/2020Desc
Sheep08/17/202011,96708/01/202008/13/202008/25/2020Desc
Loyalton08/14/202019,03208/01/202008/13/202009/06/2020Desc
Walker09/04/201922,10708/20/201909/01/201909/13/2019Asc
Camp11/08/201862,05310/23/201811/04/201811/16/2018Desc
Cascade10/08/2017404209/10/201710/04/201710/16/2017Desc
Cherokee10/08/2017340609/10/201710/04/201710/16/2017Desc
Ponderosa08/29/2017162508/06/201708/18/201709/11/2017Asc
Minerva 507/29/2017174407/01/201707/13/201708/06/2017Asc
Wall07/07/2017244106/18/201706/30/201707/12/2017Desc
Saddle09/05/2016344080/5/201608/29/201609/22/2016Asc
Table 3. The list and description of wildfire conditioning factors.
Table 3. The list and description of wildfire conditioning factors.
CategoryFactorsScale/ResolutionSource of DataReferences
TopographyAspect30 mCopernicus DEM[2]
Altitude[40]
Slope[16]
Plan curvature[14]
MeteorologicalPrecipitation800 mPRISM[2]
Maximum temperature[39]
Solar radiance4 kmNREL[41]
Windspeed100 mGlobal Wind Atlas[42]
EnvironmentalDistance to stream1:5000California State Geoportal[16]
Drought index4 kmTerra Climate[43]
Soil moisture[44]
NDVI375 mMODIS[39]
Topographic wetness index30 mCopernicus DEM[40]
AnthropologicalLand use30 mUSGS[2]
Distance to road[14]
Distance to settlement[41]
Table 4. IGR and multicollinearity results for the conditioning factors.
Table 4. IGR and multicollinearity results for the conditioning factors.
FactorIGRCollinearity Statistics
TOLVIF
Altitude0.170.175.91
Aspect0.060.981.02
Distance to stream0.060.991.01
Distance to road0.140.971.03
Distance to settlement0.030.681.47
Land use0.390.392.55
NDVI0.080.711.40
Drought index0.250.244.13
Plan curvature0.040.831.20
Precipitation0.040.392.57
Slope0.160.303.31
Soil moisture0.180.273.66
Solar0.120.412.49
Maximum temperature0.180.224.63
TWI0.160.531.90
Windspeed0.030.811.22
Table 5. Frequency ratio of wildfire conditioning factors.
Table 5. Frequency ratio of wildfire conditioning factors.
FactorClassTotal %Event %Frequency Ratio
AspectNorth1.850.000.00
Northeast10.3010.581.03
East9.745.770.59
Southeast11.3112.821.13
South12.5613.781.10
Southwest14.3116.031.12
West14.7015.381.05
Northwest13.9913.140.94
Flat11.2512.501.11
Altitude (m)13–27219.8320.511.03
272–121620.1331.091.54
1216–152519.4316.990.87
1525–177320.4616.990.83
1773–254920.1514.420.72
Distance to Road (m)0–6021.2616.350.77
60–15220.1313.460.67
152–30119.7023.401.19
301–59519.4722.121.14
595–426019.4424.681.27
Distance to Settlement (m)0–28717.9012.820.72
287–114822.3428.531.28
1148–267920.2526.281.30
2679–564620.0816.670.83
5646–24,40419.4215.710.81
Distance to stream (m)0–99719.4122.761.17
997–216020.1424.041.19
2160–357320.9822.121.05
3573–565119.7515.710.80
5651–21,19219.7215.380.78
Land useOpen water2.190.080.04
Developed, open space1.801.070.59
Developed, low intensity0.840.440.52
Developed, medium intensity0.610.220.37
Developed, high intensity0.180.060.32
Deciduous forest0.480.200.41
Evergreen forest52.0956.721.09
Mixed forest0.240.241.00
Shrub/scrub19.3418.020.93
Herbaceous9.3522.862.44
Cultivated crops10.510.010.00
Woody wetlands0.500.070.14
Emergent herbaceous wetlands1.550.020.01
Barren land0.060.000.00
Hay/pasture0.260.000.00
Drought index−1.88–−0.8319.4618.910.97
−0.84–−0.3220.3818.590.91
−0.33–−0.0819.739.290.47
−0.09–0.2121.1526.601.26
0.21–0.8919.2826.601.38
NDVI0–0.2520.0017.630.88
0.25–0.4620.0016.670.83
0.46–0.6120.0120.831.04
0.61–0.7620.0022.121.10
0.76–119.9822.761.14
Plan curvatureConcave33.8538.141.13
Flat16.1011.860.74
Convex50.0550.001.00
Precipitation (mm)13.93–25.2120.0117.310.87
25.21–33.7220.0024.681.23
33.72–43.3520.0022.761.14
43.35–67.2520.0019.550.98
67.25–116.1019.9915.710.79
Slope (degree)0–1.8019.7812.180.62
1.80–7.2220.5918.590.90
7.22–12.8919.9721.151.06
12.89–20.3719.8419.550.99
20.37–65.7519.8228.531.44
Soil moisture (mm)22.00–95.2919.7621.791.10
95.29–184.2920.1612.820.64
184.29–252.3519.9815.710.79
252.35–348.3319.5025.321.30
348.33–467.0020.6024.361.18
Solar6.14–6.3821.0927.241.29
6.38–6.5521.2327.881.31
6.55–6.7221.0019.550.93
6.72–6.8418.5012.500.68
6.84–7.2718.0512.820.71
Tmax (°C)52.07–59.0720.0612.820.64
59.07–61.8420.0115.710.78
61.84–65.6720.0118.590.93
65.67–74.1119.9730.131.51
74.11–76.0619.9522.761.14
TWI3.29–5.5218.1821.151.16
5.52–6.2520.0522.121.10
6.25–7.2720.8821.471.03
7.27–9.2021.0524.361.15
9.20–26.3119.8410.900.55
Windspeed (m/s)0.99–3.9320.0020.511.03
3.93–4.6720.0021.151.06
4.67–5.2920.0020.191.00
5.29–5.8420.0022.121.11
5.84–11.8920.0016.030.80
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Nur, A.S.; Kim, Y.J.; Lee, C.-W. Creation of Wildfire Susceptibility Maps in Plumas National Forest Using InSAR Coherence, Deep Learning, and Metaheuristic Optimization Approaches. Remote Sens. 2022, 14, 4416. https://doi.org/10.3390/rs14174416

AMA Style

Nur AS, Kim YJ, Lee C-W. Creation of Wildfire Susceptibility Maps in Plumas National Forest Using InSAR Coherence, Deep Learning, and Metaheuristic Optimization Approaches. Remote Sensing. 2022; 14(17):4416. https://doi.org/10.3390/rs14174416

Chicago/Turabian Style

Nur, Arip Syaripudin, Yong Je Kim, and Chang-Wook Lee. 2022. "Creation of Wildfire Susceptibility Maps in Plumas National Forest Using InSAR Coherence, Deep Learning, and Metaheuristic Optimization Approaches" Remote Sensing 14, no. 17: 4416. https://doi.org/10.3390/rs14174416

APA Style

Nur, A. S., Kim, Y. J., & Lee, C. -W. (2022). Creation of Wildfire Susceptibility Maps in Plumas National Forest Using InSAR Coherence, Deep Learning, and Metaheuristic Optimization Approaches. Remote Sensing, 14(17), 4416. https://doi.org/10.3390/rs14174416

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