Difference between revisions of "Halogen chemistry mechanism"

From Geos-chem
Jump to: navigation, search
(Accounting for seasalt bromide in SSA heterogeneous chemistry)
m (Overview)
 
(26 intermediate revisions by 4 users not shown)
Line 1: Line 1:
 
== Overview ==
 
== Overview ==
  
The halogen simulation present in [[GEOS-Chem v11-02#v11-02d|GEOS-Chem v11-2d]] was build atop of [[GEOS-Chem v10-01]] and brought into the main code branch at v11-01g. This code brings together published halogen developments in GEOS-Chem as described in detail by Sherwen et al [2016b], including the methyl iodide simulation [Bell et al 2002], the iodine simulation from Sherwen et al [2016a], the original bromine simulation in GEOS-Chem [Parrella et al 2012], the UCX simulation including stratospheric halogens [Eastham et al 2014], updates to bromine mechanism [Schmidt et al 2016]. Further updates were made use NASA JPL15-10 halogen cross-sections and rates by Sherwen et al [2017].
+
The halogen simulation present in [[GEOS-Chem v11-02#v11-02d|GEOS-Chem v11-02d]] (approved 12 Feb 2018) was built atop of [[GEOS-Chem v10-01]] and brought into the main code branch at v11-01g. This code brings together published halogen developments in GEOS-Chem as described in detail by Sherwen et al [2016b], including the methyl iodide simulation [Bell et al 2002], the iodine simulation from Sherwen et al [2016a], the original bromine simulation in GEOS-Chem [Parrella et al 2012] ([http://wiki.seas.harvard.edu/geos-chem/index.php/Bromine_chemistry_mechanism as described on the wiki]), the UCX simulation including stratospheric halogens [Eastham et al 2014], and updates to bromine mechanism [Schmidt et al 2016]. Further updates were made to use NASA JPL15-10 halogen cross-sections and rates [Sherwen et al. 2017].
  
 
The gas-phase chemistry follows Sherwen et al [2016b], with updates made to the heterogenous code to work with structural changes from v10-01 to v11-01.
 
The gas-phase chemistry follows Sherwen et al [2016b], with updates made to the heterogenous code to work with structural changes from v10-01 to v11-01.
  
Key additions vs. v10-01:
+
 
 +
:Key additions vs. "out of the box" v10-01:
  
 
*Addition of 21 new tracers (for further details on these tracers, please see the [http://wiki.seas.harvard.edu/geos-chem/index.php/Species_in_GEOS-Chem GEOS-Chem species wiki page]). All halogen tracers are shown below, with '''new ones in bold'''. This includes five aerosol tracers for bromine (BrSALA, BrSALC) and iodine (ISALA, ISALC, AERI).
 
*Addition of 21 new tracers (for further details on these tracers, please see the [http://wiki.seas.harvard.edu/geos-chem/index.php/Species_in_GEOS-Chem GEOS-Chem species wiki page]). All halogen tracers are shown below, with '''new ones in bold'''. This includes five aerosol tracers for bromine (BrSALA, BrSALC) and iodine (ISALA, ISALC, AERI).
Line 39: Line 40:
  
 
*Emissions of iodocarbons (CH3I, CH2I2, CH2ICl, CH2IBr) have been included following Ordonez et al [2012]. Emissions of inorganic iodine (HOI, I2) has been using the parameterisation from Carpenter et al [2013] and Macdonald et al [2014].
 
*Emissions of iodocarbons (CH3I, CH2I2, CH2ICl, CH2IBr) have been included following Ordonez et al [2012]. Emissions of inorganic iodine (HOI, I2) has been using the parameterisation from Carpenter et al [2013] and Macdonald et al [2014].
*Wet deposition of X species (HOCl, ClNO3, HI, HOI, INOy, I2OX, AERI) and dry deposition of Y species (HI, HOI, INOy, I2OX, AERI).  
+
*Wet deposition of 10 species (HOCl, ClNO3, HI, HOI, INOy, I2OX, AERI) and dry deposition of 8 species (HI, HOI, INOy, I2OX, AERI).  
  
 
=== Chlorine ===
 
=== Chlorine ===
  
Stratospheric chlorine chemistry is from UCX [Eastham et al 2014]. Tropospheric chlorine chemistry follows that described in Eastham et al [2014] and Schmidt et al [2016], with additions described in Sherwen et al [2016b]. These additions include tropospherically relevant reactions based on the JPL 10-6 compilation [Sander et al., 2011] and IUPAC [Atkinson et al., 2006]. The heterogenous reaction of N2O5 on aerosols was updated to yield
+
Stratospheric chlorine chemistry is from UCX [Eastham et al 2014]. Tropospheric chlorine chemistry follows that described in Eastham et al [2014] and Schmidt et al [2016], with additions described in Sherwen et al [2016b]. These additions include tropospherically relevant reactions based on the JPL 10-6 compilation [Sander et al., 2011] and IUPAC [Atkinson et al., 2006]. The heterogeneous reaction of N2O5 on aerosols was updated to yield
 
products of ClNO2 and HNO3  [Bertram and Thornton, 2009; Roberts et al., 2009] on sea salt and 2HNO3 on other aerosol types. Reaction probabilities are unchanged [Evans and Jacob,
 
products of ClNO2 and HNO3  [Bertram and Thornton, 2009; Roberts et al., 2009] on sea salt and 2HNO3 on other aerosol types. Reaction probabilities are unchanged [Evans and Jacob,
 
2005].
 
2005].
Line 53: Line 54:
 
=== Iodine ===
 
=== Iodine ===
  
The iodine simulation is described in Sherwen et al [2016a], with updates to coupling with bromine and chlorine described in Sherwen et al [2016b]. Notably, the cycling of iodine on aerosol now results in production of ICl and IBr, instead of I2 as in Sherwen et al [2016a].
+
The iodine simulation is described in Sherwen et al [2016a], with updates to coupling with bromine and chlorine described in Sherwen et al [2016b]. Notably, the cycling of iodine on aerosol now results in the production of ICl and IBr, instead of I2 as in Sherwen et al [2016a]. Another key change is that the sea-surface iodide now uses the value predicted by the relationship in MacDonald et al [2014] rather than Chance et al [2014].  
  
 
=== Key differences from publication ===  
 
=== Key differences from publication ===  
  
Certain changes were made during implementation of the halogens code in v10-01 into v11-02 that are not described in Sherwen et al (2016b). These include:  
+
Certain changes were made during the implementation of the halogens code in v10-01 into v11-02 that are not described in Sherwen et al (2016b). These include:  
 
*In combination with UCX [Eastham et al 2014], halogen chemistry is calculated online for both troposphere and stratosphere. In Sherwen et al [2016a, 2016b, 2107] the troposphere was calculated offline.
 
*In combination with UCX [Eastham et al 2014], halogen chemistry is calculated online for both troposphere and stratosphere. In Sherwen et al [2016a, 2016b, 2107] the troposphere was calculated offline.
*Updates to heterogeneous chemistry (see separate section with detail below).  
+
*Updates to heterogeneous chemistry (see the separate section with detail below).  
 
*Minor deposition routes included in previous publications were removed for simplicity (Br2, ICl, IBr, I2).  
 
*Minor deposition routes included in previous publications were removed for simplicity (Br2, ICl, IBr, I2).  
 
*Rates were updated following the [https://jpldataeval.jpl.nasa.gov/ JPL15-10 complication] [Sherwen et al 2017]. Rates were also updated to more consistently use JPL over alternative options (e.g. vs IUPAC).
 
*Rates were updated following the [https://jpldataeval.jpl.nasa.gov/ JPL15-10 complication] [Sherwen et al 2017]. Rates were also updated to more consistently use JPL over alternative options (e.g. vs IUPAC).
Line 67: Line 68:
 
=== Changes in heterogeneous chemistry ===
 
=== Changes in heterogeneous chemistry ===
  
Due to updates structuctural changes in GEOS-Chem (aka Flexchem making the smvgear mechanism obsolete), the heterogenous chemistry routines have been updated to allow for use of more complex heterogenous chemistry previously calculated in Calcrate.F [Eastham et al 2014, Schmidt et al. 2016, Sherwen et al 2016b].
+
Due to updates structural changes in GEOS-Chem (aka Flexchem making the smvgear mechanism obsolete), the heterogenous chemistry routines have been updated to allow for use of more complex heterogeneous chemistry previously calculated in Calcrate.F [Eastham et al 2014, Schmidt et al. 2016, Sherwen et al 2016b].
  
This work required a revamp of the heterogeneous chemistry. As a result, all heterogeneous reactions except for direct uptake (e.g. HO2) and those involving sea salt Cl- (which is not explicitly tracked) are now calculated using a second-order rate constant. This should help to make the het-chem easier to extend, as new functions have been implemented which will automatically convert a pseudo-first-order rate constant into a second order constant while performing safety checks and imposing (if requested) limits on depletion rates to avoid instability.
+
This work required a revamp of the heterogeneous chemistry. As a result, all heterogeneous reactions except for direct uptake (e.g. HO2) and those involving sea salt Cl- (which is not explicitly tracked) are now calculated using a second-order rate constant. This should help to make the het-chem easier to extend, as new functions have been implemented which will automatically convert a pseudo-first-order rate constant into a second-order constant while performing safety checks and imposing (if requested) limits on depletion rates to avoid instability.
 +
 
 +
These changes are implemented using the helper function "kIIrILtd" in gckpp_HetRates.F90. It’s a helper function which converts a pseudo-1st-order reaction to a 2nd-order one, by dividing the reaction rate by the concentration of the other reactant (kII from kI). If provided, a minimum reactant lifetime can also be provided, which has been important in the past for some chlorine activation chemistry (hence Ltd). For example, the reaction rate for ClONO2 + HCl –aerosol--> Cl2 + HNO3 is calculated as the first-order rate kI, but if we write the reaction as ClONO2 = Cl2 + HNO3 then we will get an artificial source of chlorine (while also failing to account for possible chemical depletion of HCl). The answer is to use kIIrILtd to convert kI to kII, where kII = kI / [HCl]. The overall reaction rate will be unaffected, but we can now write the “correct” reaction rate in the KPP reaction file, and will not run into the problems described above in terms of balancing Cl atoms and capturing HCl loss.
  
 
=== Inclusion of sulfur-halogen chemistry from Chen et al [2017] ===
 
=== Inclusion of sulfur-halogen chemistry from Chen et al [2017] ===
Line 86: Line 89:
  
 
In addition to changes seen on implementation of halogen chemistry in v10-01 detailed in publications [Eastham et al 2014, Schmidt et al. 2016, Sherwen et al 2016b]. Please find some details on changes to the standard simulation seen (full 1yr benchmark plots are found [http://wiki.seas.harvard.edu/geos-chem/index.php/GEOS-Chem_v11-02_benchmark_history#v11-02d inked here]).
 
In addition to changes seen on implementation of halogen chemistry in v10-01 detailed in publications [Eastham et al 2014, Schmidt et al. 2016, Sherwen et al 2016b]. Please find some details on changes to the standard simulation seen (full 1yr benchmark plots are found [http://wiki.seas.harvard.edu/geos-chem/index.php/GEOS-Chem_v11-02_benchmark_history#v11-02d inked here]).
 
==== Differences in oxidants ====
 
 
'''< insert benchmark comparisons of O3 code here - using the gridded O3 dataset from Sofen et al >'''
 
  
 
=== Source code and data files ===
 
=== Source code and data files ===
Line 99: Line 98:
 
=== Known issues ===
 
=== Known issues ===
  
*The reimplementation into GEOS-Chem v11-2 code is currently only compatible with the standard simulation (e.g. cannot be compiled with a setting of “NO_REDUCED=no”).
 
 
*CH3I is not currently calculated on-line. Instead the emissions of Bell et al [2002] are currently read in from an offline NetCDF.
 
*CH3I is not currently calculated on-line. Instead the emissions of Bell et al [2002] are currently read in from an offline NetCDF.
 
*For consistency with published work. Certain reactions that are included in the stratosphere (e.g. ClNO3+HCl, BrNO3+HCl, HOCl+HCl and HOCl+HBr) remain unconsidered in the troposphere.  
 
*For consistency with published work. Certain reactions that are included in the stratosphere (e.g. ClNO3+HCl, BrNO3+HCl, HOCl+HCl and HOCl+HBr) remain unconsidered in the troposphere.  
 
*ClNO2 production on sulfate following condensation of HCl is not considered (HCl+N2O5=(sulfate aerosol)>ClNO2) for consistency with Sherwen et al 2016b.
 
*ClNO2 production on sulfate following condensation of HCl is not considered (HCl+N2O5=(sulfate aerosol)>ClNO2) for consistency with Sherwen et al 2016b.
*The “uptake” reaction of HCl results in destruction of Cl, as there is no aerosol-phase Cl tracer.
+
*The “uptake” reaction of HCl results in the destruction of Cl, as there is no aerosol-phase Cl tracer.
 
*There are only aerosol-phase bromine tracers for sea-salt, so Bry uptake cannot occur on sulfate.
 
*There are only aerosol-phase bromine tracers for sea-salt, so Bry uptake cannot occur on sulfate.
 
*Release of IX by heterogenous iodine cycling (HOI, INOy) leads to net production of Cly and Bry as no aqueous phase Cl/Br is considered in the reaction (e.g. excess Cl-/Br- is assumed).
 
*Release of IX by heterogenous iodine cycling (HOI, INOy) leads to net production of Cly and Bry as no aqueous phase Cl/Br is considered in the reaction (e.g. excess Cl-/Br- is assumed).
 +
 +
== Configuring simulations ==
 +
 +
=== Switching on sea-salt bromide emission ===
 +
 +
Sea-salt debromination was "switched on" in Chen et al. [2016], but not in the updated heterogenous Br simulation [Schmidt et al. 2016] and subsequent papers that included this scheme [Sherwen et al 2016a,b,c,2017]. To turn on "sea-salt debromination" in v11-2, bromide is emitted with sea-salt.  This is done in v11-2d by updating the HEMCO extension settings in HEMCO_Config.rc as below:
 +
 +
:from:
 +
 +
    --> Emit Br2            :      false
 +
    --> Br2 scaling          :      1.0
 +
    --> Model sea salt Br-  :      <span style="color:red">'''false'''</span>
 +
    --> Br- mass ratio      :      2.11e-3
 +
 +
:to:
 +
 +
    --> Emit Br2            :      false
 +
    --> Br2 scaling          :      1.0
 +
    --> Model sea salt Br-  :      <span style="color:green">'''true'''</span>
 +
    --> Br- mass ratio      :      2.11e-3
 +
 +
=== Switching off HOBr + S(IV) chemistry ===
 +
 +
Chen et al [2016] describes "HOBr + S(IV)" chemistry, which is included in the standard simulation from v11-2d. Other recent GEOS-Chem halogen papers do not include this chemistry [Eastham et al 2014, Schmidt et al. 2016, Sherwen et al 2016a,b,c,2017]. To switch off this chemistry you just update in gckpp_HetRates.F90 (in KPP/Standard/, at line ~650) as shown below.
 +
 +
:from:
 +
      !--------------------------------------------------------------------
 +
      !  Get parameters for HOBr + S(IV)
 +
      !--------------------------------------------------------------------
 +
   
 +
      ! Cloud bisulfite (HSO3-) concentration [mol/l] from sulfate_mod.F
 +
      <span style="color:red">HSO3conc_Cld = State_Chm%HSO3_AQ(I,J,L) </span>
 +
   
 +
      ! Cloud sulfite (SO3--) concentration [mol/l] from sulfate_mod.F
 +
      <span style="color:red">SO3conc_Cld  = State_Chm%SO3_AQ(I,J,L)  </span>
 +
 +
:to:
 +
      !--------------------------------------------------------------------
 +
      !  Get parameters for HOBr + S(IV)
 +
      !--------------------------------------------------------------------
 +
   
 +
      ! Cloud bisulfite (HSO3-) concentration [mol/l] from sulfate_mod.F
 +
      <span style="color:green">HSO3conc_Cld = TINY(1e+0_fp)  </span>
 +
   
 +
      ! Cloud sulfite (SO3--) concentration [mol/l] from sulfate_mod.F
 +
      <span style="color:green">SO3conc_Cld = TINY(1e+0_fp) </span>
  
 
== On-going development ==
 
== On-going development ==
Line 113: Line 157:
 
== Previous issues that are now resolved ==
 
== Previous issues that are now resolved ==
  
=== Accounting for seasalt bromide in SSA heterogeneous chemistry ===
+
=== Fix for N2O5 uptake ===
 +
 
 +
<span style="color:green">'''''This update was included in [[GEOS-Chem v11-02#v11-02f|v11-02f]] (approved 17 May 2018).'''''</span>
 +
 
 +
NOTE: The bugs listed below are only present in the [[GEOS-Chem v11-02#v11-02d|v11-02d]]/[[GEOS-Chem v11-02#v11-02e|v11-02e]] development code and not in any of the halogen chemistry papers referred to above.
 +
 
 +
'''''Xuan Wang wrote:'''''
 +
 
 +
:In KPP directory, gckpp_HetRates.F90, around line 700:
 +
 
 +
      HET(ind_N2O5,  1) = kIIR1Ltd( spcVec, Ind_('N2O5'), Ind_('H2O'), &
 +
                                          HETN2O5(1.08E2_fp, 1E-1_fp))
 +
 
 +
:This changes the N2O5 uptake rate calculation from 1st order to 2nd order, but no corresponding changes in equation file are found. To fix it, change the equation file from:
 +
 
 +
      N2O5 = 2.000HNO3 :      HET(ind_N2O5,1);
 +
 
 +
:to:
 +
 
 +
      N2O5 + H2O = 2.000HNO3 : HET(ind_N2O5,1);
 +
 
 +
This bug decreases the HNO3 budget and increase N2O5.
 +
 
 +
--[[User:Bmy|Bob Yantosca]] ([[User talk:Bmy|talk]]) 15:59, 17 May 2018 (UTC)
 +
 
 +
=== Fixes in gckpp_HetRates.F90 following implementation in v11-02d ===
 +
 
 +
<span style="color:green">'''''This update was included in [[GEOS-Chem v11-02#v11-02e|v11-02e]] (approved 24 Mar 2018).'''''</span>
 +
 
 +
NOTE: The bugs listed below are only present in the [[GEOS-Chem v11-02#v11-02d|v11-02d]] development code and not in any of the halogen chemistry papers referred to above.
 +
 
 +
'''''Lei  Zhu wrote:'''''
 +
 
 +
:There are still a few bugs in gckpp_HetRates.F90 to be fixed in [[GEOS-Chem v11-02#v11-02e|v11-02d]]. Here is my list.
 +
 
 +
:1. Around line 6163 and 6171, signs of Henry’s volatility constant, CR, should be positive. This bug and its fix have been submitted to GCST by Seb.
 +
 
 +
:2. Around line 4382, 4489, and 4646, the sign of H_HOBr_E is also wrong. It should be defined as:
 +
 
 +
  REAL(fp),  PARAMETER  :: H_HOBr_E  = -6014.0
 +
 
 +
:As suggested by Tomas and Seb, we may want to retrieve Henry’s law parameters in gckpp_HetRates.F90 directly from the species database. This would prevent this kind (both 1 and 2) of bug in future.
  
'''Lei Zhu wrote:'''
+
:3. Cloud fraction (CLDF) should not be used in calculating the cloud volume (Vliq and VIce) in CLD_PARAMS function around line 6074-6075. Cloud water mixing ratio (QL and QI), by [http://acmg.seas.harvard.edu/geos/wiki_docs/geos5/GEOS5_Glossary_11_01_2011.pdf definition], is kg of cloud water per kg of dry air mass in the grid. So QL and QI contain cloud fraction information already. To fix it, change
  
:This bug was found by Xuan Wang. I tested it, and found it has a great impact on Bry budget.
+
      VLiq = CLDF(I,J,L) * QL * AD(I,J,L) / dens_h2o
 +
      VIce = CLDF(I,J,L) * QI * AD(I,J,L) / dens_h2o
  
:About the bug: [[GEOS-Chem v11-02#v11-02d|GEOS-Chem v11-02d]] doesn’t allow sea salt bromide (BrSALA and BrSALC) to be emitted as part of sea salt aerosol. Uptake of gas-phase HBr becomes the only source of bromide in sea salt aerosol. This bug has a great impact on total Bry burden and speciation through heterogeneous reactions.
+
:to
  
:To fix the bug: Turn on sea salt Br- in HEMCO_Config.rc (line 145) as default.
+
      VLiq = QL * AD(I,J,L) / dens_h2o
 +
      VIce = QI * AD(I,J,L) / dens_h2o
  
  --> Model sea salt Br-  :  true
+
:4. Lastly, we may want to reevaluate the impact of bromide from sea salt aerosol given recent updates. Now the model has additional sinks of Bry via oceanic Br+ALD2 and HOBr+S(IV). I will do some sensitivity tests on this topic.
  
--[[User:Melissa Payer|Melissa Sulprizio]] ([[User talk:Melissa Payer|talk]]) 20:18, 8 February 2018 (UTC)
+
--[[User:Melissa Payer|Melissa Sulprizio]] ([[User talk:Melissa Payer|talk]]) 19:38, 21 February 2018 (UTC)
  
 
== References ==
 
== References ==

Latest revision as of 16:51, 15 July 2019

Overview

The halogen simulation present in GEOS-Chem v11-02d (approved 12 Feb 2018) was built atop of GEOS-Chem v10-01 and brought into the main code branch at v11-01g. This code brings together published halogen developments in GEOS-Chem as described in detail by Sherwen et al [2016b], including the methyl iodide simulation [Bell et al 2002], the iodine simulation from Sherwen et al [2016a], the original bromine simulation in GEOS-Chem [Parrella et al 2012] (as described on the wiki), the UCX simulation including stratospheric halogens [Eastham et al 2014], and updates to bromine mechanism [Schmidt et al 2016]. Further updates were made to use NASA JPL15-10 halogen cross-sections and rates [Sherwen et al. 2017].

The gas-phase chemistry follows Sherwen et al [2016b], with updates made to the heterogenous code to work with structural changes from v10-01 to v11-01.


Key additions vs. "out of the box" v10-01:
  • Addition of 21 new tracers (for further details on these tracers, please see the GEOS-Chem species wiki page). All halogen tracers are shown below, with new ones in bold. This includes five aerosol tracers for bromine (BrSALA, BrSALC) and iodine (ISALA, ISALC, AERI).
CHCl3 IONO BrNO2 HCFC142b Cl2O2
CH2Cl2 IONO2 BrNO3 CFC11 INO
CH3I I2O2 CHBr3 CFC12 HBr
CH2I2 I2O3 CH2Br2 HCFC22 HCFC141b
CH2ICl I2O4 CH3Br H1211 Cl2
CH2IBr ISALA BrCl H1301 OIO
HOI ISALC HCl H2402 HOBr
I2 AERI CCl4 Cl HCFC123
IBr BrSALA CH3Cl ClO OClO
ICl BrSALC CH3CCl3 HOCl
I Br2 CFC113 ClNO3
IO Br CFC114 ClNO2
HI BrO CFC115 ClOO
  • Emissions of iodocarbons (CH3I, CH2I2, CH2ICl, CH2IBr) have been included following Ordonez et al [2012]. Emissions of inorganic iodine (HOI, I2) has been using the parameterisation from Carpenter et al [2013] and Macdonald et al [2014].
  • Wet deposition of 10 species (HOCl, ClNO3, HI, HOI, INOy, I2OX, AERI) and dry deposition of 8 species (HI, HOI, INOy, I2OX, AERI).

Chlorine

Stratospheric chlorine chemistry is from UCX [Eastham et al 2014]. Tropospheric chlorine chemistry follows that described in Eastham et al [2014] and Schmidt et al [2016], with additions described in Sherwen et al [2016b]. These additions include tropospherically relevant reactions based on the JPL 10-6 compilation [Sander et al., 2011] and IUPAC [Atkinson et al., 2006]. The heterogeneous reaction of N2O5 on aerosols was updated to yield products of ClNO2 and HNO3 [Bertram and Thornton, 2009; Roberts et al., 2009] on sea salt and 2HNO3 on other aerosol types. Reaction probabilities are unchanged [Evans and Jacob, 2005].

Bromine

The bromine simulation incorporates the work of Parrella et al [2012] (See wiki article here on incorporation in earlier versions), with updates (especially to heterogenous process) as described in detail by Schmidt et al [2016]. Minor updates to Schmidt et al [2016] are described in Sherwen et al [2016b].

Iodine

The iodine simulation is described in Sherwen et al [2016a], with updates to coupling with bromine and chlorine described in Sherwen et al [2016b]. Notably, the cycling of iodine on aerosol now results in the production of ICl and IBr, instead of I2 as in Sherwen et al [2016a]. Another key change is that the sea-surface iodide now uses the value predicted by the relationship in MacDonald et al [2014] rather than Chance et al [2014].

Key differences from publication

Certain changes were made during the implementation of the halogens code in v10-01 into v11-02 that are not described in Sherwen et al (2016b). These include:

  • In combination with UCX [Eastham et al 2014], halogen chemistry is calculated online for both troposphere and stratosphere. In Sherwen et al [2016a, 2016b, 2107] the troposphere was calculated offline.
  • Updates to heterogeneous chemistry (see the separate section with detail below).
  • Minor deposition routes included in previous publications were removed for simplicity (Br2, ICl, IBr, I2).
  • Rates were updated following the JPL15-10 complication [Sherwen et al 2017]. Rates were also updated to more consistently use JPL over alternative options (e.g. vs IUPAC).
  • Updates to VOCs (inc. oceanic acetaldehyde source) in the simulation (see “PAN updates”) will decrease Bry mixing ratios and therefore halogen impacts.
  • Uptake of HOBr by Br-/Cl-/S(IV) in aerosols and cloud droplets now treated as a "one-bulk" process following Chen et al [2017], rather than "multiple" uptake of HOBr by Br-/Cl-/S(IV) in aerosol and cloud droplets.

Changes in heterogeneous chemistry

Due to updates structural changes in GEOS-Chem (aka Flexchem making the smvgear mechanism obsolete), the heterogenous chemistry routines have been updated to allow for use of more complex heterogeneous chemistry previously calculated in Calcrate.F [Eastham et al 2014, Schmidt et al. 2016, Sherwen et al 2016b].

This work required a revamp of the heterogeneous chemistry. As a result, all heterogeneous reactions except for direct uptake (e.g. HO2) and those involving sea salt Cl- (which is not explicitly tracked) are now calculated using a second-order rate constant. This should help to make the het-chem easier to extend, as new functions have been implemented which will automatically convert a pseudo-first-order rate constant into a second-order constant while performing safety checks and imposing (if requested) limits on depletion rates to avoid instability.

These changes are implemented using the helper function "kIIrILtd" in gckpp_HetRates.F90. It’s a helper function which converts a pseudo-1st-order reaction to a 2nd-order one, by dividing the reaction rate by the concentration of the other reactant (kII from kI). If provided, a minimum reactant lifetime can also be provided, which has been important in the past for some chlorine activation chemistry (hence Ltd). For example, the reaction rate for ClONO2 + HCl –aerosol--> Cl2 + HNO3 is calculated as the first-order rate kI, but if we write the reaction as ClONO2 = Cl2 + HNO3 then we will get an artificial source of chlorine (while also failing to account for possible chemical depletion of HCl). The answer is to use kIIrILtd to convert kI to kII, where kII = kI / [HCl]. The overall reaction rate will be unaffected, but we can now write the “correct” reaction rate in the KPP reaction file, and will not run into the problems described above in terms of balancing Cl atoms and capturing HCl loss.

Inclusion of sulfur-halogen chemistry from Chen et al [2017]

Sulfur oxidation by reactive halogens, as described by Chen et al [2017], was included in the same development branch (v11-2d) of the GEOS-Chem code as the updates to the halogen simulation detailed here.

From Chen et al [2017]:

Sulfur and reactive bromine (Bry) play important roles in tropospheric chemistry and the global radiation budget. The oxidation of dissolved SO2 (S(IV)) by HOBr increases sulfate aerosol abundance and may also impact the Bry budget, but is generally not included in global climate and chemistry models. In this study, we implement HOBr + S(IV) reactions into the GEOS-Chem global chemical transport model and evaluate the global impacts on both sulfur and Bry budgets. Modeled HOBr mixing ratios on the order of 0.1–1.0 parts per trillion (ppt) lead to HOBr + S(IV) contributing to 8% of global sulfate production and up to 45% over some tropical ocean regions with high HOBr mixing ratios (0.6–0.9 ppt). Inclusion of HOBr + S(IV) in the model leads to a global Bry decrease of 50%, initiated by the decrease in bromide recycling in cloud droplets. Observations of HOBr are necessary to better understand the role of HOBr + S(IV) in tropospheric sulfur and Bry cycles.

Text S2 in the supporting document for Chen et al [2017] describes the parameterisation of HOBr + S(IV) reactions now included in GEOS-Chem.

Key differences in Simulation

In addition to changes seen on implementation of halogen chemistry in v10-01 detailed in publications [Eastham et al 2014, Schmidt et al. 2016, Sherwen et al 2016b]. Please find some details on changes to the standard simulation seen (full 1yr benchmark plots are found inked here).

Source code and data files

A monthly inventory of Iodocarbon emissions are included from Ordonez et al. [2012]. These files were the result of a project lead by Alfonso Saiz-Lopez and he has kindly allowed for them to be distributed with the main GEOS-Chem code.

These new data files are contained in the HEMCO data directory tree. For detailed instructions on how to download these data files to your disk server, please see our Downloading the HEMCO data directories wiki post.

Known issues

  • CH3I is not currently calculated on-line. Instead the emissions of Bell et al [2002] are currently read in from an offline NetCDF.
  • For consistency with published work. Certain reactions that are included in the stratosphere (e.g. ClNO3+HCl, BrNO3+HCl, HOCl+HCl and HOCl+HBr) remain unconsidered in the troposphere.
  • ClNO2 production on sulfate following condensation of HCl is not considered (HCl+N2O5=(sulfate aerosol)>ClNO2) for consistency with Sherwen et al 2016b.
  • The “uptake” reaction of HCl results in the destruction of Cl, as there is no aerosol-phase Cl tracer.
  • There are only aerosol-phase bromine tracers for sea-salt, so Bry uptake cannot occur on sulfate.
  • Release of IX by heterogenous iodine cycling (HOI, INOy) leads to net production of Cly and Bry as no aqueous phase Cl/Br is considered in the reaction (e.g. excess Cl-/Br- is assumed).

Configuring simulations

Switching on sea-salt bromide emission

Sea-salt debromination was "switched on" in Chen et al. [2016], but not in the updated heterogenous Br simulation [Schmidt et al. 2016] and subsequent papers that included this scheme [Sherwen et al 2016a,b,c,2017]. To turn on "sea-salt debromination" in v11-2, bromide is emitted with sea-salt. This is done in v11-2d by updating the HEMCO extension settings in HEMCO_Config.rc as below:

from:
   --> Emit Br2             :       false
   --> Br2 scaling          :       1.0
   --> Model sea salt Br-   :       false
   --> Br- mass ratio       :       2.11e-3
to:
   --> Emit Br2             :       false
   --> Br2 scaling          :       1.0
   --> Model sea salt Br-   :       true
   --> Br- mass ratio       :       2.11e-3

Switching off HOBr + S(IV) chemistry

Chen et al [2016] describes "HOBr + S(IV)" chemistry, which is included in the standard simulation from v11-2d. Other recent GEOS-Chem halogen papers do not include this chemistry [Eastham et al 2014, Schmidt et al. 2016, Sherwen et al 2016a,b,c,2017]. To switch off this chemistry you just update in gckpp_HetRates.F90 (in KPP/Standard/, at line ~650) as shown below.

from:
     !--------------------------------------------------------------------
     !  Get parameters for HOBr + S(IV)
     !--------------------------------------------------------------------
   
     ! Cloud bisulfite (HSO3-) concentration [mol/l] from sulfate_mod.F
     HSO3conc_Cld = State_Chm%HSO3_AQ(I,J,L) 
   
     ! Cloud sulfite (SO3--) concentration [mol/l] from sulfate_mod.F
     SO3conc_Cld  = State_Chm%SO3_AQ(I,J,L)  
to:
     !--------------------------------------------------------------------
     !  Get parameters for HOBr + S(IV)
     !--------------------------------------------------------------------
   
     ! Cloud bisulfite (HSO3-) concentration [mol/l] from sulfate_mod.F
     HSO3conc_Cld = TINY(1e+0_fp)  
   
     ! Cloud sulfite (SO3--) concentration [mol/l] from sulfate_mod.F
     SO3conc_Cld = TINY(1e+0_fp) 

On-going development

The halogen simulation will developed further in the future. Please see the GEOS-Chem model development priorities for more details.

Previous issues that are now resolved

Fix for N2O5 uptake

This update was included in v11-02f (approved 17 May 2018).

NOTE: The bugs listed below are only present in the v11-02d/v11-02e development code and not in any of the halogen chemistry papers referred to above.

Xuan Wang wrote:

In KPP directory, gckpp_HetRates.F90, around line 700:
     HET(ind_N2O5,  1) = kIIR1Ltd( spcVec, Ind_('N2O5'), Ind_('H2O'), &
                                         HETN2O5(1.08E2_fp, 1E-1_fp))
This changes the N2O5 uptake rate calculation from 1st order to 2nd order, but no corresponding changes in equation file are found. To fix it, change the equation file from:
     N2O5 = 2.000HNO3 :       HET(ind_N2O5,1);
to:
     N2O5 + H2O = 2.000HNO3 : HET(ind_N2O5,1);

This bug decreases the HNO3 budget and increase N2O5.

--Bob Yantosca (talk) 15:59, 17 May 2018 (UTC)

Fixes in gckpp_HetRates.F90 following implementation in v11-02d

This update was included in v11-02e (approved 24 Mar 2018).

NOTE: The bugs listed below are only present in the v11-02d development code and not in any of the halogen chemistry papers referred to above.

Lei Zhu wrote:

There are still a few bugs in gckpp_HetRates.F90 to be fixed in v11-02d. Here is my list.
1. Around line 6163 and 6171, signs of Henry’s volatility constant, CR, should be positive. This bug and its fix have been submitted to GCST by Seb.
2. Around line 4382, 4489, and 4646, the sign of H_HOBr_E is also wrong. It should be defined as:
  REAL(fp),  PARAMETER   :: H_HOBr_E  = -6014.0
As suggested by Tomas and Seb, we may want to retrieve Henry’s law parameters in gckpp_HetRates.F90 directly from the species database. This would prevent this kind (both 1 and 2) of bug in future.
3. Cloud fraction (CLDF) should not be used in calculating the cloud volume (Vliq and VIce) in CLD_PARAMS function around line 6074-6075. Cloud water mixing ratio (QL and QI), by definition, is kg of cloud water per kg of dry air mass in the grid. So QL and QI contain cloud fraction information already. To fix it, change
     VLiq = CLDF(I,J,L) * QL * AD(I,J,L) / dens_h2o
     VIce = CLDF(I,J,L) * QI * AD(I,J,L) / dens_h2o
to
     VLiq = QL * AD(I,J,L) / dens_h2o
     VIce = QI * AD(I,J,L) / dens_h2o
4. Lastly, we may want to reevaluate the impact of bromide from sea salt aerosol given recent updates. Now the model has additional sinks of Bry via oceanic Br+ALD2 and HOBr+S(IV). I will do some sensitivity tests on this topic.

--Melissa Sulprizio (talk) 19:38, 21 February 2018 (UTC)

References

  1. Bell, N., L. Hsu, D. J. Jacob, M. G. Schultz, D. R. Blake, J. H. Butler, D. B. King, J. M. Lobert, and E. Maier-Reimer (2002), Methyl iodide: Atmospheric budget and use as a tracer of marine convection in global models, J. Geophys. Res-Atmos., 107(D17), ACH 8-1–ACH 8-12, doi:http://dx.doi.org/10.1029/2001jd001151
  2. Eastham, S. D., D. K. Weisenstein, and S. R. H. Barrett (2014), Development and evaluation of the unified tropospheric–stratospheric chemistry extension (UCX) for the global chemistry-transport model GEOS-Chem, Atmos. Environ., 89, 52–63, doi:http://dx.doi.org/10.1016/j.atmosenv.2014.02.001
  3. Parrella, J. P. et al. (2012), Tropospheric bromine chemistry: implications for present and pre-industrial ozone and mercury, Atmos. Chem. Phys., 12(15), 6723–6740, doi:http://dx.doi.org/10.5194/acp-12-6723-2012
  4. Schmidt, J. A. et al. (2016), Modeling the observed tropospheric BrO background: Importance of multiphase chemistry and implications for ozone, OH, and mercury, J Geophys. Res-Atmos., doi:http://dx.doi.org/10.1002/2015JD024229
  5. Sherwen, T. et al. (2016a), Iodine’s impact on tropospheric oxidants: a global model study in GEOS-Chem, Atmos. Chem. Phys., 16(2), 1161–1186, doi:http://dx.doi.org/10.5194/acp-16-1161-2016
  6. Sherwen, T. et al. (2016b), Global impacts of tropospheric halogens (Cl, Br, I) on oxidants and composition in GEOS-Chem, Atmos. Chem. Phys., 16(18), 12239–12271, doi:http://dx.doi.org/10.5194/acp-16-12239-2016
  7. MacDonald, S. M., J. C. Gómez Martín, R. Chance, S. Warriner, A. Saiz-Lopez, L. J. Carpenter, and J. M. C. Plane (2014), A laboratory characterisation of inorganic iodine emissions from the sea surface: dependence on oceanic variables and parameterisation for global modelling, Atmos. Chem. Phys., 14(11), 5841–5852, doi:http://dx.doi.org/10.5194/acp-14-5841-2014
  8. Carpenter, L. J., S. M. MacDonald, M. D. Shaw, R. Kumar, R. W. Saunders, R. Parthipan, J. Wilson, and J. M. C. Plane (2013), Atmospheric iodine levels influenced by sea surface emissions of inorganic iodine, Nat. Geosci., 6(2), 108–111, doi:http://dx.doi.org/10.1038/ngeo1687
  9. Ordóñez, C., Lamarque, J.-F., Tilmes, S., Kinnison, D. E., Atlas, E. L., Blake, D. R., Sousa Santos, G., Brasseur, G., and Saiz-Lopez, A.: Bromine and iodine chemistry in a global chemistry-climate model: description and evaluation of very short-lived oceanic sources, Atmos. Chem. Phys., 12, 1423–1447, doi:http://dx.doi.org/10.5194/acp-12-1423-2012 , 2012
  10. Tomás Sherwen, Mathew John J Evans, Roberto Sommariva, Lloyd D. J. Hollis, Stephen Ball, Paul Monks, Christopher Reed, Lucy Carpenter, James D Lee, Grant Forster, Brian Bandy, Claire Reeves and William Bloss, Effects of halogens on European air-quality, Faraday Discussions, 200, 2017, DOI:http://dx.doi.org/10.1039/c7fd00026j
  11. Chen, Q., J. A. Schmidt, V. Shah, L. Jaeglé, T. Sherwen, and B. Alexander, Sulfate production by reactive bromine: Implications for the global sulfur and reactive bromine budgets, Geophys. Res. Lett., 44, 7069–7078, doi:http://dx.doi.org/10.1002/2017GL073812, 2017
  12. Bertram, T. H. and Thornton, J. A.: Toward a general parameterization of N2O5 reactivity on aqueous particles: the competing effects of particle liquid water, nitrate and chloride, Atmos. Chem. Phys., 9, 8351-8363, https://doi.org/10.5194/acp-9-8351-2009, 2009
  13. Sander, S. P., Friedl, R. R., Abbatt, J. P. D., Barker, J. R., Burkholder, J. B., Golden, D. M., Kolb, C. E., Kurylo, M. J., Moortgat, G. K., Wine, P. H., Huie, R. E., and Orkin, V. L.: Chemical kinetics and photochemical data for use in atmospheric studies, Evaluation Number 17, Tech. rep., NASA Jet Propulsion Laboratory, 2011
  14. Atkinson, R., Baulch, D. L., Cox, R. A., Crowley, J. N., Hampson, R. F., Hynes, R. G., Jenkin, M. E., Rossi, M. J., Troe, J., and IUPAC Subcommittee: Evaluated kinetic and photochemical data for atmospheric chemistry: Volume II – gas phase reactions of organic species, Atmos. Chem. Phys., 6, 3625–4055, doi:http://dx.doi.org/10.5194/acp-6-3625-2006, 2006
  15. Atkinson, R., Baulch, D. L., Cox, R. A., Crowley, J. N., Hampson, R. F., Hynes, R. G., Jenkin, M. E., Rossi, M. J., and Troe, J.: Evaluated kinetic and photochemical data for atmospheric chemistry: Volume III - gas phase reactions of inorganic halogens, Atmos. Chem. Phys., 7, 981–1191, 2007
  16. Atkinson, R., Baulch, D. L., Cox, R. A., Crowley, J. N., Hampson, R. F., Hynes, R. G., Jenkin, M. E., Rossi, M. J., Troe, J., and Wallington, T. J.: Evaluated kinetic and photochemical data for atmospheric chemistry: Volume IV - gas phase reactions of organic halogen species, J. Phys. Chem. Ref. Data, 8, 4141–4496, 2008
  17. Evans, M. J. and Jacob, D. J.: Impact of new laboratory studies of N2O5 hydrolysis on global model budgets of tropospheric nitrogen oxides, ozone, and OH, Geophys. Res. Lett., 32, L09 813, doi:http://dx.doi.org/10.1029/2005GL022469, 2005
  18. Roberts, J. M., Osthoff, H. D., Brown, S. S., Ravishankara, A. R., Coffman, D., Quinn, P., and Bates, T.: Laboratory studies of products of N2O5 uptake on Cl containing substrates, Geophys. Res. Lett., 36, doi:http://dx.doi.org/10.1029/2009GL040448, 2009