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

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

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

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

Invariant Modulation of IMF Clock Angle on the Solar Wind Energy Input into the Magnetosphere

By use of the global PPMLR Magnetohydrodynamics (MHD) model, a serial of quasisteady- state numerical simulations were conducted to examine the modulation property of the interplanetary magnetic field clock angle θ on the solar wind energy input into the magnetosphere. All the simulations can be divided into seven groups according to different criteria of solar wind conditions. For each group, 37 numerical examples are analyzed, with the clock angle varying from 0° to 360° with an interval of 10°, keeping the other solar wind parameters (such as the solar wind number density, velocity, and the magnetic field magnitude) unchanged. As expected, the solar wind energy input into the magnetosphere is modulated by the IMF clock angle. The axisymmetrical bell-shaped curve peaks at the clock angle of 180°. However, the modulation effect remains invariant with varying other solar wind conditions. The function form of such an invariant modulation is found to be sin(θ/2)2.70 + 0.25.      

Real-time global MHD simulation of the solar wind interaction with the earth’s magnetosphere

We have developed a real-time global MHD (magnetohydrodynamics) simulation of the solar wind interaction with the earth’s magnetosphere. By adopting the real-time solar wind parameters and interplanetary magnetic field (IMF) observed routinely by the ACE (Advanced Composition Explorer) spacecraft, responses of the magnetosphere are calculated with MHD code. The simulation is carried out routinely on the super computer system at National Institute of Information and Communications Technology (NICT), Japan. The visualized images of the magnetic field lines around the earth, pressure distribution on the meridian plane, and the conductivity of the polar ionosphere, can be referred to on the web site (http://www2.nict.go.jp/y/y223/simulation/realtime/).The results show that various magnetospheric activities are almost reproduced qualitatively. They also give us information how geomagnetic disturbances develop in the magnetosphere in relation with the ionosphere. From the viewpoint of space weather, the real-time simulation helps us to understand the whole image in the current condition of the magnetosphere. To evaluate the simulation results, we compare the AE indices derived from the simulation and observations. The simulation and observation agree well for quiet days and isolated substorm cases in general.     

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.     

Solar wind and its interaction with the magnetosphere: Measured parameters

The earth's magnetosphere absorbs only a minor fraction (≈ 10^?3) of the incident solar wind energy. Variations of the solar wind can often cause lively reactions in the earth's close environment. However, the physical mechanisms involved are not yet understood. It appears now that the combined action of the solar wind momentum flux, the direction of the interplanetary magnetic field as well as its fluctuations might play the dominant role. The behaviour of these parameters is governed in some characteristic way by the solar wind stream structure which reflects the condition of the solar corona and its magnetic field topology. Transients in the sun's atmosphere associated with solar activity cause reactions in the interplanetary medium which also show some typical, though very different, signatures. Taking into account the interdependence of the solar wind parameters in context with the underlying solar phenomena, we may be able to pinpoint the mechanism which controls the action of the solar wind on the magnetosphere.     

行星际磁场北向分量触发磁层扰动的可能性

本文用典型事件和统计分析论证了行星际磁场北向分量触发地磁扰动的可能性.给出了行星际电场E_y分量对磁层大尺度对流电场E_M的耦合系数.从耦合方程出发讨论了磁层对太阳风作用的响应,证明磁层不是起半波整流器作用,从而表明经典的重联理论应有所修改.      

Bow shock: Power aspects

It is clear that the primary energy source for magnetospheric processes is the solar wind, but the process of energy transfer from the solar wind into the magnetosphere, or rather, to convecting magnetospheric plasma, appears to be rather complicated. Bow shock is a powerful transformer of the solar wind kinetic energy into the gas dynamic and electromagnetic energy. A jump of the magnetic field tangential component at front crossing means that the front carries an electric current. The solar wind kinetic energy partly transforms to gas kinetic and electromagnetic energy during its passage through the bow shock front. The transition layer (magnetosheath) can use part of this energy for accelerating of plasma, but can conversely spend part its kinetic energy on the electric power generation, which afterwards may be used by the magnetosphere. Thereby, transition layer can be both consumer (sink) and generator (source) of electric power depending upon special conditions. The direction of the current behind the bow shock front depends on the sign of the IMF Bz-component. It is this electric current which sets convection of plasma in motion.     

Statistical properties of the IMF clock angle in the solar wind with northward and southward interplanetary magnetic field based on ACE observation from 1998 to 2009: Dependence on the temporal scale of the solar wind

Utilizing ACE satellite observations from 1998 to 2009, we performed the elaborate study on the properties of the clock angle θ_CA (arctan(By/Bz) (?90° to 90°) of the interplanetary magnetic field (IMF) in the solar wind at 1?AU. The solar wind with northward IMF (NW-IMF) and southward IMF (SW-IMF) are analyzed, independently. Statistical analysis shows that the solar wind with SW-IMF and NW-IMF has similar properties in general, including their durations, the IMF Bz and By components, and the IMF θ_CA. Then, the solar wind with NW-IMF (SW-IMF) is classified into five different temporal scales according to the duration of the NW-IMF (SW-IMF), i.e., very-short wind of 10–30?min, short-scale wind of 0.5–1?h, moderate-scale wind of 1–3?h, long-scale wind of 3–5?h, and super-long wind >5?h. Our analysis reveals that the IMF θ_CA has a distinct decrease with increase of the temporal scale of the solar wind. Next, the solar wind is classified into two groups, i.e., the high-speed solar wind (>450?km/s) and the low-speed solar wind (<450?km/s). Our analysis indicates that the IMF θ_CA depends highly on the solar wind speed. Statistically, high-speed solar wind tends to have larger IMF θ_CA than low-speed solar wind. The evolutions of the solar wind and IMF with the solar activity are further studied, revealing no clear solar variation of the IMF θ_CA. Finally, we analyze the monthly variation of the IMF θ_CA. Superposed epoch result strongly suggests the seasonal variation of the IMF θ_CA.     

Solar wind dependence of the Pc1 wave activity

There are a host of factors influencing the excitation of Pc1 geomagnetic pulsations, which are ULF waves in the frequency range between 0.2 and 5 Hz. We have studied carefully the dependence of the pearl-type Pc1 activity at Sodankylä, Finland (L = 5.1) on the plasma density N in front of the magnetosphere, the bulk velocity V of the solar wind, and the intensity B of the IMF. The result is as follows: high values of N and reduced values of V are favorable to appearance of Pc1, whereas the dependence of Pc1 activity on B is practically absent. We also show that the probability of Pc1 occurrence decreases with the interplanetary electric field, and increases with solar wind impact pressure and with the plasma to magnetic pressure ratio “beta”.     

The Earth’s magnetosphere response to interplanetary medium conditions on January 21–22, 2005 and on December 14–15, 2006

The Earth’s magnetosphere response to interplanetary medium conditions on January 21–22, 2005 and on December 14–15, 2006 has been studied. The analysis of solar wind parameters measured by ACE spacecraft, of geomagnetic indices variations, of geomagnetic field measured by GOES 11, 12 satellites, and of energetic particle fluxes measured by POES 15, 16, 17 satellites was performed together with magnetospheric modeling based in terms of A2000 paraboloid model. We found the similar dynamics of three particle populations (trapped, quasi-trapped, and precipitating) during storms of different intensities developed under different external conditions: the maximal values of particle fluxes and the latitudinal positions of the isotropic boundaries were approximately the same. The main sources caused RC build-up have been determined for both magnetic storms. Global magnetospheric convection controlled by IMF and substorm activity driven magnetic storm on December 14–15, 2006. Extreme solar wind pressure pulse was mainly responsible for RC particle injection and unusual January 21, 2005 magnetic storm development under northward IMF during the main phase.     

Magnetospheric Boundary Layer Structure and Dynamics as Seen From Cluster and Double Star Measurements

In this review, we discuss the structure and dynamics of the magnetospheric Low-Latitude Boundary Layer (LLBL) based on recent results from multi-satellite missions Cluster and Double Star. This boundary layer, adjacent to the magnetopause on the magnetospheric side, usually consists of a mixture of plasma of magnetospheric and magnetosheath origins, and plays an important role in the transfer of mass and energy from the solar wind into the magnetosphere and subsequent magnetospheric dynamics. During southward Interplanetary Magnetic Field (IMF) conditions, this boundary layer is generally considered to be formed as a result of the reconnection process between the IMF and magnetospheric magnetic field lines at the dayside magnetopause, and the structure and plasma properties inside the LLBL can be understood in terms of the time history since the reconnection process. During northward IMF conditions, the LLBL is usually thicker, and has more complex structure and topology. Recent observations confirm that the LLBL observed at the dayside can be formed by single lobe reconnection, dual lobe reconnection, or by sequential dual lobe reconnection, as well as partially by localized cross-field diffusion. The LLBL magnetic topology and plasma signatures inside the different sub-layers formed by these processes are discussed in this review. The role of the Kelvin-Helmholtz instability in the formation of the LLBL at the flank magnetopause is also discussed. Overall, we conclude that the LLBL observed at the flanks can be formed by the combination of processes, (dual) lobe reconnection and plasma mixing due to non-linear Kelvin-Helmholtz waves.      

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.     

Laboratory experiments and space phenomena

Two types of convection were observed in the laboratory model of the magnetosphere: viscous convection and convection due to field lines common to both the magnetosphere and artificial solar wind. With a southward field component in the solar wind, convection from the Sun is observed in the polar cap, while with a large northward component, convection is directed toward the Sun. Merging of the field lines occurs in the cleft. With the southward component, a visor appears in front of the magnetosphere boundary. The decay of the visor into small magnetic structure is observed. The formation of an induced magnetosphere with a magnetic tail is shown in the experiments of the simulated conditions near non-magnetic bodies with a plasma shell (Venus, comets). A combined induced-intrinsic magnetosphere also was investigated.     

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

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

Geomagnetic activity driven by solar wind turbulence

Within the framework of the solar wind—magnetosphere coupled system, intense perturbations in the solar wind, causing geomagnetic storms and substorms, have been widely studied by means of the so-called coupling parameters. However, remarkable variations in the geomagnetic field occur even in absence of such perturbations. In those conditions, solar wind MHD turbulence might have a role. Recent results have shown that solar wind turbulence can be described not only as a mixture of inward and outward stochastic Alfvénic fluctuations, but includes also advected structures, dominated by an excess of magnetic energy.     

Effect of the Interplanetary Electric Field on the Magnetopause From Global MHD Simulations

The north-south component B_z of the Interplanetary Magnetic Field (IMF) and solar wind dynamic pressure P_d are generally treated as the two main factors in the solar wind that determine the geometry of the magnetosphere. By using the 3D global MHD simulations, we investigate the effect of the Interplanetary Electric Field (IEF) on the size and shape of magnetopause quantitatively. Our numerical experiments confirm that the geometry of the magnetopause are mainly determined by P_dB_z, as expected. However, the dawn-dusk IEFs have great impact on the magnetopause erosion because of the magnetic reconnection, thus affecting the size and shape of the magnetopause. Higher solar wind speed with the same B_z will lead to bigger dawn-dusk IEFs, which means the higher reconnection rate, and then results in more magnetic flux removal from the dayside. Consequently, the dayside magnetopause moves inward and flank magnetopause moves outward.      

Akasofu能量耦合函数的数值检验

从太阳风-磁层能量耦合的普遍表达式出发,用34天连续的太阳风观测资料对电磁耦合机制进行了数值检验.结果表明,只有当行星际磁场有南向分量的时候,电磁耦合机制才能近似表示太阳风-磁层能量耦合过程.此时,能量输入率可以表示成p=CB_T2/3V5/3n1/3sin4（θ/2）这个函数与Akasofu能量耦合函数ε=VB2l₀2sin4（θ/2）有一定差别,但与Murayama和Hakamada,Svalgaard,Holzer和Slavin等人的结果一致.本文对影响能量耦合函数计算的几个问题从原始资料、处理方法及物理机制上进行了讨论.      

Unusually strong magnetic fields in Titan’s ionosphere: T42 case study

Observations of unusually large magnetic fields in the ionosphere indicate periods of maximum stress on Titan’s ionosphere and potentially of the strongest loss rates of ionospheric plasma. During Titan flyby T42, the observed magnetic field attained a maximum value of 37 nT between an altitude of 1200 and 1600 km, about 20 nT stronger than on any other Titan pass and close to five times greater in magnetic pressure. The strong fields occurred near the corotation-flow terminator rather than at the sub-flow point, suggesting that the flow which magnetized the ionosphere was from a direction far from corotation and possibly towards Saturn. Extrapolation of solar wind plasma conditions from Earth to Saturn using the University of Michigan MHD code predicts an enhanced solar wind dynamic pressure at Saturn close to this time. Cassini’s earlier exits from Saturn’s magnetosphere support this prediction because the Cassini Plasma Spectrometer instrument saw a magnetopause crossing three hours before the strong field observation. Thus it appears that Titan’s ionosphere was magnetized when the enhanced solar wind dynamic pressure compressed the Saturnian magnetosphere, and perhaps the magnetosheath magnetic field, against Titan. The solar wind pressure then decreased, leaving a strong fossil field in the ionosphere. When observed, this strong magnetic flux tube had begun to twist, further enhancing its strength.     

IMF北向时太阳风粒子向磁层输运的试验粒子模拟研究

以ACE卫星实时观测数据驱动的全球磁流体模拟为背景场,选取2003年10月22－24日行星际磁场（IMF）持续北向的事件,使用试验粒子方法,对太阳风粒子向磁层输运的过程进行模拟研究,分析北向IMF下太阳风粒子注入磁层过程中粒子在磁层内的空间分布和时间演化特征。IMF北向期间,进入环电流区域的粒子在晨侧区域的密度大于昏侧,且晨侧的粒子分布范围更广。背阳面磁鞘中的太阳风粒子可以通过低纬边界层进入磁层,但很难通过南北侧磁层顶进入磁层。进入磁尾的太阳风粒子聚集形成冷而密的等离子体片（CDPS）,模拟中CDPS的空间分布和密度大小与观测数据符合。在IMF长时间北向期间,磁尾的粒子数量呈现随时间增长的趋势,并存在约20 min的小幅度准周期变化和约5~6 h的较大幅度的准周期变化。