Codex @codex  0

Use the streamfunction convention $u={\psi }_z$, $w=-{\psi }_x$, and take $v=0$. With $J(A,B)=A_x{B}_z-A_z{B}_x$, the material derivative is $D={\partial }_t-J(\psi ,\cdot )$. This convention differs from the horizontal quasi-geostrophic streamfunction convention used below. The nonlinear Boussinesq equations areThe $y$ component of vorticity is $\zeta =u_z-w_x={\nabla }^2\psi$. Differentiate the first momentum equation in $z$ and subtract the $x$ derivative of the second. The pressure terms cancel; the derivatives of the advecting velocity cancel by $u_x+w_z=0$. Thus $D\zeta =-{\sigma }_x$, orConservation of total buoyancy, with background gradient $N^2(z)$, givesBoth equations retain the nonlinear advection terms.

Write $\psi =\bar{\psi }(z)+\varphi$, with ${\bar{\psi }}_z=U(z)=\bar{u}(z)$, and take the background buoyancy anomaly to be zero. The linear equations for $(\varphi ,\sigma )$ areFor a normal mode $\varphi =\hat{\psi }(z)e^{ik(x-ct)}$, the second equation gives $\hat{\sigma }=N^2\hat{\psi }/(U-c)$. Substitute this into the first to obtain the Taylor–Goldstein equationHere $k\ne 0$ and the calculation is initially away from critical levels $U=c$. The first term is the restoring effect of stable buoyancy on vertically displaced parcels. The second comes from advection of the background vorticity gradient $U^{\prime \prime }$: a displacement creates a vorticity anomaly, whose induced velocity acts back on the displacement. This is the shear-flow counterpart of a Rossby wave restoring mechanism; its sign can also permit instability, so it is not invariably stabilizing.

Rigid boundaries imply $w=-{\varphi }_x=0$, hence $\hat{\psi }(0)=\hat{\psi }(\pi H)=0$. For a real, regular normal mode, multiply the Taylor–Goldstein equation by $\hat{\psi }$, integrate, and integrate its second derivative by parts. The endpoint term vanishes, leavingFor a smooth real dispersion branch, differentiate this identity with respect to $k$. The stated stationarity of $I$ removes the term involving $d\hat{\psi }/dk$, so $2k=I_{c}dc/dk$. Direct differentiation at fixed $\hat{\psi }$ givesSince the frequency is $kc$, the group velocity is $d(kc)/dk$, and thereforeThis is variational group velocity of a Taylor-Goldstein mode. Singular critical levels or a vanishing denominator require separate limiting analysis; no stationarity proof is needed here.

For the sinusoidal profiles, choose $c=0$. In the interior,The profile ratios extend to the endpoints in this reduced equation. Thus $n=1$ always gives a nonzero real $k$, with $k_1^2=N_0^2/U_0^2$. For this normal mode,The denominator in the group velocity expression is $2H(2N_0^2/U_0^3+1/(U_0{H}^2))$. It follows thatAlthough $U$ vanishes at the endpoints, this weighted integral is finite; the vanishing of ${\hat{\psi }}_1$ supplies the required cancellation. The formula is the limiting derivative approached from the regular real branch with $c<0$, before any interior critical level is encountered.

If instead $N=N_0$ and $U=U_0$ are constant, still take $c=0$. The same vertical eigenfunctions satisfy the Taylor–Goldstein equation, now withThere is no value of $n\ge 1$ guaranteed to give real $k$ for arbitrary positive parameters: if $N_{0}H<U_0$, all these squared wavenumbers are negative. A nonzero real $k$ exists exactly when $N_{0}H>U_0$; equality gives only the degenerate $n=1$, $k=0$ limit. For a real branch the same group velocity calculation gives $c_g=k_n^2{U}_0^3/N_0^2$.

To test exact nonlinear validity, take a real normal mode $\varphi =A\sin (nz/H)\cos (kx)$ and its corresponding buoyancy perturbation $\sigma =a(z)\varphi$, where $a=N^2/(U-c)$. Becausethe quadratic vorticity self-advection is zero: $J(\varphi ,{\nabla }^2\varphi )=0$. The quadratic buoyancy self-advection, however, isThis nonlinear exactness test for a stratified streamfunction mode distinguishes the cases. For constant $N,U$, $a=N_0^2/U_0$ is constant, so both quadratic terms vanish. The background-mode terms already satisfy the linear equations, hence the constant-profile modes are exact finite-amplitude solutions of the ideal Boussinesq equations. Pressure can be recovered from the momentum equations: the vorticity equation supplies their compatibility condition in this simply connected strip.

For the sinusoidal profiles, $a=(N_0^2/U_0)\sin (z/H)$ has a nonzero derivative in the interior. For every nonzero $A$ and nonzero real $k$, $a^{\prime }\varphi {\varphi }_x$ is not identically zero; it produces a second horizontal harmonic that is absent from the proposed solution. Thus none of these nontrivial sinusoidal-profile modes is an exact finite-amplitude solution. A $k=0$ field is merely an added horizontal shear with no vertical motion and is a degenerate exception to the wave calculation, not a propagating member of it.

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.

At $K{U}_0+\Lambda =0$, the free boundary Eady edge wave has zero frequency in the ground frame. If $kϵU_0\ne 0$, a steady amplitude cannot satisfy the forced amplitude equation. With $A(0)=0$, integration givesConsequentlyThe symbol multiplying $i$ here is the horizontal wavenumber $k$, as in the PDF. Taking the real part and differentiating in $x$ givesFor $U_0=-\Lambda /K>0$, positive $ϵ,f$, and $t>0$, the meridional-velocity minima coincide with topographic maxima, and its maxima with topographic minima. This is the lower panel. The change from quadrature to opposition is the phase of a resonantly forced topographic edge wave.

The stationary relief continuously excites a boundary Rossby wave whose intrinsic propagation exactly cancels mean advection. The forcing therefore remains coherent with that wave; with no damping or radiation loss for this trapped mode, its linear amplitude grows secularly. This is a resonant topographic quasi-geostrophic wave. Unlimited growth is a conclusion of the linear idealization: eventually the small-disturbance approximation fails.

The assertion that resonance must be unsteady needs nonzero forcing. If $k=0$, then $b_x=0$ and the boundary current does not cross the relief; if $U_0=0$, the forcing also vanishes. In either case the zero disturbance is a steady solution even if $K{U}_0+\Lambda =0$. For example, choose $k=0$, $\ell \ne 0$, and $\Lambda =-K{U}_0$; this is a concrete exception to an unconditional resonance claim. When both horizontal wavenumbers vanish, $K=0$ and the stated decaying-wave ansatz itself does not apply.

Interpret ${\rho }_2$ as the ambient density: the PDF's phrase assigning it a “depth” is a typographical error. Take $0<{\rho }_1-{\rho }_2\ll {\rho }_2$, and use a common reference density in inertia while retaining the density difference in buoyancy. This is the Boussinesq approximation, with reduced gravity $g^{\prime }=g({\rho }_1-{\rho }_2)/{\rho }_{\mathrm{ref}}>0$. For the shallow water approximation, the lower-layer depth must be small compared with longitudinal variation scales, vertical accelerations must be negligible and pressure hydrostatic. The deep ambient is approximated as a hydrostatic pressure reservoir away from the nose; assume a uniform velocity over each lower-layer cross-section and neglect drag and entrainment.

Write the prismatic triangular area as $A(h)=k{h}^2$ with fixed $k>0$, so the width at height $z$ is $2kz$. The excess hydrostatic pressure force divided by reference density isThe prismatic-channel shallow water equations are thereforeFor positive depth, divide the first equation by $2kh$ and use it in the second to obtainThe coefficient matrix in variables $(h,u)$ is $\left(\begin{array}{cc}u& h/2\\ g^{\prime }& u\end{array}\right)$. Its two distinct real characteristic speeds areThus the system is strictly hyperbolic for $h>0$ and $g^{\prime }>0$, though it degenerates at zero depth. Since $dc/dh=c/(2h)$, direct substitution in the two equations givesHence $R_{\pm }=u\pm 4c$ are the Riemann invariants conserved along their respective characteristic curves. The factor four, rather than the rectangular-channel factor two, comes from the triangular relation $A\propto h^2$.

The layer equations describe long-wave adjustment behind the nose, but do not resolve the short region where the dense fluid pushes aside the ambient. A finite-height front cannot be determined by simply adjoining a vacuum solution of the one-layer momentum equation. The ambient pressure and turning flow provide additional nose physics.

In the nose frame, the material condition makes the adjacent current velocity equal to the front speed in the laboratory frame. The two interior characteristic speeds relative to that front are then $\pm c_f$: one characteristic brings information from the body and one brings a boundary condition from the nose. This is why a gravity-current front condition is needed. With a deep ambient, fixed channel shape and negligible viscosity, the local buoyancy depth supplies the velocity scale $\sqrt{g^{\prime }{h}_f}$. A dimensionless closure therefore has the formwhere a constant $F$ is an idealization for a fixed geometry and chosen nose model. Dimensional analysis identifies the form but does not determine the coefficient.

To determine an appropriate coefficient theoretically, take a control volume moving with the nose and bounded by sections far enough ahead and behind for hydrostatic pressure and approximately uniform section velocities to apply. Use the actual triangular wetted areas and integrated pressures. Conservation of fluid volume fixes the ambient return-flow velocities; integrated horizontal momentum balances the incoming and outgoing momentum fluxes against the pressure forces. A Bernoulli relation on an ambient streamline, with an explicit nose-loss assumption if needed, supplies the pressure or energy closure. The density-anomaly pressure contribution of the lower layer contains the integral $g^{\prime }k{h}_f^3/3$ derived above. Solving those balances gives $u_f^2/(g^{\prime }{h}_f)$ and hence $F^2$.

Alternatively, measured front speeds and consistently defined adjacent heights in a geometrically similar experiment can calibrate the closure. A rectangular-channel numerical coefficient must not be transferred to a triangular channel without checking those balances. The height definition and dissipative nose structure also matter; a primary study of gravity-current front conditions explains why different height conventions can lead to different inferred coefficients. No numerical value of $F$ is required for the following solution.

The Batchelor entrainment hypothesis makes the inward entrainment velocity proportional to the local plume speed:where $\alpha >0$ is the dimensionless entrainment coefficient. For an upward plume $w>0$, this is $u_e=\alpha w$. The hypothesis represents turbulent ingestion of ambient fluid; it is a closure rather than a consequence of mass conservation.

Let $A=b^2$, suppressing the common factor $\pi$ from cross-sectional integrals. In the Boussinesq approximation, use a constant reference density ${\rho }_r$ in the leading-order mass conservation equation. Dividing that equation by ${\rho }_r$ gives volume conservation:The factor two is the ratio of the perimeter $2\pi b$ to the suppressed area factor $\pi$. Density differences are retained in buoyancy, even though they are neglected in the leading inertial and volume terms.

To derive the density-deficit balance, first leave the small density differences explicit in the supplied mass equation. Because the ambient density ${\rho }_0(z)$ is time independent, the product rule givesThus the ambient mass entrained through the perimeter cancels exactly against its contribution to the volume budget, yielding the required buoyancy balance. It is understood to the accuracy of the underlying Boussinesq approximation; the volume equation is not an additional exact constraint on a fully variable-density fluid.

For later use define reduced gravity $g^{\prime }=g({\rho }_0-\rho )/{\rho }_r$ and the squared buoyancy frequency $N^2=-g{\rho }_{0z}/{\rho }_r$. With $B=A{g}^{\prime }$ and $Q=Aw$, the balance isA stably stratified ambient has $N^2>0$, so upward motion consumes the plume's positive density deficit; an unstable ambient reverses this sign.

Use the conventional kinematic plume fluxesThen $\pi Q$ is the volume flux, $\pi {\rho }_{r}Q$ the dimensional mass flux, $\pi {\rho }_{r}M$ the dimensional momentum flux, and $\pi F$ the kinematic buoyancy flux. Calling $\pi Q$ a mass flux uses a density-normalized convention: the reference density has been suppressed. If density-weighted momentum and buoyancy fluxes are desired, multiply $M$ and $F$ by the same constant ${\rho }_r$ throughout. Here the conventional kinematic definitions make the Boussinesq approximation explicit.

For an upward plume with $Q,M>0$, the inverse relations areThe volume conservation, vertical momentum conservation, and buoyancy balances from part a becomeHere the vertical momentum conservation equation uses ${\rho }_r$ in the inertia and retains $({\rho }_0-\rho )gA={\rho }_{r}B$ in the buoyancy force. The ambient has no imposed vertical momentum at entrainment.

To classify the resulting partial differential equations, work in the nonsingular primitive variables $(A,w,g^{\prime })$ with $A>0$. Subtract $w$ times the volume equation from the momentum equation and $g^{\prime }$ times the volume equation from the buoyancy equation. This givesThe principal matrix is thereforeAll three characteristic speeds are real and equal to $w$, so the model is a hyperbolic system if that term means only reality of the characteristic roots. However, for $A>0$, $(K-wI)v=0$ requires $v_w=0$, leaving only two independent eigenvectors. HenceThis weak hyperbolicity of the unsteady top-hat plume prevents a complete set of independent characteristic fields. For example, freezing $A,w$ in the homogeneous principal part gives $D_t\delta w=0$ and $D_t\delta A=-A{\partial }_z\delta w$. A Fourier perturbation of wave number $k$ therefore produces a term proportional to $kt$ in $\delta A$. Real characteristic roots alone do not provide uniform control of high-frequency perturbations. The flux-variable transformation above is invertible in the stated region, so changing to $(Q,M,F)$ does not remove this degeneracy.

For a steady upward plume with prescribed $N^2=-G^2<0$, the kinematic plume fluxes satisfyThe positive sign in the last equation expresses the gain of buoyancy flux as the plume rises through the prescribed unstable ambient. Seek a positive power law solution $Q=q{z}^p$, $M=m{z}^r$, $F=f{z}^s$, with $z>0$. Matching powers givesThe last relation gives $s=p+1$ and the first $r=2p-2$. Substitution into the middle relation yields $4p-5=2p+1$, henceThe amplitude equations areFrom the first relation $m=9q^2/(4{\alpha }^2)$, while the last two give $m^2=G^2{q}^2/16$. Taking positive amplitudes gives $m=Gq/4$, and consequentlyReconstructing the radius, speed and reduced gravity gives the top-hat plume in constant unstable stratification:As a direct check, $M_z={\alpha }^2{G}^2{z}^3/9=FQ/M$ and $F_z={\alpha }^2{G}^3{z}^3/9=G^{2}Q$, while $Q_z={\alpha }^{2}G{z}^2/3=2\alpha \sqrt{M}$.

The corresponding ambient gradient is ${\rho }_{0z}={\rho }_r{G}^2/g$. Thus $\rho (z)={\rho }_0(z)-{\rho }_r{G}^{2}z/(4g)$; for a linear ambient this is ${\rho }_0(0)+3{\rho }_r{G}^{2}z/(4g)$. The plume is locally lighter than the ambient even though both densities increase upward in this prescribed unstable stratification.

The coordinate $z$ can be measured from a virtual origin. This exact power law branch has fluxes tending to zero there and describes buoyancy generated by the maintained unstable ambient; it does not impose arbitrary nonzero source fluxes at that origin. The given idealization holds the ambient quiescent despite its static instability. Finally, the Boussinesq approximation limits the modeled height range to small relative density changes, here $G^{2}z/g\ll 1$. Within these assumptions, the displayed positive branch proves the requested existence.

In the Boussinesq approximation, replace density by a constant reference density ${\rho }_0$ in inertia and incompressible flow, retaining the small density departure only in the gravitational force. Write $\rho =\bar{\rho }(z)+{\rho }^{\prime }$ and subtract the resting hydrostatic pressure. The upward buoyancy perturbation acceleration and squared buoyancy frequency areThus stable density stratification has $N^2>0$. Let $p$ henceforth denote the pressure anomaly divided by ${\rho }_0$. The ideal, nonrotating Boussinesq equations areThe buoyancy equation follows by materially conserving the total buoyancy $b_0(z)+\sigma$, where $b_0^{\prime }=N^2$; it does not assume constant $N$.

For an $x$-periodic flow of period $P$, define $\bar{a}=P^{-1}{\int }_0^{P}adx$. An infinite-domain average with vanishing end contributions works as well. The necessary pressure condition is no mean zonal pressure gradient, $\overline{p_x}=0$: periodic pressure satisfies it, whereas a pressure containing a term proportional to $x$ does not. Use incompressible flow to write the zonal momentum equation in conservative form,Averaging its zonal derivative terms yields the exact Reynolds stress relation

For two-dimensional small disturbances about rest, the Linearized Boussinesq equations are $u_t=-p_x$, $w_t=-p_z+\sigma$, ${\sigma }_t=-N^{2}w$, and $u_x+w_z=0$. Differentiate the two momentum equations to eliminate pressure: $(u_z-w_x{)}_t=-{\sigma }_x$. Differentiating in $x$ and using $u_{xz}=-w_{zz}$ gives ${\nabla }^2{w}_t={\sigma }_{xx}$. One further time derivative givesThe equation remains valid for $N=N(z)$. No vertical derivative of $N$ was taken: ${\partial }_x^2(N^{2}w)=N^2{w}_{xx}$.

For constant $N$ and $k\ne 0$, a plane internal gravity wave hasBoth signs give the prescribed lower-boundary velocity. The radiation condition excludes energy arriving from infinity. At positive frequency, the vertical group velocity isFor a source below the fluid, energy must propagate upward, so the outgoing solution has $m<0$, even though its vertical phase velocity $\omega /m$ is downward.

To obtain this selection causally, let $T=\mu t$ and $Z=\mu z$, and replace the constant amplitude by $a(Z,T)$. The method of multiple scales treats $(z,t,Z,T)$ as independent during differentiation: the carrier derivatives hold $Z,T$ fixed, while $a_T$ holds $Z$ fixed and $a_Z$ holds $T$ fixed. Only after differentiating do we restrict to $T=\mu t$, $Z=\mu z$. Acting on the complex wave, the operators becomeWith $K^2=k^2+m^2$, the wave equation has residualThe zeroth-order term is the dispersion relation; the first-order term vanishes precisely whenIf $a_0(T)$ is the switched-on boundary amplitude, with $a_0(T)=0$ for $T<0$, the outgoing envelope is $a(Z,T)=a_0(T-Z/b)$ for $b>0$. Its characteristics carry the boundary signal into initially undisturbed fluid. If $b<0$, characteristics instead arrive from the upper half-space; zero initial disturbances and no incoming signal cannot establish the proposed boundary-forced wavetrain. This is the internal-wave envelope radiation condition.

The leading incompressible flow polarization gives $u^{\prime }=-(m/k)w^{\prime }$. Hence, to second order in amplitude and leading order in the slow modulation,There is no meridional transport here. Envelope transport gives ${\partial }_t|a{|}^2=-b{\partial }_z|a{|}^2$, so an initially vanishing mean develops asThese are the leading expressions for internal-wave momentum deposition; higher slow-modulation corrections to the stress are beyond the displayed approximation. In particular, a spatially uniform fully established wavetrain has zero stress divergence, but the front has already produced its mean flow.

The mean wave energy density per unit mass follows from ${\sigma }^{\prime }=-i{N}^2{w}^{\prime }/\omega$:Thus $\bar{u}=(k/\omega )\mathcal{E}$. In the original frame the mean kinetic energy is ${\bar{u}}^2/2=O(a^4)$. In a frame translating at $c=\omega /k$, the undisturbed velocity is $-c$, and the change in mean kinetic energy is $-c\bar{u}+O(a^4)=O(a^2)$. Indeed $\mathcal{E}-c\bar{u}=0$ at this order behind the front. Mean-flow energy is therefore essential at the same order as wave energy in the translating frame. This follows also from the frame transformation of energy, $E\mapsto E-cP$ up to the fixed background constant, and is why conservation of momentum cannot be omitted from the moving-frame conservation of energy budget.

A precise stratification–rotation analogy uses constant $N$, ideal fluid, and independence of one horizontal coordinate. For the stratified system take ${\partial }_y=0$, $v=0$, and allow prescribed forces $(F_x,0,F_z)$. The possibly nonlinear Boussinesq equations areFor the unstratified system choose right-handed coordinates $(X,Y,Z)=(z,y,-x)$, rotation $2\mathbf{\Omega }=N{\mathbf{e}}_Z$, and ${\partial }_Y=0$. Its rotating Euler equations, with the static gravitational and centrifugal contributions absorbed into the background pressure, arewhere $D_R={\partial }_t+U{\partial }_X+W{\partial }_Z$. The identification isThen $D_R=D$, $P_X=p_z$, and $P_Z=-p_x$, proving the component-by-component equivalence. In particular, buoyancy perturbation maps to transverse rotating velocity, and stratified horizontal velocity maps to rotation-axis velocity. The mapping permits arbitrary prescribed forces in the two-dimensional plane, but no transverse force: a transverse force in the rotating problem would create a source in the stratified buoyancy equation. A transverse stratified force would also invalidate $v=0$. Constant $N$ is essential to this constant-rotation equivalence.

A physical free surface introduces extra structure beyond the interior equations. After subtracting hydrostatic pressure, its zero total pressure condition gives $p=g\eta$ to first order, together with the kinematic condition ${\eta }_t=w$ for the stratified vertical displacement. The mapping takes this normal velocity into $U$, whereas an actual rotating free surface normal to its physical vertical $Z$ has normal velocity $W$ and a displacement along $Z$. The image of the stratified horizontal surface is a vertical surface in the rotating coordinates. Consequently the physical free-surface problems are not analogous, although deliberately mapped rigid boundaries or mathematical boundary data can be.

For the first prescribed force, the steady linear stratified equations are satisfied byReal parts are understood. The horizontal pressure gradient balances the force; the vertical pressure gradient is balanced by buoyancy; and $w=0$ satisfies buoyancy conservation and incompressible flow. Thus neither $u(z)$ nor $b(z)$ is fixed by the steady equations alone.

Initial rest and switching on in the unbounded fluid remove this ambiguity by causal pressure adjustment by internal waves. Work with the $e^{imz}$ coefficient and let $\hat{F}(q)$ be the horizontal Fourier transform of $f^{\prime }(x)$. For an abrupt switch at time zero, elimination in the linear equations givesAt $q=0$ the bracket cancels the apparent singularity at every finite time. The persistent pressure term is the symmetric antiderivative of the force: the inverse transform of $-i\hat{F}/q$, interpreted as a Cauchy principal value, is $f(x)-\epsilon /2$. The oscillatory part propagates away. To see why its zero-frequency part does not leave an arbitrary offset, note that $\Omega (q)\sim cq$ with $c=N/m$; its long-wave contribution consists of oppositely moving fronts. At fixed $x$, the two far-field values of $f-\epsilon /2$ are $-\epsilon /2$ and $+\epsilon /2$, whose average is zero. Nonzero horizontal wavenumbers give dispersive transients that decay locally. Smooth finite-duration switching changes the transient, not this selected steady limit. ThereforeA spatially constant pressure gauge is immaterial. The long-wave group velocity $c=N/m$ gives the adjustment scale: the elapsed time since switching must satisfy $Nt\gg 1$ and $(N/m)t\gg |x|+x_1$, and exceed the switching duration. These are asymptotic conditions, not an assertion of exact equilibration after a finite time; the full dispersion relation allows tails. The same long-wave limit gives $u_{\infty }=\epsilon m{e}^{imz}/(2N)$ if desired. Equal-amplitude outward fronts leave opposite pressure offsets and the same velocity on their two sides.