EDGAR anthropogenic emissions

From Geos-chem
Revision as of 21:14, 30 March 2014 by Meng Li (Talk | contribs)

Jump to: navigation, search

On this page we provide information about the Emissions Database for Global Atmospheric Research (EDGAR) anthropgenic emissions inventories used in GEOS-Chem.

Anthropogenic emissions

The EDGAR emissions of NOx, CO, and SOx are the standard emissions. However, the EDGAR emissions are typically overwritten by data from various regional emissions inventories listed on our Anthropogenic emissions wiki page.

The EDGAR inventory also includes emissions from ship exhaust.

--Bob Y. 12:27, 20 December 2012 (EST)

Ship emissions

For more information about the EDGAR ship emissions data files, please see the following README document: GEOS_1x1/EDGAR_200607/README

The source code routines that read the EDGAR ship emissions data from disk are:

  1. SEASCL_EDGAR_SHIP_SO2 (GeosCore/edgar_mod.F)
  2. GET_EDGAR_SHIP_SO2 (GeosCore/edgar_mod.f)

The EDGAR ship SO2 emissions are reported as per year. Routine SEASCL_EDGAR_SHIP_SO2 applies monthly scale factors (which are computed as the ratio of monthly/total ship SO2 emissions from the Corbett et al inventory to the annual EDGAR ship SO2 emissions. This is required to impose aseasonality onto the EDGAR ship SO2 emissions.

--Bob Y. 12:27, 20 December 2012 (EST)

References

  1. Olivier, J.G.J. and J.J.M. Berdowski, Global emissions sources and sinks. In: Berdowski, J., Guicherit, R. and B.J. Heij (eds.) The Climate System, pp. 33-78. A. A. Balkema Publishers/Swets & Zeitlinger Publishers, Lisse, The Netherlands., 2001
  2. A related document can be found here: EDGAR description (pdf)
  3. Aaron van Donkelaar has created a PDF document that shows some of the effects of using the EDGAR emissions.

--Bob Y. 12:31, 20 December 2012 (EST)

Use EDGAR v.4.1 in GEOS-Chem

EDGAR v4.1 emissions of CO, NOx and SO2 from 1970-2005 were prepared for GEOS-Chem model. It is available on a yearly basis with 1x1 degree spatial resolution. This inventory is not included in standard model yet but can be accessed here. It has been tested in GEOS-Chem v09-01-01 with 1-year full-chemistry benchmark for 2005.

To use EDGAR v.4.1 inventory, please download and the patch codes and data from the following URLs:

  1. http://166.111.42.46/EDGARV4.1_201111.zip
  2. http://166.111.42.46/Patch_EDGAR4_120702.zip

NOTE: for edgar_mod.f in the patch folder, please comment out line 1172 and 1173 in order to prevent double counting of international aviation for CO.

The patch folder contains 6 files: edgar_mod.f, biofuel_mod.f, input_mod.f, logical_mod.f, sulfate_mod.f, and input.geos. In input.geos, USE_EDGAR4 is a logical variable for user to choose EDGAR v.4.1 as anthropogenic emission inventory. Since biofuel emissions have been included in EDGAR v.4.1, if use EDGAR v.4.1, default biofuel emissions in GC model will be zero out automatically. If you would like to use shipping emissions from EDGAR v.4.1, also set LEDGARSHIP to True. Aircraft emissions in EDGAR v.4.1 were not implemented for NOx and SO2.


EDGAR v.4.1 emission data includes the following sectors.

EDGAR_V41_IPCC_ 1A1_ENE
Energy (fuel combustion).
EDGAR_V41_IPCC _1A2_2
Industry combustion and process emissions (fuel combustion).
EDGAR_V41_IPCC_1A3b_c_e
Ground Transport (road transport, rail transport and other transport, fuel combustion).
EDGAR_V41_IPCC_1A3d1
International shipping (fuel combustion).
EDGAR_V41_IPCC_1A3d_SHIP
Inland Waterways (fuel combustion).
EDGAR_V41_IPCC_1A4_5
Residential and other combustion (fuel combustion).
EDGAR_V41_IPCC_1B
Fuel production.
EDGAR_V41_IPCC_6A_6C
Solid Waste.
EDGAR_V41_IPCC_7
Other Sources.

Original documentation of EDGAR inventory can be found from the following URL: http://edgar.jrc.ec.europa.eu/overview.php?v=41

If you have any question, please contact Qiang Zhang (qiangzhang@tsinghua.edu.cn) at Tsinghua University .

--Bob Y. 11:37, 20 December 2012 (EST)

Implementation of EDGAR 4.2 into GEOS-Chem

EDGAR 4.2 emissions of NOx, CO, SO2, and NH3 from 1970-2008 were prepared and implemented into GEOS-Chem. Currently the implementation is based on version 9-01-02 and it will be done for version 9-02. This global emission inventory is added as a choice for users to use in GEOS-Chem simulation.

To use EDGAR 4.2 inventory, please download and the patch codes and data from the following URL: http://rain.ucis.dal.ca/upload/mli/

Please refer to the documentation for more details: http://rain.ucis.dal.ca/upload/mli/Documentation_EDGARinGEOSChem_20140328.pdf.

Previous issues that are now resolved

Potential bug in EDGAR ship emissions code

This update was tested in the 1-month benchmark simulation v9-01-02s and approved on 5 Nov 2011.

Michael Long discovered this potential bug in the EDGAR emissions code.

The current GEOS-Chem default configuration uses the ICOADS ship emissions instead of the EDGAR ship emissions. In routine READ_EMISSIONS_MENU in input_mod.f, we have a safety check...we turn off the logical switch LEDGARSHIP if switch LICOADS is turned on:

     !=================================================================
     ! Check global ship emissions options
     !=================================================================
     ! Add an ship emissions options (cklee, 6/30/09)
     ! Replace with ICOADS ship emissions
     ! User can still use regional EMEP, and ARCTAS-SO2 if he wants (phs) 
     IF ( LICOADSSHIP ) THEN
        LEDGARSHIP = .FALSE. ! ICOADS precedes EDGAR
        LSHIPSO2   = .FALSE. ! SO2-ICOADS precedes SO2-Corbett
     ENDIF    

In edgar_mod.F, the EDGAR ship emissions arrays are not allocated if the switch LEDGARSHIP == .FALSE.:

     IF ( LEDGARSHIP ) THEN

        ALLOCATE( EDGAR_SO2_SHIP( IIPAR, JJPAR ), STAT=AS )
        IF ( AS /= 0 ) CALL ALLOC_ERR( 'EDGAR_SO2_SHIP' )
        EDGAR_SO2_SHIP = 0d0

        ALLOCATE( EDGAR_NOx_SHIP( IIPAR, JJPAR ), STAT=AS )
        IF ( AS /= 0 ) CALL ALLOC_ERR( 'EDGAR_NOx_SHIP' )
        EDGAR_NOx_SHIP = 0d0

        ALLOCATE( EDGAR_CO_SHIP( IIPAR, JJPAR ), STAT=AS )
        IF ( AS /= 0 ) CALL ALLOC_ERR( 'EDGAR_CO_SHIP' )
        EDGAR_CO_SHIP = 0d0
 
     ENDIF

However, we still pass EDGAR_CO_SHIP to routine COMPUTE_EDGAR_CO even if it is not allocated:

        IF ( LEDGARCO ) THEN
           CALL COMPUTE_EDGAR_CO( YEAR, EDGAR_CO, EDGAR_CO_SHIP )
        ENDIF

Here, EDGAR_CO_SHIP is an output variable, which is supposed to return the EDGAR ship emissions at the current model resolution. This may produce undefined behavior. This may trigger an error when compiling with the check-array-bounds option (BOUNDS=yes).

Some of the more forgiving compilers may "fix" this on-the-fly and thus, you might not see any errors. However, for the sake of good programming practice, this should be fixed.

--Bob Y. 12:30, 20 December 2012 (EST)

Typos in edgar_mod.f

NOTE: This issue was resolved in GEOS-Chem v9-01-01.

In routines READ_C2H2 and READ_C2H4 of edgar_mod.f, the data directory was wrongly listed as EDGAR200607 instead of EDGAR_200607. We have now fixed these typos in v9-01-01.

--Bob Y. 12:35, 20 December 2012 (EST)