Potential vorticity

Vilhelm Bjerknes generalized Helmholtz's vorticity equation (1858) and Kelvin's circulation theorem (1869) to inviscid, geostrophic, and baroclinic fluids,[1] i.e., fluids of varying density in a rotational frame which has a constant angular speed. If we define circulation as the integral of the tangent component of velocity around a closed fluid loop and take the integral of a closed chain of fluid parcels, we obtain

${\displaystyle {\frac {DC}{Dt}}=-\oint {\frac {1}{\rho }}\nabla p\cdot \mathrm {d} \mathbf {r} -2\Omega {\frac {DA_{e}}{Dt}},}$ (1)

where ${\textstyle {\frac {D}{Dt}}}$ is the time derivative in the rotational frame (not inertial frame), ${\displaystyle C}$ is the relative circulation, ${\displaystyle A_{e}}$ is projection of the area surrounded by the fluid loop on the equatorial plane, ${\displaystyle \rho }$ is density, ${\displaystyle p}$ is pressure, and ${\displaystyle \Omega }$ is the frame's angular speed. With Stokes' theorem, the first term on the right-hand-side can be rewritten as

${\displaystyle {\frac {DC}{Dt}}=\int _{A}{\frac {\nabla \rho \times \nabla p}{\rho ^{2}}}\cdot \mathrm {d} \mathbf {A} -2\Omega {\frac {DA_{e}}{Dt}},}$(2)

which states that the rate of the change of the circulation is governed by the variation of density in pressure coordinates and the equatorial projection of its area, corresponding to the first and second terms on the right hand side. The first term is also called the "solenoid term". Under the condition of a barotropic fluid with a constant projection area ${\displaystyle A_{e}}$, the Bjerknes circulation theorem reduces to Kelvin's theorem. However, in the context of atmospheric dynamics, such conditions are not a good approximation: if the fluid circuit moves from the equatorial region to the extratropics, ${\displaystyle A_{e}}$ is not conserved. Furthermore, the complex geometry of the material circuit approach is not ideal for making an argument about fluid motions.

Carl Rossby proposed in 1939[4] that, instead of the full three-dimensional vorticity vector, the local vertical component of the absolute vorticity ${\displaystyle \zeta _{a}}$ is the most important component for large-scale atmospheric flow. Also, the large-scale structure of a two-dimensional non-divergent barotropic flow can be modeled by assuming that ${\displaystyle \zeta _{a}}$ is conserved. His later paper in 1940[5] relaxed this theory from 2D flow to quasi-2D shallow water equations on a beta plane. In this system, the atmosphere is separated into several incompressible layers stacked upon each other, and the vertical velocity can be deduced from integrating the convergence of horizontal flow. For a one-layer shallow water system without external forces or diabatic heating, Rossby showed that

${\displaystyle {\frac {D}{Dt}}\left({\frac {\zeta +f}{h}}\right)=0}$, (3)

where ${\displaystyle \zeta }$ is the relative vorticity, ${\displaystyle h}$ is the layer depth, and ${\displaystyle f}$ is the Coriolis parameter. The conserved quantity in equation (3) is later called the shallow water potential vorticity. For an atmosphere with multiple layers, with each layer having constant potential temperature, the above equation takes the form

${\displaystyle {\frac {D}{Dt}}\left({\frac {\zeta _{\theta }+f}{\Delta }}\right)=0,}$ (4)

in which ${\displaystyle \zeta _{\theta }}$ is the relative vorticity on an isentropic surface—a surface of constant potential temperature, and ${\displaystyle \Delta =-\delta p/g}$ is a measure of the weight of unit cross-section of an individual air column inside the layer.

Interpretation

Convergence and divergence of an air parcel

Equation (3) is the atmospheric equivalent to angular momentum. For example, a spinning ice skater with her arms spread out laterally can accelerate her rate of spin by contracting her arms. Similarly, when a vortex of air is broadened, it in turn spins more slowly. When the air converges horizontally, the air speed increases to maintain potential vorticity, and the vertical extent increases to conserve mass. On the other hand, divergence causes the vortex to spread, slowing down the rate of spin.

Hans Ertel generalized Rossby's work via an independent paper published in 1942.[6][7] By identifying a conserved quantity following the motion of an air parcel, it can be proved that a certain quantity called the Ertel potential vorticity is also conserved for an idealized continuous fluid. We look at the momentum equation and the mass continuity equation of an idealized compressible fluid in Cartesian coordinates:

${\displaystyle {\frac {\partial \mathbf {v} }{\partial t}}+{\frac {1}{2}}\nabla \mathbf {v^{2}} -\mathbf {v} \times (\nabla \times \mathbf {v} )+2\mathbf {\Omega } \times \mathbf {v} =-\nabla \Phi -{\frac {1}{\rho }}\nabla p,}$ (5)

${\displaystyle {\frac {D\rho }{Dt}}=-\rho \nabla \cdot \mathbf {v} ,}$ (6)

where ${\displaystyle \Phi }$ is the geopotential height. Writing the absolute vorticity as ${\textstyle \mathbf {\zeta _{a}} =\nabla \times \mathbf {v} +2\mathbf {\Omega } }$, ${\displaystyle 1/\rho }$ as ${\displaystyle \alpha }$, and then take the curl of the full momentum equation (5), we have

${\displaystyle {\frac {\partial }{\partial t}}\nabla \times \mathbf {v} -\nabla \times (\mathbf {v} \times \mathbf {\zeta _{a}} )=\nabla p\times \nabla \alpha .}$ (7)

Consider ${\displaystyle \psi =\psi (\mathbf {r} ,t)}$ to be a hydrodynamical invariant, that is, ${\textstyle {\frac {D\psi }{Dt}}}$ equals to zero following the fluid motion in question. Scalar multiplication of equation (7) by ${\displaystyle \nabla \psi }$, and note that ${\textstyle {\frac {\partial }{\partial t}}\mathbf {\zeta _{a}} ={\frac {\partial }{\partial t}}\nabla \times \mathbf {v} }$, we have

${\displaystyle \nabla \psi \cdot {\frac {\partial }{\partial t}}\mathbf {\zeta _{a}} -\nabla \psi \cdot [\nabla \times (\mathbf {v} \times \mathbf {\zeta _{a}} )]=\nabla \psi \cdot (\nabla p\times \nabla \alpha ).}$ (8)

The second term on the left-hand side of equation (8) is equal to ${\textstyle \nabla \cdot [\nabla \psi \times (\mathbf {v} \times \mathbf {\zeta _{a}} )]-(\mathbf {v} \times \zeta _{a})\cdot (\nabla \times \nabla \psi )}$, in which the second term is zero. From the triple vector product formula, we have

${\displaystyle {\begin{aligned}\nabla \psi \times (\mathbf {v} \times \mathbf {\zeta _{a}} )&=\mathbf {v} (\mathbf {\zeta _{a}} \cdot \nabla \psi )-\mathbf {\zeta _{a}} (\mathbf {v} \cdot \nabla \psi )\\&=\mathbf {v} (\mathbf {\zeta _{a}} \cdot \nabla \psi )+\mathbf {\zeta _{a}} {\frac {\partial \psi }{\partial t}},\end{aligned}}}$ (9)

where the second row is due to the fact that ${\displaystyle \psi }$ is conserved following the motion, ${\textstyle {\frac {D\psi }{Dt}}=0}$. Substituting equation (9) into equation (8) above,

${\displaystyle \nabla \psi \cdot {\frac {\partial }{\partial t}}\mathbf {\zeta _{a}} +\mathbf {v} \cdot \nabla (\mathbf {\zeta _{a}} \cdot \nabla \psi )+(\mathbf {\zeta _{a}} \cdot \nabla \psi )\nabla \cdot \mathbf {v} +\mathbf {\zeta _{a}} \cdot {\frac {\partial }{\partial t}}\nabla \psi =\nabla \psi \cdot (\nabla p\times \nabla \alpha ).}$ (10)

Combining the first, second, and fourth term in equation (10) can yield ${\textstyle {\frac {D}{Dt}}(\mathbf {\zeta _{a}} \cdot \nabla \psi )}$. Dividing by ${\textstyle \rho }$ and using a variant form of mass continuity equation,${\textstyle {\frac {1}{\rho }}\nabla \cdot \mathbf {v} =-{\frac {1}{\rho ^{2}}}{\frac {D\rho }{Dt}}={\frac {D\alpha }{Dt}}}$, equation (10) gives

${\displaystyle \alpha {\frac {D}{Dt}}(\mathbf {\zeta _{a}} \cdot \nabla \psi )+(\mathbf {\zeta _{a}} \cdot \nabla \psi ){\frac {D\alpha }{Dt}}=\alpha \nabla \psi \cdot (\nabla p\times \nabla \alpha )}$ (11)

If the invariant ${\textstyle \psi }$ is only a function of pressure ${\textstyle p}$ and density ${\textstyle \rho }$, then its gradient is perpendicular to the cross product of ${\textstyle \nabla p}$ and ${\textstyle \nabla \rho }$, which means that the right-hand side of equation (11) is equal to zero. Specifically for the atmosphere, potential temperature is chosen as the invariant for frictionless and adiabatic motions. Therefore, the conservation law of Ertel's potential vorticity is given by

${\displaystyle {\frac {D}{Dt}}PV=0.}$ (12)

the potential vorticity is defined as

${\displaystyle PV={\frac {\mathbf {\zeta _{a}} \cdot \nabla \theta }{\rho }},}$ (13)

where ${\displaystyle \rho }$ is the fluid density, ${\displaystyle \mathbf {\zeta _{a}} }$ is the absolute vorticity and ${\displaystyle \nabla \theta }$ is the gradient of potential temperature. It can be shown through a combination of the first law of thermodynamics and momentum conservation that the potential vorticity can only be changed by diabatic heating (such as latent heat released from condensation) or frictional processes.

If the atmosphere is stably stratified so that the potential temperature ${\textstyle \theta }$ increases monotonically with height, ${\textstyle \theta }$ can be used as a vertical coordinate instead of ${\textstyle z}$. In the ${\textstyle (x,y,\theta )}$ coordinate system, "density" is defined as ${\textstyle \sigma \equiv -g^{-1}\partial p/\partial \theta }$. Then, if we start the derivation from the horizontal momentum equation in isentropic coordinates, Ertel PV takes a much simpler form[8]

${\displaystyle PV_{\theta }={\frac {f+\mathbf {k} \cdot (\nabla _{\theta }\times \mathbf {v} )}{\sigma }}}$ (14)

where ${\textstyle \mathbf {k} }$ is the local vertical vector of unit length and ${\textstyle \nabla _{\theta }}$ is the 3-dimensional gradient operator in isentropic coordinates. It can be seen that this form of potential vorticity is just the continuous form of Rossby's isentropic multi-layer PV in equation (4).

One of the simplest but nevertheless insightful balancing conditions is in the form of quasi-geostrophic equations. This approximation basically states that for three-dimensional atmospheric motions that are nearly hydrostatic and geostrophic, their geostrophic part can be determined approximately by the pressure field, whereas the ageostrophic part governs the evolution of the geostrophic flow. The potential vorticity in the quasi-geostrophic limit (QGPV) was first formulated by Charney and Stern in 1960.[10] Similar to Chapter 6.3 in Holton 2004,[8] we start from horizontal momentum (15), mass continuity (16), hydrostatic (17), and thermodynamic (18) equations on a beta plane, while assuming that the flow is inviscid and hydrostatic,

${\displaystyle {\frac {D_{g}\mathbf {u_{g}} }{Dt}}+f_{0}\mathbf {k} \times \mathbf {u_{a}} +\beta y\mathbf {k} \times \mathbf {u_{g}} =0}$ (15)

${\displaystyle \nabla _{hp}\cdot \mathbf {u_{a}} +{\frac {\partial \omega }{\partial p}}=0}$ (16)

${\displaystyle {\frac {\partial \Phi }{\partial p}}=-{\frac {R\pi }{p}}\theta }$ (17)

${\displaystyle {\frac {D_{g}\theta }{Dt}}+\omega {\frac {d\theta _{0}}{dp}}={\frac {J}{c_{p}\pi }}}$ (18)

where ${\textstyle {\frac {D_{g}}{Dt}}={\frac {\partial }{\partial t}}+u_{g}{\frac {\partial }{\partial x}}+v_{g}{\frac {\partial }{\partial y}}}$ represents the geostrophic evolution, ${\textstyle \pi =(p/ps)^{R/c_{p}}}$, ${\textstyle J}$ is the diabatic heating term in ${\textstyle Js^{-1}kg^{-1}}$, ${\textstyle \Phi }$ is the geopotential height, ${\textstyle \mathbf {u_{g}} }$ is the geostrophic component of horizontal velocity, ${\textstyle \mathbf {u_{a}} }$ is the ageostrophic velocity, ${\textstyle \nabla _{hp}}$ is horizontal gradient operator in (x, y, p) coordinates. With some manipulation (see Quasi-geostrophic equations or Holton 2004, Chapter 6 for details), one can arrive at a conservation law

${\displaystyle {\frac {D_{g}}{Dt}}q=-f_{0}{\frac {\partial }{\partial p}}\left(\sigma {\frac {\kappa J}{p}}\right),}$ (19)

where ${\textstyle \sigma =-{\frac {R\pi }{p}}{\frac {d\theta _{0}}{dp}}}$ is the spatially averaged dry static stability. Assuming that the flow is adiabatic, which means ${\textstyle J=0}$, we have the conservation of QGPV. The conserved quantity ${\textstyle q}$ takes the form

${\displaystyle {q={{{1 \over f_{o}}{\nabla ^{2}\Phi }}+{f}+{{\partial \over \partial p}\left({{f_{o} \over \sigma }{\partial \Phi \over \partial p}}\right)}}},}$ (20)

which is the QGPV, and it is also known as the pseudo-potential-vorticity. Apart from the diabatic heating term on the right-hand-side of equation(19), it can also be shown that QGPV can be changed by frictional forces.

The Ertel PV reduces to the QGPV if one expand the Ertel PV to the leading order, and assume that the evolution equation is quasi-geostrophic, i.e., ${\textstyle {\frac {D}{Dt}}\approx {\frac {D_{g}}{Dt}}}$.[3] Because of this factor, one should also note that the Ertel PV conserves following air parcel on an isentropic surface and is therefore a good Lagrangian tracer, whereas the QGPV is conserved following large-scale geostrophic flow. QGPV has been widely used in depicting large-scale atmospheric flow structures, as discussed in the section #PV invertibility principle;

Apart from being a Lagrangian tracer, the potential vorticity also gives dynamical implications via the invertibility principle. For a 2-dimensional ideal fluid, the vorticity distribution controls the stream function by a Laplace operator,

${\displaystyle {\zeta ={\nabla ^{2}\Psi }},}$ (21)

where ${\displaystyle \zeta }$ is the relative vorticity, and ${\displaystyle \Psi }$ is the streamfunction. Hence from the knowledge of vorticity field, the operator can be inverted and the stream function can be calculated. In this particular case (equation 21), vorticity gives all the information needed to deduce motions, or streamfunction, thus one can think in terms of vorticity to understand the dynamics of the fluid. A similar principle was originally introduced for the potential vorticity in three-dimensional fluid in the 1940s by Kleinschmit, and was developed by Charney and Stern in their quasi-geostrophic theory.[11]

Despite theoretical elegance of Ertel's potential vorticity, early applications of Ertel PV are limited to tracer studies using special isentropic maps. It is generally insufficient to deduce other variables from the knowledge of Ertel PV fields only, since it is a product of wind (${\textstyle \zeta _{a}}$) and temperature fields (${\textstyle \theta }$ and ${\textstyle \rho }$). However, large-scale atmospheric motions are inherently quasi-static; wind and mass fields are adjusted and balanced against each other (e.g., gradient balance, geostrophic balance). Therefore, other assumptions can be made to form a closure and deduce the complete structure of the flow in question:[2]

(1) introduce balancing conditions of certain form. These conditions must be physically realizable and stable without instabilities such as static instability. Also, the space and time scales of the motion must be compatible with the assumed balance;

(2) specify a certain reference state, such as distribution of temperature, potential temperature, or geopotential height;

(3) assert proper boundary conditions and invert the PV field globally.

The first and second assumptions are expressed explicitly in the derivation of quasi-geostrophic PV. Leading-order geostrophic balance is used as the balancing condition. The second-order terms such as ageostrophic winds, perturbations of potential temperature and perturbations of geostrophic height should have consistent magnitude, i.e., of the order of Rossby number. The reference state is zonally averaged potential temperature and geopotential height. The third assumption is apparent even for 2-dimensional vorticity inversion because inverting the Laplace operator in equation (21), which is a second-order elliptic operator, requires knowledge of the boundary conditions.

For example, in equation (20), invertibility implies that given the knowledge of ${\textstyle q}$, the Laplace-like operator can be inverted to yield geopotential height ${\textstyle \Phi }$. ${\textstyle \Phi }$ is also proportional to the QG streamfunction ${\textstyle \Psi }$ under the quasi-geostrophic assumption. The geostrophic wind field can then be readily deduced from ${\textstyle \Psi }$. Lastly, The temperature field ${\textstyle \theta }$ is given by substituting ${\textstyle \Phi }$ into the hydrostatic equation (17).

• Vorticity
• Circulation (fluid dynamics)

• AMS Glossary entry
• Encyclopedia article on Potential Vorticity by Michael E. McIntyre