地球物理学报  2016, Vol. 59 Issue (3): 767-777   PDF    
太阳高能粒子(SEP)传播数值模拟中的太阳风背景场研究
魏稳稳1,2, 沈芳1, 左平兵1, 秦刚1, 杨子才1,2    
1. 中国科学院国家空间科学中心空间天气学国家重点实验室, 北京 100190;
2. 中国科学院大学, 北京 100049
摘要: 太阳高能粒子(SEP)事件是一类重要的空间天气灾害性事件,如能准确预报SEP事件,人们便可以采取必要的防护措施,保障卫星、星载设备以及航天员的安全,尽可能地降低经济损失.因此,其数值预报研究在空间天气预报研究中占有很重要的地位.SEP事件中的高能粒子在不同的时间尺度内被耀斑过程或者CME驱动的激波加速,并且在被扰动后的行星际太阳风中传输,这些过程都紧紧依赖于太阳风背景场.因此获取更加接近物理真实的太阳风背景场是模拟SEP事件的重要部分,也是提高SEP物理模式的关键因素之一.我们目前的工作基于张明等发展的SEP在行星际空间传播的模型,尝试将Parker太阳风速度解及WIND飞船观测的磁场实时数据融入模型中,研究不同的太阳风速度以及真实磁场分布对SEP在行星际空间中传播的影响.通过求解聚焦传输方程,我们的模拟结果表明:(1)快太阳风条件下,绝热冷却效应项发挥了更大的作用,使粒子能量衰减的更快,而慢太阳风对粒子的通量变化没有显著影响;(2)加入观测的磁场数据时,粒子的全向通量剖面发生了比较明显的变化,具体表现在:通量峰值推迟到达、出现多峰结构、各向异性也发生一些改变.分析表明真实磁场的极性对粒子在行星际空间中传播有着重要的影响.
关键词: 太阳高能粒子     磁场     太阳风     行星际输运    
Effects of the solar wind background field on the numerical simulation of the Solar Energetic Particle(SEP) transportation
WEI Wen-Wen1,2, SHEN Fang1, ZUO Ping-Bing1, QIN Gang1, YANG Zi-Cai1,2    
1. State Key Laboratory of Space Weather, National Space Science Center, Chinese Academy of Sciences, Beijing 100190, China;
2. University of Chinese Academy of Sciences, Beijing 100049, China
Abstract: Solar energetic particles(SEPs) pose one of the most serious hazards to spacecraft systems and constrain human activities in space. Thus, it is of importance to forecast SEP events. Several theories and numerical models are applied to simulate SEP events. Each model makes some assumptions to simplify the complex acceleration and transportation processes within such events. In general, SEP will interact with ambient solar wind and background magnetic field during transportation. It is recognized that interplanetary transport effects must be taken into account at any analysis of SEP propagation. In the previous simulation, it always assumed Parker magnetic field and fixed solar wind speed as the input parameters. However, these assumptions are too simple when compared with the real conditions. In order to get better results, it is necessary to use more accurate background conditions. Recently, we change the fixed solar wind speed into spatial-dependent speed profile based on Parker's theory, and replace the Parker magnetic field with another Parker-like magnetic field based on in situ data at 1 AU. By solving the focused transport equation with simulation of time-backward stochastic processes method, our results show that:(1) Under fast solar wind speed assumption, it is clear that the omnidirectional flux decreases faster than that for the situation with slow solar wind speed in the decay phase. We suggest that it is due to the adiabatic cooling effect. Fast solar wind speed has a significant effect on the adiabatic cooling, which leads the SEPs to lose energy more quickly during transportation. However, slow solar wind speed has less impact on the time profiles of SEP flux and anisotropy. We also compare the time profiles of SEP event observed at different observatories and energies, the results remain the same as previous;(2) When applying in situ data of magnetic field observed by WIND during different Carrington Rotations, the omnidirectional flux time profiles vary greatly, and the main results are as followings:the peak flux appears to be delayed, multi-peak occur, anisotropy also has some differences.We think it results from the magnetic field polarity, which affects the pitch angle, and, furthermore, modulates the momentum. The characteristics are similar in solar minimum and solar maximum, while the peaks seem to be more when solar activity is active. We conclude that the real magnetic field polarity may exert a significant influence during the propagation of SEP. In the future, we will try to use the real-time background conditions which obtain from MHD models in our simulations, in order to make a thorough study of the SEP propagation.
Key words: Solar energetic particle     Magnetic field     Solar wind     Interplanetary transport    
1 引言

太阳高能粒子(Solar Energetic Particles,简称SEP)事件是能量粒子的通量突然增强的事件,美国国家海洋和大气管理局(NOAA)将SEP 事件定义为能量大于10 MeV 的粒子在连续15 min以上的时间内数目超过10 pfu(cm-2·s-1·sr-1)的爆发性事件(Rodríguez-Gasén,2011).SEP 事件主要包括三种类型(Kallenrode,2003):与太阳耀斑爆发相关联的脉冲型事件,与日冕物质抛射驱动的激波相关联的缓变型事件,以及同时具有缓变型和脉冲型事件特征的混合型事件.SEP 事件会对卫星和宇航员构成严重的威胁,如能准确预报SEP 事件,人们便可以对卫星和星载设备采取必要的防护措施,同时也可以为航天器的故障分析和航天员的安全防护提供科学依据,尽可能地降低损失.然而SEP 事件产生和发展的物理机制非常复杂,目前还没有被完全解释清楚,预报SEP 事件的能力也不尽如人意(Lario,2005).因此,正确理解和恰当地动力学描述这些机制,对推进和提高空间天气预报能力是非常重要的(魏稳稳等,2015).

SEP 事件中的高能粒子在不同的时间尺度内被耀斑过程或者CME 驱动的激波加速,并且在被扰动后的行星际磁场和太阳风中传输.近几十年里,很多研究者探究了太阳风背景场对高能粒子在行星际中传播的影响.Lario等(20032008)基于Ulysses 以及ACE卫星的观测,对太阳活动高年里较低能段的高能粒子事件进行了分析,认为行星际磁力线结构会影响高能粒子的经纬度分布.并重点分析了大尺度的行星际结构对2004年9月SEP事件的影响,认为其强度和各向异性时间剖面是由一个共转的低密度、低速、低质子β值的太阳风流决定.大尺度太阳风结构对高能粒子传播的影响决定不同位置的宇宙飞船观测到的SEP事件的性质,不同于低β值区域有利于高能粒子的散射,被压缩的磁场区域会阻碍高能粒子的自由流动.磁镜效应与散射过程能够对高能粒子进行有效的约束,这些结构相对于观测者和粒子源项的位置决定SEP事件的特征(Lario et al.,2008).另外,Shen等(2008)的研究也表明磁云这种复杂结构也能对高能粒子进行有效的约束并显著影响高能粒子的通量.

近年来一些研究者采用聚焦传输方程模拟粒子在行星际空间中的传播过程,聚焦传输方程包括了许多重要的粒子传输机制:沿磁力线的流动、太阳风对流、投掷角扩散、磁场聚焦效应、垂直扩散、绝热冷却效应等.模拟结果(Ruffolo,1995;Qin et al.,2006; Zhang et al.,2009)发现,在聚焦传播方程中引入绝热冷却效应项,将增加粒子强度的衰减率.Lario等(Lario et al.,1997)研究了太阳风对能量为0.1~20 MeV的质子在行星际空间传播的影响,发现太阳风对流能使粒子更早地到达观测者,绝热冷却效应则导致粒子通量减小.

SEP 事件中的高能粒子在被扰动后的行星际磁场和太阳风中传输,这一过程紧紧依赖于太阳风和背景场(Barouch and Burlaga,1976;Lario et al.,1997;Mal and raki et al.,2007).以往粒子传输模型中大多采用理想的Parker螺旋磁场和恒定的太阳风速度,背景场的设置过于简单(Heras et al.,1995; Li et al.,2003; Zhang et al.,2009; Wang et al.,2012).因此,建立更加接近物理真实的背景场来模拟粒子在三维行星际空间中的传播是提高SEP 事件模拟结果不可或缺的部分,本文将围绕背景场模拟展开讨论,主要研究太阳风速度以及真实磁场分布对SEP在行星际空间中传播的影响.

2 模型简介

SEP在日球层中的传播是SEP研究中的一个核心问题,SEP的空间分布等都与之密切相关.研究高能粒子在行星际空间中的传输时需要求解聚焦传输方程.聚焦传输方程的本质是Fokker-Plank方程,用于描述粒子分布函数f在相空间的演化.有些学者(Kallenrode and Wibberenz,1997; Ruffolo,1995; Lario et al.,1997)通过有限差分格式对聚焦传输方程进行离散求解.Qin等(2006)Zhang等(2009)开发并发展了一种计算SEP 在三维日球层磁场中传播的模式,该模式采用后向随机微分方法解聚焦传输方程.用于模拟研究观测者处于行星际空间不同的经度、纬度以及径向距离时,SEP 通量和各向异性随时间的演化过程.

模型采用三维的聚焦传输方程来研究粒子的传播过程,其形式如下:

其中f为回旋平均分布函数,κ是垂直扩散系数,μ为投掷角的余弦值,v是粒子的速度,b 是沿着磁场方向的单位矢量,V sw是太阳风的速度,Dμμ为投掷角扩散系数,LB=(b · ∇ lnB)-1为磁场聚焦尺度,p是粒子的动量.

方程中的绝热冷却项的表达式如下:

投掷角余弦μ随时间演化的表达式为

我们根据Qin等和Zhang等(Zhang,1999; Qin et al.,2006; Zhang et al.,2009)采用的后向随机微分方法来求解聚焦传输方程.方程(1)可以改写为五个一维的后向随机微分方程:

求解随机微分方程采用后向时间的技巧,即让虚拟粒子从指定的地点r0出发,以一定的初始动量p,从时刻t=t0开始沿后向时间运动,直到粒子退出边界.我们假设所有的高能粒子都在太阳附近注入,因此,粒子分布函数在源区满足边界条件:

其中,a(Φ,θ)是一个关于日面经度Φ和日面纬度θ的函数,用来描述 SEP源区强度的空间分布,p是源区粒子动量,TcTl分别表示源区粒子注入剖面的上升和衰减时间尺度.我们以下的模拟事例中 均假设源区里面的粒子分布均匀一致,即a(Φ,θ)=1.

方程的解即为所有数值模拟过程在退出点上的平均值:

其中,x eμepe表示虚拟粒子退出边界时的参量,te是粒子退出边界的时刻,N是虚拟粒子的总量.

粒子全向通量(Omnidirectional flux)的计算公式为,各向异性(Anisotropy)指数A的定义式为,模型中用到的参数见表 1(Zhang et al.,2009),我们的模拟过程不考虑垂直扩散对粒子传播的影响,将垂直扩散系数设置为0,即κ=0.

表 1 模型中参数设定 Table 1 Model parameters used in the calculations
3 模拟结果及分析 3.1 太阳风速度对粒子传播过程的影响

根据物理性质和起源区域的不同,太阳风一般分为快太阳风(v≥550 km·s-1)和慢太阳风(v≤450 km·s-1)(Wang et al.,2000).我们分别计算了快、慢太阳风对粒子通量的影响,快太阳风对不同能量质子的传播过程以及对不同观测位置的通量剖面的影响.并且统计分析了快、慢太阳风条件下统计到的粒子的动量分布特征.模型中磁场采用Parker螺旋磁场.

为了研究发展的太阳风速度对粒子在行星际空间中传播的影响,我们将模型中沿径向恒定不变的太阳风速度 V sw=Vsw e r(Vsw=400 km·s-1)变为Parker太阳风速度解(由(6)式解出),Parker太阳风速度解的具体数值实现方法详见文献Parker(1958).从(6)式中可以看出,Parker太阳风速度解不再沿径向不变,而是与日心距离r相关.

其中Vsw满足如下方程:

图 1给出了太阳风速度接近慢太阳风的模拟结果.图 1a给出了该事例中太阳风速度随日心距离r 的变化曲线,太阳风速度在1AU处达到383 km·s-1. 以下图形中标注为慢太阳风的都是基于图 1a中的Parker慢太阳风速度解.图 1b给出了不同太阳风速度条件下100 MeV能量的质子在1AU赤道处的 全向通量对比图.黑色实线代表固定速度400 km·s-1 的模拟结果,红色虚线代表固定速度383 km·s-1的模拟结果,绿色虚线代表Parker太阳风速度解接近慢太阳风速度下的模拟结果.从图 1b中可以看出这三条通量曲线非常接近,说明慢太阳风相对固定的太阳风速度而言,对粒子的通量变化影响不大.

图 1 太阳风速度接近慢太阳风的模拟结果
(a) 太阳风速度随日心距离的变化; (b) 不同太阳风条件下能量为100 MeV的质子的全向通量对比.
Fig. 1 Simulation results under slow solar wind speed assumption
(a) Slow solar wind speed profile; (b) Comparison of 100 MeV protons flux between different solar wind speed.

图 2给出了太阳风速度接近快太阳风的模拟结果.图 2a中Parker太阳风速度解在1AU处达到644 km·s-1.以下图形中标注为快太阳风的都是基于图 2a中的Parker快太阳风速度解.图 2b给出了不同太阳风速度条件下100 MeV能量的质子在1AU赤道处的全向通量对比图.黑色实线代表固定速度400 km·s-1的模拟结果,红色虚线代表固定 速度644 km·s-1的模拟结果,绿色虚线代表Parker 太阳风速度解接近快太阳风速度下的模拟结果.从图 2b中可以看出,固定太阳风速度变大时,将增加粒子通量曲线的衰减率,而Parker快太阳风条件下,粒子通量曲线的衰减率进一步增加.通过分析我们认为这是由于快太阳风条件下,绝热冷却效应项发挥了更大的作用.在Parker太阳风速度解的模型中,绝热冷却项dp/dt的表达式如下:

图 2 太阳风速度接近快太阳风的模拟结果
(a)太阳风速度随日心距离的变化; (b) 不同太阳风条件下能量为100 MeV的质子的全向通量对比.
Fig. 2 Simulation results under fast solar wind speed assumption
(a) Fast solar wind speed profile; (b) Comparison of 100 MeV protons flux between different solar wind speed.

第一项:第二项:

从式(7)可以看出,粒子的动量随时间的函数是衰减的,绝热冷却效应会使得观测到的粒子比它在源区释放时的能量要小.相比于沿径向恒定不变的太阳风速度,在Parker太阳风速度解的条件下,绝热冷却项dp/dt右侧表达式中的方括号里多了两项,这两项在快慢太阳风条件下均为正值.在快太阳风条件下,不仅方括号里的数值增大,Vsw也超过400 km·s-1,使得表达式右侧显著变小,导致绝热冷却效应相对速度固定为400 km·s-1的太阳风明显增强.而在慢太阳风条件下,虽然方括号里的数值增大,但是Vsw小于400 km·s-1,使得表达式右侧变化不大,导致绝热冷却效应相对而言没有显著改变.另外,SEP的谱指数是小于零的,所以对于某一特定能段的粒子而言,随着时间推移,观测到的粒子通量会越来越少.因此,在快太阳风条件下,粒子的通量如图 2所示衰减的更快,粒子全向通量在后期可达到一个量级的差异.而慢太阳风条件下则如图 1所示没有显著变化.

由于在SEP的初始相时期粒子密度的梯度大,扩散效应比绝热冷却效应更重要,所以在SEP的初始相可以忽略绝热效应,在SEP的衰减时期,粒子密度空间梯度小,绝热冷却效应将起重要的作用(涂传诒等,1988).在Parker太阳风速度解的条件下,对于某一固定的径向距离r,太阳风速度增大时.相应的Vsw/r和dVsw/dr也随之增大,第一项和第二项与太阳风速度的大小呈线性关系,因此,太阳风速度越高,会造成SEP向外对流的越多,绝热能量损失也随之增大,绝热效应也就越显著.

我们又对快、慢太阳风速度条件下统计到的粒子的动量进行了统计分析.图 3给出了观测到的粒子的动量分布图.黑色实线,红色虚线分别代表Parker太阳风速度解接近快太阳风速度、慢太阳风速度时的统计结果.我们将观测到的粒子动量的最大值和最小值(pmax & pmin)等分为10个区间Δp= 0.1(pmax -pmin),然后统计每个动量区间内的粒子数占总的粒子数的百分比.表 2为不同太阳风条件下的统计结果,我们设置的目标粒子的动量为p=0.445.其中目标粒子的动量定义为观测者处最终观测到的粒子动量,通过能量计算可以得到,该能量是人为设定的,是为了便于统计分析某一个具体能量段内的粒子的传播特征、分布特征等.

图 3 快慢太阳风条件下统计到的粒子动量分布Fig. 3 Momentum distributions of observed particles

图 3表 2的结果可以看出:

表 2 不同太阳风条件下的统计结果 Table 2 Statistical results under different solar wind speed conditions

(1)从所统计到粒子的最低动量来看,它们都比目标粒子对应的动量高,说明粒子在运动过程中都有能量衰减;

(2)从所统计到粒子的最高动量来看,快太阳风条件下粒子能量的最大值更大,说明快太阳风条件下粒子能量的衰减可以更加显著;

(3)投放相同的粒子数(投放1×106个粒子)时,快太阳风条件下统计到的粒子数小于慢太阳风条件下的粒子数,也验证了上述结论;

(4)从粒子分布来看,越靠近目标粒子能量段的粒子数越多,且呈现出单调递减的趋势,说明大部分统计到的粒子是能量衰减较少的粒子.

此外,我们又模拟了快太阳风对不同能量的粒子和不同观测位置处的通量剖面的影响.图 4a4b分别给出了不同太阳风速度条件下100 MeV能量的质子在观测者位于0.5AU赤道处以及2AU、纬度80°处的全向通量对比.图 4c给出了不同太阳 风速度条件下10 MeV能量的质子在1AU赤道处的全向通量对比图.黑色实线代表固定速度为400 km·s-1 的模拟结果,红色虚线代表Parker太阳风速度解接近快太阳风速度下的模拟结果.从结果中可以得出,快太阳风对不同能量段的粒子在不同观测位置处的通量剖面有着相似的影响.

图 4 快太阳风条件下观测者处于不同位置时能量为100 MeV的质子的全向通量对比((a)和(b))以及能量为10 MeV的质子的全向通量对比(c)Fig. 4 Comparisons of 100 MeV protons flux at different locations (panels (a) and (b)) and energies (panel (c))
3.2 用实时的磁场观测数据

以往粒子传输模型中大多采用Parker螺旋磁场,Parker螺旋磁场是一个稳态模型,在太阳风是定态球对称的假定下,根据麦克斯韦方程组、磁冻结效应,以及在源表面(Rs,θs,Φs)处径向的边界条件 可以推导出Parker螺旋磁场模型,具体形式为(Parker,1958)

其中Ω为太阳自转频率,Vsw为常数太阳风速度.对于距离太阳很远的地方r≥Rs,根据磁通量守恒定律,模型中用到的是简化的Parker磁场,如表 1中所示: ,其中B0=3.54 nT是由1AU处的磁场估算得到.

然而,这种简化了的Parker磁场忽略了磁场大小和极性的变化,但是这些因素对粒子在行星际空间中的传输也可能发挥重要作用.因此,为了构造更加真实的背景磁场,我们使用了WIND卫星实时的磁场观测数据来进行背景磁场的设置.Florens等(2007)推导出了可以根据卫星实时观测数据求解磁场三分量(Brr,Bθ,BΦ)的类Parker磁场解,如下式(8)所示,我们根据该方法将WIND卫星的观测数据拓展到全空间.

假设测试粒子从太阳源表面(Rs,Φs)离开,到达空间任意位置(r,Φ),并且太阳风在每一个卡林顿周期内是不变的,则可以根据宇宙飞船的观测数据估算Rs处太阳风的参量.参数r,Rs,Φ,Φs,观测时间t,区间[t0t0+T]的开始时间t0由以下公式计算得到:

模型中使用的是WIND飞船一分钟精度的磁场数据,我们分别计算了处于太阳活动低年以及太阳活动高年的不同卡林顿自转周期(简称CR)内SEP的全向通量和各向异性的时间剖面.

太阳活动低年的模拟结果如图 5所示.其中,图 5b是根据WIND飞船在2007年5月3日至2007年5月29日处于CR2056和CR2057内的观测数据得到的黄道面内的磁场图.为了更好的可视化效果,我们将任意日心距离r处的磁场进行了标准化处理,具体标准化处理方法为:将同一条Parker螺线上,任意日心距离r处的磁场,均赋为该条螺线在1AU处的磁场数值,然后进行可视化作图.太阳坐标为(X,Y)=(rcosΦ,rsinΦ)=(0,0),地球处在黑色圆圈处,并按照顺时针方向进行背景磁场的设置.两条黑线的夹角为统计到的粒子所在源区的大致分布范围.图 5a中红色虚线为在此CR期间内根据类Parker磁场模拟得到的粒子全向通量剖面,黑色实线为原Parker磁场下的模拟结果,图 5c为相应的各向异性剖面.相比之下,我们发现在类Parker磁场下模拟得到的粒子通量峰值会到达的晚一些,全向通量剖面会出现多个峰值并且在衰减相下降的更快,而且一阶各向异性也有明显的不同.

图 5 (a) 不同行星际磁场下,对比100 MeV质子的通量; (c) 各向异性; (b) 根据WIND飞船数据得到的CR2056-CR2057内磁场Br分量图.太阳坐标为(X,Y)=(rcosΦ,rsinΦ)=(0,0),黑色圆圈为地球 Fig. 5 Distribution of 100 MeV proton omnidirectional flux (a) and anisotropy (c) at different IMF conditions during CR 2056-2057. (b) is distribution of Br based on Wind spacecraft data. The Sun is at (X,Y)=(rcosΦ, rsinΦ)=(0,0) and the Earth is at the solid black circle

由式子(10)—(12)可知,当加入磁场的观测数据后,cosβ和sinβ的大小不会改变,但是它们的符号会随着磁场的极性而发生相应的改变,从而会影响投掷角μ值的变化,并最终对动量p的值产生影响.因此,我们认为粒子通量峰值推迟到达,多峰结构以及各向异性的差异主要是由于加入的真实磁场 具有不同的极性所致.

针对通量峰值推迟到达,我们比较了Parker螺旋磁场和类Parker螺旋磁场下最大峰值附近的粒子初始能量的大小,发现Parker磁场下粒子初始能量的平均值大约为0.463,类Parker磁场下粒子初始能量的平均值大约为0.628,而我们统计的目标粒子能量为0.445.由于能量越高,衰减到目标能量所需要的时间越长,所以在类Parker磁场中由于磁场极性对粒子能量的调制,它们到达观测者所需的时间会更久,从而导致峰值推迟到达.

针对出现的多峰结构,我们在图 6中给出了第二个较大峰值到达前后不同时间段内统计到的粒子初始动量p的分布范围.黑色实线、红色虚线、绿色虚线分别代表第二个峰值到达前两天,第二个峰值期间以及第二个峰值过后两天内的粒子动量分布.从图 6可以看出,不同时间段内动量的分布存在明显差异.在第二个峰出现的时间段内(红色虚线),p=1附近出现了另一个峰值.通过对通量(为p的函数)的计算,我们发现这个峰值对该期间的粒子通量的贡献起决定性作用,并且使得该通量超过了前后两天的通量大小.通常情况下随着时间的推移,在衰减相不同时间段内统计到的大部分粒子的动量会不断增大,即类似图 6中黑色实线的峰值会不断后移,从而使得全向通量不断减小.然而,同样,由于磁场极性对粒子能量的调制,使得粒子动量出现了额外的峰值并发挥了主要作用,从而导致出现了多峰结构.

图 6 CR2056-CR2057模拟得到的不同时间段内到达观测者的粒子初始动量p的分布范围Fig. 6 Initial momentum distributions of observed particles in different intervals during CR2056-CR2057

不同于Parker磁场下的粒子各向异性,类Parker磁场下的各向异性在上升相位时小于零,而在衰减相期间,各向异性在模拟误差范围内会在零附近上下浮动.由于观测点的磁场极性保持不变,则根据各向异性指数:,以及μ= p·b/p可知,当A为正(负)时,统计到的主要是朝磁场方向同(反)向运动的粒子.因此,在上升相位期间,朝磁场相反方向运动(μ=-1)的粒子占主导.而在衰减相期间,由于磁场极性的影响以及部分粒子投掷角方向的改变,会使得各向异性在误差范围内出现波动.

太阳活动高年的模拟结果,我们以CR 1969为例展示在了图 7中.类似图 5的格式,图 7a给出了模拟得到的粒子全向通量剖面,图 7c为各向异性剖面,图 7b给出了WIND飞船在该周期内所观测到的磁场.其中两条黑色线同样代表统计到的粒子所在源区的大致范围.类似太阳活动小年的情况,从图 7a图 7c中可以看出,粒子通量峰值也比原Parker 磁场到达的晚一些,粒子的全向通量剖面出现了多 峰结构,而各向异性在误差范围内并没有太大的差异.

图 7 (a) 不同行星际磁场下,对比100 MeV质子的通量; (c) 各向异性;(b) 根据WIND飞船数据得到的CR1969的磁场Br分量Fig. 7 Distribution of 100 MeV proton omnidirectional flux (a) and anisotropy (c) at different IMF conditions during CR 1969. (b) is distribution of Br based on Wind spacecraft data

同样地,我们计算得到类Parker磁场下最大峰值附近所统计到的粒子初始能量的平均值大约为 0.562,要大于Parker磁场下粒子能量的平均值0.463,因此粒子在类Parker螺线磁场里运动到达观测者所需的时间更久,导致峰值推迟到达.但是,相比于太阳活动小年粒子的平均初始能量为0.628,该周期内统计到的粒子初始能量要更小一些,所以相比之下粒子的全向通量剖面峰值到达的也更早一些.至于更小的粒子初始能量,可能与太阳活动高年更复杂的磁场极性变化对粒子能量的调制有关.

此外,与太阳活动小年相比,在太阳活动高年里模拟得到的粒子的全向通量剖面会出现更多的峰值.图 8给出了第二、三个峰值到达时以及它们到达之前一段内观测到的粒子初始动量p的分布范围.从图 8可以看出,在这两个峰值到达之前,如黑线所示,粒子初始动量大约以0.8为峰值,而且粒子的动量分布范围更加广泛.在红线所示的第二个峰值到达期间,粒子的初始动量更加集中地分布在以0.6左右为峰值的一个小范围内,因此根据通量计算公式可知该段时间内粒子的全向通量更大.而在绿线所示的第三个峰值到达期间,粒子的初始动量仍以0.6左右为峰值,但是它们占该时间段内所统计到总粒子数的比例更小,而且来源更加广泛,这就是为什么第三个峰值要比第二个峰值低的原因.类似于太阳活动小年的情况,在太阳活动高年里磁场对粒子能量的调制会导致多峰结构,而这种更多峰值结构的出现,是由于磁场的极性出现了更多的变化,对粒子能量的调制更加复杂所致.

图 8 CR1969模拟得到的不同时间段内到达观测者的粒子初始动量p的分布范围Fig. 8 Initial momentum distributions of observed particles in different intervals during CR1969

图 7c中可以看出,类Parker磁场下的各向异性在上升相位时大于零,说明朝磁场相同方向运动(μ=1)的粒子占主导.同样地,在衰减相期间,由于磁场极性的影响以及粒子投掷角方向的改变,会使得各向异性在误差范围内出现一些波动.

从模拟结果可以看出,太阳活动高年和低年的结果虽然都包含通量峰值推迟到达、出现多峰结构、各向异性发生改变等特点,但也存在差异.由于太阳活动高年磁场比较复杂,导致通量剖面出现更多的峰,各向异性剖面也存在差异.但二者的模拟结果都充分说明了真实磁场极性对粒子在行星际空间中的传播有着重要的影响.

4 结论

本文通过求解聚焦传输方程研究了太阳风速度和背景磁场对粒子在行星际空间中传播的影响,模拟结果表明:

(1)将固定的太阳风速度改为变化的Parker太阳风速度解,从模拟结果可以看出慢太阳风对粒子的通量变化没有显著影响,而当太阳风速度加快时,粒子通量剖面在后期下降的更快.这是因为快太阳风条件下,绝热冷却效应更加显著,从而使粒子的能量衰减的更快.

(2)相对于原Parker磁场模型,加入观测的磁场数据时粒子的全向通量剖面发生了比较明显的变化:通量峰值推迟到达、出现多峰结构,各向异性也发生一些改变.通过分析,我们认为加入观测的磁场数据时,磁场有了极性,粒子的投掷角会随着磁场的极性而发生相应改变,进而也会对粒子的能量进行调制,从而导致模拟结果出现了上述的变化.无论在太阳活动低年还是高年,粒子全向通量剖面和各向异性剖面都伴随着这些改变,而太阳活动高年更加复杂的变化可能是由于磁场的极性更加复杂所致.

我们的模拟结果说明太阳风速度以及背景磁场对粒子在行星际空间中的传播有着重要的影响.以往粒子传输模型中大多采用恒定的太阳风速度和Parker螺旋磁场,背景场的设置太过理想,这可能导致无法准确模拟太阳高能粒子在行星际空间中的演化过程.因此,设置更加真实的背景场对提高我们的模拟结果有着重要的意义.未来的研究工作中我 们将尝试加入由MHD模拟(Feng et al.,20102012; Shen et al.,20112014)得到的实时变化的太阳风速度和三维磁场分量,建立更加完整、更加接近物理真实的太阳风背景场,以期得到更真实的粒子在三维行星际空间的传播物理过程,同时可以研究日冕物质抛射相互作用等对高能粒子事件形成和传播的影响(Gopalswamy et al.,2002; Shen et al.,2008; Li et al.,2012; Shen et al.,2013).

参考文献
[1] Barouch E, Burlaga L F. 1976. Three-dimensional interplanetary stream magnetism and energetic particle motion.J. Geophys. Res., 81(13):2103-2110, doi:10.1029/JA081i013p02103.
[2] Feng X S, Yang L P, Xiang C Q, et al. 2010. Three-dimensional solar WIND modeling from the sun to earth by a SIP-CESE MHD model with a six-component grid.Astrophys. J., 723(1):300-319, doi:10.1088/0004-637X/723/1/300.
[3] Feng X S, Jiang C W, Xiang C Q, et al. 2012. A data-driven model for the global coronal evolution.Astrophys. J., 758(1), doi:10.1088/0004-637X/758/1/62.
[4] Florens M S L, Cairns I H, Knock S A,et al. 2007. Data-driven solar wind model and prediction of type Ⅱ bursts. Geophys. Res. Lett., 34(4):L04104, doi:10.1029/2006GL028522.
[5] Gopalswamy N S, Yashiro S, Michałek G, et al. 2002. Interacting coronal mass ejections and solar energetic particles.Astrophys. J., 572(1):L103-L107, doi:10.1086/341601.
[6] Heras A M, Sanahuja B, Lario D, et al. 1995. Three low-energy particle events:modeling the influence of the parent interplanetary shock. Astrophys. J., 445(1):497-508, doi:10.1086/175714.
[7] Kallenrode M B. 2003. Current views on impulsive and gradual solar energetic particle events.Journal of Physics G:Nuclear and Particle Physics, 29(5):965-981, doi:10.1088/0954-3899/29/5/316.
[8] Kallenrode M B,Wibberenz G. 1997. Propagation of particles injected from interplanetary shocks:A black box model and its consequences for acceleration theory and data interpretation.J. Geophys. Res., 102(A10):22311-22334, doi:10.1029/97JA01677.
[9] Lario D, Sanahuja B, Heras A M. 1997. Modeling the interplanetary propagation of 0.1-20 MeV shock-accelerated protons. I:Effects of the adiabatic deceleration and solar wind convection. Adv. Space Res., 20(1):115-120, doi:10.1016/S0273-1177(97)00492-4.
[10] Lario D, Roelof E C, Decker R B, et al. 2003. Solar maximum low-energy particle observations at heliographic latitudes above 75 degrees.Advances in Space Research, 32(4):579-584, doi:10.1016/S0273-1177(03)00339-9.
[11] Lario D. 2005.Advances in modeling gradual solar energetic particle events.Advances in Space Research, 36(12):2279-2288, doi:10.1016/j.asr.2005.07.081.
[12] Lario D, Decker R B, Malandraki O E, et al. 2008. Influence of large-scale interplanetary structures on energetic particle propagation:September 2004 event at Ulysses and ACE. J. Geophys. Res., 113(A3):A03105, doi:10.1029/2007JA012721.
[13] Li G, Zank G P, Rice W K M. 2003. Energetic particle acceleration and transport at coronal mass ejection-driven shocks. J. Geophys. Res., 108(A2):1082, doi:10.1029/2002JA009666.
[14] Li G, Moore R, Mewaldt R A, et al. 2012.A Twin-CME scenario for ground level enhancement events.Space Science Reviews, 171(1-4):141-160, doi:10.1007/s11214-011-9823-7.
[15] Malandraki O E, Marsden R G, Tranquille C,et al. 2007. Energetic particle observations by Ulysses during the declining phase of solar cycle 23. J. Geophys. Res., 112(A6):A06111, doi:10.1029/2006JA011876.
[16] Parker E N. 1958.Dynamics of the interplanetary gas and magnetic fields.Astrophys. J., 128:664-676, doi:10.1086/146579.
[17] Qin G, Zhang M, Dwyer J R. 2006. Effect of adiabatic cooling on the fitted parallel mean free path of solar energetic particles.J. Geophys. Res., 111(A8):A08101, doi:10.1029/2005JA011512.
[18] Rodríguez-Gasén R. 2011. Modelling SEP events:latitudinal and longitudinal dependence of the injection rate of shock-accelerated protons and their flux profiles[Ph. D. thesis]. Barcelona:Universitatde Barcelona.
[19] Ruffolo D. 1995.Effect of adiabatic deceleration on the focused transport of solar cosmic rays.Astrophys. J., 442(2):861-874, doi:10.1086/175489.
[20] Shen C L, Wang Y M, Ye P Z, et al. 2008. Enhancement of solar energetic particles during a shock-magnetic cloud interacting complex structure.Solar Physics, 252(2):409-418, doi:10.1007/s11207-008-9268-7.
[21] Shen C L, Li G, Kong X, et al. 2013. Compound twin coronal mass ejections in the 2012 May 17 GLE event. Astrophys. J., 763(2):114-121, doi:10.1088/0004-637X/763/2/114.
[22] Shen F, Feng X S, Wu S T, et al. 2011. Three-dimensional MHD simulation of the evolution of the April 2000 CME event and its induced shocks using a magnetized plasma blob model.J. Geophys. Res., 116(A4):A04102, doi:10.1029/2010JA015809.
[23] Shen F, Shen C L, Zhang J, et al. 2014. Evolution of the 12 July 2012 CME from the Sun to the Earth:Data-constrained three-dimensional MHD simulations.J. Geophys. Res., 119(9):7128-7141, doi:10.1002/2014JA020365.
[24] Wang Y, Qin G, Zhang M. 2012. Effects of perpendicular diffusion on energetic particles accelerated by the interplanetary coronal mass ejection shock.Astrophys. J., 752(1):37, doi:10.1088/0004-637X/752/1/37.
[25] Wang Y M, Sheeley N R Jr, Socker D G, et al. 2000.The dynamical nature of coronal streamers.J. Geophys. Res., 105(A11):25133-25142, doi:10.1029/2000JA000149.
[26] Wei W W, Shen F, Zuo P B. 2015.Research progress on the solar energetic particle model based on magnetohydrodynamic simulation.Progress in Astronomy(in Chinese), 33(1):1-26, doi:10.3969/j.issn.1000-8349.2015.01.01.
[27] Zhang M. 1999. A Markov stochastic process theory of cosmic-ray modulation.Astrophys. J., 513(1):409-420, doi:10.1086/306857.
[28] Zhang M, Qin G, Rassoul H. 2009. Propagation of solar energetic particles in three-dimensional interplanetary magnetic fields.Astrophys. J., 692(1):109-132, doi:10.1088/0004-637X/692/1/109.
[29] 涂传诒等. 1988. 日地空间物理学.北京:科学出版社.
[30] 魏稳稳, 沈芳, 左平兵. 2015. 基于磁流体力学模拟的太阳高能粒子物理模式研究进展. 天文学进展, 33(1):1-26, doi:10.3969/j.issn.1000-8349.2015.01.01.