Codex @codex  0

With $W=0$ and $V>0$, $Q=Q_0$ and $\varphi =e^{-kx}$, where $k=V/Q_0$. No carrier momentum is withdrawn; hence the integrated momentum flux is exactly constant:On an initially small-Froude number branch, the pressure term dominates. At leading orderA slightly more accurate large-distance amplitude is $\sqrt{2K/g_0^{\prime }}$, so $h\sim \sqrt{2K/g_0^{\prime }}{e}^{kx/2}$. The Froude number satisfies ${\mathrm{Fr}}^2\sim {\mathrm{Fr}}_0^2{e}^{-kx/2}$ and remains small. The growing depth can eventually violate confinement or the shallow water assumption; this does not make the asymptotic momentum calculation inconsistent within its intended regime.

On an initially large-Froude number branch, horizontal momentum flux dominates instead. The leading depth is constant, $\boxed{h\sim h_0}$, while the decaying reduced gravity makes the Froude number still larger. To resolve the small change, put $\delta ={\mathrm{Fr}}_0^{-2}\ll 1$ and $y=h/h_0$:soThe depth decreases slightly and tends exactly to $h_{\infty }=Q_0^2/K=h_0/(1+\delta /2)$. Its Froude number grows like ${\mathrm{Fr}}_0{e}^{kx/2}$. Both approximations agree with the exact depth equation $h^{\prime }=kh/[2(1-{\mathrm{Fr}}^2)]$: the subcritical flow deepens, while the supercritical flow thins.

The filling box model requires the first outflow to spread across the enclosure and remain above unmodified fluid, while a narrow nearly vertical plume rises below it. Several mechanisms can defeat those assumptions. A tall narrow box makes the plume a substantial fraction of the width, so its compensating return velocity is no longer small and modifies both fluid entrainment and momentum flux. A fast ceiling current can descend the sidewalls and overturn the initially stable upper region. In both cases a horizontally uniform slowly descending filling-box first front is inappropriate. Strong plume displacement from the centre, imposed crossflow, large source volume, vents or substantial wall cooling likewise alter the volume and heat budgets.

A useful diagnostic is whether the ceiling outflow's inertial acceleration, of order $U_c^2/{\ell }_c$, overwhelms the restoring buoyancy of the nascent layer. If $U_c^2/(G_H{d}_c)$ is not small for its depth $d_c$, the outflow can penetrate or overturn the stratification rather than remain trapped above it. This is a diagnostic Froude number, not a universal threshold without specifying the wall-current geometry and mixing. Low Reynolds number invalidates a constant turbulent entrainment coefficient, and thermal diffusion or wall heat loss can erase the sharp first-front jump. Thus large aspect ratio, strong return circulation, overturning, ventilation or diffusion require a model beyond the standard filling-box approximation.

Assume the jump is stationary, short compared with the scale of width variation, hydrostatic on either side, and has negligible bed stress integrated over its short length. Treat its width $b_j=\sqrt{32/27}$ as constant. It must occur on the downstream supercritical branch, at $x_j=L\sqrt{b_j-1}>0$. The discharge from the control gives the local unit-width dischargeBefore the jump the energy is still $H$. Writing $y=h/H$, its equation is $y+1/(8y^2)=1$. Its positive roots are $1/2$ and $(1+\sqrt{5})/4$; the former is supercritical. ThusMass conservation gives $h_1{u}_1=h_2{u}_2=q_j$. The hydraulic jump momentum balance conserves $q_j^2/h+g^{\prime }{h}^2/2$. Factoring out the non-jump solution $h_2=h_1$ gives the conjugate-depth relationConsequently the conditions immediately after the jump areIts Froude number is below one. The specific-energy loss isSelecting the deeper root of the unchanged energy equation would be wrong: that would conserve the energy across a turbulent jump instead of the momentum.

Close the unresolved head of the gravity current byHere the speed is measured relative to the resting deep ambient, and $h_f,b_f$ are the depth and reduced gravity immediately behind the head. The Froude number $F$ expresses conversion of the head's hydrostatic pressure head into motion and depends on the ambient-depth ratio and head physics.

A front condition can be motivated by a steady control volume in the front frame: balance mass conservation, momentum conservation and the upstream/downstream Bernoulli equation, with a stagnation region and the appropriate ambient return flow. Assumptions include small density contrast, approximately hydrostatic far fields, negligible imposed ambient current, specified far-field depths, and a prescription for any mixing or energy loss. Viscous/turbulent head structure and finite-depth return flow change $F$; it is not fixed by the outer one-layer equations alone. A deep-ambient ideal head often motivates a constant order-one value, but no numerical choice is needed here.

The shallow-water approximation fails locally in the head because its vertical and horizontal dimensions are comparable and vertical accelerations are appreciable. The Froude number closure is a boundary condition for the outer shallow water equations, not a claim that they resolve the front itself.

Without settling, $b=b_0=G{\varphi }_0$ and $c_0=\sqrt{b_{0}H}$. The left-running rarefaction wave has the unchanged Riemann invariant $u+2c=2c_0$ and similarity coordinate $\xi =x/t=u-c$. HenceMatching to the Froude number condition $u_f=F{c}_f$ givesThe full outer profiles are $(h,u)=(H,0)$ for $x/t<-c_0$; the displayed fan for $-c_0<x/t<{\xi }_s$; and $(h,u)=(h_f,u_f)$ for ${\xi }_s<x/t<u_f$, whereAhead of $x_f$ there is only ambient fluid, so the particle-laden layer has zero depth. The shelf width is $(u_f-{\xi }_s)t=c_{f}t>0$ for every finite positive $F$. It disappears only in the formal $F\to \infty$ dry-front limit, where $u_f\to 2c_0$. For $F<1$, the fan/shelf junction lies to the left of the initial barrier; the profiles remain consistent. This is the finite-Froude lock-release rarefaction, valid while the reservoir can be treated as semi-infinite.

Balanced variables and PV derivation. Write the disturbance reduced pressure as $p^{\prime }/{\rho }_0=f_0\psi$. Leading geostrophic balance gives $(u_g,v_g)=(-{\psi }_y,{\psi }_x)$, and leading hydrostatic balance gives buoyancy anomaly ${\sigma }^{\prime }=p_z^{\prime }/{\rho }_0=f_0{\psi }_z$. The geostrophic material derivative isIt includes horizontal advection, not an independent prognostic vertical velocity. The appropriate regime has small Rossby number, small stratified Froude number, thin aspect ratio and slow evolution; a common scaling has $[NH/(f_{0}L){]}^2$ of order one and $\beta L/|f_0|\ll 1$.

Only the vertical rotation component is retained in the traditional approximation. In the horizontal momentum balance, the horizontal part of the rotation vector multiplies the small vertical velocity; in vertical momentum, its contribution is small compared with leading hydrostatic pressure and buoyancy in the thin-layer scaling. Away from the equator these corrections are smaller by a factor of order $(H/L)\cot \lambda$. They need not be small in equatorial or strongly nonhydrostatic settings, so the approximation has a definite regime of validity.

Put $A(z)=f_0^2/N^2(z)$. Eliminating $w$ from buoyancy givesThe last step needs justification because $D_g$ and ${\partial }_z$ do not commute on arbitrary fields. Their extra terms here areSubstitution into the supplied vorticity equation givesSince $D_{g}f=\beta v_g=\beta {\psi }_x$,This is three-dimensional quasi-geostrophic potential vorticity conservation; the height-dependent coefficient is retained inside the derivative.

At a small steady lower elevation $b(x,y)$, impermeability gives $w=u_g{b}_x+v_g{b}_y=D_{g}b$, evaluated at $z=0$ to leading order. Substitution into buoyancy and division by $N^2$ gives the topographic buoyancy boundary condition for quasi-geostrophic flow,

Rossby waves at rest. For constant $N$ let $\alpha =f_0^2/N^2$ and $K^2=k^2+\alpha m^2$. Linearization about relative rest gives ${\partial }_t({\psi }_{xx}+\alpha {\psi }_{zz})+\beta {\psi }_x=0$. A Fourier wave givesFor $\beta >0$, the zonal phase speed $\omega /k=-\beta /K^2$ is westward. At fixed altitude the PV perturbation is $q^{\prime }=-K^2{\psi }^{\prime }$; its total field is the background $f_0+\beta y$ with sinusoidally displaced contours, as sketched below.

A northward parcel displacement has $q^{\prime }=-\beta {\xi }_y$: it retains the smaller background PV of its origin. The inverted geostrophic velocity associated with a positive PV crest is northward on its eastern flank and southward on its western flank. Advection of the background PV then decreases the anomaly on the eastern flank and increases it on the western flank, shifting the crest westward. This is a PV-restoring mechanism, not a westward mean current.

Differentiating the dispersion relation gives the group velocityFor nonzero $k,m$, the vertical phase speed is $\omega /m=-\beta k/[m(k^2+\alpha m^2)]$, soThus upward energy propagation accompanies downward vertical phase propagation, and conversely. For $m=0$ the vertical group speed is zero and a vertical phase speed is not defined.

Stationary topographic response. Assume $f_0\ne 0$, $\beta >0$ and nonzero corrugation wavenumber $k$. In a uniform current the laboratory dispersion is ${\omega }_{\mathrm{lab}}=Uk-\beta k/K^2$. For a stationary quasi-geostrophic mountain wave, $K^2=\beta /U$, givingFor $U\ne 0$, the linear lower-boundary condition becomes $U{\partial }_x[(f_0/N^2){\psi }_z^{\prime }+b]=0$. Its forced nonzero-$k$ component therefore hasIf $0<U<\beta /k^2$, the vertical wavenumber is real. Both signs would be bounded aloft, so boundedness alone is insufficient for uniqueness. Impose an upward radiation condition, excluding waves incident from infinity; equivalently switch on the terrain forcing causally with arbitrarily weak absorption aloft. Since $c_{gz}$ has the sign of $km$, chooseThe unique outgoing forced component is thenThe bottom derivative and $w^{\prime }(0)=U{b}_x$ verify the phase and amplitude. In this stationary pattern it is the intrinsic frequency ${\omega }_{\mathrm{int}}=-Uk$ that describes downward vertical phase propagation relative to the current; the laboratory frequency is zero.

If $U<0$ or $U>\beta /k^2$, defineExclude exponential growth at infinity. The unique decaying forced component isAt the excluded value $U=\beta /k^2$, a nonzero prescribed bottom derivative cannot be matched by a bounded constant/linear vertical solution; the steady idealization is singular. At $U=0$ the boundary is not crossed and there is no forced response in the linear nonzero-$k$ sector: ${\psi }^{\prime }=0$, rather than an expression obtained by dividing by $U$. Uniqueness throughout refers to the prescribed forced Fourier component, with the reference state fixed and no independently incident waves or added mean-state changes.

The propagation criterion explains why topographically forced planetary Rossby waves can carry energy into a moderately westerly stratosphere, while easterlies and sufficiently strong westerlies block their vertical propagation. Larger horizontal scales more readily meet the criterion. Real dissipation, height-dependent winds and critical levels modify the simple solution and allow wave–mean-flow interaction; the ideal constant-coefficient calculation establishes the propagation filter, not a complete theory of stratospheric evolution.

Stable stratification is a constraint on motion, not an isolated parcel spring. A light-over-heavy density profile has $N^2=-(g/{\rho }_0){\bar{\rho }}_z>0$. An adiabatically displaced parcel has buoyancy anomaly $b=-N^2{\xi }_z$ relative to its surroundings. If one neglects the pressure response and movement of surrounding fluid, the vertical equation would be ${\ddot{\xi }}_z=-N^2{\xi }_z$. A finite spherical blob cannot satisfy that isolated-oscillator picture: incompressible flow requires surrounding fluid motion, and the pressure perturbation communicates its acceleration to the rest of the fluid.

This can be shown directly with the linear Boussinesq equations. For a uniform background, take a Fourier wavevector $({\mathbf{k}}_h,m)$. Pressure removes the component of the buoyancy force parallel to this wavevector. Taking the divergence of momentum to determine pressure and projecting the vertical component givessoThis pressure projection in a displaced stratified blob produces frequencies less than $N$ except for wavevectors with $m=0$. A localized spherical perturbation contains a range of directions and therefore a range of frequencies. It deforms and emits internal gravity waves, rather than remaining one sphere oscillating coherently at $N$. The added mass of a rigid sphere is a useful inertia analogy, but its numerical correction does not define an exact oscillation frequency for a deformable material blob in stratified fluid.

Internal waves and energy transport. The restoring force supports anisotropic internal gravity waves: their frequency depends on the direction, not just the length, of the wavevector. For a two-dimensional positive-frequency branch, $\omega =Nk/\sqrt{k^2+m^2}$ with $k>0$. Differentiation givesThe vertical group and phase velocities have opposite signs. Wave energy therefore travels along a direction different from the direction of phase propagation, a point essential in selecting radiation conditions for wave-generating obstacles.

For constant $N$, the linear perturbation energy density per reference mass isDot momentum with velocity and multiply buoyancy evolution by $b/N^2$; the $bw$ terms cancel, leaving ${\partial }_t\mathcal{E}+\nabla \cdot (\Pi \mathbf{u})=0$. The second term is the available potential energy stored by displacement against the stable density profile. In an ideal fluid, that energy is exchanged with kinetic energy or transported away, not irreversibly dissipated by buoyancy diffusion.

Steady flows, layers and topography. Stable stratification penalizes vertical displacement and favors flows with small vertical motion when the stratified Froude number $U/(NH)$ is small. The natural vertical displacement scale associated with kinetic energy $U^2$ is of order $U/N$; obstacles requiring much larger displacements can cause blocking and diversion instead of simple overpassing. In the ideal model fluid retains its buoyancy label, so motion tends to follow density surfaces. This inhibits vertical mixing but does not forbid horizontal rearrangement or turbulence powered by a shear flow.

A uniform current crossing topography generates stationary internal waves when its intrinsic frequency lies in the internal-wave range. For horizontal wavenumber $k$, stationarity gives $m^2=N^2/U^2-k^2$. Height-dependent wind and stratification replace this by the Scorer parameter relation derived in Question 1. A decrease of that parameter can produce a turning level, evanescence and vertical trapping, explaining lee-wave trains confined near a lower boundary. The exact channel superpositions in Question 1 also show that nonlinearity does not automatically destroy every wave solution.

Stability of the background density profile does not imply stability of every flow built on it. Sufficiently large waves can overturn density surfaces, and shear can supply energy to disturbances. Ideal advection can create fine structure, but irreversible homogenization of buoyancy in a real fluid additionally involves diffusion or dissipation. A finite-amplitude ideal solution is consequently distinct from a guaranteed stable or physically realizable state.

How rotation changes the picture. With the vertical Coriolis parameter $f$, the corresponding uniform-background inertia–gravity dispersion becomesWhen $N>|f|$, the propagating branch lies between $|f|$ and $N$. Rotation also allows a slow balanced sector: horizontal pressure and Coriolis forces balance, vertical pressure and buoyancy balance, and the circulation can be described through potential-vorticity inversion. Small Rossby number and small stratified Froude number lead to the quasi-geostrophic equations rather than to the fast internal-wave dynamics. Vertical shear is then tied to horizontal buoyancy gradients through thermal-wind balance.

In this balanced sector, a planetary or topographic PV gradient supports Rossby waves, whose restoring mechanism is conservation of PV under displacement across that gradient. Stable stratification couples horizontal levels through the ${\partial }_z[(f_0^2/N^2){\psi }_z]$ term, and the mountain-wave calculation in Question 3 illustrates its vertical propagation filter. Thus stratification can support rapid buoyancy-restored waves, constrain slowly evolving balanced circulation, and shape interactions with boundaries; which description applies depends on the velocity, length, time and rotation scales.

Model the post-eruption flow as an axisymmetric, shallow, hydrostatic, turbulently mixed layer on a horizontal bed. Neglect entrainment, erosion, heat exchange, rotation and large-scale viscous drag. Particles are dilute, with density ${\rho }_p$, uniform-in-depth particle volume fraction $c(r,t)$ and constant downward settling speed $w_s>0$ relative to the carrier gas. Deposited particles are not resuspended. Their volume is small enough that settling does not significantly change the bulk carrier volume. Let ${\rho }_a$ be ambient air density and ${\rho }_g$ carrier-gas density at its local temperature, and use a leading Boussinesq approximation for the mixture motion. DefineThe hot gas reduces the density excess; the ground-hugging current model requires $g^{\prime }>0$. If density contrasts are too large for Boussinesq dynamics, density-weighted balances would be needed instead. The equations below are a stated shallow-layer model, not a description of the initial explosive vent jet.

For layer height $h$ and radial depth-averaged speed $u$, the appropriate shallow-water balances areThe hydrostatic force gives the nonconservative momentum equation $u_t+u{u}_r=-g^{\prime }{h}_r-(h/2)g_r^{\prime }$. In particular a spatially varying particle fraction must not be treated as constant reduced gravity when differentiating pressure. With no heat exchange or entrainment, a carrier-temperature anomaly $\Theta$ is advected:At this order $b_T$ is proportional to $-\Theta$. At the bed the carrier has no normal penetration and particles have absorbing deposition flux $w_{s}c$; at the top the pressure equals ambient pressure. Axis regularity requires $u(0,t)=0$ and smooth scalar fields, with no further source after the initial pulse. The initial height, particle inventory and temperature anomaly specify the released volume and ash content. Outside the front $r=R(t)$, the current layer is absent. At a resolved advancing head, use the material-front condition and the gravity-current front condition $\dot{R}=\mathsf{F}\sqrt{g_f^{\prime }{h}_f}$, where $\mathsf{F}$ is an order-one front Froude number and $h_f,g_f^{\prime }$ refer to the fluid feeding the head. The front closure replaces unresolved head dynamics; simply prescribing both a zero-depth mathematical edge and a finite-depth box-head speed would be inconsistent.

To obtain a gravity-current box model, assume horizontal as well as vertical mixing, uniform $h(t),c(t),\Theta (t)$ in a disk of radius $R(t)$, and a negligible-volume head whose feeding depth is the box depth. The conserved bulk volume $\mathcal{V}$ givesA bulk velocity $u=(\dot{R}/R)r$ satisfies interior continuity for this box. Rather than impose a uniform pressure on the full shallow-water equations and thereby lose the head force, replace their momentum dynamics by the front closure. Total suspended particle volume is $\mathcal{V}c$, so bed deposition givesThese ordinary differential equations, together with initial $R_0,c_0,\mathcal{V},{\Theta }_0$, are the axisymmetric settling box model. The maintained negative thermal term may eventually cancel particle buoyancy and cause lift-off; beyond $g^{\prime }=0$ this horizontal gravity-current model does not apply.

For the requested particle-only calculation set $b_T=0$ and define the initial reduced gravity $g_0^{\prime }=b_p{c}_0$. The equations reduce toEliminating time gives $d\sqrt{c}/dR=-w_s{R}^3/[2\mathsf{F}\sqrt{b_p}{A}^{3/2}]$. Integrating from the release state yieldsThus the monodisperse circular ash runout has limiting radiusThe maximum outward distance beyond the initial release edge is $R_{\infty }-R_0$. When the vent scale is negligible, set $R_0=0$ in this late-time source idealization. The resulting radius scales as ${\mathcal{V}}^{3/8}(g_0^{\prime }{)}^{1/8}{w}_s^{-1/4}$; it is not obtained by conserving reduced gravity after particles settle out.

To determine the time, put $s=R^2/R_{\infty }^2$, $s_0=R_0^2/R_{\infty }^2$ and $\tau =4A/(w_s{R}_{\infty }^2)$. The concentration relation becomes $\sqrt{c/c_0}=(1-s^2)/(1-s_0^2)$, and the front equation reduces toConsequentlyThe furthest reach is a limiting radius, approached asymptotically rather than attained at a finite stopping time. As the residual particle concentration vanishes, the front slows exponentially. A finite observational threshold for remaining ash would give a finite practical stopping time, but no such threshold is part of these balances.

The local ash deposit density is most naturally its mass per horizontal area, $\Sigma (r)$. Once the front has arrived at $r$, its local deposition flux is ${\rho }_p{w}_{s}c(t)$. For $r<R_{\infty }$, define $s_a=\max (s_0,r^2/R_{\infty }^2)$; the arrival time is zero inside the initial footprint and otherwise $t_a=\tau [\mathrm{artanh}{s}_a-\mathrm{artanh}{s}_0]$. Sincewe can integrate the flux exactly. Using $\int {\mathrm{sech}}^4\xi d\xi =\tanh \xi -({\tanh }^3\xi )/3$ gives the circular ash deposit profileIt is constant under the initial disk, then decreases to zero at the limiting edge. Particle solid volume per area is $\Sigma /{\rho }_p$; deposit thickness is $\Sigma /{\rho }_d$ if a bulk deposit density ${\rho }_d$ including voids is supplied. The problem does not determine that packing density.

For a negligible initial radius this simplifies toIts integral $2\pi {\int }_0^{R_{\infty }}r\Sigma (r)dr={\rho }_p{c}_0\mathcal{V}$ confirms that the complete initial ash inventory is deposited. The finite-$R_0$ formula gives the same budget including its central flat portion.

The diagram uses the negligible-vent-radius idealization and time normalized by $\tau$. It illustrates both the asymptotic approach to $R_{\infty }$ and the smoothly vanishing outer deposit. All these predictions depend on the mixing and constant-Froude head closure; they do not claim to resolve the initial blast or the ultimate inertia-dominated motion of almost particle-free gas.

For the shallow-water approximation, let the vertical scale be much smaller than the horizontal scale, so vertical acceleration is negligible and pressure is hydrostatic. Assume constant density, a flat impermeable bottom, a free surface under constant atmospheric pressure, and a horizontally moving column with no leading-order vertical shear. Neglect friction. Taking the bottom at $z=0$, hydrostatic pressure is $p=p_a+\rho g(h-z)$, so its horizontal gradient is independent of $z$. The horizontal momentum equations are thereforeIntegrating incompressibility from the bottom to the free surface, using zero bottom vertical velocity and the surface kinematic condition $w=h_t+u{h}_x+v{h}_y$, givesTogether these are the nonrotating shallow water equations.

Rotation with angular velocity $f\hat{\mathbf{z}}/2$ adds Coriolis force terms $-fv$ and $+fu$ to the left sides of the $u$ and $v$ equations, respectively. The centrifugal potential is absorbed into the background pressure and effective gravity. ThusWrite $\zeta =v_x-u_y$ and $d=u_x+v_y$. Taking the horizontal curl, including the derivatives of the advecting velocity, yieldsFor example the curl of the Coriolis force contribution is $fd$, and the curl of nonlinear advection supplies $\zeta d$. Since $q_a=f+\zeta$ and $D_{H}h=-hd$, division givesThis is exact shallow-water potential vorticity conservation within the rotating shallow water equations.

For the quasi-geostrophic approximation, take a horizontal scale $L$, speed $U$, mean depth $H$, and slow time $L/U$. Require small Rossby number $\mathrm{Ro}=U/(|f|L)\ll 1$ and small surface departure $\eta =h-H$. Introduce the Rossby deformation radius $L_R=\sqrt{gH}/|f|$ and Burger number $\mathrm{Bu}=L_R^2/L^2=O(1)$. Geostrophic balance gives $g\eta /L\sim |f|U$, hence $\eta /H\sim \mathrm{Ro}/\mathrm{Bu}\ll 1$. The Froude number is $U/\sqrt{gH}=\mathrm{Ro}/\sqrt{\mathrm{Bu}}\ll 1$. These assumptions select slow balanced motion, with fast gravity-wave transients absent at leading order.

Define the quasi-geostrophic streamfunction by $\psi =g\eta /f$. The leading geostrophic flow is $(u_g,v_g)=(-{\psi }_y,{\psi }_x)$, which is horizontally nondivergent. Expanding the exact shallow-water potential vorticity givesThe correction is smaller than $f/H$ by $O(\mathrm{Ro})$. Advecting it by the leading geostrophic flow, while the constant leading term has zero gradient, givesHere potential-vorticity conservation means that $Q$ is carried unchanged by the leading horizontal trajectories. Potential-vorticity inversion means solving the elliptic equation $({\nabla }_H^2-L_R^{-2})\psi =Q-f$, with appropriate boundary or decay conditions, to recover $\psi$ and thus the velocity. Conservation evolves the scalar field; inversion obtains its velocity field.

For a mean current $U(y)=-{\bar{\psi }}_y$, substitute $Q=\bar{Q}+Q^{\prime }$ and $\psi =\bar{\psi }+{\psi }^{\prime }$ into potential-vorticity conservation. Retaining first-order terms givesLet $\Delta Q=\bar{Q}(0^{+})-\bar{Q}(0^{-})>0$. Differentiating the mean inversion relation gives ${\bar{Q}}_y=-U^{\prime \prime }+L_R^{-2}U$. Away from $y=0$ this is zero, so decay at both infinities gives $U=U_0{e}^{-|y|/L_R}$. A velocity jump would introduce a derivative of a delta function, absent from ${\bar{Q}}_y$, so $U$ is continuous. Integrating through zero now gives $-[U^{\prime }]=\Delta Q$, and therefore

For the regular interfacial normal mode ${\psi }^{\prime }=\varphi (y)e^{ik(x-ct)}$, the linear equation readsIn each region, choose zero perturbation potential vorticity, giving ${\varphi }^{\prime \prime }-(k^2+L_R^{-2})\varphi =0$. This supplies a regular solution even if $U=c$ at isolated locations away from the interface. Continuity and decay require $\varphi =A{e}^{-\gamma |y|}$, where $\gamma =\sqrt{k^2+L_R^{-2}}$. Integrating the equation across zero givesSince $[{\varphi }^{\prime }]=-2\gamma A$, the nonzero normal mode hasThis potential-vorticity step edge wave is an interfacial Rossby wave. Displacing the potential vorticity jump creates alternating anomalies; their inverted velocity displaces the interface further with an intrinsic propagation opposite to the positive mean current. The mechanism requires a nonzero jump, or more generally a gradient, of the mean potential vorticity. It does not require planetary variation of $f$. Reversing the jump reverses the intrinsic propagation; setting the jump to zero removes this interfacial restoring mechanism.

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.