行星际日冕物质抛射引起福布斯下降的一维随机微分模拟

倪素兰; 顾斌; 韩智伊

doi:10.7498/aps.66.139601

摘要

福布斯下降（Forbush decrease，FD）是银河宇宙线（galactic cosmic rays，GCRs）受短期剧烈太阳活动调制的重要现象之一.本文设GCRs进入由行星际日冕物质抛射（interplanetary coronal mass ejection，ICME）及其前沿激波共同形成的扰动区时，其径向扩散系数rr受抑制变为（r）rr（0 （r） 1），且抑制强度与粒子位置处的太阳风等离子体速度正相关.对任意时刻的扰动区，抑制系数（r）在激波处最小为（rsh）），并按指数规律增大，在ICME尾部归一.CME爆发时，（rsh））取全局最小值m.在扰动区向日球层外传播的过程中，（rsh））逐步恢复为1.在此基础上，根据GOES和ACE卫星观测确定模型参数，用一维随机微分方程描述GCRs在日球层内的传播，并采用倒向随机方法模拟了一个由独立Halo ICME调制GCRs引起的2005年5月15日FD事件.计算所得地面中子通量的主相、恢复相及其在CME到达地球前的增加过程，均与Oulu中子探测器观测结果一致.

关键词:

Abstract

Forbush decrease (FD) event is one of the most important short-term modulations of galactic cosmic rays (GCRs) caused by intense solar activities such as interplanetary coronal mass ejection (ICME). The modulation mechanisms of GCRs by the disturbed interplanetary magnetic fields (IMF) of ICME and the accompanying forward interplanetary shock (IP) are not clear yet. In this work, we present a one-dimensional dynamic model of the GCR barrier driven by ICME. In our model, the time dependent radial diffusion coefficient rr of GCRs is depressed to be (r)rr (0 (r) 1) as they run into the disturbed IMF. The scale factor (r) is inversely proportional to the local solar wind speed away from the Sun. Within the disturbed area at any time, (r) increases exponentially from the local minimum (rsh)) at the IP front to 1 at the end of the ICME tail. In addition, (rsh)) switches gradually from its global minimum m at the bursting of the CME to 1 as the shock moving toward the outer boundary of the heliosphere. The geometrical and dynamic parameters of the ICME and IP are derived from the observations of GOES and ACE satellites. Based on the stochastic transport theory, the one-dimensional backward stochastic differential equation (SDE) method is adopted to simulate the transport of GCRs modulated by single halo ICME. The evolution of the neutron flux at the ground is calculated according to the recently reported proton-neutron yield function. As an example, the FD event on 15 May 2005, caused by the CME event bursting on 13 May 2005, is studied and simulated. The results show that the calculated neutron flux evolution, including not only the main and recovery phases, but also the pre-enhancement before the arriving of the CME at the Earth, is consistent with the observation of Oulu neutron monitor. According to the trajectories of GCRs, it can be found that, the per-enhancement of the neutron flux is a result of the scattering by the forward IP passing 1 AU. Before the IP reaches the switch cutoff Rc, GCRs are evidently confined in the sheath between the IP and CME. After that, the GCRs will stay for longer time in the magnetic cloud of the ICME as a result of the damping of IP strength. The parameterzed one-dimensional GCRs modulation model and the SDE method, as have been confirmed by the neutron monitor observation on the Earth, can be used further to calculate and predict the GCRs fluxes of other places, such as the Mars, in the heliosphere.

Keywords:

interplanetary coronal mass ejection /
Forbush decrease /
backward stochastic differential method /
neutron flux

作者及机构信息

1.
南京信息工程大学物理系, 南京 210044;

2.
南京信息工程大学空间天气研究所, 南京 210044

通信作者: 顾斌, gubin@nuist.edu.cn

Authors and contacts

1.
Department of Physics, Nanjing University of Information Science and Technology, Nanjing 210044, China;

2.
Institute of Space Weather, Nanjing University of Information Science and Technology, Nanjing 210044, China

Corresponding author: Gu Bin, gubin@nuist.edu.cn

参考文献

[1]	Rossi B 1964 Cosmic Rays (New York: McGraw-Hill) pp110
[2]	Blasi P 2013 Astron. Astrophys. Rev. 21 70
[3]	Bothmer V, Daglis I A 2007 Space Weather- Physics and Effects (Berlin: Springer) pp103130
[4]	Le G M 2002 Ph. D. Dissertation (Beijing: Chinese Academy of Sciences) (in Chinese) [乐贵明2002 博士学位论文(北京: 中国科学院空间科学与应用研究中心)]
[5]	Potgieter M S 2013 Living Rev. Sol. Phys. 10 3
[6]	Guo W J, Zhu B Y 1990 Chin. J. Spac. Sci. 10 247 (in Chinese) [郭维吉, 朱邦耀1990 空间科学学报10 247]
[7]	Kharayat H, Prasad L, Mathpal R, Garia S, Bhatt B 2016 Sol. Phys. 291 603
[8]	Cane H V 2000 Space Sci. Rev. 93 55
[9]	Belov A, Abunin A, Abunina M, Eroshenko E, Oleneva V, Yanke V, Papaioannou A, Mavromichalaki H, Gopalswamy N, Yashiro S 2014 Sol. Phys. 289 3949
[10]	Yu X X, Lu H, Le G M, Shi F 2010 Sol. Phys. 263 223
[11]	Zhao L L, Zhang H 2016 Astrophys. J. 827 13
[12]	Le G M, Han Y B 2005 Acta Phys. Sin. 54 467 (in Chinese)[乐贵明, 韩延本2005 54 467]
[13]	Dorman L I 2005 Ann. Geophys. 23 2997
[14]	Huang Y L, Fu Y F, Chen J M, Huang G S, Liu X N 2015 Res. Astron. Astrophys. 16 82 (in Chinese) [黄寅亮, 傅元芬, 陈济民, 黄更生, 刘小宁2015 天体 16 82]
[15]	Lockwood J A 1971 Space Sci. Rev. 12 658
[16]	Jokipii J R, Kota J 1986 J. Geophys. Res. 91 2885
[17]	Li G, Webb G M, Roux J A L, Zank G P, Wiedenbeck M E 2007 Numerical Modeling of Space Plasma Flows 385 31
[18]	Li G, Webb G M, Roux J A L, Wiedenbeck M, Florinski V, Zank G P 2009 Proceedings of the 31st ICRC Td, Poland, July 715, 2009
[19]	Fisk L A 1971 J. Geophys. Res. 76 221
[20]	Hattingh M, Burger R A, Potgieter M S, Haasbroeket L J 1997 Adv. Space Res. 19 893
[21]	Gil A, Iskra K, Modzelewska R, Alania M V 2005 Adv.Space Res. 35 687
[22]	Strauss R D, Potgieter M S, Bsching I, Kopp A 2012 Astrophys. Space Sci. 339 223
[23]	Gleeson L 1968 Astrophys. J. 154 1011
[24]	Caballero-Lopez R A, Moraal H, Mcdonald F B 2004 J.Geophys. Res. 109 361
[25]	Usoskin I G, Bazilevskaya G A, Kovaltsov G A 2011 J.Geophys. Res. 116 1
[26]	Usoskin I G, Kovaltsov G A, Adriani O, Barbarino G C,Bazilevskaya G A, Bellotti R 2015 Adv. Space Res. 55 2940
[27]	Dunzlaff P, Strauss R D, Potgieter M S 2015 Comput.Phys. Commun. 192 156
[28]	Zhang M 1999 Astrophys. J. 513 409
[29]	Wawrzynczak A, Modzelewska R, Gil A 2015 J. Phys.Conf. Ser. 574 012078
[30]	Pei C, Bieber J W, Breech B, Burger R A, Clem J,Matthaeus W H 2010 J. Geophys. Res. 115 333
[31]	Luo X, Zhang M, Rassoul H K, Pogorelov N V 2011 Astrophys. J. 730 13
[32]	Luo X, Zhang M, Feng X, Mendoza-Torres J E 2013 J.Geophys. Res. 118 7517
[33]	Bobik P, Boella G, Boschini M J, Consolandi C, Della Torre S, Gervasi M, Grandi D, Kudela K, Pensotti S,Rancoita P G, Tacconi M 2012 Astrophys. J. 745 132
[34]	Parker E N 1965 Planet. Space Sci. 13 9
[35]	Schuss Z 1980 Theory and Applications of Stochastic Differential Equations (New York: John Wiley) pp1095
[36]	Wiedenbeck M E, Davis A J, Leske R A, Binna W R,Cohen C M S, Cummings A C, de Nolfo G, Israel M H,Labrador A W, Mewaldt R A, Scott L M, Stone E C,von Rosenvinge T T 2005 Proceedings of the 29th International Cosmic Ray Conference Pune, India, August 310, 2005 p227
[37]	Chih P P, Lee M A 1986 J. Geophys. Res. 91 2903
[38]	Nishida A 1982 J. Geophys. Res. 87 6003
[39]	Richardson I G, Cane H V 2010 Sol. Phys. 264 189
[40]	Thomas B T, Gall R 1984 J. Geophys. Res. 89 2991
[41]	Wang C, Du D, Richardson J D, Liu Y 2005 Proceedings of the Solar Wind 11/SOHO16, Connecting Sun and Heliosphere Conference (ESA SP-592) Whistler,Canada, June 1217, 2005 p781
[42]	Verma P L, Patel N K, Prajapati M 2014 J. Phys. Conf.Ser. 511 012057
[43]	Ahluwalia H S, Alania M V, Wawrzynczak A, Ygbuhay R C, Fikani M M 2014 Sol. Phys. 289 1763
[44]	Mishev A L, Usoskin I G, Kovalstov G A 2013 J. Geophys.Res. 118 2783
[45]	Mangeard P S, Ruffolo D, Saiz A, Madlee S, Nutaro T 2016 J. Geophys. Res. 121 7435
[46]	Papailiou M, Mavromichalaki H, Belov A, Eroshenko E,Yanke V 2012 Sol. Phys. 276 337
[47]	Feng L, Inhester B, Wei Y, Gan W Q, Zhang T L, Wang M Y 2012 Astrophys. J. 751 18

施引文献

[1]	Rossi B 1964 Cosmic Rays (New York: McGraw-Hill) pp110
[2]	Blasi P 2013 Astron. Astrophys. Rev. 21 70
[3]	Bothmer V, Daglis I A 2007 Space Weather- Physics and Effects (Berlin: Springer) pp103130
[4]	Le G M 2002 Ph. D. Dissertation (Beijing: Chinese Academy of Sciences) (in Chinese) [乐贵明2002 博士学位论文(北京: 中国科学院空间科学与应用研究中心)]
[5]	Potgieter M S 2013 Living Rev. Sol. Phys. 10 3
[6]	Guo W J, Zhu B Y 1990 Chin. J. Spac. Sci. 10 247 (in Chinese) [郭维吉, 朱邦耀1990 空间科学学报10 247]
[7]	Kharayat H, Prasad L, Mathpal R, Garia S, Bhatt B 2016 Sol. Phys. 291 603
[8]	Cane H V 2000 Space Sci. Rev. 93 55
[9]	Belov A, Abunin A, Abunina M, Eroshenko E, Oleneva V, Yanke V, Papaioannou A, Mavromichalaki H, Gopalswamy N, Yashiro S 2014 Sol. Phys. 289 3949
[10]	Yu X X, Lu H, Le G M, Shi F 2010 Sol. Phys. 263 223
[11]	Zhao L L, Zhang H 2016 Astrophys. J. 827 13
[12]	Le G M, Han Y B 2005 Acta Phys. Sin. 54 467 (in Chinese)[乐贵明, 韩延本2005 54 467]
[13]	Dorman L I 2005 Ann. Geophys. 23 2997
[14]	Huang Y L, Fu Y F, Chen J M, Huang G S, Liu X N 2015 Res. Astron. Astrophys. 16 82 (in Chinese) [黄寅亮, 傅元芬, 陈济民, 黄更生, 刘小宁2015 天体 16 82]
[15]	Lockwood J A 1971 Space Sci. Rev. 12 658
[16]	Jokipii J R, Kota J 1986 J. Geophys. Res. 91 2885
[17]	Li G, Webb G M, Roux J A L, Zank G P, Wiedenbeck M E 2007 Numerical Modeling of Space Plasma Flows 385 31
[18]	Li G, Webb G M, Roux J A L, Wiedenbeck M, Florinski V, Zank G P 2009 Proceedings of the 31st ICRC Td, Poland, July 715, 2009
[19]	Fisk L A 1971 J. Geophys. Res. 76 221
[20]	Hattingh M, Burger R A, Potgieter M S, Haasbroeket L J 1997 Adv. Space Res. 19 893
[21]	Gil A, Iskra K, Modzelewska R, Alania M V 2005 Adv.Space Res. 35 687
[22]	Strauss R D, Potgieter M S, Bsching I, Kopp A 2012 Astrophys. Space Sci. 339 223
[23]	Gleeson L 1968 Astrophys. J. 154 1011
[24]	Caballero-Lopez R A, Moraal H, Mcdonald F B 2004 J.Geophys. Res. 109 361
[25]	Usoskin I G, Bazilevskaya G A, Kovaltsov G A 2011 J.Geophys. Res. 116 1
[26]	Usoskin I G, Kovaltsov G A, Adriani O, Barbarino G C,Bazilevskaya G A, Bellotti R 2015 Adv. Space Res. 55 2940
[27]	Dunzlaff P, Strauss R D, Potgieter M S 2015 Comput.Phys. Commun. 192 156
[28]	Zhang M 1999 Astrophys. J. 513 409
[29]	Wawrzynczak A, Modzelewska R, Gil A 2015 J. Phys.Conf. Ser. 574 012078
[30]	Pei C, Bieber J W, Breech B, Burger R A, Clem J,Matthaeus W H 2010 J. Geophys. Res. 115 333
[31]	Luo X, Zhang M, Rassoul H K, Pogorelov N V 2011 Astrophys. J. 730 13
[32]	Luo X, Zhang M, Feng X, Mendoza-Torres J E 2013 J.Geophys. Res. 118 7517
[33]	Bobik P, Boella G, Boschini M J, Consolandi C, Della Torre S, Gervasi M, Grandi D, Kudela K, Pensotti S,Rancoita P G, Tacconi M 2012 Astrophys. J. 745 132
[34]	Parker E N 1965 Planet. Space Sci. 13 9
[35]	Schuss Z 1980 Theory and Applications of Stochastic Differential Equations (New York: John Wiley) pp1095
[36]	Wiedenbeck M E, Davis A J, Leske R A, Binna W R,Cohen C M S, Cummings A C, de Nolfo G, Israel M H,Labrador A W, Mewaldt R A, Scott L M, Stone E C,von Rosenvinge T T 2005 Proceedings of the 29th International Cosmic Ray Conference Pune, India, August 310, 2005 p227
[37]	Chih P P, Lee M A 1986 J. Geophys. Res. 91 2903
[38]	Nishida A 1982 J. Geophys. Res. 87 6003
[39]	Richardson I G, Cane H V 2010 Sol. Phys. 264 189
[40]	Thomas B T, Gall R 1984 J. Geophys. Res. 89 2991
[41]	Wang C, Du D, Richardson J D, Liu Y 2005 Proceedings of the Solar Wind 11/SOHO16, Connecting Sun and Heliosphere Conference (ESA SP-592) Whistler,Canada, June 1217, 2005 p781
[42]	Verma P L, Patel N K, Prajapati M 2014 J. Phys. Conf.Ser. 511 012057
[43]	Ahluwalia H S, Alania M V, Wawrzynczak A, Ygbuhay R C, Fikani M M 2014 Sol. Phys. 289 1763
[44]	Mishev A L, Usoskin I G, Kovalstov G A 2013 J. Geophys.Res. 118 2783
[45]	Mangeard P S, Ruffolo D, Saiz A, Madlee S, Nutaro T 2016 J. Geophys. Res. 121 7435
[46]	Papailiou M, Mavromichalaki H, Belov A, Eroshenko E,Yanke V 2012 Sol. Phys. 276 337
[47]	Feng L, Inhester B, Wei Y, Gan W Q, Zhang T L, Wang M Y 2012 Astrophys. J. 751 18

[1]	严豪, 丁留贯, 封莉, 顾斌. 太阳高能粒子强度与日冕物质抛射及其II型射电暴的关系. , 2024, 73(7): 079601. doi: 10.7498/aps.73.20231855
[2]	金魁, 吴颉. 高温超导体组合薄膜和相图表征高通量方法. , 2021, 70(1): 017403. doi: 10.7498/aps.70.20202102
[3]	张荣培, 王迪, 蔚喜军, 温学兵. 基于广义交替数值通量的局部间断Galerkin方法求解二维波动方程. , 2020, 69(2): 020202. doi: 10.7498/aps.69.20190613
[4]	周坤论, 丁留贯, 王智伟, 封莉. 基于射电观测的日冕物质抛射驱动激波的统计特征研究. , 2019, 68(13): 139601. doi: 10.7498/aps.68.20190223
[5]	韩伟, 周丽丹, 李富全, 王芳, 冯斌, 郑奎兴, 巩马理. 光束通量空间分布随机变化的统计分析. , 2014, 63(7): 074204. doi: 10.7498/aps.63.074204
[6]	王小林, 周朴, 马阎星, 马浩统, 许晓军, 刘泽金, 赵伊君. 基于随机并行梯度下降算法的多波长激光相干合成. , 2010, 59(8): 5474-5478. doi: 10.7498/aps.59.5474
[7]	王小林, 周朴, 马阎星, 马浩统, 许晓军, 刘泽金, 赵伊君. 基于随机并行梯度下降算法光纤激光相干合成的高精度相位控制系统. , 2010, 59(2): 973-979. doi: 10.7498/aps.59.973
[8]	王晶, 冯学尚. 日冕物质抛射引起地磁扰动的分类预报. , 2007, 56(4): 2466-2474. doi: 10.7498/aps.56.2466
[9]	吴惠彬, 张永发, 梅凤翔. 求解微分方程的Hojman方法. , 2006, 55(10): 4987-4990. doi: 10.7498/aps.55.4987
[10]	王薇, 张杰, 赵刚. 吸积盘的X射线辐射对周围星际物质的离化研究. , 2006, 55(1): 287-293. doi: 10.7498/aps.55.287
[11]	陈建文, 高鸿奕, 朱化凤, 谢红兰, 李儒新, 徐至展. 中子相衬层析成像方法. , 2005, 54(3): 1132-1135. doi: 10.7498/aps.54.1132
[12]	杨进蔚, 宋先瑛, 李旭, 张炜, 郦文忠, 王世庆. HL-1M托卡马克中的中子通量和辐射剂量当量测量. , 2005, 54(4): 1648-1652. doi: 10.7498/aps.54.1648
[13]	周云松, 阮景辉, 陈金昌, 易向东. 中子单色化新方法. , 1997, 46(5): 929-937. doi: 10.7498/aps.46.929
[14]	戴子高, 陆埮, 彭秋和. 中子星内部非奇异-奇异夸克物质的相变. , 1993, 42(8): 1210-1215. doi: 10.7498/aps.42.1210
[15]	金继荣, 金新, 韩永胜, 吉和林, 谢晖, 张治平, 姚希贤. 低通量慢中子辐照对高Tc超导体正常态电阻的影响及其机理探讨. , 1993, 42(8): 1346-1351. doi: 10.7498/aps.42.1346
[16]	金继荣, 陈武鸣, 金新, 张贻瞳, 鹿牧, 吉和林, 刘新湖, 姚希贤, 易敏瑜, 洪苍生, 万韬隃, 方俊人. 低通量慢中子对Bi系超导体的辐照效应. , 1992, 41(5): 846-850. doi: 10.7498/aps.41.846
[17]	文根旺. 量子激发态最陡下降微扰理论. , 1991, 40(9): 1388-1395. doi: 10.7498/aps.40.1388
[18]	文根旺. 简并基态最陡下降微扰理论. , 1988, 37(12): 1981-1986. doi: 10.7498/aps.37.1981
[19]	黄昆, 刘东源. 最陡下降法和多声子跃迁的多频模型. , 1985, 34(6): 709-714. doi: 10.7498/aps.34.709
[20]	李竹起, 阮景辉, 吴善令, 杨同华, 何敏, 陆挺, 成之绪, 陈桂英, 叶春堂. 用于凝聚态物质研究的旋转晶体中子飞行时间谱仪. , 1980, 29(11): 1462-1470. doi: 10.7498/aps.29.1462

计量

文章访问数: 5906
PDF下载量: 274
被引次数: 0

姓名
邮箱
手机号码
标题
留言内容
验证码

搜索

留言板