5.5. Air pollution and climate#

Cloud formation, aerosol-cloud interactions

Aims:

  1. To describe the direct effect of aerosols on climate (absorption vs. scattering of radiation).

  2. To describe the role of aerosols on cloud formation

  3. To describe the impacts of aerosols on cloud characteristics and lifetime.

Key points:

  • The negative forcing from aerosols has offset some warming from greenhouse gases in the 20th century, but this effect is decreasing due to policies to tackle air quality.

  • Aerosols can scatter and absorb radiation directly. Most aerosol species scatter solar radiation (net negative forcing).

  • Aerosol particles are needed to provide a nuclei for cloud droplets to form on

  • Droplet formation is best described using the Koehler equation, which combines the Kelvin equation (curvature effect) with the Raoult equation (solute effect)

  • An increase in CCN makes clouds brighter, causing a negative forcing.

Direct effects on climate#

Aerosol-radiation interactions (scattering vs. absorption)#

Aerosol-radiation interactions (also known as aerosol direct effects) refers to the scattering of solar radiation and absorption of solar and terrestrial radiation due to the presence of aerosol particles in the atmosphere. In general, radiative forcing due to aerosol-radiation interactions is proportional to the amount of aerosol particles in the atmosphere and their ability to extinguish radiation, which in turn is a function of particle size and wavelength.

The schematic on the left of Fig. 5.19 shows that the scattering of shortwave radiation reduces the amount of solar energy that reaches the surface, leading to a cooling of the atmosphere and surface, expressed as a net negative forcing at the TOA. The schematic on the right of Fig. 5.19 shows that an absorbing aerosol (e.g., black carbon) leads to a warming of the atmosphere because less solar energy reaches the surface and the TOA, resulting in a net positive forcing at the TOA. Note that the warmed atmosphere re-emits some of the absorbed radiation at longer wavelengths, some of which reaches the surface and partially offsets the decrease in surface temperatures.

../../_images/figure23.4.png

Fig. 5.19 Radiative influence of scattering (left panel) and absorbing aerosol particles (right panel).#

The scattering and absorption properties of different aerosol types depend on particle size (i.e., larger particles absorb more and scatter less), shape, chemical composition and mixing state. For most aerosol types, the scattering of solar radiation dominates the absorption of solar radiation, resulting in a net negative solar forcing. The major exception is Black Carbon (soot) which strongly absorbs solar radiation resulting in a positive forcing. Larger-sized aerosol particles (e.g., mineral dust) tend to strongly absorb terrestrial (longwave) radiation, which can result in a positive forcing.

Table 5.5 shows the net global mean forcing from aerosol-radiation interactions for different aerosol types and the forcing from all aerosol types combined, which is negative as most aerosol particles present in the atmosphere tend to be efficient scatters.

Table 5.5 Net global mean forcing from aerosol-radiation interactions for different aerosol types and the forcing from all aerosol types combined. Table from Boucher (2012).#

Aerosol type

Net radiative forcing (\(\mathbf{W \ m^{-2}}\))

Uncertainty range (\(\mathbf{W \ m^{-2}}\))

Fossil fuel sulphates

-0.40

-0.60 to -0.20

Fossil fuel + biofuel black carbon

+0.40

+0.05 to +0.60

Fossil fuel + biofuel organic aerosols

-0.12

-0.40 to 0.10

Biomass burning

+0.00

-0.20 to +0.20

Desert dust

-0.10

-0.30 to +0.10

All aerosol

-0.35

-0.85 to +0.15

Consider the case of scattering by a spherical particle of radius r, for which the scattering, absorption and extinction efficiency can be prescribed on the basis of theory, as a function of a dimensionless particle size parameter χ:

(5.21)#\[\chi = \frac{2 \pi r}{\lambda}\]

In general, the scattering of radiation is largest when \(\chi \simeq 1\) (i.e. the Mie scattering regime), which in the visible part of the electromagnetic spectrum is the case for those particles with radii comparable to the wavelength of visible light (ca. \(400-750 \ \mathrm{nm}\), Fig. 5.20). Note that most anthropogenic aerosol particles are efficient at scattering solar radiation because they are in the so-called “accumulation-mode” size range (\(0.1-1 \ \mathrm{\mu m}\) particle diameter).

../../_images/figure5.png

Fig. 5.20 Type of single scattering regime as a function of radiation wavelength and particle radius, assuming particles are spherical. Figure from W. Brune.#

Role of aerosol particles in cloud formation#

Cloud formation#

In the atmosphere, clouds form as a result of a combination of thermodynamic and microphysical processes. For a cloud to form, it is necessary for a large volume of moist air to become saturated. The vapour pressure of water in the parcel is referred to as \(e\), and \(e_s\) is the saturation vapour pressure of water in equilibrium with a plane surface of water at the temperature of the parcel.

The Clausius-Clapeyron equation allows us to calculate the variation of \(e_s\) with temperature.

(5.22)#\[\dv{e_s}{T} = \frac{L e_s M_v}{R T^2}\]

where \(L\) is the latent heat of vapourisation and \(M_v\) the molecular weight of water. If we assume that \(L\) is constant, independent of \(T\), then integrating the Clausius-Clapeyron equation gives:

(5.23)#\[\ln(\frac{e_s}{e_0}) = \frac{LM_v}{R T_0} - \frac{L M_v}{R T}\]

or

\[e_s = 6.11 \exp(19.83 - \frac{5417}{T})\]

where \(e_s\) is in \(\times 10^2 \ \mathrm{Pa}\) (\(\mathrm{mbar}\)) and \(T\) in Kelvin. As \(e_s\) increases exponentially with \(T\), if an air parcel cools (without changing \(e\)) the air parcel will change from being subsaturated to supersaturated.

../../_images/figure1.png

Fig. 5.21 Schematic of the saturation vapour pressure of water.#

When discussing cloud droplet formation, the saturation ratio (\(S\)) and supersaturation (\(s\)) (both unitless) are useful variables, defined as follows:

  • Saturation ratio:

    (5.24)#\[S = \frac{e}{e_s}\]

    where:

    • \(S < 1\) for a subsaturated (evaporation > condensation) environment

    • \(S = 1\) for saturation (condensation = evaporation)

    • \(S > 1\) for a supersaturated (condensation > evaporation) environment.

  • Supersaturation

    (5.25)#\[s = S-1=\frac{e}{e_s} - 1 = \mathrm{RH} - 100%\]

    with:

    • \(s = 0\) at saturation,

    • \(s < 0\) for a subsaturated environment

    • \(s > 0\) for a supersaturated environment

A parcel of air which is saturated with respect to water defines the concept of relative humidity (\(\mathrm{RH}\)), since relative humidity is defined as the mass mixing ratio \(\omega\) of a parcel at a particular temperature and pressure, compared to its saturation mixing ratio value at the same temperature and pressure:

(5.26)#\[\mathrm{RH} = \frac{\omega}{\omega_s(\rho, T)} \approx \frac{e}{e_s}\]

The equation for the saturation ratio only applies for a flat surface of pure water. Once we consider actual particles, we get into situations where curvature or the effect of a solute needs to be considered. In these cases, supersaturation needs to be considered relative to the equilibrium value of \(e\) for a particular curved surface, so we refer to \(e^*\) as the equilibrium vapour pressure, which can differ from \(e_s\) over a flat plane surface.

In the atmosphere, supersaturation can be achieved, in principle, via three basic mechanisms: radiative cooling , uplift, and mixing. Note that radiative cooling and mixing happen at constant pressure (isobaric), whereas uplift happens at constant energy (adiabatic). In the atmosphere, uplift or adiabatic expansion of ascending air is the most prevalent mechanism for a parcel of moist air to be cooled below \(T_d\) in order to reach supersaturation (see Fig. 5.22). We will see later that in an ascending parcel, the updraft velocity affects the rate of cooling and hence the amount of condensation which forms.

../../_images/figure23.2.png

Fig. 5.22 Schematic of changes in relative humidity (RH) and supersaturation for a parcel of air rising in a cloud.#

Koehler theory#

Cloud formation is best described by Koehler theory [Köhler, 1936], which combines two competing effects that most strongly determine the amount of supersaturation that a specific particle must experience in order to accumulate enough water to grow into a cloud drop. The first is a physical effect of curvature which increases the droplet’s equilibrium vapour pressure; this is called the Kelvin effect. The second is a chemical effect of the aerosol dissolving in the growing water drop and thereby reducing its equilibrium vapour pressure; this is called the solute effect or Raoult’s Law.

Kelvin effect (curvature effect)#

The Kelvin effect states that the equilibrium vapour pressure over a curved surface (e.g. a solution droplet with a finite radius \(r\)) is much greater than the saturation vapour pressure over a plane or flat surface.

The Kelvin equation is given as:

(5.27)#\[\frac{e_s(r)}{e_s(\infty)} = S = \exp(\frac{a}{r}) \qq{where} a=\frac{2 \sigma}{\rho_l R T}\]

Where \(\sigma\) is the surface tension between water and air, and is the density of the droplet evaluated at \(T\). \(R\) is the molar gas constant. At \(273.2 \ \mathrm{K}, \sigma = 0.0756 \ \mathrm{N m^{-1}}\).

Fig. 5.23 shows a plot of the Kelvin equation. Smaller particles need a higher supersaturation to grow because their surfaces are more curved and water evaporates from them more easily.

../../_images/figure3.jpg

Fig. 5.23 Vapour pressure over a curved surface for pure \(\mathrm{H_2O}\) at 298 K. Figure credit: W. Brune (Penn State University).#

Solute effect (Raoult’s Law)#

Aerosol particles are abundant in the atmosphere and some of the particles are hydrophilic and water-soluble. In terms of droplet growth, this solubility is fortunate because the saturation vapour pressure over a solution is lowered by the presence of the solute compared to that of pure water and hence growth is energetically favourable. This is described by Raoult’s Law.

Phenomenologically (and in a simplified manner) this vapour pressure reduction can be understood as follows: solute molecules (i.e. the chemical that is being dissolved in the solvent) occupy surface sites that would otherwise be occupied by water molecules (the solvent), which in turn prevents water molecules from evaporating from those sites as shown in Fig. 5.24. Consequently, a droplet containing, for instance, dissolved salt favours condensation more than would a pure water droplet of the same size. In other words, existing particles onto which water vapour condenses (cloud condensation nuclei; CCN) enable condensation at supersaturations well below those required for homogeneous nucleation (primarily because of their size) as discussed further below.

../../_images/figure23.4.png

Fig. 5.24 Schematic of the saturation vapour pressure for (a) pure water, and (b) an ideal solution. The white disks denote water molecules, and the black disks dissolved molecules. The reduction in vapour pressure from (a) to (b) corresponds to there being fewer sites from where water molecules can escape through the air–water interface. Figure from Lohmann et al. [2016].#

Raoult’s Law states that the saturation vapour pressure of an ideal solution is dependent on the vapour pressure of each chemical component and the mole fraction of the component present in the solution, so:

(5.28)#\[\frac{e^*(\infty)}{e_s(\infty)} = \frac{n_0}{n + n_0}\]

where \(e^*\) is the equilibrium vapour pressure over a solution consisting of \(n_0\) water molecules and \(n\) solute molecules.

For dilute solutions applied to droplets of size \(r\), Raoult’s Law can be approximated as:

(5.29)#\[\frac{e^*(r)}{e_s(r)} = 1 - \frac{b}{r^3} \qq{where} b = \frac{3i~m_s~M_w}{4\pi~\rho_l~M_s}\]

where \(m_s\) is the total mass of the solute, \(M_s\) and \(M_w\) are the molecular masses of the solute and water respectively, \(\rho_l\) is the density of the droplet, \(i\) and is the Van’t Hoff factor which accounts for ionic dissociation of the solute. An example is salt, \(\mathrm{NaCl}\), which splits into two ions in solution, \(\mathrm{Na^+}\) and \(\mathrm{Cl^-}\), which means that \(i = 2\) in this case.

Koehler equation: combining Kelvin equation and Raoult’s Law#

The Koehler equation combines the competing effects of the Kelvin equation and Raoult’s law as follows:

(5.30)#\[S = \frac{e^*(\infty)}{e_s(\infty)} = \exp(\frac{a}{r}) \left(1 + \frac{b}{r^3} \right)^{-1}\]

The Koehler equation describes the ratio between the equilibrium vapour pressure over a solution droplet, \(e^*(r)\) with radius \(r\) to the saturation vapour pressure over a plane surface of water \(e_s(\infty)\), with \(a\) and \(b\) as described in the Sections above.

As \(a \ll r\) and \(b \ll r^3\), expanding each term in \(S\) and neglecting higher order terms yields:

(5.31)#\[S = \frac{e^*(\infty)}{e_s(\infty)} \approx 1 + \frac{a}{r} - \frac{b}{r^3}\]

Fig. 5.25 shows a Koehler curve, i.e., the saturation ratio as a function of the droplet radius, and separate curves for the curvature and solute effects.

../../_images/figure5.png

Fig. 5.25 Exemplar Koehler curve (black line) with red dashed line indicating the contribution from Raoult’s Law and the dotted line that from the Kelvin equation. Figure from Ritchie [2017].#

From Fig. 5.25, it can also be seen that the Kelvin effect is negligible for droplets larger than \(1 \ \mathrm{\mu m}\) where all curves approach \(S = 1\).

The peak of a Koehler curve corresponds to the so-called critical saturation ratio, \(S_c\) at the critical radius , \(r_c\) (that is typically between \(0.1 \ \mathrm{\mu m}\) and \(1.0 \ \mathrm{\mu m}\)). \(S_c\) and \(r_c\) can be obtained by differentiating the Koehler equation with respect to \(r\) and setting the derivative to 0:

(5.32)#\[r_c = \left( \frac{3b}{a} \right)^\frac{1}{2} \qq{} S_c = 1 + \left( \frac{4a^3}{27b} \right)^\frac{1}{2}\]

\(S_c\) of a specific Koehler curve is the minimum saturation ratio that is required for the corresponding solution droplet to grow to cloud droplet size. Theoretically, it could grow even larger, but the growth rate becomes slower the larger the droplet, and more efficient growth mechanisms such as growth by collision-coalescence (a.k.a. accretion) with other droplets will take over. If the droplet has grown to \(r > r_c\), it is called an activated droplet. All droplets that have a critical saturation ratio \(S_c < S_{max}\) can thus be activated, where \(S_{max}\) is the maximum supersaturation reached in the cloud. \(S_{max}\) is itself a function of updraft velocity – in clouds with faster updrafts (such as storm clouds) \(S_{max}\) is higher, meaning smaller particles can be activated and more droplets are formed overall. In the context of cloud formation, the number of activated particles is always small (about \(200 \ \mathrm{cm^{-3}}\)) compared with the total number of particles in the atmosphere.

Activation of nucleation-mode aerosol particles requires very high supersaturations that do not exist in the atmosphere. In contrast, larger aerosol particles require smaller supersaturations to be activated, thus these larger particles act more readily as CCN, but it is important to remember that there are fewer of them.

Aerosol- cloud interactions#

Aerosol-cloud interactions (also known as aerosol indirect effects) refer to the modification of cloud properties due to changes in the aerosol population as summarized in Fig. 5.26.

../../_images/figure6.png

Fig. 5.26 Summary of aerosol-cloud interactions and their climate effects. Based on IPCC AR5.#

An increase in cloud condensation nuclei, CCN (i.e. a subset of the total aerosol population that can activate to form cloud droplets) number concentration will result in an increase in the number concentration of cloud droplets. These droplets compete for the available cloud liquid water content in order to grow to larger sizes. Under the assumption that the liquid water content remains constant, the increase in CCN number concentration will lead to more numerous but smaller-sized cloud droplets compared to those in an unperturbed “clean” cloud since droplet radius goes as \(r \propto \left(\frac{L}{N_d}\right)^\frac{1}{3}\) where \(L\) refers to the liquid water content (\(\mathrm{g \ m^{-3}}\)) and \(N_d\) is cloud droplet number concentration (\(\mathrm{cm^{-3}}\)).

For a cloud of thickness \(h\) (in \(\mathrm{m}\)), the cloud optical thickness, \(\tau\), is related to cloud properties as follows:

(5.33)#\[\tau \propto N_d^\frac{1}{3} L^\frac{2}{3} h\]

Therefore, for a constant liquid water content and a constant cloud thickness, \(\tau\) increases with \(N_d^\frac{1}{3}\).

A larger number of smaller-sized cloud droplets in a polluted cloud increase the total cloud droplet surface area, which results in the cloud albedo increasing (i.e. the cloud becoming brighter). A polluted cloud thus reflects more incoming solar radiation back to space and causes a negative radiative forcing since the fraction of solar radiation that is reflected back into space is proportional to the total cross-section area of droplets per projected area of cloud. This effect is referred to as 1st indirect effect or Twomey effect after [Twomey, 1974; 1977].

In general, the greater the number concentration of aerosol particles, the greater \(N_d\) but this relationship is non-linear because cloud droplets suppress supersaturation more effectively at high number concentrations, which prevents smaller-sized particles from activating. The resultant change in cloud albedo, \(\dd{A}\), due to a change in \(N_d\) is referred to as cloud albedo susceptibility [Platnick and Twomey, 1994], which can be derived by combining Twomey’s expression for optical thickness with a simple expression for cloud albedo: \(A \approx \frac{\tau}{\tau + 7}\)

(5.34)#\[\dv{A}{N_d} \approx \frac{A(1-A)}{3N_d}\]

The increase in CCN number concentration can lead to further microphysical adjustments of a cloud. One such adjustment that mostly affects boundary-layer clouds is referred to as cloud lifetime effect or Albrecht effect [Albrecht, 1989]. This refers to a slow-down of the formation of precipitation in polluted clouds because smaller droplets have a lower collision efficiency than larger droplets, and the lack or delay of precipitation prolongs the lifetime of a cloud. In more detail, in clean conditions a cloud consists of a few cloud droplets with a broader size distribution that contains at least some large particles that efficiently initiate the collision-coalescence process and form drizzle and rain (as opposed to a polluted cloud with a narrow size distribution and a lack of large particles).

A further indirect effect of aerosol on clouds that has been suggested is related to the presence of absorbing (e.g. black carbon) aerosol particles (previously referred to as semi-indirect effect). Absorbing aerosol particles may heat the lowermost atmosphere and thus suppress cumulus cloud formation (or burn off these clouds) by stabilizing the boundary layer and reducing relative humidity.