Codex @codex  0

Exchange two equal fluid elements of volume $V$ without thermal diffusion in a Boussinesq approximation. Each retains its buoyancy. If the separation is $(\delta y,\delta z)$ and the background buoyancy difference is $\Delta b=B\delta y+N^2\delta z$, the pair's gravitational potential energy changes by ${\rho }_{0}V\delta z\Delta b$. For $B<0$, $N^2>0$ and poleward-rising slope $s=\delta z/\delta y$, this is negative exactly when $0<s<-B/N^2$. The minimum at fixed $\delta y$ is attained at half the isopycnal slope. This energetic test alone is not sufficient for dynamical baroclinic instability: rotation and constraints on parcel motion must also be considered.

Baroclinic instability converts potential energy of a stratified sheared basic flow into growing disturbances. In the Eady model, coupled Boundary Rossby waves provide a useful mechanism: their induced velocities displace boundary buoyancy anomalies, and suitable wavelengths permit phase locking. A single neutral boundary wave does not by itself imply exponential instability.

When both rigid boundaries of a constant-shear Eady model have the same small slope $\alpha$, define $a=N^2\alpha /(f_0\Lambda )-1$, $M=N|k|D/|f_0|$ and $C=c/(\Lambda D)$. Zero-interior-PV modes have vertical structure $A\cosh (Mz/D)+B\sinh (Mz/D)$. The boundary conditions $(z/D-C){\hat{\psi }}_{z/D}+a\hat{\psi }=0$ at both endpoints give the displayed dispersion relation. For $a\ge 0$ its square root is real and there is no exponential baroclinic instability. The same-sign boundary-pseudomomentum criterion is necessary, not sufficient.

For constant shear $U=\Lambda z$ between impermeable horizontal walls at $0,D$, with constant buoyancy frequency and no interior potential-vorticity gradient, define $m=NK/|f_0|$. The boundary conditions yieldThe bracket is negative for $0<mD<2.399\dots$, giving baroclinic instability for nonzero zonal wavenumber. For large separation the boundary waves decouple; the lower-wave speed tends to $\Lambda /m$, matching the dispersion relation of a semi-infinite Eady edge wave.

Let $K^2=k^2+l^2$, $A=K^2+F$, $\mathcal{B}=K^2(K^2+2F)$ and use $e^{i(kx+ly-\sigma t)}$. The uniform upper-layer flow $V$ over a resting lower layer has the dispersion relationIts discriminant is $4F^2{\beta }^2{k}^2+\mathcal{B}{K}^2(K^2-2F)l^2{V}^2$. A negative value gives exponential baroclinic instability. This model includes advection across the interface-induced potential-vorticity gradients in a meridional two-layer current.

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.

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.

Let $f(y)$ denote the signed vertical component of the Coriolis parameter; it is negative in the southern hemisphere. In the steady, linear, inviscid interior, depth integration of horizontal momentum giveswhere $\mathbf{U}={\int }_{-H}^0(u,v)dz$ and $\Pi ={\int }_{-H}^{0}pdz/{\rho }_0$. Horizontal curl eliminates the pressure. Its left side is $f{\nabla }_h\cdot \mathbf{U}+\beta V$. Incompressibility and impermeable flat top and bottom imply ${\nabla }_h\cdot \mathbf{U}=w(-H)-w(0)=0$, givingThis is Sverdrup balance. It assumes steady small-amplitude motion, a Boussinesq approximation, negligible nonlinear advection, negligible Rayleigh and bottom/lateral drag in the interior, a flat fixed-depth domain and no normal top/bottom flux. The applied surface wind stress remains a forcing, even though interior friction is neglected. With cooling absent, the steady buoyancy equation additionally requires $B_z=0$, if that equation is imposed; the curl balance itself does not depend on the cooling coefficient. In the closed, $x$-independent channel below, depth-integrated $V=0$, so a nonzero wind-stress curl cannot be balanced by an inviscid interior alone.

Steady forced channel. Now take $r>0,\alpha >0$. Put $a=\pi /L$, $s=(z+h{)}_{+}=\max (z+h,0)$ and write ${\tau }_0$ for the positive wind-stress amplitude. The prescribed profiles giveThe buoyancy jumps across the sharp forcing interface in this simplified model; pressure and the circulation below are continuous. Let $P=p_y/{\rho }_0$, $D=f^2+r^2$, and $b_y^{\mathrm{top}}=-a{B}_0\sin ay/(\alpha h)$. Hydrostatic balance gives $P=C(y)+b_y^{\mathrm{top}}s$. The horizontal momentum equations and their solution areThe no-normal-flow conditions imply ${\int }_{-H}^{0}vdz=0$: this integral is independent of $y$ by incompressibility and is zero at the meridional walls. It determines the otherwise free depth-uniform pressure gradient,DefineA meridional overturning streamfunction satisfying $v=-{\Psi }_z$, $w={\Psi }_y$ isIt vanishes on all four boundaries. Together with the expressions for $b,P,u$, this solves all the steady equations. The formula permits variable $f(y)$; freezing it at a negative $f_0$ gives the simple schematic $w=aA\cos ayG(z)$. No unspecified latitude dependence is silently needed for the sketch.

Since $A>0$, upper-layer meridional transport is northward and the deep return is southward. At constant $f$, the circulation rises on the southern side and sinks on the northern side. The meridional velocity changes sign at $z=-h+h^2/(2H)$, slightly above the forcing interface, not exactly at $-h$. The required schematic follows from this derived field:

Forcing comparison. The buoyancy and wind contributions have the same vertical shape and, with these signs, the same circulation sense. Their ratio isThe printed scaling $f{\tau }_0/{\rho }_0\simeq B_0$ must be understood in magnitude, since $f<0$ and both forcing amplitudes are positive. Surface confinement supplies the small factor $h/L$, making direct buoyancy forcing negligible in this circulation model.

The missing physical mechanism in the stratified Southern Ocean is eddy-induced overturning. The wind-driven mean Deacon cell tilts isopycnals and builds available potential energy; baroclinic instability produces mesoscale eddies whose buoyancy transport largely opposes the mean overturning. The residual circulation depends on the small net balance and the diabatic buoyancy transformation. Wind energy also enters mean kinetic energy; eddies redistribute momentum through form stress and transfer energy toward dissipative scales and bottom friction. The displayed linear model already balances its wind input by imposed drag, but omits this mean-advection/eddy buoyancy budget, so the necessity of an opposing eddy circulation is a physical extension, not another mathematical constraint on the four given equations.

The potential-vorticity gradients in a meridional two-layer current are important here: besides ${\bar{q}}_{n,y}=\beta$, one has ${\bar{q}}_{1,x}=-FV$ and ${\bar{q}}_{2,x}=FV$. They must be retained when linearizing, even though the basic relative vorticity vanishes.

Let $K^2=k^2+l^2$ and $A=K^2+F$. For a normal mode the two-layer quasi-geostrophic potential vorticity amplitudes areBecause the basic velocities are northward $(V,0)$ in the two layers, their advection frequencies are $\sigma -lV$ and $\sigma$. The uniform forcing has no perturbation, $W^{\prime }=0$. The linear equations are thereforeIn particular, a perturbation's zonal velocity advects the basic interface-induced zonal potential-vorticity gradient. Substitution of the plane wave givesA nonzero disturbance exists exactly when the determinant vanishes. The requested relation for linear stability of a meridional two-layer current isEquivalently, with $\mathcal{B}=K^2(K^2+2F)$,Its discriminant isFor $K>0$, exponential baroclinic instability occurs when this discriminant is negative; otherwise the two frequencies are real. This also displays the stabilizing contribution of the planetary vorticity gradient when $k\ne 0$. As a check, $V=0$ gives the uncoupled barotropic mode and baroclinic mode frequencies $-\beta k/K^2$ and $-\beta k/(K^2+2F)$.

For an exponentially growing normal mode with $k\ne 0$, $\mathrm{Im}(kc)>0$. Such a complex $c$ cannot equal the real basic velocity anywhere, so the interior equation forces $\hat{q}=0$ at every height. Decay then fixes the unique vertical structure $C{e}^{-mz}$, and the lower-boundary condition fixes $c=-S/m$, which is real. This contradiction proves the absence of exponential instability in the semi-infinite Eady model:For $k=0$, there is no advective frequency in the linear interior equation, so an exponentially growing mode would again have $q^{\prime }=0$. The wall buoyancy equation then forces ${\psi }_z^{\prime }(0)=0$, which the nonzero decaying exponential cannot satisfy. Thus these modes do not provide a missing growing branch. This conclusion rules out exponential eigenmode growth, not every possible transient response or singular neutral potential vorticity disturbance.

With a second rigid boundary at $z=D$, the same uniform-shear flow is the Eady model. Both boundaries carry buoyancy perturbations and support Boundary Rossby waves. Their intrinsic propagation directions oppose one another relative to their respective local basic flows. At suitable wavelengths their interacting fields can phase-lock and release basic-state potential energy, producing baroclinic instability. The single-boundary system lacks that second interacting boundary wave.

For a direct comparison, write $\Lambda ={\bar{U}}_z=-S$ and $m=NK/|f_0|$. The vertical solution between the two boundaries is $A\cosh (mz)+B\sinh (mz)$, and their buoyancy conditions areEliminating $A,B$ gives the finite-depth Eady dispersion relationThe bracket is negative for $0<mD<2.399\dots$, so those modes with $k\ne 0$ have a growing and a decaying member. At short wavelengths the two edge waves interact weakly and the frequencies are real. In the semi-infinite limit, the lower-wave root tends to $c=\Lambda /m=-S/m$, recovering the neutral result above.

Let $b$ denote buoyancy in a stably stratified Boussinesq approximation. Geostrophic balance and hydrostatic balance imply the thermal wind relation $f_0{U}_z=-{\bar{b}}_y$. Consequently the specified vertical shear requiresFor northern-hemisphere $f_0>0$ and $\Lambda >0$, the mean potential temperature decreases northward and increases upward: ${\bar{\Theta }}_y=-({\Theta }_0/g)f_0\Lambda$, ${\bar{\Theta }}_z=({\Theta }_0/g)N^2$. The corresponding mean isentropes slope upward toward the north, with $dz/dy=f_0\Lambda /N^2$. Stable potential-temperature stratification does not require the actual temperature to increase upward; the dry adiabatic lapse also enters that comparison.

In the rigid-lid constant-shear baroclinic model, perturbation quasi-geostrophic potential vorticity isThe mean interior potential-vorticity gradient vanishes on this constant-$f_0$, constant-$N$, linear-shear basic state. Linear advection of that perturbation therefore becomes $({\partial }_t+\Lambda z{\partial }_x)q^{\prime }=0$, the supplied interior equation. It retains slow balanced motion and omits friction, planetary-vorticity gradients and nonlinear disturbance advection.

Substitution of ${\psi }^{\prime }=\Psi (z)\cos (ly)e^{ik(x-ct)}$ givesA possible branch has zero interior perturbation potential vorticity, soBoth supplied applications have $f_0>0$, so this equals the source's expression. For an unstable mode the nonzero imaginary part of the phase speed prevents a real critical level $\Lambda z=c$; for real speeds the displayed solution remains a valid zero-interior-PV branch away from, and smoothly through, any critical level. Taking the real part gives the physical disturbance.

The rigid lids enforce $w^{\prime }=0$. Linearized adiabatic buoyancy conservation isUse the geostrophic conventions $v^{\prime }={\psi }_x^{\prime }$ and $b^{\prime }=f_0{\psi }_z^{\prime }$. At a lid this gives $({\partial }_t+\Lambda z{\partial }_x){\psi }_z^{\prime }-\Lambda {\psi }_x^{\prime }=0$, henceThese are the no-normal-flow condition expressed through boundary buoyancy advection, not conditions of zero horizontal velocity. Boundary temperature anomalies can interact even though this solution has no interior perturbation PV, allowing baroclinic instability to extract mean available potential energy.