The tropospheric delay is one of many error sources that affect the Global Navigation Satellite System (GNSS) positioning solutions. The widely used troposphere models assume a homogeneous atmosphere so that only the zenith delay needs to be determined and is mapped through an elevation-dependent mapping function. This procedure is to reduce the computational burden and keep the positioning model full-rank. However, this assumption fails for a realistic description of the troposphere, which is always asymmetrical at a certain elevation angle, especially during a weather event when the weather conditions are very complex. These imperfectly modelled tropospheric delays may influence the positioning accuracy and integer ambiguity resolution performance. In this case, this contribution aims to investigate the effects of the model errors due to the asymmetrical troposphere on GNSS estimations. The Numerical Weather Prediction (NWP) model is applied to generate the actual ray-tracing tropospheric delay in Western Europe, and the tropospheric model errors are calculated in a normal weather condition and a weather event condition by comparing the slant delay calculated from the NWP model and the mapping function. Case studies on the same GNSS station are conducted in two weather conditions: a normal troposphere condition and a weather event with heavy rainfall. The results based on the case studies show that the troposphere in the normal weather condition is nearly homogeneous that the azimuthal-dependent discrepancies of the tropospheric delay are less than 1cm at a very low elevation angle; meanwhile, the discrepancies between different azimuthal angles can reach to more than 25cm in the weather event. A single-frequency Single Point Positioning (SPP) model and a Precise Point Positioning (PPP) model that preserves the integer property of ambiguity are chosen for studying the estimation biases caused by the troposphere model errors. It turns out that almost all horizontal positioning biases of SPP and PPP are less than 1cm in the normal weather condition; however, the scales of the horizontal and 3D biases are concentrated in 1 to 10cm in the weather event for these two models. This contribution also contains the study of the actual integer ambiguity resolution success rate in the presence of the tropospheric model errors by applying the Monte Carlo simulation, and the success rates of PPP in the normal weather condition are consistent with the theoretical values calculated with the ideal troposphere which is totally symmetrical. However, the actual success rates in the weather event are extremely low at some epochs due to the tropospheric model errors, which means that wrong fixing may occur since the theoretical values cannot take into account these model errors. Note that the horizontal tropospheric gradients are not involved in the processing, which means that an optimistic performance might be expected if the gradients are considered.

To monitor temporal variations of the Earth’s gravity field and mass transport in the Earth’s system, data from gravity recovery and climate experiment (GRACE) satellite mission and its successor GRACE Follow-On (GFO) are used. To fill in the temporal gap between these missions, other satellites’ kinematic orbits derived from GPS-based high-low satellite-to-satellite tracking data may be considered. However, it is well known that kinematic orbits are highly sensitive to various systematic errors. These errors are responsible for a non-stationary noise in the kinematic orbits, which is difficult to handle. As a result, the quality of the obtained gravity field solutions is reduced. In this research, we propose to apply an epoch-difference (ED) scheme in the context of the classical dynamic approach to gravity field recovery. Compared to the traditional undifferenced (UD) scheme, the ED scheme is able to mitigate constant or slowly varying systematic errors. To demonstrate the added value of the ED scheme, three sets of monthly gravity field solutions produced from 6 years of GRACE kinematic orbits are compared: two sets produced in-house (with the ED and UD scheme), and a set produced with the undifferenced scheme in the frame of the short-arc approach (Zehentner and Mayer-Gürr in J Geodesy 90(3):275–286, 2015. https://doi.org/10.1007/s00190-015-0872-7). As a reference, we use state-of-the-art ITSG-Grace2018 monthly gravity field solutions. A comparison in the spectral domain shows that the gravity field solutions suffer from a lower noise level when the ED scheme is applied, particularly at low-degree terms, with cumulative errors up to degree 20 being reduced by at least 20%. In the spatial domain, the ED scheme notably reduces noise levels in the mass anomalies recovered. In addition, the signals in terms of mean mass anomalies in selected regions become closer to those inferred from ITSG-Grace2018 solutions, while showing no evidence of any damping, when the ED scheme is used. We conclude that the proposed ED scheme is preferable for time-varying gravity field modeling, as compared to the traditional UD scheme. Our findings may facilitate, among others, bridging the gap between GRACE and GFO satellite mission.

This contribution implements the Kriging interpolation in predicting the tropospheric wet delays using global navigation satellite system networks. The predicted tropospheric delays can be used in strengthening the precise point positioning models and numerical weather prediction models. In order to evaluate the performances of the Kriging interpolation, a sparse network with 8 stations and a dense network with 19 stations from continuously operating reference stations (CORS) of the Netherlands are selected as the reference. In addition, other 15 CORS stations are selected as users, which are divided into three blocks: 5 stations located approximately in the center of the networks, 5 stations on the edge of the networks and 5 stations outside the networks. The zenith tropospheric wet delays are estimated at the network and user stations through the ionosphere-free positioning model; meanwhile, the predicted wet delays at the user stations are generated by the Kriging interpolation in the use of the tropospheric estimations at the network. The root mean square errors (RMSE) are calculated by comparing the predicted wet delays and estimated wet delays at the same user station. The results show that RMSEs of the stations inside the network are at a sub-centimeter level with an average value of 0.74 cm in the sparse network and 0.69 cm in the dense network. The stations on edge and outside the network can also achieve 1-cm level accuracy, which overcomes the limitation that accurate interpolations can only be attained inside the network. This contribution also presents an insignificant improvement of the prediction accuracy from the sparse network to the dense network over 1-year’s data processing and a seasonal effect on the tropospheric wet delay predictions.

The benefits of an increased number of global navigation satellite systems (GNSS) in space have been confirmed for the robustness and convergence time of standard precise point positioning (PPP) solutions, as well as improved accuracy when (most of) the ambiguities are fixed. Yet, it is still worthwhile to investigate fast and high-precision GNSS parameter estimation to meet user needs. This contribution focuses on integer ambiguity resolution-enabled Precise Point Positioning (PPP-RTK) in the use of the observations from four global navigation systems, i.e., GPS (Global Positioning System), Galileo (European Global Navigation Satellite System), BDS (Chinese BeiDou Navigation Satellite System), and GLONASS (Global’naya Navigatsionnaya Sputnikova Sistema). An undifferenced and uncombined PPP-RTK model is implemented for which the satellite clock and phase bias corrections are computed from the data processing of a group of stations in a network and then provided to users to help them achieve integer ambiguity resolution on a single receiver by calibrating the satellite phase biases. The dataset is recorded in a local area of the GNSS network of the Netherlands, in which 12 stations are regarded as the reference to generate the corresponding corrections and 21 as the users to assess the performance of the multi-GNSS PPP-RTK in both kinematic and static positioning mode. The results show that the root-mean-square (RMS) errors of the ambiguity float solutions can achieve the same accuracy level of the ambiguity fixed solutions after convergence. The combined GNSS cases, on the contrary, reduce the horizontal RMS of GPS alone with 2 cm level to GPS + Galileo/GPS + Galileo + BDS/GPS + Galileo + BDS + GLONASS with 1 cm level. The convergence time benefits from both multi-GNSS and fixing ambiguities, and the performances of the ambiguity fixed solution are comparable to those of the multi-GNSS ambiguity float solutions. For instance, the convergence time of GPS alone ambiguity fixed solutions to achieve 10 cm three-dimensional (3D) positioning accuracy is 39.5 min, while it is 37 min for GPS + Galileo ambiguity float solutions; moreover, with the same criterion, the convergence time of GE ambiguity fixed solutions is 19 min, which is better than GPS + Galileo + BDS + GLONASS ambiguity float solutions with 28.5 min. The experiments indicate that GPS alone occasionally suffers from a wrong fixing problem; however, this problem does not exist in the combined systems. Finally, integer ambiguity resolution is still necessary for multi-GNSS in the case of fast achieving very-high-accuracy positioning, e.g., sub-centimeter level.

The Gravity Recovery And Climate Experiment (GRACE) mission has achieved a quantum leap in knowledge of the Earth's gravity field. However, current gravity field solutions still cannot reach the prelaunch baseline accuracy. One of the reasons for that is the presence of colored noise in GRACE data, which is typically ignored in the classical dynamic approach to gravity field modeling. In this research, we propose to account for colored noise in the classical dynamic approach by applying the frequency-dependent data weighting (FDDW) scheme, so that enhanced estimates of gravity field solutions are produced. The monthly solutions are compared with those produced using the standard least squares adjustment without a data weighting scheme. The comparison is performed in both spectral and spatial domains, showing the positive effect of the FDDW scheme in all considered cases. For instance, the cumulative geoid height errors up to degree 96 are reduced by 18%. In the spatial domain, the FDDW scheme lowers noise level in mass changes over the oceans, Mississippi river basin, and Greenland by 20, 38, and 23%, respectively, when compared to the without a data weighting scheme. In addition, the consistency of mass changes over the Mississippi and Congo river basins with those inferred from the state-of-the-art hydrology model WaterGAP is substantially improved when the FDDW scheme is applied. These results indicate that modeling colored noise in the GRACE data allows to significantly improve the recovered monthly solutions. This finding is likely applicable also to the GRACE Follow-On mission.