1998年5月空间天气大事件的地磁场响应

地磁场与1998年5月空间大事件相对应的是5月1日至16日发生的大磁暴（k=8)。磁暴主相开始的几个小时伴随有丰富的Pc型地磁脉动，包括P c2,Pc3,Pc4等。在增暴的恢复相，甚至还有Pc5巨型脉动，在行星标磁场Bz由北向转向南向时，磁暴主相开始；南向分量达到最大值后大约2小时，地磁H分量达到最小值，恢复相开始，并且，这次磁暴与太阳风电场也存在一定的对应关系。     

基于SAMPEX卫星观测的南大西洋异常区高能质子动态分布特征

利用SAMPEX卫星1992年7月至2004年6月19～27MeV高能质子数据对南大西洋异常区的分布特征进行研究, 发现南大西洋异常区高能质子分布随高度及F_10.7的变化十分显著. 在540±25km高度上, 地磁较为平静时期南大西洋异常区高能质子微分通量随着F_10.7的增大而减小, 同时在F_10.7≥115sfu时减小趋势较为平缓. 对中等及以上磁暴进行统计分析发现, 磁暴期间南大西洋异常区高能质子微分通量和SYM-H指数的绝对值存在明显的反相关关系, 且地磁暴对南大西洋异常区高能质子微分通量存在明显的持续影响效应. 磁暴发生期间高能质子微分通量明显减少. 磁暴恢复相及其之后高能质子微分通量呈现较为显著的恢复过程.      

影响卫星故障的空间天气分析

基于美国国家地球物理数据中心(NGDC) 2384例和中国19颗卫星的263例卫星故障信息, 结合1963-2012年小时平均的多种空间环境数据, 定量分析了三种卫星故障发生期间的空间要素特征, 探讨单粒子锁定(SEU)、表面充电致静电放电(ESD)和内部深层充电所致电子引起的电磁脉冲(ECEMP)与空间天气事件的可能联系, 得出以下主要结论. (1)大部分SEU和ECEMP发生于空间天气平静时, 但在其前后3日内地磁活动达到了磁暴水平, 相对来说比例最大的发生在Dst_min之后第3日 (48～72h). (2) ESD受地磁活动和高能电子通量影响明显. 与磁暴、相对论电子通量增强事件的季节性相对应, 两分点附近ESD和ECEMP的发生率高; 93.6% 的 ESD发生前后72h内地磁活动达到磁暴水平, 故障发生时间均匀分布在 Dst_min前0～48h 和后0～24h; 54.9%的ESD 发生时处于地磁暴期(Dst <-30nT), 以-50～-30nT的小磁暴水平居多; 40.6%的ESD发生于高能电子通量高水平期(≥ 103pfu, 1pfu =1cm-2·s-1·sr-1), 81.9%的ESD发生前后72h 内高能电子通量峰值≥ 103pfu, 发生率最高时段为电子通量峰值前 48～72h. (3)高能电子对中国同步轨道卫星的SEU影响明显, 42.5% 故障发生 时高能电子通量≥ 103pfu, 故障在峰值前48～72h和峰值后48～72h 的发生概率相当, 约为23.0%. (4)同步轨道卫星SEU受太阳质子事件的影响相对较大, 22.5%的中国同步轨道卫星故障发生前后72h内发生了太阳质子事件, 季节性不明显.      

第23太阳活动周中等地磁暴行星际源的统计分析

统计了第23太阳活动周(1996--2006年)发生的183次中等强度地磁暴(-100 nT < Dst ≤ -50 nT)的行星际源,分析了中等磁暴的年分布状况以及引起中等磁暴的不同行星际结构在太阳活动周中的分布特征,同时,与强磁暴行星际源的分布状况做了对比分析,主要的统计分析结果如下. (1)共转相互作用区CIR与行星际日冕物质抛射ICME在中等磁暴中具有同等重要的作用,且在ICME中,具有磁云结构和非磁云结构的ICME在引起中等磁暴的能力方面也基本相同,但带有鞘层结构的ICME在引起中等磁暴中具有更重要的作用. (2)中等磁暴在极大年(2001年)和下降年(2003年)发生次数最多,与地磁活动的双峰年对应,在极小年(1996和2006年)发生次数最少,与地磁活动低年对应,在其他年份分布较平均. (3)中等磁暴在太阳活动极大年主要由ICME引起,在上升年和下降年CIR在其中起主要作用,且下降年基本是上升年的两倍,而对于强磁暴而言,ICME始终是最重要的行星际源.      

1998年10月18——-19日磁暴主相的行星际源分析

采用高时间分辨率的地磁指数SYM-H, 同时考虑日地连线引力平衡点(L1点)太阳风地磁效应的滞后性, 精确分析了1998年10月18---19日大磁暴主相的行星际源. 分析结果表明, 磁暴主相的行星际源仅为行星际激波和行星际日冕物质抛射之间的太阳风(Sheath), 磁云对磁暴主相没有贡献. 这个磁暴事例的研究表明, 行星际磁场南向分量与太阳风动压的乘积是影响磁暴主相发展的关键参数.      

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

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

1991年3月太阳活动引起的地磁扰动

利用北京白家疃台和马道峪台的地磁资料分析了1991年3月下旬太阳活动引起的24－27日的磁暴，从中纬记录看，这次磁暴有其特殊的形态：即急始幅度特别大，磁暴最强活动是急始后即刻开始的极区暴影响而产 生的：磁暴过程中多个时段出现脉动，其中以周期为650秒左右的脉动尤为明显；其极化方向主要表现为左旋偏振，同时夹有右旋偏振，有时则无明显规律，谱分析的结果表明，整个磁暴功率谱呈幂律P(f)=f^-2,0分布，对引起这次磁暴可能的行星际扰动进行了讨论。     

基于稀疏地磁数据的中低纬度磁暴感应地电场的算法研究

强磁暴产生的地磁感应电流是诱发电力系统灾害性事故的关键因素,预防与控制地磁感应电流的必要途径是评估其在系统中的水平.掌握磁暴感应地电场的时空分布是评估所在区域电网地磁感应电流的基础,可为研究与防治电网磁暴灾害提供参考.在中国现有稀疏分布的地磁台观测数据基础上,建立了电离层等效球面元电流系统模型,并将其与复镜像法结合,建立了磁暴感应地电场和地磁场的计算模型.利用磁暴期间实测的地磁台站数据进行了计算,算例结果与平面波理论计算结果的对比证明了算法的正确性与有效性.计算结果可直接应用于电力系统的地磁感应电流计算,为分析电力系统磁暴灾害风险提供参考.      

SABRE雷达及其观测的Pc5地磁脉动Ⅰ.统计分析结果

本文简要介绍了相干散射雷达SABRE及其观测Pc5地磁脉动的工作原理。1985—1989年观测数据的统计分析结果表明,相位向极传播脉动事件的出现率峰值在当地时间午夜附近,而向赤道传播事件的出现率峰值则在18:00LT左右。文中给出了上述两类脉动事件出现率的日变化,与地磁活动的相关,与背景电离层等离子体对流特性的关系等统计分析结果。     

X级以上耀斑与地磁效应关系研究

统计研究了2010年1月至2012年12月期间所有与耀斑爆发相伴生的日冕物质抛射(CME) 引发的地磁暴事件. 结果表明, 对于CME源区其主要分布在日面 45°E-45°W, 占总数的78.95%, 且西半球比东半球多, 即源区位于西半球的CME易产生地磁效应; X级耀斑与地磁效应的关联性更高, 60.0%的 X级耀斑在其爆发后的2～3天内观测到地磁暴, 而其他级别的耀斑与地磁效应的关联性低得多, 均不足10%; 通过对此期间日面爆发的所有X级耀斑研究分析后发现, 对于源区位于日面东经45°E-45°W 的X级耀斑, 若在其爆发过程中没有大尺度日面扰动, 则无伴生CME且后续产生地磁效应的可能性很低. 由此提出一种通过分析日面观测数据进行地磁暴预报的方法.      

Mid-Latitude Pc1, 2 Pulsations Induced by Magnetospheric Compression in the Maximum and Early Recovery Phase of Geomagnetic Storms

We investigate the properties of interplanetary inhomogeneities generating long-lasting mid-latitude Pc1, 2 geomagnetic pulsations. The data from the Wind and IMP 8 spacecrafts, and from the Mondy and Borok midlatitude magnetic observatories are used in this study. The pulsations under investigation develop in the maximum and early recovery phase of magnetic storms. The pulsations have amplitudes from a few tens to several hundred pT andlast more than seven hours. A close association of the increase (decrease) in solar wind dynamic pressure (Psw) with the onset or enhancement (attenuation or decay) of these pulsations has been established. Contrary to high-latitude phenomena, there is a distinctive feature of the interplanetary inhomogeneities that are responsible for generation of long-lasting mid-latitude Pc1, 2. It is essential that the effect of the quasi-stationary negative Bz-component of the interplanetary magnetic field on the magnetosphere extends over 4 hours. Only then are the Psw pulses able to excite the above-mentioned type of mid-latitude geomagnetic pulsations. Model calculations show that in the cases under study the plasmapause can form in the vicinity of the magnetic observatory. This implies that the existence of an intense ring current resulting from the enhanced magnetospheric convection is necessary for the Pc1, 2 excitation. Further, the existence of the plasmapause above the observation point (as a waveguide) is necessary for long-lasting Pc1 waves to arrive at the ground.      

Statistical characteristics of low latitude Pc4s: Influence of solar wind and IMF parameters

Pc4 signatures for the year 2013, extracted from geomagnetic north–south and east–west components of induction coil magnetometer (LEMI 30) from low latitude station Desalpar (DSP), operated by Institute of Seismological Research (ISR), India have been investigated vis-à-vis the prevalent interplanetary parameters (IMF) as well as the geomagnetic activity indices. A clear dominance of Pc4-5 (467 events) over Pc3 (17 events) is observed. Local time variation of Pc4 shows a peak in the noon sector in both X and Y components. Our investigations show that the dominant peak frequency is 10 mHz at low latitude region. Correlations with solar wind and IMF parameters illustrate highest occurrence of Pc4 for a solar wind speed of 300–400 km/s and average IMF B field of 3–6 nT. The amplitude of Pc4s at DSP shows an increase with increasing solar wind speed, plasma density, solar wind dynamic pressure and average B field which is also reflected in the trend of frequency variation of these pulsations. We report that IMF clock angle at low latitude does not have influence on Pc4 occurrence. Based on the characteristics of these events, detected in latitudinally distributed stations from low and mid-latitudes from northern and southern hemisphere, we infer that modes were compressional, which could be driven by K-H instability or solar wind dynamic pressure, as compressional modes can propagate to low latitude with little attenuation.     

Solar wind dependence of the Pc1 wave activity

There are a host of factors influencing the excitation of Pc1 geomagnetic pulsations, which are ULF waves in the frequency range between 0.2 and 5 Hz. We have studied carefully the dependence of the pearl-type Pc1 activity at Sodankylä, Finland (L = 5.1) on the plasma density N in front of the magnetosphere, the bulk velocity V of the solar wind, and the intensity B of the IMF. The result is as follows: high values of N and reduced values of V are favorable to appearance of Pc1, whereas the dependence of Pc1 activity on B is practically absent. We also show that the probability of Pc1 occurrence decreases with the interplanetary electric field, and increases with solar wind impact pressure and with the plasma to magnetic pressure ratio “beta”.     

低纬Pc3脉动和电离层电流体系

假设在低纬地区有一根磁力线振荡,从基本方程出发,推算出地面各点的磁场扰动情况,模拟低纬Pc3脉动。对比低纬地区Pc3脉动的观测结果,可以看出用这种理论数值模拟的结果基本上与观测事实相符,用它可以初步解释低纬Pc3脉动的频率特性和偏振特点。     

A coordinated study of field-aligned currents and Pc5 ULF waves during ejecta 1997

We show examples of long period Pc5 magnetic field pulsations near field-aligned current (FAC) regions in the high-latitude magnetosphere, observed by INTERBALL-Au, and coordinated with POLAR, GOES-9 and ground-based observations during 11 January and 11 April 1997. Identification of corresponding magnetosphere regions and subregions is provided by electrons and protons in the energy-range of 0.01–100 keV measured onboard the spacecraft. The ULF Pc5 wave occurrence is observed in both upward and downward FACs. A fairly good correlation is demonstrated between these ULF Pc5 waves and the consecutive injection of magnetosheath low energy protons. The constancy of the observed frequency peak at 1.8 mHz during quite unsteady solar wind pressure conditions could be reconciled with the surface wave mode model. The 3.1 mHz peak location area probably resembles field-line fluctuations with an interesting appearance of poloidal mode oscillation. It is suggested that the 1.3 mHz wave and its harmonic 2.6 mHz represent global compressional oscillations.     

Magnetospheric and ionospheric manifestations of interplanetary foreshocks

The geometry of a typical interplanetary shock front in the vicinity of the Earth’s orbit predicts that the leading edge of the foreshock region comes into contact with the magnetosphere a few hours ahead of geomagnetic sudden impulses (SI). There is reason to believe that the interaction of the magnetosphere with the foreshock leads to magnetic and ionospheric disturbances, which can be detected by ground-based instruments. We searched for specific precursors of SIs in data from the Scandinavian riometer network and in the short period geomagnetic pulsation data from mid-latitude magnetometers. We found that SIs were preceded by the following three features: (1) an increase in riometric absorption, (2) excitation of Pcl magnetic pulsations and (3) a spectral broadening of the Pc3 magnetic pulsations. Our observations may be useful for the study of acceleration processes in the solar wind. These observations are also of potential forecasting interest.     

Interplanetary shocks and sudden impulses during solar maximum (2000) and solar minimum (1995–1996)

In this work a study is performed on the correlation between fast forward interplanetary shock parameters at 1 Astronomical Unit and sudden impulse (SI) amplitudes in the H-component of the geomagnetic field, for periods of solar activity maximum (year 2000) and minimum (year 1995–1996). Solar wind temperature, density and speed, and total magnetic field, were taken to calculate the static pressures (thermal and magnetic) both in the upstream and downstream sides of the shocks. The variations of the solar wind parameters and pressures were then correlated with SI amplitudes. The solar wind speed variations presented good correlations with sudden impulses, with correlation coefficients larger than 0.70 both in solar maximum and solar minimum, whereas the solar wind density presented very low correlation. The parameter better correlated with SI was the square root dynamic pressure variation, showing a larger correlation during solar maximum (r = 0.82) than during solar minimum (r = 0.77). The correlations of SI with square root thermal and magnetic pressure were smaller than with the dynamic pressure, but they also present a good correlation, with r > 0.70 during both solar maximum and minimum. Multiple linear correlation analysis of SI in terms of the three pressure terms have shown that 78% and 85% of the variance in SI during solar maximum and minimum, respectively, are explained by the three pressure variations. Average sudden impulse amplitude was 25 nT during solar maximum and 21 nT during solar minimum, while average square root dynamic pressure variation is 1.20 and 0.86 nPa^1/2 during solar maximum and minimum, respectively. Thus on average, fast forward interplanetary shocks are 33% stronger during solar maximum than during solar minimum, and the magnetospheric SI response has amplitude 20% higher during solar maximum than during solar minimum. A comparison with theoretical predictions (Tsyganenko’s model corrected by Earth’s induced currents) of the coefficient of sudden impulse change with solar wind dynamic pressure variation showed excellent agreement, with values around 17 nT/nPa^1/2.     

球坐标系六片网格下三维定态行星际太阳风模拟

采用二阶MacCormack差分格式, 利用稳态的磁流体(MHD)方程组在球坐标系六片网格下模拟研究了行星际太阳风. 六片网格系统能有效避免极区奇性和网格收敛性. 迭代按径向方向推进求解, 很大程度上减少了计算量, 节约了计算时间. 内边界条件根据太阳与行星际观测确定, 比较测试了5种内边界条件, 模拟给出了1922卡林顿周的背景太阳风结构. 几种内边界条件所得模拟结果与行星际观测基本吻合. 太阳风速度采用McGregor 等的经验公式给出, 磁场由水平电流片(HCCS)模型得到, 密度和温度分别根据动量守恒和气压守恒得到, 研究表明采用这样的边界条件模拟结果最佳.      

The characteristics of sharp (small-scale) boundaries of solar wind plasma and magnetic field structures

We investigate properties of large (>20%) and sharp (<10 min) solar wind ion flux changes using INTERBALL-1 and WIND plasma and magnetic field measurements from 1996 to 1999. These ion flux changes are the boundaries of small-scale and middle-scale solar wind structures. We describe the behavior of the solar wind velocity, temperature and interplanetary magnetic field (IMF) during these sudden flux changes. Many of the largest ion flux changes occur during periods when the solar wind velocity is nearly constant, so these are mainly plasma density changes. The IMF magnitude and direction changes at these events can be either large or small. For about 55% of the ion flux changes, the sum of the thermal and magnetic pressure are in balance across the boundary. In many of the other cases, the thermal pressure change is significantly more than the magnetic pressure change. We also attempted to classify the types of discontinuities observed.