Codex @codex  0

Use pressure per unit constant density, $\pi$, with the gravitational potential and centrifugal potential absorbed into it. Strict geostrophic balance is the exact balancewith material acceleration omitted, for a homogeneous incompressible flow. Writing $f=2\Omega \ne 0$, its components are $fv={\pi }_x$, $fu=-{\pi }_y$ and ${\pi }_z=0$. Differentiating the first two equations with respect to $z$ gives $u_z=v_z=0$. Alternatively, the vector calculus identityshows, on taking the curl of the balance, that $(\mathbf{\Omega }\cdot \nabla )\mathbf{u}=0$. Thus $w_z=0$ too, proving the Taylor–Proudman theorem:The physical pressure still includes its background hydrostatic pressure; it is the reduced pressure $\pi$ that is independent of $z$. Between horizontal impermeable boundaries, the vertically constant $w$ is zero.

For boundaries $z=z_b(x,y)$ and $z=z_t(x,y)$, impermeability condition gives $w_b={\mathbf{u}}_H\cdot {\nabla }_H{z}_b$ and $w_t={\mathbf{u}}_H\cdot {\nabla }_H{z}_t$. A nearly columnar flow crossing contours of $H=z_t-z_b$ must therefore have $w_t-w_b\simeq {\mathbf{u}}_H\cdot {\nabla }_{H}H\ne 0$. Exact columnarity must give way to small corrections from ageostrophic flow. The leading horizontal geostrophic flow can remain nearly independent of $z$, and the leading reduced pressure has ${\pi }_z=0$, but $w_z$ is now nonzero and exact three-dimensional strict geostrophic balance no longer holds. For small slopes, the dominant vertical vorticity equation iswhere $q=v_x-u_y$ is relative vorticity. Equivalently, $D[(q+f)/H]/Dt\simeq 0$. Crossing geostrophic contours stretches or squashes columns: increasing $H$ produces cyclonic relative vorticity, decreasing $H$ produces anticyclonic relative vorticity. The full vorticity equation also contains tilting terms, neglected at this columnar order.

For the circular shallowing, fluid arriving from upstream has potential vorticity $f/h$. Its height is $h(1-ϵ)$ inside, so conservation givesThis uses the stated column-stretching description and initially open streamlines. To obtain the leading disturbance velocity field, use the conventional streamfunction definition $u=-{\Psi }_y$, $v={\Psi }_x$, and write $\Psi =-Uy+\psi (r)$. Then $q={\nabla }_H^2\psi =r^{-1}(r{\psi }_r{)}_r$. Regularity at $r=0$ gives ${\psi }_r=-fϵr/2$ inside. Outside, ${\psi }_r=B/r$; matching tangential velocity at $r=a$ fixes $B=-fϵa^2/2$. Choose the irrelevant additive constant so thatBoth $\psi$ and ${\psi }_r$ match at $a$, so there is no extra vorticity sheet on the rim. This is the circular topographic anticyclone. The radial expression specifies the disturbance; the complete streamfunction is $\Psi =-Uy+\psi$, since a radial function alone cannot represent the required uniform far-field current. Reversing the streamfunction sign convention reverses the displayed proportionality constant, but leaves the physical velocity unchanged.

Differentiation gives the full relative velocity field in the rotating frame:away from the narrow transition across the change of depth. The disturbance turns clockwise for $fϵ>0$ and decays as $1/r$. The horizontal field is the leading small-$ϵ$ columnar inversion: the exact variable-depth transport across a sharp rim needs higher-order corrections, rather than an exactly two-dimensional velocity continuous across a literal vertical step.

Take $U>0$ and $ϵ>0$. The greatest opposing zonal disturbance speed is $Ca$, attained at $(0,-a)$. If $Ca<U$, then $u>0$ everywhere and every particle advances in $x$, precluding closed streamlines. If $Ca>U$, the interior stagnation point $(0,-U/C)$ is a center, since the interior contours are circles about that point. The exterior stagnation point $(0,-C{a}^2/U)$ is a saddle: linearization there has real eigenvalues of opposite sign. Closed contours around the center are bounded by a separatrix through the saddle. The two points meet at the rim at onset, giving the closed-streamline threshold for a circular topographic anticyclone:At this threshold $U/(fa)=ϵ/2\ll 1$, consistent with small Rossby number. Equality is the onset, not a finite-area closed cell. Beyond onset, upstream data alone no longer determine the potential vorticity of trapped fluid; the threshold follows from the continuation of the supplied open-flow model.

The basic quasi-geostrophic streamfunction $\bar{\psi }=-\Lambda yz$ gives $\bar{u}=\Lambda z$, $\bar{v}=0$ and constant quasi-geostrophic potential vorticity $\bar{Q}=f$. DefineBecause the background has no potential-vorticity gradient, the linear potential-vorticity equation is simplyIts characteristics have constant $y,z,x-\Lambda zt$, so the general solution is $q^{\prime }=\mathcal{F}(y,z,x-\Lambda zt)$. In particular, $q^{\prime }=F(y,z)G(x-\Lambda zt)$ is the requested arbitrary separated family; one product is not the most general initial field.

Linearizing the vertical velocity requires retaining advection of ${\bar{\psi }}_z=-\Lambda y$ by $v^{\prime }={\psi }_x^{\prime }$. ThereforeThe impermeable boundary at $z=0$ imposesThis is conservation of boundary buoyancy, not merely a condition that ${\psi }_z^{\prime }$ vanish.

For zero interior potential vorticity, write ${\psi }^{\prime }=\mathrm{Re}\{\varphi (z,t)e^{i(kx+ly)}\}$. With $K_h=\sqrt{k^2+l^2}>0$, stratified quasi-geostrophic inversion givesThe decaying solution isNo time dependence has been assumed to obtain this vertical structure. Applying the boundary condition gives $-\mu \dot{A}-i\Lambda kA=0$, henceThus the Eady edge wave travels eastward with horizontal phase velocity $\omega /k=\Lambda /\mu$ for $k\ne 0$.

The propagation mechanism can be made explicit using boundary parcels. Thermal wind gives the background buoyancy gradient ${\bar{b}}_y=f{\bar{\psi }}_{zy}=-f\Lambda$. A small meridional parcel displacement $d$ produces boundary buoyancy perturbation $b^{\prime }=-d{\bar{b}}_y=f\Lambda d$. Meanwhile the decaying potential-vorticity inversion gives $b^{\prime }=f{\psi }_z^{\prime }=-f\mu {\psi }^{\prime }$ at the boundary. Consequently ${\psi }^{\prime }=-\Lambda d/\mu$ there, and the parcel equation $d_t=v^{\prime }={\psi }_x^{\prime }$ becomesA displaced buoyancy contour therefore induces a velocity field that advances that displacement eastward. No interior potential-vorticity gradient is needed because the boundary buoyancy supplies the invertibility data.

Vortex stretching is what sets the induced flow's vertical reach. In the interior, the zero-potential vorticity condition requiresBoundary isopycnal displacement decreases with height as $e^{-\mu z}$, so fluid columns are stretched or compressed over depth ${\mu }^{-1}=f/(N{K}_h)$. Conservation of potential vorticity makes this stretching balance the relative vorticity. The resulting inversion factor $b^{\prime }=-f\mu {\psi }^{\prime }$ gives the propagation speed: stronger density stratification reduces penetration and speed, while larger $f$ increases them. Zero interior potential vorticity therefore does not imply zero interior flow.

For imposed boundary vertical velocity, retain the same decaying structure and writeThe boundary relation becomes the forced Eady edge wave equationThe integrating factor $e^{i\omega t}$ gives, for ${\omega }_0\ne \omega$,The arbitrary homogeneous coefficient $B$ sets the initial disturbance. With $A(0)=0$, this is $A=C(e^{-i{\omega }_{0}t}-e^{-i\omega t})/[i(\omega -{\omega }_0)]$.

At resonance, ${\omega }_0=\omega$, integration instead yieldsAn initially undisturbed fluid has ${\psi }^{\prime }=Ct{e}^{-\mu z}\cos (kx+ly-\omega t)$. The imposed pumping reinforces the same boundary normal mode every cycle, with no dissipative loss, so the linear amplitude grows secularly. The divergence describes resonance of the undamped linear model; sufficiently large wave amplitude invalidates the assumed small-disturbance approximation.

Relative vertical vorticity is the horizontal Laplacian of the streamfunction, not its full stretched three-dimensional Laplacian. Inside the spherical potential-vorticity anomaly, the quadratic solution givesOutside, in the equatorial plane $Z=0$, let $s=(x^2+y^2{)}^{1/2}>R$ and $C=Q_0{R}^3/3$. The azimuthal velocity is $v_{\theta }={\psi }_s=C/s^2$, soThe opposite signs thus follow directly from solid-body rotation inside and the exterior velocity's rapid radial decay. The azimuthal velocity itself has the same rotational sense on both sides; its radial derivative makes the exterior relative vorticity negative.

There is no contradiction with the zero exterior PV anomaly. At $Z=0$ the exterior solution has ${\psi }_{ZZ}=C/s^3$, so $q_g+{\psi }_{ZZ}=0$. The stratification/stretching contribution compensates the horizontal relative vorticity. Inside it contributes $Q_0/3$, supplementing $2Q_0/3$ to give the prescribed excess $Q_0$.

The potential-vorticity gradients in a meridional two-layer current are important here: besides ${\bar{q}}_{n,y}=\beta$, one has ${\bar{q}}_{1,x}=-FV$ and ${\bar{q}}_{2,x}=FV$. They must be retained when linearizing, even though the basic relative vorticity vanishes.

Let $K^2=k^2+l^2$ and $A=K^2+F$. For a normal mode the two-layer quasi-geostrophic potential vorticity amplitudes areBecause the basic velocities are northward $(V,0)$ in the two layers, their advection frequencies are $\sigma -lV$ and $\sigma$. The uniform forcing has no perturbation, $W^{\prime }=0$. The linear equations are thereforeIn particular, a perturbation's zonal velocity advects the basic interface-induced zonal potential-vorticity gradient. Substitution of the plane wave givesA nonzero disturbance exists exactly when the determinant vanishes. The requested relation for linear stability of a meridional two-layer current isEquivalently, with $\mathcal{B}=K^2(K^2+2F)$,Its discriminant isFor $K>0$, exponential baroclinic instability occurs when this discriminant is negative; otherwise the two frequencies are real. This also displays the stabilizing contribution of the planetary vorticity gradient when $k\ne 0$. As a check, $V=0$ gives the uncoupled barotropic mode and baroclinic mode frequencies $-\beta k/K^2$ and $-\beta k/(K^2+2F)$.

The vertical component $q_z=\zeta -f\eta /H$ is $H$ times the first-order anomaly in shallow-water potential vorticity about rest:Thus $\mathbf{q}=q_z\hat{\mathbf{z}}$ records the combination of relative vorticity and column stretching that is conserved by the unforced linearized shallow water equations. Indeed their curl and continuity equation give $(q_z{)}_t=0$. It is the persistent forcing of the balanced component during geostrophic adjustment, while inertia-gravity waves carry away the unbalanced component.

Initially $q_z=-f{\eta }_0\mathrm{sgn}(x)/H$. In the final geostrophic flow, $u=0$, $v=g{\eta }_x/f$, and $\zeta =g{\eta }_{xx}/f$. Conservation of the linearized shallow-water potential vorticity therefore givesBoundedness on both sides and continuity of $\eta ,{\eta }_x$ at the origin determine the constants:The solution is the locally attained balanced state after outgoing inertia-gravity waves have departed; the full inviscid system need not lose its globally conserved wave energy. The signed volume flux in the $y$ direction is

For double Kelvin-wave adjustment at a depth step, put $c_{\pm }=\sqrt{g{H}^{\pm }}$ and use $R_D^{\pm }=c_{\pm }/f$ on the corresponding side of the step. The PDF's far-field wording must mean $\eta \to {\eta }_0\mathrm{sgn}(x)$, rather than ${\eta }_0$ for every $x$. In the exponentially trapped outer response, cross-step geostrophic balance gives, at $y=0$,To retain the propagation of this slowly varying trapped disturbance, the along-step acceleration in $u_t-fv=-g{\eta }_x$ must also be retained. It supplies the small cross-step ageostrophic flow. Thus the two cross-step volume fluxes areHere $d\mathrm{sgn}(x)/dx=2\delta (x)$ as a distributional derivative. Equating the fluxes and dividing by $c_{+}+c_{-}$ givesThe printed forcing is missing a factor of two. It cannot follow from the printed sign-function surface profile with the usual Dirac delta distribution normalization. Pure geostrophic balance in both horizontal directions would also omit $A_t$ and could not determine the wave evolution.

The method of characteristics with zero initial amplitude gives the corrected weak solutionBehind the eastward front, $0<x<\Delta ct$, the surface elevation at the step becomes $-{\eta }_0$, matching the upstream side. All the original cross-step transport is diverted into flow along the step once the wave has passed a fixed location: $v=0$ there, with opposing along-step velocities above and below the step. This is the complete barrier in the adjusted region, and in the long-time limit at every fixed $x$. At finite time a front still carries cross-step flux: the formulas above give $H^{\pm }{v}^{\pm }=(2{\eta }_0{c}_{+}{c}_{-}/f)\delta (x-\Delta ct)$. The discontinuous outer representation does not resolve the short-scale fast adjustment at either front.

For completeness, solving the equation literally as printed instead gives $A={\eta }_0[\Theta (x)-\Theta (x-\Delta ct)]$. Its amplitude is half that required; $\eta (x,0,t)$ retains a jump ${\eta }_0$ at the origin after passage, so that literal equation does not imply complete blocking. The corrected equation and the literal printed equation are distinguished rather than asserting an inconsistent barrier result.

Assume inviscid, unforced evolution, so shallow-water potential vorticity is conserved. If a column has initial depth $H_1$ and relative vorticity ${\zeta }_1$, moving to $H_2>H_1$ at fixed Coriolis parameter givesIts absolute vorticity is amplified by vortex stretching. A column initially without relative vorticity develops cyclonic spin: positive in the Northern Hemisphere and negative in the Southern Hemisphere. More generally the increase is in the direction of its original absolute vorticity, not necessarily the direction of $f$ if the initial spin is unusually strong.

Conservation of shallow-water potential vorticity with fixed relative vorticity givesFor ordinary planetary-vorticity-dominated columns, moving towards the equator decreases the magnitude of absolute vorticity, so the column must become shallower and spread horizontally. This is the compression counterpart of vortex stretching. For $\zeta =0$, the depth ratio is $f_2/f_1$, in either hemisphere before reaching the equator. The displayed general formula is essential if the fixed relative vorticity is large: the sign of the depth change then depends on the original absolute vorticity, and conservation cannot continue through zero absolute vorticity with positive finite depth and nonzero conserved potential vorticity.

Use a one-layer, hydrostatic, homogeneous fluid with reference depth $H_0$, no forcing or friction, small Rossby number $ϵ=U/(|f_0|L)\ll 1$, slow evolution, and leading geostrophic balance. For the beta plane, take $\beta L/|f_0|=O(ϵ)$; the relative free-surface displacement is also small, $\eta /H_0=O(ϵ)$, with the Rossby deformation radius retained at the chosen horizontal scale. This excludes the equatorial limit. If a bottom elevation $h_b$ is included, additionally require $h_b/H_0=O(ϵ)$ and stationary topography. Set $h_b=0$ for a flat bottom.

The quasi-geostrophic streamfunction satisfies $\psi =g\eta /f_0$, $u=-{\psi }_y$, $v={\psi }_x$, so relative vorticity is $\zeta ={\nabla }_h^2\psi$. Expanding the shallow-water potential vorticity using actual depth $H_0+\eta -h_b$ givesConsequently the depth-scaled anomaly, the shallow-water quasi-geostrophic potential vorticity, isThe leading material derivative uses the geostrophic flow. Define $J(a,b)=a_x{b}_y-a_y{b}_x$; then conservation yields $Q_t+J(\psi ,Q)=0$. The term $-\psi /R_D^2$ represents free-surface stretching, $\beta y$ the variation of planetary rotation, and $f_0{h}_b/H_0$ topographic stretching. Under a rigid lid or at scales much smaller than the Rossby deformation radius, the free-surface term is absent or negligible and the flat-bottom Barotropic Rossby wave model has $Q={\nabla }_h^2\psi +\beta y$.

Let $Q_b=f/H$. At sufficiently small Rossby number, retain advection of this full background shallow-water potential vorticity while dropping the higher-order advection of $\zeta /H$. For uniform-depth Sverdrup balance, this also requires relative-vorticity advection to be small against planetary-vorticity advection; small Rossby number alone does not guarantee this when $\beta$ is very small. Multiplying the resulting steady potential vorticity equation by $H$ givesDefine the topographic potential-vorticity pseudovelocityIts sign is important: it is opposite to the coefficient vector on the left of the preceding equation. Substituting the ocean transport streamfunction expression for relative vorticity yieldsThe pseudovelocity follows contours of $f/H$ and is not the fluid velocity; its units are inverse length squared per time because $\psi$ is a transport streamfunction. If $r=F=0$, the ocean transport streamfunction is constant along each connected background potential vorticity contour. Locally, at regular points,with arbitrary differentiable $\Phi$. Different disconnected components of a level set may carry different functions until boundary conditions identify them. If $f/H$ is constant over an open region, the leading equation gives no restriction there; critical points similarly require regularity and global matching rather than division by a vanishing gradient.

Let $\mathbf{\zeta }=\nabla \times \mathbf{u}$ and $\mathbf{\omega }=\mathbf{\zeta }+2\Omega {\mathbf{e}}_z$ be the relative vorticity and absolute vorticity. The nonlinear acceleration identity rewrites the momentum equation asTaking the curl removes the pressure, kinetic energy gradient and shearing-sheet tidal potential. ThusBoth $\nabla \cdot \mathbf{u}$ and $\nabla \cdot \mathbf{\omega }$ vanish. The vorticity equation is consequentlyFor the absolute-vorticity flux tensor $T_{ij}=u_i{\omega }_j-{\omega }_i{u}_j$, use the flux convention $(\nabla \cdot \mathbf{T}{)}_j={\partial }_i{T}_{ij}$, with the first index denoting transport direction. This givesThe scalar $\theta$ is a buoyancy displacement variable, so dimensional consistency gives it units of length; the printed potential-temperature terminology presupposes this normalization. Its curl source directly forces only the $y,z$ components of absolute vorticity. Hence radial vorticity conservation in a shearing sheet takes the local flux formIts volume integral is constant when the net boundary flux vanishes. However $D_t{\omega }_x=\mathbf{\omega }\cdot \nabla u_x$ still permits vortex stretching; the radial component is not generally a materially conserved scalar in three dimensions. Using divergence on the second index of the stated tensor would reverse the transport sign, so the tensor convention matters.

For the Keplerian shearing sheet background $\mathbf{U}=-3\Omega x{\mathbf{e}}_y/2$, the only spatial derivative is ${\partial }_x{U}_y=-3\Omega /2$. Since $U_x=0$ and all fields are independent of $y$, $\mathbf{U}\cdot \nabla \mathbf{U}=0$ and $\nabla \cdot \mathbf{U}=0$.

The Coriolis acceleration exactly balances the tidal term:With constant pressure and $\theta =0$, the buoyancy displacement variable equation also has zero advection, zero source $U_x$ and zero diffusion. Thereforeis a steady solution. Its absolute vorticity is ${\mathbf{\omega }}_0=(\Omega /2){\mathbf{e}}_z$, since the background relative vorticity is $-3\Omega {\mathbf{e}}_z/2$.

The vertical relative vorticity isDifferentiate the $v$-momentum equation with respect to $x$ and the $u$-momentum equation with respect to $y$, then subtract. The mixed pressure derivatives cancel, givingThe mass-conservation equation says ${\eta }_t=-h(u_x+v_y)$. Therefore the linearized shallow-water potential vorticitysatisfiesThus

The basic state has velocity fieldand buoyancy $B=N^{2}z$. For disturbances independent of $y$, the linearized equations areSubstituting a plane wave proportional to $\exp [i(kx+mz-\omega t)]$ and eliminating $p$, $v$, and $b$ gives the dispersion relationThus the requested coefficients are

An instability exists precisely when some wavenumber pair makes ${\omega }^2<0$. Since the Brunt–Väisälä frequency satisfies $N^2>0$, this is possible exactly whenThe basic relative vorticity is $\Lambda {\mathbf{e}}_z$, so its absolute vorticity is $(f+\Lambda ){\mathbf{e}}_z$. Its Ertel potential vorticity is thereforeThe instability criterion can consequently be written as $fQ<0$: the vertical absolute vorticity has the opposite sign to the planetary vorticity. This is inertial instability, approached most directly by disturbances with $|m/k|\gg 1$.

Expand both the forcing and response in the orthogonal vertical normal modes. Projection onto $P_n$ givesIf the horizontal forcing scale obeys $L_h\gg c_n/f_0$, the relative vorticity term ${\nabla }_h^2{\varphi }_n$ is small compared with the stretching term. The resulting long-wave equation isThus each mode communicates the forcing westward at its long Rossby wave speed $C_n$, with the barotropic mode fastest because $c_0$ is largest.

Potential vorticity combines vorticity with vortex stretching and density stratification. It also applies to nonrotating stratified flow through Ertel potential vorticity. In the inviscid shallow water equations,where $\zeta$ is relative vorticity and $h$ is layer thickness.

The ocean transport streamfunction convention $Hu=-{\psi }_y$, $Hv={\psi }_x$ gives relative vorticity $\zeta =H^{-1}{\nabla }^2\psi -H^{-2}\nabla H\cdot \nabla \psi$. This identity is exact for variable positive depth and does not assume that depth variations are small.

A uniform positive quasi-geostrophic potential vorticity anomaly in a sphere of stretched radius $R$ has the displayed regular, decaying streamfunction. Both it and its radial derivative match at the sphere. In its equatorial plane the horizontal relative vorticity is $2Q_0/3$ inside and $-Q_0{R}^3/(3r^3)$ outside. The exterior potential-vorticity anomaly remains zero because its vertical stretching contribution compensates that relative vorticity. A stretched sphere corresponds to a flattened ellipsoid in physical coordinates when $N\gg |f|$.