Codex @codex  0

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.

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.

The phase velocity is the motion of a constant-wave phase surface. For a one-dimensional plane wave $u=a\cos (kx-\omega t+{\varphi }_0)$ with real $k\ne 0$ and real $\omega$, fixing $kx-\omega t+{\varphi }_0$ givesEach wave crest of this monochromatic travelling wave advances at that signed phase velocity. Its magnitude is the phase speed $|c_p|=\lambda /T$, where $\lambda =2\pi /|k|$ is the wavelength and $T=2\pi /|\omega |$ is the temporal period when $\omega \ne 0$. Choosing $\omega >0$ makes positive $k$ mean rightward phase velocity and negative $k$ mean leftward phase velocity. Reversing both signs leaves the real plane wave and its phase velocity unchanged. At $k=0$ there are no spatial crests to track, so $\omega /k$ is undefined; a finite limit as $k\to 0$ is a limiting phase speed, not the speed of a spatially uniform oscillation.

In several dimensions, a plane wave has wave phase $\mathbf{k}\cdot \mathbf{x}-\omega t+{\varphi }_0$. Constant-wave phase planes obey $\mathbf{k}\cdot \mathbf{v}=\omega$. Their normal phase velocity isTangential motion within a plane does not move that wavefront, which is why the normal direction fixes the vector convention. In an anisotropic medium the group velocity need not point along this wavevector. Likewise the horizontal phase velocity of a tilted wavefront is its crossing speed along a horizontal line, rather than a component of the normal phase velocity.

For the nondispersive wave branch $\omega =ck$, a wave crest and a wave envelope both move at $c$. For deep-water surface gravity waves, $\omega =\sqrt{gk}$ gives $c_p=\sqrt{g/k}$ and group velocity $c_g=c_p/2$. With the illustrative choices $g=9.81\mathrm{m}{\mathrm{s}}^{-2}$ and $\lambda =1\mathrm{m}$, the model predicts $c_p\simeq 1.25\mathrm{m}{\mathrm{s}}^{-1}$ and $c_g\simeq 0.625\mathrm{m}{\mathrm{s}}^{-1}$. These are calculated values for negligible surface tension and sufficiently large $kh$, not measurements. A wave crest moves forward through the wave packet, so timing a crest does not time the packet's journey.

There is a direct laboratory meaning to this distinction. MIT's 2005 surface-wave lab uses a long water tank, a controllable hydraulic paddle, and an absorbing beach. Resistance probes convert submerged length into voltage; calibration converts their readouts into free surface elevation versus time. Two probes a known distance apart, backed by a ruler and stopwatch, let students time a wave crest and separately time a feature of the wave packet. The lab records depth and varies frequency and wave amplitude.^[ref]

For a clean measurement of phase velocity, fit the periodic portion of both records or follow one identifiable wave crest. A phase difference is only known modulo $2\pi$: probe spacing below half the shortest relevant wavelength, or continuity while sweeping frequency, avoids assigning the wrong whole-cycle delay. Then compare the timed phase velocity with $\lambda /T$. Tracking pepper grains instead measures fluid-particle motion, which has a different meaning from either phase velocity or group velocity.

The Klein-Gordon packet velocities example makes the distinction especially sharp: its phase velocity exceeds $c$, yet its narrow wave envelope moves below $c$. A carrier wave labels successive oscillations; its speed alone does not establish information velocity. The front velocity depends on the full causal model, including angular frequencies far outside the carrier's narrow range.

Choose $k,m>0$, $K=(k^2+m^2{)}^{1/2}$ and the positive angular frequency $\omega =Nk/K$. Then the phase velocity is $Nk(k,m)/K^3$ and the group velocity is $N(m^2,-km)/K^3$. Their dot product vanishes and their sum is horizontal with magnitude $N/K$. Draw the group velocity from the tip of the phase velocity: the two form the perpendicular sides of a right triangle. The horizontal hypotenuse is their sum, not the horizontal phase velocity $\omega /k$.