Data Differentiations (GPS Observation Equations and Equivalence Properties)

Data differentiations are methods of combining GPS data (of the same type) measured at different stations. For the convenience of later discussions, tidal effects and relativ-istic effects are considered corrected before forming the differences. The original code, phase and Doppler observables as well as their standardised combinations can be rewritten as (cf. Eqs. 6.44-6.47)

tmpD-83_thumbtmpD-84_thumb

wheretmpD-85_thumbis the index of frequency f, subscript i is the index of station number and superscript k is the id number of satellite.


Single Differences

Single difference (SD) is the difference formed by data observed at two stations on the same satellite as

tmpD-87_thumb

where O is the original observable, andtmpD-88_thumbare two id number of the stations.

Supposing the original observables have the same variance oftmpD-89_thumbthen the single difference observable has a variance oftmpD-90_thumbConsidering Eqs. 6.89-6.92, one has

tmpD-94_thumb

wheretmpD-95_thumbis the time differentiation oftmpD-96_thumbare the differenced ionospheric and tropospheric effects at the two stations related to the satellite k , respectively.

The most important property of single differences is that the satellite clock error terms in the model are eliminated. However, it should be emphasised that the satellite clock error, which implicitly affects the computation of satellite position, still has to be carefully considered. Ionospheric and tropospheric effects are reduced through difference forming, especially for those stations that are not very far away from each other. Because of the identical mathematical models of the station clock errors and ambiguities, not all clock and ambiguity parameters can be resolved in the single difference equations of Eqs. 6.94-6.96.

For the original observable vector of station i1 and i2,

tmpD-99_thumb

the single differences

tmpD-100_thumb

can be formed by a linear transformation

tmpD-101_thumb

Where common satellitestmpD-102_thumbare observed, E is an identity matrix that has the size of the observed satellite number; in the above example the size is 3 x 3. The covariance matrix of the single differences is then

tmpD-104_thumb

i.e., the weight matrix is

tmpD-105_thumb

That is, the single differences are un-correlated observables in the case of a single baseline. C in Eq. 6.97 is a general form, so C is denoted bytmpD-106_thumbis the number of commonly viewed satellites.

Single differences can be formed for any baselines as long as the two stations have common satellites in sight. However, the baselines should be a set of "independent" ones. The most-used methods are to form the radial baselines or traverse baselines. Supposing the stations’ id vector istmpD-107_thumband the baseline between station i1 and i2 is denoted bytmpD-108_thumbthen the radial baselines could be formed, e.g., bytmpD-109_thumb

tmpD-110_thumband the traverse baselines could be formed, e.g., bytmpD-111_thumb tmpD-112_thumbStation i1 is called a reference station and is freely selectable. In some cases, a mixed radial and traverse baselines have to be formed such as, e.g., bytmpD-113_thumb

tmpD-114_thumbSometimes the baselines have to be formed by several groups, and therefore several references have to be selected. A method of forming an independent and optimal baseline network will be discussed Sects. 9.1 and 9.2.

In case three stations are used to measure the GPS data, the original observable vector of station i1, i2 and i3 is

tmpD-124_thumb

where n is the commonly observed satellite number. The single differences of the baseline (i, j) are

tmpD-125_thumb

If the baselines are formed in a radial way, i.e., baselines are formed astmpD-126_thumband tmpD-127_thumbthen one has

tmpD-130_thumbtmpD-131_thumb

If the baselines are formed in a traverse way, i.e., baselines are formed astmpD-132_thumb andtmpD-133_thumbthen one has

tmpD-136_thumb

It is obvious that the single differences are correlated if the station numbers are more than two. And the correlation depends on the ways the baselines are formed. Therefore, a general covariance formula of the single differences of a network is not possible to be derived. Furthermore, the commonly viewed satellite number n could be different from baseline to baseline, so the formulation of the covariance matrix could be more complicated.

A baseline-wise processing of the GPS data of a network by using single differences is equivalent to an omission of the correlation between the baselines.

Double Differences

Double differences are formed between two single differences related to two observed satellites as

tmpD-137_thumb

where k1 and k2 are the two id numbers of the satellites. Supposing the original observables have the same variance oftmpD-138_thumbthen the double differenced observables have a variance oftmpD-139_thumbConsidering Eqs. 6.89-6.92, one has

tmpD-142_thumb

wheretmpD-143_thumbare the differenced ionospheric and tropospheric effects at the two stations related to the two satellites, respectively. For the ionosphere-free combined observables (denoted by j = 4 for distinguishing), the ionospheric error terms have vanished from above equations.

The most important property of the double differences is that the clock error terms in the equation (model) are completely eliminated. It should be emphasised that the clock error, which implicitly affects the computation of the position of the satellite, still has to be carefully considered. Ionospheric and tropospheric effects are reduced greatly through difference forming, especially for those stations that are not very far away from each other. Double differenced Doppler directly describes the geometry change. Double differenced ambiguities can be denoted by

tmpD-145_thumb

The original ambiguities used in Eq. 6.103 are for convenience in case of reference satellite changing.

For the single difference observable vector

tmpD-146_thumb

the double differences

tmpD-147_thumb

can be formed by a linear transformation

tmpD-148_thumb

where E is an identity matrix of size m x m, I is a 1 vector of size m (all elements of the vector are 1), m is the number of formed double differences, and m = n – 1. The cova-riance matrix of the double differences is then

tmpD-149_thumb

For single and double differences

tmpD-150_thumb

the linear transformation matrixtmpD-151_thumband the covariance matrix can be obtained by

tmpD-153_thumb

For the general case of

tmpD-155_thumb

it is obvious that the general transformation matrixtmpD-156_thumband the related covariance matrix can be represented as

tmpD-158_thumb

wheretmpD-159_thumbis an m x m matrix whose elements are all 1, and the weight matrix has the form of

tmpD-161_thumb

where n = m +1. Equation 6.118 can be verified by an identity matrix test (i.e., P • cov(DD(0))= E).

In the case of three stations, supposing n common satellitestmpD-162_thumbare viewed, then the single and double differences can be written as

tmpD-164_thumb

Then one has the transformation and covariance

tmpD-165_thumb

Because of the dependency of the cov(SD) on the baselines forming, cov(DD) is also dependent on the baselines forming. A baseline-wise processing of a network GPS data using double differences is equivalent to an omission of the correlation between the baselines.

Triple Differences

Triple differences are formed between two double differences related to the same stations and satellites at the two adjacent epochs as

tmpD-166_thumb

where t1 and t2 are two adjacent epochs. Supposing the original observables have the same variance oftmpD-167_thumbthen the triple differenced observables have a variance oftmpD-168_thumb Considering Eqs. 6.102-6.104, one has

tmpD-171_thumb

Ionospheric and tropospheric effects are eliminated. If there are no cycle slips during the time, the term of Eq. 6.124 is zero. Therefore, triple differences of Eq. 6.122 can also be used as a check for the cycle slips. Through triple difference forming, the systematic cycle slip turns out to be an effect like an outlier.

The most important property of the triple differences is that only the geometric changing is left in the models. Triple differences of Doppler describe the acceleration of the position.

For double differences

tmpD-172_thumb

one has

tmpD-173_thumb

where

tmpD-174_thumb

Then the related covariance matrix can be represented as

tmpD-175_thumb

where 1tmpD-176_thumbis the double difference transformation matrix of two epochs. Because double differences are independent epoch wise,tmpD-177_thumbis a diagonal matrix oftmpD-178_thumb

tmpD-182_thumb

It is notable that the triple differences formed by epochs (t1, t2) are correlated to the differences formed by epochs (t0, t1) and (t1, t2). Such a correlation makes a sequential processing of the triple difference data very complicated. Sequentially using the above covariance formula indicates an omission of the correlation related to the previous epoch and the next epoch.

Taking the correlation between the baselines into account, an exact correlation description of the triple differences of a GPS network turns out to be very complicated.

Next post:

Previous post: