Oxidation of low-molecular-weight organic compounds in cloud droplets: development of the Jülich Aqueous-phase Mechanism of Organic Chemistry (JAMOC) in CAABA/MECCA (version 4.5.0)
- 1Institute of Energy and Climate Research: Troposphere (IEK-8), Forschungszentrum Jülich GmbH, Jülich, Germany
- 2Atmospheric Chemistry Department, Max Planck Institute for Chemistry, Mainz, Germany
Correspondence: Simon Rosanka (firstname.lastname@example.org)
The Jülich Aqueous-phase Mechanism of Organic Chemistry (JAMOC) is developed and implemented in the Module Efficiently Calculating the Chemistry of the Atmosphere (MECCA; version 4.5.0). JAMOC is an explicit in-cloud oxidation scheme for oxygenated volatile organic compounds (OVOCs), suitable for global model applications. It is based on a subset of the comprehensive Cloud Explicit Physico-chemical Scheme (CLEPS; version 1.0). The phase transfer of species containing up to 10 carbon atoms is included, and a selection of species containing up to 4 carbon atoms reacts in the aqueous phase. In addition, the following main advances are implemented: (1) simulating hydration and dehydration explicitly; (2) taking oligomerisation of formaldehyde, glyoxal, and methylglyoxal into account; (3) adding further photolysis reactions; and (4) considering gas-phase oxidation of new outgassed species. The implementation of JAMOC in MECCA makes a detailed in-cloud OVOC oxidation model readily available for box as well as for regional and global simulations that are affordable with modern supercomputing facilities. The new mechanism is tested inside the box model Chemistry As A Boxmodel Application (CAABA), yielding reduced gas-phase concentrations of most oxidants and OVOCs except for the nitrogen oxides.
Aqueous-phase chemistry in cloud droplets differs significantly from gas-phase chemistry, mainly due to enhanced photolysis based on scattering effects within cloud droplets (Bott and Zdunkowski, 1987; Mayer and Madronich, 2004), faster reaction rates, and ion reactions that do not occur in the gas phase (Herrmann, 2003; Epstein and Nizkorodov, 2012). Moreover, conversion of nitrogen monoxide (NO) to nitrogen dioxide (NO2) by peroxy radicals (RO2) essentially does not take place in aqueous droplets because NO is insoluble. In the aqueous phase, oxygenated volatile organic compounds (OVOCs) are mainly oxidised during the daytime by the hydroxyl radical (OH) and by the nitrate radical (NO3) during the nighttime (Herrmann et al., 2015). Even though ozone (O3) is not very soluble, it can be taken up into cloud droplets where it is destroyed by
for which the superoxide anion () is in equilibrium with its conjugate acid, the hydroperoxyl radical (HO2). This indicates that the in-cloud O3 destruction is sensitive to in-cloud OVOC oxidation. Lelieveld and Crutzen (1990) have already proposed that clouds can influence HOx () and NOx (), resulting in regional changes of up to 40 % in particular locations, being subject to cloud processing. At the tropics and mid-latitudes, Liang and Jacob (1997) suggest that clouds may reduce O3 by 3 % in summer. By changing the gas-phase oxidant budgets, clouds can indirectly influence the formation of secondary organic aerosols (SOAs). Within cloud droplets, OVOC oxidation additionally can lead to the formation and destruction of SOA precursors, and clouds can act as SOA sources (Blando and Turpin, 2000). Further modelling studies suggest that clouds may contribute to the SOA formation on a par with gas-phase sources (Ervens et al., 2011; Lin et al., 2012; Ervens, 2015). By scattering, SOAs are known to influence the aerosol optical depth (AOD), leading to a reduction in NO2 photolysis (Tie et al., 2005). In addition, SOAs may act as cloud condensation nuclei (CCN) (Andreae and Rosenfeld, 2008), affecting cloud properties. An increased formation of SOAs would thus influence tropospheric HOx and O3 chemistry.
When performing global modelling studies, it is thus desirable to include the in-cloud oxidation of OVOCs. However, compared to gas-phase chemistry, knowledge of aqueous-phase chemistry still suffers from large uncertainties and most global models only include very limited representations. Most global models include only the uptake of a few soluble compounds, their acid–base equilibria, and the oxidation of sulfur dioxide (SO2) by ozone (O3) and hydrogen peroxide (H2O2) (Ervens, 2015, their Table 1). The explicit oxidation of OVOCs is currently not considered in any global model, with one exception – though limited to species containing one carbon atom (Tost et al., 2006). Mouchel-Vallon et al. (2017) recently presented the Cloud Explicit Physico-chemical Scheme (CLEPS; version 1.0), a complex new oxidation scheme coupled to the gas-phase Master Chemical Mechanism (MCM; version 3.3.1; Jenkin et al., 2015). However, their comprehensive mechanism is targeted for box-model applications and is not suitable for global model applications due to its complexity.
In this study, the in-cloud OVOC oxidation scheme Jülich Aqueous-phase Mechanism of Organic Chemistry (JAMOC) is presented and implemented into the chemistry mechanism Module Efficiently Calculating the Chemistry of the Atmosphere (MECCA). Here, JAMOC's representation of organic chemistry is based on CLEPS and is thus an addition to MECCA's existing aqueous-phase chemical mechanism. Therefore, JAMOC needs to be selected by the user upon compilation of MECCA's chemical mechanism. A visualisation of this procedure can be found in MECCA's user manual available in the archived model code (caaba_mecca_manual.pdf). The modular structure of MECCA allows it to be connected to different base models, e.g. to the Chemistry As A Boxmodel Application (CAABA) by Sander et al. (2019) or to the global ECHAM/MESSy Atmospheric Chemistry Model (EMAC) by Jöckel et al. (2010). In this combination, the proposed mechanism closes the gap between box models and global model applications. In addition to the new aqueous-phase OVOC chemistry, MECCA also contains the gas-phase Mainz Organic Mechanism (MOM; Sander et al., 2019) with an extensive oxidation scheme for isoprene (Taraborrelli et al., 2009, 2012; Nölscher et al., 2014), monoterpenes (Hens et al., 2014), and aromatics (Cabrera-Perez et al., 2016). VOCs are oxidised by OH, O3, and NO3, whereas RO2 reacts with HO2, NOx, and NO3 and undergoes self- and cross-reactions (Sander et al., 2019).
The mechanism of JAMOC is described in Sect. 2, followed by a short description of its implications in the box model CAABA (Sect. 3). Global implications are analysed in our companion paper (Rosanka et al., 2021a), and the mechanism's importance for global models simulating extreme pollution events is addressed by Rosanka et al. (2020). Modelling uncertainties are discussed in Sect. 4 before drawing final conclusions in Sect. 5.
The detailed mechanism CLEPS includes 850 aqueous-phase reactions, focusing on the oxidation of species containing up to four carbon atoms (Mouchel-Vallon et al., 2017). Since our target is to simulate OVOC chemistry inside the global model EMAC (Rosanka et al., 2021a), using such a large mechanism is not feasible. Therefore, we have developed the reduced mechanism JAMOC. Only a selection of species containing up to four carbon atoms is considered in the aqueous phase. The gas-phase oxidation of the most abundant hydrocarbons (e.g. methane, isoprene) leads to many highly soluble organic species with one or two carbon atoms (e.g. formaldehyde, methanol, glyoxal). In order to properly represent these degradation products, JAMOC includes the aqueous-phase oxidation of all species containing one and two carbon atoms treated in CLEPS. Even though isoprene (C5H8), the biogenic VOC emitted the most (Guenther et al., 2012), is not soluble, the representation of its oxidation products containing more than two carbon atoms (i.e. methylglyoxal, methacrolein, methyl vinyl ketone) is desirable for global model applications. Therefore, the oxidation of species containing three carbon atoms in JAMOC focuses on the representation of the aqueous-phase oxidation of methylglyoxal and its aqueous-phase oxidation products (e.g. pyruvic acid). Additionally, its aqueous-phase sources from acetone, hydroxy acetone, isopropanol, hydroperoxide, and isopropyl hydroperoxide are included. The aqueous-phase oxidation of species containing four carbon atoms in JAMOC is limited to methacrolein (MACR) and methyl vinyl ketone (MVK). The phase transfer of species containing up to 10 carbon atoms is included so that their wet deposition can be represented in global model applications (i.e. by using EMAC; see Rosanka et al., 2021a). In order to reduce the stiffness of the ordinary differential equation (ODE) system and the required computational demand, the representation of organic radicals (see Sect. 2.7) is simplified. It is assumed that the following reactions occur instantly and are not explicitly represented in the ODE system if they are the only fate of the respective radical: (1) the O2 addition to alkyl radicals; (2) the HO2 elimination of α-hydroxyperoxyl; and (3) the carbon bond scission or 1,2-hydrogen shift of alkoxyl radicals. In addition to the chemistry from CLEPS, JAMOC includes (1) explicit hydration and dehydration; (2) oligomerisation of formaldehyde, glyoxal, and methylglyoxal as an in-cloud SOA source; (3) further aqueous-phase photolysis reactions; and (4) the gas-phase photo-oxidation of new outgassed species. The complete aqueous-phase mechanism represents the phase transfer of 368 species, 68 equilibria (acid–base and hydration–dehydration), 402 reactions, and 27 aqueous-phase photolysis reactions. In the gas phase, 1 photolysis and 18 OH oxidation reactions are added to MOM. A list detailing the complete mechanism is available in the archived model code.
This section provides a general overview of the developed mechanism. For completeness, short summaries of CLEPS are provided if no significant difference exists between both mechanisms. Figures 1 and 2 give a graphical representation of all parts of the developed mechanism, using glyoxal and oxalic acid as examples.
2.1 Inorganic chemistry
The inorganic chemistry for the proposed mechanism is very similar to the inorganic chemistry of the standard aqueous-phase mechanism used in EMAC (Tost et al., 2007; Jöckel et al., 2016). In this standard mechanism, the major aqueous-phase O3 sink, the reaction with , is represented as
In JAMOC, this aqueous-phase O3 chemistry is updated to the mechanism proposed by Staehelin et al. (1984) with corrections from Staehelin and Hoigné (1985), in which the O3 destruction by is represented as given in Reaction (R1).
2.2 Uptake of gaseous species into cloud droplets
The mass transfer of species between the gas and the aqueous phase is described following Schwartz (1986) (see Sander, 1999; Tost et al., 2006). The explicit bidirectional phase transfer of 45 carbon-containing species, which explicitly react in the aqueous phase, is considered (indicated in pink in Figs. 1 and 2). In this model framework, Henry's law constants are mainly taken from Sander (2015), Burkholder et al. (2015), and sources therein. In order to account for the hydration of aldehydes (for more details see Sect. 2.3), a distinction is made between the effective Henry's law constant (H*) and the intrinsic Henry's law constant (H). The latter is calculated by
where Khyd is the ratio between the forward and reverse kinetic rate constant of the hydration equilibrium (see Reaction R3). Table 1 gives an overview of the hydration constants and the effective Henry's law constants, including the resulting intrinsic Henry's law constants, for all aldehydes. The temperature dependencies of the intrinsic Henry's law constants are assumed to be the same as for the effective constants. The accommodation constant (α) is known for a few species; if unknown, the standard EMAC estimate of 0.1 is used. In addition to the phase transfer of all species that explicitly react in the aqueous phase, the phase transfer of all soluble MOM species containing up to 10 carbon atoms is represented in order to allow their removal by wet deposition in global models (i.e. by using EMAC; see Rosanka et al., 2021a). A list summarising all Henry's law and accommodation constants is available in the archived model code.
2.3 Hydration of carbonyls
Gem-diols are formed when aldehydes (carbonyl compounds) hydrate:
In the new mechanism, 12 carbonyl species undergo hydration (indicated with blue arrows in Fig. 1). The monohydrate of glyoxal (dihydroxyacetaldehyde) undergoes additional hydration to form its dihydrate (1,1,2,2-ethanetetrol). Pseudo-first-order rate constants for the hydration and dehydration are mainly obtained from the literature (e.g. Doussin and Monod, 2013). In the case of formyldioxidanyl and hydroperoxyacetaldehyde, the pseudo-first-order rate constants are assumed to be the same as for formaldehyde and glycolaldehyde, respectively.
The typical lifetime of a warm cloud droplet can be several minutes, but their typical evaporation timescale is less than 100 s (Jarecka et al., 2013). Following the dehydration constants presented by Doussin and Monod (2013), the dehydration of some gem-diols can be slower than the typical cloud droplet evaporation timescale. Additionally, their rapid transfer across the phases is expected to affect the gas-phase concentration of gem-diols, for which no other significant source is known. This process could be an important removal of gem-diols from the aqueous phase, without yielding the original aldehyde. Therefore, their outgassing is considered for use with the models representing evaporating clouds like the EMAC model (following Sect. 2.2). However, their Henry's law constants are unknown. Thus, estimates are obtained at 25 ∘C using the bond method (Meylan and Howard, 1991) from the United States Environmental Protection Agency Estimation Programs Interface (EPI) Suite (United States Environmental Protection Agency (US EPA), 2012). An overview of all estimated effective Henry's law constants is given in Table 2.
In CLEPS, acyl peroxy radicals (RC(O)(OO)) are assumed to be in a hydration–dehydration equilibrium similarly to their parent aldehydes (Mouchel-Vallon et al., 2017). However, experimental results by Villalta et al. (1996) show that in the case of peroxyacetyl radicals (CH3C(O)(OO)), no equilibrium exists. Instead, hydrolysis takes place, likely yielding acetic acid (CH3CO2H) and HO2. It is thus assumed that all acyl peroxy radicals undergo hydrolysis following Reaction (R4), with a reaction rate constant of 7.0×105 M−1 s−1, as proposed by Villalta et al. (1996).
2.4 Acid dissociation
The dissociation of acids is taken into account following
which is indicated in green in Fig. 2. The acidity constants (Ka) for most of the one-carbon-, two-carbon-, and three-carbon-containing acids taken into account in JAMOC are known from the literature (Rumble, 2020). If unknown, the acidity constants are used as proposed by Mouchel-Vallon et al. (2017). The dissociation and association rate constants are selected such that the equilibrium between dissociation and association is reached quickly, while still avoiding numerical stiffness problems in the numerical integrator.
2.5 Oxidation by OH, NO3, and other oxidants
In JAMOC, OH and NO3 are the main oxidants taken into account. Reactions of OVOCs with oxidants are treated as proposed by Mouchel-Vallon et al. (2017). Organic compounds may react in three different ways with OH radicals (Herrmann et al., 2015), each indicated in orange in Figs. 1 and 2. They form an alkyl radical following H abstraction:
If the organic compound contains a double bond, OH addition is favoured.
With anions like carboxylates, electron transfer takes place.
When available, rate constants are obtained from the literature. If unavailable, the rate constant for the H abstraction is estimated based on the structure–activity relationship (SAR) by Doussin and Monod (2013), which for carboxylate compounds is extended to account for the electron transfer as described by Mouchel-Vallon et al. (2017). In all cases, branching ratios are obtained from the SAR with simplifications by Mouchel-Vallon et al. (2017).
During the nighttime, OH radical concentrations are low and, due to missing photolysis, NO3 radicals are considered the main nighttime oxidant. Similar to CLEPS, JAMOC only considers the H abstraction leading to alkyl radicals for NO3 reactions (Herrmann et al., 2015):
For most species containing one or two carbon atoms, rate constants are obtained from the literature. In contrast to OH, no SAR is available for the H abstraction by NO3. Therefore, rate constants are obtained from the similar criteria described by Mouchel-Vallon et al. (2017). Due to missing branching ratios from the literature, branching ratios are assumed to be the same as for the H abstraction by OH.
In addition to reactions of organic compounds with OH and NO3, reactions with other oxidants are implemented when available from the literature. The oxidants considered here are , O3, H2O2, , and sulfur-containing oxidants ( and ). For all oxidation reactions, reaction rates and branching ratios are either taken from the literature or as proposed by Mouchel-Vallon et al. (2017).
The formation of oligomers within the atmospheric aqueous phase is known to be a source of SOAs. Even though Tan et al. (2009) suggest that the formation of oligomers becomes increasingly important for aerosol water, where precursor concentrations are found to be higher, Lin et al. (2012) have demonstrated that SOA formation from cloud processing is globally important. Therefore, JAMOC includes self- and cross-reactions leading to oligomers for formaldehyde, glyoxal, and methylglyoxal. The oligomerisation of formaldehyde is implemented following Hahnenstein et al. (1995), in which the methanediol formed from hydrolysis (see Sect. 2.3) reacts with itself and the dimer formed from this self-reaction. Ervens and Volkamer (2010) studied the oligomerisation of glyoxal. Here, glyoxal and its hydrates react with the monohydrate to form three oligomers (indicated in green in Fig. 1). The oligomerisation of methylglyoxal is assumed to follow the same mechanisms as for glyoxal. However, only the monohydrate of methylglyoxal is taken into account in this mechanism, leading to only two oligomers. Each oligomer is assumed to react with OH, leading to HO2, with reaction rate constants that are double for the corresponding (hydrated) monomer due to an increased number of abstractable H atoms.
2.7 Organic radicals
Organic radicals are generally treated following Mouchel-Vallon et al. (2017). Alkyl radicals can either form oligomers via self- and cross-reactions (e.g. Lim et al., 2013; Ervens et al., 2015) or undergo O2 addition:
As proposed by Mouchel-Vallon et al. (2017), it is assumed that O2 addition is the fastest pathway, due to high O2 concentrations following fast O2 saturation in cloud droplets (Ervens, 2015). Thus, oligomers formed from the self- and cross-reactions of alkyl radicals are not considered in JAMOC.
Peroxyl radicals generally undergo self- or cross-reactions forming short-lived tetroxides that quickly decompose (von Sonntag and Schuchmann, 1997). Due to limited computation resources, only self-reactions are taken into account. Mouchel-Vallon et al. (2017) propose three similarity criteria for the decomposition of tetroxides depending on the peroxyl radical: (1) for β-peroxycarboxylic acids (RC(OO)C(=O)(OH)) experimental results from Schuchmann et al. (1985) are generalised, (2) β-hydroxyperoxyl radicals () are represented according to Piesiak et al. (1984), and (3) β-oxoperoxyl radicals () are treated based on Zegota et al. (1986) and Poulain et al. (2010). If some products are unknown, branching ratios of the known products are rescaled to 100 % in order to preserve mass. The peroxyl radicals undergo HO2 elimination (von Sonntag, 1987) if the hydroxyl moiety is in the alpha position (α-hydroxyperoxyl).
The generalised corresponding rate constants are used as proposed by Mouchel-Vallon et al. (2017, Table 3), which are based on the work of von Sonntag (1987). In CLEPS, peroxyl radicals additionally undergo elimination when reacting with OH− (Zegota et al., 1986; Mouchel-Vallon et al., 2017).
In order to decrease the number of reactions and due to the fast HO2 elimination, this elimination is not considered explicitly in JAMOC.
Acyl peroxy radicals (RC(O)(OO)) are treated like peroxyl radicals, as described in Monod et al. (2007), but only form alkoxyl radicals. Peroxyl radicals that have not explicitly been discussed so far are treated following Monod et al. (2007) (Mouchel-Vallon et al., 2017).
Mouchel-Vallon et al. (2017) suggest that alkoxyl radicals (RO) undergo either a carbon bond scission (Hilborn and Pincock, 1991) if the neighbouring carbon atom is oxygenated (Reaction R13) or a 1,2-hydrogen shift (DeCosta and Pincock, 1989) if the neighbouring carbon atom is not oxygenated (Reaction R14).
In general, the photolysis of some organic compounds (e.g. organic peroxides, pyruvic acid) competes with other oxidation pathways (see Sect. 2.5) and can be a source of OH. In Rosanka et al. (2021a), a global tropospheric in-cloud OH budget is presented. When using JAMOC, EMAC predicts that about 40 % of all in-cloud OH is produced from the photolysis of a selection of organic compounds. However, Fenton chemistry is not considered by Rosanka et al. (2021a), and the relative contribution is therefore expected to be overestimated. The photolysis of glyoxal and oxalic acid is indicated in orange in Figs. 1 and 2. The number of photolytic reactions known from the literature, of which some are implemented in CLEPS (Mouchel-Vallon et al., 2017), is limited. In JAMOC, the photolysis of additional compounds is taken into account. This includes the photolysis of oxalic acid ((COOH)2), which is implemented following Yamamoto and Back (1985) using the ultraviolet absorption spectrum presented in Back (1984). If available, additional photolysis reactions are implemented following Sander et al. (2014). In order to account for scattering effects within cloud droplets (Ruggaber et al., 1997), an enhancement factor of 2.33, the same as that used in EMAC's standard aqueous-phase mechanism for the photolysis of H2O2 (Tost et al., 2007; Jöckel et al., 2016), is applied to each gas-phase photolysis rate.
2.9 Gas-phase oxidation of new species
Oxalic acid was not represented in the gas-phase mechanism (i.e. in MOM). The gas-phase oxidation of oxalic acid by OH and its photolysis are implemented in order to realistically represent oxalic acid in the gas phase. Similarly to the implementation in the aqueous phase, the photolysis of oxalic acid is implemented following Yamamoto and Back (1985) and Back (1984). All gem-diols (see Sect. 2.3) formed from hydration are transferred to the gas phase and oxidised by OH (indicated in orange in Figs. 1 and 2). All OH oxidation reaction rates are estimated following the description of Sander et al. (2019).
The implications of the developed mechanism are tested by comparing it to the minimum in-cloud oxidation scheme available in CAABA/MECCA and EMAC. The minimum mechanism only includes the uptake of a few soluble compounds, their acid–base equilibria, and the oxidation of SO2 by O3 and H2O2 (Jöckel et al., 2006). This minimal mechanism is thus representative of most global models (Ervens, 2015). For both mechanisms, an air parcel is simulated in CAABA, taking the same conditions into account: the air parcel is simulated during summer at a mid-latitude with a constant temperature of 278 K and relative humidity of 100 %. Table 3 provides a selection of initial mixing ratios and emission fluxes of gas-phase species treated in MOM. The initial conditions are a modified version of the scenario used by Taraborrelli et al. (2009). Within the air parcel, a stable cloud droplet population is simulated with a radius of 20 µm and a liquid water content of 0.3 g m−3. Both simulations are intended as a sensitivity study of JAMOC. Therefore, CAABA is initialised at 00:00 UTC and simulates the air parcel for 5 d in total. A realistic cloud event with a cloud droplet lifetime of 1 h using CAABA is presented in Rosanka et al. (2021a). In addition, Rosanka et al. (2021a) study the implications of JAMOC on a global scale using EMAC.
Figure 3 gives an overview of the temporal development of the total mixing ratios (gas + aqueous phase) for a selection of species during the simulated daily cycles of 5 d. Comparing the new and the minimum mechanisms, it becomes clear that the newly developed mechanism has a significant impact on most trace gases. With the explicit oxidation of many OVOCs in the aqueous phase, the mixing ratio of the sum of all OVOCs explicitly reacting in JAMOC (∑OVOCs; see Eq. A1 in Appendix A) is significantly reduced. This reduction is a combined effect from (1) the in-cloud oxidation of these OVOCs and (2) their dampened gas-phase production. In the gas-phase, most OVOCs are formed by secondary production (e.g. oxidation of primarily emitted VOCs). The decrease in the main VOC oxidant (i.e. OH) leads to reduced oxidation of primarily emitted VOCs resulting in a reduced gas-phase OVOC formation. The calculated diurnal cycles of OH, HO2, NOx, and O3 are similar for both mechanisms and differ mainly in the absolute mixing ratios calculated. When JAMOC is used, HO2 partitions into the cloud droplets, whereas NO stays in the gas phase due to its low solubility (Jacob, 1986; Lelieveld and Crutzen, 1990). This results in substantial changes in the NOx–HOx relation, resulting in reduced OH formation from its second-most-important atmospheric gas-phase source:
Overall, this results in reduced HOx and elevated NOx mixing ratios. In addition, lower HO2 mixing ratios lead to a reduced removal of NOx by the formation of nitric acid (HNO3) and peroxynitric acid (HNO4). Within the cloud droplet, is in equilibrium with its conjugated base HO2. Higher in-cloud HO2 concentrations, caused by mass transfer and in-cloud OVOC oxidation, consequently lead to an increased destruction of O3 via Reaction (R1). This results in an enhanced uptake of O3 into the cloud droplet and an increased importance of cloud droplets as O3 sinks.
The impact of the newly proposed mechanism is consistent with earlier box-model studies. The reduction in OVOCs is similar to the findings given in Mouchel-Vallon et al. (2017) when using CLEPS. In contrast, the reduction in methylglyoxal differs since, in CLEPS, gas-phase methylglyoxal mixing ratios first increase and later decrease during the modelled cloud event of Mouchel-Vallon et al. (2017). This difference is most likely linked to the usage of the intrinsic Henry's law constant and the explicit representation of the methylglyoxal hydration–dehydration in JAMOC. In contrast to Mouchel-Vallon et al. (2017), CAABA predicts a reduction in OH levels. However, this reduction in OH is in line with other modelling studies predicting a similar reduction in gas-phase OH during cloud events (Tilgner et al., 2013). It is important to keep in mind that in Mouchel-Vallon et al. (2017), a different cloud event is simulated, including different initial conditions and a different emission scenario. In their study, the cloud forms after a certain time period, whereas in CAABA the cloud is present the whole time.
The uncertainties associated with the present kinetic model are mainly attributed to (1) assumptions and simplifications in the aqueous-phase mechanism and (2) missing sinks of key oxidants. Each possible uncertainty is discussed in this section.
In general, aqueous-phase kinetics data suffer from many large uncertainties compared to the data available for the gas phase. In the development of the implemented in-cloud oxidation scheme JAMOC, some assumptions are made that introduce modelling uncertainties. If rate constants are unknown, estimates are taken from Mouchel-Vallon et al. (2017). These are based on a structure–activity relationship (SAR) for the H abstraction by OH for dissolved carbonyls and carboxylic acids considered in this study (Doussin and Monod, 2013). However, it is expected that the uncertainty in the estimated rate constants is low since Doussin and Monod (2013) report that when evaluated using experimental data, their estimates were within ±20 % for 58 % of the calculated rate constants. Also the up-scaling of branching ratios to conserve mass further influences the predictions of VOC oxidation. The mechanism should be updated with rate constants and branching ratios as soon as experimental results become available. The increased concentration and burden of certain organic acids heavily depend on the chemistry and solubility of some gem-diols. For example, the gas-phase oxidation of the methylglyoxal monohydrate leads to the formation of pyruvic acid. The gas-phase production of pyruvic acid therefore depends on the mass transfer of this specific monohydrate. In the current implementation, the Henry's law constants for all gem-diols are estimated. For the methylglyoxal monohydrate, the estimated values range from 3.5×103 to 2.4×104 M atm−1.
Phase transfer of soluble VOCs into cloud droplets is considered in JAMOC even when their oxidation is not explicitly represented (see Sect. 2.2). This allows their removal from the atmosphere by rain-out when JAMOC is connected to a global model (e.g. using EMAC; see Rosanka et al., 2021a). Arakaki et al. (2013) point out that by not taking the oxidation of all dissolved organic carbon (DOC) into account, aqueous-phase OH concentrations might be overestimated. Based on observational estimates, they suggest a general scavenging rate constant of M−1 s−1 for all DOC. If each DOC species reacts with OH, the gas-phase concentration would be reduced, further influencing gas-phase VOC concentrations and the overall oxidation capacity. Implementing the DOC oxidation, suggested by Arakaki et al. (2013), for every scavenged DOC species would increase the aqueous-phase mechanism by more than 280 reactions, which is almost a doubling of the proposed organic mechanism. Within the scope of this study, it is thus computationally not feasible to include this additional OH sink. Currently, the model runtime increases from 4.3 s for EMACs minimum in-cloud oxidation scheme to 6.5 s for the newly proposed mechanism JAMOC.
Reducing the model uncertainties introduced by estimates of Henry's law constants of gem-diols and missing in-cloud DOC oxidation is outside the scope of this study due to the uncertainties' complexity. Model representation of the latter is expected to influence the oxidation rate of VOCs in the cloud droplets and aerosols.
In this study, the new in-cloud oxidation scheme of soluble VOCs JAMOC is developed and implemented into MECCA. This mechanism is suitable for global model applications and based on the box-model mechanism CLEPS proposed by Mouchel-Vallon et al. (2017). The mechanism considers the phase transfer of OVOCs containing up to 10 carbon atoms. For a selection of OVOCs containing up to 4 carbon atoms, their acid–base and/or hydration–dehydration equilibria and their reactions with OH, NO3, and other oxidants (if available) are explicitly represented. Additionally, the gas-phase photo-oxidation of gem-diols and oxalic acid is implemented into the gas-phase mechanism MOM. Finally, JAMOC is tested within the CAABA box model.
The proposed mechanism leads to a significant reduction in OVOCs and an overall reduction in important oxidants. These findings are in line with other box-model studies and demonstrate the importance of in-cloud chemistry in atmospheric chemistry. By not taking the in-cloud oxidation of OVOCs into account, global models will tend to overestimate the levels of OVOCs and atmospheric oxidants. A complete analysis on the importance of JAMOC at a global scale is presented in Rosanka et al. (2021a). In future studies, the modular implementation of JAMOC, with the necessary adjustments, will allow its application to aerosol water.
In Fig. 3, the mixing ratios of the sum of all the OVOCs explicitly reacting in JAMOC (∑OVOCs) are shown. In these cases, ∑ OVOCs is defined as follows:
The current version of the CAABA/MECCA model code is available as a community model in the code repository at https://gitlab.com/RolfSander/caaba-mecca (last access: 25 May 2021, Sander, 2021a), published under the GNU General Public License (http://www.gnu.org/copyleft/gpl.html, last access: 23 April 2021).
The exact version of the CAABA/MECCA model (version 4.5.0) developed in this paper and used in each simulation presented in this paper is archived at Zenodo (http://doi.org/10.5281/zenodo.4707938; Sander, 2021b). All future versions of CAABA/MECCA will be made available at https://doi.org/10.5281/zenodo.4707937.
The archived model code includes a list of all chemical reactions including rate constants and references (caaba/manual/meccanism.pdf), a list of all Henry's law and accommodation constants (caaba/tools/chemprop/chemprop.pdf), and a user manual (caaba/manual/caaba_manual_manual.pdf). For further information and updates, the CAABA/MECCA web page at http://www.mecca.messy-interface.org (last access: 23 April 2021) can be consulted.
SR and DT developed the chemical mechanism. The chemical mechanism was reviewed by RS. SR, DT, and RS implemented the mechanism into MECCA. The results were discussed by all co-authors. The manuscript was prepared by SR with the help of all co-authors.
The authors declare that they have no competing interests.
The work described in this paper has received funding from the Initiative and Networking Fund of the Helmholtz Association through the project Advanced Earth System Modelling Capacity (ESM). The content of this paper is the sole responsibility of the authors, and it does not represent the opinion of the Helmholtz Association, and the Helmholtz Association is not responsible for any use that might be made of the information contained. The authors gratefully acknowledge the Earth System Modelling (ESM) project for funding this work by providing computing time on the ESM partition of the supercomputer JUWELS at the Jülich Supercomputing Centre (JSC).
This research has been supported by the Initiative and Networking Fund of the Helmholtz Association through the project Advanced Earth System Modelling Capacity (ESM) (grant no. DB001549).
The article processing charges for this open-access publication were covered by the Forschungszentrum Jülich.
This paper was edited by Christoph Knote and reviewed by two anonymous referees.
Andreae, M. and Rosenfeld, D.: Aerosol–cloud–precipitation interactions. Part 1. The nature and sources of cloud-active aerosols, Earth-Sci. Rev., 89, 13–41, https://doi.org/10.1016/j.earscirev.2008.03.001, 2008. a
Arakaki, T., Anastasio, C., Kuroki, Y., Nakajima, H., Okada, K., Kotani, Y., Handa, D., Azechi, S., Kimura, T., Tsuhako, A., and Miyagi, Y.: A General Scavenging Rate Constant for Reaction of Hydroxyl Radical with Organic Carbon in Atmospheric Waters, Environ. Sci. Technol., 47, 8196–8203, https://doi.org/10.1021/es401927b, 2013. a, b
Blando, J. D. and Turpin, B. J.: Secondary organic aerosol formation in cloud and fog droplets: a literature evaluation of plausibility, Atmos. Environ., 34, 1623–1632, https://doi.org/10.1016/S1352-2310(99)00392-1, 2000. 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: https://jpldataeval.jpl.nasa.gov/pdf/JPL_Publication_15-10.pdf (last access: 26 May 2021), 2015. a, b
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
DeCosta, D. P. and Pincock, J. A.: Control of product distribution by Marcus type electron-transfer rates for the radical pair generated in benzylic ester photochemistry, J. Am. Chem. Soc., 111, 8948–8950, https://doi.org/10.1021/ja00206a045, 1989. a
Doussin, J.-F. and Monod, A.: Structure–activity relationship for the estimation of OH-oxidation rate constants of carbonyl compounds in the aqueous phase, Atmos. Chem. Phys., 13, 11625–11641, https://doi.org/10.5194/acp-13-11625-2013, 2013. a, b, c, d, e, f
Epstein, S. A. and Nizkorodov, S. A.: A comparison of the chemical sinks of atmospheric organics in the gas and aqueous phase, Atmos. Chem. Phys., 12, 8205–8222, https://doi.org/10.5194/acp-12-8205-2012, 2012. a
Ervens, B. and Volkamer, R.: Glyoxal processing by aerosol multiphase chemistry: towards a kinetic modeling framework of secondary organic aerosol formation in aqueous particles, Atmos. Chem. Phys., 10, 8219–8244, https://doi.org/10.5194/acp-10-8219-2010, 2010. a, b, c
Ervens, B., Turpin, B. J., and Weber, R. J.: Secondary organic aerosol formation in cloud droplets and aqueous particles (aqSOA): a review of laboratory, field and model studies, Atmos. Chem. Phys., 11, 11069–11102, https://doi.org/10.5194/acp-11-11069-2011, 2011. a
Ervens, B., Renard, P., Tlili, S., Ravier, S., Clément, J.-L., and Monod, A.: Aqueous-phase oligomerization of methyl vinyl ketone through photooxidation – Part 2: Development of the chemical mechanism and atmospheric implications, Atmos. Chem. Phys., 15, 9109–9127, https://doi.org/10.5194/acp-15-9109-2015, 2015. a
Guenther, A. B., Jiang, X., Heald, C. L., Sakulyanontvittaya, T., Duhl, T., Emmons, L. K., and Wang, X.: The Model of Emissions of Gases and Aerosols from Nature version 2.1 (MEGAN2.1): an extended and updated framework for modeling biogenic emissions, Geosci. Model Dev., 5, 1471–1492, https://doi.org/10.5194/gmd-5-1471-2012, 2012. a
Hahnenstein, I., Albert, M., Hasse, H., Kreiter, C. G., and Maurer, G.: NMR Spectroscopic and Densimetric Study of Reaction Kinetics of Formaldehyde Polymer Formation in Water, Deuterium Oxide, and Methanol, Ind. Eng. Chem. Res., 34, 440–450, https://doi.org/10.1021/ie00041a003, 1995. 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 HOx radicals in a boreal forest, Atmos. Chem. Phys., 14, 8723–8747, https://doi.org/10.5194/acp-14-8723-2014, 2014. a
Herrmann, H., Schaefer, T., Tilgner, A., Styler, S. A., Weller, C., Teich, M., and Otto, T.: Tropospheric Aqueous-Phase Chemistry: Kinetics, Mechanisms, and Its Coupling to a Changing Gas Phase, Chem. Rev., 115, 4259–4334, https://doi.org/10.1021/cr500447k, 2015. a, b, c
Hilborn, J. W. and Pincock, J. A.: Rates of decarboxylation of acyloxy radicals formed in the photocleavage of substituted 1-naphthylmethyl alkanoates, J. Am. Chem. Soc., 113, 2683–2686, https://doi.org/10.1021/ja00007a049, 1991. a
Jacob, D. J.: Chemistry of OH in remote clouds and its role in the production of formic acid and peroxymonosulfate, J. Geophys. Res.-Atmos., 91, 9807–9826, https://doi.org/10.1029/JD091iD09p09807, 1986. a
Jarecka, D., Grabowski, W. W., Morrison, H., and Pawlowska, H.: Homogeneity of the Subgrid-Scale Turbulent Mixing in Large-Eddy Simulation of Shallow Convection, J. Atmos. Sci., 70, 2751–2767, https://doi.org/10.1175/JAS-D-13-042.1, 2013. a
Jöckel, P., Tost, H., Pozzer, A., Brühl, C., Buchholz, J., Ganzeveld, L., Hoor, P., Kerkweg, A., Lawrence, M. G., Sander, R., Steil, B., Stiller, G., Tanarhte, M., Taraborrelli, D., van Aardenne, J., and Lelieveld, J.: The atmospheric chemistry general circulation model ECHAM5/MESSy1: consistent simulation of ozone from the surface to the mesosphere, Atmos. Chem. Phys., 6, 5067–5104, https://doi.org/10.5194/acp-6-5067-2006, 2006. 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
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
Jülich Supercomputing Centre: JUWELS: Modular Tier-0/1 Supercomputer at the Jülich Supercomputing Centre, J. Large-Scale Res. Fac., 5, A171, https://doi.org/10.17815/jlsrf-5-171, 2019.
Lim, Y. B., Tan, Y., and Turpin, B. J.: Chemical insights, explicit chemistry, and yields of secondary organic aerosol from OH radical oxidation of methylglyoxal and glyoxal in the aqueous phase, Atmos. Chem. Phys., 13, 8651–8667, https://doi.org/10.5194/acp-13-8651-2013, 2013. a
Lin, G., Penner, J. E., Sillman, S., Taraborrelli, D., and Lelieveld, J.: Global modeling of SOA formation from dicarbonyls, epoxides, organic nitrates and peroxides, Atmos. Chem. Phys., 12, 4743–4774, https://doi.org/10.5194/acp-12-4743-2012, 2012. a, b
Mayer, B. and Madronich, S.: Actinic flux and photolysis in water droplets: Mie calculations and geometrical optics limit, Atmos. Chem. Phys., 4, 2241–2250, https://doi.org/10.5194/acp-4-2241-2004, 2004. a
Monod, A., Chevallier, E., Jolibois, R. D., Doussin, J., Picquet-Varrault, B., and Carlier, P.: Photooxidation of methylhydroperoxide and ethylhydroperoxide in the aqueous phase under simulated cloud droplet conditions, Atmos. Environ., 41, 2412–2426, https://doi.org/10.1016/j.atmosenv.2006.10.006, 2007. a, b
Mouchel-Vallon, C., Deguillaume, L., Monod, A., Perroux, H., Rose, C., Ghigo, G., Long, Y., Leriche, M., Aumont, B., Patryl, L., Armand, P., and Chaumerliac, N.: CLEPS 1.0: A new protocol for cloud aqueous phase oxidation of VOC mechanisms, Geosci. Model Dev., 10, 1339–1362, https://doi.org/10.5194/gmd-10-1339-2017, 2017. a, b, c, d, e, f, g, h, i, j, k, l, m, n, o, p, q, r, s, t, u, v, w
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
Piesiak, A., Schuchmann, M. N., Zegota, H., and von Sonntag, C.: β-Hydroxyethylperoxyl radicals: a study of the γ-radiolysis andpulse radiolysis of ethylene in oxygenated aqueous solutions, Z. Naturforsch., 39, 1262–1267, 1984. a
Poulain, L., Katrib, Y., Isikli, E., Liu, Y., Wortham, H., Mirabel, P., Calvé, S. L., and Monod, A.: In-cloud multiphase behaviour of acetone in the troposphere: Gas uptake, Henry’s law equilibrium and aqueous phase photooxidation, Chemosphere, 81, 312–320, https://doi.org/10.1016/j.chemosphere.2010.07.032, 2010. a
Rosanka, S., Franco, B., Clarisse, L., Coheur, P.-F., Wahner, A., and Taraborrelli, D.: Organic pollutants from tropical peatland fires: regional influences and its impact on lower stratospheric ozone, Atmos. Chem. Phys. Discuss. [preprint], https://doi.org/10.5194/acp-2020-1130, in review, 2020. a
Rosanka, S., Sander, R., Franco, B., Wespes, C., Wahner, A., and Taraborrelli, D.: Oxidation of low-molecular-weight organic compounds in cloud droplets: global impact on tropospheric oxidants, Atmos. Chem. Phys., 21, 9909–9930, https://doi.org/10.5194/acp-21-9909-2021, 2021a. a, b, c, d, e, f, g, h, i, j
Rosanka, S., Sander, R., Wahner, A., and Taraborrelli, D.: Model output from CAABA/MECCA obtained during the development of JAMOC [Data set], Jülich DATA, https://doi.org/10.26165/JUELICH-DATA/SD9F6B, 2021b. a
Ruggaber, A., Dlugi, R., Bott, A., Forkel, R., Herrmann, H., and Jacobi, H.-W.: Modelling of radiation quantities and photolysis frequencies in the aqueous phase in the troposphere, Atmos. Environ., 31, 3137–3150, https://doi.org/10.1016/S1352-2310(97)00058-7,1997. a
Rumble, J. R. (Ed.): CRC Handbook of Chemistry and Physics, 101st edn., CRC Press, Boca Raton, FL, 2020. a
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, R., Baumgaertner, A., Cabrera-Perez, D., Frank, F., Gromov, S., Grooß, J.-U., Harder, H., Huijnen, V., Jöckel, P., Karydis, V. A., Niemeyer, K. E., Pozzer, A., Riede, H., Schultz, M. G., Taraborrelli, D., and Tauer, S.: The community atmospheric chemistry box model CAABA/MECCA-4.0, Geosci. Model Dev., 12, 1365–1385, https://doi.org/10.5194/gmd-12-1365-2019, 2019. a, b, c, d
Schuchmann, M. N., Zegota, H., and von Sonntag, C.: Acetateperoxyl radicals, : a study on the γ-radiolysis andpulse radiolysis of acetate in oxygenated aqueous solutions, Z. Naturforsch. Pt. B, 40, 215–221, 1985. a
Schwartz, S. E.: Mass-Transport Considerations Pertinent to Aqueous Phase Reactions of Gases in Liquid-Water Clouds, in: Chemistry of Multiphase Atmospheric Systems, edited by: Jaeschke, W., Springer Berlin Heidelberg, Berlin, Heidelberg, 415–471, 1986. a
Staehelin, J. and Hoigné, J.: Decomposition of ozone in water in the presence of organic solutes acting as promoters and inhibitors of radical chain reactions, Environ. Sci. Technol., 19, 1206–1213, https://doi.org/10.1021/es00142a012, 1985. a
Staehelin, J., Buehler, R. E., and Hoigné, J.: Ozone decomposition in water studied by pulse radiolysis. 2. Hydroxyl and hydrogen tetroxide (HO4) as chain intermediates, J. Phys. Chem., 88, 5999–6004, https://doi.org/10.1021/j150668a051, 1984. a
Tan, Y., Perri, M. J., Seitzinger, S. P., and Turpin, B. J.: Effects of Precursor Concentration and Acidic Sulfate in Aqueous Glyoxal–OH Radical Oxidation and Implications for Secondary Organic Aerosol, Environ. Sci. Technol., 43, 8105–8112, https://doi.org/10.1021/es901742f, 2009. 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
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
Tie, X., Madronich, S., Walters, S., Edwards, D. P., Ginoux, P., Mahowald, N., Zhang, R., Lou, C., and Brasseur, G.: Assessment of the global impact of aerosols on tropospheric oxidants, J. Geophys. Res.-Atmos., 110, D03204, https://doi.org/10.1029/2004JD005359, 2005. a
Tilgner, A., Bräuer, P., Wolke, R., and Herrmann, H.: Modelling multiphase chemistry in deliquescent aerosols and clouds using CAPRAM3.0i, J. Atmos. Chem., 70, 221–256, https://doi.org/10.1007/s10874-013-9267-4, 2013. a
Tost, H., Jöckel, P., Kerkweg, A., Sander, R., and Lelieveld, J.: Technical note: A new comprehensive SCAVenging submodel for global atmospheric chemistry modelling, Atmos. Chem. Phys., 6, 565–574, https://doi.org/10.5194/acp-6-565-2006, 2006. a, b
Tost, H., Jöckel, P., Kerkweg, A., Pozzer, A., Sander, R., and Lelieveld, J.: Global cloud and precipitation chemistry and wet deposition: tropospheric model simulations with ECHAM5/MESSy1, Atmos. Chem. Phys., 7, 2733–2757, https://doi.org/10.5194/acp-7-2733-2007, 2007. a, b
von Sonntag, C. and Schuchmann, H.-P.: Peroxyl Radicals in Aqueous Solutions, in: The Chemistry of Free Radicals: Peroxyl Radicals, Wiley, New York, 1997. a
Zegota, H., Schuchmann, M. N., Schulz, D., and von Sonntag, C.: Acetonylperoxyl radicals, CH3COCH2O2: A study on the γ-radiolysis and pulse radiolysis of acetone in oxygenated aqueoussolutions, Z. Naturforsch., 41, 1015–1022, 1986. a, b