The SOA/VOC/NOx system: an explicit model of secondary organic aerosol formation

International audienceOur current understanding of secondary organic aerosol (SOA) formation is limited by our knowledge of gaseous secondary organics involved in gas/particle partitioning. The objective of this study is to explore (i) the potential for products of multiple oxidation steps contributing to SOA, and (ii) the evolution of the SOA/VOC/NOx system. We developed an explicit model based on the coupling of detailed gas-phase oxidation schemes with a thermodynamic condensation module. Such a model allows prediction of SOA mass and speciation on the basis of first principles. The SOA/VOC/NOx system is studied for the oxidation of 1-octene under atmospherically relevant concentrations. In this study, gaseous oxidation of octene is simulated to lead to SOA formation. Contributors to SOA formation are shown to be formed via multiple oxidation steps of the parent hydrocarbon. The behaviour of the SOA/VOC/NOx system simulated using the explicit model agrees with general tendencies observed during laboratory chamber experiments. This explicit modelling of SOA formation appears as a useful exploratory tool to (i) support interpretations of SOA formation observed in laboratory chamber experiments, (ii) give some insights on SOA formation under atmospherically relevant conditions and (iii) investigate implications for the regional/global lifetimes of the SOA

Modeling SOA formation from the oxidation of intermediate volatility <i>n</i>-alkanes

The chemical mechanism leading to SOA formation and ageing is expected to be a multigenerational process, i.e. a successive formation of organic compounds with higher oxidation degree and lower vapor pressure. This process is here investigated with the explicit oxidation model GECKO-A (Generator of Explicit Chemistry and Kinetics of Organics in the Atmosphere). Gas phase oxidation schemes are generated for the C₈–C₂₄ series of <i>n</i>-alkanes. Simulations are conducted to explore the time evolution of organic compounds and the behavior of secondary organic aerosol (SOA) formation for various preexisting organic aerosol concentration (<i>C</i>_OA). As expected, simulation results show that (i) SOA yield increases with the carbon chain length of the parent hydrocarbon, (ii) SOA yield decreases with decreasing <i>C</i>_OA, (iii) SOA production rates increase with increasing <i>C</i>_OA and (iv) the number of oxidation steps (i.e. generations) needed to describe SOA formation and evolution grows when <i>C</i>_OA decreases. The simulated oxidative trajectories are examined in a two dimensional space defined by the mean carbon oxidation state and the volatility. Most SOA contributors are not oxidized enough to be categorized as highly oxygenated organic aerosols (OOA) but reduced enough to be categorized as hydrocarbon like organic aerosols (HOA), suggesting that OOA may underestimate SOA. Results show that the model is unable to produce highly oxygenated aerosols (OOA) with large yields. The limitations of the model are discussed

SCIAMACHY formaldehyde observations: constraint for isoprene emission estimates over Europe?

Formaldehyde (HCHO) is an important intermediate compound in the degradation of volatile organic compounds (VOCs) in the troposphere. Sources of HCHO are largely dominated by its secondary production from VOC oxidation, methane and isoprene being the main precursors in unpolluted areas. As a result of the moderate lifetime of HCHO, its spatial distribution is determined by reactive hydrocarbon emissions. We focus here on Europe and investigate the influence of the different emissions on HCHO tropospheric columns with the CHIMERE chemical transport model in order to interpret the comparisons between SCIAMACHY and simulated HCHO columns. Europe was never specifically studied before for these purposes using satellite observations. The bias between measurements and model is less than 20% on average. The differences are discussed according to the errors on the model and the observations and remaining discrepancies are attributed to a misrepresentation of biogenic emissions. This study requires the characterisation of: (1) the model errors and performances concerning formaldehyde. The errors on the HCHO columns, mainly related to chemistry and mixed emission types, are evaluated to 2&amp;times;10&lt;sup&gt;15&lt;/sup&gt; molecule/cm&lt;sup&gt;2&lt;/sup&gt; and the model performances evaluated using surface measurements are satisfactory (~13%); (2) the observation errors that define the needs in spatial and temporal averaging for meaningful comparisons. Using SCIAMACHY observations as constraint for biogenic isoprene emissions in an inverse modelling scheme reduces their uncertainties by about a factor of two in region of intense emissions. The retrieved correction factors for the isoprene emissions range from a factor of 0.15 (North Africa) to a factor of 2 (Poland, the United Kingdom) depending on the regions

Explicit modelling of SOA formation from α-pinene photooxidation: sensitivity to vapour pressure estimation

The sensitivity of the formation of secondary organic aerosol (SOA) to the estimated vapour pressures of the condensable oxidation products is explored. A highly detailed reaction scheme was generated for α-pinene photooxidation using the Generator for Explicit Chemistry and Kinetics of Organics in the Atmosphere (GECKO-A). Vapour pressures (P^(vap)) were estimated with three commonly used structure activity relationships. The values of P^(vap) were compared for the set of secondary species generated by GECKO-A to describe α-pinene oxidation. Discrepancies in the predicted vapour pressures were found to increase with the number of functional groups borne by the species. For semi-volatile organic compounds (i.e. organic species of interest for SOA formation), differences in the predicted Pvap range between a factor of 5 to 200 on average. The simulated SOA concentrations were compared to SOA observations in the Caltech chamber during three experiments performed under a range of NO_x conditions. While the model captures the qualitative features of SOA formation for the chamber experiments, SOA concentrations are systematically overestimated. For the conditions simulated, the modelled SOA speciation appears to be rather insensitive to the P^vap estimation method

Explicit modeling of volatile organic compounds partitioning in the atmospheric aqueous phase

The gas phase oxidation of organic species is a multigenerational process involving a large number of secondary compounds. Most secondary organic species are water-soluble multifunctional oxygenated molecules. The fully explicit chemical mechanism GECKO-A (Generator of Explicit Chemistry and Kinetics of Organics in the Atmosphere) is used to describe the oxidation of organics in the gas phase and their mass transfer to the aqueous phase. The oxidation of three hydrocarbons of atmospheric interest (isoprene, octane and α-pinene) is investigated for various NOx conditions. The simulated oxidative trajectories are examined in a new two dimensional space defined by the mean oxidation state and the solubility. The amount of dissolved organic matter was found to be very low (yield less than 2% on carbon atom basis) under a water content typical of deliquescent aerosols. For cloud water content, 50% (isoprene oxidation) to 70% (octane oxidation) of the carbon atoms are found in the aqueous phase after the removal of the parent hydrocarbons for low NOx conditions. For high NOx conditions, this ratio is only 5% in the isoprene oxidation case, but remains large for α-pinene and octane oxidation cases (40% and 60%, respectively). Although the model does not yet include chemical reactions in the aqueous phase, much of this dissolved organic matter should be processed in cloud drops and modify both oxidation rates and the speciation of organic species

Explicit modeling of volatile organic compounds partitioning in the atmospheric aqueous phase

The gas phase oxidation of organic species is a multigenerational process involving a large number of secondary compounds. Most secondary organic species are water-soluble multifunctional oxygenated molecules. The fully explicit chemical mechanism GECKO-A (Generator of Explicit Chemistry and Kinetics of Organics in the Atmosphere) is used to describe the oxidation of organics in the gas phase and their mass transfer to the aqueous phase. The oxidation of three hydrocarbons of atmospheric interest (isoprene, octane and α-pinene) is investigated for various NOx conditions. The simulated oxidative trajectories are examined in a new two dimensional space defined by the mean oxidation state and the solubility. The amount of dissolved organic matter was found to be very low (yield less than 2% on carbon atom basis) under a water content typical of deliquescent aerosols. For cloud water content, 50% (isoprene oxidation) to 70% (octane oxidation) of the carbon atoms are found in the aqueous phase after the removal of the parent hydrocarbons for low NOx conditions. For high NOx conditions, this ratio is only 5% in the isoprene oxidation case, but remains large for α-pinene and octane oxidation cases (40% and 60%, respectively). Although the model does not yet include chemical reactions in the aqueous phase, much of this dissolved organic matter should be processed in cloud drops and modify both oxidation rates and the speciation of organic species

Explicit modeling of organic chemistry and secondary organic aerosol partitioning for Mexico City and its outflow plume

The evolution of organic aerosols (OA) in Mexico City and its outflow is investigated with the nearly explicit gas phase photochemistry model GECKO-A (Generator of Explicit Chemistry and Kinetics of Organics in the Atmosphere), wherein precursor hydrocarbons are oxidized to numerous intermediate species for which vapor pressures are computed and used to determine gas/particle partitioning in a chemical box model. Precursor emissions included observed C3-10 alkanes, alkenes, and light aromatics, as well as larger &lt;i&gt;n&lt;/i&gt;-alkanes (up to C25) not directly observed but estimated by scaling to particulate emissions according to their volatility. Conditions were selected for comparison with observations made in March 2006 (MILAGRO). The model successfully reproduces the magnitude and diurnal shape for both primary (POA) and secondary (SOA) organic aerosols, with POA peaking in the early morning at 15–20 μg m&lt;sup&gt;−3&lt;/sup&gt;, and SOA peaking at 10–15 μg m&lt;sup&gt;−3&lt;/sup&gt; during mid-day. The majority (&amp;ge;75%) of the model SOA stems from reaction products of the large &lt;i&gt;n&lt;/i&gt;-alkanes, used here as surrogates for all emitted hydrocarbons of similar volatility, with the remaining SOA originating mostly from the light aromatics. Simulated OA elemental composition reproduces observed H/C and O/C ratios reasonably well, although modeled ratios develop more slowly than observations suggest. SOA chemical composition is initially dominated by &amp;delta;-hydroxy ketones and nitrates from the large alkanes, with contributions from peroxy acyl nitrates and, at later times when NOx is lower, organic hydroperoxides. The simulated plume-integrated OA mass continues to increase for several days downwind despite dilution-induced particle evaporation, since oxidation chemistry leading to SOA formation remains strong. In this model, the plume SOA burden several days downwind exceeds that leaving the city by a factor of &gt;3. These results suggest significant regional radiative impacts of SOA

Study of the unknown HONO daytime source at a European suburban site during the MEGAPOLI summer and winter field campaigns

International audienceNitrous acid measurements were carried out during the MEGAPOLI summer and winter field campaigns at SIRTA observatory in Paris surroundings. Highly variable HONO levels were observed during the campaigns, ranging from 10 ppt to 500 ppt in summer and from 10 ppt to 1.7 ppb in winter. Significant HONO mixing ratios have also been measured during daytime hours, comprised between some tenth of ppt and 200 ppt for the summer campaign and between few ppt and 1 ppb for the winter campaign. Ancillary measurements, such as NOx , O3 , photolysis frequencies, meteorological parameters (pressure, temperature, relative humidity , wind speed and wind direction), black carbon concentration , total aerosol surface area, boundary layer height and soil moisture, were conducted during both campaigns. In addition, for the summer period, OH radical measurements were made with a CIMS (Chemical Ionisation Mass Spectrometer). This large dataset has been used to investigate the HONO budget in a suburban environment. To do so, calculations of HONO concentrations using PhotoStationary State (PSS) approach have been performed, for daytime hours. The comparison of these calculations with measured HONO concentrations revealed an underestimation of the calculations making evident a missing source term for both campaigns. This unknown HONO source exhibits a bell-shaped like average diurnal profile with a maximum around noon of approximately 0.7 ppb h−1 and 0.25 ppb h−1 , during summer and winter respectively. This source is the main HONO source during daytime hours for both campaigns. In both cases, this source shows a slight positive correlation with J (NO2) and the product between J (NO2) and soil moisture. This original approach had, thus, indicated that this missing source is photolytic and might be heterogeneous occurring at ground surface and involving water content available on the ground. Published by Copernicus Publications on behalf of the European Geosciences Union. 2806 V. Michoud et al.: Study of the unknown HONO daytime sourc

Secondary organic aerosol formation and composition from the photo-oxidation of methyl chavicol (estragole)

The increasing demand for palm oil for uses in biofuel and food products is leading to rapid expansion of oil palm agriculture. Methyl chavicol (also known as estragole and 1-allyl-4-methoxybenzene) is an oxygenated biogenic volatile organic compound (VOC) that was recently identified as the main floral emission from an oil palm plantation in Malaysian Borneo. The emissions of methyl chavicol observed may impact regional atmospheric chemistry, but little is known of its ability to form secondary organic aerosol (SOA). The photo-oxidation of methyl chavicol was investigated at the European Photoreactor chamber as a part of the atmospheric chemistry of methyl chavicol (ATMECH) project. Aerosol samples were collected using a particle into liquid sampler (PILS) and analysed offline using an extensive range of instruments including; high-performance liquid chromatography mass spectrometry (HPLC-ITMS), high-performance liquid chromatography quadrupole time-of-flight mass spectrometry (HPLC-QTOFMS) and Fourier transform ion cyclotron resonance mass spectrometry (FTICR-MS). The SOA yield was determined as 18 and 29% for an initial VOC mixing ratio of 212 and 460 ppbv (parts per billion by volume) respectively; using a VOC:NOx ratio of ~5:1. In total, 59 SOA compounds were observed and the structures of 10 compounds have been identified using high-resolution tandem mass spectrometry. The addition of hydroxyl and/or nitro-functional groups to the aromatic ring appears to be an important mechanistic pathway for aerosol formation. This results in the formation of compounds with both low volatility and high O:C ratios, where functionalisation rather than fragmentation is mainly observed as a result of the stability of the ring. The SOA species observed can be characterised as semi-volatile to low-volatility oxygenated organic aerosol (SVOOA and LVOOA) components and therefore may be important in aerosol formation and growth

Exploration of the influence of environmental conditions on secondary organic aerosol formation and organic species properties using explicit simulations: development of the VBS-GECKO parameterization

Atmospheric chambers have been widely used to study secondary organic aerosol (SOA) properties and formation from various precursors under different controlled environmental conditions and to develop parameterization to represent SOA formation in chemical transport models (CTMs). Chamber experiments are however limited in number, performed under conditions that differ from the atmosphere and can be subject to potential artefacts from chamber walls. Here, the Generator for Explicit Chemistry and Kinetics of Organics in the Atmosphere (GECKO-A) modelling tool has been used in a box model under various environmental conditions to (i) explore the sensitivity of SOA formation and properties to changes on physical and chemical conditions and (ii) develop a volatility basis set (VBS)-type parameterization. The set of parent hydrocarbons includes n-alkanes and 1-alkenes with 10, 14, 18, 22 and 26 carbon atoms, α-pinene, β-pinene and limonene, benzene, toluene, o-xylene, m-xylene and p-xylene. Simulated SOA yields and their dependences on the precursor structure, organic aerosol load, temperature and NOx levels are consistent with the literature. GECKO-A was used to explore the distribution of molar mass, vaporization enthalpy, OH reaction rate and Henry's law coefficient of the millions of secondary organic compounds formed during the oxidation of the different precursors and under various conditions. From these explicit simulations, a VBS-GECKO parameterization designed to be implemented in 3-D air quality models has been tuned to represent SOA formation from the 18 precursors using GECKO-A as a reference. In evaluating the ability of VBS-GECKO to capture the temporal evolution of SOA mass, the mean relative error is less than 20&thinsp;% compared to GECKO-A. The optimization procedure has been automated to facilitate the update of the VBS-GECKO on the basis of the future GECKO-A versions, its extension to other precursors and/or its modification to carry additional information.</p