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

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

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

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

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

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

行星际起伏向磁层顶的输运

时间尺度为分钟数量级的太阳风速度和行星际磁场大幅度扰动实际上始终存在于行星际空间的。这些扰动一直传输到紧贴磁层边界面外侧的区域。它们在磁鞘等离子体和磁层顶的相互作用过程中可能起很重要的作用。行星际起伏中的磁场分量在通过地球弓激波时首先经历一次跳跃,然后一部分扰动被带到磁层边界面处。在边界面附近磁场扰动幅度被大大地放大了。弓激波上游的太阳风条件控制了放大因子。本文所作的数值模拟研究结果表明,如果上游有大幅度的扰动,在边界面附近就有大幅度的Alfven起伏的磁场分量。当上游磁场接近垂直于日地联线时,放大因子变得相当大,而且放大因子随上游的等离子体β值和/或Alfven马赫数的增加而增加。上游各向异性对放大因子的影响不大。在磁层边界附近存在大幅度起伏表明这里不存在稳定的片流。      

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

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

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

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

1994年2月21日行星际激波引起的磁暴

利用Imp-8,Geotail和Goes-6等卫星资料,研究了1994年2月21日0900UT到达地球磁层的行星际激波引起的磁暴期间,从太阳风向磁层传输能量的有关问题.结果指出:（1）南向行星际磁场（IMF）的长持续时间不是太阳风向磁层输能的必要条件.南北振荡的,较强IMF也能产生显著的能量传输;（2）行星际扰动磁场通过弓激波和磁层顶后扰动磁能增加,增幅将近5倍;（3）在磁层内扰动磁场的Bz分量在1×10-4Hz附近显著被吸收.这一低频扰动磁场可能是磁暴期间导致氧离子和质子等环电流粒子向内扩散并被加速的原因之一.     

磁云边界层中的重联慢激波观测分析

在Petschek模型中，排空区边界处的一对慢激波是能量耗散的重要机制.已有大量行星际空间的Petschek型磁场重联排空区观测事件被报道，但是只有少量的排空区边界处观测到了慢激波.针对一例位于磁云边界层中的Petschek型磁场重联排空区观测事件，在排空区靠近磁云一侧边界处证认了一例慢激波.激波跃变层两侧的磁场和等离子体参数满足Rankine-Hugoniot关系，且激波上下游的中间马赫数均小于1，上游的慢马赫数为2.94（>1），下游的慢马赫数为0.65（<1），符合慢激波的观测特征.磁云内部的等离子体β值很低，局地阿尔芬速度高，同时磁云边界层中可能发生丰富的磁场重联活动，这可能是磁云前边界处慢激波形成的原因.      

磁层顶旋转间断的稳定性

本文应用三层模式和ISEE卫星观测资料,讨论了磁层顶旋转间断的稳定性。结果表明:(1)在磁层顶旋转间断中可以激发一种不稳定性。随着波数k增大,不稳定性增长率也将增加。(2)当行星际磁场为北向时,磁层顶旋转间断是稳定的;当行星际磁场逐渐变为南向时,不稳定性增长率将迅速增加。(3)当太阳风速度较大时,不稳定性增长率相应地也较大。(4)当行星际磁场为南向时,随着行星际磁场与磁层顶切平面交角的增大,不稳定性的增长率也迅速增加。      

太阳风——磁层相互作用的磁流体力学数值模拟研究

磁层位于地球空间的最外层, 太阳风与磁层的相互作用是空间天气变化因果链中承上启下的关键环节, 是揭示地球空间天气基本规律的关键科学课题. 地球空间由于时变、多成分、多自由度的关联相互作用,使得传统的理论分析变得非常困难. 数值模拟作为近几十年发展起来的一个新的研究手段,对地球空间的理论和应用研究产生了深刻的影响. 国际上磁层的全球MHD数值模拟工作开始于20世纪70年代末, 最初的研究局限于二维空间. 由于磁层内在的三维特性, 20世纪80年代三维MHD数值模拟工作兴起. 本文简要说明了三维全球磁层MHD (磁流体力学)研究的特点及现状, 给出了三维全球磁层模型的基本框架, 综述了行星际激波与磁层相互作用、大尺度电流体系、重联电压和越极电位、磁层顶K-H不稳定性等方面的太阳风--磁层相互作用的MHD数值模拟的研究进展.      

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.     

Small-scale deformation of the bow shock

The prediction of the bow shock location is a proof of our understanding of the processes governing the solar wind – magnetosphere interaction. However, the models describing the bow shock location as a function of upstream parameters are based on a statistical processing of bow shock crossings observed by a single spacecraft. Such crossings locate the bow shock in motion, i.e., in a non-equilibrium state and this fact can be a source of significant errors. We have carefully analyzed a long interval of simultaneous observations of the bow shock and magnetopause and another interval of bow shock observations at two well-separated points. Our results suggest that often a small-scale deformation of the bow shock front due to magnetosheath fluctuations is the most appropriate interpretation of observations. Since the low-frequency magnetosheath variations exhibit largest amplitudes, a simultaneous bow shock displacement over a distance of 10–15 R_E can be observed. We suggest that bow shock models can be probably improved if the tilt angle would be implemented as a parameter influencing the bow shock location in high latitudes.     

Electrodynamic Coupling in the Solar Wind-Magnetosphere-Ionosphere System

This paper presents a brief summary of our recent work based on global MHD simulations of the Solar wind-Magnetosphere-Ionosphere (SMI) system with emphasis on the electrodynamic coupling in the system. The main conclusions obtained are summarized as follows. (1) As a main dynamo of the SMI system, the bow shock contributes to both region 1 Field-Aligned Current (FAC) and cross-tail current. Under strong interplanetary driving conditions and moderate Alfven Mach numbers, the bow shock's contribution may exceed more than fifty percent of the total of either region 1 or cross-tail currents. (2) In terms of more than 100 simulation runs with due southward Interplanetary Magnetic Field (IMF), we have found a combined parameter f = E_swP_swM_A-1/2 (E_sw, P_sw, and M_A are the solar wind electric field, ram pressure, and Alfven Mach number, respectively): both the ionospheric transpolar potential and the magnetopause reconnection voltage vary linearly with f for small f, but saturate for large f. (3) The reconnection voltage is approximately fitted by sin3/2θ_IMF/2, where θ_IMF is the IMF clock angle. The ionospheric transpolar potential, the voltage along the polar cap boundary, and the electric fields along the merging line however defined they may be, respond differently to θ_IMF, so it is not justified to take them as substitutes for the reconnection voltage.      

A stationary bow shock model for plasmas: The spherical blunt obstacle problem

We present an analytic model of a stationary bow shock which describes the interaction between a supermagnetosonic ambient wind and an obstacle with spherical-like frontal shape. We develop expressions for the bow shock’s geometry and the physical properties of the plasma sheath as functions of the upstream conditions. The solution is limited to magnetic fields parallel to the upstream velocity. The model allows to use any value of the upstream alfvenic and sonic Mach numbers and the polytropic index (γ$\gamma$), pointing out the influence of γ$\gamma$ for the magnetosheath compression and the bow shock shape. When both Mach numbers are small, the upstream magnetic field intensity affects also the bow shock shape. We compare our results with other models finding important consistencies. We also compare our results with in-situ data, we fund a reasonable qualitative agreement; however, it seems that our model underestimates the magnetosheath size.     

Plasma flow variations and energetic protons upstream of the earth’s bow shock: A statistical study

Foreshock is a special region located upstream of the Earth’s bow shock characterized by the presence of various plasma waves and fluctuations caused by the interaction of the solar wind plasma with particles reflected from the bow shock or escaping from the magnetosphere. On the other hand, foreshock fluctuations may modify the bow shock structure and, being carried through the magnetosheath, influence the magnetopause. During the years 1995–2000, the INTERBALL-1 satellite made over 10,000 hours of plasma and energetic particles measurements in the solar wind upstream of the Earth’s bow shock. We have sorted intervals according to the level of solar wind ion flux fluctuations and/or according to the flux of back-streaming energetic protons. An analysis of connection between a level of ion flux fluctuations and fluxes of high-energy protons and their relation to the IMF orientation is presented.     

Global Hybrid Simulation of Magnetic Reconnection in the Magnetosheath

A three-dimensional (3-D) global hybrid simulation is carried out for the generation and structure of magnetic reconnection in the magnetosheath due to interaction of an interplanetary Tangential Discontinuity (TD) with the bow shock and magnetosphere. Runs are performed for solar wind TDs possessing different initial half-widths. As the TD propagates through the bow shock toward the magnetopause, it is greatly narrowed by a two-step compression processes, a "shock compression' followed by a subsequent ``convective compression'. In cases with a relatively thin solar wind TD, 3-D patchy reconnection is initiated in the transmitted TD, forming magnetosheath flux ropes. Multiple components of ion particles are present in the velocity distribution in the magnetosheath merging, accompanied by ion heating. For cases with a relatively wide initial TD, a dominant single X-line appears in the subsolar magnetosheath after the transmitted TD is narrowed. A shock analysis is performed for the detailed structure of magnetic reconnection in the magnetosheath. Rotational Discontinuity (RD)/Time-Dependent Intermediate Shock (TDIS) are found to dominate the reconnection layer, which and some weak slow shocks are responsible for the ion heating and acceleration.      

Forward and reverse CIR shocks at 4–5 AU: Ulysses

In this article, we study fast shocks at CIR boundaries during an extended interval of 15 consecutive major high speed solar wind streams in 1992–1993. Ulysses was 4–5 AU from the sun. The Abraham-Schrauner shock normal method and the Rankine-Hugoniot relations were used to determine fast shock directions and speeds. Out of 33 potential CIR shocks, 14 were determined to be fast forward shocks (FSs) and 14 were fast reverse shocks (RSs). Of the remaining 5 events, 2 were forward waves and 3 were reverse waves. CIR edges at latitudes below ∼30^o were, for the most part, bounded by fast magnetosonic shocks. The forward shocks were generally quasi-perpendicular (average θ_nBo = 67^o). The reverse shocks were more oblique (average θ_nBo = 52^o), but they extended to all angles. Both FSs and RSs had magnetosonic Mach numbers ranging from 1 to 5 or 6. The average Mach numbers were 2.4 and 2.6 for FSs and RSs, respectively. The shock Mach numbers were noted to generally decrease with increasing latitude. The non-shock events or waves were noted to occur preferentially at high (∼−30° to −35°) heliolatitudes where stream-stream interactions were presumably weaker. These results are consistent with expectations, indicating the general accuracy of the Abraham-Schrauner technique.     

Dayside reconnection during IMF northward: A possible foreshock effect

The northward and southward orientation of the interplanetary magnetic field (IMF) is usually considered as providing the external boundary conditions in the solar wind interaction with the Earth's magnetopause but it is the magnetic field in the magnetosheath that interacts with the Earth's magnetic field. In this paper, we consider the possibility that the wave activity in the foreshock region may affect the magnetic field orientation in the magnetosheath with time scales that might be geomagnetically effective. If magnetosheath magnetic field becomes disturbed on plasma streamlines which are connected to the quasi-parallel bow shock and foreshock, the magnetic field orientation on the inner magnetosheath may differ significantly from the undisturbed IMF. We present a model of dayside reconnection which may occur when the IMF northward and illustrate its effects on the erosion of the magnetopause.