Testing the three-dimensional IRI-SIRMUP-P mapping of the ionosphere for disturbed periods

This paper describes the three-dimensional (3-D) electron density mapping of the ionosphere given as output by the assimilative IRI-SIRMUP-P (ISP) model for three different geomagnetic storms. Results of the 3-D model are shown by comparing the electron density profiles given by the model with the ones measured at two testing ionospheric stations: Roquetes (40.8°N, 0.5°E), Spain, and San Vito (40.6°N, 17.8°E), Italy. The reference ionospheric stations from which the autoscaled foF2 and M(3000)F2 data as well as the real-time vertical electron density profiles are assimilated by the ISP model are those of El Arenosillo (37.1°N, 353.3°E), Spain, Rome (41.8°N, 12.5°E), and Gibilmanna (37.9°N, 14.0°E), Italy. Overall, the representation of the ionosphere made by the ISP model is better than the climatological representation made by only the IRI-URSI and the IRI-CCIR models. However, there are few cases for which the assimilation of the autoscaled data from the reference stations causes either a strong underestimation or a strong overestimation of the real conditions of the ionosphere, which is in these cases better represented by only the IRI-URSI model. This ISP misrepresentation is mainly due to the fact that the reference ionospheric stations covering the region mapped by the model turn out to be few, especially for disturbed periods when the ionosphere is very variable both in time and in space and hence a larger number of stations would be required. The inclusion of new additional reference ionospheric stations could surely smooth out this concern.     

An automatic quality factor for Autoscala foF2 values

We propose a new parameter for quality evaluation of ionogram traces reconstructed by Autoscala. This parameter efficiently assesses the reliability of the automatic interpretation of ionospheric characteristics. Based on an extensive analysis of the data, the parameter values are statistically associated with the accuracy of f_oF2 data automatically scaled by Autoscala. Therefore, Autoscala will be improved by providing f_oF2 accuracy as supplementary output information.     

Application of Autoscala to ionograms recorded by the VIPIR ionosonde

In November 2008, the ionosonde station at Boulder, Colorado, USA (40.0°N; 105.3°W) became the host of a new ionosonde (VIPIR, Vertical Incidence Pulsed Ionospheric Radar) developed and built by Scion Associates.     

Multi-beam sounding ionograms in the polar cap region: Absorption induced by proton precipitations

The simulation of the multi-beam ionograms in the polar cap region, assessing absorption effect is performed. It is reasonable to distinguish among four different mechanisms responsible for absorption: regular absorption due to solar UV illumination, absorption associated with energetic particles precipitation, absorption connected with X-rays flare and additional absorption in Auroral oval area. In this paper the absorption attributed to proton precipitations is envisaged. The computational model of the high-latitude ionosphere with irregularities oriented to application for the high frequency wave propagation problem was elaborated (Zaalov et al., 2005). A number of the quasi-vertical ionograms in the polar cap region were simulated on the basis of this model. A well-known algorithm (Sauer and Wilkinson, 2008) is applied for the absorption effects calculation. The simulated high-latitude ionograms with the absorption effect and the measured ionograms exhibit quite a good resemblance. This paper illustrates the importance of the understanding and taking into account the absorption effect in the presence of the various structural features in the polar ionosphere (in particular, patches of enhanced electron density) in interpreting ionosonde data.     

A new inversion algorithm for HF sky-wave backscatter ionograms

HF sky-wave backscatter sounding system is capable of measuring the large-scale, two-dimensional (2-D) distributions of ionospheric electron density. The leading edge (LE) of a backscatter ionogram (BSI) is widely used for ionospheric inversion since it is hardly affected by any factors other than ionospheric electron density. Traditional BSI inversion methods have failed to distinguish LEs associated with different ionospheric layers, and simply utilize the minimum group path of each operating frequency, which generally corresponds to the LE associated with the F₂ layer. Consequently, while the inversion results can provide accurate profiles of the F region below the F₂ peak, the diagnostics may not be so effective for other ionospheric layers. In order to resolve this issue, we present a new BSI inversion method using LEs associated with different layers, which can further improve the accuracy of electron density distribution, especially the profile of the ionospheric layers below the F₂ region. The efficiency of the algorithm is evaluated by computing the mean and the standard deviation of the differences between inverted parameter values and true values obtained from both vertical and oblique incidence sounding. Test results clearly manifest that the method we have developed outputs more accurate electron density profiles due to improvements to acquire the profiles of the layers below the F₂ region. Our study can further improve the current BSI inversion methods on the reconstruction of 2-D electron density distribution in a vertical plane aligned with the direction of sounding.     

Electron density profile calculation technique for Autoscala ionogram analysis

An electron density profile model with free parameters is introduced. Initially the parameters are calculated on the basis of the ionospheric characteristics automatically obtained from the ionograms by Autoscala and considering the helio-geophysical conditions. The technique used to adjust the free parameters to the particular ionograms recorded is presented.     

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.     

Histogram-based ionogram displays and their application to autoscaling

A simple method is described for displaying and auto scaling the basic ionogram parameters foF2 and h’F2 as well as some additional layer parameters from digital ionograms. The technique employed is based on forming frequency and height histograms in each ionogram. This technique has now been applied specifically to ionograms produced by the IPS5D ionosonde developed and operated by the Australian Space Weather Service (SWS). The SWS ionograms are archived in a cleaned format and readily available from the SWS internet site. However, the method is applicable to any ionosonde which produces ionograms in a digital format at a useful signal-to-noise level. The most novel feature of the technique for autoscaling is its simplicity and the avoidance of the mathematical imaging and line fitting techniques often used. The program arose from the necessity to display many days of ionogram output to allow the location of specific types of ionospheric event such as ionospheric storms, travelling ionospheric disturbances and repetitive ionospheric height changes for further investigation and measurement. Examples and applications of the method are given including the removal of sporadic E and spread F.     

The effect of collisions in ionogram inversion

The results of this paper demonstrate that the effect of collisions on the group refraction index is small, when the ordinary ray is considered. If, however, in order to improve the performance of a system for automatic interpretation of ionograms, the information contained in ordinary and extraordinary traces is combined, the effect of collisions between the electrons and neutral molecules should be taken into account for the extraordinary ray. The magnitude of these differences is generally very small and must be compared with the resolution in the virtual vertical height of the ionosonde, resolution which is typically of the order of few kilometers.     

Validation of University of New Brunswick Ionospheric Modeling Technique with ionosonde TEC estimation over South Africa

For more than a decade, ionospheric research over South Africa has been carried out using data from ionosondes geographically located at Madimbo (28.38°S, 30.88°E), Grahamstown (33.32°S, 26.50°E), and Louisvale (28.51°S, 21.24°E). The objective has been modelling the bottomside ionospheric characteristics using neural networks. The use of Global Navigation Satellite System (GNSS) data is described as a new technique to monitor the dynamics and variations of the ionosphere over South Africa, with possible future application in high frequency radio communication. For this task, the University of New Brunswick Ionospheric Modelling Technique (UNB-IMT) was applied to compute midday (10:00 UT) GNSS-derived total electron content (GTEC). GTEC values were computed using GNSS data for stations located near ionosondes for the years 2002 and 2005 near solar maximum and minimum, respectively. The GTEC was compared with the midday ionosonde-derived TEC (ITEC) measurements to validate the UNB-IMT results. It was found that the variation trends of GTEC and ITEC over all stations are in good agreement and show a pronounced seasonal variation for the period near solar maximum, with maximum values (∼80 TECU) around autumn and spring equinoxes, and minimum values (∼22 TECU) around winter and summer. Furthermore, the residual ΔTEC = GTEC − ITEC was computed. It was evident that ΔTEC, which is believed to correspond to plasmaspheric electron content, showed a pronounced seasonal variation with maximum values (∼20 TECU) around equinoxes and minimum (∼5 TECU) around winter near solar maximum. The equivalent ionospheric and total slab thicknesses were also computed and comprehensively discussed. The results verified the use of UNB-IMT as one of the tools for future ionospheric TEC research over South Africa.