Multi-constellation single-frequency ionospheric-free precise point positioning with low-cost receivers

We analyze the observation quality, assess the performance and identify the constraints of quadruple-constellation single-frequency ionospheric-free precise point positioning (SF-IF PPP) with low-cost receivers. It is revealed that low-cost receivers with patch antennas exhibit lower C/N0 records and a weaker elevation dependence of C/N0 than the high-grade equipment. The results demonstrate that low-cost receivers can offer code measurements with similar noise compared to high-grade receivers providing that the multipath effect is eliminated. Regarding positioning performance, it is shown how SF-IF PPP for the high-grade receiver converges approximately two times faster than for the low-cost receiver with a patch antenna. It is confirmed that an application of a survey-grade antenna instead of the patch one noticeably enhances the performance of low-cost receiver SF-IF PPP. The study also reveals that the multipath effect is a dominant factor that constrains the performance of SF-IF PPP with low-cost receivers.


Introduction and motivation
An increasing interest in the precise positioning with lowcost receivers is recognized. In one of the first research in this area, Wisniewski et al. (2013) assessed the feasibility of precise positioning with the LEA-6T low-cost receiver. This was followed by the studies on single-baseline realtime kinematic (RTK) by Teunissen (2016, 2017), who investigated positioning with dual-system single-frequency (SF) and dual-frequency (DF) observations, respectively, and Garrido-Carretero et al. (2019), who showed the performance of SF-RTK with u-blox NEO-M8P. A centimeter-level positioning of multi-constellation RTK with Skytraq S2525F8-BD, u-blox ZED F9P and EVK-M8T receivers during high ionospheric activity and for long-baselines was reported in Teunissen 2019, 2020), respectively.
However, there is still a need to evaluate the accuracy level that may be reached with low-cost receivers in singlefrequency precise point positioning (SF-PPP) mode. This is warranted by a high demand for autonomous positioning and the fact that many devices are equipped with GNSS chipsets that offer observations on a single-frequency band only.
The most important challenge for SF-PPP is the elimination of the ionospheric delay. This error may be handled by using ionospheric models and global ionospheric maps. Respective research on this topic for low-cost receiver SF-PPP is found in Gill et al. (2017), use of real-time ionospheric products provided, e.g., by the Centre National d'Etudes Spatiales is found in Wang et al. (2020), and models such as ionospheric-weighted, ionospheric-constrained, and ionospheric-free (IF) were used in . The recently developed Smart-PPP for smartphone observations was addressed in Wang et al. (2021). Single-frequency ionospheric-corrected code observations of low-cost receivers were used in PPP by Nie et al. (2019). However, the method proposed by the authors combines SF measurements with conventional IF-DF linear combination (LC) and thus is suitable for DF receivers only.
Another option for handling ionospheric delay in SF positioning is to employ a combination of phase and code observations that was originally known as GRoup And PHase Ionospheric Correction (GRAPHIC) (Yunck 1996). The quality of code observations is a dominant factor that drives the positioning performance based on this linear combination. Therefore, a growing precision of code pseudoranges offered by modernized and emerging constellations is one of the aspects that induced us to draw special attention to this single-frequency ionospheric-free PPP model called SF-IF PPP. SF-IF PPP model was, until now, applied only to the observations of selected constellations, e.g., GPS + GLONASS as shown in (Cai et al. 2013), GPS + Galileo that was analyzed in (Bahadur and Nohutcu 2020) or single constellation, e.g., BDS (Wanninger and Beer 2015). Moreover, existing studies were based on the observations collected by high-grade receivers and antennas and with such equipment a decimeter-to sub-decimeter-level accuracy depending on the session length was reported (Bahadur and Nohutcu 2021).
Such progress motivates to assess, for the first time, the performance of quadruple-constellation SF-IF PPP of the low-cost receivers against high-grade ones and to identify the constraints of SF-PPP with low-cost receivers. An insight into the quality of GNSS observations of current low-cost receivers is provided. Furthermore, the observation quality and positioning results are analyzed in the light of the multipath impact, and the potential benefit from applying a survey-grade antenna instead of a low-cost one is verified.
In the next section, the functional model of SF-IF PPP is presented. The following section describes the data collection and the experiment design. Next, the observation quality analysis is given. Afterward, the performance of low-cost receiver SF-IF PPP against that of a high-grade receiver is assessed, and the impact of the multipath effect on the positioning is investigated. Finally, the conclusions are provided in the last section.

Methodology of single-frequency ionospheric-free precise point positioning
The observation equations for phase Φ and code P undifferenced signals between the satellite s and receiver r, and generalized to frequency i read: in which s r,0 denotes the geometric range between the satellite and a priori position of the station in meters; c represents the speed of light in meters per second; dt r and dt s correspond to the receiver and satellite clock corrections in seconds, respectively; i is the signal wavelength in meters on selected frequency i ; I s r,i and T s r are the slant ionospheric and tropospheric delays in meters, respectively;N s r,i is the integer ambiguity of the phase observable in cycles; B r,i and B s i are the receiver and satellite phase delays in meters, while b r,i and b s i correspond to the receiver and satellite (1) code pseudorange delays in meters, respectively; finally, s r,P and s r,Φ denote the noise of code and phase observations, correspondingly.
The availability of SF phase and code observations creates an opportunity to employ GRAPHIC IF-LC. The combination takes advantage of the fact that phase and code pseudorange observations transmitted on the same frequency band are subject to an ionospheric delay of the same magnitude but the opposite sign. However, it may also cause a model rank deficiency with single-epoch observations when using a least-squares adjustment (Sterle et al. 2015). The GRAPHIC combination applied to undifferenced observations is expressed as follows: where mh s r and mw s r are coefficients of the hydrostatic and non-hydrostatic tropospheric mapping functions, respectively; ZHD r and ZWD r represent zenith hydrostatic and non-hydrostatic tropospheric delays; e s r denotes the line-ofsight vector from the user receiver to the satellite; X r refers to the vector of the corrections to a priori coordinates of the station in meters;dt r is the receiver clock offset that combines receiver clock corrections and the temporally variable receiver hardware delays in seconds, A s r is the non-integer ambiguity term that aggregates one half of the phase integer ambiguity and satellite-specific but constant biases. Finally, s r,SF−IF denotes observational noise of SF-IF LC.
The final system of observation equations for combined quadruple-constellation signals that handles different for each system time scales reads: where GNSS constellations are denoted with the following superscripts: GPS (G), GLONASS (R), Galileo (E) or BDS (C) and dt R−G r , dt E−G r , dt C−G r state for GPS-GLONASS, GPS-Galileo, and GPS-BDS receiver inter system clock biases, respectively.
The vector of unknown estimates contains corrections to the a priori coordinates of the station x r , y r , z r , , and tropospheric nonhydrostatic zenith delay ZWD r .

Data collection
The dataset was acquired on November 17-19, 2020 andJune, 16-20, 2021 in the vicinity of UWM Olsztyn campus. GNSS measurements were collected with 5 s sampling interval and two collocated zero-baselines (ZB) (Fig. 1). One used two u-blox C099-F9P evaluation kits equipped with ZED F9P chipsets. An antenna splitter connected these to the same u-blox ANN-MB patch antenna having a circular ground plane. Another u-blox receiver of the same type and Trimble Alloy shared Leica LEIAR10 antenna. The antennas were mounted on the oblique roof with some obstacles. Therefore, a non-negligible multipath effect is expected. Since the inter-antenna distance was 50 cm, the observational conditions may be considered as the same for all receivers. The employed low-cost receivers acquire GPS L1 C/A and L2C, Galileo E1-B/C and E5b, BeiDou B1I and B2I, as well as GLONASS L1 and L2 signals. Since the SF-IF model only uses signals from the first frequency, the observations transmitted on other frequency bands have not been analyzed. During the first experiment, the Trimble receiver operated with enabled Everest™ multipath rejection technology that takes advantage of advanced digital signal processing within the satellite signal tracking loops. During the experiment conducted in June 2021, this option was disabled after two days of data collection to evaluate its potential advantages.

Signal strength
Since carrier-to-noise density ratio (C/N0) is prone to, inter alia, the losses and gains of signal strength driven by receiver hardware, this parameter is used as the first indicator that characterizes the quality of GNSS signals of lowcost receivers against high-grade ones. Figure 2 presents the frequency distribution histograms of C/N0s for the signals on the first frequency band for the dataset collected in November 2020. The histograms show a comparable C/N0 distribution for all receiver + antenna sets. In all cases, C/N0s exhibit the highest concentration in 41-45 dB-Hz bin for GPS, Galileo, and GLONASS and in 36-40 dB-Hz bin for BDS. For all constellations except BDS, the high-grade receiver achieved the highest frequency in 41-45 dB-Hz. A significantly higher number of Galileo observations with the strongest gain that fit 46-50 dB-Hz bin for the high-grade receiver as compared to the low-cost ones is also reported.
A more detailed view on the C/N0 is provided by Table 1 where differences between the C/N0 of the corresponding signals logged by Trimble Alloy + LEIAR10 and two u-blox receivers are presented. As expected, the low-cost receiver with a patch antenna always exhibits lower signal strength than the high-grade receiver. This is, for example, mirrored in between Trimble Alloy + LEIAR10 and u-blox ZED F9P + ANN-MBaverage C/N0 differences reaching  3.7 dB-Hz and 3.1 dB-Hz for Galileo and GLONASS observations, respectively. In the case of GPS and BDS constellations, these differences are not so significant and reach up to 0.8 dB-Hz. What is interesting, however, is that for GPS, BDS and GLONASS constellations, u-blox ZED F9P + LEIAR10 shows slightly higher values of average C/ N0 as compared to Trimble Alloy + LEIAR10. However, it is recalled that this low-cost receiver shares the survey-grade antenna with Trimble Alloy, and the latter receiver powers the shared antenna by a splitter.
What transpires from Fig. 3 is that C/N0 of the signals collected by Trimble Alloy + LEIAR10 demonstrate a strong elevation dependence that is typical for high-grade GNSS equipment but unusual for low-cost chipsets and antennas embedded in smartphones (Li and Geng 2019;Paziewski et al. 2021Paziewski et al. , 2019. The C/N0s for ZED F9P + ANN-MB exhibit the lowest elevation dependence when compared to both receivers equipped with the LEIAR10 antenna. The figure also shows how the patch antenna provides the observations of low-elevated satellites with C/N0s of a much noisier nature as compared to both receivers with a survey-grade antenna. Interestingly, also a low-cost receiver with a surveygrade antenna may be affected by occasional drops of the C/N0 for low-elevated satellites. A clear decline in the C/ N0 at high elevations for the low-cost receiver with a patch antenna is also worthy of notice.

Observation noise
To assess the impact of the noise and multipath effect on code observations, code-minus-carrier phase (CMC) linear combination is used, which for the satellite s and receiver r, is expressed as follows: where M s r,Φ and M s r,P denote the multipath effect in the unit of meters for phase and code observations, respectively.
CMC LC highlights the combined impact of noise and multipath of code observations but is also subject to ambiguities and satellite/receiver code and phase biases that are considered constant over time, as well as a doubled ionospheric delay that, in turn, is obviously a time-variant factor. All these unwanted effects were properly handled to isolate code noise and multipath as shown in (de Bakker et al. 2012). Figure 4 shows the example time series of CMC LC and the corresponding C/N0 records for GPS PRN 12 that are typical for all satellites. An observer may notice that surprisingly, GPS CMC LCs for low-cost receivers demonstrate less noise manifested by a random component of the time series compared to the high-grade receiver. This finding holds true only for the signals with relatively high signal strength, i.e., C/N0 over approximately 35 dB-Hz. On the contrary, CMC values that are characterized with low C/N0 are much more disturbed for low-cost receivers when compared to the highgrade one and exhibit a high magnitude of the low-frequency component. Such effects are clearly seen at the beginning and at the end of the observational arcs when low-elevated signals are acquired. Such behavior of low-cost CMC time series may be attributed to the multipath effect. A less noisy nature of the random component of CMC LC time series for low-cost receivers compared to Trimble Alloy may result  (Park et al. 2017;Odolinski and Teunissen 2020). However, at this time, we are not informed about specific signal processing algorithms that are implemented in the employed low-cost receivers. In Fig. 5, standard deviations (STD) of CMC LC as a function of signal strength are shown. The figure confirms a much noisier nature of code observations with C/N0 lower than approximately 35 dB-Hz for low-cost receivers compared to the high-grade one. This finding is especially true for the receiver equipped with a path antenna. The figure also reveals that code observations of the high-grade receiver are much less dependent on C/N0 than those of low-cost equipment. Table 2 reports CMC LC STDs that are distinguished between these related to all of the acquired observations and those related to the dataset with excluded observations characterized with low C/N0 records (below 35 dB-Hz).
As expected, the CMC STDs for the low-cost receiver with a patch antenna always exhibit the highest values. What was also found is an identical performance of highgrade receivers and the low-cost receiver with a surveygrade antenna in terms of GPS L1 code observation noise expressed by CMC STD equaled to 32 cm. Table 2 also reveals that Trimble Alloy significantly outperforms lowcost receivers in terms of Galileo E1 code accuracy. For the former receiver CMC STD equals to 25 cm, whereas for low-cost receivers, these statistics reach 37 cm and 41 cm, respectively. Such results are in agreement with C/N0 frequency distribution histograms given in Fig. 2 where numerous Galileo observations with the highest C/N0 logged by the high-grade receiver were shown. The opposite situation occurs for code noise of BDS and GLONASS. Surprisingly, in these cases, the low-cost receiver with a survey-grade antenna outperforms high-grade equipment. However, this is consistent with a slightly higher average C/N0 for the former receiver than the latter, reported in Table 1. One should note that BDS-2 code observations are subject to elevation-dependent variations of code pseudoranges that may gain the CMC STDs (Wanninger and Beer 2015). Overall, the results are in line with the findings of Odolinski and Teunissen (2019), who reported a code noise fitting the range of 43-49 cm depending on the constellation for u-blox EVK-M8T + ANN-MS.
Further inspection of the statistics given in Table 2 reveals a noticeable reduction of CMC STDs for low-cost receivers after the elimination of low-gain signals. Such high enhancement is not observed for the high-grade receiver, which was however expected if we recall the results of CMC STD as a function of C/N0 illustrated in Fig. 5. Taking as an example the GPS L1 observations, after the elimination of low-gain signals, Trimble Alloy CMC STD improved only by a 1 cm, whereas for the low-cost receivers with a surveygrade antenna and a patch antenna, these improvements reached 4 cm and 8 cm, respectively. The outcomes confirm that low-cost GNSS equipment is highly prone to the multipath effect that mostly affects low-elevated and therefore low-gain signals, and may exhibit a competitive performance to that of high-grade grade equipment providing that multipath affected observations were eliminated or corrected.

Phase and code noise analysis with zero-baseline setup
Double-differenced (DD) observations of zero-baseline (ZB) setup built of two ZED F9P receivers connected by a splitter to the same ANN-MB patch antenna are used to assess the noise of phase observations. Phase STDs of low-cost receivers that were normalized to the undifferenced level are presented in Table 3. The results show a precision of phase observations comparable to highgrade receivers and an agreement with the findings of Odolinski and Teunissen (2019). For all except the GLONASS constellation, the phase STD did not exceed 2.9 mm. The noise of GLONASS phase observations is higher and equals to 4.1 mm. Nonetheless, the noise of phase observations has a marginal impact on SF-IF PPP performance since it is approximately two orders lower as compared to that of code observations.
More importantly, the ZB setup allows the separation of the observational noise from the multipath effect since the latter affects the observations of both receivers of the ZB setup in the same manner and thus is reduced in double-differencing. The results given in Table 3 confirm a high impact of multipath on low-cost equipment since STDs of code noise are noticeably lower as compared to those reported in Table 2, which correspond to the combined impact of noise and multipath. U-blox ZED F9P + ANN-MB code STDs reduce from 39 cm, 41 cm, 42 cm, and 48 cm for GPS, Galileo, BDS, and GLONASS to 29 cm, 32 cm, 32 cm, and 35 cm, respectively, when the multipath effect has been eliminated in ZB setup.

Performance assessment of SF-IF PPP with low-cost receivers
In this section, the details of the processing strategy are provided and the positioning performance assessment of quadruple-constellation SF-IF PPP is presented. A full dataset was processed using processing windows of different lengths.

Processing strategy
Correction models that are typically considered in the multiconstellation PPP were employed. Specifically final MGEX orbits and clock corrections were used (Prange et al. 2020). The MGEX bias products are applied to account for satellite code biases and to correct satellite clock offsets containing ionospheric-free satellite code hardware delays for singlefrequency measurements. A priori hydrostatic tropospheric delays are modeled using the Saastamoinen model and Global Mapping Function (Boehm et al. 2006;Saastamoinen 1973). Other corrections to eliminate effects such as relativity, solid earth tides, phase wind-up effect, and antenna phase center offset and variations with igs14.atx files are also accounted for according to the conventional models (Kouba and Héroux 2001).
The parameter estimation is performed using the Kalman filter. Spectral density values of process noise equal to 10 2 m 2 s −1 , 10 −7 m 2 s −1 , and 10 −9 m 2 s −1 for receiver clock offset, receiver inter system clock biases, and nonhydrostatic zenith tropospheric delays, respectively, were adopted in random walk processes. Static coordinates and  float ambiguity parameters are considered as constant values. The zenith-referenced STDs of code observations are adopted, and taking advantage of the observation assessment results, an elevation-dependent weighting scheme based on the cosecant function is used. The data were processed with 5 s sampling interval and a cutoff angle of 10 degrees.

Positioning results
The performance of quadruple-constellation SF-IF PPP was evaluated in terms of positioning errors and convergence time. Positioning errors are computed by referring the coordinate estimates to the benchmark coordinates. Convergence time is defined here as the period required to achieve and keep a position for at least 5 min with a three-dimensional (3D) error lower than the adopted thresholds of 1 m or 1.5 m. Figure 6 shows a snapshot positioning performance of low-cost receivers against a high-grade receiver. The figure presents example time series of quadruple-constellation SF-IF PPP positioning errors for static daily solutions on November 17-18, 2020. What transpires from the figure is a noticeably slower convergence for both low-cost receivers when compared to the high-grade one. A clear benefit from the application of survey-grade antenna instead of the patch one is also reported, as u-blox ZED F9P + LEIAR10 gives lower positioning errors and faster convergence when compared to ZED F9P + ANN-MB. Figure 7 presents SF-IF PPP 3D positioning root-meansquare (RMS) errors that were averaged over the full dataset. The Trimble receiver results are distinguished between those that are related to the period with an enabled Ever-est™ multipath rejection technology and those that refer to the period after disabling the function. The figure shows how the 3D positioning accuracy converges to the level of 0.5 m after approximately 4-5 h depending on the employed receiver + antenna set. The full dataset results confirm the snapshot of daily solutions given in Fig. 6, as Trimble Alloy converges noticeably faster than both low-cost receivers. A convergence enhancement from the application of the Everest™ algorithm in the Trimble receiver is worthy of attention.
A detailed evaluation of the positioning results begins with the convergence time assessment. Figure 8 provides a focus on the average over the full dataset 3D RMS positioning errors that were obtained when the filter was reinitialized every 1 h. The most evident outperformance of high-grade receivers over the low-cost ones is observed during the first 45 min of data processing. Suppose we focus on the first 15-20 min. In that case, we will also notice the benefit of applying the survey-grade antenna to low-cost receiver positioning and the advantage from Trimble's Everest™ technology. Nonetheless, after approximately a 1 h initialization period, all receivers reach a comparable accuracy level. Table 4 summarizes the performance of the convergence time. Statistics prove how the high-grade receiver converges faster than the low-cost ones. Trimble Alloy with enabled Everest™ requires 20.8 min to converge to the defined 3D accuracy threshold of 1 m. The corresponding indicators for ZED F9Ps with a patch or a survey-grade antenna are  The solutions obtained using a filter reinitialization every 6 h were averaged over the full dataset more than twice as high. An advantage from the Everest™ algorithm is also proved. The convergence time extends for almost one-quarter from 20.8 up to 25.7 min, when this option is disabled in the Trimble receiver.
The statistics for a less challenging 3D positioning error threshold of 1.5 m aim to highlight the divergences between the receivers during the filter initialization period. Trimble Alloy with Everest™ is the best performing receiver, whereas ZED F9P with a patch antenna is the worst one. The latter needs 37.8 min to reach the defined 3D positioning accuracy level, whereas it is only 8.3 min for the former.
To further investigate the inequalities in the positioning performance among the employed receivers, the empirical RMSs of positioning errors are presented in Table 5. The statistics correspond to the accuracy level that may be reached after 15 min, 45 min, and 6 h of data processing. For the first two intervals, the positioning accuracy is expected to be significantly lower than that of daily solutions. In the latter case, the positioning accuracy may converge to subdecimeter level, as shown in Fig. 6 and past studies (Cai et al. 2013). However, such long site occupations are rarely expected in real-life applications.
Overall, Table 5 shows that the positioning errors are at a dm-level depending on the adopted processing window and receiver. As we may read from the table, Trimble Alloy with Everest™ always outperforms both low-cost receiver sets and Trimble Alloy with a disabled Everest™ in terms of positioning accuracy. The Trimble Alloy performance corresponds to the previous results reported for high-grade receivers in Bahadur and Nohutcu (2021), which showed an 3D RMS of 0.5-0.7 m depending on the satellite products used, and also confirmed the better performance of multiconstellations solution over single-, dual-or triple-constellation solutions.
Moreover, the statistics manifest significant divergences in positioning performance between high-grade and lowcost receivers during the initialization period. Taking 15 min long processing window as an example, averaged RMSs of positioning errors for Trimble Alloy with enabled Everest™ are of 47.1 cm, 68.7 cm, and 49.0 cm for North, East, and Up components, respectively. As expected, a lower accuracy is reached by the Trimble Alloy with the disabled Ever-est™ and in this case RMSs are of 45.0 cm, 76.6 cm, and 64.0 cm, for the respective coordinate components. ZED F9P + LEIAR10 antenna gives comparable to the high-grade receiver RMSs for North and East components of 46.0 cm, 73.4 cm, respectively, but over three times higher RMS for the vertical one, that is of 190.8 cm. This between-receiver divergence reaches an even higher magnitude when the high-grade equipment is compared to the low-cost receiver  with a patch antenna. Such outcomes prove how the multipath effect, to which the low-cost receivers are especially prone, propagates into estimates and, as a consequence, produces biased position time series and deteriorates position convergence.
Inspecting Table 5 again, the reader may notice that after a few hours, the position converges to the 3D positioning RMS of about 3-4 dm for all receivers and therefore, the important disagreements between the receivers disappear.

Insight into the observation residuals
SF-IF observation residuals are examined to understand the poorer positioning performance of low-cost receivers when compared to the high-grade ones. An insight into the example residuals of GPS PRN 12 is given as they are typical for all satellites. Figure 9 shows the time series of residuals in gray color for June 16-17, 2021 and satellite elevation in a blue color when the Trimble receiver operated with enabled Everest™ algorithm. The residual time series demonstrate two components, a high-frequency random component that is most likely driven only by observation noise and a lowfrequency one related to unmodeled errors that remain in the SF-IF observations, which is predominately the multipath effect (Sterle et al. 2015). The red line denotes the residuals that were filtered with a moving average to expose the lowfrequency component. The STDs given in Fig. 9 are defined from the low-frequency component of the residuals.
The hypothesis that the multipath effect is mainly responsible for the poor performance of low-cost receiver positioning may be verified by an analysis of day-by-day repeatability of the residual time series. It is evident that the multipath effect and antenna phase-center variation error exhibit repeatability since the same measuring environment is provided each sidereal day by the repetition of GPS orbits (Xue et al. 2021). Indeed, a consistency between the low-frequency components of the residuals of consecutive days for the low-cost receivers that are given in Fig. 9 is discovered.
A high impact of the multipath is also confirmed by the change in the nature of the Trimble residuals after the deactivation of the Everest™ multipath rejection algorithm. Such residuals, shown in Fig. 10, are subject to evident low-frequency components characterized with STDs of 0.13-0.14 m that are about 60% higher as compared to the period when the Everest™ was enabled (c.f. Figure 9). A noticeable drop in STD of the filtered residuals driven by the multipath mitigation algorithm implies that there is still some room for positioning improvement, providing that the multipath is carefully handled.
It also follows from Fig. 9 that the estimated observation residuals for both low-cost receivers exhibit a lower noise that is represented by a random component of the time series and a higher impact of low-frequency effects that is likely driven by multipath as compared to high-grade receivers. Such results are consistent with code-minus-carrier LC given in Fig. 4.
If we further investigate Fig. 9, we discover that an application of the survey-grade antenna to the low-cost receiver did not significantly change the nature of residual time series as they still exhibit a noticeable low-frequency component. This finding is supported by the example STDs of SF-IF residuals given in the figure, as they are always higher than the Trimble receiver. Therefore, the low-frequency component of the residual time series is attributed to the way how Red lines represent the residuals filtered with a moving average filter and corresponded to low-frequency unmodeled effects. Satellite elevation is shown in blue. The first and second panels correspond to the residuals for Trimble Alloy + LEIAR10 with enabled Everest™. The U-blox ZED F9P + LEIAR10 residuals are given in panels three and four, respectively. ZED F9P + ANN-MB residuals are shown in fifth and sixth panels, correspondingly receivers handle the multipath. It may be concluded that the application of a survey-grade antenna instead of a patch one, enhances the PPP performance of low-cost receivers to some extent that is in agreement with the past studies on RTK (Odolinski and Teunissen 2016). However, it requires advanced multipath mitigation algorithms to reach the positioning performance close to that of high-grade receivers. Given that theoretically, phase multipath may reach only up to a quarter of a signal wavelength, it becomes clear that the multipath of code observations should be a subject of further investigations with regard to the low-cost receiver SF-PPP.

Conclusions
First, it was shown that the low-cost receiver with a patch antenna exhibits lower records of the C/N0 than the highgrade receiver. Then, a weaker elevation dependence of C/N0 for a patch antenna as compared to the high-grade one was proved. The results also showed how the low-cost antenna provides the observations of low-elevated satellites with C/N0s of much noisier nature than the receivers with a survey-grade antenna and exhibits a drop in C/N0 at high elevations.
The observation assessment revealed that low-cost receivers may offer code measurements characterized with a similar noise compared to high-grade receivers providing that the multipath effect has been eliminated. A slight outperformance of the low-cost receivers over the high-grade one in terms of code noise was noticed for GPS and GLO-NASS constellations when low-gain signals were eliminated. Moreover, the ZB setup experiment proved a high precision of low-cost receiver phase observations corresponding to high-grade receivers.
Assessing the positioning performance in a non-negligible multipath environment, it was then demonstrated how SF-IF PPP for the high-grade receiver converges approximately two times faster than for the low-cost receiver with a patch antenna. Analysis of the positioning errors confirmed a significant outperformance of the high-grade equipment when compared to the low-cost one during the initialization period. The tests demonstrated how after approximately 6 h the low-cost receivers may reach a comparable to highgrade ones accuracy level characterized with 3D RMS of 0.3-0.4 m. It was also shown how the application of a survey-grade antenna instead of a patch one might enhance the low-cost receiver positioning by improving the convergence time and reducing positioning errors.
Finally, SF-IF observation residuals were investigated to understand the poorer positioning performance of the lowcost receivers when compared to the high-grade ones. It was revealed that the multipath effect is a significant part of the observation residuals and a dominant factor that constrains the performance of SF-IF PPP with low-cost receivers. It is anticipated that an application of advanced multipath mitigation algorithms similar to Trimble's Everest™ may move the positioning performance of low-cost receivers towards the level close to that of high-grade receivers.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http:// creat iveco mmons. org/ licen ses/ by/4. 0/.