New Approach Evaluating Peatland Fires in Indonesian Factors

: Until 2018, the El Niño–Southern Oscillation (ENSO) was used as an explanation for ﬁres in Indonesia’s peatlands. However, when the 2019 ﬁres occurred independently of El Niño, more suitable indicators and methods were required to (a) analyze, (b) evaluate and (c) forecast peatland ﬁres. In this study, we introduced the OLR–MC index—one of the rain-related indices derived from OLR (outgoing longwave radiation) in MC (maritime continent) area in Indonesia. This index showed stronger correlation with active peatland ﬁres than the conventional ENSO index, and is likely to be able to respond to heat and dry weather supposed to be under climate-change conditions. We then analyzed peatland ﬁres in the top six ﬁre years from 2002 to 2018 and showed that peatland ﬁres occurred in three stages—surface ﬁre, shallow peatland ﬁre and deep peatland ﬁre. To explain each stage, we proposed a one-dimensional groundwater level (GWL) prediction model (named as MODEL-0). MODEL-0 predicts GWL from daily rainfall. Analysis using MODEL-0 showed the GWL thresholds for the three ﬁre stages were between -300 mm and -500 mm; peatland ﬁre activities during the three ﬁre stages were dependent on these GWL values. The validity of MODEL-0 was shown by comparison with the measured values of GWL in the top three ﬁre years.


Introduction
Indonesia has some of the world's highest rates of deforestation and forest degradation, the principle drivers of which are agricultural expansion and wood extraction in combination with an increased incidence of fire [1]. An average annual clearance rate of 10 × 10 3 km 2 during the 1980s has increased to an average of 20 × 10 3 km 2 per year since 1996. The largest tree-cover loss since 2001 was 24.2 × 10 3 km 2 in 2016 [2].
As fire has been employed as a cheap land-clearance tool, uncontrolled fires from various sources have entered adjacent forests or plantation estates, and burnt both the forest biomass and, in peatland areas, underlying peat [1]. Analysis of recent land-cover change has underscored the close link between land management and fire activity [3]. Most burn scars were located along forest edges (i.e., disturbed forest), usually in close proximity to human settlements that provided a source of ignition. More remote, intact forests were unaffected by fire, even during extended El Niño-Southern Oscillation (ENSO) drought [1].
ENSO events or droughts also play a role in exacerbating fire occurrence and severity [1]. The worst fires were associated with strong ENSO episodes, including the fires of 1972-73, 1982-83 (the 'Great Fire of Borneo'), 1994, 1997-98, 2002 and 2006. The 1982-83 fires were estimated to have burnt 32 × 10 3 km 2 , The risk of fire on peatland is increased greatly by drainage, which lowers the water table, exposing a greater volume of dry peat to combustion. This effect is demonstrated by a study of the fire regime in the former mega rice project (MRP) area on peatland in southern Central Kalimantan [11]. From 1973 to 1996, fires affected 24% of a 4.5 × 10 3 -km 2 Block C area in MRP. Following multiple fires in peat swamp forest, the numbers of tree species and individual trees, saplings and seedlings within the secondary vegetation are greatly reduced and, at the highest levels of degradation, succession back to forest is diverted to a retrogressive succession to communities dominated by ferns with very few or no trees [12,13]. Non-forested areas in Kalimantan had shorter average fire return intervals (FRI) than Sumatra (13 years versus 40 years), with ferns/low shrub areas burning most frequently [14].
Peatland fires make a major contribution to emissions of greenhouse gases, fine particulate matter and aerosols, thus contributing to climate change as well as presenting a problem for human health [1]. The devastating 1997-1998 Indonesian fires were among the largest peak emission events in the recorded history of fires in equatorial Southeast Asia [15,16]. According to a recent report on intense fires in 2015, peat fires in Indonesia during July-October 2015 released about 2 Tg of carbon into the atmosphere, 81% of which was in the form of carbon dioxide (CO 2 ); 16% carbon monoxide (CO) and 2.3% methane (CH 4 ) [17]. Fine particulate matter (PM2.5) emissions from fires across Sumatra and Kalimantan (Borneo) during September-October 2015 were estimated at 7.33 Tg [18]. The level of major air-pollution gases and particulate matter (PM10) at Palangkaraya near the mega rice project was already reported [19].
Smoldering wildfire in peatlands is one of the largest and longest-lasting fire phenomena on Earth [20]. Losses of carbon from above-ground biomass (AGB) usually greatly exceed those from soil, but in regions with carbon-rich substrates-e.g., peat swamp forests-combustion losses from below-ground stocks can be much greater than AGB losses [1]. For this reason, much peat research is carried out to control fires in peatland and to clarify the background conditions of multiple fires in peatland.
Comprehensive peat study clarified the combustion properties of peat, the thermal characteristics of peat fire, the temperature distribution in peat layers during fire, the pattern of peat fire fronts and the peat fire spreading speed [21].
Experiments on smoldering and flaming conditions of peat demonstrated the successful piloted flaming ignition of peat soil with moisture up to 100 wt% under external radiation, indicating that flames may spread rapidly on peatland before transitioning to a conventional smoldering peat fire [20].
The upward and downward spread of smoldering peat fire was analyzed using a one-dimensional computational model. Results showed that smoldering combustion can spread over peat layers with very high moisture content (MC) (MC > 250%) if the layer is thin and located below a thick, drier layer [22].
In this study, we discuss two issues related to peatland fires in Indonesia. First, we introduce one of existing indicators to explain active fires in 2019 occurred regardless of ENSO-neutral. The correlation between the OLR-MC and the fire occurrence is analyzed. The OLR-MC is one of indices issued by JMA and a rain-related index derived from outgoing longwave radiation (OLR) in Indonesia (MC: maritime continent, area from Kalimantan to west Papua).
Second, we propose a one-dimensional groundwater level (GWL) prediction model (MODEL-0) to evaluate peat fire activities. MODEL-0 needs daily rainfall as input data to calculate the diurnal change of GWL. Therefore, rainfall data from satellite observations, meteorological observatory data and estimated values from OLR are compared.
Finally, we analyze peatland fire activities during the top six fire years from 2002 to 2018 to find their active fire periods. To explain fire activities during the fire stages, we used MODEL-0 to identify the relationship with GWL. Figure 1 shows the peatland distribution and study area in Kalimantan (Borneo). Palangkaraya is the capital of Central Kalimantan is located in the southern part of the province (2.21 • S, 113.92 • E-hereafter, a simpler expression without the degree sign will be used, e.g., 2.21-2.23 S). Peatlands in Kalimantan cover about 57,600 km 2 and are equivalent in area to that of Sumatra. Central Kalimantan alone has 30,100 km 2 of peatlands. The distance from the nearest coastline is about 100 km and the average altitude is only around 10 m. The MRP was built on tropical swamp forest areas on the eastern and southern sides of Palangkaraya. Before the disturbance, the tropical swamp forest could hold enough water to stay wet even in the dry season. However, the newly constructed 4000-km-long MRP canal built for irrigation facilitated not only illegal logging, but also a loss of water through drainage from most of the peatlands in the MRP area. These disturbances are the main reasons for intense and recurrent fire activity in the MRP area.

Study Area
We refer to the study area as 'MRP+' in this report. MRP+ covers the MRP area and its vicinity (a part of the Sebangau National Park) (1.75-3.5 S, 113.5-115.0 E). MRP is divided into five blocks: Block A to Block E, as shown in Figure 1. Their actual boundaries are defined mainly by rivers. In this report, approximate boundaries defined by latitude and longitude lines are used for convenience. MRP+ was chosen as the study area simply because MRP+ area has major peatland in Central Kalimantan and is one of the highest hotspot density areas in Indonesia [23].   Information  for  Resource  Management System, https://firms2.modaps.eosdis.nasa.gov/download/(Accessed 6 June 2020)). We use only the spatial and temporal (latitude, longitude and acquisition date and time) hotspot data in this study. The number of daily hotspots is used to identify fire activities and the fire year.

OLR-MC Index
The OLR-MC (rain-related index derived from outgoing longwave radiation (OLR), MC: maritime continent) data were obtained from the Japan Meteorological Agency (https://www.data.jma.go.jp/gmd/cpd/db/diag/2018/index/html/soiolru/index_html_soiolru_2018.htm l (Accessed 6 June 2020)). The OLR-MC index monitors convective activity in Indonesia (5 N-5 S, 110-135 E) and used to identify a relationship with fire occurrence in Kalimantan. One study showed a high negative correlation between OLR and rainfall in South East Asia [9].

Rainfall Data
Some rainfall data are used in this report to compensate missing data, mainly in 2019. daily rainfall data were needed to estimate the daily change in groundwater level (GWL).
Generally speaking, Indonesian meteorological records lack data. In particular, there was much missing data on rainfall in 2019. But several Japanese researchers, including us, continue to collect rainfall data in collaboration with Indonesian researchers, enabling about 40 years (from 1978 to 2018) rainfall data to be recorded. Long-term data are useful for calculating average rainfall during dry and wet seasons, for example.

Hotspot (Fire) Data
Eighteen years of hotspot (HS, temperature anomaly) data detected by moderate resolution imaging spectroradiometer (MODIS) on the Terra and Aqua satellites are used to evaluate fires in boreal forests. MODIS HS data collected from 2002 to 2019 were obtained from NASA FIRMS (Fire Information for Resource Management System, https://firms2.modaps.eosdis.nasa.gov/download/ (Accessed 6 June 2020)). We use only the spatial and temporal (latitude, longitude and acquisition date and time) hotspot data in this study. The number of daily hotspots is used to identify fire activities and the fire year.

OLR-MC Index
The OLR-MC (rain-related index derived from outgoing longwave radiation (OLR), MC: maritime continent) data were obtained from the Japan Meteorological Agency (https://www.data.jma.go.jp/ gmd/cpd/db/diag/2018/index/html/soiolru/index_html_soiolru_2018.html (Accessed 6 June 2020)). The OLR-MC index monitors convective activity in Indonesia (5 N-5 S, 110-135 E) and used to identify a relationship with fire occurrence in Kalimantan. One study showed a high negative correlation between OLR and rainfall in South East Asia [9].

Rainfall Data
Some rainfall data are used in this report to compensate missing data, mainly in 2019. daily rainfall data were needed to estimate the daily change in groundwater level (GWL).
Generally speaking, Indonesian meteorological records lack data. In particular, there was much missing data on rainfall in 2019. But several Japanese researchers, including us, continue to collect rainfall data in collaboration with Indonesian researchers, enabling about 40 years (from 1978 to 2018) Remote Sens. 2020, 12, 2055 5 of 16 rainfall data to be recorded. Long-term data are useful for calculating average rainfall during dry and wet seasons, for example.
(2) Satellite rainfall data Daily rainfall data provided by Japan Aerospace Exploration Agency (JAXA) [26] were used to estimate groundwater level during the fire (dry) season. We downloaded rainfall data by setting a target rectangle region in Palangkaraya (see a black dotted line rectangle in Figure 1).
(3) Estimated rainfall using OLR Using the above-mentioned OLR-MC index, we analyzed the relationship between OLR values and rainfall data. The OLR value at point O (2.21 S, 114 E, 'O' in Figure 1) was read from the OLR contour map.
Daily values and daily long-term mean values were estimated from (3-1) NOAA daily (non-interpolated) OLR [27] and (3-2) NOAA Interpolated OLR [28,29], respectively. daily values and daily long-term mean values were estimated from a correlation equation between OLR values and weather rainfall data at TRMS.

Ground Water Level (GWL) Data
Ten-year (1997-2006) GWL measurements [30] were averaged to develop a GWL simulation model. GWL was measured near point U (2.32 S, 113.9 E, 'U' in Figure 1) on the upper catchment of the Sebangau River about 800 m from the forest edge and 2 km from the river along the logging railway on the right bank of the Sebangau River. The depth of peat near the measurement point is approximately 3 m [30]. Unfortunately, the measuring equipment was stolen in 2014, so the data location in 2015 and 2019 is near 'U' (ca. 0.9 km northeast of 'U').

ENSO Data
One of ENSO (El Niño-Southern Oscillation) indices, Niño 3.4, is obtained from the National Weather Service, Climate Prediction Center, NOAA [31].

Ground Water Level (GWL) Simulation (MODEL 0)
We used one-dimensional GWL model [32,33]. The water balance equation is expressed as [32]: where ∆S is the change in storage of groundwater, R is rainfall, Q is the sum of surface runoff, transpiration, interflow and surface detention, E is evapotranspiration. The daily change of GWL was calculated by the following formula: where ∆t = one day, Q = 0, E is assumed as 8 mm day −1 then, (GWL of the day) = (GWL of the previous day)+ R-8 where E (evapotranspiration) is assumed as 8 mm day −1 by considering rate of GWL decrease derived from GWL data (1997)(1998)(1999)(2000)(2001)(2002)(2003)(2004)(2005)(2006) [30]. We ignore the value of Q because all component of Q are soil surface processes (Q = 0) [32].

Annual Hotspot (Fire) Trend in Kalimantan
The total annual number of hotspots (HSs) in each of the 18 years (2002-2019) is summarized in Table 2, which is arranged in descending order of the annual number of HSs in the MRP+ region (vertical row named as '(1) HS MRP+'). The annual number of HSs detected in Kalimantan Island (including Malaysia and Brunei) is also added in vertical row named as '(2) HS Kalimantan' in Table 2. The difference in the number of HSs-(2) HS Kalimantan in Table 2-between 2014 (top 7th) and 2012 (top 8th)-is large or ca. 13,000. The top seven years from Rank 1 (2015) to Rank 7 (2014) are called 'fire year' in this report. The study area MRP+ is 32,400 km 2 in size, which accounts for about 6% of Kalimantan Island (544,150 km 2 ). Nevertheless, about 20% (see the ratio (1)/(2) in Table 2) of all Kalimantan fires have occurred in the study area. This suggests that peatland fires existing in MRP+ contribute greatly. The total number of top seven fires in MRP+ was 77,446, about 80% of the total number (97,423) for 18 years. Most fires in MRP+ occurred in the top seven fires.

OLR-MC Index
OLR-MC index is introduced as an index for predicting peatland fire occurrence in the study area, MRP+. Figure where HSs is the number of HSs, OLR is the OLR-MC index.
The relatively high value of decision coefficient (R 2 = 0.84) suggests that there was a good correlation between the number of HSs and the OLR-MC index over the last 18 years, including 2019.

Niño 3.4 Index
The most referenced fire index in Indonesia is Niño 3.4 index, one of the ENSO indices. Niño 3.4 is introduced here to identify the relationship with the number of HSs in study area MRP+.

Niño 3.4 Index
The most referenced fire index in Indonesia is Niño 3.4 index, one of the ENSO indices. Niño 3.4 is introduced here to identify the relationship with the number of HSs in study area MRP+.
where HSs is the number of HSs, and Niño is the Niño 3.4 index. The number of observed HSs and estimated HSs using Equation (5) where HSs is the number of HSs, and Niño is the Niño 3.4 index. The number of observed HSs and estimated HSs using Equation (5) (Niño 3.4 index = −0.04) in 2019 were 7334 and 2242, respectively. This large underestimated HS value and the relatively low value of decision coefficient (R 2 = 0.60)-compared with that of (R 2 = 0.84) for OLR-MC index-suggests that the Niño 3.4 index will work for fire occurrence related to ENSO, but will not be suitable for fires like those in 2019 supposed to be under climate-change conditions.
The OLR-MC index shows a stronger correlation with active peatland fires than the conventional ENSO index, and is likely to be able to respond to future fire forecasts under climate-change conditions.

Comparison of OLR-MC and Niño 3.4
The relationship between the number of HSs and OLR-MC index is shown in Figure 4a. Figure 4a shows a good fitting of the regression equation (Equation (4) The relationship between the number of HSs and OLR-MC index is shown in Figure 4a. Figure  4a shows a good fitting of the regression equation (Equation (4)

Estimated Rainfall using OLR
The OLR-MC index already showed a good correlation with HS in Section 4.2.1. A report [34] showed there are negative correlations between OLR and rainfall over equatorial regions (0-10 S). Therefore, we compared the daily average rainfall data from 1978 to 2018 measured at TRMS ('A' in

Estimated Rainfall using OLR
The OLR-MC index already showed a good correlation with HS in Section 4.2.1. A report [34] showed there are negative correlations between OLR and rainfall over equatorial regions (0-10 S). Therefore, we compared the daily average rainfall data from 1978 to 2018 measured at TRMS ('A' in  Figure 5 shows the relationship between averaged daily rainfall data (TRMS) and smoothed OLR. A total of 161 data (DN =170-330) were plotted in Figure 5 daily rainfall was estimated using the linear correlation equation: daily rainfall (mm) = 43.881 − 0.16986 OLR (R 2 = 0.52965) (6) obtained from the OLR contour map. OLR values were smoothed by the 5-day running mean values. Figure 5 shows the relationship between averaged daily rainfall data (TRMS) and smoothed OLR. A total of 161 data (DN =170-330) were plotted in Figure 5 daily rainfall was estimated using the linear correlation equation: daily rainfall (mm) = 43.881 − 0.16986 OLR (R 2 = 0.52965) (6)

Rainfall during Dry Season
Rainfall trend during dry season in Palangkaraya are shown in Figure 6. The dry season usually begins in June and ends in October, with the rainy season from November to May. As daily rainfall tends to decrease gradually from June 19-Day Number (DN)=170)-we make two accumulated rainfall curves, one for the 41 years average rainfall from 1978 to 2018 and another for the top six fire years excluding 2019. Rainfall trend during dry season in Palangkaraya are shown in Figure 6. The dry season usually begins in June and ends in October, with the rainy season from November to May. As daily rainfall tends to decrease gradually from June 19 -Day Number (DN)=170) -we make two accumulated rainfall curves, one for the 41 years average rainfall from 1978 to 2018 and another for the top six fire years excluding 2019. Over those 41 years, the average daily rainfall was 4 mm day −1 in the dry season, 10 mm day −1 in the rainy season and 5 mm day −1 in the transition period from DN = 170. The 30-year average rainfall estimated by OLR shows a similar trend with weather station data. As a result of this comparison, we may say the estimated rainfall using OLR can be applied to various analyses because the margin for error is not so large.
Two sets of average rainfall data of the top six fire years, excluding 2019, are also shown in Figure  6. The 41 years average daily rainfall for the top six fire years is 2 mm day −1 in the dry season, 10 mm day −1 in the rainy season and 2.4 mm day −1 in the transition period. Rainfall for the top six fire years in the dry season is half that of 41-year average. Figure 6 clearly shows severe drought conditions. The JAXA-averaged rainfall for the top six fire years also shows a similar trend with weather station data. From this comparison, we may say the JAXA rainfall data will be applicable for various analyses because the margin for error is not so large.  Over those 41 years, the average daily rainfall was 4 mm day −1 in the dry season, 10 mm day −1 in the rainy season and 5 mm day −1 in the transition period from DN = 170. The 30-year average rainfall estimated by OLR shows a similar trend with weather station data. As a result of this comparison, we may say the estimated rainfall using OLR can be applied to various analyses because the margin for error is not so large.

Fire Occurrence Trends in Fire Years
Two sets of average rainfall data of the top six fire years, excluding 2019, are also shown in Figure 6. The 41 years average daily rainfall for the top six fire years is 2 mm day −1 in the dry season, 10 mm day −1 in the rainy season and 2.4 mm day −1 in the transition period. Rainfall for the top six fire years in the dry season is half that of 41-year average. Figure 6 clearly shows severe drought conditions. The JAXA-averaged rainfall for the top six fire years also shows a similar trend with weather station data. From this comparison, we may say the JAXA rainfall data will be applicable for various analyses because the margin for error is not so large.  Figure 7, the fire occurrence tendency of each fire year can be grasped. In Figure 7, straight lines are drawn on each accumulated hotspot (HS) curves of the fire year. These lines are useful to find fire activity, active fire periods (with start and end dates) and duration. Fire rate (HS day −1 ) derived from the gradient of the straight line is used to evaluate fire activities.  Figure 8 shows the average fire occurrence trend of the fire year from June 20 (DN = 170) and estimated GWL using MODEL-0. Additional straight lines on an average accumulated HS curve are drawn to define several fire stages. These stages are inherently derived from fire occurrence trends on the average accumulated HS curve. Here, we could define a few fire stages by gradients of several straight lines and named them Stage 1, Stage 2 and Stage 3, which are shown in Figure 8. We follow the stage numbering (Stage 1, Stage 2 and Stage 3) in the previous report [21]. Fires in Stage 1 (fire rate ca. 10 HS day −1 ) start from the middle of July and end in late July. Then, fires in Stage 2 (fire rate ca. 70 HS day −1 ) begin from the middle of August and last until late August. Finally, most active fires in Stage 3 (fire rate ca. 160 HS day −1 ) occur from the beginning of September and last until the end of October. However, the fire will continue until early November when the full rainy season begins.  Figure 8. We follow the stage numbering (Stage 1, Stage 2 and Stage 3) in the previous report [21].

Fire Occurrence Trends in Fire Years
Fires in Stage 1 (fire rate ca. 10 HS day −1 ) start from the middle of July and end in late July. Then, fires in Stage 2 (fire rate ca. 70 HS day −1 ) begin from the middle of August and last until late August. Finally, most active fires in Stage 3 (fire rate ca. 160 HS day −1 ) occur from the beginning of September and last until the end of October. However, the fire will continue until early November when the full rainy season begins.
Estimated GWL using MODEL-0 in Figure 8 was obtained by using average rainfall of top fire year (TRMS). Equation (3) ((GWL of the day) = (GWL of the previous day) + R-8) is used to estimate daily GWL. Calculation starts from June 19. Initial GWL on June 19 is set to 0 mm from averaged groundwater level [30]. GWL of next day (June 20) is +3.3 mm (=0 + R-8 (where R (rainfall on June 20) is 11.3 mm from average TRMS rainfall data). Like this manner, daily GWL is estimated. GWL becomes under −500 mm on September 1 (DN = 244). (MODEL-0) 4.5.1. Simulation (MODEL-0) Results for Fire Years Figure 8 shows the average fire occurrence trend of the fire year from June 20 (DN = 170) and estimated GWL using MODEL-0. Additional straight lines on an average accumulated HS curve are drawn to define several fire stages. These stages are inherently derived from fire occurrence trends on the average accumulated HS curve. Here, we could define a few fire stages by gradients of several straight lines and named them Stage 1, Stage 2 and Stage 3, which are shown in Figure 8. We follow the stage numbering (Stage 1, Stage 2 and Stage 3) in the previous report [21]. Fires in Stage 1 (fire rate ca. 10 HS day −1 ) start from the middle of July and end in late July. Then, fires in Stage 2 (fire rate ca. 70 HS day −1 ) begin from the middle of August and last until late August. Finally, most active fires in Stage 3 (fire rate ca. 160 HS day −1 ) occur from the beginning of September and last until the end of October. However, the fire will continue until early November when the full rainy season begins.

Hotspot (Fire) Simulation
Estimated GWL using MODEL-0 in Figure 8 was obtained by using average rainfall of top fire year (TRMS). Equation  The fire activity of the three fire stages was dependent on the GWL depth, with thresholds of ca. −300 mm and ca. -500 mm. The three fire stages are classified into a surface fire, a shallow peat fire and a deep peatland fire [21]. This is because the lower the GWL, the lower the moisture content of the peat and the more easily the peat can be burned [35,36]. The onset of peatland fire is confirmed by measuring air pollutant gases and particulate matter from peat burning. In fact, increasing trends in carbon monoxide (CO), sulfur dioxide (SO 2 ), and particulate matter (PM10) were observed from mid-August 2002 [19]. Increased SO 2 from around 5 x10 −6 g m −3 in early August to 85.8 × 10 −6 g m −3 in mid-October was thought to be primarily due to peat fires, as there are no other major sources of SO 2 (such as industry) nearby. CO and PM10 also began to increase in early August. These increases are due to the low-temperature, underground peatland fire (smoldering combustion under low levels of oxygen). The increase in CO, SO 2 and PM10 at the same time can be considered one of the characteristics of peatland fires. Figure 9 shows the relationship between fire trends and GWL in 2015 (the top fire year, HSs = 19,936, Niño 3.4 = 2.12 (very strong), OLR-MC = −1.9). In Figure 9, GWL measured values are plotted in order to evaluate the validity of MODEL-0 (GWL estimation model). From Figure 9, we could see the following: Fires in 2015 became active (ca. 160 HS day −1 ) from mid-August when the GWL was below −300 mm. In early September, after the GWL fell below −500 mm, the fire became very intense (ca. 370 HS day −1 ) and continued until late October. The fires in 2015 were weakening due to the rainy and hazy conditions in September, but the fires continued because of the low GWL (under −500 mm). Figure 10 shows the relationship between fire trends and GWL in 2006 (the top second fire year, HSs = 11,754, Niño 3.4 = 0.66 (weak), OLR-MC = −1.7). In Figure 10, GWL measured values are plotted in order to evaluate the validity of MODEL-0.

Simulation (MODEL-0) Results for 2015 (top) & 2006 (2nd)
From Figure 10, we could see the following: Fires in 2006 became active (ca. 70 HS day −1 ) from mid-August when the GWL was below −300 mm. In mid-September, after the GWL fell below −500 mm, the fire became very intense (ca. 240 HS day −1 ) and continued until mid-November. The fires in 2006 were weakening due to the rainfalls of 19 mm and 9.6 mm in September and October, but the fires continued because of the low GWL (under −500 mm).
19,936, Niño 3.4 = 2.12 (very strong), OLR-MC = −1.9). In Figure 9, GWL measured values are plotted in order to evaluate the validity of MODEL-0 (GWL estimation model). From Figure 9, we could see the following: Fires in 2015 became active (ca. 160 HS day −1 ) from mid-August when the GWL was below −300 mm. In early September, after the GWL fell below −500 mm, the fire became very intense (ca. 370 HS day −1 ) and continued until late October. The fires in 2015 were weakening due to the rainy and hazy conditions in September, but the fires continued because of the low GWL (under −500 mm).   From Figure 10, we could see the following: Fires in 2006 became active (ca. 70 HS day −1 ) from mid-August when the GWL was below −300 mm. In mid-September, after the GWL fell below −500 mm, the fire became very intense (ca. 240 HS day −1 ) and continued until mid-November. The fires in 2006 were weakening due to the rainfalls of 19 mm and 9.6 mm in September and October, but the fires continued because of the low GWL (under −500 mm).  Figure 11, GWL measured values are plotted in order to evaluate the validity of MODEL-0. Here, we adjusted the GWL measured for comparison with estimated GWL. The GWL measured was high until late July, probably because of the river's high water level.
From Figure 11, we can see the following: Fires in 2019 became active (ca. 40 HS day −1 ) from late July when the GWL was below −300 mm. In late August, after the GWL fell below −500 mm, the fire became very intense (ca. 270 HS day −1 ) and continued until mid-November. The large difference between GWL measured values and the estimated GWL is mainly due to a different measuring position. However, the tendency of the GWL decrease is almost the same. The peatland fires in 2019 weakened due to the heavy rainfall of 44 mm in August, but they continued because the rain could not sufficiently moisten the peat (GWL was around −450 mm). A few days after the heavy rain, GWL began to decline again, and fires were very active when GWL fell to ca. −500 mm. The 2019 case demonstrates how effective MODEL-0 was in forecasting GWL. In particular, it is very important to predict the occurrence of deep peatland fires (fire stage 3), and MODEL-0 could do this for the 2019 case, if applied.  Figure 11 shows the relationship between fire trends and GWL in 2019 (the third top fire year, non-El Niño year, HSs = 11,070, Niño 3.4 = −0.04 (neutral), OLR-MC = −1.9). In Figure 11, GWL measured values are plotted in order to evaluate the validity of MODEL-0. Here, we adjusted the GWL measured for comparison with estimated GWL. The GWL measured was high until late July, probably because of the river's high water level.
From Figure 11, we can see the following: Fires in 2019 became active (ca. 40 HS day −1 ) from late July when the GWL was below −300 mm. In late August, after the GWL fell below −500 mm, the fire became very intense (ca. 270 HS day −1 ) and continued until mid-November. The large difference between GWL measured values and the estimated GWL is mainly due to a different measuring position. However, the tendency of the GWL decrease is almost the same. The peatland fires in 2019 weakened due to the heavy rainfall of 44 mm in August, but they continued because the rain could not sufficiently moisten the peat (GWL was around −450 mm). A few days after the heavy rain, GWL began to decline again, and fires were very active when GWL fell to ca. −500 mm. The 2019 case demonstrates how effective MODEL-0 was in forecasting GWL. In particular, it is very important to predict the occurrence of deep peatland fires (fire stage 3), and MODEL-0 could do this for the 2019 case, if applied. weakened due to the heavy rainfall of 44 mm in August, but they continued because the rain could not sufficiently moisten the peat (GWL was around −450 mm). A few days after the heavy rain, GWL began to decline again, and fires were very active when GWL fell to ca. −500 mm. The 2019 case demonstrates how effective MODEL-0 was in forecasting GWL. In particular, it is very important to predict the occurrence of deep peatland fires (fire stage 3), and MODEL-0 could do this for the 2019 case, if applied.

Discussion
Many researchers have revealed peat-burning characteristics as well as the relationship with groundwater levels, carbon dioxide emissions, air pollution and ENSO. However, little research was carried out on: (1) predictive indicators other than ENSO for active peatland fires, and (2) the groundwater level (GWL) prediction model (MODEL-0) to evaluate the onset of active peatland fire, as described in this report.

Predictive Indicators Other than ENSO Indices for Active Peatland Fires
Until 2019, most peatland fires were associated with an El Niño event. Especially active peatland fires in 1982-83, 1997-98 and 2015 occurred under very strong El Niño conditions [1] However, active peatland fires broke out in 2019, unrelated to the El Niño phenomenon (see Figures 2 and 3) and could be linked to climate change. Preliminarily, we introduced the OLR-MC index (see 2.2.2) and examined the relationship with the number of HSs and OLR-MC index readings. As a result, the regression equation between fires (the number of HSs) and OLR-MC index readings showed a relatively high decision coefficient (R 2 = 0.84). This suggested that there was a good correlation between the number of HSs and the OLR-MC index over the past 18 years, including 2019 (see Figure 2, Equation (4)).
On the other hand, a regression equation between fires (the number of HSs) and Niño 3.4 (one of the ENSO indices) had a relatively low decision coefficient (R 2 = 0.60). In addition, there was a large difference between the number of observed HSs and estimated HSs from the regression equation regarding Niño 3.4 (index = −0.04 in 2019). The values were 7334 and 2242, respectively (see Figure 3, Equation (5)).
From the above results, we may conclude that an OLR index, such as OLR-MC, will be useful for analyzing future active fires related to both El Niño and climate change. In addition, this report compared various rainfall data. Results showed that the remote sensing and satellite data provided by JAXA and the estimated data using OLR are effective for rainfall analysis, such as that presented in this report.

Groundwater Level (GWL) Prediction Model ((MODEL-0))
Peatland fire could be classified into three stages, namely surface fire, surface peatland fire and deep peatland fire [21]. In this report, the authors proposed a one-dimensional groundwater level (GWL) prediction model (MODEL-0) to clarify the relationship between peatland fire activity and GWL. MODEL-0 calculates daily changes in GWL using daily precipitation (for simplicity, dry season evapotranspiration is assumed as constant). GWL is important for knowing the dry state of peat.
From the average number of fire hotspots in the top six fire years (excluding 2019) and the daily change in GWL trend calculated from the average daily rainfall, peatland fires are divided into three stages as GWL penetrates deeper. The fire activity of the three stages was dependent on the GWL depth, with their threshold values being ca. −300 mm and ca. −500 mm. Thus, MODEL-0 clarified the relationship between the fire activity of the three fire stages and GWL (see Figure 8).
Peatland fire analysis using MODEL-0 was successfully applied to the top three peatland fires in 2005, 2006 and 2019 (see Figures 9-11). MODEL-0 was able to estimate the groundwater level relatively well and could explain the active peatland fire occurrence depending on the groundwater level. In particular, MODEL-0 was also effective in the 2019 fire analysis, which was a neutral El Niño year, so it can be said that the possibility of responding to extreme weather supposed to be under climate change conditions was demonstrated.

Conclusions
We analyzed recent active peatland fires in the MRP+ region in Central Kalimantan from 2002 to 2019. The results allowed us to make the following conclusions: (1) OLR-MC (monitoring convective activity in Indonesia using a rain-related index derived from the OLR) index is introduced as one of the fire indices for predicting peatland fire occurrence. The OLR-MC index showed a stronger correlation with active peatland fires than the conventional ENSO index and is likely to be able to respond to future fire forecasts supposed to be under climate change conditions (Figures 2 and 3).  Figure 11). The 2019 case demonstrated how effective MODEL-0 was in forecasting GWL. In other words, the peatland fire situation (fire stages 1-3) can be ascertained by the GWL value derived from MODEL-0. In particular, it is important to predict the occurrence of deep peatland fires (fire stage 3) and MODEL-0 could predict it. (7) A method for estimating rainfall using the OLR distribution map was developed. Rainfall at any arbitrary location can be estimated with this method. MODEL-0 can be applied in any target region to determine the fire stages from the estimated GWL ( Figures 5 and 6). Funding: This research is partially funded by Ministry of Agriculture, Forestry and Fisheries of Japan grant to CIFOR.: Project named "Enhancing climate-resilient livelihoods in boreal and tropical high carbon forests and peatlands".