Search

Polycyclic aromatic hydrocarbons (PAHs) are carcinogenic air pollutants. The dispersion of PAHs in the atmosphere is influenced by gas–particle partitioning and chemical loss. These processes are closely interlinked and may occur at vastly differing timescales, which complicates their mathematical description in chemical transport models. Here, we use a kinetic model that explicitly resolves mass transport and chemical reactions in the gas and particle phases to describe and explore the dynamic and non-equilibrium interplay of gas–particle partitioning and chemical losses of PAHs on soot particles. We define the equilibration timescale τeq of gas–particle partitioning as the e-folding time for relaxation of the system to the partitioning equilibrium. We find this metric to span from seconds to hours depending on temperature, particle surface area, and the type of PAH. The equilibration time can be approximated using a time-independent equation, τeq≈1kdes+kads, which depends on the desorption rate coefficient kdes and adsorption rate coefficient kads, both of which can be calculated from experimentally accessible parameters. The model reveals two regimes in which different physical processes control the equilibration timescale: a desorption-controlled and an adsorption-controlled regime. In a case study with the PAH pyrene, we illustrate how chemical loss can perturb the equilibrium particulate fraction at typical atmospheric concentrations of O3 and OH. For the surface reaction with O3, the perturbation is significant and increases with the gas-phase concentration of O3. Conversely, perturbations are smaller for reaction with the OH radical, which reacts with pyrene on both the surface of particles and in the gas phase. Global and regional chemical transport models typically approximate gas–particle partitioning with instantaneous-equilibration approaches. We highlight scenarios in which these approximations deviate from the explicitly coupled treatment of gas–particle partitioning and chemistry presented in this study. We find that the discrepancy between solutions depends on the operator-splitting time step and the choice of time step can help to minimize the discrepancy. The findings and techniques presented in this work not only are relevant for PAHs but can also be applied to other semi-volatile substances that undergo chemical reactions and mass transport between the gas and particle phase.

Atmospheric aerosols can exist in amorphous semi-solid or glassy phase states whose viscosity varies with atmospheric temperature and relative humidity. The temperature and humidity dependence of viscosity has been hypothesized to be predictable from the combination of a water–organic binary mixing rule of the glass transition temperature, a glass-transition-temperature-scaled viscosity fragility parameterization, and a water uptake parameterization. This work presents a closure study between predicted and observed viscosity for sucrose and citric acid. Viscosity and glass transition temperature as a function of water content are compiled from literature data and used to constrain the fragility parameterization. New measurements characterizing viscosity of sub-100 nm particles using the dimer relaxation method are presented. These measurements extend the available data of temperature- and humidity-dependent viscosity to −28 ∘C. Predicted relationships agree well with observations at room temperature and with measured isopleths of constant viscosity at ∼107 Pa s at temperatures warmer than −28 ∘C. Discrepancies at colder temperatures are observed for sucrose particles. Simulations with the kinetic multi-layer model of gas–particle interactions suggest that the observed deviations at colder temperature for sucrose can be attributed to kinetic limitations associated with water uptake at the timescales of the dimer relaxation experiments. Using the available information, updated equilibrium phase-state diagrams (-80∘C

Secondary organic aerosols (SOA) account for a substantial fraction of air particulate matter, and SOA formation is often modeled assuming rapid establishment of gas–particle equilibrium. Here, we estimate the characteristic timescale for SOA to achieve gas–particle equilibrium under a wide range of temperatures and relative humidities using a state-of-the-art kinetic flux model. Equilibration timescales were calculated by varying particle phase state, size, mass loadings, and volatility of organic compounds in open and closed systems. Model simulations suggest that the equilibration timescale for semi-volatile compounds is on the order of seconds or minutes for most conditions in the planetary boundary layer, but it can be longer than 1 h if particles adopt glassy or amorphous solid states with high glass transition temperatures at low relative humidity. In the free troposphere with lower temperatures, it can be longer than hours or days, even at moderate or relatively high relative humidities due to kinetic limitations of bulk diffusion in highly viscous particles. The timescale of partitioning of low-volatile compounds into highly viscous particles is shorter compared to semi-volatile compounds in the closed system, as it is largely determined by condensation sink due to very slow re-evaporation with relatively quick establishment of local equilibrium between the gas phase and the near-surface bulk. The dependence of equilibration timescales on both volatility and bulk diffusivity provides critical insights into thermodynamic or kinetic treatments of SOA partitioning for accurate predictions of gas- and particle-phase concentrations of semi-volatile compounds in regional and global chemical transport models.

The effect of relative humidity (RH) on the chemical composition of secondary organic aerosol (SOA) formed from low-NOx toluene oxidation in the absence of seed particles was investigated. SOA samples were prepared in an aerosol smog chamber at x toluene SOA. The volatility distributions of the SOA compounds, estimated from the distribution of molecular formulas using the molecular corridor approach, confirmed that low-NOx toluene SOA became more volatile on average under high-RH conditions. In contrast, the effect of RH on SOA mass loading was found to be much smaller for high-NOx toluene SOA. The observed increase in the oligomer fraction and particle mass loading under dry conditions were attributed to the enhancement of condensation reactions, which produce water and oligomers from smaller compounds in low-NOx toluene SOA. The reduction in the fraction of oligomeric compounds under humid conditions is predicted to partly counteract the previously observed enhancement in the toluene SOA yield driven by the aerosol liquid water chemistry in deliquesced inorganic seed particles.

We report the first measurements of HO2 uptake coefficients, γ, for secondary organic aerosol (SOA) particles and for the well-studied model compound sucrose which we doped with copper(II). Above 65 % relative humidity (RH), γ for copper(II)-doped sucrose aerosol particles equalled the surface mass accommodation coefficient α = 0.22 ± 0.06, but it decreased to γ = 0.012 ± 0.007 upon decreasing the RH to 17 %. The trend of γ with RH can be explained by an increase in aerosol viscosity and the contribution of a surface reaction, as demonstrated using the kinetic multilayer model of aerosol surface and bulk chemistry (KM-SUB). At high RH the total uptake was driven by reaction in the near-surface bulk and limited by mass accommodation, whilst at low RH it was limited by surface reaction. SOA from two different precursors, α-pinene and 1,3,5-trimethylbenzene (TMB), was investigated, yielding low uptake coefficients of γ

Diffusion of organic vapours to the surface of aerosol or cloud particles is an important step for the formation and transformation of atmospheric particles. So far, however, a database of gas phase diffusion coefficients for organic compounds of atmospheric interest has not been available. In this work we have compiled and evaluated gas phase diffusivities (pressure-independent diffusion coefficients) of organic compounds reported by previous experimental studies, and we compare the measurement data to estimates obtained with Fuller's semi-empirical method. The difference between measured and estimated diffusivities are mostly < 10%. With regard to gas-particle interactions, different gas molecules, including both organic and inorganic compounds, exhibit similar Knudsen numbers (Kn) although their gas phase diffusivities may vary over a wide range. This is because different trace gas molecules have similar mean free paths in air at a given pressure. Thus, we introduce the pressure-normalised mean free path, λP ≈ 100 nm atm, as a near-constant generic parameter that can be used for approximate calculation of Knudsen numbers as a simple function of gas pressure and particle diameter to characterise the influence of gas phase diffusion on the uptake of gases by aerosol or cloud particles. We use a kinetic multilayer model of gas-particle interaction to illustrate the effects of gas phase diffusion on the condensation of organic compounds with different volatilities. The results show that gas phase diffusion can play a major role in determining the growth of secondary organic aerosol particles by condensation of low-volatility organic vapours.

Soot particles are the most efficient light absorbing aerosol species in the atmosphere, playing an important role as a driver of global warming. Their climate effects strongly depend on their mixing state, which significantly changes their light absorbing capability and cloud condensation nuclei (CCN) activity. Therefore, knowledge about the mixing state of soot and its aging mechanism becomes an important topic in the atmospheric sciences. The size-resolved (30–320 nm diameter) mixing state of soot particles in polluted megacity air was measured at a suburban site (Yufa) during the CAREBeijing 2006 campaign in Beijing, using a volatility tandem differential mobility analyzer (VTDMA). Particles in this size range with non-volatile residuals at 300 °C were considered to be soot particles. On average, the number fraction of internally mixed soot in total soot particles (Fin), decreased from 0.80 to 0.57 when initial Dp increased from 30 to 320 nm. Further analysis reveals that: (1) Fin was well correlated with the aerosol hygroscopic mixing state measured by a CCN counter. More externally mixed soot particles were observed when particles showed more heterogeneous features with regard to hygroscopicity. (2) Fin had pronounced diurnal cycles. For particles in the accumulation mode (Dp at 100–320 nm), largest Fin were observed at noon time, with "apparent" turnover rates (kex → in) up to 7.8% h−1. (3) Fin was subject to competing effects of both aging and emissions. While aging increases Fin by converting externally mixed soot particles into internally mixed ones, emissions tend to reduce Fin by emitting more fresh and externally mixed soot particles. Similar competing effects were also found with air mass age indicators. (4) Under the estimated emission intensities, actual turnover rates of soot (kex → in) up to 20% h−1 were derived, which showed a pronounced diurnal cycle peaking around noon time. This result confirms that (soot) particles are undergoing fast aging/coating with the existing high levels of condensable vapors in the megacity Beijing. (5) Diurnal cycles of Fin were different between Aitken and accumulation mode particles, which could be explained by the faster growth of smaller Aitken mode particles into larger size bins. To improve the Fin prediction in regional/global models, we suggest parameterizing Fin by an air mass aging indicator, i.e., Fin = a + bx, where a and b are empirical coefficients determined from observations, and x is the value of an air mass age indicator. At the Yufa site in the North China Plain, fitted coefficients (a, b) were determined as (0.57, 0.21), (0.47, 0.21), and (0.52, 0.0088) for x (indicators) as [NOz]/[NOy], [E]/[X] ([ethylbenzene]/[m,p-xylene]) and ([IM] + [OM])/[EC] ([inorganic + organic matter]/[elemental carbon]), respectively. Such a parameterization consumes little additional computing time, but yields a more realistic description of Fin compared with the simple treatment of soot mixing state in regional/global models.