Difference between revisions of "Biomass burning emissions"

From Geos-chem
Jump to: navigation, search
(GFED3 monthly data)
(update FIRECAM details)
(47 intermediate revisions by 3 users not shown)
Line 1: Line 1:
 
This page describes the options for biomass burning emissions in GEOS-Chem.
 
This page describes the options for biomass burning emissions in GEOS-Chem.
  
== GFED3 ==
+
== GFED4 ==
  
=== GFED3 monthly data ===
+
The [[GFED4 biomass burning emissions|GFED4 biomass burning emission inventory]] is now the current default biomass burning dataset for GEOS-Chem (introduced in [[GEOS-Chem v10-01|v10-01]], June 2015).  The current GFED4 archive contains data from 1997 through 2016.
  
GFED3 is availalble for use in [[GEOS-Chem v9-01-02]] and higher versions.  Prasad Kasibhatla and the GEOS-Chem Support Team implemented GFED3 into [[GEOS-Chem v9-01-02]].
+
GFED4 has since superseded the obsolete [[GFED3 biomass burning emissions|GFED3]] and [[GFED2 biomass burning emissions|GFED2]] inventories.
  
For more information about the GFED3 data set, please see the following:
+
--[[User:Bmy|Bob Yantosca]] ([[User talk:Bmy|talk]]) 19:18, 2 January 2019 (UTC)
 
+
#[http://www.globalfiredata.org/ GFED web page].
+
#[http://www.falw.vu/~gwerf/GFED/ Guido van der Werf's GFED3 web page].
+
#[http://www.atmos-chem-phys-discuss.net/10/16153/2010/acpd-10-16153-2010.html ''van der Werf et al'' (2010)].  This is the  reference paper for the GFED3 inventory.
+
 
+
'''''[[User:Psk9|Prasad Kasibhatla]] wrote:'''''
+
 
+
<blockquote>I have updated the GFED3 database to include data for 2011.  Please note that GPCP precip data (one of the input datasets needed for GFED3) was not available for all of 2011 - the ERA-Iterim precip  (http://climatedataguide.ucar.edu/reanalysis/era-interim) was therefore used for the second half of 2011.</blockquote>
+
<blockquote>There were also small isolated changes prior to 2011 in this version of GFED3. I have therefore reprocessed the entire 1997-2011 GFED3 dataset for GEOS-Chem.</blockquote>
+
 
+
==== GFED3 in GEOS-Chem v10-01 and newer versions ====
+
 
+
In [[GEOS-Chem v10-01]] and newer versions, the GFED3 emissions of SO2 are read with the [[HEMCO|HEMCO emissions component]].  We have created [http://ftp.as.harvard.edu/gcgrid/data/ExtData/HEMCO/GFED3/v2014-10/README new GFED3 data files] (in [http://ferret.wrc.noaa.gov/noaa_coop/coop_cdf_profile.html COARDS-compliant netCDF format]) for use with HEMCO.  These new data files are contained in the [[HEMCO data directories|HEMCO data directory tree]].  For detailed instructions on how to download these data files to your disk server, please see our [[HEMCO_data_directories#Downloading_the_HEMCO_data_directories|''Downloading the HEMCO data directories'' wiki post]].
+
 
+
==== GFED3 in GEOS-Chem versions prior to v10-01 ====
+
 
+
GFED3 data is available from 1997-2011.  In GEOS-Chem v9-02 and older versions, the GFED3 data files are stored in binary punch format at this directory tree:
+
 
+
ftp ftp.as.harvard.edu
+
cd gcgrid/geos-chem/data/GEOS_NATIVE/GFED3_201212/
+
 
+
--[[User:Bmy|Bob Y.]] 14:02, 31 March 2011 (EDT)<br>
+
--[[User:Melissa Payer|Melissa Payer]] 15:03, 21 February 2012 (EST)
+
 
+
The data is read from disk by module <tt>GeosCore/gfed3_biomass_mod.F</tt>.
+
 
+
--[[User:Bmy|Bob Y.]] 15:12, 13 February 2015 (EST)
+
 
+
=== GFED3 daily and 3-hourly emissions ===
+
 
+
==== Original implementation ====
+
 
+
'''''This update was tested in the 1-month benchmark simulation [[GEOS-Chem v9-01-03 benchmark history#v9-01-03l|v9-01-03l]] and approved on 18 May 2012.'''''
+
 
+
'''''[mailto:psk9@duke.edu Prasad Kasibhatla] wrote:'''''
+
 
+
:I have updated GEOS-Chem to incorporate the capabilty to use GFED3 daily and 3-hourly fire fractions, and to regrid on the fly from the GFED3 native resolution (0.5x0.5) to the model resolution.
+
 
+
:NOTES:
+
 
+
:*Daily and 3hourly fractions are only available for 2003-2010. For methodology, see
+
::: Mu, M., J.T. Randerson, G.R. van der Werf, L. Giglio, P. Kasibhatla, D. Morton, G.J. Collatz, R.S. DeFries,  E.J. Hyer, E.M. Prins, D.W.T. Griffith, D. Wunch, G.C. Toon, V. Sherlock, and P.O. Wennberg. 2011. ''Daily and 3-hourly variability in global fire emissions and consequences for atmospheric model predictions of carbon monoxide''. <u>Journal of Geophysical Research-Atmospheres</u>. '''116''': D24303. doi:10.1029/2011JD016245.
+
:*Since fire fraction are not available as a function of type of fire type (as is the case for dry matter burnt), the same fractions are applied to all 6 fire types (ag waste, deforestation, extratropical forest, peat, savanna, and woodland).
+
:*I took the opportunity to rewrite <tt>gfed3_biomass_mod.F</tt> which mainly involved cleaning up code but with one big change - regridding directly from the GFED3 native resolution to the model resolution by adapting the MAP_A2A regriddng package mentioned in Bob's Model Engineer [[GEOS-Chem_Newsletter_Winter_2012#Other_development_tasks|Report]].
+
:*As a result of the change to the regridding, the GFED3 data file archive has to be updated to bpch files at the native resolution. I have these ready and will supply to [[GEOS-Chem Support Team|GC support team]]. Note that these require the creation of new directory structure since they do not belong in the <tt>GEOS_1x1</tt> data directory.
+
:*GFED3 for 2011 is likely to be available in late spring-early summer. An added data layer for fire emissions that may be missed by the burnt area algorithm is currently being developed.
+
:*There are occasional cells where the burnt dry matter seems to unrealistically high. For example in Dec 2009, the cell centered at 13.25 degree lat, 123.75 degree lon has DM = 13054 g DM/m<sup>2</sup>/month. The next highest grid box that month is about 723 g DM/m<sup>2</sup>/month. The GFED team is aware of this and is investigating.
+
 
+
--[[User:Melissa Payer|Melissa Payer]] 15:03, 21 February 2012 (EST)
+
 
+
'''''Update 3/9/12:'''''
+
 
+
'''''[mailto:psk9@duke.edu Prasad Kasibhatla] wrote:'''''
+
 
+
:I have updated the GFED3 patch file to:
+
 
+
:#revise the regridding-on-the-fly algorithm  - the orginal MAP_A2A algorithm results in a bit of incorrect smearing when regridding grid-average fluxes.
+
:#add code so that the GFED3 code works correctly for nested grid simulations
+
:#correct a small error - in a few tropical forest boxes, woodland (rather than tropical forest) emissions factros were being applied.
+
 
+
--[[User:Melissa Payer|Melissa Payer]] 16:13, 9 March 2012 (EST)
+
 
+
==== Bug fix for 3-hourly GFED3 emissions ====
+
 
+
'''''This update was tested in the 1-month benchmark simulation [[GEOS-Chem_v9-02_benchmark_history#v9-02l|v9-02l]] and approved on 26 Jun 2013.'''''
+
'''''This update is included in Adjoint [[GEOS-Chem_Adjoint_v35 | v35a]].'''''
+
 
+
'''''[mailto:skoplitz@fas.harvard.edu Shannon Koplitz] wrote:'''''
+
 
+
:I think there might be a bug in gfed3_biomass_mod.F when trying to run with 3-hourly GFED emissions. Running with 3-hourly emissions from the first day of a month seems to work fine, but trying to start on any other day appears to cause a problem. This didn't happen when running with daily emissions and I think it might have to do with how the monthly and 3-hourly emissions are coupled.
+
 
+
'''''[mailto:ckeller@seas.harvard.edu Christoph Keller] wrote:'''''
+
 
+
:I think I know what's going on. Since the same 3-hourly fractions are used throughout the entire month, the fractions are only read once per month.
+
 
+
:To construct the 3-hourly filename, the <tt>YYYYMMDD</tt> string is used (In gfed3_biomass_mod.F, line 670), which is hardcoded to have the day set to 1 (this happens on line 542). However, the <tt>TAU0</tt> value that is passed to <tt>READ_BPCH2_GFED3</tt> is constructed using the current day (<tt>DD</tt>), which can differ from 1 (<tt>DD</tt> is defined on line 505, the corresponding <tt>TAU0</tt> on line 668). So <tt>TAU0</tt> does not match the time stamp in the 3-hourly file if the starting day differs from the first of the month, and you get the error in <tt>READ_BPCH2_GFED3</tt> that you observe.
+
 
+
:If you hardcode <tt>DD</tt> to 1 in <tt>GET_TAU0</tt> on line 668 I would expect that the error disappears. So just use
+
 
+
TAU0 = GET_TAU0( MM, 1, YYYY, IT3H )
+
 
+
:Instead of
+
 
+
TAU0 = GET_TAU0( MM, DD, YYYY, IT3H )
+
 
+
--[[User:Melissa Payer|Melissa Payer]] 15:02, 16 April 2013 (EDT)
+
 
+
=== Update GFED3 emissions to 2011 ===
+
 
+
'''''This update was tested in the 1-month benchmark simulation [[GEOS-Chem_v9-02_benchmark_history#v9-02m|v9-02m]] and approved on 30 Jul 2013.'''''
+
'''''This update is included in Adjoint [[GEOS-Chem_Adjoint_v35 | v35d]].'''''
+
 
+
'''''[mailto:psk9@duke.edu Prasad Kasibhatla] wrote:'''''
+
 
+
:I have updated the GFED3 database to include data for 2011.  Please note that GPCP precip data (one of the input datasets needed for GFED3) was not available for all of 2011 - the ERA-Iterim precip (http://climatedataguide.ucar.edu/reanalysis/era-interim) was therefore used for the second half of 2011.
+
 
+
:There were also small isolated changes prior to 2011 in this version of GFED3. I have therefore reprocessed the entire 1997-2011 GFED3 dataset for GEOS-Chem, now stored in at Harvard in the:
+
 
+
  ftp//ftp.as.harvard.edu/gcgrid/data/GEOS_NATIVE/GFED3_201212/ 
+
 
+
:directory structure. 
+
 
+
:If you wish to use this updated GFED3 data, please also modify the following line in the module file <tt>GeosCore/gfed3_biomass_mod.F</tt>, from:
+
 
+
  CALL GFED3_AVAILABLE( YYYY, 1997, 2010 )
+
 
+
:to
+
 
+
  CALL GFED3_AVAILABLE( YYYY, 1997, 2011 )
+
 
+
:Please note that this is the last update for GFED3.  Plans are underway for GFED4, whcih will ultimately replace the pre-2011 GFED3 dataset and go forward in time post-2011.
+
 
+
'''''[[User:Bmy|Bob Yantosca]] wrote:'''''
+
 
+
:Please note, we are placing this data into the <tt>GEOS_NATIVE/GFED3_201212</tt> subdirectory so that it will not conflict with the existing GFED3 data in <tt>GEOS_NATIVE/GFED3_201203</tt>.  By not overwriting the data, we will be able to replicate GEOS-Chem results made with prior model versions using the data in GFED3_201203.
+
 
+
:We will add support for the updated GFED3 data to the standard GEOS-Chem along with other emissions updates (slated for [[GEOS-Chem_v9-02]]), and we will run our customary 1-month benchmark to test the results.  In the meantime, if you would like to use the new GFED3 data for 2011, then please download the data to your location and make the above-mentioned update in gfed3_biomass_mod.F in your research code.
+
 
+
--[[User:Melissa Payer|Melissa Payer]] 18:02, 18 December 2012 (EST)
+
 
+
=== GFED3 data files and GAMAP ===
+
 
+
If you want to read the GFED3 dry matter files with GAMAP, please make sure you add the following lines to these files:
+
 
+
In <tt>IDL/gamap2/input_files/diaginfo.dat</tt>, add the line:
+
 
+
    8000 GFED3-BB                                GFED3 biomass burning
+
 
+
In <tt>IDL/gamap2/input_files/tracerinfo.dat</tt> files, add the lines:
+
 
+
#==============================================================================
+
# OFFSET=8000: GFED3 dry matter emissions
+
#==============================================================================
+
AGW_DM  AGW dry matter emissions            0.e0  1    8091  1.0E+00 g/m2/month
+
DEF_DM  DEF dry matter emissions            0.e0  1    8092  1.0E+00 g/m2/month
+
FOR_DM  FOR dry matter emissions            0.e0  1    8093  1.0E+00 g/m2/month
+
PET_DM  PET dry matter emissions            0.e0  1    8094  1.0E+00 g/m2/month
+
SAV_DM  SAV dry matter emissions            0.e0  1    8095  1.0E+00 g/m2/month
+
WDL_DM  WDL dry matter emissions            0.e0  1    8096  1.0E+00 g/m2/month
+
 
+
These modifications will be made in GAMAP v2-16 and higher versions.
+
 
+
--[[User:Bmy|Bob Y.]] 16:37, 30 January 2012 (EST)
+
  
 
== FINNv1 ==
 
== FINNv1 ==
  
<div style="color: #aa0000; background: #eeeeee;border: 3px solid red; padding: 1em; margin: auto; width: 90%; ">'''<p>NOTE: At present, the FINN biomass emissions inventory has been added to the non-standard SEAC4RS research version of GEOS-Chem.  We plan to add the FINN biomass emissions into the standard GEOS-Chem code via the [[HEMCO|HEMCO emissions component]], most likely in version [[GEOS-Chem v10-01|v10-01]]. (Bob Yantosca, 13 May 2014)</p>'''</div>
+
The [[FINNv1 biomass burning emissions|FINNv1 biomass burning emission inventory]] has been added to [[GEOS-Chem v10-01]] via the [[HEMCO|HEMCO emissions component]]. This inventory may be used to replace [[GFED4 biomass burning emissions|GFED4]] for research purposes.
 
+
"The Fire INventory from NCAR version 1.0 (FINNv1) provides daily, 1-km resolution, global estimates of the trace gas and particle emissions from open
+
burning of biomass, which includes wildfire, agricultural fires, and pre- scribed burning and does not include biofuel use and trash burning" (Wiedinmyer et al., 2011). For more information about this inventory, see:<br>
+
Wiedinmyer, C., S. K. Akagi, R .J. Yokelson, L. K. Emmons, J. A. Al-Saadi, J. J. Orlando, and A. J. Soja. ''The Fire INventory from NCAR (FINN): a high resolution global model to estimate the emissions from open burning''. Geosci. Model Dev., 4, 625–641, 2011. [http://www.geosci-model-dev.net/4/625/2011/]
+
 
+
=== Implementation in GEOS-Chem ===
+
Jenny Fisher (U. Wollongong) and Min Huang (JPL) are working with Christine Wiedinmyer to implement FINNv1 into the standard version of GEOS-Chem. This implementation will include a number of features, many of which are improvements over the non-standard implementations currently used (see below):
+
#Improved speciation for GEOS-Chem, including specialty simulations like SOA
+
#Compatibility with new emissions module (including use of netcdf, high-resolution inputs, online regridding, etc.)
+
#GFED-like online computation of emissions for most species based on emissions of a single species for each land type (significantly reducing file storage requirements)
+
#Ability to use daily or monthly emissions
+
 
+
If you have questions or suggestions, please contact Jenny Fisher (jennyf@uow.edu.au) and Min Huang (Min.Huang@jpl.nasa.gov).
+
--[[User:Bmy|Bob Y.]] 10:54, 13 May 2014 (EDT)
+
 
+
'''UPDATE (July 2014)'''
+
 
+
A new version of FINNv1 emissions has been implemented in GEOS-Chem v9-02 and is currently being incorporated into HEMCO for inclusion into v10-01. The new implementation starts from CO2 emissions (gridded to high resolution from original FINN files) for 6 different land types. CO2 emissions are converted to emissions of other gas phase and aerosol species using emission ratios provided by Christine Wiedinmyer. Emissions of non-methane organic compounds are derived from a single NMOC species (also computed from ratios to CO2) using speciation factors (also provided by Christine Wiedinmyer). These calculations are performed online in GEOS-Chem, then summed over land types and regridded to model resolution. FINN can be used with daily or monthly resolution, and emissions will be available with a lag of approximately 6 months behind real time (in special cases, near-real-time emissions can also be provided).
+
 
+
This implementation includes support for many species, including those used in specialised GEOS-Chem simulations. Some of these species are not currently emitted in publicly available releases of GEOS-Chem, and interested users will be responsible for adding relevant emissions-related code (i.e. defining IDBs, etc.). No changes to the FINN code itself should be necessary to use these species. The current list of FINN species includes:
+
#Species currently emitted by biomass burning in standard versions for full chemistry and/or specialty simulations ''(no code changes necessary)'': CO2, CO, CH4, NO, SO2, OC, BC, NH3, ACET, ALD2, ALK4, BENZ, C2H2, C2H4, C2H6, C3H8, CH2BR2, CH2O, CH3BR, GLYC, GLYX, HAC, MEK, MGLY, PRPE, TOLU, XYLE
+
#Species either not currently emitted or emitted but not by biomass burning in standard versions ''(code changes necessary)'': ACTA, CH3CN, CH3I, DMS, EOH, HCN, HCOOH, HNO2, ISOP, MACR, MNO3, MOH, MVK, R4N2, RCHO, LIMO, MTPA, MTPO
+
#Species currently lumped online in the FINN code: ETBENZ (with TOLU), STYR (with TOLU), TMB (with XYLE), APINE (with MTPA), BPINE (with MTPA), CARENE (with MTPA)
+
#Species/groupings that don't currently exist but that are defined in FINN inputs for possible implementation by interested users: AROM, FUR, ROH, RCOOH, SESQ
+
 
+
For the original lumping of VOCs into GEOS-Chem species, see [http://wiki.seas.harvard.edu/geos-chem/images/FINN_final_lumping.pdf FINN_final_lumping.pdf]
+
 
+
Potential future updates include more detailed treatment of terpenes and of secondary species like PAN. If you have an interest in working on emissions of these species, please contact me (jennyf@uow.edu.au) and/or Christine (christin@ucar.edu).
+
 
+
--Jenny Fisher, 4 June 2014
+
 
+
=== Current non-standard implementations ===
+
As of IGC6 (May 2013), FINNv1 has been implemented in non-standard versions of GEOS-Chem by at least 3 groups. These implementations allow users to use the ASCII emissions files (or a converted gridded netcdf equivalent) available directly from http://bai.acd.ucar.edu/Data/fire/. If you are interested in working with FINNv1 emissions immediately, please contact Jingqiu Mao (Jingqiu.Mao@noaa.gov), Min Huang (Min.Huang@jpl.nasa.gov), or Jenny Fisher (jennyf@uow.edu.au).
+
 
+
--[[User:Bmy|Bob Y.]] 10:54, 13 May 2014 (EDT)
+
 
+
== GFED2 ==
+
 
+
'''''NOTE: Once the GFED3 daily and hourly biomass burning emissions are installed into GEOS-Chem v9-01-03, then GFED2 will be removed.'''''
+
 
+
=== GFED2 monthly data thru 2008 ===
+
 
+
As of July 2009, the GFED2 monthly data is now available for years 1997-2008.  The data is available in bpch format from the following directory structure:
+
 
+
ftp ftp.as.harvard.edu
+
cd pub/geos-chem/data/GEOS_1x1/GFED2_200601/
+
 
+
For more information about the data, please see the following README file: [ftp://ftp.as.harvard.edu/pub/geos-chem/data/GEOS_1x1/GFED2_200601/README GEOS_1x1/GFED2_200601/README]
+
 
+
You will need to make a small modification in routine GFED2_COMPUTE_BIOMASS (in <tt>gfed2_biomass_mod.f</tt>) in order to update the maximum available year of data.  Look for the lines:
+
 
+
      ! Availability of MONTHLY data
+
      !-------------------------------
+
      ELSE IF ( LGFED2BB ) THEN
+
       
+
        CALL GFED2_AVAILABLE( YYYY, 1997, 2007 )
+
 
+
and change the 2007 to 2008:
+
 
+
      ! Availability of MONTHLY data
+
      !-------------------------------
+
      ELSE IF ( LGFED2BB ) THEN
+
       
+
        !-----------------------------------------------------------
+
        ! Prior to 7/8/09:
+
        ! GFED2 2008 monthly data is now available (bmy, 7/8/09)
+
        !CALL GFED2_AVAILABLE( YYYY, 1997, 2007 )
+
        !-----------------------------------------------------------
+
        CALL GFED2_AVAILABLE( YYYY, 1997, 2008 )
+
 
+
You may also download the updated version of <tt>gfed2_biomass_mod.f</tt> from:
+
 
+
ftp ftp.as.harvard.edu
+
cd pub/geos-chem/patches/v8-02-02
+
get gfed2_biomass_mod.f
+
 
+
--[[User:Bmy|Bob Y.]] 10:32, 24 February 2010 (EST)
+
 
+
=== GFED2 8-day emissions ===
+
 
+
The GFED2 8-day emissions are now included in the standard model as a user option as of [[GEOS-Chem v8-02-01]]. The data files for years 2001-2007 may be downloaded from:
+
 
+
ftp ftp.as.harvard.edu
+
cd pub/geos-chem/data/GEOS_1x1/GFED2_8day_200712/
+
 
+
For more information on the GFED2 8-day emissions data, please see the README file: [ftp://ftp.as.harvard.edu/pub/geos-chem/data/GEOS_1x1/GFED2_8day_200712/README GEOS_1x1/GFED2_8day_200712/README]
+
 
+
'''''[mailto:ray.nassar@utoronto.ca Ray Nassar] wrote:'''''
+
 
+
:The [GFED2 8-day emissions] code is working and I am using it for my work.  Overall, the results I get for CO and ozone are very similar to my monthly GFED runs but this changes if one focuses on a  specific region over with a short time scale.
+
 
+
:{Here are} ... [http://acmg.seas.harvard.edu/geos/wiki_docs/emissions/GFED2_compare.pdf Some of {my comparisons}.]  This is by no means a complete analysis but more of an initial assessment.
+
 
+
:In my opinion, the move to 8-day emissions is a definite improvement, although as I said, the  differences are fairly localized and really only matter for time scales shorter than one month.
+
 
+
--[[User:Bmy|Bob Y.]] 10:32, 24 February 2010 (EST)
+
 
+
== FLAMBE biomass emissions ==
+
 
+
<div style="color: #aa0000; background: #eeeeee;border: 3px solid red; padding: 1em; margin: auto; width: 90%; ">'''<p>NOTE: At this time FLAMBE has not been implemented into the mainline standard GEOS-Chem, but it has been used for the NRT-ARCTAS codes.  Therefore probably only a small fraction of GEOS-Chem users will be concerned with these emissions at this time. (Bob Yantosca, 13 May 2014)</p>'''</div>
+
 
+
=== Data issues ===
+
  
'''''[mailto:jafisher@fas.harvard.edu Jenny Fisher] wrote:'''''
+
--[[User:Melissa Payer|Melissa Sulprizio]] 10:21, 18 February 2015 (EST)
  
:Sorry to revisit the old FLAMBE issues, but we are seeing some strange results in the fire emissions, and I am trying to verify with Ed Hyer that we are processing the data correctly.
+
== QFED ==
  
:To get a handle on how we are processing the file, I have been looking in Philippe's directory <tt>~phs/IDL/dvpt/flambe/</tt>
+
The [[QFED_biomass_burning_emissions|QFED biomass burning emission inventory]] has been added to [[GEOS-Chem v10-01]] via the [[HEMCO|HEMCO emissions component]]. This inventory may be used to replace [[GFED4 biomass burning emissions|GFED4]] for research purposes.
  
:I took a look at some of the data files in the data/ directory in there. In the column that corresponds to carbon emissions (column 11, or 10 in IDL accounting), I see values that range from 4.95e4 to 765e4. Hyer tells me that in his version of the files, these range from 4.95 to 7650 (i.e. 4 orders of magnitude different). We are assuming these emissions are in g/m2, which they clearly aren't at our values.
+
--[[User:Melissa Payer|Melissa Sulprizio]] ([[User talk:Melissa Payer|talk]]) 18:27, 14 September 2016 (UTC)
  
:Do either of you know if I am looking at old files, and if what we have looks more like what he has? Or are we processing things 4 orders of magnitude too large?? Or, are we assuming a different unit on the emissions when we actually process them? I can't seem to find raw flambe data files anywhere besides Philippe's development directory...
+
== GFAS ==
  
'''''[mailto:plesager@seas.harvard.edu Philippe Le Sager] wrote:'''''
+
The [[GFAS biomass burning emissions|GFAS biomass burning emission inventory]] was added in [[GEOS-Chem 12#12.2.0|GEOS-Chem 12.2.0]] via the [[HEMCO|HEMCO emissions component]]. This inventory may be used to replace [[GFED4 biomass burning emissions|GFED4]] for research purposes.
  
:We had a couple of problems with the data. The one that gave me headache was the difference between forecast and analysis data, which can be mixed in some files. There was the issue of double counting fires seen by both GOES and MODIS, and we found a problem with unit.  It was Kg instead of g as advertised in the ppt [that was given to us] (the only documentation for the data).  Ed did correct the files and back processed the data on its server.
+
--[[User:Melissa Payer|Melissa Sulprizio]] ([[User talk:Melissa Payer|talk]]) 12:19, 4 October 2018 (UTC)
  
:I have a "new" directory in my test data. You probably look at the old set of data. The new one is in the "new" subdirectory:
+
== FIRECAM tool ==
  
    /home/phs/IDL/dvpt/flambe/data/new/
+
I developed the [https://globalfires.earthengine.app/view/firecam FIRECAM tool] for end-users to quickly compare fire emissions estimates from five global inventories (GFEDv4s, GFASv1.2, FINNv1.5, QFEDv2.5r1, and FEERv1.0-G1.2) for a given study region. 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 [https://sites.google.com/view/firecam/home FIRECAM website] and [https://doi.org/10.1016/j.rse.2019.111557 our paper] for more information. The paper describes the methodological differences among these inventories and metrics we developed to diagnose spatial biases and uncertainties. We recommend end-users to use multiple inventories in model runs, if possible. As we show in our paper, the choice of emissions inventory can significantly impact model results. If this is not feasible, please take the time to see how much emissions estimates vary in your study region using the FIRECAM tool.
  
:You can also still get the data at the following website:
+
--[[User:tianjialiu|Tianjia Liu]] ([[User talk:tianjialiu|talk]]) 15:30, 21 August 2019 (EDT)
  
    http://www.nrlmry.navy.mil/aerosol_web/arctas_flambe/data_hourly/
+
== Previous issues that are now resolved ==
  
:Finally here are some totals we check w/ Maria and Jingqiu:
+
=== Fix bug in biomass burning emissions when GFED is turned off ===
  
    Data for one set of satellites was still given Kg.  Now fixed, I have the following
+
<span style="color:green">'''''This update was included in [[GEOS-Chem v11-02#v11-02f|v11-02f]] (approved 17 May 2018).'''''</span>
    total for one day all over the world, assuming 3 land types. Seems reasonable, no?
+
+
    NOx (Tg N) assuming ALL is savanna / trop forest / extraTrop forest:
+
      0.080809369    0.063615890      0.10316091
+
+
    CO (Tg) assuming ALL is savanna / trop forest / extraTrop forest:
+
      4.6422410      7.5160097      7.8107557
+
+
    SMOKE (Tg) assuming ALL is savanna / trop forest / extraTrop forest:
+
      0.47039473      0.47039473      0.47039473
+
  
--[[User:Bmy|Bob Y.]] 15:07, 18 February 2009 (EST)
+
'''''Corey Trujillo (U. Colorado) wrote:'''''
  
=== Obtaining the FLAMBE data ===
+
:I've been learning [[GEOS-Chem v11-01]] and have encountered two anomalies (one a bug, one a potential bug) when I turn off GFED.  These also appear to persist in v11-02d.
  
==== Raw data ====
+
:The first issue is that the model defaults to using BOND biomass burning if GFED (or FINN) is off. This is a bit strange, especially since BOND biomass burning would then be turned on if a user turned off GFED and FINN but turned on QFED - this would lead to double biomass burning emissions for BC and OC. This can be avoided by defining a separate emission extension for BOND biomass burning in HEMCO_Config.rc, i.e. replacing this:
The FLAMBE biomass burning data files are available at the following archive: http://www.nrlmry.navy.mil/aerosol_web/arctas_flambe/data_hourly/.
+
  
[[User:Bmy|Bob Yantosca]] has some scripts that can be used to download the data (contact him for more info)These are:
+
  (((.not.GFED.or.FINN
 +
  0 BOND_BM_BCPI  $ROOT/BCOC_BOND/v2014-07/Bond_biomass.nc BC  2000/1-12/1/0 C xy kg/m2/s BCPI 70    2 1
 +
  0 BOND_BM_BCPO  -                                        -    -            - - -      BCPO 71    2 1
 +
  0 BOND_BM_OCPI  $ROOT/BCOC_BOND/v2014-07/Bond_biomass.nc OC  2000/1-12/1/0 C xy kg/m2/s OCPI 72    2 1
 +
  0 BOND_BM_OCPO  -                                        -    -            - -  -      OCPO 73    2 1
 +
  0 BOND_BM_POG1  -                                        -    -            - -  -      POG1 74/76 2 1
 +
  0 BOND_BM_POG2  -                                        -    -            - -  -      POG2 74/77 2 1
 +
  ))).not.GFED.or.FINN
  
;sleepFlambe: Perl script which issues a [[Downloading GEOS-Chem source code and data#Using wget to download files|Unix <tt>wget</tt>]] command to download one day (24 hourly files) of FLAMBE data to disk.
+
:with something like this:
;wrap_flambe.pro: IDL batch file which acts as a driver for looping over all 24 hours of FLAMBE data for a given date.
+
;read_flambe.pro: IDL program which reads each FLAMBE "raw" data file and saves it to [http://acmg.seas.harvard.edu/gamap/doc/Chapter_6.html#6.2 GEOS-Chem bpch format].
+
  
==== Processed data ====
+
  (((BOND_BBURN
The processed FLAMBE data (for input to GEOS-Chem) are kept on the Harvard data archive in the directories:
+
  0 BOND_BM_BCPI  $ROOT/BCOC_BOND/v2014-07/Bond_biomass.nc BC  2000/1-12/1/0 C xy kg/m2/s BCPI 70    2 1
 +
  0 BOND_BM_BCPO  -                                        -    -            - -  -      BCPO 71    2 1
 +
  0 BOND_BM_OCPI  $ROOT/BCOC_BOND/v2014-07/Bond_biomass.nc OC  2000/1-12/1/0 C xy kg/m2/s OCPI 72    2 1
 +
  0 BOND_BM_OCPO  -                                        -    -            - -  -      OCPO 73    2 1
 +
  0 BOND_BM_POG1  -                                        -    -            - -  -      POG1 74/76 2 1
 +
  0 BOND_BM_POG2  -                                        -    -            - -  -      POG2 74/77 2 1
 +
  )))BOND_BBURN
  
ftp ftp.as.harvard.edu
+
:The second issue is that if someone is to use the BOND biomass burning emissions, either intentionally or not, the emissions are not showing up in the BC-BIOB diagnostic. Instead, that diagnostic is all zeros, even though there is still clearly biomass burning emissions occurring in the model. Figures demonstrating this for a one month simulation of the standard code for July 2013, 4x5, GEOSFP, can be found here.
cd /pub/geos-chem/NRT-ARCTAS/flambe/YYYY/MM
+
  
where <tt>YYYY/MM</tt> are the year and month of the data.  Currently we have the full data archive from 2008 and 2009, with some months of 2010.
+
:There is thus a bug in the BC-BIOB diagnostic — has anyone else encountered this or have a suggestion for a fix?
  
--[[User:Bmy|Bob Y.]] 14:32, 10 February 2010 (EST)
+
'''''Therese Carter (MIT) wrote:'''''
  
== Duncan et al "Seasonal" Biomass Emissions ==
+
:A question came up recently about v11-01 defaulting to using BOND biomass burning emissions if QFED is switched on and both FINN and GFED are switched off. This can be fairly easily fixed by adding the following code in red:
  
See ''Bey et al'' [2001].
+
  (((.not.GFED.or.FINN<span style="color:red">.or.QFED2</span>
 +
  0 BOND_BM_BCPI  $ROOT/BCOC_BOND/v2014-07/Bond_biomass.nc BC  2000/1-12/1/0 C xy kg/m2/s BCPI 70    2 1
 +
  0 BOND_BM_BCPO  -                                        -    -            - -  -      BCPO 71    2 1
 +
  0 BOND_BM_OCPI  $ROOT/BCOC_BOND/v2014-07/Bond_biomass.nc OC  2000/1-12/1/0 C xy kg/m2/s OCPI 72    2 1
 +
  0 BOND_BM_OCPO  -                                        -    -            - -  -      OCPO 73    2 1
 +
  0 BOND_BM_POG1  -                                        -    -            - -  -      POG1 74/76 2 1
 +
  0 BOND_BM_POG2  -                                        -    -            - -  -      POG2 74/77 2 1
 +
  ))).not.GFED.or.FINN<span style="color:red">.or.QFED2</span>
  
== Duncan et al "Interannual" Biomass Emissions ==
+
--[[User:Bmy|Bob Yantosca]] ([[User talk:Bmy|talk]]) 16:02, 17 May 2018 (UTC)
  
 
== References ==
 
== References ==
Line 330: Line 91:
 
#Mu, M., J.T. Randerson, G.R. van der Werf, L. Giglio, P. Kasibhatla, D. Morton, G.J. Collatz, R.S. DeFries,  E.J. Hyer, E.M. Prins, D.W.T. Griffith, D. Wunch, G.C. Toon, V. Sherlock, and P.O. Wennberg, ''Daily and 3-hourly variability in global fire emissions and consequences for atmospheric model predictions of carbon monoxide'', <u>Journal of Geophysical Research-Atmospheres</u>, '''116''', D24303, doi:10.1029/2011JD016245, 2011.
 
#Mu, M., J.T. Randerson, G.R. van der Werf, L. Giglio, P. Kasibhatla, D. Morton, G.J. Collatz, R.S. DeFries,  E.J. Hyer, E.M. Prins, D.W.T. Griffith, D. Wunch, G.C. Toon, V. Sherlock, and P.O. Wennberg, ''Daily and 3-hourly variability in global fire emissions and consequences for atmospheric model predictions of carbon monoxide'', <u>Journal of Geophysical Research-Atmospheres</u>, '''116''', D24303, doi:10.1029/2011JD016245, 2011.
 
#Nassar, R., J. A. Logan, I. A. Megretskaia, L. T. Murray, L. Zhang, and D. B. A. Jones, ''Analysis of tropical tropospheric ozone, carbon monoxide and water vapor during the 2006 El Niño using TES observations and the GEOS-Chem model'', <u>J. Geophys. Res.</u>, '''114''', D17304, doi:10.1029/2009JD011760, 2009. [http://acmg.seas.harvard.edu/publications/2009/nassar_2009.pdf PDF]
 
#Nassar, R., J. A. Logan, I. A. Megretskaia, L. T. Murray, L. Zhang, and D. B. A. Jones, ''Analysis of tropical tropospheric ozone, carbon monoxide and water vapor during the 2006 El Niño using TES observations and the GEOS-Chem model'', <u>J. Geophys. Res.</u>, '''114''', D17304, doi:10.1029/2009JD011760, 2009. [http://acmg.seas.harvard.edu/publications/2009/nassar_2009.pdf PDF]
#van der Werf, G., 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)'', <u>Atmos. Chem. Phys. Discuss.</u>, '''10''', 16153-16230, 2010. [http://www.atmos-chem-phys-discuss.net/10/16153/2010/acpd-10-16153-2010.html PDF]
+
#Liu, T., L.J. Mickley, R.S. DeFries, M.E. Marlier, M.F. Khan, M.T. Latif, and A. Karambelas, ''Diagnosing spatial uncertainties and relative biases in global fire emissions inventories: Indonesia as regional case study'', <u>Remote Sens. Environ.</u>, '''237''', 111557. [https://doi.org/10.1016/j.rse.2019.111557 Article]
 +
#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'', <u>Earth Sys. Sci. Data</u>, '''9''', 697-720, 2017. [https://doi.org/10.5194/essd-9-697-2017 Article]
 +
#Wiedinmyer, C., S.K. Akagi, R.J. Yokelson, L.K. Emmons, J.J. Orlando, and A.J. Soja, ''The Fire INventory from NCAR (FINN): a high resolution global model to estimate the emissions from open burning'', <u>Geosci. Model Dev.</u>, '''4''', 625–641, 2011. [https://doi.org/10.5194/gmd-4-625-2011 Article]
 +
#Kaiser, J.W., A. Heil, M.O. Andreae, A. Benedetti, N. Chubarova, L. Jones, J.J. Morcrette, M. Razinger, M.G. Schultz, M. Suttie, and G.R. van der Werf, ''Biomass burning emissions estimated with a global fire assimilation system based on observed fire radiative power'', <u>Biogeosciences</u>, '''9''', 527–554, 2012. [https://doi.org/10.5194/bg-9-527-2012 Article]
 +
#Darmenov, A.S. and A. da Silva, ''The Quick Fire Emissions Dataset (QFED) - Documentation of versions 2.1, 2.2, and 2.4, NASA Technical Report Series on Global Modeling and Data Assimilation, Volume 32'', 2013. [http://citeseerx.ist.psu.edu/viewdoc/summary?doi=10.1.1.406.7724 Article]
 +
#Ichoku, C. and L. Ellison, ''Global top-down smoke-aerosol emissions estimation using satellite fire radiative power measurements'', <u>Atmos. Chem. Phys.</u>, '''14''', 6643–6667, 2014. [https://doi.org/10.5194/acp-14-6643-2014 Article]
  
 
--[[User:Melissa Payer|Melissa Payer]] 15:10, 21 February 2012 (EST)
 
--[[User:Melissa Payer|Melissa Payer]] 15:10, 21 February 2012 (EST)

Revision as of 07:30, 12 December 2019

This page describes the options for biomass burning emissions in GEOS-Chem.

GFED4

The GFED4 biomass burning emission inventory is now the current default biomass burning dataset for GEOS-Chem (introduced in v10-01, June 2015). The current GFED4 archive contains data from 1997 through 2016.

GFED4 has since superseded the obsolete GFED3 and GFED2 inventories.

--Bob Yantosca (talk) 19:18, 2 January 2019 (UTC)

FINNv1

The FINNv1 biomass burning emission inventory has been added to GEOS-Chem v10-01 via the HEMCO emissions component. This inventory may be used to replace GFED4 for research purposes.

--Melissa Sulprizio 10:21, 18 February 2015 (EST)

QFED

The QFED biomass burning emission inventory has been added to GEOS-Chem v10-01 via the HEMCO emissions component. This inventory may be used to replace GFED4 for research purposes.

--Melissa Sulprizio (talk) 18:27, 14 September 2016 (UTC)

GFAS

The GFAS biomass burning emission inventory was added in GEOS-Chem 12.2.0 via the HEMCO emissions component. This inventory may be used to replace GFED4 for research purposes.

--Melissa Sulprizio (talk) 12:19, 4 October 2018 (UTC)

FIRECAM tool

I developed the FIRECAM tool for end-users to quickly compare fire emissions estimates from five global inventories (GFEDv4s, GFASv1.2, FINNv1.5, QFEDv2.5r1, and FEERv1.0-G1.2) for a given study region. 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 and our paper for more information. The paper describes the methodological differences among these inventories and metrics we developed to diagnose spatial biases and uncertainties. We recommend end-users to use multiple inventories in model runs, if possible. As we show in our paper, the choice of emissions inventory can significantly impact model results. If this is not feasible, please take the time to see how much emissions estimates vary in your study region using the FIRECAM tool.

--Tianjia Liu (talk) 15:30, 21 August 2019 (EDT)

Previous issues that are now resolved

Fix bug in biomass burning emissions when GFED is turned off

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

Corey Trujillo (U. Colorado) wrote:

I've been learning GEOS-Chem v11-01 and have encountered two anomalies (one a bug, one a potential bug) when I turn off GFED. These also appear to persist in v11-02d.
The first issue is that the model defaults to using BOND biomass burning if GFED (or FINN) is off. This is a bit strange, especially since BOND biomass burning would then be turned on if a user turned off GFED and FINN but turned on QFED - this would lead to double biomass burning emissions for BC and OC. This can be avoided by defining a separate emission extension for BOND biomass burning in HEMCO_Config.rc, i.e. replacing this:
  (((.not.GFED.or.FINN
  0 BOND_BM_BCPI   $ROOT/BCOC_BOND/v2014-07/Bond_biomass.nc BC   2000/1-12/1/0 C xy kg/m2/s BCPI 70    2 1
  0 BOND_BM_BCPO   -                                        -    -             - -  -       BCPO 71    2 1
  0 BOND_BM_OCPI   $ROOT/BCOC_BOND/v2014-07/Bond_biomass.nc OC   2000/1-12/1/0 C xy kg/m2/s OCPI 72    2 1
  0 BOND_BM_OCPO   -                                        -    -             - -  -       OCPO 73    2 1
  0 BOND_BM_POG1   -                                        -    -             - -  -       POG1 74/76 2 1
  0 BOND_BM_POG2   -                                        -    -             - -  -       POG2 74/77 2 1
  ))).not.GFED.or.FINN
with something like this:
  (((BOND_BBURN
  0 BOND_BM_BCPI   $ROOT/BCOC_BOND/v2014-07/Bond_biomass.nc BC   2000/1-12/1/0 C xy kg/m2/s BCPI 70    2 1
  0 BOND_BM_BCPO   -                                        -    -             - -  -       BCPO 71    2 1
  0 BOND_BM_OCPI   $ROOT/BCOC_BOND/v2014-07/Bond_biomass.nc OC   2000/1-12/1/0 C xy kg/m2/s OCPI 72    2 1
  0 BOND_BM_OCPO   -                                        -    -             - -  -       OCPO 73    2 1
  0 BOND_BM_POG1   -                                        -    -             - -  -       POG1 74/76 2 1
  0 BOND_BM_POG2   -                                        -    -             - -  -       POG2 74/77 2 1
  )))BOND_BBURN
The second issue is that if someone is to use the BOND biomass burning emissions, either intentionally or not, the emissions are not showing up in the BC-BIOB diagnostic. Instead, that diagnostic is all zeros, even though there is still clearly biomass burning emissions occurring in the model. Figures demonstrating this for a one month simulation of the standard code for July 2013, 4x5, GEOSFP, can be found here.
There is thus a bug in the BC-BIOB diagnostic — has anyone else encountered this or have a suggestion for a fix?

Therese Carter (MIT) wrote:

A question came up recently about v11-01 defaulting to using BOND biomass burning emissions if QFED is switched on and both FINN and GFED are switched off. This can be fairly easily fixed by adding the following code in red:
  (((.not.GFED.or.FINN.or.QFED2
  0 BOND_BM_BCPI   $ROOT/BCOC_BOND/v2014-07/Bond_biomass.nc BC   2000/1-12/1/0 C xy kg/m2/s BCPI 70    2 1
  0 BOND_BM_BCPO   -                                        -    -             - -  -       BCPO 71    2 1
  0 BOND_BM_OCPI   $ROOT/BCOC_BOND/v2014-07/Bond_biomass.nc OC   2000/1-12/1/0 C xy kg/m2/s OCPI 72    2 1
  0 BOND_BM_OCPO   -                                        -    -             - -  -       OCPO 73    2 1
  0 BOND_BM_POG1   -                                        -    -             - -  -       POG1 74/76 2 1
  0 BOND_BM_POG2   -                                        -    -             - -  -       POG2 74/77 2 1
  ))).not.GFED.or.FINN.or.QFED2

--Bob Yantosca (talk) 16:02, 17 May 2018 (UTC)

References

  1. Duncan, B.N., et al., Interannual and Seasonal Variability of Biomass Burning Emissions Constrained by Satellite Observations, J. Geophys. Res., 108(D2), 4040, doi:10.1029/2002JD002378, 2003. PDF
  2. Hyer, E., FLAMBE Biomass Burning emissions for ARCTAS, 2008. PDF
  3. Lobert, J. M., W. C. Keene, J. A. Logan, and R. Yevich, Global chlorine emissions from biomass burning: the reactive chlorine emissions inventory, J. Geophys. Res., 8, 2999-3014, 2008.
  4. Mu, M., J.T. Randerson, G.R. van der Werf, L. Giglio, P. Kasibhatla, D. Morton, G.J. Collatz, R.S. DeFries, E.J. Hyer, E.M. Prins, D.W.T. Griffith, D. Wunch, G.C. Toon, V. Sherlock, and P.O. Wennberg, Daily and 3-hourly variability in global fire emissions and consequences for atmospheric model predictions of carbon monoxide, Journal of Geophysical Research-Atmospheres, 116, D24303, doi:10.1029/2011JD016245, 2011.
  5. Nassar, R., J. A. Logan, I. A. Megretskaia, L. T. Murray, L. Zhang, and D. B. A. Jones, Analysis of tropical tropospheric ozone, carbon monoxide and water vapor during the 2006 El Niño using TES observations and the GEOS-Chem model, J. Geophys. Res., 114, D17304, doi:10.1029/2009JD011760, 2009. PDF
  6. Liu, T., L.J. Mickley, R.S. DeFries, M.E. Marlier, M.F. Khan, M.T. Latif, and A. Karambelas, Diagnosing spatial uncertainties and relative biases in global fire emissions inventories: Indonesia as regional case study, Remote Sens. Environ., 237, 111557. 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
  8. Wiedinmyer, C., S.K. Akagi, R.J. Yokelson, L.K. Emmons, J.J. Orlando, and A.J. Soja, The Fire INventory from NCAR (FINN): a high resolution global model to estimate the emissions from open burning, Geosci. Model Dev., 4, 625–641, 2011. Article
  9. Kaiser, J.W., A. Heil, M.O. Andreae, A. Benedetti, N. Chubarova, L. Jones, J.J. Morcrette, M. Razinger, M.G. Schultz, M. Suttie, and G.R. van der Werf, Biomass burning emissions estimated with a global fire assimilation system based on observed fire radiative power, Biogeosciences, 9, 527–554, 2012. Article
  10. Darmenov, A.S. and A. da Silva, The Quick Fire Emissions Dataset (QFED) - Documentation of versions 2.1, 2.2, and 2.4, NASA Technical Report Series on Global Modeling and Data Assimilation, Volume 32, 2013. Article
  11. Ichoku, C. and L. Ellison, Global top-down smoke-aerosol emissions estimation using satellite fire radiative power measurements, Atmos. Chem. Phys., 14, 6643–6667, 2014. Article

--Melissa Payer 15:10, 21 February 2012 (EST)