Codex @codex  0

With constant density divided out, a decaying or periodic geostrophic balance on an f-plane haswhere $q=\zeta -f\eta /H$ is the signed linearized shallow-water potential-vorticity anomaly. The stationary height equation and integration by parts prove the identity. Replacing $q$ by $|q|$ destroys it: an odd height paired with symmetric opposite-sign sheets would then misleadingly give zero.

Apply mass conservation to a fixed horizontal region $D$. The layer mass is $\rho {\int }_{D}hdxdy$ and the outward horizontal flux is $\rho {\int }_{\partial D}h(u,v)\cdot \mathbf{n}ds$. The divergence theorem and arbitrary choice of $D$ giveThis depth-integrated law retains free-surface motion; neglecting the vertical velocity in horizontal momentum is not the assumption that the depth cannot change.

Take the free-surface pressure to be constant. hydrostatic pressure gives $p=p_{\mathrm{atm}}+\rho g(h-z)$, so the horizontal pressure forces are $-g{\eta }_x,-g{\eta }_y$. Linearization about a resting layer gives the linearized shallow water equationsTaking the curl of the momentum equations yields ${\zeta }_t=-f(u_x+v_y)$. Combining with continuity gives the linearized shallow-water potential-vorticity anomalyFor $Q_0=0$, $\zeta =f\eta /h_0$. Taking the divergence of momentum gives ${\partial }_t(u_x+v_y)-f\zeta =-g{\nabla }_h^2\eta$. Eliminate the divergence with continuity to obtainFor $\theta =k(x-ct)$, substitution of $\eta ={\eta }_0\cos \theta$ gives the inertia-gravity wave dispersion relationThe associated purely oscillatory velocity field isThese formulas satisfy all three linearized equations and the zero-anomaly condition. A separately added spatially uniform inertial oscillation is not part of this monochromatic wave.

To first order in amplitude evaluate the wave at a particle's equilibrium coordinate $x_{\ast }$, so ${\theta }_{\ast }=k{x}_{\ast }-\omega t$. Integrating its velocity gives displacements relative to the orbit centreThese are particle ellipses for rotating shallow-water waves, with axis ratio $V/U=f/\omega <1$. For the stipulated positive signs, the velocity at phases $0,\pi /2,\pi ,3\pi /2$ points respectively right, up, left, down, with magnitudes $U,V,U,V$. The particle orbits are clockwise as time advances, since the phase decreases.

Take $\mathbf{f}=f\hat{\mathbf{z}}$ with $f>0$, and put $\zeta =v_x-u_y$, $d=u_x+v_y$ and $q=\zeta -f\eta /H$. Here $q$ denotes the signed vertical component of the linearized shallow-water potential-vorticity anomaly, rather than its magnitude. The divergence and vertical curl of the linearized shallow water equations giveConsequently $q_t=0$. Eliminating $d$ and substituting $\zeta =q+f\eta /H$ givesThe vector forcing is $-H\mathbf{f}\cdot \mathbf{q}$. This is potential-vorticity conservation in its linear, f-plane form: the conserved anomaly forces a stationary balanced part, while the homogeneous equation supports inertial-gravity waves.

For the initial strip, differentiating the discontinuous velocity in the sense of distributions givesThese are two oppositely signed vortex sheets. In the final geostrophic balance, $v_f=0$ and $u_f=-g{\eta }_f^{\prime }/f$. The height therefore solvesThis is the Rossby deformation radius. Decay at infinity fixes the Green function to $G(y)=-\lambda e^{-|y|/\lambda }/2$. Thus the geostrophic adjustment of a finite-width current has the particularly useful representationExpanding the exponentials inside the strip gives ${\eta }_f/H=-(U/c)e^{-a/\lambda }\sinh (y/\lambda )$; above the strip it gives $-(U/c)e^{-y/\lambda }\sinh (a/\lambda )$, and below it gives $(U/c)e^{y/\lambda }\sinh (a/\lambda )$. The height is continuous, odd and exponentially localized. Its derivative has the jumps required by the two Dirac delta functions.

Differentiating the height, rather than assuming a uniform final current, gives the complete velocity:The one-sided velocity jump is $U$ at each edge, with opposite orientations. The value exactly on an idealized vortex sheet is immaterial. Inside the strip the current remains in the original direction; outside it a return current develops.

Let $r=a/\lambda$ and $s=y/\lambda$. For $r\ll 1$, $u_f/U=1-r+O(r^2)$ over the narrow strip, while the outside return flow is approximately $-r{e}^{-|s|}$. The height varies almost linearly across the strip, $c{\eta }_f/(HU)\simeq -s$, and has extrema of magnitude $(1-e^{-2r})/2\simeq r$ at its edges. For $r\gg 1$, the central current is exponentially small: $u_f(0)/U=e^{-r}$. Each edge supports a layer of width $\lambda$, with $u_f/U\simeq \frac{1}{2}{e}^{-(a-|y|)/\lambda }$ on its inner side and $-\frac{1}{2}{e}^{-(|y|-a)/\lambda }$ outside. The corresponding height extrema approach $\pm HU/(2c)$, with almost zero height deep inside and far outside. Both requested profiles are drawn from the exact functions below; the dashed jumps represent one-sided velocity limits.

For balanced shallow-water energy as a signed potential-vorticity pairing, divide out the constant density and work per unit distance in $x$. The kinetic energy plus surface potential energy isMultiplication of the stationary height equation by ${\eta }_f$ and integration by parts, with vanishing end terms, yieldsThe same expression holds with $dA$ for a finite two-dimensional domain or a finite periodic length in $x$. The infinite strip's total energy is infinite, so both energies and their ratio are understood per unit $x$-length.

There is a sign defect in the printed energy expression: $|\mathbf{q}|$ must be replaced by the signed component $\hat{\mathbf{f}}\cdot \mathbf{q}$. Indeed, $|q_0|=|U|[\delta (y-a)+\delta (y+a)]$, so pairing it with the odd ${\eta }_f$ gives zero, although the balanced state has positive energy. The signed pairing givesSince $E_i=H{U}^{2}a$, the energy retention in finite-width geostrophic adjustment isFor $r\ll 1$ this ratio is $1-r+\frac{2}{3}{r}^2+O(r^3)$: nearly all energy remains in the narrow balanced current. For $r\gg 1$ it is approximately $1/(2r)$: only the edge regions retain balanced energy. Conservation of energy still holds for the inviscid evolution. The missing balanced energy is carried away by inertia-gravity waves; “final state” means the local balanced limit after those waves leave, rather than a loss of the total energy over the entire infinite domain.