scholarly journals Modelling of the photooxidation of toluene: conceptual ideas for validating detailed mechanisms

2003 ◽  
Vol 3 (1) ◽  
pp. 89-106 ◽  
Author(s):  
V. Wagner ◽  
M. E. Jenkin ◽  
S. M. Saunders ◽  
J. Stanton ◽  
K Wirtz ◽  
...  

Abstract. Toluene photooxidation is chosen as an example to examine how simulations of smog-chamber experiments can be used to unravel shortcomings in detailed mechanisms and to provide information on complex reaction systems that will be crucial for the design of future validation experiments. The mechanism used in this study is extracted from the Master Chemical Mechanism Version 3 (MCM v3) and has been updated with new modules for cresol and g-dicarbonyl chemistry. Model simulations are carried out for a toluene-NOx experiment undertaken at the European Photoreactor (EUPHORE). The comparison of the simulation with the experimental data reveals two fundamental shortcomings in the mechanism: OH production is too low by about 80%, and the ozone concentration at the end of the experiment is over-predicted by 55%. The radical budget was analysed to identify the key intermediates governing the radical transformation in the toluene system. Ring-opening products, particularly conjugated g-dicarbonyls, were identified as dominant radical sources in the early stages of the experiment. The analysis of the time evolution of radical production points to a missing OH source that peaks when the system reaches highest reactivity. First generation products are also of major importance for the ozone production in the system. The analysis of the radical budget suggests two options to explain the concurrent under-prediction of OH and over-prediction of ozone in the model: 1) missing oxidation processes that produce or regenerate OH without or with little NO to NO2 conversion or 2) NO3 chemistry that sequesters reactive nitrogen oxides into stable nitrogen compounds and at the same time produces peroxy radicals. Sensitivity analysis was employed to identify significant contributors to ozone production and it is shown how this technique, in combination with ozone isopleth plots, can be used for the design of validation experiments.

2002 ◽  
Vol 2 (4) ◽  
pp. 1217-1259 ◽  
Author(s):  
V. Wagner ◽  
M. E. Jenkin ◽  
S. M. Saunders ◽  
J. Stanton ◽  
K. Wirtz ◽  
...  

Abstract. Toluene photooxidation is chosen as an example to examine how simulations of smog-chamber experiments can be used to unravel shortcomings in detailed mechanisms and to provide information on complex reaction systems that will be crucial for the design of future validation experiments. The mechanism used in this study is extracted from the Master Chemical Mechanism Version 3 (MCM v3) and has been updated with new modules for cresol and g-dicarbonyl chemistry. Model simulations are carried out for a toluene-NOx experiment undertaken at the European Photoreactor (EUPHORE). The comparison of the simulation with the experimental data reveals two fundamental shortcomings in the mechanism: OH production is too low by about 80%, and the ozone concentration at the end of the experiment is over-predicted by 55%. The radical budget was analysed to identify the key intermediates governing the radical transformation in the toluene system. Ring-opening products, particularly conjugated g-dicarbonyls, were identified as dominant radical sources in the early stages of the experiment. The analysis of the time evolution of radical production points to a missing OH source that peaks when the system reaches highest reactivity. First generation products are also of major importance for the ozone production in the system. The analysis of the radical budget suggests two options to explain the concurrent under-prediction of OH and over-prediction of ozone in the model: 1) missing oxidation processes that produce or regenerate OH without or with little NO to NO2 conversion or 2) NO3 chemistry that sequesters reactive nitrogen oxides into stable nitrogen compounds and at the same time produces peroxy radicals. Sensitivity analysis was employed to identify significant contributors to ozone production and it is shown how this technique, in combination with ozone isopleth plots, can be used for the design of validation experiments.


2011 ◽  
Vol 11 (10) ◽  
pp. 28851-28894 ◽  
Author(s):  
Y. Kanaya ◽  
A. Hofzumahaus ◽  
H.-P. Dorn ◽  
T. Brauers ◽  
H. Fuchs ◽  
...  

Abstract. A photochemical box model constrained by ancillary observations was used to simulate OH and HO2 concentrations for three days of ambient observations during the HOxComp field campaign held in Jülich, Germany in July 2005. OH and HO2 levels, observed by four and three instruments, respectively, were fairly well reproduced to within 33% by a base model run (Regional Atmospheric Chemistry Mechanism with updated isoprene chemistry adapted from Master Chemical Mechanism ver. 3.1) with high R2 values (0.72–0.97) over a range of isoprene (0.3–2 ppb) and NO (0.1–10 ppb) mixing ratios. Adding isomerization of isoprene peroxy radicals to the model increased OH and HO2 by 43% and 48% on average. Although these are still only 15% and 21% higher than the observations made by one of the instruments, larger overestimations (>60%) occurred with respect to the observations made by the other three instruments, suggesting that the rates of the isomerization were not readily supported by the ensemble of radical observations. These model runs tend to underestimate observed OH reactivity which may be explained by unmeasured hydrocarbon species. By selecting hydrocarbon types to be added to the model in amounts that accounted for the missing fractions of observed OH reactivity, the gaps between HOx observations and model results with and without isomerization could be individually diminished to within uncertainty levels. In this case, however, the HO2/OH ratio rose on addition of hydrocarbons and diverged from observations. In the case where we used modeled HO2(*), taking into account the sensitivity toward speciated RO2 (organic peroxy) radicals, as recently reported from one of the participating instruments in the HO2 measurement mode, the model's overestimation for HO2 became evident (by factors of more than 1.8). These results strongly indicated that more loss processes for peroxy radicals were necessary to explain the observations. One of the measurement days was characterized by low isoprene concentrations (~0.5 ppb) and OH reactivity that was well explained by the observed species, especially before noon. For this selected period, as opposed to the general behavior, the model tended to underestimate HO2 (and HO2(*)) with respect to observations made by the three instruments. We found that this tendency is associated with high NOx concentrations, suggesting that some HO2 production or HO2 regeneration processes under high NOx conditions were being overlooked; this might require revision of ozone production regimes.


2012 ◽  
Vol 12 (5) ◽  
pp. 2567-2585 ◽  
Author(s):  
Y. Kanaya ◽  
A. Hofzumahaus ◽  
H.-P. Dorn ◽  
T. Brauers ◽  
H. Fuchs ◽  
...  

Abstract. A photochemical box model constrained by ancillary observations was used to simulate OH and HO2 concentrations for three days of ambient observations during the HOxComp field campaign held in Jülich, Germany in July 2005. Daytime OH levels observed by four instruments were fairly well reproduced to within 33% by a base model run (Regional Atmospheric Chemistry Mechanism with updated isoprene chemistry adapted from Master Chemical Mechanism ver. 3.1) with high R2 values (0.72–0.97) over a range of isoprene (0.3–2 ppb) and NO (0.1–10 ppb) mixing ratios. Daytime HO2(*) levels, reconstructed from the base model results taking into account the sensitivity toward speciated RO2 (organic peroxy) radicals, as recently reported from one of the participating instruments in the HO2 measurement mode, were 93% higher than the observations made by the single instrument. This also indicates an overprediction of the HO2 to OH recycling. Together with the good model-measurement agreement for OH, it implies a missing OH source in the model. Modeled OH and HO2(*) could only be matched to the observations by addition of a strong unknown loss process for HO2(*) that recycles OH at a high yield. Adding to the base model, instead, the recently proposed isomerization mechanism of isoprene peroxy radicals (Peeters and Müller, 2010) increased OH and HO2(*) by 28% and 13% on average. Although these were still only 4% higher than the OH observations made by one of the instruments, larger overestimations (42–70%) occurred with respect to the OH observations made by the other three instruments. The overestimation in OH could be diminished only when reactive alkanes (HC8) were solely introduced to the model to explain the missing fraction of observed OH reactivity. Moreover, the overprediction of HO2(*) became even larger than in the base case. These analyses imply that the rates of the isomerization are not readily supported by the ensemble of radical observations. One of the measurement days was characterized by low isoprene concentrations (∼0.5 ppb) and OH reactivity that was well explained by the observed species, especially before noon. For this selected period, as opposed to the general behavior, the model tended to underestimate HO2(*). We found that this tendency is associated with high NOx concentrations, suggesting that some HO2 production or regeneration processes under high NOx conditions were being overlooked; this might require revision of ozone production regimes.


2019 ◽  
Author(s):  
Michelle L. Lew ◽  
Pamela S. Rickly ◽  
Brandon P. Bottorff ◽  
Sofia Sklaveniti ◽  
Thierry Léonardis ◽  
...  

Abstract. Reactions of the hydroxyl (OH) and peroxy radicals (HO2 and RO2) play a central role in the chemistry of the atmosphere. In addition to controlling the lifetimes of many trace gases important to issues of global climate change, OH radical reactions initiate the oxidation of volatile organic compounds (VOCs) which can lead to the production of ozone and secondary organic aerosols in the atmosphere. Previous measurements of these radicals in forest environments characterized by high mixing ratios of isoprene and low mixing ratios of nitrogen oxides (NOx) have shown serious discrepancies with modeled concentrations. These results bring into question our understanding of the atmospheric chemistry of isoprene and other biogenic VOCs under low NOx conditions. During the summer of 2015, OH and HO2 radical concentrations as well as total OH reactivity were measured using Laser-Induced Fluorescence - Fluorescence Assay by Gas Expansion (LIF-FAGE) techniques as part of the Indiana Radical, Reactivity and Ozone Production Intercomparison (IRRONIC). This campaign took place in a forested area near the Indiana University, Bloomington campus characterized by high mixing ratios of isoprene and low mixing ratios of NOx. Supporting measurements of photolysis rates, VOCs, NOx, and other species were used to constrain a zero-dimensional box model based on the Regional Atmospheric Chemistry Mechanism (RACM2) and the Master Chemical Mechanism (MCM). Using an OH chemical scavenger technique, the study revealed the presence of an interference with the LIF-FAGE measurements of OH that increased with both ambient concentrations of ozone and temperature. Subtraction of the interference resulted in measured OH concentrations that were in better agreement with model predictions, although the model still underestimated the measured concentrations, likely due to an underestimation of the concentration of NO at this site. Measurements of HO2 radical concentrations during the campaign included a fraction of isoprene-based peroxy radicals (HO2* = HO2 + αRO2) and were found to agree with model predictions. On average, the measured reactivity was consistent with that calculated from measured OH sinks to within 20 %, with modeled oxidation products accounting for the missing reactivity, although significant missing reactivity (approximately 40 % of the total measured reactivity) was observed on some days.


2011 ◽  
Vol 11 (21) ◽  
pp. 10779-10790 ◽  
Author(s):  
M. M. Galloway ◽  
A. J. Huisman ◽  
L. D. Yee ◽  
A. W. H. Chan ◽  
C. L. Loza ◽  
...  

Abstract. We present first-generation and total production yields of glyoxal, methylglyoxal, glycolaldehyde, and hydroxyacetone from the oxidation of isoprene, methyl vinyl ketone (MVK), and methacrolein (MACR) with OH under high NOx conditions. Several of these first-generation yields are not included in commonly used chemical mechanisms, such as the Leeds Master Chemical Mechanism (MCM) v. 3.2. The first-generation yield of glyoxal from isoprene was determined to be 2.1 (±0.6)%. Inclusion of first-generation production of glyoxal, glycolaldehyde and hydroxyacetone from isoprene greatly improves performance of an MCM based model during the initial part of the experiments. In order to further improve performance of the MCM based model, higher generation glyoxal production was reduced by lowering the first-generation yield of glyoxal from C5 hydroxycarbonyls. The results suggest that glyoxal production from reaction of OH with isoprene under high NOx conditions can be approximated by inclusion of a first-generation production term together with secondary production only via glycolaldehyde. Analogously, methylglyoxal production can be approximated by a first-generation production term from isoprene, and secondary production via MVK, MACR and hydroxyacetone. The first-generation yields reported here correspond to less than 5% of the total oxidized yield from isoprene and thus only have a small effect on the fate of isoprene. However, due to the abundance of isoprene, the combination of first-generation yields and reduced higher generation production of glyoxal from C5 hydroxycarbonyls is important for models that include the production of the small organic molecules from isoprene.


2006 ◽  
Vol 6 (8) ◽  
pp. 2193-2214 ◽  
Author(s):  
Z. L. Fleming ◽  
P. S. Monks ◽  
A. R. Rickard ◽  
D. E. Heard ◽  
W. J. Bloss ◽  
...  

Abstract. Peroxy radical (HO2+ΣRO2) measurements, using the PEroxy Radical Chemical Amplification (PERCA) technique at the North Atlantic Marine Boundary Layer EXperiment (NAMBLEX) at Mace Head in summer 2002, are presented and put into the context of marine, boundary-layer chemistry. A suite of other chemical parameters (NO, NO2, NO3, CO, CH4, O3, VOCs, peroxides), photolysis frequencies and meteorological measurements, are used to present a detailed analysis of the role of peroxy radicals in tropospheric oxidation cycles and ozone formation. Under the range of conditions encountered the peroxy radical daily maxima varied from 10 to 40 pptv. The diurnal cycles showed an asymmetric shape typically shifted to the afternoon. Using a box model based on the master chemical mechanism the average model measurement agreement was 2.5 across the campaign. The addition of halogen oxides to the model increases the level of model/measurement agreement, apparently by respeciation of HOx. A good correlation exists between j(HCHO).[HCHO] and the peroxy radicals indicative of the importance of HCHO in the remote atmosphere as a HOx source, particularly in the afternoon. The peroxy radicals showed a strong dependence on [NO2] with a break point at 0.1 ppbv, where the radicals increased concomitantly with the reactive VOC loading, this is a lower value than seen at representative urban campaigns. The HO2/(HO2+ΣRO2) ratios are dependent on [NOx] ranging between 0.2 and 0.6, with the ratio increasing linearly with NOx. Significant night-time levels of peroxy radicals were measured up to 25 pptv. The contribution of ozone-alkenes and NO3-alkene chemistry to night-time peroxy radical production was shown to be on average 59 and 41%. The campaign mean net ozone production rate was 0.11±0.3 ppbv h-1. The ozone production rate was strongly dependent on [NO] having linear sensitivity (dln(P(O3))/dln(NO)=1.0). The results imply that the N(O3) (the in-situ net photochemical rate of ozone production/destruction) will be strongly sensitive in the marine boundary layer to small changes in [NO] which has ramifications for changing NOx loadings in the European continental boundary layer.


2015 ◽  
Vol 15 (21) ◽  
pp. 31247-31286
Author(s):  
L. K. Whalley ◽  
D. Stone ◽  
B. Bandy ◽  
R. Dunmore ◽  
J. F. Hamilton ◽  
...  

Abstract. Near-continuous measurements of OH reactivity in the urban background atmosphere of central London during the summer of 2012 are presented. OH reactivity behaviour is seen to be broadly dependent on airmass origin with the highest reactivity and the most pronounced diurnal profile observed when air had passed over central London to the East, prior to measurement. Averaged over the entire observation period of 26 days, OH reactivity peaked at ~ 27 s−1 in the morning with a minimum of ~ 15 s−1 during the afternoon. A maximum OH reactivity of 116 s−1 was recorded on one day during morning rush hour. A detailed box model using the Master Chemical Mechanism was used to calculate OH reactivity, and was constrained with an extended measurement dataset of volatile organic compounds (VOCs) derived from GC-FID and a two-dimensional GC instrument which included heavier molecular weight (up to C12) aliphatic VOCs, oxygenated VOCs and the biogenic VOCs of α pinene and limonene. Comparison was made between observed OH reactivity and modelled OH reactivity using (i) a standard suite of VOC measurements (C2-C8 hydrocarbons and a small selection of oxygenated VOCs) and (ii) a more comprehensive inventory including species up to C12. Modelled reactivities were lower than those measured (by 33 %) when only the reactivity of the standard VOC suite was considered. The difference between measured and modelled reactivity was improved, to within 15 %, if the reactivity of the higher VOCs (≥ C9) was also considered, with the reactivity of the biogenic compounds of α pinene and limonene and their oxidation products almost entirely responsible for this improvement. Further improvements in the model's ability to reproduce OH reactivity (to within 6 %) could be achieved if the reactivity and degradation mechanism of unassigned two-dimensional GC peaks were estimated. Neglecting the contribution of the higher VOCs (≥ C9) (particularly α pinene and limonene) and model-generated intermediates worsened the agreement between modelled and observed OH concentrations (by 41 %) and the magnitude of in situ ozone production calculated from the production of RO2 was significantly lower (60 %). This work highlights that any future ozone abatement strategies should consider the role that biogenic emissions play alongside anthropogenic emissions in influencing London's air quality.


2021 ◽  
Vol 21 (3) ◽  
pp. 2125-2147
Author(s):  
Lisa K. Whalley ◽  
Eloise J. Slater ◽  
Robert Woodward-Massey ◽  
Chunxiang Ye ◽  
James D. Lee ◽  
...  

Abstract. Measurements of OH, HO2, complex RO2 (alkene- and aromatic-related RO2) and total RO2 radicals taken during the integrated Study of AIR Pollution PROcesses in Beijing (AIRPRO) campaign in central Beijing in the summer of 2017, alongside observations of OH reactivity, are presented. The concentrations of radicals were elevated, with OH reaching up to 2.8×107moleculecm-3, HO2 peaking at 1×109moleculecm-3 and the total RO2 concentration reaching 5.5×109moleculecm-3. OH reactivity (k(OH)) peaked at 89 s−1 during the night, with a minimum during the afternoon of ≈22s-1 on average. An experimental budget analysis, in which the rates of production and destruction of the radicals are compared, highlighted that although the sources and sinks of OH were balanced under high NO concentrations, the OH sinks exceeded the known sources (by 15 ppbv h−1) under the very low NO conditions (<0.5 ppbv) experienced in the afternoons, demonstrating a missing OH source consistent with previous studies under high volatile organic compound (VOC) emissions and low NO loadings. Under the highest NO mixing ratios (104 ppbv), the HO2 production rate exceeded the rate of destruction by ≈50ppbvh-1, whilst the rate of destruction of total RO2 exceeded the production by the same rate, indicating that the net propagation rate of RO2 to HO2 may be substantially slower than assumed. If just 10 % of the RO2 radicals propagate to HO2 upon reaction with NO, the HO2 and RO2 budgets could be closed at high NO, but at low NO this lower RO2 to HO2 propagation rate revealed a missing RO2 sink that was similar in magnitude to the missing OH source. A detailed box model that incorporated the latest Master Chemical Mechanism (MCM3.3.1) reproduced the observed OH concentrations well but over-predicted the observed HO2 under low concentrations of NO (<1 ppbv) and under-predicted RO2 (both the complex RO2 fraction and other RO2 types which we classify as simple RO2) most significantly at the highest NO concentrations. The model also under-predicted the observed k(OH) consistently by ≈10s-1 across all NOx levels, highlighting that the good agreement for OH was fortuitous due to a cancellation of missing OH source and sink terms in its budget. Including heterogeneous loss of HO2 to aerosol surfaces did reduce the modelled HO2 concentrations in line with the observations but only at NO mixing ratios <0.3 ppbv. The inclusion of Cl atoms, formed from the photolysis of nitryl chloride, enhanced the modelled RO2 concentration on several mornings when the Cl atom concentration was calculated to exceed 1×104atomscm-3 and could reconcile the modelled and measured RO2 concentrations at these times. However, on other mornings, when the Cl atom concentration was lower, large under-predictions in total RO2 remained. Furthermore, the inclusion of Cl atom chemistry did not enhance the modelled RO2 beyond the first few hours after sunrise and so was unable to resolve the modelled under-prediction in RO2 observed at other times of the day. Model scenarios, in which missing VOC reactivity was included as an additional reaction that converted OH to RO2, highlighted that the modelled OH, HO2 and RO2 concentrations were sensitive to the choice of RO2 product. The level of modelled to measured agreement for HO2 and RO2 (both complex and simple) could be improved if the missing OH reactivity formed a larger RO2 species that was able to undergo reaction with NO, followed by isomerisation reactions reforming other RO2 species, before eventually generating HO2. In this work an α-pinene-derived RO2 species was used as an example. In this simulation, consistent with the experimental budget analysis, the model underestimated the observed OH, indicating a missing OH source. The model uncertainty, with regards to the types of RO2 species present and the radicals they form upon reaction with NO (HO2 directly or another RO2 species), leads to over an order of magnitude less O3 production calculated from the predicted peroxy radicals than calculated from the observed peroxy radicals at the highest NO concentrations. This demonstrates the rate at which the larger RO2 species propagate to HO2, to another RO2 or indeed to OH needs to be understood to accurately simulate the rate of ozone production in environments such as Beijing, where large multifunctional VOCs are likely present.


2011 ◽  
Vol 11 (4) ◽  
pp. 10693-10720 ◽  
Author(s):  
M. M. Galloway ◽  
A. J. Huisman ◽  
L. D. Yee ◽  
A. W. H. Chan ◽  
C. L. Loza ◽  
...  

Abstract. We present first-generation and total production yields of glyoxal, methylglyoxal, glycolaldehyde, and hydroxyacetone from the oxidation of isoprene, methyl vinyl ketone (MVK), and methacrolein (MACR) with OH under high NOx conditions. Several of these first-generation yields are not included in commonly used chemical mechanisms, such as the Leeds Master Chemical Mechanism (MCM) v. 3.1. Inclusion of first-generation production of glyoxal, glycolaldehyde and hydroxyacetone from isoprene and methylglyoxal from MACR greatly improves performance of an MCM based model during the initial part of the experiments. In order to further improve performance of the MCM based model, higher generation glyoxal production was reduced by lowering the first-generation yield of glyoxal from C5 carbonyls. The results suggest that glyoxal production from reaction of OH with isoprene under high NOx conditions can be approximated by inclusion of a first-generation production term together with secondary production only via glycolaldehyde. Analogously, methylglyoxal production can be approximated by a first-generation production term from isoprene, and secondary production via MVK, MACR and hydroxyacetone. The first-generation yields reported here correspond to less than 5% of the total oxidized yield from isoprene and thus only have a small effect on the fate of isoprene. However, due to the abundance of isoprene, the combination of first-generation yields and reduced higher generation production of glyoxal from C5 carbonyls is important for models which include the production of the small organic molecules from isoprene.


2014 ◽  
Vol 14 (8) ◽  
pp. 3789-3800 ◽  
Author(s):  
T. P. Riedel ◽  
G. M. Wolfe ◽  
K. T. Danas ◽  
J. B. Gilman ◽  
W. C. Kuster ◽  
...  

Abstract. Nitryl chloride (ClNO2) is produced at night by reactions of dinitrogen pentoxide (N2O5) on chloride containing surfaces. ClNO2 is photolyzed during the morning hours after sunrise to liberate highly reactive chlorine atoms (Cl·). This chemistry takes place primarily in polluted environments where the concentrations of N2O5 precursors (nitrogen oxide radicals and ozone) are high, though it likely occurs in remote regions at lower intensities. Recent field measurements have illustrated the potential importance of ClNO2 as a daytime Cl· source and a nighttime NOx reservoir. However, the fate of the Cl· and the overall impact of ClNO2 on regional photochemistry remain poorly constrained by measurements and models. To this end, we have incorporated ClNO2 production, photolysis, and subsequent Cl· reactions into an existing master chemical mechanism (MCM version 3.2) box model framework using observational constraints from the CalNex 2010 field study. Cl· reactions with a set of alkenes and alcohols, and the simplified multiphase chemistry of N2O5, ClNO2, HOCl, ClONO2, and Cl2, none of which are currently part of the MCM, have been added to the mechanism. The presence of ClNO2 produces significant changes to oxidants, ozone, and nitrogen oxide partitioning, relative to model runs excluding ClNO2 formation. From a nighttime maximum of 1.5 ppbv ClNO2, the daytime maximum Cl· concentration reaches 1 × 105 atoms cm−3 at 07:00 model time, reacting mostly with a large suite of volatile organic compounds (VOC) to produce 2.2 times more organic peroxy radicals in the morning than in the absence of ClNO2. In the presence of several ppbv of nitrogen oxide radicals (NOx = NO + NO2), these perturbations lead to similar enhancements in hydrogen oxide radicals (HOx = OH + HO2). Neglecting contributions from HONO, the total integrated daytime radical source is 17% larger when including ClNO2, which leads to a similar enhancement in integrated ozone production of 15%. Detectable levels (tens of pptv) of chlorine containing organic compounds are predicted to form as a result of Cl· addition to alkenes, which may be useful in identifying times of active Cl· chemistry.


Sign in / Sign up

Export Citation Format

Share Document