GFED4 biomass burning emissions

From Geos-chem
Jump to: navigation, search

GFED4 biomass burning emissions are now included in GEOS-Chem v10-01 and higher versions.

Source code and data

GFED4 was implemented into GEOS-Chem v10-01 as an extension of the HEMCO emissions component. The source code files may be found in the GEOS-Chem source code folders:

File Description
HEMCO/Extensions/hcox_gfed_mod.F90 Module for computing GFED4 biomass emissions
HEMCO/Extensions/hcox_gfed_include_gfed4.H Emission factors for GFED4 species (inlined for computational efficiency)

We have created new GFED4 data files (in COARDS-compliant netCDF format) for use with HEMCO. These new data files are contained in the HEMCO data directory tree. For detailed instructions on how to download these data files to your disk server, please see our Downloading the HEMCO data directories wiki post.

--Bob Y. (talk) 15:53, 17 June 2015 (UTC)

Data visualization

I uploaded GFEDv4s data (1997-2018) to Google Earth Engine and created an online app to visualize GFEDv4s emissions. The GFEDv4s app plots annual and monthly time series of burned area or a given species for custom regions: (1) global, (2) basis region, (3) country, (4) pixel, and (5) custom polygon using lon/lat coordinates. Emissions data are also further subdivided by the six GFEDv4s land use and land cover types: SAVA, DEFO, BORF, TEMF, PEAT, and AGRI.

To compare GFEDv4s with other global fire emissions inventories (GFASv1.2, FINNv1.5, QFEDv2.5r1, and FEERv1.0-G1.2), please use the FIRECAM tool. Currently, FIRECAM supports six species: CO2, CO, CH4, OC, BC, and PM2.5. FIRECAM compares the five inventories at an aggregated 0.5deg x 0.5deg spatial resolution for the 2003-2018 time period. Please see the FIRECAM website for more information.

--Tianjia Liu (talk) 14:53, 21 August 2019 (EDT)

Discussion and updates

Update to GFED 4.1

This update was validated with 1-month benchmark simulation v11-01d and 1-year benchmark simulation v11-01d-Run1. This version was approved on 12 Dec 2015.

Prasad Kasibhatla wrote:

Please see attached. This latest version GFED 4.1s is being released July 3, 2015.
The monthly GFED files that go into GEOS-Chem have to be updated (though updates are very minor). In addition, from 2003 onwards there is now daily and diurnal info available in the dataset - this should be added to GEOS-Chem.

--Melissa Sulprizio (talk) 13:25, 13 July 2015 (UTC)

Discussion following 1-month benchmark v10-01h

GFED4 biomass burning emissions were validated in GEOS-Chem v10-01h. The following discussion ensued.

Prasad Kasibhatla wrote:

I have been looking at the CO biomass emission ratio plot in the 10-01h benchmark - there are a couple of regions where things look a bit different than what I expect from the emissions I compute offline at native resolution. Would it be possible for you to send me plots of absolute CO emissions from biomass burning from v10-01h and v10-01f_UCX that went into the ratio plot. And a difference plot?

Melissa Sulprizio wrote:

Please see the attached image for maps of CO biomass emissions from v10-01f_UCX (upper left) and v10-01h (upper right). The absolute difference (lower left) and percent difference (lower right) maps are also included.
V10-01 CO biomass emission diffs.gif

Prasad Kasibhatla wrote:

I computed CO emissions offline for GFED3 for July 2011 and GFED4 for July 2013 (corresponding to what was used in the v10-01f_UCX and v10-01h benchmark runs) and reproduced the emissions plot you sent.
Don't worry about the labels on my plot - my panels correspond to your panels. As you can see your patterns match patterns I compute - but absolute values of emissions that you are getting from the model are higher than what I calculate offline. I have not been able to figure out if I am making a mistake in my offline calculation or if there is a problem in the geoschem implementation.
Obviously, we need to figure this out. Perhaps Christoph and I can focus on one or two pixels and resolve what is going on?
GFED4 CO emissions from prasad.png
For both the GFED3 and GFED4 offline plots, I accumulated dry matter burned in various categories into the generic 1x1 grid, computed emissions in molecules/cm2 on the 1x1 grid using category-specific emission factors, and then regridded to the 4x5 grid using ctm_regridh with the /per_unit_area option.

Christoph Keller wrote:

Thanks, this may explain it. In v10-01e and newer, dry matter burned and vegetation mask are first regridded onto the simulation grid, and emissions are computed on that grid.
I looked into this for version v10-01e, and the changes in GFED-3 emissions (compared to pre-HEMCO versions) are attributable to the regridding order. In v10-01h, we didn’t change anything in GFED-3, and GFED-4 uses the same code as GFED-3 except for the updated dry matter emissions fields.

--Melissa Sulprizio 13:38, 27 March 2015 (EDT)


  1. GFED website:
  2. Akagi, S. K., R. J. Yokelson, C. Wiedinmyer, M. J. Alvarado, J. S. Reid, T. Karl, J. D. Crounse, and P. O. Wennberg, Emission factors for open and domestic biomass burning for use in atmospheric models, Atmos. Chem. Phys., 11, 4039-4072, 2011. Article
  3. Andreae, M. O., P. Merlet, Emission of trace gases and aerosols from biomass burning, Glob. Biogeochem. Cycles, 15, Issue 4, pages 955-966, December 2001. Article
  4. Giglio, L., J. T. Randerson, and G. R. van der Werf, Analysis of daily, monthly, and annual burned area using the fourth-generation global fire emissions database (GFED4)", J. Geophys. Res, Biogeosciences., 118, Issue 1, 317-328, March 2013. Article
  5. Randerson, J. T., Y. Chen1, G. R. van der Werf, B. M. Rogers, and D. C. Morton, Global burned area and biomass burning emissions from small fires, J. Geophys. Res, Biogeosciences, 117, Issue G4, December 2012. Article
  6. van der Werf, G. R, J. T. Randerson, L. Giglio, G. J. Collatz, M. Mu, P. S. Kasibhatla, D. C. Morton, R. S. DeFries, Y. Jin, and T. T. van Leeuwen, Global fire emissions and the contribution of deforestation, savanna, forest, agricultural, and peat fires (1997-2009) Atmos. Chem. Phys., 10, 11707-11735, 2010. Article
  7. van der Werf, G. R., J. T. Randerson, L. Giglio, T. T. van Leeuwen, Y. Chen, B. M. Rogers, M. Mu, M. J. E. van Marle, D. C. Morton, G. J. Collatz, R. J. Yokelson, and P. S. Kasibhatla, Global fire emissions estimates during 1997–2016, Earth Sys. Sci. Data, 9, 697-720, 2017. Article

--Bob Y. 10:56, 17 March 2015 (EDT)

Previous issues that are now resolved

Update GFED4 emission factors

This update was included in v11-02f (approved 17 May 2018).

Katie Travis wrote:

I was looking through the GFED4 emission factors in v11, and realized they appear out of date from is currently on the website ( This table appears to have been posted last June.
I put the table online into the HEMCO format using Bob’s code. I also updated the burned area categories in the code to match GFED4 to avoid future confusion. Since the GFED3 burned area categories are different as you say, I was not sure how to deal with a few things like Hg0 etc that are needed for the model but not standard for GFED. I cc’d Emily since she also made some changes for her PAN paper that I am not sure are properly represented anymore.

Prasad Kasibhatla wrote:

When I used GFED3 emission factors with GFED4 burned area files, I did something pretty simple. Here is the mapping as I remember:
  GFED4 Burned Area Category                GFED3 Emission Factor Category

  SAVA                                      SAV
  PEAT                                      PET
  BORF                                      FOR
  AGRI                                      AGW
  TEMP                                      FOR
  DEFO                                      DEF
Note in GFED4 implementation of GFED3 emission factors, there is no need of Humid Tropical Forest mask. In GFED3 implementation of GFED3 emission factors, GFED3 burned area category DEF used DEF emission factor for Humid Tropical Forest, and WDL emission factor outside Humid Tropical Forest.
So, now GFED4 emission factors corresponding to GFED4 burned area categories should be used. But the GFED3 code should not be changed. But you seem to gave changed the GFED3 part of the code?
In fact, we should now just get rid of the GFED3 option.

Katie Travis wrote:

I translated Hg0, (and HAC and GLYX, which are not in GFED4), from GFED3 to GFED4 using Prasad’s translation. I stripped out any calls to GFED3 in the code.

--Bob Yantosca (talk) 15:55, 17 May 2018 (UTC)

Fixed unit conversion error on 18 Jun 2015

On 18 Jun 2015, we fixed a unit conversion error in the GFED4 data files that were contained in the HEMCO data directory GFED4/v2015-06/:

Please note that this correction is now obsolete. GFED4/v2015-06 is superceded by the new GFED4 release located at GFED4/v2015-10.

--Lizzie Lundgren (talk) 16:55, 13 December 2017 (UTC)

Zhen Qu wrote:

The biomass burning emission seems to be smaller by about 3 magnitude with the conversion error.

Christoph Keller wrote:

Yes, this sounds about right. The GFED-4 netCDF files that we use in GEOS-Chem are directly generated from the original GFED-4 data. These data changed units between the provisional and the official GFED 4 release (gDM/m2/month to kgDM/m2/month). I didn’t realize that and the GFED-4 files were hence a factor of 1000 too low.

The GFED4 data files with the correct units were uploaded to the HEMCO data directory tree at 4PM EDT on 18 Jun 2015. If you obtained GFED4 data prior to 18 Jun 2015, we ask you to re-download this data.

To download the data, you can use the hemco_data_download installer, as described on our Downloading the HEMCO data directories wiki post. Or if you would like to get these files directly, you can use this FTP command:

ftp "*.nc"

We apologize for the inconvenience.

--Bob Y. (talk) 20:11, 18 June 2015 (UTC)

Currently unresolved issues

None at this time.

--Bob Y. 14:27, 27 March 2015 (EDT)