Incorrect modeling of troposphere delays is one of the major error sources for space geodetic techniques such as Global Navigation Satellite Systems (GNSS) or Very Long Baseline Interferometry (VLBI). Over the years, many approaches have been devised which aim at mapping the delay of radio waves from zenith direction down to the observed elevation angle, so-called mapping functions. This paper contains a new approach intended to refine the currently most important discrete mapping function, the Vienna Mapping Functions 1 (VMF1), which is successively referred to as Vienna Mapping Functions 3 (VMF3). It is designed in such a way as to eliminate shortcomings in the empirical coefficients b and c and in the tuning for the specific elevation angle of 3 • . Ray-traced delays of the ray-tracer RADIATE serve as the basis for the calculation of new mapping function coefficients. Comparisons of modeled slant delays demonstrate the ability of VMF3 to approximate the underlying ray-traced delays more accurately than VMF1 does, in particular at low elevation angles. In other words, when requiring highest precision, VMF3 is to be preferable to VMF1. Aside from revising the discrete form of mapping functions, we also present a new empirical model named Global Pressure and Temperature 3 (GPT3) on a 5 • × 5 • as well as a 1 • × 1 • global grid, which is generally based on the same data. Its main components are hydrostatic and wet empirical mapping function coefficients derived from special averaging techniques of the respective (discrete) VMF3 data. In addition, GPT3 also contains a set of meteorological quantities which are adoptedas they stand from their prede-B Daniel Landskron

The Vienna VLBI and Satellite Software (VieVS) is state-of-the-art Very Long Baseline Interferometry (VLBI) analysis software for geodesy and astrometry. VieVS has been developed at Technische Universität Wien (TU Wien) since 2008, where it is used for research purposes and for teaching space geodetic techniques. In the past decade, it has been successfully applied on Very Long Baseline Interferometry (VLBI) observations for the determination of celestial and terrestrial reference frames as well as for the estimation of celestial pole offsets, universal Time (UT1-UTC), and polar motion based on least-squares adjustment. Furthermore, VieVS is equipped with tools for scheduling and simulating VLBI observations to extragalactic radio sources as well as to satellites and spacecraft, features which proved to be very useful for a variety of applications. VieVS is now available as version 3.0 and we do provide the software to all interested persons and institutions. A wiki with more information about VieVS is available at http://vievswiki.geo.tuwien.ac.at/.

Missing or incorrect consideration of azimuthal asymmetry of troposphere delays is a considerable error source in space geodetic techniques such as Global Navigation Satellite Systems (GNSS) or Very Long Baseline Interferometry (VLBI). So-called horizontal troposphere gradients are generally utilized for modeling such azimuthal variations and are particularly required for observations at low elevation angles. Apart from estimating the gradients within the data analysis, which has become common practice in space geodetic techniques, there is also the possibility to determine the gradients beforehand from different data sources than the actual observations. Using ray-tracing through Numerical Weather Models (NWMs), we determined discrete gradient values referred to as GRAD for VLBI observations, based on the standard gradient model by Chen and Herring (J Geophys Res 102(B9):20489-20502, 1997. https://doi.org/10.1029/97JB01739) and also for new, higher-order gradient models. These gradients are produced on the same data basis as the Vienna Mapping Functions 3 (VMF3) (Landskron and Böhm in J Geod, 2017. https://doi.org/10.1007/s00190-017-1066-2), so they can also be regarded as the VMF3 gradients as they are fully consistent with each other. From VLBI analyses of the Vienna VLBI and Satellite Software (VieVS), it becomes evident that baseline length repeatabilities (BLRs) are improved on average by 5% when using a priori gradients GRAD instead of estimating the gradients. The reason for this improvement is that the gradient estimation yields poor results for VLBI sessions with a small number of observations, while the GRAD a priori gradients are unaffected from this. We also developed a new empirical gradient model applicable for any time and location on Earth, which is included in the Global Pressure and Temperature 3 (GPT3) model. Although being able to describe only the systematic component of azimuthal asymmetry and no short-term variations at all, even these empirical a priori gradients slightly reduce (improve) the BLRs with respect to the estimation of gradients. In general, this paper addresses that a priori horizontal gradients are actually more important for VLBI analysis than previously assumed, as particularly the discrete model GRAD as well as the empirical model GPT3 are indeed able to refine and improve the results.

Abstract. In Global Navigation Satellite System (GNSS) tomography, precise information about the tropospheric water vapor distribution is derived from integral measurements like ground-based GNSS slant wet delays (SWDs). Therefore, the functional relation between observations and unknowns, i.e., the signal paths through the atmosphere, have to be accurately known for each station–satellite pair involved. For GNSS signals observed above a 15∘ elevation angle, the signal path is well approximated by a straight line. However, since electromagnetic waves are prone to atmospheric bending effects, this assumption is not sufficient anymore for lower elevation angles. Thus, in the following, a mixed 2-D piecewise linear ray-tracing approach is introduced and possible error sources in the reconstruction of the bended signal paths are analyzed in more detail. Especially if low elevation observations are considered, unmodeled bending effects can introduce a systematic error of up to 10–20 ppm, on average 1–2 ppm, into the tomography solution. Thereby, not only the ray-tracing method but also the quality of the a priori field can have a significant impact on the reconstructed signal paths, if not reduced by iterative processing. In order to keep the processing time within acceptable limits, a bending model is applied for the upper part of the neutral atmosphere. It helps to reduce the number of processing steps by up to 85 % without significant degradation in accuracy. Therefore, the developed mixed ray-tracing approach allows not only for the correct treatment of low elevation observations but is also fast and applicable for near-real-time applications.

No abstract