Near Earth Vortices Driving of Field Aligned Currents Based on Magnetosphere Multiscale and Swarm Observations

A long-standing mystery in the study of Field-Aligned Currents (FACs) has been that:how the currents are generated and why they appear to be much stronger at high altitudes than in the ionosphere.Here we present two events of magnetotail FACs observed by the Magnetospheric Multiscale Spacecraft (MMS) on 1st July and 14th July 2016,to show how the Substorm Current Wedges (SCW) were formed.The results show that particles were transferred heading towards the Earth during the expansion phase of substorms. The azimuthal flow formed clockwise (counter-clockwise) vortex-like motion,and then generated downward (upward) FACs on the tailward/poleward side of the distorted field with opposite vorticity on their Earthward/equatorward side.We also analyzed the Region 1 FACs observed by the Earth Explorer Swarm spacecraft on 1st July 2016 and found that they were associated with FACs observed by MMS,although differing by a factor of 10.This difference suggests that either there was the closure of the currents at altitudes above 500 km or the currents were not strictly parallel to B and closed at longitudes away from where they were generated.      

大尺度场向电流的控制因素

将ＩＳＥＥ－１和ＩＳＥＥ－２飞般测量的地球内磁层场向电流作为行星际磁场与极区地磁活动水平（由ＡＬ表征）的函数。发现大约７５％的场向电流发生在行星际磁场南向时,其余２５％发生在行星际磁场由南向转到北向的半小时内。而且,发生在ＡＬ〈－１５０ｎＴ的场向电流也大约是７５％。场向电流的强度和密度随行星际磁场南向分量及ＡＬ的负值增加而增大。由此而得出结论,内磁层场向电流的产生主要是由行星际磁场控制的,是太阳风     

Controlling factors of Region 2 field-aligned current and its relationship to the ring current: Model results

One essential component of magnetosphere and ionosphere coupling is the closure of the ring current through Region 2 field-aligned current (FAC). Using the Comprehensive Ring Current Model (CRCM), which includes magnetosphere and ionosphere coupling by solving the kinetic equation of ring current particles and the closure of the electric currents between the two regions, we have investigated the effects of high latitude potential, ionospheric conductivity, plasma sheet density and different magnetic field models on the development of Region 2 field-aligned currents, and the relationship between R2 FACs and the ring current. It is shown that an increase in high latitude potential, ionospheric conductivity or plasma sheet density generally results in an increase in Region 2 FACs’ intensity, but R2 FACs display different local time and latitudinal distributions for changes in each parameter due to the different mechanisms involved. Our simulation results show that the magnetic field configuration of the inner magnetosphere is also an important factor in the development of Region 2 field-aligned current. More numerical experiments and observational results are needed in further our understanding of the complex relationship of the two current systems.     

Hemisphere-coupled modeling of nighttime medium-scale traveling ionospheric disturbances

Nighttime medium-scale traveling ionospheric disturbances (MSTIDs), which have tilted frontal structures in the midlatitude ionosphere, are investigated by the midlatitude ionosphere electrodynamics coupling (MIECO) model in this study. It has been proposed that the electrodynamic coupling between the E and F regions plays an important role in generating MSTIDs within a few hours. An intriguing aspect of MSTIDs is that they were simultaneously observed at magnetic conjugate locations in the Northern and Southern Hemispheres. In order to study the hemisphere-coupled electrodynamics, the MIECO model has been upgraded to consist of two simulation domains for both hemispheres in which the electrostatic potential is solved by considering electrodynamics in both hemispheres. The simultaneous occurrence of MSTIDs at the magnetic conjugate stations has clearly been reproduced when the F-region neutral wind satisfies the unstable condition in both hemispheres and a sporadic-E layer is given only at the Northern (summer) Hemisphere. Even if the unstable condition is satisfied in the summer hemisphere, an unfavorable F-region neutral wind in the winter hemisphere largely suppresses the growth of MSTIDs in both hemispheres.     

On the effect of near-equatorial thunderstorms on the global distribution of ionospheric potential

We develop our earlier attempts to perform an indirect quantitative examination of the hypothesis that electric currents flowing up from thunderstorms to the ionosphere (also known as Wilson currents) charge the ionosphere to a large positive potential with respect to the Earth. First, we take the electrostatic potential arising from the interaction of the solar wind with the Earth’s magnetosphere derived from an experimental data-based model of the high-latitude field-aligned currents. We then obtain the global distribution of ionospheric potential, utilizing a thin shell model, based on integration along field lines of the current continuity equation with a realistic model of ionospheric conductivity. Next, we include additional upward currents to simulate the effect of the three main thunderstorm regions over equatorial Asia/Oceania, Africa and the Americas. We compare the local time variation of the eastward electric field in the ionosphere produced by these three equatorial sources separately, and seek to understand the substantial differences between them. Finally, we examine the variation with local time of the eastward electric field in the ionosphere at low latitudes.     

Remote sensing of the ring current using energetic neutral atoms

Energetic neutral atom (ENA) images of the storm-time ring current obtained from the ISEE-1 spacecraft provide information for a “zero-order” global model of the energetic ion distribution. With the assumption of isotropic pressure and magnetostatic, non-convective pressure balance, the global system of electrical currents driven by the ion pressure can be calculated using Euler potentials for the divergenceless current density. Radial pressure gradients drive azimuthal currents, and azimuthal pressure gradients drive radial currents. The radial currents cause current lines in the inner magnetosphere to close in the ionosphere, forming a “partial” ring current. The intensities and locations of these field-aligned currents driven into and out of the ionosphere resemble those of the observed Region 2 current system, but not all observed properties of the Region 2 system are reproduced by the “zero-order” model.     

场向电流随亚暴位相的变化

利用ISEE1和2卫星测量的磁场数据,计算了电离层中的场向电流。依据每个场向电流事件所伴随的亚暴位相,分别计算了一区和二区场向电流强度、密度及电流片厚度在亚暴成长相、膨胀相和恢复相的平均值及中间值。其结果,从成长相到膨胀相,一区和二区场向电流的强度和密度增加,从膨胀相到恢复相,其值减小。平均说来,一区的电流强度约是二区的1.4倍。电流片厚度的变化在上述期间内与电流强度及密度的变化趋势相反。     

Investigation of the planar midlatitude ionospheric trend under different levels of solar activity

A better understanding of the ionosphere through accurate mathematical models is no doubt a crucial element. This study focuses on the challenging problem of building a model representing the complex structure of the midlatitude ionosphere. Previous studies have shown that a regional planar model is suitable in representing the total electron content (TEC) trend in the midlatitude ionosphere in both hemispheres. In this study, the planar trend model for 12 non-overlapping northern hemisphere regions in three groups of geographically near 4 regions is further investigated under different levels of solar activity; low, moderate and high. To that end, the coefficients of the model are estimated in the least squares sense using total electron content values from global ionospheric maps (GIMs) for the years 2009, 2012 and 2014. Subsequently, these coefficients are used to reconstruct estimated TEC maps which are then compared with actual GIM-TEC by investigating their difference in normalized $L_2$ norm squared sense. The regional planar trend model provides a particularly successful representation in the years 2012 and 2014 for which the solar activity level is the dominant factor determining the TEC trend. Under low solar activity conditions of 2009, other factors such as ocean currents, temperature variations and meteorological phenomena are suspected to have a considerable effect in some regions depending on their geographic location and on seasonal trends in those regions. As an example, studies show that under the influence of the Pacific Decadal Oscillation (PDO) and Siberian High (SH), a significant cooling trend between 2004 and 2018 in autumn is observed in Eurasia, which, in conjunction with the low solar activity levels, may be related to the deviations from the actual GIM-TEC in 2009 in these regions. As solar radiation increases, however, such bottom-side forcings are masked in 2012 and 2014 and these deviations are no longer observed.     

Field-aligned plasmaspheric flows at moderate latitudes

In a model of the plasmasphere, coupled time-dependent continuity, momentum and energy equations are solved for thermal O⁺, H⁺ and electrons. The field-aligned mass flow coupling and thermal coupling of conjugate ionsopheres via the protonosphere are studied. For solstice conditions, thermal coupling between conjugate hemispheres gives rise to very strong upward flows of O⁺ in the topside ionosphere of the summer hemisphere at the time of sunrise in the conjugate (winter) ionosphere; a less marked effect (but with downward flow) occurs in the summer ionosphere at winter sunset. In addition, there are strong upward and downward flows of O⁺ at local sunrise and sunset, respectively, in both hemispheres. At both L = 1.5 and L = 3, the 24-hour time-integrated interhemispheric H⁺ flux is in the summer - winter direction. At L = 1.5 its magnitude is in good agreement with the magnitude of the time- integrated field-aligned plasma (O⁺ + H⁺) flux at 1000 km altitude; there is no such agreement at L = 3.     

火星空间大尺度电流体系分布特征

在MHD模型基础上, 对火星空间环境的电流分布进行了模拟. 结果表明, 火星空间存在着弓激波电流、磁堆积区电流、电离层电流和磁尾电流. 弓激波电流在激波曲面上均由北向南自成体系, 电流密度在弓激波顶区域较大. 在向阳面磁堆积区边界电流与电离层电流彼此耦合形成完整的回路, 在背阳面磁堆积区边界电流与磁尾中心电流片耦合形成完整的回路.     

火星磁场产生的原因及其分布

构建了一个可以得到火星赤道面上磁场分布的模型. 模型根据卫星观测数据, 提出了火星电离层、磁层顶和磁尾电流片上都各自通有电流的假设. 由电流的连续性条件可知, 这三种背景条件下的电流之间满足一定关系, 即火星磁层顶上的总电流是电离层上的总电流与磁尾电流片上的总电流之和. 这些电流产生的磁场与太阳风磁场共同构成了火星赤道面上磁场分布. 通过计算发现, 采用这种磁场模型得到的结果与目前卫星所观测的结果以及与采用其他方法得到的结果符合得较好.     

Magnetic field disturbances in the mid-altitude cusp

Magnetic field disturbances and hot particles in the cusp as seen mainly by both the IMAP-3 magnetometer and PROMICS -3 spectrometer aboard the INTERBALL-AURORAL PROBE (AU) (perigee 4 Re, inclination 65 deg) are discussed. Orbits from March, 1997 are analysed accordingly. The INTERBALL-AU transits cross the middle of the cusp from low to high latitude from pre-noon to post-noon hours 11–13 MLT. Although the crossings are not exactly meridional, they reach 80–81 ILAT and permitting detection of the cusp, its subregions and their equatorward and poleward boundaries. Magnetic field disturbances are analysed for field-aligned currents (FACs) during different IMF conditions. Regions with structured forms are emphasised, in particular where a more intense current is concentrated. The following cases are compared: (a) the most intense current in the boundary cusp, adjacent to the ring current (March 7, southward IMF); (b) FACs are concentrated completely inside the cusp (March 13, southward -downward IMF) and (c) FACs with pronounced intensities equatorward from the cusp proper (March 11, duskward IMF).     

Electric currents from thunderstorms to the ionosphere during a solar cycle: Quasi-static modeling of the coupling mechanism

Conduction and displacement currents, and their sum the Maxwell current, generated over a thunderstorm (TS) with recurrent lightning discharges are investigated theoretically. The aim is to study better the influence of different factors on these currents, which form the link between thunderstorms and the ionosphere in the global atmospheric electrical circuit. The factors studied concern the thunderstorm characteristics (the charge separation current, and the lightning discharge parameters), as well as the atmospheric and cloud conductivity. Some of these factors may show long-term changes with the 11-year solar cycle, possibly realized through an inverse dependence of the cosmic ray flux on solar activity. Earlier investigations have suggested that the lightning-related charge redistribution and subsequent relaxation, rather than the high intensity current, is mainly the source of the energy coupled to the ionosphere. With respect to this, a quasi-electrostatic analytical model is proposed, based on Maxwell’s equations. The currents are generated by a TS modeled as a positive vertical dipole with charges which are first accumulated and then destroyed by lightning. Our computations show that the mean and peak values of the conduction and total Maxwell currents to the ionosphere depend significantly on the charge moment change. The mean currents are also sensitive to the reduction of the conductivity in thunderclouds. Small variations of the stratospheric conductivity (20% at geomagnetic latitude 40° and 40–50% at 55°) with the solar activity do not influence the currents to the ionosphere very much.     

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

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

Sq（Y）的季节变化和场向电流

本文对中国东部地磁台站链地磁静日变化东向分量进行了分析,有结果表明,在冬、夏季清晨,南北半球之间存在着电离层发电机驱动的场向电流,其方向是从夏季半球到冬季半球;在黄昏,也可能存在方向相反的电流。     

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.     

Occurrence of medium-scale travelling ionospheric disturbances identified in the Tasman international geospace environment radar observations

The occurrence characteristics of medium-scale travelling ionospheric disturbances (MSTIDs) were investigated using the Tasman International Geospace Environment Radar (TIGER). From the occurrence study of sea echoes, we found two maxima, one pre-noon and the other after noon. They are less obvious with increase of magnetic activities, and more obvious when B_z is northwards. It is suggested that this maxima were related to fore- and after-noon maxima in the distribution of net field-aligned currents flowing from the magnetosphere to the ionosphere, and that these two regions were sources of atmospheric gravity waves (AGWs) due to enhancement of Hall conductivities in the ionosphere. The Lorentz force is suggested to be a possible mechanism for the excitation of MSTIDs in the dayside ionosphere.     

Adaptation of the bottomside electron density profile of the IRI model to data of topside radiosounding

A method is proposed for reconstructing the electron density profiles N(h) of the IRI model from ionograms of topside satellite sounding of the ionosphere. An ionograms feature is the presence of traces of signal reflection from the Earth's surface. The profile reconstruction is carried out in two stages. At the first stage, the N(h) –profile is calculated from the lower boundary of the ionosphere to the satellite height (total profile) by the method presented in this paper using the ionogram. In this case, the monotonic profile of the topside ionosphere is calculated by the classical method. The profile of the inner ionosphere is represented by analytical functions, the parameters of which are calculated by optimization methods using traces of signal reflection, both from the topside ionosphere and from the Earth. At the second stage, the profile calculated from the ionogram is used to obtain the key parameters: the height of the maximum hmF2 of the F2 layer, the critical frequency foF2, the values of B0 and B1, which determine the profile shape in the F region in the IRI model. The input of key parameters, time of observation, and coordinates of sounding into the IRI model allows obtaining the IRI-profile corrected to real experimental conditions. The results of using the data of the ISIS-2 satellite show that the profiles calculated from the ionograms and the IRI profiles corrected from them are close to each other in the inner ionosphere and can differ significantly in the topside ionosphere. This indicates the possibility of obtaining a profile in the inner ionosphere close to the real distribution, which can significantly expand the information database useful for the IRTAM (IRI Realmax Assimilative Modeling) model. The calculated profiles can be used independently for local ionospheric research.     

The solar wind control of the equatorial ionosphere dynamics during geomagnetic storms

The interplanetary magnetic field, geomagnetic variations, virtual ionosphere height h′F, and the critical frequency foF2 data during the geomagnetic storms are studied to demonstrate relationships between these phenomena. We study 5-min ionospheric variations using the first Western Pacific Ionosphere Campaign (1998–1999) observations, 5-min interplanetary magnetic field (IMF) and 5-min auroral electrojets data during a moderate geomagnetic storm. These data allowed us to demonstrate that the auroral and the equatorial ionospheric phenomena are developed practically simultaneously. Hourly average of the ionospheric foF2 and h′F variations at near equatorial stations during a similar storm show the same behavior. We suppose this is due to interaction between electric fields of the auroral and the equatorial ionosphere during geomagnetic storms. It is shown that the low-latitude ionosphere dynamics during these moderate storms was defined by the southward direction of the B_z-component of the interplanetary magnetic field. A southward IMF produces the Region I and Region II field-aligned currents (FAC) and polar electrojet current systems. We assume that the short-term ionospheric variations during geomagnetic storms can be explained mainly by the electric field of the FAC. The electric fields of the field-aligned currents can penetrate throughout the mid-latitude ionosphere to the equator and may serve as a coupling agent between the auroral and the equatorial ionosphere.     

Equatorial F-region ionization differences between March and September, 1979

Ionospheric total electron content and the F-region maximum electron density at a number of stations in the equatorial region, during the recent solar activity maximum period 1979 to 1980, show significant differences between the two equinoctial periods. Ionization during the month of March is higher than in September, irrespective of the station location both in northern and southern hemispheres, and in different longitude sectors. The observed pattern is compared with those predicted by different models, in particular with one of the authors which includes processes such as ionization production, loss, electrodynamic drifts, winds and global composition changes involved in the equatorial ionosphere. It is found that a change in the neutral composition is primarily responsible for the observed F-region density differences between March and September.