投稿时间:2014-08-11最后修改时间: 2014-12-07 .基金项目：国家自然科学基金(41204120, 41304130, 41474141);中央高校自主科研项目(2042014kf0251, 2042014kf0269)和中国博士后科学基金(2013M542051)联合资助.作者简介：沈晓飞,男,1992年生,2014年毕业于武汉大学获学士学位,现为北京大学研究生/div>

　　作者简介：沈晓飞,男,1992年生,2014年毕业于武汉大学获学士学位,现为北京大学研究生/div通讯作者：倪彬彬,男,1978年生,2005年毕业于武汉大学获博士学位,现为武汉大学空间物理系教授,主要从事磁层物理、波粒相互作用及空间天气效应与建模等方面的研究,E-mail:.cn：日冕物质抛射(Coronal Mass Ejection,简称CME)和共转相互作用区(Corotating Interaction Region,简称CIR)是造成日地空间行星际扰动和地磁扰动的两个主要原因,提供了地球磁暴的主要驱动力,进而显著影响地球空间环境.为深入研究太阳风活动及受其主导影响的地磁活动的时间分布特征,本文对大量太阳风参数及地磁活动指数的数据进行了详细分析.首先,采用由NASA OMNIWeb提供的太阳风参数及地磁活动指数的公开数据,通过自主编写matlab程序对第23太阳活动周期(1996-01-01—2008-12-31)的数据包括行星际磁场指数等主要的地磁指数进行统计分析,建立了包括269个CME事件和456个CIR事件列表的数据库.采用事例分析法和时间序列叠加法分别对两类太阳活动的四个重要太阳风参数(IMFst指数最小值确定了第23太阳活动周期内的355个孤立地磁暴事件,并以st指数最小值为标准将这些磁暴进一步分类为145个弱磁暴、123个中等磁暴、70个强磁暴、12个剧烈磁暴和5个巨大磁暴.最后,采用时间序列叠加法对不同强度磁暴的太阳风参数和地磁指数进行统计分析.统计分析表明,对于CME事件,线性拟合斜率一般为负,这可作为辨别CME和CIR事件的一种有效方法.从平均意义上讲,相较于CIR事件,CME事件有更大的南向IMF.一般情况下,CME事件有更大的可能性驱动极强地磁暴.总体而言,对于不同强度的地磁暴,st指数衰减的速度变快.CME和CIR事件以及其各自驱动的地磁暴事件有着很多不同,因此,需要将CME事件驱动的磁暴及CIR事件驱动的磁暴分开研究.建立CME、CIR事件及地磁暴的数据库以及获取的统计分析结果,将为深入研究地球磁层等离子体片、辐射带及环电流对太阳活动的响应特征提供有利的帮助.A statistical analysis of solar wind parameters and geomagnetic indices for the Solar Cycle 231. School of Physical Science and Technology, Wuhan University, Wuhan 430072, China;

　　太阳风是由源于太阳的等离子体流以及随着等离子体流一起流动的行星际磁场构成的.太阳风经过地球时，地磁场阻挡太阳风，并形成地球弓激波、磁鞘和磁层顶.太阳风和地球磁层的耦合导致在后者区域产生电流、电场，并驱动磁层中大尺度对流，造成磁层顶位形变化，驱使着地球磁层的内部变化，产生磁暴或者亚暴；也影响着磁层内部各电流体系及各区域空间的分布状况，如环电流、辐射带等(苑顺周，2011).

　　地磁活动是受太阳风影响很大的空间天气中的一种重要现象.地磁环境发生变化时，特别是巨大磁暴期间，电离层、高层大气和高能粒子环境都会出现强烈的响应，空间环境变得非常恶劣，会对各种航天器以及宇航员安全构成极大威胁.日冕物质抛射(CME)和共转相互作用区(CIR)是影响地磁活动的 主要因素(Heberet al.，1999；Cane and Richardson，2003;Zhanget al.，2008). 其中CME事件是大尺度的太阳活动，主要发生在太阳活动极大期，持续时间较短，但经常伴有磁暴急始，且强度很大；与此相反的是，CIR事件往往较温和，主要发生在太阳活动极小期，持续时间较长，很少出现磁暴急始，且强度较小，很少出现Dst值小于-100 nT的情况.地磁暴的产生与CME和CIR事件有着密切关系(Turneret al.，2009)；而地磁暴的产生又能对地球磁层的粒子动力学过程以及等离子体波扰动产生显著的影响(Reeveset al.，2003)，并随着时间的长短表现出不同时间尺度的变化特征.因此，深入、细致地了解太阳风以及地磁活动的时变特性对全面、系统地认知地球磁层对太阳活动的响应特性以及太阳风-磁层的粒子、能量耦合过程具有重要的价值.

　　本文通过对第23太阳活动周期(1996—2008)太阳风参数以及地磁指数的详细分析，构造可靠的CME事件、CIR事件以及地球磁暴的数据库，进而细致、定量地研究多个重要的太阳风参数与地磁水平指数的统计与平均特性.2 数据与分析方法2.1 数据说明

　　本文用于研究的太阳风及地磁活动的数据均来源于美国NASA的 OMNIWeb公开数据(主要分析的参量与指数包括行星际磁场Bz分量(IMFBz)、太阳风质子密度(Nsw)、太阳风速度(Vsw)、太阳风动压(Pdyn=2×10－6NswV2sub>

　　sw)、地磁Dst、AE以及Kp指数.采用数据的时间分辨率为1 h.2.2 数据分析方法

　　本文主要采用事例分析法和时间序列叠加法详细分析1996—2008年的三类重要事件：日冕物质抛射、共转相互作用区以及孤立地磁暴.

　　基于事例分析法，对第23太阳活动周期所有CME事件、CIR事件和孤立地磁暴事件进行判别，确定各个事件的起始与结束时间，建立CME、CIR以及磁暴事件的完整列表，绘制每个事件重要太阳风与地磁参数的时间变化图，并分析具体事例参数间的潜在联系.

　　时间序列叠加方法则对某一类事件按照相同标准选择时间零点，晚于时间零点记为正，早于时间零点记为负.例如，将发生于2002年5月23日10时16分选为时间零点，则2002年5月24日10时16分记为+1天，2002年5月22日10时16分记为-1天，以此类推.然后将待分析的参量数据依据相同的时间序列描述，在感兴趣的代表性时间序列长度内对各个参数的变化特征及平均特性进行详细的量化分析，以求获取一般性的时变规律.

　　本文在分析CME及CIR事件时选取事件的起始时间为时间零点，CME及CIR事件的起始时间依据行星际激波的到达时间或行星际磁场南向开始时间进行选取，结束时间依据太阳风参数及地磁指数恢复平静的时间进行选取(Chenet al.，2012).其优势在于可以获取CME及CIR事件随时间变化的一般性规律，为研究和辨识CME及CIR事件提供依据.在分析地磁暴时选取Dst指数达到最小值的时刻为时间零点，其优势在于便于对地磁暴进行分类，获取不同类型地磁暴的发生规律，也便于研究地磁暴期间Dst指数与其他地磁指数的关系.3 分析结果3.1 CME与CIR事件举例及其磁暴

　　中左图给出一个CME事件及其对应的磁暴现象.CME事件如图中两条竖虚线年第269天(9月26日)1142UT，持续时间约为1.5天.而磁暴发生的时间更长，其结束时间一般选择为

　　st指数基本恢复到暴前水平的时刻.在这期间，太阳风参数与地磁指数有以下的变化：(1)初始有一个较大的南向IMF

　　分量，约为-17.7 nT，发生在第268天0200UT.对应着IMF

　　最终在0 nT附近小幅度变化.太阳风速度降低缓慢，太阳风质子密度与动压降低速度相对更快，说明这是一个密度增强主导的动压脉冲事件；(2)对应着南向IMF

　　st指数最小值为-207 nT，这是一个剧烈磁暴，发生在第268天0900UT，对应着较小南向IMF

　　指数最大可达8.3，发生在第268天0600UT.磁暴恢复相持续大约1天半的时间；AE峰值后出现较多极值幅度为几百nT的振荡，表明亚暴的持续发生；

　　Examples of CME and CIR events and associated geomagnetic storms

　　图 1右图给出一个CIR事件及其对应的磁暴现象.如图中两条竖虚线年第84天(即3月24日)0000UT，持续时间约达4天.在这期间，太阳风参数与地磁指数表现出 以下的变化：(1)南向IMFBz值更小，仅为-5.8 nT，发生在第80天0700UT，但是By与Bz分量均表现出更频繁的振荡.太阳风速度从397 km·s-1增加到653 km·s-1，于第81天2300UT达到最大值，随后缓慢降低.太阳风质子密度与动压变化也非常缓慢，最大质子密度为15.5 cm-3，而太阳风动压不超过6 nPa；(2)Dst指数最小值仅为-66 nT，发生在第81天2200UT，属于一个中等磁暴，而且没有磁暴急始现象.AE指数最大值为844 nT，同样发生在第80天1900UT；AE指数振荡比较明显，说明亚 暴事件比较频繁.Kp指数峰值仅为5.7，并在Kp=3 附近保持很长时间，对应着持续一天多的磁暴恢复相.引发地磁暴强度较弱和Kp指数恢复平静时间较长都是CIR事件激发磁暴的典型特征(Borovsky and Denton，2006；Chenet al.，2012).3.2 CME与CIR事件完整列表

　　基于3.1节展示的对CME与CIR事件的判断，同时依据Borovsky等总结出的CME与CIR事件特征的主要差别(Borovsky and Denton，2006)，本文分析了1996—2008年太阳风参数以及地磁指数的OMNIWeb数据库，共判定269个CME事件和456个CIR事件及其起始与结束时间.下面对这些CME与CIR事件进行分析.3.3 CME与CIR事件的太阳风参数统计分析3.3.1 时间序列叠加法分析

　　本文采用时间序列叠加方法对269个CME事件和456个CIR事件(少数CME事件及CIR事件因跨越年份或数据缺失严重等原因未纳入分析)对应的IMFBz、太阳风速度(Vsw)、太阳风质子密度(Nsw)、太阳风动压(Pdyn)、Dst指数、AE指数及Kp指数7个参量进行分析.选取CME或CIR事件的起始时间为时间零点，对于CME事件，时间范围选取为起始时间前一天及后两天，共三天数据；对于CIR事件，时间范围选取为起始时间前一天及后四天，共五天数据.

　　图 2左右两列图分别是1996—2008年所有CME事件和CIR事件的太阳风参数(IMFBz分量、太阳风速度、太阳风质子密度、太阳风动压)采用时间序列叠加分析得到的统计结果.其中黑色“+” 号表示每个CME或CIR事件的数据点，红线表示采用时间序列叠加分析后得到的算术平均时变曲线.从统计意义上讲，CME事件比CIR事件表现出 跨越度更大的IMFBz、太阳风速度与太阳风动压数值区间，而CIR事件对应的太阳风质子密度数值区 间更大.比如，有的CME事件对应着>

　　1000 km·s-1太阳风速 度，而CIR事件的太阳风速度普遍-1； 另外，太阳风动压峰值大于30 nPa的CME事件甚多于相应的CIR事件.CME事件的太阳风质子密度均不超过60 cm-3，而CIR事件的质子密度峰值可达100 cm-3.时间序列叠加后的平均结果显示，CME与CIR事件对应的IMFBz分量相差不大.虽然CME与CIR事件对应的平均太阳风速度在时间零点都表现出显著的增大，但前者的变化更迅速，且在零点后几个小时内随时间明显降低，后者的变化相对平缓，且在零点后1~2天时间内缓慢增大.另外，CME事件对应的平均动压变化比CIR事件明显，但CIR事件对应的平均质子密度比CME事件略强.图 2Fig. 2图 21996—2008年所有CME及CIR事件对应的四个主要太阳风参数

　　对于CME事件，起始时间后各参数均很快达到最值；而CIR事件，各参数达到最值的时间往往较长.结合表 1，可以获知，从平均意义上讲，CME事件有更大的南向IMFBz分量、太阳风动压Pdyn、AE指数、Kp指数与更小的Dstmin，说明CME事件往往比CIR事件更加剧烈，激发的磁暴与亚暴也更强(Dentonet al.，2006)；而CIR事件的太阳风速度和太阳风质子密度值更大，则是因为CIR事件常与太阳风质子高速流有关.3.3.2Nsw/Pdyn斜率的线给出的时间序列叠加分析结果，本文对1996—2008年所有CME及CIR事件的

　　Nsw与Pdyn参数进行了相关分析：用Nsw除以Pdyn，并对该曲线随时间的变化趋势进行线性拟合获取直线斜率，拟合区域为时间零点后整个CME或CIR事件的持续时间.图 4给出第23太阳活动周期所有CME和CIR事件Nsw/Pdyn线性拟合斜率的分布，其中红色“+”号代表CME事件，蓝色“*”代表CIR事件.分析发现，CME事件和CIR事件发生期间的N

　　sub

　　曲线线性拟合斜率的显著差异体现了CME、CIR发生期间太阳风速度不同的变化特征.CME期间，太阳风速度先在较短时间内大幅上升，然后缓慢下降，下降持续时间比上升时间长，对应到

　　拟合曲线斜率为正.CIR期间，由于高速质子流的影响，太阳风速度持续增强并维持较长时间，使得

　　结果表明，太阳风速度的整体变化趋势是辨别CME和CIR事件的有效、合理参量.

　　Scatter plot of the linear fitting slope of

　　corresponding to all CME and CIR events of 1996—2008

　　确定地磁暴的发生必须满足两个基本条件：(1)具有较大的地磁扰动强度，即

　　st指数最小值小于-30 nT；(2)有明显的起始相、主相和恢复相.其中，孤立地磁暴具有重要的研究价值.本文着重于研究孤立地磁暴事件，其判断标准是：(1)

　　st时变曲线比较平滑或可用平滑曲线)起始时间和结束时间对应的Dst值接近于磁静状态.通过对1996—2008年的Dst指数曲线图逐一进行目测分析，最终确定355个孤立地磁暴事件，平均每年约27个.地磁暴强度分类标准有很多种，本文采用Loewe等给出的利用Dstmin的分类标准(Loewe and Prlss，1997

　　st衰减的速度变快，并且恢复相的持续时间变长.但总体而言，对于不同强度的地磁暴，

　　st指数的变化呈现出一定的相似性.对于其他六个参数，也有类似结果.另外，需要说明的是：由于巨大磁暴样本非常少，叠加平均后仍有较大波动.

　　Superposed epoch analysis results of all isolated geomagnetic storms of 1996—2008

　　CME事件和CIR事件是造成日地空间行星际扰动和地磁扰动的两个主要原因，提供了地球磁暴的主要驱动力，进而显著影响了地球空间环境.在深入研究1996—2008年(第23太阳活动周期)太阳风活动及受其影响的地磁活动的时间分布特征后，本文得到以下主要结论：

　　(1)第23太阳活动周期共发生269个CME事件与456个CIR事件.从统计意义上讲，CME事件比CIR事件表现出跨越范围更大的IMF

　　、太阳风速度与太阳风动压数值区间，而CIR事件对应的太阳风质子密度数值区间更大.

　　dyn曲线线性拟合斜率k值为正，99.8%的CIR事件对应的Nsw/Pdyn曲线线性拟合斜率k值为负，体现了CME、CIR发生期间太阳风速度不同的变化特征，因此太阳风速度的整体变化趋势是辨别CME与CIR事件的有效、合理参量.(4)1996—2008年共有355个孤立地磁暴.随着地磁暴强度的增强，Dst衰减的速度变快，并且恢复相的持续时间变长.但总体而言，对于不同强度的地磁暴，Dst指数的变化呈现出一定的相似性.本文的研究专注于第23太阳活动周期，对第24太阳活动周期的研究将在后续的工作中展开.本文工作集中统计分析了期间的CME事件、CIR事件及孤立地磁暴事件特征，对于太阳风参数及地磁指数表现出来的统计性特征或差异性并未深入探讨.下一步的工作也将包括通过AE-

　　Dst双指数间的比较考察磁层亚暴与磁暴之间的关联性，定量分析驱动不同强度磁暴的太阳风条件及参数特征，以及将CME事件和CIR事件引起的磁暴分立开分别进行更细致、更本质的分析.另外，把太阳风以及地磁活动的详细分析与地球磁层的粒子动力学过程以及波场时空分布数据系统地结合起来，将大大有助于更深入地认识地球磁层包括等离子体片、环电流、辐射带对太阳活动的时变响应特性，从而更有效、合理地建立太阳风-磁层-电离层的全球耦合模型.致谢感谢UCLA大学的Tung-Shin Tsu博士以及NASA的Lan Jian博士给予的帮助与讨论.