Articles | Volume 18, issue 4
Biogeosciences, 18, 1375–1393, 2021
Biogeosciences, 18, 1375–1393, 2021

Research article 23 Feb 2021

Research article | 23 Feb 2021

Intraseasonal variability of greenhouse gas emission factors from biomass burning in the Brazilian Cerrado

Intraseasonal variability of greenhouse gas emission factors from biomass burning in the Brazilian Cerrado
Roland Vernooij1, Marcos Giongo2, Marco Assis Borges3, Máximo Menezes Costa3, Ana Carolina Sena Barradas3, and Guido R. van der Werf1 Roland Vernooij et al.
  • 1Department of Earth Sciences, Faculty of Science, Vrije Universiteit Amsterdam, Amsterdam, the Netherlands
  • 2Center for Environmental Monitoring and Fire Management (CEMAF), Federal University of Tocantins, Gurupi, Brazil
  • 3Chico Mendes Institute for Conservation of Biodiversity (ICMBio), Rio da Conceição, Brazil

Correspondence: Roland Vernooij (


Landscape fires, often referred to as biomass burning (BB), emit substantial amounts of (greenhouse) gases and aerosols into the atmosphere each year. Frequently burning savannas, mostly in Africa, Australia, and South America are responsible for over 60 % of total BB carbon emissions. Compared to many other sources of emissions, fires have a strong seasonality. Previous research has identified the mitigation potential of prescribed fires in savanna ecosystems; by burning cured fuels early in the dry season when landscape conditions still provide moist buffers against fire spread, fires are in general smaller, patchier, and less intense. While it is widely accepted that burned area (BA) and the total carbon consumed are lower when fires are ignited early in the dry season, little is known about the intraseasonal variability of emission factors (EFs). This is important because potentially, higher EFs in the early dry season (EDS) could offset some of the carbon benefits of EDS burning. Also, a better understanding of EF intraseasonal variability may improve large-scale BB assessments, which to date rely on temporally static EFs. We used a sampling system mounted on an unmanned aerial vehicle (UAV) to sample BB smoke in the Estação Ecológica Serra Geral do Tocantins in the Brazilian states of Tocantins and Bahia. The protected area contains all major Cerrado vegetation types found in Brazil, and EDS burning has been implemented since 2014. Over 800 smoke samples were collected and analysed during the EDS of 2018 and late dry season (LDS) of 2017 and 2018. The samples were analysed using cavity ring-down spectroscopy, and the carbon balance method was used to estimate CO2, CO, CH4, and N2O EFs. Observed EF averages and standard deviations were 1651 (±50) g kg−1 for CO2, 57.9 (±28.2) g kg−1 for CO, 0.97 (±0.82) g kg−1 for CH4, and 0.096 (±0.174) g kg−1 for N2O. Averaged over all measured fire prone Cerrado types, the modified combustion efficiency (MCE) was slightly higher in the LDS (0.961 versus 0.956), and the CO and CH4 were 10 % and 2.3 % lower in the LDS compared to the EDS. However, these differences were not statistically significant using a two-tailed t test with unequal variance at a 90 % significance level. The seasonal effect was larger in more wood-dominated vegetation types. N2O EFs showed a more complex seasonal dependency, with opposite intraseasonal trends for savannas that were dominated by grasses versus those with abundant shrubs. We found that the N2O EF for the open Cerrado was less than half the EF suggested by literature compilations for savannas. This may indicate a substantial overestimation of the contribution of fires in the N2O budget. Overall, our data imply that in this region, seasonal variability in greenhouse gas emission factors may offset only a small fraction of the carbon mitigation gains in fire abatement programmes.

1 Introduction

Landscape fires emit large amounts of greenhouse gases and aerosols, which significantly impact atmospheric chemistry and biogeochemical cycles on local to global scales (Andreae and Merlet, 2001; Reid et al., 2005; van der Werf et al., 2017). The primary greenhouse gases emitted from biomass burning are carbon dioxide (CO2), methane (CH4), and nitrous oxide (N2O). Over the period 1997–2016, average annual emissions from landscape fires were 7.3 Pg CO2, 16 Tg CH4, and 0.9 Tg N2O according to the Global Fire Emissions Database (GFED4s, Van der Werf et al., 2017). Tropical savannas accounted for the majority of these global landscape fire emissions with 4.9 Pg CO2, 6 Tg CH4, and 0.6 Tg N2O. South American savannas on average accounted for about 10 % of the global fire-related carbon emissions in savannas, corresponding to 6.5 % of the total fire-related carbon emissions over this period. In general, biomass burning CO2 emissions are compensated for by regrowth of vegetation after the fire (Beringer et al., 2007; Landry and Matthews, 2016). Therefore, fires only impact long-term atmospheric CO2 concentrations when regrowth does not take place (e.g. following deforestation and tropical peatland fires) or if fire regimes become more severe and result in regional carbon sources. Although the Intergovernmental Panel on Climate Change (IPCC) recommends national accounting of indirect GHGs (greenhouse gases) like CO (carbon monoxide) and NMHCs (non-methane hydrocarbons) from biomass burning (Goodwin et al., 2019), carbon monoxide is generally not considered in discussions on emission abatement schemes (Australian Government – Department of the Environment and Energy, 2018; Cook et al., 2015; Lipsett-Moore et al., 2018). Like CH4, CO reacts with atmospheric OH radicals and is eventually oxidized to CO2 (Crutzen and Zimmermann, 1991; Daniel and Solomon, 1998). The depletion of OH radicals by enhanced CO concentrations leads to an increase in the atmospheric lifetime of CH4 and the formation of ozone (O3) (Crutzen and Zimmermann, 1991; Fry et al., 2012; Sudo and Akimoto, 2007). Albeit poorly understood, CO and CH4 indirectly affect clouds by altering the abundance of oxidants which convert SO2 into sulfate (Penner et al., 2006). Therefore, CO can be viewed as an indirect greenhouse gas, more potent than CO2 (Daniel and Solomon, 1998; Myhre et al., 2013).

Emission factors (EFs) are used to quantify the conversion of the total amount of carbon and other elements in the consumed fuel to emissions of various trace gases and aerosols. They are often reported in grams per kilogram of dry biomass consumed. Biomass burning EFs are derived from laboratory, ground-based, and aircraft measurements and have been reported for a large number of chemical species and vegetation types (Akagi et al., 2011; Andreae, 2019; Andreae and Merlet, 2001). The modified combustion efficiency (MCE), defined as the amount of carbon emitted as CO2 divided by the amount of emitted carbon in CO2 and CO combined, is often used as an indication of the relative contribution of flaming and smouldering combustion (Akagi et al., 2011). The MCE ranges from about 0.65 in smouldering peat fires to values close to one for highly efficiently oxidizing grass fires. The negative correlation of the MCE with EFs for incomplete combustion products such as CH4, non-methane hydrocarbons (NMHCs), and carbonaceous particulate matter (CPM) (Hoffa et al., 1999; Urbanski, 2013) makes it a useful metric for emission estimations.

A substantial amount of research has been conducted to understand what environmental factors affect the EFs for greenhouse gases (e.g. Chen et al., 2010; Korontzi et al., 2003a; Urbanski, 2014). While the drivers of variability in CO and CH4 EFs have received considerable attention, relatively little is known about the BB contribution to the N2O budget. N2O is formed through the oxidation of HCN and NH3, in which the reaction of HCN through NCO is the dominant pathway. The N2O EF is strongly dependent on the C : N ratio in the fuel (Lobert and Warnatz, 1993) as well as the temperature and partial pressure of oxygen during combustion (Kilpinen and Hupa, 1991; Winter et al., 1999a). CH4 is formed during incomplete oxidation of biomass, with higher EF when fuels are relatively moist (Chen et al., 2010) or when fuels are densely packed (Bertschi et al., 2003; Urbanski, 2013). While some of the drivers of variability in these EFs are qualitatively known, large-scale studies have so far relied on biome-average estimates due to the lack of quantitative information, thus ignoring spatio-temporal variability within biomes (Van Leeuwen and Van der Werf, 2011).

The Cerrado in South America consists of a mosaic of grasslands, shrublands, and forests. The biome covers roughly 24 % of Brazil, as well as smaller parts of Paraguay and Bolivia (Klink and Machado, 2005). Vegetation dynamics and distribution in wild Cerrado areas are primarily determined by water availability, soil type, and fire history (Pivello, 2011). The Cerrado can be categorized based on the abundance of woody species, ranging from campo limpo (open grassland), campo sujo (grassland with sparse presence of shrubs), campo rupestre (rock field), parque cerrado (grass- or shrub-dominated with scattered trees), and cerrado típico (tree-dominated with scattered shrubs and a grass understory) to cerradão (forest savanna). Forested landscapes, such as gallery forests which tend to directly line the river and are found in riparian zones within the Cerrado, are particularly fire-sensitive (Ribeiro and Walter, 1998). Humid grasslands are also found here, consisting of gleysols that remain flooded in the rain season and are typically covered with grass and sparse palm trees. Fires have a dominant role in limiting the proportion of trees in the Cerrado, and fire frequencies generally range from 3 to 8 years (Fidelis et al., 2018). Natural fires in the Cerrado are caused by lightning and mostly occur at the beginning and end of the dry season. Anthropogenic fires, lit for example for cattle-ranging pasture improvement, typically occur around the middle of the dry season in July–August (Pivello, 2011; Ramos-Neto and Pivello, 2000; Schmidt et al., 2018). Fire intensity is a key landscaping factor that can also feed back on the long-term vegetation state; high-intensity fires limit tree cover and promote the formation of open grassland, which in turn promotes higher fire frequency (Miranda et al., 2009; Oliveras et al., 2013; Staver et al., 2011a).

Under seasonally dry conditions, Cerrado grass species dry out and senesce, leading to standing dead fuel accumulation (Fidelis and Fernanda, 2013). Local practices have relied on prescribed burning in the past, and scientific research showed that fire has a key role in maintaining the Cerrado's high biodiversity. However, until the first integrated fire management approach for some protected areas was launched in 2014, a “zero-fire” policy had been maintained in the Brazilian Cerrado for decades (Durigan et al., 2016; Schmidt et al., 2018). While fire suppression strategies can be effective as a tool to enhance carbon sequestration and total carbon stocks (Murphy et al., 2010; Staver et al., 2011a), keeping fire out of the Cerrado altogether potentially leads to a sharp decline in biodiversity through the loss of light-demanding savanna species (Abreu et al., 2017; Durigan et al., 2016). In larger continuous landscapes, fire suppression strategies have led to a shift towards more high-intensity late-dry-season (LDS) fires which are more difficult to suppress. Frequent, high-intensity fires can cause long-term losses of soil nitrogen and phosphorous (Kauffman et al., 1994), which in turn decreases the total amount of carbon that is sequestered by net primary productivity. This may in time alter the carbon sink capacity of frequently burning savanna grasslands (Pellegrini et al., 2018). To combat the rise of intense LDS fires, it is important to look for alternative fire management strategies. Somewhat ironically, fire exclusion experiments have thus shown that well-managed fire intervals and intensities are vital for sustaining biodiversity in fire-prone savanna systems (Abreu et al., 2017; Durigan et al., 2016; Scott et al., 2012).

Given that fire exclusion and thus a fire-free Cerrado is hardly possible or beneficial, previous studies have suggested the potential for prescribed burning in the early dry season (EDS) as an alternative to devastating LDS fires (Fidelis et al., 2018). Cerrado fuels dry at different rates under different landscape conditions (e.g. slopes versus valley bottoms) as the EDS progresses, resulting generally in smaller, more patchy, less intense fires (Rissi et al., 2017). Networks of strategically placed EDS fires can be used to reinforce natural (e.g. riparian corridors) and built (e.g. road) barriers, thereby reducing the risk of more intense, extensive LDS wildfires. For this reason, prescribed EDS burning is suggested as a climate mitigation strategy in the savanna (Anderson et al., 2015; Lipsett-Moore et al., 2018; Penman et al., 2011; Russell-Smith et al., 2013).

Africa and South America combined collectively account for about 65 %–77 % of the fire-prone savannas, and carbon schemes may provide incentivized alternatives for delivering less destructive savanna fire regimes as shown in Australia (Lipsett-Moore et al., 2018; Maraseni et al., 2016; Russell-Smith et al., 2013). Wildfire emissions are the product of fire extent, fuel load, combustion completeness, and EFs for the emitted species (Seiler and Crutzen, 1980). EDS fires have been demonstrated to be smaller and more patchy (Oliveira et al., 2015; Price et al., 2012) and result in lower combustion completeness (Yates et al., 2015). Total fuel consumption is therefore lower. However, through more incomplete combustion under more humid fuel conditions, higher CH4 EFs offset some of the emission gains from reduced fuel consumption (Hoffa et al., 1999; Ito and Penner, 2004; Korontzi, 2005; Van Leeuwen and Van der Werf, 2011; Yokelson et al., 2011). Understanding and quantifying the intraseasonal variability in EFs is therefore essential to assess the implications of natural and human-induced fire regime shifts.

In this study we have used a novel UAV-based (unmanned aerial vehicle) approach to sample smoke from landscape fires during three field campaigns, covering different parts of the dry season and various fire-prone Cerrado vegetation types. Our main objective was to assess the spatio-temporal variability in EFs for the main greenhouse gases associated with BB. With this knowledge we are in a better position to understand the carbon mitigation potential of savanna fire management, and these findings may improve the representation of EFs in large-scale fire databases.

2 Methodology

2.1 Study area

The Estação Ecológica Serra Geral do Tocantins (hereafter referred to as EESGT) is a protected area located in the Brazilian states of Tocantins and Bahia (Fig. 1a). Covering  700 000 ha, it is one of Brazil's largest ecological stations, a type of strictly protected area established to preserve representative samples of the different biomes in Brazil. EESGT used to be one of the most frequently burning protected areas in the Cerrado. On average, about 30 % of the protected area burned each year (Fidelis et al., 2018). Since 2014, prescribed EDS burning has been used within EESGT as a tool to reduce the negative impacts from uncontrolled, high-intensity LDS fires. The strategy focusses on creating a mosaic of smaller areas with different fire histories and hence varied fuel loads (Barradas et al., 2020). Over the years 2015–2018, the strategy resulted in an increase in the number of fires but a decrease in average fire size and total burned area (Barradas, 2017; Fidelis et al., 2018; Schmidt et al., 2018). This is in contrast to other protected areas in the Cerrado, which have implemented limited prescribed burning and experienced an increase in burned area over the same period. For example, 78 % of nearby Chapada dos Veadeiros and 85 % of Reserva Natural Serra do Tombador were burned in 2017, inciting calls for wider implementation of EDS fire management (Fidelis et al., 2018).

Figure 1(a) Location of the Estação Ecológica Serra Geral do Tocantins in the Cerrado biome in the Brazilian state of Tocantins. (b) Vegetation types in the Estação Ecológica Serra Geral do Tocantins (Franke et al., 2018; Orozco-Filho, 2017) with the locations of the measurements.

EESGT has a semi-arid to tropical climate with two well-defined seasons: hot and dry in May to September and wet and cooler in the rainy season between October and April (Fig. 2a). With an average annual rainfall of around 1400–1500 mm, EESGT is somewhat wetter than the Cerrado average of 1300 mm (Seplan, 2003). The area is dominated by nutrient-poor, deep arenaceous quartz soils and has a high floristic biodiversity. All the major Cerrado vegetation types are represented in the ecological station, but the area is dominated by open grasslands (campo limpo and campo sujo) and vegetation of open savanna (cerrado ralo and cerrado típicosensu stricto) (Fig. 1b).

Figure 2Seasonality and interannual variability of (a) daily burned area (BA; MCD64A1 C6; Giglio et al., 2018) as well as monthly rainfall and soil moisture, averaged over the 2013–2018 period. The prescribed-burning season and repression season are hatched. (b) Early-dry-season (EDS, before 1 July) and late-dry-season (LDS, after 1 July) annual burned area and active fire detection (AFD, MOD14A1v6/MYD14A1v6; Giglio and Justice, 2015) over the 2007–2018 period.


The fire season in EESGT roughly lasts from May until October and peaks around September (Fig. 2a). In the EDS, managers apply fires during a “safe-burning window” which depends on the vegetation type, fuel conditions, and weather. Typically, prescribed EDS fires are applied in the afternoon and extinguish after sunset as the temperature and wind speed drop and relative humidity (RH) increases. Managers actively suppress intense LDS fires to protect vulnerable vegetation and surrounding communities (Barradas et al., 2018).

2.2 Measurement campaigns

We carried out one EDS and two LDS measurement campaigns. During the 2017 LDS (23 September–11 October), fires were ignited between 09:30 and 18:00, and air temperature ranged from 25 to 38 C. No relative humidity measurements were taken during this campaign. In the EDS of 2018 (16–30 June), fires were ignited between 12:00 and 18:00. The air temperature during this period, as measured prior to the fire at an elevation of 15 m, was between 31 and 36 C with an average relative humidity of 18 %. During the 2018 LDS (23 September–11 October) samples were collected from prescribed and two non-prescribed fires, in which we sampled smoke from fires in gallery forest and humid grasslands during a LDS fire repelling effort. Although the LDS campaign in 2018 occurred after the first rainfall of the season, which came early in 2018, fire intensities appeared to be much higher than during the EDS campaign (Fig. 3). The temperature ranged from 37 to 42 C with an average relative humidity of 13 %. More information about the number of measurements taken and vegetation type coverage during each campaign is listed in Table 1.

Figure 3Typical post-fire images showing the much smaller impact of EDS fires, in this case in June (a), compared to LDS fires in September (b).


Table 1Number of samples and ancillary information about the field campaigns.

Download Print Version | Download XLSX

2.3 Sampling strategy

We filled single-polypropylene fitted Tedlar bags (SKC, type 232-01) with fresh smoke using a UAV-based (DJI, Matrice 100) sampling system. Most of the samples were taken 15–20 m above the fire, with the height increasing with the intensity of the fire. Our sampling system consisted of a container mounted on top of the UAV which contained four Tedlar sample bags. We filled 1 L bags with ±0.7 L of smoke, which took 35 s for each bag, creating a 35 s averaged mixture of trace gases in the bag. The sample inlet was located on the top of the UAV and fitted with a 60 µm sintered porous metal filter. During the sampling period, the system logged time, GPS coordinates, pressure, temperature, and relative humidity measured at the UAV.

Most sampled fires were ignited by the EESGT rangers using a drip torch to start a fire line of at least 30 m. We started sampling when the fire front had advanced 10–20 m. We positioned the UAV to capture a mixture of the fast-ascending flaming combustion products and the smouldering products that were generated upwind from the flaming fire front. While the majority of the fires sampled were prescribed burns, we also sampled several non-prescribed LDS fires. These fires were most likely escaped pasture fires or poaching fires, given that lightning did not occur during our LDS campaign. We sampled both EDS and LDS fires over various vegetation types with the time since the last fire ranging between 2 and 5 years (Table 1).

2.4 Smoke analysis

We used cavity ring-down spectroscopy to analyse CO2, CO, CH4, and N2O concentrations from the sample bags. After sampling, the Tedlar bags were kept in a dark environment and analysed within 12 h. This was done in order to minimize the oxidation of CO by OH radicals inside the bags. According to Meyer et al. (2012) and our own tests, CO2, CO, N2O, and CH4 concentrations are stable in the Tedlar bags for extended periods under these conditions. The samples were measured for 20 s at a flow rate of 1.3 L min−1 using a CO2 and CH4 analyser (Los Gatos Research, Microportable Gas Analyzer) followed by 20 s at a flow rate of 0.25 L min−1 using a CO and N2O analyser (Aeris Technologies, Pico series); see Table 2. Measurement of the trace gas concentration in the bags was based on the 10 s average concentration following a 10 s initial flushing period. Before each fire, we filled four “background” samples at 15 m altitude. The average concentration of these background samples was subtracted from those in the plumes to get the excess mixing ratio (EMR) in the sample bags. Variability between the background samples during a single day was smaller than 5 %. Both analysers were calibrated before and after each campaign using certified standard calibration gas (Table 2). No significant calibration drift was observed during the campaigns.

Table 2Description of analysis equipment used.

Download Print Version | Download XLSX

2.5 EF calculation

We converted the EMR (sample minus background concentration) in the bags to EFs for CO2, CO, and CH4 in grams of emitted species per kilogram of dry matter burned, following the carbon mass balance method (Urbanski, 2013; Yokelson et al., 1999):

(1) EF i = F c × MM i AM c × C i C total × 1000 g kg - 1 ,

where EFi is the emission factor of species i and Fc is the carbon content of the fuel by weight fraction. In this study, we used 48 % for grassland–savanna and humid grasslands and 50 % for gallery forest, based on carbon content measurements from different Cerrado vegetation types by Susott et al. (1996). MMi is the molecular mass of species i divided by the atomic mass of carbon (AMc). Ci is the number of moles of carbon emitted in species i; Ctotal is the total number of moles of emitted carbon. Because we did not measure non-methane hydrocarbons (NMHCs) and carbonaceous particulate matter (CPM), these fractions were estimated based on ratios from savanna-burning literature. The total mass of emitted CPM was estimated to be 7 % of the emitted mass of CO (Andreae, 2019), with carbon accounting for 70 % of the CPM mass (Reid et al., 2005). The total amount of carbon in NMHCs was estimated to be 3.5 times the ER(CH4/CO2) (emission ratio) based on common ratios for savanna fires (Andreae, 2019). We did not consider residual ash in our calculations, which can represent significant amounts of carbon (Jones et al., 2019). Although this is common practice in EF calculations, leaving out ash may lead to overestimation of carbon emissions (Surawski et al., 2016). To calculate the EFs for N2O, we used Eq. (2) described by Andreae and Merlet (2001). This method uses the emission ratio (ER(i/y)) of the species i to a relatively inert, co-emitted carbonaceous species y.

(2) EF i = ER ( i / y ) × MM i MM y × EF y

We used CO2 as the co-emitted reference gas following earlier work (Hao et al., 1991; Hurst et al., 1994a; Surawski et al., 2015). Although CO is also sometimes used for this purpose due to its low background variability (Meyer et al., 2012), based on previous continuous emission measurements, we found N2O to be more closely correlated with CO2. Calculation using CO as a co-emitted reference gas for the N2O EF on average would lead to N2O EFs that are 2.4 % higher. EFs were calculated for each bag separately, and we partitioned the bags into different seasonal, vegetation type, and fire history classes (Sect. 2.6). To get the weighted-average EF for these classes, we calculated EFs over the cumulative EMR of the respective trace gas species from all samples in the class. Samples with low overall trace gas concentrations thus have low impact on the weighted-average EF.

To assess the intraseasonal effect of the emissions on radiative forcing (RF), we calculated the CO2 equivalent (eq.) EF based on the EFs weighted by the 100-year global warming potentials (GWP100) of the emitted species. CH4 and N2O have a GWP100 of 34 and 298, respectively, when climate–carbon feedback mechanisms are included (Myhre et al., 2013). CO is usually not included, but through its removal of hydroxide it leads to a longer lifetime of CH4, which is a precursor for O3, and eventually oxidizes to CO2 resulting in 1.57 g CO2 per gram oxidized CO (Goodwin et al., 2019). Therefore we have also taken CO into account. Estimates of the indirect CO GWP100 vary from 1.8 (Fry et al., 2012) to 5.4 when taking into account primary and secondary aerosol effects on clouds (Shindell et al., 2009). We used a GWP100 for CO of 2.2, which is on the conservative side of these estimates and does not include the effect of the oxidized CO2, since this is assumed to be compensated for by regrowth.

2.6 Spatial analysis and upscaling

All samples were geolocated using the coordinates of the UAV. This location was used to tag the samples with vegetation type and the number of years since the previous fire derived from satellite data. Most of the plumes were sampled close to the fire, but we manually checked this information with satellite BA data to avoid mismatches due to plume advection. To calculate the fire history of the burned vegetation, we used BA data based on the 30 m Landsat Thematic Mapper (TM), Enhanced Thematic Mapper (ETM), and Operational Land Imager (OLI) from the Instituto Nacional de Pesquisas Especiais (INPE) (Melchiori et al., 2015). The dataset uses consecutive Landsat scenes to detect changes in the Normalized Difference Vegetation Index (dNDVI) and Normalized Burn Ratio (dNBR) for BA classification. The BA classification is manually validated in the field, and thresholds in the algorithm were optimized for EESGT as described by Barradas et al. (2018). The number of years since the last fire was determined based on the location of the sample and the Landsat 30 m burn scars of the last years.

For the vegetation classification described in Table 1, we used maps created by the University of Brasilia (Fig. 1b), which were derived from 5 m RapidEye multispectral imagery (Orozco-Filho, 2017). The classification is based on spectral characterization of the different vegetation types and distinguishes the following Cerrado classes sampled by this study: campo limposujo (open grassland; 0 %–5 % tree cover), cerrado sensu stricto ralo (open Cerrado; 5 %–20 %), cerrado sensu stricto típico (typical Cerrado; 20 %–50 %), cerrado denso (dense Cerrado; >50%), gallery forest (continuous canopy), and riparian zones (sparse palm trees in wetlands). The classification matched well with our field observations during the campaigns, but we did not validate the map formally. It should be noted that the fractional tree cover (FTC) classification in the RapidEye map generally leads to higher FTC values compared to the MODIS-based (Moderate Resolution Imaging Spectroradiometer) vegetation continuous fields dataset (MCD44Bv6, DiMiceli et al., 2015) or the Landsat-based rescaling of the MCD44Bv6 dataset (Sexton et al., 2013). Hence, care should be taken with spatial extrapolation of these vegetation classes using different FTC products.

The weighted-average emission factor (EF) for vegetation types of combined Cerrado in EESGT was calculated through Eq. (3), in which n is the number of vegetation types; BAi is the burned area based on the aforementioned EESGT-optimized INPE BA (Barradas et al., 2018) for the vegetation class (i) over the years 2013 to 2018; and BAtot is the total burned area over the same period (Fig. 4).

(3) EF = i = 1 n EF i × BA i BA tot

Since we lack data on the fuel load and combustion completeness, we weighed the EFs by the percentage of BA in the different classes (Fig. 4). Further, given that we do not have measurements of EFs of dense Cerrado, BA of the dense Cerrado was accounted for as typical Cerrado. As the BA composition of EDS fires primarily depends on management considerations, both seasons were weighed by the total averaged BA composition.

Figure 4Partitioning of the burned area over the EDS (before 1 July) (left columns) and LDS (after 1 July) (right columns) for the various vegetation types.


3 Results

The weighted-average EFs for the different vegetation types, as well as the EF for vegetation of combined Cerrado, are listed in Table 3. Since the introduction of prescribed LDS burning in EESGT in 2014, the proportion of BA before 1 July has gradually increased (Fig. 2b). This has been the case for all dominant fire-prone vegetation types found in EESGT. As the samples were unevenly distributed over the different vegetation types (Table 1), we had to account for the sample bias in vegetation type to compare EDS and LDS EFs, which is our main objective. To obtain a seasonal weighted-average emission factor (EF) for Cerrado vegetation, we therefore weighted the different Cerrado vegetation class EFs by their contribution to the fires in EESGT. Over the 2013–2018 timeframe, the distribution of BA over the different fire-prone ecosystems (vegetation types) most common in EESGT was approximately 23 % in open grassland, 42 % in open Cerrado, 28 % in typical Cerrado, and 7 % in dense Cerrado.

Table 3Weighted-mean EFs (g kg−1) for various vegetation types of EDS and LDS fires. The standard error of the mean (SEM) is given in brackets.

a Using a two-tailed t test with unequal variance, the difference is statistically significant (p<0.1). b N2O weighted-average EFs and SEMs are given for samples with (>15 mol) additional carbon.

Download Print Version | Download XLSX

3.1 Intraseasonal variability

Although the variability within individual fires (we collected several samples from each fire), vegetation types, and campaigns was high, the difference between the EDS- and LDS-averaged CO and CH4 EFs was limited (Fig. 5). The MCE increased slightly from a weighted average of 0.957 in the EDS to 0.963 in the LDS. When considering individual vegetation types, more efficient combustion in the LDS campaigns is apparent. For example, the difference between the LDS and EDS when averaged over all vegetation types (−15 % for CO and −13 % for CH4) is more pronounced when focusing on more shrub-dominated areas (open Cerrado). For example, CO and CH4 EFs were 18 % and 21 % lower in the LDS for vegetation of typical Cerrado (Table 3). As a result of the large spread in EFs and a limited number of samples in some vegetation types, only the slight differences in open grasslands and the 14 % and 34 % increases in N2O EF for open Cerrado and typical Cerrado, respectively, were statistically significant using a two-tailed t test with unequal variance at a 90 % significance level.

Figure 5EFs (g kg−1) in the EDS and LDS as well as the EFs from savanna measurements used in the Andreae (2019) EF compilation. The green diamond represents the arithmetic mean, and the red cross represents the EMR-weighted-mean value. The purple dot represents the value that is used in GFED for savanna fires.


Campaign-averaged N2O EFs were 0.105 g kg−1 in the EDS and 0.123 g kg−1 in the LDS. However, internal variability within the campaigns was high with standard deviations of 0.183 g kg−1 in the EDS and 0.263 g kg−1 in the LDS. In Table 3, N2O EFs are reported for samples with enhanced carbon concentrations of over 15 mol (as explained in Sect. 4.4), in order to minimize propagation of measurement error in the standard error of the mean. Though not significantly altering the weighted average, this improved the significance of the found relationships. In Figs. 5–7 the green diamond represents the arithmetic mean, and the red cross represents the EMR-weighted mean. Measurements more than 1.5 times the interquartile range (IQR) above the upper or below the lower quartile are presented as outliers (open circles). Whiskers represent the outermost values within 1.5 times the IQR of the respective quartiles. The third boxplot represents the spread in EFs from different studies on BB EFs in savannas, and the value that is currently used in large-scale emission assessments. If we investigate the N2O EF intraseasonal variability within the vegetation type classes, we find opposite trends (Table 3). In the open grasslands (campo limpocampo sujo), the weighted-average N2O EF in the EDS was more than double the N2O EF in the LDS. In the open Cerrado (cerrado ralo) and typical Cerrado (cerrado típico), however, the weighted-average N2O EFs were 14 % and 34 % higher in the LDS.

Figure 6EFs (g kg−1) of CO, CH4, and N2O for the various vegetation types. The green diamond represents the arithmetic mean, and the red cross represents the EMR-weighted mean. The purple dot represents the values that are used in GFED for “savanna”, “peat”, and “tropical deforestation” fires, respectively.


Figure 7EFs (g kg−1) for CO, CH4, and N2O for samples of open grassland for different periods since the last fire. The green diamond represents the arithmetic mean, and the red cross represents the EMR-weighted mean.


3.2 EF variability in vegetation type and fire history

We found no significant differences in the MCE, CO EF, and CH4 EF between the EMR-averaged values of the different Cerrado vegetation types, despite substantial differences in tree cover density (Fig. 6). The samples we took over gallery forest contained much higher EFs for CO and CH4, indicating more smouldering combustion. The N2O EF was found to be positively correlated with tree cover and was a factor of 5 times higher in the gallery forest compared to savanna vegetation.

Fire history had some effect on the burning efficiency. We found a decrease in the CO EF and CH4 EF (and thus increase in the MCE) with increasing time between fires ranging from 2 to 4 years in samples from the open grasslands (Fig. 7). Although the measurements in typical Cerrado did not cover the entire fire-frequency span, the available data suggested no significant relation between EFs and the years since the last fire in vegetation of both open Cerrado and typical Cerrado (not shown).

3.3 GWP variability between EDS and LDS fire

Figure 8 shows the cumulative CO2 eq. of the respective gases, based on a 100-year time span. Overall, CO2 eq. emissions per kilogram of dry fuel in the Cerrado were 8.2 % lower in the LDS compared to the EDS. The difference between EDS and LDS CO2 eq. can largely be contributed to somewhat more efficient combustion in the LDS, which is partially compensated for by a higher N2O EF. The black error bar represents the propagation to the net CO2 eq. emissions of the combined standard error of the mean of all species; 12 % to 50 % of this error comes from the propagation of the uncertainty in N2O EFs. Even without taking aerosol effects into account, the indirect radiative forcing due to CO made up a significant portion (45 %–65 %) of total CO2 eq. emissions.

Figure 8CO2 equivalents using GWP with a 100-year horizon and including indirect atmospheric effects for various fire types. The black error bar represents the propagation of the standard error of the mean (SEM) of the combined CO2 eq. emissions.


4 Discussion

4.1 Difference in EFs between EDS and LDS fires

Korontzi et al. (2003b) found that the seasonal curing cycle affected the MCE in prescribed-burn plots in southern African savannas. This intraseasonal shift would limit or even cancel climate benefits of prescribed EDS burning. They found that for “dambo” grasslands, EFs for reduced species were strongly correlated with the percentage of green grass in the fuel. This percentage decreases as grasses cure over the course of the dry season. A similar trend was found by Yokelson et al. (2011) when comparing EF measurements for EDS fires in Mexico to LDS African savanna measurements. Direct measurements taken during the West Arnhem Land Fire Abatement Project (WALFA) in northern Australia, however, showed no significant intraseasonal fluctuation in both CH4 and N2O EFs (Hurst et al., 1994b; Meyer et al., 2012). Measurements taken in Zambian miombo woodlands did not show significant intraseasonal MCE fluctuation either (Hoffa et al., 1999).

In this study we measured EFs during lower-intensity fires in the EDS as well as higher-intensity LDS fires, all in the same region. Although we also found some intraseasonal difference, the decrease of EFs for CO (−15 %) and CH4 (−13 %) was small compared to the −68 % (for CO) and −81 % (for CH4) change found for African grassland fires (Korontzi et al., 2003a). In addition, intraseasonal variability was smaller compared to the variability within EDS or LDS campaigns, and the difference was not statistically significant (p<0.1). The average N2O EF over the samples of combined Cerrado showed a slight increase over the season, though stronger and opposing intraseasonal trends were found in the individual vegetation classes. Meyer et al. (2012) also found opposing intraseasonal N2O EF trends for different vegetation types. While the formation process of N2O is often linked to combustion characteristics (Kilpinen and Hupa, 1991; Meyer et al., 2012; Winter et al., 1999a), we did not find a significant correlation of the N2O EF with the MCE. Overall, MCE values were higher than the average MCE values derived from CO2 and CO EFs for savanna and grassland fires in Andreae (2019) but within the range of previous measurements from Cerrado vegetation (Ferek et al., 1998; Ward et al., 1991). Over all Cerrado vegetation types combined, the weighted-average CH4 EF slightly declined over the dry season.

We conducted the EDS experiments in June when the majority of the prescribed burning takes place (Fig. 2a). Although the LDS measurements in 2018 were taken after the first rains, conditions were still hotter and dryer than during the EDS, and the combustion completeness appeared to be higher (Fig. 3). No fuel moisture measurements were done during the 2018 campaigns, but co-located measurements from 2017 showed limited drying occurring from June to September, with respective average fuel moisture content (FMC) declining from 63.8 % to 55.4 % for live grass and 11.7 % to 7.2 % for dead grass (Santos et al., 2020). Larger differences may be expected earlier in the EDS period of March–May (N'Dri et al., 2018), when the FMC and live-to-dead fuel ratio are even higher (Santos et al., 2020). During these months, when humidity is still very high, prescribed-burning efforts focus on the protection of vulnerable ecosystems such as peatlands and gallery forests, as well as areas around homes and farmlands, but total BA is limited. Additional measurements in the very start of the dry season (March–May) should confirm whether EFs increase for these fires. Rissi et al. (2017) measured fuel characteristics, the rate of spread, flame height, fire intensity (kW m−1), and combustion completeness in campo sujo vegetation (<20 % tree cover) for prescribed burns in May, July, and October. Although the spread in fire intensity between fires was higher in the late season, they found no significant differences in these characteristics between the July and August treatments. Fire intensity was best explained by fuel build-up (Rissi et al., 2017); this is consistent with the MCE increase we found between 2 and 4 years of fuel build-up in vegetation of open grassland (Fig. 7). The finding that the number of years since the last burn did not significantly affect the combustion efficiency after 4 years is consistent with the results from Govender et al. (2006). However, we only found this correlation in open grassland with annual grasses leading to accumulation of easily combustible dead biomass.

4.2 Variability in CO and CH4 EFs

According to our results, there was no significant difference in CO and CH4 EFs between the dominant savanna vegetation types in EESGT: campo limpocampo sujo, cerrado ralo, and cerrado típico. Overall, the weighted-average CO and CH4 EFs for these fuel types of combined savanna were lower than most of the existing literature on savanna fires (Akagi et al., 2011; Andreae, 2019) (Fig. 5). The discrepancy with the literature is particularly strong for CH4 as shown in Fig. 9, where the individual CH4 EF measurements are plotted as a function of the MCE measured for the Cerrado vegetation types. Results from other studies, plotted as the study averages, are shown based on the individual papers included in Andreae (2019). The averaged EFs were rather similar between EDS and LDS campaigns, but within each campaign, the EFs varied substantially. The shift in the LDS towards a steeper slope of the CH4 EF–MCE linear regression in Fig. 9 may be an indication of a shift toward more combustion of woody fuels (Van Leeuwen and Van der Werf, 2011). Although the lower CH4 EF found in this study can partially be explained by on average higher MCE values in our plots, the CH4 EFs were much lower than average CH4 EFs from savanna literature studies with the same MCE. Within the total range of variability, the slopes of the linear regression we found for both EDS and LDS campaigns were significantly less steep compared to the regression slope based on previous measurements of savanna vegetation CH4 EFs. This is to some degree surprising given that the relation between the MCE and CH4 is thought to be well understood. In part, the lower slope comes from a larger number of earlier observations in the 0.90–0.95 MCE range; in the higher MCE ranges our results deviate less from earlier work. This may indicate that there is more variability in fire processes between different savanna types than previously reported. Also compared to earlier measurements from Cerrado vegetation the CH4 EFs were low; Ferek et al. (1998) found an average CH4 EF of 3.7 g kg−1 and Ward et al. (1992) found CH4 EFs ranging from 1 to 1.6 g kg−1. This indicates that more research is needed over ideally a larger range of Cerrados and regions to understand what drives this variability.

Figure 9Relation between the MCE and CH4 EF for all EDS and LDS samples from Cerrado vegetation fires (i.e. excluding humid grasslands and gallery forest samples). Existing savanna measurements are shown using the study-average values in the Andreae (2019) database.


The difference between EDS and LDS weighted-average CH4 EFs is partly the result of a larger spread and high concentration of residual smouldering combustion (RSC) samples in the LDS (Fig. 5). Although the CH4 EF was lower in the LDS (−13 %, Table 3), the overall spread of CH4 EFs in the LDS fires was higher than during EDS fires. Moreover, during the EDS, high CH4 EFs are mostly found in samples with low overall trace gas EMRs (Fig. 10), meaning their impact on the EMR-weighted average was small. An explanation for the increased spread of CH4 EFs in the LDS when the relative humidity was lower may be the effects of more complete combustion of grasses and fine fuels on one hand and an increased share of RSC-prone woody fuels in the fuel mixture leading to a higher CH4 EF on the other hand (Bertschi et al., 2003; Hoffa et al., 1999). These fuels typically contain more moisture in the EDS and are densely packed; therefore, they are more likely to burn in the LDS when humidity is low (Akagi et al., 2011; Eck et al., 2013; van Leeuwen et al., 2014). This is also observed in Australian savannas, where combustion completeness of woody debris was found to be twice as high in the LDS compared to EDS fires (Yates et al., 2015).

Figure 10CH4 EFs as a function of relative humidity based on measurements on the UAV at the time of sampling. The size of the dots represents the ΔCH4 EMR (ppm) in the sample and therefore depicts the relative contribution to the weighted mean. The black diamonds show the weighted-average CH4 EF for each 5 % relative humidity bin. The black line represents the standard error of the class' mean.


Savanna areas with higher tree cover had slightly higher EFs for N2O. Furthermore, there was an opposite intraseasonal trend in N2O EFs from grass-dominated campo limpocampo sujo (−55 % from EDS to LDS) and shrub-dominated cerrado ralo (+14 %) and cerrado típico (+34 %). Winter et al. (1999b) found N2O EFs to be closely correlated with the nitrogen content of the fuel. Susott et al. (1996) and Ward et al. (1992) measured the dry-weight carbon and nitrogen content of various fractions of savanna fuels. For the Cerrado, they analysed dead and living grass, dicots, litter, leaves, and various woody-debris fractions for the most fire-prone Cerrado classes studied in this paper. While they found that carbon content in living grasses was only slightly higher compared to dead grasses, nitrogen content in living grass was on average more than double the content of dead grass. They also found that nitrogen content of leaf, litter, and dicot fractions increased in more woody vegetation types. The nitrogen content of coarse woody debris tends to decrease with the size of the debris. The opposite intraseasonal trends in N2O EFs may therefore be related to a seasonal shift in vegetation types that burn. Many shrubs and trees in EESGT are deciduous or semi-deciduous and drop all or part of their leaves throughout the dry season. This creates a fire-prone, nitrogen-rich litter layer that burns mostly in the LDS fires. In the open grasslands however, where leaf litter is not as significant to the fuel mixture, the ratio of dead versus living grasses increases, which could reduce the nitrogen content of the fuel (Yokelson et al., 2011). The decline found in N2O EF from open grasslands that have not burned for some years (Fig. 7) may thus be related to the increased dead-to-live grass ratio of the fuel mixture as found by Santos et al. (2020). Whether this is indeed the explanation for the opposite intraseasonal trends in N2O emission factors requires future campaigns which include measurements of fuel load, combustion completeness, and nitrogen content over the whole season.

During the LDS, fires can escape into the peatlands and gallery forests lining the rivers. Many EF measurements in the savanna biome are conducted in research plots that are representative of the typical savanna vegetation. These plots, therefore, do not include EFs of these fine-scale landscape features. For this reason, we assessed them separately and have not included them in the Cerrado weighted averages. Intense fires will only occur in these vegetation types in the LDS, when fuels are relatively dry and the groundwater table is low. Late-wet-season management fires in these vulnerable vegetation types are used to reduce moribund fuels. Because we only took a few samples from gallery forest (26 samples) and humid grassland (15 samples), more research is needed in these vegetation types. Based on our measurements in EESGT and the relatively high N : C ratio of these ecosystems as described in literature, the N2O EF of 0.2 g kg−1 currently applied in emission databases both for gallery forest (tropical forest) and humid grasslands (savanna) (Akagi et al., 2011) is likely a significant underestimate for both vegetation types, respectively.

4.3 Uncertainties

The main uncertainties associated with calculating fire-averaged EFs from field measurements include representativeness of the measurements taken related to the sampling strategy, measurement uncertainties, and assumptions based on other literature to represent factors not measured but required to compute EFs.

4.3.1 Sampling strategy

Given the unpredictable nature of fires and difficulties to move around during a spreading fire in a protected area without many roads, we tackled each fire differently. We could not standardize the strategy with regards to sampling head, back, and sideway-propagating fires. Especially in the LDS fires, it was difficult to take many samples over the fast-moving fire front. Therefore, sideway-propagating fires may be overrepresented in the dataset. According to Surawski et al. (2015) based on wind tunnel experiments and Wooster et al. (2011) based on experimental field burns, fire spread mode affects EFs with, in general, a slightly lower MCE occurring in head fires. Compared to earlier studies, we have taken a much larger number of samples thus lowering biases. To better calculate the representative mean requires better-contained fires that are easier to access and continuous sampling at various locations.

During the LDS, fires were predominantly sampled from 11:00 to 16:00 when temperatures are highest. However, these LDS fires generally last for multiple days, and measurements taken during the night and early morning are underrepresented in the dataset. Diurnal fluctuations in temperature, wind, and humidity may cause these fires to behave more similarly to EDS fires during these times. Even though the amount of carbon consumed during those times is presumably lower than during the day, future efforts could shed light on the diurnal cycle of EFs.

An additional source of uncertainty stems from a potential bias related to sampling of RSC conditions. If the sampling period overlapped with the fire duration including the RSC, as was often the case for grasslands, derived EMR values are likely to have been representative. However, as the RSC may persist also after we stopped sampling, especially in more woody fuels, EFs of predominantly RSC products such as CO and CH4 may be underestimated using our sampling strategy. In Fig. 4 the difference between the arithmetic mean (green triangle) and the weighted mean (red square) represents the effect of weighing the bags by the excess mixing ratio. In most cases, the difference is small, suggesting that the total contribution of the RSC is limited. This is consistent with Ward et al. (1992), who measured BB emissions in Cerrado vegetation. They found that over 97 % of the total carbon released was emitted during the flaming phase. The relatively low significance of the RSC in grass-dominated savannas was also found for experiments in the Kruger National Park, South Africa (Cofer et al., 1996; Wooster et al., 2011).

While the role of the RSC in these grass-dominated ecosystems is thus thought to be small, the significance of the RSC in areas with more woody fuel may be higher (Bertschi et al., 2003; Christian et al., 2007; Hao et al., 1991). With prescribed-fire management, dead organic matter and woody carbon stocks may increase over time (Oliveras et al., 2013; Pivello, 2011; Veenendaal et al., 2018). For long-term emission abatement potential, it is therefore important to understand how these changes in fuel composition affect EFs.

4.3.2 EF calculations and assumptions

Ideally, EF calculations are based on measurements of all emitted carbon-containing species. This allows for the direct conversion of emission ratios to EFs per unit of burned fuel. We did not measure non-methane hydrocarbons (NMHCs) and carbonaceous particulate matter (CPM). When combined, these can account for a significant portion of the total carbon emitted. To account for this, we have made assumptions for the CNMHC/CCH4 ratio (3.5, Andreae, 2019), the PM2.5/ CO mass ratio (0.07), and the carbon fraction of PM2.5 (0.70), based on Andreae (2019) and Reid et al. (2005). This adds an additional 0.4 %–2.7 % C from NMHCs and 0.5 %–1.9 % C from CPM to the total carbon balance. Most studies only include carbonaceous trace gases in the total carbon. However, leaving out part of the carbonaceous emissions artificially increases the EFs of the measured species. This inflation is proportional to the carbon that is not accounted for and will likely be in the 1 %–5 % range (Akagi et al., 2011; Yokelson et al., 2013). EFs for both NMHCs and PM are negatively correlated with combustion efficiency (Hoffa et al., 1999; Yokelson et al., 2013). Therefore, the overestimation of EF would be slightly larger in the EDS compared to the LDS. As the N2O EF is coupled to a carbonaceous co-emitted species, in our case CO2, this inflation will also affect the N2O EF.

Another source of uncertainty is the carbon content of the fuel. EFs scale linearly with this fraction, and we used 48 % for vegetation of typical Cerrado and humid grasslands and 50 % for gallery forests (Susott et al., 1996). The carbon content in humid grasslands is based on the assumption that no peat, which has a higher carbon content of ∼56 % (Susott et al., 1996), was combusted in the fire. Had we made other assumptions, for example of 45 % (Andreae, 2019; Andreae and Merlet, 2001) or 50 % (Akagi et al., 2011; Urbanski, 2014), our EF estimates would have been 6 % lower to 4 % higher in types of typical Cerrado and humid grasslands and 10 % lower or equal in gallery forest. This scaling does not affect the spatial and temporal patterns we found.

4.3.3 CO2 eq. calculations and assumptions

Finding a useful metric to assess the direct and indirect impact on RF and climate is challenging, as mechanisms and timeframes of the impact often differ between studies (Fuglestvedt et al., 2010). Atmospheric impact may also depend on geographic location, injection height, or atmospheric conditions (Daniel and Solomon, 1998; Fry et al., 2012). There is substantial uncertainty in GWP (e.g. ±40 % for CH4 GWP100), dominated by uncertainty in the actual GWP for CO2 (i.e. the denominator of the GWP ratio) and inclusion of indirect effects (Myhre et al., 2013). We used an indirect GWP100 for CO of 2.2, i.e. taking into account the CH4 and O3 effects but not considering primary and secondary aerosol effects. Including these effects would increase the effect of CO by 50 % for only primary and 140 % for primary and secondary aerosol effects (Shindell et al., 2009). Due to the short atmospheric lifetimes of CO (2–3 months) and CH4 (12.4 years), using the short-term GWP20 would lead to a ∼3 times higher impact for CO and ∼2.5 times higher impact for CH4 (Myhre et al., 2013). Since we assume sequestration of atmospheric CO2 upon regrowth, the GWP100 values we used for CO (2.2) and CH4 (32) do not include the GWP of CO2 from methane oxidation, which would add roughly (1×MMCO2MMi) to the GWPi, depending on the considered timeframe.

In the savanna biome, fires typically occur frequently with fire return times strongly dependent on the amount of rainfall, hence productivity (Bistinas et al., 2014; Govender et al., 2006; Staver et al., 2011b). As the vegetation recovers after a fire, atmospheric CO2 is captured during photosynthesis, thus balancing CO2 emissions during the fire. This net-zero emission for CO2 is true for stable savanna systems with rapid regrowth, while forest CO2 emissions from fires take longer to be compensated for. For peat underlying humid grasslands, however, some of these emissions might be attributed to carbon that was stored over thousands of years. These carbon stocks will not regenerate at a rate that is relevant to current climate change. As peat layers are still moist in the EDS, the ratio of short-carbon-cycle aboveground fuel to long-carbon-cycle peat may be seasonally dependent. Also, in the case of deforestation, CO2 uptake does not balance out the loss in biosphere carbon stocks due to the fire. Based on our measurements, we cannot conclude whether peat from the soil underlying the humid grasslands contributed to the fuel mixture. If we would not assume CO2 uptake, CO2 eq. EFs would be 453 % and 297 % higher for gallery forest and peat, respectively. Our assumptions to calculate the climate impact of these fires may therefore be seen as conservative and are only valid for stable systems.

4.4 N2O EF uncertainty

N2O EFs were significantly lower than the 0.20 g kg−1 value that is currently used in GFED based on Akagi et al. (2011) and the 0.21 g kg−1 value for savanna in Andreae and Merlet (2001). However, the values we find are more in line with other savanna measurements from South America (0.05–0.07 g kg−1; Hao et al., 1991; Susott et al., 1996), Australia (0.07–0.12 g kg−1; Hurst et al., 1994a; Meyer et al., 2012; Surawski et al., 2015), and Africa (0.16 g kg−1; Cofer et al., 1996). The high average N2O EF in the Akagi et al. (2011) and Andreae and Merlet (2001) databases may partially be linked to the use of stainless-steel sample containers in older studies, leading to N2O formation in the sample container (Muzio and Kramlich, 1988). Due to the low concentrations and small departure from background conditions, N2O is notoriously difficult to measure. Figure 10 shows that many EFs were negative. This occurs when concentrations in the smoke samples were below background concentrations. Although N2O is destroyed in flaming combustion (Winter et al., 1999a, b), and negative emissions are thus theoretically possible, we expect it is more likely to be a measurement error. We found extremely high and low EFs in samples with low overall EMRs. The normal Gaussian distribution pattern in Fig. 11 indicates high measurement uncertainty at low smoke concentrations. The positive and negative tails of this Gaussian error partially balance out, and their weight is low relative to higher concentration measurements. Therefore, the effects of this error on the weighted-average EFs should be limited. Still, a degree of caution is advised while dealing with N2O EFs. The relative error in the 2017 campaign was higher than in the 2018 campaigns due to improvements in the algorithms used to stabilize the CO and N2O sensor implemented after the first campaign. When comparing the same dataset based on vegetation type, a clear shift of the average N2O EF can be found (Fig. 11b). For vegetation types with a low number of measurements or cumulative smoke signal, the large spread reflects much higher uncertainty.

Figure 11N2O EFs plotted against the cumulative EMR of the carbonaceous trace gases in the sample based on (a) all Cerrado measurements in the three campaigns. The black diamonds represent the averages of each 5 mol C bin. (b) Combined EDS and LDS measurements in open grasslands and of vegetation of open Cerrado and typical Cerrado. The dotted lines and numbers on the right represent the weighted-average N2O EFs over all campaigns.


4.5 Limitations of the study

The findings of this study have to be seen in light of some limitations. Field measurements take place in an uncontrolled environment. This means wind conditions vary, possibly affecting the temperature and combustion efficiency of the fire and the type of fuel it consumes. Many processes happen at once during a fire, making it challenging to obtain a representative EF for all stages of the process. Future research will focus on further improving the UAV-based measurement methodology to avoid possible biases as discussed in Sect. 4.2.1. We used generalized vegetation classes based on remote sensing, although we lacked fuel measurements to substantiate or nuance this classification. Although the number of samples taken is substantially higher than earlier campaigns, the sample size for individual categories of “vegetation class” and “years without fire” is in some cases small, meaning we could not always disentangle all different combinations of classes. Measuring more fires, covering a larger geographical area, and adding fuel and wind speed measurements could provide further insights for the variability we found.

5 Conclusions

We obtained over 800 fresh smoke samples in different Cerrado vegetation types, during three fieldwork campaigns at various stages of the fire season. EFs of CO2, CO, CH4, and N2O were calculated from the difference between sample bag and background concentrations based on the carbon mass balance method. Weighted-average EFs over the vegetation of combined Cerrado in the EESGT study region for CO, CH4, and N2O were 48, 0.78, and 0.11 g kg−1, respectively, in the early dry season. In the late dry season, weighted-average EFs were 41 g kg−1 for CO (−15 % compared to early dry season), 0.68 g kg−1 for CH4 (−13 %), and 0.12 g kg−1 for N2O (+17 %). Apart from the intraseasonal N2O EF decrease in grasslands and increase in typical Cerrado, we did not find major intraseasonal EF differences that were statistically significant (p<0.1). Some variability was explained by vegetation type and fire history in open grasslands, whereas relative humidity only had a minor impact on variability. While we found some evidence pointing towards more efficient combustion in the LDS, the difference in weighted-average EFs over the campaigns was low, while the variability during each campaign was substantial. Our findings thus imply that the effectiveness of carbon mitigation in fire abatement programmes is not significantly impacted by intraseasonal changes in EFs for the fieldwork site and length of fire season sampled.

Overall, EFs for CO and CH4 were 36 % and 72 % lower than EFs found in previous studies in the Cerrado and savanna fires in general. The lower CH4 EFs compared to previous studies were not fully explained by the higher MCE but rather by a reduction in the steepness of the slope of the linear regression of CH4 EFs as a function of the MCE. We found that in our study region, N2O EFs for Cerrado vegetation were approximately half the value used in large-scale emission assessments. Uncertainties for N2O measurements are high, especially in low-concentration samples. However, these lower EFs are also found in more recent savanna studies and could indicate a substantial overestimation of the contribution of fires in the N2O budget in global databases. Seasonal effects of N2O EFs were opposite for grass fuels contrasted with more shrub-dominated vegetation types. Finally, our findings indicate that accounting for CO should be considered in carbon schemes. While not a direct greenhouse gas, it has a significant effect on fire radiative forcing through its indirect effect on the CH4 and O3 concentration.

Data availability

Measurement data are available upon request.

Author contributions

RV, GRvdW, and MG designed the study. RV and MAB conducted the experiments. RV, MMC, and GRvdW participated in data analysis and/or interpretation. RV wrote the paper. RV, GRvdW, and ACSB edited the paper.

Competing interests

The authors declare that they have no conflict of interest.


The Chico Mendes Institute for Conservation of Biodiversity (ICMBio) and the Center for Environmental Monitoring and Fire Management (CEMAF) led the fieldwork, which would not have been possible without the rangers working at EESGT. Also, Alan Silva, Eduardo Ganassoli Neto, Micael Moreira, and Jader Nunes Cachoeira were indispensable for providing all logistics related to the fieldwork. We thank Robert Yokelson and Martin Wooster for valuable discussions on emission factor calculation. Finally, we wish to thank Anja Hoffmann for connecting the authors.

Financial support

This research has been supported by the Netherlands Organisation for Scientific Research (NWO) (Vici scheme research programme, grant no. 016.160.324).

Review statement

This paper was edited by Anja Rammig and reviewed by two anonymous referees.


Abreu, R. C. R., Hoffmann, W. A., Vasconcelos, H. L., Pilon, N. A., Rossatto, D. R., and Durigan, G.: The biodiversity cost of carbon sequestration in tropical savanna, Sci. Adv., 3, 1–8,, 2017. 

Akagi, S. K., Yokelson, R. J., Wiedinmyer, C., Alvarado, M. J., Reid, J. S., Karl, T., Crounse, J. D., and Wennberg, P. O.: Emission factors for open and domestic biomass burning for use in atmospheric models, Atmos. Chem. Phys., 11, 4039–4072,, 2011. 

Anderson, R., Beatty, R., Russell-Smith, J., and van der Werf, G. R.: The global potential of indigenous fire management: findings of the regional feasibility assessments, UNU-IAS, Tokyo, available at: (last access: 19 February 2021), 2015. 

Andreae, M. O.: Emission of trace gases and aerosols from biomass burning – an updated assessment, Atmos. Chem. Phys., 19, 8523–8546,, 2019. 

Andreae, M. O. and Merlet, P.: Emission of trace gases and aerosols from biomass burning, Biogeochemistry, 15, 955–966,, 2001. 

Australian Government, Department of the Environment and Energy: Explanatory Statement – Carbon Credits (Carbon Farming Initiative) Act 2011, Australian government, Canbera, available at: Statement/Text (last access: 19 February 2021), 2018. 

Barradas, A. C. S.: A gestão do fogo na Estação Ecológica Serra Geral do Tocantins, Brasil. Mestrado Profissional em Biodiversidade em Unidades de Conservação, master thesis, Instituto de Pesquisas Jardim Botânico do Rio de Janeiro, 2017. 

Barradas, A. C. S., Borges, M. A., and Costa, M. M.: Plano de manejo integrado do fogo – Estação Ecológica Serra Geral do Tocantins, Rio da Conceição, ICMBIO, Rio da Conceição, 2018. 

Barradas, A. C. S., Borges, M. A., Costa, M. M., and Ribeiro, K. T.: Paradigmas da Gestão do Fogo em Áreas Protegidas no Mundo e o Caso da Estação Ecológica Serra Geral do Tocantins, Biodiversidade Bras., 10, 71–86,, 2020. 

Beringer, J., Hutley, L. B., Tapper, N. J., and Cernusak, L. A.: Savanna fires and their impact on net ecosystem productivity in North Australia, Glob. Change Biol., 13, 990–1004,, 2007. 

Bertschi, I., Yokelson, R. J., Ward, D. E., Babbitt, R. E., Susott, R. A., Goode, J. G., and Hao, W. M.: Trace gas and particle emissions from fires in large diameter and belowground biomass fuels, J. Geophys. Res.-Atmos., 108, 1–12,, 2003. 

Bistinas, I., Harrison, S. P., Prentice, I. C., and Pereira, J. M. C.: Causal relationships versus emergent patterns in the global controls of fire frequency, Biogeosciences, 11, 5087–5101,, 2014. 

Chen, L.-W. A., Verburg, P., Shackelford, A., Zhu, D., Susfalk, R., Chow, J. C., and Watson, J. G.: Moisture effects on carbon and nitrogen emission from burning of wildland biomass, Atmos. Chem. Phys., 10, 6617–6625,, 2010. 

Christian, T. J., Yokelson, R. J., Jr, A. C., Griffith, D. W. T., Alvarado, E. C., Santos, C., Gomes, T., Neto, S., Veras, C. A. G., and Hao, W. M.: The tropical forest and fire emissions experiment: Trace gases emitted by smoldering logs and dung from deforestation and pasture fires in Brazil, J. Geophys. Res., 112, 1–14,, 2007. 

Cofer, W. R., Levine, J. S., Winstead, E. L., Cahoon, D. R., Sebacher, D. I., Pinto, P., and Stocks, B. J.: Source compositions of trace gases released during African savanna fires, J. Geophys. Res., 101, 23597–23602, 1996. 

Cook, G. D., Liedloff, A. C., and Murphy, B. P.: Towards a methodology for increased carbon sequestration in dead fuels through implementation of less severe fire regimes in savannas, in: Carbon Accounting and Savanna Fire Management, edited by: Murphy, B. P., Edwards, A. C., Meyer, M., and Russell-Smith, J., CSIRO Publishing, Clayton South, Australia, 321–326, 2015. 

Crutzen, P. J. and Zimmermann, P. H.: The changing photochemistry of the troposphere, Tellus, 43, 136–151,, 1991. 

Daniel, J. S. and Solomon, S.: On the climate forcing of carbon monoxide, J. Geophys. Res.-Atmos., 103, 13249–13260,, 1998. 

DiMiceli, C., Carroll, M., Sohlberg, R., Kim, D., Kelly, M., and Townshend, J.: MOD44B MODIS/Terra Vegetation Continuous Fields Yearly L3 Global 250m SIN Grid V006 [Data set], NASA EOSDIS Land Processes DAAC,, 2015. 

Durigan, G., Ratter, J. A., State, P., and Box, P. O.: The need for a consistent fire policy for Cerrado conservation, J. Appl. Ecol., 53, 11–15,, 2016. 

Eck, T. F., Holben, B. N., Reid, J. S., Mukelabai, M. M., Piketh, S. J., Torres, O., Jethva, H. T., Hyer, E. J., Ward, D. E., Dubovik, O., Sinyuk, A., Schafer, J. S., Giles, D. M., Sorokin, M., Smirnov, A., and Slutsker, I.: A seasonal trend of single scattering albedo in southern African biomass-burning particles: Implications for satellite products and estimates of emissions for the world's largest biomass-burning source, J. Geophys. Res.-Atmos., 118, 6414–6432,, 2013. 

Ferek, R. J., Reid, J. S., Hobbs, P. V., Blake, D. R., and Liousse, C.: Emission factors of hydrocarbons, halocarbons, trace gases and particles from biomass burning in Brazil, J. Geophys. Res., 103, 32107–32118,, 1998. 

Fidelis, A. and Fernanda, M.: Above- and below-ground biomass and carbon dynamics in Brazilian Cerrado wet grasslands, J. Veg. Sci., 24, 356–364,, 2013. 

Fidelis, A., Alvarado, S., Barradas, A., Pivello, V., Fidelis, A., Alvarado, S. T., Barradas, A. C. S., and Pivello, V. R.: The Year 2017: Megafires and Management in the Cerrado, Fire, 1, 49,, 2018. 

Franke, J., Barradas, A. C. S., Borges, M. A., Menezes Costa, M., Dias, P. A., Hoffmann, A. A., Orozco Filho, J. C., Melchiori, A. E., and Siegert, F.: Fuel load mapping in the Brazilian Cerrado in support of integrated fire management, Remote Sens. Environ., 217, 221–232,, 2018. 

Fry, M. M., Naik, V., West, J. J., Schwarzkopf, M. D., Fiore, A. M., Collins, W. J., Dentener, F. J., Shindell, D. T., Atherton, C., Bergmann, D., Duncan, B. N., Hess, P., MacKenzie, I. A., Marmer, E., Schultz, M. G., Szopa, S., Wild, O., and Zeng, G.: The influence of ozone precursor emissions from four world regions on tropospheric composition and radiative climate forcing, J. Geophys. Res.-Atmos., 117, 1–16,, 2012. 

Fuglestvedt, J. S., Shine, K. P., Berntsen, T., Cook, J., Lee, D. S., Stenke, A., Skeie, R. B., Velders, G. J. M., and Waitz, I. A.: Transport impacts on atmosphere and climate: Metrics, Atmos. Environ., 44, 4648–4677,, 2010. 

Giglio, L. and Justice, C. O.: MOD14A1 MODIS/Terra Thermal Anomalies/Fire Daily L3 Global 1km SIN Grid V006, NASA EOSDIS,, 2015. 

Giglio, L., Boschetti, L., Roy, D. P., Humber, M. L., and Justice, C. O.: The Collection 6 MODIS burned area mapping algorithm and product, Remote Sens. Environ., 217, 72–85,, 2018. 

Goodwin, J., Gillenwater, M., Romano, D., and Radunsky, K.: Precursors and indirect emissions, in: 2019 Refinement to the 2006 IPCC Guidelines for National Greenhouse Gas Inventories – General Guidance and Reporting, Volume 1, edited by: Gómez, D. and Irving, W., Institute for Global Environmental Strategies (IGES) for the IPCC, available at: (last access: 10 December 2020), 2019. 

Govender, N., Trollope, W. S. W., and Van Wilgen, B. W.: The effect of fire season, fire frequency, rainfall and management on fire intensity in savanna vegetation in South Africa, J. Appl. Ecol., 43, 748–758,, 2006. 

Hao, W. M., Scharffe, D., Lob, J. M., and Crutzen, P. J.: Emissions of N20 from the Burning of Biomass in an Experimental System, Geophys. Res. Lett., 18, 999–1002, 1991. 

Hoffa, E. A., Ward, D. E., Hao, W. M., Susott, R. A., and Wakimoto, R. H.: Seasonality of carbon emissions from biomass burning in a Zambian savanna, J. Geophys. Res., 104, 13841–13853,, 1999. 

Hurst, D. F., Griffith, D. W. T., Carras, J. C., Williams, D. J., and Fraser, P. J.: Measurements of Trace Gases Emitted by Australian Savanna Fires During the 1990 Dry Season, J. Atmos. Chem., 18, 33–56,, 1994a. 

Hurst, D. F., Griffith, D. W. T., and Cook, G. D.: Trace gas emissions from biomass burning in tropical Australian savannas, J. Geophys. Res., 99, 16441,, 1994b. 

Ito, A. and Penner, J. E.: Global estimates of biomass burning emissions based on satellite imagery for the year 2000, J. Geophys. Res.-Atmos., 109, 1–18,, 2004. 

Jones, M. W., Santín, C., van der Werf, G. R., and Doerr, S. H.: Global fire emissions buffered by the production of pyrogenic carbon, Nat. Geosci., 12, 742–747,, 2019. 

Kauffman, J. B., Cummings, D. L., and Ward, D. E.: Relationships of Fire, Biomass and Nutrient Dynamics along a Vegetation Gradient in the Brazilian Cerrado, J. Ecol., 82, 519–531,, 1994. 

Kilpinen, P. and Hupa, M.: Homogeneous N2O chemistry at fluidized bed combustion conditions: A kinetic modeling study, Combust. Flame, 85, 94–104,, 1991. 

Klink, C. A. and Machado, R. B.: Conservation of the Brazilian Cerrado, Conserv. Biol., 19, 707–713,, 2005. 

Korontzi, S.: Seasonal patterns in biomass burning emissions from southern African vegetation fires for the year 2000, Glob. Change Biol., 11, 1680–1700,, 2005. 

Korontzi, S., Justice, C. O., and Scholes, R. J.: Influence of timing and spatial extent of savanna fires in southern Africa on atmospheric emissions, J. Arid Environ., 395–404,, 2003a. 

Korontzi, S., Ward, D. E., Susott, R. A., Yokelson, R. J., Justice, C. O., Hobbs, P. V., Smithwick, E. A. H., and Hao, W. M.: Seasonal variation and ecosystem dependence of emission factors for selected trace gases and PM2.5 for southern African savanna fires, J. Geophys. Res.-Atmos., 108, 4758,, 2003b. 

Landry, J.-S. and Matthews, H. D.: Non-deforestation fire vs. fossil fuel combustion: the source of CO2 emissions affects the global carbon cycle and climate responses, Biogeosciences, 13, 2137–2149,, 2016. 

Lipsett-Moore, G. J., Wolff, N. H., and Game, E. T.: Emissions mitigation opportunities for savanna countries from early dry season fire management, Nat. Commun., 9, 1–8,, 2018. 

Lobert, J. M. and Warnatz, J.: Emissions from the Combustion Process in Vegetation, in: Fire in the Environment: The Ecological, Atmospheric, and Climatic Importance of Vegetation Fires, edited by: Crutzen, P. J. and Goldammer, J. G., John Wiley and Sons, Hoboken, NJ, USA, 15–38, 1993. 

Maraseni, T. N., Smith, K. R., Griffiths, G., and Apan, A.: Savanna burning methodology for fire management and emissions reduction: a critical review of influencing factors, Carbon Balance Manag., 11, 25,, 2016. 

Melchiori, A. E., Setzer, A. W., Morelli, F., Libonati, R., Cândido, P. D. A., and De Jesús, S. C.: A Landsat-TM/OLI algorithm for burned areas in the Brazilian Cerrado – preliminary results A Landsat-TM/OLI algorithm for burned areas in the Brazilian Cerrado – preliminary results, VII International Conference on Forest Fire Research, edited by: Viegas, D. X., 2014,, 2015. 

Meyer, C. P., Cook, G. D., Reisen, F., Smith, T. E. L., Tattaris, M., Russell-Smith, J., Maier, S. W., Yates, C. P., and Wooster, M. J.: Direct measurements of the seasonality of emission factors from savanna fires in northern Australia, J. Geophys. Res.-Atmos., 117, D20305,, 2012. 

Miranda, H. S., Sato, M. N., Neto, W. N., and Aires, F. S.: Fires in the cerrado, the Brazilian savanna, in: Tropical Fire Ecology, edited by: Mason, J., Springer, Berlin and Heidelberg, Germany, 427–450, 2009. 

Murphy, B. P., Russell-Smith, J., and Prior, L. D.: Frequent fires reduce tree growth in northern Australian savannas: implications for tree demography and carbon sequestration, Glob. Change Biol., 16, 331–343,, 2010. 

Muzio, L. J. and Kramlich, J. C.: An artifact in the measurement of N2O from combustion sources, Geophys. Res. Lett., 15, 1369–1372,, 1988. 

Myhre, G., Shindell, D., Bréon, F.-M., Collins, W., Fuglestvedt, J., Huang, J., Koch, D., Lamarque, J.-F., Lee, D., Mendoza, B., Nakajima, T., Robock, A., Stephens, G., Takemura, T., and Zhang, H.: Anthropogenic and Natural Radiative Forcing, in: Climate Change 2013: The Physical Science Basis. Contribution of Working Group I to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change, edited by: Stocker, T. F., Qin, D., Plattner, G.-K., Tignor, M., Allen, S. K., Boschung, J., Nauels, A., Xia, Y., Bex, V., and Midgley, P. M., Cambridge University Press, Cambridge, United Kingdom and New York, NY, USA, 659 pp., 2013. 

N'Dri, B., Soro, T. D., Gignoux, J., Dosso, K., Koné, M., Dri, J. K. N., Koné, N. A., and Barot, S.: Season affects fire behavior in annually burned humid savanna of West Africa, Fire Ecol., 14, 5,, 2018. 

Oliveira, S. L. J., Campagnolo, M. L., Price, O. F., Edwards, A. C., Russell-Smith, J., and Pereira, J. M. C.: Ecological implications of fine-scale fire patchiness and severity in tropical savannas of northern Australia, Fire Ecol., 11, 10–31,, 2015. 

Oliveras, I., Meirelles, S. T., Hirakuri, V. L., Freitas, C. R., Miranda, H. S., and Pivello, V. R.: Effects of fire regimes on herbaceous biomass and nutrient dynamics in the Brazilian savanna, Int. J. Wildl. Fire, 22, 368–380,, 2013. 

Orozco-Filho, J. C.: Avaliação do uso da abordagem orientada-objeto com imagens de alta resolução RapidEye na classificação das fitofisionomias do cerrado, Universidade de Brasília, Brasilia, available at: (last access: 19 February 2021), 2017. 

Pellegrini, A. F. A., Ahlström, A., Hobbie, S. E., Reich, P. B., Nieradzik, L. P., Staver, A. C., Scharenbroch, B. C., Jumpponen, A., and Anderegg, W. R. L.: Fire frequency drives decadal changes in soil carbon and nitrogen and ecosystem productivity, Nature, 553, 194–198,, 2018. 

Penman, T. D., Christie, F. J., Andersen, A. N., Bradstock, R. A., Cary, G. J., Henderson, M. K., Price, O., Tran, C., Wardle, G. M., Williams, R. J., and York, A.: Prescribed burning: how can it work to conserve the things we value?, Int. J. Wildl. Fire, 20, 721,, 2011. 

Penner, J. E., Quaas, J., Storelvmo, T., Takemura, T., Boucher, O., Guo, H., Kirkevåg, A., Kristjánsson, J. E., and Seland, Ø.: Model intercomparison of indirect aerosol effects, Atmos. Chem. Phys., 6, 3391–3405,, 2006. 

Pivello, V. R.: The use of fire on the cerrado and amazonian rainforrests, Fire Ecol., 7, 24–39,, 2011. 

Price, O. F., Russell-Smith, J., and Watt, F.: The influence of prescribed fire on the extent of wildfire in savanna landscapes of western Arnhem Land, Australia, Int. J. Wildl. Fire, 21, 297–305,, 2012. 

Ramos-Neto, M. B. and Pivello, V. R.: Lightning Fires in a Brazilian Savanna National Park: Rethinking Management Strategies, Environ. Manage., 26, 675–684,, 2000. 

Reid, J. S., Koppmann, R., Eck, T. F., and Eleuterio, D. P.: A review of biomass burning emissions part II: intensive physical properties of biomass burning particles, Atmos. Chem. Phys., 5, 799–825,, 2005. 

Ribeiro, J. F. and Walter, B. M. T.: As principais Fitofisionomias do Bioma Cerrado, in: Cerrado: ecologia e flora, Embrapa cerrados, Brasilia, 152–212, 1998. 

Rissi, M. N., Baeza, J., Gorgone-barbosa, E., Zupo, T., and Fidelis, A.: Does season affect fire behaviour in the Cerrado?, Int. J. Wildl. Fire, 26, 427–433,, 2017. 

Russell-Smith, J., Cook, G. D., Cooke, P. M., Edwards, A. C., Lendrum, M., Meyer, C. P., and Whitehead, P. J.: Managing fire regimes in north Australian savannas: Applying Aboriginal approaches to contemporary global problems, Front. Ecol. Environ., 11, e55–e63,, 2013. 

Santos, M. M., Batista, A. C., Silva, A. D. P. da, Ganassoli Neto, E., Barradas, A. C. S., and Giongo, M.: Characterization and Dynamics of Surface Fuel of Cerrado Grassland in Jalapão Region – Tocantins, Brazil, Floresta, 51, 127,, 2020. 

Schmidt, I. B., Moura, L. C., Ferreira, M. C., Eloy, L., Sampaio, A. B., Dias, P. A., and Berlinck, C. N.: Fire management in the Brazilian savanna: First steps and the way forward, J. Appl. Ecol., 55, 2094–2101,, 2018. 

Scott, K., Setterfield, S. A., Douglas, M. M., Parr, C. L., Schatz, J. O. N., and Andersen, A. N.: Does long-term fire exclusion in an Australian tropical savanna result in a biome shift? A test using the reintroduction of fire, Austral Ecol., 37, 693–711,, 2012. 

Seiler, W. and Crutzen, P. J.: Estimates of Gross and Net Fluxes of Carbon Between the biosphere and the athmosphere from biomass burning, Clim. Change, 2, 207–247, 1980. 

Seplan: Zoneamento Ecológico Econômico de Tocantins, available at: (last access: 10 December 2020), 2003. 

Sexton, J. O., Song, X. P., Feng, M., Noojipady, P., Anand, A., Huang, C., Kim, D. H., Collins, K. M., Channan, S., DiMiceli, C., and Townshend, J. R. G.: Global, 30-m resolution continuous fields of tree cover: Landsat-based rescaling of MODIS Vegetation Continuous Fields with lidar-based estimates of error, Int. J. Digit. Earth, 6, 427–448,, 2013. 

Shindell, D. T., Faluvegi, G., Koch, D. M., Schmidt, G. A., Unger, N., and Bauer, S. E.: Improved Attribution of Climate Forcing to Emissions, Interactions, 326, 716–718, 2009. 

Staver, A. C., Archibald, S., and Levin, S. A.: The Global Extent and Determinants of Savanna and Forest as Alternative Biome States, Science, 334, 230–232,, 2011a. 

Staver, A. C., Archibald, S., and Levin, S.: Tree cover in sub-Saharan Africa: Rainfall and fire constrain forest and savanna as alternative stable states, Ecology, 92, 1063–1072,, 2011b. 

Sudo, K. and Akimoto, H.: Global source attribution of tropospheric ozone: Long-range transport from various source regions, J. Geophys. Res.-Atmos., 112, D12302,, 2007. 

Surawski, N. C., Sullivan, A. L., Meyer, C. P., Roxburgh, S. H., and Polglase, P. J.: Greenhouse gas emissions from laboratory-scale fires in wildland fuels depend on fire spread mode and phase of combustion, Atmos. Chem. Phys., 15, 5259–5273,, 2015. 

Surawski, N. C., Sullivan, A. L., Roxburgh, S. H., Meyer, C. P. M., and Polglase, P. J.: Incorrect interpretation of carbon mass balance biases global vegetation fire emission estimates, Nat. Commun., 7, 1–5,, 2016. 

Susott, R. A., Olbu, G., Baker, S. P., Ward, D. E., Kauffman, J. B., and Shea, R. W.: Carbon, hydrogen, nitrogen, and thermogravimetric analysis of tropical ecosystem biomass, in: Biomass Burning and Global Change: Remote sensing, modeling and inventory Development and Biomass Burning in Africa, edited by: Levine, J. S., MIT Press, Cambridge, MA, USA, 350–360, 1996. 

Urbanski, S.: Forest Ecology and Management Wildland fire emissions, carbon, and climate: Emission factors, Forest Ecol. Manage., 317, 51–60,, 2014. 

Urbanski, S. P.: Combustion efficiency and emission factors for wildfire-season fires in mixed conifer forests of the northern Rocky Mountains, US, Atmos. Chem. Phys., 13, 7241–7262,, 2013. 

van der Werf, G. R., Randerson, J. T., Giglio, L., van Leeuwen, T. T., Chen, Y., Rogers, B. M., Mu, M., van Marle, M. J. E., Morton, D. C., Collatz, G. J., Yokelson, R. J., and Kasibhatla, P. S.: Global fire emissions estimates during 1997–2016, Earth Syst. Sci. Data, 9, 697–720,, 2017. 

van Leeuwen, T. T. and van der Werf, G. R.: Spatial and temporal variability in the ratio of trace gases emitted from biomass burning, Atmos. Chem. Phys., 11, 3611–3629,, 2011. 

van Leeuwen, T. T., van der Werf, G. R., Hoffmann, A. A., Detmers, R. G., Rücker, G., French, N. H. F., Archibald, S., Carvalho Jr., J. A., Cook, G. D., de Groot, W. J., Hély, C., Kasischke, E. S., Kloster, S., McCarty, J. L., Pettinari, M. L., Savadogo, P., Alvarado, E. C., Boschetti, L., Manuri, S., Meyer, C. P., Siegert, F., Trollope, L. A., and Trollope, W. S. W.: Biomass burning fuel consumption rates: a field measurement database, Biogeosciences, 11, 7305–7329,, 2014. 

Veenendaal, E. M., Torello-Raventos, M., Miranda, H. S., Sato, N. M., Oliveras, I., Langevelde, F., Van, Asner, G. P., Lloyd, J., and Lloyd, J.: On the relationship between fire regime and vegetation structure in the tropics, New Phytol., 218, 153–166,, 2018. 

Ward, D. E., Setzer, A. W., Kaufman, Y. J., and Rasmussen, R. A.: Characteristics of smoke emissions from biomass fires of the Amazon region – BASE-A experiment, in: Global Biomass Burning – Atmospheric, Climatic and Biospheric Implications, edited by: Levine, J. S., The MIT press, Cambridge, MA, USA, 394–401, 1991. 

Ward, D. E., Susott, R. A., Kauffman, J. B., Babbitt, R. E., Cummings, D. L., Dias, B., Holben, B. N., Kaufman, Y. J., Rasmussen, R. A., and Setzer, A. W.: Smoke and fire characteristics for cerrado and deforestation burns in Brazil: BASE-B Experiment, J. Geophys. Res.-Atmos., 97, 14601–14619,, 1992. 

Winter, F., Wartha, C., and Hofbauer, H.: NO and N2O formation during the combustion of wood, straw, malt waste and peat, Bioresour. Technol., 70, 39–49,, 1999a. 

Winter, F., Wartha, C., and Hofbauer, H.: The Relative Importance of Radicals on the N2O and NO Formation and Destruction Paths in a Quartz CFBC, J. Energy Resour. Technol., 121, 131–136,, 1999b. 

Wooster, M. J., Freeborn, P. H., Archibald, S., Oppenheimer, C., Roberts, G. J., Smith, T. E. L., Govender, N., Burton, M., and Palumbo, I.: Field determination of biomass burning emission ratios and factors via open-path FTIR spectroscopy and fire radiative power assessment: headfire, backfire and residual smouldering combustion in African savannahs, Atmos. Chem. Phys., 11, 11591–11615,, 2011. 

Yates, C. P., Russell-Smith, J., Murphy, B. P., Desailly, M., Evans, J., Legge, S., Lewis, F., Lynch, D., and Edwards, A. C.: Fuel accumulation, consumption and fire patchiness in the lower rainfall savanna region, in: Carbon Accounting and Savanna Fire Management, edited by: Murphy, B. P., Edwards, A. C., Meyer, M., and Russell-Smith, J., CSIRO Publishing, Clayton South, Australia, 115–127, 2015. 

Yokelson, R. J., Goode, J. G., Ward, D. E., Susott, R. A., Babbitt, R. E., Wade, D. D., Bertschi, I., Griffith, D. W. T., and Hao, W. M.: Emissions of formaldehyde, acetic acid, methanol, and other trace gases from biomass fires in North Carolina measured by airborne Fourier transform infrared spectroscopy, J. Geophys. Res.-Atmos., 104, 30109–30125,, 1999. 

Yokelson, R. J., Burling, I. R., Urbanski, S. P., Atlas, E. L., Adachi, K., Buseck, P. R., Wiedinmyer, C., Akagi, S. K., Toohey, D. W., and Wold, C. E.: Trace gas and particle emissions from open biomass burning in Mexico, Atmos. Chem. Phys., 11, 6787–6808,, 2011.  

Yokelson, R. J., Burling, I. R., Gilman, J. B., Warneke, C., Stockwell, C. E., de Gouw, J., Akagi, S. K., Urbanski, S. P., Veres, P., Roberts, J. M., Kuster, W. C., Reardon, J., Griffith, D. W. T., Johnson, T. J., Hosseini, S., Miller, J. W., Cocker III, D. R., Jung, H., and Weise, D. R.: Coupling field and laboratory measurements to estimate the emission factors of identified and unidentified trace gases for prescribed fires, Atmos. Chem. Phys., 13, 89–116,, 2013. 

Short summary
We used drones to measure greenhouse gas emission factors from fires in the Brazilian Cerrado. We compared early-dry-season management fires and late-dry-season fires to determine if fire management can be a tool for abating emissions. Although we found some evidence of increased CO and CH4 emission factors, the seasonal effect was smaller than that found in previous studies. For N2O, the third most important greenhouse gas, we found opposite trends in grass- and shrub-dominated areas.
Final-revised paper