外辐射带动态演化的粒子模拟研究

运用计算机实验方法即三维电磁粒子模拟方法初步研究了暴时扰动情况下外辐射带粒子环境的动态演化特性。模拟计算了暴时辐射带粒子环境的演化情况。模拟计算结果显示，高能质子、电子注入午辐射滞后，两者一方面沿磁力线做弹跳运动，向高纬扩展，部分注入质子和电子沉降于南北两极区域；另一方面，注入质子和电子还经历顺时针和逆时针方向的漂移运动，粒子能量越大，漂移速度越快。暴时多次注入引起整个辐射带粒子能量的大幅增强；粒子注入颗次和强度越大，辐射带粒子通量增幅越大。本项研究为开发研制完整的辐射带动态模式积累了有益的经验。     

FEM code simulation of the magnetospheric proton fluxes

We report a study of the numeric solution to the diffusive transport equation for energetic protons magnetically trapped in the Earth's equatorial magnetosphere. The analysis takes into account the pertinent physical processes in this region, including deceleration of protons by Coulomb collisional interactions with free and bound electrons, the charge exchange process, cosmic ray albedo neutron decay source, and electric and magnetic radial diffusion. These results were obtained using the Finite Element Method with magnetic moment and geomagnetic L-shell as free variables. Steady state boundary conditions were imposed at L=1 as zero distribution function and at L=7 with proton distribution function extracted from ATS 6 satellite observations. The FEM-code yields unidirectional proton flux in the energy range of 0.1–1000 MeV at the equatorial top of the geomagnetic lines, and the results are found to be in satisfactorily agreement with the empirical NASA AP-8 model proton flux within the energy range of 0.5–100 MeV. Below 500 keV, the empirical AP-8 model proton fluxes are several orders of magnitude greater than those computed with the FEM-code at L<3. This discrepancy is difficult to explain by uncertainties of boundary spectrum parameters or transport coefficients.     

Solar cosmic rays in the near Earth space and the atmosphere

Solar energetic particles (SEPs) constitute a distinct population of energetic charged particles, which can be often observed in the near Earth space. SEP penetration into the Earth’s magnetosphere is a complicated process depending on particle magnetic rigidity and geomagnetic field structure. Particles in the several MeV energy range can only access to periphery of the magnetosphere and the polar cap regions, while the GeV particles can arrive at equatorial latitudes. Solar protons with energies higher than 100 MeV may be observed in the atmosphere above ∼30 km, and those with energies more than 1 GeV may be recorded even at the sea level. There are some observational evidences of SEP influence on atmospheric processes. Intruding into the atmosphere, SEPs affect middle atmosphere odd-nitrogen and ozone chemistry. Since spatial and temporal variations of SEP fluxes in the near Earth space are controlled by solar activity, SEPs may present an important link between solar activity and climate. The paper outlines dynamics of SEP fluxes in the near Earth space during the last decades. This can be useful for tracing relationship between SEPs and atmospheric processes.     

Variation of energetic particle precipitation with local magnetic time

The detailed study of the precipitation of magnetospheric particles into the atmosphere is complicated by the rather complex spatial configuration of the precipitation region and its variability with geomagnetic activity. In this paper we will introduce polar oval coordinates and apply them to POES observations of 30 keV to 2.5 MeV electrons and comparable protons to illustrate the dependence of particle precipitation on local time and geomagnetic activity. These coordinates also allow an easy separation of the spatial precipitation patterns of solar and magnetospheric particles. The results indicate that (a) the spatial precipitation pattern of energetic magnetospheric electrons basically follows the pattern of the field parallel Birkeland currents up to MeV energies and (b) at least in the mesosphere the influence of magnetospheric electrons is comparable to the one of solar electrons. Implications for modeling of atmospheric chemistry will be sketched.     

Comparison of Antarctic riometer radio wave absorption and THEMIS mission energetic electron fluxes

Simultaneous observations of in situ plasma properties in the tail of the Earth’s magnetosphere and of ground based instruments, lying on the same geomagnetic field lines, have recently proved to yield significant new results. In most cases magnetosphere ionosphere interactions during the night-time northern hemisphere conditions are studied. Here, observations of energetic electrons in the tail of the Earth’s magnetosphere made by the THEMIS mission satellites are compared with auroral radio wave absorption determined by riometers in the Antarctic for sunlit conditions. Days for which satellites and riometers are connected by the same geomagnetic field line are selected using a geomagnetic field model. The six days analysed show clear associations between fluxes and absorptions in some cases. However, these do not necessarily correspond to conjugacy intervals. Hours of positive associations are 1.65 times those for negative associations, all hours and days considered (1.42–3.6 on five days and 0.58 on the other day). These computations are assumed appropriate since the footprints of the satellites used approximately follow corrected geomagnetic parallels for all six days studied. The use of a finer parameterization of geomagnetic models to determine conjugacy may be needed.     

用数值模拟的方法研究磁暴主相期间影响环电流分布特征的要素

环电流是距离地心2～7 Re的带电粒子围绕地球西向漂移形成的.环电流的增强将引起全球磁场的降低,反映了地磁暴的强度.磁暴主相期间,对流电场驱动等离子体片中的能量粒子经历E×B漂移与俘获注入环电流,进入损失锥的粒子沉降到大气中.本文采用磁暴主相期间环电流离子分布的模型,结合上述因素研究不同离子能量下对流电场对环电流离子通量的直接影响,以及强弱对流电场下环电流能量离子投掷角的变化,并从物理上阐述造成此种通量分布特性的原因.      

A new inversion method of plasma density distribution of plasmasphere in the geomagnetic equatorial plane from IMAGE data

The plasma density distribution of plasmasphere in the geomagnetic equatorial plane can help us study the magnetosphere like plasmasphere, ionosphere and their kinetics. In this paper, we introduce a new inversion method, GE-ART, to calculate the plasma density distribution in the geomagnetic equatorial plane from the Extreme Ultraviolet (EUV) data of IMAGE satellite under the assumption that the plasma density is constant along each geomagnetic field line. The new GE-ART algorithm was derived from the traditional Algebraic Reconstruction Techniques (ART) in Computed Tomography (CT) which was different from the several existing methods. In this new method, each value of the EUV image data was back-projected evenly to the geomagnetic field lines intersected by this EUV sight. A 3-D inversion matrix was produced by the contributions of all the voxels contained in the plasmasphere covered by the EUV sensor. That is, we considered that each value of the EUV image data was relative to the plasma densities of all the voxels passed through by the corresponding EUV radiation, which is the biggest difference to all the existing inversion methods. Finally, the GE-ART algorithm was evaluated by the real EUV data from the IMAGE satellite.     

质子地磁垂直截止刚度在2004年11月7-8日两个时刻处的数值模拟

利用单粒子轨道理论, 在T04和IGRF2000磁场模式建立的磁层模型基础上, 应用四阶龙格库塔方法, 模拟计算宇宙线带电质子在地球磁层中的运动以及沿天顶方向入射到达地球磁场内的某一特定位置, 得到了距地球表面450km高度处全球质子垂直截止刚度在2004年11月7-8日中两个时刻的计算值. 根据计算得到的684个不同位置处的截止刚度值, 分析了同一时刻地磁垂直截止刚度随磁纬和磁经的变化. 与此同时, 模拟计算了相同时刻下磁纬为30°, 磁经为0°, 45°, 90°, 135°, 180°, 225°, 270°, 315°处, 地磁有效垂直截止刚度随高度的变化情况. 结果表明, 在磁纬30°处, 其质子有效垂直截止刚度随距离地心高度的变化沿向阳处向磁尾处方向, 变化越来越缓慢.      

磁扰期间的三维漂移壳分离

基于磁层粒子动力学理论,首先对比了计算漂移壳分离的引导中心法和磁力线追踪法,计算表明两种方法的计算结果一致.然后分别采用T89c和T96磁层磁场模式,用磁力线追踪法数值计算了不同初始位置(≤9Re)、不同初始投掷角、不同Kp指数和不同太阳风压力下,带电粒子的漂移壳分离.计算结果揭示了漂移壳分离随初始位置、投掷角、Kp指数和太阳风压力的变化.其具体特征如下. (1)随着径向距离的增大,漂移壳分离效应愈加显著,由正午出发的粒子将被稳定捕获,而午夜出发的径向距离≥7Re的部分大投掷角粒子将沿磁层顶逃逸. (2)正午出发的粒子,漂移到午夜时其漂移壳随投掷角减小向外排列;午夜出发的粒子,漂移到正午时其漂移壳随投掷角增大排列; 90°投掷角粒子在磁赤道面的漂移壳沿着磁场等值线排列. (3)漂移壳分离随Kp指数和太阳风压力增大变得显著,且随这两种扰动参数的变化特征和趋势是基本相似的.      

太阳风里高能电子自发辐射的Langmuir波

本文讨论太阳风里太阳耀班高能电子产生的Langmuir波的自发辐射。理论估计得出,在时间尺度γ_km-1内,Langmuir波自发辐射的电场幅值约为10-3-10-2mV/m(依赖于高能电子速度分布的具体形式),这里γ_km为在波数k处的峰增长率。此理论结果比飞船在太阳风里的观测值低2-3个数量级。因此认为,太阳风里自发辐射产生的Langmuir波辐射是可以忽略的。      

Do ring current ions set the outer plasmasphere in motion?

GEOS-2 experiments have observed on a regular basis the transition from teh plasmasphere to the trough revealing a more tenuous and energetic population. This transition occured typically between 18 and 21 h local time under quiet geomagnetic conditions. On the plasmasphere side of the transition, large electric fields of several mV/m, mainly in the duskward direction, were observed. The present statistical results, over one year of data, assess clearly the mean magnitude and direction of these electric fields in relation to the boundary between these two plasma regimes under quiet geomagnetic conditions. The possibility of ring current energetic ions loosing part of their energy to the plasmasphere and accelerating this plasma sunward will be discussed.     

Structure and properties of the Earth's plasmasphere

The widely used concept of the plasmapause as the last closed electric field equipotential in the equatorial plane of the magnetosphere is oversimplified. The field aligned plasma motions are of substantial importance in the plasmapause formation and should be taken into account. Distributions of the main plasma parameters measured from the Prognoz-5 satellite are presented. The diurnal variations of the plasmapause height and the plasmasphere thermal properties are discussed.     

Energetic neutral atom (ENA) and charged particle periodicities in Saturn’s magnetosphere

The Magnetospheric Imaging Instrument (MIMI) on the Cassini spacecraft has observed energetic neutral atoms (ENA) and charged particles at Saturn from mid-2004 to the present. The particles often but not always reveal striking periodic behavior that seems to depend on the type of particle and spacecraft location. When subjected to a Lomb periodogram analysis, energetic electrons (>150 keV) exhibited strong frequency peaks near 10.80 h (the nominal or “base” period of Saturn kilometric radiation) during 2006–2008, but essentially no periodicity during 2005. The electron periodograms also show pronounced “double” frequency peaks in 2007 and 2008. Energetic protons (3–26 keV) show strong peaks near the same period for 2005–2007, but none for 2008. Oxygen ions at the same energies display strong peaks for 2005 and 2006, but not for 2007 and 2008. By projecting the ENA images onto Saturn’s equatorial plane or onto a plane perpendicular to the equatorial plane and then summing the data in the appropriate dimension, “strip” images can be constructed from which a time history can be derived. These time histories of ENA emissions are also subjected to a Lomb periodogram analyses. The energetic hydrogen neutrals (20–50 keV) exhibited periodic behavior only during 2007, while energetic oxygen neutrals (64–144 keV) displayed a strong SKR-like period in 2005 and 2006 but not for 2007 or 2008. Some of this behavior may be due to changing spacecraft aspect relative to the ENA emissions, and some of it may be real. This periodic behavior may be consistent with a rotating anomaly that “flashes” brightly in the midnight-to-dawn sector once per 10.8 h, with the flash parameters depending on particle species and energy.     

Towards multi-dimensional space weather modeling for energetic oxygen ions in the earth''s inner magnetosphere: Equilibrium configuration

We report the first 3+1 dimensional model development for energetic atomic oxygen ions in the Earth's radiation belts. Energetic Oxygen ions cans be supplied to the Earth's Inner magnetosphere from the sun (as a component of solar wind and solar energetic particles), from anomalous cosmic rays, and from acceleration processes acting on ionospheric atomic oxygen ions. We have built a multi-dimensional oxygen ion model in the following free parameters: geomagnetic L-shell, the magnetic moment, the second adiabatic invariant, and the discrete charge state number. Quiet time, steady state oxygen ion distributions have been obtained numerically from an assumed outer radiation zone boundary condition at L=7, average values of the radial diffusion coefficients, and standard values for the exospheric neutral densities due to the MSIS-86 upper atmosphere and exosphere neutral thermal particle density model. Average distributions of free electrons in the plasmasphere were also assumed with a mean plasmapause location just beyond L=4. We included the six lowest ionic charge states of atomic oxygen (¹⁶O) based on an existing charge exchange cross section compilation by Spjeldvik and Fritz (1978). Computed oxygen ion distributions include the resulting equilibrium structure of energy oxygen ions between 10 KeV and 100 MeV.     

Plasma convection model and generation of longitudinal currents in the earth's magnetosphere

The aim of this paper is to investigate processes in the magnetosphere and in particular the problems of the interaction of the solar wind with the Earth's magnetic field to produce large-scale convection, electric fields and longitudinal currents in the magnetosphere. The investigation is carried out in the frame of magnetic hydrodynamics. The reason for such an approach can be found in /1/. When calculating the transfer coefficients, the Böhm approximation is used, i.e. it is considered that the plasma in the near-equatorial part of the magnetosphere (quasiplanar geometry is used in the problem for simplification) is sufficiently turbulent that the condition ωτ ≈ 1 is valid (ω is the Larmor frequency of electrons, τ is effective time between two Quasi-collisions). The main subjects of investigation in this paper are the input near the equatorial boundary layer and the plasma layer of the magnetosphere tail.     

High energy particle dispersion events observed by interball-1 and -2

Time period from October 1996 until January 1998 was checked on high energy resolution DOK2 energetic particle instrument measurements on Interball-1 and Interball-2 for the ion (> 20 keV) dispersive events (EDIS) with the exclusion of Interball-1 orbit parts in the tail. A variety of energy dispersive events, both in ion and electron spectra with different duration is found in the auroral regions, in the outer magnetosphere and near the cusp. While EDIS were observed in all sectors of MLT, the best conditions for their observation were in the afternoon local time. The characteristics of dispersive events observed by DOK2 are consistent with their explanation by the gradient-curvature drift of particles from the injection point(s) in the night local time sector given in Lutsenko at al., 2000a, b.     

On the transport of cosmic rays in the distant heliosheath

This paper discusses the transport of energetic charged particles through a sectored magnetic field in distant regions of the inner heliosheath. As the plasma flow slows down on approach to the stagnation point on the heliopause, the distance between the folds of the current sheet decreases to the point where it becomes comparable to the cyclotron radius of a cosmic ray particle. Under these conditions a particle can effectively drift across the stack of magnetic sectors with a speed comparable with the particle’s velocity. For a random distribution of current sheet separation distances, a diffusive transport across the stack of sectors occurs instead. The proposed mechanism could have contributed to unusually high intensities of galactic cosmic rays measured by Voyager 1 in the heliosheath during 2009–2010.     

1998年9月24日地磁暴前宇宙线各向异性特性分析

分析了日本Nagoya 宇宙线闪烁体望远镜30°, 49°, 64° 倾角的东、西、南、北方向探测数据的变化特点, 运用小波分析方法定性地探讨了磁暴前后宇宙线南北、东西各向异性的变化特征. 研究发现, 当发 生大地磁暴时, 地面宇宙线强度的各向异性特征将发生非常大的变化, 这种变化一般在磁暴发生前10~20 h就开始出现. 当描述这种各向异性特征的各向异性指数的小波系数变化达到一定阈值时, 就可能有大地磁暴发生.      

极光西向涌浪的二维模式(线性部分)

极光西向涌浪(WTS)被定义为磁层亚暴开始的重要标志, 它是伴随带电粒子沉降过程在极光带电离层出现的特定的极光运动形态.本文给出了极光西向涌浪的二维动力学模式, 亚暴开始, 伴随粒子沉降在电离层产生离子密度梯度(电导率梯度), 在背景电场中激发出低频漂移波, 它的传播给出极光西向涌浪运动的主要特性.