Earth and Space Sciences · Book 4 · Bachelor Year 2

University Earth and Space Sciences — Year 2

University Earth and Space Sciences — Year 2 · Bachelor Year 2

21Atmospheric Thermodynamics

On 16 April 2026, at 12 UTC, a balloon rose from the weather station of Norman, Oklahoma, as balloons rise from more than eight hundred stations at that hour every day. The instrument it carried radioed back a vertical profile that a forecaster reads at a glance: air near the ground almost saturated, a thin warm lid six hundred metres up, and above it a temperature falling fast for ten kilometres. Lifted through the lid, the surface air would have found itself warmer than its surroundings over most of the troposphere, with some two thousand joules per kilogram to spend on rising. This chapter learns to read that energy from a sounding: how much water air holds, how a rising parcel cools before and after it condenses, how the diagrams of meteorologists turn areas into energies, and how cloud droplets form at all.

You already know

From the physics books: saturation vapour pressure, relative humidity, dew point, latent heat, the Clausius–Clapeyron relation and its integrated form; nucleation and supercooling; the first law, adiabatic changes and the heat capacities of an ideal gas. From the chemistry books: partial pressures, Raoult’s law, surface tension and the Kelvin equation. From the school volumes: the condensation level and its rule of 125 m125\,\mathrm{m} per degree, stable and unstable air, the condensation nucleus, the warm-rain and ice-crystal processes. From the Year 1 volume: the standard atmosphere, the environmental and dry adiabatic lapse rates (Γd=9.8 K/km\Gamma_d = 9.8\,\mathrm{K}/\mathrm{km}), the air parcel, the potential temperature, static stability, temperature inversions and the hypsometric equation.

Throughout the chapter: Rd=R/Md=287.0 J kg−1 K−1R_d = R/M_d = 287.0\,\mathrm{J}\,\mathrm{kg}^{-1}\,\mathrm{K}^{-1} for dry air of molar mass 28.97 g/mol28.97\,\mathrm{g}/\mathrm{mol}; Rv=R/Mw=461.5 J kg−1 K−1R_v = R/M_w = 461.5\,\mathrm{J}\,\mathrm{kg}^{-1}\,\mathrm{K}^{-1} for water vapour, of molar mass 18.0153 g/mol18.0153\,\mathrm{g}/\mathrm{mol}; ε=Mw/Md=0.622\varepsilon = M_w/M_d = 0.622; cp=72Rd=1004.5 J kg−1 K−1c_p = \tfrac72 R_d = 1004.5\,\mathrm{J}\,\mathrm{kg}^{-1}\,\mathrm{K}^{-1} (the Year 1 volume, with the ratio of heat capacities 1.4); g=9.81 m s−2g = 9.81\,\mathrm{m}\,\mathrm{s}^{-2}; and the latent heat of vaporization at the triple point, L=2.501×106 J/kgL = 2.501 \times 10^{6}\,\mathrm{J}/\mathrm{kg}.

21.1 Moist air

Air is a mixture of dry air and water vapour, both close to ideal gases. The physics books measure its water by the relative humidity and the dew point; a rising parcel needs measures that do not change when it expands.

Definition 21.1 (Vapour pressure; mixing ratio; specific humidity)

The vapour pressure ee of moist air is the partial pressure of its water vapour. Its mixing ratio is r=mv/mdr = m_v/m_d, the mass of water vapour per unit mass of dry air in the same volume, and its specific humidity is q=mv/(md+mv)q = m_v/(m_d + m_v), the mass of vapour per unit mass of moist air. Both are usually given in grams per kilogram.

Proposition 21.2 (Humidity variables)

In moist air at the total pressure pp with the vapour pressure ee,

r=εep−e,q=r1+r=εep−(1−ε)e.r = \frac{\varepsilon e}{p - e}, \qquad q = \frac{r}{1 + r} = \frac{\varepsilon e}{p - (1 - \varepsilon)e} .

At saturation, e=es(T)e = e_s(T) and r=rs(T,p)r = r_s(T, p), the saturation mixing ratio. Neither rr nor qq changes when a parcel expands or is compressed without gaining or losing water.

Proof. In the same volume VV at the temperature TT, each gas obeys its own ideal-gas law with its partial pressure: mv=eV/(RvT)m_v = eV/(R_vT) and md=(p−e)V/(RdT)m_d = (p - e)V/(R_dT), so r=(Rd/Rv) e/(p−e)r = (R_d/R_v)\,e/(p - e), and Rd/Rv=Mw/Md=εR_d/R_v = M_w/M_d = \varepsilon. Dividing mvm_v by md+mvm_d + m_v gives q=r/(1+r)q = r/(1 + r). Expansion changes VV, pp and ee but not the two masses. ∎

Definition 21.3 (Virtual temperature)

The virtual temperature TvT_v of moist air is the temperature that dry air would need, at the same pressure, to have the same density:

Tv=T 1+r/ε1+r≈T (1+0.61 r).T_v = T\,\frac{1 + r/\varepsilon}{1 + r} \approx T\,(1 + 0.61\,r) .

The proof that this temperature does what it says is a line: the density of moist air is (md+mv)/V=(p−e)/(RdT)+e/(RvT)(m_d + m_v)/V = \bigl(p - e\bigr)/(R_dT) + e/(R_vT), which, with e=pr/(ε+r)e = pr/(\varepsilon + r), equals p/(RdTv)p/(R_dT_v). Water vapour is lighter than dry air, so moist air is less dense than dry air at the same temperature and pressure, as if it were a little warmer.

Example 21.4 (The air at Norman)

At the ground the sounding gives p=971.1 hPap = 971.1\,\mathrm{hPa}, T=19.0 ∘CT = 19.0\,{}^{\circ}\mathrm{C} and a dew point of 18.5 ∘C18.5\,{}^{\circ}\mathrm{C}. The vapour pressure is the saturation vapour pressure at the dew point, e=es(18.5 ∘C)=21.5 hPae = e_s(18.5\,{}^{\circ}\mathrm{C}) = 21.5\,\mathrm{hPa} (Proposition 21.5 below). Then r=0.622×21.5/(971.1−21.5)=14.1 g/kgr = 0.622 \times 21.5/(971.1 - 21.5) = 14.1\,\mathrm{g}/\mathrm{kg}, q=13.9 g/kgq = 13.9\,\mathrm{g}/\mathrm{kg}, and Tv−T=0.61×0.0141×292.15 K=2.5 KT_v - T = 0.61 \times 0.0141 \times 292.15\,\mathrm{K} = 2.5\,\mathrm{K}: this air is as light as dry air two and a half degrees warmer.

21.2 Saturation vapour pressure

Over a flat surface of pure water, the vapour pressure at equilibrium depends only on the temperature. The physics books derive its law, the Clausius–Clapeyron relation, and integrate it with a constant latent heat; for water vapour, ML/R=L/RvM L/R = L/R_v.

Proposition 21.5 (Saturation vapour pressure)

With the latent heat LL held constant, the saturation vapour pressure over water is

es(T)=etexp⁡ ⁣[LRv(1Tt−1T)],et=6.117 hPa,Tt=273.16 K,LRv=5419 K,e_s(T) = e_t \exp\!\Bigl[\frac{L}{R_v}\Bigl(\frac{1}{T_t} - \frac{1}{T}\Bigr)\Bigr], \qquad e_t = 6.117\,\mathrm{hPa}, \quad T_t = 273.16\,\mathrm{K}, \quad \frac{L}{R_v} = 5419\,\mathrm{K},

where (Tt,et)(T_t, e_t) is the triple point of water.

Corollary 21.6 (Seven per cent per kelvin)

The relative rate of increase of the saturation vapour pressure is dln⁡es/dT=L/(RvT2)\mathrm{d}\ln e_s/\mathrm{d}T = L/(R_v T^2): 7.3 %7.3\,\% per kelvin at 0 ∘C0\,{}^{\circ}\mathrm{C}, 6.3 %6.3\,\% at 20 ∘C20\,{}^{\circ}\mathrm{C}, 5.9 %5.9\,\% at 30 ∘C30\,{}^{\circ}\mathrm{C}.

The constant latent heat is the formula’s only flaw: the latent heat falls from 2500.9 kJ/kg2500.9\,\mathrm{kJ}/\mathrm{kg} at the triple point to 2453.5 kJ/kg2453.5\,\mathrm{kJ}/\mathrm{kg} at 20 ∘C20\,{}^{\circ}\mathrm{C}, so the formula runs ahead of the measured values as the temperature rises, by 1.2 %1.2\,\% at 20 ∘C20\,{}^{\circ}\mathrm{C} and 2.5 %2.5\,\% at 30 ∘C30\,{}^{\circ}\mathrm{C} (figure below). An empirical formula fitted to measurements, es=6.112 hPa×exp⁡[17.67 t/(t+243.5)]e_s = 6.112\,\mathrm{hPa} \times \exp[17.67\,t/(t + 243.5)] with tt in degrees Celsius, stays within a few parts in a thousand; data products use it, but the proposition is the one that can be differentiated and understood, and this chapter keeps it.

Below 0 ∘C0\,{}^{\circ}\mathrm{C}, water can stay liquid in clouds, supercooled, and two saturation vapour pressures exist: over liquid water and over ice. The second is lower, because the latent heat of sublimation exceeds that of vaporization by the latent heat of fusion. Air saturated with respect to supercooled droplets is supersaturated with respect to ice, by 10 %10\,\% at −10 ∘C-10\,{}^{\circ}\mathrm{C} and 22 %22\,\% at −20 ∘C-20\,{}^{\circ}\mathrm{C}: the last section returns to that difference.

Left: the saturation vapour pressure over water by , by the empirical formula and as measured. Right: the ratio of the saturation vapour pressures over supercooled water (empirical formula) and over ice (IAPWS sublimation equation).
Left: the saturation vapour pressure over water by Proposition 21.5, by the empirical formula and as measured. Right: the ratio of the saturation vapour pressures over supercooled water (empirical formula) and over ice (IAPWS sublimation equation).

21.3 Rising moist air

A parcel that rises without exchanging heat cools at the dry adiabatic lapse rate, keeps its potential temperature, and, if it does not condense, keeps its mixing ratio. Its vapour pressure e=pr/(ε+r)e = pr/(\varepsilon + r) falls with the pressure, but its saturation vapour pressure falls faster, with its temperature. The two meet at a level that the school volumes found by a rule.

Definition 21.7 (Lifting condensation level)

The lifting condensation level (LCL) of a parcel is the level at which it becomes saturated when it is lifted dry adiabatically. It is the condensation level of the school volumes, the height of the flat bases of the clouds that form by lifting.

Proposition 21.8 (Height of the lifting condensation level)

A parcel at the temperature TT with the dew point TdT_d reaches saturation about

zLCL≈T−TdΓd−Γdew,Γdew=gRvTd2LRdT≈1.8 K/km,z_{\text{LCL}} \approx \frac{T - T_d}{\Gamma_d - \Gamma_{\text{dew}}}, \qquad \Gamma_{\text{dew}} = \frac{g R_v T_d^2}{L R_d T} \approx 1.8\,\mathrm{K}/\mathrm{km},

above its starting level: about 125 m125\,\mathrm{m} for each degree of dew-point depression.

Proof. The parcel’s temperature falls at Γd\Gamma_d. Its dew point is the temperature at which its vapour pressure ee would saturate; with rr constant, ee is proportional to pp, so dln⁡e=dln⁡p=−g dz/(RdT)\mathrm{d}\ln e = \mathrm{d}\ln p = -g\,\mathrm{d}z/(R_dT) by the hydrostatic balance. By Proposition 21.5, dln⁡e=L dTd/(RvTd2)\mathrm{d}\ln e = L\,\mathrm{d}T_d/(R_vT_d^2). Hence dTd/dz=−gRvTd2/(LRdT)=−Γdew\mathrm{d}T_d/\mathrm{d}z = -gR_vT_d^2/(LR_dT) = -\Gamma_{\text{dew}}: the dew point falls too, but more slowly. The gap T−TdT - T_d closes at the rate Γd−Γdew\Gamma_d - \Gamma_{\text{dew}}. With T≈Td≈288 KT \approx T_d \approx 288\,\mathrm{K}, Γdew=1.8 K/km\Gamma_{\text{dew}} = 1.8\,\mathrm{K}/\mathrm{km}, and 1/(9.8−1.8)1/(9.8 - 1.8) kilometres per kelvin is 125 m125\,\mathrm{m}. ∎

At Norman the dew-point depression was 0.5 K0.5\,\mathrm{K}: the surface air would have condensed after 63 m63\,\mathrm{m} of ascent, at 964 hPa964\,\mathrm{hPa}. The full computation, which follows the dry adiabat and the saturation mixing ratio without linearizing, gives 64 m64\,\mathrm{m}.

Above the LCL the parcel condenses as it rises. Condensation releases latent heat, and the parcel cools more slowly.

Definition 21.9 (Saturated adiabatic lapse rate; moist adiabat)

The saturated adiabatic lapse rate Γs\Gamma_s is the rate at which the temperature of a saturated parcel falls as it rises, its vapour condensing to keep it saturated. The curve it follows on a diagram of temperature against pressure is a moist adiabat. When the condensed water is assumed to fall out at once, the ascent is called pseudo-adiabatic.

Proposition 21.10 (Saturated adiabatic lapse rate)

For a pseudo-adiabatic ascent, with rs=rs(T,p)r_s = r_s(T, p),

Γs=g 1+LrsRdTcp+L2rsRvT2.\Gamma_s = g\,\frac{1 + \dfrac{L r_s}{R_d T}}{c_p + \dfrac{L^2 r_s}{R_v T^2}} .

It tends to Γd=g/cp\Gamma_d = g/c_p when rsr_s tends to zero, and is smaller otherwise.

Proof. Per unit mass of dry air, the first law for a parcel that receives the latent heat of the vapour it condenses, −L drs-L\,\mathrm{d}r_s, reads cp dT−(RdT/p) dp=−L drsc_p\,\mathrm{d}T - (R_dT/p)\,\mathrm{d}p = -L\, \mathrm{d}r_s; with the hydrostatic balance (RdT/p) dp=−g dz(R_dT/p)\,\mathrm{d}p = -g\,\mathrm{d}z, it becomes cp dT+g dz+L drs=0c_p\,\mathrm{d}T + g\,\mathrm{d}z + L\,\mathrm{d}r_s = 0. With rs≈εes/pr_s \approx \varepsilon e_s/p, drs=rs (dln⁡es−dln⁡p)=rs(L dT/(RvT2)+g dz/(RdT))\mathrm{d}r_s = r_s\,(\mathrm{d}\ln e_s - \mathrm{d}\ln p) = r_s\bigl(L\,\mathrm{d}T/(R_vT^2) + g\,\mathrm{d}z/(R_dT)\bigr). Substituting and collecting the terms in dT\mathrm{d}T and in dz\mathrm{d}z gives −dT/dz=Γs-\mathrm{d}T/\mathrm{d}z = \Gamma_s. ∎

The rate depends on how much water there is to condense. At 1000 hPa1000\,\mathrm{hPa} it is 6.5 K/km6.5\,\mathrm{K}/\mathrm{km} at 0 ∘C0\,{}^{\circ}\mathrm{C}, 5.2 K/km5.2\,\mathrm{K}/\mathrm{km} at 10 ∘C10\,{}^{\circ}\mathrm{C}, 4.2 K/km4.2\,\mathrm{K}/\mathrm{km} at 20 ∘C20\,{}^{\circ}\mathrm{C} and 3.4 K/km3.4\,\mathrm{K}/\mathrm{km} at 30 ∘C30\,{}^{\circ}\mathrm{C}; high in the troposphere, at 200 hPa200\,\mathrm{hPa} and −60 ∘C-60\,{}^{\circ}\mathrm{C}, it is back to 9.6 K/km9.6\,\mathrm{K}/\mathrm{km}, almost the dry rate. A moist adiabat is therefore steep near warm ground and bends towards the dry adiabats aloft.

Definition 21.11 (Equivalent potential temperature)

The equivalent potential temperature θe\theta_e of a parcel is the temperature it would have if it were lifted until all its vapour had condensed and fallen out, then brought down dry adiabatically to 1000 hPa1000\,\mathrm{hPa}. A good approximation is

θe≈θexp⁡ ⁣(LrcpTLCL),\theta_e \approx \theta \exp\!\Bigl(\frac{L r}{c_p T_{\text{LCL}}}\Bigr),

with θ\theta the potential temperature and TLCLT_{\text{LCL}} the temperature at the parcel’s lifting condensation level.

The potential temperature is conserved in dry ascent only; θe\theta_e is conserved, to this approximation, in dry and saturated ascent alike, because the heat it counts as already released is exactly the heat condensation releases. It labels a moist adiabat as θ\theta labels a dry one. The Norman surface air has θ=292.15×(1000/971.1)2/7=294.6 K\theta = 292.15 \times (1000/971.1)^{2/7} = 294.6\,\mathrm{K} and θe=294.6exp⁡[2.501×106×0.0141/(1004.5×291.5)]=332 K\theta_e = 294.6 \exp[2.501 \times 10^6 \times 0.0141/(1004.5 \times 291.5)] = 332\,\mathrm{K}.

21.4 Thermodynamic diagrams

Definition 21.12 (Atmospheric sounding)

An atmospheric sounding is a vertical profile of the pressure, temperature, humidity and wind of the atmosphere above a place, measured by an instrument carried aloft, usually a radiosonde under a balloon.

A sounding is not quite vertical, because the balloon drifts with the wind, nor quite instantaneous, because it takes several minutes to reach the upper troposphere; and the main network launches only twice a day. Its numbers become useful once they are drawn on a chart where a rising parcel can be followed with a ruler.

Definition 21.13 (Thermodynamic diagram; skew-T log-p diagram; tephigram)

A thermodynamic diagram is a chart whose coordinates are functions of the pressure and the temperature, on which isobars, isotherms, dry adiabats, moist adiabats and lines of constant saturation mixing ratio are printed, and on which a sounding is plotted as two curves, its temperature and its dew point. On a skew-T log-p diagram the vertical coordinate is −ln⁡p-\ln p and the horizontal one is T+kln⁡(p0/p)T + k\ln(p_0/p), so that isobars are horizontal and isotherms are straight lines slanting up to the right. A tephigram uses the temperature and the entropy per unit mass cpln⁡θc_p\ln\theta as coordinates, so that isotherms and dry adiabats are perpendicular straight lines.

Proposition 21.14 (Areas are energies)

On a skew-T log-p diagram, the area enclosed between two curves T1(p)T_1(p) and T2(p)T_2(p) over the pressures pb>pap_b > p_a is proportional to ∫papbRd (T1−T2) dln⁡p\int_{p_a}^{p_b} R_d\,(T_1 - T_2)\,\mathrm{d}\ln p, the energy per unit mass that Proposition 21.19 below gives to a parcel following T1T_1 through an environment T2T_2.

Proof. Take Y=ln⁡(p0/p)Y = \ln(p_0/p) and X=T+kYX = T + kY. The Jacobian of (T,Y)↦(X,Y)(T, Y) \mapsto (X, Y) is 1, so areas in the (X,Y)(X, Y) plane equal areas in the (T,Y)(T, Y) plane, where the area between two curves is ∫(T1−T2) dY=∫(T1−T2) dln⁡p\int (T_1 - T_2)\,\mathrm{d}Y = \int (T_1 - T_2)\,\mathrm{d}\ln p taken from pap_a to pbp_b with the sign fixed. Multiplying by RdR_d and by the scale of the chart gives the energy. ∎

Method 21.15 (Lifting a parcel on a diagram)

  1. Plot the sounding’s temperature and dew point against pressure.
  2. From the parcel’s temperature, follow the dry adiabat upward; from its dew point, follow the line of saturation mixing ratio. Where they meet is the LCL.
  3. From the LCL, follow the moist adiabat upward.
  4. Where the parcel’s curve lies to the right of the sounding’s temperature, the parcel is warmer than its surroundings; the levels where the two curves cross are the levels of Definition 21.17 below.

21.5 Stability and convective energy

The Year 1 volume compared the lapse rate Γ\Gamma of a layer with Γd\Gamma_d. A saturated parcel cools at Γs\Gamma_s instead, and a third case appears between the two rates.

Definition 21.16 (Conditional instability)

A layer whose lapse rate lies between the saturated and the dry adiabatic lapse rates, Γs<Γ<Γd\Gamma_s < \Gamma < \Gamma_d, is in conditional instability: stable for unsaturated parcels, unstable for saturated ones.

At Norman the temperature fell by 37.0 K37.0\,\mathrm{K} between 850 hPa850\,\mathrm{hPa} (1489 m1489\,\mathrm{m}) and 500 hPa500\,\mathrm{hPa} (5725 m5725\,\mathrm{m}), at 8.7 K/km8.7\,\mathrm{K}/\mathrm{km}: conditional instability over four kilometres. Whether the instability is released depends on where a parcel from below becomes saturated and then warmer than its surroundings.

Definition 21.17 (Level of free convection; equilibrium level)

The level of free convection (LFC) of a lifted parcel is the level above which, saturated, it becomes warmer than its surroundings and rises freely. The equilibrium level (EL) is the level above the LFC where its temperature again equals that of its surroundings.

The 12 UTC sounding of Norman, Oklahoma, on 16 April 2026, on a skew-T log-p diagram, with the surface parcel lifted along its dry adiabat to the LCL and its moist adiabat above. Red: the positive area (CAPE); blue: the negative area below the LFC (CIN), both computed with the temperature. Saturation mixing ratios in grams per kilogram, labelled at 600\, hPa.
The 12 UTC sounding of Norman, Oklahoma, on 16 April 2026, on a skew-T log-p diagram, with the surface parcel lifted along its dry adiabat to the LCL and its moist adiabat above. Red: the positive area (CAPE); blue: the negative area below the LFC (CIN), both computed with the temperature. Saturation mixing ratios in grams per kilogram, labelled at 600 hPa600\,\mathrm{hPa}.

Definition 21.18 (Convective available potential energy; convective inhibition)

The convective available potential energy (CAPE) of a parcel is the work of the buoyancy force on unit mass of it between its LFC and its EL,

CAPE=∫zLFCzELg Tv,parcel−Tv,envTv,env dz,\text{CAPE} = \int_{z_{\text{LFC}}}^{z_{\text{EL}}} g\,\frac{T_{v,\text{parcel}} - T_{v,\text{env}}}{T_{v,\text{env}}}\,\mathrm{d}z ,

and its convective inhibition (CIN) is the same integral, negative, from its starting level to its LFC.

Proposition 21.19 (Buoyancy on a diagram; the parcel’s speed)

  1. In terms of pressure, CAPE=Rd∫pELpLFC(Tv,parcel−Tv,env) dln⁡p\text{CAPE} = R_d \int_{p_{\text{EL}}}^{p_{\text{LFC}}} (T_{v,\text{parcel}} - T_{v,\text{env}})\,\mathrm{d}\ln p: on a skew-T log-p diagram, it is proportional to the area between the parcel’s curve and the sounding.
  2. A parcel that leaves the LFC at rest and keeps all the work of its buoyancy reaches the EL with the vertical speed wmax⁡=2 CAPEw_{\max} = \sqrt{2\,\text{CAPE}}.

Proof. 1. The parcel and its surroundings share the pressure, so the density ratio is the ratio of the virtual temperatures and the buoyancy force per unit mass is g(Tv,parcel−Tv,env)/Tv,envg(T_{v,\text{parcel}} - T_{v,\text{env}})/T_{v,\text{env}}. The environment is hydrostatic: dz=−RdTv,env dln⁡p/g\mathrm{d}z = -R_dT_{v,\text{env}}\,\mathrm{d}\ln p/g. Substituting cancels gg and Tv,envT_{v,\text{env}}; the area follows from Proposition 21.14. 2. With bb the buoyancy force per unit mass, dw/dt=b\mathrm{d}w/\mathrm{d}t = b and w=dz/dtw = \mathrm{d}z/ \mathrm{d}t give w dw=b dzw\,\mathrm{d}w = b\,\mathrm{d}z; integrating from the LFC, where w=0w = 0, to the EL gives w2/2=CAPEw^2/2 = \text{CAPE}. ∎

Example 21.20 (The energy in the Norman sounding)

Lift the surface air of Example 21.4. It saturates at 964 hPa964\,\mathrm{hPa}, 64 m64\,\mathrm{m} up, and follows a moist adiabat. Below 900 hPa900\,\mathrm{hPa} the sounding is itself nearly saturated and the parcel stays a few tenths of a degree colder than its surroundings; then comes the inversion between 908 hPa908\,\mathrm{hPa} and 889 hPa889\,\mathrm{hPa}, where the temperature rises by 4.1 K4.1\,\mathrm{K} in 178 m178\,\mathrm{m} and the parcel finds itself 5 K5\,\mathrm{K} too cold. It becomes warmer than its surroundings only at 732 hPa732\,\mathrm{hPa}, 2.4 km2.4\,\mathrm{km} above the ground, its LFC, and stays warmer up to its EL at 236 hPa236\,\mathrm{hPa}, 10.9 km10.9\,\mathrm{km} above sea level, just below the tropopause the sounding reports at 235 hPa235\,\mathrm{hPa}. On a fine grid of pressures the areas of the skew-T figure are CAPE =2100 J/kg= 2100\,\mathrm{J}/\mathrm{kg} and CIN =−200 J/kg= -200\,\mathrm{J}/\mathrm{kg} with the temperature, and 2230 J/kg2230\,\mathrm{J}/\mathrm{kg} and −130 J/kg-130\,\mathrm{J}/\mathrm{kg} with the virtual temperature: the moisture of the parcel adds its own buoyancy. Parcel theory gives wmax⁡=2×2100=65 m/sw_{\max} = \sqrt{2 \times 2100} = 65\,\mathrm{m}/\mathrm{s}.

Forecasters’ tables call a CAPE up to 1500 J/kg1500\,\mathrm{J}/\mathrm{kg} positive, one between 1500 and 2500 J/kg2500\,\mathrm{J}/\mathrm{kg} large and one above 2500 J/kg2500\,\mathrm{J}/\mathrm{kg} extreme; values above 1000 J/kg1000\,\mathrm{J}/\mathrm{kg} are common in thunderstorm environments, and rare ones exceed 5000 J/kg5000\,\mathrm{J}/\mathrm{kg}. CAPE is largest in the warm season, especially in late spring. The lifted index, the sounding’s temperature at 500 hPa500\,\mathrm{hPa} minus that of the lifted parcel, sums the same information in one number: at Norman, −18.4−(−8.9)=−9.5 K-18.4 - (-8.9) = -9.5\,\mathrm{K}, in the class from −8-8 to −10-10 that the forecasters’ tables call extreme instability. None of these values is a threshold above which storms must form.

The CIN is the price of admission. The warm lid kept the surface air from rising on its own: a parcel would have to be pushed through it, by a front or by the slow ascent of the air ahead of a disturbance aloft, or warmed from below by the afternoon sun until it was buoyant all the way. And the speed wmax⁡w_{\max} is an upper bound: a real updraft carries its condensed water, which weighs it down, and mixes with the drier air around it, which cools it by evaporation. The equilibrium level remains a good guide to where the storm tops spread out into an anvil.

A thunderstorm over Nebraska, 20 April 1985: its updraft spreads out into an anvil where it reaches its equilibrium level, near the tropopause. (Stephen Corfidi, NOAA/NWS, public domain, via Wikimedia Commons.)
A thunderstorm over Nebraska, 20 April 1985: its updraft spreads out into an anvil where it reaches its equilibrium level, near the tropopause. (Stephen Corfidi, NOAA/NWS, public domain, via Wikimedia Commons.)

21.6 Cloud microphysics

A parcel at its LCL is saturated over a flat water surface. A droplet is not flat: the physics books show that a new phase does not appear at equilibrium, but only through nucleation, and the chemistry books that the vapour pressure over a curved surface exceeds that over a flat one (the Kelvin equation). With the constant AA of Proposition 21.22 below, a droplet of pure water 2 nm2\,\mathrm{nm} across is in equilibrium only at the saturation ratio exp⁡(2.1/2)=2.9\exp(2.1/2) = 2.9, a supersaturation of 190 %190\,\%, far above the fractions of a per cent at which particles activate: droplets form on particles.

Definition 21.21 (Supersaturation; cloud condensation nucleus)

The supersaturation of air is s=e/es−1s = e/e_s - 1, the excess of its vapour pressure over the saturation vapour pressure above a flat surface of pure water, usually given in per cent. A cloud condensation nucleus (CCN) is an aerosol particle on which a cloud droplet forms at the small supersaturations of clouds. It is the condensation nucleus of the school volumes.

Proposition 21.22 (Köhler curve)

A solution droplet of diameter DD grown on a dry particle of diameter DdD_d and hygroscopicity κ\kappa is in equilibrium with the saturation ratio

S(D)=D3−Dd3D3−Dd3(1−κ) exp⁡(AD),A=4σMwRTρw,S(D) = \frac{D^3 - D_d^3}{D^3 - D_d^3(1 - \kappa)}\,\exp\Bigl(\frac{A}{D}\Bigr), \qquad A = \frac{4\sigma M_w}{R T\rho_w},

with σ\sigma the surface tension and ρw\rho_w the density of water (admitted). The curve S(D)S(D) is the Köhler curve. For D≫DdD \gg D_d,

s=S−1≈AD−κDd3D3,s = S - 1 \approx \frac{A}{D} - \frac{\kappa D_d^3}{D^3},

which is maximal at the critical diameter and supersaturation

Dc=3κDd3A,sc=4A327 κDd3.D_c = \sqrt{\frac{3\kappa D_d^3}{A}}, \qquad s_c = \sqrt{\frac{4A^3}{27\,\kappa D_d^3}} .

Proof of the approximation. The exponential is the Kelvin factor; the fraction is the water activity, a version of Raoult’s law in which κ\kappa measures how much water the dissolved matter binds. For D≫DdD \gg D_d and A≪DA \ll D, exp⁡(A/D)≈1+A/D\exp(A/D) \approx 1 + A/D and the fraction is 1−κDd3/D31 - \kappa D_d^3/D^3 to first order; the product, to first order, is 1+A/D−κDd3/D31 + A/D - \kappa D_d^3/D^3. Setting the derivative −A/D2+3κDd3/D4-A/D^2 + 3\kappa D_d^3/D^4 to zero gives Dc2=3κDd3/AD_c^2 = 3\kappa D_d^3/A, and then sc=A/Dc−A/(3Dc)=2A/(3Dc)s_c = A/D_c - A/(3D_c) = 2A/(3D_c), which is the formula. ∎

With σ=0.072 J/m2\sigma = 0.072\,\mathrm{J}/\mathrm{m}^{2} at 298.15 K298.15\,\mathrm{K}, A=2.1 nmA = 2.1\,\mathrm{nm}. Ammonium sulfate has κ=0.61\kappa = 0.61; sodium chloride is more hygroscopic still (κ=1.28\kappa = 1.28), and atmospheric particles mostly lie between 0.1 and 0.9. A dry particle of 50 nm50\,\mathrm{nm} of ammonium sulfate has Dc=0.33 µmD_c = 0.33\,\text{µ}\mathrm{m} and sc=0.42 %s_c = 0.42\,\%; one of 100 nm100\,\mathrm{nm}, 0.93 µm0.93\,\text{µ}\mathrm{m} and 0.15 %0.15\,\%. The approximation differs from the exact maximum by less than 1 %1\,\%.

Köhler curves of ammonium sulfate particles (= 0.61) of three dry diameters at 298.15\, K, with their critical points (dots). In air at 0.5\,\% supersaturation, the 50\, nm and 100\, nm particles grow without limit; the 30\, nm particle stays a haze droplet.
Köhler curves of ammonium sulfate particles (κ=0.61\kappa = 0.61) of three dry diameters at 298.15 K298.15\,\mathrm{K}, with their critical points (dots). In air at 0.5 %0.5\,\% supersaturation, the 50 nm50\,\mathrm{nm} and 100 nm100\,\mathrm{nm} particles grow without limit; the 30 nm30\,\mathrm{nm} particle stays a haze droplet.

The curve divides the particles. Below its peak a droplet is in stable equilibrium: if it grows, the air is no longer supersaturated enough for it and it shrinks back; such haze droplets swell and shrink with the humidity. Once the supersaturation of the air exceeds scs_c, the droplet is activated: past the peak, it needs less and less supersaturation as it grows, and it grows for as long as there is vapour. A rising parcel cools, its supersaturation climbs, and the largest and most soluble particles activate first; the droplets they make then take up the vapour and cap the supersaturation. Over remote oceans there are about 110 CCN per cubic centimetre active at 0.4 %0.4\,\%, over remote continents almost twice as many, and in polluted air about ten times more. The same water shared among more droplets makes them smaller.

Remark 21.23 (Why condensation alone does not make rain)

A droplet grows by the diffusion of vapour towards it. The flux through a sphere around a droplet of radius rr is 4πrDv Δρv4\pi r D_v\,\Delta\rho_v, with DvD_v the diffusivity of vapour and Δρv\Delta\rho_v the excess of vapour density far away, so its mass 43πr3ρw\tfrac43\pi r^3\rho_w grows at a rate proportional to rr: r dr/dtr\,\mathrm{d}r/\mathrm{d}t is constant and r2r^2 grows linearly in time. Large droplets grow more slowly than small ones, and a population growing by condensation alone becomes narrower, not broader. Raindrops form when droplets of different sizes, falling at different speeds, collide and coalesce.

Above the 0 ∘C0\,{}^{\circ}\mathrm{C} level a cloud is made of supercooled droplets until ice nucleating particles start ice crystals. Then the difference of the right panel of the saturation figure acts: air saturated over the droplets is 10 %10\,\% supersaturated over ice at −10 ∘C-10\,{}^{\circ}\mathrm{C}, tens of times the supersaturations at which the droplets activated. The crystals grow by deposition at the expense of the droplets, which evaporate (the Bergeron process, the ice-crystal process of the school volumes), and they collect droplets that freeze on them (riming).

The two roads from vapour to precipitation: activation and growth of droplets, then collision and coalescence (top); ice crystals that grow at the expense of supercooled droplets, rime, and melt as they fall into warmer air (bottom).
The two roads from vapour to precipitation: activation and growth of droplets, then collision and coalescence (top); ice crystals that grow at the expense of supercooled droplets, rime, and melt as they fall into warmer air (bottom).

In the field — A morning launch

A radiosonde released at dawn from a weather station (AI-generated illustration).

At each upper-air station the technician prepares the balloon and the radiosonde in the hour before 00 and 12 UTC, the times at which more than eight hundred stations of the world observe together. The instrument, of 60 to 80 g80\,\mathrm{g}, hangs below the balloon and reports the pressure, the temperature, the humidity and its satellite position every second. The balloon rises about 300 m300\,\mathrm{m} a minute, for two hours or more, past 35 km35\,\mathrm{km} before it bursts, and may drift more than 300 km300\,\mathrm{km}. Forecasters plot the profile on a skew-T diagram, and on days when severe weather is expected some stations launch extra balloons.

Exercises

Data for the exercises, unless stated otherwise: es(T)=6.117 hPa×exp⁡[5419 K×(1/273.16 K−1/T)]e_s(T) = 6.117\,\mathrm{hPa} \times \exp[5419\,\mathrm{K} \times (1/273.16\,\mathrm{K} - 1/T)]; ε=0.622\varepsilon = 0.622, Rd=287.0 J kg−1 K−1R_d = 287.0\,\mathrm{J}\,\mathrm{kg}^{-1}\,\mathrm{K}^{-1}, Rv=461.5 J kg−1 K−1R_v = 461.5\,\mathrm{J}\,\mathrm{kg}^{-1}\,\mathrm{K}^{-1}, cp=1004.5 J kg−1 K−1c_p = 1004.5\,\mathrm{J}\,\mathrm{kg}^{-1}\,\mathrm{K}^{-1}, L=2.501×106 J/kgL = 2.501 \times 10^{6}\,\mathrm{J}/\mathrm{kg}, g=9.81 m s−2g = 9.81\,\mathrm{m}\,\mathrm{s}^{-2}, Γd=9.8 K/km\Gamma_d = 9.8\,\mathrm{K}/\mathrm{km}.

Exercise 21.1 ★

At 850 hPa850\,\mathrm{hPa} the temperature is 10.0 ∘C10.0\,{}^{\circ}\mathrm{C} and the dew point 5.0 ∘C5.0\,{}^{\circ}\mathrm{C}. (a) Find the vapour pressure and the saturation vapour pressure. (b) Find the relative humidity. (c) Find the mixing ratio and the specific humidity, in grams per kilogram.

Solution

Solution of Exercise 21.1.

(a) e=es(278.15 K)=8.73 hPae = e_s(278.15\,\mathrm{K}) = 8.73\,\mathrm{hPa} and es(283.15 K)=12.32 hPae_s(283.15\,\mathrm{K}) = 12.32\,\mathrm{hPa}. (b) 8.73/12.32=70.9 %8.73/12.32 = 70.9\,\%. (c) r=0.622×8.73/(850−8.73)=6.45 g/kgr = 0.622 \times 8.73/(850 - 8.73) = 6.45\,\mathrm{g}/\mathrm{kg} and q=r/(1+r)=6.41 g/kgq = r/(1 + r) = 6.41\,\mathrm{g}/\mathrm{kg}.

Exercise 21.2 ★

(a) By how many per cent per kelvin does ese_s rise at 0 ∘C0\,{}^{\circ}\mathrm{C}, 20 ∘C20\,{}^{\circ}\mathrm{C} and 30 ∘C30\,{}^{\circ}\mathrm{C}? (b) Compute ese_s at 20 ∘C20\,{}^{\circ}\mathrm{C} and at 25 ∘C25\,{}^{\circ}\mathrm{C}, and their ratio; compare with 1.06351.063^5. (c) The measured value at 20 ∘C20\,{}^{\circ}\mathrm{C} is 23.39 hPa23.39\,\mathrm{hPa}. By how much does the formula err, and in which direction? Why?

Solution

Solution of Exercise 21.2.

(a) dln⁡es/dT=5419/T2\mathrm{d}\ln e_s/\mathrm{d}T = 5419/T^2: 7.26 %7.26\,\% per kelvin at 273.15 K273.15\,\mathrm{K}, 6.31 %6.31\,\% at 293.15 K293.15\,\mathrm{K}, 5.90 %5.90\,\% at 303.15 K303.15\,\mathrm{K}. (b) es(20 ∘C)=23.66 hPae_s(20\,{}^{\circ}\mathrm{C}) = 23.66\,\mathrm{hPa} and es(25 ∘C)=32.26 hPae_s(25\,{}^{\circ}\mathrm{C}) = 32.26\,\mathrm{hPa}; ratio 1.363, against 1.0635=1.3571.063^5 = 1.357. The rate is a derivative of ln⁡es\ln e_s, so over 5 K5\,\mathrm{K} the factor is exp⁡(5×0.0620)=1.363\exp(5 \times 0.0620) = 1.363, with 0.0620 the rate at the mid-point 22.5 ∘C22.5\,{}^{\circ}\mathrm{C}; compounding 6.3 %6.3\,\% five times is close but not exact. (c) 23.66/23.39=1.01223.66/23.39 = 1.012: the formula is 1.2 %1.2\,\% too high, because it keeps the latent heat of the triple point, while the latent heat decreases as the temperature rises; the formula’s slope is too steep, and its error grows with the distance from 0 ∘C0\,{}^{\circ}\mathrm{C}.

Exercise 21.3 ★

At noon the surface air is at 30 ∘C30\,{}^{\circ}\mathrm{C} with a dew point of 14 ∘C14\,{}^{\circ}\mathrm{C}. (a) At what height above the ground do cumulus bases form? (b) What is the temperature of a parcel there, by the dry adiabat? (c) What is its dew point there, with Γdew=1.8 K/km\Gamma_{\text{dew}} = 1.8\,\mathrm{K}/\mathrm{km}? Comment. (d) Later the air warms to 33 ∘C33\,{}^{\circ}\mathrm{C}, the dew point unchanged. Where are the cloud bases now?

Solution

Solution of Exercise 21.3.

(a) 125×(30−14)=2000 m125 \times (30 - 14) = 2000\,\mathrm{m}. (b) 30−9.8×2.0=10.4 ∘C30 - 9.8 \times 2.0 = 10.4\,{}^{\circ}\mathrm{C}. (c) 14−1.8×2.0=10.4 ∘C14 - 1.8 \times 2.0 = 10.4\,{}^{\circ}\mathrm{C}: temperature and dew point meet, as they must at the level of saturation; the rule is this calculation. (d) 125×(33−14)=2375 m125 \times (33 - 14) = 2375\,\mathrm{m}: the bases rise as the afternoon warms the ground air.

Exercise 21.4 ★

Read the skew-T diagram of the Norman sounding. (a) Up to which pressure is the air near the ground saturated, or nearly so? (b) Where is the inversion, and how does it show? (c) Read the LCL, the LFC and the EL of the surface parcel. (d) How do the dew point and the temperature compare just above the inversion, and what does that say about this air?

Solution

Solution of Exercise 21.4.

(a) The dew-point curve lies on the temperature curve from the ground to about 908 hPa908\,\mathrm{hPa}: the lowest half-kilometre is saturated, or nearly. (b) Between about 908 and 889 hPa889\,\mathrm{hPa} the temperature curve turns to the right, to warmer values, as the pressure decreases: the temperature rises with height. (c) LCL at about 964 hPa964\,\mathrm{hPa}, a few tens of metres up; LFC at about 730 hPa730\,\mathrm{hPa}; EL at about 235 hPa235\,\mathrm{hPa}. (d) At the top of the inversion the dew point is 11 K11\,\mathrm{K} below the temperature, and 20 K20\,\mathrm{K} below it at 850 hPa850\,\mathrm{hPa}: warm, very dry air lies on top of cool, moist air. Its relative humidity is low, so it did not acquire its warmth by condensation; air that sinks warms by compression at the dry adiabatic rate while keeping its mixing ratio, which would give such a layer.

Exercise 21.5 ★★

Tropical air at 1000 hPa1000\,\mathrm{hPa} has T=30.0 ∘CT = 30.0\,{}^{\circ}\mathrm{C} and r=20 g/kgr = 20\,\mathrm{g}/\mathrm{kg}. (a) Find its virtual temperature exactly and with the approximation T(1+0.61 r)T(1 + 0.61\,r). (b) Find the densities of this air and of dry air at the same temperature and pressure. (c) By how much would dry air have to be warmed to be as light? (d) Why does this matter for convection over warm seas?

Solution

Solution of Exercise 21.5.

(a) Tv=303.15×(1+0.020/0.622)/1.020=306.76 KT_v = 303.15 \times (1 + 0.020/0.622)/1.020 = 306.76\,\mathrm{K}; the approximation gives 303.15×(1+0.61×0.020)=306.85 K303.15 \times (1 + 0.61 \times 0.020) = 306.85\,\mathrm{K}. (b) ρ=p/(RdTv)=105/(287.0×306.76)=1.1358 kg/m3\rho = p/(R_dT_v) = 10^5/(287.0 \times 306.76) = 1.1358\,\mathrm{kg}/\mathrm{m}^{3}, against 105/(287.0×303.15)=1.1494 kg/m310^5/(287.0 \times 303.15) = 1.1494\,\mathrm{kg}/\mathrm{m}^{3} for dry air: 1.2 %1.2\,\% lighter. (c) By Tv−T=3.6 KT_v - T = 3.6\,\mathrm{K}. (d) A rising parcel needs only to be lighter than its surroundings: over warm seas the moisture of the air adds a buoyancy equivalent to several kelvin, as much as a sizable excess of temperature.

Exercise 21.6 ★★

Use Proposition 21.10 at 1000 hPa1000\,\mathrm{hPa}, with rsr_s from Proposition 21.2. (a) Compute rsr_s and Γs\Gamma_s at 0 ∘C0\,{}^{\circ}\mathrm{C}, 20 ∘C20\,{}^{\circ}\mathrm{C} and 30 ∘C30\,{}^{\circ}\mathrm{C}. (b) How much warmer than a dry parcel is a saturated parcel at 20 ∘C20\,{}^{\circ}\mathrm{C} after rising 1 km1\,\mathrm{km}, if its lapse rate stays the same? (c) Why does the difference shrink higher up?

Solution

Solution of Exercise 21.6.

(a) ese_s = 6.11, 23.66 and 43.54 hPa43.54\,\mathrm{hPa}; rs=0.622 es/(1000−es)r_s = 0.622\,e_s/(1000 - e_s) = 3.82, 15.07 and 28.31 g/kg28.31\,\mathrm{g}/\mathrm{kg}; Γs\Gamma_s = 6.48, 4.20 and 3.44 K/km3.44\,\mathrm{K}/\mathrm{km}. (b) The dry parcel cools by 9.8 K9.8\,\mathrm{K}, the saturated one by 4.2 K4.2\,\mathrm{K}: the saturated parcel is 5.6 K5.6\,\mathrm{K} warmer. (c) Higher up the air is colder, rsr_s is smaller, little vapour is left to condense, and Γs\Gamma_s returns towards Γd\Gamma_d.

Exercise 21.7 ★★

The Norman surface parcel (p=971.1 hPap = 971.1\,\mathrm{hPa}, T=19.0 ∘CT = 19.0\,{}^{\circ}\mathrm{C}, Td=18.5 ∘CT_d = 18.5\,{}^{\circ}\mathrm{C}) has its LCL at 18.35 ∘C18.35\,{}^{\circ}\mathrm{C}. (a) Compute its mixing ratio, its potential temperature and its equivalent potential temperature. (b) Lifted, it is at −8.9 ∘C-8.9\,{}^{\circ}\mathrm{C} at 500 hPa500\,\mathrm{hPa}, saturated. Compute its θ\theta and θe\theta_e there. (c) Which of the two is conserved?

Solution

Solution of Exercise 21.7.

(a) e=es(291.65 K)=21.52 hPae = e_s(291.65\,\mathrm{K}) = 21.52\,\mathrm{hPa}, r=0.622×21.52/(971.1−21.52)=14.1 g/kgr = 0.622 \times 21.52/(971.1 - 21.52) = 14.1\,\mathrm{g}/\mathrm{kg}; θ=292.15×(1000/971.1)0.2857=294.6 K\theta = 292.15 \times (1000/971.1)^{0.2857} = 294.6\,\mathrm{K}; θe=294.6exp⁡[2.501×106×0.0141/(1004.5×291.5)]=294.6×1.128=332.3 K\theta_e = 294.6 \exp[2.501 \times 10^6 \times 0.0141/(1004.5 \times 291.5)] = 294.6 \times 1.128 = 332.3\,\mathrm{K}. (b) es(264.25 K)=3.13 hPae_s(264.25\,\mathrm{K}) = 3.13\,\mathrm{hPa}, rs=0.622×3.13/(500−3.13)=3.92 g/kgr_s = 0.622 \times 3.13/(500 - 3.13) = 3.92\,\mathrm{g}/\mathrm{kg}; θ=264.25×20.2857=322.1 K\theta = 264.25 \times 2^{0.2857} = 322.1\,\mathrm{K}; the parcel is saturated, so TLCL=TT_{\text{LCL}} = T and θe=322.1exp⁡[2.501×106×0.00392/(1004.5×264.25)]=334.3 K\theta_e = 322.1 \exp[2.501 \times 10^6 \times 0.00392/(1004.5 \times 264.25)] = 334.3\,\mathrm{K}. (c) θ\theta rose by 27.5 K27.5\,\mathrm{K}, warmed by the latent heat of the 10.2 g/kg10.2\,\mathrm{g}/\mathrm{kg} of vapour that condensed; θe\theta_e changed by 0.6 %0.6\,\%, within the accuracy of its approximate formula: it is conserved.

Exercise 21.8 ★★

A parcel is warmer than its surroundings by 4 K4\,\mathrm{K} between 700 hPa700\,\mathrm{hPa}, its LFC, and 300 hPa300\,\mathrm{hPa}, its EL; ignore the virtual-temperature correction. (a) Compute its CAPE. (b) Compute wmax⁡w_{\max}. (c) Which class does the CAPE fall into? (d) On a skew-T diagram, what shape has the positive area, and why is its area what was computed?

Solution

Solution of Exercise 21.8.

(a) CAPE=287.0×4×ln⁡(700/300)=287.0×4×0.847=973 J/kg\text{CAPE} = 287.0 \times 4 \times \ln(700/300) = 287.0 \times 4 \times 0.847 = 973\,\mathrm{J}/\mathrm{kg}. (b) 2×973=44.1 m/s\sqrt{2 \times 973} = 44.1\,\mathrm{m}/\mathrm{s}. (c) Between 1 and 1500 J/kg1500\,\mathrm{J}/\mathrm{kg}: positive, but not large. (d) A band of constant horizontal width between two curves; by Proposition 21.14 its area is proportional to ∫ΔT dln⁡p=4ln⁡(700/300)\int \Delta T\,\mathrm{d}\ln p = 4 \ln(700/300), whatever the slant of the isotherms.

Exercise 21.9 ★★

Take A=2.1 nmA = 2.1\,\mathrm{nm} and ammonium sulfate, κ=0.61\kappa = 0.61. (a) Find the critical diameter and supersaturation of a dry particle of 80 nm80\,\mathrm{nm}. (b) What is the smallest dry particle activated in air at 0.3 %0.3\,\% supersaturation? At 0.1 %0.1\,\%? (c) An updraft raises the supersaturation from 0.1 %0.1\,\% to 0.3 %0.3\,\%: what happens to the number of droplets, and to their size for the same liquid water?

Solution

Solution of Exercise 21.9.

(a) Dc=3×0.61×803/2.1 nm=668 nm=0.668 µmD_c = \sqrt{3 \times 0.61 \times 80^3/2.1}\,\mathrm{nm} = 668\,\mathrm{nm} = 0.668\,\text{µ}\mathrm{m}; sc=2A/(3Dc)=2×2.1/(3×668)=0.00210=0.210 %s_c = 2A/(3D_c) = 2 \times 2.1/(3 \times 668) = 0.00210 = 0.210\,\%. (b) Solving sc=4A3/(27κDd3)s_c = \sqrt{4A^3/(27\kappa D_d^3)} for DdD_d: Dd=(4A3/(27κsc2))1/3D_d = \bigl(4A^3/(27\kappa s_c^2)\bigr)^{1/3}, 63 nm63\,\mathrm{nm} at 0.3 %0.3\,\% and 131 nm131\,\mathrm{nm} at 0.1 %0.1\,\%. (c) Every particle between 63 and 131 nm131\,\mathrm{nm} now activates as well: more droplets, and, for the same liquid water, smaller ones.

Exercise 21.10 ★★★

The Norman surface parcel had a CIN of about −130 J/kg-130\,\mathrm{J}/\mathrm{kg} with the virtual temperature. (a) With what vertical speed would it have to leave the ground to reach its LFC by its own momentum? (b) The inversion top is at 889.1 hPa889.1\,\mathrm{hPa} and 20.5 ∘C20.5\,{}^{\circ}\mathrm{C}. To what temperature must the surface air at 971.1 hPa971.1\,\mathrm{hPa} warm for a parcel lifted dry adiabatically to be no colder than the inversion top? (c) If its dew point stayed at 18.5 ∘C18.5\,{}^{\circ}\mathrm{C}, where would its cloud base be? (d) Name two ways in which the atmosphere removes such a lid without the sun.

Solution

Solution of Exercise 21.10.

(a) w=2×130=16.1 m/sw = \sqrt{2 \times 130} = 16.1\,\mathrm{m}/\mathrm{s}, impossible for air at rest near the ground. (b) The potential temperature of the inversion top is 293.65×(1000/889.1)0.2857=303.7 K293.65 \times (1000/889.1)^{0.2857} = 303.7\,\mathrm{K}; a surface parcel with this θ\theta has T=303.7×(971.1/1000)0.2857=301.2 KT = 303.7 \times (971.1/1000)^{0.2857} = 301.2\,\mathrm{K}, 28.0 ∘C28.0\,{}^{\circ}\mathrm{C}. (c) 125×(28.0−18.5)=1190 m125 \times (28.0 - 18.5) = 1190\,\mathrm{m} above the ground, above the old inversion top. (d) A front, which lifts the air on its slope; the slow ascent of the air ahead of a disturbance aloft, which cools the lid itself.

Exercise 21.11 ★★★

(a) Show from Proposition 21.10 that Γs→Γd\Gamma_s \to \Gamma_d as rs→0r_s \to 0. (b) Show that when Lrs/(RdT)≫1Lr_s/(R_dT) \gg 1, Γs→gRvT/(LRd)\Gamma_s \to gR_vT/(LR_d), and evaluate it at 303 K303\,\mathrm{K}. (c) Compare with Γdew\Gamma_{\text{dew}} of Proposition 21.8 when Td=TT_d = T, and explain why the two limits coincide. (d) Is the limit of (b) reached in the atmosphere? Use the value Γs=3.4 K/km\Gamma_s = 3.4\,\mathrm{K}/\mathrm{km} at 30 ∘C30\,{}^{\circ}\mathrm{C} and 1000 hPa1000\,\mathrm{hPa}.

Solution

Solution of Exercise 21.11.

(a) When rs→0r_s \to 0, the numerator tends to gg and the denominator to cpc_p. (b) Write L2rs/(RvT2)=cp×(Lrs/(RdT))×LRd/(RvTcp)L^2r_s/(R_vT^2) = c_p \times \bigl(Lr_s/(R_dT)\bigr) \times LR_d/(R_vTc_p), with LRd/(RvTcp)=5.1LR_d/(R_vTc_p) = 5.1 at 303 K303\,\mathrm{K}: when Lrs/(RdT)≫1Lr_s/(R_dT) \gg 1, also L2rs/(RvT2)≫cpL^2r_s/(R_vT^2) \gg c_p, and Γs→g (Lrs/(RdT))/(L2rs/(RvT2))=gRvT/(LRd)=9.81×461.5×303/(2.501×106×287.0)=1.91 K/km\Gamma_s \to g\,(Lr_s/(R_dT))/(L^2r_s/(R_vT^2)) = gR_vT/(LR_d) = 9.81 \times 461.5 \times 303/(2.501 \times 10^6 \times 287.0) = 1.91\,\mathrm{K}/\mathrm{km}. (c) With Td=TT_d = T, Γdew=gRvT/(LRd)\Gamma_{\text{dew}} = gR_vT/(LR_d), the same. In that limit the latent heat is so large that the parcel can condense almost none of its vapour without warming far above its saturation point: it keeps its mixing ratio and follows the line of constant mixing ratio, along which the temperature falls at Γdew\Gamma_{\text{dew}}. (d) At 30 ∘C30\,{}^{\circ}\mathrm{C} and 1000 hPa1000\,\mathrm{hPa}, Lrs/(RdT)=0.81Lr_s/(R_dT) = 0.81, not much greater than 1, and L2rs/(RvT2)=4180L^2r_s/(R_vT^2) = 4180, four times cpc_p: Γs=3.4 K/km\Gamma_s = 3.4\,\mathrm{K}/\mathrm{km}, still far from 1.9 K/km1.9\,\mathrm{K}/\mathrm{km}. The limit is not reached.

Exercise 21.12 ★★★

Over supercooled water and over ice the saturation vapour pressures are 2.868 and 2.599 hPa2.599\,\mathrm{hPa} at −10 ∘C-10\,{}^{\circ}\mathrm{C}, 1.257 and 1.032 hPa1.032\,\mathrm{hPa} at −20 ∘C-20\,{}^{\circ}\mathrm{C}. (a) Find the supersaturation over ice of air saturated over water at both temperatures. (b) Compare with the supersaturations at which droplets activate. (c) Explain why a cloud of supercooled droplets that receives a few ice crystals turns to ice. (d) Why is the effect weaker near 0 ∘C0\,{}^{\circ}\mathrm{C}?

Solution

Solution of Exercise 21.12.

(a) 2.868/2.599=1.1042.868/2.599 = 1.104 and 1.257/1.032=1.2181.257/1.032 = 1.218: 10.4 %10.4\,\% and 21.8 %21.8\,\%. (b) Droplets activate at a few tenths of a per cent of supersaturation, tens of times less. (c) The crystals grow fast by deposition; they take vapour from the air, which falls below saturation over the droplets; the droplets evaporate and their vapour goes to the ice (the Bergeron process), until the cloud is glaciated. (d) At 0 ∘C0\,{}^{\circ}\mathrm{C} the two saturation vapour pressures are equal: they meet at the triple point, and the ratio tends to one near it.

Problem 21.1

Weekend problem — Reading the morning sounding

The table gives the 12 UTC sounding of Norman, Oklahoma, on 16 April 2026, at the levels used here: pressure, height above sea level, temperature, dew point, and the temperature of the surface parcel lifted as in Method 21.15, computed by integrating Proposition 21.10 in small steps. The sounding reports the tropopause at 235.3 hPa235.3\,\mathrm{hPa} and 10 954 m10\,954\,\mathrm{m}. Use the data of the exercises and the temperatures (not the virtual temperatures) unless told otherwise.

pp (hPa\mathrm{hPa})zz (m\mathrm{m})TT (∘C{}^{\circ}\mathrm{C})TdT_d (∘C{}^{\circ}\mathrm{C})TparcelT_{\text{parcel}} (∘C{}^{\circ}\mathrm{C})
971.134519.018.519.0
925.076317.317.316.9
907.792416.416.416.2
889.1110220.59.515.4
850.0148918.6−1.8-1.813.7
700.031104.5−12.2-12.26.1
500.05725−18.4-18.4−36.6-36.6−8.9-8.9
400.07357−30.4-30.4−55.6-55.6−20.4-20.4
300.09349−42.7-42.7−66.3-66.3−37.1-37.1
250.010560−49.8-49.8−69.8-69.8−48.3-48.3
200.012006−52.3-52.3−78.3-78.3−61.8-61.8

Part I — The air near the ground.

  1. Find the vapour pressure of the surface air, to one decimal.
  2. Find its mixing ratio, to one decimal.
  3. Find its specific humidity.
  4. Find its relative humidity.
  5. By how much does its virtual temperature exceed its temperature?
  6. Find the height of its LCL above the ground with the 125 m125\,\mathrm{m} rule, and the pressure there (take a mean virtual temperature of 294.7 K294.7\,\mathrm{K} for the layer).
  7. Find its potential temperature.

Part II — The lid and the layers above.

  1. At 925 hPa925\,\mathrm{hPa}, is the parcel warmer or colder than its surroundings, and by how much?
  2. Find the lapse rate between 907.7 and 889.1 hPa889.1\,\mathrm{hPa}. What kind of layer is this?
  3. By how much is the parcel colder than its surroundings at 889.1 hPa889.1\,\mathrm{hPa}?
  4. Find the relative humidity at 850 hPa850\,\mathrm{hPa}. Where could such air have come from?
  5. Find the mean lapse rate between 850 and 500 hPa500\,\mathrm{hPa} and classify the layer.
  6. Find the lifted index and its class.

Part III — The energy.

  1. Compute the excess ΔT=Tparcel−T\Delta T = T_{\text{parcel}} - T at each level of the table.
  2. Find the LFC by interpolating ΔT\Delta T linearly in ln⁡p\ln p, to the nearest hectopascal, and its height by interpolating zz the same way.
  3. Find the EL and its height in the same way. Compare with the tropopause.
  4. Compute the CAPE by the trapezium rule in ln⁡p\ln p, from the LFC to the EL through the levels of the table.
  5. Which class does it fall into?
  6. Compute the CIN in the same way, from the ground to the LFC.
  7. A computation on a grid of 1 hPa1\,\mathrm{hPa} gives a CAPE of 2100 J/kg2100\,\mathrm{J}/\mathrm{kg} with the temperatures and 2230 J/kg2230\,\mathrm{J}/\mathrm{kg} with the virtual temperatures. Why do the three values differ?

Part IV — The storm the air could feed.

  1. What vertical speed does parcel theory give with the CAPE of question 17?
  2. At half that speed throughout, how long would the air take to rise from the LFC to the EL?
  3. Name two effects that make a real updraft slower.
  4. To sum up: what updraft speed does parcel theory predict for a storm fed by this morning’s surface air?
Solution

Solution of Problem 21.1.

1. e=es(291.65 K)=21.5 hPae = e_s(291.65\,\mathrm{K}) = 21.5\,\mathrm{hPa}. 2. r=0.622×21.5/(971.1−21.5)=14.1 g/kgr = 0.622 \times 21.5/(971.1 - 21.5) = 14.1\,\mathrm{g}/\mathrm{kg}. 3. q=14.1/1.0141=13.9 g/kgq = 14.1/1.0141 = 13.9\,\mathrm{g}/\mathrm{kg}. 4. es(292.15 K)=22.2 hPae_s(292.15\,\mathrm{K}) = 22.2\,\mathrm{hPa}; 21.5/22.2=97 %21.5/22.2 = 97\,\%. 5. 0.61×0.0141×292.15=2.5 K0.61 \times 0.0141 \times 292.15 = 2.5\,\mathrm{K}. 6. 125×0.5=63 m125 \times 0.5 = 63\,\mathrm{m}; p=971.1exp⁡[−9.81×62.5/(287.0×294.7)]=964 hPap = 971.1 \exp[-9.81 \times 62.5/(287.0 \times 294.7)] = 964\,\mathrm{hPa}. 7. θ=292.15×(1000/971.1)0.2857=294.6 K\theta = 292.15 \times (1000/971.1)^{0.2857} = 294.6\,\mathrm{K}. 8. 16.9−17.3=−0.4 K16.9 - 17.3 = -0.4\,\mathrm{K}: colder. 9. (20.5−16.4)/(1102−924)=4.1/178(20.5 - 16.4)/(1102 - 924) = 4.1/178: the temperature rises by 23 K/km23\,\mathrm{K}/\mathrm{km}, a lapse rate of −23 K/km-23\,\mathrm{K}/\mathrm{km}. An inversion: absolutely stable, a lid. 10. 15.4−20.5=−5.1 K15.4 - 20.5 = -5.1\,\mathrm{K}. 11. es(−1.8 ∘C)=5.4 hPae_s(-1.8\,{}^{\circ}\mathrm{C}) = 5.4\,\mathrm{hPa} and es(18.6 ∘C)=21.7 hPae_s(18.6\,{}^{\circ}\mathrm{C}) = 21.7\,\mathrm{hPa}: 25 %25\,\%. Warm and dry, it can be air that sank and warmed by compression, or air heated over dry ground upwind; it is not the moist air below. 12. 18.6+18.4=37.0 K18.6 + 18.4 = 37.0\,\mathrm{K} over 5725−1489=4236 m5725 - 1489 = 4236\,\mathrm{m}, 8.73 K/km8.73\,\mathrm{K}/\mathrm{km}: between Γs\Gamma_s and Γd\Gamma_d, conditionally unstable. 13. −18.4−(−8.9)=−9.5 K-18.4 - (-8.9) = -9.5\,\mathrm{K}: extreme instability (−8-8 to −10-10). 14. From the ground up: 0, −0.4-0.4, −0.2-0.2, −5.1-5.1, −4.9-4.9, +1.6+1.6, +9.5+9.5, +10.0+10.0, +5.6+5.6, +1.5+1.5, −9.5-9.5 K\mathrm{K}. 15. Between 850 (−4.9-4.9) and 700 hPa700\,\mathrm{hPa} (+1.6+1.6), at the fraction 4.9/6.5=0.7544.9/6.5 = 0.754: ln⁡p=ln⁡850+0.754ln⁡(700/850)\ln p = \ln 850 + 0.754 \ln(700/850), p=734 hPap = 734\,\mathrm{hPa}; z=1489+0.754×1621=2711 mz = 1489 + 0.754 \times 1621 = 2711\,\mathrm{m}, 2.37 km2.37\,\mathrm{km} above the ground. 16. Between 250 (+1.5+1.5) and 200 hPa200\,\mathrm{hPa} (−9.5-9.5), at the fraction 1.5/11.0=0.1361.5/11.0 = 0.136: p=242.5 hPap = 242.5\,\mathrm{hPa} and z=10560+0.136×1446=10 757 mz = 10560 + 0.136 \times 1446 = 10\,757\,\mathrm{m}, about 200 m200\,\mathrm{m} below the reported tropopause (10 954 m10\,954\,\mathrm{m}). 17. The segments (mean ΔT\Delta T times ln⁡\ln of the pressure ratio): 0.80×0.0474=0.0380.80 \times 0.0474 = 0.038; 5.55×0.3365=1.8675.55 \times 0.3365 = 1.867; 9.75×0.2231=2.1769.75 \times 0.2231 = 2.176; 7.80×0.2877=2.2447.80 \times 0.2877 = 2.244; 3.55×0.1823=0.6473.55 \times 0.1823 = 0.647; 0.75×0.0305=0.0230.75 \times 0.0305 = 0.023; sum 6.995 K6.995\,\mathrm{K}. CAPE=287.0×6.995=2010 J/kg\text{CAPE} = 287.0 \times 6.995 = 2010\,\mathrm{J}/\mathrm{kg}, to three figures. 18. Between 1500 and 2500 J/kg2500\,\mathrm{J}/\mathrm{kg}: large. 19. −0.2×0.0486−0.3×0.0189−2.65×0.0207−5.0×0.0450−2.45×0.1467=−0.655-0.2 \times 0.0486 - 0.3 \times 0.0189 - 2.65 \times 0.0207 - 5.0 \times 0.0450 - 2.45 \times 0.1467 = -0.655 K\mathrm{K}; CIN=287.0×(−0.655)=−188 J/kg\text{CIN} = 287.0 \times (-0.655) = -188\,\mathrm{J}/\mathrm{kg}. 20. The table keeps only eleven levels: between them the trapezium rule replaces the curves by straight segments in ln⁡p\ln p and misses their bends. The 1 hPa1\,\mathrm{hPa} grid follows the curves. The virtual temperature adds the buoyancy of the parcel’s water vapour, more abundant in the saturated parcel than in the dry air around it: a larger CAPE, a smaller CIN. 21. 2×2010=63 m/s\sqrt{2 \times 2010} = 63\,\mathrm{m}/\mathrm{s}. 22. 10757−2711=8046 m10757 - 2711 = 8046\,\mathrm{m} at 63.4/2=31.7 m/s63.4/2 = 31.7\,\mathrm{m}/\mathrm{s}: 254 s254\,\mathrm{s}, about four minutes. 23. The weight of the condensed water the updraft carries; the mixing in of drier air around it (entrainment), which cools it as droplets evaporate. 24. Parcel theory gives the updraft of a storm fed by this air a speed of about 63 m s−163\,\mathrm{m}\,\mathrm{s}^{-1}, an upper bound.

Terms defined in this chapter

See all 957 terms in the glossary