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.     

Analysis of characteristic rhythms of heliogeomagnetic activity and their roles in synchronization of rhythms of biological objects

Rhythms of heliogeomagnetic activity and their dynamics in the solar activity (SA) cycles are analyzed. Heliogeomagnetic rhythms are the mixture of a periodic signal and noise component, which differs from a periodic signal by a substantial suddenness. Therefore, the filtering methods optimized for smoothness were used to analyze these rhythms and to separate rhythm periodic components. Short-period rhythms of the Kp index of geomagnetic activity and Wolf numbers (W) have been revealed. The rhythm dynamics in SA cycles has been indicated. Clearly defined near-seven-day, semi-seven-day, and other components, which are harmonics of the solar rotation period, are present in the rhythms of the Kp index in contrast to W. A comparison of the Kp and Dst rhythms indicates that the seasonal and annual variations in these indices are substantially similar. It has been indicated that the rhythms of heliogeomagnetic activity could be the external synchronizers of biological rhythms, having integrated into the endogenous temporal structure of these rhythms in the process of evolution, and that the dependence of the population macroparameters on W is much weaker than on Kp.     

Propagation of internal gravity waves to the altitudes of the ionospheric <Emphasis Type="Italic">D</Emphasis> and dynamo regions in the seismic region (Kamchatka): Preliminary results

A spectral analysis of simultaneous diurnal variations in the E _z component of the quasi-static electric field in the near-Earth atmosphere, VLF radio noise, and the horizontal component of the geomagnetic field, observed at Kamchatka in September 1999, has been performed. These geophysical parameters are indirectly used to study wave processes in the near-Earth atmosphere and in the ionospheric D and dynamo regions within the band of periods of internal gravity waves (T = 0.5?3.5 h). The correlation method in the frequency region is used to analyze the interrelation between the wave processes in these atmospheric regions. The power cross-spectra of various pairs of geophysical parameters have been studied depending on meteorological, seismic, and geomagnetic activities. It is shown that the oscillations in the power spectra in the T ～ 1–1.5 h band of periods are caused by the sources of internal gravity waves in the near-Earth atmosphere and by the remote sources above the dynamo region of the ionosphere within the T ～ 1.5–3 h band of periods.     

Spectra of atmospheric waves in the dynamo region of the ionosphere at Kamchatka

A spectral analysis of the diurnal variations in the geomagnetic field horizontal component, observed at Kamchatka and Barrow polar observatory in September–October 1999, has been performed. The complete set of oscillations of thermal tidal atmospheric waves with T = 24, 12, 8, and 4 h has been detected in the variation spectral power (Sq) at Kamchatka, and only the fundamental harmonic with T = 24 h has been distinguished at Barrow. The above periods vary in both directions relative to stable maximums during strong geomagnetic disturbances. The relative spectral intensity at subharmonics also vary toward the fundamental harmonic with a period of 24 h. In the frequency band 0.5–3 h (IGW periods), the maximal intensity in the background spectra is observed at T ～ 2 h and increases by an order of magnitude with increasing geomagnetic activity at both Kamchatka and Barrow. A day before earthquakes, the intensity of this maximum is below the rms background values, and the spectra widen toward the region of periods shorter than 2 h. A similar effect was previously observed in the power spectra of the diurnal variations in the quasistatic electric field and VLF noise, simultaneously measured in September–October 1999.     

The relation between solar activity and orientation of the solar wind electric field relative to the Earth’s magnetic moment

The relation of the Kp index of geomagnetic activity to the solar wind electric field (E _SW) and the projection of this field onto the geomagnetic dipole has been estimated. An analysis indicated that the southward component of the IMF vector (B _z < 0) is the main geoeffective parameter, as was repeatedly indicated by many researchers. The presence of this component in any combinations of the interplanetary medium parameters is responsible for a high correlation between such combinations and geomagnetic activity referred to by the authors of different studies. Precisely this field component also plays the main role in the relation between the Kp index and the relative orientation of E _SW and the Earth’ magnetic moment.     

Effect of geomagnetic disturbances on the operation of railroad automated mechanisms and telemechanics

Geomagnetic variations generate electric currents in long conductors such as high-voltage lines, pipelines, and telecoms cables. The aim of our work is to study the possible effect of geomagnetic disturbances on the operation of automated systems and telemechanics of a midlatitude railroad based on the data on the malfunctions and breakdowns registered in 2004 on the East Siberian railroad (VSZhD). It has been obtained that the total daily duration of malfunctions and breakdowns (T) during disturbed periods is controlled by geomagnetic activity. When a peak of geomagnetic activity is reached during a storm, T increases about three times. Moreover, a correlation between T and the local index of geomagnetic activity (A), measured at Podkamennaya Tunguska Siberian observatory, is high during disturbed periods. Specifically, the correlation coefficient (K) is equal to 0.83 and 0.71 for the strongest two storms of 2004 that occurred in July 17–August 2 and November 4–18, respectively.     

Regular Variations of Intensity of the Westward Auroral Electrojet according to the Geomagnetic <Emphasis Type="Italic">AL</Emphasis> Index

Diurnal and seasonal variations of the geomagnetic AL index are studied. It is found in disturbed days that the mode of AL diurnal variation depends on the angle between the Sun–Earth line prolongation in the direction towards the magnetotail and the plane of geomagnetic equator; on quiet days, AL depends on the angle of attack between the geomagnetic axis and the Earth–Sun line. Seasonal AL variations are characterized by annual variations with summer maximum and semiannual variations with equinoctial maxima. It is shown that the semiannual AL variations can be described by a simplified model of plasma convection in the magnetotail based on a plasma electron cooling mechanism.     

Comparative analysis of atmospheric and ionospheric variability by measurements of temperature in the mesopause region and peak electron density <Emphasis Type="Italic">NmF</Emphasis>2

The results of studying the atmospheric and ionospheric variability in the region of Eastern Siberia are presented. The analysis involved data on the atmosphere temperature at mesopause heights (Tm) and vertical sounding data on the peak electron density (NmF2). The data on temperature were obtained by spectrometric observations of the hydroxyl molecule emission (band ОН (6-2), 834.0 nm, maximum emission height ~87 km). The analysis covers the period from 2008 to 2015. Seasonal and year-to-year variations in the variability of Tm and NmF2 were studied and compared in different time periods: day-to-day variations (T > 24 h), tidal variations (8 h ≤ T ≤ 24 h), and variations with periods of internal gravity waves (T < 8 h). Both common features and distinctions in the behavior of the analyzed parameters have been found, and their possible physical causes are analyzed.     

Magnetospheric effects of cosmic rays. 1. Long-term changes in the geomagnetic cutoff rigidities for the stations of the global network of neutron monitors

Vertical geomagnetic cutoff rigidities are obtained for the stations of the global network of neutron monitors via trajectory calculations for each year of the period from 1950 to 2020. Geomagnetic cutoff rigidities are found from the model of the Earth’s main field International Geomagnetic Reference Field (IGRF) for 1950–2015, and the forecast until 2020 is provided. In addition, the geomagnetic cutoff rigidities for the same period are obtained by Tsyganenko model T89 (Tsyganenko, 1989) with the average annual values of the Kp-index. In each case, the penumbra is taken into account in the approximation of the flat and power spectra of variations of cosmic rays. The calculation results show an overall decrease in geomagnetic cutoff rigidities, which is associated with the overall decrease and restructuring of the geomagnetic field during the reporting period, at almost all points.     

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.     

Relation of the <Emphasis Type="Italic">F</Emphasis>2 region to stratospheric parameters: Dependence on magnetic activity in two solar cycles

The dependence of the correlation coefficient r(h, fo) between the stratospheric parameter h(100) and critical frequency foF2 revealed in the data of two solar cycles (1979–1989 and 1990–2000) on geomagnetic activity is analyzed. It is shown that the character of the r(h, fo) dependence on limitation on the Ap geomagnetic index is the same in both cycles but depends on the time of day and solar activity level for the given year. It is also found that there is a considerable difference in the absolute values of r(h, fo) between two cycles.     

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.     

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.     

Geomagnetic disturbances and earthquakes in the Vrancea zone

Based on the example of the Vrancea zone of concentrated seismicity, it is shown how the stress-strain state of the medium responds to a disturbance of the geomagnetic field. Geomagnetic conditions are examined in relation to earthquakes in the Vrancea zone in the period 1988–1996. It is established that the seismic energy release in the Vrancea zone is associated with differences (“gradients”) in the H component of the geomagnetic field. Such a gradient preceding earthquakes is shown to be the midnight polar substorm and the degree of its mid-latitude effect. The time interval from the maximum of the substorm development to a shock (τ, h) is directly related to the focal depth. The seismic characteristics K _en and h (km) are demonstrated to be related to morphological features of the substorm development, namely, its duration T (min), intensity, and background. Differences in the duration of polar substorms before crustal (shallow) and deep earthquakes are revealed. Morphological features of the spectrum of geomagnetic variations preceding the seismic energy release are established.     

Spatial and time variations in the <Emphasis Type="Italic">foF</Emphasis>2 (night)/<Emphasis Type="Italic">foF</Emphasis>2 (day) ratio: Specification of some effects

A series of important aspects of the time and spatial variations in the critical frequency ratio fo(night)/fo(day) considered by one of the authors [Danilov, 2007, 2008] is specified. The list of the Eastern-Hemisphere stations, for which an analysis of the above indicated ratio was performed, is completed. The available stations of the Western Hemisphere are considered. It is shown that the character of the variations in the fo(night)/fo(day) ratio is independent of limitations imposed on the Ap index of magnetic activity for the analyzed days. The consideration of the long-term variations in geomagnetic activity using the annual mean value of the Ap index does not influence the principal conclusions of the study, although in some cases changes statistical characteristics of the obtained behavior of the fo(night)/fo(day) ratio after 1980. A comparison of the results, based on the interpretation of the trends of the fo(night)/fo(day) value with the changes in the stratosphere-ionosphere coupling from the 1980s to the 1990s obtained by the authors earlier, confirms the assumption that there occurs a systematic change in the zonal wind in the upper atmosphere.     

The geomagnetic field at the Proterozoic-Paleozoic boundary and lower-mantle plumes

Data on the amplitude of variations in the direction of the geomagnetic field and the frequency of reversals in the Vendian-Cambrian are presented. It has been established from these data that (a) distributions of variations in the direction of the geomagnetic field S _p are bimodal (modes 9° and 11°); (b) the maximum of the average amplitude S _p takes place by 5–10 Myr later than the Vendian-Cambrian boundary; (c) S _p tends to increase as plume epicenters are approached; and (d) the plume formation is more often confined to intervals with different frequencies of geomagnetic reversals than to the interval of a stable state of the geomagnetic field without reversals (Vendian hyperchron). The listed features of the geomagnetic field behavior are repeated near all boundaries of geological eras of the Phanerozoic.     

Dependence of geomagnetic disturbances on extreme values of the solar wind <Emphasis Type="Italic">E</Emphasis><Subscript><Emphasis Type="Italic">y</Emphasis></Subscript> component

The dependence of the zonal geomagnetic indices (AE, Ap, Kp, Kn, and Dst) on the solar wind parameters (the electric field E _y component, dynamic pressure P _d and IMF irregularity σB) has been studied for two types of events: magnetic clouds and high-speed streams. Based on the empirical relationships, it has been established that the AE, Ap, Kp, and Kn indices are directly proportional to the E _y value at E _y < 12 mV m^?1 and are inversely proportional to this value at E _y > 12 mV m^?1 for the first-type events. On the contrary, the dependence of Dst on E _y is monotonous nonlinear. A linear dependence of all geomagnetic indices on E _y is typical of the second-type events. It has been indicated that the specific features of geoeffectiveness of magnetic clouds and high-speed solar wind streams are caused by the dependence of the electric field potential across the polar cap on the electric field, solar wind dynamic pressure, and IMF fluctuations.     

Empirical Model of the Location of the Main Ionospheric Trough

The empirical model of the location of the main ionospheric trough (MIT) is developed based on an analysis of data from CHAMP satellite measured at the altitudes of ~350–450 km during 2000–2007; the model is presented in the form of the analytical dependence of the invariant latitude of the trough minimum Φ_m on the magnetic local time (MLT), the geomagnetic activity, and the geographical longitude for the Northern and Southern Hemispheres. The time-weighted average index Kp(τ), the coefficient of which τ = 0.6 is determined by the requirement of the model minimum deviation from experimental data, is used as an indicator of geomagnetic activity. The model has no limitations, either in local time or geomagnetic activity. However, the initial set of MIT minima mainly contains data dealing with an interval of 16–08 MLT for Kp(τ) < 6; therefore, the model is rather qualitative outside this interval. It is also established that (a) the use of solar local time (SLT) instead of MLT increases the model error no more than by 5–10%; (b) the amplitude of the longitudinal effect at the latitude of MIT minimum in geomagnetic (invariant) coordinates is ten times lower than that in geographical coordinates.     

The effect of baric disturbances in the atmosphere on microseismic processes in the crust

The effect of baric variations of different origins on characteristics of seismic noise is analyzed in the frequency range 0.03–20 Hz. Long period variations in atmospheric pressure caused by cyclones, whose period T ranges from half a day to a few days, are shown to increase the microseismic background amplitude by two to four times in the frequency range 0.03–1 Hz (the coefficient of linear correlation between time variations in the amplitude and atmospheric pressure is K = 0.65 at a significance level of r = 0.95). Short-period baric variations with T ～ 5–30 min associated with the passage of cold fronts lead to a tenfold increase in the microseismic background amplitude in the frequency range 4–8 Hz (K = 0.67 at r = 0.95). In this case, disturbances of seismic background are recorded for 20–60 min after the passage of an atmospheric front and display an exponential drop in the amplitude. In distinction to cyclones, an atmospheric front increases the number of impulsive microseismic events of the resonance type.     

Precursor times of abnormal <Emphasis Type="Italic">b</Emphasis>-values prior to mainshocks

Seismic observations exhibit the presence of abnormal b-values prior to numerous earthquakes. The time interval from the appearance of abnormal b-values to the occurrence of mainshock is called the precursor time. There are two kinds of precursor times in use: the first one denoted by T is the time interval from the moment when the b-value starts to increase from the normal one to the abnormal one to the occurrence time of the forthcoming mainshock, and the second one denoted by T _p is the time interval from the moment when the abnormal b-value reaches the peak one to the occurrence time of the forthcoming mainshock. Let T* be the waiting time from the moment when the abnormal b-value returned to the normal one to the occurrence time of the forthcoming mainshock. The precursor time, T (usually in days), has been found to be related to the magnitude, M, of the mainshock expected in a linear form as log(T)?=?q?+?rM where q and r are the coefficient and slope, respectively. In this study, the values of T, T _p , and T* of 45 earthquakes with 3?≤?M?≤?9 occurred in various tectonic regions are compiled from or measured from the temporal variations in b-values given in numerous source materials. The relationships of T and T _p , respectively, versus M are inferred from compiled data. The difference between the values of T and T _p decreases with increasing M. In addition, the plots of T*/T versus M, T* versus T, and T* versus T-T* will be made and related equations between two quantities will be inferred from given data.