Inconsistency of magnetic field and plasma velocity variations in the distant plasma sheet: Violation of the “frozen-in” criterion?

Measurements of the magnetic field and low energy plasma by the GEOTAIL spacecraft have been used to study the relationship between variations of the plasma velocity and of the magnetic field in the distant (100–200 R_E) and middle (40–80 R_E) tail. The analysis was carried out separately for the tail lobes and the plasma sheet. It is shown that the absolute values of the magnetic field and plasma velocity, as well as their corresponding components (V_X and B_X, V_Y and B_Y, V_Z and B_Z), are linearly connected in the tail lobes. In the plasma sheet, however, the plasma velocity and the magnetic field do not seem to be related to one another. The distant plasma sheet seems to be in a regime of turbulence. The diffusion coefficients estimated from our data set of the velocity parameters in the plasma sheet are in good agreement with the theoretical predictions of Antonova and Ovchinnikov (1996, 1999).     

IMF effect on ionospheric trough occurrence at equinoxes

Previous observations have shown that there is a relationship between the F region trough and both B_z and B_y components of the interplanetary magnetic field (IMF). Since IMF governs the polar cap convection, we investigate here if this relationship can be explained by means of polar cap convection. The study is limited to equinox seasons. The poleward and equatorward edges of the trough are determined from satellite tomographic observations and their locations are plotted in magnetic coordinates together with the convection pattern given by Papitashvili and Rich [Papitashvili, V.O., Rich, F.J. High-latitude ionospheric convection models derived from DMSP ion drift observations and parameterized by the IMF strength and direction. J. Geophys. Res. 107, 2002, doi:10.1029/2001JA000264] using IMF measurements coincident with trough observations. The results indicate a close relationship between the troughs and convection. Most of the troughs are seen within the dusk cell and the pattern of trough observations rotates with the convection pattern, when B_y changes its sign. More dayside troughs are observed when B_z is negative than in the opposite case, i.e. fast convective flow favours the dayside trough occurrence. Nightside troughs are observed more frequently when B_y is negative. In both evening and morning sectors the trough is situated close to the edges of convection cells, which partly contradicts previous results showing that the troughs are associated with the convection reversal. It is concluded that plasma convection has an important role in trough generation, although the effect of a strong electric field and other mechanisms like precipitation certainly have a role of their own.     

The model of turbulent plasma sheet during IMF Bz > 0

Theory of the plasma sheet with medium-scale developed turbulence gives the possibility to explain the main processes of plasma sheet bifurcation and theta-aurora formation during IMF B_z > 0. The model suggests that during IMF B_z > 0 small bulge structure in the plasma sheet center is formed. The polarization of the bulge due to dawnward electron motion and duskward ion motion decreases the large-scale electric field in the bulge region. The decrease of the large-scale field in the conditions of constant coefficient of diffusion leads to the bulge growth. The results of plasma sheet bifurcation and theta-aurora formation modelling are presented.     

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.     

Large geomagnetic storms of extreme solar event periods in solar cycle 23

During extreme solar events such as big flares or/and energetic coronal mass ejections (CMEs) high energy particles are accelerated by the shocks formed in front of fast interplanetary coronal mass ejections (ICMEs). The ICMEs (and their sheaths) also give rise to large geomagnetic storms which have significant effects on the Earth’s environment and human life. Around 14 solar cosmic ray ground level enhancement (GLE) events in solar cycle 23 we examined the cosmic ray variation, solar wind speed, ions density, interplanetary magnetic field, and geomagnetic disturbance storm time index (D_st). We found that all but one of GLEs are always followed by a geomagnetic storm with D_st  −50 nT within 1–5 days later. Most(10/14) geomagnetic storms have D_st index  −100  nT therefore generally belong to strong geomagnetic storms. This suggests that GLE event prediction of geomagnetic storms is 93% for moderate storms and 71% for large storms when geomagnetic storms preceded by GLEs. All D_st depressions are associated with cosmic ray decreases which occur nearly simultaneously with geomagnetic storms. We also investigated the interplanetary plasma features. Most geomagnetic storm correspond significant periods of southward B_z and in close to 80% of the cases that the B_z was first northward then turning southward after storm sudden commencement (SSC). Plasma flow speed, ion number density and interplanetary plasma temperature near 1 AU also have a peak at interplanetary shock arrival. Solar cause and energetic particle signatures of large geomagnetic storms and a possible prediction scheme are discussed.     

Interconnection of high-latitude and low-latitude boundary layers when IMF BY is dominant

The antiparallel merging model places the location of the reconnection region for a dominant interplanetary magnetic field (IMF) B_Y at high latitudes at the dayside magnetopause and predicts that the low-latitude boundary layer (LLBL) is located on open field lines of the magnetospheric flanks. Interball-1 data obtained in the wide local time range near the low-latitude magnetopause makes it possible to analyze the LLBL plasma population and to find a link between possible reconnection at high latitudes and LLBL occurrence. We found that no boundary layer was observed in the regions which have no topological connection with the merging site. All cases of LLBL observations are located downstream from a specific boundary. This boundary coincides with the first magnetospheric field line touching the reconnection region and can be located in a wide local time region depending on the instant IMF direction. Even the LLBL on closed field lines shows the tendency to be concentrated in the vicinity of this boundary. Thus we show that all types of observed LLBLs are linked to reconnection sites predicted by the antiparallel merging model.     

INTERBALL statistical study of ion flow fluctuations in the plasma sheet

We use ion distribution measurements with CORALL instrument on-board the INTERBALL/Tail spacecraft to study plasma flows in the mid-tail (−9> X> −27 R_E) plasma sheet. Three velocity components computed every 2 minutes exhibit two types of velocity variations: Earthward bursty bulk flows (BBFs) and random flow fluctuations. Their properties are in a good agreement with the observations of the ISEE-2 spacecraft (Borovsky et al., 1997). The INTERBALL/Tail spacecraft configuration favors measurements of V_z component, in contrast to previous experiments in which only V_x and V_y were measured reliably. In the outer part of the plasma sheet V_y and V_z fluctuations were close to each other (variances σ(V_y) and σ(V_z) were about 160 and 110 km/s, respectively), but in the inner part at the dusk flank amplitude of V_y fluctuations increased and was 2 times higher than that of V_z component. This asymmetry of fluctuations should be taken into account during modern theoretical analysis and simulations.     

Current collection through the transport of electrons across magnetic field lines

Collection of electrons by a long conducting cylinder in a flowing plasma is studied by means of numerical simulations. The plasma flow simulates the relative motion between a spacecraft and plasma. The sheath structures and the levels of electron current collections for the cases with and without an ambient magnetic field ( ) are studied. It is found that for the flow perpendicular to the magnetic field, the current is considerably enhanced depending on the relative drift velocity. In the case of a non-zero magnetic field perpendicular to the cylinder axis, the potential structure is a two-dimensional double layer with dimensions L L_|, where L and L_| are the dimensions perpendicular and parallel to , respectively. L is found to be the current limiting radius given by the Parker-Murphy model. For the flow along , the electron current is found to be smaller than that for the flow perpendicular to . This is explained in terms of the potential structures.     

Collisionless reconnection in the magnetotail

Recent developments regarding collisionless reconnection in current sheets with a finite normal magnetic field component (B_z) are reviewed. In 2-D x, z configurations the ion tearing mode is stabilized by the electron compressibility. When the y dependence is included, cross-field current instabilities can be excited. Of these, the drift kink mode appears to be particularly important. 3-D electromagnetic particle simulations indicate that this mode can act as the precursor to the growth of tearing modes and subsequent reconnection.     

The collective gyration of a heavy ion cloud in a magnetized plasma

In both the ionospheric barium injection experiments CRIT I and CRIT II, a long-duration oscillation was seen with a frequency close to the gyro frequency of barium and a time duration of about one second. A model for the phenomenon which was proposed for the CRIT 1 experiment is here compared to the results from CRIT II which made a much more complete set of measurements. The model follows the motion of a low-β ion cloud through a larger ambient plasma. The internal field of the model is close to antiparallel to the injection direction v_i but slightly tilted towards the self-polarization direction E_P = −v_ixB. As the ions move across the magnetic field, the space charge is continuously neutralized by magnetic-field aligned electron currents from the ambient ionosphere, drawn by the divergence in the perpendicular electric field. These currents give a perturbation of the magnetic field related to the electric field perturbation by ΔE/ΔB ≈ V_a. The model predictions agree quite well with the observed vector directions, field strengths, and decay times of the electric and magnetic fields in CRIT II. The possibility to extend the model to the active region, where the ions are produced in this type of self-ionizing injection experiments, is discussed.     

两种用于NS-DSMC耦合方法的连续失效参数对比

NS-DSMC（Navier Stokes-Direct Simulation Monte Carlo）耦合方法是计算连续-稀薄跨流域流动的主要方法,应用过程中如何确定连续流域和稀薄流域的界面是此方法的关键问题之一,界面位置通常通过连续失效参数来判定.为合理选择连续失效参数,对目前广泛使用的两种连续失效参数Kn_Q和B参数进行了理论上的分析和比较,表明虽然他们建立的出发点不同,但在数学形式上具有一定相似性.通过圆柱绕流问题的数值试验进一步验证了两种连续失效参数具有一定对应关系.Kn_Q和B两种连续失效参数在NS-DSMC耦合方法中应用效果相差不大,但Kn_Q 的阈值0.05适用范围较广,而B参数的阈值对不同流动问题会有变化.      

Interplanetary shocks and geomagnetic activity during solar maximum (2000) and solar minimum (1995–1996)

Plasma and magnetic field parameter variations through fast forward interplanetary shocks were correlated with the peak geomagnetic activity index Dst in a period from 0 to 3 days after the shock, during solar maximum (2000) and solar minimum (1995–1996). Solar wind speed (V) and total magnetic field (B_t) were the parameters with higher correlations with peak Dst index. The correlation coefficients were higher during solar minimum (r² = 56% for V and 39% for B_t) than during solar maximum (r² = 15% for V and 12% for B_t). A statistical distribution of geomagnetic activity levels following interplanetary shocks was obtained. It was observed that during solar maximum, 36% and 28% of interplanetary shocks were followed by intense (Dst  −100 nT) and moderate (−50  Dst < −100 nT) geomagnetic activity, whereas during solar minimum 13% and 33% of the shocks were followed by intense and moderate geomagnetic activity. It can be concluded that the upstream/downstream variations of V and B_t through the shocks were the parameters better correlated with geomagnetic activity level, and during solar maximum a higher relative number of interplanetary shocks can be followed by intense geomagnetic activity than during solar minimum. One can extrapolate, for forecasting goals, that during a whole solar cycle a shock has a probability of around 50% to be followed by intense/moderate geomagnetic activity.     

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”.     

AL index dependence on the solar wind input revisited

Estimates of the geomagnetic indices made with real-time solar wind measurements form the basis of many space weather forecast techniques. We analyze 20 years of hourly AL and OMNI solar wind data to determine geomagnetic importance of various solar wind and IMF parameters. Besides the solar wind driver of primary importance (VB_s), the first-order contributions, significantly increasing the quality of the model are: solar wind velocity, 2 h of solar wind prehistory and 1 h of AL history. The factors of secondary importance, marginally contributing to overall statistical quality, are IMF B_y, solar wind density, and IMF fluctuations. The dynamic pressure is geomagnetically effective only if the pressure is lower than the average. Modelling of the same data set with an artificial neural network (ANN) confirmed our selection of important factors. Statistically the ANN model was just marginally better than our analytic expression . The AU index dependence is principally different from AL in several respects; therefore modelling of the AE composite index is physically misleading.     

行星际磁场对磁尾场向电流发生率的影响

利用Cluster卫星的磁场和等离子体探测数据, 研究了行星际磁场(IMF)时钟角(clock angle) Φ和锥角(cone angle) θ对磁尾等离子体片边界层(PSBL)区场向电流发生率的影响. 当时钟角Φ >0时, 磁尾场向电流 的发生率较高, 这表明磁尾场向电流的发生与昏向太阳风条件更为密切; 当 90°<|Φ|<180°时, 场向电流的发生率较高, 这表明 场向电流的发生与南向IMF更为密切. 当锥角θ <30°时(即IMF与 日地连线夹角较小时)场向电流的发生率较低. 而当θ> 30°时, 场向电流在90°<|Φ|<180°的情况下发生率明显增大, 这说明南向IMF情况下, 场向电流发生率明显增大. 但是当|Φ|<90°时 (北向IMF情况下), 尽管θ很大, 场向电流的发生率并未明显增大. 当θ>70°时, 且在140°< < i>Φ<160°的行星 际磁场条件下, 磁尾等离子体片边界层区场向电流的发生率最大.      

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.     

近地尾向流和磁场偶极化——TC-1卫星观测结果

2004年10月12日, 在01:30---04:30 UT期间, 位于向阳侧磁层顶附近的Geotail卫星探测到行星际磁场为持续南向. 此太阳风条件驱动了一个小磁暴, Sym-H指数在04:12 UT达到最小值-33nT. 在磁暴主相期间, AE指数维持在较高的水平, 其最大值达400nT. 02:00---03:00 UT期间, TC-1卫星在近地磁尾(-10.6, 3.2, -0.1)Re处观测到明显的亚暴膨胀相特征和磁场偶极化过程. 在偶极化前1min, 有较强的(v_x<-100 km/s)持续时间超过3min的尾向流发生. 分析发现该尾向流具有低温、高密度和沿磁场流动的特点, 这说明尾向流具有来源于电离层风的特征. 尾向流期间, TC-1观测的磁场分量Bx和总的磁场强度增加, 磁倾角减小, 磁场结构变成非偶极型, 说明尾向流对磁场结构有一定的影响, 文中尝试给出了相应的物理解释. 观测表明, 该事例中的近地磁尾尾向流可能对磁场偶极化过程的发生有重要意义.      

镜模波识别方法研究及其在火星磁鞘中的应用

镜模波是温度各向异性等离子体中的一种波动结构,根据磁场和离子分布及波动特性可以进行识别.本文对比了只使用磁场数据与同时使用磁场及离子数据两种识别方法,分析了两类方法的特点.只使用磁场数据的方法基于磁场强度变化大、方向沿背景磁场的特征,通常使用磁场强度的波动幅度ΔB/|B|以及磁场变化方向与背景磁场的夹角θ_min,θ_max作为参数;同时使用磁场及粒子数据的方法利用的是磁场纵波特性、总压平衡和波动在等离子体坐标系下静止的特征.使用两种方法对MAVEN卫星在火星磁鞘内的数据进行识别,结果表明在某些情况下,只使用磁场数据会导致对镜模波的误判.通过研究改变上述参数阈值时识别结果的变化,发现当θ_min> 40°,θ_max < 40°,ΔB/|B|> 80%时,只用磁场数据可取得较好的识别效果.      

Auroral signatures of transient processes in the outer magnetosphere

We present an analysis of sporadic and recurrent injections of magnetospheric ions in the midnight auroral oval during substorms and of the associated ionospheric ion outflows. The source of plasma sheet precipitating ions is determined using a simple method, based on the measured relation between the ion inverse velocity and time (l = v × t). This method is applied here to two typical passes of the Interball-Auroral (IA) satellite at distances of 3 R_E above the auroral regions. Substorm related ion injections are shown to be mainly due to time of flight effects. In contrast with particle trajectory computations (Sauvaud et al., 1999), the inverse velocity method does not require magnetic and electric field models and can thus be used systematically for the detection of time of flight dispersed ion structures (TDIS). This allowed us to build a large database of TDIS events and to perform a statistical analysis of their spatial distribution. For the cases presented here the source region of the injected ions is found at radial distances from 18 to 30 R_E near the equatorial magnetosphere. At Interball altitudes ( 3 R_E), ion injections detected at the poleward boundary of the nighside auroral oval are associated with shear Alfvén waves superimposed over large-scale quasi-static current structures. We show that the most poleward TDIS are collocated with a large outflow of ionospheric H⁺ and O⁺ displaying pitch-angle distributions peaked in the pitch-angle range 90°–120°. These ions are thus accelerated perpendicularly to the magnetic field not only in the main auroral acceleration region but also up to at least 3 R_E. The expanding auroral bulge thus constitutes a significant source of H⁺ and O⁺ ions for the mid-tail magnetosphere.     

Statistical characteristics of low latitude Pc4s: Influence of solar wind and IMF parameters

Pc4 signatures for the year 2013, extracted from geomagnetic north–south and east–west components of induction coil magnetometer (LEMI 30) from low latitude station Desalpar (DSP), operated by Institute of Seismological Research (ISR), India have been investigated vis-à-vis the prevalent interplanetary parameters (IMF) as well as the geomagnetic activity indices. A clear dominance of Pc4-5 (467 events) over Pc3 (17 events) is observed. Local time variation of Pc4 shows a peak in the noon sector in both X and Y components. Our investigations show that the dominant peak frequency is 10 mHz at low latitude region. Correlations with solar wind and IMF parameters illustrate highest occurrence of Pc4 for a solar wind speed of 300–400 km/s and average IMF B field of 3–6 nT. The amplitude of Pc4s at DSP shows an increase with increasing solar wind speed, plasma density, solar wind dynamic pressure and average B field which is also reflected in the trend of frequency variation of these pulsations. We report that IMF clock angle at low latitude does not have influence on Pc4 occurrence. Based on the characteristics of these events, detected in latitudinally distributed stations from low and mid-latitudes from northern and southern hemisphere, we infer that modes were compressional, which could be driven by K-H instability or solar wind dynamic pressure, as compressional modes can propagate to low latitude with little attenuation.