Difference between revisions of "Biomass burning emissions"

Revision as of 19:32, 10 February 2010

GFED2

GFED2 monthly data thru 2008

As of July 2009, the GFED2 monthly data is now available for years 1997-2008. The new 2008 data is available in bpch format at:

ftp ftp.as.harvard.edu
cd pub/geos-chem/data/GEOS_1x1/GFED2_200601/2008

You will need to make a small modification in routine GFED2_COMPUTE_BIOMASS (in gfed2_biomass_mod.f) 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 gfed2_biomass_mod.f from:

ftp ftp.as.harvard.edu
cd pub/geos-chem/patches/v8-02-02
get gfed2_biomass_mod.f

--Bob Y. 12:44, 13 July 2009 (EDT)

04-Apr-2008

Here is an update on the status of the GFED2 8-day emissions:

Ray Nassar (ray.nassar@utoronto.ca) 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} ... 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.

NOTE: GFED2 8-day emissions are now included in the standard model as a user option as of GEOS-Chem v8-02-01.

--Bob Y. 09:52, 9 July 2009 (EDT)

FLAMBE biomass emissions

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.

Data issues

Jenny Fisher (jafisher@fas.harvard.edu) wrote:

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.

To get a handle on how we are processing the file, I have been looking in Philippe's directory ~phs/IDL/dvpt/flambe/

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.

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...

Philippe Le Sager (plesager@seas.harvard.edu) wrote:

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.

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:

   /home/phs/IDL/dvpt/flambe/data/new/

You can also still get the data at the following website:

   http://www.nrlmry.navy.mil/aerosol_web/arctas_flambe/data_hourly/

Finally here are some totals we check w/ Maria and Jingqiu:

   Data for one set of satellites was still given Kg.  Now fixed, I have the following
   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

--Bob Y. 15:07, 18 February 2009 (EST)

Obtaining the FLAMBE data

The FLAMBE files are available at: http://www.nrlmry.navy.mil/aerosol_web/arctas_flambe/data_hourly/.

Bob Yantosca has some scripts that can be used to download the data in his NRT_ARCTAS/perl subdirectory. These are:

sleepFlambe

Perl script which issues a Unix wget command to download one day (24 hourly files) of FLAMBE data to disk.

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 GEOS-Chem bpch format.

The FLAMBE data are kept on the Harvard data archive in the directories:

ftp ftp.as.harvard.edu
cd /pub/geos-chem/NRT-ARCTAS/flambe/YYYY/MM

where YYYY/MM are the year and month of the data. Currently we have the full data archive from 2008 and 2009, with some months of 2010.

--Bob Y. 14:32, 10 February 2010 (EST)

Duncan et al "Seasonal" Biomass Emissions

Duncan et al "Interannual" Biomass Emissions