Use of varying shell heights derived from ionosonde data in calculating vertical total electron content (TEC) using GPS – New method

The dispersive nature of the ionosphere makes it possible to measure its total electron content (TEC). Thus Global Positioning System, which uses dual-frequency radio signals, is an ideal system to measure TEC. When data from an ionosonde situated in polar region was observed, the height of an approximated thin shell of electrons (shell height) used in GPS studies was seen not to be fixed but rather changing with time. Here we introduce a new method in which we included the varying shell heights derived from the ionosonde to map the slant total electron content from GPS to obtain a more precise vertical total electron content of the ionosphere contrary to some previous methods which used fixed shell heights. In this paper we also compared the ionosonde derived TEC with the GPS derived vertical TEC (vTEC) values. These GPS vTEC values were obtained from GPS slant TEC (sTEC) measurements using both fixed shell height and varying shell heights (from ionosonde measurements). For the polar regions, the varying shell height approach produced better results than the fixed shell height and compared to exponential function, Chapman function seems to be a better function to model the topside ionosphere.     

Evaluation of the IRI-2007 model options for the topside electron density

The International Reference Ionosphere (IRI) 2007 provides two new options for the topside electron density profile: (a) a correction of the IRI-2001 model, and (b) the NeQuick topside formula. We use the large volume of Alouette 1, 2 and ISIS 1, 2 topside sounder data to evaluate these two new options with special emphasis on the uppermost topside where IRI-2001 showed the largest discrepancies. We will also study the accurate representation of profiles in the equatorial anomaly region where the profile function has to accommodate two latitudinal maxima (crests) at lower altitudes but only a single maximum (at the equator) higher up. In addition to IRI-2001 and the two new IRI-2007 options we also include the Intercosmos-based topside model of Triskova, Truhlik, and Smilauer [Triskova, L., Truhlik, V., Smilauer, J. An empirical topside electron density model for calculation of absolute ion densities in IRI. Adv. Space Res. 37 (5), 928–934, 2006] (TTS model) in our analysis. We find that overall IRI-2007-NeQ gives the best results but IRI-2007-corrected provides a more realistic representation of the altitudinal–latitudinal structure in the equatorial anomaly region. The applicability of the TTS model is limited by the fact that it is not normalized to the F2 peak density and height.     

Symmetry and asymmetry of ionospheric weather at magnetic conjugate points for two midlatitude observatories

Variations of the ionospheric weather W-index for two midlatitude observatories, namely, Grahamstown and Hermanus, and their conjugate counterpart locations in Africa are studied for a period from October 2010 to December 2011. The observatories are located in the longitude sector, which has consistent magnetic equator and geographic equator so that geomagnetic latitudes of the line of force are very close to the corresponding geographic latitudes providing opportunity to ignore the impact of the difference of the gravitational field and the geomagnetic field at the conjugate points on the ionosphere structure and dynamics. The ionosondes of Grahamstown and Hermanus provide data of the critical frequency (foF2), and Global Ionospheric Maps (GIM) provide the total electron content (TECgps) along the magnetic field line up to the conjugate point in the opposite hemisphere. The global model of the ionosphere, International Reference Ionosphere, extended to the plasmasphere altitude of 20,200 km (IRI-Plas) is used to deliver the F2 layer peak parameters from TECgps at the magnetic conjugate area. The evidence is obtained that the electron gas heated by day and cooled by night at the summer hemisphere as compared with the opposite features in the conjugate winter hemisphere testifies on a reversal of plasma fluxes along the magnetic field line by the solar terminator. The ionospheric weather W-index is derived from NmF2 (related with foF2) and TECgps data. It is found that symmetry of W-index behavior in the magnetic conjugate hemispheres is dominant for the equinoxes when plasma movement along the magnetic line of force is imposed on symmetrical background electron density and electron content. Asymmetry of the ionospheric storm effects is observed for solstices when the plasma diffuse down more slowly into the colder winter hemisphere than into the warmer summer hemisphere inducing either plasma increase (positive phase) or decrease (negative phase of W-index) in the ionospheric and plasmaspheric plasma density.     

On the use of NeQuick topside option in IRI-2007

The latest version of IRI includes various options for the computation of the topside electron density profile. One of the possible choices is based on NeQuick model. Its inclusion in IRI has been made transferring all the formulations used in NeQuick model. In details, an Epstein layer function is used to describe the electron density profile and the topside shape is controlled by an empirical parameter, connected to the NeQuick F2 bottomside thickness parameter, B2bot. It is computed also in this IRI topside option in order to maintain self-consistency with its original formulation. This paper analyses the possibility of using the IRI bottomside parameters for this option and its impact on the profile and TEC. The case of experimental peak values given as input is also analysed.     

Evaluation of the IRI-2016 and NeQuick electron content specification by COSMIC GPS radio occultation,ground-based GPS and Jason-2 joint altimeter/GPS observations

We examined performance of two empirical profile-based ionospheric models, namely IRI-2016 and NeQuick-2, in electron content (EC) and total electron content (TEC) representation for different seasons and levels of solar activity. We derived and analyzed EC estimates in several representative altitudinal intervals for the ionosphere and the plasmasphere from the COSMIC GPS radio occultation, ground-based GPS and Jason-2 joint altimeter/GPS observations. It allows us to estimate a quantitative impact of the ionospheric electron density profiles formulation in several altitudinal intervals and to examine the source of the model-data discrepancies of the EC specification from the bottom-side ionosphere towards the GPS orbit altitudes. The most pronounced model-data differences were found at the low latitude region as related to the equatorial ionization anomaly appearance. Both the IRI-2016 and NeQuick-2 models tend to overestimate the daytime ionospheric EC and TEC at low latitudes during all seasons of low solar activity. On the contrary, during high solar activity the model results underestimated the EC/TEC observations at low latitudes. We found that both models underestimated the EC for the topside ionosphere and plasmasphere regions for all levels of solar activity. For low solar activity, the underestimated EC from the topside ionosphere and plasmasphere can compensate the overestimation of the ionospheric EC and, consequently, can slightly decrease the resulted model overestimation of the ground-based TEC. For high solar activity, the underestimated EC from the topside ionosphere and plasmasphere leads to a strengthening of the model underestimation of the ground-based TEC values. We demonstrated that the major source of the model-data discrepancies in the EC/TEC domain comes from the topside ionosphere/plasmasphere system.     

Reconstruction of topside density profile by using the topside sounder model profiler and digisonde data

To improve the accuracy of the real time topside electron density profiles given by the Digisonde software a new model-assisted technique is used. This technique uses the Topside Sounder Model (TSM), which provides the plasma scale height (H_s), O⁺–H⁺ transition height (H_T), and their ratio Rt = H_s/H_T, derived from topside sounder data of Alouette and ISIS satellites. The Topside Sounder Model Profiler (TSMP) incorporates TSM and uses the model quantities as anchor points in construction of topside density (Ne) profiles. For any particular location, TSMP calculates topside Ne profiles by specifying the values of foF2 and hmF2. In the present version, TSMP takes the F2 peak characteristics – foF2, hmF2, and the scale height at hmF2 – from the Digisonde measurements. The paper shows results for the Digisonde stations Athens and Juliusruh. It is found that the topside scale height used in Digisonde reconstruction is less than that extracted from topside sounder profiles. Rough comparison of their bulk distributions showed that they differ by an average factor of 1.25 for locations of Athens and Juliusruh. When the Digisonde scale heights are adjusted by this factor, the reconstructed topside profiles are close to those provided by TSM. Compared with CHAMP reconstruction profiles in two cases, TSMP/Digisonde profiles show lower density between 400 and 2000 km.     

Collocated ionosonde and dense GPS/GLONASS network measurements of midlatitude MSTIDs

To analyze midlatitude medium-scale travelling ionospheric disturbances (MSTIDs) over Kazan (55.5°N, 49°E), Russia, the sufficiently dense network of GNSS receivers (more than 150 ground-based stations) were used. For the first time, daytime MSTIDs in the form of their main signature (band structure) on high-resolution two-dimensional maps of the total electron content perturbation (TEC maps) are compared with ionosonde data with a high temporal resolution. For a pair of events, a relationship between southwestward TEC perturbations and evolution of F2 layer traces was established. So F2 peak frequency varied in antiphase to TEC perturbations. The ionograms show that during the movement of plasma depletion band (overhead ionosonde) the F2 peak frequency is the highest, and vice versa, for the plasma enhancement band, the F2 peak frequency is the lowest. One possible explanation may be a greater inclination of the radio beam from the vertical during the placement of a plasma enhancement band above the ionosonde, as evidenced by the absence of multiple reflections and the increased occurrence rate of additional cusp trace. Another possible explanation may be the redistribution of the electron content in the topside ionosphere with a small decrease in the F peak concentration of the layer with a small increase in TEC along the line-of-sight. Analysis of F2 peak frequency variation shows that observed peak-to-peak values of TEC perturbation equal to 0.4 and 1 TECU correspond to the values of ΔN/N equal to 13% and 28%. The need for further research is evident.     

Four-peak longitudinal distribution of the equatorial plasma bubbles observed in the topside ionosphere: Possible troposphere tide influence

In this paper we consider an idea of the troposphere tide influence on the character of the longitudinal variations in the distribution of the equatorial plasma bubbles (EPBs) observed in the topside ionosphere. For this purpose, the obtained EPB longitudinal patterns were compared with the thermosphere and ionosphere characteristics having the prominent “wave-like” longitudinal structures with wave number 4, which are uniquely associated with the influence of the troposphere DE3 tides. The characteristics of the equatorial mass density anomaly (EMA), equatorial ionization anomaly (EIA), zonal wind and pre-reversal E?×?B drift enhancement (PRE) were used for comparison. The equinox seasons during high solar activity were under consideration. It was obtained that the longitudinal patterns of the EMA and zonal wind show the surprising similarity with the EPB distributions (R???0.8, R???0.72). On the other hand, the resemblance with the ionosphere characteristics (EIA, PRE) is rather faint (R???0.37, R???0.12). It was shown that the thermosphere zonal winds are the most possible transfer mediator of the troposphere DE3 tide influence. The most successful moment for the transfer of the troposphere DE3 tide energy takes place in the beginning of the EPB production, namely, during the seed perturbation development.     

Neural network prediction of the topside electron content over the Euro-African sector derived from Swarm-A measurements

This study presents the first prediction results of a neural network model for the vertical total electron content of the topside ionosphere based on Swarm-A measurements. The model was trained on 5 years of Swarm-A data over the Euro-African sector spanning the period 1 January 2014 to 31 December 2018. The Swarm-A data was combined with solar and geomagnetic indices to train the NN model. The Swarm-A data of 1 January to 30 September 2019 was used to test the performance of the neural network. The data was divided into two main categories: most quiet and most disturbed days of each month. Each category was subdivided into two sub-categories according to the Swarm-A trajectory i.e. whether it was ascending or descending in order to accommodate the change in local time when the satellite traverses the poles. Four pairs of neural network models were implemented, the first of each pair having one hidden layer, and the second of each pair having two hidden layers, for the following cases: 1) quiet day-ascending, 2) quiet day-descending, 3) disturbed day-ascending, and 4) disturbed day-descending. The topside vertical total electron content predicted by the neural network models compared well with the measurements by Swarm-A. The model that performed best was the one hidden layer model in the case of quiet days for descending trajectories, with RMSE = 1.20 TECU, R = 0.76. The worst performance occurred during the disturbed descending trajectories where the one hidden layer model had the worst RMSE = 2.12 TECU, (R = 0.54), and the two hidden layer model had the worst correlation coefficient R = 0.47 (RMSE = 1.57).In all cases, the neural network models performed better than the IRI2016 model in predicting the topside total electron content. The NN models presented here is the first such attempt at comparing NN models for the topside VTEC based on Swarm-A measurements.     

A new ionospheric tomography model combining pixel-based and function-based models

Considering the limitations of current single pixel-based and function-based computerized ionospheric tomography (CIT) models, this paper proposes a new tomography model – COMBI, which combines these two models. COMBI model is able to reconstruct the three dimensional distribution of electron density with fewer parameters, and easy to compute, as well as very convenient to use. Through experiments with simulated data and measured data, it is verified that the new COMBI model not only can better describe refine structure of ionospheric electron density, but also is superior to these two pixel-based and function-based CIT models in application.     

Linkage of the ionospheric peak electron density and height deduced from the topside sounding data

Topside sounding electron density profiles are analyzed to explore interrelations of the F2 layer critical frequency and the peak height for a representative set of conditions provided by ISIS1, ISIS2, IK19 and Cosmos-1809 satellites for the period of 1969–1987. The foF2 and hmF2 are delivered with exponential extrapolation of electron density profile to zero of its 1st derivative. It is shown that the linear regression exists between foF2 and hmF2 under different conditions. The linkage between the two parameters amended to the empirical model of the peak height [Gulyaeva, T.L., Bradley, P.A., Stanislawska, I., Juchnikowski, G. Towards a new reference model of hmF2 for IRI. Adv. Space Res. 42, 666–672, doi:10.1016/j.asr.2008.02.021, 2008] results in an empirical model of the both foF2 and hmF2 expressed by superposition of functions in terms of local-time, season, geodetic longitude, modified dip latitude and solar activity. For the solar activity we use a proxy Fsp index averaged from the mean solar radio flux F10.7s for the past 81 days (3 solar rotations) and F10.7 value for 1 day prior the day of observation. Impact of geomagnetic activity is not discernible with the topside sounding data due to mixed positive and negative storm-time effects. Appreciable differences have been revealed between IRI-CCIR predictions and outcome of the new model which might be attributed to the different techniques of the peak electron density and height derivation, different epochs and different global distribution of the source data as well as the different mathematical functions involved in the maps and the model presentation.     

Propagation of plasma bubbles observed in Brazil from GPS and airglow data

Equatorial spread-F is a common occurrence in the equatorial ionosphere that is associated with large variations in plasma density that often cause scintillation and interference in communication signals. These events are known to result from Rayleigh–Taylor instability, but the day-to-day variability of their occurrence is not well understood. The triggering mechanism of plasma depletions is still a matter of debate, but may be linked to gravity waves that under favorable conditions propagate to the middle atmosphere. Understanding the triggering of ESF was the focus of the SpreadFEx campaign near Brasilia, Brazil in 2005. The campaign provided co-located airglow and GPS observations to study the onset of plasma depletions and their evolution as they traversed the region. Comparisons between the 630.0 nm airglow data and GPS data demonstrate the ability of the compact dual frequency GPS array to detect the plasma bubbles and retrieve reliable propagation characteristics of the depletions. In this case study, a plasma depletion was detected and moved over the array at velocities of 85–110 m/s, slowing as it moved towards the east. Correlation of consecutive airglow images gives consistent estimates of the eastward drift over the same time period. Mapping the airglow data to the GPS line-of-sight geometry allows direct comparison and reveals a resolvable westward tilt of the plasma depletion that may be due to vertical shear. The uniqueness of this study is the ability to resolve locally the characteristics of the plasma depletion without relying on assumptions about the mapping of the depletion along magnetic field lines to large latitudinal distances. It presents new information for understanding ESF development and the development of depletions strong enough to produce scintillation.     

TEC derived from some GPS stations in Nigeria and comparison with the IRI and NeQuick models

Total electron content (TEC) measured simultaneously using Global Positioning System (GPS) ionospheric monitors installed at some locations in Nigeria during the year 2011 (Rz = 55.7) was used to study the diurnal, seasonal, and annual TEC variations. The TEC exhibits daytime maximum, seasonal variation and semiannual variations. Measured TEC were compared with those predicted by the improved versions of the International Reference Ionosphere (IRI) and NeQuick models. The models followed the diurnal and seasonal variation patterns of the observed values of TEC. However, IRI model produced better estimates of TEC than NeQuick at all locations.     

Comparison of model and experimental ionospheric parameters at high latitudes

Modern use and study of the auroral region needs to attract a wider class of models for describing conditions of radio wave propagation in the ionosphere. In this paper the possibilities of the International Reference Ionosphere model, well-proven and widespread in the mid-latitudes, are investigated in the high latitude zone. Model and measured values of the critical frequency foF2 for two mid-latitude stations (Juliusruh and Goosebay) and four high-latitude ones (Loparskaya, Sodankyla, Sondrestrom, Thule) are compared. Deviations of medians, variations from day to day and solar activity trends seemed to be similar for both areas. This similarity is irrespective of the RZ12 index. Special attention is paid to the TEC parameter and its determination using 6 versions of models, a new version of the model IRI2010 (IRI-Plas) among them. It is shown that the IRI-Plas model significantly improves the definition of TEC in contrast to the versions of IRI2007 and the new model NeQuick. The use of the median of the experimental equivalent slab thickness, together with the current values of the TEC, increases by a factor of two the agreement between calculated and measured foF2 values as compared with the variations from day to day. This allows foF2 to be defined in near-real time.     

Mapping seasonal trends of electron temperature in the topside ionosphere based on DEMETER data

The diurnal, seasonal and latitudinal variations of the electron temperature in the Earth‘s topside ionosphere during relatively low solar activity period of 2005 – 2008 are investigated. In order to examine seasonal variations and morphology of the topside ionospheric plasma temperature, CNES micro-satellite DEMETER ISL data are used. Presented study is oriented on the dataset gathered in 2005 and 2008. Within conducted analysis, global maps of electron temperature for months of equinoxes and solstices have been developed. Furthermore, simultaneous studies on two-dimensional time series based on DEMETER measurements and predictions obtained with the IRI-2012 model supply examination of the topside ionosphere during recent deep solar minimum. Comparison with the IRI-2012 model reveals discrepancies between data and prediction, that are especially prominent during the periods of very low solar activity.     

South American regional ionospheric maps computed by GESA: A pilot service in the framework of SIRGAS

SIRGAS (Geocentric Reference Frame for the Americas) is an international enterprise of the geodetic community that aims to realize the Terrestrial Reference Frame in the America’s countries. In order to fulfill this commitment, SIRGAS manages a network of continuously operational GNSS receivers totalling around one hundred sites in the Caribbean, Central, and South American region. Although the network was not planed for ionospheric studies, its potential to be used for such a purpose was recently recognized and SIRGAS started a pilot experiment devoted to establish a regular service for computing and releasing regional vertical TEC (vTEC) maps based on GNSS data. Since July, 2005, the GESA (Geodesia Espacial y Aeronomía) laboratory belonging to the Facultad de Ciencias Astronómicas y Geofísicas of the Universidad Nacional de La Plata computes hourly maps of vertical Total Electron Content (vTEC) in the framework of the SIRGAS pilot experiment. These maps exploit all the GNSS data available in the South American region and are computed with the LPIM (La Plata Ionospheric Model). LPIM implements a de-biasing procedure that improves data calibration in relation to other procedures commonly used for such purposes. After calibration, slant TEC measurements are converted to vertical and mapped using local-time and modip latitude. The use of modip latitude smoothed the spatial variability of vTEC, especially in the South American low latitude region and hence allows for a better vTEC interpolation. This contribution summarizes the results obtained by GESA in the framework of the SIRGAS pilot experiment.     

Medium-scale traveling ionospheric disturbances (MSTID) modeling using a dense German GPS network

The MSTIDs are wave-like perturbations of the ionospheric plasma, which cause the most common ionospheric disturbances in mid-latitude regions. Generally the MSTIDs have velocities of several hundred meters per second and wavelengths of several hundred kilometers. The wave-like effect of the MSTID is one of the main obstacles for accurate interpolation of ionospheric corrections in a medium-scale reference GPS network. In this paper we show a new method of detecting and modeling MSTIDs using dense German GPS network. The between-epoch single difference ionospheric delays from a medium scale dense GPS network are used to estimate the parameter of the MSTID e.g. amplitude, wavelength and velocity. The efficiency of the approach is tested with data from about 320 GPS stations in and near Germany. A MSTID wave moving from east to west across Germany was observed at September 27 in 2009. Its wavelength is about 302 km, with a period of ∼7 min and velocity of about 700 m/s.     

Effects of X2-class solar flare events on ionospheric GPS-TEC and radio waves over Brazilian sector

In this investigation, we present and discuss the effects of 6 X2-class solar flare events in the ionospheric F region over Brazilian sector that occurred during 2013 to 2015. For this investigation, we present the vertical total electron content (VTEC) observations from nearly 120 Global Positioning System (GPS) receivers all over the Brazilian sector for each event. Also, ionospheric sounding observations obtained in São José dos Campos (23.2°S, 45.9°W, dip latitude 17.6°S; hereafter referred to as SJC), under the southern crest of the equatorial ionospheric anomaly (EIA), Brazil, are presented. The observations show that the greatest TEC impact occurs with the EUV fluxes increases lasting for more than one hour and when the solar active region is located close to the solar disc center. We present a detailed study of the efficiency of the EUV flux with wavelengths ranging from 0.1 to 190?nm for the F region ionization. The largest increase of ΔTEC occurs below the magnetic equator line, covering mainly the central, northeast, southeast and south regions, which includes the equatorial ionospheric anomaly (EIA) region. The ionograms show partial or total fade out in the echoes traces observed causing blackouts of radio signals of up to 60?min, which can have serious consequences to technological systems of public and private agencies around Brazilian sector. This study can help to better understand the effects of solar flares in the ionospheric F region.     

An analysis of the topside ionosphere parameters based on the long-duration Irkutsk incoherent scatter radar measurements

The topside ionosphere parameters are studied based on the long-duration Irkutsk incoherent scatter radar (52.9N, 103.3E) measurements conducted in September 2005, June and December 2007. As a topside ionosphere parameter we chose the vertical scale height (VSH) related to the gradient of the electron density logarithm above the peak height. For morphological studies we used median electron density profiles. Besides the median behavior we also studied VSH disturbances (deviations from median values) during the magnetic storm of September 11th 2005. We compared the Irkutsk incoherent scatter radar data with the Millstone Hill and Arecibo incoherent scatter radar observations, the IRI-2007 prediction (using the two topside options) and VSH derived from the Irkutsk DPS-4 Digisonde bottomside measurements.