Implementation of RETRO Anthropogenic Emissions
From Geos-chem
On this page we shall provide information about the RETRO anthropogenic VOC emissions inventory, and its implementation into GEOS-Chem.
Contents |
Overview
Dear colleagues:
Updating the anthropogenic VOC emissions was identified as a priority update for GEOS-Chem at the previous GEOS-Chem Scientific and Users' Meeting. We have recently completed an update which allows RETRO anthropogenic VOC emissions to be read into the model. The RETRO inventory has improved temporal and spatial resolution and a more inclusive suite of speciated NMVOCs than many of the inventories currently employed. The new code allows for input of 12 anthropogenic VOCs, with a notable exception of acetone.
A description of RETRO and a comparison with the existing inventories is posted here. We have submitted this update to the GEOS-Chem support team for inclusion in the mainline model. It will remain an optional inventory (rather than the default) until there is some opportunity for comparison against observations.
Best regards,
Wes Reinhart
University of Minnesota
Implementation details
Here is a timeline of the implementation of RETRO emissions into GEOS-Chem.
| Date | Update |
|---|---|
| 29 Mar 2011 | The RETRO code was ported from v8-02-01 to v8-03-01 with minor changes. Please refer to the above documentation for details. |
| 20 May 2011 | Updated documentation with full-chemistry simulation results and several minor edits. |
| 17 Aug 2011 | RETRO emissions added into standard GEOS-Chem code and tested with 1-month benchmark simulation v9-01-02k. Benchmark was approved by Dylan Millet and Wes Reinhart. |
| 22 Aug 2011 | Two minor errors in the RETRO implementation were discovered by Dylan Millet. This update was tested in the 1-month benchmark simulation v9-01-02o and approved on 27 Sep 2011. |
--Bob Y. 16:21, 22 August 2011 (EDT)
References
- Implementation of RETRO Anthropogenic Emissions into GEOS-Chem, Wes Reinhart @ U. Minnesota
- Bolshcer, M, et al., RETRO Deliverable D1-6, RETRO Documentation, 2007.
- Pulles, T, et al., Assessment of Global Emissions from Fuel Combustion in the Final Decades of the 20th Century, TNO report A-R0132/B, 2007.
--Bob Y. 16:54, 24 March 2011 (EDT)
Previous issues that are now resolved
Double-counting of residential biofuel
This update was tested in the 1-month benchmark simulation v9-01-02o and approved on 27 Sep 2011.
Dylan Millet wrote:
- I know this timing is bad, but I found a couple issues to be corrected with the retro code.
- The retro anthropogenic emissions include residential biofuel, so we need to shut off biofuel for the relevant VOCs when using retro.
- Two small bug corrections in retro_mod.f
- For biofuel_mod.f, I implemented a correction that should work in biofuel_mod.f, starting from the 9-1-1 version of the module. In routine BIOFUEL_BURN: we have to make these changes:
USE LOGICAL_MOD, ONLY : LRETRO ! (dbm, 8/18/2011)
- and
! Biofuel emissions in [molec/cm3/s]
BIOFUEL(N,I,J) = BIOFUEL_KG(N,I,J) *
& ( 6.023d23 / MOLWT(N) ) /
& ( 365d0 * 86400d0 * BOXVL(I,J,1) )
!-----------------------------------------------------------
! RETRO anthropogenic emissions include residential biofuel,
! so set this to zero if we are using RETRO for anthro VOCs
! (dbm, 8/18/2011)
!-----------------------------------------------------------
IF ( LRETRO ) THEN
IF ( (NN == IDTC2H2) .or. (NN == IDTC2H4) .or.
& (NN == IDTC2H6) .or. (NN == IDTPRPE) .or.
& (NN == IDTC3H8) .or. (NN == IDTALK4) .or.
& (NN == IDTCH2O) .or. (NN == IDTALD2) .or.
& (NN == IDTACET) .or. (NN == IDTMEK ) .or.
& (NN == IDTBENZ) .or. (NN == IDTTOLU) .or.
& (NN == IDTXYLE) ) THEN
BIOFUEL(N,I,J) = 0.d0
ENDIF
ENDIF
- We also have to make two small changes in module retro_mod.F. In routine EMISS_RETRO we need to comment out this code:
!------------------------------------------------------------------------------ ! Prior to 8/18/11: ! This shouldn't be here (dbm, 8/18/2011) ! ! Set base year to 2000 ! BASEYEAR = 2000 ! ! CALL GET_ANNUAL_SCALAR( 71, BASEYEAR, SCALEYEAR, SC) ! ALK4(:,:) = ALK4(:,:) * SC(:,:) !-----------------------------------------------------------------------------
- And we need to make this change in routine GET_RETRO_ANTHRO:
IF ( N == IDTALK4 ) THEN
RETRO = RETRO_ALK4(I,J)
!------------------------------------------------------------------------------
! ! Partition RETRO ketones into 75% acetone and 25% MEK (dbm, 8/18/2011)
! ELSE IF ( N == IDTACET ) THEN
! RETRO = RETRO_ACET(I,J)
! ELSE IF ( N == IDTMEK ) THEN
! RETRO = RETRO_MEK(I,J)
!------------------------------------------------------------------------------
ELSE IF ( N == IDTACET ) THEN
RETRO = 0.75d0*RETRO_MEK(I,J) ! RETRO ketones --> 75% ACET
ELSE IF ( N == IDTMEK ) THEN
RETRO = 0.25d0*RETRO_MEK(I,J) ! RETRO ketones --> 25% MEK
These fixes will be added into GEOS-Chem v9-01-02.
--Bob Y. 16:31, 22 August 2011 (EDT)
Further fix for double-counting of residential biofuel
This update was tested in the 1-month benchmark simulation v9-01-02s and approved on 5 Nov 2011.
Chris Homles wrote:
- The fix for double-counting home biofuel emissions with RETRO is *incorrect* as currently described on the wiki here.
- Maybe you already discovered this, but the code needs to update NN (transported tracer ID) within the loop. With this change the BIOFUELS emissions for RETRO species (e.g. ALK4, ACET, MEK...) should be zero. It looks like they are not zero in the v9-01-2o benchmark runs.
--Melissa Payer 12:12, 20 October 2011 (EDT)
RETRO anthropogenic ethane emissions are too low
This update was tested in the 1-month benchmark simulation v9-01-03i and approved on 26 Mar 2012.
Bess Corbitt (Harvard) found that the ethane emissions in RETRO inventory are too low, as evidenced by this plot:
Also, Cyndi Whaley (of U. Toronto) has created a PDF plot comparing GEOS-Chem (v8-02-01 and v9-01-01) to ground-based FTIR measurements for 2006. (The measurements come from Kim Strong's group at U. Toronto.) As you can see, ethane v9-01-02 is much lower than both the measurements and than the previous model version.
The emissions used in Yaping Xiao's offline ethane simulation are much better. Dylan Millet suggested replacing the ethane emissions in RETRO with off-line simulation values from Xiao et al. (2008). This has been implemented in GEOS-Chem v9-01-03.
--Bob Y. 15:41, 23 November 2011 (EST)
Update 3/22/12: Emily Fischer found that that ethane emissions in the EPA/NEI2005, Streets, and EMEP regional inventories are also too low. Therefore, we will replace these with Yaping Xiao's offline ethane simulation emissions as well.
--Melissa Payer 15:18, 23 March 2012 (EDT)
