FAOSTAT estimates of greenhouse gas emissions from biomass and peat fires

The Global Fire Emissions Database (GFED3) and the FAOSTAT Emissions database, containing estimates of greenhouse gas (GHG) emissions from biomass burning and peat fires, are compared. The two datasets formed the basis for several analyses in the fifth Assessment Report of the Intergovernmental Panel on Climate Change (IPCC AR5), and thus represent a critical source of information for emissions inventories at national, regional and global level. The two databases differ in their level of computational complexity in estimating emissions. While both use the same burned area information from remote sensing, estimates of available biomass are computed in GFED3 at tier 3 using a complex dynamic vegetation model, while they are computed in FAOSTAT using default, tier 1 parameters from the Intergovernmental Panel on Climate Change (IPCC). Over the analysis period 1997–2011, the two methods were found to produce very similar global GHG emissions estimates for each of the five GFED aggregated biomass fire classes: i) Savanna; ii) Woodland; iii) Forest; iv) Deforestation; v) Peatlands; with total emissions ranging 6–8 Gt CO2eq yr-1. The main differences between the two datasets were found with respect to peat fires, with FAOSTAT showing a lower 1997–1998 peak in emissions compared with GFED3, within an otherwise good agreement for the rest of the study period, when limited to the three tropical countries covered by GFED. Conversely, FAOSTAT global emissions from peat fires, including both boreal and tropical regions, were several times larger than those currently estimated by GFED3. Results show that FAOSTAT activity data and emission estimates for biomass fires offer a robust alternative to the more sophisticated GFED data, representing a valuable resource for national GHG inventory experts, especially in countries where technical and institutional constraints may limit access, generation and maintenance of more complex methodologies and data.

show that FAOSTAT activity data and emission estimates for biomass fires offer a robust alternative to the more sophisticated GFED data, representing a valuable resource for national GHG inventory experts, especially in countries where technical and institutional constraints may limit access, generation and maintenance of more complex methodologies and data.

Introduction
Biomass fires emit greenhouse gases (GHGs) in the atmosphere, where they affect the interannual variability and growth rate of CO 2 , CH 4 and N 2 O and other trace gases . Anthropogenic emissions of non-CO 2 gases from peat and forest fires were 300 million tCO 2 eq yr -1 on average over the period 2001-2010 , with prescribed burning of savanna adding 200 million tCO 2 eq yr -1 (Tubiello et al. 2014). Slash and burn practices associated with deforestation were furthermore responsible for over 1.5 billion tCO 2 eq yr -1 as direct CO 2 gas .
Specific fire emission categories are significant in terms of their overall regional impacts. For instance, in Africa, prescribed burning of savanna represents 25 % of total agricultural emissions (Tubiello et al. 2014). In South-east Asia, fires in drained tropical peatland, often linked to palm oil cultivation, can lead to significant emissions in specific years, with impacts seen at large regional and even global scales (Hayasaka et al. 2014).
Fire emissions linked to human-induced land-use and land use change activities, such as prescribed burning of savanna, deforestation and peatland degradation, need to be reported by countries to the UN Framework Convention on Climate Change (UNFCCC) . Reporting obligations include compilation and submission of national GHG inventories within both National Communications and Biennial Update Reports. The Guidelines for National Greenhouse Gas Inventories, developed by the Intergovernmental Panel on Climate Change (IPCC 2006) allow countries to estimate and report data at various tiers of complexity, from simpler (Tier 1) to more complex (Tiers 2 and 3) approaches.
The simplest Tier 1 approach consists in the multiplication of burned area by default, static parameters describing the biomass available for combustion per hectare, and by GHG emission factors, differentiated by broad geographical, climatic or vegetation characteristics. Tier 2 approaches use the same method, but employ country specific, rather than default parameters. Tier 3 approaches allow for the use of more complex models to derive the necessary parameters, including simulation of biophysical variables at fine temporal and geospatial scales.
Most published estimation methods for quantifying GHG emissions from biomass fires are Tier 3 approaches (e.g., Stroppiana et al. 2010). Experts in developing countries however, including those tasked with implementing national GHG inventories, often face technical and institutional constraints that may limit access, generation and maintenance of data using the complex estimation methods linked to these higher tiers. Understanding differences in and availability of emission data estimated at differing levels of complexity is therefore important in order to plan and implement sustainable national GHG inventory systems.
To this end, we compared data from two available methodologies from the recent IPCC AR5 Report: the Global Fire Emissions Database version 3.1 (GFED3) ( Van der Werf et al. 2010), used in WGI chapter 6 on the carbon cycle (Ciais et al. 2013); and the FAOSTAT Emissions database (Tubiello et al. 2013), used in WGIII chapter 11 on Agriculture, Forestry, and Other Land Uses (AFOLU) .
While both databases rely on remote-sensing information on burned area by land cover class, GFED3 is a Tier 3 approach, employing a dynamic vegetation model to determine fuel biomass consumption values at fine temporal and spatial scales. By contrast, FAOSTAT applies a Tier 1 approach using default vegetation parameters from the 2006 IPCC Guidelines (Tubiello et al. 2014). The aim of this work is to conduct a first robust assessment of differences and similarities between the two databases. While the results of such a comparison are of general scientific interest, given the above-discussed importance of biomass fires emissions at global and regional scales, they also provide specific, valuable information to national GHG inventory compilers, especially in developing countries, towards assessing the costs and the benefits of moving from simple approaches, of easy applicability and sustainability, to more sophisticated modelling tools requiring complex data processing and management.

Materials and methods
The 2006 IPCC Guidelines present a generic equation for the estimation of biomass and peat fires emissions, based on the bottom-up approach first proposed by Seiler and Crutzen (1980): and & G ef,x,v is the emission factor for gas x and for the vegetation type v.
Altogether, M B,v and C f,v are referred to as fuel biomass consumption values. They can be derived either from field data or via models, while emission factors are derived from experimental data. The following sub-sections specify the input data used in Eq. (1) by the two methods considered herein.

Burned area and land cover classes
Both GFED3 and FAOSTAT employ the GFED Burned Area dataset version 3.1 (GFED3-BA)  as input for burned areas in Eq. (1) above. GFED3-BA is largely based on the Moderate Resolution Imaging Spectroradiometer (MODIS) Direct Broadcast Burned Area Monthly L3 Global 500 m product (MCD64A1) (Giglio et al. 2009), with time coverage extended to the pre-MODIS era back to June 1996, by harmonizing MODIS data with those from other sensors, namely the Tropical Rainfall Measuring Mission (TRMM) Visible and Infrared Scanner (VIRS) and the Along-Track Scanning Radiometer (ATSR) family of sensors. The GFED Burned Area Datasets provide the longest consistent global time series of observed burned area data available in the literature. The GFED3-BA covers the period June 1997-March 2012, with monthly maps at 0.5-degree spatial resolution; the more recent GFED4-BA (Giglio et al. 2013) covers the period June 1995-March 2014, with a 0.25-degree spatial resolution and daily/monthly temporal resolution. Each data value represents the total area burned in a given cell over 1 month. An additional data layer contains the share of the total burned area by land cover class, following the University of Maryland (UMD) classification system (Hansen et al. 2000), identified through the MODIS Land Cover Type Yearly L3 Global 500 m product (MCD12Q1) (Friedl et al. 2010). The UMD classification contains the following six land cover classes: i) Savanna; ii) Woody savanna; iii) Open shrublands; iv) Closed shrublands; v) Grassland; vi) Forests. Yearly burned area maps are produced in FAOSTAT by extracting monthly GFED-BA data and then aggregating in yearly composites for each UMD land cover class.
Emissions from peatland fires are also computed in both GFED3 and FAOSTAT, in an attempt to quantify the additional burning of the peat soil component underlying a given burned area and land cover type. The association of a given burned area to peat fire is made independently from the MODIS MCD12Q1 land cover information, as peat soil is not a land cover type. To this end, GFED3 uses a map of ecosystem types (Olson et al. 2001), as a proxy to associate burned areas to peat fires, but limited to Brunei, Indonesia and Malaysia. FAOSTAT applies instead the Harmonized World Soil Database (HWSD) (FAO/IIASA/ISRIC/ISSCAS/JRC 2012), a robust proxy for characterizing peat soils globally (IPCC 2006).

Fuel biomass consumption values
Estimates of fuel biomass consumption values represented the main difference between GFED and FAOSTAT. GFED employs a revised version of the Carnegie-Ames-Stanford-Approach (CASA) dynamic vegetation growth model (Potter et al. 1993) to estimate available biomass amounts as a function of climate and soil parameters, at each pixel and over time for the period of interest. Fuel biomass consumption values and their associated GHG emissions are estimated in monthly intervals and were then aggregated to yearly values. By contrast, FAOSTAT utilizes IPCC default annual fuel biomass consumption values from IPCC (2006), differentiated by different vegetation, soil and climate types (see Table 1), but constant over time (Tubiello et al. 2014). Peatland coefficients were taken from the IPCC Wetlands Supplement (IPCC 2013).
In order to match the IPCC classes with the UMD land cover classes used by GFED3-BA, and thus assign appropriate IPCC fuel biomass consumption values to each pixel, we determined the agro-ecological zones and ecosystem types associated to specific IPCC fuel biomass consumption values, based on IPCC references, and then assigned these values using geo-referenced maps. Specifically, IPCC fuel biomass consumption values for forest subcategories (i.e., boreal, temperate, tropical forest) were assigned to burned forest area pixels by using the FAO Global Ecological Zones map (FAO 2012). For pixels characterized by other vegetation types (e.g., savanna, woody savanna, grassland), tropical and non-tropical IPCC parameter values were assigned using the JRC-IPCC Climate Zones map (JRC 2010).

Emission factors
GFED3 and FAOSTAT use very similar GHG emissions factors for Eq. (1). The GFED3 values, derived from Andreae and Merlet (2001) and Christian et al. (2003), are consistent with the default IPCC 2006 values used by FAOSTAT. For peatlands, for which no 2006 IPCC emission factors exist, FAOSTAT applied those used by GFED3 (Table 1).

GHG emissions
GFED3 emissions are computed by classifying fires into six land cover classes: i) Deforestation; ii) Forest; iii) Woodland; iv) Savanna and Grassland; v) Peatland; and vi) Agriculture. While the underlying, GFED-BA burned area data follow the six land cover UMD classes, as discussed, the GFED emissions data are partitioned using UMD classes and ancillary datasets, such as the MOD44 MODIS continuous vegetation field product and the already mentioned ecosystems map. By contrast, FAOSTAT emissions classes are solely based on the UMD land cover classification, and aggregated in line with IPCC guidelines and UNFCCC reporting requirements. FAOSTAT emissions thus include the following six classes: i) Humid Tropical Forest; ii) Other Forest; iii) All Savanna (subdivided in savanna, woody savanna, closed and open shrublands, and grasslands); iv) Peatland; and v) Agriculture.
For the purpose of this work, it was necessary to re-aggregate FAOSTAT emission classes into those used by GFED. The re-aggregation performed was rather straightforward (Table 2) and allowed for direct comparison between the two emissions datasets. The overall mapping of UMD-based FAOSTAT into the GFED emission classes deforestation and peat fires nonetheless required application of common sense and expert judgement, corroborated by repeated testing. This is because in GFED, fire emissions are associated to deforestation and peat fires based on more than land cover classification. Specifically, in GFED the deforestation fire emission class is determined by using additional MODIS fire products, and the peat fires emission class is determined by using additional information on soil water status, obtained via complex hydrological modelling. In the end, FAOSTAT emission estimates from fire in humid tropical forests were used to compare to the GFED deforestation class, and FAOSTAT emissions estimates from peatlands were used to compare to the GFED peat fires class.
The comparison between GFED and FAOSTAT emissions data was limited to CH 4 , N 2 O and CO 2 gases, as these are those included in national GHG inventories. We did not compare emissions for the agricultural fires class, since FAOSTAT does not use GFED-BA products to estimate them, but rather the simpler IPCC Tier 1 approach, based on total harvested area by crop from national statistics. Comparisons over the remaining five classes were furthermore performed using CO 2 eq units rather than considering the three single gases CO 2 , CH 4 and N 2 O, since trends in each component gas are fully correlated to their equivalent total. This is because the emission factor for each gas is a constant multiplied by the same underlying biomass combustion value. Finally, a downscaling factor of 10 was applied to all GFED3 burned area data, in order to reconcile data from the GFED online database with those published in van der Werf et al. (2010) (L. Giglio, personal communication).
The reader is advised that while we discuss herein results based on the GFED3-BA products in order to maintain coherence with published GFED emissions data and with the data used in the IPCC AR5, the FAOSTAT GHG emission estimates currently distributed by FAO (Tubiello et al. 2015) are already based on the new GFED4-BA dataset (Giglio et al. 2013).

Results and discussion
Yearly GHG emissions from biomass burning produced by GFED3 and FAOSTAT were compared over the period 1997-2011, at both global scale and grid cell level. This paper reports mainly on the results of the comparisons performed at global scale, with some additional information on the grid cell level analysis. The overall GHG estimates were found to be similar between approaches, i.e., not statistically different for all the classes considered (Table 3, Fig. 1).
GFED3 and FAOSTAT estimates indicate global total emissions over the study period at 7.18 ± 1.37 Pg CO 2 eq yr -1 and 6.92 ± 0.77 Pg CO 2 eq yr -1 , respectively, with moderate positive correlation, i.e., with a Pearson's correlation coefficient ρ = 0.53 (p < .05) (Table 3, Fig. 1). The larger variability of GFED versus FAOSTAT data observed at the global level, i.e., coefficient of variation (CV) = 19 % vs. CV = 11 %, respectively, reflects the use of the CASA model by GFED to derive variable biomass combustion values, as opposed to the static default IPCC values used by FAOSTAT. Further investigation confirmed this feature for the other land cover classes, The overall good agreement over the study period 1997-2011 between global GFED and FAOSTAT total emission estimates was also evident maintained for single fire emission classes (Table 3, Fig. 2). For deforestation, global GHG estimates with GFED and FAOSTAT were 1.48 ± 0.62 Pg CO 2 eq yr -1 and 1.66 ± 0.37 Pg CO 2 eq yr -1 , respectively, with good positive correlation, i.e. ρ = 0.71 (p < .05). For forest, estimates were 1.02 ± 0.35 Pg CO 2 eq yr -1 and 1.01 ± 0.17 Pg CO 2 eq yr -1 , respectively, also with good positive correlation, i.e. ρ = 0.81 (p < .05). For savanna, estimates were 3.15 ± 0.34 Pg CO 2 eq yr -1 and 2.86 ± 0.51 Pg CO 2 eq yr -1 , respectively, with mild positive correlation, i.e. ρ = 0.46 (p < .1). For woodland, estimates were 1.16 ± 0.20 Pg CO 2 eq yr -1 and 1.20 ± 0.13 Pg CO 2 eq yr -1 , respectively, with good positive correlation, i.e. ρ = 0.82 (p < .05). Finally for peat fires, only estimates for Brunei-Darussalam, Malaysia, and Indonesia were compared, since GFED3 only provides data for these countries. In this geographical region, estimates were 0.37 ± 0.69 Pg CO 2 eq yr -1 and 0.21 ± 0.30 Pg CO 2 eq yr -1 , respectively, with correlation ρ = 0.99 (p < .2). FAOSTAT peat fires emissions were in general in good agreement with GFED3, although peaks in El Niño years 1997, 2005 and 2008 were higher for GFED3. This implies that the hydrological analysis method employed by GFED3 associates larger peat biomass available for burning in years with pronounced drought, such as those associated with the El Niño phenomena-while the IPCC methods employed by FAOSTAT result in constant values irrespective of soil drought status.
Importantly, it should be noted that FAOSTAT emission estimates for peat fires apply globally, while emissions estimates by GFED only cover three tropical countries: Brunei-Darussalam, Malaysia, and Indonesia. Global estimates of peat fires emissions in FAOSTAT, including both tropical and boreal zones, were 1.00 ± 0.37 Pg CO 2 eq yr -1 over the period This results indicate that global emissions from peat fires, once both boreal and tropical fires are included, could be several times larger than those currently estimated by GFED (Fig. 2). While the above mean emission values by class were found to be not statistically different, differences in CVs between GFED and FAOSTAT were apparent for several emission classes (Table 3). This is related to the inherent higher interannual variability of biomass used in GFED compared to the static values used in FAOSTAT. This mechanism may explain differences in CV found for forest (34 and 16 %, respectively), woodland (17 and 11 %, respectively) and deforestation (42 and 22 %, respectively), although we did not have access to the underlying biomass data to verify quantitatively this conclusion. For peat fire emission estimates CVs were 189 and 132 %, respectively.
The agreement between GFED and FAOSTAT at regional scale was assessed producing geo-spatial maps with grid cell-level analyses of GFED and FAOSTAT data over the period 1997-2011, with respect to single emission classes (Fig. 3). The correlation at pixel level was analysed using the Pearson correlation coefficient, while the difference between the two datasets was assessed using the normalized mean absolute error (NMAE), as suggested by Willmott and Matsuura (2005). The results suggest that the significant consistency between the two methodologies that was found at global level is also maintained at regional scale, given the same burned area datasets used as input by the two methodologies, and despite the use of a complex dynamic biomass predictor in GFED versus the simpler, static FAOSTAT approach described earlier. Relatively lower correlations and higher difference between the data series at pixel level were found in mixed transition areas among different ecosystems, probably due to the different procedures adopted to map the original UMD land cover classes used by the GFED-BA dataset into the classes GFED used to group emissions ( Table 2).
Results of the comparisons undertaken herein confirm the original analysis of GFED on the global distribution of GHG emissions from biomass fires (Fig. 4), also summarized in several chapters of IPCC AR5. Global total fire emissions were in the range of 6-8 Gt CO 2 eq yr -1 over the period 2001-2010. Of this, savanna fires contributed about one-third of total emissions, or about 2 Gt CO 2 eq yr -1 . Deforestation, woodland and forest fires emissions contributed the remaining two-thirds. Only about one-third of these emissions constitute anthropogenic =0.64 forcing of the atmosphere and are to be included in national GHG inventories. These include about 1.7 Gt CO 2 eq yr -1 (of which 1.5 Gt CO 2 eq yr -1 as CO2 gas, 0.06 Gt CO 2 eq yr -1 as N 2 O, and 0.13 Gt CO 2 eq yr -1 as CH 4 ) from deforestation fires, and 1 GtCO 2 eq yr -1 (of which 0.77 Gt CO 2 eq yr -1 as CO 2 , 0.03 Gt CO 2 eq yr -1 as N 2 O, and 0.2 Gt CO 2 eq yr -1 as CH 4 ) from peat fires, while non-CO 2 emissions from savanna and other forest fires amount in total to some 0.3 Gt CO 2 eq yr -1 (of which 0.15 Gt CO 2 eq yr -1 as N 2 O, and 0.15 Gt CO 2 eq yr -1 as CH 4 ) Tubiello et al. 2014).

Conclusions
This work investigates differences between GHG emissions from biomass fires, comparing the FAOSTAT estimates, obtained using a default Tier 1 approach according to the IPCC  estimates. The data estimated with the two methods were shown to have close correspondence with respect to all the land cover classes considered. This comparison exercise will be repeated with FAOSTAT and GFED4 emission data. Furthermore, in order to improve access to FAOSTAT historical and updated data, we will examine the possibility to employ directly the MCD45A1 burned areas dataset from the NASA MODIS programme instead of GFED processed data.
In conclusion, this study shows that the FAOSTAT Tier 1 approach, using the default widelyapplicable IPCC methodology, produces in general estimates of emissions from biomass burning in line with the results of the widely cited GFED3 database, which are obtained through a complex vegetation growth modelling approach (Tier 3) much more demanding in terms of input data and know-how. This is very important information for inventory compilers from an operational point of view. In many contexts, technical constraints or lack of data limit application and maintenance of higher Tiers methodologies. This is especially the case of developing countries, which starting with 2014 are required to regularly produce estimations of their GHG emissions and submit Biennial Update Report (BUR) to the UNFCCC. Our results show that in these contexts the Tier 1 approach can be considered a reliable option for use in national GHG inventories. FAO makes emissions data available free of charge at www.faostat.fao.org.