Extension, development, and evaluation of the representation of the OH-initiated dimethyl sulfide (DMS) oxidation mechanism in the Master Chemical Mechanism (MCM) v3.3.1 framework

Abstract

Understanding dimethyl sulfide (DMS) oxidation can help us constrain its contribution to Earth's radiative balance. Following the discovery of hydroperoxymethyl thioformate (HPMTF) as a DMS oxidation product, a range of new experimental chamber studies have since improved our knowledge of the oxidation mechanism of DMS and delivered detailed chemical mechanisms. However, these mechanisms have not undergone formal intercomparisons to evaluate their performance.

This study aimed to synthesise the recent experimental studies and develop a new, near-explicit, DMS mechanism, through a thorough literature review. A simple box model was then used with the mechanism to simulate a series of chamber experiments and evaluated through comparison with four published mechanisms. Our modelling shows that the mechanism developed in this work outperformed the other mechanisms on average when compared to the experimental chamber data, having the lowest fractional gross error for 8 out of the 14 DMS oxidation products studied. A box model of a marine boundary layer was also run, demonstrating that the deviations in the mechanisms seen when comparing them against chamber data are also prominent under more atmospherically relevant conditions.

Although this work demonstrates the need for further experimental work, the mechanism developed in this work has been evaluated against a range of experiments, which validate the mechanism and reduce the bias from individual experiments. Our mechanism provides a good basis for a near-explicit DMS oxidation mechanism that would include other initiation reactions (e.g. halogens) and can be used to compare the performance of reduced mechanisms used in global models.

1 Introduction

Dimethyl sulfide (${\mathrm{CH}}_{\mathrm{3}}{\mathrm{SCH}}_{\mathrm{3}}$; DMS) emissions are the largest natural source of sulfur in the atmosphere . It is formed from phytoplankton when they undergo physiological stress and, to a smaller extent, from terrestrial vegetation . The presence of a sulfur atom in the molecule leads to rather complex oxidation mechanisms, much more complex than that of similar sized hydrocarbons (e.g. ). The earliest study on DMS $\mathrm{OH}$-initiated oxidation dates back to the 1970s , the decade before DMS was postulated as being involved in a global homeostatic cycle later termed the CLAW hypothesis . This (commonly refuted) hypothesis that warmer temperatures would cause phytoplankton to emit more DMS, resulting in higher concentrations of cloud condensation nuclei (CCN) and an increase in cloud formation, effectively counteracting global warming, has since led to an abundance of research . performed a comprehensive overview of the oxidation of DMS, capturing the major literature up to ca. 2006. The review from highlights the major features of DMS oxidation: (1) $\mathrm{OH}$ is the principle oxidant of DMS and can initiate oxidation via $\mathrm{H}$-atom abstraction at a methyl group or $\mathrm{OH}$ addition to the sulfur atom; (2) the nitrate radical (${\mathrm{NO}}_{\mathrm{3}}$) primarily reacts via $\mathrm{H}$-atom abstraction; and (3) halogen atoms ($\mathrm{Cl}$, $\mathrm{Br}$, and $\mathrm{I}$) and halogen oxides ($\mathrm{ClO}$, $\mathrm{BrO}$, and $\mathrm{IO}$) can participate in $\mathrm{H}$-atom abstraction, halogen-atom addition, and $\mathrm{O}$-atom addition to the sulfur atom to form dimethyl sulfoxide (${\mathrm{CH}}_{\mathrm{3}}\mathrm{S}={\mathrm{OCH}}_{\mathrm{3}}$; DMSO). Based on these initial oxidation mechanisms, a range of oxidation products are possible.

In spite of the extensive previous work on DMS, several aspects of its oxidation have remained highly uncertain, and new surprises have emerged in the last few years. One of the major oxidation products of DMS is the methylthiomethyl peroxy radical (${\mathrm{CH}}_{\mathrm{3}}{\mathrm{SCH}}_{\mathrm{2}}{\mathrm{O}}_{\mathrm{2}}$; MTMP). Based on ab initio calculations, suggested that this could undergo atmospheric autoxidation and generate hydroperoxymethyl thioformate ($\mathrm{O}={\mathrm{CHSCH}}_{\mathrm{2}}\mathrm{OOH}$; HPMTF), a mechanistic pathway that was not considered previously in the review. The existence of HPMTF in the atmosphere was conclusively established thanks to global aircraft observations . Additionally, recent work suggests that the formation and deposition of a major product of DMS oxidation, methanesulfonic acid (${\mathrm{CH}}_{\mathrm{3}}{\mathrm{SO}}_{\mathrm{3}}\mathrm{H}$; MSA), are not accurately modelled . MSA has been known to contribute to cloud condensation nuclei (CCN) due to its low volatility and hygroscopicity , making the accurate modelling of the compound important to understand the role of DMS oxidation in cloud formation.

New updates to our understanding of DMS oxidation could have important impacts on our understanding of the role of DMS in the climate system. DMS sets the natural sulfur background in most parts of the world and, as has been described previously, could play an important role in CCN. looked at the effects of updates to the gas- and aqueous-phase DMS oxidation on radiative forcing using the CESM2 (CAM6-chem) model. They found that updates to gas-phase chemistry, including HPMTF, led to significant changes in the pre-industrial aerosol burden and so the change in radiative forcing from the pre-industrial period to the present day. However, these changes were counteracted when they accounted for updates to the aqueous-phase chemistry of DMS.

The mechanism of was a significant improvement from the standard chemistry used in CESM2 (CAM6-chem) for DMS (similar for other CMIP6-era Earth system models, such as ). However, the DMS mechanism itself was never evaluated against experimental chamber data. Instead, the authors evaluated the performance of the scheme through the comparison of measured and simulated species related to DMS, an approach that can be affected by emission bias. used as their emissions of DMS the sea surface DMS climatology. Whilst this is widely regarded as a reference DMS climatology, the simulation of atmospheric DMS by was significantly high biased compared to the median of aircraft- and ground-based observations.

The best way to improve our understanding of DMS oxidation, and by proxy its impacts on the climate system, is through the development of comprehensive mechanisms evaluated against experimental data and the incorporation of these mechanisms into complex models. There have been a number of experiments recently, and a number of detailed mechanisms have been developed. , , and performed simulation chamber experiments under a wide range of conditions (such as mixing ratios of reagents, photolysis environments, temperature, and humidity) and modelled their experiments using near-explicit mechanisms that included the HPMTF pathway along with adjustments to the DMS subset of the Master Chemical Mechanism (MCM v3.3.1). In the case of the study, their mechanism was tied to the mechanism, which led to more deviations from the MCM. The MCM is often used as a starting point for mechanism development; however, the MCM DMS scheme suffers from a number of problems. Firstly, unlike the other VOCs simulated by the MCM (alkanes, alkenes, aromatics, and oxygenates), the DMS scheme never underwent evaluation against chamber experiments. Secondly, the MCM DMS scheme is rather outdated. It fails to account for the autoxidation of ${\mathrm{CH}}_{\mathrm{3}}{\mathrm{SCH}}_{\mathrm{2}}{\mathrm{O}}_{\mathrm{2}}$ and the sulfur chemistry updates from the literature captured in the most recent NASA Panel report .

Moreover, when mechanisms have been developed, they have generally been developed based on one or a small subset of experimental chamber studies. This is especially true for the recent studies focused on DMS that have solely developed mechanisms based on individual chamber studies. Given that the settings of a chamber (its volume, material, and its inputs of gases) vary significantly, it is unlikely that each independent chamber study has developed a DMS mechanism under the same set of conditions. There have also been no studies that synthesise and intercompare these detailed mechanisms.

In this study, we use the MCM DMS mechanism as a template to develop a near-explicit gas-phase DMS mechanism, focusing on the $\mathrm{OH}$-initiated chemistry. The newly developed mechanism was compared with recently developed mechanisms that have also, largely, been based on the MCM. We intercompared the results of mechanisms reported in the literature with each other and with numerical simulations of the corresponding array of chamber experiments they were derived from. In this paper, we outline the mechanism generated in this study and the other mechanisms used. We outline the chamber experiments that all mechanisms were evaluated under along with our evaluation metrics and the performance of the mechanisms. We detail the results for some key DMS oxidation products and discuss the key atmospheric implications and conclusions of the study.

2 Evaluated mechanisms

In this study, four published DMS mechanisms were compared to each other and to results from experimental chamber studies. The mechanisms were extracted from the primary studies as , , , and the MCM v3.3.1 (; hereafter Ye, Jernigan, Shen, and MCM). In addition to the four mechanisms from the literature that were evaluated, we also developed a new DMS mechanism for this study. A comparison of the major reactions included in the five mechanisms evaluated in this paper is given in Fig. .

Figure 1

The major reactions of the DMS oxidation pathway for the five mechanisms discussed in this work, with the colours of the arrows representing the mechanisms that included the respective reactions: Jernigan (blue), Ye (red), MCM (purple), Shen (orange), and the mechanism developed in this work (yellow). Only the sulfur products are included. Note that not every reaction that was included in the mechanisms is shown (see Table A1 in the Appendix and Sect. S1 in the Supplement for full details).

[Figure omitted. See PDF]

Only the sulfur reactions from the mechanisms were compared in this study. In addition to the MCM, the Ye mechanism included 12 reactions of the HPMTF pathway. The Jernigan mechanism added only five new reactions and adjusted four existing reactions. Most of the changes involved adding a simplified HPMTF pathway, along with changes to the methylthiomethyl peroxy (${\mathrm{CH}}_{\mathrm{3}}{\mathrm{SCH}}_{\mathrm{2}}{\mathrm{O}}_{\mathrm{2}}$; MTMP) radical reaction with ${\mathrm{RO}}_{\mathrm{2}}$ and the DMSO reaction with an $\mathrm{OH}$ radical. These changes to the MCM were based on the mechanism file included in the data in the Supplement of the paper, which differed slightly from the table of reactions given in their Supplement. Note that the study also included a mechanism with additional adjustments to the MCM, based on theoretical calculations of the reactions following the OH-initiated oxidation of HPMTF (multigenerational mechanism). As this mechanism was outperformed by the simpler, first-generational mechanism, the first-generational mechanism was used in this study and is referred to as the “Jernigan mechanism”. The modelling results of the multigenerational model are included in the Supplement (Figs. S3, S5, S7, S9, and S11). The Shen mechanism was based on the mechanism developed in and as such has the largest deviation from MCM; 25 sulfur reactions were added and 4 were adjusted (compared to the MCM). In addition, there were 19 non-sulfur reactions added, adjusted, or removed; however, these reactions were not included as part of the Shen mechanism used in this study but only in the auxiliary file for the experiment (this is explained in further detail in Sect. 3.1). The reactions adjusted or added to the MCM by the different mechanisms are given in the Supplement (Sect. S1).

The mechanism developed in this study was based on a thorough literature review to update and improve the DMS mechanism in the MCM. To determine which rate constants should be used in the model, the same construction methodology as the MCM was used. This methodology prioritises evaluated data (such as the NASA Panel report and IUPAC ) followed by published experimental data, structure–activity relationships, and theoretical calculations. In some cases where rate constants had not been experimentally determined, we manually adjusted them to improve the performance of the mechanism in the chamber studies. One exception was the decomposition of ${\mathrm{CH}}_{\mathrm{3}}{\mathrm{SO}}_{\mathrm{2}}$, for which the rate constant for the reaction had been experimentally determined, but there was no consensus in the literature (described in more detail in Sect. 5.4), so it was also adjusted for this study. Overall, 73 reactions were added and 21 reactions were adjusted from the base MCM mechanism. A more detailed explanation of the construction of this mechanism, along with tables of the reactions added and adjusted, is included in the Appendix.

3 Experiments studied

The experiments studied in this work were conducted under a range of different experimental conditions and by extension model input parameters. These experimental conditions are summarised in Table . To evaluate the mechanisms against the different experiments, a consistent approach was taken to deal with the different methods of modelling the experiments used by the authors of the papers.

Table 1

A summary of the experimental conditions of the experiments studied in this paper . The ${\mathrm{RO}}_{\mathrm{2}}$ and $\mathrm{OH}$ radical concentrations were found through the box modelling of the experiments (using the mechanism developed in this work).

Experiment
		experiment 1	experiment 2a
Temp (K)	290	295	295	298	263
Chamber (m${}^{\mathrm{3}}$)	0.34	7.5	7.5	0.6	26.1
$\mathrm{OH}$ source${}^{\ast }$	${\mathrm{H}}_{\mathrm{2}}{\mathrm{O}}_{\mathrm{2}}$	$\mathrm{HONO}$	${\mathrm{H}}_{\mathrm{2}}{\mathrm{O}}_{\mathrm{2}}$	TME + ${\mathrm{O}}_{\mathrm{3}}$	${\mathrm{O}}_{\mathrm{3}}$ (+ ${\mathrm{H}}_{\mathrm{2}}\mathrm{O}$)
Avg. $\mathrm{OH}$ (cm${}^{-\mathrm{3}}$)	$\mathrm{2.9}\times {\mathrm{10}}^{\mathrm{7}}$	$\mathrm{5.9}\times {\mathrm{10}}^{\mathrm{6}}$	$\mathrm{1.5}\times {\mathrm{10}}^{\mathrm{6}}$	$\mathrm{1.4}\times {\mathrm{10}}^{\mathrm{6}}$	$\mathrm{5.7}\times {\mathrm{10}}^{\mathrm{6}}$
Avg. ${\mathrm{RO}}_{\mathrm{2}}$ (cm${}^{-\mathrm{3}}$)	$\mathrm{2.0}\times {\mathrm{10}}^{\mathrm{11}}$	$\mathrm{1.6}\times {\mathrm{10}}^{\mathrm{7}}$	$\mathrm{4.7}\times {\mathrm{10}}^{\mathrm{8}}$	$\mathrm{1.8}\times {\mathrm{10}}^{\mathrm{9}}$	$\mathrm{3.2}\times {\mathrm{10}}^{\mathrm{8}}$
RH (%)	–	1	1	$<$ 0.5	70
DMS (ppb)	15 000	72.8	82	10	0.6
NO (ppb)	–	50	–	–	–
NO${}_{\mathrm{2}}$ (ppb)	–	90	–	–	–
${\mathrm{H}}_{\mathrm{2}}{\mathrm{O}}_{\mathrm{2}}$ (ppb)	25 000	–	1500	–	–
HONO (ppb)	–	90	–	–	–
CO (ppb)	–	–	–	–	120
O${}_{\mathrm{3}}$ (ppb)	–	–	–	23	125
VOC photolysis	No	Yes (300–400 nm)	Yes (300–400 nm)	No	No
Duration (h)	0.5	2	5	20	6

${}^{\ast }$ Apart from during the experiment, $\mathrm{OH}$ was formed through the photolysis of the given precursor.

3.1 Modelling of the experiments

The zero-dimensional box model BOXMOX was used in this work , which is an open-source wrapper for the Kinetic PreProcessor (KPP) software.

When first modelling an experiment that had been modelled by the respective authors, the input parameters and mechanism from that paper were used. This allowed us to directly compare our modelled output to their modelled outputs and ensure we were able to correctly replicate the chemical system. Figure  shows our replication of the output from these papers was in generally excellent agreement (compare the solid black and dashed coloured lines). In replicating the modelled outputs from the papers, some of our model outputs deviated by up to 16 % at the end of the experiment; however, the larger deviations tended to be for the minor species and in all cases were well within the deviations between the different mechanisms evaluated, providing confidence that we were able to faithfully simulate using a unified framework, the different experimental chamber studies from the , , and papers.

Figure 2

Reproduction of the modelling results from the chamber experiments simulated in this study. The modelling results of the experiments from the authors of , , and are shown in dashed coloured lines (obtained via request or through data in the supplements attached to the respective papers); our modelling using the same mechanism and input files as given is shown in solid black lines.

[Figure omitted. See PDF]

It is worth noting that special treatment was needed in modelling the experiment 1 and experiment. In the case of the experiment 1, the input concentration of ${\mathrm{NO}}_{\mathrm{2}}$ was increased to reflect the measured concentration, which resulted in less DMS being consumed. See the Supplement (Sects. S2 and S3) for details.

When comparing the mechanisms to the experiments, we only wanted to compare the effects of changing the gas-phase sulfur reactions. To do this, we kept the model inputs and mechanism consistent between the simulations of an experiment, with the exception of the sulfur reactions, which were adjusted between the mechanisms. and also adjusted non-sulfur reactions, such as the reaction between methyl peroxy (${\mathrm{CH}}_{\mathrm{3}}{\mathrm{O}}_{\mathrm{2}}$) and $\mathrm{OH}$ radicals. These non-sulfur reactions were kept in the auxiliary mechanism of an experiment to keep our modelling of an experiment from a paper consistent with the modelling in that paper by the authors. In addition, reactions for the loss from dilution, wall loss, and, in the case of , heterogeneous wall reactions forming DMSO and dimethyl sulfone (${\mathrm{CH}}_{\mathrm{3}}{\mathrm{SO}}_{\mathrm{2}}{\mathrm{CH}}_{\mathrm{3}}$; DMSO2) were also included by us in the auxiliary mechanisms for each experiment. In the auxiliary mechanism for the experiment, the tetramethylethylene ($({\mathrm{CH}}_{\mathrm{3}}{)}_{\mathrm{2}}\mathrm{C}=\mathrm{C}({\mathrm{CH}}_{\mathrm{3}}{)}_{\mathrm{2}}$; TME) subset of the MCM was included along with the adjustments from the authors. Links to these auxiliary mechanisms can be found in the “Data availability” section.

3.2 Comparison of the mechanisms

After the input parameters and auxiliary mechanisms for the experiments had been determined, the different mechanisms were compared to the chamber studies. This was done for all experiments and is included in the Supplement (Sect. S4), with the experiment results shown here as an example.

The experiment was performed in a low-${\mathrm{NO}}_{\mathrm{x}}$ environment with $\mathrm{OH}$ produced from the ozonolysis of tetramethylethylene ($({\mathrm{CH}}_{\mathrm{3}}{)}_{\mathrm{2}}\mathrm{C}=\mathrm{C}({\mathrm{CH}}_{\mathrm{3}}{)}_{\mathrm{2}}$; TME). The seven products measured in the experiment are shown in Fig. , with the modelling outputs of the Jernigan, Shen, Ye, and MCM mechanisms, along with the mechanism developed in this study.

Figure 3

Comparison of the products measured in the experiment (solid grey lines) to our modelling results using the Jernigan (dotted–dashed blue lines), Ye (dashed red lines), MCM (dotted–dashed purple lines), and Shen mechanisms (dotted orange lines), along with the mechanism developed in this work (dotted–dashed yellow lines).

[Figure omitted. See PDF]

As shown in Fig.  (and Figs. S3–S12 in the Supplement), the performance of the DMS mechanisms deviated greatly from each other and the experiment(s), especially in the low- ${\mathrm{NO}}_{\mathrm{x}}$ conditions of the experiment. Some products, such as thioperformic acid ($\mathrm{S}=\mathrm{CHOOH}$, TPA), were not included by most mechanisms, while others, such as carbonyl sulfide ($\mathrm{OCS}$), were produced in much lower concentrations by some of the mechanisms than was observed experimentally.

3.3 Evaluation metrics used

The mechanisms in this paper are evaluated on the basis of three metrics: fractional gross error (FGE), modified mean bias (MMB), and correlation ($\mathit{\rho }$). These metrics are normalised and thus independent of the units used or the relative intensities of the species. This makes it easier to compare across different experiments, which use different units, and a large range of precursor concentrations (the initial DMS mole fractions used in the experiments ranged from 0.6–15 000 ppb). The measured and modelled concentrations over the time steps of the experiment were utilised for these calculations, and the modelled outputs were interpolated to the observation time steps, with any gaps in the observations not included in the analysis.

The MMB (given in the equation below) is a normalised version of mean bias. This metric provides the bias between the model and the observations within a range of $-\mathrm{2}$ (negative bias) to 2 (positive bias), with 0 being ideal, where the model matches the observation.

1$$\mathrm{MMB}={\displaystyle \frac{\mathrm{2}}{n}}\sum _i{\displaystyle \frac{M_i-O_i}{M_i+O_i}},$$ where $\mathit{M}=(M_i{)}_{\mathrm{0}\le i\le n}$ is a vector of modelled values and $\mathit{O}=(O_i{)}_{\mathrm{0}\le i\le n}$ is a vector of the observed (experimental) values.

The FGE is the normalised version of the mean error. The FGE measures the error in the model within a range of 0 to 2, with 0 being ideal, where the model matches the observation. 2$$\mathrm{FGE}={\displaystyle \frac{\mathrm{2}}{n}}\sum _i\left|{\displaystyle \frac{M_i-O_i}{M_i+O_i}}\right|$$

Finally, for correlation, as the data were not normally distributed, Spearman's rank correlation coefficient was used. Spearman's rank correlation coefficient, referred to as just “correlation” or $\mathit{\rho }$ in this study, is a measure of how linearly correlated the ranks of the measured and modelled values are. The values range from $-\mathrm{1}$ (negatively correlated) to 1 (positively correlated), with 1 being perfectly correlated. 3$$\mathit{\rho }={\displaystyle \frac{\mathrm{cov}\left(R\right(\mathit{O}),R(\mathit{M}\left)\right)}{{\mathit{\sigma }}_{R\left(\mathit{O}\right)}{\mathit{\sigma }}_{R\left(\mathit{M}\right)}}},$$ where $R\left(\mathit{M}\right)$ is the rank of the modelled values and $R\left(\mathit{O}\right)$ is the rank of the observed (experimental) values.

4 Intercomparison and evaluation of recent DMS mechanisms

The modified mean bias (MMB), fractional gross error (FGE), and correlation ($\mathit{\rho }$) for the species in all the experiments we modelled (with all the mechanisms) can be found in the Supplement (Sect. S4); however, the outcomes of the study are summarised in Fig. 4.

Figure  shows how the mechanisms perform for each species. In the case where more than one experiment measured a certain species, the metrics for that species were averaged between the different experiments. For example, as DMSO was measured by all the experiments, to obtain the average MMB for DMSO for the Shen mechanism, the MMB from the DMSO modelled in each of the experiments with the Shen mechanism were averaged. The individual MMBs for DMSO found for each experiment modelled with the Shen mechanism are also shown in Fig. , with different symbols corresponding to different experiments. This was done for all species and all mechanisms.

Figure 4

The average fractional gross error, correlation, and modified mean bias of the five mechanisms (Jernigan, Shen, Ye, MCM, and this work) for each product found in the experiments .

[Figure omitted. See PDF]

We note that none of the statistical metrics we calculated are perfect descriptors for the performance of a mechanism against experimental data. However, the use of these metrics provides a succinct and quantitative way for the evaluation of the mechanisms to be performed. An idealised mechanism would have an average FGE and MMB of 0 for all compounds and a $\mathit{\rho }$ of 1. Figure  demonstrates that none of the mechanisms are ideal, with the performance of the mechanisms differing over the different compounds. The metrics from the individual experiments (shown as the fainter symbols) also show a large range between different experiments, even for the same mechanism. This spread is demonstrated with methyl thioformate (${\mathrm{CH}}_{\mathrm{3}}\mathrm{SCHO}$; MTF), where the Jernigan mechanism shows a negative bias of around $-\mathrm{1}$ when modelling the experiment and a positive bias of 1.6 for the experiment. Although the average bias of the Jernigan mechanism for MTF of 0.16 is close to the ideal bias of 0, the deviations from the experiments are captured in the larger average error of 0.96. This demonstrates that although the MMB can be used to assess the performance of a mechanism and provides more information than the FGE, an average bias of close to zero may not be indicative of a mechanism performing well in replicating the production of a compound; the range of the bias from different experiments should be taken into account. However, the average FGE can help summarise the overall performance of each mechanism.

For 8 of the 14 species included in this study, the mechanism developed in this work had the lowest average FGE. For 12 species, the mechanism from this study is among the two mechanisms with the lowest error. The average FGE for the mechanism was 0.70, with most species having an error lower than 1. The MCM has the highest or equal highest average error for eight of the species, with four of those species not being formed by the MCM (resulting in an FGE of 2). The high errors from the MCM tended to come from its poor performance for the experiments conducted in low-${\mathrm{NO}}_{\mathrm{x}}$ conditions. These experiments were where HPMTF formation dominates, a pathway missing from the MCM.

The large spread in the error and bias between the mechanisms demonstrates that the adjustments between the mechanisms, sometimes as few as 9 or 12 reactions in the case of the Jernigan and Ye mechanisms, are important in the modelling of these experiments. Since the mechanisms being adjusted are compared to only one set of experiments, they tend to perform better for their own experiment compared to the others. This could be due to the experimental conditions impacting which reactions are important or the different products measured in the experiments. One such example is the Shen mechanism having the largest average error for ${\mathrm{SO}}_{\mathrm{2}}$ (0.76). The Shen mechanism underestimates the ${\mathrm{SO}}_{\mathrm{2}}$ in the and experiments, resulting in an average bias of $-\mathrm{0.73}$. Although ${\mathrm{SO}}_{\mathrm{2}}$ was expected as a major product, it was not measured in the experiment, and as such the mechanism could not have been evaluated for its performance of modelling ${\mathrm{SO}}_{\mathrm{2}}$. This highlights the need for mechanism development to include a range of mechanisms and experiments, as is done in this study.

Spearman's rank correlation coefficient ($\mathit{\rho }$) can be used to assess the correlation between the mechanisms and the experiments, with two caveats. This form of correlation is a measure of whether the observed values and modelled values are both increasing or decreasing during the same time step. However, in one time step, if the increase/decrease in the observed value due to noise is larger than the actual increase/decrease in the concentration of a compound, this will affect the $\mathit{\rho }$ found. The significant noise in the concentration of some products in the experiment, the $\mathrm{OCS}$ concentration in the experiment, and the ${\mathrm{SO}}_{\mathrm{2}}$ concentration in the experiment 2a contribute to lower $\mathit{\rho }$. Additionally, as the and experiments are steady-state experiments, once the experiment reaches a steady state, small deviations in the experiment can result in lower $\mathit{\rho }$ values. However, the reduction in $\mathit{\rho }$ due to noise and the experiments reaching a steady state will affect the performance of all mechanisms similarly, and the range in correlation found between the mechanisms for each compound can still be used to assess the performance of a mechanism. For 8 of the 14 compounds, the mechanism developed in this work has the highest or equal highest correlation, with the mechanism having the highest or second-highest correlation for 10 compounds.

In an ideal world, a developed mechanism could approach the “perfect” FGE and MMB of 0 and $\mathit{\rho }$ of 1. Deviations from the ideal can be attributed to uncertainties and unknowns in the rate constants and reactions of the mechanism, although that would assume that the experiments themselves represent the “truth”. In reality, there are uncertainties in the concentrations of the products, especially in the case of the low concentrations measured in the experiments and the difficulty in determining the sensitivities of the species measured. The compounds measured using chemical ionisation mass spectrometry (CIMS), such as HPMTF and ${\mathrm{CH}}_{\mathrm{3}}{\mathrm{SO}}_{\mathrm{4}}{\mathrm{NO}}_{\mathrm{2}}$, are affected by the largest uncertainties, with including a 50 % relative standard deviation to the species they measured by I${}^{\mathrm{-}}$ CIMS. In the case of , in addition to the large uncertainties in the experimentally determined sensitivity of HPMTF, thioperformic acid ($\mathrm{S}=\mathrm{CHOOH}$; TPA) and methanesulfinic acid (${\mathrm{CH}}_{\mathrm{3}}{\mathrm{SO}}_{\mathrm{2}}\mathrm{H}$; MSIA) sensitivities were determined by calculating the species binding energy to iodine and then comparing them to the binding energy of HPMTF to then scale the experimentally determined HPMTF sensitivity. These uncertainties again emphasise the importance of comparing multiple experiments from different sources when developing and evaluating a mechanism, such as in this study.

5 Discussion of key products

We now focus on a subset of DMS oxidation products. These products (DMSO, HPMTF, MSA, and ${\mathrm{SO}}_{\mathrm{2}}$) were chosen as they are found in field studies and were modelled differently by the mechanisms.

5.1 Dimethyl sulfoxide

Firstly, we evaluate the performance of mechanisms in simulating DMSO. DMSO is a primary oxidation product of DMS (formed from both $\mathrm{OH}$-addition and halogen oxide reactions; ). The modelled DMSO from most mechanisms was similar, with the exception of the Shen mechanism. The Shen mechanism is based on the mechanism, which uses the explicit mechanism for the $\mathrm{OH}$ addition to DMS; the addition of $\mathrm{OH}$ to DMS is reversible, forming ${\mathrm{CH}}_{\mathrm{3}}\mathrm{S}\left(\mathrm{OH}\right){\mathrm{CH}}_{\mathrm{3}}$, which can react with ${\mathrm{O}}_{\mathrm{2}}$ irreversibly to form HODMSO2. These reactions, along with their recommended rate constants from the NASA Panel report, are included in Table .

Table 2

The reversible $\mathrm{OH}$-addition reaction of DMS, along with the addition of ${\mathrm{O}}_{\mathrm{2}}$, with rate constants from the NASA Panel report.

Reaction				Rate constant
1	DMS $+\mathrm{OH}$	$\to$	${\mathrm{CH}}_{\mathrm{3}}\mathrm{S}\left(\mathrm{OH}\right){\mathrm{CH}}_{\mathrm{3}}$	$\mathrm{3}\times {\mathrm{10}}^{-\mathrm{31}}\times {\frac{T}{\mathrm{298}}}^{-\mathrm{6.24}}\times \left[M\right]$
2	${\mathrm{CH}}_{\mathrm{3}}\mathrm{S}\left(\mathrm{OH}\right){\mathrm{CH}}_{\mathrm{3}}+{\mathrm{O}}_{\mathrm{2}}$	$\to$	HODMSO2	$\mathrm{8.5}\times {\mathrm{10}}^{-\mathrm{13}}$
3	${\mathrm{CH}}_{\mathrm{3}}\mathrm{S}\left(\mathrm{OH}\right){\mathrm{CH}}_{\mathrm{3}}$	$\to$	DMS $+\mathrm{OH}$	$\frac{\mathrm{3}\times {\mathrm{10}}^{-\mathrm{31}}\times {\frac{T}{\mathrm{298}}}^{-\mathrm{6.24}}\times \left[M\right]}{\mathrm{9.6}\times {\mathrm{10}}^{-\mathrm{27}}\times e^{\mathrm{5376}/T}}$

The MCM, along with the Jernigan and Ye mechanisms, combines the three reactions from Table  into one, using the combined rate constant from IUPAC, $\mathrm{2.2}\times {\mathrm{10}}^{-\mathrm{12}}$ cm${}^{\mathrm{3}}$ molecules${}^{-\mathrm{1}}$ s${}^{-\mathrm{1}}$ at 298 K and 1 atm . In the paper, they use the same combined rate constant; however, they use it for the forward (reversible) reaction, referencing the MCM. This combined rate constant is smaller than the forward reaction recommended by the NASA Panel report ($\mathrm{7.4}\pm \mathrm{3.3}\times {\mathrm{10}}^{-\mathrm{12}}$ cm${}^{\mathrm{3}}$ molecules${}^{-\mathrm{1}}$ s${}^{-\mathrm{1}}$ at 1 atm and 298 K), as the combined reaction takes into account the backward reaction. The rate constant used by for the backward reaction ($\mathrm{2.3}\times {\mathrm{10}}^{\mathrm{6}}$ s${}^{-\mathrm{1}}$ at 298 K) is from , which is smaller than the backward reaction from the NASA Panel report ($\mathrm{1.1}\times {\mathrm{10}}^{\mathrm{7}}$ s${}^{-\mathrm{1}}$). However, due to the fast reaction of ${\mathrm{CH}}_{\mathrm{3}}\mathrm{S}\left(\mathrm{OH}\right){\mathrm{CH}}_{\mathrm{3}}$ with O${}_{\mathrm{2}}$, the smaller forward reaction rate constant used by (and the Shen mechanism) results in less DMSO being produced, which is why less DMSO is formed via the Shen mechanism.

5.2 Hydroperoxymethyl thioformate

Recent global modelling points to HPMTF being a major oxidation product of DMS that was unaccounted for until very recently. The major uncertainties surrounding the modelling of HPMTF are the first isomerisation ($\mathrm{H}$ shift of ${\mathrm{CH}}_{\mathrm{3}}{\mathrm{SCH}}_{\mathrm{2}}{\mathrm{O}}_{\mathrm{2}}$), photolysis of HPMTF, the reactions of ${\mathrm{CH}}_{\mathrm{3}}{\mathrm{SCH}}_{\mathrm{2}}{\mathrm{O}}_{\mathrm{2}}$, and the uptake of HPMTF onto aerosol surfaces . HPMTF was measured in the experiment, the experiment 2a, and the experiment; the observed HPMTF along with the HPMTF modelled by the various mechanisms is shown in Fig. .

Figure 5

Comparison of modelled and measured HPMTF in the experiment 2a, the experiment, and the experiment. The measured HPMTF is shown in grey, and modelled HPMTF, using the various mechanisms, is shown using the same colours and line styles as Fig. .

[Figure omitted. See PDF]

In the mechanism developed in this work, the temperature-dependent rate constant was used for the first $\mathrm{H}$ shift ($\mathrm{0.06}\pm \mathrm{0.02}$ s${}^{-\mathrm{1}}$ at 298 K). This rate constant was used as it was both measured directly and is temperature-dependent; however, it is smaller than the other rate constants measured at 298 K . The HPMTF formed is sensitive to this reaction. In the experiment, our new mechanism outputs the second-lowest concentration of HPMTF, forming 0.2 % more HPMTF than the Jernigan mechanism. However, the Jernigan mechanism includes a larger rate constant for the ${\mathrm{RO}}_{\mathrm{2}}$ reaction of ${\mathrm{CH}}_{\mathrm{3}}{\mathrm{SCH}}_{\mathrm{2}}{\mathrm{O}}_{\mathrm{2}}$, which reduces the HPMTF formed for that experiment.

The rate constant used by is the rate constant recommended by the NASA Panel report for the ${\mathrm{CH}}_{\mathrm{3}}{\mathrm{SCH}}_{\mathrm{2}}{\mathrm{O}}_{\mathrm{2}}$ self-reaction ($\mathrm{1.0}\pm \mathrm{0.6}\times {\mathrm{10}}^{-\mathrm{11}}$ cm${}^{\mathrm{3}}$ molecules${}^{-\mathrm{1}}$ s${}^{-\mathrm{1}}$). All the other mechanisms use the rate constant from the MCM v3.3.1 ($\mathrm{3.7}\times {\mathrm{10}}^{-\mathrm{12}}$ cm${}^{\mathrm{3}}$ molecules${}^{-\mathrm{1}}$ s${}^{-\mathrm{1}}$), which uses the same self-reaction rate constant recommended by the NASA Panel report. Since the MCM uses a pooled ${\mathrm{RO}}_{\mathrm{2}}$ concentration instead of explicitly representing the ${\mathrm{RO}}_{\mathrm{2}}$ reactions, the rate constant they use for ${\mathrm{RO}}_{\mathrm{2}}$ reactions is double the geometric mean of the self-reaction rate constant of the species in question and the self-reaction rate constant for ${\mathrm{CH}}_{\mathrm{3}}{\mathrm{O}}_{\mathrm{2}}$ at 298 K, $\mathrm{3.5}\pm \mathrm{1.5}\times {\mathrm{10}}^{-\mathrm{13}}$ cm${}^{\mathrm{3}}$ molecules${}^{-\mathrm{1}}$ s${}^{-\mathrm{1}}$ . This method is used as ${\mathrm{CH}}_{\mathrm{3}}{\mathrm{O}}_{\mathrm{2}}$ is generally the major ${\mathrm{RO}}_{\mathrm{2}}$ reacting; in the case of these experiments, ${\mathrm{CH}}_{\mathrm{3}}{\mathrm{SCH}}_{\mathrm{2}}{\mathrm{O}}_{\mathrm{2}}$ makes up 1 %–30 % of the ${\mathrm{RO}}_{\mathrm{2}}$ pool, whereas ${\mathrm{CH}}_{\mathrm{3}}{\mathrm{O}}_{\mathrm{2}}$ makes up 40 %–90 % of the ${\mathrm{RO}}_{\mathrm{2}}$ pool (based on our modelling results). use a higher rate constant for the ${\mathrm{CH}}_{\mathrm{3}}{\mathrm{SCH}}_{\mathrm{2}}{\mathrm{O}}_{\mathrm{2}}$ reaction with pooled ${\mathrm{RO}}_{\mathrm{2}}$, which is why the Jernigan mechanism forms less HPMTF than the other mechanisms in their experiment (due to the shorter lifetime of ${\mathrm{CH}}_{\mathrm{3}}{\mathrm{SCH}}_{\mathrm{2}}{\mathrm{O}}_{\mathrm{2}}$). As the experiment involved the ozonolysis of TME, it had a higher concentration of total ${\mathrm{RO}}_{\mathrm{2}}$ compared to the experiment 2a and the experiment by 383 % and 563 %, respectively. The larger ${\mathrm{RO}}_{\mathrm{2}}$ concentration meant that the different rate constant used for the ${\mathrm{RO}}_{\mathrm{2}}$ reaction of ${\mathrm{CH}}_{\mathrm{3}}{\mathrm{SCH}}_{\mathrm{2}}{\mathrm{O}}_{\mathrm{2}}$ was more significant for the experiment.

Another source of difference between the mechanisms is the rate constant used for the $\mathrm{OH}$-initiated oxidation of HPMTF. The Shen mechanism uses the smallest rate constant, $\mathrm{1.4}\times {\mathrm{10}}^{-\mathrm{12}}$ cm${}^{\mathrm{3}}$ molecules${}^{-\mathrm{1}}$ s${}^{-\mathrm{1}}$, which was based on the computational paper by . The rate constant used, 1.4 (0.27–$\mathrm{2.4})\times {\mathrm{10}}^{-\mathrm{11}}$ cm${}^{\mathrm{3}}$ molecules${}^{-\mathrm{1}}$ s${}^{-\mathrm{1}}$, was based on the best fit to their experiment; however, that fit is dependent on the reactions forming HPMTF. used $\mathrm{1.0}\times {\mathrm{10}}^{-\mathrm{11}}$ cm${}^{\mathrm{3}}$ molecules${}^{-\mathrm{1}}$ s${}^{-\mathrm{1}}$ based on , who found their rate constant through the best fit to observations, and the assumption that the rate constant will be similar to the rate constant measured for methyl thioformate (${\mathrm{CH}}_{\mathrm{3}}\mathrm{SCHO}$, MTF) due to structural similarities. The rate constant used in this work, $\mathrm{1.75}\times {\mathrm{10}}^{-\mathrm{11}}$ cm${}^{\mathrm{3}}$ molecules${}^{-\mathrm{1}}$ s${}^{-\mathrm{1}}$, is an average of the rate constant obtained by ($\mathrm{2.1}\pm \mathrm{0.1}\times {\mathrm{10}}^{-\mathrm{11}}$ cm${}^{\mathrm{3}}$ molecules${}^{-\mathrm{1}}$ s${}^{-\mathrm{1}}$, found by looking at the decay of HPMTF after adding NO) and the value of the best fit from the study.

The importance of the rate constant used for the isomerisation of ${\mathrm{HOOCH}}_{\mathrm{2}}{\mathrm{SCH}}_{\mathrm{2}}{\mathrm{O}}_{\mathrm{2}}$ into HPMTF is dependent on the other reactions of ${\mathrm{HOOCH}}_{\mathrm{2}}{\mathrm{SCH}}_{\mathrm{2}}{\mathrm{O}}_{\mathrm{2}}$. Other than isomerisation, the two reactions included in most mechanisms are the reaction with $\mathrm{NO}$ and the reaction with ${\mathrm{HO}}_{\mathrm{2}}$. In addition, in our mechanism we included the reactions with ${\mathrm{RO}}_{\mathrm{2}}$ and ${\mathrm{NO}}_{\mathrm{3}}$. For the isomerisation of ${\mathrm{HOOCH}}_{\mathrm{2}}{\mathrm{SCH}}_{\mathrm{2}}{\mathrm{O}}_{\mathrm{2}}$, both the Shen mechanism and the mechanism from this work use the rate constant calculated by , which at 298 K is a factor of 317 smaller than calculated in the study (used in the Ye mechanism) and a factor of 5.6 smaller than calculated in the study. The calculation by for the first $\mathrm{H}$ shift agrees well with the measured rate constant by at 298 K, which is why the calculation for the second $\mathrm{H}$ shift was chosen for this work. However, the smaller rate constant used in this work results in the reaction of ${\mathrm{HOOCH}}_{\mathrm{2}}{\mathrm{SCH}}_{\mathrm{2}}{\mathrm{O}}_{\mathrm{2}}$ with ${\mathrm{HO}}_{\mathrm{2}}$ and, to a lesser extent, ${\mathrm{RO}}_{\mathrm{2}}$, reducing the amount of HPMTF formed.

The rate constants used in the mechanism developed in this work seem to model the HPMTF from the and experiments the best. However, apart from the MCM, it has the highest fractional gross error for the experiment 2a of all the mechanisms studied, although an uncertainty of around 50 % was included for the HPMTF measured in that experiment.

5.3 Methane sulfonic acid

MSA is measured in two experiments, namely the experiment and the experiment 1. In these experiments, the modelled MSA from the mechanism developed in this work came from two different reactions involving ${\mathrm{CH}}_{\mathrm{3}}{\mathrm{SO}}_{\mathrm{3}}$. In the experiment, nearly all of the MSA that was modelled came from the reaction with ${\mathrm{HO}}_{\mathrm{2}}$; however, in the experiment 1, MSA came from the reaction of ${\mathrm{CH}}_{\mathrm{3}}{\mathrm{SO}}_{\mathrm{3}}$ with DMS.

Apart from the mechanism developed in this work, for the experiment 1, the other mechanisms produced a modified mean bias of around $-\mathrm{2}$ for MSA (meaning almost no MSA formed). The paper, which involved a review and an evaluation of a DMS mechanism, included a few reactions forming MSA where the ${\mathrm{CH}}_{\mathrm{3}}{\mathrm{SO}}_{\mathrm{3}}$ radical abstracted hydrogen from different species. mention that type of reaction as key for MSA formation; the MCM already includes the ${\mathrm{CH}}_{\mathrm{3}}{\mathrm{SO}}_{\mathrm{3}}$ and ${\mathrm{HO}}_{\mathrm{2}}$ reaction, which is the source of MSA in the experiment and discussed below. The estimated rate constants for ${\mathrm{CH}}_{\mathrm{3}}{\mathrm{SO}}_{\mathrm{3}}$ $+$ R-H from were based on the bond dissociation energy of the relevant bond between hydrogen and the $\mathrm{H}$ donor. We included all the MSA-forming reactions from the paper; however, for this experiment, the model was only sensitive to the ${\mathrm{CH}}_{\mathrm{3}}{\mathrm{SO}}_{\mathrm{3}}$ reaction with DMS. The rate constant for this reaction was increased by a factor of 2.1 until the ratio of sulfuric acid (${\mathrm{H}}_{\mathrm{2}}{\mathrm{SO}}_{\mathrm{4}}$) to MSA was the same as measured in the experiment (as ${\mathrm{H}}_{\mathrm{2}}{\mathrm{SO}}_{\mathrm{4}}$, measured as sulfate in this experiment, is the other major fate of the ${\mathrm{CH}}_{\mathrm{3}}{\mathrm{SO}}_{\mathrm{3}}$ radical). Not only did the addition of this reaction to the mechanism explain the formation of MSA, but it also accounted for the total amount of DMS reacted with, which was underestimated by the other mechanisms (Fig. ). However, the formation of MSA is dependent on the formation of ${\mathrm{CH}}_{\mathrm{3}}{\mathrm{SO}}_{\mathrm{3}}$ in the model, which depends on ${\mathrm{CH}}_{\mathrm{3}}{\mathrm{SO}}_{\mathrm{2}}$ reactions. The rate constant of ${\mathrm{CH}}_{\mathrm{3}}{\mathrm{SO}}_{\mathrm{2}}$ decomposition was increased to 6 s${}^{-\mathrm{1}}$ (at 295 K) in our mechanism, which is discussed in more detail in the following section.

Figure  shows that following these adjustments, the mechanism from this work is able to reproduce the major products (MSA, ${\mathrm{SO}}_{\mathrm{2}}$, and ${\mathrm{H}}_{\mathrm{2}}{\mathrm{SO}}_{\mathrm{4}}$) along with the DMS lost during the experiment. Again, this method of tuning is dependent on the accuracy of the measurements of the different species and the rate constants of other reactions and is not necessarily accurate. However, it does indicate that more studies should be conducted on the ${\mathrm{CH}}_{\mathrm{3}}{\mathrm{SO}}_{\mathrm{2}}$ and ${\mathrm{CH}}_{\mathrm{3}}{\mathrm{SO}}_{\mathrm{3}}$ reactions, which include large uncertainties, as the modelled DMS oxidation products are sensitive to them.

Figure 6

The DMS, ${\mathrm{SO}}_{\mathrm{2}}$, MSA, and ${\mathrm{H}}_{\mathrm{2}}{\mathrm{SO}}_{\mathrm{4}}$ measured in the experiment 1, along with the modelling outputs from the different mechanisms (Jernigan, Ye, MCM, Shen, and the mechanism developed in this work). The experimental DMS shown has not been corrected for dilution, which was included in the modelled DMS (explained in more detail in Sect. S3 in the Supplement).

[Figure omitted. See PDF]

Our mechanism overestimates MSA in the experiment, resulting in a modified mean bias of 1.09. The source for the overestimation is likely uncertainties in the DMSO and MSIA reactions with $\mathrm{OH}$ at 263 K, along with the reaction of ${\mathrm{CH}}_{\mathrm{3}}{\mathrm{SO}}_{\mathrm{3}}$ and ${\mathrm{HO}}_{\mathrm{2}}$, indicating a need for further temperature-dependent experiments. The modelling of MSA for the experiment is discussed in further detail in the Supplement (Sect. S5).

5.4 Sulfur dioxide

The modelling of sulfur dioxide (${\mathrm{SO}}_{\mathrm{2}}$) is complex as it forms from a range of different reactions, which are dependent on the conditions of the experiments. The first column of Fig.  shows the ${\mathrm{SO}}_{\mathrm{2}}$ measured by the different experiments and the model output from the mechanisms. The second column of the figure shows the rate of formation of ${\mathrm{SO}}_{\mathrm{2}}$ by the major reactions from the mechanism developed in this study.

Figure 7

The column on the left shows the ${\mathrm{SO}}_{\mathrm{2}}$ formed in the various experiments compared to the modelled ${\mathrm{SO}}_{\mathrm{2}}$ from the different mechanisms. The column on the right shows the different rates of ${\mathrm{SO}}_{\mathrm{2}}$ formation modelled from our mechanism for those experiments; the reactants of those reactions are included in the legend. Only the major ${\mathrm{SO}}_{\mathrm{2}}$-forming reactions are included.

[Figure omitted. See PDF]

Our mechanism generally performs similarly to the Ye and Jernigan mechanisms for the total ${\mathrm{SO}}_{\mathrm{2}}$ formed in the experiments modelled. The MCM overestimates the ${\mathrm{SO}}_{\mathrm{2}}$ formed in the experiment, which is due to HPMTF being a major product in this experiment and a product missing in the MCM. The Shen mechanism tends to underestimate ${\mathrm{SO}}_{\mathrm{2}}$ in all experiments apart from the experiment. The largest deviations from the other mechanisms are in the low-${\mathrm{NO}}_{\mathrm{x}}$ experiments where HPMTF is a major product; this is mostly due to the Shen mechanism having the smallest rate constant for the reaction of HPMTF and $\mathrm{OH}$ radicals, of which ${\mathrm{SO}}_{\mathrm{2}}$ is a secondary product.

Figure  demonstrates that the formation of ${\mathrm{SO}}_{\mathrm{2}}$ from the decomposition of ${\mathrm{CH}}_{\mathrm{3}}{\mathrm{SO}}_{\mathrm{2}}$ is a major reaction for all of the experiments. The Ye and Jernigan mechanisms use the same temperature-dependent rate constant as the MCM (0.4 s${}^{-\mathrm{1}}$ at 298 K), whereas the Shen mechanism uses a larger temperature-dependent rate constant (7.0 s${}^{-\mathrm{1}}$ at 298 K) for the decomposition of ${\mathrm{CH}}_{\mathrm{3}}{\mathrm{SO}}_{\mathrm{2}}$. The rate constant used by the MCM is lower than the experimental upper bound estimated by at 1 s${}^{-\mathrm{1}}$ at 100–660 Torr and 300 K. determined the barrier for ${\mathrm{CH}}_{\mathrm{3}}{\mathrm{SO}}_{\mathrm{2}}$ dissociation via velocity map imaging to be 14 kcal mol${}^{-\mathrm{1}}$ and calculated a high-pressure ($P=\mathrm{\infty }$) rate constant of $\mathrm{1}\times {\mathrm{10}}^{\mathrm{3}}$ s${}^{-\mathrm{1}}$. Using a UCCSD(T)//UCCSD-level calculation, also calculated an energy barrier of 14 kcal mol${}^{-\mathrm{1}}$ but calculated a rate constant of $\mathrm{2}\times {\mathrm{10}}^{\mathrm{3}}$ s${}^{-\mathrm{1}}$. More recently, in their box modelling, used a rate constant of 20 s${}^{-\mathrm{1}}$ to replicate their experiments. Due to the wide range of rate constants estimated and the sensitivity of the reaction in the experiment 1, the rate constant of the decomposition of ${\mathrm{CH}}_{\mathrm{3}}{\mathrm{SO}}_{\mathrm{2}}$ was adjusted in our mechanism until the formation of MSA (and the loss of DMS) in the model matched the experiment (6 s${}^{-\mathrm{1}}$ at 295 K).

Although a larger rate constant was used for the decomposition of ${\mathrm{CH}}_{\mathrm{3}}{\mathrm{SO}}_{\mathrm{2}}$ compared to the other mechanisms in the study, Fig.  shows that this tuning did not seem to negatively affect the modelling of the experiments. The range of experiments we analysed show the multiple pathways through which ${\mathrm{SO}}_{\mathrm{2}}$ is formed and that our mechanism performs well in all the conditions studied.

6 Atmospheric implications

6.1 Marine boundary layer run

We have shown that in simulating the chamber studies, there was a large range in the performance of the mechanisms applied. Although some of the above experiments were performed to simulate realistic, marine conditions, high concentrations of DMS or $\mathrm{OH}$ precursors result in different conditions from those found in the marine environment. To compare the mechanisms to marine conditions, and determine the atmospheric implications for the divergence amongst the mechanisms, the remote marine boundary layer box model run by was used. The same input parameters as were used, with the exception of photolysis parameters, where the zenith angle was used to obtain the photolysis rate constants based on the $l$, $m$, and $n$ photolysis parameters used by the MCM . The runs were performed over 8 cloud-free days, in low-${\mathrm{NO}}_{\mathrm{x}}$ conditions (around 10 ppt) and with a zenith angle of zero during solar noon. The planetary boundary layer height (based on ) and temperature were varied throughout the day. The input parameters for the run are included in the Supplement (Sect. S6).

The results of the marine boundary layer modelling, shown in Fig. , demonstrate that there is still a significant spread between the different mechanisms under these more atmospherically relevant conditions.

Figure 8

The distribution of DMS oxidation products from the marine boundary layer run for each mechanism based on the average concentration of the products over the last 2 d of the run. The data from this figure are included as percentages in the Supplement (Table S8).

[Figure omitted. See PDF]

All mechanisms show that ${\mathrm{SO}}_{\mathrm{2}}$ is the major product formed, but the range in the fraction of ${\mathrm{SO}}_{\mathrm{2}}$ varies from 0.32 to 0.75 (more than a factor of 2). This result is important. Using our new mechanism, which we have demonstrated performs best against the range of experimental chamber studies evaluated, we show that there are much more diverse sets of products formed under atmospheric conditions than most mechanisms would predict. Our calculations imply that the use of the mechanism would result in HPMTF being the major gas-phase oxidation product of DMS. Our results are in best agreement with the results of , but we note that more detailed observational studies would be required to determine if this wider spread of DMS oxidation products simulated with our mechanism is also seen in reality. We note recent reports of significant amounts of DMSO2 , which the MCM in particular does not predict would form under the conditions investigated but that we calculate would account for approximately 2.6 % of DMS oxidation products.

7 Conclusions

The oxidation of DMS is complex but key in understanding the climate impacts of the major natural source of sulfur in the atmosphere. In this work, we used the MCM v3.3.1 DMS oxidation scheme as a template to further develop, performed a comprehensive evaluation against an array of recent DMS chamber experiments, and benchmarked the ability of recently proposed DMS mechanisms to simulate this array of experimental data. Basic statistical metrics were applied to determine the ability of our new mechanism alongside the existing mechanisms to simulate the experimental data. Based on an analysis of these statistics, we concluded that our new mechanism shows greater overall skill in simulating DMS oxidation than the other mechanisms studied. The worst-performing mechanism overall was the MCM, mostly due to the lack of the HPMTF pathway.

However, this work shows that more experimental work needs to be done to reduce the uncertainty in some of the key reactions involved in DMS oxidation. This is especially the case for the rate constants that we adjusted as the model was sensitive to them, but they had not been explored experimentally or computationally, such as ${\mathrm{ROONO}}_{\mathrm{2}}$ decomposition reactions and direct $\mathrm{OCS}$ and DMSO2 formation from the reactions of HPMTF and DMSO with $\mathrm{OH}$, respectively. Additionally, although the decomposition of ${\mathrm{CH}}_{\mathrm{3}}{\mathrm{SO}}_{\mathrm{2}}$ has been experimentally determined previously, there is no consensus in the literature for the decomposition rate constant , and more experiments should be done to constrain the reaction. Finally, the modelling of MSA in the experiment indicated that further experiments exploring the rate constants for the reactions of DMSO and MSIA with $\mathrm{OH}$ radicals at lower temperatures, along with the reaction of ${\mathrm{CH}}_{\mathrm{3}}{\mathrm{SO}}_{\mathrm{3}}$ and ${\mathrm{HO}}_{\mathrm{2}}$, could improve the modelling of MSA in the marine environment.

Additionally, $\mathrm{OCS}$ and $\mathrm{OCHSOH}$ ($\mathrm{O}=\mathrm{CHSOH}$) represent major products of our modelling of marine conditions. These products stem from the $\mathrm{OH}$-initiated oxidation of HPMTF, a pathway that mostly includes structure–activity relationships and theoretical calculations from the paper. As our mechanism is sensitive to these reactions in marine conditions, further experimental studies should be performed to constrain them.

Photolysis reactions remain a major source of uncertainty due to the lack of experimental data. The photolysis rate constants in our mechanism are based on structure–activity relationships of compounds that do not contain sulfur, with the exception of the ${\mathrm{CH}}_{\mathrm{3}}\mathrm{SCHO}$ absorption cross-section used for the photolysis of the carbonyl group. Although the experiments included UV lamps with wavelengths between 300–400 nm, which allowed some evaluation of the photolysis reactions in our mechanism, further experiments exploring the photolysis of DMS oxidation products should be performed.

This paper also highlights the importance of intercomparison studies. By evaluating a mechanism across experiments that include a range of conditions, it reduces the importance of systematic uncertainties in the experiments and ensures the mechanism is robust over a wider range of conditions. Future experiments in different, marine conditions (including reactions with halogens) measuring a wide range of products would be useful to further constrain the DMS mechanism. To increase the ease of modelling these experiments in future intercomparison studies, these studies should include the model input files (representing the experimental parameters) in the Supplement, along with the mechanism files.

Appendix A The mechanism developed in this work

The mechanism developed in this study was based on a thorough literature review to update and improve the DMS mechanism in the MCM. To determine which rate constants should be used in the model, we used the same construction methodology as the MCM . The full mechanism is given in Table A1, with a description of the development of this mechanism included here.

In this methodology, evaluated experimental data took priority. The NASA Panel report and IUPAC provide these evaluations on the experimental data; however, the NASA Panel report provides a more recent review, and as such was relied on during this study. In this mechanism, 9 of the MCM reactions were updated to the NASA Panel report current recommendations, and 13 reactions were added from the report.

When evaluated experimental data were not provided, published experimental data were used. An additional nine reactions were either adjusted or added with experimental values for rate constants. In the case where there were multiple experiments that recorded a rate constant (three reactions in this mechanism), the experimental values were either averaged or evaluated to find the rate constant to use. An average that was not used was the ${\mathrm{CH}}_{\mathrm{3}}{\mathrm{SCH}}_{\mathrm{2}}{\mathrm{O}}_{\mathrm{2}}$ $\mathrm{H}$ shift (for the HPMTF pathway). In that case, the rate constant from was used due to it being a more direct measurement and temperature-dependent (this decision is discussed in more detail in Sect. 5.2).

When there were no experimental data to base rate constants on, structure–activity relationships (SARs) or estimates were used. The MCM provides literature-based SARs related to carbon-based chemistry; however, these SARs do not take into account sulfur chemistry. A comprehensive DMS mechanism paper by considers sulfur rate experiments, ab initio calculations, and bond dissociation energies to estimate rate constants that had not been experimentally determined. However, the MCM is based on more recent experiments and includes temperature dependence. To decide between the estimates of and the SARs of the MCM, we used the following methodology:

• If MCM and used a similar rate constant at 298 K, the MCM value was used as it includes temperature dependence.

• If MCM and used different rate constants, then two outcomes are possible.

  • If has a sulfur-based reasoning for their rate constant, their value was used.

  • Otherwise, the MCM value was used, as it is based on more recent literature.

If there were no appropriate SAR or estimates to use, a theoretical rate constant was used. This was the case for seven reactions. The only theoretical rate constant used that was not for the HPMTF pathway was for the methanesulfinic acid (MSIA) and ozone (${\mathrm{O}}_{\mathrm{3}}$) reaction from . The other rate constants where theoretical studies were applied were the second $\mathrm{H}$ shift forming HPMTF and the reactions following the formation of ${\mathrm{HOOCH}}_{\mathrm{2}}\mathrm{S}$ from the reaction between HPMTF and $\mathrm{OH}$ radicals. Only the major reactions following the HPMTF and $\mathrm{OH}$ reaction based on the theoretical rate constants calculated by were included in the mechanism. For the second $\mathrm{H}$ shift, we had to choose between three theoretical papers: , , and . note that their rate constants for both $\mathrm{H}$ shifts leading to HPMTF are smaller than (by factors of 51 and 317 at 298 K). say that this difference is mainly attributed to the different computational methods and consider their calculations to be more accurate. In addition, the rate constant they calculate for the first $\mathrm{H}$ shift (0.058 s${}^{-\mathrm{1}}$ at 298 K) agrees well with direct measurements from , which is why the rate constant was used in this work for the second $\mathrm{H}$ shift instead of the calculations by or .

In addition to the above, some reactions were adjusted in this work to better fit the chamber experiments. This was done as these rate constants or branching ratios had not been experimentally determined; however, the model was sensitive to them. found theoretically that carbonyl sulfide ($\mathrm{OCS}$) can form from the decomposition of ${\mathrm{HOOCH}}_{\mathrm{2}}\mathrm{SCO}$ (formed from the reaction of HPMTF with $\mathrm{OH}$). However, their predicted branching ratio of direct $\mathrm{OCS}$ formation was 3 % at 298 K, albeit with an uncertainty factor of at least 3. Although some $\mathrm{OCS}$ was formed in our mechanism from the reactions of ${\mathrm{HOOCH}}_{\mathrm{2}}\mathrm{S}$, the branching ratio of $\mathrm{OCS}$ formed from the initial HPMTF reaction with $\mathrm{OH}$ was increased until the modelled $\mathrm{OCS}$ matched the $\mathrm{OCS}$ observed in the experiment. This branching ratio was found to be 9 %, which is within the upper limit of uncertainty calculated by .

Table A1

The sulfur reactions used in the mechanism developed in this work along with their source (if not from the MCM) and reference.

	Reaction			Rate constant	Comment	Ref.
1	${\mathrm{SO}}_{\mathrm{2}}$ + $\mathrm{O}$	=	${\mathrm{SO}}_{\mathrm{3}}$	(KinfT $\times$ K0T$\times$ [M]) $\times \left({\mathrm{0.6}}^{\text{Kb}}\right)/$	Evaluated data	${}^{\mathrm{a}}$
				$(\text{KinfT}+\text{K0T}\times \text{[M]})$
2	${\mathrm{SO}}_{\mathrm{2}}$ + $\mathrm{OH}$	=	${\mathrm{HSO}}_{\mathrm{3}}$	(KinfT2 $\times$ K0T2$\times$ [M])$\times \left({\mathrm{0.6}}^{\text{Kb2}}\right)/$	Evaluated data	${}^{\mathrm{a}}$
				$(\text{KinfT2}+\text{K0T2}\times \text{[M]})$
3	${\mathrm{HSO}}_{\mathrm{3}}$	=	${\mathrm{HO}}_{\mathrm{2}}$ + ${\mathrm{SO}}_{\mathrm{3}}$	$\mathrm{1.3}\times {\mathrm{10}}^{-\mathrm{12}}\times e^{-\mathrm{330}/\text{T}}\times$ [${\mathrm{O}}_{\mathrm{2}}$]	MCM	${}^{\mathrm{b}}$
4	${\mathrm{SO}}_{\mathrm{3}}$	=	SA	$\mathrm{8.5}\times {\mathrm{10}}^{-\mathrm{41}}\times e^{\mathrm{6540}/\text{T}}\times$[${\mathrm{H}}_{\mathrm{2}}\mathrm{O}$]${}^{\mathrm{2}}$	Evaluated data	${}^{\mathrm{a}}$
5	${\mathrm{SO}}_{\mathrm{3}}$ + ${\mathrm{NO}}_{\mathrm{2}}$	=	${\mathrm{NSO}}_{\mathrm{5}}$	$\mathrm{1.00}\times {\mathrm{10}}^{-\mathrm{19}}$	Evaluated data	${}^{\mathrm{a}}$
6	DMS + ${\mathrm{NO}}_{\mathrm{3}}$	=	${\mathrm{CH}}_{\mathrm{3}}{\mathrm{SCH}}_{\mathrm{2}}{\mathrm{O}}_{\mathrm{2}}$ + ${\mathrm{HNO}}_{\mathrm{3}}$	$\mathrm{1.9}\times {\mathrm{10}}^{-\mathrm{13}}\times e^{\mathrm{530}/\text{T}}$	Evaluated data	${}^{\mathrm{a}}$
7	DMS + $\mathrm{OH}$	=	${\mathrm{CH}}_{\mathrm{3}}{\mathrm{SCH}}_{\mathrm{2}}{\mathrm{O}}_{\mathrm{2}}$	$\mathrm{1.2}\times {\mathrm{10}}^{-\mathrm{11}}\times e^{-\mathrm{280}/\text{T}}$	Evaluated data	${}^{\mathrm{a}}$
8	DMS + $\mathrm{O}$	=	${\mathrm{CH}}_{\mathrm{3}}\mathrm{SO}$ + ${\mathrm{CH}}_{\mathrm{3}}{\mathrm{O}}_{\mathrm{2}}$	$\mathrm{1.3}\times {\mathrm{10}}^{-\mathrm{11}}\times e^{\mathrm{410}/\text{T}}$	Evaluated data	${}^{\mathrm{a}}$
9	DMS + $\mathrm{OH}$	=	${\mathrm{CH}}_{\mathrm{3}}\mathrm{S}\left(\mathrm{OH}\right){\mathrm{CH}}_{\mathrm{3}}$	[M]$\times \mathrm{3}\times {\mathrm{10}}^{-\mathrm{31}}\times (\text{T}/\mathrm{298}{)}^{-\mathrm{6.24}}$	Evaluated data	${}^{\mathrm{a}}$
10	${\mathrm{CH}}_{\mathrm{3}}\mathrm{S}\left(\mathrm{OH}\right){\mathrm{CH}}_{\mathrm{3}}$	=	HODMSO2	$\mathrm{8.5}\times {\mathrm{10}}^{-\mathrm{13}}\times$[${\mathrm{O}}_{\mathrm{2}}$]	Evaluated data	${}^{\mathrm{a}}$
11	${\mathrm{CH}}_{\mathrm{3}}\mathrm{S}\left(\mathrm{OH}\right){\mathrm{CH}}_{\mathrm{3}}$	=	${\mathrm{CH}}_{\mathrm{3}}\mathrm{SOH}$ + ${\mathrm{CH}}_{\mathrm{3}}{\mathrm{O}}_{\mathrm{2}}$	$\mathrm{5.0}\times {\mathrm{10}}^{\mathrm{5}}$	Estimate	${}^{\mathrm{c}}$
12	${\mathrm{CH}}_{\mathrm{3}}\mathrm{S}\left(\mathrm{OH}\right){\mathrm{CH}}_{\mathrm{3}}$	=	DMS + $\mathrm{OH}$	$\text{[M]}\times \mathrm{3}\times {\mathrm{10}}^{-\mathrm{31}}\times (\text{T}/\mathrm{298}{)}^{-\mathrm{6.24}}/$	Evaluated data	${}^{\mathrm{a}}$
				$(\mathrm{9.6}\times {\mathrm{10}}^{-\mathrm{27}}\times e^{\mathrm{5376}/\text{T}})$
13	${\mathrm{CH}}_{\mathrm{3}}\mathrm{SOH}$ + $\mathrm{OH}$	=	${\mathrm{CH}}_{\mathrm{3}}\mathrm{SO}$	$\mathrm{5.00}\times {\mathrm{10}}^{-\mathrm{11}}$	Estimate	${}^{\mathrm{d}}$
14	${\mathrm{CH}}_{\mathrm{3}}\mathrm{SOH}$ + ${\mathrm{O}}_{\mathrm{3}}$	=	${\mathrm{CH}}_{\mathrm{3}}{\mathrm{O}}_{\mathrm{2}}$ + ${\mathrm{HO}}_{\mathrm{2}}$ + ${\mathrm{SO}}_{\mathrm{2}}$	$\mathrm{2.00}\times {\mathrm{10}}^{-\mathrm{12}}$	Fit to data	${}^{\mathrm{e}}$
15	${\mathrm{CH}}_{\mathrm{3}}{\mathrm{SCH}}_{\mathrm{2}}{\mathrm{O}}_{\mathrm{2}}$ + ${\mathrm{HO}}_{\mathrm{2}}$	=	${\mathrm{CH}}_{\mathrm{3}}{\mathrm{SCH}}_{\mathrm{2}}\mathrm{OOH}$	KRO2HO2 $\times \mathrm{0.387}$	MCM	${}^{\mathrm{b}}$
16	${\mathrm{CH}}_{\mathrm{3}}{\mathrm{SCH}}_{\mathrm{2}}{\mathrm{O}}_{\mathrm{2}}$ + $\mathrm{NO}$	=	${\mathrm{CH}}_{\mathrm{3}}{\mathrm{SCH}}_{\mathrm{2}}\mathrm{O}$ + ${\mathrm{NO}}_{\mathrm{2}}$	$\mathrm{4.9}\times {\mathrm{10}}^{-\mathrm{12}}\times e^{\mathrm{260}/\text{T}}$	MCM	${}^{\mathrm{b}}$
17	${\mathrm{CH}}_{\mathrm{3}}{\mathrm{SCH}}_{\mathrm{2}}{\mathrm{O}}_{\mathrm{2}}$ + ${\mathrm{NO}}_{\mathrm{2}}$	=	${\mathrm{CH}}_{\mathrm{3}}{\mathrm{SCH}}_{\mathrm{2}}{\mathrm{OONO}}_{\mathrm{2}}$	$\mathrm{9.20}\times {\mathrm{10}}^{-\mathrm{12}}$	Experiment	${}^{\mathrm{f}}$
18	${\mathrm{CH}}_{\mathrm{3}}{\mathrm{SCH}}_{\mathrm{2}}{\mathrm{O}}_{\mathrm{2}}$ + ${\mathrm{NO}}_{\mathrm{3}}$	=	${\mathrm{CH}}_{\mathrm{3}}{\mathrm{SCH}}_{\mathrm{2}}\mathrm{O}$ + ${\mathrm{NO}}_{\mathrm{2}}$	KRO2NO3	MCM	${}^{\mathrm{b}}$
19	${\mathrm{CH}}_{\mathrm{3}}{\mathrm{SCH}}_{\mathrm{2}}{\mathrm{O}}_{\mathrm{2}}$	=	${\mathrm{CH}}_{\mathrm{3}}{\mathrm{SCH}}_{\mathrm{2}}\mathrm{O}$	$\mathrm{3.74}\times {\mathrm{10}}^{-\mathrm{12}}\times$[${\mathrm{RO}}_{\mathrm{2}}$]$\times \mathrm{0.8}$	MCM	${}^{\mathrm{b}}$
20	${\mathrm{CH}}_{\mathrm{3}}{\mathrm{SCH}}_{\mathrm{2}}{\mathrm{O}}_{\mathrm{2}}$	=	${\mathrm{CH}}_{\mathrm{3}}{\mathrm{SCH}}_{\mathrm{2}}\mathrm{OH}$	$\mathrm{3.74}\times {\mathrm{10}}^{-\mathrm{12}}\times$ [${\mathrm{RO}}_{\mathrm{2}}$]$\times \mathrm{0.1}$	MCM	${}^{\mathrm{b}}$
21	${\mathrm{CH}}_{\mathrm{3}}{\mathrm{SCH}}_{\mathrm{2}}{\mathrm{O}}_{\mathrm{2}}$	=	${\mathrm{CH}}_{\mathrm{3}}\mathrm{SCHO}$	$\mathrm{3.74}\times {\mathrm{10}}^{-\mathrm{12}}\times$ [${\mathrm{RO}}_{\mathrm{2}}$] $\times \mathrm{0.1}$	MCM	${}^{\mathrm{b}}$
22	${\mathrm{CH}}_{\mathrm{3}}{\mathrm{SCH}}_{\mathrm{2}}{\mathrm{O}}_{\mathrm{2}}$	=	${\mathrm{HOOCH}}_{\mathrm{2}}{\mathrm{SCH}}_{\mathrm{2}}{\mathrm{O}}_{\mathrm{2}}$	$\mathrm{2.39}\times {\mathrm{10}}^{\mathrm{9}}\times e^{-\mathrm{7278}/\text{T}}$	Experiment	${}^{\mathrm{g}}$
23	${\mathrm{CH}}_{\mathrm{3}}{\mathrm{SCH}}_{\mathrm{2}}{\mathrm{OONO}}_{\mathrm{2}}$	=	${\mathrm{CH}}_{\mathrm{3}}{\mathrm{SCH}}_{\mathrm{2}}{\mathrm{O}}_{\mathrm{2}}$ + ${\mathrm{NO}}_{\mathrm{2}}$	$\mathrm{0.134}$	SAR (${\mathrm{CH}}_{\mathrm{3}}{\mathrm{SO}}_{\mathrm{4}}{\mathrm{NO}}_{\mathrm{2}}$)	${}^{\mathrm{h}}$
24	HODMSO2 + $\mathrm{NO}$	=	DMSO2 + ${\mathrm{HO}}_{\mathrm{2}}$ + ${\mathrm{NO}}_{\mathrm{2}}$	$\mathrm{5.0}\times {\mathrm{10}}^{-\mathrm{12}}$	SAR	${}^{\mathrm{c}}$
25	HODMSO2	=	DMSO + ${\mathrm{HO}}_{\mathrm{2}}$	$\mathrm{8.90}\times {\mathrm{10}}^{\mathrm{10}}\times e^{-\mathrm{6040}/\text{T}}$	MCM	${}^{\mathrm{b}}$
26	${\mathrm{CH}}_{\mathrm{3}}{\mathrm{SCH}}_{\mathrm{2}}\mathrm{OOH}$ + $\mathrm{OH}$	=	${\mathrm{CH}}_{\mathrm{3}}\mathrm{SCHO}$ + $\mathrm{OH}$	$\mathrm{7.03}\times {\mathrm{10}}^{-\mathrm{11}}$	MCM	${}^{\mathrm{b}}$
27	${\mathrm{CH}}_{\mathrm{3}}{\mathrm{SCH}}_{\mathrm{2}}\mathrm{OOH}$	=	${\mathrm{CH}}_{\mathrm{3}}{\mathrm{SCH}}_{\mathrm{2}}\mathrm{O}$ + $\mathrm{OH}$	J(“CH3OOH”)	MCM	${}^{\mathrm{b}}$
28	${\mathrm{CH}}_{\mathrm{3}}{\mathrm{SCH}}_{\mathrm{2}}\mathrm{O}$	=	${\mathrm{CH}}_{\mathrm{3}}\mathrm{S}$ + $\mathrm{HCHO}$	KDEC	MCM	${}^{\mathrm{b}}$
29	${\mathrm{CH}}_{\mathrm{3}}{\mathrm{SCH}}_{\mathrm{2}}\mathrm{OH}$ + $\mathrm{OH}$	=	${\mathrm{CH}}_{\mathrm{3}}\mathrm{SCHO}$ + ${\mathrm{HO}}_{\mathrm{2}}$	$\mathrm{2.78}\times {\mathrm{10}}^{-\mathrm{11}}$	MCM	${}^{\mathrm{b}}$
30	${\mathrm{CH}}_{\mathrm{3}}\mathrm{SCHO}$	=	${\mathrm{CH}}_{\mathrm{3}}\mathrm{S}$ + $\mathrm{CO}$ + ${\mathrm{HO}}_{\mathrm{2}}$	J(“MTF”)	Experiment/SAR	${}^{\mathrm{b},\mathrm{i}}$
31	${\mathrm{CH}}_{\mathrm{3}}\mathrm{SCHO}$ + $\mathrm{OH}$	=	${\mathrm{CH}}_{\mathrm{3}}\mathrm{S}$ + $\mathrm{CO}$	$\mathrm{1.23}\times {\mathrm{10}}^{-\mathrm{11}}$	Experiment	${}^{\mathrm{i},\mathrm{j}}$
32	DMSO2 + $\mathrm{OH}$	=	DMSO2O2	$\mathrm{1.0}\times {\mathrm{10}}^{-\mathrm{14}}$	Estimate	${}^{\mathrm{c}}$
33	DMSO + $\mathrm{OH}$	=	MSIA + ${\mathrm{CH}}_{\mathrm{3}}{\mathrm{O}}_{\mathrm{2}}$	$\mathrm{6.1}\times {\mathrm{10}}^{-\mathrm{12}}\times e^{\mathrm{800}/\text{T}}\times \mathrm{0.9}$	MCM	${}^{\mathrm{b},\mathrm{k}}$
34	DMSO + $\mathrm{OH}$	=	DMSO2 + ${\mathrm{HO}}_{\mathrm{2}}$	$\mathrm{6.1}\times {\mathrm{10}}^{-\mathrm{12}}\times e^{\mathrm{800}/\text{T}}\times \mathrm{0.1}$	Estimate	${}^{\mathrm{k}}$
35	DMSO + $\mathrm{O}$	=	${\mathrm{SO}}_{\mathrm{2}}$ + ${\mathrm{CH}}_{\mathrm{3}}{\mathrm{O}}_{\mathrm{2}}$ + ${\mathrm{CH}}_{\mathrm{3}}{\mathrm{O}}_{\mathrm{2}}$	$\mathrm{2.0}\times {\mathrm{10}}^{-\mathrm{12}}\times e^{\mathrm{440}/\text{T}}$	Evaluated data	${}^{\mathrm{a}}$
36	DMSO + ${\mathrm{NO}}_{\mathrm{3}}$	=	DMSO2 + ${\mathrm{NO}}_{\mathrm{2}}$	$\mathrm{2.90}\times {\mathrm{10}}^{-\mathrm{13}}$	Evaluated data	${}^{\mathrm{a}}$
37	${\mathrm{CH}}_{\mathrm{3}}\mathrm{S}$ + ${\mathrm{NO}}_{\mathrm{2}}$	=	${\mathrm{CH}}_{\mathrm{3}}\mathrm{SO}$ + $\mathrm{NO}$	$\mathrm{3}\times {\mathrm{10}}^{-\mathrm{11}}\times e^{\mathrm{240}/\text{T}}$	Evaluated data	${}^{\mathrm{a}}$
38	${\mathrm{CH}}_{\mathrm{3}}\mathrm{S}$ + ${\mathrm{O}}_{\mathrm{3}}$	=	${\mathrm{CH}}_{\mathrm{3}}\mathrm{SO}$	$\mathrm{1.5}\times {\mathrm{10}}^{-\mathrm{12}}\times e^{\mathrm{360}/\text{T}}$	Evaluated data	${}^{\mathrm{a}}$
39	${\mathrm{CH}}_{\mathrm{3}}\mathrm{S}$	=	${\mathrm{CH}}_{\mathrm{3}}\mathrm{SOO}$	$\mathrm{1.20}\times {\mathrm{10}}^{-\mathrm{16}}\times e^{\mathrm{1580}/\text{T}}\times$ [${\mathrm{O}}_{\mathrm{2}}$]	MCM	${}^{\mathrm{b}}$
40	DMSO2O2 + ${\mathrm{HO}}_{\mathrm{2}}$	=	DMSO2OOH	KRO2HO2 $\times \mathrm{0.387}$	MCM	${}^{\mathrm{b}}$
41	DMSO2O2 + $\mathrm{NO}$	=	DMSO2O + ${\mathrm{NO}}_{\mathrm{2}}$	KRO2NO	MCM	${}^{\mathrm{b}}$
42	DMSO2O2 + ${\mathrm{NO}}_{\mathrm{3}}$	=	DMSO2O + ${\mathrm{NO}}_{\mathrm{2}}$	KRO2NO3	MCM	${}^{\mathrm{b}}$
43	DMSO2O2	=	${\mathrm{CH}}_{\mathrm{3}}{\mathrm{SO}}_{\mathrm{2}}\mathrm{CHO}$	$\mathrm{2.00}\times {\mathrm{10}}^{-\mathrm{12}}\times$ [${\mathrm{RO}}_{\mathrm{2}}$] $\times \mathrm{0.2}$	MCM	${}^{\mathrm{b}}$
44	DMSO2O2	=	DMSO2O	$\mathrm{2.00}\times {\mathrm{10}}^{-\mathrm{12}}\times$ [${\mathrm{RO}}_{\mathrm{2}}$] $\times \mathrm{0.6}$	MCM	${}^{\mathrm{b}}$
45	DMSO2O2	=	DMSO2OH	$\mathrm{2.00}\times {\mathrm{10}}^{-\mathrm{12}}\times$ [${\mathrm{RO}}_{\mathrm{2}}$] $\times \mathrm{0.2}$	MCM	${}^{\mathrm{b}}$
46	MSIA + $\mathrm{OH}$	=	${\mathrm{CH}}_{\mathrm{3}}{\mathrm{SO}}_{\mathrm{2}}$	$\mathrm{9.00}\times {\mathrm{10}}^{-\mathrm{11}}$	Evaluated data	${}^{\mathrm{a}}$
47	MSIA + ${\mathrm{NO}}_{\mathrm{3}}$	=	${\mathrm{CH}}_{\mathrm{3}}{\mathrm{SO}}_{\mathrm{2}}$ + ${\mathrm{HNO}}_{\mathrm{3}}$	$\mathrm{1.0}\times {\mathrm{10}}^{-\mathrm{13}}$	Estimate	${}^{\mathrm{c}}$
48	MSIA + ${\mathrm{O}}_{\mathrm{3}}$	=	MSA	$\mathrm{1.79}\times {\mathrm{10}}^{-\mathrm{22}}$	Theory	${}^{\mathrm{l}}$
49	${\mathrm{CH}}_{\mathrm{3}}\mathrm{SO}$ + ${\mathrm{NO}}_{\mathrm{2}}$	=	${\mathrm{CH}}_{\mathrm{3}}{\mathrm{O}}_{\mathrm{2}}$ + ${\mathrm{SO}}_{\mathrm{2}}$ + $\mathrm{NO}$	$\mathrm{1.20}\times {\mathrm{10}}^{-\mathrm{11}}\times \mathrm{0.25}$	MCM	${}^{\mathrm{b}}$
50	${\mathrm{CH}}_{\mathrm{3}}\mathrm{SO}$ + ${\mathrm{NO}}_{\mathrm{2}}$	=	${\mathrm{CH}}_{\mathrm{3}}{\mathrm{SO}}_{\mathrm{2}}$ + $\mathrm{NO}$	$\mathrm{1.20}\times {\mathrm{10}}^{-\mathrm{11}}\times \mathrm{0.75}$	MCM	${}^{\mathrm{b}}$

Table A1

Continued.

	Reaction			Rate constant	Comment	Ref.
51	${\mathrm{CH}}_{\mathrm{3}}\mathrm{SO}$ + ${\mathrm{O}}_{\mathrm{3}}$	=	${\mathrm{CH}}_{\mathrm{3}}{\mathrm{O}}_{\mathrm{2}}$ + ${\mathrm{SO}}_{\mathrm{2}}$	$\mathrm{4.00}\times {\mathrm{10}}^{-\mathrm{13}}$	MCM	${}^{\mathrm{b}}$
52	${\mathrm{CH}}_{\mathrm{3}}\mathrm{SO}$	=	${\mathrm{CH}}_{\mathrm{3}}{\mathrm{SOO}}_{\mathrm{2}}$	$\mathrm{3.12}\times {\mathrm{10}}^{-\mathrm{16}}\times e^{\mathrm{1580}/\text{T}}\times$ [${\mathrm{O}}_{\mathrm{2}}$]	MCM	${}^{\mathrm{b}}$
53	${\mathrm{CH}}_{\mathrm{3}}\mathrm{SOO}$ + $\mathrm{NO}$	=	${\mathrm{CH}}_{\mathrm{3}}\mathrm{SO}$ + ${\mathrm{NO}}_{\mathrm{2}}$	$\mathrm{1.10}\times {\mathrm{10}}^{-\mathrm{11}}$	MCM	${}^{\mathrm{b}}$
54	${\mathrm{CH}}_{\mathrm{3}}\mathrm{SOO}$ + ${\mathrm{NO}}_{\mathrm{2}}$	=	${\mathrm{CH}}_{\mathrm{3}}\mathrm{SO}$ + ${\mathrm{NO}}_{\mathrm{3}}$	$\mathrm{2.20}\times {\mathrm{10}}^{-\mathrm{11}}$	MCM	${}^{\mathrm{b}}$
55	${\mathrm{CH}}_{\mathrm{3}}\mathrm{SOO}$	=	${\mathrm{CH}}_{\mathrm{3}}{\mathrm{O}}_{\mathrm{2}}$ + ${\mathrm{SO}}_{\mathrm{2}}$	$\mathrm{5}$	Fit to data	${}^{\mathrm{m}}$
56	${\mathrm{CH}}_{\mathrm{3}}\mathrm{SOO}$	=	${\mathrm{CH}}_{\mathrm{3}}\mathrm{S}$	$\mathrm{1.5}\times {\mathrm{10}}^{\mathrm{5}}$	Fit to data	${}^{\mathrm{m}}$
57	${\mathrm{CH}}_{\mathrm{3}}\mathrm{SOO}$ + ${\mathrm{HO}}_{\mathrm{2}}$	=	${\mathrm{CH}}_{\mathrm{3}}\mathrm{SOOH}$	$\mathrm{4.0}\times {\mathrm{10}}^{-\mathrm{12}}$	SAR	${}^{\mathrm{c}}$
58	${\mathrm{CH}}_{\mathrm{3}}\mathrm{SOOH}$ + $\mathrm{OH}$	=	${\mathrm{CH}}_{\mathrm{3}}\mathrm{SOO}$	$\mathrm{3.68}\times {\mathrm{10}}^{-\mathrm{13}}\times e^{\mathrm{635}/\text{T}}$	SAR ($\mathrm{ROOH}$)	${}^{\mathrm{n}}$
59	${\mathrm{CH}}_{\mathrm{3}}\mathrm{SOOH}$	=	${\mathrm{CH}}_{\mathrm{3}}\mathrm{SO}$ + $\mathrm{OH}$	J(“CH3OOH”)	SAR (${\mathrm{CH}}_{\mathrm{3}}\mathrm{OOH}$)	${}^{\mathrm{b}}$
60	DMSO2OOH + $\mathrm{OH}$	=	${\mathrm{CH}}_{\mathrm{3}}{\mathrm{SO}}_{\mathrm{2}}\mathrm{CHO}$ + $\mathrm{OH}$	$\mathrm{1.26}\times {\mathrm{10}}^{-\mathrm{12}}$	MCM	${}^{\mathrm{b}}$
61	DMSO2OOH + $\mathrm{OH}$	=	DMSO2O2	$\mathrm{3.60}\times {\mathrm{10}}^{-\mathrm{12}}$	MCM	${}^{\mathrm{b}}$
62	DMSO2OOH	=	DMSO2O + $\mathrm{OH}$	J(“CH3OOH”)	MCM	${}^{\mathrm{b}}$
63	DMSO2O	=	${\mathrm{CH}}_{\mathrm{3}}{\mathrm{SO}}_{\mathrm{2}}$ + $\mathrm{HCHO}$	KDEC	MCM	${}^{\mathrm{b}}$
64	${\mathrm{CH}}_{\mathrm{3}}{\mathrm{SO}}_{\mathrm{2}}\mathrm{CHO}$ + $\mathrm{OH}$	=	${\mathrm{CH}}_{\mathrm{3}}{\mathrm{SO}}_{\mathrm{2}}$ + $\mathrm{CO}$	$\mathrm{1.78}\times {\mathrm{10}}^{-\mathrm{12}}$	MCM	${}^{\mathrm{b}}$
65	${\mathrm{CH}}_{\mathrm{3}}{\mathrm{SO}}_{\mathrm{2}}\mathrm{CHO}$	=	${\mathrm{CH}}_{\mathrm{3}}{\mathrm{SO}}_{\mathrm{2}}$ + $\mathrm{CO}$ + ${\mathrm{HO}}_{\mathrm{2}}$	J(“MTF”)	SAR (MTF)	${}^{\mathrm{b},\mathrm{i}}$
66	DMSO2OH + $\mathrm{OH}$	=	${\mathrm{CH}}_{\mathrm{3}}{\mathrm{SO}}_{\mathrm{2}}\mathrm{CHO}$ + ${\mathrm{HO}}_{\mathrm{2}}$	$\mathrm{5.23}\times {\mathrm{10}}^{-\mathrm{13}}$	MCM	${}^{\mathrm{b}}$
67	DMSO2OH + $\mathrm{OH}$	=	DMSO2O	$\mathrm{1.40}\times {\mathrm{10}}^{-\mathrm{13}}$	MCM	${}^{\mathrm{b}}$
68	${\mathrm{CH}}_{\mathrm{3}}{\mathrm{SO}}_{\mathrm{2}}$ + ${\mathrm{O}}_{\mathrm{3}}$	=	${\mathrm{CH}}_{\mathrm{3}}{\mathrm{SO}}_{\mathrm{3}}$	$\mathrm{3.00}\times {\mathrm{10}}^{-\mathrm{13}}$	MCM	${}^{\mathrm{b}}$
69	${\mathrm{CH}}_{\mathrm{3}}{\mathrm{SO}}_{\mathrm{2}}$	=	${\mathrm{CH}}_{\mathrm{3}}{\mathrm{O}}_{\mathrm{2}}$ + ${\mathrm{SO}}_{\mathrm{2}}$	$\mathrm{1.04}\times {\mathrm{10}}^{\mathrm{15}}\times e^{-\mathrm{9673}/\text{T}}$	Fit to data	${}^{\mathrm{b},\mathrm{h}}$
70	${\mathrm{CH}}_{\mathrm{3}}{\mathrm{SO}}_{\mathrm{2}}$	=	${\mathrm{CH}}_{\mathrm{3}}{\mathrm{SO}}_{\mathrm{2}}{\mathrm{O}}_{\mathrm{2}}$	$\mathrm{1.03}\times {\mathrm{10}}^{-\mathrm{16}}\times e^{\mathrm{1580}/\text{T}}\times$ [${\mathrm{O}}_{\mathrm{2}}$]	MCM	${}^{\mathrm{b}}$
71	${\mathrm{CH}}_{\mathrm{3}}{\mathrm{SO}}_{\mathrm{2}}$ + ${\mathrm{NO}}_{\mathrm{2}}$	=	${\mathrm{CH}}_{\mathrm{3}}{\mathrm{SO}}_{\mathrm{3}}$ + $\mathrm{NO}$	$\mathrm{2.20}\times {\mathrm{10}}^{-\mathrm{12}}$	Evaluated data	${}^{\mathrm{a}}$
72	${\mathrm{CH}}_{\mathrm{3}}{\mathrm{SO}}_{\mathrm{2}}$ + $\mathrm{OH}$	=	MSA	$\mathrm{5.0}\times {\mathrm{10}}^{-\mathrm{11}}$	Estimate	${}^{\mathrm{c}}$
73	${\mathrm{CH}}_{\mathrm{3}}{\mathrm{SOO}}_{\mathrm{2}}$ + ${\mathrm{HO}}_{\mathrm{2}}$	=	${\mathrm{CH}}_{\mathrm{3}}{\mathrm{SO}}_{\mathrm{2}}$ + $\mathrm{OH}$	KAPHO2 $\times \mathrm{0.44}$	MCM	${}^{\mathrm{b}}$
74	${\mathrm{CH}}_{\mathrm{3}}{\mathrm{SOO}}_{\mathrm{2}}$ + ${\mathrm{HO}}_{\mathrm{2}}$	=	${\mathrm{CH}}_{\mathrm{3}}\mathrm{SOOOH}$	KAPHO2 $\times \mathrm{0.41}$	MCM	${}^{\mathrm{b}}$
75	${\mathrm{CH}}_{\mathrm{3}}{\mathrm{SOO}}_{\mathrm{2}}$ + ${\mathrm{HO}}_{\mathrm{2}}$	=	MSIA + ${\mathrm{O}}_{\mathrm{3}}$	KAPHO2 $\times \mathrm{0.15}$	MCM	${}^{\mathrm{b}}$
76	${\mathrm{CH}}_{\mathrm{3}}{\mathrm{SOO}}_{\mathrm{2}}$ + $\mathrm{NO}$	=	${\mathrm{CH}}_{\mathrm{3}}{\mathrm{SO}}_{\mathrm{2}}$ + ${\mathrm{NO}}_{\mathrm{2}}$	$\mathrm{1.00}\times {\mathrm{10}}^{-\mathrm{11}}$	MCM	${}^{\mathrm{b}}$
77	${\mathrm{CH}}_{\mathrm{3}}{\mathrm{SOO}}_{\mathrm{2}}$ + ${\mathrm{NO}}_{\mathrm{2}}$	=	${\mathrm{CH}}_{\mathrm{3}}{\mathrm{SOO}}_{\mathrm{2}}{\mathrm{NO}}_{\mathrm{2}}$	$\mathrm{1.20}\times {\mathrm{10}}^{-\mathrm{12}}\times (\text{T}/\mathrm{300}{)}^{-\mathrm{0.9}}$	MCM	${}^{\mathrm{b}}$
78	${\mathrm{CH}}_{\mathrm{3}}{\mathrm{SOO}}_{\mathrm{2}}$ + ${\mathrm{NO}}_{\mathrm{3}}$	=	${\mathrm{CH}}_{\mathrm{3}}{\mathrm{SO}}_{\mathrm{2}}$ + ${\mathrm{NO}}_{\mathrm{2}}$	KRO2NO3 $\times \mathrm{1.74}$	MCM	${}^{\mathrm{b}}$
79	${\mathrm{CH}}_{\mathrm{3}}{\mathrm{SOO}}_{\mathrm{2}}$	=	${\mathrm{CH}}_{\mathrm{3}}\mathrm{SO}$	$\mathrm{9.10}\times {\mathrm{10}}^{\mathrm{10}}\times e^{-\mathrm{3560}/\text{T}}$	MCM	${}^{\mathrm{b}}$
80	${\mathrm{CH}}_{\mathrm{3}}{\mathrm{SOO}}_{\mathrm{2}}$	=	${\mathrm{CH}}_{\mathrm{3}}{\mathrm{SO}}_{\mathrm{2}}$	$\mathrm{1.00}\times {\mathrm{10}}^{-\mathrm{11}}\times$ [${\mathrm{RO}}_{\mathrm{2}}$] $\times \mathrm{0.7}$	MCM	${}^{\mathrm{b}}$
81	${\mathrm{CH}}_{\mathrm{3}}{\mathrm{SOO}}_{\mathrm{2}}$	=	MSIA	$\mathrm{1.00}\times {\mathrm{10}}^{-\mathrm{11}}\times$ [${\mathrm{RO}}_{\mathrm{2}}$] $\times \mathrm{0.3}$	MCM	${}^{\mathrm{b}}$
82	${\mathrm{CH}}_{\mathrm{3}}{\mathrm{SO}}_{\mathrm{3}}$ + ${\mathrm{HO}}_{\mathrm{2}}$	=	MSA	$\mathrm{5.00}\times {\mathrm{10}}^{-\mathrm{11}}$	MCM	${}^{\mathrm{b}}$
83	${\mathrm{CH}}_{\mathrm{3}}{\mathrm{SO}}_{\mathrm{3}}$	=	${\mathrm{CH}}_{\mathrm{3}}{\mathrm{O}}_{\mathrm{2}}$ + ${\mathrm{SO}}_{\mathrm{3}}$	$\mathrm{3.34}\times {\mathrm{10}}^{\mathrm{13}}\times e^{-\mathrm{9946}/\text{T}}$	Experiment	${}^{\mathrm{b},\mathrm{o}}$
84	${\mathrm{CH}}_{\mathrm{3}}{\mathrm{SO}}_{\mathrm{3}}$ + MSIA	=	MSA + ${\mathrm{CH}}_{\mathrm{3}}{\mathrm{SO}}_{\mathrm{2}}$	$\mathrm{1.0}\times {\mathrm{10}}^{-\mathrm{13}}$	Estimate	${}^{\mathrm{c}}$
85	${\mathrm{CH}}_{\mathrm{3}}{\mathrm{SO}}_{\mathrm{3}}$ + $\mathrm{HCHO}$	=	MSA + ${\mathrm{HO}}_{\mathrm{2}}$ + $\mathrm{CO}$	$\mathrm{1.6}\times {\mathrm{10}}^{-\mathrm{15}}$	Estimate	${}^{\mathrm{c}}$
86	${\mathrm{CH}}_{\mathrm{3}}{\mathrm{SO}}_{\mathrm{3}}$ + $\mathrm{HONO}$	=	MSA + ${\mathrm{NO}}_{\mathrm{2}}$	$\mathrm{6.6}\times {\mathrm{10}}^{-\mathrm{16}}$	Estimate	${}^{\mathrm{c}}$
87	${\mathrm{CH}}_{\mathrm{3}}{\mathrm{SO}}_{\mathrm{3}}$ + ${\mathrm{H}}_{\mathrm{2}}{\mathrm{O}}_{\mathrm{2}}$	=	MSA + ${\mathrm{HO}}_{\mathrm{2}}$	$\mathrm{3.0}\times {\mathrm{10}}^{-\mathrm{16}}$	Estimate	${}^{\mathrm{c}}$
88	${\mathrm{CH}}_{\mathrm{3}}{\mathrm{SO}}_{\mathrm{3}}$ + ${\mathrm{CH}}_{\mathrm{3}}\mathrm{OOH}$	=	MSA + ${\mathrm{CH}}_{\mathrm{3}}{\mathrm{O}}_{\mathrm{2}}$	$\mathrm{3.0}\times {\mathrm{10}}^{-\mathrm{16}}$	Estimate	${}^{\mathrm{c}}$
89	${\mathrm{CH}}_{\mathrm{3}}{\mathrm{SO}}_{\mathrm{3}}$ + ${\mathrm{CH}}_{\mathrm{3}}\mathrm{OH}$	=	MSA + ${\mathrm{HO}}_{\mathrm{2}}$ + $\mathrm{HCHO}$	$\mathrm{1.0}\times {\mathrm{10}}^{-\mathrm{16}}$	Estimate	${}^{\mathrm{c}}$
90	${\mathrm{CH}}_{\mathrm{3}}{\mathrm{SO}}_{\mathrm{3}}$ + DMS	=	MSA + ${\mathrm{CH}}_{\mathrm{3}}{\mathrm{SCH}}_{\mathrm{2}}{\mathrm{O}}_{\mathrm{2}}$	$\mathrm{1.45}\times {\mathrm{10}}^{-\mathrm{13}}$	Estimate	${}^{\mathrm{c},\mathrm{h}}$
91	${\mathrm{CH}}_{\mathrm{3}}{\mathrm{SO}}_{\mathrm{2}}{\mathrm{O}}_{\mathrm{2}}$ + ${\mathrm{HO}}_{\mathrm{2}}$	=	${\mathrm{CH}}_{\mathrm{3}}{\mathrm{SO}}_{\mathrm{2}}\mathrm{OOH}$	KAPHO2 $\times \mathrm{0.41}$	MCM	${}^{\mathrm{b}}$
92	${\mathrm{CH}}_{\mathrm{3}}{\mathrm{SO}}_{\mathrm{2}}{\mathrm{O}}_{\mathrm{2}}$ + ${\mathrm{HO}}_{\mathrm{2}}$	=	${\mathrm{CH}}_{\mathrm{3}}{\mathrm{SO}}_{\mathrm{3}}$ + $\mathrm{OH}$	KAPHO2 $\times \mathrm{0.44}$	MCM	${}^{\mathrm{b}}$
93	${\mathrm{CH}}_{\mathrm{3}}{\mathrm{SO}}_{\mathrm{2}}{\mathrm{O}}_{\mathrm{2}}$ + ${\mathrm{HO}}_{\mathrm{2}}$	=	MSA + ${\mathrm{O}}_{\mathrm{3}}$	KAPHO2 $\times \mathrm{0.15}$	MCM	${}^{\mathrm{b}}$
94	${\mathrm{CH}}_{\mathrm{3}}{\mathrm{SO}}_{\mathrm{2}}{\mathrm{O}}_{\mathrm{2}}$ + $\mathrm{NO}$	=	${\mathrm{CH}}_{\mathrm{3}}{\mathrm{SO}}_{\mathrm{3}}$ + ${\mathrm{NO}}_{\mathrm{2}}$	$\mathrm{1.00}\times {\mathrm{10}}^{-\mathrm{11}}$	MCM	${}^{\mathrm{b}}$
95	${\mathrm{CH}}_{\mathrm{3}}{\mathrm{SO}}_{\mathrm{2}}{\mathrm{O}}_{\mathrm{2}}$ + ${\mathrm{NO}}_{\mathrm{2}}$	=	${\mathrm{CH}}_{\mathrm{3}}{\mathrm{SO}}_{\mathrm{4}}{\mathrm{NO}}_{\mathrm{2}}$	$\mathrm{1.20}\times {\mathrm{10}}^{-\mathrm{12}}\times (\text{T}/\mathrm{300}{)}^{-\mathrm{0.9}}$	MCM	${}^{\mathrm{b}}$
96	${\mathrm{CH}}_{\mathrm{3}}{\mathrm{SO}}_{\mathrm{2}}{\mathrm{O}}_{\mathrm{2}}$ + ${\mathrm{NO}}_{\mathrm{3}}$	=	${\mathrm{CH}}_{\mathrm{3}}{\mathrm{SO}}_{\mathrm{3}}$ + ${\mathrm{NO}}_{\mathrm{2}}$	KRO2NO3 $\times \mathrm{1.74}$	MCM	${}^{\mathrm{b}}$
97	${\mathrm{CH}}_{\mathrm{3}}{\mathrm{SO}}_{\mathrm{2}}{\mathrm{O}}_{\mathrm{2}}$	=	${\mathrm{CH}}_{\mathrm{3}}{\mathrm{SO}}_{\mathrm{2}}$	$\mathrm{3.01}\times {\mathrm{10}}^{\mathrm{10}}\times e^{-\mathrm{3560}/\text{T}}$	MCM	${}^{\mathrm{b}}$
98	${\mathrm{CH}}_{\mathrm{3}}{\mathrm{SO}}_{\mathrm{2}}{\mathrm{O}}_{\mathrm{2}}$	=	${\mathrm{CH}}_{\mathrm{3}}{\mathrm{SO}}_{\mathrm{3}}$	$\mathrm{1.00}\times {\mathrm{10}}^{-\mathrm{11}}\times$ [${\mathrm{RO}}_{\mathrm{2}}$] $\times \mathrm{0.7}$	MCM	${}^{\mathrm{b}}$
99	${\mathrm{CH}}_{\mathrm{3}}{\mathrm{SO}}_{\mathrm{2}}{\mathrm{O}}_{\mathrm{2}}$	=	MSA	$\mathrm{1.00}\times {\mathrm{10}}^{-\mathrm{11}}\times$ [${\mathrm{RO}}_{\mathrm{2}}$] $\times \mathrm{0.3}$	MCM	${}^{\mathrm{b}}$
100	${\mathrm{CH}}_{\mathrm{3}}\mathrm{SOOOH}$ + $\mathrm{OH}$	=	${\mathrm{CH}}_{\mathrm{3}}{\mathrm{SOO}}_{\mathrm{2}}$	$\mathrm{9.00}\times {\mathrm{10}}^{-\mathrm{11}}$	MCM	${}^{\mathrm{b}}$
101	${\mathrm{CH}}_{\mathrm{3}}\mathrm{SOOOH}$	=	${\mathrm{CH}}_{\mathrm{3}}{\mathrm{SO}}_{\mathrm{2}}$ + $\mathrm{OH}$	J(“CH3OOH”)	MCM	${}^{\mathrm{b}}$
102	${\mathrm{CH}}_{\mathrm{3}}{\mathrm{SOO}}_{\mathrm{2}}{\mathrm{NO}}_{\mathrm{2}}$ + $\mathrm{OH}$	=	MSIA + ${\mathrm{NO}}_{\mathrm{2}}$	$\mathrm{1.00}\times {\mathrm{10}}^{-\mathrm{11}}$	MCM	${}^{\mathrm{b}}$
103	${\mathrm{CH}}_{\mathrm{3}}{\mathrm{SOO}}_{\mathrm{2}}{\mathrm{NO}}_{\mathrm{2}}$	=	${\mathrm{CH}}_{\mathrm{3}}{\mathrm{SOO}}_{\mathrm{2}}$ + ${\mathrm{NO}}_{\mathrm{2}}$	$\mathrm{0.134}$	SAR (${\mathrm{CH}}_{\mathrm{3}}{\mathrm{SO}}_{\mathrm{4}}{\mathrm{NO}}_{\mathrm{2}}$)	${}^{\mathrm{h}}$

Table A1

Continued.

	Reaction			Rate constant	Comment	Ref.
104	MSA + $\mathrm{OH}$	=	${\mathrm{CH}}_{\mathrm{3}}{\mathrm{SO}}_{\mathrm{3}}$	$\mathrm{2.24}\times {\mathrm{10}}^{-\mathrm{14}}$	MCM	${}^{\mathrm{b}}$
105	${\mathrm{CH}}_{\mathrm{3}}{\mathrm{SO}}_{\mathrm{2}}\mathrm{OOH}$ + $\mathrm{OH}$	=	${\mathrm{CH}}_{\mathrm{3}}{\mathrm{SO}}_{\mathrm{2}}{\mathrm{O}}_{\mathrm{2}}$	$\mathrm{3.60}\times {\mathrm{10}}^{-\mathrm{12}}$	MCM	${}^{\mathrm{b}}$
106	${\mathrm{CH}}_{\mathrm{3}}{\mathrm{SO}}_{\mathrm{2}}\mathrm{OOH}$	=	${\mathrm{CH}}_{\mathrm{3}}{\mathrm{SO}}_{\mathrm{3}}$ + $\mathrm{OH}$	J(“CH3OOH”)	MCM	${}^{\mathrm{b}}$
107	${\mathrm{CH}}_{\mathrm{3}}{\mathrm{SO}}_{\mathrm{4}}{\mathrm{NO}}_{\mathrm{2}}$ + $\mathrm{OH}$	=	${\mathrm{CH}}_{\mathrm{3}}{\mathrm{SO}}_{\mathrm{2}}{\mathrm{O}}_{\mathrm{2}}$ + ${\mathrm{HNO}}_{\mathrm{3}}$	$\mathrm{3.60}\times {\mathrm{10}}^{-\mathrm{13}}$	MCM	${}^{\mathrm{b}}$
108	${\mathrm{CH}}_{\mathrm{3}}{\mathrm{SO}}_{\mathrm{4}}{\mathrm{NO}}_{\mathrm{2}}$	=	${\mathrm{CH}}_{\mathrm{3}}{\mathrm{SO}}_{\mathrm{2}}{\mathrm{O}}_{\mathrm{2}}$ + ${\mathrm{NO}}_{\mathrm{2}}$	$\mathrm{0.134}$	Fit to experiment	${}^{\mathrm{h}}$
109	${\mathrm{HOOCH}}_{\mathrm{2}}{\mathrm{SCH}}_{\mathrm{2}}{\mathrm{O}}_{\mathrm{2}}$	=	HPMTF + $\mathrm{OH}$	$\mathrm{6.1}\times {\mathrm{10}}^{\mathrm{11}}\times e^{-\mathrm{9.5}\times {\mathrm{10}}^{\mathrm{3}}/\text{T}+\mathrm{1.1}\times {\mathrm{10}}^{\mathrm{8}}/{\text{T}}^{\mathrm{3}}}$	Theory	${}^{\mathrm{p}}$
110	${\mathrm{HOOCH}}_{\mathrm{2}}{\mathrm{SCH}}_{\mathrm{2}}{\mathrm{O}}_{\mathrm{2}}$ + $\mathrm{NO}$	=	${\mathrm{HOOCH}}_{\mathrm{2}}{\mathrm{SCH}}_{\mathrm{2}}\mathrm{O}$ + ${\mathrm{NO}}_{\mathrm{2}}$	$\mathrm{4.9}\times {\mathrm{10}}^{-\mathrm{12}}\times e^{\mathrm{260}/\text{T}}$	SAR (${\mathrm{CH}}_{\mathrm{3}}{\mathrm{SCH}}_{\mathrm{2}}{\mathrm{O}}_{\mathrm{2}}$)	${}^{\mathrm{b}}$
111	${\mathrm{HOOCH}}_{\mathrm{2}}{\mathrm{SCH}}_{\mathrm{2}}{\mathrm{O}}_{\mathrm{2}}$ + ${\mathrm{HO}}_{\mathrm{2}}$	=	${\mathrm{HOOCH}}_{\mathrm{2}}{\mathrm{SCH}}_{\mathrm{2}}\mathrm{OOH}$	KRO2HO2 $\times \mathrm{0.387}$	SAR (${\mathrm{CH}}_{\mathrm{3}}{\mathrm{SCH}}_{\mathrm{2}}{\mathrm{O}}_{\mathrm{2}}$)	${}^{\mathrm{b}}$
112	${\mathrm{HOOCH}}_{\mathrm{2}}{\mathrm{SCH}}_{\mathrm{2}}{\mathrm{O}}_{\mathrm{2}}$ + ${\mathrm{NO}}_{\mathrm{3}}$	=	${\mathrm{HOOCH}}_{\mathrm{2}}{\mathrm{SCH}}_{\mathrm{2}}\mathrm{O}$ + ${\mathrm{NO}}_{\mathrm{2}}$	KRO2NO3	SAR (${\mathrm{CH}}_{\mathrm{3}}{\mathrm{SCH}}_{\mathrm{2}}{\mathrm{O}}_{\mathrm{2}}$)	${}^{\mathrm{b}}$
113	${\mathrm{HOOCH}}_{\mathrm{2}}{\mathrm{SCH}}_{\mathrm{2}}{\mathrm{O}}_{\mathrm{2}}$	=	${\mathrm{HOOCH}}_{\mathrm{2}}{\mathrm{SCH}}_{\mathrm{2}}\mathrm{O}$	$\mathrm{3.74}\times {\mathrm{10}}^{-\mathrm{12}}\times$ [${\mathrm{RO}}_{\mathrm{2}}$] $\times \mathrm{0.8}$	SAR (${\mathrm{CH}}_{\mathrm{3}}{\mathrm{SCH}}_{\mathrm{2}}{\mathrm{O}}_{\mathrm{2}}$)	${}^{\mathrm{b}}$
114	${\mathrm{HOOCH}}_{\mathrm{2}}{\mathrm{SCH}}_{\mathrm{2}}{\mathrm{O}}_{\mathrm{2}}$	=	${\mathrm{HOOCH}}_{\mathrm{2}}{\mathrm{SCH}}_{\mathrm{2}}\mathrm{OH}$	$\mathrm{3.74}\times {\mathrm{10}}^{-\mathrm{12}}\times$ [${\mathrm{RO}}_{\mathrm{2}}$] $\times \mathrm{0.1}$	SAR (${\mathrm{CH}}_{\mathrm{3}}{\mathrm{SCH}}_{\mathrm{2}}{\mathrm{O}}_{\mathrm{2}}$)	${}^{\mathrm{b}}$
115	${\mathrm{HOOCH}}_{\mathrm{2}}{\mathrm{SCH}}_{\mathrm{2}}{\mathrm{O}}_{\mathrm{2}}$	=	HPMTF	$\mathrm{3.74}\times {\mathrm{10}}^{-\mathrm{12}}\times$ [${\mathrm{RO}}_{\mathrm{2}}$] $\times \mathrm{0.1}$	SAR (${\mathrm{CH}}_{\mathrm{3}}{\mathrm{SCH}}_{\mathrm{2}}{\mathrm{O}}_{\mathrm{2}}$)	${}^{\mathrm{b}}$
116	${\mathrm{HOOCH}}_{\mathrm{2}}{\mathrm{SCH}}_{\mathrm{2}}\mathrm{O}$	=	${\mathrm{HOOCH}}_{\mathrm{2}}\mathrm{S}$ + $\mathrm{HCHO}$	KDEC	SAR (${\mathrm{CH}}_{\mathrm{3}}{\mathrm{SCH}}_{\mathrm{2}}\mathrm{O}$)	${}^{\mathrm{b}}$
117	HPMTF + $\mathrm{OH}$	=	${\mathrm{HOOCH}}_{\mathrm{2}}\mathrm{S}$ + $\mathrm{CO}$	$\mathrm{1.75}\times {\mathrm{10}}^{-\mathrm{11}}\times \mathrm{0.91}$	Experiment	${}^{\mathrm{h},\mathrm{j},\mathrm{k}}$
118	HPMTF + $\mathrm{OH}$	=	$\mathrm{OH}$ + $\mathrm{HCHO}$ + $\mathrm{OCS}$	$\mathrm{1.75}\times {\mathrm{10}}^{-\mathrm{11}}\times \mathrm{0.09}$	Fit to data	${}^{\mathrm{h},\mathrm{j},\mathrm{k}}$
119	HPMTF	=	${\mathrm{HOOCH}}_{\mathrm{2}}\mathrm{S}$ + ${\mathrm{HO}}_{\mathrm{2}}$ + $\mathrm{CO}$	J(“MTF”)	SAR (MTF)	${}^{\mathrm{b},\mathrm{i}}$
120	HPMTF	=	${\mathrm{OCH}}_{\mathrm{2}}\mathrm{SCHO}$ + $\mathrm{OH}$	J(“CH3OOH”)	SAR (${\mathrm{CH}}_{\mathrm{3}}\mathrm{OOH}$)	${}^{\mathrm{b}}$
121	$\mathrm{OCHSCHO}$ + $\mathrm{OH}$	=	$\mathrm{OCS}$ + $\mathrm{CO}$ + ${\mathrm{HO}}_{\mathrm{2}}$	$\mathrm{2.6}\times {\mathrm{10}}^{-\mathrm{11}}$	Theory	${}^{\mathrm{k}}$
122	$\mathrm{OCHSCHO}$	=	$\mathrm{OCS}$ + ${\mathrm{HO}}_{\mathrm{2}}$ + $\mathrm{CO}$ + ${\mathrm{HO}}_{\mathrm{2}}$	J(“MTF”)	SAR (MTF)	${}^{\mathrm{b},\mathrm{i}}$
123	${\mathrm{HOOCH}}_{\mathrm{2}}{\mathrm{SCH}}_{\mathrm{2}}\mathrm{OH}$ + $\mathrm{OH}$	=	HPMTF + ${\mathrm{HO}}_{\mathrm{2}}$	$\mathrm{2.78}\times {\mathrm{10}}^{-\mathrm{11}}$	SAR (${\mathrm{CH}}_{\mathrm{3}}{\mathrm{SCH}}_{\mathrm{2}}\mathrm{OH}$)	${}^{\mathrm{b}}$
124	${\mathrm{HOOCH}}_{\mathrm{2}}{\mathrm{SCH}}_{\mathrm{2}}\mathrm{OOH}$ + $\mathrm{OH}$	=	${\mathrm{HOOCH}}_{\mathrm{2}}{\mathrm{SCH}}_{\mathrm{2}}{\mathrm{O}}_{\mathrm{2}}$	$\mathrm{2}\times \mathrm{3.68}\times {\mathrm{10}}^{-\mathrm{13}}\times e^{\mathrm{635}/\text{T}}$	SAR ($\mathrm{ROOH}$)	${}^{\mathrm{n}}$
125	${\mathrm{HOOCH}}_{\mathrm{2}}{\mathrm{SCH}}_{\mathrm{2}}\mathrm{OOH}$	=	${\mathrm{HOOCH}}_{\mathrm{2}}{\mathrm{SCH}}_{\mathrm{2}}\mathrm{O}$ + $\mathrm{OH}$	J(“CH3OOH”)	SAR (${\mathrm{CH}}_{\mathrm{3}}\mathrm{OOH}$)	${}^{\mathrm{b}}$
126	$\mathrm{OCS}$ + $\mathrm{O}$	=	$\mathrm{CO}$ + $\mathrm{SO}$	$\mathrm{2.1}\times {\mathrm{10}}^{-\mathrm{11}}\times e^{-\mathrm{2200}/\text{T}}$	Evaluated data	${}^{\mathrm{a}}$
127	$\mathrm{OCS}$ + $\mathrm{OH}$	=	$\mathrm{SO}$ + $\mathrm{OH}$	$\mathrm{7.2}\times {\mathrm{10}}^{-\mathrm{14}}\times e^{-\mathrm{1070}/\text{T}}$	Evaluated data	${}^{\mathrm{a}}$
128	$\mathrm{SO}$	=	${\mathrm{SO}}_{\mathrm{2}}$ + $\mathrm{O}$	$\mathrm{1.6}\times {\mathrm{10}}^{-\mathrm{13}}\times e^{-\mathrm{2280}/\text{T}}\times$ [${\mathrm{O}}_{\mathrm{2}}$]	Evaluated data	${}^{\mathrm{a}}$
129	$\mathrm{SO}$ + ${\mathrm{O}}_{\mathrm{3}}$	=	${\mathrm{SO}}_{\mathrm{2}}$	$\mathrm{3.4}\times {\mathrm{10}}^{-\mathrm{12}}\times e^{-\mathrm{1100}/\text{T}}$	Evaluated data	${}^{\mathrm{a}}$
130	$\mathrm{SO}$ + ${\mathrm{NO}}_{\mathrm{2}}$	=	${\mathrm{SO}}_{\mathrm{2}}$ + $\mathrm{NO}$	$\mathrm{1.40}\times {\mathrm{10}}^{-\mathrm{11}}$	Evaluated data	${}^{\mathrm{a}}$
131	$\mathrm{SO}$ + $\mathrm{OH}$	=	${\mathrm{SO}}_{\mathrm{2}}$ + ${\mathrm{HO}}_{\mathrm{2}}$	$\mathrm{2.6}\times {\mathrm{10}}^{-\mathrm{11}}\times e^{\mathrm{330}/\text{T}}$	Evaluated data	${}^{\mathrm{a}}$
132	${\mathrm{HOOCH}}_{\mathrm{2}}\mathrm{S}$ + ${\mathrm{O}}_{\mathrm{3}}$	=	${\mathrm{HOOCH}}_{\mathrm{2}}\mathrm{SO}$	$\mathrm{1.5}\times {\mathrm{10}}^{-\mathrm{12}}\times e^{\mathrm{360}/\text{T}}$	SAR (${\mathrm{CH}}_{\mathrm{3}}\mathrm{S}$)	${}^{\mathrm{a}}$
133	${\mathrm{HOOCH}}_{\mathrm{2}}\mathrm{S}$ + ${\mathrm{NO}}_{\mathrm{2}}$	=	${\mathrm{HOOCH}}_{\mathrm{2}}\mathrm{SO}$ + $\mathrm{NO}$	$\mathrm{3.0}\times {\mathrm{10}}^{-\mathrm{11}}\times e^{\mathrm{240}/\text{T}}$	SAR (${\mathrm{CH}}_{\mathrm{3}}\mathrm{S}$)	${}^{\mathrm{a}}$
134	${\mathrm{HOOCH}}_{\mathrm{2}}\mathrm{S}$	=	${\mathrm{HOOCH}}_{\mathrm{2}}\mathrm{SOO}$	$\mathrm{1.20}\times {\mathrm{10}}^{-\mathrm{16}}\times e^{\mathrm{1580}/\text{T}}\times$ [${\mathrm{O}}_{\mathrm{2}}$]	SAR (${\mathrm{CH}}_{\mathrm{3}}\mathrm{S}$)	${}^{\mathrm{b}}$
135	${\mathrm{HOOCH}}_{\mathrm{2}}\mathrm{SOO}$	=	TPA + ${\mathrm{HO}}_{\mathrm{2}}$	$\mathrm{7.13}\times {\mathrm{10}}^{-\mathrm{31}}\times {\text{T}}^{\mathrm{14.02}}\times e^{-\mathrm{2556}/\text{T}}$	Theory	${}^{\mathrm{k}}$
136	${\mathrm{HOOCH}}_{\mathrm{2}}\mathrm{SOO}$	=	${\mathrm{HOOCH}}_{\mathrm{2}}\mathrm{S}$	$\mathrm{1.5}\times {\mathrm{10}}^{\mathrm{5}}$	SAR (${\mathrm{CH}}_{\mathrm{3}}\mathrm{SOO}$)	${}^{\mathrm{m}}$
137	${\mathrm{HOOCH}}_{\mathrm{2}}\mathrm{SOO}$	=	${\mathrm{SO}}_{\mathrm{2}}$ + $\mathrm{HCHO}$ + $\mathrm{OH}$	$\mathrm{5}$	SAR (${\mathrm{CH}}_{\mathrm{3}}\mathrm{SOO}$)	${}^{\mathrm{m}}$
138	TPA + $\mathrm{OH}$	=	$\mathrm{OCS}$ + $\mathrm{OH}$	$\mathrm{5}\times {\mathrm{10}}^{-\mathrm{11}}\times \mathrm{0.14}$	Theory	${}^{\mathrm{k}}$
139	TPA + $\mathrm{OH}$	=	$\mathrm{OCHSOH}$ + $\mathrm{OH}$	$\mathrm{5}\times {\mathrm{10}}^{-\mathrm{11}}\times \mathrm{0.86}$	Theory	${}^{\mathrm{k}}$
140	$\mathrm{OCHSOH}$ + $\mathrm{OH}$	=	$\mathrm{OCS}$ + $\mathrm{OH}$	$\mathrm{1.40}\times {\mathrm{10}}^{-\mathrm{12}}$	Theory	${}^{\mathrm{k}}$
141	$\mathrm{OCHSOH}$	=	$\mathrm{CO}$ + ${\mathrm{HO}}_{\mathrm{2}}$ + $\mathrm{SO}$ + ${\mathrm{HO}}_{\mathrm{2}}$	J(“MTF”)	SAR (MTF)	${}^{\mathrm{b},\mathrm{i}}$
142	${\mathrm{HOOCH}}_{\mathrm{2}}\mathrm{SO}$ + ${\mathrm{O}}_{\mathrm{3}}$	=	${\mathrm{SO}}_{\mathrm{2}}$ + $\mathrm{HCHO}$ + $\mathrm{OH}$	$\mathrm{4.00}\times {\mathrm{10}}^{-\mathrm{13}}$	SAR (${\mathrm{CH}}_{\mathrm{3}}\mathrm{SO}$)	${}^{\mathrm{b}}$
143	${\mathrm{HOOCH}}_{\mathrm{2}}\mathrm{SO}$ + ${\mathrm{NO}}_{\mathrm{2}}$	=	${\mathrm{SO}}_{\mathrm{2}}$ + $\mathrm{HCHO}$ + $\mathrm{OH}$ + $\mathrm{NO}$	$\mathrm{1.2}\times {\mathrm{10}}^{-\mathrm{11}}$	SAR (${\mathrm{CH}}_{\mathrm{3}}\mathrm{SO}$)	${}^{\mathrm{b}}$
144	TPA	=	$\mathrm{OCS}$ + ${\mathrm{HO}}_{\mathrm{2}}$ + $\mathrm{OH}$	J(“CH3OOH”)	SAR (${\mathrm{CH}}_{\mathrm{3}}\mathrm{OOH}$)	${}^{\mathrm{b}}$
145	${\mathrm{HOOCH}}_{\mathrm{2}}{\mathrm{SCH}}_{\mathrm{2}}\mathrm{OH}$	=	${\mathrm{OCH}}_{\mathrm{2}}{\mathrm{SCH}}_{\mathrm{2}}\mathrm{OH}$ + $\mathrm{OH}$	J(“CH3OOH”)	SAR (${\mathrm{CH}}_{\mathrm{3}}\mathrm{OOH}$)	${}^{\mathrm{b}}$
146	${\mathrm{OCH}}_{\mathrm{2}}{\mathrm{SCH}}_{\mathrm{2}}\mathrm{OH}$	=	${\mathrm{HOCH}}_{\mathrm{2}}\mathrm{S}$ + $\mathrm{HCHO}$	KDEC	SAR (${\mathrm{CH}}_{\mathrm{3}}{\mathrm{SCH}}_{\mathrm{2}}\mathrm{O}$)	${}^{\mathrm{b}}$
147	${\mathrm{HOCH}}_{\mathrm{2}}\mathrm{S}$ + ${\mathrm{O}}_{\mathrm{3}}$	=	${\mathrm{HOCH}}_{\mathrm{2}}\mathrm{SO}$	$\mathrm{1.5}\times {\mathrm{10}}^{-\mathrm{12}}\times e^{\mathrm{360}/\text{T}}$	SAR (${\mathrm{CH}}_{\mathrm{3}}\mathrm{S}$)	${}^{\mathrm{a}}$
148	${\mathrm{HOCH}}_{\mathrm{2}}\mathrm{S}$ + ${\mathrm{NO}}_{\mathrm{2}}$	=	${\mathrm{HOCH}}_{\mathrm{2}}\mathrm{SO}$ + $\mathrm{NO}$	$\mathrm{3.0}\times {\mathrm{10}}^{-\mathrm{11}}\times e^{\mathrm{240}/\text{T}}$	SAR (${\mathrm{CH}}_{\mathrm{3}}\mathrm{S}$)	${}^{\mathrm{a}}$
149	${\mathrm{HOCH}}_{\mathrm{2}}\mathrm{S}$	=	${\mathrm{HOCH}}_{\mathrm{2}}\mathrm{SOO}$	$\mathrm{1.20}\times {\mathrm{10}}^{-\mathrm{16}}\times e^{\mathrm{1580}/\text{T}}\times$ [${\mathrm{O}}_{\mathrm{2}}$]	SAR (${\mathrm{CH}}_{\mathrm{3}}\mathrm{S}$)	${}^{\mathrm{b}}$
150	${\mathrm{HOCH}}_{\mathrm{2}}\mathrm{SOO}$	=	${\mathrm{HOCH}}_{\mathrm{2}}\mathrm{S}$	$\mathrm{1.5}\times {\mathrm{10}}^{\mathrm{5}}$	SAR (${\mathrm{CH}}_{\mathrm{3}}\mathrm{SOO}$)	${}^{\mathrm{m}}$
151	${\mathrm{HOCH}}_{\mathrm{2}}\mathrm{SOO}$	=	${\mathrm{SO}}_{\mathrm{2}}$ + $\mathrm{HCHO}$ + ${\mathrm{HO}}_{\mathrm{2}}$	$\mathrm{5}$	SAR (${\mathrm{CH}}_{\mathrm{3}}\mathrm{SOO}$)	${}^{\mathrm{m}}$
152	${\mathrm{HOCH}}_{\mathrm{2}}\mathrm{SO}$ + ${\mathrm{O}}_{\mathrm{3}}$	=	${\mathrm{SO}}_{\mathrm{2}}$ + $\mathrm{HCHO}$ + ${\mathrm{HO}}_{\mathrm{2}}$	$\mathrm{4.00}\times {\mathrm{10}}^{-\mathrm{13}}$	SAR (${\mathrm{CH}}_{\mathrm{3}}\mathrm{SO}$)	${}^{\mathrm{b}}$
153	${\mathrm{HOCH}}_{\mathrm{2}}\mathrm{SO}$ + ${\mathrm{NO}}_{\mathrm{2}}$	=	${\mathrm{SO}}_{\mathrm{2}}$ + $\mathrm{HCHO}$ + ${\mathrm{HO}}_{\mathrm{2}}$ + $\mathrm{NO}$	$\mathrm{1.2}\times {\mathrm{10}}^{-\mathrm{11}}$	SAR (${\mathrm{CH}}_{\mathrm{3}}\mathrm{SO}$)	${}^{\mathrm{b}}$
154	${\mathrm{OCH}}_{\mathrm{2}}\mathrm{SCHO}$	=	$\mathrm{HCHO}$ + $\mathrm{OCS}$ + ${\mathrm{HO}}_{\mathrm{2}}$	KDEC	SAR (${\mathrm{CH}}_{\mathrm{3}}{\mathrm{SCH}}_{\mathrm{2}}\mathrm{O}$)	${}^{\mathrm{b}}$

References: ${}^{\mathrm{a}}$ NASA Panel report . ${}^{\mathrm{b}}$ MCM v3.3.1 . ${}^{\mathrm{c}}$ . ${}^{\mathrm{d}}$ . ${}^{\mathrm{e}}$ . ${}^{\mathrm{f}}$ . ${}^{\mathrm{g}}$ . ${}^{\mathrm{h}}$ This work. ${}^{\mathrm{i}}$ . ${}^{\mathrm{j}}$ . ${}^{\mathrm{k}}$ . ${}^{\mathrm{l}}$ . ${}^{\mathrm{m}}$ . ${}^{\mathrm{n}}$ . ${}^{\mathrm{o}}$ . ${}^{\mathrm{p}}$ .Rate constants: KinfT $=\mathrm{4.1}\times {\mathrm{10}}^{-\mathrm{14}}\times (\text{T}/\mathrm{298}{)}^{\mathrm{1.8}}$, K0T $=\mathrm{1.8}\times {\mathrm{10}}^{-\mathrm{33}}\times (\text{T}/\mathrm{298}{)}^{\mathrm{2}}$, Kb $=(\mathrm{1}+({\text{log}}_{\mathrm{10}}(\text{K0T}\times \text{[M]}/\text{KinfT}){)}^{\mathrm{2}}{)}^{-\mathrm{1}}$, KinfT2 $=\mathrm{1.7}\times {\mathrm{10}}^{-\mathrm{12}}\times (\text{T}/\mathrm{298}{)}^{\mathrm{0.2}}$, K0T2 $=\mathrm{2.9}\times {\mathrm{10}}^{-\mathrm{31}}\times (\text{T}/\mathrm{298}{)}^{-\mathrm{4.1}}$, Kb2 $=(\mathrm{1}+({\text{log}}_{\mathrm{10}}(\text{K0T2}\times \text{[M]}/\text{KinfT2}){)}^{\mathrm{2}}{)}^{-\mathrm{1}}$, KRO2HO2 $=\mathrm{2.91}\times {\mathrm{10}}^{-\mathrm{13}}\times e^{\mathrm{1300}/\text{T}}$, KRO2NO3 $=\mathrm{2.3}\times {\mathrm{10}}^{-\mathrm{12}}$, J(“CH3OOH”) = photolysis rate for ${\mathrm{CH}}_{\mathrm{3}}\mathrm{OOH}$, KDEC $=\mathrm{1.0}\times {\mathrm{10}}^{\mathrm{6}}$, J(“MTF”) = photolysis rate for MTF (${\mathrm{CH}}_{\mathrm{3}}\mathrm{SCHO}$), KRO2NO $=\mathrm{2.7}\times {\mathrm{10}}^{-\mathrm{12}}\times e^{\mathrm{360}/\text{T}}$, KAPHO2 $=\mathrm{5.2}\times {\mathrm{10}}^{-\mathrm{13}}\times e^{\mathrm{980}/\text{T}}$.

The decomposition of ${\mathrm{ROONO}}_{\mathrm{2}}$ products formed from our mechanism (${\mathrm{CH}}_{\mathrm{3}}{\mathrm{SCH}}_{\mathrm{2}}{\mathrm{OONO}}_{\mathrm{2}}$, ${\mathrm{CH}}_{\mathrm{3}}{\mathrm{SOO}}_{\mathrm{2}}{\mathrm{NO}}_{\mathrm{2}}$, and ${\mathrm{CH}}_{\mathrm{3}}{\mathrm{SO}}_{\mathrm{4}}{\mathrm{NO}}_{\mathrm{2}}$; MSPN) was adjusted to improve the fit to the experiment 1. The decomposition rate constants for these ${\mathrm{ROONO}}_{\mathrm{2}}$ products in the MCM were not experimentally determined; however, the experiment is sensitive to these reactions. The decomposition rate constant for ${\mathrm{CH}}_{\mathrm{3}}{\mathrm{SO}}_{\mathrm{4}}{\mathrm{NO}}_{\mathrm{2}}$ was increased until the modified mean bias was around zero when compared to the measured product in the experiment 1. The same decomposition rate was used for ${\mathrm{CH}}_{\mathrm{3}}{\mathrm{SCH}}_{\mathrm{2}}{\mathrm{OONO}}_{\mathrm{2}}$ (which was not included as a species in the MCM) and ${\mathrm{CH}}_{\mathrm{3}}{\mathrm{SOO}}_{\mathrm{2}}{\mathrm{NO}}_{\mathrm{2}}$. Although these adjustments may not be realistic and should be adjusted further if they are experimentally determined, due to the low concentration of ${\mathrm{NO}}_{\mathrm{x}}$ in the marine environment these ${\mathrm{ROONO}}_{\mathrm{2}}$ species are not considered important for modelling DMS in the environment.

One exception to our methodology was the decomposition of ${\mathrm{CH}}_{\mathrm{3}}{\mathrm{SO}}_{\mathrm{2}}$. Although the rate constant of ${\mathrm{CH}}_{\mathrm{3}}{\mathrm{SO}}_{\mathrm{2}}$ was experimentally determined to be less than 1 s${}^{-\mathrm{1}}$ , both and used higher rate constants of 7.0 and 20 s${}^{-\mathrm{1}}$ at 298 K, respectively, to model their experiments. Additionally, calculated a decomposition rate constant of $\mathrm{2}\times {\mathrm{10}}^{\mathrm{3}}$ s${}^{-\mathrm{1}}$. Due to the large range of rate constants used and the sensitivity of the experiment 1 to the reaction, in this work the rate constant was set to the best fit. To do this, initially the estimated rate constant of the reaction between ${\mathrm{CH}}_{\mathrm{3}}{\mathrm{SO}}_{\mathrm{3}}$ and DMS from was increased by a factor of 2.1 for the ratio of sulfuric acid (${\mathrm{H}}_{\mathrm{2}}{\mathrm{SO}}_{\mathrm{4}}$; measured as sulfate) to methanesulfonic acid (MSA) modelled to be the same as found by the experiment 1 ($\mathrm{1.45}\times {\mathrm{10}}^{-\mathrm{13}}$ cm${}^{\mathrm{3}}$ molecules${}^{-\mathrm{1}}$ s${}^{-\mathrm{1}}$). Then, the rate constant of the decomposition of ${\mathrm{CH}}_{\mathrm{3}}{\mathrm{SO}}_{\mathrm{2}}$ was adjusted until the ${\mathrm{SO}}_{\mathrm{2}}$-to-MSA ratio modelled by our mechanism was the same as measured by the experiment 1 (6 s${}^{-\mathrm{1}}$ at 295 K). Both these adjustments are described in more detail in Sect. 5. These adjustments do not necessarily provide correct rate constants; however, they can be used to improve the fit to the chamber studies until more experiments are performed. Finally, the temperature dependence of ${\mathrm{CH}}_{\mathrm{3}}{\mathrm{SO}}_{\mathrm{2}}$ decomposition (and ${\mathrm{CH}}_{\mathrm{3}}{\mathrm{SO}}_{\mathrm{3}}$) was included by using the same temperature dependence as the MCM.

Table A1 lists all the sulfur reactions included in our mechanism along with the sources of those reactions. Similarly to the MCM, ${\mathrm{H}}_{\mathrm{2}}\mathrm{O}$, ${\mathrm{O}}_{\mathrm{2}}$, and ${\mathrm{CO}}_{\mathrm{2}}$ are not included as products or reactants in the mechanism. In the case that ${\mathrm{H}}_{\mathrm{2}}\mathrm{O}$ or ${\mathrm{O}}_{\mathrm{2}}$ are reactants, their concentration is included in the rate constant of the reaction (in molecules cm${}^{-\mathrm{3}}$). The mechanism assumes oxygen is a bath gas, and, as such, the formations of ${\mathrm{CH}}_{\mathrm{3}}$ and ${\mathrm{CH}}_{\mathrm{3}}{\mathrm{SCH}}_{\mathrm{2}}$ radicals are included as ${\mathrm{CH}}_{\mathrm{3}}{\mathrm{O}}_{\mathrm{2}}$ and ${\mathrm{CH}}_{\mathrm{3}}{\mathrm{SCH}}_{\mathrm{2}}{\mathrm{O}}_{\mathrm{2}}$, respectively, due to their fast reaction with oxygen. Additionally, instead of including individual ${\mathrm{RO}}_{\mathrm{2}}$ reactions, the ${\mathrm{RO}}_{\mathrm{2}}$ radicals formed in the model run are lumped into a total ${\mathrm{RO}}_{\mathrm{2}}$ concentration, which is included in the rate constant calculations of ${\mathrm{RO}}_{\mathrm{2}}$ reactions. The variable $M$ in the rate constants of Table A1 is the total bath gas concentration.

In some cases, acronyms (such as DMSO, DMSO2, and MSA) are included instead of the chemical formulas. These acronyms are given alongside their structural formula in Fig. . The exceptions to this are SA (${\mathrm{H}}_{\mathrm{2}}{\mathrm{SO}}_{\mathrm{4}}$, sulfuric acid), DMSO2O2 (${\mathrm{CH}}_{\mathrm{3}}{\mathrm{SO}}_{\mathrm{2}}{\mathrm{CH}}_{\mathrm{2}}{\mathrm{O}}_{\mathrm{2}}$), DMSO2OOH (${\mathrm{CH}}_{\mathrm{3}}{\mathrm{SO}}_{\mathrm{2}}{\mathrm{CH}}_{\mathrm{2}}\mathrm{OOH}$), DMSO2OH (${\mathrm{CH}}_{\mathrm{3}}{\mathrm{SO}}_{\mathrm{2}}{\mathrm{CH}}_{\mathrm{2}}\mathrm{OH}$), and DMSO2O (${\mathrm{CH}}_{\mathrm{3}}{\mathrm{SO}}_{\mathrm{2}}{\mathrm{CH}}_{\mathrm{2}}\mathrm{O}$).

The photolysis rate constants included in Table 1A, J(“MTF”) and J(“CH3OOH”), are wavelength-dependent and calculated for each experiment. In the case of J(“MTF”), the absorption cross-section of methyl thioformate (${\mathrm{CH}}_{\mathrm{3}}\mathrm{SCHO}$; MTF) measured by was used, with the quantum yield from the MCM (based on ${\mathrm{C}}_{\mathrm{3}}{\mathrm{H}}_{\mathrm{7}}\mathrm{CHO}$). The photolysis rate constant J(“CH3OOH”) was based on the absorption cross-section and quantum yield used by the MCM.

Code and data availability

The model input files for all the simulated experiments and output files for the model runs using the mechanism developed for this work can be found through Apollo (10.17863/CAM.101652; ).

The supplement related to this article is available online at: https://doi.org/10.5194/acp-24-3329-2024-supplement.

Author contributions

LSDJ developed the mechanism and ran the box model simulations under the supervision of ATA and CG. LSDJ and ATA wrote the paper. All authors were involved in helpful discussions and contributed to the paper.

Competing interests

The contact author has declared that none of the authors has any competing interests.

Disclaimer

Publisher’s note: Copernicus Publications remains neutral with regard to jurisdictional claims made in the text, published maps, institutional affiliations, or any other geographical representation in this paper. While Copernicus Publications makes every effort to include appropriate place names, the final responsibility lies with the authors.

Acknowledgements

The authors would like to thank Qing Ye, Matthew B. Goss, Jesse H. Kroll, Jiali Shen, Xu-Cheng He, Christopher M. Jernigan, and Timothy H. Bertram for giving additional information and data to aid in the modelling of their experiments. Lorrie S. D. Jacob acknowledges Cambridge Australia Scholarships and the Cambridge Trust for stipend and tuition support.

Financial support

This research has been supported by the Natural Environment Research Council (grant no. NE/W009412/1).

Review statement

This paper was edited by Marc von Hobe and reviewed by two anonymous referees.