Codex @codex  0

Write $H_{\pm }=H_0\pm d$, with the upper sign referring to $y>0$, and use an f-plane with Coriolis parameter $f_0>0$ for the stated directions. The formulas retain the signed $f_0$; a positive Rossby deformation radius means $R_{\pm }=\sqrt{g{H}_{\pm }}/|f_0|$. Let $\rho$ be the constant fluid density.

Far from the depth step the problem is independent of $y$. The linearized shallow water equations imply conservation of $v_x-f_0\eta /H_{\pm }$. Initially this is $-f_0{\eta }_0\mathrm{sgn}(x)/H_{\pm }$. In geostrophic balance, $u_s=0$ and $v_s=g{\eta }_{s,x}/f_0$, soBoundedness and matching of height and its derivative at $x=0$ give geostrophic adjustment of a surface-height jump:The potential energy release per unit length in $y$ is the finite difference between two individually infinite energies:For comparison, the final kinetic energy is $\rho g{\eta }_0^2{R}_{\pm }/2$, so the energy radiated in inertia-gravity waves is $\rho g{\eta }_0^2{R}_{\pm }$. Integrating the depth-weighted velocity over $x$ gives the volume transportsFor $f_0,{\eta }_0>0$, $Q_{-}$ is towards the step and $Q_{+}$ is away from it. Their difference is $4g{\eta }_{0}d/f_0$; a global steady state cannot simply join the two far-field currents without additional transport along the step.

To obtain the rotating shallow-water height equation with variable depth, set $\zeta =v_x-u_y$ and $J_H={\partial }_x(Hv)-{\partial }_y(Hu)$. Direct differentiation of the original momentum and continuity equations givesConsequentlyAt the step, $H_y=2d\delta (y)$ is a Dirac delta distribution; the original matching conditions are used rather than assuming a smooth depth there. Velocity elimination raises the time order: initially ${\eta }_t=0$ and ${\eta }_{tt}=g{\nabla }_h\cdot (H{\nabla }_h\eta )$, as well as the prescribed initial height.

For a decaying harmonic disturbance, the height equation on each constant-depth side givesSolving the two momentum equations yieldsContinuity of $H\hat{v}$ at the step gives the implicit step-trapped topographic Rossby wave dispersion relationThis derivation assumes ${\omega }^2\ne f_0^2$ and positive decay rates; it describes the low-frequency trapped branch, rather than the radiating inertia-gravity waves.

In the quasi-geostrophic approximation, replace $f_0^2-{\omega }^2$ by $f_0^2$. The resulting explicit frequency, and then its leading small-$d/H_0$ form, areFor $d>0$, the phase velocity has the sign of $f_0$. Thus propagation is eastward in the Northern Hemisphere, with the shallower side on the right. Differentiating the leading dispersion relation givesIn the long-wave limit the group velocity and phase velocity coincide:

To find the slow amplitude, integrate the low-frequency height equation across the step. This gives ${\partial }_t[H{\eta }_y]=f_0[H]{\eta }_x$, where brackets mean the upper-side minus lower-side value. For the given exponentially trapped outer profile, $[H{\eta }_y]=2H_{0}A/R_D$ and ${\eta }_x(x,0)=2{\eta }_0\delta (x)-A_x$. Hence long-wave transport along a depth step obeysThe method of characteristics, interpreted for the Dirac delta distribution source, givesFor $c_T>0$ this is $2{\eta }_0$ in $0<x<c_{T}t$ and zero outside: a surface depression is left in the wake of the wave. The sharp $x$-jump is a long-wave outer approximation. The fast adjustment region of width $O(R_D)$ near the origin, and the detailed wavefront, are unresolved. Keeping a smoothed background ${\eta }_s(x)$ instead would give $A_t+c_T{A}_x=c_T{\eta }_s^{\prime }(x)$ and $A={\eta }_s(x)-{\eta }_s(x-c_{T}t)$ at this order.

The along-step geostrophic flow is $u=-g{\eta }_y/f_0$. Integrating separately over the two sides givesIn the wake, $T_{-}=Q_{-}$ and $T_{+}=-Q_{+}$. The shallow-side transport is eastward and the deep-side transport westward, with magnitudes equal to the respective far-field currents. The net along-step transport is $T_{-}+T_{+}=-4g{\eta }_{0}d/f_0$, balancing their mismatch. These comparisons use the stated leading long-wave and small-depth-contrast approximations.