We developed operators to assimilate Global Navigation Satellite System (GNSS) Zenith Total Delays (ZTDs) and horizontal delay gradients into a numerical weather model. In this study we experiment with refractivity fields derived from the Global Forecast System (GFS) available with a horizontal resolution of 0.5°. We begin our investigations with simulated observations. In essence, we extract the tropospheric parameters from the GFS analysis, add noise to mimic observation errors and assimilate the simulated observations into the GFS 24h forecast valid at the same time. We consider three scenarios: (1) the assimilation of ZTDs (2) the assimilation of horizontal delay gradients and (3) the assimilation of both ZTDs and horizontal delay gradients. The impact is measured by utilizing the refractivity fields. We find that the assimilation of the horizontal delay gradients in addition to the ZTDs improves the refractivity field around 800 hPa. When we consider a single station there is a clear improvement when horizontal delay gradients are assimilated in addition to the ZTDs because the horizontal delay gradients contain information that is not contained in the ZTDs. On the other hand, when we consider a dense station network there is not a significant improvement when horizontal delay gradients are assimilated in addition to the ZTDs because the horizontal delay gradients do not contain information that is not already contained in the ZTDs. Finally, we replace simulated by real observations, that is, tropospheric parameters from a Precise Point Positioning solution provided with the G-Nut/Tefnut software, in order to show that the GFS 24h forecast is indeed improved when GNSS horizontal delay gradients are assimilated in addition to GNSS ZTDs; for the considered station (Potsdam, Germany) and period (June and July, 2017) we find an improvement in the retrieved refractivity of up to 4%.

The Benchmark data set collected within the European COST Action ES1206 has aimed to support the development and validation of advanced Global Navigation Satellite System (GNSS) tropospheric products, in particular high-resolution zenith delays and tropospheric gradients. In this work we utilize this unique data set to show that the interpolation of GNSS Zenith Wet Delays (ZWDs) can be improved by utilizing tropospheric gradients. To do this we first prove the concept with simulated observations, that is, zenith delays and tropospheric gradients derived from a Numerical Weather Model. We show how tropospheric gradients can be converted to ZWD gradients. Then the ZWD gradients together with the ZWDs at selected reference stations are used in an inverse distance weighting interpolation scheme to estimate the ZWD at some target station. For a station configuration with an average station distance of 50 km in Germany and a period of two months (May and June 2013), we find an improvement of 20% in interpolated ZWDs when tropospheric gradients are taken into account. Next, we replace the simulated by real observations, that is, zenith delays and tropospheric gradients from a Precise Point Positioning (PPP) solution provided with the G-Nut/Tefnut analysis software. Here we find an improvement of 10% in interpolated ZWDs when tropospheric gradients are taken into account.

We developed a new strategy for a synchronous generation of real-time (RT) and near real-time (NRT) tropospheric products. It exploits the precise point positioning method with Kalman filtering and backward smoothing, both supported by real-time orbit and clock products. The strategy can be optimized for the latency or the accuracy of NRT production. In terms of precision, it is comparable to the traditional NRT network solution using deterministic models in the least-square adjustment. Both RT and NRT solutions provide a consistent set of tropospheric parameters such as zenith total delays, horizontal tropospheric gradients and slant delays, all with a high resolution and optimally exploiting all observations from available GNSS multi-constellations. As the new strategy exploits RT processing, we assessed publicly precise RT products and results of RT troposphere monitoring. The backward smoothing applied for NRT solution, when using an optimal latency of 30 min, reached an improvement of 20% when compared to RT products. Additionally, multi-GNSS solutions provided more accurate (by 25%) tropospheric parameters, and the impact will further increase when constellations are complete and supported with precise models and products. The new strategy is ready to replace our NRT contribution to the EUMETNET EIG GNSS Water Vapour Programme (E-GVAP) and effectively support all modern multi-GNSS tropospheric products.

The tropospheric delays estimated from the Global Navigation Satellite System (GNSS) have been proven to be an efficient product for monitoring variations of water vapor, which plays an important role in meteorology applications. The operational GNSS water vapor monitoring system is currently based on the Global Positioning System (GPS) and GLObal NAvigation Satellite System(GLONASS) dual-frequency observations. The Galileo satellite navigation system has been evolving continuously, and on 11 February 2019, the constellation reached 22 active satellites, achieving a capability of standalone Precise Point Positioning (PPP) and tropospheric estimation that is global in scope. This enhancement shows a 37% improvement if the precision of the Galileo-only zenith tropospheric delay, while we may anticipate further benefits in terms of tropospheric gradients and slant delays in the future if an optimal multi-constellation and multi-frequency processing strategy is used. First, we analyze the performance of the multi-frequency troposphere estimates based on the PPP raw observation model by comparing it with the standard ionosphere-free model. The performance of the Galileo-only tropospheric solution is then validated with respect to GPS-only solution using 48 globally distributed Multi-GNSS Experiment (MGEX) stations. The averaged bias and standard deviations are −0.3 and 5.8 mm when only using GPS satellites, respectively, and 0.0 and 6.2 mm for Galileo, respectively, when compared to the International GNSS Service (IGS) final Zenith Troposphere Delay(ZTD) products. Using receiver antenna phase center corrections from the corresponding GPS dual-frequency observations does not affect the Galileo PPP ambiguity float troposphere solutions. These results demonstrate a comparable precision achieved for both Galileo-only and GPS-only ZTD solutions, however, horizontal tropospheric gradients, estimated from standalone GPS and Galileo solutions, still show larger discrepancies, mainly due to their being less Galileo satellites than GPS satellites. Including Galileo E1, E5a, E5b, and E5 signals, along with their proper observation weighting, show the benefit of multi-frequency observations, further improving the ZTD precision by 4% when compared to the dual-frequency raw observation model. Overall, the presented results demonstrate good prospects for the application of multi-frequency Galileo observations for the tropospheric parameter estimates.