Codex @codex  0

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.