|
|
(37 intermediate revisions by 2 users not shown) |
Line 1: |
Line 1: |
− | This page contains information about the methane (CH4) simulation in GEOS-Chem.
| + | __FORCETOC__ |
| + | '''''[[Aerosol-only simulation|Previous]] | [[CO2 simulation|Next]] | [[Guide to GEOS-Chem simulations]]''''' |
| + | #[[Simulations using KPP-built mechanisms|Simulations using KPP-built mechanisms (carbon, fullchem, Hg)]] |
| + | #[[Aerosol-only simulation]] |
| + | #<span style="color:blue">'''CH4 simulation'''</span> |
| + | #[[CO2 simulation]] |
| + | #[[Metals simulation]] |
| + | #[[POPs simulation]] |
| + | #[[Tagged CO simulation]] |
| + | #[[Tagged O3 simulation]] |
| + | #[[TransportTracers simulation]] |
| + | |
| | | |
| == Overview == | | == Overview == |
− | | + | This page contains information about the methane (CH4) simulation in GEOS-Chem. |
− | === Authors and collaborators ===
| + | |
− | * [http://jdmaasakkers.github.io/ Bram Maasakkers] ''(Harvard)''
| + | |
− | * [[User:Melissa Payer|Melissa Sulprizio]] ''(Harvard)''
| + | |
− | * [https://scholar.harvard.edu/jsheng Jianxiong Sheng] ''(Harvard)''
| + | |
− | * [http://www.prism.gatech.edu/~yzhang425/ Yuzhong Zhang] ''(Harvard)''
| + | |
− | * [https://eps.harvard.edu/people/tia-scarpelli Tia Scarpelli] ''(Harvard)''
| + | |
− | * [https://alexjturner.github.io Alex Turner] ''(UC Berkeley)''
| + | |
− | | + | |
− | Previous authors include:
| + | |
− | * Kevin Wecht ''(formerly @ Harvard)''
| + | |
− | * Christopher Pickett-Heaps ''(formerly @ Harvard)''
| + | |
− | * Jerome Drevet ''(formerly @ EPFL)''
| + | |
− | * James Wang ''(formerly @ Harvard)''
| + | |
| | | |
| === CH4 simulation user groups === | | === CH4 simulation user groups === |
| | | |
| See [[Carbon_Cycle_Working_Group#Current_Carbon_Cycle_Projects_.28please_add_yours.21.29|this post on the ''Carbon Cycle Working Group'' wiki page]]. | | See [[Carbon_Cycle_Working_Group#Current_Carbon_Cycle_Projects_.28please_add_yours.21.29|this post on the ''Carbon Cycle Working Group'' wiki page]]. |
− |
| |
− | == Recent science updates ==
| |
− |
| |
− | The following updates were recently added into the GEOS-Chem CH4 simulation:
| |
− |
| |
− | === Updated CH4 simulation in v11-02 and later ===
| |
− |
| |
− | <span style="color:darkorange">'''''This update is slated for inclusion in [[GEOS-Chem v11-02]].'''''</span>
| |
− |
| |
− | Bram Maasakkers and Melissa Sulprizio have updated the CH4 simulation in [[GEOS-Chem v11-02]]. These updates include the following emission updates (brought in via HEMCO):
| |
− |
| |
− | #[http://edgar.jrc.ec.europa.eu/ EDGAR v4.3.2] global anthropogenic emissions
| |
− | #Gridded US EPA emissions ([http://pubs.acs.org/doi/pdf/10.1021/acs.est.6b02878 Maasakkers et al. 2016])
| |
− | #Canada and Mexico oil & gas emissions ([http://www.sciencedirect.com/science/article/pii/S1352231017301164 Sheng et al. 2017])
| |
− | #Daily [[QFED biomass burning emissions]]
| |
− | #JPL WetCHARTs emissions ([https://www.geosci-model-dev.net/10/2141/2017/gmd-10-2141-2017.pdf Bloom et al. 2017])
| |
− | #Global geological seep emissions (Maasakkers et al., in prep)
| |
− | #Global lake emissions (Maasakkers et al., in prep)
| |
− |
| |
− | and code updates:
| |
− |
| |
− | #CH4 loss by Cl (Maasakkers et al., in prep)
| |
− | #GOSAT and TCCON observation operators (modified from the adjoint code)
| |
− |
| |
− | The following data files are required for these updates
| |
− |
| |
− | ftp://ftp.as.harvard.edu/gcgrid/data/ExtData/HEMCO/CH4/v2017-10/ (entire directory)
| |
− |
| |
− | ftp://ftp.as.harvard.edu/gcgrid/data/ExtData/HEMCO/MASKS/v2014-07/CONUS_Mask.001x001.nc
| |
− |
| |
− | ftp://ftp.as.harvard.edu/gcgrid/data/ExtData/HEMCO/MASKS/v2014-07/Canada_Mask.Extended_NAgrid.001x001.nc
| |
− |
| |
− | ftp://ftp.as.harvard.edu/gcgrid/data/ExtData/HEMCO/MASKS/v2014-07/Mexico_Mask.Extended_NAgrid.001x001.nc
| |
− |
| |
− | A few caveats:
| |
− | *EDGAR v4.3.2 lumps oil, gas, and coal emissions in the files available on their website. Ideally, we would like to separate these sources, so we are looking into that. In addition, Tia Scarpelli at Harvard is developing a global oil and gas CH4 emissions inventory which will be included in the future. Until we resolve this, the EDGAR v4.3.2 lumped oil/gas/coal emissions are added to CH4_OIL in the tagged CH4 simulation and CH4_GAS, CH4_COL only include emissions from Canada, US, and Mexico.
| |
− | *EDGAR v4.3.2 rice emissions by default are annual emissions. For our CH4 simulations, we generated monthly EDGAR rice emissions by applying seasonality from the [http://onlinelibrary.wiley.com/doi/10.1002/2016GB005381/abstract Zhang et al. (2016)]. We had to preprocess the EDGAR files to do this and haven't yet figured out a way to apply the seasonality via HEMCO -- it may require writing an extension. If you're interested in using our EDGAR seasonal rice emissions, please let us know and we can send you those files.
| |
− | *We have updated rice emissions from Zhang et al. (2016), but we don't have permission to share those data yet. We used those rice emissions in Maasakkers et al. (in prep) to apply seasonality to the EDGAR v4.3.2 rice emissions.
| |
− | *We recommend not use the lakes emissions (they are turned off by default in HEMCO_Config.rc). There is some concern it leads to double-counting when used together with WetCHARTS. It’s something we’ll look into but just using WetCHARTS for now is a cleaner approach.
| |
− | *There’s a small issue with incorporating offshore GEPA emissions right now (as those fall outside of the US mask. We’re working on resolving that.
| |
− |
| |
− | --[[User:Melissa Payer|Melissa Sulprizio]] ([[User talk:Melissa Payer|talk]]) 17:34, 16 November 2017 (UTC)
| |
− |
| |
− | === Update CH4 emissions to EDGAR v4.2 ===
| |
− |
| |
− | The CH4 simulation now uses emissions from the EDGAR v4.2 inventory.
| |
− |
| |
− | ==== In GEOS-Chem v10-01 and newer versions ====
| |
− |
| |
− | In [[GEOS-Chem v10-01]] and newer versions, EDGAR v4.2 CH4 emissions are read with the [[HEMCO|HEMCO emissions component]]. We have created [http://ftp.as.harvard.edu/gcgrid/data/ExtData/HEMCO/CH4/v2014-09/README new emission data files for the CH4 simulation] (in [[Preparing data files for use with HEMCO|COARDS-compliant netCDF format]]) for use with HEMCO. These new data files are contained in the [[HEMCO data directories|HEMCO data directory tree]]. For detailed instructions on how to download these data files to your disk server, please see our [[HEMCO_data_directories#Downloading_the_HEMCO_data_directories|''Downloading the HEMCO data directories'' wiki post]].
| |
− |
| |
− | --[[User:Bmy|Bob Y.]] 13:23, 3 March 2015 (EST)
| |
− |
| |
− | ==== In GEOS-Chem versions prior to v10-01 ====
| |
− |
| |
− | <span style="color:green">'''''This update was tested in the 1-month benchmark simulation [[GEOS-Chem_v9-02_benchmark_history#v9-02m|v9-02m]] and approved on 30 Jul 2013. This update is included in Adjoint [[GEOS-Chem_Adjoint_v35 | v35d]].'''''</span>
| |
− |
| |
− | Kevin Wecht has prepared EDGAR v4.2 emissions of methane for 2004-2008. The data files are available at:
| |
− |
| |
− | ftp://ftp.as.harvard.edu/pub/geos-chem/data/GEOS_0.5x0.666_NA/CH4_201305/
| |
− |
| |
− | ftp://ftp.as.harvard.edu/pub/geos-chem/data/GEOS_2x2.5/CH4_201305/
| |
− |
| |
− | ftp://ftp.as.harvard.edu/pub/geos-chem/data/GEOS_4x5/CH4_201305/
| |
− |
| |
− | You can download these directories with anonymous FTP or the Wget utility. For instructions, please see [http://acmg.seas.harvard.edu/geos/doc/man/chapter_2.html#2.4 Chapter 2.4, Downloading the GEOS-Chem shared data directories] in the GEOS-Chem Online User's Guide.
| |
− |
| |
− | Modifications to subroutines <tt>RICE_EMIS</tt> and <tt>ASEASONAL_ANTHRO_EMIS</tt> in <tt>global_ch4_mod.F</tt> include changing the directory names and updating CYEAR and XTAU as follows:
| |
− |
| |
− | WRITE( CYEAR, '(i4)' ) GET_YEAR()
| |
− | XTAU = GET_TAU0( 1, 1, GET_YEAR() )
| |
− |
| |
− | IF ( GET_YEAR() .LT. 2004 ) THEN
| |
− | CYEAR='2004'
| |
− | XTAU = GET_TAU0( 1, 1, 2004 )
| |
− | ENDIF
| |
− | IF ( GET_YEAR() .GT. 2008 ) THEN
| |
− | CYEAR='2008'
| |
− | XTAU = GET_TAU0( 1, 1, 2008 )
| |
− | ENDIF
| |
− |
| |
− | --[[User:Melissa Payer|Melissa Payer]] 16:53, 2 July 2013 (EDT)
| |
− |
| |
− | == Validation ==
| |
| | | |
| == References == | | == References == |
Line 118: |
Line 32: |
| #Wang, J.S., J.A. Logan, M.B. McElroy, B.N. Duncan, I.A. Megretskaia, and R.M. Yantosca, ''A 3-D model analysis of the slowdown and interannual variability in the methane growth rate from 1988 to 1997'', <u>Global Biogeochem. Cycles</u>, '''18''', GB3011, doi:10.1029/2003GB002180, 2004. | | #Wang, J.S., J.A. Logan, M.B. McElroy, B.N. Duncan, I.A. Megretskaia, and R.M. Yantosca, ''A 3-D model analysis of the slowdown and interannual variability in the methane growth rate from 1988 to 1997'', <u>Global Biogeochem. Cycles</u>, '''18''', GB3011, doi:10.1029/2003GB002180, 2004. |
| | | |
− | == Previous issues that are now resolved ==
| + | ---- |
− | | + | '''''[[Aerosol-only simulation|Previous]] | [[CO2 simulation|Next]] | [[Guide to GEOS-Chem simulations]]''''' |
− | The following bugs and/or issues with the code have now been resolved:
| + | |
− | | + | |
− | === Now call INIT_GLOBAL_CH4 from input_mod.F ===
| + | |
− | | + | |
− | <span style="color:green">'''''This update was validated in the 1-month benchmark simulation [[GEOS-Chem v10-01 benchmark history#v10-01c|v10-01c]] and approved on 29 May 2014.'''''</span>
| + | |
− | | + | |
− | Prior to [[GEOS-Chem v10-01]], routine <tt>INIT_GLOBAL_CH4</tt> (in module <tt>GeosCore/global_ch4_mod.F</tt>) was called from the main driver routine <tt>main.F</tt> at the start of the simulation.
| + | |
− | | + | |
− | In order to make the code more compatible with our [[Grid-Independent GEOS-Chem]] development, we now call <tt>INIT_TAGGED_CH4</tt> from routine <tt>GIGC_INIT_EXTRA</tt> (in module <tt>GeosCore/input_mod.F</tt>).
| + | |
− | | + | |
− | --[[User:Bmy|Bob Y.]] 16:59, 30 May 2014 (EDT)
| + | |
− | | + | |
− | === Minor fixes to CH4 simulation for MERRA meteorology ===
| + | |
− | | + | |
− | <span style="color:green">'''''These updates were validated in the 1-month benchmark simulation [[GEOS-Chem v10-01 benchmark history#v10-01c|v10-01c]] and approved on 29 May 2014.'''''</span>
| + | |
− | | + | |
− | We have made a couple of fixes to the CH4 simulation module (<tt>GeosCore/global_ch4_mod.F</tt>) for compatibility with the [[MERRA]] met fields. At lines 1114, and again at line 1503, there were #if blocks that looked like this:
| + | |
− | | + | |
− | #if defined( GEOS_5 ) || defined( GEOS_FP )
| + | |
− | | + | |
− | but these should have been:
| + | |
− | | + | |
− | #if defined( GEOS_5 ) || defined( GEOS_FP ) || defined( MERRA )
| + | |
− | | + | |
− | Furthermore, in data directory <tt>GEOS_4x5/CH4_201305/wetlands</tt>, we have symbolically linked files:
| + | |
− | | + | |
− | GWETTOP.geos5.4x5.YYYY.bpch --> GWETTOP.merra.4x5.YYYY.bpch
| + | |
− | TSKIN.geos5.4x5.YYYY.bpch --> TSKIN.merra.4x5.YYYY.bpch
| + | |
− | | + | |
− | where <tt>YYYY = { 2004, 2005, 2006, 2007, 2008, 2009, 2010 }</tt>.
| + | |
− | | + | |
− | These issues were discovered with the [[GEOS-Chem Unit Tester]].
| + | |
− | | + | |
− | --[[User:Bmy|Bob Y.]] 16:59, 30 May 2014 (EDT)
| + | |
− | | + | |
− | == Unresolved issues ==
| + | |
− | | + | |
− | None reported so far.
| + | |
− | | + | |
− | --[[User:Bmy|Bob Y.]] 14:07, 30 March 2012 (EDT)
| + | |
− | | + | |
− | == Obsolete information ==
| + | |
− | | + | |
− | These sections pertain to code that has been removed from the most recent versions of GEOS-Chem. We shall leave this information here for reference.
| + | |
− | | + | |
− | === Notes about the methane simulation ===
| + | |
− | | + | |
− | [[Image:Obsolete.jpg]]
| + | |
− | | + | |
− | <span style="color:red">'''''In [[GEOS-Chem v10-01]] and newer versions, CH4 emissions are implemented via the [[HEMCO|HEMCO emissions component]].'''''</span>
| + | |
− | | + | |
− | :1. '''annual_met.pro:''' The methane simulation uses monthly and annual average skin temperature and soil moisture values from the met fields. These are calculated offline using IDL ("annual_met.pro"). This file can be found on the ftp site in the methane data directory with the LPJ wetland emission scheme (e.g., "GEOS_0.5x0.666_NA/CH4_201305/wetlands"). The script should work as is for 2004-2011, however we don't currently have the met data for 2012 so you may need to modify "global_ch4_mod.F" if you're interested in 2012 or 2013.
| + | |
− | ::: '''Subroutine "WETLAND_EMIS"
| + | |
− | YEAR = GET_YEAR()
| + | |
− | #if defined( GEOS_5 )
| + | |
− | IF ( YEAR .LT. 2004 ) YEAR = 2004
| + | |
− | IF ( YEAR .GT. 2011 ) YEAR = 2011
| + | |
− | | + | |
− | ::: '''Subroutine "RICE_EMIS"
| + | |
− | ! Get annual and monthly mean soil wetness from GEOS
| + | |
− | ! One file contains both monthly and annual mean GWETTOP
| + | |
− | YEAR = GET_YEAR()
| + | |
− | #if defined( GEOS_5 )
| + | |
− | IF ( YEAR .LT. 2004 ) YEAR = 2004
| + | |
− | IF ( YEAR .GT. 2011 ) YEAR = 2011
| + | |
− | | + | |
− | :2. '''Seasonality of Rice Emissions:''' Seasonality in the rice emissions is imposed through annual and monthly average soil moisture values. This works reasonably well for major rice growing areas like South/East Asia where rice paddies are rain-fed, but it fails in other regions. For example, California rice paddies are irrigated and crops are grown in the summer, precisely when the natural soil is the driest. You can enter the "global_ch4_mod.F", subroutine "RICE_EMIS", and impose any seasonality you want on the annual emission rates from EDGAR.
| + | |
− | | + | |
− | :3. '''Potential Global Mean Bias:''' The methane simulation may not simulate global average concentrations accurately. That is, you may have a mean bias of up to +/-30ppbv when you compare the simulation to observations. A global mean bias means that the model has a slight imbalance between the sources and sinks relative to the world. You may want to consider removing this mean bias before analyzing your model concentrations.
| + | |
− | | + | |
− | --[[User:Alex Turner|Alex Turner]] 11:01, 24 January 2014 (EST)
| + | |
− | | + | |
− | === Fixes for wetland emissions ===
| + | |
− | | + | |
− | [[Image:Obsolete.jpg]] | + | |
− | | + | |
− | <span style="color:red">'''''In [[GEOS-Chem v10-01]] and newer versions, CH4 emissions are implemented via the [[HEMCO|HEMCO emissions component]].'''''</span>
| + | |
− | | + | |
− | <span style="color:green">'''''These issues were resolved during the [[GEOS-Chem_v9-02#Period_of_public_comment|public comment period]] following the provisional release of GEOS-Chem v9-02 and were included in the official release of [[GEOS-Chem v9-02]] (03 Mar 2014). This update is included in Adjoint [[GEOS-Chem_Adjoint_v35 | v35j]].'''''</span>
| + | |
− | | + | |
− | '''''Kevin Wecht''''' wrote:
| + | |
− | | + | |
− | :I have some updates/bug fixes to wetland emissions in the methane simulation. Both changes impact global_ch4_mod.F, subroutine WETLAND_EMIS.
| + | |
− | | + | |
− | ::1. Christopher Pickett-Heaps added a feature that zeroes wetland emissions if snow covers the ground. This was never incorporated into the standard code. It is included in the attached file.
| + | |
− | | + | |
− | ::2. Parameters in the subroutine are not correct. The following values are the correct and are included in the attached file:
| + | |
− | MOIST_SCALE = 0.205
| + | |
− | EMIT_FACT = 0.018
| + | |
− | | + | |
− | --[[User:Melissa Payer|Melissa Sulprizio]] 16:11, 23 January 2014 (EST)<br>
| + | |
− | --[[User:Bmy|Bob Y.]] 10:25, 4 April 2014 (EDT)
| + | |
− | | + | |
− | === Error reading GWETTOP fields ===
| + | |
− | | + | |
− | [[Image:Obsolete.jpg]]
| + | |
− | | + | |
− | <span style="color:red">'''''In [[GEOS-Chem v10-01]] and newer versions, CH4 emissions are implemented via the [[HEMCO|HEMCO emissions component]].'''''</span>
| + | |
− | | + | |
− | <span style="color:green">'''''This update was tested in the 1-month benchmark simulation [[GEOS-Chem_v9-01-03_benchmark_history#v9-01-03h|v9-01-03h]] and approved on 09 Mar 2012.'''''</span>
| + | |
− | | + | |
− | There is a file read error in the current CH4 simulation. We will fix this issue in [[GEOS-Chem v9-01-03]].
| + | |
− | | + | |
− | In RICE_EMIS (global_ch4_mod.F), XTAU is reset to the wrong value when GEWTTOP is read. Kevin Wecht provided a fix that accounts for different availabilities of GEOS-4 and GEOS-5 met fields. Prior to setting FILENAME for GWETTOP , add the following:
| + | |
− | | + | |
− | YEAR = GET_YEAR()
| + | |
− | #if defined( GEOS_5 )
| + | |
− | IF ( YEAR .LT. 2004 ) YEAR = 2004
| + | |
− | IF ( YEAR .GT. 2009 ) YEAR = 2009
| + | |
− | #elif defined( GEOS_4 )
| + | |
− | IF ( YEAR .LT. 2000 ) YEAR = 2000
| + | |
− | IF ( YEAR .GT. 2006 ) YEAR = 2006
| + | |
− | #endif
| + | |
− | WRITE( CYEAR, '(i4)' ) YEAR
| + | |
− | | + | |
− | :Also change:
| + | |
− | | + | |
− | XTAU = GET_TAU0( 1, 1, GET_YEAR() )
| + | |
− | and
| + | |
− | XTAU = GET_TAU0( GET_MONTH(), 1, GET_YEAR() )
| + | |
− | | + | |
− | :to:
| + | |
− | | + | |
− | XTAU = GET_TAU0( 1, 1, YEAR )
| + | |
− | and
| + | |
− | XTAU = GET_TAU0( GET_MONTH(), 1, YEAR )
| + | |
− | | + | |
− | --[[User:Melissa Payer|Melissa Payer]] 18:06, 2 February 2012 (EST)
| + | |
This page contains information about the methane (CH4) simulation in GEOS-Chem.