Historical and Future Emissions

From Geos-chem
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)


  1. Leibensperger, E.M., L.J. Mickley, D.J. Jacob, W.-T. Chen, J.H. Seinfeld, A. Nenes, P.J. Adams, D.G. Streets, N. Kumar, and D. Rind, Climatic effects of 1950-2050 changes in US anthropogenic aerosols - Part 1: Aerosol trends and radiative forcing, Atmos. Chem. Phys., 12, 3333-3348, 2012.

--Bob Y. 14:30, 27 February 2015 (EST)

Bug fix for regridding of historical emissions

NOTE: This issue was resolved in GEOS-Chem v9-02e. Several other regridding issues were also corrected in this version.

This code was removed in GEOS-Chem v10-01. All emissions are now handled by the HEMCO emissions component.

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)

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

RCP future emission scenarios

In GEOS-Chem v10-01 and newer versions

In GEOS-Chem v10-01, the RCP emissions are implemented via the HEMCO emissions component. For more information, please see our RCP future emission scenarios wiki page.

--Bob Y. 14:28, 27 February 2015 (EST)

In versions prior to GEOS-Chem v10-01

GEOS-Chem v9-02 and prior versions read the RCP future emissions scenarios from binary punch files. For more information about this implementation, please see our RCP future emission scenarios wiki page.

--Bob Y. 14:28, 27 February 2015 (EST)