Research Article  Open Access
Patrick Henkel, Christoph Günther, "Joint L/CBand Code and Carrier Phase Linear Combinations for Galileo", International Journal of Navigation and Observation, vol. 2008, Article ID 651437, 8 pages, 2008. https://doi.org/10.1155/2008/651437
Joint L/CBand Code and Carrier Phase Linear Combinations for Galileo
Abstract
Linear code combinations have been considered for suppressing the ionospheric error. In the Lband, this leads to an increased noise floor. In a combined L and Cband (5010–5030 MHz) approach, the ionosphere can be eliminated and the noise floor reduced at the same time. Furthermore, combinations that involve both code and carrierphase measurements are considered. A new Lband codecarrier combination with a wavelength of 3.215 meters and a noise level of 3.92 centimeters is found. The double difference integer ambiguities of this combination can be resolved by extending the system of equations with an ionospherefree L/Cband code combination. The probability of wrong fixing is reduced by several orders of magnitude when Cband measurements are included. Carrier smoothing can be used to further reduce the residual variance of the solution. The standard deviation is reduced by a factor 7.7 if Cband measurements are taken into account. These initial findings suggest that the combined use of L and Cband measurements, as well as the combined code and phase processing are an attractive option for precise positioning.
1. Introduction
The integer ambiguity resolution of carrierphase measurements has been simplified by the consideration of linear combinations of measurements at multiple frequencies. Early methods were the threecarrier ambiguity resolution (TCAR) method introduced by Forssell et al. [1], as well as the cascade integer resolution (CIR) developed by Jung et al. [2]. The weighting coefficients of threefrequency phase combinations are designed either to eliminate the ionosphere at the price of a rather small wavelength or to reduce the ionosphere only by a certain amount with the advantage of a larger wavelength.
The systematic search of all possible GPS L1L2 widelane combinations has been performed by Cocard and Geiger [3]. The L1L2 linear combination of maximum wavelength ( m) amplifies the ionospheric error by a factor . Collins gives an overview of reduced ionosphere L1L2 combinations with wavelengths up to cm ( widelane) in [4].
The authors have extended this work to threefrequency (3F) Galileo combinations (E1E5aE5b) in [5]. A 3F widelane combination with a wavelength of m and an ionospheric suppression of dB was found. Furthermore, a 3F narrowlane combination with a wavelength of cm could reduce the ionospheric error by as much as dB. Sets of linear carrierphase combinations that are robust against residual biases were studied in [6]. The integer ambiguities of the linear combinations can be estimated by the leastsquares ambiguity decorrelation adjustment (LAMBDA) algorithm developed by Teunissen [7]. The method includes an integer transformation which can also be used to determine optimum sets of linear combinations [8].
In this paper, the authors used code and carrierphase measurements in the linear combinations for obtaining ionospheric elimination, large wavelengths, and a low noise level at the same time. The E5a and E5b code measurements are of special interest due to their large bandwidth ( MHz) and their low associated CramerRao bound of cm [9]. The Cband phase measurements are particularly interesting due to their small wavelength and their thus reduced phase noise. The properties of codecarrier linear combinations are optimized by including both L and Cband measurements. The cost function is defined as the ratio of half the wavelength and the noise level of the ionospherefree codecarrier combination. It is called combination discrimination and it is a measure of the radius of the decision regions expressed in units given by the standard deviation of the noise. The LBand signals of Galileo are defined in the GalileoICD [10]. The Cband signals are foreseen in a band between MHz to MHz [11]. The signal propagation and tracking characteristics in the Cband have been analyzed by Irsigler et al. [12]. The larger frequencies result in an additional free space loss of dB that has to be compensated by a larger transmit power.
The paper is organized as follows: the next section introduces the design of codecarrier linear combinations. The underlying tradeoff between a low noise level and strong ionospheric reduction turns out to be controlled by the weighting coefficients of E5a/E5b code measurements.
In Section 3, codecarrier linear combinations are computed in a way that include both L and Cband measurements. An ionospherefree codeonly combination is determined that benefits from a 4.5 times lower noise level than a pure Lband combination. Furthermore, a pure Lband ionospherefree codecarrier combination with a wavelength of 3.215 m and a noise standard deviation of 3.9 cm is found. The combined use of the two reduces the probability of wrong fixing of the latter solution by 9 orders of magnitude with respect to a pure Lband solution.
The use of Cband measurements for ionospherefree carrier smoothing is discussed in Section 4: an ionospherefree codecarrier combination of arbitrary wavelength is smoothed by a pure phase combination. The low noise level of Cband measurements provides a linear combination that benefits from an dB lower noise level as compared to the equivalent Lband combination.
2. CodeCarrier Linear Combinations
Linear combinations of carrierphase measurements are constructed to increase the wavelength (widelane), suppress the ionospheric error, and to simplify the integer ambiguity resolution. The properties of the linear combinations can be improved by including weighted code measurements into the pure phase combinations. Figure 1 shows a three frequency (3F) linear combination where the phase measurements are weighted by , , , and the code measurements are scaled by , , . The weighting coefficients are generally restricted by a few conditions: first, the geometry should be preserved, that is Moreover, the superposition of ambiguities should be an integer multiple of a common wavelength , that is which can be split into three sufficient conditionswith denoting the space of integers. These integer constraints are rewritten to obtain the weighting coefficients
Mixed codecarrier combinations weight the phase part by and the code part by . The border cases are pure phase () and pure code () combinations. The parameter has a significant impact on the properties of the linear combination, and it is optimized later in this section. Replacing the weighting coefficients in by (4) yields the wavelength of the codecarrier combinationThe generalized widelane criterion is given for by . Equivalently, it can be expressed as a function of , and asThe linear combination scales the ionospheric error byThe thermal noise of the elementary carrier phase measurements is assumed Gaussian with the standard deviation given by Kaplan and Hegarty [13]where denotes the loop bandwidth, the carriertonoise ratio, and the predetection integration time. The overall noise contribution of the linear combination is written aswith being the noise variance of the code measurements. Table 1 shows the CramerRao bound (CRB) for some Galileo signals as derived by Hein et al. [9]. A DLL bandwidth of Hz has been assumed. The 4 MHz receiver bandwidth for E1 has been chosen to avoid sidelobe tracking.

For E1, E5a, E5b, E6 phase measurements, the wavelength scaling of can be neglected due to the close vicinity of the frequency bands. However, it plays a major role when CBand measurements are included.
Figure 2 shows the benefit of the code contribution to the (E1), (E5b), (E5a) linear combination: a slight increase in noise level results in a considerable reduction of the ionospheric error. The phase weighting has been fixed to so that , , , and are uniquely determined. The E5b and E5a code weights are adapted continuously and the ionosphere is eliminated in the border caseE1 code measurements have not been taken into account due to the increased noise level but might be included with a low weight.
The combination discrimination—measured by the ratio of half the wavelength and the noise level —is proposed as a cost function to select linear combinations due to its independence of the geometry. It is shown for multiple ionospherefree codecarrier combinations in Figure 3. The strong dependency on the phase weighting suggests an optimization with respect to this parameter. Note that the legend refers to the elementary wavelengths which have to be scaled by .
The computation of the optimum takes again only E5a/ E5b code measurements into account as the benefit of the E1 code measurement is negligible (). The notation is simplified bywith , and implicitly given by (5), (7), and (9). The E5a/E5b code weights are determined from the ionospherefree and geometrypreserving conditions as
The combination discrimination becomes from (5), (11), and (12): Setting the derivative to zero yields the optimum weightingwhich is independent of both and . Table 2 contains the weighting coefficients and characteristics of the codecarrier combinations shown in Figure 3.

Figure 4 shows the benefit of adaptive code and phase weighting for the codecarrier linear combination with (E1), (E5b), and (E5a). Obviously, the wavelength increases linearly with and the ionosphere can be eliminated with any . The noise amplification depends on the level of ionospheric reduction: a linear increase can be observed near the pure phase combination, while the increase becomes negligible for the ionospherefree combination. Thus, the combination discrimination of the ionospherefree codecarrier combination is increased by almost the same factor as is risen.
3. CBand Aided CodeCarrier Linear Combinations
The MHz wide CBand () has been reserved for Galileo. The higher frequency range has a multitude of advantages and drawbacks: an additional free space loss of dB occurs which has to be compensated by a larger transmit power. The ionospheric delay is approximately times lower than in the E1 band. The small wavelength of complicates direct ambiguity resolution but results in an approximately times lower standard deviation of phase noise. Moreover, the CBand offers additional degrees of freedom for the design of linear combinations.
3.1. Reduced Noise IonosphereFree CodeOnly Combinations
The design of three frequency codeonly combinations that preserve geometry and eliminate ionospheric errors is characterized by one degree of freedom used for noise minimization. The weighting coefficients are derived from the geometry preserving and ionospherefree constraints in (1), (7) asMinimization of yieldsIonospherefree codeonly combinations with more than three frequencies are obtained by a multidimensional derivative. Table 3 shows that the pure Lband E1E5bE5a combination is characterized by a noise level of . If the E5 signal is received with full bandwidth, the CRB is reduced to cm but the number of degrees of freedom is reduced by one so that the noise level of the E1E5 combinations are lightly increased (Table 4). These combinations will play a role in conjunction with codecarrier combinations.


The Cband is split into bands of MHz bandwidth centered at and allows a significant reduction of the noise level. Note that the noise of any contributing elementary combination is reduced by a weighting coefficient smaller than one.
3.2. Joint L/CBand Widelane Combinations
Codecarrier linear combinations can also include both Lband and Cband measurements. Therefore, (1)–(7) are extended to include the additional measurements. The weighting coefficients and are computed such that the discriminator output of (13) is maximized for a given set of integer coefficients .
Table 5 contains ionospherefree joint L/Cband codecarrier widelane combinations The E1E5 pure Lband combination benefits from a noise level of only which simplifies the resolution of the integer ambiguities. In contrast to the codeonly combinations, the use of the fullbandwidth E5 signal is advantageous compared to separate E5a and E5b measurements. The Cband offers no benefit for these wavelengths. In the last row of Table 5, a linear combination with a pure Lband code and pure Cband phase part is described. The combination discrimination equals but the noise level is also increased to .

Figure 5 shows the tradeoff between wavelength and noise level for joint L/Cband ionospherefree linear codecarrier combinations with and . The E1E5 combination is of special interest but the maximum combination discrimination is obtained for a joint L/Cband combination.
3.3. Joint L/CBand Narrowlane Combinations
There exists a large variety of joint codecarrier narrowlane combinations where Cband measurements help to reduce the noise substantially. Figure 6 shows the tradeoff between wavelength and noise level for and . For , the consideration of Cband measurements reduces the noise level by a factor of compared to a pure Lband combination (Table 6).

3.4. Reliability of Ambiguity Resolution
The integer ambiguity resolution is based on the linear combination of four different variable types: doubledifference measurements for eliminating clock errors and satellite/receiver biases; multifrequency combinations for suppressing the ionosphere; code and carrier phase measurements for reducing the noise level; and finally, L/Cband combinations for noise and discrimination characteristics.
Two joint L/Cband codecarrier ionospherefree combinations are chosen for realtime (single epoch) ambiguity resolution. The , combination of Table 5 and one further combination of Table 4. The double difference (DD) ionospherefree combinations are modeled aswithand the DD geometry matrix , the baseline and the integer ambiguities . The doubledifferenced troposphere is assumed to be negligible or known a priori (e.g., from an accurate continued fraction model).
Note that the troposphere has the same impact on all geometrypreserving combinations and does not affect the optimization of the mixed codecarrier combinations. The noise vector is Gaussian distributed, that is,where denotes the Kronecker product. models the linear combination induced correlation and includes the correlation due to double difference measurements from visible satellites. The standard deviation of the most critical ambiguity estimate can be written asand the probability of wrong fixing follows asIn the following analysis, the location of the reference station is at N, E with a baseline length of km.
Figure 7 shows the benefit of Cband measurements for integer ambiguity fixing. If the E1E5 pure Lband combination is used as second combination in (17), the failure rate varies between and due to its poor noise characteristics. The use of two additional Cband measurements reduces the maximum probability of wrong fixing to . For three Cband frequencies, the failure rate is at most which corresponds to a gain of to orders of magnitude compared to the pure Lband combination.
The reliability of ambiguity resolution can be further improved by using the LAMBDA method of Teunissen [7]. The float ambiguity estimates are decorrelated by an integer transformation and the ambiguity covariance matrix is written as with the decomposition into a lower triangular matrix and a diagonal matrix . The probability of wrong fixing of the sequential bootstrapping estimator is given by Teunissen [14] aswith . It represents a lower bound for the success rate of the integer leastsquare estimator and is depicted in Figure 8. Obviously, the use of joint L/Cband linear combinations reduces the probability of wrong fixing by several orders of magnitude compared to pure Lband combinations.
3.5. Accuracy of Baseline Estimation
After integer ambiguity fixing, the baseline is reestimated from (17). The covariance matrix of the baseline estimate in local coordinates is given bywith the rotation matrix . Figure 9 shows the achievable horizontal and vertical accuracies for the two optimized joint L/Cband combinations.
The pure Lband combinations in the first row of Tables 4 and 5 have been again selected as reference scenario. It can be observed that the use of joint L/Cband linear combinations enables a slight improvement in position estimates compared to the significant benefit for ambiguity resolution.
4. Joint L/CBand Carrier Smoothed Carrier
Ionospherefree codecarrier linear combinations are characterized by a noise level that is one to two orders of magnitude larger than of the underlying carrierphase measurements (Table 5). Both noise and multipath of the codecarrier combinations can be reduced by the smoothing filter of Hatch [15] which is shown in Figure 10. The upper input can be an ionospherefree codecarrier combination of arbitrary wavelength. The lower input is a pure ionospherefree phase combination that is determined by three conditions: the first ensures that the geometry is preserved, the second eliminates the ionosphere, and the third minimizes the noise, that is,
Note that the superposition of ambiguities of the pure phase combination is not necessarily an integer number of a common wavelength. The respective ambiguities are not affected by the low pass filter and do not occur in the smoothed output due to different signs in the addition to (Figure 10).
Table 7 shows an ionospherefree E1E5aE5b phase combination that increases the noise level by a factor . However, the low noise level of Cband measurements suggests the use of the second combination with . In this case, the noise level is not only reduced by smoothing but also by the coefficients of the pure phase combination.

The variance of the smoothed combination is given bywith the lowpass filtered noise (e.g., Konno et al. [16])and the smoothing time . Setting (27) into (26), and using the definition of a geometric series yields For long smoothing times, only the low noise of the joint L/C pure carrierphase combination remains (Table 7).
5. Conclusions
In this paper, new joint L/Cband linear combinations that include both code and carrierphase measurements have been determined. The weighting coefficients are selected such that the ratio between wavelength and noise level is maximized. An ionospherefree Lband combination (IFL) could be found at a wavelength of m with a noise level of cm.
The combination of L and Cband measurements reduces the noise level of ionospherefree codeonly combinations by a factor compared to pure Lband combinations. This increases the reliability of an ambiguity resolution option for the IFL combination by 9 orders of magnitude.
The residual variance of the noise can be further reduced by smoothing. An L/Cband carrier combination can smooth the noise with a residual variance below the Lband phase noise variance. The smoothed solution can either be used directly or can be used to resolve the narrowlane ambiguities. The variance is basically the same in both cases. The resolved ambiguities, however, provide instantaneous independent solutions.
References
 B. Forssell, M. MartinNeira, and R. A. Harris, “Carrier phase ambiguity resolution in GNSS2,” in Proceedings of the 10th International Technical Meeting of the Satellite Division of the Institute of Navigation (ION GPS '97), vol. 2, pp. 1727–1736, Kansas City, Mo, USA, September 1997. View at: Google Scholar
 J. Jung, P. Enge, and B. Pervan, “Optimization of cascade integer resolution with three civil frequencies,” in Proceedings of the 13th International Technical Meeting of the Satellite Division of the Institute of Navigation (ION GPS '00), Salt Lake City, Utah, USA, September 2000. View at: Google Scholar
 M. Cocard and A. Geiger, “Systematic search for all possible Widelanes,” in Proceedings of the 6th International Geodetic Symposium on Satellite Positioning, Columbus, Ohio, USA, March 1992. View at: Google Scholar
 P. Collins, “An overview of gps interfrequency carrier phase combinations,” Technical Memorandum, Geodetic Survey Division, University of New Brunswick, Ottawa, Ontario, Canada, 1999. View at: Google Scholar
 P. Henkel and C. Günther, “Three frequency linear combinations for Galileo,” in Proceedings of the 4th Workshop on Positioning, Navigation and Communication (WPNC '07), pp. 239–245, Hannover, Germany, March 2007. View at: Google Scholar
 P. Henkel and C. Günther, “Sets of robust fullrank linear combinations for widearea differential ambiguity fixing,” in Proceedings of the 2nd ESA Workshop on GNSS Signals and Signal Processing, Noordwijk, The Netherlands, April 2007. View at: Google Scholar
 P. Teunissen, “Leastsquares estimation of the integer ambiguities,” 1993, Invited lecture, Section IV, Theory and Methodology, IAG General Meeting, Beijing, China. View at: Google Scholar
 P. J. G. Teunissen, “On the GPS widelane and its decorrelating property,” Journal of Geodesy, vol. 71, no. 9, pp. 577–587, 1997. View at: Publisher Site  Google Scholar
 G. Hein, J. Godet, J. Issler et al., “Status of galileo frequency and signal design,” in Proceedings of the 15th International Technical Meeting of the Satellite Division of the Institute of Navigation (ION GPS '02), pp. 266–277, Portland, Ore, USA, September 2002. View at: Google Scholar
 “Galileo Open Service SignalinSpace ICD,” http://www.galileoju.com. View at: Google Scholar
 European Radiocommunications Office, “The European Table of Frequency Allocations and Utilisations in the Frequency Range 9 kHz to 1000 GHz,” 2007, ERC Report 25, p. 132. View at: Google Scholar
 M. Irsigler, G. Hein, B. Eissfeller et al., “Aspects of Cband satellite navigation: signal propagation and satellite signal tracking,” in Proceedings of the European Navigation Conference (ENC GNSS '02), Kopenhagen, Denmark, May 2002. View at: Google Scholar
 E. Kaplan and C. Hegarty, Understanding GPS: Principles and Applications, Artech House, London, UK, 2nd edition, 2006.
 P. J. G. Teunissen, “Success probability of integer GPS ambiguity rounding and bootstrapping,” Journal of Geodesy, vol. 72, no. 10, pp. 606–612, 1998. View at: Publisher Site  Google Scholar
 R. R. Hatch, “A new threefrequency, geometryfree, technique for ambiguity resolution,” in Proceedings of the 19th International Technical Meeting of the Satellite Division of the Institute of Navigation (ION GNSS '06), vol. 1, pp. 309–316, Fort Worth, Tex, USA, September 2006. View at: Google Scholar
 H. Konno, S. Pullen, J. Rife, and P. Enge, “Evaluation of two types of dualfrequency differential GPS techniques under anomalous ionosphere conditions,” in Proceedings of the National Technical Meeting of the Institute of Navigation (NTM '06), vol. 2, pp. 735–747, Monterey, Calif, USA, January 2006. View at: Google Scholar
Copyright
Copyright © 2008 Patrick Henkel and Christoph Günther. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.