The community atmospheric chemistry box model CAABA/MECCA-4.0

Sander, Rolf; Baumgaertner, Andreas; Cabrera-Perez, David; Frank, Franziska; Gromov, Sergey; Grooß, Jens-Uwe; Harder, Hartwig; Huijnen, Vincent; Jöckel, Patrick; Karydis, Vlassis A.; Niemeyer, Kyle E.; Pozzer, Andrea; Riede, Hella; Schultz, Martin G.; Taraborrelli, Domenico; Tauer, Sebastian

doi:https://doi.org/10.5194/gmd-12-1365-2019

Articles | Volume 12, issue 4

https://doi.org/10.5194/gmd-12-1365-2019

Special issue:

https://doi.org/10.5194/gmd-12-1365-2019

Articles | Volume 12, issue 4

Model description paper

05 Apr 2019

Model description paper |

| 05 Apr 2019

The community atmospheric chemistry box model CAABA/MECCA-4.0

Rolf Sander, Andreas Baumgaertner, David Cabrera-Perez, Franziska Frank, Sergey Gromov, Jens-Uwe Grooß, Hartwig Harder, Vincent Huijnen, Patrick Jöckel, Vlassis A. Karydis, Kyle E. Niemeyer, Andrea Pozzer, Hella Riede, Martin G. Schultz, Domenico Taraborrelli, and Sebastian Tauer

Abstract

We present version 4.0 of the atmospheric chemistry box model CAABA/MECCA that now includes a number of new features: (i) skeletal mechanism reduction, (ii) the Mainz Organic Mechanism (MOM) chemical mechanism for volatile organic compounds, (iii) an option to include reactions from the Master Chemical Mechanism (MCM) and other chemical mechanisms, (iv) updated isotope tagging, and (v) improved and new photolysis modules (JVAL, RADJIMT, DISSOC). Further, when MECCA is connected to a global model, the new feature of coexisting multiple chemistry mechanisms (PolyMECCA/CHEMGLUE) can be used. Additional changes have been implemented to make the code more user-friendly and to facilitate the analysis of the model results. Like earlier versions, CAABA/MECCA-4.0 is a community model published under the GNU General Public License.

Download & links

Article (PDF, 3408 KB)

Supplement (25206 KB)

Download & links

Article (3408 KB)
Supplement (25206 KB)
BibTeX
EndNote

How to cite

How to cite.

Dates

Received: 08 Aug 2018 – Discussion started: 03 Sep 2018 – Revised: 11 Mar 2019 – Accepted: 12 Mar 2019 – Published: 05 Apr 2019

1 Introduction

MECCA (Module Efficiently Calculating the Chemistry of the Atmosphere) is an atmospheric chemistry module that contains a comprehensive chemical mechanism with tropospheric and stratospheric chemistry of both the gas and the aqueous phases. For the numerical integration, MECCA uses the KPP (Kinetic PreProcessor) software (Sandu and Sander, 2006).

To apply the MECCA chemistry to atmospheric conditions, MECCA must be connected to a base model via the MESSy (Modular Earth Submodel System) interface (Jöckel et al., 2010). This base model can be a complex 3-dimensional model but it can also be a simple box model. CAABA (Chemistry As A Boxmodel Application) is such a box model, simulating the atmospheric environment in which the MECCA chemistry takes place.

A full description of CAABA/MECCA has already been published elsewhere (Sander et al., 2005; R. Sander et al., 2011). Here, we only present new features that have been implemented after version 3.0. Section 2 describes all changes related to the chemical mechanism of MECCA. In Sect. 3 we show several new options for calculating photolysis rate coefficients in the model. Section 4 presents new features that are only useful when MECCA is coupled to a global (3-dimensional) base model.

2 The chemical mechanism MECCA

MECCA is a chemistry submodel that contains a comprehensive atmospheric reaction mechanism. In addition to the basic HO_x (OH+HO₂), NO_x (NO+NO₂), and CH₄ chemistry, it also includes nonmethane volatile organic compounds (NMVOCs), halogens (Cl, Br, I), sulfur (S), and mercury (Hg) chemistry. Recent extensions of MECCA are presented in the following sections.

2.1 The Mainz Organic Mechanism (MOM)

The MOM is the default oxidation mechanism for volatile organic compounds (VOCs) in MECCA. The current MOM mechanism is a further development of the versions used by Lelieveld et al. (2016) and Cabrera-Perez et al. (2016). It includes developments from Taraborrelli et al. (2012), Hens et al. (2014), and Nölscher et al. (2014). MOM chemistry has been used by Mallik et al. (2018) to study oxidation processes in the Mediterranean atmosphere. Figure 1 shows all 43 emitted species that are treated by MOM. These species are alkanes and alkenes up to four carbon atoms, ethyne (acetylene), two nitriles, isoprene, 2-methyl-3-buten-2-ol (MBO), five monoterpenes, and nine aromatics. Most of the oxidation scheme is explicit. Lumping is used for some isomers with similar properties, e.g., the MOM species “LXYL” presents the sum of o-, m-, and p-xylene. All lumped species are marked by the prefix “L” in their names. The full mechanism includes about 600 species and 1600 reactions. A list of all chemical reactions, including rate coefficients and references, is available in the Supplement (meccanism.pdf).

https://www.geosci-model-dev.net/12/1365/2019/gmd-12-1365-2019-f01

Figure 1Emitted VOCs treated by MOM.

Download

The mechanism for the isoprene oxidation was developed starting from MIM2 (Taraborrelli et al., 2009), which is a reduction of MCM v3.1 (Master Chemical Mechanism; Rickard and Pascoe, 2009; Jenkin et al., 1997). The major mechanisms that regenerate OH under low-NO_x conditions are included. OH addition to the unsaturated isoprene hydroperoxides has been implemented yielding, entirely, epoxydiols and OH according to Paulot et al. (2009). The Z-1,4- and Z-4,1-ISOPO2 isomers undergo 1,6-H-shifts as originally proposed by Peeters et al. (2009). In MOM, the corresponding rate coefficients are those computed by Taraborrelli et al. (2012), and the 66 % yields of isoprene-derived hydroperoxyenals (HPALDs) are according to Nölscher et al. (2014). For the non-HPALD-yielding channel, the corresponding mechanisms proposed by Peeters et al. (2014) and Jenkin et al. (2015) have been included but in a simplified manner. The estimated photo-induced cascade of reactions produces substantial amounts of OH (see Sect. 2.1.5). Finally, methacrolein (MACR) oxidation has been implemented according to Orlando et al. (1999), except for the fate of the methylvinyl radical. The rate of the 1,4-H-shift for the MACRO2 radical is now calculated using the expression reported by Crounse et al. (2012).

Oxidation of the two important terpenes, α-pinene and β-pinene, is based on MCM (Jenkin et al., 2000). However, important modifications following the theoretical work of L. Vereecken have been implemented with some simplifications (Vereecken et al., 2007; Nguyen et al., 2009; Vereecken and Peeters, 2012; Capouet et al., 2008). For instance, minor channels of the OH- and O₃-initiated oxidation are neglected.

Aromatics (benzene, toluene, xylenes) are oxidized in the mechanism by Cabrera-Perez et al. (2016), which is to a large extent a reduction of the corresponding MCM (Jenkin et al., 2003; Bloss et al., 2005). Photolysis of ortho-nitrophenols yielding HONO has been added according to Bejan et al. (2006) and Chen et al. (2011). Finally, reactions of phenyl peroxy radicals with NO₂ yielding NO₃ have been added, consistent with Jagiella and Zabel (2007).

Oxidation of VOCs by O₃ and NO₃ is similar to that in MCM. The oxidation by OH, however, significantly differs from MCM treatment and therefore is detailed in the next section.

2.1.1 VOC reactions with OH

Reactions of OH with organic molecules can be either H abstraction or OH addition. If available, experimental rate coefficients are preferred and taken mostly from the IUPAC kinetic data evaluation (Atkinson et al., 2006; Wallington et al., 2018). Unmeasured rate coefficients for the C₁ to C₅ species are estimated with a site-specific structure–activity relationship (SAR) similar to MCM, based on the work of Atkinson (1987) and Kwok and Atkinson (1995). The base rate coefficients for OH addition to double bonds are taken from the more recent SAR by Peeters et al. (2007). For the C₆ to C₁₁ closed-shell species, the MCM rate coefficients are retained. It is worth noting that the SAR-estimated ones have no temperature dependence and are only given at 298 K. The effect of neighboring groups is expressed by substituent factors and is differentiated by functional group. Most substituent factors by Kwok and Atkinson (1995) are updated or calculated ex novo by computing the relative rate coefficient of OH with the simplest VOC bearing the substituent relative to the one of its parent compound (Table 1). A clear limitation of this approach is that for OH addition, no substituent effect on the branching ratios is considered. No rigorous evaluation of the SAR has been conducted and the estimation uncertainty is expected to be in the same range as for the SAR used by MCM.

Table 1Structure–activity relationship (SAR) parameters and substituent factors in MOM, largely based on Kwok and Atkinson (1995) for H abstraction and on Peeters et al. (2007) for OH addition, unless noted otherwise. Most base rate constants and substituent factors are updated with data from Atkinson et al. (2006). Original values for the substituent factors given by Kwok and Atkinson (1995) are listed in parentheses. All rate constants refer to reactions with OH.

^a There is a sign error in Kwok and Atkinson (1995) who present the value $\exp (- 696 K / T)$ instead of exp (696 K∕T). ^b Median value from the range calculated by Vereecken and Peeters (2001).

Download Print Version | Download XLSX

The general formulae for H abstraction by OH are

\begin{array}{l} (1) & k ({CH}_{3} X) & = k_{p} \cdot F (X), \\ (2) & k ({CH}_{2} X Y) & = k_{s} \cdot F (X) \cdot F (Y), \\ (3) & k (CH X Y Z) & = k_{t} \cdot F (X) \cdot F (Y) \cdot F (Z), \end{array}

where k_p, k_s, and k_t are the group rate coefficients for the hydrogens on the primary, secondary, and tertiary carbon atoms, respectively, and F(X) is the factor for the substituent X.

The SAR for OH addition to (poly)alkenes is based on the hypothesis that the site-specific rate coefficient depends solely on the stability of the radical product (Peeters et al., 2007). Thus, rate coefficients for the formation of primary, secondary, and tertiary radicals are derived from the high-pressure limits for ethene, 2-butene, and 2,3-dimethyl-2-butene, respectively. It is worth noting that for the tertiary radical formation, Peeters et al. (2007) solely used the rate coefficient for 2,3-dimethyl-2-butene and not that for 2-methyl-2-butene minus that for the secondary radical.

2.1.2 RO₂ reactions with NO_x and NO₃

Reactions with NO are the dominant sink for RO₂ under polluted conditions. The RO₂-size independent MCM rate coefficient is used with the exception of CH₃O₂ and CH₃CH₂O₂, for which the IUPAC recommendations are followed (Atkinson et al., 2006). In general, the two possible reaction channels are considered:

\begin{array}{l} (R1) & {RO}_{2} + NO & \to (1 - α) \times (RO + {NO}_{2}) \\ (R2) & \to α \times {RONO}_{2} \end{array}

with α being the yield for the formation of alkyl nitrates, which curb tropospheric ozone production. Acyl RO₂ do not form nitrates. The CH₃ONO₂ yield is calculated according to Butkovskaya et al. (2012) with a reduction according to Flocke et al. (1998). The CH₃CH₂ONO₂ yield is calculated according to Butkovskaya et al. (2010). For all other peroxy radicals, the corresponding alkyl nitrate yields are calculated with the relationship by Arey et al. (2001), which depends on temperature, pressure, and molecular size. However, the latter is not represented by the number of carbon atoms but by the number of heavy atoms (excluding the −OO moiety) according to Teng et al. (2015). The oxygen atom in β-carbonyl RO₂ is not counted. Due to disagreement in the literature, no dependence of α on the degree of RO₂ substitution (primary, secondary, and tertiary) is considered. Reduction factors for β- and γ-carbonyl RO₂ are derived from Praske et al. (2015) and for bicyclic RO₂ from aromatics are derived from Elrod (2011). As an example, Fig. 2 shows the predicted variable yield for the nitrate of the secondary hydroxy butyl peroxy radical.

https://www.geosci-model-dev.net/12/1365/2019/gmd-12-1365-2019-f02

Figure 2Temperature- and pressure-dependent nitrate yield for the secondary hydroxybutyl peroxy radical obtained and calculated by MOM. A constant yield of about 10 % (“Old model”) is used by MCM.

Download

Formation and decomposition of many peroxy nitrates is considered. The equilibria of acyl peroxy nitrates with their parent RO₂ are represented as in MCM but the JPL (Jet Propulsion Laboratory) kinetic data (Burkholder et al., 2015) are used. Only three alkyl peroxy nitrates – CH₃O₂NO₂, CH₃CH₂O₂NO₂, and CH₃COCH₂O₂NO₂ – are represented. The equilibrium reactions for the latter are taken from Tyndall et al. (2001), Sehested et al. (1998), and Kirchner et al. (1999). Reactions of peroxy radicals with NO₃ all produce the corresponding alkoxy radical and NO₂:

\begin{matrix} (R3) & {RO}_{2} + {NO}_{3} \to RO + {NO}_{2} + O_{2} . \end{matrix}

The temperature-independent rate coefficient of $k (C_{2} H_{5} O_{2} + {NO}_{3}) = 2.3 \times 10^{- 12}$ ${cm}^{- 3} s^{- 1}$ is used for all RCH₂O₂. For acyl peroxy radicals, an enhancement factor of $k ({CH}_{3} C (O) OO + {NO}_{3}) / k (C_{2} H_{5} O_{2} + {NO}_{3}) = 1.74$ is calculated based on the peroxy acetyl radical.

2.1.3 RO₂ reactions with HO_x

HO₂ reactions are often competitive with NO reactions of peroxy radicals. The former reactions are known to proceed via three channels:

\begin{array}{l} (R4) & {RO}_{2} + {HO}_{2} & \to RO + OH + O_{2}, \\ (R5) & \to ROOH + O_{2}, \\ (R6) & \to ROH + O_{3}, \end{array}

of which only the first is a radical-propagating channel. Alkyl peroxy radicals cannot have the O₃ channel and their rate coefficient is calculated as a function of the number of carbons according to the fitting formula provided by Saunders et al. (2003) and Boyd et al. (2003). The branching ratios of the OH channel for β-carbonyl, alkoxy, and bicyclic peroxy radicals are taken from Dillon and Crowley (2008), Orlando and Tyndall (2012), and Birdsall et al. (2010), respectively. A 10 % OH yield for reactions of β-hydroxyl peroxy radicals is taken from the isoprene oxidation study of Liu et al. (2013), which is consistent with the results of Groß (2013) and Paulot et al. (2009). The HO₂ reaction of the simplest acyl peroxy radical (CH₃CO₃) has unique branching ratios as determined by direct OH and O₃ measurements (Groß et al., 2014). For all other acyl peroxy radicals, the kinetic data for β-hydroxy acyl peroxy radicals, e.g., HOCH₂CO₃, are taken from Groß (2013) with the rate coefficient having the temperature dependence as recommended by IUPAC.

There is laboratory evidence for a non-negligible reaction of CH₃O₂ with OH (Bossolasco et al., 2014):

\begin{matrix} (R7) & {CH}_{3} O_{2} + OH \to {CH}_{3} O + {HO}_{2} . \end{matrix}

The lower limit of the rate coefficient $1.4 \times 10^{- 10}$ ${cm}^{- 3} s^{- 1}$ reported by Bossolasco et al. (2014) is used in MOM. This is consistent with the revised experimental value by the same laboratory (Assaf et al., 2016). The major reaction channel involving HO₂ elimination represents (80±20) % and is set as the only channel (Assaf et al., 2017). The other possible channels are very uncertain and are therefore not included.

2.1.4 RO₂ permutation reactions

Atkinson et al. (2006)Glover and Miller (2005)Orlando and Tyndall (2012)Wallington et al. (2018)Orlando and Tyndall (2012)Orlando and Tyndall (2012)Orlando and Tyndall (2012)Saunders et al. (2003)Saunders et al. (2003)Saunders et al. (2003)

Table 2Second-order rate constants k^2nd for permutation reactions (in ${cm}^{- 3} s^{- 1}$ ). Here, $k_{{CH}_{3} O_{2}} = 1.03 E - 13 \times \exp (365 K / T)$ ${cm}^{- 3} s^{- 1}$ is for the self reaction of CH₃O₂.

Download Print Version | Download XLSX

The self and cross reactions of organic peroxy radicals are treated according to the permutation reaction formalism in MCM (Jenkin et al., 1997). Every organic peroxy radical reacts in a pseudo-first-order reaction with a rate coefficient that is expressed as

\begin{matrix} (4) & k^{1 st} = 2 \times \sqrt{k_{{RO}_{2}} \times k_{{CH}_{3} O_{2}}} \times [{RO}_{2}], \end{matrix}

where $k_{{RO}_{2}}$ is the second-order rate coefficient of the self reaction of the organic peroxy radical, $k_{{CH}_{3} O_{2}}$ is the second-order rate coefficient of the self reaction of CH₃O₂, and [RO₂] is the sum of the concentrations of all organic peroxy radicals. The formalism is a simplification of the approach by Madronich and Calvert (1990) under the assumption that the dominant co-reactant of RO₂ is CH₃O₂. The value of $k_{{CH}_{3} O_{2}}$ is taken from the IUPAC recommendations. Expressions for $k_{{RO}_{2}}$ distinguish acyl from alkyl peroxy radicals. The latter are differentiated by the degree and kind of substituents close to the −OO moiety. The rate expressions (Table 2) are not from MCM, except for β-hydroxyl radicals, and have a temperature dependence.

2.1.5 Photo-induced reactions

The enhanced photolysis of carbonyl nitrates from isoprene is implemented according to Barnes et al. (1993) and Müller et al. (2014). The enhancement is applied to the photolysis rate coefficients (j values) of nitrooxyacetone (NOA), nitrooxyacetaldehyde (NO3CH2CHO), lumped nitrates of methyl ethyl ketone (LMEKNO3), nitrates of MVK (methyl vinyl ketone) and MACR, and unsaturated C₅-nitrooxyaldehyde from the isoprene + NO₃ reaction.

Keto–enol tautomerization of aldehydes induced by light absorption is implemented based on data for acetaldehyde (Clubb et al., 2012). The enols are in equilibrium with the corresponding aldehydes by HCOOH catalysis (da Silva, 2010). Formic acid is then produced upon reaction of the enols with OH similarly to the simplest enol (So et al., 2014). Vinyl alcohol is also produced in the photolysis of propanal.

Photolysis of HPALDs is according to Peeters et al. (2014) and Jenkin et al. (2015) and the subsequent photolysis of the resulting carbonyl enols HVMK (hydroxy vinyl methyl ketone) and HMAC (3-hydroxy-2-methyl-acrolein) is treated according to Nakanishi et al. (1977) and Messaadia et al. (2015).

Nitrophenols undergo photolysis yielding HONO, according to Bejan et al. (2006) and Chen et al. (2011), and assumed co-products being cyclic ketenes. However, the OH-formation channel (Cheng et al., 2009; Vereecken et al., 2016) is not implemented.

Conjugated unsaturated dialdehydes (like butenedial and 2-methyl-butenedial resulting from the oxidation of isoprene and aromatics) undergo photolysis based on Xiang et al. (2007) and MCM. Only the major channel, CO loss, is considered, and the j values are scaled with j(NO₂). The ketenes from photolysis of hydroperoxyacetyl-conjugated unsaturated aldehydes from isoprene, conjugated unsaturated dialdehydes and nitrophenols undergo photodissociation yielding CO and an excited Criegee intermediate. The j value is assumed to be the same as that for MVK with a unity quantum yield.

2.2 Other chemical mechanisms

In addition to the native chemistry mechanism of MECCA (available in the file gas.eqn), several other independent mechanisms are now provided as well. The chemical mechanisms CB05BASCOE and MOZART from the Copernicus Atmosphere Monitoring Service project (CAMS 42) and the Jülich Atmospheric Mechanism (JAM002) have been converted to KPP format and introduced into MECCA. It is also possible to use our previous simple mechanism MIM1 (Jöckel et al., 2016). In addition, chemical mechanisms extracted and downloaded from the MCM web page can be converted with a script that makes them compatible with CAABA/MECCA. All mechanisms are suitable for stratospheric as well as tropospheric calculations. They all include the chemistry of chlorine, bromine, and isoprene. They differ in the treatment of terpenes. MIM1 has no terpenes at all. CB05BASCOE and MOZART include terpenes as a lumped species. Only MCM, MOM, and JAM002 treat some terpenes individually, e.g., pinene. The JAM002 mechanism is larger than CB05BASCOE and MOZART but small compared to MECCA with MOM. The very detailed MCM is the largest of all. More information about the chemical mechanisms is provided in the following sections.

https://www.geosci-model-dev.net/12/1365/2019/gmd-12-1365-2019-f03

Figure 3Intercomparison of the MOM (black), CB05BASCOE (red), MOZART (green), MIM1 (blue), MCM (magenta), and JAM002 (cyan) mechanisms. The simulations represent the boundary layer over the Amazon forest. They start on 1 August at midnight and last for 5 days. Temperature, pressure, and relative humidity are set to 301 K, 101 325 Pa, and 70 %, respectively. The model is initialized with 2 nmol mol⁻¹ isoprene (C₅H₈), 500 pmol mol⁻¹ of terpenes (MOM: 100 pmol mol⁻¹ of α-pinene, β-pinene, camphene, carene, and sabinene each; CB05BASCOE and MOZART: lumped terpenes; MIM1: no terpenes; MCM and JAM002: 200 pmol mol⁻¹ α-pinene and 300 pmol mol⁻¹ β-pinene), and 100 pmol mol⁻¹ PAN. During the model simulation, emissions of NO are set to $3.3 \times 10^{- 9}$ ${cm}^{- 2} s^{- 1}$ (Taraborrelli et al., 2009).

Download

2.2.1 CB05BASCOE

The CB05BASCOE scheme (Huijnen et al., 2016) is a merge of a tropospheric and stratospheric chemistry scheme. The tropospheric chemistry is based on the Carbon Bond mechanism 2005 (CB05, Yarwood et al., 2005). Here, a lumping approach is adopted for organic species by defining a separate tracer species for specific types of functional groups. The scheme has been modified and extended to include an explicit treatment of C₁ to C₃ species (Williams et al., 2013), SO₂, dimethyl sulfide (DMS), methyl sulfonic acid (MSA), and ammonia (NH₃), as described by Huijnen et al. (2010). The reaction rates follow the recommendations given in either the JPL or IUPAC evaluation (Burkholder et al., 2015; Wallington et al., 2018). The stratospheric chemistry is based on that from the BASCOE (Belgian Assimilation System for Chemical ObsErvations) system (Errera et al., 2008) and is labeled “sb15b”. This chemical scheme merges the reaction lists developed by Errera and Fonteyn (2001) to produce short-term analyses, with the list included in the SOCRATES 2-D model for long-term studies of the middle atmosphere (Brasseur et al., 2000; Chabrillat and Fonteyn, 2003). The list of species includes all the ozone-depleting substances and greenhouse gases necessary for multidecadal simulations of the couplings between dynamics and chemistry in the stratosphere, as well as the reservoir and short-lived species necessary for a complete description of stratospheric ozone photochemistry. Gas-phase and heterogeneous reaction rates are taken from the JPL evaluations 17 and 18 (S. P. Sander et al., 2011; Burkholder et al., 2015). The merged reaction mechanism includes 99 species interacting through 211 gas-phase and 10 heterogeneous reactions. Details regarding its implementation and evaluation within the ECMWF Integrated Forecasting System (IFS) are given by Huijnen et al. (2016). The complete mechanism of CB05BASCOE (species and equations) can be found in the directory mecca/eqn/cb05bascoe/ in the Supplement.

2.2.2 MOZART

The tropospheric chemistry in MOZART (Model of OZone And Related Tracers) is based on the MOZART-3 mechanism by Kinnison et al. (2007). It includes additional species and reactions from MOZART-4 (Emmons et al., 2010) and further updates from the Community Atmosphere Model with interactive chemistry, referred to as CAM4-chem (Lamarque et al., 2012). The chemical mechanism includes an updated isoprene oxidation scheme and a better treatment of volatile organic compounds, with lumped species to represent large alkanes, alkenes, and aromatic compounds as well as their oxidation products. Overall, it includes the degradation of C₁, C₂, C₃, C₄, C₅, C₇, and C₁₀ species. The heterogeneous chemistry in the troposphere is implemented according to the corresponding module from CB05BASCOE. MOZART includes the extended stratospheric chemistry discussed by Kinnison et al. (2007) with further updates from CAM4-chem (Lamarque et al., 2012; Tilmes et al., 2016). This includes detailed gas-phase halogen chemistry of chlorine and bromine. The stratospheric chemistry accounts for heterogeneous processes on liquid sulfate aerosols and polar stratospheric clouds, following the approach of Considine et al. (2000). The complete mechanism of MOZART (species and equations) can be found in the directory mecca/eqn/mozart/ in the Supplement. Overall, the MOZART mechanism includes 117 gas-phase species, 65 photolysis reactions, and 247 gas-phase reactions. Rate coefficients are taken from the JPL recommendations (Sander et al., 2006; S. P. Sander et al., 2011).

2.2.3 JAM002

Version 2 of the Jülich Atmospheric Mechanism (JAM002) has been implemented in the ECHAM-HAMMOZ chemistry–climate model (Schultz et al., 2018). It is a blend of the stratospheric chemistry scheme of the Whole Atmosphere Chemistry Climate Model (WACCM; Kinnison et al., 2007) and version 4 of the tropospheric MOZART model (Emmons et al., 2010). The combined chemistry scheme of WACCM and MOZART has been enhanced with a detailed representation of the oxidation of isoprene following the Mainz Isoprene Mechanism 2 (MIM2; Taraborrelli et al., 2009), and by adding a few primary volatile organic compounds and their oxidation chains. The isoprene oxidation scheme includes recent discoveries of 1,6 H-shift reactions (Peeters et al., 2009), the formation of epoxide (Paulot et al., 2009), and the photolysis of HPALDs (Wolfe et al., 2012). Some of the reaction products and rates were taken from MCM (Jenkin et al., 2015). Radical–radical reactions have been substantially revised since Emmons et al. (2010). In contrast to MCM, JAM002 does not use a radical pool but instead follows the pathways of peroxy radical reactions with HO₂, CH₃O₂, and CH₃COO₂ (peroxy acetyl) as explicitly as possible. Inorganic tropospheric chemistry considers ozone, NO, NO₂, NO₃, N₂O₅, HONO, HNO₃, HNO₄, HCN, CO, H₂, OH, HO₂, H₂O₂, NH₃, chlorine and bromine species, SO₂, and oxygen atoms. The complete mechanism of JAM002 (species and equations) can be found in the directory mecca/eqn/jam/ in the Supplement. In total, JAM002 contains 246 species and 733 reactions, including 142 photolysis reactions. Detailed information can be found in Schultz et al. (2018).

2.2.4 Master Chemical Mechanism (MCM)

The MCM describes in detail the tropospheric degradation of more than a hundred VOCs (Jenkin et al., 1997; Saunders et al., 2003). It is widely used as the reference mechanism for modeling studies of atmospheric processes. Although the standard organic chemistry mechanism in MECCA (MOM, described above) is sufficient for many model applications, a more explicit mechanism can be necessary when studying specific VOCs. For example, the fate of limonene (C₁₀H₁₆) emitted from boreal forests is not included in the standard MECCA mechanism. To use MCM reactions inside MECCA, the new tool xmcm2mecca has been added. It converts a subset extracted from the MCM web page (http://mcm.leeds.ac.uk/MCM; last access: 26 March 2019) to a KPP equation file that is compatible with MECCA. The user manual provides a detailed description of this new tool.

Table 3Simplified example list of species with overall interaction coefficients (OICs). The full mechanism includes all species; the skeletal mechanisms s1, s2, and s3 only include species above a certain OIC threshold. Targets with OIC =1 are always included. The color coding of the skeletal mechanism (used also in Fig. 4) shows in which mechanism a species occurs. For example, orange is used for species that are included in the full mechanism and in s1 but not in s2 and s3.

Download Print Version

https://www.geosci-model-dev.net/12/1365/2019/gmd-12-1365-2019-f04

Figure 4Skeletal reduction of terpene chemistry in the MOM reaction scheme (only C₁₀ species are shown here). Vertex colors and OIC values correspond to those in Table 3; only the green and yellow species are kept in the reduced mechanism.

Download

2.3 Mechanism intercomparison

Having several mechanisms implemented in the same modeling system enables mechanism intercomparison studies under exactly the same conditions. This approach ensures that any resulting differences come from the chemical mechanism, not from any other parts of the model. We have performed such an intercomparison for MOM, CB05BASCOE, MOZART, MIM1, JAM002, and a comparable subset of MCM. Details about these model runs and the results for all species are available in the testsuite/cams directory in the Supplement. Some representative results are shown in Fig. 3. All mechanisms show a very similar decay of the initial isoprene because they all use similar rate constants for the main reactions of isoprene with ozone, OH and NO₃. In contrast, the results for the terpenes differ. In CB05BASCOE and MOZART, the rate constants for the lumped terpenes are taken from α-pinene. In the other mechanisms, β-pinene (and other terpenes) are considered individually. Since β-pinene reacts with ozone much slower than α-pinene, the explicit treatment of β-pinene in the mechanism leads to a slower decay of the terpenes than in the lumped mechanisms. With respect to peroxyacetyl nitrate (PAN), CB05BASCOE especially shows very different values during the first day of the simulation. The calculated diurnal cycles of ozone, OH, and NO₂ are similar for all mechanisms but their absolute values vary; highest concentrations are produced by MIM1 and MCM, and the lowest by CB05BASCOE and JAM002. MOM and MOZART are in between.

It is interesting to compare our results to a mechanism intercomparison study conducted about 10 years ago by Emmerson and Evans (2009), who partially used predecessors of the mechanisms in our code. They found significant differences for both PAN and isoprene. Using the present-day versions of the mechanisms, we still see differences for PAN but very similar results for isoprene.

2.4 Skeletal mechanism reduction

In the area of fuel combustion research, chemical models require highly complex mechanisms to describe ignition, flame propagation, and other properties. In order to save computing time, several methods have been developed to create a simplified chemical mechanism (called skeletal mechanism), which still produce similar results as the full mechanism (e.g., Tomlin and Turányi, 2013). One of these methods is DRGEP (Directed Relation Graph with Error Propagation), which was introduced by Pepiot-Desjardins and Pitsch (2008) and implemented into the MARS (Mechanism Automatic Reduction Software) model by Niemeyer et al. (2010) and Niemeyer and Sung (2011). The DRGEP code from MARS has been implemented in CAABA/MECCA, making the skeletal reduction method available for atmospheric chemistry mechanisms. The most important quantities of DRGEP are briefly explained below; full details can be found in Niemeyer et al. (2010).

Targets:

Important chemical species, for which the skeletal mechanism has to produce similar results as the full mechanism.

Sample points:

A set of environmental conditions (temperature, pressure, concentrations of chemical species) simulated by the chemistry model.

Interaction coefficients (DIC, PIC, OIC):

The importance of chemical species in a mechanism is defined in terms of several interaction coefficients. The direct interaction coefficient (DIC) describes the importance of one species for another, based on its normalized contribution to production/consumption rates through reactions involving both species. Then, a graph search calculates a path interaction coefficient (PIC) based on the product of direct interaction coefficients along the path from target to species, where nodes represent species and weighted directed edges represent DICs. Finally, the overall interaction coefficient (OIC) is the maximum of all PICs between target and species. It is calculated for all sample points and expressed as a value between 0 (unimportant) and 1 (important). For targets, OIC =1 by definition. OIC values are only calculated for the full mechanism.

Error δ_skel:

A normalized value describing the error when using a skeletal mechanism instead of the full mechanism. A skeletal mechanism is suitable if δ_skel<1 for all targets and sample points. To allow individual weighting, the calculation of δ_skel depends on a target threshold AbsTol and a maximum acceptable relative tolerance RelTol, which are defined for all targets:

\begin{matrix} (5) & δ_{skel} = |\frac{max (x_{skel}, AbsTol)}{max (x_{full}, AbsTol)} - 1| / RelTol, \end{matrix}

where x_full and x_skel are the mixing ratios calculated with the full and the skeletal mechanism, respectively.

OIC threshold ε_ep:

A chemical species is considered important if OIC(species) >ε_ep. The final ε_ep calculated by DRGEP is the maximum value for which δ_skel<1 still holds.

To test the skeletal mechanism generation, we chose HCHO, HO₂, NO, O₃, and OH as targets, allowing a relative tolerance of RelTol=20 % for mixing ratios above a threshold of AbsTol=1 pmol mol⁻¹. Sample points were extracted from a global atmospheric chemistry simulation with a setup similar to that presented by Lelieveld et al. (2016). The chemical compositions were taken from several boxes at two altitudes (at the surface and at about 1 km). As we want the skeletal mechanism to perform well not only at typical concentrations of the targets but also when they are very high or very low, we picked boxes where the targets reach their minimum, average, or maximum concentrations, respectively. This resulted in the generation of 30 sample points (5 targets times (min/ave/max) times 2 altitudes), covering a wide range of values. The full mechanism contained the complete set of species from MOM (Sect. 2.1). To illustrate the mechanism, the subset describing terpene chemistry is shown in Fig. 4. The importance (OIC values) of a few selected species is shown in Table 3. Three skeletal mechanisms (s1, s2, s3) were generated, reducing the number of species from 663 in the full mechanism to 462, 429, and 411, respectively. The number of reactions was reduced from 2091 to 1444, 1320, and 1262, respectively. The third skeletal mechanism (s3) was rejected because it did not fulfill the criterion δ_skel<1. Results obtained with the full mechanism and with s2 were compared in a global simulation, as described below in Sect. 4.2.

2.5 Kinetic and isotope tagging

We have updated the sub-submodel MECCA-TAG (Gromov et al., 2010), which had been introduced in version 3.0 of CAABA. Several improvements to the kinetic tagging technique were implemented. These new features include the following.

Selectable composition transfer mode. Depending on the research question, prescribed-, molecular-, or element-weighted composition transfer may be selected. These modes determine the shares with which each reactant contributes to the products in the tagged chemical reactions: according to user-specified weightings, proportional to the reacting molecules count, or following the given element (e.g., C or H) content, respectively. Whilst the latter mode is intrinsic to isotope tagging, the others may be used for custom tagging configurations, e.g., product yield calculations.
Diagnostics for unaccounted production or loss of elemental composition. MECCA-TAG optionally adds passive diagnostic species to the tagged reactions with unbalanced transfer of the element of interest. This helps to quantify the amount of atoms the chemical mechanism receives from or loses to “nothing”, including the isotope composition of such mass-balance violations.
The new “class shifting” tagging mode. This mode allows for migration of molecules between the tagging classes in specified reactions, which allows for quantifying various exchange processes in the mechanism. For instance, one can distinguish oxidation generations: in reactions with given oxidants the products become “promoted” to the tagging class of the next oxidation generation. Another application of “class shifting” is quantifying the efficiency of recycling chains. In essence, such is the “online” implementation of the approach similar to that of Lehmann (2004), with the number of tagging classes defining the maximum of the recycling sequences it is possible to follow.

The range of MECCA-TAG applications was extended with new tagging configurations.

Radiocarbon configurations. These facilitate simulating the ¹⁴C content in a desired set of species, including the routines for calculating abundances using conventional units like pMC (percent Modern Carbon).
Hydrogen isotope chemistry. Now MECCA-TAG allows for tracing pathways of H transfer between the species in the mechanism. Furthermore, D∕H isotope chemistry (including relevant kinetic isotope effects for HO_x and C₁–C₂ chemistry) are included. The configuration and calculations of the composition transfer were extended with the possibility to specify isotope branching ratios necessary for the consistent D∕H kinetics simulations. Both H transfer and D∕H chemistry are currently evaluated in stratospheric setups of CAABA (Frank et al., 2018).
O₂ clumped isotope chemistry. Simulation of nonstochastic distributions of ¹⁸O¹⁸O and ¹⁷O¹⁸O isotopologues (Δ₃₆ and Δ₃₅ signatures) resulting from O(³P)-mediated temperature-dependent isotope exchange kinetics.

There are also some changes in the implementation and software requirements. There is no “doubling” mode anymore for evaluating the results of the optimized tagging. Performing kinetic tagging of the chemical mechanism with MECCA-TAG requires the Free Pascal Compiler (fpc, https://www.freepascal.org/, last access: 26 March 2019, version ≥2.6) at the time the xmecca script is run. The sub-submodel files are located in the mecca/tag/ directory of the distribution. The directory mecca/tag/cfg/ contains tagging configuration control files (*.cfg). The option to tag a newly created chemical mechanism is available in the xmecca script (also via batch files). Further details about the MECCA-TAG code development can be found in the file mecca/tag/CHANGELOG within the CAABA distribution.

3 Photolysis

CAABA contains several submodels that provide photolysis rate coefficients j, also called “j values”. The simple submodels READJ and SAPPHO have already been described by R. Sander et al. (2011). READJ has not changed since version 3.0. SAPPHO photolysis rates can now be scaled using a common enhancement factor “efact” for all photolysis rates. This has, for instance, been used to simulate the very bright conditions within a cloud top (Heue et al., 2014). The updated and new photolysis submodels JVAL and RADJIMT are described in the sections below.

3.1 JVAL

The submodel JVAL inside the CAABA/MECCA model calculates j values using the method of Landgraf and Crutzen (1998). It was first updated to the version described by Sander et al. (2014) and then additional changes were made. Many new photolysis reactions have been added, most of them related to either species from the MOM mechanism (CH3NO3, CH3O2NO2, CH3ONO, CH3O2, HCOOH, C2H5NO3, NOA, MEKNO3, BENZAL, HOC6H4NO2, CH3COCO2H, IPRCHO2HCO, C2H5CHO2HCO, C3H7CHO2HCO, PeDIONE24, PINAL2HCO) or organic halogen compounds (CF2ClCFCl2, CH3CFCl2, CF3CF2Cl, CF2ClCF2Cl, CHF2Cl, CHCl3, CH2Cl2). Besides, bug fixes were necessary regarding incorrect temperature dependencies of the ozone and OCS cross sections in the input data.

3.2 RADJIMT

RADJIMT is a new submodel that provides dissociation and ionization rates due to absorption of light and energetic photoelectrons in the mesosphere and thermosphere (see Table 4). It is part of the upper atmosphere extension of MESSy initially described by Baumgaertner et al. (2013), which was partly based on the implementations from the middle and upper atmosphere model CMAT2 (Harris, 2001; Dobbin, 2005; Dobbin and Aylward, 2008). For upper atmosphere simulations with CAABA, MECCA was extended by the relevant chemical species (electrons and ions) and reactions (labeled %Up in gas.eqn). For the respective literature sources, see meccanism.pdf in the Supplement.

Table 4New upper atmosphere reactions for which RADJIMT provides j values.

Download Print Version | Download XLSX

Photodissociation and photoionization due to the absorption of solar X-ray, EUV (extreme ultraviolet), and UV radiation are calculated using fluxes from the SOLAR2000 empirical model (Tobiska et al., 2000), the GLOW model (Solomon et al., 1988), and data presented by Henke et al. (1993) and Fennelly and Torr (1992). Relative partitioning between the possible products of the ionization process are based on the model of Strickland and Meier (1982) and Fuller-Rowell (1993).

For solar zenith angles larger than 75^∘, the atmospheric column of each absorbing species is calculated using an approximation of the Chapman grazing incidence function (Smith III and Smith, 1972).

Reaction enthalpies in kJ mol⁻¹ (exothermic chemical heating) are provided as a product of the relevant chemical reactions when “set enthalpy = y” is defined in the MECCA batch file. Radiative heating and cooling is also calculated by the submodel (variable “heatrates”).

As an example, we have performed simulations with CAABA using the MECCA and RADJIMT submodels. The mechanism was created using the batch file mtchem.bat, which selects reactions of the upper atmosphere labeled %Up. The model setup in caaba_mtchem.nml was used: the temperature was kept constant at 195 K and the pressure was set to 0.5 Pa (approximately 85 km). The model starts on 1 January. Chemical species were initialized using the values provided by Brasseur and Solomon (2005) in their Tables A.6.1 and A.6.2. The default time step length of 20 min was used. For MECCA and RADJIMT, the default settings were used. Model-calculated mixing ratios for a few selected species are shown in Fig. 5. A comprehensive set of plots is available in radjimt_mixrat.pdf and radjimt_jvalues.pdf in the Supplement.

https://www.geosci-model-dev.net/12/1365/2019/gmd-12-1365-2019-f05

Figure 5Model-calculated mixing ratios from an upper atmosphere simulation with MECCA and RADJIMT: diurnal cycles for 4 January (after 3 days of spinup) for the Equator (black) and a latitude of 50^∘ N (red). Time is in hours with local noon at 12. See Sect. 3.2 for further details.

Download

3.3 DISSOC

The new MESSy submodel DISSOC is based on the photolysis scheme by Meier et al. (1982). Briefly, it calculates a table of the so-called enhancement factor, which is basically the ratio of the actinic flux at a specific location to the solar irradiance at the top of the atmosphere. The enhancement factor depends on the pressure level, solar zenith angle, and wavelength. Input data are the solar irradiance at the top of the atmosphere, absorption cross sections, and ozone and oxygen profiles. For the implementation into global models, the input profiles are allowed to be latitude dependent, which increases the dimensions of the enhancement factor table from 3 to 4. Photolysis rates are calculated from the tabulated enhancement factor as a wavelength integral over the product with the absorption cross sections. The calculation is formulated in spherical geometry, such that it can be also applied to zenith angles above 90^∘. Rayleigh scattering is calculated based on Nicolet et al. (1982). Absorption cross sections are taken from the current JPL recommendations (Burkholder et al., 2015).

The code was first implemented by Lary and Pyle (1991) and coupled to a stratospheric chemistry-box model (Müller et al., 1994). Becker et al. (2000) improved the treatment of the diffuse actinic flux and corrected an implementation error of Meier et al. (1982). The extension to the use of multiple latitudes was introduced within the development of the model CLaMS (McKenna et al., 2002). The possibility to calculate diurnally averaged photolysis rates was introduced for the simplified fast chemistry setup used in multi-annual CLaMS simulations (Pommrich et al., 2014).

In the current configuration, DISSOC determines the photolysis rates for 38 photolysis reactions that are primarily of relevance in stratospheric chemistry. A standard setup contains 36 pressure levels, 18 latitude bins, and 28 solar zenith angle bins (of which 8 are above 90^∘). Of the 203 standard wavelength intervals between 116 and 850 nm, typically only the 159 intervals above 175 nm are used for tropospheric and stratospheric applications.

4 MECCA in the MESSy modeling system

Apart from using MECCA inside the CAABA box model, it is also possible to connect MECCA chemistry to a trajectory or global 3-dimensional model via the MESSy infrastructure (Jöckel et al., 2010, 2016). Recent developments of MECCA shown in this section are related to its implementation inside MESSy.

4.1 TRAJECT

The TRAJECT submodel by Riede et al. (2009) allows for simulations of atmospheric chemistry along precalculated Lagrangian trajectories. For this purpose, the air parcel simulated by CAABA is moved through space and time along a trajectory taken from an external input file, while simulating atmospheric photochemistry with MECCA and JVAL. More generally, TRAJECT allows us to prescribe physical boundary conditions for CAABA box model simulations. A typical application is the simulation of atmospheric trajectories (balloon measurements or backward trajectories). However, laboratory conditions (e.g., in a flow reactor) can also be prescribed. The previous TRAJECT version, described by R. Sander et al. (2011), has been updated. The output is now more consistent with the trajectory input file, as physical information is now written out beginning with the first time step instead of the second. In general, an integration time step of chemical kinetics is always performed with the physical parameters given for the end of the time step. In that way, the mixing ratios written out at the end of a time step are consistent with the physical conditions at that point. Also, solar zenith angle and local time at the end of a time step are now consistent with the given longitude and latitude for that trajectory point.

In addition to the trajectory input file, an external input file with j values for NO₂ can be used to scale all j values with the factor

\begin{matrix} (6) & jfac = \frac{j ({NO}_{2}, external)}{j ({NO}_{2}, JVAL)} . \end{matrix}

To facilitate the analysis of the scaling impact, jfac is now written to output. Scaling thresholds have been implemented to prevent artifacts that would occur when j(NO₂,JVAL) is very small and the calculation of jfac approaches a division by zero.

4.2 PolyMECCA/CHEMGLUE

In a standard global model simulation, the MESSy submodel MECCA contains one chemical mechanism that is used for all grid boxes. This ensures a consistent chemistry simulation from the surface to the upper atmosphere. However, in some cases, it may be preferable to allow different mechanisms in different boxes, e.g., terpene chemistry only in the troposphere and ion chemistry only in the mesosphere.

With the script xpolymecca, several independent chemical MECCA mechanisms can be produced. The first mechanism has the name “mecca”, as usual. Additional mechanisms are labeled with a three-digit suffix. For example, the code of mechanism 2 is contained in messy_mecca002_kpp.f90 and related files.

To select an appropriate mechanism at each point in space and time, the MESSy submodel CHEMGLUE has been written. The name of the submodel was chosen because CHEMGLUE can also glue together different chemical mechanisms at the border where a chemical species is included in one mechanism but not in the other. CHEMGLUE defines the new channel object “meccanum”, which contains the mechanism number for each grid point. These values can either be selected statically, e.g., depending on the model level number or the sea–land fraction mask. Alternatively, a dynamic (time-dependent) selection based on chemical or meteorological variables is possible, e.g., pressure, temperature, or the concentrations of ozone or isoprene.

Note that even when different boxes of a global model simulation use different chemistry mechanisms, the set of tracers contains all species from all mechanisms for all boxes.

The implementation ensures binary identical results when one chemical mechanism (“mecca”) is replaced by two identical copies of it (“mecca” and “mecca002”).

https://www.geosci-model-dev.net/12/1365/2019/gmd-12-1365-2019-f06

Figure 6Results of the global comparison between the FULL, POLY, and SKEL mechanisms (see Sect. 4.2 for details). Shown are surface mixing ratios of ozone (a, c, e) and isoprene (b, d, f) at the end of the simulation, i.e., after one month. Panels (a) and (b) show results obtained with the FULL chemistry mechanism. Panels (c) and (d) compare POLY to FULL, and panels (e) and (f) compare SKEL to FULL.

Download

For a more realistic test, we created two different chemical mechanisms for organics. In the first mechanism, only the oxidation of methane is considered, and all nonmethane hydrocarbons are neglected. The second (FULL) contains the full set of MOM (Sect. 2.1) reactions. CHEMGLUE selects the second mechanism whenever the mixing ratios of organics are above a threshold (isoprene >100 pmol mol⁻¹, α-pinene >100 pmol mol⁻¹, or toluene >10 pmol mol⁻¹). To investigate how much CPU time can be saved and how much the simplification affects the results, we have performed global test simulations based on the ECHAM5/MESSy atmospheric chemistry (EMAC) model by Jöckel et al. (2016). The horizontal resolution was T42 ( ${2.8}^{\circ} \times {2.8}^{\circ}$ ), with 47 vertical levels. Starting on 1 January 2009, 1 month was simulated. To facilitate the intercomparison between the simulations, the feedback of chemistry on the meteorology was switched off. Three different chemical schemes were tested.

FULL. Full MOM chemistry was activated throughout the atmosphere.
POLY. PolyMECCA/CHEMGLUE switches between the full MOM chemistry and the methane-only chemistry as described above.
SKEL. The skeletal mechanism s2 as described in Sect. 2.4 was activated throughout the atmosphere.

The CPU usage for the POLY and SKEL simulations are 62 % and 65 % of the FULL simulation, respectively. Results are shown in Fig. 6. Overall, the agreement between the simulations is quite good, considering that the simplified mechanisms neglect many reactions.

4.3 CHEMPROP

Chemical properties of the species in the reaction mechanism are needed at many locations in the model, e.g., molar masses (M), Henry's law constants (H), accommodation coefficients (α), acidity constants (K_A), and ion charge numbers (z). These values have so far been stored at different locations in the code (gas.tex, messy_cmn_gasaq.f90, and elsewhere). Because maintaining data that are spread over several source files is tedious and error-prone, the new CHEMPROP database has been created, which stores all values centrally in the ASCII table messy_main_tracer_chemprop.tbl. MECCA (and other submodels) can access these chemical property data via MESSy tracer containers, as described by Jöckel et al. (2008).

5 Further changes

The new subroutines dilute and dilute_once dilute the concentrations of chemicals in an air parcel by mixing it with unperturbed air. This can, for example, be used for modeling chemistry in an expanding volcanic plume or smog plume. An alternative usage for these subroutines is the simulation of the flow in and out of a reaction chamber (e.g., van Eijck et al., 2013).
A new functionality has been implemented for the external initialization of chemical species from a netCDF file: if the time axis of the input file contains more than one point, the time values are used to interpolate mixing ratios at model start time. This is convenient for bundling several initializations into one file, for instance to initialize several CAABA simulations from different points along a trajectory with recorded mixing ratios (see also Sect. 4.1). If the time axis of the input file contains only one point, the mixing ratios are read into CAABA, regardless of the time value.
We extended CAABA with parameters to optionally control the output step frequency (output_step_freq) and the output synchronization frequency (output_sync_freq). The first variable sets the frequency at which values are written to the output. A value of output_step_freq =α skips α−1 time steps and writes only every αth time step to the output file. The second variable controls the output synchronization. Data are buffered for output_sync_freq time steps before they are written to the output files. Both parameters enable the user to carry out very long box model simulations without being constrained by machine I/O performance, and they can individually regulate the output file size. A high value of output_sync_freq has a positive effect on performance. However, in case of machine failure, buffered output steps are lost.
The treatment of humidity has been improved. Now specific as well as relative humidity (RH) is available throughout CAABA, and can be interconverted with generic conversion functions. Of the two, specific humidity is the more robust variable for humidity because the definition of RH can be based on either partial pressure or on specific humidity (Jacobson, 1999). There are various parameterizations for saturation water vapor pressure, and RH can be defined over liquid surface even below 0 ^∘C if supercooling is allowed. Functions that use humidity as input (concentration of air, conversion between humidity, and water vapor concentration) now use the unambiguous specific humidity. If necessary, it is derived from relative humidity taking all of the above considerations into account.
For better model time control, two boolean namelist parameters have been introduced: l_groundhogday=T repeats a diurnal cycle while l_freezetime=T repeats a certain point in time, effectively freezing the solar zenith angle.
The selection of various chemical species to define steady state has been simplified to allow for more flexibility in the criteria. The progress towards the defined steady state is now logged during CAABA runtime. Artifacts by species' concentrations close to zero are now prevented.
Several shell scripts have been converted to python (xcaaba.py, multirun.py, montecarlo.py). They use the netcdf4 interface and don't depend on the availability of the NetCDF operators (e.g., ncks) anymore. Currently, the python scripts are in beta testing. In future versions, they will replace the current tcsh scripts.
Model results can now be visualized with the python script caabaplot.py using matplotlib. The previously used ferret scripts are still included but not actively supported anymore.
Complex reaction mechanisms can be interpreted as graphs, with species representing vertices and reactions representing edges. To visualize and analyze these graphs, the “graph-tool” software by Peixoto (2014) can now be used. For example, Fig. 4 was created with graph-tool.
Rate coefficients have been updated to the latest JPL recommendations (Burkholder et al., 2015) and recent laboratory studies. A complete list of chemical reactions, rate coefficients, and references is available in the Supplement (meccanism.pdf).
The kinetic preprocessor KPP (Sandu and Sander, 2006) performs the numerical integration of the chemical reaction mechanism. It has been updated to the latest version 2.2.3, which contains a number of small fixes throughout the code (http://people.cs.vt.edu/~asandu/Software/Kpp, last access: 26 March 2019).
The python scripts check_eqntags.py and check_eqns.pl check the internal consistency of the chemical mechanism.
Details of all new features have been added to the updated user manual, which now also includes an index. Additional minor bug fixes can be found in the CHANGELOG file.

6 Summary and outlook

We have presented the current version of the atmospheric chemistry module MECCA-4.0, which includes several new features: skeletal mechanism reduction, the MOM chemical mechanism for organic compounds, optional inclusion of reactions from MCM and other chemical mechanisms, updated isotope tagging, and improved and new photolysis modules. When MECCA is connected to a global model, PolyMECCA and CHEMGLUE allow coexisting multiple chemistry mechanisms. CAABA/MECCA is now available to the research community.

Based on the model development described in this paper, our current and upcoming goals are the following. (For work in progress, initials of the principal investigators are shown in parentheses.)

reduce complex mechanisms to a size suitable for global model simulations (RS, KEN),
perform a chemistry module intercomparison including CB05BASCOE and MOZART within a global chemistry modeling framework (Huijnen et al., 2019),
evaluate MOM chemistry and its effect on secondary aerosol formation (AP),
compare MOM chemistry to measurements obtained during the recent AQABA field campaign (HH),
advance our understanding of the role of organic compounds on the tropospheric ozone and HO_x budgets (DT),
compare model results with studies at the SAPHIR chamber (DT),
investigate the multiphase chemical pathways leading to organic acids and aerosols (DT),
simulate stratospheric isotope H exchanges between CH₄ and H₂O (SG),
implement additional photolysis modules (e.g., CLOUDJ, TUV) and compare the resulting j values (HH),
parallelize to distribute independent (e.g., Monte Carlo or sensitivity) box model simulations on multiple cores (HH),
study the impact of aromatic compounds on atmospheric chemistry (RS, manuscript in preparation).

Code and data availability

The CAABA/MECCA model code is available as a community model published under the GNU General Public License (http://www.gnu.org/copyleft/gpl.html, last access: 26 March 2019). The model code can be found in the Supplement. In addition to the complete code, a list of chemical reactions – including rate coefficients and references (meccanism.pdf) – and a user manual (caaba_mecca_manual.pdf) are available in the manual directory of the Supplement. For further information and updates, the MECCA web page at http://www.mecca.messy-interface.org (last access: 26 March 2019) can be consulted.

Supplement

The supplement related to this article is available online at: https://doi.org/10.5194/gmd-12-1365-2019-supplement.

Author contributions

RS develops and maintains the CAABA/MECCA software. AB provided RADJIMT. DCP provided the aromatic chemistry mechanism of MOM. FF added code to control the model output. JUG provided DISSOC and helped with its implementation in MESSy. SG provided the MECCA-TAG sub-submodel. HH and ST provided code for the inclusion of the MCM reaction schemes. PJ contributed to several model development projects (MESSy modeling system, PolyMECCA, CHEMPROP, CHEMGLUE) and maintains the interfaces to ensure that the modules are not only compatible with the box model but also with the 3-D models. VH contributed through initiating the provision of CAMS chemistry models for inclusion in MECCA, and for generation of the CB05BASCOE merged chemical mechanism. VAK integrated CB05BASCOE and MOZART into MECCA. KEN provided code for the skeletal mechanism generation. AP contributed to several model development projects (MESSy modeling system, scenarios for skeletal mechanism generation, MOM, CAMS, PolyMECCA testing). HR provided an update of the TRAJECT submodel. MGS provided JAM002, and DT provided MOM.

Competing interests

The authors declare that they have no conflict of interest.

Acknowledgements

We thank Simon Chabrillat and Idir Bouarar for provision of the mechanisms BASCOE and MOZART, respectively. We also thank Tim Butler for contributing the diagnostic tool check_eqns.pl. Duy Cai added some photolysis reactions to JVAL. Tilo Fytterer and Stefan Versick discovered and reported the temperature dependence bug of the ozone and OCS photolyses in JVAL. Vincent Huijnen acknowledges funding from the Copernicus Atmosphere Monitoring Service (CAMS). Model results were plotted using the open-source software matplotlib (https://matplotlib.org, last access: 26 March 2019) and ferret (http://ferret.pmel.noaa.gov, last access: 26 March 2019).

The article processing charges for this open-access
publication were covered by the Max Planck Society.

Review statement

This paper was edited by Slimane Bekki and reviewed by two anonymous referees.

References

Arey, J., Aschmann, S. M., Kwok, E. S. C., and Atkinson, R.: Alkyl nitrate, hydroxyalkyl nitrate, and hydroxycarbonyl formation from the NO_x-air photooxidations of C₅-C₈ n-alkanes, J. Phys. Chem. A, 105, 1020–1027, https://doi.org/10.1021/jp003292z, 2001. a

Assaf, E., Song, B., Tomas, A., Schoemaecker, C., and Fittschen, C.: Rate constant of the reaction between CH₃O₂ radicals and OH radicals revisited, J. Phys. Chem. A, 120, 8923–8932, https://doi.org/10.1021/acs.jpca.6b07704, 2016. a

Assaf, E., Sheps, L., Whalley, L., Heard, D., Tomas, A., Schoemaecker, C., and Fittschen, C.: The reaction between CH₃O₂ radicals and OH radicals: Product yields and atmospheric implications, Environ. Sci. Technol., 51, 2170–2177, https://doi.org/10.1021/acs.est.6b06265, 2017. a

Atkinson, R.: A structure-activity relationship for the estimation of rate constants for the gas-phase reactions of OH radicals with organic compounds, Int. J. Chem. Kinet., 19, 799–828, https://doi.org/10.1002/kin.550190903, 1987. a

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, https://doi.org/10.5194/acp-6-3625-2006, 2006. a, b, c, d

Barnes, I., Becker, K. H., and Zhu, T.: Near UV absorption-spectra and photolysis products of difunctional organic nitrates – possible importance as NO_x reservoirs, J. Atmos. Chem., 17, 353–373, https://doi.org/10.1007/BF00696854, 1993. a

Baumgaertner, A. J. G., Jöckel, P., Aylward, A. D., and Harris, M. J.: Simulation of particle precipitation effects on the atmosphere with the MESSy model system, in: Climate and Weather of the Sun-Earth System (CAWSES), Highlights from a Priority Program, edited by: Lübken, F.-J., 301–316, Springer Verlag, Berlin, 2013. a

Becker, G., Grooß, J.-U., McKenna, D. S., and Müller, R.: Stratospheric photolysis frequencies: Impact of an improved numerical solution of the radiative transfer equation, J. Atmos. Chem., 37, 217–229, https://doi.org/10.1023/A:1006468926530, 2000. a

Bejan, I., Abd El Aal, Y., Barnes, I., Benter, T., Bohn, B., Wiesen, P., and Kleffmann, J.: The photolysis of ortho-nitrophenols: a new gas phase source of HONO, Phys. Chem. Chem. Phys., 8, 2028–2035, https://doi.org/10.1039/B516590C, 2006. a, b

Birdsall, A. W., Andreoni, J. F., and Elrod, M. J.: Investigation of the role of bicyclic peroxy radicals in the oxidation mechanism of toluene, J. Phys. Chem. A, 114, 10655–10663, https://doi.org/10.1021/jp105467e, 2010. a

Bloss, C., Wagner, V., Jenkin, M. E., Volkamer, R., Bloss, W. J., Lee, J. D., Heard, D. E., Wirtz, K., Martin-Reviejo, M., Rea, G., Wenger, J. C., and Pilling, M. J.: Development of a detailed chemical mechanism (MCMv3.1) for the atmospheric oxidation of aromatic hydrocarbons, Atmos. Chem. Phys., 5, 641–664, https://doi.org/10.5194/acp-5-641-2005, 2005. a

Bossolasco, A., Faragó, E. P., Schoemaecker, C., and Fittschen, C.: Rate constant of the reaction between CH₃O₂ and OH radicals, Chem. Phys. Lett., 593, 7–13, https://doi.org/10.1016/j.cplett.2013.12.052, 2014. a, b

Boyd, A. A., Flaud, P.-M., Daugey, N., and Lesclaux, R.: Rate constants for RO₂+HO₂ reactions measured under a large excess of HO₂, J. Phys. Chem. A, 107, 818–821, https://doi.org/10.1021/JP026581R, 2003. a

Brasseur, G., Smith, A., Khosravi, R., Huang, T., Walters, S., Chabrillat, S., and Kockarts, G.: Natural and human-induced perturbations in the middle atmosphere: A short tutorial, in: Atmospheric Science Across the Stratopause, Geophysical Monograph 123, edited by: Siskind, D. E., Eckermann, S. D., and Summers, M. E., Am. Geophys. Union, Washington, 7–20, https://doi.org/10.1029/GM123p0007, 2000. a

Brasseur, G. P. and Solomon, S.: Aeronomy of the Middle Atmosphere: Chemistry and Physics of the Stratosphere and Mesosphere, Springer Verlag, https://doi.org/10.1007/1-4020-3824-0, 2005. a

Burkholder, J. B., Sander, S. P., Abbatt, J., Barker, J. R., Huie, R. E., Kolb, C. E., Kurylo, M. J., Orkin, V. L., Wilmouth, D. M., and Wine, P. H.: Chemical Kinetics and Photochemical Data for Use in Atmospheric Studies, Evaluation No. 18, JPL Publication 15-10, Jet Propulsion Laboratory, Pasadena, available at: http://jpldataeval.jpl.nasa.gov (last access: 26 March 2019), 2015. a, b, c, d, e

Butkovskaya, N., Kukui, A., and Le Bras, G.: Pressure and temperature dependence of ethyl nitrate formation in the C₂H₅O₂+NO reaction, J. Phys. Chem. A, 114, 956–964, https://doi.org/10.1021/jp910003a, 2010. a

Butkovskaya, N., Kukui, A., and Le Bras, G.: Pressure and temperature dependence of methyl nitrate formation in the CH₃O₂ + NO reaction, J. Phys. Chem. A, 116, 5972–5980, https://doi.org/10.1021/jp210710d, 2012. a

Cabrera-Perez, D., Taraborrelli, D., Sander, R., and Pozzer, A.: Global atmospheric budget of simple monocyclic aromatic compounds, Atmos. Chem. Phys., 16, 6931–6947, https://doi.org/10.5194/acp-16-6931-2016, 2016. a, b

Capouet, M., Müller, J.-F., Ceulemans, K., Compernolle, S., Vereecken, L., and Peeters, J.: Modeling aerosol formation in alpha-pinene photo-oxidation experiments, J. Geophys. Res., 113D, https://doi.org/10.1029/2007JD008995, 2008. a

Chabrillat, S. and Fonteyn, D.: Modelling long-term changes of mesospheric temperature and chemistry, Adv. Space Res., 32, 1689–1700, https://doi.org/10.1016/S0273-1177(03)90464-9, 2003. a

Chen, J., Wenger, J. C., and Venables, D. S.: Near-ultraviolet absorption cross sections of nitrophenols and their potential influence on tropospheric oxidation capacity, J. Phys. Chem. A, 115, 12235–12242, https://doi.org/10.1021/jp206929r, 2011. a, b

Cheng, S.-B., Zhou, C.-H., Yin, H.-M., Sun, J.-L., and Han, K.-L.: OH produced from o-nitrophenol photolysis: A combined experimental and theoretical investigation, J. Chem. Phys., 130, 234 311, https://doi.org/10.1063/1.3152635, 2009. a

Clubb, A. E., Jordan, M. J. T., Kable, S. H., and Osborn, D. L.: Phototautomerization of acetaldehyde to vinyl alcohol: a primary process in UV-irradiated acetaldehyde from 295 to 335 nm, J. Phys. Chem. Lett., 3, 3522–3526, https://doi.org/10.1021/jz301701x, 2012. a

Considine, D. B., Douglass, A. R., Connell, P. S., Kinnison, D. E., and Rotman, D. A.: A polar stratospheric cloud parameterization for the global modeling initiative three-dimensional model and its response to stratospheric aircraft, J. Geophys. Res., 105D, 3955–3973, https://doi.org/10.1029/1999JD900932, 2000. a

Crounse, J. D., Knap, H. C., Ørnsø, K. B., Jørgensen, S., Paulot, F., Kjaergaard, H. G., and Wennberg, P. O.: Atmospheric fate of methacrolein. 1. peroxy radical isomerization following addition of OH and O₂, J. Phys. Chem. A, 116, 5756–5762, https://doi.org/10.1021/jp211560u, 2012. a

da Silva, G.: Carboxylic acid catalyzed keto-enol tautomerizations in the gas phase, Angew. Chem.-Int. Ed., 122, 7685–7687, https://doi.org/10.1002/ange.201003530, 2010. a

Dillon, T. J. and Crowley, J. N.: Direct detection of OH formation in the reactions of HO₂ with CH₃C(O)O₂ and other substituted peroxy radicals, Atmos. Chem. Phys., 8, 4877–4889, https://doi.org/10.5194/acp-8-4877-2008, 2008. a

Dobbin, A. L.: Modelling studies of possible coupling mechanisms between the upper and middle atmosphere, PhD thesis, Department of Physics and Astronomy, University College London, available at: http://discovery.ucl.ac.uk/1444625 (last access: 26 March 2019), 2005. a

Dobbin, A. L. and Aylward, A. D.: A three-dimensional modelling study of the processes leading to mid latitude nitric oxide increases in the lower thermosphere following periods of high geomagnetic activity, Adv. Space Res., 42, 1576–1585, https://doi.org/10.1016/j.asr.2008.03.004, 2008. a

Elrod, M. J.: Kinetics study of the aromatic bicyclic peroxy radical + NO reaction: overall rate constant and nitrate product yield measurements, J. Phys. Chem. A, 115, 8125–8130, https://doi.org/10.1021/jp204308f, 2011. a

Emmerson, K. M. and Evans, M. J.: Comparison of tropospheric gas-phase chemistry schemes for use within global models, Atmos. Chem. Phys., 9, 1831–1845, https://doi.org/10.5194/acp-9-1831-2009, 2009. a

Emmons, L. K., Walters, S., Hess, P. G., Lamarque, J.-F., Pfister, G. G., Fillmore, D., Granier, C., Guenther, A., Kinnison, D., Laepple, T., Orlando, J., Tie, X., Tyndall, G., Wiedinmyer, C., Baughcum, S. L., and Kloster, S.: Description and evaluation of the Model for Ozone and Related chemical Tracers, version 4 (MOZART-4), Geosci. Model Dev., 3, 43–67, https://doi.org/10.5194/gmd-3-43-2010, 2010. a, b, c

Errera, Q. and Fonteyn, D.: Four-dimensional variational chemical assimilation of CRISTA stratospheric measurements, J. Geophys. Res., 106D, 12253–12265, https://doi.org/10.1029/2001JD900010, 2001. a

Errera, Q., Daerden, F., Chabrillat, S., Lambert, J. C., Lahoz, W. A., Viscardy, S., Bonjean, S., and Fonteyn, D.: 4D-Var assimilation of MIPAS chemical observations: ozone and nitrogen dioxide analyses, Atmos. Chem. Phys., 8, 6169–6187, https://doi.org/10.5194/acp-8-6169-2008, 2008. a

Fennelly, J. A. and Torr, D. G.: Photoionization and photoabsorption cross sections of O, N₂, O₂, and N for aeronomic calculations, Atom. Data Nucl. Data Tables, 51, 321–363, https://doi.org/10.1016/0092-640X(92)90004-2, 1992. a

Flocke, F., Atlas, E., Madronich, S., Schauffler, S. M., Aikin, K., Margitan, J. J., and Bui, T. P.: Observations of methyl nitrate in the lower stratosphere during STRAT: implications for its gas phase production mechanisms, Geophys. Res. Lett., 25, 1891–1894, https://doi.org/10.1029/98GL01417, 1998. a

Frank, F., Jöckel, P., Gromov, S., and Dameris, M.: Investigating the yield of H₂O and H₂ from methane oxidation in the stratosphere, Atmos. Chem. Phys., 18, 9955–9973, https://doi.org/10.5194/acp-18-9955-2018, 2018. a

Fuller-Rowell, T. J.: Modeling the solar cycle change in nitric oxide in the thermosphere and upper mesosphere, J. Geophys. Res., 98A, 1559–1570, https://doi.org/10.1029/92JA02201, 1993. a

Glover, B. G. and Miller, T. A.: Near-IR cavity ringdown spectroscopy and kinetics of the isomers and conformers of the butyl peroxy radical, J. Phys. Chem. A, 109, 11191–11197, https://doi.org/10.1021/jp054838q, 2005. a

Gromov, S., Jöckel, P., Sander, R., and Brenninkmeijer, C. A. M.: A kinetic chemistry tagging technique and its application to modelling the stable isotopic composition of atmospheric trace gases, Geosci. Model Dev., 3, 337–364, https://doi.org/10.5194/gmd-3-337-2010, 2010. a

Groß, C. B. M.: Kinetische Studien zur OH-Bildung über die Reaktionen von HO₂ mit organischen Peroxyradikalen, PhD thesis, Johannes Gutenberg-Universität, 2013. a, b

Groß, C. B. M., Dillon, T. J., Schuster, G., Lelieveld, J., and Crowley, J. N.: Direct kinetic study of OH and O₃ formation in the reaction of CH₃C(O)O₂ with HO₂, J. Phys. Chem. A, 1, 974–985, https://doi.org/10.1021/jp412380z, 2014. a

Harris, M. J.: A new coupled middle atmosphere and thermosphere general circulation model: Studies of dynamic, energetic and photochemical coupling in the middle and upper atmosphere, PhD thesis, Department of Physics and Astronomy, University College London, available at: http://search.proquest.com/docview/1758611641 (last access: 26 March 2019), 2001. a

Henke, B. L., Gullikson, E. M., and Davis, J. C.: X-Ray Interactions: Photoabsorption, Scattering, Transmission, and Reflection at E = 50–30,000 eV, Z = 1–92, Atom. Data Nucl. Data Tables, 54, 181–342,https://doi.org/10.1006/adnd.1993.1013, 1993. a

Hens, K., Novelli, A., Martinez, M., Auld, J., Axinte, R., Bohn, B., Fischer, H., Keronen, P., Kubistin, D., Nölscher, A. C., Oswald, R., Paasonen, P., Petäjä, T., Regelin, E., Sander, R., Sinha, V., Sipilä, M., Taraborrelli, D., Tatum Ernest, C., Williams, J., Lelieveld, J., and Harder, H.: Observation and modelling of HO_x radicals in a boreal forest, Atmos. Chem. Phys., 14, 8723–8747, https://doi.org/10.5194/acp-14-8723-2014, 2014. a

Heue, K.-P., Riede, H., Walter, D., Brenninkmeijer, C. A. M., Wagner, T., Frieß, U., Platt, U., Zahn, A., Stratmann, G., and Ziereis, H.: CARIBIC DOAS observations of nitrous acid and formaldehyde in a large convective cloud, Atmos. Chem. Phys., 14, 6621–6642, https://doi.org/10.5194/acp-14-6621-2014, 2014. a

Huijnen, V., Williams, J., van Weele, M., van Noije, T., Krol, M., Dentener, F., Segers, A., Houweling, S., Peters, W., de Laat, J., Boersma, F., Bergamaschi, P., van Velthoven, P., Le Sager, P., Eskes, H., Alkemade, F., Scheele, R., Nédélec, P., and Pätz, H.-W.: The global chemistry transport model TM5: description and evaluation of the tropospheric chemistry version 3.0, Geosci. Model Dev., 3, 445–473, https://doi.org/10.5194/gmd-3-445-2010, 2010. a

Huijnen, V., Flemming, J., Chabrillat, S., Errera, Q., Christophe, Y., Blechschmidt, A.-M., Richter, A., and Eskes, H.: C-IFS-CB05-BASCOE: stratospheric chemistry in the Integrated Forecasting System of ECMWF, Geosci. Model Dev., 9, 3071–3091, https://doi.org/10.5194/gmd-9-3071-2016, 2016. a, b

Huijnen, V., Pozzer, A., Arteta, J., Brasseur, G., Bouarar, I., Chabrillat, S., Christophe, Y., Doumbia, T., Flemming, J., Guth, J., Josse, B., Karydis, V. A., Marécal, V., and Pelletier, S.: Quantifying uncertainties due to chemistry modeling – evaluation of tropospheric composition simulations in the CAMS model, Geosci. Model Dev. Discuss., https://doi.org/10.5194/gmd-2018-331, in review, 2019. a

Jacobson, M. Z.: Fundamentals of atmospheric modeling, Cambridge University Press, Cambridge, 1999. a

Jagiella, S. and Zabel, F.: Reaction of phenylperoxy radicals with NO₂ at 298 K, Phys. Chem. Chem. Phys., 9, 5036–5051, https://doi.org/10.1039/B705193J, 2007. a

Jenkin, M., Saunders, S. M., and Pilling, M. J.: The tropospheric degradation of volatile organic compounds: A protocol for mechanism development, Atmos. Environ., 31, 81–104, https://doi.org/10.1016/S1352-2310(96)00105-7, 1997. a, b, c

Jenkin, M. E., Shallcross, D. E., and Harvey, J. N.: Development and application of a possible mechanism for the generation of cis-pinic acid from the ozonolysis of α- and β-pinene, Atmos. Environ., 34, 2837–2850, https://doi.org/10.1016/S1352-2310(00)00087-X, 2000. a

Jenkin, M. E., Saunders, S. M., Wagner, V., and Pilling, M. J.: Protocol for the development of the Master Chemical Mechanism, MCM v3 (Part B): tropospheric degradation of aromatic volatile organic compounds, Atmos. Chem. Phys., 3, 181–193, https://doi.org/10.5194/acp-3-181-2003, 2003. a

Jenkin, M. E., Young, J. C., and Rickard, A. R.: The MCM v3.3.1 degradation scheme for isoprene, Atmos. Chem. Phys., 15, 11433–11459, https://doi.org/10.5194/acp-15-11433-2015, 2015. a, b, c

Jöckel, P., Kerkweg, A., Buchholz-Dietsch, J., Tost, H., Sander, R., and Pozzer, A.: Technical Note: Coupling of chemical processes with the Modular Earth Submodel System (MESSy) submodel TRACER, Atmos. Chem. Phys., 8, 1677–1687, https://doi.org/10.5194/acp-8-1677-2008, 2008. a

Jöckel, P., Kerkweg, A., Pozzer, A., Sander, R., Tost, H., Riede, H., Baumgaertner, A., Gromov, S., and Kern, B.: Development cycle 2 of the Modular Earth Submodel System (MESSy2), Geosci. Model Dev., 3, 717–752, https://doi.org/10.5194/gmd-3-717-2010, 2010. a, b

Jöckel, P., Tost, H., Pozzer, A., Kunze, M., Kirner, O., Brenninkmeijer, C. A. M., Brinkop, S., Cai, D. S., Dyroff, C., Eckstein, J., Frank, F., Garny, H., Gottschaldt, K.-D., Graf, P., Grewe, V., Kerkweg, A., Kern, B., Matthes, S., Mertens, M., Meul, S., Neumaier, M., Nützel, M., Oberländer-Hayn, S., Ruhnke, R., Runde, T., Sander, R., Scharffe, D., and Zahn, A.: Earth System Chemistry integrated Modelling (ESCiMo) with the Modular Earth Submodel System (MESSy) version 2.51, Geosci. Model Dev., 9, 1153–1200, https://doi.org/10.5194/gmd-9-1153-2016, 2016. a, b, c

Kinnison, D. E., Brasseur, G. P., Walters, S., Garcia, R. R., Marsh, D. R., Sassi, F., Harvey, V. L., Randall, C. E., Emmons, L., Lamarque, J. F., Hess, P., Orlando, J. J., Tie, X. X., Randel, W., Pan, L. L., Gettelman, A., Granier, C., Diehl, T., Niemeier, U., and Simmons, A. J.: Sensitivity of chemical tracers to meteorological parameters in the MOZART-3 chemical transport model, J. Geophys. Res., 112D, https://doi.org/10.1029/2006JD007879, 2007. a, b, c

Kirchner, F., Mayer-Figge, A., Zabel, F., and Becker, K. H.: Thermal stability of peroxynitrates, Int. J. Chem. Kinet., 31, 127–144, https://doi.org/10.1002/(SICI)1097-4601(1999)31:2<127::AID-KIN6>3.0.CO;2-L, 1999. a

Kwok, E. S. C. and Atkinson, R.: Estimation of hydroxyl radical reaction rate constants for gas-phase organic compounds using a structure-reactivity relationship: an update, Atmos. Environ., 29, 1685–1695, https://doi.org/10.1016/1352-2310(95)00069-B, 1995. a, b, c, d, e

Lamarque, J.-F., Emmons, L. K., Hess, P. G., Kinnison, D. E., Tilmes, S., Vitt, F., Heald, C. L., Holland, E. A., Lauritzen, P. H., Neu, J., Orlando, J. J., Rasch, P. J., and Tyndall, G. K.: CAM-chem: description and evaluation of interactive atmospheric chemistry in the Community Earth System Model, Geosci. Model Dev., 5, 369–411, https://doi.org/10.5194/gmd-5-369-2012, 2012. a, b

Landgraf, J. and Crutzen, P. J.: An efficient method for online calculations of photolysis and heating rates, J. Atmos. Sci., 55, 863–878, https://doi.org/10.1175/1520-0469(1998)055<0863:AEMFOC>2.0.CO;2, 1998. a

Lary, D. J. and Pyle, J. A.: Diffuse radiation, twilight, and photochemistry – I, J. Atmos. Chem., 13, 373–406, https://doi.org/10.1007/BF00057753, 1991. a

Lehmann, R.: An algorithm for the determination of all significant pathways in chemical reaction systems, J. Atmos. Chem., 47, 45–78, https://doi.org/10.1023/B:JOCH.0000012284.28801.B1, 2004. a

Lelieveld, J., Gromov, S., Pozzer, A., and Taraborrelli, D.: Global tropospheric hydroxyl distribution, budget and reactivity, Atmos. Chem. Phys., 16, 12477–12493, https://doi.org/10.5194/acp-16-12477-2016, 2016. a, b

Liu, Y. J., Herdlinger-Blatt, I., McKinney, K. A., and Martin, S. T.: Production of methyl vinyl ketone and methacrolein via the hydroperoxyl pathway of isoprene oxidation, Atmos. Chem. Phys., 13, 5715–5730, https://doi.org/10.5194/acp-13-5715-2013, 2013. a

Madronich, S. and Calvert, J. G.: Permutation reactions of organic peroxy radicals in the troposphere, J. Geophys. Res., 95D, 5697–5715, https://doi.org/10.1029/JD095ID05P05697, 1990. a

Mallik, C., Tomsche, L., Bourtsoukidis, E., Crowley, J. N., Derstroff, B., Fischer, H., Hafermann, S., Hüser, I., Javed, U., Keßel, S., Lelieveld, J., Martinez, M., Meusel, H., Novelli, A., Phillips, G. J., Pozzer, A., Reiffs, A., Sander, R., Taraborrelli, D., Sauvage, C., Schuladen, J., Su, H., Williams, J., and Harder, H.: Oxidation processes in the eastern Mediterranean atmosphere: evidence from the modelling of HO_x measurements over Cyprus, Atmos. Chem. Phys., 18, 10825–10847, https://doi.org/10.5194/acp-18-10825-2018, 2018. a

McKenna, D. S., Grooß, J.-U., Günther, G., Konopka, P., Müller, R., Carver, G., and Sasano, Y.: A new chemical Lagrangian model of the stratosphere (CLaMS) 2. Formulation of chemistry scheme and initialization, J. Geophys. Res., 107, 4256, https://doi.org/10.1029/2000JD000113, 2002. a

Meier, R. R., Anderson Jr., D. E., and Nicolet, M.: Radiation field in the troposphere and stratosphere from 240–1000 nm – I: General analysis, Planet. Space Sci., 30, 923–933, https://doi.org/10.1016/0032-0633(82)90134-9, 1982. a, b

Messaadia, L., Dib, G. E., Ferhati, A., and Chakir, A.: UV-visible spectra and gas-phase rate coefficients for the reaction of 2,3-pentanedione and 2,4-pentanedione with OH radicals, Chem. Phys. Lett., 626, 73–79, https://doi.org/10.1016/j.cplett.2015.02.032, 2015. a

Müller, J.-F., Peeters, J., and Stavrakou, T.: Fast photolysis of carbonyl nitrates from isoprene, Atmos. Chem. Phys., 14, 2497–2508, https://doi.org/10.5194/acp-14-2497-2014, 2014. a

Müller, R., Peter, T., Crutzen, P. J., Oelhaf, H., Adrian, G. P., von Clarmann, T., Wegner, A., Schmidt, U., and Lary, D.: Chlorine chemistry and the potential for ozone depletion in the Arctic stratosphere in the winter of 1991/92, Geophys. Res. Lett., 21, 1427–1430, https://doi.org/10.1029/94GL00465, 1994. a

Nakanishi, H., Morita, H., and Nagakura, S.: Electronic structures and spectra of the keto and enol forms of acetylacetone, B. Chem. Soc. Jpn., 50, 2255–2261, https://doi.org/10.1246/bcsj.50.2255, 1977. a

Nguyen, T. L., Peeters, J., and Vereecken, L.: Theoretical study of the gas-phase ozonolysis of β-pinene (C₁₀H₁₆), Phys. Chem. Chem. Phys., 11, 5643–5656, https://doi.org/10.1039/b822984h, 2009. a

Nicolet, M., Meier, R. R., and Anderson Jr., D. E.: Radiation field in the troposphere and stratosphere from 240–1000 nm – II. Numerical analysis, Planet. Space Sci., 30, 935–983, https://doi.org/10.1016/0032-0633(82)90135-0, 1982. a

Niemeyer, K. E. and Sung, C.-J.: On the importance of graph search algorithms for DRGEP-based mechanism reduction methods, Combust. Flame, 158, 1439–1443, https://doi.org/10.1016/j.combustflame.2010.12.010, 2011. a

Niemeyer, K. E., Sung, C.-J., and Raju, M. P.: Skeletal mechanism generation for surrogate fuels using directed relation graph with error propagation and sensitivity analysis, Combust. Flame, 157, 1760–1770, https://doi.org/10.1016/j.combustflame.2009.12.022, 2010. a, b

Nölscher, A., Butler, T., Auld, J., Veres, P., Muñoz, A., Taraborrelli, D., Vereecken, L., Lelieveld, J., and Williams, J.: Using total OH reactivity to assess isoprene photooxidation via measurement and model, Atmos. Environ., 89, 453–463, https://doi.org/10.1016/j.atmosenv.2014.02.024, 2014. a, b

Orlando, J. J. and Tyndall, G. S.: Laboratory studies of organic peroxy radical chemistry: an overview with emphasis on recent issues of atmospheric significance, Chem. Soc. Rev., 41, 6294–6317, https://doi.org/10.1039/C2CS35166H, 2012. a, b, c, d, e

Orlando, J. J., Tyndall, G. S., and Paulson, S. E.: Mechanism of the OH-initiated oxidation of methacrolein, Geophys. Res. Lett., 26, 2191–2194, https://doi.org/10.1029/1999GL900453, 1999. a

Paulot, F., Crounse, J. D., Kjaergaard, H. G., Kürten, A., St. Clair, J. M., Seinfeld, J. H., and Wennberg, P. O.: Unexpected epoxide formation in the gas-phase photooxidation of isoprene, Science, 325, 730–733, https://doi.org/10.1126/science.1172910, 2009. a, b, c

Peeters, J., Boullart, W., Pultau, V., Vandenberk, S., and Vereecken, L.: Structure-activity relationship for the addition of OH to (poly)alkenes: Site-specific and total rate constants, J. Phys. Chem. A, 111, 1618–1631, https://doi.org/10.1021/jp066973o, 2007. a, b, c, d

Peeters, J., Nguyen, T. L., and Vereecken, L.: HO_x radical regeneration in the oxidation of isoprene, Phys. Chem. Chem. Phys., 11, 5935–5939, https://doi.org/10.1039/B908511D, 2009. a, b

Peeters, J., Müller, J.-F., Stavrakou, T., and Nguyen, V. S.: Hydroxyl radical recycling in isoprene oxidation driven by hydrogen bonding and hydrogen tunneling: the upgraded LIM1 mechanism, J. Phys. Chem. A, 118, 8625–8643, https://doi.org/10.1021/jp5033146, 2014. a, b

Peixoto, T. P.: The graph-tool python library, Figshare, https://doi.org/10.6084/m9.figshare.1164194, 2014. a

Pepiot-Desjardins, P. and Pitsch, H.: An efficient error-propagation-based reduction method for large chemical kinetic mechanisms, Combust. Flame, 154, 67–81, https://doi.org/10.1016/j.combustflame.2007.10.020, 2008. a

Pommrich, R., Müller, R., Grooß, J.-U., Konopka, P., Ploeger, F., Vogel, B., Tao, M., Hoppe, C. M., Günther, G., Spelten, N., Hoffmann, L., Pumphrey, H.-C., Viciani, S., D'Amato, F., Volk, C. M., Hoor, P., Schlager, H., and Riese, M.: Tropical troposphere to stratosphere transport of carbon monoxide and long-lived trace species in the Chemical Lagrangian Model of the Stratosphere (CLaMS), Geosci. Model Dev., 7, 2895–2916, https://doi.org/10.5194/gmd-7-2895-2014, 2014. a

Praske, E., Crounse, J. D., Bates, K. H., Kurtén, T., Kjaergaard, H. G., and Wennberg, P. O.: Atmospheric fate of methyl vinyl ketone: peroxy radical reactions with NO and HO₂, J. Phys. Chem. A, 119, 4562–4572, https://doi.org/10.1021/jp5107058, 2015. a

Rickard, A. and Pascoe, S.: The Master Chemical Mechanism (MCM), available at: http://mcm.leeds.ac.uk (last access: 26 March 2019), 2009. a

Riede, H., Jöckel, P., and Sander, R.: Quantifying atmospheric transport, chemistry, and mixing using a new trajectory-box model and a global atmospheric-chemistry GCM, Geosci. Model Dev., 2, 267–280, https://doi.org/10.5194/gmd-2-267-2009, 2009. a

Sander, R., Kerkweg, A., Jöckel, P., and Lelieveld, J.: Technical note: The new comprehensive atmospheric chemistry module MECCA, Atmos. Chem. Phys., 5, 445–450, https://doi.org/10.5194/acp-5-445-2005, 2005. a

Sander, R., Baumgaertner, A., Gromov, S., Harder, H., Jöckel, P., Kerkweg, A., Kubistin, D., Regelin, E., Riede, H., Sandu, A., Taraborrelli, D., Tost, H., and Xie, Z.-Q.: The atmospheric chemistry box model CAABA/MECCA-3.0, Geosci. Model Dev., 4, 373–380, https://doi.org/10.5194/gmd-4-373-2011, 2011. a, b, c

Sander, R., Jöckel, P., Kirner, O., Kunert, A. T., Landgraf, J., and Pozzer, A.: The photolysis module JVAL-14, compatible with the MESSy standard, and the JVal PreProcessor (JVPP), Geosci. Model Dev., 7, 2653–2662, https://doi.org/10.5194/gmd-7-2653-2014, 2014. a

Sander, S. P., Friedl, R. R., Golden, D. M., Kurylo, M. J., Moortgat, G. K., Keller-Rudek, H., Wine, P. H., Ravishankara, A. R., Kolb, C. E., Molina, M. J., Finlayson-Pitts, B. J., Huie, R. E., and Orkin, V. L.: Chemical Kinetics and Photochemical Data for Use in Atmospheric Studies, Evaluation Number 15, JPL Publication 06-2, Jet Propulsion Laboratory, Pasadena, CA, available at: http://jpldataeval.jpl.nasa.gov (last access: 26 March 2019), 2006. a

Sander, S. P., Abbatt, J., Barker, J. R., Burkholder, J. B., Friedl, R. R., Golden, D. M., Huie, R. E., Kolb, C. E., Kurylo, M. J., Moortgat, G. K., Orkin, V. L., and Wine, P. H.: Chemical Kinetics and Photochemical Data for Use in Atmospheric Studies, Evaluation No. 17, JPL Publication 10-6, Jet Propulsion Laboratory, Pasadena, available at: http://jpldataeval.jpl.nasa.gov (last access: 26 March 2019), 2011. a, b

Sandu, A. and Sander, R.: Technical note: Simulating chemical systems in Fortran90 and Matlab with the Kinetic PreProcessor KPP-2.1, Atmos. Chem. Phys., 6, 187–195, https://doi.org/10.5194/acp-6-187-2006, 2006. a, b

Saunders, S. M., Jenkin, M. E., Derwent, R. G., and Pilling, M. J.: Protocol for the development of the Master Chemical Mechanism, MCM v3 (Part A): tropospheric degradation of non-aromatic volatile organic compounds, Atmos. Chem. Phys., 3, 161–180, https://doi.org/10.5194/acp-3-161-2003, 2003. a, b, c, d, e

Schultz, M. G., Stadtler, S., Schröder, S., Taraborrelli, D., Franco, B., Krefting, J., Henrot, A., Ferrachat, S., Lohmann, U., Neubauer, D., Siegenthaler-Le Drian, C., Wahl, S., Kokkola, H., Kühn, T., Rast, S., Schmidt, H., Stier, P., Kinnison, D., Tyndall, G. S., Orlando, J. J., and Wespes, C.: The chemistry–climate model ECHAM6.3-HAM2.3-MOZ1.0, Geosci. Model Dev., 11, 1695–1723, https://doi.org/10.5194/gmd-11-1695-2018, 2018. a, b

Sehested, J., Christensen, L. K., Nielsen, O. J., Bilde, M., Wallington, T. J., Schneider, W. F., Orlando, J. J., and Tyndall, G. S.: Atmospheric chemistry of acetone: Kinetic study of the ${CH}_{3} C (O) {CH}_{2} O_{2} + NO / {NO}_{2}$ reactions and decomposition of CH₃C(O)CH₂O₂NO₂, Int. J. Chem. Kinet., 30, 475–489, https://doi.org/10.1002/(SICI)1097-4601(1998)30:7<475::AID-KIN4>3.0.CO;2-P, 1998. a

Smith III, F. L. and Smith, C.: Numerical evaluation of Chapman's grazing incidence integral ch (X, χ), J. Geophys. Res., 77, 3592–3597, https://doi.org/10.1029/JA077i019p03592, 1972. a

So, S., Wille, U., and da Silva, G.: Atmospheric chemistry of enols: a theoretical study of the vinyl alcohol + OH + O₂ reaction mechanism, Environ. Sci. Technol., 48, 6694–6701, https://doi.org/10.1021/es500319q, 2014. a

Solomon, S. C., Hays, P. B., and Abreu, V. J.: The auroral 6300 Å emission – Observations and modeling, J. Geophys. Res., 93, 9867–9882, https://doi.org/10.1029/JA093iA09p09867, 1988. a

Strickland, D. J. and Meier, R. R.: A photoelectron model for the rapid computation of atmospheric excitation rates, Tech. Rep. ADA122871, Naval Research Lab Washington DC, available at: https://apps.dtic.mil/dtic/tr/fulltext/u2/a122871.pdf (last access: 26 March 2019), 1982. a

Taraborrelli, D., Lawrence, M. G., Butler, T. M., Sander, R., and Lelieveld, J.: Mainz Isoprene Mechanism 2 (MIM2): an isoprene oxidation mechanism for regional and global atmospheric modelling, Atmos. Chem. Phys., 9, 2751–2777, https://doi.org/10.5194/acp-9-2751-2009, 2009. a, b, c

Taraborrelli, D., Lawrence, M. G., Crowley, J. N., Dillon, T. J., Gromov, S., Groß, C. B. M., Vereecken, L., and Lelieveld, J.: Hydroxyl radical buffered by isoprene oxidation over tropical forests, Nat. Geosci., 5, 190–193, https://doi.org/10.1038/NGEO1405, 2012. a, b

Teng, A. P., Crounse, J. D., Lee, L., St. Clair, J. M., Cohen, R. C., and Wennberg, P. O.: Hydroxy nitrate production in the OH-initiated oxidation of alkenes, Atmos. Chem. Phys., 15, 4297–4316, https://doi.org/10.5194/acp-15-4297-2015, 2015. a

Tilmes, S., Lamarque, J.-F., Emmons, L. K., Kinnison, D. E., Marsh, D., Garcia, R. R., Smith, A. K., Neely, R. R., Conley, A., Vitt, F., Val Martin, M., Tanimoto, H., Simpson, I., Blake, D. R., and Blake, N.: Representation of the Community Earth System Model (CESM1) CAM4-chem within the Chemistry-Climate Model Initiative (CCMI), Geosci. Model Dev., 9, 1853–1890, https://doi.org/10.5194/gmd-9-1853-2016, 2016. a

Tobiska, W. K., Woods, T., Eparvier, F., Viereck, R., Floyd, L., Bouwer, D., Rottman, G., and White, O. R.: The SOLAR2000 empirical solar irradiance model and forecast tool, J. Atmos. Sol.-Terr. Phy., 62, 1233–1250, https://doi.org/10.1016/S1364-6826(00)00070-5, 2000. a

Tomlin, A. S. and Turányi, T.: Mechanism reduction to skeletal form and species lumping, in: Cleaner Combustion, edited by: Battin-Leclerc, F., Simmie, J. M., and Blurock, E., Springer Verlag, Berlin, 447–466, https://doi.org/10.1007/978-1-4471-5307-8_17, 2013. a

Tyndall, G. S., Cox, R. A., Granier, C., Lesclaux, R., Moortgat, G. K., Pilling, M. J., Ravishankara, A. R., and Wallington, T. J.: The atmospheric chemistry of small organic peroxy radicals, J. Geophys. Res., 106D, 12157–12182, https://doi.org/10.1029/2000JD900746, 2001. a

van Eijck, A., Opatz, T., Taraborrelli, D., Sander, R., and Hoffmann, T.: New tracer compounds for secondary organic aerosol formation from β-caryophyllene oxidation, Atmos. Environ., 80, 122–130, https://doi.org/10.1016/J.ATMOSENV.2013.07.060, 2013. a

Vereecken, L. and Peeters, J.: H-atom abstraction by OH-radicals from (biogenic) (poly)alkenes: C–H bond strengths and abstraction rates, Chem. Phys. Lett., 333, 162–168, https://doi.org/10.1016/S0009-2614(00)01347-6, 2001. a

Vereecken, L. and Peeters, J.: A theoretical study of the OH-initiated gas-phase oxidation mechanism of β-pinene (C₁₀H₁₆): first generation products, Phys. Chem. Chem. Phys., 14, 3802–3815, https://doi.org/10.1039/c2cp23711c, 2012. a

Vereecken, L., Müller, J.-F., and Peeters, J.: Low-volatility poly-oxygenates in the OH-initiated atmospheric oxidation of α-pinene: impact of non-traditional peroxyl radical chemistry, Phys. Chem. Chem. Phys., 9, 5241–5248, https://doi.org/10.1039/b708023a, 2007. a

Vereecken, L., Chakravarty, H. K., Bohn, B., and Lelieveld, J.: Theoretical study on the formation of H- and O-atoms, HONO, OH, NO, and NO₂ from the lowest lying singlet and triplet states in ortho-nitrophenol photolysis, Int. J. Chem. Kinet., 48, 785–795, https://doi.org/10.1002/kin.21033, 2016. a

Wallington, T. J., Ammann, M., Cox, R. A., Crowley, J. N., Herrmann, H., Jenkin, M. E., McNeill, V., Mellouki, A., Rossi, M. J., and Troe, J.: IUPAC Task group on atmospheric chemical kinetic data evaluation: Evaluated kinetic data, available at: http://iupac.pole-ether.fr (last access: 26 March 2019), 2018. a, b, c

Williams, J. E., van Velthoven, P. F. J., and Brenninkmeijer, C. A. M.: Quantifying the uncertainty in simulating global tropospheric composition due to the variability in global emission estimates of Biogenic Volatile Organic Compounds, Atmos. Chem. Phys., 13, 2857–2891, https://doi.org/10.5194/acp-13-2857-2013, 2013. a

Wolfe, G. M., Crounse, J. D., Parrish, J. D., Clair, J. M. S., Beaver, M. R., Paulot, F., Yoon, T. P., Wennberg, P. O., and Keutsch, F. N.: Photolysis, OH reactivity and ozone reactivity of a proxy for isoprene-derived hydroperoxyenals (HPALDs), Phys. Chem. Chem. Phys., 14, 7276–7286, https://doi.org/10.1039/c2cp40388a, 2012. a

Xiang, B., Zhu, L., and Tang, Y.: Photolysis of 4-oxo-2-pentenal in the 190–460 nm region, J. Phys. Chem. A, 111, 9025–9033, https://doi.org/10.1021/jp0739972, 2007. a

Yarwood, G., Rao, S., Yocke, M., and Whitten, G.: Updates to the carbon bond chemical mechanism: CB05, Final report RT-04-00675, United States Environmental Protection Agency, available at: http://www.camx.com/files/cb05_final_report_120805.aspx (last access: 26 March 2019), 2005. a

Articles

Short summary

We present the atmospheric chemistry box model CAABA/MECCA which now includes a number of new features: skeletal mechanism reduction, the MOM chemical mechanism for volatile organic compounds, an option to include reactions from the Master Chemical Mechanism (MCM) and other chemical mechanisms, updated isotope tagging, improved and new photolysis modules, and the new feature of coexisting multiple chemistry mechanisms. CAABA/MECCA is a community model published under the GPL.

The community atmospheric chemistry box model CAABA/MECCA-4.0

2.1 The Mainz Organic Mechanism (MOM)

2.1.1 VOC reactions with OH

2.1.2 RO2 reactions with NOx and NO3

2.1.3 RO2 reactions with HOx

2.1.4 RO2 permutation reactions

2.1.5 Photo-induced reactions

2.2 Other chemical mechanisms

2.2.1 CB05BASCOE

2.2.2 MOZART

2.2.3 JAM002

2.2.4 Master Chemical Mechanism (MCM)

2.3 Mechanism intercomparison

2.4 Skeletal mechanism reduction

2.5 Kinetic and isotope tagging

3.1 JVAL

3.2 RADJIMT

3.3 DISSOC

4.1 TRAJECT

4.2 PolyMECCA/CHEMGLUE

4.3 CHEMPROP

2.1.2 RO₂ reactions with NO_x and NO₃

2.1.3 RO₂ reactions with HO_x

2.1.4 RO₂ permutation reactions