Modelling the formation and composition of secondary organic aerosol from ?- and ?-pinene ozonolysis using MCM v3

International audienceThe formation and detailed composition of secondary organic aerosol (SOA) from the gas phase ozonolysis of ?- and ?-pinene has been simulated using the Master Chemical Mechanism version 3 (MCM v3), coupled with a representation of gas-to-aerosol transfer of semivolatile and involatile oxygenated products. A kinetics representation, based on equilibrium absorptive partitioning of ca. 200 semivolatile products, was found to provide an acceptable description of the final mass concentrations observed in a number of reported laboratory and chamber experiments, provided partitioning coefficients were increased by about two orders of magnitude over those defined on the basis of estimated vapour pressures. This adjustment is believed to be due, at least partially, to the effect of condensed phase association reactions of the partitioning products. Even with this adjustment, the simulated initial formation of SOA was delayed relative to that observed, implying the requirement for the formation of species of much lower volatility to initiate SOA formation. The inclusion of a simplified representation of the formation and gas-to-aerosol transfer of involatile dimers of 22 bi- and multifunctional carboxylic acids (in addition to the absorptive partitioning mechanism) allowed a much improved description of SOA formation for a wide range of conditions. The simulated SOA composition recreates certain features of the product distributions observed in a number of experimental studies, but implies an important role for multifunctional products containing hydroperoxy groups (i.e. hydroperoxides). This is particularly the case for experiments in which 2-butanol is used to scavenge OH radicals, because [HO2]/[RO2] ratios are elevated in such systems. The optimized mechanism is used to calculate SOA yields from ?- and ?-pinene ozonolysis in the presence and absence of OH scavengers, and as a function of temperature

Modelling the formation and composition of secondary organic aerosol from α- and β-pinene ozonolysis using MCM v3

The formation and detailed composition of secondary organic aerosol (SOA) from the gas phase ozonolysis of α- and β-pinene has been simulated using the Master Chemical Mechanism version 3 (MCM v3), coupled with a representation of gas-to-aerosol transfer of semivolatile and involatile oxygenated products. A kinetics representation, based on equilibrium absorptive partitioning of ca. 200 semivolatile products, was found to provide an acceptable description of the final mass concentrations observed in a number of reported laboratory and chamber experiments, provided partitioning coefficients were increased by about two orders of magnitude over those defined on the basis of estimated vapour pressures. This adjustment is believed to be due, at least partially, to the effect of condensed phase association reactions of the partitioning products. Even with this adjustment, the simulated initial formation of SOA was delayed relative to that observed, implying the requirement for the formation of species of much lower volatility to initiate SOA formation. The inclusion of a simplified representation of the formation and gas-to-aerosol transfer of involatile dimers of 22 bi- and multifunctional carboxylic acids (in addition to the absorptive partitioning mechanism) allowed a much improved description of SOA formation for a wide range of conditions. The simulated SOA composition recreates certain features of the product distributions observed in a number of experimental studies, but implies an important role for multifunctional products containing hydroperoxy groups (i.e. hydroperoxides). This is particularly the case for experiments in which 2-butanol is used to scavenge OH radicals, because [HO₂]/[RO₂] ratios are elevated in such systems. The optimized mechanism is used to calculate SOA yields from α- and β-pinene ozonolysis in the presence and absence of OH scavengers, and as a function of temperature

Sensitivities of the absorptive partitioning model of secondary organic aerosol formation to the inclusion of water

The predicted distribution of semi-volatile organic components between the gaseous and condensed phase as a function of ambient relative humidity and the specific form of the partitioning model used has been investigated. A mole fraction based model, modified so as not to use molar mass in the calculation, was found to predict identical RH dependence of component partitioning to that predicted by the conventional mass-based partitioning model which uses a molar mass averaged according to the number of moles in the condensed phase. A recently reported third version of the partitioning model using individual component molar masses was shown to give significantly different results to the other two models. Further sensitivities to an assumed pre-existing particulate loading and to parameterised organic component non-ideality are explored and shown to contribute significantly to the variation in predicted secondary organic particulate loading

Protocol for the development of the Master Chemical Mechanism, MCM v3 (Part A): tropospheric degradation of non-aromatic volatile organic compounds

Kinetic and mechanistic data relevant to the tropospheric degradation of volatile organic compounds (VOC), and the production of secondary pollutants, have previously been used to define a protocol which underpinned the construction of a near-explicit Master Chemical Mechanism. In this paper, an update to the previous protocol is presented, which has been used to define degradation schemes for 107 non-aromatic VOC as part of version 3 of the Master Chemical Mechanism (MCM v3). The treatment of 18 aromatic VOC is described in a companion paper. The protocol is divided into a series of subsections describing initiation reactions, the reactions of the radical intermediates and the further degradation of first and subsequent generation products. Emphasis is placed on updating the previous information, and outlining the methodology which is specifically applicable to VOC not considered previously (e.g. <font face='Symbol' >a</font>- and <font face='Symbol' >b</font>-pinene). The present protocol aims to take into consideration work available in the open literature up to the beginning of 2001, and some other studies known by the authors which were under review at the time. Application of MCM v3 in appropriate box models indicates that the representation of isoprene degradation provides a good description of the speciated distribution of oxygenated organic products observed in reported field studies where isoprene was the dominant emitted hydrocarbon, and that the <font face='Symbol' >a</font>-pinene degradation chemistry provides a good description of the time dependence of key gas phase species in <font face='Symbol' >a</font>-pinene/NO_X photo-oxidation experiments carried out in the European Photoreactor (EUPHORE). Photochemical Ozone Creation Potentials (POCP) have been calculated for the 106 non-aromatic non-methane VOC in MCM v3 for idealised conditions appropriate to north-west Europe, using a photochemical trajectory model. The POCP values provide a measure of the relative ozone forming abilities of the VOC. Where applicable, the values are compared with those calculated with previous versions of the MCM

Evaluation of detailed aromatic mechanisms (MCMv3 and MCMv3.1) against environmental chamber data

International audienceA high quality dataset on the photo-oxidation of benzene, toluene, p-xylene and 1,3,5-trimethylbenzene has been obtained from experiments in the European Photoreactor (EUPHORE), a large outdoor environmental reaction chamber. The experiments were designed to test sensitive features of detailed aromatic mechanisms, and the dataset has been used to evaluate the performance of the Master Chemical Mechanism Version 3 (MCMv3). An updated version (MCMv3.1) was constructed based on recent experimental data, and details of its development are described in a companion paper. The MCMv3.1 aromatic mechanisms have also been evaluated using the EUPHORE dataset. Significant deficiencies have been identified in the mechanisms, in particular: 1) an over-estimation of the ozone concentration, 2) an under-estimation of the NO oxidation rate, 3) an under-estimation of OH. The use of MCMv3.1 improves the model-measurement agreement in some areas but significant discrepancies remain

Dilatancy, Jamming, and the Physics of Granulation

Granulation is a process whereby a dense colloidal suspension is converted into pasty granules (surrounded by air) by application of shear. Central to the stability of the granules is the capillary force arising from the interfacial tension between solvent and air. This force appears capable of maintaining a solvent granule in a jammed solid state, under conditions where the same amount of solvent and colloid could also exist as a flowable droplet. We argue that in the early stages of granulation the physics of dilatancy, which requires that a powder expand on shearing, is converted by capillary forces into the physics of arrest. Using a schematic model of colloidal arrest under stress, we speculate upon various jamming and granulation scenarios. Some preliminary experimental results on aspects of granulation in hard-sphere colloidal suspensions are also reported.Comment: Original article intended for J Phys Cond Mat special issue on Granular Materials (M Nicodemi, Ed.

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

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-NO_x 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 NO₂ conversion or 2) NO₃ 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

Estimation of rate coefficients and branching ratios for gas-phase reactions of OH with aromatic organic compounds for use in automated mechanism construction

Reaction with the hydroxyl (OH) radical is the dominant removal process for volatile organic compounds (VOCs) in the atmosphere. Rate coefficients for the reactions of OH with VOCs are therefore essential parameters for chemical mechanisms used in chemistry transport models, and are required more generally for impact assessments involving estimation of atmospheric lifetimes or oxidation rates for VOCs. A structure–activity relationship (SAR) method is presented for the reactions of OH with aromatic organic compounds, with the reactions of aliphatic organic compounds considered in the preceding companion paper. The SAR is optimized using a preferred set of data including reactions of OH with 67 monocyclic aromatic hydrocarbons and oxygenated organic compounds. In each case, the rate coefficient is defined in terms of a summation of partial rate coefficients for H abstraction or OH addition at each relevant site in the given organic compound, so that the attack distribution is defined. The SAR can therefore guide the representation of the OH reactions in the next generation of explicit detailed chemical mechanisms. Rules governing the representation of the reactions of the product radicals under tropospheric conditions are also summarized, specifically the rapid reaction sequences initiated by their reactions with O2

The CRI v2.2 reduced degradation scheme for isoprene

The reduced representation of isoprene degradation in the Common Representative Intermediates (CRI) mechanism has been systematically updated, using the Master Chemical Mechanism (MCM v3.3.1) as a reference benchmark, with the updated mechanism being released as CRI v2.2. The complete isoprene degradation mechanism in CRI v2.2 consists of 186 reactions of 56 closed shell and free radical species, this being an order of magnitude reduction in size compared with MCM v3.3.1. The chemistry initiated by reaction with OH radicals, NO3 radicals and ozone (O3) is treated. An overview of the updates is provided, within the context of reported kinetic and mechanistic information. The revisions mainly relate to the OH-initiated chemistry, which tends to dominate under atmospheric conditions, although these include updates to the chemistry of products that are also generated from the O3- and NO3-initiated oxidation. The revisions have impacts in a number of key areas, including recycling of HOx and NOx. The performance of the CRI v2.2 isoprene mechanism has been compared with those of the preceding version (CRI v2.1) and the reference MCM v3.3.1 over a range of relevant conditions, using a box model of the tropical forested boundary layer. In addition, tests are carried out to ensure that the performance of MCM v3.3.1 remains robust to more recently reported information. CRI v2.2 has also been implemented into the STOCHEM chemistry-transport model, with a customized close-variant of CRI v2.2 implemented into the EMEP MSC-W chemistry-transport model. The results of these studies are presented and used to illustrate the global-scale impacts of the mechanistic updates on HOx radical concentrations