Hypsometric equation

Atmospheric equation in meteorology

The hypsometric equation, also known as the thickness equation, relates an atmospheric pressure ratio to the equivalent thickness of an atmospheric layer considering the layer mean of virtual temperature, gravity, and occasionally wind. It is derived from the hydrostatic equation and the ideal gas law.

Formulation

The hypsometric equation is expressed as:[1] h = z 2 z 1 = R T v ¯ g ln ( p 1 p 2 ) , {\displaystyle h=z_{2}-z_{1}={\frac {R\cdot {\overline {T_{v}}}}{g}}\,\ln \left({\frac {p_{1}}{p_{2}}}\right),} where:

  • h {\displaystyle h} = thickness of the layer [m],
  • z {\displaystyle z} = geometric height [m],
  • R {\displaystyle R} = specific gas constant for dry air,
  • T v ¯ {\displaystyle {\overline {T_{v}}}} = mean virtual temperature in Kelvin [K],
  • g {\displaystyle g} = gravitational acceleration [m/s2],
  • p {\displaystyle p} = pressure [Pa].

In meteorology, p 1 {\displaystyle p_{1}} and p 2 {\displaystyle p_{2}} are isobaric surfaces. In radiosonde observation, the hypsometric equation can be used to compute the height of a pressure level given the height of a reference pressure level and the mean virtual temperature in between. Then, the newly computed height can be used as a new reference level to compute the height of the next level given the mean virtual temperature in between, and so on.

Derivation

The hydrostatic equation:

p = ρ g z , {\displaystyle p=\rho \cdot g\cdot z,}

where ρ {\displaystyle \rho } is the density [kg/m3], is used to generate the equation for hydrostatic equilibrium, written in differential form:

d p = ρ g d z . {\displaystyle dp=-\rho \cdot g\cdot dz.}

This is combined with the ideal gas law:

p = ρ R T v {\displaystyle p=\rho \cdot R\cdot T_{v}}

to eliminate ρ {\displaystyle \rho } :

d p p = g R T v d z . {\displaystyle {\frac {\mathrm {d} p}{p}}={\frac {-g}{R\cdot T_{v}}}\,\mathrm {d} z.}

This is integrated from z 1 {\displaystyle z_{1}} to z 2 {\displaystyle z_{2}} :

p ( z 1 ) p ( z 2 ) d p p = z 1 z 2 g R T v d z . {\displaystyle \int _{p(z_{1})}^{p(z_{2})}{\frac {\mathrm {d} p}{p}}=\int _{z_{1}}^{z_{2}}{\frac {-g}{R\cdot T_{v}}}\,\mathrm {d} z.}

R and g are constant with z, so they can be brought outside the integral. If temperature varies linearly with z (e.g., given a small change in z), it can also be brought outside the integral when replaced with T v ¯ {\displaystyle {\overline {T_{v}}}} , the average virtual temperature between z 1 {\displaystyle z_{1}} and z 2 {\displaystyle z_{2}} .

p ( z 1 ) p ( z 2 ) d p p = g R T v ¯ z 1 z 2 d z . {\displaystyle \int _{p(z_{1})}^{p(z_{2})}{\frac {\mathrm {d} p}{p}}={\frac {-g}{R\cdot {\overline {T_{v}}}}}\int _{z_{1}}^{z_{2}}\,\mathrm {d} z.}

Integration gives

ln ( p ( z 2 ) p ( z 1 ) ) = g R T v ¯ ( z 2 z 1 ) , {\displaystyle \ln \left({\frac {p(z_{2})}{p(z_{1})}}\right)={\frac {-g}{R\cdot {\overline {T_{v}}}}}(z_{2}-z_{1}),}

simplifying to

ln ( p 1 p 2 ) = g R T v ¯ ( z 2 z 1 ) . {\displaystyle \ln \left({\frac {p_{1}}{p_{2}}}\right)={\frac {g}{R\cdot {\overline {T_{v}}}}}(z_{2}-z_{1}).}

Rearranging:

z 2 z 1 = R T v ¯ g ln ( p 1 p 2 ) , {\displaystyle z_{2}-z_{1}={\frac {R\cdot {\overline {T_{v}}}}{g}}\ln \left({\frac {p_{1}}{p_{2}}}\right),}

or, eliminating the natural log:

p 1 p 2 = e g R T v ¯ ( z 2 z 1 ) . {\displaystyle {\frac {p_{1}}{p_{2}}}=e^{{\frac {g}{R\cdot {\overline {T_{v}}}}}\cdot (z_{2}-z_{1})}.}

Correction

The Eötvös effect can be taken into account as a correction to the hypsometric equation. Physically, using a frame of reference that rotates with Earth, an air mass moving eastward effectively weighs less, which corresponds to an increase in thickness between pressure levels, and vice versa. The corrected hypsometric equation follows:[2] h = z 2 z 1 = R T v ¯ g ( 1 + A ) ln ( p 1 p 2 ) , {\displaystyle h=z_{2}-z_{1}={\frac {R\cdot {\overline {T_{v}}}}{g(1+A)}}\cdot \ln \left({\frac {p_{1}}{p_{2}}}\right),} where the correction due to the Eötvös effect, A, can be expressed as follows: A = 1 g ( 2 Ω u ¯ cos ϕ + u ¯ 2 + v ¯ 2 r ) , {\displaystyle A=-{\frac {1}{g}}\left(2\Omega {\overline {u}}\cos \phi +{\frac {{\overline {u}}^{2}+{\overline {v}}^{2}}{r}}\right),} where

  • Ω {\displaystyle \Omega } = Earth rotation rate,
  • ϕ {\displaystyle \phi } = latitude,
  • r {\displaystyle r} = distance from Earth center to the air mass,
  • u ¯ {\displaystyle {\overline {u}}} = mean velocity in longitudinal direction (east-west), and
  • v ¯ {\displaystyle {\overline {v}}} = mean velocity in latitudinal direction (north-south).

This correction is considerable in tropical large-scale atmospheric motion.

See also

  • Barometric formula
  • Vertical pressure variation

References

  1. ^ "Hypsometric equation - AMS Glossary". American Meteorological Society. Retrieved 12 March 2013.
  2. ^ Ong, H.; Roundy, P.E. (2019). "Nontraditional hypsometric equation". Q. J. R. Meteorol. Soc. 146 (727): 700–706. arXiv:2011.09576. doi:10.1002/qj.3703.