Codex @codex  0

For constant buoyancy frequency, density forcing $r_0{e}^{-\mu z}\cos kx\cos {\omega }_{0}t$ gives two frequencies $\omega =\pm {\omega }_0$. The quasi-geostrophic streamfunction satisfies ${\lambda }_{\omega }^2=(N^2/f_0^2)(k^2+\beta k/\omega )$ and $C_{\omega }=ig\mu r_0/(2{\rho }_0{f}_0\omega )$. With zero streamfunction at the bottom, a positive ${\lambda }_{\omega }^2$ selects the decaying solution $C_{\omega }(e^{-\mu z}-e^{-{\lambda }_{\omega }z})/({\mu }^2-{\lambda }_{\omega }^2)$. For $\beta ,k>0$, the westward branch propagates vertically when ${\omega }_0<\beta /k$; a radiation condition selects positive vertical group velocity, not merely bounded amplitude. Coincident exponential exponents are handled by the finite limit $-C_{\omega }z{e}^{-\mu z}/(2\mu )$.

The quasi-geostrophic approximation isolates slowly evolving, nearly balanced motions in a rotating, stably stratified fluid. Ocean eddies and mid-latitude atmospheric disturbances often have this character: a leading Coriolis force–pressure balance suppresses rapid horizontal acceleration, while hydrostatic balance suppresses rapid vertical acceleration. Small ageostrophic flow nevertheless supplies the vertical motion that changes the leading vorticity. The approximation filters fast inertia-gravity waves while retaining nonlinear transport, vertical shear and Rossby waves.

Use a reference density ${\rho }_0$ and remove the background hydrostatic pressure. Let $\pi =p^{\prime }/{\rho }_0$ and let the upward buoyancy perturbation be $b=-g{\rho }^{\prime }/{\rho }_0$. A vertically varying reference density stratification has $N^2(z)>0$. Under the Boussinesq approximation, the inviscid and adiabatic equations areHere $D/Dt={\partial }_t+{\mathbf{u}}_H\cdot {\nabla }_H+w{\partial }_z$. Density variations are small compared with ${\rho }_0$ but retained in the vertical buoyancy force. This is a useful local, incompressible description; a compressible atmosphere needs the corresponding weighted or pressure-coordinate version.

Let horizontal and vertical lengths be $L,H$, speed be $U$, and take the slow time $L/U$. Take $f_0>0$ for definiteness; use $|f_0|$ in the scale estimates in the opposite hemisphere. The principal small parameter is the Rossby number $\mathrm{Ro}=U/(f_{0}L)\ll 1$ at a latitude where $f_0\ne 0$. For the hydrostatic, strongly stratified regime, $H/L\ll 1$ and $f_0/N\ll 1$. Retaining both horizontal vorticity and vortex stretching gives the Burger numberThe geostrophic balance pressure scale is $f_{0}UL$ and the buoyancy perturbation scale is $B=f_{0}UL/H$. Hence $B/(N^{2}H)=\mathrm{Ro}/\mathrm{Bu}\ll 1$: buoyancy perturbations produce small displacements relative to the background stratification scale. The ageostrophic flow is $O(\mathrm{Ro}U)$ and incompressibility then gives $w=O(\mathrm{Ro}UH/L)$. Vertical acceleration is smaller than hydrostatic terms by order ${\mathrm{Ro}}^2(f_0/N{)}^2$ in this scaling. Weak forcing or diffusion can be added at the slow order, but neither is present here.

Define the quasi-geostrophic streamfunction by $\psi =\pi /f_0$. Leading geostrophic balance and hydrostatic balance implyDifferentiating horizontally and vertically gives the thermal-wind balanceThus vertical shear is determined by horizontal buoyancy gradients. The homogeneous Taylor–Proudman theorem is recovered only when those gradients vanish. Stable stratification allows sloping density surfaces and different geostrophic flow at different levels, without abandoning leading balance.

The Prandtl ratio of scales makes that modification quantitative. Horizontal and vertical terms in balanced potential-vorticity inversion have sizes $\psi /L^2$ and $(f_0^2/N^2)\psi /H^2$. Equality givesFor constant $N$, $Z=(N/|f_0|)z$ makes the inversion operator an ordinary three-dimensional Laplacian. Equivalently, an interior balanced disturbance of horizontal size $L$ naturally penetrates a depth of order $|f_0|L/N$. With $N\gg |f_0|$ this is a thin, vertically sheared structure, rather than the arbitrarily tall homogeneous columns suggested by strict Taylor–Proudman theorem behaviour. This is an aspect-ratio statement, not the diffusivity Prandtl number.

To allow a weak meridional variation of the Coriolis parameter, use the beta plane $f=f_0+\beta y$, with $\beta L/f_0=O(\mathrm{Ro})$, or equivalently $\beta L^2/U=O(1)$. The small correction to the leading balance is dynamically important on the slow time. Let $\xi ={\nabla }_H^2\psi$ be vertical relative vorticity and let $D_g={\partial }_t+{\mathbf{u}}_g\cdot {\nabla }_H$. At the first nontrivial order, the vertical vorticity equation and buoyancy equation reduce toRelative-vorticity stretching and tilting are higher order here; planetary-vorticity stretching $f_0{w}_z$ must remain. Eliminate $w$ to expose a materially conserved balanced scalar. Because ${\mathbf{u}}_{g,z}=(-b_y/f_0,b_x/f_0)$, the commutator term ${\mathbf{u}}_{g,z}\cdot {\nabla }_{H}b$ is exactly zero. For $N=N(z)$ this givesAdding the two equations proves the quasi-geostrophic potential-vorticity equationAn additive constant such as $f_0$ changes no dynamics. The horizontal streamfunction advection bracket can express this as $Q_t+{\psi }_x{Q}_y-{\psi }_y{Q}_x=0$. Although small departures from balance generated this evolution equation, the leading geostrophic flow itself advects $Q$ nonlinearly. The stretching term couples horizontal motion to vertical displacement of density surfaces.

This gives the practical advection-and-inversion description. Start with interior three-dimensional quasi-geostrophic potential vorticity, advect it with ${\mathbf{u}}_g$, and recover $\psi$ fromThe stratified quasi-geostrophic inversion is elliptic for positive $N^2$. It must be supplied with lateral and vertical boundary conditions and any circulation or mean-flow constraints; interior potential vorticity alone is not sufficient. At rigid horizontal lids, $w=0$ makes boundary buoyancy satisfy $D_{g}b=0$, and $b=f_0{\psi }_z$ supplies the vertical derivative data for inversion. Boundary buoyancy perturbations can therefore support flow even when the interior anomaly is zero. Once $\psi$ is found, its horizontal derivatives give velocity, its vertical derivative gives buoyancy, and the thermodynamic equation diagnoses $w$. This balanced evolution excludes the independent fast-wave initial data present in the full equations.

Finally, a gradient in $f$ supplies a restoring mechanism even in an otherwise uniform fluid. Linearizing about rest with constant $N$, the quasi-geostrophic potential-vorticity equation becomesSubstitution of $\psi \propto e^{i(kx+ly+mz-\omega t)}$ gives the Rossby wave dispersionFor $\beta >0$, its zonal phase speed is westward relative to a resting background. A meridionally displaced parcel retains potential vorticity, so it acquires a relative-vorticity anomaly that induces the velocity tending to return the disturbance. Background shear and stratification modify the full potential-vorticity gradient; they can alter propagation or permit instability. On a constant f-plane the planetary contribution vanishes, but gradients of relative vorticity, boundary buoyancy or topography can still support balanced wave motion. Close enough to the equator that $f_0$ vanishes, this mid-latitude quasi-geostrophic approximation and its scaling must be replaced.

The basic quasi-geostrophic streamfunction $\bar{\psi }=-\Lambda yz$ gives $\bar{u}=\Lambda z$, $\bar{v}=0$ and constant quasi-geostrophic potential vorticity $\bar{Q}=f$. DefineBecause the background has no potential-vorticity gradient, the linear potential-vorticity equation is simplyIts characteristics have constant $y,z,x-\Lambda zt$, so the general solution is $q^{\prime }=\mathcal{F}(y,z,x-\Lambda zt)$. In particular, $q^{\prime }=F(y,z)G(x-\Lambda zt)$ is the requested arbitrary separated family; one product is not the most general initial field.

Linearizing the vertical velocity requires retaining advection of ${\bar{\psi }}_z=-\Lambda y$ by $v^{\prime }={\psi }_x^{\prime }$. ThereforeThe impermeable boundary at $z=0$ imposesThis is conservation of boundary buoyancy, not merely a condition that ${\psi }_z^{\prime }$ vanish.

For zero interior potential vorticity, write ${\psi }^{\prime }=\mathrm{Re}\{\varphi (z,t)e^{i(kx+ly)}\}$. With $K_h=\sqrt{k^2+l^2}>0$, stratified quasi-geostrophic inversion givesThe decaying solution isNo time dependence has been assumed to obtain this vertical structure. Applying the boundary condition gives $-\mu \dot{A}-i\Lambda kA=0$, henceThus the Eady edge wave travels eastward with horizontal phase velocity $\omega /k=\Lambda /\mu$ for $k\ne 0$.

The propagation mechanism can be made explicit using boundary parcels. Thermal wind gives the background buoyancy gradient ${\bar{b}}_y=f{\bar{\psi }}_{zy}=-f\Lambda$. A small meridional parcel displacement $d$ produces boundary buoyancy perturbation $b^{\prime }=-d{\bar{b}}_y=f\Lambda d$. Meanwhile the decaying potential-vorticity inversion gives $b^{\prime }=f{\psi }_z^{\prime }=-f\mu {\psi }^{\prime }$ at the boundary. Consequently ${\psi }^{\prime }=-\Lambda d/\mu$ there, and the parcel equation $d_t=v^{\prime }={\psi }_x^{\prime }$ becomesA displaced buoyancy contour therefore induces a velocity field that advances that displacement eastward. No interior potential-vorticity gradient is needed because the boundary buoyancy supplies the invertibility data.

Vortex stretching is what sets the induced flow's vertical reach. In the interior, the zero-potential vorticity condition requiresBoundary isopycnal displacement decreases with height as $e^{-\mu z}$, so fluid columns are stretched or compressed over depth ${\mu }^{-1}=f/(N{K}_h)$. Conservation of potential vorticity makes this stretching balance the relative vorticity. The resulting inversion factor $b^{\prime }=-f\mu {\psi }^{\prime }$ gives the propagation speed: stronger density stratification reduces penetration and speed, while larger $f$ increases them. Zero interior potential vorticity therefore does not imply zero interior flow.

For imposed boundary vertical velocity, retain the same decaying structure and writeThe boundary relation becomes the forced Eady edge wave equationThe integrating factor $e^{i\omega t}$ gives, for ${\omega }_0\ne \omega$,The arbitrary homogeneous coefficient $B$ sets the initial disturbance. With $A(0)=0$, this is $A=C(e^{-i{\omega }_{0}t}-e^{-i\omega t})/[i(\omega -{\omega }_0)]$.

At resonance, ${\omega }_0=\omega$, integration instead yieldsAn initially undisturbed fluid has ${\psi }^{\prime }=Ct{e}^{-\mu z}\cos (kx+ly-\omega t)$. The imposed pumping reinforces the same boundary normal mode every cycle, with no dissipative loss, so the linear amplitude grows secularly. The divergence describes resonance of the undamped linear model; sufficiently large wave amplitude invalidates the assumed small-disturbance approximation.

With the free-surface elevation denoted by $\zeta$, the rotating shallow water equations over a fixed bottom areThe momentum gradient is that of the surface elevation, not depth alone; bottom slope enters through the depth and continuity equation. Taking ${\partial }_x$ of meridional momentum minus ${\partial }_y$ of zonal momentum gives the vorticity equationThe terms quadratic in velocity gradients combine into ${\omega }_r{\nabla }_H\cdot {\mathbf{u}}_H$; the constant Coriolis parameter contributes $f{\nabla }_H\cdot {\mathbf{u}}_H$. Continuity gives $D_{H}h/Dt=-h{\nabla }_H\cdot {\mathbf{u}}_H$. Dividing the two equations proves the exact shallow-water potential vorticity invariant

For the quasi-geostrophic approximation, let $U_s$ and $\ell$ be representative velocity and horizontal length scales. Require the Rossby number $\mathrm{Ro}=U_s/(|f|\ell )\ll 1$, a shallow layer, and slow evolution on time $\ell /U_s$. Put $L_R=\sqrt{g{h}_{00}}/|f|$ and assume $\mathrm{Bu}=(L_R/\ell {)}^2$ is order one. Leading momentum then gives geostrophic balance, and a quasi-geostrophic streamfunction can be chosen asIts surface-height scale satisfies $|\zeta |/h_{00}=O(\mathrm{Ro}/\mathrm{Bu})\ll 1$. Take the bottom-height ratio of the same small order, with gentle slopes, so $\eta =\zeta -b$ obeys $|\eta |/h_{00}\ll 1$. Relative vorticity is order $\mathrm{Ro}f$, and expanding the exact potential vorticity givesBecause the leading geostrophic velocity is horizontally nondivergent, ageostrophic advection of the anomaly is one order smaller. At the first evolving order the shallow-water quasi-geostrophic potential vorticity equation isThe constant reference $f/h_{00}$ has zero gradient and drops out. These assumptions also require all topographically induced velocities to retain a small Rossby number; the small bottom-height condition alone is not a full justification of quasi-geostrophic balance.

The uniform eastward current has ${\psi }_0=-Uy$ up to an additive constant, so ${\zeta }_0=-(fU/g)y$. With $b_0=-(fU/g)y$, surface and bottom slopes are equal and $h=h_{00}$. Its relative vorticity vanishes and $Q_0=0$. This sloping background is used locally where its surface and bottom departures remain within the small-height regime.

Add the compact ridge and set $\psi =-Uy+\varphi (x)$. Then $u=U$, $v={\varphi }^{\prime }(x)$, and $Q$ depends only on $x$. The steady equation becomes $U{Q}_x=0$. Uniform incoming potential vorticity fixes that constant to zero, giving the cosine-ridge shallow-water geostrophic responseLet $a=\pi L/2$, $\lambda =L_R^{-1}$ andThe forcing is even, so the decaying response is even. A cosine particular solution inside the ridge leaves an even homogeneous term; outside, retain only the decaying exponential. Since the forcing is bounded and contains no delta function, both $\varphi$ and ${\varphi }^{\prime }$ are continuous at the edges. Write $\varphi =A\cos (x/L)+D\cosh (x/L_R)$ inside and $\varphi =B{e}^{-(|x|-a)/L_R}$ outside. Matching at $x=a$ givesSolving yields $D=A(L_R/L)e^{-a/L_R}$, henceA globally decaying homogeneous solution that is continuously differentiable is zero, so this matched solution is unique.

A steady streamline is a constant-$\psi$ curve. If its upstream ordinate is $y_{\infty }$, it therefore hasFor $f>0$, $\varphi$ is positive: the flow first turns toward increasing $y$, passes the crest at its largest deflection, then returns to its upstream ordinate. Reversing the sign of $f$ reverses the deflection. For $L=3L_R$, the center deflection is $A[1+e^{-3\pi /2}/3]/U$ and each tail decays on scale $L_R$. The figure shows several streamlines and the ridge support; the deflection extends beyond the ridge.

One compatible realization of the initial condition has $h(x)=h_{00}-\hat{b}(x)$ and $\zeta =b_0(y)$. Its stated zonal velocity then satisfies $u(x)h(x)=U{h}_{00}$, so the initial mass flux is constant and $v=0$. The data given only far upstream do not uniquely determine the initial height everywhere; this is the natural mass-flux-compatible choice. It is not fully balanced over the ridge: $fu(x)$ differs from $-g{\zeta }_y=fU$. The resulting transverse acceleration initiates geostrophic adjustment and launches fast inertia-gravity waves.

Initially the relative vorticity is zero, but exact shallow-water potential vorticity is $f/[h_{00}-\hat{b}(x)]$, not the uniform upstream value. That localized anomaly is materially transported downstream, approximately at $U$ in the leading quasi-geostrophic approximation; it is not annihilated by adjustment. Fresh incoming parcels of uniform PV subsequently supply the ridge region, whose balanced potential-vorticity inversion approaches the steady solution above. The fast waves can carry unbalanced energy away on the inertia-gravity time scale, short compared with the advective time at small Rossby number. In an open domain this can leave the local balanced response; an ideal closed or reflecting domain need not settle monotonically, because undamped waves can return. Radiation, advection of the initial PV anomaly, and any weak damping must be distinguished from irreversible loss of PV.

Let $S(z)=d{\rho }_s/dz$, $N^2=-gS/{\rho }_0>0$, and define the quasi-geostrophic streamfunction by $\psi =\tilde{p}/({\rho }_0{f}_0)$. Use horizontal scale $L$, vertical scale $D$, velocity scale $U$ and advective time $L/U$. The required quasi-geostrophic approximation has small Rossby number $\epsilon =U/(|f_0|L)$, $\beta L/|f_0|=O(\epsilon )$, and Burger number $N^2{D}^2/(f_0^2{L}^2)=O(1)$. Hydrostatic, thin-layer scaling has $D/L\ll 1$; density variations are small compared with ${\rho }_0$, as required by the Boussinesq approximation. The leading horizontal flow is nondivergent, with vertical velocity $w=O(\epsilon UD/L)$. Retain forcing of the same slow order as density advection, $r=O(U\tilde{\rho }/L)$; stronger forcing need not preserve this balanced scaling.

Leading geostrophic balance and the hydrostatic approximation giveWrite $D_g/Dt={\partial }_t-{\psi }_y{\partial }_x+{\psi }_x{\partial }_y$. The density equation becomesTaking the curl of the horizontal momentum equations at the next order and using continuity gives the stretching equationSubstitute the expression for $w$. Vertical differentiation commutes with the geostrophic material derivative acting on ${\psi }_z$ here: the extra Jacobian is $J({\psi }_z,{\psi }_z)=0$. Since $N$ depends only on height, this yields the diabatically forced quasi-geostrophic potential vorticity equationThe minus sign is appropriate to the specified density forcing; a positive upward buoyancy forcing is $R=-gr/{\rho }_0$ and gives $D_{g}q/Dt=f_0{\partial }_z(R/N^2)$.

Now take constant $N$ and linearize about rest. Put $s=f_0^2/N^2$ and $S_0=g{f}_0\mu r_0/({\rho }_0{N}^2)$. The forced equation isFor the harmonic heating response of quasi-geostrophic flow, represent the two traveling Fourier components with frequencies ${\omega }_{+}={\omega }_0$ and ${\omega }_{-}=-{\omega }_0$:Each component of the forcing has amplitude $S_0/2$, soEquivalently its vertical structure obeysThese are the time-periodic forced responses; unspecified initial data may add free Rossby waves. Selecting no incoming waves describes forcing switched on causally and then its persistent harmonic part.

For $\beta >0$, the eastward component always has ${\lambda }_{+}^2>0$. Let ${\lambda }_{+}>0$. The growing exponential would introduce an unbounded velocity and pressure disturbance at infinity, so discard it. The bottom condition ${\varphi }_{+}(0)=0$ then givesIf $\mu ={\lambda }_{+}$, this expression has the regular limit $-C_{+}z{e}^{-\mu z}/(2\mu )$, rather than a physical singularity.

For ${\omega }_0>\beta /k$, the westward component also has positive ${\lambda }_{-}^2$. Choose ${\lambda }_{-}>0$ and use the same decaying expression with minus subscripts, including the same coincident-exponent limit. Thus both components are vertically evanescent above the heating.

For $0<{\omega }_0<\beta /k$, putBoth homogeneous solutions $e^{\pm imz}$ are bounded, so boundedness is insufficient. A free Rossby wave with phase $kx+mz-\omega t$ hasThe radiation condition demands energy propagation away from the heating, upward into $z>0$. For $\beta ,k>0$ that selects positive $m$, hence $e^{imz}$ in the westward component. The resulting solution isIt satisfies ${\varphi }_{-}(0)=0$ and, away from the decaying forcing, ${\varphi }_{-}^{\prime }\sim im{\varphi }_{-}$. Its far-field amplitude does not decay in this nondissipative model: continuous forcing supplies an outgoing wave train. Thus the low-frequency westward component propagates vertically while the eastward component remains evanescent.

At the threshold ${\omega }_0=\beta /k$, the westward vertical wavenumber is zero. Excluding the linearly growing homogeneous solution and retaining the bottom condition gives ${\varphi }_{-}=C_{-}(e^{-\mu z}-1)/{\mu }^2$. This bounded limiting response tends to a constant and has zero vertical group velocity; demanding both bottom zero and decay to zero at infinity would overdetermine this threshold problem. For $\beta \le 0$, the branch labels and propagation inequality must be reconsidered; the comparison stated here assumes the usual positive beta-plane gradient.

The hydrostatic buoyancy anomaly is $b^{\prime }=f_0{\psi }_z$, and its leading adiabatic evolution is $D_g{b}^{\prime }/Dt+N_0^{2}w=0$. The small vertical velocity advects the background stratification; vertical advection of the perturbation buoyancy is higher order. ConsequentlyThe impermeability condition at a stationary lower boundary $z=\alpha y$ is $w=\alpha v$. Using $v_g={\psi }_x$ gives the topographic buoyancy boundary condition for quasi-geostrophic flowFor application at the reference plane, the actual small-height condition is $|\alpha |L/D\ll 1$. Slopes scaled by $D/L$ are of Rossby-number order, so this topographic contribution is retained alongside the small ageostrophic vertical velocity. The dimensional slope should not be compared to the Rossby number without also specifying the aspect ratio.

The basic quasi-geostrophic streamfunction is $\overline{\psi }=-\Lambda zy$. It gives $U=\Lambda z$, basic buoyancy gradient ${\overline{b}}_y=-f_0\Lambda$, and spatially constant interior quasi-geostrophic potential vorticity. Linearizing its conservation law therefore givesThe printed equation omits the shear-advection term. This cannot be omitted for general disturbances: an interior anomaly proportional to $\cos [k(x-\Lambda zt)]$ is materially conserved but not independent of time at a fixed location. The intended zero-interior-PV Eady model modes do have $q^{\prime }=0$ and satisfy both expressions; the ensuing mode calculation is valid in that subspace.

At either boundary, linearizing the geostrophic material derivative includes $v^{\prime }{\overline{\psi }}_{zy}=-\Lambda {\psi }_x^{\prime }$. LetThe two boundary equations areThese are exactly the respective bottom and top slope corrections to the rigid-boundary buoyancy condition for quasi-geostrophic waves.

Assume periodicity horizontally and first consider $q^{\prime }=0$. Multiply its defining equation by ${\psi }_x^{\prime }$ and integrate over the layer. The two horizontal Laplacian terms vanish by integration by parts; the remaining volume term after one vertical integration is an $x$-derivative. ThereforeFor nonzero $A_0,A_D$, the boundary equations giveThe top advection term integrates to zero. Subtracting proves conservation of the boundary pseudomomentum of a sloping Eady layer:The original PDF uses squared vertical derivatives and the upper slope $\gamma$ in the first denominator; these details are corrupted in the converted TeX. For arbitrary nonzero interior PV the same calculation instead gives $d\mathcal{P}/dt=(f_{0}s{)}^{-1}\int {\psi }_x^{\prime }{q}^{\prime }dV$, which need not vanish. For example, at an instant take ${\psi }^{\prime }=\cos kx+(z/D{)}^2\sin kx$, independent of $y$; the zonal average of ${\psi }_x^{\prime }{q}^{\prime }$ is $-ks/D^2$. Thus the boundary-only conservation statement also needs the zero-interior-PV qualification.

For an exponentially growing smooth normal mode, this qualification is automatic: the interior equation is $ik(\Lambda z-c)\hat{q}=0$, and a phase speed with nonzero imaginary part cannot equal the real basic velocity. Therefore $\hat{q}=0$. If $A_0{A}_D<0$, the conserved boundary quadratic form is definite, so a nonzero exponentially growing boundary amplitude is impossible. Instability consequently requiresfor $\Lambda \ne 0$. When one $A$ vanishes, the corresponding boundary derivative is simply advected and must vanish for a growing mode. The other boundary then yields a real single-edge-wave speed; if both vanish, the zero-interior-PV Neumann problem has no nonzero growing mode. These degenerate cases are neutral and should not be handled by dividing by zero.

For equal slopes, put $\zeta =z/D$, $C=c/(\Lambda D)$, $M=N_0|k|D/|f_0|$, and $a=N_0^2\alpha /(f_0\Lambda )-1$. The normal-mode interior equation is ${\Phi }_{\zeta \zeta }-M^2\Phi =0$ for the zero-interior-PV modes. ThusThe coefficient equations areSetting their determinant to zero, without dividing by $C$ or $a$, givesso the Eady model with parallel sloping boundaries hasFor positive $k,f_0$, $M$ is the printed $\mu$, and $a$ is the printed $\tilde{\alpha }$. The expression is even in the sign of that vertical scale. For $a\ge 0$ the radicand is at least $1/4$, so both speeds are real and these modes have no exponential growth. For equal slopes the earlier necessary condition is merely $a^2>0$, which is satisfied in many stable cases and is not sufficient. At $a=0$ the speeds are $0$ and $\Lambda D$; the invariant with nonzero denominators cannot be used, but the determinant still correctly gives neutral modes. No claim excluding transient amplification is needed.

For the Boussinesq approximation, subtract the resting background buoyancy $N^{2}z$ and denote its perturbation by $b$. With constant nonzero $f$ and constant $N>0$, the leading geostrophic flow and hydrostatic balance define the quasi-geostrophic streamfunctionThe leading horizontal flow is nondivergent. At the next order, its relative vertical vorticity $q_g={\nabla }_H^2\psi$ is stretched by the small vertical velocity, while the background buoyancy is advected vertically. Thus the curl and buoyancy equations giveThese follow from $D{q}_a/Dt=-q_a{\nabla }_H\cdot \mathbf{u}$ with $q_a\simeq f$ and ${\nabla }_H\cdot {\mathbf{u}}_a=-w_z$, retaining the geostrophic advection of the first-order anomalies. Vertical advection of those anomalies is of higher quasi-geostrophic order.

Differentiate the buoyancy equation in $z$. It is important not to assume without justification that ${\partial }_z$ commutes with $D_g$. Their commutator on $b$ isby hydrostatic balance and thermal wind. Hence $D_g{b}_z=-N^2{w}_z$. Adding $f/N^2$ times this equation to the vorticity equation eliminates $w_z$, yieldingThis is the quasi-geostrophic potential-vorticity equation on an $f$-plane. The additive $f$ makes $Q$ equal to the first-order normalized Ertel PV of total buoyancy $N^{2}z+b$; using a PV anomaly instead removes that constant and does not change the dynamics.

For stretched coordinates for quasi-geostrophic inversion, putThen $(f^2/N^2){\partial }_z^2={\partial }_Z^2$, and the potential-vorticity relation becomesStable stratification makes this operator elliptic. Uniform $N$ and nonzero $f$ are essential to this simple constant-coefficient stretching.

Use the streamfunction convention $u={\psi }_z$, $w=-{\psi }_x$, and take $v=0$. With $J(A,B)=A_x{B}_z-A_z{B}_x$, the material derivative is $D={\partial }_t-J(\psi ,\cdot )$. This convention differs from the horizontal quasi-geostrophic streamfunction convention used below. The nonlinear Boussinesq equations areThe $y$ component of vorticity is $\zeta =u_z-w_x={\nabla }^2\psi$. Differentiate the first momentum equation in $z$ and subtract the $x$ derivative of the second. The pressure terms cancel; the derivatives of the advecting velocity cancel by $u_x+w_z=0$. Thus $D\zeta =-{\sigma }_x$, orConservation of total buoyancy, with background gradient $N^2(z)$, givesBoth equations retain the nonlinear advection terms.

Write $\psi =\bar{\psi }(z)+\varphi$, with ${\bar{\psi }}_z=U(z)=\bar{u}(z)$, and take the background buoyancy anomaly to be zero. The linear equations for $(\varphi ,\sigma )$ areFor a normal mode $\varphi =\hat{\psi }(z)e^{ik(x-ct)}$, the second equation gives $\hat{\sigma }=N^2\hat{\psi }/(U-c)$. Substitute this into the first to obtain the Taylor–Goldstein equationHere $k\ne 0$ and the calculation is initially away from critical levels $U=c$. The first term is the restoring effect of stable buoyancy on vertically displaced parcels. The second comes from advection of the background vorticity gradient $U^{\prime \prime }$: a displacement creates a vorticity anomaly, whose induced velocity acts back on the displacement. This is the shear-flow counterpart of a Rossby wave restoring mechanism; its sign can also permit instability, so it is not invariably stabilizing.