Chemistry Working Group

From Geos-chem
Revision as of 22:07, 30 August 2010 by Jmao (Talk | contribs) (yield of isoprene nitrates)

Jump to: navigation, search

Oxidants and Chemistry Working Group

All users interested in the GEOS-Chem chemistry scheme and associated processes (photolysis, heterogeneous, deposition) are encouraged to subscribe to the chemistry email list (click on the link in the contact information section below).

Contact information

Current GEOS-Chem Chemistry Projects (please add yours!)

User Group Description Contact Person Date Added
Caltech/Harvard Examining the effects of isoprene nitrates and epoxide on tropical ozone and OH using SHADOZ and OMI measurements Fabien Paulot
Jingqiu Mao
April 28, 2009
Caltech/Harvard Examining the HCHO yield, OH recycling and ozone production with INTEX-A observations from new isoprene chemistry Jingqiu Mao
Fabien Paulot
April 28, 2009
NIA / LaRC Tropospheric ozone over East Asia: Ozonesonde observations and modeling analysis Yiqiang Zhang
Hongyu Liu
June 29, 2010

Current GEOS-Chem Chemistry Issues (please add yours!)

Isoprene chemistry

I've created a page with some of the recent literature on isoprene chemistry. Please add more papers as they come along! ( MJE Leeds)

isoprene

yield of isoprene nitrates

Fabien Paulot found a problem in current chemistry scheme.In v8-02-01 and beyond, isoprene nitrates are produced twice: one through channel A and one through 10% loss in channel B. This makes the loss of NOx larger than it should be (18.7% vs. 10%) and also reduces the yield of MVK/MACR/CH2O by about 13%.

A  453 2.70E-12  0.0E+00    350 1 B   0.00     0.     0.         
       5.00E+00  0.0E+00      0 0     0.00     0.     0.         
      RIO2          +     NO                                              
=0.900NO2           +0.900HO2           +0.340IALD          +0.340MVK     
+0.220MACR          +0.560CH2O          +                   +  
          
A  453 2.70E-12  0.0E+00    350 1 A   0.00     0.     0.         
       5.00E+00  0.0E+00      0 0     0.00     0.     0.         
      RIO2          +     NO                                              
=1.000HNO3          +                   +                   +             

So it should be corrected as (no channel A):

A  453 2.70E-12  0.0E+00    350 1 0   0.00     0.     0.         
       5.00E+00  0.0E+00      0 0     0.00     0.     0.         
      RIO2          +     NO                                              
=0.900NO2           +0.900HO2           +0.340IALD          +0.340MVK     
+0.220MACR          +0.560CH2O          +                   +      

D  453 2.70E-12  0.0E+00    350 1 A   0.00     0.     0.         
       5.00E+00  0.0E+00      0 0     0.00     0.     0.         
      RIO2          +     NO                                              
=1.000HNO3          +                   +                   +  

--J Mao. 18:04, 30 Aug 2010 (EDT)

HO2 + CH2O

Scheme does not contain the HO2 + CH2O --> Adduct reaction (MJE Leeds)

Hermans, I., et al. (2005), Kinetics of alpha-hydroxy-alkylperoxyl radicals in oxidation processes. HO2 center dot-initiated oxidation of ketones/aldehydes near the tropopause, Journal of Physical Chemistry A, 109(19), 4303-4311.

According to this paper, this reaction is significant when Temperature is below 220K.

--J Mao. 15:00, 10 Aug 2009 (EDT)

near-IR photolysis of HNO4

1. Since FastJX already takes this into account with cross section data at 574nm, we do not need to redo this in calcrate.f. We can therefore comment out this entire IF block:

        !---------------------------------------------------------------------
        ! Prior to 10/27/09:
        ! FastJX has taken near-IR photolysis into account with
        ! cross section at 574nm, so we don't need to add 1e-5 anymore.
        ! According to Jimenez et al., "Quantum yields of OH, HO2 and
        ! NO3 in the UV photolysis of HO2NO2", PCCP, 2005, we also
        ! changed the branch ratio from 0.67(HO2)/0.33(OH) to 0.95/0.05
        ! This will put most weight of near-IR photolysis on HO2 channel.
        ! (jmao, bmy, 10/27/09)
        !
        !!==============================================================
        !! HARDWIRE addition of 1e-5 s-1 photolysis rate to 
        !! HNO4 -> HO2+NO2 to account for HNO4 photolysis in near-IR -- 
        !! see Roehl et al. 'Photodissociation of peroxynitric acid in 
        !! the near-IR', 2002. (amf, bmy, 1/7/02)
        !!
        !! Add NCS index to NKHNO4 for SMVGEAR II (gcc, bmy, 4/1/03)
        !!==============================================================
        !IF ( NKHNO4(NCS) > 0 ) THEN
        !
        !   ! Put J(HNO4) in correct spot for SMVGEAR II
        !   PHOTVAL = NKHNO4(NCS) - NRATES(NCS)
        !   NKN     = NKNPHOTRT(PHOTVAL,NCS)
        !
        !   DO KLOOP=1,KTLOOP
        !      RRATE(KLOOP,NKN)=RRATE(KLOOP,NKN) + 1d-5
        !   ENDDO
        !ENDIF
        !---------------------------------------------------------------------


2. We need to change the branch ratio of HNO4 photolysis in ratj.d. Change these lines from:

13 HNO4       PHOTON     OH         NO3                  0.00E+00  0.00     33.3  HO2NO2 
14 HNO4       PHOTON     HO2        NO2                  0.00E+00  0.00     66.7  HO2NO2 

to:

13 HNO4       PHOTON     OH         NO3                  0.00E+00  0.00      5.0  HO2NO2 
14 HNO4       PHOTON     HO2        NO2                  0.00E+00  0.00     95.0  HO2NO2

This is based on Jimenez et al. (Quantum yields of OH, HO2 and NO3 in the UV photolysis of HO2NO2, PHYSICAL CHEMISTRY CHEMICAL PHYSICS, 2005) shows that HO2 yield should be 0.95 and OH yield should be 0.05 for wavelength above 290nm.

This way all the near-IR photolysis will have most weight on HO2 channel(Stark et al., Overtone dissociation of peroxynitric acid (HO2NO2): Absorption cross sections and photolysis products, JOURNAL OF PHYSICAL CHEMISTRY A, 2008).

This update has now been added to the chemistry mechanism documentation file.

--J Mao. 11:00, 26 Oct 2009 (EDT)

Documentation

--Bob Y. 15:41, 27 October 2009 (EDT)