Modeling PSBL high speed ion beams observed by Cluster and Double Star

On October 8, 2004, the Cluster and Double Star spacecraft crossed the near-Earth (12–19 R_E) magnetotail neutral sheet during the recovery phase of a small, isolated substorm. Although they were separated in distance by ∼7 R_E and in time by ∼30 min, both Cluster and Double Star observed steady, but highly structured Earthward moving >1000 km/s high speed H⁺ beams in the PSBL. This paper utilizes a global magnetohydrodynamic (MHD) simulation driven by Wind spacecraft solar wind input to model the large-scale structure of the PSBL and large-scale kinetic (LSK) particle tracing calculations to investigate the similarities and differences in the properties of the observed beams. This study finds that the large-scale shape of the PSBL is determined by the MHD configuration. On smaller scales, the LSK calculations, in good qualitative agreement with both Cluster and Double Star observations, demonstrated that the PSBL is highly structured in both time and space, on time intervals of less than 2 min, and spatial distances of the order of 0.2–0.5 R_E. This picture of the PSBL is different from the ordered and structured region previously reported in observations.     

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.     

On high-altitude electric fields in the auroral downward current region and their coupling to ionospheric electric fields

The downward field-aligned current region plays an active role in magnetosphere–ionosphere coupling processes associated with aurora. A quasi-static electric field structure with a downward parallel electric field forms at altitudes between 800 km and 5000 km, accelerating ionospheric electrons upward, away from the auroral ionosphere. Other phenomena including energetic ion conics, electron solitary waves, low-frequency wave activity, and plasma density cavities occur in this region, which also acts as a source region for VLF saucers. Results are presented from high-altitude Cluster observations with particular emphasis on the characteristics and dynamics of quasi-static electric field structures. These, extending up to altitudes of at least 4–5 Earth radii, appear commonly as monopolar or bipolar electric fields. The former occur at sharp boundaries, such as the polar cap boundary whereas the bipolar fields occur at softer boundaries within the plasma sheet. The temporal evolution of quasi-static electric field structures, as captured by the pearls-on-a-string configuration of the Cluster spacecraft, indicates that the formation of electric field structures and of ionospheric plasma density cavities are closely coupled processes. A related feature of the downward current is a broadening of the current sheet with time, possibly related to the depletion process. Preliminary studies of the coupling of electric fields in the downward current region, show that small-scale structures are typically decoupled from the ionosphere, similar to what has been found for the upward current region. However, exceptions are also found where small-scale electric fields couple perfectly between the ionosphere and Cluster altitudes. Recent FAST results indicate that the degree of coupling differs between sheet-like and curved structures, and that it is typically partial. The electric field coupling further depends on the current–voltage relationship, which is highly non-linear in the downward current region, and still unrevealed, as to its specific form.     

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.     

2001年9月15日Cluster卫星对电流片的观测研究

2001年9月15日0430-0515 UT期间,Cluster卫星多次穿越磁尾电流片.由FGM、CIS等仪器获得了电流片磁场,粒子速度等数据变化情况,并得到了磁尾高速流的两次明显反转.本文采用GSM坐标系.利用求磁场空间梯度张量的方法获得了越尾电流,其电流密度的峰值为28nA/m2.并对电流片在垂直方向上的摆动和在晨昏方向的波动现象进行了分析.数据显示此时电流片为薄电流片并有一个变薄的趋势,其厚度大约为0.2-0.3 Re.磁场重联时地向流与尾向流均超过了1000 km/s,并测得了电流片的法线方向和运动速度,从而得到了磁尾电流片的结构和运动情况.      

中性电流片的形成和环形日冕瞬变

在双极背景场下,光球层反向磁通量的喷发将会在新老磁场之间形成中性电流片.本文从理想磁流体方程组出发,考虑磁场和日冕等离子体的相互作用,对上述电流片的形成过程进行了数值研究.结果表明,对亚音速喷发,将由里向外形成四个区域:(1)由喷发物质直接形成的低温,高密度日珥,位于最里层;(2)紧挨抛射日珥的低温稀疏区;(3)喷发物质和日冕物质向中性电流片集中形成的高温.高密度物质环;(４)在环的周围,由快磁声波形成的,密度略比日冕背景为高的前鞘区.上述结构与典型的环形日冕瞬变的观测特征相符.由此表明双极背景场下反向磁通量的喷发可能是触发这类瞬变的重要机制.      

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.     

Three Dimensional Simulation of the Steady Solar Wind in Carrington Rotation 1935

A three-dimensional MHD simulation is conducted to study the steady solar wind in Carrington Rotation (CR) 1935 by using the three-dimensional numerical magnetohydrodynamic (MHD) model introduced by Feng et al. The numerical results demonstrate that the neutral current sheet has two peaks and two valleys, which is consistent with the result of PFSS model at Wilcox Solar Observatory (WSO). The obtained proton number density at 2.5 Rs is of the same order of magnitude as the result estimated from K-coronal brightness during the CRs 1733-1742 in 1983made by Wei et al. The radial velocity profile along heliocentric distance is consistent with that of low solar wind speed deduced by Sheeley and Wang et al. However, it is not able to reproduce the fast-speed flow in coronal holes and slow solar wind in streamers because of oversimplified energy equation adopted in our model. Future efforts must be made to remedy this deficiency.      

引导场对无碰撞电流片低频电磁不稳定性的影响

通过数值求解无碰撞电流片中可压缩磁流体力学模型下得到的一般形式的色散关系,讨论了无碰撞电流片中引导场对低频电磁波不稳定性的影响.结果表明,平衡态磁场中的引导场,对于三维扰动传播的波不稳定性有很强的影响.(1)在电流片中间平面上(z=0),无引导场时,没有不稳定性发生,但若存在引导场,不稳定性便发生,并随着引导场的增强,不稳定性明显增强,不稳定的波模可能是低混杂模.(2)在中间平面附近(z=0.2),电流片是不稳定的.随着引导场的增强,不稳定性增长率明显地增强,不稳定的波模从平行和反平行两个方向传播变为反平行方向一个方向传播,并且是斜传播的,具有低频哨声模或低混杂模的特征.(3)在电流片边缘附近(z=0.8),引导场对不稳定的波模和增长率没有明显影响,不稳定的波模都是准平行的哨声波.      

Comparison of BBfs Observed by Single and Multi Satellites

This paper, using the data of three Cluster satellites, compares the observations of Bursty Bulk Flow (BBF) by single satellite with those by multi satellites. The results indicate that there exists remarkable difference between observations of BBF by single satellite and multi satellites. The observations of BBF by a single satellite depend on its position relative to the flow channel. The difference is caused by the localization characteristics of fast flows in the plasma sheet, and can lead to diverging views about substorm and causal relations among substorm phenomena.      

Transient fast reconnection in dynamic current sheets with anomalous resistivity

Two-dimensional compressible magnetohydrodynamic simulations of current sheet dynamics under the influence of multiple anomalous resistivity areas and slight asymmetries are presented. Following induced tearing and multiple coalescence, a plasmoid is formed and accelerated. Dominant X-points drive the dynamical evolution and lead to transient occurrence of a Petschek-like reconnection geometry. The dependence of current density extrema, plasmoid bulk velocity and maximum reconnection rate on the Lundquist number is examined.     

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.     

Universal scaling laws for fully-developed magnetic field turbulence near and far upstream of the Earth’s bow shock

We analyze the multifractal scaling of the modulus of the interplanetary magnetic field near and far upstream of the Earth’s bow shock, measured by Cluster and ACE, respectively, from 1 to 3 February 2002. The maximum order of the structure function is carefully estimated for each time series using two different techniques, to ensure the validity of our high-order statistics. The first technique consists of plotting the integrand of the pth order structure function, and the second technique is a quantitative method which relies on the power-law scaling of the extreme events. We compare the scaling exponents computed from the structure functions of magnetic field differences with the predictions obtained by the She–Lévêque model of intermittency in anisotropic magnetohydrodynamic turbulence. Our results show a good agreement between the model and the observations near and far upstream of the Earth’s bow shock, rendering support for the modelling of universal scaling laws based on the Kolmogorov phenomenology in the presence of sheet-like dissipative structures.     

Application of an MHD simulation to the study of substorms

The substorm mechanism is studied by the numerical solutions obtained from a resistive magnetohydrodynamic (MHD) simulation. After a southward turning of the interplanetary magnetic field (IMF), the simulation results reproduce observed features of the growth phase. The numerical solutions show that the plasma sheet thinning during the growth phase is formed under the dynamic balance between the flux pileup from the midtail and the flux removal toward the dayside controlled by the convection in the magnetosphere-ionosphere (M-I) coupling system. After the growth phase, dipolarization is generated in the near-earth tail accompanied by a plasma injection into the inner magnetosphere, the formation of plasmoid in the midtail, and the enhancement of the nightside field-aligned currents (FACs). The direct cause of this onset is the state (phase space) transition of the convection system from a thinned state to a dipolarized state associated with a self-organization in the nonlinear system.     

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.     

Plasma sheet oscillations and their relation to substorm development: Cluster and double star TC1 case study

We examined two consecutive plasma sheet oscillation and dipolarization events observed by Cluster in the magnetotail, which are associated with a pseudo-breakup and a small substorm monitored by the IMAGE spacecraft. Energy input from the solar wind and an associated enhancement of the cross-tail current lead to current sheet thinning and plasma sheet oscillations of 3–5 min periods, while the pseudo-breakups occur during the loading phase within a spatially limited area, accompanied by a localized dipolarization observed by DSP TC1 or GOES 12. That is, the so-called “growth phase” is a preferable condition for both pseudo-breakup and plasma sheet oscillations in the near-Earth magnetotail. One of the plasma sheet oscillation events occurs before the pseudo-breakup, whereas the other takes place after pseudo-breakup. Thus there is no causal relationship between the plasma sheet oscillation events and pseudo-breakup. As for the contribution to the subsequent small substorm, the onset of the small substorm took place where the preceding plasma sheet oscillations can reach the region.     

日冕冲浪形成的磁流体动力学模拟

应用二维时变可压缩磁流体动力学模拟,数值研究了双极-单极磁场中电阻撕裂模不稳定性引起的磁场重联过程,用于模拟日冕冲浪的形成.结果表明,在包含有三区——双极场、电流片和单极场的磁静力平衡初态下,双极场和单极场中的磁力线将会直接重联,磁场演变成鞭状(whip)结构.由弯曲磁力线支撑的等离子体团向上运动到最高位置后,逐渐下落和弥散.等离子体团上升速度可达到0.10v_A(v_A为双极场中的Alfv'én速度).模拟结果证实日冕冲浪的形成可能与双极-单极场中的磁场重联密切相关.      

Three-dimensional global simulation of multiple ICMEs’ interaction and propagation from the Sun to the heliosphere following the 25–28 October 2003 solar events

This study performs simulations of interplanetary coronal mass ejection (ICME) propagation in a realistic three-dimensional (3D) solar wind structure from the Sun to the Earth by using the newly developed hybrid code, HAFv.2+3DMHD. This model combines two simulation codes, Hakamada–Akasofu–Fry code version 2 (HAFv.2) and a fully 3D, time-dependent MHD simulation code. The solar wind structure is simulated out to 0.08 AU (18 Rs) from source surface maps using the HAFv.2 code. The outputs at 0.08 AU are then used to provide inputs for the lower boundary, at that location, of the 3D MHD code to calculate solar wind and its evolution to 1 AU and beyond. A dynamic disturbance, mimicking a particular flare’s energy output, is delivered to this non-uniform structure to model the evolution and interplanetary propagation of ICMEs (including their shocks). We then show the interaction between two ICMEs and the dynamic process during the overtaking of one shock by the other. The results show that both CMEs and heliosphere current sheet/plasma sheet were deformed by interacting with each other.     

Structure and kinetic properties of slow-mode shocks associated with magnetic reconnection in the near-Earth magnetotail

We study the structure and kinetic properties of slow-mode shocks near the plasma sheet boundary layer (PSBL) associated with magnetic reconnection by Cluster observation. The presence of slow-mode shocks is confirmed by traditional Rankine–Hugoniot (RH) analysis and Monte-Carlo shock fitting method. The Walén analysis, applied to the tailward flow associated with slow-mode shocks, also supports that plasma was accelerated across a Petschek-type slow-mode shock connected to the diffusion region. Back-streaming ions were observed on the shock layer, and cold ions were accelerated and heated by slow-mode shocks. In addition, whistler and electrostatic solitary waves were observed around the slow-mode shocks. These waves might be excited by the observed field-aligned electron beams near the shocks.     

Cluster observations showing the indication of the formation of a modified-two-stream instability in the geomagnetic tail

This study presents several observations of the Cluster spacecraft on September 24, 2003 around 15:10 UT, which show necessary prerequisites and consequences for the formation of the so-called modified-two-stream instability (MTSI). Theoretical studies suggest that the plasma is MTSI unstable if (1) a relative drift of electrons and ions is present, which exceeds the Alfvèn speed, and (2) this relative drift or current is in the cross-field direction. As consequences of the formation of a MTSI one expects to observe (1) a field-aligned electron beam, (2) heating of the plasma, and (3) an enhancement in the B-wave spectrum at frequencies in the range of the lower-hybrid-frequency (LHF). In this study we use prime parameter data of the CIS and PEACE instruments onboard the Cluster spacecraft to verify the drift velocities of ions and electrons, FGM data to calculate the expected LHF and Alfvèn velocity, and the direction of the current. The B-wave spectrum is recorded by the STAFF instrument of Cluster. Finally, a field aligned beam of electrons is observed by 3D measurements of the IES instrument of the RAPID unit. Observations are verified using a theoretical model showing the build-up of a MTSI under the given circumstances.