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.     

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

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

Empirical model of the TEC response to the geomagnetic activity over the North American region

The paper presents an empirical model of the total electron content (TEC) response to the geomagnetic activity described by the K_p-index. The model is built on the basis of TEC measurements covering the region of North America (50°W–150°W, 10°N–60°N) for the period of time between October 2004 and December 2009. By using a 2D (latitude-time) cross-correlation analysis it is found that the ionospheric response to the geomagnetic activity over the considered geographic region and at low solar activity revealed both positive and negative phases of response. The both phases of the ionospheric response have different duration and time delay with respect to the geomagnetic storm. It was found that these two parameters of the ionospheric response depend on the season and geographical latitude. The presence of two phases, positive and negative, of the ionospheric response imposed the implementation of two different time delay constants in order to properly describe the two different delayed reactions. The seasonal dependence of the TEC response to geomagnetic storms is characterized by predominantly positive response in winter with a short (usually ∼5–6 h) time delay as well as mainly negative response in summer with a long (larger than 15 h) time delay. While the TEC response in March and October is more close to the winter one the response in April and September is similar to the summer one.     

On the mid-latitude ionospheric storm association with intense geomagnetic storms

The bulk association between ionospheric storms and geomagnetic storms has been studied. Hemispheric features of seasonal variation of ionospheric storms in the mid-latitude were also investigated. 188 intense geomagnetic storms (Dst ≤ 100 nT) that occurred during solar cycles 22 and 23 were considered, of which 143 were observed to be identified with an ionospheric storm. Individual ionospheric storms were identified as maximum deviations of the F2 layer peak electron density from quiet time values. Only ionospheric storms that could clearly be associated with the peak of a geomagnetic storm were considered. Data from two mid-latitude ionosonde stations; one in the northern hemisphere (i.e. Moscow) and the other in the southern hemisphere (Grahamstown) were used to study ionospheric conditions at the time of the individual geomagnetic storms. Results show hemispheric and latitudinal differences in the intensity and nature of ionospheric storms association with different types of geomagnetic storms. These results are significant for our present understanding of the mechanisms which drive the changes in electron density during different types of ionospheric storms.     

Investigation of ionospheric response to two moderate geomagnetic storms using GPS–TEC measurements in the South American and African sectors during the ascending phase of solar cycle 24

The responses of the ionospheric F region using GPS–TEC measurements during two moderate geomagnetic storms at equatorial, low-, and mid-latitude regions over the South American and African sectors in May 2010, during the ascending phase of solar cycle 24, are investigated. The first moderate geomagnetic storm studied reached a minimum Dst value of −64 nT at 1500 UT on 02 May 2010 and the second moderate geomagnetic storm reached a minimum Dst value of −85 nT at 1400 UT on 29 May 2010. In this paper, we present vertical total electron content (VTEC) and phase fluctuations (in TECU/min) from Global Positioning System (GPS) observations from the equatorial to mid-latitude regions in the South American and African sectors. Our results obtained during these two moderate geomagnetic storms from both sectors show significant positive ionospheric storms during daytime hours at the equatorial, low-, and mid-latitude regions during the main and recovery phases of the storms. The thermospheric wind circulation change towards the equator is a strong indicator that suggests an important mechanism is responsible for these positive phases at these regions. A pre-storm event that was observed in the African sector from low- to the mid-latitude regions on 01 May 2010 was absent in the South American sector. This study also showed that there was no generation or suppression of ionospheric irregularities by storm events. Therefore, knowledge about the suppression and generation of ionospheric irregularities during moderate geomagnetic storms is still unclear.     

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).     

Spread-F occurrence during geomagnetic storms near the southern crest of the EIA in Argentina

This work presents, for the first time, the analysis of the occurrence of ionospheric irregularities during geomagnetic storms at Tucumán, Argentina, a low latitude station in the Southern American longitudinal sector (26.9°S, 294.6°E; magnetic latitude 15.5°S) near the southern crest of the equatorial ionization anomaly (EIA). Three geomagnetic storms occurred on May 27, 2017 (a month of low occurrence rates of spread-F), October 12, 2016 (a month of transition from low to high occurrence rates of spread-F) and November 7, 2017 (a month of high occurrence rates of spread-F) are analyzed using Global Positioning System (GPS) receivers and ionosondes. The rate of change of total electron content (TEC) Index (ROTI), GPS Ionospheric L-band scintillation, the virtual height of the F-layer bottom side (h'F) and the critical frequency of the F2 layer (foF2) are considered. Furthermore, each ionogram is manually examined for the presence of spread-F signatures.The results show that, for the three events studied, geomagnetic activity creates favorable conditions for the initiation of ionospheric irregularities, manifested by ionogram spread-F and TEC fluctuation. Post-midnight irregularities may have occurred due to the presence of eastward disturbance dynamo electric fields (DDEF). For the May storm, an eastward over-shielding prompt penetration electric field, (PPEF) is also acting. A possibility is that the PPEF is added to the DDEF and produces the uplifting of the F region that helps trigger the irregularities. Finally, during October and November, strong GPS L band scintillation is observed associated with strong range spread-F (SSF), that is, irregularities extending from the bottom-side to the topside of the F region.     

Dayside ionospheric uplift during strong geomagnetic storms as detected by the CHAMP,SAC-C,TOPEX and Jason-1 satellites

Ionosphere response to severe geomagnetic storms that occurred in 2001–2003 was analyzed using data of global ionosphere maps (GIM), altimeter data from the Jason-1 and TOPEX satellites, and data of GPS receivers on-board CHAMP and SAC-C satellites. This allowed us to study in detail ionosphere redistribution due to geomagnetic storms, dayside ionospheric uplift and overall dayside TEC increase. It is shown that after the interplanetary magnetic field turns southward and intensifies, the crests of the equatorial ionization anomaly (EIA) travel poleward and the TEC value within the EIA area increases significantly (up to ∼50%). GPS data from the SAC-C satellite show that during the main phase of geomagnetic storms TEC values above the altitude of 715 km are 2–3 times higher than during undisturbed conditions. These effects of dayside ionospheric uplift occur owing to the “super-fountain effect” and last few hours while the enhanced interplanetary electric field impinged on the magnetopause.     

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.     

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.     

中低纬电离层不同季节地磁暴响应的事例分析

利用中国中低纬台站漠河（53.5°N,122.3°E）、北京（40.3°N,116.2°E）、武汉（30.5°N,114.2°E）和三亚（18.3°N,109.6°E）的电离层观测数据,对比分析了4个台站电离层参数在2015年不同季节4个地磁扰动事件期间的变化特征.结果表明,4个磁暴事件期间电离层的响应特征并不完全一致,有着明显的季节特征,春季、夏季和秋季电离层以负相扰动为主,冬季以正相扰动为主.分析发现,中性成分O/N₂的降低与电离层负相扰动有关,但三亚地区的负相扰动还与扰动发电机电场相关.正相扰动的机制在不同事件中并不相同,穿透电场可能是引起春季磁暴事件期间电离层短时正暴效应的原因,而冬季长时间的正暴效应则是扰动电场和中性风共同作用的结果.      

Responses of the African equatorial ionization anomaly (EIA) to some selected intense geomagnetic storms during the maximum phase of solar cycle 24

This study investigates the morphology of the GPS TEC responses in the African Equatorial Ionization Anomaly (EIA) region to intense geomagnetic storms during the ascending and maximum phases of solar cycle 24 (2012–2014). Specifically, eight intense geomagnetic storms with Dst ≤ ?100 nT were considered in this investigation using TEC data obtained from 13 GNSS receivers in the East African region within 36–42°E geographic longitude; 29°N–10°S geographic latitude; ± 20°N magnetic latitude. The storm-time behavior of TEC shows clear positive and negative phases relative to the non-storm (median) behavior, with amplitudes being dependent on the time of sudden commencement of the storm and location. When a storm starts in the morning period, total electron content increases for all stations while a decrease in total electron content is manifested for a storm that had its sudden commencement in the afternoon period. The TEC and the EIA crest during the main phase of the storm is significantly impacted by the geomagnetic storm, which experiences an increase in the intensity of TEC while the location and spread of the crest usually manifest a poleward expansion.     

The impact of large solar events on the total electron content of the ionosphere at mid latitudes

Ionospheric disturbances associated with solar activity may occur via two basic mechanisms. The first is related to the direct impact on the ionosphere of EUV photons from a flare, and the second by prompt electric field penetration into the magnetosphere during geomagnetic storms. In this paper we examine the possibility that these two mechanisms may have an impact at mid latitudes by calculating the total electron content (TEC) from GPS stations in Mexico during several large X-ray flares. We have found that indeed large, complex flares, which are well located, may affect the mid latitude ionosphere. In fact, in the solar events of July 14, 2000 and April 2001 storms, ionospheric disturbances were observed to increase up to 138 and 150 TECu, respectively, due to the influence of EUV photons. Also, during the solar events of July 2000, April 2001, Halloween 2003, January 2005 and December 2006, there are large ionospheric disturbances (up to 393 TECu in the Halloween Storms), due to prompt penetration electric field, associated with CME producing geomagnetic storm.     

Effects of CME and CIR induced geomagnetic storms on low-latitude ionization over Indian longitudes in terms of neutral dynamics

This paper presents the response of the ionosphere during the intense geomagnetic storms of October 12–20, 2016 and May 26–31, 2017 which occurred during the declining phase of the solar cycle 24. Total Electron Content (TEC) from GPS measured at Indore, Calcutta and Siliguri having geomagnetic dips varying from 32.23°N, 32°N and 39.49°N respectively and at the International GNSS Service (IGS) stations at Lucknow (beyond anomaly crest), Hyderabad (between geomagnetic equator and northern crest of EIA) and Bangalore (near magnetic equator) in the Indian longitude zone have been used for the storms. Prominent peaks in diurnal maximum in excess of 20–45 TECU over the quiet time values were observed during the October 2016 storm at Lucknow, Indore, Hyderabad, Bangalore and 10–20 TECU for the May 2017 storm at Siliguri, Indore, Calcutta and Hyderabad. The GUVI images onboard TIMED spacecraft that measures the thermospheric O/N₂ ratio, showed high values (O/N₂ ratio of about 0.7) on October 16 when positive storm effects were observed compared to the other days during the storm period. The observed features have been explained in terms of the O/N₂ ratio increase in the equatorial thermosphere, CIR-induced High Speed Solar Wind (HSSW) event for the October 2016 storm. The TEC enhancement has also been explained in terms of the Auroral Electrojet (AE), neutral wind values obtained from the Horizontal Wind Model (HWM14) and equatorial electrojet strength from magnetometer data for both October 2016 and May 2017 storms. These results are one of the first to be reported from the Indian longitude sector on influence of CME- and CIR-driven geomagnetic storms on TEC during the declining phase of solar cycle 24.     

基于因果卷积与LSTM网络的电离层总电子含量预报

电离层总电子含量（TEC）不仅是分析电离层形态的关键参数之一,同时为导航及定位等空间应用系统消除电离层附加时延提供重要支撑。由于电离层TEC的时空变化特征,本文融合因果卷积和长短时记忆网络,以太阳活动指数F_10.7、地磁活动指数Dst和电离层TEC历史数据作为特征输入,构建深度学习模型,实现提前24 h预报电离层TEC。进一步利用2005－2013年连续9年的CODE TEC数据,全面评估了模型在北京站（40°N,115°E）、武汉站（30.53°N,114.36°E）和海口站（20.02°N,110.38°E）的预报性能。结果显示不同太阳活动条件下三个站的TEC值与真实测量值的相关系数都大于0.87,均方根误差大都集中在0～1 TECU以内,且模型预报精度与纬度、太阳、地磁活动程度、季节变化相关。与仅由长短时记忆网络构成的预报模型相比,本实验模型均方根误差降低了15%,为电离层TEC预报模型的实际应用提供了参考。      

基于深度学习递归神经网络的电离层总电子含量经验预报模型

利用行星际太阳风参数与太阳活动指数、地磁活动指数、电离层总电子含量格点化地图数据,首次基于一种能处理时间序列的深度学习递归神经网络（Recurrent Neural Network,RNN）,建立提前24h的单站电离层TEC预报模型.对北京站（40°N,115°E）的预测结果显示,RNN对扰动电离层的预测误差低于反向传播神经网络（Back Propagation Neural Network,BPNN）0.49～1.46TECU,将太阳风参数加入预报因子模型后对电离层正暴预测准确率的提升可达16.8%.RNN对2001和2015年31个强电离层暴预报的均方根误差比BPNN低0.2TECU,将太阳风参数加入RNN模型可使31个事件的平均预报误差降低0.36～0.47TECU.研究结果表明深度递归神经网络比BPNN更适用于电离层TEC的短期预报,且在预报因子中加入太阳风数据对电离层正暴的预报效果有明显改善.      

Analysis of ionospheric TEC response to solar and geomagnetic activities at different solar activity stages

Studying the relationship of total electron content (TEC) to solar or geomagnetic activities at different solar activity stages can provide a reference for ionospheric modeling and prediction. On the basis of solar activity indices, geomagnetic activity parameters, and ionospheric TEC data at different solar activity stages, this study analyzes the overall variation relationships of solar and geomagnetic activities with ionospheric TEC, the characteristics of the quasi-27-day periodic oscillations of the three variables at different stages, and the delayed TEC response of solar activity by conducting correlation analysis, Butterworth band-pass filtering, Fourier transform, and time lag analysis. The following results are obtained. (1) TEC exhibits a significant linear relationship with solar activity at different solar activity stages. The correlation coefficients |R| are arranged as follows: |R|_EUV > |R|_F10.7 > |R|sunspot number. No significant linear relationship exists between TEC and geomagnetic activity parameters (|R| < 0.35). (2) TEC, solar activity indices, and geomagnetic activity parameters have a period of 10.5 years. The maximum amplitudes of the Fourier spectrum for TEC and solar activity indices are nearly 27 days and those of geomagnetic activity parameters are nearly 27 and 13.5 days. (3) The deviations of the quasi-27-day significant periodic oscillation of TEC and solar activity indices are consistent. (4) No evident relationship exists between the quasi-27-day periodic oscillation of TEC and geomagnetic activity parameters. (5) The delay time of TEC for the 10.7 cm solar radio flux and extreme ultraviolet is always consistent, whereas that for sunspot number varies at each stage.     

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.     

2009年7月22日日全食期间电离层参量的变化

利用多个电离层垂测站的数据和IGS-TEC数据资料, 结合日地空间环境指数, 分析了2009年7月22日日全食期间中国地区电离层参量(反射回波最低频率fmin及f0F2和TEC)的变化特征. 结果表明, 日食发生后fmin迅速降低, 日食结束后fmin迅速恢复到正常水平; 在食甚时刻附近, f0F2和TEC出现明显的降低, 显示了明显的光食效应. 日食结束后5~6 h, f0F2和TEC出现不同程度的正扰动, 在驼峰区更明显; 日食结束后9~10 h, f0F2和TEC出现较显著的负扰动. 由于此次日食发生时伴随着中等强度的磁暴和低纬电场穿透等空间天气事件, 给此次日食电离层效应的深入分析带来很大困难.      

Traveling ionospheric disturbances observed at South African midlatitudes during the 29–31 October 2003 geomagnetically disturbed period

This paper presents traveling ionospheric disturbances (TIDs) observations from GPS measurements over the South African region during the geomagnetically disturbed period of 29–31 October 2003. Two receiver arrays, which were along two distinct longitudinal sectors of about 18°-20° and 27°-28° were used in order to investigate the amplitude, periods and virtual propagation characteristics of the storm induced ionospheric disturbances. The study revealed a large sudden TEC increase on 28 October 2003, the day before the first of the two major storms studied here, that was recorded simultaneously by all the receivers used. This pre-storm enhancement was linked to an X-class solar flare, auroral/magnetospheric activities and vertical plasma drift, based on the behaviour of the geomagnetic storm and auroral indices as well as strong equatorial electrojet. Diurnal trends of the TEC and foF2 measurements revealed that the geomagnetic storm caused a negative ionospheric storm; these parameters were depleted between 29 and 31 October 2003. Large scale traveling ionospheric disturbances were observed on the days of the geomagnetic storms (29 and 31 October 2003), using line-of-sight vertical TEC (vTEC) measurements from individual satellites. Amplitude and dominant periods of these structures varied between 0.08–2.16 TECU, and 1.07–2.13 h respectively. The wave structures were observed to propagate towards the equator with velocities between 587.04 and 1635.09 m/s.