ALBERT

Description: 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.

Description: In this paper, we present results of the second reprocessing of all data from 1996 to 2014 from all stations in International Association of Geodesy (IAG) Reference Frame Sub-Commission for Europe (EUREF) Permanent Network (EPN) as performed at the Geodetic Observatory Pecný (GOP). While the original goal of this research was to ultimately contribute to the realization of a new European Terrestrial Reference System (ETRS), we also aim to provide a new set of GNSS (Global Navigation Satellite System) tropospheric parameter time series with possible applications to climate research. To achieve these goals, we improved a strategy to guarantee the continuity of these tropospheric parameters and we prepared several variants of troposphere modelling. We then assessed all solutions in terms of the repeatability of coordinates as an internal evaluation of applied models and strategies and in terms of zenith tropospheric delays (ZTDs) and horizontal gradients with those of the ERA-Interim numerical weather model (NWM) reanalysis. When compared to the GOP Repro1 (first EUREF reprocessing) solution, the results of the GOP Repro2 (second EUREF reprocessing) yielded improvements of approximately 50 and 25 % in the repeatability of the horizontal and vertical components, respectively, and of approximately 9 % in tropospheric parameters. Vertical repeatability was reduced from 4.14 to 3.73 mm when using the VMF1 mapping function, a priori ZHD (zenith hydrostatic delay), and non-tidal atmospheric loading corrections from actual weather data. Raising the elevation cut-off angle from 3 to 7° and then to 10° increased RMS from coordinates' repeatability, which was then confirmed by independently comparing GNSS tropospheric parameters with the NWM reanalysis. The assessment of tropospheric horizontal gradients with respect to the ERA-Interim revealed a strong sensitivity of estimated gradients to the quality of GNSS antenna tracking performance. This impact was demonstrated at the Mallorca station, where gradients systematically grew up to 5 mm during the period between 2003 and 2008, before this behaviour disappeared when the antenna at the station was changed. The impact of processing variants on long-term ZTD trend estimates was assessed at 172 EUREF stations with time series longer than 10 years. The most significant site-specific impact was due to the non-tidal atmospheric loading followed by the impact of changing the elevation cut-off angle from 3 to 10°. The other processing strategy had a very small or negligible impact on estimated trends.

Description: An analysis of processing settings impact 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 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 with the G-Nut/Tefnut software. The impact of the gradient mapping function, elevation cut-off angle, GNSS constellation 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 prediction models. Generally, all the solutions in the post-processing mode provided a robust tropospheric gradient estimation with a clear relation to real weather conditions. The quality of tropospheric gradient estimates in real-time mode mainly depends on the actual quality of the real-time orbits and clocks. Best results were achieved using the 3° elevation angle cut-off and a combined GPS+GLONASS constellation. Systematic effects of up to 0.3mm were observed in estimated tropospheric gradients when using different gradient mapping functions which depend on the applied observation elevation-dependent weighting. While the latitudinal troposphere tilting causes a systematic difference in the north gradient component on a global scale, large local wet gradients pointing to a direction of increased humidity cause systematic differences in both gradient components depending on the gradient direction.

Description: 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%.

Description: 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.

Description: In this paper, we present results of the 2nd reprocessing of all data from 1996 to 2014 from all stations in the European GNSS permanent network as performed at the Geodetic Observatory Pecný (GOP). While the original goal of this research was to ultimately contribute to new realization of the European terrestrial reference system, we also aim to provide a new set of GNSS tropospheric parameter time series with possible applications to climate research. To achieve these goals, we enhanced a strategy to guarantee the continuity of these tropospheric parameters and we prepared several variants of troposphere modelling. We then assessed all solutions in terms of the repeatability of coordinates as an internal evaluation of applied models and strategies, and in terms of zenith tropospheric delays (ZTD) and horizontal gradients with those of ERA-Interim numerical weather model (NWM) reanalysis. When compared to the GOP Repro1 solution, the results of the GOP Repro2 yielded improvements of approximately 50 % and 25 % in the repeatability of the horizontal and vertical components, respectively, and of approximately 9 % in tropospheric parameters. Vertical repeatability was reduced from 4.14 mm to 3.73 mm when using the VMF1 mapping function, a priori ZHD, and non-tidal atmospheric loading corrections from actual weather data. Raising the elevation angle cut-off from 3° to 7° and then to 10° increased RMS from coordinates’ repeatability, which was then confirmed by independently assessing GNSS tropospheric parameters with the NWM reanalysis. The assessment of tropospheric horizontal gradients with respect to the ERA-Interim revealed a strong sensitivity of estimated gradients to the quality of GNSS antenna tracking performance. This impact was demonstrated at the Mallorca station, where gradients systematically grew up to 5 mm during the period between 2003 and 2008, before this behaviour disappeared when the antenna at the station was changed.