Dependence of geomagnetic activity during magnetic storms on the solar wind parameters for different types of streams: 3. Development of storm

This paper is a continuation of (Nikolaeva et al., 2011, 2012) and it analyzes the development of the main phase of 190 magnetic storms with Dst ≤ −50 nT depending on the type of source in the solar wind (magnetic clouds, MC; corotating interaction regions, CIR; Ejecta; Sheath before them, Sh_E; Sheath before MC, Sh_MC; all Sheath regions before ICME, Sh_E + Sh_MC; all ICME, MC + Ejecta; and an indeterminate type of solar wind stream, IND).     

Dependence of geomagnetic activity during magnetic storms on the solar wind parameters for different types of streams

The dependence of the maximal values of the |Dst| and AE geomagnetic indices observed during magnetic storms on the value of the interplanetary electric field (E _y ) was studied based on the catalog of the large-scale solar wind types created using the OMNI database for 1976–2000 [Yermolaev et al., 2009]. An analysis was performed for eight categories of magnetic storms caused by different types of solar wind streams: corotating interaction regions (CIR, 86 storms); magnetic clouds (MC, 43); Sheath before MCs (Sh_MC, 8); Ejecta (95); Sheath (Sh_E, 56); all ICME events (MC + Ejecta, 138); all compression regions Sheaths before MCs and Ejecta (Sh_MC + Sh_E, 64); and an indeterminate type of storm (IND, 75). It was shown that the |Dst| index value increases with increasing electric field E _y for all eight types of streams. When electric fields are strong (E _y > 11 mV m⁻¹), the |Dst| index value becomes saturated within magnetic clouds MCs and possibly within all ICMEs (MC + Ejecta). The AE index value during magnetic storms is independent of the electric field value E _y for almost all streams except magnetic clouds MCs and possibly the compressed (Sheath) region before them (Sh_MC). The AE index linearly increases within MC at small values of the electric field (E _y < 11 mV m⁻¹) and decrease when these fields are strong (E _y > 11 mV m⁻¹). Since the dynamic pressure (Pd) and IMF fluctuations (σB) correlate with the E _y value in all solar wind types, both geomagnetic indices (|Dst| and AE) do not show an additional dependence on Pd and IMF δB. The nonlinear relationship between the intensities of the |Dst| and AE indices and the electric field E _y component, observed within MCs and possibly all ICMEs during strong electric fields E _y , agrees with modeling the magnetospheric-ionospheric current system of zone 1 under the conditions of the polar cap potential saturation.     

Does the duration of the magnetic storm recovery phase depend on the development rate in its main phase? 2. A new method

In contrast to our previous work (Yermolaev et al., 2015), in which we used the magnetic storm recovery phase duration, the exponential time of the recovery phase of magnetic storms generated by three interplanetary driver types (CIR, Sheath, and ICME) is introduced in the present work. The dependence of these times on the storm development rate |Dst_min|/ΔT (where ΔT is the storm main phase duration) is studied. A similar physical result has been achieved despite the different data analysis method used: the times of the storm recovery and development rates correlate for storms induced by CIR and Sheath compression regions, and any relation between these parameters is absent for storms induced by ICME.     

Development of substorm bulges during storms of different interplanetary origins

Different solar wind structures are observed: magnetic clouds (MC), recurrent streams (RS), and regions of their interaction with undisturbed solar wind (Sheath and CIR). Three of these structures, Sheath, CIR, MC, are the sources of geomagnetic storms. We have searched for distinctions in the development of substorm bulges occurring during geomagnetic storms connected with the MC, Sheath and CIR. Solar wind parameters were taken from the Wind spacecraft and the auroral bulge parameters were obtained from the Ultra Violet Imager onboard Polar spacecraft. We determined the dimensions of the auroral bulges, the poleward aurora propagation, and the onset latitude of auroral bulges. It is shown that auroral bulges “geometry” is different for the examined types of storms. In consequence, the ratio between longitudinal and latitudinal sizes is also different.     

Prediction of maximal daily average values of relativistic electron fluxes in geostationary orbit during the magnetic storm recovery phase

The relation of the maximal daily average values of the relativistic electron fluxes with an energy higher than 2 MeV, obtained from the measurements on GOES geostationary satellites, during the recovery phase of magnetic storms to the solar wind parameters and magnetospheric activity indices has been considered. The parameters of Pc5 and Pi1 geomagnetic pulsations and the relativistic electron fluxes during the prestorm period and the main phase of magnetic storms have been used together with the traditional indices of geomagnetic activity (A _E, K _p, D _st). A simple model for predicting relativistic electron fluxes has been proposed for the first three days of the magnetic storm recovery phase. The predicted fluxes of the outer radiation belt relativistic electrons well correlate with the observed values (R ∼ 0.8–0.9).     

Erosion of the inner magnetosphere during geomagnetic storms

Using the empirical magnetic field model dependent on the Dst index and solar wind dynamic pressure, we calculated the behaviour of the contour B = B_s in the equatorial plane of the magnetosphere where B_s is the magnetic field in the subsolar point at the magnetopause. The inner domain of the magnetosphere outlined by this contour contains the bulk of geomag-netically trapped particles. During quiet time the boundary of the inner magnetosphere passes at the distance ∼10R_E at noon and at ∼7R_E at midnight. During very intense storms this distance can be reduced to 4–5 R_E for all MLT. The calculation results agree well with the satellite measurements of the magneto-pause location during storms. The ionospheric projection of the B = B_s contour calculated with the Euler potential technique is close to the equatorward edge of the auroral oval.     

磁暴主相多阶发展行星际起因的初步分析

本文利用1998～2006年与磁云有关的80起中强磁暴(Dst*≤-50 nT),对其主相期间不同发展阶数磁暴的行星际起因进行了统计分析.重点研究了鞘区磁场单独作用、磁云本体单独作用、鞘区与磁云共同作用以及其他复杂行星际结构在磁暴主相多阶发展中的相对重要性,并对导致磁暴主相增加一阶的行星际起因做了初步分析.统计结果表明...     

Magnetic storms and their effects in the lower ionosphere: Differences in storms of various types

During magnetic storms (MS’s) in the ionospheric D region, changes in the electron density and corresponding effects on radiowave propagation are observed. The differences in manifestations of MS’s in the lower ionosphere are mainly caused by the time and spatial differences in precipitations of energetic electrons. It is shown that the observed differences in the effects of storms in the D region are related to the differences in the corresponding types of MS’s determined by the observed fluxes of energetic electrons (E ∼ 0.1–2 MeV) at L ≈ 3–8. The storm types are identified by changes in the geomagnetic ap and AE indices and the ap/Dst and AE/Dst ratios during the recovery phase of a storm.     

The size of the auroral belt during magnetic storms

Using the auroral boundary index derived from DMSP electron precipitation data and the Dst index, changes in the size of the auroral belt during magnetic storms are studied. It is found that the equatorward boundary of the belt at midnight expands equatorward, reaching its lowest latitude about one hour before Dst peaks. This time lag depends very little on storm intensity. It is also shown that during magnetic storms, the energy of the ring current quantified with Dst increases in proportion to L_e^–3, where L_e is the L-value corresponding to the equatorward boundary of the auroral belt designated by the auroral boundary index. This means that the ring current energy is proportional to the ion energy obtained from the earthward shift of the plasma sheet under the conservation of the first adiabatic invariant. The ring current energy is also pronortional to E_mag, the total magnetic field energy contained in the spherical shell bounded by L_e and L_eq, where L_eq corresponds to the quiet-time location of the auroral precipitation boundary. The ratio of the ring current energy E_R to the dipole energy E_mag is typically 10%. The ring current leads to magnetosphere inflation as a result of an increase in the equivalent dipole moment.     

Variations in the amplitude of the chandler wobble

It is shown that within the framework of the Kolmogorov model the “energy” of the pole E(t) = x ₁² + x ₂² can be interpreted as a Markovian process. The exact analytical expression has been obtained for the density of the conditional probability of the quantity E(t) and the problem of the first passage time of the process E(t) has been analyzed. It was shown that the available data on the swing of the function E(t) are not at variance with the Kolmogorov model and a short-period drop of the amplitude of the Chandler wobble in the early 20th century fits this model at Q = 50–200 too; values of Q > 350 are less reasonable.     

A superposed epoch analysis of geomagnetic storms

A superposed epoch analysis of geomagnetic storms has been undertaken. The storms are categorised via their intensity (as defined by the Dst index). Storms have also been classified here as either storm sudden commencements (SSCs) or storm gradual commencements (SGCs, that is all storms which did not begin with a sudden commencement). The prevailing solar wind conditions defined by the parameters solar wind speed (v_sw), density (_sw) and pressure (P_sw) and the total field and the components of the interplanetary magnetic field (IMF) during the storms in each category have been investigated by a superposed epoch analysis. The southward component of the IMF, appears to be the controlling parameter for the generation of small SGCs (-100 nT< minimum Dst\leq-50 nT for\geq4 h), but for SSCs of the same intensity solar wind pressure is dominant. However, for large SSCs (minimum Dst\leq-100 nT for \geq4 h) the solar wind speed is the controlling parameter. It is also demonstrated that for larger storms magnetic activity is not solely driven by the accumulation of substorm activity, but substantial energy is directly input via the dayside. Furthermore, there is evidence that SSCs are caused by the passage of a coronal mass ejection, whereas SGCs result from the passage of a high speed/ slow speed coronal stream interface. Storms are also grouped by the sign of B_z during the first hour epoch after the onset. The sign of B_z at t=+1 h is the dominant sign of the B_z for \sim24 h before the onset. The total energy released during storms for which B_z was initially positive is, however, of the same order as for storms where B_z was initially negative.     

Variations in the ULF index of daytime geomagnetic pulsations during recurrent magnetic storms

A new index of wave activity (ULF index) is applied to analyze daytime magnetic pulsations in the Pc5 range (f = 2–7 mHz) during ten successive recurrent magnetic storms (CIR (corotating interaction region) storms) of 2006. The most intense daytime geomagnetic Pc5 pulsations on the Earth’s surface in all phases of CIR storms are predominantly observed in the pre-noon sector at latitudes higher than 70°, while those in CME storms (storms initiated by coronal mass ejection (CME)) are observed at latitudes lower than 70°. A comparison of wave activity during CIR and CME storms has shown that the amplitude of Pc5 pulsations in CIR storms is much smaller than that in CME storms and the spectrum maximum is observed at lower frequencies and higher latitudes. At the same time, the mechanism of ULF wave generation during both types of magnetic storms seems to be similar, namely, resonance of magnetic field lines due to the development of the Kelvin-Helmholtz instability caused by an approach of a high-velocity solar wind stream to the Earth’s magnetosphere. Since resonance oscillations are excited only in the closed magnetosphere, the higher-latitude position of the Pc5 pulsation intensity maximum in CIR storms points to larger dimensions of the daytime magnetosphere during CIR storms as compared to CME storms.     

南北地震带两次典型强震周边自然电场变化及机理探讨

根据场源不同可将地电场E分为大地电场E_T和自然电场E_SP。空间S_q电流系和固体潮一般被认为是大地电场E_T的起源;自然电场E_SP源于地下介质的物理化学作用,其局部性变化相对稳定。2013年南北地震带相继发生了芦山M_S7.0地震（30.3°N,103.0°E）和岷县、漳县M_S6.6地震（34.5°N,104.2°E）。两次地震周边约400 km范围内的地电场台站（成都、汉王）数据表明震前自然电场E_SP出现明显的小幅度突跳或大幅度跃变,同时在时间上具有准同步性,其他地电场台站也有类似变化。基于大地电场岩体裂隙水（电荷）渗流（移动）模型,对两次地震周边自然电场E_SP的变异机理进行探讨,尝试解析其物理过程。     

Properties of electric turbulence in the polar cap ionosphere

Small-scale (scales of ∼0.5–256 km) electric fields in the polar cap ionosphere are studied on the basis of measurements of the Dynamics Explorer 2 (DE-2) low-altitude satellite with a polar orbit. Nineteen DE-2 passes through the high-latitude ionosphere from the morning side to the evening side are considered when the IMF z component was southward. A rather extensive polar cap, which could be identified using the ɛ-t spectrograms of precipitating particles with auroral energies, was formed during the analyzed events. It is shown that the logarithmic diagrams (LDs), constructed using the discrete wavelet transform of electric fields in the polar cap, are power law (μ ∼ s ^α). Here, μ is the variance of the detail coefficients of the signal discrete wavelet transform, s is the wavelet scale, and index α characterizes the LD slope. The probability density functions P(δE, s) of the electric field fluctuations δE observed on different scales s are non-Gaussian and have intensified wings. When the probability density functions are renormalized, that is constructed of δE/s ^γ, where γ is the scaling exponent, they lie near a single curve, which indicates that the studied fields are statistically self-similar. In spite of the fact that the amplitude of electric fluctuations in the polar cap is much smaller than in the auroral zone, the quantitative characteristics of field scaling in the two regions are similar. Two possible causes of the observed turbulent structure of the electric field in the polar cap are considered: (1) the structure is transferred from the solar wind, which is known to have turbulent properties, and (2) the structure is generated by convection velocity shears in the region of open magnetic field lines. The detected dependence of the characteristic distribution of turbulent electric fields over the polar cap region on IMF B _y and the correlation of the rms amplitudes of δE fluctuations with IMF B _z and the solar wind transfer function (B _y ² + B _z ²)^1/2sin(θ/2), where θ is the angle between the geomagnetic field and IMF reconnecting on the dayside magnetopause when IMF B _z < 0, together with the absence of dependence on the IMF variability are arguments for the second mechanism.     

The Polar Cap (PC) indices: Relations to solar wind parameters and global magnetic activity

The polar geomagnetic activity resulting from solar wind–magnetosphere interactions can be characterized the Polar Cap (PC) indices, PCN and PCS. PC index values are derived from polar magnetic variations calibrated on a statistical basis such that the index approximate values in units of mV/m of the interplanetary “geo-effective” (or “merging”) electric field (E_M) conveyed by the solar wind. The timing and amplitude relations of the PC index to solar wind plasma and magnetic field parameters are reported. The solar wind effects are parameterized in terms of the geo-effective electric field (E_M) and the dynamical pressure (P_DYN). The PC index has a delayed and damped response to E_M variations and display saturation-like effects for E_M values exceeding 10 mV/m. Steady or slowly varying levels of solar wind dynamical pressure have little or no impact on the PC index above the effects related to E_M for which the solar wind velocity is also a factor. Sharp increases in the dynamical pressure generate impulsive variations in the PC index comprising a initial negative impulse of 5–10 min duration followed by a positive impulse lasting 10–20 min. Typical amplitudes of both the negative and the positive impulses are 0.2–0.5 units. A sharp decrease in the pressure produces the inverse sequence of pulses in the PC index. Auroral substorm activity represented by the AL index level has a marked influence on the average PC/E_M level at the transition from very quiet (AL0 nT) to disturbed conditions while more or less disturbed conditions (AL<−100 nT) have no systematic effect on the average PC/E_M values. At distinct substorm events the PC/E_M ratio has a minimum (0.8) in the pre-onset phase at around 20 min before substorm onset. The average ratio gradually increases in the expansion phase to reach a maximum value (1.1) at around 40 min after substorm onset (or 20 min after the largest (negative) peak in AL). At substorm recovery during the next 2 h the PC/E_M ratio decreases. Finally, we report on the application of polar magnetic variations to model the disturbance storm time (Dst) index development during magnetic storms by using the PC index as a source function to quantify the energy input to the ring current representing accumulated storm energy and characterized by the Dst index.     

Experimental evidence for direct penetration of ULF waves from the solar wind and their possible effect on acceleration of radiation belt electrons

The event of March 12–19, 2009, when a moderately high-speed solar wind stream flew around the Earth’s magnetosphere and carried millihertz ultralow-frequency (ULF) waves, has been analyzed. The stream caused a weak magnetic storm (D _{st min} = −28 nT). Since March 13, fluxes of energetic (up to relativistic) electrons started increasing in the magnetosphere. Comparison of the spectra of ULF oscillations observed in the solar wind and magnetosphere and on the Earth’s surface indicated that a stable common spectral peak was present at frequencies of 3–4 mHz. This fact is interpreted as evidence that waves penetrated directly from the solar wind into the magnetosphere. Possible scenarios describing the participation of oscillations in the acceleration of medium-energy (E > 0.6 MeV) and high-energy (E > 2.0 MeV) electrons in the radiation belt are discussed. Based on comparing the event with the moderate magnetic storm of January 21–22, 2005, we concluded that favorable conditions for analyzing the interaction between the solar wind and the magnetosphere are formed during a deep minimum of solar activity.     

Seasonal variations in the occurrence of geomagnetic storms

Seasonal variations in the onset of magnetic storms are investigated. For the purposes of this study storms have been defined as events in which Dst falls below -50 nT for at least four consecutive hours. The storms have been classified as either storm sudden commencements (SSCs; storms initiated by a sudden commencement) or as storm gradual commencements (SGCs; all other storms). It is found that the semi-annual variation of magnetic activity is reflected in the occurrence statistics of SGC events only, indicative that the solar wind origin is different for SSCs and SGCs. It is suggested that the heliospheric latitude model of seasonal magnetic activity is relatively ineffective in modulating the previously observed seasonal variations in the occurrence of magnetic storms.     

Specific features of development of two-dimensional convection in a rotating compressible conducting fluid

An axisymmetric model of convection in a rotating cylinder in an external uniform magnetic field has been considered. In the considered model, the meridional circulation is created by a nonuniform rotation of the lower boundary relative to the other boundaries. In the considered model, the time of formation of the stationary regime in the magnetic field considerably increases if the vertical density (compressibility) inhomogeneity is taken into account for Ekman numbers of E = E _M = 3 × 10⁻³. This example shows that the compressibility of a medium should be taken into account in the convection and dynamics of the magnetic field when the magnetohydrodynamics of the Earth is analyzed.