Air mass (astronomy)

Last updated

In astronomy, air mass or airmass is a measure of the amount of air along the line of sight when observing a star or other celestial source from below Earth's atmosphere ( Green 1992 ). It is formulated as the integral of air density along the light ray.

Contents

As it penetrates the atmosphere, light is attenuated by scattering and absorption; the thicker atmosphere through which it passes, the greater the attenuation. Consequently, celestial bodies when nearer the horizon appear less bright than when nearer the zenith. This attenuation, known as atmospheric extinction, is described quantitatively by the Beer–Lambert law.

"Air mass" normally indicates relative air mass, the ratio of absolute air masses (as defined above) at oblique incidence relative to that at zenith. So, by definition, the relative air mass at the zenith is 1. Air mass increases as the angle between the source and the zenith increases, reaching a value of approximately 38 at the horizon. Air mass can be less than one at an elevation greater than sea level; however, most closed-form expressions for air mass do not include the effects of the observer's elevation, so adjustment must usually be accomplished by other means.

Tables of air mass have been published by numerous authors, including Bemporad (1904), Allen (1973), [1] and Kasten & Young (1989).

Definition

The absolute air mass is defined as: where is volumetric density of air. Thus is a type of oblique column density.

In the vertical direction, the absolute air mass at zenith is:

So is a type of vertical column density.

Finally, the relative air mass is:

Assuming air density to be uniform allows removing it from the integrals. The absolute air mass then simplifies to a product: where is the average density and the arc length of the oblique and zenith light paths are:

In the corresponding simplified relative air mass, the average density cancels out in the fraction, leading to the ratio of path lengths:

Further simplifications are often made, assuming straight-line propagation (neglecting ray bending), as discussed below.

Calculation

Plots of air mass using various formulas Viewing angle and air masses.svg
Plots of air mass using various formulas

Background

The angle of a celestial body with the zenith is the zenith angle (in astronomy, commonly referred to as the zenith distance ). A body's angular position can also be given in terms of altitude , the angle above the geometric horizon; the altitude and the zenith angle are thus related by

Atmospheric refraction causes light entering the atmosphere to follow an approximately circular path that is slightly longer than the geometric path. Air mass must take into account the longer path ( Young 1994 ). Additionally, refraction causes a celestial body to appear higher above the horizon than it actually is; at the horizon, the difference between the true zenith angle and the apparent zenith angle is approximately 34 minutes of arc. Most air mass formulas are based on the apparent zenith angle, but some are based on the true zenith angle, so it is important to ensure that the correct value is used, especially near the horizon. [2]

Plane-parallel atmosphere

When the zenith angle is small to moderate, a good approximation is given by assuming a homogeneous plane-parallel atmosphere (i.e., one in which density is constant and Earth's curvature is ignored). The air mass then is simply the secant of the zenith angle :

At a zenith angle of 60°, the air mass is approximately 2. However, because the Earth is not flat, this formula is only usable for zenith angles up to about 60° to 75°, depending on accuracy requirements. At greater zenith angles, the accuracy degrades rapidly, with becoming infinite at the horizon; the horizon air mass in the more realistic spherical atmosphere is usually less than 40.

Interpolative formulas

Many formulas have been developed to fit tabular values of air mass; one by Young & Irvine (1967) included a simple corrective term: where is the true zenith angle. This gives usable results up to approximately 80°, but the accuracy degrades rapidly at greater zenith angles. The calculated air mass reaches a maximum of 11.13 at 86.6°, becomes zero at 88°, and approaches negative infinity at the horizon. The plot of this formula on the accompanying graph includes a correction for atmospheric refraction so that the calculated air mass is for apparent rather than true zenith angle.

Hardie (1962) introduced a polynomial in : which gives usable results for zenith angles of up to perhaps 85°. As with the previous formula, the calculated air mass reaches a maximum, and then approaches negative infinity at the horizon.

Rozenberg (1966) suggested which gives reasonable results for high zenith angles, with a horizon air mass of 40.

Kasten & Young (1989) developed [3]

which gives reasonable results for zenith angles of up to 90°, with an air mass of approximately 38 at the horizon. Here the second term is in degrees.

Young (1994) developed in terms of the true zenith angle , for which he claimed a maximum error (at the horizon) of 0.0037 air mass.

Pickering (2002) developed where is apparent altitude in degrees. Pickering claimed his equation to have a tenth the error of Schaefer (1998) near the horizon. [4]

Atmospheric models

Interpolative formulas attempt to provide a good fit to tabular values of air mass using minimal computational overhead. The tabular values, however, must be determined from measurements or atmospheric models that derive from geometrical and physical considerations of Earth and its atmosphere.

Nonrefracting spherical atmosphere

Atmospheric effects on optical transmission can be modelled as if the atmosphere is concentrated in approximately the lower 9 km. Airmass geometry.png
Atmospheric effects on optical transmission can be modelled as if the atmosphere is concentrated in approximately the lower 9 km.

If atmospheric refraction is ignored, it can be shown from simple geometrical considerations (Schoenberg 1929, 173) that the path of a light ray at zenith angle through a radially symmetrical atmosphere of height above the Earth is given by or alternatively, where is the radius of the Earth.

The relative air mass is then:

Homogeneous atmosphere

If the atmosphere is homogeneous (i.e., density is constant), the atmospheric height follows from hydrostatic considerations as:[ citation needed ] where is the Boltzmann constant, is the sea-level temperature, is the molecular mass of air, and is the acceleration due to gravity. Although this is the same as the pressure scale height of an isothermal atmosphere, the implication is slightly different. In an isothermal atmosphere, 37% (1/e) of the atmosphere is above the pressure scale height; in a homogeneous atmosphere, there is no atmosphere above the atmospheric height.

Taking , , and gives . Using Earth's mean radius of 6371 km, the sea-level air mass at the horizon is

The homogeneous spherical model slightly underestimates the rate of increase in air mass near the horizon; a reasonable overall fit to values determined from more rigorous models can be had by setting the air mass to match a value at a zenith angle less than 90°. The air mass equation can be rearranged to give matching Bemporad's value of 19.787 at  = 88° gives   631.01 and   35.54. With the same value for as above,   10,096 m.

While a homogeneous atmosphere is not a physically realistic model, the approximation is reasonable as long as the scale height of the atmosphere is small compared to the radius of the planet. The model is usable (i.e., it does not diverge or go to zero) at all zenith angles, including those greater than 90° (see § Homogeneous spherical atmosphere with elevated observer ). The model requires comparatively little computational overhead, and if high accuracy is not required, it gives reasonable results. [5] However, for zenith angles less than 90°, a better fit to accepted values of air mass can be had with several of the interpolative formulas.

Variable-density atmosphere

In a real atmosphere, density is not constant (it decreases with elevation above mean sea level. The absolute air mass for the geometrical light path discussed above, becomes, for a sea-level observer,

Isothermal atmosphere

Several basic models for density variation with elevation are commonly used. The simplest, an isothermal atmosphere, gives where is the sea-level density and is the density scale height. When the limits of integration are zero and infinity, the result is known as Chapman function. An approximate result is obtained if some high-order terms are dropped, yielding ( Young 1974 , p. 147),

An approximate correction for refraction can be made by taking ( Young 1974 , p. 147) where is the physical radius of the Earth. At the horizon, the approximate equation becomes

Using a scale height of 8435 m, Earth's mean radius of 6371 km, and including the correction for refraction,

Polytropic atmosphere

The assumption of constant temperature is simplistic; a more realistic model is the polytropic atmosphere, for which where is the sea-level temperature and is the temperature lapse rate. The density as a function of elevation is where is the polytropic exponent (or polytropic index). The air mass integral for the polytropic model does not lend itself to a closed-form solution except at the zenith, so the integration usually is performed numerically.

Layered atmosphere

Earth's atmosphere consists of multiple layers with different temperature and density characteristics; common atmospheric models include the International Standard Atmosphere and the US Standard Atmosphere. A good approximation for many purposes is a polytropic troposphere of 11 km height with a lapse rate of 6.5 K/km and an isothermal stratosphere of infinite height ( Garfinkel 1967 ), which corresponds very closely to the first two layers of the International Standard Atmosphere. More layers can be used if greater accuracy is required. [6]

Refracting radially symmetrical atmosphere

When atmospheric refraction is considered, ray tracing becomes necessary ( Kivalov 2007 ), and the absolute air mass integral becomes [7] where is the index of refraction of air at the observer's elevation above sea level, is the index of refraction at elevation above sea level, , is the distance from the center of the Earth to a point at elevation , and is distance to the upper limit of the atmosphere at elevation . The index of refraction in terms of density is usually given to sufficient accuracy (Garfinkel 1967) by the Gladstone–Dale relation

Rearrangement and substitution into the absolute air mass integral gives

The quantity is quite small; expanding the first term in parentheses, rearranging several times, and ignoring terms in after each rearrangement, gives ( Kasten & Young 1989 )

Homogeneous spherical atmosphere with elevated observer

Air mass for elevated observer in homogeneous spherical atmosphere HomogSphElevObsAM.png
Air mass for elevated observer in homogeneous spherical atmosphere

In the figure at right, an observer at O is at an elevation above sea level in a uniform radially symmetrical atmosphere of height . The path length of a light ray at zenith angle is ; is the radius of the Earth. Applying the law of cosines to triangle OAC, expanding the left- and right-hand sides, eliminating the common terms, and rearranging gives

Solving the quadratic for the path length s, factoring, and rearranging,

The negative sign of the radical gives a negative result, which is not physically meaningful. Using the positive sign, dividing by , and cancelling common terms and rearranging gives the relative air mass:

With the substitutions and , this can be given as

When the observer's elevation is zero, the air mass equation simplifies to

In the limit of grazing incidence, the absolute airmass equals the distance to the horizon. Furthermore, if the observer is elevated, the horizon zenith angle can be greater than 90°.

Maximum zenith angle for elevated observer in homogeneous spherical atmosphere HomogSphElevObsZmax.png
Maximum zenith angle for elevated observer in homogeneous spherical atmosphere

Nonuniform distribution of attenuating species

Atmospheric models that derive from hydrostatic considerations assume an atmosphere of constant composition and a single mechanism of extinction, which isn't quite correct. There are three main sources of attenuation ( Hayes & Latham 1975 ): Rayleigh scattering by air molecules, Mie scattering by aerosols, and molecular absorption (primarily by ozone). The relative contribution of each source varies with elevation above sea level, and the concentrations of aerosols and ozone cannot be derived simply from hydrostatic considerations.

Rigorously, when the extinction coefficient depends on elevation, it must be determined as part of the air mass integral, as described by Thomason, Herman & Reagan (1983). A compromise approach often is possible, however. Methods for separately calculating the extinction from each species using closed-form expressions are described in Schaefer (1993) and Schaefer (1998). The latter reference includes source code for a BASIC program to perform the calculations. Reasonably accurate calculation of extinction can sometimes be done by using one of the simple air mass formulas and separately determining extinction coefficients for each of the attenuating species (Green 1992, Pickering 2002).

Implications

Air mass and astronomy

Atmospheric transmittance across the electromagnetic spectrum Atmospheric electromagnetic transmittance or opacity.jpg
Atmospheric transmittance across the electromagnetic spectrum

In optical astronomy, the air mass provides an indication of the deterioration of the observed image, not only as regards direct effects of spectral absorption, scattering and reduced brightness, but also an aggregation of visual aberrations, e.g. resulting from atmospheric turbulence, collectively referred to as the quality of the "seeing". [8] On bigger telescopes, such as the WHT ( Wynne & Worswick 1988 ) and VLT ( Avila, Rupprecht & Beckers 1997 ), the atmospheric dispersion can be so severe that it affects the pointing of the telescope to the target. In such cases an atmospheric dispersion compensator is used, which usually consists of two prisms.

The Greenwood frequency and Fried parameter, both relevant for adaptive optics, depend on the air mass above them (or more specifically, on the zenith angle).

In radio astronomy the air mass (which influences the optical path length) is not relevant. The lower layers of the atmosphere, modeled by the air mass, do not significantly impede radio waves, which are of much lower frequency than optical waves. Instead, some radio waves are affected by the ionosphere in the upper atmosphere. Newer aperture synthesis radio telescopes are especially affected by this as they “see” a much larger portion of the sky and thus the ionosphere. In fact, LOFAR needs to explicitly calibrate for these distorting effects (van der Tol & van der Veen 2007; de Vos, Gunst & Nijboer 2009), but on the other hand can also study the ionosphere by instead measuring these distortions ( Thidé 2007 ).

Air mass and solar energy

Solar irradiance spectrum above atmosphere and at surface Solar Spectrum.png
Solar irradiance spectrum above atmosphere and at surface

In some fields, such as solar energy and photovoltaics, air mass is indicated by the acronym AM; additionally, the value of the air mass is often given by appending its value to AM, so that AM1 indicates an air mass of 1, AM2 indicates an air mass of 2, and so on. The region above Earth's atmosphere, where there is no atmospheric attenuation of solar radiation, is considered to have "air mass zero" (AM0).

Atmospheric attenuation of solar radiation is not the same for all wavelengths; consequently, passage through the atmosphere not only reduces intensity but also alters the spectral irradiance. Photovoltaic modules are commonly rated using spectral irradiance for an air mass of 1.5 (AM1.5); tables of these standard spectra are given in ASTM G 173-03. The extraterrestrial spectral irradiance (i.e., that for AM0) is given in ASTM E 490-00a. [9]

For many solar energy applications when high accuracy near the horizon is not required, air mass is commonly determined using the simple secant formula described in § Plane-parallel atmosphere .

See also

Notes

  1. Allen's air mass table was an abbreviated compilation of values from earlier sources, primarily Bemporad (1904).
  2. At very high zenith angles, air mass is strongly dependent on local atmospheric conditions, including temperature, pressure, and especially the temperature gradient near the ground. In addition low-altitude extinction is strongly affected by the aerosol concentration and its vertical distribution. Many authors have cautioned that accurate calculation of air mass near the horizon is all but impossible.
  3. The Kasten and Young formula was originally given in terms of altitude as in this article, it is given in terms of zenith angle for consistency with the other formulas.
  4. Pickering (2002) uses Garfinkel (1967) as the reference for accuracy. (p. 22, continuation of footnote 39 beginning on p. 21)
  5. Although acknowledging that an isothermal or polytropic atmosphere would have been more realistic, Janiczek & DeYoung (1987) used the homogeneous spherical model in calculating illumination from the Sun and Moon, with the implication that the slightly reduced accuracy was more than offset by the considerable reduction in computational overhead.
  6. The notes for Reed Meyer's air mass calculator describe an atmospheric model using eight layers and using polynomials rather than simple linear relations for temperature lapse rates.
  7. See Thomason, Herman & Reagan (1983) for a derivation of the integral for a refracting atmosphere.
  8. Observing tips: air mass and differential refraction retrieved 15 May 2011.
  9. ASTM E 490-00a was reapproved without change in 2006.

Related Research Articles

<span class="mw-page-title-main">Centripetal force</span> Force directed to the center of rotation

A centripetal force is a force that makes a body follow a curved path. The direction of the centripetal force is always orthogonal to the motion of the body and towards the fixed point of the instantaneous center of curvature of the path. Isaac Newton described it as "a force by which bodies are drawn or impelled, or in any way tend, towards a point as to a centre". In Newtonian mechanics, gravity provides the centripetal force causing astronomical orbits.

<span class="mw-page-title-main">Spherical coordinate system</span> Coordinates comprising a distance and two angles

In mathematics, a spherical coordinate system is a coordinate system for three-dimensional space where the position of a given point in space is specified by three real numbers: the radial distancer along the radial line connecting the point to the fixed point of origin; the polar angleθ between the radial line and a polar axis; and the azimuthal angleφ as the angle of rotation of the radial line around the polar axis. (See graphic re the "physics convention".) Once the radius is fixed, the three coordinates (r, θ, φ), known as a 3-tuple, provide a coordinate system on a sphere, typically called the spherical polar coordinates.

<span class="mw-page-title-main">Relative density</span> Ratio of two densities

Relative density, also called specific gravity, is a dimensionless quantity defined as the ratio of the density of a substance to the density of a given reference material. Specific gravity for solids and liquids is nearly always measured with respect to water at its densest ; for gases, the reference is air at room temperature. The term "relative density" is preferred in SI, whereas the term "specific gravity" is gradually being abandoned.

<span class="mw-page-title-main">Navier–Stokes equations</span> Equations describing the motion of viscous fluid substances

The Navier–Stokes equations are partial differential equations which describe the motion of viscous fluid substances. They were named after French engineer and physicist Claude-Louis Navier and the Irish physicist and mathematician George Gabriel Stokes. They were developed over several decades of progressively building the theories, from 1822 (Navier) to 1842–1850 (Stokes).

<span class="mw-page-title-main">Cylindrical coordinate system</span> 3-dimensional coordinate system

A cylindrical coordinate system is a three-dimensional coordinate system that specifies point positions by the distance from a chosen reference axis (axis L in the image opposite), the direction from the axis relative to a chosen reference direction (axis A), and the distance from a chosen reference plane perpendicular to the axis (plane containing the purple section). The latter distance is given as a positive or negative number depending on which side of the reference plane faces the point.

In atmospheric science, geostrophic flow is the theoretical wind that would result from an exact balance between the Coriolis force and the pressure gradient force. This condition is called geostrophic equilibrium or geostrophic balance. The geostrophic wind is directed parallel to isobars. This balance seldom holds exactly in nature. The true wind almost always differs from the geostrophic wind due to other forces such as friction from the ground. Thus, the actual wind would equal the geostrophic wind only if there were no friction and the isobars were perfectly straight. Despite this, much of the atmosphere outside the tropics is close to geostrophic flow much of the time and it is a valuable first approximation. Geostrophic flow in air or water is a zero-frequency inertial wave.

<span class="mw-page-title-main">Projectile motion</span> Motion of launched objects due to gravity

Projectile motion is a form of motion experienced by an object or particle that is projected in a gravitational field, such as from Earth's surface, and moves along a curved path under the action of gravity only. In the particular case of projectile motion on Earth, most calculations assume the effects of air resistance are passive.

<span class="mw-page-title-main">Stokes parameters</span> Set of values that describe the polarization state of electromagnetic radiation

The Stokes parameters are a set of values that describe the polarization state of electromagnetic radiation. They were defined by George Gabriel Stokes in 1852, as a mathematically convenient alternative to the more common description of incoherent or partially polarized radiation in terms of its total intensity (I), (fractional) degree of polarization (p), and the shape parameters of the polarization ellipse. The effect of an optical system on the polarization of light can be determined by constructing the Stokes vector for the input light and applying Mueller calculus, to obtain the Stokes vector of the light leaving the system. They can be determined from directly observable phenomena. The original Stokes paper was discovered independently by Francis Perrin in 1942 and by Subrahamanyan Chandrasekhar in 1947, who named it as the Stokes parameters.

<span class="mw-page-title-main">Euler–Bernoulli beam theory</span> Method for load calculation in construction

Euler–Bernoulli beam theory is a simplification of the linear theory of elasticity which provides a means of calculating the load-carrying and deflection characteristics of beams. It covers the case corresponding to small deflections of a beam that is subjected to lateral loads only. By ignoring the effects of shear deformation and rotatory inertia, it is thus a special case of Timoshenko–Ehrenfest beam theory. It was first enunciated circa 1750, but was not applied on a large scale until the development of the Eiffel Tower and the Ferris wheel in the late 19th century. Following these successful demonstrations, it quickly became a cornerstone of engineering and an enabler of the Second Industrial Revolution.

<span class="mw-page-title-main">Chapman function</span>

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

<span class="mw-page-title-main">Multiple integral</span> Generalization of definite integrals to functions of multiple variables

In mathematics (specifically multivariable calculus), a multiple integral is a definite integral of a function of several real variables, for instance, f(x, y) or f(x, y, z).

<span class="mw-page-title-main">Toroidal coordinates</span>

Toroidal coordinates are a three-dimensional orthogonal coordinate system that results from rotating the two-dimensional bipolar coordinate system about the axis that separates its two foci. Thus, the two foci and in bipolar coordinates become a ring of radius in the plane of the toroidal coordinate system; the -axis is the axis of rotation. The focal ring is also known as the reference circle.

<span class="mw-page-title-main">Shallow water equations</span> Set of partial differential equations that describe the flow below a pressure surface in a fluid

The shallow-water equations (SWE) are a set of hyperbolic partial differential equations that describe the flow below a pressure surface in a fluid. The shallow-water equations in unidirectional form are also called (de) Saint-Venant equations, after Adhémar Jean Claude Barré de Saint-Venant.

In mathematics, the cylindrical harmonics are a set of linearly independent functions that are solutions to Laplace's differential equation, , expressed in cylindrical coordinates, ρ (radial coordinate), φ (polar angle), and z (height). Each function Vn(k) is the product of three terms, each depending on one coordinate alone. The ρ-dependent term is given by Bessel functions (which occasionally are also called cylindrical harmonics).

The air mass coefficient defines the direct optical path length through the Earth's atmosphere, expressed as a ratio relative to the path length vertically upwards, i.e. at the zenith. The air mass coefficient can be used to help characterize the solar spectrum after solar radiation has traveled through the atmosphere.

<span class="mw-page-title-main">Vacuum airship</span> Hypothetical airship that is evacuated rather than filled with a lighter-than-air gas

A vacuum airship, also known as a vacuum balloon, is a hypothetical airship that is evacuated rather than filled with a lighter-than-air gas such as hydrogen or helium. First proposed by Italian Jesuit priest Francesco Lana de Terzi in 1670, the vacuum balloon would be the ultimate expression of lifting power per volume displaced.

Blade element momentum theory is a theory that combines both blade element theory and momentum theory. It is used to calculate the local forces on a propeller or wind-turbine blade. Blade element theory is combined with momentum theory to alleviate some of the difficulties in calculating the induced velocities at the rotor.

An axial fan is a type of fan that causes gas to flow through it in an axial direction, parallel to the shaft about which the blades rotate. The flow is axial at entry and exit. The fan is designed to produce a pressure difference, and hence force, to cause a flow through the fan. Factors which determine the performance of the fan include the number and shape of the blades. Fans have many applications including in wind tunnels and cooling towers. Design parameters include power, flow rate, pressure rise and efficiency.

The planar reentry equations are the equations of motion governing the unpowered reentry of a spacecraft, based on the assumptions of planar motion and constant mass, in an Earth-fixed reference frame.

References