Sectional Aerosol Research Articles

MATCH-SALSA – Multi-scale Atmospheric Transport and CHemistry model coupled to the SALSA aerosol microphysics model – Part 1: Model description and evaluation

Abstract. We have implemented the sectional aerosol dynamics model SALSA (Sectional Aerosol module for Large Scale Applications) in the European-scale chemistry-transport model MATCH (Multi-scale Atmospheric Transport and Chemistry). The new model is called MATCH-SALSA. It includes aerosol microphysics, with several formulations for nucleation, wet scavenging and condensation. The model reproduces observed higher particle number concentration (PNC) in central Europe and lower concentrations in remote regions. The modeled PNC size distribution peak occurs at the same or smaller particle size as the observed peak at four measurement sites spread across Europe. Total PNC is underestimated at northern and central European sites and accumulation-mode PNC is underestimated at all investigated sites. The low nucleation rate coefficient used in this study is an important reason for the underestimation. On the other hand, the model performs well for particle mass (including secondary inorganic aerosol components), while elemental and organic carbon concentrations are underestimated at many of the sites. Further development is needed, primarily for treatment of secondary organic aerosol, in terms of biogenic emissions and chemical transformation. Updating the biogenic secondary organic aerosol (SOA) scheme will likely have a large impact on modeled PM2.5 and also affect the model performance for PNC through impacts on nucleation and condensation.

Near-road modeling and measurement of cerium-containing particles generated by nanoparticle diesel fuel additive use.

Cerium oxide nanoparticles (nCe) are used as a fuel-borne catalyst in diesel engines to reduce particulate emissions, yet the environmental and human health impacts of the exhaust particles are not well understood. To bridge the gap between emission measurements and ambient impacts, size-resolved measurements of particle composition and mass concentration have been performed in Newcastle-upon-Tyne, United Kingdom, where buses have used an nCe additive since 2005. These observations show that the noncrustal cerium fraction thought to be associated with the use of nCe has a mass concentration ∼ 0.3 ng m(-3) with a size distribution peaking at 100-320 nm in aerodynamic diameter. Simulations with a near-roadway multicomponent sectional aerosol dynamic model predict that the use of nCe additives increases the number concentration of nuclei mode particles (<50 nm in diameter) while decreasing the total mass concentration. The near-road model predicts a downwind mass size distribution of cerium-containing particles peaking at 150 nm in aerodynamic diameter, a value similar to that measured for noncrustal cerium in Newcastle. This work shows that both the emission and atmospheric transformation of cerium-containing particles needs to be taken into account by regional modelers, exposure scientists, and policymakers when determining potential environmental and human health impacts.

A computational and experimental study of coflow laminar methane/air diffusion flames: Effects of fuel dilution, inlet velocity, and gravity

The influences of fuel dilution, inlet velocity, and gravity on the shape and structure of laminar coflow CH4–air diffusion flames were investigated computationally and experimentally. A series of nitrogen-diluted flames measured in the Structure and Liftoff in Combustion Experiment (SLICE) on board the International Space Station was assessed numerically under microgravity (μg) and normal gravity (1g) conditions with CH4 mole fraction ranging from 0.4 to 1.0 and average inlet velocity ranging from 23 to 90cm/s. Computationally, the MC-Smooth vorticity–velocity formulation was employed to describe the reactive gaseous mixture, and soot evolution was modeled by sectional aerosol equations. The governing equations and boundary conditions were discretized on a two-dimensional computational domain by finite differences, and the resulting set of fully coupled, strongly nonlinear equations was solved simultaneously at all points using a damped, modified Newton’s method. Experimentally, flame shape and soot temperature were determined by flame emission images recorded by a digital color camera. Very good agreement between computation and measurement was obtained, and the conclusions were as follows. (1) Buoyant and nonbuoyant luminous flame lengths are proportional to the mass flow rate of the fuel mixture; computed and measured nonbuoyant flames are noticeably longer than their 1g counterparts; the effect of fuel dilution on flame shape (i.e., flame length and flame radius) is negligible when the flame shape is normalized by the methane flow rate. (2) Buoyancy-induced reduction of the flame radius through radially inward convection near the flame front is demonstrated. (3) Buoyant and nonbuoyant flame structure is mainly controlled by the fuel mass flow rate, and the effects from fuel dilution and inlet velocity are secondary.

An experimental and computational study of soot formation in a coflow jet flame under microgravity and normal gravity

Upon the completion of the Structure and Liftoff in Combustion Experiment (SLICE) in March 2012, a comprehensive and unique set of microgravity coflow diffusion flame data was obtained. This data covers a range of conditions from weak flames near extinction to strong, highly sooting flames, and enabled the study of gravitational effects on phenomena such as liftoff, blowout and soot formation. The microgravity experiment was carried out in the Microgravity Science Glovebox (MSG) on board the International Space Station (ISS), while the normal gravity experiment was performed at Yale utilizing a copy of the flight hardware. Computational simulations of microgravity and normal gravity flames were also carried out to facilitate understanding of the experimental observations. This paper focuses on the different sooting behaviors of CH4 coflow jet flames in microgravity and normal gravity. The unique set of data serves as an excellent test case for developing more accurate computational models. Experimentally, the flame shape and size, lift-off height, and soot temperature were determined from line-of-sight flame emission images taken with a color digital camera. Soot volume fraction was determined by performing an absolute light calibration using the incandescence from a flame-heated thermocouple. Computationally, the MC-Smooth vorticity–velocity formulation was employed to describe the chemically reacting flow, and the soot evolution was modeled by the sectional aerosol equations. The governing equations and boundary conditions were discretized on an axisymmetric computational domain by finite differences, and the resulting system of fully coupled, highly nonlinear equations was solved by a damped, modified Newton’s method. The microgravity sooting flames were found to have lower soot temperatures and higher volume fraction than their normal gravity counterparts. The soot distribution tends to shift from the centerline of the flame to the wings from normal gravity to microgravity.

Impact of radiatively interactive dust aerosols in the NASA GEOS‐5 climate model: Sensitivity to dust particle shape and refractive index

AbstractThe radiative effects of Saharan dust aerosols are investigated in the NASA GEOS‐5 atmospheric general circulation model. A sectional aerosol microphysics model (CARMA) is run online in GEOS‐5. CARMA treats the dust aerosol lifecycle, and its tracers are radiatively coupled to GEOS‐5. A series of AMIP‐style simulations are performed, in which input dust optical properties (particle shape and refractive index) are varied. Simulated dust distributions for summertime Saharan dust compare well to observations, with best results found when the most absorbing dust optical properties are assumed. Dust absorption leads to a strengthening of the summertime Hadley cell circulation, increased dust lofting to higher altitudes, and a strengthening of the African easterly jet, resulting in increased dust atmospheric lifetime and farther northward and westward transport. We find a positive feedback of dust radiative forcing on emissions, in contrast with previous studies, which we attribute to our having a relatively strong longwave forcing caused by our simulating larger effective particle sizes. This longwave forcing reduces the magnitude of midday net surface cooling relative to other studies, and leads to a nighttime warming that results in higher nighttime wind speeds and dust emissions. The radiative effects of dust particle shape have only minor impact on transport and emissions, with small (~5%) impact on top of atmosphere shortwave forcing, in line with previous studies, but relatively more pronounced effects on shortwave atmospheric heating and surface forcing (~20% increase in atmospheric forcing for spheroids). Shape effects on longwave heating terms are of order ~10%.

Numerical Modeling of an RF Argon–Silane Plasma with Dust Particle Nucleation and Growth

A one-dimensional numerical model and simulation results are presented for a capacitively-coupled radio frequency parallel-plate argon-silane dusty plasma. The model includes self-consistently coupled numerical modules, including a plasma fluid model, a sectional aerosol model, and a simple chemistry model to predict rates of particle nucle- ation and surface growth. Operating conditions considered include 13.56 MHz frequency, 100 mTorr pressure, a 4-cm electrode gap, gas flow through the top electrode with a 30:1 ratio of argon to silane, and applied radio frequency voltage amplitude of either 100 or 250 V. In the higher voltage case two lobes of relatively large particles are formed by ion drag, while fresh nucleation occurs in the void between these lobes. It is shown that the reason that fresh nucleation occurs in the void involves an interplay among several coupled phenomena, including nanoparticle transport, the plasma potential profile, and trapping of silicon hydride anions that drive nucleation in this system.

Modal Bin Hybrid Model: A surface area consistent, triple‐moment sectional method for use in process‐oriented modeling of atmospheric aerosols

AbstractA triple‐moment sectional (TMS) aerosol dynamics model, Modal Bin Hybrid Model (MBHM), has been developed. In addition to number and mass (volume), surface area is predicted (and preserved), which is important for aerosol processes and properties such as gas‐to‐particle mass transfer, heterogeneous reaction, and light extinction cross section. The performance of MBHM was evaluated against double‐moment sectional (DMS) models with coarse (BIN4) to very fine (BIN256) size resolutions for simulating evolution of particles under simultaneously occurring nucleation, condensation, and coagulation processes (BINx resolution uses x sections to cover the 1 nm to 1 µm size range). Because MBHM gives a physically consistent form of the intrasectional distributions, errors and biases of MBHM at BIN4‐8 resolution were almost equivalent to those of DMS at BIN16–32 resolution for various important variables such as the moments Mk (k: 0, 2, 3), dMk/dt, and the number and volume of particles larger than a certain diameter. Another important feature of MBHM is that only a single bin is adequate to simulate full aerosol dynamics for particles whose size distribution can be approximated by a single lognormal mode. This flexibility is useful for process‐oriented (multicategory and/or mixing state) modeling: Primary aerosols whose size parameters would not differ substantially in time and space can be expressed by a single or a small number of modes, whereas secondary aerosols whose size changes drastically from 1 to several hundred nanometers can be expressed by a number of modes. Added dimensions can be applied to MBHM to represent mixing state or photochemical age for aerosol mixing state studies.

Representation of nucleation mode microphysics in a global aerosol model with sectional microphysics

Abstract. In models, nucleation mode (1 nm &lt; Dp &lt; 10 nm) particle microphysics can be represented explicitly with aerosol microphysical processes or can be parameterized to obtain the growth and survival of nuclei to the model's lower size boundary. This study investigates how the representation of nucleation mode microphysics impacts aerosol number predictions in the TwO-Moment Aerosol Sectional (TOMAS) aerosol microphysics model running with the GISS GCM II-prime by varying its lowest diameter boundary: 1 nm, 3 nm, and 10 nm. The model with the 1 nm boundary simulates the nucleation mode particles with fully resolved microphysical processes, while the model with the 10 nm and 3 nm boundaries uses a nucleation mode dynamics parameterization to account for the growth of nucleated particles to 10 nm and 3 nm, respectively. We also investigate the impact of the time step for aerosol microphysical processes (a 10 min versus a 1 h time step) to aerosol number predictions in the TOMAS models with explicit dynamics for the nucleation mode particles (i.e., 3 nm and 1 nm boundary). The model with the explicit microphysics (i.e., 1 nm boundary) with the 10 min time step is used as a numerical benchmark simulation to estimate biases caused by varying the lower size cutoff and the time step. Different representations of the nucleation mode have a significant effect on the formation rate of particles larger than 10 nm from nucleated particles (J10) and the burdens and lifetimes of ultrafine-mode (10 nm &amp;amp;leq; Dp &amp;amp;leq; 70 nm) particles but have less impact on the burdens and lifetimes of CCN-sized particles. The models using parameterized microphysics (i.e., 10 nm and 3 nm boundaries) result in higher J10 and shorter coagulation lifetimes of ultrafine-mode particles than the model with explicit dynamics (i.e., 1 nm boundary). The spatial distributions of CN10 (Dp &amp;amp;geq; 10 nm) and CCN(0.2%) (i.e., CCN concentrations at 0.2% supersaturation) are moderately affected, especially CN10 predictions above ~ 700 hPa where nucleation contributes most strongly to CN10 concentrations. The lowermost-layer CN10 is substantially improved with the 3 nm boundary (compared to 10 nm) in most areas. The overprediction in CN10 with the 3 nm and 10 nm boundaries can be explained by the overprediction of J10 or J3 with the parameterized microphysics, possibly due to the instantaneous growth rate assumption in the survival and growth parameterization. The errors in CN10 predictions are sensitive to the choice of the lower size boundary but not to the choice of the time step applied to the microphysical processes. The spatial distribution of CCN(0.2%) with the 3 nm boundary is almost identical to that with the 1 nm boundary, but that with the 10 nm boundary can differ more than 10–40% in some areas. We found that the deviation in the 10 nm simulations is partly due to the longer time step (i.e., 1 h time step used in the 10 nm simulations compared to 10 min time step used in the benchmark simulations), but, even with the same time step, the 10 nm cutoff showed noticeably higher errors than the 3 nm cutoff. In conclusion, we generally recommend using a lower diameter boundary of 3 nm for studies focused on aerosol indirect effects but down to 1 nm boundary for studies focused on CN10 predictions or nucleation.

Anthropogenic contribution to cloud condensation nuclei and the first aerosol indirect climate effect

Atmospheric particles influence the climate indirectly by acting as cloud condensation nuclei (CCN). The first aerosol indirect radiative forcing (FAIRF) constitutes the largest uncertainty among the radiative forcings quantified by the latest IPCC report (IPCC2007) and is a major source of uncertainty in predicting climate change. Here, we investigate the anthropogenic contribution to CCN and associated FAIRF using a state-of-the-art global chemical transport and aerosol model (GEOS-Chem/APM) that contains a number of advanced features (including sectional particle microphysics, online comprehensive chemistry, consideration of all major aerosol species, online aerosol–cloud–radiation calculation, and usage of more accurate assimilated meteorology). The model captures the absolute values and spatial distributions of CCN concentrations measured in situ around the globe. We show that anthropogenic emissions increase the global mean CCN in the lower troposphere by ∼60–80% and cloud droplet number concentration by ∼40%. The global mean FAIRF based on GEOS-Chem/APM is −0.75 W m−2, close to the median values of both IPCC2007 and post-IPCC2007 studies. To the best of our knowledge, this is the first time that a global sectional aerosol model with full online chemistry and considering all major aerosol species (including nitrate, ammonium, and second organic aerosols) has been used used to calculate FAIRF.

Microphysical simulations of large volcanic eruptions: Pinatubo and Toba

AbstractSimulations of stratospheric clouds from eruptions ranging in size from the 1991 eruption of Mount Pinatubo to that of Toba 74,000 years ago have been completed using a 3D microphysical sectional aerosol model advectively coupled to a general circulation model with prognostic chemistry (Whole Atmosphere Community Climate Model/Community Aerosol and Radiation Model for Atmospheres). For Pinatubo, properties of the aerosol cloud peak within the ranges derived from observations in the Northern Hemisphere, but reduce faster than observed, and a general low bias is found in the Southern Hemisphere. These biases could be reduced by adding aerosol radiative coupling, a quasi‐biennial oscillation, and the Cerro Hudson eruption to the model. Simulations of eruptions 10 times and 100 times larger than Pinatubo suggest burdens and Aerosol Optical Depth increase less than linearly (a 100‐fold injection increase produces a 20‐fold AOD increase) due to particle growth and sedimentation, consistent with previous work that also found the radiative forcings from large eruptions to be self‐limiting. Global‐averaged AOD remains elevated for 1, 2, and 4 years, respectively, for the three simulated eruptions. The inclusion of van der Waals forces in our coagulation scheme increases peak effective radius and reduces peak AOD by about 10–20%, with bigger effects for larger eruptions. Our simulations find peak mode size to vary by up to an order of magnitude and mode width to vary by up to 50%, suggesting that two‐moment modal models may not accurately capture the evolving size distribution. These simulations suggest the value of including van der Waals forces in the coagulation scheme and sectional size distributions in climate models.

GEM-AQ/EC, an on-line global multi-scale chemical weather modelling system: model development and evaluation of global aerosol climatology

Abstract. A global air quality modeling system GEM-AQ/EC was developed by implementing tropospheric chemistry and aerosol processes on-line into the Global Environmental Multiscale weather prediction model – GEM. Due to the multi-scale features of the GEM, the integrated model, GEM-AQ/EC, is able to investigate chemical weather at scales from global to urban domains. The current chemical mechanism is comprised of 50 gas-phase species, 116 chemical and 19 photolysis reactions, and is complemented by a sectional aerosol module CAM (The Canadian Aerosol Module) with 5 aerosols types: sulphate, black carbon, organic carbon, sea-salt and soil dust. Monthly emission inventories of black carbon and organic carbon from boreal and temperate vegetation fires were assembled using the most reliable areas burned datasets by countries, from statistical databases and derived from remote sensing products of 1995–2004. The model was run for ten years from from 1995–2004 with re-analyzed meteorology on a global uniform 1° × 1° horizontal resolution domain and 28 hybrid levels extending up to 10 hPa. The simulating results were compared with various observations including surface network around the globe and satellite data. Regional features of global aerosols are reasonably captured including emission, surface concentrations and aerosol optical depth. For various types of aerosols, satisfactory correlations were achieved between modeled and observed with some degree of systematic bias possibly due to large uncertainties in the emissions used in this study. A global distribution of natural aerosol contributions to the total aerosols is obtained and compared with observations.

Evaluation of the sectional aerosol microphysics module SALSA implementation in ECHAM5-HAM aerosol-climate model

Abstract. We present the implementation and evaluation of a sectional aerosol microphysics module SALSA within the aerosol-climate model ECHAM5-HAM. This aerosol microphysics module has been designed to be flexible and computationally efficient so that it can be implemented in regional or global scale models. The computational efficiency has been achieved by minimising the number of variables needed to describe the size and composition distribution. The aerosol size distribution is described using 10 size classes with parallel sections which can have different chemical compositions. Thus in total, the module tracks 20 size sections which cover diameters ranging from 3 nm to 10 μm and are divided into three subranges, each with an optimised selection of processes and compounds. The implementation of SALSA into ECHAM5-HAM includes the main aerosol processes in the atmosphere: emissions, removal, radiative effects, liquid and gas phase sulphate chemistry, and the aerosol microphysics. The aerosol compounds treated in the module are sulphate, organic carbon, sea salt, black carbon, and mineral dust. In its default configuration, ECHAM5-HAM treats aerosol size distribution using the modal method. In this implementation, the aerosol processes were converted to be used in a sectional model framework. The ability of the module to describe the global aerosol properties was evaluated by comparing against (1) measured continental and marine size distributions, (2) observed variability of continental number concentrations, (3) measured sulphate, organic carbon, black carbon and sea-salt mass concentrations, (4) observations of aerosol optical depth (AOD) and other aerosol optical properties from satellites and AERONET network, (5) global aerosol budgets and concentrations from previous model studies, and (6) model results using M7, which is the default aerosol microphysics module in ECHAM5-HAM. The evaluation shows that the global aerosol properties can be reproduced reasonably well using a coarse resolution of 10 sections in size space. The simulated global aerosol budgets are within the range of previous studies. Surface concentrations of sulphate and carbonaceous species have an annual mean within a factor of two of the observations. The simulated sea-salt concentrations reproduce the observations within a factor of two, apart from the Southern Ocean over which the concentrations are within a factor of five. Regionally, AOD is in a relatively good agreement with the observations (within a factor of two). At mid-latitudes the observed AOD is captured well, while at high-latitudes as well as in some polluted and dust regions the modelled AOD is significantly lower than observed. Regarding most of the investigated aerosol properties, the SALSA and the modal aerosol module M7 perform comparably well against observations. However, SALSA reproduces the observed number concentrations and the size distribution of CCN sized particles much more accurately than M7, and is therefore a good choice for aerosol-cloud interaction studies in global models. Our study also shows that when activation type nucleation in the boundary layer is included, the observed concentration of particles under 50 nm in diameter are reproduced much better compared to when only binary nucleation in the free troposphere is assumed.

Microphysical simulations of sulfur burdens from stratospheric sulfur geoengineering

Abstract. Recent microphysical studies suggest that geoengineering by continuous stratospheric injection of SO2 gas may be limited by the growth of the aerosols. We study the efficacy of SO2, H2SO4 and aerosol injections on aerosol mass and optical depth using a three-dimensional general circulation model with sulfur chemistry and sectional aerosol microphysics (WACCM/CARMA). We find increasing injection rates of SO2 in a narrow band around the equator to have limited efficacy while broadening the injecting zone as well as injecting particles instead of SO2 gas increases the sulfate burden for a given injection rate, in agreement with previous work. We find that injecting H2SO4 gas instead of SO2 does not discernibly alter sulfate size or mass, in contrast with a previous study using a plume model with a microphysical model. However, the physics and chemistry in aircraft plumes, which are smaller than climate model grid cells, need to be more carefully considered. We also find significant perturbations to tropospheric aerosol for all injections studied, particularly in the upper troposphere and near the poles, where sulfate burden increases by up to 100 times. This enhanced burden could have implications for tropospheric radiative forcing and chemistry. These results highlight the need to mitigate greenhouse gas emissions rather than attempt to cool the planet through geoengineering, and to further study geoengineering before it can be seriously considered as a climate intervention option.

Intercomparison of modal and sectional aerosol microphysics representations within the same 3-D global chemical transport model

Abstract. In the most advanced aerosol-climate models it is common to represent the aerosol particle size distribution in terms of several log-normal modes. This approach, motivated by computational efficiency, makes assumptions about the shape of the particle distribution that may not always capture the properties of global aerosol. Here, a global modal aerosol microphysics module (GLOMAP-mode) is evaluated and improved by comparing against a sectional version (GLOMAP-bin) and observations in the same 3-D global offline chemistry transport model. With both schemes, the model captures the main features of the global particle size distribution, with sub-micron aerosol approximately unimodal in continental regions and bi-modal in marine regions. Initial bin-mode comparisons showed that the current values for two size distribution parameter settings in the modal scheme (mode widths and inter-modal separation sizes) resulted in clear biases compared to the sectional scheme. By adjusting these parameters in the modal scheme, much better agreement is achieved against the bin scheme and observations. Annual mean surface-level mass of sulphate, sea-salt, black carbon (BC) and organic carbon (OC) are within 25% in the two schemes in nearly all regions. Surface level concentrations of condensation nuclei (CN), cloud condensation nuclei (CCN), surface area density and condensation sink also compare within 25% in most regions. However, marine CCN concentrations between 30° N and 30° S are systematically 25–60% higher in the modal model, which we attribute to differences in size-resolved particle growth or cloud-processing. Larger differences also exist in regions or seasons dominated by biomass burning and in free-troposphere and high-latitude regions. Indeed, in the free-troposphere, GLOMAP-mode BC is a factor 2–4 higher than GLOMAP-bin, likely due to differences in size-resolved scavenging. Nevertheless, in most parts of the atmosphere, we conclude that bin-mode differences are much less than model-observation differences, although some processes are missing in these runs which may pose a bigger challenge to modal schemes (e.g., boundary layer nucleation and ultra-fine sea-spray). The findings here underline the need for a spectrum of complexity in global models, with size-resolved aerosol properties predicted by modal schemes needing to be continually benchmarked and improved against freely evolving sectional schemes and observations.

A study of new particle formation in the marine boundary layer over the central Arctic Ocean using a flexible multicomponent aerosol dynamic model

ABSTRACTEnhancement of number concentrations of particles with sizes less than 25 nm diameter has been frequently observed in the boundary layer over the central Arctic Ocean during summer. The sectional aerosol dynamics model for Marine Aerosol Formation (MAFOR) was applied to evaluate the capability of different nucleation mechanisms to reproduce nucleation events observed during three expeditions (1996, 2001 and 2008) onboard the Swedish icebreaker Oden. Model calculations suggest that a source rate of a condensable organic vapour (OV) of about 2×105 cm−3 s−1 is required to reproduce observed growth of nucleation mode particles. Nucleation rates predicted by the newly proposed combined sulphuric acid nucleation mechanism, which best described new particle formation in the Arctic, ranged from 0.04 to 0.1 cm−3 s−1. This mechanism additively combines ion-mediated nucleation and cluster activation, and treats condensation of OV without correction of the Kelvin effect. In several events, the simultaneous number enhancement of particles in the 20–50 nm size range remained unexplained by the nucleation mechanisms. This lends support to alternative theories such as the fragmentation of marine gels (≈200–500 nm diameter in size) by physical or chemical processes.

Microphysical simulations of new particle formation in the upper troposphere and lower stratosphere

Abstract. Using a three-dimensional general circulation model with sulfur chemistry and sectional aerosol microphysics (WACCM/CARMA), we studied aerosol formation and microphysics in the upper troposphere and lower stratosphere (UTLS) as well as the middle and upper stratosphere based on three nucleation schemes (two binary homogeneous schemes and an ion-mediated scheme related to one of the binary schemes). Simulations suggest that ion-mediated nucleation rates in the UTLS are 25 % higher than its related binary scheme, but that the rates predicted by the two binary schemes vary by two orders of magnitude. None of the nucleation schemes is superior at matching the limited observations available at the smallest sizes. However, it is found that coagulation, not nucleation, controls number concentration at sizes greater than approximately 10 nm. Therefore, based on this study, processes relevant to atmospheric chemistry and radiative forcing in the UTLS are not sensitive to the choice of nucleation schemes. The dominance of coagulation over other microphysical processes in the UTLS is consistent with other recent work using microphysical models. Simulations using all three nucleation schemes compare reasonably well to observations of size distributions, number concentration across latitude, and vertical profiles of particle mixing ratio in the UTLS. Interestingly, we find that we need to include Van der Waals forces in our coagulation scheme to match the UTLS aerosol concentrations. We conclude that this model can reasonably represent sulfate microphysical processes in the UTLS, and that the properties of particles at atmospherically relevant sizes appear to be insensitive to the details of the nucleation scheme. We also suggest that micrometeorites, which are not included in this model, dominate the aerosol properties in the upper stratosphere above about 30 km.

Application of an enhanced PAH growth model to soot formation in a laminar coflow ethylene/air diffusion flame

A recently developed chemical kinetic scheme for C 2 fuel combustion with PAH growth has been implemented in a parallelized coflow flame solver. The reaction mechanism has been developed to include almost all reasonably well-established reaction classes for aromatic ring formation and soot particle precursor molecular weight growth. The model has recently been validated for zero- and one-dimensional premixed flame systems [N.A. Slavinskaya, P. Frank, Combust. Flame 156 (2009) 1705–1722] and has now been updated and extended to a sooting ethylene/air diffusion flame in the coflow geometry. Updates to the mechanism reflect the latest advances in the literature and address numerical stiffness that was present in diffusion flame systems. The chemical kinetic mechanism has been coupled to a sectional aerosol dynamics model for soot growth, considering PAH-based inception and surface condensation, surface chemistry (growth and oxidation), coagulation, and fragmentation. The sectional model predicts the soot aggregate number density and the number of primary particles per aggregate in each section, so as to yield information on particle size distribution and structure. Flame simulation data for the present mechanism is compared to data computed using two other reaction schemes [J. Appel, H. Bockhorn, M. Frenklach, Combust. Flame 121 (2000) 122–136; N.M. Marinov, W.J. Pitz, C.K. Westbrook, A.M. Vincitore, M.J. Castaldi, S.M. Senkan, Combust. Flame 114 (1998) 192–213]. The computed data are also compared to numerous experimental data sets. Whereas the fuel oxidation chemistry in all three mechanisms are essentially the same, the PAH growth pathways vary considerably. It is shown that soot concentrations on the wings of the flame (where soot formation is dominated by surface chemistry) can be predicted with two of the three mechanisms. However, only the present mechanism with its enhanced PAH growth routes can also predict the correct order of magnitude of soot volume fraction in the low-sooting, inception-dominated, central region of the flame. In applying this chemical mechanism, the parameter α, which describes the portion of soot surface sites that are available for chemical reaction, has been reduced to a theoretically acceptable range, thus improving the quality of the model.

Formation of secondary aerosols over Europe: comparison of two gas-phase chemical mechanisms

Abstract. The impact of two recent gas-phase chemical kinetic mechanisms (CB05 and RACM2) on the formation of secondary inorganic and organic aerosols is compared for simulations of PM2.5 over Europe between 15 July and 15 August 2001. The host chemistry transport model is Polair3D of the Polyphemus air-quality platform. Particulate matter is modeled with a sectional aerosol model (SIREAM), which is coupled to the thermodynamic model ISORROPIA for inorganic species and to a module (MAEC) that treats both hydrophobic and hydrophilic species for secondary organic aerosol (SOA). Modifications are made to the gas-phase chemical mechanisms to handle the formation of SOA. In order to isolate the effect of the original chemical mechanisms on PM formation, the addition of reactions and chemical species needed for SOA formation was harmonized to the extent possible between the two gas-phase chemical mechanisms. Model performance is satisfactory with both mechanisms for speciated PM2.5. The monthly-mean difference of the concentration of PM2.5 is less than 1 μg m−3 (6%) over the entire domain. Secondary chemical components of PM2.5 include sulfate, nitrate, ammonium and organic aerosols, and the chemical composition of PM2.5 is not significantly different between the two mechanisms. Monthly-mean concentrations of inorganic aerosol are higher with RACM2 than with CB05 (+16% for sulfate, +11% for nitrate, and +10% for ammonium), whereas the concentrations of organic aerosols are slightly higher with CB05 than with RACM2 (+22% for anthropogenic SOA and +1% for biogenic SOA). Differences in the inorganic and organic aerosols result primarily from differences in oxidant concentrations (OH, O3 and NO3). Nitrate formation tends to be HNO3-limited over land and differences in the concentrations of nitrate are due to differences in concentration of HNO3. Differences in aerosols formed from aromatic SVOC are due to different aromatic oxidation between CB05 and RACM2. The aromatic oxidation in CB05 leads to more cresol formation, which then leads to more SOA. Differences in the aromatic aerosols would be significantly reduced with the recent CB05-TU mechanism for toluene oxidation. Differences in the biogenic aerosols are due to different oxidant concentrations (monoterpenes) and different particulate organic mass concentrations affecting the gas-particle partitioning of SOA (isoprene). These results show that the formulation of a gas-phase chemical kinetic mechanism for ozone can have significant direct (e.g., cresol formation) and indirect (e.g., oxidant levels) effects on PM formation. Furthermore, the incorporation of SOA into an existing gas-phase chemical kinetic mechanism requires the addition of reactions and product species, which should be conducted carefully to preserve the original mechanism design and reflect current knowledge of SOA formation processes (e.g., NOx dependence of some SOA yields). The development of chemical kinetic mechanisms, which offer sufficient detail for both oxidant and SOA formation is recommended.

A new flexible multicomponent model for the study of aerosol dynamics in the marine boundary layer

A new sectional aerosol dynamics model, MAFOR, was developed with the focus to study nucleation in the marine boundary layer. Novel aspects of the model are (1) flexibility in the treatment of gas phase chemistry, (2) treatment of liquid phase chemistry, which can be extended according to needs and (3) simultaneous calculation of number and mass concentration distributions of a multicomponent aerosol as functions of time. Comparison with well-documented aerosol models (MONO32 and AEROFOR), a comprehensive data set on gas phase compounds, aerosol size distribution and chemical composition obtained during theA OE-96 (Arctic Ocean Expedition, 1996)was used to evaluate the model. Dimethyl sulphide decay during advection of an air parcel over the Arctic pack ice was well captured by the applied models and predicted concentrations of gaseous sulphuric acid and methane sulphonic acid range up to 1.0 × 106 cm−3 and 1.8 × 106 cm−3, respectively. Different nucleation schemes were implemented in MAFOR which allow the simulation of new particle formation. Modelled nucleation rates from sulphuric acid nucleation via cluster activation were up to 0.21 cm−3 s−1 while those from ion-mediated nucleationwere below 10−2 cm−3 s−1. Classical homogeneous binary and ternary nucleation theories failed to predict nucleation over the central Arctic Ocean in summer.

Modeling of Oxidation-Driven Soot Aggregate Fragmentation in a Laminar Coflow Diffusion Flame

In this study, three different oxidation-driven soot aggregate fragmentation models with 1:1, 2:1, and 10:1 fragmentation patterns are developed and implemented into a laminar coflow ethylene/air diffusion flame, together with a pyrene-based soot model and a sectional aerosol dynamics model. It is found that the average degree of particle aggregation (n p ) in the soot oxidation region is not correctly predicted if oxidation-driven aggregate fragmentation is neglected; whereas the incorporation of aggregate fragmentation significantly improves the n p prediction in the soot oxidation region. Similar results are obtained using the 1:1 and 2:1 fragmentation patterns. However, as the pattern ratio increases to 10:1, appreciable difference in the predicted n p is observed. As the pattern ratio becomes larger, the fragmentation effect diminishes and the predicted n p approaches that of the original model neglecting fragmentation.