Ionospheric F-region response to geomagnetic disturbances

The F2-region reaction to geomagnetic storms usually called as an ionospheric storm is a rather complicated event. It consists of so called positive and negative phases, which have very complicated spatial and temporal behavior. The main morphological features of ionospheric storms and the main processes governing their behavior were understood at the end of the 1900s and described in a series of review papers. During the recent decade there were many publications dedicated to the problem of ionospheric storms. In this paper a concept of ionospheric storm morphology and physics formulated at the end of the 1990s is briefly summarized and the most interesting results obtained in the 2000s are described. It is shown that the main features of the studies of the previous decade were: the use of GPS TEC data for analyzing the ionospheric storm morphology, attraction of sophisticated theoretical models for studying the processes governing ionospheric behavior in disturbed conditions, and accent to analysis of ionospheric behavior during prominent events (very strong and great geomagnetic storms). Also a special attention was paid to the pre-storm enhancements in foF2 and TEC.     

Ionospheric precursors observed at low latitudes around the time of koyna earthquake

On December 11, 1967 at 05:21 LT, an immense earthquake of magnitude 6.7 struck Koyna, the Indian province of Maharashtra. Its epicenter was located at geographic latitude 17.37°N and longitude 73.75°E with depth of about 3 km. Ground based measurements show variation in the critical frequency of ionospheric F2 layer (foF2) before and after the shock. In the present study the behavior of F2-region of ionosphere has been examined over the equatorial and low latitudinal region ionosphere during the month of December 1967 around the time of Koyna earthquake. For this purpose, the ionospheric data collected with the help of ground-based ionosondes installed at Hyderabad (located close to the earthquake epicenter) Ahmedabad, Trichirapulli, Kodaikanal and Trivendrum have been utilized. The upper and lower bound of Interquartile range (IRQ) are constructed to monitor the variations in foF2 other than day-to-day and diurnal pattern for finding the seismo-ionospheric precursors. Some anomalous electron density variations are observed between post midnight hours to local pre-noon hours at each station. These anomalies are strongly time dependent and appeared a couple of days before the main shock. The period considered in this study comes under the quiet geomagnetic conditions. Hence, the observed anomalies (which are more than the usual day-to-day variability) over all stations are likely to be associated with this imminent earthquake. The possible mechanism to explain these anomalies is the effect of seismogenic electric field generated just above the surface of earth within the earthquake preparation zone well before the earthquake due to emission of radioactive particles and then propagated upward, which perturbs the F-region ionosphere.     

Recent Progresses on Ionospheric Climatology Investigations

The ionosphere varies over multiple time scales, which are classified into two categories: the climatology and weather variations. In this national report, we give a brief summary of recent progresses in ionospheric climatology with focus on (1) the seasonal variations, (2) solar cycle effects, and (3) empirical modeling of the ionosphere. The seasonal variations of the ionosphere have been explored in many works to give a more detailed picture with regional and global features at various altitudes by analyzing the observation data from various sources and models. Moreover, a series of studies reported the response of the ionosphere to solar cycle variations, which revealed some novel and detailed features of solar activity dependence of ionospheric parameters at different altitudes. These investigations have improved our understanding on the states of the ionosphere and underlying fundamental processes, provided clues to future studies on ionospheric weather, and guided ionospheric modeling, forecasting and related applications.     

F2 layer response to geomagnetic disturbances in Eastern Asia under the low solar activity

In this paper, the peculiarities of ionospheric response to geomagnetic disturbances observed at the decay and minimum of solar activity (SA) in the period 2004–2007 are investigated with respect to different geomagnetic conditions. Data from ionospheric stations and results of total electron content (TEC) measurements made at the network of GPS ground-based receivers located within the latitude–longitude sector (20–70°N, 90–160°Е) are used in this study. Three groups of anomalous ionospheric response to geomagnetic disturbances have been observed during low solar activity. At daytime, the large-scale traveling ionospheric disturbances (LSTIDs) could generally be related to the main phase of magnetic storm. Quasi-two-days wavelike disturbances (WLDs) have been also observed in the main phase independent of the geomagnetic storm intensity. Sharp electron density oscillations of short duration (OSD) occurred in the response to the onset of both main and recovery phases of the magnetic storm in the daytime at middle latitudes. A numerical model for ionosphere–plasmasphere coupling was used to interpret the occurrence of LS TIDs. Results showed that the LSTIDs might be associated with the unexpected lifting of F2 layer to the region with the lower recombination rate by reinforced meridional winds that produces the increase of the electron density in the F2 layer maximum.     

Statistical properties of the most powerful solar and heliospheric disturbances

We present and discuss here the first version of a data base of extreme solar and heliospheric events. The data base contains now 87 extreme events mostly since 1940. An event is classified as extreme if one of the three critical parameters passed a lower limit. The critical parameters were the X-ray flux (parameter R), solar proton flux (parameter S) and geomagnetic disturbance level (parameter G). We find that the five strongest extreme events based on four variables (X-rays SEP, Dst, Ap) are completely separate except for the October 2003 event which is one the five most extreme events according to SEP, Dst and Ap. This underlines the special character of the October 2003 event, making it unique within 35 years. We also find that the events based on R and G are rather separate, indicating that the location of even extreme flares on the solar disk is important for geomagnetic effects. We also find that S = 3 events are not extreme in the same sense as R > 3 and G > 3 events, while S = 5 events are missing so far. This suggests that it might be useful to rescale the classification of SEP fluxes.     

Ionospheric wave signature of the American solar eclipse on 21 August 2017 in Europe

A total solar eclipse occurred on 21 August 2017, with the path of totality starting over the North Pacific Ocean, crossing North-America and ending over the Mid-Atlantic Ocean slightly North of the equator. As a result, a partial solar eclipse was observed as far away as the Western Europe. The ionospheric observatory in Dourbes, Belgium, was right on the edge of the partial eclipse and was exposed for a very short period of only few minutes just before the local sunset. High-resolution ionospheric measurements were carried out at the observatory with collocated digital ionosonde and GNSS receivers. The data analysis revealed a clear wave-like pattern in the ionosphere that can be seen arriving before the local onset of the eclipse. The paper details the analysis and provides a possible explanation of the observed phenomenon.     

Ionospheric disturbances detected by MEXART

The radio telescope MEXART was developed to make observations of interplanetary scintillation (IPS) produced by large scale disturbances associated with solar events. In this work it is shown that on occasion there are disturbances in the ionosphere that are related with these events and which cannot only contaminate the IPS but actually be the main contribution to the observed oscillations. This was the case of the event of 15 December 2006 observed by MEXART, which presented clear scintillation. The total electron content (TEC) of the ionosphere above Mexico was calculated for the same period. It was found that the variations in TEC were associated with the scintillations detected by MEXART.     

Solar activity variations of ionosonde measurements and modeling results

The time series of hourly electron density profiles N(h) obtained at several mid-latitude stations in Europe have been used to obtain N(h) profiles on a monthly basis and to extract both the expected bottomside parameters and a proxy of the ionospheric variability as functions of time and height. With these data we present advances on a “Local Model” technique for the parameters B0 and B1, its applicability to other ionospheric stations, to other bottomside ionospheric parameters, and to modeling the time/height variability of the profile. The Local Model (LM) is an empirical model based on the experimental results of the solar activity dependence of the daily and seasonal behavior of the above parameters. The LM improves the IRI-2001 prediction of the B0 and B1 by factor of two at mid-latitudes. Moreover, the LM can be used to simulate other ionospheric parameters and to build mean N(h) profiles and the deviations from them. The modeling of both the average N(h) profiles and their deviations is an useful tool for ionospheric model users who want to know both the expected patterns and their deviations.     

Ionospheric heating due to solar flares as measured by SROSS-C2 satellite

The data on thermal fluctuations of the topside ionosphere have been measured by Retarding Potential Analyser (RPA) payload aboard the SROSS-C2 satellite over the Indian region for half of the solar cycle (1995–2000). The data on solar flare has been obtained from National Geophysical Data Center (NGDC) Boulder, Colorado (USA) and other solar indices (solar radio flux and sunspot number) were download from NGDC website. The ionospheric electron and ion temperatures show a consistent enhancement during the solar flares. The enhancement in the electron temperature is 28–92% and for ion temperature it is 18–39% compared to the normal day’s average temperature. The enhancement of ionospheric temperatures due to solar flares is correlated with the variation of sunspot and solar radio flux (F10.7cm). All the events studied in the present paper fall in the category of subflare with almost same intensity. The ionospheric electron and ion temperatures enhancement have been compared with the IRI model values.     

Comparison of peak characteristics of F2 ionospheric layer over Tehran region at a low solar activity period with IRI-2001 and IRI-2007 models predictions

In the present work values of peak electron density (NmF2) and height of F2 ionospheric layer (hmF2) over Tehran region at a low solar activity period are compared with the predictions of the International Reference Ionosphere models (IRI-2001 and IRI-2007). Data measured by a digital ionosonde at the ionospheric station of the Institute of Geophysics, University of Tehran from July 2006 to June 2007 are used to perform the calculations. Formulations proposed by  and  are utilized to calculate the hmF2. The International Union of Radio Science (URSI) and International Radio Consultative Committee (CCIR) options are employed to run the IRI-2001 and IRI-2007 models. Results show that both IRI-2007 and IRI-2001 can successfully predict the NmF2 and hmF2 over Tehran region. In addition, the study shows that predictions of IRI-2007 model with CCIR coefficient has closer values to the observations. Furthermore, it is found that the monthly average of the percentage deviation between the IRI models predictions and the values of hmF2 and NmF2 parameters are less than 10% and 21%, respectively.     

Ionospheric reflection of the magnetic activity described by the index η

Differences in the external part of the vertical geomagnetic component point to the existence of local inhomogeneities in the magnetosphere or the ionosphere. Usually used magnetic indices are not sufficient to express the state of ionosphere, the common used global Kp index derived in the three-hour interval does not indicate much more rapidly changes appearing in ionosphere. Magnetic index η reflects ionospheric disturbances when other indices show very quiet conditions. Data of ionospheric characteristics (foE, foEs, h’E, h’F2) during 28-day long quiet day conditions (Kp = 0–2) in 2004 were analyzed. The correlations between strong local disturbances in ionosphere during very quiet days and high values of magnetic index η were found. The most sensitive to magnetic influence – ionospheric E layer data (foE characteristic) – reaches median deviations up to (+0.8 MHz and −0.8 MHz) during very low magnetic activity (Kp = 0–1). The high peaks (2–2.7) of the magnetic index η correlate in time with large local median deviations of foE. Such local deviations can suggest local inhomogeneities (vertical drifts) in the ionosphere. The correlation in space is not trivial. The strong peak of η is situated between the positive and negative deviations of foE. Additional observation is connected with correlation in time of the high η value with the negative median deviations of h’F2 (in some cases up to −90 km). The analysis was based on one-minute data recorded at each of 20 European Magnetic Observatories working in the INTERMAGNET network and from 19 ionosondes for 2004. Ionospheric data are sparse in time and in space in opposite to the magnetic data. The map of the magnetic indices can suggest the behavior of ionospheric characteristics in the areas where we have no data.     

Low-latitude geomagnetic response to the interplanetary conditions during very intense magnetic storms

The variations in the horizontal and declination components of the geomagnetic field in response to the interplanetary shocks driven by fast halo coronal mass ejections, fast solar wind streams from the coronal hole regions and the dynamic pressure pulses associated with these events are studied. Close association between the field-aligned current density (j_∥) and the fluctuations in the declination component (ΔD_ABG) at Alibag is found for intense storm conditions. Increase in the dawn-dusk interplanetary electric field (E_y) and ΔD_ABG are generally in phase. However, when the magnetospheric electric field is directed from dusk to dawn direction, a prominent scatter occurs between the two. It is suggested that low-latitude ground magnetic data may serve as a proxy for the interplanetary conditions in the solar wind.     

Trends in foF2 and the 24th solar activity cycle

The trends in foF2 are analyzed based on the data of Juliusruh and Boulder ionospheric stations. It is shown that using the traditional solar activity index F10.7 leads to an impossible trend in foF2 when the data for the 24th solar activity cycle are included into the analysis. It is assumed that the F10.7 index does not describe correctly the solar ultraviolet radiation variations in that cycle. A correction of this index using the Rz (sunspot number) and Ly (intensity of the Lyman-α line in the solar spectrum) is performed, and it is shown that in that case reasonable values of the foF2 trends are obtained.     

Global characteristics of ionospheric electric fields and disturbances during the first hours of magnetic storms

The ionospheric plasma density can be significantly disturbed during magnetic storms. In the conventional scenario of ionospheric storms, the negative storm phases with plasma density decreases are caused by neutral composition changes, and the positive storm phases with plasma density increases are often related to atmospheric gravity waves. However, recent studies show that the global redistribution of the ionospheric plasma is dominated primarily by electric fields during the first hours of magnetic storms. In this paper, we present the measurements of ionospheric disturbances by the DMSP satellites and GPS network during the magnetic storm on 6 April 2000. The DMSP measurements include the F region ion velocity and density at the altitude of ∼840 km, and the GPS receiver network provides total electron content (TEC) measurements. The storm-time ionospheric disturbances show the following characteristics. The plasma density is deeply depleted in a latitudinal range of ∼20° over the equatorial region in the evening sector, and the depletions represent plasma bubbles. The ionospheric plasma density at middle latitudes (20°–40° magnetic latitudes) is significantly increased. The dayside TEC is increased simultaneously over a large latitudinal range. An enhanced TEC band forms in the afternoon sector, goes through the cusp region, and enters the polar cap. All the observed ionospheric disturbances occur within 1–5 h from the storm sudden commencement. The observations suggest that penetration electric fields play a major role in the rapid generation of equatorial plasma bubbles and the simultaneous increases of the dayside TEC within the first 2 h during the storm main phase. The ionospheric disturbances at later times may be caused by the combination of penetration electric fields and neutral wind dynamo process.     

Ionospheric slab thickness – Analysis, modelling and monitoring

A review of the climatological and storm-time behaviour of the ionospheric slab thickness is presented based on long-time observations at a European mid-latitude site, Dourbes (50.1°N, 04.6°E), and on published results from other studies. An operational electron density and slab thickness monitoring system, established to provide real-time characterisation of the local ionospheric dynamics, is outlined together with some exemplary results.     

Ionospheric climate index as a driving parameter for the NeQuick model

Ionospheric delay is one of the significant error sources for global navigation satellite system (GNSS) positioning. GNSSs broadcast the coefficients of the ionospheric model to correct ionospheric delay for single-frequency users. A modified three-dimensional model (NeQuick G) based on the NeQuick climatological model is adopted for Galileo users. The NeQuick G model uses the effective ionization level (Az) instead of the sunspot number as the driving parameter. In this study, we introduce the ionospheric climate index (ICI) as a new driving parameter for the NeQuick model. In comparison, the ICI-driven NeQuick model has a better performance than the Az-driven NeQuick G model at both low and high latitudes. In addition, only one GNSS station at low latitudes is required to calculate the ICI, which would save maintenance costs and improve the efficiency of updating the broadcast coefficients. This model has potential application value for future upgrades of Galileo’s ionospheric broadcast model.     

Predicting indices of the ionosphere response to solar activity for the ascending phase of the 25th solar cycle

The international reference ionosphere, IRI, and its extension to plasmasphere, IRI-Plas, models require reliable prediction of solar and ionospheric proxy indices of solar activity for nowcasting and forecasting of the ionosphere parameters. It is shown that IRI prediction errors could increase for the F2 layer critical frequency foF2 and the peak height hmF2 due to erroneous predictions of the ionospheric global IG index and the international sunspot number SSN1 index on which IRI and IRI-Plas models are built. Regression relation is introduced to produce daily SSN1 proxy index from new time series SSN2 index provided from June 2015, after recalibration of sunspots data. To avoid extra errors of the ionosphere model a new solar activity prediction (SAP) model for the ascending part of the solar cycle SC25 is proposed which expresses analytically the SSN1 proxy index and the 10.7-cm radio flux F10.7 index in terms of the phase of the solar cycle, Φ. SAP model is based on monthly indices observed during the descending part of SC24 complemented with forecast of time and amplitude for SC25 peak. The strength of SC25 is predicted to be less than that of SC24 as shown with their amplitudes for eight types of indices driving IRI-Plas model.     

Behavior of F2-layer parameters and solar activity indices in the 24th cycle

The analysis of the behavior of the critical frequency foF2 during the 24th solar activity cycle (Danilov and Konstantinova, 2020a, c) is prolonged for two more months and the nighttime hours. In addition to the Rz and Ly-α indices used in the aforementioned papers for correction of the F10.7 index during the 24th cycle, the commonly used Mg II index is added. The results confirm the previous conclusions on the existence of the “vague” period with chaotic behavior of foF2 and the recovery of the negative trend in foF2 after 2008–2010. A comparison of the F10.7 index with three other SA indices (Ly-α, Rz, and Mg II) for the 22nd, 23rd, and 24th SA cycles is performed. It is shown that the relationship between F10.7 and other indices is close in the 22nd and 23rd cycles but differs from that in the 24th cycle. The corrected values of F10.7 in the 24th cycle are proposed for analysis of ionospheric trends during that cycle.     

Some ionospheric storm effects at equatorial and low latitudes

In this paper, the response of the equatorial and low latitude ionosphere to three intense geomagnetic storms occurred in 2002 and 2003 is reported. For that, critical frequency of F2-layer foF2 and the peak height hmF2 hmF2 for the stations Jicamarca (11.9°S), Ascension Is (7.92°S) and Tucuman (26.9°S) are used. The results show a “smoothing” of the Equatorial Anomaly structure during the development of the storms. Noticeable features are the increases in foF2 before the storm sudden commencement (SC) at equatorial latitudes and the southern crest of the Equatorial Anomaly. In some cases nearly simultaneous increases in foF2 are observed in response to the storm, which are attributed to the prompt electric field. Also, positive effects observed at equatorial and low latitudes during the development of the storm seem to be caused by the disturbance dynamo electric field due to the storm-time circulation. Increases in foF2 above the equator and simultaneous decreases in foF2 at the south crest near to the end of a long-duration main phase are attributed to equatorward-directed meridional winds. Decreases in foF2 observed during the recovery phase of storms are believed to be caused by composition changes. The results indicate that the prompt penetration electric field on the EA is important but their effect is of short lived. More significant ionospheric effects are the produced by the disturbance dynamo electric field. The role of storm-time winds is important because they modify the “fountain effect” and transport the composition changes toward low latitudes.     

Mid-latitude amplitude scintillation of GPS signals and GPS performance slips

Degradation of transionospheric radio signals and operation failures during ionospheric disturbances constitute a crucial factor of space weather influence on radio engineering satellite systems performance. We found that during the main phase of strong magnetic storms in 2000–2003 when the auroral oval expands into mid-latitudes, its southern boundary develops a region with intense small-scale electron density irregularities. Such irregularities may cause strong amplitude scintillations of GPS signals at both GPS operating frequencies. The another consequence of it was significant random GPS signal phase fluctuations, breaking-down of signal tracking, and sharp increasing of GPS positioning errors as a result.