Codex @codex  0

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.

Use the meridional overturning streamfunction conventionDifferentiating geostrophic and hydrostatic balance gives the thermal wind relation $f_0{\overline{u}}_z=(g/{\rho }_0){\overline{\rho }}_y$. Differentiate it in time and use the mean momentum and density equations. The resulting diagnostic elliptic boundary value problem isFor the transformed circulation it becomes the Eliassen equation for residual circulation, with signs appropriate to the convention above:At the rigid meridional walls, impermeability requires each streamfunction to be constant along each wall. Set those constants to zero, choosing the solution with no imposed net vertical throughflow. Since the prescribed density flux vanishes at both walls, the same conditions apply to ${\mathcal{X}}^{\ast }$.

There is a genuine data limitation: prescribing $F$ alone does not prescribe $R$. The term $-f_0{R}_{yz}$ remains in both circulation equations. To give the usual explicit step-flux circulation, assume vanishing meridional momentum-flux divergence, $R_y=0$, with no independently imposed circulation. The formulas below are conditional on that assumption, not a unique consequence of the density flux alone.

Let $l=\pi /L$, $\kappa =N_{0}l/|f_0|$ and $A_0=-1/S>0$, using the flux unit in which its lower plateau is $-1$. The bounded solution has $\mathcal{X}=\sin (ly)\chi (z)$. Its vertical equation isAt vertical infinity require bounded Eulerian flow and vanishing residual flow. Requiring the Eulerian vertical velocity to vanish also at lower infinity would be inconsistent with the continuing lower-region flux divergence. At $z=0$, continuity of $\chi$ and ${\chi }^{\prime }$ follows because the forcing has a step but no delta function. The step-flux residual circulation in a stratified channel is thereforeHere ${\mathcal{X}}^{\ast }=\sin (ly){\chi }^{\ast }$, because $F/S=A_0\sin (ly)$ below zero and vanishes above it. The residual streamfunction decays at both infinities, whereas the Eulerian streamfunction tends to $-A_0\sin (ly)$ below and to zero above.

For $z\ne 0$, the Eulerian velocity isThe flow descends near the $y=0$ wall and rises near $y=L$, with a negative meridional return flow concentrated near $z=0$. Below the dissipation level its vertical flow balances the continuing divergence of eddy density transport. There is no independently added throughflow, but opposing vertical transports at lower infinity are required by this idealized forcing.

Away from the discontinuity the residual meridional velocity equals the Eulerian value, whileThus the residual mean circulation forms two oppositely rotating cells, confined within vertical distance $O(1/\kappa )$ of the termination of the wave flux. Its streamfunction jumps by $-A_0\sin (ly)$ at zero. Distributionally this supplies the horizontal return transportThe Dirac delta sheet is a consequence of the discontinuous imposed flux, not an extra boundary condition. Smoothing the flux termination replaces it by a thin finite return current. Streamlines must not be drawn as continuously crossing the jump without that transport.

The physical Eulerian mean flow is continuous. In the transformed momentum equation, the sheet contribution to $f_0{v}^{\ast }$ cancels the sheet in the flux divergence, leaving the smooth mean accelerationThus termination of the vertical Eliassen–Palm flux produces a finite column-integrated wave drag, redistributed over depth by the residual mean circulation.

For completeness, the general bounded response can be calculated without setting $R_y$ to zero. Expand $R_{yz}$ in wall-compatible sine modes, with coefficients $B_n(z)$, $l_n=n\pi /L$ and ${\kappa }_n=N_0{l}_n/|f_0|$. Add to either displayed streamfunction the same correctionfor decaying momentum-flux forcing. The Green function sign follows from $({\partial }_z^2-{\kappa }_n^2)e^{-{\kappa }_n|z-z^{\prime }|}=-2{\kappa }_n\delta (z-z^{\prime })$.

An explicit counterexample to uniqueness is the additional momentum flux $R=ϵ[1-\cos (2ly)]e^{-a|z|}$, with $a>0$. It vanishes at both walls and leaves the prescribed density flux unchanged. For $a\ne 2\kappa$ it addsto both circulations. This nonzero bounded, wall-impermeable solution verifies that the original mean equations need an additional momentum-flux specification to select a unique sketch.

Interpret strict ellipticity as uniform ellipticity up to the boundary, as required for the global elliptic boundary value problem. Mere pointwise positive definiteness in the interior is insufficient: $L=x^2{D}^2$ on $(0,1)$ has trivial homogeneous Dirichlet kernel but $Lu=1$ has no $C^2$ solution up to $x=0$.

Set $X=\{u\in C^{2,\alpha }(\overline{\Omega }):u{|}_{\partial \Omega }=0\}$ and $Y=C^{0,\alpha }(\overline{\Omega })$. The Fredholm alternative for an elliptic Dirichlet problem says $L:X\to Y$ is a Fredholm operator of index zero. Its null space $N$ is finite dimensional, its range is closed, and the range codimension equals $\dim N$. Hence $\boxed{N=\{0\}⟺L:X\to Y\text{ is invertible}}$; in this case every forcing has a unique solution.

Otherwise there are nonzero homogeneous solutions and finitely many compatibility conditions: $Lu=f$ is solvable exactly when $\ell (f)=0$ for every continuous linear functional on $Y$ vanishing on $L(X)$. When solvable, its solutions are $u_0+N$. This Fredholm solvability condition can be expressed with formal-adjoint null solutions for smoother coefficients. The functional formulation avoids differentiating the merely $C^{0,\alpha }$ coefficients.

The preceding estimate gives constants $A,B>0$ such that, whenever $\Vert w{\Vert }_{H^2}\le R$,Choose $R>0$ so that $B{R}^2\le 1/2$, and then choose $r>0$ so that $Ar\le R/2$. If $\Vert f{\Vert }_2\le r$, the closed ball $B_R(0)$ is mapped into itself.

For $w,z\in B_R(0)$, factor the cubic gradient nonlinearity asThe same Sobolev embedding theorem and the Holder inequality implyThe elliptic estimate therefore yieldsShrinking $R$ further makes $C^{\prime }{R}^2<1$, so $\Phi$ is a contraction mapping of the closed ball. This ball is complete because $H^2(U)\cap H_0^1(U)$ is a Banach space. The contraction mapping theorem gives a fixed point $u=\Phi (u)$, and its defining equation isThus the nonlinear elliptic boundary value problem has a solution for sufficiently small $\Vert f{\Vert }_2$.

Potential-vorticity inversion recovers a balanced velocity field from potential vorticity and specified boundary and global data. In two-dimensional vortex dynamics this is a Poisson equation; shallow-water quasi-geostrophic inversion and stratified quasi-geostrophic inversion use different elliptic boundary value problems. Advance the materially conserved scalar, then invert again to recover the evolving flow. Boundary buoyancy can induce flow with zero interior anomaly. The balance relation and boundary conditions are essential: unrestricted gravity waves are additional degrees of freedom. McIntyre's potential-vorticity account explains the general invertibility principle.

Suppose the inverse $S:C^{0,\alpha }\to C_0^{2,\alpha }$ of a linear elliptic boundary value problem is bounded by $K$, and a nonlinear map satisfiesThen $T(v)=S(Q(v)+f)$ maps a sufficiently small closed ball into itself and is a contraction mapping when $\delta$ and $\Vert f{\Vert }_{C^{0,\alpha }}$ are sufficiently small. The Banach fixed-point theorem gives a solution of $Lu=Q(u)+f$ with zero Dirichlet boundary condition, unique within that ball. Invertibility follows from the Fredholm alternative for an elliptic Dirichlet problem when the homogeneous null space is zero.

The three-dimensional quasi-geostrophic potential vorticity relation defines an anisotropic elliptic boundary value problem. Interior potential vorticity and boundary buoyancy perturbation $b=f_0{\psi }_z$, together with lateral and global conditions, determine the quasi-geostrophic streamfunction up to its allowed constant. If two solutions have the same data, their difference has zero elliptic source and homogeneous boundary data. Integration by parts gives zero integral of $|{\nabla }_h\psi {|}^2+(f_0^2/N^2)|{\psi }_z{|}^2$, proving uniqueness modulo a constant in a connected domain when the boundary terms vanish. Eady edge waves show that zero interior source alone is insufficient.