Potential vorticity ($q$)

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.

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.

In two-dimensional vortex dynamics, solve the Poisson equation for the streamfunction with the prescribed boundary conditions. On the plane its fundamental solution is $(2\pi {)}^{-1}\log |\mathbf{x}|$, so differentiating the displayed inverse gives the planar Biot-Savart kernel. The harmonic function ${\psi }_h$ represents boundary and global-flow data. An additive constant changes no velocity; circulation around holes and a uniform flow are additional degrees of freedom if their boundary or far-field values are not prescribed. The same inverse acts independently on each level in nonrotating layerwise-two-dimensional vortex dynamics.

The shallow-water quasi-geostrophic potential vorticity relation is a modified Helmholtz equation for the quasi-geostrophic streamfunction. Its inverse with specified boundary conditions recovers horizontal geostrophic flow and surface elevation $\eta =f_0\psi /g$. The Rossby deformation radius $L_D=\sqrt{gH}/|f_0|$ is the range of the balanced response: for a horizontal Fourier component, $\hat{\psi }=-\hat{q}/(|{\mathbf{k}}_h{|}^2+L_D^{-2})$, after subtracting background potential vorticity.

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.

The Prandtl ratio compares horizontal and vertical scales in rotating density stratification. In the three-dimensional quasi-geostrophic potential vorticity operator, horizontal curvature scales as $\psi /L^2$ and vertical stretching as $(f_0^2/N^2)\psi /H^2$. Comparable terms therefore require $H/L\sim |f_0|/N$. For constant $N$, the stretched coordinate $Z=(N/|f_0|)z$ makes the operator isotropic. This aspect ratio concerns rotational and stratification effects, not the viscous-to-thermal diffusivity Prandtl number.

For constant nonzero Coriolis parameter and positive constant buoyancy frequency, the anisotropic quasi-geostrophic potential vorticity operator becomes the ordinary three-dimensional Laplacian in the displayed coordinates. Its decaying whole-space inverse is $\psi (\mathbf{X})=-\int [Q({\mathbf{X}}^{\prime })-f]/(4\pi |\mathbf{X}-{\mathbf{X}}^{\prime }|)d^3{X}^{\prime }$. Boundary or far-field data are needed to exclude arbitrary harmonic additions.

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|$.

The spatial gradient of potential vorticity controls the anomaly created when a parcel is displaced. The linear potential-vorticity equation contains ${\mathbf{u}}^{\prime }\cdot \nabla \bar{q}$. A gradient of the Coriolis parameter contributes planetary vorticity, but gradients of relative vorticity and layer stretching can be equally important.

A potential vorticity budget has the formThe material derivative follows the fluid motion and $S_q$ collects forcing or dissipation. With $S_q=0$ it becomes potential-vorticity conservation. For a quasi-geostrophic streamfunction, horizontal advection is $J(\psi ,q)={\psi }_x{q}_y-{\psi }_y{q}_x$.

In unforced adiabatic inviscid motion, potential vorticity is conserved following a parcel:The linearization about a basic state is $({\partial }_t+\bar{\mathbf{u}}\cdot \nabla )q^{\prime }+{\mathbf{u}}^{\prime }\cdot \nabla \bar{q}=0$. This separates advection of a disturbance from its generation by displacement across a basic potential-vorticity gradient.

If quasi-geostrophic potential vorticity initially has gradient $\beta$ across a channel of width $W$ and is homogenized, a height-independent change of zonal velocity obeys $\delta Q=-\delta U_y$. Requiring $\delta U(0)=\delta U(W)=0$ gives uniform final potential vorticity equal to its original channel mean, $\delta U=\beta y(y-W)/2$, and channel-mean momentum change $-\beta W^2/12$. The negative change for $\beta >0$ is wave drag associated with mixing, not a loss of total momentum in an isolated unforced system.

In zonally averaged barotropic quasi-geostrophic approximation flow with zero mean meridional transport, the Taylor identity for quasi-geostrophic flux relates the meridional potential vorticity flux to the convergence of Reynolds stress. A local downgradient closure $\overline{v^{\prime }{q}^{\prime }}=-\kappa {\overline{Q}}_y$, with $\kappa \ge 0$, therefore gives westward mean acceleration where the background gradient is positive. Rossby wave breaking can stir PV into fine filaments; coarse-graining or weak diffusion produces the inferred mixing. Globally the momentum-flux convergence integrates to boundary stresses, so localized westward changes need compensating transport or momentum changes elsewhere in an isolated system.

Ertel potential vorticity is the scalar product of absolute vorticity with the gradient of a materially conserved stratifying scalar. In Boussinesq flow with buoyancy $b$,It is materially conserved by inviscid, adiabatic motion.

For dry atmospheric thermodynamics, choose potential temperature as the materially conserved scalar in Ertel's theorem. The displayed density-normalized scalar has units $\mathrm{K}{\mathrm{m}}^2\mathrm{k}{\mathrm{g}}^{-1}{\mathrm{s}}^{-1}$. Its conservation in adiabatic inviscid motion depends on conservative body forces and a composition-fixed equation of state. It differs in normalization from the Boussinesq approximation convention using buoyancy without the variable density factor. Large-scale balanced weather can be organized by its distribution, with potential-vorticity inversion additionally requiring balance and boundary information.

For a composition-fixed equation of state, let $\dot{\theta }=D\theta /Dt$ and let $\mathbf{F}$ be nonconservative force per unit mass. The vorticity and scalar-gradient evolution equations give the source terms in the title. Spatially varying heating projected along absolute vorticity and the curl of friction can generate or destroy Ertel potential vorticity. Heating that changes potential temperature does not automatically change PV everywhere: its spatial gradient and orientation matter. Unresolved irreversible mixing can also alter the parcelwise conservation picture.

On potential-temperature surfaces, $q=({\zeta }_{\theta }+f)/{\sigma }_{\theta }$, where ${\sigma }_{\theta }=-(1/g)\partial p/\partial \theta$ is isentropic mass density. In adiabatic frictionless flow it is materially conserved. Diagnosing a velocity field also requires balance and boundary data.

For an adiabatic, inviscid parcel, reduced isentropic mass thickness requires reduced absolute vertical vorticity; increased mass thickness requires the reverse. A rigid-upper-boundary geometric cartoon with layer depth $H-h_m$ gives an anticyclonic tendency while crossing a mountain in the Northern Hemisphere. Actual isentropes need not behave like a rigid lid. Blocking, meridional displacement, heating and friction must be assessed before applying that cartoon to real terrain.

${\sigma }_{\theta }=-(1/g)\partial p/\partial \theta$ measures mass per horizontal area per unit potential-temperature interval under hydrostatic balance. It is positive in a statically stable atmosphere.

Heating over land often favours a deeper cyclonic thermal low, while cooling favours a shallower anticyclonic thermal high. Diabatic heating and surface friction violate the assumptions of simple parcel PV conservation; a vertical profile is needed for a quantitative statement.

Spatially varying drag has a vorticity tendency beyond simple damping of existing relative vorticity. The cross-gradient term depends on wind direction and the gradient of surface resistance. Increasing roughness on moving from ocean to land does not alone specify a universal cyclonic or anticyclonic sign. Seasonal heating also changes stratification and isentropic mass thickness, so it cannot generally be replaced by adiabatic PV conservation.

An inviscid fluid with conservative body force, a materially conserved stratifying scalar $\alpha$, and an equation of state making $\nabla \rho ,\nabla p,\nabla \alpha$ coplanar conserves its Ertel potential vorticity. The stretching of absolute vorticity per unit density cancels the deformation of the scalar gradient; the remaining baroclinic triple product vanishes. In an incompressible density-stratified ideal fluid, choosing the advected density as $\alpha$ gives that cancellation directly.

Shallow-water potential vorticity is materially conserved by inviscid, unforced shallow water equations. Linearizing about rest with depth $H_0$ gives an anomaly proportional to

For inviscid nonrotating shallow water equations, columnar horizontal flow has material vorticity equation $D\zeta /Dt=-\zeta {\nabla }_h\cdot {\mathbf{u}}_h$ and depth equation $Dh/Dt=-h{\nabla }_h\cdot {\mathbf{u}}_h$. Thus $\zeta /h$ is conserved, not $\zeta$ itself. Squashing a column to much smaller depth decreases its vertical vorticity. This permits a late draining flow to approximate an irrotational $u_{\theta }=C/r$ profile outside the drain, despite nonzero initial solid-body vorticity. At finite nonzero depth the residual vorticity need not be exactly zero. In a rotating frame, replace $\zeta$ by absolute vorticity $\zeta +f$ in the conserved ratio. Kelvin's circulation theorem remains compatible with these changes because it follows a material loop, not a fixed spatial circle.

For an upstream irrotational current on an f-plane in a rigid layer of height $h$, a circular region of height $h(1-ϵ)$ has the displayed relative vorticity if its fluid retains upstream potential vorticity $f/h$. To leading order in the height contrast, take $u=-{\Psi }_y$, $v={\Psi }_x$ and $\Psi =-Uy+\psi$. Solving the Poisson equation for $\psi$ gives $\psi =-C{r}^2/2$ inside and $\psi =C{a}^2[\log (a/r)-1/2]$ outside, where $C=fϵ/2$. The disturbance is clockwise for positive $fϵ$. Variable-depth transport and the transition across the rim introduce higher-order corrections to this leading columnar inversion.

For the circular topographic anticyclone in a positive uniform current $U$, the opposing disturbance speed is at most $Ca=fϵa/2$. Below the displayed threshold $u$ remains positive everywhere, so no streamline can close. Above it an interior stagnation point at $(0,-U/C)$ is a center, while the exterior point $(0,-C{a}^2/U)$ is a saddle. At equality the two meet at the rim. This is the onset criterion; once closed contours form, the potential vorticity of trapped fluid is not determined solely by its present upstream conditions.

For a steady, slowly varying rotating channel, negligible transverse velocity and zero shallow-water potential vorticity give $u_y=f$. Cross-channel geostrophic balance gives $g{h}_y=-fu$, producing the displayed profiles. The specific energy $E=h+u^2/(2g)$ is constant across the channel, and $E+H$ is constant along it. Zero potential vorticity is not the value in an ordinary finite-depth reservoir at rest in the rotating frame; it requires vanishing absolute vorticity or a deep-reservoir approximation.

At a hydraulic control, the specific-energy function of a zero-potential-vorticity rotating channel flow is stationary with respect to wall depth at fixed discharge and width. Writing $v=u_1+fb/2$ gives $v(v+fb/2)=g{h}_1$ and $Q_c=b{v}^3/g$, hence the displayed discharge. The fully wetted positive-flow solution requires $fb<\sqrt{2g{h}_1}$; at equality the other wall dries and the first wall velocity vanishes. Beyond that limit the assumed two-wet-wall profile must be replaced.

For a constant-depth resting layer on an f-plane, the depth-scaled linear potential vorticity anomaly isIt satisfies $q_t=0$ under the linearized shallow water equations. The physical perturbation of $(f+\zeta )/(H+\eta )$ is $q/H$. The signed scalar $q$, not its absolute value, is used in potential-vorticity inversion and balanced energy pairings.

For small Rossby number and small free-surface displacement, shallow-water potential vorticity reduces, up to a constant factor and additive constant, to relative vorticity plus background PV gradients and the stretching term $-\psi /L_D^2$.

For a uniform upstream current whose free-surface and bottom slopes compensate, steady quasi-geostrophic approximation flow with uniform incoming potential vorticity obeys the displayed inversion equation. Let $\hat{b}=ϵ\cos (x/L)$ on $|x|<a=\pi L/2$ and zero outside. With $C=(fϵ/h_0)/(L^{-2}+L_R^{-2})$, the decaying even response is $\varphi =C[\cos (x/L)+(L_R/L)e^{-a/L_R}\cosh (x/L_R)]$ inside and $\varphi =C(L_R/L)e^{-a/L_R}\cosh (a/L_R)e^{-(|x|-a)/L_R}$ outside. Continuity of $\varphi$ and ${\varphi }^{\prime }$ fixes the coefficients. Its streamlines satisfy $y=y_{\infty }+\varphi (x)/U$, so the response extends beyond the compact ridge through exponentially decaying balanced tails.

On a constant f-plane with finite Rossby deformation radius, a geostrophic uniform current $(U,0)$ has streamfunction $-Uy$ and a sloping free surface. The full shallow-water quasi-geostrophic potential vorticity is $Q=f+{\nabla }^2\Psi -\Psi /L_R^2+fb/H$, so a flat-bottom background has the displayed nonconstant value. Omitting its gradient is not a harmless streamfunction gauge choice. A uniform-current, constant-PV inversion must explicitly supply a compensating background, external balance, or a rigid-lid limit.

In a shallow-water potential-vorticity inversion with the background gradient explicitly compensated, uniform PV gives $({\nabla }^2-L_R^{-2}){\psi }^{\prime }=-fb/H$. The displayed Gaussian function solves it for $b=ϵ[4(1-r^2/a^2)+a^2/L_R^2]e^{-r^2/a^2}$. This is not, without that compensation, a steady solution of the standard finite-depth constant-f current problem: its residual is $U{\psi }_x^{\prime }/L_R^2$. Closed streamlines can occur while Rossby number and surface displacement remain small.

For total streamfunction $\Psi =-Uy+A{e}^{-(x^2+y^2)/a^2}$, the maximum opposing zonal disturbance speed is $\sqrt{2/e}|A|/a$. Below the displayed threshold the zonal velocity never changes sign and streamlines cannot close. Above it two stagnation points appear on the symmetry axis: one has definite Hessian and is a center, the other has indefinite Hessian and is a saddle. Closed contours surround the center. Equality is a degenerate onset, not a finite-area closed-cell regime.

Homogenizing the potential vorticity over a finite constant bottom slope creates the odd anomaly displayed here while leaving the exterior unchanged. Inverting $({\partial }_y^2-L_R^{-2})\psi =\Delta Q$ with decay at both infinities and matching $\psi ,{\psi }_y$ gives a central return current and flanking jets. The deformation radius sets the smoothing length. For $L_R\ll a$, the central surface follows the bottom slope to leading order, but edge corrections and decaying exterior flow are essential.

A positive jump $\Delta Q$ of shallow-water quasi-geostrophic potential vorticity supports a current $U=U_0{e}^{-|y|/L_R}$ with $U_0=L_R\Delta Q/2$. A disturbance has structure $\varphi =A{e}^{-\sqrt{k^2+L_R^{-2}}|y|}$. Integrating the linear potential-vorticity conservation equation across the interface gives $(U_0-c)[{\varphi }_y]+\Delta Q\varphi (0)=0$, hence the displayed dispersion relation. The wave propagates intrinsically against the mean current because the displaced PV interface induces the velocity that advects it; no planetary beta is required.

The barotropic deformation radius is the horizontal scale at which free-surface gravity and rotation balance. It contributes $L_D^{-2}$ to the denominator of the one-layer Rossby-wave dispersion relation.