CH4 simulation

NOTE: Page under construction

This page contains information about the methane (CH4) simulation in GEOS-Chem.

Overview

The CH4 simulation in GEOS-Chem ...

Authors and collaborators

• Christopher Pickett-Heaps (Harvard)
• Kevin Wecht (Harvard)
• Jerome Drevet (formerly @ EPFL)
• James Wang (formerly @ Harvard)

CH4 simulation user groups

User Group	Personnel	Projects
Harvard University	Christopher Pickett-Heaps Kevin Wecht	Methane sources in the Arctic
U. Wollongong	Chris Chan Miller	Methane sources in tropical Australia

Planned updates

Update methane emissions and add capability for North America nested grid

Note: This update is slated to be incorporated in GEOS-Chem v9-01-03 or later.

Kevin Wecht wrote:

I have some GEOS-Chem methane updates that I'd like to add to the standard code.

1. GFED emissions - I've updated GFED-2 and GFED-3 biomass burning to include methane emissions.
2. Methane simulation - I've updated emissions and error catches to make the simulation compatible with the nested simulation over North America. I added linearized CH4 chemistry from GMI output (CH4 loss frequencies generated by Lee Murray). These changes should be made to the data directory:

a. Replace Carbon_soil, Carbon_litter, and Wetfrac files in data directory (for all horizontal resolutions)

b. Add 2010-2011 TSKIN and GWETTOP files to data directory (for all horizontal resolutions)

c. Add nested north american OH files to OH directory

d. Add biofuel emissions to data directory (for all horizontal resolutions)

e. Add all emission files to data directory for nested North America simulation

--Melissa Payer 17:35, 2 February 2012 (EST)

Known issues

Error reading GWETTOP fields

Melissa Payer discovered a file read error in the current CH4 simulation. We will fix this issue in GEOS-Chem v9-01-03.

Melissa Payer wrote:

I got the offline CH4 simulation rice emissions to work. It turns out that XTAU is being reset [to a the wrong value] when GWETTOP is read. To fix this, I added the following code in RICE_EMIS prior to setting the soil wetness FILENAME:

     ! Determine year for soil wetness files (mpayer, 1/31/12)
     IF ( GET_YEAR() < 2004 ) THEN
        CYEAR  = '2004'
        CYEAR2 = 2004
     ELSEIF (GET_YEAR() > 2009 ) THEN
        CYEAR  = '2009'
        CYEAR2 = 2009
     ELSE
        WRITE( CYEAR, '(i4)' ) GET_YEAR()
        CYEAR2 = GET_YEAR()
     ENDIF

I also changed:

     XTAU = GET_TAU0( 1, 1, GET_YEAR() )

to:

     XTAU = GET_TAU0( 1, 1, CYEAR2 )

For both annual and monthly mean GWETTOP. That seems to work and it reads the file for the proper simulation year.

--Bob Y. 10:29, 1 February 2012 (EST)