Codex @codex  0

For an initially resting shallow water layer whose height rises linearly from $-h$ to $h$ across $(-L,L)$, let $R_d=\sqrt{gH}/|f|$ and $\alpha =L/R_d$. Conserved linear shallow-water potential vorticity gives $R_d^2{\eta }_s^{\prime \prime }-{\eta }_s=-{\eta }_0$. Boundedness, oddness, and continuity of height and slope give ${\eta }_s/h=x/L-e^{-\alpha }\sinh (x/R_d)/\alpha$ inside the ramp, and ${\eta }_s/h=\mathrm{sgn}(x)[1-(\sinh \alpha /\alpha )e^{-|x|/R_d}]$ outside. The balanced velocity is $u_s=0$, $v_s=g{\eta }_s^{\prime }/f$. Narrow ramps spread over the Rossby deformation radius; wide ramps retain their initial height except for narrow matching layers.

A constant-depth layer initially at rest with a step in its free surface adjusts by radiating inertia-gravity waves. Conservation of linearized shallow-water potential vorticity leaves the bounded balanced profile ${\eta }_s={\eta }_0\mathrm{sgn}(x)(1-e^{-|x|/R_D})$, where $R_D=\sqrt{gH}/|f|$ is the Rossby deformation radius. The released potential energy per unit transverse length is $3\rho g{\eta }_0^2{R}_D/2$; the final kinetic energy is $\rho g{\eta }_0^2{R}_D/2$, leaving $\rho g{\eta }_0^2{R}_D$ in radiated waves.

In nearly columnar flow between rigid boundaries, a geostrophic contour is a contour of column height divided by the Coriolis parameter. If relative vorticity is negligible, conservation of shallow-water potential vorticity constrains leading geostrophic flow to these contours. Motion across them stretches or squashes columns, generating relative vorticity; for constant $f$, the contours are simply contours of $H$.

Let $H=h_{00}+\zeta$ and write the horizontal shallow water equations on the constant f-plane asKeeping only terms linear in the disturbance about rest gives the linearized shallow water equationsDefine the horizontal divergence $\delta =u_x+v_y$ and relative vorticity $q=v_x-u_y$. Taking respectively the divergence and vertical curl of the momentum equation yieldsTogether with ${\zeta }_t=-h_{00}\delta$, this provesThe nonlinear shallow-water potential vorticity is $Q=(f+q)/H$. The nonlinear vertical vorticity equation gives $D(f+q)/Dt=-(f+q)\delta$, while mass conservation gives $DH/Dt=-H\delta$; their ratio is therefore materially conserved. Expanding about rest,Advection of its constant background contributes nothing at first order, explaining why the displayed linear anomaly is conserved at each fixed position.

A further time derivative eliminates $q$ and $\zeta$ from the divergence equation:For a plane wave with horizontal wavevector $(k,l)$, write $K=\sqrt{k^2+l^2}$. More directly, substituting the plane wave into the three linear equations gives a coefficient matrix with determinant proportional to $\omega [{\omega }^2-f^2-c_0^2{K}^2]$. Thus the full linear rotating shallow-water dispersion relation has one balanced branch and two wave branches:The latter are inertia-gravity waves. The positive branch starts at $f$ with zero slope and tends to the straight line $c_{0}K$; the negative branch is its reflection. The requested sketch, including the stationary branch, is

For the positive-frequency branch, the magnitudes of the phase velocity and group velocity areThe group velocity vector is $c_0^2(k,l)/\omega$; the negative branch reverses its direction. Near $\omega =f$, namely $K{R}_D\ll 1$ with the Rossby deformation radius $R_D=c_0/f$, $c_p\sim f/K$ is large while $c_g\sim c_0^{2}K/f$ is small. The $K=0$ limit is an inertial oscillation, for which a spatial phase speed is not defined. Far above $f$, both speeds approach $c_0$, the shallow water gravity-wave speed, with $c_p>c_0$ and $c_g<c_0$ at finite $K$.

On the zero-frequency branch, geostrophic balance givesFor a nonconstant plane wave, $\hat{q}=-(g{K}^2/f)\hat{\zeta }$: relative vorticity and elevation have opposite signs. They can both be nonzero and spatially structured, whereas the horizontal divergence vanishes identically. A constant elevation is also a zero-frequency solution, with zero velocity and zero relative vorticity.

For the initially resting ridge, set ${\zeta }_0(x)=ϵh_{00}{e}^{-x^2/L^2}$ and assume an unbounded domain, uniform in $y$. Initial potential vorticity fixesIn the balanced remainder, insert $q_s=(g/f){\zeta }_s^{\prime \prime }$ to obtainThe decaying Green's function of this modified Helmholtz equation gives the Geostrophic adjustment of a Gaussian height ridge explicitly:The initial outward pressure force creates a cross-ridge current; the Coriolis acceleration turns that current and builds the along-ridge geostrophic flow. The remainder is radiated inertia-gravity waves, not a dissipative disappearance of energy.

For precision, use the Fourier transform $\hat{\zeta }(k)=\int \zeta (x)e^{-ikx}dx$. Initial rest implies ${\zeta }_t(x,0)=0$, and the conserved anomaly givesConsequentlyThese formulas satisfy both initial conditions. They also show that the waves carry no linear potential vorticity anomaly: $q=(f/h_{00})(\zeta -{\zeta }_0)$ and $\delta =-{\zeta }_t/h_{00}$ reconstruct the remaining fields, with $u$ and $v$ obtained from their derivatives. On the infinite line, dispersive wave packets propagate away and decay locally, leaving the balanced ridge and its oppositely directed flanking currents. Conservation of energy still holds globally. In a reflecting finite basin, standing waves may persist, so a permanent pointwise relaxation is not guaranteed without dissipation.

If $L\gg R_D$, most of the initial elevation remains balanced: ${\zeta }_s\simeq {\zeta }_0$ with small smoothing corrections. If $L\ll R_D$, the balanced ridge spreads to width $R_D$, and away from the original narrow coreThus the final central height is much smaller than the initial one and the wave component is substantial. The integral of the balanced elevation equals the initial integral, since the Green's function has unit integral; the oscillatory component has zero net added volume.

Potential-vorticity inversion is the diagnostic recovery of the flow from its advected vortical scalar, using the balance relation and specified boundary conditions. Its usefulness is a separation of tasks: transport the scalar by the current velocity, invert to find the new velocity, and repeat. The inverse is spatially nonlocal even though the transport law is local. The four models differ chiefly in the operator being inverted and the information needed at boundaries.

In two-dimensional vortex dynamics, planar incompressibility permits $\mathbf{u}=(-{\psi }_y,{\psi }_x)$, and vertical vorticity is $\zeta ={\nabla }^2\psi$. The curl of inviscid momentum givesThere is no vortex stretching in strictly planar motion, so vorticity is materially conserved. The two-dimensional potential-vorticity inversion is the Poisson equation ${\nabla }^2\psi =\zeta$. On the plane, with suitable integrability and far-field behavior, its Green's function is $(2\pi {)}^{-1}\log |\mathbf{x}-{\mathbf{x}}^{\prime }|$, givingDifferentiation gives the planar Biot-Savart kernel. Here ${\psi }_h$ is a harmonic function fixed by boundary, circulation and far-field data. A rigid impermeable wall is a streamline, so $\psi$ is constant on it; constants or circulations around distinct holes require specification. The vorticity alone cannot distinguish two flows differing by a permitted harmonic contribution, such as a uniform current on an unbounded plane. A constant added to $\psi$ has no physical effect.

In nonrotating layerwise-two-dimensional vortex dynamics, strong stable density stratification constrains motion to almost horizontal layers. In the ideal independent-layer limit, vertical velocity vanishes at leading order, $(u,v)=(-{\psi }_y,{\psi }_x)$ andhold separately at each fixed $z$. Height is a parameter: the horizontal inverse Laplacian contains no vertical derivative. With background buoyancy $N^{2}z$, the leading Ertel potential vorticity is $Q=\mathbf{\omega }\cdot \nabla b\simeq N^2\zeta$, so the materially advected normalized scalar is $Q/N^2\simeq \zeta$. The same horizontal Poisson equation is inverted layer by layer, using each layer's boundary data. This is an asymptotic constrained model, not arbitrary three-dimensional Euler motion with vertical velocity set to zero: its validity requires small Froude number and controlled vertical scales. At finite buoyancy frequency, the weak vertical motion and tilted isopycnals restore coupling. The inversion may be independent by layer while the reconstructed fields vary with height. The horizontal Euler equations for an inviscid fluid determine the pressure gradient, while hydrostatic balance relates its vertical derivative to buoyancy.

In shallow-water quasi-geostrophic inversion, begin with the exact shallow-water potential vorticity $(f+\zeta )/h$. Linearizing its anomaly about depth $H$ gives, after multiplying by $H$, $\zeta -f_0\eta /H$. For small Rossby number, geostrophic balance gives $\psi =g\eta /f_0$, while $\zeta ={\nabla }_h^2\psi$. Keeping the beta plane background gradient givesThus $q-\beta y$ is inverted with the modified Helmholtz equation $({\nabla }_h^2-L_D^{-2})\psi =q-\beta y$, under the specified boundary and far-field conditions. The velocity follows from $\psi$ and surface elevation from $\eta =f_0\psi /g$. A horizontal Fourier mode of the anomaly satisfiesThe Rossby deformation radius therefore limits the reach of the balanced response. Layer stretching, represented by $-\psi /L_D^2$, is absent from ordinary planar vorticity inversion. Unlike the pure Poisson equation, the finite-deformation operator also responds to the constant part of $\psi$, which specifies mean elevation rather than an arbitrary streamfunction offset. Specified bottom topography adds a known source to $q$; it is subtracted along with the background before inversion.

In stratified quasi-geostrophic inversion, continuously varying vertical displacement replaces the shallow layer's single stretching term. For positive $N^2(z)$ and constant nonzero $f_0$, the three-dimensional quasi-geostrophic potential vorticity isGeostrophic balance and hydrostatic balance give velocity $(-{\psi }_y,{\psi }_x)$ and buoyancy perturbation $b=f_0{\psi }_z$. Invert the anisotropic elliptic boundary value problem using interior $q$, prescribed top and bottom buoyancy as Neumann boundary conditions, and the lateral, far-field and global-flow data. Boundary buoyancy evolves by its own material conservation when the boundary is impermeable. It cannot be discarded: the Eady edge wave in the preceding solution has zero interior potential vorticity anomaly but nonzero boundary buoyancy and a nonzero induced flow.

For constant $N,f_0$, the stretched coordinates for quasi-geostrophic inversion $Z=Nz/|f_0|$ turn the inversion operator into the three-dimensional Laplacian. Its infinite-space Green's function is $-1/(4\pi |\mathbf{X}-{\mathbf{X}}^{\prime }|)$ in $(x,y,Z)$ coordinates. Hence a horizontal disturbance of scale $L$ communicates vertically over scale $|f_0|L/N$, explicitly coupling layers. In a bounded domain with homogeneous difference data, multiplying the difference of two inversions by its streamfunction and applying integration by parts givesThis proves uniqueness up to a constant when the relevant boundary terms vanish; a mean or pressure reference fixes the remaining constant. Compatible boundary fluxes are needed for a purely Neumann boundary-value problem.

Across all four models, inversion provides the induced flow, while scalar transport supplies time evolution. The independent horizontal inverse in the first two models becomes a finite-deformation inverse in shallow water and a vertically coupled elliptic inverse in continuous stratification. The principle applies within the stated balance or constrained model. In the unrestricted Boussinesq equations or shallow water equations, gravity waves and divergent motion can share the same potential vorticity, so the vortical scalar alone does not reconstruct every possible flow.

With the free-surface elevation denoted by $\zeta$, the rotating shallow water equations over a fixed bottom areThe momentum gradient is that of the surface elevation, not depth alone; bottom slope enters through the depth and continuity equation. Taking ${\partial }_x$ of meridional momentum minus ${\partial }_y$ of zonal momentum gives the vorticity equationThe terms quadratic in velocity gradients combine into ${\omega }_r{\nabla }_H\cdot {\mathbf{u}}_H$; the constant Coriolis parameter contributes $f{\nabla }_H\cdot {\mathbf{u}}_H$. Continuity gives $D_{H}h/Dt=-h{\nabla }_H\cdot {\mathbf{u}}_H$. Dividing the two equations proves the exact shallow-water potential vorticity invariant

For the quasi-geostrophic approximation, let $U_s$ and $\ell$ be representative velocity and horizontal length scales. Require the Rossby number $\mathrm{Ro}=U_s/(|f|\ell )\ll 1$, a shallow layer, and slow evolution on time $\ell /U_s$. Put $L_R=\sqrt{g{h}_{00}}/|f|$ and assume $\mathrm{Bu}=(L_R/\ell {)}^2$ is order one. Leading momentum then gives geostrophic balance, and a quasi-geostrophic streamfunction can be chosen asIts surface-height scale satisfies $|\zeta |/h_{00}=O(\mathrm{Ro}/\mathrm{Bu})\ll 1$. Take the bottom-height ratio of the same small order, with gentle slopes, so $\eta =\zeta -b$ obeys $|\eta |/h_{00}\ll 1$. Relative vorticity is order $\mathrm{Ro}f$, and expanding the exact potential vorticity givesBecause the leading geostrophic velocity is horizontally nondivergent, ageostrophic advection of the anomaly is one order smaller. At the first evolving order the shallow-water quasi-geostrophic potential vorticity equation isThe constant reference $f/h_{00}$ has zero gradient and drops out. These assumptions also require all topographically induced velocities to retain a small Rossby number; the small bottom-height condition alone is not a full justification of quasi-geostrophic balance.

The uniform eastward current has ${\psi }_0=-Uy$ up to an additive constant, so ${\zeta }_0=-(fU/g)y$. With $b_0=-(fU/g)y$, surface and bottom slopes are equal and $h=h_{00}$. Its relative vorticity vanishes and $Q_0=0$. This sloping background is used locally where its surface and bottom departures remain within the small-height regime.

Add the compact ridge and set $\psi =-Uy+\varphi (x)$. Then $u=U$, $v={\varphi }^{\prime }(x)$, and $Q$ depends only on $x$. The steady equation becomes $U{Q}_x=0$. Uniform incoming potential vorticity fixes that constant to zero, giving the cosine-ridge shallow-water geostrophic responseLet $a=\pi L/2$, $\lambda =L_R^{-1}$ andThe forcing is even, so the decaying response is even. A cosine particular solution inside the ridge leaves an even homogeneous term; outside, retain only the decaying exponential. Since the forcing is bounded and contains no delta function, both $\varphi$ and ${\varphi }^{\prime }$ are continuous at the edges. Write $\varphi =A\cos (x/L)+D\cosh (x/L_R)$ inside and $\varphi =B{e}^{-(|x|-a)/L_R}$ outside. Matching at $x=a$ givesSolving yields $D=A(L_R/L)e^{-a/L_R}$, henceA globally decaying homogeneous solution that is continuously differentiable is zero, so this matched solution is unique.

A steady streamline is a constant-$\psi$ curve. If its upstream ordinate is $y_{\infty }$, it therefore hasFor $f>0$, $\varphi$ is positive: the flow first turns toward increasing $y$, passes the crest at its largest deflection, then returns to its upstream ordinate. Reversing the sign of $f$ reverses the deflection. For $L=3L_R$, the center deflection is $A[1+e^{-3\pi /2}/3]/U$ and each tail decays on scale $L_R$. The figure shows several streamlines and the ridge support; the deflection extends beyond the ridge.

One compatible realization of the initial condition has $h(x)=h_{00}-\hat{b}(x)$ and $\zeta =b_0(y)$. Its stated zonal velocity then satisfies $u(x)h(x)=U{h}_{00}$, so the initial mass flux is constant and $v=0$. The data given only far upstream do not uniquely determine the initial height everywhere; this is the natural mass-flux-compatible choice. It is not fully balanced over the ridge: $fu(x)$ differs from $-g{\zeta }_y=fU$. The resulting transverse acceleration initiates geostrophic adjustment and launches fast inertia-gravity waves.

Initially the relative vorticity is zero, but exact shallow-water potential vorticity is $f/[h_{00}-\hat{b}(x)]$, not the uniform upstream value. That localized anomaly is materially transported downstream, approximately at $U$ in the leading quasi-geostrophic approximation; it is not annihilated by adjustment. Fresh incoming parcels of uniform PV subsequently supply the ridge region, whose balanced potential-vorticity inversion approaches the steady solution above. The fast waves can carry unbalanced energy away on the inertia-gravity time scale, short compared with the advective time at small Rossby number. In an open domain this can leave the local balanced response; an ideal closed or reflecting domain need not settle monotonically, because undamped waves can return. Radiation, advection of the initial PV anomaly, and any weak damping must be distinguished from irreversible loss of PV.

Write $c=\sqrt{gH}$ and $R_d=c/|f|$, with $f\ne 0$. For the profiles below take $f>0$, so $\alpha =L/R_d$ agrees with the stated parameter; changing the sign of $f$ reverses the balanced meridional velocity but not the height profile. The linearized shallow water equations areTaking the horizontal curl gives $(v_x-u_y{)}_t=-f(u_x+v_y)=f{\eta }_t/H$. Hence the signed linear shallow-water potential vorticity anomaly is conserved:The initial data are independent of $y$, so the solution remains independent of $y$ and $v_x=f(\eta -{\eta }_0)/H$. Differentiating continuity in time and using the zonal momentum equation gives

For a horizontal Fourier mode with wavenumber $k$, put ${\Omega }_k^2=f^2+c^2{k}^2$. The solution splits into a balanced component and inertia-gravity waves:These waves have frequencies $|\omega |\ge |f|$, phase velocity $\omega /k$, and group velocity $c^{2}k/\omega$, whose magnitude is less than $c$. They have zero perturbation potential vorticity and carry the unbalanced part away. On the unbounded line, the localized difference between the initial and balanced profiles disperses, leaving a local steady limit. No damping is needed for this local geostrophic adjustment; in a closed periodic domain, undamped waves would persist and a literal pointwise steady limit would generally not exist.

The bounded steady height satisfies the modified Helmholtz equationOddness removes the even homogeneous solution in the central region. The exterior solution is bounded and tends to the appropriate plateau. Matching both height and its first derivative at $x=\pm L$ gives the geostrophic adjustment of a finite-width height ramp:For example, writing the interior solution as $hx/L+B\sinh (x/R_d)$, matching its slope to the decaying exterior at $L$ gives $B=-h{e}^{-\alpha }/\alpha$.

The steady meridional momentum equation gives $u_s=0$, and the zonal equation gives the geostrophic balance $f{v}_s=g{\eta }_s^{\prime }$. ThusThe height is odd and monotone, and the meridional current is even, with maximum $v_s(0)=gh(1-e^{-\alpha })/(fL)$ for $h>0$.

For $\alpha \ll 1$, the initial ramp is narrow compared with the Rossby deformation radius. To leading order,Near the origin the height rises with slope $h/R_d$, rather than $h/L$; the current has width $R_d=L/\alpha$ and peak $gh\alpha /(fL)$. For $\alpha \gg 1$, the central height remains nearly $hx/L$ and the central velocity nearly $gh/(fL)$. Around each end of the ramp there is a matching layer of width $R_d=L/\alpha$, with height correction of order $h/\alpha$; the exterior current decays over the same width. At $x=\pm L$ the current is approximately half its central value. In both limits the zonal current is zero. These are the broad rotational adjustment of a narrow disturbance and the weak edge adjustment of a broad disturbance, respectively.

Use the energy normalization specified for the height energy: the linear total energy density is $(u^2+v^2)/2+g{\eta }^2/(2H)$. The linear equations implyFor $\alpha \ll 1$, let the potential energy loss mean the positive difference between the initial and final values, per unit transverse length. Although both absolute potential energies diverge on the infinite line, their difference is finite. The leading step approximation givesThe finite ramp changes this leading expression by a relative $O(\alpha )$ correction. The initial kinetic energy is zero and the final gain isThe remainder $\Delta V-\Delta T\simeq g{h}^2{R}_d/H$ is energy carried away by the outgoing inertia-gravity waves. This explains why the ratio is less than one without invoking dissipation: conservation of energy includes the waves, not only the final balanced flow.

Full equations and the balanced scaling. Let $H=h_{00}$, free-surface displacement $\eta$ and layer thickness $h=H+\eta -b$. The inviscid rotating shallow water equations areHere gradients are horizontal and $b$ is a fixed bottom height. For horizontal scales $L$, velocities $V$ and time scale $L/V$, acceleration is $V^2/L$, whereas the Coriolis term is $|f|V$. A small Rossby number $\mathrm{Ro}=V/(|f|L)$ therefore gives leading geostrophic balance. DefineThe last length is the Rossby deformation radius. Geostrophy gives $\eta \sim |f|VL/g$ and consequentlyThe usual distinguished quasi-geostrophic limit takes the Burger number $L_R^2/L^2$ of order one, small surface displacement, and bottom amplitude $b/H=O(\mathrm{Ro})$. The ageostrophic velocity is smaller than the geostrophic velocity by a factor of order $\mathrm{Ro}$, but its small divergence supplies the first nonzero thickness tendency. Motions are slow compared with $|f{|}^{-1}$, and the shallow-layer approximation also requires depth small compared with horizontal scale. Other deformation-radius limits are possible only if the small-height condition is retained.

Taking the vertical curl and combining it with thickness continuity gives exact shallow-water potential vorticity conservation,Expand to first order in relative vorticity, displacement and relief; multiplication by $H$ removes the irrelevant constant depth factor. To quasi-geostrophic accuracy,Potential-vorticity conservation transports this scalar with the geostrophic flow. Potential-vorticity inversion solves the elliptic problemwith the required boundary or far-field conditions, then differentiates $\psi$ to recover velocity and surface height. The deformation term suppresses the far-field influence over distances comparable with $L_R$; at scales much shorter than $L_R$ the inversion approaches the two-dimensional Poisson problem.

The full shallow-water system evolves two velocity components and layer thickness independently, with three first-order time equations. The quasi-geostrophic system evolves one scalar; balanced velocity and height are diagnostic consequences of inversion. In particular it filters the two fast inertia–gravity branches, whose resting-layer dispersion is ${\omega }^2=f^2+gH|\mathbf{k}{|}^2$, and does not describe their fast geostrophic adjustment. Reducing the time order restricts admissible initial data; it does not turn an arbitrary full shallow-water state into a quasi-geostrophic state automatically.

The intended Gaussian inversion, and its required background qualification. Set $r^2=x^2+y^2$ and ${\psi }^{\prime }=A{e}^{-r^2/a^2}$. Direct differentiation givesConsequently, in an inversion with uniform PV and its background gradient explicitly compensated, the printed relief is balanced byThis is the Gaussian topographic response with compensated uniform potential vorticity. In that model the perturbation PV is identically zero, and any steady velocity field advects this uniform PV trivially.

There is a genuine uniform-current obstruction to constant shallow-water QG potential vorticity in the literal constant-$f$, finite-$L_R$ formulation. A geostrophic background current $(U,0)$ has total streamfunction $\Psi =-Uy+{\psi }^{\prime }$ and a background free-surface slope. Over flat bottom its PV is $\overline{Q}=f+Uy/L_R^2$, not uniform. For the displayed Gaussian and printed relief, the two disturbance terms cancel in inversion but leaveFor $Uϵ\ne 0$ and finite $L_R$ this residual is nonzero; for example, it is nonzero at $x=a/\sqrt{2},y=0$. Thus the claimed Gaussian is not a steady solution of the standard unforced constant-$f$ equations as literally stated. Omitting the mean contribution to the stretching term is not a legitimate gauge choice.

A precise intended repair is to supply a background PV gradient $-Uy/L_R^2$ from an explicitly compensating weak bottom slope or beta-plane background, so the current really does have uniform PV. Another repair is the rigid-lid limit $L_R\to \infty$. In the uncompensated finite-$L_R$ current model, a Gaussian can instead solve steady equations if the $a^2/L_R^2$ term is removed from the relief: then ${\nabla }^2{\psi }^{\prime }+fb/H=0$ and $Q=f-\Psi /L_R^2$, so $J(\Psi ,Q)=0$. Its PV is not uniform. These repairs distinguish the intended calculation from the false literal claim.

Smallness and closed streamlines of the intended response. Its maximum disturbance speed isIts relative vorticity is of order $|fϵ|/H$, and its surface displacement has maximum $|{\eta }^{\prime }|=|ϵ|a^2/L_R^2$. The relief satisfies $|b|\le |ϵ|(4+a^2/L_R^2)$. Sufficient conditions for quasi-geostrophic consistency are thereforeThe second condition limits the background geostrophic height change across the disturbance scale. The last controls disturbance Rossby number, surface height and bottom relief, including positive layer depth. For $a$ comparable with $L_R$ it reduces to $|ϵ|\ll H$. It does not require disturbance velocity to be small relative to $U$.

Take $U>0$ without loss of horizontal orientation. The total streamfunction is $\Psi =-Uy+A{e}^{-r^2/a^2}$, withIf $\sqrt{2/e}|A|/a<U$, then $u>0$ everywhere and no streamline can close. In terms of the relief amplitude,is the strict no-closed-cell condition. Above this closed-streamline threshold for a Gaussian disturbance in uniform flow, stagnation points occur at $x=0$ on the side opposing the current. The equation $|y|e^{-y^2/a^2}=|U|a^2/(2|A|)$ has two roots: the inner one is a center because both Hessian eigenvalues of $\Psi$ have the same sign; the outer one is a saddle. Closed contours surround the center and are bounded by a separatrix. Equality gives a degenerate onset, with no finite-area closed cell yet. If $U=0$, every nontrivial Gaussian response instead has circular closed streamlines.

Closed streamlines are compatible with small Rossby number: their threshold can be reached with $|ϵ|/H$ of order the already small background Rossby number. They do not alone invalidate quasi-geostrophy. They do, however, disconnect the cell from upstream PV data. If only the upstream PV is specified, the PV in such a cell depends on its formation history and cannot generally be assigned from the upstream value alone. Globally uniform initial PV would remain uniform in the ideal compensated model, but mixing, dissipation and forcing history can select a different physical trapped circulation. These qualifications are separate from the constant-$f$ inconsistency demonstrated above.

Choose a fixed vertical datum so the physical free surface is at $H+\zeta$, with $H$ constant; its depth above the bottom is $h=H+\zeta -b$. The gradient of the constant datum does not enter the pressure force. With $D_H={\partial }_t+u{\partial }_x+v{\partial }_y$, the rotating shallow water equations areThe fluid is hydrostatic with depth-independent horizontal velocity at this order. The pressure-gradient force uses the free-surface slope, not the layer-depth slope: $\nabla \zeta =\nabla h+\nabla b$.

Let $q=v_x-u_y$ be relative vertical vorticity, $q_a=f+q$ the absolute vorticity, and $\delta =u_x+v_y$. Taking the horizontal curl of the two momentum equations and expanding the advective derivatives givesMass conservation is $D_{H}h=-h\delta$. The two stretching factors cancel, proving the exact shallow-water potential vorticity conservation lawThe denominator is $h$, the actual column thickness that changes under horizontal convergence; an arbitrary surface elevation $\zeta$ does not obey that column mass equation. There is no $w{\partial }_z$ term because the reduced layer quantities have no vertical dependence. The vertical motion has already supplied the layer-thickness equation, rather than remaining an independent coordinate of this model.

Now take $H=h_{00}$ and assume $|q|\ll |f|$, $|\zeta |,|b|\ll H$. To first order,For a small Rossby number, leading geostrophic flow satisfies $fv=g{\zeta }_x$ and $fu=-g{\zeta }_y$. ThusFor $f\ne 0$, the barotropic deformation radius is $L_R=\sqrt{gH}/|f|$. Multiplication of the first-order PV by $H$ givesThis is the required shallow-water quasi-geostrophic potential vorticity, including the bottom contribution.

Write $S=fϵ/H$ and $L=L_R$. The initial state has $\psi =0$, so ${\mathcal{Q}}_i=f+(f/H)b$. The specified homogenization sets ${\mathcal{Q}}_f=f$ in $|y|<a$ and leaves it unchanged outside. Hence the requested change isFor $f>0$ this is a descending straight segment from $+Sa$ to $-Sa$, with zero exterior values and jumps at the two edges. Endpoint values do not affect the inversion. It is the anomaly, not the full final PV profile, that appears as the source below.

For the zonal mean, $v=0$, $u=-{\psi }^{\prime }(y)$, and potential-vorticity inversion requireswith decaying disturbance at both infinities. The odd source makes the unique decaying solution odd. A particular interior solution is $S{L}^{2}y$, so writeBoth $\psi$ and ${\psi }^{\prime }$ must be continuous at $y=\pm a$: a jump would introduce a delta function or its derivative into a source that has only finite jumps. Matching at $a$ givesand thereforeDifferentiate the interior streamfunction to obtainThe exterior profiles needed for the sketches are $u=(D/L)e^{-(|y|-a)/L}$ and $\zeta =(f/g)\mathrm{sgn}(y)D{e}^{-(|y|-a)/L}$.

For $L\gg a$, put $A=a/L$, $Y=y/L$, so $|Y|\le A$. Uniform Taylor expansion givesConsequentlyFor $f>0$ the leading interior current is westward and approximately parabolic.