Abstract. For more than 2 decades the IGS (International GNSS Service) ionosphere associated analysis centers (IAACs) have provided global maps of the vertical total electron content (VTEC). In general, the representation of a 2-D or 3-D function can be performed by means of a series expansion or by using a discretization technique. While in the latter case, pixels or voxels are usually chosen for a spherical function such as VTEC, for a series expansion spherical harmonics (SH) are primarily used as basis functions. The selection of the best suited approach for ionosphere modeling means a trade-off between the distribution of available data and their possibility of representing ionospheric variations with high resolution and high accuracy. Most of the IAACs generate global ionosphere maps (GIMs) based on SH expansions up to the spectral degree n=15 and provide them with a spatial resolution of 2.5∘×5∘ with respect to the latitudinal and longitudinal directions, respectively, and a temporal sampling interval of 2 h. In recent years, it has frequently been claimed that the spatial resolution of the VTEC GIMs has to be increased to a spatial resolution of 1∘×1∘ and to a temporal sampling interval of about 15 min. Enhancing the grid resolution means an interpolation of VTEC values for intermediate points but with no further information about variations in the signal. n=15 in the SH case, for instance, corresponds to a spatial sampling of 12∘×12∘. Consequently, increasing the grid resolution concurrently requires an extension of the spectral content, i.e., to choose a higher SH degree value than 15. Unlike most of the IAACs, the VTEC modeling approach at Deutsches Geodätisches Forschungsinstitut der Technischen Universität München (DGFI-TUM) is based on localizing basis functions, namely tensor products of polynomial and trigonometric B-splines. In this way, not only can data gaps be handled appropriately and sparse normal equation systems be established for the parameter estimation procedure, a multi-scale representation (MSR) can also be set up to determine GIMs of different spectral content directly, by applying the so-called pyramid algorithm, and to perform highly effective data compression techniques. The estimation of the MSR model parameters is finally performed by a Kalman filter driven by near real-time (NRT) GNSS data. Within this paper, we realize the MSR and create multi-scale products based on B-spline scaling, wavelet coefficients and VTEC grid values. We compare these products with different final and rapid products from the IAACs, e.g., the SH model from CODE (Berne) and the voxel solution from UPC (Barcelona). In contrast to the abovementioned products, DGFI-TUM's products are based solely on NRT GNSS observations and ultra-rapid orbits. Nevertheless, we can conclude that the DGFI-TUM's high-resolution product (“othg”) outperforms all products used within the selected time span of investigation, namely September 2017.

The Kalman filter (KF) is widely applied in (ultra) rapid and (near) real-time ionosphere modeling to meet the demand on ionosphere products required in many applications extending from navigation and positioning to monitoring space weather events and naturals disasters. The requirement of a prior definition of the stochastic models attached to the measurements and the dynamic models of the KF is a drawback associated with its standard implementation since model uncertainties can exhibit temporal variations or the time span of a given test data set would not be large enough. Adaptive methods can mitigate these problems by tuning the stochastic model parameters during the filter run-time. Accordingly, one of the primary objectives of our study is to apply an adaptive KF based on variance component estimation to compute the global Vertical Total Electron Content (VTEC) of the ionosphere by assimilating different ionospheric GNSS measurements. Secondly, the derived VTEC representation is based on a series expansion in terms of compactly supported B-spline functions. We highlight the morphological similarity of the spatial distributions and the magnitudes between VTEC values and the corresponding estimated B-spline coefficients. This similarity allows for deducing physical interpretations from the coefficients. In this context, an empirical adaptive model to account for the dynamic model uncertainties, representing the temporal variations of VTEC errors, is developed in this work according to the structure of B-spline coefficients. For the validation, the differential slant total electron content (dSTEC) analysis and a comparison with Jason-2/3 altimetry data are performed. Assessments show that the quality of the VTEC products derived by the presented algorithm is in good agreement, or even more accurate, with the products provided by IGS ionosphere analysis centers within the selected periods in 2015 and 2017. Furthermore, we show that the presented approach can be applied to different ionosphere conditions ranging from very high to low solar activity without concerning time-variable model uncertainties, including measurement error and process noise of the KF because the associated covariance matrices are computed in a self-adaptive manner during run-time.

The ionosphere is an ionized part of the upper atmosphere, spanning from 60 to around 1,000 km in altitude (Hargreaves, 1992). It arises mainly due to the photoionization effects from the solar extreme ultraviolet (EUV) radiation and charged energetic-particle precipitation (Kivelson & Russell, 1995). Generally, the ionosphere is strongly coupled with the thermosphere (Astafyeva, 2019). The latter supplies the neutral particles that can be ionized, and plays a crucial role in the interplay between the production (source) and recombination (loss) processes. The ionosphere affects the propagation of the Global Navigation Satellite System (GNSS) signals by introducing frequency-dependent delays. Unlike the neutral atmosphere, which can cause errors in navigation and positioning in the order of several meters, ionospheric effects can yield uncertainties of up to  E 100 m (e.g., Hernández-Pajares et al., 2011;Petit & Luzum, 2010). Ionospheric delays are inversely related to the square of carrier frequency, and directly proportional to electron density integrated along the ray path (e.g., Goss et al., 2019;Hobiger & Jakowski, 2017).Electron density distribution in the ionosphere strongly depends on altitude and can be divided into several layers, originally identified from ionograms: the D-layer (60-90 km altitude), E-layer (90-130 km), and F-layer (above 130 km), which can be subdivided into F1 and F2 layers (e.g., Astafyeva, 2019). The dominant contribution to electron density profiles comes from the peak of the F2 layer, generally located between

The ionosphere is one of the largest error sources in GNSS (Global Navigation Satellite Systems) applications and can cause up to several meters of error in positioning. Especially for single-frequency users, who cannot correct the ionospheric delay, the information about the state of the ionosphere is mandatory. Dual- and multi-frequency GNSS users, on the other hand, can correct the ionospheric effect on their observations by linear combination. However, real-time applications such as autonomous driving or precision farming, require external high accuracy corrections for fast convergence. Mostly, this external information is given in terms of grids or coefficients of the vertical total electron content (VTEC). Globally distributed GNSS stations of different networks, such as the network of the International GNSS Services (IGS), provide a large number of multi-frequency observations which can be used to determine the state of the ionosphere. These data are used to generate Global Ionosphere Maps (GIM). Due to the inhomogeneous global distribution of GNSS real-time stations and especially due to the large data gaps over oceanic areas, the global VTEC models are usually limited in their spatial and spectral resolution. Most of the GIMs are mathematically based on globally defined radial basis functions, i.e., spherical harmonics (SH), with a maximum degree of 15 and provided with a spatial resolution of 2.5 ° × 5 ° in latitude and longitude, respectively. Regional GNSS networks, however, offer dense clusters of observations, which can be used to generate regional VTEC solutions with a higher spectral resolution. In this study, we introduce a two-step model (TSM) comprising a global model as the first step and a regional model as the second step. We apply polynomial and trigonometric B-spline functions to represent the global VTEC. Polynomial B-splines are used for modelling the finer structures of VTEC within selected regions, i.e., the densification areas. The TSM provides both, a global and a regional VTEC map at the same time. In order to study the performance, we apply the developed approach to hourly data of the global IGS network as well as the EUREF network of the European region for St. Patrick storm in March 2015. For the assessment of the generated maps, we use the dSTEC analysis and compare both maps with different global and regional products from the IGS Ionosphere Associated Analysis Centers, e.g., the global product from CODE (Berne, Switzerland) and from UPC (Barcelona, Spain), as well as the regional maps from ROB (Brussels, Belgium). The assessment shows a significant improvement of the regional VTEC representation in the form of the generated TSM maps. Among all other products used for comparison, the developed regional one is of the highest accuracy within the selected time span. Since the numerical tests are performed using hourly data with a latency of one to two hours, the presented procedure is seen as an intermediate step for the generation of high precision regional real-time corrections for modern applications.

Global ionosphere maps (GIM) computed from dual-frequency GNSS measurements have been widely used for monitoring ionosphere as well as providing ionospheric corrections in Space Geodesy since 1998. This work focuses on a comprehensive study of the influence of time resolution on GIM performance. One and a half solar cycle of the IGS GIM with higher time resolution and accuracy (the UPC-IonSAT Quarter-of-an-hour time resolution Rapid GIM, AKA uqrg) has been taken as baseline to downsample them to all possible sub-daily temporal resolutions. The performance of the resulting GIMs has been assessed by directly comparing with external Vertical Total Electron Content (VTEC) measurements from Jason altimeters over oceanic regions. In order to perform a complete assessment and analysis of involved GIMs, the influence of geographical position, solar and geomagnetic activity were also taken into account during more than one solar cycle. In addition, to have a clear view at the smaller time resolutions, a more accurate assessment, the dSTEC test based on external GNSS measurements not used in the GIM generation, was also done during two solstice and two equinox days in 2015 over continental regions. The assessment shows that discrepancy among GIMs with different time resolutions becomes more apparent at low latitudes and also at high solar-geomagnetic activity. The results also suggest that the accuracy for GIMs with time resolution smaller or equal to 60 minutes is consistent