Performance analysis of frequency-mixed PPP-RTK using low-cost GNSS chipset with different antenna configurations

Low-cost Global Navigation Satellite System (GNSS) devices offer a cost-effective alternative to traditional GNSS systems, making GNSS technology accessible to a wider range of applications. Nevertheless, low-cost GNSS devices often face the challenges in effectively capturing and tracking satellite signals, which leads to losing the observations at certain frequencies. Moreover, the observation peculiarities of low-cost devices are in contradistinction to those of traditional geodetic GNSS receivers. In this contribution, a low-cost PPP-RTK model that considers the unique characteristics of different types of measurements is developed and its performance is fully evaluated with u-blox F9P receivers equipped with three distinctive antenna configurations: vertical dipole, microstrip patch, and helix antennas. Several static and kinematic experiments in different scenarios are conducted to verify the effectiveness of the proposed method. The results indicate that the mixed-frequency PPP-RTK model outperforms the traditional dual-frequency one with higher positioning accuracy and fixing percentage. Among the three low-cost antennas tested, the vertical dipole antenna demonstrates the best performance under static conditions and shows a comparable performance as geodetic antennas with a positioning accuracy of 0.02 m, 0.01 m and 0.07 m in the east, north, and up components, respectively. Under low-speed kinematic scenarios, the helix antenna outperforms the other two with a positioning accuracy of (0.07 m, 0.07 m, 0.34 m). Furthermore, the helix antenna is also proved to be the best choice for vehicle navigation with an ambiguity fixing rate of over 95% and a positioning accuracy of (0.13 m, 0.14 m, 0.36 m).


Introduction
Global Navigation Satellite System (GNSS) has been widely used in geodesy relying on geodetic receivers, which exhibit commendable capabilities in terms of multipath suppression, noise reduction, and observation quality.However, the widespread adoption of geodetic receivers is hindered by their high cost, bulky size, substantial power consumption, and inconvenient portability in the mass market.Thanks to the rapid advancements in manufacturing technology, GNSS receivers are undergoing a transformative shift towards miniaturization, lightweight design, and affordability (Amami et al., 2014;Hofmann-Wellenhof et al., 2012).Many commercial companies such as Broadcom (www. broad com. com), u-blox (www.u-blox. com), and unicorecomm (www. unico recomm. com) have developed low-cost GNSS chipsets for as little as 100-200 dollars.This has made GNSS with low-cost receivers be a new research hotspot and of great potential for emerging massive and automatic applications (Mongrédien et al., 2016;Odolinski & Teunissen, 2016, 2017).
Numerous studies investigated the observation characteristics of low-cost receivers, including miniaturized GNSS boards like u-blox, entry-level NovAtel devices, and smartphones (Caldera et al., 2016;Chen & Qin, 2012;Li & Geng, 2019;Wanninger & Heßelbarth, 2020).The research indicated that the low-cost devices had higher noise level and lower carrier-to-noise-density ratio compared to the geodetic ones.Additionally, low-cost devices are characterized by severe systematic errors, such as the multipath (Gogoi et al., 2018;Zhang et al., 2019b), interfrequency bias (Håkansson et al., 2017), and inter-system bias (Mi et al., 2019).To better access these systematic errors, Zhang et al. (2021) developed an easily implementable method based on the ionosphere-corrected and geometry-free model, which is effective regardless the number of satellites and frequencies in most scenarios.Moreover, the quality of observations greatly relies on the receiving antenna.Cina and Piras (2015) reported that the combination of low-cost GNSS receivers with a geodetic antenna could improve their performance.Hamza et al. (2021) and Amami (2022) tested different types of antennas and illustrated their advantages and limitations.However, their studies focused mainly on the assessment of noise, and more detailed analysis is required.
As to positioning capability, several studies were conducted to evaluate the performance of low-cost GNSS receivers.Takasu andYasuda (2008, 2009) were the pioneers in evaluating the Real-Time Kinematic (RTK) performance and demonstrated that RTK method using the low-cost GNSS receivers could achieve comparable performances with geodetic grade receivers under favorable circumstances.The similar conclusions were drawn by Jo et al. (2013), Cina and Piras (2015), and Garrido-Carretero et al. (2019).Unlike the RTK algorithm, Precise Point Positioning (PPP) technology can achieve centimeter-level positioning accuracy using only a single GNSS receiver.Gill et al. (2017) assessed the single-frequency PPP performance with the ionospheric products provided by the Centre National d'Etudes Spatiales (CNES).Following this, the ionosphere-weighted PPP model (Li et al., 2019) to identify various situations with the different uncertainties of ionospheric constraints was addressed.Wen et al. (2020) and Wang et al. (2021) utilized the observations from smart devices to conduct PPP.Furthermore, Li and Geng (2022) enhanced the ambiguity resolution method for Android data.Considering the coexistence of single-and dual-frequency observations on low-cost GNSS devices, Nie et al. (2020) employed the singlefrequency ionosphere-corrected code measurements and the conventional dual-frequency ionosphere-free code and phase measurements, shorterning PPP convergence time to a few minutes.
Both RTK and PPP technologies using low-cost devices have their own pros and cons.To achieve rapid and accurate positioning for a wide range of users, researchers proposed a combined approach called PPP-RTK (Li et al., 2011;Teunissen et al., 2010;Wübbena et al., 2005).The PPP-RTK technique has the merits of flexibility, high accuracy, extensive coverage, and enhanced privacy protection, making it as a promising tool with innovative opportunities for emerging mass applications.Notably, scholars have explored PPP-RTK using low-cost devices.Nadarajah et al. (2018) investigated the PPP-RTK performance with low-cost single-frequency U-blox M8 receivers and found that 5cm accuracy could be achieved within 9 min using the data with a 30-s sampling interval.Li et al. (2022) demonstrated the feasibility of PPP-RTK using android GNSS raw measurements with an external antenna, revealing the potential for instantaneous decimeter-level positioning accuracy under low-speed dynamic conditions.However, there are still some challenges in PPP-RTK with low-cost devices.Firstly, data loss often occurs due to the hardware constraints of low-cost devices.Effectively utilizing all available data becomes crucial.Secondly, the disparities in observation characteristics between low-cost and high-end devices render the traditional stochastic model inadequate, necessitating the determination of appropriate weight ratios.Additionally, there are various types of low-cost antennas in market, such as helix antennas, microstrip patch antennas, and vertical dipole antennas, which have different performances depending on the manufacturing process and physical structure (Amami et al., 2014;Bancroft, 2019).Therefore, further exploration of low-cost devices' PPP-RTK performance is urgently needed.
In this contribution, we implement PPP-RTK with lowcost devices and evaluate their performance.A mixedfrequency PPP-RTK model and a modified stochastic model considering the prior knowledge of different types of measurements are proposed.In addition, the signal characteristics and positioning performance of different low-cost devices are compared and analyzed.After this introduction, the PPP-RTK system with low-cost devices is first described.Then the experimental sets and the processing strategies are detailed.Hereafter, the effectiveness of the method proposed in this study and the performance of different low-cost devices are evaluated.Finally, conclusions are summarized.

Methodology
In this section, the implementation of the mixed-frequency PPP-RTK model and the cascade ambiguity resolution strategy are introduced first.A stochastic model considering the prior information on low-cost GNSS devices is then detailed.

Mixed-frequency PPP-RTK model
The performances of low-cost GNSS devices are compromised in some functions, such as preamplifier, Digital Signal Processing (DSP), and Low Noise Amplifier (LNA), resulting in the loss of observations, particularly on the second band in certain periods.To address this limitation, this paper proposes a mixed-frequency PPP-RTK model that leverages all available observation information to improve positioning performance.As illustrated in Fig. 1, the model comprises two main components: the server and the client.The former provides the precise atmospheric corrections derived from ambiguity-fixed PPP solutions with the precise orbit, clock, and Uncalibrated where ρ s r,i is the geometric distance between the receiver r and the satellite s; t r and t s denote clock offset at the receiver and satellite side respectively in seconds scaled by the light speed in vacuum c ; I s r,1 refers to the slant ionospheric delay at the first frequency linked to other frequencies by γ i = 2 i / 2 1 with the wavelength i ; the zenith tropospheric delay Z r is mapped to the slant direc- tion by the function m s r ; the integer ambiguity N s r,i is in cycles scaled by the wavelength i ; b r,i and B r,i denote the receiver hardware biases in pseudorange and carrier phase observations, respectively, while b s i and B s i are satellite counterparts; e s r,i and ε s r,i represent measurement noise of pseudorange and carrier phase, respectively.Other errors such as the satellite and receiver antenna Phase Center Offsets (PCOs) and Phase Center Variations (PCVs), the relativistic effects, the tidal loadings, and the phase wind-up effect are corrected by the existing models (Kouba, 2009).
Considering the coexistence of single-, dual-, multifrequency GNSS observations using low-cost devices, the uncombined PPP-RTK model with the atmospheric constraints can be expressed as: (1) Phase Delays (UPD) products (Li et al., 2021;Psychas & Verhagen, 2020).The latter establishes an uncombined observation model (Geng et al., 2020;Zhang et al., 2019a) that incorporates single-, dual-, and multi-frequency observations.The rapid ambiguity resolution in the mixed-frequency mode is then performed with the assistance of the atmospheric corrections provided by the server component.For a receiver r which tracks a satellite s of system sys at frequency i, the GNSS pseudorange measurement ( P s r,i ) and carrier phase measurement ( L s r,i ) can be described as: (3) where E(•) denotes the expectation operator; p s r,i and l s r,i represent the observed-minus-computed code and phase observations after the removal of the atmospheric corrections f aug P s r,i and f aug L s r,i , respectively; H (the representative of p s r , l s r , N s r , . . . ) is a column vec- tor whose number of rows depends on the number of frequencies k; X u and L u represent the user coordinates to be estimated and the design matrix after linearization, respectively; I V is a column vector with all elements of 1; I and O are identity matrix and zero matrix, respec- tively.The satellite clock offsets t s have been corrected with the International GNSS Service (IGS) precise products.t G r denotes the receiver clock offset absorbing the liner combination of receiver hardware biases b r,IF12 .N s r refers to the vector of the float ambiguities ( N s r,i ), which can recover its integer feature by deducting the hardware delay d r,i with the assistance of UPD products.Since the tropospheric hydrostatic delay can be corrected by the existing models, only the wet component Z r,w remains unknown.f aug P s r,i and f aug L s r,i represent the precise atmospheric corrections, which contain the tropospheric wet delay ( Z r,w ) and ionospheric delay with hardware delay at the receiver side ( I s r,1 ).δI s r,1 and δZ r,w denote the residual of the slant ionospheric delay on the first frequency and the tropospheric zenith wet delay, respectively.The numbers 1, 2 represent the two frequencies used in the satellite clock offset estimation process.However, if the frequency of the observation does not correspond to these frequencies, the inter-frequency bias b F must be considered.Moreover, the receiver inter-sys- tem bias S sys -G r is introduced for each system except for Global Positioning System (GPS) under the assumption that the multi-GNSS code observations share the same receiver clock.Thus, the parameters X to be estimated can be written as follows: where X u , t G r are estimated as a white noise model, and b F , S sys -G r , N s r are estimated as a random walk model.By correcting atmospheric products, the prior knowledge of δI s r,1 and δZ r,w are set to zero with high precision, which can be expressed as: Note that the ionospheric corrections absorb the hardware delays, which brings additional biases to the user side and needs to be carefully considered.The satelliteinduced bias remains consistent in both the server and user ends, thereby will not affect the positioning.However, the receiver-induced bias varies with stations.Once only single-frequency observations are used, this bias could be absorbed by the estimated receiver clock and ambiguity parameters.While for dual-frequency or mixed-frequency processing, an additional (5) With the assistance of precise atmospheric corrections, a rapid ambiguity resolution can be achieved.For the combined use of single-, dual-, and multi-frequency ambiguities, a cascade mixed-frequency ambiguity resolution strategy is conducted in Fig. 2. The float ambiguities derived from the above PPP-RTK model will be divided into three categories: the multi-frequency one, the dual-frequency one, and the single-frequency one.When there are more than two ambiguities, the Extra-Wide-Lane (EWL) and Wide-Lane (WL) ambiguities are formulated and fixed using the Least-Square AMBiguity Decorrelation Adjustment (LAMBDA) method (Temiissen, 1995) in sequence.When there are two ambiguities, only WL ambiguities are considered.The integer EWL ambiguities will be taken as virtual observations in the parameter estimation for accelerating the search process of WL ambiguities, and integer WL ambiguities will also serve as virtual observations for N1 ambiguity determination.Then N1 ambiguity resolution will be conducted also using the LAMBDA method.Note that multi-and dual-frequency ambiguities that are not fixed in WL mode will be removed in the next step.For the single-frequency ambiguities, they are fixed directly after correcting the UPD.To improve the success rate of ambiguity resolution, some quality control measures are employed to select the proper candidates for fixing, wherein the cutoff elevation, the minimum carrier-tonoise-density ratio C N 0 , and the maximum fractional part of ambiguity are set as the threshold.If the number of ambiguities after N1 ambiguity quality control is less than four or LAMBDA search fails, float solutions will be output.

Modified stochastic model considering the prior measurement accuracy
As the functional model describes the relationship between GNSS raw measurements and the parameters to be estimated, the observation stochastic model depicts the expectation and probability distribution of the random errors of measurements by defining a variance matrix, which is of great importance for parameter estimation and integer ambiguity resolution (Teunissen, 2018;Wang et al., 2002).
When low-cost GNSS devices are used for positioning in complex urban scenes, the traditional stochastic model using the elevation weighting scheme is no longer applicable as the satellites with high elevation may even  be blocked.In contrast, the carrier-to-noise-density ratio ( C N 0 ), which is closely associated with the propaga- tion path of the signal, is a better choice for defining the variance matrix.Moreover, the GNSS signal strength depends greatly on the type of receiver and antenna, which may vary significantly when equipped with different low-cost devices.In this way, to better account for the actual situation in defining the variance matrix, a modified stochastic model is given below: where C 0 and C 1 are the fitting parameters determined based on the prior observation information, and σ is the measurement precision determined by C N 0 with the unit in dB•Hz.
In order to obtain the appropriate fitting parameters C 0 and C 1 , the long-term observations of different types of GNSS devices are collected and tested in advance.The zero-baseline (or ultra-short baseline) method is applied to evaluate the noise of pseudorange observation in the condition that the satellite orbit error and clock offset error are eliminated, and the effect of atmospheric delay is alleviated (Amiri-simkooei & Tiberius, 2007;Pirazzi et al., 2017).The formula to evaluate the noise level of code observations is shown as follows: where r 1 is the reference station, and denotes the sin- gle-difference operator.The receiver clock offset is equal to the mean value of all satellites' code single-difference residual ( P s r 1 r 2 ,f ) at the current epoch.After subtracting the receiver clock offset from the code single-difference residual, the code noise e s r 1 r 2 ,f is acquired.After the computation of the measurement noise, the relationship between the noise and C N 0 is addressed.We divide the C N 0 into sections with an interval of 1 dB•Hz and calculate the mean value of the code noise at a specific section.As shown in Eq. ( 9), the least squares method is used to estimate the fitting parameters C 0 , C 1 .( 7) where σ refers to the residual of code observation noise, which equals the real code observation noise ( σ ) minus the first term of Taylor expansion ( σ 0 ).L 0 and L 1 repre- sent the coefficient of the parameters C 0 and C 1 respec- tively after linearization.C 0 0 and C 0 1 are the approximate values of the parameters to be estimated ( C 0 and C 1 ), respectively.
The fitting curves and results of helix antenna are taken as examples and shown in Fig. 3 and Table 1.Galileo and BDS of Fig. 3 represent Galileo navigation satellite system and BeiDou Navigation Satellite System, respectively.In this study, the weight of code observations was determined first, followed by the weight of phase observations, which was determined based on the empirical ratio of 100:1 of phase measurement accuracy to code measurement accuracy.

Data collection and experiment design
In this section, several experiments were designed for the verification and performance evaluation of the proposed PPP-RTK model using low-cost GNSS devices.The detailed information on GNSS devices including GNSS boards and antennas is shown in Table 2.A geodetic GNSS board (Septentrio PolarRx5S) and a geodetic-grade antenna (Trimble Zephyr 2), named SEPT + Trimble, was employed as a benchmark.Three types of low-cost GNSS devices: u-blox ZED-F9P GNSS board together with Vertical Dipole Antenna (VDA), Microstrip Patch Antenna (MPA), and Helix Antenna (HA) named as F9P + VDA, F9P + MPA, and F9P + HA, respectively were tested and compared.In kinematic experiments, a terminal that consists of a low-cost GNSS board (u-blox ZED-F9P) together with and more affordable in the mass market.
To evaluate the positioning performance, experiments were conducted in both static and kinematic modes.The static experiment was carried out from GPS Time (GPST) 12:00 to 22:00 on Day of Year (DOY) 334 (November 30th), 2021, with the experimental devices placed on the roof of a building, as shown in subgraph a of Fig. 4. The surrounding environment, including trees and tall buildings, as shown in subgraph b of Fig. 4, could potentially disrupt GNSS signals.
In addition to the static experiment, two kinematic tests were conducted to validate the low-cost PPP-RTK method.Experiment A took place on the playground from GPST 8:35 to 9:00 on DOY 017 (January 17th), 2022.Experiment B is a vehicle experiment conducted in a semi-urban environment from GPST 6:50 to 7:20 on DOY 12 (January 12th), 2022.The subgraphs a, b and c in Fig. 5 depict the experimental equipment and the trajectory of the two tests, respectively.For the kinematic experiments, the tightly coupled multi-GNSS RTK/INS solutions calculated by commercial Inertial Explorer (IE) 8.9 software (NovAtel, 2018) using the raw observation of the tactical-level IMU and GNSS terminals were taken as the reference coordinates.
Table 3 provides a detailed processing strategy for the PPP-RTK scheme.The multi-GNSS observations (GPS, Galileo, and BDS) with a sampling interval of 1 s were employed.In addition, the precise products from the German Research Centre for Geosciences (GFZ) were used to correct the orbit and clock offset errors, while the satellite phase biases were computed from a set of globally distributed stations using open-source software called GREAT-UPD (Li et al., 2021).Figure 6 shows the reference station network, which generates the atmospheric corrections for the clients.

Results and discussion
In this section, we first evaluate the observation quality of different low-cost devices in terms of data availability, C/N 0 , noise level, and cycle slip.Then the performances of PPP-RTK with different low-cost devices in static and kinematic conditions are investigated.

Observation characteristics of different low-cost GNSS devices
To figure out the observation characteristics of various low-cost GNSS devices, the data availability, C/N 0 , observation noise, and cycle slip rate with different lowcost devices are analyzed successively.
Firstly, the availability of observations is analyzed to explore the sensitivity and stability of the different lowcost GNSS devices.Figure 7 shows the availability of the code and phase observations on the second band of GNSS (e.g., GPS L2, BDS B2I, Galileo E5a), represented by the red and blue lines, respectively.Noting that the Septentrio GNSS receiver can track BDS Geosynchronous Earth Orbit (GEO) satellites from C01 to C05, while they are unavailable for the u-blox F9P due to the differences in inner acquisition algorithms and thresholds.For low-cost antennas, the observation continuity and integrity of some satellites, such as G26, E04 and E24, are inferior compared to the Trimble antenna.Overall, the analysis suggests that the performance in tracking satellites is largely determined by the GNSS boards used, while the continuity and integrity Table 4 reveals the percentage of single-frequency data in different scenarios.Notably, in the case of the campus playground, which is surrounded by dense trees and buildings, the single-frequency data accounts for the largest proportion.Conversely, in the vehicle experiment, where the car traversed a road section with a wide field of vision for nearly 20 min, the proportion of single-frequency data is comparatively smaller.Among the tested GNSS terminals, SEPT + Trimble exhibits the smallest proportion of single-frequency data, followed by F9P + Trimble, while the remaining low-cost terminals show similar proportions of singlefrequency data.This indicates that low-cost boards and the environment significantly influence data tracking capabilities.(Amami, 2022;Huang, 2022).It is worth noting that the C/N 0 of HA is higher than that of the other two low-cost antennas and is even comparable to that of Trimble when the elevation is below 45°.It suggests that the HA's high C/N 0 value may be attributed to its stronger tracking characteristics for low-elevation satellites.

Table 2 The information of GNSS devices used in this study
Note that the letters "D" and "H" represent the diameter and height of the antenna, respectively.The "±" sign indicates a range of values.For example, ( 40 Figure 9 demonstrates the code and phase noises of four GNSS terminals on the first band of GPS (L1), where the code noise is evaluated by the multipath combination (Kleusberg & Teunissen, 1996), and phase noise is calculated by the cubic difference between epochs (Amiri-Simkooei & Tiberius, 2007).The SEPT + Trimble device demonstrates superior performance, exhibiting smaller code and phase noise Root Mean Square (RMS) values of 0.302 and 0.005 m, respectively, compared to the three low-cost GNSS devices.For the low-cost GNSS devices, the RMS values for code noise of F9P + VDA, F9P + MPA, and F9P + HA are 0.3, 1.6, and 1.0 m, respectively, while RMSs for the phase noise of all three antennas are comparable, about 0.007 m.Note that the systematic and periodic deviations are observed in the code observations of F9P + MPA and F9P + HA, which lead to higher RMS values.One can conclude that under static conditions microstrip patch antenna is most vulnerable to noise and multipath interference, while vertical dipole antenna is of the best multipath suppression ability among the three low-cost GNSS antennas.
Figure 10 displays the code and phase observation noise series of different GNSS devices in real-world kinematic scenarios.In such situations, GNSS signals are particularly susceptible to interference or interruption, leading to multiple outliers in GNSS observations.This problem will be more pronounced when low-cost antennas are employed (Kaplan & Hegarty, 2017).As a result, the RMS of code noise for SEPT + Trimble and F9P + Trimble is 1.6 m, while it rises to 2.7 m, 2.8 m, and 1.8 m for F9P + VDA, F9P + MPA, and F9P + HA, respectively.In general, the noise level of phase measurement for the four terminals is comparable.
Moreover, the cycle slip rate, which is calculated as the percentage of the epochs with cycle slip over the total epochs, is also a key indicator for assessing carrier phase measurement.Table 5 shows the number of cycle slips and cycle slip rates of different GNSS devices under the static and the kinematic scenarios.As expected, cycle slips are more likely to occur on low-cost GNSS boards under kinematic conditions.For instance, the u-blox F9P receiver with a Trimble or a helix antenna exhibits a cycle slip rate of around 25‰ in the kinematic scenario, which is 3.6 times larger than that of SEPT + Trimble.Notably, among the three low-cost antennas, the helix antenna demonstrates the best performance.Additionally, Table 6 demonstrates the performance ranking of different low-cost devices in terms of C/N 0 , noise level, and cycle slip rate.The results show that the helix antenna exhibits the best performance in terms of C/N 0 and cycle slip rate.As for the noise level, the ver- tical dipole antenna and the microstrip patch antenna perform slightly better under static and kinematic conditions, respectively.

PPP-RTK performance in static experiments
Figure 11 illustrates the positioning error series and satellite numbers (NSAT) using a low-cost u-blox receiver with a microstrip patch antenna.It is evident that our PPP-RTK method achieves rapid and consistent centimeter-level positioning accuracy for most of the 7-h period.In contrast, traditional dual-frequency PPP-RTK is prone to the occlusions from the surrounding environment, leading to several reductions in the number of available satellites and reconvergence processes as depicted in Fig. 4. The performance degradation is alleviated using our proposed mixed-frequency PPP-RTK system, which demonstrates superior accuracy and stability.Table 7 further presents the fixed rate and positioning accuracy.The 3D-RMS with our new method is reduced from 0.53 m of the traditional method to 0.07 m, an improvement by 86.8%.In the traditional dual-frequency model, if only the single-frequency signal of a satellite is received, the satellite will be abandoned.Limited by the hardware performance of low-cost equipment, there are few observations available in some epochs with the traditional dual-frequency model.By contrast, the proposed mixedfrequency method can fully use all available measurements, and thus outperforms the traditional approach.
Compared with a survey-grade device, the positioning series with the new PPP-RTK method for three types of low-cost GNSS terminals (F9P + VDA, F9P + MPA, F9P + HA) is further presented in Fig. 12.The number of available satellites (NSAT), position dilution of precision (PDOP), and fixing status of each epoch is also given in this figure.It is worth noting that the low-cost GNSS devices can achieve the similar positioning accuracy as the geodetic one when the number of satellites is more than ten.However, the positioning performance significantly degrades when the number of the satellites decreases to about five, as illustrated in the second red box of Fig. 12     experience much more cycle slips than the geodetic one, and the F9P receiver with an MPA antenna is the most susceptible to cycle slips.This suggests that frequent cycle slips are the primary factor that affects the positioning performance of low-cost devices.
Table 9 shows the fixed rate and the positioning accuracy in the east, north, and up components.Among the four types of GNSS terminals, SEPT + Trimble and F9P + VDA perform best with a horizontal accuracy of 0.04 m and vertical accuracy of 0.07 m.F9P + HA with the positioning accuracy of (0.05 m, 0.03 m, 0.14 m) is slightly inferior to them, while F9P + MPA with the positioning accuracy of (0.16 m, 0.14 m, 0.50 m) exhibits the worst performance in the fixed rate.This may be attributed to the frequent cycle slip and large noise level of F9P + MPA.

PPP-RTK performance in kinematic experiments
Figure 13 illustrates the positioning error series and data availability for selected satellites in a kinematic scenario with a low-cost u-blox receiver and a low-cost helix antenna.The graph clearly demonstrates the effectiveness of our PPP-RTK method, as it consistently achieves rapid centimeter-level positioning accuracy for most of the 15-min period.In contrast, the traditional dual-frequency PPP-RTK method can only use the satellites with dual-frequency data, leading to several reductions in the    10 further presents the fixed rate and positioning accuracy.
The 3D-RMS of our new method is reduced from 2.94 m of the traditional method to 0.46 m, an improvement by 84.3%. Figure 14 depicts the positioning error series of all five GNSS terminals in kinematic Experiment A, and Table 11 reports their fixed rates and position accuracy.The SEPT + Trimble and F9P + Trimble perform best among all the terminals with a fixed rate of more than 99%.It illustrates that with the geodetic Trimble antenna the low-cost u-blox GNSS board can achieve the same positioning performance as the geodetic GNSS board Septentrio.Among the three low-cost GNSS devices, F9P + HA is the top performer with a positioning accuracy of (0.07 m, 0.07 m, 0.34 m), owing to its lower cycle slip rate and higher gain for low-elevation satellites.In contrast, F9P + VDA, which performs best under static conditions, exhibits the worst performance in kinematic scenarios.We find that F9P + VDA exhibits a high cycle slip rate of 35.8‰ in the initial period, which impacts its positioning performance.Our results suggest that the combination of the    F9P receiver and helix antenna, which can achieve geodetic-level positioning performance, is the most suitable for kinematic applications.Furthermore, Fig. 15 presents the position error series and satellite numbers of different devices in Experiment B. The number of available satellites (NSAT) of F9P + HA is also presented in the figure.It is worth noting that the number of available satellites for other devices is comparable to that of F9P + HA.During the first 18 min, the number of available GNSS satellites varies between 10 and 20, which ensures a continuous and reliable positioning.The F9P + HA configuration can achieve comparable positioning accuracy to that of SEPT + Trimble and F9P + Trimble.However, the positioning results obtained with F9P + MPA and F9P + VDA configurations are comparatively inferior.Then, the vehicle runs into a semiurban environment with several trees, tall buildings, and overpasses where the signal tracking becomes discontinuous and the NSAT drops frequently.Because of the frequent signal interruptions, the positioning accuracy of different devices, especially low-cost equipment, decreases.Notably, SEPT + Trimble and F9P + Trimble consistently achieve decimeter-level positioning accuracy, even in unfavorable observation conditions.However, the three low-cost devices exhibit poor performance.This phenomenon highlights the susceptibility of low-cost antennas to environmental interference.
Table 12 summarizes the fixed rate and positioning accuracy achieved by the different GNSS terminals during the first 18 min in a GNSS favorable scenario.SEPT + Trimble and F9P + Trimble exhibit the best performance, achieving centimeter-level accuracy with a fixed rate exceeding 99%.F9P + HA, which has the best performance in Experiment A, also outperforms the others in Experiment B with a horizontal accuracy of 0.2 m and vertical accuracy of 0.4 m.In contrast, F9P + VDA and F9P + MPA demonstrate poor positioning accuracy, despite their fixed rates of above 90%, due to the high incidence of wrongly fixed solutions.Hence, the helix antenna is the best choice among the three types of antennas under favorable kinematic scenarios.

Conclusions
This contribution proposed a modified PPP-RTK algorithm for the low-cost GNSS devices and evaluated its performance in different scenarios using different low-cost GNSS devices equipped with vertical dipole,  13 Positioning errors series and data availability for selected satellites using a low-cost u-blox receiver with a helix antenna.In the positioning errors series on the left panel, the red and blue points represent the proposed method and the traditional one used, respectively.In the data availability on the right panel, the red and blue points represent the data on the first band and the second band, respectively microstrip patch, and helix antennas, respectively.The observation quality of different low-cost devices in terms of data availability, C/N 0 , noise level, and cycle slip as well as the PPP-RTK positioning performance were investigated.
The experiment results illustrate that the ability to track satellites is determined by the GNSS board, while the continuity and integrity of the observations highly depend on the antenna.The helix antenna exhibits superior performance in C/N 0 compared to the other two low-cost antennas, approaching the level of a survey-grade Trimble antenna when the satellite elevation is below 45°.The microstrip patch antenna has the highest cycle slip rate under static conditions, while the vertical dipole antenna has the highest cycle slip rate under kinematic conditions.Furthermore, the level of measurement noise is also significantly affected by the type of GNSS antenna.In static mode, the code and phase noises of the three low-cost antennas range from 0.3-1.6 to 0.007 m, respectively, whereas in low-speed kinematic scenarios, they range from 1.8-2.8 to 0.011-0.034m, respectively.In vehicle experiments, the noise level increases to 3.9 m and 0.3 m for code and phase noise, respectively.
Compared to the traditional PPP-RTK model, the proposed method shows a better performance with the positioning error decreasing from 0.53 to 0.07 m.Both the static and low-speed kinematic experiments demonstrate that centimeter-level positioning accuracy can be achieved with low-cost devices.The experiment results indicate that the vertical dipole antenna with the positioning accuracy of (0.02 m, 0.01 m, 0.07 m) performs best under static conditions, and the helix antenna exhibits the best performance with the positioning accuracy of (0.07 m, 0.07 m, 0.34 m) under the low-speed kinematic scenario among the three types of low-cost antennas.Additionally, the helix antenna is the best choice for vehicle navigation with a fixed rate of over 95%, a horizontal accuracy of 0.2 m, and a vertical accuracy of 0.4 m under favorable observation conditions.

Fig. 3
Fig. 3 Fitting curves of helix antenna.The blue triangle and the orange line denote the raw mean value of noise and the fitting curve, respectively

Figure 8
Figure 8 illustrates the C/N 0 values on the first band of GPS (L1) for four types of GNSS terminals: SEPT + Trimble, F9P + VDA, F9P + MPA, and F9P + HA, represented by red, blue, green, and orange, respectively.A box diagram is drawn to demonstrate the distribution of the C/N 0 for every 15° interval.The C/N 0 values are rising with the increase of the elevation angles.When the elevation is larger than 45°, the C/N 0 values of the four types of GNSS terminals exceed 40 dB•Hz.Owing to the high gain of the Trimble antenna, the SEPT + Trimble shows the best performance in signal intensity(Amami, 2022;Huang, 2022).It is worth noting that the C/N 0 of HA is higher than that of the other two low-cost antennas and is even comparable to that of Trimble when the elevation is below 45°.It suggests that the HA's high C/N 0 value may be attributed to its stronger tracking characteristics for low-elevation satellites.

Fig. 4 Fig. 5
Fig. 4 Experimental devices and surrounding environment in the static experiment.a depicts the experimental devices; b represents the surrounding environment in the static experiment

Fig. 6
Fig.6The reference station network.The red triangles and blue dots represent the reference and user stations, respectively

Fig. 11
Fig. 11Positioning error series and satellite numbers using a low-cost u-blox receiver with a microstrip patch antenna.The red and blue points represent the proposed method and the traditional one, respectively

Fig. 12
Fig. 12 Positioning errors of different GNSS terminals in the static experiment Fig. 13Positioning errors series and data availability for selected satellites using a low-cost u-blox receiver with a helix antenna.In the positioning errors series on the left panel, the red and blue points represent the proposed method and the traditional one used, respectively.In the data availability on the right panel, the red and blue points represent the data on the first band and the second band, respectively

Fig. 14
Fig. 14 Positioning errors of different GNSS terminals in experiment A

Fig. 15
Fig. 15 Positioning errors of different GNSS terminals in experiment A

Table 1
The fitting result of helix antenna . Moreover,

Table 3
Processing strategies for PPP-RTK Observation availability of the different GNSS devices

Table 4
Percentage of single-frequency data in different scenarios

Table 5
Cycle slip rate of different GNSS terminals

Table 6
Data quality ranking of different low-cost GNSS antennas

Table 7
Fixed rate and positioning accuracy of the new method and the traditional one

Table 8
Cycle slip rate of different GNSS terminals

Table 9
Fixed rate and positioning accuracy of different GNSS terminals in static Experiment

Table 10
Fixed rate and positioning accuracy of the new method and the traditional one

Table 11
Fixed rate and positioning accuracy of different GNSS terminals in Experiment A

Table 12
Fixed rate and positioning accuracy of different GNSS terminals in favorable scenario of Experiment B