Test of Math Mode
Actually, both math and LaTeX2html. Taken from
RefBendDelayCalc but with the figures left out.
Another test:
and the
old way.
The ALMA Delay Server
The ALMA delay server is a software module which is responsible for distributing the delay correction for each antenna to the various bits of hardware in the antenna and to the correlator. The delay correction is defined as:
where
is some arbitrary offset to make sure all delays are positive,
is a per antenna constant representing the electronic delay,
is the geometric delay derived by making a call to CALC, and
is the atmospheric delay correction (which can be gotten from CALC but is (as of 2008/03/18) not used and set to zero. In the long term the atmospheric delay will be derived using ATM to calculate the zenith delay coupled to an appropriate mapping function.
Overview
The delay experienced by an incoming signal due to its propagation through the Earth's atmosphere is given by:
where
is the path through and n is the refractive index of the atmosphere. Since
is very nearly unity for the Earth's atmosphere one normally uses the "refractivity" (
) instead of the index of refraction. Refractivity and refractive index are related as follows:
Furthermore, the atmospheric delay can be separated into contributions due to the dry and wet atmosphere:
where
is the contribution due to dry air while
is the contribution due to wet air. In general
and
are parameterized in terms of a zenith contribution to the delay which is dependent upon local atmospheric conditions (
) and a "mapping function" (
) which relate delays at an arbitrary elevation angle
to that at the zenith:
Since the elevation angle
is the unrefracted source elevation, refraction effects are included in the mapping functions
. In the following I describe calculations of
and
.
Zenith Delay
The contribution to the atmospheric delay at the zenith (
) is a measure of the integrated refractivity of the atmosphere at the zenith. In general, the refractivity of moist air at microwave frequencies depends upon the permanent and induced dipole moments of the molecular species that make up the atmosphere. The primary species that make up the dry atmosphere, nitrogen and oxygen, do not have permanent dipole moments, so contribute to the refractivity via their induced dipole moments. Water vapour does have a permanent dipole moment. Permanent dipole moments contribute to the refractivity as
, while induced dipole moments contribute as
, where
is the pressure and
is the temperature of the species.
A simple parameterization of the refractivity at the zenith is given by the Smith-Weintraub equation (Smith & Weintraub 1953):
where
and
are the partial pressures due to dry and wet air, T is the temperature of the atmosphere, and
,
, and
are constants. The pressures and temperature are usually taken to be equivalent to the local ambient values at the antenna station under study. The dry and wet air refractivities are then given by:
The dry air contribution to this refractivity (
) is primarily due to oxygen and nitrogen, and is nearly in hydrostatic equilibrium. Therefore,
does not depend upon the detailed behaviour of dry air pressure and temperature along the path through the atmosphere, and can be derived based on local atmospheric temperature and pressure measurements. The wet air refractivity (
) can be inferred from local water vapour radiometry measurements. Also, one can use an atmospheric model to derive the total atmospheric refractivity using an atmospheric model (such as ATM) with local atmospheric conditions as input.
Mapping Function
The simplest form for the mapping function (
), which relates the delay at an arbitrary elevation angle
to that at the zenith, is given by the plane-parallel approximation for the Earth's atmosphere:
This simple form is in fact inadequate, which led Marini (1972) to consider corrections to this simple functional form which accounted for the Earth's curvature. Assuming an exponential atmospheric profile where the atmospheric refractivity varies exponentially with height above the antenna, Marini (1972) developed a continued fraction form for the mapping function:
where I include only the first three terms in the continued fraction. A slight modification to the Marini (1972) continued fraction functional form which forces
at the zenith replaces the even numbered
terms (i.e. the second, fourth, sixth, etc.) with
Chao (1974) introduced this modification by truncating the Marini (1972) form to include only two terms.
A more generalized continued-fractional form for the mapping function was developed by Yan and Ping (1995):
...where...
...is the "normalized effective zenith argument" of function which includes the "normalized effective height" of the atmosphere (
) defined as:
For an exponentially-distributed atmosphere:
Normally,
(i.e. start the integration from the ground).
The constants
, etc. in the continued fraction forms above are generally derived from analytic fits to ray-tracing results either for standard atmospheres or for observed atmospheric profiles based on radiosonde measurements. The mapping functions derived in Niell (1996) and Davis (1985) are derived in this way.
A physically more correct mapping function has been derived by Lanyi (1984). Unlike previous mapping functions, Lanyi does not fully separate the dry and wet contributions to the delay, which is a more physically correct approximation. It is based on an ideal model atmosphere whose temperature is constant from the surface to the inversion layer
, then decreases linearly with height at rate
(the lapse rate) from
to the tropopause height
, then is assumed to be constant above
. This mapping function is designed then to be a semi-analytic approximation to the atmospheric delay integral that retains an explicit temperature profile that can be determined using meteorological measurements. The mapping function is expanded as a second-order polynomial in
and
, plus the largest third-order term). It is nonlinear in
and
. It also contains terms which couple
and
, thus including terms which arise from the bending of the signal path through the atmosphere. The functional form for the atmospheric delay in this Lanyi (1984) model is given by:
where
where
With standard values of
,
,
(appropriate for mid-latitudes), and
,
.
The dry, wet, and bending contributions are expressed in terms of moments of the refractivity. The bending terms are evaluated for the ideal model atmosphere and thus give the dependence of the delay on the four parameters
,
,
, and
. Therefore, the Lanyi (1984) model relies upon accurate surface meteorological measurements at the time of the observations to which the delay model is applied.
Antenna Height Correction to Total Atmospheric Delay
In the calculation of the zenith atmospheric delay at an antenna it is assumed that the atmospheric properties (P, T, RH) are the values measured at the focal plane of the antenna. For example, in VLBI each station has a set of associated weather measurements which are used to calculate
. For a clustered array like the VLA or ALMA, the affect of the differences in antenna focal plane height above some reference point need to be accounted for.
For the VLA (not
EVLA), CALC was not used to calculate the atmospheric delay. The antenna height correction was incorporated with a simple atmospheric delay correction by correcting for the path difference between each VLA antenna and a reference point at the center of the array. For the VLA case, the extra atmospheric path due to a difference in antenna height above the center-of-the-array reference point (
, in ns) is given by:
where
is the atmospheric refractivity,
is the atmospheric zenith delay calculated using the VLA weather station (which is located near the center of the array),
is the geometric w of the antenna (in ns). The first term is the antenna height correction to the zenith delay, while the second term is a simple atmospheric delay correction. For
EVLA, CALC will be used to calculate both geometric and atmospheric delay. We believe (though have not confirmed) that CALC also calculates the antenna height correction (first term in the equation above) given antenna heights relative to the reference point at the center of the array. ALMA will need to include this antenna height correction term.
A simple estimate of the magnitude of the antenna height difference correction at the zenith can be gotten by assuming that the pressure P changes linearly with height. Then 52 cm of additional antenna height (the current difference in height between the two ATF antennas) out of a total atmospheric height of 8 km would correspond to:
where I have assumed P = 1053 mb. The dry term zenith atmospheric delay changes approximately like 2.3 mm/mb of pressure change. A pressure change of 0.068 mb corresponds to approximately 156 micron of path difference. This is consistent with alternate back-of-the-envelope calculations of this quantity.
Differential Excess Atmospheric Delay Between Two Antennas
NOTE: The following is just an aside. Since CALC or any other analysis of the atmospheric delay at an antenna calculates the total integrated delay along the path of observation, the differential delay between two antennas is accounted for in any differencing calculations done during baseline determination.
The differential delay induced in an interferometer by a horizontally stratified troposphere results from the difference in zenith angle of the source at the antennas. Thompson, Moran, and Swensen (2001), pp. 516-518 discuss the atmospheric delay induced along an interferometer baseline. The excess path length is given by:
where
is the refractivity at the Earth's surface,
is the height above the Earth's surface,
is the atmospheric scale height,
is the length coordinate along the direction to the source, and E is the antenna elevation while observing the source. Note that refraction is neglected. One can relate
,
,
, and
as follows (see Figure 13.4 in TMS, page 517; reproduced below) using the cosine rule on the triangle formed by
,
, and
:
Solving for
and using elevation rather than zenith angle yields:
For the triangle which is formed by sides
,
, and the side which is equal to
, we can write:
Since
and
(the height of the troposphere),
. Since
,
. The equation for
in terms of
,
, and
then becomes:
(Thanks to Dick Thompson for filling-in some of the details of this calculation).
We can now write the expression for
as follows:
Since
, the second term in the equation above can be expanded with a Taylor series so that:
Integration yields:
Writing this equation in terms involving
, the excess path length
becomes:
Taking the derivative of
with respect to E and multiplying this derivative by the baseline length
divided by
yields the atmospheric differential delay between two antennas separated by baseline
:
...where
is in m,
is in km, and
is in km. Note that once must calculate
using a suitable atmospheric model which uses measurements of the local atmospheric pressure, temperature and relative humidity to derive the resultant differential residual delay.
--
PatrickMurphy - 2009-01-21