## 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