2001年3月2日磁通量传输事件特性的研究

2001年3月2日11:00 至11:15 UT 期间,Cluster Ⅱ在南半球极尖区晨侧附近磁鞘内探测到3个通量传输事件（简称FTEs）. 本文利用Cluster Ⅱ星簇4颗卫星观测到的磁场和等离子体资料研究了这些通量传输事件的磁场形态和粒子特征. 并利用它们探测到的空间磁场梯度资料由安培定律直接求出星簇所在区域的电流分布. 结果指出：（1）BY占优势的行星际磁场结构在磁层顶的重联可以在极尖区附近发生;（2）FTEs通量管形成初期内外总压差和磁箍缩应力不一定平衡,达到平衡有一发展过程;（3）FTEs通量管截面在L M平面内的线度约为1.89RE;（4）FTEs通量管中等离子体主要沿轴向场方向流动,整个通量管以慢于背景等离子体的速度沿磁层顶向南向尾运动;（5）FTEs通量管中不仅有轴向电流,也存在环向电流. 轴向电流基本沿轴向磁场方向流动. 轴向和环向电流在管内均呈体分布,因而轴向电流产生的环向磁场接近管心时不断减小到零,而环向电流生成的轴向场则不断增大到极值;（6）在通量管的磁鞘部分观测到磁层能量粒子流量的增强,这表明通量管通过磁层顶将磁鞘和磁层内部连通起来了.     

2001年1月26日高纬磁层顶通量管事件的观测研究

2001年1月26日11:10～11:40UT, ClusterⅡ卫星簇位于午后高纬磁鞘边界层和磁鞘区,此 时行星际磁场Bz为南向. 本文对在此期间观测到的多次磁通量管事件作了详细的研究 ,获得一系列的新发现：（1）高纬磁鞘边界层磁通量管的出现具有准周期性,周期约为78s ,比目前已知的磁层顶向阳面FTE的平均周期（8～11min）小得多. （2）这些通量管都具有 强的核心磁场;其主轴多数在磁场最小变化方向,少数在中间变化方向,有些无法用PAA判 定其方向（需要用电流管PAA确定）,这与卫星穿越通量管的相对路径有关. （3）每个事件 都存在很好的HT参考系,在HT参考系中这些通量管是准定常态结构;所有通量管都沿磁层顶 表面运动,速度方向大体相同,都来自晨侧下方. 通量管的径向尺度为1～2RE, 与通 常的FTE通量管相当. （4）起源于磁层的强能离子大体上沿着管轴方向由磁层向磁鞘运动; 起源于太阳风的热等离子体沿管轴向磁层传输. 通量管为太阳风等离子体向磁层输运和磁层 粒子向行星际空间逃逸提供了通道. （5）每个通量管事件都伴随有晨昏电场的反转,该电 场为对流电场.     

2004年3月18日高纬磁层顶多重通量管事件分析

2004-03-18 23:10～23:50 UT期间，“双星（Double Star）”探测一号卫星（TC 1）在向阳面磁层顶高纬晨侧由内向外穿越磁层顶，其时TC_1的GSM坐标为 (75RE, -55RE, -54RE), RE为地球半径.穿越过程中TC_1观测到了8个通量管和1个磁通量传输事件（FTEs）.在此期间Cluster星簇位于向阳面太阳风内，其GSM坐标为(180RE, -31RE, -62RE)，其4颗卫星监测到行星际磁场（IMF）的BZ分量持续南向，BＹ有较大的负值.本文的研究表明：TC_1观测到的前7个通量管具有准周期重现性，周期大约是1～4 min，明显小于以前所观测到的FTEs的平均周期(8～11 min)；所有的通量管都具有较强的核心场.本文分别使用最小方差分析法（MVA）和Grad_Shafranov反演方法（GSR）对通量管的轴向进行了分析和对比，发现所有的通量管主轴基本沿晨昏向，结果显示GSR方法在轴向分析上比MVA优越.本文使用GSR方法对通量管的磁场结构进行了分析，恢复出了通量管的磁场在卫星穿越面的结构图；此外，本文还对这次多重通量管事件进行了deHoffmann Teller(HT)分析，结果表明，所有通量管大致朝南极方向运动，均来源于向日面低纬区域.这说明它们可能起源于向日面低纬区，由该区的磁场分量重联产生.     

近尾无碰撞多重X线重联的观测

利用Cluster星簇sc3卫星2001年10月1日09:46~09:50UT时段磁场和等离子数据,探测到近尾尾向传播的两个磁通量绳,时间间隔为26 s.通量绳结构的尾侧存在高速地向流,地侧有高速的尾向流.同时观测到了尾向流的重联X线源区和地向流的重联出流区,直接地观测到了近尾无碰撞多重X线重联.     

木卫二尾迹的数值模拟

木卫二在木星磁层中运动,与快速共转的木星磁场相互作用形成木星极光中木卫二的尾迹,此尾迹反映了连续扰动的磁通量管对木星电离层电流注入的过程.本文应用细丝理论对此过程进行数值模拟,未扰动的上游磁通量管经过大约720 s的时间与木卫二作用形成角度约为4°且相对于木卫二静止的拉伸磁通量管,随后在磁张力的作用下加速运动.在随后大约2160 s的加速时间内,磁通量管的速度逐渐与背景磁力线接近,向木星电离层注入的电流也由9×10⁵A逐渐衰减至0.以1/e初始电流处认为是可观测的木卫二尾迹限度,则木卫二的尾迹仅大约有一至两个经度,这与哈勃望远镜观测到的木卫二尾迹大小大致相同.木卫二尾迹与木卫一尾迹相比,输送的能量约为2.3×10¹³J,仅有木卫一的1/10,导致尾迹长度远小于木卫一的尾迹长度.此外,模拟过程中发现木卫二下游的压力槽对磁通量管运动的影响很小,可以忽略.     

孤立磁拱的约束平衡位形

本文讨论分层大气中孤立约束磁通量拱的二维静平衡位形．通量拱中的磁场初步假定为无力场，通过利用有限元计算方法，得到这个非线性自由边界问题的数值解．结果清楚地表明磁拱的内外场的耦合关系，并讨论了大气压力分布，尤其是大气标高对通量拱位形的影响．     

基于COMSOL Multiphysics的长白山天池多物理场耦合响应研究

采用COMSOL Multiphysics建立长白山天池2D多物理场模型,模拟长白山天池流场、自然电位分布及其磁场响应.研究结果表明:天池和大气降水入渗补给地下水,地下水在重力作用下向下流动时,产生电位差,形成自然电场,进而产生磁场.在火山通道外侧,电流密度方向由地表指向地下,通道周边地区磁通量密度呈现高值;在火山通道/断层破碎带内,向下运移的地下水遇到致密岩体后向上补给天池或涌出地表成泉,该过程产生向上运移的电流,电流密度方向由地下指向地表,磁通量密度呈现低值.研究长白山天池多物理场响应,有助于进一步理解长白山天池地区地下水循环过程及电磁异常出现的机理,亦可用于指导火山地区勘查、监测与开发.     

内磁层电流体系的地磁效应

对Ⅱ区场向电流及其伴随的部分环电流和电离层电流组成的内磁层三维电流体系（PRFI电流系）的磁场效应进行了数值计算.这一三维电流体系在中低纬度地面产生的磁场呈现出特殊的纬度分布:X分量几乎不随纬度变化,Y分量随纬度增高近似呈线性变化.这些特征明显不同于对称环电流的磁场分布特征（X∝ cosφ,φ是纬度,Y=0）,也不同于DP2、S_q、L等电流体系的磁场分布特征.利用这一特征我们可以从地磁台子午链观测到的磁场扰动中分离出PRFI电流系的贡献.用1989年3月磁暴的实例检验了上述模型,观测结果与理论结果符合得很好.分析结果还表明,最大的D_st指数并不一定对应着最强的对称环电流.     

伴随着高速流的磁尾电流片拍动观测研究

磁尾电流片在磁尾动力学过程中起着重要作用.卫星观测表明磁尾电流片经常处于拍动状态.但磁尾电流片拍动的特性和产生机制至今仍然没有被完全弄清楚.本文主要利用欧洲空间局Cluster卫星数据,研究一个伴随高速离子流的电流片拍动事件.该电流片拍动事件具有很强的周期性.拍动的周期约是2min,磁场振荡幅度约为20nT.能量电子和离子的通量具有周期性增强和减弱的特征.电流密度X和Y分量也具有周期性的振荡,并且振荡周期与磁场振荡周期一致.通过对粒子流速矢量与电流矢量的分析,发现粒子运动具有涡旋的特征.因此可以推断,该磁尾电流片的拍动不是由磁尾等离子体片高速流产生的,而是与局地等离子体不稳定性有关.     

通量传输事件磁场的理论计算及其与观测的比较

根据涡旋诱发重联理论,对通量传输事件（FTEs）磁场分布特性作了计算.结果表明,卫星测到的FTEs的不同磁场分布形态,是取决于通量管的运动方向及卫星穿越通量管的部位.在北半球,当通量管由低纬向高纬（由南向北）直向运动时,不论卫星通过什么部位,绝大多数情况下观测到先正后负的B_x,变化（即正FTE）,个别部位观测到先负后正的B_x变化（即反FTE）;B_z是单峰分布形式,表现为V型、倒V型或是U型和倒U型.当通量管在x方向有正或负速度分量即斜向运动时,大部分部位测到的B_x呈不规则变化,B_z表现为双极分布.与61个FTEs的观测实例作了对比,理论计算与观测符合得较好.     

Multiple Flux Rope Events at the High-Latitude Magnetopause: Cluster/Rapid Observation on 26 January, 2001

Cluster measurements of the cusp and high latitude magnetopause boundary on 26 January, 2001 confirm that the cusp is a dynamic region full of energetic charged particles and turbulence. An energetic ion layer at high-latitudes beyond and adjacent to the duskside magnetopause exists when the Interplanetary Magnetic Field (IMF) has a southward orientation. Multiple energetic ion flux bursts were observed in the energetic ion layer. Each energetic ion flux burst was closely related to a magnetic flux rope. The axes of the flux ropes lie in the direction pointing duskward/tailward and somewhat upward. An intense axis-aligned current flows inside the ropes, with the current density reaching ∼10⁻⁸ A/m². The main components of the energetic ions are protons, helium and CNO ions, which originate from the magnetosphere, flowing out into the magnetosheath along the axis of the flux ropes. The velocity of the magnetosheath thermal plasma relative to the deHoffman-Teller (DHT) frame is found to be basically along the axis of the flux ropes also, but towards the magnetosphere. These flux ropes seem to be produced somewhere away via magnetic reconnection and move at similar DHT velocities passing over the spacecraft. These observations further confirm that the high-latitude magnetopause boundary region plays an important role in the solar wind-magnetopause coupling.     

Skin-layer of the eruptive magnetic flux rope in large solar flares

The analysis of observations of large solar flares made it possible to propose a hypothesis on existence of a skin-layer in magnetic flux ropes of coronal mass ejections. On the assumption that the Bohm coefficient determines the diffusion of magnetic field, an estimate of the skin-layer thickness of ~10⁶ cm is obtained. According to the hypothesis, the electric field of ~0.01–0.1 V/cm, having the nonzero component along the magnetic field of flux rope, arises for ~5 min in the surface layer of the eruptive flux rope during its ejection into the upper corona. The particle acceleration by the electric field to the energies of ~100 MeV/nucleon in the skin-layer of the flux rope leads to their precipitation along field lines to footpoints of the flux rope. The skin-layer presence induces helical or oval chromospheric emission at the ends of flare ribbons. The emission may be accompanied by hard X-ray radiation and by the production of gamma-ray line at the energy of 2.223 MeV (neutron capture line in the photosphere). The magnetic reconnection in the corona leads to a shift of the skin-layer of flux rope across the magnetic field. The area of precipitation of accelerated particles at the flux-rope footpoints expands in this case from the inside outward. This effect is traced in the chromosphere and in the transient region as the expanding helical emission structures. If the emission extends to the spot, a certain fraction of accelerated particles may be reflected from the magnetic barrier (in the magnetic field of the spot). In the case of exit into the interplanetary space, these particles may be recorded in the Earth’s orbit as solar proton events.     

The self-similar,non-linear evolution of rotating magnetic flux ropes

We study, in the ideal MHD approximation, the non-linear evolution of cylindrical magnetic flux tubes differentially rotating about their symmetry axis. Our force balance consists of inertial terms, which include the centrifugal force, the gradient of the axial magnetic pressure, the magnetic pinch force and the gradient of the gas pressure. We employ the “separable” class of self-similar magnetic fields, defined recently. Taking the gas to be a polytrope, we reduce the problem to a single, ordinary differential equation for the evolution function. In general, two regimes of evolution are possible; expansion and oscillation. We investigate the specific effect rotation has on these two modes of evolution. We focus on critical values of the flux rope parameters and show that rotation can suppress the oscillatory mode. We estimate the critical value of the angular velocity _crit, above which the magnetic flux rope always expands, regardless of the value of the initial energy. Studying small-amplitude oscillations of the rope, we find that torsional oscillations are superimposed on the rotation and that they have a frequency equal to that of the radial oscillations. By setting the axial component of the magnetic field to zero, we study small-amplitude oscillations of a rigidly rotating pinch. We find that the frequency of oscillation is inversely proportional to the angular velocity of rotation ; the product being proportional to the inverse square of the Alfvén time. The period of large-amplitude oscillations of a rotating flux rope of low beta increases exponentially with the energy of the equivalent 1D oscillator. With respect to large-amplitude oscillations of a non-rotating flux rope, the only change brought about by rotation is to introduce a multiplicative factor greater than unity, which further increases the period. This multiplicative factor depends on the ratio of the azimuthal speed to the Alfvén speed. Finally, considering interplanetary magnetic clouds as cylindrical flux ropes, we inquire whether they rotate. We find that at 1 AU only a minority do. We discuss data on two magnetic clouds where we interpret the presence in each of vortical plasma motion about the symmetry axis as a sign of rotation. Our estimates for the angular velocities suggest that the parameters of the two magnetic clouds are below critical values. The two clouds differ in many respects (such as age, bulk flow speed, size, handedness of the magnetic field, etc.), and we find that their rotational parameters reflect some of these differences, particularly the difference in age. In both clouds, a rough estimate of the radial electric field in the rigidly rotating core, calculated in a non-rotating frame, yields values of the order mV m⁻¹.     

Model of quiescent prominence with the helical magnetic field

A new exact analytical solution of the magnetohydrostatic problem describes the equilibrium of a solitary, dense-cool solar filament maintained against the gravity by magnetic force in hot solar corona at heights up to 20–40 Mm. The filament is assumed to be uniform along the axis (the translation symmetry). The magnetic field of the filament has the helical structure (magnetic flux rope) with a typical strength of a few Gauss in the region of minimal temperature (about 4000 K). The model can be applied to the quiescent prominence of both normal and inverse magnetic polarity.

     

New method for determining central axial orientation of flux rope embedded within current sheet using multipoint measurements

A new method for determining the central axial orientation of a two-dimensional coherent magnetic flux rope (MFR) via multipoint analysis of the magnetic-field structure is developed. The method is devised under the following geometrical assumptions: (1) on its cross section, the structure is left-right symmetric; (2) the projected structure velocity is vertical to the line of symmetry. The two conditions can be naturally satisfied for cylindrical MFRs and are expected to be satisfied for MFRs that are flattened within current sheets. The model test demonstrates that, for determining the axial orientation of such structures, the new method is more efficient and reliable than traditional techniques such as minimum-variance analysis of the magnetic field, Grad-Shafranov (GS) reconstruction, and the more recent method based on the cylindrically symmetric assumption. A total of five flux transfer events observed by Cluster are studied using the proposed approach, and the application results indicate that the observed structures, regardless of their actual physical properties, fit the assumed geometrical model well. For these events, the inferred axial orientations are all in excellent agreement with those obtained using the multi-GS reconstruction technique.     

Dependence of the Occurrence of Coronal Mass Ejections on the Initial Length of the Eruptive Prominence

Geomagnetism and Aeronomy - A model of the eruption of a magnetic flux rope with ends rigidly fixed in the photosphere is analyzed. Long and short flux ropes exhibit different scenarios of...