Search

In June 2014, a number of meteotsunamis were detected in the Mediterranean and Black Sea area. These meteotsunamis were initiated by a unique high-altitude dynamical system which was initially originated above Spain and traveled across the Mediterranean Sea towards Black Sea and Turkey. Meteotsunamis unlike tsunamis driven by strong earthquakes are local events, and their formation has different mechanism. Atmospheric internal gravity waves (IGWs) are one of the main known sources of meteotsunamis (e.g. Vilibic et al. in Pure Appl Geophys 165:2169–2195, 2008). The synoptic system produced short-lived and small-scale atmospheric pressure perturbations which drifted with the jet stream-like bubbles and generated tsunami-like waves in the open waters. The bubbles with typical dimensions 15–60 km continuously form and collapse in the atmosphere at altitudes of 3–6 km. Such a “boiled” atmosphere generated IGWs propagating both downward, where they produced meteotsunamis (presumably under Proudman resonance condition) and upward into the ionosphere, with following dissipation and excitation of plasma density perturbations. One of the few experimental techniques, which can monitor perturbations of the ionization within the lower ionosphere, uses long-wave probing by very low and low frequency (VLF/LF) radio signals. To study the ionospheric disturbances observed during the chain of meteotsunamis affecting the Mediterranean Sea, we used VLF/LF data collected in South Europe by “The International Network for Frontier Research on Earthquake Precursors”. By applying the spectral analysis method to the anomalous VLF/LF signals, it was found that revealed periods of the signal variations were from 10 to 40–70 min in different stations, which are in the range of the atmospheric pressure oscillations and the meteotsunami events. These periods also correspond to the periods of IGWs.

We analyze retrospectively/prospectively the transient variations of six different physical parameters in the atmosphere/ionosphere during the M7.8 and M7.3 earthquakes in Nepal, namely: 1) outgoing longwave radiation (OLR) at the top of the atmosphere (TOA); 2) GPS/TEC; 3) the very-low-frequency (VLF/LF) signals at the receiving stations in Bishkek (Kyrgyzstan) and Varanasi (India); 4) Radon observations; 5) Atmospheric chemical potential from assimilation models; and; 6) Air Temperature from NOAA ground stations. We found that in mid-March 2015, there was a rapid increase in the radiation from the atmosphere observed by satellites. This anomaly was located close to the future M7.8 epicenter and reached a maximum on April 21–22. The GPS/TEC data analysis indicated an increase and variation in electron density, reaching a maximum value during April 22–24. A strong negative TEC anomaly in the crest of EIA (Equatorial Ionospheric Anomaly) occurred on April 21, and a strong positive anomaly was recorded on April 24, 2015. The behavior of VLF-LF waves along NWC-Bishkek and JJY-Varanasi paths has shown abnormal behavior during April 21–23, several days before the first, stronger earthquake. Our continuous satellite OLR analysis revealed this new strong anomaly on May 3, which was why we anticipated another major event in the area. On May 12, 2015, an M7.3 earthquake occurred. Our results show coherence between the appearance of these pre-earthquake transient’s effects in the atmosphere and ionosphere (with a short time-lag, from hours up to a few days) and the occurrence of the 2015 M7.8 and M7.3 events. The spatial characteristics of the pre-earthquake anomalies were associated with a large area but inside the preparation region estimated by Dobrovolsky-Bowman. The pre-earthquake nature of the signals in the atmosphere and ionosphere was revealed by simultaneous analysis of satellite, GPS/TEC, and VLF/LF and suggest that they follow a general temporal-spatial evolution pattern that has been seen in other large earthquakes worldwide.

Diagnostics of the outer envelopes of the Earth through global satellite navigational systems and through the equally global network of ground receivers and very-low-frequency transmitters make it possible to explore fairly complex processes of lithosphere–ionosphere interactions. The effects of earthquakes and tsunamis on the ionosphere are an example of these processes. A brief review of the current advances in studying these processes is presented.

Natural disasters, such as earthquakes, tsunamis, volcanic eruptions, tropical cyclones (typhoons in the Pacific, hurricanes in America), often lead to enormous human casualties, causing great damage to the human environment. This is due to the unpredictable nature of disasters. Investigation of the precursor signs of the development of such dangerous events (including an assessment of their power) allows you to inform the relevant services in advance of a natural disaster. This is often a very urgent task. The experimental and theoretical study of the electromagnetic response of the ionosphere to the development of natural catastrophic events is an important component of furthering disaster preparedness. This includes improvement of methods for event precursor isolation for use in warning and development systems. In this, validation of lithosphere-ionosphere communication models are very important. Therefore, in this work, we consider the amplitude anomalies of VLF signals recorded at a station in Moscow (Russia) during the passage of an atmospheric/ionospheric seismic wave from an earthquake in Chile, which epicenter is removed at a distance of ~ 15 thousand km.The work is performed according to the Russian Government Program of Competitive Growth of Kazan Federal University.

Modeling propagation of VLF electromagnetic beams in the waveguide earth-ionosphere (WGEI) is of a great importance because variation in the characteristics of these waves is an effective instrument for diagnostics the influences on the ionosphere from above (Sun-Solar Wind-Magnetosphere-Ionosphere), from below (the most powerful meteorological, seismogenic and other sources in the lower atmosphere and lithosphere/Earth, such as hurricanes, earthquakes, tsunamis etc.), from inside the ionosphere (strong thunderstorms and lightning discharges) and even from the far space (such as gamma-flashes, cosmic rays etc.). Thus, VLF became one of the most universal instrument for monitoringthe Space Weather in the direct sense of this term, i.e. the state of the Sun-Earth space and the ionosphere as it is particularly determined by all possible relatively powerful sources, wherever they are placed. This paper is devoted mostly to modelling VLF electromagnetic beam propagation in the WGEI. We present a new tensor impedance method for modelling propagation of electromagnetic beams (TIMEB) in a multi-layered/inhomogeneous waveguide. Suppose that such a waveguide, i.e. WGEI, possesses the gyrotropy and inhomogeneity with a thick cover layer placed above the waveguide. Note a very useful and attractive feature of the proposed TIMEB method: in spite of a large thickness of the waveguide cover layer, the proposed effective impedance approach reflects an impact of such a cover on the electromagnetic (EM) waves, which propagate in the waveguide. This impedance approach can be applied for EM waves/beams in layered gyrotropic/anisotropic active media in very wide frequency range, from VLF to optics. Moreover, this approach can be applied to calculations of EM waves/beams propagation in the media of an artificial origin such as metamaterial microwave or optical waveguides. The results of the modelling the propagation of VLF beams in the WGEI are included. The qualitative comparison between the theory and experimental observation of increasing losses of VLF waves in the WGEI is discussed. The new proposed method and its further development allows the comparison with the results of the future rocket experiment. This method allows to model (i) excitation of the VLF modes in the WGEI and their excitation by the typical VLF sources, such as radio wave transmitters and lightning discharges and (ii) leakage of VLF waves/beams into the upper ionosphere/magnetosphere.

An experimental study of the phase and amplitude observations of sub‐ionospheric very low and low frequency (VLF/LF) signals is performed to analyze the response of the lower ionosphere during the 21 August 2017 total solar eclipse in the United States of America. Three different sub‐ionospheric wave paths are investigated. The length of the paths varies from 2,200 to 6,400 km, and the signal frequencies are 21.4, 25.2, and 40.75 kHz. The two paths cross the region of the total eclipse, and the third path is in the region of 40–60% of obscuration. None of the signals reveal any noticeable amplitude changes during the eclipse, while negative phase anomalies (from −33° to −95°) are detected for all three paths. It is shown that the effective reflection height of the ionosphere in low and middle latitudes is increased by about 3–5 km during the eclipse. Estimation of the electron density change in the lower ionosphere caused by the eclipse, using linear recombination law, shows that the average decrease is by 2.1 to 4.5 times.

In this paper we analysed sudden phase anomalies (SPAs) of VLF–LF signals recorded at Graz (Austria), Birr (Ireland) and Moscow (Russia) stations during two strong solar flares in September 2017. The first X-class 9.3 flare occurred on 6 September at 12:02 UT, and the second X-class 8.2 flare was observed on 10 September 2017 at 16:06 UT. Data from seven transmitters in a frequency range between 20 and 45 kHz are used for the analysis. The SPAs were observed in all middle-latitude paths (differently orientated) with path lengths from 350 to 7000 km. Solar X-ray burst data were taken from GOES satellite observations in the wavelength band of 0.05–0.4 nm. It was found that (i) the amplitude of SPAs in different paths varies from 10 to 282∘, and (ii) the correlation between the amplitudes of SPAs, the lengths of paths and the signal frequency is weak. The change in effective height of reflection due to lowering of the reflecting layer during the flares was found to be about 12 km for the first event and about 9 km for the second event. Spectral analysis of the X-ray and LF data, filtered in the range between 5 s and 16 min, showed that the LF signal spectra are very similar to X-ray spectra. Maxima of both X-ray and LF spectra are in 2–16 min interval.

The analyzed amplitude and phase variations in electromagnetic VLF and LF signals at 20–45 kHz, received in Moscow, Graz (Austria), and Sheffield (UK) during the total solar eclipse of March 20, 2015, are considered. The 22 analyzed paths have lengths of 200—6100 km, are differently oriented, and cross 40–100% occultation regions. Fifteen paths crossed the region where the occultation varied from 40 to 90%. Solar eclipse effects were found only on one of these paths in the signal phase (–50°). Four long paths crossed the 90–100% occultation region, and signal amplitude and phase anomalies were detected for all four paths. Negative phase anomalies varied from–75° to–90°, and the amplitude anomalies were both positive and negative and were not larger than 5 dB. It was shown that the effective height of the ionosphere varied from 6.5 to 11 km during the eclipse.

The disturbances of very low and low frequency signals in the lower mid-latitude ionosphere caused by magnetic storms, proton bursts and relativistic electron fluxes are investigated on the basis of VLF–LF measurements obtained in the Far East and European networks. We have found that magnetic storm (−150 < Dst < −100 nT) influence is not strong on variations of VLF–LF signals. The anomalies with negative amplitude were registered during the main and recovery phases for several magnetic storms (mainly for three northernmost paths). The correlation between VLF–LF signals and geomagnetic activity is rather weak even for these paths (≈ 12–18%). Also, the correlation between magnetic activity and VLF signal variations recorded onboard the DEMETER satellite is not found. The significant influence of outer-zone particles (energetic particle sensor on board/Geostationary Operational Environmental Satellite (GOES) measurements) on the VLF–LF signal variations is found for almost half of the sub-ionospheric paths.

International audience; The paper presents the results of a joint analysis of ground-based and satellite observations of very low-frequency and low-frequency (VLF/LF) signals during periods of strong seismic activity in the region of Kuril Islands and Japan in 2004–2010. Ground and satellite data was processed using a method based on the difference between the real signal in nighttime and that of a model. The results of the analysis show a good correlation between ground-based and satellite data for several cases of strong (M ⩾ 6.8) earthquakes.

The subionospheric data from a Japanese very low frequency/low frequency (VLF/LF) receiving station at Moshiri, Hokkaido, are used to detect the response of the lower ionosphere to the tsunami triggered by the 2011 Tohoku earthquake. Disturbances in the phase and amplitude of VLF signals propagating from the transmitter in Hawaiian Islands are observed during the tsunami wave passage, and these effects in the ionosphere are compared to the in situ sea-level global positioning system (GPS) measurements near Japan. The frequency of the maximum spectral amplitude both for the VLF and GPS data is found to be in the range of periods of 8–50 min, which is likely to correspond to the period of the internal gravity waves generated by the tsunami.

We report on specific fluctuations in phase and amplitude of VLF signals that correlate both spatially and temporally with the passage of the tsunamis recorded by the Deep-ocean Assessments and Reporting of Tsunamis bottom pressure stations. Measurements from the VLF/LF receiver sited in Petropavlovsk-Kamchatsky and sensor buoys placed throughout the Pacific Ocean at great distances (Hawaii and Japan) from the epicenter are consistent with the hypothesis that the ocean tsunami following the Chile earthquake on 27 February 2010 radiated internal gravity waves which propagated through the lower ionosphere.

In this paper we develop a new method for the analysis of excitation and propagation of planetary electromagnetic waves (PEMW) in the ionosphere of the Earth. The nonlinear system of equations for PEMW, valid for any height, from D to F regions, including intermediate altitudes between D and E and between E and F regions, is derived. In particular, we have found the system of nonlinear one-fluid MHD equations in the β-plane approximation valid for the ionospheric F region (Aburjania et al., 2003a, 2005). The series expansion in a "small" (relative to the local geomagnetic field) non-stationary magnetic field has been applied only at the last step of the derivation of the equations. The small mechanical vertical displacement of the media is taken into account. We have shown that obtained equations can be reduced to the well-known system with Larichev–Reznik vortex solution in the equatorial region (see e.g. Aburjania et al., 2002). The excitation of planetary electromagnetic waves by different initial perturbations has been investigated numerically. Some means for the PEMW detection and data processing are discussed.

Data from the Pacific network of VLF receivers have been used to study the response of the lower ionosphere to the January 2011 Mt. Kirishima (South Japan) volcanic eruption. A major explosive eruption occurred on January 27, which was preceded by several small eruptions. Perturbations of nighttime subionospheric VLF signals have been detected on the day of the first small eruption on January 18 (UT) with the maximum observed about 1.5 h after the eruption. The nighttime signal remained disturbed during the subsequent pre-eruptive and eruptive activity of Mt. Kirishima. The daytime perturbations were not observed. The frequency of the maximum spectral amplitude was found to be in the range of periods of 6–30 min, which corresponds to the periods of internal gravity waves. These results suggest that the observed VLF ionospheric effects can possibly be produced by the penetration of gravity waves caused by the volcanic activity into the ionosphere.

Two methods for diagnosing ionospheric disturbances using VLF signals received on the DEMETER satellite are considered in connection with two cases of high seismic activity: the method for analyzing a signal reception zone in order to find large-scale spatial variations and the residual signal method used to reveal time variations. Based on an analysis of the VLF signal reception zone, performed in connection with the catastrophic earthquake near Sumatra, it has been found out that the signal amplitude decreased during a month before the earthquake. As a result, it has been concluded that the zone of ionospheric disturbances extends for several thousand kilometers. It has been indicated that the residual signal method showed good agreement between the ground and satellite data when high seismic activity near Japan was analyzed.

Data from the European network of very low/ low frequency (VLF/LF) receivers has been used to study the response of the lower ionosphere to the earthquake of magnitude 5.5 in Vrancea area on November 22, 2014. Negative amplitude anomalies have been observed during 3 days before the earthquake and two days after, on the LF (45.9 kHz) signal passed above the seismic area. No perturbations have been found for the same signal in control paths during this period. Other possible influences both from above and below which can produce perturbations in the ionosphere have been taken into consideration.

The analysis of the VLF signals radiated by ground transmitters and received on board of the French DEMETER satellite, reveals a drop of the signals (scattering spot) connected with the occurrence of large earthquakes. The extension of the "scattering spots" zone is large enough (1000–5000 km) and, probably, it increases with the magnitude of the "relative" earthquake. A possible model to explain the phenomenology, based on the acoustic gravity waves and the ionosphere turbulence, is proposed. The method of diagnostics applied to this study has the advantage to be a global one due to the world wide location of the powerful VLF transmitters and of the satellite reception. However, a specific disadvantage exists because the method requires rather a long time period of analysis due to the large longitudinal displacements among the successive satellite orbits. At the moment, at least, one month seems to be necessary.

Variations of the LF signal from observations in 6 receiving stations during 2000–2004 are analyzed. Five stations are located in Japan (Kochi, Moshiri, Kasugai, Shimizu and Chiba) and one station is in Russia (Petropavlovsk–Kamchatski). All the stations are equipped with the identical OMNIPAL receivers. JJY (40 kHz) transmitter is situated in Fukushima prefecture (Japan). Sensitivity of the phase and amplitude of the LF signal to ionosphere and atmosphere perturbations associated with variations of the magnetic and meteorological conditions and solar activity is studied. Correlation of phase and amplitude variations of the LF signal with D st -index, outer-zone particles (protons and electrons) with high-pitch angle is established for all the stations. It is found that LF anomalies are typical for the main phase of magnetic storm if the sudden commencement happens in the night time. The similarity of phase and amplitude anomalies to the structure of the simultaneously observed Pi3 geomagnetic pulsations is found for all the stations. Correlation between variations of atmosphere pressure and amplitude of the LF signal is established. Correlation between variations of humidity for the stations located on the coast of Japan and amplitude of the LF signal is found.

The variations of phase and amplitude of LF signal propagating from the transmitter JJY (40 kHz, Fukushima prefecture, Japan) are analyzed. The receiver is situated in Petropavlovsk-Kamchatski (Russia). The time interval from July1, 2004 till January 24, 2005 is included in the examination. This interval is characterized by quiet seismic conditions up to the beginning of November, but rather strong seismic activity occurs in November and December not far from Hokkaido (Japan) and in the region of northern Kuril Islands. Seismic activity in a zone of sensitivity of our wave paths is determined by three series of earthquakes with M = 5.6–7.1. Nighttime “bay-like” phase and amplitude anomalies of the LF signal begin some days before the first earthquake and continue till the last earthquake of this series ( M = 6.2, November 11, 2004). Most strong and obvious anomalies appear in the amplitude of the LF signal. After that, during about one fortnight a period of seismic calm is observed. Then there is a following anomaly of the LF signal before two earthquakes ( M = 7.1, 28.11.04 and M = 6.8, 6.12.04). The next LF anomaly is observed before the earthquakes on December, 18 and 21 ( M = 5.8 and M = 5.6). During the whole period of seismic activity a significant shift in terminator times is also evident. The spectrum of LF seismo-induced anomalies shows a clear increase for the periods of 10–25 min.

We present a general concept of mechanisms of preseismic phenomena in the atmosphere and ionosphere. After short review of observational results we conclude: 1. Upward migration of fluid substrate matter (bubble) can lead to ousting of the hot water/gas near the ground surface and cause an earthquake (EQ) itself in the strength-weakened area; 2. Thus, time and place of the bubble appearance could be random values, but EQ, geochemistry anomaly and foreshocks (seismic, SA and ULF electromagnetic ones) are casually connected; 3. Atmospheric perturbation of temperature and density could follow preseismic hot water/gas release resulting in generation of atmospheric gravity waves (AGW) with periods in a range of 6–60min; 4. Seismo-induced AGW could lead to modification of the ionospheric turbulence and to the change of over-horizon radio-wave propagation in the atmosphere, perturbation of LF waves in the lower ionosphere and ULF emission depression at the ground.

Very low and low frequency (VLF/LF) data recorded in the Far Eastern stations Petropavlovsk-Kamchatsky (158.92° E, 53.15° N), Yuzhno-Sakhalinsk (142.75° E, 46.95° N) and Yuzhno-Kurilsk (145.861° E, 44.03° N) are investigated to study the meteorological effects in the lower ionosphere. The results demonstrate the sensitivity of the VLF/LF signals to the variations of atmospheric pressure, humidity, wind velocity and temperature, and the VLF/LF record at the station of Yuzhno-Kurilsk is found to be most sensitive to those variations of atmospheric parameters. The region under consideration is characterized by high winter cyclonic activity in mid-latitudes and strong summer and autumn typhoon activity in low latitudes. VLF/LF signal variations during eight tropical cyclones (TCs) with different intensity are considered. Negative nighttime anomalies in the signal amplitude that are most probably caused by TC activity are found for six events. Those anomalies are observed during 1–2 days when TCs move inside the sensitivity zones of the subionospheric paths. Perturbations of the VLF signal observed during two TCs can be caused by both the TC influence and seismic activity, but no correlation between TC intensity and magnitude of the signal anomalies is found. Spectral analysis of the typhoon-induced disturbed signals revealed the fluctuations with time periods in the range of 7–16 and 15–55 min that corresponds to the range of internal gravity waves periods.

We examined changes in the spectral composition of the low frequency (LF) subionospheric signals from the NRK transmitter (37.5 kHz) in Iceland that were received in Bari (Italy) relative to the earthquake that occurred in L’Aquila on April 6, 2009. In our previous studies, we have reported the occurrence of preseismic night-time anomalies using observations from three receivers located in Bari, Graz (Austria) and Moscow (Russia). The strongest anomalies in the signals were observed in the NRK-Bari propagation path during the period 5-6 days before the L’Aquila earthquake, as well as during the series of aftershocks. During this period, similar very low frequency (VLF)/LF amplitude anomalies were also observed along several other propagation paths that crossed the L’Aquila seismogenic zone. Spectral analysis of the LF signals filtered in the frequency range 0.28 mHz to 15 mHz shows differences in the spectra for seismo-disturbed days when compared to those for either quiet or geomagnetically disturbed days. These spectral anomalies, which are only observed in the propagation path between NRK and Bari, contain signals with periods of about 10 min to 20 min. These periodic signals are absent both in the spectra of the undisturbed signals for the control paths, and in the spectra of the signals received during geomagnetic storms. The same changes in the spectral composition were observed in the analysis of LF (40 kHz) signals from the JJY transmitter in Japan that were received in Petropavlovsk-Kamchatsky (Russia) during the occurrence of three strong earthquakes with M ≥7.0. The results of this study support the theoretical prediction that the possible mechanism for energy penetration from the origin of an earthquake through the atmosphere and into the ionosphere is based on the excitation and upward propagation of internal gravity waves.

In this paper we present results obtained from more than two years continuous very low frequency (VLF) measurements between a network of VLF transmitters and receivers. The focus is on VLF amplitude and phase variations in the Earth lithosphere-ionosphere cavity with the scientific objective to characterise long term trends possibly related to planetary waves and climatological factors. After considering the nominal diurnal and seasonal behaviour of the individual paths in the VLF waveguide, we distinguish between natural and artifical disturbances on different temporal and spatial scales. We conclude that it is too early to speak of a clear trend, but VLF investigations together with complementary measurements can be a useful tool in long term environmental and climatological studies on medium spatial scales.

By using the network observation of subionospheric VLF/LF signals in Japan and in Russia, we have found a significant ionospheric perturbation prior to the recent 2011 March 11 Japan earthquake (EQ) in the off-sea of the Tohoku area, which was an exceptionally huge plate-type EQ. A remarkable anomaly (with decrease in the nighttime amplitude and also with enhancement in dispersion) has been detected on March 5 and 6 on the propagation path from the NLK transmitter (Seattle, USA) to Chofu (together with Kochi and Kasugai), and also we have observed the corresponding VLF anomaly during a prolonged period of March 1–6, with minima in the nighttime amplitude on March 3 and 4 on the path from JJI transmitter (Miyazaki, Kyushu) to Kamchatka, Russia. Используя наблюдения распространения СДВ/ДВ-радиоволн над Тихим океаном на японской и российской сети станций, удалось обнаружить значительное возмущение ионосферы, предшествовавшее последнему мощному землетрясению в Японии 11.03.2011 г. Эпицентр землетрясения находился в море, в области Тохоку, а само событие относится к исключительно мощным землетрясениям, связанным с перемещением тектонических плит. Явно выраженная аномалия (уменьшение ночной амплитуды сигнала при увеличении ее дисперсии) была обнаружена 5 и 6 марта на трассе распространения от передатчика NLK (Сиэтл, США) к наблюдателю в Чофу, Япония (аналогичные явления – на трассах распространения в Кочи и Кацугаи). Аналогичная длительная аномалия в СДВ-распространении регистрировалась с 1 по 6 марта с минимальной ночной амплитудой 3 и 4 марта на трассе от передатчика JJI (Миязаки, Кюсю) до Камчатки, Россия. Використовуючи спостереження поширення СДВ/ДВ-радіохвиль над Тихим океаном на японській і російській мережі станцій, вдалося виявити значне збурення іоносфери, що сталося перед останнім потужним землетрусом у Японії 11.03.2011 р. Епіцентр землетрусу знаходився в морі, в області Тохоку, а сама подія відноситься до виключно потужних землетрусів, пов’язаних з переміщенням тектонічних плит. Явно виражена аномалія (зменшення нічної амплітуди сигналу при збільшенні її дисперсії) було виявлено 5 та 6 березня на трасі від передавача NLK (Сіетл, США) до спостерігача в Чофу, Японія (аналогічні явища – на трасах поширення до Кочі й Кацугаї). Аналогічну тривалу аномалію в СДВ-поширенні реєстрували з 1 по 6 березня з мінімальною нічною амплітудою 3 і 4 березня на трасі від передавача JJI (Міязакі, Кюсю) до Камчатки, Росія.

We studied six different parameters derived from three different radiophysical techniques used to monitor area over the central Italy around the time of L'Aquila M6.3 earthquake of 6 April 2009. Namely VLF monitoring, ground based ionospheric sounding and GPS vertical TEC. The parameters derived from these data are: variations of amplitude of VLF signal on different propagation passes, critical frequency foF2 variations, cross-correlation coefficient for different pairs of ground based ionosondes, vertical GPS TEC time series, local variability index for the network of GPS receivers, and GPS TEC map. High congruence of techniques in time and space is found.

In our earlier papers we have found the effect of VLF transmitter signal depression over epicenters of the large earthquakes from observation on the French DEMETER satellite that can be considered as new method of global diagnostics of seismic influence on the ionosphere. At present paper we investigate a possibility VLF signal-ionospheric turbulence interaction using additional characteristic of VLF signal-spectrum broadening. This characteristic is important for estimation of the interaction type: linear or nonlinear scattering. Our main results are the following: – There are two zones of increased spectrum broadening, which are centered near magnetic latitudes Φ=±10° and Φ=±40°. Basing on the previous case study research and ground ionosonde registrations, probably it is evidence of nonlinear (active) scattering of VLF signal on the ionospheric turbulence. However occurrence rate of spectrum broadening in the middle-latitude area is higher than in the near-equatorial zone (~15–20% in comparison with ~100% in former area) that is probably coincides with the rate of ionospheric turbulence. – From two years statistics of observation in the selected 3 low-latitude regions and 1 middle-latitude region inside reception area of VLF signal from NWC transmitter we find a correlation of spectrum broadening neither with ion-cyclotron noise (f=150–500 Hz), which possibly means poor representation of the turbulence by the noise due to its mixture with natural ELF emission (which correlates with whistler), nor with magnetic storm activity. – We find rather evident correlation of ion-cyclotron frequency noise, VLF signal depression and weak correlation of spectrum broadening with seismicity in the middle-latitude region over Japan. But in the low-latitude regions we do not find such a correlation. Statistical decrease of VLF signal supports our previous case study results. However rather weak spectrum broadening-seismicity statistical correlation means probably that passive scattering prevails upon nonlinear (active) one.

The results of the monitoring of three VLF/LF signals collected in Petropavlovsk station (Kamchatka, Russia) and one VLF signal collected on board of the DEMETER French satellite are presented. Two periods of the seismic activity occurred in the Japan-Kamchatka area during November–December 2004 and July–September 2005 were investigated and the earthquakes with M≥6.0 in the Japan-Kamchatka area, located inside one or more of the third Fresnel zones of the three radio paths were considered. The ground data were analysed using residual signal of phase dP or of amplitude dA, defined as the difference between the signal and the average of few quiet days (±5 days) immediately preceding or following the current day. Also the satellite data were processed by a method based on the difference between the real signal and the reference one, but in order to obtain this last signal it was necessary to construct previously a model of the signal distribution over the selected area. The method consists: (a) in averaging all the data available in the considered region over a period characterized by low level seismicity, regardless of the global disturbances, in particular, of the magnetic activity; (b) in computing a polynomial expression for the surface as a function of the longitude and the latitude. The model well describes the real data in condition of their completeness and in absence of magnetic storms or seismic forcing. In the quoted periods of seismic activity clear anomalies both in the ground and in satellite data were revealed. The influence of the geomagnetic activity cannot to be excluded, but the seismic forcing seems more probable.

By using network observation of subionospheric VLF (very low frequency)/LF (low frequency) signals in Japan and in Russia, we have found a significant ionospheric perturbation prior to the recent 2011 March 11 Japan earthquake (EQ) which occurred at sea proximate to the Tohoku area on the main island (Honshu) of Japan was an exceptionally huge plate-type EQ. A remarkable anomaly (with a decrease in the nighttime amplitude and also with enhancement in dispersion) was detected on March 5 and 6 along the propagation path from the NLK (Seattle, USA) transmitter to Chofu (together with Kochi and Kasugai). We also have observed the corresponding VLF anomaly during a prolonged period of March 1 - 6, with minima in the nighttime amplitude on March 3 and 4 along the path from JJI (Miyazaki, Kyushu) to Kamchatka, Russia. This ionospheric perturbation has been discussed extensively with respect to its reliability. (1) How abnormal is this VLF/LF propagation anomaly? (2) What was the temporal evolution of terminator times? (3) Were there any solar-terrestrial effects (especially the effect from geomagnetic storms) on the VLF/LF propagation anomaly? (4) The effect of any other EQs and foreshock activities on the VLF/LF anomaly? (5) Were there any correlations with other related phenomena? Finally, (6) are there any other examples of a VLF/LF propagation anomaly for oceanic EQs? We then compared the temporal properties of ionospheric perturbations for this EQ with those of a huge number of inland EQs and compared the corresponding spatial scale with the former result of the same oceanic 2004 Sumatra EQ with nearly the same magnitude. Finally, the generation mechanism of those seismo-ionospheric perturbations is briefly discussed.