Process based Modelling of Chemical and Physical Aerosol Properties Relevant for Climate and Health

Atmospheric aerosol particles have substantial influence on climate and air quality. However, the anthropogenic influence on the atmospheric aerosol is still poorly known. This limits the understanding of past and future climate changes. Additionally, both epidemiological and toxicological studies indicate adverse health effects of inhaled aerosol particles. In order to study the effect of atmospheric processes on the particle properties relevant for climate and health, two models were developed and implemented. The first is a 2D-Lagrangian model for Aerosol Dynamics, gas phase CHEMistry and radiative transfer (ADCHEM), which treats the dispersion in the vertical and horizontal direction perpendicular to air mass trajectories. The second model is a kinetic multilayer model for Aerosol Dynamics, gas and particle phase chemistry in laboratory CHAMber environments (ADCHAM). With ADCHAM it is possible to study process based formation and evaporation of secondary organic aerosol particles, and mass transfer limitations and reactions within the particle phase. ADCHEM was used to quantify the anthropogenic influence from the city of Malmö (280 000 inhabitants) in southern Sweden. In Malmö and a few tens of kilometres downwind, the primary particle emissions have a large influence on the particle number concentration. However, more than 2 hours downwind Malmö, the anthropogenic particle mass contribution is dominated by secondary ammonium nitrate. To quantify the direct and indirect climate impact of urban aerosol emissions, the secondary aerosol formation which changes the optical and hygroscopic properties of the primary soot particles, needs to be addressed in future measurements and process modelling. ADCHAM was used to simulate different laboratory chamber experiment, with focus on potential influential but poorly known processes for secondary organic aerosol properties, formation and evaporation rates in the atmosphere (i.e. oligomerization, organic salt formation, salting-out effects, oxidation of organic compounds in the particle phase and mass transfer limitations in the particle phase). The model results reveal that formation of small amounts of low-volatile and long lived oligomers, which accumulate in the particle surface layers, can effectively prevent the evaporation of more volatile compounds. This can significantly prolong the lifetime of SOA in the atmosphere

Experimental investigation into the volatilities of highly oxygenated organic molecules (HOMs)

Secondary organic aerosol (SOA) forms a major part of the tropospheric submicron aerosol. Still, the exact formation mechanisms of SOA have remained elusive. Recently, a newly discovered group of oxidation products of volatile organic compounds (VOCs), highly oxygenated organic molecules (HOMs), have been proposed to be responsible for a large fraction of SOA formation. To assess the potential of HOMs to form SOA and to even take part in new particle formation, knowledge of their exact volatilities is essential. However, due to their exotic, and partially unknown, structures, estimating their volatility is challenging. In this study, we performed a set of continuous flow chamber experiments, supported by box modelling, to study the volatilities of HOMs, along with some less oxygenated compounds, formed in the ozonolysis of α-pinene, an abundant VOC emitted by boreal forests. Along with gaseous precursors, we periodically injected inorganic seed aerosol into the chamber to vary the condensation sink (CS) of low-volatility vapours. We monitored the decrease of oxidation products in the gas phase in response to increasing CS, and were able to relate the responses to the volatilities of the compounds. We found that HOM monomers are mainly of low volatility, with a small fraction being semi-volatile. HOM dimers were all at least low volatility, but probably extremely low volatility; however, our method is not directly able to distinguish between the two. We were able to model the volatility of the oxidation products in terms of their carbon, hydrogen, oxygen and nitrogen numbers. We found that increasing levels of oxygenation correspond to lower volatilities, as expected, but that the decrease is less steep than would be expected based on many existing models for volatility, such as SIMPOL. The hydrogen number of a compound also predicted its volatility, independently of the carbon number, with higher hydrogen numbers corresponding to lower volatilities. This can be explained in terms of the functional groups making up a molecule: high hydrogen numbers are associated with, e.g. hydroxy groups, which lower volatility more than, e.g. carbonyls, which are associated with a lower hydrogen number. The method presented should be applicable to systems other than α-pinene ozonolysis, and with different organic loadings, in order to study different volatility ranges.Secondary organic aerosol (SOA) forms a major part of the tropospheric submicron aerosol. Still, the exact formation mechanisms of SOA have remained elusive. Recently, a newly discovered group of oxidation products of volatile organic compounds (VOCs), highly oxygenated organic molecules (HOMs), have been proposed to be responsible for a large fraction of SOA formation. To assess the potential of HOMs to form SOA and to even take part in new particle formation, knowledge of their exact volatilities is essential. However, due to their exotic, and partially unknown, structures, estimating their volatility is challenging. In this study, we performed a set of continuous flow chamber experiments, supported by box modelling, to study the volatilities of HOMs, along with some less oxygenated compounds, formed in the ozonolysis of alpha-pinene, an abundant VOC emitted by boreal forests. Along with gaseous precursors, we periodically injected inorganic seed aerosol into the chamber to vary the condensation sink (CS) of low-volatility vapours. We monitored the decrease of oxidation products in the gas phase in response to increasing CS, and were able to relate the responses to the volatilities of the compounds. We found that HOM monomers are mainly of low volatility, with a small fraction being semi-volatile. HOM dimers were all at least low volatility, but probably extremely low volatility; however, our method is not directly able to distinguish between the two. We were able to model the volatility of the oxidation products in terms of their carbon, hydrogen, oxygen and nitrogen numbers. We found that increasing levels of oxygenation correspond to lower volatilities, as expected, but that the decrease is less steep than would be expected based on many existing models for volatility, such as SIM-POL. The hydrogen number of a compound also predicted its volatility, independently of the carbon number, with higher hydrogen numbers corresponding to lower volatilities. This can be explained in terms of the functional groups making up a molecule: high hydrogen numbers are associated with, e.g. hydroxy groups, which lower volatility more than, e.g. carbonyls, which are associated with a lower hydrogen number. The method presented should be applicable to systems other than alpha-pinene ozonolysis, and with different organic loadings, in order to study different volatility ranges.Peer reviewe

Aerosol mass yields of selected biogenic volatile organic compounds– a theoretical study with nearly explicit gas-phase chemistry

In this study we modeled secondary organic aerosol (SOA) mass loadings from the oxidation (by O-3, OH and NO3) of five representative biogenic volatile organic compounds (BVOCs): isoprene, endocyclic bond-containing monoterpenes (alpha-pinene and limonene), exocyclic double-bond compound (beta-pinene) and a sesquiterpene (beta-caryophyllene). The simulations were designed to replicate an idealized smog chamber and oxidative flow reactors (OFRs). The Master Chemical Mechanism (MCM) together with the peroxy radical autoxidation mechanism (PRAM) were used to simulate the gas-phase chemistry. The aim of this study was to compare the potency of MCM and MCM + PRAM in predicting SOA formation. SOA yields were in good agreement with experimental values for chamber simulations when MCM + PRAM was applied, while a stand-alone MCM underpredicted the SOA yields. Compared to experimental yields, the OFR simulations using MCM + PRAM yields were in good agreement for BVOCs oxidized by both O-3 and OH. On the other hand, a stand-alone MCM underpredicted the SOA mass yields. SOA yields increased with decreasing temperatures and NO concentrations and vice versa. This highlights the limitations posed when using fixed SOA yields in a majority of global and regional models. Few compounds that play a crucial role (> 95% of mass load) in contributing to SOA mass increase (using MCM + PRAM) are identified. The results further emphasized that incorporating PRAM in conjunction with MCM does improve SOA mass yield estimation.Peer reviewe

Aerosol ageing in an urban plume - implication for climate

The climate effects downwind of an urban area resulting from gaseous and particulate emissions within the city are as yet inadequately quantified. The aim of this work was to estimate these effects for Malmo city in southern Sweden (population 280 000). The chemical and physical particle properties were simulated with a model for Aerosol Dynamics, gas phase CHEMistry and radiative transfer calculations (ADCHEM) following the trajectory movement from upwind of Malmo, through the urban background environment and finally tens and hundreds of kilometers downwind of Malmo. The model results were evaluated using measurements of the particle number size distribution and chemical composition. The total particle number concentration 50 km (similar to 3 h) downwind, in the center of the Malmo plume, is about 3700 cm(-3) of which the Malmo contribution is roughly 30%. Condensation of nitric acid, ammonium and to a smaller extent oxidized organic compounds formed from the emissions in Malmo increases the secondary aerosol formation with a maximum of 0.7-0.8 mu gm(-3) 6 to 18 h downwind of Malmo. The secondary mass contribution dominates over the primary soot contribution from Malmo already 3 to 4 h downwind of the emission sources and contributes to an enhanced total surface direct or indirect aerosol shortwave radiative forcing in the center of the urban plume ranging from -0.3 to -3.3 Wm(-2) depending on the distance from Malmo, and the specific cloud properties

alpha-Pinene Autoxidation Products May Not Have Extremely Low Saturation Vapor Pressures Despite High O:C Ratios

COSMO-RS (conductor-like screening model for real solvents) and three different group-contribution methods were used to compute saturation (subcooled) liquid vapor pressures for 16 possible products of ozone-initiated alpha-pinene autoxidation, with elemental compositions C10H16O4-10 and C20H30O10-12. The saturation vapor pressures predicted by the different methods varied widely. COSMO-RS predicted relatively high saturation vapor pressures values in the range of 10(-6) to 10(-10) bar for the C10H16O4-10 "monomers", and 10(-11) to 10(-16) bar for the C20H30O10-12 "dimers". The group-contribution methods predicted significantly (up to 8 order of magnitude) lower saturation vapor pressures for most of the more highly oxidized monomers. For the differs, the COSMO-RS predictions were within the (wide) range spanned by the three group-contribution methods. The main reason for the discrepancies between the methods is likely that the group-contribution methods do not contain the necessary parameters to accurately treat autoxidation products containing multiple hydroperoxide, peroxy acid or peroxide functional groups, which form intramolecular hydrogen bonds with each other. While the COSMO-RS saturation vapor pressures for these systems may be overestimated, the results strongly indicate that despite their high O:C ratios, the volatilities of the autoxidation products of alpha-pinene (and possibly other atmospherically relevant alkenes) are not necessarily extremely low. In other words, while autoxidation products are able to, adsorb onto aerosol particles, their evaporation back into the gas phase cannot be assumed to be negligible, especially from the smallest nanometer-scale particles. Their observed effective contribution to aerosol particle growth may therefore involve rapid heterogeneous reactions (reactive uptake) rather than effectively irreversible physical absorption.Peer reviewe

Coupling an aerosol box model with one-dimensional flow : a tool for understanding observations of new particle formation events

Field observations of new particle formation and the subsequent particle growth are typically only possible at a fixed measurement location, and hence do not follow the temporal evolution of an air parcel in a Lagrangian sense. Standard analysis for determining formation and growth rates requires that the time-dependent formation rate and growth rate of the particles are spatially invariant; air parcel advection means that the observed temporal evolution of the particle size distribution at a fixed measurement location may not represent the true evolution if there are spatial variations in the formation and growth rates. Here we present a zero-dimensional aerosol box model coupled with one-dimensional atmospheric flow to describe the impact of advection on the evolution of simulated new particle formation events. Wind speed, particle formation rates and growth rates are input parameters that can vary as a function of time and location, using wind speed to connect location to time. The output simulates measurements at a fixed location; formation and growth rates of the particle mode can then be calculated from the simulated observations at a stationary point for different scenarios and be compared with the 'true' input parameters. Hence, we can investigate how spatial variations in the formation and growth rates of new particles would appear in observations of particle number size distributions at a fixed measurement site. We show that the particle size distribution and growth rate at a fixed location is dependent on the formation and growth parameters upwind, even if local conditions do not vary. We also show that different input parameters used may result in very similar simulated measurements. Erroneous interpretation of observations in terms of particle formation and growth rates, and the time span and areal extent of new particle formation, is possible if the spatial effects are not accounted for.Peer reviewe

Implementation of the sectional aerosol module SALSA2.0 into the PALM model system 6.0 : Model development and first evaluation

Urban pedestrian-level air quality is a result of an interplay between turbulent dispersion conditions, background concentrations, and heterogeneous local emissions of air pollutants and their transformation processes. Still, the complexity of these interactions cannot be resolved by the commonly used air quality models. By embedding the sectional aerosol module SALSA2.0 into the large-eddy simulation model PALM, a novel, high-resolution, urban aerosol modelling framework has been developed. The first model evaluation study on the vertical variation of aerosol number concentration and size distribution in a simple street canyon without vegetation in Cambridge, UK, shows good agreement with measurements, with simulated values mainly within a factor of 2 of observations. Dispersion conditions and local emissions govern the pedestrian-level aerosol number concentrations. Out of different aerosol processes, dry deposition is shown to decrease the total number concentration by over 20 %, while condensation and dissolutional increase the total mass by over 10 %. Following the model development, the application of PALM can be extended to local-and neighbourhood-scale air pollution and aerosol studies that require a detailed solution of the ambient flow field.Peer reviewe

Modelling non-equilibrium secondary organic aerosol formation and evaporation with the aerosol dynamics, gas- and particle-phase chemistry kinetic multilayer model ADCHAM

We have developed the novel Aerosol Dynamics, gas- and particle-phase chemistry model for laboratory CHAMber studies (ADCHAM). The model combines the detailed gas-phase Master Chemical Mechanism version 3.2 (MCMv3.2), an aerosol dynamics and particle-phase chemistry module (which considers acid-catalysed oligomerization, heterogeneous oxidation reactions in the particle phase and non-ideal interactions between organic compounds, water and inorganic ions) and a kinetic multilayer module for diffusion-limited transport of compounds between the gas phase, particle surface and particle bulk phase. In this article we describe and use ADCHAM to study (1) the evaporation of liquid dioctyl phthalate (DOP) particles, (2) the slow and almost particle-size-independent evaporation of alpha-pinene ozonolysis secondary organic aerosol (SOA) particles, (3) the mass-transfer-limited uptake of ammonia (NH3) and formation of organic salts between ammonium (NH4+) and carboxylic acids (RCOOH), and (4) the influence of chamber wall effects on the observed SOA formation in smog chambers. ADCHAM is able to capture the observed alpha-pinene SOA mass increase in the presence of NH3(g). Organic salts of ammonium and carboxylic acids predominantly form during the early stage of SOA formation. In the smog chamber experiments, these salts contribute substantially to the initial growth of the homogeneously nucleated particles. The model simulations of evaporating alpha-pinene SOA particles support the recent experimental findings that these particles have a semi-solid tar-like amorphous-phase state. ADCHAM is able to reproduce the main features of the observed slow evaporation rates if the concentration of low-volatility and viscous oligomerized SOA material at the particle surface increases upon evaporation. The evaporation rate is mainly governed by the reversible decomposition of oligomers back to monomers. Finally, we demonstrate that the mass-transfer-limited uptake of condensable organic compounds onto wall-deposited particles or directly onto the Teflon chamber walls of smog chambers can have a profound influence on the observed SOA formation. During the early stage of the SOA formation the wall-deposited particles and walls themselves serve as an SOA sink from the air to the walls. However, at the end of smog chamber experiments the semi-volatile SOA material may start to evaporate from the chamber walls. With these four model applications, we demonstrate that several poorly quantified processes (i.e. mass transport limitations within the particle phase, oligomerization, heterogeneous oxidation, organic salt formation, and chamber wall effects) can have a substantial influence on the SOA formation, lifetime, chemical and physical particle properties, and their evolution. In order to constrain the uncertainties related to these processes, future experiments are needed in which as many of the influential variables as possible are varied. ADCHAM can be a valuable model tool in the design and analysis of such experiments