Geosynchronous magnetic field response to the large and fast solar wind dynamic pressure change

We present a comparison of large and sharp solar wind dynamic pressure changes, observed by several spacecraft, with fast disturbances in the magnetospheric magnetic field measured by the GOES-8, 9 and 10 geosynchronous satellites. Almost 400 solar wind pressure changes in the period 1996–2003 were selected for this study. Using the large statistics we confirmed that increases (decreases) in the dynamic pressure always results in increases (decreases) in the magnitude of geosynchronous B_z component. The amplitude of the geosynchronous B_z response strongly depends on the location of observer relative to the noon meridian, from the value of solar wind pressure before the disturbance arriving and firstly – from the amplitude of the pressure change.     

Comparison of the magnetic field before the subsolar magnetopause with the magnetic field in the solar wind before the bow shock

Crossings of the magnetopause near the subsolar point are analyzed using data of THEMIS mission. Variations of the magnetic field near magnetopause measured by one of THEMIS satellites are studied and compared with simultaneous measurements in the solar wind by another THEMIS satellite. The time delay of the solar wind arrival at the subsolar point of the magnetopause is taken into account. 30 and 90 s averaging of the magnetic field in the magnetosheath is produced. The results of averaging are compared with the results of measurements in the solar wind before the bow shock and foreshock. It is shown, that _Bx$B_x$ component of the magnetic field near magnetopause is near to zero, which supports the possibility to consider the magnetopause as the tangential discontinuity. Comparatively good correlation of _By$B_y$ component in the solar wind and near the magnetopause is observed. The correlation of _Bz$B_z$ component near the magnetopause and IMF is practically absent, the sign of the _Bz$B_z$ near the subsolar point does not coincide with the sign of IMF _Bz$B_z$ in ∼$\sim$30% cases.     

Statistical properties of the IMF clock angle in the solar wind with northward and southward interplanetary magnetic field based on ACE observation from 1998 to 2009: Dependence on the temporal scale of the solar wind

Utilizing ACE satellite observations from 1998 to 2009, we performed the elaborate study on the properties of the clock angle θ_CA (arctan(By/Bz) (?90° to 90°) of the interplanetary magnetic field (IMF) in the solar wind at 1?AU. The solar wind with northward IMF (NW-IMF) and southward IMF (SW-IMF) are analyzed, independently. Statistical analysis shows that the solar wind with SW-IMF and NW-IMF has similar properties in general, including their durations, the IMF Bz and By components, and the IMF θ_CA. Then, the solar wind with NW-IMF (SW-IMF) is classified into five different temporal scales according to the duration of the NW-IMF (SW-IMF), i.e., very-short wind of 10–30?min, short-scale wind of 0.5–1?h, moderate-scale wind of 1–3?h, long-scale wind of 3–5?h, and super-long wind >5?h. Our analysis reveals that the IMF θ_CA has a distinct decrease with increase of the temporal scale of the solar wind. Next, the solar wind is classified into two groups, i.e., the high-speed solar wind (>450?km/s) and the low-speed solar wind (<450?km/s). Our analysis indicates that the IMF θ_CA depends highly on the solar wind speed. Statistically, high-speed solar wind tends to have larger IMF θ_CA than low-speed solar wind. The evolutions of the solar wind and IMF with the solar activity are further studied, revealing no clear solar variation of the IMF θ_CA. Finally, we analyze the monthly variation of the IMF θ_CA. Superposed epoch result strongly suggests the seasonal variation of the IMF θ_CA.     

Hybrid simulations of the O ion escape from Venus: Influence of the solar wind density and the IMF x component

As an initial effort to study the evolution of the Venus atmosphere, the influence of the solar wind density and the interplanetary magnetic field (IMF) x component (the x-axis points from Venus towards the Sun) on the O⁺ ion escape rate from Venus is investigated using a three-dimensional quasi-neutral hybrid (HYB-Venus) model. The HYB-Venus model is first applied to a case of the high-density (100 cm⁻³) solar wind interaction with Venus selected from the Pioneer Venus Orbiter observations to demonstrate its capability for the study. Two sets of simulations with a wide range of solar wind densities and different IMF x components are then performed. It is found that the O⁺ ion escape rate increases with increasing solar wind density. The O⁺ ion escape rate saturates when the solar wind density becomes high (above 100 cm⁻³). The results also suggest that the IMF x component enhances the O⁺ ion escape rate, given a fixed IMF component perpendicular to the x-axis. Finally, the results imply a higher ion loss rate for early-Venus, when solar conditions were dramatically different.     

On the possibility to check the magnetosphere’s model by CR: The strong geomagnetic storm in November 2003

We report on a study of cosmic ray cutoff rigidity variations during the strong geomagnetic storm of 18–24 November 2003. We employed the most recent Tsyganenko magnetospheric model to represent the very strong disturbed Magnetosphere. We used this magnetic field for the cosmic ray trajectory calculations to determine the geomagnetic cutoff rigidity throughout this period of severe geomagnetic disturbance. We determine the cutoff rigidity changes during this period by two methods, by trajectory calculations and by the spectrographic global survey method. The values of the change of cutoff rigidities obtained by two different methods are correlated with the D_st and interplanetary magnetic field and plasma parameters and result in correlation coefficients in the range 0.63–0.84 for the various cosmic ray stations. The result of this study indicates that the most significant contributions to the cutoff rigidity changes are due to D_st variation although the influence of solar wind density and B_z and B_y components of IMF variations is significant.     

Symmetric or asymmetric energy transfer from Interplanetary Coronal Mass Ejections to the magnetosphere depending on the solar dipole

The annual behaviour of monthly number of hours spent by the Earth in domains of either positive or negative B_y component of the interplanetary magnetic field (IMF) was studied. We used the hourly OMNI data in the cases of K_p > 3. The study was confined to the ascending phases of the four recent sunspot cycles when Interplanetary Coronal Mass Ejections (ICMEs) dominate among the sources of geoeffectiveness. Definite differences were found between the annual variations of the hourly sums. When the solar dipole is opposite to the terrestrial one, the sums exhibit the the combined effect of Rosenberg–Coleman and Russell–McPherron effects. Thus, in the geomagnetically active hours the negative B_y dominates during the first half of the year and the positive B_y dominates during the second half of the year. However, these effects can not be detected in the occurrence of the negative and positive GSM B_y values when the solar and terrestrial dipoles are parallel. In this case one can see polarity-independent semiannual variations instead of the polarity-dependent opposite annual variations. It is well-known that the B_y component modulates the energy transfer from the solar wind to the magnetosphere causing marked asymmetries in magnetospheric convective flow patterns at high latitudes. Our results hint that the occurrences of these asymmetries related to the ICMEs depend on the solar dipole cycle. In the antiparallel years one of them dominates during half a year causing asymmetric energy transfer to the magnetosphere. In the parallel years the occurrences of the two kind of asymmetries are equal on monthly time scale, thus the energy transfer is symmetric within a monthly and yearly time interval.     

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

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

Correction of B2bot for NeQuick during low solar activity at Hainan station

This work is a continuation of the previous article and it focuses on low solar activity and modeling effort. NeQuick model uses Epstein layer formalism to model each part of the profile. We study the diurnal and seasonal variations of B2_bot, ΔB2 (B2_best − B2_NeQuick2) and R (B2_best/B2_{NeQuick 2}) at Hainan station during low solar activity. The results show it is possible to improve the B2_bot parameter of the NeQuick model at that region during low solar activity. Then, we use a function ?(t) with LT in different seasons to correct the B2_bot formula of NeQuick 2. The correction shows that (1) By the correction formula, the B2_bot of NeQuick is improved. The maximum standard deviation is improved for 9 km. (2) The correction formula is more effective in summer than in equinox and winter and performs better during early morning hours than during the rest of the day.     

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.     

Determination of reconnected flux via remote sensing

Magnetic reconnection is one of the most fundamental processes in the magnetosphere. We present here a simple method to determine the essential parameters of reconnection such as reconnected flux and location of the reconnection site out of single spacecraft data via remote sensing. On the basis of a time-dependent reconnection model, the dependence of the reconnected flux on the magnetic field z-component B_z is shown. The integral of B_z over time is proportional to the reconnected flux and depends on the distance between the reconnection site and the actual position where B_z is measured. This distance can be estimated from analysis of magnetic field B_z data. We apply our method to Cluster measurements in the Earth’s magnetotail.     

Statistical survey on sawtooth events,SMCs and isolated substorms

Solar wind driving can cause a variety of different responses in the magnetosphere. Strong and steady driving during geomagnetic storms may result in sawtooth events. Strong to moderate driving may be followed by either sawtooth events or steady magnetospheric convection (SMC) events. Lower solar wind energy input typically leads to the formation of isolated non-storm substorms. This study uses superposed epoch analysis to reveal the typical properties of these three event groups as well as their similarities and differences. We use IMF and solar wind parameters, as well as ground-based indices (AL, SYM-H, ASY-H, PCN) to examine the level of solar wind driving and its response in the magnetosphere. Our results show that sawtooth events are associated with the strongest ionospheric activity. The subgroups of events during constant solar wind _EY$E_Y$ show that the key difference between the events is the average solar wind speed. Particularly, the high activity during sawtooth events is driven by high solar wind speed, while the lowest average speed during the SMCs may explain the lack of substorm activity during the steady convection periods.     

Dependence of the 27-day variation of cosmic rays on the global magnetic field of the Sun

We show that the higher range of the heliolongitudinal asymmetry of the solar wind speed in the positive polarity period (A > 0) than in the negative polarity period (A < 0) is one of the important reasons of the larger amplitudes of the 27-day variation of the galactic cosmic ray (GCR) intensity in the period of 1995–1997 (A > 0) than in 1985–1987 (A < 0). Subsequently, different ranges of the heliolongitudinal asymmetry of the solar wind speed jointly with equally important corresponding drift effect are general causes of the polarity dependence of the amplitudes of the 27-day variation of the GCR intensity. At the same time, we show that the polarity dependence is feeble for the last unusual minimum epoch of solar activity 2007–2009 (A < 0); the amplitude of the 27-day variation of the GCR intensity shows only a tendency of the polarity dependence. We present a three dimensional (3-D) model of the 27-day variation of GCR based on the Parker’s transport equation. In the 3-D model is implemented a longitudinal variation of the solar wind speed reproducing in situ measurements and corresponding divergence-free interplanetary magnetic field (IMF) derived from the Maxwell’s equations. We show that results of the proposed 3-D modeling of the 27-day variation of GCR intensity for different polarities of the solar magnetic cycle are in good agreement with the neutron monitors experimental data. To reach a compatibility of the theoretical modeling with observations for the last minimum epoch of solar activity 2007–2009 (A < 0) a parallel diffusion coefficient was increased by ∼40%.     

Features of the 11-year variation of galactic cosmic rays in different periods of solar magnetic cycles

Variations of galactic cosmic ray intensity have been studied based on the neutron monitors and interplanetary magnetic field experimental data for different ascending and descending epochs of solar activity. The dependence of the diffusion coefficient on the cosmic ray particles rigidity R is stronger in the maxima epoch than in the minima epoch of solar activity. For the period of 1977–1981 (qA > 0) the diffusion coefficient for the minimum epoch is, χ_min ∝ R^{0.7 ± 0.04} and for the maximum χ_max ∝ R^{1.3 ± 0.05}; for the period of 1987–1990 (qA < 0), χ_min ∝ R^{0.8 ± 0.05} and χ_max ∝ R^{1.1 ± 0.04}. The exponents ν_y and ν_z of the power spectral density of the B_y and B_z components of the IMF in the region of the frequencies (10⁻⁶– 4 × 10⁻⁶) Hz are larger for the minimum epoch of 1987 (ν_y ≈ 2.0 and ν_z ≈ 1.93) than for the maximum epoch of 1990 (ν_y ≈ 1.43 and ν_z ≈ 1.27).     

Moderate geomagnetic storms of January 22–25, 2012 and their influences on the wave components in ionosphere and upper stratosphere-mesosphere regions

Moderate geomagnetic storms occurred during January 22–25, 2012 period. The geomagnetic storms are characterized by different indices and parameters. The SYM-H value on January 22 increased abruptly to 67 nT at sudden storm commencement (SSC), followed by a sharp decrease to −87 nT. A second SSC on January 24 followed by a shock on January 25 was also observed. These SSCs before the main storms and the short recovery periods imply the geomagnetic storms are CME  -driven. The sudden jump of solar wind dynamic pressure and IMF _Bz$B_z$ are also consistent with occurrence of CMEs. This is also reflected in the change in total electron content (TEC) during the storm relative to quiet days globally. The response of the ionospheric to geomagnetic storms can also be detected from wave components that account for the majority of TEC variance during the period. The dominant coherent modes of TEC variability are diurnal and semidiurnal signals which account upto 83% and 30% of the total TEC variance over fairly exclusive ionospheric regions respectively. Comparison of TEC anomalies attributed to diurnal (DW1) and semidiurnal (SW2) tides, as well as stationary planetary waves (SPW1) at 12 UTC shows enhancement in the positive anomalies following the storm. Moreover, the impact of the geomagnetic storms are distinctly marked in the daily time series of amplitudes of DW1, SW2 and SPW1. The abrupt changes in amplitudes of DW1 (5 TECU) and SW2 (2 TECU) are observed within 20°S–20°N latitude band and along 20°N respectively while that of SPW1 is about 3 TECU. Coherent oscillation with a period of 2.4 days between interplanetary magnetic field and TEC was detected during the storm. This oscillation is also detected in the amplitudes of DW1 over EIA regions in both hemispheres. Eventhough upward coupling of quasi two day wave (QTDWs) of the same periodicity, known to have caused such oscillation, are detected in both ionosphere and upper stratosphere, this one can likely be attributed to the geomagnetic storm as it happens after the storm commencement. Moreover, further analysis has indicated that QTDWs in the ionosphere are strengthened as a result of coherent oscillation of interplanetary magnetic field with the same frequency as QTDWs. On the otherhand, occurrences of minor SSW and geomagnetic storms in quick succession complicated clear demarcation of attribution of the respective events to variability of QTDWs amplitudes over upper stratosphere.     

Dynamics of the convection in the inner magnetosphere by observations of the diffuse aurora,SAR arc and ionospheric drift

The first results of the comparison of subauroral luminosity dynamics in 557,7 and 630,0 nm emission with simultaneous measurements of the ionospheric drift in the F2 region with a digisonde DPS-4 at the Yakutsk meridian (CGMC: 55–60N, 200°E) at Kp = 2–6 are presented. It is shown from the analysis of individual events that during the magnetospheric convection intensification after the turn of the IMF B_z – component to the south the equatorward extension of diffuse aurora takes place. At the same time the westward ionospheric drift velocity increases both in the diffuse aurora region and much equatorward of it due to the occurrence of the northward polarization electric field. We suppose that the generation of polarization field can be associated with the development of the region 2 FAC during the intensification of magnetospheric convection. The comparison of ground-based observations with measurements of the plasma drift aboard the DMSP-F15 satellite has been carried out.     

Energy conversion of the flare due to direct electric fields from the sheared reconnection

In this paper we present a new mechanism of the main energy conversion of the solar flare. Since a flare inducing prominence (flux tube) rises V_z ? 300 km s⁻¹, the plasmas below it cannot continuously eject with Alfvén speeds of V_A = 3000 km s⁻¹ but probably with V_z ≈ ±100 km s⁻¹. Plasma up and downflows with V_A will within a short duration be blocked between the chromosphere where reconnected flux tubes are piling up, and the slowly rising flux rope. Hence the Petschek slow shock mechanism is difficult to be realized as a major energy converting mechanism.     

Energy dependence of the rigidity spectrum of Forbush decrease of galactic cosmic ray intensity

We show that rigidity spectrum of Forbush decrease (Fd) of galactic cosmic ray (GCR) intensity in September 9–23, 2005 clearly depends on energy. We calculated rigidity spectrum of the Fd based on the neutron monitors and Nagoya muon telescope channels’ data divided in three groups according to their cut off rigidities. We found that temporal changes of rigidity spectrum exponent γ are approximately similar for all cut off rigidity groups, but γ values are the larger the higher are cut off rigidities. We conclude that rigidity spectrum of Fd is hard for lower energy range and is soft for the higher energy range. We believe that an energy dependence of the power law rigidity spectrum of Fd is observed owing to the preferential convection–diffusion mechanism during Fd in September 9–23, 2005. It is a reflection of an influence of the temporal changes of the structure of the interplanetary magnetic field (IMF) turbulence in different range of frequency f during Fd. Particularly, a decisive role in formation of the character of the rigidity spectrum belongs to the changes of the exponent ν of the power spectral density (PSD) of the IMF turbulence (PSD ∝ f^−ν). The exponent ν is greater for high frequency region of the IMF turbulence (responsible for scattering of low rigidity particles of GCR), than for low frequency region of the IMF turbulence (being responsible for scattering of higher rigidity particles). Also, we challenge to estimate an existence of slab/2D structure of solar wind turbulence during the Fd in September 9–23, 2005 based on the distribution of average turbulence energy among the IMF’s components.     

Contribution of geometry of interaction between interplanetary and terrestrial magnetic fields into global magnetospheric state and geomagnetic activity

The paper presents results of our study of dependence of geomagnetic activity from geoeffective parameters taking into account mutual orientation of the interplanetary magnetic field, electric field of the solar wind and geomagnetic moment. We attract a reconnection model elaborated by us made allowance for changes of geometry of the solar wind–magnetosphere interaction during annual and diurnal motions of the Earth. We take as our data base the interplanetary magnetic field and solar wind velocity measured at 1 a.u. at ecliptic plane for the period of 1963–2005 and K_p, D_st, a_m indices. Taken as a whole a geoeffective parameter suggested by us explains 95% of observed variations of the indices. Changes of the geometric factor determined by mutual orientation of the solar wind electric field and geomagnetic moment explain larger than 75% of observed statistical variations of D_st and a_m indices. Based on our results we suggest a new explanation of semi-annual and UT variation of geomagnetic activity.