Codex @codex  0

A useful stipulated dispersive medium obeys the Klein-Gordon equation $u_{tt}-c^2{u}_{xx}+{\omega }_0^{2}u=0$, with $c>0$ and ${\omega }_0>0$. Its positive-angular frequency dispersion relation isHere ${\omega }_0$ is the gap at zero wavenumber, not the central angular frequency ${\omega }_c=\omega (k_0)$ of a wave packet; the parameter in the preceding formula is $A={\omega }_0/c$. For $k>0$, the phase velocity satisfies $v_p>c$, the group velocity satisfies $0<v_g<c$, and $v_p{v}_g=c^2$. In this stipulated medium $c$ is the limiting propagation speed of the model; it need not mean the vacuum speed of light. The same square-root dispersion relation appears in waveguide phase and group velocities, with a waveguide cutoff frequency playing the role of ${\omega }_0$.^[ref]

For the original dimensionless choice $c={\omega }_0=k_0=1$, a carrier wave has phase velocity $v_p=\sqrt{2}$, while a narrow wave packet has leading group velocity $v_g=1/\sqrt{2}$. After time $t=8$, the carrier's zero-wave phase reference has moved to $x=8\sqrt{2}\simeq 11.31$, whereas the leading wave envelope centre is at $x=8/\sqrt{2}\simeq 5.66$. A carrier wave crest therefore moves through the wave envelope; it does not carry an identifiable piece of the envelope with it.

The origin-to-curve secant on the dispersion diagram gives phase velocity, while its local tangent gives group velocity. As $k\to 0^{+}$, $v_p\to \infty$ and $v_g\to 0$; at $k=0$ the phase velocity itself is undefined. As $k\to \infty$, both approach $c$, from opposite sides. On the same positive-angular frequency branch with $k<0$, both signed velocities are negative and their product remains $c^2$. These conventions distinguish direction from phase speed. Reversing both $k$ and $\omega$ represents the same real plane wave; changing just one sign changes its direction.

For an original finite-bandwidth check, take equal-amplitude plane waves at $k_1=0.8$ and $k_2=1.2$, still with $c={\omega }_0=1$. The exact two-tone wave envelope speed iswhile its carrier wave has phase velocity $(\sqrt{1+0.8^2}+\sqrt{1+1.2^2})/2\simeq 1.421337$. Their product is exactly $1$: rationalizing the difference gives $v_{\mathrm{beat}}=2/({\omega }_1+{\omega }_2)$. Nevertheless $v_{\mathrm{beat}}\ne {\omega }^{\prime }(1)$; the derivative describes the limit of vanishing separation, not an arbitrary pair of tones. The two-tone wave envelope repeats indefinitely, so it has no isolated first arrival and cannot measure front velocity.

The long step-trapped topographic Rossby wave has equal phase velocity and group velocity, $c_T=fd{R}_D/H_0$. For an outer surface profile $\eta ={\eta }_0\mathrm{sgn}(x)-A(x,t)e^{-|y|/R_D}$, cross-step matching gives $A_t+c_T{A}_x=2c_T{\eta }_0\delta (x)$. With $A(x,0)=0$, its weak solution is $A={\eta }_0[\mathrm{sgn}(x)-\mathrm{sgn}(x-c_{T}t)]$. This describes the slowly varying outer response; the fast adjustment region at the origin is unresolved.

Use a slit light source, a rectangular transparent block on plain paper, pencil, ruler and protractor. Trace the block outline, keep it fixed, and mark widely separated points on incident and emergent light rays. Remove the block and join the entry and exit points to reconstruct its internal light ray. Draw the surface normal vector at entry; measure $i$ outside and $r$ inside from that normal line. Choose a path exiting the opposite parallel face rather than a side edge. Repeat for a range of $i$ with the same block and source colour; record $i,r,\sin i,\sin r$ and repeat-reading spread. Plot $\sin i$ vertically against $\sin r$ horizontally: Snell's law predicts slope $n_{\mathrm{block}}/n_{\mathrm{air}}$, approximately the block's refractive index. The emergent light ray should be parallel to the incident one. Beam width, pencil marks, protractor resolution, block movement and colour spread limit precision; a narrow-colour source reduces spread from optical dispersion. These are measured angular readings, not a direct timing of phase velocity.^[ref]

A normal incidence trial checks the drawn normal and the absence of bending, but the ratio $\sin i/\sin r$ is then $0/0$. Small nonzero angles also make that ratio sensitive to reading errors. Compare several trials through the slope rather than relying on one nearly axial pair. This rectangular-block measurement complements the existing total internal reflection experiment, which probes a different boundary condition.

Write $\beta =D/{\rho }_w$ and $\alpha =m/{\rho }_w$ in the dispersion relation of a floating elastic plate. A stationary point of squared phase velocity satisfies $\beta k^4(3+2\alpha k)=g(1+2\alpha k)$. If plate inertia is neglected, its solution is $k=(g/(3\beta ){)}^{1/4}$. A steadily translating load can resonate only when its speed matches a phase velocity; this threshold differs from the group-velocity minimum of a flexural-gravity wave.

The zero-phase velocity ansatz $w=(\mathrm{sech}z{)}^k(\tanh z{)}^{1-k}$ solves the Taylor–Goldstein equation away from its critical level of an internal gravity wave exactly when $J=k(1-k)$. Let $y=\tanh z$ and $a=1-k$. Then $w^{\prime }/w=a/y-y$ and $w^{\prime \prime }/w=a(a-1)/y^2-a-1+2y^2$, so the differential-equation residual divided by $w$ is $[J-k(1-k)](y^{-2}-1)$. The curve has maximum $1/4$ at $k=1/2$. The endpoint $k=0$ does not decay at infinity; fractional powers at the critical level require a branch and regularity qualification.

A nondispersive wave has constant phase velocity on a propagation branch. For $\omega =ck$, its group velocity is also $c$. A frequency relation with a nonzero offset, such as $\omega =ck+b$, has constant group velocity but a phase velocity $c+b/k$ depending on wavenumber; constant group velocity alone is therefore not sufficient for nondispersion.

At normal incidence, the incident light ray follows the surface normal vector. Snell's law gives ${\theta }_t=0$ for positive refractive indices, so the transmitted light ray does not bend. Nevertheless, its phase velocity and wavelength change if the refractive index changes. A single axial light ray therefore cannot measure apparent depth: neighboring light rays and their directions supply that information. A reflected component can still be present; no bending of the transmitted light ray does not imply complete transmission.

With no forcing, the linear dispersion relation has positive phase velocity $\sigma k^2/(\beta {\alpha }_1)$ and group velocity three times this. A decaying permanent wave with $\xi =x+ct$ instead obeys $\sigma (A^{\prime }{)}^2=A^2(\beta {\alpha }_{1}c+{\alpha }_{2}A/3)$. For positive coefficients a nonzero homoclinic wave requires $c>0$ and is negative, travelling upstream. Its inverse width is $k_s=\frac{1}{2}\sqrt{\beta {\alpha }_{1}c/\sigma }$. The opposite directions reflect the nonlinear amplitude-dependent propagation speed.

In units of elastohydrodynamic penetration length and inverse forcing frequency, $y_t=-y_{xxxx}$. Prescribing $y(0,t)=y_0\cos t$, zero bending moment $y_{xx}(0,t)=0$, and decay at infinity givesThe two components have phase velocities $-1/S$ and $1/C$. Their unequal spatial attenuation produces a nonreciprocal shape cycle even though the endpoint actuation is reciprocal.

Let $H=h_{00}+\zeta$ and write the horizontal shallow water equations on the constant f-plane asKeeping only terms linear in the disturbance about rest gives the linearized shallow water equationsDefine the horizontal divergence $\delta =u_x+v_y$ and relative vorticity $q=v_x-u_y$. Taking respectively the divergence and vertical curl of the momentum equation yieldsTogether with ${\zeta }_t=-h_{00}\delta$, this provesThe nonlinear shallow-water potential vorticity is $Q=(f+q)/H$. The nonlinear vertical vorticity equation gives $D(f+q)/Dt=-(f+q)\delta$, while mass conservation gives $DH/Dt=-H\delta$; their ratio is therefore materially conserved. Expanding about rest,Advection of its constant background contributes nothing at first order, explaining why the displayed linear anomaly is conserved at each fixed position.

A further time derivative eliminates $q$ and $\zeta$ from the divergence equation:For a plane wave with horizontal wavevector $(k,l)$, write $K=\sqrt{k^2+l^2}$. More directly, substituting the plane wave into the three linear equations gives a coefficient matrix with determinant proportional to $\omega [{\omega }^2-f^2-c_0^2{K}^2]$. Thus the full linear rotating shallow-water dispersion relation has one balanced branch and two wave branches:The latter are inertia-gravity waves. The positive branch starts at $f$ with zero slope and tends to the straight line $c_{0}K$; the negative branch is its reflection. The requested sketch, including the stationary branch, is

For the positive-frequency branch, the magnitudes of the phase velocity and group velocity areThe group velocity vector is $c_0^2(k,l)/\omega$; the negative branch reverses its direction. Near $\omega =f$, namely $K{R}_D\ll 1$ with the Rossby deformation radius $R_D=c_0/f$, $c_p\sim f/K$ is large while $c_g\sim c_0^{2}K/f$ is small. The $K=0$ limit is an inertial oscillation, for which a spatial phase speed is not defined. Far above $f$, both speeds approach $c_0$, the shallow water gravity-wave speed, with $c_p>c_0$ and $c_g<c_0$ at finite $K$.

On the zero-frequency branch, geostrophic balance givesFor a nonconstant plane wave, $\hat{q}=-(g{K}^2/f)\hat{\zeta }$: relative vorticity and elevation have opposite signs. They can both be nonzero and spatially structured, whereas the horizontal divergence vanishes identically. A constant elevation is also a zero-frequency solution, with zero velocity and zero relative vorticity.

For the initially resting ridge, set ${\zeta }_0(x)=ϵh_{00}{e}^{-x^2/L^2}$ and assume an unbounded domain, uniform in $y$. Initial potential vorticity fixesIn the balanced remainder, insert $q_s=(g/f){\zeta }_s^{\prime \prime }$ to obtainThe decaying Green's function of this modified Helmholtz equation gives the Geostrophic adjustment of a Gaussian height ridge explicitly:The initial outward pressure force creates a cross-ridge current; the Coriolis acceleration turns that current and builds the along-ridge geostrophic flow. The remainder is radiated inertia-gravity waves, not a dissipative disappearance of energy.

For precision, use the Fourier transform $\hat{\zeta }(k)=\int \zeta (x)e^{-ikx}dx$. Initial rest implies ${\zeta }_t(x,0)=0$, and the conserved anomaly givesConsequentlyThese formulas satisfy both initial conditions. They also show that the waves carry no linear potential vorticity anomaly: $q=(f/h_{00})(\zeta -{\zeta }_0)$ and $\delta =-{\zeta }_t/h_{00}$ reconstruct the remaining fields, with $u$ and $v$ obtained from their derivatives. On the infinite line, dispersive wave packets propagate away and decay locally, leaving the balanced ridge and its oppositely directed flanking currents. Conservation of energy still holds globally. In a reflecting finite basin, standing waves may persist, so a permanent pointwise relaxation is not guaranteed without dissipation.

If $L\gg R_D$, most of the initial elevation remains balanced: ${\zeta }_s\simeq {\zeta }_0$ with small smoothing corrections. If $L\ll R_D$, the balanced ridge spreads to width $R_D$, and away from the original narrow coreThus the final central height is much smaller than the initial one and the wave component is substantial. The integral of the balanced elevation equals the initial integral, since the Green's function has unit integral; the oscillatory component has zero net added volume.

Let $p$ be kinematic pressure, with the resting hydrostatic part subtracted, and let $\sigma$ be the buoyancy perturbation about background buoyancy $N^{2}z$. The nonrotating inviscid Boussinesq equations, with no buoyancy diffusion, areHere the third equation is mass conservation for incompressible flow, and the last expresses buoyancy conservation following a parcel. Linearization about rest removes the products of disturbances. For a plane internal gravity wave proportional to $e^{i(kx+mz-\omega t)}$, the resulting Linearized Boussinesq equations becomeWith $K=(k^2+m^2{)}^{1/2}$ and $k\ne 0$, the first and third equations give $\hat{p}=-m\omega \hat{w}/k^2$. Substitution in vertical momentum and buoyancy gives the dispersion relationSince the velocity is the time derivative of the displacement field, $\hat{\mathbf{u}}=-i\omega \hat{\mathbf{\xi }}$. Thus the complete buoyancy polarization of a plane internal gravity wave isBoth velocity and displacement are perpendicular to the wavevector; pressure and velocity are in quadrature with buoyancy and displacement.

For upward/rightward phase velocity, choose $k,m>0$ and $\omega =Nk/K$. Differentiating this dispersion relation gives the group velocity, while the vector phase velocity points along the wavevector:Their dot product is zero and ${\mathbf{c}}_p+{\mathbf{c}}_g=(N/K,0)$. Draw ${\mathbf{c}}_g$ from the tip of ${\mathbf{c}}_p$: the right-triangle geometry of internal-wave velocities has horizontal hypotenuse $N/K$. This hypotenuse is distinct from the horizontal phase velocity $\omega /k$.

At an instant with real complex amplitude $\hat{\sigma }=A$, write $\varphi =kx+mz-\omega t$. The spatial fields areThe constant-phase lines of an internal gravity wave slope down to the right, with $dz/dx=-k/m$. Buoyancy maxima lie a quarter wavelength from pressure maxima; the oscillating displacement and velocity run along the phase lines. The group velocity points down and right along them, despite upward/rightward motion of the wave phase.

Let $P=\overline{{\xi }_{x}u+{\zeta }_{x}w}$, using horizontal averaging over a wavelength. This is the displacement pseudomomentum of an internal gravity wave, whose sign convention is opposite to $k\overline{E}/\omega$. For complex coefficients, $\overline{ab}=\frac{1}{2}\mathrm{Re}(\hat{a}{\hat{b}}^{\ast })$. The above wave polarization consequently givesSince $c_{gz}=-\omega m/K^2$, these obey $\overline{{\zeta }_{x}p}=c_{gz}P$. With weak kinematic viscosity, each leading monochromatic velocity component satisfies ${\nabla }^{2}u=-K^{2}u$, ${\nabla }^{2}w=-K^{2}w$. The dissipative term is therefore $-\nu K^{2}P$, so the supplied averaged balance reduces toFor constant $k,m,N$ and a steady wave envelope, $c_{gz}{P}_z=-\nu K^{2}P$. As $P$ is proportional to squared wave amplitude, viscous attenuation of an internal gravity wave givesHere $c_{gz}<0$, so amplitude decreases below a source at $z_s$: the source is above the receiving level. The same conclusion follows directly because downward group velocity carries energy from the source downward. Upward phase velocity does not locate the source.

For a slowly varying mean shear flow $U(z)$, the local intrinsic frequency is $\hat{\omega }=\omega -kU(z)$. A local WKB approximation uses ${\hat{\omega }}^2=N^2{k}^2/(k^2+m(z{)}^2)$, with varying vertical wavenumber and group velocity. The appropriate transported quantity is wave action, or horizontal wave pseudomomentum $kE/\hat{\omega }$; its vertical flux has viscous loss, and the changing polarization factors must be retained. In particular, $P/|\hat{\sigma }{|}^2$ is no longer constant.

Near a critical level of an internal gravity wave, $\hat{\omega }\to 0$, so $|m|\sim N|k|/|\hat{\omega }|$ and $|c_{gz}|\sim {\hat{\omega }}^2/(N|k|)$. The damping per vertical distance, of order $\nu K^2/|c_{gz}|$, grows rapidly. The waves are strongly attenuated over a narrow altitude range and transfer momentum to the mean flow there or before reaching the critical level. The inviscid WKB approximation eventually fails; viscosity and possibly internal-wave breaking regulate the small scales.

For a homogeneous lossless elastic solid, insert the plane wave $u_i=\mathrm{Re}\{U{a}_i{e}^{i(\mathbf{k}\cdot \mathbf{x}-\omega t)}\}$ into the equations of linear elasticity. Put $\mathbf{k}=\kappa \mathbf{n}$, $|\mathbf{n}|=1$, and $v=\omega /\kappa$. The resulting eigenvalue problem isThe minor and major symmetries make this acoustic tensor a real symmetric matrix. For a stable material it is positive-definite, so its three positive eigenvalues give three real positive phase velocities and therefore three slownesses $s_a=1/v_a$. Coincident eigenvalues are possible, so there need not be three distinct speeds. For distinct eigenvalues the displacement polarization are orthogonal; in an anisotropic solid they are generally quasi-longitudinal and quasi-transverse rather than exactly longitudinal and transverse. Stability, or at least strong ellipticity in elasticity, is required for the assertion about positive real speeds.

The elastic slowness surface is the locus of the three vectors $\mathbf{p}=\mathbf{n}/v_a(\mathbf{n})=\mathbf{k}/\omega$. Its equation isFor each wavefront normal there are three radial sheets, which may meet at degeneracies. To obtain the energy flux directly, multiply the elastic equation by ${\dot{u}}_i$. Forthe symmetry of stress gives ${\partial }_{t}W={\partial }_j({\sigma }_{ij}{\dot{u}}_i)$. Thus the instantaneous elastic-wave energy flux is $F_j=-{\sigma }_{ij}{\dot{u}}_i$. With a real branch polarization $\mathbf{a}$, harmonic averaging givesNow differentiate $c_{ijkl}{p}_j{p}_l{a}_k=\rho a_i$ along a regular slowness surface sheet and contract with $a_i$. The terms containing $d\mathbf{a}$ cancel by symmetry and the original eigenvalue equation, leavingEvery tangent $d\mathbf{p}$ is consequently perpendicular to $⟨\mathbf{F}⟩$: the energy flux is normal to the slowness surface. At a degenerate point this statement is made on the limiting branches, rather than assuming a unique normal.

For completeness, equality of the mean kinetic and strain energies gives $⟨W⟩=\rho {\omega }^2|U{|}^2|\mathbf{a}{|}^2/2$. Differentiating the acoustic tensor equation in the wavevector variables similarly givesThus the group velocity is the elastic energy velocity. The homogeneous elastic dispersion relation is degree one in $\mathbf{k}$, so differentiating $\omega (\lambda \mathbf{k})=\lambda \omega (\mathbf{k})$ at $\lambda =1$ gives $\mathbf{p}\cdot \mathbf{g}=1$.

The elastic wave surface is the unit-time locus of these group velocities. From an impulsive point source, the propagating wavefronts at time $t$ are its sheets scaled to $\mathbf{x}=t\mathbf{g}$. A particular source may fail to excite a mode in some directions, and its amplitudes require the source polarization and geometrical spreading; the wave surface determines the arrival geometry. It does not specify all near-field motion between the radiative fronts. Geometrically it is the envelope of the planes $\mathbf{p}\cdot \mathbf{x}=1$: tangent variation gives $d\mathbf{p}\cdot \mathbf{x}=0$, so $\mathbf{x}$ is the appropriately normalized normal to the slowness surface.

In a symmetry-plane cross-section write $\mathbf{p}=r(\theta )\mathbf{n}$, $\mathbf{n}=(\cos \theta ,\sin \theta )$ and $\mathbf{t}=(-\sin \theta ,\cos \theta )$. The two envelope conditions giveThe numerator is the signed curvature numerator of the slowness surface cross-section. At a simple inflection point, $h=v+v^{\prime \prime }=0$ and $h^{\prime }\ne 0$. Then ${\mathbf{g}}^{\prime \prime }=h^{\prime }\mathbf{t}$ and ${\mathbf{g}}^{\prime \prime \prime }=h^{\prime \prime }\mathbf{t}+2h^{\prime }{\mathbf{t}}^{\prime }$ are linearly independent. These are precisely the local conditions for an ordinary cusp. The slowness curvature and wave-surface cusps relation explains why a smooth nonconvex slowness surface can generate a cusped wavefront and multiple ray arrivals.

The example below uses a positive orthotropic stiffness in normalized units: $\rho =1$, $C_{11}=C_{22}=C_{33}=4$, $C_{44}=C_{55}=1.2$, $C_{66}=1$, and all normal off-diagonal stiffnesses zero. In the symmetry plane $n_3=0$, its two in-plane acoustic tensor eigenvalues are $v_{\pm }^2=(5\pm \sqrt{5+4\cos 4\theta })/2$, and its out-of-plane shear speed is $\sqrt{1.2}$. The in-plane shear sheet is nonconvex and produces the marked cusps; the construction uses an actual stable elastic tensor.

Use $p$ for disturbance pressure divided by the constant inertial density. The ideal, nonrotating Linearized Boussinesq equations areHere $N^2>0$ is a prescribed time-independent buoyancy frequency squared. Taking the scalar product of momentum with velocity gives ${\partial }_t(|\mathbf{u}{|}^2/2)=-\nabla \cdot (p\mathbf{u})+\sigma w$, using incompressibility. Multiplying the buoyancy equation by $\sigma /N^2(z)$ gives ${\partial }_t[{\sigma }^2/(2N^2)]=-\sigma w$. No spatial derivative of $N$ is needed in this cancellation. ThusThis is positive disturbance energy per unit inertial density; its flux is pressure work.

For constant $N$, write a plane internal gravity wave with wave vector $\mathbf{K}=({\mathbf{k}}_h,m)$ and phase $\theta ={\mathbf{k}}_h\cdot {\mathbf{x}}_h+mz-\omega t$. Let $k_h=|{\mathbf{k}}_h|>0$ and $K^2=k_h^2+m^2$. In complex amplitudes, incompressibility and horizontal momentum giveConsequently $\hat{p}=-\omega m\hat{w}/k_h^2$. Substitution into vertical momentum, $-i\omega \hat{w}=-im\hat{p}+\hat{\sigma }$, yieldsThe same dispersion relation follows directly from the displacement equation: its two terms become $K^2{\omega }^2\hat{\zeta }$ and $-N^2{k}_h^2\hat{\zeta }$. The positive-frequency branch is $\omega =N{k}_h/K$; the other branch is its sign reversal.

Choose real vertical-displacement amplitude $A$. The internal-wave polarization can be writtenVelocity and displacement lie along the phase planes, perpendicular to $\mathbf{K}$, and are a quarter cycle out of phase. Buoyancy is opposite in sign to vertical displacement. In the requested two-dimensional sketch take $k<0$, $m<0$ and $\omega >0$: the phase velocity $\omega \mathbf{K}/K^2$ points down and left. Then $\mathbf{e}=(-m/k,1)$ points up and left, and positive $p$ accompanies positive $w$ and negative $u$. Negative-pressure bands have the opposite velocity. Multiplication by pressure therefore transports energy up and left in both signs of band.

Differentiating the positive-frequency dispersion relation gives the group velocityIn particular $\mathbf{K}\cdot {\mathbf{c}}_g=0$. Since $|\mathbf{e}{|}^2=K^2/k_h^2$, the wavelength averages areUsing $N^2={\omega }^2{K}^2/k_h^2$ in the last expression provesFor $k,m<0$ the group velocity points up and left, perpendicular to the down-left phase velocity. The pressure sign pattern, rather than the direction of moving crests, determines the energy-flux direction shown below.

For the moving-boundary problem, first take the fluid to occupy $z>0$ with the boundary near $z=0$. An outgoing boundary-forced internal gravity wave must have $c_{gz}>0$. With the specified $k>0$ and $\omega >0$, the radiation condition therefore selectsIf instead the fluid occupies $z<0$, radiation requires $c_{gz}<0$ and the positive root for $m$. The fluid half-space is needed to select a unique sign; the following formulas use the upper half-space.

For slowly modulated internal-wave boundary forcing, put $\zeta =A(T,Z)e^{i\theta }$, $T=\mu t$, $Z=\mu z$. Acting on $A$ by ${\partial }_t\mapsto -i\omega +\mu {\partial }_T$ and ${\partial }_z\mapsto im+\mu {\partial }_Z$, the displacement equation has leading term $(K^2{\omega }^2-N^2{k}^2)A$. Its order-$\mu$ term isThus characteristics carry the boundary amplitude into the fluid:Real parts are understood. To the same leading order,Slow-envelope corrections to these polarization formulas are order $\mu$. At $T=2$, the amplitude is $ϵ$ for $0\le Z\le c_{gz}$, falls smoothly across $c_{gz}<Z<2c_{gz}$, and is zero beyond $2c_{gz}$. This translating ramp is also shown in the figure.

The mean boundary work per unit horizontal area is $\overline{W}={\rho }_{00}\overline{p(0)w(0)}$ at quadratic disturbance order. The small boundary displacement permits evaluation at the unperturbed plane. From the flux calculation,The finite causal wave train has energy integralDifferentiation, or integration of the local energy law with zero flux beyond the front, proves

Finally, distinguish the oscillatory velocity from the second-order Eulerian mean $\overline{u}$. The displacement polarizations giveThe supplied mean-flow relation therefore yields ${\overline{u}}_t=(k/\omega ){\overline{E}}_t+O({\mu }^2)$ at quadratic wave order. Equivalently, the horizontal momentum flux is $\overline{uw}=(k/\omega )c_{gz}\overline{E}$, and envelope-energy transport gives the same acceleration through ${\overline{u}}_t=-{\partial }_z\overline{uw}$.

In the frame moving horizontally at $c=\omega /k$, the boundary pattern is stationary and the initial mean current is $-c$. Subtract that background's kinetic energy before integrating over the infinite half-space. The internal-wave work in a phase-speed frame is thenIts rate of change isThe mean-current cross term supplies the quadratic energy loss in this moving frame; the lab-frame mean kinetic energy alone is fourth order in wave amplitude. The equality is at the leading slow and quadratic orders used throughout the forcing calculation.

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.