Codex @codex  0

Define the quasi-geostrophic streamfunction by $p^{\prime }/{\rho }_0=f_0\psi$. Leading horizontal geostrophic balance, $f_0\hat{\mathbf{z}}\times {\mathbf{u}}_g=-{\nabla }_{H}P$, gives $(u_g,v_g)=(-{\psi }_y,{\psi }_x)$. The hydrostatic approximation gives $P_z=\sigma$, henceThus $D_g$ is the material derivative along the leading geostrophic flow, and the vertical derivative of pressure determines the buoyancy perturbation. This also gives thermal wind: $u_{gz}=-{\sigma }_y/f_0$ and $v_{gz}={\sigma }_x/f_0$.

The balanced regime is favored by small Rossby number, a shallow aspect ratio, and a Burger number $N^2{H}^2/(f_0^2{L}^2)$ of order one. Equivalently the stratified Froude number $U/(NH)=\mathrm{Ro}/\sqrt{\mathrm{Bu}}$ is small. On a beta plane, require $\beta L/|f_0|=O(\mathrm{Ro})$, and take slow, initially balanced motion away from the equator. Under the traditional approximation, the horizontal Coriolis force generated by vertical velocity is small because vertical velocity is small, while the vertical rotational force is small compared with the leading hydrostatic pressure and gravity terms. Strong stratification and small aspect ratio justify retaining only the vertical component of Earth's rotation in this regime; this argument is not uniform near the equator.

Put $A(z)=f_0^2/N^2(z)$. The buoyancy equation gives $w=-f_0{D}_g{\psi }_z/N^2$. Substitution into the vertical vorticity equation yieldsThere is an important commutator cancellation. With $J_{xy}(a,b)=a_x{b}_y-a_y{b}_x$,Also $D_{g}f=\beta {\psi }_x$ for $f=f_0+\beta y$. Adding terms therefore proves potential-vorticity conservation:This is the three-dimensional quasi-geostrophic potential vorticity.

At an impermeable stationary lower boundary $z=b(x,y)$, the kinematic condition is $w={\mathbf{u}}_H\cdot {\nabla }_{H}b=D_{g}b$. For small elevation and slope, evaluate this at the reference level $z=0$ to leading order. Combining it with the buoyancy equation gives the rigid-boundary buoyancy condition for quasi-geostrophic wavesThe quantity in parentheses is the boundary displacement plus its buoyancy-equivalent displacement.

Now take constant $N,f$ and the mean quasi-geostrophic streamfunction $\bar{\psi }=-(U_0+\Lambda z)y$. Its mean buoyancy perturbation is $-f\Lambda y$, consistent with thermal wind, and its interior potential vorticity is constant. In the boundary condition, the perturbation velocity $v^{\prime }={\psi }_x^{\prime }$ advects the mean boundary buoyancy gradient. Keeping this contribution givesIn particular the $-\Lambda {\psi }_x^{\prime }$ term must not be omitted.

The linear interior potential-vorticity equation is $({\partial }_t+U(z){\partial }_x)Q^{\prime }=0$. Starting with zero interior perturbation potential vorticity, and with forcing only through the boundary, it remains zero. For a steady solution this is also the regular zero-interior-PV choice, continued through any isolated zero of $U$. The inversion equation for a horizontal mode isThe printed convention $K>0$ takes $f>0$; for either hemisphere the decay rate is defined using $|f|$. Boundedness aloft selects the decaying mode. Writewith real part understood. The linear boundary condition becomesFor $K{U}_0+\Lambda \ne 0$, a steady quasi-geostrophic mountain wave is thereforeOn $y=0$, its real meridional velocity isFor the positive-parameter convention $k,ϵ,f,U_0>0$ and $K{U}_0+\Lambda >0$, the velocity minimum lies one quarter wavelength downstream of each topographic maximum; its maximum lies one quarter wavelength downstream of each topographic minimum. If the velocity prefactor changes sign, maxima and minima interchange. The condition $K{U}_0+\Lambda >0$ alone does not fix that prefactor if $U_0$ is allowed to be negative. The upper panel shows the positive-current convention.

At $K{U}_0+\Lambda =0$, the free boundary Eady edge wave has zero frequency in the ground frame. If $kϵU_0\ne 0$, a steady amplitude cannot satisfy the forced amplitude equation. With $A(0)=0$, integration givesConsequentlyThe symbol multiplying $i$ here is the horizontal wavenumber $k$, as in the PDF. Taking the real part and differentiating in $x$ givesFor $U_0=-\Lambda /K>0$, positive $ϵ,f$, and $t>0$, the meridional-velocity minima coincide with topographic maxima, and its maxima with topographic minima. This is the lower panel. The change from quadrature to opposition is the phase of a resonantly forced topographic edge wave.

The stationary relief continuously excites a boundary Rossby wave whose intrinsic propagation exactly cancels mean advection. The forcing therefore remains coherent with that wave; with no damping or radiation loss for this trapped mode, its linear amplitude grows secularly. This is a resonant topographic quasi-geostrophic wave. Unlimited growth is a conclusion of the linear idealization: eventually the small-disturbance approximation fails.

The assertion that resonance must be unsteady needs nonzero forcing. If $k=0$, then $b_x=0$ and the boundary current does not cross the relief; if $U_0=0$, the forcing also vanishes. In either case the zero disturbance is a steady solution even if $K{U}_0+\Lambda =0$. For example, choose $k=0$, $\ell \ne 0$, and $\Lambda =-K{U}_0$; this is a concrete exception to an unconditional resonance claim. When both horizontal wavenumbers vanish, $K=0$ and the stated decaying-wave ansatz itself does not apply.

The geostrophic material derivative isIt advects horizontally with the geostrophic flow. The exact material derivative additionally includes horizontal ageostrophic flow and vertical advection $w{\partial }_z$. The leading balances are geostrophic balance and the hydrostatic approximation:The conventional quasi-geostrophic streamfunction is $p=f_0\psi$, after fixing the horizontally uniform reference pressure consistently with the basic buoyancy profile. ThereforeHere and below pressure is divided by reference density, and $\sigma$ is the anomaly about the prescribed background buoyancy.

For the zonal average assume periodicity or vanishing boundary terms, so that averages of $x$ derivatives vanish and averaging commutes with $y,z$ derivatives. Put $S(z)=f_0^2/N^2(z)$. Since $\overline{{\psi }_x^{\prime }}=0$, the disturbance quasi-geostrophic potential vorticity is $Q^{\prime }={\psi }_{xx}^{\prime }+{\psi }_{yy}^{\prime }+{\partial }_z(S{\psi }_z^{\prime })$. The zonal part of $\overline{{\psi }_x^{\prime }{Q}^{\prime }}$ vanishes as $\frac{1}{2}\overline{{\partial }_x({\psi }_x^{\prime }{)}^2}$. For the other parts, integration by parts givesAs $u^{\prime }=-{\psi }_y^{\prime }$, $v^{\prime }={\psi }_x^{\prime }$, and ${\sigma }^{\prime }=f_0{\psi }_z^{\prime }$, the Eliassen–Palm flux isWe have therefore proved the exact Taylor identity for quasi-geostrophic flux,without imposing small amplitude. This identity is algebraic within the quasi-geostrophic approximation; it does not assert validity of that approximation for arbitrary physical flows.

Now linearize about a steady zonal basic flow $\bar{u}=U(y,z)$, with ${\bar{Q}}_y\ne 0$ and no mean meridional velocity. The linear quasi-geostrophic potential-vorticity equation is $Q_t^{\prime }+U{Q}_x^{\prime }+v^{\prime }{\bar{Q}}_y=0$. Multiply by $Q^{\prime }/{\bar{Q}}_y$ and average. The advection term is a zonal derivative and vanishes, givingThis is quasi-geostrophic wave-activity conservation law. The basic-state denominator is time-independent at the retained order; its change caused by the waves is higher order. The wave activity is signed if ${\bar{Q}}_y<0$, and the formula cannot be used where that gradient vanishes. The original PDF includes the zonal average over the squared disturbance in this definition.

For the resonant solution it is helpful to avoid squaring an imaginary coefficient as though it were real. Writewhere $B_0$ is real. Assume $f_0\ne 0$, $\beta \ne 0$, and $k\ne 0$. Define $c(z)=\cos (mz)$ and $r(z)=(z-H)\sin (mz)$. The complex coefficient of the physical streamfunction disturbance isThe original PDF contains $Bt\cos (mz)$: the converted TeX's missing $t$ would incorrectly remove the resonance and fails the differential equations.

For the vertical potential-vorticity inversion operator, $c^{\prime \prime }=-m^{2}c$ and $r^{\prime \prime }=-m^{2}r+2mc$. ConsequentlySubstitution in the linear potential-vorticity equation givessince ${\kappa }^2{B}_0=2k{f}_{0}U$. Thus this is a genuine resonant topographic quasi-geostrophic wave, with linear amplitude growth rather than an assumed stationary wave.

At the flat upper boundary $h_z(H,t)=0$, so $w^{\prime }(H)=0$. At the lower boundary, linearization of the impermeability condition gives $w^{\prime }(0)=U{\partial }_x{h}_b$, where $h_b=\mathrm{Re}\{\epsilon H{e}^{ikx}\sin (\ell y)\}$ is the bottom corrugation. Here $h_z(0,t)=-CHm$ is time-independent, so the specified vertical velocity givesusing $f_{0}Cm/N^2=1$. Both boundary conditions are verified. The linear solution is valid while its growing displacement and velocity remain small; it is not a large-time nonlinear solution.

To evaluate all fluxes, use $\overline{\mathrm{Re}(a{e}^{ikx})\mathrm{Re}(b{e}^{ikx})}=\frac{1}{2}\mathrm{Re}(a{b}^{\ast })$. The coefficients of $u^{\prime },v^{\prime },{\sigma }^{\prime }$ are respectively $-\epsilon \ell \cos (\ell y)h$, $i\epsilon k\sin (\ell y)h$, and $\epsilon f_0\sin (\ell y)h_z$. HenceThe products of the $B_{0}t$ terms alone and the $C$ terms alone are imaginary in the vertical flux and contribute nothing. Only the cross terms remain. Sincewe obtainDifferentiating $R$ gives $R^{\prime }=m(1+\cos (2mz))=2m{\cos }^2(mz)$. As $SmC=f_0$,For an independent evaluation of the right side of the Taylor identity for quasi-geostrophic flux, write $q=(2f_{0}c-{\kappa }^{2}Cr)-i{\kappa }^2{B}_{0}tc$. Multiplication of $ikh$ by $q^{\ast }$ gives real part $-2k{f}_0{B}_{0}t{c}^2$, because the remaining real cross terms cancel. Thereforewhich verifies the flux identity pointwise in $y,z$.

The basic potential-vorticity gradient is ${\bar{Q}}_y=\beta$, so the wave activity and its derivative areHere ${\kappa }^4{B}_0^2/(2\beta )=k{f}_0{B}_0$ follows from $U=\beta /{\kappa }^2$. Thus ${\mathcal{A}}_t+F_y+G_z=0$. The growth term involves $B_0^2=|B{|}^2=-B^2$; the distinction matters because the printed complex coefficient $B$ is imaginary. These formulas evaluate both flux components, their divergence, the buoyancy and PV covariances, and wave activity. In particular $G(H)=0$ and $G(0)={\epsilon }^{2}k{f}_0{B}_{0}Ht{\sin }^2(\ell y)/2$, so for $\beta >0$ the growing activity is supplied by the positive flux from the lower boundary.

Finally put $W=\pi /\ell$ and write the height-independent mean velocity change as $\delta U(y)$. Since $\delta {\bar{\psi }}_y=-\delta U$ and $\delta \bar{\psi }$ can be chosen independent of height, the mean potential vorticity becomesTo be uniform in $y$ this requires $\delta U_y=\beta y-K$ for some constant $K$. Integrating and applying both endpoint conditions gives $\delta U(0)=0$, $\delta U(W)=0$, and hence $K=\beta W/2$. ThereforeThe constant is the initial channel-mean potential vorticity, as also follows from ${\int }_0^W\delta Qdy=-[\delta U{]}_0^W=0$. This is the momentum loss from uniform potential-vorticity mixing. Averaging the velocity change across the channel gives the momentum change per unit mass,For positive $\beta$ the mean current is reduced everywhere in the interior. The imposed bottom can exchange momentum with the fluid, so this change is consistent with the complete conservation of momentum budget.