Difference between revisions of "Halogen chemistry mechanism"

From Geos-chem
Jump to: navigation, search
m (References)
(Key differences from publications)
Line 62: Line 62:
 
*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 JPL15-10 complication rates [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 JPL15-10 complication rates [Sherwen et al 2017]. Rates were also updated to more consistently use JPL over alternative options (e.g. vs IUPAC).
*Updates to VOCs in the simulation (see “PAN updates”) will decrease Bry mixing ratios and therefore halogen impacts.
+
*Updates to VOCs in the simulation (see [http://wiki.seas.harvard.edu/geos-chem/index.php/PAN “PAN updates”]) will decrease Bry mixing ratios and therefore halogen impacts.
  
 
=== Changes in heterogeneous chemistry ===
 
=== Changes in heterogeneous chemistry ===

Revision as of 17:59, 10 January 2018

Overview

The halogen simulation present in 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 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:

  • Addition of 21 new tracers. 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 X species (HOCl, ClNO3, HI, HOI, INOy, I2OX, AERI) and dry deposition of Y 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 heterogenous 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 production of ICl and IBr, instead of I2 as in Sherwen et al [2016a].

Key differences from publications

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:

  • In combination with UCX [Eastham et al 2014], halogen chemistry is calculated online for both troposphere and stratosphere. In Sherwen et al [2016, 2016b] the troposphere was calculated offline.
  • Updates to heterogeneous chemistry (see 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 rates [Sherwen et al 2017]. Rates were also updated to more consistently use JPL over alternative options (e.g. vs IUPAC).
  • Updates to VOCs in the simulation (see “PAN updates”) will decrease Bry mixing ratios and therefore halogen impacts.

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].

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.

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 benchmark plots are found here: <insert link>).

Differences in oxidants

< insert benchmark comparisons of O3 code here - using the gridded O3 dataset from Sofen et al >

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

  • 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.
  • 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 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).

Previous issues that are now resolved

None yet

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