Next Article in Journal
Secondary Succession in the Tropical Lowland Rainforest Reduced the Stochasticity of Soil Bacterial Communities through the Stability of Plant Communities
Next Article in Special Issue
Characterizing Global Fire Regimes from Satellite-Derived Products
Previous Article in Journal
Spatial Distribution and Regulating Factors of Soil Nutrient Stocks in Afforested Dump of Pingshuo Opencast Coalmine, China
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Burned Area Detection Using Multi-Sensor SAR, Optical, and Thermal Data in Mediterranean Pine Forest

1
Department of Geomatics Engineering, Hacettepe University, Ankara 06800, Turkey
2
Department of Geomatics Engineering, Zonguldak Bulent Ecevit University, Zonguldak 67100, Turkey
3
Department of Architecture, Cukurova University, Adana 01250, Turkey
4
Department of Geomatics Engineering, Istanbul Technical University, Maslak, Istanbul 34469, Turkey
5
Department of GIS and Remote Sensing, University of Tabriz, Tabriz 5166616471, Iran
6
Remote Sensing Laboratory, University of Tabriz, Tabriz 5166616471, Iran
7
Department of Architecture and Building Engineering, Tokyo Institute of Technology, Yokohama 226-8502, Japan
8
Department of Geomatic Engineering, Yildiz Technical University, Istanbul 34210, Turkey
*
Author to whom correspondence should be addressed.
Submission received: 31 December 2021 / Revised: 15 February 2022 / Accepted: 16 February 2022 / Published: 18 February 2022

Abstract

:
Burned area (BA) mapping of a forest after a fire is required for its management and the determination of the impacts on ecosystems. Different remote sensing sensors and their combinations have been used due to their individual limitations for accurate BA mapping. This study analyzes the contribution of different features derived from optical, thermal, and Synthetic Aperture Radar (SAR) images to extract BA information from the Turkish red pine (Pinus brutia Ten.) forest in a Mediterranean ecosystem. In addition to reflectance values of the optical images, Normalized Burn Ratio (NBR) and Land Surface Temperature (LST) data are produced from both Sentinel-2 and Landsat-8 data. The backscatter of C-band Sentinel-1 and L-band ALOS-2 SAR images and the coherence feature derived from the Interferometric SAR technique were also used. The pixel-based random forest image classification method is applied to classify the BA detection in 24 scenarios created using these features. The results show that the L-band data provided a better contribution than C-band data and the combination of features created from Landsat LST, NBR, and coherence of L-band ALOS-2 achieved the highest accuracy, with an overall accuracy of 96% and a Kappa coefficient of 92.62%.

Graphical Abstract

1. Introduction

World forests provide environmental, social, and ecological benefits as well as their marketing values. They affect the regional and global climate by means of biological, chemical, and physical processes that influence atmospheric composition, hydrologic cycle, and planetary energetics [1]. In addition to the forests’ benefits to the climate, they help all living creatures, especially humankind, by filtering the water, purifying the air, and preventing erosion [2]. Furthermore, forests are home to numerous plants and animals, and they supply fundamental natural resources from timber and food to medicinal plants [3]. However, wildfires affect around 350 million hectares of land annually [4,5], and these fires result in the loss of huge amounts of forest cover. The fundamental causes of forest fires can be categorized into: (i) human activities such as land use land cover (LULC) changes, campfires, smoking, etc., and (ii) natural conditions such as lightning, climatic conditions, vegetation status, and topography [6,7,8]. Over the past five years, the world has witnessed a large number of wildfires, and the responses of the increased fire activities are a rise in fire intensity, changes in climate, reduced moisture availability (dryer summer seasons), higher greenhouse gas emissions, and a loss of various plant and animal species [9,10]. Moreover, forest fires have huge impacts on the economy of some countries [11].
After a fire, determining the total affected area, the grade of the damage, and appropriate methods for vegetation rehabilitation are the primary goals of forest management [6]. Spatio-temporal mapping of the burned areas (BAs) is crucial for ecological and environmental impact assessment, determination of economic losses, and nowadays, modeling climatic influences of the burning of biomass [12,13]. Therefore, fire management agencies and local governments need timely, accurate/reliable, and cost-effective fire mapping systems, especially under the current situation of the increasing frequency and severity of wildfire events [4,14]. Concerning the BA mapping, Koutsias and Karteris [15] introduced three basic approaches depending on the scale level of measurements, called micro-scale, meso-scale, and macro-scale. The micro-scale level refers to in situ human-made measurements, which are highly accurate but expensive and time-consuming. The meso-scale level represents using aerial photography at a scale ranging from 1:10,000 to 1:25,000. At the macro-scale level, satellite-based remote sensing data with various spatial resolutions starting from 10 meters are used. Satellite imageries with their high spatial and temporal (revisit time) resolutions play an important role in BA mapping, especially in large extent areas with low accessibility. Furthermore, satellite data make it possible to generate continuous and long-term BA databases [16].
Considering the literature, the main remote sensing techniques for BA mapping include supervised and unsupervised classification methods, spectral indices, spectral mixture analysis, logistic regression modeling, principal component analysis, and recently, machine learning [7,13,15,17,18,19,20]. The datasets used in these methods are generally bi-temporal images or single post-fire images. In the bi-temporal approach, the different images obtained pre- and post-fire are used to reveal the differences in the spectral response of the BA. The most commonly used process of the bi-temporal approach includes two steps: (i) spectral index extraction, where popular ones are the Normalized Difference Vegetation Index (NDVI) and Normalized Burn Ratio (NBR), and (ii) different images between pre- and post-fire indices, namely, dNDVI and dNBR [11,21,22,23,24]. Apart from NDVI and NBR, the Mid-Infrared Burned Index (MIRBI), Global Environmental Monitoring Index (GEMI), Burned Area Index Modified (BAIM), and Char Soil Index (CSI) are among the other indices used for BA mapping [25,26,27,28]. Specifically, dNBR has been used for burn-severity analysis as the best measure [11,29,30,31]. Mallinis et al. [30] compared different spectral indices derived from Sentinel-2 and Landsat-8 for severity analysis in a Mediterranean pine forest. It was concluded that dNBR has a high relationship with field data and Sentinel-2 provided slightly better accuracy of classification. Achour et al. [31] applied dNBR and dNBR-based indices calculated from Sentinel-2 and Landsat-8 over Tunisia. Similarly, dNBR and the relativized version of it determined good results. It was also concluded that Sentinel-2 was found to be slightly more effective than Landsat-8 data for BA mapping.
Previous studies revealed the ability of optical, thermal, and radar remote sensing data to be used in BA mappings [4,9,32,33,34,35]. The Moderate-Resolution Imaging Spectrometer (MODIS), The Advanced Very High-Resolution Radiometer (AVHRR), and Visible Infrared Imaging Radiometer Suite (VIIRS) are the preliminary sensors used for BA mapping, whose spatial resolutions range between 300 m and 1 km. Although these sensors have very high temporal resolution (up to two or three revisits in a day), they have low spatial resolutions to present the exact extent of the BA, especially for small extent BAs [33,34,35,36,37,38]. On the other hand, with their relatively high spatial resolutions but low temporal resolutions, Landsat satellites and Sentinel-2 provide high-quality images for BA mapping at a satisfactory extent [39,40]. All of the above instruments provide optical and/or thermal images, and one major drawback of optical and thermal remote sensing is the presence of cloud cover, smoke, and/or haze, which is an important limitation for many areas, such as in tropical regions. Therefore, researchers investigated the capability of Synthetic Aperture Radar (SAR) and Interferometric SAR (InSAR) data for mapping the BA since the active microwave has the ability to penetrate through the cloud coverage and it is less affected by weather conditions [4,32,41,42,43]. Moreover, some studies considered a multi-sensor approach using SAR and optical/thermal data for BA mapping [4,33,44,45]. Philipp and Levick [10] compared Sentinel-2 and Sentinel-1 mapping results of fire scar over tropical savannas. Stroppiana et al. [33] integrated backscatter of C-band ENVISAT data with Landsat TM data to extract burned and unburned regions. Zhang et al. [4] used backscatter, coherence, and reflectance of Sentinel-2 data for the BA mapping. Liu et al. [23] combined the Land Surface Temperature (LST) of Landsat-8 with the spectral bands to identify the BAs.
The Mediterranean is very vulnerable to fire risk and about 57% of Turkey’s forests are under this risk. Forests in Turkey cover an area of about 22.9 million ha, which is 29% of the country. Coniferous forest covers 48%, broadleaved covers 32%, and mixed forest covers the remaining 20%. Pinus brutia Ten. is the largest distributed coniferous tree species that covers about 22.74% of the total amount of forest in Turkey [46]. As a coniferous tree species native to the Mediterranean, Pinus brutia Ten. (Red pine) is located in the western part of the Mediterranean and Pinus Halepensis Mill. (Aleppo pine) is located in the eastern part. Together, they cover an area of approximately 6.8 million hectares in total [47]. They are distributed at the coastal zones where they provide high ecological characteristics. These species are preferred in post-fire afforestation owing to their resistance to dry conditions and their ability to regenerate in a short time. They also have economic importance such as pulp production [46,47,48]. However, they are also the species most exposed to forest fires in Europe [49]. According to data from the General Directorate of Forestry (OGM), in Turkey, the number of wildfires has increased from around 1000 to around 3500 per year since the 1940s [50]. It is seen that the majority (~90%) of the fires that occurred in 2020 and 2019 were human-induced. In 2019, one fire greater than 1000 ha and one of four fires larger than 1000 ha in 2020 occurred in Izmir Region. The number of fires is higher between August and October compared to other months [51,52]. In July and August 2021, 1700 km2 of forest in Turkey’s Mediterranean Region was affected by a series of more than two hundred wildfires, meaning the worst-ever wildfire season in the country’s history [53]. Therefore, it is crucial for Turkey to determine the accurate extent of the forest fires for post-fire rehabilitation. The major aim of this study is to map BA using multi-sensor data such as SAR data, optical data, and thermal data with the random forest (RA) classifier. SAR data include backscatter images from ALOS-2 and Sentinel-1, and interferometric coherence from ALOS-2. As optical data, reflective bands and spectral indices of Sentinel-2 and Landsat-8 were implemented, and LST images from the Landsat-8 Thermal Infrared Sensor (TIRS) were utilized. Twenty-four different scenarios were determined based on these datasets, and as a sub-aim, we investigated the contribution of different sensor data on the accuracy of BA mapping with the use of the RF classifier. To the best of our knowledge, so far, no study has considered such a rich database including SAR backscatter, interferometric coherence, optical data, and thermal data to investigate their effectiveness in BA mapping. In addition to the BA mapping, the burn severity maps were generated from both Sentinel-2- and Landsat-8-derived NBR to reveal the grades of the fire damage, and to compare the results of both sensors.

2. Materials and Methods

2.1. Study Area

The fire event considered in this study occurred on 18 August 2019 in the district of İzmir metropolitan city (Figure 1). The region where the forest fire was located has a Mediterranean climate. The forest is dominated by coniferous Turkish red pine (Pinus brutia Ten.) and maquis vegetation cover [54]. The elevation ranges from 140 to 975 m. Approximately 5000 ha was affected by the fire [24,54].

2.2. Multi-Sensor Data

In this study, a multi-sensor approach was proposed to determine the post-fire BA detection. For the analysis, visible, near infrared (NIR), and thermal bands of optical and thermal sensors and multi-frequency SAR data were used (Table 1). The optical satellites, Sentinel-2 and Landsat-8, of pre- and post-fire images were acquired from the archives of ESA (https://scihub.copernicus.eu/) (accessed on 1 September 2021) and NASA (https://earthexplorer.usgs.gov/) (accessed on 1 September 2021). Similarly, pre- and post-fire images of C-band Sentinel-1(https://scihub.copernicus.eu/) (accessed on 1 September 2021)and L-band ALOS-2 data were used as SAR images. As the time difference between pre- and post-event acquisitions of ALOS-2 data includes a long period, the same acquisition period was also preferred for the other datasets. Different features from optical, thermal, and SAR data were obtained and their contribution to the BA extraction was tested with the RF classification method. A generalized workflow of the study is presented in Figure 2.

2.3. Burn Severity Spectral Indices of Optical Images

Remote sensing-based burn severity indices have been developed and used due to their simple computation and direct application [55]. Spectral indices based on the Near Infrared (NIR) and Short-Wave Infrared (SWIR) bands, specifically the Normalized Burn Ratio (NBR) and its bi-temporal approach, such as the differenced Normalized Burn Ratio (dNBR), have been identified as one of the optimal burn severity measures [56,57].
In this study, to determine the burn severity of the selected region, NBR and dNBR were calculated by using pre- and post-fire Landsat-8 OLI and Sentinel-2 MS data (Table 2). BA severity can be easily classified by the dNBR index. The dNBR is obtained by using the bi-temporal difference of the preprocessed NBR images [29].
Hypothetically, the values for dNBR are between −2 and +2, where the values for BAs are ranging from 0.10 to 1.35 and the unburned areas from −0.10 to +0.10 in most areas, but the minimum and maximum values can be different based on the study area (Table 3) [58]. These ranges are flexible and scene-pair-dependent, and shifts in thresholds (±0.1) are possible. dNBR values of less than about −0.55, or greater than about +1.35 may occur, but usually are not considered burned. Rather, they are likely anomalies caused by miss-registration, clouds, or other factors not related to real land cover differences.

2.4. LST Extraction

LST images were obtained from the Landsat-8 TIR-1 band using the Mono Window Algorithm (MWA) proposed by Qin et al. [60]. As the Land Surface Emissivity (LSE) model in LST algorithm, we considered the NDVI threshold-based model introduced by Sobrino et al. [61] since a previous LST accuracy assessment study provided the highest results with Landsat-8 by this combination [62]. A difference LST (dLST) is created using the difference between pre- and post-fire Landsat-8 data.

2.5. Preprocessing of SAR Images

The Interferometric Wide (IW) mode Sentinel-1 and Strip Map (SM) mode ALOS-2 data were acquired in Single-Look Complex (SLC) image format. For each data point, the backscatter value was obtained after calibration, filtering by the 5 × 5 Lee filter and geometric correction with the Shuttle Radar Topography Mission (SRTM) digital elevation model (30 m). The calibration step converts the digital number of the pixels ( D N i ) into the backscattering coefficient ( σ p q 0 ). In the below equations, K is the calibration constant and σ p q 0 is the backscattering coefficient, also called Sigma Nought. The backscattering coefficients which have linear units were converted into decibel ( σ d B 0 ) units using Equation (3). A differential image was created to identify the BAs using pre-event backscatter ( σ p r e f i r e ) and post-event backscatter ( σ p o s t f i r e ) values for each polarimetry. All images were resampled to a common geometric resolution of 30 m.
σ p q 0 = K · D N i 2
σ p q 0 = | D N i | 2 A i 2
σ d B 0 = 10 · l o g 10 ( σ p q 0 )
d σ d B 0 = σ p r e f i r e σ p o s t f i r e
As the SAR data were collected in the complex image format, the SLC images have both phase and amplitude information. Using two co-registered SLC images, a cross-correlation factor between the pairs can be extracted with using Equation (5):
γ = | S 1 S 2 * | ( S 1 S 1 * S 2 S 2 * )
where γ represents coherence, S1 and S2 represent complex matrices of the backscatter coefficient of the pre- and post-SAR images, * indicates complex conjugation, ‹› indicates spatial averaging, and | | indicates the absolute value [63,64]. The coherence values range between zero (low coherence) and one (high coherence). Low coherence indicates a high degree of variation between repeat-pass SAR data. Coherence of Sentinel-1 data was not used further because its values were also very low for non-burned regions. The open-source tools of the Sentinel Application Platform [65] software were used for the preprocessing steps of the SAR data.

2.6. Land Cover Classification

The supervised pixel-based RF classification approach, which represents robust performance [66], was used to detect land cover classes, namely, burned forest, unburned forest, agriculture, and bare land. The general information of RF has been provided extensively in several studies [67,68]. In brief, the RF approach is mainly a bagging algorithm that adds the randomness feature to the model.
The RF algorithm is based on the creation of many decision trees, combining them to obtain more accurate and reliable results. It uses a series of binary decision trees and it branches each node by using the best among the randomly selected variables at each node. The Information Gain Rate criterion and the Gini Index are the attribute selection measures most frequently used to induce the decision tree. We chose the Gini Index, which measures the impurity of an attribute in relation to the classes. For a given T training set, it randomly selects a case (pixel) and determines the class that it belongs to. In this work, the optimal number of the trees in the RF classification was applied as 100. This number was determined based on the trial-and-error method, and only this number was used in each scenario to ensure the consistency of the research.
The initialization of a supervised classification process requires composite images and training samples (polygons). The training polygons selected for four classes as BAs, unburned areas, agriculture, and bare land in the composite images were used to classify the images. Training (n = 66) and testing (n = 50) samples were collected randomly on high-resolution images of Google Earth. The number of training pixels were 965, 3660, 546, and 371 for burned forest, unburned forest, agriculture, and bare land, respectively. The number of test pixels were 678, 2344, 388, and 242 for burned forest, unburned forest, agriculture, and bare land, respectively.
A total of 24 scenarios were described to perform the contribution of different features to extract the fire BAs (Table 4). For this purpose, each feature was added as an additional band to the multispectral bands and their contribution was examined. As a last step, only contributing features were added to identify the contribution of these features to the BA analysis. The resulting BA maps were evaluated with an error matrix using test samples [68]. The EnMap-Box integrated to QGIS was used for the classification and evaluations [69,70,71].
Accuracy assessments of all of the classification results were calculated to determine the performance of the selected sensors and created scenarios by using the error matrix.
Overall accuracy (OA), user’s accuracy (UA), and producer’s accuracy (PA), the Kappa coefficient [72], and the F1-score were calculated to assess the performance of the random forest classification. The F1-score provides an accuracy that combines the UA and PA [73]. Google Earth images were used to collect a validation dataset to determine the accuracy of the classification results.

3. Results

The contribution of different features derived from multi-sensors was evaluated with the comparison of the accuracy analysis (Figure 3 and Figure 4). In general, all scenarios provided an OA higher than 91%, a Kappa value higher than 83.9%, and an F1 value higher than 83%. In the first twelve scenarios, post-fire Sentinel-2 with VIS and NIR bands were combined with features derived from Sentinel-2 (NBR2020, dNBR), Landsat-8 (dLST), and SAR images (Table 4). The dNBR improved the accuracy better than the index of S-2 NBR2020. The backscattering coefficients of ALOS-2 and Sentinel-1 could not improve the accuracy, but coherence of ALOS-2 enhanced the accuracy by 1.64%. Including features of optic and SAR images, the OA and Kappa were improved by 3.78% and 6.62%, respectively, by S12.
Similar to Sentinel-2, the dNBR improved the accuracy better than L-8 NBR2020 when combined with Landsat-8. Both dB values of ALOS-2 and Sentinel-1 could not contribute to the S13. However, similar to Sentinel-2, Landsat-8 combined coherence improved the results. The combination of all data (S24) provided an increase of 3.2% in the OA and 5.54% in the Kappa values. Including dLST and coherence of ALOS-2, the results slightly improved the UA from 92.85% to 93.23% and decreased the PA from 97.64% to 97.49%. Landsat-8 increased the UA of BA by 3.8%. S24 provided the highest PA and F1 among all scenarios (Figure 5).
The highly accurate classification results which were provided by the S12 and S24 are mapped in Figure 5. When comparing the contributions of the Sentinel-2 and Landsat-8 combined datasets, it showed that the results have similarity and the two sensors confirmed each other. The highest accuracy obtained with S24 indicates that the BA was found to be about 5742 ha. About 97% of the BA was covered with Turkish Red pine forest. The results of the burn severity analysis showed that the BA is mostly covered with moderate- to high-severity levels and Sentinel-2 data provided higher severity regions (3107.87 ha) than Landsat-8 data (1292.42 ha). The moderate–high-severity levels of Sentinel-2 (2066.87 ha) are closer to the result of Landsat-8 (2447.33 ha) (Figure 6).
The severity map was analyzed with a slope map that was derived from ALOS World 3D [74]. The data have 30 m spatial resolution. The overlapped region was divided to five classes of slope, as gentle (0°–5°), low (5°–15°), medium (15°–25°), high (25°–35°), and steep (>35°). In places where the slope is close to flat, it was seen that the severity of the fire was low, and as the slope increases to a medium degree, the severity of the fire also increased. On the contrary, it was determined that the burn severity decreased from moderate to steep slopes in both cases (Figure 7).

4. Discussions

The spectral indices of optical images have been widely used for BA extraction. In addition to different indices derived from individual satellites [75], the combination of backscatter of SAR and indices derived from optical data [33], the combination of backscatter, coherence, and reflectance of Sentinel-2 [4], and the combined approach of LST of Landsat-8 with the spectral bands improved BA mapping [23]. While other studies indicated success of coherence data, Philipp and Levick [10] failed to add a contribution compared to backscatter data.
In this study, different from the above-mentioned studies, spectral bands of optical data, spectral indices, and LST were combined with SAR features to improve the accuracy of BA detection. Among all of the features, dLST and ALOS-2 individually contributed more than the other features. Only using NBR2020 (S2 and S14) presented almost similar results to the results of Sentinel-2 (S1) and Landsat-8 (S13). The contribution of NBR2019 could be seen when dNBR was added to the reflectance data. While considering the polarimetry, it was noticed that each σ d B 0 (bakscatter difference) of polarimetric data could not contribute to the classification accuracy and the OA was decreased.
In general, the results of Sentinel-2 (S1 to S12) and L-8 (S13 to S24) showed that Landsat-8 provided slightly higher accuracies than Sentinel-2. Considering the class of burned forest, the accuracies were close to each other, which is similar to the spatial distribution pattern of burned areas in two products (S12 and S24). The results showed that the scenario combined with Sentinel-2 (S12) detected a slightly higher BA (5813 ha) than the scenario combined with Landsat-8 (S24) (5742 ha). The BA of Turkish Red pine was about 98%, whereas 10% was degraded forest in both results. The rest is composed of agriculture (~55 ha in S12 and ~50 ha in S24) and maquis (~59 ha in S12 and ~60 ha in S24). Ngadze et al. [40] compared RF classified Sentinel-2 and Landsat-8 images over two areas and showed that Sentinel-2 had higher accuracies on one site while Landsat-8 accuracy was higher on the other site. They report that the reason for their result may be due to the acquisition periods of the satellite images. In our case, there was also a time difference between the two satellite image acquisitions. Another reason might also be due to the up-sampling of Sentinel-2 images to Landsat-8 resolution. When comparing bands, the thermal band added scenarios (S4 and S16) had higher contributions than scenarios with a combination of NIR and SWIR bands (S1 and S13), as seen in both Sentinel-2 and Landsat-8 datasets.
Differences between wavelengths showed that L-band ALOS-2 data illustrated a higher change after the fire compared to C-band Sentinel-1 data. The value of σ d B 0 represents the sum of the intensity collected by the sensor after backscattering. The γ (coherence) presents the direct measurement of distributed scatters obtained by the phase information using SAR interferometry [9]. As coherence is widely affected by the temporal decorrelation pre- and post-fire, the SAR pair should provide low coherence at the BA compared to the unburned forest area. However, Sentinel-1 did not show a low-level coherence value at the burned region and it was not included in the processing. On the other hand, ALOS-2 data provided a high level of coherence loss, which added a contribution to discriminate the BA. Tanase et al. [9] showed that the sensitivity of L-band ALOS-1 data is higher than other bands when using almost one and a half month acquisition periods over the pine forest of Spain. In our study, even though the temporal acquisition periods were almost one year, the coherence of L-band ALOS-2 for both cross- and co-polarized data indicated a sensitivity to determine BA over the Turkish red pine-dominated forest region. Sentinel-1 showed limitations for the extraction of BA over a one-year image pair acquisition. Compared to the short-wavelength C-band, the longer wavelength L-band indicated higher coherence even when the temporal baseline was about one year. However, the backscatter of ALOS-2 could not improve the accuracies when it was included.
The burn severity analysis using dNBR showed that Sentinel-2 resulted in a higher level of severity compared to the Landsat-8 data. The differences between the two severity maps might be due to the differences in image acquisition periods and the spatial resolutions [30]. The burn severity and its distribution on different slope levels showed that Sentinel-2 and Landsat-8 provided similar trends: 42% and 41.8% of burn severity were located over the medium slope using Sentinel-2 and Landsat-8, respectively.
The study area was also analyzed by the previous studies that focused on merely optical datasets. Nasery and Kalkan [76] used Sentinel-2 data with dNBR and dSAVI indices to assess the burn severity. The results indicated that the dNBR index provided better discrimination of severity levels than dSAVI. However, the detected BA using both indices (dNBR ~ 6910 ha, dSAVI ~ 6700 ha) was not closer to the reference area (~5000 ha) presented in the previous studies [24,54].
Atak and Tonyaloglu [54] evaluated five indices, namely, dNDVI, dARVI, dBAI, dNBR, and dNBR2, calculated from the pre- and post-fire images of L-8. An object-based classification was applied to detect burned and unburned regions using S-2 data as the reference data. The results indicated that dNDVI achieved the closest area (~6204 ha) that was detected, however dNBR2 and dARVI provided the best results of classification accuracy, with an OA of 84%, a UA of 90%, and a PA of 90% in both indices.
Masshadi and Alganci [24] evaluated four indices, namely dNDVI, dNBR, RdNBR, and RBR, for burn severity mapping derived from both Sentinel-2 and Landsat-8 datasets. The RdNBR index presented the highest accuracies in both Landsat-8 (89%) and Sentinel-2 (93%) data with respect to the active fire dataset of NASA. However, the estimated areas were not close to the declared area as the RdNBR results from Sentinel-2 and Landsat-8 were ~7979.4 and ~6669.6 ha, respectively. The study also applied SVM and RF classification approaches on reflectance images of Sentinel-2 and Landsat-8 to extract the BA. The results showed that SVM estimated a larger BA than the RF approach. SVM provided better accuracies (OA = 89%, Kappa = 0.87) compared to RF (OA = 86%, Kappa = 0.85). However, RF better-estimated the BA using Landsat-8 (~5722.5 ha) than Sentinel-2 (~5790.8 ha) data.
Compared to the previous studies over the same region, this study provided a higher OA (98.92%) with a multi-sensor approach using the RF method. Masshadi and Alganci [24] estimated a BA slightly closer to the reference area; however, the OA and Kappa accuracies were lower than in this study. This might be due to the different acquisition periods of the optical data. Additionally, the selected and distributed numbers of training and test pixels were not same. In this study, in addition to the major portion of the BA, a small BA located in the south-west direction was also extracted (Figure 5). The previous studies, conducted over the same test site, merely used spectral indices of Sentinel-2 for severity analysis [76] and spectral indices of Sentinel-2 with object-based image classification [54]. In addition, a comparison was applied using spectral indices of Sentinel-2 and Landsat-8 data and pixel-based classification of spectral bands with SVM and RF methods [24]. As to these previous methods, our study demonstrated that the additional features from different satellite sensors can improve the accuracy of BA extraction using the pixel-based RF method. Concerning the other multi-sensor-based studies about BA mappings, the combination of multispectral and thermal bands of optical data [23], integration of optical and SAR data [33], SAR backscatter, coherence map, and optical spectral information [4] all showed the advantages of multi-sensor analysis. Moreover, in this study, not only the thermal band of optical satellites but also the coherence map of SAR data played an important role when they were combined with multispectral bands of optical data for accurate BA detection. Additionally, spatial analysis was applied to the burn severity map to show the spatial distribution of severity levels at different slopes.

5. Conclusions

This study analyzed the performance of different features which were extracted from optical and SAR data over a BA in a Mediterranean coniferous forest. The analysis showed that the use of a thermal band to extract the BA contributed more than the NIR and SWIR bands of optical data. The results also indicated that Sentinel-2 and Landsat-8 provided compatible results and they can be used as alternatives of each other.
The long-wavelength L-band SAR data provided a better contribution than C-band SAR data when the InSAR-based coherence feature was considered. The extraction of the BA was affected by the long temporal resolution between image acquisitions. L-band ALOS-2 was more sensitive to the changes, and the loss of coherence more clearly showed the changes between pre-fire and post-fire.
The study highlighted the contribution of the combined use of LST and coherence for mapping burnt regions in Mediterranean forests. In addition to optical VNIR bands, we used features derived from optical images, namely NBR, LST, and SAR data, as coherence information, which have not been evaluated all together for forest fire burn analysis in the literature.
Further analysis will focus on features derived from full polarimetric data in combination with optical features to identify the forest reforestation using multi-temporal analysis. The new L-band NISAR mission is expected to contribute to this issue with its short revisit time and open data policy.

Author Contributions

Conceptualization, S.A., C.B., A.S., F.B.B. and F.B.S.; methodology, S.A., C.B., A.S., F.B.B. and F.B.S.; analysis and validation, S.A., C.B., A.S. and F.B.B.; writing—original draft preparation, S.A., A.S. and F.B.B.; writing—review and editing, S.A., C.B., A.S., F.B.B., S.K., M.M. and F.B.S.; visualization, C.B. 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

The Copernicus Sentinel-1 and Sentinel-2 are openly available on the website of https://scihub.copernicus.eu/, accessed on 1 September 2021. Landsat-8 and ALOS World 3D data used in this study are openly available in the websites of NASA (https://earthexplorer.usgs.gov/, accessed on 1 September 2021) and JAXA (https://www.eorc.jaxa.jp/ALOS/en/dataset/aw3d30/aw3d30_e.htm, accessed on 1 September 2021), respectively.

Acknowledgments

The authors would like to thank Japan Aerospace Exploration Agency (JAXA), which provided the ALOS-2 data through the JAXA working group and ALOS World 3D data. The authors would like to thank United States Geological Survey (USGS) and the National Aeronautics and Space Administration (NASA) for providing Landsat-8 images at free of cost. Besides, the authors would like to thank European Space Agency (ESA) for providing Sentinel-1 and Sentinel-2 images freely available.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Table A1. Accuracy analysis of each scenario. The best result of each accuracy metric is written in bold.
Table A1. Accuracy analysis of each scenario. The best result of each accuracy metric is written in bold.
ScenariosFeaturesOA (%)Kappa (%)F1 Accuracy (%)
S1S-2 (B2, B3, B4, B8, B12)91.6584.9684.33
S2S-2 (B2, B3, B4, B8, B12), NBR202091.7085.0884.59
S3S-2 (B2, B3, B4, B8, B12), dNBR92.9987.3386.17
S4S-2 (B2, B3, B4, B8, B12), dLST93.4688.1387.62
S5 S - 2   ( B 2 ,   B 3 ,   B 4 ,   B 8 ,   B 12 ) ,   A - 2 γ H V 93.2087.5886.31
S6 S - 2   ( B 2 ,   B 3 ,   B 4 ,   B 8 ,   B 12 ) ,   A - 2 γ H H 93.2987.6886.36
S7 S - 2   ( B 2 ,   B 3 ,   B 4 ,   B 8 ,   B 12 ) ,   A - 2 d σ H V 0 91.2184.2183.70
S8 S - 2   ( B 2 ,   B 3 ,   B 4 ,   B 8 ,   B 12 ) ,   A - 2 d σ H V 0 91.0583.9083.19
S9 S - 2   ( B 2 ,   B 3 ,   B 4 ,   B 8 ,   B 12 ) ,   S - 1 d σ V V 0 91.1384.0783.59
S10 S - 2   ( B 2 ,   B 3 ,   B 4 ,   B 8 ,   B 12 ) ,   S - 1 d σ V H 0 91.5484.7584.03
S11 S - 2   ( B 2 ,   B 3 ,   B 4 ,   B 8 ,   B 12 ) ,   dNBR ,   dLST ,   A - 2 d σ H V 0 95.3791.4990.31
S12 S - 2   ( B 2 ,   B 3 ,   B 4 ,   B 8 ,   B 12 ) ,   dNBR ,   dLST ,   A - 2 d σ H H 0 95.4391.5890.39
S13L-8 (B2, B3, B4, B5, B7)92.8087.0886.93
S14L-8 (B2, B3, B4, B5, B7), NBR202092.8387.1387.12
S15L-8 (B2, B3, B4, B5, B7), dNBR93.3588.0087.43
S16L-8 (B2, B3, B4, B5, B7), dLST94.4189.8689.12
S17 L - 8   ( B 2 ,   B 3 ,   B 4 ,   B 5 ,   B 7 ) ,   A - 2 γ H V 94.6690.2189.23
S18 L - 8   ( B 2 ,   B 3 ,   B 4 ,   B 5 ,   B 7 ) ,   A - 2 γ H H 94.3089.5988.54
S19 L - 8   ( B 2 ,   B 3 ,   B 4 ,   B 5 ,   B 7 ) ,   A - 2 d σ H V 0 92.6186.7886.60
S20 L - 8   ( B 2 ,   B 3 ,   B 4 ,   B 5 ,   B 7 ) ,   A - 2 d σ H V 0 92.4286.4686.29
S21 L - 8   ( B 2 ,   B 3 ,   B 4 ,   B 5 ,   B 7 ) ,   S - 1 d σ V V 0 92.7787.0686.90
S22 L - 8   ( B 2 ,   B 3 ,   B 4 ,   B 5 ,   B 7 ) ,   S - 1 d σ V H 0 92.7787.0486.91
S23 L - 8   ( B 2 ,   B 3 ,   B 4 ,   B 5 ,   B 7 ) ,   dNBR ,   dLST ,   A - 2 d σ H V 0 95.9292.4791.61
S24 L - 8   ( B 2 ,   B 3 ,   B 4 ,   B 5 ,   B 7 ) ,   dNBR ,   dLST ,   A - 2 d σ H H 0 96.0092.6291.51
Table A2. Accuracy analysis of each classes in each scenarios. The best result of each class is written in bold.
Table A2. Accuracy analysis of each classes in each scenarios. The best result of each class is written in bold.
BurnedForestBarelandAgriculture
ScenariosUAPAF1UAPAF1UAPAF1UAPAF1
S185.7196.4690.7799.2394.0396.5685.1572.4278.2761.5485.9571.72
S286.1796.4691.0299.2393.8696.4787.3572.9479.4960.2387.6071.38
S392.8597.6495.1899.2494.8897.0186.2177.3281.5260.0086.7870.95
S490.3097.4993.7699.1695.2297.1585.1578.3581.6169.2389.2677.98
S587.5896.7691.9499.1796.3397.7388.6472.4279.7267.6486.3675.86
S687.2296.6191.6798.9596.7697.8490.7370.6279.4268.4086.7876.50
S786.1396.1790.8799.0593.5696.2387.7870.3678.1157.5788.0269.61
S884.2696.3189.8899.2393.7396.4086.6068.3076.3758.8286.7870.12
S985.5196.6190.7299.0593.4396.1687.4670.1077.8357.9787.1969.64
S1085.5196.6190.7299.,1594.0796.5486.5271.1378.0860.3585.5470.77
S1193.2397.4995.3199.3597.6598.4990.2681.1985.4875.1790.0881.95
S1293.0097.9495.4099.3597.7098.5291.7480.1585.5674.8390.9182.09
S1391.1597.2094.0899.1993.6496.3491.2282.9986.9158.6888.0270.41
S1491.1597.2094.0899.1493.6096.2992.4982.4787.1958.6589.6770.92
S1594.0297.3595.6599.0694.5096.7290.6582.4786.3759.2888.4370.98
S1693.7597.3595.5199.3495.7397.5090.3384.2887.2066.3689.6776.27
S1792.1196.4694.2499.1796.8998.0189.4783.2586.2571.8286.3678.42
S1892.5496.9094.6799.0496.4297.7190.4082.4786.2567.6585.5475.55
S1991.7696.9094.2699.3293.3996.2689.9282.7386.1757.3388.8469.69
S2090.9097.2093.9499.4593.1796.2191.0781.4485.9956.2589.2669.01
S2191.3997.0594.1399.3293.5296.3391.0183.5187.1057.9988.4370.05
S2290.9796.6193.7199.3293.6496.4090.4883.2586.7158.9088.8470.84
S2394.9597.0595.9999.0997.9998.5490.6184.5487.4778.8590.9184.45
S2494.7097.4996.0899.2798.1298.6991.1184.5487.7078.3489.6783.62

References

  1. Bonan, G.B. Forests and climate change: Forcings, feedbacks, and the climate benefits of forests. Science 2008, 320, 1444–1449. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Li, Q. Forest Bathing: How Trees Can Help You Find Health and Happiness; Viking Press: New York, NY, USA, 2018; ISBN 0525559868. [Google Scholar]
  3. WWF. Available online: https://www.worldwildlife.org/industries/responsible-forestry (accessed on 18 December 2021).
  4. Zhang, Q.; Ge, L.; Zhang, R.; Metternicht, G.I.; Du, Z.; Kuang, J.; Xu, M. Deep-learning-based burned area mapping using the synergy of Sentinel-1&2 data. Remote Sens. Environ. 2021, 264, 112575. [Google Scholar] [CrossRef]
  5. Schroeder, W.; Oliva, P.; Giglio, L.; Quayle, B.; Lorenz, E.; Morelli, F. Active fire detection using Landsat-8/OLI data. Remote Sens. Environ. 2016, 185, 210–220. [Google Scholar] [CrossRef] [Green Version]
  6. Palandjian, D.; Gitas, I.Z.; Wright, R. Burned area mapping and post-fire impact assessment in the Kassandra peninsula (Greece) using Landsat TM and Quickbird data. Geocarto Int. 2009, 24, 193–205. [Google Scholar] [CrossRef]
  7. Sertel, E.; Alganci, U. Comparison of pixel and object-based classification for burned area mapping using SPOT-6 images. Geomatics Nat. Hazards Risk 2016, 7, 1198–1206. [Google Scholar] [CrossRef] [Green Version]
  8. Widodo, J.; Sulaiman, A.; Awaluddin, A.; Riyadi, A.; Nasucha, M.; Perissin, D.; Sri Sumantyo, J.T. Application of SAR Interferometry Using ALOS-2 PALSAR-2 Data as Precise Method to Identify Degraded Peatland Areas Related to Forest Fire. Geosciences 2019, 9, 484. [Google Scholar] [CrossRef] [Green Version]
  9. Tanase, M.A.; Santoro, M.; Wegmüller, U.; de la Riva, J.; Pérez-Cabello, F. Properties of X-, C-and L-band repeat-pass interferometric SAR coherence in Mediterranean pine forests affected by fires. Remote Sens. Environ. 2010, 114, 2182–2194. [Google Scholar] [CrossRef]
  10. Philipp, M.B.; Levick, S.R. Exploring the potential of C-Band SAR in contributing to burn severity mapping in tropical savanna. Remote Sens. 2020, 12, 49. [Google Scholar] [CrossRef] [Green Version]
  11. Xulu, S.; Mbatha, N.; Peerbhay, K. Burned Area Mapping over the Southern Cape Forestry Region, South Africa Using Sentinel Data within GEE Cloud Platform. ISPRS Int. J. Geo-Inf. 2021, 10, 511. [Google Scholar] [CrossRef]
  12. Viedma, O.; Moreno, J.M.; Güngöroglu, C.; Cosgun, U.; Kavgacı, A. Recent land-use and land-cover changes and its driving factors in a fire-prone area of southwestern Turkey. J. Environ. Manag. 2017, 197, 719–731. [Google Scholar] [CrossRef]
  13. Elhag, M.; Yimaz, N.; Bahrawi, J.; Boteva, S. Evaluation of optical remote sensing data in burned areas mapping of Thasos Island, Greece. Earth Syst. Environ. 2020, 4, 813–826. [Google Scholar] [CrossRef]
  14. Thompson, M.P.; MacGregor, D.G.; Dunn, C.J.; Calkin, D.E.; Phipps, J. Rethinking the wildland fire management system. J. For. 2018, 116, 382–390. [Google Scholar] [CrossRef] [Green Version]
  15. Koutsias, N.; Karteris, M. Burned area mapping using logistic regression modeling of a single post-fire Landsat-5 Thematic Mapper image. Int. J. Remote Sens. 2000, 21, 673–687. [Google Scholar] [CrossRef]
  16. Pereira, A.A.; Libonati, R.; Rodrigues, J.A.; Nogueira, J.; Santos, F.L.M.; Oom, D.; Sanches, W.; Alvarado, S.T.; Pereira, J. Multi-sensor, active fire-supervised, one-class burned area mapping in the brazilian savanna. Remote Sens. 2021, 13, 4005. [Google Scholar] [CrossRef]
  17. Sali, M.; Piaser, E.; Boschetti, M.; Brivio, P.A.; Sona, G.; Bordogna, G.; Stroppiana, D. A Burned Area Mapping Algorithm for Sentinel-2 Data Based on Approximate Reasoning and Region Growing. Remote Sens. 2021, 13, 2214. [Google Scholar] [CrossRef]
  18. Matcı, D.K.; Avdan, U. Comparative analysis of unsupervised classification methods for mapping burned forest areas. Arab. J. Geosci. 2020, 13, 711. [Google Scholar] [CrossRef]
  19. Smiraglia, D.; Filipponi, F.; Mandrone, S.; Tornato, A.; Taramelli, A. Agreement index for burned area mapping: Integration of multiple spectral indices using Sentinel-2 satellite images. Remote Sens. 2020, 12, 1862. [Google Scholar] [CrossRef]
  20. Quintano, C.; Fernández-Manso, A.; Fernández-Manso, O.; Shimabukuro, Y.E. Mapping burned areas in Mediterranean countries using spectral mixture analysis from a uni-temporal perspective. Int. J. Remote Sens. 2006, 27, 645–662. [Google Scholar] [CrossRef]
  21. Kasischke, E.S.; French, N.H.F.; Harrell, P.; Christensen, N.L., Jr.; Ustin, S.L.; Barry, D. Monitoring of wildfires in boreal forests using large area AVHRR NDVI composite image data. Remote Sens. Environ. 1993, 45, 61–71. [Google Scholar] [CrossRef]
  22. Santana, N.C.; de Carvalho Júnior, O.A.; Gomes, R.A.T.; Guimarães, R.F. Burned-area detection in Amazonian environments using standardized time series per pixel in MODIS data. Remote Sens. 2018, 10, 1904. [Google Scholar] [CrossRef] [Green Version]
  23. Liu, S.; Wang, S.; Chi, T.; Wen, C.; Wu, T.; Wang, D. An improved combined vegetation difference index and burn scar index approach for mapping cropland burned areas using combined data from Landsat 8 multispectral and thermal infrared bands. Int. J. Wildl. Fire 2020, 29, 499–512. [Google Scholar] [CrossRef]
  24. Mashhadi, N.; Alganci, U. Determination of forest burn scar and burn severity from free satellite images: A comparative evaluation of spectral indices and machine learning classifiers. Int. J. Environ. Geoinform. 2021, 8, 488–497. [Google Scholar] [CrossRef]
  25. Smith, A.M.S.; Drake, N.A.; Wooster, M.J.; Hudak, A.T.; Holden, Z.A.; Gibbons, C.J. Production of Landsat ETM+ reference imagery of burned areas within Southern African savannahs: Comparison of methods and application to MODIS. Int. J. Remote Sens. 2007, 28, 2753–2775. [Google Scholar] [CrossRef]
  26. Martín, M.P.; Gómez, I.; Chuvieco, E. Burnt Area Index (BAIM) for burned area discrimination at regional scale using MODIS data. For. Ecol. Manag. 2006, 234, S221. [Google Scholar] [CrossRef]
  27. Pinty, B.; Verstraete, M.M. GEMI: A non-linear index to monitor global vegetation from satellites. Vegetatio 1992, 101, 15–20. [Google Scholar] [CrossRef]
  28. Trigg, S.; Flasse, S. An evaluation of different bi-spectral spaces for discriminating burned shrub-savannah. Int. J. Remote Sens. 2001, 22, 2641–2647. [Google Scholar] [CrossRef]
  29. Miller, J.D.; Thode, A.E. Quantifying burn severity in a heterogeneous landscape with a relative version of the delta Normalized Burn Ratio (dNBR). Remote Sens. Environ. 2007, 109, 66–80. [Google Scholar] [CrossRef]
  30. Mallinis, G.; Mitsopoulos, I.; Chrysafi, I. Evaluating and comparing sentinel 2A and landsat-8 operational land imager (OLI) spectral indices for estimating fire severity in a Mediterranean pine ecosystem of Greece. GISci. Remote Sens. 2018, 55, 1–18. [Google Scholar] [CrossRef]
  31. Achour, H.; Toujani, A.; Trabelsi, H.; Jaouadi, W. Evaluation and comparison of Sentinel-2 MSI, Landsat 8 OLI, and EFFIS data for forest fires mapping. Illustrations from the summer 2017 fires in Tunisia. Geocarto Int. 2021, 36, 1–20. [Google Scholar] [CrossRef]
  32. Ban, Y.; Zhang, P.; Nascetti, A.; Bevington, A.R.; Wulder, M.A. Near real-time wildfire progression monitoring with Sentinel-1 SAR time series and deep learning. Sci. Rep. 2020, 10, 1322. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Stroppiana, D.; Azar, R.; Calò, F.; Pepe, A.; Imperatore, P.; Boschetti, M.; Silva, J.; Brivio, P.A.; Lanari, R. Integration of optical and SAR data for burned area mapping in Mediterranean Regions. Remote Sens. 2015, 7, 1320–1345. [Google Scholar] [CrossRef] [Green Version]
  34. Filipponi, F. Exploitation of sentinel-2 time series to map burned areas at the national level: A case study on the 2017 italy wildfires. Remote Sens. 2019, 11, 622. [Google Scholar] [CrossRef] [Green Version]
  35. Chuvieco, E.; Lizundia-Loiola, J.; Pettinari, M.L.; Ramo, R.; Padilla, M.; Tansey, K.; Mouillot, F.; Laurent, P.; Storm, T.; Heil, A. Generation and analysis of a new global burned area product based on MODIS 250 m reflectance bands and thermal anomalies. Earth Syst. Sci. Data 2018, 10, 2015–2031. [Google Scholar] [CrossRef] [Green Version]
  36. Giglio, L.; Boschetti, L.; Roy, D.P.; Humber, M.L.; Justice, C.O. The Collection 6 MODIS burned area mapping algorithm and product. Remote Sens. Environ. 2018, 217, 72–85. [Google Scholar] [CrossRef]
  37. Fernández-Manso, A.; Quintano, C. A synergetic approach to burned area mapping using maximum entropy modeling trained with hyperspectral data and VIIRS hotspots. Remote Sens. 2020, 12, 858. [Google Scholar] [CrossRef] [Green Version]
  38. Otón, G.; Ramo, R.; Lizundia-Loiola, J.; Chuvieco, E. Global detection of long-term (1982–2017) burned area with AVHRR-LTDR data. Remote Sens. 2019, 11, 2079. [Google Scholar] [CrossRef] [Green Version]
  39. Roy, D.P.; Huang, H.; Boschetti, L.; Giglio, L.; Yan, L.; Zhang, H.H.; Li, Z. Landsat-8 and Sentinel-2 burned area mapping-A combined sensor multi-temporal change detection approach. Remote Sens. Environ. 2019, 231, 111254. [Google Scholar] [CrossRef]
  40. Ngadze, F.; Mpakairi, K.S.; Kavhu, B.; Ndaimani, H.; Maremba, M.S. Exploring the utility of Sentinel-2 MSI and Landsat 8 OLI in burned area mapping for a heterogenous savannah landscape. PLoS ONE 2020, 15, e0232962. [Google Scholar] [CrossRef]
  41. Moreira, A.; Prats-Iraola, P.; Younis, M.; Krieger, G.; Hajnsek, I.; Papathanassiou, K.P. A tutorial on synthetic aperture radar. IEEE Geosci. Remote Sens. Mag. 2013, 1, 6–43. [Google Scholar] [CrossRef] [Green Version]
  42. Tariq, A.; Shu, H.; Li, Q.; Altan, O.; Khan, M.R.; Baqa, M.F.; Lu, L. Quantitative Analysis of Forest Fires in Southeastern Australia Using SAR Data. Remote Sens. 2021, 13, 2386. [Google Scholar] [CrossRef]
  43. Donezar, U.; De Blas, T.; Larrañaga, A.; Ros, F.; Albizua, L.; Steel, A.; Broglia, M. Applicability of the multitemporal coherence approach to sentinel-1 for the detection and delineation of burnt areas in the context of the copernicus emergency management service. Remote Sens. 2019, 11, 2607. [Google Scholar] [CrossRef] [Green Version]
  44. Verhegghen, A.; Eva, H.; Ceccherini, G.; Achard, F.; Gond, V.; Gourlet-Fleury, S.; Cerutti, P.O. The potential of Sentinel satellites for burnt area mapping and monitoring in the Congo Basin forests. Remote Sens. 2016, 8, 986. [Google Scholar] [CrossRef] [Green Version]
  45. Belenguer-Plomer, M.A.; Tanase, M.A.; Chuvieco, E.; Bovolo, F. CNN-based burned area mapping using radar and optical data. Remote Sens. Environ. 2021, 260, 112468. [Google Scholar] [CrossRef]
  46. General Directorate of Forestry, OGM 2020. ISBN 978-605-7599-68-1. Available online: https://www.ogm.gov.tr/tr/ormanlarimiz/Turkiye-Orman-Varligi (accessed on 28 December 2021). (In Turkish)
  47. Mauri, A.; Di Leo, M.; de Rigo, D.; Caudullo, G. Pinus Halepensis and Pinus Brutia in Europe: Distribution, Habitat, Usage and Threats; San-Miguel-Ayanz, J., de Rigo, D., Caudullo, G., Houston Durrant, T., Mauri, A., Eds.; European Atlas of Forest Tree Species; Publications Office of the EU: Luxembourg, 2016; p. e0166b8.
  48. Kaynaş, B.Y.; Öney Birol, S. Afforestation of post-fire Pinus brutia Ten forests Effects of Eucalyptus camaldulensis Dehnh P brutia and Pinus pinea L leaf extracts on cell division in P brutia seeds. Allelopathy J. 2020, 49, 125–140. [Google Scholar] [CrossRef]
  49. Chambel, M.R.; Climent, J.; Pichot, C.; Ducci, F. Chapter 5: Mediterranean pines (Pinus halepensis Mill. and brutia Ten.). In Forest Tree Breeding in Europe; Managing Forest Ecosystems 25; Springer: Dordrecht, The Netherlands, 2013; pp. 229–265. [Google Scholar] [CrossRef]
  50. Dailysabah. Available online: https://www.dailysabah.com/opinion/op-ed/fires-floods-mucilage-whats-happening-in-turkey (accessed on 20 December 2021).
  51. San-Miguel-Ayanz, J.; Durrant, T.; Boca, R.; Maianti, P.; Liberta, G.; Artes-Vivancos, T.; Oom, D.; Branco, A.; de Rigo, D.; Ferrari, D.; et al. Forest Fires in Europe, Middle East and North Africa 2019, EUR 30402 EN; Publications Office of the European Union: Luxembourg, 2020; ISBN 978-92-76-23209-4. [CrossRef]
  52. San-Miguel-Ayanz, J.; Durrant, T.; Boca, R.; Maianti, P.; Liberta, G.; Artes-Vivancos, T.; Oom, D.; Branco, A.; de Rigo, D.; Ferrari, D. Forest Fires in Europe, Middle East and North Africa 2020, EUR 30862 EN; Publications Office of the European Union: Luxembourg, 2021; ISBN 978-92-76-42351-5. [CrossRef]
  53. Reuters. Available online: https://www.reuters.com/world/middle-east/fire-near-turkish-power-plant-under-control-local-mayor-2021-08-04/ (accessed on 20 December 2021).
  54. Atak, B.K.; Tonyaloğlu, E.E. Evaluating spectral indices for estimating burned areas in the case of Izmir/Turkey. Eurasian J. For. Sci. 2020, 8, 49–59. [Google Scholar] [CrossRef]
  55. Veraverbeke, S.; Lhermitte, S.; Verstraeten, W.W.; Goossens, R. Evaluation of pre/post-fire differenced spectral indices for assessing burn severity in a Mediterranean environment with Landsat Thematic Mapper. Int. J. Remote Sens. 2011, 32, 3521–3537. [Google Scholar] [CrossRef] [Green Version]
  56. Miller, J.D.; Knapp, E.E.; Key, C.H.; Skinner, C.N.; Isbell, C.J.; Creasy, R.M.; Sherlock, J.W. Calibration and validation of the relative differenced Normalized Burn Ratio (RdNBR) to three measures of fire severity in the Sierra Nevada and Klamath Mountains, California, USA. Remote Sens. Environ. 2009, 113, 645–656. [Google Scholar] [CrossRef]
  57. Veraverbeke, S.; Lhermitte, S.; Verstraeten, W.W.; Goossens, R. The temporal dimension of differenced Normalized Burn Ratio (dNBR) fire/burn severity studies: The case of the large 2007 Peloponnese wildfires in Greece. Remote Sens. Environ. 2010, 114, 2548–2563. [Google Scholar] [CrossRef] [Green Version]
  58. Key, C.H.; Benson, N.C. Landscape Assessment: Remote Sensing of Severity, the Normalized Burn Ratio and Ground Measure of Severity, the Composite Burn Index; FIREMON: Fire Effects Monitoring and Inventory System, General Technical Report, RMRS-GTR-164-CD:LA1-LA51; USDA Forest Service, Rocky Mountain Research Station: Ogden, UT, USA, 2005.
  59. Key, C.H.; Benson, N.C. Measuring and remote sensing of burn severity. In Proceedings of the Joint Fire Science Conference and Workshop; University of Idaho and International Association of Wildland Fire: Boise, ID, USA, 1999; Volume 2, p. 284. [Google Scholar]
  60. Qin, Z.; Karnieli, A.; Berliner, P. A mono-window algorithm for retrieving land surface temperature from Landsat TM data and its application to the Israel-Egypt border region. Int. J. Remote Sens. 2001, 22, 3719–3746. [Google Scholar] [CrossRef]
  61. Sobrino, J.A.; Jiménez-Muñoz, J.C.; Sòria, G.; Romaguera, M.; Guanter, L.; Moreno, J.; Plaza, A.; Martínez, P. Land surface emissivity retrieval from different VNIR and TIR sensors. IEEE Trans. Geosci. Remote Sens. 2008, 46, 316–327. [Google Scholar] [CrossRef]
  62. Sekertekin, A.; Bonafoni, S. Land surface temperature retrieval from Landsat 5, 7, and 8 over rural areas: Assessment of different retrieval algorithms and emissivity models and toolbox implementation. Remote Sens. 2020, 12, 294. [Google Scholar] [CrossRef] [Green Version]
  63. Hagberg, J.O.; Ulander, L.M.H.; Askne, J. Repeat-pass SAR interferometry over forested terrain. IEEE Trans. Geosci. Remote Sens. 1995, 33, 331–340. [Google Scholar] [CrossRef]
  64. Abdikan, S.; Bayik, C.; Ustuner, M.; Sanli, F.B. Repeat-pass interferometric and backscatter analysis of X-band PAZ satellite–first results. Int. Arch. Photogramm. Remote Sens. Spat. Inf. Sci. 2020, 43, 253–258. [Google Scholar] [CrossRef]
  65. SNAP v8.0. Available online: http://step.esa.int (accessed on 14 December 2021).
  66. Shao, Z.; Sumari, N.S.; Portnov, A.; Ujoh, F.; Musakwa, W.; Mandela, P.J. Urban sprawl and its impact on sustainable urban development: A combination of remote sensing and social media data. Geo-Spat. Inf. Sci. 2021, 24, 241–255. [Google Scholar] [CrossRef]
  67. Breiman, L. Random forests. Mach. Learn. 2001, 45, 5–32. [Google Scholar] [CrossRef] [Green Version]
  68. Ok, A.O.; Akar, O.; Gungor, O. Evaluation of random forest method for agricultural crop classification. Eur. J. Remote Sens. 2012, 45, 421–432. [Google Scholar] [CrossRef]
  69. Van der Linden, S.; Rabe, A.; Held, M.; Jakimow, B.; Leitão, P.J.; Okujeni, A.; Schwieder, M.; Suess, S.; Hostert, P. The EnMAP-Box—A toolbox and application programming interface for EnMAP data processing. Remote Sens. 2015, 7, 11249–11266. [Google Scholar] [CrossRef] [Green Version]
  70. EnMAP-Box Developers. EnMAP-Box 3—A QGIS Plugin to Process and Visualize Hyperspectral Remote Sensing Data. 2019. Available online: https://enmap-box.readthedocs.io (accessed on 8 December 2021).
  71. QGIS. Available online: http://www.qgis.org (accessed on 12 December 2021).
  72. Foody, G.M. Status of land cover classification accuracy assessment. Remote Sens. Environ. 2002, 80, 185–201. [Google Scholar] [CrossRef]
  73. Ustuner, M.; Balik Sanli, F. Polarimetric Target Decompositions and Light Gradient Boosting Machine for Crop Classification: A Comparative Evaluation. ISPRS Int. J. Geo-Inf. 2019, 8, 97. [Google Scholar] [CrossRef] [Green Version]
  74. Tadono, T.; Ishida, H.; Oda, F.; Naito, S.; Minakawa, K.; Iwamoto, H. Precise global DEM generation by ALOS PRISM. ISPRS Ann. Photogramm. Remote Sens. Spat. Inf. Sci. 2014, 2, 71–76. [Google Scholar] [CrossRef] [Green Version]
  75. Quintano, C.; Fernández-Manso, A.; Fernández-Manso, O. Combination of Landsat and Sentinel-2 MSI data for initial assessing of burn severity. Int. J. Appl. Earth Obs. Geoinf. 2018, 64, 221–225. [Google Scholar] [CrossRef]
  76. Nasery, S.; Kalkan, K. Burn area detection and burn severity assessment using Sentinel 2 MSI data: The case of Karabağlar district, İzmir/Turkey. Turk. J. Geosci. 2020, 1, 72–77. [Google Scholar]
Figure 1. The location of the forest fire. BA is seen in the black patch on the Sentinel-2 image (NIR-G-B) presented in the map on the right.
Figure 1. The location of the forest fire. BA is seen in the black patch on the Sentinel-2 image (NIR-G-B) presented in the map on the right.
Forests 13 00347 g001
Figure 2. Workflow of the study.
Figure 2. Workflow of the study.
Forests 13 00347 g002
Figure 3. Accuracies of the scenarios (%) (the details are presented in Appendix A Table A1).
Figure 3. Accuracies of the scenarios (%) (the details are presented in Appendix A Table A1).
Forests 13 00347 g003
Figure 4. User accuracy (UA) and producer accuracy (PA) of each class (the details are presented in Appendix A Table A2).
Figure 4. User accuracy (UA) and producer accuracy (PA) of each class (the details are presented in Appendix A Table A2).
Forests 13 00347 g004
Figure 5. The highest classification accuracy results of scenarios using Sentinel-2 (S12) and Landsat-8 (S24).
Figure 5. The highest classification accuracy results of scenarios using Sentinel-2 (S12) and Landsat-8 (S24).
Forests 13 00347 g005
Figure 6. Burn severity maps derived from Sentinel-2 and Landsat-8.
Figure 6. Burn severity maps derived from Sentinel-2 and Landsat-8.
Forests 13 00347 g006
Figure 7. Spatial analysis of burn severity over slope.
Figure 7. Spatial analysis of burn severity over slope.
Forests 13 00347 g007
Table 1. Multi-sensor remote sensing satellite data specifications used in this study.
Table 1. Multi-sensor remote sensing satellite data specifications used in this study.
Data TypeSatellite MissionAcquisition
Date(Pre–Post)
Image Channel (Band or Polarization)Central Wavelength or FrequencySpatial Resolution (m)
OpticalSentinel-24 March 2019
13 March 2020
Blue (B2)492.4 nm10
Green (B3)559.8 nm10
Red (B4)664.6 nm10
NIR (B8)832.8 nm10
Narrow NIR (B8A)864.7 nm20
SWIR (B12)2202.4 nm20
Landsat-816 March 2019
2 March 2020
Blue (B2)482 nm30
Green (B3)561 nm30
Red (B4)655 nm30
NIR (B5)865 nm30
SWIR 2 (B7)2200 nm30
ThermalLandsat-816 March 2019
2 March 2020
TIR-1 (B10)10,800 nm100 (Resampled to 30)
SARSentinel-110 March 2019
4 March 2020
VV (C band)5.405 GHz2.33 × 13.91 (r × az)
VH (C band)5.405 GHz2.33 × 13.91 (r × az)
ALOS-210 March 2019
8 March 2020
HH (L band)1.2 GHz4.29 × 3.41 (r × az)
HV (L band)1.2 GHz4.29 × 3.41 (r × az)
Table 2. Burn severity spectral indices.
Table 2. Burn severity spectral indices.
Spectral IndexLandsat-8 OLI EquationSentinel-2 MSI EquationReference
NBR(B5 − B7)/(B5 + B7)(B8A − B12)/(B8A + B12)Key and Benson [50]
dNBR(NBRpre-NBRpost)(NBRpre-NBRpost)Miller and Thode [49]
Table 3. Burn severity levels obtained by calculating dNBR, proposed by [59].
Table 3. Burn severity levels obtained by calculating dNBR, proposed by [59].
Severity LeveldNBR Range (Not Scaled)
Unburned<−0.1
Low Severity0.1–0.26
Moderate Low Severity0.27–0.43
Moderate High Severity0.44–0.65
High Severity>0.66
Table 4. The scenarios created for the comparison. S-2 = post-fire Sentinel-2 data, L-8 = Landsat-8, A-2 = ALOS-2, S-1 = Sentinel-1, γ   = coherence, d σ 0 = backscatter difference.
Table 4. The scenarios created for the comparison. S-2 = post-fire Sentinel-2 data, L-8 = Landsat-8, A-2 = ALOS-2, S-1 = Sentinel-1, γ   = coherence, d σ 0 = backscatter difference.
ScenariosFeaturesScenariosFeatures
S1S-2 (B2, B3, B4, B8, B12)S13L-8 (B2, B3, B4, B5, B7)
S2S-2 (B2, B3, B4, B8, B12), S-2 NBR2020S14L-8 (B2, B3, B4, B5, B7), L-8 NBR2020
S3S-2 (B2, B3, B4, B8, B12), dNBRS15L-8 (B2, B3, B4, B5, B7), dNBR
S4S-2 (B2, B3, B4, B8, B12), dLSTS16L-8 (B2, B3, B4, B5, B7), dLST
S5 S - 2   ( B 2 ,   B 3 ,   B 4 ,   B 8 ,   B 12 ) ,   A - 2 γ H V S17 L - 8   ( B 2 ,   B 3 ,   B 4 ,   B 5 ,   B 7 ) ,   A - 2 γ H V
S6 S - 2   ( B 2 ,   B 3 ,   B 4 ,   B 8 ,   B 12 ) ,   A - 2 γ H H S18 L - 8   ( B 2 ,   B 3 ,   B 4 ,   B 5 ,   B 7 ) ,   A - 2 γ H H
S7 S - 2   ( B 2 ,   B 3 ,   B 4 ,   B 8 ,   B 12 ) ,   A - 2 d σ H V 0 S19 L - 8   ( B 2 ,   B 3 ,   B 4 ,   B 5 ,   B 7 ) ,   A - 2 d σ H V 0
S8 S - 2   ( B 2 ,   B 3 ,   B 4 ,   B 8 ,   B 12 ) ,   A - 2 d σ H V 0 S20 L - 8   ( B 2 ,   B 3 ,   B 4 ,   B 5 ,   B 7 ) ,   A - 2 d σ H V 0
S9 S - 2   ( B 2 ,   B 3 ,   B 4 ,   B 8 ,   B 12 ) ,   S - 1 d σ V V 0 S21 L - 8   ( B 2 ,   B 3 ,   B 4 ,   B 5 ,   B 7 ) ,   S - 1 d σ V V 0
S10 S - 2   ( B 2 ,   B 3 ,   B 4 ,   B 8 ,   B 12 ) ,   S - 1 d σ V H 0 S22 L - 8   ( B 2 ,   B 3 ,   B 4 ,   B 5 ,   B 7 ) ,   S - 1 d σ V H 0
S11 S - 2   ( B 2 ,   B 3 ,   B 4 ,   B 8 ,   B 12 ) ,   dNBR ,   dLST ,   A - 2 d σ H V 0 S23 L - 8   ( B 2 ,   B 3 ,   B 4 ,   B 5 ,   B 7 ) ,   dNBR ,   dLST ,   A - 2 d σ H V 0
S12 S - 2   ( B 2 ,   B 3 ,   B 4 ,   B 8 ,   B 12 ) ,   dNBR ,   dLST ,   A - 2 d σ H H 0 S24 L - 8   ( B 2 ,   B 3 ,   B 4 ,   B 5 ,   B 7 ) ,   dNBR ,   dLST ,   A - 2 d σ H H 0
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Abdikan, S.; Bayik, C.; Sekertekin, A.; Bektas Balcik, F.; Karimzadeh, S.; Matsuoka, M.; Balik Sanli, F. Burned Area Detection Using Multi-Sensor SAR, Optical, and Thermal Data in Mediterranean Pine Forest. Forests 2022, 13, 347. https://0-doi-org.brum.beds.ac.uk/10.3390/f13020347

AMA Style

Abdikan S, Bayik C, Sekertekin A, Bektas Balcik F, Karimzadeh S, Matsuoka M, Balik Sanli F. Burned Area Detection Using Multi-Sensor SAR, Optical, and Thermal Data in Mediterranean Pine Forest. Forests. 2022; 13(2):347. https://0-doi-org.brum.beds.ac.uk/10.3390/f13020347

Chicago/Turabian Style

Abdikan, Saygin, Caglar Bayik, Aliihsan Sekertekin, Filiz Bektas Balcik, Sadra Karimzadeh, Masashi Matsuoka, and Fusun Balik Sanli. 2022. "Burned Area Detection Using Multi-Sensor SAR, Optical, and Thermal Data in Mediterranean Pine Forest" Forests 13, no. 2: 347. https://0-doi-org.brum.beds.ac.uk/10.3390/f13020347

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