Understanding magnetotail current sheet meso-scale structures using MHD simulations

Recent Cluster observations have strongly supported the existence of meso-scale structure in the magnetotail current sheet. In our study, a magnetohydrodynamic simulation event study exhibited current sheet behavior comparable to that seen in the Cluster observations. Geotail and DoubleStar observations also show that the simulation is providing a realistic representation of the magnetosphere during the period of interest; that is, when the current sheet evidently becomes bifurcated. The magnetohydrodynamic simulation allows us to place the local observations into a global contest. It shows that the observations can be explained in terms of localized reconnection tailward of the Cluster location and the formation of a flux rope nearby. The simulation also features wave-like structure across the current sheet.     

Tearing,coalescence and fragmentation processes in solar flare current sheet and drifting pulsating structures

The paper presents a summary of results from two different simulations which study the tearing, coalescence and fragmentation of current sheets, the associated production of energetic electrons and of plasma waves from these electrons which could explain drifting pulsation structures observed at radio wavelengths. Using a 2.5-D particle-in-cell (PIC) model of the current sheet it is shown that due to the tearing mode instability the current sheet tears into plasmoids and these plasmoids later on coalesce into larger ones. During these processes electrons are accelerated and they produce observable electromagnetic waves. Furthermore, the 3-D PIC model with two current sheets extended in the electric current direction shows their fast fragmentation associated with the exponential dissipation of the free magnetic field energy. An example of the drifting pulsating structure which is considered to be a radio signature of the above mentioned processes in solar flares is shown.     

Kinetic theory of the current sheath. I. On polarization of an equilibrium current sheath

The basic theme of this paper is the investigation of the polarizing electric fields caused by anisotropy of thermal plasma of a current sheath. The paper provides a generalization of the equilibrium distribution function of the current sheath (Harris function) to the case of temperature anisotropy of plasma along and across the sheath. It appears that solution of the electroneutral equilibrium represents a narrow class of solutions and is true only for certain relations between the parameters of the problem. In the general case, the plasma of a stationary current sheath is polarized. The structure of the current sheath of a magnetosphere tail is investigated with regard to the effect of polarization, and profiles of the polarization electric fields are obtained. This field should be taken into account in the study of current sheath stability.     

Evaluation of substorm models with Cluster observations of plasma flow reversal in the magnetotail

We evaluate two prevailing substorm models with an event of plasma flow reversal from tailward to Earthward detected by Cluster at the downstream distance of ∼19 R_E in the magnetotail during a substorm on August 22, 2001. We use the unique capability of Cluster measurements in determining gradients to examine the associated current density, Lorentz force, and current dissipation/dynamo term. In association with plasma flow reversal, it is found that (1) there was no clear quadrupole magnetic perturbation signature, (2) the x-component of the Lorentz force did not change sign, (3) the y-component of the product of the current density and the electric field was occasionally negative indicative of a dynamo effect, and (4) the timing sequence of flow reversal from the Cluster configuration did not match tailward motion of a single plasma flow source. These observations are consistent with the near-Earth initiation model for substorms with multiple current disruption sites moving progressively tailward near the late stage of substorm expansion.     

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.     

Fluxes and nuclear abundances of cosmic rays inside the magnetosphere using a transmission function approach

At 1 AU and outside the Earth’s magnetosphere, the relative abundances to protons for He (He/p), C (C/p) and Fe (Fe/p) nuclei were calculated using the observation data of AMS-01 (for p and He) and HEAO-3 (for C and Fe) above 0.8 GeV/nucleon. In addition, the transmission function (TF) for the GCR propagation inside the magnetosphere was evaluated using the IGRF and T96 (introduced by Tsyganenko and Stern) models to obtain permitted and forbidden trajectories inside the magnetosphere. The TF allowed one to derive the primary He-nuclei fluxes in the same geomagnetic regions of AMS-01 observations. These fluxes were found in good agreement with the observation data. Furthermore inside the magnetosphere in addition to the flux of helium, it allowed one to obtain those of the primary p, C, and Fe nuclei and the relative abundances of He, C and Fe nuclei to protons from the same observation data of AMS-01 and HEAO-3 above ≈0.8 GeV/nucleon. Up to a geomagnetic latitude of ≈45.84°, the relative isotopic abundances were found to depend on the mass number I_isot and, on average, range from a factor ≈2.31 up to ≈3.35 larger than those outside the magnetosphere at 1 AU. Thus, the magnetospheric isotopic/nuclear relative abundances differ from those inside the solar cavity and those in the interstellar space. The usage of the TF approach can allow one to determine the nuclear abundances in the magnetosphere at any geomagnetic latitude and, thus, any orbit, provided that the CR spectra are determined at 1 AU.     

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.     

Outstanding issues in understanding the dynamics of the inner plasma sheet and ring current during storms and substorms

The paper deals with five selected issues of the dynamical coupling of the near-Earth plasma sheet and magnetosphere, (1) substorm initiation, (2) dipolarization, (3) pressure release of the outer magnetosphere via the auroral energy conversion process, (4) magnetization of the very high beta plasma assembling at the inner edge of the tail, and (5) penetration of energetic particles into the ring current below L 4. One outstanding and strongly debated subject is not discussed here, the origin of the substorm current wedge. The main conclusions (or personal preferences) are: (1) the substorm is initiated by formation of a near-Earth neutral line; (2) dipolarization occurs through magnetic flux transport by the earthward reconnection flow and not by current diffusion; (3) the auroral energy conversion process, the “auroral pressure valve”, contributes substantially to the pressure release during the substorms; (4) high beta ( 10) plasma breaks up into smaller scale blobs under slow magnetization; and (5) deep and prolonged penetration of hot plasma sheet plasma into the middle magnetosphere produces currents and electric fields which lead to the growth of the storm-time ring current.     

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.     

Geomagnetic effects on mid-latitude railways: A statistical study of anomalies in the operation of signaling and train control equipment on the East-Siberian Railway

Intense geomagnetically induced currents (GIC) can hamper rail traffic by disturbing signaling and train control systems. GIC threats have been a concern for technological systems at high-latitude locations due to geomagnetic disturbances driven by substorm expansion electrojet or convection electrojet intensifications. However, other geomagnetic storm processes such as storm sudden commencement (SSC) and geomagnetic pulsations can also cause GIC concerns for technological systems. We present in this paper the first evidence based on statistical data for links between geomagnetic disturbances and faulty operations (anomalies) in the functioning of railway automatics and telemetry. We analyze anomalies of automatic signaling and train control equipment which occurred in 2004 on the East-Siberian Railway (corrected geomagnetic latitude _m = 46–51°N and longitude λ_m = 168–187°E). Our results reveal a seasonal effect in the number of anomalies per train similar to the one observed in geomagnetic activity (Kp, Ap, Dst indices). We also found an increase by a factor of 3 in the total duration of daily anomalies during intense geomagnetic storms (local geomagnetic index specific to Siberian Observatory A_max > 30), with a significant correlation between the daily sum of durations of anomalies with geomagnetic activity. Special attention was paid to failures not related to recognized technical malfunctions. We found that the probability of these failures occurring in geomagnetically disturbed periods was 5–7 times higher than the average anomaly occurrence.     

Relative contributions of terrestrial and solar wind ions in the plasma sheet

A major uncertainty concerning the origins of plasma sheet ions is due to the fact that terrestrial H⁺ can have similar fluxes and energies as H⁺ from the solar wind. The situation is especially ambiguous during magnetically quiet conditions (AE < 60γ) when H⁺ typically contributes more than 90% of the plasma sheet ion population. In this study we examine that problem using a large data set obtained by the ISEE-1 Plasma Composition Experiment. The data suggest that one component of the H⁺ increases in energy with increasing activity, roughly in proportion to $\text{1}\text{4}$ the energy of the He⁺⁺, whereas the other H⁺ component has about the same energy at all activity levels, as do the O⁺ and the He⁺. If we can assume that the H⁺ of solar wind origin on the average has about the same energy-per-nucleon as the He⁺⁺, which is presumably almost entirely from the solar wind, then the data imply that as much as 20–30% of the H⁺ can be of terrestrial origin even during quiet conditions.     

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.     

Temporal variation of the arterial pressure in healthy young people and its relation to geomagnetic activity in Mexico

We present a study of the temporal behavior of the systolic (SBP) and diastolic (DBP) blood pressure for a sample of 51 normotensive, healthy volunteers, 18 men and 33 women with an average age of 19 years old in Mexico City, Mexico, during April and May, 2008. We divided the data by sex along the circadian rhythm. Three geomagnetic storms occurred during the studied time-span. The strongest one, a moderate storm, is attributed to a coronal hole border that reached the Earth. The ANOVA test applied to the strongest storm showed that even though we are dealing with a moderate geomagnetic storm, there are statistically significant responses of the blood pressure. The superposed epoch analysis during a three-day window around the strongest storm shows that on average the largest changes occurred for the SBP. Moreover, the SBP largest increases occurred two days before and one day after this storm, and women are the most sensitive group as they present larger SBP and DBP average changes than men. Finally, given the small size of the sample, we cannot generalize our results.     

Ballooning perturbations in the inner magnetosphere of the Earth: Spectrum,stability and eigenmode analysis

We investigate the generation of ballooning perturbations in the inner magnetosphere of the Earth in the dipole model of the geomagnetic field taking into account ionospheric boundary conditions. The ionosphere is considered as a thin layer with finite conductivity. The eigenmode spectrum is discrete and consists of Alfvén, slow magnetosonic, flute and incompressible modes. Their interaction depends on ionospheric conductivity. The decay rate is small in noon and night sectors and large in dawn and dusk sectors. The lowest stability threshold α/γ ≈ 4.25 is determined by flute modes.     

Modeling of ionospheric parameter variations in East Asia during the moderate geomagnetic disturbances

The results of modeling of ionospheric disturbances observed in the East Asian region during moderate storms are presented. The numerical model for ionosphere–plasmasphere coupling developed at the ISTP SB RAS is used to interpret the data of observations at ionospheric stations located in the longitudinal sector of 90–130°E at latitudes from auroral zone to equator. There is obtained a reasonable agreement between measurements and modeling results for winter and equinox. In the summer ionosphere, at the background of high ionization by the solar EUV radiation in the quiet geomagnetic period the meridional thermospheric wind strongly impacts the electron concentration in the middle and auroral ionosphere. The consistent calculations of the thermospheric wind permit to obtain the model results which are closer to summer observations. The actual information about the space-time variations of thermosphere and magnetosphere parameters should be taken into account during storms.     

Response of the ionosphere to super geomagnetic storms: Observations and modeling

The relative importance of the main drivers of positive ionospheric storms at low-mid latitudes is studied using observations and modeling for the first time. In response to a rare super double geomagnetic storm during 07–11 November 2004, the low-mid latitude (17°–48°N geomag. lat.) ionosphere produced positive ionospheric storms in peak electron density (NmF2) in Japan longitudes (≈125°–145°E) on the day of main phase (MP1) onset (06:30 LT) and negative ionospheric storms in American longitudes (≈65°–120°W) on the following day of MP1 onset (13:00–16:00 LT). The relative effects of the main drivers of the positive ionospheric storms (penetrating daytime eastward electric field, and direct and indirect effects of equatorward neutral wind) are studied using the Sheffield University Plasmasphere Ionosphere Model (SUPIM). The model results show that the penetrating daytime (morning–noon) eastward electric field shifts the equatorial ionisation anomaly crests in NmF2 and TEC (total electron content) to higher than normal latitudes and reduces their values at latitudes at and within the anomaly crests while the direct effects of the equatorward wind (that reduce poleward plasma flow and raise the ionosphere to high altitudes of reduced chemical loss) combined with daytime production of ionisation increase NmF2 and TEC at latitudes poleward of the equatorial region; the later effects can be major causes of positive ionospheric storms at mid latitudes. The downwelling (indirect) effect of the wind increases NmF2 and TEC at low latitudes while its upwelling (indirect) effect reduces NmF2 and TEC at mid latitudes. The net effect of all main drivers is positive ionospheric storms at low-mid latitudes in Japan longitude, which qualitatively agrees with the observations.     

The effect of the longitudinal propagating waves on the electron acceleration in the reconnecting current sheet

Based upon the most efficient electron acceleration near the midplane of 3D non-neutral driven reconnecting current sheet (RCS) and the electrostatic wave excitation by the drift Maxwellian distribution of electrons in Vlasov simulation, we assume that the electrostatic waves mainly propagate opposite to the reconnecting electric field and investigated how these waves affect the electron acceleration. The main results are: (1) when the electron’s velocity equals to the phase speed of the waves, they will be trapped and have the different accelerating characteristics from the untrapped electrons through solving the momentum equations of electrons analytically; (2) the test particle simulations further prove that the number of the energetic electrons decreases with the increasing intensity of unstable waves, and the distribution of the energetic electrons takes on the double power-law.     

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.     

Longitudinal variations of geomagnetic and ionospheric parameters in the Northern Hemisphere during magnetic storms according to multi-instrument observations

We present a joint analysis of longitude-temporal variations of ionospheric and geomagnetic parameters at middle and high latitudes in the Northern Hemisphere during the two severe magnetic storms in March and June 2015 by using data from the chains of magnetometers, ionosondes and GPS/GLONASS receivers. We identify the fixed longitudinal zones where the variability of the magnetic field is consistently high or low under quiet and disturbed geomagnetic conditions. The revealed longitudinal structure of the geomagnetic field variability in quiet geomagnetic conditions is caused by the discrepancy of the geographic and magnetic poles and by the spatial anomalies of different scales in the main magnetic field of the Earth. Variations of ionospheric parameters are shown to exhibit a pronounced longitudinal inhomogeneity with changing geomagnetic conditions. This inhomogeneity is associated with the longitudinal features of background and disturbed structure of the geomagnetic field. During the recovery phase of a storm, important role in dynamics of the mid-latitude ionosphere may belong to wave-like thermospheric disturbances of molecular gas, propagating westward for several days. Therefore, it is necessary to extend the time interval for studying the ionospheric effects of strong magnetic storms by a few days after the end of the magnetospheric source influence, while the disturbed regions in the thermosphere continues moving westward and causes the electron density decrease along the trajectories of propagation.     

Optical observations of ion tail structures and velocity flows in comet Halley during the period of the five spacecraft encounters

Observations of the distribution and evolution of a number of the major constituents of the neutral coma (CN, C₂, CH, O, H, Na) of Comet Halley were made during two observing periods, each of 3 weeks duration, from the Table Mountain Observatory, California. The first period was pre-perihelion, in late November/December 1985. The second period, from Feb 28 to March 22 1986, covered the five close spacecraft encounters with Halley, and when ICE flew some 20 M Km upstream of Halley. Sodium emission was recorded in early Dec 1985 from the near-nuclear region at a heliocentric distance of 1.4 AU, an observation confirmed with the UCL Doppler Imaging system. The CN coma could be detected to an outer diameter of more than 4M Km in Dec 1985, and 5 – 6M Km in early March 1986, allowing the production of heavy cometary pick-up ions to be estimated. Observations of the cometary ion coma (H₂O⁺ and CO⁺ ions) showed considerable variability from day to day, particularly during the period of the spacecraft encounters. These observations have been used, in conjuction with the neutral coma data, to map the flow field of cometary ions. In early Dec. 1985, Halley developed a traditional “type I” ion tail, which persisted until late April 1986. It has also been possible to evaluate the ion flow fields within the narrow core of the ion tail, and in the surrounding diffuse, low density, regions populated by pick-up and extracted cometary ions, and by slowed solar wind ions. Tail disconnection events were observed on several occasions, particularly between the VEGA 2 and GIOTTO encounters, and with a highly spectacular event on March 19 1986.