Codex @codex  0

The hydrostatic buoyancy anomaly is $b^{\prime }=f_0{\psi }_z$, and its leading adiabatic evolution is $D_g{b}^{\prime }/Dt+N_0^{2}w=0$. The small vertical velocity advects the background stratification; vertical advection of the perturbation buoyancy is higher order. ConsequentlyThe impermeability condition at a stationary lower boundary $z=\alpha y$ is $w=\alpha v$. Using $v_g={\psi }_x$ gives the topographic buoyancy boundary condition for quasi-geostrophic flowFor application at the reference plane, the actual small-height condition is $|\alpha |L/D\ll 1$. Slopes scaled by $D/L$ are of Rossby-number order, so this topographic contribution is retained alongside the small ageostrophic vertical velocity. The dimensional slope should not be compared to the Rossby number without also specifying the aspect ratio.

The basic quasi-geostrophic streamfunction is $\overline{\psi }=-\Lambda zy$. It gives $U=\Lambda z$, basic buoyancy gradient ${\overline{b}}_y=-f_0\Lambda$, and spatially constant interior quasi-geostrophic potential vorticity. Linearizing its conservation law therefore givesThe printed equation omits the shear-advection term. This cannot be omitted for general disturbances: an interior anomaly proportional to $\cos [k(x-\Lambda zt)]$ is materially conserved but not independent of time at a fixed location. The intended zero-interior-PV Eady model modes do have $q^{\prime }=0$ and satisfy both expressions; the ensuing mode calculation is valid in that subspace.

At either boundary, linearizing the geostrophic material derivative includes $v^{\prime }{\overline{\psi }}_{zy}=-\Lambda {\psi }_x^{\prime }$. LetThe two boundary equations areThese are exactly the respective bottom and top slope corrections to the rigid-boundary buoyancy condition for quasi-geostrophic waves.

Assume periodicity horizontally and first consider $q^{\prime }=0$. Multiply its defining equation by ${\psi }_x^{\prime }$ and integrate over the layer. The two horizontal Laplacian terms vanish by integration by parts; the remaining volume term after one vertical integration is an $x$-derivative. ThereforeFor nonzero $A_0,A_D$, the boundary equations giveThe top advection term integrates to zero. Subtracting proves conservation of the boundary pseudomomentum of a sloping Eady layer:The original PDF uses squared vertical derivatives and the upper slope $\gamma$ in the first denominator; these details are corrupted in the converted TeX. For arbitrary nonzero interior PV the same calculation instead gives $d\mathcal{P}/dt=(f_{0}s{)}^{-1}\int {\psi }_x^{\prime }{q}^{\prime }dV$, which need not vanish. For example, at an instant take ${\psi }^{\prime }=\cos kx+(z/D{)}^2\sin kx$, independent of $y$; the zonal average of ${\psi }_x^{\prime }{q}^{\prime }$ is $-ks/D^2$. Thus the boundary-only conservation statement also needs the zero-interior-PV qualification.

For an exponentially growing smooth normal mode, this qualification is automatic: the interior equation is $ik(\Lambda z-c)\hat{q}=0$, and a phase speed with nonzero imaginary part cannot equal the real basic velocity. Therefore $\hat{q}=0$. If $A_0{A}_D<0$, the conserved boundary quadratic form is definite, so a nonzero exponentially growing boundary amplitude is impossible. Instability consequently requiresfor $\Lambda \ne 0$. When one $A$ vanishes, the corresponding boundary derivative is simply advected and must vanish for a growing mode. The other boundary then yields a real single-edge-wave speed; if both vanish, the zero-interior-PV Neumann problem has no nonzero growing mode. These degenerate cases are neutral and should not be handled by dividing by zero.

For equal slopes, put $\zeta =z/D$, $C=c/(\Lambda D)$, $M=N_0|k|D/|f_0|$, and $a=N_0^2\alpha /(f_0\Lambda )-1$. The normal-mode interior equation is ${\Phi }_{\zeta \zeta }-M^2\Phi =0$ for the zero-interior-PV modes. ThusThe coefficient equations areSetting their determinant to zero, without dividing by $C$ or $a$, givesso the Eady model with parallel sloping boundaries hasFor positive $k,f_0$, $M$ is the printed $\mu$, and $a$ is the printed $\tilde{\alpha }$. The expression is even in the sign of that vertical scale. For $a\ge 0$ the radicand is at least $1/4$, so both speeds are real and these modes have no exponential growth. For equal slopes the earlier necessary condition is merely $a^2>0$, which is satisfied in many stable cases and is not sufficient. At $a=0$ the speeds are $0$ and $\Lambda D$; the invariant with nonzero denominators cannot be used, but the determinant still correctly gives neutral modes. No claim excluding transient amplification is needed.

Let $V=4\pi b^3/3$ and $B=Vg({\rho }_0-\rho )/{\rho }_{\mathrm{ref}}$. The four balances describe, respectively, fluid entrainment through the spherical surface at normal speed $\alpha u$, the increase of vertical momentum due to buoyancy, conservation of total buoyancy in a uniform ambient, and motion of the thermal centre. Entrained ambient fluid initially has no vertical momentum, so the momentum equation includes the inertia needed to accelerate it. This model neglects added mass and drag; these effects would alter the numerical similarity coefficients.

With $|\rho -{\rho }_0|\ll {\rho }_0$, the Boussinesq approximation replaces inertial density by ${\rho }_0$ while retaining the density difference in buoyancy. Take ${\rho }_{\mathrm{ref}}={\rho }_0$ to this accuracy. ThenFor a point source, $b(0)=z(0)=0$ and the initial integrated impulse is zero. Therefore $b=\alpha z$ and $Vu=B_{0}t$. Integrating $b^3\dot{b}=3\alpha B_{0}t/(4\pi )$ yieldsThe mass density follows from conserved buoyancy:Thus the radius and rise height grow as $t^{1/2}$, the rise speed decays as $t^{-1/2}$, and the density deficit decays as $t^{-3/2}$. In height coordinates, put $K=4\pi {\alpha }^3/3$ and $C=\sqrt{B_0/(2K)}$; then $V=K{z}^3$ and $u=C/z$.

If a reference density different from the inertial ambient density is retained, $d(Vu)/dt=({\rho }_{\mathrm{ref}}/{\rho }_0)B_0$: replace $B_0$ in the kinematic formulas by $({\rho }_{\mathrm{ref}}/{\rho }_0)B_0$, while retaining the original $B_0$ in the density-deficit formula. This keeps reference-density normalization separate from the physical acceleration.

The point-source spherical thermal similarity is singular at $t=0$ and cannot satisfy the small-density-difference condition arbitrarily close to its formal origin. A finite-radius, zero-impulse regularization has $b=b_s+\alpha z$ and $b^4=b_s^4+3\alpha B_0{t}^2/(2\pi )$, with $u=B_{0}t/V$. This displays the required initial data rather than attributing the singularity to a physical infinite initial speed.

In stable density stratification, ambient density decreases with height. A rising buoyant thermal encounters progressively lighter surroundings and entrains denser fluid from below, so its net buoyancy decreases. It reaches neutral buoyancy, overshoots through its remaining momentum, and eventually stops rising, spreads as a buoyant intrusion or excites internal gravity waves. The finite rise is a consequence of this restoring effect, not just drag.

Total buoyancy has dimensions $[B_0]=L^4{T}^{-2}$ and buoyancy frequency has $[N]=T^{-1}$. For an ideal point release with no other dimensional source length or time, dimensional analysis therefore givesThe dimensionless coefficient depends on entrainment coefficient, shape and the momentum model. A finite source radius or impulse supplies additional dimensionless parameters; dimensional analysis cannot determine the coefficient by itself.

Within precisely the spherical model, $b=\alpha z$ and the evolving thermal buoyancy satisfies $dB/dz=-N^{2}V=-K{N}^2{z}^3$, so $B(z)=B_0-K{N}^2{z}^4/4$. Neutrality occurs at $z_n^4=4B_0/(K{N}^2)$. For $P=Vu$, the momentum equation and $\dot{z}=P/V$ giveThe first nonzero-height stopping point is $z_{\max }^4=8B_0/(K{N}^2)$, so $z_{\max }=2^{1/4}{z}_n$. This illustrates both the dimensional scale and the distinction between neutral height and maximum rise; the precise coefficients rely on the specified ideal fluid entrainment and impulse model.

The cell conservation equation equates the change of cell number density $n$ to minus the divergence of its total flux. The flux $n\mathbf{u}$ is advection by the water, $n{V}_s⟨\mathbf{p}⟩$ is the mean swimming flux, and $-\mathbf{D}\cdot \nabla n$ is the dispersal flux. The tensor $\mathbf{D}$ describes effective translational diffusion, whereas $D_r$ in part (i) describes rotational diffusion; they have different dimensions. This conservation law neglects birth, death and sources at the walls.

The equation $\nabla \cdot \mathbf{u}=0$ expresses incompressible flow. In the momentum equation, $\rho$ is the reference suspension density, $\rho [{\partial }_t\mathbf{u}+(\mathbf{u}\cdot \nabla )\mathbf{u}]$ is material acceleration times density, $-\nabla p_e$ is the pressure force after subtracting the pure-water hydrostatic pressure, and $\mu {\nabla }^2\mathbf{u}$ is the viscous force for constant dynamic viscosity $\mu$ in a Newtonian fluid. Since $n$ counts cells per volume, $g^{\prime }$ in $-g^{\prime }n{\mathbf{e}}_z$ is excess downward weight per cell; it is not an acceleration. The small cell-fluid density difference permits the Boussinesq approximation: retain that difference in buoyancy while using constant $\rho$ elsewhere. Diluteness permits neglect of cell-cell interactions, changes of viscosity, and active cell stresses. Any separate settling speed is assumed negligible or already incorporated in the specified mean drift. The light strength and swimming speed are uniform.

Consider the intended long-channel bulk approximation. At the sidewalls impose zero normal cell flux and the no-slip boundary condition. The closed ends require zero total fluid flux through each horizontal section, so ${\int }_{-a}^{a}wdx=0$; localized turning regions are omitted from the bulk equations. The mean concentration is ${\int }_{-a}^{a}ndx=2a{n}_0$. The final paragraph below distinguishes these bulk conditions from exact closed-end cell conservation.

For $\mathbf{u}=w(x){\mathbf{e}}_z$, the vorticity is $-w^{\prime }(x){\mathbf{e}}_y$. Hence $\epsilon =-w^{\prime }/(2D_r)$, and the leading horizontal mean swimming direction is $K_0{\mathbf{e}}_x$. With $\mathbf{D}=DI$ and $n=n(x)$, steady cell conservation equation says that $J_x=V_s{K}_{0}n-D{n}^{\prime }$ is constant. Its wall value is zero, so $n^{\prime }=\beta n$ with $\beta =V_s{K}_0/D$. The phototactic concentration in a vertical channel is thereforeThis result neglects the $O({\epsilon }^2)$ correction to horizontal swimming. For positive phototaxis, $\beta >0$ and cells accumulate towards the illuminated wall $x=a$.

The steady inertial terms vanish for $w=w(x)$, since $\mathbf{u}\cdot \nabla =w{\partial }_z$. The horizontal momentum equation gives ${\partial }_x{p}_e=0$. The vertical momentum equation then requires a constant excess pressure gradient $G={\partial }_z{p}_e$ and gives $\mu w^{\prime \prime }=G+g^{\prime }{n}_1{e}^{\beta x}$. Twice integrating producesSet $b=a\beta$. The two no-slip boundary conditions give $A=-C\sinh b/a$ and $B=-C\cosh b-G{a}^2/(2\mu )$. The zero fluid-flux condition follows by integrating this expression:Consequently the remaining constants areThe resulting velocity can also be written without separate polynomial constants asHere $G$ is the zero-flux pressure gradient in a cell-driven channel. It is negative for $b>0$: pressure decreases upwards, creating the upward pressure force that balances the excess cell weight and drives the return flow. Its magnitude is determined by zero net fluid transport, not by prescribing pressure independently. In the zero-light limit, the continuous result is $n=n_0$, $G=-g^{\prime }{n}_0$ and $w=0$.

For the requested large-$b$ sketch, cells occupy a boundary layer of thickness ${\beta }^{-1}$ near $x=a$. Put $\xi =x/a$ and $W=g^{\prime }{n}_0{a}^2/(\mu b)$. The exact dimensionless profile above isAway from the illuminated-wall layer it tends to $(1+\xi )(1-3\xi )/2$. There is broad upward return flow for $-1<\xi <1/3$ and downward flow on the illuminated side. The outer maximum is $2W/3$ at $\xi =-1/3$. Near $\xi =1$, the exponential term is approximately $2e^{-b(1-\xi )}$ and restores the no-slip boundary condition from the negative outer velocity to zero at the wall. The finite-$b$ profile shown below retains that layer. The small-vorticity assumption must still be enforced, for example by sufficiently weak excess weight or a sufficiently dilute suspension; $b\gg 1$ alone does not ensure $|w^{\prime }|/(2D_r)\ll 1$.

There is a genuine global qualification to the bulk construction. For an exactly steady channel with impermeable ends as well as impermeable sidewalls, the closed-channel cell-flux constraint also requires $\int J_{z}dx=0$ at every height. In this bulk ansatz the first-order orientation closure givessince integration by parts, $n^{\prime }=\beta n$ and $w(\pm a)=0$ imply $\int n{w}^{\prime }=-\beta \int nw$. Moreover, multiply $\mu w^{\prime \prime }=G+g^{\prime }n$ by $w$ and integrate. The wall values and zero fluid flux givefor a nontrivial profile. Thus zero axial cell flux generally fails unless the extra factor $1+V_s\beta K_1/(2D_r)$ vanishes, or the flow is trivial. A complete steady closed-channel solution generically needs axial concentration variation. The explicit profile solves the stated bulk equations, sidewall conditions, concentration normalization and zero fluid-flux condition; it is not, for arbitrary parameters, an exact globally closed steady cell distribution. Ignoring end regions gives the intended bulk approximation, but does not by itself prove the stronger global existence claim.

The geodynamo operates in the electrically conducting liquid of Earth's outer core. buoyancy from thermal and compositional convection supplies motion, while the Coriolis force organizes rotating flow. Rotation together with radial stratification permits kinetic helicity, typically with opposite signs in the two hemispheres. In mean-field language this allows an alpha effect that converts poloidal and toroidal components into one another. A large-scale alpha-squared dynamo is therefore a plausible regeneration mechanism even without dominant differential rotation.

The growth condition compares this regeneration with the sum of molecular and turbulent magnetic diffusion on the core scale. A spherical model must impose the appropriate magnetic matching at the core boundaries and permit spatially varying, generally anisotropic alpha and diffusivity tensors. Rotation selects an axis, and stratification and boundaries select positions, so homogeneous isotropic turbulence is only a local idealization. Such a model can favor a predominantly axial dipole, but the preferred parity and threshold must be found from its global eigenproblem rather than deduced from the plane-wave formula alone. If mean shear is significant, alpha–Omega regeneration can also contribute.

The first-order-smoothing helicity formula is a useful guide to the sign and magnitude of alpha, not an automatically controlled numerical closure for the core. The important convective motions need not have small eddy magnetic Reynolds number or short enough correlation time; scale separation can also be limited. Furthermore, the Lorentz force modifies the flow and its correlations once the field grows. Saturation therefore requires nonlinear feedback, including dynamo quenching and magnetic-helicity constraints where relevant.

The observed field varies and occasionally reverses. A constant-coefficient growing kinematic eigenmode does not itself explain these reversals or select a saturated amplitude. Spatial and temporal variations, competing modes, and nonlinear flow–field interactions must be included. mean-field dynamo theory supplies a physical regeneration mechanism and a testable threshold, while the geodynamo requires a rotating, bounded, nonlinear realization of that mechanism.

Equations and conventions. Let ${\rho }_0$ be a constant reference density and define total buoyancy acceleration by $\sigma =-g(\rho -{\rho }_0)/{\rho }_0$. With reduced pressure $\Pi =p/{\rho }_0+gz$, the nonrotating ideal Boussinesq equations areDensity departures are small compared with ${\rho }_0$, so they can be neglected in inertia and mass continuity but retained in the gravitational force. The reference stratification must also change density only weakly over the depth of interest. incompressible flow is already assumed; the Boussinesq approximation additionally restricts the size of the density variations, rather than following automatically from incompressible flow. An ideal fluid here has neither viscosity nor buoyancy diffusion.

A stationary reference profile $\bar{\sigma }(z)$ has buoyancy frequencyStable stratification has $N^2>0$: an upward-displaced parcel retains its old buoyancy and has negative buoyancy relative to its new surroundings.

Vorticity equation and linearization. For motion in the $(x,z)$ plane, choose the streamfunction by $u={\psi }_z$, $w=-{\psi }_x$. The $y$-component of vorticity is then $u_z-w_x={\nabla }^2\psi$, where ${\nabla }^2={\partial }_x^2+{\partial }_z^2$. Taking the $y$-component of the curl of momentum givesThe minus sign in the buoyancy curl follows from $(\nabla \times \sigma \hat{z}{)}_y=-{\sigma }_x$.

Write $\psi =\bar{\psi }(z)+{\psi }^{\prime }$ with ${\bar{\psi }}_z=\bar{u}(z)$ and $\sigma =\bar{\sigma }(z)+{\sigma }^{\prime }$. The linear material derivative is $D_t={\partial }_t+\bar{u}{\partial }_x$. The perturbation equations areThe shear-curvature term arises because $w^{\prime }$ advects the basic vorticity ${\bar{u}}_z$. Apply $D_t$ to the first equation and use the second; $D_t$ commutes with ${\partial }_x$ and does not differentiate the height-only coefficients. Therefore

For a stationary nonzero horizontal wavenumber $k$, the buoyancy equation gives $\hat{\sigma }=N^2\hat{\psi }/\bar{u}$. Substitution into the vorticity equation givesHere ${\ell }^2$ is the Scorer parameter. Where $m^2>0$, the vertical structure is oscillatory; where $m^2<0$, it is evanescent. A rigid flat lower boundary gives $\hat{\psi }(0)=0$ for $k\ne 0$. In an upper evanescent region the growing solution is excluded by boundedness. When propagation persists to infinity, an appropriate radiation condition distinguishes outgoing from incoming energy for a forced problem.

If $\bar{u}=U_0+Sz$ with $S>0$, its curvature term vanishes. Bounded $N^2$ implies ${\ell }^2=N^2/{\bar{u}}^2\to 0$ aloft. Thus every fixed nonzero $k$ eventually lies in an evanescent region. If stratification is sufficiently strong near the bottom, a lower oscillatory region is separated from the upper evanescent region by a turning level of a stationary internal wave. Waves reflect at that level and at the rigid lower boundary, giving vertically trapped modes for the horizontal wavenumbers that satisfy both conditions. In a slowly varying example the phase condition is approximately ${\int }_0^{z_t}mdz=(n-1/4)\pi$, from a lower node and the turning-point connection; the key point is trapping by the decrease in the Scorer parameter, not a physical lid.

Rigid channel modes. With constant $U=\bar{u}>0$ and $N$, the vertical nodes require $m_n=n\pi /H$, $n=1,2,\dots$. Stationarity then requiresThe specified strict interval admits positive $k_n^2$ for exactly $n=1,2$. Choosing $k_1,k_2>0$, every real superposition of these nonzero-wavenumber modes isThe two independent complex constants encode the amplitudes and horizontal phases; negative wavenumbers are supplied by complex conjugation, not additional independent constants.

Exact finite amplitude. Put $\lambda =N^2/U^2$. Both modes and every superposition satisfy ${\nabla }^2{\psi }^{\prime }=-\lambda {\psi }^{\prime }$. For the total streamfunction $\Psi =Uz+{\psi }^{\prime }$, the total buoyancy is $\sigma =N^{2}z+(N^2/U){\psi }^{\prime }=(N^2/U)\Psi$. A stationary velocity $({\Psi }_z,-{\Psi }_x)$ is tangent to its streamfunction contours, so $D\sigma /Dt=0$ exactly. Moreover the nonlinear self-advection of ${\nabla }^2{\psi }^{\prime }$ vanishes because it is proportional to ${\psi }^{\prime }$ itself. The remaining vorticity equation iswhich is an identity. Both rigid boundaries remain impermeable, and the curl-free residual in momentum determines a pressure. Thus the displayed two-mode family consists of exact finite-amplitude Boussinesq solutions, with no small-amplitude restriction on the constants in these equations. Very large amplitudes can overturn the stratification or leave the physical Boussinesq regime; exact existence does not establish stability.

The two-constant statement concerns wave disturbances with $k\ne 0$. If $k=0$ is literally included, arbitrary stationary horizontal shear adjustments ${\psi }^{\prime }=F(z)$ and height-only buoyancy changes are also possible. They are changes of the basic state and are not constrained by the wave equation obtained after division by $k$. Excluding that mean-flow sector is necessary for the intended count.

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.