2. 中国科学院近地空间环境重点实验室, 合肥 230026;
3. 安徽蒙城地球物理国家野外科学观测研究站, 合肥 230026
2. Key Laboratory of Geospace Environment, Chinese Academy of Sciences, Hefei 230026, China;
3. Mengcheng National Geophysical Observatory, Hefei 230026, China
在中低纬台站经常观测到地磁场水平分量大幅度减小,并且大约持续数天,这种现象叫做地磁暴(Gonzalez et al., 1994).地磁暴是一类灾害性空间 天气事件,它能引起磁层、电离层和中高层大气等的剧烈变化,从而对航天器、输电系统等产生各种危害.一般地,地磁暴的强度可用AE指数、Kp指数以及Dst指数等来描述.根据Dst指数最小值的不同,地磁暴可分为弱磁暴(Dstmin≤-30 nT)、中等磁暴(Dstmin≤-50 nT)以及强磁暴(Dstmin≤-100 nT)等.
一般认为,地磁暴是具有南向分量的行星际磁场与地球磁层发生相互作用而形成(Gonzalez et al., 1994).因而,可能伴随行星际磁场南向分量的行星际日冕物质抛射、激波以及共转相互作用区等被认为是地磁暴的主要行星际源(Gonzalez et al., 1994,1999; Zhang et al., 2003; Wang et al., 2002b; Xue et al., 2005; Gonzalez et al., 2007; Zhang et al., 2007b;Echer et al., 2008; Gonzalez et al., 2011; Echer et al., 2013).近年特别是进入23太阳活动周以来,随着观测技术的进步,人们开始对太阳活动和近地空间的环境进行长期的监测.利用这些观测数据,人们对地磁暴特别是强地磁暴的行星际源进行了深入、广泛的讨论.根据太阳和行星际局地观测,Zhang等(2007a,2007b)统计分析了1996—2005年间Dstmin≤-100 nT的强地磁暴的行星际和太阳源,发现大部分(87%)的强地磁暴主要由日冕物质抛射(Coronal Mass Ejection,CME)及其相关行星际结构引起,而共转相互作用区(Corotating Interaction Regions,CIR)则引起了约13%的强地磁暴.对1996—2008年的中等地磁暴(-100 nT<Dstmin≤-50 nT)的行星际源分析,Echer等(2013)发现共转相互作用区和太阳风高速流(HSSs)是中等地磁暴的主要驱动源.进一步,人们根据23太阳周完整的行星际和地磁暴观测数据,对太阳活动不同强度时的地磁暴强度及起因进行了研究,发现在太阳低年,共转相互作用区是引起地磁暴的主要起因,而在太阳高年,行星际日冕物质抛射(Interplanetary Coronal Mass Ejection,ICME)及其相关结构则是引起地磁暴的主要因素(Zhang et al., 2007a; Echer et al., 2013).
现有的研究工作已对上一太阳周的地磁暴事件进行了较完整的分析和讨论. 2008年开始,太阳进入了一个活动极弱的时期.2008年和2009年中,分别出现了266天和262天的无黑子日(李可军等,2010).同时,人们也开始注意到23太阳周后地磁活动性急剧减弱,并根据不同的观测资料对地磁活动性减弱的原因进行了讨论(Richardson,2013; Gopalswamy et al., 2014),发现可能的原因是在该段时间内日冕物质抛射的数目和速度(Richardson,2013)、太阳风的磁场强度(Gopalswamy et al., 2014)等均明显下降.同时,2007—2012年间地磁暴及其行星际源的完整列表仍没有建立,对这段时期内地磁暴的完整统计分析研究仍缺乏.本文将根据WIND卫星和ACE卫星的局地行星际磁场以及太阳风等离子体参数等的观测数据,对2007—2012年间的所有中等以上(Dstmin≤-50 nT)地磁暴及其行星际源进行统计分析.第2节将给出地磁暴事例的选取及发生率随太阳周的变化.第3节将介绍地磁暴行星际源的认定方法、典型事例及统计结果.最后,将在第4节中给出结论和讨论. 2 地磁暴事件的选取及其发生率年变化
根据国际数据中心(http://wdc.kugi.kyoto-u.ac.jp/dstdir/)提供的Dst指数数据,我们挑选出 了2007年至2012年所有的Dst指数峰值小于-50 nT(Dstmin≤-50 nT)的中等以上地磁暴事件.图 1给出了这段时间内Dst指数随时间的变化.图中的菱形符号对应于我们所选取的事件及其峰值大小.在此期间,国际数据中心一共记录到51次Dstmin≤-50 nT的中等以上地磁暴,其中有9次Dstmin≤-100 nT的强磁暴,没有Dstmin≤-200 nT的特大地磁暴.2007年至2012年中,最强的地磁暴发生于2011年10月25日,其Dst指数峰值为-132 nT.
为了比较,选取了上一活动周类似的区间:1996年至2001年.图 2给出了太阳黑子数随时间的变化,图中阴影区域(b)为本文所研究的时间段,阴影区域(a)为上次类似时间段(1996年至2001年).可以看到,在本文研究的区间中(图 2中阴影区域(b)),黑子数远小于上一类似时期,这表明在该段时期内太阳活动性显著减弱.同时,对地磁暴产生数的对比发现,2007年至2012年间的地磁暴数目显著减小.在1996年至2001年间,发生了48次Dstmin≤-100 nT的强磁暴事件(Zhang et al., 2007b),其中Dstmin≤-200 nT的地磁暴发生了10次,最大地磁暴的Dst指数峰值为-387 nT(2001年3月31日).对比结果表明,2007年开始的这个太阳活动极小期及24周的太阳活动上升期,地 磁活动显著减弱,这与以前所得结果一致(Gopalswamy et al., 2014).
图 3给出了这段时间地磁暴发生数随年度的变化.从图 3可以看到,地磁暴的发生数呈现明显的太阳活动周变化,其发生数与太阳黑子数(图 3中的红色“﹡”符号)的演化规律基本一致.2007年和2008年,每年各发生了5次中等以上地磁暴事件.而在太阳活动极低的2009年,只有一个中等以上地磁暴发生.随着太阳活动的增强,从2010年开始,地磁暴的数目快速增加.特别的,2012年共记录到中等以上地磁暴19次.图 3中蓝色部分给出了Dstmin≤-100 nT的强地磁暴数目随年份的变化.可以看到,在太阳活动较弱的2007年至2010年,没有发生Dstmin≤-100 nT的强地磁暴.随着太阳黑子数的迅速增加,2011年和2012年分别发生了3次和6 次Dstmin≤-100 nT的强磁暴.这表明地磁活动的强弱可能受太阳活动强弱的直接影响.
我们利用WIND卫星和ACE卫星的局地行星际磁场及太阳风等离子体观测数据认证了2007年至2012年间的51次地磁暴事件的行星际源.在Zhang等(2007a)以及Echer等(2013)的工作中,地磁暴对应的行星际源被详细地分为行星际日冕物质抛射、具有前导激波的行星际日冕物质抛射、激波鞘区、共转相互作用区、高速太阳风等.在本文中,我们把地磁暴的行星际源简单分为与日冕物质抛射相关的现象、与共转相互作用区相关现象及其他等三类.表 1给出了这些事件的列表和行星际源认证结果.同时,我们的在线列表(http://space.ustc.edu.cn/dreams/geostorms2007/(保持更新))还进一步给出了这些事件的相关图像.下面,将对每一类事件的典型事例进行详细说明.
行星际日冕物质抛射是行星际局地观测所探测到的特殊结构,一般被认为是日冕物质抛射的行星际观测表现.行星际日冕物质抛射的可能观测特征包括增强的磁场、高能粒子通量的降低、增强的磁场、不同寻常的电离态以及双向电子流等(Wang et al., 2002b; Jian et al., 2006a).快速的日冕物质抛射可能驱动激波,进而与太阳风相互作用形成由激波、激波鞘区和行星际日冕物质抛射所组成的结构.图 4给出了一个由与日冕物质抛射相关结构引起的地磁暴的典型事例.图中灰色阴影区域为可能的行 星际日冕物质抛射区域.该行星际日冕物质抛射于2010年5月28日19 ∶ 50UT到达地球,2010年5月29日12 ∶ 10UT结束.在该区域中,磁场强度明显增强且平滑旋转、质子温度和等离子β低且存在明显的双向电子流,这是一个典型的行星际日冕物质抛射观测.同时,该行星际日冕物质抛射驱动激波于2010年5月28日01 ∶ 50UT(如图中红色虚线所示)被WIND卫星观测到.从Dst指数随时间的演化可以看出,该激波的到来引起了Dst指数的突然增加,即地磁暴急始的开始.随后,伴随着鞘区和ICME所携带的南向行星际磁场的到来,Dst指数缓慢下降,开始了磁暴的主相,并于2010年5月29日14 ∶ 00到达其峰值-85 nT.
一般认为,地球附近观测到的行星际激波由日冕物质抛射驱动.但是,一些局地探测到的行星际激波可能并没有观测到随后的行星际日冕物质抛射结构(Gopalswamy et al., 2009).这并不表明这类激波不是由日冕物质抛射驱动,可能的原因有:(1)激波宽度远大于日冕物质抛射本体,这导致我们的局地观测只能探测到激波而探测不到日冕物质抛射本体;(2)日冕物质抛射在传播过程中发生偏转(Shen et al., 2011a; Wang et al., 2004; Gopalswamy et al., 2009; Wang et al., 2014),而导致其本体没有被局地卫星探测到.激波到达地球压缩行星际磁场 也可能引起地磁暴,图 5给出一个只由激波引起的地 磁暴事例:2012年1月24日事例.Shen等(2013,2014)报道了这次激波事件及其驱动的日冕物质抛射,他们发现这个激波由爆发于日面(15°W,40°N)的一次全晕状日冕物质抛射引起.但由于日冕物质抛射本体的宽度较小或可能的偏转,WIND卫星未观测到日冕物质抛射本体对应的行星际日冕物质抛射,而只观测到了这次具有较宽宽度的激波.由于该 类事件中的激波也是由日冕物质抛射所驱动,因而 在本文中我们仍认为该类事件为与日冕物质抛射相关的地磁暴.
基于卫星的局地观测数据,多个日冕物质抛射在行星际空间中形成复杂结构也已被广泛地报道(Wang et al., 2003a,2002a,2003b,2003c; Lugaz and Farrugia, 2014),在本文中,我们没有单独区分这种复杂结构,而简单归为与日冕物质抛射相关结构类.因而,与日冕物质抛射相关结构包括单个行星际日冕物质抛射、行星际日冕物质抛射及其驱动激波结构以及日冕物质抛射相互作用复杂结构等三小类,这三类均与太阳爆发活动相关. 3.2 与共转相互作用区相关的地磁暴
共转相互作用区是高速太阳风和低速太阳风相互作用形成的区域(Jian et al., 2006b).图 6给出 了一次由共转相互作用区引起地磁暴的典型事例:2010 年5月2日地磁暴事件.在该次事件中,共转相互作用区于2010年5月2日03 ∶ 00UT至16 ∶ 00UT 间被局地卫星观测到.从图 6可以看到,该区间为低速太阳风和高速太阳风的过渡区间,在该区间内太阳风速度缓慢上升.而高速太阳风与低速太阳风的相互作用导致了该区域内的磁场、等离子密度和温度的增强.本次共转相互作用区所携带的南向行星际磁场使得Dst指数从约2010年5月2日09 ∶ 00UT 开始下降,并于2010年5月2日18 ∶ 00UT达到其峰值-66 nT.
在本文中,其他不与日冕物质抛射或者共转相互作用区相关的行星际源均归为其他类.图 7给出了一个事例:2010年12月28日地磁暴事件.从图中可以看出,该次地磁暴开始时间约为2010年12 月28日11 ∶ 00UT,并于2010年12月28日17 ∶ 00UT 到达其峰值-50 nT.在该段时间附近,行星际观测无法找到明显的结构.从图 7中可以看到,在该段时间中太阳风中存在一些小尺度的磁场扰动.这些扰动中的南向行星际磁场是形成这次地磁暴事件的主因.
表 2给出了对所有事件的行星际源的统计结果.从表中可以看到,在2007年至2012年这段时间内,大部分(65%)的地磁暴由与日冕物质抛射相关的行星际结构引起,而由共转相互作用区引起的地磁暴占总地磁暴的31%.特别的,所有Dstmin≤-100 nT的强磁暴均由与日冕物质抛射相关的行星际结构引起.而在1996年至2001年间的,共转相互作用区引起了3次Dstmin≤-100 nT的强地磁暴 事件(Zhang et al., 2007a).这表明,在这段太阳活 动较弱的时期内,共转相互作用区的地磁效应也被削弱.同时,日冕物质抛射相关结构也没能引起Dstmin≤-200 nT的地磁暴,这也表明日冕物质抛射本身的地磁效应相对也显著减弱.图 8给出了不同行星际源引起地磁暴数目随年份的变化.从图中可看到,在太阳下降期的2007年和2008年,大部分地磁暴均由共转相互作用区引起.特别的,在太阳极低的2009年,唯一地磁暴由共转相互作用区引起.而随着太阳活动的增强,2010年开始与日冕物质抛 射相关的行星际结构则成了地磁暴的主要行星际源.
本文建立了最近一段太阳活动极低及24周太阳活动上升期(2007年至2012年)所有Dstmin≤-50 nT的中等以上地磁暴事件列表,并根据行星际观测资料认证了这些地磁暴的行星际源.在此基础上,我们讨论了这些地磁暴及其行星际源的特征和随太阳周的变化规律,并与上一周类似太阳活动 时期的已有研究结果进行了初步的对比.主要结论包括:
(1)2007年至2012年间,共发生地磁暴51次,其中Dstmin≤-100 nT的强地磁暴9次,没有Dstmin≤-200 nT的超强地磁暴产生,远小于上一类似太阳活动期间(1996年至2001年)的地磁暴数,这与已有的最近太阳活动低年和上升期的地磁活动水平弱的结果一致(Gopalswamy et al., 2014; Kilpua et al., 2014).进一步,我们发现地磁暴发生的次数与太阳黑子数正相关,且Dstmin≤-100 nT的强地磁暴只在上升期和高年间发生,这与以前的结果一致.
(2)在该段时间内,约65%的地磁暴由与日冕物质抛射相关结构引起,31%的地磁暴由共转相互作用区引起,这与以前的观测结果类似.同时,在该段时间内所有Dstmin≤-100 nT的强地磁暴均由与日冕物质抛射相关结构引起.与上一个类似太阳活 动区间的比较发现,与日冕物质抛射相关结构没有 引起Dstmin≤-200 nT的超强地磁暴,共转相互作用区也没能引起Dstmin≤-100 nT的强地磁暴.该结果表明,在该段太阳活动极弱的时期内,共转相互作用区和日冕物质抛射的地磁效应弱,没能产生强的地磁暴.对于共转相互作用区,一个可能的解释是在该段时间内背景太阳风速度、行星际磁场特别是其南向分量强度等均显著降低(Jian et al., 2011; Kilpua et al., 2014),这可能导致共转相互作用区的磁场强度等均减小,故而地磁效应也相应减弱.最近,Gopalswamy等(2014)发现引起24周前4年的地磁暴的行星际日冕物质抛射的磁场强度、太阳风速度等均小于23周相似时间段,这可能是日冕物质抛射相关结构地磁效应减弱的直接原因.
(3)在太阳下降期(2007年,2008年)和低年(2009年),共转相互作用区是引起地磁暴的主要行星际源.而在上升期和高年,与日冕物质抛射相关的行星际结构是引起地磁暴的主要因素,这与以前的观测结果一致.
需要指出的是,在这段时间内,有约18%(9/51)的地磁暴由持续时间小于12 h的行星际磁结构引起.在本文中,由于这些结构具有行星际日冕物质抛射的特征或者包含在共转相互作用区中,因而我们仍归为与日冕物质抛射相关或者与共转相互作用区相关的地磁暴类.但是,最近研究也提出太阳风中的小结构也可能来源于冕流区域或者在太阳风中产生(Yu et al., 2013).因而,利用多点、多种类的观测资料,对这些引起地磁暴的小尺度规则结构进行详细分析,将为我们更进一步了解可能的地磁暴行星际源及这些小结构的形成、传播具有重要的意义.
日冕物质抛射相互作用对地磁效应的影响也已被研究(Wang et al., 2003a,2002a,2003b,2003c; Lugaz et al., 2005; Xue et al., 2005; Shen et al., 2011b).在我们分析的事件中可能存在相互作用的事件为7个,约占所有与日冕物质抛射相关事件的21%,这一比例略小于Zhang 等(2007b)的结果.这可能是因为在该段时间中太阳活动较弱,因而连续爆发日冕物质抛射进而发生相互作用的可能性降低.同时,存在相互作用的由日冕物质抛射事件引起的地磁暴也没有显著增强.需要说明的是,现在的日冕物质抛射相互作用只通过局地数据认证,这可能存在一定的偏差.因而,结合日球层成像观测数据对真实存在日冕物质抛射相互作用事件的地磁效应的分析也需要深入的进行.
[1] | "Echer E, Gonzalez W D, Tsurutani B T. 2008. Interplanetary conditions leading to superintense geomagnetic storms (Dst≤-250 nT) during solar cycle 23. Geophysical Research Letters, 35(6): L06S03, doi: 10.1029/2007GL031755. |
[2] |
Echer E, Tsurutani B T, Gonzalez W D. 2013. Interplanetary origins of moderate (-100 nT |
[3] | Gonzalez W D, Joselyn J A, Kamide Y, et al. 1994. What is a geomagnetic storm? Journal of Geophysical Research, 99(A4): 5771-5792. |
[4] | Gonzalez W D, Tsurutani B T, Clúa de Gonzalez A L. 1999. Interplanetary origin of geomagnetic storms. Space Science Reviews, 88(3-4): 529-562. |
[5] | Gonzalez W D, Echer E, Clúa Gonzalez A L, et al. 2007. Interplanetary origin of intense geomagnetic storms (Dst<-100 nT) during solar cycle 23. Geophysical Research Letters, 34(6): L06101, doi: 10.1029/2006GL028879. |
[6] | Gonzalez W D, Echer E, Tsurutani B T, et al. 2011. Interplanetary origin of intense, superintense and extreme geomagnetic storms. Space Science Reviews, 158(1): 69-89. |
[7] | Gopalswamy N, Mäkelä P, Xie H, et al. 2009. CME interactions with coronal holes and their interplanetary consequences. Journal of Geophysical Research, 114(A3): A00A22, doi: 10.1029/2008JA013686. |
[8] | Gopalswamy N, Akiyama S, Yashiro S, et al. 2014. Anomalous expansion of coronal mass ejections during solar cycle 24 and its space weather implications. Geophysical Research Letters, 41(8): 2673-2680. |
[9] | Jian L, Russell C T, Luhmann J G, et al. 2006a. Properties of interplanetary coronal mass ejections at one AU during 1995—2004. Solar Physics, 239(1-2): 393-436. |
[10] | Jian L, Russell C T, Luhmann J G, et al. 2006b. Properties of stream interactions at one AU during 1995—2004. Solar Physics, 239(1-2): 337-392. |
[11] | Jian L K, Russell C T, Luhmann J G. 2011. Comparing solar minimum 23/24 with historical solar wind records at 1 AU. Solar Physics, 274(1-2): 321-344. |
[12] | Kilpua E K J, Luhmann J G, Jian L, et al. 2014. Why have geomagnetic storms been so weak during the recent solar minimum and the rising phase of cycle 24? Journal of Atmospheric and Solar-Terrestrial Physics, 107: 12-19, doi: 10.1016/j.jastp.2013.11.001. |
[13] | Li K J, Feng W, Liang H F. 2010. The abnormal 24th solar cycle-The first complete solar cycle of the new millennium. Scientia Sinica Phys., Mech. & Astron. (in Chinese), 40(10): 1293-1301. |
[14] | Lugaz N, Farrugia C J. 2014. A new class of complex ejecta resulting from the interaction of two CMEs and its expected geoeffectiveness. Geophysical Research Letters, 41(3): 769-776. |
[15] | Lugaz N, Manchester IV W B, Gombosi T I. 2005. Numerical simulation of the interaction of two coronal mass ejections from sun to earth. Astrophysical Journal, 634(1): 651-662. |
[16] | Richardson I G. 2013. Geomagnetic activity during the rising phase of solar cycle 24. Journal of Space Weather and Space Climate, 3(A8), doi: 10.1051/swsc/2013031. |
[17] | Shen C L, Wang Y M, Gui B, et al. 2011a. Kinematic evolution of a slow CME in corona viewed by STEREO-B on 8 October 2007. Solar Physics, 269(2): 389-400. |
[18] | Shen C L, Wang Y M, Pan Z H, et al. 2013. Full halo coronal mass ejections: Do we need to correct the projection effect in terms of velocity? Journal of Geophysical Research: Space Physics, 118(11): 6858-6865, doi: 10.1002/2013JA018872. |
[19] | Shen C L, Wang Y M, Pan Z H, et al. 2014. Full-halo coronal mass ejections: arrival at the earth. Journal of Geophysical Research: Space Physics, 119(7): 5107-5116, doi: 10.1002/2014JA020001. |
[20] | Shen F, Feng X S, Wang Y M, et al. 2011b. Three-dimensional MHD simulation of two coronal mass ejections' propagation and interaction using a successive magnetized plasma blobs model. Journal of Geophysical Research, 116(A9), doi: 10.1029/2011JA016584. |
[21] | Wang Y M, Shen C L, Wang S, et al. 2004. Deflection of coronal mass ejection in the interplanetary medium. Solar Physics, 222(2): 329-343. |
[22] | Wang Y M, Wang B Y, Shen C L, et al. 2014. Deflected propagation of a coronal mass ejection from the corona to interplanetary space. Journal of Geophysical Research: Space Physics, 119(7): 5117-5132, doi: 10.1002/2013JA-19537. |
[23] | Wang Y M, Wang S, Ye P Z. 2002a. Multiple magnetic clouds in interplanetary space. Solar Physics, 211(1-2): 333-344. |
[24] | Wang Y M, Ye P Z, Wang S, et al. 2002b. A statistical study on the geoeffectiveness of Earth-directed coronal mass ejections from March 1997 to December 2000. Journal of Geophysical Research, 107(A11): SSH 2-1-SSH 2-9. |
[25] | Wang Y M, Ye P Z, Wang S. 2003a. Multiple magnetic clouds: Several examples during March-April 2001. Journal of Geophysical Research, 108(A10): 1370, doi: 10.1029/2003JA009850. |
[26] | Wang Y M, Ye P Z, Wang S, et al. 2003b. Theoretical analysis on the geoeffectiveness of a shock overtaking a preceding magnetic cloud. Solar Physics, 216(1-2): 295-310. |
[27] | Wang Y M, Ye P Z, Wang S, et al. 2003c. An interplanetary cause of large geomagnetic storms: Fast forward shock overtaking preceding magnetic cloud. Geophysical Research Letters, 30(13): 31-33. |
[28] | Xue X H, Wang Y M, Ye P Z, et al. 2005. Analysis on the interplanetary causes of the great magnetic storms in solar maximum (2000—2001). Planetary and Space Science, 53(4): 443-457. |
[29] | Yu W, Farrugia C J, Lugaz N, et al. 2013. A statistical analysis of properties of small transients in the solar wind 2007—2009: STEREO and Wind observations. Journal of Geophysical Research: Space Physics, 119(2): 689-708, doi: 10.1002/2013JA019115. |
[30] | Zhang J, Dere K P, Howard R A, et al. 2003. Identification of solar sources of major geomagnetic storms between 1996 and 2000. Astrophysical Journal, 582(1): 520-533. |
[31] | Zhang J, Richardson I G, Webb D F, et al. 2007a. Correction to ""Solar and interplanetary sources of major geomagnetic storms (Dst≤-100 nT) during 1996—2005"". Journal of Geophysical Research, 112(A12): A12103, doi: 10.1029/2007JA012891. |
[32] | Zhang J, Richardson I G, Webb D F, et al. 2007b. Solar and interplanetary sources of major geomagnetic storms (Dst≤-100 nT) during 1996—2005. Journal of Geophysical Research, 112(A10): A10102, doi: 10.1029/2007JA012321." |
[33] | 李可军, 冯雯, 梁红飞. 2010. 异常的第24太阳活动周——新千年的第一个完整的太阳活动周. 中国科学: 物理学 力学 天文学, 40(10): 1293-1301. |