Building capacity for estimating fire emissions from tropical peatlands; a worked example from Indonesia

Tropical peatlands are globally significant in the terrestrial carbon cycle as they are comprised of a large forest carbon sink and a large peat carbon store—both of which can potentially be exchanged with the atmosphere on decadal time frames. Greenhouse gas emissions from fire-disturbance and development of tropical peatlands over the last few decades, and the potential for ongoing emissions, highlights the need for policy to slow or halt emissions and to activate mechanisms to sequester carbon through restoration of degraded peatlands. The UN REDD + scheme provides a means for developing countries to receive payments for avoided deforestation and forest degradation, but the steps to achieve REDD+ compliance are rigorous and the details required can be a barrier to activating benefits—especially for peatlands where repeated cycles of fire interrupt forest recovery and create a range of recovery classes. Therefore, to improve estimates of peat fire emissions and of carbon balance of tropical peatlands, the biomass and combustion factor parameters need to be developed and applied according to forest recovery stage. In this study we use published activity data from the extensive 1997 fires in the peatlands of Indonesian Borneo to detail a transparent and accountable way to estimate and report emissions from tropical peatland fires. This example for estimating and reporting emissions is provided to assist REDD+ countries to efficiently develop their capacity for improving emissions estimates from fire-impacted tropical peatlands.

gained national and international significance as a mechanism for addressing climate change.For these reasons the Food and Agriculture Organisation (FAO) of the United Nations (UN) has recently declared that improving the assessment of GHG emissions from peatlands is a global strategic priority 8 .Improving the assessment of GHG emissions from peatlands will allow for the development of policy that establishes incentives to achieve peatland restoration at a country scale, particularly through schemes such as REDD + .Several studies over the last two decades have documented the rapid degradation of tropical peatlands, with clearing and drainage often followed by frequent fires and development of agriculture and plantations over large areas of former peat swamp forest 9,10 .The opportunity for conserving and restoring tropical peatlands is a priority for climate change mitigation as well as biodiversity conservation 11 and water resource protection.
International pressure on Indonesia to reduce emissions from peat fires is ongoing, yet the magnitude of emissions released to the atmosphere is uncertain.This uncertainty is partly due to different approaches to account for emissions from the diversity of peatland landscapes that have evolved from the initial fire occurrence, that may or may not have followed initial logging.Peat swamp forests (PSF) degraded by logging and fire subsequently become either secondary forests, or plantations on peat, or repeatedly fire-impacted degraded peatlands.Intact or primary peat swamp forests store more carbon aboveground and in the peat layer than burnt and recovering forests or repeatedly burnt and degraded forest (Fig. 1), where the latter are generally dominated by shrubs and grasses and have a shallower peat layer than less disturbed forests 12 .
Recognizing the varying outcomes of either recovery or continued degradation of peatland landscapes following disturbance is important for carbon and emission accounting, as forest condition sets the starting point for selecting the most appropriate parameters for emissions calculations.However most published studies rely on the default parameters provided in the Guideline of the Intergovernmental Panel on Climate Change IPCC 13 that apply to primary forests, and rarely adopt parameters that are more appropriate to the burning of already degraded peat swamp forests.Moreover, as we mentioned above, highly cited publications on the amount of emissions released from peat fires in Indonesia have applied simplified, non-IPCC methods for estimating emissions, and have also not distinguished between primary and already degraded peatlands in making their emissions estimates 7,14 .
Initial methodologies for estimating emissions from biomass burning were published by the IPCC in 1996 15 .The parameters required for estimating emissions are the activity data (AD, or area burnt) and emission factors (EF)-Eq.( 1).
The EF is estimated as the mass of the fuel available for combustion (M), multiplied by the fuel combustion factor (Cf) and a gas specific emission factor (G i ).Cf refers to how much fuel is combusted (or released to the atmosphere) and G i refers to the amount of i greenhouse gas (i.e., CO, CH 4 , N 2 O) emitted per kg of fuel burnt-Eq.(2).
In the case of peatlands, separate EFs are required for aboveground and peat layers-as both are combusted in fires (Fig. 1), thus information on mass of peat (M PEAT ) and AGB and their corresponding Cf and Gi will be required for the estimates.
Mass of peat ( M PEAT ) combusted in fires is estimated as peat bulk density (BD) multiplied by the depth of burnt peat (h)-Eq.( 3).
(1) Emissions are reported in CO 2 -equivalent and global warming potential (GWP) for each of the reported gases (CH 4 , N 2 O, CO) is also required.
Thus, to estimate emissions for i gas from peat, we need: We also need to estimate emissions for i gas from Above Ground Biomass (AGB): Thus, the total emissions will be the sum of emissions from peat and AGB for all reported i gases (e.g., CO 2 , CO, N 2 O, CH 4 )-Eq.( 6): As we mentioned above, using the same M for primary or intact forests and for degraded peatlands will either over or underestimate emissions (Fig. 1).Therefore, emissions estimates should be stratified by peatland cover classes (or fire frequencies)-Eq.(7).
As research advances and as technologies progress revision of Eq. ( 4) and (5) parameters i.e., AGB, BD, h, Gi, Cf, acts to improve emissions estimates and reduce uncertainties.One good example of continuous updates in the emission equation parameters is for peat combustion factor (Cf) of Eq. ( 4).The widely used IPCC 2014 Supplement 13 states on page 2.39: "For all organic soil fires, the default combustion factor is 1.0, since the assumption is that all fuel is combusted (Yokelson et al., 1997)".Thus a single study is used for the peat related emissions estimates worldwide.However, on a detailed examination of the study referred too 16 , there is no mention of a Cf .Thus a great number of peat emissions estimates in reports and publications around the world are based on a default Cf that is not supported by the literature.The IPCC has been alerted to this oversight and is considering an update to default parameters for calculating emissions from peat fires.The worked example presented here applies peat combustion factors (Cf), that have been stratified by fire frequency, according to the authors recent peer-reviewed study 17 .
For biomass burning, the United Nations Framework Convention on Climate Change (UNFCCC) encourage non-Annex I Parties to provide information on anthropogenic emissions by source for CO 2 , CH 4 , N 2 O (Decision 17/CP.8,annex, paragraph14), CO, NO x , and non-methane VOCs (Decision 17/CP.8,annex, paragraph16).Following from this requirement, a further challenge for developing countries is that many of the parameters required for emission estimates are derived from scientific publications that are only available at considerable cost via journal subscriptions or paywall access to the required articles.Also, parameters for some compounds are not available, such as GWP for NOx 18 .For many gases, such as non-methane volatile organic compounds (NMVOC) or NOx-an abbreviation for nitric oxide (NO) and nitrogen dioxide (NO 2 ), specific G i are not easy to identify for staff responsible for the national level reporting, who may not have a chemistry background.Thus, on the one hand, it is a good practice to report emissions as comprehensively as possible, but on the other hand many countries lack the experience and capacity to do so, or elect not to report all gases.
The UNFCCC established the REDD+ scheme as a means for developing countries to receive result-based payments for reducing emissions from deforestation and forest degradation.To achieve REDD+ compliance each country is required to establish a forest reference level (FRL) or forest reference emissions level (FREL) and a process for monitoring, verification and reporting of annual forest carbon stock changes and emissions.The uptake and overall success of the REDD+ scheme requires that each country develops expertise in GHG emissions reporting-a significant step given the complexities of the policy, the underlying science, and the many steps to implementation.For these reasons capacity building is an essential component of climate change mitigation and adaptation initiatives that are promoted through UNFCCC mechanisms.Countries where emissions from biomass burning is one of the major sources of national emissions need to have appropriately skilled and experienced staff, thus training and capacity building is essential.Emissions reporting also relies on relevant peer-reviewed publications to provide a transparent and defensible scientific basis for parameters selected and applied in emissions calculations.
As an example of the complex nature of emissions reporting, the Government of Indonesia (GoI) did not report peat fire emissions for its first FREL submitted in 2016, due to a high level of uncertainty in the parameters used in the estimates, including those as the Appendix 19 .Emissions from peat fires were included in their second FREL submission in 2022 due to improved capacities and due to newly published parameters for peat combustion factors (Cf PEAT ) 17 .This is a significant development for the improvement of global emissions estimates from peatlands.
In this study we provide a simple example of estimating emissions from peat fires using the IPCC methodology in a transparent and accountable way.We stratified our peat swamp forests into three Peatland cover classes: (1) Primary or secondary long-unburnt peat swamp forests, (1st fire); (2) secondary peat swamp forests burnt in 2-3 consecutive fires (2nd-3rd fires), and (3) degraded peatlands, frequently burnt, shrub-dominated (> 3 fires).We combined data from various sources (details below) to separate emissions by fire frequencies.While our intention was to report emissions for all major GHGs as encouraged by the UNFCCC, a review of the literature for G i and GWP revealed that we can estimate emissions for only four major GHGs (CO 2 , CO, CH 4 and N 2 O). ( 4)

Methods
Study rationale and approach.To describe a worked example for the calculation of GHG emissions from major landscape fires in tropical peatlands we used published activity data (i.e.area burnt by peat landcover class) from the 1997 fires in peatlands of Central Kalimantan, Indonesia 7 .The 1997 fires were unprecedented in the history of Indonesia, burning through 729,500 million hectares of peatlands (Fig. 2).We also chose the 1997 fires because of their wide coverage in the scientific literature, including fire emissions estimates equivalent to between 13 and 40% of the mean annual global carbon emissions accounted from fossil fuels 7 .
The methodology applied in this study estimated emissions separately for peat, and for the aboveground biomass, using the emissions equations (Eq.3-7) as described above (a Tier 2 IPCC method).
Datasets.AGB and peat bulk density parameters were derived from three datasets based on field measurements in peat swamp forests of Kalimantan 17,20,21 according to the categories of fire frequency as described above.The bulk density and biomass parameters were bootstrapped for 10,000 iterations and the average and 95% CI values were used in the calculations.
The other parameters such as h, Gi, Cf were derived from the literature (Table 1).
Estimating mass and EF for peat and AGB.As described above in Eq. ( 3), peat mass (M PEAT ) was estimated from peat bulk density (BD) and the depth (h) of peat burnt, while peat EF was estimated from peat mass (M PEAT ) and a combustion factor (Cf PEAT ) according to fire frequency (Eq.4) and (Table 2).AGB was estimated as the sum of biomass components-trees, litter, deadwood and pyrogenic carbon, pending data availability (Table 3).
Developing fire area map.A fire area map for the 1997 fires in Kalimantan was developed using QGIS

Results and discussion
GHG emissions from peat layer and aboveground biomass.Emissions were higher from the 1st fire compared to more frequent fire categories, reflecting the greater depth of peat burnt in first fire (Table 4).Total emissions decreased with increased fire frequency as the depth of peat burnt decreased.Emissions from AGB followed the same pattern as for peat layer with the highest emissions from long unburnt forests (Table 5).

Carbon released from 1997 peat fires-comparison of estimates reported in Page et al. 2002 with the estimates from this study.
As mentioned in the introduction, the total burnt area for Kalimantan from the study of 7 was applied to each of the peatland cover class reported (Table 6), and emission equation parameters derived in this study were applied to estimate total emissions for the peat layer (Table 3) and for AGB (Table 4).Our emissions results for the study area are 0.14 Gt C, half the reported 0.24-0.28Gt C in 7 .The 50% lower emissions estimate can be attributed to the following factors: (1) changes in Cf PEAT : the 7 study used Cf PEAT = 1, assuming that 100% of peat up to 50 cm depth was emitted to the atmosphere over the whole of burnt peatlands, which is likely an overestimation and an oversimplification of the heterogeneity of burnt peatland; (2) depth of peat burnt: 50 cm of peat burnt was used in 7 -later studies 14,23 proposed a shallower depth of peat burnt-33 cm-the value used in our estimates; (3) AGB: the study by 7 used a single AGB carbon density of 250 t C ha −1 (i.e.500 t ha −1 )-applied to primary PSF, logged over PSF and bushlands alike.We note that AGB of 500 t ha −1 reflects highly productive forests and our dataset, which included data from multiple sources with 110 data points had only two forests with AGB reported above 500 t ha −1 .Surely, logged over PSF and bushlands would not have 500 t ha −1 stored in AGB (Fig. 1).
The average emission estimates in our study of 1.92 × 10 −7 Gt C ha −1 (i.e., the total emission of 0.14 Gt C divided by the total area burnt of 726,782.1 ha), demonstrates that the 1997 fires released between 0.47 Gt C, an intermediate estimate (1.92 × 10 −7 Gt C ha −1 × 2,441,000 ha) and an upper estimate 1.31 Gt C (1.92 × 10 −7 Gt C ha −1 × 6,804,668 ha), where 2,441,000 ha and 6,804,668 ha are the peat area burnt over the entire Indonesia extracted from 7 .This is a half of the 0.81-2.57Gt C emissions previously reported.Our estimates match well with the top-down estimates of GFEDv4 of 0.5 Gt C released from the 1997 fires.While GFEDv4 estimates produced a reasonable assessment of Indonesian 1997 peat fires emissions, it may be technically challenging for some countries to use GFEDv4.The example emissions calculations shown in this study are relatively simple, can be completed in an excel spreadsheet, and also do not require an in depth knowledge of rather complex and lengthy IPCC methodologies.
To better enable the development of climate change mitigation actions we believe that building the capacity of countries that aspire to REDD+ payments, through international collaborations with countries more experienced in emission reporting, is a priority for the UNFCCC COP meetings to address.This study can help to guide communities from REDD+ countries, and anyone interested more generally in emissions reporting, to make informed IPCC-method compliant estimates.

Conclusions
In this study we demonstrate the application of IPCC compliant methodology to estimate and report GHG emissions (CO 2 , CO, CH 4 , N 2 O) from tropical peatland fires impacting both aboveground and peat carbon pools.The method is transparent and accounts for a range of fire regimes observed in former peat swamp forest areas over recent decades-to account for either different stages of forest recovery or ongoing chronic fire-disturbance cycles.This nuancing of GHG reporting is especially important in degraded PSF areas because of their potential contribution to either carbon storage (sequestration), or to GHG emissions.This example is provided to assist countries aspiring to participate in the REDD+ scheme to more efficiently navigate the compliance requirements for estimating and reporting emissions from tropical peatlands.
The recently published peat combustion factor data 17 has been used in Indonesia's 2nd FREL submitted to the UNFCCC in early 2022, and it has been assessed by UNFCCC Technical experts.In addition, Indonesia's national GHG inventory reporting is also being updated following the use of the new data in the 2nd FREL report.
We suggest that a more detailed calculation of peat fire emissions according to fire frequency class is achievable in future reporting.It will require detailed activity data to differentiate the areas that have been burnt once, twice, or multiple times which would require resources and high technology/capacity to run the process at the national scale.
We hope that more mechanisms will be developed to increase capacity in developing countries for REDD+ compliance so that they can lead climate change mitigation actions at local to regional and to global levels.Table 6.Comparison of the estimates of carbon released from 1997 peat fires using the emission data from Tables 4 and 5 and as reported in 7 .Total study area and damage in class are derived from Tables 1 and 2 of 7 .

Peatland cover class
Total study area, 10

Figure 1 .
Figure 1.A schematic diagram of the relationship between fire frequencies and peat swamp forests Aboveground Biomass and peat depth.

3. 3 (
Open Source Geospatial Foundation) by overlaying the NOAA-AVHRR Hotspot 1998 layer with the Mawas

Figure 2 .
Figure 2. Area of peat swamp forests burnt in 1997 fires.Map created using the Free and Open Source QGIS.

Table 1 .
Parameters required to estimate emissions from peat and AGB and their sources.

biomass, AGB Peat To estimate CO 2 -eq AGB Cf AGB G i_AGB h BD Cf PEAT G i_PEAT GWP
25rnt area 1997, the Natural Forest Cover 1996, and the peat distribution and administration boundary for Kalimantan.The 1996 Natural Forest Cover map was provided by the Ministry of Environment and Forestry, and the Peat Distribution Map was provided by the Ministry of Agriculture.The Mawas Burnt Area for 1997 was extracted from25.

Table 2 .
Parameters required to estimate M PEAT and EF PEAT and their sources.Values are means (95%CI).

Table 3 .
Parameters required to estimate M ABG and EF AGB and their sources.Values are means (95%CI).