Simulating the formation of secondary organic aerosol from the photooxidation of aromatic hydrocarbons

The formation and composition of secondary organic aerosol (SOA) from the photooxidation of benzene, p-xylene, and 1,3,5-trimethylbenzene has been simulated using the Master Chemical Mechanism version 3.1 (MCM v3.1) coupled to a representation of the transfer of organic material from the gas to particle phase. The combined mechanism was tested against data obtained from a series of experiments conducted at the European Photoreactor (EUPHORE) outdoor smog chamber in Valencia, Spain. Simulated aerosol mass concentrations compared reasonably well with the measured SOA data only after absorptive partitioning coefficients were increased by a factor of between 5 and 30. The requirement of such scaling was interpreted in terms of the occurrence of unaccounted-for association reactions in the condensed organic phase leading to the production of relatively more nonvolatile species. Comparisons were made between the relative aerosol forming efficiencies of benzene, toluene, p-xylene, and 1,3,5-trimethylbenzene, and differences in the OH-initiated degradation mechanisms of these aromatic hydrocarbons. A strong, nonlinear relationship was observed between measured (reference) yields of SOA and (proportional) yields of unsaturated dicarbonyl aldehyde species resulting from ring-fragmenting pathways. This observation, and the results of the simulations, is strongly suggestive of the involvement of reactive aldehyde species in association reactions occurring in the aerosol phase, thus promoting SOA formation and growth. The effect of NO, concentrations on SOA formation efficiencies (and formation mechanisms) is discussed.

Simulating the detailed chemical composition of secondary organic aerosol formed on a regional scale during the TORCH 2003 campaign in the southern UK

Following on from the companion study (Johnson et al., 2006), a photochemical trajectory model (PTM) has been used to simulate the chemical composition of organic aerosol for selected events during the 2003 TORCH (Tropospheric Organic Chemistry Experiment) field campaign. The PTM incorporates the speciated emissions of 124 nonmethane anthropogenic volatile organic compounds (VOC) and three representative biogenic VOC, a highly-detailed representation of the atmospheric degradation of these VOC, the emission of primary organic aerosol (POA) material and the formation of secondary organic aerosol (SOA) material. SOA formation was represented by the transfer of semi and non-volatile oxidation products from the gas-phase to a condensed organic aerosol-phase, according to estimated thermodynamic equilibrium phase-partitioning characteristics for around 2000 reaction products. After significantly scaling all phase-partitioning coefficients, and assuming a persistent background organic aerosol (both required in order to match the observed organic aerosol loadings), the detailed chemical composition of the simulated SOA has been investigated in terms of intermediate oxygenated species in the Master Chemical Mechanism, version 3.1 ( MCM v3.1). For the various case studies considered, 90% of the simulated SOA mass comprises between ca. 70 and 100 multifunctional oxygenated species derived, in varying amounts, from the photooxidation of VOC of anthropogenic and biogenic origin. The anthropogenic contribution is dominated by aromatic hydrocarbons and the biogenic contribution by alpha-and beta-pinene (which also constitute surrogates for other emitted monoterpene species). Sensitivity in the simulated mass of SOA to changes in the emission rates of anthropogenic and biogenic VOC has also been investigated for 11 case study events, and the results have been compared to the detailed chemical composition data. The role of accretion chemistry in SOA formation, and its implications for the results of the present investigation, is discussed.

Simulating regional scale secondary organic aerosol formation during the TORCH 2003 campaign in the southern UK

A photochemical trajectory model has been used to simulate the chemical evolution of air masses arriving at the TORCH field campaign site in the southern UK during late July and August 2003, a period which included a widespread and prolonged photochemical pollution episode. The model incorporates speciated emissions of 124 nonmethane anthropogenic VOC and three representative biogenic VOC, coupled with a comprehensive description of the chemistry of their degradation. A representation of the gas/aerosol absorptive partitioning of ca. 2000 oxygenated organic species generated in the Master Chemical Mechanism (MCM v3.1) has been implemented, allowing simulation of the contribution to organic aerosol (OA) made by semi- and non-volatile products of VOC oxidation; emissions of primary organic aerosol (POA) and elemental carbon (EC) are also represented. Simulations of total OA mass concentrations in nine case study events (optimised by comparison with observed hourly-mean mass loadings derived from aerosol mass spectrometry measurements) imply that the OA can be ascribed to three general sources: (i) POA emissions; (ii) a '' ubiquitous '' background concentration of 0.7 mu g m(-3); and (iii) gas-to-aerosol transfer of lower volatility products of VOC oxidation generated by the regional scale processing of emitted VOC, but with all partitioning coefficients increased by a species-independent factor of 500. The requirement to scale the partitioning coefficients, and the implied background concentration, are both indicative of the occurrence of chemical processes within the aerosol which allow the oxidised organic species to react by association and/or accretion reactions which generate even lower volatility products, leading to a persistent, non-volatile secondary organic aerosol (SOA). The contribution of secondary organic material to the simulated OA results in significant elevations in the simulated ratio of organic carbon (OC) to EC, compared with the ratio of 1.1 assigned to the emitted components. For the selected case study events, [OC]/[EC] is calculated to lie in the range 2.7-9.8, values which are comparable with the high end of the range reported in the literature.

Oxidation of biogenic and water-soluble compounds in aqueous and organic aerosol droplets by ozone: a kinetic and product analysis approach using laser Raman tweezers

The results of an experimental study into the oxidative degradation of proxies for atmospheric aerosol are presented. We demonstrate that the laser Raman tweezers method can be used successfully to obtain uptake coeffcients for gaseous oxidants on individual aqueous and organic droplets, whilst the size and composition of the droplets is simultaneously followed. A laser tweezers system was used to trap individual droplets containing an unsaturated organic compound in either an aqueous or organic ( alkane) solvent. The droplet was exposed to gas- phase ozone and the reaction kinetics and products followed using Raman spectroscopy. The reactions of three different organic compounds with ozone were studied: fumarate anions, benzoate anions and alpha pinene. The fumarate and benzoate anions in aqueous solution were used to represent components of humic- like substances, HULIS; a alpha- pinene in an alkane solvent was studied as a proxy for biogenic aerosol. The kinetic analysis shows that for these systems the diffusive transport and mass accommodation of ozone is relatively fast, and that liquid- phase di. ffusion and reaction are the rate determining steps. Uptake coe. ffcients, g, were found to be ( 1.1 +/- 0.7) x 10(-5), ( 1.5 +/- 0.7) x 10 (-5) and ( 3.0 - 7.5) x 10 (-3) for the reactions of ozone with the fumarate, benzoate and a- pinene containing droplets, respectively. Liquid- phase bimolecular rate coe. cients for reactions of dissolved ozone molecules with fumarate, benzoate and a- pinene were also obtained: k(fumarate) = ( 2.7 +/- 2) x 10 (5), k(benzoate) = ( 3.5 +/- 3) x 10 (5) and k(alpha-pinene) = ( 1-3) x 10(7) dm(3) mol (-1) s (- 1). The droplet size was found to remain stable over the course of the oxidation process for the HULIS- proxies and for the oxidation of a- pinene in pentadecane. The study of the alpha- pinene/ ozone system is the first using organic seed particles to show that the hygroscopicity of the particle does not increase dramatically over the course of the oxidation. No products were detected by Raman spectroscopy for the reaction of benzoate ions with ozone. One product peak, consistent with aqueous carbonate anions, was observed when following the oxidation of fumarate ions by ozone. Product peaks observed in the reaction of ozone with alpha- pinene suggest the formation of new species containing carbonyl groups.

Mechanisms for the formation of secondary organic aerosol components from the gas-phase ozonolysis of alpha-pinene

Gas-phase ozonolysis of alpha-pinene was studied in static chamber experiments under 'OH-free' conditions. A range of multifunctional products-in particular low-volatility carboxylic acids-were identified in the condensed phase using gas chromatography coupled to mass spectrometry after derivatisation. The dependence of product yields on reaction conditions (humidity, choice of OH radical scavengers, added Criegee intermediate scavengers, NO2 etc.) was investigated to probe the mechanisms of formation of these products; additional information was obtained by studying the ozonolysis of an enal and an enone derived from alpha-pinene. On the basis of experimental findings, previously suggested mechanisms were evaluated and detailed gas-phase mechanisms were developed to explain the observed product formation. Atmospheric implications of this work are discussed.

Chemical ageing and transformation of diffusivity in semi-solid multi-component organic aerosol particles

Recent experimental evidence underlines the importance of reduced diffusivity in amorphous semi-solid or glassy atmospheric aerosols. This paper investigates the impact of diffusivity on the ageing of multi-component reactive organic particles representative of atmospheric cooking aerosols. We apply and extend the recently developed KM-SUB model in a study of a 12-component mixture containing oleic and palmitoleic acids. We demonstrate that changes in the diffusivity may explain the evolution of chemical loss rates in ageing semi-solid particles, and we resolve surface and bulk processes under transient reaction conditions considering diffusivities altered by oligomerisation. This new model treatment allows prediction of the ageing of mixed organic multi-component aerosols over atmospherically relevant time scales and conditions. We illustrate the impact of changing diffusivity on the chemical half-life of reactive components in semisolid particles, and we demonstrate how solidification and crust formation at the particle surface can affect the chemical transformation of organic aerosols.

Chemical ageing and transformation of diffusivity in semi-solid multi-component organic aerosol particles

Recent experimental evidence underlines the importance of reduced diffusivity in amorphous semi-solid or glassy atmospheric aerosols. This paper investigates the impact of diffusivity on the ageing of multi-component reactive organic particles approximating atmospheric cooking aerosols. We apply and extend the recently developed KMSUB model in a study of a 12-component mixture containing oleic and palmitoleic acids. We demonstrate that changes in the diffusivity may explain the evolution of chemical loss rates in ageing semi-solid particles, and we resolve surface and bulk processes under transient reaction conditions considering diffusivities altered by oligomerisation. This new model treatment allows prediction of the ageing of mixed organic multi-component aerosols over atmospherically relevant timescales and conditions. We illustrate the impact of changing diffusivity on the chemical half-life of reactive components in semi-solid particles, and we demonstrate how solidification and crust formation at the particle surface can affect the chemical transformation of organic aerosols.

The AeroCom evaluation and intercomparison of organic aerosol in global models

This paper evaluates the current status of global modeling of the organic aerosol (OA) in the troposphere and analyzes the differences between models as well as between models and observations. Thirty-one global chemistry transport models (CTMs) and general circulation models (GCMs) have participated in this intercomparison, in the framework of AeroCom phase II. The simulation of OA varies greatly between models in terms of the magnitude of primary emissions, secondary OA (SOA) formation, the number of OA species used (2 to 62), the complexity of OA parameterizations (gas-particle partitioning, chemical aging, multiphase chemistry, aerosol microphysics), and the OA physical, chemical and optical properties. The diversity of the global OA simulation results has increased since earlier AeroCom experiments, mainly due to the increasing complexity of the SOA parameterization in models, and the implementation of new, highly uncertain, OA sources. Diversity of over one order of magnitude exists in the modeled vertical distribution of OA concentrations that deserves a dedicated future study. Furthermore, although the OA / OC ratio depends on OA sources and atmospheric processing, and is important for model evaluation against OA and OC observations, it is resolved only by a few global models. The median global primary OA (POA) source strength is 56 Tg a−1 (range 34–144 Tg a−1) and the median SOA source strength (natural and anthropogenic) is 19 Tg a−1 (range 13–121 Tg a−1). Among the models that take into account the semi-volatile SOA nature, the median source is calculated to be 51 Tg a−1 (range 16–121 Tg a−1), much larger than the median value of the models that calculate SOA in a more simplistic way (19 Tg a−1; range 13–20 Tg a−1, with one model at 37 Tg a−1). The median atmospheric burden of OA is 1.4 Tg (24 models in the range of 0.6–2.0 Tg and 4 between 2.0 and 3.8 Tg), with a median OA lifetime of 5.4 days (range 3.8–9.6 days). In models that reported both OA and sulfate burdens, the median value of the OA/sulfate burden ratio is calculated to be 0.77; 13 models calculate a ratio lower than 1, and 9 models higher than 1. For 26 models that reported OA deposition fluxes, the median wet removal is 70 Tg a−1 (range 28–209 Tg a−1), which is on average 85% of the total OA deposition. Fine aerosol organic carbon (OC) and OA observations from continuous monitoring networks and individual field campaigns have been used for model evaluation. At urban locations, the model–observation comparison indicates missing knowledge on anthropogenic OA sources, both strength and seasonality. The combined model–measurements analysis suggests the existence of increased OA levels during summer due to biogenic SOA formation over large areas of the USA that can be of the same order of magnitude as the POA, even at urban locations, and contribute to the measured urban seasonal pattern. Global models are able to simulate the high secondary character of OA observed in the atmosphere as a result of SOA formation and POA aging, although the amount of OA present in the atmosphere remains largely underestimated, with a mean normalized bias (MNB) equal to −0.62 (−0.51) based on the comparison against OC (OA) urban data of all models at the surface, −0.15 (+0.51) when compared with remote measurements, and −0.30 for marine locations with OC data. The mean temporal correlations across all stations are low when compared with OC (OA) measurements: 0.47 (0.52) for urban stations, 0.39 (0.37) for remote stations, and 0.25 for marine stations with OC data. The combination of high (negative) MNB and higher correlation at urban stations when compared with the low MNB and lower correlation at remote sites suggests that knowledge about the processes that govern aerosol processing, transport and removal, on top of their sources, is important at the remote stations. There is no clear change in model skill with increasing model complexity with regard to OC or OA mass concentration. However, the complexity is needed in models in order to distinguish between anthropogenic and natural OA as needed for climate mitigation, and to calculate the impact of OA on climate accurately.

Airborne measurements of the spatial distribution of aerosol chemical composition across Europe and evolution of the organic fraction

The spatial distribution of aerosol chemical composition and the evolution of the Organic Aerosol (OA) fraction is investigated based upon airborne measurements of aerosol chemical composition in the planetary boundary layer across Europe. Sub-micron aerosol chemical composition was measured using a compact Time-of-Flight Aerosol Mass Spectrometer (cToF-AMS). A range of sampling conditions were evaluated, including relatively clean background conditions, polluted conditions in North-Western Europe and the near-field to far-field outflow from such conditions. Ammonium nitrate and OA were found to be the dominant chemical components of the sub-micron aerosol burden, with mass fractions ranging from 20--50% each. Ammonium nitrate was found to dominate in North-Western Europe during episodes of high pollution, reflecting the enhanced NO_x and ammonia sources in this region. OA was ubiquitous across Europe and concentrations generally exceeded sulphate by 30--160%. A factor analysis of the OA burden was performed in order to probe the evolution across this large range of spatial and temporal scales. Two separate Oxygenated Organic Aerosol (OOA) components were identified; one representing an aged-OOA, termed Low Volatility-OOA and another representing fresher-OOA, termed Semi Volatile-OOA on the basis of their mass spectral similarity to previous studies. The factors derived from different flights were not chemically the same but rather reflect the range of OA composition sampled during a particular flight. Significant chemical processing of the OA was observed downwind of major sources in North-Western Europe, with the LV-OOA component becoming increasingly dominant as the distance from source and photochemical processing increased. The measurements suggest that the aging of OA can be viewed as a continuum, with a progression from a less oxidised, semi-volatile component to a highly oxidised, less-volatile component. Substantial amounts of pollution were observed far downwind of continental Europe, with OA and ammonium nitrate being the major constituents of the sub-micron aerosol burden. Such anthropogenically perturbed air masses can significantly perturb regional climate far downwind of major source regions.

The gas-phase ozonolysis of unsaturated volatile organic compounds in the troposphere

The gas-phase reactions of ozone with unsaturated hydrocarbons are significant sources of free radical species (including (OH)-O-center dot) and particulate material in the Earth's atmosphere. In this tutorial review, the kinetics, products and mechanisms of these reactions are examined, starting with a discussion of the original mechanism proposed by Criegee and following with a summary presentation of the complex, free radical-mediated reactions of carbonyl oxide (Criegee) intermediates. The contribution of ozone-terpene reactions to the atmospheric burden of secondary organic aerosol material is also discussed from the viewpoint of the formation of non-volatile organic acid products from the complex chemistry of ozone with alpha-pinene. Throughout the article, currently accepted understanding is supported through the presentation of key experimental results, and areas of persistent or new uncertainty are highlighted.

A study on the sensitivities of simulated aerosol optical properties to composition and size distribution using airborne measurements

We present a flexible framework to calculate the optical properties of atmospheric aerosols at a given relative humidity based on their composition and size distribution. The similarity of this framework to climate model parameterisations allows rapid and extensive sensitivity tests of the impact of uncertainties in data or of new measurements on climate relevant aerosol properties. The data collected by the FAAM BAe-146 aircraft during the EUCAARI-LONGREX and VOCALS-REx campaigns have been used in a closure study to analyse the agreement between calculated and measured aerosol optical properties for two very different aerosol types. The agreement achieved for the EUCAARI-LONGREX flights is within the measurement uncertainties for both scattering and absorption. However, there is poor agreement between the calculated and the measured scattering for the VOCALS-REx flights. The high concentration of sulphate, which is a scattering aerosol with no absorption in the visible spectrum, made the absorption measurements during VOCALS-REx unreliable, and thus no closure study was possible for the absorption. The calculated hygroscopic scattering growth factor overestimates the measured values during EUCAARI-LONGREX and VOCALS-REx by ∼30% and ∼20%, respectively. We have also tested the sensitivity of the calculated aerosol optical properties to the uncertainties in the refractive indices, the hygroscopic growth factors and the aerosol size distribution. The largest source of uncertainty in the calculated scattering is the aerosol size distribution (∼35%), followed by the assumed hygroscopic growth factor for organic aerosol (∼15%), while the predominant source of uncertainty in the calculated absorption is the refractive index of organic aerosol (28–60%), although we would expect the refractive index of black carbon to be important for aerosol with a higher black carbon fraction.

Aerosol scattering and absorption during the EUCAARI-LONGREX flights of the Facility for Airborne Atmospheric Measurements (FAAM) BAe-146: can measurements and models agree?

Scattering and absorption by aerosol in anthropogenically perturbed air masses over Europe has been measured using instrumentation flown on the UK’s BAe-146-301 large Atmospheric Research Aircraft (ARA) operated by the Facility for Airborne Atmospheric Measurements (FAAM) on 14 flights during the EUCAARI-LONGREX campaign in May 2008. The geographical and temporal variations of the derived shortwave optical properties of aerosol are presented. Values of single scattering albedo of dry aerosol at 550 nm varied considerably from 0.86 to near unity, with a campaign average of 0.93 ± 0.03. Dry aerosol optical depths ranged from 0.030 ± 0.009 to 0.24 ± 0.07. An optical properties closure study comparing calculations from composition data and Mie scattering code with the measured properties is presented. Agreement to within measurement uncertainties of 30% can be achieved for both scattering and absorption,but the latter is shown to be sensitive to the refractive indices chosen for organic aerosols, and to a lesser extent black carbon, as well as being highly dependent on the accuracy of the absorption measurements. Agreement with the measured absorption can be achieved either if organic carbon is assumed to be weakly absorbing, or if the organic aerosol is purely scattering and the absorption measurement is an overestimate due to the presence of large amounts of organic carbon. Refractive indices could not be inferred conclusively due to this uncertainty, despite the enhancement in methodology compared to previous studies that derived from the use of the black carbon measurements. Hygroscopic growth curves derived from the wet nephelometer indicate moderate water uptake by the aerosol with a campaign mean f (RH) value (ratio in scattering) of 1.5 (range from 1.23 to 1.63) at 80% relative humidity. This value is qualitatively consistent with the major chemical components of the aerosol measured by the aerosol mass spectrometer, which are primarily mixed organics and nitrate and some sulphate.

Kinetic multi-layer model of gas-particle interactions in aerosols and clouds (KM-GAP): linking condensation, evaporation and chemical reactions of organics, oxidants and water

We present a novel kinetic multi-layer model for gas-particle interactions in aerosols and clouds (KM-GAP) that treats explicitly all steps of mass transport and chemical reaction of semi-volatile species partitioning between gas phase, particle surface and particle bulk. KM-GAP is based on the PRA model framework (Pöschl-Rudich-Ammann, 2007), and it includes gas phase diffusion, reversible adsorption, surface reactions, bulk diffusion and reaction, as well as condensation, evaporation and heat transfer. The size change of atmospheric particles and the temporal evolution and spatial profile of the concentration of individual chemical species can be modelled along with gas uptake and accommodation coefficients. Depending on the complexity of the investigated system, unlimited numbers of semi-volatile species, chemical reactions, and physical processes can be treated, and the model shall help to bridge gaps in the understanding and quantification of multiphase chemistry and microphysics in atmo- spheric aerosols and clouds. In this study we demonstrate how KM-GAP can be used to analyze, interpret and design experimental investigations of changes in particle size and chemical composition in response to condensation, evaporation, and chemical reaction. For the condensational growth of water droplets, our kinetic model results provide a direct link between laboratory observations and molecular dynamic simulations, confirming that the accommodation coefficient of water at 270 K is close to unity. Literature data on the evaporation of dioctyl phthalate as a function of particle size and time can be reproduced, and the model results suggest that changes in the experimental conditions like aerosol particle concentration and chamber geometry may influence the evaporation kinetics and can be optimized for eðcient probing of specific physical effects and parameters. With regard to oxidative aging of organic aerosol particles, we illustrate how the formation and evaporation of volatile reaction products like nonanal can cause a decrease in the size of oleic acid particles exposed to ozone.

Kinetic multi-layer model of gas-particle interactions in aerosols and clouds (KM-GAP): linking condensation, evaporation and chemical reactions of organics, oxidants and water

We present a novel kinetic multi-layer model for gas-particle interactions in aerosols and clouds (KMGAP) that treats explicitly all steps of mass transport and chemical reaction of semi-volatile species partitioning between gas phase, particle surface and particle bulk. KMGAP is based on the PRA model framework (P¨oschl-Rudich- Ammann, 2007), and it includes gas phase diffusion, reversible adsorption, surface reactions, bulk diffusion and reaction, as well as condensation, evaporation and heat transfer. The size change of atmospheric particles and the temporal evolution and spatial profile of the concentration of individual chemical species can be modeled along with gas uptake and accommodation coefficients. Depending on the complexity of the investigated system and the computational constraints, unlimited numbers of semi-volatile species, chemical reactions, and physical processes can be treated, and the model shall help to bridge gaps in the understanding and quantification of multiphase chemistry and microphysics in atmospheric aerosols and clouds. In this study we demonstrate how KM-GAP can be used to analyze, interpret and design experimental investigations of changes in particle size and chemical composition in response to condensation, evaporation, and chemical reaction. For the condensational growth of water droplets, our kinetic model results provide a direct link between laboratory observations and molecular dynamic simulations, confirming that the accommodation coefficient of water at 270K is close to unity (Winkler et al., 2006). Literature data on the evaporation of dioctyl phthalate as a function of particle size and time can be reproduced, and the model results suggest that changes in the experimental conditions like aerosol particle concentration and chamber geometry may influence the evaporation kinetics and can be optimized for efficient probing of specific physical effects and parameters. With regard to oxidative aging of organic aerosol particles, we illustrate how the formation and evaporation of volatile reaction products like nonanal can cause a decrease in the size of oleic acid particles exposed to ozone.

Airborne measurements of trace gases and aerosols over the London metropolitan region

The Emissions around the M25 motorway (EM25) campaign took place over the megacity of London in the United Kingdom in June 2009 with the aim of characterising trace gas and aerosol composition and properties entering and emitted from the urban region. It featured two mobile platforms, the UK BAe-146 Facility for Airborne Atmospheric Measurements (FAAM) research aircraft and a ground-based mobile lidar van, both travelling in circuits around London, roughly following the path of the M25 motorway circling the city. We present an overview of findings from the project, which took place during typical UK summertime pollution conditions. Emission ratios of volatile organic compounds (VOCs) to acetylene and carbon monoxide emitted from the London region were consistent with measurements in and downwind of other large urban areas and indicated traffic and associated fuel evaporation were major sources. Sub-micron aerosol composition was dominated by secondary species including sulphate (24% of sub-micron mass in the London plume and 29% in the non-plume regional aerosol), nitrate (24% plume; 20% regional) and organic aerosol (29% plume; 31% regional). The primary sub-micron aerosol emissions from London were minor compared to the larger regional background, with only limited increases in aerosol mass in the urban plume compared to the background (~12% mass increase on average). Black carbon mass was the major exception and more than doubled in the urban plume, leading to a decrease in the single scattering albedo from 0.91 in the regional aerosol to 0.86 in the London plume, on average. Our observations indicated that regional aerosol plays a major role on aerosol concentrations around London, at least during typical summertime conditions, meaning future efforts to reduce PM levels in London must account for regional as well as local aerosol sources.