地球弓激波前的低频磁流体波

一、引言 在地球弓激波前存在着低频磁流体波。这种低频磁流体波是太阳风在地球弓激波上的反射粒子和太阳风粒子之间相互作用产生的。根据人造卫星的观测资料可以得到,在地球弓激波前,Pc3-4脉动频率范围内的低频磁流体波的主频率和行星际磁场强度     

内日球太阳风的激波加热

利用HeliosA,B飞船1974年至1980年的太阳风探测资料,分析了不同速度间隔太阳风质子温度径向变化指数在太阳不同活动期的变化,以及不同太阳活动期间内日球行星际激波强度分布的变化．结果指出,在0．3-1．0AU区间行星际激波可能是太阳风加热的一个重要因素,这一因素在太阳活动高年可能起着主要的作用．激波MHD数值模拟也从量的方面表明激波加热太阳风的有效性．     

电离层电导对地球磁层顶和舷激波尺度的影响

本文在如下假定下分析电离层电导对地球磁层顶和舷激波尺度的影响：（1）对电离层采用球壳近似,Pedersen电导Σ_P均匀,Hall电导为零;（2）地磁偶极矩处于正南方向,行星际磁场（IMF）只有南向分量（B_z<0）.磁层顶和舷激波的尺度分别由它们与GSE坐标系三个轴的交点,即日下点、晨昏侧翼点和南北顶点的地心距离表征.对给定的太阳风条件、B_z和Σ_P,通过三维全球MHD模拟获得系统的准定态.结果表明,在大约1～5 S范围内,Σ_P值显著影响磁层顶和舷激波的尺度,而在该范围之外则几乎没有影响.随着Σ_P的增加,磁层顶和舷激波整体向外扩张,前者的扩张程度低于后者,以至磁鞘区的范围扩大.磁层顶的侧翼点的位置随Σ_P的变化与B_z的幅度有关：在弱南向IMF情况下磁层顶的侧翼点随Σ_P的增加向内移动,而在强南向IMF情况下则向外移动.上述结果表明,在构建磁层顶和舷激波的经验模型时,有必要计入电离层电导的影响.     

行星际激波的传播及其相应的地磁效应

基于 2.5 维理想磁流体力学（Magnetohydrodynamic,MHD）方程组分析了行星际激波在日球层子午面内的传播过程及其相应的地磁效应.日球层电流片（Heliospheric Current Sheet,HCS）-日球层等离子体片（Heliospheric Plasma Sheet,HPS）对于行星际激波的传播具有一定的阻碍作用.当行星际激波相对于HCS 倾斜传播时,相对于扰动源位于HCS 异侧的激波强度较同侧的明显减弱.局地激波面的法线（或形状）对通过激波阵面的磁力线发生偏转的程度和方向起决定性作用.沿激波传播方向其为准平行激波,磁场偏转程度较小,而其两侧部分则为斜激波,磁场偏转程度较大.位于HCS-HPS 位置处的波前形成凹槽,磁力线偏转程度明显加强.行星际激波对磁场的偏转效应是其驱动地磁暴的重要机制,而且地磁效应的强度与地球相对于HCS 的角距离Δθ_p有明显关系.数值模拟结果表明：任何行星际激波,Δθ_p=0°处均无法形成较大强度的地磁效应;沿HCS 传播的行星际激波,地磁效应最强的区域位于HCS 两侧;相对于HCS 倾斜传播的行星际激波,地磁效应最强的区域位于HCS 异侧.     

地球远磁尾中的磁场重联——（Ⅱ）局部驱动力的影响

基于可压缩磁流体动力学模型,数值研究了由太阳风引起的局部驱动力对地球远磁尾中磁场重联的影响.结果表明,在远磁尾等离子体片中将发生强迫磁场重联,并形成磁岛和等离子体团.形成磁岛的特征时间很大于流动撕裂模不稳定性引起磁岛非线性饱和的特征时间.磁岛宽度随着磁Reynolds数S的增大而减小,随着尾瓣中等离子体压力与磁压之比值β_∞的降低而减小.认为太阳风引起的局部驱动力对地球远磁尾等离子体片中磁场重联的影响,可能不如流动撕裂模不稳定性那样显著.     

等离子体片离子与太阳风及地磁条件的关联研究

本文根据搭载于Cluster卫星的CIS/CODIF和RAPID仪器的观测数据,统计研究了等离子体片中的H⁺、O⁺离子在磁暴期间的时间变化特性,及其对太阳风条件的响应.观测结果表明：(1) 磁暴开始前,O⁺离子（0~40 keV）数密度保持在较低水平.随着磁暴的发展,O⁺数密度缓慢上升,其峰值出现在D_st极小值附近;H⁺离子(0~40 keV)数密度在磁暴开始之前的较短时间迅速增加并达到峰值,在磁暴开始之后迅速降低,并在整个主相和恢复相期间保持在相对较低水平.更高能量的离子则在磁暴开始后迅速增多,并在低能O⁺离子达到峰值之前达到峰值.因此我们推测磁暴初期从等离子体片注入环电流的主要是H⁺离子,主相后期O⁺离子可能扮演更为重要的角色.(2)在地磁活动时期,太阳风密度和动压强与等离子体片中的H⁺、O⁺数密度存在一定相关性.等离子体片中的H⁺离子对北向IMF B_z较为敏感,而IMF B_z南向条件下更有利于太阳风参数对等离子体片中O⁺数密度的影响.在地磁活动平静期,太阳风条件对等离子体片中的离子没有明显影响.     

太阳风压力系数的研究

利用全球磁流体力学(MHD)的模拟结果,研究了太阳风压力系数与上游太阳风参数和日下点磁层顶张角的相关性.在识别出日下点附近磁层顶位置后,通过拟合得到日下点附近的磁层顶张角.在考虑上游太阳风中的磁压和热压以及磁层顶外侧的太阳风动压的情况下,计算了太阳风压力系数.通过分析行星际磁场不同方向时太阳风动压在日地连线上与磁压和热压的转化关系,详细研究了太阳风参数和日下点磁层顶张角对太阳风压力系数的影响,得到以下相关结论:(1) 在北向行星际磁场较大(B_z≥5 nT)时,磁层顶外侧磁压占主导,南向行星际磁场时磁层顶外侧热压占主导;(2) 太阳风压力系数随着行星际磁场的增大而增大,随着行星际磁场时钟角的增大而减小;并且在行星际磁场大小和其他太阳风条件相同时,北向行星际磁场时的太阳风压力系数要大于南向行星际磁场时的;北向行星际磁场时,太阳风压力系数随着太阳风动压的增大而减小,南向行星际磁场时,太阳风压力系数随着太阳风动压的增大而增大;以上结论是对观测结果的扩展;(3) 最后,我们还发现太阳风压力系数随着日下点磁层顶张角的增大而增大.     

耀斑引发的激波初始速度对激波到达时间预测的影响

利用Hakamada-Akasofu-Fry运动学太阳风模型模拟了1981到1985卡林顿周的48个太阳耀斑事件激发的激波到达地球的时间.结果表明,对模式输入参数之一的激波初始速度进行调整,可以使模拟结果和实际观测基本一致.通过对发生在该时期多个事件的统计分析,分别得到日球表面东、西两半球耀斑爆发对应的激波初始速度调整...     

高密度极尖区的形成原因

极尖区是太阳风进入磁层的一个重要窗口,极尖区密度是反映这一物理过程的重要参量,通常情况下极尖区密度约为1~10 cm^-3,但有时卫星会观测到密度大于40 cm^-3的极尖区,本文称之为高密度极尖区.我们分析了Cluster卫星2001—2009年的观测数据,在470个极尖区穿越中找到28个高密度极尖区穿越事件并进行了统计研究,分析了高密度极尖区事件的形成原因,进而讨论了太阳风高效地进入极尖区的外部条件.结果表明：距正午的距离（｜MLT-12｜）较小,太阳风的密度高,低纬有磁层顶磁重联发生以及正偶极倾角都是观测到高密度极尖区事件的有利条件,并且当同时满足上述4个条件时,高密度极尖区事件发生率为100%;而低纬磁层顶磁重联以及大的正偶极倾角被认为是太阳风高效地进入极尖区的重要条件.这些研究结果有助于我们更进一步地理解太阳风进入极尖区的物理机制.     

冕洞与太阳风高速流关系的统计研究:到达时间、持续时间和峰值强度

冕洞是太阳风高速流的源区.当冕洞出现在中低纬区域时,太阳风高速流会扫过地球并引发地球空间环境扰动,如地磁暴和高能电子暴等.在太阳活动周下降年和低年,这种类型的扰动占据主导地位.因此,冕洞高速流的到达时间、峰值时间、峰值强度和持续时间等,是空间天气预报的重要内容.本文基于2010年5月到2016年12月的SDO/AIA太阳极紫外图像以及1AU处ACE和WIND卫星的太阳风观测数据,确定了160个冕洞-太阳风高速流事件,定量计算了他们的特征参数,包括冕洞与太阳风高速流的开始时间、峰值时间、峰值强度和结束时间,分析了各个特征参数的分布规律,对冕洞-高速流之间的关系进行了统计研究,并提出了一种新的预报方法,为基于冕洞成像观测的太阳风高速流的精准预报提供了依据.     

Analytical description of electric currents in the magnetosheath region

Volume currents in the magnetosheath region are calculated within the framework of a new analytical model. Magnetic field structure in the region is found, satisfying boundary conditions on the bow shock and the magnetopause, and then volume currents are calculated using the Maxwell equation. Surface bow shock and magnetopause currents are calculated, too. Free parameters of the model are interplanetary magnetic field, Mach number of the solar wind flow, distances to the bow shock and to the magnetopause, and field compression at the magnetopause.     

地球弓激波的旋转非对称性

通过对太阳风-磁层-电离层系统的全球MHD模拟,研究地球弓激波相对日地连线的旋转非对称性．模拟限于太阳风速度沿日地连线、地球磁偶极矩和行星际磁场（IMF）与日地连线垂直的简单情况．模拟结果表明,即便对于IMF强度为零的情况,弓激波相对日地连线也不具备旋转对称性质：终端面（晨昏子午面）及其向阳侧的弓激波截线的东西宽度大于南北宽度（约9%~11％）,终端面尾侧的弓激波截线东西宽度小于南北宽度（约8％）．在存在IMF的情况下,弓激波的位形同时受到磁层顶的形状和快磁声波速度各向异性的影响．磁层顶向外扩张并沿IMF方向拉伸,且其扩张和拉伸程度随IMF由北转南而增强．在磁鞘中,垂直于磁场方向的快磁声波速度高于平行方向．因此,磁层顶拉伸方向与快磁声波速度最大方向垂直,它们对弓激波位置的效应恰好相反;弓激波的最终形状取决于何种效应占据主导地位．对于终端面尾侧,快磁声波速度的各向异性起主导作用,弓激波截线沿IMF垂直方向的宽度大于平行方向．对于终端面及其向阳侧,弓激波截线的形状与IMF取向有关：在准北向或晨昏向IMF情况下,弓激波截线沿IMF垂直方向的宽度仍大于平行方向;在准南向IMF情况下,弓激波截线沿IMF垂直方向的宽度小于平行方向的．鉴于弓激波形状同IMF取向之间的密切关系,我们提议以IMF为基准方向,提取弓激波截线的平行半宽度R_b∥和垂直半宽度R_b⊥作为尺度参数．这些尺度参数和通常引入的弓激波截线的东西半宽度y_b和南北半宽度z_b相比,更为合理地表征了弓激波的几何性质．模拟结果表明,在终端面上,y_b/z_b和R_b∥/R_b⊥在IMF各向同性取向下的统计平均值均低于1,与观测得到的结论一致.     

Interaction of the solar wind with Mars from Mars Global Surveyor MAG/ER observations

The observations of the magnetometer/electron reflectometer (MAG/ER) investigation onboard the Mars Global Surveyor (MGS) have greatly contributed to improve our understanding of the interaction of the solar wind with Mars. These observations established conclusively that a global dynamo-generated magnetic field does not exist at Mars, and that the interaction with solar wind is of the atmospheric type. This article reviews the most important results obtained from MGS MAG/ER on the study of two major features in the Mars solar wind interaction. The first feature is the occurrence of large-amplitude, highly coherent waves at the proton cyclotron frequency in the region upstream from the Martian bow shock. The second feature is the magnetic pileup boundary (MPB), a well-defined plasma boundary inside of which the planetary exospheric ions outnumber the solar wind ions. The study of these two elements is crucial to characterize the properties of the Martian exosphere. In addition, the occurrence of an MPB at comets and Venus reveals common processes to all these unmagnetized atmospheric bodies in spite of their different physical nature and characteristic scales.     

Effects of the interplanetary magnetic field clock angle on the shape of bow shock

Using the global magnetohydrodynamics(MHD) simulation model, we investigated the effects of the interplanetary magnetic field(IMF) clock angle on the shape of bow shock, including its rotational asymmetry and subsolar point. For general northward IMF( z component Bz 0), the rotational symmetry of the bow shock is broken by the effects of fast magnetosonic Mach number(Mms), and the cross-sectional line of the bow shock is an ellipse with the semi-major axis along the direction perpendicular to the IMF. The ratio or D-value between semi-major and semi-minor axis can be used to illustrate the extent of asymmetry of the bow shock. On the basis of the multiple parameters fitting, we obtain the changing relationship of both semi-axes with the clock angle and the distance away from the Earth. For general southward IMF(Bz 0), the cross sectional line of the bow shock is highly asymmetrical under the multiple effects of magnetopause and Mms. The effects of IMF clock angle on subsolar point depend mainly on the changing subsolar point of magnetopause as an obstacle. The distance of subsolar point of bow shock from the Earth increases with the increasing IMF clock angle for Bz 0, and decreases with the increasing IMF clock angle for Bz 0.     

南向行星际磁场事件与磁暴关系的研究

利用172-182年IMP-8飞船的太阳风观测资料和相应地磁活动性指数Dst和AE,研究了43个南向行星际磁场事件期间太阳风和磁层的耦合问题. 与这43个事件对应的地磁暴是中等的和强的磁暴(Dst＜-50nT). 结果表明:(1) 在43个事件中有11个(约占25.6髎)紧随激波之后,18个处于激波下游流场中(占42髎),其余14个(占33髎)和激波没有关连. 绝大多数事件都伴有太阳风动压和总磁场强度的增加;(2) 当行星际晨昏向电场强度EI＞-4mV/m时,只引起磁亚暴,对Dst指数没有明显影响. 仅当EI＜-5mV/m时,磁亚暴和磁暴才会同时出现;(3) 太阳风动压的增加会增强能量向环电流的输入,但不是密度和速度单独起作用,而是以PK=ρV2的组合形式影响能量的输入;(4) 虽然行星际磁场(IMF)南向分量BZ对太阳风和磁层的耦合起着关键作用,但IMF的BX和BY分量相对于BZ的大小对太阳风向磁层的能量传输也有一定影响. 当BX、BY相对BZ较大时能量耦合加强.     

The Influence of Parameters of the Interplanetary Medium and Magnetosheath Boundaries on the Correlation Coefficient between the Ion Flux Measured in the Solar Wind and the Magnetosheath

In this paper, the correlation coefficient between the ion fluxes in the solar wind and the magnetosheath is analyzed with the use of data of two satellites of the THEMIS mission and the THEMIS/Spektr-R satellites obtained in 2008 and 2011?2014, respectively. We have distinguished the conditions in which a high level of correlation between the measurements in the solar wind and the magnetosheath is observed, i.e., the correlation coefficient exceeds 0.7. As key factors, we consider both direct parameters of the solar wind, such as the density, the magnetic field magnitude, the magnetosonic Mach number, and the ratio β of the thermal pressure to the magnetic, and a more general factor—the type of large-scale structure of the solar wind. In addition, the effect of the satellite location in the magnetosheath relative to its boundaries—the bow shock and the magnetopause—on the correlation level is considered. It has been shown that, in roughly one third of cases, the plasma structures of the solar wind undergo a strong modification at the bow shock and in the magnetosheath, which results in a low correlation level corresponding to a correlation coefficient of less than 0.5; a high correlation level is observed in half of cases, i.e., the plasma structures are weakly disturbed. It has been determined that (1) the low correlation level in the magnetosheath behind quasi-perpendicular bow shock is more often observed near the magnetopause than in region just behind the bow shock, (2) the probability of observations of a high correlation level is independent of the profile shape of the quasi-perpendicular bow shock, and (3) the high correlation is more probable for the events corresponding to the solar wind of the Corotating Interaction Region (CIR) type than for those with the other solar wind types observed in the considered period.     

Analytical model of the near-Earth magnetopause according to the data of the Prognoz and Interball satellite data

Based on the magnetopause observations near the Earth by the Prognoz/Interball satellites in 1972–2000, the empirical model of this boundary has been proposed, and the magnetopause behavior at different parameters of the oncoming solar wind has been studied. For the first time, it has been detected that the Earth’s magnetopause is compressed by ∼5% in the direction perpendicular to the plane including the vectors of the solar wind velocity and IMF. At the same time, any dependence of the subsolar magnetopause position on the IMF B _z component has not been revealed in the Progrnoz/Interball data. The proposed magnetopause model can be used to model the position and shape of the near-Earth bow shock.     

Shape and position of Earth’s bow shock near-lunar orbit based on ARTEMIS data

Earth’s bow shock is the result of interaction between the supersonic solar wind and Earth’s magnetopause. However, data limitations mean the model of the shape and position of the bow shock are based largely on near-Earth satellite data. The model of the bow shock in the distant magnetotail and other factors that affect the bow shock, such as the interplanetary magnetic field (IMF) B_y, remain unclear. Here, based on the bow shock crossings of ARTEMIS from January 2011 to January 2015, new coefficients of the tail-flaring angle α of the Chao model (one of the most accurate models currently available) were obtained by fitting data from the middle-distance magnetotail (near-lunar orbit, geocentric distance -20R_E>X>-50R_E). In addition, the effects of the IMF B_y on the flaring angle α were analyzed. Our results showed that: (1) the new fitting coefficients of the Chao model in the middle-distance magnetotail are more consistent with the observed results; (2) the tail-flaring angle α of the bow shock increases as the absolute value of the IMF B_y increases. Moreover, positive IMF B_y has a greater effect than negative IMF B_y on flaring angle. These results provide a reference for bow shock modeling that includes the IMF B_y.