行星际激波导致的磁尾等离子片中ULF波动事件

磁层中的超低频（ULF）波动在太阳风和磁层之间的能量输运过程中具有重要作用.ULF波动主要发生在内磁层，且内磁层中ULF波动影响粒子的加速及沉降，而在夜侧磁层尤其是磁尾等离子片中观测到的ULF波动比较少.基于中国自主磁层探测卫星TC-1的观测数据，发现了两例行星际激波导致的磁尾中心等离子片中ULF波动事件，并发现这两例ULF事例都包含很强的环向模驻波分量，与以往THEMIS卫星报道的同类事件观测特征相符.根据ULF波的观测特征，分析了这两例ULF波动的可能触发机制.研究结果有助于深入理解磁层对行星际激波的全球响应.      

夜侧极区电离层对流对行星际激波的响应

利用超级双子极光雷达网（Super Dual Aurora Radar Network,SuperDARN）高频雷达、北半球IMAGE地磁台链以及南极中山站的极光观测数据,研究电离层对流对2012年7月14日一个行星际激波扰动事件的响应.在18:10UT行星际激波到达地球并与磁层相互作用触发地磁急始和磁层亚暴,SuperDARN雷达观测到北半球夜侧极区电离层对流显著增强,观测视野覆盖黄河站的Hankasalmi雷达观测到从激波到达地球至18:33UT,电离层F层出现剧烈扰动,雷达回波数明显增多,并出现局部对流速度反转现象.18:33UT之后,观测到F层出现三块速度高达600m·s-1的逆阳运动不规则体.而与Hankasalmi雷达地磁共轭的南半球Kerguelen雷达探测到的回波主要来自E层,回波数量几乎无变化,但是Kerguelen雷达观测视野内的中山站全天空光学成像仪观测到极光活动显著增强.南北半球夜侧电离层观测结果的差异,主要是由于它们分别处于极夜和极昼.      

2004年12月13日IMF北向期间极区亚暴电离层电动力学特征

利用KRM地磁反演方法, 结合北半球中高纬度地磁台站数据, 研究了2004年12月13日行星际磁场北向期间发生的亚暴事件, 极区电离层电动力学参量(电流矢量、等效电流函数以及电势)的分布特征. 结果表明, 在该亚暴膨胀相起始后, 午夜之前西向电集流急剧增强, 且等效电流体系表现为夜侧双涡, 同时伴随夜侧增强的南向电场. 由于极弱的直接驱动过程, 卸载过程引起的电离层效应得到清楚显示. 卸载过程在膨胀相期间起绝对主导性作用. 同时, 夜侧电导率的增强是电集流区域电流急剧增强的主要原因.      

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

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

地球弓激波的三维模拟

利用磁流体动力学(MHD)全球模拟结果,根据弓激波的跃变特性确定出弓激波位置,建立了一个新的综合考虑了快磁声马赫数、太阳风动压、行星际磁场强度以及磁层顶曲率半径的弓激波三维位型模型.将新模型与以往模型的模拟结果进行比较发现,新的弓激波全球模型结果可靠,解决了部分现有模型不能描述弓激波三维位型的问题.研究结果表明,在行星际磁场北向时,随着快磁声马赫数的增大,弓激波日下点距离减小,但是在行星际磁场南向时,快磁声马赫数的变化对弓激波日下点距离影响不大;弓激波位型在赤道面与子午面上存在明显的不对称性,而且随着行星际磁场的转向,这种非对称性也会发生相应改变;行星际磁场南向,B_z值较小时,子午面内弓激波位型已经不是简单的抛物线,出现了明显的类似于极尖区磁层顶的凹陷变化区.      

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

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

磁鞘等离子体密度统计建模

利用理想磁流体LFM模型的模拟数据,基于非参数统计方法对2004年11月14日03:00UT-07:00UT磁暴恢复相期间磁鞘等离子体平均密度进行建模.分析磁鞘等离子体平均密度与上游太阳风参数、行星际磁场参数及地磁扰动参数的统计关系,建立基于数据降维的经验模型.结果表明,电离层扰动强度因子、太阳风-磁层耦合强度因子和日地空间因果链耦合强度因子是影响磁鞘等离子体平均密度的三个主要方面.磁暴恢复相期间电离层上行离子是磁层环电流和磁尾等离子体的重要离子来源.建模分析过程表明,利用经验模型对空间物理过程开展建模,数据的严重多重共线性通常会导致模型的精度较差.而利用SIR和LPR建立的磁鞘等离子体平均密度随相关参数变化的经验模型可以有效解决该问题,具有较好的预测精度,统计特征显著.      

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

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

电离层电场对环电流变化响应的时延

本文讨论了强磁暴期间磁层环电流能量变化率与电离层电场变化之间的联系.STARE和SABRE雷达资料表明,电离层对环电流变化响应的主要特点是:(1)在磁地方时午后区,响应的时延达最大(约1—2小时),场强以指数形式增加;在其它时区内,无系统的增强过程,仅观测到较大的、有明显涨落的电场值.(2)STARE(70.2°N)和SABRE(65.8°N)测到的电场变化往往具有相反的趋势.(3)在STARE视场内,环电流能量变化达极大后,较低纬(70.2°N)上的电场值经常大于较高纬(71.8°N)上的值.分析结果表明,磁暴期间磁层-电离层耦合过程中,环电流起着重要作用.      

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

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

Propagation of interplanetary shocks into the Earth’s magnetosphere

This paper is devoted to the study of propagation of disturbances caused by interplanetary shocks (IPS) through the Earth’s magnetosphere. Using simultaneous observations of various fast forward shocks by different satellites in the solar wind, magnetosheath and magnetosphere from 1995 till 2002, we traced the interplanetary shocks into the Earth’s magnetosphere, we calculated the velocity of their propagation into the Earth’s magnetosphere and analyzed fronts of the disturbances. From the onset of disturbances at different satellites in the magnetosphere we obtained speed values ranging from 500 to 1300 km/s in the direction along the IP shock normal, that is in a general agreement with results of previous numerical MHD simulations. The paper discusses in detail a sequence of two events on November 9th, 2002. For the two cases we estimated the propagation speed of the IP shock caused disturbance between the dayside and nightside magnetosphere to be 590 km/s and 714–741 km/s, respectively. We partially attributed this increase to higher Alfven speed in the outer magnetosphere due to the compression of the magnetosphere as a consequence of the first event, and partially to the faster and stronger driving interplanetary shock. High-time resolution GOES magnetic field data revealed a complex structure of the compressional wave fronts at the dayside geosynchronous orbit during these events, with initial very steep parts (10 s). We discuss a few possible mechanisms of such steep front formation in the paper.     

Ballooning perturbations in the inner magnetosphere of the Earth: Spectrum,stability and eigenmode analysis

We investigate the generation of ballooning perturbations in the inner magnetosphere of the Earth in the dipole model of the geomagnetic field taking into account ionospheric boundary conditions. The ionosphere is considered as a thin layer with finite conductivity. The eigenmode spectrum is discrete and consists of Alfvén, slow magnetosonic, flute and incompressible modes. Their interaction depends on ionospheric conductivity. The decay rate is small in noon and night sectors and large in dawn and dusk sectors. The lowest stability threshold α/γ ≈ 4.25 is determined by flute modes.     

ISEE-1 and 2 observations of field-aligned currents in the distant midnight magnetosphere

Magnetic field measurements obtained in the nightside magnetosphere by the co-orbiting ISEE-1 and 2 spacecraft have been examined for signatures of field-aligned currents (FAC). Such currents are found on the boundary of the plasma sheet both when the plasma sheet is expanding and when it is thinning. Plasma sheet boundary layer current structure and substorm associated dynamics can be determined using the two spacecraft, although for slow traversals of the FAC sheet the spatial/temporal ambiguity is still an issue. We often find evidence for the existence of waves on the plasma sheet boundary, leading to multiple crossings of the FAC sheet. At times the boundary layer FAC sheet orientation is nearly parallel to the X-Z GSM plane, suggesting ‘protrusions’ of plasma sheet into the lobes. The boundary layer current polarity is, as expected, into the ionosphere in the midnight to dawn local time sector, and outward near dusk. Current sheet thicknesses and velocities are essentially independent of plasma sheet expansion or thinning, having typical values of 1500 km and 20–40 km/s respectively. Characteristic boundary layer current densities are about 10 nanoamps per square meter.     

How do interplanetary shock impact angles control the size of the geoeffective magnetosphere?

In this paper, we investigate temporal and spatial magnetosphere response to the impact of interplanetary (IP) shocks with different inclinations and speeds on the Earth’s magnetosphere. A data set with more than 500 IP shocks is used to identify positive sudden impulse (SI⁺) events as expressed by the SuperMAG partial ring current index. The SI⁺ rise time (RT), defined as the time interval between compression onset and maximum SI⁺ signature, is obtained for each event. We use RT and a model suggested by Takeuchi et al. (2002) to calculate the geoeffective magnetospheric distance (GMD) in the shock propagation direction as a function of shock impact angle and speed for each event. GMD is a generalization of the geoeffective magnetosphere length (GML) suggested by Takeuchi et al. (2002), defined from the subsolar point along the X line toward the tail. We estimate statistical GMD and GML values which are then reported for the first time. We also show that, similarly to well-known results for RT, the highest correlation coefficient for the GMD and impact angle is found for shocks with high speeds and small impact angles, and the faster and more frontal the shock, the smaller the GMD. This result indicates that the magnetospheric response depends heavily on shock impact angle. With these results, we argue that the prediction and forecasting of space weather events, such as those caused by coronal mass ejections, will not be accurately accomplished if the disturbances’ angles of impact are not considered as an important parameter within model and observation scheme capabilities.     

基于Cluster观测的磁尾场向电流在南北半球投影位置的分布

利用Cluster四颗卫星的磁场探测数据计算磁尾场向电流并投影到极区电离层,研究其投影位置在南北半球的分布规律,统计过程中去除了强磁暴（磁暴主相Dst<–100 nT）期间的场向电流事件。结果显示:磁尾场向电流事件在极区投影位置的纬度分布具有明显的南北半球不对称性,北半球为单峰结构,南半球为双峰结构。在北半球投影到较低纬度（<64°）的场向电流事件数目明显多于南半球,并且所能达到的最低纬度更低;在南半球投影到较高纬度（>74°）的场向电流事件数目明显多于北半球,并且所能达到的最高纬度更高。地磁平静条件下（|AL|<100 nT）,磁尾场向电流密度随磁地方时（MLT）呈递增趋势,这一结果与低高度卫星在极区对I区场向电流的探测结果符合很好。研究结果表明,磁尾场向电流投影位置的纬度分布呈现出明显的南北不对称性,这与南北半球磁尾场向电流的空间分布以及磁层中磁场结构具有密切关系。      

Some peculiarities of longitudinal distribution of proton fluxes at high latitudes

Dynamical features of proton fluxes at high and middle latitudes were studied based on data measured by Sun-synchronous low-altitude (800 km height) polar-orbiting three NOAA series satellites: POES 15, 16, and 17 during the geomagnetic storm on January, 21–22, 2005. Data from three satellites that passed the Northern hemisphere along different MLTs allow reconstructing the longitudinal distribution of the proton fluxes. Measurements of protons with energies of 30–80 keV and 80–240 keV (the ring current energy range) by 0- and 90-detectors were used to evaluate and compare the longitudinal asymmetry of proton flux distribution measured in the regions equatorward and poleward of the isotropic boundary. It was found that during all the phases of the geomagnetic storm distribution of the maximum flux of precipitating protons (0-detector data) is sufficiently asymmetric. The maximal flux position along MLT is moving from pre-midnight sector in quiet time to post-midnight one before and during SSC and moving back during recovery phase. The longitudinal distribution of precipitation maxima demonstrates the local increase in afternoon sector (approximately at 13:30 MLT) and decrease in the dusk one during SSC. These features are evident consequence of the magnetosphere compression. To identify the origin of the particles, the locations of maximum fluxes have been projected to the magnetosphere. It was determined that during geomagnetic storm main and recovery phases maximum fluxes were measured at latitudes poleward of the isotropic boundary. To evaluate the trapped particle flux asymmetry, the particle fluences (90-detector data) were calculated along the satellite orbit from L = 2 to the isotropic boundary. The total fluences of trapped particles calculated along the satellite orbit show regular asymmetry between dusk and dawn during main and recovery phases. The maximal intensity of proton fluxes of both investigated populations located poleward and equatorward of the isotropic boundary is achieved during SSC. The total flux measured during crossing the anisotropic region can be considered as a proxy for ring current injection rate.     

Spatial distribution of the magnetosheath ion flux

The magnetosheath plays a crucial role in solar wind-magnetosphere interaction because it is the magnetosheath magnetic field and plasma that interact with the magnetopause and magnetosphere, not the unshocked solar wind. We are presenting ion flux measurement statistics at both the dawn and dusk flanks of the magnetosheath and their comparison with a gasdynamic magnetosheath model. The study is based on three years of INTERBALL-1 measurements supported by simultaneous WIND solar wind and magnetic field observations. Statistical processing has shown (1) the limitations of the gasdynamic model, (2) the conditions favorable for the creation of a plasma depletion layer adjacent to the flank magnetopause, (3) strong dawn-dusk asymmetry of the ion fluxes, and (4) an evidence for the presence of a slow mode front adjacent to the magnetopause.     

地球磁场对太阳风扰动响应的研究

分别对行星际激波、太阳风动压增大事件和减小事件的地球磁场响应进行了比较. 分析结果表明, 同步轨道磁场对太阳风扰动在向阳面产生较强的正响应, 在背阳面 响应较弱且有时会出现负响应, 地磁指数SYM-H对太阳风扰动的响应为正响应. 同时还得出, 向阳侧同步轨道磁场响应幅度d Bz与地磁指数响应幅度d SYM-H、上下游动压均方差均具有较好的相关性. 地磁指数响应幅度与同步轨道磁场响应幅度相关关系在激波和动压增大事件中具有一致性, 动压减小事件出 现明显差异, 这说明激波和动压增大事件在影响地球磁场方面具有某种共性.      

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

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

磁暴期间电离层电磁离子回旋波的地方时分布差异

利用Swarm卫星的高精度（50 Hz）磁场观测数据,对2015年3月16—25日磁暴期间中纬度电离层电磁离子回旋（EMIC）波时空分布特征进行了研究.结果表明:晨侧EMIC波事件数与昏侧大致相当,午前时段明显多于子夜前时段.昏侧EMIC波高发生率与等离子体羽状结构有关,晨侧EMIC波高发生率与太阳风动压增强及稠密冷等离子体有关.晨侧-正午前EMIC波频率高于昏侧-子夜前,表明源区位置以及离子成分占比存在地方时差异.昏侧事件大多发生在早期恢复相,晨侧事件大多发生在晚期恢复相,晨-昏两侧的时间差异源于磁暴期间高能离子西向漂移所需时间及等离子体层顶位置的地方时差异.磁暴期间,EMIC波以H+波和He+为主,其中H+波主要分布在06:00 MLT—10:00 MLT（磁地方时）扇区,He+波主要分布在18:00 MLT—22:00 MLT扇区.在磁暴主相期间没有出现H+带波,但是出现He+-O+双波段EMIC波,表明磁暴主相期间环电流高浓度氧离子对H+带EMIC波具有抑制作用.