dered at 100% r.h. Because the droplet should evaporate in such conditions, both thermophoresis and diffusiophoresis should push away the aerosol particles. Therefore, in our opinion, no scavenging should be observed due to phoretic forces in the considered conditions.
Bae et al.,  considered the thermophoretic and diffusiophoretic contribution to the below-cloud scavenging of aerosol (range 10–3 μm - 102 μm), by assuming a lognormal raindrop size distribution. The collection efficiency due to phoresis was evaluated, varying the air relative humidity between 50% and 90%, the ambient temperature, and the difference of temperature between air and raindrop surface (1˚C, 3˚C, 5˚C). The results of the model shown in Table 3 and Figures 4 and 5 are of the paper evidence positive scavenging coefficient for both thermophoresis and diffusiophoresis. The relative contribution of both thermophoresis and diffusiophoresis turns out to be about 80%, when particle diameter is between 0.1 and 1 μm. Examining the data reported in Tables 2 and 3 of the paper, we disagree with the conclusions of the authors. In particular, as the driving force of the diffusiophoretic process is positive (), the droplets evaporate, and by including hydrodynamic Stefan flow, the diffusiophoresis pushes the particles away from the droplets . Therefore thermophoresis and diffusiophoresis act in a reverse way, not in the same way, as claimed by the authors, who write that “as r.h. becomes lower, diffusionphoresis is more effective because particles move to raindrop faster due to rapid evaporation of raindrops, and thus the difference of water vapour density becomes larger”.
In order to evaluate the relative importance of phoresis in the scavenging process, it is also useful to consider published papers on wet scrubbers, devices which can remove even particles smaller than a few micrometers from indoor or industrial gases, by bringing them into contact with liquid droplets.
Sparks and Pilat , Pilat and Prem [27,28] calculated particle collection efficiencies for a single droplet (d = 100 μm) and for spray droplet scrubbers, considering inertial impaction, diffusiophoresis and thermophoresis, showing that positive diffusiophoresis (condensation) could substantially improve the collection efficiency of sub-micron particles. It was also seen that the effect of phoretic forces becomes greater as the droplet size is increased, because of the longer time necessary to come into equilibrium with the gas. The model assumed a potential flow for the gas flow profile around the droplet moving at settling velocity and a stagnant gas layer surrounding the droplet across which exists a vapour and temperature gradient. However, according Wang et al., , at the considered Reynolds number, the viscous nature of flow of air past drops is still of great importance and the temperature gradients considered by the above authors are not typically found in the atmosphere.
Horst  considered the scavenging of aerosol from gas for a condensing steam/air system on a cold wall, assuming that the steam/air system was saturated everywhere. Following Waldmann and Schmidt theory, he found that the particles are carried in the direction of water vapour flux at a velocity: V = 103 mw (mw, mass flux of vapour, g·cm–2·s–1; V, cm·s–1), two orders of magnitude higher than thermophoretic velocity.
Viswanathan  performed a similar numerical study, for a range of collector radii between 50 and 4000 mm and particle size range between 0.05 and 10 mm, to include particles in the “Greenfield gap”, by considering Stokes flow. The results for a water temperature of 10˚C and a gas temperature ranging between 20˚C and 95˚C show an increase of particle collection efficiency of several orders of magnitude, which is most prominent at lower Reynolds number.
It is worth noting that both Pilat and Prem [27,28], and Viswanathan  considered the temperatures of the drop and of the air to be constant (Tair > Tdrop), with thermophoresis and diffusiophoresis consequently behaving in the same way. This assumption does not always reflect the atmospheric situation. The relative importance of phoretic phenomena is also highlighted both theoretically and experimentally from Vasudevan et al.  in the removal of particulate from fuel gas in spray scrubbers.
3. Laboratory and Fields Studies on Drop Scavenging, In-Cloud and Below-Clouds
3.1. In-Cloud Scavenging
It should be noted that the distinction between in-cloud and below-cloud scavenging is not absolute, as precipitation which forms at higher altitude, falling into a cloud layer, leads to scavenging even within the cloud. Here, the review focuses on only a few published papers, which consider phoretic forces in the in-cloud scavenging process.
The phenomena of thermoand diffusiophoresis were predicted by Stefan  and subsequently observed by Facy . Facy was the first to suggest that diffusiophoresis may contribute to the scavenging of the atmospheric aerosol during cloud droplet growth. He photographed growing and evaporating water drops in an aerosol-laden atmosphere and showed that a dust-free space surrounded the evaporating drop, due to diffusiophoretic forces set up by the flux of water vapour away from the drop. Conversely, a growing drop encouraged the migration of aerosol to the surface of the drop.
Goldsmith et al.,  conducted a cloud-chamber experiment in which salt droplets were grown up to 10 μm in radius in the atmosphere of a radioactive aerosol (0.01 μm - 0.1 μm). They assessed the absolute efficiency of the diffusiophoretic mechanism in capturing the submicron particles by the growing droplets. Concerning the role of diffusiophoresis in the atmosphere, the authors concluded that the efficiency of diffusiophoresis in scavenging particles from the atmosphere should be low, for each condensation cycle inside a cloud, compared with any other processes. However, the overall efficiency is dependent on the number of condensation-evaporationcondensation cycles taking place in natural clouds. It should be pointed out that in the performed experiments thermophoresis and diffusiophoresis were present simultaneously, and the results showed a prevalence of diffusiophoresis with respect to thermophoresis by considering sub-micron particles.
Phoretic forces could play a role in the discrepancy observed in mixed clouds between the concentration of ice nuclei and ice particles. There are abundant observations of enhanced ice formation in regions where cloud droplets are evaporating in cumuliform, stratiform, and wave clouds [49,61-67]. Moreover, it has been speculated that contact nucleation of aerosol on the external surface of the supercooled droplets may be responsible for this “evaporating freezing”. On considering the contribution of phoretic effects, it should be concluded that thermophoresis is prevalent with respect to diffusiophoresis.
Cooper and Vali  studied simple layer clouds. The patterns in the development of ice in these clouds suggest that ice originates in association with the initial condensation process, near the upwind edge of the cloud. Since continuous ice production does not occur beyond that region, the ice development can be attributed to condensation-freezing and contact nucleation. In this case a prevalence of diffusiophoresis during water vapour condensation should be present.
These contradictory results was possibly overcome by Durant and Shaw , whose laboratory experiments showed that the freezing temperature of an evaporating supercooled drop will suddenly become higher, once an immersed ice nucleus contacts the internal drop surface (“evaporation freezing”). That is to say, when a particle is immersed and very near or touching the surface of a supercooled droplet, it is a better ice nucleus than when it is far from the surface. The difference in the two modes of nucleation was found to be between 4˚C and 5˚C for all the materials tested. Therefore the notion of contact nucleation should be generalized to include surface crystallization from particles contacting a supercooled drop from the inside-out, as well as from the outside-in. In conclusion, if internal contact nucleation were the prevalent freezing mechanism in the above reported cases of evaporating cloud droplets, the hypothesis that thermophoresis was prevalent over diffusiophoresis should decline.
Ladino et al.,  performed experiments in the CLINCH chamber, using pure water droplets of radii between 12.8 and 20.0 μm, allowed to fall freely, and to collide in laminar flow with particles having radii between 0.05 and 0.33 μm. Evaporation of the droplets was taken into account since the relative humidity inside the chamber was below 100%. Even if the data are very scattered, for droplet size at 12.8 μm and particle radii in the range from 0.05 to 0.33 μm, a decrease of the collection efficiency can be observed, higher than the Wang et al.,  theoretical model, with a minimum at about 0.25 μm. The trend is partially in agreement with previous experimental studies by Lai et al., , performed with much higher drop diameter. In addition, in experiments with a fixed aerosol radius (0.24 μm), Ladino et al.,  found a decrease of collection efficiency with increasing droplet size (12.8 - 20.0 μm) as predicted by theory , although values were smaller by up to one order of magnitude (r.h. = 90%). This could mean that for evaporating droplets diffusiophoresis decreases the scavenging.
3.2. Below-Cloud Scavenging
Considering laboratory experiments, the different experimental conditions (temperature, relative humidity, droplet and aerosol diameter, etc.), make it difficult to compare results. Several laboratory measurements of capture efficiency have been made, disregarding the relative humidity present in the experimental apparatus, making it impossible to evaluate the contribution of thermodiffusiophoresis [71-73].
In some cases the humidity in the experimental apparatus has been measured. In a laboratory experiments, Pranesha and Kamra  considered particle size (1.9 μm - 6.4 μm) and 3.6 mm - 4.8 mm diameter drops. The collection efficiency increased as the particle size increased for a fixed drop size, and decreased sharply with increasing drop size for each particle size. Although the humidity in the aerosol chamber was observed to be between 35% to 50%, the authors claim that evaporation from the surface of the drops was expected to be negligible, and consequently phoretic effects were not considered.
Ebert et al.,  performed a field experiment in which aerosol (radius varying from 0.19 to 1.8 μm) was used to measure particle scavenging processes in the atmosphere during real precipitation events. Even if the r.h. varied between 65% - 100% during the experiments, no evaluation was made concerning the effect of phoresis.
The first laboratory experiment considering all the scavenging mechanisms (Brownian diffusion, hydrodynamic, phoretic and electrical effects) was performed by Wang and Pruppacher . Water drops (150 to 2500 mm equivalent radius) were allowed to fall through an aerosol chamber in sub-saturated air. Although their data show good agreement with the efficiency calculations of Grover et al.,  as far as the dependence on the drop size was concerned, their measurements were for only one particle size (0.25 μm radius; 23% r.h.). The collection efficiency showed a minimum for 200 μm drop radius and decreased for drop with radius higher than 500 mm.
The authors write: “despite the agreement found between theory and experiment, it appears necessary to test the theory of Grover for other temperature and humidity conditions and for additional drop and particle sizes”.
Leong et al.,  investigated the effect of particle size (0.58 μm - 3.2 μm radius) on scavenging for evaporating cloud drops at 30% r.h., with mean radii of 55 to 93 μm. Good agreement between experiments and theory was found for larger particles scavenged by inertial impaction, but not for the smaller ones also scavenged by phoretic forces .
In an effort to evaluate the contribution of phoretic forces in the scavenging process, Wang et al.,  performed experiments assuming a quasi steady-process and comparing the collection efficiency obtained with a water droplet evaporating (500 μm radius) and an oil droplet with very low vapour pressure, in the same experimental conditions. The data obtained indicated that for 500 μm diameter droplets an efficiency of 1.6 × 10–3 can be attributed to phoretic forces.
Concerning field measurements, the scavenging coefficient is usually calculated from measurements of change in aerosol size distribution with different rainfall in different environments. For very small and very large particles, there is prevalently an agreement with theoretical studies [19,25,38]. However, theoretical parameterizations prevalently underestimate observed L values by one to two orders of magnitude for particles in the 0.1 μm - 1 μm radius range and by one order of magnitude for particles smaller than 0.1 mm, compared to the available field measurements .
To avoid experimental biases, it should be guaranteed that in the course of precipitation the air at the measuring station is not replaced by another mass with a different aerosol loading, i.e. no significant convective activity should be present and no frontal passage should occur over the time of interest. This requirement is not always met in field experiments. In addition, a change of r.h. during the event may alter the aerosol size distribution due to hygroscopic effects. Since L is very sensitive function of the aerosol diameter, a change of r.h. will result in a change of the depletion rate. Precipitation is hardly even a constant process. Large fluctuations in the rainfall rate and the drop size distribution are often observed. These parameters influence L directly. If the above requirements are fulfilled, the differences between theoretical and experimental data could depend on theoretical reasons, i.e. on the fact that phoretic forces acting upon submicron particles are underestimated by theory, on vertical turbulent diffusion, or on the dependence of the scavenging coefficient on aerosol chemical composition [80-82].
Davenport and Peters  analysed wash-out coefficients for moderate rains (1.1 mm·h–1 to 3.0 mm·h–1) over sampling times of 5 to 6.5 hours. Based on data collected during the events, the authors claimed that theoretical scavenging coefficients agree quite well with experimental results in the coarse aerosol range. For particle in the range 0.007 μm - 0.022 μm and 0.4 μm - 1 μm, the model which consider inertial impaction, interception, Brownian diffusion, and phoresis, predicted one or two orders of magnitude less than field measurements. No available experimental data concerned aerosol in the 0.02 μm - 0.4 μm range. The results of Davenport and Peters are rather surprising, as they did not find any dependence on aerosol radius at all in their scavenging coefficients in the measured ranges. One may speculate that scavenging effects were not properly distinguished from transport phenomena.
Moreover, Radke et al.,  measured the scavenging of aerosol particle (natural “man made” as well as plume and volcanic aerosol, d > 0.01 μm) mainly by showertype precipitation. Despite large variations in the nature of the aerosol particles and precipitations, the scavenging collection efficiencies as a function of particle size showed marked similarity. For particles ³ 1 μm in diameter, where inertial impaction dominates scavenging, measurements were in good agreement with the theoretical results of Dana and Hales , Wang et al.,  and the laboratory measurements of Lai et al., . For submicron particles the scavenging collection efficiencies were generally much higher than values theoretically predicted for scavenging due to Brownian and phoretic effects.
Schumann , Volken and Schumann  performed a field study on Mt-Rigi (central Switzerland), measuring the below-cloud scavenging of aerosol particle during 7-year winter precipitations. The scavenging coefficient showed a dependence on particle diameter.
For particle smaller than about 0.6 μm, the value of L was 10–5 s–1 on average, increasing to values around 3 × 10–5 s–1 for particles 1 μm < d < 3 μm. A small local maximum near 1 μm was also visible. The unexpected result was the absence of an evident minimum of the scavenging coefficient. Based on data collected during the events (raindrop size distribution, relative humidity), and collection kernels from Flossmann , the theoretical scavenging coefficients calculated agreed quite well with experimental results in the coarse aerosol range, while for submicron particles they tended to be 20 - 40 times smaller. An attempt was made to explain the differences, by considering the influence of r.h., the particle shape, electric effects, and phoretic forces.
Laakso et al.,  studied the dependence of the scavenging coefficient for particle diameters of 10 - 500 nm and for rain intensity 0 - 20 mm·h–1 at a boreal forest site in Southern Finland. When the rain intensity increased from 0.5 to 9 mm·h–1, the corresponding scavenging coefficient changed from 1 × 10–5 to 4 × 10–5 s–1. The scavenging coefficient showed a very little size dependence, with a minimum in the 0.1 μm - 1 μm range. In the accumulation size range, the results were higher than those suggested by theories presented by Slinn , Schumann  and McGann and Jennings . To explain the differences, the authors refer to reasons suggested by Volken and Schumann .
Chate and Pranesha  noted that the collection of atmospheric particles during thunderstorm rain by the mechanisms of Brownian diffusion, directional interception, inertial impaction and electrostatic charge effects cannot adequately explain the depletion in particle number concentrations observed during these field experiments. They suggested that the phenomenon of thermophoresis, diffusiophoresis, coagulation and condensational growth of hygroscopic particles may be important during thunderstorm rain as removal processes for atmospheric particulate.
In conclusion, it seems that even parameterizations which take into account the additional collection processes due to thermophoresis, diffusiophoresis, and electrostatic forces, are not able to explain the discrepancies between field and observed scavenging coefficient in the “Greenfield zone”.
The review of the published papers on the role of phoretic forces in atmospheric scavenging of aerosol can be summarized as follows:
• Theoretical models agree that Brownian diffusion appears to dominate scavenging of aerosol in the size r < 0.1 μm range, while inertial impaction dominates scavenging in the size r > 1 μm range. Thus, there is a minimum collection efficiency for particles in the approximate range of 0.1 μm - 1 μm, where phoretic and electric effects are felt. The exact depth, width and position of this minimum depend on the properties of the aerosols and hydrometeors, and on ambient conditions.
• By considering in-cloud scavenging, where scavenging theories assume a thermal equilibrium between the drop and the environment, evaporation being considered the only factor determining the temperature gradient, thermophoretic and diffusiophoretic forces (by including Stefan flow) act in the opposed direction. Considering as an example an evaporating droplet, the thermophoretic force should push the aerosol towards the droplet, and diffusiophoresis should push the aerosol away from the droplet. Disagreement exists in both laboratory experiments and the theoretical model on the prevalent force. Some authors evaluate thermophoresis as prevalent for particles of radius less than 1 μm, while diffusiophoresis is prevalent for particles with higher radius.
• Considering below-cloud scavenging, the proposed theoretical models follow different approaches. Chate , Andronache et al., , and Bae et al.,  assume a constant temperature for droplets and air, neglecting the variation of droplet temperature due to growth or evaporation. In our opinion, the above authors in the theoretical approach consider only thermophoresis and pure diffusiophoresis, by neglecting Stefan flow. Considering an evaporating droplet, the movement of the aerosol due to diffusiophoresis process was considered towards the droplet, while by including the Stefan flow the aerosol should be rejected.
• Studies on the scavenging process in spray water scrubbers, based on the mechanisms valid in the atmosphere, highlight a prevalence of diffusiophoresis with respect to thermophoresis.
• Concerning field measurements, for very small and very large particles, there is generally an agreement with theoretical studies. However, significant discrepancies between observations and theoretical estimates have been found for sub-micron particles in the “Greenfield gap”, with measured values being one to two orders of magnitude higher than predicted. This discrepancies remain even the contribution of phoretic forces is considered.
• In conclusion, more laboratory experiments under controlled conditions and the collection of field data under different rain conditions and other environmental conditions are required, to evaluate the contribution of phoretic forces to the scavenging processes.