Technical Article • July/August 2009
Making Sense of GPS InterSignal CorrectionsSatellite Calibration Parameters in Legacy and Modernized Ionosphere Correction AlgorithmsFigure 1 Because satellite designers cannot precisely eliminate small variations in delays among the various signal paths within the spacecraft, GPS signals don’t emerge from the satellite antennas at exactly the same time or from the same location. If not accounted for, these delay offsets would produce navigation fix errors for GPS users. In this article, the authors derive and explain the modernized ionosphere dual and singlefrequency correction algorithms that will resolve these intersignal variations. They also explain how to modernize any new algorithm, how a single ephemeris models the satellite location when each Lband signal may have separate departure points, and why that one ephemeris might point to a location behind the satellite.
Share via: Slashdot Technorati Twitter Facebook
The downloadable PDF (above) contains bonus material not available in the print edition. Appendices and other information on this subject is available at the bottom of the page. The Block IIRM satellites add the new L2C, L1M, and L2M signals on the existing Lband carriers. The Block IIF will also transmit these signals as well as a new L5 signal on the recently established L5 carrier. Signals/ranging/modulation codes all refer to the various GPS broadcast waveforms that can be present on the L1, L2, or L5 carriers. These codes perform three functions: code division multiple access (CDMA), processing/antijam gain, and indication of transmission time by providing each chip a unique timetag. The “new” algorithm, which we refer to as “the modernized ionospherefree pseudorange algorithm,” contains a mix of new parameters, intersignal corrections (ISCs), and the legacy scaled group delay differential parameter T_{GD}. This article describes why a new algorithm is needed, and how the new and old parameters are combined. With the U.S. Air Force set to begin broadcasting this fall the first of the new CNAV navigation messages on the L2 civil signal transmitted by the Block IIRM satellites, an understanding of these issues will be important for GPS receiver designers, GNSS signal simulator manufacturers, and end users who require highprecision results from their GPS equipment. In our presentation, we will first derive the new algorithm by using ISGPS200’s definition of delay error to model the GPS satellite equipment delays. Then, working through the various paragraphs on how to perform single and dualfrequency corrections, we can reconstruct the new equations by analyzing the modernized SV model. Having this modernized model allows one to also derive new results, such as a modernized version of the popular alternative ionosphere correction algorithm — the ionosphere pseudorange difference algorithm and understand why the new algorithm offers a more accurate correction.
Background If no other significant frequencydependent errors exist, two ranging measurements on different Lband carrier frequencies allow the ionosphere error to be eliminated. In general, until now only Department of Defense (DoD) receivers that track both the L1 P(Y) and the L2 P(Y) signals could make such a correction. (Some civilian receivers can track the DoD P(Y) code without a crypto key when the signal to noise ratio is sufficiently “high.”) With modernized GPS, however, civilian users will finally have access to a second frequency, the new L2C and/or L5 signals. For DoD users, the new L1M and L2M military signals will provide better cryptography as well as improved accuracy. Moreover, because a pair of M codes will reside on different carriers, dualfrequency ionosphere corrections capability can be maintained on modernized GPS military user equipment. Because we cannot precisely eliminate small variations in SV equipment delays among the various signal paths within the satellite, these new signals — as well as the legacy signals — don’t exactly emerge from the satellite antenna at the same time or from the same location. If not accounted for, these delay offsets would produce GPS navigation fix errors for dual and singlefrequency users. Consequently, ISGPS200 provides one ISC for each signal on each Lband. So, a Block IIRM SV with L1 CA, L1P(Y), L1M, L2C, L2P(Y), and L2M, has six ISCs. A GPS followon Block IIF satellite that also includes L5I and L5Q (in phase and in quadrature signals, respectively) will have two additional ISCs. The legacy T_{GD} parameter is really the ISC parameter for L2P(Y) scaled by a constant, which we will discuss later. We will use the following subscript notation when a pair of Lband measurements is being discussed: Li,x will denote one carrier frequency with signal x, and Lj,z will be used for the second Lband carrier frequency and the signals on that carrier. So for L1, Li=L1, and x=CA, P(Y), or M codes. For L2, Lj=L2, and z=L2C, P(Y), or M. For L5, Lj=L5, and z=L5I or L5Q. In all of the dual frequency correction algorithms, ISGPS200 and ISGPS705 use the symbol γ_{ij} to represent the ratio of (f_{Li}/f_{Lj})^{2} of the Lband carrier frequencies squared. If γ is not subscripted, one can assume that it is γ_{12}. To understand where the small SV equipment–delay errors come from, consider putting a picture to the words in paragraph 3.3.1.7 of ISGPS200 in order to build up a model of the modernized SV as shown in Figure 1 (above right). According to paragraph 3.3.1.7, the SVequipment group delay for each signal is the amount of time it takes the signal to start out from the common clock, travel through each code generator, modulator, transmitter, tri or quadraplexor, and finally emerge from the satellite’s antenna. Thus, the total delay consists of an electrical portion and an antenna portion. Note that the SV antenna is effectively a twoelement array, made up of an innerring that produces a broader beam and an outer ring that produces a narrower beam. As discussed in the articles by C. Choi and by G. Mader and F. Czopek listed in the Additional Resources section near the end of this article, the two rings are effectively phased 180 degrees apart so that the narrower beam is subtracted from the wider beam. As a result, the power across the surface of the Earth over the 13degree half angle of the transmitted signal is approximately constant. Because the SV antenna is an array, it can have gain, phase, and group delay variations across the beam, although for a properly designed antenna array, any angular variations of delay and phase would be small compared to the total signalinspace (SIS) error budget. Choi’s article notes the distinct locations of the phase centers for each Lband — all signals within an Lband sharing the same phase center, but it does not discuss the group delay characteristics of the SV antenna, nor does it discuss the possibility that each Lband signal may have its own distinct group delay center. As noted earlier, ISC_{Li,x} is the difference of the transit delay through the SV (including antenna) for L1P(Y) minus the transit delay for the x^{th} or z^{th} signal on the Li or Lj carrier. (It is important to point out that the ISC values, one for each signal, are not affected by the common clock error modeled by a second order polynomial with the coefficients a_{f0}, a_{f1}, and a_{f2} broadcast in subframe 1, as noted in the center of Figure 1, because this error is common to all signals and won’t contribute to delay differences.) The ISC values may age with time or change when redundant components are swapped in as the satellite ages. Each satellite has its own set of ISC values. The legacy parameter T_{GD} is really ISC_{L2P(Y)} scaled by (1γ). The ISCs for civil users are broadcast in Message #30 of the CNAV dataset, and the Mcode ISCs will be present in the MNAV dataset. The values are based on contractormeasured data. For T_{GD}, onorbit data is also used, as discussed in the article by B. Wilson et alia listed in Additional Resources.
Problem and Simplified Solution We should note that, strictly speaking, ISGPS200 defines the ISCs in terms of timetag differences, not delay differences, and it takes a substantial amount of effort and verbiage to prove that the delay notation we use is equivalent to the timetag notation used in reference 1. Also of interest is the fact that, historically, ISGPS200 was written for engineers who were building entire GPS receivers; thus, many of the compensation algorithms are stated in terms of the delay lock loop’s code phase readout of the time of transmission as measured by the receiver. Only later is the final pseudorange formed as the user time minus time of transmission scaled by the speed of light, c. For our full derivation of the new algorithm with all of the exhaustive details, readers can download an extended version of this article using the PDF link at the top of this story. Although ISGPS200 can be tedious to read, a tremendous number of insights can be gleamed from it. (As an aside, it would be nice if a team of authors could write a companion document to ISGPS200 that describes why certain decisions were made, and to provide derivations, so that valuable knowledge will not be lost as people retire.) For the complete, extended version of this story, including figures, equations, and the following topics, please download the PDF of the article, above.
Acknowlegments
Additional Resources Additional Information» julyaug09tetewskyappendices.pdf Copyright © 2017 Gibbons Media & Research LLC, all rights reserved. 
