Abstract. Initial objectives and design of the Benchmark campaign organized within the European COST Action ES1206

Abstract. 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 intertechnique 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.

Abstract. An analysis of processing settings impacts on estimated tropospheric gradients is presented. The study is based on the benchmark data set collected within the COST GNSS4SWEC action with observations from 430 Global Navigation Satellite Systems (GNSS) reference stations in central Europe for May and June 2013. Tropospheric gradients were estimated in eight different variants of GNSS data processing using precise point positioning (PPP) with the G-Nut/Tefnut software. The impacts of the gradient mapping function, elevation cut-off angle, GNSS constellation, observation elevation-dependent weighting and real-time versus post-processing mode were assessed by comparing the variants by each to other and by evaluating them with respect to tropospheric gradients derived from two numerical weather models (NWMs). Tropospheric gradients estimated in post-processing GNSS solutions using final products were in good agreement with NWM outputs. The quality of high-resolution gradients estimated in (near-)real-time PPP analysis still remains a challenging task due to the quality of the real-time orbit and clock corrections. Comparisons of GNSS and NWM gradients suggest the 3∘ elevation angle cut-off and GPS+GLONASS constellation for obtaining optimal gradient estimates provided precise models for antenna-phase centre offsets and variations, and tropospheric mapping functions are applied for low-elevation observations. Finally, systematic errors can affect the gradient components solely due to the use of different gradient mapping functions, and still depending on observation elevation-dependent weighting. A latitudinal tilting of the troposphere in a global scale causes a systematic difference of up to 0.3 mm in the north-gradient component, while large local gradients, usually pointing in a direction of increasing humidity, can cause differences of up to 1.0 mm (or even more in extreme cases) in any component depending on the actual direction of the gradient. Although the Bar-Sever gradient mapping function provided slightly better results in some aspects, it is not possible to give any strong recommendation on the gradient mapping function selection.

GPS tomography has been investigated since 2000 as an attractive tool for retrieving the 3D field of water vapour and wet refractivity. However, this observational technique still remains a challenging task that requires improvement of its methodology. This was the purpose of this study, and for this, GPS data from the Australian Continuously Operating Research Station (CORS) network during a severe weather event were used. Sensitivity tests and statistical cross-comparisons of tomography retrievals with independent observations from radiosonde and radio-occultation profiles showed improved results using the presented methodology. The initial conditions, which were associated with different time-convergence of tomography inversion, play a critical role in GPS tomography. The best strategy can reduce the normalised root mean square (RMS) of the tomography solution by more than 3 with respect to radiosonde estimates. Data stacking and pseudo-slant observations can also significantly improve tomography retrievals with respect to non-stacked solutions. A normalised RMS improvement up to 17% in the 0-8 km layer was found by using 30 min data stacking, and RMS values were divided by 5 for all the layers by using pseudo-observations. This result was due to a better geometrical distribution of mid-and low-tropospheric parts (a 30% coverage improvement). Our study of the impact of the uncertainty of GPS observations shows that there is an interest in evaluating tomography retrievals in comparison to independent external measurements and in estimating simultaneously the quality of weather forecasts. Finally, a comparison of multi-model tomography with numerical weather prediction shows the relevant use of tomography retrievals to improving the understanding of such severe weather conditions. Global Positioning System (GPS) tomography considers the use of slant-integrated estimates, wet delays, or corresponding water vapour content estimated from the data records of ground-based GPS stations to respectively retrieve the 3D field of wet refractivity or water vapour density, as introduced by [1,2]. Comparisons of tomography retrievals with other techniques (e.g., those using a water vapour radiometer, radiosonde, raman lidar, or atmospheric emitted radiance interferometer) and with numerical weather models have shown relevant results and an encouraging understanding of meteorological situations ([3-17]) The resolution and configuration/geometry of the network of GPS stations are critical parameters with which to obtain the best scenario for applying GPS tomography to retrieve water vapour density or wet refractivity. These fields can be ideally retrieved for meteorological applications with a horizontal resolution of a few kilometres, a vertical resolution of~500 m in the lower troposphere, and a vertical resolution of~2 km in the upper troposphere, with a time resolution of 15 to 5 min. However, to obtain this remarkable geometrical resolution, data from a dense homogeneously distributed network of GPS stations (e.g., a netw...

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.

Abstract. An extensive validation of line-of-sight tropospheric Slant Total Delays (STD) from Global Navigation Satellite Systems (GNSS), ray-tracing in Numerical Weather Prediction Models (NWM) fields and microwave Water Vapour Radiometer (WVR) is presented. Ten GNSS reference stations and almost two months of data from 2013, including severe weather events, entered the comparison. Seven institutions delivered their STDs based on GNSS observations processed using five software and eleven strategies. STDs from NWM ray-tracing came from three institutions using three different NWM models. Results show generally a very good mutual agreement among all solutions from all the techniques. The mean bias (over all stations) between the GNSS solution selected as reference, which did not use post-fit residuals in STDs, and all other GNSS solutions without post-fit residuals is −0.6 mm for STDs scaled in the zenith direction, and the corresponding mean standard deviation is 3.7 mm. Standard deviations of comparisons between GNSS a NWM ray-tracing solutions are typically 10 mm ± 2 mm (scaled in the zenith direction), depending on the NWM model and the particular station considered. When comparing GNSS versus WVR STDs, standard deviations reached 12 mm ± 2 mm, as scaled in zenith direction. Moreover, the influence of adding raw GNSS post-fit residuals, as well as residuals screened out of systematic effects, to STDs was studied. It was found that adding raw post-fit residuals always led to lower quality of GNSS STDs while the situation was not that straightforward after the post-fit residuals cleaning.

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.