Accuracy assessment of the National Geodetic Survey’s OPUS-RS utility

OPUS-RS is a rapid static form of the National Geodetic Survey’s On-line Positioning User Service (OPUS). Like OPUS, OPUS-RS accepts a user’s GPS tracking data and uses corresponding data from the U.S. Continuously Operating Reference Station (CORS) network to compute the 3-D positional coordinates of the user’s data-collection point called the rover. OPUS-RS uses a new processing engine, called RSGPS, which can generate coordinates with an accuracy of a few centimeters for data sets spanning as little as 15 min of time. OPUS-RS achieves such results by interpolating (or extrapolating) the atmospheric delays, measured at several CORS located within 250 km of the rover, to predict the atmospheric delays experienced at the rover. Consequently, standard errors of computed coordinates depend highly on the local geometry of the CORS network and on the distances between the rover and the local CORS. We introduce a unitless parameter called the interpolative dilution of precision (IDOP) to quantify the local geometry of the CORS network relative to the rover, and we quantify the standard errors of the coordinates, obtained via OPUS-RS, by using functions of the form

Geocenter variations derived from a combined processing of LEO- and ground-based GPS observations

GNSS observations provided by the global tracking network of the International GNSS Service (IGS, Dow et al. in J Geod 83(3):191–198, 2009) play an important role in the realization of a unique terrestrial reference frame that is accurate enough to allow a detailed monitoring of the Earth’s system. Combining these ground-based data with GPS observations tracked by high-quality dual-frequency receivers on-board low earth orbiters (LEOs) is a promising way to further improve the realization of the terrestrial reference frame and the estimation of geocenter coordinates, GPS satellite orbits and Earth rotation parameters. To assess the scope of the improvement on the geocenter coordinates, we processed a network of 53 globally distributed and stable IGS stations together with four LEOs (GRACE-A, GRACE-B, OSTM/Jason-2 and GOCE) over a time interval of 3 years (2010–2012). To ensure fully consistent solutions, the zero-difference phase observations of the ground stations and LEOs were processed in a common least-squares adjustment, estimating all the relevant parameters such as GPS and LEO orbits, station coordinates, Earth rotation parameters and geocenter motion. We present the significant impact of the individual LEO and a combination of all four LEOs on the geocenter coordinates. The formal errors are reduced by around 20% due to the inclusion of one LEO into the ground-only solution, while in a solution with four LEOs LEO-specific characteristics are significantly reduced. We compare the derived geocenter coordinates w.r.t. LAGEOS results and external solutions based on GPS and SLR data. We found good agreement in the amplitudes of all components; however, the phases in x- and z-direction do not agree well.     

Removing attitude-related variations in the line-of-sight for kinematic GPS positioning

We are presenting a method for removing attitude-related variations (ATTRVs) in dynamic 1 Hz GPS positioning. The ATTRVs are separated into vertical and horizontal components. These result from the translational and rotational motions of the vehicle, which is a marine research vessel in our case. We have developed new observation equations that use corrected pseudoranges and carrier phases to account for ATTRVs. In the present contribution, we are only focusing on the vertical signals. These modeled ATTRVs are included as corrections in the line-of-sight (LoS) to each GPS satellite in the RINEX data sets using a new software called RNXATTCOR. Precise IGS sp3-orbits are used as inputs together with precise lever arm coordinates of the onboard GPS antennas, observations from the marine inertial navigation system and a priori 3D position of the vessel. The corrected RINEX data sets are then processed using kinematic processing or sequential processing in Precise Point Positioning (PPP) mode. The reduction of the standard deviation from a running mean in the final processed GPS time series is about 95%. The current method is being proposed for marine geodesy science applications such as along-track calibration/validation of altimetry products, and also for applications related to offshore kinematic high precision GPS positioning such as drilling, offshore platforms stability, marine pipeline positioning, earthquake monitoring and tsunami detection.     

A double-differenced cycle slip detection and repair method for GNSS CORS network

The difficulty to detect and repair cycle slip of carrier phase measurements is a key limit for continuously high accuracy of GNSS positioning and navigation services. We propose an automated cycle slip detection and repair method for data preprocessing of a CORS network. The method jointly uses double-differenced (DD) geometry-free (GF) combination and ionospheric-free observation corrected for the computed geometrical distance (IF-OMC) to estimate the cycle slips in dual-frequency observations. The DD GF combination, which is only affected by the ionospheric residual, can be used to detect cycle slips with high reliability except for special pairs such as (77, 60) on GPS L1/L2 frequencies. The detection principle of the IF-OMC observable is such that there is a large discontinuity related to the previous epoch when cycle slips occur at the present epoch. The disadvantages of these two combinations can be overcome employing the proposed detection method. The cycle slip pair (77, 60) has no effect on the GF combination, while a change of 14.65 m is derived from GPS L1/L2 observations using the IF-OMC algorithm. Using pre-determined station coordinates as precise values, we found that the accuracy of the DD IF-OMC combination was 18 mm for a 200-km CORS baseline. Therefore, cycle slips in dual-frequency observations can be correctly and uniquely determined using DD GF and IF-OMC equations. The proposed method was verified by adding simulated cycle slips in observations collected from the CORS network under a quiet ionosphere and shown to be effective. Moreover, the method was assessed with observations made during intense ionospheric activity, which generated extensive cycle slips. The results show that the algorithm can detect and repair all cycle slips apart from two exceptions relating to long data gaps.     

CORS在线定位用户服务系统

首先对基于CORS的在线定位用户服务系统功能进行了较细致的分析,将其概括为7项基本功能和4项扩展功能。在此基础上提出了系统架构方案,将系统划分为3个部件和20个子部件,并对各组成部分的功能作了较详细的阐述。最后,针对用户数据上传与存储、基线解算和基线向量网平差等核心功能,提出了较具体的设计方案,并对软件实现与系统集成方法进行了阐述。     

全球导航卫星系统的新进展

本文综合介绍了于 2 0 0 4年 9月 2 1日至 2 4日在美国加州举行的“全球导航卫星系统” 2 0 0 4年年会(GNSS2 0 0 4 )会议的主要议题 ,并对其中我们可能关切的方面进行了重点介绍。①美国的GPS连续运行站网(CORS)。CORS由美国大地测量局 (NGS)主持运行。用户可以通过NGS网络 ,获得用户的GPS待定点相邻的CORS站 (三个以上 )的GPS相应载波相位和码距 ,以支持用户的GPS准实时或后处理定位。NGS也可以为用户通过网络提供GPS定位计算服务 ,这一服务可以在用户提供待定点的观测资料后的几个小时内完成 ,称为NGS的在线GPS定位服务。CORS目前在美国已有 5 0 0余个站。②GPS系统的进展。GPSⅡR型卫星从体形和功能方面都比较优秀 ,使GPS卫星在轨的位置误差显著降低 ,测距精度提高近一倍 ,目前GPSⅡR型卫星截止至 2 0 0 4年 1月 1日时有 9颗在轨。③利用L1 ,L2频道的GPS空基增强系统 (WAAS)。在美国大部分地区WAAS系统的水平精度可达 1～ 2m ,垂直精度可达 2～ 3m。④GPS信号的重构。美国已发展了一种高度逼真的和适应各种情况的虚拟GPS信号系统 ,这种虚拟发射装置可以是陆基的 ,空基的 ,或者星基的。GPS接收机可以利用这一虚拟的GPS信号进行精密定位。⑤Galileo卫星导航系统运行的准备工作。欧洲空间局已经重新和?     

Site-specific multipath characteristics of global IGS and CORS GPS sites

The site-specific multipath characteristics of 217 Global Positioning System (GPS) sites worldwide were analyzed using the variability of the post-fit phase residuals. Because the GPS satellite constellation returns to the same configuration in a sidereal day (23 h 56 min 4 s), the multipath repeats on that period. However, daily GPS position estimates are usually based on the solar day. When several days of GPS data are processed, this steady change in the orientation of the satellite constellation with respect to the station manifests itself in the form of patterns in the post-fit phase residuals which shift by 3 min 56 s per day. It was found that the mean root mean square of the time-shifted post-fit phase residuals is highly dependent on the GPS antenna type. The conclusions derived from the analysis of the time-shifted post-fit residuals were verified by performing a cross-correlation of the post-fit residuals across many days for selected sites.     

Accuracy analysis for the NAD 83 geodetic reference system

The North American Datum of 1983 (NAD 83) provides horizontal coordinates for more than 250,000 geodetic stations. These coordinates were derived by a least squares adjustment of existing terrestrial and space-based geodetic data. For pairs of first order stations with interstation distances between 10km and 100km, therms discrepancy between distances derived fromNAD 83 coordinates and distances derived from independentGPS data may be suitably approximated by the empirical rulee=0.008 K^0.7 where e denotes therms discrepancy in meters and K denotes interstation distance in kilometers. For the same station pairs, therms discrepancy in azimuth may be approximated by the empirical rule e=0.020 K^0.5. Similar formulas characterize therms discrepancies for pairs involving second and third order stations. Distance and orientation accuracies, moreover, are well within adopted standards. While these expressions indicate that the magnitudes of relative positional accuracies depend on station order, absolute positional accuracies are similar in magnitude for first, second, and third order stations. Adjustment residuals reveal a few local problems with theNAD 83 coordinates and with the weights assigned to certain classes of observations.     

W0: an estimate in the Finnish Height Datum N60, epoch 1993.4, from twenty-five GPS points of the Baltic Sea Level Project

Based upon a data set of 25 points of the Baltic Sea Level Project, second campaign 1993.4, which are close to mareographic stations, described by (1) GPS derived Cartesian coordinates in the World Geodetic Reference System 1984 and (2) orthometric heights in the Finnish Height Datum N60, epoch 1993.4, we have computed the primary geodetic parameter W ₀(1993.4) for the epoch 1993.4 according to the following model. The Cartesian coordinates of the GPS stations have been converted into spheroidal coordinates. The gravity potential as the additive decomposition of the gravitational potential and the centrifugal potential has been computed for any GPS station in spheroidal coordinates, namely for a global spheroidal model of the gravitational potential field. For a global set of spheroidal harmonic coefficients a transformation of spherical harmonic coefficients into spheroidal harmonic coefficients has been implemented and applied to the global spherical model OSU 91A up to degree/order 360/360. The gravity potential with respect to a global spheroidal model of degree/order 360/360 has been finally transformed by means of the orthometric heights of the GPS stations with respect to the Finnish Height Datum N60, epoch 1993.4, in terms of the spheroidal “free-air” potential reduction in order to produce the spheroidal W ₀(1993.4) value. As a mean of those 25 W ₀(1993.4) data as well as a root mean square error estimation we computed W ₀(1993.4)=(6 263 685.58 ± 0.36) kgal × m. Finally a comparison of different W ₀ data with respect to a spherical harmonic global model and spheroidal harmonic global model of Somigliana-Pizetti type (level ellipsoid as a reference, degree/order 2/0) according to The Geodesist's Handbook 1992 has been made. Received: 7 November 1996 / Accepted: 27 March 1997     

GPS sidereal filtering: coordinate- and carrier-phase-level strategies

Multipath error is considered one of the major errors affecting GPS observations. One can benefit from the repetition of satellite geometry approximately every sidereal day, and apply filtering to help minimize this error. For GPS data at 1 s interval processed using a double-difference strategy, using the day-to-day coordinate or carrier-phase residual autocorrelation determined with a 10-h window leads to the steadiest estimates of the error-repeat lag, although a window as short as 2 h can produce an acceptable value with > 97% of the optimal lag’s correlation. We conclude that although the lag may vary with time, such variation is marginal and there is little advantage in using a satellite-specific or other time-varying lag in double-difference processing. We filter the GPS data either by stacking a number of days of processed coordinate residuals using the optimum “sidereal” lag (23 h 55 m 54 s), and removing these stacked residuals from the day in question (coordinate space), or by a similar method using double-difference carrier-phase residuals (observational space). Either method results in more consistent and homogeneous set of coordinates throughout the dataset compared with unfiltered processing. Coordinate stacking reduces geometry-related repeating errors (mainly multipath) better than carrier-phase residual stacking, although the latter takes less processing time to achieve final filtered coordinates. Thus, the optimal stacking method will depend on whether coordinate precision or computational time is the over-riding criterion.     

CGCS2000框架维持方法分析

针对目前国内各GNSS应用在归算到2000国家大地坐标系（CGCS2000）时常采用强制符合及拟稳平差方法存在的问题及局限性,利用国家测绘地理信息局收集的全国近1 800个连续运行参考站（Continuously Operating Reference Stations,CORS）连续观测1个月的数据进行数据处理得到的基线解,分别基于强基准拟稳平差方式及顾及板块运动的速度场归算方法获得测站在CGCS2000下的坐标,分析两者的差异及误差分布情况,指出对于省级CORS站纳入CGCS2000的情况,采用强制符合拟稳平差方法所基于的基准站必须考虑板块的运动趋势进行CGCS2000归算,否则平差结果与CGCS2000最大差异将达到dm量级。     

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.     

Variance-covariance estimation of GPS Networks

Summary It is quite easy to estimate the variance-covariance (VCV) matrix for single session surveys or local networks, but difficult where these local networks are combined together to form a regional network. Our main aim is to develop an appropriate VCV model to combine all the different types of networks, either global, regional or local. By careful estimation and combination of the individual VCVs of the local networks, we can form a unique VCV for local, regional and global networks. Different techniques are used to derive appropriate models for the variancecovariance components of the Global Positioning System (GPS) networks. The VCV models were estimated using homogeneous and heterogeneous data. The variance-covariance components are empirically derived using (a) the covariance of the observations of homogeneous data, (b) a combination of the covariance of the observationsP ^–1 and the covariance of the signal componentsC _ss (for either homogeneous and/or heterogeneous data), (c) only the variances are used to determine the variancecovariance, their covariances being zeros. We compare the solutions of the VCV developed for homogeneous and/or heterogeneous data with other developed VCVs. It was observed that the derived VCV model for the combined homogeneous and/or heterogeneous data of case (b) gives the best estimates in all cases.     

Accuracy of GPS-derived relative positions as a function of interstation distance and observing-session duration

 Ten days of GPS data from 1998 were processed to determine how the accuracy of a derived three-dimensional relative position vector between GPS antennas depends on the chord distance (denoted L) between these antennas and on the duration of the GPS observing session (denoted T). It was found that the dependence of accuracy on L is negligibly small when (a) using the `final' GPS satellite orbits disseminated by the International GPS Service, (b) fixing integer ambiguities, (c) estimating appropriate neutral-atmosphere-delay parameters, (d) 26 km ≤ L ≤ 300 km, and (e) 4 h ≤T ≤ 24 h. Under these same conditions, the standard error for the relative position in the north–south dimension (denoted S _n and expressed in mm) is adequately approximated by the equation S _n =k _n /T ^0.5 with k _n =9.5 ± 2.1 mm · h^0.5 and T expressed in hours. Similarly, the standard errors for the relative position in the east–west and in the up-down dimensions are adequately approximated by the equations S _e =k _e /T ^0.5 and S _u =k _u /T ^0.5, respectively, with k _e =9.9 ± 3.1 mm · h^0.5 and k _u =36.5 ± 9.1 mm · h^0.5. Received: 5 February 2001 / Accepted: 14 May 2001     

Combined processing of the IGS and the IGEX network

 Since the beginning of the International Global Navigation Satellite System (GLONASS) Experiment, IGEX, in October 1998, the Center for Orbit Determination in Europe (CODE) has acted as an analysis center providing precise GLONASS orbits on a regular basis. In CODE's IGEX routine analysis the Global Positioning System (GPS) orbits and Earth rotation parameters are introduced as known quantities into the GLONASS processing. A new approach is studied, where data from the IGEX network are combined with GPS observations from the International GPS Service (IGS) network and all parameters (GPS and GLONASS orbits, Earth rotation parameters, and site coordinates) are estimated in one processing step. The influence of different solar radiation pressure parameterizations on the GLONASS orbits is studied using different parameter subsets of the extended CODE orbit model. Parameterization with three constant terms in the three orthogonal directions, D, Y, and X (D = direction satellite–Sun, Y = direction of the satellite's solar panel axis), and two periodic terms in the X-direction, proves to be adequate for GLONASS satellites. As a result of the processing it is found that the solar radiation pressure effect for the GLONASS satellites is significantly different in the Y-direction from that for the GPS satellites, and an extensive analysis is carried out to investigate the effect in detail. SLR observations from the ILRS network are used as an independent check on the quality of the GLONASS orbital solutions. Both processing aspects, combining the two networks and changing the orbit parameterization, significantly improve the quality of the determined GLONASS orbits compared to the orbits stemming from CODE's IGEX routine processing. Received: 10 May 2000 / Accepted: 9 October 2000     

Validating ocean tide loading models using GPS

Ocean tides cause periodic deformations of the Earths surface, also referred to as ocean tide loading (OTL). Tide-induced displacements of the Earths crust relying on OTL models are usually taken into account in GPS (Global Positioning System) data analyses. On the other hand, it is also possible to validate OTL models using GPS analyses. The following simple approach is used to validate OTL models. Based on a particular model, instantaneous corrections of the site coordinates due to OTL are computed. Site-specific scale factors, f, for these corrections are estimated in a standard least-squares adjustment process of GPS observations together with other relevant parameters. A resulting value of f close to unity indicates a good agreement of the model with the actual site displacements. Such scale factors are computed for about 140 globally distributed IGS (International GPS Service) tracking sites. Three OTL models derived from the ocean tide models FES95.2.1, FES99, and GOT00.2 are analyzed. As expected, the most reliable factors are estimated for sites with a large loading effect. In general, the scaling factors have a value close to unity and no significant differences between the three ocean tide models could be observed. It is found that the validation approach is easy to apply. Without requiring much additional effort for a global and self-consistent GPS data analysis, it allows detection of general model misfits on the basis of a large number of globally distributed sites. For detailed validation studies on OTL models, the simultaneous estimation of amplitudes and phases for the main contributing partial tides within a GPS parameter adjustment process would provide more detailed answers.     

浅谈CORS站观测墩高量测方法

目前,利用多基站网络RTK技术建立的连续运行卫星定位服务综合系统(Continuous Operational Reference System,缩写为CORS)已成为各省GPS应用的发展重点。CORS系统是卫星定位技术、计算机网络技术、数字通讯技术等高新科技多方位、深度的产物。观测墩高度的测量精度将直接影响CORS站系统高程精度,本文结合LNCORS站观测墩高测量方法及在准备阶段研究成果,阐述CORS站观测墩高测量方法及精度分析。     

Combined Earth orientation parameters based on homogeneous and continuous VLBI and GPS data

The CONT02 campaign is of great interest for studies combining very long baseline interferometry (VLBI) with other space-geodetic techniques, because of the continuously available VLBI observations over 2 weeks in October 2002 from a homogeneous network. Especially, the combination with the Global Positioning System (GPS) offers a broad spectrum of common parameters. We combined station coordinates, Earth orientation parameters (EOPs) and troposphere parameters consistently in one solution using technique- specific datum-free normal equation systems. In this paper, we focus on the analyses concerning the EOPs, whereas the comparison and combination of the troposphere parameters and station coordinates are covered in a companion paper in Journal of Geodesy. In order to demonstrate the potential of the VLBI and GPS space-geodetic techniques, we chose a sub-daily resolution for polar motion (PM) and universal time (UT). A consequence of this solution set-up is the presence of a one-to-one correlation between the nutation angles and a retrograde diurnal signal in PM. The Bernese GPS Software used for the combination provides a constraining approach to handle this singularity. Simulation studies involving both nutation offsets and rates helped to get a deeper understanding of this singularity. With a rigorous combination of UT1–UTC and length of day (LOD) from VLBI and GPS, we showed that such a combination works very well and does not suffer from the systematic effects present in the GPS-derived LOD values. By means of wavelet analyses and the formal errors of the estimates, we explain this important result. The same holds for the combination of nutation offsets and rates. The local geodetic ties between GPS and VLBI antennas play an essential role within the inter-technique combination. Several studies already revealed non-negligible discrepancies between the terrestrial measurements and the space-geodetic solutions. We demonstrate to what extent these discrepancies propagate into the combined EOP solution.     

Non-linear motions of Australian geodetic stations induced by non-tidal ocean loading and the passage of tropical cyclones

We investigate daily and sub-daily non-tidal oceanic and atmospheric loading (NTOAL) in the Australian region and put an upper bound on potential site motion examining the effects of tropical cyclone Yasi that crossed the Australian coast in January/February 2011. The dynamic nature of the ocean is important, particularly for northern Australia where the long-term scatter due to daily and sub-daily oceanic changes increases by 20–55 % compared to that estimated using the inverted barometer (IB) assumption. Correcting the daily Global Positioning System (GPS) time series for NTOAL employing either a dynamic ocean model or the IB assumption leads to a reduction of up to 52 % in the weighted scatter of daily coordinate estimates. Differences between the approaches are obscured by seasonal variations in the GPS precision along the northern coast. Two compensating signals during the cyclone require modelling at high spatial and temporal resolution: uplift induced by the atmospheric depression, and subsidence induced by storm surge. The latter dominates ( \(>\) 135 %) the combined net effect that reaches a maximum of 14 mm, and 10 mm near the closest GPS site TOW2. Here, 96 % of the displacement is reached within 15 h due to the rapid transit of cyclones and the quasi-linear nature of the coastline. Consequently, estimating sub-daily NTOAL is necessary to properly account for such a signal that can be 3.5 times larger than its daily-averaged value. We were unable to detect the deformation signal in 2-hourly GPS processing and show that seasonal noise in the Austral summer dominates and precludes GPS detection of the cyclone-related subsidence.     

Mapping GPS positional errors using spatial linear mixed models

Nowadays, GPS receivers are very reliable because of their good accuracy and precision; however, uncertainty is also inherent in geospatial data. Quality of GPS measurements can be influenced by atmospheric disturbances, multipathing, synchronization of clocks, satellite geometry, geographical features of the observed region, low broadcasting coverage, inadequate transmitting formats, or human or instrumental unknown errors. Assuming that the scenario and technical conditions that can influence the quality of GPS measurements are optimal, that functional and stochastic models that process the signals to a geodetic measurement are correct, and that all the GPS observables are taken in the same conditions, it is still possible to estimate the positional errors as the difference between the real coordinates and those measured by the GPS. In this paper, three spatial linear mixed models, one for each axis, are used for modelling real-time kinematic GPS accuracy and precision, of a multiple-reference-station network in dual-frequency with carrier phase measurements. Along the paper, the proposed models provide an estimate of the “accuracy” in terms of bias defined as the difference between real coordinates and measured coordinates after being processed and “precision” through the standard errors of the estimated differences. This is done using ten different transmitting formats. Mapping and quantifying these differences can be interesting for users and GPS professionals. The performance of these models is illustrated by mapping positional error estimates within the whole region of Navarre, Spain. Sampled data have been taken in 54 out of the 211 geodetic vertex points of this region. Maps show interesting error patterns depending on transmitting formats, the different axes, and the geographical characteristics of the region. Higher differences are found in regions with bad broadcasting coverage, due to the presence of mountains and high degree of humidity.