亚暴起始位置和膨胀相持续时间对行星际磁场Bz分量的响应

磁层亚暴是太阳风–磁层–电离层耦合过程中的重要爆发性事件,其特性受太阳风参数的影响很大。本文利用对IMAGE卫星在2000 － 2005年观测到的4193个亚暴起始事件,统计研究了在不同的行星际磁场（IMF）B_z 条件下亚暴起始位置和膨胀相持续时间。结果表明,南向IMF发生的亚暴比北向IMF下发生的亚暴要多。南向IMF条件下亚暴AE指数最大值的平均值基本上>600 nT,并有随南向IMF持续时间增大而增大的趋势。北向IMF条件下亚暴AE指数最大值的平均值基本上<500 nT,并有随北向IMF持续时间增大而减小的趋势。亚暴的起始磁纬度基本上位于65° － 70°之间。当南向IMF或北向IMF的持续时间增大,超过80 min时,北半球的亚暴起始磁纬度会降低。亚暴起始磁地方时大部分位于22:15 － 23:15 MLT之间。但整体分布比较分散,显示不出特别清晰的随IMF B_z持续时间变化的趋势。相比于南向的IMF,北向IMF期间发生亚暴的平均膨胀相持续时间增大了将近10 min,表明南向IMF期间,亚暴强度虽然较大,但其膨胀相持续时间较短,亚暴能量释放和耗散的速度更快。      

太阳风扰动对极光电急流和环电流指数的影响

利用WIND卫星的太阳风观测数据和地磁活动指数, 研究了太阳风扰动对环电流SYM-H指数, 西向极光电急流AL指数和东向极光电急流AU指数的影响. 结果表明, 太阳风动压增长和减少能够同步或延迟引起地磁活动指数的强烈扰动, 其包括环电流指数的上升, 西向极光电急流指数的下降和东向极光电急流指数的上升. 太阳风动压的突然剧烈增加还能够触发超级亚暴和大磁暴. 太阳风动压脉冲引起的地磁效应具有复杂的表现形式, 这说明太阳风动压脉冲的地磁效应不仅与太阳风动压脉冲大小和持续时间有关, 还与磁层本身所处的状态有关. 时间尺度较长, 消耗能量较大的磁暴只有大的持续时间较长的太阳风动压脉冲才能激发.      

低轨道中等能量电子对地磁扰动的响应

低轨道高度上能量电子通量变化与地磁扰动程度密切相关.利用我国资源2号(ZY-2)03星空间环境监测分系统在轨工作期间所获得的能量电子探测数据,以及美国NOAA-15,NOAA-16,NOAA-17三颗卫星中等能量电子探测器自1998年以来积累的太阳同步轨道中等能量电子探测数据,结合地磁活动观测数据,对低轨道高度上中等能量电子对地磁扰动的响应特性进行了统计分析.结果表明,该区域的中等能量电子通量在磁暴、磁层亚暴期间有显著增强,增幅大小与地磁活动程度呈正相关关系,强磁暴期间增幅可达一个数量级左右,在响应时间上存在电子通量变化滞后于磁扰的时间特性.      

2021年10月11日地磁暴对两座变电站GIC的影响

近年来中国相继监测到地磁暴侵害电网、铁路轨道电路和油气管道系统产生的地磁感应电流（Geomagnetically Induced Current,GIC）数据,但是目前实测的GIC数据还相对较少。根据2021年10月9日日冕物质抛射事件（CME）产生的Kp指数为6的地磁扰动（Geomagnetic Disturbance,GMD）数据,500 kV阿拉坦变电站（48.7°N,116.8°E）和上河变电站（33.4°N,119.2°E）及输电系统的参数,分析了2021年10月11日地磁暴期间在两座变电站监测到的GIC数据以及输电系统参数对GIC量值的影响。结果表明:地磁暴在500 kV上河变电站产生的GIC比阿拉坦变电站GIC量值相对较大。分析结果说明,在这次磁暴事件中,输电线路导线电阻是影响变电站GIC的主要因素。      

TC-1和Geotail对亚暴过程中近地等离子体片磁扰动的联合观测

以2004年9月28日02:53:20 UT的亚暴为例, 通过TC-1在磁尾约12.5 R_e 和Geotail卫星在近地磁尾等离子体片约8～9 R_e的联合观测, 研究亚暴触发过程中近地磁尾等离子体片中等离子体波动特征. 结果表明, 亚暴触发区是近地磁尾中心等离子体片中较小的一个区域, 在亚暴触发区中低混杂不稳定性在近地磁尾等离子体片中存在, 准垂直传播的低混杂波发生在亚暴触发过程中, 而亚暴触发过程中近地磁尾等离子体片外边界区内的磁场偶极化信号和扰动都非常微弱. 在亚暴触发和亚暴膨胀相过程中出现了多次具有不同特征的磁场偶极化现象.      

行星际南向磁场事件与强磁暴

利用１９７８－１９８８年期间的太阳风和地磁资料对行星际磁场（ＩＭＦ）南向分量Ｂｓ事件（即Ｂｓ〉１０ｎＴ及其所驱动的错向电场ＶＢｓ〉５ｍＶ／ｍ、持续时间△Ｔ〉３ｈ的事件）与弱磁暴（Ｄｓｔ≤－１００ｎＴ）关系进行了分析。结果表明,１００％的Ｂｓ事件能能引起磁暴的发生,但其中只有８４％为强磁暴;强磁暴的发生都与较强的ＩＭＦ Ｂｓ活动密切相关,但只有６８％的强磁共伴随Ｂｓ事件而发生;Ｂｓ事件与强磁暴并不是     

第23太阳活动周太阳极紫外辐射与磁暴相关性分析

统计第23个太阳活动周内中等及以上强度（Dst_min<-50nT）的磁暴事件,线性拟合分析磁暴主相DDst_min和达到DDst_min前一个表征太阳极紫外辐射强度的F_10.7之间的相关性.结果表明:随着太阳极紫外辐射增强,DDst_min<-50nT的磁暴出现的总数增多,在弱、中等和强太阳极紫外辐射条件下,其数量分别为56,84和85;随着太阳极紫外辐射增强,强磁暴（-200nT ≤ Dst_min<-100nT）和大磁暴（Dst_min<-200nT）发生的数量和相对发生率呈增长趋势,尤其是大磁暴数目（1,4,12）和相对发生率（1.79%,4.76%,14.12%）明显呈增长趋势;大磁暴（|Dst_min|）与太阳极紫外辐射（F_10.7）之间存在中度正相关关系,其相关系数为0.532,并且主要体现在大磁暴（|Dst_min|）与强太阳极紫外辐射（F_10.7）之间的中度正相关性,其相关系数为0.582.大磁暴与强太阳极紫外辐射之间的相关性可为空间天气预报提供参考依据.      

基于Cluster观测的磁尾场向电流在南北半球投影位置的分布

利用Cluster四颗卫星的磁场探测数据计算磁尾场向电流并投影到极区电离层,研究其投影位置在南北半球的分布规律,统计过程中去除了强磁暴（磁暴主相Dst<–100 nT）期间的场向电流事件。结果显示:磁尾场向电流事件在极区投影位置的纬度分布具有明显的南北半球不对称性,北半球为单峰结构,南半球为双峰结构。在北半球投影到较低纬度（<64°）的场向电流事件数目明显多于南半球,并且所能达到的最低纬度更低;在南半球投影到较高纬度（>74°）的场向电流事件数目明显多于北半球,并且所能达到的最高纬度更高。地磁平静条件下（|AL|<100 nT）,磁尾场向电流密度随磁地方时（MLT）呈递增趋势,这一结果与低高度卫星在极区对I区场向电流的探测结果符合很好。研究结果表明,磁尾场向电流投影位置的纬度分布呈现出明显的南北不对称性,这与南北半球磁尾场向电流的空间分布以及磁层中磁场结构具有密切关系。      

太阳活动低年广州地区地磁扰动与电离层闪烁的统计特征分析

基于肇庆地磁台的地磁监测数据和广州气象卫星地面站建立的华南地区GPS电离层闪烁监测网的监测数据, 统计分析了2008年7月至2010年7月太阳活动低年期间广州地区地磁扰动与电离层闪烁的关系. 用肇庆台地磁水平分量H的变化量换算出肇庆地磁指数K, 以此来代表广州地区地磁扰动情况.分析结果表明, 磁暴/强地磁扰动对广州地区电离层闪烁的发生总体表现为抑制作用, 电离层闪烁主要发生在低K值期间, 而在K ≥ 4时电离层闪烁的发生呈下降趋势. 电离层闪烁发生率随季节和地磁活动的变化规律表现在, 春季的弱闪烁发生率、夜间中等以上闪烁发生率和夏季中等以上闪烁的发生率明显与地磁活动指数K相关, 即随$K$指数的增大而减小; 在秋季和冬季闪烁发生率与K指数变化无明显关系. 同时还综合分析了地磁与太阳活动的变化对电离层活动的影响, 广州地区闪烁主要发生在太阳活动较低的磁静日期间.      

2004年12月13日IMF北向期间极区亚暴电离层电动力学特征

利用KRM地磁反演方法, 结合北半球中高纬度地磁台站数据, 研究了2004年12月13日行星际磁场北向期间发生的亚暴事件, 极区电离层电动力学参量(电流矢量、等效电流函数以及电势)的分布特征. 结果表明, 在该亚暴膨胀相起始后, 午夜之前西向电集流急剧增强, 且等效电流体系表现为夜侧双涡, 同时伴随夜侧增强的南向电场. 由于极弱的直接驱动过程, 卸载过程引起的电离层效应得到清楚显示. 卸载过程在膨胀相期间起绝对主导性作用. 同时, 夜侧电导率的增强是电集流区域电流急剧增强的主要原因.      

镜模波识别方法研究及其在火星磁鞘中的应用

镜模波是温度各向异性等离子体中的一种波动结构,根据磁场和离子分布及波动特性可以进行识别.本文对比了只使用磁场数据与同时使用磁场及离子数据两种识别方法,分析了两类方法的特点.只使用磁场数据的方法基于磁场强度变化大、方向沿背景磁场的特征,通常使用磁场强度的波动幅度ΔB/|B|以及磁场变化方向与背景磁场的夹角θ_min,θ_max作为参数;同时使用磁场及粒子数据的方法利用的是磁场纵波特性、总压平衡和波动在等离子体坐标系下静止的特征.使用两种方法对MAVEN卫星在火星磁鞘内的数据进行识别,结果表明在某些情况下,只使用磁场数据会导致对镜模波的误判.通过研究改变上述参数阈值时识别结果的变化,发现当θ_min> 40°,θ_max < 40°,ΔB/|B|> 80%时,只用磁场数据可取得较好的识别效果.      

Diagnosing the storm-time ring current and other magnetospheric plasma domains using relativistic electron data

Ring current ions and relativistic electrons simultaneously measured on board MOLNIYA-1 are analyzed in comparison with the ground-based magnetometer data for the period of a strong magnetic storm (|D_st|max≈230 nT). Injection of >500 keV electrons into the slot region (L≈3) near equatorial plane is occurred on time scale ≈1 hour, when, during the magnetic storm maximum, the extreme low-latitude position of auroral electrojets is reached and ring current becomes more symmetrical. Positions of both the ring current maximum and electron intensity maximum (L_max) are consistent to our previous result: |D_st|_max = 2.75 • 10⁴/L⁴_max. An extreme storm-time low-latitude position of the west electrojet center (for amplitudes of |D_st|_max up to 600 nT) is shown to be in a good consistence with this empirical dependence. It is supposed the trapped radiation boundary collapses down to L≈L_max in the course of the storm main phase.     

Prediction tests by using “ISF” method for the geomagnetic disturbances

A so-called “ISF” prediction method for geomagnetic disturbances caused by solar wind storms blowing to the Earth is suggested. The method is based on a combined approach of solar activity, interplanetary scintillation (I) and geomagnetic disturbance observations during the period 1966–1982 together with the dynamics of solar wind storm propagation (S) and fuzzy mathematics (F). It has been used for prediction tests for 37 geomagnetic disturbance events during the descending solar activity phase 1984–1985, and was presented in 33rd COSPAR conference. Here, it has been improved by consideration of the three dimensional propagation characteristics of each event, the search for the best radio source and the influence of the southward components of interplanetary magnetic fields on the geomagnetic disturbances. It is used for prediction tests for 24 larger geomagnetic disturbance events that produced space anomalies during the period 1980–1999. The main results are: (1) for the onset time of the geomagnetic disturbance, the relative error between the observation, T_obs, and the prediction, T_pred, ΔT_pred/T_obs  10% for 45.8% of all events, 30% for 78.3% and >30% for only 21.7%; (2) for the magnetic disturbance magnitude, the relative error between the observation, ∑K_p,obs, and the prediction, ∑K_p,pred, Δ∑K_p,pred/∑K_p,obs  10% for 41.6% of all events, 30% for 79% and 45% for 100%. This shows that the prediction method described here has encouraging prospects for improving predictions of large geomagnetic disturbances in space weather events.     

Novel technology for detecting atmospheric disturbances using GPS. Instantaneous response of the ionosphere to a sudden commencement of the strong magnetic storms

We developed a new technology for global detection of atmospheric disturbances, on the basis of phase measurements of the total electron content (TEC) using the international GPS network. Temporal dependencies of TEC are obtained for a set of spaced receivers of the GPS network simultaneously for the entire set of visible satellites. These series are subjected to filtering in the selected range of oscillation periods using known algorithms for spatiotemporal analysis of signals. An “instantaneous” ionospheric response to the sudden commencement of the strong magnetic storms were detected. On the dayside of the Earth the largest value of the net response amplitude was found to be of order 0.4 TECU (2–3 % of the background TEC value), and the delay with respect to the SSC in mid-latitudes was about 360 s. In higher latitudes the delay goes up to 15 min. On the nightside these values are 0.2 TECU and 30 min, respectively. The velocity of the traveling disturbance from the middle to high latitudes on the dayside as well as from the dayside to the nightside was about 10–20 km/s.     

Storm-time formation of a relativistic electron belt and some relevant phenomena in other magnetospheric plasma domains

An empirical formula relating the strength of a storm given by its |Dst|_max with the L-coordinate of the peak of storm-injected relativistic electrons is one of a few well-confirmed quantitative relations found in the magneto-spheric physics. We successively extended a dataset of the formula’s basic storms with several events of high Dst-amplitude up to the highest observed |Dst|_max = 600 nT. Possible applying of the formula to the predicting of the ring-current plasma-pressure distribution and the lowest westward electrojet position for a storm are discussed. We have also analyzed the 2000–2001 years’ data on relativistic electrons from our instruments installed on EXPRESS-A (geosynchronous orbit; E_e = 0.8–6 MeV), Molniya-3 (h = 500 × 40 000 km, i = 63°; E_e = 0.8–5.5 MeV) and GLONASS (h = 20 000 km, i = 64°; E_e  l MeV) along with other correlated measurements: GOES series (E_e > 2 MeV), geomagnetic indices (Dst, AE, AL) and interplanetary parameters (solar wind, IMF). The goal is to investigate which outer conditions are most responsible for the high/low output of the storm-injected relativistic electrons. For the geosynchronous orbit, two factors are found as the necessary condition of the highest electron output: high and long-lasting substorm activity on a storm recovery phase and high velocity of solar wind. On the contrary, extremely low substorm activity surely observed during whole the storm recovery phase constitutes a sufficient condition of the non-increased after-storm electron intensity. For the first time found cases of the increased after-storm electron intensity observed at the inner L-shells with no simultaneously seen increase in the geosynchronous distances are presented.     

宇宙线强度变化与磁扰K类型

本文把1966—1983年期间发生的679个地磁暴进行了分类,利用统计方法分析了各类磁暴发生前后宇宙线强度的变化特征.突发急始脉冲发生后,宇宙线强度没有出现显著的变化;缓始型暴发生后,宇宙线强度出现Forbush下降,但下降幅度很小;急始型暴发生后,宇宙线强度出现十分明显的Forbush下降.当把急始型暴按K指数大小和持续时间分为5种类型,发现它们伴随的宇宙线Forbush下降是不一样的,其下降幅度随磁暴的增强而加大,下降的速率随磁暴的增强而加快,扰动的持续时间随磁暴的减弱而增加.      

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.     

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.     

Initial phase of mid-latitude aurora during strong geomagnetic storms

Analysing the initial mid-latitude aurora phase during strong geomagnetic storms we found that the initial phase of the mid-latitude aurorae observed at 630 nm emission during the strong geomagnetic storms on March 24, 1991, April 6, 2000, October 30 and November 20, 2003 is characterized by a short (∼1 h) wave-like disturbance. This disturbance corresponds to the beginning of main phase of the magnetic storms. The marked effect of the mid-latitude aurorae is analyzed using data on magnetosphere and ionosphere conditions in observation periods. The features of the dynamics of the 630 nm emission intensity and its connection with the dynamics of magnetospheric–ionospheric structures are considered. Possible excitation mechanisms of the atomic oxygen emission (630 nm) during these disturbances are discussed.