Statistical ionospheric E layer properties measured with the Cyprus Digisonde and comparisons with IRI predictions

This paper reports the diurnal, seasonal, and long term variability of the E layer critical frequency (foE) and peak height (hmE) derived from Digisonde measurements from 2009 to 2016 at the low-middle latitude European station of Nicosia, Cyprus (geographical coordinates: 35°N, 33°E, geomagnetic lat. 29.38°N, I = 51.7°). Manually scaled monthly median values of foE and hmE are compared with IRI-2012 predictions with a view to assess the predictability of IRI. Results show that in general, IRI slightly overestimates foE values both at low and high solar activity. At low solar activity, overestimations are mostly limited to 0.25?MHz (equivalent electron density, 0.775?×?10³?el/m^?3) but can go as high as 0.5?MHz (equivalent electron density, 3.1?×?10³?el/m^?3, during noon) around equinox. In some months, underestimations, though sporadic in nature, up to 0.25?MHz are noted (mostly during sunrise and sunset). At high solar activity, a similar pattern of over-/underestimation is evident. During the entire period of study, over-/under estimations are mostly limited to 0.25?MHz. In very few cases, these exceed 0.25?MHz but are limited to 0.5?MHz. Analysis of hmE reveals that: (1) hmE remains almost constant during ±2 to ±4?h around local noon, (2) hmE values are higher in winter than in spring, summer and autumn, (3) there are two maxima near sunrise and sunset with a noontime minimum in between. During the entire period of study, significant differences between observed hmE and the IRI predictions have been noted. IRI fails to predict hmE and outputs a constant value of 110?km, which is higher than most of the observed values. Over- and under estimations range from 3 to 13?km and from 0 to 3?km respectively.     

The response of the ionosphere over Ilorin to some geomagnetic storms

The effects of some geomagnetic storms on the F2 layer peak parameters over Ilorin, Nigeria (Lat. 8:53°N, Long. 4.5°E, dip angle, −2.96°) have been investigated. Our results showed that the highest intensity of the noon bite-out occurred during the March equinox and lowest during the June Solstice on quiet days. Quiet day NmF2 disturbances which appeared as a pre-storm enhancement, but not related to the magnetic storm event that followed were observed at this station. These enhancements were attributed to the modification of the equatorial electric field as a result of injection of the Auroral electric field to the low and equatorial ionosphere. For disturbed conditions, the morphology of the NmF2 on quiet days is altered. Daytime and nighttime NmF2 and hmF2 enhancements were recorded at this station. Decreases in NmF2 were also observed during the recovery periods, most of which appeared during the post-noon period, except the storm event of May 28–29. On the average, enhancements in NmF2 (i.e. Positive phases) are the prominent features of this station. Observations from this study also indicate that Dst, Ap and Kp which have been the most widely used indices in academic research in describing the behavior of geomagnetic storms, are not sufficient for storm time analysis in the equatorial and low latitude ionosphere.     

A different method to update monthly median hmF2 values

The height, h_mF2, and the electron density, N_mF2, of the F2 peak are key model parameters to characterize the actual state of the ionosphere. These parameters, or alternatively the propagation factor, M3000F2, and the critical frequency, f_oF2, of the F2 peak, which are related to h_mF2 and N_mF2, are used to anchor the electron density vertical profile computed with different models such as the International Reference Ionosphere ( Bilitza, 2002), as well as for radio propagation forecast purposes. Long time series of these parameters only exist in an inhomogeneous distribution of points over the surface of Earth, where dedicated instruments (typically ionosondes) have been working for many years. A commonly used procedure for representing median values of the aforementioned parameters all over the globe is the one recommended by the ITU-R ( ITU-R, 1997). This procedure, known as the Jones and Gallet mapping technique, was based on ionosondes measurements gathered from 1954 to 1958 by a global network of around 150 ionospheric stations (  and ). Even though several decades have passed since the development of that innovative work, only few efforts have been dedicated to establish a new mapping technique for computing h_mF2 and N_mF2 median values at global scale or to improve the old method using the increased observational database. Therefore, in this work three different procedures to describe the daily and global behavior of the height of the F2 peak are presented. All of them represent a different and simplified method to estimate h_mF2 and are based on different mathematical expressions. The advantages and disadvantages of these three techniques are analyzed, leading to the conclusion that the recommended procedure to represent h_mF2 is best characterized by a Spherical Harmonics expansion of degree and order equal to 15, since the differences between the h_mF2 values obtained with the Jones and Gallet technique and those obtained using the abovementioned procedure are of only 1%.     

An investigation of ionospheric F region response in the Brazilian sector to the super geomagnetic storm of May 2005

In this paper, we have investigated the responses of the ionospheric F region at equatorial and low latitude regions in the Brazilian sector during the super geomagnetic storm on 15–16 May 2005. The geomagnetic storm reached a minimum Dst of −263 nT at 0900 UT on 15 May. In this paper, we present vertical total electron content (vTEC) and phase fluctuations (in TECU/min) from Global Positioning System (GPS) observations obtained at Belém, Brasília, Presidente Prudente, and Porto Alegre, Brazil, during the period 14–17 May 2005. Also, we present ionospheric parameters h’F, hpF2, and foF2, using the Canadian Advanced Digital Ionosonde (CADI) obtained at Palmas and São José dos Campos, Brazil, for the same period. The super geomagnetic storm has fast decrease in the Dst index soon after SSC at 0239 UT on 15 May. It is a good possibility of prompt penetration of electric field of magnetospheric origin resulting in uplifting of the F region. The vTEC observations show a trough at BELE and a crest above UEPP, soon after SSC, indicating strengthening of nighttime equatorial anomaly. During the daytime on 15 and 16 May, in the recovery phase, the variations in foF2 at SJC and the vTEC observations, particularly at BRAZ, UEPP, and POAL, show large positive ionospheric storm. There is ESF on the all nights at PAL, in the post-midnight (UT) sector, and phase fluctuations only on the night of 14–15 May at BRAZ, after the SSC. No phase fluctuations are observed at the equatorial station BELE and low latitude stations (BRAZ, UEPP, and POAL) at all other times. This indicates that the plasma bubbles are generated and confined on this magnetically disturbed night only up to the low magnetic latitude and drifted possibly to west.     

Towards a new reference model of hmF2 for IRI

A numerical model of the peak height of the F2 layer, hmF2_top, is derived from the topside sounding database of 90,000 electron density profiles for a representative set of conditions provided by ISIS1, ISIS2, IK19 and Cosmos-1809 satellites for the period of 1969–1987. The model of regular hmF2 variations is produced in terms of local time, season, geomagnetic latitude, geodetic longitude and solar radio flux. No geomagnetic activity trends were discernible in the topside sounding data. The nighttime peak of hmF2_top evident for mid-latitudes disappears near the geomagnetic equator where a maximum of hmF2_top occurs at sunset hours when it can exceed 500 km at solar maximum. The hmF2 given by the IRI exceeds hmF2_top at the low solar activities. The hmF2_top, obtained by extrapolation of the first derivative of the topside profile to zero shows saturation similar to foF2 the greater the solar activity. The proposed model differs from hmF2 given by IRI based on M(3000)F2 to hmF2 conversion by empirical relationships in terms of foF2, foE and R12 with these quantities mapped globally by the ITU-R (former CCIR) from ground-based ionosonde data. The differences can be attributed to the different techniques of the peak 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.     

Local ionospheric electron density profile reconstruction in real time from simultaneous ground-based GNSS and ionosonde measurements

The purpose of the LIEDR (local ionospheric electron density profile reconstruction) system is to acquire and process data from simultaneous ground-based total electron content (TEC) and digital ionosonde measurements, and subsequently to deduce the vertical electron density distribution above the ionosonde’s location. LIEDR is primarily designed to operate in real time for service applications and, for research applications and further development of the system, in a post-processing mode. The system is suitable for use at sites where collocated TEC and digital ionosonde measurements are available. Developments, implementations, and some preliminary results are presented and discussed in view of possible applications.     

Validation of the Neustrelitz Global Model according to the low latitude ionosphere

Empirical modeling including empirical model for the total electron content (TEC) is important for the study of the ionosphere and practical applications. In this paper goodness of new Neustrelitz Global Model (then NGM) at low latitudes is studied. The NGM model includes such parameters as the maximal electron density (NmF2) and altitude of the maximum (hmF2). As of today, besides NGM there are several empirical models for NmF2 and hmF2. Therefore, a comparison of these parameters of the NGM model, not only with the experimental data, but also with two versions of the International Reference Ionosphere (the IRI model): IRI2001 and IRI-Plas would be instructive. Because the NGM model incorporates special factor describing the equatorial anomaly, the comparison in lower latitude areas is particularly interesting. As one can see from the presented example of the data from low latitude stations located in the northern and southern hemispheres near the Greenwich meridian, the NGM model may have certain advantages over the IRI model versions. In particular, NGM TEC is preferable regardless of solar activity level while NGM NmF2 is only preferable under high solar activity conditions. Next, NGM equivalent slab thickness of the ionosphere: τ(NGM) = TEC(NGM)/NmF2(NGM) has been calculated and tested to answer the question whether τ(NGM) can be used as a proxy of the slab thickness of the ionosphere for an empirical modeling. The answer is positive for the near equatorial stations and periods of high solar activity, and under such conditions predicted τ(NGM) can be used for deriving NmF2 from the experimental values of TEC(CODE) in real time.     

The equatorial ionization anomaly at the topside F region of the ionosphere along 75°E

Electron density measured by the Indian satellite SROSS C2 at the altitude of ∼500 km in the 75°E longitude sector for the ascending half of the solar cycle 22 from 1995 to 1999 are used to study the position and density of the equatorial ionization anomaly (EIA). Results show that the latitudinal position and peak electron density of the EIA crest and crest to trough ratios of the anomaly during the 10:00–14:00 LT period vary with season and from one year to another. Both EIA crest position and density are found to be asymmetric about the magnetic equator and the asymmetry depends on season as well as the year of observation, i.e., solar activity. The latitudinal position of the crest of the EIA and the crest density bears good positive correlation with F_10.7 and the strength of the equatorial electrojet (EEJ).