太阳风向磁层输入电磁能量研究

利用全球磁流体力学（MHD）模拟结果,通过确立包含磁层顶的太阳风流线内边界来识别三维磁层顶位形,并以极尖区位置作为磁层顶日侧与夜侧的分界线,在此基础上定量研究了不同条件下穿过磁层顶向磁层内输入的电磁能量. 研究发现,磁层顶的能量传输与太阳风条件密切相关,磁重联是控制电磁能量传输的重要机制. 结果表明,当IMF（行星际磁场）南向时,极尖区后方的磁尾附近存在电磁能输入最大值,当IMF北向时,电磁能输入最大值发生在极尖区附近;南向IMF条件下,在IMF强度增大或太阳风密度增大时,磁层顶电磁能传输的电磁能量比北向IMF条件时增加更显著. 太阳风通过调节磁层顶面积间接影响到磁层顶能量传输大小. 研究还发现,北向IMF与南向IMF条件下穿过磁层顶的电磁能输入的比值范围约为10%～30%,此比值一定程度上反映了北、南方向IMF与地磁场磁重联效率的比值.      

极端太阳风条件下的磁层顶位形

基于极端太阳风条件下的三维MHD数值模拟数据, 构建了一种极端太阳风条件下的三维非对称磁层顶位形模型. 所提出的模型考虑了行星际南向磁场(IMF) B_z日下点距离侵蚀的饱和效应, 太阳风动压B_d对磁层顶张角影响的饱和效应, 赤道面、昼夜子午面磁层顶的不对称性以及极尖区的内凹结构和内凹中心的移动, 并利用Levenberg-Marquart多参量非线性拟合方法拟合了模型参数. 数值模拟研究表明, 在极端太阳风条件下, 随B_d增大, 磁层顶日下点距离减小, 磁层顶磁尾张角几乎不变; 随南向(IMF)B_z增大, 磁层顶日下点距离略有减小, 磁层顶磁尾张角减小, 极尖区内凹中心向低纬移动. 通过对2010年8月1日太阳风暴事件验证发现, 本文所建立的模型能够描述极端太阳风条件下的三维磁层顶位形.      

基于MHD模拟数据的正午午夜子午面磁层顶位形研究

通过分析太阳风-磁层-电离层系统的三维全球磁流体力学(MHD)模型的计算数据, 给出了正午-午夜子午面磁层顶位形的定量模型. 分析表明, 正午-午夜子午面磁层顶位形可以用文献[3]提出的基于卫星观测数据的、描述赤道面磁层顶位形的函数来描述. 与赤道面磁层顶不同, 正午-午夜子午面磁层顶位形更为复杂. 在忽略极尖区(cusp)的简化条件下, 磁层顶位形仍需利用两条曲线来拟合. 太阳风动压Dp与行星际磁场分量Bz是控制磁层顶位形的主要因素. 行星际磁场为北向时, 磁场增强, 日下点距离r0增大; 行星际磁场为南向时, 磁场增强, 磁层顶日下点距离r0减小. 整体而言, 行星际磁场分量Bz由南转北时, r0增大, 且Bz对r0的影响减弱. 太阳风动压Dp是控制磁层顶日下点的主要因素, Dp增大, r0减小. 磁层顶位形的另一个参数磁层顶磁尾张角α, 随着行星际磁场南向分量增强而增大, 即磁层顶张开程度更加显著, 更多的磁通量由向阳侧传输到夜侧; Dp增大, α略增大, 这意味着Dp对磁通量由日侧向夜侧的传输也有一定的贡献.      

磁鞘快速流引起的磁层顶凹陷事件

最近研究表叽磁层顶凹陷对磁层-电离层耦合具有重要作用.但是,磁层顶凹陷现象的确认需要多颗卫星的联合观测,目前为止报道的磁层顶凹陷事例非常少.本文利用THEMIS 5颗卫星的联合观测结果,分析了一例由磁鞘快速流引起的磁层顶凹陷事件.2007年7月21日10:00UT-10:45UT,位于日下点磁层顶附近的THEMIS卫星...     

三维试验粒子轨道法在磁层粒子全球输运中的应用

根据磁层粒子动力学理论, 通过偶极磁场模型验证利用三维试验粒子轨道方法模拟近地球区(r < 8R_e)带电粒子运动特征的可靠性. 在此基础上, 以太阳风和磁层相互作用的全球MHD模拟结果为背景, 利用三维试验粒子轨道方法, 对非磁暴期间南向行星际磁场背景下太阳风离子注入磁层的情形进行数值模拟, 并对北向行星际磁场背景下太阳风离子注入极尖区以及内磁层的几种不同情形进行了单粒子模拟. 模拟结果反映了南向和北向行星际磁场离子向磁层的几种典型输入过程, 揭示出行星际磁场南向时太阳风粒子在磁层内密度分布的晨昏不对称性以及其在磁鞘和磁层内的大致分布, 并得出统计规律. 模拟结果与理论预测和观测结论相一致, 且通过数值模拟发现, 行星际磁场北向时靠近极尖区附近形成的非典型磁镜结构对于能量粒子经由极尖区注入环电流区域过程有重要的影响和作用.      

磁层顶通量传输事件轴向统计分析

磁层顶通量传输事件（Flux Transfer Event,FTE）与磁重联相关,其典型特征为磁场法向分量的双极变化.在不同FTE模型里,FTE结构可能为重联的通量管、由多X线重联形成的闭合磁通量绳或者由单X线重联形成的开放磁场环,从而在磁层顶有不同的整体位形.使用一种新的轴向分析方法,对Cluster在一个向阳面磁层顶穿越季观测到的505个FTE进行统计研究.结果表明:在磁层顶中低纬度的侧翼,大多数FTE轴向均为沿磁层磁力线方向即南北方向,少数FTE轴向沿着不同于磁层磁力线方向的东西方向;在高纬磁层顶,大多数FTE轴向沿东西方向,少数FTE轴向沿着磁层磁力线方向即南北方向.这些统计特征有助于重新认识FTE的全球形态.      

磁层顶穿越事件自动识别算法

磁层顶是太阳风与磁层进行质量、动量、能量交换的关键区域.磁层顶穿越事件(MCEs)可通过对卫星探测到的粒子能谱和磁场数据图进行人工分析的方式来识别.因太阳风动压和行星际磁场的易变性,位于磁层顶附近的卫星经过长期观测可能会经历成千上万次的磁层顶穿越.人工分析的方法工作量巨大,而且识别速度慢.本文发展了一种新的日下点附近M...     

磁层顶低混杂漂移不稳定性及等离子体幔、磁层对流的形成机制

磁层顶低混杂漂移不稳定性的理论和观测使我们可以提出一个新的磁层对流驱动模式,为了解释磁层对流的形成、磁层顶厚度等一系列磁层现象,已经提出了三种磁层模式,Dunge提出互联模式,认为行星际磁场磁力线与地磁场磁力线在磁层顶前部相互联接起来,磁层顶为一旋转间断面,太阳风粒子可直接通过磁层顶进入磁层内,虽然这一     

基于能量色散特征计算磁层顶重联位置的方法研究

磁层顶磁场重联是太阳风向磁层输入能量的主要方式.重联如何触发一直是空间物理研究的难点,其机制仍然有待深入研究.由于卫星穿越磁层顶时,很难恰好穿越重联发生的区域,因此难以观测到重联的触发条件.本文利用THEMIS卫星观测,确立了反演磁层顶重联点的方法.当重联刚开始发生时,卫星能够观测到离子的能量色散特征,可利用其计算卫星到重联发生位置的距离.沿着磁力线模型追踪该距离即可反演出磁层顶发生重联的位置.与其他方法进行了对比分析,结果显示本文方法比其他方法具有更高的精度.      

磁层顶附近流场剪切度的确定

磁层顶附近的流场剪切度与磁层顶附近能量转换的程度有关.很多磁层顶数值模拟用到流场剪切度这个输入参数,但一直是假设的.本文利用Cluster多卫星同时观测数据及独特的时空分辨功能,采用线性插值和重心坐标的方法确定了磁层顶附近流场的剪切度.通过对晨侧和昏侧磁层顶及附近磁层磁鞘流场剪切度的真实空间分布的研究结果表明,在平静的太阳风条件和地磁条件下磁层顶附近流场剪切度有时也很大,可达每百公里330 km/s的相对速度差.但在很多情况下流场是弱剪切的,在上千公里的距离上只有每秒几十公里的相对速度差.本文确定流场剪切度的方法可以推广用来确定任一位置的流场剪切度.      

The magnetopause erosion and the magnetosheath magnetic field penetration into the dayside magnetosphere

The earthward displacement of the magnetopause observed during a southward IMF (or the magnetopause erosion) and its dependence on the solar wind plasma and magnetic field parameters is studied by investigating data of about 30 magnetopause crossings by the ISEE 1 and 2 spacecraft. It is shown that the magnetopause erosion may be explained by a depression of the magnetic field intensity in the dayside magnetosphere caused by the penetration of the magnetosheath magnetic field (component perpendicular to the reconnection line) into the magnetosphere. The penetration coefficient (the ratio of the intensity of the penetrated field to the intensity of the magnetosheath magnetic field) is estimated and found to equal approximately 1.     

土星磁层顶的Kelvin-Helmholtz不稳定性和磁流体表面波

本文采用可压缩磁流体力学近似并选取Voyager-1的观测数据,讨论了土星磁层顶的Kelvin-Helmholtz不稳定性和磁流体表面波.计算表明,土星磁层随星体自转同磁鞘流之间的速度差,可以在地方时午前、中午和距中午不远的午后磁层顶激发不稳定性;午前和午后区有很强的非对称性.计算所得的波的特性与观测一致.      

磁层顶对磁层磁场的影响

太阳风与地磁场相互作用形成的磁层顶对磁层内磁场有重要影响。本文假定地磁场为偶极子,太阳风为理想导体,在太阳风与磁层的边界上满足磁场法向分量为零的边界条件。采用最小二乘法求得磁层顶电流在磁层内产生的磁场的球谐系数。从计算结果可以看出磁层顶对磁层磁场的影响。结果表明,向阳面的磁场、中性点、极光区的位置与形状与实际观测比较接近;磁尾磁场与实际观测相差较远,原因是没有加上磁尾片电流。文中还给出了太阳风与地磁轴交不同角度时的磁层磁场的计算结果。      

A Theoretic Interpretation of Movement of the Cusp Equatorward Boundary

A preliminary model is proposed to describe quantitatively the position and movement of cusp equatorward boundary. This integrated model, consisting of an empirical model of the magnetopause and a compressed dipolar model of Open/Closed field line, connects quantitatively the solar wind conditions, subsolar magnetopause and cusp equatorward boundary. It is shown that the increasing solar wind dynamic pressure and the increasing southward Interplanetary Magnetic Field (IMF) component drive the magnetopause to move inward and the cusp equatorward. This model is adopted to interpret quantitatively the cusp movement of August 14, 2001 observed by Cluster. The results show that the subsolar magnetopause moved earthward from 10.7 He to 9.0 Re during the period of 002300-002800 UT, and correspondingly the cusp equatorward boundary shifted equatorward. The observations of Cluster C1 and C4 show the cusp equatorward boundary that Cluster Cl and C4 were crossing during same interval moved equatorward by 4.6°. The cusp equatorward boundary velocity computed in the theoretical model (10.7km/s) is in good agreement with the observed value (9.4km/s) calculated from the data of CIS of Cluster C4 and C1.     

行星际磁场北向时磁层顶区磁场重联的全球模式

在对背阳面磁层顶区局域磁场重联模拟的基础上提出了一个行星际磁场北向时磁层顶磁场重联的全球模式。行星际磁场北向时碰层顶磁场重联导致近地尾瓣的能量被输送到远磁尾,太阳风能量不在磁尾储存,向阳面磁层顶变厚,磁层受到一系列扰动。      

Magnetopause poleward of the cusp: Comparison of plasma and magnetic signature of the boundary for southward and northward directed interplanetary magnetic field

A coherent data set of high-latitude dayside magnetopause encounters by old (Heos 2, Hawkeye, Prognoz 7, 8) and new (Polar, Interball Tail, Cluster) spacecraft is needed to build a realistic model of the magnetopause (MP) including an indentation in the cusp. In building such a coherent data set a caution is necessary as the dayside magnetopause at high-latitudes may be less clearly defined than in the case of observations at low latitudes. It is due to expected presence of bundles of newly-reconnected magnetic field lines forming an extended boundary layer on the magnetosheath (MS) side of the magnetopause in the cusp region. Moreover, numerical magnetohydrodynamic (MHD) models of the solar wind-magnetosphere interaction predict that under northward interplanetary magnetic field (IMF) an additional thin current sheet should form inside the magnetopause at high latitudes on the dayside (e.g., Wu, 1983; Palmroth et al., 2001). Such a thin currect sheet is absent in empirical magnetosphere models. This internal current sheet, if a real one, may be mistaken for the magnetopause if magnetic field data are only taken into account and/or plasma data are unavailable. The Interball-Tail orbit allows for a full transition of magnetopause boundary layers at high-latitudes. We compare plasma and magnetic field signatures of the magnetopause poleward of the cusp for southward and northward IMF. The distance between the magnetic signature of the magnetopause (the current layer) and a cold and laminarly antisunward flowing MS plasma (so called free-flow MS) was found to be 0.5 to 1 R_E, at least. These observations were made under nominal solar wind of v350 km/s and p_dyn=1 to 4 nPa. We also observed several transient magnetic field reversals in the cusp related to pulses of solar wind dynamic pressure and/or the IMF discontinuity arrival. These transient reversals occurred at the same distance to the model MP as well defined full MP crossing, so most probably they represent just short encounters with the magnetopause current layer. Our analysis suggests that an indentation of the magnetopause with a subtle structure dependent on the local magnetic shear would explain and allow to predict the magnetic configuration in the high-altitude cusp.     

On the modular approach of a 3D model of the system magnetosheath–magnetosphere

This article presents some preliminary features of a new self consistent model of the system magnetosheath–magnetosphere, recently developed in the Institute of Mechanics, Sofia, Bulgaria. The flow in the magnetosheath is governed by 3D ideal gas-dynamic equations. The positions and the shapes of the bow shock and the magnetopause are calculated iteratively as a part of the solution. These surfaces are essentially three-dimensional (generally non-axially-symmetric). The self-consistency between the regions is ensured via the boundary conditions. The magnetopause cusp indentations are formed, influencing essentially the magnetosheath flow. Prediction of the position and the shape of the bow shock for different values of the sonic Mach number are derived. Distribution of some flow parameters in the magnetosheath is presented. 3D numerical finite element model, calculating the field due to the magnetopause currents for an arbitrary magnetopause geometry, is used in the magnetosphere. The fields due to the current systems inside the magnetosphere(cross-tail current, ring current, and Birkeland current) are taken from the Tsyganenko empirical model. The magnetopause surface is calculated from the requirement the outside gas-dynamic pressure to be balanced by the magnetic pressure inside. The magnetosphere model can be viewed as an improved version of the empirical model but with more realistic magnetopause form and shielding field. Not a final but a beta version is used in this approach. The final model version as well the model details will be presented elsewhere.