Chapman function

A Chapman function describes the integration of atmospheric absorption along a slant path on a spherical Earth, relative to the vertical case. It applies to any quantity with a concentration decreasing exponentially with increasing altitude. To a first approximation, valid at small zenith angles, the Chapman function for optical absorption is equal to

Graph of ch(x, z)

$\sec(z),\$

where z is the zenith angle and sec denotes the secant function.

The Chapman function is named after Sydney Chapman, who introduced the function in 1931.

Definition

In an isothermal model of the atmosphere, the density ${\textstyle \varrho (h)}$ varies exponentially with altitude ${\textstyle h}$ according to the Barometric formula:

$\varrho (h)=\varrho _{0}\exp \left(-{\frac {h}{H}}\right)$,

where ${\textstyle \varrho _{0}}$ denotes the density at sea level (${\textstyle h=0}$) and ${\textstyle H}$ the so-called scale height. The total amount of matter traversed by a vertical ray starting at altitude ${\textstyle h}$ towards infinity is given by the integrated density ("column depth")

$X_{0}(h)=\int _{h}^{\infty }\varrho (l)\,\mathrm {d} l=\varrho _{0}H\exp \left(-{\frac {h}{H}}\right)$.

For inclined rays having a zenith angle ${\textstyle z}$, the integration is not straight-forward due to the non-linear relationship between altitude and path length when considering the curvature of Earth. Here, the integral reads

$X_{z}(h)=\varrho _{0}\exp \left(-{\frac {h}{H}}\right)\int _{0}^{\infty }\exp \left(-{\frac {1}{H}}\left({\sqrt {s^{2}+l^{2}+2ls\cos z}}-s\right)\right)\,\mathrm {d} l$,

where we defined ${\textstyle s=h+R_{\mathrm {E} }}$ (${\textstyle R_{\mathrm {E} }}$ denotes the Earth radius).

The Chapman function ${\textstyle \operatorname {ch} (x,z)}$ is defined as the ratio between slant depth ${\textstyle X_{z}}$ and vertical column depth ${\textstyle X_{0}}$. Defining ${\textstyle x=s/H}$, it can be written as

$\operatorname {ch} (x,z)={\frac {X_{z}}{X_{0}}}=\mathrm {e} ^{x}\int _{0}^{\infty }\exp \left(-{\sqrt {x^{2}+u^{2}+2xu\cos z}}\right)\,\mathrm {d} u$.

Representations

A number of different integral representations have been developed in the literature. Chapman's original representation reads

$\operatorname {ch} (x,z)=x\sin z\int _{0}^{z}{\frac {\exp \left(x(1-\sin z/\sin \lambda )\right)}{\sin ^{2}\lambda }}\,\mathrm {d} \lambda$.

Huestis developed the representation

$\operatorname {ch} (x,z)=1+x\sin z\int _{0}^{z}{\frac {\exp \left(x(1-\sin z/\sin \lambda )\right)}{1+\cos \lambda }}\,\mathrm {d} \lambda$,

which does not suffer from numerical singularities present in Chapman's representation.

Special cases

For ${\textstyle z=\pi /2}$ (horizontal incidence), the Chapman function reduces to

$\operatorname {ch} \left(x,{\frac {\pi }{2}}\right)=x\mathrm {e} ^{x}K_{1}(x)$.

Here, ${\textstyle K_{1}(x)}$ refers to the modified Bessel function of the second kind of the first order. For large values of ${\textstyle x}$, this can further be approximated by

$\operatorname {ch} \left(x\gg 1,{\frac {\pi }{2}}\right)\approx {\sqrt {{\frac {\pi }{2}}x}}$.

For ${\textstyle x\rightarrow \infty }$ and ${\textstyle 0\leq z<\pi /2}$, the Chapman function converges to the secant function:

$\lim _{x\rightarrow \infty }\operatorname {ch} (x,z)=\sec z$.

In practical applications related to the terrestrial atmosphere, where ${\textstyle x\sim 1000}$, ${\textstyle \operatorname {ch} (x,z)\approx \sec z}$ is a good approximation for zenith angles up to 60° to 70°, depending on the accuracy required.

See also

• Air mass
• Atmospheric physics
• Ionosphere

External links

• Chapman function at Science World
• Smith, F. L.; Smith, Cody (1972). "Numerical evaluation of Chapman's grazing incidence integral ch(X,χ)". J. Geophys. Res. 77 (19): 3592–3597. Bibcode:1972JGR....77.3592S. doi:10.1029/JA077i019p03592.