Applying the Tropical Peatland Combustion Algorithm to Landsat-8 Operational Land Imager (OLI) and Sentinel-2 Multi Spectral Instrument (MSI) Imagery

This study establishes a new technique for peatland fire detection in tropical environments using Landsat-8 and Sentinel-2. The Tropical Peatland Combustion Algorithm (ToPeCAl) without longwave thermal infrared (TIR) (henceforth known as ToPeCAl-2) was tested on Landsat-8 Operational Land Imager (OLI) data and then applied to Sentinel-2 Multi Spectral Instrument (MSI) data. The research is aimed at establishing peatland fire information at higher spatial resolution and more frequent observation than from Landsat-8 data over Indonesia’s peatlands. ToPeCAl-2 applied to Sentinel-2 was assessed by comparing fires detected from the original ToPeCAl applied to Landsat-8 OLI/Thermal Infrared Sensor (TIRS) verified through comparison with ground truth data. An adjustment of ToPeCAl-2 was applied to minimise false positive errors by implementing pre-process masking for water and permanent bright objects and filtering ToPeCAl-2’s resultant detected fires by implementing contextual testing and cloud masking. Both ToPeCAl-2 with contextual test and ToPeCAl with cloud mask applied to Sentinel-2 provided high detection of unambiguous fire pixels (>95%) at 20 m spatial resolution. Smouldering pixels were less likely to be detected by ToPeCAl-2. The detected smouldering pixels from ToPeCAl-2 applied to Sentinel-2 with contextual testing and with cloud masking were only 35% and 56% correct, respectively; this needs further investigation and validation. These results demonstrate that even in the absence of TIR data, an adjusted ToPeCAl algorithm (ToPeCAl-2) can be applied to detect peatland fires at 20 m resolution with high accuracy especially for flaming. Overall, the implementation of ToPeCAl applied to cost-free and available Landsat-8 and Sentinel-2 data enables regular peatland fire monitoring in tropical environments at higher spatial resolution than other satellite-derived fire products.


Introduction
Tropical peatland is a significant source of global carbon [1,2]. The state of tropical peatland is mainly disturbed by conversion of peatland forest to other uses, especially for plantations. Peatland drainage [3] and peatland fires [4] have contributed to accelerating atmospheric carbon increases consequential climate change [5,6]. Smouldering peatland fires burn organic soil and thick fuels (e.g., decaying leaves, tree branches and trunks), which release more carbon and smoke with greater particulate concentration than flaming fires [7,8]. Drought, unmanaged drainage, and land use change from natural peatland habitat to agricultural/plantation areas have led to high flammability conditions of degraded peatlands [9][10][11][12][13][14][15][16]. Peatland fires together with agriculture, forestry, and other land use are the main contributing sources (about 51%) to greenhouse gases emissions in Indonesia [17]. Figure 1 shows the estimated annual emission of peatland fires during 1997-2019 in the Equatorial Asian region including Indonesia and neighbouring countries [18]. The peak emissions from peatland fires, which occurred in 1997,2002,2006,2009,2015, and 2019, corresponded to the warm episodes of El Niño Southern Oscillation (ENSO) and the positive modes of Indian Ocean Dipole, as shown in the Oceanic Niño Index (ONI) and Dipole Mode Index (DMI) recorded by Climate Prediction Center and Physical Sciences Laboratory of National Oceanic and Atmospheric Administration, respectively [19,20]. In general, positive modes of these climate indices were associated with extreme droughts and contributed to Indonesian fires [21,22]. Beside the influences of extreme drought, mainly, peatland fires in Indonesia have anthropogenic origins [23][24][25]. Therefore, it is important to develop a comprehensive monitoring of peatland fires for control and minimisation of carbon emission. www.globalfiredata.org/data.html) and the Oceanic Niño Index (ONI) and Dipole Mode Index (DMI) using 3 month running average of sea surface temperature anomalies during November-December-January period. Climate indices data were accessed at https://origin.cpc.ncep.noaa.gov/products/analysis_monitoring/ensostuff/ONI_v5.php and https://psl.noaa.gov/gcos_wgsp/Timeseries/Data/dmi.had.long.data.
Remote sensing satellites have been widely used for biomass burning detection. Research on remote sensing satellite applications, particularly for peatland fire detection and its assessment in Indonesia has been conducted from medium spatial resolution to low spatial resolution satellite sensors and from Unpiloted Aerial Vehicle (UAV) as shown in Table 1. Remote sensing satellites have been widely used for biomass burning detection. Research on remote sensing satellite applications, particularly for peatland fire detection and its assessment in Indonesia has been conducted from medium spatial resolution to low spatial resolution satellite sensors and from Unpiloted Aerial Vehicle (UAV) as shown in Table 1. Table 1. List of peatland fire detection research using medium to low spatial resolution satellite data and unpiloted aircraft in Indonesia (OLI is Operational Land Imager; TIRS is Thermal Infrared Sensor; BIRD is Bi-spectral InfraRed Detection; HSRS is HotSpot Recognition System; TET-1 is Technology Experiment Carrier; Suomi NPP is Suomi National Polar-orbiting Partnership; SWIR is shortwave infrared; TIR is longwave thermal infrared; MWIR is mid-wave infrared; VIIRS is Visible Infrared Imaging Radiometer Suite; UAV is Unpiloted Aerial Vehicle).  16 30 Peatland fire detection at nighttime & daytime [26,27] BIRD HSRS MWIR,  Peatland fire detection at daytime [28] TET-1 TET-1 MWIR 2-5 160 Peatland fire detection at daytime [29] Suomi NPP VIIRS SWIR 0.5 750

Satellite/Aircraft
Peatland fire detection at nighttime [30] MWIR, TIR 0.5 375 Assessment of active fire [31] UAV Tau2 TIR TIR as needed 0.185 Peatland fire detection at daytime [32] Cost-free satellite data with moderately frequent observations have potential to be used in regular peatland fires information collection. The cost-free and available medium spatial resolution (20-30 m) e.g., Sentinel-2 Multi Spectral Instrument (MSI) and Landsat-8 Operational Land Imager (OLI)/ Thermal Infrared Sensor (TIRS) regularly observe the Earth in the morning at about 10 am every five and 16 days, respectively. At lower spatial resolution (375 m) VIIRS is also available free of cost and provides the highest spatial resolution global active fire product for daytime and night-time among low spatial resolution cost-free active fire products. However, because fire products are not available for Sentinel-2, there exists an opportunity for the development of a detection algorithm for this sensor. A fire detection algorithm, particularly for peatland fires in tropical environments, e.g., Indonesia, called Tropical Peatland Combustion Algorithm (ToPeCAl), has been established using aerosol, shortwave infrared (SWIR) and longwave thermal infrared (TIR) spectra and applied to Landsat-8 OLI/TIRS daytime [27].
The possibility of merging fire information from medium spatial resolution data with low spatial and high temporal resolution data, Visible Infrared Imaging Radiometer Suite (VIIRS) 375 m, was undertaken via a comparison with ToPeCAl fire information from Landsat-8 OLI/TIRS and resulted in accuracies up to 90% for fires in Indonesia's peatland [31]. Meanwhile, many of the combined fire products from middle spatial resolution sensor satellites are specifically designed for burnt area mapping [33][34][35][36]. A global active fire algorithm from middle spatial resolution data has been developed by [37][38][39][40] but is still limited for specific tropical peatland fires environments [26,27]. Moreover, research on spectral comparison between middle spatial resolution data, e.g., Landsat-8 and Sentinel-2, in many regions, e.g., conterminous United States, Australia, Southern Africa, Europe, and China showed good relationships [41][42][43][44][45]. However, this relationship needs to be explored over equatorial tropical environments.
Higher spatial and temporal information of peatland fires is potentially possible by applying ToPeCAl to Sentinel-2 MSI data. However, Sentinel-2 MSI only has SWIR for fire detection, while ToPeCAl originally used SWIR and TIR in Landsat-8 OLI/TIRS. An approach of using SWIR from medium spatial resolution data for fire detection has been explored by [37][38][39][40]. It may be possible to apply ToPeCAl without TIR to Sentinel-2 MSI for peatland fires detection in tropical environments.

Landsat-8 and Sentinel-2A/B
The National Aeronautics and Space Administration (NASA) and the Department of Interior United States Geological Survey (USGS) developed and launched the Landsat 8 satellite on 11 February 2013 with two science instruments, OLI and TIRS. Nine spectral bands, which span from the visible (VIS), the near infrared (NIR), the SWIR and the cirrus band at 30 m spatial resolution and 15 m for the panchromatic band, are collected in OLI. TIRS collects two bands of TIR spectra at 100 m spatial resolution. Landsat-8 with its near-polar and Sun-synchronous orbit provides a repeat earth observation every 16-days at an equatorial crossing time of 10.00 ± 15 min in a 185 km swath [46,47]. In this study, Landsat-8 data were accessed from https://earthexplorer.usgs.gov/.
Sentinel-2 satellites were developed by the European Commission and the European Space Agency (ESA). Sentinel-2A and 2B were launched on 23 June 2015 and 7 March 2017, respectively. The on-board MSI has 13 spectral bands ranging from the VIS, the NIR, and the SWIR at different spatial resolutions (10, 20 and 60 m). The orbit of Sentinel-2 satellites is Sun-synchronous with observation times close to Landsat local overpass time at 10.30 (descending orbit). The swath width of Sentinel-2 MSI is 290 km with a frequent revisit of five days at the equator from Sentinel-2A and Sentinel-2B. The Sentinel-2 MSI product is provided in tiles which cover 10,000 km 2 . In a single day of acquisition, the Sentinel-2 MSI product may consist of one or more tiles in a certain area [48]. Sentinel-2 data were accessed from https://scihub.copernicus.eu/dhus/#/home. The bands used for applying ToPeCAl to Landsat-8 OLI/TIR and Sentinel-2 MSI and background mapping, e.g., water, cloud, and bright object are listed in Table 2.  Nine pairs of Landsat-8 and Sentinel-2 acquired on the same day were used in this study (Table  4) to explore the second objective. The pairs of Landsat-8 and Sentinel-2 were used to undertake a spectral reflectance comparison prior to applying and assessing the application of ToPeCAl-2 to  Nine pairs of Landsat-8 and Sentinel-2 acquired on the same day were used in this study (Table 4) to explore the second objective. The pairs of Landsat-8 and Sentinel-2 were used to undertake a spectral reflectance comparison prior to applying and assessing the application of ToPeCAl-2 to Sentinel-2 MSI; this constituted step 2 in the analysis. Data pairs with fire features analysed from the RGB image Table 4. List of Landsat-8 and Sentinel-2 data for applying and assessing application of ToPeCAl-2 to images over Central Kalimantan, South Kalimantan, West Kalimantan, Riau, Jambi, and South Sumatera. There are four pairs of Landsat-8 and Sentinel-2 images contained fires in the same areas on the day of acquisition. Descending orbit Landsat-8 images were acquired at 2-12 min after Sentinel-2 image acquisition. The state of fire features was identified from the Red-Green-Blue (RGB) image composites of SWIR-2, SWIR-1, and Red bands of Sentinel-2 and Landsat-8 images.

Tropical Peatland Combustion
Algorithm with Shortwave Infrared (SWIR) and Thermal Infrared (TIR) Figure 3 illustrates the principle of Planck's Law which describes spectral radiance simulation using the MODerate resolution atmospheric TRANSmission (MODTRAN-5) for peatland with limited vegetation under assumption of ground albedo of 0.1 [52] and emissivity of 0.9 [53] at given temperatures of 300 K, 360K, 450 K, 500 K, 600 K, 700 K, and 1000 K under a tropical atmosphere and using a rural model with visibility of 3 km and hazy conditions. Previous studies observed that peatlands combustion temperature occurred at 360-500 K for smouldering, 500-700 K for smouldering but close to flaming, and at 700-1000 K for flaming [26,53,54]. The atmospheric windows where remote sensing of the biosphere best operates (light blue shading in Figure 3) shows the radiant energy of flaming peatlands (>700 K) clearly rising in SWIR region (1.6-2.2 µm) and reaching its peak in MWIR region (3-4 µm). Smouldering peatland spectral radiance falls in 360-500 K area, showing less energy in SWIR than MWIR. The presence of cool fires can also be observed in the Longwave Infrared (LWIR) which shows higher radiance than in SWIR. However, smouldering fires close to flaming (500-700 K) have less energy at LWIR compared with MWIR. The background temperature (~300 K) has a more significant emission in LWIR region. Thus, the wavelengths with high potential to observe peatland combustion using satellite remote sensing are in the 1.6-2.2 µm, 3-4 µm and 8-12 µm zones. However, reflected solar energy in SWIR may interfere with the detection of flaming at daytime, while the background temperature may affect the detection of cool fire in LWIR [55]. with an assumption of emissivity of 0.9 over burning peatland at 300 K, 360 K, 450 K, 500 K, 700 K and 1,000 K under a tropical atmosphere and using a rural model. The atmospheric transmittance is presented in light blue background, while water vapour and gas absorption is portrayed in grey. This graph assumes 100% of the pixels corresponding to objects are at the nominated temperature and that might be unrealistic for some flaming at 30 × 30 m pixels, but it may be realistic for smouldering pixels. The graph is not definitive-it gives suggestions as to where to observe peatland fires using remote sensing.
To examine the radiance profile of peatland fire from Landsat-8 OLI/TIRS, [27] analysed the top of atmosphere (TOA) radiance for peatland fires locations based on ground truth. The ground truth data were collected on the same day and up to two days after Landsat-8 observation over Central Kalimantan peatlands on 28 September 2018 at path/row of 118/062. The ground truth data were collected by helicopter from Disaster Management Agency, drone imaging for inaccessible locations, and by direct measurement for accessible peatland fires [27]. Figure 4 presents examples of areas of interest in burnt zones and smouldering peatland for extracting TOA radiances from Landsat-8 based on drone imagery in Central Kalimantan. Drone imagery in Figure 4a observed burnt areas on the same day with Landsat-8 observation on 28 September 2018. The burnt area is shown in bright brown tones from the drone image and also from surface pictures which were taken by using digital camera (Figure 4b). There was no smoke emanating from burnt areas and no smoky atmosphere in the area at the time these images were taken. Meanwhile, Figure 4c shows smouldering peatland which is shown in dark areas mixed with white tones, indicating ash as the product of smouldering [7,8,56]. The drone imagery of smouldering peatlands was taken two days after Landsat-8 observed the area on 28 September 2018. The surface picture of a part of smouldering peatlands (Figure 4d) was taken at the same time as drone imagery ( Figure 4c). Smoky atmosphere existed in the areas when the field measurements were taken and some smoke was still emanating from the smouldering peatlands ( Figure 4d). The accessible smouldering peatland areas were measured for surface temperature resulting in an average of 452 K by using a handheld IRtek IR150 dual-beam laser infrared thermometer (Figure 4e) while the average surface temperature of burnt area was about 315 K. with an assumption of emissivity of 0.9 over burning peatland at 300 K, 360 K, 450 K, 500 K, 700 K and 1000 K under a tropical atmosphere and using a rural model. The atmospheric transmittance is presented in light blue background, while water vapour and gas absorption is portrayed in grey. This graph assumes 100% of the pixels corresponding to objects are at the nominated temperature and that might be unrealistic for some flaming at 30 × 30 m pixels, but it may be realistic for smouldering pixels. The graph is not definitive-it gives suggestions as to where to observe peatland fires using remote sensing.
To examine the radiance profile of peatland fire from Landsat-8 OLI/TIRS, [27] analysed the top of atmosphere (TOA) radiance for peatland fires locations based on ground truth. The ground truth data were collected on the same day and up to two days after Landsat-8 observation over Central Kalimantan peatlands on 28 September 2018 at path/row of 118/062. The ground truth data were collected by helicopter from Disaster Management Agency, drone imaging for inaccessible locations, and by direct measurement for accessible peatland fires [27]. Figure 4 presents examples of areas of interest in burnt zones and smouldering peatland for extracting TOA radiances from Landsat-8 based on drone imagery in Central Kalimantan. Drone imagery in Figure 4a observed burnt areas on the same day with Landsat-8 observation on 28 September 2018. The burnt area is shown in bright brown tones from the drone image and also from surface pictures which were taken by using digital camera (Figure 4b). There was no smoke emanating from burnt areas and no smoky atmosphere in the area at the time these images were taken. Meanwhile, Figure 4c shows smouldering peatland which is shown in dark areas mixed with white tones, indicating ash as the product of smouldering [7,8,56]. The drone imagery of smouldering peatlands was taken two days after Landsat-8 observed the area on 28 September 2018. The surface picture of a part of smouldering peatlands (Figure 4d) was taken at the same time as drone imagery (Figure 4c). Smoky atmosphere existed in the areas when the field measurements were taken and some smoke was still emanating from the smouldering peatlands ( Figure 4d). The accessible smouldering peatland areas were measured for surface temperature resulting in an average of 452 K by using a handheld IRtek  The graph in Figure 5 shows the profile of TOA radiances of Landsat-8 that were extracted from multiple pixels at a total of 1.62 ha of flaming, 2.43 ha of smouldering and 1.62 ha of burnt areas. Each The graph in Figure 5 shows the profile of TOA radiances of Landsat-8 that were extracted from multiple pixels at a total of 1.62 ha of flaming, 2.43 ha of smouldering and 1.62 ha of burnt areas. Each of flaming, smouldering, and burnt area covered at least >80-90% of the pixel areas. The area of interest of smouldering and burnt areas were imaged by drone, while for flaming based features on large fire areas from 7-6-4 RGB image composites of Landsat-8 on 28 September 2018 in Central Kalimantan. The profile for flaming shows high spectral radiance in all bands of Landsat-8 (band 1-7, band-10) compared with smouldering and burnt area radiance. Smouldering showed a higher radiance at NIR (band-5) than burnt areas due to higher scattering from smoky surfaces and also some un-burnt vegetation still existing in the ongoing smouldering peatland which reflected more NIR than in the burnt area. In SWIR bands, burnt area reflectance was higher than smouldering area in band-6, but shown a similar value at band-7. The burnt area with less vegetation than in smouldering peatland may have lower surface moisture which leads to a higher reflection of SWIR in burnt areas. In TIR, smouldering peatland shows a slightly higher radiance than burnt area. Observation of the spectral plots suggests that for separating flaming, smouldering, and burnt area, taking a ratio of band-7 and band-6 may be productive. The spectral signatures in Figure 5 were acquired for guiding peatland fire validation from Landsat-8 images for some of locations which were not accessible. These profiles supported the development of ToPeCAl [27].
In addition, previous tropical peatland fire studies using field experiments and using remote sensing data were reviewed when selecting samples of peatland fires. In respect to acquisition times the low speed of peatland fire propagation permits sampling of smouldering areas on different days; in our case the time difference was up to 2 days after Landsat-8 acquisition. Low surface peatland fire propagation speed was reported at 3.8 cm/h (0.9 m/day) on average, or 6.5 cm/h (1.6 m/day) on maximum through field experiments in Central Kalimantan [54]. Meanwhile, the authors of [29] showed active fire propagation using TET-1 data, lasting for about 30 days, at 4 m/day (minimum) to 1,903 m/day (maximum) in Central Kalimantan. From personal communication with local staff of the Indonesian Disaster Management Agency, who accompanied ground truth work on 28 September-1 October 2018, peatland fires can last more than one week if no attempt to extinguish is made. Smouldering is a very persistent type of combustion, which is relatively easy to commence and difficult to suppress [7,8]. Once ignition has occurred, smouldering can exist for long periods (days, months or years). When peatland is ignited, the smouldering combustion dominates the flaming combustion [7,8] and its propagation develops both laterally and vertically [54,56]. Smouldering and flaming combustion can lead to one other; however, when the oxygen supply increases, as occurs with an increase in wind velocity, combustion becomes faster, temperature increases and a flaming state may eventuate [7,8].
ToPeCAl was first developed using SWIR and TIR spectral bands for peatland fire detection with a classification of clear and smoky/hazy tropical atmospheric conditions from the aerosol band (0.433-0.453 µm) of Landsat-8 over Indonesia [27]. The TOA reflectance of aerosol band, the TOA reflectance of SWIR bands (1.566-1.651 µm and 2.107-2.294 µm), and the TOA brightness temperature of TIR band (10.61-11.19 µm) were used in ToPeCAl. TOA reflectance for fire detection was introduced in previous studies using Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER) and Landsat-8 by [37][38][39][40]. ToPeCAl uses TIR and SWIR to enhance the detection of smouldering and flaming combustion of tropical peatland fires [27]. ToPeCAl applied to Landsat-8 OLI/TIRS provided 82% accuracy compared with ground truth data with omission errors in small smouldering areas of non-homogeneous pixels less than one pixel of Landsat-8 and commission errors corresponding to bright building rooftops. Elimination of the urban area (building rooftops) from the analysis using a landuse map generated from Landsat-8 increased the accuracy to 93% for ToPeCAl fire mapping [27]. However, false detection may also occur for smouldering detection resulting from ToPeCAl, as the brightness temperature of smouldering peatlands showed slight differences with burnt areas ( Figure 5).  In ToPeCAl [27], the potential combustion areas are first identified by use of the Shortwave Infrared Combustion Index (SICI ρ ), the ratio of TOA reflectance of SWIR-2 at the central wavelength of 2.2 µm (ρ 2.2 ) to SWIR-1 at the central wavelength of 1.6 µm (ρ 1.6 ). SICI ρ values greater than 1 are classified as potential combustion areas, while the non-combustion areas, e.g., vegetation, water, urban area, and cloud, were masked using SICI ρ < 1. Then, flaming and smouldering pixels are identified within the potential combustion areas using TOA reflectance of SWIR-2 (ρ 2.2 ) and TOA brightness temperature of TIR (BT 10.9 ) on relatively clear or smoky/hazy atmospheric conditions based on TOA reflectance of the aerosol band (ρ 0.44 ). The ToPeCAl formulation is given by equations 1−12 [27]: Peatland combustion types in relatively-clear atmospheric conditions: Peatland combustion types in smoky/hazy atmosphere condition: AND ρ 0.44 ≥ 0.27 (8) AND 0.11 ≥ ρ 2.2 ≤ 0.32 AND BT 10.9 ≥ 297 (smouldering) (9) 0.32 < ρ 2.2 < 0.47 AND BT 10.9 ≥ 297 (mix of flaming and smouldering) (10) AND SICI ρ ≤ 1 AND ρ 2.2 ≥ 0.47 AND BT 10.9 ≥ 303 (flaming for close to saturated pixel of SWIR radiance) (12) In this study, the ToPeCAl algorithm with SWIR and TIR is termed ToPeCAl-1, while ToPeCAl without TIR band is termed ToPeCAl-2. An initial comparison of ToPeCAl-1 and ToPeCAl-2 was conducted through visual analysis on the ToPeCAl fires product compared with 6-5-4 (SWIR-1, NIR, Red) RGB images of Landsat-8. In this visual analysis, the attention was to check what fire objects resulting from ToPeCAl-2 were not detected by ToPeCAl-1; that is, the false positives of ToPeCAl-2. From the 6-5-4 RGB images of Landsat-8, most of the false positive of ToPeCAl-2 fires occurred in the edge of cloud, permanent bright objects such as urban features or factory rooftops and mining areas ( Table 5). This is indicative of false positive fire detection using ToPeCAl-2 for sensors with no TIR, e.g., Landsat-8 OLI and Sentinel-2 MSI. Table 5. Visual comparison on false positive fire pixels of ToPeCAl with SWIR and TIR (ToPeCAl-1) and ToPeCAl without TIR (ToPeCAl-2) applied to Landsat-8 images.

ToPeCAl-1 ToPeCAl-2 Recommendation for Eliminating the False Positive Pixels on ToPeCAl-2
Edge of cloud no yes Contextual test or cloud masking Permanent bright objects less more permanent bright object masking

Tropical Peatland Combustion Algorithm without TIR (ToPeCAl-2)
ToPeCAl-1 formulation was adjusted to perform ToPeCAl-2. The general flowchart is shown in Figure 6. The adjustment of ToPeCAl-2 included: • Pre-processing for water body and permanent bright objects masking, • A variation of ToPeCAl-2 with contextual test for comparing candidate fire pixels with background pixels [37,38], • An alternative variation to ToPeCAl-2 using cloud masking.

Tropical Peatland Combustion Algorithm without TIR (ToPeCAl-2)
ToPeCAl-1 formulation was adjusted to perform ToPeCAl-2. The general flowchart is shown in Figure 6. The adjustment of ToPeCAl-2 included: • Pre-processing for water body and permanent bright objects masking, • A variation of ToPeCAl-2 with contextual test for comparing candidate fire pixels with background pixels [37,38], • An alternative variation to ToPeCAl-2 using cloud masking. Figure 6. Flowchart of general methods used for ToPeCAl-2 compared with ToPeCAl-1. Two methods were assessed, i.e., ToPeCAl-2 with contextual test and ToPeCAl-2 with cloud masking. F is flaming; FS is mix of flaming and smouldering; S is smouldering. TOA is top of atmosphere; ρ is reflectance; BT is brightness temperature; SWIR is shortwave infrared; TIR is thermal infrared.

Water Masking
Water and permanent bright objects are outliers or potential outliers in the analysis for fire detection; thus, they were masked, firstly before applying ToPeCAl-2. Normalised Difference Water Index (NDWI) and Modified Normalised Difference Water Index (MNDWI) are widely used for water body mapping from middle spatial resolution data, e.g., Landsat-8 and Sentinel-2 [43,50,57]. NDWI uses green ( . ) and NIR ( . ) reflectance while MNDWI uses green and SWIR ( . ) Figure 6. Flowchart of general methods used for ToPeCAl-2 compared with ToPeCAl-1. Two methods were assessed, i.e., ToPeCAl-2 with contextual test and ToPeCAl-2 with cloud masking. F is flaming; FS is mix of flaming and smouldering; S is smouldering. TOA is top of atmosphere; ρ is reflectance; BT is brightness temperature; SWIR is shortwave infrared; TIR is thermal infrared.

Water Masking
Water and permanent bright objects are outliers or potential outliers in the analysis for fire detection; thus, they were masked, firstly before applying ToPeCAl-2. Normalised Difference Water Index (NDWI) and Modified Normalised Difference Water Index (MNDWI) are widely used for water body mapping from middle spatial resolution data, e.g., Landsat-8 and Sentinel-2 [43,50,57]. NDWI uses green (ρ 0.56 ) and NIR (ρ 0.865 ) reflectance while MNDWI uses green and SWIR (ρ 1.61 ) reflectance in their formulation (Equation (13)). We found that NDWI can map water bodies well in clear atmospheric conditions, but for some images with many fires and smoky conditions, NIR does not perform well due to smoke aerosol scattering; thus, MNDWI with SWIR wavelengths which are able to be better detected through smoke was used to map the water bodies. Both NDWI and MNDWI were applied in this study to overcome the variety of smoke conditions over multiple fire images. According to [50], NDWI values within 0 to 0.1 and MNDWI values within 0.2 to 0.35 can map water bodies well. In this study, NDWI > 0.1 and MNDWI > 0.35 was used to map water bodies. Some cloud shadows which have the same spectral character with water also were mapped using these thresholds (Equation (13)).

Permanent Bright Object Masking
Previously, in the application of ToPeCAl-1 an urban GIS map was used for masking potential false positive pixels in urban areas [27]. However, multitemporal high spatial resolution images from Google Earth confirmed that the false positive pixels were still found on building or factory rooftops outside the urban area and these objects are relatively permanent, e.g., building/factory and mining areas which existed through several years.
These permanent bright reflective objects can be more effectively eliminated to minimise the false positive pixels through the use of TOA reflectance of SWIR-2 (ρ 2.2 ). The mean value of ρ 2.2 > 0.2 was applied by [38] on multitemporal data of Landsat-8 OLI in co-located cloud free pixels in the previous 6 months (176 days) to eliminate false positives over permanent bright objects. However, it is rare to obtain cloud free pixels from multitemporal data of Landsat-8 for every 16-days in tropical environments. Thus, some adjustments were applied in this study to map the permanent bright objects with the following procedure: • Use of multitemporal Sentinel-2 MSI images instead of Landsat-8 OLI images to acquire more images containing cloud free pixels. This is a result of the higher temporal resolution of Sentinel-2 acquisition compared with Landsat-8. The cloud free pixels of the multitemporal images from Level-1C Sentinel-2 MSI for permanent bright object mapping were obtained via the Google Earth Engine platform [58]. According to [48], geometric image quality performance for multitemporal Sentinel-2 20 m bands was 6 m (2σ). To match with the fire product of ToPeCAl-2 applied to Landsat-8, the product of permanent bright objects from SWIR-2 band of Sentinel-2 was resampled to 30 m.

•
The median values of ρ 2.2 as applied by [40] instead of the mean values as used by [38]; thus, in this study, the median values of ρ 2.2 > 0.18 were used. This condition needs to exist in the first half-year (January-June) and in the second half-year (July-December) either in the year of data acquisition, or in the previous year. This seasonal period was selected to gain increased possibilities of obtaining cloud free images which were not based on the period of seasons in Indonesia. Otherwise, cloud free images were usually found in dry season (April-September) and barely any cloud free data in wet season (October-March). Permanent bright objects should also have relatively stable and high SWIR reflectance in any period. Thus, these conditions can be expressed as: AND ((ρ 2.2 at tt (jul−dec, y) OR ρ 2.2 at t (jul−dec, y−1) ) > 0.18 (15) then TRUE as bright permanent objects. Where t (jan-june, y) and t (jan-june, y−1) are time period in January-June in on-going year and in previous year; t (jul-dec, y) and t (jul-dec, y−1) are time period in July-December in on-going year and in previous year.

ToPeCAl-2 with Contextual Test
The form of ToPeCAl-2 for relatively clear sky and smoky/hazy atmospheres is similar to ToPeCAl-1 but without BT 10.9 . Unambiguous detection is for flaming fires (Equation (16)). For detection of flaming close to saturated pixel, SICI ρ values in ToPeCAl-1 (Equations (6) and (12)) were expressed as ≤1 with conditions for ρ 2.2 and BT 10.9 . This SICI ρ ≤ 1 resulted from ρ 1.6 or ρ 2.2 ≥ 1 and with a condition of ρ 1.6 ≥ ρ 2.2 . The actual values of SICI ρ for close to saturated pixel were not explicitly stated in ToPeCAl-1 (Equations (6) and (12)), but because TIR does not exist in ToPeCAl-2 then it should be rewritten as ≥0.9 (Equation (17)). Meanwhile additional candidate fire pixels are from a mix of flaming and smouldering (FS) and smouldering pixels (Equations (18) and (19)).
Unambiguous fire pixels on clear sky or smoky/hazy atmosphere: (16) Fire pixels at close to saturation for flaming pixel on clear sky or smoky/hazy atmosphere: Candidate fire pixels for FS pixels on clear sky or smoky/hazy atmosphere: Candidate fire pixels for smouldering pixels on clear sky or smoky/hazy atmosphere: A contextual algorithm was implemented to compare the candidate fire pixels from Equations (18) and (19) with background pixels. However, the unambiguous fire pixel (flaming) from Equation (16) and the close to saturated flaming pixels from Equation (17) were retained as peatland fire pixels. The contextual algorithm was sourced from [38] with a modification particularly on using the mean and standard deviation of the ratio of ρ 2.2 and ρ 1.6 rather than the ratio of ρ 2.2 and ρ 0.865 as used in [38]. The peatland fire pixels from ToPeCAl-2 were derived from the ratio of ρ 2.2 to ρ 1.6 and the reflectance of SWIR-2 (ρ 2.2 ), thus the statistical calculations of the background pixels were also analysed based on these spectra.
The contextual search was set to a 61 × 61 pixel window similar to [37,38]. Valid background pixels ideally should not be contaminated with outlier pixels. However, in tropical peatland environments the background is mostly heterogeneous, potentially resulting in a non-normal distribution statistic of the background values. Dense cloud might also exist and also the fire pixels including candidate fire pixels in 61 × 61 pixel windows. These features can also influence the background statistics. In this study, the cloud mask product from Landsat-8 was not chosen as it mapped slightly larger cloud areas which misclassified some clear area as clouds. This result also was reported by [51] who suggested to use the Cloud Detection Algorithm Generation (CDAG) to perform better cloud masking. Thus, CDAG with single band of Red reflectance (ρ 0.65 ) was applied to map clouds in this study. The formulation of cloud mask is given in Equation (20) [51].
Previous studies [37,38,40] only masked water/cloud shadow and unambiguous fire pixels from background analyses. In this study, permanent bright objects, dense cloud pixels and all of fire pixels, including candidate pixels from the population of the background were also masked. An example is shown in Figure 7c, but with no bright permanent object pixels existing in this sample. The histogram and population statistics of the ratio of reflectance of band-7 to band-6 and the reflectance of band-7 resulting from the image data corresponding to perform better cloud masking. Thus, CDAG with single band of Red reflectance (ρ . ) was applied to map clouds in this study. The formulation of cloud mask is given in Equation (20) [51].
Previous studies [37,38,40] only masked water/cloud shadow and unambiguous fire pixels from background analyses. In this study, permanent bright objects, dense cloud pixels and all of fire pixels, including candidate pixels from the population of the background were also masked. An example is shown in Figure 7c, but with no bright permanent object pixels existing in this sample. The histogram and population statistics of the ratio of reflectance of band-7 to band-6 and the reflectance of band-7 resulting from the image data corresponding to Figure   The same RGB image as (a) overlaid with water/cloud shadow mask (blue) and obvious fire pixels (yellow) as per-previous studies [37,38,40]; saturated flaming pixels and permanent bright objects did not exist in this sample; (c) The same with (b) but with additional masks for cloud (light green) and candidate fire pixels (pink); (d) The histogram and population statistics of the ratio reflectance band-7 ( . ) to band-6 ( . ) and reflectance of band-7 of background pixels which excluded water/cloud shadow and unambiguous fire pixels as shown spatially in (b); (e) The histogram and population statistics of the same variable with (d) but with additional masks of cloud The same RGB image as (a) overlaid with water/cloud shadow mask (blue) and obvious fire pixels (yellow) as per-previous studies [37,38,40]; saturated flaming pixels and permanent bright objects did not exist in this sample; (c) The same with (b) but with additional masks for cloud (light green) and candidate fire pixels (pink); (d) The histogram and population statistics of the ratio reflectance band-7 (ρ 2.2 ) to band-6 (ρ 1.6 ) and reflectance of band-7 of background pixels which excluded water/cloud shadow and unambiguous fire pixels as shown spatially in (b); (e) The histogram and population statistics of the same variable with (d) but with additional masks of cloud and candidate fire pixels. The ratio of band-7/band-6 and reflectance of band-7 data (d,e) do not follow a normal distribution according the Anderson-Darling goodness of fit statistic at confidence level of 99.7% (3σ), corresponding to a heterogeneous background, but data in (e) is more normal then (d).
After selecting valid background pixels, the candidate fire pixels from Equations (18) and (19) were tested by using a contextual algorithm (Equation (21)). The fire pixels that passed this test were merged with fire pixels resulting from Equations (16) and (17) as the peatland fires pixels.
Contextual test for candidate fire pixels: are the mean and standard deviation of the ratio ρ 2.2 to ρ 1.6 respectively, while ρ 2.2 and σρ 2.2 are the mean and standard deviation of ρ 2.2 respectively from 61 × 61 pixel window centred on each potential fire pixel. The values of 0.8 and 0.08 in Equation (21) were sourced from [38] to avoid false alarms in unusually homogeneous regions [37].

ToPeCAl-2 with Cloud Masking
ToPeCAl-2 with cloud masking was another approach to eliminate false positives in cloud, dense cloud and the edge of cloud, instead of the contextual test approach. In this second approach, ToPeCAl-2 fires resulted from Equation (16) to Equation (20) were filtered with cloud masking, particularly for FS and smouldering pixels, to remove false positives. However, the cloud mask data was not produced from Equation (21); rather it was obtained from cloud mask product of Level-2 pixel-qa of Landsat-8. The cloud mask from pixel-qa provided a larger cloud area which assisted in removing the false positive pixels in the edge of cloud.
3.2.5. Pixel-Based Comparison of ToPeCAl-2 Applied to Landsat-8 Operational Land Imager (OLI) with ToPeCAl-1 Applied to Landsat-8 OLI/Thermal Infrared Sensor (TIRS) A pixel-based comparison of 15 Landsat-8 images was conducted between the fire products of ToPeCAl-2 with the fire products of ToPeCAl-1. The reference fire data was ToPeCAl-1 fires. The comparison of peatland fires was calculated for three classes of peatland combustions and one combined class of fire. The three classes of peatland combustions from ToPeCAl are flaming (F), smouldering (S), and mix of flaming and smouldering (FS). The one fire class (Fire) was created by integrating the three peatland combustion classes into one class. Statistical analyses generated the probability of detection (POD), false positive (commission error) and false negative (omission error). POD was calculated by the ratio of true positive fires between ToPeCAl-2 and ToPeCAl-1 to the total fires produced in the reference data (ToPeCAl-1). Commission error was counted from the ratio of the ToPeCAl-2 fires which were not detected in ToPeCAl-1 to the total ToPeCAl-2 fires, while omission error can be computed as 1-POD.
In addition to ToPeCAl-1 Landsat-8 fires as the reference, others image data, e.g., Google Earth and PlanetScope were used to confirm any indication of fires or false detections. The PlanetScope data were accessed from Planet explorer at www.planet.com.

Sentinel-2 MSI Data Pre-Processing
The TOA reflectance of Sentinel-2 MSI data were obtained from Level-1C. Some of the Sentinel-2 listed in Table 4 consisted of more than one tile for the corresponding Landsat-8 scene, thus the data of the Sentinel-2 tiles were mosaicked with nearest neighbor resampling method to preserve the original values of data. In the overlap zones, the images with less cloud cover were selected as the primary data. The tiles were mosaicked for each band of aerosol (ρ 0.443 ), green (ρ 0.560 ), NIR (ρ 0.865 ), SWIR-1(ρ 1.610 ) and SWIR-2 (ρ 2.190 ). The spatial resolution of aerosol and green bands was resampled to 20 m to match with spatial resolution of the SWIR bands to produce fire pixels. Pre-processing to mask out water and permanent bright objects was conducted by using Equation (13) to Equation (15). A general flowchart, which is presented in Figure 6, was also applied to Sentinel-2 to produce fire pixels with an initial investigation of sensor comparison for both Landsat-8 OLI and Sentinel-2 MSI and is described in Section 3.3.2.

Sentinel-2A MSI and Landsat-8 OLI Spectral Reflectance Sensor Comparison
In order to apply ToPeCAl-2 developed from Landsat-8 OLI data, an initial comparison of the difference between spectral reflectance from Landsat-8 OLI and Sentinel-MSI at band-1, SWIR-1 and SWIR-2 was conducted. Figure 8 shows that the spectral response function of Sentinel-2A/2B MSI slightly varied compared with Landsat-8 OLI at aerosol band, SWIR-1 and SWIR-2 bands. The spectral response data was sourced from ESA [59,60]. original values of data. In the overlap zones, the images with less cloud cover were selected as the primary data. The tiles were mosaicked for each band of aerosol (ρ . ), green (ρ . ), NIR (ρ . ), SWIR-1(ρ . ) and SWIR-2 (ρ . ). The spatial resolution of aerosol and green bands was resampled to 20 m to match with spatial resolution of the SWIR bands to produce fire pixels. Pre-processing to mask out water and permanent bright objects was conducted by using Equation (13) to Equation (15). A general flowchart, which is presented in Figure 6, was also applied to Sentinel-2 to produce fire pixels with an initial investigation of sensor comparison for both Landsat-8 OLI and Sentinel-2 MSI and is described in Section 3.3.2.

Sentinel-2A MSI and Landsat-8 OLI Spectral Reflectance Sensor Comparison
In order to apply ToPeCAl-2 developed from Landsat-8 OLI data, an initial comparison of the difference between spectral reflectance from Landsat-8 OLI and Sentinel-MSI at band-1, SWIR-1 and SWIR-2 was conducted. Figure 8 shows that the spectral response function of Sentinel-2A/2B MSI slightly varied compared with Landsat-8 OLI at aerosol band, SWIR-1 and SWIR-2 bands. The spectral response data was sourced from ESA [59,60]. Pairs of Sentinel-2A/2B MSI pixels and Landsat-8 OLI pixels on the same day of acquisition also exhibiting fire features in the same area in both data were compared at aerosol, SWIR-1 and SWIR-2 bands. All bands used in ToPeCAl were set at the same spatial resolution of 20 m using the nearest neighbor resampling method in projection of UTM/WGS84. In this study, a valid pixel selection for comparing TOA reflectance of Sentinel-2 and Landsat-8 was built from the same data acquisition utilising criteria which were provided in Scene Classification (SC) of Sentinel-2 from Level-2A through Sen2Cor and pixel-qa of Landsat-8 from Level-2. The valid pixels were from vegetation and non-vegetation class based on the SC data.
In addition to the cloud masking, the physical environment close to cloud generates cloud adjacency effects or 3-D radiative effects. The cloud adjacency effects cause a brighter pixels around the edge of clouds and an increased high-biased aerosol optical depth thus contributing to an increase in observed reflectance [61][62][63][64]. The parallax error due to a difference of earth viewing angles observed from Sentinel-2 spectral bands also was mentioned as a possible displacement of cloud mask [61]. In this study, an additional 40 m (2 pixels) buffer from the cloud mask was generated and applied to masking the cloud adjacency effects in the remaining imagery. The cloud mask was merged from medium and high confidence cloud from Sentinel-2 data.
Pairs of Sentinel-2A/2B MSI pixels and Landsat-8 OLI pixels on the same day of acquisition also exhibiting fire features in the same area in both data were compared at aerosol, SWIR-1 and SWIR-2 bands. All bands used in ToPeCAl were set at the same spatial resolution of 20 m using the nearest neighbor resampling method in projection of UTM/WGS84. In this study, a valid pixel selection for comparing TOA reflectance of Sentinel-2 and Landsat-8 was built from the same data acquisition utilising criteria which were provided in Scene Classification (SC) of Sentinel-2 from Level-2A through Sen2Cor and pixel-qa of Landsat-8 from Level-2. The valid pixels were from vegetation and non-vegetation class based on the SC data.
In addition to the cloud masking, the physical environment close to cloud generates cloud adjacency effects or 3-D radiative effects. The cloud adjacency effects cause a brighter pixels around the edge of clouds and an increased high-biased aerosol optical depth thus contributing to an increase in observed reflectance [61][62][63][64]. The parallax error due to a difference of earth viewing angles observed from Sentinel-2 spectral bands also was mentioned as a possible displacement of cloud mask [61].
In this study, an additional 40 m (2 pixels) buffer from the cloud mask was generated and applied to masking the cloud adjacency effects in the remaining imagery. The cloud mask was merged from medium and high confidence cloud from Sentinel-2 data.
Meanwhile, the mask data from pixel-qa of Landsat-8 were cloud, cloud shadow, water, and high confidence cirrus. All of the mask data selected from SC data and pixel-qa data were merged as one mask and then applied to Sentinel-2A/2B MSI and Landsat-8 OLI images.
The valid pixels of TOA reflectance data of Sentinel-2A/2B and Landsat-8 OLI in aerosol, SWIR-1 and SWIR-2 bands were then statistically compared and analysed using t-tests to determine if there was a significant difference between the means of Sentinel-2 and Landsat-8 TOA reflectance. Ordinary least squares (OLS) regression analysis was also conducted to explore the relationships using the slope and intercept coefficient values between the TOA reflectance of the two sensors. The goodness of fit and significance of the OLS regressions were defined by the coefficient of determination and statistic p-value.

Implementation of ToPeCAl-2 Applied to Sentinel-2 MSI
ToPeCAl-2 was applied to original values of TOA reflectance of Sentinel-2 MSI and to the transformed values of TOA reflectance of Sentinel-2 MSI resulting from the OLS regression equations between Landsat-8 and Sentinel-2 at band aerosol, SWIR-1 and SWIR-2. As tested in Landsat-8 OLI, ToPeCAl-2 applied to Sentinel-2 was also adjusted by implementing two approaches, i.e., ToPeCAl-2 with contextual test and ToPeCAl-2 with cloud masking.
The formulation of ToPeCAl-2 with contextual test was based on to Equation (16) to Equation (21). However, the cloud mask was not generated from Equation (20). The cloud product from SC map of Sentinel-2 was used for background analysis in the contextual test and for removing false positive pixels from FS and smouldering pixels in the second approach of ToPeCAl-2 adjustment, ToPeCAl-2 with cloud masking.
The masking data used for calculating ToPeCAl-2 was different from the masks data for analysing the spectral relationship between Sentinel-2 and Landsat-8 in the Section 3.3.2. The mask data used in ToPeCAl-2 application was high confidence cloud from SC map with additional buffer zone of five pixels (100 m) from outside high confidence cloud polygon. According to [65], the high and medium confidence cloud has above 65% and 35% of cloud probability, respectively. The buffer created around dense cloud was also aimed to mask the cloud adjacency effects close to dense cloud. The medium confidence cloud, which was not adjacent to high confidence cloud, was not masked out from the calculation of ToPeCAl-2 in order to detect possible fire pixels under the thin cloud.

Pixel-Based Comparison of ToPeCAl-2 Applied to Sentinel-2 MSI with ToPeCAl-1 Applied to Landsat-8 OLI/TIRS
The comparison of fire products of ToPeCAl-2 applied to Sentinel-2 MSI were conducted on a pixel-basis with the fire products of ToPeCAl-1 applied to Landsat-8 OLI/TIRS as reference. The statistical comparisons conducted were POD, false positive (commission error), and false negative (omission error). Potential geometric error existing between both Sentinel-2 and Landsat-8, which is about 38 m (2σ) [66] and peatland fires propagation zone during the difference in time of observation between Landsat-8 and Sentinel-2 were considered for analysing the zones of true positive fires. According to [29], the peatland fires propagation can be approximated at 1.3 m/min. This value was then multiplied by the time difference between acquisition of Landsat-8 and Sentinel-2 to estimate the potential spread of peatland fires between the observations. Thus, fire pixels located in the zone of geometric error plus the zone of peat fire propagation from the true positive fires were counted as spatially related true positive fire pixels. Others fire features from higher spatial resolution images, e.g., Google Earth and PlanetScope were also used to confirm any indication of fires or false detections.
Remote Sens. 2020, 12, x FOR PEER REVIEW 20 of 37 IOE resulted from ToPeCAl-2 with cloud masking (2%). Detail statistical values can be found in Appendix- Table A1. Figure 9b shows the comparison of ToPeCAl-1 and ToPeCAl-2 from two approaches in three classes of peatland fires. Both of F and FS pixels resulting from ToPeCAl-2 with either contextual test or cloud masking provided high POD (93-99%). Meanwhile, the accuracy of detection of smouldering pixels from ToPeCAl-2 with contextual test was less than from ToPeCAl-2 with cloud masking. The POD of smouldering pixels from ToPeCAl-2 with contextual test was only 35% with IOE of 65%. Referring to the result shown in Figure 9a, the low POD of ToPeCAl-2 with contextual test was related to the poorer detection of smouldering pixels, as shown in Figure 9b. Detailed statistical values of this comparison can be found in Appendix-Table A2.  Figure 10 shows an example comparing the results of ToPeCAl-2 from the two approaches compared to ToPeCAl-1 in Riau's peatland. There are obvious fires features visually shown in orange to yellow tones in the false composite of RGB Landsat-8 image (Figure 10a) and Figure 10b shows the fire areas resultant from ToPeCAl-1 as the fire reference. ToPeCAl-2 with contextual test (Figure 10c) resulted in less true positive (TP) pixels than from ToPeCAl-2 with cloud masking (Figure 10d). The TP from ToPeCAl-2 with contextual test which are shown in red colour were general spatially correlated with flaming pixels from ToPeCAl-1 (Figure 10b), while the TP from ToPeCAl-2 with cloud masking (Figure 10d) corresponds to all the fires in ToPeCAl-1. Figure 10e is the natural RGB composite of PlanetScope data with spatial resolution of 3.7 m about 11 min before Landsat-8 observation on 29 August 2019. Fire features cannot be displayed from the natural RGB image of PlanetScope, but smoke as an indication of fires, shown from this RGB image, confirmed peatland fires detected by ToPeCAl-1 and ToPeCAl-2. An enlarged area (red square) of burning peatlands on PlanetScope image was made and overlaid with fires resulting from ToPeCAl-2 with contextual test (Figure 10f) and ToPeCAl-2 with cloud masking (Figure 10g). Most of FN pixels from ToPeCAl-2 with contextual test (Figure 10f) were in the area of burning peatlands with smoke features in Figure 10h. Meanwhile, Figure 10i shows PlanetScope image at the same area with Figure 10h after three days on 1 September 2019. Some of the burnt areas in dark tones were covered by smoke which could be sourced from smouldering peatlands or surrounding areas. However, based on the experience observed during ground truth and drone mapping in Central Kalimantan, the burnt peatlands were still emitting smoke after 2-3 days [27].  Figure 9b shows the comparison of ToPeCAl-1 and ToPeCAl-2 from two approaches in three classes of peatland fires. Both of F and FS pixels resulting from ToPeCAl-2 with either contextual test or cloud masking provided high POD (93-99%). Meanwhile, the accuracy of detection of smouldering pixels from ToPeCAl-2 with contextual test was less than from ToPeCAl-2 with cloud masking. The POD of smouldering pixels from ToPeCAl-2 with contextual test was only 35% with IOE of 65%. Referring to the result shown in Figure 9a, the low POD of ToPeCAl-2 with contextual test was related to the poorer detection of smouldering pixels, as shown in Figure 9b. Detailed statistical values of this comparison can be found in Appendix A- Table A2. Figure 10 shows an example comparing the results of ToPeCAl-2 from the two approaches compared to ToPeCAl-1 in Riau's peatland. There are obvious fires features visually shown in orange to yellow tones in the false composite of RGB Landsat-8 image (Figure 10a) and Figure 10b shows the fire areas resultant from ToPeCAl-1 as the fire reference. ToPeCAl-2 with contextual test (Figure 10c) resulted in less true positive (TP) pixels than from ToPeCAl-2 with cloud masking (Figure 10d). The TP from ToPeCAl-2 with contextual test which are shown in red colour were general spatially correlated with flaming pixels from ToPeCAl-1 (Figure 10b), while the TP from ToPeCAl-2 with cloud masking (Figure 10d) corresponds to all the fires in ToPeCAl-1. Figure 10e is the natural RGB composite of PlanetScope data with spatial resolution of 3.7 m about 11 min before Landsat-8 observation on 29 August 2019. Fire features cannot be displayed from the natural RGB image of PlanetScope, but smoke as an indication of fires, shown from this RGB image, confirmed peatland fires detected by ToPeCAl-1 and ToPeCAl-2. An enlarged area (red square) of burning peatlands on PlanetScope image was made and overlaid with fires resulting from ToPeCAl-2 with contextual test (Figure 10f) and ToPeCAl-2 with cloud masking (Figure 10g). Most of FN pixels from ToPeCAl-2 with contextual test (Figure 10f) were in the area of burning peatlands with smoke features in Figure 10h. Meanwhile, Figure 10i shows PlanetScope image at the same area with Figure 10h after three days on 1 September 2019. Some of the burnt areas in dark tones were covered by smoke which could be sourced from smouldering peatlands or surrounding areas. However, based on the experience observed during ground truth and drone mapping in Central Kalimantan, the burnt peatlands were still emitting smoke after 2-3 days [27].

Relationship Between Top of Atmosphere (TOA) Reflectance of Landsat-8 and Sentinel-2 for Aerosol and SWIR Bands
The number of pixel pairs for TOA reflectance values comparison between Landsat-8 and Sentinel-2 ranged between two million and 27.7 million. The aerosol band showed high variability in regression slope values with the lowest value at 0.14 in Jambi and South Sumatera sites and the highest value at 1.07 in a Central Kalimantan site. The data pairs with weak relationships for the aerosol bands also experienced lower R 2 values in the SWIR-1 and SWIR-2 bands compared with the higher R 2 data for other regressions. Each data pair with higher R 2 values (greater than 0.70 at SWIR-2) were merged into one data set with 58.7 million pairs of pixels and analysed with OLS regression. The merged data provided 0.66, 0.78, and 0.83 for R 2 at aerosol, SWIR-1, and SWIR-2 bands, respectively. The results of OLS these regression analyses were significant at p-values < 0.001. The statistical values can be found in Appendix A- Table A3. Table 6 summaries the t-test statistic resulting from the mean differences between TOA reflectance of Landsat-8 OLI and TOA reflectance of Sentinel-2 MSI with transformation and no transformation at aerosol, SWIR-1 and SWIR-2 bands. In general, the Sentinel-2 TOA reflectance values with transformation and no transformation were slightly greater than Landsat-8 OLI. The mean difference of TOA reflectance from the two sensors with no transformation at aerosol, SWIR-1, and SWIR-2 bands were 0.0045, 0.0171, and 0.0086, respectively with p-values < 0.001 indicating statistical significance. Meanwhile, the mean differences of the transformed TOA reflectance values of Sentinel-2 MSI sourced from the merged OLS model and the local OLS model with the TOA reflectance of OLI were lower than the mean difference resulting from Sentinel-2 TOA reflectance with no transformation. Table 6. The mean differences of TOA reflectance of Sentinel-2 Multi Spectral Imager (MSI) and Landsat-8 OLI at aerosol, SWIR-1, and SWIR-2 bands. Total data (n) were compiled from three data sets representing study sites in Central Kalimantan and South Kalimantan (CK-SK), Jambi and South Sumatera (JSS), and Riau. * Resulted from t-test with significant different at p-values < 0.001.

OLS Model
Site Acquisition Date n (Pixel)

Implementation and Comparison of ToPeCAl-2 Applied to Sentinel-2 MSI
As mentioned in Table 4 only four of the nine pairs of data sets were found with fire features on both Sentinel-2 and Landsat-8 images. ToPeCAl-2 with contextual test and ToPeCAl-2 with cloud masking were applied to more than 128 million pixels of Sentinel-2 MSI at 20 m spatial resolution. Figure 11a shows the performance of ToPeCAl-2 fires resultant from Sentinel-2 MSI compared with the ToPeCAl-1 fires of Landsat-8 OLI/TIRS. Based on the three OLS transformation models, the highest POD resulted from no transformation value for both ToPeCAl-2 with contextual test and ToPeCAl-2 with cloud mask. Overall, the POD of ToPeCAl-2 with cloud mask (88%) was higher than ToPeCAl-2 with contextual test (62%). However, ToPeCAl-2 with cloud mask showed a higher independent commission error (ICE) and a lower independent omission error (IOE) than ToPeCAl-2 with contextual test. The ICE of ToPeCAl-2 with cloud masking and contextual test from the no transformation model were 21% and 10%, while the IOE were 12% and 38%, respectively. Detailed statistical values can be found in Appendix A- Table A4.
Remote Sens. 2020, 12, x FOR PEER REVIEW 23 of 37 with contextual test. The ICE of ToPeCAl-2 with cloud masking and contextual test from the no transformation model were 21% and 10%, while the IOE were 12% and 38%, respectively. Detailed statistical values can be found in Appendix- Table A4.
To distinguish the origin of classes of fire in ICE and IOE pixels, Figure 11b presents the ToPeCAl fires comparison shown in Figure 11a in the three peatland fires classes, i.e., flaming, mix of flaming and smouldering (FS) and smouldering. Flaming pixels were well detected (≥95%) in ToPeCAl-2 from either ToPeCAl-2 with contextual test or ToPeCAl-2 with cloud masking. About 76-78% of FS pixels from ToPeCAl-2 with contextual test agreed with ToPeCAl-1, while 86-97% of FS pixels were found from ToPeCAl-2 with cloud masking. Smouldering pixels were less well-detected (32-56%) from both ToPeCAl-2 with highest detection from ToPeCAl-2 with cloud masking (56%). The ICE shows the highest error percentage in flaming or FS pixel, while IOE was high in smouldering pixel. Detail statistical values can be found in Appendix -Table A5.
In general, Figure 11 shows a lower smouldering pixel detection and a higher FS or flaming detection from ToPeCAl-2. ToPeCAl-2 with cloud masking showed a higher peatland fire detection for overall fires, but with higher false positive pixels in FS or flaming fire and lower false negative pixels in smouldering fire. Whilst ToPeCAl-2 with contextual test was likely to have lower false positive pixels in FS fire but, higher false negative in smouldering fire.  Figure 12 shows an example of the comparison of ToPeCAl-2 fires resulted from Sentinel-2 image with ToPeCAl-1 fires resulted from Landsat-8 image in Central Kalimantan. The Sentinel-2 image was acquired nine minutes before the Landsat-8 image. The RGB composite of SWIR-2, SWIR-1 and Red bands was applied to Landsat-8 image (Figure 12a) and to the Sentinel-2 image ( Figure  12b) to enhance the fire features under heavy smoke conditions; fires cannot be seen in the natural RGB composite as shown in Figure 12c. Figure 12d,e show fires resultant from ToPeCAl-2 with contextual test and ToPeCAl-2 with could masking applied to Sentinel-2 image, respectively. True positive (TP) fire pixels are shown in red colour in most of the fire areas. A boundary of spatially related-false positive fires (RFP) is shown in black bordered squares and spatially related-false negative fires (RFN) in cyan bordered squares shown adjacent to the true positive fires. Some of these RFP/RFN were within the distance of the geometric error from both Sentinel-2 and Landsat-8 and within the boundary of peatland fires propagation which was based on estimates from the previous research [29] at 1.3 m/m horizontally. These RFP and RFN pixels were counted as the true positive pixels. In Figure 12d, some ToPeCAl-2 fires with contextual test remained as false negative (FN) To distinguish the origin of classes of fire in ICE and IOE pixels, Figure 11b presents the ToPeCAl fires comparison shown in Figure 11a in the three peatland fires classes, i.e., flaming, mix of flaming and smouldering (FS) and smouldering. Flaming pixels were well detected (≥95%) in ToPeCAl-2 from either ToPeCAl-2 with contextual test or ToPeCAl-2 with cloud masking. About 76-78% of FS pixels from ToPeCAl-2 with contextual test agreed with ToPeCAl-1, while 86-97% of FS pixels were found from ToPeCAl-2 with cloud masking. Smouldering pixels were less well-detected (32-56%) from both ToPeCAl-2 with highest detection from ToPeCAl-2 with cloud masking (56%). The ICE shows the highest error percentage in flaming or FS pixel, while IOE was high in smouldering pixel. Detail statistical values can be found in Appendix A- Table A5.
In general, Figure 11 shows a lower smouldering pixel detection and a higher FS or flaming detection from ToPeCAl-2. ToPeCAl-2 with cloud masking showed a higher peatland fire detection for overall fires, but with higher false positive pixels in FS or flaming fire and lower false negative pixels in smouldering fire. Whilst ToPeCAl-2 with contextual test was likely to have lower false positive pixels in FS fire but, higher false negative in smouldering fire. Figure 12 shows an example of the comparison of ToPeCAl-2 fires resulted from Sentinel-2 image with ToPeCAl-1 fires resulted from Landsat-8 image in Central Kalimantan. The Sentinel-2 image was acquired nine minutes before the Landsat-8 image. The RGB composite of SWIR-2, SWIR-1 and Red bands was applied to Landsat-8 image (Figure 12a) and to the Sentinel-2 image (Figure 12b) to enhance the fire features under heavy smoke conditions; fires cannot be seen in the natural RGB composite as shown in Figure 12c. Figure 12d,e show fires resultant from ToPeCAl-2 with contextual test and ToPeCAl-2 with could masking applied to Sentinel-2 image, respectively. True positive (TP) fire pixels are shown in red colour in most of the fire areas. A boundary of spatially related-false positive fires (RFP) is shown in black bordered squares and spatially related-false negative fires (RFN) in cyan bordered squares shown adjacent to the true positive fires. Some of these RFP/RFN were within the distance of the geometric error from both Sentinel-2 and Landsat-8 and within the boundary of peatland fires propagation which was based on estimates from the previous research [29] at 1.3 m/m horizontally. These RFP and RFN pixels were counted as the true positive pixels. In Figure 12d, some ToPeCAl-2 fires with contextual test remained as false negative (FN) pixels in blue colour while from ToPeCAl-2 with cloud masking some of these FN pixels were classified as true positive pixels (red colour) in Figure 12e. An enlargement of fires area which classified as FN from ToPeCAl-2 with contextual test and TP from ToPeCAl-2 with cloud masking is presented in the PlanetScope image with natural colour composite (Figure 12f) on the same day with Sentinel-2 and Landsat-8. PlanetScope images have a higher spatial resolution (3.7 m) than Sentinel-2 but SWIR is not available in PlanetScope. However, it can show smoke as an indication of fires, as illustrated in Figure 12f.
Remote Sens. 2020, 12, x FOR PEER REVIEW 24 of 37 pixels in blue colour while from ToPeCAl-2 with cloud masking some of these FN pixels were classified as true positive pixels (red colour) in Figure 12e. An enlargement of fires area which classified as FN from ToPeCAl-2 with contextual test and TP from ToPeCAl-2 with cloud masking is presented in the PlanetScope image with natural colour composite (Figure 12f) on the same day with Sentinel-2 and Landsat-8. PlanetScope images have a higher spatial resolution (3.7 m) than Sentinel-2 but SWIR is not available in PlanetScope. However, it can show smoke as an indication of fires, as illustrated in Figure 12f. Another example in smaller peatland fires (about 104 ha) is presented in Figure 13. Peatland fires started with bluish smoke emanating from ground in the false colour composite of Sentinel-2 image (Figure 13a) and pixels detected as fire from ToPeCAl-2 with contextual test and ToPeCAl-2 with cloud masking in Figure 13b,c, respectively. The FN areas (blue colour) in Figure 13b,c correspond to the dark area with smoke shown by PlanetScope image on Figure 13d. The PlanetScope image was acquired three minutes after Sentinel-2 acquisition. Some smoke emanates from dark burning Another example in smaller peatland fires (about 104 ha) is presented in Figure 13. Peatland fires started with bluish smoke emanating from ground in the false colour composite of Sentinel-2 image (Figure 13a) and pixels detected as fire from ToPeCAl-2 with contextual test and ToPeCAl-2 with cloud masking in Figure 13b,c, respectively. The FN areas (blue colour) in Figure 13b,c correspond to the dark area with smoke shown by PlanetScope image on Figure 13d. The PlanetScope image was acquired three minutes after Sentinel-2 acquisition. Some smoke emanates from dark burning peatland and is one of the characteristics of smouldering peatlands which were observed in the field. The RFP and RFN fire pixels showed adjacent to TP areas resultant from ToPeCAl-2 corresponded to smoke features from the ground in PlanetScope image (Figure 13d). The RFN/RFP pixels adjacent to TP fires were classified as FS and smouldering pixels from ToPeCAl-1 of Landsat-8. The FS and smouldering pixels surrounding unambiguous fire pixels were found by [54] in their field experiment in Central Kalimantan's peatland. Peatland fires (smouldering or FS fires) were found at 20-30 m distance surrounding the front line of the burning area due to lateral propagation of smouldering fires at <20 cm depth as described in Figure 13e. This scheme for peatland fires propagation supports the results of ToPeCAl fires detection.
Remote Sens. 2020, 12, x FOR PEER REVIEW 25 of 37 distance surrounding the front line of the burning area due to lateral propagation of smouldering fires at <20 cm depth as described in Figure 13e. This scheme for peatland fires propagation supports the results of ToPeCAl fires detection.  [54] from field experiment in Central Kalimantan's peatland after ignition and during surface burning in the first stage; the propagation were also in depth spread in the next stage.

Comparison and Adjustment of ToPeCAl-2 to Landsat-8 OLI
The comparison and adjustment of ToPeCAl without TIR (ToPeCAl-2) was conducted in order to explore the potential of this algorithm for the detection of tropical peatland fires using Sentinel-2 MSI data. An initial test was accomplished by implementing ToPeCAl-2 compared with ToPeCAl-1 In this data pairs of Sentinel-2 and Landsat-8, the RFP and RFN pixels were within geometric error (38 m) and propagation zones (12 m) which about 50 m from TP pixels (magenta square line); (c) As per Figure (b) but from ToPeCAl-2 with cloud mask; (d) PlanetScope image in natural colour composite showed smoke emanating from fires areas from in the same area and in three minutes after Sentinel-2 s acquisition; (e) Surface peatland fires propagation in lateral (horizontal) direction according to [54] from field experiment in Central Kalimantan's peatland after ignition and during surface burning in the first stage; the propagation were also in depth spread in the next stage.

Comparison and Adjustment of ToPeCAl-2 to Landsat-8 OLI
The comparison and adjustment of ToPeCAl without TIR (ToPeCAl-2) was conducted in order to explore the potential of this algorithm for the detection of tropical peatland fires using Sentinel-2 MSI data. An initial test was accomplished by implementing ToPeCAl-2 compared with ToPeCAl-1 on the same source of Landsat-8 data. In this study, ground truth to validate the ToPeCAl-2 fires performance were not conducted; however, ToPeCAl-1 was validated by ground truth data and provided 82% accuracy which can be improved by removing permanent bright objects [27]. Although ToPeCAl-1 has SWIR and TIR to support peatland fire mapping, the spatial resolution of TIR in Landsat-8 is quite low to observe small peatland fires, often in mixed pixels which contain heterogeneous objects. However, smouldering pixels resultant from ToPeCAl-1, which are located adjacent to or surrounding unambiguous fire, may be considered as the true smouldering peatland. Smouldering is a very persistent type of combustion, which is relatively easy to commence and difficult to suppress [7] and its residual can last for days after flaming fires pass over [56]. Smouldering and flaming peatland fires can lead to one other, however, when the oxygen supply increases, as occurs with an increase in wind velocity, combustion becomes faster, temperature increases and a flaming state may eventuate [7,8]. Other than that, from field experiment, [54] found that smouldering peatlands were found in surrounding or adjacent un-slashed areas to the burning peatland. The information gained from laboratory and field experiments [7,8,49,52] for peatland fire characteristics support the interpretation of peatland fires resultant from satellite.
In addition, for confirming of ToPeCAl-2 fire detections, some verification occurred by observing the state of smoke emanating from fire areas which are shown in natural RGB images from PlanetScope on the same day with Landsat-8/Sentinel-2 as shown in Figures 10, 12 and 13. However, some of the smouldering or FS pixels detected by ToPeCAL-2 were difficult to be verified using specific smoke features in a generally smoky atmosphere, which can influence the quality of radiometry of PlanetScope RGB images. Further investigation on the radiometry of smoke features can potentially support the validation of ToPeCAl-2 for mapping peatland fires.
In this study, the adjustments of ToPeCAl were performed by adding a pre-processing for water masking, permanent bright object masking and performing some filtering methods on ToPeCAl-2 results, i.e., contextual test and cloud masking. ToPeCAl-2 with contextual test examining FS and smouldering fire pixels was compared to their background statistics using mean and standard deviation in 61 × 61 pixel window. In this contextual test, high confidence fire pixels were adopted from previous literature [37,38,40] by using statistical values of mean+3σ or 99.7% of confidence level for active fire detection. The fire pixels that passed this contextual test were highly likely to be active fire or flaming fire. Peatland fires particularly smouldering fire has lower combustion power and low temperature which is significantly different with flaming fires. A lower confidence level, e.g., 95% (mean+2σ) or 68% (mean+σ) for the background statistic can be tested to cover low confidence fires such as smouldering peatland fires. However, ToPeCAl-2 depends only on SWIR for detecting peatland fires, thus it may be possible to generate more false positive smouldering pixels by lowering the confidence level, for example, in burnt areas. Additional analysis such as smoke detection and analysis for smouldering peatland can be added to confirm smouldering fire pixels. In this study, ToPeCAl-2 with contextual test applied to Landsat-8 OLI provided low detection (55%) of fires compared with ToPeCAl-1 fires. From the results comparison using peatland fires types derived from ToPeCAl-1, the missed detection of ToPeCAl-2 with contextual test was for smouldering fire at about 45%. Flaming and FS fire pixels were well detected by ToPeCAl-2 with contextual test (≥93%).
A second approach was conducted by implementing cloud mask to ToPeCAl-2 fires. Initial assessment of ToPeCAl-2 fires compared with RGB image of Landsat-8 using visual analysis showed that false positive pixels occurred in cloud and at the edge of cloud. These false positive pixels were from FS and smouldering fire pixels detected by ToPeCAl-2. Thus, the cloud mask was only applied to FS and smouldering pixels, while unambiguous fire pixels (flaming and close to saturated flaming) were classified as true fire pixels. Dense cloud with highly reflective aerosols tends to brighten the cloud [67] which might be classified as false detection in ToPeCAl-2. The uncertainty in observation on the edge of cloud is commonly termed "cloud adjacency effects" or "3-D radiative effects". These can increase the reflectance around the edge of cloud with an increase of aerosol properties [62][63][64]. The cloud adjacency effect has also exhibited as a problem in other optical remote sensing applications, e.g., ocean monitoring [68] and atmospheric monitoring [62,64]. The results of ToPeCAl-2 with cloud mask provided higher percentage of peatland fire detection (98%) because only FS and smouldering pixels occurring in cloud areas were removed as fire pixels. Lower counts of false positive and false negative pixel resulted from this model. Most of FS and smouldering pixels that passed cloud masking were located in areas surrounding true positive fire pixels. However, the smouldering pixels resultant from ToPeCAl-2 may include false positive pixels in burnt areas as they have similar SWIR reflectance as shown in Figure 5.
Although most permanent bright objects pixels were excluded in ToPeCAl-2 by implementing bright object masking, it is still possible to have some of them generate false positive pixels in ToPeCAl-2 s results. Undetected bright objects were potentially due to the mixture of reflectance with other dominant objects (e.g., soil and vegetation in one-pixel area) and some of the rooftops had lower SWIR reflectance than the threshold used (<0.18). Most of the mining areas identified as false positive pixels also exhibited mixed reflectance with some dark soils, thus resulting in lower reflectance than the threshold for bright object masking. Some miss-detection of mining areas was possible due to their typical landscape which often consists of open-cut quarries which generate deep shadows leading to a decrease in reflectance. However, the state of the particular mining landscapes in this research requires further checking with very high spatial resolution data (for example Geo-eye, World-View 4 or SkySat). The removal of false positives related to mining areas and industrial gas flares can also be assisted via the use of other data sources, such as the record of location of registered mining areas and industrial gas flares from government agencies.

Relationships between TOA Reflectance of Sentinel-2 and Landsat-8 at Aerosol and SWIR bands
There was a significant difference of the mean values between TOA reflectance of Sentinel-2 and Landsat-8 at aerosol, SWIR-1, and SWIR-2 bands in the study areas. The mean differences of TOA reflectance from the two sensors at SWIR bands found in this study were slightly higher than the mean difference resulting from previous research in the Conterminous United States and southern Africa [41,45]. There was a lack of consistent relationship found in aerosol band which was also reflected at SWIR-1 and SWIR-2 bands. Conversely, good relationships of TOA reflectance at SWIR-1 and SWIR-2 bands were associated with a good relationship at aerosol band at other locations and times. The data resampling from lower to higher spatial resolution applied to aerosol band of Sentinel-2 and Landsat-8 could possibly influence the poorer relationship between both sensors at the aerosol band. Moreover, in this study sensor to sensor relationship analysis was conducted based on a linear empirical model, whereas the pre-launch relative spectral response of Landsat-8 and Sentinel-2, at mid spatial resolution, particularly in the aerosol band (Figure 8), shows obvious sensitivity differences of signal across the wavelength range. The atmospheric conditions, e.g., moist or dry, also influence the spectral response of SWIR [60]. The variability in coefficients from OLS modelling, suggests that corrections from TOA reflectance of Sentinel-2 to Landsat-8 should not be applied in ToPeCAl, as the modelling exhibits too much uncertainty. However, the subtle differences in band sensitivities between Sentinel-2 and Landsat-8 might be significant for other applications.
The atmospheric conditions in the tropics can be an important factor which influences the relationship of TOA reflectance between Sentinel-2 and Landsat-8, although the environments, season, and geographic location may also affect the relationships. Another aspect to consider is the cloud product of Sentinel-2 data which is generated from Sen2COR [65]. This might generate uncertainty in providing a clear image due to a lack of thermal data that can impact the approach and efficacy of cloud masking [61,69], particularly in tropical atmospheres. The influence of terrestrial environments, such as vegetation types, seasons, and an improvement of cloud masking product in tropical environments can be investigated further for relating the TOA reflectance of Sentinel-2 and Landsat-8.

The Application of ToPeCAl-2 to Sentinel-2 Data
The results of ToPeCAl-2 were not sensitive to the effects of TOA reflectance correction from Sentinel-2 to Landsat-8. Thus, ToPeCAl-2 with no transformation value can be applied to Sentinel-2 MSI. General results for peatland fires pattern resultant from ToPeCAl-2 applied to Landsat-8 were also found in ToPeCAl-2 applied to Sentinel-2. ToPeCAl-2 applied to Sentinel-2 with cloud masking provided higher detection for overall peatland fires with POD of 88%, ICE of 21% and IOE of 12%. Meanwhile, ToPeCAl-2 with a contextual test provided lower POD (62%), but with lower ICE (10%) and higher IOE (38%). The results also showed that smouldering fires were less likely to be detected from ToPeCAl-2 applied to Sentinel-2 data, but flaming pixels were well detected. Although ToPeCAl-2 with cloud masking resulted lower IOE than ToPeCAl-2 with contextual test, especially for smouldering pixels, validation still needs to be investigated for those smouldering pixels which may contain false detection in burnt areas. Meanwhile, the false positive pixels from ToPeCAL-2 might be the result from another fire at a different time of observation or un-detected bright object pixels. The POD values from the comparison of two fires produced from Sentinel-2 and Landsat-8 were presented in order to illustrate the strong similarity of fire pattern and the existence of fire for a certain time difference. Thus, they can be utilised together to fill temporal gaps in detailed fire spatial information.
The choice of using ToPeCAl-2 with contextual test or ToPeCAl-2 with cloud masking is dependent on the purpose for obtaining the peatland fires information. Users who are dealing with fire extinguishment may utilise ToPeCAl-2 fires from contextual test which provided high confidence of flaming detection. For users who work in peatland rehabilitation or carbon loss inventory from peatland fires, then the fire product from ToPeCAl-2 with cloud masking can be considered. However, smouldering resulting from ToPeCAl-2 especially from cloud masked approach still needs further validation. One should also consider the advantages and the limitation of both ToPeCAl-2 models.
The application of ToPeCAl without TIR to Sentinel-2 MSI and relying on SWIR to detect fires in tropical peatland environments is possible, but may present with false positives related to bright objects and warm soil due to the contamination of sun radiation on SWIR reflectance during daytime.
Smouldering peatlands can also be confirmed by observation of smoke over smouldering surfaces; a process which can be explored in further research. Overall, without MIR or TIR, ToPeCAl has a weaker detectability of smouldering peatland which is important for carbon emission estimations. This constraint introduces a future challenge to provide a sensor, or constellation of sensors incorporating MIR/TIR with similar, or better, temporal and spatial resolution specifications to Sentinel-2 to improve peatland fire monitoring; for example the Copernicus High Spatio-Temporal Resolution Land Surface Temperature Mission -LSTM mission which is planned by European Space Agency for future earth observation in 2025 [70].

Conclusions
A new technique for peatland fires detection in tropical peatland environments has been established using Landsat-8 OLI and Sentinel-2 MSI data in this study. ToPeCAl-2 (ToPeCAl without TIR) applied to Sentinel-2 MSI generates reliable results for fire detection at 20 m resolution especially for flaming and mix of flaming and smouldering with less smouldering detection. Furthermore, it can be combined with the use of ToPeCAl-1 to fill in the peatland fires information gaps at Landsat-8 spatial resolution (30 m) for tropical environment conditions, e.g., Indonesia. ToPeCAl-2 and ToPeCAl-1 has similar capability on detecting flaming fires (>95%) in tropical peatland environments. Smouldering peatlands detected by ToPeCAl-2 were only 35-56% correct, which still may include false detection in warm soil or burnt areas. ToPeCAl-2 with contextual test provides less detection of smouldering pixel than ToPeCAl-2 with cloud-masked approach. To further extend these smouldering results, especially from ToPeCAl with cloud-masked, requires further investigation and ground validation.
Adjustments to eliminate the false positive fires pixels were applied to ToPeCAl-2 to enhance its performance. It is important for tropical peatland fires management, with limited resources, to support a priority to extinguish actual biomass fires, which extend over large areas, rather than send resources to incorrectly identified fires.
In relating peatland fires products from Sentinel-2 and Landsat-8 at different times certain aspects were also examined, including the geometric errors between the two sensors, spectral response, and the propagation of fires between images captured at different times. The subtle differences in the spectral response of Landsat-8 and Sentinel-2 were less significant for fire detection, possibly due to specific local atmospheric conditions which can interfere with TOA reflectance in tropical environments.
There is a high potential to apply ToPeCAl-2 to Sentinel-2 MSI for operational peatland fires monitoring especially for flaming conditions, together with ToPeCAl-1 applied to Landsat-8 OLI/TIRS, leading to a higher temporal monitoring frequency than currently available from a single sensor alone. An integration with lower spatial resolution fire products from VIIRS 375 m data will also complement peatland fires monitoring on daily basis. Table A3. Landsat-8 OLI TOA reflectance as a function of Sentinel-2 MSI TOA reflectance at aerosol, SWIR-1, and SWIR-2 bands. Intercept and slope were derived from the ordinary least squares (OLS) regression of total data (n). The coefficients of determination (R 2 ) were significant at p-values < 0.001. The data sites were Central Kalimantan and South Kalimantan (CK-SK), Jambi and South Sumatera (JSS), and Riau.     F  6301  4101  1892  2209  2817  2319  498  95  17  5  FS  3195  19,150  9852  9298  5234  3299  1935  89  36  11  S  ----25,