Burned Area Mapping Using Multi-Temporal Sentinel-2 Data by Applying the Relative Differenced Aerosol-Free Vegetation Index (RdAFRI)

Assessing the development of wildfire scars during a period of consecutive active fires and smoke overcast is a challenge. The study was conducted during nine months when Israel experienced massive pyro-terrorism attacks of more than 1100 fires from the Gaza Strip. The current project strives at developing and using an advanced Earth observation approach for accurate post-fire spatial and temporal assessment shortly after the event ends while eliminating the influence of biomass burning smoke on the ground signal. For fulfilling this goal, the Aerosol-Free Vegetation Index (AFRI), which has a meaningful advantage in penetrating an opaque atmosphere influenced by biomass burning smoke, was used. On top of it, under clear sky conditions, the AFRI closely resembles the widely used Normalized Difference Vegetation Index (NDVI), and it retains the same level of index values under smoke. The relative differenced AFRI (RdAFRI) set of algorithms was implemented at the same procedure commonly used with the Relative differenced Normalized Burn Ratio (RdBRN). The algorithm was applied to 24 Sentinel-2 Level-2A images throughout the study period. While validating with ground observations, the RdAFRI-based algorithms produced an overall accuracy of 90%. Furthermore, the RdAFRI maps were smoother than the equivalent RdNBR, with noise levels two orders of magnitude lower than the latter. Consequently, applying the RdAFRI, it is possible to distinguish among four severity categories. However, due to different cloud cover on the two consecutive dates, an automatic determination of a threshold level was not possible. Therefore, two threshold levels were considered through visual inspection and manually assigned to each imaging date. The novel procedure enables calculating the spatio-temporal dynamics of the fire scars along with the statistics of the burned vegetation species within the study area.


Introduction
In recent years, wildfires have become serious human and environmental concerns for several reasons. Most importantly, they pose a threat to human life, flora, and fauna, as well as properties and economic losses. Smoke, which consists of a mixture of fine particulate matter and gases, can cause health problems (e.g., [1]). Wildfires are considered as one of the worst ecological disturbances to long-term records of vegetation phenology, since land-cover alterations are the basis for the understanding of the biological responses to climate change, such as impacting carbon emissions, at regional to continental where ρ is the reflectance value of the indicated spectral band-red, NIR, and SWIR2 (around 2.1 µm). The different approach behind these two indices is raised since the red and the SWIR bands have a different response to the wildfire consequences [27,28]. The NDVI relies on the chlorophyll content in the red band and therefore indicates the state of the live vegetation in contrast to bare surfaces, e.g., a burned area. At the same time, the NBR is strongly related to water content in soils and vegetation as characterized by the SWIR band [29]. Several studies concluded that the performance of SWIR-derived indices is better than those who build upon the red and NIR bands (e.g., [23]). NBR has recently turned to be a standard method for assessing burn severity using different Earth observation satellites, e.g., Landsat [30][31][32], MODIS [24,33], ASTER, [24,34], and Sentinel-2 [10]. Based on the NBR, the multi-date change detection between the pre-(t0) and the post-burned area (t1) is formulated as the difference (delta, d, dNBR) between these areas, based on the NBR [26]: Additionally, an improvement of dNBR accuracy and removing heteroscedasticity from the index distribution in terms of a relative scale was named Relative dNBR (RdNBR) [29,35]: Near real time of consecutive fires might raise another difficulty, since the biomass burning smoke often contaminates the atmosphere. The NBR, involving the SWIR band, can penetrate an opaque atmospheric column influenced by biomass burning smoke, without the need for explicit correction for the aerosol effect [36,37].
Similar to the NBR, the Aerosol-Free Vegetation Index (AFRI) is also based on the correlation between the visible-red and the SWIR2 band [15,38]: Unlike NBR, the SWIR2 band is multiplied by 0.5. This empirical constant was proven to correlate the SWIR2 and the red bands in detecting vegetation under clear sky conditions [15]. Therefore, on top of the AFRI's main advantage in penetrating a contaminated atmosphere by smoke and smog, it closely resembles the NDVI when no smoke exists. Since the values of both indices are almost identical, it is possible to assess vegetation conditions without using the visible bands that are more sensitive to atmospheric scattering.
Wildfire can be either nature-or human-induced [39]. Worldwide, the majority of the fires have been set by people by traditional land clearing for hunting and shifting cultivation [40]. Wildfires also exist in well-developed countries. For example, in the United States, about 85% of wildland fires are caused by human activities, such as campfires left unattended, debris burning, equipment use and malfunctions, discarded cigarettes, and intentional acts of arson [41,42]. Another type of fire arsons is known as pyro-terrorism: action intentionally to cause damage to property, burn crops, destroy infrastructures, and frighten the civilian population [43].
The Kit Terror is a name referring to Palestinian protest activity that, among other actions, flew an enormous number of kites and balloons, carrying Molotov cocktails across the Israel-Gaza Strip perimeter fence to set fire in the Israeli side of the border. The kite and balloon attacks, which had begun on 11 April 2018, continued non-stop, day by day, during the Mediterranean dry seasons-spring, summer, and fall, 2018-and faded at the beginning of the rainy season. As a result, more than 1100 fires were started, destroying tens of square kilometers of agricultural fields, nature reserves, and forests. The Kit Terror was the driver to establish a continuous monitoring program for counting and mapping the burned scars in the areas, which were mostly forests and woodlands owned and maintained by the Jewish National Fund (JNF).
Assessing the development of fire scars along with reliable landscape observation, during a period of consecutive active fires and smoke overcast, is a challenge. Therefore, the current project strives at developing and using an advanced approach for accurate post-fire assessment while eliminating the influence of biomass burning smoke on the ground signal. To fulfill this goal, the study took advantage of the relatively high revisit time of the S2 constellation and the ability of the AFRI to penetrate the smoky atmosphere. The study presents an alternative index to the NDVI and NBR for detecting fire scars and quantify burn severity. The landscape is characterized by low hilly terrain. The typical soils are either alluvial, loessial, or sandy soils. The area contains diverse types of planted and natural trees and shrubs, as listed in Table 1. Although the study vegetation is termed 'forests', the trees are rather sparse, savanna-like, spread within understory of grasses and shrubs. All species were mapped by the JNF and are stored in a geographic information system (GIS) database.  The landscape is characterized by low hilly terrain. The typical soils are either alluvial, loessial, or sandy soils. The area contains diverse types of planted and natural trees and shrubs, as listed in Table 1. Although the study vegetation is termed 'forests', the trees are rather sparse, savanna-like, spread within understory of grasses and shrubs. All species were mapped by the JNF and are stored in a geographic information system (GIS) database.

Sentinel-2 Data Collection
The S2 constellation is an Earth observation mission that acquires optical images over land and coastal waters. There are five S2 data products available on the Copernicus Open Access Hub [44]: Level-0, Level-1A, Level-1B, Level-1C, and Level-2A. The Level-0 products correspond to raw images still on board compressed. The Level-1A products are raw images after decompression. The Level-1B are radiometrically corrected products, whereas the Level-1C are orthorectified images providing top of atmosphere reflectance. Finally, the Level-2A are orthorectified products providing surface reflectance and basic pixel classification (including classes for different types of clouds) [45]. In this study, 25 Level-2A products were selected with cloud coverage inferior to 15%, from 6 April 2018 to 22 December 2018, in order to monitor the study area during the kite and balloon attacks period. The NIR band (B8) and the SWIR2 (B12), at 10 and 20 m spatial resolution, respectively, were used for calculating the relevant spectral indices.

Ground Observations
Throughout the period of the Kit Terror, from 11 April to 22 December 2018, daily ground inspections were performed by the district JNF's foresters, rangers, and firefighters. The affected areas were accurately mapped immediately after the active fire was extinguished using a designated mobile application. As a result, a two-class (burned/unburned) digital GIS database of fire scars was produced ( Figure 2). This database was used as ground truth for validating the Earth observation's algorithms.
Remote Sens. 2020, 12, x 6 of 21 The S2 Level-2A product contains vector cloud and cirrus masks that are created as a product of the atmospheric correction. These masks are used in this proposed methodology in order to compute the relative dAFRI only at the cloud-free pixels. The flowchart of the proposed classification algorithm is detailed in Figure 3. The indicator , ( ) represents the pixel located at row and column in the i-th image. The AFRI of the i-th image is computed from Equation (5), only for the cloud-free pixels, where ρ and ρ are the surface reflectance available at bands 8 and 12, respectively. Band 12 is resampled to 10 m in order to preserve the spatial resolution of band 8. The RdAFRI(i − 1, i) is computed only if the pixel ( , ) of the (i − 1)-th image is also cloud-free. If not, the relative differenced AFRI is derived by the i-th and (i − 2)-th images, i.e., RdAFRI = RdAFRI(i − 2, i). Finally, the pixel ( , ) is considered a burned surface if RdAFRI is greater than a fixed threshold. In this study, two thresholds

AFRI-Based Algorithms Development
Similar to the procedures developed for the NBR-based set of algorithms, the most intuitive burned area-mapping indicator consists of an absolute change detection methodology obtained subtracting a post-fire AFRI image from a pre-fire AFRI image to derive the differenced AFRI (dAFRI): Then, the dAFRI, for two successive images collected in this study, can be formulated as: where (i) indicates the i-th image in the database. The dAFRI(i − 1, i) can present problems in the cases with low vegetation values for the image taken at (i − 1): the absolute change will be small, and the index will not be able to detect the burned area. In order to avoid this issue, similarly to the work proposed in [29,35], the relative differenced AFRI (RdAFRI) was defined as: Positive RdAFRI(i − 1, i) values represent a decrease in vegetation cover, while negative values represent an increase in vegetation cover. Therefore, the burned area can be easily identified when the relativized delta AFRI is greater than a threshold that has been chosen a priori.
The S2 Level-2A product contains vector cloud and cirrus masks that are created as a product of the atmospheric correction. These masks are used in this proposed methodology in order to compute the relative dAFRI only at the cloud-free pixels. The flowchart of the proposed classification algorithm is detailed in Figure 3. The indicator p j,k (i) represents the pixel located at row j and column k in the i-th image. The AFRI of the i-th image is computed from Equation (5), only for the cloud-free pixels, where ρ NIR and ρ SWIR2 are the surface reflectance available at bands 8 and 12, respectively. Band 12 is resampled to 10 m in order to preserve the spatial resolution of band 8. The RdAFRI(i − 1, i) is computed only if the pixel ( j, k) of the (i − 1)-th image is also cloud-free. If not, the relative differenced AFRI is derived by the i-th and (i − 2)-th images, i.e., RdAFRI = RdAFRI(i − 2, i). Finally, the pixel ( j, k) is considered a burned surface if RdAFRI is greater than a fixed threshold. In this study, two thresholds were used: THR 1 for the relative differenced AFRI obtained by the i-th and (i − 1)-th images, and THR 2 for the relative differenced AFRI obtained by the i-th and (i − 2)-th images. With no loss of importance, the proposed methodology can be extended in order to manage more multi-temporal images, i.e., p j,k (i − 3), p j,k (i − 4), etc. In order to avoid that seasonal vegetation changes may be falsely detected as a burned area, a limited time expansion should be considered. The choice of the thresholds THR 1 and THR 2 is crucial in order to label a pixel as a burned or unburned area. An automated method, already used for wildfire burning and severity detection using Sentinel-2 images [21], is the Otsu threshold. This methodology distinguishes between background and foreground in imagery by creating two classes with minimal intraclass variance [46]. Alternatively, the value of the thresholds can also be chosen manually. In the current study, results were presented using both Otsu and manual thresholds.
Since the proposed classification algorithm is based on the RdAFRI values compared to a threshold, it is an applicative practice to quantify the noise degree of the relative difference AFRI images. For instance, the level of noise, denoted hereafter as a noise index (NI) of an image, can be computed from the following equation: where σ 2 j,k is the variance of the pixel ( j, k) computed in a neighborhood as a 3 × 3 matrix, N is the total number of pixels, and m and n are the number of rows and columns, respectively. THR2 is crucial in order to label a pixel as a burned or unburned area. An automated method, already used for wildfire burning and severity detection using Sentinel-2 images [21], is the Otsu threshold. This methodology distinguishes between background and foreground in imagery by creating two classes with minimal intraclass variance [46]. Alternatively, the value of the thresholds can also be chosen manually. In the current study, results were presented using both Otsu and manual thresholds. Since the proposed classification algorithm is based on the RdAFRI values compared to a threshold, it is an applicative practice to quantify the noise degree of the relative difference AFRI images. For instance, the level of noise, denoted hereafter as a noise index (NI) of an image, can be computed from the following equation: where , is the variance of the pixel ( , ) computed in a neighborhood as a 3 × 3 matrix, is the total number of pixels, and and are the number of rows and columns, respectively.

Results
Throughout the study area, more than 1100 fires were counted by the JNF's foresters, rangers, and firefighters. Ground observations revealed that most of the wildfires could be classified as 'surface fire', consuming mostly the understory grasses, shrubs, litter, and woody material lying on the ground [47]. Occasionally, the surface fire reached and injured the trunks of mature trees. In eucalyptus plots, damages were also observed in the lower foliage ( Figure 4A). Pinus halepensis (Aleppo pine), the dominant pine species in the study area, is very sensitive to heat due to its relative thin bark and high amount of resin [48,49]. Therefore, many pines did not survive the surface fire ( Figure 4B). Only the minority of the fires were characterized as 'canopy fires'. In other species, such as acacia, mesquites, and tamarisk, the foliage was damaged, but the trees were recovered soon after the fire.
For comparing the performance of the NBR, NDVI, and AFRI indices, a section of the S2 Level-2A image obtained on 10 July 2018, close to the Gaza Strip border, is presented in Figure 5. The transects represent a common situation when smoke, at different intensities, covers a variety of substrates-cultivated, bare soils, fire scars, and more. The true-color composite image (RGB = 0.665, 0.56, 0.49 µm) shows the open fire (light-orange hue), the biomass burning smoke (white hue), as well as burned scars that are a few days old (dark surfaces). The three indices were produced at 10-m spatial resolution, according to Equations (1), (2), and (5), from the surface reflectance values along Remote Sens. 2020, 12, 2753 9 of 22 a cross-section of 2771 m ( Figure 6, line A-A). This line was selected since it passes cultivated fields, bare soil, and was overcast by light smoke that characterized the entire region along the whole study period. This line is subdivided into several segments. From pixel 0 to 50 (and similarly from pixel 180 to 190) over the agricultural field where no smoke exists, the AFRI values accurately mimic those of the NDVI, but the NBR values are significantly lower. From pixel 50 to 180, over the bare soil, AFRI values are somehow higher than those of the NDVI. However, the NBR values are negative and much lower. From pixel 180 to 318, under the smoke, the AFRI values of the crops remain at the same high level as in the smoke-clear section, while both the NDVI and NBR produce low values.

Results
Throughout the study area, more than 1100 fires were counted by the JNF's foresters, rangers, and firefighters. Ground observations revealed that most of the wildfires could be classified as 'surface fire', consuming mostly the understory grasses, shrubs, litter, and woody material lying on the ground [47]. Occasionally, the surface fire reached and injured the trunks of mature trees. In eucalyptus plots, damages were also observed in the lower foliage ( Figure 4A). Pinus halepensis (Aleppo pine), the dominant pine species in the study area, is very sensitive to heat due to its relative thin bark and high amount of resin [48,49]. Therefore, many pines did not survive the surface fire ( Figure 4B). Only the minority of the fires were characterized as 'canopy fires'. In other species, such as acacia, mesquites, and tamarisk, the foliage was damaged, but the trees were recovered soon after the fire. For comparing the performance of the NBR, NDVI, and AFRI indices, a section of the S2 Level-2A image obtained on 10 July 2018, close to the Gaza Strip border, is presented in Figure 5. The transects represent a common situation when smoke, at different intensities, covers a variety of substrates-cultivated, bare soils, fire scars, and more. The true-color composite image (RGB = 0.665, 0.56, 0.49 µm) shows the open fire (light-orange hue), the biomass burning smoke (white hue), as well as burned scars that are a few days old (dark surfaces). The three indices were produced at 10-m spatial resolution, according to Equations (1), (2), and (5), from the surface reflectance values along a cross-section of 2771 m ( Figure 6, line A-A). This line was selected since it passes cultivated fields, bare soil, and was overcast by light smoke that characterized the entire region along the whole study period. This line is subdivided into several segments. From pixel 0 to 50 (and similarly from pixel 180 to 190) over the agricultural field where no smoke exists, the AFRI values accurately mimic those of the NDVI, but the NBR values are significantly lower. From pixel 50 to 180, over the bare soil, AFRI values are somehow higher than those of the NDVI. However, the NBR values are negative and much lower. From pixel 180 to 318, under the smoke, the AFRI values of the crops remain at the same high level as in the smoke-clear section, while both the NDVI and NBR produce low values. Figure 7A shows the AFRI and NBR profiles along line B-B ( Figure 5). This line was selected since it passes burned forest and cultivated fields. As already described in Figure 6, it can be noticed that the AFRI has a similar profile to NBR but with higher and positive values ( Figure 7A). Furthermore, the latter index glides the zero value for the entire transect, expect in the crop area (from pixel 300 to 330). Figure 7B,C depict the profiles of dAFRI and dNBR and RdAFRI and RdNBR computed by Equations      Figure 7A shows the AFRI and NBR profiles along line B-B ( Figure 5). This line was selected since it passes burned forest and cultivated fields. As already described in Figure 6, it can be noticed that the AFRI has a similar profile to NBR but with higher and positive values ( Figure 7A). Furthermore, the latter index glides the zero value for the entire transect, expect in the crop area (from pixel 300 to 330). Figure 7B,C depict the profiles of dAFRI and dNBR and RdAFRI and RdNBR computed by Equations (6), (3), (8), and (4), respectively, between 5 July 2018 and 25 July 2018. Figure 8 shows the false-color composite (RGB = 0.842, 0.665, 0.56 µm), the RdAFRI, and the RdNBR images for a portion of the study area obtained on 15 June 2018. The reference of the difference is the image acquired on 5 June 2018. It can be noticed that both RdAFRI and RdNBR emphasize the severity of the fire scars. The dark areas in the false-color composite are similar to the bright regions in the RdAFRI and RdNBR. However, in the RdNBR image, many pixels out of the burned area appear to have high index values that mislead any burned area interpretation. On the other hand, the RdAFRI image is smoother. Consequently, it is easier to select the threshold value (manually or from a systematics methodology) that flags a pixel as burned or not. By way of example, Figure 9 depicts the classification of the fire scars obtained by four levels of severity that are commonly used [29]. Table 2 details the values of NI for RdAFRI and RdNBR images as calculated from Equation (9) for selected dates and in the area with the lower burned severity index, i.e., RdAFRI <5.8. High values of NI mean high discrepancy between the index values with respect to their average values. Evidently, the noise degree in the RdNBR images is two orders of magnitude greater than the same calculated in the RdAFRI image.  systematics methodology) that flags a pixel as burned or not. By way of example, Figure 9 depicts the classification of the fire scars obtained by four levels of severity that are commonly used [29].   (9) for selected dates and in the area with the lower burned severity index, i.e., RdAFRI <5.8. High values of NI mean high discrepancy between the index values with respect to their average values. Evidently, the noise degree in the RdNBR images is two orders of magnitude greater than the same calculated in the RdAFRI image.  Remote Sens. 2020, 12, x 11 of 21 systematics methodology) that flags a pixel as burned or not. By way of example, Figure 9 depicts the classification of the fire scars obtained by four levels of severity that are commonly used [29].   (9) for selected dates and in the area with the lower burned severity index, i.e., RdAFRI <5.8. High values of NI mean high discrepancy between the index values with respect to their average values. Evidently, the noise degree in the RdNBR images is two orders of magnitude greater than the same calculated in the RdAFRI image.         Table 3 shows that the total burned area, within the JNF forests, is 7.02 km 2 , which corresponds to 13.79% of the total JNF area. It is worth mentioning that in some cases, fires were set in areas that have already been set on fire and were previously counted as scars. For example, the fire on 10 July 2018 covered an area of 0.38 km 2 , but only 0.15 km 2 were detected in new areas that were not burned before. The values of the thresholds THR1 and THR2 were chosen manually by visual inspection of the 24 S2 images.
For demonstration, Figure 11 presents the results of the proposed classification algorithm, i.e., AFRI, obtained for three successive S2 images, on 11 May, 21 May, and 5 June 2018. For each date, the false-color image (RGB = 0.842, 0.665, 0.56 µm) and the corresponding burned map image are provided. 23 12 December  Table 3 shows that the total burned area, within the JNF forests, is 7.02 km 2 , which corresponds to 13.79% of the total JNF area. It is worth mentioning that in some cases, fires were set in areas that have already been set on fire and were previously counted as scars. For example, the fire on 10 July 2018 covered an area of 0.38 km 2 , but only 0.15 km 2 were detected in new areas that were not burned before. The values of the thresholds THR1 and THR2 were chosen manually by visual inspection of the 24 S2 images.
For demonstration, Figure 11 presents the results of the proposed classification algorithm, i.e., AFRI, obtained for three successive S2 images, on 11 May, 21 May, and 5 June 2018. For each date, the false-color image (RGB = 0.842, 0.665, 0.56 µm) and the corresponding burned map image are provided. Blue pixels indicate clouds, whereas pink pixels indicate burned areas (only within the JNF study polygons, delineated by the yellow lines). It can be noticed that some burned areas can be identified comparing two successive false-color images (black area).
As described in Section 2.3, a field campaign was conducted on 23 July 2018 in order to monitor the forest in the proximity of the Gaza Strip. The result of this campaign is vector data that locates all the burned territories inside the JNF areas, as defined in Figure 1. These ground-truth vector data were used to realize an accuracy assessment of the results obtained on 25 July 2018 by the proposed classification algorithm. Around 5000 points have been randomly distributed within all JNF areas, where each area has several points proportional to its relative area. The confusion matrix [50] that summarizes the classification performance is shown in Table 4. The first column of the matrix indicates that 750 points have been classified by the algorithm as burned areas. Between them, 563 have been Figure 11. False-color (RGB = 0.842, 0.665, 0.56 µm) images with the corresponding burned maps obtained by the proposed AFRI algorithm from 11 May 2018 to 5 June 2018. Blue pixels indicate clouds, whereas pink pixels indicate burned areas (only within the JNF study polygons, delineated by the yellow lines). It can be noticed that some burned areas can be identified comparing two successive false-color images (black area).
As described in Section 2.3, a field campaign was conducted on 23 July 2018 in order to monitor the forest in the proximity of the Gaza Strip. The result of this campaign is vector data that locates all the burned territories inside the JNF areas, as defined in Figure 1. These ground-truth vector data were used to realize an accuracy assessment of the results obtained on 25 July 2018 by the proposed classification algorithm. Around 5000 points have been randomly distributed within all JNF areas, where each area has several points proportional to its relative area. The confusion matrix [50] that summarizes the classification performance is shown in Table 4. The first column of the matrix indicates that 750 points have been classified by the algorithm as burned areas. Between them, 563 have been detected as "Burned areas" also by the ground-truth vector data, whereas 187 have been classified as "Not burned". The positive predictive precision is 75%, and the negative predictive value is 93%, whereas the total accuracy is 90%. Table 5 shows the statistics of the affected vegetation species in the JNF forest areas on 22 December 2018. For example, 16.36% of the eucalyptus trees were burned.  Figure 12 presents the spatial distribution of the accumulated fire scars (red pixels) in the JNF forests as for 22 December 2018 obtained by the proposed classification algorithm. Figure 13A,B is zoomed in to show the areas in more detail.   In order to compare the classification results obtained by the methodology presented in Section 2.4 with respect to an equivalent algorithm based on RdNBR, the Otsu threshold method was applied. A smaller area (5.3 km 2 ) was considered with respect to the JNF forest presented in Figure 1. The area presents a uniform vegetation background (Figure 14). Among all 24 dates in the previous analysis, only 3 dates were selected: 5 June, 15 June, and 25 June 2018. The main reason is that in this area, the majority of the fire scars appears only on these 3 dates; therefore, the results could be compared with the ground measurement. Figure 14A-C show the ground measurements of the accumulated fire scars obtained by using the RdNBR-based algorithm and the RdAFRI-based algorithm, respectively. Since the RdNBR has high values even if ρ NIR approaches ρ SWIR2 (Figure 8), we have processed the classification results of the RdNBR-based algorithm with a majority filter tool in ArcGIS in order to replace pixel classification based on the majority of their contiguous neighboring pixels ( Figure 14D).
Finally, Tables 6-8 summarize the classification performances of the three methods. Around 1500 points have been randomly distributed within the study area of Figure 14. It can be noticed that the results obtained by RdAFRI are better than those obtained by RdNBR: the positive predictive precision is 86% and 69%, the negative predictive value is 87% and 84%, whereas the total accuracy is 87% and 80%, respectively. On the other hand, the RdNBR-based algorithm with a filter tool performs similarly to the proposed algorithm: the positive predictive precision is 78%, the negative predictive value is 86%, and the total accuracy is 84%.   Finally, Tables 6-8 summarize the classification performances of the three methods. Around 1500 points have been randomly distributed within the study area of Figure 14. It can be noticed that the results obtained by RdAFRI are better than those obtained by RdNBR: the positive predictive precision is 86% and 69%, the negative predictive value is 87% and 84%, whereas the total accuracy is 87% and 80%, respectively. On the other hand, the RdNBR-based algorithm with a filter tool performs similarly to the proposed algorithm: the positive predictive precision is 78%, the negative predictive value is 86%, and the total accuracy is 84%.

Discussion
Wildfire is a complex issue, since it is related to changes in ecosystems, climate, land-use and land-cover, and management practices, and it also has many socio-economic implications [51,52]. Therefore, maps of fire scars provide valuable information for studying forestry, agriculture, and pedology, as well as climate change. Since Earth observation data are considered to be the most informative means for quantifying the scars, a large variety of space systems have been used for this purpose. Selecting the most appropriate satellite depends on several criteria. Obviously, for the continental to global scales, the large-swath satellites (e.g., NOAA-AVHRR, MODIS) are more suitable. In contrast, for a regional scale, high spatial resolution and narrow swath are preferable. The second criterion is the temporal resolution required for the specific application. The number of fires, burned areas, as well as derived environmental consequences, such as the amount of CO 2 emissions, are usually summed for a monthly, seasonally, or yearly period using the large-swath satellites. However, for fast response and accurate mapping (e.g., for insurance assessment), more frequent images from high-resolution satellites (e.g., Sentinel-2, VENµS, RapidEye, etc.) are essential. Thirdly, the most meaningful spectral index for enhancing the burned signal and differentiating severity levels needs to be identified. Concerning the study objectives, S2 was found to be a suitable satellite due to its characteristics-high spatial resolution of 10 m, high temporal resolution of about five days, and the SWIR bands for computing the AFRI.
AFRI was initially developed as an approved vegetation index. While NDVI is based on the red band that is more sensitive to atmospheric aerosols, AFRI relies upon the SWIR band(s) that enable the radiation to penetrate a polluted atmospheric column. Empirical linear relationships between the red and the SWIR2 band reveal that ρ red ≈ 0.5ρ SWIR2 . Consequently, AFRI performs similar to NDVI, producing identical values under clear-sky conditions. However, under an opaque atmosphere influenced by biomass burning smoke, the NDVI values drop. Still, the AFRI continuously detects the real index values without the need for explicit correction for the aerosol effect.
In the current study, AFRI and its related algorithms were successfully used for delineating post-fire scars. While the AFRI values resembled the NDVI values when no smoke exists, the NBR values were significantly lower and even produced negative values over bare soils. Thoughtfully, the difference between positive and negative values enables the detection of the fire scars. In fact, beneath the smoke, the values of the NBR were as low as those of the NDVI. At the same time, the AFRI always kept high index values, disregarding even the light smoke that characterized the region during the nine months of the Kit Terror. The dAFRI and dNBR showed identical behavior both over unburned and burned areas. On the other hand, even if both successfully detected the burned area, the RdAFRI and RdNBR presented different behavior. Along the entire study area, the RdAFRI appeared to be smoother than the RdNBR, which was characterized by many noisy pixels. Mathematically, this phenomenon occurs when the NBR was computed on the pre-fire images-i.e., when the denominator of Equation (4) attains zero. As a result, some areas can be erroneously considered as fire scars, even if the dNBR value is low. Since the ρ SWIR2 is multiplied by a factor of 0.5, this behavior does not occur in the RdAFRI index profile. The noise levels of the images produced by the RdAFRI algorithm were two orders of magnitude lower than the images produced by RdNBR. Consequently, by applying the RdAFRI, it is possible to identify manual thresholds for distinguishing among the four severity categories.
During the eight months of this study, the fire scars appeared in small areas, i.e., around 2% maximum of the whole forest area. The images of the entire area presented a different type of background (e.g., diverse types of vegetation). Consequently, the "automatic" Otsu threshold did not succeed in accurately classifying these small burned patches correctly concerning the whole study area. On the other hand, Otsu's methodology worked very well if applied to images that circumscribe the burned areas, where the background was more uniform. In order to have an overview of the fire scars in the entire JNF forest, we opted for manual thresholds. These values were not unique for all the dates. The main reason is that the RdAFRI, but also the RdNBR, was sensible to the clouds over the pre and post-burned area. Even if the cloud mask product was used in order to avoid applying the algorithm at cloudy pixels, it appeared that many cloudy areas were not correctly flagged by the mask and were analyzed by the algorithm. Therefore, in the corresponding high cloud coverage images, a higher threshold was required in order to locate the burned area pixels concerning the cloud pixels correctly. Therefore, an accurate cloud mask product is a key factor for the accuracy of burned area mapping algorithms. For companion, the proposed algorithm was also implemented with the cloud screening and atmospheric correction procedure provided by the MACCS-ATCOR Joint Algorithm (MAJA) processor [43] that obtained similar results. In this regard, we admit that the main drawback of the proposed procedure (and also of the RdNBR) is determining the threshold level. This disadvantage can be overcome by more accurate filed calibration [35,53].

Conclusions
The study was derived by the need to map a sequence of more than 1100 wildfires that last nine months, partially over forest areas. Obviously, Earth observation is the preferable means for monitoring fire scars. Sentinel-2 was selected due to its considerable high spatial and temporal resolution, as well as suitable spectral bands. Naturally, two obstacles prevent accurate scar mapping. The first is the biomass burning smoke that continuously overcast the study area. This issue was solved by using the SWIR bands that can penetrate fine-mode atmospheric polluted aerosols. The AFRI spectral index, which uses the 2.1 µm band, is an atmospheric-resistant index that can observe vegetation through the smoke. Furthermore, its values closely resemble the NDVI values under clear-sky conditions, and it retains the same level of index values under smoke. In the latter regard, RdAFRI has proven to assess burned areas accurately as the equivalent RdNBR index with the advantage of not requiring sieved or filtering tools to remove erroneously fire scars pixel. Ground validation of the RdAFRI values reveals an overall accuracy of 90%. The RdAFRI enables one to distinguish among four severity categories; however, due to different cloud cover in two consecutive dates, an automatic determination of a threshold level was not possible. Therefore, two threshold levels were considered through visual inspection and manually assigned to each imaging date. It is concluded that the S2 AFRI-derived is the desirable means for calculating and mapping the spatio-temporal dynamics of the fire scars along with the statistics of the burned vegetation species. RdAFRI was found applicable, particularly when the area is affected by many fires that broke out at different places across the region and contaminated the air away from the fire center.