1 / 47

国家重点基础研究发展计划 973 计划 课题 四 2011 年度总结报告

国家重点基础研究发展计划 973 计划 课题 四 2011 年度总结报告. 项目名称：日地空间天气预报的物理基础与模式 研究课题名称：磁层空间暴的能量传输与释放过程研究 课题成员：沈超，傅绥燕，曹晋滨，陈耿雄， 肖池阶，王慧，邓晓华，李毅 , 张永存等. 报 告 提 纲. 一、 2011 年研究计划 和目标 二、 年度计划任务完成情况 三、 2011 年度 研究工作的主要进展 四、总结. 一、 2011 年研究计划和目标. 计划研究内容 :

morgana-arlais
Download Presentation

国家重点基础研究发展计划 973 计划 课题 四 2011 年度总结报告

An Image/Link below is provided (as is) to download presentation Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author. Content is provided to you AS IS for your information and personal use only. Download presentation by click this link. While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server. During download, if you can't get a presentation, the file might be deleted by the publisher.

E N D

Presentation Transcript

  1. 国家重点基础研究发展计划973计划课题四2011年度总结报告国家重点基础研究发展计划973计划课题四2011年度总结报告 项目名称：日地空间天气预报的物理基础与模式 研究课题名称：磁层空间暴的能量传输与释放过程研究 课题成员：沈超，傅绥燕，曹晋滨，陈耿雄， 肖池阶，王慧，邓晓华，李毅, 张永存等

  2. 报 告 提 纲 一、 2011年研究计划 和目标 二、年度计划任务完成情况 三、 2011年度研究工作的主要进展 四、总结

  3. 一、 2011年研究计划和目标 计划研究内容: (1) 收集国内外相关的卫星和地面观测台站对磁层空间暴（亚暴-磁暴-粒子暴）的探测数据，并选择一些典型事件进行个例分析；自主开发具有原创性的多点卫星探测数据分析新方法。 (2) 利用双星，Cluster，Themis以及其它卫星的观测数据，进行统计和个例研究，综合分析各种行星际条件下，尤其强CME-磁层耦合期间，磁层顶三维磁重联X-线方向的分布规律，以及FTEs事件中磁通量管的结构特性；运用自主开发的多点卫星数据分析方法研究磁层顶边界层的全球拓扑结构及其内部结构； (3) 深入分析辐射带粒子动力学规律，研究磁层高能粒子的来源、加速和耗散机制；探讨磁层相对论电子的物理加速和扩散机制。

  4. 预期目标: • 建立磁层空间暴典型事件数据集；提出原创性的多点卫星探测数据分析新方法；初步分析磁层顶边界层的全球拓扑结构及其内部结构；初步认识磁层粒子暴物理机理。 • 完成学术论文8-12 篇。

  5. 二、年度计划任务完成情况 2011年度主要完成成果： 1. 提出原创性的多点卫星探测数据分析新方法，能够应用于任意颗卫星、任意卫星位形的多点探测数据； 2 . 第一次由探测分析得出地球磁层顶以及极尖区的几何结构和曲率半径的空间分布； 3 . 首次发现太阳风动压和IMF不同步现象，连续磁重联触发先于瞬时重联；Cluster 在高纬区探测到了反平行磁重联过程。

  6. 4 . 统计分析得出磁尾电流片大空间尺度结构分布（ AGU流行榜）； 5 . 曹晋滨等建立了第一个磁尾等离子体片内边界的全球模型； 6 . 磁尾等离子体片高速流与偶极化锋面特性统计分析; 7 . 磁层亚暴的多次分立爆发特征分析;

  7. 8. 亚极光区极化流产生机制及其与亚暴的关系; 9. 磁暴环电流离子分布的双星中性原子成像探测研究； 10．辐射带高能电子演化的变分辛积分算法研究 11. 初步建立辐射带高能粒子动态演化数值模式，具有良好的稳定性和准确度。

  8. 2011年已发表SCI论文共计18篇；在印和接受论文6篇。2011年已发表SCI论文共计18篇；在印和接受论文6篇。 获得空间探测仪器发明专利4项。 地球空间双星探测计划获2010年度国家科技进步奖一等奖（主要完成人：刘振兴等）； 2011年傅绥燕荣获国际宇航科学院通讯院士。 课题4已经按计划完成了2011年度课题研究任务，实现了预期目标。

  9. 三. 2011年度研究工作的主要进展

  10. 1．基于任意卫星星簇形态的多点探测来确定物理量空间梯度的方法1．基于任意卫星星簇形态的多点探测来确定物理量空间梯度的方法 物理量空间梯度 我们成功运用新方法分析THEMIS三点卫星探测，Cluster三点等离子体探测，计算了三点星簇平面法向的电流密度和等离子体流的涡度。

  11. K-H waves at flank magnetopause 左图 2001年11月20日磁层顶昏侧剪切流场激发的K-H涡旋的电流密度和流场涡度的分布

  12. 2. 地球磁层顶以及极尖区的几何结构和曲率半径的空间分布 第一次通过磁场曲率的直接测量来得到中等太阳风动压条件下晨昏子午面向阳侧磁层顶表面的曲率半径的空间分布。 直接证实了极尖区磁力线凹点的存在。 在磁层顶附近磁力线的最小曲率半径都大约是2Re，无论是极尖区太阳一侧还是磁尾一侧。 越靠近赤道面，磁层顶的曲率半径越大；当接近日下点时，磁层顶的曲率半径和其地心距离相当。

  13. 磁层顶曲面的几何结构 [Shen et al.,2011] 磁层顶曲面的拓扑位形 磁层顶曲面曲率-天顶角关系

  14. 图4 磁层顶边界层曲率半径随极向角的空间分布特性

  15. 图5 磁层顶边界层拓扑结构示意图

  16. 3. 磁层顶磁重联过程及其对太阳风变化的响应 综合Cluster和“双星计划”TC-1等卫星的协同观测数据，发现在日下点存在反平行磁重联过程； 首次发现太阳风动压和IMF不同步现象，动压先于IMF到达磁层顶；连续磁重联触发先于瞬时重联； 这些结果为建立向阳面磁层顶三维磁重联的磁场结构和物理模型提供更加细致的观测基础。

  17. 4、低纬Pi2脉动产生机制研究 Pi2脉动作为磁层亚暴的主要特征之一，一直是国内外学者所关心的研究课题。我们利用THEMIS同时分别位于等离子体层内、外三颗卫星的观测资料，研究了两个伴随中等强度亚暴而发生的Pi2事件的产生和传播机制。 结果表明，卫星和地面观测到的Pi2脉动之间的相关性高度依赖卫星相对等离子体顶的位置，而等离子体层的空腔共振是低纬Pi2脉动的主要形成机制。

  18. 图2. P3（蓝线），P4（红线），P5（绿线）和ASC台站在磁地方时-磁壳参数L值坐标系中的分布。实心圆点代表相应卫星穿过等离子体层顶时的位置。而粗短线表示对应于两个Pi2事件的时间区间（0015 ～ 0040 UT）。 图1. AU/AL指数以及ASC和HER两个地面台站分别记录到的两个Pi2事件。

  19. 4 . 磁尾电流片大空间尺度结构分布统计分析 基于Cluster卫星簇在2001、2003、2004年期间在距地心15-19个地球半径（RE）范围内对磁尾磁场的多点同步观测，本文在GSM坐标系下统计性地研究了地球磁尾电流片中心处的磁场结构。

  20. GSM，磁场及其分量的统计分布 磁场几何结构的空间分布 磁场几何结构的频次统计分布 －Bz的出现概率，及其曲率方向的分布

  21. 电流片法向 倾斜电流片出现概率的分布 80％的倾斜电流片中， 磁力线倾角 是小于450 这表明倾斜电流片中的磁力线 结构大多是相互滑移的。 电流片的分类

  22. 最小曲率半径Rc,min与中性片半厚度h • 在方位角~140°-~200°, 即磁地方时~21:00-01:00, Rc,min(0.6-1RE) 与h (0.2-0.4RE) 相对较小 • 而在晨昏两翼Rc,min(3-5RE) 和h (0.7-1.6RE) 相对较大。 • 薄电流片（h<1000km)容易在方位角~140°-~200°范围内出现。

  23. [Shen et al.,AG.2008c;Rong et al,PSS2010 ] 电流片中心处，电流密度的方向 电流密度的平均强度为 4~8nA/m2,基本沿电流片 切向方向，由晨侧“西向” 流向昏侧方向。 尤其是Y>10RE.的区域

  24. 统计结果表明，磁尾电流片存在明显的晨昏不对称性。在磁地方时～21:00–～01:00范围内，电流片普遍较薄，负Bz出现概率较高。这表明，在该区域范围内磁尾能量容易通过重联或电流中断等活动得到释放，相应若将该电流片区域映射到极区电离层，那么对应亚暴极光活动在磁地方时～21:00–～01:00范围内的发生概率也会相对较高。这与IMAGE-FUV的统计观测是相一致的。统计结果表明，磁尾电流片存在明显的晨昏不对称性。在磁地方时～21:00–～01:00范围内，电流片普遍较薄，负Bz出现概率较高。这表明，在该区域范围内磁尾能量容易通过重联或电流中断等活动得到释放，相应若将该电流片区域映射到极区电离层，那么对应亚暴极光活动在磁地方时～21:00–～01:00范围内的发生概率也会相对较高。这与IMAGE-FUV的统计观测是相一致的。 极光亚暴随磁地方时的分布（IMAGE-FUV ）(引自Frey et al., 2004).

  25. 5. 磁尾等离子体片内边界的全球模型 等离子体片的地向渗透一直是空间天气研究的重要内容。2011年曹晋滨等人利用TC-1卫星2004－2006年的离子探测数据，研究了等离子体片内边界环绕地球的分布，以及随地磁扰动Kp指数的变化规律，并建立了第一个磁尾等离子体片内边界的全球模型。

  26. 平静时期和活动时期等离子体内边界赤道面分布平静时期和活动时期等离子体内边界赤道面分布

  27. 6 . 磁尾等离子体片高速流与偶极化锋面特性统计分析 利用搭载于Cluster四颗卫星上的磁场和等离子体观测仪器，对2001和2002两年发生在磁尾等离子体片中的高速流事件中的磁场变化进行了统计研究。统计结果表明，在高速流前端，伴随着等离子体整体速度的增加，磁场Bz分量出现先短暂减小然后剧烈增加的现象。

  28. Cluster卫星观测到的219个高速流事件统计图

  29. 7.磁层亚暴的多次分立爆发特征分析 利用THEMIS五颗卫星的数据，同时结合THEMIS地面台站的极光和地磁观测，研究了2008年3月1日UT 08:30-10:00时间段的一个中等强度亚暴事件。发现不同的卫星在不同的时间观测到了不同数量的偶极化事件，这一系列的偶极化的触发位置有一个明显的尾向移动，而每一次偶极化事件又有一个地向的传输过程。这些观测结果和地面极光的增亮过程可以很好的联系起来。同时，在一次次偶极化事件中，我们会看到离子和电子的能通量都有明显的增加，离子和电子有地向注入的过程。综上所述，我们观测到了亚暴过程中局地瞬时的磁能释放过程。

  30. 亚暴过程中磁尾一系列偶极化事件演化过程示意图亚暴过程中磁尾一系列偶极化事件演化过程示意图 [Ma, Shen, et al., in submission, 2011]

  31. 8. 磁暴环电流离子分布的双星中性原子成像探测研究 分析双星计划TC-2NUADU仪器第一手的中性原子成像数据，开展环电流离子分布形态以及动态演化的研究。探测数据分析结果显示，磁暴期间南向的行星际磁场（IMF）分量在离子从磁尾向内磁层注入和随后的环电流增长过程中起着关键的作用；能量越高的离子越不容易进入小L值区域被捕获。利用环电流离子稳态模型模拟计算环电流区的ENA通量，与NUADU的观测数据进行对比分析；模拟计算ENA通量空间分布与NUADU观测ENA图像一致。

  32. 磁暴环电流区中性原子成像探测与模拟结果比较磁暴环电流区中性原子成像探测与模拟结果比较 图5 2004-11-07 磁暴（Kp=8，Dst = -88nT ）期间，于23:08:49UT，r = 4.18 RE，NUADU观测ENA图像（左列）和模拟ENA图像（右列）的对比

  33. 图5 2004-11-07 磁暴（Kp=8，Dst = -88nT ）期间，于23:49:19UT，r = 2.95 RE，NUADU观测ENA图像（左列）和模拟ENA图像（右列）的对比

  34. 环电流离子的地方时分布具有能量无差异特性，不同能量环电流离子的注入区域、漂移路径和损失过程基本相同。 环电流离子的地方时分布具有能量无差异特性，不同能量环电流离子的注入区域、漂移路径和损失过程基本相同。 • 主相期间观测的ENA通量在方位角100度附近达到最大，模拟图像通量最大方位角位置与观测图像不同。 观测与模拟比较研究 – 地方时分布特性 2004-11-07 23:08:49UT，主相，Kp=8

  35. +z 4524 km S1 11.25º T13 T14 T15 26343 km D1 E2 T8 D2 S2 A B A1 D 8287 km T11 C2 T12 +y 60º B1 O C B2 6400 km A2 C1 Earth 600 km Exobase • 模拟图像具有更宽的仰角分布范围，通量最大值方向相对远离地球边界线。 • 模拟计算采用的地冕中性层模型在低高度区域低估了中性层密度。 • ENA电离损失：（2004-11-08 03:53:37UT） 观测与模拟比较研究 – 高度/纬度分布特性

  36. 9．辐射带高能电子演化的变分辛积分算法研究9．辐射带高能电子演化的变分辛积分算法研究 变分辛积分方法保持系统的辛结构，其数值计算结果在整个积分时间里都保证误差在一定的小范围内。在长时间计算中，辛积分方法比传统的RK4方法更有优越性。计算结果清晰显示RK4方法所得粒子轨迹发生扭曲，而变分辛积分方法在总时间内保持误差在一定范围内。事实上，这种辛积分这种特性在任意长时间的计算里都能保持。在将来波粒子相互作用的研究中，我们将采用更精确的磁场与电场模型。另外磁层中屏蔽电场也将考虑。

  37. 图（3） 考虑对流电场与共转电场，4KeV质子在磁层中的“香蕉”轨道。（a）为RK4计算结果，（b）为辛积分计算结果。

  38. 10. 磁暴环电流离子动态演化数值模式 依据环电流粒子暴的动态演化的福克-普朗克方程，初步编制了环电流粒子动态演化计算程序，具有良好的稳定性和准确度。此数值模式能够描述各能量段（中能和高能）环电流粒子在各种磁层电场作用下的全球空间分布，能谱特征以及投掷角分布。数值模拟分析了若干种典型情况。

  39. 环电流离子全球演化的动力学模拟 磁暴注入离子(能量100Kev, 投掷角90度)在内磁层的磁飘移运动

  40. 四. 总结 2011年度主要完成成果： 1. 提出原创性的多点卫星探测数据分析新方法； 2 . 由探测分析得出地球磁层顶以及极尖区的几何结构和曲率半径的空间分布； 3 . 首次发现太阳风动压和IMF不同步现象，连续磁重联触发先于瞬时重联；Cluster 在高纬区探测到了反平行磁重联过程。

  41. 4 .磁尾电流片/等离子体片结构:统计分析得出磁尾电流片大空间尺度结构分布,建立了第一个磁尾等离子体片内边界的全球模型； 5 .磁层亚暴多时空演化规律:统计分析磁尾等离子体片高速流与偶极化锋面特性,分析了磁层亚暴的多次分立爆发特征; 6. 磁暴环电流离子分布的双星中性原子成像探测研究； 7．开发了辐射带高能电子演化的变分辛积分算法，初步建立环电流离子全球演化的动力学数值模式。

  42. 2011年已发表SCI论文共计18篇；获得空间探测仪器发明专利4项。2011年已发表SCI论文共计18篇；获得空间探测仪器发明专利4项。 地球空间双星探测计划获2010年度国家科技进步奖一等奖（主要完成人：刘振兴等）；2011年傅绥燕荣获国际宇航科学院通讯院士。 课题4已经按计划完成了2011年度课题研究任务，实现了预期目标。

  43. 五、2012年度研究计划 2012年度研究计划如下： (1) 利用双星, Cluster， THEMIS和ACE等多点卫星全球联合探测数据，对磁层亚暴进行大量的个例分析以及统计分析研究；分析近地磁尾等离子体片结构，建立经验模式；分析近地磁尾电流片的不稳定性及效应，粒子加速与加热机制；对等离子体片中高速流与磁结构的关联以及演化问题开展统计研究；着重开展卫星资料的联合观测，注重新的现象和物理过程的发现，为建立新的理论模型提供依据；分析研究亚暴在近地磁尾蓄积磁能，触发以及能量释放的整体演化过程； (2)分析研究磁暴期间，尤其超强磁暴期间，磁层-电离层系统中的电场和电流体系全球三位分布及演化，概括磁层磁暴多尺度特性及演化规律，揭示磁暴-亚暴物理关系； (3) 研究相对论电子事件与地球磁暴和亚暴的关系；在充分分析辐射带粒子暴机理的基础上，建立辐射带粒子暴动力学方程，编制数值模拟程序。

  44. 预期目标： 建立磁尾动力过程的若干模式，为揭示亚暴机制建立重要基础； 获得磁层-电离层系统中的电场和电流体系全球三维分布及演化规律； 认识磁暴-亚暴关系物理图像； 改进辐射带粒子暴模型，数值模拟粒子暴的全球发展过程，服务于航天活动空间天气安全保障。 完成学术论文8-12 篇。

  45. 谢 谢！

More Related