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).

The earth atmosphere's scale height is about 8.5km, as can be confirmed from this diagram of air pressure p by altitude h: At an altitude of 0, 8.5, and 17 km, the pressure is about 1000, 370, and 140 hPa, respectively.

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^[1][2]

${\displaystyle H={\frac {kT}{mg}}}$ 

or equivalently

${\displaystyle H={\frac {RT}{Mg}}}$ 

where:

• k = Boltzmann constant = 1.38 x 10⁻²³ J·K⁻¹
• R = gas constant
• T = mean atmospheric temperature in kelvins = 250 K^[3] 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/s²)

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 ρ 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:

${\displaystyle {\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

${\displaystyle \rho ={\frac {MP}{RT}}}$ 

Combining these equations gives

${\displaystyle {\frac {dP}{P}}={\frac {-dz}{\frac {kT}{mg}}}}$ 

which can then be incorporated with the equation for H given above to give:

${\displaystyle {\frac {dP}{P}}=-{\frac {dz}{H}}}$ 

which will not change unless the temperature does. Integrating the above and assuming P₀ is the pressure at height z = 0 (pressure at sea level) the pressure at height z can be written as:

${\displaystyle P=P_{0}\exp \left(-{\frac {z}{H}}\right)}$ 

This translates as the pressure decreasing exponentially with height.^[4]

In Earth's atmosphere, the pressure at sea level P₀ averages about 1.01×10⁵ Pa, the mean molecular mass of dry air is 28.964 u and hence m = 28.964 × 1.660×10⁻²⁷ = 4.808×10⁻²⁶ kg. As a function of temperature, the scale height of Earth's atmosphere is therefore H/T = k/mg = (1.38/(4.808×9.81))×10³ = 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/m³ at sea level to 0.5³ = .125 g/m³ 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 ρ₀ roughly equal to 1.2 kg m⁻³
• 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 follow.

• Venus: 15.9 km^[5]
• Earth: 8.5 km^[6]
• Mars: 11.1 km^[7]
• Jupiter: 27 km^[8]
• Saturn: 59.5 km^[9]

• Titan: 21 km^[10]

• Uranus: 27.7 km^[11]
• Neptune: 19.1–20.3 km^[12]
• Pluto: ~50 km^[13]

Scale height for a thin disk

 

A schematic depiction of the force balance in a gas disk around a central object, e.g., a star.

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:

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

where:

• G = Gravitational constant ≈ 6.674 x 10⁻¹¹ m³·kg⁻¹·s⁻²
• 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
• ${\displaystyle \rho }$  = the gas mass density in the disk

In the thin disk approximation, ${\displaystyle z\ll r}$  and the hydrostatic equilibrium equation is

${\displaystyle {\frac {dP}{dz}}\approx -{\frac {GM_{*}\rho z}{r^{3}}}}$ 

To determine the gas pressure, one can use the ideal gas law:

${\displaystyle P={\frac {\rho kT}{\bar {m}}}}$ 

with:

• T = the gas temperature in the disk, where the temperature is a function of r, but independent of z
• ${\displaystyle {\bar {m}}}$  = the mean molecular mass of the gas

Using the ideal gas law and the hydrostatic equilibrium equation, gives:

${\displaystyle {\frac {d\rho }{dz}}\approx -{\frac {GM_{*}{\bar {m}}\rho z}{kTr^{3}}}}$ 

which has the solution

${\displaystyle \rho =\rho _{0}\exp \left(-\left({\frac {z}{h_{D}}}\right)^{2}\right)}$ 

where ${\displaystyle \rho _{0}}$  is the gas mass density at the midplane of the disk at a distance r from the center of the star and ${\displaystyle h_{D}}$  is the disk scale height with

${\displaystyle h_{D}={\sqrt {\frac {2kTr^{3}}{GM_{*}{\bar {m}}}}}\approx 0.0306{\sqrt {\frac {\left(T/100\ {\rm {K}}\right)\left(r/1\ {\rm {au}}\right)^{3}}{\left(M_{*}/{\rm {M}}_{\odot }\right)\left({\bar {m}}/2\ {\rm {amu}}\right)}}}\ \ {\rm {au}}}$ 

with ${\displaystyle {\rm {M}}_{\odot }}$  the solar mass, ${\displaystyle {\rm {au}}}$  the astronomical unit and ${\displaystyle {\rm {amu}}}$  the atomic mass unit.

As an illustrative approximation, if we ignore the radial variation in the temperature, ${\displaystyle T}$ , we see that ${\displaystyle 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, ${\displaystyle 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.^[14][15][16] 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: ^[16]

${\displaystyle \rho (r,z)=\rho _{0}(r)\exp \left(-\left({\frac {z}{h_{D}}}\right)^{2}\right)-\rho _{\rm {cut}}(r)\left[1-\exp \left(-\left({\frac {z}{h_{D}}}\right)^{2}\right)\right]\ }$ 

where the cut-off density ${\displaystyle \rho _{\rm {cut}}}$  has the form

${\displaystyle \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

• ${\displaystyle \mu _{0}}$  is the permeability of free space
• ${\displaystyle \sigma _{D}}$  is the electrical conductivity of the disk
• ${\displaystyle B_{z}}$  is the magnetic flux density of the poloidal field in the ${\displaystyle z}$  direction
• ${\displaystyle \Omega _{*}}$  is the rotational angular velocity of the central object (if the poloidal magnetic field is independent of the central object then ${\displaystyle \Omega _{*}}$  can be set to zero)
• ${\displaystyle \Omega _{K}}$  is the keplerian angular velocity of the disk at a distance ${\displaystyle r}$  from the central object.

These formulae give the maximum height, ${\displaystyle H_{B}}$ , of the magnetized disk as

${\displaystyle H_{B}=h_{D}{\sqrt {\ln \left(1+\rho _{0}/\rho _{\rm {cut}}\right)}}}$ ,

while the e-folding magnetic scale height, ${\displaystyle h_{B}}$ , is

${\displaystyle h_{B}=h_{D}{\sqrt {\ln \left(1+{\frac {1-1/e}{1/e+\rho _{\rm {cut}}/\rho _{0}}}\right)}}\ .}$ 

See also

• Time constant

References

1. "Glossary of Meteorology - scale height". American Meteorological Society (AMS).
2. "Pressure Scale Height". Wolfram Research.
3. "Daniel J. Jacob: "Introduction to Atmospheric Chemistry", Princeton University Press, 1999".
4. (PDF). Archived from the original (PDF) on 2011-07-16.
5. "Venus Fact Sheet". NASA. Retrieved 28 September 2013.
6. "Earth Fact Sheet". NASA. Retrieved 28 September 2013.
7. "Mars Fact Sheet". NASA. Retrieved 28 September 2013.
8. . NASA. Archived from the original on 13 October 2011. Retrieved 28 September 2013.
9. . NASA. Archived from the original on 18 August 2011. Retrieved 28 September 2013.
10. Justus, C. G.; Aleta Duvall; Vernon W. Keller (1 August 2003). "Engineering-Level Model Atmospheres For Titan and Mars". International Workshop on Planetary Probe Atmospheric Entry and Descent Trajectory Analysis and Science, Lisbon, Portugal, October 6–9, 2003, Proceedings: ESA SP-544. ESA. Retrieved 28 September 2013.
11. "Uranus Fact Sheet". NASA. Retrieved 28 September 2013.
12. "Neptune Fact Sheet". NASA. Retrieved 28 September 2013.
13. "Pluto Fact Sheet". NASA. Retrieved 2020-09-28.
14. Lovelace, R.V.E.; Mehanian, C.; Mobarry, C. M.; Sulkanen, M. E. (September 1986). "Theory of Axisymmetric Magnetohydrodynamic Flows: Disks". Astrophysical Journal Supplement. 62: 1. Bibcode:1986ApJS...62....1L. doi:10.1086/191132. Retrieved 26 January 2022.
15. Campbell, C. G.; Heptinstall, P. M. (August 1998). "Disc structure around strongly magnetic accretors: a full disc solution with turbulent diffusivity". Monthly Notices of the Royal Astronomical Society. 299 (1): 31. Bibcode:1998MNRAS.299...31C. doi:10.1046/j.1365-8711.1998.01576.x.
16. Liffman, Kurt; Bardou, Anne (October 1999). "A magnetic scaleheight: the effect of toroidal magnetic fields on the thickness of accretion discs". Monthly Notices of the Royal Astronomical Society. 309 (2): 443. Bibcode:1999MNRAS.309..443L. doi:10.1046/j.1365-8711.1999.02852.x.