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.     

A comparative study of ionospheric IRIEup and ISP assimilative models during some intense and severe geomagnetic storms

Three-dimensional (3-D) electron density matrices, computed in the Mediterranean area by the IRI climatological model and IRIEup and ISP nowcasting models, during some intense and severe geomagnetic-ionospheric storms, were ingested by the ray tracing software tool IONORT, to synthesize quasi-vertical ionograms. IRIEup model was run in different operational modes: (1) assimilating validated autoscaled electron density profiles only from a limited area which, in our case, is the Mediterranean sector (IRIEup_re(V) mode); (2) assimilating electron density profiles from a larger region including several stations spread across Europe: (a) without taking care of validating the autoscaled data in the assimilation process (IRIEup(NV)); (b) validating carefully the autoscaled electron density profiles before their assimilation (IRIEup(V)).The comparative analysis was carried out comparing IRI, IRIEup_re(V), ISP, IRIEup(NV), and IRIEup(V) foF2 synthesized values, with corresponding foF2 measurements autoscaled by ARTIST, and then validated, at the truth sites of Roquetes (40.80°N, 0.50°E, Spain), San Vito (40.60°N, 17.80°E, Italy), Athens (38.00°N, 23.50°E, Greece), and Nicosia, (35.03°N, 33.16°E, Cyprus). The outcomes demonstrate that: (1) IRIEup_re(V), performs better than ISP in the western Mediterranean (around Roquetes); (2) ISP performs slightly better than IRIEup_re(V) in the central part of Mediterranean (around Athens and San Vito); (3) ISP performance is better than the IRIEup_re(V) one in the eastern Mediterranean (around Nicosia); (4) IRIEup(NV) performance is worse than the IRIEup(V) one; (5) in the central Mediterranean area, IRIEup(V) performance is better than the IRIEup_re(V) one, and it is practically the same for the western and eastern sectors.Concerning the overall performance, nowcasting models proved to be considerably more reliable than the climatological IRI model to represent the ionosphere behaviour during geomagnetic-ionospheric storm conditions; ISP and IRIEup(V) provided the best performance, but neither of them has clearly prevailed over the other one.     

A comparative sporadic-E layer study between two mid-latitude ionospheric stations

Hourly systematic measurements of the highest frequency reflected by the sporadic-E layer (foEs) recorded from January 1976 to June 2009 at the ionospheric stations of Rome (Italy, 41.8°N, 12.5°E) and Gibilmanna (Italy, 37.9°N, 14.0°E) were considered to carry out a comparative study between the sporadic E layer (Es) over Rome and Gibilmanna. Different statistical analysis were performed taking into account foEs observations near the periods of minimum and maximum solar activity. The results reveal that: (1) independently from the solar activity, Es develops concurrently over extended regions in space, instead of being a spatially limited layer which is transported horizontally by neutral winds over a larger area; especially during summer months, when an Es layer is present at Rome, there is a high probability that an Es layer is also present over Gibilmanna, and vice versa; (2) Es layer lifetimes of 1–5 h were found; in particular, Es layers with lifetimes of 5 h both over Gibilmanna and Rome are observed with highest percentages of occurrence in summer ranging between 80% and 90%, independently from the solar activity; (3) latitudinal effects on Es layer occurrence emerge mostly for low solar activity during winter, equinoctial, and summer months, when Es layers are detected more frequently over Gibilmanna rather than Rome; (4) when the presence of an Es layer over Rome and Gibilmanna is not simultaneous, Es layer appearance both over Rome and Gibilmanna confirms to be a locally confined event, because drifting phenomena from Rome to Gibilmanna or vice versa have not been emphasized.     

Power variation analysis of echo signals from ionospheric reflectors

A series of Power Virtual Height measurements (PVH) of radio echoes reflected from the ionosphere were acquired at a given frequency during the period 3–22 January 2008 with the purpose of studying the slow fading variations through time of the ionospheric channel. To obtain PVH data, an ionospheric vertical sounding system was suitably adapted to work at a single fixed frequency. PVH measurements were recorded between two routine ionospheric vertical soundings, providing a data type that enables evaluation of fading fluctuation through time. The time stability of the ionospheric layers is determined by analyzing the level of the received signal power within a chosen threshold. In this paper the fading behaviour and its characteristics are described, considering only temporal periodicity above 0.5 s. In a further analysis a relation is demonstrated between the recorded fading and the time stability of the signal within a fixed interval of values.     

Fading in the HF ionospheric channel and the role of irregularities

It is well known that the ionosphere affects radio wave propagation especially in the high frequency (HF) range. HF radio waves reflected by the ionosphere can reach considerable distances, often with changes in amplitude, phase, and frequency. The ionosphere is a dispersive in frequency and time, bi-refractive, absorbing medium, in which multipath propagation due to traveling irregularities is very frequent. The traveling irregularities undulate the reflecting ionospheric layer, introducing variations in signal amplitude (fading). In this multipath time variant channel fading is mainly considered, even though it is not the sole effect. Echo signals from a single reflection, as in ionospheric vertical sounding (VIS) techniques, are affected by a certain degree of variability even in quiet ionospheric conditions. In this work the behavior of the ionospheric channel is studied and characterized by observing the power variation of received echoes using the VIS technique. Multipath fading was analyzed quantifying the power variation of the signal echo due to irregularities on a temporal scale from 0.5 to 256 s. An experimental set-up derived from an ionosonde was implemented and the analysis was performed employing a special numerical algorithm operating off-line on the acquired time sequence of the signal. The gain-loss of the irregularity shapes are determined in some special cases.     

Short-term forecasting regional model to predict M(3000)F2 over the European sector: Comparisons with the IRI model during moderate,disturbed, and very disturbed geomagnetic conditions

The hourly measurements of M(3000)F2 (M(3000)F2_meas) and the hourly quiet-time values of M(3000)F2 (M(3000)F2_QT) relative to the ionospheric observatories of Poitiers, Lannion, Dourbes, Slough, Rome, Juliusruh, Kaliningrad, Uppsala, Lyckesele, Sodankyla, and Kiruna as well as the hourly time-weighted accumulation series derived from the geomagnetic planetary index ap (ap(τ)), were considered during the period January 1957–December 2003 and used for the development of 11 short-term forecasting local models (STFLM) of M(3000)F2.     

Time of flight measurements over a radio link from Uppsala to Bruntingthorpe and their application to testing predictions methods that approximate the ray tracing technique

Time of flight measurements (TOF) over the radio link between Uppsala (Tx: 59.9°N, 17.6°E) and Bruntingthorpe (Rx: 52.5°N, 1.1°W) have been performed every 2 min at six frequencies (4.637, 6.954, 8.008, 10.391, 11.118, and 14.364 MHz) during the period November 2006–January 2008. Such measurements have been compared with the TOF provided by three prediction methods that approximate the ray tracing technique: IRI-95, SIRM&BR_D, and ICEPAC. The root mean square deviation (rms) between TOF monthly median measurements and TOF monthly median predictions and the differences (DP) between the length of the median and predicted ray path have been calculated. The results, which are presented in terms of rms and DP for different seasons and different time periods, have indicated that the approximate methods are inadequate and that for more accurate predictions ray tracing techniques should be applied.     

The COMPLEIK subroutine of the IONORT-ISP system for calculating the non-deviative absorption: A comparison with the ICEPAC formula

The present paper proposes to discuss the ionospheric absorption, assuming a quasi-flat layered ionospheric medium, with small horizontal gradients. A recent complex eikonal model (Settimi et al., 2013b) is applied, useful to calculate the absorption due to the ionospheric D-layer, which can be approximately characterized by a linearized analytical profile of complex refractive index, covering a short range of heights between h₁ = 50 km and h₂ = 90 km. Moreover, Settimi et al. (2013c) have already compared the complex eikonal model for the D-layer with the analytical Chapman’s profile of ionospheric electron density; the corresponding absorption coefficient is more accurate than Rawer’s theory (1976) in the range of middle critical frequencies. Finally, in this paper, the simple complex eikonal equations, in quasi-longitudinal (QL) approximation, for calculating the non-deviative absorption coefficient due to the propagation across the D-layer are encoded into a so called COMPLEIK (COMPLex EIKonal) subroutine of the IONORT (IONOspheric Ray-Tracing) program ( Azzarone et al., 2012). The IONORT program, which simulates the three-dimensional (3-D) ray-tracing for high frequencies (HF) waves in the ionosphere, runs on the assimilative ISP (IRI-SIRMUP-P) discrete model over the Mediterranean area ( Pezzopane et al., 2011). As main outcome of the paper, the simple COMPLEIK algorithm is compared to the more elaborate semi-empirical ICEPAC formula (Stewart, undated), which refers to various phenomenological parameters such as the critical frequency of E-layer. COMPLEIK is reliable just like the ICEPAC, with the advantage of being implemented more directly. Indeed, the complex eikonal model depends just on some parameters of the electron density profile, which are numerically calculable, such as the maximum height.     

Oblique-incidence ionospheric soundings over Central Europe and their application for testing now casting and long term prediction models

After a first oblique-incidence ionospheric sounding campaign over Central Europe performed during the period 2003–2004 over the radio links between Inskip (UK, 53.5°N, 2.5°W) and Rome (Italy, 41.8°N, 12.5°E) and between Inskip and Chania (Crete, 35.7°N, 24.0°E), new and more extensive analysis of systematic MUF measurements from January 2005 to December 2006 have been performed. MUF measurements collected during moderately disturbed days (17 ? Ap ? 32), disturbed days (32 < Ap ? 50) and very disturbed days (Ap > 50), have been used to test the long term prediction models (ASAPS, ICEPAC and SIRM&LKW), and the now casting models (SIRMUP&LKW and ISWIRM&LKW). The performances of the different prediction methods in terms of r.m.s are shown for selected range of geomagnetic activity and for each season.     

Occurrence of sporadic-E layer over the ionospheric station of Rome: Analysis of data for thirty-two years

Hourly systematic measurements of the highest frequency reflected by the sporadic-E layer (foEs) recorded at the Rome ionospheric observatory (Italy, 41.8° N, 12.5° E), were considered during the period January 1976–December 2007, to calculate the percentage of occurrence of sporadic-E layer with frequencies foEs greater than a given threshold value f_T, P(foEs > f_T).