MEGAN v2.1 biogenic emissions

From Geos-chem
Jump to: navigation, search

GEOS-Chem currently uses MEGAN v2.1 plus Guenther 2012 biogenic emissions via HEMCO.

On this page, we describe the MEGAN v2.1 biogenic emissions, which were used in GEOS-Chem v8-02-04 through GEOS-Chem v9-02.


Mike Barkley has written a document entitled Description of MEGAN biogenic VOC emissions in GEOS-Chem. Please see this document for a full description of the MEGAN v2.1 emissions.

For more information, please follow this discussion about the implementation of MEGAN 2.1:

Mike Barkley wrote:

I've read in the latest GEOS-Chem news letter that your working on terpene emissions from MEGAN. I've got the latest v2.1 emission factors for these species (+isoprene & MBO) already included in my model code (v8-01-04) which I drive with MODIS LAI collection 5 (a supposedly much improved version) & which is divided by the fraction of vegetated area of each grid cell.
Be aware that the current implementation of the emissions in Geos-Chem is 'wrong' as it takes the LAI into account twice: once in the canopy model & once in the seasonality factor (although somehow the global total comes out okay ~500 Tg isoprene year). I've spoken to Alex about this, if the LAI is included in the canopy model, then the factor: 0.49d0 * LAI / SQRT( 1.d0 + 0.2d0 * LAI*LAI ) should not be applied as well.
I now run the emissions with two different algorithms:
  1. Using the existing canopy model but with the megan module restructured & with new leaf-age & temperature algorithms.
  2. Using the new PCEEA approach (which doesn't use a canopy model) outlined in Gunther 2006 which now takes into account the light dependency of monoterpene emissions.
From these two parameterizations the 10 year (1996-2006) global averages using ND46 at 4x5 & GEOS-4 met data are:
  1. = 518 ± 26 Tg isoprene year (1sd)
  2. = 403 ± 21 Tg isoprene year (1sd)
Obviously, these are quite large differences which is why I've hesitated about putting them forward to the standard code. I've done some validation at a three tropical sites, at 2 sites the algorithms give pretty much the same and at the other there are about ~30%-ish differences in the emissions. With the PCEEA algorithm the emissions are much closer in total to the estimates of Muller et al (2008) which were ~410 Tg isoprene year. It makes a big difference to the HCHO columns!
What I am going to do is put an optional flag in the input.geos file to turn the PCEEA algorithm on/off as required. I've got to tidy all this code up, including the LAI data sets but once I've done this I can send it you & Harvard, as it needs checking thoroughly.
--Dbm 13:10, 17 June 2009 (EDT)

Yuhang Wang wrote:

Jun Nam in my group has done some updates with isoprene emissions in GEOS-Chem. He's using MEGAN v2.1 with PCEEA (and w/o canopy model). The biggest problem he's got so far is when we switched from GEOS-4 to GEOS-5 met fields. Isoprene emissions are much higher. Apparently, surface temperature in GEOS-5 is too high (for 2006; we are using the current release of GEOS-Chem and GEOS-5 met field). He's sent emails to Steve Pawson but didn't get a response yet. So he used NCEP surface temperature instead and the global isoprene emission rate is 562 Tg C/year.
Can you set up a thread on wiki for this discussion so that Jun can provide more details if anyone wants to know.
--Dbm 11:00, 18 June 2009 (EDT)

Jun Nam wrote:

GEOS-5 has higher temperature than NCEPv2 reanalysis for 2006. Global isoprene emission is 658 Tg C (GEOS-5) vs. 562 Tg C (NCEPv2). For this analysis, I used PCEEA algorithm (Guenther et al., 2006) and MEGAN V2.04 (most recent version). I've checked high HCHO regions in 5 continents (except australia) and the surface temperature difference (in July afternoon) is 2-3 degrees (temporal variation agrees very well). As a result, activity factor due to temperature (gamma_t) is about 30-45% higher using GEOS-5 than NCEPv2.
(...) I'm using version 2.0 of the MEGAN LAIv (LAIv data described in Guenther et al. 2006), which is in netcdf_30min format available in MEGAN download portal (CDP) input folder.

Mike Barkley wrote:

I've done a quick comparison between the emissions resulting from GEOS-4 and GEOS-5 met data for the year 2004 only, using the PECCA model approach:
   GEOS-4 ==> 2004  ::  362 Tg C  =  410 Tg isoprene
   GEOS-5 ==> 2004  ::  338 Tg C  =  383 Tg isoprene
which is about ~7% difference. This difference maybe due to variations in LAI as well as in the meteorology.
--Dbm 14:56, 24 June 2009 (EDT)
There seems to be some contradiction here between what Jun is finding for GEOS-5 2006 (658 TgC) and what Mike is finding for GEOS-5 2004 (338 TgC).

Dylan Millet wrote:

Michael Barkley has packaged his updates to MEGAN in GEOS-Chem into a beta/non-standard version of the model (v8-2-1), and is very kindly making it available. There's a description above, but briefly it drives emissions of isoprene, terpenes, and MBO with MODIS LAI collection 5 and uses updated emission factors and emission algorithms (MEGAN v2.1). It also includes an option to use the canopy model currently implemented in GEOS-Chem or an alternative approach (PCEEA) described by Guenther et al. [ACP 2006].
For anyone interested in testing/using this code, you can find it as a gzipped tarball here:
Once you unpack it you'll find a README detailing changes from standard code.
If you do use it please be sure to get in touch with Mike (and also let me know, so I can see how much of a demand there is) - and if you use it in a paper it would be appropriate to offer Mike co-authorship.
--Dbm 17:54, 10 July 2009 (EDT)

GEOS-Chem v8-02-04 is the first "standard" code version to include the MEGAN v2.1 emissions. Please see the Overview section below for full details.

Speciated biogenic emissions

The MEGAN v2.1 inventory contains speciated emissions for several secondary organic aerosols. Please see this post on our Secondary Organic Aerosols wiki page for more information.

--Bob Y. 14:02, 19 March 2010 (EDT)

Met field requirements for MEGAN

The current MEGAN v2.1 emissions code computes a 10-day running average of the surface temperature field. At the start of a GEOS-Chem simulation, the MEGAN emissions code will attempt to read the previous 10 days's worth of surface temperature data. If you do not have the met field data files that contain the surface temperature field on disk for the 10 days prior to the start of your GEOS-Chem simulation, then GEOS-Chem will die with an error..

P.S. The now-obsolete MEGAN 2.0 emissions code used to keep a 15-day running average of surface temperature.

--Bob Y. 10:12, 8 March 2011 (EST)

New Update Required in MEGAN algorithm

Chris Chan Miller found that the isoprene emissions in MEGAN do not capture the diurnal variability of isoprene observed in the Michigan state forest. This is because the meteorological variables input to GEOS-Chem (temperature and photosynthetically active radiation) are not interpolated between the 3-h average surface data read in to megan_mod.f. This can be fixed by including a CALL INTERP statement in main.f for I-3 met fields. Chris is working on updating this in his working code and will send necessary updates to Bob to be included in the standard code.

--Emarais 16:32, 4 April 2012 (EDT)


  1. Barkley, M., Description of MEGAN biogenic VOC emissions in GEOS-Chem, 2010. PDF
  2. Buermann, W., Wang, Y.J., Dong, J.R., Zhou, L.M., Zeng, X.B., Dickinson, R.E., Potter, C.S., and Myneni, R.B.: Analysis of a multiyear global vegetation leaf area index data set, J. Geophys. Res., 107, 4646, doi:10.1029/2001JD000975, 2002.
  3. Guenther, A., Baugh, B., Brasseur, G., Greenberg, J., Harley, P., Klinger, L., Serca, D., and Vierling, L.: Isoprene emission estimates and uncertainties for the Central African EXPRESSO study domain, J. Geophys. Res., 104, 30625-30639, 1999.
  4. Guenther, A., Karl, T., Harley, P., Wiedinmyer, C., Palmer, P.I., and Geron, C.: Estimates of global terrestrial isoprene emissions using MEGAN (Model of Emissions of Gases and Aerosols from Nature), Atmos. Chem. Phys., 6, 3181-3210, 2006.
  5. Guenther, A., and C. Wiedinmyer, User's guide to the Model of Emissions of Gases and Aerosols from Nature (MEGAN), Version 2.01, 2007.
  6. Guenther, A. B., Jiang, X., Heald, C. L., Sakulyanontvittaya, T., Duhl, T., Emmons, L. K., and Wang, X.: The Model of Emissions of Gases and Aerosols from Nature version 2.1 (MEGAN2.1): an extended and updated framework for modeling biogenic emissions, Geosci. Model Dev., 5, 1471-1492, doi:10.5194/gmd-5-1471-2012, 2012.
  7. Millet, D.B., Jacob, D.J., Boersma, K.F., Fu, T.M., Kurosu, T.P., Chance, K., Heald, C.L., and Guenther, A.: Spatial distribution of isoprene emissions from North America derived from formaldehyde column measurements by the OMI satellite sensor, J. Geophys. Res., 113, D02307, doi:10.1029/2007JD008950, 2008. PDF
  8. Mueller, J.-F., et al. Global isoprene emissions estimated using MEGAN, ECMWF analyses and a detailed canopy environment model, Atmos. Chem. Phys., 8, 1329-1341, 2008.
  9. Myneni, R. B., et al., Large seasonal swings in leaf area of Amazon rainforests, Proceedings of the National Academy of Sciences, 104(12), 4820{4823, doi:10.1073/pnas.0611338104, 2007.
  10. Palmer, P.I., Abbot, D.S., Fu, T.M., Jacob, D.J., Chance, K., Kurosu, T.P., Guenther, A., Wiedinmyer, C., Stanton, J.C., Pilling, M.J., Pressley, S.N., Lamb, B., and Sumner, A.L.: Quantifying the seasonal and interannual variability of North American isoprene emissions using satellite observations of the formaldehyde column, J. Geophys. Res., 111, D12315, doi:10.1029/2005JD006689, 2006. PDF
  11. Sakulyanontvittaya, T., T. Duhl, C. Wiedinmyer, D. Helmig, S. Matsunaga, M. Potosnak, J. Milford, and A. Guenther, Monoterpene and Sesquiterpene Emission Estimates for the United States, Environ. Sci. Technol., 42(5), 1623{1629, doi:10.1021/es702274e, 2008.

--Bob Y. 09:39, 3 March 2010 (EST)