Skip to main content

Carbon dioxide and particulate emissions from the 2013 Tasmanian firestorm: implications for Australian carbon accounting



Uncontrolled wildfires in Australian temperate Eucalyptus forests produce significant smoke emissions, particularly carbon dioxide (CO2) and particulates. Emissions from fires in these ecosystems, however, have received less research attention than the fires in North American conifer forests or frequently burned Australian tropical savannas. Here, we use the 2013 Forcett–Dunalley fire that caused the first recorded pyrocumulonimbus event in Tasmania, to understand CO2 and particulate matter (PM2.5) emissions from a severe Eucalyptus forest fire. We investigate the spatial patterns of the two emissions using a fine scale mapping of vegetation and fire severity (50 m resolution), and utilising available emission factors suitable for Australian vegetation types. We compare the results with coarse-scale (28 km resolution) emissions estimates from Global Fire Emissions Database (GFED) to determine the reliability of the global model in emissions estimation.


The fine scale inventory yielded total CO2 emission of 1.125 ± 0.232 Tg and PM2.5 emission of 0.022 ± 0.006 Tg, representing a loss of 56 t CO2 ha−1 and 1 t PM2.5 ha−1. The CO2 emissions were comparable to GFED estimates, but GFED PM2.5 estimates were lower by a factor of three. This study highlights the reliability of GFED for CO2 but not PM2.5 for estimating emissions from Eucalyptus forest fires. Our fine scale and GFED estimates showed that the Forcett–Dunalley fire produced 30% of 2013 fire carbon emissions in Tasmania, and 26–36% of mean annual fire emissions for the State, representing a significant single source of emissions.


Our analyses highlight the need for improved PM2.5 emission factors specific to Australian vegetation, and better characterisation of fuel loads, particularly coarse fuel loads, to quantify wildfire particulate and greenhouse gas emissions more accurately. Current Australian carbon accountancy approach of excluding large wildfires from final GHG accounts likely exaggerates Tasmania’s claim to carbon neutrality; we therefore recommend that planned and unplanned emissions are included in the final national and state greenhouse gas accounting to international conventions. Advancing these issues is important given the trajectory of more frequent large fires driven by anthropogenic climate change.


Fire plays an important role in the functioning of many terrestrial ecosystems globally and affects climate via the release of greenhouse gases (GHGs) and aerosols in smoke. Emerging evidence suggests that climate change is causing worsening fire weather, longer fire seasons and more intense wildfires globally [1]. Frequent and intense fires have the potential to release enormous quantities of greenhouse gases, thereby exacerbating climate change in a positive feedback process. Carbon dioxide (CO2) contributes the largest proportion of total wildfire smoke emissions (90% of carbon emissions) and is therefore an important driver of radiative forcing [2]. CO2 is assimilated by plants in subsequent growing seasons post-fire; however, frequent fires and changing climate may limit the ability of ecosystems to recover from the fires, resulting in net positive CO2 emissions [3]. Another important product of wildfire combustion is particulate emission which accounts for < 5% of total carbon emissions [4]. Smoke particles affect climate in complex and poorly understood ways causing both short term regional climate cooling due to regional haze formation [5], somewhat analogous to volcanic eruptions [6], and also atmospheric warming, affecting precipitation patterns [7]. Particulates (especially PM2.5, the fraction of particles with a diameter < 2.5 µm) have an important and demonstrable harmful effects on human health, including worsened respiratory symptoms, exacerbation of respiratory and cardiovascular diseases, and premature mortality from cardiovascular complications [8].

These issues are well illustrated by fire activity in Australian temperate forests that have experienced increased fire danger due to extreme fire weather conditions, with resultant lengthening of fire seasons earlier into spring months, associated with climate change [9]. Further, the recent 2019–2020 Black Summer fires in south-eastern Australia are historically unprecedented and most likely exacerbated by climate change [10,11,12]. Analyses involving remote sensing of atmospheric chemistry suggest that the Black Summer fires emitted 715 Tg of CO2 [13], in broad agreement with a bootstrapped emissions estimate of c 670 Tg [3]. It is estimate that 0.3–1.1 Tg of smoke particles were injected into the stratosphere by these fires [14]. Associated particulate pollution from the 2019–2020 fires is estimated to have caused premature death of 429 people and caused nearly 2 billion Australian dollars in health costs [15]. The emissions for the 2019–2020 season are estimated to be 80 times higher than the average fire season apparent in the satellite record [16], highlighting the importance of understanding the impacts of wildfires on GHG emissions.

Despite their capacity to pollute the atmosphere, there are surprisingly few studies of carbon and particulate emission from individual Australian fires. Savanna fires in northern Australia have received the greatest attention, motivated by interest in landscape carbon abatement programs, e.g., [17,18,19]. In temperate Eucalyptus forests, the majority of the studies are based on emissions from prescribed fires, e.g., [20,21,22] with a few exceptions involving laboratory measurements, e.g., [23] or wildfires, e.g., [24, 25]. Particulate emissions from Australian fires still remain largely unexplored, with one study conducted from prescribed fires in south-eastern Australia [26] and a second study on Black Summer fires [14]. A consequence of this limited inquiry is that global analyses often extrapolate these few studies to the entire Australian continent, or use gaseous and particulate emission coefficients from other biomes globally, especially North America, or both. For instance, a frequent source of emissions data is the Global Fire Emissions Database (GFED), which is the most widely used global emissions inventory and has also been critical in assessing the global and regional burden of mortality due to PM2.5 pollution from landscape fires [27, 28].

Accurate estimation of carbon emissions is important for a complete understanding of regional and national carbon accounts. Emissions from Australian wildfires are accounted for in the national GHG accounting to the Intergovernmental Panel on Climate Change (IPCC); however, very large fires are attributed as natural disturbances, so they are excluded in the final fire-related emissions estimation [29]. This approach likely affects the claim of ‘carbon neutrality’ by the state of Tasmania given a spate of large wildfires that have burned around 25% of the island since 1990.

Accurate assessments of particulate emissions are essential for quantifying the exposure of populations to smoke pollution, and in assessing the trade-offs in health impacts from prescribed fires and wildfires [30]. Beyond substantial health costs, particulates have, as aforementioned, a demonstrable harmful impact on human health with regard to cardiovascular and respiratory complications [27]. Particulate emission estimation is also important from a climate perspective because of their influence on haze and cloud dynamics that affects atmospheric chemistry and radiative balance at regional and hemispherical scales [31].

The January 2013 Forcett–Dunalley fire presents an ideal model system to understand smoke emissions from a single, intense fire in a southeast Australian temperate Eucalyptus forest. This fire is notable because it generated a pyrocumulonimbus (PyroCb)—a fire-induced thunderstorm that almost destroyed a small town [32], with the fire burning 25,950 ha of natural Eucalyptus forests, Eucalyptus plantations and agricultural lands. Approximately 55% of the area burnt as high-very high severity, under the influence of extreme weather and dry fuels in the landscape, coupled with a conducive undulating terrain that amplified the fire intensity, estimated to reach c. 68,000 kW m−1 [33]. The PyroCb from the fire was the first record for the island state, although it is becoming increasingly common across eastern Australia and in North America, likely due to climate change [14].

In this study, we test the hypothesis that CO2 and PM2.5 emissions from a single intense wildfire that were estimated from an existing geographically coarse-scale global model is closely correlated with estimates from a purpose-built local model using spatially high-resolution inputs. We then explore how global and local scale wildfire emission estimates can improve regional and national carbon accounting approaches and thereby shape the understanding of carbon ‘costs’ of wildfires. Building on prior analyses of the fire [33] and the bootstrapped emissions analysis of Bowman et al. [3], we: (1) use fine-scale mapping of vegetation and fire severity to map the spatial distribution of CO2 and fine particulate matter (PM2.5) emissions from the fire using the original model by Seiler and Crutzen [34]; (2) compare the spatial distribution and total emissions of the two pollutants between the basic model and the global GFED model to determine the effect of geographical resolution of fire severity and vegetation mapping on emissions estimation; (3) compare daily emissions estimates between the two inventories during the days of concurrently recorded fire activity (3–18 January 2013); and (4) contextualise the Forcett–Dunalley emissions and determine by how much the emissions contribute to overall wildfire and GHG emissions in Tasmania. This study is limited to estimation of CO2; estimation of additional gaseous species such as methane and nitrous oxides could in future be scaled beyond CO2 and expressed as CO2-equivalent emissions.


Study area

The Forcett–Dunalley fireground has a cool temperate climate with annual rainfall of 700–1000 mm, mean daily temperature of 17 °C in summer and 9 °C in winter, and elevation rising from sea-level to 600 m above sea level (Fig. 1b). Native Eucalyptus forests, and Pinus and Eucalyptus plantations are found within the area, with the dry Eucalyptus forest as the most dominant vegetation type (Fig. 1c). The fire occurred from 3 to 18 January 2013 on the Forestier and Tasman Peninsulas in the south-east of Tasmania, the southernmost island state of Australia (Fig. 1a). The fire was ignited possibly from a smouldering stump from an unextinguished campfire. The fire burnt under varying fire weather conditions, topography and fuel characteristics leading to spatial variability in fire severity within the fireground (Fig. 1d). By the time of containment, the fire had burnt approximately 20,200 ha of the 25,950 ha fireground, mostly affecting native vegetation and rural lands (Fig. 1c). A detailed description of the fire and associated broader environmental conditions have been provided in Ndalila et al. [33].

Fig. 1
figure 1

Adapted from Ndalila et al. [33]

Location of the Forcett–Dunalley fireground in SE Tasmania: a Annual rainfall (in mm) and elevation (in m) across Tasmania and the location of major fires in the 2013 fire season including Forcett–Dunalley (1). b Elevation and mean annual rainfall across the Forestier and Tasman Peninsulas, derived from Worldclim dataset [35]. The location of Dunalley township is indicated on the map. c Dominant vegetation in the Forestier and Tasman Peninsulas based on TASVEG 3.0, an integrated vegetation map of Tasmania. d Fire severity patterns within the fireground.

Data preparation

Emission factors (EFs) in our study represent the total mass of a series of gaseous or particulate species emitted per mass of dry fuel burnt. In order to calculate total emissions from biomass burning over a defined area, emission factors are multiplied by the mass of fuel consumed, in a relationship defined by Eq. 1 [34]. The equation incorporates emissions factors (EFs) for the emitted gases and particulates, in addition to the standard estimates of area burnt, fuel loads and the fraction of fuel consumed. A grid covering the extent of the fire perimeter, with 50 m-resolution grid cells, was used for the emissions analysis.

$$\mathrm{E}i =\mathrm{ A }(x) *\mathrm{ FL }(x) *\mathrm{ CC }*\mathrm{ EF}i$$

where, Ei is mass (in g) of emitted species i; A is area burned (in m2) at grid cell x; FL is total fuel load (in kg m−2) at grid cell x; CC is combustion completeness (or the fraction of consumed fuel, 0–1 scale); and EFi is the emission factor (in g kg −1) of the chemical species i.

Our study focused on the spatiotemporal variability of CO2 and PM2.5 emissions given their crucial role in regulating the earth’s carbon and energy budget, and the latter influencing human health; as such, other CO2-equivalent gases (methane and nitrous oxide) were not considered. Likewise, the choice of PM2.5 over PM10 was guided by the fact that for biomass combustion emissions, PM2.5 makes up the majority of PM10, and is more damaging to human health than PM10. The pollutant can penetrate the lungs and be transported to other organs through the bloodstream and trigger reactions such as bronchitis, asthma attacks, cardiovascular diseases and premature mortality [36].

Area burnt records, obtained from Tasmania Fire Service, included unburnt patches within the fire perimeter. These were excluded from the analysis so that only burnt areas remained, covering 20,200 ha of the perimeter. Fuel load estimates representative of all vegetation within the perimeter were absent except for one site that was sampled after the fire from paired burnt-unburnt plots [37]. We therefore adopted fuel load estimates (in t ha−1 dry matter) across Tasmania and from literature on southern Eucalyptus forests of Australia (Table 1). Since a large variability of fuel loads existed across different regions in Australia, emissions calculation involved a bootstrapping of all available ranges of fuel load within each vegetation class to account for the uncertainties propagated by fuel loads. Fuel loads in Table 1 have been stratified into fine (diameter < 0.6 cm) and coarse woody debris (CWD, diameter > 0.6 cm), where fine fuels represent surface to elevated fuels (e.g., litter, standing herbs, grass and fine twigs), while CWD represents fallen twigs, branch wood, logs and stumps.

Table 1 The variability of fuel load (in t ha−1 of dry matter) within the general southern Australia Eucalyptus forests

Since CWD fuel loads identified in the literature included outlying extreme values in the native forests, some of which were obtained following logging operations and included exaggerated coarse debris, we decided to use variability of mean values of fine and coarse fuels within each vegetation type to limit the influence of these outlying values. We conducted 100 simulations where within each run, all grid cells for a given vegetation type were assigned the same random fuel load value drawn from a uniform distribution from the available range of mean fuel load values (Table 1). For example, in any one simulation, all cells within the dry forest class were assigned a similar fine fuel value between 9 and 21 t ha−1, and CWD value between 16 and 74 t ha−1, with the values changing for every simulation so that at the end, 100 emissions estimates are produced. The fuel load values were converted to kg m−2 and aggregated to 2500 m2 to harmonise all analyses at a 50 m × 50 m grid cell scale.

Combustion completeness (fraction of fuel burnt) was determined based on a combination of previous fire severity mapping for this study area [33] and field measurements of fuel consumption in prescribed and wildfires in Eucalyptus-dominated forests in Tasmania and south-eastern Australia (Table 2). We chose these data sources to estimate fuel consumption because field measurements of consumption after the Forcett–Dunalley fire were largely lacking. We partitioned fuel consumption according to severity classes mapped from the Forcett–Dunalley fire based on the assumption that areas with high fire severity have most (or all) of the fine, coarse dead fuels and canopy burnt while for areas that burnt under mild severity, a lower fraction of the fuel mass is consumed (Table 2). The CWD combustion estimates in Table 2 concur with woody fuel consumption estimates reported by Hollis et al. [38] in two high-severity fires: the Kilmore East fire and the Pickering Brook fire. Fire patchiness, which is usually incorporated in the estimation of combustion efficiency, was assumed in this study to be accounted for by the high spatial resolution of the severity mapping. Therefore, patchiness at a resolution below that of the pixel dimensions was not considered. No differentiation in fuel consumption is made between different woody vegetation classes (native or plantation forest). We acknowledge the lack of site-specific fuel consumption also introduces uncertainties in estimation of emissions [39, 40].

Table 2 Estimates of consumed biomass per fuel size class and fire severity (dNBR) class for native (dry and wet Eucalyptus forests) and plantation (Pinus and Eucalyptus) forests, obtained from previous field measurements of native forests in Tasmania and mainland Australia

Lastly, we adopted emission factors for CO2 and PM2.5 from literature based on lab analysis and previous prescribed burning campaigns in southern Australian Eucalyptus forests (Table 3). Emissions factors have not been partitioned into different vegetation classes because estimates are lacking in most classes found in the study area.

Table 3 Emission factors (in g kg−1) for CO2 and PM2.5 for fine and coarse fuels as used in Southern Australian Eucalyptus-dominated landscapes

Spatiotemporal distribution of emissions

The spatial distribution of emissions was determined by combining the aforementioned model variables in Eq. 1 [34] using R version 3.6.1 [46] and ArcGIS 10.3 [47]. The fine scale approach (using 50 m grid resolution) followed a schematic workflow (Fig. 2), which includes the mentioned input variables in Eq. 1. A feature of this analysis is the use of detailed fire severity information and vegetation mapping to estimate emissions. Maps of the spatial distribution of emissions of both CO2 and PM2.5 were produced, where estimates in each grid cell were totals from emissions values for both fine and coarse fuels. Total emission for each pollutant was determined for each of the 100 runs by summing values from all grid cells. We then obtained a bootstrapped mean and standard deviation of total emissions across the runs.

Fig. 2
figure 2

Systematic flowchart of emissions analysis from the Forcett–Dunalley fire, with inputs obtained from available geospatial datasets, previous field assessments and literature

A daily variation of these emissions was determined by intersecting the final emissions map with the fire progression isochrones and summing emissions contained within each temporal polygon. It is worth noting that at the start of the fire, the fire spread polygons were available at sub-daily intervals but as the fire progressed, the time interval between available boundary mapping increased to day(s). We therefore aggregated emissions from sub-daily resolution to daily progressions by combining all emissions for each day.

Comparison with GFED inventory

To assess the reliability of a global emissions model (GFED) in situations of unavailability of site-specific fire data, we compared the spatial and temporal variability of CO2 and PM2.5 emissions between the above fine scale analysis and the GFED4 inventory for January 2013. GFED4 is an industry-standard global emissions model that provides 3-hourly, daily and monthly estimates of 42 emissions species from across the globe at 0.25° (~ 28 km) spatial resolution from the year 1997 [48]. GFED is based on a Carnegie–Ames–Stanford Approach (CASA) biogeochemical model that simulates carbon fluxes from satellite-based observations of vegetation, weather, area burnt and combustion completeness. A full description of the model is provided in van der Werf et al. [48].

We downloaded two gridded datasets (combusted dry matter (DM), and the area burnt layer for January 2013) from the GFED website [49], and multiplied the variables with recommended GFED emission factors for temperate forests (12.9 and 1647 g kg −1 for PM2.5 and CO2 respectively). The result was a spatial map of the two emissions for the entire Tasmania, and a monthly estimate for January 2013 for specific cells that represent the Forcett–Dunalley fireground. These monthly estimates were partitioned into daily emissions by using a daily fraction file that contains the contribution of each grid cell to the total emissions. The daily and spatial variation of the resulting maps from the fine scale and GFED inventories were quantitatively and visually compared to determine the effect of the geographic resolution of fire severity and vegetation mapping on emissions. It should be acknowledged that the daily GFED estimates were only available for 3–14 January, which coincide with the duration of MODIS thermal hotspots data available for the study area. It is therefore likely that the burnt area layer was obtained from a combination of spectral reflectance of burnt area and thermal hotspot data, the latter of which is adopted in GFED4 to represent small fires that would have been missed in previous GFED versions.

We then validated the two emissions inventories using FullCAM simulation of carbon emission (which can be converted to CO2 via 3.67 factor) over the Forcett–Dunalley fireground. FullCAM is a modelling interface used in Australian GHG accounting of the land sector [50], and can simulate fire emissions as an event by feeding in carbon flux estimates from combustion of forest debris and live biomass. Major emission outputs of the model include methane, nitrous oxide and carbon. To determine carbon emissions within FullCAM, we used input parameters values recommended in Surawski et al. [50] for wildfire events with fire intensities of > 7000 kW m−1 in which trees have not been killed.

Contextualising emissions in Tasmania

To gauge the relative contributions of the Forcett–Dunalley fire (that included a significant PyroCb event) to typical annual fire emissions in the state, we compared the Dunalley emissions with the mean fire emission estimates for Tasmania for the period 1997–2020 (the period of the available GFED record). First, we merged GFED estimates across the different vegetation types in Tasmania to produce an annual emission estimate for the above period. Since the GFED emissions were available as carbon emissions, for comparison with estimates from Dunalley fire, we converted GFED’s carbon emissions estimates to CO2 (using 3.67 conversion factor). The percentage of Forcett–Dunalley emissions was then estimated relative to: (1) the total 2013 fire emissions across the state, and (2) mean annual fire emissions for the state. We then examined Tasmania’s fire emissions relative to the state wide carbon (GHG) emissions budget, in order to quantify the effect of excluding severe fires from GHG accounting under the assumption that the fires are natural disturbances and beyond human control.


Spatial distribution of emissions

From the fine scale emissions inventory, total CO2 emissions were 1.125 ± 0.232 million tonnes (or 1.125 ± 0.232 Tg), translating to 55.7 t ha−1 of CO2 released from the 20,200-ha burnt area (Table 4). PM2.5 emissions reached 0.022 ± 0.006 Tg and 1.1 t ha−1 when normalized by area burnt. Carbon dioxide emissions varied across the fireground, reaching 33 tonnes per 50 m resolution grid cell, while the PM2.5 emission peaked at 0.72 tonnes (Fig. 3). It is worth noting that the spatial patterns of both CO2 and PM2.5 are identical because they are based on the same amount of consumed fuel per unit area, but only differ in their respective emissions factors. In both pollutants, the highest emissions were in the south-southwest of the fireground, characterized by the highest fire severity classes (Fig. 1d). These areas also coincided with a large flaming zone in the classified infra-red linescan map for 4 January (see Additional file 1: Fig. S1) which was associated with elevated fire weather.

Table 4 Total CO2 and PM2.5 emission, and emissions standardized by burnt area from the Forcett–Dunalley fire
Fig. 3
figure 3

Spatial distribution of CO2 and PM2.5 emissions (in tonnes per 50 m grid cell) from the Forcett–Dunalley fire as a bootstrapped mean of total emissions per grid cell, from the 100 simulations. Note the similarity in emissions patterns for the two emissions

Overall, the dry forest contributed the highest proportion (77–79%) of total CO2 and PM2.5 emissions respectively, while the wet forests contributed approximately 10% of both emissions (Fig. 4). This reflects the greater proportion of dry forests in the area burned at higher intensity, although the highest variance was in the wet forest and Pinus (softwood) plantation (coefficient of variation of ~ 31% CO2 and ~ 40% PM2.5 for both vegetation classes), with only a few areas burning intensely. The emissions variability for the dry forest was around 26–32% for the two pollutants respectively, while the Eucalyptus (hardwood) plantation displayed the lowest variability, at 20–28% for the two pollutants respectively (Fig. 4).

Fig. 4
figure 4

Bootstrapped mean and variability of total emissions from the different vegetation types found within the fireground. a represents CO2 emissions and b PM2.5 emissions

Model comparison

The fine scale estimation, that incorporated detailed fire severity and vegetation mapping, had a better characterisation of the spatial variability of both emission types than GFED (compare Figs. 3 and 5). Nonetheless, GFED detected the area with the highest emissions, with an added advantage of providing a synoptic view of several fires burning across Tasmania. A comparison of total CO2 and PM2.5 emissions between the two inventories revealed comparable emissions estimates, especially for CO2 (Table 4). The fine scale analysis produced total CO2 emissions (and range) of 1.125 Tg (0.893–1.357 Tg) compared to GFED’s estimate of 0.822 Tg which is 73% (range of 65–92%) of the CO2 emissions estimate from the fine scale inventory. However, for PM2.5, GFED reported much lower emissions of 0.006 Tg relative to 0.022 ± 0.006 Tg from the fine scale analysis, representing 30% (24–41%) of the emissions estimate in the fine scale inventory. Per-hectare emissions were comparable but lower for GFED, with 36 t ha−1 for CO2 and 0.3 t ha−1 for PM2.5 (Table 4). It’s worth noting that the area burnt estimate from GFED was approximately 22,851 ha, which is similar to the area estimated by the fine scale analysis (20,200 ha).

Fig. 5
figure 5

Spatial distribution of CO2 and PM2.5 emissions (in tonnes per 28-km grid cell) from several fires in mainland Tasmania, including the Forcett–Dunalley fire (red polygon) for the entire January 2013 from GFED4 analysis

The GFED estimates for the study area were only available until 14 January 2013 and during this period, temporal variability of the two emissions showed similar trends between the fine scale and GFED inventories (Fig. 6). These trends were significantly correlated (r = 0.99, p < 0.05), albeit emission estimates from GFED were always lower than the fine scale analysis. The 4 January had the highest emissions of all days, a day notable for the formation a pyrocumulonimbus (PyroCb). Emissions then drastically declined on 5–6 January and subsequently stabilised at lower values till containment of the fire.

Fig. 6
figure 6

Daily variability of CO2 and PM2.5 emissions from the Forcett–Dunalley fire between the fine scale (FS) and GFED inventories. a, b Represent CO2 variability while c, d show PM2.5 variability for each of the inventories. The error bars represent the standard deviation values around the mean of bootstrapped total daily emissions. 4 January is the day of the PyroCb occurrence

Overall, burnt area mapping from GFED closely aligned with area estimates from the fine scale inventory; total emissions for CO2 were comparable in both inventories; with the models capable of capturing the temporal evolution of CO2 and PM2.5 emissions. However, validation of both inventories using FullCAM simulation over the Forcett–Dunalley fireground yielded approximately 38.6 t ha−1 of carbon emission (or 142 t ha−1 of CO2), which is more than twice the estimates from both inventories.

Fire emissions in Tasmania

Wildfire-derived carbon dioxide emissions and area burnt across Tasmanian fires revealed an interannual variability (Fig. 7a and b), both showing a similar trend where more emissions were produced with an increased area of unplanned fire (correlation of 0.925). Further, correlation for all fires combined (both planned and unplanned) was 0.891 although emissions from planned fires were negatively correlated with area burnt (r = − 0.203), suggesting that increased planned fire area slightly reduces CO2 emissions. Conversely, the wildfire emissions trends do not correlate with Tasmania’s GHG (CO2-equivalent) accounts (Fig. 7c), which show a sharp decline in GHG emissions in 2012 and a stable reduction in the after-years (to being net carbon sink from 2013) despite a spate of large Tasmanian fires in 2013, 2016 and 2019. It is worth noting that fire emissions for the period January-March 2019 are missing from the GFED record, a period characterised by extensive wildfires. It is likely that fire emissions for year 2019 are considerably underestimated.

Fig. 7
figure 7

Time series of carbon emissions across Tasmania for the period 1990–2019. a interannual variability of area burnt within the state; b variability of total annual wildfire emissions based on the available GFED record; and c interannual variability of GHG (CO2-equivalent) emissions according to the State’s Greenhouse Gas Inventory for 2019 that includes the period 1990–2019

From the wildfire-related estimates in Fig. 7b, the Forcett–Dunalley fire represented 28% (almost a third) of fire emissions in Tasmania during the 2013 fires, and 36% and 26% of mean annual fire emissions (3.12 Tg CO2) for the period 1997–2020, based on fine scale and GFED estimates respectively.


This study adopted a ‘bottom-up’ emissions methodology to quantify CO2 and PM2.5 emissions from the 3–18 January 2013 Forcett–Dunalley fire in south-eastern Tasmania. We show that total CO2 and PM2.5 emissions from the fine scale analysis reached 1.125 ± 0.232 Tg and 0.022 ± 0.006 Tg respectively. A comparison of the fine scale (50 m) analysis that uses local fuel and fire severity estimates, and a coarse scale global emissions model GFED (0.25 degrees or ~ 28 km) showed that GFED had a good agreement with the fine-scale analysis regarding total CO2 emissions but not PM2.5 emissions. Naturally, fine scale analysis had more detailed spatial patterns of both emissions than GFED. Validation of the emissions estimates using the FullCAM model yielded 142 t CO2 ha−1 (> 2 times the estimates from both inventories), suggesting that further refinement of FullCAM is important, especially the parameters used in calibrating the model (e.g., debris pool) which are subject to large uncertainties [50].

Other wildfire emissions

A comparison of Forcett–Dunalley fire emissions with other Australian temperate fires showed similarities with some fires and considerable differences with other fires (see Additional file 2: Table S1). For example, the per-hectare CO2 estimate from this study was 55.7 t CO2 ha−1 whereas Volkova et al. [43] reported emission of 105 t CO2 ha−1 from a wildfire in a long-unburnt dry shrubby Eucalyptus forest in Victoria. However, our values are comparable to those reported by these authors from the areas within that wildfire that were previously fuel-reduced (42 t ha−1 of CO2). The 2003 Canberra fire produced 20.2 Tg of CO2 emissions based on the Australian FullCAM model [51], translating to approximately 78 t CO2 ha−1 from the 260,000 ha-fire size, assuming no unburnt patches. However, other studies have reported carbon emissions estimates of 40 M tonnes (or 40 Tg) from the same fire [52]; it is likely that CO2 emissions from that fire exceeded 400 t ha−1 given that CO2 emission are 3.67 times more than carbon emission.

Previous studies in Australia have shown high agreement between GFED and other models/field observations in CO2 emissions e.g., Paton-Walsh et al. [24]. This is despite GFED treating vegetation types, particularly Eucalyptus forests and woodlands, and fire behaviour in south-eastern Australia as the same as those found in the temperate biomes in Northern Hemisphere. The overall good performance of GFED’s CO2 estimates in this study also likely reflects an improved detection of smaller fires in GFED4 compared to previous versions of GFED [48].

Per-hectare estimates for PM2.5 in this study (1.1 t ha−1) were inconsistent with emissions estimates from other Australian temperate fires (Additional file 2: Table S1). For example, Reisen et al. [26] reported emissions of 73.7–163.9 kg ha−1 (0.07–0.16 t ha−1) from prescribed fires in Victorian Eucalyptus forests while another Tasmanian study reported PM2.5 emissions of 7789 tonnes (or 6.9 t ha−1) from a high-intensity regeneration fire in a southern Tasmanian native forest [53]. It should be noted that there is paucity of data on PM2.5 emission from temperate Australian forest fires; most of the studies have instead focused on PM2.5 concentration in urban airsheds for air quality purposes, involving a mix of emission sources. Beyond Australia, western US wildfires between 2011 and 2015 were estimated to have emitted 1530 Gg (1.53 Tg) of PM2.5 annually [54]. Similar to our study, the authors report that the emissions were three times higher than the estimates from the US national inventory. Further, in another study, the GFED3 PM2.5 emission estimate across contiguous US was lower by a factor of eight compared to the national emissions inventory [55], revealing a likely systematic underestimation of PM emission across jurisdictions.

Deficiencies in current fire emissions approaches

The discrepancy in GFED modelling in this study was the lower PM2.5 emissions by a factor of three, likely due to lower emissions factors (EFs) used for PM2.5 within GFED (12.9 g kg−1). These EFs do not accurately reflect temperate Eucalyptus-dominated fuels in Australia, as they are averaged across the temperate biome globally. One of the main differences significantly affecting emissions amongst the temperate biomes is fire behaviour. For example, compared to other biomes, Australian forests and woodlands typically have a higher biomass of sclerophyllous leaves and bark, which burn intensely and support short-long distance transport and spotting of embers that spread landscape fire [56]. Eucalyptus fuels have lower rates of decomposition (and therefore low/absent duff layer [57]) compared to northern hemisphere conifer/boreal forests that have a more-developed duff layer that supports smouldering combustion and can contribute up to 50–74% of fuel consumption [58]. An upward revision of PM2.5 EFs to 16.9–38.8 g kg−1 [26] is therefore recommended to better accommodate typical fuels within these Australian ecosystems.

The accuracy of bottom-up approaches (such as the above inventories) that adopt fuel consumption estimates in emissions estimations has been a topic of debate relative to the more accurate top-down approaches that use satellite observations to directly estimate emissions within the atmospheric column [59,60,61]. Despite these limitations, two previous carbon emissions studies on the recent Australian Black Summer fires using top-down and bottom-up approaches revealed comparable CO2 estimates between the two methods [3, 13]. This highlights the importance of validating emissions estimations with diverse methods, including satellite and on-ground observations, to reduce the inherent uncertainties.

Smoke emissions analyses are constrained by the quality and representativeness of data on fuel types, requiring greater sampling of a broader range of vegetation and fuels [20]. Field protocols should include detailed inventories of vegetation characteristics, e.g., Prior et al. [62] and measurement of fuel loads across all fuel components, ranging from subsurface to overstorey fuels, and from fine to woody fuels. To date, coarse woody debris (CWD) estimation, being the less studied fuel component than fine fuels, is the most common source of emissions uncertainties in temperate Australian landscapes. This is because CWD is influenced in different regions by among other factors, the disturbance history (past fire or logging activities), forest age, and site productivity [18, 63]. More field inventories across Australia and particularly in Tasmania where there has been scarcity of fuel load data [45] are needed to provide confidence in emissions estimates.

Fire behaviour modelling in Australia has shifted from an emphasis on fine fuel loads, to a more realistic determination of fuel hazard scores across fuel types; nonetheless, we contend that there remains a need for accurate fine and coarse fuel load measurements to underpin fire emissions analysis [64]. These inventories could make use of recent technologies such as LiDAR to increase the accuracy of fuel estimation, especially the amount of coarse woody debris, within a forest. Previous research has shown that carbon losses from forest regeneration burns are around 200 t ha−1 [65]. However, the relationship between forest harvesting and likelihood of uncontrolled fires, that would cause higher carbon emissions than if native forests were unharvested, is highly controversial and demands further research [66, 67]. Another important knowledge gap concerns the comparative assessment of particulate and carbon emissions and associated costs of fuel management burns, post-logging (or regeneration) burns and wildfires. Previous research into health economics suggests the public health cost of both fuel management burns and wildfires can be substantial [68].

Fire severity scales with fuel consumption, with high-severity fires typically associated with high consumption of vegetation; however, the general lack of empirical fuel consumption data can introduce variability in total emissions, despite the availability of fire severity information. This was evident in the spectral signatures (from satellite observations) in grassland areas of the Forcett–Dunalley fireground which exhibited very high severities despite their very low fuel loads and minimal biological impact. Fuel consumption estimates in this study were inferred from a few studies on temperate Eucalyptus forests (Table 2). Therefore, there is need to improve data collection of fuel consumption during wildland fires (supplemented by remote sensing), and measurement of residence time of flaming and smouldering to partition emissions into the different combustion stages. Although these attributes can be inferred from lab experiments, variability in fuel size, especially coarser fuels are difficult to accurately characterise in the lab [69]. There is also a need to clearly establish a quantitative link between severity measurements and fuel consumption for better applicability of fire severity data in future emissions studies.

Greenhouse gas accounting

Estimates of emissions from wildfires are of increasing interest given their contribution to climate change. Indeed, emissions from Australian wildfires are accounted for in the national GHG accounting to the Intergovernmental Panel on Climate Change, however, what constitutes a wildfire and a human-caused fire in the accounting is subject to debate and a number of pragmatic and often poorly justified ‘rules’. For example, the Australian Government accounting uses a burned area threshold (that is 16,950 ha in Tasmania) and fire emissions threshold (2 standard deviations above the mean of gross annual fire emissions) to exclude large fires or fire years, with the assumption that the fires were not human-caused and therefore are under no human control [29]. These statistically large fires are therefore attributed as natural disturbances and are excluded in the final carbon accounting, in the same way post-logging regeneration fires are excluded. It is therefore likely that the Forcett–Dunalley fire (with a burnt area of > 20,000 ha was excluded based on these criteria despite it being anthropogenically-caused. While there is some logic to this reasoning, there is uncertainty as to how to treat severe wildfires, such as the Dunalley disaster, that are human-caused, are exacerbated by anthropogenic climate change, burn over a highly human-modified landscape, and are subject to intensive human control efforts, yet they exceed the above threshold for defining anthropogenic fires.

Although, it is commendable that from the year 2019, the Australian government can report to IPCC on fire emissions within the ‘natural disturbance’ provision [51], we recommend inclusion of all emissions from large, human-caused fires as well as post-logging burns at state and national levels in the final accounting, to prevent situations where net carbon credits are claimed despite insufficient fire management. Current accounting approaches can potentially lead to perverse outcomes where carbon neutrality could be claimed by reducing the extent of planned fires that are an important tool in mitigating uncontrolled bushfire and reducing emissions (Fig. 7). Current arrangements therefore provide disincentives to effective wildfire management to reduce carbon emissions from large wildfires and post-logging fires that ultimately exacerbate climate change. Furthermore, the national policy is inconsistent because in north Australian savannas, there are carbon emissions abatement programs which reward pre-emptive early dry season burning to limit the high smoke emissions associated with late season burning [70].

Tasmanian government’s GHG reporting reveals that since 2012, forestry-related activities (LULUCF) have counteracted anthropogenic non-forestry GHG emissions [71, 72], with an average removal of − 9.17 Tg between the years 2012–2019, and an increased carbon sequestration from − 5.920 Tg in 2012 to − 10.04 Tg in 2019. These estimates seem impressive; however, they are unaffected by major wildfires such as Dunalley disaster that according to the GFED model, accounted for one third of the state’s annual fire emissions. If severe-fire emissions were incorporated in the forestry-related GHG accounting for 2013 (− 10.952 Tg in forest land), Dunalley CO2 emissions (1.125 Tg) could have reduced forest land CO2 sequestration (or removal) by 10%. These results suggest that if wildfire emissions are included, then Tasmania may not be actually achieving carbon neutrality.

An important consideration in the understanding and accounting of carbon emissions is the influence of climate change on, and feedbacks with, fire regimes. In the GHG accounting across many national jurisdictions, the emitted carbon from wildfires is assumed to be assimilated by forests in the following growing seasons via tree growth, and therefore carbon uptake post-fire can be substantial. However, it is not clear how the regrowth and carbon sequestration can be relied upon in a changing hotter or drier climate. For instance, a warming earth has increased the vulnerability of ecosystems to frequent and intense fires, which in turn emit large quantities of emissions, thereby creating a positive feedback loop where forests are converted to a treeless state [73]. This calls for more investigation using diverse tools ranging from experiments, observations and models, to understand the complex interactions between climate, ecosystem structure and fire dynamics.


This study quantified CO2 and PM2.5 emissions from the January 2013 Forcett–Dunalley fire using two standard emissions inventories. We report the release of approximately 1.125 ± 0.232 Tg of CO2 and 0.022 ± 0.006 Tg of PM2.5 into the atmosphere using a basic model that incorporated local fuel attributes. We investigated the reliability of a global model GFED4 in emissions estimation assuming the absence of field data. Our findings show that both the fine scale and GFED inventories produced comparable estimates for CO2, although PM2.5 estimates were lower by a factor of three for GFED. We therefore show that GFED was able to produce reliable emissions estimates within the limits of emissions uncertainties, although the model did not accurately capture the spatial distribution of the two emissions. By contextualising these estimates with wildfire emissions and overall GHG accounting in Tasmania, we show that the fire injected approximately 30% of fire emissions during the 2013 fire season, and represented 25–34% of mean annual fire emissions from the state. These findings showed the influence of the extreme fire event to overall carbon balance for the state, although the Forcett–Dunalley fire appears to have been excluded from the state and national carbon accounting due to the criteria that excludes natural disturbances fires. Such exclusions could have a major influence on a national or local jurisdiction’s claim of carbon neutrality. This analysis also investigated knowledge gaps in emissions quantification in Australian temperate Eucalyptus forests. We show that fuel attributes, especially the amount of coarse wood fuels within a forest stand, and the fraction of fuel consumed, contributed the most to uncertainties in emissions estimates. More accurate fine-scale analyses demand improved data on fuel types and their emission factors.

Availability of data and materials

The datasets supporting the conclusions of this article are available in the following open-source databases: Burned area (fire history) records were obtained from Tasmania’s LISTmap (, last accessed: 4 March 2022). Two gridded datasets used in GFED4s emissions estimation (combusted dry matter and burned area) were downloaded from GFED website (, last accessed: 4 March 2022). Time series data on Tasmania’s carbon accounting was extracted from the State and Territory Greenhouse Gas Inventory 2019, covering the period 1990–2019. (, last accessed: 4 March 2022).


  1. Jolly WM, Cochrane MA, Freeborn PH, Holden ZA, Brown TJ, Williamson GJ, Bowman DMJS. Climate-induced variations in global wildfire danger from 1979 to 2013. Nat Commun. 2015;6:7537.

    Article  CAS  Google Scholar 

  2. Urbanski S. Wildland fire emissions, carbon, and climate: emission factors. For Ecol Manage. 2014;317:51–60.

    Article  Google Scholar 

  3. Bowman DMJS, Williamson GJ, Price OF, Ndalila MN, Bradstock RA. Australian forests, megafires and the risk of dwindling carbon stocks. Plant Cell Environ. 2021;44(2):347–55.

    Article  CAS  Google Scholar 

  4. Reid JS, Koppmann R, Eck TF, Eleuterio DP. A review of biomass burning emissions part II: intensive physical properties of biomass burning particles. Atmos Chem Phys. 2005;5(3):799–825.

    Article  CAS  Google Scholar 

  5. Bowman DMJS, Balch JK, Artaxo P, Bond WJ, Carlson JM, Cochrane MA, et al. Fire in the earth system. Science. 2009;324(5926):481–4.

    Article  CAS  Google Scholar 

  6. McCormick MP, Thomason LW, Trepte CR. Atmospheric effects of the Mt Pinatubo eruption. Nature. 1995;373:399–404.

    Article  CAS  Google Scholar 

  7. Bond TC, Doherty SJ, Fahey DW, Forster PM, Berntsen T, DeAngelo BJ, Flanner MG, Ghan S, Kärcher B, Koch D, et al. Bounding the role of black carbon in the climate system: a scientific assessment. J Geophys Res Atmos. 2013;118(11):5380–552.

    Article  CAS  Google Scholar 

  8. Johnston FH, Hanigan IC, Henderson SB, Morgan GG, Portner T, Williamson GJ, Bowman DMJS. Creating an integrated historical record of extreme particulate air pollution events in Australian cities from 1994 to 2007. J Air Waste Manage Assoc. 2011;61(4):390–8.

    Article  CAS  Google Scholar 

  9. Di Virgilio G, Evans JP, Blake SAP, Armstrong M, Dowdy AJ, Sharples J, McRae R. Climate change increases the potential for extreme wildfires. Geophys Res Lett. 2019;46(14):8517–26.

    Article  Google Scholar 

  10. Canadell JG, Meyer CP, Cook GD, Dowdy A, Briggs PR, Knauer J, Pepler A, Haverd V. Multi-decadal increase of forest burned area in Australia is linked to climate change. Nat Commun. 2021;12(1):6921.

    Article  CAS  Google Scholar 

  11. van Oldenborgh GJ, Krikken F, Lewis S, Leach NJ, Lehner F, Saunders KR, van Weele M, Haustein K, Li S, Wallom D, et al. Attribution of the Australian bushfire risk to anthropogenic climate change. Nat Hazards Earth Syst Sci. 2021;21(3):941–60.

    Article  Google Scholar 

  12. Nolan RH, Bowman DMJS, Clarke H, Haynes K, Ooi MKJ, Price OF, Williamson GJ, Whittaker J, Bedward M, Boer MM, et al. What do the Australian black summer fires signify for the global fire crisis? Fire. 2021;4(4):97.

    Article  Google Scholar 

  13. van der Velde IR, van der Werf GR, Houweling S, Maasakkers JD, Borsdorff T, Landgraf J, Tol P, van Kempen TA, van Hees R, Hoogeveen R, et al. Vast CO2 release from Australian fires in 2019–2020 constrained by satellite. Nature. 2021;597(7876):366–9.

    Article  CAS  Google Scholar 

  14. Peterson DA, Fromm MD, McRae RHD, Campbell JR, Hyer EJ, Taha G, Camacho CP, Kablick GP, Schmidt CC, DeLand MT. Australia’s black summer pyrocumulonimbus super outbreak reveals potential for increasingly extreme stratospheric smoke events. NPJ Clim Atmos Sci. 2021;4(1):38.

    Article  Google Scholar 

  15. Johnston FH, Borchers-Arriagada N, Morgan GG, Jalaludin B, Palmer AJ, Williamson GJ, Bowman DMJS. Unprecedented health costs of smoke-related PM2.5 from the 2019–20 Australian megafires. Nat Sustain. 2021;4(1):42–7.

    Article  Google Scholar 

  16. Mallapaty S. Australian bush fires belched out immense quantity of carbon. Nature. 2021;597(7877):459–60.

    Article  CAS  Google Scholar 

  17. Hurst DF, Griffith DWT. Trace gas emissions from biomass burning in tropical Australian savannas. J Geophys Res Atmos. 1994;99(D8):16441–56.

    Article  CAS  Google Scholar 

  18. Russell-Smith J, Murphy BP, Meyer CP, Cook GD, Maier S, Edwards AC, Schatz J, Brocklehurst P. Improving estimates of savanna burning emissions for greenhouse accounting in northern Australia: limitations, challenges, applications. Int J Wildland Fire. 2009;18(1):1–18.

    Article  CAS  Google Scholar 

  19. Meyer CP, Cook GD, Reisen F, Smith TEL, Tattaris M, Russell-Smith J, Maier SW, Yates CP, Wooster MJ. Direct measurements of the seasonality of emission factors from savanna fires in northern Australia. J Geophys Res Atmos. 2012.

    Article  Google Scholar 

  20. Volkova L, Weston C. Redistribution and emission of forest carbon by planned burning in Eucalyptus obliqua (L. Hérit.) forest of south-eastern Australia. For Ecol Manage. 2013;304:383–90.

    Article  Google Scholar 

  21. Paton-Walsh C, Smith TEL, Young EL, Griffith DWT, Guérette ÉA. New emission factors for Australian vegetation fires measured using open-path Fourier transform infrared spectroscopy—part 1: methods and Australian temperate forest fires. Atmos Chem Phys. 2014;14(20):11313–33.

    Article  CAS  Google Scholar 

  22. Guérette EA, Paton-Walsh C, Desservettaz M, Smith TEL, Volkova L, Weston CJ, Meyer CP. Emissions of trace gases from Australian temperate forest fires: emission factors and dependence on modified combustion efficiency. Atmos Chem Phys. 2018;18(5):3717–35.

    Article  CAS  Google Scholar 

  23. Surawski NC, Sullivan AL, Meyer CP, Roxburgh SH, Polglase PJ. Greenhouse gas emissions from laboratory-scale fires in wildland fuels depend on fire spread mode and phase of combustion. Atmos Chem Phys. 2015;15(9):5259–73.

    Article  CAS  Google Scholar 

  24. Paton-Walsh C, Emmons LK, Wilson SR. Estimated total emissions of trace gases from the Canberra Wildfires of 2003: a new method using satellite measurements of aerosol optical depth & the MOZART chemical transport model. Atmos Chem Phys. 2010;10(12):5739–48.

    Article  CAS  Google Scholar 

  25. Lawson SJ, Cope M, Lee S, Galbally IE, Ristovski Z, Keywood MD. Biomass burning at Cape Grim: exploring photochemistry using multi-scale modelling. Atmos Chem Phys. 2017;17(19):11707–26.

    Article  CAS  Google Scholar 

  26. Reisen F, Meyer CP, Weston CJ, Volkova L. Ground-based field measurements of PM2.5 emission factors from flaming and smoldering combustion in eucalypt forests. J Geophys Res Atmos. 2018;123(15):8301–14.

    CAS  Google Scholar 

  27. Johnston FH, Henderson SB, Chen Y, Randerson JT, Marlier M, DeFries RS, Kinney P, Bowman DMJS, Brauer M. Estimated global mortality attributable to smoke from landscape fires. Environ Health Perspect. 2012;120(5):695–701.

    Article  Google Scholar 

  28. Marlier ME, DeFries RS, Voulgarakis A, Kinney PL, Randerson JT, Shindell DT, Chen Y, Faluvegi G. El Niño and health risks from landscape fire emissions in southeast Asia. Nat Clim Chang. 2013;3(2):131–6.

    Article  CAS  Google Scholar 

  29. Commonwealth of Australia: national inventory report volume 2. Australian Government Department of Industry, Science, Energy and Resources. 2020. Accessed 12 July 2021.

  30. Williamson GJ, Bowman DMJS, Price OF, Henderson SB, Johnston FH. A transdisciplinary approach to understanding the health effects of wildfire and prescribed fire smoke regimes. Environ Res Lett. 2016;11:125009.

    Article  CAS  Google Scholar 

  31. Langmann B, Duncan B, Textor C, Trentmann J, van der Werf GR. Vegetation fire emissions and their impact on air pollution and climate. Atmos Environ. 2009;43(1):107–16.

    Article  CAS  Google Scholar 

  32. Ndalila MN, Williamson GJ, Fox-Hughes P, Sharples J, Bowman DMJS. Evolution of a pyrocumulonimbus event associated with an extreme wildfire in Tasmania, Australia. Nat Hazards Earth Syst Sci. 2020;20(5):1497–511.

    Article  Google Scholar 

  33. Ndalila MN, Williamson GJ, Bowman DMJS. Geographic patterns of fire severity following an extreme eucalyptus forest fire in southern Australia: 2013 Forcett–Dunalley fire. Fire. 2018;1(3):40.

    Article  Google Scholar 

  34. Seiler W, Crutzen PJ. Estimates of gross and net fluxes of carbon between the biosphere and the atmosphere from biomass burning. Clim Chang. 1980;2(3):207–47.

    Article  CAS  Google Scholar 

  35. Hijmans RJ, Cameron SE, Parra JL, Jones PG, Jarvis A. Very high resolution interpolated climate surfaces for global land areas. Int J Climatol. 2005;25(15):1965–78.

    Article  Google Scholar 

  36. Reid CE, Brauer M, Johnston FH, Jerrett M, Balmes JR, Elliott CT. Critical review of health impacts of wildfire smoke exposure. Environ Health Perspect. 2016;124(9):1334–43.

    Article  Google Scholar 

  37. Murphy BP, Prior LD, Cochrane MA, Williamson GJ, Bowman DMJS. Biomass consumption by surface fires across Earth’s most fire prone continent. Glob Chang Biol. 2019;25(1):254–68.

    Article  Google Scholar 

  38. Hollis JJ, Anderson WR, McCaw WL, Cruz MG, Burrows ND, Ward B, Tolhurst KG, Gould JS. The effect of fireline intensity on woody fuel consumption in southern Australian eucalypt forest fires. Aust For. 2011;74(2):81–96.

    Article  Google Scholar 

  39. Knorr W, Lehsten V, Arneth A. Determinants and predictability of global wildfire emissions. Atmos Chem Phys. 2012;12(15):6845–61.

    Article  CAS  Google Scholar 

  40. Weise DR, Wright CS. Wildland fire emissions, carbon and climate: characterizing wildland fuels. For Ecol Manage. 2014;317:26–40.

    Article  Google Scholar 

  41. Hollis JJ, Matthews S, Ottmar RD, Prichard SJ, Slijepcevic A, Burrows ND, Ward B, Tolhurst KG, Anderson WR, Gould JS. Testing woody fuel consumption models for application in Australian southern eucalypt forest fires. For Ecol Manage. 2010;260(6):948–64.

    Article  Google Scholar 

  42. O'Loughlin EM, Cheney NP, Burns J. The Bushrangers experiment: hydrological response of a eucalypt catchment to fire. In: National Symposium on Forest Hydrology (The Institution of Engineers, Australia: Canberra). 1982;132–138.

  43. Volkova L, Meyer CM, Murphy S, Fairman T, Reisen F, Weston C. Fuel reduction burning mitigates wildfire effects on forest carbon and greenhouse gas emission. Int J Wildland Fire. 2014;23(6):771–80.

    Article  CAS  Google Scholar 

  44. Environment Australia. NPI emissions estimation technique manual for aggregated emissions from prescribed burning and wildfires. Environment Australia. 1999.

  45. Roxburgh S, Volkova L, Surawski N, Meyer M, Weston C. Review of fuel loads, burn efficiencies, emissions factors and recovery functions used to estimate greenhouse gas emissions and removals associated with wildfire on temperate forested lands. CSIRO. 2015.

    Article  Google Scholar 

  46. R Core Team. R: a language and environment for statistical computing. R foundation for statistical computing, Vienna, Austria. 2019. Accessed 4 March 2022.

  47. Environmental Systems Research Institute (ESRI). ArcGIS desktop release 10.3. Redlands, California. 2015.

  48. van der Werf GR, Randerson JT, Giglio L, van Leeuwen TT, Chen Y, Rogers BM, Mu MQ, van Marle MJE, Morton DC, Collatz GJ, et al. Global fire emissions estimates during 1997–2016. Earth Syst Sci Data. 2017;9(2):697–720.

    Article  Google Scholar 

  49. Global Fire Emissions Database Version 4.1 (GFEDv4). Accessed 4 Mar 2022.

  50. Surawski NC, Sullivan AL, Roxburgh SH, Cook GD. Review of FullCAM forest fire event parameters with recommendations supported by a literature review. CSIRO client report EP 28061232 submitted to the Department of Climate Change and Energy Efficiency. 2012.

  51. Commonwealth of Australia. Estimating greenhouse gas emissions from bushfires in Australia’s temperate forests: focus on 2019–20: Department of Industry, Science, Energy and Resources. 2020.

  52. Mitchell RM, O’Brien DM, Campbell SK. Characteristics and radiative impact of the aerosol generated by the Canberra firestorm of January 2003. J Geophys Res Atmos. 2006;111(D2):D02204.

    Article  Google Scholar 

  53. Meyer CP, Reisen F, Keywood MD, Crumeyrolle S. Impacts of smoke from regeneration burning on air quality in the Huon Valley, Tasmania. Melbourne: CSIRO; 2011.

    Google Scholar 

  54. Liu X, Huey LG, Yokelson RJ, Selimovic V, Simpson IJ, Müller M, Jimenez JL, Campuzano-Jost P, Beyersdorf AJ, Blake DR, et al. Airborne measurements of western U.S. wildfire emissions: comparison with prescribed burning and air quality implications. J Geophys Res Atmos. 2017;122(11):6108–29.

    Article  CAS  Google Scholar 

  55. Larkin NK, Raffuse SM, Strand TM. Wildland fire emissions, carbon, and climate: U.S. emissions inventories. For Ecol Manage. 2014;317:61–9.

    Article  Google Scholar 

  56. Cruz M, Gould J, Hollis J, McCaw W. A hierarchical classification of wildland fire fuels for Australian vegetation types. Fire. 2018;1(1):13.

    Article  Google Scholar 

  57. Gould J, Cruz M. Australian fuel classification: stage II. Ecosystem sciences and climate adaption flagship. Canberra: CSIRO; 2012.

    Google Scholar 

  58. Bertschi I, Yokelson RJ, Ward DE, Babbitt RE, Susott RA, Goode JG, Hao WM. Trace gas and particle emissions from fires in large diameter and belowground biomass fuels. J Geophys Res Atmos. 2003.

    Article  Google Scholar 

  59. Wang J, Christopher SA, Nair US, Reid JS, Prins EM, Szykman J, Hand JL. Mesoscale modeling of Central American smoke transport to the United States: 1. “Top-down” assessment of emission strength and diurnal variation impacts. J Geophys Res Atmos. 2006.

    Article  Google Scholar 

  60. Kaiser JW, Heil A, Andreae MO, Benedetti A, Chubarova N, Jones L, Morcrette JJ, Razinger M, Schultz MG, Suttie M, et al. Biomass burning emissions estimated with a global fire assimilation system based on observed fire radiative power. Biogeosciences. 2012;9(1):527–54.

    Article  CAS  Google Scholar 

  61. Ichoku C, Ellison L. Global top-down smoke-aerosol emissions estimation using satellite fire radiative power measurements. Atmos Chem Phys. 2014;14(13):6643–67.

    Article  CAS  Google Scholar 

  62. Prior LD, Murphy BP, Williamson GJ, Cochrane MA, Jolly WM, Bowman DMJS. Does inherent flammability of grass and litter fuels contribute to continental patterns of landscape fire activity? J Biogeogr. 2017;44(6):1225–38.

    Article  Google Scholar 

  63. Volkova L, Bi H, Murphy S, Weston CJ. Empirical estimates of aboveground carbon in open eucalyptus forests of south-eastern Australia and its potential implication for national carbon accounting. Forests. 2015;6(10):3395–411.

    Article  Google Scholar 

  64. Volkova L, Roxburgh SH, Surawski NC, Meyer CP, Weston CJ. Improving reporting of national greenhouse gas emissions from forest fires for emission reduction benefits: an example from Australia. Environ Sci Policy. 2019;94:49–62.

    Article  Google Scholar 

  65. Slijepcevic A. Loss of carbon during controlled regeneration burns in Eucalyptus obliqua forest. Tasforests. 2001;13:281–90.

    Google Scholar 

  66. Bowman DMJS, Williamson GJ, Gibson RK, Bradstock RA, Keenan RJ. Reply to: logging elevated the probability of high-severity fire in the 2019–20 Australian forest fires. Nat Ecol Evol. 2022.

    Article  Google Scholar 

  67. Lindenmayer DB, Zylstra P, Kooyman R, Taylor C, Ward M, Watson JEM. Logging elevated the probability of high-severity fire in the 2019–20 Australian forest fires. Nat Ecol Evol. 2022.

    Article  Google Scholar 

  68. Borchers-Arriagada N, Bowman DMJS, Price O, Palmer AJ, Samson S, Clarke H, Sepulveda G, Johnston FH. Smoke health costs and the calculus for wildfires fuel management: a modelling study. Lancet Planet Health. 2021;5(9):e608–19.

    Article  Google Scholar 

  69. Urbanski SP. Combustion efficiency and emission factors for wildfire-season fires in mixed conifer forests of the northern Rocky Mountains, US. Atmos Chem Phys. 2013;13(14):7241–62.

    Article  CAS  Google Scholar 

  70. Edwards A, Archer R, De Bruyn P, Evans J, Lewis B, Vigilante T, Whyte S, Russell-Smith J. Transforming fire management in northern Australia through successful implementation of savanna burning emissions reductions projects. J Environ Manage. 2021;290: 112568.

    Article  Google Scholar 

  71. Commonwealth of Australia. National greenhouse accounts 2019. State and territory greenhouse gas inventories: annual emissions. Australian Government Department of Industry, Science, Energy and Resources. 2021. Accessed 22 Nov 2021.

  72. Tasmanian Climate Change Office. Tasmanian greenhouse gas emissions report 2021. Department of Premier and Cabinet. 2021. Accessed 03 May 2022.

  73. Bowman DMJS, Kolden CA, Abatzoglou JT, Johnston FH, van der Werf GR, Flannigan M. Vegetation fires in the Anthropocene. Nat Rev Earth Environ. 2020;1(10):500–15.

    Article  Google Scholar 

  74. Bresnehan SJ. An assessment of fuel characteristics and fuel loads in the dry sclerophyll forests of south-east Tasmania. Hobart: University of Tasmania; 2003.

    Google Scholar 

  75. Watson P. Fuel load dynamics in NSW vegetation. Part 1: forests and grassy woodlands. Centre for Environmental Risk Management of Bushfires, University of Wollongong. 2011.

  76. Hollis JJ, Matthews S, Anderson WR, Cruz MG, Burrows ND. Behind the flaming zone: predicting woody fuel consumption in eucalypt forest fires in southern Australia. For Ecol Manage. 2011;261(11):2049–67.

    Article  Google Scholar 

  77. Woldendorp G, Keenan RJ. Coarse woody debris in Australian forest ecosystems: a review. Austral Ecol. 2005;30(8):834–43.

    Article  Google Scholar 

  78. Bowman DMJS, Murphy BP, Neyland DLJ, Williamson GJ, Prior LD. Abrupt fire regime change may cause landscape-wide loss of mature obligate seeder forests. Glob Chang Biol. 2014;20(3):1008–15.

    Article  Google Scholar 

  79. Volkova L, Bi HQ, Hilton J, Weston CJ. Impact of mechanical thinning on forest carbon, fuel hazard and simulated fire behaviour in Eucalyptus delegatensis forest of south-eastern Australia. For Ecol Manage. 2017;405:92–100.

    Article  Google Scholar 

  80. Newnham G, Opie K, Leonard J. A methodology for State-wide mapping annual fuel load and bushfire hazard in Queensland. CSIRO: EP175130. CSIRO. 2017.

  81. Leonard S. Predicting sustained fire spread in Tasmanian native grasslands. Environ Manage. 2009;44(3):430.

    Article  Google Scholar 

Download references


We are grateful to Lynda Prior for the provision of some datasets for the study, and the Tasmania Fire Service for their end-user support through BNHCRC funding. We are also thankful to John Hunter for the discussions on carbon accounting protocols within Australia. The support of the University of Tasmania to the main author is also acknowledged.


This research was funded by the Australian Research Council Linkage Program (Grant No. LLP130100146) and the Bushfire and Natural Hazards Cooperative Research Centre (BNHCRC). The funding agencies had no role in the design of the study and collection, analysis, and interpretation of data and in writing the manuscript.

Author information

Authors and Affiliations



MNN, GJW, and DMJSB designed the research; MNN analysed the data, with some additional analysis provided by GJW who also developed the code for emissions analysis. MNN prepared the paper with contributions from DMJSB and GJW. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Mercy N. Ndalila.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1: Figure S1.

Spatiotemporal progression of combustion. Spatiotemporal progression of combustion during the early days of the fire, from classification of infrared linescan imagery obtained from a Victoria DELWP aircraft. The 4 January displayed dynamic fire behaviour of all the days during the fire. The original 20-cm resolution imagery has been resampled after classification to fit the 50-m resolution of the analysis.

Additional file 2: Table S1.

Comparison of total emissions (in Tg) and per-hectare emissions (in t ha−1) among wildfires in Australia. Comparison of total emissions (in Tg) and per-hectare emissions (in t ha−1) among wildfires in Australia. Burnt area estimates (BA; in ha) for each fire event are indicated in brackets. CO2-equivalent (CO2-e) emissions are totals from CO2, methane and nitrous oxide emissions. The estimate for the Forcett–Dunalley fire (this study) has also been compared with estimates from the FullCAM model that is used in Australia for national GHG accounting.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Ndalila, M.N., Williamson, G.J. & Bowman, D.M.J.S. Carbon dioxide and particulate emissions from the 2013 Tasmanian firestorm: implications for Australian carbon accounting. Carbon Balance Manage 17, 7 (2022).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: