Codex @codex  0

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.

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.

Inertial instability occurs when displaced fluid parcels amplify their displacement by conserving absolute momentum in a rotating flow. For a uniformly stratified basic state on an f-plane, a necessary condition is that the Ertel potential vorticity have the opposite sign to the Coriolis parameter.

In the limiting model of infinitely strong stable density stratification, small Froude number and fixed vertical scales suppress vertical displacement. Each horizontal level evolves by two-dimensional vortex dynamics; $z$ labels independent copies of the horizontal Poisson equation. The leading normalized Ertel potential vorticity is $Q/N^2=\zeta$. This is a constrained asymptotic model: finite buoyancy frequency restores vertical motion and coupling, and arbitrarily short vertical scales invalidate the independent-layer limit. McIntyre's discussion of potential-vorticity inversion explains the limiting independent-layer picture.

For a hydrostatic, stably stratified atmosphere, define isentropic mass density ${\sigma }_{\theta }=-(1/g)\partial p/\partial \theta >0$ and relative vertical vorticity evaluated on an isentrope, ${\zeta }_{\theta }$. Then isentropic potential vorticity isIt is the hydrostatic isentropic form of dry Ertel potential vorticity, $q_E={\mathbf{\omega }}_a\cdot \nabla \theta /\rho$. For inviscid, adiabatic motion with fixed composition and conservative body forces, $D{q}_E/Dt=0$, so parcels carry their PV along isentropes. Stable stratification makes the coordinate transformation well defined.

Heating and friction spoil simple material conservation. More explicitly, with $\dot{\theta }=D\theta /Dt$ and nonconservative force per mass $\mathbf{F}$,Thus the distribution is especially useful away from strongly diabatic clouds and frictional boundary layers. Explaining a velocity field by potential-vorticity inversion additionally requires an appropriate balanced-flow approximation and boundary information; a PV value alone does not specify the entire flow. ^[ref]

The isentropic potential vorticity is Ertel potential vorticity expressed on constant potential temperature surfaces:Here ${\zeta }_{\theta }$ is the vertical curl of horizontal velocity differentiated on an isentropic surface, and ${\sigma }_{\theta }$ is mass per unit area per unit potential-temperature interval. For smooth inviscid adiabatic motion, $Dq/Dt=0$. Conserved anomalies are advected and deformed, explaining much balanced large-scale weather variability. Diagnosing velocity from PV additionally needs balance, stratification and boundary information. Strong latent/radiative heating, friction and irreversible mixing generate or modify PV, limiting a purely conservative explanation.^[ref]

Choose $y$ northward and $z$ upward in the Northern Hemisphere. In a locally linear Boussinesq approximation, write the background buoyancy as $\bar{b}=b_{\ast }+By+N^{2}z$, with horizontal buoyancy gradient $B<0$ and positive squared buoyancy frequency $N^2$. The stable density stratification resists a purely vertical displacement, but the isopycnals slope upward poleward:The geopotential surfaces are horizontal in this local approximation. Their inclination differs from that of the isopycnals, which is the relevant baroclinicity.

Exchange two equal small fluid elements, each of volume $V$, between positions separated by $(\delta y,\delta z)$. The adiabatic parcel exchange has no thermal diffusion, so each parcel retains its buoyancy. The background buoyancy difference is $\Delta \bar{b}=B\delta y+N^2\delta z$. Because the variable-density contribution to gravitational potential energy is $-{\rho }_{0}Vbz$, the change in the pair's gravitational potential energy isFor a poleward-rising exchange, set $\delta z=s\delta y$ with $s>0$. ThenThe trajectory lies between a geopotential surface and an isopycnal. The poleward-rising parcel starts more buoyant than the parcel it replaces, and the equatorward-sinking parcel starts less buoyant: the buoyancy does positive work on both. This releases available potential energy into kinetic energy, allowing eddies to amplify if the full rotating dynamics supports the displacement and dissipation does not remove the gain.

A trajectory steeper than the isopycnals, $s>s_{\rho }$, instead costs gravitational potential energy; an equatorward-rising trajectory also costs it. Motion exactly along a geopotential surface or an isopycnal gives zero change in this ideal exchange. At fixed $\delta y$, the largest release occurs at $s=s_{\rho }/2$, since $N^{2}s(s-s_{\rho })$ is minimized there.

This is an energetic possibility argument for baroclinic instability, not a sufficient normal mode instability criterion. In particular, unrestricted parcel motion must not be assumed in a rotating fluid: conservation of angular momentum can impose an additional restoring effect. Symmetric instability has a separate Ertel potential vorticity criterion; ordinary Eady instability can occur without violating that criterion.

The quasi-geostrophic potential-vorticity equation is the leading balanced form of potential-vorticity conservation for an inviscid rotating fluid with stable density stratification, small Rossby number and an adiabatic process for each fluid element. Its parent conservation law is Ertel potential vorticity conservation. Horizontal geostrophic balance, vertical hydrostatic balance, and the beta-plane approximation reduce it to advection of one scalar by the leading horizontal geostrophic flow.

One can see its origin by combining the vertical vorticity equation with the thermal equation. Set $D_g={\partial }_t+{\mathbf{u}}_g\cdot {\nabla }_h$, ${\mathbf{u}}_g=(-{\psi }_y,{\psi }_x)$, and ${\zeta }_g={\nabla }_h^2\psi$. The quasi-geostrophic approximation givesHere $b$ is the buoyancy perturbation relative to the reference stable density stratification, $N=N(z)$ its buoyancy frequency, and $w$ the small vertical velocity. The stretching of absolute vorticity in the first equation is offset by the change of stratification represented by the second. More precisely,The last dot product is zero because thermal-wind balance makes ${\partial }_z{\mathbf{u}}_g$ perpendicular to ${\nabla }_{h}b$. Adding therefore yieldsThe constant $f_0$ has no effect on the tendency, but its inclusion makes the connection with absolute vorticity clear. Vertical advection of $q$ is beyond the retained order; the vertical velocity has already supplied the stretching and thermodynamic terms.

For a simple balanced weather forecast, the essential cycle is diagnose, advect, invert:

This is an elliptic boundary value problem, not a local algebraic diagnosis. Its solution supplies the next geostrophic flow and, through $b=f_0{\psi }_z$, the next thermal field. Repeat the cycle to move the balanced weather systems forward.

Interior potential vorticity alone is insufficient: boundary thermal anomalies can induce flow even with zero interior anomaly, as in the Eady model. The method filters fast inertia-gravity waves and is useful for large-scale extratropical weather, while tropical dynamics, convection, moisture, friction and diabatic heating require additional physics. If vertical motion is wanted, the quasi-geostrophic omega equation supplies a further diagnostic calculation.

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.

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

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.

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.