Past exam of the geophysical fluid dynamics course of the University of Oxford / 2023 / 1

Choose $y$ northward and $z$ upward in the Northern Hemisphere. In a locally linear Boussinesq approximation, write the background buoyancy as $\bar{b}=b_{\ast }+By+N^{2}z$, with horizontal buoyancy gradient $B<0$ and positive squared buoyancy frequency $N^2$. The stable density stratification resists a purely vertical displacement, but the isopycnals slope upward poleward:The geopotential surfaces are horizontal in this local approximation. Their inclination differs from that of the isopycnals, which is the relevant baroclinicity.

Exchange two equal small fluid elements, each of volume $V$, between positions separated by $(\delta y,\delta z)$. The adiabatic parcel exchange has no thermal diffusion, so each parcel retains its buoyancy. The background buoyancy difference is $\Delta \bar{b}=B\delta y+N^2\delta z$. Because the variable-density contribution to gravitational potential energy is $-{\rho }_{0}Vbz$, the change in the pair's gravitational potential energy isFor a poleward-rising exchange, set $\delta z=s\delta y$ with $s>0$. ThenThe trajectory lies between a geopotential surface and an isopycnal. The poleward-rising parcel starts more buoyant than the parcel it replaces, and the equatorward-sinking parcel starts less buoyant: the buoyancy does positive work on both. This releases available potential energy into kinetic energy, allowing eddies to amplify if the full rotating dynamics supports the displacement and dissipation does not remove the gain.

A trajectory steeper than the isopycnals, $s>s_{\rho }$, instead costs gravitational potential energy; an equatorward-rising trajectory also costs it. Motion exactly along a geopotential surface or an isopycnal gives zero change in this ideal exchange. At fixed $\delta y$, the largest release occurs at $s=s_{\rho }/2$, since $N^{2}s(s-s_{\rho })$ is minimized there.

This is an energetic possibility argument for baroclinic instability, not a sufficient normal mode instability criterion. In particular, unrestricted parcel motion must not be assumed in a rotating fluid: conservation of angular momentum can impose an additional restoring effect. Symmetric instability has a separate Ertel potential vorticity criterion; ordinary Eady instability can occur without violating that criterion.

Let ${\mathbf{U}}_h=\Gamma z{\mathbf{u}}_h(x,y)$ denote the actual horizontal velocity field. With constant nonzero Coriolis parameter $f_0$, geostrophic balance and hydrostatic balance, using pressure per reference density $\pi$, giveDifferentiate the first equation in $z$ to obtain thermal-wind balance:Taking the dot product with ${\mathbf{U}}_h$ provesThus the horizontal velocity is tangent to the horizontal buoyancy contours at every height. For a fluid whose thermal equation of state makes horizontal buoyancy contours coincide with horizontal temperature contours, there is no horizontal advection across those contours. In a dry terrestrial atmosphere, use potential temperature as the materially conserved thermal variable. A salinity contribution to buoyancy would require separate treatment before interpreting the result as a statement about temperature.

The usual statement about net meridional heat transport also follows explicitly. Introduce a horizontal streamfunction $\Psi$ such that ${\mathbf{u}}_h=(-{\Psi }_y,{\Psi }_x)$. Integration of thermal-wind balance givesand henceFor a zonally periodic function $\Psi$, or vanishing endpoint contributions, the zonal integral is zero. Likewise, the anomalous eddy buoyancy flux isThe required correlation between meridional velocity and thermal anomaly is absent. A nonzero vertical phase tilt, as in the next part, permits that correlation.

Qualification of the wording. The result does not say that the local heat flux vector $b{\mathbf{U}}_h$ vanishes, or that an arbitrary open-domain integral vanishes without boundary conditions. For a concrete counterexample to that literal reading, take ${\mathbf{u}}_h=(0,1)$ and $b=f_0\Gamma x+N^{2}z$. These satisfy the displayed thermal-wind balance, but $vb=\Gamma z(f_0\Gamma x+N^{2}z)$ is generally nonzero. What the hypotheses alone guarantee is zero transport across horizontal buoyancy contours; the vanishing zonally integrated heat transport uses the stated endpoint condition.

Take eastward $x$, upward $z$, positive wavenumbers $k,K$, and a real instantaneous amplitude $A$. A common growing factor does not affect the phase geometry. Put $C=\cosh (Kz)$ and $S=\sinh (Kz)$, using the hyperbolic cosine and hyperbolic sine. The prescribed streamfunction iswhereFollowing a fixed phase givesTherefore the streamfunction pattern tilts westward with height. The lower boundary leads the upper boundary eastward.

In the conventional quasi-geostrophic streamfunction normalization, hydrostatic balance gives $b^{\prime }=f_0{\psi }_z^{\prime }$. ThereforeFor the phase calculation, multiplication by the constant $f_0$ is immaterial. Define a continuous phase ${\delta }_b=atan2(\alpha C,S)$; this avoids a spurious discontinuity at $z=0$. ThenThe buoyancy perturbation pattern tilts eastward with height, opposite to the streamfunction pattern. At $z=0$, ${\psi }^{\prime }=A\cos (kx)$ while $b^{\prime }=-f_{0}AK\alpha \sin (kx)$, so the two fields are a quarter wavelength apart there.

Both perturbation amplitudes are largest at the boundaries. Indeed,which increase with $|z|$. The physical reason is particularly clear through stratified quasi-geostrophic inversion: the Eady model has no interior background potential-vorticity gradient, and a growing normal mode has zero interior potential vorticity anomaly. The independent thermal anomalies reside at the rigid upper and lower boundaries, where meridional motion advects the basic horizontal buoyancy gradient. Each Eady edge wave induces an interior balanced flow that decays away from its boundary on scale $K^{-1}$. Their interacting, phase-locked combination produces the growing baroclinic instability. Boundary confinement of the thermal sources, rather than a rigid-wall requirement ${\psi }^{\prime }=0$, explains the amplitude pattern.

For a sketch with a quarter-wavelength difference between the boundary streamfunction patterns, chooseAn upper positive streamfunction crest is then ${\lambda }_x/4$ west of its lower counterpart, where ${\lambda }_x=2\pi /k$. The buoyancy perturbation crests slope the other way. The quarter-wavelength boundary assumption is used for the streamfunction; differentiating in $z$ fixes the buoyancy perturbation and does not generally give it the same boundary phase separation.

The connection to poleward heat transport can be checked without relying on the sketch. The meridional geostrophic flow and buoyancy perturbation areAverage over one zonal wavelength, using equal averages $1/2$ for squared sine and cosine and zero for their product:In the Northern Hemisphere this is positive: warm anomalies preferentially move northward and cold anomalies southward. For the same eastward-increasing basic thermal wind in the Southern Hemisphere, $f_0<0$ and the basic buoyancy gradient reverses; the negative northward eddy buoyancy flux is again poleward. For thermal expansion coefficient ${\alpha }_T>0$ and specific heat capacity at constant pressure $c_p$, the corresponding sensible heat flux in a thermal Boussinesq approximation is ${\rho }_0{c}_p\overline{v^{\prime }{T}^{\prime }}={\rho }_0{c}_p\overline{v^{\prime }{b}^{\prime }}/(g{\alpha }_T)$. The growing amplitude makes this flux grow as $A^2$. The tilt supplies the correlation absent in an untilted disturbance and reduces the background available potential energy.

The low-latitude regions of Earth absorb more incoming solar radiation than they emit as outgoing thermal radiation, whereas high latitudes have a net annual radiation deficit. The obliquity of solar incidence and the distribution of daylight create a strong latitude dependence of absorbed solar energy; emitted thermal radiation varies less strongly because terrestrial temperature contrasts are smaller. The terrestrial atmosphere and ocean therefore carry energy from the tropical surplus toward the polar deficits. Otherwise the low latitudes would keep warming and the high latitudes cooling until a different local radiative equilibrium was reached.

For a steady, zonally averaged energy balance, let $R(\varphi )$ be net downward radiation at the top of the terrestrial atmosphere and $H(\varphi )$ the total northward energy transport across latitude $\varphi$. A narrow spherical band and conservation of energy giveThe poleward direction is northward in the Northern Hemisphere and southward in the Southern Hemisphere. Seasonal internal energy storage modifies the instantaneous balance, but not the annual need for poleward redistribution.

Several circulation structures contribute in addition to the transient eddies produced by baroclinic instability of the middle latitudes:

The Brewer-Dobson circulation additionally redistributes energy in the stratosphere, with tropical ascent and extratropical descent, although most total atmospheric heat transport occurs lower down. An Eulerian mean middle-latitude Ferrel cell is thermally indirect and can oppose the poleward transport by eddies; listing every meridional circulation as a direct poleward heat carrier would be incorrect.