Guide To Gps Positioning Ebook
Guide To Gps Positioning Ebook
VXAAAAA/
by
Norman Beck Demitris Delikaraoglou Alfred Kleusberg Edward J. Krakiwsky Gerard Lachapelle Richard B. Langley Mete Nakiboglu Klaus-Peter Schwarz James M. Tranquilla Petr Vanicek David Wells
First printing, December 1986 Second printing, with corrections, May 1987
All rights reserved. No part of this publication may be reproduced, stored in a retrieval system, or transmitted, in any form, or by any means, electronic, mechanical, photocopying, recording or otherwise, without the prior permission of the copyright owner.
ISBN: 0-920-114-73-3
Orders: Additional copies may be obtained from Canadian Institute of Surveying and Mapping Box 5378 Postal Station F Ottawa, Ontario Canada K2C 3J1 Telephone: (613)224-9851 Price: $35.00 in Canadian funds (shipping extra) Questions and correspondence: Canadian GPS Associates Box 3184 Postal Station B Fredericton, New Brunswick Canada E3A5G9
Recommended citation and reference: Citation: [Wells et al., 1986] Reference: Wells, D.E., N. Beck, D. Delikaraoglou, A. Kleusberg, EJ. Krakiwsky, G. Lachapelle, R.B. Langley, M. Nakiboglu, K.P. Schwarz, J.M. Tranquilla and P. Vanic'ek (1986). Guide to GPS Positioning. Canadian GPS Associates, Fredericton, N.B., Canada.
Printed in Fredericton, New Brunswick, Canada by the University of New Brunswick Graphic Services. Cover designed by Richard B. Langley.
INTRODUCTION
Goals, content, and format
The Guide to GPS Positioning is a self-contained introduction to the Global Positioning System, designed to be used in any of the following three ways: as a self-study guide, as lecture notes for formal post-secondary education courses, or as hand-out material to support short-course and seminar presentations on GPS. The intended audience includes practicing surveyors, hydrographers, engineers, geophysicists, geologists, geographers, oceanographers, space scientists, and professionals in information management, transportation, forestry, and agriculture. In short, anyone interested in the potential of GPS to provide accurate, inexpensive, consistent positioning. The content has been intentionally kept on the conceptual level as much as possible. The Guide is primarily an introduction, rather than an explication of the details of the latest research results concerning GPS. A minimum of mathematics is used, with each term hopefully fully explained. On the other hand, we attempt to describe all the concepts required to understand and use GPS, and we believe that even those actively using GPS will find this Guide to be of value.
The structure of this guide is modular in both a macro and a micro sense. At the macro level the fifteen chapters in this Guide can be grouped into four units: A description of GPS (chapters 1 to 5) GPS data collection and processing (chapters 6 to 10) GPS applications (chapters 11 to 15) GPS receivers (Appendix A). These four units can be used together to provide a complete overview of GPS, or can be used independently from each other for more specialized audiences. At the micro level, within each chapter, each topic is dealt with in two pages: a graphic or summary page illustrating the concept or listing the main points, and a text page providing a detailed description of the topic. In adopting this format, we followed two earlier GPS publications, Everyman's Guide to Satellite Navigation, by Steven D. Thompson, and The WM GPS Primer, by Rene Scherrer (see Appendix C for details). Our 2-page topic modules are designed to be as independent from each other as possible, so that a specialized presentation on GPS can be assembled by selecting only the appropriate topic modules required.
0.00
TATION
position of receiver R = (X,Y,Z) T (x,y,z) position of satellite r - R| | distance between satellite andreceiver
p =
<$>
di P=c =
carrier phase observation (in cycles) carrier phase observation (in length units) time delay observation pseudo-range observation (in length units) range error due to inaccurate ephemerides satellite clock error receiver clock error range error due to ionospheric refraction range error due to tropospheric refraction integer carrier phase ambiguity
dt dT
trop
A
TABLE OF CONTENTS
0.00
Chapter 1
LOO 1.01 1.02 1.03 1.04 1.05 1.06 1.07 1.08 1.09 1.10
Chapter 2
2.00 2.01 2.02 2.03 2.04 2.05 2.06 2.07 2.08 2.09
Chapter 3
3.00 3.01 3.02 3.03 3.04 3.05 3.06
Chapter 7
7.00 7.01 7.02 7.03 7.04 7.05 7.06 7.07 7.08 7.09 7.10 7.11 7.12 7.13 7.14 7.15 7.16 7.17 7.18 7.19
Chapter 8
8.00 8.01 8.02 8.03 8.04 8.05 8.06 8.07 8.08 8.09 8.10 8.11
Observation equations Pseudo-range observations Carrier beat phase observations Continuous carrier phase Linear combinations of observations Between-epoch (Doppler) single differences Between-receiver single differences Between-satellite single differences Receiver-time double differences Receiver-satellite double differences Receiver-satellite-time triple difference Correlation between observations
Chapter 9
9.00 9.01 9.02 9.03 9.04 9.05 9.06 9.07 9.08 9.09 9.10 9.11 9.12 9.13 9.14 9.15
Chapter 10
10.00 10.01 10.02 10.03 10.04 10.05 10.06 10.07 10.08 10.09 10.10 10.11 10.12 10.13 10.14 10.15 10.16 10.17 10.18 10.19 10.20
Solutions
Solutions Positioning problems Main characteristics of static and kinematic positioning Classification of types of solutions Point positioning using pseudo-range Satellite multi-ranging Least-squares solution for point positioning Relative positioning using pseudo-ranges Relative positioning using phase measurements GPS phase measurement differencing Network solutions Network solution and orbit improvement Weighted short arc solution Absolute kinematic positioning Relative kinematic positioning Combined kinematic point positioning Postmission smoothing Constrained solutions Finding antenna height above the reference ellipsoid Finding the geoid height: accuracy and difficulty Typical geoid shape
I.00 2.00 3.00 4.00 5.00 6.00 7.00 8.00 9.00 10.00 II.00 12.00 13.00 14.00 15.00 16.00 17.00 18.00 19.00 20.00 21.00 22.00 23.00 24.00 25.00 26.00 27.00
Which GPS receiver? CMA-786 (C-set) JMR SatTrak Astrolab II ISTAC-SERIES GPS Positioner Model 2002 ISTAC-SERIES Marine Position Sensor MPS-1 JLR-4000 GPS navigator LTN-700 Macrometer V-1000 Macrometer II T-set MX 1100 GPS Navigator MX 4400 Eagle Mini-Ranger Norstar 1000 PS 8400 Collins Navcore I TR5S GTT-2000 SEL GPS receiver GPS monitor station receiver PAHRS Model 1 TI-4100 NAVSTAR Navigator 4000A GPS Locator, and 4000S GPS Surveyor LORAN-GPS 10X WM101 Other GPS receivers
CHAPTER ONE
Astronomical positioning Basic concept of satellite point positioning Ranges and directions to satellites Ranges to satellites Relative satellite positioning
1.01
POINT
Towards Conventional t z CT International Origin (CIO)
p (point to be positioned)
POSITION VECTOR
R CT =(X CT
CTZCT
Greenwich Meridian
CT
"Equator
This particular coordinate system is the CONVENTIONAL TERRESTRIAL SYSTEM (CT)
Canadian GPS Associates, May 1987
POINT POSITIONING
By positioning we understand the determination of positions of stationary or moving objects [Vanic'ek and Krakiwsky, 1986]. These positions can be determined:
(1) With respect to a well-defined coordinate system usually by three coordinate values. (Here we assume that the 'well-defined' coordinate system is itself positioned and oriented with respect to the earth.)
(2) With respect to another point, taking one point as the origin of a local coordinate system.
The first mode of positioning is known as point positioning, the second as relative positioning. (Some people call relative positioning, differential positioning.) It does not matter, conceptually, which coordinate system is used for point positioning. Nevertheless, it is usual to employ a geocentric system, i.e., a system whose origin coincides with the centre of mass (C) of the earth, as shown here. Therefore, the term absolute positioning is sometimes used instead of point positioning, especially in kinematic applications.
If the object to be positioned is stationary, we speak of static positioning. When the object is moving, we speak of kinematic positioning. Both positioning modes can be and are being used for either kind of positioning. Both modes and both kinds are dealt with in this guide. We note that static positioning is used in surveying; kinematic in navigation. They differ in some important aspects, as we shall see later.
RELATIVE POSITIONING
Relative positioning using terrestrial methods is simpler than point positioning, particularly if the two points are intervisible. Many relative positioning techniques exist, based on a variety of physical and geometrical concepts. Relative positioning is the basic mode used in surveying practice. Any local coordinate system may be used for relative positioning. When survey instruments are used for relative positioning, it is then natural to use the Local Astronomical (LA) system, shown here. As we shall see later for relative positioning from space, other systems are being used. Suppose now that the position R1 of point Pj is known in one coordinate system and the interstation vector AR 12 (AX, AY, AZ) is determined in the same coordinate system. Then the vector equation
=R
gives us the position of P2 in the same coordinate system. If AR12 has been determined in a different coordinate system (*) then it has to be transformed into the same (as Rj) system by the following transformation:
1.03
Provided all
AR^
NETWORK OF POINTS
Often it is advantageous to connect a whole family of points into a network. This mode is normally used for stationary, properly-monumented points. This is because of the necessity to re-occupy the individual points of the network. When these points are positioned together, connected together by means of redundant links, they make a geometrically stronger configuration.
In a network, each link may be considered as an interstation vector and each pair of adjacent points as a pair of points to be positioned relatively with respect to each other. If all the interstation vectors are known derived from observations then the network can serve as a medium for transmitting positions from one end to the other, as shown here. From this point of view, the network mode can be regarded as an incremental mode. As such, it naturally suffers from error accumulation, and care must be taken to curb the systematic error propagation. In surveying practice, networks are the most common mode of positioning. The points linked together by a network are called control points.
For accurate positioning, we have always relied on relative positioning techniques. Relative positioning in the network mode has served as the basic mode of geodetic positioning. Since terrestrially-based optical instruments (transit, theodolite, level, electro-optical distance measuring device, etc.) have been used exclusively to observe the interstation vectors, the links between (adjacent) points have been confined to lines of sight. A horizontal geodetic network is a network for which the horizontal coordinates, latitude and longitude, of the control points in the network are determined as accurately as possible, whereas heights are known only approximately. A geodetic levelling network is comprised of vertical control points (bench marks) whose heights are known as accurately as possible and whose horizontal positions are only very weak.
1.04
AND RELATIVE
/ CONFIDENCE ELLIPSOID
ELLIPSE
R
Relative Error in Relative Position
A R 21
'21
AR 2 i
It is customary among some people to use a proportional error presentation instead of confidence regions. This representation is obtained by dividing the standard (position) error in the desired direction by the distance from the origin of the coordinate system used. Proportional accuracy can be defined for point positioning as well as for relative positioning: one metre error in the geocentric position of a point represents a relative error of 1 m divided by the radius of the earth (6.371 x 106m) and is equal to 0.16 parts per million (0.16 ppm). To get a proportional accuracy of 0.16 ppm in a relative position of two points 10 km apart, their relative position would have to be known with o = 1.6 mm.
The proportional errors of relative positions of adjacent points in the present horizontal networks in Canada are at best about 10 ppm (first-order control). For levelling networks, the relative accuracy is about one order of magnitude higher. Both horizontal and levelling networks are now being redefined, re-adjusted, and generally improved to become more compatible with new positioning capabilities. In North America, the new horizontal network contains many Transit (Doppler) satellite-determined positions that have been included to strengthen the geometry, location, and orientation of the network.
Another concept which has been frequently used to describe accuracy is the concept of Dilution of Precision (DOP). The DOP is, in essence, a root-sum-squared measure of the size of the confidence region.
Both the DOP as well as confidence regions stem from a more general concept of a covariance matrix, which in this context is an assembly of variances (squares of standard deviations) and covariances of individual coordinates of a set of points of simultaneous interest. These coordinates may be one or several dimensional. For example, a set consisting of one point in three dimensions has a 3 by 3 covariance matrix associated with it. Different kinds of DOPS are nothing but the square roots of the traces of different submatrices of the covariance matrix, divided by the measurement standard deviation, whereas the semi-axes of the confidence regions are related to the eigenvalues of this matrix [Vanicek and Krakiwsky, 1986].
1.05
.TUM;
GEOCENTRIC REFERENCE ELLIPSOID
If the geoid is replaced with a biaxial ellipsoid, one can define geometrical heights h, also called heights above the (reference) ellipsoid, which have very little use in practice. However, their attraction is that they can be obtained directly from three-dimensional Cartesian coordinates of the point of interest, provided the location and orientation of the reference ellipsoid in the Cartesian coordinate system is known [Vanicek and Krakiwsky, 1986]. The equation that links the two kinds of heights is shown in the figure, where the elevation N of the geoid above the reference ellipsoid is usually called the geoidal height. Globally, the absolute value of N with respect to the best fitting geocentric reference ellipsoid is almost everywhere smaller than 100 metres.
Whereas the use of the ellipsoid as a reference for heights is impractical for many tasks, its use as a reference for horizontal coordinates, latitude cp and longitude A,, is widespread. This is why such reference ellipsoids are commonly called horizontal datums. Many horizontal datums (both geocentric and non-geocentric) are used throughout the world. Clearly, if we were happy describing positions by Cartesian coordinates in a system such as, for example, the Conventional Terrestrial system, the use of both vertical and horizontal datums could be abandoned. This is, conceptually, the case with positioning by satellites. However, in practice, satellite-determined Cartesian | coordinates must usually be transformed to (p, K, and H. This transformation, as stated above, requires knowledge of: the location and orientation with respect to the Cartesian coordinate system of the | horizontal datum selected, the height of the geoid above the horizontal datum.
1.06
Local Zenith
Greenwich Meridian
Local Meridian
Astronomical! Longitude
Astronomical Latitude
This figure does not illustrate a coordinate system in the strict sense, only a concept
Canadian GPS Associates, May 1987
ASTRONOMICAL POSITIONING
To accurately measure point positions on the surface of the earth has been a human endeavour since time immemorial. The only technique available before the advent of satellites was optical astronomy. The basic idea of astronomical positioning is to find the position of local zenith among stars, as shown here. From the relative position of local zenith with respect to the surrounding stars (whose positions, in turn, are already known), the astronomical latitude and astronomical longitude are derived. The equations for doing it can be found in Vanicek and Krakiwsky [1986].
The main problem with astronomically-determined 'positions' is that they have to be transformed into a well-defined coordinate system, to be useful for accurate positioning. This task requires a good knowledge of the earth's gravity field. Without this knowledge, the distortions in astronomically derived 'positions' may reach many hundreds of metres. Nevertheless, such 'positions' may still be quite useful for rough positioning or navigation on land, in the air, or at sea. With good knowledge of the gravity field, the errors in positions after the transformation are of the order of 10 metres, i.e., the proportional accuracy of those positions is of the order of 1 ppm.
Today, astronomical positioning is not used much for precise positioning. Rather, the astronomical 'positions' are used to study the earth's gravity field; they provide an easy tool to evaluate the deflection of the vertical, the slope of the geoid with respect to the selected reference ellipsoid.
1.07
jth Satellite
MEASURED
1 in Antenna :-
tti
SOUGHT
Ri = r j - e j p j i
R. is the position vector of the i th antenna
r * is the position vector of the j th satellite * 6. p. is the range vector between the two.
Canadian GPS Associates, May 1987
To predict accurately the position of the satellite rj(t) in time is a rather difficult proposition. The task of predicting the ephemeris which is the proper name of rj as a function of time calls for a specialized knowledge of satellite dynamics which historically belongs in the field of celestial mechanics. The ephemerides are usually determined and predicted (in time) by the operators of the satellite system. It is, however, possible for the user of the system to improve the given ephemerides. It is also possible for the users to establish an independent network of tracking stations, preferably a global one, and compute the ephemerides in a postmission mode. It should be clear that if the antenna changes its position Rj in time (kinematic positioning), then R i has to be continuously re-evaluated. If the antenna is stationary, then the re-evaluation should always lead to the same result, and we will have a redundancy of results. In many applications the position of the antenna itself is not of primary interest. It is the position of the survey marker, the centre of a hydrographic launch, etc., that is to be determined. For this reason, the offset (eccentricity) of the antenna with respect to the derived point must be established or monitored.
Photograph
RA
Right Ascension This particular coordinate system is the RIGHT ASCENSION SYSTEM (RA)
Canadian GPS Associates, May 1987
It is worth noting, however, that optical direction measurements remain an important means of determining the positions of satellites (satellite tracking). Catalogues of the thousands of objects in orbit around the earth are maintained by NASA, and other agencies, based in part on photographs of these objects from known ground stations, against the star background.
1.09
r1, r 2 , r
rJ-R,||= P
Canadian GPS Associates, May 1987
j =1,2,3
RANGING TO SATELLITES
Much better results have been obtained from observing ranges. When all three components of the vector ep are not observed simultaneously, a minimum of three non-coplanar ranges have to be observed to get the required position, as shown here. If Rj is the position vector of a stationary antenna, then the three (or more) noncoplanar ranges can be measured at different times and hence may (or may not) pertain to one and the same satellite. If the antenna is moving, then to get a position at any instant of time, often called the fix, the three (or more) ranges have to be measured simultaneously. This infers that, for ranging in a kinematic point positioning mode, at least three satellites have to be 'visible' to the antenna at the same time. Ranging can be done using different parts of the electromagnetic spectrum. The visible part of the spectrum is used in laser ranging; NAVSTAR/GPS uses L-band frequencies (1.2 to 1.6 GHz), whereas the Transit (Doppler) system uses VHP and UHF radio frequencies (150 and 40| MHz). Generally, laser ranging is more accurate than radio ranging, but it is less practical because of its limitations due to visibility, portability, and excessive cost of equipment. Laser ranging has an accuracy of about 5 cm (one sigma level) from three month's worth of ranging [Smith et al., 1985]. The corresponding accuracy from GPS cannot be assessed because no such experiment has yet been conducted, to the best of our knowledge. It is thought, however, that an accuracy of a few decimetres could be eventually achieved.
RELATIVE SATELLITE
1.10
Conceptually:
= r J - ej '- ej2pj2
AR12=R2-R1=
& Canadian GPS Associates, May 1987
The effects limiting the accuracy of radio positioning are: the atmospheric propagation delay, drift in the satellite and possibly the receiver clocks, and error in satellite ephemerides. Some of these effects can be suppressed, or even eliminated, if relative positions, rather than point positions, are sought. Put in other words: it is possible to determine relative positions, and thus also geometrical configurations of points, much more accurately than it is to locate these points on the earth's surface.
This is true for both ranging and range differencing. It is even more the case when simultaneous, rather than sequential, observations to each satellite are used. The relative positioning mode was used extensively with the Transit system under the name of translocation. The best proportional accuracies currently obtainable in relative positioning are 5 to 10 ppm, when using the Transit system, and 0.1 to 2.0 ppm, when using GPS. The relative positioning ranging mode can be used for both static and kinematic kinds of positioning.
CHAPTER TWO
2.00
IITIONE
LAND APPLICATIONS
MILITARY APPLICATIONS
2.01
TIONS IN
Geodynamics
Deformation Monitoring
Control Networks
Canadian GPS Associates, April 1986
cadastral surveying geodetic control local deformation monitoring global deformation monitoring.
Surveying for cadastral purposes normally requires an accuracy in relative position of 10~4. This accuracy can be easily achieved utilizing GPS observations. However, for a wide application of GPS in this field, completely automated processing, if possible on site, is required. It is expected that this feature will be available in GPS receivers within the next few years.
Geodetic control can involve high accuracy networks, from which lower order networks are established. The required accuracies in relative positioning are of the order of 5 x 10'6 to 1 x 10"6 for distances of 20 to 100 km. This accuracy level is routinely achieved by postprocessing GPS carrier phase observations with standard software. Lower order geodetic control (for mapping, for example) can also be established economically by GPS.
Local deformation monitoring (mining subsidence, structural deformation) requires accuracies of 1 mm to 1 cm over distances of up to a few kilometres. The presently achievable accuracy for these applications is limited by uncertainties in the variation of the electrical centre of the GPS antennas and signal distortions due to the reflective environment of the antennas. Further problems arise from limited satellite visibility due to signal shadowing in a typical industrial environment.
Global deformation monitoring (plate tectonics, etc.) requires accuracies of 10"7 to 10"8 over intercontinental distances. The main differences compared with the previously discussed applications lie in the necessity for an elaborate modelling of the GPS satellite orbits, atmospheric signal propagation delays, and other biases.
LAND APPLICATIONS IN
2.02
NAVIGATION
MONITORING
Monitoring the location and movement of vehicles can be achieved if they are equipped with automated transmitters in addition to the GPS receiver. The position determined by the GPS receiver/processor can be transmitted to a central site for display. This application can be of significant importance for controlling and dispatching boxcars, taxi fleets, truck fleets, etc.
2.03
Shoal surveying
Canadian GPS Associates, April 1986
The positioning requirements for field work in physical oceanography can be met by GPS. Carrier phase measurements provide additionally precise ship's velocity, which is needed in ocean current studies.
E APPL
2.05
| Bathymetry
In aerial photogrammetry, sensor (camera) navigation accuracy requirements of several tens of metres can be met easily by GPS. High precision postmission sensor positioning (and sensor orientation) by GPS can replace aerotriangulation and thus render ground control points superfluous. Positioning accuracy requirements for this application vary between 0.5 m and 25 m for different mapping scales.
Airborne laser profiling can be used for direct digital mapping of the topography (digital terrain model) if the sensor (laser) position is known with an accuracy of 0.5 to 1 m vertically and several metres horizontally. GPS is expected to provide this position accuracy in postmission analysis.
Roughly the same positioning requirements have to be met for airborne gravity and gravity gradiometry. In these applications, GPS measurements can additionally determine the sensor velocity needed for the Eotvos reduction of the gravity data.
Airborne laser bathymetry and radar imaging require less accurate sensor positions which can be easily determined by GPS measurements.
2.06
Forestry spraying
En route navigation
Nonprecision approach
In other airborne applications (forestry and crop spraying, etc.), where navigation is not required for aviation safety but for the deployment of freight, GPS can easily meet the navigation accuracy requirements.
SPACE APPLICATIONS
Space applications of GPS consist mainly of positioning and navigation of other space vehicles carrying geophysical and/or geodetic sensors. Usually, these satellites have low orbits and, therefore, the measurement geometry is similar to terrestrial applications. Typical examples are satellite remote sensing and radar altimetry. Satellite positions obtained from GPS measurements can be used to improve and/or simplify the orbit computations for these space vehicles or even replace the orbit integration by discrete orbit positions.
2.08
Mountain climbing
RECREATIONAL APPLICATIONS
It is expected that the cost of GPS receivers will continue to decrease. At a price level of several hundred dollars, a GPS receiver will be affordable for non-professional users. At the same time, advances in microelectronics will have reduced the volume, weight, and complexity of receivers. In this case, recreational activities will provide a huge market for these low cost, portable, and easy-to-use receivers. Low to medium accuracy positions (tens to hundreds of metres) will be used by yacht and boat owners, hikers and bikers, safari participants, and many others.
/mip / ^ ^^
^|0S /i )
MILITARY
*
2.09
POSITIONING AND NAVIGATION en route navigation low-level navigation target acquisition photo reconnaissance
remotely operated vehicles
i i
updating inertial navigation systems sensor emplacement missile guidance command and control
NUCLEAR DETONATIONS
Canadian GPS Associates, April 1986
MILITARY APPLICATONS
The Global Positioning System is primarily designed for real-time military positioning. Military applications encompass airborne, marine, and land navigation. GPS will serve as a stand-alone system and as a part of integrated navigation systems. This figure shows a list of the main military applications. In addition, the GPS satellites carry sensors to detect and monitor nuclear detonations.
CHAPTER THREE
XT
3.00
IN
TEI
5TRIAL
THE PAST Optical Systems Early Radio Ranging and Direction Finding Systems THE PRESENT Transit Satellite Positioning Principle of Doppler Positioning Transit Doppler Measurements ARGOS Global Positioning System Satellite Laser Ranging
Very Long Baseline Interferometry New Technology Accuracy vs. Distance GPS and Other Systems
THE FUTURE
ED IN
3.01
Positional Astronomy
cflS
Ofi 00
fe
^2 O o HH ^
oo
oooo oo
H eu
i/> oe
01)
S So,
W
C3 U C3 V *^ ; -*;
s. ^ - 1 0 c5 **
**
,2 s S o
s w o a
a o w s
H cc
M
C*
^ _rt ^ PC 'S PC T2 g 2 ^! *- 5* ** S3
iH O
OPTICAL SYSTEMS
Optical systems are those systems which utilize the visible part of the electromagnetic spectrum. This definition would encompass astronomical positioning using a theodolite or sextant and also laser ranging. However, we will not consider these techniques in this section. Conventional astronomical positioning has been discussed in Chapter 1, and laser ranging will be. discussed later in this chapter.
Optical systems range in accuracy and complexity from simple binoculars to 60 cm tracking cameras. They have been used in the past for both satellite tracking (i.e., orbit determination) and terrestrial positioning but only the tracking cameras provided geodetically-useful results [Henriksen, 1977]. A tracking camera records a photographic image or a series of images of a satellite, either from a flashing beacon or via reflected sunlight, together with the images of background stars. The time of exposure is also precisely recorded. The positions of the satellite images on the processed film are subsequently measured with respect to stars in the field of view. If the positions of the stars are known, then the topocentric right ascension and declination of the satellite can be determined. Three such determinations combined with the known coordinates of the camera are sufficient for determining the orbit of the satellite. Additional observations can be used to improve the orbit's accuracy.
Alternatively, if the orbit of the satellite is known, it is possible to determine the geocentric coordinates of the camera. Two techniques have been used. In the geometric method, two or more cameras simultaneously record satellite images and the camera positions are then determined by triangulation. In the dynamic method, observations are not necessarily synchronized but precise knowledge of the satellite orbit is required. Several three-dimensional networks were established using these techniques [Henriksen, 1977]. One of them was the BC-4 worldwide satellite triangulation network, a network of 45 stations whose coordinates were established through observations of the balloon satellite PAGEOS, using Wild BC-4 cameras. The relative positions of the stations were determined with an estimated accuracy of about 5 m. Greater accuracy in position determination using satellite photography is not readily achievable due to a number of error sources. Chief among these are scintillations of the satellite due to atmospheric turbulence and distortions in the photographic emulsion. Consequently, in geodesy, satellite photography has been superseded by other techniques.
Optical systems are still an important source of information for satellite orbit determination and surveillance. They figure prominently, for example, in the U.S. Air Force's Spacetrack System and its Ground-based Electro-Optical Deep Space Surveillance (GEODSS) System.
M
PQ
OS
C/5
PQ
C3 ** 2 ^ S * a* u
00
to
so
oo
O 10
t> o
(=J
u
a -
ss ^
8 a
W
The GRARR system relies on a transponder in the satellite to respond to an S-band signal transmitted from an antenna on the ground. Superimposed on the carriers by phase modulation are 6 to 8 frequencies between 8 Hz and 500 kHz. The relative phases of the transmitted and received modulations are used to infer the range to the satellite. The multiple modulation frequencies are used to remove the inherent cycle ambiguity. The phase rate or Doppler shift of the received carrier is also measured. The range accuracy of GRARR is about 5 m.
Minitrack, or more fully Prime Minitrack, was NASA's first tracking system. It consists of a pair of crossed interferometers which measure the relative phase of a nominally 136 (formerly 108) MHz signal transmitted by suitably-equipped satellites as received at the two ends of the interferometer. The same principle is used in the U.S. Naval Space Surveillance System and in very long baseline interferometry (see later in this chapter). With Minitrack, the altitude and azimuth of a satellite can be determined to about 20".
The SECOR system, which was operated between 1962 and 1970, consisted of four ground stations each transmitting modulated signals on 420.9 MHz to a transponder on certain satellites. The transponders returned the signal on 449 and 224.5 MHz. By measuring the relative phase of the transmitted and received 585.533 kHz modulation, the range to the satellite could be determined. An additional three modulation frequencies provided ambiguity resolution. Range accuracies of 3 to 10 m were obtained.
3.04
Receives orbit data Measures integrated Doppler shifts Computes his position
3.05
The classical explanation of the effect is that the observer receives more wave crests per second, i.e., the frequency is increased if the source and the observer are moving closer together, whereas fewer wave crests per second are received, i.e., the frequency is decreased, if the source and the observer are moving farther apart. If the relative speed of the source and observer is much less than the speed of light, then the received frequency is given approximately as f r f s ( l - 1/cdS/dt) where f s is the frequency at the source, c is the speed of light, and S the distance or range between the source and the observer; dS/dt is the range rate.
Returning to the train at the level crossing, you may have noticed that the closer you are to the track, the faster the change in pitch of the horn. And even if you could not see or feel the train, you can tell when it passes the crossing (the point of closest approach) by noting the instant when the pitch of the horn is mid-way between the high and low extremes (fs). Therefore by monitoring the frequency of the received sound as the train passes and knowing its assumed constant speed, you can establish your position in a two-dimensional coordinate system where the x-axis, say, runs along the track and the y-axis runs perpendicular to it. The origin may be assigned arbitrarily. This is the principle of Doppler positioning.
In the case of a Transit satellite (or any other satellite for that matter), the position of a receiver can be established by continuously recording the Doppler shift of the received signals (or the number of cycles of the Doppler frequency which is a more precisely obtained observable). Subsequently these data are combined with accurate coordinates of the satellite to determine the position of the receiver. As with the passage of a train, a single satellite pass can provide at most only two coordinates of the receiver's position. Whereas this may be satisfactory for navigation at sea where the height above the reference ellipsoid is approximately known, three-dimensional positioning requires observing multiple satellite passes.
1IT 5UREMENT,
3.06
~ 80 ppm
25 ppm
10
12
14
3.07
NOAAj
Platform i (1=1,4000)
ODlOD
ARGOS
Another satellite system which uses the Doppler effect for positioning is Argos [Service Argos, 1978], a cooperative project of the French Centre National d'Etudes Spatiales (CNES), NASA, and the U.S. National Oceanic and Atmospheric Administration (NOA A). However, in contrast to the Transit system, the transmitters are operated by the users and the receivers are in the satellites. An Argos transmitter on an instrumented 'platform' of some sort (oceanographic or navigation buoy, radiosonde balloon, remote weather station, etc.), periodically emits a 401.65 MHz signal carrying information from the platform's sensors. One of two passing U.S. TIROS/N-class weather satellites picks up this signal and records its Doppler shift along with the sensor data. These data are subsequently played back when the satellite is in range of one of three tracking stations: Wallops Island, Virginia; Gilmore Creek, Alaska; or Lannion, France. The tracking stations relay their data to NOAA's National Environmental Satellite Service in Suitland, Maryland, where they are sorted and then passed on to the Argos Data Processing Centre at CNES in Toulouse, France. CNES computes the position and velocity of the platform from the recorded Doppler shifts. This information along with th sensor data is conveyed to the operator of the platform by telex or letter. The two-dimensional (latitude and longitude) positions and velocities are claimed to be accurate to 3 km and 0.5 m s"1 (3a), respectively. Up to 4000 platforms requiring location service can be handled by the Argos system, assuming these are uniformly distributed over the earth's surface.
The Argos receiving system was first implemented on the prototype TIROS/N spacecraft, orbitted in 1978. Subsequent TIROS/N satellies, NOAA-6 through NOAA-10, have also carried the Argos Data Collection System.
Position determination from signals uplinked to a satellite is also utilized in the COSPAS-SARSAT search and rescue system [McPherson, 1981].
3.08
CONTROL SYSTEM Time synchronization Orbit prediction Data injection Satellite health monitoring
Canadian GPS Associates, May 1987
3.09
LAS:
p= cCT2-T1)/2|
\
1
r
>r a*"
PHOTO1 OTVYO
1 CL
J
TAPE RECORDER
The first SLR systems began operating in the early 1960s. At the present time there are over 25 fixed, precision SLR systems worldwide. In addition, NASA has developed and deployed 8 MOBLAS-class mobile systems and 4 TLRS easily transportable systems; a further two transportable systems, MTLRS-1 and -2, have been built in Europe.
Most SLR systems now use short pulse length (100 ps) neodymium YAG lasers with an output at 532 nm and a power of a few ml per pulse. Some systems still use the more powerful (several joules per pulse) but wider pulse width (several nanoseconds) ruby lasers with output at 694 nm. Most SLR systems have separate telescopes for transmitting and receiving; some use a single telescope. The receiving telescopes used range in size from 28 cm diameter to a metre or more. The receiving sensitivities of systems range from a single photon to more than 100 photons. The best single shot range accuracy is about 2 cm r.m.s.
SLR data are usually analysed using the dynamical approach in which data from many stations spanning periods from about 5 to 30 days are combined. Station positions, orbital elements, and various biases are determined simultaneously. The lengths of global intersite baselines have been determined to precisions of 3 to 5 cm. The SLR networks also provide accurate values of the earth rotation parameters.
Over 15 retroreflector-equipped satellites have been placed in earth orbit. Of these, the most fruitful have been the dedicated LAGEOS and Starlette satellites. Recently, Japan launched the third fully reflecting satellite, the Experimental Geodetic Payload, 'Ajisi.' Laser ranging to the earth's natural satellite, the moon, is also carried out by a few specially equipped stations. Reflectors were deployed on the lunar surface by the crews of Apollo 11, 14, and 15 and remotely by two Soviet Luna landers.
Suggested further reading: Degnan, J.J. (1985). "Satellite laser ranging: Current status and future prospects." IEEE Transactions on Geoscience and Remote Sensing, Special Issue on satellite geodynamics, Vol. GE-23, No. 4, pp. 398-413.
3.10
3RY L(
INT:
ASELINE
,- x
QUASAR|
B-s = c T I
xv>
sv
H-MASER CLOCK
H-MASER CLOCK
The primary observable of geodetic VLB I is the group delay, the difference in arrival times of the quasar signal wavefronts at the radio telescopes. In principle, the delay can be measured in the correlation process by noting the time offset between a pair of tape recordings required to achieve maximum correlation. The phase (delay) of the correlation function and its time rate of change, the delay rate, are also measured. In practice, the group delay is obtained from measurements of the phase delay at different frequencies.
The primary component of the measured group delays is the geometric delay, T = (B-s) / c, where B is the baseline vector connecting two radio telescopes, s is the unit vector in the direction of the radio source, and c is the speed of light. From observations of a dozen or more radio sources during a nominal 24-hour session, the three components of the baseline vector can be retrieved. A number of biases in the data must be carefully measured or modelled.
The radio frequencies used for geodetic VLB I are in the microwave portion of the electromagnetic spectrum. The Mark-Ill VLBI system [Clark et al., 1985] uses 14 frequencies at S band (2.2 - 2.3 GHz) and X band (8.2 - 8.6 GHz); the two bands being used in order to correct for ionospheric delay (see Chapter 9).
t=1
.J
a I a I
_a
3.12
ONE - WAY SYSTEMS (military) GPS = NAVSTAR Proposed in the early 1970s 11 prototype satellites (1978-86) 28 production satellites (1989-95) TRANSIT = NNSS Operational (1964-94) Typically 4-6 usable satellites Impact on control surveys
(civilian)
GEOSTAR (U.S. COMMERCIAL) Electronic packages piggy-backed on communication satellites Low cost (U.S.$450 + subscription) Planned U.S. coverage by late 1980s 100 000 users by 1990
NAVSAT (ESA PROPOSAL) Two-way GPS Regional control centres Deployment not yet funded
Canadian GPS Associates, May 1987
Civilian systems need not impose radio silence on the user or require the satellites to operate autonomously during a crisis. Consequently two-way systems are possible. The complexity of the system can then be transferred from the satellites and the receivers to the control segment. Two of the proposed two-way systems are Geostar and NAVSAT. Geostar requires two-way communication between the user and the satellites. NAVSAT requires continuous two-way communication between the satellites and the control stations. These systems are described on the following pages of this chapter.
3.13
GEOSYNCHRONOUS COMMUNICATION
CENTRAL COMPUTER
GEOSTAR
Geostar is one of several satellite-based positioning systems [Anon., 1986] that have been proposed by commercial interests. The concept of the system has evolved since 1983 when it was first proposed. The current design [O'Neill, 1984] calls for two geosynchronous satellites to relay two-way signals between a central computational facility, 'Geostar Central,' and user terminals. Spread-spectrum signals are transmitted to one satellite which retransmits the signals to the user terminals. The user terminals transmit identically-coded signals back through both satellites to the central facility using a time-division multiple-access scheme, adding a unique digital identification code for security and billing purposes. The two-way travel times of the signals via each of the two satellites are then determined for each of the user terminals. By combining the two ranges with height information obtained from a digital terrain map, the two-dimensional coordinates of a land- or sea-based user can be determined as often as required. Accuracy is improved by combining data from a fixed reference station within 200 km of a user terminal. The system also has a message capability. Messages can be sent from the central facility to users or vice versa. For aircraft in flight, height information can be supplied by the pilot or through automatic interrogation of an encoding altimeter. The two satellites would cover the continental U.S. and peripheral regions. Four additional satellites could provide essentially worldwide coverage (except at high latitudes). The system is estimated to have sufficient capacity to serve a user population of 50 million user terminals. As all of the positioning computations are performed at the central facility, the user terminals can be kept very simple. A terminal consists of a transceiver to relay the radio signals and processing and message storage integrated circuits. Messages and position information are displayed on a liquid crystal panel and outgoing messages entered on a calculator-like keypad. Hand-held versions of the terminal will be powered by penlight batteries, weigh about half a kilogram, and could cost as little as $US 450 (1984 dollars) if produced in quantities of the order of 100 000.
Initial tests of the Geostar concept were carried out in the fall and winter of 1983/84 using 'pseudo-satellites' on mountain tops in the vicinity of Lake Tahoe, Nevada. Those tests, using pedestrians, automobiles, and aircraft, demonstrated a positioning accuracy of about one metre. The first true satellite test (of message capability from user terminals to Geostar Central) began in March 1986 with the launch of GTE Spacenet's GStar 2 communication satellite. The satellite carries a one-way Geostar transponder called LinkOne. Unfortunately, both the primary and backup LinkOne L-band receivers failed in May 1986 [Anon., 1986].
3.14
NAVSAT
-12 hr orbit W inclined 63 to equatorial plane
39 500 km
^""lllimilllllluilllllluilNIIIJIIIIIIIIIUIINIIIIIIMIIIIIII'l'1
24 hr geostationary orbit; 36 500 km altitude ORBITS TYPE co(o) # SATS (*) elliptic 63.4 0 270 2 elliptic 63.4 120 270 2 elliptic 63.4 240 270 2 circular 0 0 0 6
Diagram patterned after one by ESA Canadian GPS Associates, May 1987
PHASES() 180,240 168,205.5 154.5, 192 0, 60, 120, 180, 240, 300
3.15
The near-term trend in satellite positioning will be toward smaller, cheaper, and more easily-operated receivers. A hand-held GPS receiver not much larger than a standard deck of playing cards is under development by Rockwell Collins for the U.S. Department of Defense. The single-channel receiver with a pop-up antenna will make use of the latest gallium arsenide very high speed, very large scale integrated circuit technology. We foresee a natural evolution toward even smaller, more accurate receivers. Eventually, we may even have the 'wrist locator,' a cheap ($10 in current funds), accurate (1 mm) device that would be as ubiquitous as today's electronic wrist watches [Vanicek et al., 1983]. Such devices will influence our future daily lives in ways that are presently impossible to predict.
CHAPTER FOUR
4.00
SATELLITES GPS satellite Prototype satellite constellation Final 21-satellite constellation Satellite ground tracks GPS ground stations GPS SIGNAL STRUCTURE MEASUREMENTS Propagation and ranging concepts Pseudo-range and carrier beat phase measurements RECEIVER COMPONENTS AND FUNCTIONS POSITIONS Static, kinematic, and relative positioning Geometry and accuracy Dilution of precision (DOP) Constrained solutions IMPACT U.S. Federal Radionavigation Plan How many GPS users will there be ? Main features of GPS Present and future GPS characteristics
Canadian GPS Associates, April 1986
Signals. The next figure quickly reviews the concepts behind the rather complicated GPS signal structure.
Measurements. The next five figures review propagation and range-measuring concepts, and introduce the two types of GPS measurements pseudo-ranges and carrier beat phase measurements. Receivers. The next two figures describe the components and functions of a typical GPS receiver.
Positions. The next four figures present the concepts involved in static and kinematic positioning; relative positioning; multi-receiver positioning; and relative kinematic positioning.
Accuracy. The next three figures introduce the relationship between the GPS satellite constellation geometry and positioning accuracy; the Dilution of Precision (DOP) concept; and the effect of constrained solutions.
Impact. So far the material in this chapter has previewed the remainder of the Guide. In this final set of five figures, we discuss the policies and effects of the United States Government Federal Radionavigation Plan; the possible size of the GPS user community; the features of GPS that will form the basis for its impact on the positioning professions and on society at large; and a speculative look at the future of high-technology positioning capabilities.
4.01
GPS SATELLITE
So far three generations of GPS satellites have been planned. These are known as the 'Block I,' 'Block II,' and 'Block IIP satellites.
Eleven Block I satellites were constructed and launched between 1978 and 1985. These are the prototype satellites presently in use.
Construction of the 28 Block II satellites is nearing completion. Launch of 21 of these will signal the beginning of the GPS operational era.
Some Block II GPS satellites will weigh 845 kg (upon insertion into final orbit). Other Block II GPS satellites will be heavier, having a NUDET (NUclear DETection) module on board. Block II satellites are designed to fulfill their function for 7.5 years. This figure shows what these satellites look like in orbit.
(a) Receive and store information transmitted by the 'control segment,' i.e., by the operators of the system. (b) Do limited data processing on board by means of its own microprocessor. (c) Maintain very accurate time by means of several oscillators (two cesium, two rubidium) carried on board. (d) Transmit information to the user by means of various signals. (e) Manoeuvre by means of thrusters controlled by the system operators.
The activities of the satellites are powered by batteries charged by solar panels which when deployed have an area of 7.25 m2.
The Block III GPS satellites are presently at the design stage. It is planned that they will be available for launch, to replace the Block II generation, beginning in the late 1990s.
(N O
.9.9
ta ts fc s
'cd 3
-^ c^
00
CO
* $ %
CM
CSJ^^Or^OCS^H^-I^HOO
8
H
.2
(NoocNcS'4c^r^voa\o
oooooooooo*-HcOTf'<frin t^-r^t^r-oooooooooooooo
a (3
id
o o o ^H o o
OO OO OO OO O O
t^ r^ r^ t^ oo oo ON ON ON ON ON ON
oooo CO Tt Tj- IO OO OO OO OO
ON ON ON ON
co 13 g < tj I
oo S <. o jo
ON ON T-H ON OO CO t^ CS
T-H O CS T-H
^H CS CO rf
1H
T-H
OO
TI
T-H
a
Er o DH O
1I
1(
1I
T(
T-H
T-H
ON CN oo rr^ in t^ en
ON ON ON VO
T-H CO C<l
VO
OO
ON
uu<u
<D
43 G S
o <o x>
4.03
11
10 N A n 9 v
p
R N
12 13 11 9 8 6
0 6 12 Fredericton, April 28, 1986
8 6
4
S T A
18 24 Hour (UT)
Presently (April 28, 1986) in Fredericton, New Brunswick, this five-hour period falls between 3 am and 9 am Universal Time. We see from the figure that the sequence of satellite risings is as follows:
orbit position C-5 rises (NAVSTAR 11, PRN 3 ) orbit position C-3 rises ( NAVSTAR 4, PRN 8 ) orbit position A-3 rises ( NAVSTAR 3, PRN 6 ) orbit position C-2 rises (NAVSTAR 8, PRN 11 ) orbit position A-2 rises ( NAVSTAR 6, PRN 9 ) orbit position A-l rises ( NAVSTAR 10, PRN 12 ) orbit position C-l rises ( NAVSTAR 9, PRN 13 ).
Due to the orbit configuration, orbit position C-5 sets below the horizon about ten minutes before positions A-l and C-l rise. Thus the maximum number of visible satellites is six out of the seven now available. We see that this occurs during the period from the rise of A-l and C-l, and the set of A-3.
This situation differs from one part of the world to another.
Because the GPS satellites circle the earth once every 12 hours, configurations repeat themselves every twelve hours (successively in opposite hemispheres). The orbital periods are 12 hours of sidereal time, and due to the difference between sidereal and solar time, the satellites appear about four minutes earlier each day. Thus each month the coverage period will advance by two hours, or roughly 30 minutes per week.
q ^f*
Due to the delays resulting from the space shuttle Challenger disaster, shuttle flights will not begin again before early 1988. The shuttle program was designed to replace supposedly more expensive rocket launching of satellites into orbit, not only for NASA but for the U.S. military as well. With the problems and delays in the shuttle program, a return to rocket launches appears likely. The USAF has contracted four design studies for conventional rocket launching GPS (and other) satellites. Rocket launches may start in late 1988. Three scenarios are shown here. If shuttle-only launches are used, it may take as long as four years to place the full GPS constellation of 18 satellites in orbit. If a combination of shuttle and rocket launches is used, the most optimistic date for achievement of the full constellation is late 1989. The most likely guess, using both shuttle and rocket launches, is early 1991.
For a continuous three-dimensional navigation capability, at least 4 satellites have to be everywhere visible simultaneously at all times. Four satellites enable the user to eliminate unknown time offset between satellite and user oscillators. Similarly, for two-dimensional navigation capability, 3 simultaneously visible satellites are necessary. Dates when these capabilities may be reached are shown here.
GPS can be used in the navigation mode, three dimensional or two dimensional, even before these dates with the understanding that the coverage will not be continuous in time. This temporal discontinuity affects stationary applications less dramatically. It results only in the necessity of spending more time 'observing' the satellites thus making the position determinations more expensive.
'STELLATION IN
4.05
18 satellites, plus 3 active spares 6 orbit planes 12 hour period 20,000 km height (almost) full coverage (24 hours per day everywhere in the world)
The satellites are in circular orbits approximately 20 000 kilometres above the earth's surface. This is more than three times the earth's radius. At this height, the satellite revolves around the earth once every twelve hours. For an observer on the surface of the earth, any particular satellite will be above the horizon for about five of these twelve hours. At various times of day, and at various locations, the number of satellites above the horizon may vary between a minimum of four and a maximum of ten.
The coverage provided by this '18+3' constellation is not perfect at certain localities, and at certain times of day, when only four satellites are available, the satellite geometry provides much worse positioning accuracy than usual. Such occurrences are brief (lasting less than 10 minutes usually). They are referred to as outages.
Constellations proposed earlier included: 24 satellites in three orbital planes of eight satellites each. This would have provided fewer outages than the current proposal, but was shelved when budget cuts reduced the number of satellites to 18; 18 satellites in three orbital planes of 6 satellites each. This resulted in more serious outages than the one currently adopted.
vi +4
a o
V)
H t=4 ^ .J H H
-2 "3
-4^
C O
ca 00
U .J
* OO
O\
4.07
03:00 to 05:00 UT
P numbers
Block II Satellites
Canadian GPS Associates, April 1986
Note that the 4 minutes per day regression of the satellites appears as an eastward shift of the tracks by approximately 1.
These plots were prepared using program MacSat, available from the Geodetic Research Laboratory, Department of Surveying Engineering, University of New Brunswick.
4.08
iATELLITE
rises at 4:18
sets at 4:07
rises at 4:25
rises at 3:25
270
sets at 4:29
These plots correspond to the track plots of the previous figure. Both plots are for Fredericton, April 28, 1986, 03:00 to 05:00 UT. The top plot is for the present constellation. The bottom plot assumes that the final 18-satellite constellation was available for the same day and time.
These plots represent a handy tool to judge the satellite configuration at any given instant. For instance, the top plot shows that between the setting of PRN 3 at 4:07, and the rise of PRN 12 at 4:18 the configuration is poor; the remaining four satellites are bunched together (just past the middle of their two-hour tracks) towards the northwest, and all are at a relatively high elevation. However, once PRNs 12 and 13 have risen but remain at low elevations, and the other satellites have progressed further along their tracks, particularly PRNs 6 and 8 moving away from the bunching region, the configuration improves considerably. A comparison between the two plots indicates the limits of the present constellation. The top plot shows that all seven of the prototype satellites are visible during the two hour period, and the bottom plot shows that only six of the 18-satellite constellation would be visible for the same period. However, the distribution of these six is far better, leading to a better configuration throughout the period, than is the case for the top plot. Also, satellite 4 rises before satellite 5 sets, so that there would always be at least five satellites visible. The conclusion is obvious what counts is not only how many satellites you can see but also how they are distributed spatially around you.
These plots were prepared using program MacGEPSAL, available from the Geodetic Research Laboratory, Department of Surveying Engineering, University of New Brunswick.
ON
rt
All five stations are monitor stations, tracking all GPS signals for use in controlling the satellites and predicting their orbits. The tracking is done by means of two-frequency receivers equipped with cesium oscillators. Meteorological data are also collected to allow for the most accurate evaluation of tropospheric delays. Positions of these monitor stations are known to a very high accuracy. Three of the stations (Ascension, Diego Garcia, and Kwajalein) are capable of transmitting data up to the satellites, including new ephemerides, clock corrections, and other broadcast message data, and command telemetry. One station (Colorado Springs) is the Master Control Station.
Tracking data from the monitor stations are transmitted to the Master Control Station for processing. This processing involves the computation of satellite ephemerides and satellite clock corrections. The master station is also responsible for controlling orbital corrections when any satellite strays too far from the assigned position. The Master Control Station would also initiate the necessary manoeuvres to have a dead satellite replaced by a spare.
4.10
Similar to P-code, with secret generating equation MESSAGE 1500 bit sequence of 1 steps at 50 bps
PHASE MODULATION TECHNIQUE
= "Multiply" Signal by -1
Canadian GPS Associates, May 1987
Both carriers carry the broadcast satellite message, a low frequency (50 Hz) stream of data designed to inform the user about the health and position of the satellite. These data can be decoded by the receiver and used for positioning in real time.
4.11
WAVE,
electromagnetic wave described as a function of elapsed time and distance travelled
receiver
at a fixed position
at a fixed time
phase = kx x = phase / k = phase A, / 2n
The distance travelled by the wave between the satellite and the receiver may be determined by either of two methods. If both the frequency and the exact instant of transmission from the satellite are known, then a measurement of elapsed time (between signal transmission and reception) allows one to compute the distance by
distance = velocity of propagation x elapsed time . This pseudo-range technique requires accurate time information.
The other method of distance determination is to measure the number of phase cycles over the transmission path. This can be done by mixing the incoming signal with the known signal frequency. Then, at any given instant, the propagation path distance is given from phase = kx = 2n/K x x = phase x 7J2n .
This method requires that the frequency of the satellite signal be accurately known and that the phase of the signal can be accurately tracked.
4.12
Measure At = 2p / c = two way time of travel Calculate p = c At / 2 One clock used to measure At Single user system
ONE-WAY RANGING (GPS)
^X/V/WAA
One-way ranging, as used in GPS, is more complicated. Each terminal, transmitter and receiver, is controlled by a separate clock. The transmitter clock generates the signal. The receiver clock detects when the signal arrives. Both clocks must keep the same time. An error in the synchronization of the two clocks of 1 jis creates an error in range of 30f metres. Since it is practically impossible to keep the two clocks perfectly synchronized physically, it is usually done mathematically. In general each clock will run at its own rate, keeping its own time. However, if the relationship between the two time bases defined by the clocks is known, then they can still be said to be synchronized. The problem, of course, is to find the relationship between the two time bases. GPS works only because we can assume that all GPS satellite clocks are synchronized (mathematically, not physically) among themselves. This is achieved by careful monitoring of differences in satellite clock timing by GPS ground station monitors, and including the resulting mathematical corrections in the satellite message broadcast to all users. Users can then assume that all the simultaneous GPS ranges measured by their receiver are related to the same (ficticious) clock at the satellite end of the range measurements. As we will see, synchronizing the receiver clock to this common fictitious satellite clock is then a 'piece of cake' (more or less).
4.13
CONTINUOUS WAVE
PULS]
at t 0 attj att
CONTINUOUS WAVE at transmitter
at receiver
att o
One of the disadvantages of the pulsed ranging is the sparse information content pulses arrive only at discrete intervals rather than continuously. GPS pseudoranging shares this sparseness: code epochs occur at discrete intervals rather than continuously.
Representation of the continuous wave one-way ranging concept is more complicated. The transmitter is now an oscillator that continuously generates a wave whose phase varies both in time and in space. We now need two sets of 'snapshots': (a) of the phase variation in time, viewed from a fixed point in space, and (b) of the phase variation in space, viewed at a fixed epoch in time. We assume the transmitter and receiver oscillators are perfectly synchronized in frequency and phase. In this case, plots of type (a) at both transmitter and receiver will be identical. To simplify the concept, consider plots of type (b) immediately after turning on the CW transmitter: i.e., 'snapshots' at three epochs as for the pulsed system. The range measurement made at t2 is given by the difference between the phase of the incoming signal received from the transmitter, and the phase of the internal oscillator in the receiver. The total phase difference involved in the range measurement includes the number of full wavelengths, or cycles of phase, between the transmitter and receiver,in the snapshot at t2. However, all we can directly measure is the fractional phase difference, between 0 and 360, or within one cycle of the signal. We must determine the integer cycle count some way other than by directly measuring it. This cycle ambiguity problem is the most serious drawback to CW range measuring systems. The main advantage is the continuous information density the CW wave is always available for measuring, not only at discrete epochs as for pulses.
H tt
V_____J
2
o
1
Q
PSEUDO-RANGE MEASUREMENTS
A GPS receiver can make basically only two kinds of measurements: pseudo-range, and carrier beat phase.
Pseudo-range is the time shift required to line up a replica of the code generated in the receiver with the received code from the satellite multiplied by the speed of light. Ideally the time shift is the difference between the time of signal reception (measured in the receiver time frame) and the time of emission (measured in the satellite time frame). In fact, the two time frames will be different, which introduces a bias into the measurement. These biased time delay measurements are thus referred to as pseudo-ranges. Pseudo-ranges are measured by a correlation detector that controls a delay lock loop which maintains the alignment (correlation peak) between the receiver-generated code replica and the actual incoming code. In simple terms, the code correlation concept involves generating a replica of the code sequence (the sequence of +1 and -1 values) within the receiver, and aligning this replica in time (correlate) with the incoming signal.
The pseudo-range measurement is then the delay which must be added to epochs of the receiver clock to keep the replica and incoming code aligned (correlated).
A rule of thumb for the precision of pseudo-range measurements (that is, the precision with which the correlation peak can be maintained) is 1% of the period between successive code epochs. For the P-code, successive epochs ^are 0.1 microsecond apart, implying a measurement precision of 1 nanosecond. When multiplied by the speed of light, this implies a range measurement precision of 30 centimetres. For the C/A-code, the numbers are ten times less precise, or a range measurement precision of 3 metres.
y2 = A2 cos2(cot + (|>) = A2[l + cos(2(Dt + 2<|>)]/2 . Since A(t) is the sequence of +1 and -1 values representing the code. A(t)2 = A2 is always equal to +1 and may be dropped. The resulting signal, y , is then pure carrier, but at twice the original frequency. Because the wavelength of the carrier is much shorter than the wavelength of either of the codes, the precision of carrier beat phase measurements is much higher than the precision of code pseudo-ranges. For the GPS LI carrier signal, the wavelength is about 20 cm. Using the rule of thumb that phase measurements can be made to about 1 % of the wavelength, this implies a precision of 2 mm.
The two main disadvantages of carrier beat phase measurements both involve the cycle ambiguity problem.
Obtaining the initial number of integer cycles of the carrier between the satellite and the receiver is very difficult, if not impossible. One way out is to take measurements which can be assumed to have the same (unknown) initial cycle ambiguity, and subtract one from another. More about this later. Maintaining an integer cycle count as the satellite-to-receiver range changes with time is something most good quality GPS receivers can do most of the time. However, for a variety of reasons, such as a noisy signal or a shaded antenna, any receiver will suffer cycle slips, or the loss of a coherent integer cycle count. In many cases painstaking postprocessing permits the detection and correction of cycle slips. However, the possibility of cycle slips restricts the use of carrier beat phase measurements for real-time applications.
4.16
electronics
Hardware
RF components
antennas
algorithms I
The hardware segment includes the electronic circuitry, RF components and antennas, and any peripheral equipment required to operate the receiver. The main features in this segment must include ruggedness, portability, field reliability, and ease of operation.
The software segment includes the obvious programming and data processing needed to convert the GPS messages into useful positioning or navigation information. This programming must allow adequate user interaction in order to take full advantage of the many features of GPS. It must also be adapted to ready field use and should be designed to provide useful system status and progress messages to the operator. Also included in the software segment are the programs developed independent of the GPS receiver which can evaluate such factors as satellite availability and precision confidence levels. Automatically included in the user equipment package is the technology which went into the various hardware and software components. This technology spans every area related to GPS, such as receiver design improvement, analysis and modelling of the various antenna and propagation effects and their incorporation into the data processing software, development of reliable communication link systems for short- to long-range differential GPS operation, and monitoring of the developing trends in equipment cost and performance.
4.17
GENE]
RF Section
Micro processor
Recording Device
Power Supply
Recording device: Either tape recorders or floppy disks are utilized to record the observations and other usefiil information extracted from the received signal. Power supply: Most receivers need low voltage DC power; only a few require AC power.
oo
TI
Tfr
03
E
H
U
a o
o a> ^ .S a> g
4^ P^
JS
ii
O <a
s =3 s
> o
O> *
"2 03
PM
-fi <> o
*& .
^^
vi & % O -
a *^ !5 H ^ 0
W.T
pQ
H t=l
a
c*
U
*
t=4
.o
H < H
u S o -*- t2 *
s s ^'S
43 Jg
18,
5S
OJ
X {=(
2 O t=1
&
pfi Cfr*
x^^
o o
S^TJ C3 O
b e^
< 03
S'g a
-"S 0) 4)
5
Q t=j
H
fe
>SS
S *-r S !2 fl ** h
4^
a> t gj
h*
g -o <
&g ^
fe
O U O) U
U *O
'S u
2"*
M t5 e
?| S
S .2 a) ^ *-
^ E ^ fi o S'T
** *S
v v
o ^ a
!
< 2 o
4.19
RELATIVE PO,
c ))
Receiver 2
AX, AY, AZ
Receiver 1
When solving for relative positions AX, AY, AZ, these errors are eliminated or reduced
RELATIVE POSITIONING
When a higher accuracy is required, then relative positioning must be used. In this operational mode, two antennas are positioned with respect to each other and two receivers (at the ends of the observed baseline) are operated simultaneously. The higher accuracy achieved in this mode occurs because several errors burdening the observed ranges are either identical or at least similar at both ends of the baseline. These errors then get either eliminated or, at least, significantly reduced when relative positions are sought. One particularly intriguing and promising mode of relative positioning, first suggested by Ben Remondi of the U.S. National Geodetic Survey, is relative semi-kinematic positioning. The idea is to use one stationary instrument and rove around with the other instrument. As long as no cycle slips occur in either receiver, accuracies of better than one cycle (20 cm) of the carrier beat phase signal can be continuously maintained in the relative positions between the stationary and roving receiver. This suggestion has two implications: kinematic applications may benefit from the much higher carrier measurement accuracies, rather than being limited to code measurement accuracies, and a much wider spectrum of GPS applications has been opened up. Setting out surveys, and aerotriangulation without ground control are two that are under investigation at present.
MULTI
4.20
MULTI-RECEIVER POSITIONING
Accuracy 0f results can be significantly improved when several receivers are deployed in a network positioning mode. A network is a geometrically stronger configuration than a single baseline because of redundancy the individual baselines connected into the network have to satisfy the conditions imposed by the geometry. This redundancy is used to control the effect of different errors, random as well as systematic, contained in the observations. We note that even with only two receivers available, it is preferable to connect desired baselines into a network, thus improving the accuracy of the positions.
When several receivers are deployed, one is faced with unique logistic problems concerning the pattern in which these are moved around and observing periods at individual stations. Under these circumstances, it is advisable to pay particular attention to the optimization of observing schedules to achieve the best accuracy with the most economical means.
(N Tf
(=)
t=l
t=l
4.22
Positioning accuracy
a = DOP a,
Geometry (Dilution of Precision)
Measurement accuracy
GOOD GDOP (ideal case) one satellite overhead 3 on horizon, 120 apart in azimuth
The effect of satellite configuration geometry is expressed by the dilution of precision (DOP) factor, which is the ratio of the positioning accuracy to the measurement accuracy, or
a = DOP a o ,
where aQ is the measurement accuracy (standard deviation), and
Where does the DOP come from? Recall the confidence ellipsoid discussed in Chapter 1. The standard deviation in one coordinate (say height) is represented by the distance from the centre to the surface of the ellipsoid along the local vertical direction (the height axis). According to the above equation, it is also equal to VDOP GQ . For horizontal positioning, in confidence ellipsoid terms, we want to express the 'size' of the horizontal ellipse: one reasonable measure is square root of the sum of squares of the two axes of the horizontal ellipse. This is HDOP-aQ. In general, any DOP is equivalent to the square root of the sum of the squares of the confidence region axes corresponding to the parameters being assessed.
4.23
GDOP
r
Fredericton D April 28, 1986 5 degree mask
Hours (UT)
I
PDOP
HDOP
TDOP
VDOP
DILUTION OF PRECISION
The dilution of precision, or DOP, is a measure of the geometrical 'strength' of the GPS satellite configuration. This will change with time, as the satellites travel along their orbits, and from place to place, since the satellite configuration is position-dependent. This figure shows five different DOPs plotted for Fredericton for April 28, 1986. The present constellation of seven prototype GPS satellites was used, with a mask angle of 5. These DOP plots correspond to the satellite configuration described in Figures 4.03,4.07, and 4.08.
We want the DOP to be as small as possible. Assuming a measurement accuracy of 10 m, and a DOP value of 5, our positioning accuracy is 50 m. If the DOP value is close to unity, then our positioning accuracy is close to our 10m measurement accuracy (a most blessed state!).
What can these DOP plots tell us? First we note that the GDOP plot has a spike at about 04:00. Recalling Figure 4.08, this occurs when PRN 3 sets at 4:07 and ends when PRNs 12 rises ten minutes or so later. What is the problem during this period? We note that PDOP, TDOP, and VDOP all contain the same spike, but HDOP does not. Because both TDOP (time) and VDOP (height) contain the spike, a common configuration weakness is suspected: when all satellites are too high in elevation, the solutions for time and for height become highly correlated, meaning that neither is well determined. Recalling the discussion of Figure 4.08, that appears to be the case here. Note that if our objective was horizontal positioning, we need not be concerned about this spike.
This figure was prepared with the help of program MacGEPSAL, available from the Geodetic Research Laboratory, Department of Surveying Engineering, University of New Brunswick.
'S
H P
CONSTRAINED SOLUTIONS
This figure has two purposes: it continues the discussion of the GDOP spike identified in Figure 4.23; it demonstrates the effect of a constrained solution on this spike.
The top right plot is a blow-up of the GDOP plot for Fredericton from Figure 4.23, during the period of the spike. The top left plot is a contour map of GDOP variations in space, for the instant the spike reaches its peak (at 4:08 UT). Note that to the north and east of Fredericton, the high GDOP value does not occur. Referring to Figure 4.08, we can see why. To the^north, PRN 12 will already have risen. To the east, PRN 3 will not yet have set.
What is a constrained solution? Let us say we have a reasonably good idea of our height, for example, and that our primary interest is in horizontal positioning. Then we could hold our 'known' height fixed and solve only for latitude and longitude (and time offset). However, a better idea is not to hold height fixed absolutely, but to constrain it to not vary more than a specified amount (equal to the uncertainty in our knowledge of the height) from the * known' value. The effect of constraining height (or any other parameter or parameters) is to improve the overall geometry the constrained height acts rather like an added 'observation'.
The bottom two plots show what happens to our spike when we know our height with a standard deviation of 5 metres and use this fact by constraining the solution for height. The GDOP spike has been cut almost in half, from just over 10 to just over 5. This improvement appears to be over the whole of the GDOP map area on the left-hand plots.
This figure was prepared with the help of program MacGEPSAL, available from the Geodetic Research Laboratory, Department of Surveying Engineering, University of New Brunswick.
c -,
A A
<S fS
<**>
^ fS
<**>
<n fS
:=>d jfe g
u
HH
J2 ,2
*j CQ -** CQ
l 1 J{ * 1111 S ^
-* CQ
'So s ^^ ^ cC '*
Q
C9
x
00 9s
* O tf>
HH HH^I H^
< z
^w HW
*
<S 9\ Tt ^\ t^ 9\
<D
2
03 'O
ON
9\
OS
&\
t9\
oo * ja o\ t3
^
08
OO XT
tt
c
QJ
a o
VO
O O
O
^ 0
O O
O
O O
10
O O
4J
oo
r- <s ^H oo
T3 5
*p
08
H CA
>H
W
C^
U * .ts A ^ 5
S -.
fl
O ^^
WD
PUi
c
C3
Dfi
<*
^ -g
S C5
O
-^ C
a
3 .3 'o o
H-l
I ! 8 5j 3 1 O H > H B
g i Ci V3
2 03 o CA
-1^ *
o w
<
2 o
a
IMPACT OF EXI
immediate impact
ON LAND Transit ISS (horizontal) ISS (vertical) EDM (engineering) EDM (cadastral) Theodolite Loran-C (nav) ON THE SEA Transit Loran-C Omega Syledis, etc. IN THE AM INS Omega VOR TACAN ILS/MLS X X X
4.26
X X X
X X
X X X X
If and when GPS receiver costs drop below the costs of owning and operating equipment for competing positioning systems, then GPS will erode the user community for those systems. It is worth noting that the prudent navigator will always prefer to have at least two different systems to mutually check each other's performance. However, a system like GPS, under complete military control, will have to overcome initial skepticism before it will be unequivocably accepted where safety of navigation has paramount importance.
MANY UJ
4.27
500 000
400 000
300 000
200 000
100 000
MILITARY GEODETIC
1975
1980
1985
1990
1995 2000
4.28
ACCURACY
Single Rx
Differential
5 m (3 m) 20 m (5 m) P-code 8 m (5 m) 50 m (50 m) C/A-code 1 ppm (0.01 ppm) Carrier (static) 10 cm (< 10 cm) Carrier (kinematic)
Receiver costs (high now, low later) Operating costs (no intervisibility requirement)
UNIFIED COORDINATE SYSTEM
Worldwide, all-weather, continuously available GPS and local geodetic coordinates simply related Forms the basis of land information systems
Widespread use of GPS by positioning professionals to replace or supplement existing techniques in society at large for many new applications
Canadian GPS Associates, April 1986
Low cost. No one would claim that GPS receivers are low cost today (1986). However, the real markets for GPS receivers await launching of the operational constellation. It is expected that at that time volume production will bring the cost of receivers down dramatically. However, receiver capital costs are not the only costs involved in positioning applications. The fact that GPS removes the need for mtervisibility between stations decreases the number of stations that must be established and permits the selection of locations that are convenient logistically, rather than for mtervisibility alone.
Unified coordinate system. This is perhaps the feature of GPS which will have the broadest impact. If (because of the first two features mentioned) GPS becomes widely adopted for positioning, then it would form the basis for unifying data from heterogenous sources in position-related information systems.
The consequence of these features is that GPS will likely be used not only by positioning professionals, such as surveyors and mappers, but also by large segments of the general public for new position-related tasks which were not previously feasible. The advent of a new emphasis in our society, that of the 'position based society,' may well result.
/mifa J Vd^^
/ ")
TO Jji
4.29
c? C iu)
1986
COST
1990
Ideal
$10
wrist locator
POWER SUPPLY
15 min (1 min)
milliseconds
Yes
0.01 ppm (few cm)
Yes
ACCURACY
1 mm
9
Canadian GPS Associates, April 1986
9 9
Some of the trends behind these speculations are as follows: Several manufacturers have developed specialized GPS Very-Large-ScaleIntegrated circuits (VLSI chips) which dramatically reduce the receiver parts count (and hence cost and power consumption). The ultimate goal of a standard GPS receiver-on-a-chip may well be achieved by 1990. A GPS receiver constructed for as little as $500 was predicted as early as 1978. While this low cost is a prerequisite for GPS success in the consumer market, positioning professionals will need and want a complete package, including warranty, servicing, data recording hardware, real-time and postprocessing software, and training. The receiver hardware cost will likely be a minor component in such a package. The U.S. Department of Defense has contracted the development of hand-held receivers for military use. Commercial exploitation of this development is sure to follow. It is becoming standard to provide carrier beat phase measurements in GPS receivers or to use them at least internally in the receiver. The initial military design specification that a GPS receiver need not be capable of tracking more than the required minimum of four satellites is being replaced with receivers capable of tracking up to the maximum possible number of visible satellites. Dual frequency receivers not dependent on access to the P-code are expected to become widely available over the next few years. Real-time communication links, and receivers capable of using them, are beginning to emerge, making possible real-time relative GPS positioning.
CHAPTER FIVE
5.00
ORBIT DESCRIPTION IDEAL WORLD coordinate systems forces acting on satellites Kepler's laws Keplerian orbital elements position and velocity ORBIT DESCRIPTION REAL WORLD non-central earth gravity third body, tides, solar radiation, drag
ORBIT DETERMINATION ephemeris prediction postcomputed ephemerides active control system concept
Orbital information comes in two flavours: predicted and postcomputed. The methods, sources, and accuracies of each type are discussed. In particular the orbital information broadcast by the GPS satellites is described in detail.
Recommended further reading:
Broughton, D.W., R.I. Esai and P.J. Farrow (1984). Brouwer, D., and G.M. Clemence (1961). Buffett, B.A. (1985). Delikaraoglou, D. and R.R. Steeves (1985). Henriksen, S.W. (Ed.) (1977). Kaula, W.M. (1966). Mertikas, S.P., D. Delikaraoglou and R. Santerre (1986). Nakiboglu, S.M., B. Buffett, K.P. Schwarz, EJ. Krakiwsky, B. Wanless (1984). OToole, J.W. (1976). Remondi, B.W. (1985d). Russell, S.S. and J.H. Schaibly (1978). Steeves, R.R., D. Delikaraoglou and N. Beck (1986). Swift, E.R. (1985). van Dierendonk, A.J., S.S.Russell, E.R. Kopitzke and M. Birnbaum (1978).
5.01
INFORMATION
POINT POSITIONING
^phemeris \
rue orbit dr = satellite position error dp = equivalent range error dR = resulting station position error RELATIVE POSITIONING RULE OF THUMB
'Jphemeris
*True orbit
db
Accurate orbital data is essential for the reduction of the observations collected by the GPS receivers in order to obtain the user's position, which is usually the final product sought in GPS positioning. In the case of point position computations, there is clearly a stricter requirement for precise orbital information since any error in the ephemeris data will be directly transferred into an equivalent range error, and hence it will bias the final position results in a systematic way. For relative positioning, the effects of any uncertainties in the ephemeris data on the final position results are somewhat alleviated because, for baselines which are short compared with the 21000 km GPS satellite altitudes, satellite orbit errors tend to cancel out. The equation that defines the vector error db in the baseline vector b, as a function of the vector errors di^ in the positions of the i satellites used to determine b is
with pj denoting the range to the ith satellite, and Cj denoting the unit vector to the
\th satellite. Clearly db depends not only on the magnitudes of p, b, and dr but also on their directions. As well, the baseline solution depends on measurements from several satellites over an extended time period, and the satellite position errors di^ involved will, in general, be different. The effect of one particular satellite orbit error, say drk, is given by
db e'kPk = b div Taking the length of the vectors in this equation, we obtain the rule of thumb shown in the figure. This rule of thumb is an upper bound for the magnitude ||db|| of the baseline error as a function of the ephemeris error. Evidently, on baselines a few tens of kilometres long, the orbit error has a very small impact: accuracy of a few centimetres can be obtained with present techniques. On longer baselines, of the order of several hundreds of kilometres, the impact can be large enough to be a major obstacle at present for the general use of GPS for precise geodetic and geodynamic applications.
20 2
1
10 100 0.1
10
5.02
APPROXIMATE
a Calgary-Fredericton baseline, April 28, 1986
Calgary Frederic
4
Canadian GPS Associates, April 1986
UT hours6
The actual computation of GPS positions requires accurate orbital information. Other tasks involved in the use of GPS also require knowledge of satellite positions, but at a much lower accuracyperhaps only within a few hundred kilometres. Two of these tasks are: Planning the use of GPS. Knowing when (what part of the day) to use GPS receivers requires orbital information. During the present period of incomplete coverage, when the satellites are visible only during short 'windows' each day, this is particularly important. Signal acquisition. The first thing any GPS receiver must do is acquire the satellite signals to be used. If nothing whatever is known about the satellites to be tracked, then the receiver must perform a 'cold start' and do an 'all sky search/ This involves selecting each of the 32 possible C/A-codes, and searching over all possible Doppler shifted frequencies for each signal, attempting to lock on to any one of them. Provided that at least one satellite is visible, this search may take as long as 30 minutes. This time can be reduced considerably if the receiver knows which satellites are visible, which limits the search to only a few of the 32 possible codes. If the satellite positions are known as well, then estimates of the Doppler shift can be predicted, limiting the frequency band over which the search must be performed. In this way, the search can be reduced to a few minutes. Once one GPS satellite signal has been acquired, the message from that satellite provides the receiver with the approximate location of all of the other GPS satellites ('almanac' data), so that their signals can be acquired very rapidly. The approximate satellite positions computed for planning or receiver aiding are known as satellite alerts. For planning purposes, there are several useful graphical representations of the results: a linear plot against time of the elevation and azimuth (with respect to a specified location) for each satellite. a polar plot, in which elevations and azimuths as a function of time for each satellite, with respect to a specified location, are predicted and plotted. a track plot, in which the satellite subpoint (the point where the line joining the satellite to the geocentre pierces the earth's surface) is plotted on a map. a visibility plot, against time of day, of the number of satellites which are visible from a specified location, or of the period during which each satellite is visible from a specified location (shown here). a Dilution of Precision (DOP) plot, against time of day, showing the changing geometrical strength of the visible GPS satellite constellation (shown here). Area contours of DOPs at the specific tune epoch are also used.
5.03
.ATELLITE SELECTION
SATELLITE SELECTION
Using the present GPS constellation of seven satellites, it is possible that as many as six satellites will be above the observer's horizon (be visible) at once. With the full operational GPS constellation, as many as ten satellites may be visible at any one time. Most GPS receivers do not have the capability of tracking all visible satellites at once. In many cases, only four satellites can be tracked simultaneously.
The question then arises: How to select the 'best* satellites to be tracked? The usual criterion to be used is to choose those satellites which will result in the most accurately determined receiver positions. One measure of how accurately positions are likely to be determined is the relevant Dilution of Precision factor, or DOP. The satellite selection problem then reduces to finding those satellites which provide the lowest DOP.
The most straightforward use of GPS is for point positioning in four dimensions (latitude, longitude, height, and time). In this case the appropriate DOP is the Geometrical Dilution of Precision, or GDOP. For the case of a receiver capable of tracking four satellites simultaneously, it has been shown that GDOP is related to the volume V of the spatial tetrahedron iformed by the unit vectors from the receiver to four satellites. GDOP is inversely proportional to this volume V, hence the satellite selection that maximizes V also minimizes GDOP. This has provided a simple algorithm that can be used for selection of the best available satellites to be tracked. This algorithm is simple enough to be programmed into software contained within the receiver itself.
However, for receivers capable of tracking more than four satellites simultaneously, but still less than the maximum number of visible satellites, this simple algorithm does not apply. Similarly if GDOP is not the most appropriate DOP for the GPS application (for example, for two-dimensional horizontal positioning, or for time transfer applications), then again this simple algorithm.will not necessarily yield the best satellite selection. Finally for relative positioning and network surveying applications, the satellite selection problem is also more complicated. It has been found that the optimal satellite selection for point positioning will also be the optimal selection for relative positioning, when the baselines are short. However for long baselines (say 5000 km), this is no longer the case. The optimal satellite selection changes with time, as satellites rise and set, and the spatial tetrahedra which they form change in shape. One strategy is to switch satellites whenever a new optimal selection occurs. However, there are other criteria which may be more important, such as maintaining continuity of tracking the same satellites for longer periods. In these cases, the satellite selection problem has not yet been solved.
MOTION OF iS AROUND
AN
THE IDEAL EARTH: Perfectly uniform spherical shape No atmosphere
5.04
mass of the earth mass of the satellite universal gravitational constant distance of the satellite from the geocentre force of the attraction
F
GPS
Geostationary
PRACTICAL LIMITATIONS:
a perfectly radial gravitational field; and no atmosphere. After the satellite is launched, it begins its unrestrained revolution around the earth. This motion in a radial gravitational field is governed by Newton's gravitational law: F =GMm/r2
where F is the force of attraction, M is the mass of the earth, m is the mass of the satellite, r is the distance of the satellite from the geocentre, and G is the universal gravitational constant.
The practical implication of Newton's gravitational law is that gravity never becomes zero. No matter where a satellite is in space, there always is a force pulling it toward the earth, although the earth's pull may become very small compared to that of other bodies in deep space.
5.05
3. transfer orbit
4. final orbit
shuttle takes satellite up to parking orbit Payload Assist Module spins satellite up to 70 rpm and boosts satellite into elliptical transfer orbit
1 2 3-4 5 6-7 8 9 10 11 12 13 14 15
rocket rocket shuttle rocket shuttle rocket rocket rocket shuttle shuttle rocket shuttle rocket
Feb89 May 89 Jun89 Aug89 Sept 89 Nov89 2D Feb90 May 90 May 90 May 90 Aug90 Oct90 Nov90
16-17 18 19 20 21 22 23 24 25 26 27 28
Achievement of worldwide 2D coverage by Nov 1989 assumes that four of the present prototype satellites are still functioning. Achievement of 3D coverage by February 1991 assumes that two of the present satellites are still functioning.
5.06
!Y STEM,
reenwich Meridian
Equator Satellite motion is described in RA-system Terrestrial positioning refers to CT-system (^ ,X2,X3) The two are related by : Greenwich apparent sidereal time (0) Polar motion components (Xp,yp)
Canadian GPS Associates, May 1987
COORDINATE SYSTEMS
The dynamics of satellites is governed by Newton's law of motion which states that the acceleration of the centre of mass of a free body is proportional to the force applied to it. The proportionality constant here is the mass of the body. The position, velocity, acceleration, and force components should refer to a special coordinate system called inertial. An inertial coordinate system, loosely speaking, is one which is either stationary or in uniform motion (with no accleration) in space. Our first task in describing satellite motion, therefore, is to define an appropriate inertial coordinate system. We conveniently adopt the right ascension (RA) system, which fulfills our reference frame requirements for short time periods. The orientation of the RA-system relative to the fixed stars changes by about 1" per year, due to the forced periodic motion of the earth's rotation axis and the plane of the ecliptic. However, this motion, called precession and nutation, is known very well, and the orientation of the RA-system at different epochs may be traced back.
The RA-system is not a convenient reference frame for terrestrial positioning purposes. For the latter, we require a coordinate system which is 'fixed to the solid earth' so that a stationary point on earth will have time-invariant coordinates. The conventional terrestrial system (CT-system), which is also known as the average terrestrial system, is one such system. The first axis of the CT-system passes through the intersection of the Greenwich meridian and the equatorial plane. The third axis of the CT-system is defined by the Conventional International Origin (CIO), which is the average position of the earth's rotation pole for the years 1900 to 1905. The second axis is orthogonal to the first and third axes. The RA-system and CT-system are related as
R = R2(-xp) R^) R3(8) r.
Here, 0 is Greenwich apparent sideral time; x and y are the components of polar motion; R and r are the CT-system and RA-system coordinates, respectively; and Rj (i= 1,2,3) are rotation matrices.
Until July 1985, the realization of the CT-system used to describe the GPS broadcast ephemeris was the World Geodetic System 1972 (WGS-72). From July 1985 to December 1986 the WGS-84 system is used for all ephemeris prediction computations, but the broadcast ephemeris information is then converted back to WGS-72. Since January 1987 the WGS-84 system has been used throughout.
5.07
ACID ELLIT
orbit SUN
MOON
F r
Geocentric gravitational attraction of the earth Non-central gravitational attraction of the earth Third body effects of the moon and the sun Solar radiation pressure
Atmospheric drag force Tidal forces (not shown) Magnetic forces (not shown)
Other forces
Fa
(a) gravitational attraction of the earth, (b) gravitational attractions of the sun, moon, and planets (called the third body effects), (c) atmospheric drag force, (d) solar radiation pressure, both direct and albedo effects, (e) magnetic forces, and (f) the variable part of the earth's gravitational field arising from tidal and other deformations of the solid earth and oceans.
Among the forces listed above, the most dominant is the gravitational attraction of the earth. The earth's gravitational attraction may be further separated into two parts:
(a) the main part, called the central part, is the radial gravitational attraction, and (b) the non-central part, which is the rest.
The central part of the gravitational attraction is three orders of magnitude larger than the non-central part and all the other forces combined. Therefore, the main characteristics of satellite motion are determined by the central gravitational field, while the remaining forces introduce some small variations into the motion. These small variations and forces causing them are called the perturbations and perturbing forces, respectively.
5.08
1L
THE ORBIT IS AN ELLIPSE WITH ONE OF THE FOCI LOCATED AT THE EARTH'S CENTRE OF MASS
To get lower inclinations requires "parking orbit" with second stage launch at equator crossing => heavy, complex and expensive satellites
5.09
Apogee
Perigee
This law states that the speed of a satellite in an elliptical orbit is not constant. Practical use can be made of this law. It is desirable for spy satellites to travel low and fast over their target area consequently, they are placed in highly elliptical orbits with perigee over the target area. It is desirable for communication satellites to have as wide a coverage area as possible, and for as long as possible. Most communications satellites are placed in geostationary orbits over the equator, in which they appear to be stationary in the sky. Unfortunately, this does not provide very good coverage at high latitudes. Therefore other communications satellites are placed into highly elliptical orbits with apogees over their coverage areas. They are not geostationary, but (from Kepler's second law) they travel very slowly near apogee, and hence can be used for most of their orbital period. The lifetime of a satellite before it decays into the atmosphere is determined by how quickly it is slowed down by atmospheric drag. This slowing effect is proportional to the satellite's speed, and is felt mostly at perigee, where the atmosphere is densest and the speed highest (Kepler's second law again). Hence the initial perigee height determines the lifetime of the satellite. Typical satellites and their lifetimes are:
Satellite
Perigee fkm)
Apogee fknri
Lifetime
15 3 12 1000 106
5.10
THE
(Orbital Period)2 (Semi-major Axis)'
Constant
JS 10
w*
Geosynchronous
" 10
J-
:AL IMPLICATION
Two satellites with the same semi-major axis length have equal orbital periods regardless of eccentricity
Satellite 2 Satellite 1 Geocentre
It is a direct consequence of this law that the GPS satellites, which have an orbital radius equal to about four earth radii, have a longer orbital period (about 12 hours) than other commonly used geodetic satellites, like Transit or LAGEOS, which are in much lower orbits.
IT
5.11
Satellite
Perigee
Vernal Equinox
X
The stationarity of orbital plane, perigee, size, and shape of the orbit, together with constant period, lead to a realization that the Keplerian motion is fully describable in terms of six parameters, of which only one is a function of time.
The choice of these parameters is not unique. A particular set of parameters called the Keplerian elements is commonly used in satellite positioning. The definitions of these are as follows:
(a) Right ascension of the ascending node (Q). The geocentric angle between the nodal directions and vernal equinox measured in the equatorial plane. (b) Inclination (i). The angle between the equatorial and orbital planes. (c) Argument of perigee (co). The angle between the nodal and perigee directions measured in the orbital plane. (d) Semi-major axis of the elliptical orbit (a). (e) Eccentricity of the orbit (e). (f) An element describing the position of the satellite on the orbital ellipse. One of the anomalies is normally used.
The Keplerian elements, Q and i, define the orientation of the orbital plane in space; co defines the location of the perigee on the orbit; and a and e define the size and shape of the orbit. The choice of Keplerian elements as the constants of motion is inappropriate in some cases. For example, for circular orbits we have e=0, and co and / are undefined.
ATELLITE
)SITION IN TIME
5.12
q=
a cos E - a e a V l - e2 'sinE
0
a q = ( 1-encos E )
- sinE - e2 cos E 0
position = R
and
velocity = R
where
The orbital period T of the satellite in terms of n is, as we have already seen,
Satellite position and velocity vectors in the RA-system can now be expressed in terms of Keplerian elements and time, as shown in the figure.
5.13
Satellite at t + At
p At
ong track Radial Satellite at t
r = R 3 (9)R x q q.
The topocentric range and range-rate of a satellite from a tracking station with geodetic latitude (p, longitude A,, and ellipsoidal height h are readily obtained as
ap/at=(ar/ao T
where,
P = Hr-R||
/P,
aq / at = n a /(1-e cosE)
ee2 sin2(p)1/2
N = ae /(I is the radius of curvature of the earth ellipsoid in the prime vertical direction; ae and ee are the semi-major axis and eccentricity of the earth
ellipsoid; and a6/at = CDC 0.729 211 514 7 x 10~4 rad/sec is the earth's rotation rate.
Various satellite observations, like range, Doppler effect, and phase, may be expressed in terms of p and/or dp/at.
5.1
10
Hours
1S
OUT-OF-PLANE PERTURBATION
ALONG-TRACK PERTURBATION
200^
Canadian GPS Associates, May 1987
Non-central part of the earth's gravitational attraction. The density distribution within the earth departs considerably from radial symmetry.
Consequently, the gravitational acceleration varies with latitude and longitude in addition to the geocentric distance. The non-central gravitational force is primarily due to the equatorial bulge of the earth. Also, the earth at the south pole is slightly more depressed than at the north pole resulting in a 'pear shape.' There are other lesser mass anomalies contributing to the non-central gravitational attraction. The attraction force attenuates rapidly with altitude. Therefore, only a few terms in the earth's gravitational potential are sufficient for modelling GPS orbits.
The equatorial bulge leads to 2 main perturbations of the ideal Keplerian orbit: Rotation of the orbital plane in space in a direction opposite to the satellite's eastwest motion. This motion depends on inclination according to cos i. Thus an orbit with less than 90 inclination will precess to the west. For a truly polar obit (90 inclination), this effect is zero, and for an equatorial orbit (approaching 0 inclination), a maximum. This perturbation can be modelled by }-dot. Rotation of the major axis in the orbital plane. This perturbation depends on inclination according to (1-5 cos2i), resulting in rotation in a direction which is forward for equatorial orbits, zero for an orbital inclination of 63.4, and backward for polar orbits. This perturbation can be modelled by co-dot.
The perturbations of GPS orbital arcs of four hours' duration, due to the non-central earth's gravitation effects are: Effect of Effect of Orbital element Equatorial bulge (m) Higher order terms (m)
a e i Q co + M
20 5 5 3 40
This figure shows the growth of the effect of the higher-order terms over 12 hours.
5.15
HI]
Oil
EFFECT
THE DIRECT ATTRACTION OF THE MOON IS THE MOST SIGNIFICANT THIRD BODY EFFECT.
ALTHOUGH THE SUN IS MORE MASSIVE, IT IS ALSO MUCH FARTHER AWAY AND, THEREFORE, ITS EFFECT IS LESS THAN THAT OF THE MOON. PLANETARY EFFECTS FOR GPS ARE NEGLIGIBLE.
INDIRECT EFFECT
THIRD BODY DEFORMS EARTH (I.E., TIDES) WHICH IN TURN CHANGE THE EARTH'S GRAVITATIONAL FORCE AT THE SATELLITE.
The perturbations of GPS orbital arcs of four hours' duration, due to the lunar third body effects are:
Orbital element
a e i Q
CG + M
Since the value of M* / |r*|3 for the sun is 0.46 that for the moon, the solar third body effect is about half that of the moon. For the planets, the masses M* are far smaller than for the sun, and distances |r*| much greater than for the moon; the planetary values for M* / |r*|3 are so small that their third-body effects on GPS satellites are negligible.
Tides. The gravitational attractions of the moon and the sun also have an indirect effect on satellite motion. The solid earth and the oceans undergo significant tidal deformations at primarily the semi-diurnal and diurnal frequencies. The tidal deformations cause a periodic re-arrangement of the earth's mass which in turn alters the earth's gravitational attraction. This potential change depends on the deformation characteristics of the earth, and it is closely related to the third-body potential. These effects are insignificant for GPS satellites over orbital periods.
>LAR RAJ
5.16
ACCELERATION DEPENDS ON SOLAR RADIATION AT SATELLITE (SHADOW EFFECT) SATELLITE MASS AND SURFACE AREA EXPOSED TO SUN WHETHER SATELLITE IS IN SHADOW OR NOT TO GET 2 m GPS ORBITS, RADIATION MUST BE TAKEN INTO ACCOUNT
ATMOSPHERIC DRAG ACCELERATION DEPENDS ON ATMOSPHERIC DENSITY SATELLITE CROSS SECTIONAL AREA AND MASS SATELLITE SPEED EFFECT IS PRIMARILY ON APOGEE HEIGHT; PERIGEE IS CONSTANT. AT THE HEIGHT OF GPS SATELLITES, ATMOSPHERIC DENSITY IS SO LOW AS TO MAKE DRAG NEGLIGIBLE.
Canadian GPS Associates, April 1986
Orbital element
a e i ft cfl + M
Atmospheric drag. The atmospheric drag is a non-conservative force which arises from the friction between the satellite surface and the surrounding atmosphere. The consequence of frictional dissipation is that the satellite loses energy which results in secular (non-recoverable) changes in a and M. The atmospheric drag force per unit mass on a satellite varies according to p v2 A / m where p is the atmospheric density at the satellite, v is the satellite's velocity relative to the atmosphere, A is the satellite cross sectional area, and m is the satellite mass. The atmospheric density at GPS satellite altitudes, ~20 000 km, is so small that the drag effects are negligible.
5.17
p (t)
+>( EDIT
i CORRECTIONS)
| SMOOTH" ^_
KALMAN FILTER
OUTPUT)
Canadian GPS Associates, April 1986
EPHEMERIS PREDICTION
The broadcast ephemeris information transmitted by the GPS satellites is computed and controlled by the GPS Master Control Segment. Ephemeris prediction is a two-step process: (a) an off-line least-squares fit to produce a reference ephemeris, and (b) an on-line first-order correction using a Kalman filter and additional measurements provides current estimates of the satellite states which, in turn, are used to predict future states.
The reference ephemeris is an initial estimate of the satellite trajectory, computed from about one week of data from the five GPS monitor sites (Colorado Springs, Diego Garcia, Ascension, Kwajalein, Hawaii). It is used to permit accurate linearization of the orbital model in the Kalman filter. Errors in the reference ephemeris of order 100 m would result in linearization errors well within the GPS ephemeris prediction specification of 1 m per day. Linear perturbations about this reference ephemeris are computed by the on-line Kalman filter. The reference ephemeris also is used to compute the time dependent state transition matrix (partial derivatives) required by the Kalman filter. The on-line ephemeris prediction is driven by pseudo-range and integrated Doppler measurements (nominally at 6 sec intervals) from each satellite in view at each monitor station. After editing and time tag correction, these measurements are compared with equivalent quantities computed using the reference ephemeris and misclosures formed. These misclosures are corrected for various known biases, such as ionospheric delay, tropospheric refraction, relativistic effects, and satellite and monitor station antenna phase centre offsets. The range measurements are then smoothed to give a set of smoothed pseudo-range data every 15 minutes for each satellite-station pair. The integrated Doppler measurements are not smoothed but sampled every 15 minutes. These data are processed by a Kalman filter estimator to produce estimates of the satellite position and velocity states. The output of the Kalman filter is 6 orbital element perturbations, the satellite clock frequency offset and drift rate, 3 solar radiation pressure parameters per satellite, monitor station clock bias and frequency offset, tropospheric scaling bias for all monitor stations, and 3 earth rotation component residual states.
The orbital element perturbations are used to correct the reference ephemeris, and an extrapolated ephemeris is predicted for upload to the satellite. Similarly, the satellite clock states are propagated throughout the next prediction span, and a time polynomial representation is prepared for upload. Predicted ephemeris errors of about 1 m radially, 7 m along track, and 3 m cross track are expected, once GPS is fully operational. These predictions are reformatted into GPS message parameters and uploaded into the satellite memory at least once a day. The normal upload contains sufficient data for 26-hour orbits. Special 6-hour uploads are used to support testing requirements.
5.18
PRECISE EPHEMERIS
IL NETWORKS IN
POSTCOMPUTED EPHEMERIDES
The predicted broadcast ephemeris will not always satisfy the needs of all users. Some will use receivers which do not decode the broadcast message (for example, so called Modeless' receivers). Some will want more accurately determined orbits than are provided by the broadcast ephemeris. The alternative to the broadcast ephemeris is a postcomputed ephemeris. Several sources for postcomputed ephemerides are presently available, or have been proposed.
The so-called precise ephemeris is computed by the U.S.Naval Surface Weapons Center (NSWC). This orbital determination operation is similar to the operation used to determine the reference trajectories for the broadcast ephemeris computation. However, data from the five GPS monitor stations are augmented by data from sites in Australia, Seychelles, England, and Argentina. A more extensive and international network is planned for the future. The U.S. National Geodetic Service has announced that they will offer reformatted precise ephemeris data to civilian users, under two standard GPS orbit data product subscription services: one designed for large computers and magnetic tape distribution, and the other for small computers and telecommunications distribution. FORTRAN 77 routines to handle the data will also be provided. The proposed subscription costs would be $US 100 per week, or $US 2000 for six months.
The U.S. Geological Survey has contracted with the University of Texas to produce postcomputed ephemerides, beginning in September 1986, based on data from three stations ^Austin, Texas; Westford, Massachusetts; and Mojave, California). A Canadian station (Yellowknife) is to be added to the network in early 1987. A German station (Wettzell) may also be added. The initial operations are for six months, but may be continued for up to three years. TI 4100 receivers, controlled by IBM PC/XT computers and cesium clocks, are used at each station. Orbits are computed at the University of Texas. Alternative orbital predictions are planned by the Canadian Geodetic Survey.
Three dual frequency Macrometer receivers have been in operation since mid-1985 at the POLARIS stations at Westford, Massachusetts; Fort Davis, Texas; and Richmond, Florida. Orbits based on this data are believed to be computed using the PEP program at MIT.
Litton Aeroservices offer a commercial GPS orbit data product. The data source for this orbit has not been revealed by Litton.
Real-time data distribution will depend on cost-effective mobile satellite communication services being available. This is expected early in the 1990s.
This conceptual diagram illustrates the various data and control communications which may be involved.
5.20
T EPHEMERIi
M0
An e
o
o
CO
Rate of inclination
us Corrections to argument of latitude Corrections to orbital radius
, c IS
t oe
Corrections to inclination
Ephemeris reference time
Each broadcast ephemeris parameter set is intended for use only during the one hour period to which they refer. Extrapolation of the broadcast ephemeris well beyond this period leads to an exponential error growth which can affect significantly even differential positioning accuracies, especially for large station separations. The actual accuracy of the broadcast ephemeris information has not been adequately assessed to date. Reported results of several GPS campaigns have implied (from comparisons of GPS derived station coordinates at known locations) an accuracy at the level of 20 to 50 metres. However, such estimates are affected by the 'freshness' of the ephemeris and the area of application.
5.21
Y
node
@t 0 bs
node
@t 0 bs
-toe)
Y
GASTweelc ^^^^^^^^^^
GASTpbs "H^
(tnbs~ tweek)
zero meridian
@ tweek Q0(broadcast)
GASTobs= GASTweek+
where &0bs is the angle between the vernal equinox and the ascending node (the right ascension of the ascending node) at tobs, and GASTobs is the angle between the vernal equinox and the Greenwich meridian (Greenwich Apparent Sidereal Time) at tobs. However the broadcast ephemeris does not supply values for Q and GAST at tobs. This complicates things slightly. As shown in B of the figure, the value we use for Q is for the ephemeris reference time t^. Since we are given the rate of change of right ascension Q-dot, we can use
"obs = "oe +
Q dot
Q-dot (the orbit plane precession rate) is typically a few millidegrees per hour. New values for Q-dot and t^ are given hourly in the broadcast ephemeris.
As shown in C of the figure, the value we use for GAST is for the beginning of the week (Saturday/Sunday midnight in the GPS time system), tweek. Since a value
coe is approximately 15 degrees per hour (the earth's rotation rate). The times tobs and toe are the elapsed time since tweek, so we can take it as having a zero value. We must still obtain values for Q^ and GASTweek (actually their difference). As indicated in D of the figure, the right ascension parameter in the broadcast ephemeris is QQ(broadcast) = QQe - GASTweek, and not, as might be expected, equal to Q^ alone. It is therefore mislabelled as a right ascension. Making all the above substitutions (and setting tweek = 0), we finally obtain
5.22
IT DEi
16 PARAMETERS NEW VALUES EACH HOUR
C
uc' Cus
rc'
rs
ic' is
Satellite
Reference
Perigee
*Y.
Ascending Node
The ascending node. This is the intersection of the orbital plane, and the equatorial plane, as the satellite passes from the southern hemisphere into the northern hemisphere.
Perigee. This is point at which the satellite most closely approaches the earth. The angle between the ascending node and perigee, measured in the orbital plane.
The reference position. This is position of the satellite at the reference time toe (the time for which MQ, iQ, and Q,Q are given). The satellite position. This is what we want to determine. It is separated from perigee by the true anomaly /, and from the ascending node by the argument of latitude u = co + /.
In order to transform the orbit parameters to CT-system coordinates, three rotations between four coordinate systems are used:
The satellite position is defined in the plane of the satellite orbit in the (x lfl , y1", z" f ) coordinate system, with x1" toward the satellite, y'" in the orbital plane, and z'" normal to the orbital plane.
A rotation of the (x ltf , y !tl , z l f t ) system about the z lfl -axis (3-axis) from the satellite to the node (through the argument of latitude u) produces the (x", y", z") system, with x" aligned to the ascending node, y" in the orbital plane, and z" = z it?
A rotation of the (x lf , y", z") system about the x"-axis (1-axis) from the orbital plane to the equatorial plane (through the inclination angle i) produces the (x', y1, z1) system, with x1 = x", y1 now in the equatorial plane, and z1 aligned to the earth's rotation axis.
A final rotation of the (x1, y1, z1) system about the z'-axis (3-axis) from the ascending node to the Greenwich meridian (through the longitude of the ascending node at the observation epoch, A, u ) leads to coordinates (x, y, z) of the satellite in the CT-system. These coordinates can be computed for any time during a given observation span.
5.23
iATELLITE
CONSTANTS (these values must be used) gravitation constant |i = 3.986008 x 1014 m3/ sec2 mean earth rotation rate coe = 7.292115147 x 10"5rad/sec n =3.1415926535898 TRUE ANOMALY / k time since reference epoch t k = t - t oe mean anomaly at t k M k = M 0 + (VjT/ Va" 3 +An) t k solve Kepler's equation iteratively for eccentric anomaly Ek M k = Ek- ejsin Ek -1 true anomaly /. / k = tan" cos Ek- e JIG OF LATITUDE u. , RADIUS rv INCLINATION i
Jt
&L
k= 0 + 4+
k=
CLcos 2
l\s
the argument of perigee co, the true anomaly f^ and the correction terms with coefficients Cuc and Cug.
Compute the orbit radius r^. using the expression for elliptical radial distance, and correction terms with coefficients Crc and Crs.
Compute the orbital plane inclination i^, using the inclination iQ, at reference time toe the linear change in inclination since reference time, and correction terms with coefficients C^c and C|s. Compute the longitude of the ascending node A^, by adding the right ascension parameter, &o(broadcast), the change in Greenwich Apparent Sidereal Time between the beginning of the week and the reference time toe, and the change in longitude of the ascending node since the reference time t(oe> Compute the CT-system coordinates by applying the three rotations (through iu, and A,^,) described previously.
CHAPTER SIX
SIGNAL STRUCTURE
6.00
Why is the GPS signal so complicated? Definitions SIGNAL Code modulation technique Code generation Message format and content
What leaves the satellite antenna ? GPS signal tolerances Speculations on Selective Availability
FROM SATELLITE TO RECEIVER
SIGNAL STRUCTURE
In this chapter we examine the structure of the signals sent out by the GPS satellites. Many of the strengths of GPS depend on various features of this signal structure. We start by considering why the structure is so complicated. In order to describe the signal, we need to understand some unfamiliar terms, so these are defined next.
The GPS signal consists of a carrier upon which three modulations are impressed: the C/A-code, the P-code, and the broadcast message. The method of modulating the carrier by binary biphase modulations is described. The two codes are pseudo-random noise codes. These are generated using a device called a tapped feedback switch register. The operation of this device, and how it leads to the C/Aand P-codes, is considered next. The message format and content conclude our discussion of the various GPS signal components.
The methods by which these signal components are combined to form the signal which actually leaves a GPS satellite are considered by tracing the signal through a simplified schematic of the satellite's inner workings. The specifications for the officially-stated tolerances on various aspects of the signal generation are reviewed.
The policy of the U.S. Department of Defense is to intentionally degrade the accuracy obtainable from GPS by civilian users, once the GPS system is fully deployed. This degradation, officially termed Selective Availability is reviewed, and some speculations are made on what parts of the signal generation might be involved. Finally, one of the features of the signal as it is received on the ground, is its Doppler shift, due to the relative motion between satellite and receiver. We examine the features of the Doppler shift and its derivatives.
Complete details of the GPS signal structure are specified in the Interface Control Document [Rockwell, 1984]. A detailed review of the signal structure is given by Spilker [1978].
REAL-TIME POSITIONING
Simultaneous measurements from many satellites how to identify signal from each satellite? Unambiguous range measurements signal modulation to get time delay Real-time satellite positions broadcast message containing ephemeris
High accuracy velocities (10 cm/s) and positions (10 m). Continuous centimetre-per-second-level velocities require Doppler shift measurements on a centimetre-wavelength (microwave) carrier. Continuous metre-level positions require metre-wavelength modulations, such as the 10 MHz P-code. Metre-level positions also require compensation for ionospheric refraction, most easily accomplished by making dual-frequency corrections.
Pollution of the electromagnetic spectrum. New systems like GPS do not easily gain permission to use parts of the electromagnetic spectrum. For military use deliberate electromagnetic pollution (jamming) is something to worry about. The modern communication-technology solution to both problems is to spread the signal over a wide bandwidth, in such a way that it can share the same part of the spectrum with many other users and,at the same time, be impervious to jamming. GPS uses such a spread-spectrum technique, implemented by the choice of modulating code, also used for the range measurements. Mixture of military and civilian users. GPS will have both the military Precise Positioning Service (PPS) using the P-code, and the civilian Standard Positioning Service (SPS) using the C/A-code.
6.02
PSEUDO-1 NOISE CODE A binary sequence (string of O's and lfs) Has noise-like properties Maximum autocorrelation at zero lag
BINARY BIPHASE MODULATIONS
Modulations on a constant frequency carrier Two phase states: normal and mirror image, modelled by y = A(t) cos( cot + <j> ), where A(t) is a sequence of values which are either +1 representing normal phase state (binary 0) -1 representing mirror image state (binary 1) Phase shifts restricted to 0 and 180 degrees Transition from binary 0 to 1 or 1 to 0 results in 180 degree phase shift Sequence of all 1's or all O's results in no phase shift SPREAD SPECTRUM Characteristic of a signal which is spread over a frequency band much wider than the minimum bandwidth needed to transmit the information being sent.
Canadian GPS Associates, May 1987
DEFINITIONS
Three unfamiliar terms, necessary to understand GPS codes, are defined here.
Pseudo-random noise code. A code is a system for representing information, together with rules for using the system. Many codes nowadays are binary codes (strings of O's and 1's), because they are handled easily and naturally by computers. GPS codes are binary codes. In information theory, a quantity is noise if it is not wanted, or does not contain useful information. Random noise has some specific statistical properties; in particular, the autocorrelation function (plot of the scalar product of the noise sequence with a delayed copy of itself) is zero, except at zero delay. Random noise is assumed to arise from unpredictable causes, and therefore would be subject to the rules needed to specify a code. Pseudo-random noise, on the other hand, is predictable and may contain information (namely, a reading of the transmitter clock). The pseudo-random noise codes used in GPS are generated by a carefully specified algorithm, and precisely identical codes are generated independently in both the satellite and receiver in order to track the signal. The term pseudo-random noise is used to indicate that the code has some random noise-like properties, in particular the autocorrelation is small (but not zero) except at zero delay.
Binary biphase modulations. This is the method used to implement the pseudo-random noise code in GPS signals. An unmodulated (pure) carrier signal carries no information and no code. Many kinds of modulations are possible: the frequency, amplitude, or phase of the carrier can be changed. Phase modulations are used for GPS. Since a binary code is used, there must be two states of the phase modulations. For GPS these two states are the normal state, representing a binary 0, and the mirror image state, representing a binary 1. The normal state can be modelled by multiplying the unmodulated carrier by +1 (no change). The mirror image state can be modelled by multiplying the unmodulated carrier by -1 (phase reversal to obtain the mirror image state). A code transition from 0 to 1 (normal to mirror image) or from 1 to 0 (mirror image to normal) each involves a phase reversal or 180 phase shift.
Spread spectrum. An unmodulated carrier has a bandwidth approaching zero (its spectrum looks like a single spike). The modulations needed to communicate information result in a broader bandwidth. If the modulations result in a bandwidth much broader than the minimum required to transmit the infomation, then the signal has a spread spectrum. The reasons for using a spread spectrum signal are to combat interference and ensure communication security. For positioning, however, spread spectrum signals also permit more accurate range measurements. The signal spreading method used by GPS, called direct sequence spreading, [Dixon, 1976] uses a pseudo-random noise code and binary biphase modulations.
CO
H
t=l
H ^
g.
Code chips. The P-code is at fQ, and the C/A-code at f Q /10. Here we illustrate the concept with code chips at fQ. The values are either 0 or 1, and are called chips rather than bits to indicate that they do not carry data.
Code states. A code chip value of 0 corresponds to a code state of +1 (normal state). A code chip value of 1 corresponds to a code state of -1 (mirror image state).
Modulated carrier. The carrier is modulated by multiplying it by the code states. Each transition from +1 to -1 or from -1 to +1 leads to a 180 phase shift of the carrier.
What happens when a GPS code is multiplied by a replica of itself? Since the code can be represented by a sequence of +1 and -1 values, when we multiply the two copies together, the states of the resulting signal will be the multiples of the individual states of each of the two copies. When the code and its replica are synchronized (lined up in time), the corresponding states of each copy will also be synchronized: that is, +1 states will occur simultaneously in both copies, as will -1 states. Then the states of the resulting multiplied signal will be either (+1)2 = +1 or (-1)2 = +1. In other words, only the normal state remains, and the code modulation has been removed from the signal. On the other hand, if the original code and its replica are not synchronized, then the modulation will not be removed from the resulting multiplied signal. Detecting the presence of modulations in this resulting signal permits adjustment of the time delay between the original code and its replica until the modulation is no longer present in the multiple and synchronization has been achieved. These features of the GPS code are crucial to the design of GPS receivers.
6.04
Polynomial = 1 + x + x
10
Output
;th
th 0 (i+1)
C/A-
Gl =
10
Output
Canadian GPS Associates, May 1987
CODE GENERATION
The pseudo-random binary sequences of the C/A-code and the P-code are generated using devices called 'tapped feedback shift registers' [Golumb, 1967]. The operation of these devices is described and applied to the generation of GPS codes.
A shift register is a set of connected one-bit storage cells (here labelled 1 to 10). When a clock pulse is applied to the register the content of each cell is shifted one bit to the right. Hence the content of cell 10 is 'read out' as output. The special properties of such shift registers depends on how information is 'read in' to cell 1.
A tapped feedback shift register is a shift register for which the input (the value to be placed in the leftmost cell) is determined by the state of one or more of the other cells in the register. Here we show that the binary sum of the states of cells 3 and 10 are used to determine the input to cell 1. Hence, if cells 3 and 10 have different states (one is 1 and the other 0), a 1 will be read into cell 1 on the next clock pulse. If cells 3 and 10 have the same state (both either 1 or 0), then a 0 will be read into cell 1. A shorthand way of denoting this particular design is by the polynomial 1 + x3 + x10. As illustrated, if we start with every cell containing a binary 1, 12 clock pulses later we will have the contents 0010001110. The next clock pulse will take the 1 in cell 3 and the 0 in cell 10, add them to get a binary 1, and place this in cell 1. Meanwhile, all other bits have shifted along one cell to the right, and the 0 in cell 10 becomes the next bit in the output sequence.
The C/A-code is generated by two 10-bit tapped feedback shift registers. One of these, Gl, is the same as the 1 + x3 + x10 register already described. The other shift register, G2, has a polynomial representation 1 + x2 + x3 + x6 + x8+ x9 + x10. That is, cells 2,3,6,8,9, and 10 are tapped and binary added to get the new input value for cell 1. In this case, the output comes not from cell 10 but from a second set of taps. Various pairs of these second taps are binary added, and the result binary added to the output of Gl. That becomes the C/A-code. The tapping shown here is for the C/A-code for the satellite with PRN 1. Various alternative pairs of taps are used to generate the complete set of 32 PRN C/A-codes. This family of codes is known as Gold codes [Gold, 1967], and has the property that any two have a very low cross correlation (are nearly orthogonal).
The P-code generation follows the same principles, except four 10-cell shift registers are used. Two registers are combined to produce the XI code, which is 15 345 000 chips long and repeats every 1.5 seconds; and two registers are combined to produce the X2 code, which is 15 345 037 chips long. The XI and X2 codes can be combined with 37 different delays on the X2 code to produce 37 different one-week segments of the P-code, five of which are associated with ground stations, and each of the remaining 32 with a different satellite.
One Master frame contains all the information in all pages of subframes 4 and 5, and consists of 25 complete frames. One Master frame takes 12.5 minutes to transmit.
Each message bit lasts for 20 milliseconds. During this period, the complete C/A-code (1023 chips at 1.023 Mbps) will have repeated exactly 20 times. A divide-by-20 counter on the C/A-code output provides the timing for the message bits. During one message bit 20 460 C/A-code chips, 204 600 P-code chips, and 31 508 400 cycles of the LI carrier will have occurred. Each message subframe lasts for six seconds. During this period, the XI component of the P-code will have repeated exactly four times.
6.06
SUBRAME 1 Flags (L2 code & data; week #; satellite accuracy and health) Age of data Satellite clock correction coefficients SUBFRAMES 2 & 3
Orbit parameters
SUBFRAME 4 Almanac for satellites 25-32 (pages 2,3,4,5,7,8,9,10) Ionospheric model, and UTC data (page 18) Antispoof flag - 32 satellites (page 25) Satellite configuration - 32 satellites (page 25) Health of satellites 25-32 (page 25)
Reserved (pages 1,6,11,12,16,19,20,21,22,23,24) Spares (pages 13,14,15) Special messages (page 17)
SUBFRAME 5
Almanac for satellites 1-24 (pages 1-24) Health of satellites 1-24 (page 25)
Subframe 4, so far, has sensible information on only 10 of its 25 pages. These pages contain an ionospheric model, UTC data, flags for each satellite indicating whether 'antispoofing' is on (the Y-code, or encrypted version of the P-code, replacing the P-code) and, if more than 24 satellites are in orbit, the almanac data and health for the satellites in excess of 24.
Subframe 5 contains the almanac data and health status for the first 24 satellites in orbit. The almanac data, a crude description of the satellite orbit, is used to determine where every satellite in the constellation is, and then to acquire signals from satellites which are above the observer's horizon but not yet being tracked. Hence, once one satellite is being tracked, acquisition of signals from all other satellites is relatively easy.
At the beginning of each 6-second subframe, are two special words called the telemetry word (TLM), and the hand-over word (HOW). The HOW word contains the so-called 'Z-count', and changes every 6 seconds. The TLM word changes when uploads or other satellite operations are underway.
The Z-count is the number of 1.5-second XI code repetitions since the beginning of the week (Saturday/Sunday midnight in the GPS system time scale). The Z-count thus ranges from 0 at the start of the week to 403 200 just before the end of the week. If we know the Z-count, we know exactly where in the week-long P-code to attempt to acquire and track the P-code.
The HOW word contains a number which, when multiplied by 4, gives the Z-count of the next 6-second subframe. It also contains the subframe identification (1 to 5), and a warning flag indicating when the subframes are not exactly aligned with the XI code (this is expected to be an extremely rare event).
The TLM word contains a fixed 8-bit synchronization pattern, and a 14-bit anc jw which contains, when appropriate [van Dierendonck et YTAUV/AI WV/AAMUJII.O, wiivn appivjpnaut/ ivaii A *'""*** f***f* ' message al., 1978], the status nf imlnaHino Hata frn th= catpllite the status of uploading data to the satellite diagnostic messages messages such as the Z-count of roll momentum dump.
M.
6.07
LI (154 F)
C/A-CODE (F/10)
MESSAGE
P-CODE (F)
L2 (120 F)
1
OSCILLATOR (F)
Mod 2 add
A B Mixer A + B Sum
All components of the signals leaving the satellites are coherently generated from the same fundamental frequency, f Q = 10.23 MHz (point A on the figure). Two carrier frequencies are used: LI with fj = 154 f 0 = 1575.42 MHz, and a wavelength of about 20 cm (point B), L2 with f 2 = 120 f 0 = 1227.60 MHz, and a wavelength of about 25 cm (point C). Three phase modulations are imposed on these carriers: The C/A-code C(t) has a chip rate of f Q 710 = 1.023 Mbps = f{ I 1540, a wavelength of about 300 m, a length of 1023 bits, and a period of 1 ms (point D). A different C/A-code is assigned to each GPS satellite. Only a 90 phase shifted version (quadrature component) of the LI carrier (point E) is modulated by the C/A-code. The P-code P(t) has a chip rate of f Q = 10.23 Mbps = fl I 154 = f 2 / 120, a wavelength of about 30 m, a length of 235 469 592 765 000 bits, and a period of 38 weeks (266 days), 9 hours, 45 minutes and approximately 55.5 seconds (point F). A different one-week segment of the P-code is assigned to each of up to 32 GPS satellites. Both the LI and L2 carriers are modulated by the P-code. The data sequence D(t) has a bit rate of f Q / 204 600 = 50 bps = ft / 31 508 400 = f 2 / 24 552 000 (point G). One frame of D(t) has a length of 1500 bits, and hence a period of 30 seconds. The complete signal leaving the satellite antenna (point H) can be represented by
AcC(t)D(t) sin(2;cf 1t+(|)c)+Ap P(t)D(t) cos(27if jt+4> J + Ap P(t)D(t) cos(27cf2t+(|> 2)
where the first component is the C/A-code signal (point J), the second component is the LI P-code signal (point K), and the third component is the L2 P-code signal (point L). The code and data phase modulation sequences can be modelled as sequences of 1 amplitude states. The data-and-code combinations C(t) D(t) (point M) and P(t) D(t) (point N) are modulo 2 additions (e.g., C(t) D(t) = +1 if C(t) = D(t), and C(t) D(t) = -1 if C(t) * D(t)). It is important to stress again that these modulations can be stripped from the carrier by multiplying them by a replica of themselves (or even by squaring them) in other words, a double phase reversal recovers the original signal.
^ "\
L TOL
Ll (154 F)
C/A-CODE (F/10)
MESSAGE
P-CODE (F)
the equipment group delay specification (Section 3.3.1.7) is worded: The delay between the L-band radiated output of a specific SV (measured at the antenna phase center) and the output of that SV's on-board
frequency source consists of a bias term and an uncertainty. The bias term is of no concern to the user since it is included in clock correction parameters relayed in the NAV data, and is therefore accounted for by the user computations of system time. The uncertainty (variation) of this delay as 'well as the difference between the LI andL2 delays are: 0 the effective uncertainty of the group delay shall not exceed 3 ns (2 a) 0 the group delay differential between the radiated LI and L2 signals consists of random plus bias (mean differential) components. For a given navigation pay load redundancy configuration, the absolute value of the mean differential delay shall not exceed 15 ns. The random variations about the mean shall not exceed 3 ns (2 a).
the signal coherence specification (Section 3.3.1.8) is worded: All transmitted signals for a particular SV shall be coherently derived from the same on-board frequency standard; all digital signals shall be clocked in coincidence with the PRN transitions for the P-signal and occur at the P-signal transition speed. On the LI channel the data transitions of the two modulating signals shall be such that the average time difference between the transitions does not exceed 10ns (2 a).
PECULATION;
6.09
AVAILABILITY
C/A-CODE (F/10)
OSCILLATOR (F)
Two options for what we can assume about how SA will be implemented are: We can assume SA will affect both C/A-code receivers and P-code receivers. In this case, the SA must follow some algorithm that P-code receivers can overcome by using some special device or information. We can assume that SA will affect only C/A-code receivers, and thus be implemented by adding some kind of *dither' to one or more parts of the C/A-code on LI Ac C(t) D(t) sin This figure shows where three kinds of dither might be inserted: 1 - The message data D*(t) transmitted together with the C/A-code is changed (dithered) relative to the 'true' message data D(t) transmitted together with the P-code. However, the message data is changed only once per hour. Hence any implementation of SA this way would have to be represented by 'constant' errors in some of the message parameters over that hour. 2 - The clocking of the combined data and code on the SPS signal C(t*) D(t*) is dithered in time relative to the 'true' clocking for the C(t) D(t) combined data and code. However, the P-code and C/A-code transitions are supposed to be simultaneous within 10 ns. This would place a limit on SA well below 100 m . 3 - The phase noise (|>c * on the SPS carrier is increased. However, the phase noise (and quadrature) of the P-code and C/A-code carriers is supposed to be less than 5 degrees of phase. This is an infinitesimal size compared to the phase noise which would be required to implement 100 m SA.
V * V
CHAPTER SEVEN
ANTENNAS AND
7.00
VERS Basic concepts Code and code-free channels Continuous, sequencing, and multiplexing receivers Single vs. dual frequency receivers Frequency standards and clocks Available receiver types Receiver cost trends
We conclude this chapter with a classification of some of the receivers currently available and look at the trend in their cost.
7.01
induced current
direction of propagation
current flow
The incident electromagnetic wave is characterized by its amplitude (how strong the signal is), its frequency (time dependence), its phase (referenced to how far the wave has travelled from its source), and its polarization (orientation in space). When the wave is incident on the surface of an electrical conductor, a current is induced to flow in the conductor which then sets up its own electromagnetic field. The boundary condition (a statement of physical laws which describes the interaction between electromagnetic energy and electrically conducting materials) requires that the total tangential electric field (incident field plus induced field) on the surface of an ideal electrical conductor must be zero. The induced field is caused by the induced current, hence there is a direct relationship between the incident electric field and the induced current. If the receiving antenna is connected to a circuit, this induced current will flow through it. Since the current is directly responsive to the incident electromagnetic field, any information contained in the field will also be present in the induced current, hence coded (modulated) information can be extracted from the current in the receiver.
The relationship between the induced current and the incident field is a complex function of frequency, amplitude, and phase distribution of the incident wave, as well as the shape, size, and orientation of the antenna element.
7.02
Monopole
Quadrifilar helix
Microstrip
Spiral helices
All GPS system antennas are required to be circularly polarized. In most cases, the receiver preamplifier is built as an integral part of the antenna baseplate in an attempt to minimize signal loss and noise introduction.
The monopole/dipole configurations usually require a groundplane for proper operation. They are single frequency antennas and have the advantage of small element size and simplicity of construction. The groundplane size can be used as a design parameter to reduce multipath interference.
The quadrifilar or volute antennas are single frequency devices which are more complicated in construction than the dipoles, are difficult to adjust (at time of manufacture) for proper element phasing and polarization, and are not azimuthally symmetric. They do, however, have good gain and pattern coverage characteristics. They do not require a groundplane, although one may be used for rear lobe control. Some designs incorporate a small pad of electrically lossy material for a baseplate so as to reduce the rear lobe radiation.
Microstrip antennas are the most rugged and simple in construction. They may be either single or dual frequency antennas with a very low profile which makes
them ideal for airborne applications. The main disadvantage of this class of
antennas is the relatively low gain, although this can be largely offset by low noise preamplifiers. Spiral helices offer the advantage of dual frequency operation at both LI and L2. The antenna has good gain and pattern coverage characteristics, although its high profile puts it at some disadvantage with respect to its phase characteristics. Like the volute design, the tapered helices exhibit some azimuthal asymmetry which is partially offset by careful ground orientation.
7.03
Multipath reflection
Earth horizon
LET PHAS
TANGENT SPHERE METHOD
7.04
Actual equiphase contour Apparent phase centre for each observer position
Antenna
Rotation axis
[ Apparent ] nt 'centre VJ
TRENDS IN ANTENNA
7.05
The term channel may also be applied to the transmission medium (the space between the satellite and the antenna) or the complete circuit comprising the satellite transmitter, transmission medium, and the path through the receiver that a particular signal follows.
From the point of view of the broadcast message, GPS is a spread-spectrum communications system [Ziemer and Peterson, 1985; Jordan, 1985]. A spreadspectrum system is one that takes a signal of small bandwidth yW, where y 1, and converts it (spreads it) into a signal of large bandwidth, W, by modulating the information waveform with a wideband waveform, such as pseudo-random noise. For GPS, the message bandwidth (null-to-null in power spectrum) is 100 Hz, and the P-code bandwidth, W, is about 20 MHz, hence, y is 5 x 10~6 or, for the C/Acode, W is about 2 MHz and yis 5 x 10"5. The wideband signal is transmitted through the channel, received, and then compressed back (despread) into the original small bandwidth signal by demodulating the signal with the identical wideband waveform. Spread spectrum systems have the advantage of being secure since the receiver must have access to the wideband waveform used at the transmitter. They also possess some immunity from jamming and, of course, are capable of high resolution ranging.
The signals transmitted by the satellites are detected and decoded by the hardware and software channels in a GPS receiver. These channels can be mechanised in different forms depending on the measurement type being sought and whether a message-decoding capability is required.
7.07
Bandpass filter passes a selected band of frequencies Frequency multiplier used to obtain a higher frequency harmonically related to oscillator frequency
A B
AB
Mixer multiplies two signals together to generate additional signals with frequencies equal to the sum and difference of the original frequencies
RECEIVER RF COMPONENTS
In order that a reader with a limited background in the workings of radio receivers may understand how the different types of GPS receiver channels work, we will briefly describe the operation of some common radio frequency (RF) components. Signal processing in a microwave radio receiver is generally more easily handled if the band of signals arriving at the antenna are translated to a common lower frequency band. This translation is achieved by combining the incoming signals centred on a frequency of f c with a pure harmonic signal of frequency f Q from an oscillator in a device known as a mixer. The mixer functions as a mathematical multiplier of the two input signals. The output of an ideal mixer is a superposition of two bands of signals: one centred on a frequency of f c .+ f Q and one centred on f c - fQ, called the up-conversion and down-conversion components, respectively. A bandpass filter is used to select the down-conversion component. All frequencies outside of the band centred on f c - f Q are rejected. The central frequency, f c - f Q , is known as an intermediate frequency (IF) or beat frequency. A GPS receiver in general will have a series of IF stages. In order to pass the P-code pseudo-random noise modulation, the bandwidth of a filter near the front end of a GPS receiver must be about 20 MHz. Once the code modulation has been removed by despreading, and as long as the Doppler shift of the received signals has been annulled by phase tracking the received carrier, the width of a filter in a GPS receiver can be of the order of the bandwidth of the message (about 100 Hz). Initial signal acquisition may require a bandwidth of several kiloHertz.
If the signal of interest is a slowly varying direct current (DC) signal or is composed of a band of frequencies centred on DC, then a lowpass filter may be used to remove all higher frequencies, such as those of extraneous noise received by the antenna or generated within the receiver. Frequency oscillators generate signals at a precise single frequency or, for a voltage-controlled oscillator, over a narrow range of frequencies. To obtain other frequencies that are harmonically-related to the oscillator frequency, a frequency multiplier (or divider) is used. The output of a frequency multiplier is a signal with a frequency which is an integer multiple of the input frequency. Similarly, a frequency divider yields the input frequency divided by an integer. By combining multipliers and dividers, a range of frequencies may be synthesized from a single input frequency.
oo
CORRELATION CHANNEL
Most of the receivers that have been developed are code-correlating receivers which are composed of one or more correlation channels. It is the only type of receiver that can decode the message. This block diagram of a correlation channel is intended to portray the concepts involved, and the detailed design of the channel will vary from receiver to receiver. Links to the receiver's microprocessor are not shown. A correlation receiver consists of two major blocks: a code tracking loop and a carrier tracking loop. The purpose of the code tracking loop is to provide a measurement of the pseudo-range either from the P- or C/A-code and at the same time to 'despread' the signal so that the message may be subsequently retrieved. In the loop, a code generator replicates the code that is imposed on the signal at the satellite. The start of the code is offset in time to account for the propagation delay between the satellite and the receiver. The correct delay is achieved by crosscorrelating the received signal with the receiver-generated code. If the two codes are not aligned, a correction signal is applied to the clock controlling the receiver code generator until a match is obtained. Once the receiver code is aligned with the incoming code, it stays locked to it. This type of tracking loop is called a delay lock loop. There are various ways of implementing a delay lock loop such as the 'full-time' and the 'tau-dither' early-late tracking loops [Ziemer and Peterson, 1985]. When the tracking loop is locked, a reading of the delayed code clock is essentially a reading of the satellite clock. The offset of the code clock from the receiver clock, which is assumed to be keeping GPS time, is the pseudo-range. Once the code tracking loop is locked, the PRN code can be removed from the satellite signal simply by mixing it with the locally generated one and bandpass filtering the resultant signal. The despread signal then passes to the carrier tracking loop which demodulates the satellite message by aligning the phase of the channel's local oscillator frequency with the phase of the intermediate or beat frequency. This action is achieved by controlling the frequency of a voltage-controlled oscillator. If the phase of the oscillator signal is not correct, this is detected in the demodulator and a correction signal is then applied to the oscillator. Once the oscillator is locked to the satellite signal it will continue to follow the variations in the phase of the carrier as the range to the satellite changes. The carrier beat phase observable is obtained in principle simply by counting the elapsed cycles and by measuring the fractional phase of the locked local oscillator signal. Most implementations of carrier tracking use the Costas loop method [Ziemer and Peterson, 1985]. Some use the squaring method. All carrier tracking loops suffer from a 180 phase ambiguity which must be resolved using the message structure. Temporarily-incorrect choices of phase can lead to so-called 'half cycle slips'. The message bits are decoded using bit synchronization and a data detection filter.
7.09
SQUARING
Phase inversions (code steps)
wwwww
Output
SQUARING CHANNEL
When the Global Positioning System was being designed, it was not widely realised that the satellite signals could be used without explicit knowledge of either the P- or C/A-codes. Prof. Charles Counselman, at the Massachusetts Institute of Technology, and Peter MacDoran, then at the Jet Propulsion Laboratory, independently developed code-free receiver designs. Both designs evolved into commercial products. In this section, we discuss the squaring channel used in Counselman's Macrometer receiver; in the next section, we discuss the code phase channel used in MacDoran's ISTAC receiver.
A squaring channel provides only the carrier phase observable. Neither pseudo-range measurements nor the satellite message is available. In principle, the operation of a squaring channel is very simple. In common with other GPS receivers, the incoming signal is first translated to an IF signal. The carrier, or rather the beat frequency carrier, is obtained simply by squaring this signal. Any phase inversions hi the IF signal due to the codes or message are thereby removed. This happens because a phase inversion can be considered as a change in the amplitude of the IF signal from +1 to -1 or from -1 to +1 and the instantaneous amplitude is therefore either plus or minus one. Squaring the signal results in a signal with a constant amplitude of unity, although with a frequency equal to twice the original. However, the phase of this signal is easily related to the phase of the original carrier. Of course, in the squaring process both the codes and the message are lost.
In practice, a channel in a squaring receiver would be more complex than that described here. Of particular concern, perhaps even more so than in a codecorrelating receiver, are the noise levels on the signals, since squaring the signals also squares the noise. This results in a signal-to-noise ratio loss which must be made as small as possible. Squaring can also be used in the carrier tracking loop of a code-correlating receiver to recover the beat carrier which is subsequently used to demodulate the message.
The ambiguity of the time delay measurements must be resolved before they can be utilized to determine the precise coordinates of the receiver's antenna. One technique that has been adopted to resolve the ambiguity uses the unambiguous Doppler shift of the received code frequency [Buennagel et al., 1984]. The frequency of the code can be measured using, say, a Fast Fourier Transform technique, from which the position of the receiver car be established to about 100 m. The C/A-code ambiguity can then be resolved. The P-code ambiguity can subsequently be resolved by interpolating within the C/A-code modulation.
T-H
r^
JS +
C3
&=!
E#3
&d >
Code-free
< P ^
#
0>
{==1
H ^ "*^ 1
{=11 H S ^^ ^5
T3
o
<
s J
w
r>
1)
o y
S -o 5 o > o 2
pO
^.& >>g -o *> o* a> ^ S a a 5 w a 3 2 | -S a-2 S a g & o 1 "* 8 OJ CO _fi vQ ^lx** S a> S > ^ O w to o * ^ 5f ^ o
* a-e,
S3 D S
<! *5
s
L*
' a>
V
a o
^*
o fe * 0 V ^ a, JS * cu
"s ^
M \
0 S2 z
U
^^
N4
^H P^
saovxuvAdv
saoviMVAdvsia
externally synchrony
S3
a,
N-l
sU *- S S a> C
a>
-*
S J3
SQUARING
7.12
CONTINUOUS
Four or more dedicated hardware channels Each channel tracks a single satellite signal Inter-channel biases must be calibrated May enjoy a signal-to-noise advantage Redundancy capability
SWITCHING
One or more hardware switching channels Each channel samples more than one signal Software more complex Potentially more difficult to keep track of phase In principle, cheaper than continuous receiver
HYBRID
A continuous-tracking receiver has four or more dedicated hardware channels. Each channel tracks a single satellite and maintains continuous code and/or phase lock on the signal. Each channel is controlled and sampled by the receiver's microprocessor with input/output operations being performed fast enough so that tracking is not disturbed. Continuous-tracking receivers may enjoy a signal-to-noise advantage over switching receivers in that the satellite signal is continuously available and may be more frequently sampled. A further advantage is a potential redundancy capability. Should one of the hardware channels fail, it may still be possible to obtain sufficient data to determine a position. One disadvantage of a multichannel receiver is that the differences in signal path delay in the channels, the so-called inter-channel biases, must be well calibrate^.
A switching receiver has one or more hardware channels. Each channel samples more than one satellite signal. Code and/or carrier tracking for the individual signals is performed in software (or, more usually, firmware) by the microprocessor. As a result, greater demands are placed on the microprocessor in a switching receiver and its programming is necessarily more complex; hardware complexity is exchanged for software complexity. The term pseudo-channel is sometimes used to identify software tracking channels. Depending on the application and the environment, a switching receiver may be more susceptible to cycle slips than a continuous-tracking receiver. If the cost of hardware components is a significant factor in determining the selling price of a receiver, then in principle, the cost of a switching receiver should be cheaper than a continuous receiver.
It is possible to combine continuous and switching channels in the same receiver. Such a receiver is termed a hybrid.
7.13
1
2
D D D D D D D D D D D D
D D D D D D 0 D D D D D D D D D D D D D D 0 0 D
D, D D D D D D D D D D 0
signals
signals
t
Canadian GPS Associates, April 1986
A multiplexing channel is one for which the sequencing time to sample all satellites assigned to the channel is equal to 20 milliseconds, the period of one bit in the satellite message. The sampling can be arranged so that no message bit boundary is spanned by any tracking interval. In this way, the messages from all satellites phase tracked by the channel can be read simultaneously. A multiplexing channel can be used to obtain both LI and L2 data by alternating between the frequencies every 20 ms. The TI4100 has a single multiplexing hardware channel that is used to sample both LI and L2 signals from up to four satellites.
If a channel switches between signals at a rate which is asynchronous with the message bit rate, the channel is referred to as a sequencing channel. A fast sequencing channel is one which takes the order of seconds to sequence through the signals. A slow sequencing channel takes minutes to hours. A single sequencing channel would lose bits in a particular satellite message during those intervals spent sampling the signals from other satellites. Consequently, sequencing receivers may have an extra hardware channel just for data reading. Alternatively, sufficient orbital data must be accumulated before the receiver starts the tracking cycle for real-time positioning. A disadvantage of a sequencing channel used to track both LI and L2 signals is that since the measurements on the two frequencies are not simultaneous, the ionospheric correction may be degraded.
a o P* P* a *i M CQ
OJD C
a> P^
^ a
P-J
1 % ^
-tt J5
CM O
O-e a
IBS
T3 g g O W W
^T 0)
5 s gj
S a
O w
-"
U Z
p ^
o'd HZ MO
t=4
W < A V -. fi
s-S 3 ^
L- .,
C3
o 5 03
a>
a C3
11
i c^
va "^
1-s
saoviNVAav
saoviMVAavsia
7.15
CLOC:
QUARTZ CRYSTAL frequency generated by vibrating crystal with applied electric field (piezoelectric effect) inexpensive, compact, and robust RUBIDIUM VAPOUR CELL based on 6.834 682 605 GHz resonance frequency of a hyperfine energy transition of rubidium 87 atoms least expensive of the atomic clocks, ^$5,000 CESIUM BEAM TUBE based on 9.192 631 770 GHz resonance frequency of a hyperfine energy transition of cesium 133 atoms resonance frequency used for the definition of the SI (atomic) second most accurate (frequency most precisely set) of the atomic clocks HYDROGEN MASER based on 1.420 405 757 68 GHz resonance frequency of a hyperfine energy transition of atomic hydrogen most stable of the atomic clocks
GUI
o O
oo
fe
m
H t=l
H
cs
UIUI
o O 'O
Frequency domain stability is usually expressed by the one-sided spectral density of frequency fluctuations, S (f), or the one-sided spectral density of phase fluctuations, S(f), where f is the Fourier frequency offset from VQ. The two quantities are related by Sy(f) = (f2/vQ2)S^(f). These spectral densities are invariant when, with ideal equipment, the frequency v(t) is multiplied or divided.
A more intuitive measure of frequency stability may be the time-domain two-sample (Allan) variance of the fractional frequency variations. If we define the fractional average frequency over an interval T at time t as
then the Allan variance is given by where the angular brackets denote an infinite time average. The Allan variances of typical quartz crystal, rubidium vapour cell, cesium beam tube, and hydrogen maser oscillators are shown in this figure [Hellwig, 1979]. The typical drifts of 1 part in 1010 and 1 part in 1012 per day for the quartz crystal and rubidium oscillators respectively, have been removed.
H t=l
HH
Good oscillator stability is particularly important in kinematic applications if the receiver clock is not updated at each solution epoch. Such will be the case during outages in the 18-satellite constellation. A precision oscillator will enable navigation during a brief outage with minimal degradation in position accuracy. An atomic frequency standard, such as a rubidium oscillator, may be advantageous in these situations. A precision oscillator will also extend the useful period of navigation for the prototype GPS constellation.
7.18
CODE
CORRELATION
CODE PHASE
CODE
CODE
Macrometer V-1000 E
EDO SatTrak
Norstar 1000
Canadian GPS Associates, April 1986
The only code phase receiver currently available on the market is the ISTAC-SERIES GPS Positioner Model 2002. It operates on the LI signals only, however, there are plans to extend its capability to two frequencies. Single and dual-frequency versions of the Macrometer squaring channel receivers are available. The pioneering single channel version is known as the Macrometer V-1000; the dual-frequency model is known as the Macrometer II. The selection of a specific receiver type or model should be predicated on its intended use. The physical environment, accuracy and reliability required, and the dynamics of the moving platform (if used) are parameters of major importance in the selection process.
7.19
U.S.$1,000,000
GPS Operational
$100,000
$10,000
$1000
$500
$100
1975 1980
I I 1985 1990 1995 2000
CHAPTER EIGHT
OBSERVATION EQUATIONS
8.00
OBSERVATION TYPES
Between-epoch single differences Between-receiver single differences Between-satellite single differences Receiver-time double differences Receiver-satellite double differences Receiver-satellite-time triple differences
OBSERVATION EQUATIONS
In this chapter we first describe the two main GPS observation equations: the equation modelling code pseudo-ranges, and the equation modelling carrier beat phase measurements. We will see that GPS observations can always be thought of as some kind of biased range. Both kinds of observation equations contain linear bias terms modelling the effects of receiver clock offsets, satellite clock offsets, ionospheric and tropospheric delays. The carrier beat phase observation equation contains an additional linear bias term the initial cycle ambiguity term.
We then describe various linear combinations (single, double, and triple differences) of these observation equations. These linear combinations are often used in processing GPS data. This is because various linear combinations reduce or eliminate the effects of some of the linear biases in the observation equations: those due to receiver and satellite clock offsets, and initial cycle ambiguities. Differenced observations also tend to reduce the effects of orbital errors, and the effects of ionospheric and tropospheric delays. However, this reduction is more effective the closer the separation of the simultaneously observing receivers.
In developing the observation equations, where necessary we shall use subscripts to indicate the receiver station and superscripts to indicate the transmitting satellite. These indices will be dropped when there is no risk of confusing the meaning of each quantity.
8.01
-ERVATION,
Satellite clock
GPS time
Receiver clock
P=
PSEUDO-RANGE OBSERVATIONS
The basic pseudo-range observable, dx, is the difference between the time of transmission (in the satellite time scale t) and the time of arrival (in the receiver time scale T) of a particular signal transmitted by the satellite. In addition to the satellite and receiver time scales, there is a more or less ideal time scale called GPS time T. For the pseudo-range, we can write:
dT = T(Tb)-t(Ta) .
- T we have
-(*.-
-t(ta)]-[VT(tb)]
Ideally, the first term is the travel time of the signal which, when multiplied by c = 299 792 458 m/secthe speed of light in a vacuumis equal to the true range p to the satellite, ignoring atmospheric effects. The second term represents the offset dt of the satellite clock from GPS time. The third term represents the offset dT of the receiver clock from GPS time. When ionospheric and tropospheric delays are taken into account, the complete pseudo-range equation takes the form
p = c-dT = p + c-(dt - dT) + dion +
We are primarily interested in receiver coordinates. Where are they in this equation? In fact, they are hidden in the true range p. If the satellite coordinates are perfectly known, then p = ||r - R|| . However, r is generally not known perfectly, so we add a term, dp, to account for the effect of ephemeris errors and get p = ||r - R|| + dp . Then the observation equation becomes
p = ||r-R|| + dp+c-(dt-dT) +
which is in the form / = f(R,r,b), where b is a vector of terms we may wish to consider as biases (any or all of dp, c-(dt - dT), dion, and d^. ). Ways of dealing with these bias terms are explained in Chapter 9. Finally, we must consider effects due to measurement noise and unmodelled influences, both of which we assume to be random. These are included in the model by adding a residual term to the right-hand side of the observation equation.
EAT PHASE
8.02
GPS Oscillator
where <|> is the phase of the signal transmitted by the kth satellite at time t, and fy is the phase at the ith receiver at reception time T (both expressed in cycles).
The carrier beat phase model can be reformulated using two relations The phase/frequency relationship which is valid for highly stable oscillators over a short time interval is
and obtain
Then
T-t = dt-dT and the new carrier beat phase model is $ = - (f/c).p - f-(dt - dT) - (f/c).(-dion +
8.03
PHASE
-/ \->al I ivl t. l will
.AAAAAAAA
Receiver lock-on
-X .^VtV^lVtVI
Received Carrier
VxdA 1 Id
AAAAAAAA-^t =t
fraction of I cycle
u >% U CM o
JD
E a
could also be called the cycle unknown. As long as the receiver maintains continuous phase lock during an observing session (that is, as long as the Int(<|); tQ, t) count is maintained), there is only one ambiguity per satellite/receiver pair. However, if there are breaks in phase lock, several ambiguities or cycle slips per
satellite may result. Ways of dealing with this problem will be discussed in Chapter 9. What is then observed by the receiver? It is <l>measured = Fr(<W + Int(4>; to* *)
Hence <t>total = ^mQ^ re<4 + N(tQ). The equation for the instantaneous carrier beat phase observable for one satellite, one receiver, and one instant of time may be written as
(total = - <f/c>P -f'(dt - dT> - (f/cH-dion + <W + N '
where the phase difference between the satellite and receiver oscillators are now implied by the term f-(dt - dT) for the satellite and receiver clock errors.
Multiplying by the wavelength A, = c/f, and defining
~ " '^measured'
one can obtain the corresponding carrier phase equation (in length units) <D = p + c-(dt - dT) -.- X.N - dion + d^, which is directly comparable with the pseudo-range equation, except for the ambiguity term (and the sign of the ionospheric term). This clearly indicates that except for the initial ambiguity, N, O can be thought of as a range, and its observation equation is almost the same as that of the biased range. The biases for this kind of measurement are the same as for pseudoranges, except that the ambiguity term has been added.
8.04
BETWEEN EPOCHS
BETWEEN RECEIVERS
O O
Using GPS data, both point and relative positioning are possible. However, because of the uncertainties in the satellite positions, clock behaviour, and propagation delays, 'absolute* point positioning will likely be obtained with accuracies of only a few metres. For most positioning requirements for geodetic and geodynamic applications, it will be necessary to use GPS in a relative mode. Of the error sources affecting the GPS signals, satellite orbit, satellite clock, receiver clock, and atmospheric propagation errors will exhibit some correlation among signals received at several stations simultaneously tracking the same several satellites. The goal of the differential process is to take advantage of these correlations to improve the accuracy of relative positions. Using the code or carrier phase measurements for relative positioning consists of taking differences between measurements. In this way, the effect of the various errors, common to the measurements being differenced, are removed or greatly reduced.
GPS measurements can be differenced between receivers, between satellites, and between epochs, or combinations thereof. Many different differencing combinations are possible. Presently, the most commonly used approach is to difference in the above order, i.e., first between receivers, then between satellites, and lastly between epochs. The notation we use for taking these differences is intended to be mnemonic: A denotes differences between two receivers (two vertices on the bottom) V denotes differences between two satellites (two vertices on the top) 8 denotes differences between two epochs.
8.05
Static receiver
Kinematic receiver
we obtain
8O = Sp + c-(Sdt-SdT) -8d-ion
Note that this observation is free of any cycle ambiguity.
Such observations can be expressed in the form / = f(R, r, b), where b is a bias vector containing whatever terms we wish to consider as biases; R and r are contained in 8p = ||r - R||epoch 2 - II1" - RHepoch i and / = So ^ observed phase change.
The processing technique required to correctly interpret the data is determined by the relationship among the epochs over which the change of phase is measured:
If the Doppler count is recorded for a short interval, the counter reset to zero, and sometime later, repeatedly, another Doppler count recorded, we obtain Intermittently Integrated Doppler (IID) measurements. Using this technique, much potentially useful data is ignored. A second approach is to reset the counter to zero after each counting (integration) interval, but immediately start counting again, without losing any cycles. This type of data are Consecutive Doppler Counts (CDC). These measurements can be interpreted as range changes over the integration intervals. The third approach is to simply sum the consecutive Doppler counts throughout the satellite pass to obtain Continuously Integrated Doppler (CID) measurements. These measurements can be interpreted as the range changes since some initial (signal lock-on) epoch. If the initial range is considered as a range unknown (ambiguity), the Doppler observation can be considered to be a difference of two biased range observations.
8.06
/ ~ v
'receiver 2
''receiver 1
=(> receiver 2
to obtain
" ''receiver 1
= Ap - c-AdT + Ad ion
tropospheric delays, and AdT is a differential correction for the receiver clock errors.
Similarly, if the two stations are collecting carrier beat phase data, a between-receiver single difference carrier phase observable can be derived as l AO = Ap - c-AdT + X-AN - Adiorl + Adtrop
Between-receiver single differences remove or greatly reduce the effects of errors associated with the satellites: satellite clock errors, and to a great extent orbit errors and atmospheric delays as well (if the baseline lengths are short compared with the 20 000 km altitude of the GPS satellites).
WEEN-SATELLITE LE
8.07
V(
satellite 2
satellite 1
VP
is given by Vp = Vp
The corresponding model for the between-satellite single difference carrier phase observable is
VO = Vp + c-Vdt + A,.VN-Vd ion + The between-satellite single difference observable is free from receiver clock errors. If the goal of a particular GPS observation campaign is eventually to form double differences, the between satellite differencing could be done at each individual receiver at the time of observation.
iLE
:EIVER=TIME
8.08
ATELLITE
8.09
VAd^ + VAd^
VA<D = VA p + A,-VAN - VA dion + VAdtrop
& Canadian GPS Associates, May 1987
Consecutive Doppler
DOPPLER
Correlation is Removed
GENERAL
trop
any combination of the differencing operations we have considered (A, V, or 8) may be represented by an operator (matrix) D, which is composed entirely of the elements +!,-!, and 0 operating on a vector / of 'observed ranges' p or O.
If the covariance matrix of the vector of ranges / is Cj, then the covariance matrix of the differenced range, D /, is
This provides a straightforward (but perhaps computationally inconvenient) method of accounting for the mathematical correlations introduced by differencing.
Physical correlations are more difficult to handle. Physical correlation may be spatial or temporal in nature, or a combination of both. Propagation anomalies (atmospheric inhomogeneities) are one significant source of both spatial and temporal physical correlations. To completely account for the correlations from just this one source would require a realistic model for the atmospheric parameters.
Proper accounting for the correlated nature of differenced data is necessary in order to make the differenced model equivalent to that involving the uncorrelated undifferenced data. However, correlation matrices are often computationally inconvenient. Hence, of necessity, differenced observables are often treated as uncorrelated. Excellent results have been obtained even when correlations have been neglected.
CHAPTER NINE
9.00
SES SATELLITE DEPENDENT Orbit representation biases Satellite clock model biases
STATION DEPENDENT Receiver clock biases Station coordinates
OBSERVATION DEPENDENT Ionospheric delay Tropospheric delay Carrier beat phase ambiguity
Residual Biases Cycle Slips Multipath Antenna Phase Centre Movement Random Observation Error
Canadian GPS Associates, May 1987
(a) Satellite biases consist of biases in the satellite ephemeris (e.g., the satellite is not where the GPS broadcast data message or other orbital information tell us it is), and biases in models for the satellite clocks supplied in the broadcast message (e.g., the satellite clocks, even with the broadcast message models, are not perfectly synchronized to GPS time). These biases are thought to be uncorrelated between satellites. They affect both code and carrier beat phase measurements equally, and they depend on the number and the location of the tracking stations providing data for orbital determination, the orbital force model used, and the satellite geometry. (b) Station biases usually consist of receiver clock biases and, for non-positioning types of GPS applications, such as time transfer and orbital tracking, of biases induced by uncertainties in the coordinates of the stations.
(c) Observation dependent biases include those associated with the signal propagation and other biases dependent on the observation type, such as, for instance, ambiguity biases inherent in the carrier beat phase observables.
The effect of biases is removed, or at least suppressed, by an attempt to model them. They are assumed to have functional relations with a variety of arguments such as time, position, temperature, etc.
Beside biases, the accuracy of positions and/or time obtained by GPS is dependent on two general influences: the geometric strength of the satellite configuration being observed, and the errors affecting the measurements themselves plus the remnant from the biases after the main effects have been modelled out. The former have already been discussed in Chapter 5, where it,was shown that the measurement errors propagate into the position proportionally to the various Dilution of Precision factors. Errors from each of the sources will have complicated spectral characteristics and other properties, and there will be correlations between some of these errors. However, at this stage in GPS development, the error models are usually limited to the simple approach of predicting typical standard deviations of uncorrelated equivalent range errors from each error source.
9.01
ALL TYPES OF BIASES ARE FELT AS RANGE MEASUREMENT ERRORS UNLESS MODELLED OR ELIMINATED
The range bias can be quite sizeable, depending on what hardware and what information is used to arrive at the observed biased range. Individual contributions to the range bias may reach the following maximum values:
orbital bias
300 000 m (decreases to 10 m if broadcast correction is used) 10 m to 100 m (depends on the type of receiver oscillator) 80 m for today's broadcast ephemerides (will decrease to 10 m when operational tracking network is in place) 150 m at the horizon (decreases to 50 m at the zenith) 20 m at 10 above horizon (decreases to 2 m at the zenith) anything.
Even though the individual contributions to the range bias would generally have varying signs (so that their sum would be smaller than the sum of their absolute values), it would clearly be impossible to work with these biased ranges. Thus an attempt must be made to remove the above biases by modelling or by differencing.
9.02
TIMING BIASES
SIZE OF BIASES
Satellite clocks 1 ms => 300 km without broadcast message corrections 30 ns => 10 m using broadcast message corrections
Receiver clocks
Synchronization with each other and with GPS time Short term stability (noise) Long term stability (drift and higher order terms)
OPTIONS FOR MODELLING BIASES Assume negligible and ignore Apply broadcast message correction Assume independent bias at each epoch for each clock Between satellite differencing 0 eliminates receiver clock bias Between receiver differencing 0 eliminates satellite clock bias Estimating independent clock biases at each epoch 0 same result as differencing Assume correlation between biases at different epochs Use time polynomial Treat as weighted parameters
Canadian GPS Associates, May 1987
where tQ is some reference epoch, and aQ is the satellite clock time offset, a is the
9.03
IT BIASEi
TTFTpl TT^v TT A (fa ~\ii i
Radial
Along track
Out of plane
SIZE
ORBIT BIASES
Satellite epherneris errors are the most difficult to deal with. Clocks would be improved, or their undesirable effects be eliminated or greatly removed, for instance, by simultaneous tracking by two or more receivers and the use of differenced observables. Ephemeris errors, on the other -hand, would require a better estimation of the orbits, a process which is hampered by insufficient knowledge of the forces acting on the satellites since often these forces cannot be measured directly and adequately from the ground tracking sites. At present, using the orbital data provided by the broadcast ephemerides, satellite positions can be computed with a typical accuracy of about 20 m, with occasional errors reaching up to 80 m. When the future Operational Tracking Network is in place, it is anticipated that the quality of the broadcast ephemerides will be improved to the point that satellite positions will most likely be accurate to about 5 m, with occasional errors up to 10 m. The main approach to solving the problems created by the imperfect modelling of the physical phenomena responsible for these forces has been usually directed toward their source; that is, toward developing more accurate models to begin with. Alternatively, it may be possible to use parametric models of various complexities for the ephemeris biases (i.e., with parameters that can be adjusted as part of an orbit estimation process).
9.04
Assume the given ephemeris is perfect, and hold the orbit fixed.
OPTION 2 ~ RELAX THE ORBIT
1) Assume geometrical biases, and estimate one to six bias parameters per arc
2) Assume dynamic model, use a force model and six initial conditions 3) Assume "free orbit", and estimate independent orbital biases at every epoch
OPTION 3 - DIFFERENCE THE DATA (RATHER THAN MODELLING THE ORBI
For example, between-receiver differencing reduces the effect of 20 m orbital error on baseline vector solution to 1 ppm or less.
Canadian GPS Associates, May 1987
An alternative approach is to assume that the satellite ephemeris errors manifest themselves as geometrical biases, such as translation, rotation, etc., of the given orbital arc with respect to the true orbit. In this case, one to six orbital bias parameters per arc are included in the model depending on the length of the arc. Usually three biases expressed in a reference frame attached to the satellite with the axes pointing in the along-, across-, and out-of-plane directions are used for short arcs.
A physically more meaningful approach is the inclusion in the model of initial conditions (e.g., six Keplerian elements) of the satellite orbit based on an assumed force model. The third alternative is to assume in the model independent orbital biases at every observation epoch. This is the free orbit approach.
A priori weights to the orbit bias model parameters are usually applied in all three of these approaches.
In relative positioning, the effects of certain systematic orbital errors can be removed or greatly reduced by differencing the observations from two or more sites observing simultaneously the same satellites. This is a common approach often used to date as an alternative to rigorous orbital modelling. As has already been indicated, in this approach, a 20 m orbital error, for instance, results in about 1 ppm error in a baseline vector solution. Any estimation of orbital biases using difference observations requires-that some information on orbital biases remains in the differenced observations. This means that GPS networks extending over large areas are more suitable than ones extending over smaller areas.
9.05
1PHE
l >~s
03
Dispersion curve
CO. CO
Ionospheric fluctuations
Note: co = 2 n f
Canadian GPS Associates, May 1987
terms of v as v = v + P 3 v / 9p. Corresponding to the phase velocity is the phase index of refraction n = v / c, and corresponding to the group velocity is the group index of refraction n = v / c. We can express n in terms of n^ as
The form of IL for the ionosphere is n^ 1 - a N / f2, where N is the ionospheric electron density in electrons/m3 at the time and place we are evaluating ru and a is a constant (see, e.g., Bomford [1980]). Hence n 1 + a N / f2.
o
The change in path length due to the ionosphere is Ap = J(n-l) ds, where the integration is along the propagation path. GPS code signals are dependent upon the group index of refraction, so the ionospheric group delay is Ap = J( n g ~l) ds = + a NT / f2. Phase measurements, on the other hand, are governed by the phase index of refraction, so the ionospheric phase delay is Ap = J(n -1) ds = - a NT / f 2 . Here NT is the total electron content along the propagation path in electrons/m2. NT varies from 1016 to 1018 electrons/m2. If Ap is in metres and f in Hz, then a 40.28. If we let dion = a NT / f2, the effect on observed code ranges is +dion, and the effect on observed carrier phase ranges is -dion. The magnitude of dion for GPS varies from tens of centimetres to tens of metres, depending on NT.
9.06
SIZE
OF BIASES Extremal value at zenith ^ 50 m Extremal value at horizon ~ 3 x 50 m Diurnal variation: day ^ five times night Annual variation: November ~ four times July Sunspot variation: at max ~ four times at min
FOR MODELLING BIAS Dual frequency data - use dispersion dion (LI) = [p (LI) -p (L2)] {f^ /(fj - { )}
single frequency between-receiver differences 0 reduces common ionospheric effect 0 residual effect: regional baselines - typically 1 ppm continental baselines - not as good
Canadian GPS Associates, April 1986
IONOSPHERIC MODELLING
At GPS frequencies, the ionospheric range effect may vary from more than 150 m (at sunspot maximum periods, midday, satellite near horizon) to less than 5 m (sunspot minimum, night, satellite at zenith). Since the ionospheric effect is frequency dependent, LI and L2 measurements can be compared to estimate the effect as
dp;on(Ll) = [p(Ll) - p(L2)][f22/(f22 . fj2)] .
A similar approach is used to correct the observed carrier beat phase measurements
by f 1/f2 N(L2)]
where N(L1), N(L2) are cycle ambiguities on the LI and L2 data. As long as the observed phases of die LI and L2 signals are continuously counted (i.e., there are no cycle slips), the last term on the right-hand side of the equation remains constant (but unknown). Hence a time series of the linear combination of the LI and L2 phases reflects mainly the variations of the ionospheric delay of the LI signals.
The dual frequency correction removes most of the ionospheric effect in the code or carrier phase measurements. The residual effect, however, may still be significant for some applications, in particular for mid-afternoon observations, and more so during a maximum solar activity cycle.
The above technique obviously can be only applied to two-frequency data. Users of single-frequency receivers must use ionospheric models or rely on the ionospheric model parameters transmitted by the satellites (cf. Chapter 6). However, these types of models are still in the early stages of development and are unlikely to be more than 75% effective so that for an ionospheric effect of 50 m, for instance, residual errors of up to 15m still remain. In fact, a study using actual data [Lachapelle and Wade, 1982] has shown that using the current broadcast ionospheric correction model actually gave worse results than not applying the correction at all.
9.07
20 metres
Ignore Use standard vertical refractivity profile models (a) with standard atmosphere (b) driven by surface weather data (c) with scale parameter estimated (d) both (b) and (c) Use local vertical refractivity profile model driven by surface weather data Water vapour radiometer observations for wet component Ray tracing (not yet attempted)
TROPOSPHERIC EFFECTS
Refraction in the neutral atmospherewhich includes the troposphere and other regions up to 80 km in altitudeis essentially independent of frequency over the entire radio spectrum. Unlike the ionosphere, the troposphere is not dispersive for frequencies below 30 GHz, so that group and phase delay are the same. Refraction in the neutral atmosphere may be conveniently separated into 'dry' and 'wet' components. The dry component can be approximated by
DTC = 2.27 10'3 P0 where DTC is the dry term range contribution in the zenith direction, and PQ is the surface pressure (in mbar). For an average atmospheric pressure of 1013 mbar, this corresponds to a range error of 2.3 m. The dry term contains approximately 90% of the total zenith range error, and may be estimated from surface pressure data with an accuracy of about 0.2%.
The wet component, on the other hand, depends on the atmospheric conditions all along the signal path. These conditions are not necessarily well correlated with surface conditions. There are various models available to date to model for this component taking into account such factors as the water vapour content, temperature, altitude, and the elevation angle of the signal path.
To date, the available tropospheric models, even with real-time meteorological data, appear to reduce the combined tropospheric effects by only 92% to 95% depending on the amount of atmospheric information available to the user. Users with good knowledge of the atmospheric conditions in the vicinity of their operations will have a smaller residual error than users who must rely on approximate or average data.
For differential observations, the differences of the tropospheric effects at the two ends of a baseline are also important. The available tropospheric models can be used to calculate this relative effect given the meteorological conditions at each site. In this case, one may choose the simplest approach, assuming that any tropospheric effects can be removed or greatly reduced by just differencing the observations from the two sites. This indeed may be a viable approach for most favourable cases, especially when the baselines are short. To achieve centimetre level accuracy, however, water vapour radiometers probably will need be used at both ends of the baselines to infer the integrated water vapour content along the lines of sight to the satellites. This is because the atmospheric water vapour content is likely to be quite variable in time and inhomogeneous spatially in many locations. This time variability is comparable to the time scale of the measurements at a given site. Typically, the wet component will be about 30 cm or more for observations at 20 elevation. If only surface measurements are used to estimate the total tropospheric correction, the uncertainty typically will be about 3.5 cm.
BEAT PHASE
9.08
An integer multiple of the carrier wavelength ~ 20 cm for code correlation ~ 10 cm for squaring As long as receiver maintains phaselock (absence of cycle slips) same value for all epochs different value for each satellite/receiver pair
OPTIONS FOR MODELLING BIAS
Estimate as a real number leave as real number choose nearest integer, and hold fixed Use combination of carrier and code measurements to estimate Search likely integer values for "best" values Difference between epochs to eliminate bias
Canadian GPS Associates, May 1987
The problem is further compounded by the occurrence of frequent cycle slips. While cycle slips have to be detected and rectified before the final parameter estimation, cycle ambiguity is normally estimated together with the rest of the unknown parameters.
9.09
RESOLUTION
"3-2
* ^^
0
Vf
*5
-5
AMBIGUITY RESOLUTION
One technique, usable with receivers, such as the TI-4100, which reconstruct the carrier with the knowledge of the pseudo-random noise (PRN) modulation code, relies on the synergistic combination of the pseudo-range and the carrrier beat phase measurements to a satellite to resolve the cycle ambiguities. It can be shown (e.g., Delikaraoglou et al. [1986]) that for each observation epoch where P-code pseudo-range and carrier phase measurements are available on the LI and L2 frequencies, this combination allows to estimate through a recursive algorithm the cycle ambiguities on LI and L2 as
N(L1) = 4.091 f x tj - 3.967 f 2 ^ - <|>(L1) + 5N(L1) N(L2) = 3.967 f x tj - 4.091 f 2 \^ - <|>(L2) + 8N(L2) with similar expressions for 8N(L1) and 8N(L2). In practice, the success of an ambiguity resolution of the desired subcycle level hinges on the capability of reducing the effect of the instrumental errors, especially those associated with multipath errors in the code measurements. Usually such limitations can be overcome by long-term averaging leading to average ambiguity estimates with formal errors of the order of 1 to 3 cycles.
The code/phase combination over the long-term averaging process can be conveniently complemented and/or checked against alternative algorithms which search for die integer cycle ambiguities on purely geometrical principles. For instance, it can be shown (e.g., Delikaraoglou et al. [1986]) that if one assumes some approximate values for the sought integer cycle ambiguities, integer corrections to them can be found as solutions to a line equation, i.e., (f2 / fx)5N(Ll) - 8N(L2) = x ,
in the two-dimensional real space coordinated by 8N(L1) and 8N(L2), whose slope is exactly known (i.e., equal to f 2 / fj = 120/154), whereas its intercept x is affected by instrumental and other errors. The best solution (because of the errors, the equation does not generally have an exact solution when SN(L1) and 8N(L2) are both integers) can be found through a search on an integer grid, (shown in the figure) for the point on this grid that lies closest to the straight line defined above.
9.10
CYCLE SLIPS
-500-r
0 -600 -i-
X
10 20 30 Elapsed time (minutes)
-700 -0
CYCLE SLIPS
When a satellite signal is obstructed in any way, it can no longer be tracked. When signal lock is resumed, the fractional part of the measured phase would still be the same as if tracking had been maintained. The integer number of cycles, however, exhibits a discontinuity or cycle slip. Such occurrences can be quite frequent, as shown in the figure.
There are various possible approaches to dealing with this problem. A common approach is to hold the station positions fixed (in a preadjustment), and edit the data manually by inspecting the corresponding residual series for breaks. This has been proven to work, but it can be a tedious and time consuming process. Another approach is to model the data with a piecewise continuous polynomial for each satellite. To implement this approach, however, the data have to be examined to find the times of slip occurrences (breaks) first. This again is quite tedious and, in addition, still some manual editing at the few cycles level may be necessary.
An alternative approach is to use between-receivers, satellites, and time triple difference observables to determine station locations. Once reasonably good values of station coordinates have been obtained, a search can be carried out automatically on the triple difference residuals to identify discontinuities in the double differences which form the specific triple differences [Remondi, 1985d]. Triple differences are particularly convenient to use in this connection, since they are not affected by receiver or satellite clock errors. Once such large residuals have been isolated over the whole epoch, one can easily determine which satellite is responsible for the integer cycle slips by looking at the triple difference residuals of the different satellite pair combinations.
MULTIFATH I
9.11
mmmimmsm
Coupling between antenna and image
Canadian GPS Associates, April 1986
9.12
MULTIPATH SIGNALS
Direct signal
n_TL_n
Interfering multipath signal Rejected multipath signal
n__TLTL_n
Time
Canadian GPS Associates, April 1986
9.13
DUE MOVEMENT
Phase centre
where the first exponential contains the time dependence and the second exponential contains the phase information due to the geometry of the source antenna, the excitation current, and the placement of the antenna and observer in a reference coordinate system. Distance measurements may be performed either by measuring the transmission time of the signal or by doing a signal phase measurement. In a phase reference system, a single phase measurement is made and the path distance is calculated from
R =(O-
n = (O - cc)v/co
where O is the indicated phase length of the path, and R is to be determined. Of course, R is a function of the geometry of the antenna and observer positions, i.e.,
R = z cosG - d sinO sin<|) - r , where the true unknown distance (i.e., from the coordinate origin to the observation point) is given by r. Thus
(O - CC)V/CG = z cos0 - d sin6 sin<|> - r /. r = zcos9-dsin0sinf-(O-cc)v/CQ . The usual assumption in radiopositioning is that the unknown distance is given simply by
r = - Ov/co , which may be interpreted as the radius of a spherical wave centred at the origin (r=0), that is, the so-called phase centre of the wave coincides with the coordinate origin. The effect of the additional terms may be interpreted as moving the phase origin away from the coordinate origin by an amount which is a function of the element dimensions and the observer-source orientation (9,<|)).
^M^
9.14
fit*/ v^^x
OBSERVATION ERROR;
RULE OF THUMB
WAVELENGTH 300 m 30 m 20 cm 3m 30 cm 2 mm
ERROR SOURCE
SATELLITE ORBIT CLOCK PROPAGATION IONO TROPO RECEIVER MULTIPATH CLOCK NOISE
TYPICAL ERROR
WHAT TO DO ?
difference difference diff. + 2 freq. model & diff. antenna siting model in soln
OBSERVATION ERRORS
The observed time delays to be converted to ranges will have some additional errors associated with them. These are due to the limitations of the receiver's electronics and are of a random character. As a general rule, they are proportional to the wavelength of the signal and consequently have significantly different magnitudes (standard deviations) for C/A-code ranging, P-code ranging, and carrier ranging, as shown in the table.
In addition to these truly observational errors, the random errors in the model have to be taken into account. The model errors are the residual errors that remain in the ranges and in the satellite positions after the biases enumerated in this chapter have been modelled out. These errors are generally much more damaging than the true observational errors, and no effort must be spared to reduce them.
Most of the model errors can be eliminated, or at least drastically reduced, by differencing the observed ranges between satellites, between receivers, and between times. Differencing either eliminates some of the errors because they are common to the differenced ranges (clock errors), or reduces the rest of the errors because of their spatial and temporal coherence (propagation errors, orbital errors). This reduction is more likely to be successful if the two simultaneously ranging receivers are closer together than if they are farther apart. The multipath errors can be neither modelled nor differenced out. The only protection against a multipath error is a careful siting of the antenna with the aim of avoiding reflections off surrounding structures.
9.15
OPTIONS FOR MODELLING BIAS Assume perfectly known (a possible case for non-positioning applications such as time transfer and orbit determination) Assume known with some uncertainty (by varying the uncertainty assumed, this approach is appropriate for all positioning and non-positioning applications)
Assume unknown (the usual case for positioning applications)
EFFECTS ON DETERMINATION
Continental tracking network (to obtain 3 m orbits over the region, relative station coordinates must be known to 0.5 ppm) Global tracking network
?
Each metre of error in station coordinate leads to typically 3 ns time transfer error
Canadian GPS Associates, April 1986
STATION COORDINATES
For non-positioning applications, i.e., orbital determination and time transfer, the positions of antennas (stations) must be assumed known either perfectly or with some uncertainty. The latter case is theoretically more appropriate because the ground positions are never known perfectly and, as such, should always be allowed to adjust even in the estimation of non-position parameters.
It should be obvious that the orbits can be predicted more accurately the better the ground positions are known. Thus, for instance, if an accuracy of 3 metres is desired in the determination of satellite positions, relative station coordinates must be known to 0.5 ppm. In addition, the spacing of tracking stations should be as wide as possible. For the time transfer, an uncertainty of one metre in the coordinates of the two stations induces an error of 3 ns in the transferred time.
Incidentally, in positioning, the accuracy of a priori knowledge of the sought coordinates also has an effect on the finally estimated coordinates. This effect gets smaller and smaller with time spent on tracking.
CHAPTER TEN
SOLUTIONS
10.00
OLUTION,
POSITIONING PROBLEMS Main characteristics of static and kinematic positioning Classification of types of solutions POINT POSITIONING USING PSEUDO-RANGES
Satellite multi-ranging Least-squares solution for point positioning
RELATIVE POSITIONING Using pseudo-ranges Using phase measurements GPS phase measurement differencing
NETWORK SOLUTIONS Orbit improvement and network solution Weighted short-arc solution KINEMATIC POSITIONING Absolute Relative Combined POSTMISSION SMOOTHING CONSTRAINED SOLUTIONS
Canadian GPS Associates, April 1986
Solutions 10.0
SOLUTIONS
In this chapter, we look at alternative approaches to processing the observation equations discussed in Chapter 8, taking into account the biases described in Chapter 9, in order to obtain resulting positions for the tracking receivers.
We begin by considering the kinds of positioning problems that may be encountered, and the processing methods which are appropriate in each case. Generally, these problems can be divided into static and kinematic cases, in which the receivers can be considered to be stationary during the GPS observation campaign, or not, respectively.
The two approaches of point positioning and relative positioning are discussed, in which either a single receiver is used or at least two receivers are used simultaneously. Relative positioning leads naturally into a discussion of network determination.
Kinematic positioning can be accomplished by means of several optional approaches involving either single or multiple GPS receivers or the use of additional observations, for example, from inertial navigation systems. The effect of GPS satellite constellation outages on kinematic positioning, and the role of constraints in overcoming them, complete this chapter.
10.01
LE:
POINT
RELATIVE
NETWORK
Solutions 10.1
POSITIONING PROBLEMS
The three basic positioning problems to be solved are:
point positioning a one point problem; relative positioning a two point problem; and positioning a network of points a many point problem.
Further, in reality these points may not be moving and thus are considered static, or they may be moving and thus are considered kinematic. Examples of the latter are a moving vehicle on land, an aircraft, or a vessel at sea.
10.02
STATIC:
One, two, or several stationary points Results not needed in real time High redundancy, reliability and accuracy
KINEMATIC:
One moving point Position relative to geocentre or Position relative to a second point Results needed in real time Low redundancy and thus lower accuracy Reliability required, thus combination and integration needed Improvement of results by postmission smoothing
Solutions 10.2
10.03
Solutions 10.3
u;
POINT
10.04
Orbit
Solutions 10.4
where the observed range p is equal to the range p expressed as a function of the coordinates of the point (R) and the satellite coordinates (r). The ionospheric and tropospheric corrections to the range are applied and assumed correct. The clock bias c(dt - dT) is due to the time offset between satellite and receiver clocks, and is left in the solution as a nuisance parameter to be solved for along with the position vector R of the point in question.
10.05
MULTI
SCALAR EQUATION
P.
LIMITATIONS
Any error in r^ will result in an error in R Propagation errors will cause errors in r 0
Canadian GPS Associates, April 1986
Solutions 10.5
SATELLITE MULTI-RANGING
The mathematical model in conceptual form is
p = f(R,r,(dt-dT)) ,
where the three coordinates in R and (dt - dT) are the four unknowns, while p is the observed quantity. It is clear from this figure that at least three p's are needed to solve for R, and an additional p is needed so that a solution can be made for (dt-dT) as well. Note that this solution is only possible if the satellite position (r) is assumed known. This means that the ephemeris error is assumed to be zert. This assumption may degrade the accuracy of the solution.
In practice, dt is given a priori via the broadcast message, therefore the solution actually is made only for dT. The corresponding redundancy v of the solution is v = n-4 , where n is the total number of pseudo-ranges observed. A measure of the goodness of the accuracy of the solution is given by means of the covariance matrix C corresponding to R and dT.
In order to show how C is computed, we must now introduce the linearized version of the model f given above, namely:
where w is a vector of numbers obtained by evaluating f with approximate values of the position R and observed values of p; v is the vector of n residuals corresponding to the n pseudo-ranges, while AR and AT are the design matrices corresponding to the two sets of unknowns dR and dT. Note the dimensions of the various vectors and matrices. For n equal to 4, we have the unique solution mentioned above. Also needed is the covariance matrix C of the pseudo-ranges; it is simply a diagonal matrix consisting of variances for each range.
10.06
POINT
MODEL
ITIONING
p = f{R, r, (dt - dT)}
v = (A A,J dR
R
LINERIZED MODEL
NORMAL EQUATIONS
+W
nxl N
nx4
SOLUTION
dR dT
N" X u
R
'R,T
C = R
,2 T
a:
Canadian GPS Associates, April 1986
Solutions 10.6
and
The covariance for the three coordinates of the point and time bias as follows:
r CR I CRT 1
I 3x3
4x4 | CjR
|
|
3x1
Cj
I
I
CR.T-N-I- I - - - - - i - - - - - I
L 1x3 I 1x1 J
where CRy and C jR are cross-covariance matrices, and the matrices of interest are
T ax2
CR =
3x3
L symmetric and
CT = a-r2, 1x1
which define the variances of the four unknowns. These four variances are used to define the various dilution of precision (DOP) measures. The above equations are summarized in this figure.
10.07
USING PSEUDO
ARy = Rj-Rj I
Relative positions A R y are determined more accurately than "absolute" positions RI and RJ , since several sources of error will be common at both receivers, and will cancel.
Canadian GPS Associates, April 1986
Solutions 10.7
10.08
;NTS
Adt = 0
(code)
Solutions 10.8
where two sets of nuisance parameters need be solved. The relative receiver clock bias, AdT, remains as the satellite clock error, dt, cancels out. The effects of ephemerides and refraction errors are almost identical at each receiver, thus, as a consequence of differencing, they are greatly reduced. Note, however, the presence of the carrier beat phase ambiguity AN. Hence, in relative positioning, with carrier phase observations, the inherent number of unknowns is three for AR, one for AdT, and one for the ambiguity parameter a minimum of five equations are needed.
10.09
PHASE MEASUREMENT
SATELLITE p
Epoch t1
poch ti SATELLITE q
t2
RECEIVER i
Instantaneous (Single) Difference over Receiver (A)
RECEIVER j
(Double) Difference Over Satellite
(Triple) Difference Over Time (5)
8VAO
o>:(t 2 )
A<D
AO
Solutions 10.9
The next model in this sequence is the so-called double difference model: VAO = f(AR, r; AN) . Note that the relative receiver clock bias cancels.
The final model in the sequence is the so-called triple difference model: 8VAO = f(AR,r). Note that the ambiguity parameter has been eliminated.
One complication not mentioned above is that there is a price to pay for not having to solve for the ambiguity parameter. The covariance matrix for double and triple differences is non-diagonal, and the computation of its elements is given in Remondi [1985d]. A comforting note is that the solution obtained when using only a diagonal matrix is not adversely affected [Goad and Remondi, 1984]. Another disadvantage of using triple differences is that it is not possible to take advantage of the integer values of the ambiguity parameters. In some cases, however, it is not even possible to enforce the integer value; for example, for long baselines in which the atmospheric conditions are different. An important advantage of using triple differences is that cycle slips can be easily identified, however, when such slips occur, perfectly good undifferenced observations require rejection because the *partner* observation in the difference is missing. This lack of continuity can be resolved successfully during preprocessing.
10.10
AR
ARzs
AR
Provided all
ARy
Solutions 10.10
NETWORK SOLUTIONS
GPS network processing techniques are still in their infancy, and at the time of writing this Guide, it is not possible to distinguish between all possibilities and state the definite way to perform this task. Nevertheless, various researchers have followed one of two procedures depending upon the aim and extent of the network. When dealing with local networks of a few kilometres, but not larger than, say, 100 km, receiver coordinates are the main unknowns with orbit biases possibly appearing as nuisance parameters. The solution for corrections to orbital elements of GPS satellites is left to networks of regional or continental extent. Coming back to local networks, there exist two approaches that have been used, namely: (1) Reduce the observation campaign to simultaneous data observed from pairs of stations (baselines), and then combine these baselines into a network. (2) Take the GPS observations observed simultaneously at all stations directly into a network adjustment where all the coordinates of the network are present as unknowns. The baseline approach is rigorous if there are no common observations and biases between baselines. In this figure, an array of interconnected baselines is shown. If one considers one baseline as the basic unit, then the basic model for a three-dimensional adjustment is the following:
where AR- is the quasi-observed value of the coordinate difference with covariance matrix C^. Note, the role played by R- and Rj is that of unknowns to be solved for. This setup is akin to a levelling network, but in this context we get threedimensional information instead of just height (one dimensional) information. Another possibility is to perform a condition adjustment on the ARs. This means that for each independent vector polygon the sum of the vectors AR must be a zero vector. With the adjusted ARs, one is then able to compute the coordinates of any point, as indicated in this figure.
10.11
ORBIT
true orbit
estimated orbit
Solutions 10.11
Regional networks are those having an extent greater than 100 km and, in some applications, may be continental in extent. When station separations are of this magnitude, then it becomes both necessary and possible to improve the orbit along with the solution for the coordinates of the points in the network.
There are basically two approaches that have been developed for this purpose. One is the so-called orbit bias approach, and the second is the short-arc state vector approach. The mathematical model for both of these approaches can be written as follows: f! (R, r, b, 1) = 0, q, Cb f 2 (r,z) = 0, Cz ,
where we have a combination of models to solve. The first model contains the receiver coordinates R, the satellite positions r, some biases b (e.g., clock), and observables 1 which can be ranges and/or phase measurements. The second model represents a reparameterization of the rs by z a vector of six Keplerian orbital elements. C denotes a covariance matrix.
>LUT:
,r,b,l) = 0 ,
f (r,Z 0 ) = 0 ,C
0
T-
10. 12
{r
= B Z0
ADVANTAGES:
Solutions 10.12
10.13
^9
Geocentre
Solutions 10.13
with redundancy v equal to 0 if four pseudo-ranges to four satellites are observed in order to solve for the three coordinates of R and the one clock bias parameter. Clearly, the goodness of the solution depends upon the geometrical distribution of the satellites.
10.14
Solutions 10.14
The mathematical model for this mode is identical with the static relative positioning model. For pseudo-ranges
Ap = f(AR, AdT), where AR is the unknown relative position, AdT is the unknown clock bias of the moving receiver relative to reference receiver, and Ap is the observed receiver pseudo-range difference. Clearly, observations to four satellites are needed for a unique solution for the unknowns. The relative model using the more accurate phase measurements has the following form:
10.15
POSIT:
COMBINED KINEMATIC POSITIONING
Solutions 10.15
where the first model is the standard point positioning model for pseudo-ranges, while the second contains a time differenced phase (change in phase) measurement yielding information on the change in position between contiguous epochs of pseudo-range determinations. Note the presence of the clock drift parameter, 8dT.
In essence, the above model is a linkage of the discrete position determinations resulting from the kinematic absolute positioning mode described above. The solution of the above equations can be done by sequential least squares. Kinematic relative positioning, using a combination of pseudo-ranges and phase (or Doppler) measurements, makes use of the following two models: Ap = fj(AR, AdT)
10.16
TIME
faQ
iiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiimiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiitiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiniiiii
TIME
Canadian GPS Associates, April 1986
Solutions 10.16
POSTMISSION SMOOTHING
As already stated, kinematic positioning involves the determination of the trajectory of a moving vehicle. The amount of information (observations) on the trajectory increases with time. Positions computed in the past, however, do not benefit from the information collected in the present, and certainly do not benefit from the data to be collected in the future.
The key to being able to go backward in time and improve the results is whether the discrete time epochs have been linked in some way. This linkage has a very specific meaning, and that meaning is exemplified by those models f, discussed earlier, that involved the 8 operator. These included combined point positioning, combined relative positioning, and integrated point and relative positioning. What we had in these models was the ability to predict the position at ti+1, and determine (filter) the position at tk+1. This task can be performed with the well-known Kalman (or Bayes) filter equations. The literature also contains the necessary equations that allows one to back solve (reverse) the filter equations and obtain updated (smoothed and improved) estimates of positions in the past from information that was collected since its filtered position was determined. Clearly, one would expect the variances corresponding to the smoothed solutions to be smaller than the filtered positions.
10.17
>LUTION;
WHAT OPTIONS DO WE HAVE ? Unconstrained solution for 3D position and time Height-constrained solution for 2D position and time Time-constrained solution for 3D position Height-and-time-constrained solution for 2D position
WHY [STRAIN OUR SOLUTION ? Improved geometrical strength improved position fix accuracy
WHY NOT CONSTRAIN OUR SOLUTION ? Must know "true" value of constrained parameter using wrong value may lead to much worse fix accuracy than the unconstrained solution TIME CONSTRAINTS Useful constraints are 10 ns or smaller Only time observatories keep 10 ns time May be useful to "coast" through outage Cesium clock may keep 10 ns for one day Rubidium clock may keep 100 ns for one day Crystal clock may keep 10 jus for one day
CONSTRAINTS Useful constraints are 5 metres (15 ns) or smaller can we know our height to 5 m in a launch ?
Canadian GPS Associates, Apnl 1986
Solutions 10.17
CONSTRAINED SOLUTIONS
The geometry for a four-dimensional fix will be worse than that for a fix where we actually know one or more of the solution parameters in advance, and take this knowledge into account. We can use a technique known as weighted constraints we let our solution vary from the assumed 'known' values of the weighted parameters only slightly. Constrained solutions are particularly useful to help the solution 'coast' through a rough period (that is, an outage in the geometry). However, there can be a price to pay: if the constrained parameter is weighted at an inappropriate value ('wrong' by much more than the weight constraint would indicate), then the remaining solution will be adversely affected. Height constraints are a useful approach, for example, in the case of a conventional hydrographic survey from launches. Height constraints are quite effective at the 5 metre level, and 'knowing' our height to 5 metres, when at sea level in a survey launch, seems quite feasible. This figure shows how effective this can be in eliminating satellite outages. In the following three figures, we consider this c^se study, and look at how feasible it is to constrain height to five metres realistically.
Solutions 10.18
Height (H) of the antenna above sea level: This can be measured directly or indirectly in port to an accuracy corresponding to a standard deviation GH smaller than 10 centimetres. However, it will change with varying ship loading, according to the characteristics of the vessel.
Height (N) of the geoid above the reference ellipsoid: The geoidal height, N, may reach as much as 100 metres. The best geoidal height information available today is good to about 2 metres (one sigma level) in an absolute sense with local relief being more accurate.
Roll and/or pitch angle (co): If we wish to account for roll, it either must be modelled in the solution, or has to be monitored; here we assume it can be monitored to an accuracy high enough so that we do not have to worry about the effect of o .
Tide and sea state (TS): These have to be predicted from models. The combined magnitude may reach many metres. However, if good models are used for the prediction, its accuracy would be some 20 centimetres (one sigma level).
Semi-permanent sea surface topography (SST): This is a rather difficult quantity to get hold of. Estimates exist [Lisitzin, 1974; Levitus, 1982] but it is somewhat uncertain just how good they are. Probably <*SST = 20 cm would again be a fair assessment. The magnitude of SST is below 2 metres so that if the SST is considered to equal to zero, OSST would have to be raised to some 80 centimetres.
The deflection of the vertical (0): This is normally much smaller than one minute of arc and can be safely neglected.
ON
C O
5 CUu* sL-OO
o o
so *- g
l-< *PN
C C
^ S^ Pu. g ffi U
ss
S ,2
M 10
H
t==]
^ ^ 0>
o>
ON O>
^ W3
t-
2 Dfi
fi _,
OA D
S E o o
o>
ii
>
o
H
{==)
A L
S E
10
o
10
10
a o P* CA
SO O
C fc* ,
S S
M ^
Solutions 10.19
Geoidal heights expressed as high order and degree spherical harmonic expansion: At least two sets of global potential coefficients that come to mind would satisfy the accuracy requirements: The Ohio State University (OSU) 180 by 180 model [Rapp, 1983] contains 32 761 coefficients. The evaluation of one N value from the OSU model takes about 5 minutes on an HP 9816 computer, once the coefficients are in the memory. A low degree expansion like the GEM 10B [Lerch et al., 1981] will not necessarily satisfy the above requirements, though its 1369 coefficients would require much less storage space. Geoidal heights expressed on a regular geographical grid: Gravimetric geoids, such as the Geodetic Survey of Canada geoid [Lachapelle, 1977] or the UNB geoid [VaniSek et al., 1986], can be precomputed on a regular grid, say 10' by 10'. Linear interpolation (very fast) between the grid nodes would suffice to give an adequately accurate N anywhere within the boundaries of the solution. The number of grid values (36 per 1 by 1 cell) to be stored would depend on the area covered by the navigating ship.
Satellite altimetry expressed on a regular geographical grid: Gridded OSU SEASAT altimetry [Rapp, 1982] (heights of short period mean sea level above a reference ellipsoid) on a 10' by 10' grid and NASA 7.5* by 7.5' gridded SEASAT/GEOS-3 altimetry [Marsh et al., 1984] are available globally. These data sets would be really better suited for the antenna height determination because they automatically take care of at least part of the SST.
CN O
Solutions 10.20
CHAPTER ELEVEN
STATIC APPLICATION;
11.00
STATIC APPLICATIONS: Geographical positioning applications Impact of GPS on control networks Height and geoid determination applications
RESULTS: Ottawa test network Western Canada triangle test Alberta township corner Urban tests in Quebec City Gimli array test Long baseline test Gaithersberg test
Static means that the points being positioned have no perceptible motion associated with them relative to surrounding points. All geographical positioning applications fall into this category in which the goal is to attach coordinates to points or features of interest.
Quasi-static means that while there is perceptible relative motion among the points being positioned, it is generally so slow that it would not be noticed during a GPS survey campaign (a few days to weeks), but may become evident between campaigns (months to decades apart). Geodynamics and crustal movement monitoring applications fall into this category.
GEOGRAPHICAL POSITIONING Attaching coordinates to objects of interest: e.g., control monument, gravity data point, census address, etc.
ONE RECEIVER + "ACTIVE CONTROL SYSTEM" Relative coordinates COMPARISON TO OTHER TECHNIQUES
IMPACT ON LAND SURVEYING
With the advent of GPS, accurate coordinates will be easily obtained for a variety of geographical applications. In addition to fulfilling the usual requirement of accurate coordinates for survey control monuments, GPS can provide equally-precise coordinates for other objects of interest, such as gravity data points, water sample points, and addresses for census purposes.
In the near term, the greatest impact of GPS in geographical positioning will be on the surveying industry. For many jobs, a pair of GPS receivers can replace a conventional total station instrument and greatly exceed its capabilities. With GPS, there is no intervisibility requirement and practically no limit to the distance between receivers. Furthermore, no reference azimuth is required as the coordinates obtained are relative three-dimensional coordinates of one point with respect to the other.
In addition to their high accuracy, GPS surveys will be carried out more cheaply than conventional surveys. For some types of surveys, this is even true at the present time with only the prototype constellation available. No line cutting or survey tower construction is required for GPS observations. It may be necessary to clear a small area surrounding a GPS antenna in order that the sky not be obscured, however, this is generally neither time consuming nor expensive. As a result, survey practice will be considerably changed with the adoption of GPS technology.
A further benefit accruing from the use of GPS is that all coordinates are in the same unique system. This characteristic considerably simplifies the use of position information in different applications and makes the creation of land information data bases that much easier.
GPS will have a particularly significant impact on survey control networks. This impact will be discussed in detail below.
11.02
WHAT INFORMATION DO COORDINATES CONTAIN? Intrinsically - zero information Information depends on attaching something to the coordinates
WHAT IS A CONTROL NETWORK ? A set of monuments with attached coordinates
WHAT IS IT USED FOR ? Transferring coordinates to other features Boundary markers Map sheet corners Highways, power lines, etc.
WHY ARE CONTROL NETWORKS INDEPENDENT FROM FEATURES OF REAL INTEREST ? Limitations of control network survey and maintenance techniques
11.03
= ORTHOMETRIC HEIGHT
'/Topography
N=h-H
( Medium wavelength geoid ^ ( Detailed geoid ^ I (based on gravity anomalies)! [(includes terrain effects))
Canadian GPS Associates, May 1987
GPS is a three-dimensional positioning system which can naturally provide heights and/or height differences above a selected ellipsoid. On the other hand, orthometric heights can be measured accurately by geometric levelling techniques. However, this is a slow and expensive process. If orthometric heights are known at satellite positioned stations, a direct and accurate measure of the geoid can be obtained by the following equation (see the figure):
N=h-H.
Geoidal undulations (N) can also be estimated by gravimetric and other methods. Until now, it has been difficult to obtain reliable accuracy estimates for these methods. A comparison of gravimetrically-derived and satellite-derived undulations provides a reliable accuracy estimate for these quantities. For instance, this technique was used, with Transit-derived geoidal undulations, to assess the accuracy of the Canadian geoid during the late 1970s [Lachapelle, 1979], Due to the low accuracy (50 to 100 cm) of the Transit-derived undulations, it was not possible to estimate die accuracy of the gravimetric geoid better than within the above range. The use of GPS-derived undulations nowadays permits an assessment of gravimetric undulation differences with an accuracy of the order of 5 to 10 cm. Recent tests performed in Canada (e.g., Schwarz et al. [1985]) show that the accuracy of the gravimetric geoidal height differences is currently of the order of 10 cm in the areas tested. This is very important as the various methods of computing the geoid (e.g., Lachapelle [1979], Schwarz et al. [1985], Vanicek et al. [1986]) can be more fully analysed and refined to produce a better result.
This, in turn, will lead to the combination of GPS-derived ellipsoidal heights and available geoidal undulations to derive orthometric heights in areas where geometric levelling operations are impractical. If the above gravimetric geoid accuracy of 10 cm can be confirmed over most or all of the country, it would become possible in the near future to establish orthometric heights with an accuracy of the order of 10 to 20 cm using GPS. This accuracy could be further improved by additional land and satellite-based gravimetric data and more refined geoid estimation methods.
11.04
QUASAR
QUASAR
VLBI Monitor
sites
GEODYNAMICS APPLICATIONS
The accuracy of baseline lengths determined from GPS carrier phase observations has already reached the level of 0.1 ppm. A relative error of 0.1 ppm translates to 3 cm on a baseline of, say, 300 km. This accuracy level implies that satellite orbital errors have been reduced to about 2 m. Such improved orbits can be obtained by simultaneously tracking the satellites from a network of VLB I sites, as well as from the sites of interest. With the ability to measure baselines of several hundred kilometres to an accuracy of a few centimetres, changes in baselines resulting from crustal plate motion, such as that along tectonic faults, can be monitored. GPS is therefore becoming an important tool for geophysicists attempting to understand better how the earth works.
The motions associated with plate tectonics are very small, generally of the order of a few centimetres per year. Nevertheless, such motion can be measured if a network of GPS receiver sites is established in the region of interest and the network is remeasured occasionally over a suitably long period of time. Such networks could be set up across fault zones, in areas of subsidence, or across spreading zones. The Jet Propulsion Laboratory (JPL) under the auspices of NASA has undertaken a campaign to measure the complex tectonic motion of parts of Mexico and the Caribbean using GPS [Thornton et al., 1986]. Another application of GPS in the study of geodynamics is the international campaign for monitoring the crustal motions in Iceland.
If an earthquake occurs in a region being monitored by GPS, it should be possible to accurately measure the co- and postseismic deformation. Such measurements could provide valuable clues in understanding the mechanism of earthquakes, tectonic or other. The use of GPS to monitor the bulging of volcanos preceding eruption has also been suggested [King et al., 1985]. Although such monitoring can be done using conventional surveying techniques, it may be awkward, time consuming, and more importantly unsafe. GPS data from individual receivers placed on the flank and in the crater of a volcano could be telemetered to a central site at some distance where the volcano's deformation could be monitored safely.
11.05
MOVEMENT MONITORING
GPS, thanks to its three-dimensional accuracy, can be used for a wide variety of man-made movement monitoring tasks, ranging from surface subsidence caused by sub-surface mining to hydro-electric power dam deformation. Three specific applications are described herein.
Hydro-Electric Dam Deformation: This type of structure is constantly being deformed due to the load exerted on the earth's crust by the water behind the dam. In a large hydro-electric complex, the continuous monitoring of horizontal and vertical deformations may be required for safety reasons. GPS may be the most cost effective tool for this task. Firstly, it is capable tf meeting the stringent accuracy requirements (1.0 to 0.1 ppm). Secondly, an array of GPS receivers could be installed permanently at proper sites and controlled remotely. In distant areas, the data could be regularly transmitted by means of an appropriate telemetry link to a processing centre for reduction and analysis.
Offshore Platform Subsidence: This GPS application has arisen in connection with the subsidence of petroleum exploration platforms in the North Sea which is due to the extraction of oil and gas [Collins, 1986]. The rate of subsidence can typically range between 10 and 50 cm per year. The platforms are periodically surveyed, say, every month to every six months. Each survey is conducted over a period of several days to insure a satisfactory short-term repeatability. The accuracy requirements being very stringent, orbital and multipath effects can be significant. Nearby platforms, assumed stable, are often used to obtain a higher accuracy through relative positioning.
Surface Subsidence: This potential application arises in connection with subsurface mining, water, or oil withdrawal. The use of GPS can be very cost effective in such a case, since no interstation visibility is required. An array of GPS receivers can be deployed and moved quickly between sites of interest. Once the GPS constellation is complete, an array of n receivers could monitor (n-1) x 8 sites per working day (of 10 to 12 hours). Weather dependance would be decreased significantly as compared to classical survey methods.
VO
H H
Macrometer V-1000 results: The carrier beat phase observations were adjusted using three different programs: Macrometries' software for baseline vector solutions, UNB software for baseline solutions, and UNB software for network solutions (the last two being the third generation UNB software package DIPOP, used in two of its modes). All three results agreed at the 1 ppm level. The network solution is the most reliable one since it combines all observations at all stations to derive a unique set of station coordinates. Thus nonsystematic errors due to neglected ionospheric propagation delays (single frequency observations!) and other sources are reduced, at least partially.
Texas Instruments TI-4100 results: The dual frequency observations were adjusted at UNB in the network solution mode, also using DIPOP. The adjustment was done separately for the LI carrier beat phases only and the ionospheric delay free linear combination of LI and L2 phases. Both results agreed at the 1 ppm level.
Comparison of V-1000 and TI-4100: A comparison of the two network solutions (dual frequency result for TI-4100) showed an agreement at the 1 ppm level for the baseline lengths and height differences. The horizontal coordinate differences between the two solutions exhibited a significant rotation of about 0.5", thought to be due to orbital errors.
H m pq H
Each daily collected data set was processed independently using program PHASER [Goad, 1985]. The broadcast ephemerides were utilized. Carrier phase measurements at 30 second intervals were used. For the 3 days during which dual frequency observations were taken, LI only, L2 only, and ionospherically corrected (L1/L2) measurements were successively processed to isolate the effect of the ionosphere. The stations utilized were not tied to existing networks, and the daily solutions were compared with the average L1/L2 solution based on the three days of dual frequency observations. In general, the single frequency results are consistent to within 2 to 3 ppm in terms of maximum scatter, while the dual frequency result consistency is of the order of 1 ppm. The effect of the ionosphere is thus of the order of 1 ppm. This is significant for accurate applications. The observations took place during the middle of the night when the effect of the ionosphere is expected to be minimal. Also, the sunspot activity during 1985 was near a minimum. A maximum will be reached in 1991. At that time, the effect of the ionosphere could reach several parts per million.
oo
y)
11.09
'\
Canadian GPS Associates, April 1986
1 1000 m
A 16-station network located in the Quebec City area, with baselines ranging from 200 m to 2 km, was established in 1984 and re-observed in 1985 using successively the Macrometer V-1000 and TI-4100 receivers. The purpose of this project was to assess the capability of GPS to meet the Quebec survey specifications for this application under a wide range of climatic and multipath environments [Moreau et al. 1985; Moreau 1986]. The Macrometer V-1000 test was conducted in January 1984. Three receivers were used to obtain 4 quasi-independent baselines each day. The field observations were completed in 4 days. The data were processed using both Beutler's [1984] software and that developed by the contractor; the agreement between the results was of the order of 1 mm, which was considered fully satisfactory. The external consistency of the GPS coordinate differences with the local network was of the order of 4 mm (r.m.s.) horizontally and 5 mm vertically.
The TI-4100 test was conducted in March 1985. Again, three receivers were used to obtain 4 quasi-independent baselines each day. The field observations were completed in 7 days. The data were processed by the sub-contractor. Comparison of the GPS results with the local network revealed some inconsistencies on a few baselines which were subsequently rejected from the final adjustment. After rejection of these baselines, the average agreement with the local coordinates was of the order of 1 cm. Overall, the results were not as satisfactory as the Macrometer V-1000 results.
The general conclusion is that GPS will become a cost effective tool for urban densification once the complete satellite constellation is available.
^-_^'
ON OS
X I tfl LU
tt LU
H
pq
t=J
cn O
0)
HI
"3
o
_J
CD U. <
o o
Most results that have been obtained with GPS have been on relatively short baselines ranging from less than a kilometre in length to one or two hundred kilometres. However, several tests have been carried out with baselines as long as several thousand kilometres. Such long baseline tests are useful for separating and understanding the various biases affecting GPS observations and for accurately modelling the orbital biases of the satellites.
One such long baseline test was the Spring 1985 High-precision Baseline Test organized by JPL. Between 29 March and 5 April, dual-frequency GPS receivers were operated at 10 sites in the continental U.S. Six of the sites have permanent VLB I radio telescopes near by. Another site (Mammoth Lakes) has been visited by a mobile VLB I station. All 10 sites were occupied by Texas Instruments TI-4100 receivers. At two sites (Mojave and Big Pine), JPL's SERIES-X receivers were also operated, and at Westford, Richmond, and Fort Davis, U.S. Air Force Geophysical Laboratory dual-frequency receivers (the forerunner of the Macrometer II) were additionally operated. The receivers were used with a variety of frequency standards: hydrogen masers, cesium and rubidium frequency standards, and receiver internal oscillators. Water vapour radiometers were operated at three of the sites in California: Hat Creek, Mojave, and Big Pine.
A number of research groups have analysed the data from the test. Preliminary results of the analyses of data from all three types of receiver indicate agreement in baseline lengths with those determined by VLB I to about 0.1 ppm [University of Texas, 1986]. In order to achieve this accuracy, the satellite orbits were adjusted using the data from three or more of the stations. Use of the broadcast ephemerides alone resulted in baseline errors up to 2 ppm.
11.12
SEMI-KINEMATIC GPS
Eleven baselines One hour total duration Initial baseline in 7 minutes 45-105 seconds for other baselines RESULTS: Worst case 3.5 cm away from reference Typically 1.2 cm away from reference
Canadian GPS Associates, April 1986
CHAPTER TWELVE
12.00
ULT,
Kinematic positioning models Comparison of kinematic models Real-time vs postmission processing Accuracy and instrumentation Range of applications
Land vehicle test traverse Land vehicle positioning results Land vehicle velocity results
SHIPBORNE APPLICATIONS Hydrographic survey test area Hydrographic survey test results 3-D seismic surveys
AIRCRAFT
Laser bathymetry Aerotriangulation without ground control Airborne gravimetry Helicopter applications
PRINCIPLE OF INEMATIC
ABSOLUTE KINE]
SITIONING
12.01
Kinematic positioning refers to applications in which the trajectory of a moving object is determined. In the applications considered here, the moving object is a GPS receiver carried by a road vehicle, a ship, or an aircraft. Slow kinematic carriers, like glaciers and tectonic plates, have been treated in Chapter 11 and will not be considered here. Kinematic positioning is closely related to navigation. In the first application, the positioning aspects are emphasized, while in the second application, the guidance aspects, i.e., speed and heading, are of primary interest. In addition, real-time results are required in navigation, while postmission results are quite acceptable in many kinematic positioning applications. Although the two areas overlap to some extent, the above distinctions are useful from an application point of view.
Kinematic positioning with respect to an implied coordinate system, as, e.g., the conventional terrestrial system, will be called absolute kinematic positioning. An alternative term often used is kinematic point positioning. Positioning with respect to one or several receivers at known positions, which track the same radiopositioning signals as the moving receiver, will be called relative kinematic positioning. An alternative term often used is differential kinematic positioning. Real-time relative kinematic positioning requires a radio link between the stationary and the moving receiver. The principle of absolute and relative kinematic positioning is shown in this figure. Because of the simultaneity of the measurements and the resulting reduction of satellite errors, relative positioning is more accurate than absolute positioning. This should be understood such that the vector AR can be determined more accurately than the vector R.
GPS measurements used for absolute and relative kinematic positioning are pseudo-ranges and phase measurements or a combination of both. GPS measurements can also be combined with observations coming from other sensors. The resulting models will be briefly reviewed in the next section.
KINEMATIC
12.02
80 = f {5R, 8dT}
2
v =
(v , 8v }
EELATEVE
KINEMATIC
PHASE ONLY :
COMBINED:
A8<D = g. {A8R,A8dT} 2
g and g as above
8v>
gl and g 2 as above
X 2
INTEGRATED:
81 v =g.
In the absolute mode, pseudo-range kinematic positioning uses pseudo-range measurements p to at least four satellites to determine the receiver position R and the clock bias dT. In the relative mode, measurements are differenced over receivers, and p, R, dT are replaced by Ap, AR, AdT.
Phase kinematic positioning is intrinsically a relative method [Remondi, 1985b]. It uses phase measurements to at least four satellites from both receivers to
determine coordinate differences with respect to the known initial position of the moving receiver and the clock drift A8dT between the two receivers.
Combined kinematic positioning uses a minimum of four pseudo- ranges p and four phase measurements O to determine the receiver positions Rj and Ri+1
the receiver clock bias dT, and its drift 5dT [Kleusberg, 1986]. Absolute and relative positioning are handled as outlined above.
Integrated kinematic positioning adds a dynamic component g3 to the measurements, resulting in the determination of a practically continuous trajectory [Wong and Schwarz, 1983]. In case of an INS, for instance, velocity measurements l y are available about every 50 milliseconds, and thus the kinematic point density improves dramatically. This has considerable advantages if precise interpolation is required.
12.03
GPS ALONE: PSEUDO-RANGE KINEMATIC POSITIONING No vessel dynamics Trajectory as string of independent position fixes Equivalent to static case Absolute mode makes for simple system Relative mode reduces errors PHASE KINEMATIC POSITIONING No vessel dynamics Trajectory as sum of 8 vectors Sequential least squares Relative mode reduces errors Problems with cycle slips COMBINED KINEMATIC POSITIONING Similar to phase method trajectory as combination of position fixes and 5 vectors INTEGRATED SYSTEMS: INTEGRATED KINEMATIC POSITIONING Full vessel dynamics (INS) Practically continuous trajectory Kalman filter and smoother for dynamic system Best accuracy and reliability Complex system
Canadian GPS Associates, April 1986
All models, with the exception of the integrated one, use GPS measurements only. This makes for a simple measurement system which is usually also advantageous in terms of costs. Compared to the more elaborate integrated systems, however, it has drawbacks in terms of accuracy, reliability, and interpolation.
In the simplest case, absolute pseudo-range positioning, only one receiver is needed. The trajectory is modelled by a string of independent point fixes which currently are seconds apart. The vessel dynamics are not taken into account, and simplifying assumptions are necessary for interpolation. The accuracy is basically the same worldwide and is sufficient for most navigation purposes. It is, however, not adequate for many geodetic and geophysical applications. The simplicity of the measurement system is thus paid for in terms of accuracy and reliability. The other absolute positioning models (combined, integrated) will give a smoother but not a more accurate trajectory.
The accuracy is considerably better for relative pseudo-range positioning. Although the modelling technique is the same, the elimination of a number of common error sources improves the accuracy of the results. However, the accuracy now becomes dependent on the distance between stationary and moving receivers, and the simplicity of the original configuration is gone.
An even better accuracy is possible with relative phase positioning. Differenced phase measurements provide a measure for the receiver position change between epochs. Since the noise of the phase measurements is much lower than that of pseudo-ranges, the accuracy of this positioning method is potentially higher. It requires, however, the absence of cycle slips.
Relative combined positioning goes one step further in terms of modelling and operational reliability. Combining pseudo-range and phase measurements by a filtering mechanism makes the detection of large cycle slips possible as long as they are not correlated with errors in the pseudo-range measurements.
Finally, relative integrated positioning has the potential to provide higher accuracy, improved reliability and precise dynamic information, albeit at the price of more system complexity and higher costs. The addition of an inertial system, for instance, gives a more accurate interpolation between GPS fixes, helps to eliminate cycle slips, increases relative positioning accuracy, provides precise velocity information, and improves overall system reliability.
12.04
VERSUi POSTMIi
REAL-TIME PROCESSING: ADVANTAGES Results available immediately (navigation, guidance) No data storage problems
DISADVANTAGES Less accurate Measurements often not stored (for postmission computation) Blunder detection difficult Real-time link needed for relative methods
1ING
POSTMISSION COMPUTATION: ADVANTAGES More accurate Blunder detection easier Detailed analysis of residuals possible Model identification No real-time data link needed
DISADVANTAGES Results cannot be used in the field Data storage problems
Canadian GPS Associates, May 1987
12.05
INSTRUMENTATION
MODEL ABSOLUTE
INSTRUMENTATION REQUIREMENTS
20-30 15-20
15-20 15-20
1 receiver
1 receiver plus inertial system
RELATIVE
Pseudo-range
Combined
5-10
2 receivers 2 receivers
0.5-1.5 0.2-0.3
Integrated
The first important conclusion is that adding relative position information, either in the form of phase differences or in the form of inertial velocity measurements, does not change the absolute positioning accuracy to any great extent. The reason for this is obvious. The major error sources in absolute kinematic positioning are, at present, orbital errors and ionospheric refraction effects. Both are not eliminated using carrier phase or inertial system measurements. The increase in accuracy for the combined and integrated method using a C/A-code receiver comes from a smoothing of the noisy data by the information provided by the inertial system or carrier phases, not from an elimination of the above errors. The same effect is present in the P-code receiver results but is so small that it does not change the accuracy estimates.
The dramatic improvement in accuracy achieved by relative kinematic positioning shows the potential of these methods for geodetic applications. It is remarkable that, at present, no real difference between P-code and C/A-code results exists. This will change, however, once SPS is degraded. The accuracy estimates for the combined kinematic method have been obtained under the assumption that cycle slips can be satisfactorily eliminated from the data. Although higher accuracies have been predicted by several authors, test results at present do not confirm these predictions. The results for the integrated kinematic method are taken from a simulation study and need further confirmation.
12.06
10m
HYDROGRAPHIC SURVEYING AIRBORNE BATHYMETRY
3-D SEISMIC AEROTRIANGULATION GRAVITY VECTOR MAPPING SYSTEM AIRBORNE GRAVIMETRY AIRBORNE LASER
1m
RANGE OF APPLICATIONS
Although the full impact of GPS will only become apparent in a decade or so, many new and interesting applications are envisaged right now and are being actively pursued by research groups in many countries. The following overview is not complete. It concentrates on work that is currently in progress in Canada. Following this figure, the applications will be subdivided into three major groups. In the first group, low accuracy applications will be treated which require positioning at the 20 to 30 m (lo) level and can be handled by absolute kinematic positioning. In the second group, medium accuracy applications in the 5 to 10 m (la) range will be treated which can make use of the relative pseudo-range model. In the third group, high accuracy applications in the 0.5 to 2.0 m range will be discussed which require the use of either relative phase, relative combined, or relative integrated kinematic positioning.
Low accuracy applications will be sufficient for most conventional marine and airborne navigation tasks and will simplify to a considerable degree the positioning methods used in airborne geophysical surveys. New applications are emerging in the land vehicle mode where the positioning and guidance of truck fleets, ambulances, and taxis pose interesting problems. These applications will become economical once the price of C/A-code receivers drops to a few thousand dollars.
In medium accuracy applications, GPS will most likely replace a number of hydrographic survey systems currently in use. Launch positioning and bathymetry are two of them. Since real-time results are usually required, the development of a reliable radio link is important for these applications. In high accuracy applications, GPS, in the phase, combined, and integrated mode, will give rise to a variety of dedicated systems which will solve problems which could not be adequately handled by conventional methods. Among them are aerotriangulation without ground control, airborne gravimetry and gradiometry, and direct gravity and geoid determination from a moving land vehicle. It will also facilitate a number of applications for which present survey technology is marginal or cumbersome, such as, e.g., three-dimensional seismic and airborne laser profiling. The decision whether to use GPS alone or the integrated method will, to a large extent, depend on how well the cycle slip problem can be solved.
LAND
LOW ACCURACY APPLICATIONS
12.07
Stationary receiver
Station to _\ be positioned
Level surface
Canadian GPS Associates, April 1986
In low accuracy applications, the receiver is one component of a system which integrates the real-time positioning capability of GPS with stored map information to either display the current position of the vehicle on an electonic map or to guide the vehicle to a specific location along a traffic optimized path. Realization of such a system is dependent on low cost receivers and is currently envisaged for ambulances, truck fleets, taxis and, perhaps, railway cars.
High accuracy applications are currently discussed for precise positioning and for gravity field determination. In the first application, GPS will be used in the relative phase mode. The survey operation will be similar to that used in inertial surveying, i.e., the vehicle will start from a known position, drive to the unknown position, and stop there for a short period of time before driving on to the next unknown position [Remondi, 1985a; 1985b]. Gravity vector mapping is possible by using an integrated INS/GPS. In this application, the relative combined GPS mode is used to feed precise position information to the INS. The inertial sensors can then be used to detect changes in the gravity field and to plot a continuous profile of gravity, deflections of the vertical, and the geoid [Schwarz, 1986].
12.08
' 43 km
CONTROL Regular zero velocity updates at two minute intervals 15 second-order control points
KINEMATIC CHECK POINTS Number : 360 Accuracy : 0.3 to 0.8 m (la) in position 1 to 3 cm/s (la ) in velocity
Canadian GPS Associates, May 1987
^j
iRx _-^^
^.M^
LAND VEH:
LJo
IB
12.09
ERROR IN <P K h
RESULTS Mean (m) RMS (m) -4.0 6.2 6.6 0.4 0.9 -2.0 6.4 9.5 13.2 6.6 6.3 7.2
9 'k h 9 A, h
12.10
3 o
1? js <+^
m/s
0
UT 03:00 -0.1 cm/s east = -0.8 cm/s 2.0 cm/s east = 3.2 cm/s
16:30
17:00
17:30
2.9 cm/s east = 9.3 cm/s 7.6 cm/s east =13.4 cm/s
12.11
APPLICATIONS:
ACCURACY REQUIREMENTS
Navigation in open waters 1 km Navigation in coastal channel 50-100 m Hydrographic surveying 2-5 m Three-dimensional seismic surveys 2-5 m Gravity surveys <10 cm s"1
PROBLEMS:
SHIPBORNE APPLICATIONS
Marine uses of GPS cover the whole range of accuracies achievable with the system. In this section, only shipborne applications will be treated, while airborne marine applications will be discussed in the next section. One problem common to most shipborne applications is the mounting of the antenna high above the centre of gravity of the vessel. This results in large antenna accelerations which may result in frequent cycle slips. Thus, the use of the phase method in this environment depends on an operational solution to the cycle slip problem. In addition, the antenna movements make it difficult to determine a precise trajectory of the vessel. Therefore, the high accuracies that can be achieved in land vehicles and aircraft require additional instrumentation in shipborne applications.
Low accuracy^ applications of GPS in the marine environment comprise most navigation tasks in open waters. For this, GPS in the absolute positioning mode is an ideal system. Its accuracy and global availability make it superior to the systems currently used for this purpose, and the development of rather inexpensive C/A-code receivers will soon make it very economical as well.
Medium accuracy applications of GPS using the relative pseudo-range mode are of considerable interest to hydrographic surveying. Real-time accuracies of 5 m or better are required in real time. Tests have shown that the relative pseudo-range mode is marginal in this case and that the use of phase measurements is required. This would result in the difficulties mentioned above or require the use of a simple inertial strapdown system to smooth the noise in the relative pseudo-range measurements.
High accuracy applications are, at present, mainly needed in marine geophysical exploration. Positioning accuracies of 2 to 5 m (lo) are required for three-dimensional seismic surveying, while velocity accuracies of 10 cm s"1 or better are needed in marine gravity surveys to keep the error in the Eotvos correction below 1 mGal.
12.12
44-35
44-20 64-00
The tests in 1983 were conducted to assess both absolute and relative pseudorange positioning with fast switching single channel receivers. In this case, the reference trajectory was obtained from a Syledis radionavigation system, whose accuracy was estimated to be in the 3 to 5 m range. The baseline length, i.e., the distance between stationary and moving receivers, was about 10 km on the first day and about 50 km on the next two days. The total time available for each test was again about four hours.
12.13
ULT;
Number 4-
of visible satellites
Cfl
^ 1
100-r
m
0
O
-100-U
100-r
m
d
d
0>
-100-L
12
13
14 15 16 UT hour of day
The equipment for the second test [Lachapelle et al., 1984b] consisted of two Tl-4100 multiplexing receivers, performing dual frequency P-code pseudo-range and phase measurements on up to four satellites. Pseudo-range measurements, recorded at intervals of 3 to 6 seconds, were used for the tests. Results showed that absolute pseudo-range positioning using the precise positioning system (P-code) is possible at the 15 m (1 a) level given good geometry for three or four satellites. Under similar circumstances, the 5 to 10 m (1 a) range was achieved for relative pseudo-range positioning. Another important result from these tests is that the use of the LI pseudo-range only produces better differential results than ionospherically corrected data. This implies that similar results can be achieved with C/A-code receivers. This was confirmed by later tests.
12.14
-D SEISMIC
PRINCIPLE hydrophone array tail buoy
seismic source
seismic reflector
BOUNDARY DELINEATION
L is 25m to 100m
Canadian GPS Associates, April 1986
AIRCRAFT
12.15
Mounting of antenna
AIRCRAFT APPLICATIONS
Airborne applications of GPS again cover the whole range of achievable accuracies. However, much interest is at present directed towards high accuracy methods at the level of 1 m (1 a) or below because they open the way to new applications in photogrammetry, geodesy, and geophysics which could not be considered up to now. Among them are aerotriangulation without ground control, airborne gravimetry and gradiometry, and airborne laser profiling.
From an operational point of view, fixed wing applications and helicopter applications should be distinguished. Fixed wing aircraft operate in a rather stable dynamic environment and are well suited for most of the above applications which require constant flight altitude and regularity of procedure. Cycle slip problems may occur during takeoffs and turning manoeuvres. On the other hand, helicopters experience a more varied dynamic environment and, therefore, may be more prone to cycle slip problems. Their flexibility in terms of hovering, landing, manoeuvering and velocity, results in a wide range of GPS positioning applications.
In both cases, the mounting of the antenna on the aircraft requires special consideration. On fixed wing aircraft, it is usually installed on the upper part of the fuselage in a location which minimizes shading and multipath effects. Permanent modifications to the aircraft are usually made. The problem is more difficult for helicopters. If the antenna is installed under the main rotor, signal interference will inevitably cause frequent losses of the GPS signal. On large helicopters, it may be possible to install the antenna on the top of the main rotor shaft, above the rotor. However, these types of helicopters are relatively expensive to operate. On smaller types of helicopters, the upper extremity of the vertical tail fin may be well suited for a GPS antenna due to its location above the vertical tail rotor of the helicopter. The site is horizon clear on many types of helicopters and losses of signal can be avoided altogether. The modification of the vertical tail fin is, however, relatively complex.
12.16
LASER BATHYMETRY
Hsi
VJ
ocean bottom
li
depth = - Hsi
ACCURACY REQUIRED :
METHOD :
15 m horizontally
LASER BATHYMETRY
One of the airborne systems that has been developed for remote depth sounding of shallow waters is the Lidar bathymeter. It consists of a dual frequency laser with one frequency chosen such that it is reflected at the surface of the water, while the other penetrates the water and is reflected at the ocean floor. The water depth thus can be derived as the difference of the two elevations measured in the airplane. To locate the measured depth profile, the horizontal position of the aircraft is needed as a function of time. Although navigation systems are presently available which will give the required accuracy of 15 to 20 m, the GPS was tested for this application because of its availability in remote areas. The experiment took place in 1983 in the Lake Huron test area of the Canadian Hydrographic Survey [Lachapelle et al., 1984a]. Two TI-4100 were used, one installed in a DC-3 aircraft, while the second was set up about 100 km away. The accuracy of the results was assessed using four transponders of a Del Norte Trisponder System (UHF), which provided reference positions with an accuracy of 1 to 3 m (1 o) depending on geometry. Aircraft speeds during the test were between 210 and 300 km h'1. Results for absolute pseudo-range positioning using the precise positioning system (P-code) were 10 to 15 m (lo) for the horizontal coordinates and 20 to 25 m in height when using ionospherically corrected data. When using LI data only, i.e., when not applying the ionospheric correction, the noise level dropped but biases of up to 50 m showed in the results. In the relative pseudo-range mode, which eliminates a large part of the ionospheric propagation effects, results using LI only were usually better giving standard errors of 10 m or better for the horizontal coordinates and 15 m for height. In this case, biases were almost completely eliminated by using the differential technique.
ADVANTAGES
Homogeneous control No targetting or access problems Free choice of flight patterns Instantaneous positioning (repeated flights)
12.18
\ J
/ /
lower density
AIRBORNE GRAVIMETRY
Gravity measurements are one of the major data sources for geodesy and are also an important reconnaissance tool used in geophysics. Terrestrial gravity measurements are done at discrete points. Although very accurate, the measuring procedure is slow and rather difficult to use in remote and mountainous areas. Attempts to make the method airborne have been made, therefore, more than once during the last 30 years but have so far not met with success, at least not in the fixed wing aircraft mode. This failure was due to at least three causes: The difficulty to position the aircraft with an accuracy of 2 m in height, the difficulty to obtain the aircraft velocity to 10 cm s"1, and the difficulty to separate gravitational and inertial accelerations. GPS has the potential to provide solutions to the first two problems.
In order to secure 1 mGal accuracy for the gravity anomalies, the orthometric height of the aircraft must be known to better than 3 m, while the horizontal position is required to within 50 m. The latter requirement can be easily met by any of the absolute GPS methods. The height accuracy is not as easy to achieve. Since the orthometric height is composed of the ellipsoidal height obtained from GPS and the geoidal height, errors in both components must be considered. In general, the geoidal height can be computed with an accuracy of 1 to 3 m (1 a). This means that the ellipsoidal height is required with an accuracy of about 2 m (1 a). This accuracy can be achieved by the relative phase and the relative combined method assuming that cycle slips can be eliminated. If that is not possible, only the relative integrated method will give the necessary accuracy. The additional measurements can be either pressure altimeter data or, again, inertial data. Aircraft velocities have to be known to 10 cm s"1 or better in order to keep the error in the Eotvos correction below 1 mGal. This is achievable with time differenced phase data or through the use of an integrated system. The third problem mentioned above can be rigorously solved by airborne gradiometry. However, a practically sufficient solution is perhaps also possible by airborne gravimetry. It requires the smoothing of high aircraft accelerations and the separation of gravitational and inertial effects by a priori statistical information. The resulting smoothed gravity data would be sufficient for geodetic purposes and may have enough short wavelength information for reconnaissance gravity surveys. Major advantages of the method would be the rapid data acquisition even in remote areas, the coverage of large areas, and the possibility of flying gravity and aeromagnetic surveys on the same mission.
12.19
APPLICATIONS
Laser Profiling
HELICOPTER APPLICATIONS
The major advantage of using helicopters instead of fixed wing aircraft is their greater flexibility. The hovering mode, the capability for quick manoeuvres, and the small space required for landing opens up a whole new range of applications which are not possible with fixed wing aircraft. The hovering mode, for instance, is used to position objects or sites where landing is impossible; examples are positioning of buoys and objects in heavily forested areas. The manoeuvering flexibility makes detailed laser profiling possible in areas of rugged topography. Landing flexibility makes establishment or resurvey of large numbers of control points in forest-free regions very cost effective. Many helicopter-borne inertial surveying techniques currently utilized should be replaceable by GPS techniques once GPS satellites are available on a continuous basis. A series of tests was conducted in 1984-85 to assess the accuracy of helicopters for buoy positioning and Loran-C calibration. A Bell 206B helicopter and two TI-4100 receivers were used for the tests. The GPS antenna was installed on the upper extremity of the vertical tail fin. The tests in 1984 suffered from frequent loss of the GPS signal due to the tracking software used. They indicated, however, that latitude and longitude can be determined with an accuracy of about 10m when using the relative pseudo-range method while hovering. Tests with the improved Phase II tracking software, conducted in 1985, showed clearly that the signal loss problem had been overcome. For control point positioning and laser profiling, accuracies at the half metre level are needed. This requires the use of phase data or an integrated system approach. While the relative phase or the relative combined method will be adequate for the first application, an attitude reference and thus an integrated system is needed in the second application. Test results are not available yet, but extrapolating the results obtained in road tests shows that the potential to achieve these accuracies is certainly there.
CHAPTER THIRTEEN
13.00
DIFFERENTIAL GPS REQUIREMENTS COMMUNICATION LINK REQUIREMENTS What factors must be considered when designing a real-time differential GPS data communication link?
GPS
COMMUNICATION METHODS What signal propagation modes are available? SYSTEM OPTIONS What options are actually available? How do we choose one? SYSTEM COMPONENTS What do we need at the monitor? What do we need at the remote user site?
Canadian GPS Associates, May 1987
The chapter concludes by describing the basic system components of a differential GPS communication link at both the monitor and remote sites.
13.01
GPS REQUIREMENT;
STATIC APPLICATIONS
May use postprocessing Data may be buffered or stored Communication need not be in real time Link equipment may not be power or range limited
KINEMATIC APPLICATIONS
Real-time processing desirable Strict correction message timing, bandwidth, and content Allowed processing time is determined by the vessel dynamics and the control response time Link equipment must be designed to tight power and range specifications Very high data link reliability required
In static applications, the field user is either not moving or is able to remain stationary for as long as desired. This would be the case in determining a precise land position. Under these conditions, the GPS correction may be applied after all the measurements have been taken and recorded (postprocessing) or the correction message may be temporarily stored or buffered for later transmission and processing. For this reason, there is no need for the communication channel to remain in operation at all times. Also, there will usually be more flexibility in the available transmitter power (and, hence, range) for stationary sites.
In kinematic applications, the field user may be an aircraft, ship, or other mobile platform. Postprocessing is not always possible in these cases since the GPS data are being used at least partially for real-time navigation. Under these conditions, the correction message must follow strict length and timing specifications, and the receiver processing time (and message updating interval) are determined by the vessel dynamics and control response time. The communication link equipment must be carefully selected or designed to rigid power and range requirements with emphasis on reliability and redundancy.
13.02
Range
Type of Service
Duty Period
Propagation Constraints
Dynamic Constraints
Cost
Number of Users
Range: What is the operational distance between the monitor station and the user? What is the tradeoff between channel bandwidth, range, and equipment cost? Type of Service: We may select from a variety of communication methods and system configurations, each having its own constraints, advantages, and disadvantages. Duty Period: The communication link may be required to operate continuously, on demand, or only at specified intervals. The duty cycle will be determined by the positioning accuracy requirements at hand. Kinematic applications, such as hydrographic surveys, may require continuous operation, whereas static applications may only require periodic fixing. The duty period may place additional constraints on the communication link data rate. Propagation Constraints: Different methods of radio propagation present different limitations in terms of cost, reliability, range of operation, and equipment complexity. Dynamic Constraints: The degree of movement of the user will have a great influence on the possible equipment choices. Available equipment space, shock loading, receiver tracking capability, and antenna beam coverage are factors which must be considered. The motion of the user will also determine the duty cycle of the communication link. Cost: What is the tradeoff between system complexity, reliability, and unit cost? What is the operational cost of a unit failure? How many units are in the system? Will new technology introduce cost benefits? Will new technology make the present system obsolete, or can the system be expanded easily to take advantage of new components without sacrificing existing units? Service and Maintainability: What are the ambient operating conditions? Will the equipment require specially trained operators or maintenance technicians? What is the expected service lifetime of the units? How and where will units be repaired or replaced? Are the units interchangeable either in whole or in part? Can the monitor be operated unattended? Number of Users: How many users will the system support? Are the user units fully interchangeable at all ranges? Can the system be easily expanded to accommodate additional users? Can two or more such systems be operated side-by-side without interference? Can the fixed monitor station be easily adjusted to serve different networks of users?
13.03
word 1
word 2
message
16 MESSAGE TYPES
1 2 3 4 5 6 7 8 9 10 11 12 13-15 16
Differential corrections Delta differential corrections Station parameters Surveying Constellation health Null frame Beacon health Pseudolite health High rate differential corrections P-code differential corrections C/A-code L1,L2 differential corrections Health message To be decided Special messages
The actual correction messages contain a wide variety of information of use to the surveyor or navigator. These messages relay extensive information on monitor station position and health, satellite constellation health, and corrections to be applied to the direct GPS observations. The most commonly used message will be Type 1 which contains the actual pseudo-range correction. The Type 2 message, for example, will be used to augment the Type 1 data in the event mat new orbital or clock parameters are introduced by the monitor station.
13.04
Preamble (8 bits) Message type which follows (4 bits) Station identification (12 bits) Parity bits (6 bits)
Modified z-count or timing (13 bits) Sequence number (3 bits) Message length in words (5 bits) Station health (3 bits) Parity (6 bits)
13.05
TYFE 1 MESSAGE
Corrections to GPS pseudo-ranges in the form
M ^TT^TIJ
<1
[Fill / Parityj
TYPE 1 MESSAGE
This is the primary message type which provides the pseudo-range corrections. The message contains the correction term, as well as the range rate correction, to enable application of corrections at any time. Timing information for this adjustment is taken from the second word in the message header. The message also contains the satellite identification number for which these corrections apply, satellite health summary, age of the data, and the estimated user range error. The type 1 message is transmitted separately for each GPS satellite.
13.06
TYPE
Corrections to Type 1 data in the form
f Fill /Parity}
Always used to augment a Type 1 message Does not contain satellite ID,heaIth or data age is used only when the broadcast message parameters being used by the monitor station change. It will correct for the difference between the old and new broadcast messages, until such time as the remote has begun using the new set also. The remote uses "age of data" in Type 1 as a flag to use Type 2.
Canadian GPS Associates, April 1986
TYPE 2 MESSAGE
The reference station will transmit a Type 2 message whenever it starts using new GPS navigation data, i.e., when the station switches from one set of orbital and clock parameters to another. If the user does not switch to the new orbital and clock parameters at the same time, the Type 1 data which he has received will be useless to him. Therefore, the Type 2 data gives corections to the Type 1 data for use by those who have not yet switched to the new orbital and clock parameters.
13.07
ionospheric bounce
ground wave
x
remote user
COMMUNICATION METHODS
Many differential applications may allow a choice of the method of communication, and in other cases the choices may be extremely limited. The principal characteristics of several methods are as follows:
Ground wave propagation: At low frequencies, the range may extend to several hundred kilometres. The wave is 'bound' to the earth-air interface and may require large transmitter antennas.
Sky wave propagation: Includes direct and tropospheric refracted waves. At HF frequencies, the range may extend to a few hundred kilometres, and at higher frequencies the range becomes reduced to near horizon. Ionospheric bounce: Service is sporadic and varies from day to night, as well as from season to season. The bounce distances tend to be much larger than the typical differential GPS ranges, and no service is available at short ranges of under a few hundred kilometres.
Meteor trail bounce: Provides burst service at continuous ranges out to several hundred kilometres.
Communication satellite relay: Power limitations of the satellite restrict this usage to remote users who can mount steerable directional antennas. The communication satellite must be geosynchronous so as to provide continuous operation. True mobile satellite communication service is anticipated by the year 2000.
13.08
1YSTEM
Piggyback
Frequency Hopping
Satellite Broadcasting
Direct Dedicated Radio Service
Meteorite Scattering
l l
Pseudolite
SYSTEM OPTIONS
The implementation of the communication channel link allows further considerations as follows.
Piggyback operation: Can the differential correction message be added to some presently used radio service?
Frequency hopping: A technique employed to avoid some propagation limitations. Such systems may be range-limited, and the remote units may not be interchangeable. The system may be complex and expensive but operationally sound in certain applications.
Satellite broadcasting: Not practical for direct mobile service, but may be used as part of a split link.
Direct dedicated radio service: Establishing and licensing a radio service for differential message handling independent of existing radio services.
Meteor scattering: Commercial systems are available but may require a costly transmitting station to serve the region.
Pseudolite service: Provides useful features in some applications (e.g., terminal air navigation), and may be practical as the final leg of a long-range system assuming that interference will not occur with normal GPS operation.
For each of the choices above, we must consider whether the system will be designed to offer real-time or postprocessed correction. Further, the system may be designed to provide direct correction messages to the user, or an intermediate relay link may be used.
^M^
13.09
SYSTEM
GPS Receiver
Computer Link Transmitter
GPS Receiver
Computer Link
Monitor Station
Remote User
SYSTEM COMPONENTS
Monitor Station
The monitor station consists of a GPS receiver, a computer to calculate the pseudo-range and range rate errors and to compose the correction message, a modem, and a radio frequency (RF) transmitter with its associated antenna assembly and coupler. The computer/modem should include the provision for an error correction scheme for the RF link. This is in addition to any error correction or checking schemes which may be built into the GPS message.
Application requirements may dictate that the monitor station be easily transported by, e.g., a helicopter, and operated unmanned. It may also be necessary to select a radio transmitter which can be easily adjusted to serve different user networks on different frequencies.
Remote User
The remote user equipment consists of a GPS receiver, a computer for decoding the correction message and incorporating it with the GPS message, a modem, and a radio link receiver with its associated antenna assembly. The emphasis of the communication equipment selection criteria should be on reliability, simplicity of operation, and servicability. Ideally, all user units should be fully interchangeable. Unit cost for the communication equipment should be low enough to use spares to guarantee backup, particularly where equipment failure would involve very costly mission failure.
Application requirements may place further restrictions on equipment choices in order to meet ambient and dynamic performance standards (e.g., marine spray, arctic, airborne, small launch).
CHAPTER FOURTEEN
14.00
1URVEY BE
GEOMETRICAL DESIGN
SELECTION OF OBSERVING ACCURACY
LOGISTIC DESIGN
OPERATION ISSING SOFTWARE
OUTAGES CONSTRAINTS
>ESIGN IN
PARAMETER
14.01
SURVEY BE;
PARAMETERS
Station geometry
Satellite geometry
Receivers
Time considerations
number of stations configuration of stations baselines number of satellites configuration of satellites mask angle type(s) of receiver to be used number of receivers available configurations for each session weather observations to be made ? manual, remote control, or preprogrammed ? time of day (ionospheric effects) observing window and duration sampling rate (preprocessing) code (C/A or P) or carrier or both LI or L2 or both datum ranges or range differences network or baseline by baseline broadcast or precise ephemerides bias modelling number of personnel method of transportation leased or owned equipment
Whereas there are infinitely many choices of configurations, both of the ground (receiver) stations as well as satellites, and of observing times, for practical reasons, only a finite number of many combinations of configurations, observing times, and other design parameters can be tried. The design is then reduced to a trial and error process. Once a selection of a set of design parameters is made, elements of the N matrix can be evaluated with relative ease. Under those circumstances, the criterion of max Tr(N) instead of min TrCN"1) would lead to the same optimal selection of design parameters and, therefore, to the same design.
In the rest of this chapter, we shall treat three particular subproblems to the above general problem of design. We shall show the choices available for controlling the accuracy of ranging/range differencing. We shall discuss the concepts involved in the ground station configuration design. We shall describe the principles used in selecting optimal satellite configurations.
14.02
Centroid
DATUM DEFINITION
When constructing the design matrix for the survey network it is necessary to describe how to position and orient the network on the earth's surface. This procedure is often known as datum definition. (Note that this is a looser use of the term datum than its mathematical definition, which is that a datum is a particular coordinate surface.) This can be done in several ways: holding one point fixed (placing an infinite weight on the known position of the point), holding one point properly weighted (placing the appropriate weight Px on the one known position), holding no point fixed and using a generalized inverse to invert the (singular) N, holding several points fixed, holding several points properly weighted by Px. The first three approaches use the minimum necessary number of positional constraints and are therefore known as minimum constraint adjustments. The last two belong to a family of overconstrained adjustments. Other techniques also exist.
All minimum constraint adjustments yield identical results for the network configuration (relative positions among the points). It is said that network configuration is invariant with respect to the choice of the minimal set of constraints. The placement of the network on the earth's surface, though, would generally be different. Also, the relative confidence ellipsoids will be the same for all minimum constraint adjustments, whereas absolute confidence ellipsoids will vary. Different situations are shown in the figure in a two-dimensional fashion; the concepts are analogous, however, in three dimensions.
It should be emphasized that the only statistically rigorous way of selecting the datum is to take all the known positions into consideration by putting appropriate weights on them. If a position of only one point is known this includes the case when only a very approximate position (good to only one kilometre) of the network is available and if we are not particularly interested in learning where exactly the network is located, in an absolute sense, then any minimum constraint adjustment will do the job. This should be borne in mind when realizing that GPS determined relative positions are inherently superior to the GPS point positions. The concepts involved in the foothold technique (used, for example, in the Active Control System), will be treated separately in a subsequent section.
14.03
Strong configuration
Weaker configuration
No geometrical redundancy
Canadian GPS Associates, April 1986
GEOMETRICAL DESIGN
Often the survey design problem to be solved includes the geometry of both the satellites to be used and the network of ground stations. The aim then is to design the geometrically strongest combined configuration. Clearly, if time and expense are of no importance and receivers capable of tracking any number of satellites are available, then the use of all visible satellites and all possible baselines that can be created within the network will result in the strongest configuration and will give the best accuracies of relative positions. This is, however, never the case and limitations must be considered. The number of satellites to be tracked simultaneously is limited by the number of receiver channels, real or simulated. Most receivers are capable of tracking four satellites simultaneously, and the selection of the best four satellites is thus the most often posed problem.
Baseline length and orientation affect the accuracy of estimated baseline components relative positions. The selection of baselines within the network can either reduce or enhance geometrical redundancy, i.e., the number of conditions that sums of baselines must satisfy to make a consistent collection (so that all loops made of individual baselines close to zero).
To reduce biases introduced by incorrect modelling of tropospheric delays, the mask angle (vertical angle cutoff) during tracking should be kept reasonably large (8 to 15). Observed values of weather parameters are likely to reduce further suspected biases in tropospheric modelling. If the utmost accuracy is desired, then water vapour radiometers may have to be used. All these precautions are, of course, much more important for long baselines where range differencing is less likely to difference away biases present in the observations.
Selection of the appropriate time of day for the observations reduces the effect of the ionosphere, and keeping the antenna away from reflective structures reduces the danger of multipath problems. Finally, when carrier ranging is employed, potential accuracy decreases with every new range unknown that has to be solved for. Hence, switching to a marginally more favourable set of satellites for a brief period of time would not be a good strategy because it would result in decreased accuracy. Reasonable and practical combinations of choices should be tested in constructing the N matrices. We note that the selection of observing period will have the most significant effect on the internal structure of N. The rest of the variables will affect only the scale aQ2 of N.
ATELLITE SELECTION
14.05
Good configuration
Bad configuration
b dr = p e db
formulated for a single baseline b, where e is the mean unit vector in the direction of the satellite. For several satellites we get
either small or tend to be perpendicular to Zj p. Cj, which in turn, for the optimal satellite configuration, will tend to be directed perpendicularly upwards from the baseline.
Needless to say,we have to make sure that all the satellites in the selected configuration are visible simultaneously to all the tracking receivers.
E DESIGN
TT^vTF^ (TN 17/^*^3 TivvT
14.06
Step 1 Session 1
Step 2 Session 2
Step 3 Session 2
Step 4 Session 3
Step 5 Session 3
Step 6 Session 4
where the economy comes from the fact that AN will have relatively few non-zero elements. For instance, removal of one baseline will affect only two points resulting in AN having only a 6 by 6 non-zero submatrix imbedded in a zero matrix. The addition of a baseline renders a similar result. The use of an incremental matrix AN cuts down the number of operations necessary to evaluate the individual elements of N. Since each observation gives an independent contribution to all the affected elements of N, the number of operations needed to assemble each element may be quite large. Denoting by E the number of contemplated sessions (receiver location permutations), by N the number of observations during each session, by S the number of tracked satellites, and by R the number of simultaneously operating receivers, we have to deal with the total of either ENS(R-l) independent single differences, or EN(S-1)(R-1) independent double differences, or E(N-1)(S-1)(R-1) independent triple differences. This represents, for instance, 10 800 observed double differences even for a relatively modest case of 20 sessions of 60 observations each using four satellites and four receivers.
When the stepwise updating of N is used, then for the deletion of one baseline only EN(S-l) observations of double differences have to be considered in assembling AN.
This figure is based on the use of three receivers to position six points. We start with a three-station setup (step 1); try two alternative expansions to a fourth station (steps 2 and 3), and decide that die former is better; try two alternative expansions to a fifth station (steps 4 and 5), and decide that the latter is better; and finally expand to the sixth station.
14.07
Session 2
Session 3
Session 1
Session 3
From the practical point of view, it is usually desirable to minimize the cost of field work while adhering to the prescribed accuracy. The components in this minimization problem are all or a selection of the following:
personnel transportation means and their deployment number and type of receivers use of own or rented receivers use of remotely controlled or preprogrammed receivers.
The cost associated with the individual components generally varies from company to company, from place to place, and from job to job. It is, therefore, impossible to give a generally valid prescription here, and a solution should be sought for each case separately.
Mathematically, it is possible to work out the exact solution using linear and quadratic programming. Software for these so-called economic problems is commercially available. Most of the time, an approximate solution works out satisfactorily.
The choice of economically appropriate software for processing the collected observations should be considered in the logistic design.
14.08
CASE
Master station
FOOTHOLD OPERATION
In the foothold mode of operation, when a permanent master station is operated and maintained by a service organization, the survey design would have a different character. The service organization may or may not be equipped to design the survey procedures for the user. If not, then it would be up to the user/surveyor to do it himself.
The foothold mode can be used with either one or several roving receivers. In either case, the master station serves as the datum point whose coordinates would be weighted relatively highly, because the accuracy of its position would be quite high.
The distance of roving receiver(s) from the master station would normally be much longer than distances between roving receivers. As a result, the integrity of the relative positions of roving receivers will be quite high, whereas their location with respect to the master station (in metric units) would be less well known.
If the user must carry out the design by himself, the steps to be followed are analogous to those described earlier.
14.09
Not restricted to GPS observations Solutions for station coordinates, satellite orbit, force model and nuisance parameters Mainframe computer
BUILT-IN RECEIVER SOFTWARE
Dedicated to operate on output of a particular receiver type Baseline and/or network mode Personal computer
3rd PARTY SOFTWARE Not restricted to particular receiver type Baseline and/or network mode Mainframe, mini, or personal computer
PROCESSING SOFTWARE
Even though with proper formulation the results are invariant with respect to the processing mode (ranges, range differences, double differences, triple differences, etc.), there are still a number of other features that deserve a thought from the point of view of selecting the appropriate software for processing the collected data. First, the available software varies considerably as far as bias modelling is concerned: orbital biases, clock biases, range biases, atmospheric biases. Secondly, not all software is able to process data in the network mode. If only individual baselines can be processed, then a separate program has to be used subsequently to adjust the baselines within the network. Modelling of tropospheric and ionospheric delays may also vary significantly. Some software allows input of actual weather parameters; others do not. Data rate reduction (normal point formation) is treated differently in different programs and so are cycle slips. Lastly, different pieces of software provide different kinds of output, particularly in the error estimates. Not all software is available for all categories of computers. If only a personal computer is available, then the choices may be fairly limited. Also computer time consumption may have to be taken into account, since the performance of different programs varies appreciably.
OUTAGES
The geometry provided by GPS satellites changes continually, as the satellites move in their orbits, and as the user moves from place to place over the surface of the earth. For certain times and at certain places the geometry can become very bad. Such cases are referred to as outages. An outage is the occurrence (in time and space) of GPS satellite geometry so bad that the GPS position accuracy is worse than some predefined level. Full information about the accuracy of a GPS position is contained in the covariance matrix associated with that position. It is often useful to separate the accuracy estimate into two factors: that due to range measurement errors, and that due to geometry, as represented in a simple way by the Dilution of Precision (DOP). This single number describes the geometrical strength of the GPS satellite constellation at any point in time and space, and is defined as the factor by which GPS range errors should be multiplied to yield the corresponding positioning error (see Chapter 4). Outages are usually defined in terms of DOP values. For example, we may specify that an outage occurs whenever GDOP exceeds 6 that is, whenever the r.m.s. solution error exceeds six times the range measurement error. For range errors of, say, 2 metres, this would yield an r.m.s. solution error of 12 metres. If we take the full covariance matrix information into account, we may find each of the position coordinates and time to have to quite different accuracies. But as a simple measure of the effect of geometry, the DOP value is quite useful.
Once the full GPS constellation is in place, outages will still occur, but briefly (typically for less than 10 minutes, once per day), and only in specified locations on the surface of the earth. This figure shows the occurrence in time and space of outages near Baffin Island for the 18 + 3 constellation. In this case, the outage is specified in terms of PDOP. Outages depend strongly on mask angle: larger mask angles will cause more serious outages, both in extent and duration. For the mask angle of 5, considered in this figure, outages last typically for about 15 minutes, and are of limited extent. At a mask angle of 20, this outage expands to cover most of North America. Outages are of primary concern in navigation and real-time point positioning, where during an outage other positioning systems must be used. In surveying, they are a minor nuisance, causing only a forced extension of time needed for field observations thus making the GPS survey more expensive.
This figure was prepared with the help of program MacGEPSAL, available from the Geodetic Research Laboratory, Department of Surveying Engineering, University of New Brunswick.
,FFEGT OF HEIGHT
1 onr^Or^ /f*^f*\T^ /f^Tr^
TT TTTr^TT/^"^ TT TT^Tn
14.11
HDOP at Halifax, Nova Scotia New Year's Day 1990 Mask angle 10 No constraints
24 Hour of day
5CH O
HDOP at Halifax, Nova Scotia" New Year's Day 1990 Mask angle 10 Height constraint of 5 metres.
Q ffl
0 0
Canadian GPS Associates, April 1986
12
i i 18 24 Hour of day
Investigating the effect of various constraints on expected results for the remaining parameters is an important aspect of GPS survey design. Since we will never know our position coordinates or time perfectly ahead of time, it is important that we use a technique known as weighted constraints, discussed earlier. We still solve for the constrained parameters, but we let our solution vary from the assumed 'known' value only slightly. For example, if we assume that we know the height to 5 metres, we can instruct the GPS software to prevent the height value in the solution from varying more than 5 metres from the 'known' value we specify for height going into the solution.
What do we gain from a constrained solution? Improved geometrical strength, and enhanced redundancy, and hence improved accuracy for the solution for the remaining, unconstrained, parameters. The tighter the constraint (that is, the more accurately we know the constrained parameters) the greater this improvement in geometrical strength. What do we risk by using a constrained solution? This depends on whether we actually know the constrained parameters as well as we have told the algorithm we do. For example, if we constrain the height to 5 metres, but use a 'known' value for height which is wrong by 50 metres, we will have a much worse horizontal position than if we had assumed we knew nothing about the height and had solved for it. How much worse depends on the geometry again, but would typically be of the same order as our height error50 metres. The trick is to match the constraint to the accuracy with which we know the constrained parameter. This figure shows the results of a GPS hydrographic survey planning exercise in which the application of height constraints to improve GPS geometry was investigated. The top plot in this figure shows the variation in HDOP at Halifax, Nova Scotia, on New Year's Day 1990, assuming the 18-satellite constellation and a 10 mask angle. Two outages are evident, where HDOP approaches a value of 10. The bottom plot shows the result for the same case, but with a 5 metre height constraint applied. This constraint eliminates the outage. We conclude that, if the shipboard GPS receiver ellipsoid height can be predetermined to 5 metres (requiring accurate geoid maps, estimates of tide and swell, height of antenna above the water level, and the effects of roll and pitch), then constraining it will improve the geometrical strength enough to eliminate outages.
This figure was prepared with the help of program MacGEPSAL, available from the Geodetic Research Laboratory, Department of Surveying Engineering, University of New Brunswick.
CHAPTER FIFTEEN
L Ai
15.00
STANDARD;
GPS survey standards Field procedure standards Data processing standards
Perhaps the next few pages, although not providing an exhaustive description of GPS surveying, will at least provide a sketch of the procedures to be followed and some of the criteria which may soon make up the GPS surveying standards.
15.01
PLANNING
DESIGN
PLANNING
Until the cost of GPS receivers drops and a full satellite constellation is available, careful planning of GPS surveys is critical for a successful project [HothemetaL, 1984]. The initial planning stage involves the project layout and network design which are primarily dependent on accuracy requirements, satellite constellation, and geometry (as already discussed in a previous chapter). Consideration must be given to the number of receivers, redundancy, time per station, number of stations per day, observation window, and whether one station will be used as a base for the whole project or a leap frog approach will be used.
Accuracy requirements will determine whether or not an L1/L2 receiver is required and whether or not precise ephemerides need be arranged for. Choice of receiver will affect the type of recording system and media which will be required. Codeless receivers will require that receivers are physically synchronized often during a project. Other receiver characteristics which affect logistics are whether power requirements are AC, 12 or 24 v DC, whether the equipment is portable and ruggedized or merely transportable, and what antenna cable lengths are available.
Terrain and station accessibility will determine whether car, truck, or helicopter is used to transport the receivers and the number of personnel required.
Planning cannot be finalized until site reconnaissance is completed. Sites must be inspected for vehicle access, cable length requirement, obstructions, foliage and, if necessary, the monument itself must be found. While receivers are expensive and the operating window is short, few operations can afford surprises during the observation campaign.
15.02
FIEL1
BEFORE LEAVING BASE Receiver synchronization (if required) EN ROUTE TO OBSERVING STATIONS Maintain power to receiver SETTING UP AT OBSERVING SITES Antenna - set-up, centre, measure phase centre height Receiver - initialize, select satellites (if required) Recorder - ensure data is being recorded Perform any required calibrations OPERATION AT OBSERVING SITE
Monitor receiver operation Collect temperature, humidity, and pressure data Enter the following in the log book Weather Receiver problems Operator errors Obstruction diagram
COMPLETION OF OBSERVING SESSION
Remove and label data tape (disk, etc.) Repeat calibration checks (if possible) Next site or return to field office
Canadian GPS Associates, April 1986
The daily operations in the field office, wherever it may be, include [Beck, 1985]: synchronization of receivers when codeless techniques are used, charging of batteries, fuelling generators and vehicles, and initializing data storage media. Before heading out to the stations, receivers are often placed in standby mode to warm up the internal oscillator. Once at the site, the antenna is centred above the station on a tripod and the height to the antenna phase centre is measured. All cables are connected and the receivers are initialized so that visible satellites are acquired. For receivers that do not track all satellites in view, it is important that all receivers used in the project are simultaneously tracking the same satellites. When tracking has begun, the operator must ensure that the recording device is functioning properly and that both measurements and broadcast ephemeris (if applicable) are recorded. The operator then monitors tracking performance by watching receiver signal quality indicators. Switching satellites and making pressure, temperature, and relative humidity measurements may also be required. Log books are maintained to record any operator errors, receiver problems, tracking problems, obstruction diagrams, and weather data.
When the predetermined amount of data has been collected, the antenna height and centring is checked, the antenna taken down, and the equipment moved to the next site or back to the field office or camp.
15.03
mo
IN THE FIELD
Data verification
PREPROCESSING
Smoothing Editing (cycle slips, etc.) Ambiguity resolution (first step) Ephemens representation
ADJUSTMENT
Ambiguity resolution (second step) Baseline relative positioning computation Three-dimensional network adjustment
DATA PROCESSING
The first stage of data processing is transferring the data to (usually) diskette for backup and/or further processing. There is a trend toward using transportable microcomputers for the data processing which allows processing for most applications to be completed in the field (where no external ephemeris source is needed). Re-observation can be scheduled before demobilizing as required by the processing output. Data is verified to ensure data quantity and quality. Preprocessing also includes measurement editing and smoothing, such as correction or flagging of cycle slips and differencing of phase data between stations, satellites, and epochs as required. The collected broadcast ephemeris must be decoded and interpreted; precise ephemeris can be obtained from the Canadian Geodetic Survey or the U.S. National Geodetic Service but at best two weeks later. In either case, the ephemeris representation is usually altered by, e.g., fitting a curve through consecutive broadcast segments to form a contiguous representation or by computing satellite coordinates for each satellite at each transmission epoch. One approach to the adjustment stage of the processing is as follows. Accounting for horizontal and vertical antenna offsets and weather data, three-dimensional geometric relative positions using non-differenced or differenced phase and range observations are computed. The processing scheme and software choice can vary depending on accuracy objectives. The results of this single baseline, or session (multi-baseline) or network (multi-session) processing, consisting of coordinate differences and associated covariance matrices are then input into a minimally constrained three-dimensional network adjustment program in order to generate reliable statistics and detect any remaining errors. An alternative adjustment approach is to input the preprocessed observations from a multi-session network campaign directly into a network adjustment package. This approach forces the adjusted coordinates to be consistent over all observations in the campaign, rather than only on a session by session basis as above. In this way, biases which are strongly deterministic with one session but to some extent random between sessions (e.g., orbital biases and ionospheric biases) will be able to contribute to the estimation of a more realistic network covariance matrix. This latter approach is taken by DIPOP, the third-generation differential GPS positioning program package developed at the University of New Brunswick for the Geodetic Survey of Canada. Data analysis includes checking whether the internal consistency of results and the results of repeated baselines meet the desired accuracy. Residuals are examined and the variance factor tested to ensure proper weighting The final processing step is to present the results, and to archive all text, diskettes, tapes, etc., in an orderly manner.
&
wK
11
1 . 4 50
I 1 1 |
I 1 1
NEW SPECIFICATIONS MUST ADDRESS 1 Various accuracy requirements 1 New parameters 1 Receiver and software calibration 1
1 1
No published government standards for GPS surveys formally exist yet. New parameters will require definition in these forthcoming standards. Some of the issues are [Hothem, 1985]:
GPS provides three-dimensional capabilities requiring ties to bench marks satellite geometry and ephemeris accuracy constraints need definition geoid slopes and accuracies need definition different measurement types and accuracies are available baseline, session, or network solutions are possible redundancy can come from the number of occupations per stations, repeat baseline observations, measuring antenna height twice.
Even when standards are defined, two major problems will remain: instrumentation calibration/compatibility, and software calibration. To compare results of instrumentation or combination of receivers and to evaluate realism of statistics and results generated by GPS processing software and to ensure compatibility and confidence in results generated from different software packages, three-dimensional calibration networks akin to current EDM baselines are being established, and test data sets are being provided for software calibration.
At least 100 GPS survey sets are currently being used in Canada and the U.S. alone. The rapidly increasing number of receivers is imposing a need to establish surveying standards for both field and office procedures. The following pages address some of these but are by no means near complete. GPS surveying standards will evolve as experience with GPS grows.
FIELD
HOUR:
15.05
INSTRUMENT CALIBRATION
Number of horizontal control points Number of bench marks Minimum number of receivers Number of occupations per station (redundancy) Minimum & maximum station separations Connections between stations Observing period per point LI only, or LI & L2 Meteorological data acquisition
Very preliminary values, for the above parameters as a function of accuracy, can be obtained from USNGS for use as a guideline.
15.06
ING
Ephemeris source and accuracy Measurement rejection criteria Maximum acceptable measurement rejection rate
POSTPROCESSING ANALYSIS STANDARDS
After computing vectors or baselines, any baselines observed more than once would have to meet repeatability criteria for each of the three-dimensional vector components as a function of baseline length. A minimally constrained three-dimensional adjustment would reveal whether realistic observation weights were used and what the standard errors (precision) of the three-dimensional vector component are. A constrained adjustment using existing control coordinates accounting for scale, orientation, and undulation relationships between the GPS satellite and the control network datums would provide proper integration of the GPS relative positioning results (obviously at a level corresponding to the control network accuracy). Specifications for the display and storage of results and variance-covariance information will also be required.
15.07
ESTABLISHING CONTROL POINTS AS REQUIRED, USING GPS RELATIVE POSITIONING TECHNIQUES, IS BECOMING MORE COST EFFECTIVE
THAN MAINTAINING MONUMENTATION
EVENTUALLY
MONUMENTATION
The need for stable expensive long term monumentation at often difficult to access sites might soon disappear except in urban areas where GPS satellites are not readily visible and for very high precision geodynamic studies networks. With GPS, re-establishing points is becoming less expensive then monumenting. Because there is no intervisibility requirement, new points established should be accessible by land transport and should be placed where they will be most useful.
Azimuth reference marks may need be added for local survey orientation. These can be established using GPS, from star observations using a theodolite, or by observations using a gyrotheodolite.
Continuously tracking GPS monitor stations, as proposed by the Texas State Department of Highways and the Canadian Geodetic Survey, could satisfy most monumentation requirements and provide additional benefits to the surveyor [Steevesetal., 1986]: no physical access to control points required no descriptions of control points required may be able to use fewer receivers automatic access to latest control point coordinates single datum automatic integration of information into the national framework.
APPENDIX A
A-l
i.iL^vt=xlLijJl v IL^rJ
Ff F1FVF1
APPLICATION Real time or postmission Static or dynamic Differential or point positioning Accuracy requirement TYPE
One or two frequencies Code dependant or codeless Multichannel, switching, hybrid C/A- or P-code Fully integrated system or sensor only
FEATURES Carrier phase Data logging Physical characteristics and portability Power requirements Software
the application for which it will be used; the accuracy requirement; the type of signal processing required; the necessary features and options required to suit the application; and, of course, cost.
The more applications a receiver will fulfill, the greater its flexibility, the fancier the features and options, the greater the cost. Real-time applications may only require displayed positions, velocities, and waypoint navigation, whereas postmission applications obviously require data logging devices and interface ports. Static receivers for point positioning or time transfer do not require a fast sequential or multi-channel receiver configuration, whereas one mounted in a helicopter might. Defining the accuracy requirement to 50 m, 15 m, 5 m or submetre determines whether differential techniques are required, C/A- or L1/L2 P-code is needed, and whether carrier beat phase measurements must be available. The environment in which the receiver will be used dictates the antenna type and mount, receiver size, portability, and degree of ruggedness. Sometimes only a GPS receiver sensor, as part of a multi-sensor integrated navigation system, is required, as opposed to a complete self-contained receiver system. GPS receiver technology is rapidly evolving. What follows is essentially a catalogue of commercially available GPS receiver products. The information has been gleaned from company brochures, manuals, phone conversations, and proceedings. (Where no reference is given, information was obtained from brochures.) Due to the dissimilar type and quantity of information available, the discussion format for each receiver varies.
DISCLAIMER: An attempt has been made to include most known, commercially available GPS receivers; however, some may have been overlooked or insufficient data was available. Trademarks used belong to the respective manufacturers. Specifications, prices, availability, and material which follows are for information only and are subject to change without notice. For specific and current information, GPS receiver manufacturers should be contacted. No endorsement is hereby intended or implied.
A-2
MANUFACTURER:
2 channel LI, C/A-code, multiplex 1 channel tracks up to 8 satellites pseudo-range and pseudo-range rate 2nd channel to acquire satellites, collect almanac and ephemeris data 1/2 ATR short Standard interface ports for external sensor aiding and relative positioning 3-D navigation function Optional low profile antenna (0.8 x 13 x 13 cm)
AVAILABILITY:
1986 $ to be determined
Canadian GPS Associates, May 1987
CMA-786 (C-SET)
Manufacturer: Canadian Marconi Company 415LeggettDrive P.O. Box 13330 Kanata, Ontario Canada K2K2B2
Canadian Marconi is taking advantage of experience gained in producing the military CMA-774 (D-Set) P- and C/A-code L1/L2 medium dynamics receiver in its development of a commercial GPS receiver initially intended for air navigation applications [Canadian Marconi Company, 1985]. What makes this navigation receiver different from any other is its dual channel C/A-code LI multiplexing architecture. One of the channels is dedicated to tracking all of the up to 8 acquired satellites to generate pseudo-ranges and Doppler measurements. The second channel performs utility functions including: satellite acquisition (and re-acquisition as necessary); data bit synchronization and almanac and ephemeris data collection.
The advantages of tracking all visible satellites, as opposed to the minimum of four required for three-dimensional positioning, are: no need to acquire new satellites during satellite shading; no requirement for GDOP calculation; redundant measurements and therefore improvement in navigation accuracy.
In its final design, the CMA-786 will accept (optionally) external aiding through various standard interfaces (ARINC 426 I/O, RS-232C and RS-422) from sensors such as altimeter, Omega, and Doppler and a data link to support differential mode of operation. Operational features include automatic signal acquisition and tracking and typical three-dimensional waypoint navigation functions.
Physically, the CMA-786 is housed in a 1/2 ATR short unit weighing 8 kg, drawing 50 W maximum from 24 VDC or 115 VAC. Full temperature range components are optional. The CDU can be an off-the-shelf CMA-823. Two choices of antenna are available; either a 4 cm diameter by 13 cm helix cylinder type for ship or helicopter, or a crossed-slot low profile 8 x l 3 x l 3 c m antenna for commercial aircraft.
A-3 I
MANUFACTURER:
4 satellite sequencing on 1 LI, C/A-code channel Single piece (30 x 40 x 17 cm) ruggedized receiver/processor, CDU, and thermal printer Operates at -10C to +50C 11-15 VDC (1 hour internal battery) No ephemeris or measurement data output 30 waypoint navigation Compact (6 x 35 cm) antenna Usable in relative positioning mode (static) 30 m single point mode (r.m.s.) 5 m relative positioning static (after 5 min.) 0.1 m/s (r.m.s.) speed
AVAILABILITY:
JMR SatTrak
Manufacturer: EDO Canada Ltd. 8-6320 llth Street S.E. Calgary, Alberta Canada T2H2L7
As of early 1986, EDO Canada Ltd. had available a portable GPS positioning instrument transportable in a single package (including antenna) of 30 x 40 x 17 cm weighing about 10 kg. This LI C/A-code single channel sequential receiver is for low dynamic and static positioning applications. The system functions on its own internal 12 VDC battery or can be connected to an external power supply.
The displayed (or printed) output consists primarily of two-dimensional or three-dimensional positions, speeds, waypoint information and time. The best four satellites to track are automatically selected.
In recent tests, 30 m r.m.s. accuracies for both static and dynamic positioning were achieved [Ayers and Gonthier, 1986]. By taking 5 minute averages at two static sites and performing a translocation, three-dimensional relative positioning accuracies of 5 m r.m.s. were achieved.
A-4
1TROLAB II
APPLICATION:
Low dynamics navigation, time transfer, research and development training set
1 channel 4 satellite sequencing LI, C/A-code Interfaces DECCA, ARGO, Syledis, Mini-Ranger, Speed Log and Gyro 48 x 14 x 44 cm (rack mounted) 18 kg 19-32 VDC, 110-220 VAC ( 30 V) Waypoint navigation Relative positioning capability (option) Resident user software
AVAILABILITY:
90 days delivery
Canadian GPS Associates, April 1986
ASTROLAB II
Manufacturer:
Interstate Electronics Corporation (IEC) 1001 E. Ball Road Anaheim, California 92805 U.S.A.
IEC has had many years of GPS experience, primarily for the military in support of the Trident I and II programs. Their experience won them the contract for development of GPS hardware for a new range application program, involving differential GPS and integrated low cost INS.
In 1984, IEC had available commercially the model 4200, a C/A-code LI time transfer and positioning receiver. Since that time they have acquired additional GPS technology and now have available the Astrolab II, again a C/A-code LI sequential receiver for time transfer, low dynamics navigation and relative positioning. This rack mount non-ruggedized unit is intended as a training set for users wishing to learn about GPS through experience. The system has both IEEE 488 and RS-423 interfaces which allow interfacing to print and plot devices along with surface radio-location systems, such as DECCA, ARGO, Syledis, and Mini-Ranger, with available software packages. It also has an optional speed-log and gyro interface. Software also allows for relative positioning, controlling input/output to radio links. Several waypoint navigation functions are included.
The 48 x 14 x 44 cm, 18 kg receiver is tolerant of a wide range of power inputs: 19-32 VDC or 110, 220 V (30 V) AC.
A-5
MANUFACTURER: APPLICATION:
FEATURES: LI, all visible satellites, codeless pseudo-range Consists of 3 pieces Antenna/receiver (7 kg) Rubidium clock interface unit with internal power (18 kg) Data analyser/recorder/computer (7 kg) External ephemeris source required Not restricted to GPS satellites Accuracy: 10 cm in 15 min 30 cm in 3 min
AVAILABILITY: 45 days delivery $US 57 000 per unit complete with all processing software L1/L2 and carrier phase processing options may be available in 1987
Canadian GPS Associates, May 1987
SERIES is an acronym for Satellite Emission Range Inferred Earth Surveying. The technology originated at the Jet Propulsion Laboratory (JPL) and evolved to a configuration of Modeless' pseudo-ranging in 1980 [MacDoran et al. 1985a], which will not be affected by the U.S. Department of Defense (USDoD) selective availability policies. The first commercial version of the product, the Model 1991 Land Surveyor marketed by GEO/HYDRO of Rockville, Maryland, was introduced in 1985. Since then, a more portable and efficient version, the GPS Positioner Model 2002, has been introduced. Both systems are best used in relative positioning modes. The Model 2002, in its three piece portable configuration of a total of 32 kg, can achieve accuracies at the 30 to 10 cm level in 3 to 15 minutes, respectively (with good geometry). Like other Modeless' receivers, an external ephemeris is required. ISTAC makes orbit determination software available along with an ephemeris service. The unique features of the ISTAC products are that they track all visible GPS satellites without knowledge of their orbit, and that the receiver could be tuned to track other satellites, such as GLONASS satellites. The system configuration is different from most GPS receivers, consisting of the antenna/receiver assembly (23 x 42 cm), connected to the clock interface unit which houses the rubidium atomic clock and timing electronics connected to the field portable computer/data analyser/recorder (46 x 33 x 13 cm) for data collection and verification. This computer can also be used in the field for baseline postprocessing using ISTAC supplied software.
MANUFACTURER: APPLICATIONS:
ISTAC, Inc. Low dynamics relative positioning for marine environment (e.g., seismic surveys, inland harbour navigation)
FEATURES:
Under development
As Model 2002 but L1/L2 for marine application < 10 knots Communication cross-link required 1 second update interval Accuracy: 5 m @ 500 km 10 m @ 1000 km
AVAILABILITY:
Using the same receiver technology as employed in the Model 2002, the two frequency SERIES marine positioning sensor (MPS-1) is under development for real-time marine relative positioning applications at speeds of < 10 knots [MacDoran et al., 1984; MacDoran et al., 1985b]. In order to achieve real-time position, communication cross-links must be supplied between the receivers. Accuracies of 5 m to 10 m are expected at receiver separation distances of 500 to 1000 km, respectively. Applications include offshore survey vessel location, seismic vessel and streamer location, and harbour navigation.
JLRGPS NAVIGATOR
MANUFACTURER:
A-7
APPLICATION:
FEATURES: C/A-code 1 channel LI, sequential receiver 30 waypoint memory Cross-track error, distance, bearing, time to go, etc. No raw data available Outputs positions, time, velocities 1 second cycle rate Antenna: 35 x 6 cm; < 1 kg; -25C to +70C Receiver/Processor/Display: 32 x 25 x 12 cm; 6 kg; 0C to +50C; 110/115/220 VAC or 24 VDC, < 30 watts Optional interface to remote display, printer, plotter Accuracy: (as at April 1986) 30 m r.m.s. position 0.1 m/sec r.m.s. speed
AVAILABILITY:
SCAN 32 ooo
Canadian GPS Associates, May 1987
JRC has been developing a GPS receiver since 1979. They completed the prototype JLR-4000 in 1982 and announced its availability in 1984. This sequential C/A-code, LI, compact, single channel receiver is designed solely for low dynamics navigation of, e.g., marine vessels and automobiles [Japan Radio Corporation, 1985; Tomioka and Yamada, n.d.]. It is designed to provide 30 m (r.m.s.) position accuracy and 0.1 m/s r.m.s. velocity accuracy each second with the current GPS constellation. This information is available to the control display unit of the receiver and optionally to a remote display, colour plotter and printer. Navigation functions, such as 30 waypoint storage, distance, bearing, time to go, cross track error, arrival alarm, etc., are included. Because it is solely for navigation, broadcast message data and measurement data are not available for output. The compact receiver/processor/display unit is powered from external 10-40 VDC or 110/115/220/230 VAC consuming on average less than 30 watts. It is operational in non-hostile environments at temperatures of 0C to 50C. The 6 cm diameter 35 cm quadrifilar helix omnidirectional antenna can withstand temperatures of -25C to +70C. The maximum antenna cable length is 50 m.
In Canada, the JLR-4000 is available through the Canadian Marconi Company, 2442 Trenton Avenue, Montreal, Quebec, Canada H3P 1Y9.
A-8
LT
:ANUFACTURER;
APPLICATION:
Receiver/processor -15C to +70C 115 VAC (400 Hz) Access to data and measurements Interfaces: ARINC 575 to CDU ARINC 429 auxiliary RS-232 for testing, etc. Marine/land applications options Ruggedized antenna/preamp assembly Various power supply options Differential GPS interface
AVAILABILITY:
LTN-700
Manufacturer: Litton Aero Products 6101 Condor Drive Moorpark, California 93021 U.S.A.
The Litton Aero Products LTN-700 GPS Navigation Set, designed primarily for the commercial aviation market, has been available since 1984. Although it is a GPS sensor, it is usually sold as a system consisting of the LI C/A-code, 4 satellite fast sequencing receiver/processor, control display unit and antenna. What makes this receiver unique is how it sequences through the satellite signals [Sturza, 1984]. During acquisition mode, the receiver sequentially acquires the satellites (through a slow sequencing scheme). Then during tracking a fast sequencing scheme is implementedeach 20 ms is divided into five, 4 ms intervals. Carrier tracking of four satellites and code tracking on one (of the four) satellite takes place every 20 milliseconds. A complete cycle of code and carrier tracking of 4 satellites requires 80 ms, easily allowing the receiver to update navigation outputs every 1 second. This scheme also allows for simultaneous data demodulation.
The system has also been designed to be used as a developmental tool allowing both raw and processed measurements and data to be accessed via RS-232 ports. Optional configurations allow the system to be used in marine and land navigation applications including a differential GPS interface, ruggedized antenna/ preamplifier assembly, and various power supplies.
The next Litton Aero Products GPS air navigation receiver system will incorporate a two-card receiver system to enable * all-in-view' satellite tracking.
A-9
MANUFACTURER:
APPLICATION:
FEATURES: 6 channel LI only Carrier phase squaring (codeless) No broadcast ephemeris decoding Both 12 VDC and 115 VAC (350 watts) required 91 kg without generator 3 mm antenna phase centre accuracy Synchronization of receivers required Ephemeris available through contract Antenna/preamp operable at -50C to +70C AVAILABILITY: No longer for sale Still widely used (services can be leased) Minimum configuration includes: 2 field units 1 processing system and software 1 GOES satellite time receiver 1 time interval counter
Canadian GPS Associates, May 1987
MACROMETER V-1000
Manufacturer:
Litton Aero Service 8100 Westpark Drive P.O. Box 1939 Houston, Texas 77251-1939 U.S.A.
Early versions of Macrometer technology were introduced to the survey community in 1982. Since the January 1983 tests by the (U.S.) Federal Geodetic Control Committee (FGCC), the Macrometer V-1000 has been used extensively for precise relative surveying. In spite of its size (58 x 56 x 64 cm) and weight (73 kg), exclusive of the antenna (18 kg), the system is transportable and has been successfully used in both small 4 wheel vehicles and helicopters. Because the units use codeless carrier phase 'squaring' techniques, no broacast ephemeris data decoding is done and no clock synchronization using the GPS signal is possible. This scheme requires that the receivers be synchronized to each other and also to UTC via, e.g., a GOES (Geostationary Operational Environmental Satellite) satellite time receiver, and that power is maintained continuously to the internal oscillators. Ephemeris for postprocessing is available within 3 to 5 days under contract with Litton Aero Service or within several weeks through the U.S. National Geodetic Survey (USNGS).
The antenna pre-amplifier design is such that the unit can be operated between -50C and +70C The antenna phase centre uncertainty is estimated to be 2 to 3 mm regardless of orientation, and the attached large ground plane is effective in reducing multipath.
The data are recorded in the non-volatile bubble memory unit of the receiver unit and transferred to a small data cartridge after the recording session, which can last from 10 minutes to 6 hours depending on accuracy requirements and geometry. The data are processed postmission with the Aero Service-supplied processing system and software. Macrometer is a registered trademark.
A-10 II
As V-1000 plus L2 frequency 3 portable units per system 27 kg (less antenna) 12 VDC (200 watts) 1 to 2 ppm currently achievable using provided ephemeris 10 min to 6 hrs depending on desired accuracy
AVAILABILITY:
$US 240 000 Approximately 6 months delivery Complete system 2 Macrometer II portable field systems 1 data processing system and software 1 GOES satellite time receiver 1 HP time interval counter
GPS Receivers A. 10
MACROMETER II
Manufacturer:
Litton Aero Service SlOOWestparkDrive P.O. Box 1939 Houston, Texas 77251-1939 U.S.A.
The Macrometer II, conceptually identical with the V-1000, has been redesigned to consist of three portable waterproof units weighing a total of 27 kg. This dual frequency version of Macrometer technology allows for ionospheric calibration and high accuracy relative positioning on long baselines even after the U.S. Department of Defense implementation of denial of P-code access. Operationally, the Macrometer II is used like the V-1000 except an AC generator is no longer required, just a 12 V battery (200 watts) will suffice.
Accuracies of 1 to 2 pjHn have been readily achieved with the Macrometer system, as indicated in the many papers written by a variety of users.
In the future, Litton Aero Service, in conjunction with Litton Aero Products, will be introducing the Mini-Mac, a compact lightweight hybrid C/A-code carrier phase 'squaring' receiver system. Not much else is known about it at this time.
Macrometer II and Mini-Mac are trademarks of Aero Service Division, Western Geophysical Company of America.
A-ll
MANUFACTURER:
Magnavox
GPS receiver card to fit in DEC VT-103 terminal including PDP-11/23 computer Controlled environment operation Receiver physical specs defined by VT-103 Antenna: -40 C to +70C (land or marine configuration available) Position, pseudo-range, ephemeris can be output to an RS-232 port for recording, printing, etc.
AVAILABILITY:
$US 50 000
GPS Receivers A. 11
T-SET
Manufacturer: Magnavox 2829 Maricopa Street Torrance, California 90503 U.S.A.
Magnavox has been involved in the military GPS program for about 20 years. The first GPS receiver which they offered commercially was the Magnavox T-set GPS Navigator. The initial versions were sequencing receivers; now they are a 5 channel continuous tracking LI C/A-code system. The GPS receiver itself is contained on one card made to fit in a DEC VT-103 terminal with a PDP 11/23 computer. Thus, environmental, physical, and electrical characteristics are defined by the VT-103 terminal specifications, which means that the operation must be confined to controlled environments with conventional AC power. The land or sea antenna/preamp units (1.3 kg or 4.5 kg, respectively) operate at -40C to +70C.
Four of the five channels are designed to continuously track satellites, while the fifth is used for inter-channel bias calibration, ephemeris data gathering, and locking onto new satellites. Positions and speeds are output to the display, while processed and measured pseudo-ranges and Doppler and ephemeris data are available for output on RS-232 ports.
A-12 1
MANUFACTURER:
Magnavox
Combined GPS/Transit receiver system 2 channel, 4 satellite, LI C/A-code tracking Sequences through 4 satellites and updates position every 1.6 seconds Physically similar to MX 1100 series navigator Transit receiver systems MX 1100 Transit systems can be retrofit with GPS addition Combination Transit/GPS antenna or two antenna configuration available
AVAILABILITY: New $US 45 000 Retrofit $US 25 000
Canadian GPS Associates, April 1986
The MX 1100 GPS Navigator combines GPS and Transit capability (at least until GPS is fully operational) in the familiar MX 1100 Transit receiver package. The two receivers operate concurrently displaying navigation results, with the system picking the optimal one. On the GPS receiver card, one of the two LI C/A-code channels sequences through each of four satellites every 1.6 seconds, while the other searches for new satellites and gathers data. Positions are updated every 1.6 seconds.
The MX 1100 series navigators can be retrofitted with the additional GPS receiver system. A configuration of either a combination Transit/GPS antenna or two separate antennas is available.
Until the full GPS constellation is available, the combination GPS/Transit receiver provides better accuracy than a Transit-only system, and more continuous navigation than a GPS-only system.
A-13
MANUFACTURER:
Magnavox
APPLICATION;
2 channel, LI C/A-code Output of selectable groups of processed and raw GPS data; 2 RS-422 input/output ports 8-state Kalman filter Automatic best constellation selection Position updated every 1.2 s May be interfaced to Transit or Loran-C receiver and/or heading and speed sensors for dead reckoning Receiver: 10 - 28 VDC (20 W); -20C to +50C; waterproof; 35 x 33 x 11 cm; 7.2 kg; rack-mount option Antenna: -40C to +50C; 33 x 21 cm; 4.5 kg Demonstrated accuracy: 15 m r.m.s static: 5.1 m (half hour avg.) velocity: 0.05 m/s
AVAILABILITY: 1986 SCAN 25 000
Canadian GPS Associates, Apnl 1986
OPS Receivers A. 13
MX 4400
Manufacturer:
The MX 4400 GPS Positioning and Navigation System is Magnavox's dedicated, self-contained GPS positioning receiver. It shares the technology that has been developed for other receivers in the Magnavox GPS product line. The MX 4400 is a 2 channel, LI C/A-code tracking receiver similar to that in the MX 1100. It has been designed to meet the needs of a broad spectrum of scientific, military, and commercial positioning and navigation requirements. The receiver is housed in a compact, rugged, waterproof enclosure and is available in a self-supporting or rack-mounted configuration. The front panel has a two-line back-lit display and a simple, straightforward keypad for receiver control. The receiver can display position updates every 1.2 s. Two RS-422 input/output ports are provided. These may be used to output selectable groups of processed and raw data and in the future to input differential pseudo-range and delta pseudo-range corrections in RTCM 108 signal formats. A 5 MHz signal from an external atomic frequency standard may be supplied for positioning with fewer than four satellites. The MX 4400 has software for waypoint navigation and may be interfaced to heading and/or speed sensors for dead reckoning during GPS outages. An advertised future option is a capability to interface to a Transit or Loran-C receiver.
A-14
EAGL]
INI-
MANUFACTURER: APPLICATION:
Motorola Inc. Air, land, and sea static, kinematic, and relative positioning
FEATURE;
4 channel, LI C/A-code/carrier tracking sensor 1 second update rate RS-232 for CDU (DCE) RS-232 for relative positioning data link (DTE) 18-32 VDC or 10-17 VDC (30 W max) Receiver: 1/4 ATR, 30 x 18 x 6 cm; 4.5 kg; -20C to +55C Antenna/preamp: 11 x 11 x 5 cm; 1.4 kg; -40C to +85C Options: 1 sec time mark 429 ARINC interface
AVAILABILITY: November 1986 $US 19 000
Canadian GPS Associates, April 1986
EAGLE MINI-RANGER
Manufacturer: Motorola Inc. 2100 E.Elliot Road Tempe, Arizona U.S.A.
Motorola is producing a 4 channel, LI, C/A-code sensor for ground, sea, and air, and static, kinematic and relative positioning applications [Motorola Inc., 1986a; 1986b]. The compact receiver (30 x 18 x 6 cm) comes with two RS-232 ports to facilitate differential operations. Relative positioning data, including satellite identification, satellite correction factors, measurement quality, and time, is updated for output every second through the data link port.
The 4.5 kg receiver comes with two DC power supply options and is operable at -20C to +55C. The 1.4 kg antenna/preamp unit operates at -40C to +70C.
A-15
MANUFACTURER:
APPLICATION:
FEATURES: 5 or (optional) 7 channel LI C/A-code operation 1 sec update rate > 1000 waypoints, > 1000 events Mass memory unit stores up to 8 hrs of 30 second data 4 lines of 40 character display (menu driven) 4 configurable serial I/O ports Phase, pseudo-range, positions, etc. output to any or all ports 13 x 35 x 30 cm; 24 VDC, 110/220 VAC optional -40C to +60C 8kg AVAILABILITY: Mid-1987, 7 channel end of 1987 $CAN 79 950 plus $13 995 for 7 channel set
Canadian GPS Associates, May 1987
GPS Receivers A. 15
NORSTAR 1000
Manufacturer:
Norstar Instruments Ltd. Second floor, 319 - 2nd Avenue S.W. Calgary, Alberta Canada T2P OC5
By the end of 1986, Norstar Instruments Ltd. expects to have completed development of the NORSTAR 1000, a ruggedized, compact, and portable (or optionally rack mountable) 5 (or optionally 7) channel, C/A-code receiver for point or relative (dynamic and static) positioning [Norstar Instruments, 1986]. One feature of this system is that it is menu driven and that it allows inspection and modification of many parameters including, for example, the standard deviation of the navigation filter parameters. As well, the three I/O ports can each be configured for data logging various raw and processed values to several peripherals, such as mass memory units, simultaneously. The receiver allows for the recording of code, carrier phase measurements, and broadcast ephemeris to either the mass memory unit (8 hours of data at 30 seconds) or to any I/O port for post-processing applications, such as precise geodetic relative positioning.
A-16
Prakla-Seismos AG
APPLICATION:
FEATURES: 1 channel 4 satellite multiplex C/A-code Waypoint navigation functions to local or remote displays Accepts input from log, gyro, modem Output capability to print, plot, or record devices Position averaging for static Manual or automatic relative mode operation Consists of Antenna unit: 20 x 40 cm, 4 kg; -20C to +70C Receiver unit (not for outside use): 0C to +55C CDU: 25 x 15 x 18 cm (8 kg)
End of 1986
Canadian GPS Associates, April 1986
GPS Receivers A. 16
PS 8400
Manufacturer:
Prakla-Seismos AG
Buchholzer Str. 100 D-3000 Hannover 51 Federal Republic of Germany
Prakla-Seismos, in a joint venture with Kongsberg Vapenfabrik of Norway, has already developed and tested two PS 8400F (functional) sequential C/A-code receivers for both static and dynamic applications to meet their respective requirements in geophysical exploration and offshore positioning [Stiller, 1985]. Kinematic positioning tests [Gerlach, 1985] yielded typical C/A-code accuracies of better than 30 m. By using relative positioning techniques for static positioning, 2
(49 x 25 x 27 cm) operates at 0C to +50C. Standard navigation and waypoint functions are available. Through an external modem and communication link, relative positioning can be done automatically in real time. In cooperation with the Geodetic Institute of the University of Bonn, Prakla-Seismos is also developing a GPS receiver for precise relative geodetic
positioning.
A-17
MANUFACTURER: APPLICATION:
Rockwell International
Timing, navigation, relative positioning
4 satellite, 1 channel, LI C/A-code sequential GPS sensor Complete satellite and positioning cycle every 1 second Three configurations Timing Navigation Positioning/differential Two RS-232 ports 1 for control/display 1 for processed and raw data out 20 x 19 x 12 cm, 3 kg, 0C to 55C 12-40 VDC (30 watts)
AVAILABILITY: $US 18 000
Canadian GPS Associates, April 1986
OPS Receivers A. 17
Rockwell International Collins Industrial GPS Products Cedar Rapids, Iowa 52498 U.S.A.
During the NAVSTAR GPS Phase IIB Full Scale Development program of the Join Program Office, Rockwell-Collins developed a family of user equipment, for a wide range of applications, consisting of one, two, and five channel sets [Krishnamurti et al., n.d.]. One-channel sets were designed for manpack and tank configurations; two-channel sets for helicopters and aircraft carriers; and five-channel configuration for submarines and aircraft, such as the F-16 and B-52. Drawing from these development experiences, Avionics Group of Rockwell International has produced the Collins NAVCORE I GPS navigation receiver (sensor) for the commercial market. The NAVCORE I, LI, C/A-code receiver sequences through 4 satellites in one second providing completely updated navigation solutions at up to once per second. Functionally, three GPS receiver units can be derived from the single GPS sensor hardware module: a timing receiver, a navigation receiver, and a position/differential receiver. The timing configuration outputs static position and time and time mark only (with respect to UTC) once per second. The navigation receiver configuraton outputs navigation solutions once per second, and GDOP information and timing every two minutes. The positioning/differential configuration is the most flexible, able to output everything the other sets do plus raw pseudo-range and delta range measurements and ephemeris data, some at variable rates. Physically, the sensor is available in one 20 x 19 x 12 cm or two half height modules with interconnect cables, weighing about 3 kg, requiring 12 to 40 VDC (30 Watts) power, and operating at 0C to +55C. The two RS-232 ports provide (a) control/display I/O, and (b) output of raw or processed data. A one pulse per second timing signal is also provided. Because the NAVCORE I is a GPS sensor only, necessary peripheral equipment must include a control display unit and antenna. Other commercial GPS products are under development.
A-18
MANUFACTURER:
Sercel-France
APPLICATION:
FEATURES:
5 channel, LI C/A-code Code, carrier phase measurements, and broadcast data available through 2 RS-232 ports 2- or 3-dimensional navigation Non-portable configuration Receiver: 37 x 53 x 45 cm; 30 kg; +9C to +40C, 24 DC Antenna/preamp: 5 kg, -20C to +50C Monitor (video): 28 x 25 x 43 cm; 11 kg; 0C to +40C, 220 VAC Keyboard: 6 x 37 x 26 cm; 1 kg; 0C to +40C
AVAILABILITY: 450,000 FF or $US 67 500
Canadian GPS Associates, April 1986
GPS Receivers A. 18
TR5S
Manufacturer: Contact: Sercel-France Sercel Electronics of Canada Ltd. Bay 1,1715 27th Avenue N.E. Calgary, Alberta Canada T2E7E1
Although the TR5S is not ruggedized or very portable, with its keyboard and video monitor configuration, it is suitable for situations where space and power are not major considerations. The receiver is most useful as a tool to evaluate the GPS capabilities in real time and postmission with the raw code and carrier phase measurements, broadcast data, and recorded processed positioning results. Data are available at a 0.6 second rate.
The TR5S is intended for low dynamics of 100 km/hr with accelerations of 0.5 g. For precise surveying applications, the antenna-preamp assembly includes a ground plane. The 30 kg receiver, of dimensions 37 x 53 x 45 cm, operates on 24 VDC at temperatures between +9C to +40C; the 5 kg antenna operates at temperatures between -20C to +50C; the 11 kg video monitor (28 x 25 x 43 cm) and the 1 kg keyboard (6 x 37 x 26 cm) both operate at 0C to +40C.
A-19
lANUFACTURER;
D
LICATI0N;
FEATURES:
4 channel, LI C/A-code GPIB interface (or RS-232 option) for position and ephemeris output Auto start and acquisition function 15 waypoints navigation functions LCD display Antenna: 7 x 28 cm; 1 kg; -20C to +75C Receiver unit: 16 x 21 x 30 cm; 7 kg; -10C to +50C
9-32 VDC (75 W) AVAILABILITY: Early 1987
GPS Receivers A. 19
GTT-2000
Manufacturer:
The Sony Corporation has already developed a time transfer receiver and a navigation receiver. Development of a receiver for precise geodetic applications has also begun. The GTT-2000 navigation receiver [Sony Corporation, n.d.] will be available commercially in early 1987. The 7 kg compact 16 x 21 x 30 cm receiver unit is available with three power supply options: 9-16 VDC, 18-32 VDC, and 9-32 VDC. It operates at -10C to +50C. The 1 kg, 7 cm diameter by 28 cm antenna operates at -20C to +75C. GPIB or RS-232 interface ports are available for computed position outputs.
A-20
;L
:ANUFACTURER;
'(TFIVFTR
fVssxlL^IL V JL>iJLvv
4 channel, LI C/A-code Receiver: 16 x 14 x 40 cm Antenna: 10 x 10 x 2 cm plus preamp Weight: 8 kg Power: 12/24 VDC (15 watts) For land navigation, auxilliary sensor package compensates for GPS signal loss
AVAILABILITY;
1989
< $5000
SEL is developing a 4 channel, LI, C/A-code GPS receiver for air, land, and marine applications. Only preliminary specifications are available. The compact, lightweight system is expected to sell for less than $5000 when it is made available in early 1989. For land navigation applications, an auxilliary sensor package will be incorporated to compensate for GPS signal loss due to buildings and trees.
A joint venture between SEL (Germany), TRT/LMT (France), Elmer (Italy), and Magnavox (U.S.A.) is producing an L1/L2 5 channel P-code set for use at high dynamics to be available in 1989.
A-21
[ANUFACTURER:
Down converter chassis for up to 6 dual receiver chassis, and test transmitter (STI 5012C) in double rack configuration Each dual receiver can track all L1/L2, C/A- and P-code combinations of measurements for 2 satellites continuously Supports anti-spoof acquisition and tracking Performance:
Measurement
pseudo-range
Mode
C/A
Resolution (m)
2.3 0.23
.0007
pseudo-range
Doppler
P
P
The first commercially available GPS receiver (1980) was the STI-5010, an L1/L2, C/A- and P-code slow sequencing receiver. To cycle through 4 satellites took about 5 minutes. The receiver unit alone required about 30 cm of rack space. The pseudo-range measurements were performed by external counters. The receiver was controlled and the position was computed by an external computer. The 5010 is excellent for experimenting with GPS, providing all raw measurements and data along with intermediate frequencies.
Because of the 5010 performance, STI was awarded the contract to produce the monitor station receivers (MSR) for the GPS control stations. Each MSR consists of two racks of equipment: 1 down converter chassis provides a signal for up to 6 dual receiver chassis; an STI-5012C test transmitter is also included to provide a calibration signal at regular intervals. Each receiver can track L1/L2, C/A-, P-, and future Y-codes to make pseudo-range and Doppler measurements.
Other STI GPS products include the STI-502, a time transfer receiver, the STI-5012C GPS test transmitter and, towards the end of 1986, GPS receiver modules (with configurations containing from 1 to 8 channels) for the OEM market.
A-22
MODEL 1
MANUFACTURER: Starnav Corporation
APPLICATIONS: Low to medium dynamics navigation Heading and attitude FEATURES: 5 channel, LI C/A-code Phase, pseudo-range and rate and satellite message output Kalman filter for position, velocity, and time Least-squares estimator for attitude and heading Waypoint navigation and route management 3-element antenna: 38 x 38 x 18 cm; 22.7 kg; -40C to +50C 2-piece receiver: 115 VAC (90 W) or 28 VDC (110 W); -20C to +50C signal processor: 27 x 23 x 32 cm; 13.5 kg display and keyboard: 32 x 22 x 28 cm; 9 kg Expected accuracies: position: 30 m worst case velocity: 0.1 m/s heading: 0.3 roll, pitch: 1 AVAILABILITY: Product announcement October 1986
SCAN 25 ooo
PAHRS MODEL 1
Manufacturer:
Starnav Corporation 148 Colonnade Road Suite 201 Ottawa, Ontario Canada K2E7R4
The Starnav Corporation, founded by a former Canadian Marconi employee, is developing a unique GPS receiver, the PAHRS Model 1. This receiver, in addition to providing the usual real-time position and velocity information, can also determine the attitude and heading of the vehicle, vessel, or aircraft to which the antenna is attached. This capability is achieved through the use of a novel three-element heated antenna which in principle works like an interferometer. The receiver switches between the three elements and code and phase tracks the signals from each element in the usual fashion. By suitably combining the outputs of the trackers, the microprocessor computes the heading, roll, and pitch of the antenna and hence of the vehicle. A Kalman filter is used to improve the raw, real-time estimates. The anticipated heading accuracy is 0.3, and roll and pitch accuracy isl. PAHRS is designed to replace a conventional mechanical gyrocompass. In addition to the antenna, PAHRS consists of a signal processing unit and a control and display unit. The size of the Model 1 signal processor is 27 x 23 x 32 cm and the weight is 13.5 kg. The control and display unit is 32 x 22 x 28 cm in size and weighs 9 kg. Starnav may reconfigure future models to reduce the size and weight of the units. PAHRS can run off 115 VAC or 28 VDC power and consumes about 100W. Software will be provided for waypoint navigation and route management.
A-23
MANUFACTURER:
Static to medium dynamics single point or relative positioning Precise geodetic relative positioning
FEATURES:
L1/L2, P-code and C/A-code multiplex receiver 8 pseudo-channels: 4 LI, 4 L2 Pseudo-range, carrier phase and message output Real-time navigation functions with 9 waypoints Ruggedized; 22-32 VDC (90 W); -20C to +50C Receiver unit: 37 x 45 x 21 cm; 24.5 kg Antenna unit: 30 x 17 cm dia.; 2 kg Optional recorder: 21 x 36 x 30 cm; 7 kg
AVAILABILITY:
The TI-4100, available since 1982, is the only commercially available C/A-code, P-Code, L1/L2 receiver. Because of these features and its single channel 4 satellite multiplex design, the TI 4100 is the most versatile receiver available today [Texas Instruments, 1986]. Successfully used and proven applications for the receiver include [Henson and Montgomery, 1984]: 3-D navigation with 4 satellites at stationary, low, and medium dynamics in ships, aircraft and on land with position accuracies better than 15m with L1/L2 P-code and velocity accuracies at 0.15 m/s to 2.00 m/s for medium dynamics. 3-D relative positioning (4 satellites) in aircraft, ships, and on land at distances up to hundreds of kilometres with accuracies of 2 to 5 m. With 3 satellites and altitude aiding 2-D accuracies of about 15m. With 3 satellites and external clock aiding 3-D positioning at about 15m. With 2 satellites and altitude and time aiding again about 15m 2-D positioning. 3-D geodetic precise relative positioning with accuracies from 1:100 000 to better than 1:1 000 000 depending on software used to postprocess the data and amount and quality of simultaneous data collected. Currently, when using C/A-ctde only, accuracies of 20 to 30 m are achieved when ionospheric refraction is minimal.
The ruggedized (-20C to +50C) portable or rack mountable TI 4100 system consists of a 24.5 kg, 37 x 45 x 21 cm receiver/processor with hand-held control and display unit (CDU); a 2 kg, 30 cm x 17 cm diameter antenna; an optional 7 kg dual cassette recorder (21 x 36 x 20 cm); and a 22 to 32 VDC power source (90 W) and cables. Both CDU and recorder (replacable by computer) are interfaced through RS-232 ports.
/
The receiver makes available to the recorder port the time tagged code, carrier beat phase, and carrier phase rate measurements along with satellite navigation messages, and positions and velocities. Various navigation functions are included, such as cross-track error, distance to go, etc., with storage for up to 9 waypoints.
TI also has available SATPLAN software which computes and graphically displays satellite rise and set times, geometry, etc., and GEOMARK for determining 3-D baselines from carrier phase data for geodetic applications through postmission processsing. Both programs run on the TI PC.
A-24
URVEYOR
MANUFACTURER: Trimble Navigation
APPLICATION: 4000A - low dynamics navigation, time transfer 4000S - plus precise relative positioning FEATURES: 4 channel, LI C/A-code 20 kg, 17 x 45 x 48 cm, rack mounted (4000S in protective casing) 0C to +50C 115/230 VAC or 24 VDC (4000A) or 12 VDC (4000S) Outputs ephemeris and positions, etc., pseudoranges. 4000S also outputs ^carrier phase and integrated Doppler Automatic search and acquisition 1.5 kg compact antenna, -30C to +75C 4000S includes ground plane Computer control and recording interface RS-232C (option on 4000A)
AVAILABILITY:
4000A - $US 25 000 1 month delivery 4000S - $US 44 000 2 months delivery Includes postprocessing software
Canadian GPS Associates, May 1987
Trimble Navigation 585 North Mary Avenue Sunnyvale, California 94086 U.S.A.
As of March 1986, Trimble had more GPS receiver models available to the civilian market at low cost than any other manufacturer. In 1984, Trimble introduced the 4000A, a 4 channel, LI, C/A-code navigator and time transfer receiver; a rack mountable 24 VDC, 115/230 VAC 50 watt receiver with a compact antenna which operates at -30C to +75C. Updated positions and velocities are available for output to the integral CDU or RS-232 I/O port at a 10 second rate. At the end of 1985, the 4000S was introduced. Enhanced, with more powerful processor and software features, the 4000S allows for recording of time tagged carrier phase, pseudo-range, and integrated Doppler measurements and broadcast ephemeris, and can be used for precise relative positioning in either postprocessing (currently available) mode or real time (planned option) [Trimble Navigation, 1985; 1986]. At least two receivers make up a survey system, each consisting of an encased 4000S, a recording system (usually a battery operable personal computer), an external 12 VDC power supply (60 watts), and the antenna and cable. Each receiver system weighs a total of less than 40 kg (without external power supply). For precise three-dimensional relative accuracies of about 1 cm + 2 ppm of the baseline distance of 1 to 100 km, at least 30 minutes of good continuous data from optimal satellite geometry is required. The postprocessing is usually done on the same computer used for recording the data, requiring about 15 minutes of processing for 30 minutes of simultaneous carrier phase data.
The accuracy of relative static positioning from pseudo-ranges is 2 m. 'The receiver is also capable of dynamic positioning to 200 knots.
Trimble has recently upgraded the capabilities of the 4000S and renamed it the 4000SX.
A-25
LORAN-GPS
MANUFACTURER:
APPLICATION:
FEATURES: LORAN-GPS combination 2 channel, LI C/A-code multiplexed GPS receiver Each CDU 25 x 12 x 7 cm Receiver/processor: 11 x 25 x 30 cm Weight with antenna: 7 kg 30 watts at 12 or 24 VDC -10C to +60C 100 waypoint navigation 1 second update rate
AVAILABILITY:
$US 16 000 to $US 25 000, depending on number of CDUs, software, with/without LORAN
Canadian GPS Associates, May 1987
LORAN-GPS 10X
Manufacturer:
Trimble Navigation 585 North Mary Avenue Sunnyvale, California 94086 U.S.A.
To increase continuous navigation capabilities of stand-alone GPS or LORAN, Trimble has produced a marine navigation system to combine the advantages of both GPS and LORAN in the LORAN-GPS 10X receiver [Braisted et al., n.d.b]. The two channel, LI, C/A-code multiplexed receiver portion of this system provides position and velocity updates at 1 second intervals, whereas the LORAN portion provides an update at 6 second intervals when GPS is not operational. The compact receiver/processor (11 x 25 x 30 cm) and two CDUs (25 x 12 x 7 cm) consume 30 watts at 12 or 24 VDC and weigh 7 kg with antenna.
Trimble Navigation also produces a GPS sensor for low dynamics applications [Braisted et al., n.d.a]. The Model 400, a 2 channel, LI, C/A-code sequential receiver is packaged in a 1/4 ATR short ARINC-sized package for interfacing to existing display and keyboard modules via dual RS-232 or ARINC 429 interfaces.
A-26
MANUFACTURER:
FUTURE OPTIONS: Meteorological sensor interface Modem interface L2 'codeless' carrier phase
AVAILABILITY: Includes postprocessing software $US 80 000/unit Spring 1986
Canadian GPS Associates, May 1987
WM101
Manufacturer: WM Satellite Survey Company c/o Wild Heerbrugg CH-9435 Heerbrugg Switzerland
Wild Heerbrugg and Magnavox have joined forces as the WM Satellite Survey Company and are offering as their first product the WM101 surveying receiver [Scherrer, 1985]. This is the first commercially available GPS receiver product in which the power supply, recorder, control display unit, and receiver/processor have been combined into a single compact robust unit weighing only 14 kg in a package slightly larger than a briefcase (50 x 40 x 17 cm). This, and the only additional piece required in the fieldthe 1.5 kg antenna of 20 cm heightmakes this system very portable.
Both antenna and receiver are waterproof and operable at temperatures of -25C to +55C and -40C to +70C, respectively. The antenna cable length options vary from 10 to 120 metres. The 4 channel, LI, C/A-code GPS receiver can be used for low dynamics navigation but is designed primarily for precise relative positioning by postprocessing the data. Up to 6 satellites are tracked from which pseudo-ranges, broadcast ephemeris, and carrier phase data can be recorded on cassette tape. To compute the three-dimensional vector between stations, WM supplies a software package, called PoPS, which runs on a personal computer.
Accuracies of 1 cm plus 2 ppm are expected using the broadcast ephemeris. WM plans to add codeless L2 carrier tracking along with a meteorological sensor interface and real-time baseline determination software as options sometime in the future.
TH:
A-27
other positioning receivers time transfer receivers receivers for military use
SERIES, SERIES-X Boojum Universities: New Brunswick Delft Bonn Stuttgart Leeds
For development and research purposes, receivers have been developed using off-the-shelf components at several universities, including the University of New Brunswick, and Delft, Leeds, Bonn, and Stuttgart universities.
Manufacturers of GPS receivers not represented in this Appendix or those introducing new GPS receivers are invited to send product information to:
Canadian GPS Associates Post Office Box 3184, Postal Station B Fredericton, New Brunswick Canada E3A 5G9.
The new information will be included in later editions of this publication.
APPENDIX B
Between-epoch difference The difference between two complete carrier beat phase measurements made by the same receiver on the same signal (same satellite, same frequency), but at different time epochs. Between-frequency difference The instantaneous difference between (or, more generally, any other linear combination involving) the complete carrier beat phase measurements made by the same receiver observing signals from the same satellite at two (or more) different frequencies. Between-receiver difference The instantaneous difference in the complete carrier beat phase measurement made at two receivers simultaneously observing the same received signal (same satellite, same frequency).
Between-satellite difference The instantaneous difference in the complete carrier beat phase measurement made by the same receiver observing two satellite signals simultaneously (same
frequency).
Binary pulse code modulation Pulse modulation using a string (code) of binary numbers. This coding is usually represented by ones and zeros with definite meanings assigned to them, such as changes in phase or direction of a wave [Dixon, 1976].
Binary biphase modulation Phase changes on a constant frequency carrier of either 0 or 180 (to represent binary 0 or 1, respectively). These can be modelled by y = A(t) cos (cot + <()), where the amplitude function A(t) is a sequence of +1 and -1 values (to represent 0 and 180 phase changes, respectively) [Dixon, 1976].
Carrier A radio wave having at least one characteristic (e.g., frequency, amplitude, phase) which may be varied from a known reference value by modulation [Bowditch, 1981, Vol. II].
Carrier frequency The frequency of the unmodulated fundamental output of a radio transmitter [Bowditch, 1981, Vol. II].
Carrier beat phase The phase of the signal which remains when the incoming Doppler-shifted satellite carrier signal is beat (the difference frequency signal is generated) with the nominally-constant reference frequency generated in the receiver.
Carrier beat phase ambiguity The uncertainty in the initial measurement, which biases all measurements in an unbroken sequence. The ambiguity consists of three components
i-
PJ is the fractional initial phase in the satellite (both due to various contributions to phase bias, such as unknown clock phase, circuit delays, etc.), and NjJ is an integer cycle bias in the initial measurement.
Channel A channel of a GPS receiver consists of the radio frequency and digital hardware, and the software, required to track the signal from one GPS satellite at one of the two GPS carrier frequencies.
Chip The minimum time interval of either a zero or a one in a binary pulse code.
C/A-code The standard (Coarse/Acquisition, or Clear/Access) GPS codea sequence of 1023 pseudo-random binary biphase modulations on the GPS carrier at a chip rate of 1.023 MHz, thus having a code repetition period of one millisecond.
Complete instantaneous phase measurement A measurement of carrier beat phase which includes the integer number of cycles of carrier beat phase since the initial phase measurement. See Fractional instantaneous phase measurement. Correlation-type channel A GPS receiver channel which uses a delay lock loop to maintain an alignment (correlation peak) between the replica of the GPS code generated in the receiver and the incoming code. Cycle slip A cycle slip is a discontinuity of an integer number of cycles in the measured carrier beat phase resulting from a temporary loss-of-lock in the carrier tracking loop of a GPS receiver. Delay lock The technique whereby the received code (generated by the satellite clock) is compared with the internal code (generated by the receiver clock) and the latter shifted in time until the two codes match. Delay lock loops can be implemented in several ways, for example, tau dither and early-minus-late gating [Spilker, 1978]. Delta pseudo-range The difference between two carrier beat phase measurements, made coincidentally with (code) pseudo-range epochs. Differenced measurements see Between-epoch difference; Between-frequency difference; Betweenreceiver difference; Between-satellite difference. Many combinations of differences are possible. Which differences, and their order, should be specified in describing a processing method (for example, Receiver-satellite double differences).
Differential positioning see Relative positioning Dilution of precision (DOP) A description of the purely geometrical contribution to the uncertainty in a kinematic position fix, given by the expression
DOP =
where A is the design matrix for the solution (dependent on satellite/receiver geometry). The DOP factor depends on the parameters of the position fix solution. Standard terms in the case of kinematic GPS are: GDOP (three position coordinates plus clock offset in the solution) PDOP (three coordinates) HDOP (two horizontal coordinates) VDOP (height only) TDOP (clock offset only), and HTDOP (horizontal position and tune).
Doppler shift The apparent change in frequency of a received signal due to the rate of change of the range between the transmitter and receiver. See Carrier beat phase. Dynamic positioning See Kinematic positioning Fast switching channel A switching channel with a sequence time short enough to recover (through software prediction) the integer part of the carrier beat phase. Fractional instantaneous phase measurement A measurement of the carrier beat phase which does not include any integer cycle count. It is a value between zero and one cycle. See Complete instantaneous phase measurement. Frequency band A range of frequencies in a particular region of the electromagnetic spectrum [Wells, 1974]. Frequency spectrum The distribution of amplitudes as a function of frequency of the constituent waves in a signal [Wells, 1974].
Handover word The word in the GPS message that contains time synchronization information for the transfer from the C/A-code to the P-code [Milliken and Zoller, 1980]. Independent baselines Baselines determined from independent observing sessions.
Independent observing sessions Sessions for which any common biases affecting the observations can be ignored.
Ionospheric refraction A signal travelling through the ionosphere (which is a nonhomogeneous and dispersive medium) experiences a propagation time different from that which would occur in a vacuum. Phase advance depends on electron content and affects carrier signals. Group delay depends on dispersion in the ionosphere as well, and affects signal modulation (codes). The phase advance and group delay are of the same magnitude but opposite sign [Davidson et al., 1983].
Interferometry See Relative positioning
Kinematic positioning Kinematic positioning refers to applications in which a trajectory (of a ship, ice field, tectonic plate, etc.) is determined.
Lane The area (or volume) enclosed by adjacent lines (or surfaces) of zero phase of either the carrier beat phase signal, or of the difference between two carrier beat phase signals. On the earth's surface, a line of zero phase is the locus of all points for which the observed value would have an exact integer value for the complete instantaneous phase measurement. In three dimensions, this locus becomes a surface. L-band The radio frequency band extending from 390 MHz to (nominally) 1550 MHz [Bowditchl981,Vol.II].
Multipath error An error resulting from interference between radiowaves which have travelled between the transmitter and the receiver by two paths of different electrical lengths [Bowditch, 1981, Vol. II]. Multichannel receiver A receiver containing many channels.
Multiplexing channel A receiver channel which is sequenced through a number of satellite signals (each from a specific satellite and at a specific frequency) at a rate which is synchronous with the satellite message bit-rate (50 bits per second, or 20 milliseconds per bit). Thus one complete sequence is completed in a multiple of 20 milliseconds. Observing session The period of time over which GPS data is collected simultaneously by two or more receivers.
Outage The occurrence in time and space of a GPS Dilution of Precision value exceeding a specified maximum.
Phase lock The technique whereby the phase of an oscillator signal is made to become a smoothed replica of the phase of a reference signal by first comparing the phases of the two signals and then using the resulting phase difference signal to adjust the reference oscillator frequency to eliminate phase difference when the two signals are next compared [Bowditch, 1981, Vol. II]. The smoothing time span occurs over approximately the inverse of the bandwidth. Thus a 40 hertz loop bandwidth implies an approximately 25 millisecond smoothing time constant. Phase observable See Carrier beat phase
P-code The Precise (or Protected) GPS codea very long (about 1014 bit) sequence of pseudo-random binary biphase modulations on the GPS carrier at a chip rate of 10.23 MHz which does not repeat itself for about 267 days. Each one-week segment of the P-code is unique to one GPS satellite and is reset each week.
Precise positioning service (PPS) The highest level of dynamic positioning accuracy that will be provided by GPS, based on the dual frequency P-code [U.S. DoD/DOT, 1982]. Pseudolite The ground-based differential GPS station which transmits a signal with a structure similar to that of an actual GPS satellite [Kalafus, 1984]. Pseudo-random noise (PRN) code Any of a group of binary sequences that exhibit noise-like properties, the most important of which is that the sequence has a maximum autocorrelation, at zero lag [Dixon, 1976].
Pseudorange The time shift required to align (correlate) a replica of the GPS code generated in the receiver with the received GPS code, scaled into distance by the speed of light. This time shift is the difference between the time of signal reception (measured in the receiver time frame) and the time of emission (measured in the satellite time frame).
Pseudo-range difference See Carrier beat phase Receiver channel See Channel
Restart capability The property of a sequential processing computer program, that data can be processed rigorously in a sequence of computer runs, rather than only in one long run.
S-code See C/A-code
Satellite constellation The arrangement in space of the complete set of satellites of a system like GPS. Satellite configuration The state of the satellite constellation at a specific time, relative to a specific user or set of users. Simultaneous measurements Measurements referred to time frame epochs which are either exactly equal or else so closely spaced in time that the time misalignment can be accommodated by correction terms in the observation equation, rather than by parameter estimation.
Slow switching channel A switching channel with a sequencing period which is too long to allow recovery of the integer part of the carrier beat phase.
Spread spectrum systems A system in which the transmitted signal is spread over a frequency band much wider than the minimum bandwidth needed to transmit the information being sent [Dixon, 1976].
Squaring-type channel A GPS receiver channel which multiplies the received signal by itself to obtain a second harmonic of the carrier, which does not contain the code modulation.
Standard positioning service (SPS) The level of kinematic positioning accuracy that will be provided by GPS based on the single frequency C/A-code [U.S. DoD/DOT, 1982].
Static positioning Positioning applications in which the positions of points are determined, without regard for any trajectory they may or may not have. Switching channel A receiver channel which is sequenced through a number of satellite signals (each from a specific satellite and at a specific frequency).
source, converted into range units, assuming that error source is uncorrelated with all other error sources [Martin, 1980].
Z-count word The GPS satellite clock time at the leading edge of the next data subframe of the transmitted GPS message (usually expressed as an integer number of 1.5 second periods) [van Dierendock et al., 1978].
APPENDIX C
GPS INFORMATION SOURCES This appendix contains sources of information on GPS of interest to those who wish to learn more about GPS or to use it in practice. Sources in three categories are listed: publications complementary to this Guide, information sources on the current status of GPS, and additional GPS services available from Canadian GPS Associates.
OTHER PUBLICATIONS
Everyman's Guide to Satellite Navigation A booklet designed to acquaint managers with principles involved in using GPS for navigation. August 1985. 66 pages. Single copy free. Available from Steven D. Thompson ARINC Research Corp. 2551 Riva Road Annapolis, MD 21401, U.S.A. Telephone: (301) 266-4650 The WM GPS Primer Information required for basic understanding of GPS satellite surveying, highlighting the WM101 receiver and its postprocessing software, PoPS. 1985. 116 pages. Single copy free. Available from Rene Scherrer (3388) WM Satellite Survey Company Wild Heerbrugg Geodesy Division CH-9435 Heerbrugg, Switzerland Surveying with GPS by R.W. King, E.G. Masters, C. Rizos, A. Stolz, and J. Collins. Reference book on GPS designed for the surveying community. November 1985. 128 pages. $AUS20. ISBN 0-85839-042-6. Available as Monograph No. 9 School of Surveying The University of New South Wales P.O. Box 1 Kensington, NSW, Australia 2033
Positioning with GPS 1985 Proceedings of the first international symposium on precise positioning with GPS, Rockville, Md, April 1985. Volumes I and II, 946 pages. $US 24 the set. Available as Publication N/CG 17X2 National Geodetic Information Center National Ocean Service, NOAA Rockville, MD 20852, U.S.A. Telephone: (301) 443-8316
Proceedings of the Fourth International Geodetic Symposium on Satellite Positioning Austin, Tx, May 1986.
Volumes I and II, 1510 pages. $US 50 the set. Available from Attn: Ms Barbara Becker Applied Research Laboratories University of Texas at Austin PO Box 8029 Austin, TX 78713-8029, U.S.A. Telephone: (512) 835-3466
Federal Radionavigation Plan Official U.S. plan describing the impact of GPS on other publicly supported radionavigation systems, and specifying the civilian availability of GPS. Issued approximately every two years. 1984 edition available as DoT-TSC-RFPA-84-8 Accession number AD-A151 295 National Technical Information Service Springfield, VA 22161, U.S.A.
STATUS OF GPS
GPS satellite status and health NAVSTAR/GPS Operational Control System Attn: Rufus N. Biggs, 2SCS/DOO Falcon Air Force Station Colorado Springs, CO 80914-5000, U.S.A.
Recorded daily operations report (303) 550-2099 Mission operations controller (live) (303) 550-2363 or 2477 Telexed or mailed satellite status/health reports (Notice Advisory to NAVSTAR Users (NANU)) sent to those on AIG 5825 distribution list
GPS satellite clock behaviour and related GPS information U.S. Naval Observatory Washington, DC 20392-5100, U.S.A.
Dial-up on-line computer files Requires terminal/computer/modem set to full duplex, 7 data bits, even parity, upper case. 300 or 1200 baud: (202) 653-1079; international CCITT.21: (202) 653-1095 2400 baud or CCITT.21 or CCITT.22 bis: (202) 653-1783. Files of interest: @GPS: GPS status information @GPSSY: GPS system description and more permanent information
USNO Series 4 bulletins contain information on the status of GPS including GPS time steering.
Approximate GPS orbit information (can be in error by tens of km) NASA Prediction Bulletin subscriptions Control Center Support Section Code 513.2 Project Operations Branch NASA Goddard Space Flight Center Greenbelt, MD 20771, U.S.A.
Request should specify NASA satellite identification numbers (see our Figure 4.02, column 5) of satellites for which prediction bulletins are desired.
Accurate GPS orbit information NGS: Precise orbital positions and velocities for GPS satellites are available to the U.S. public from the U.S. National Geodetic Survey. Satellite ephemerides are scheduled to be available within four weeks after the orbital data are collected. Contact Chief National Geodetic Information Center, N/CG17 National Geodetic Survey National Ocean Service, NOAA Rockville, MD 20852, U.S.A. Telephone: (301)443-8775
Litton Aero Service: Aero Service, a division of Western Geophysical and a member of the Litton Resources Group, offers a commercial orbit data service utilizing data obtained at its tracking network stations in the United States. Contact Mr. Jim Cain Manager, GPS Services Aero Service 8100 Westpark, Postal Box 1939 Houston, TX 77251, U.S.A. Telephone: (713) 784-5800
Civil access to the GPS Precise Positioning Service Assistant Secretary of Defense, Command, Control, Communications and Intelligence Attention: Colonel Phillip J. Baker The Pentagon Room3D174 Washington, DC 30201-3040, U.S.A.
Performance of GPS geodetic receivers Larry Hothem, Test Coordinator Federal Geodetic Control Committee 6001 Executive Boulevard Room 305/C1 Rockville, MD 20852, U.S.A. Telephone: (301)443-8171
Interface Control Document ICD-GPS-200 Official description of the GPS signal structure and specifications. September 1984. Available from Mail Code SN-65 Rockwell International Corporation Satellite Systems Division 12214 Lakewood Boulevard
Downey,CA 90241, U.S.A.
Courses tailored to specialized audiences Presented anywhere in the world 8.5" x 11" transparencies of figures from this volume SCAN 500 per complete set of 262. Contact
Canadian GPS Associates Box 3184 Postal StationB Fredericton, New Brunswick Canada E3A 5G9
APPENDIX D
BIBLIOGRAPHY
Bibliography D.I
APPENDIX D BIBLIOGRAPHY
anon. (1986). "FCC resolves issue of satellite position fixing technology." Aviation Week and Space Technology, Vol. 124, No. 24, p. 115.
Ayers, H. and M.D. Gonthier (1986). "JMR SatTrak GPS receiver." GPS seminar, 13-14 March, Qu6bec City, Quebec, Canada
Bauersima, I. (1983). "NAVSTAR/Global Positioning System (GPS) III: Erdvermessung durch radiointerferometriche Satellitenbeobachtungen." Miteilungen der Satellitenbeobachtungstation Zimmenvald Nr. 12, Druckerei der Universitat Bern, Switzerland.
Beck, N. (1985). "Practical aspects of satellite positioning." Presented at seminar on Impact of GPS on Surveying Practice in Ontario, Erindale College, Mississauga, Ontario, Canada, November.
Beutler, G. (1984). "GPS carrier phase difference observation software." Department of Surveying Engineering Technical Memorandum TM-2, University of New Brunswick, Fredericton, N.B., Canada.
Beutler, G., D.A. Davidson, R.B. Langley, R. Santerre, P. Vanicek and D.E. Wells (1984). "Some theoretical and practical aspects of geodetic positioning using carrier phase difference observations of GPS satellites." Department of Surveying Engineering Technical Report No. 109, University of New Brunswick, Fredericton, N.B., Canada
Bomford, G. (1980). Geodesy. 4th ed. Oxford University Press, Oxford.
Bbssler, J.D., C.C. Goad and P.L. Bender (1980). "Using the Global Positioning System (GPS) for geodetic positioning." Bulletin Geodesique, 54, pp. 553-563.
Bowditch, N. (1981). "Useful tables, calculations, glossary of marine navigation." Vol. II, The American Practical Navigator, Publication No. 9, DMAHTC, Washington, DC, U.S.A.
Braisted, P., R. Eschenbach and R. Sutherland (n.d.a). "$1,000 GPS receiver technology A fifth generation GPS receiver." Trimble Navigation, Sunnyvale, Ca, U.S.A.
Bibliography D.2
Braisted, R., R. Eschenbach and A. Tivori (n.d.b). "Combining LORAN and GPSThe best of both worlds." Trimble Navigation, Sunnyvale, Ca, U.S.A. Broughton, D.W., R.I. Essai and PJ. Farrow (1984). "Principles of satellite orbits." In: Global Civil Satellite Navigation Systems, Collected papers of NAV'84, Royal Institute of Navigation, London, U.K., May, pp. 3-6.
Brouwer, D. and G.M. Clemence (1961). Methods of Celestial Mechanics. Academic Press, NY, U.S.A.
Buennagel, L.A., P.P. MacDoran, R.E. Neilan, DJ. Spitzmesser and L.E. Young (1984). "Satellite Emission Range Inferred Earth Survey (SERIES) project: Final report on research and development phase, 1979 to 1983." Jet Propulsion Laboratory Publication 84-16, Pasadena, Ca, U.S.A., March. Buffett, B.A. (1985). "A short-arc orbit determination for the Global Positioning System." M.Sc. thesis, Division of Surveying Engineering, The University of Calgary, Calgary, Alberta, Canada.
Canadian Marconi Company (1985). "CMC capabilities and experience." CMC document No. 1162-1042, December, Kanata, Ontario, Canada.
Cannon, M.E., K.P. Schwarz and R.V.C. Wong (1986). "Kinematic positioning with GPSAn analysis of road tests." Proceedings of the Fourth International Geodetic Symposium on Satellite Positioning, University of Texas at Austin, Austin, Tx, U.S.A., April, pp. 1251-1268.
Clark, T.A., B.E. Corey, J.L. Davis and others (1985). "Precision geodesy using the Mark-Ill very-long-baseline interferometry system." IEEE Transactions, Geoscience and Remote Sensing, Special Issue on satellite geodynamics, Vol. GE-23, No. 4, pp. 438-449. Clynch, J.R. and D.S. Coco (1986). "Error characteristics of high quality geodetic GPS measurements: Clocks, orbits and propagation effects." Proceedings of the Fourth International Geodetic Symposium on Satellite Positioning, University of Texas at Austin, Austin, Tx, U.S.A., April, pp. 539-556.
Collins, J. (1986). "Measuring platform subsidence using GPS satellite surveying." Proceedings of Colloquium IV, Land, Sea and Space Today's Survey Challenge, CPA, CHA, Lake Louise, Alberta, Canada, April, pp. 377-382. Available from the Canadian Institute of Surveying and Mapping, Ottawa, Canada.
Bibliography D.3
Counselman, C.C. and S. Gourevitch (1981). "Miniature interferometer terminals for earth surveying: Ambiguity and multipath with Global Positioning System." IEEE Transactions, Geoscience and Remote Sensing, Vol. GE-19, No. 4, pp. 244-252.
Davidson, D., D. Delikaraoglou, R. Langley, B. Nickerson, P. Vanicek and D. Wells (1983). "Global Positioning System differential positioning simulations." Department of Surveying Engineering Technical Report No. 90, University of New Brunswick, Fredericton, N.B., Canada, March.
Degnan, J.J. (1985). "Satellite laser ranging: Current status and future prospects." IEEE Transaction, Geoscience and Remote Sensing, Special Issue on satellite geodynamics, Vol. GE-23, No. 4, pp. 398-413.
Delikaraoglou, D. and R.R. Steeves (1985). "The impact of VLBI and GPS on geodesy in Canada." Proceedings of the First International Symposium on Precise Positioning with the Global Positioning System, Rockville, Md, U.S.A., April, Vol. II, pp. 743-752.
Delikaraoglou, D., D.J. McArthur, P. H6roux and N. Beck (1986). "Improvements in the second generation GPS software of the Canadian Geodetic Survey." Proceedings of the Fourth International Geodetic Symposium on Satellite Positioning, University of Texas at Austin, Austin, Tx, U.S.A., April, pp.821-839.
Dixon, R.C. (1976). Spread Spectrum Systems. John Wiley, New York, U.S.A.
European Space Agency (1985). "A global civil navigation satellite system concept." ESA submission to ICAO and IMO working groups on navigation systems. Fruehauf, H.J. (1982). "Precision oscillators and their role and performance in navigation systems." Proceedings of IEEE PLANS '82, pp. 289-301.
Gerlach, B.E. (1985). "Positioning with a GPS pseudorange receiverFunctional model test results." Proceedings, Inertial, Doppler and GPS Measurements for National and Engineering Surveys, Munich, F.R.G., July.
Goad, C.C. (1985). "Precise relative position determination using Global Positioning System carrier phase measurements in a nondifference mode." Proceedings of the First International Symposium on Precise Positioning with the Global Positioning System, Rockville, Md, U.S.A., April, Vol. I, pp. 347-356.
Bibliography D.4
Goad, C.C. and B.W. Remondi (1984). "Initial relative positioning results using the Global Positioning System." Bulletin Geodesique, Vol. 58, No. 2, pp. 193-210. Gold, R. (1967). "Optimal binary sequences for spread spectrum multiplexing." IEEE Transactions on Information Theory, pp. 619-621. Goldfarb, J.M. and K.P. Schwarz (1985). "Kinematic positioning with an integrated INS-differential/GPS." Proceedings of the First International Symposium on Positioning with the Global Positioning System, Rockville, Md, U.S.A., April, Vol. II, pp. 757-772. Golumb, S.W. (1967). Shift Register Sequences. Holden-Day, San Francisco, Ca. Hatch, R. and K. Larson (1985). "MAGNET-4100 GPS survey program processing techniques and test results." Proceedings of the First International Symposium on Precise Positioning with the Global Positioning System, Rockville, Md, U.S.A., April, Vol. I, pp. 285-298. Hellwig, H. (1979). "Microwave time and frequency standards." Radio Science, Vol. 14, No. 4, pp. 561-572. Henriksen, S.W. (ed.) (1977). National Geodetic Satellite Program. Vols. 1 and 2, NASA SP-365, American Geophysical Union, Washington, DC, U.S.A. Henson, D. and B. Montgomery (1984). "TI 4100 Navstar Navigator test results." Presented at the 40th Annual Meeting of The (U.S.) Institute of Navigation, Cambridge, Ma, U.S.A., June. Hothem, L.D. (1985). "Proposed revisions to geodetic survey standards and preliminary specifications for geodetic surveys using relative positioning GPS techniques." National Geodetic Survey, version of 27 October, Rockville, Md, U.S.A. Hothem, L.D., C.C. Goad and B.W. Remondi (1984). "GPS satellite surveyingPractical aspects." The Canadian Surveyor, Vol. 38, No. 3, pp. 177-192. Japan Radio Corporation (1985). "JLR-4000 GPS navigator instruction manual." Tokyo, Japan. Johnson, C.R., P.W. Ward, M.D. Turner and S.D. Roemerman (1981). "Applications of a multiplexed GPS user set." In: Global Positioning System. Papers published in Navigation, reprinted by The (U.S.) Institute of Navigation, Vol.11, 1984, pp. 61-77.
Bibliography D.5
Jordan, E.D. (ed.) (1985). Reference Data for Engineers: Radio, Electronics, Computer, and Communications. 7th ed., Howard W. Sams and Co., Inc., Indianapolis, U.S.A.
Jorgensen, P.S. (1980). "NAVSTAR/Global Positioning System 18-satellite constellation." In: Global Positioning System. Papers published in Navigation, reprinted by The (U.S.) Institute of Navigation, Vol.II, 1984, pp. 1-12.
Kalafus, R.M. (1984). "RTCM SC-104 progress on differential GPS standards." Presented at the 40th Annual Meeting of The (U.S.) Institute of Navigation, Cambridge, Ma, U.S.A., June.
Kalafus, R.M., J. Vilcans and N. Knable (1983). "Differential operation of NAVSTARGPS." Navigation, Vol. 30, pp. 187-204.
Kaula, W.M. (1966). Theory of Satellite Geodesy. Blaisdell, Waltham, Ma, U.S.A.
King, R.W., R.I. Abbot, Y. Bock and C.C. Counselman (1985). "Deformation measurement by radio interferometry with GPS." EOS, Transactions of the AGU, Vol. 66, No. 46, p. 854 (abstract only).
Kleusberg, A. (1986). "Kinematic relative positioning using GPS code and carrier beat phase observations." Marine Geodesy, 10, pp. 167-184.
Kleusberg, A., S.H. Quek, D.E. Wells and J. Hagglund (1985a). "Comparison of INS and GPS ship velocity determination." Proceedings of the Third International Symposium on Inertial Technology for Surveying and Geodesy, Banff, Alberta, September. Division of Surveying Engineering, UCSE Report No. 60005, University of Calgary, Calgary, Alberta, Canada, 1986. Kleusberg, A., S.H. Quek, D.E. Wells, G. Lachapelle and J. Hagglund (1986). "GPS relative positioning techniques for moving platforms." Proceedings of the Fourth International Geodetic Symposium on Satellite Positioning, University of Texas at Austin, Austin, Tx, U.S.A., April, pp. 1299-1310. Kleusberg, A., R.B. Langley, R. Santerre, P. Vani2ek, D.E. Wells and G. Beutler (1985b). "Comparison of survey results from different types of GPS receivers." Proceedings of the First International Symposium on Positioning with GPS, Rockville, Md, U.S.A., April, pp. 579-592.
Krishnamurti, G., S.A. Harshbarger and T.N. Smith (n.d.). "The design and performance of GPS Phase II user equipment navigation software." Rockwell International, Cedar Rapids, la, U.S.A.
Bibliography D.6
Lachapelle, G. (1977). "Estimation of disturbing potential components using a combined integral formulae and collocation approach." In: Collected Papers of the Geodetic Survey of Canada, Surveys and Mapping Branch, Ottawa, Ontario, Canada.
Lachapelle, G. (1979). "Comparison of Doppler-derived and gravimetric geoid undulations in North America." Proceedings of the Second International Geodetic Symposium on Satellite Positioning, Austin, Tx, U.S.A., April, pp.653-670. Lachapelle, G. (1986). "GPSCurrent capabilities and prospects for multipurpose offshore surveying." Proceedings of FIG XVIII International Congress, Toronto, Ontario, Canada, June, Commission 8, pp. 86-101.
Lachapelle, G. and E. Cannon (1986). "Single and dual frequency GPS results for baselines of 10 to 500 km." The Canadian Surveyor, Vol. 40, No. 2, pp. 173-183.
Lachapelle, G. and R.L. Wade (1982). "NAVSTAR/GPS single point positioning." Technical papers of the 42nd Annual Meeting of the American Congress on Surveying and Mapping, Denver, Co, U.S.A., March, pp. 603-609.
Lachapelle, G., J. Lethaby and M. Casey (1984a). "Airborne single point and differential GPS navigation for hydrographic bathymetry." The Hydrographic Journal (U.K.), 37(2):200-208.
Lachapelle, G., J. Hagglund, H. Jones and M Eaton (1984b). "Differential GPS marine navigation." Proceedings of the IEEE PLANS '84, San Diego, Ca, U.S.A., November, pp. 245-255. Lachapelle, G., J. Hagglund, W. Falkenberg, P. Bellemare, M. Casey and M. Eaton (1986). "GPS land kinematic positioning experiences." Proceedings of the Fourth International Geodetic Symposium on Satellite Positioning, University of Texas at Austin, Austin, Tx, U.S.A., April, pp. 1327-1344. Langley, R.B., G. Beutler, D. Delikaraoglou, B. Nickerson, R. Santerre, P. Vanic'ek and D.E. Wells (1984). "Studies in the application of the Global Positioning System to differential positioning." Department of Surveying Engineering Technical Report No. 108, University of New Brunswick, Fredericton, N.B., Canada, September.
Leeman, R.W. and L.B. Fletcher (1985). "GPSA user's critique." Proceedings of the First International Symposium on Precise Positioning with the Global Positioning System, Rockville, Md, U.S.A., April, Vol. II, pp. 613-624.
Bibliography D.7
Lerch, F.J., B. Putney, S. Klosko and C. Wagner (1981). "Goddard earth models for oceanographic applications (GEM 10B and 10C)." Marine Geodesy, No. 5, pp. 145-198. Levitus, S. (1982). "Climatological atlas of the world ocean." NOAA Professional Paper No. 13, U.S. Department of Commerce, Washington, DC, U.S.A. Lindlohr, W. and D.E. Wells (1985). "GPS design using undifferenced carrier beat phase observations." Manuscripta Geodaetica, Vol. 10, No. 3, pp. 255-295.
Lisitzin, E. (1974). Sea Level Changes. Elsevier Scientific Publications Co., Amsterdam, The Netherlands.
MacDoran, P.P., J.H. Whitcomb and R.B. Miller (1984). "Codeless GPS positioning offers sub-meter accuracy." Sea Technology, October. Also in CSTG Bulletin No. 7, December, pp. 192-194.
MacDoran, P.P., R.B. Miller, L.A. Buennagel and J.H. Whitcomb (1985a). "Codeless systems for positioning with NAVSTAR-GPS." Proceedings of the First International Symposium on Precise Positioning with the Global Positioning System, Rockville, MD, April, Vol. I, pp. 181-190.
MacDoran, P.P., J.H. Whitcomb, R.B. Miller and L.A. Buennagel (1985b). "Codeless GPS systems for marine positioning and navigation." 17th Annual OTC, Houston, Tx, U.S.A., May (OTC 4994). Mader, G.L. (1986). "Decimeter level aircraft positioning using GPS carrier phase measurements." Proceedings of the Fourth International Symposium on Satellite Positioning, University of Texas at Austin, Austin, Tx, U.S.A., April, pp.1311-1326. Marsh, J.G., A.C. Bremner, B.D. Beckley and T.V. Martin (1984). "Global mean sea surface based on the SEASAT altimeter data." Preprint.
Martin, E.H. (1980). "GPS user equipment error models." In: Global Positioning System. Papers published in Navigation, reprinted by The (U.S.) Institute of Navigation, pp. 109-118.
McArthur, D., N. Beck, K. Lockhead and D. Delikaraoglou (1985). "Precise relative positioning with the Macrometer V-1000 interferometric surveyor: Experiences at the Geodetic Survey of Canada." Proceedings of the First International Symposium on Precise Positioning with the Global Positioning System, Rockville, Md, U.S.A., April, Vol. I, pp. 521-532.
Bibliography D.8
McDonald, K.D. and R.L. Greenspan (1985). "A survey of GPS satellite system alternatives and their potential for precise positioning." Proceedings of the First International Symposium on Precise Positioning with the Global Positioning System, Rockville, Md, U.S.A., April, Vol. I, pp. 457-472.
McPhersbn, W.R. (1981). "SARSAT applications." Canadian Aeronautics and Space Journal Vol. 27, No. 2, pp. 118-126.
Merrifield, J.T. (1985). "Fifth-generation GPS receiver aims at commercial market." Aviation Week and Space Technology, 19 August, pp. 123,125.
Mertikas, S. (1985). "Differential Global Positioning System navigation: A geometrical analysis." Department of Surveying Engineering Technical Report 95, University of New Brunswick, Fredericton, N.B., Canada.
Mertikas, S.P., D. Delikaraoglou and R. Santerre (1986). "Alert program for NAVSTAR Global Positioning System, Transit, LAGEOS, and Starlette satellites." Department of Surveying Engineering Technical Report No. 85, University of New Brunswick, Fredericton, N.B., Canada
Milliken, RJ. and CJ. Zoller (1980). "Principle of operation of NAVSTAR and system characteristics." In: Global Positioning System. Papers published in Navigation, reprinted by The (U.S.) Institute of Navigation, Vol.I, 1980, pp. 3-14.
Moreau, R. (1986). "The Quebec 1985 TI 4100 test." Proceedings of the Fourth International Geodetic Symposium on Satellite Positioning, University of Texas at Austin, Austin, Tx, U.S.A., April, pp. 1059-1072.
Moreau, R., J.G. Leclerc, B. Labrecque, G. Beutler, R.B. Langley and R. Santerre (1985). "The Quebec 1984 Macrometer test." Proceedings of the First International Symposium on Precise Positioning with the Global Positioning System, Rockville, Md, U.S.A., April, Vol. II, pp. 557-566. Morgan, V.G. (1983). "The challenge of precisely positioning a 3D seismic survey." Navigation, 30(3):261-272. Motorola Inc. (1986a). "Preliminary specification GPS Core receiver (Rev. C)." March.
Motorola Inc. (1986b). "Eagle Mini-Ranger Global Positioning System receiver." Equipment brochure.
Bibliography D.9
Nakiboglu, S.M., B. Buffett, K.P. Schwarz, EJ. Krakiwsky and B. Wanless (1984). "A multi-station approach to GPS orbital improvement and precise positioning." Proceedings, Proceedings of IEEE PLANS '84, San Diego, Ca, November, pp. 153-162.
NATO (1983). "Method of expressing navigation accuracies." North Atlantic Treaty Organization Standardization Agreement 4278.
Norstar Instruments (1986). "Norstar 1000 operational executive summary (preliminary)." Calgary, Alberta, Canada, March.
O'Neill, O.K. (1984). "The Geostar system." In: Global Civil Satellite Navigation Systems, Collected papers of NAV'84 conference, Royal Institute of Navigation, London, U.K., May, pp. 339-345.
OToole, J.W. (1976). "The CELEST computer program for computing satellite orbits." U.S. Naval Surface Weapons Center TR 3565, Dahlgren, Va, U.S.A.
Radio Technical Commission for Marine Services (1985). "Recommendations of Special Committee 104 [Chaired by R.M. Kalafus] differential NAVSTAR/GPS service." RTCM, P.O. Box 19087, Washington, DC 20036, U.S.A. Rapp, R.H. (1982). "A global atlas of sea surface heights based on the adjusted SEASAT altimeter data." Department of Geodetic Science and Surveying Technical Report No. 333, The Ohio State University, Columbus, Oh, U.S.A.
Rapp, R.H. (1983). "The development of the January 1983 1 x 1 mean free-air anomaly data tape." Internal Report, Department of Geodetic Science and Surveying, The Ohio State University, Columbus, Oh, U.S.A.
Remondi, B.W. (1985a). "Performing centimeter accuracy relative surveys in seconds using GPS carrier phase." Proceedings of the First International Symposium on Precise Positioning with the Global Positioning System, Rockville, Md, U.S.A., April, Vol. II, pp. 789-798. Remondi, B.W. (1985b). "Performing centimeter-level surveys in seconds with GPS carrier phase: Initial results." NOAA Technical Memorandum, NOS-NGS 45, Rockville, Md, U.S.A. Also ^Navigation, Vol. 32, No. 4, pp. 386-400.
Remondi, B.W. (1985c). "Modeling the GPS carrier phase for geodetic applications." Proceedings of the First International Symposium on Precise Positioning with the Global Positioning System, Rockville, Md, U.S.A., April, Vol, I, pp. 325-336.
Bibliography D.10
Remondi, B.W. (1985d). "Using the Global Positioning System (GPS) phase observable for relative geodesy: Modeling, processing, and results." Ph.D. dissertation, Center for Space Research, The University of Texas at Austin, Austin, Tx, U.S.A., May. Available from National Geodetic Information Center, N/CG17X2, Rockville, Md 20852, U.S.A. Rockwell (1984). "NAVSTAR GPS space segment/navigation user interfaces." Interface Control Document ICD-GPS-200, Rockwell International Space Systems Group, Satellite Systems Divison, dated 9/25/1984. Russell, S.S. and J.H. Schaibly (1978). "Control segment and user performance." In: Global Positioning System. Papers published in Navigation, reprinted by The (U.S.) Institute of Navigation, Vol. I, 1980, pp. 74-80.
Scherrer, R. (1985). "The WM GPS primer." WM Satellite Survey Company, Heerbrugg, Switzerland.
Schwarz, K.P. (1983). "Kalman filtering and optimal smoothing." In: Papers for the CIS Adjustment and Analysis Seminars, ed. EJ. Krakiwsky, Canadian Institute of Surveying, Ottawa, Ontario, Canada. Schwarz, K.P. (1986). "Geoid profiles from an integration of GPS-satellites and inertia! data." Paper presented at the International Symposium on the Definition of the Geoid, Florence, Italy, May. Schwarz, K.P., C.S. Eraser and P.C. Gustafson (1984a). "Aerotriangulation without ground control." International Archives of Photogrammetry and Remote Sensing, 25, Part Al, Rio de Janeiro, Brazil, June.
Schwarz, K.P., M.G. Sideris and R. Forsberg (1985). "Orthometric heights without levelling." Presented at the Annual Meeting of the Canadian Institute of Surveying, Edmonton, Alberta, Canada.
Schwarz, K.P., R.V.C. Wong, J. Hagglund and G. Lachapelle (1984b) "Marine positioning with a GPS-aided inertial navigation system." Proceedings of the National Technical Meeting of the Institute of Navigation, San Diego, Ca, U.S.A., January.
Scott, V.D. and J.G. Peters (1983). "A standardized exchange format for NAVSTAR GPS geodetic data." Applied Research Laboratories, The University of Texas at Austin, Austin, Tx, U.S.A., March.
Bibliography D.ll
Scott, V.D., J.R. Clynch and J.G. Peters (1986). "A proposed standardized exchange format for NAVSTAR GPS geodetic data." Preliminary paper, Applied Research Laboratories, University of Texas at Austin, Austin, Tx, U.S.A.
Service Argos (1978). "Argos user's guide." Centre National d'Etudes Spatiale, Toulouse, France.
Smith, D.E., D.C. Christodoulidis, R. Kolenkiewicz, P.J. Dunn, S.M. Klosko, M.H. Torrence, S. Fricke and S. Blackwell (1985). "A global geodetic reference frame from LAGEOS ranging (SL5.1AP)." Journal of Geophysical Research, Vol. 90, No. Bll, pp. 9221-9233. Sony Corporation (n.d.). "Instruction manual for GTT-2000 GPS receiver." (Preliminary.)
Spilker, JJ. (1978). "GPS signal structure and performance characteristics." In: Global Positioning System. Papers published in Navigation, reprinted by The (U.S.) Institute of Navigation, Vol. I, 1980, pp. 29-54.
Stansell, T.A. (1978). The Transit Navigation Satellite System. Magnavox Government and Industrial Electronics Company, Torrence, Ca, U.S.A.
Steeves, R.R., D. Delikaraoglou and N. Beck (1986). "Development of a Canadian Active Control System using GPS." Proceedings of the Fourth International Geodetic Symposium on Satellite Positioning, University of Texas at Austin, Austin, Tx, U.S.A., April, pp. 1189-1203. Stiffler, JJ. (1966). telecommunications." Vol. V of Space Technology, NASA. Stiller, A. (1985). "Development of civilian GPS receivers for the Federal Republic of Germany for different applications." Deutsche Forschungs- und Versuchsanstalt fur Luft und Raumfahrt, E.V. Sturza, M.A. (1984). "Commercial aviation GPS navigation set architecture." The National Technical Meeting of The (U.S.) Institute of Navigation, San Diego, Ca, U.S.A., January.
Swift, E.R. (1985). "NSWC's GPS orbit/clock determination system." Proceedings of the First International Symposium on Precise Positioning with the Global Positioning System, Rockville, Md, U.S.A., April, Vol. I, pp. 51-62.
Texas Instruments (1986). "Product specifications for TI-4100 Navstar Navigator." Lewisville, Ky, U.S.A,, January.
Bibliography D.I2
Thompson, S.D. (1985). Everyman's Guide to Satellite Navigation. ARINC Research Corp.
Thornton, C.L., W.G. Melbourne and T.H. Dixon (1986). "NASA GPS-based geodetic program in Mexico and the Caribbean." Proceedings of the Fourth International Geodetic Symposium on Satellite Positioning, University of Texas at Austin, Austin, Tx, U.S.A., April, pp. 1003-1018.
Tomioka, G. and K. Yamada (n.d.). "GPS navigator." Japan Radio Corproation, Tokyo, Japan. Trimble Navigation (1985). "Trimble 4000S GPS surveyor." Rev. 11/27/85, Sunnyvale, Ca, U.S.A. Trimble Navigation (1986). "GPS surveyor system." Sunnyvale, Ca, U.S.A., February.
U.S. Department of Defense/Department of Transportation (1982). Federal Radionavigation Plan. Vols. I-IV, March.
U.S. Department of Defense/Department of Transportation (1984). Federal Radionavigation Plan. Final report: March 1982-December 1984. U.S. Department of Defense document #DoD-4650.4 and U.S. Department of Transportation document #DoT-TSC-RSPA-84-8. Available from NTIS, Document AD-A151 295.
U.S. Federal Geodetic Control Committee (1986). "Proposed geometric geodetic survey standards and specifications for geodetic surveys using GPS relative positioning techniques." Version of 1 September, Rockville, Md, U.S.A.
University of Texas (1986). Proceedings of the Fourth International Geodetic Symposium on Satellite Positioning. University of Texas at Austin, Austin, Tx, U.S.A., April, 2 volumes.
van Dierendonck, A.J., S.S. Russell, E.R. Kopitzke and M. Birnbaum (1978). "The GPS navigation message." In: Global Positioning System. Papers published in Navigation, reprinted by The (U.S.) Institute of Navigation, Vol. I, 1980, pp. 55-73.
Vanicek, P. and E.J. Krakiwsky (1986). Geodesy: The Concepts. 2nd rev. ed., North-Holland, Amsterdam, The Netherlands.
Vanic'ek, P., R.B. Langley, D.E. Wells and D. Delikaraoglou (1984). "Geometrical aspects of differential GPS positioning." Bulletin Geodesique, 58, pp. 37-52.
Bibliography D.13
Vanicek, P., G. Beutler, A. Kleusberg, R.B. Langley, R. Santerre, and D.E. Wells (1985). "DIPOP differential positioning program package for the Global Positioning System." Department of Surveying Engineering Technical Report No. 115, University of New Brunswick, Fredericton, N.B., Canada.
Vantfek, P., D.E. Wells, A. Chrzanowski, A.C. Hamilton, R.B. Langley, J.D. McLaughlin and B.G. Nickerson (1983). "The future of geodetic networks." Proceedings of the International Association of Geodesy Symposia, Symposium d, The Future of Terrestrial and Space Methods for Positioning, Hamburg, FRG, August, Vol. 2, pp. 372-379. Vanicek, P., A. Kleusberg, R-G Chang, H. Fashir, N. Christou, M. Hofman, T. Kling and T. Arsenault (1986). "The Canadian geoid." Final contract report for the Geodetic Survey of Canada; Energy, Mines and Resources Canada, Ottawa, Ontario, Canada. Vessot, R.F.C. (1976). "Frequency and time standards." In: Methods of Experimental Physics, Vol. 12, Pt C, ed. M.L. Meeks, Academic Press, N.Y., NY, U.S.A. Wanless, B. (1985). "A prototype multi-station, multi-pass satellite data reduction program." M.Sc. thesis, Division of Surveying Engineering, The University of Calgary, Calgary, Alberta, Canada. Ward, P. (1981). "An inside view of pseudorange. and delta pseudorange measurements in a digital NAVSTAR GPS receiver." Presented at the ITC/USA/'81 International Telemetering Conference, San Diego, Ca, U.S.A., October. Wells, D.E. (1974). "Doppler satellite control." Department of Surveying Engineering Technical Report 29, University of New Brunswick, Fredericton, N.B., Canada.
Wells, D.E., P. Vanicek and D. Delikaraoglou (1981). "The application of NAVSTAR/GPS to geodesy in Canada: Pilot study." Department of Surveying Engineering Technical Report 76, University of New Brunswick, Fredericton, N.B., Canada.
Wells, D.E., A. Kleusberg, S.H. Quek, J. McCullough and J. Hagglund (1985). "Precise ship's velocity from GPS: Some test results." Proceedings of the First International Symposium on Precise Positioning with the Global Positioning System, Rockville, Md, U.S.A., April, Vol. II, pp. 799-808.
Bibliography D.14
Wong, R.V.C. and KP. Schwarz (1983). "Dynamic positioning with an integrated GPS-INS." Division of Surveying Engineering UCSE Report 30003, University of Calgary, Calgary, Alberta, Canada.
Yiu, K.P., R. Crawford and R. Eschenbach (1982). "A low-cost GPS receiver for land navigation." Navigation, Journal of The (U.S.) Institute of Navigation, Vol. 29, No. 3, pp. 204-220. Also in Global Positioning System. Papers published in Navigation, reprinted by The (U.S.) Institute of Navigation, Vol. II, 1984, pp. 44-60.
Ziemer, R.E. and R.L. Peterson (1985). Digital Communications and Spread Spectrum Systems. Macmillan Publishing Co., N.Y., NY, U.S.A.
ISBN:
0-920-114-73-3