Lecture: Turbulent Heat Exchange

1. Introduction: The Breakdown of Pure Radiative Equilibrium

In our previous lectures on the Gray Gas model, we derived the analytical temperature profile for an atmosphere in pure Radiative Equilibrium. While mathematically elegant, that solution led to a severe physical paradox at the boundary:

T_g = T_e \left( 1 + \frac{\tau_\infty}{2} \right)^{1/4} \quad > \quad T_{air}(0) = T_e \left[ \frac{1 + \tau_{\infty}}{2} \right]^{1/4}

The ground (\(T_g\)) is strictly hotter than the air immediately above it (\(T_{air}(0)\)). This creates an incredibly steep, superadiabatic lapse rate near the surface.

The "Dancing Test Tube" & The Stove Analogy: As discussed in our fluid dynamics examples, heating a fluid aggressively from below creates buoyancy. It is physically impossible to maintain a system where a surface is vastly hotter than the fluid touching it without inducing motion. A parcel of hot air right at the surface is lighter than the cooler air above it, causing it to accelerate upward. This means a purely radiative atmosphere is dynamically unstable.

To restore stability and build a realistic climate model, we must introduce a new physical mechanism that rapidly transports this excess thermal energy away from the surface: Turbulent Heat Exchange.

2. The Nature of Atmospheric Turbulence and Mass Conservation

In the Planetary Boundary Layer, the flow of air is chaotic and swirling. Molecular diffusion is far too slow to transport significant heat over planetary scales. Instead, friction against the rough surface of the planet combined with intense surface heating breaks the air flow into macroscopic eddies.

The Dust Analogy: Consider how dust from the ground reaches the air. Wind blowing parallel to the surface only moves dust laterally. For dust to reach the upper air, there must be a vertical wind component (\(w\)). However, to satisfy mass conservation, air moving upward must be balanced by air moving downward.

If an equal mass of air moves up and down, why is there a net flux of dust (or heat)? The answer lies in concentration. The air moving up carries the high surface concentration, while the pristine air moving down carries zero dust. A net flux exists as long as the upward-moving parcels carry a different property value than the downward-moving parcels.

3. The Basic Formulation: Kinematic Heat Flux (\(\overline{w'T'}\))

To quantify this transport, we use Reynolds Decomposition. We split atmospheric variables into a time-averaged mean component (overbar) and a turbulent fluctuation (prime):

w = \bar{w} + w' \] \[ T = \bar{T} + T'

Near the surface over flat terrain, the mean vertical velocity is essentially zero (\(\bar{w} \approx 0\)). However, the instantaneous vertical transport of temperature is \(wT\). The time average of this product leaves us with the covariance of the fluctuations:

\text{Kinematic Heat Flux} = \overline{w'T'}

Physical Meaning:

A parcel of air heated by the ground becomes buoyant and rises (\(w' > 0\)). Because it is hotter than average, its temperature fluctuation is positive (\(T' > 0\)). The product \(w'T'\) is positive.
Conversely, to replace that rising air, a cooler parcel must sink (\(w' < 0\)). Because it is cooler than average (\(T' < 0\)). A negative times a negative is positive.

Both rising hot eddies and sinking cold eddies result in a net upward transport of heat.

4. The Bulk Aerodynamic Formula

Global Climate Models (GCMs) cannot resolve individual microscopic eddies. We parameterize \(\overline{w'T'}\) using large-scale, "bulk" variables. The magnitude of vertical mixing depends on the macroscopic gradient and the vigor of mechanical turbulence (driven by the mean horizontal wind speed, \(U\)):

\overline{w'T'} \approx C_H \cdot U \cdot (T_g - T_{air})

Here, \(C_H\) is the Bulk Transfer Coefficient for Heat. It is a dimensionless empirical constant that depends on the aerodynamic roughness of the planet's surface. Over a smooth ocean, \(C_H \approx 1 \times 10^{-3}\). Over a rough forest canopy, it can be up to \(1 \times 10^{-2}\).

5. Sensible and Latent Heat Fluxes

To incorporate this into our planetary energy balance, we must convert the kinematic flux (units of \([K \cdot m/s]\)) into an energetic flux (units of \([W/m^2]\)).

A. Sensible Heat Flux (\(SH\))

We multiply the kinematic flux by the volumetric heat capacity of the air (\(\rho c_p\)):

SH = \rho c_p \overline{w'T'} = \rho c_p C_H U (T_g - T_{air})

This is the physical transfer of dry thermal energy.

B. Latent Heat Flux (\(LH\))

On planets with condensable liquids, eddies transport vapor exactly as they transport heat. Evaporation absorbs the latent heat of vaporization (\(L_v\)). Swapping temperature for specific humidity (\(q\)):

LH = \rho L_v C_E U (q_g - q_{air})

Where \(C_E\) is the moisture exchange coefficient, and \(q_g\) is the saturation humidity at the ground.

6. Closing the Surface Energy Balance

At the surface of a real planet, net radiation (\(R_{net}\)) is balanced by turbulence and ground conduction (\(G\)):

R_{net} = SH + LH + G

Turbulence is so efficient that it completely obliterates the pure radiative surface discontinuity. These fluxes siphon energy out of the ground and dump it into the lower atmosphere, tightly coupling \(T_g\) and \(T_{air}\).

7. Surface Properties, Energy Partitioning, and the Bowen Ratio

To formalize how a surface partitions available energy, meteorologists use the Bowen Ratio (\(B\)):

B = \frac{SH}{LH}

Deserts (Sensible Heat Extreme, \(B \ge 5.0\)): With no moisture (\(q_g \approx 0\)), \(LH \approx 0\). To balance intense solar radiation via \(SH\), the temperature gradient \((T_g - T_{air})\) must be massive. This forces the desert sand to spike dramatically in temperature.
Forests (Latent Heat Extreme, \(B \approx 0.2 \text{ to } 0.5\)): Trees act as biological straws, guaranteeing high \(q_g\). Evaporating water absorbs massive energy, leaving little for \(SH\). Consequently, the required \((T_g - T_{air})\) remains small, keeping forests cool.
Urban Areas (Heat Island Effect, \(B \approx 1.5 \text{ to } 3.0\)): Paving removes moisture availability (lowering \(LH\)) while concrete provides high thermal conductivity (raising \(G\)). The lack of evaporative cooling generates the Urban Heat Island.
Oceans (The Infinite Moisture Buffer, \(B \approx 0.1\)): The surface is always at complete saturation. \(LH\) heavily dominates the turbulent exchange.

Quantitative Example: Forcing a Humidity Change Over the Ocean

Because the ocean provides infinite water, the boundary layer is heavily moisture-laden. Let's estimate the continuous energy required to alter the marine relative humidity (\(RH\)) by just 10%.

Using our Latent Heat Flux formula, where specific humidity \(q_{air} = RH \cdot q_g\), the required change in flux is: \[ \Delta LH = \rho L_v C_E U (\Delta RH \cdot q_g) \] Assuming standard marine conditions (\(\rho = 1.2 \text{ kg/m}^3\), \(L_v = 2.5 \times 10^6 \text{ J/kg}\), \(C_E = 10^{-3}\), wind speed \(U = 5 \text{ m/s}\), and saturation humidity \(q_g \approx 0.01 \text{ kg/kg}\)): \[ \Delta LH = (1.2) \times (2.5 \times 10^6) \times (10^{-3}) \times (5) \times (0.1 \times 0.01) \] \[ \Delta LH = 1.2 \times 2500 \times 5 \times 0.001 = \mathbf{15 \text{ W/m}^2} \]
To put this \(15 \text{ W/m}^2\) into perspective, the global radiative forcing caused by doubling atmospheric \(CO_2\) is only about \(3.7 \text{ to } 4.0 \text{ W/m}^2\). It would take roughly four times the energy of a doubled-\(CO_2\) atmosphere just to continuously shift the marine boundary layer's relative humidity by a mere 10%!