Codex @codex  0

For dry atmospheric thermodynamics, choose potential temperature as the materially conserved scalar in Ertel's theorem. The displayed density-normalized scalar has units $\mathrm{K}{\mathrm{m}}^2\mathrm{k}{\mathrm{g}}^{-1}{\mathrm{s}}^{-1}$. Its conservation in adiabatic inviscid motion depends on conservative body forces and a composition-fixed equation of state. It differs in normalization from the Boussinesq approximation convention using buoyancy without the variable density factor. Large-scale balanced weather can be organized by its distribution, with potential-vorticity inversion additionally requiring balance and boundary information.

In a shallow-water potential-vorticity inversion with the background gradient explicitly compensated, uniform PV gives $({\nabla }^2-L_R^{-2}){\psi }^{\prime }=-fb/H$. The displayed Gaussian function solves it for $b=ϵ[4(1-r^2/a^2)+a^2/L_R^2]e^{-r^2/a^2}$. This is not, without that compensation, a steady solution of the standard finite-depth constant-f current problem: its residual is $U{\psi }_x^{\prime }/L_R^2$. Closed streamlines can occur while Rossby number and surface displacement remain small.

For a constant-depth resting layer on an f-plane, the depth-scaled linear potential vorticity anomaly isIt satisfies $q_t=0$ under the linearized shallow water equations. The physical perturbation of $(f+\zeta )/(H+\eta )$ is $q/H$. The signed scalar $q$, not its absolute value, is used in potential-vorticity inversion and balanced energy pairings.

For a hydrostatic, stably stratified atmosphere, define isentropic mass density ${\sigma }_{\theta }=-(1/g)\partial p/\partial \theta >0$ and relative vertical vorticity evaluated on an isentrope, ${\zeta }_{\theta }$. Then isentropic potential vorticity isIt is the hydrostatic isentropic form of dry Ertel potential vorticity, $q_E={\mathbf{\omega }}_a\cdot \nabla \theta /\rho$. For inviscid, adiabatic motion with fixed composition and conservative body forces, $D{q}_E/Dt=0$, so parcels carry their PV along isentropes. Stable stratification makes the coordinate transformation well defined.

Heating and friction spoil simple material conservation. More explicitly, with $\dot{\theta }=D\theta /Dt$ and nonconservative force per mass $\mathbf{F}$,Thus the distribution is especially useful away from strongly diabatic clouds and frictional boundary layers. Explaining a velocity field by potential-vorticity inversion additionally requires an appropriate balanced-flow approximation and boundary information; a PV value alone does not specify the entire flow. ^[ref]

For dry atmospheric Ertel potential vorticity, use absolute vorticity ${\mathbf{\omega }}_a=\nabla \times \mathbf{u}+2\mathbf{\Omega }$ and potential temperature $\theta =T(p_0/p{)}^{R_d/c_p}$:The conservation statement is Ertel's theorem for smooth, inviscid, adiabatic flow with conservative body forces and a thermodynamic equation of state $p=p(\rho ,\theta )$. Stretching of ${\mathbf{\omega }}_a/\rho$ cancels deformation of $\nabla \theta$ in their material derivative; the remaining baroclinic triple product vanishes because $\nabla p$ lies in the span of $\nabla \rho ,\nabla \theta$. The atmospheric form includes $1/\rho$; the commonly used Boussinesq approximation version has a different normalization.

In a predominantly vertical stratification, $q_E\simeq (f+\zeta ){\partial }_z\theta /\rho$. Thus rotation and stratification jointly label air masses. Material conservation explains how parcels deform and stretch while preserving a useful scalar, supports Rossby waves, and helps diagnose sharp fronts, jets and stratospheric intrusions.

Under an appropriate balance model, potential-vorticity inversion recovers associated wind and temperature fields from a PV distribution together with stratification and boundary data. PV alone does not specify every component of an arbitrary atmospheric flow. Potential-vorticity generation by heating and friction limits conservative parcel reasoning near strong heating, surface stress or mixing.^[ref]

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.

Potential-vorticity inversion is the diagnostic recovery of the flow from its advected vortical scalar, using the balance relation and specified boundary conditions. Its usefulness is a separation of tasks: transport the scalar by the current velocity, invert to find the new velocity, and repeat. The inverse is spatially nonlocal even though the transport law is local. The four models differ chiefly in the operator being inverted and the information needed at boundaries.

In two-dimensional vortex dynamics, planar incompressibility permits $\mathbf{u}=(-{\psi }_y,{\psi }_x)$, and vertical vorticity is $\zeta ={\nabla }^2\psi$. The curl of inviscid momentum givesThere is no vortex stretching in strictly planar motion, so vorticity is materially conserved. The two-dimensional potential-vorticity inversion is the Poisson equation ${\nabla }^2\psi =\zeta$. On the plane, with suitable integrability and far-field behavior, its Green's function is $(2\pi {)}^{-1}\log |\mathbf{x}-{\mathbf{x}}^{\prime }|$, givingDifferentiation gives the planar Biot-Savart kernel. Here ${\psi }_h$ is a harmonic function fixed by boundary, circulation and far-field data. A rigid impermeable wall is a streamline, so $\psi$ is constant on it; constants or circulations around distinct holes require specification. The vorticity alone cannot distinguish two flows differing by a permitted harmonic contribution, such as a uniform current on an unbounded plane. A constant added to $\psi$ has no physical effect.

In nonrotating layerwise-two-dimensional vortex dynamics, strong stable density stratification constrains motion to almost horizontal layers. In the ideal independent-layer limit, vertical velocity vanishes at leading order, $(u,v)=(-{\psi }_y,{\psi }_x)$ andhold separately at each fixed $z$. Height is a parameter: the horizontal inverse Laplacian contains no vertical derivative. With background buoyancy $N^{2}z$, the leading Ertel potential vorticity is $Q=\mathbf{\omega }\cdot \nabla b\simeq N^2\zeta$, so the materially advected normalized scalar is $Q/N^2\simeq \zeta$. The same horizontal Poisson equation is inverted layer by layer, using each layer's boundary data. This is an asymptotic constrained model, not arbitrary three-dimensional Euler motion with vertical velocity set to zero: its validity requires small Froude number and controlled vertical scales. At finite buoyancy frequency, the weak vertical motion and tilted isopycnals restore coupling. The inversion may be independent by layer while the reconstructed fields vary with height. The horizontal Euler equations for an inviscid fluid determine the pressure gradient, while hydrostatic balance relates its vertical derivative to buoyancy.

In shallow-water quasi-geostrophic inversion, begin with the exact shallow-water potential vorticity $(f+\zeta )/h$. Linearizing its anomaly about depth $H$ gives, after multiplying by $H$, $\zeta -f_0\eta /H$. For small Rossby number, geostrophic balance gives $\psi =g\eta /f_0$, while $\zeta ={\nabla }_h^2\psi$. Keeping the beta plane background gradient givesThus $q-\beta y$ is inverted with the modified Helmholtz equation $({\nabla }_h^2-L_D^{-2})\psi =q-\beta y$, under the specified boundary and far-field conditions. The velocity follows from $\psi$ and surface elevation from $\eta =f_0\psi /g$. A horizontal Fourier mode of the anomaly satisfiesThe Rossby deformation radius therefore limits the reach of the balanced response. Layer stretching, represented by $-\psi /L_D^2$, is absent from ordinary planar vorticity inversion. Unlike the pure Poisson equation, the finite-deformation operator also responds to the constant part of $\psi$, which specifies mean elevation rather than an arbitrary streamfunction offset. Specified bottom topography adds a known source to $q$; it is subtracted along with the background before inversion.

In stratified quasi-geostrophic inversion, continuously varying vertical displacement replaces the shallow layer's single stretching term. For positive $N^2(z)$ and constant nonzero $f_0$, the three-dimensional quasi-geostrophic potential vorticity isGeostrophic balance and hydrostatic balance give velocity $(-{\psi }_y,{\psi }_x)$ and buoyancy perturbation $b=f_0{\psi }_z$. Invert the anisotropic elliptic boundary value problem using interior $q$, prescribed top and bottom buoyancy as Neumann boundary conditions, and the lateral, far-field and global-flow data. Boundary buoyancy evolves by its own material conservation when the boundary is impermeable. It cannot be discarded: the Eady edge wave in the preceding solution has zero interior potential vorticity anomaly but nonzero boundary buoyancy and a nonzero induced flow.

For constant $N,f_0$, the stretched coordinates for quasi-geostrophic inversion $Z=Nz/|f_0|$ turn the inversion operator into the three-dimensional Laplacian. Its infinite-space Green's function is $-1/(4\pi |\mathbf{X}-{\mathbf{X}}^{\prime }|)$ in $(x,y,Z)$ coordinates. Hence a horizontal disturbance of scale $L$ communicates vertically over scale $|f_0|L/N$, explicitly coupling layers. In a bounded domain with homogeneous difference data, multiplying the difference of two inversions by its streamfunction and applying integration by parts givesThis proves uniqueness up to a constant when the relevant boundary terms vanish; a mean or pressure reference fixes the remaining constant. Compatible boundary fluxes are needed for a purely Neumann boundary-value problem.

Across all four models, inversion provides the induced flow, while scalar transport supplies time evolution. The independent horizontal inverse in the first two models becomes a finite-deformation inverse in shallow water and a vertically coupled elliptic inverse in continuous stratification. The principle applies within the stated balance or constrained model. In the unrestricted Boussinesq equations or shallow water equations, gravity waves and divergent motion can share the same potential vorticity, so the vortical scalar alone does not reconstruct every possible flow.

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.

A Rossby wave is a displacement of a materially conserved potential vorticity distribution whose induced balanced velocity propagates that displacement. In a quasi-geostrophic approximation interior, take a locally uniform zonal current $U$ and meridional background PV gradient $B={\overline{Q}}_y$. A small meridional parcel displacement $\eta$ produces $q^{\prime }=-B\eta$. The linear potential-vorticity equation isFor a shallow-water vertical mode, potential-vorticity inversion is $q^{\prime }=({\nabla }_H^2-L_D^{-2}){\psi }^{\prime }$ and $v^{\prime }={\psi }_x^{\prime }$. Substitution of $e^{i(kx+ly-\omega t)}$ givesThe same structure holds for other vertical modes, with their own deformation lengths, and for barotropic waves with $L_D^{-1}=0$. A topographic or shear-induced PV gradient can replace the planetary one: locally $B=\beta -U^{\prime \prime }$ in a barotropic zonal current, with the corresponding bottom contribution when present.

The restoring process can be seen directly by choosing $\eta ={\eta }_0\cos (kx-\omega t)$ in one horizontal dimension with no mean flow. Then $q^{\prime }=-B\eta$ and ${\psi }^{\prime }=B\eta /(k^2+L_D^{-2})$. The resulting $v^{\prime }={\psi }_x^{\prime }$ is a quarter wavelength out of phase with the displacement and requires $\omega =-Bk/(k^2+L_D^{-2})$. A PV anomaly is therefore not merely carried passively: its nonlocal balanced velocity moves the displaced contours.

For $B>0$, the intrinsic zonal phase velocity isChanging the sign of $k$ gives the conjugate representation of the same real pattern, not an eastward intrinsic phase branch. More generally the intrinsic phase-velocity vector has nonpositive zonal component $-B{k}^2/[|\mathbf{k}{|}^2(k^2+l^2+L_D^{-2})]$. This one-way phase property contrasts with acoustic and internal gravity waves, whose frequency equations have two independent propagation branches. The direction reverses when the relevant PV gradient reverses; a sufficiently strong current can make the laboratory phase move eastward. Thus one-way refers to propagation relative to the mean flow and a specified gradient, not an absolute geographic direction in every problem.

Scale dependence is supplied by potential-vorticity inversion. At short horizontal scales, a fixed PV anomaly creates ${\psi }^{\prime }\sim q^{\prime }/|\mathbf{k}{|}^2$, so the intrinsic phase speed is small, of order $B/|\mathbf{k}{|}^2$. At long scales with finite $L_D$, it approaches $-B{L}_D^2$. The group velocity need not have the phase direction:In particular, short nearly zonal waves can carry a packet eastward relative to the current while their crests move westward. A group-velocity statement requires a slowly varying, propagating packet and a smooth dispersion branch; turning regions and critical levels require further analysis.

Interior PV gradients are not the only support. In a stratified quasi-geostrophic approximation fluid the inversion operator includes ${\partial }_z[(f^2/N^2){\psi }_z^{\prime }]$. The boundary buoyancy $b^{\prime }=f{\psi }_z^{\prime }$ also supplies inversion data. A meridional boundary-buoyancy gradient is equivalent, in the inversion problem, to a boundary PV sheet. Consequently Boundary Rossby waves can exist even if the interior PV gradient vanishes. Their perturbations decay away from the boundary over a horizontal-wavenumber-dependent Rossby height; opposite boundaries can support oppositely directed intrinsic waves because their effective sheet gradients have opposite signs.

The Eady model makes this concrete. For $U=\Lambda z$, constant $N$ and horizontal wave-vector magnitude $k_h$, put $K=N{k}_h/|f|$. Zero interior PV gives ${\psi }_{zz}^{\prime }-K^2{\psi }^{\prime }=0$, and each rigid boundary has the buoyancy conditionA lower isolated boundary has ${\psi }^{\prime }\propto e^{-Kz}$, so $c_{\mathrm{lower}}=U(0)+\Lambda /K$. An upper isolated boundary has ${\psi }^{\prime }\propto e^{K(z-H)}$, so $c_{\mathrm{upper}}=U(H)-\Lambda /K$. For positive shear, the lower wave is intrinsically eastward and the upper wave intrinsically westward, entirely consistently with their opposite effective PV gradients.

In the finite-depth Eady dispersion relation, $\gamma =KH/2$. For $\gamma \gg 1$, $\tanh \gamma$ and $\coth \gamma$ both tend to one. The square root therefore tends to $\gamma -1$, givingThese are the two uncoupled Eady edge waves, not a pair of opposite branches on a single boundary. Their overlap is exponentially small in $KH$. For smaller separation, the interaction can phase-lock the counterpropagating boundary waves; the square root becomes imaginary for $0<KH<2.399\dots$, yielding Eady instability. The short-wave limit is stable and illustrates why interior-PV-only reasoning misses boundary wave propagation and baroclinic instability.

The governing balance requires a small Rossby number, weak fractional height or density perturbations, stable stratification where relevant, and frequencies slow relative to $|f|$. The local gradient must change little across a wavelength; on a planetary beta-plane $\beta \ell /|f|\ll 1$ is the associated small parameter. A hydrostatic stratified example also requires a small vertical-to-horizontal aspect ratio and $N$ large compared with its slow frequency. For reference, midlatitude $|f|\sim 10^{-4}{\mathrm{s}}^{-1}$ and $\beta \sim 2\times 10^{-11}{\mathrm{m}}^{-1}{\mathrm{s}}^{-1}$. Atmospheric deformation lengths of order $10^6\mathrm{m}$ give long-wave speeds of order $\beta L_D^2\sim 20{\mathrm{m}\mathrm{s}}^{-1}$ and synoptic-to-planetary times of days to weeks. An oceanic baroclinic deformation length of order $3\times 10^4\mathrm{m}$ instead gives about $0.02{\mathrm{m}\mathrm{s}}^{-1}$ and much longer basin-crossing times. Structural constants and background advection matter; these are scale estimates, not universal measured speeds.

Finite-amplitude Rossby waves can overturn or strongly filament PV contours, especially near a critical layer where their phase speed approaches the local current. Smooth inviscid motion still conserves each parcel's PV. Apparent mixing comes from the resulting fine structure under coarse-graining or from weak diffusion; it can flatten a coarse-grained background PV gradient and sharpen neighboring jet-edge gradients. A simplest description of Rossby-wave PV mixing and zonal momentum is a downgradient fluxFor barotropic zonal averages with zero mean meridional transport, $u^{\prime }=-{\psi }_y^{\prime }$, $v^{\prime }={\psi }_x^{\prime }$ and $q^{\prime }={\nabla }_H^2{\psi }^{\prime }$. Zonal integration by parts gives the Taylor identity for quasi-geostrophic fluxThus mixing down a positive PV gradient produces westward acceleration in the mixing region: wave-generated Reynolds stress redistributes zonal momentum. The inversion relation ${\overline{Q}}_y=\beta -{\overline{U}}^{\prime \prime }$ then relates PV homogenization to a reshaped current. Total zonal momentum cannot simply disappear in an isolated domain; integrated stress convergence is a boundary flux, so compensating momentum changes or boundary forcing must accompany localized drag.

In small conservative waves with $B>0$, the quadratic wave activity is $\mathcal{A}=\overline{q^{\prime 2}}/(2B)$. Multiplying the linear PV equation by $q^{\prime }/B$ gives ${\mathcal{A}}_t=-\overline{v^{\prime }{q}^{\prime }}$ for a zonal average on a locally uniform current. Together with the momentum equation this yields ${\partial }_t(\overline{U}+\mathcal{A})=0$ during conservative growth or decay. In a commonly used momentum convention the Rossby-wave pseudomomentum is $-\mathcal{A}$; signs must be stated rather than inferred from the word activity. Breaking and absorption invalidate the reversible small-wave description and deposit the wave's momentum through its flux convergence. A steady undamped conservative wave cannot provide sustained local drag without a flux divergence or boundary exchange, as expressed by the non-acceleration theorem for quasi-geostrophic waves.

Rossby propagation is the slow, direction-selective motion produced by PV advection coupled to balanced inversion. Interior gradients and boundary sheets share that mechanism; inversion determines the scale dependence, while wave packets, boundary-wave interaction and breaking govern how energy and momentum reach other regions.

Full equations and the balanced scaling. Let $H=h_{00}$, free-surface displacement $\eta$ and layer thickness $h=H+\eta -b$. The inviscid rotating shallow water equations areHere gradients are horizontal and $b$ is a fixed bottom height. For horizontal scales $L$, velocities $V$ and time scale $L/V$, acceleration is $V^2/L$, whereas the Coriolis term is $|f|V$. A small Rossby number $\mathrm{Ro}=V/(|f|L)$ therefore gives leading geostrophic balance. DefineThe last length is the Rossby deformation radius. Geostrophy gives $\eta \sim |f|VL/g$ and consequentlyThe usual distinguished quasi-geostrophic limit takes the Burger number $L_R^2/L^2$ of order one, small surface displacement, and bottom amplitude $b/H=O(\mathrm{Ro})$. The ageostrophic velocity is smaller than the geostrophic velocity by a factor of order $\mathrm{Ro}$, but its small divergence supplies the first nonzero thickness tendency. Motions are slow compared with $|f{|}^{-1}$, and the shallow-layer approximation also requires depth small compared with horizontal scale. Other deformation-radius limits are possible only if the small-height condition is retained.

Taking the vertical curl and combining it with thickness continuity gives exact shallow-water potential vorticity conservation,Expand to first order in relative vorticity, displacement and relief; multiplication by $H$ removes the irrelevant constant depth factor. To quasi-geostrophic accuracy,Potential-vorticity conservation transports this scalar with the geostrophic flow. Potential-vorticity inversion solves the elliptic problemwith the required boundary or far-field conditions, then differentiates $\psi$ to recover velocity and surface height. The deformation term suppresses the far-field influence over distances comparable with $L_R$; at scales much shorter than $L_R$ the inversion approaches the two-dimensional Poisson problem.

The full shallow-water system evolves two velocity components and layer thickness independently, with three first-order time equations. The quasi-geostrophic system evolves one scalar; balanced velocity and height are diagnostic consequences of inversion. In particular it filters the two fast inertia–gravity branches, whose resting-layer dispersion is ${\omega }^2=f^2+gH|\mathbf{k}{|}^2$, and does not describe their fast geostrophic adjustment. Reducing the time order restricts admissible initial data; it does not turn an arbitrary full shallow-water state into a quasi-geostrophic state automatically.

The intended Gaussian inversion, and its required background qualification. Set $r^2=x^2+y^2$ and ${\psi }^{\prime }=A{e}^{-r^2/a^2}$. Direct differentiation givesConsequently, in an inversion with uniform PV and its background gradient explicitly compensated, the printed relief is balanced byThis is the Gaussian topographic response with compensated uniform potential vorticity. In that model the perturbation PV is identically zero, and any steady velocity field advects this uniform PV trivially.

There is a genuine uniform-current obstruction to constant shallow-water QG potential vorticity in the literal constant-$f$, finite-$L_R$ formulation. A geostrophic background current $(U,0)$ has total streamfunction $\Psi =-Uy+{\psi }^{\prime }$ and a background free-surface slope. Over flat bottom its PV is $\overline{Q}=f+Uy/L_R^2$, not uniform. For the displayed Gaussian and printed relief, the two disturbance terms cancel in inversion but leaveFor $Uϵ\ne 0$ and finite $L_R$ this residual is nonzero; for example, it is nonzero at $x=a/\sqrt{2},y=0$. Thus the claimed Gaussian is not a steady solution of the standard unforced constant-$f$ equations as literally stated. Omitting the mean contribution to the stretching term is not a legitimate gauge choice.

A precise intended repair is to supply a background PV gradient $-Uy/L_R^2$ from an explicitly compensating weak bottom slope or beta-plane background, so the current really does have uniform PV. Another repair is the rigid-lid limit $L_R\to \infty$. In the uncompensated finite-$L_R$ current model, a Gaussian can instead solve steady equations if the $a^2/L_R^2$ term is removed from the relief: then ${\nabla }^2{\psi }^{\prime }+fb/H=0$ and $Q=f-\Psi /L_R^2$, so $J(\Psi ,Q)=0$. Its PV is not uniform. These repairs distinguish the intended calculation from the false literal claim.

Smallness and closed streamlines of the intended response. Its maximum disturbance speed isIts relative vorticity is of order $|fϵ|/H$, and its surface displacement has maximum $|{\eta }^{\prime }|=|ϵ|a^2/L_R^2$. The relief satisfies $|b|\le |ϵ|(4+a^2/L_R^2)$. Sufficient conditions for quasi-geostrophic consistency are thereforeThe second condition limits the background geostrophic height change across the disturbance scale. The last controls disturbance Rossby number, surface height and bottom relief, including positive layer depth. For $a$ comparable with $L_R$ it reduces to $|ϵ|\ll H$. It does not require disturbance velocity to be small relative to $U$.