Scale height



In atmospheric, earth, and planetary sciences, a scale height, usually denoted by the capital letter H, is a distance (vertical or radial) over which a physical quantity decreases by a factor of e (the base of natural logarithms, approximately 2.718).

Scale height used in a simple atmospheric pressure model
For planetary atmospheres, scale height is the increase in altitude for which the atmospheric pressure decreases by a factor of e. The scale height remains constant for a particular temperature. It can be calculated by

$$H = \frac{k_\text{B}T}{mg}$$ or equivalently $$H = \frac{RT}{Mg}$$ where:
 * kB = Boltzmann constant =
 * R = gas constant
 * T = mean atmospheric temperature in kelvins = 250 K for Earth
 * m = mean mass of a molecule (units kg)
 * M = mean mass of one mol of atmospheric particles = 0.029 kg/mol for Earth
 * g = acceleration due to gravity at the current location (m/s2)

The pressure (force per unit area) at a given altitude is a result of the weight of the overlying atmosphere. If at a height of z the atmosphere has density &rho; and pressure P, then moving upwards an infinitesimally small height dz will decrease the pressure by amount dP, equal to the weight of a layer of atmosphere of thickness dz.

Thus: $$\frac{dP}{dz} = -g\rho$$ where g is the acceleration due to gravity. For small dz it is possible to assume g to be constant; the minus sign indicates that as the height increases the pressure decreases. Therefore, using the equation of state for an ideal gas of mean molecular mass M at temperature T, the density can be expressed as $$\rho = \frac{MP}{RT}$$

Combining these equations gives $$\frac{dP}{P} = \frac{-dz}{\frac{k_\text{B}T}{mg}}$$ which can then be incorporated with the equation for H given above to give: $$\frac{dP}{P} = - \frac{dz}{H}$$ which will not change unless the temperature does. Integrating the above and assuming P0 is the pressure at height z = 0 (pressure at sea level) the pressure at height z can be written as: $$P = P_0\exp\left(-\frac{z}{H}\right)$$

This translates as the pressure decreasing exponentially with height.

In Earth's atmosphere, the pressure at sea level P0 averages about $1.01 Pa$, the mean molecular mass of dry air is 28.964 u and hence m = 28.964 &times; $1.66$ = $4.808 kg$. As a function of temperature, the scale height of Earth's atmosphere is therefore H/T = k/mg = (1.38/(4.808&times;9.81))×103 = $29.26 m/K$. This yields the following scale heights for representative air temperatures.


 * T = 290 K, H = 8500 m
 * T = 273 K, H = 8000 m
 * T = 260 K, H = 7610 m
 * T = 210 K, H = 6000 m

These figures should be compared with the temperature and density of Earth's atmosphere plotted at NRLMSISE-00, which shows the air density dropping from 1200 g/m3 at sea level to 0.53 = 0.125 g/m3 at 70 km, a factor of 9600, indicating an average scale height of 70/ln(9600) = 7.64 km, consistent with the indicated average air temperature over that range of close to 260 K.

Note:
 * Density is related to pressure by the ideal gas laws. Therefore, density will also decrease exponentially with height from a sea level value of &rho;0 roughly equal to 1.2 kg m&minus;3
 * At heights over 100 km, an atmosphere may no longer be well mixed. Then each chemical species has its own scale height.
 * Here temperature and gravitational acceleration were assumed to be constant but both may vary over large distances.

Planetary examples
Approximate atmospheric scale heights for selected Solar System bodies:


 * Venus: 15.9 km
 * Earth: 8.5 km
 * Mars: 11.1 km
 * Jupiter: 27 km
 * Saturn: 59.5 km
 * Titan: 21 km
 * Uranus: 27.7 km
 * Neptune: 19.1–20.3 km
 * Pluto: ~50 km

Scale height for a thin disk
For a disk of gas around a condensed central object, such as, for example, a protostar, one can derive a disk scale height which is somewhat analogous to the planetary scale height. We start with a disc of gas that has a mass which is small relative to the central object. We assume that the disc is in hydrostatic equilibrium with the z component of gravity from the star, where the gravity component is pointing to the midplane of the disk:

$$\frac{dP}{dz} = - \frac{GM_*\rho z}{(r^2+z^2)^{3/2}}$$

where:


 * G = Gravitational constant ≈
 * r = the radial cylindrical coordinate for the distance from the center of the star or centrally condensed object
 * z = the height/altitude cylindrical coordinate for the distance from the disk midplane (or center of the star)
 * M* = the mass of the star/centrally condensed object
 * P = the pressure of the gas in the disk
 * $$\rho$$ = the gas mass density in the disk

In the thin disk approximation, $$ z \ll r $$ and the hydrostatic equilibrium equation is $$\frac{dP}{dz} \approx - \frac{GM_*\rho z}{r^3}$$

To determine the gas pressure, one can use the ideal gas law: $$ P = \frac{\rho k_\text{B} T}{\bar{m}} $$ with:
 * T = the gas temperature in the disk, where the temperature is a function of r, but independent of z
 * $$\bar{m}$$ = the mean molecular mass of the gas

Using the ideal gas law and the hydrostatic equilibrium equation, gives: $$\frac{d\rho}{dz} \approx - \frac{GM_* \bar{m}\rho z}{k Tr^3}$$ which has the solution $$ \rho = \rho_0 \exp\left(-\left(\frac{z}{h_D}\right)^2 \right) $$ where $$ \rho_0 $$ is the gas mass density at the midplane of the disk at a distance r from the center of the star and $$ h_D $$ is the disk scale height with

$$ h_D = \sqrt{\frac{2kTr^3}{GM_* \bar{m}}} \approx 0.0306 \sqrt{\frac{\left(T/100 \ \text{K}\right)\left(r/1 \text{ au} \right)^3}{\left(M_* / M_\odot\right)\left(\bar{m}/2 \text{ amu} \right)}} \ \text{ au} $$ with $$ M_\odot $$ the solar mass, $$ \text{au} $$ the astronomical unit and $$ \text{amu} $$ the atomic mass unit.

As an illustrative approximation, if we ignore the radial variation in the temperature, $$ T $$, we see that $$ h_D \propto r^{3/2} $$ and that the disk increases in altitude as one moves radially away from the central object.

Due to the assumption that the gas temperature in the disk, T, is independent of z, $$ h_D $$ is sometimes known as the isothermal disk scale height.

Disk scale height in a magnetic field
A magnetic field in a thin gas disk around a central object can change the scale height of the disk. For example, if a non-perfectly conducting disk is rotating through a poloidal magnetic field (i.e., the initial magnetic field is perpendicular to the plane of the disk), then a toroidal (i.e., parallel to the disk plane) magnetic field will be produced within the disk, which will pinch and compress the disk. In this case, the gas density of the disk is:

$$ \rho(r, z) = \rho_0(r) \exp \left(- \left(\frac{z}{h_D}\right)^2 \right) - \rho_\text{cut}(r) \left[1 - \exp \left(- \left(\frac{z}{h_D}\right)^2 \right) \right] $$ where the cut-off density $$ \rho_\text{cut} $$ has the form $$ \rho_{\rm cut}(r) = (\mu_0 \sigma_D r)^2 \left(\frac{B_z^2}{\mu_0}\right) \left(\frac{\Omega_*}{\Omega_K} - 1 \right)^2 $$ where
 * $$ \mu_0 $$ is the permeability of free space
 * $$ \sigma_D $$ is the electrical conductivity of the disk
 * $$ B_z $$ is the magnetic flux density of the poloidal field in the $$ z $$ direction
 * $$ \Omega_* $$ is the rotational angular velocity of the central object (if the poloidal magnetic field is independent of the central object then $$ \Omega_* $$ can be set to zero)
 * $$ \Omega_K $$ is the keplerian angular velocity of the disk at a distance $$ r $$ from the central object.

These formulae give the maximum height, $$ H_B $$, of the magnetized disk as

$$ H_B = h_D \sqrt{\ln\left(1 + \rho_0/\rho_{\rm cut} \right)} ,$$ while the e-folding magnetic scale height, $$ h_B $$, is $$ h_B = h_D \sqrt{\ln\left(1 + \frac{1 - 1/e}{1/e + \rho_\text{cut}/\rho_0} \right)} \. $$