Time delay of ionospheric TEC storms to geomagnetic storms and pre-storm disturbance events in East Asia

Using the TEC data at Beijing (39.61°N, 115.89°E)/Yakutsk (62.03°N, 129.68°E) stations of East Asia regions and relevant geomagnetic data from 2010 to 2017, we have studied the time delay of ionospheric storms to geomagnetic storms and compare it with our previous results of Taoyuan (25.02°N, 121.21°E) station (Zhang et al., 2020). The data shows a well-known local time dependence of the time delay, and seasonal dependences are different at these stations. In addition, there is no correlation between the time delay and the magnetic storm intensity /solar activity, except the time delay of negative storms has weakly negative dependence on the solar activity. Comparing with the results of Taoyuan station which is located at EIA region in East Asia, we find that the time delay increases nonlinearly as the latitude decreases due to different ionospheric backgrounds at these places. Moreover, the pre-storm disturbance events are found to have similar statistical characteristics as the pre-storm enhancement in Europe middle latitudes (Bure?ová and La?tovi?ka, 2007). By subtracting the common features of the pre-storm disturbance events, we preliminarily infer that auroral activity might be main driver of the pre-storm disturbance events.     

Statistical study of the time delay of ionospheric TEC storms to geomagnetic storms in Taoyuan,Taiwan

We have studied the time delay of ionospheric storms to geomagnetic storms at a low latitude station Taoyuan (25.02°N, 121.21°E), Taiwan using the Dst and TEC data during 126 geomagnetic storms from the year 2002 to 2014. In addition to the known local time dependence of the time delay, the statistics show that the time delay has significant seasonal characteristics, which can be explained within the framework of the seasonal characteristics of the ionospheric TEC. The data also show that there is no correlation between the time delay and the intensity of magnetic storms. As for the solar activity dependence of the time delay, the results show that there is no relationship between the time delay of positive storms and the solar activity, whereas the time delay of negative storms has weakly negative dependence on the solar activity, with correlation coefficient −0.41. Especially, there are two kinds of extreme events: pre-storm response events and long-time delay events. All of the pre-storm response events occurred during 15–20 LT, manifesting the Equator Ionospheric Anomaly (EIA) feature at Taoyuan. Moreover, the common features of the pre-storm response events suggest the storm sudden commencement (SSC) and weak geomagnetic disturbance before the main phase onset (MPO) of magnetic storms are two main possible causes of the pre-storm response events. By analyzing the geomagnetic indices during the events with long-time delay, we infer that this kind of events may not be caused by magnetic storms, and they might belong to ionospheric Q-disturbances.     

Ionospheric disturbances under low solar activity conditions

The paper is focused on ionospheric response to occasional magnetic disturbances above selected ionospheric stations located at middle latitudes of the Northern and Southern Hemisphere under extremely low solar activity conditions of 2007–2009. We analyzed changes in the F2 layer critical frequency foF2 and the F2 layer peak height hmF2 against 27-days running mean obtained for different longitudinal sectors of both hemispheres for the initial, main and recovery phases of selected magnetic disturbances. Our analysis showed that the effects on the middle latitude ionosphere of weak-to-moderate CIR-related magnetic storms, which mostly occur around solar minimum period, could be comparable with the effects of strong magnetic storms. In general, both positive and negative deviations of foF2 and hmF2 have been observed independent on season and location. However positive effects on foF2 prevailed and were more significant. Observations of stormy ionosphere also showed large departures from the climatology within storm recovery phase, which are comparable with those usually observed during the storm main phase. The IRI STORM model gave no reliable corrections of foF2 for analyzed events.     

Validation of the STORM model in IRI-2001 at a high latitude station

The International Reference Ionosphere IRI-2001 model contains geomagnetic activity dependence based on an empirical storm time ionospheric correction (STORM model). An extensive validation of the STORM model for the middle latitude region has been performed. In this paper the ability of the STORM model to predict foF2 values at high latitudes is analyzed. For this, ionosonde data obtained at Base Gral. San Martin (68.1°S, 293°E) are compared with those obtained by the IRI-2001 model with or without storm correction during four geomagnetic storms that occurred in 2000 (Rz₁₂ = 117) and 2001 (Rz₁₂ = 111). The results show that predicted values with the STORM model follow the behaviour of foF2 experimental data better than without the STORM model. The relative deviation between measured and predicted foF2 reaches values of up to 24% and 43% with and without the STORM model in IRI-2001, during the main phase of the storms. In order to explain increases of electron density that occurred prior to the storm onset and also decreases of electron density observed during the first part of the recovery of the storm, possible physical mechanisms are discussed.     

Geomagnetic pulsation over conjugate locations during geomagnetic storms and substorm

Geomagnetic data collected during magnetic storm over magnetically conjugate pair (according to IGRF 2000) of high latitude stations viz., Maitri (70° 45′ S, 11° 42′ E) and Tromso (69° 40′ N, 18° 56′ E) reveal that amplitudes of Pc6 pulsation characteristically differ. The amplitude obtained from horizontal magnetic field for the Pc6 pulsation frequency range between 0.6 and 1.6 mHz significantly differs in time corresponding to peak amplitude. The relative differences in its time of occurrence found to gradually increase around initial phase of storm and remain exactly out of phase at peak amplitude of storm. Thence, it is found to be in agreement in phase gradually until storm unwinds. This indicates that simultaneous amplitude of Pc6 pulsation at conjugate pair of stations and its time of occurrence could be a key factor to infer storm arrivals somewhat prior to its peak effects. The emphasize remains on prediction of storm arrivals only by utilizing ground based magnetometer observations. However, it is necessary to understand differences on the basis of weak, moderate, strong, and super strong cases and more exactly how they behave along the line of magnetic Meridian. Nevertheless, the analysis implies that geo-effective magnetic ejecta/clouds/CIRs/sheaths/CMEs/ICMEs giving rise to geomagnetic storm can be predicted ahead of its peak effects by having magnetometer data over conjugate locations.     

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.     

Large geomagnetic storms of extreme solar event periods in solar cycle 23

During extreme solar events such as big flares or/and energetic coronal mass ejections (CMEs) high energy particles are accelerated by the shocks formed in front of fast interplanetary coronal mass ejections (ICMEs). The ICMEs (and their sheaths) also give rise to large geomagnetic storms which have significant effects on the Earth’s environment and human life. Around 14 solar cosmic ray ground level enhancement (GLE) events in solar cycle 23 we examined the cosmic ray variation, solar wind speed, ions density, interplanetary magnetic field, and geomagnetic disturbance storm time index (D_st). We found that all but one of GLEs are always followed by a geomagnetic storm with D_st  −50 nT within 1–5 days later. Most(10/14) geomagnetic storms have D_st index  −100  nT therefore generally belong to strong geomagnetic storms. This suggests that GLE event prediction of geomagnetic storms is 93% for moderate storms and 71% for large storms when geomagnetic storms preceded by GLEs. All D_st depressions are associated with cosmic ray decreases which occur nearly simultaneously with geomagnetic storms. We also investigated the interplanetary plasma features. Most geomagnetic storm correspond significant periods of southward B_z and in close to 80% of the cases that the B_z was first northward then turning southward after storm sudden commencement (SSC). Plasma flow speed, ion number density and interplanetary plasma temperature near 1 AU also have a peak at interplanetary shock arrival. Solar cause and energetic particle signatures of large geomagnetic storms and a possible prediction scheme are discussed.     

宇宙线变化特征在强地磁暴预报中的应用

利用宇宙线中子探测数据定性分析了地面宇宙线多台站之间的相互联系以及大磁暴与宇宙线之间的响应关系. 以Irkutsk和Oulu宇宙线台站为例, 运用小波去噪技术提高数据的稳定性. 结果表明, 相同世界时条件下, 两站宇宙线通量相关性在事件发生时较高; 而相同地方时条件下, 相关性则在平静期较高. 进一步采用相同地方时条件对不同宇宙线台站的通量在平静期和扰动期的相对变化进行分析, 选取2004年7月强地磁暴典型事例进行直观分析, 发现大地磁暴前Irkutsk和Oulu台站的宇宙线相对通量发生明显差异, 可以尝试作为强地磁暴宇宙线先兆特征. 通过对2001年3月至2005年5月的强磁暴和中强磁暴进行统计, 得到与强地磁暴相关的适当宇宙线相对差异阈值. 将得到的阈值对2005年9月至2011年12月所有强磁暴及中强磁暴进行验证, 总成功率达到87.5%, 误报率为35.7%, 结果较好.      

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.     

Relationship between substorm activity and the interplanetary medium parameters during the main phase of strong magnetic storms

In the present paper dependences of substorm activity on the solar wind velocity and southward component (Bz) of interplanetary magnetic field (IMF) during the main phase of magnetic storms, induced by the CIR and ICME events, is studied. Strong magnetic storms with close values of Dst_min?≈??100?±?10?nT are considered. For the period of 1979–2017 there are selected 26 magnetic storms induced by the CIR and ICME (MC?+?Ejecta) events. It is shown that for the CIR and ICME events the average value of the AE index (AE_aver) at the main phase of magnetic storm correlates with the solar wind electric field. The highest correlation coefficient (r?=?0.73) is observed for the magnetic storms induced by the CIR events. It is found that the AE_aver for magnetic storms induced by ICME events, unlike CIR events, increases with the growth of average value of the southward IMF Bz module. The analysis of dependence between the AE_aver and average value of the solar wind velocity (Vsw_aver) during the main phase of magnetic storm shows that in the CIR events, unlike ICME, the AE_aver correlates on the Vsw_aver.     

Moderate geomagnetic storms of January 22–25, 2012 and their influences on the wave components in ionosphere and upper stratosphere-mesosphere regions

Moderate geomagnetic storms occurred during January 22–25, 2012 period. The geomagnetic storms are characterized by different indices and parameters. The SYM-H value on January 22 increased abruptly to 67 nT at sudden storm commencement (SSC), followed by a sharp decrease to −87 nT. A second SSC on January 24 followed by a shock on January 25 was also observed. These SSCs before the main storms and the short recovery periods imply the geomagnetic storms are CME  -driven. The sudden jump of solar wind dynamic pressure and IMF _Bz$B_z$ are also consistent with occurrence of CMEs. This is also reflected in the change in total electron content (TEC) during the storm relative to quiet days globally. The response of the ionospheric to geomagnetic storms can also be detected from wave components that account for the majority of TEC variance during the period. The dominant coherent modes of TEC variability are diurnal and semidiurnal signals which account upto 83% and 30% of the total TEC variance over fairly exclusive ionospheric regions respectively. Comparison of TEC anomalies attributed to diurnal (DW1) and semidiurnal (SW2) tides, as well as stationary planetary waves (SPW1) at 12 UTC shows enhancement in the positive anomalies following the storm. Moreover, the impact of the geomagnetic storms are distinctly marked in the daily time series of amplitudes of DW1, SW2 and SPW1. The abrupt changes in amplitudes of DW1 (5 TECU) and SW2 (2 TECU) are observed within 20°S–20°N latitude band and along 20°N respectively while that of SPW1 is about 3 TECU. Coherent oscillation with a period of 2.4 days between interplanetary magnetic field and TEC was detected during the storm. This oscillation is also detected in the amplitudes of DW1 over EIA regions in both hemispheres. Eventhough upward coupling of quasi two day wave (QTDWs) of the same periodicity, known to have caused such oscillation, are detected in both ionosphere and upper stratosphere, this one can likely be attributed to the geomagnetic storm as it happens after the storm commencement. Moreover, further analysis has indicated that QTDWs in the ionosphere are strengthened as a result of coherent oscillation of interplanetary magnetic field with the same frequency as QTDWs. On the otherhand, occurrences of minor SSW and geomagnetic storms in quick succession complicated clear demarcation of attribution of the respective events to variability of QTDWs amplitudes over upper stratosphere.     

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.     

Seasonal variations of GPS derived TEC at three different latitudes of the southern hemisphere during geomagnetic storms

Data from the archive of the International GNSS Services (IGS) were used to study the seasonal variations of Total Electron Content (TEC) over three stations located at different latitudes in the southern hemisphere during the geomagnetic storms of 11 January, 6 April, 8 June, and 13 October 2000, representing storms that occurred in summer, autumn equinox, winter and spring equinox, respectively. The percentage TEC deviation with respect to reference values differs substantially from season to season. A strong seasonal anomaly and clear equinoctial asymmetry in TEC response to the storms were observed. Weak and short-lived positive TEC deviations as well as strong and long-lasting negative trends were observed in summer storm during the main and recovery phases respectively over the high and low latitudes whereas in winter storm, the highest positive TEC deviations was recorded during the main phase over the entire latitudes. TEC enhancement dominated all the stations during the autumn (March) equinox storm while TEC depletion was majorly observed during the spring (September) equinox. All these variations find their explanations in the thermospheric composition change and circulation. Future work with direct or modeled measurement of atomic Oxygen to molecular Nitrogen ratio (O/N₂), large number of storms and other possible factors such as variations in storm’s intensity and local time dependence of the storm onset is expected to validate the observations in this study.     

Latitudinal variation of F2-region response to geomagnetic disturbance

The effect of geomagnetic storms on the F2 region was studied by calculating the deviation, ΔfoF2, of foF2 during 40 magnetic storms, ranging from moderate (Dst < −50 nT) to very intense (Dst < −200 nT) of the 21st solar cycle. In order to study the variation of storm-time foF2 with latitude, season and storm strength, ionosonde data were obtained from eight stations spanning a latitudinal range of +60–−60°. The stations chosen lay in a narrow longitudinal range of 140–151°, so that local time difference between the stations is practically negligible. The features exhibited by positive and negative phases were essentially different. The storm time ΔfoF2 clearly exhibited a latitudinal variation and this variation were found to be coupled with the seasonal variation. As for the variation with storm intensity, though ΔfoF2 was found to vary even between two storms of almost equal intensity, the amplitude of a positive or negative phase, |ΔfoF2max| showed a distinct upper limit for each intensity category of storms.     

2000年7月和2003年10月大磁暴期间东亚地区中低纬电离层的GPS TEC的响应研究

利用武汉电离层观象台研制的GPS TEC的现报方法及现报系统,对东亚地区GPS台网的观测数据进行处理分析,特别对2000年7月14-18日和2003年10月28日至11月1日两次特大磁暴期间的数据进行了对比考察,文中分析了两次磁暴间的电离层响应,得到对应不同磁暴时段电离层TEC的不同变化情况,着重揭示了TEC赤道异常峰的压缩和移动以及赤道异常随时间的压缩—反弹—恢复的过程,并结合高纬电离层的部分响应机制进行了说明,结果显示,两次磁暴期的电离层响应表现出了各自不同的特点,从而反映出因季节变化引起的高纬电离层暴时能量注入的不同而造成的全球性电离层扰动的不同形态,由此看出,磁暴期间电离层TEC的变化直接与太阳扰动发生的时间及其对高纬电离层的耦合有关,若短时期内连续发生多次磁暴,则电离层反应更加复杂,不能简单地当做单一磁暴叠加处理。     

Relation of decline characteristics of 2–4.6 MeV protons in SEP events to solar wind parameters

The shape of the particle flux decline in solar energetic particle (SEP) events is of particular importance in understanding the propagation of energetic particles in the interplanetary medium. Power-law time profiles indicate the dominance of diffusive propagation, whereas exponential-law decline emphasizes convection transport and adiabatic deceleration. Values obtained theoretically for the decay time in the latter case are reasonably close to the fitted slopes in nearly half of all events when the solar wind speed stays nearly constant. Dependencies of characteristic decay time τ and spectral index γ on environmental plasma parameters are considered. Parts of exponential-law declines when solar wind speed: (a) decreases with time, (b) is constant, and (c) increases with time through the interval are analyzed separately. Both average values and dispersions of size distributions of τ for these three groups markedly differ in accordance with theoretical expectations.     

Some effects in the upper atmosphere during geomagnetic storms

This paper examines the response of the high latitude ionosphere–thermosphere system during two intense geomagnetic storms. For that, data taken by instruments on board Dynamic Explorer 2 taken at heights of the F2-layer are used. These results represent a comparison of simultaneous measurements of storm disturbances in gas composition, electron density and temperature in common local time sectors. Documented are an increase in electron temperature and a decrease in electron density; increases both in electron temperature and electron density; and the correlation between electron density decreases and increases in the ratio N₂/O. It is noticed that the decrease in electron density is sometimes due to an increase in the molecular nitrogen density N₂ and not always is attributed both to the increase in N₂ density and the simultaneous decrease in the atomic oxygen density.     

Variations in cosmic ray cutoff rigidities during the great geomagnetic storm of November 2004

A very strong interplanetary and magnetospheric disturbance observed on 7–13 November 2004 can be regarded as one of the strongest events during the entire period of space observations. In this paper we report on the studies of cosmic ray cutoff rigidity variations during 7–13 November 2004 showing how storm conditions can affect the direct cosmic ray access to the inner magnetosphere. Effective cutoff rigidities have been calculated for selected points on the ground by tracing trajectories of cosmic ray particles through the magnetospheric magnetic field of the “storm-oriented” Tsyganenko 2003 model. Cutoff rigidity variations have also been determined by the spectrographic global survey method on the basis of experimental data of the neutron monitor network. Relations between the calculated and experimental cutoff rigidities and the geomagnetic Dst-index and interplanetary parameters have been investigated. Correlation coefficients between the cutoff rigidities obtained by the trajectory tracing method and the spectrographic global survey method have been found to be in the limits 0.76–0.89 for all stations except the low-latitude station Tokyo (0.35). The most pronounced correlation has been revealed between the cutoff rigidities that exhibited a very large variation of ∼1–1.5 GV during the magnetic storm and the Dst index.     

Qualitative estimation of magnetic storm efficiency in producing relativistic electron flux in the Earth’s outer radiation belt using geomagnetic pulsations data

It is well known that during many but not all of the geomagnetic storms enhanced fluxes of high-energy electrons are observed in the outer radiation belt. Here we examine relativistic (>2 MeV) electron fluxes measured by GOES at the synchronous orbit and on-ground observations of two types of ULF pulsations during 30 magnetic storms occurred during 1996–2000. To characterize the effectiveness of the chosen magnetic storms in producing relativistic electron fluxes, following to (Reeves, G.D., McAdams, K.L., Friedel, R.H.W., O’Brien, T.R. Acceleration and loss of relativistic electrons during geomagnetic storms. Geophys. Res. Lett. 30, doi:10.1029/2002GL016513, 2003), we calculate a ratio of the maximum daily-averaged electron flux measured during the recovery phase, to the mean pre-storm electron flux. A storm is considered an effective one if its ratio exceeds 2. We compare behavior of Pi1 and Pc5 geomagnetic pulsations during effective and non-effective storms and find a tendency for a storm efficiency to be higher when the mid-latitude Pi1 pulsations are observed for a long time during the magnetic storm main phase. We note also that the prolonged powerful Pc5 pulsation activity during the recovery phase of a magnetic storm is the necessary condition for the storm effectiveness. To interpret the found dependences, we suggest that there are two prerequisites for generating relativistic electron populations during a storm: (1) the availability of seed electrons in the magnetosphere, and Pi1 emissions are indicators of the mid-energy electron interaction with the ionosphere and (2) acceleration of the seed electrons to MeV energies, and interaction of electrons with the MHD wave activity in the Pc5 range is one of the most probable mechanisms proposed in the literature for this purpose.