Multi-station observation of ionospheric irregularities over South Africa during strong geomagnetic storms

This paper presents results pertaining to the response of the mid-latitude ionosphere to strong geomagnetic storms that occurred from 31 March to 02 April 2001 and 07–09 September 2002. The results are based on (i) Global Positioning Systems (GPSs) derived total electron content (TEC) variations accompanying the storm, (ii) ionosonde measurements of the ionospheric electrodynamic response towards the storms and (iii) effect of storm induced travelling ionospheric disturbances (TIDs) on GPS derived TEC. Ionospheric data comprising of ionospheric TEC obtained from GPS measurements, ionograms, solar wind data obtained from Advanced Composition Explorer (ACE) and magnetic data from ground based magnetometers were used in this study. Storm induced features in vertical TEC (VTEC) have been obtained and compared with the mean VTEC of quiet days. The response of the mid-latitude ionosphere during the two storm periods examined may be characterised in terms of increased or decreased level of VTEC, wave-like structures in VTEC perturbation and sudden enhancement in hmF2 and h′F. The study reveals both positive and negative ionospheric storm effects on the ionosphere over South Africa during the two strong storm conditions. These ionospheric features have been mainly attributed to the travelling ionospheric disturbances (TIDs) as the driving mechanism for the irregularities causing the perturbations observed. TEC perturbations due to the irregularities encountered by the satellites were observed on satellites with pseudo random numbers (PRNs) 15, 17, 18 and 23 between 17:00 and 23:00 UT on 07 September 2002.     

Magnetospheric disturbances associated with the 13 December 2006 solar flare and their ionospheric effects over North-East Asia

We present an observational study of magnetospheric and ionospheric disturbances during the December 2006 intense magnetic storm associated with the 4В/Х3.4 class solar flare. To perform the study we utilize the ground data from North–East Asian ionospheric and magnetic observatories (60–72°N, 88–152°E) and in situ measurements from LANL, GOES, Geotail and ACE satellites. The comparative analysis of ionospheric, magnetospheric and heliospheric disturbances shows that the interaction of the magnetosphere with heavily compressed solar wind and interplanetary magnetic field caused the initial phase of the magnetic storm. It was accompanied by the intense sporadic E and F2 layers and the total black-out in the nocturnal subauroral ionosphere. During the storm main phase, LANL-97A, LANL 1994_084, LANL 1989-046 and GOES_11 satellites registered a compression of the dayside magnetosphere up to their orbits. In the morning–noon sector the compression was accompanied by an absence of reflections from ionosphere over subauroral ionospheric station Zhigansk (66.8°N, 123.3°E), and a drastic decrease in the F2 layer critical frequency (foF2) up to 54% of the quite one over subauroral Yakutsk station (62°N, 129.7°E). At the end of the main phase, these stations registered a sharp foF2 increase in the afternoon sector. At Yakutsk the peak foF2 was 1.9 time higher than the undisturbed one. The mentioned ionospheric disturbances occurred simultaneously with changes in the temperature, density and temperature anisotropy of particles at geosynchronous orbit, registered by the LANL-97A satellite nearby the meridian of ionospheric and magnetic measurements. The whole complex of disturbances may be caused by radial displacement of the main magnetospheric domains (magnetopause, cusp/cleft, plasma sheet) with respect to the observation points, caused by changes in the solar wind dynamic pressure, the field of magnetospheric convection, and rotation of the Earth.     

The solar wind control of the equatorial ionosphere dynamics during geomagnetic storms

The interplanetary magnetic field, geomagnetic variations, virtual ionosphere height h′F, and the critical frequency foF2 data during the geomagnetic storms are studied to demonstrate relationships between these phenomena. We study 5-min ionospheric variations using the first Western Pacific Ionosphere Campaign (1998–1999) observations, 5-min interplanetary magnetic field (IMF) and 5-min auroral electrojets data during a moderate geomagnetic storm. These data allowed us to demonstrate that the auroral and the equatorial ionospheric phenomena are developed practically simultaneously. Hourly average of the ionospheric foF2 and h′F variations at near equatorial stations during a similar storm show the same behavior. We suppose this is due to interaction between electric fields of the auroral and the equatorial ionosphere during geomagnetic storms. It is shown that the low-latitude ionosphere dynamics during these moderate storms was defined by the southward direction of the B_z-component of the interplanetary magnetic field. A southward IMF produces the Region I and Region II field-aligned currents (FAC) and polar electrojet current systems. We assume that the short-term ionospheric variations during geomagnetic storms can be explained mainly by the electric field of the FAC. The electric fields of the field-aligned currents can penetrate throughout the mid-latitude ionosphere to the equator and may serve as a coupling agent between the auroral and the equatorial ionosphere.     

Effects of TADs on the F region of the mid-latitude ionosphere during an intense geomagnetic storm

Based on observations of two ionosondes at Wuhan and Kokubunji, this paper presents effects of TADs on the daytime mid-latitude ionosphere during the intense geomagnetic storm on March 31, 2001. During a positive ionospheric storm, the start of the enhancement of the foF2 (F2 peak plasma frequency) at Wuhan lags that at Kokubunji by 15 min, which corresponds to the time interval of traveling atmospheric disturbances (TADs’) propagation from Kokubunji to Wuhan. Associated with the uplifting of the hmF2 (height of F2 peak) caused by TADs, it is observed by the two ionosondes that the F1 cusp becomes better developed. Therefore, during a geomagnetic storm, TADs originating from the auroral oval may have a strong influence on the shape of the electron density profile in the F1 region ionosphere at middle latitudes. It is highly likely that TADs are responsible for the evolution of the F1 cusp.     

Ionospheric response to solar and magnetospheric protons during January 15–22, 2005: EAGLE whole atmosphere model results

We present an analysis of the ionosphere and thermosphere response to Solar Proton Events (SPE) and magnetospheric proton precipitation in January 2005, which was carried out using the model of the entire atmosphere EAGLE. The ionization rates for the considered period were acquired from the AIMOS (Atmospheric Ionization Module Osnabrück) dataset. For numerical experiments, we applied only the proton-induced ionization rates of that period, while all the other model input parameters, including the electron precipitations, corresponded to the quiet conditions. In January 2005, two major solar proton events with different energy spectra and proton fluxes occurred on January 17 and January 20. Since two geomagnetic storms and several sub-storms took place during the considered period, not only solar protons but also less energetic magnetospheric protons contributed to the calculated ionization rates. Despite the relative transparency of the thermosphere for high-energy protons, an ionospheric response to the SPE and proton precipitation from the magnetotail was obtained in numerical experiments. In the ionospheric E layer, the maximum increase in the electron concentration is localized at high latitudes, and at heights of the ionospheric F2 layer, the positive perturbations were formed in the near-equatorial region. An analysis of the model-derived results showed that changes in the ionospheric F2 layer were caused by a change in the neutral composition of the thermosphere. We found that in the recovery phase after both solar proton events and the enhancement of magnetospheric proton precipitations associated with geomagnetic disturbances, the TEC and electron density in the F region and in topside ionosphere/plasmasphere increase at low- and mid-latitudes due to an enhancement of atomic oxygen concentration. Our results demonstrate an important role of magnetospheric protons in the formation of negative F-region ionospheric storms. According to our results, the topside ionosphere/plasmasphere and bottom-side ionosphere can react to solar and magnetospheric protons both with the same sign of disturbances or in different way. The same statement is true for TEC and foF2 disturbances. Different disturbances of foF2 and TEC at high and low latitudes can be explained by topside electron temperature disturbances.     

Peculiarities of the nighttime winter foF2 increase over Irkutsk

An analysis of properties and peculiarities of the nighttime winter foF2 increases (NWI) in the East Siberia is made on data of ionospheric station Irkutsk in the periods 1958–1992 and 2002–2009 and the empirical model of the F2 layer critical frequency under the geomagnetic quiet conditions deduced from these data (model Q-F2). It is revealed, that the NWI is the stable regularity of the quiet ionosphere over Irkutsk. The amplitude of the NWI (the difference between maximum and minimum foF2 values at night hours) is the greatest in December–January and nearly the same at low and middle solar activity. It is a peculiarity of the quiet ionosphere in the East Siberia. Maximum in night foF2 under quiet geomagnetic conditions is observed mainly after midnight (02-04 LT) and is shifted to predawn hours as solar activity increases. At low solar activity the quiet ionosphere at ∼02–04 LT shows the following properties: (a) the fluctuations of foF2 and hmF2 are in the reverse correlation but this dependence is weak; (b) very strong fluctuations of foF2 (|δfoF2| > 30%) occur seldom (∼4% of events) and almost all of them are positive; an example of very strong fluctuations of foF2 up to 60% can be an extreme increase in the foF2 on 19.12.2008; (c) the very strong enhancements of foF2 in the NWI maximum can be observed at the low geomagnetic activity, they occur more often during substorms but very seldom during geomagnetic storms. Possible reasons of these properties of NWI are discussed.     

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.     

Effects of geomagnetic storms on the bottomside ionospheric F region

The geomagnetic storm is a complex process of solar wind/magnetospheric origin. The variability of the ionospheric parameters increases substantially during geomagnetic storms initiated by solar disturbances. Various features of geomagnetic storm act at various altitudes in the ionosphere and neutral atmosphere. The paper deals with variability of the electron density of the ionospheric bottomside F region at every 10 km of altitude during intense geomagnetic storms with attention paid mainly to the distribution of the F1 region daytime ionisation. We have analysed all available electron density profiles from some European middle latitude stations (Chilton, Pruhonice, Ebro, Arenosillo, Athens) for 36 events that occurred in different seasons and under different levels of solar activity (1995–2003). Selected events consist of both depletion and increase of the F2 region electron density. For European higher middle and middle latitude the F1 region response to geomagnetic storm was found to be negative (decrease of electron density) independent on the storm effect on the F2 region. For lower middle latitude the F1 response is weaker and less regular. Results of the analysis also show that the maximum of the storm effect may sometimes occur below the height of the maximum of electron density (NmF2).     

Deformation of the ionosphere structure during the space weather events on October–November 2003

We present the spatial maps of the ionosphere–plasmasphere slab thickness τ (ratio of the vertical total electron content, TEC, to the F-region peak electron density, NmF2) during the intense ionospheric storms of October–November 2003. The model-assisted technology for estimate of the upper boundary of the ionosphere, hup, from the slab thickness components in the bottomside and topside ionosphere – eliminating the plasmasphere contribution of τ – is applied at latitudes 35° to 70°N and longitudes −10° to 40°E, from the data of 20 observatories of GPS-TEC and ionosonde networks, for selected days and hours of October and November 2003. The daily–hourly values of NmF2, hmF2 and TECgps are used as the constrained parameters for the International Reference Ionosphere extended to the plasmasphere, IRI-Plas, during the ionospheric quiet days, positive and negative storm phases for estimate of τ and hup. Good correlation has been found between the slab thickness and the upper boundary of the ionosphere for the intense ionospheric storms at October–November 2003. During the negative phase of the ionospheric storm, when the ionospheric plasma density is exhausted, the nighttime upper boundary of the ionosphere is greatly uplifted towards the magnetosphere tail, while the daytime upper boundary of the ionosphere is reduced below 500 km over the Earth.     

Updating of the IRI-1990 with respect to substorms

The method of updating the IRI-1990 model for the effects of large scale travelling ionospheric disturbances (LSTIDs) of the F-region have been established as a result of the statistical analysis of VI sounding data of European- Asian sector for 380 substorms. The effect of a LSTID occurring during a substorm is modelled as an additive correction to the quiet IRI foF2 and hmF2 values which depends on the time since the onset of the substorm and on the maximum AE-index value reached in it. In case of a sequence of substorms the additive corrections ΔfoF2 and ΔhmF2 are taking into account. The method for night-time ionosphere is described in detail, the listing of the FORTRAN-77 computer program is available on the floppy disk. The method is workable as correction to any ionospheric model.     

Variability of the ionosphere over Irkutsk at low solar activity

Statistical and spectral analyses are performed to investigate variations of two ionosphere F2 layer key parameters, the critical frequency (foF2) and the peak height (hmF2), that were measured over Irkutsk (52.5°N, 104.0°E) from December 2006 to January 2008 under solar minimum. The analyses showed that both parameters contain quasi-harmonic oscillations with periods of T_n = 24/n hours (n = 1–7), among which the diurnal (n = 1) and semidiurnal (n = 2) ones are the strongest. Seasonal variations are explored of mean and median values, spectrum, amplitude, and phase of the diurnal and semidiurnal components of foF2 and hmF2.     

Ionospheric response to the 3 August 2010 geomagnetic storm at mid and mid-high latitudes

The global ionospheric response to the geomagnetic storm occurred of 3 August 2010 is studied in terms of the ionospheric parameter foF2. Data from three longitudinal sectors (Asia/Pacific, Europe/Africa and America) are considered. Some new aspects of the storm time ionospheric behavior are revealed. Results of the analysis show that the main ionospheric effects of the storm under consideration are: (a) prior to the storm, Japanese, Australian and American stations show increases in foF2, irrespective of the local time. (b) During the main phase, the stations of mid latitudes of the American sector show positive disturbances (in the pre-dusk hours), which subsequently change to negative. (c) During the recovery phase of the magnetic storm long-duration positive disturbances are observed at mid-low latitudes of the African chain. Also positive disturbances are observed in the Australian sector. In the European sector long-duration negative disturbances are seen at mid-high latitudes during the last part of the recovery phase while at mid-low latitudes a positive disturbance is seen, followed by a negative disturbance. In general, the ionospheric storm effects show a clear hemispheric asymmetry.     

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.     

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.     

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.     

Effects observed in the equatorial and low latitude ionospheric F-region in the Brazilian sector during low solar activity geomagnetic storms and comparison with the COSMIC measurements

The main objective of the present investigation has been to compare the ionospheric parameters (NmF2 and hmF2) observed by two ground-based ionospheric sounders (one at PALMAS- located near the magnetic equator and the other at Sao Jose dos Campos-located in the low-latitude region) in the Brazilian sector with that by the satellite FORMOSAT-3/COSMIC radio occultation (RO) measurements during two geomagnetic storms which occurred in December 2006 and July 2009. It should be pointed out that in spite of increasing the latitude (to 10°) and longitude (to 20°) around the stations; we had very few common observations. It has been observed that both the peak electron density (NmF2) and peak height (hmF2) observed by two different techniques (space-borne COSMIC and ground-based ionosondes) during both the geomagnetic storm events compares fairly well (with high correlation coefficients) at the two stations in the Brazilian sector. It should be pointed out that due to equatorial spread F (ESF) in the first storm (December 2006) and no-reflections from the ionosphere during nighttime in the second storm (July 2009), we had virtually daytime data from the two ionosondes.     

An evaluation of the IRI-2007 storm time model at low latitude stations

This paper discusses the ability of the International Reference Ionosphere IRI-2007 storm time model to predict foF2 ionospheric parameter during geomagnetic storm periods. Experimental data (based on availability) from two low latitude stations: Vanimo (geographic coordinates, 2.7 °S, 141.3 °E, magnetic coordinates, 12.3 °S, 212.50 °E) and Darwin (geographic coordinates, 12.45 °S, 130.95 °E, magnetic coordinates, 22.9 °S, 202.7 °E) during nine storms that occurred in 2000 (Rz₁₂ = 119), 2001(Rz₁₂ = 111) and 2003 (Rz₁₂ = 64) are compared with those obtained by the IRI-2007 storm model. The results obtained show that the percentage deviation between the experimental and IRI predicted foF2 values during these storm periods is as high as 100% during the main and recovery phases. Based on the values of “relative deviation module mean” (RDMM) obtained (i.e. between 0.08 and 0.60), it is observed that there is a reasonable to poor agreement between measured foF2 values and the IRI-storm model prediction values during main and recovery phases of the storms under investigation. As a result, in addition to other studies that have been carried out from different sectors, more studies are required to be carried out. This will enable IRI community to improve on the present performance of the model. In general the IRI-storm model predictions follow normal trend of the foF2 measured values but does not reproduce well the measured values.     

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.     

Responses of ionospheric foF2 to geomagnetic activities in Hainan

Responses of low-latitude ionospheric critical frequency of F2 layer to geomagnetic activities in different seasons and under different levels of solar activity are investigated by analyzing the ionospheric foF2 data from DPS-4 Digisonde in Hainan Observatory during 2002–2005. The results are as follows: (1) the response of foF2 to geomagnetic activity in Hainan shows obvious diurnal variation except for the summer in low solar activity period. Generally, geomagnetic activity will cause foF2 to increase at daytime and decrease at nighttime. The intensity of response of foF2 is stronger at nighttime than that at daytime; (2) seasonal dependence of the response of foF2 to geomagnetic activity is very obvious. The negative ionospheric storm effect is the strongest in summer and the positive ionospheric storm effect is the strongest in winter; (3) the solar cycle has important effect on the response of foF2 to geomagnetic activity in Hainan. In high solar activity period, the diurnal variation of the response of foF2 is very pronounced in each season, and the strong ionospheric response can last several days. In low solar activity period, ionospheric response has very pronounced diurnal variation in winter only; (4) the local time of geomagnetic activities occurring also has important effect on the responses of foF2 in Hainan. Generally, geomagnetic activities occurred at nighttime can cause stronger and longer responses of foF2 in Hainan.     

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.