Omega equation

Derivation

The derivation of the ${\displaystyle \omega }$ equation is based on the vertical component of the vorticity equation, and the thermodynamic equation. The vertical vorticity equation for a frictionless atmosphere may be written using pressure as the vertical coordinate:

Here ${\displaystyle \xi }$ is the relative vorticity, ${\displaystyle V}$ the horizontal wind velocity vector, whose components in the ${\displaystyle x}$ and ${\displaystyle y}$ directions are ${\displaystyle u}$ and ${\displaystyle v}$ respectively, ${\displaystyle \eta }$ the absolute vorticity ${\displaystyle \xi +f}$, ${\displaystyle f}$ is the Coriolis parameter, ${\displaystyle \omega =\frac{dp}{dt}}$ the material derivative of pressure ${\displaystyle p}$, ${\displaystyle \hat{k}}$ is the unit vertical vector, ${\displaystyle \mathrm{\nabla }}$ is the isobaric Del (grad) operator, ${\displaystyle \left(\xi \frac{\mathrm{\partial }\omega }{\mathrm{\partial }p}-\omega \frac{\mathrm{\partial }\xi }{\mathrm{\partial }p}\right)}$ is the vertical advection of vorticity and ${\displaystyle \hat{k}\cdot \mathrm{\nabla }\omega \times \frac{\mathrm{\partial }V}{\mathrm{\partial }p}}$ represents the "tilting" term or transformation of horizontal vorticity into vertical vorticity.

The thermodynamic equation may be written as:

where ${\displaystyle k\equiv \left(\frac{\mathrm{\partial }Z}{\mathrm{\partial }p}\right)\frac{\mathrm{\partial }}{\mathrm{\partial }p}\ln \theta }$, in which ${\displaystyle Q}$ is the heating rate (supply of energy per unit time and unit mass), ${\displaystyle C_{\text{p}}}$ is the specific heat of dry air, ${\displaystyle R}$ is the gas constant for dry air, ${\displaystyle \theta }$ is the potential temperature and ${\displaystyle \varphi }$ is geopotential ${\displaystyle (gZ)}$.

The ${\displaystyle \omega }$ equation (1) is obtained from equation (2) and (3) by casting both equations in terms of geopotential Z, and eliminating time derivatives based on the physical assumption that thermal wind imbalance remains small across time, or d/dt(imbalance) = 0. For the first step, the relative vorticity must be approximated as the geostrophic vorticity: $${\displaystyle \xi =\frac{g}{f}{\mathrm{\nabla }}^{2}Z}$$

Expanding the final "tilting" term in (2) into Cartesian coordinates (although we will soon neglect it), the vorticity equation reads:

Differentiating (4) with respect to ${\displaystyle p}$ gives:

Taking the Laplacian (${\displaystyle {\mathrm{\nabla }}^2}$) of (3) gives:

Adding (5) to g/f times (6), substituting ${\displaystyle gk=\sigma }$, and approximating horizontal advection with geostrophic advection (using the Jacobian formalism) gives:

Equation (7) is now a diagnostic, linear differential equation for ${\displaystyle \omega }$, which can be split into two terms, namely ${\displaystyle {\omega }_1}$ and ${\displaystyle {\omega }_2}$, such that:

and

where ${\displaystyle {\omega }_1}$ is the vertical velocity attributable to all the flow-dependent advective tendencies in Equation (8), and ${\displaystyle {\omega }_2}$ is the vertical velocity due to the non-adiabatic heating, which includes the latent heat of condensation, sensible heat fluxes, radiative heating, etc. (Singh & Rathor, 1974). Since all advecting velocities in the horizontal have been replaced with geostrophic values, and geostrophic winds are nearly nondivergent, neglect of vertical advection terms is a consistent further assumption of the quasi-geostrophic set, leaving only the square bracketed term in Eqs. (7-8) to enter (1).

Interpretation

Equation (1) for adiabatic ${\displaystyle {\omega }_{QG}}$ is used by meteorologists and operational weather forecasters to anticipate where upward motion will occur on synoptic charts. For sinusoidal or wavelike motions, where Laplacian operators act simply as a negative sign, and the equation's meaning can be expressed with words indicating the sign of the effect: Upward motion is driven by positive vorticity advection increasing with height (or PVA for short), plus warm air advection (or WAA for short). The opposite signed case is logically opposite, for this linear equation.

In a location where the imbalancing effects of adiabatic advection are acting to drive upward motion (where ${\displaystyle {\omega }_{QG}<0}$ in Eq. 1), the inertia of the geostrophic wind field (that is, its propensity to carry on forward) is creating a demand for decreasing thickness ${\displaystyle \frac{-\mathrm{\partial }Z}{\mathrm{\partial }p}}$ in order for thermal wind balance to continue to hold. For instance, when there is an approaching upper-level cyclone or trough above the level in question, the part of ${\displaystyle {\omega }_{QG}}$ attributable to the first term in Eq. 1 is upward motion needed to create the increasingly cool air column that is required hypsometrically under the falling heights. That adiabatic reasoning must be supplemented by an appreciation of feedbacks from flow-dependent heating, such as latent heat release. If latent heat is released as air cools, then an additional upward motion will be required based on Eq. (9) to counteract its effect, in order to still create the necessary cool core. Another way to think about such a feedback is to consider an effective static stability that is smaller in saturated air than in unsaturated air, although a complication of that view is that latent heating mediated by convection need not be vertically local to the altitude where cooling by ${\displaystyle {\omega }_{QG}}$ triggers its formation. For this reason, retaining a separate Q term like Equation (9) is a useful approach.