Klobuchar, NeQuickG, BDGIM, GLONASS, IRI-2016, IRI-2012, IRI-Plas, NeQuick2, and GEMTEC Ionospheric Models: A Comparison in Total Electron Content and Positioning Domains.

Global navigation satellite systems (GNSS) provide a great data source about the ionosphere state. These data can be used for testing ionosphere models. We studied the performance of nine ionospheric models (Klobuchar, NeQuickG, BDGIM, GLONASS, IRI-2016, IRI-2012, IRI-Plas, NeQuick2, and GEMTEC) both in the total electron content (TEC) domain-i.e., how precise the models calculate TEC-and in the positioning error domain-i.e., how the models improve single frequency positioning. The whole data set covers 20 years (2000-2020) from 13 GNSS stations, but the main analysis involves data during 2014-2020 when calculations are available from all the models. We used single-frequency positioning without ionospheric correction and with correction via global ionospheric maps (IGSG) data as expected limits for errors. Improvements against noncorrected solution were as follows: GIM IGSG-22.0%, BDGIM-15.3%, NeQuick2-13.8%, GEMTEC, NeQuickG and IRI-2016-13.3%, Klobuchar-13.2%, IRI-2012-11.6%, IRI-Plas-8.0%, GLONASS-7.3%. TEC bias and mean absolute TEC errors for the models are as follows: GEMTEC--0.3 and 2.4 TECU, BDGIM--0.7 and 2.9 TECU, NeQuick2--1.2 and 3.5 TECU, IRI-2012--1.5 and 3.2 TECU, NeQuickG--1.5 and 3.5 TECU, IRI-2016--1.8 and 3.2 TECU, Klobuchar-1.2 and 4.9 TECU, GLONASS--1.9 and 4.8 TECU, and IRI-Plas-3.1 and 4.2 TECU. While TEC and positioning domains differ, new-generation operational models (BDGIM and NeQuickG) could overperform or at least be at the same level as classical empirical models.

B2 Thickness Parameter Response to Equinoctial Geomagnetic Storms.

The thickness parameters that most empirical models use are generally defined by empirical relations related to ionogram characteristics. This is the case with the NeQuick model that uses an inflection point below the F2 layer peak to define a thickness parameter of the F2 bottomside of the electron density profile, which is named B2. This study is focused on the effects of geomagnetic storms on the thickness parameter B2. We selected three equinoctial storms, namely 17 March 2013, 2 October 2013 and 17 March 2015. To investigate the behavior of the B2 parameter before, during and after those events, we have analyzed variations of GNSS derived vertical TEC (VTEC) and maximum electron density (NmF2) obtained from manually scaled ionograms over 20 stations at middle and low latitudes of Asian, Euro-African and American longitude sectors. The results show two main kinds of responses after the onset of the geomagnetic events: a peak of B2 parameter prior to the increase in VTEC and NmF2 (in ~60% of the cases) and a fluctuation in B2 associated with a decrease in VTEC and NmF2 (~25% of the cases). The behavior observed has been related to the dominant factor acting after the CME shocks associated with positive and negative storm effects. Investigation into the time delay of the different measurements according to location showed that B2 reacts before NmF2 and VTEC after the onset of the storms in all the cases. The sensitivity shown by B2 during the studied storms might indicate that experimentally derived thickness parameter B2 could be incorporated into the empirical models such as NeQuick in order to adapt them to storm situations that represent extreme cases of ionospheric weather-like conditions.

On the link between the topside ionospheric effective scale height and the plasma ambipolar diffusion, theory and preliminary results.

Over the years, an amount of models relying on effective parameters were implemented in the challenging issue of the topside ionosphere description. These models are based on different analytical functions, but all of them depend on a parameter called effective scale height, that is deduced from topside electron density measurements. As their names state, they are effective in reproducing the topside electron density profile only when applied to the analytical function used to derive them. Then, in principle, they do not have any physical meaning. It is the goal of this paper to mathematically link the effective scale height modeled through the Epstein layer to the vertical scale height theoretically deduced from the plasma ambipolar diffusion theory. Firstly, effective and theoretical scale heights are linked through a mathematical relation by showing that they tend to each other in the topside ionosphere. Secondly, their connection is preliminarily demonstrated by calculating effective scale height values from the entire COSMIC/FORMOSAT-3 radio occultation dataset. Thirdly, a possible connection between the vertical gradient of the topside scale height (as obtained by COSMIC/FORMOSAT-3 satellites) and the electron temperature (as obtained by ESA Swarm B satellite) is studied by highlighting corresponding similarities in the diurnal, seasonal, solar activity, and latitudinal variability.