Omega equation explained

The omega equation is a culminating result in synoptic-scale meteorology. It is an elliptic partial differential equation, named because its left-hand side produces an estimate of vertical velocity, customarily expressed by symbol

\omega

, in a pressure coordinate measuring height the atmosphere. Mathematically,

\omega=

dp
dt
, where

{d\overdt}

represents a material derivative. The underlying concept is more general, however, and can also be applied[1] to the Boussinesq fluid equation system where vertical velocity is

w=

dz
dt
in altitude coordinate z.

Concept and summary

Vertical wind is crucial to weather and storms of all types. Even slow, broad updrafts can create convective instability or bring air to its lifted condensation level creating stratiform cloud decks. Unfortunately, predicting vertical motion directly is difficult. For synoptic scales in Earth's broad and shallow troposphere, the vertical component of Newton's law of motion is sacrificed in meteorology's primitive equations, by accepting the hydrostatic approximation. Instead, vertical velocity must be solved through its link to horizontal laws of motion, via the mass continuity equation. But this presents further difficulties, because horizontal winds are mostly geostrophic, to a good approximation. Geostrophic winds merely circulate horizontally, and do not significantly converge or diverge in the horizontal to provide the needed link to mass continuity and thus vertical motion.

The key insight embodied by the quasi-geostrophic omega equation is that thermal wind balance (the combination of hydrostatic and geostrophic force balances above) holds throughout time, even though the horizontal transport of momentum and heat by geostrophic winds will often tend to destroy that balance. Logically, then, a small non-geostrophic component of the wind (one which is divergent, and thus connected to vertical motion) must be acting as a secondary circulation to maintain balance of the geostrophic primary circulation. The quasi-geostrophic omega

\omegaQG

is the hypothetical vertical motion whose adiabatic cooling or warming effect (based on the atmosphere's static stability) would prevent thermal wind imbalance from growing with time, by countering the balance-destroying (or imbalance-creating) effects of advection. Strictly speaking, QG theory approximates both the advected momentum and the advecting velocity as given by the geostrophic wind.

In summary, one may consider the vertical velocity that results from solving the omega equation as that which would be needed to maintain geostrophy and hydrostasy in the face of advection by the geostrophic wind.[2]

The equation reads:

where

f

is the Coriolis parameter,

\sigma

is related to the static stability,

Vg

is the geostrophic velocity vector,

\zetag

is the geostrophic relative vorticity,

\phi

is the geopotential,
2
\nabla
H
is the horizontal Laplacian operator and

\nablaH

is the horizontal del operator.[3] Its sign and sense in typical weather applications is: upward motion is produced by positive vorticity advection above the level in question (the first term), plus warm advection (the second term).

Derivation

The derivation of the

\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

\xi

is the relative vorticity,

V

the horizontal wind velocity vector, whose components in the

x

and

y

directions are

u

and

v

respectively,

η

the absolute vorticity

\xi+f

,

f

is the Coriolis parameter,

\omega=

dp
dt
the material derivative of pressure

p

,

\hatk

is the unit vertical vector,

\nabla

is the isobaric Del (grad) operator,

\left(\xi

\partial\omega
\partialp

-\omega

\partial\xi
\partialp

\right)

is the vertical advection of vorticity and

\hatk\nabla\omega x

\partialV
\partialp

represents the "tilting" term or transformation of horizontal vorticity into vertical vorticity.[4]

The thermodynamic equation may be written as:where

k\equiv\left(

\partialZ
\partialp

\right)

\partial
\partialp

ln\theta

, in which

Q

is the heating rate (supply of energy per unit time and unit mass),

Cp

is the specific heat of dry air,

R

is the gas constant for dry air,

\theta

is the potential temperature and

\phi

is geopotential

(gZ)

.

The

\omega

equation is obtained from equation and 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:\xi = \frac\nabla^2 Z

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

Differentiating with respect to

p

gives:

Taking the Laplacian (

\nabla2

) of gives:

Adding to g/f times, substituting

gk=\sigma

, and approximating horizontal advection with geostrophic advection (using the Jacobian formalism) gives:

Equation is now a diagnostic, linear differential equation for

\omega

, which can be split into two terms, namely

\omega1

and

\omega2

, such that:andwhere

\omega1

is the vertical velocity attributable to all the flow-dependent advective tendencies in Equation, and

\omega2

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

\omegaQG

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,[5] 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

\omegaQG<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
-\partialZ
\partialp
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

\omegaQG

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

\omegaQG

triggers its formation. For this reason, retaining a separate Q term like Equation (9) is a useful approach.[6]

References

  1. Davies. Huw. 2015. The Quasigeostrophic Omega Equation: Reappraisal, Refinements, and Relevance.. Monthly Weather Review. 143. 1. 3–25. 10.1175/MWR-D-14-00098.1. 2015MWRv..143....3D. free.
  2. Book: Holton, James. An Introduction to Dynamic Meteorology. Elsevier Academic Press. 2004. 0123540151.
  3. Holton, J.R., 1992, An Introduction to Dynamic Meteorology Academic Press, 166-175
  4. Singh & Rathor, 1974, Reduction of the Complete Omega Equation to the Simplest Form, Pure and Applied Geophysics, 112, 219-223
  5. Web site: Quasi-Geostrophic Omega Equation Lab. METEd, CoMET program. 10 November 2019.
  6. Nie. Ji. Fan. Bowen. 2019-06-19. Roles of Dynamic Forcings and Diabatic Heating in Summer Extreme Precipitation in East China and the Southeastern United States. Journal of Climate. 32. 18. 5815–5831. 10.1175/JCLI-D-19-0188.1. 0894-8755. 2019JCli...32.5815N. free.

External links