Effect of Moon phases in riometer absorption and in the ionospheric and geomagnetic parameters

Variations in the frequency of occurrence of riometer absorption, minimum frequency of reflection of the ionospheric F layer, minimum height, and height of maximum electron density of the ionospheric F layer near the solar minimum have been studied. Application of the superposed epoch technique has detected the Moon phase effect on these ionospheric parameters. This effect was: three events per day in the occurrence of riometer absorption, 0.056 MHz in the minimum frequency of reflection of the F layer, and 2.6 and 6.7 km, in the change of the minimum height of reflection and height of reflection from the region with maximum electron density of the ionospheric F layer, respectively. The lunar tide action changes the ionospheric conductivity and, thus, influences the current systems of the magnetosphere. Through changes of magnetospheric currents, the Moon phase effect is exhibited in the Ap and Dst indices and is 4.3 and 4.25 nT, respectively.     

Variations of Geomagnetic Cosmic Ray Thresholds and Their Latitudinal Behavior in the Period of Solar Disturbance in September 2005

The period of interplanetary, geomagnetic and solar disturbances of September 7–15, 2005, is characterized by two sharp increases of solar wind velocity to 1000 km/s and great Dst variation of the geomagnetic field (~140 nT). The time variations of theoretical and experimental geomagnetic thresholds observed during this strong geomagnetic storm, their connection with solar wind parameters and the Dst index, and the features of latitudinal behavior of geomagnetic thresholds at particular times of the storm were studied. The theoretical geomagnetic thresholds were calculated with cosmic ray particle tracing in the magnetic field of the disturbed magnetosphere described by Ts01 model. The experimental geomagnetic thresholds were specified by spectrographic global survey according to the data of cosmic ray registration by the global station network.     

Dynamics of the ring current–magnetotail currents relationships during geomagnetic storms of different intensity

The results of studying the intensity of fluxes of 30–80 keV ions from the data of measurements of the NOAA (POES) sun-synchronous satellites during geomagnetic storms of different intensity are presented. For 15 geomagnetic storms with |Dst|_max from ~37 to ~422 nT, the storm-time maximum ion fluxes in the near-equatorial region (trapped particles) and at high latitudes (precipitating particles) have been considered. It is shown that the maximum fluxes of trapped particles, which are considered a ring-current proxy, increase with the storm power. In this case, if a smooth growth of fluxes is recorded for storms with |Dst|_max < 250 nT in the near-equatorial region, a significantly steeper growth of fluxes of trapped particles is observed when storm power increases during storms with |Dst|_max > 250 nT. This may be evidence of both an increasing of the contribution of the ring current relative to magnetotail currents to the development of high-intensity storms and to a nonlinear link between the ring current and ion fluxes at low altitudes in the near-equatorial region. Despite large variations in fluxes of precipitating particles in the polar region above the boundary of isotropization, a decreasing tendency, as a whole, in fluxes of these particles is observed with increasing the storm intensity. This is the evidence of the effect of saturation of magnetotail currents and of an increase in the relative role of the ring current during strong magnetic storms.     

Estimation of storm-time level of day-side wave geomagnetic activity using a new <Emphasis Type="Italic">ULF</Emphasis> index

The level of wave geomagnetic activity in the morning and daytime sectors of auroral latitudes during strong magnetic storms with Dst _min varying from ?100 to ?150 nT in 1995–2002 have been studied using a new ULF index of wave activity proposed in [Kozyreva et al., 2007]. It has been detected that daytime Pc5 pulsations (2–6 mHz) are most intense during the main phase of a magnetic storm rather than during the recovery phase as was considered previously. It has been indicated that morning geomagnetic pulsations during the substorm recovery phase mainly contribute to daytime wave activity. The appearance of individual intervals with the southward IMF B _z component during the magnetic storm recovery phase results in increases in the ULF index values.     

Daytime quasiperiodic geomagnetic pulsations during the recovery phase of the strong magnetic storm of May 15, 2005

Intense quasimonchromatic geomagnetic pulsations with a period of ～15 min, observed on the Earth’s surface in the near-noon sector at the beginning of the recovery phase of a very strong (Dst _min = ?260 nT) magnetic storm of May 15, 2005, are analyzed. The variations were registered at auroral latitudes only in the X field component, and wave activity shifted into the postnoon sector of the polar cap an hour later; in this case pulsations were observed in the X and Y field components. Within the magnetosphere the source of magnetic pulsations could be the surface waves on the magnetopause caused by the pulse of the solar wind magnetic pressure. Geomagnetic pulsations in the polar cap, observed in phase at different latitudes, could apparently reflect quasiperiodic variations in the NBZ system of field-aligned currents. Such variations can originate due to the series of pulsed reconnections in the postnoon outer cusp at large (～20 nT) positive B _z values and large (about ?40 nT) negative values of IMF B _x .     

An Empirical Model for the Ionospheric Global Electron Content Storm-Time Response

By analyzing the variations of global electron content (GEC) during geomagnetic storm events, the ratio “GEC/GECQT” is found to be closely correlated with geomagnetic Kp index and time weighted Dst index, where GECQT is the quiet time reference value. Moreover, the GEC/GECQT will decrease with the increase of the solar flux F10.7 index. Furthermore, we construct a linear model for storm-time response of GEC. Eighty-two storm events during 1999–2011 were utilized to calculate the model coefficients, and the performance of the model was tested using data of 8 storm events in 2012 by comparing the outputs of the model with the observed GEC values. Results suggest that the model can capture the characteristics of the GEC variation in response to magnetic storms. The component describing the solar activity influence shows a counteracting effect with the geomagnetic activity component; and the influence of Kp index causes an increase of GEC, while the time weighted Dst index causes a decrease of GEC.     

Global empirical model of critical frequency of the ionospheric <Emphasis Type="Italic">F</Emphasis>2-layer for quiet geomagnetic conditions

Using the foF2 database obtained from satellites and ground-based ionospheric stations, we have constructed a global empirical model of the critical frequency of the ionospheric F2-layer (SDMF2—Satellite and Digisonde Data Model of the F2 layer) for quiet geomagnetic conditions (Kp < 3). The input parameters of this model are the geographical coordinates, UT, day, month, year, and the integral index F10.7 (day, τ = 0.96) of solar activity for a given day. The SDMF2 model was based on the Legendre method for the spatial expansion of foF2 monthly medians to 12 in latitude and 8 in longitude of spherical harmonics. The resulting spatial coefficients have been expanded by the Fourier method in three spherical harmonics with respect to UT. The effect of the saturation of critical frequency of the ionospheric F2-layer at high solar activity was described in the SDMF2 model by foF2 as a logarithmic function of F10.7 (day, τ = 0.96). The difference between the SDMF2 and IRI models is a maximum at low solar activity as well as in the Southern Hemisphere and in the oceans. The testing on the basis of ground-based and satellite data has indicated that the SDMF2 model is more accurate than the IRI model.     

Dynamics of the ionospheric electric currents and auroral luminosity boundaries during strong magnetic storms

The regularities in the southward drift of the ionospheric current centers and luminosity boundaries during strong magnetic storms of November 2003 and 2004 (with Dst ≈ ?400 and ?470 nT, respectively) are studied based on the global geomagnetic observations and TV measurements of auroras. It has been indicated that the eastward and westward electrojets in the dayside and nightside sectors simultaneously shift equatorward to minimal latitudes of Φ _min ^° ～53°–55°. It has been obtained that the Φ _min ^° latitude decreases with increasing negative values of Dst, IMF B _z component, and westward electric field strength in the solar wind. The dependence of the electrojet equatorward shift velocity (V _av) on the rate of IMF B _z variations (ΔB _z /Δt) has been determined. It is assumed that the electrojet dynamics along the meridian is caused by a change in the structure of the magnetosphere and electric fields in the solar wind and the Earth’s magnetosphere.     

Peculiarities of the outer radiation belt dynamics during the strong magnetic storm of May 15, 2005

The dynamics of energetic electrons (E _e =0.17–8 MeV) and protons (E _p =1 MeV) of the outer radiation belt during the magnetic storm of May 15, 2005, at high (GOES-10 and LANL-84 geosynchronous satellites) and low (Meteor-3M polar satellite) altitudes is analyzed. The data have been compared to the density, plasma velocity, solar wind, and magnetic field measurements on the ACE satellite and geomagnetic disturbances. During the magnetic storm main phase, the nighttime boundary of the region of trapped radiation and the center of westward electrojet shifted to L ～ 3. Enhancements of only low-energy electrons were observed on May 15, 2005. The belt of relativistic electrons with a maximum at L ～ 4 was formed during the substorm, the amplitude of which reached its maximum at ～0630 UT on May 16. The results are in good agreement with the regularity relating the position of a maximum of the new relativistic electron belt, boundaries of the trapped radiation region, and extreme low-latitude position of westward electrojet center to the Dst variation amplitude.     

Ionospheric disturbances during the severe magnetic storm of November 7–10, 2004

Results of the study of the behavior of the F ₂ region and topside ionosphere during the magnetic storm on November 7–10, 2004, which was a superposition of two sequent Severe magnetic disturbances (Kp = 9–) are presented. The observations were conducted by the incoherent scatter radar at Kharkov. Considerable effects of a negative ionospheric disturbance are registered, including a decrease in the electron density in the F ₂-layer maximum by a factor of 6–7 and of the total electron content up to a height of 1000 km by a factor of 2, a lifting up of the ionospheric F ₂ layer by 300 km at night and by 150–180 km in the daytime, unusual nighttime heating of the plasma with an increase of the ion and electron temperatures up to 2000 and 3000 K, respectively, and a decrease in the relative density of hydrogen ions N(H⁺)/N _e by a factor of up to 3.5 because of the emptying of the magnetic flux tube passing over Kharkov. The effects usually observed in the high-latitude ionosphere, including the coherent echoes, are detected during the main phase of the storm. The results obtained manifest a shift of the large-scale structures of the high-latitude ionosphere (the auroral oval, main ionospheric trough, hot zone, etc.) down to latitudes close to the latitude of the Kharkov radar.     

Properties of solar activity and ionosphere for solar cycle 25

Based on the known forecast of solar cycle 25 amplitude (Rz _max ≈ 50), the first assessments of the shape and amplitude of this cycle in the index of solar activity F10.7 (the magnitude of solar radio flux at the 10.7 cm wavelength) are given. It has been found that (F10.7)_max ≈ 115, which means that it is the lowest solar cycle ever encountered in the history of regular ionospheric measurements. For this reason, many ionospheric parameters for cycle 25, including the F2-layer peak height and critical frequency (hmF2 and foF2), will be extremely low. For example, at middle latitudes, typical foF2 values will not exceed 8–10 MHz, which makes ionospheric heating ineffective in the area of upper hybrid resonance at frequencies higher than 10 MHz. The density of the atmosphere will also be extremely low, which significantly extends the lifetime of low-orbit satellites. The probability of F-spread will be increased, especially during night hours.     

Daytime geomagnetic disturbances at high latitudes during a strong magnetic storm of June 21–23, 2015: The storm initial phase

The high-latitude geomagnetic effects of an unusually long initial phase of the largest magnetic storm (SymH ~–220 nT) in cycle 24 of the solar activity are considered. Three interplanetary shocks characterized by considerable solar wind density jumps (up to 50–60 cm^–3) at a low solar wind velocity (350–400 km/s) approached the Earth’s magnetosphere during the storm initial phase. The first two dynamic impacts did not result in the development of a magnetic storm, since the IMF Bz remained positive for a long time after these shocks, but they caused daytime polar substorms (magnetic bays) near the boundary between the closed and open magnetosphere. The magnetic field vector diagrams at high latitudes and the behaviour of high-latitude long-period geomagnetic pulsations (ipcl and vlp) made it possible to specify the dynamics of this boundary position. The spatiotemporal features of daytime polar substorms (the dayside polar electrojet, PE) caused by sudden changes in the solar wind dynamic pressure are discussed in detail, and the singularities of ionospheric convection in the polar cap are considered. It has been shown that the main phase of this two-stage storm started rapidly developing only when the third most intense shock approached the Earth against a background of large negative IMF Bz values (to–39 nT). It was concluded that the dynamics of convective vortices and the related restructing of the field-aligned currents can result in spatiotemporal fluctuations in the closing ionospheric currents that are registered on the Earth’s surface as bay-like magnetic disturbances.     

Scatter of <Emphasis Type="Italic">hmF</Emphasis>2 values as an indicator of trends in the thermospheric dynamics

The variability degree of the F ₂-layer height, hmF2, from the 1950s–1960s to the 1990s has been analyzed based on the vertical sounding data for a series of midlatitude ionospheric stations. It has been found that the scatter of the hmF2 values (standard deviation) abruptly increases from the earlier decades to the later ones. This increase is more evident in the spring period of the year and is independent of geomagnetic activity. An increase in the scatter of hmF2 apparently indicates systematic changes (trends) in the thermospheric dynamics, the existence of which was suggested in the recent publications of the authors.     

Properties of the <Emphasis Type="Italic">F</Emphasis>2-layer critical frequency median in the nocturnal subauroral ionosphere during low and moderate solar activity

Based on an analysis of data from the European ionospheric stations at subauroral latitudes, it has been found that the main ionospheric trough (MIT) is not characteristic for the monthly median of the F2-layer critical frequency (foF2), at least for low and moderate solar activity. In order to explain this effect, the properties of foF2 in the nocturnal subauroral ionosphere have been additionally studied for low geomagnetic activity, when the MIT localization is known quite reliably. It has been found that at low and moderate solar activity during night hours in winter, the foF2 data from ionospheric stations are often absent in the MIT area. For this reason, a model of the foF2 monthly median, which was constructed from the remaining data of these stations, contains no MIT or a very weakly pronounced MIT.     

Great geomagnetic storms in 1841–1870 according to the data from the network of Russian geomagnetic observatories

Great magnetic storms (geomagnetic index C9 is ≥8 for St. Petersburg, which can correspond to Kp ≥ 8 or Dst < ?200 nT), registered from 1841 to 1870 at the St. Petersburg, Yekaterinburg, Barnaul, Nerchinsk, Sitka, and Beijing (at the Russian embassy) observatories are analyzed. A catalog of intensive magnetic storms during this period, which includes solar cycles 9–11, has been compiled. The statistical characteristics of great magnetic storms during this historical period have been obtained. These results indicate that high solar activity played a decisive role in the generation of very intense magnetic storms during the considered period. These storms are characterized by only one peak in a solar cycle, which was registered in the years of the cycle minimum (or slightly earlier): the number of great geomagnetic storms near the solar activity maximum was twice as large as the number of such storms during less active periods. A maximum in September–October and an additional maximum in February are observed in the annual distribution of storms. In addition, the storm intensity inversely depends on the storm duration.     

Effect of zonal E × B plasma drift on electron density in the low-latitude ionospheric <Emphasis Type="Italic">F</Emphasis> region at a solar activity maximum near vernal equinox

The variations in the density of the ionospheric F2 layer maximum (NmF2) under the action of the zonal plasma drift perpendicularly to the magnetic (B) and electric (E) fields in the direction geomagnetic west-geomagnetic east have been studied using the three-dimensional nonstationary theoretical model of electron and ion densities (N _e and N _i ) and temperatures (T _e and T _i ) in the low-latitude and midlatitude ionospheric F region and plasmasphere. The method of numerical calculations of N _e , N _i , T _e , and T _i , including the advantages of the Lagrangian and Eulerian methods, is used in the model. A dipole approximation of the geomagnetic field (B), taking into account the non-coincidence of the geographic and geomagnetic poles and differences between the positions of the Earth’s and geomagnetic dipole centers, is accepted in the calculations. The calculated NmF2 and altitudes of the F2 layer maximum (hmF2) have been compared with these quantities measured at 16 low-latitude ionospheric sounding stations during the geomagnetically quiet period October 11–12, 1958. This comparison made it possible to correct the input model parameters: the NRLMSISE-00 model [O], the meridional component of the neutral wind velocity according to the HWW90 model, and the meridional component of the equatorial plasma drift due to the electric field specified by the empirical model. It has been indicated that the effect of the zonal E × B plasma drift on NmF2 can be neglected under daytime conditions and changes in NmF2 and hmF2 under the action of this drift are insignificant under nighttime conditions north of 25° and south of ?26° geomagnetic latitude. The effect of the zonal E × B plasma drift on NmF2 and hmF2 is most substantial in the nightside ionosphere approximately from ?20° to 20° geomagnetic latitude, and the neglect of this drift results in an up to 2.4-fold underestimation of NmF2. The found dependence of the effect of the zonal E × B plasma drift on NmF2 and hmF2 on geomagnetic latitude is related to the longitudinal asymmetry of B, asymmetry of the neutral wind about the geomagnetic equator, and changes in the meridional E × B plasma drift at a change in geomagnetic longitude.     

Effect of the IMF azimuthal component on the longitudinal position of the westward electrojet intensity maximum during strong geomagnetic storms

The distribution of the ionospheric currents during the geomagnetic storms of November 20–21, 2003, November 7–8, 2004, and November 9–10, 2004, depending on the IMF B _y component, has been studied based on the data from the global network of magnetic stations. It has been indicated that, during geomagnetic disturbances, the westward electrojet intensity maximum is localized in the evening sector at IMF B _y < 0 and in the morning sector at IMF B _y > 0. The region of the westward electrojet intensity maximum shifts to morning hours with increasing positive B _y values. Thus, the IMF azimuthal component forms not only the magnetospheric convection pattern during magnetic storms but is also responsible for the longitudinal position of ionospheric structures.     

Formation Mechanisms of the Spring–Autumn Asymmetry of the Midlatitudinal <Emphasis Type="Italic">NmF</Emphasis>2 under Daytime Quiet Geomagnetic Conditions at Low Solar Activity

Formation mechanism of the spring–autumn asymmetry of the F2-layer peak electron number density of the midlatitudinal ionosphere, NmF2, under daytime quiet geomagnetic conditions at low solar activity are studied. We used the ionospheric parameters measured by the ionosonde and incoherent scatter radar at Millstone Hill on March 3, 2007, March 29, 2007, September 12, 2007, and September 18, 1984. The altitudinal profiles of the electron density and temperature were calculated for the studied conditions using a one-dimensional, nonstationary, ionosphere–plasmasphere theoretical model for middle geomagnetic latitudes. The study has shown that there are two main factors contributing to the formation of the observed spring–autumn asymmetry of NmF2: first, the spring–autumn variations of the plasma drift along the geomagnetic field due to the corresponding variations in the components of the neutral wind velocity, and, second, the difference between the composition of the neutral atmosphere under the spring and autumn conditions at the same values of the universal time and the ionospheric F2-layer peak altitude. The seasonal variations of the rate of O⁺(⁴S) ion production, which are associated with chemical reactions with the participation of the electronically excited ions of atomic oxygen, does not significantly affect the studied NmF2 asymmetry. The difference in the degree of influence of O⁺(⁴S) ion reactions with vibrationally excited N₂ and O₂ on NmF2 under spring and autumn conditions does not significantly change the spring–autumn asymmetry of NmF2.     

Statistical Analysis and Modeling of the Local Ionospheric Critical Frequency: A Mid-latitude Single-Station Model for Use in Forecasting

The hourly values of the F-layer critical frequency from the ionospheric sounder in Dourbes (50.1°N, 4.6°E) during the time interval from 1957 to 2010, comprising five solar cycles, were analyzed for the effects of the solar activity. The hourly time series were reduced to hourly monthly medians which in turn were used for fitting a single station foF₂ monthly median model. Two functional approaches have been investigated: a statistical approach and a spectral approach. The solar flux F_10.7 is used to model the dependence of foF₂ on the solar activity and is incorporated into both models by a polynomial expression. The statistical model employs polynomial functions to fit the F-layer critical frequency while the spectral model is based on spectral decomposition of the measured data and offers a better physical interpretation of the fitting parameters. The daytime and nighttime foF₂ values calculated by both approaches are compared during high and low solar activity. In general, the statistical model has a slightly lower uncertainty at the expense of the larger number of fitting parameters. However, the spectral approach is superior for modeling the periodic effects and performs better when comparing the results for high and low solar activity. Comparison with the International Reference Ionosphere (IRI 2012) shows that both local models are better at describing the local values of the F-layer critical frequency.     

Position of projections of the nightside auroral oval equatorward and poleward edges in the magnetosphere equatorial plane

The position of the auroral oval poleward and equatorward boundary projections on the equatorial plane in the nightside MLT sector during magnetically quiet periods (|AL| < 200 nT, |Dst| < 10 nT) has been determined. The oval boundary positions were determined according to the precipitation model developed at Polar Geophysical Institute (http://apm.pgia.ru/). The isotropy of the averaged plasma pressure and the experimentally confirmed balance of pressures during the nighttime have been taken into account. The morphological mapping method has been used to map the oval poleward and equatorward edges without the use of any magnetic field model on the assumption that the condition of magnetostatic equilibrium is valid. Ion pressures at ionospheric altitudes and in the equatorial plane have been compared. It has been shown that the auroral oval equatorward boundary in the midnight sector is localized at geocentric distances of ~7 R_E, which is in good agreement with the position of the energetic particle injection boundary in the equatorial plane. The oval poleward edge is localized at the ~10 R_E geocentric distance, which is in good agreement with the position of the equatorward boundary of the region with a high turbulence level in the Earth’s magnetosphere plasma sheet.