1. Introduction
Air quality models are employed by policy-makers to formulate emission control strategies with a view to combatting photochemical air pollution, particulate matter and acid rain in Europe, Asia and North America. Chemical mechanisms are essential components of air quality models because ozone (O
3) is a secondary pollutant that is not emitted into the atmosphere. Ozone is formed by photochemical reactions in the sunlit atmosphere from emitted precursors: oxides of nitrogen (NO
x) and volatile organic compounds (VOCs). Hence air quality policies and strategies, and the models that address them, have to focus on the NO
x and VOC precursors. Tackling elevated ozone levels is a major policy activity because ozone is an important air pollutant that at elevated levels damages human health and vegetation and also contributes to climate change [
1].
The chemical mechanisms in air quality models incorporate information on the chemical kinetics and pathways that transform the primary emitted precursor pollutants into secondary pollutants, particularly ozone and suspended particulate matter. It is essential that these chemical mechanisms faithfully represent the actual chemistry of the real world if the derived policies are to deliver the required improvements in air quality. If the chemical mechanisms contain inadequately characterised representations of important atmospheric chemistry processes, then the policy predictions may underestimate the emissions reductions required or may overstate them, causing unnecessary implementation costs. How are chemical mechanisms to be evaluated? Perhaps more importantly, how are chemical mechanisms to be chosen for particular air quality policy modelling applications from the plethora of available chemical mechanisms? Chemical mechanism selection has been a difficult issue for regulatory policy development since the early days [
2]. These concerns still exist [
3,
4]. These difficulties are compounded by the often limited choice of chemical mechanisms available within the widely used air quality modelling systems because chemical mechanisms are not offered on a plug-and-play basis as with meteorological parameterisations, for example. Where differences in policy-relevant predictions have been found by substituting different chemical mechanisms as in, for example [
5], there is then the difficulty of proving which predictions are correct.
In previous chemical mechanism intercomparison studies, the impact of different chemical mechanisms for air quality policy formulation have been investigated under North American [
6] and European [
7] conditions. Here the focus is more global in context, rather than regional, as policy attention has shifted to the background troposphere. This change in emphasis reflects the shift in the attention of policy-makers to address the intercontinental transport of ozone [
8,
9,
10,
11], the issue of the policy-relevant ozone background across the United States of America [
12,
13,
14,
15,
16,
17,
18,
19,
20,
21] the global burden of disease resulting from air pollution [
22,
23] and the linkage between air quality and climate change [
11,
24,
25,
26,
27,
28,
29].
There are several aspects of chemical mechanism development that could be compared and assessed in an intercomparison. In this study, attention is focussed on the representation of the oxidation and degradation of the organic compounds and its relationship to the photochemical formation of ozone and hydroxyl (OH) radicals. The representation of these processes presents a formidable challenge to mechanism developers because of their complexity, because of the myriad of organic compounds that need to be considered and because of the limited nature of current understanding. Mechanism developers inevitably have to rely on approximations, assumptions and in many cases, important processes have to be neglected altogether. To enable a fair intercomparison of the representation of the oxidation and degradation chemistry of organic compounds in each chemical mechanism, a common basis has to be established so that the main differences between the chemical mechanisms are highlighted and brought to the fore. Attention is given to the photochemical ozone production rates and the OH number densities in the base case and their responses to incremental additions of organic compounds.
2. Methodology
A zero-dimensional box model was set up to provide a framework for the intercomparison and evaluation of the thirteen chemical mechanisms. The focus of the intercomparison was on O3 and OH radicals, together with other reactive free radical species. In the paragraphs below, the background environmental conditions underpinning the intercomparison are described. Attention is then directed to the implementation of the chemical mechanisms themselves. Finally, a description is given of the chemical development in the box model and how the mixing ratios of the box model species were constrained to the background environmental conditions.
The formulation of the box model was based on the Photochemical Trajectory Model (PTM), the details of which are given elsewhere [
30]. In this application, wet and dry deposition, exchange with the free troposphere and emission processes have been switched off, leaving the complete focus on the chemical development of photochemical ozone and that of the hydroxyl (OH) radical species that drives it.
Table 1 summarises the details of the thirteen chemical mechanisms studied, together with their literature references. They varied marked in complexity from the highly detailed and explicit Master Chemical Mechanism (MCM) v3.3.1 to the highly condensed Carbon Bond and SAPRC mechanisms. Our focus was principally on the background troposphere and the regional and global chemistry-transport models used in policy formulation and assessment on the regional and hemispheric scales. Each of the major chemical mechanisms studied here has undergone frequent updating and no attempt has been made to keep track of these. There is no implication underpinning our choice of mechanism and version that we consider them to be of more policy importance than any other. For most chemical mechanisms, there are a range of versions spanning many years of development and, again, our choice of version is completely arbitrary.
The mechanisms were, however, not implemented as published but were harmonised to minimise the influence of publication date which spanned four decades from RADM to CRIv2.2. The first harmonisation addressed the so-called ‘inorganic chemistry’. This set of chemical reactions establishes the fast photochemical balance and involves the hydroxyl (OH), hydroperoxy (HO
2) and oxygen (O
1D and O
3P) atoms and their reactions with nitric oxide (NO), nitrogen dioxide (NO
2), ozone (O
3), water vapour, carbon monoxide (CO), hydrogen (H
2) and sulphur dioxide (SO
2). The ‘inorganic chemistry’ provided with each mechanism was removed and replaced with a set of 49 chemical reaction pathways and rate coefficients, together with their temperature and pressure dependences, taken from IUPAC (
http://iupac.pole-ether.fr/) [
42].
The second harmonisation involved the photolysis rate coefficients. The photolysis rate coefficients provided with each mechanism were replaced with a standard set taken from the MCM website (
http://mcm.york.ac.uk/MCM) [
33]. The final harmonisation involved the rate coefficients for the formation and decomposition of the peroxyacyl nitrates and for the reactions of the methyl peroxy (CH
3O
2) radical. Again, these were removed and replaced with a standard set of rate coefficients, together with their pressure and temperature dependences taken from the MCM website (
http://mcm.york.ac.uk/MCM) [
33].
It is understood that these harmonisation steps may well move the mechanisms away from the conditions and chemical regimes under which they were developed by their originators. This was considered inevitable. As a result, the performance of the mechanisms may be different from that if no changes had been made. To draw attention to the harmonisation and standardisation steps, the names of the mechanisms have been printed in italics to indicate that they have not been implemented as originally developed.
Any intercomparison of chemical mechanisms needs input data on background environmental conditions to set up an appropriate chemical regime to frame the evaluation. In this study, output has been taken from a global Lagrangian chemistry-transport model (STOCHEM-CRI) [
43] and was used to provide realistic mixing ratios for a number of trace gases in the box model. These trace gases, of which there are thirty in all see
Table 2, included closed-shell ozone precursors and reaction products with atmospheric lifetimes of the order of minutes and longer. Free radical species have much shorter lifetimes and were not set up in the same way but allowed to establish their own levels based on the time-dependent photochemical activity.
The background environmental data were taken from the STOCHEM-CRI model for a grid of 108 points at 20° intervals, covering a region between 60° N and 60° S, near-ground, 1st July and are summarised in
Table 2. For further details of the STOCHEM-CRI model [
43]. Not all of these species were represented in each chemical mechanism and so care had to be taken not to introduce bias into the model results with any chemical mechanism through the choice of background environmental conditions. Since some mechanisms relied heavily on lumped VOC emission surrogates, VOC coverage was limited to common VOCs and to those explicitly or near-explicitly represented in all chemical mechanisms. The background data therefore necessarily represented the minimum or lowest common denominator approach so as not the introduce bias between the highly detailed and highly reduced chemical mechanisms.
In a box model, a differential equation of the form
was set up for each model species, i, where c
i is the concentration of species, i, in the box, P
i is its production rate from emissions, chemistry and boundary conditions and l
i is the first order loss coefficient arising from chemistry, deposition and boundary conditions. In a constrained box model, the above differential equation is modified by the addition of a flux, F
i, to the right hand side of the equation so that the rate of change of the constrained species, i, as listed in
Table 2, remains zero and its concentration remains constant at the constrained value, c
i*:
The Gear’s method automatic numerical integrator FACSIMILE [
44] returns the flux, F
i, required at the end of each time step to maintain the concentrations of each species at its constraint. This flux, which is diurnally-varying through the time-dependent photolysis rates, is then integrated over a time period of 5 days to give the time-integrated production or loss flux (depending on its sign) for that species. The output of the constrained box model is therefore the time-integrated production or loss fluxes for each of the thirty constrained species in
Table 2. Particular focus was given in this study to the time-integrated ozone production flux, P
O3.
For the many species without constrained values in
Table 1, these species would reach some form of local instantaneous photochemical steady state and their time varying concentrations would be set by equations of the form of (1) above. Their concentrations were averaged over the 5-day time period and provided another component of the output from the constrained box model. Particular focus was given to the 5-day average hydroxyl radical number density, [
OH], in this study.
3. Testing and Evaluating the Constrained Box Model
The constrained box model was set up with each of the thirteen reduced chemical mechanisms and each set of 108 sets of background conditions and was integrated for five days. For each set of conditions and each chemical mechanism, the P
O3 and [
OH] values were noted, and these comprised the base case results for each chemical mechanism. P
O3 for MCMv3.3.1 varied from −0.5 ppt per hour to 11 ppb per hour over the 108 sets of conditions, with an average of 1.1 ppb per hour. The P
O3 values were closely similar between the different chemical mechanisms and closely tracked each other across the different sets of conditions. Scatter plots were constructed between the base case P
O3 results for each reduced chemical mechanism versus the MCMv3.3.1 results. Linear regression equations were fitted through each scatter plot, revealing excellent correlations with R
2 values consistently greater than 0.99, slopes that were close to unity and intercepts that were statistically indistinguishable from zero. P
O3 values closely overlapped and average intermechanism differences from the MCMv3.3.1 were within the range −0.6% to +0.4% across all 108 sets of background conditions and are plotted out in
Figure 1.
[
OH] for the MCMv3.3.1 varied from 0.3 to 3.6 × 10
6 molecule cm
−3, with an average of 1.6 × 10
6 molecule cm
−3. Again, [
OH] values were closely similar between the different chemical mechanisms. However, there was much more variability between the chemical mechanisms and sets of conditions compared to P
O3. Scatter plots were constructed between [
OH] for each chemical mechanism versus MCMv3.3.1 and linear regression equations were fitted through the scatter plots. Correlations were good with R
2 values generally greater than 0.98. The slopes and intercepts were generally indistinguishable from zero. Average intermechanism differences were generally in the range −4% to +8%, see
Figure 2, (RADM results excluded).
This close agreement in PO3 and [OH] between the chemical mechanisms confirmed the view that the standardisation and harmonisation steps had not introduced significant bias between the chemical mechanisms.
4. Comparing the Impact of Organic Compounds on Background Tropospheric Chemistry
4.1. Overview
For the purposes on this intercomparison, we consider that organic compounds exert two major impacts on tropospheric chemistry. In the first, the hydroxyl radical driven degradation of an organic compound leads to the formation of hydroperoxy and organic peroxy radicals (RO2) which oxidise nitric oxide (NO) to nitrogen dioxide (NO2), and hence producing ozone (O3 This impact is most important in urban areas and in the polluted atmospheric boundary layer where there are abundant sources of nitrogen oxides (NOx = NO + NO2). However, in the background troposphere, the subject of this intercomparison, it is of diminished policy importance because background tropospheric chemistry is generally NOx-limited.
In the second major impact of organic compounds on tropospheric chemistry, the hydroxyl radical driven degradation of organic compounds changes the sources and sinks of the major free radical species: hydroxyl (OH), hydroperoxy (HO2) and organic peroxy (RO2) radicals. Since background tropospheric chemistry is generally free radical limited then these changes may be of some significance. Changes to the sources and sinks of the hydroxyl radical have an important impact on tropospheric chemistry because of the vast numbers of trace gases whose sinks are controlled by OH radical driven oxidation and degradation. This second major impact of organic compounds on tropospheric chemistry is therefore of immediate policy interest and is the focus of this intercomparison.
4.2. Methodology
The methodology chosen to address the impact of organic compounds on the major free radical species involved adding increments of each organic compound and following the changes in hydroxyl radical number densities ([
OH]) and the rates of oxidation of organic compounds (ROC). Experiments were performed with 100 ppt (1 ppt = 1 part in 10
12 parts of air) increments of each organic compound from 100 ppt up to 500 ppt at each of the 108 background locations. This methodology gave four differences (Δ): 100 ppt vs. 200 ppt, 200 ppt vs. 300 ppt and so on. For each difference, the following hydroxyl radical incremental reactivity (IROC) was defined and estimated:
where [
OH] is the five-day average OH number density in molecule cm
−3 and ROC is the time-averaged oxidation rate of the organic compound in ppb per hour.
To illustrate the methodology, let us consider an arbitrary example employing ethylene (ethene) (C2H4) as the organic compound, the most remote of the 108 background locations and the MCMv3.3.1 mechanism. With an additional 100 ppt increment of ethylene added, the time-averaged ethylene oxidation rate increased by 2.62 ppt per hour and increased to 5.01 ppt per hour with 200 ppt added. At the same time, the five-day averaged OH number density declined from 4.91 × 105 to 4.73 × 105 molecule cm−3. On this basis, IRC2H4100 ppt vs. 200 ppt was found to be (4.73–4.91) × 105 / (5.01–2.61) × 10−3 molecule cm−3 per ppb hour−1, that is to say: IRC2H4100ppt vs. 200ppt was −7.5 × 106 molecule cm−3 per ppb hour−1. When averaging over all 108 background locations, IRC2H4100 ppt vs. 200 ppt was found to be −7.2 × 106 and this decreased to −6.9 × 106, −6.7 × 106, −6.5 × 106 and −6.5 × 106 molecule cm−3 per ppb hour−1 for 200 ppt vs. 300 ppt, 300 ppt vs. 400 ppt and 400 ppt vs. 500 ppt, respectively. The decline from −7.2 to −6.5 × 106 molecule cm−3 per ppb hour−1 illustrated that the hydroxyl radical responses were not accurately linear in ethylene increments and that there was some non-linearity in the reacting system amounting to ±10%. Overall, it was concluded that using the MCMv3.3.1 mechanism, hydroxyl radical number densities declined with increasing ethylene additions, with an ‘initial’, (that is to say, 100 ppt vs. 200 ppt) IRC2H4 of −7.2 × 106 molecule cm−3 per ppb hour−1.
This same methodology was then applied to all thirteen organic compounds and all thirteen chemical mechanisms. The results for all the alkanes in each chemical mechanism are presented in
Table 3, man-made alkenes in
Table 4, biogenic alkenes and dialkenes in
Table 5 and aromatics in
Table 6. A rapid overview of
Table 3,
Table 4,
Table 5 and
Table 6 revealed that all the organic compounds studied gave negative initial hydroxyl radical incremental reactivities. Generally speaking, incremental reactivities were smallest for the alkane class in
Table 3 and largest for the aromatics in
Table 6. There were, however, significant divergences between the incremental reactivities found by the different chemical mechanisms. Understanding these divergences is the main aim of this chemical mechanism intercomparison.
4.3. Assessing the Impact of the Alkanes on Hydroxyl Radicals
Table 3 presents the initial hydroxyl radical incremental reactivities for the three simplest alkanes using the thirteen chemical mechanisms. Incremental reactivities were the most negative for n-butane and least negative for ethane in almost all of the mechanisms. Averages and standard deviations of the estimates are presented in
Table 3. There was a spread of values for each organic compound and, generally speaking, the values for MCMv3.3.1 were close to the mechanism average. It was concluded that most of the mechanisms were able to account reasonably accurately for the impact of alkane increments on hydroxyl radical number densities. However, incremental reactivities did span over a factor of two between
CBM4 and
RACM-2.
The mechanism for the hydroxyl radical-driven oxidation of ethane (C
2H
6) laid out in the MCMv3.3.1 involves 46 chemical species, 119 chemical and photochemical reactions producing 21 stable product species including alkyl hydroperoxides, nitrates, alcohols, aldehydes and peroxyacyl nitrates. The complete mechanism represents a cascade system through the organic peroxy radicals, ethyl peroxy, acetyl peroxy, methyl peroxy and hydroperoxyl, after the initial attack of OH on ethane:
There is therefore a single reaction which is unique to ethane and a core of 118 reactions that are used in many other VOC oxidation mechanisms. Almost all of the 119 reactions have some impact on hydroxyl radical number densities because they may be sources and sinks of free radicals or convert one free radical into another. Overall, however, the oxidation of ethane is a net sink for hydroxyl radical, as shown by the negative hydroxyl radical incremental reactivities in
Table 3.
The twelve reduced or condensed mechanisms, see
Table 3, appear to have reduced or condensed markedly the detail of ethane oxidation into a handful of reactions without losing any of the chemical complexity and all confirm that ethane oxidation is a net sink for hydroxyl radicals. In the
CB05 mechanism, for example, the reaction of hydroxyl radicals with ethane (ETHA) is represented by:
where the generic organic peroxy radical XO2 and XO2N are shared by the degradation schemes of several hydrocarbons and the non-stoichiometric terms: 0.991 and 0.009 represent the net of the peroxy radical cascade system. In the
SAPRC-99 mechanism, the oxidation of ethane is represented with a single generic peroxy radical product, RO2_R as:
At the other end of the reduced mechanism spectrum, the CRIv2.2 mechanism uses 29 species and 68 chemical reactions to describe ethane oxidation in a manner that is accurately consistent with that in the MCMv3.3.1.
The initial hydroxyl radical incremental reactivities for propane, see
Table 3, were about three times more negative than those for ethane. Increments of propane therefore showed a significantly greater propensity to deplete 5-day mean hydroxyl radical number densities compared to those of ethane. This increased the propensity results from the longer organic peroxy radical cascade between propyl peroxy and hydroperoxyl. The MCMv3.3.1 oxidation scheme for propane employs 90 chemical species and 256 chemical reactions whilst that for the CRIv2.2 employs 49 and 129, respectively.
The reduced mechanisms gave a reasonably accurate account of the increased incremental hydroxyl radical reactivity for propane relative to ethane, showing that they were able to account for the significant increase in complexity in moving from ethane to propane. This increased complexity arises not only from the increased length of the peroxy radical cascade but also from the formation of relatively unreactive ketones such as acetone, compared with the highly reactive and photochemically-labile aldehydes such as acetaldehyde formed in ethane oxidation.
Initial hydroxyl radical incremental reactivities for n-butane are considerably more negative than those of propane and ethane, see
Table 3. There is a noticeably increased range in the available values for n-butane compared with those for ethane which amounts to a factor three. The MCMv3.3.1 mechanism represents the oxidation of n-butane using 184 species in 550 reactions, representing a dramatic increase in complexity over and above the oxidation of propane. This follows on from the increased number of isomeric peroxy radicals and stable products, such as ketones and hydroperoxides, and the increased length of the peroxy radical cascade system.
Before moving on from the simple alkanes, we note that the initial hydroxyl radical incremental reactivities for n-butane, propane and ethane have the relative ratios: 1:0.85:0.2. The corresponding ratios of their OH+ alkane rate coefficients are: 1:0.47:0.11, whereas the lengths of their oxidation mechanisms in the MCMv3.3.1 are: 1:0.47:0.21, showing no simple relationship with the incremental reactivities.
It therefore appears that in constructing mechanisms for alkane oxidation from smog chamber data or by compiling near-explicit mechanisms from individual elementary reactions, it has been possible to describe the impacts of alkanes on both ozone and hydroxyl radicals in ways that are mutually consistent. We have been able to show here that the descriptions of the impacts on photochemical ozone production rates and on 5-day mean hydroxyl radical number densities are reasonably consistent across the mechanisms. Examples will be found later where this degree of mutual consistency between photochemical ozone production rates and mean hydroxyl radical number densities is not so clearly demonstrated.
4.4. Assessing the Impact of the Alkenes on Hydroxyl Radicals
The hydroxyl radical incremental reactivities for the man-made alkenes are presented in
Table 4. The incremental reactivities for ethylene and propylene (propene) are similar in almost all the mechanisms (except for
CBM4 and
RADM) and appear to be three times more negative compared with those for trans but-2-ene. The relative ratios of the incremental reactivities for ethylene propylene trans, but-2-ene were found to be 1.0:1.0:0.3 compared to 1.0:3.2:7.8 for their OH rate coefficients. Since the alkenes are all highly reactive on the 5-day timescale of this study, it is not surprising that the OH rate coefficients give little indication of the relative hydroxyl radical incremental reactivities. In comparison, the relative numbers of reactions in the MCMv3.3.1 oxidation schemes were 1.0:2.1:1.5, respectively.
Almost all the mechanisms indicated a strong depletion of hydroxyl radical number densities driven by additional increments of ethylene, see
Table 4, except for the
RADM mechanism, the oldest mechanism studied, which appears to be an outsider in this instance. The remaining mechanisms showed excellent agreement and a relatively small range about the MCMv3.3.1 value. The similarity between the incremental reactivities found for ethylene and propylene in
Table 4 was somewhat surprising, as was the observations that this close similarity applied to almost all of the mechanisms studied.
There appeared to be a significantly large range in the estimates of the hydroxyl radical incremental reactivities for trans but-2-ene, one of the most reactive of all organic compounds. Despite this exceptional OH-reactivity, its incremental hydroxyl radical incremental reactivity was found to be significantly less than that for either ethylene or propylene, see
Table 4. The range in the estimates was also significantly larger than for the other alkenes, with the MCMv3.3.1 value towards the lower end of the range of estimates.
All three alkenes are expected to generate similar first-generation reaction products, namely; ethylene oxidation generates two formaldehyde molecules, propylene generates one molecule of formaldehyde and one of acetaldehyde and trans but-2-ene generates two molecules of acetaldehyde. Hence the less negative incremental reactivity and its greater uncertainty range of the estimates for trans but-2-ene are, at first sight, somewhat surprising. We explain these discrepancies as symptoms of the difficulty mechanism developers have in reconciling both photochemical ozone production and free radical sources and sinks, in the face of incomplete and conflicting understanding of the elementary reactions involved in free radical production and loss in the OH-oxidation of trans but-2-ene.
4.5. Assessing the Impact of the Biogenic Alkenes and Alkadienes on Hydroxyl Radicals
Table 5 presents the hydroxyl radical incremental reactivities for the biogenic alkenes and alkadienes. The biogenic alkenes have been separated from the man-made alkenes because of the different environments in which their atmospheric reactions become of policy significance. The biogenic alkenes are emitted in largely rural environments where NO
x levels are low, compared to the man-made alkenes that are emitted mainly in urban environments with associated high NO
x levels. Smog chamber studies have been largely performed under high NO
x conditions and so the preparation of mechanisms for biogenic alkenes involves significantly greater extrapolation to conditions appropriate to the atmosphere. It is possible therefore that OH oxidation mechanisms for biogenic alkenes are necessarily somewhat more uncertain than for man-made alkenes.
The hydroxyl radical incremental reactivities for isoprene in
Table 5 found with the different mechanisms span a factor of two range from −4.9 for
CBM-4 to −9.8 for
MOZART-4, with an average of −7.8 ± 1.4. The MCMv3.3.1 value lies at the lower end of the 1-σ confidence range with the CRIv2.2 value close to the middle of the confidence range. Overall, almost all of the mechanisms gave incremental reactivities for isoprene that were more negative than the MCMv3.3.1 value, with only
CBM-4 returning a less negative value.
In contrast, hydroxyl radical incremental reactivities for α-pinene were significantly less negative than those for isoprene and significantly more accurately defined as indicated by the 1-σ confidence range. For this organic compound, the difference between the MCMv3.3.1 and CRIv2.2 mechanisms was small and significant. The MOZART-4 mechanism gave the most negative reactivity whilst CB6 gave the least negative.
Hydroxyl radical incremental reactivities for β-pinene were the most negative of all the three biogenic alkenes and alkadienes. However, coverage by the available mechanisms was poor and results were limited to those from the MCMv3.3.1 and CRIv2.2 mechanisms. Both mechanisms made β-pinene between two and three times more negative than α-pinene.
4.6. Assessing the Impact of the Acetylene and the Aromatics on Hydroxyl Radical Number Densities
Incremental reactivities for the acetylene (ethyne) are presented in
Table 6. Acetylene is generally considered unreactive because of its low reactivity with OH. Estimates of its hydroxyl radical incremental reactivity are available from six mechanisms and they demonstrate that acetylene increments have a relatively weak influence on hydroxyl radical number densities. Incremental reactivities from
CB6 and
RACM-2 lie within the ranges of those from the CRI and MCM mechanisms.
The incremental reactivities estimated for the aromatics in
Table 6 present a complex picture which is difficult to resolve. Starting with benzene, we find that estimates span a wide range from
MELCHIOR-2, the least negative, to
SAPRC-99 the most negative. The CRIv2.2 estimate was found to be significantly less negative than the mechanism-average and the MCMv3.3.1 estimate significantly more negative. The
RACM-2 and
SAPRC-07 estimates were close to the average.
There was no agreement concerning the hydroxyl radical incremental reactivities for toluene (methylbenzene) and o-xylene (1:2-dimethylebenzene). There was a factor of three range between the CRIv2.2 and the MCMv3.3.1 estimates which was unusual since both mechanisms are generally close. All mechanisms gave estimates that were less negative for both toluene and o-xylene compared to the MCMv3.3.1 and many were less negative than the CRIv2.2 mechanism. The MCMv3.3.1 mechanism therefore appeared to be an outsider for both toluene and o-xylene.
5. Discussions and Conclusions
This intercomparison has taken twelve chemical mechanisms and compared how they treat VOC oxidation and degradation and its relationship to the photochemical formation of ozone and hydroxyl radicals. Representation of VOC chemistry represents a formidable challenge for mechanism developers because of the complexity involved and because of the sheer number of organic compounds and chemical reactions that require treatment. Mechanism developers have to rely on approximations, assumptions and, in many cases, important processes have to be neglected. Here, we have looked in some detail at the incremental responses of hydroxyl radicals to incremental additions of a range of organic compounds under conditions appropriate to the background atmosphere.
Most of the chemical mechanisms studied demonstrated increasingly negative incremental hydroxyl radical reactivities with increasing carbon numbers for the alkanes ethane, propane and n-butane. That is to say, increasing alkane oxidation rates led to increased hydroxyl radical depletion, with the depletion caused by n-butane > propane > ethane.
Hydroxyl radical incremental reactivities for the simple alkenes ethylene and propylene, were reasonably consistent across the chemical mechanisms studied. However, this consistent representation did not extend to trans but-2-ene where reactivity estimates spanned a range of a factor of five. Incremental reactivities were reasonably well-defined for isoprene which was encouraging in view of its importance to background tropospheric chemistry.
The most serious discrepancies emerging from this study were found with the aromatics toluene and o-xylene and with the Master Chemical Mechanism. These discrepancies are in part explained by the Master Chemical Mechanism being an assembly of reported information on elementary reactions which has not been designed for explaining aromatic + NO
x smog chamber data accurately. For the Master Chemical Mechanism to do this, an additional hydroxyl radical source had to be included [
45]. Because the reduced mechanisms effectively included this additional hydroxyl radical source, they necessarily gave significantly less negative incremental hydroxyl radical reactivities for toluene and o-xylene in this intercomparison.
There is a general level of uncertainty in the understanding and representation of aromatic degradation chemistry in chemical mechanisms and their ability to represent photochemical ozone formation relative to hydroxyl radical depletion in smog chamber experiments. In view of the lack of agreement between the mechanisms, it is difficult to comment on the relative incremental reactivities for the aromatics in the real world. Further laboratory and smog chamber studies will be required to resolve the differences between the Master Chemical Mechanism and the condensed and reduced smog chamber mechanisms.