Application of ray-traced tropospheric slant delays to geodetic VLBI analysis

The correction of tropospheric influences via so-called path delays is critical for the analysis of observations from space geodetic techniques like the very long baseline interferometry (VLBI). In standard VLBI analysis, the a priori slant path delays are determined using the concept of zenith delays, mapping functions and gradients. The a priori use of ray-traced delays, i.e., tropospheric slant path delays determined with the technique of ray-tracing through the meteorological data of numerical weather models (NWM), serves as an alternative way of correcting the influences of the troposphere on the VLBI observations within the analysis. In the presented research, the application of ray-traced delays to the VLBI analysis of sessions in a time span of 16.5 years is investigated. Ray-traced delays have been determined with program RADIATE (see Hofmeister in Ph.D. thesis, Department of Geodesy and Geophysics, Faculty of Mathematics and Geoinformation, Technische Universität Wien. http://resolver.obvsg.at/urn:nbn:at:at-ubtuw:1-3444, 2016) utilizing meteorological data provided by NWM of the European Centre for Medium-Range Weather Forecasts (ECMWF). In comparison with a standard VLBI analysis, which includes the tropospheric gradient estimation, the application of the ray-traced delays to an analysis, which uses the same parameterization except for the a priori slant path delay handling and the used wet mapping factors for the zenith wet delay (ZWD) estimation, improves the baseline length repeatability (BLR) at 55.9% of the baselines at sub-mm level. If no tropospheric gradients are estimated within the compared analyses, 90.6% of all baselines benefit from the application of the ray-traced delays, which leads to an average improvement of the BLR of 1 mm. The effects of the ray-traced delays on the terrestrial reference frame are also investigated. A separate assessment of the RADIATE ray-traced delays is carried out by comparison to the ray-traced delays from the National Aeronautics and Space Administration Goddard Space Flight Center (NASA GSFC) (Eriksson and MacMillan in http://lacerta.gsfc.nasa.gov/tropodelays, 2016) with respect to the analysis performances in terms of BLR results. If tropospheric gradient estimation is included in the analysis, 51.3% of the baselines benefit from the RADIATE ray-traced delays at sub-mm difference level. If no tropospheric gradients are estimated within the analysis, the RADIATE ray-traced delays deliver a better BLR at 63% of the baselines compared to the NASA GSFC ray-traced delays.     

A proposed GPS method with multi-antennae and single receiver

A method based on multi-antennae linked to a common GPS receiver is proposed. The goal of the technique is to improve height determination for baselines a few kilometres in length. The advantage of this technique resides in the elimination of relative clock parameters in the between-antenna single difference observations. Because single difference observations are free of clock errors more geometrical strength remains to determine the baseline components. This statement is valid as long as intercable biases can be carefully calibrated. For millimetre height determination, the intercable calibration must be done at the same level of accuracy. Under this assumption it is shown that in general the height standard deviation improves by a factor of about three compared to standard GPS data processing. With the proposed method, the effect of relative tropospheric zenith delay errors becomes a bit smaller (in absolute value), compared to standard data processing. To absorb this error, a relative tropospheric zenith delay parameter may be estimated. Even with this additional parameter in the solution the height standard deviation remains two times smaller than the results of standard processing techniques (without tropospheric zenith delay parameter), and at least five times smaller than in the results obtained from standard processing including one tropospheric zenith delay parameter.     

区域GPS网对夏季降雨过程中大气可降水量的探测

利用GPS网的观测资料,通过GAMIT软件求得5个测站对流层天顶总延迟,进而求出各测站对流层湿延迟;利用湿延迟与大气可降水量之间的转换关系得到各测站的大气可降水量。将所得GPS-PWV值与同时段探空资料所得的大气可降水量以及地表实际降水量进行对比分析,结果表明:GPS-PWV值与探空资料所得的PWV值比较相符;在降水前后,GPS-PWV有比较明显的变化,降水一般出现在GPS-PWV值迅速增加的4-6h内;实际降水量峰值与GPS-PWV增量大小也有较强的相关性。     

Forecast Vienna Mapping Functions 1 for real-time analysis of space geodetic observations

The Vienna Mapping Functions 1 (VMF1) as provided by the Institute of Geodesy and Geophysics (IGG) at the Vienna University of Technology are the most accurate mapping functions for the troposphere delays that are available globally and for the entire history of space geodetic observations. So far, the VMF1 coefficients have been released with a time delay of almost two days; however, many scientific applications require their availability in near real-time, e.g. the Ultra Rapid solutions of the International GNSS Service (IGS) or the analysis of the Intensive sessions of the International VLBI Service (IVS). Here we present coefficients of the VMF1 as well as the hydrostatic and wet zenith delays that have been determined from forecasting data of the European Centre for Medium-Range Weather Forecasts (ECMWF) and provided on global grids. The comparison with parameters derived from ECMWF analysis data shows that the agreement is at the 1 mm level in terms of station height, and that the differences are larger for the wet mapping functions than for the hydrostatic mapping functions and the hydrostatic zenith delays. These new products (VMF1-FC and hydrostatic zenith delays from forecast data) can be used in real-time analysis of geodetic data without significant loss of accuracy.     

Demonstration of sub-meter GPS orbit determination and 1.5 parts in 108 three-dimensional baseline accuracy

Differential tracking of theGPS satellites in high-earth orbit provides a powerful relative positioning capability, even when a relatively small continental U.S. fiducial tracking network is used with less than one-third of the fullGPS constellation. To demonstrate this capability, we have determined baselines of up to2000 km in North America by estimating high-accuracyGPS orbits and ground receiver positions simultaneously. The2000 km baselines agree with very long baseline interferometry(VLBI) solutions at the level of1.5 parts in10 ⁸ and showrms daily repeatability of0.3–2 parts in10 ⁸. The orbits determined for the most thoroughly trackedGPS satellites are accurate to better than1 m. GPS orbit accuracy was assessed from orbit predictions, comparisons with independent data sets, and the accuracy of the continental baselines determined along with the orbits. The bestGPS orbit strategies included data arcs of at least one week, process noise models for tropospheric fluctuations, estimation ofGPS solar pressure coefficients, and combined processing ofGPS carrier phase and pseudorange data. For data arcs of two weeks, constrained process noise models forGPS dynamic parameters significantly improved the solutions.     

对流层延迟对GPS大地高的影响

分析了对流层延迟的基本特性及常用的处理措施,通过实例表明了不同对流层的处理方法对大地高精度有较大影响:当基线两端高差较大时,即使基线较短,其气象条件也有所差异,仅仅通过差分或模型改正仍不能消除对流层延迟的影响,残余的对流层延迟仍然影响大地高的解算精度,通过对残余对流层延迟进行估计能够明显改善大地高的精度。     

Tropospheric parameters: combination studies based on homogeneous VLBI and GPS data

The combination of tropospheric parameters derived from different space-geodetic techniques has not been of large interest in geodesy so far. However, due to the high correlation between station coordinates and tropospheric parameters, the latter should not be neglected in combinations. This paper deals with the comparison and combination of tropospheric parameters derived from global positioning system (GPS) and very long baseline interferometry (VLBI) observations stemming from a 15-day campaign of continuous VLBI observations in 2002 (CONT02). The observation data of both techniques were processed homogeneously to avoid systematic differences between the solutions. We compared the tropospheric estimates of GPS and VLBI at eight co-location sites and found a very good agreement in the temporal behavior of the tropospheric zenith path delays (ZPD), reflected by correlation factors up to 0.98. Following this, a combination of the tropospheric parameters was performed. We demonstrate that the combination of tropospheric parameters leads to a stabilization of combined station networks. This becomes visible in the improvement of the repeatabilities of the station height components. Furthermore, the potential use of independent data from water vapor radiometers (WVRs) to validate space-technique-derived tropospheric parameters was investigated. Correlation coefficients of 0.95 or better were estimated between the tropospheric parameters of WVR and GPS or VLBI. Additionally, the utility of the tropospheric parameters for validation of local tie vectors was investigated. Both tropospheric zenith delays and tropospheric gradients were found to be very suitable to validate the height component and the horizontal components of the local tie, respectively.     

GPS inferred geocentric reference frame for satellite positioning and navigation

Accurate geocentric three dimensional positioning is of great importance for various geodetic and oceanographic applications. While relative positioning accuracy of a few centimeters has become a reality using Very Long Baseline Interferometry (VLBI), the uncertainty in the offset of the adopted coordinate system origin from the geocenter is still believed to be of the order of one meter. Satellite Laser Ranging (SLR) is capable of determining this offset to better than10 cm, though, because of the limited number of satellites, this requires a long arc of data. The Global Positioning System (GPS) measurements provide a powerful alternative for an accurate determination of this origin offset in relatively short period of time. Two strategies are discussed, the first utilizes the precise relative positions predetermined byVLBI, where as the second establishes a reference frame by holding only one of the tracking sites longitude fixed. Covariance analysis studies indicate that geocentric positioning to an accuracy of a few centimeters can be achieved with just one day of preciseGPS pseudorange and carrier phase data.     

The TropGrid2 standard tropospheric correction model

TropGrid2 is a new version of a tropospheric model that is based on climatology and provides tropospheric propagation delay corrections for standard positioning users without temperature, pressure and humidity measurements. Zenith hydrostatic and wet delays are modeled as special harmonic functions taking seasonal and diurnal variations into consideration. The grid-point values are height-reduced and can be interpolated horizontally to the user position. The database used to derive this model consists of more than 9 years of 3D numerical weather fields of the NOAA NCEP GDAS weather model. We validated this standard model using 10 years of GPS-derived zenith path delays at 290 International GPS service reference stations. The gridded version is accurate at a level of 3.8 cm (root mean square in zenith direction) on global average; the average long-term bias is ?0.3 cm. The standard deviations computed by the model turn out to be slightly too pessimistic for almost all stations under investigation, in contrast to the site-specific version, which is only marginally (1 mm) more accurate on global scale.     

精密单点定位中卫星星历对天顶对流层延迟估计的影响

为了分析不同卫星星历对天顶对流层延迟估计的影响,本文选取不同的卫星星历产品分别进行静态精密单点定位试验,估计天顶对流层延迟,并与IGS发布的天顶对流层延迟产品相比。结果表明,采用最终星历、快速星历和超快星历实测部分时,天顶对流层延迟的平均RMS值分别为4.5mm、4.3mm和4.6mm,估计精度一致。而采用超快星历外推部分时,平均RMS值为6.3mm,估计精度略低。     

GPS超长基线解算的误差特性与精度分析

基于单基站的超长基线定位技术在地壳形变监测、高精度授时等领域具有广泛应用,但仍有诸多因素制约着超长基线解算精度。从观测方程出发,利用单差观测值对长（超长）基线(146～1 724 km)解算中的卫星轨道误差、对流层延迟误差、地球潮汐误差和相位缠绕误差等误差特性进行了详细分析。分析结果表明,当基线小于500 km时广播星历误差可忽略不计;超过500 km时需要采用精密星历,同时需要考虑地球潮汐误差的影响;利用参数估计法同时估计基线两端的天顶对流层延迟误差可获得1～2 cm精度;相位缠绕误差对基线小于2 000 km的解算影响可忽略。基于估计天顶对流层延迟的方法解算了5条长(超长)基线（146 km、491 km、837 km、 1 043 km和1 724 km）。实验结果表明,当基线小于500 km时,采用广播星历可获得水平方向优于0.05 m、高程方向优于0.08 m的定位精度;当基线小于2 000 km时,采用超快速精密星历可获得水平方向优于0.025 m、高程方向优于0.055 m的定位精度。解算的初始收敛时间随着基线长度增加而缩短。     

顾及映射函数误差的对流层延迟两步估计法

利用VLBI、GNSS等空间大地测量数据估计对流层延迟时,通常不考虑映射函数的误差,但这种处理策略最终会对估计结果造成不可忽略的影响。首先,根据多年的气象资料,分析了映射函数误差随高度角的变化特性,然后构建了映射函数误差模型,进而提出了一种顾及映射函数误差的对流层延迟两步估计法。试验表明,本文提出的方法可有效削弱斜路径延迟残差,并在一定程度上改善对流层湿延迟的估计精度。     

A modified mix-differenced approach for estimating multi-GNSS real-time satellite clock offsets

A modified mixed-differenced approach for estimating multi-GNSS real-time clock offsets is presented. This approach, as compared to the earlier presented mixed-differenced approach which uses epoch-differenced and undifferenced observations, further adds a satellite-differenced process. The proposed approach, based on real-time orbit products and a mix of epoch-differenced and satellite-differenced observations to estimate only satellite clock offsets and tropospheric zenith wet delays, has fewer estimated parameters than other approaches, and thus its implementing procedure is efficient and can be performed and extended easily. To obtain high accuracy, the approach involves three steps. First, the high-accuracy tropospheric zenith wet delay of each station is estimated using mixed-differenced carrier phase observations. Second, satellite clock offset changes between adjacent epochs are estimated using also mixed-differenced carrier phase observations. Third, the satellite clock offsets at the initial epoch are estimated using satellite-differenced pseudorange observations. Finally, the initial epoch clock results and clock offset changes are concatenated to obtain the clock results of the current epoch. To validate the real-time satellite clock results, multi-GNSS post-processing clock products from IGS ACs were selected for comparison. From the comparison, the standard deviations of the GPS, GLONASS, BeiDou and Galileo systems clock results are approximately 0.1–0.4 ns, except for the BeiDou GEO satellites. The root mean squares are about 0.4–2.3 ns, which are similar to those of other international real-time products. When the clock estimates were assessed based on a pseudo-kinematic PPP procedure, the positioning accuracies in the East, North and Up components reach 5.6, 5.5 and 7.6 cm, respectively, which meet the centimeter level and are comparable to the application of other products.     

Modeling tropospheric wet delays with national GNSS reference network in China for BeiDou precise point positioning

During past decades, precise point positioning (PPP) has been proven to be a well-known positioning technique for centimeter or decimeter level accuracy. However, it needs long convergence time to get high-accuracy positioning, which limits the prospects of PPP, especially in real-time applications. It is expected that the PPP convergence time can be reduced by introducing high-quality external information, such as ionospheric or tropospheric corrections. In this study, several methods for tropospheric wet delays modeling over wide areas are investigated. A new, improved model is developed, applicable in real-time applications in China. Based on the GPT2w model, a modified parameter of zenith wet delay exponential decay wrt. height is introduced in the modeling of the real-time tropospheric delay. The accuracy of this tropospheric model and GPT2w model in different seasons is evaluated with cross-validation, the root mean square of the zenith troposphere delay (ZTD) is 1.2 and 3.6 cm on average, respectively. On the other hand, this new model proves to be better than the tropospheric modeling based on water-vapor scale height; it can accurately express tropospheric delays up to 10 km altitude, which potentially has benefits in many real-time applications. With the high-accuracy ZTD model, the augmented PPP convergence performance for BeiDou navigation satellite system (BDS) and GPS is evaluated. It shows that the contribution of the high-quality ZTD model on PPP convergence performance has relation with the constellation geometry. As BDS constellation geometry is poorer than GPS, the improvement for BDS PPP is more significant than that for GPS PPP. Compared with standard real-time PPP, the convergence time is reduced by 2–7 and 20–50% for the augmented BDS PPP, while GPS PPP only improves about 6 and 18% (on average), in horizontal and vertical directions, respectively. When GPS and BDS are combined, the geometry is greatly improved, which is good enough to get a reliable PPP solution, the augmentation PPP improves insignificantly comparing with standard PPP.     

Multi-technique comparison of tropospheric zenith delays derived during the CONT02 campaign

In October 2002, 15 continuous days of Very Long Baseline Interferometry (VLBI) data were observed in the Continuous VLBI 2002 (CONT02) campaign. All eight radio telescopes involved in CONT02 were co-located with at least one other space-geodetic technique, and three of them also with a Water Vapor Radiometer (WVR). The goal of this paper is to compare the tropospheric zenith delays observed during CONT02 by VLBI, Global Positioning System (GPS), Doppler Orbitography Radiopositioning Integrated by Satellite (DORIS) and WVR and to compare them also with operational pressure level data from the European Centre for Medium-Range Weather Forecasts (ECMWF). We show that the tropospheric zenith delays from VLBI and GPS are in good agreement at the 3–7 mm level. However, while only small biases can be found for most of the stations, at Kokee Park (Hawaii, USA) and Westford (Massachusetts, USA) the zenith delays derived by GPS are larger by more than 5 mm than those from VLBI. At three of the four DORIS stations, there is also a fairly good agreement with GPS and VLBI (about 10 mm), but at Kokee Park the agreement is only at about 30 mm standard deviation, probably due to the much older installation and type of DORIS equipment. This comparison also allows testing of different DORIS analysis strategies with respect to their real impact on the precision of the derived tropospheric parameters. Ground truth information about the zenith delays can also be obtained from the ECMWF numerical weather model and at three sites using WVR measurements, allowing for comparisons with results from the space-geodetic techniques. While there is a good agreement (with some problems mentioned above about DORIS) among the space-geodetic techniques, the comparison with WVR and ECMWF is at a lower accuracy level. The complete CONT02 data set is sufficient to derive a good estimate of the actual precision and accuracy of each geodetic technique for applications in meteorology.     

利用地基GPS资料分析成都地区大气可降水量

利用成都地区GPS地震监测网的观测资料,通过Bernese软件求解获得CHDU、JYAN、PIXI、QLAI、RENS、ZHJI测站的天顶总延迟和湿延迟,运用湿延迟与大气可降水量之间的转换关系得到各测站大气可降水量,与站点地面实际降雨量及SONDE大气可降水量进行了比较分析。结果表明：地基GPS监测大气可降水量变化与地面实际降雨量有很强的相关性,与SONDE大气可降水量的平均差值为0．43mm,均方偏差为2．56mm。多站点GPS监测大气可降水量联合分析,对于研究区域大气可降水量的变化有一定的意义。     

A case study of using Raman lidar measurements in high-accuracy GPS applications

This paper investigates the impact of rapid small-scale water vapor fluctuations on GPS height determination. Water vapor measurements from a Raman lidar are used for documenting the water vapor heterogeneities and correcting GPS signal propagation delays in clear sky conditions. We use data from four short observing sessions (6 h) during the VAPIC experiment (15 May–15 June 2004). The retrieval of wet delays from our Raman lidar is shown to agree well with radiosonde retrievals (bias and standard deviation (SD) were smaller than 1 and 2.8 mm, respectively) and microwave radiometers (from two different instruments, bias was 6.0/−6.6 mm and SD 1.3/3.8 mm). A standard GPS data analysis is shown to fail in accurately reproducing fast zenith wet delay (ZWD) variations. The ZWD estimates could be improved when mean post-fit phase residuals were removed. Several methodologies for integrating zenith lidar observations into the GPS data processing are also presented. The final method consists in using lidar wet delays for correcting a priori the GPS phase observations and estimating a scale factor for the lidar wet delays jointly with the GPS station position. The estimation of this scale factor allows correcting for a mis-calibration in the lidar data and provides in the same way an estimate of the Raman lidar instrument constant. The agreement of this constant with an independent determination using radiosonde data is at the level of 1–4%. The lidar wet delays were derived by ray-tracing from zenith pointing measurements: further improvement in GPS positioning is expected from slant path lidar measurements that would properly account for water vapor anisotropy.     

Ray-traced slant factors for mitigating the tropospheric delay at the observation level

Three-dimensional ray tracing through a numerical weather model has been applied to a global precise point positioning (PPP) campaign for modeling both the elevation angle- and azimuth-dependence of the tropospheric delay. Rather than applying the ray-traced slant delays directly, the delay has been parameterized in terms of slant factors, which are applied in a similar manner to traditional mapping functions, but which can account for the azimuthal asymmetry of the delay. Five strategies are considered: (1) Vienna Mapping Functions 1 (VMF1) and estimation of a residual zenith delay parameter; (2) VMF1, estimation of a residual zenith delay and estimation of two tropospheric gradient parameters; (3) three-dimensional ray-traced slant factors and estimation of a residual zenith delay; (4) using only ray-traced slant factors and no estimation of any tropospheric parameters and; (5) using both ray-traced slant factors and estimating a residual zenith delay and two tropospheric gradient parameters. The use of the ray-traced slant factors (solution 3) showed a 3.8% improvement in the repeatability of the up component when compared to the assumption of a symmetric atmosphere (solution 1), while the estimation of two tropospheric gradient parameters gave the best results showing an 7.6% improvement over solution 1 in the up component. Solution 4 performed well in the horizontal domain, allowing for sub-centimeter repeatability but the up component was degraded due to deficiencies in the modeling of the zenith delay, particularly for stations located at equatorial latitudes. The magnitude of the differences in the mean coordinates between solution 2 and solution 3, and the strong correlation with the differences between the north component and the ray-traced gradients (coefficient of correlation of 0.83), as well as the impact of observation geometry on the gradient solution indicate that the use of the ray-traced slant factors could have an implication on the realization of reference frames. The estimated tropospheric products from the PPP solutions were compared to those derived from ray tracing. For the zenith delay, a root mean square (RMS) of 5.4 mm was found, while for the gradient terms, a correlation coefficient of 0.46 for the N–S and 0.42 for the E–W was found for the north–south and east–west components, suggesting that there are still important differences in the gradient parameters which could be due to either errors in the NWM or to non-tropospheric error sources leaking into the PPP-estimated gradients.     

Using external tropospheric corrections to improve GNSS positioning of hot-air balloon

High accurate global navigation satellite systems (GNSS) require to correct a signal delay caused by the troposphere. The delay can be estimated along with other unknowns or introduced from external models. We assess the impact of the recently developed augmentation tropospheric model on real-time kinematic precise point positioning (PPP). The model is based on numerical weather forecast and thus reflects the actual state of weather conditions. Using the G-Nut/Geb software, we processed GNSS and meteorological data collected during the experiment using a hot-air balloon flying up to an altitude of 2000 m. We studied the impacts of random walk noise setting of zenith total delay (ZTD) on estimated parameters and the mutual correlations, the use of external tropospheric corrections, the use of data from a single or dual GNSS constellation and the use of Kalman filter and backward smoothing processing methods. We observed a significant negative correlation of the estimated rover height and ZTD which depends on constraining ZTD estimates. Such correlation caused a degraded performance of both parameters when estimated simultaneously, in particular for a single GNSS constellation. The impact of ZTD constraining reached up to 50-cm differences in the rover height. Introducing external tropospheric corrections improved the PPP solution regarding: (1) shortened convergence, (2) better overall robustness, particularly, in case of degraded satellite geometry, (3) less adjusted parameters with lower correlations. The numerical weather model-driven PPP resulted in 9–12- and 5–6-cm uncertainties in the rover altitude using the Kalman filter and the backward smoothing, respectively. Compared to standard PPP, it indicates better performance by a factor of 1–2 depending on the availability of GNSS constellations, the troposphere constraining and the processing strategy.     

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.