Quality assessment of GPS rapid static positioning with weighted ionospheric parameters in generalized least squares

Precise GPS positioning requires the processing of carrier-phase observations and fixing integer ambiguities. With increasing distance between receivers, ambiguity fixing becomes more difficult because ionospheric and tropospheric effects do not cancel sufficiently in double differencing. A popular procedure in static positioning is to increase the length of the observing session and/or to apply atmospheric (ionospheric) models and corrections. We investigate the methodology for GPS rapid static positioning that requires just a few minutes of dual-frequency GPS observations for medium-length baselines. Ionospheric corrections are not required, but the ionospheric delays are treated as pseudo-observations having a priori values and respective weights. The tropospheric delays are reduced by using well-established troposphere models, and satellite orbital and clock errors are eliminated by using IGS rapid products. Several numerical tests based on actual GPS data are presented. It is shown that the proposed methodology is suitable for rapid static positioning within 50–70 km from the closest reference network station and that centimeter-level precision in positioning is feasible when using just 1 min of dual-frequency GPS data.     

Time-Relative Positioning with a Single Civil GPS Receiver

Time-relative positioning is a recent method for processing GPS phase observations. The operational method undertaken in this paper consists of the following steps: first, recording phase observations at a station of known coordinates; second, moving the GPS receiver to an unknown station (which can be located up to a few hundred meters away, dependint on what type of transportation – e. g., walking, motorcycle – is available) while continuously observing carrier phases; and, third, recording phase observations at a second station of unknown coordinates with a single GPS receiver. To obtain the position of the unknown station relative to the first (known) station, the processing method uses combined observations taken at two different epochs and two different stations with the same receiver. For this reason, the errors that vary between two epochs must be taken into account in an appropriate way, especially errors in satellite clock corrections and ephemerides, and errors related to tropospheric and ionospheric delays. Ionospheric modeling using IONEX files (the ionospheric maps calculated by the International GPS Service) was also tested to correct L1 phase observations. This method has been used to calculate short vectors with an accuracy of a few centimeters (for a processing interval of 30 s) using a single civil GPS receiver. ? 2001 John Wiley & Sons, Inc.     

Assessment of ZTD derived from ECMWF/NCEP data with GPS ZTD over China

The accuracy and feasibility of computing the zenith tropospheric delays (ZTDs) from data of the European Center for Medium-Range Weather Forecasts (ECMWF) and the United States National Centers for Environmental Prediction (NCEP) are studied. The ZTDs are calculated from ECMWF/NCEP pressure-level data by integration and from the surface data with the Saastamoinen model method and then compared with the solutions measured from 28 global positioning system (GPS) stations of the Crustal Movement Observation Network of China (CMONOC) for 1 year. The results are as follows: (1) the error of the integration method is 1–3 cm less than that of the Saastamoinen model method. The agreement between the ECMWF ZTD and GPS ZTD is better than that between NCEP ZTD and GPS ZTD; (2) the bias and root mean square difference (RMSD), especially the latter, have a seasonal variation, and the RMSD decreases with increasing altitude while the variation with latitude is not obvious; and (3) when using the full horizontal resolution of 0.5° × 0.5° of the ECMWF meteorological data in place of a reduced 2.5° × 2.5° grid, the mean RMSD between GPS and ECMWF ZTD decreases by 4.5 mm. These results illuminated the accuracy and feasibility of computing the tropospheric delays and establishing the ZTD prediction model over China for navigation and positioning with ECMWF and NCEP data.     

电离层高阶项改正对参考框架实现及测站坐标的影响分析

根据电离层高阶项改正模型,确定了高阶项延迟对不同GPS观测量的影响量。根据IERS协议2010推荐的最新模型,对全球均匀分布的104个IGS基准站数据进行了重新处理,研究了电离层高阶项延迟(二、三阶项)对坐标参考框架实现及测站坐标的影响,量化了不同地磁模型下的高阶项改正影响变化。结果表明,电离层延迟高阶项改正对参考框架原点有较大影响,犣方向的平移可达20mm,犡、犢方向的平移所受的影响相对较小,大部分维持在5mm以下;电离层高阶项延迟会引起测站坐标变化的区域性偏移现象;地磁模型的不同会导致参考框架原点犣方向有着约10mm的变化,在局部区域也会引起显著的测站坐标变化差异。最后对电离层高阶项改正的影响结果进行了分析与讨论。     

Combination of GNSS and SLR observations using satellite co-locations

Satellite Laser Ranging (SLR) observations to Global Navigation Satellite System (GNSS) satellites may be used for several purposes. On one hand, the range measurement may be used as an independent validation for satellite orbits derived solely from GNSS microwave observations. On the other hand, both observation types may be analyzed together to generate a combined orbit. The latter procedure implies that one common set of orbit parameters is estimated from GNSS and SLR data. We performed such a combined processing of GNSS and SLR using the data of the year 2008. During this period, two GPS and four GLONASS satellites could be used as satellite co-locations. We focus on the general procedure for this type of combined processing and the impact on the terrestrial reference frame (including scale and geocenter), the GNSS satellite antenna offsets (SAO) and the SLR range biases. We show that the combination using only satellite co-locations as connection between GNSS and SLR is possible and allows the estimation of SLR station coordinates at the level of 1–2 cm. The SLR observations to GNSS satellites provide the scale allowing the estimation of GNSS SAO without relying on the scale of any a priori terrestrial reference frame. We show that the necessity to estimate SLR range biases does not prohibit the estimation of GNSS SAO. A good distribution of SLR observations allows a common estimation of the two parameter types. The estimated corrections for the GNSS SAO are 119 mm and −13 mm on average for the GPS and GLONASS satellites, respectively. The resulting SLR range biases suggest that it might be sufficient to estimate one parameter per station representing a range bias common to all GNSS satellites. The estimated biases are in the range of a few centimeters up to 5 cm. Scale differences of 0.9 ppb are seen between GNSS and SLR.     

Analysis of long-range network RTK during a severe ionospheric storm

The network-based GPS technique provides a broad spectrum of corrections to support RTK (real-time kinematic) surveying and geodetic applications. The most important among them are the ionospheric corrections generated in the reference network. The accuracy of these corrections depends upon the ionospheric conditions and may not always be sufficient to support ambiguity resolution (AR), and hence accurate GPS positioning. This paper presents the analyses of the network-derived ionospheric correction accuracy under extremely varying – quiet and stormy – geomagnetic and ionospheric conditions. In addition, the influence of the correction accuracy on the instantaneous (single-epoch) and on-the-fly (OTF) AR in long-range RTK GPS positioning is investigated, and the results, based on post-processed GPS data, are provided. The network used here to generate the ionospheric corrections consists of three permanent stations selected from the Ohio Continuously Operating Reference Stations (CORS) network. The average separation between the reference stations was ∼200 km and the test baseline was 121 km long. The results show that, during the severe ionospheric storm, the correction accuracy deteriorates to the point when the instantaneous AR is no longer possible, and the OTF AR requires much more time to fix the integers. The analyses presented here also outline the importance of the correct selection of the stochastic constraints in the rover solution applied to the network-derived ionospheric corrections.     

Atmospheric loading corrections at the observation level in VLBI analysis

This paper investigates whether in very long baseline interferometry (VLBI) analysis atmospheric loading corrections should be applied a priori at the observation level or whether it is sufficient to correct for atmospheric loading effects a posteriori by adding constant values per session to the estimated station coordinates. Simulated observations at single stations corresponding to the precise point positioning approach of global navigation satellite systems show that the atmospheric loading effect can be fully recovered by a posteriori corrections, i.e., the height differences between both approaches stay well below 1 mm. However, real global VLBI network solutions with sessions from 1984 to 2008 reveal that the effect of neglected atmospheric loading corrections at the stations is distributed to the other stations in the network, thus resulting in station height differences between solutions with observation level and with a posteriori corrections which can be as large as 10 mm and a ‘damping’ effect of the corrections. As soon as the terrestrial reference frame and the corresponding coordinate time series are determined, it would be conceptually wrong to apply atmospheric loading corrections at the VLBI stations. We recommend the rigorous application of atmospheric loading corrections at the observation level to all stations of a VLBI network because the seven parameters for translation, rotation, and in particular the network-scale of VLBI networks are significantly affected.     

GPS tomography in the polar cap: comparison with ionosondes and in situ spacecraft data

Tomographic 4D reconstructions of ionospheric anomalies appearing in the high-latitude polar cap region are compared with plasma density measurements by digital ionosonde located near the north magnetic pole at Eureka station and with in situ plasma measurements on-board DMSP spacecraft. The moderate magnetic storm of 14–17 October 2002 is taken as an example of a geomagnetic disturbance which generates large-scale ionospheric plasma anomalies at mid-latitudes and in the polar cap region. Comparison of the GPS tomographic reconstructions over Eureka station with the ionosonde measurements of the F layer peak densities indicates that the GPS tomography correctly predicts the time of arrival and passage of the ionospheric tongue of ionization over the magnetic pole area, although the tomographic technique appears to under-estimate the value of F peak plasma density. Comparison with the in situ plasma measurements by the DMSP SSIES instruments shows that the GPS tomography correctly reproduces the large-scale spatial structure of ionospheric anomalies over a wide range of latitudes from mid-latitudes to the high-latitude polar cap region, though the tomographic reconstructions tend to over-estimate the density of the topside ionosphere at 840 km DMSP orbit. This study is essential for understanding the quality and limitations of the tomographic reconstruction techniques, particularly in high-latitude regions where GPS TEC measurements and other ionospheric data sources are limited.     

Numerical simulation of troposphere-induced errors in GPS-derived geodetic time series over Japan

Troposphere-induced errors in GPS-derived geodetic time series, namely, height and zenith total delays (ZTDs), over Japan are quantitatively evaluated through the analyses of simulated GPS data using realistic cumulative tropospheric delays and observed GPS data. The numerical simulations show that the use of a priori zenith hydrostatic delays (ZHDs) derived from the European Centre for Medium-Range Weather Forecasts (ECMWF) numerical weather model data and gridded Vienna mapping function 1 (gridded VMF1) results in smaller spurious annual height errors and height repeatabilities (0.45 and 2.55 mm on average, respectively) as compared to those derived from the global pressure and temperature (GPT) model and global mapping function (GMF) (1.08 and 3.22 mm on average, respectively). On the other hand, the use of a priori ZHDs derived from the GPT and GMF would be sufficient for applications involving ZTDs, given the current discrepancies between GPS-derived ZTDs and those derived from numerical weather models. The numerical simulations reveal that the use of mapping functions constructed with fine-scale numerical weather models will potentially improve height repeatabilities as compared to the gridded VMF1 (2.09 mm against 2.55 mm on average). However, they do not presently outperform the gridded VMF1 with the observed GPS data (6.52 mm against 6.50 mm on average). Finally, the commonly observed colored components in GPS-derived height time series are not primarily the result of troposphere-induced errors, since they become white in numerical simulations with the proper choice of a priori ZHDs and mapping functions.     

Comparison of a GPS-defined global reference frame with ITRF2000

The Global Positioning System is a constellation of 24–28 satellites, which can be used to define a global terrestrial reference frame. Daily offsets between a GPS defined frame and ITRF2000 have been estimated using more than a decade of GPS observations from 1990–2001. A linear fit to the full span of data shows agreement between the two frames at the level of –1 ppb and –0.1 ppb/year for scale, 5 mm and 0 mm/year for the X component of center of mass, –2 mm and –3 mm/year for the Y component, and 4 mm and 6 mm/year for the Z component. GPS is a viable tool for defining the global reference frame either alone, or in combination with other geodetic techniques. Electronic Publication     

Precise ionosphere-free single-frequency GNSS positioning

Ionospheric delays can be efficiently eliminated from single-frequency data using a combination of carrier phases and code ranges. Unfortunately, GPS and GLONASS ranges are relatively noisy which can limit the use of the positioning method. Nevertheless, position standard deviations are in the range of 6–8 cm (horizontal) and 7–9 cm (3d) obtained from diurnal data batches from selected IGS reference stations can be further reduced to 2–3 cm (3d) for weekly smoothed averages. GPS data sets collected in Ghana (Africa) reveal a typical level of 10 cm of deviation that must be anticipated under average conditions. Looking at the future of GNSS, the European Galileo system will, in contrast to GPS, provide the broadband signal E5 that is by far less affected by multipath thus providing rather precise range measurements. Simulated processing runs featuring both high ionospheric and tropospheric delay variations show a 3d position precision of 4 cm even for a data batch as short as just 1 h, whereas GPS L1/Galileo E1 performance is close to 13 cm for the same data set.     

GPS precision as a function of session duration and reference frame using multi-point software

As an aid to survey design, we used data acquired from three European continuous GPS networks to test the precision of position estimates from static observations as a function of the length of the observing session and the number and distribution of reference stations. Our criterion was the weighted RMS of estimates over 31 days with respect to coordinates determined from 24-h sessions over a 2-year period. With a single reference station, a precision of 3 mm horizontal and 10 mm vertical could be achieved reliably only for session lengths of 3 h or longer and baselines less than 200 km. If four or more reference stations are used, these levels of precision were usually achieved with sessions as short as 2 h. With sessions 6 h or longer and four or more reference stations, the precision is typically 1–2 mm in horizontal and about 3–5 mm in vertical. Increasing the number of reference stations further provides only marginal improvement. Although there is some variation in precision in 4-station networks with the choice of reference stations, the dependence on distance and geometric distribution is weak.     

Absolute phase center corrections of satellite and receiver antennas

Results of the estimation of azimuth-dependent phase center variations (PCVs) of GPS satellite antennas using global GPS data are presented. Significant variations of up to ±3–4 mm that are demonstrated show excellent repeatability over eight years. The application of the azimuthal PCVs besides the nadir-dependent ones will lead to a further reduction in systematic antenna effects. In addition, the paper focuses on the benefit of a possible transition from relative to absolute PCVs. Apart from systematic changes in the global station coordinates, one can expect the GPS results to be less dependent on the elevation cut-off angle. This, together with the significant reduction of tropospheric zenith delay biases between GPS and VLBI, stands for an important step toward more consistency between different space geodetic techniques.     

Short note: Crustal deformation in the key stone network detected by satellite laser ranging

The paper presents the results of crustal deformation, as evidenced by changed station coordinates, in the Tokyo metropolitan area detected by the satellite laser ranging (SLR) technique. The coordinates of two Key Stone SLR stations, Tateyama and Kashima, were determined from 4 weeks of orbital arcs of the LAGEOS-1 and LAGEOS-2 satellites with respect to 16 SLR stations kept fixed in the ITRF2000 reference frame. The station coordinates were calculated using the NASA GEODYN-II orbital program. The orbital RMS-of-fit for both satellites was 16 mm. The standard deviation of the estimated positions was 3 mm. A jump of about 5 cm in the baseline length between the Kashima and Tateyama stations was detected in June–August 2000 by VLBI and GPS techniques. This work confirms this crustal deformation as determined by SLR and vice versa. Analysis of coordinates of these stations shows that this effect was caused by a 4.5-cm displacement of the Tateyama station in the north-east direction. The change in the vertical component was not significant.     

Single GPS receiver time-relative positioning with loop misclosure corrections

Time-relative positioning makes use of observations taken at two different epochs and stations with a single global positioning system (GPS) receiver to determine the position of the unknown station with respect to the known station. The limitation of this method is the degradation over time of the positioning accuracy due to the temporal variation of GPS errors (ionospheric delay, satellite clock corrections, satellite ephemerides, and tropospheric delay). The impact of these errors is significantly reduced by adding to the one-way move from the known to the unknown station, a back move to the known station. A loop misclosure is computed from the coordinates obtained at the known station at the beginning and at the end of the loop, and is used to correct the coordinates of the unknown station. The field tests, presented in this paper, show that using the loop misclosure corrections, time-relative positioning accuracy can be improved by about 60% when using single frequency data, and by about 40% with dual frequency data. For a 4-min processing interval (an 8-min loop) and a 95% probability level, errors remain under 20 cm for the horizontal components and 36 cm for the vertical component with single frequency data; and under 11 cm for the horizontal components and 29 cm for the vertical component with dual frequency data.     

Assessment of ECMWF-derived tropospheric delay models within the EUREF Permanent Network

The Global Positioning System (GPS) observations from the EUREF Permanent Network (EPN) are routinely analyzed by the EPN analysis centers using a tropospheric delay modeling based on standard pressure values, the Niell Mapping Functions (NMF), a cutoff angle of 3° and down-weighting of low elevation observations. We investigate the impact on EPN station heights and Zenith Total Delay (ZTD) estimates when changing to improved models recommended in the updated 2003 International Earth Rotation and Reference Systems Service (IERS) Conventions, which are the Vienna Mapping Functions 1 (VMF1) and zenith hydrostatic delays derived from numerical weather models, or the empirical Global Mapping Functions (GMF) and the empirical Global Pressure and Temperature (GPT) model. A 1-year Global Positioning System (GPS) data set of 50 regionally distributed EPN/IGS (International GNSS Service) stations is processed. The GPS analysis with cutoff elevation angles of 3, 5, and 10° revealed that changing to the new recommended models introduces biases in station heights in the northern part of Europe by 2–3 mm if the cutoff is lower than 5°. However, since large weather changes at synoptic time scales are not accounted for in the empirical models, repeatability of height and ZTD time series are improved with the use of a priori Zenith Hydrostatic Delays (ZHDs) derived from numerical weather models and VMF1. With a cutoff angle of 3°, the repeatability of station heights in the northern part of Europe is improved by 3–4 mm.     

The corridor correction concept

Atmospheric delays are contributors to the GNSS error budget in precise GNSS positioning that can reduce positioning accuracy considerably if not compensated appropriately. Both ionospheric and tropospheric delay corrections can be determined with help of reference stations in active GNSS networks. One approach to interpolate these error terms to the user’s location that is employed in Germany’s SAPOS network is the determination of area correction parameters (ACP, German: “Fl?chenkorrekturparameter—FKP”). A 2D interpolation scheme using data from at least 3 reference stations surrounding the rover is employed. A modification of this method was developed which only makes use of as few as 2 reference stations and provides 1D linear correction parameters along a “corridor” in which the user’s rover is moving. We present the results of a feasibility study portraying results from use of corridor correction parameters for precise RTK-like positioning. The differences to the reference coordinates (3D) attained in average for 1 h of data employing selected network nodes in Germany are between 0.8 and 2.0 cm, which compares well with the traditional area correction method that yields an error of 0.7 up to 1.1 cm.     

Low-degree earth deformation from reprocessed GPS observations

Surface mass variations of low spherical harmonic degree are derived from residual displacements of continuously tracking global positioning system (GPS) sites. Reprocessed GPS observations of 14 years are adjusted to obtain surface load coefficients up to degree n _max = 6 together with station positions and velocities from a rigorous parameter combination. Amplitude and phase estimates of the degree-1 annual variations are partly in good agreement with previously published results, but also show interannual differences of up to 2 mm and about 30 days, respectively. The results of this paper reveal significant impacts from different GPS observation modeling approaches on estimated degree-1 coefficients. We obtain displacements of the center of figure (CF) relative to the center of mass (CM), Δr _CF–CM, that differ by about 10 mm in maximum when compared to those of the commonly used coordinate residual approach. Neglected higher-order ionospheric terms are found to induce artificial seasonal and long-term variations especially for the z-component of Δr _CF–CM. Daily degree-1 estimates are examined in the frequency domain to assess alias contributions from model deficiencies with regard to satellite orbits. Finally, we directly compare our estimated low-degree surface load coefficients with recent results that involve data from the Gravity Recovery and Climate Experiment (GRACE) satellite mission.     

A grid-based tropospheric product for China using a GNSS network

Tropospheric delay accounts for one source of error in global navigation satellite systems (GNSS). To better characterize the tropospheric delays in the temporal and spatial domain and facilitate the safety-critical use of GNSS across China, a method is proposed to generate a grid-based tropospheric product (GTP) using the GNSS network with an empirical tropospheric model, known as IGGtrop. The prototype system generates the GTPs in post-processing and real-time modes and is based on the undifferenced and uncombined precise point positioning (UU-PPP) technique. GTPs are constructed for a grid form (\(2.0{^{\circ }}\times 2.5{^{\circ }}\) latitude–longitude) over China with a time resolution of 5 min. The real-time GTP messages are encoded in a self-defined RTCM3 format and broadcast to users using NTRIP (networked transport of RTCM via internet protocol), which enables efficient and safe transmission to real-time users. Our approach for GTP generation consists of three sequential steps. In the first step, GNSS-derived zenith tropospheric delays (ZTDs) for a network of GNSS stations are estimated using UU-PPP. In the second step, vertical adjustments for the GNSS-derived ZTDs are applied to address the height differences between the GNSS stations and grid points. The ZTD height corrections are provided by the IGGtrop model. Finally, an inverse distance weighting method is used to interpolate the GNSS-derived ZTDs from the surrounding GNSS stations to the location of the grid point. A total of 210 global positioning system (GPS) stations from the crustal movement observation network of China are used to generate the GTPs in both post-processing and real-time modes. The accuracies of the GTPs are assessed against with ERA-Interim-derived ZTDs and the GPS-derived ZTDs at 12 test GPS stations, respectively. The results show that the post-processing and real-time GTPs can provide the ZTDs with accuracies of 1.4 and 1.8 cm, respectively. We also apply the GTPs in real-time kinematic GPS PPP, and the results show that the convergence time of the PPP solutions is shortened. These results confirm that the GTPs can act as an efficient information source to augment GNSS positioning over China.     

电离层延迟高阶项对GPS精密单点定位精度影响分析

随着GPS卫星轨道、钟差及各种误差修正模型的不断精化,静态精密单点定位(PPP)定位精度达到mm级,进行电离层延迟高阶项较小量级的误差改正研究,对改进PPP数据处理策略具有重要的参考价值。本文利用分布在不同地理纬度的5个IGS跟踪站3天的观测数据,对比分析了电离层延迟二阶项、三阶项对GPS观测值精度及静态PPP定位精度的影响。分析结果表明,电离层延迟二阶项、三阶项对GPS观测值精度的影响分别为cm级和mm级,对低纬度地区PPP定位精度的影响大于3 mm,但对中高纬度的测站观测值、定位精度的影响比低纬度地区小很多。