Ionospheric Time-delay over Akure Using Global Positioning System Observations

Ionospheric time delay (VΔt) variability using Global Positioning System (GPS) data over Akure (7.15°N, 5.12°E), Nigeria, has been studied. The observed variability of VΔt in comparison to older results of vertical total electron content (TEC) across similar regions has shown equivalent signatures. Higher monthly mean values of VΔt (MVΔt) were observed during daytime as compared to nighttime (pre- and post-midnight) hours in all months. The highest MVΔt observed in September during daytime hours range between ~6 and ~21 ns (~1.80 and ~6.30 m) and at post-midnight, they are in the range of ~1 to ~6 ns (~0.3 to ~1.80 m). The possible mechanisms responsible for this variability were discussed. Seasonal VΔt were investigated as well.     

IRI model in the problem of predicting ionospheric radio-wave propagation under conditions of high solar activity

The paper presents the results of an analysis of the correspondence between model representations the monthly mean diurnal dependences of critical frequency and vertical profiles of plasma frequency at local noon at IZMIRAN station for the middle months of the four seasons of 2014, the year of the maximum solar activity in the current 24th cycle. It is shown that in general the IRI model reliably describes the daily variation of foF2, and the smallest discrepancy is achieved when its basic input parameter is given by the ionospheric index of solar activity IG12. An exception is April, for which there is a fundamental discrepancy with the model both in the daily variation of the critical frequency foF2 and in the N _e (h)-profile for local noon time. For this month, the inadequacy in the model representation of the vertical distribution of the electron density turned out to be very significant in the calculation of the MUF: the relative error can reach 20%. The simulation results are confirmed by data from oblique-incidence ionospheric radio sounding.     

Modeling of Physical Processes by Analysis of Hard X-Ray and Microwave Radiations in the Solar Flare of November 10, 2002

For electron acceleration during solar flares, it is very important to determine the pitch-angle and energy dependences of the electron distribution function. At present, this cannot be done directly from observations. Therefore, it is necessary to perform a numerical simulation of the propagation of accelerated electrons in the magnetic field of the flare loop (loops) and calculate the X-ray and radio emissions. For the solar flare of November 10, 2002, we have obtained qualitative and quantitative agreements of modeled X-ray and radio maps with the RHESSI satellite and Nobeyama Radioheliograph data. We have determined the flare model parameters that agree with observations. The pitch-angle anisotropy of electrons determined by highly directional functions of the S(α) = cos⁸(α) type, the energy spectrum consist of two electron populations, the low-energy part of the spectrum up to an energy of break of 350 keV is characterized by a power law with the exponent δ₁ = 2.7–2.9, and the energy spectrum is more rigid above 420 keV (δ₂ = 2–2.3).     

The seismic wave absorption in the crust and upper mantle in the vicinity of the Kislovodsk seismic station

The Q-factor estimates of the Earth’s crust and upper mantle as the functions of frequency (Q(f)) are obtained for the seismic S-waves at frequencies up to ~35 Hz. The estimates are based on the data for ~40 earthquakes recorded by the Kislovodsk seismic station since 2000. The magnitudes of these events are M_W > 3.8, the sources are located in the depth interval from 1 to 165 km, and the epicentral distances range from ~100 to 300 km. The Q-factor estimates are obtained by the methods developed by Aki and Rautian et al., which employ the suppression of the effects of the source radiation spectrum and local site responses in the S-wave spectra by the coda waves measured at a fixed lapse time (time from the first arrival). The radiation pattern effects are cancelled by averaging over many events whose sources are distributed in a wide azimuthal sector centered at the receiving site. The geometrical spreading was specified in the form of a piecewise-continuous function of distance which behaves as 1/R at the distances from 1 to 50 km from the source, has a plateau at 1/50 in the interval from 50–70 km to 130–150 km, and decays as \({\raise0.7ex\hbox{$1$} \!\mathord{\left/ {\vphantom {1 {\sqrt R }}}\right.\kern-\nulldelimiterspace} \!\lower0.7ex\hbox{${\sqrt R }$}}\) beyond 130–150 km. For this geometrical spreading model and some of its modifications, the following Q-factor estimates are obtained: Q(f) ~ 85f^0.9 at the frequencies ranging from ~1 to 20 Hz and Q(f) ~ 75f^1.0 at the frequencies ranging from ~1 to 35 Hz.     

Long-period irregular pulsations under the conditions of a quiet magnetosphere

Simultaneous observations of high-latitude long-period irregular pulsations at frequencies of 2.0–6.0 mHz (ipcl) and magnetic field disturbances in the solar wind plasma at low geomagnetic activity (Kp ~ 0) have been studied. The 1-s data on the magnetic field registration at Godhavn (GDH) high-latitude observatory and the 1-min data on the solar wind plasma and IMF parameters for 2011–2013 were used in an analysis. Ipcl (irregular pulsations continuous, long), which were observed against a background of the IMF Bz reorientation from northward to southward, have been analyzed. In this case other solar wind plasma and IMF parameters, such as velocity V, density n, solar wind dynamic pressure P = ρV² (ρ is plasma density), and strength magnitude B, were relatively stable. The effect of the IMF Bz variation rate on the ipcl spectral composition and intensity has been studied. It was established that the ipcl spectral density reaches its maximum (~10–20 min) after IMF Bz sign reversal in a predominant number of cases. It was detected that the ipcl average frequency (f) is linearly related to the IMF Bz variation rate (ΔBz/Δt). It was shown that the dependence of f on ΔBz/Δt is controlled by the α = arctan(By/Bx) angle value responsible for the MHD discontinuity type at the front boundary of magnetosphere. The results made it possible to assume that the formation of the observed ipcl spectrum, which is related to the IMF Bz reorientation, is caused by solar wind plasma turbulence, which promotes the development of current sheet instability and surface wave amplification at the magnetopause.     

Modeling of the <Emphasis Type="Italic">D</Emphasis>-region behavior in the polar ionosphere during the solar flare of April 5, 2004

On the basis of the ion chemistry theoretical model, the impact of a powerful solar flare on variations in the ion composition and electron density in the D region of the polar ionosphere is considered. Good agreemnt between the model profiles of the electron density N _e (h) and the experimental data obtained during the flare by the partial reflection method is found. It is shown that the decrease in the effective recombination coefficient observed during disturbances is explained by the depletion of the relative content of the rapidly recombining complex ion clusters.     

Behavior of Parameters of Enhancements in the Nighttime Electron Concentration in the Ionospheric <Emphasis Type="Italic">F</Emphasis>2 Layer

We have analyzed the behavior of the F2 layer parameters during nighttime periods of enhanced electron concentration by the results of vertical sounding of the ionosphere carried out with five-minute periodicity in Almaty (76°55′ E, 43°15′ N) in 2001–2012. The results are obtained within the frameworks of the unified concept of different types of ionospheric plasma disturbances manifested as variations in the height and half-thickness of the layer accompanied by an increase and decrease of N _m F2 at the moments of maximum compression and expansion of the layer. A good correlation is found between height h _Am , which corresponds to the maximum increase, and layer peak height h _m F, while h _Am is always less than h _m F. The difference between h _Am and h _m F linearly increases with increasing h _m F. Whereas the difference is ~38 km for h _m F = 280 km, it is ~54 km for h _m F = 380 km. Additionally, the correlation is good between the increase in the electron concentration in the layer maximum ΔN _m and the maximum enhancement at the fixed height ΔN; the electron concentration enhancement in the layer maximum is about two to three times lower than its maximum enhancement at the fixed height.     

Reaction of the high-latitude lower ionosphere to solar proton events from observations in the ELF range

The reaction of the lower ionosphere to the solar proton events that occurred in 2011–2012 is studied in this paper based on the results of measurements of the propagation velocity and the E _z /H _τ ratio of the low-frequency electromagnetic pulses (atmospherics) in the ELF range at the high-latitude observatories Lovozero and Barentsburg. With numerical modeling methods, it is shown that horizontal local irregularities of the lower ionosphere conductivity profile could be a cause of the splashes in the E _z /H _τ ratio observed in the experiment during the solar proton event of March 7, 2012, which was a unique event in both the proton flux value and energy.     

Attenuation coefficients of Rayleigh and <Emphasis Type="Italic">Lg</Emphasis> waves

Analysis of the frequency dependence of the attenuation coefficient leads to significant changes in interpretation of seismic attenuation data. Here, several published surface-wave attenuation studies are revisited from a uniform viewpoint of the temporal attenuation coefficient, denoted by χ. Theoretically, χ( f) is expected to be linear in frequency, with a generally non-zero intercept γ?=?χ(0) related to the variations of geometrical spreading, and slope dχ/df = π/Q _e caused by the effective attenuation of the medium. This phenomenological model allows a simple classification of χ( f) dependences as combinations of linear segments within several frequency bands. Such linear patterns are indeed observed for Rayleigh waves at 500–100-s and 100–10-s periods, and also for Lg from ~2 s to ~1.5 Hz. The Lg χ( f) branch overlaps with similar linear branches of body, Pn, and coda waves, which were described earlier and extend to ~100 Hz. For surface waves shorter than ~100 s, γ values recorded in areas of stable and active tectonics are separated by the levels of \(\gamma _{D} \approx 0.2 \times 10^{-3}\) s^???1 (for Rayleigh waves) and 8 ×10^???3 s^???1 (for Lg). The recently recognized discrepancy between the values of Q measured from long-period surface waves and normal-mode oscillations could also be explained by a slight positive bias in the geometrical spreading of surface waves. Similarly to the apparent χ, the corresponding linear variation with frequency is inferred for the intrinsic attenuation coefficient, χ _i , which combines the effects of geometrical spreading and dissipation within the medium. Frequency-dependent rheological or scattering Q is not required for explaining any of the attenuation observations considered in this study. The often-interpreted increase of Q with frequency may be apparent and caused by using the Q-based model of attenuation and following preferred Q( f) dependences while ignoring the true χ( f) trends within the individual frequency bands.     

Response of the night aurora to a negative sudden impulse

Data from the meridian scanning photometers of the NORSTAR network and all-sky cameras of the THEMIS network were used for a detailed study of the response of night auroras to the sharp decrease of the solar wind dynamic pressure on September 28, 2009. The decrease in dynamic pressure was accompanied by a corresponding depression of the magnetic field in the SYM-H index and the origin of a negative sudden impulse (SI) with a duration of 5–8 min and amplitude of 150–200 nT in the horizontal component of the magnetic field at stations of the night sector of the auroral zone. The magnetic impulse was preceded by a long calm magnetic period, although the IMF Bz-component was negative for ~1.5 hour before the SI ^–. The commencement of the SI ^–, which was determined by variations in the magnetic field at ~0650 UT, was accompanied by a sharp increase in the intensity of discrete forms of polar auroras in the midnight sector of the auroral zone and their fast propagation to the pole. Approximately 6–8 min after the SI ^–, the auroral intensity in the emissions, which were excited by the fluxes of precipitated electrons and protons, quickly began to decrease in the night sector. Analysis of the optical observations showed the two-stage character of the response of the night auroras to the SI ^– in the considered event: first, fast movement of the discrete aurora forms to the pole with a significant increase in their intensity, and a further fast decrease in auroral intensity with a delay of ~6–8 min relative to the SI ^–. The possible reasons for such aurora behavior are discussed.     

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.     

Features of Radiation and Propagation of Seismic Waves in the Northern Caucasus: Manifestations in Regional Coda

Based on the Anapa (ANN) seismic station records of ~40 earthquakes (M_W > 3.9) that occurred within ~300 km of the station since 2002 up to the present time, the source parameters and quality factor of the Earth’s crust (Q(f)) and upper mantle are estimated for the S-waves in the 1–8 Hz frequency band. The regional coda analysis techniques which allow separating the effects associated with seismic source (source effects) and with the propagation path of seismic waves (path effects) are employed. The Q-factor estimates are obtained in the form Q(f) = 90 × f 0.7 for the epicentral distances r < 120 km and in the form Q(f) = 90 × f1.0 for r > 120 km. The established Q(f) and source parameters are close to the estimates for Central Japan, which is probably due to the similar tectonic structure of the regions. The shapes of the source parameters are found to be independent of the magnitude of the earthquakes in the magnitude range 3.9–5.6; however, the radiation of the high-frequency components (f > 4–5 Hz) is enhanced with the depth of the source (down to h ~ 60 km). The estimates Q(f) of the quality factor determined from the records by the Sochi, Anapa, and Kislovodsk seismic stations allowed a more accurate determination of the seismic moments and magnitudes of the Caucasian earthquakes. The studies will be continued for obtaining the Q(f) estimates, geometrical spreading functions, and frequency-dependent amplification of seismic waves in the Earth’s crust in the other regions of the Northern Caucasus.     

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.     

Use of the index of TEC vertical variation disturbance in studying ionospheric effects of the Chelyabinsk meteorite

The results of an analysis of the ionospheric effects accompanying fall of the Chelyabinsk meteorite on February 15, 2013 are presented using a method of calculating the index of the disturbance of total electron content vertical variations (Wtec) according to data from the GPS receiver network. A substantial increase (by a factor of 2–3) in the Wtec index with a duration of ~1.5 h was observed in the studied region after the main height explosion accompanying the meteorite fall at 0320 UT. The ionospheric response in Wtec was most significant statistically registered at the radio rays “receiver–satellite” for the GPS located southward from the place of explosion.     

Disturbances in the topside ionosphere during solar flares

The results of studying the ionospheric response to solar flares, obtained from the data of the GPS signal observations and incoherent scatter radars and as a result of the model calculations, are presented. It is shown that, according to the GPS data, a flare can cause a decrease in the electron content at altitudes of the topside ionosphere (h > 300 km). Similar effects of formation of a negative disturbance in the ionospheric F region were also observed during the solar flares of May 21 and 23, 1967, with the Arecibo incoherent scatter radar. The mechanism by which negative disturbances appear in the topside ionosphere during solar flares has been studied in this work based on the theoretical model of the ionosphere-plasmasphere coupling. It has been indicated that the formation of the electron density negative disturbance in the topside ionosphere is caused by an intense removal of O⁺ ions into the overlying plasmasphere under the action of an abrupt increase in the ion production rate and thermal expansion of the ionospheric plasma.     

Oxygen cyclotron harmonic waves observed using Van Allen Probes

Fine structured multiple-harmonic electromagnetic emissions at frequencies around the equatorial oxygen cyclotron harmonics are observed by Van Allen Probe A outside the core plasmasphere(L~5) off the magnetic equator(MLAT~.7.5°)during a geomagnetic storm. We find that the multiple-harmonic emissions have power spectrum density(PSD) peaks during 2–8equatorial oxygen gyroharmonics( f ~ n fO+, n=2–8), while the fundamental mode(n=1) is absent, implying that the harmonic waves are generated near the equator and propagate into the observation region. Additionally, these electromagnetic emissions are linearly polarized. Different from the equatorial noise emission that propagates considerably obliquely, these emissions have moderate wave normal angles(approximately 40°–60°), which predominately increase as the harmonic number increases.Considering their frequency and wave normal angle characteristics, it is suggested that these multiple-harmonic emissions play an important role in the dynamic variation of radiation belt electrons.     

The <Emphasis Type="Italic">Q</Emphasis>-factor estimates for the crust and upper mantle in the vicinity of Sochi and Anapa (North Caucasus)

The regularities in the radiation and propagation of seismic waves in the regions of the North Caucasus are analyzed for estimating the ground motion parameters during the probable future strong earthquakes. Based on the records of the regional earthquakes with magnitudes M_W ~ 3.9–5.6 within epicentral distances up to ~300 km obtained during the period of digital measurements at the Sochi and Anapa seismic stations, the Q-factors in the vicinities of these sites are estimated at ~55 f^0.9 and ~90f^0.7, respectively. The estimates were obtained by the coda normalization method developed by Aki, Rautian, and other authors. This method is based on the phenomenon of suppression of the earthquake (source) effects and local (site) responses by coda waves in the S-wave spectra. The obtained Q-factor estimates can be used for forecasting the ground shaking parameters for the future probable strong earthquakes in the North Caucasus in the vicinities of Sochi and Anapa.     

Intersector disbalance of the large-scale open solar magnetic field fluxes,active and passive boundaries,and solar-terrestrial extrastorms

The dynamics (from rotation to rotation) of the absolute values of the large-scale open solar magnetic field fluxes in the four-sector field structure has been considered for the first time, using CRs 2032–2035 in July–October 2005 as examples. An important role of the ratio of the fluxes at the eastern and western sector boundaries (Φ _E /Φ _W ) is confirmed. As in the cases of the two-sector structure, Φ _E /Φ _W > 1 is typical of active rigidly corotating boundaries with intense sunspot formation, flares, and interplanetary and geomagnetic disturbances. A remarkable property of the considered structure was the presence of a rapidly increasing flux in an initially narrow sector and the flux interaction with a stable rigidly corotating sector in the zone of the main active longitudes, which caused an unexpectedly strong geoeffective long-range action of flares near the corresponding active boundary.     

Role of solar activity in formation of the anomalous El Nin’o current

Monthly indices of Southern Atmospheric Oscillation (SOI) and corresponding Wolf numbers, geoeffective solar flares, magnetic AE indices as well as daily average values of the southward component of the interplanetary magnetic field (IMF B _z) and data on the wind characteristics at Antarctic stations Vostok, Leningradskaya, and Russkaya are analyzed. It is shown that a sharp decrease in the SOI indices, which corresponds to the beginning of El Nin’o (ENSO), is preceded one or two months before by a 20% increase in the monthly average Wolf numbers. In warm years of Southern Atmospheric Oscillation a linear relationship is observed between the SOI indices and the number of geoeffective solar flares with correlation coefficients p < ?0.5. It is shown that in warm years a change in the general direction of the surface wind to anomalous at the above stations is preceded one or two days before by an increase in the daily average values of IMF B _z. An increase in the SOI indices is preceded one or two months before by a considerable increase in the monthly average values of the magnetic AE indices.     

Activation of seismicity in Central and South Asia after the Makran earthquakes: Possible acceleration of preparation of large seismic events in the Tien Shan region

Two zones of seismicity (ten events with M _w = 7.0–7.7) stretching from Makran and the Eastern Himalaya to the Central and EasternTien Shan, respectively, formed over 11 years after the great Makran earthquake of 1945 (M _w = 8.1). Two large earthquakes (M _w = 7.7) hit theMakran area in 2013. In addition, two zones of seismicity (M ≥ 5.0) occurred 1–2 years after theMakran earthquake in September 24, 2013, stretching in the north-northeastern and north-northwestern directions. Two large Nepal earthquakes struck the southern extremity of the “eastern” zone (April 25, 2015, M _w = 7.8 and May 12, 2015, M _w = 7.3), and the Pamir earthquake (December 7, 2015, M _w = 7.2) occurred near Sarez Lake eastw of the “western” zone. The available data indicate an increase in subhorizontal stresses in the region under study, which should accelerate the possible preparation of a series of large earthquakes, primarily in the area of the Central Tien Shan, between 70° and 79° E, where no large earthquakes (M _w ≥ 7.0) have occurred since 1992.