Study of three dimensional structure of the fast convection flow in the plasma sheet by MHD simulations on the basis of spontaneous fast reconnection model

Three dimensional structure of the fast convection flow in the plasma sheet is examined using magnetohydrodynamic (MHD) simulations on the basis of spontaneous fast reconnection model. The fast flow observed in the near-Earth magnetotail is one of the key phenomena in order to understand the causal relationship between magnetic substorm and magnetic reconnection. In this paper, we focus on this earthward fast flow in the near-Earth magnetotail. Our previous studies have shown that the fast reconnection produces the Alfvénic fast reconnection outflow and drastic magnetic field dipolarization in the finite extent. In this paper, the results of our simulations are compared with those of the in-situ observations in the geomagnetotail. They have consistent temporal profiles of the plasma quantities. It is suggested that the fast convection flows are caused by spontaneous fast reconnection.     

Particle simulation study of substrom triggering with a southward IMF

This paper reports the spatial and temporal development of bursty bulk flows (BBFs) created by reconnection as well as current disruptions (CDs) in the near-Earth tail using our 3-D global electromagnetic (EM) particle simulation with a southward turning interplanetary magnetic field (IMF) in the context of the substorm onset. Recently, observations show that BBFs are often accompanied by current disruptions for triggering substorms. We have examined the dynamics of BBFs and CDs in order to understand the timing and triggering mechanism of substorms. As the solar wind with the southward IMF advances over the Earth, the near-Earth tail thins and the sheet current intensifies. Before the peak of the current density becomes maximum, reconnection takes place, which ejects particles from the reconnection region. Because of earthward flows the peak of the current density moves toward Earth. The characteristics of the earthward flows depend on the ions and electrons. Electrons flow back into the inflow region (the center of reconnection region), which provides current closure. Therefore the structure of electron flows near the reconnection region is rather complicated. In contrast, the ion earthward flows are generated far from the reconnection region. These earthward flows pile up near the Earth. The ions mainly drift toward the duskside. The electrons are diverted toward the dawnside. Due to the pile-up, dawnward current is generated near Earth. This dawnward current dissipates rapidly with the sheet current because of the opposite current direction, which coincides with the dipolarization in the near-Earth tail. At this time the wedge current may be created in our simulation model. This simulation study shows the sequence of the substorm dynamics in the near-Earth tail, which is similar to the features obtained by multisatellite observations. Identification of the timing and mechanism of triggering substorm onset requires further studies in conjunction with observations.     

Theoretical properties of offset bipolar electric field solitary structures in space plasmas

The bipolar electric field solitary (EFS) structures have been frequently observed in the near Earth plasma regions, such as auroral zone, magnetopause, cusp regions, and magneto-tail. Sometimes these structures are observed as offset bipolar structures. In this paper, the properties of the offset bipolar EFS structures parallel to the magnetic field are studied with an ion fluid model in a cylindrical symmetry by considering electrostatic condition. The model results show that the offset bipolar EFS structures can develop from both ion-acoustic waves and ion cyclotron waves, and propagate along the magnetic field line in the space plasmas if plasma satisfies some conditions. The offset bipolar EFS structures can have both polarities. It will be first negative pulse and then positive pulse if the initial electric field E₀ < 0 or reverse in order if E₀ > 0. The amplitude of the offset bipolar EFS structures first decreases and then increases with the wave propagation velocity. Some results from our model are consistent with the observations.     

Cluster observations of earthward propagating plasmoid and flux ropes in the near-tail during the course of a substorm event

We analyze the magnetic structures in the near-tail at Xgsm = −17.5 Re on September 19, 2003 by Cluster. During the course of a substorm event, the earthward propagating plasmoid and flux ropes in the near-tail are observed. The earthward propagating plasmoid is associated with the bipolar Bz and By signatures. The two flux ropes are embedded within the earthward plasma flows, which might be referred to the population as ‘‘BBF-type’’ flux ropes. The first flux rope diameter is about 0.7 Re and duration based upon the Bz signature is ∼20 s, while the second one diameter is about 1.4 Re and duration is ∼30 s. The earthward propagating plasmoid and flux ropes could have influence upon the dipolarization and injection in inner magnetosphere. The Cluster observations of earthward propagating plasmoid and flux ropes can be interpreted as strong evidence for multiple X-lines. Our observations are consistent with that multiple plasmoids or flux ropes are formed repeatedly and ejected tailward in the course of geomagnetically active time.     

有限β等离子体中密度和磁场不均匀驱动的动理学Alfven波

在分析有限β等离子体中的密度、磁场不均匀引起的漂移波不稳定性的基础上，剖析了漂移波不稳定性对动理学Alfven波激发的作用．动理学理论能正确地处理有限拉莫半径效应和波粒共振相互作用，本文根据带电粒子在电磁场中的运动特性，采用Vlasov方程描述离子运动，运用漂移动理学方程对电子运动进行描述．密度不均匀和磁场不均匀对产生漂移不稳定性的对比分析表明：在有限β等离子体中，密度不均匀比磁场不均匀更易激发漂移不稳定性，且密度不均匀激发漂移不稳定性中的能量转换和转移更为强烈．这种能量的转换为动理学Alfven波的激发提供了物理基础．所得数值解表明：动理学Alfven波在磁层中能广泛地被激发产生，特别是在磁层空间的极尖区、磁层顶和等离子体片边界层等具有明显的不均匀性区域中更容易被激发产生．本文的研究结果进一步表明动理学Alfven波对磁层空间中能量传输具有重要作用．     

Standing Alfvén waves with m ≫ 1 in a dipole magnetosphere with moving plasma and aurorae

The structure of standing Alfvén waves with large azimuthal wave numbers (m ? 1) is studied in a dipole model of the magnetosphere with rotating plasma. In the direction across magnetic shells the structure of such waves is determined by their dispersion associated with curvature of geomagnetic field lines and corresponds to the travelling wave localized between toroidal and poloidal resonant surfaces. In projection into the ionosphere (along geomagnetic field lines) this structure is similar to the structure of a discrete auroral arc. The azimuthal structure of an auroral arc is similar to azimuthal structure of Alfvén waves with m ∼ 100. Possible interaction mechanisms between the Alfvén waves and energetic electron fluxes forming auroral arcs are discussed.     

Multi-point observations of earthward fast flow in the plasma sheet by virtual satellites located in the MHD simulation domain

Time profiles of some physical values in earthward fast flows in the plasma sheet are observed at three dimensionally different positions by employing virtual satellites located in the three-dimensional magnetohydrodynamic simulation domain, and these simulations are done on the basis of the spontaneous fast reconnection model. In the spontaneous fast reconnection evolution, the width of the flow channel is narrow in the dawn-dusk direction, and it does not spread until the plasma collides with the magnetic loop. The enhancements in B_z and V_x are larger at the center of the fast flow channel than those at its dawn and dusk edges, reflecting the differences in the reconnection rate in the diffusion region. The enhancement in V_x is shorter near the plasma sheet boundary layer than that near the neutral sheet, reflecting the changes in the thickness of the flow channel.     

Alfvén波在任意流场、磁场位形中的传播

本文研究具有任意偏振和非单色的Alfvén波, 在任意流场、磁场位形且有粘、热各向异性等离子体流介质中的传播.采用WKB近似可以获得它的波幅矢量的张量表达式, 它是介质密度、速度、磁场、Alfvén波速、热各向异性和粘性的函数.文中对简单情形进行了讨论, 结果表明:本文结果是对普遍情形的一种简明的综合和概括, 便于研究三维传播问题.      

Acceleration of charged particles by fluctuating and steady electric fields in a X-type magnetic field

Release of stored magnetic energy via particle acceleration is a characteristic feature of astrophysical plasmas. Magnetic reconnection is one of the mechanisms for releasing energy from magnetized plasmas. Collisionless magnetic reconnection could provide both the energy release mechanism and the particle accelerator in space plasmas. Here we studied particle acceleration when fluctuating (in-time) electric fields are superposed on an static X-type magnetic field in collisionless hot solar plasma. This system is chosen to mimic the reconnective dissipation of a linear MHD disturbance. Our results are compared to particle acceleration from constant electric field superposed on an X-type magnetic field. The constant electric field configuration represents the effects of steady state magnetic reconnection. Time evolution of ion and electron distributions are obtained by numerically integrating particle trajectories. The frequencies of the electric field represent a turbulent range of waves. Depending on the frequency and amplitude of the electric field, electrons and ions are accelerated to different degrees and have energy distributions of bimodal form consisting of a lower energy part and a high energy tail. For frequencies (ω in dimensioless units) in the range 0.5 ? ω ? 1.0 a substantial fraction (20%–30%) of the proton distribution is accelerated to gamma-ray producing energies. For frequencies in the range 1 ? ω ? 100.0 the bulk of the electron distribution is accelerated to hard X-ray producing energies. The acceleration mechanism is important for solar flares and solar noise storms but it could be applicable to all collisionless astrophysical plasmas.     

Electron acceleration and parallel electric fields due to kinetic Alfvén waves in plasma with similar thermal and Alfvén speeds

We investigate electron acceleration due to shear Alfvén waves in a collissionless plasma for plasma parameters typical of 4–5R_E radial distance from the Earth along auroral field lines. Recent observational work has motivated this study, which explores the plasma regime where the thermal velocity of the electrons is similar to the Alfvén speed of the plasma, encouraging Landau resonance for electrons in the wave fields. We use a self-consistent kinetic simulation model to follow the evolution of the electrons as they interact with a short-duration wave pulse, which allows us to determine the parallel electric field of the shear Alfvén wave due to both electron inertia and electron pressure effects. The simulation demonstrates that electrons can be accelerated to keV energies in a modest amplitude sub-second period wave. We compare the parallel electric field obtained from the simulation with those provided by fluid approximations.     

磁扰动和磁静时近地等离子体片中脉冲电场的分布

Ｇｅｏｔａｉｌ卫星的电场数据被用于分析近地磁尾等离子体片中电场在磁扰动（Ｄｓｔ＜－２５ｎＴ）和磁静时（Ｄｓｔ＞－２５ ｎＴ的统计分布．结果表明，伴随着地向高速离子流，在Ｘ＞－１６Ｒｅ以内区域出现强电场（高达 ５—８ ｍＶ／ｍ）．磁扰动期间强电场的幅值较磁静时大，并且出现在更靠近地球的位置．较强和较靠近地球的强电场与磁扰动时更薄的等离子体片和更接近地球的等离子体片内边界相联系．观测结果意味着磁扰动期间的亚暴可能更有效地将高能粒子注射到环电流中．这对磁暴和亚暴的关系问题的解决有重要意义．     

Generation mechanism of the whistler-mode waves in the plasma sheet prior to magnetic reconnection

The whistler-mode waves and electron temperature anisotropy play a key role prior to and during magnetic reconnection. On August 21, 2002, the Cluster spacecrafts encountered a quasi-collisionless magnetic reconnection event when they crossed the plasma sheet. Prior to the southward turning of magnetospheric magnetic field and high speed ion flow, the whistler-mode waves and positive electron temperature anisotropy are simultaneously observed. Theoretic analysis shows that the electrons with positive temperature anisotropy can excite the whistler-mode waves via cyclotron resonances. Using the data of particles and magnetic field, we estimated the whistler-mode wave growth rate and the ratio of whistler-mode growth rate to wave frequency. They are 0.0016f_ce (Electron cyclotron frequency) and 0.0086f_ce, respectively. Therefore the whistler-mode waves can grow quickly in the current sheet. The combined observations of energetic electron beams and waves show that after the southward turning of magnetic field, energetic electrons in the reconnection process are accelerated by the whistler-mode waves.     

Electromagnetic ion cyclotron resonance heating in the VASIMR

Plasma physics has found an increasing range of practical industrial applications, including the development of electric spacecraft propulsion systems. One of these systems, the Variable Specific Impulse Magnetoplasma Rocket (VASIMR) engine, both applies several important physical processes occurring in the magnetosphere. These processes include the mechanisms involved in the ion acceleration and heating that occur in the Birkeland currents of an auroral arc system. Auroral current region processes that are simulated in VASIMR include lower hybrid heating, parallel electric field acceleration and ion cyclotron acceleration. This paper will focus on using a physics demonstration model VASIMR to study ion cyclotron resonance heating (ICRH). The major purpose is to provide a VASIMR status report to the COSPAR community. The VASIMR uses a helicon antenna with up to 20 kW of power to generate plasma. This plasma is energized by an RF booster stage that uses left hand polarized slow mode waves launched from the high field side of the ion cyclotron resonance. The present setup for the booster uses 2–4 MHz waves with up to 20 kW of power. This process is similar to the ion cyclotron heating in tokamaks, but in the VASIMR the ions only pass through the resonance region once. The rapid absorption of ion cyclotron waves has been predicted in recent theoretical studies. These theoretical predictions have been supported with several independent measurements in this paper. The ICRH produced a substantial increase in ion velocity. Pitch angle distribution studies show that this increase takes place in the resonance region where the ion cyclotron frequency is equal to the frequency on the injected RF waves. Downstream of the resonance region the perpendicular velocity boost should be converted to axial flow velocity through the conservation of the first adiabatic invariant as the magnetic field decreases in the exhaust region of the VASIMR. In deuterium plasma, 80% efficient absorption of 20 kW of ICRH input power has been achieved. No evidence for power limiting instabilities in the exhaust beam has been observed.     

Study of waves in the magnetotail region with cluster and DSP

The study of the neutral sheet is of fundamental importance in understanding the dynamics of the Earth’s magnetosphere. From the earliest observation of the magnetotail, it has been found that the neutral sheet frequently appears to be in motion due to changing solar wind conditions and geomagnetic activity. Multiple crossings of the neutral sheet by spacecraft have been attributed to a flapping motion of the neutral sheet in the north–south direction, a wavy profile either along the magnetotail or the dawn–dusk direction. Cluster observations have revealed that the flapping motions of the Earth’s magnetotail are of internal origin and that kink-like waves are emitted from the central part of the tail and propagate toward the tail flanks. This flapping motion is shown here to propagate at an angle of ∼45° with x_GSM. A possible assumption that the flapping could be created by a wake travelling away from a fast flow in the current sheet is rejected. Other waves in the magnetotail are found in the ULF range. One conjunction event between Cluster and DoubleStar TC1 is presented where all spacecraft show ULF wave activity at a period of approximately 5 min during fast Earthward flow. These waves are shown to be Kelvin–Helmholtz waves on the boundaries of the flow channel. Calculations show that the conversion of flow energy into magnetic energy through the Kelvin–Helmholtz instability can contribute to a significant part of flow breaking between Cluster and DoubleStar TC1.     

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.     

Substorms in the inner plasma sheet

Thin Current Sheets (TCS) are regularly formed prior to substorm breakup, even in the near-Earth plasma sheet, as close as the geostationary orbit. A self-consistent kinetic theory describing the response of the plasma sheet to an electromagnetic perturbation is given. This perturbation corresponds to an external forcing, for instance caused by the solar wind (not an internal instability). The equilibrium of the configuration of this TCS in the presence of a time varying perturbation is shown to produce a strong parallel thermal anisotropy (T T) of energetic electrons and ions (E>50keV) as well as an enhanced diamagnetic current carried by low energy ions (E<50keV). Both currents tend to enhance the confinement of this current sheet near the magnetic equator. These results are compared with data gathered by GEOS-2 at the geostationary orbit, where the magnetic signatures of TCS, and parallel anisotropics are regularly observed prior to breakup. By ensuring quasi-neutrality everywhere we find, when low frequency electromagnetic perturbations are applied, that although the magnetic field line remains an equipotential to the lowest order in T_e/T_i, a field-aligned potential drop exists to the next order in (T_e/T_i). Thus the development of a TCS implies the formation of a field-aligned potential drop ( few hundred volts) to ensure the quasi-neutrality everywhere. For an earthward directed pressure gradient, a field-aligned electric field, directed towards the ionosphere, is obtained, on the western edge of the perturbation (i.e. western edge of the current sheet). Thus field aligned beams of electrons are expected to flow towards the equatorial region on the western edge of the current sheet. We study the stability of these electron beams and show that they are unstable to “High Frequency” (HF) waves. These “HF” waves are regularly observed at frequencies of the order of the proton gyrofrequency (f_H+) just before, or at breakup. The amplitude of these HF waves is so large that they can produce a strong pitch-angle diffusion of energetic ions and a spatial diffusion that leads to a reduction of the diamagnetic current. The signature of a fast ion diffusion is indeed regularly observed during the early breakup; it coincides with the sudden development of large amplitude transient fluctuations, ballooning modes, observed at much lower frequencies (fH+). These results suggest that the HF waves, generated by field-aligned electron beams, provide the dissipation which is necessary to destabilize low frequency (ballooning) modes.     

Heating heavy ions in the polar corona by collisionless shocks: A one-dimensional simulation

Recently a new model for explaining the observations of preferential heating of heavy ions in the polar solar corona was proposed ( and ). In that model the ion energization mechanism is the ion reflection off supercritical quasi-perpendicular collisionless shocks in the corona and the subsequent acceleration by the motional electric field E = −V × B/c. The mechanism of heavy ion reflection is based on ion gyration in the magnetic overshoot of the shock. The acceleration due to the motional electric field is perpendicular to the magnetic field, giving rise to large temperature anisotropy with T_⊥ ? T_∥, in agreement with SoHO observations. Such a model is tested here by means of a one dimensional test particle simulation where ions are launched toward electric and magnetic profiles representing the shock transition. We study the dynamics of O⁵⁺, as representative of coronal heavy ions for Alfvénic Mach numbers of 2–4, as appropriate to solar corona. It is found that O⁵⁺ ions are easily reflected and gain more than mass proportional energy with respect to protons.     

太阳风中的动力论Alfven激波

本文指出因动力论Alfven波的存在非线性作用，从而能引起波连续陡化;这种情况下考虑频散效应，理论上应该形成中间激波.针对这一物理思想，在二元流体模型下，引进离子声反常阻尼，数值解出中间激波的结构.     

Alfven波在微电离大气中的衰减特性研究

利用简单的偶极子地磁场模型以及大气电子密度和电导率模式, 分析地面产生的磁扰动以Alfven波的模式传播到近地空间区域. 这种地面的磁扰动可能干扰近地空间卫星对空间磁扰动的观测. 通过对地面磁扰动Alfven波模式1000 km高度内的衰减情况进行模拟, 认为在近地空间采用地磁偶极子模型是合理的. 由于衰减随扰动频率的增大而急剧增强, 分析还得到了近地卫星能够探测到地面磁扰动的最大频率. 计算结果表明, Alfven波的衰减主要集中在高度50 km以下, 这个区域内的大气电导率极其微弱, 使Alfven波的传播受到极大衰减. 0.4 Hz以下的Alfven波沿磁力线传播到1000 km高度后衰减结为原来扰动幅度的千分之一, 因此频率在0.4 Hz以下的Alfven波可能会干扰低轨卫星探测磁场脉动.      

Height and time variability of planetary wave activity

The height–season and year-to-year regularities of parameters of first and second spatial harmonics determine the structure of the stratosphere and mesosphere circulation and its variability. In the period 1992–2002 at heights 0–55 km, the amplitudes and phases of the first and second spatial harmonics in the field of temperature, geopotential height, zonal and meridional wind were calculated by the method of harmonic decomposition. Dispersion (standard or mean square deviation) of their day-to-day and year-to-year variations was calculated by their wavelength constants. Height and season patterns of variability have been estimated. The difference in height–longitude variability for wave numbers m = 1 and 2 has been discovered. At the same time, the intensity of wave disturbances for m = 1 is less than for m = 2 excluding the polar areas, where a significant variability appears at the heights 0–55 km. There is also a tendency for the intensity of year-to-year variations to decrease in comparison with day-to-day variations. In cold and warm periods the amplitude of perturbation waves with m = 2 both for day-to-day and year-to-year variations is greater than for waves with m = 1. Transient height areas in the interval of 20–30 km are more distinct for day-to-day variations of polar area.