Alberto García-Rigo

Global prediction of the vertical total electron content of the ionosphere based on GPS data

[1] Although vertical total electron content (VTEC) forecasting is still an open subject of research, the use of predictions of the ionospheric state at a scale of several days is an area of increased interest. A global VTEC forecast product for two days ahead, which is based exclusively on actual Global Positioning System (GPS) data, has been developed in the frame of the International Global Navigation Satellite Systems (GNSS) Service (IGS) Ionospheric Working Group (IGS Iono‐WG). The UPC ionospheric VTEC prediction model is based on the Discrete Cosine Transform (DCT), which is widely used in image compression (for instance, in JPEG format). Additionally, a linear regression module is used to forecast the time evolution of each of the DCT coefficients. The use of the DCT coefficients is justified because they represent global features of the whole two‐ dimensional VTEC map/image. Also, one can therefore introduce prior information affecting the VTEC, for instance, smoothness or the distribution of relevant features in different directions. For this purpose, the use of a long time series of final/rapid UPC VTEC maps is required. Currently, the UPC Predicted product is being automatically generated in test mode and is made available through the main IGS server for public access. This product is also used to generate two days ahead preliminary combined IGS Predicted product. Finally, the results presented in this work suggest that the two days ahead UPC Predicted product could become an official IGS product in the near future.

Distribution and mitigation of higher‐order ionospheric effects on precise GNSS processing

Higher-order ionospheric effects (I2+) are one of the main limiting factors in very precise Global Navigation Satellite Systems (GNSS) processing, for applications where millimeter accuracy is demanded. This paper summarizes a comprehensive study of the I2+ effects in range and in GNSS precise products such as receiver position and clock, tropospheric delay, geocenter offset, and GNSS satellite position and clock. All the relevant higher-order contributions are considered: second and third orders, geometric bending, and slant total electron content (dSTEC) bending (i.e., the difference between the STEC for straight and bent paths). Using a realistic simulation with representative solar maximum conditions on GPS signals, both the effects and mitigation errors are analyzed. The usage of the combination of multifrequency L band observations has to be rejected due to its increased noise level. The results of the study show that the main two effects in range are the second-order ionospheric and dSTEC terms, with peak values up to 2 cm. Their combined impacts on the precise GNSS satellite products affects the satellite Z coordinates (up to +1 cm) and satellite clocks (more than ±20 ps). Other precise products are affected at the millimeter level. After correction the impact on all the precise GNSS products is reduced below 5 mm. We finally show that the I2+ impact on a Precise Point Positioning (PPP) user is lower than the current uncertainties of the PPP solutions, after applying consistently the precise products (satellite orbits and clocks) obtained under I2+ correction

Consistency of seven different GNSS global ionospheric mapping techniques during one solar cycle

In the context of the International GNSS Service (IGS), several IGS Ionosphere Associated Analysis Centers have developed different techniques to provide global ionospheric maps (GIMs) of vertical total electron content (VTEC) since 1998. In this paper we present a comparison of the performances of all the GIMs created in the frame of IGS. Indeed we compare the classical ones (for the ionospheric analysis centers CODE, ESA/ESOC, JPL and UPC) with the new ones (NRCAN, CAS, WHU). To assess the quality of them in fair and completely independent ways, two assessment methods are used: a direct comparison to altimeter data (VTEC-altimeter) and to the difference of slant total electron content (STEC) observed in independent ground reference stations (dSTEC-GPS). The main conclusion of this study, performed during one solar cycle, is the consistency of the results between so many different GIM techniques and implementations.

A linear scale height Chapman model supported by GNSS occultation measurements

GNSS radio occultations allow the vertical sounding of the Earths atmosphere, in particular the ionosphere. The physical observables estimated with this technique permit to test theoretical models of the electron density such as, for example, the Chapman and the Vary-Chap models. The former is characterized by a constant scale height, whereas the latter considers a more general function of the scale height with respect to height. We propose to investigate the feasibility of the Vary-Chap model where the scale height varies linearly with respect to height. In order to test this hypothesis, the scale height data provided by radio occultations from a receiver on board a low Earth orbit (LEO) satellite, obtained by iterating with a local Chapman model at every point of the topside F2 layer provided by the GNSS satellite occultation, are fitted to height data by means of a linear least squares fit (LLS). Results, based on FORMOSAT-3/COSMIC GPS occultation data inverted by means of the Improved Abel transform inversion technique (which takes into account the horizontal electron content gradients), show that the scale height presents a more clear linear trend above the F2 layer peak height, hm, which is in good agreement with the expected linear temperature dependence. Moreover, the parameters of the linear fit obtained during four representative days for all seasons, depend significantly on local time and latitude, strongly suggesting that this approach can significantly contribute to build realistic models of the electron density directly derived from GNSS occultation data.

Impact and implementation of higher-order ionospheric effects on precise GNSS applications

High precision Global Navigation Satellite Systems (GNSS) positioning and time transfer require correcting signal delays, in particular higher-order ionospheric (I2+) terms. We present a consolidated model to correct second- and third-order terms, geometric bending and differential STEC bending effects in GNSS data. The model has been implemented in an online service correcting observations from submitted RINEX files for I2+ effects. We performed GNSS data processing with and without including I2+ corrections, in order to investigate the impact of I2+ corrections on GNSS products. We selected three time periods representing different ionospheric conditions. We used GPS and GLONASS observations from a global network and two regional networks in Poland and Brazil. We estimated satellite orbits, satellite clock corrections, Earth rotation parameters, troposphere delays, horizontal gradients, and receiver positions using a global GNSS solution, Real-Time Kinematic (RTK) and Precise Point Positioning (PPP) techniques. The satellite-related products captured most of the impact of I2+ corrections, with the magnitude up to 2 cm for clock corrections, 1 cm for the along- and cross-track orbit components, and below 5 mm for the radial component. The impact of I2+ on troposphere products turned out to be insignificant in general. I2+ corrections had limited influence on the performance of ambiguity resolution and the reliability of RTK positioning. Finally, we found that I2+ corrections caused a systematic shift in the coordinate domain that was time- and region-dependent, and reached up to -11 mm for the North component of the Brazilian stations during the most active ionospheric conditions. .

Electron density extrapolation above F2 peak by the linear Vary‐Chap model supporting new Global Navigation Satellite Systems‐LEO occultation missions

The new radio-occultation (RO) instrument on-board the future EUMETSAT Polar System 2nd Generation (EPS-SG) satellites, flying at a height of 820 km, is primarily focusing on neutral atmospheric profiling. It will also provide an opportunity for RO ionospheric sounding, but only below impact heights of 500 km, in order to guarantee a full data gathering of the neutral part. This will leave a gap of 320 km, which impedes the application of the direct inversion techniques to retrieve the electron density profile. To contribute to overcome this challenge, we have looked for new ways (accurate and simple) of extrapolating the electron density (also applicable to other Low-Earth-Orbiting, LEO, missions like CHAMP): a new Vary-Chap Extrapolation Technique (VCET). VCET is based on the scale height behaviour, linearly dependent on the altitude above hmF2 . This allows the electron density profile extrapolation for impact heights above its peak height (this is the case for EPS-SG), up to the satellite orbital height. VCET has been assessed with more than 3700 complete electron density profiles obtained in 4 representative scenarios of FORMOSAT-3/COSMIC occultations, in solar maximum and minimum conditions, and geomagnetically disturbed conditions, by applying an updated Improved Abel Transform Inversion technique to dual-frequency GPS measurements. It is shown that VCET performs much better than other classical Chapman models, with 60% of occultations showing relative extrapolation errors below 20%, in contrast with conventional Chapman model extrapolation approaches with 10% or less of the profiles with relative error below 20%.

Ionospheric Effects on GNSS Performance

This paper presents the features of the MONITOR project. This project initiated by ESA / ESTEC aims to increase the knowledge of the ionospheric effects and its impact on GNSS systems during active periods of solar activity. It includes the deployment of a set of GNSS-based ionospheric monitoring receivers worldwide distributed, the development of specific analysis software tools some of them integrated on a common platform, others distributed providing products routinely and a measurement campaign which will last beyond the peak of the current solar cycle.

The Barcelona ionospheric mapping function (BIMF) and its application to northern mid-latitudes

A simple way of improving the Global Navigation Satellite Systems (GNSS) slant ionospheric correction from Vertical Total Electron Content (VTEC) models is presented. In many GNSS applications, a mapping function is required to convert from VTEC, which may be provided in Global Ionospheric Maps (GIMs), to Slant TEC (STEC). Typical approaches assume a single ionospheric shell with constant height, which is unrealistic, especially for low-elevation signals. To reduce the associated conversion error, we propose the Barcelona Ionospheric Mapping Function and its first implementation at northern mid-latitudes (BIMF-nml). BIMF is based on a climatic prediction of the distribution of the topside vertical electron content fraction of VTEC (hereinafter µ2). BIMF is convenient to be applied since no external data are required in practice. To evaluate its performance, we use as independent reference the STEC difference (so-called dSTEC) values directly measured from mid-latitude dual-frequency Global Positioning System (GPS) receivers that have not been used in the computation of the VTEC GIMs under assessment. It is shown that the use of BIMF improves the GIM STEC estimation compared to the single-layer assumptions. This is the case for the mapping functions used by the International GNSS Service (IGS) and Satellite-Based Augmentation Systems (SBAS). This improvement is valid not only for the UPC GIMs, up to 15% for the year 2014, but especially for the GIMs of other analysis centers, such as those produced by CODE and JPL, up to 32 and 29%, respectively.

GPS as a solar observational instrument: Real-time estimation of EUV photons flux rate during strong, medium, and weak solar flares

In this manuscript, the authors show how the Global Navigation Satellite Systems, GNSS (exemplified in the Global Positioning System, GPS), can be efficiently used for a very different purpose from that for which it was designed as an accurate Solar observational tool, already operational from the open global GPS measurements available in real-time, and with some advantages regarding dedicated instruments onboard spacecraft. The very high correlation of the solar extreme ultraviolet (EUV) photon flux rate in the 26–34 mm spectral band, obtained from the solar EUV monitor instrument onboard the SOHO spacecraft during Solar flares, is shown with the GNSS solar flare activity indicator (GSFLAI). The GSFLAI is defined as the gradient of the ionospheric vertical total electron content rate versus the cosine of the Solar zenith angle in the day hemisphere (which filters out nonsolar over ionization), and it is measured from data collected by a global network of dual frequency GPS receivers (giving in this way continuous coverage). GSFLAI for 60 X class flares, 320 M class flares, and 300 C class flares, occurred since 2001, were directly compared with the EUV solar flux rate data to show existing correlations. It was found that the GSFLAI and EUV flux rate present the same linear relationship for all classes of flares, not only the strong and medium intensity ones, X and M class, as in previous works, but also for the weakest C class solar flares, which is a remarkable result.

Electron density extrapolation above F2 peak by the linear Vary-Chap model supporting new Global Navigation Satellite Systems-LEO occultation missions: ELECTRON DENSITY EXTRAPOLATION

The new radio-occultation (RO) instrument on board the future EUMETSAT Polar System-Second Generation (EPS-SG) satellites, flying at a height of 820 km, is primarily focusing on neutral atmospheric profiling. It will also provide an opportunity for RO ionospheric sounding, but only below impact heights of 500 km, in order to guarantee a full data gathering of the neutral part. This will leave a gap of 320 km, which impedes the application of the direct inversion techniques to retrieve the electron density profile. To overcome this challenge, we have looked for new ways (accurate and simple) of extrapolating the electron density (also applicable to other low-Earth orbiting, LEO, missions like CHAMP): a new Vary-Chap Extrapolation Technique (VCET). VCET is based on the scale height behavior, linearly dependent on the altitude above hmF2. This allows extrapolating the electron density profile for impact heights above its peak height (this is the case for EPS-SG), up to the satellite orbital height. VCET has been assessed with more than 3700 complete electron density profiles obtained in four representative scenarios of the Constellation Observing System for Meteorology, Ionosphere, and Climate (COSMIC) in the United States and the Formosa Satellite Mission 3 (FORMOSAT-3) in Taiwan, in solar maximum and minimum conditions, and geomagnetically disturbed conditions, by applying an updated Improved Abel Transform Inversion technique to dual-frequency GPS measurements. It is shown that VCET performs much better than other classical Chapman models, with 60% of occultations showing relative extrapolation errors below 20%, in contrast with conventional Chapman model extrapolation approaches with 10% or less of the profiles with relative error below 20%.