In GNSS data processing, the station height, receiver clock and tropospheric delay (ZTD) are highly correlated to each other. Although the zenith hydrostatic delay of the troposphere can be provided with sufficient accuracy, zenith wet delay (ZWD) has to be estimated, which is usually done in a random walk process. Since ZWD temporal variation depends on the water vapor content in the atmosphere, it seems to be reasonable that ZWD constraints in GNSS processing should be geographically and/or time dependent. We propose to take benefit from numerical weather prediction models to define optimum random walk process noise. In the first approach, we used archived VMF1-G data to calculate a grid of yearly and monthly means of the difference of ZWD between two consecutive epochs divided by the root square of the time lapsed, which can be considered as a random walk process noise. Alternatively, we used the Global Forecast System model from National Centres for Environmental Prediction to calculate random walk process noise dynamically in real-time. We performed two representative experimental campaigns with 20 globally distributed International GNSS Service (IGS) stations and compared real-time ZTD estimates with the official ZTD product from the IGS. With both our approaches, we obtained an improvement of up to 10% in accuracy of the ZTD estimates compared to any uniformly fixed random walk process noise applied for all stations.

The tropospheric delay is one of the major error sources in precise point positioning (PPP), affecting the accuracy and precision of estimated coordinates and convergence time, which raises demand for a reliable tropospheric model, suitable to support PPP. In this study, we investigate the impact of three tropospheric models and mapping functions regarding position accuracy and convergence time. We propose a routine to constrain the tropospheric estimates, which we implemented in the in-house developed real-time PPP software. We take advantage of the high spatial resolution (4 × 4 km²) numerical weather prediction Weather Research and Forecasting (WRF) model and near real-time GNSS data combined by the least-squares collocation estimation to reconstruct the tropospheric delays. We also present mapping functions calculated from the WRF model using the ray-tracing technique. The performance tests are conducted on 14 Polish EUREF Permanent Network (EPN) stations during 3 weeks of different tropospheric conditions: calm, standard and severe. We consider six GNSS data processing variants, including two commonly used variants using a priori ZTD and mapping functions from UNB3m and VMF1-FC models, one with a priori ZTD and mapping functions calculated directly from WRF model and three variants using the aforementioned mapping functions but with ZTD model based on GNSS and WRF data used as a priori troposphere and to constrain tropospheric estimates. The application of a high-resolution GNSS/WRF-based ZTD model and mapping functions results in the best agreement with the official EPN coordinates. In both static and kinematic modes, this approach results in an average reduction of 3D bias by 20 and 10 mm, respectively, but an increase of 3D SDs by 1.5 and 4 mm, respectively. The application of high-resolution tropospheric model also shortens the convergence time, for example, for a 10 cm convergence level, from 67 to 58 min for the horizontal components and from 79 to 63 min for the vertical component.

An extensive validation of line-of-sight tropospheric slant total delays (STD) from Global Navigation Satellite Systems (GNSS), ray tracing in numerical weather prediction model (NWM) fields and microwave water vapour radiometer (WVR) is presented. Ten GNSS reference stations, including collocated sites, and almost 2 months of data from 2013, including severe weather events were used for comparison. Seven institutions delivered their STDs based on GNSS observations processed using 5 software programs and 11 strategies enabling to compare rather different solutions and to assess the impact of several aspects of the processing strategy. STDs from NWM ray tracing came from three institutions using three different NWMs and ray-tracing software. Inter-techniques evaluations demonstrated a good mutual agreement of various GNSS STD solutions compared to NWM and WVR STDs. The mean bias among GNSS solutions not considering post-fit residuals in STDs was −0.6 mm for STDs scaled in the zenith direction and the mean standard deviation was 3.7 mm. Standard deviations of comparisons between GNSS and NWM ray-tracing solutions were typically 10 mm ± 2 mm (scaled in the zenith direction), depending on the NWM model and the GNSS station. Comparing GNSS versus WVR STDs reached standard deviations of 12 mm ± 2 mm also scaled in the zenith direction. Impacts of raw GNSS post-fit residuals and cleaned residuals on optimal reconstructing of GNSS STDs were evaluated at inter-technique comparison and for GNSS at collocated sites. The use of raw post-fit residuals is not generally recommended as they might contain strong systematic effects, as demonstrated in the case of station LDB0. Simplified STDs reconstructed only from estimated GNSS tropospheric parameters, i.e. without applying post-fit residuals, performed the best in all the comparisons; however, it obviously missed part of tropospheric signals due to non-linear temporal and spatial variations in the troposphere. Although the post-fit residuals cleaned of visible systematic errors generally showed a slightly worse performance, they contained significant tropospheric signal on top of the simplified model. They are thus recommended for the reconstruction of STDs, particularly during high variability in the troposphere. Cleaned residuals also showed a stable performance during ordinary days while containing promising information about the troposphere at low-elevation angles.

The Global Navigation Satellite System (GNSS) is commonly recognized by its all-weather capability. However, observations depend on atmospheric conditions which requires the induced tropospheric delay to be estimated as an unknown parameter. In the following study, we investigate the impact of intense weather events on GNSS estimates. GNSS slant total delays (STD) in Precise Point Positioning technique (PPP) strategy were calculated for stations in southwest Poland in a 56 days period covering several heavy precipitation cases. The corresponding delays retrieved from Weather Research and Forecasting (WRF) model by a ray-tracing technique considered only gaseous parts of the atmosphere. The discrepancies are correlated with rain rates and cloud type products from remote sensing platforms. Positive correlation is found as well as GNSS estimates tend to be systematically larger than modeled delays. Mean differences mapped to the zenith direction are showed to vary between 10 mm and 30 mm. The magnitude of discrepancies follows the intensity of phenomena, especially for severe weather events. Results suggest that effects induced by commonly neglected liquid and solid water terms in the troposphere modeling should be considered in precise GNSS applications for the atmosphere monitoring. The state-of-art functional model applied in GNSS processing strategies shows certain deficits. Estimated tropospheric delays with gradients and post-fit residuals could be replaced by a loosely constrained solution without loss of quality.