The loose combination (LC) and the tight combination (TC) are two different models in the combined processing of four global navigation satellite systems (GNSSs). The former is easy to implement but may be unusable with few satellites, while the latter should cope with the inter-system bias (ISB) and is applicable for few tracked satellites. Furthermore, in both models, the inter-frequency bias (IFB) in the GLObal NAvigation Satellite System (GLONASS) system should also be removed. In this study, we aimed to investigate the performance difference of ambiguity resolution and position estimation between these two models simultaneously using the single-frequency data of all four systems (GPS + GLONASS + Galileo + BeiDou Navigation Satellite System (BDS)) in three different environments, i.e., in an open area, with surrounding high buildings, and under a block of high buildings. For this purpose, we first provide the definition of ISB and IFB from the perspective of the hardware delays, and then propose practical algorithms to estimate the IFB rate and ISB. Thereafter, a comprehensive performance comparison was made between the TC and LC models. Experiments were conducted to simulate the above three observation environments: the typical situation and situations suffering from signal obstruction with high elevation angles and limited azimuths, respectively. The results show that in a typical situation, the TC and LC models achieve a similar performance. However, when the satellite signals are severely obstructed and few satellites are tracked, the float solution and ambiguity fixing rates in the LC model are dramatically decreased, while in the TC model, there are only minor declines and the difference in the ambiguity fixing rates can be as large as 30%. The correctly fixed ambiguity rates in the TC model also had an improvement of around 10%. Once the ambiguity was fixed, both models achieved a similar positioning accuracy.

The multipath effect is a crucial error source caused by the environment around the station and cannot be eliminated or mitigated by differential algorithms. Theoretically, the maximum value for the carrier phase is a quarter the wavelength, i.e., about 4.8 cm for the GPS L1 signal. Considering the increasing demands of high-precision applications, the multipath error has become a major factor affecting the accuracy and reliability of GPS millimeter-level data processing. This paper proposes a multi-point hemispherical grid model (MHGM) to mitigate the multipath effect. In this method, the hemisphere centered on each station is divided into a grid, and the multipath error at the station is estimated based on the parameterization of the grid points. The double-differenced (DD) observed-minus-calculated (OMC) values on some previous days are treated as the observation values to model the present multipath error. Contrary to the present methods which rely much on the platform of data collection and processing, MHGM can be potentially applied to GPS data processing with the existing hardware and software. Experiments in high-multipath and low-multipath environments are designed by mounting a baffle or not. The experimental results show that MHGM is effective in mitigating the multipath effect. When using data from the previous day, an average improvement of about 63.3% in the RMS of DD OMC can be made compared with that without correction, and this is basically consistent with the sidereal filtering (SF) method which is 63.0%. Furthermore, the effectiveness of the above two methods is better than that of the empirical site model (ESM). The kinematic positioning results are also basically consistent with the statistical results of the RMS values of DD OMC. Historical data from more than one day can more explicitly and effectively model the MHGM. Furthermore, compared with the SF method, the MHGM can be used not only to mitigate the multipath error, but also to orientate the sources of the multipath error around the station, and give guidance in the physical elimination of these sources.

The multipath error is considered to be the most limiting factor for high precision positioning applications. The sidereal filtering (SF) method can be used to mitigate the multipath error in the observation domain, and it has been successfully applied in the multipath mitigation in global positioning systems (GPS) and regional BeiDou navigation satellite systems (BDS2). However, there are few reports on the SF method in other systems. The performance of the SF method relies on the explicit orbit repeat periods of satellites in diverse systems or even different types of constellations. It is therefore inconvenient to utilize the SF method for multi-GNSS multipath error mitigation. Alternatively, a space domain multipath error reduction method, which establishes the multi-point hemispherical grid model (MHGM) using the residuals of the double-differenced carrier phase observations in the ambiguity-fixed period, has been modified. It is an integrated model for multi-GNSS, without considering the diversity of different systems and constellations. To compare the performance of MHGM and SF from a multi-GNSS point of view, the determination method of orbit repeat periods via the broadcast ephemerides is summarized, and the SF method is extended to the global BeiDou navigation satellite system (BDS3) and Galileo navigation satellite system. Further test results show that the performance of MHGM and SF are comparable from the perspective of root mean squares (RMS) and the power spectrum analysis of double-differenced residuals, as well as the static positioning results. This implies that the space domain MHGM can obtain similar correction effects as the SF method in the observation domain, but the former is more flexible for modeling with various systems’ data. In addition, the established MHGM using the data of multi orbit periods demonstrates a better performance compared with that of only one orbit period, and an average improvement of 13.1% in the RMS of the double-differenced residuals can be achieved.

Considering the influence of the ionosphere, troposphere, and other systematic errors on double-differenced ambiguity resolution (AR), we present an optimal triple-frequency code-phase combination determination method driven by both the model and the real data. The new method makes full use of triple-frequency code measurements (especially the low-noise of the code on the B3 signal) to minimize the total noise level and achieve the largest AR success rate (model-driven) under different ionosphere residual situations (data-driven), thus speeding up the AR by directly rounding. With the triple-frequency Beidou Navigation Satellite System (BDS) data collected at five stations from a continuously-operating reference station network in Guangdong Province of China, different testing scenarios are defined (a medium baseline, whose distance is between 20 km and 50 km; a medium-long baseline, whose distance is between 50 km and 100 km; and a long baseline, whose distance is larger than 100 km). The efficiency of the optimal code-phase combination on the AR success rate was compared with that of the geometry-free and ionosphere-free (GIF) combination and the Hatch-Melbourne-Wübbena (HMW) combination. Results show that the optimal combinations can always achieve better results than the HMW combination with B2 and B3 signals, especially when the satellite elevation angle is larger than 45°. For the wide-lane AR which aims to obtain decimeter-level kinematic positioning service, the standard deviation (STD) of ambiguity residuals for the suboptimal combination are only about 0.2 cycles, and the AR success rate by directly rounding can be up to 99%. Compared with the HMW combinations using B1 and B2 signals and using B1 and B3 signals, the suboptimal combination achieves the best results in all baselines, with an overall improvement of about 40% and 20%, respectively. Additionally, the STD difference between the optimal and the GIF code-phase combinations decreases as the baseline length increases. This indicates that the GIF combination is more suitable for long baselines. The proposed optimal code-phase combination determination method can be applied to other multi-frequency global navigation satellite systems, such as new-generation BDS, Galileo, and modernized GPS.