Kinetic multi-layer model of aerosol surface and bulk chemistry (KM-SUB): the influence of interfacial transport and bulk diffusion on the oxidation of oleic acid by ozone

We present a novel kinetic multi-layer model that explicitly resolves mass transport and chemical reaction at the surface and in the bulk of aerosol particles (KM-SUB). The model is based on the PRA framework of gas–particle interactions (P¨oschl et al., 5 2007), and it includes reversible adsorption, surface reactions and surface-bulk exchange as well as bulk diffusion and reaction. Unlike earlier models, KM-SUB does not require simplifying assumptions about steady-state conditions and radial mixing. The temporal evolution and concentration profiles of volatile and non-volatile species at the gas-particle interface and in the particle bulk can be modeled along with surface 10 concentrations and gas uptake coefficients. In this study we explore and exemplify the effects of bulk diffusion on the rate of reactive gas uptake for a simple reference system, the ozonolysis of oleic acid particles, in comparison to experimental data and earlier model studies. We demonstrate how KM-SUB can be used to interpret and analyze experimental data from laboratory stud15 ies, and how the results can be extrapolated to atmospheric conditions. In particular, we show how interfacial transport and bulk transport, i.e., surface accommodation, bulk accommodation and bulk diffusion, influence the kinetics of the chemical reaction. Sensitivity studies suggest that in fine air particulate matter oleic acid and compounds with similar reactivity against ozone (C=C double bonds) can reach chemical lifetimes of 20 multiple hours only if they are embedded in a (semi-)solid matrix with very low diffusion coefficients (10−10 cm2 s−1). Depending on the complexity of the investigated system, unlimited numbers of volatile and non-volatile species and chemical reactions can be flexibly added and treated with KM-SUB. We propose and intend to pursue the application of KM-SUB 25 as a basis for the development of a detailed master mechanism of aerosol chemistry as well as for the derivation of simplified but realistic parameterizations for large-scale atmospheric and climate models

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

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

Coupling aerosol surface and bulk chemistry with a kinetic double layer model (K2-SUB): oxidation of oleic acid by ozone

We present a kinetic double layer model coupling aerosol surface and bulk chemistry (K2-SUB) based on the PRA framework of gas-particle interactions (Poschl-Rudich-Ammann, 2007). K2-SUB is applied to a popular model system of atmospheric heterogeneous chemistry: the interaction of ozone with oleic acid. We show that our modelling approach allows de-convoluting surface and bulk processes, which has been a controversial topic and remains an important challenge for the understanding and description of atmospheric aerosol transformation. In particular, we demonstrate how a detailed treatment of adsorption and reaction at the surface can be coupled to a description of bulk reaction and transport that is consistent with traditional resistor model formulations. From literature data we have derived a consistent set of kinetic parameters that characterise mass transport and chemical reaction of ozone at the surface and in the bulk of oleic acid droplets. Due to the wide range of rate coefficients reported from different experimental studies, the exact proportions between surface and bulk reaction rates remain uncertain. Nevertheless, the model results suggest an important role of chemical reaction in the bulk and an approximate upper limit of similar to 10(-11) cm(2) s(-1) for the surface reaction rate coefficient. Sensitivity studies show that the surface accommodation coefficient of the gas-phase reactant has a strong non-linear influence on both surface and bulk chemical reactions. We suggest that K2-SUB may be used to design, interpret and analyse future experiments for better discrimination between surface and bulk processes in the oleic acid-ozone system as well as in other heterogeneous reaction systems of atmospheric relevance

Turbulent Diffusion and Turbulent Thermal Diffusion of Aerosols in Stratified Atmospheric Flows

The paper analyzes the phenomenon of turbulent thermal diffusion in the Earth atmosphere, its relation to the turbulent diffusion and its potential impact on aerosol distribution. This phenomenon was predicted theoretically more than 10 years ago and detected recently in the laboratory experiments. This effect causes a non-diffusive flux of aerosols in the direction of the heat flux and results in formation of long-living aerosol layers in the vicinity of temperature inversions. We demonstrated that the theory of turbulent thermal diffusion explains the GOMOS aerosol observations near the tropopause (i.e., the observed shape of aerosol vertical profiles with elevated concentrations located almost symmetrically with respect to temperature profile). In combination with the derived expression for the dependence of the turbulent thermal diffusion ratio on the turbulent diffusion, these measurements yield an independent method for determining the coefficient of turbulent diffusion at the tropopause. We evaluated the impact of turbulent thermal diffusion to the lower-troposphere vertical profiles of aerosol concentration by means of numerical dispersion modelling, and found a regular upward forcing of aerosols with coarse particles affected stronger than fine aerosols.Comment: 19 pages, 10 figure

Dominant aerosol processes during high-pollution episodes over Greater Tokyo

This paper studies two high-pollution episodes over Greater Tokyo: 9 and 10 December 1999, and 31 July and 1 August 2001. Results obtained with the chemistry-transport model (CTM) Polair3D are compared to measurements of inorganic PM2.5. To understand to which extent the aerosol processes modeled in Polair3D impact simulated inorganic PM2.5, Polair3D is run with different options in the aerosol module, e.g. with/without heterogeneous reactions. To quantify the impact of processes outside the aerosol module, simulations are also done with another CTM (CMAQ). In the winter episode, sulfate is mostly impacted by condensation, coagulation, long-range transport, and deposition to a lesser extent. In the summer episode, the effect of long-range transport largely dominates. The impact of condensation/evaporation is dominant for ammonium, nitrate and chloride in both episodes. However, the impact of the thermodynamic equilibrium assumption is limited. The impact of heterogeneous reactions is large for nitrate and ammonium, and taking heterogeneous reactions into account appears to be crucial in predicting the peaks of nitrate and ammonium. The impact of deposition is the same for all inorganic PM2.5. It is small compared to the impact of other processes although it is not negligible. The impact of nucleation is negligible in the summer episode, and small in the winter episode. The impact of coagulation is larger in the winter episode than in the summer episode, because the number of small particles is higher in the winter episode as a consequence of nucleation.Comment: Journal of Geophysical Research D: Atmospheres (15/05/2007) in pres

Influence of organic films on the evaporation and condensation of water in aerosol

Uncertainties in quantifying the kinetics of evaporation and condensation of water from atmospheric aerosol are a significant contributor to the uncertainty in predicting cloud droplet number and the indirect effect of aerosols on climate. The influence of aerosol particle surface composition, particularly the impact of surface active organic films, on the condensation and evaporation coefficients remains ambiguous. Here, we report measurements of the influence of organic films on the evaporation and condensation of water from aerosol particles. Significant reductions in the evaporation coefficient are shown to result when condensed films are formed by monolayers of long-chain alcohols [C(n)H((2n+1))OH], with the value decreasing from 2.4 × 10(−3) to 1.7 × 10(−5) as n increases from 12 to 17. Temperature-dependent measurements confirm that a condensed film of long-range order must be formed to suppress the evaporation coefficient below 0.05. The condensation of water on a droplet coated in a condensed film is shown to be fast, with strong coherence of the long-chain alcohol molecules leading to islanding as the water droplet grows, opening up broad areas of uncoated surface on which water can condense rapidly. We conclude that multicomponent composition of organic films on the surface of atmospheric aerosol particles is likely to preclude the formation of condensed films and that the kinetics of water condensation during the activation of aerosol to form cloud droplets is likely to remain rapid

Transfer Functions and Penetrations of Five Differential Mobility Analyzers for Sub-2 nm Particle Classification

The transfer functions and penetrations of five differential mobility analyzers (DMAs) for sub-2 nm particle classification were evaluated in this study. These DMAs include the TSI nanoDMA, the Caltech radial DMA (RDMA) and nanoRDMA, the Grimm nanoDMA, and the Karlsruhe-Vienna DMA. Measurements were done using tetra-alkyl ammonium ion standards with mobility diameters of 1.16, 1.47, and 1.70 nm. These monomobile ions were generated by electrospray followed by high resolution mobility classification. Measurements were focused at an aerosol-to-sheath flow ratio of 0.1. A data inversion routine was developed to obtain the true transfer function for each test DMA, and these measured transfer functions were compared with theory. DMA penetration efficiencies were also measured. An approximate model for diffusional deposition, based on the modified Gormley and Kennedy equation using an effective length, is given for each test DMA. These results quantitatively characterize the performance of the test DMAs in classifying sub-2 nm particles and can be readily used for DMA data inversion

Transport and diffusion in the embedding map

We study the transport properties of passive inertial particles in a $2-d$ incompressible flows. Here the particle dynamics is represented by the $4-d$ dissipative embedding map of $2-d$ area-preserving standard map which models the incompressible flow. The system is a model for impurity dynamics in a fluid and is characterized by two parameters, the inertia parameter $\alpha$, and the dissipation parameter $\gamma$. We obtain the statistical characterisers of transport for this system in these dynamical regimes. These are, the recurrence time statistics, the diffusion constant, and the distribution of jump lengths. The recurrence time distribution shows a power law tail in the dynamical regimes where there is preferential concentration of particles in sticky regions of the phase space, and an exponential decay in mixing regimes. The diffusion constant shows behaviour of three types - normal, subdiffusive and superdiffusive, depending on the parameter regimes. Phase diagrams of the system are constructed to differentiate different types of diffusion behaviour, as well as the behaviour of the absolute drift. We correlate the dynamical regimes seen for the system at different parameter values with the transport properties observed at these regimes, and in the behaviour of the transients. This system also shows the existence of a crisis and unstable dimension variability at certain parameter values. The signature of the unstable dimension variability is seen in the statistical characterisers of transport. We discuss the implications of our results for realistic systems.Comment: 28 pages, 14 figures, To Appear in Phys. Rev. E; Vol. 79 (2009

Applicability of condensation particle counters to measure atmospheric clusters

This study presents an evaluation of a pulse height condensation particle counter (PH-CPC) and an expansion condensation particle counter (E-CPC) in terms of measuring ambient and laboratory-generated molecular and ion clusters. Ambient molecular cluster concentrations were measured with both instruments as they were deployed in conjunction with an ion spectrometer and other aerosol instruments in Hyytiälä, Finland at the SMEAR II station between 1 March and 30 June 2007. The observed cluster concentrations varied and ranged from some thousands to 100 000 cm -3. Both instruments showed similar (within a factor of ~5) concentrations. An average size of the detected clusters was approximately 1.8 nm. As the atmospheric measurement of sub 2-nm particles and molecular clusters is a challenging task, we conclude that most likely we were unable to detect the smallest clusters. Nevertheless, the reported concentrations are the best estimates to date for minimum cluster concentrations in a boreal forest environment