Cloud Processing of Gases and Aerosols in Air Quality Modeling

The representations of cloud processing of gases and aerosols in some of the current state-of-the-art regional air quality models in North America and Europe are reviewed. Key processes reviewed include aerosol activation (or nucleation scavenging of aerosols), aqueous-phase chemistry, and wet deposition/removal of atmospheric tracers. It was found that models vary considerably in the parameterizations or algorithms used in representing these processes. As an emerging area of research, the current understanding of the uptake of water soluble organics by cloud droplets and the potential aqueous-phase reaction pathways leading to the atmospheric secondary organic aerosol (SOA) formation is also reviewed. Sensitivity tests using the AURAMS model have been conducted in order to assess the impact on modeled regional particulate matter (PM) from: (1) the different aerosol activation schemes, (2) the different below-cloud particle scavenging algorithms, and (3) the inclusion of cloud processing of water soluble organics as a potential pathway for the formation of atmospheric SOA. It was found that the modeled droplet number concentrations and ambient PM size distributions were strongly affected by the use of different aerosol activation schemes. The impact on the modeled average ambient PM mass concentration was found to be limited in terms of averaged PM 2.5 concentration (~a few percents) but more significant in terms of PM 1.0 (up to 10 percents). The modeled ambient PM was found to be moderately sensitive to the below-cloud particle scavenging algorithms, with relative differences up to 10% and 20% in terms of PM 2.5 and PM 10 , respectively, when using the two different algorithms for the scavenging coefficient (Λ) 568 corresponding to the lower and upper bounds in the parameterization for Λ. The model simulation with the additional cloud uptake and processing of water-soluble organic gases was shown to improve the evaluation statistics for modeled PM 2.5 OA compared to the IMPROVE network data, and it was demonstrated that the cloud processing of water-soluble organics can indeed be an important mechanism in addition to the traditional secondary organic gas uptake to the particle organic phase.


Introduction
Clouds play an active role in the processing and cycling of chemical tracers in the atmosphere.Gases and aerosol particles can enter cloud droplets through absorption/condensation (of soluble gases) and activation and impact scavenging (of aerosol particles) [1,2].Once inside the cloud droplets these tracers can dissolve, dissociate, and undergo chemical reactions [3].In particular, it has long been recognized that a large portion of the atmospheric particulate sulphate, which contributes to a significant fraction of the total PM mass, is produced in cloud via aqueous-phase oxidation.Past studies have estimated that the aqueous-phase oxidation reactions account for up to 80% of the total production of sulphate globally and are responsible for roughly 50% of the global sulphate burden (e.g., [4,5]).At a regional scale, about 30% of the ambient particulate sulphate concentration in eastern North America can be attributed to the aqueous-phase oxidation production (e.g., [6]).On the other hand, clouds also act as an effective sink for gases and aerosols through precipitation and wet deposition.As a large fraction of clouds globally are non-precipitating [7,8], a significant portion of the cloud processed tracers, often physically and chemically modified, will be released back to the atmosphere.Another aspect of cloud processing is its role in the vertical transport and redistribution of tracers, through the updraft and downdraft often associated with convective clouds [9][10][11] and through precipitation evaporation which brings the tracers processed in cloud to lower levels (closer to the ground).
Cloud processing of gases and aerosols has always been an important aspect in air quality modeling.In the late 1980s and early 1990s, relating to acid rain issues, there were numerous modeling studies investigating the role of cloud processing in producing atmospheric ambient sulphate (e.g., [8,[12][13][14][15][16][17]).There have been a number of studies dedicated to the role of cloud chemistry processing on tropospheric photochemistry (e.g., [18][19][20]).In the context of cloud-aerosol interaction and climate change issues, modeling studies (mostly general circulation models, or GCMs) have been conducted to examine the impact of cloud processing on the global sulphur cycle and on aerosol direct and indirect radiative forcing (e.g., [4,[21][22][23][24][25][26]).More recently, there has been an increasing interest in looking at aerosol-cloud interaction in the context of feedback to dynamics and weather systems using fully coupled models (e.g., [27,28]).There is a wide range of sophistication and completeness in the representation of the various cloud processes in these models.Cloud processing of gases and aerosols remains to be a major source of uncertainty in current air quality models.
Just as particulate sulphate contributes to a significant fraction of the total atmospheric PM mass, particulate organic matter is also a major contributor [29].In recent years, the air quality modeling community has noted the poor performance, shared by almost all models, in predicting atmospheric organic aerosol mass (OA).Modeled OA, with the traditional secondary organic aerosol (SOA) formation mechanism, mostly based on dry or nearly dry smog chamber data, tends to be much lower than the observations.Recently there is increasing evidence that reactions involving organic compounds in the aqueous phase may be an important pathway in the formation of SOA in the atmosphere.Ambient oxygenated organic aerosols with higher organic mass to organic carbon (OM/OC) ratios than the secondary organic aerosol (SOA) formed in smog chambers from volatile organic compound (VOC) precursors have been observed [30], suggesting that our current understanding of ambient organic aerosol chemical processes may be incomplete.It has been postulated that cloud and wet aerosol processing of water soluble organic gases could be a missing source of SOA.The water soluble organic gases are very likely of secondary formation (e.g., glyoxal, methylglyoxal, glycolaldehyde).In traditional SOA formation, precursors produce gas-phase semi-volatile products, in only a few oxidation steps, which can then partition to the organic phase.In contrast, aqueous SOA precursors tend to be smaller molecules because they must be water soluble.When oxidation occurs on molecules with low carbon number, the products can have very high OM/OC ratios (e.g., oxalic acid OM/OC = 3.8; oligomers of methylglyoxal, glycolaldehyde, phenol OM/OC = 2.0 − 2.6).The aqueous phase is unique compared to the gas phase as it concentrates water soluble species.Water soluble aldehyde species enter hydrolysis equilibria creating diols.These diol species react with OH to form peroxy radicals which can decompose efficiently to form an acid group.The OM/OC ratio of the water soluble organic aerosol fraction in aged ambient aerosol has been measured in several studies and ranges from 1.6 to 2.5 ([31-33]).The OM/OC ratio for ambient aerosol has been analyzed with the aerosol mass spectrometer (AMS) data using positive matrix factorization (PMF).The AMS PMF hydrocarbon-like organic factor (HOA) and lowest volatility oxygenated organic factor (OOA-1) factors typically span the lower and upper limit and are typically ~1.3 and ~2.4,respectively [30].Overall, aqueous oxidation chemistry can help explain the high OM/OC ratios for the water soluble organic fraction of aerosol.
Furthermore, ambient particle size distributions for organic species have been found in the droplet-mode size range (centered at diameter 0.7 µm), similar to the droplet-mode sulphate size range [34,35].Sorooshian et al. [36] measured cloud droplet residuals and observed a strong correlation between oxalate and sulphate.Also, Crahan et al. [37] measured two and three times more air equivalent concentrations of sulphate and oxalate in cloud water as compared to ambient below cloud air, respectively.Several modeling studies have simulated the formation of low-molecular weight dicarboxylic acids from in-cloud chemistry [38][39][40].Current knowledge indicates that dicarboxylic acids are only produced in the aqueous phase.Sorooshian et al. [36] obtained a good agreement for the relative production of sulphate to oxalate from in-cloud chemistry between the calculation using a cloud parcel model (constrained by measurements of liquid water content and precursor VOCs) and observed ratios of sulphate to oxalate.Although the mechanisms for SOA formation via aqueous phase are now being explored, there have been a limited number of modeling studies (at regional scale) attempting to assess the impact of these potential additional SOA production pathways on model predictions [39,41,42].In-cloud processing of organic gases also has the potential to alter the vertical distribution of ambient OA as cloud processing occurs at cloud level.Using a regional chemical transport model, Carlton et al. [41] showed large improvement in OA concentration bias and correlation for a flight designed to investigate clouds.
In this paper we review the various approaches and algorithms in representing cloud processing of gas and aerosols in some of the state-of-the-art regional air quality models in North America and Europe.The review will be focused on some of the key processes, such as aerosol activation (Section 2), aqueous-phase chemistry, including the potential pathways for SOA formation via aqueous phase (Section 3), and wet deposition (Section 4).We will explore the variability in the approaches taken in these current air quality models in each of these areas, and assess, wherever possible, the impact of the different approaches on the modeled regional air pollutants based on existing work and sensitivity model runs conducted for this study.Conclusions and recommendations are provided at the end.

Aerosol Activation (Nucleation Scavenging)
The cloud processing of gases and aerosols starts from aerosol activation.Aerosol particles, particularly those that are water-soluble or that contain water-soluble components, can serve as cloud condensation nuclei (CCN).Under suitable supersaturation conditions, condensation of water vapor on these CCN can proceed rapidly and lead to the formation of cloud droplets, i.e., aerosol activation process (or droplet nucleation).Aerosol activation is a key process in aerosol-cloud interaction.It controls the cloud droplet number concentration, which has direct implication on cloud optical properties and precipitation formation; it can affect droplet acidity which may in turn affect the uptake of gases and aqueous phase reactions; it also determines where the aerosol mass addition, due to in-cloud production (e.g., of sulphate), will reside after cloud evaporation and, hence, the cloud processed aerosol size spectrum, which will affect aerosol activation in subsequent cloud cycles (e.g., [43,44]).The modification of aerosol size spectrum will have an impact on other aerosol dynamic and physical processes, such as coagulation, sedimentation, and removal (dry and wet) due to the size dependant nature of these processes.
It has long been known that the nucleation of cloud droplets at the base of clouds is controlled primarily by the aerosol particles that behave as CCN and by the rate of cooling of the air parcel [3,[45][46][47].The activation theory, commonly expressed in a form of the Köhler equation is well established [3,47].One can determine, based on this equation, the equilibrium supersaturation needed for a CCN of given size and physical and chemical properties (e.g., surface tension and solubility) to reach its critical size for growth by water condensation and continue to grow beyond this point even if the supersaturation decreases, i.e., for this CCN to be activated.In the real world situation for a population of aerosol particles the competition between rate of cooling of the rising air parcel and the rate of water condensation on particles results in a maximum supersaturation within the cloud several meters above the base.Under equilibrium, the particles with their critical supersaturation lower than the maximum supersaturation will be activated while the others will remain as interstitial particles, although there are conditions when kinetic limitations may prevent some of the particles from growing to their critical sizes even with their critical supersaturation lower than the maximum supersaturation in the updraft [48,49].The critical supersaturation of a particle is largely influenced by the amount of material it contains that will dissolve in the solution droplet and the degree to which it dissociates (e.g., [3]).Because of their generally lower hygroscopicity in water, organic components appear to have a reduced impact on the nucleation of cloud droplets relative to the inorganic component in terms of the Raoult effect (e.g., [50]).Organic material may enhance activation by reducing the surface tension (e.g., [51][52][53]) or inhibit activation by lowering the rate of accommodation of water molecules by the droplet surface (e.g., [54,55]).Readers are referred to McFiggans et al. [56] for an in-depth review on the effect of physical and chemical aerosol properties on warm cloud droplet activation.
The representation of aerosol activation in 3-D large scale (global and regional) models however is challenging due to the scale difference between the actual process and the model spatial and temporal resolution.It is impractical, if not impossible, to carry out an explicit first principle calculation of aerosol activation in these models.Rather, the process will have to be parameterized.The parameterization varies from primarily empirically based, e.g., linking droplet number concentration to either aerosol mass or number concentration (e.g., [57,58]), to more mechanistically or physically based, often relying on a parcel model concept (e.g., [59,60]).Most of the attention on the representation of aerosol activation and droplet spectrum has been within the climate modeling community as it poses a critical uncertainty in assessing the aerosol indirect effect on climate forcing.Studies have shown that there can be significant variability in modeled droplet number concentrations from these different activation parameterizations under certain conditions (e.g., [24,25]).There is a tendency within the climate modeling community to move away from the more empirically based parameterizations to more physically based mechanistic approaches for representing the aerosol activation processes.On the other hand, these mechanistic approaches have not reduced the uncertainty range of the indirect aerosol effect due to new uncertainties in characterizing size-dependent aerosol mass, composition, and mixing state, and the required parameters for the mechanistically based algorithms [24][25][26]61,62].Within the air quality modeling community, the attention so far has mainly been focused on the prediction of aerosol mass (e.g., PM 2.5 and PM 10 mass) rather than the prediction of aerosol size distribution.However, as alluded to earlier, aerosol activation (or nucleation scavenging), which primarily affects aerosol size distribution, can also have an impact on aerosol mass: it directly affects the wet removal of aerosols in the case of precipitation production ("rain-out", see Section 4) and indirectly the dry removal of aerosols through the impact on the size distribution of cloud-processed aerosols.
Table 1 summarizes the various representations of aerosol activation (or nucleation scavenging) in some of the current regional air quality models.It is not intended to be an exclusive or exhaustive list but these models are representative of the current state of regional air quality modeling in North America and Europe.Amongst the models listed, only a few have a specific representation of the aerosol activation process.In part, how this process is represented in a model is constrained by its aerosol size representation.For example, it is much more difficult to consider partial activation of a mode with a modal representation of aerosol size spectrum without mode splitting-the activation and the subsequent in-cloud production of secondary aerosol tend to modify aerosol size spectrum creating a so-called "Hoppel minimum" [63,64].With a sectional approach, the division of size sections (or resolution) may also dictate the sophistication of the representation of the activation process.Zhang et al. [65] conducted a study dedicated to the impact of aerosol size representation on modeling aerosol-cloud interactions, and demonstrated the limitation of modal approach in representing aerosol size distribution when aqueous-phase chemistry is included.The treatment for aerosol activation in these models is widely varied, from explicit mechanistic representation (e.g., GATOR) to simply assuming all aerosol particles are incorporated in droplets in presence of clouds (e.g., the 2-bin version of CAMx).The GATOR model [66][67][68][69] has the most detailed treatment of the process amongst all, by integrating explicit growth (by water condensation) equations for size resolved (sectional) aerosols.A mechanistic parameterization of Abdul-Razzak and Ghan for multi-mode [59] and sectional [70] aerosols is used in WRF-CHEM [71][72][73].AURAMS [74][75][76], in its current version, employs an empirical parameterization simply relating cloud droplet number concentration to that of aerosol particles with an upper limit of 375 cm −3 [77]; the critical diameter for activation is obtained by counting down aerosol numbers from the top end of the size bins until the total droplet number is reached [74].CMAQ [78,79] with its modal representation of aerosol size distribution assumes all particles in accumulation mode are activated while the Aitken mode particles remain as interstitial [80].For CAMx ( [81,82]) a fixed activation diameter of 0.7 µm is imposed when using a 10-bin configuration while in a 2-bin configuration all aerosol particles are assumed to exist in cloud water in presence of cloud (resolved cloud only).There is no specific representation of aerosol activation in STEM [83,84], it is assumed that once the cloud forms the number of sulphate aerosols activated is the same as the number of cloud droplets [85].For most of the current European regional air quality models (CHIMERE, COSMO-MUSCAT, EMEP), aerosol activation is only considered as part of the wet removal by precipitation (or in-cloud scavenging) and represented by a first order scavenging rate, parameterized by precipitation rate and an experimentally estimated scavenging efficiency for Aitken, accumulation, and coarse mode, respectively [86,87].The LOTOS-EUROS model does not represent nucleation scavenging in its present version [88].Bulk, including primary pathways of sulphate and nitrate production in cloud water, equilibrium for mass transfer between gas and aqueous phase [14], based on Chameides and Davis [98], Chameides [99], and Jacob [100].
Sulphate produced in cloud due to aqueousphase oxidation is re-injected into gas phase as hydrometeors evaporate [14], which can then be partitioned to aerosol phase, e.g., using SCAPE-II [84].To illustrate the impact of aerosol activation on modeled regional aerosols, sensitivity runs using AURAMS were carried out for 1 July-31 August 2004, employing two different aerosol activation schemes: the empirical parameterization of Jones et al. [77] as in the current AURAMS version (see Table 1) and a newly implemented mechanistically-based Abdul-Razzak and Ghan scheme [70].In contrast to the Jones scheme, which links droplet number concentration empirically to aerosol number concentration only, the Abdul-Razzak and Ghan scheme establishes a parameterized relationship (based on the Kohler theory) between the maximum supersaturation reached by an air parcel (under given updraft) and an effective critical supersaturation which can be evaluated from the critical supersaturation of the individual sections/bins.For this particular implementation, the standard deviation of updraft, used as the characteristic updraft in the calculation of aerosol activation [113,114], is parameterized as a function of liquid water content [115], and the aerosol wet size at critical supersaturation is used as a first guess for evaluating droplet size dependent parameters.A separate manuscript is under preparation on more extensive testing and evaluation (with available aircraft and satellite observations) of the aerosol activation schemes in AURAMS.Other processes represented in the off-line AURAMS CTM include emissions from surface and from elevated sources, horizontal and vertical advection, vertical diffusion, gas-phase, aqueous-phase, and inorganic heterogeneous chemistry, secondary organic particle formation, dry and wet deposition, and particle nucleation, condensation, coagulation, and sedimentation [74][75][76].

WRF-CHEM
The simulations were carried out at two horizontal resolutions: 42-km, over a domain covering most of the continental US and Canada, and 15-km, over a smaller domain covering eastern U.S. and south-eastern Canada.The simulation at 42-km resolution provided time-varying boundary conditions for the 15-km resolution simulation.The period coincides with one of the major air quality field campaign conducted over eastern North American and the North Atlantic (ICARTT 2004, [116]) and was chosen for a comparative model evaluation study [117].The details of model set-up are given in [117].
Figure 1 shows the modeled droplet number concentration for a 24-h period (10 August 2004) at the 1235 m model level from the base run (with Jones scheme) and the sensitivity run (with Abdul-Razzak & Ghan scheme).The differences in modeled droplet number concentration from the two different aerosol activation schemes are evident.One immediate observation from Figure 1 is that the Jones scheme seems to result in a more homogeneous droplet number concentration field, which is in part due to its capping at 375 cm −3 , while the Abdul-Razzak & Ghan scheme results in significant higher peak values and more in-homogeneity corresponding to the variation in updraft.Although the aim here is to demonstrate the variability in modeled droplet number concentration due to different activation schemes rather than evaluation, the cap imposed by the Jones scheme may be too low, particularly under continental polluted conditions, as much higher cloud droplet number concentrations have been reported from both aircraft and surface-based high-elevation measurements (e.g., [118,119]).As the primary concern in air quality modeling community has been mass concentration of particulate matter (e.g., PM 2.5 ), modeled surface concentration of sulphate 2.5 and sulphate 1.0 from the two activation schemes are compared in Figure 2. As seen, while the difference in aerosol activation has limited impact on modeled sulphate 2.5 (a few percents), the impact is considerably more significant on modeled sulphate 1.0 (up to 10 percents).This is due to the shift in the activation diameter from the difference in the activation scheme: the Abdul-Razzak & Ghan scheme tends to activate more aerosols (i.e., smaller activation radius) compared to the Jones scheme causing the cloud processed aerosol size distribution to shift towards smaller sizes.This is illustrated in Figure 3 which shows the modeled averaged size distribution of sulphate mass at two grid locations coinciding with two monitoring sites (Dolly Sods in West Virginia and Presque Isle in Maine) within the IMPROVE network.It is seen that, compared to the Jones scheme used in the base-case simulation, the use of Abdul-Razzak & Ghan scheme leads to an increase in mass at smaller sizes and a decrease at larger sizes with the switchover at around 0.5 µm.With the increase and decrease in mass over the size spectrum largely cancelling out each other, the overall impact on the bulk aerosol mass is considerably reduced.The simple sensitivity test shown above provides an indication of the kind of variability in modeled droplet number concentrations one can expect from how the aerosol activation is represented in the model.This is of particular importance for looking at the feedback to radiation through cloud and aerosol optical properties and to cloud microphysics, e.g., precipitation production and subsequent aerosol activation, in fully coupled models (see Zhang [120] for a detailed review of online coupled meteorology and chemistry models; note that while AURAMS is an off-line model, i.e., it does not include the feedback from aerosol to cloud microphysics and dynamics, the change in aerosol size distribution as a result of change in droplet number may affect the diagnosed supersaturation in subsequent aerosol activation).It also demonstrates the strong influence on modeled aerosol size spectrum from the activation parameterization, although the impact on overall modeled bulk aerosol mass is less significant.It needs to be pointed out, however, that the sensitivity shown here could well be influenced by the model parameterizations of the various relevant processes (such as coagulation, wet and dry deposition, and gravitational settling, which may be impacted by the difference in modeled aerosol size distribution due to different activation schemes) and by the aerosol size resolution used (for example, greater impact on aerosol size distribution from the different activation schemes may be expected with higher size/bin resolution).In addition, it is important to note that, while the aerosol activation seems to have limited impact on PM 2.5 mass from this sensitivity demonstration, its impact on PM 1.0 is much more significant; the latter is relevant to ambient measurements of aerosols as some of the current continuous speciated PM measurement technologies experience significant decrease in transmission efficiency for particles greater than 1 µm (e.g., aerosol mass spectrometer: [121]).Finally, the variability tested here (the two different activation schemes) does not cover the full range of variability in the parameterization of the various factors affecting aerosol activation, such as influence of organic surfactants (e.g., [122,123]), updraft velocity (e.g., [113,123]), and entrainment (e.g., [124]).

Aqueous-Phase Chemistry
Aqueous-phase chemistry usually refers to several processes, including mass transfer between the gas phase and aqueous phase (solution droplets), i.e., gas phase diffusion to/from aqueous droplets, surface accommodation, hydrolysis/dissociation of certain dissolved species, diffusion in aqueous solution, and chemical reactions amongst various species in aqueous solution [3].The nature of these processes and the pollutants involved may be quite different for concentrated solution, e.g., aqueous aerosols or haze particles, and dilute solution, e.g., cloud droplets.Since this article deals with cloud processing of gases and aerosols, we will focus on the aqueous-phase chemistry in dilute solution droplets (or close to ideal solution) [125].
Driven by acid deposition issues, a significant body of work in the 80's has been dedicated to the development and testing of reaction mechanisms for atmospheric aqueous-phase chemistry in cloud and rain including the mass transfer aspect (e.g., [96,[98][99][100][126][127][128][129][130]).The mechanisms generally cover oxygen-hydrogen chemistry, carbonate chemistry, chlorine chemistry, nitrite and nitrate chemistry, organic chemistry focused on formaldehyde oxidation, and sulphur chemistry [96].By conducting a formal sensitivity analysis on a comprehensive atmospheric aqueous-phase chemistry mechanism, Pandis and Seinfeld [96] found that, by comparing the sensitivity of various aqueous-phase species concentrations with respect to each of the individual reactions in the mechanism, the dominant reaction in the system is the oxidation of S(IV) to S(VI) by H 2 O 2 followed by other S(IV) oxidation pathways, e.g., catalyzed oxidation and oxidation by O 3 .Reactions that produce or consume H 2 O 2 are important as they influence the S(IV)-to-S(VI) conversion, and aqueous phase reactions are also important pathways for the production of formic acid through the reaction of dissolved formaldehyde with OH.The mechanisms developed in the 80s formed the basis for the aqueous-phase chemistry mechanism used in the current air quality models.
In more recent years, the atmospheric chemistry community has again turned its attention to aqueous-phase processes to look for possible pathways contributing to atmospheric SOA formation (e.g., [34,131]).There is increasing evidence that cloud processing of water soluble organics may indeed be playing a role in the formation of atmospheric SOA.These missing pathways from the current atmospheric models may at least partly explain the under-prediction of organic aerosol mass commonly found with these models (e.g., [132][133][134][135][136]).
In the following, we will first review the representations of cloud aqueous-phase chemistry in current regional air quality models; we will then discuss some of the recent studies on the possible cloud aqueous-phase oxidation pathways leading to atmospheric SOA formation and some preliminary assessments of their potential impact on modeled aerosol mass on a regional scale.

Representation of Aqueous-Phase Chemistry in the Current Regional Air Quality Models
The last two columns in Table 1 summarize the treatment of aqueous phase chemistry and its feedback to atmospheric gases and aerosols in some of the current regional air quality models.As seen, all models have some representation of in-cloud sulphate production through aqueous-phase oxidation.The various representations can be grouped into three categories as far as aqueous-phase chemistry mechanism is concerned: (1) models employing more complete aqueous-phase chemistry mechanism, such as GATOR, STEM, CAMx (with CMU aqueous-phase chemistry) and WRF-CHEM (PNNL version); (2) models employing more condensed mechanisms focusing on aqueous-phase sulphur chemistry only, such as AURAMS, CAMx (with RADM aqueous-phase chemistry), CMAQ, CHIMERE, and EMEP unified model; and (3) models employing a highly parameterized representation of sulphur conversion in aqueous-phase, such as COSMO-MUSCAT and LOTUS-EUROS.The more complete mechanisms include not only the sulphur chemistry but also reactions involving other groups, e.g., the oxygen-hydrogen, chlorine, carbonyl.The ADOM (acid deposition and oxidant model, [13]) aqueous-phase chemistry mechanism, used by AURAMS, and the RADM (regional acid deposition model, [12]) aqueous-phase chemistry mechanism, used in CMAQ and CAMx, are very similar.Both include the mass transfer and aqueous-phase dissociation reactions (for SO 2 , O 3 , H 2 O 2 , organic peroxides, HNO 3 , NH 3 , and CO 2 ) and S(IV)-to-S(VI) oxidation by H 2 O 2 , O 3 , organic peroxides, and O 2 (catalyzed by iron and manganese).Comparing with the more complete mechanism, this sulphur oxidation-focused mechanism does not resolve the additional sources and sinks of peroxy radicals and peroxides in aqueous phase (through the oxygen-hydrogen chemistry) which will affect aqueous-phase sulphur oxidation.However if the system is not oxidant-limited this omission may not have as much impact as far as sulphur oxidation is concerned.On the other hand, the impact on tropospheric photo-oxidant chemistry from the aqueous-phase oxygen-hydrogen chemistry and the aqueous-phase sink for formaldehyde can be significant [18,20].The approach taken by LOTUS-EUROS and COSMO-MUSCAT is to represent the aqueous-phase sulphur oxidation by a single first order reaction with a fixed base rate constant modulated by cloud cover and relative humidity [88,111].Although this approach may work well in some cases, it may not be as robust given the use of a single fixed base rate constant which in effect will impose a certain level of oxidants which may or may not be representative of any given region and/or season.
Almost all models use a bulk representation for aqueous-phase chemistry.The variable size resolution model (VSRM, [92]) is used in CAMx (as an option) and WRF-CHEM, by which, depending on a set of criteria, the aqueous-phase chemistry is solved either in bulk or in two bins (sections) separating the droplets formed on larger, more alkaline particles from those formed on smaller, more acidic particles.Past studies (e.g., [64,[137][138][139]) have shown enhanced sulphate production from a size-resolved representation of aqueous-phase chemistry compared to a bulk representation.The difference is mainly due to the inhomogeneity of pH amongst cloud droplets and the strong pH dependency of the aqueous-phase SO 2 -O 3 reaction rate.Kreidenweis et al. [64] showed that, for a parcel-model comparison case, nearly three times more sulphate is produced through the ozone oxidation pathway from the size-resolved models than from the bulk models.However, sulphate production through the H 2 O 2 oxidation pathway is less from the size-resolved models, and the net enhancement in sulphate production from the size-resolved calculation is about 15% for the case considered.Barth [140] looked into the importance of droplet size representation from the cloud photochemistry perspective and found that, for species such as formaldehyde and formic acid, model prediction is very sensitive to the droplet size representation since the fate of these species in aqueous phase is highly dependent on pH.
Other variability in aqueous-phase chemistry representation amongst the various models includes the treatment of mass transfer between gaseous and aqueous phase and cloud water pH.Almost all the models with explicit aqueous-phase chemistry mechanism, whether more complete or more condensed, assume Henry's law equilibrium between gaseous and aqueous phase concentrations for species involved in mass transfer processes.In contrast, AURAMS employs a kinetic approach for mass transfer between gas and aqueous phase (except for O 3 and CO 2 ) using the Fuchs and Sutugin [141] formulation for diffusion coefficients.Under most of the atmospheric conditions, the Henry's law equilibrium can be established very quickly for moderately and less soluble species.For highly soluble species such as HNO 3 the equilibrium partitioning between gas and aqueous-phase may not be achieved easily.As to the cloud water pH, most of the models that do treat the pH-dependent dissociation and oxidation reactions diagnose pH (or hydrogen ion concentration) from charge balance, while EMEP unified model uses a fixed pH of 4.5, a value which may be representative of an average condition found in Europe but may not be representative at a finer resolution (i.e., at a local scale).Even for the diagnosed pH, the model may not have the representation of the full set of base cations and therefore pH may be underestimated (i.e., more acidic) from an incomplete charge balance.
Models do vary in how the aqueous-phase mass production (bulk) feeds back to aerosols after cloud evaporation.Depending on the aerosol size representation in a model, the bulk mass increment due to the aqueous-phase S(IV)-to-S(VI) oxidation may be added to the aerosol size spectrum in various ways.For three of the sectional models (AURAMS, GATOR, and WRF-CHEM PNNL version), the bulk mass increment is distributed to the aerosol size spectrum according to the ratio of the liquid water content in each aerosol size sections to the total (bulk) liquid water content, although each model has a different way of determining the liquid water in each size sections.In AURAMS, the bulk cloud liquid water (provided by the meteorological driver model) is distributed evenly to all activated aerosol particles and the liquid water in each size sections is determined by the aerosol number in each activated or partially activated sections [74]; in GATOR, cloud (bulk) microphysics is solved first for both sub-grid scale and resolved clouds, and the water vapor evaporated/sublimated from the bulk microphysics is then re-condensed/re-deposited onto aerosol particles along with other aerosol-hydrometeor and hydrometeor-hydrometeor interactions to produce aerosol/hydrometeor spectrum [69]; in WRF-CHEM (PNNL version) a Gamma distribution of droplets and a maximum correlation between droplet size and dry particle size are assumed, i.e., the largest activated aerosol particles are associated with the largest droplets, [65].For the CMAQ model with modal representation, the bulk mass due to aqueous-phase oxidation is added to the accumulation mode only, maintaining total number and geometric standard deviation of the mode [80].In STEM, the sulphate produced in cloud due to aqueous-phase oxidation is re-injected into gas phase as hydrometeors evaporate [14], which can then be partitioned to the aerosol phase, e.g., using a thermodynamics module for inorganic gas-particle equilibrium partitioning, SCAPE-II (Simulating Composition of Atmospheric Particles at Equilibrium, [142][143][144]) [84].For those models with very coarse size resolution for the aerosol size spectrum, e.g., two bins (modes), the bulk mass addition may be added to one or two of the bins (modes) depending on the size cut.As discussed in Section 2 above, the modification of aerosol size spectrum from cloud processing (particularly, aerosol activation and sulphate production in cloud) and its representation will have a direct impact on modeled aerosol optical properties (important for direct aerosol radiative forcing) and on aerosol hygroscopicity (important for subsequent cloud processing).The difference in cloud processed aerosol size distribution due to the different model representations will also impact some of the size dependant aerosol microphysical and dynamic processes, which will in turn have an impact on modeled aerosol mass concentrations.In a fully coupled chemistry-dynamic system, the modified aerosol size spectrum (both mass and composition) through cloud processing will feed back to subsequent cloud microphysics which will, in turn, impact aqueous-phase chemistry, although the majority of current regional air quality models do not have the representation of this feedback mechanism.
Although it is outside the scope of this review, it is worth pointing out that the largest uncertainty by far, in modeling cloud aqueous-phase chemistry, may be from modeled cloud, including timing, location, and liquid water content.Through their sensitivity analysis, Pandis and Seinfeld [96] demonstrated that the liquid water content is one of the most important parameters of the aqueous-phase chemistry system.Particularly, the aqueous-phase sulphate mass production is very sensitive to the liquid water content.Modeling clouds in large scale models is extremely challenging due to the spatial and temporal scales of the processes involved and their feedback to the dynamic system.Zhang et al. [145] conducted an evaluation of modeled cloud properties against in-situ measurements for air quality application, and found that the model's ability to reproduce observed clouds depends on the model resolution, the particular microphysics parameterization employed, and the type of clouds modeled.The uncertainty in modeled in-cloud production of aerosol mass is expected to be large due to the uncertainty in the model predicted cloud and liquid water content.

Mechanism for Cloud Processing of Organic Gases
The hypothesis is that some of the gaseous organic species will be taken up by droplets, and some of them will undergo irreversible reactions (oxidations) and form products of low volatility which after cloud droplets evaporate will be partitioned in particle phase.
Similar to aqueous phase inorganic chemistry, the uptake of gaseous organic species by droplets will be affected by a number of processes (steps): gas phase diffusion, surface accommodation (or transfer across the interface), attainment of aqueous-phase equilibrium (solubility and hydrolysis), and possible reactions in aqueous phase.Gas-phase diffusion is usually fast and there is not a vast variation amongst different species.Surface accommodation, described by an accommodation (or sticking) coefficient, is also unlikely to be a rate limiting step (unless the accommodation coefficient is much smaller than 0.01).The remaining key factors affecting the uptake of an organic gaseous species would be solubility and hydrolysis of a given species and its possible reactions in the aqueous phase.Mathematically, the uptake of gaseous organic species may be formulated as follows [141,[146][147][148][149]: where C g is the gas-phase concentration, γ is a parameterized uptake coefficient (E2), A is the cloud liquid surface area per volume of air, v is the gas-phase molecular velocity, K H is the effective Henry's law coefficient, R is the universal gas constant, T is the temperature, a is the aqueous drop radius, D g is the gas-phase diffusion coefficient, k rxn is the pseudo-first order reaction rate coefficient in the aqueous phase, D L is the liquid phase diffusion rate coefficient, and α is the accommodation coefficient.The three terms in E2, from left to right, describe the gas-phase diffusion, surface accommodation, and aqueous-phase reactions, respectively.
Henry's law coefficients describe the solubility or dissolution of an organic (or inorganic) gas in aqueous solution.Hydration equilibrium coefficients can be combined with Henry's law coefficients to form effective Henry's law coefficients (K H ). Table 2 summarizes the K H for some water soluble organic species from existing work as listed in the 4th column of the table.Interestingly, glycolaldehyde has a larger Henry's law coefficient compared to glyoxal; however, glyoxal has a very large hydration coefficient resulting in glyoxal having an order of magnitude larger effective K H than glycolaldehyde.In looking at equation (E2), the 3rd term includes the product of K H , k rxn and D L .However, if one considers that k rxn and D L have physical upper limits, then one can calculate a minimum K H for which uptake becomes important.Gelencser and Varga [149] estimated a threshold effective Henry's law coefficient of 103 M/atm, above which mass transfer to droplets becomes important.A review of the structure-activity relationships to estimate K H [150] based on specific functional groups in a molecule revealed that, in general, monofunctional alcohols, peroxides, ethers, aldehydes, ketones, organic nitrates, nitro compounds, and esters do not meet the threshold of 10 3 M/atm.Aromatic alcohols and organic acids do meet the general rule of thumb.Many bifunctional compounds meet the threshold and need to be further assessed; however, species with ether and ester groups, ether and ketone groups, two ether groups, two ketone groups, two ester groups, two nitrate groups, two halide groups, nitro-aromatics and ketone groups, and nitro-aromatics and alcohol groups did not meet the criteria.All bifunctional species with at least one aldehyde group met the criteria.For trifunctional species, the majority met the criteria with the exception of species with all halide groups.Once an organic species is dissolved in solution, it can undergo irreversible oxidation reactions.For organic species, OH-initiated oxidation is likely the most important pathway [161].Dutot et al. [62] developed a neural network method to estimate aqueous OH-rate coefficients for organic compounds.Their training database was based largely on compilation of rate coefficients in [163].Chemical groups considered include aliphatic bonds, unsaturated bonds and functional groups including alcohols, carboxylic acids, carboxalate ions, ketones, aldehydes, esters, ethers and aromatic functional groups, but not hydroperoxide and nitrate groups because of insufficient data.Also, Monod et al. [164] and Monod et al. [165] developed structure-activity relationships (SAR) for aqueous OH-reactivity, similar to gas-phase SAR methods developed by Atkinson [166].Table 3 lists the OH-initiated aqueous oxidation rate coefficients for several C2 compounds with different functional groups based on Monod et al. [165].An examination of Table 3 shows that the OH reaction rate is faster for ethanol, the unhydrated form of acetaldehyde, and the glyoxalate, CHOCOO-relative to the ethane reaction rate.The hydrated acetaldehyde, hydrated glyoxal, glycolaldehyde all have ratios relative to ethane between 0.2 and 1.0.Acetone, methyl peroxide, acetic acid, oxalic acid dissociation ions have values between 0.01 and 0.2.Thus, on inspection, aldehydes, alcohols and acid groups adjacent to carbon chains do not vastly affect the H-abstraction from the adjacent carbon chain.The non-dissociated dicarboxylic acid, oxalic acid, does have a vastly slower reaction rate relative to ethane, 8.8 × 10 −4 .Also noted from Table 3, anions usually react faster than their corresponding acids with OH because of different mechanisms (electron transfer compared to H-abstraction) [167].Based on recent literature, the following summarizes some of the key atmospheric soluble organic species that can undergo aqueous OH-initiated oxidation with the potential to form SOA:  In cloud droplets, during daylight hours, glyoxal will react with OH to irreversibly form glyoxalic acid which, in turn, can form oxalic acid ( [169]; Table S1 of [125] for reaction mechanism).At cloud relevant glyoxal concentrations (10 −5 to 10 −4 M glyoxal), the predicted mass yields of oxalic acid are 140% related to the change in aqueous glyoxal concentration [125].In wet aerosol, irreversible reactions with OH help produce high molecular weight oligomers (with lower volatility than oxalic acid) via radical-radical reactions. At lower methylglyoxal concentrations typical of cloud concentrations (dissolved methylglyoxal ~µM order), the predicted temporal evolution of methylglyoxal products (oxalic acid, pyruvic acid) and total water soluble organic carbon compared quite well to the measurements in an UV illuminated reaction vessel (aqueous scheme with 19 reactions specific to methylglyoxal oxidation, Table 1 in [170]).After irradiation conditions representative of one cloud cycle ([OH] ~ 4 × 10 −12 M for 10 min), the pyruvic acid mass yield per mass of methylglyoxal reacted was 108%.After 90 min of irradiation, acetic acid reached a maximum concentration, pyruvic acid was decreasing and oxalic acid was increasing.After 210 min of irradiation, oxalic acid reached a maximum, acetic acid was slowly declining and pyruvic acid was a minor product (Figure 2 in [170]).All of these acid species have effective Henry's law coefficients higher than 10 3 M/atm and thus remain, in part, in the cloud phase.As the cloud droplets evaporate, some of the (semi-)volatile species will return to gas phase while other (semi-)volatile species may remain bound in particle phase due to oligomerization reactions (refer to Section 3. .While formaldehyde does affect the aqueous OH budget, it is not likely a critical precursor to in-cloud SOA production. Acetone has an intermediate reactivity and effective Henry's law coefficient and, thus, has the potential to make a small contribution in terms of the in-cloud production of formaldehyde, methylglyoxal, hydroxyacetone and acetic acid [171]. Methyl peroxide has two major in-cloud reaction pathways, namely, OH reaction to form the methyl peroxy radical, CH 3 O 2 , and reaction with HSO 4 − to form SO 4 2− , H + and methanol.Since the cloud phase is depleted in oxides of nitrogen relative to the gas-phase, then the methyl peroxy radicals can recombine to reform organic peroxides.While methyl peroxide does affect the aqueous radical budget, it is not likely a critical precursor to in-cloud SOA production. Liu et al. [172] and Zhang et al. [173] proposed cloud-relevant radical chemistry mechanisms for methacrolein (MACR) and methyl vinyl ketone (MVK) initiated by OH-oxidation.The mechanisms were successful at predicting the first-generation oxidation products of MACR and MVK for laboratory experimental conditions, but high carbon number organic acids were poorly predicted.It was postulated that the decomposition of higher oligomers contributed to the oxalate ion.El Haddad et al. [174] investigated the fate of MACR under cloud conditions and measured an aqueous SOA yield of 2-12% relative to a change in the initial MACR aqueous concentration.These aqueous in-cloud SOA yields are not directly comparable to the cloud-free isoprene SOA yield measured (1-3% relative to a change in gas-phase precursor concentration) in low relative humidity smog chambers [175].The Henry's law coefficients for MVK and MACR are relatively low, 41 M/atm and 6.5 M/atm, respectively ( [153], 298K), and thus it is unlikely that these carbonyl compounds will be significant contributors to atmospheric SOA.The current aqueous yields were also measured at higher initial MACR molarities (mM range typical of wet aerosol) and after long reaction times (22 h).Further laboratory studies are recommended at ambient aqueous MACR molarities.Chemical transport modelling sensitivity studies with and without the aqueous pathway are required to assess its importance on regional atmospheric scales over forest ecosystems.Glycolaldeyde is a water-soluble species with K H > 4 × 10 4 M/atm [161] and has been observed in the µM concentration range in cloud water [176].In the cloud-phase, glycolaldehyde is expected to oxidize to glyoxal and glycolic acid, which can, in turn, oxidize further to form glyoxalic acid and oxalic acid [177].At early times (<25 min), glycolic acid is the dominant product species and SOA yields are 50-60% relative to the change in aqueous glycolaldehyde concentration.After 100 min, oxalic acid is the dominant product and the SOA yield dropped to a relatively constant value of 30% relative to the change in aqueous glycolaldehyde (see Table 1 in [177] for reaction mechanism).Given the large gas-phase yields of glycolaldehyde from MVK (~70% by mass) and large K H , it is possible that glycolaldehyde could contribute to aqueous-phase SOA production.Glycolaldehyde is also efficiently produced in the gas-phase from the oxidation of anthropogenic ethane and is emitted directly from the smouldering stage of biomass burning [178]. Cloud processing of phenolic species (phenol, guaiacol, syringol) has been shown to produce SOA with mass yields near unity relative to a change in aqueous precursor concentration [179].The lifetime of phenols is a couple of hours at midday summertime conditions for both gas-phase and aqueous-phase OH-initiated oxidation.Henry's law coefficient (at cloud temperatures of 5 °C and liquid water content of 0.3 g/m 3 ) predicts a phase partitioning of 3% phenol (K H ~ 5.5 × 10 3 M/atm), 4% guaiacol (K H ~ 7.2 × 10 3 M/atm), and 15% syringol (K H ~ 2.7 × 10 4 M/atm), by mass, in the cloud-phase.The OM/OC ratios measured are in the range 2.2-2.6 [179] which corresponds to the oxidation states observed for the aerosol mass spectrometer lowest volatility factor in the positive matrix factorization analysis (termed OOA-1).The SOA mass spectra show evidence for the formation of dimers and higher oligomers via C-O and C-C linkages of phenoxyl radicals [179].Prominent mass spectra peaks observed in the cloud processing of wood combustion are consistent with the peaks in the phenolic OH-initiated aqueous oxidation product mass spectra.The reader is also referred to a recent, more detailed review focusing on many aspects of aqueous phase SOA chemistry [180].
Chen et al. [39] coupled a near explicit gas-phase chemical mechanism (Caltech Atmospheric Chemistry Mechanism, CACM) to a cloud chemistry mechanism.In their study, 14 C 1 -C 3 organic species are explicitly treated and partitioned between the gas and cloud phase based on individual species properties.In addition, 58 semi-volatile species are lumped into 11 surrogate species of which 4 are considered hydrophilic and reactive enough to participate in aqueous cloud oxidation chemistry (47 aqueous reactions, Table 1 in [39]).Two of these surrogates are formed from aromatic oxidation chemistry and two are from monoterpene oxidation chemistry.Aqueous-phase OH-initiated reactions for these 4 surrogates are considered and each reaction forms one non-volatile product with unity yield.Three of the multi-functional products have organic acid groups and, thus, these species can enter dissociation equilibrium.A 48-h averaged box model simulation for a rural scenario showed an average 32% increase in PM 2.5 SOA due to cloud processing for 3-h per day.A 48-h averaged box model simulation for an urban scenario showed an average 11% increase in PM 2.5 SOA due to cloud processing for 3-h per day.
The most detailed aqueous-phase mechanism (CAPRAM v3.0i, C 1 -C 5 chemistry, [181]) has been reduced to a version potentially applicable to regional chemical transport models, CAPRAM-RED [161].The number of aqueous phase species decreased from 380 in the full mechanism to 130 in the reduced mechanism (40% reduction in computation times).The averaged percentage differences calculated between the two mechanisms were below 5% for most important organic species (ketones, alcohols, aldehydes, organic acids).A 2-D version of the COSMO-MUSCAT modeling system [106] was developed to compare the CAPRAM-RED mechanism to a simpler inorganic only aqueous mechanism for a prescribed cloud case with two different emission scenario tests (urban and remote).The 2-D test simulations with the reduced organic mechanism and the inorganic only mechanism showed large differences in oxidant levels for the urban case (up to a factor of 4 lower gas-phase OH mixing ratio with the reduced organic mechanism) and smaller differences for the remote case.For the remote test case, oxalate and glycolate were the dominant organic acids predicted downwind of the cloud, whereas, for the urban case, oxalate and pyruvate were the dominant organic acids predicted.

Impact of Cloud Processing of Organic Gases on SOA Mass Concentrations at Regional Scales (Upper Limit Assessment)
There are only a few studies which attempted to assess the atmospheric implications of cloud processing of organic gases on a regional scale, all with the use of the Community Multiscale Air Quality (CMAQ) model [39,42,169].Carlton et al. [41] assumed Henry's law equilibrium for glyoxal and methylglyoxal (lumped together) to the cloud-phase followed by aqueous-phase OH-initiated oxidation to form non-volatile SOA mass (4% molar yield).More recent laboratory and box model studies suggest the 4% molar yield is a conservative estimate [125].Nevertheless, Carlton et al. [41] did observe significant improvement in both the SOA bias and model vs. measurement correlation coefficient for comparisons where aircraft measurements were impacted by clouds and model results adequately predicted cloud location.CMAQv4.7 was also run with a continental U.S. domain (36 km grid spacing) over an annual cycle for 2003.When averaged by month across 97 ground sites, the temporal correlation did show a significant improvement, but there still remained a negative model bias for the surface SOA mass concentration during the summer months [42].
The Chen et al. [39] mechanism discussed above was also incorporated into CMAQ and run for a short, 2-day summertime 2004 period, but with the detailed gas-phase Cal-Tech Atmospheric Chemistry Mechanism and 47 reaction aqueous-phase scheme.On average, over an eastern U.S. domain, surface SOA predictions increased by 9% for a 2-day August 2004 period.Results also indicated that the maximum contribution of SOA from irreversible aqueous phase organic reactions was 0.6 µg/m 3 for several one-hour averaged grid-cell locations.
Recently, an aqueous-phase organic oxidation scheme was placed inside the global model, TM4-ECPL, to study the oxalic acid budget [182].An empirical reaction scheme (based on [38,169,183,184]) was employed assuming Henry's law equilibrium for formic acid, acetic acid, glyoxal, methylglyoxal and glycolaldehyde.Model predictions of oxalate at rural locations were consistent with predictions in terms of both bias and correlation.Most oxalate was formed from in-cloud processing rather than wet aerosol; 66% of oxalate was removed by deposition processes and 34% by further aqueous oxidation.The oxalate production originated largely from gas-phase isoprene oxidation.The condensed phase net source of oxalate in conjunction with a global mean turnover time against deposition of about 5 days, maintained oxalate's global tropospheric burden at 0.2-0.3Tg (i.e., 0.05-0.1 Tg-C) which is about 5-9% of model-calculated water soluble organic carbon burden, though the model prediction of total organic carbon was lower by approximately a factor of 2 relative to observations.
To get a broad idea of how the potential cloud processing of organic gases may impact model predicted organic aerosol component on a regional scale over North America, a sensitivity run was conducted using AURAMS.AURAMS uses the ADOM-II gas-phase chemical mechanism [185,186] and three water soluble organic species were selected to partition to the cloudphase: MGLY (C 2 -C 3 dicarbonyl), DIAL (larger dicarbonyls, from aromatic oxidation) and CRES (from aromatic oxidation and emission).A group contribution method to estimate effective Henry's law coefficients at 25 °C [150] yielded K H > 1 × 10 3 M/atm for these three surrogate species (model values listed in Table 2 footnotes).A pseudo-first order uptake was assumed for the calculation of mass transfer to the modeled cloud LWC (assuming 10 µm diameter drops), following (E2).Gas-phase diffusion constants and molecular velocities were calculated at 288 K (in the range D g = 1.0 × 10 −5 to 1.2 × 10 −5 m 2 /s, v = 268-316 m/s).Model K H values were calculated depending on the model grid cell temperature using enthalpy's of vaporization (see Table 2 footnotes).An accommodation coefficient of 0.023 was estimated for MGLY [40], and 0.05 for DIAL and CRES [39].Aqueous reactions were assumed to form non-volatile SOA mass with a yield of unity.A maximum pseudo-first order OHinitiated aqueous reaction rate coefficient of 10 −4 s −1 was assumed and scaled by the cosine of the solar zenith angle.The aqueous SOA was not assumed to oxidize further.These assumptions likely represent an upper limit for cloud organic oxidation chemistry to form SOA. A significant difference between our base case and sensitivity run would suggest that computationally efficient, reduced mechanisms for cloud organic chemistry should be tested and implemented into CTMs, such as the mechanisms proposed by Deguillaume et al. [161] and Lim et al. [125] (see Section 3.2.1 above).A small difference between runs would suggest that organic aerosol produced by cloud uptake and cloud aqueous chemistry is small compared to organic aerosol produced from traditional SOA formation mechanisms.
For the simulation, hourly anthropogenic point, area, and mobile emission files were prepared for the CTM using the 2005 Canadian and 2005 US emissions inventories and version 2.3 of the SMOKE emissions processing system (Sparse Matrix Operator Kernel Emission, [187]).Biogenic emissions were calculated on-line by AURAMS using BEIS version 3.09, the Biogenic Emissions Landcover Database (BELD3) vegetation data set, (30 tree species and 20 crop species used for Canada), and the modeled meteorological fields (temperature and irradiance).Biogenic VOC emissions are speciated into four groups: isoprene (ISOP), monoterpenes (PINE), sesquiterpenes (SESQ) and "other VOCs".For the soluble organics species partitioned to cloud phase in this sensitivity test, isoprene is the sole relevant biogenic precursor species.Its standard emission fluxes (defined at 30 °C and photosynthetically active radiation, PAR, at 1000 μmol/m 2 /s) from oak, populus, maple and fir trees were 26,000, 21,000, 38, and 150 μgC/m 2 /h, respectively, (product of basal emission rate, leaf area index and leaf mass density); the actual emission flux in the model is modulated by prognostic temperature and radiation (irradiance).
Similar to the setup for the AURAMS sensitivity run on aerosol activation in Section 2, a 2-month simulation with the simple implementation of the additional SOA formation from cloud processing of water soluble organics was performed for the period of 1 July to 31 August 2004 at two horizontal grid resolutions (42-km and 15-km).Figure 4a is the modeled 2-month averaged PM 2.5 SOA mass concentration from the base-case simulation, and Figure 4b shows the difference between the sensitivity run (with the cloud-phase non-volatile SOA production) and the base case.It is seen that the modeled PM 2.5 SOA is enhanced quite significantly by the inclusion of in-cloud production of SOA (up to 40%).There is a clear indication of biogenic origin of the modeled SOA, as has been observed in other studies by Chen et al. [39] and Hennigan et al. [188].Shown in Figure 4c is a model-observation comparison of PM 2.5 organic aerosol (OA) component at available IMPROVE network sites in terms of averages over the period.Speciated PM 2.5 measurements from filter samples (24-h, one in 3 days) are available at 54 sites over the simulation period.Table 4 summarizes the PM 2.5 OA performance statistics for the two simulations as compared to the IMPROVE network data.It is evident that the sensitivity run including a simple representation of aqueous organic chemistry mechanisms improves predictions in terms of the bias and root mean square error as well as the slope for PM 2.5 OA (refer to Figure 4c, model vs. measurement OA mass concentration correlation).Although these results need to be taken in with caution as the representation is very simplified, designed to give an upper bound, and subject to many uncertainties, they do imply that the potential impact on modeled SOA from the cloud processing of water soluble organics can be significant.Further study is definitely warranted.

Uncertainties in Mechanisms for Cloud Processing of Organic Gases
A key uncertainty related to forming SOA from cloud processing arises from the current lack of understanding of the processes that occur as cloud droplets evaporate.As water evaporates from cloud droplets, solute concentrations increase which can result in complex radical and non-radical chemistry:  Radical chemistry: Volkamer et al. [189] showed that photochemistry enhanced SOA formation from glyoxal uptake to wet aerosol.They observed that in the presence of UV light, the rate of volume growth for glyoxal was 2-3 orders of magnitude larger as compared to dark conditions.Observed SOA formation rates in laboratory experiments were shown to scale linearly with aqueous phase OH.In evaporating wet particles, the uptake of glyoxal and their subsequent OH-oxidation allows the accumulation of alkyl radicals in the particle phase [190].These alkyl radicals either react with O 2 or recombine with other alkyl radicals.In a drying particle, the concentration of alkyl radicals increase which favors the recombination reactions.This leads to an increase in the molecular size of products ( [125,183]).Guzman et al. [190] also observed that the alkyl radical + O 2 reaction rate coefficients were considerably smaller in the aqueous phase due to hydrogen bonding and/or a H 2 O cage around the alkyl radicals inhibiting O 2 attack.Aqueous particles are believed to have greater OH concentrations compared to large cloud drops due to faster OH diffusion rates to a particle (λ  1/a in first term of E2 where a is the radius), thus small wet particles may have greater alkyl radical production which favors the self-bimolecular termination reactions. Non-radical chemistry: As droplets evaporate, glyoxal and methylglyoxal can undergo self-oligomerization by hemiacetal formation and/or aldol condensation and the sensitivity of this pathway to acidity is still uncertain [189,191,192]. Oligomer formation through amino acid catalyzed and ammonium ion catalyzed aldol condensation [184,193] has been observed in concentrated acetaldehyde solutions representative of wet aerosol. Reactions of glyoxal and methylglyoxal with amino acids, methylamine and (NH 4 ) 2 SO 4 have also been observed in wet aqueous aerosol and evaporating cloud drops.Their results are consistent with rapid imine formation followed by the formation of nitrogen-containing oligomers, methylimidazole and dimethylimidazole products [184,[192][193][194][195][196][197]. The mechanism of organosulphate formation in evaporating drops is uncertain.Organosulphate (sulphate esters) can form by (1) the acid catalyzed esterification reaction of an alcohol (hydrated aldehyde group) with sulfuric acid followed by dehydration or by (2) radical-radical reactions.Recent studies have shown that organosulphate yields for precursor alkenes vary with OH radical concentrations (i.e., sunlight conditions) and provide strong evidence for the radicalradical mechanism [195,198].The OH • radicals can convert (1) bisulphate ion, HSO 4 − , to the SO 4 − radical, (2) sulphuric acid to the HSO 4 radical and, (3) RH (alkenes, aldehydes) to the R alkyl radical.The alkyl radical can then react with either the SO 4 − radical or the HSO 4 radical to produce organosulphates.With a time constant of several hours for radical-radical reactions, these reactions are more plausible in the atmosphere compared to the slower time scales for esterification [170].Fundamental aqueous-phase kinetic rate coefficients are clearly needed for all the above processes under ambient aqueous solution conditions from bulk (reaction-vessel) laboratory experiments. The reversibility of oligomer formation is uncertain in terms of extent of reactions and timescales.
There are several studies that suggest photochemically formed oligomers are irreversible whereas oligomers formed in the absence of radicals are reversibly formed [125,183].Hastings et al. [199] observed oligomer formation and no significant evaporation of dehydrated glyoxal when drying solutions.Loeffler et al. [192] suggest that oligomer formation during cloud evaporation and in wet aerosol is fast compared to evaporation of the surface dehydrated glyoxal and de-oligomerization is slow allowing oligomers to stabilize.Schwier et al. [200] also observed no reversible process of glyoxal oligomers upon dilution over several hours.Other studies suggest the process is only partially reversible.Future experiments on evaporating drops are clearly warranted.The weight of evidence suggests that glyoxal oligomers formed through acid catalyzed pathways are reversible based on thermodynamic considerations, but kinetically the glyoxal oligomers appear to be stable, at least for timescales of hours [125]. The impact of hydrophobic surfactants on the evaporating cloud droplet surface is an additional uncertainty related to surface accommodation/evaporation of water soluble gases.Biogenic fatty-acid coatings of the kind reported by Tervahattu et al. [201] appear to be more effective at inhibiting mass transport across the air-water interface than shorter, more soluble organic surfactants.Both short-chain and longer surfactants can affect the rates of surfacial chemical reactions; for the reactions studied to date, this effect is primarily through enhancing solvation of reagents at the surface.For a comprehensive review of atmospheric aerosol organic films the reader is referred to Donaldson and Vaida [202]. Our mechanistic understanding of the gas-phase formation of multi-functional organic compounds is limited by the availability of fundamental laboratory kinetic and product studies:  Of particular interest is aromatic ring fragmentation as this reaction channel is known to produce dicarbonyls in high gas-phase yields.Dicarbonyls can hydrate efficiently in solution to form diols and, thus, have high effective Henry's law coefficients.
 Our ability to predict Henry's law coefficients and OH-oxidation rate coefficients from structure-activity relations is limited for some functional groups such as organic nitrates [150,156,158].

Wet Deposition
Wet deposition refers to the processes of tracer scavenging and transport by precipitation.Atmospheric tracers (gaseous and particulate) can enter hydrometeors in various ways.Particles can be incorporated into hydrometeors through either nucleation scavenging (i.e., serving as CCN or ice nuclei, see Section 2) or impaction scavenging (i.e., being collected, through collision, by hydrometeors including cloud droplets, raindrops, and snow crystals).The incorporation of gaseous pollutants in hydrometeors is mainly through diffusion and subsequent aqueous-phase chemistry processes (see Section 3).The wet deposition (or removal) is contributed by both tracer transfer from cloud droplets to raindrops due to autoconversion (precipitation production), sometimes referred to as "rain-out", and the direct impact scavenging of aerosol particles and soluble gases by falling hydrometeors, sometimes referred to as "below-cloud" scavenging or "wash-out".Tracers captured in precipitation are removed from the atmosphere once the precipitation reaches the ground.However, part or all of the precipitation may evaporate before reaching the ground, in which case the tracers may be released back to the atmosphere from the evaporating hydrometeors.This can be an effective mechanism for tracer redistribution in the vertical.
In this section, we will mainly focus on the below-cloud scavenging aspect.Scavenging schemes for particles and gases by liquid and solid precipitation in some of the current North American and European regional air-quality models (same as listed in Table 1) are summarized in Table 5 and will be discussed below.The treatment of tracer release due to precipitation evaporation and the relative importance of in-and below-cloud scavenging to total wet deposition are also briefly discussed.
Table 5. Summary of below-cloud gas and aerosol scavenging schemes in some of the North American and European regional air quality models (see Table 1 for developers and references for these models).

Model
Scavenging schemes aerosol by rain aerosol by snow and ice gas by rain gas by snow

CMAQ
Accumulation and coarse mode: completely absorbed by cloud and rain water.Aitken mode: treated as interstitial aerosol subject to scavenging by precipitation, [79].

N/A
Mass transfer for species participating cloud chemistry and equilibrium for other species [79].
Same as for rain but with a lower limit for E, and different formulas for D.
Same as for rain but with different formulas for D and V t .
Mass-transfer rates for most species; equilibrium for O 3 and CO 2 .[83].

N/A
N/A : scavenging coefficient (s −1 ); d: aerosol diameter; D: droplet diameter; D g : diffusivity for gaseous species; E: droplet-aerosol collection efficiency; K c : Mass transfer coefficient; N Sh : the Sherwood number; P: precipitation rate; V t : droplet terminal velocity; τ washout , washout time representing the amount of time required to remove all of the water from the cloud volume at the specified precipitation rate; Re and Sc; respectively, the Reynolds and Schmidt numbers of drops.

Below-Cloud Particle Scavenging
Below-cloud particle scavenging processes depend on the net action of various forces influencing the relative motions of aerosol particles and hydrometeors.Ideally, both aerosol particles and hydrometers need to be size-resolved in air quality models in order to quantify the below-cloud particle removal.Such an approach is used in GATOR in which the detailed coagulation processes between size-resolved aerosol particles and hydrometeors are treated explicitly [69].Both liquid phase and solid phase hydrometeors are considered in GATOR.Despite the explicit treatment of the coagulation process, the uncertainties can be significant due to the parameterizations used for the various collection efficiency mechanisms (e.g., Brownian diffusion enhancement, turbulent inertial motion, turbulent shear, gravitational settling, thermophoresis, diffusiophoresis, and electric charge) and for droplet gravitational settling velocities [211].
Due to the computational constraint, the majority of the air quality models reviewed here choose to use bulk precipitation information, although the aerosol particles are size-resolved to some degree in these models, either with a sectional approach (e.g., AURAMS, CAMx, STEM, WRF-CHEM/MOSAIC, CHIMERE) or a modal approach (CMAQ, WRF-CHEM/MADE/SORGAM, EMEP unified model, COSMO-MUSCAT, LOTOS-EUROS).In many of these models the below-cloud particle removal by rain is assumed to follow a first-order decay and quantified by a parameter called scavenging coefficient, Λ(d), parameterized as a function of precipitation intensity (P), droplet-aerosol collection efficiency (E), mean droplet size (D) or mean droplet setting velocity (V t ) (see the second column in Table 5).A recent detailed review on size-resolved Λ(d) by rain [211] suggested that the use of a mean (or representative) droplet size in place of a full droplet spectra may be acceptable because the sensitivity of Λ(d) to droplet size spectrum is smaller compared to the sensitivity of Λ(d) to the collection efficiency parameterization.The study also recommended to use the parameterization for E which includes the considerations for such collection mechanisms as thermophoresis, diffusiophoresis, and electrostatic forces, in additional to the primary collection mechanisms (e.g., Brownian diffusion, interception, and inertial impaction).This is in consideration that even with the enhanced collection efficiency, the theoretical Λ(d) values were still much lower than those derived from field observations.However, a subsequent numerical sensitivity study following [211] suggests that the large apparent underestimation of Λ(d) by the theoretical formulation, compared to the field derived values, can be mostly explained by the fact that the field derived scavenging coefficients are inherently affected by the vertical turbulent diffusion, which is not part of the scavenging process (hence not accounted for by the theoretical formulation) but does impact the field observations [212].The study therefore concludes that the current theoretical framework for Λ(d) can provide a reasonable approximation of below-cloud aerosol particle scavenging by rain in size-resolved aerosol transport models as long as the vertical diffusion process is represented in the models [212].
Several of the air-quality models reviewed here choose to ignore snow scavenging of aerosol particles, while a few others (e.g., AURAMS, CAMx) use a similar formulation for Λ(d) as for particle scavenging by rain but with a different parameterization for the collector (hydrometeors) size (D) and collection efficiency (E) for scavenging of aerosol particles by snow and ice.At present, there has not been a systematic study on the magnitude of differences between different snow-scavenging formulations.It is expected that the uncertainty in the parameterization of particle scavenging by snow may be greater than that by rain due to the difficulty in quantifying the non-spherical shape, orientation, and size spectrum of snow and ice particles [213][214][215], although this needs to be confirmed through a detailed assessment of the uncertainties of existing snow scavenging parameterizations and associated parameters.Furthermore, the relative importance of particle scavenging by rain vs. by snow has not been assessed, though intuitively the scavenging by snow may be more efficient.
Besides the two approaches discussed above (i.e., the explicit coagulation approach and the scavenging coefficient approach), two models (CMAQ, STEM) treat the precipitation scavenging process through a diagnostic cloud column model that treats cloud dynamics (entrainment/detrainment), aqueous-phase chemistry, and wet deposition [94,208] for certain species (sulphate, nitrate, ammonium).
To get a sense of how sensitive the modeled aerosol mass concentration is to the use of different Λ(d) formulas, the AURAMS base-case simulation at 15-km resolution as described in Section 2 was repeated for two days (August 9 and 10, 2004) using two different Λ(d) schemes, Andronache et al. [216] and Mircea et al. [217].As shown in [211] these two scavenging schemes encompass the extent of the variation in the current Λ(d) parameterizations, with the Andronache scheme giving highest Λ(d) and the Mircea scheme the lowest; hence, this pair of runs will provide a ballpark for the model uncertainty (in terms of modeled PM mass) due to the variability in the current Λ(d) formulations.Note that the scavenging process is more efficient for large particles dominating PM 10 mass than for fine particles determining the PM 2.5 mass.Thus, the differences in predicted mass concentrations from using the different scavenging coefficient parameterizations are expected to be larger for PM 10 than for PM 2.5 on short (hourly to daily) time scales.This assumption is confirmed by sensitivity test results from AURAMS.As seen from Figure 5b,c, the relative difference in the modeled daily mean mass concentrations between the two schemes are up to 10% for PM 2.5 and up to 20% for PM 10 in areas with moderate to heavy precipitation.It is noticed that the patterns in Figure 5b,c are slightly different from the precipitation distribution pattern in Figure 5a, e.g., the largest differences in Figure 5b,c are not necessarily at areas with largest amount of precipitation.This is because, for areas with heavy precipitation, if the majority of the PM mass has been removed rapidly and new particles have not replenished in time, the differences in mass concentrations from using different schemes would become smaller at later times.This is consistent with results from [211] where it was shown that the difference in bulk mass from using different Λ(d) schemes increased with time during the first several hours and then decrease with time at later times under heavy rain conditions.

Below-Cloud Gas Scavenging
The simplest approach for treating the below-cloud gas scavenging by rain in air quality models is perhaps the 'scavenging coefficient' approach (Λ(g), with g denoting different gaseous species), where the depletion of atmospheric gaseous concentration is assumed to follow a first-order decay (irreversible scavenging), similar to the Λ approach for particle scavenging discussed above.In AURAMS, this approach is used only for very soluble species (e.g., HNO 3 , NH 3 , H 2 O 2 , ROOH) while in CAMx, this approach is used for all gaseous species.
The Λ(g) approach does not require information on aqueous phase concentrations.Instead, Λ(g) is parameterized as a function of the gaseous diffusivity, precipitation rate, mean droplet size and settling velocity.The scavenging process involves the diffusion of gases to droplets and the subsequent aqueous-phase reactions.Thus, it is expected that Λ(g) depends more on total droplet surface area than on the precipitation rate, which is confirmed in a numerical study by Zhang et al. [218] for moderately to highly soluble species (SO 2 , HNO 3 , NH 3 , H 2 O 2 ) under weak precipitation conditions.Unless the parameterized mean droplet size can properly represent the total droplet surface area under various precipitation conditions, it is likely that parameterizing Λ(g) using a single droplet size without the consideration of a full droplet spectra will introduce large errors in modeling the scavenging for soluble species.The same study also showed that Λ(g) for several soluble species (such as SO 2 and H 2 O 2 ) can vary by one order of magnitude depending on conditions such as droplet pH value and concentrations of other chemical species, due to the aqueous-phase reactions.Therefore parameterizing Λ(g) only as a function of physical parameters of precipitation (e.g., precipitation rate, mean droplet size) may cause large uncertainties.Given the same precipitation amount, uncertainties in the Λ(g) approach are expected to be larger in weak precipitation, which has a large amount of smaller droplets and hence large droplet surface area, than in strong precipitation.
The second relatively simple approach is the Henry's law equilibrium approach which assumes that the gas-and aqueous-phase concentrations are in equilibrium always (i.e., reversible scavenging).This approach is used in several models reviewed here for some or all gaseous species (see the third column in Table 5 for details).Theoretically, only weakly soluble species can reach equilibrium with small size droplets [208,219].For example, equilibrium assumption should work well for species like O 3 and CO 2 , but not for very soluble species.Deviation from Henry's law equilibrium can be significant for moderately to very soluble species (e.g., SO 2 , H 2 O 2 , HNO 3 ), even for very small droplets [218,220].Errors in this approach are largest for very soluble species and for large droplets (or heavy precipitation).
The theoretically correct and more complex approach is the kinetic mass transfer approach (also reversible scavenging) for those species whose gas-liquid phase equilibrium cannot be established readily.This approach is of course much more computationally intensive, having to integrate kinetic equations.Uncertainties in this approach can arise from the representation of droplet spectra and the modeled aqueous-phase concentrations.This approach is used in CMAQ, STEM, WRF-CHEM for some soluble species.
The impact on model prediction of gaseous pollutant concentrations due to the different representations of below-cloud scavenging of gases has not been assessed.In theory, errors from using the kinetic approach should be smaller than those from using other approaches; however, this does not guarantee better model results due to its complexity and the uncertainty in the parameters for describing the kinetic uptake.Considering that the below-cloud scavenging of gases is much more effective compared to the below-cloud scavenging of aerosol particles (except for very large ones), further detailed studies are needed to quantify the uncertainties from this process and to identify the best approaches for use in air quality models.
Gas scavenging by snow is ignored in most air-quality models.A scavenging coefficient approach is used in CAMx, and only a few gases (e.g., HNO 3 , NH 3 ) are considered in AURAMS and STEM.In AURAMS, the scavenging of these species by snow is treated similarly to the scavenging of H 2 SO 4 vapor which is based on an analogy to water vapor deposition rate on snow and ice [74].Apparently, knowledge of gas scavenging by snow is very limited and more studies are needed both in field and laboratories and through targeted experiments.

"Rain-Out" vs. "Wash-Out" and Precipitation Evaporation
As mentioned in the beginning of this section, wet deposition and removal of atmospheric tracers consist of both the below-cloud scavenging and removal, or "wash-out", and the removal of tracers in cloud water (through any of the processes: nucleation scavenging, impact scavenging of interstitial aerosol particles, gas diffusion and condensation on droplets, and in-cloud production) due to precipitation production (auto-conversion), or "rain-out".Most of the models reviewed here do consider both wet deposition/removal processes with the exception of LOTIS-EUROS, which does not currently consider the rain-out component [88].GATOR has a full, size resolved representation for the hydrometeor-hydrometeor and hydrometeor-aerosol microphysics, and hence the autoconversion (or precipitation production) and the tracer transfer (from cloud droplets to rain droplets) are represented explicitly.The approaches taken in AURAMS and WRF-CHEM are similar, i.e., using the bulk autoconversion rate (or precipitation production rate) as the tracer removal rate [73,74].Both CMAQ and STEM treat rain-out removal within their respective diagnostic 1-D cloud models [94,208].In CAMx, the rain-out process is treated as accretion of cloud droplets (or collection of cloud droplets by precipitation) with a scavenging coefficient parameterized by mean droplet diameter, droplet fall speed, collection efficient, and hydrometeor number density [91].CHIMERE, COSMO-MUSCAT, and the EMEP unified model share the same algorithm for in-cloud scavenging and removal, i.e., using a first order scavenging rate parameterized by precipitation rate and prescribed scavenging efficiency (see Table 1).The relative importance of the rain-out mechanism vs. the wash-out mechanism differs for particles and gases.For gases, the wash-out process is expected to be more important than rain-out because of higher tracer concentrations at lower levels compared to at cloud levels.However, the removal of atmospheric aerosol particles through rain-out is more important (as far as mass is concerned) than wash-out due to the low efficiency of impact scavenging (except for very large particles; see for example the Λ(d)-vs.-sizeplot in [211]: the scavenging coefficient is several orders of magnitude smaller for particles in 0.1-1 µm range than for particles in the 10 µm and up range).In contrast, with an average activation diameter around 0.1 µm, nucleation scavenging is much more efficient for particles in accumulation range.In addition, tracer mass removed through rain-out would also include the fraction produced via aqueous-phase oxidation; for example, a past study has indicated that a considerable portion (about 50%) of wet deposited sulphate can be attributed to in-cloud production [74].
Another process that will have an impact on wet deposition and removal of atmospheric tracers is precipitation evaporation.It can play an important role in redistributing tracer mass in the vertical.As rain falls through unsaturated air, some or all may evaporate before reaching the ground, and consequently tracers present in the rain (due to in-cloud or below-cloud scavenging or in-cloud production) may be released back into the atmosphere closer to the ground.Whether and how the tracer release during precipitation evaporation is treated can have a significant impact on model predictions of ambient tracer concentrations and wet deposition.The raindrops may evaporate completely or only partially before hitting the ground.Depending on the phase (i.e., dissolved in solution or as solid core) and volatility of the tracers in the raindrops, they may be released to the gas phase (through equilibrium mass transfer) as the raindrops evaporate or remain in the drops until complete evaporation has occurred.Amongst the models reviewed here, only two do treat the precipitation evaporation process.GATOR has explicit treatment for the hydrometeor evaporation/sublimation during settling by accounting for vapor exchange and change in hydrometeor surface temperature [69].Aerosol cores are retained in droplets until complete evaporation.AURAMS simply assumes that the rate for the release of tracers to the atmosphere equals to that of the bulk precipitation evaporation.As discussed in [74] this simple approach is likely to result in overestimation of the release of particles since they may have a tendency to remain in the rain droplets until complete evaporation.A sensitivity test by turning off the release of particles from precipitation evaporation showed a significant impact on the modeled ambient concentration and wet deposition flux of sulphate [74].Although the majority of models reviewed here do not explicitly consider the evaporation process, the effects of evaporation on gases may be implicitly represented in some of the models using equilibrium partition between gas and aqueous phase (rain water).
Similar to the case of modeling aqueous-phase chemistry, large uncertainty in modeling precipitation scavenging can arise from the modeled precipitation and other relevant microphysical parameters, such as autoconversion and precipitation evaporation.Our confidence in these modeled microphysical parameters is relatively low.

Conclusions
In this paper we reviewed the representations of cloud processing of gases and aerosols in some of the current regional air quality models in North America and Europe.We have focused on a number of key processes: aerosol activation (or droplet nucleation), aqueous-phase chemistry, and wet deposition (removal) by precipitation, and we have mostly focused on warm cloud processes.Through some simple sensitivity tests we attempted to illustrate some of the potential impact on modeled regional air quality from the variability in current representations for several processes, namely aerosol activation, potential droplet aqueous-phase pathways for atmospheric secondary organic aerosol formation, and the parameterization of below-cloud scavenging of aerosol particles.The following are some of the findings.

Aerosol Activation
Although aerosol activation or nucleation scavenging of aerosols is an important process in aerosol-cloud interaction and has been an active area of research within the climate modeling community, it has not been a process attracting significant attention within the air quality modeling community.The approaches to represent this process in the current regional air quality models are varied.A few of the sectional models have explicit representations, from mechanistically based parameterizations to simple parameterizations of using an empirical formula or a fixed activation diameter.For models with modal or broad sectional (2-4 sections or modes) representation for the aerosol size spectrum, it is usually assumed that all the accumulation mode aerosol particles, or all particles in the accumulation mode range, are activated and incorporated in droplets in the presence of cloud.The sensitivity model runs conducted, for a two-month period over eastern North America, with two different activation parameterization schemes, showed a significant impact on the modeled droplet number concentration and on the averaged size distribution of the modeled ambient particulate matter.The impact on the modeled average (over the 2-month period) ambient PM 2.5 concentration is limited (~a few percents), but the impact on the modeled average PM 1.0 concentration is considerably higher (up to 10 percents).Also, with the increasing interest in the area of air quality-climate interaction and the potential model applications in this area, there is a need to re-examine the representation of aerosol activation in current air quality models.

Inorganic Aqueous-Phase Chemistry
Almost all of the regional air quality models reviewed here have some representation for the aqueous-phase oxidation pathways leading to the production of sulphate in cloud.The models differ in chemistry mechanisms, from more complete atmospheric aqueous-phase chemistry, to sulphur oxidation focused chemistry, to highly parameterized single first-order reaction representation.Almost all models use a bulk approach for the aqueous-phase chemistry; a few have an option to use a variable-size-resolution-model approach to allow either a bulk or, when necessary, a two-bin representation in order to separate the droplets formed on larger, more alkaline particles from those formed on smaller, more acidic particles.Models also differ in how cloud water pH is determined, which will have an impact on aqueous-phase production of sulphate under certain conditions.As pointed out, though not reviewed here, the largest uncertainty in modeled aqueous-phase sulphate production may arise from the uncertainty in modeled cloud and liquid water content.
To date, the emphasis with regard to aqueous-phase feedback to aerosols is on the aerosol bulk mass within the air quality community.A few of the current sectional models do consider size distributed mass addition over the aerosol size spectrum from the aqueous-phase production, whereas this process is difficult to resolve in other air quality models due to their particular (or lack of) aerosol size representations.As discussed, the combination of aerosol activation and aqueous-phase secondary aerosol production has profound impact on the cloud processed aerosol size distribution and hence on aerosol optical properties.Again, this area may deserve some attention within the air quality modeling community.

Cloud Processing of Water-Soluble Organics
The cloud processing of water soluble organics and its potential impact on atmospheric SOA formation is a new area of active research within the atmospheric chemistry community.Given the large mass of biogenic and anthropogenic VOCs emitted into the atmosphere compared to the SO 2 emissions, it is important to assess the role of aqueous-phase production of SOA.Based on the current literature, in most cases, the liquid phase reaction term (in E2) is the rate limiting step in the uptake of water soluble organic gases.Thus, water soluble gas uptake to cloud droplets depends critically on available liquid water content, species K H , species liquid reaction rate, and species liquid diffusion coefficient.Existing studies have assessed the importance of the effective K H in determining the potential for a given species to partition to the aqueous phase and have suggested a minimum value of 10 3 M/atm.A recent review of the structure activity relationships (needed for the calculation of effective K H ) has revealed that most bifunctional species have K H greater than the 10 3 M/atm threshold, especially if one of the functional groups is an aldehyde or a peroxide.The weight of evidence from recent laboratory studies suggests that during the daytime the radical reactions dominate cloud organic chemistry, largely OH-initiated oxidation converting aldehyde groups to carboxylic acid groups.Of course, the details of the daytime radical chemistry (propagation vs. termination reactions, role of reservoir cycling) remain to be determined for atmospherically-relevant cloud conditions.
The results from the AURAMS sensitivity run, designed to give an upper-limit estimation of the impact of potential cloud droplet processing of atmospheric water soluble organic gases on SOA formation, suggest that indeed water soluble organic gas uptake to clouds and subsequent processing can be an important mechanism in addition to the traditional secondary organic gas uptake to the particle organic phase.The simple test with the additional pathway seems to result in improved evaluation statistics for modeled PM 2.5 OA compared to the IMPROVE network data for a 2-month period, e.g., reduced negative model bias and RMSE and improved slope for the linear correlation between model and observation.It should be noted that aqueous production of SOA occurs under different conditions than traditional SOA partitioning, such as high liquid water contents (clouds, high relative humidity).This will be particularly true in locations with large isoprene emissions and high cloud water contents such as the southeastern US.Overall, aqueous chemistry can also help explain the high OM/OC ratios and dicarboxylic acid measurements in ambient particles.

Wet Deposition
The majority of the models reviewed here use the scavenging coefficient (Λ) approach for below-cloud aerosol scavenging by rain or snow.The AURAMS sensitivity test for a two day period showed that the modeled daily ambient concentrations under rain conditions can differ by up to 10% for PM 2.5 and by up to 20% for PM 10 , using two different theoretical Λ parameterizations (corresponding to the lower and upper bounds).This suggests that the uncertainties in below-cloud aerosol scavenging process are smaller than many other processes considered in the model.This is because particle mass is dominated by large particles and existing theoretical scavenging coefficient formulas agree with each other quite well for particles larger than 3 m.It is also due to the fact that particles in the size range of 0.1-2 μm have lowest sensitivity to the scavenging process.The current Λ formulas used in the air quality models for snow scavenging are similar to those for rain scavenging.There is a lack of detailed and systematic studies to assess the uncertainty in the parameterization for scavenging by snow.Given the same precipitation amount, snow scavenging is expected to remove more particles from the atmosphere than rain scavenging due to the larger surface area of snow particles; and it is also expected that the uncertainties in current snow scavenging parameterizations will be larger than those in rain scavenging.
Various methods have been used for gas scavenging by rain, including the Λ approach (irreversible scavenging), Henry's law equilibrium (reversible scavenging), and kinetic mass transfer approach (also reversible scavenging).In some cases different approaches are used in the same model for different species (depending on the physical and chemical properties of the species).The uncertainties in gas scavenging by precipitation are known to be large, and further assessment of the uncertainties by comparing the different approaches within the same model framework would be worth pursuing.Scavenging (or uptake) of gases by snow is also an area needing further study both in the context of deposition and as part of snow chemistry and processing.composition measurements.Two anonymous reviewers provided valuable comments and suggestions which are much appreciated.

Figure 2 .
Figure 2. Impact of aerosol activation on modeled ambient sulphate mass: top-left panel-base-case (with Jones scheme) surface sulphate 2.5 concentration averaged over the simulation period; top-right panel-relative difference (in percent) in the modeled average surface sulphate 2.5 concentrations using the different activation schemes; bottom panels are the same as the top ones except for sulphate 1.0 .

Figure 3 .
Figure 3. Modeled average PM sulphate size distribution at two IMPROVE sites: Dolly Sods, WV (left) and Presque Isle, ME (right); black lines-base case (Jones scheme), red lines-sensitivity (Abdul-Razzak & Ghan scheme), blue lines-relative difference between the sensitivity and the base-case runs.

2 . 3
for wet aerosol chemistry). Hydrated formaldehyde can react with OH to form formic acid with a rate coefficient of 1.1 × 10 −9 M −1 •s −1 .The non-hydrated form of HCHO can also react with HSO 3 − to form HOCH 2 SO 3 − + OH − .In turn, the HOCH 2 SO 3 − can react with OH to form formic acid, HO 2

Figure 5 :
Figure 5: Modeled daily mean precipitation (a) and percentage difference in PM 2.5 (b) and PM 10 (c) mass concentrations from using two different scavenging schemes, Androche (A) and Mircea (M), for 10 August 2004.The percentage differences are calculated as ((M−A)/A) × 100.

Table 1 .
Summary of model representations of some of the in-cloud processes in a number of current regional air quality models in North America and Europe.
* All sizes referred here are dry sizes.

Table 2 .
Effective Henry's law coefficients for the dissolution of several organic species.

Table 3 .
OH-initiated aqueous oxidation of representative organic compounds with different functional groups (298 K).