Regularities in changes of the daytime and nighttime electric fields in the ionosphere during a substorm

Dynamics of changes of the electric fields in the daytime and nighttime sectors of the high-latitude ionosphere is studied. The following facts are demonstrated: (a) With the onset of active phases of substorms, a decrease occurs in the value of the electric potential difference in the daytime sector. It is especially abrupt with the onset of the expansion phase, which is related to a decrease and the end of reconnection at the dayside magnetopause. (b) With the onset of active phases of substorms, an increase in the value of the electric potential difference in the nighttime sector occurs. It is sharply intensified with the onset of the expansion phase, this fact being caused by the intensification of the reconnection in the nighttime magnetotail. (c) There exists a relation, close to a functional dependence, between the values of the electric potential differences in the nighttime sector and | AL | index.     

A Magnetospheric Substorm and the Motion of the Magnetopause

Using a single event as an example, we make an analysis of the time development of a substorm and estimate its influence on the motion of the low-latitude boundary of the magnetosphere. To this end, we compare the data on plasma and magnetic field obtained by five spacecraft (WIND, INTERBALL-1, GEOTAIL, GOES-8, and GOES-9) with measurements made by ground-based stations. It is shown that the release of energy of the geomagnetic tail begins from a disruption of the current sheet near the Earth. The high-speed plasma stream that transfers a magnetic flux to the Earth and can have an effect on the magnetic field configuration near the Earth is detected later. Almost simultaneously with a substorm onset a series of magnetopause crossings has been detected by the INTERBALL-1 satellite on the evening side of the low-latitude magnetosphere. In this paper we consider some of possible causes of this motion of the magnetosphere boundary, including variations of parameters of the solar wind, Kelvin-Helmholtz instability, and substorm processes. It is shown that fast motions of the magnetopause are detected almost simultaneously with field variations in the near magnetotail of the Earth and geomagnetic pulsations Pi2 on ground-based stations. A sufficiently high degree of correlation (K = 0.67) between the amplitude of Pi2 pulsations and the amplitude of magnetic field variations near the magnetopause is probably indicative of the connection of short-term motions of the magnetosphere boundary with the tail current disruption and the process of formation of a substorm current wedge.__________Translated from Kosmicheskie Issledovaniya, Vol. 43, No. 4, 2005, pp. 248–259.Original Russian Text Copyright © 2005 by Nikolaeva, Parkhomov, Borodkova, Klimov, Nozdrachev, Romanov, Yermolaev.     

Some features of solar cosmic ray penetration into the Earth’s magnetosphere

We compared fluxes of the 1–100 MeV solar energetic particles (SEP) measured in the interplanetary medium (ACE) and in the magnetosphere (Universitetsky-Tatiana, POES—in polar caps, and GOES-11—at geosynchronous orbit) during several SEP events of 2005–2006. Peak intensities of the SEP fluxes inside and outside the magnetosphere were compared for each event. It is shown that observed inside-outside difference depends mainly on direction of interplanetary magnetic field (IMF), on degree of the SEP anisotropy (pitch-angle distribution) in IMF, and on distance of the dayside magnetopause from the Earth.     

The Fields of Magnetospheric Current Systems on the Earth's Surface in an Interval of the Magnetic Storm Studied under the International Program on Solar–Terrestrial Physics

Using the magnetic storm in January 1997 as an example, we examined the possibilities to employ the magnetospheric field T96 [1, 2] and the dynamic paraboloid model PM of the magnetosphere [3] for modeling the D _st variation. We have revealed the necessity to refine the results of normalizing the free parameters of the model T96 according to the solar wind parameters. The contributions to the D _st variation of magnetic fields of basic large-scale magnetospheric current systems (the field DCF on the magnetopause, the field DR of the ring current, and the field DT in the magnetotail) are estimated for different phases of the storm from model calculations. Possible causes of a discrepancy between the results of modeling D _st using the T96 and PM models are discussed. Special emphasis is made on the ratios of contributions into the D _st variation of the fields of the magnetotail and the ring current in the main phase of magnetic storms and on the contributions to D _st of the fields of various current systems at the recovery phase.     

Geosynchronous magnetopause crossings on October 29–31, 2003

During the period October 29–31, 2003, geosynchronous magnetopause crossings (GMC) have been identified based on the magnetic data of the GOES series spacecraft and plasma data of the LANL series spacecraft. It is shown that most of the time the size of the dayside magnetosphere was highly decreased under the effect of very high pressure associated with high velocities and densities of the solar wind plasma, as well as high negative values of the B_z component of the interplanetary magnetic field (IMF). For tens of hours the subsolar magnetopause was deep inside the geosynchronous orbit. During the main phase and at the maximum of the strong geomagnetic storms that occurred in the period under consideration, the dayside magnetosphere was characterized by a strong dawn-dusk asymmetry, so that its size in the postnoon sector considerably exceeded the size in the pre-noon sector. The geomagnetic disturbances in the morning on October 30 and 31, 2003 were accompanied by global magnetospheric pulsations with periods of 5–10 min and high amplitude (up to 0.8 R_E).Translated from Kosmicheskie Issledovaniya, Vol. 42, No. 6, 2004, pp. 574–584.Original Russian Text Copyright © 2004 by Dmitriev, Suvorova.     

Verification of the Accuracy of the Magnetopause Empirical Models Using the INTERBALL-1 Satellite Data

A statistical analysis of the shape and location of the magnetopause according to the INTERBALL-1 satellite data for the period 1995–1997 is carried out. The instants of crossing the magnetosphere boundaries obtained by the plasma and magnetic data are compared with computations based on three empirical models, namely, Petrinec and Russel, 1996; Shue et al., 1997; and Shue et al., 1998. The state of the interplanetary medium (dynamic pressure of the solar wind plasma P _d and the B _z component of the interplanetary magnetic field) was determined by the measurements onboard the WIND spacecraft. We estimate the accuracy of the considered models for different groups of boundary crossings: single, multiple with small duration (less than 40 min), and multiple with large duration (more than 40 min). It is demonstrated that the small-scale motions of the boundary (<1R _E) are observed more often in the dayside magnetosphere, especially near the cusp region. Large-scale boundary oscillations (>1R _E) are more common in the tail region of the magnetosphere, namely, its flanks. Various models give similar results: about 50% of all events have deviations by more than 1R _E from the model locations. In some cases, the deviation of the measured location of the magnetosphere boundary from the model prediction may be as large as 5–6R _E for all three models considered, the actual boundary being more often located nearer to the Earth than the result of model computations. The best model is that of Shue et al., 1998, but it does not differ significantly from the other models.     

Test-Particle Trajectories in a “Sheared” Stationary Field: Newton–Lorentz and First Order Drift Numerical Simulations1

We study a magnetic field distribution that is nonuniform and sheared like in tangential discontinuities. This distribution is an input parameter for the numerical integration of the equations of motion of the test-particle and of its guiding center. Two different electric field distributions are alternatively tested. In the first case, the electric field is uniform and constant like the electric field prescribed in the large-scale, steady-state reconnection models. The numerical solution shows that in this case the test-particle is trapped within the discontinuity into a region where (i) B goes to zero or (ii) the magnetic vector becomes exactly parallel to the electric field. In the second case, we consider an electric field, which is nonuniform. Its components are computed such that the zero order (or electric) drift is everywhere perpendicular to the discontinuity surface and its value is conserved throughout the simulation. In this case the numerically integrated trajectory of the test-particle penetrates the discontinuity for any angle of shear of B. Direct comparison between exact (Newton–Lorentz) and approximated (first order drift) numerical solutions shows that the mathematical singularities of the latter do not correspond to any physical singularity of the exact equation of motion of the particle.     

Why does the total pressure on the subsolar magnetopause differ from the solar wind dynamic pressure?

Based on analysis of MHD equations and the results of numerical simulation in the magneto-sheath it is demonstrated that the total pressure on the magnetopause differs from the solar wind dynamic pressure in the majority of cases. From the equation of motion it follows that the total pressure is reduced due to deflection from the Sun-Earth line. At the same time, it increases because of formation of a magnetic barrier. This result is consistent with experimentally observed expansion of the magnetosphere for the radial direction of the interplanetary magnetic field, when no magnetic barrier is formed. In this paper we compare the behavior of pressure along the Sun-Earth line for the northward and radial interplanetary field, using the results of numerical MHD simulation and observational data from THEMIS. In the isotropic MHD approximation, the difference between the total pressure on the subsolar magnetopause at northern and radial IMFs does not exceed 10–12 percent. However, in the anisotropic approximation this difference increases up to 15–20 percent. The results of anisotropic modeling well agree with observed averaged profiles of pressure components in the subsolar magnetosheath.     

Turbulent fluctuations of plasma and magnetic field parameters in the magnetosheath and the low-latitude boundary layer formation: Multisatellite observations on March 2, 1996

The results of simultaneous analysis of plasma and magnetic field characteristics measured on the INTERBALL/Tail Probe, WIND and Geotail satellites on March 2, 1996, are presented. During these observations the INTERBALL/Tail Probe crossed the low-latitude boundary layer, and the WIND and Geotail satellites measured the solar wind’s and magnetosheath’s parameters, respectively. The plasma and magnetic field characteristics in these regions have been compared. The data of the Corall, Electron, and MIF instruments on the INTERBALL/Tail Probe satellite are analyzed. Fluctuations of the magnetic field components and plasma velocity in the solar wind and magnetosheath, measured onboard the WIND and Geotail satellites, are compared. The causes resulting in appearance of plasma jet flows in the low-latitude boundary layer are analyzed. The amplitude of magnetic field fluctuations in the magnetosheath for a studied magnetosphere boundary crossing is shown to exceed the magnetic field value below the magnetopause near the cusp. The possibility of local violation of pressure balance on the magnetopause is discussed, as well as penetration of magnetosheath plasma into the magnetosphere, as a result of magnetic field and plasma flux fluctuations in the magnetosheath.     

Dust Plasma Dynamics in the Geomagnetic Tail

The results of numerical modeling of dust plasma dynamics in the geomagnetic tail are presented for the southward and northward orientations of the interplanetary magnetic field in the model with reconnection of magnetic field lines of geomagnetic and interplanetary magnetic fields. It is shown that in reconnection regions the compression shock waves arise, at the fronts of which the dust plasma fluxes with higher-than-background density are generated.     

Relationship between ionospheric electric fields and magnetic activity indices

The relations between electric fields in the daytime and nighttime sectors of the polar ionosphere and magnetic activity indices of auroral region (AL) and northern polar cap (PCN) are studied. It is found that the above relations do exist and are described by: a) equations U _pc⁽¹⁾ (kV) = 27.62 + 21.43PCN with a correlation coefficient R = 0.87 and U _pc⁽¹⁾ (kV) = 4.06 + 49.21PCN - 6.24 PCN² between the difference in the electric potentials across the polar cap in the daytime sector U _pc⁽¹⁾ and PCN and b) regression equation U _pc⁽²⁾ (kV) = 23.33 + 0.08|AL| with R = 0.86 between the difference in the electric potentials across the polar cap in the nighttime sector U _pc⁽²⁾ and |AL|. It is shown that: a) it is possible to use the AL and PCN indices for real-time diagnostics of instantaneous values of the electric fields in the daytime and nighttime sectors of the polar ionosphere in the process of a substorm development; b) at the expansion phase of a substorm, due to calibration of PCN values by the values of the solar wind electric field E _sw, the PCN index does not feel the contribution of the western electrojet and, accordingly, the contribution of the nighttime ionospheric electric field U _pc⁽²⁾, governed by the reconnection in the magnetospheric tail.     

Activation of the tail open part during the magnetospheric storm

In each polar cap (PC) we mark out “old PC” observed during quiet time before the event under consideration, and “new PC” that emerges during the substorm framing the old one and expanding the PC total area. Old and new PCs are the areas for the magnetosphere old and new tail lobes, respectively. The new lobe variable magnetic flux Ψ₁ is usually assumed to be active, i.e. it provides the electromagnetic energy flux (Poynting flux) ɛ′ transport from solar wind (SW) into the magnetosphere. The old lobe magnetic flux Ψ₂ is supposed to be passive, i.e. it remains constant during the disturbance and does not participate in the transporting process which would mean the old PC electric field absolute screening from the convection electric field created by the magnetopause reconnection. In fact, screening is observed, but far from absolute. We suggest a model of screening and determine its quantitative characteristics in the selected superstorm. The coefficient of a screening is the β = Ψ₂/Ψ₀₂, where Ψ₀₂ = const is open magnetic flux through the old PC measured prior to the substorm, and Ψ₂ is variable magnetic flux through the same area measured during the substorm. We consider three various regimes of disturbance. In each, the coefficient β decreased during the loading phase and increased at the unloading phase, but the rates and amplitudes of variations exhibited a strong dependence on the regime. We interpreted decrease in β as a result of involving the old PC magnetic flux Ψ₂, which was considered to be constant earlier, in the Poynting flux ɛ′ transport process from solar wind into the magnetosphere. Transport process weakening at the subsequent unloading phase creates increase in β. Estimates showed that coefficient β during each regime and the computed Poynting flux ɛ′ varied manifolds. In general, unlike the existing substorm conception, the new scenario describes an unknown earlier of tail lobe activation process during a substorm growth phase that effectively increases the accumulated tail energy for the expansion and recovery phases.     

Possible Origin of Ultralow-Frequency Whistlers

The sonographic analysis of records of ultralow-frequency emissions recorded by the induction magnetometer at the Mirnyi observatory (Antarctica) in 1981–1985 has revealed the presence of a special class of signals in the frequency band 0.25–5 Hz having a characteristic dispersion reminiscent of the dynamic spectra of LHR-whistlers observed in the VLF band. The ULF whistlers are observed, as a rule, at morning and evening hours of local time at moderate values of the K _p-index (0 < K _p < 4). The analysis of the frequency–time behavior of observed signals shows that the canalized propagation of short broadband pulses as magnetosonic waves in the layered medium, for example, in the magnetospheric current sheet can serve as a possible cause of the dispersion. Such sporadic phenomena on the magnetopause as microreconnections or FTE-events may be sources of these pulses.     

Pressure balance on the magnetopause near the subsolar point according to observational data of the THEMIS project satellites

An analysis of the pressure balance on the magnetopause near the subsolar point has been made for 18 crossings of the magnetopause by the THEMIS project satellites under magneto-quiet conditions. Dynamic and static pressures of plasma are determined, as well as magnetic pressure in the magnetosheath, and magnetic and plasma static pressure inside the magnetosphere. Variations of the total pressure have been studied in the case when one satellite is located inside the magnetosphere and another one stays in the magnetosheath near the magnetopause. It is demonstrated that for 18 investigated events the condition of pressure balance at the subsolar point is valid on average with an accuracy of 7%, within measurement errors and under applicability of the approximation of anisotropic magnetic hydrodynamics to collisionless plasma of the magnetosheath and magnetosphere.     

Modeling the Modification of the Nighttime High-Latitude F-Region by Powerful HF Radiowaves

A mathematical model of the high-latitude F-region, taking the ionospheric plasma convection into account, is used for modeling a response of the auroral F-region to irradiation by powerful high-frequency (HF) radiowaves. The model enables us to calculate the time variations in the profiles of the electron concentration, the velocity of positive ions, and the ion and electron temperatures in the part of a magnetic-flux tube moving over an artificial ionospheric heater under the action of a convective electric field. The modeling was carried out for a part of the magnetic flux tube intersecting the F-layer irradiated by the Norwegian heater at Tromsö when it is located near the midnight magnetic meridian. The calculations were made for the equinox conditions under the high solar and low geomagnetic activity. The results of our modeling show that substantial variations in the profiles of the electron temperature, the velocity of positive ions, and the electron concentration can be produced in the nighttime high-latitude F-layer due to HF heating. The perturbation caused by a 20 s rectangular pulse should exist for about 20 min at the level of the F-layer maximum. The disturbed plasma volume can leave the region irradiated by the heater and move away for over 500 km within the above-mentioned time period.     

Structure of the earth’s magnetosheath at small fluctuations in the interplanetary magnetic field direction

Spatial structure of the magnetosheath of the Earth was studied under the conditions when no sharp (more than 40° during 5 min) changes in the interplanetary magnetic field direction were observed. On the basis of 24 flights of the Interball-1 satellite through the magnetosheath, it is found that three regions differing from each other by parameters of the field and plasma can be observed in the magnetosheath under the above-indicated conditions. These regions also differ from the solar wind region before front of the Earth’s magnetospheric bow shock. Empirical distributions of parameters were studied in each region. Taking into account the influence of the interplanetary magnetic field direction on the processes in the magnetosheath, the cases of quasi-perpendicular and quasi-parallel shock waves were considered separately. The study showed that the distribution of parameters in the selected regions (in the solar wind before front of the bow shock, in the magnetosheath behind the bow shock (post-shock), in the region of the magnetosheath with minimal fluctuations in the field, and in the inner magnetosheath) differ from each other at any interplanetary magnetic field direction.     

Low-energy ion fluxes and magnetic field on the inner boundary of the heliosphere

A comparison of temporal profiles of low-energy ion intensity and magnetic field magnitude in different periods of solar activity in the outer heliosphere is carried out using the data of the Voyager 1 and Voyager 2 spacecraft. It is shown that temporal, spectral, and statistical characteristics of particle fluxes and magnetic field in the heliospheric regions before and after the terminal shock in 2002–2008 had similar dynamics in different hemispheres. This similarity allowed one to assume that, in the region of the inner heliospheric boundary, a quasistable spatial structure existed moving together with the terminal shock in accordance with the solar wind pressure, as well as, probably, under the action of the interstellar medium. It was revealed that the spatial dimensions of most details of this structure are less on Voyager 2, which, probably, is due to variation of the solar activity level, difference in latitude of spacecraft disposition, and also the influence of the interstellar magnetic field.     

Calculation of the Initial Magnetic Field for Mercury’s Magnetosphere Hybrid Model

Several types of numerical models are used to analyze the interactions of the solar wind flow with Mercury’s magnetosphere, including kinetic models that determine magnetic and electric fields based on the spatial distribution of charges and currents, magnetohydrodynamic models that describe plasma as a conductive liquid, and hybrid models that describe ions kinetically in collisionless mode and represent electrons as a massless neutralizing liquid. The structure of resulting solutions is determined not only by the chosen set of equations that govern the behavior of plasma, but also by the initial and boundary conditions; i.e., their effects are not limited to the amount of computational work required to achieve a quasi-stationary solution. In this work, we have proposed using the magnetic field computed by the paraboloid model of Mercury’s magnetosphere as the initial condition for subsequent hybrid modeling. The results of the model have been compared to measurements performed by the Messenger spacecraft during a single crossing of the magnetosheath and the magnetosphere. The selected orbit lies in the terminator plane, which allows us to observe two crossings of the bow shock and the magnetopause. In our calculations, we have defined the initial parameters of the global magnetospheric current systems in a way that allows us to minimize paraboloid magnetic field deviation along the trajectory of the Messenger from the experimental data. We have shown that the optimal initial field parameters include setting the penetration of a partial interplanetary magnetic field into the magnetosphere with a penetration coefficient of 0.2.     

Formation of a Polarization Jet during the Expansion Phase of a Substorm: Results of Ground-Based Measurements

The ground-based polarization jet measurements at the Yakutsk ionosphere station (L= 3.0) for the years 1989–1991 (110 events) are compared with variations of the AE-index and with parameters of the local magnetic activity. It is shown that polarization jet development in the near midnight sector can occur during a period of no longer than 10 min on the expansion phase of a substorm. The formation of the polarization jet is accompanied by a specific magnetic field variation corresponding in shape to a fast passage of the Harang Discontinuity above the station. Statistical data are given on ground level observations of the polarization jet, which are close to those measured from satellites. The mean delay (averaged over the full data bank) between the onset of a substorm with AE 500 nT and the moment of the polarization jet appearance at L= 3.0 is equal to 0.5 h near midnight and to 1.0 – 1.5 h in the evening sector. Estimations show that the duration of the polarization jet formation when energetic ions are injected into the Harang Discontinuity region above the ground station can last for about 10 min, and during this time the Harang Discontinuity can be shifted to the west. This is in qualitative agreement with the described observations.     

Solar wind inhomogeneity front inclination effect on properties of front-caused long-period geomagnetic pulsations

On the basis of data of two networks of Canadian stations and also of extra- and intra-magnetospheric satellites, daytime long-period geomagnetic pulsations related to sudden impulses of the dynamic pressure of the solar wind (SW) are studied. The influence of SW parameters, interplanetary magnetic field (IMF), and geomagnetic activity on the propagation direction, polarization, and amplitude of pulsations is discussed. It is shown that at arrival front of the solar wind inhomogeneity at the place of its tangency, surface oscillations within the range of Pc5 geomagnetic pulsations are excited on the magnetopause, and they run away from the tangency point to the nighttime side with increasing amplitude and opposite polarization. The pulsation properties and the position of the running-away point are explained by the mechanism of their excitation on the magnetopause by the inclined front of the inhomogeneity and also by the Kelvin-Helmholtz instability. Increases in SW density observed ahead of the shock front were able to cause pulsation excitation onsets prior to the sudden storms commencement (SSC) front arrival. The observed increase in geomagnetic activity after SSC could change the direction of pulsation propagation from anti-sunward to sunward. The analysis of oscillation spectra made it possible to assume that pulsations with a frequency of the order of 2.5 mHz are of a global character, they are not related to oscillations in SW and are excited by sharp SSC fronts.