Historical and Future Emissions

From Geos-chem
Revision as of 18:50, 3 May 2013 by Chris Holmes (Talk | contribs) (RCP future emission scenarios)

Jump to: navigation, search

On this page we shall provide information about historical and future emissions inventories.

Historical emission inventories of SO2, NOx, BC, and POA


This update was tested in the 1-month benchmark simulation v9-01-03h and approved on 09 Mar 2012.

Leibensperger et al. (2011) implemented historical emission inventories of SO2, NOx, BC, and POA in order to simulate the evolution aerosol concentrations between 1950 and 2000 in the US. Global emissions are available from 1850 for BC/POA and 1890 for SO2/NOx. Emissions are from the following sources:

  • SO2/NOx between 1890-1990 - EDGAR HYDE (van Aardenne et al., 2001)
  • SO2/NOx 2000 - EDGAR 3.2 FT (Olivier and Berdowski, 2001: Global emissions sources and sinks in "The Climate System")
  • BC/POA between 1850-2000 - Bond et al. 2007

Future emissions of SO2, NOx, BC, and POA (2010 to 2050) were calculated by applying regional scale factors to emissions for year 2000. Scale factors follow the IPCC A1B scenario as implemented by Streets et al. (2004).

These emissions were implemented into the standard code in v9-01-03h. For more information, contact Eric Leibensperger.

--Melissa Payer 15:51, 9 January 2013 (EST)

Bug fix for regridding of historical emissions

In v9-01-03 there is a bug in the regridding of historical emissions in carbon_mod.F. The following lines:

           ! Cast to REAL*8 and resize
           CALL TRANSFER_2D( GENERIC_1x1(:,:,1), FD2D )

should be replaced with:

           ! Cast from REAL*4 to REAL*8
           GEN8_1x1 = GENERIC_1x1

           ! Regrid to current model resolution [kg/yr]
           INGRID  => GEN8_1x1(:,:)
           CALL DO_REGRID_A2A( LLFILENAME, I1x1, J1x1-1, 
    &                          INGRID,     FD2D, PERAREA=1,
    &                          netCDF=.TRUE.                )

           ! Free pointer
           NULLIFY( INGRID )

The above changes should be made everywhere that historical emissions are read in carbon_mod.F. The following USE statements and definitions will also need to be added at the top of subroutine ANTHRO_CARB_TBOND.


     CHARACTER(LEN=255)            :: LLFILENAME
     REAL*8, POINTER               :: INGRID (:,:) => NULL()
     REAL*8, TARGET                :: GEN8_1x1(I1x1,J1x1-1)
     REAL*4                        :: GENERIC_1x1(I1x1,J1x1-1)

This bug was fixed in v9-02e along with several other fixes for regridding.

--Melissa Payer 18:14, 14 February 2013 (EST)

RCP future emission scenarios

Chris Holmes (UC Irvine) has implemented the RCP anthropogenic emission inventories into GEOS-Chem. Please see this document for a description of the RCP emission scenarios, implementation into GEOS-Chem, and comparison to existing GEOS-Chem emission inventories.

This updated is slated for inclusion in GEOS-Chem v9-02.

--Melissa Payer 15:51, 9 January 2013 (EST)

The large difference in O3 and OH seen over the oceans when comparing RCP emissions to the prior GEOS-Chem emissions is due to a bug in the ICOADS-AMVER ship emission inventory. See Section 3.3 of this document. The ICOADS-AMVER inventory should be 5.0 Tg(N)/yr, but was actually 1.5 Tg(N)/yr in GEOS-Chem (v9-01-03 and earlier). Correcting this bug dramatically reduces the differences in simulated O3 and OH concentrations over the oceans and makes the CH4 lifetime in the RCP about 5% longer than in the base GEOS-Chem simulation.

--Chris Holmes 14:50, 3 May 2013 (EDT)