A time-frequency analysis of the thermal radiation background anomalies caused by large earthquakes: A case study of the Wenchuan 8.0 earthquake

Strong earthquakes have an impact on the regional thermal radiation background, which has been both observed and confirmed. This effect produces anomalies in the thermal radiation background (TRBA) and increases the difficulty of extracting a thermal radiation anomaly (TRA) that is associated with an earthquake occurring during the same time period. The extraction and identification of such anomalies has been ignored by previous studies. In this study, we investigate the time-frequency analysis (TFA) method, together with the wavelet filtering of the Daubechies method and the relative power spectrum analysis of the Fourier Welch method to extract and analyse the TRBA caused by the 2008 Wenchuan Ms 8.0 earthquake and the TRA of the 2013 Minxian-Zhangxian Ms 6.6 earthquake using data concerning the brightness temperature of a black body (TBB) from the Fengyun-2 series of geostationary meteorological satellites developed by China. The result showed that this method can effectively extract and analyse the TRBA caused by the Wenchuan earthquake and the TRA of the Minxian-Zhangxian earthquake form a complex background environment. Furthermore, we discussed the impact of the earthquake on the TRBA and segmented the process. The impact is mainly reflected by three aspects; the characteristic period of the TRBA changes, the TRBA occurs at the same time every year, which is identical to the time at which the earthquake anomaly occurred, and the impact process is in stages. We also summarized the correlation between the characteristic parameters of a TRA and the regional thermal radiation background, geography, and climatic factors.     

Crustal deformation characteristics of Sichuan-Yunnan region in China on the constraint of multi-periods of GPS velocity fields

In order to obtain deformation parameters in each block of Sichuan-Yunnan Region (SYG) in China by stages and establish a dynamic model about the variation of the strain rate fields and the surface expansion in this area, we taken the Global Positioning System (GPS) sites velocity in the region as constrained condition and taken advantage of the block strain calculation model based on spherical surface. We also analyzed the deformation of the active blocks in the whole SYG before and after the Wenchuan earthquake, and analyzed the deformation of active blocks near the epicenter of the Wenchuan earthquake in detail. The results show that, (1) Under the effects of the carving from India plate and the crimping from the potential energy of Tibetan Plateau for a long time, there is a certain periodicity in crustal deformation in SYG. And the period change and the earthquake occurrence have a good agreement. (2) The differences in GPS velocity fields relative Eurasian reference frame shows that the Wenchuan earthquake and the Ya'an earthquake mainly affect the crustal movement in the central and southern part of SYG, and the average velocity difference is about 4–8?mm/a for the Wenchuan earthquake and 2–4?mm/a for the Ya'an earthquake. (3) For the Wenchuan earthquake, the average strain changed from 10 to 20 nanostrian/a before earthquake to 40–50 nanostrian/a after the earthquake, but before and after the Ya’an earthquake, the strain value increased from about 15 nanostrian/a to about 30 nanostrian/a. (4) The Wenchuan earthquake has changed the strain parameter of each active block more or less. Especially, the Longmen block and Chengdu block near the epicenter. The research provides fundamental material for the study of the dynamic mechanism of the push extrusion from the north-east of the India plate and the crimp from Qinghai Tibet Plateau, and it also provides support for the study of crustal stress variation and earthquake prediction in Sichuan Yunnan region.     

Statistical analysis of the ionospheric ion density recorded by DEMETER in the epicenter areas of earthquakes as well as in their magnetically conjugate point areas

Results of a statistical variation of total ion density observed in the vicinity of epicenters as well as around magnetically conjugated points of earthquakes are presented in this paper. Two data sets are used: the ion density measured by DEMETER during about 6.5?years and the list of strong earthquakes (M_W?≥?4.8) occurring globally during this period (14,764 earthquakes in total). First of all, ionospheric perturbations with 23–120?s observation time corresponding to spatial scales of 160–840?km are automatically detected by a software (64,287 anomalies in total). Second, it is checked if a perturbation could be associated either with the epicenter of an earthquake or with its magnetically conjugated point (distance?<?1500?km and time?<?15?days before the earthquake). The index Kp?<?3 is also considered in order to reduce the effect of the geomagnetic activity on the ionosphere during this period. The results show that it is possible to detect variations of the ionospheric parameters above the epicenter areas as well as above their conjugated points. About one third of the earthquakes are detected with ionospheric influence on both sides of the Earth. There is a trend showing that the perturbation length increases as the magnitude of the detected EQs but it is more obvious for large magnitude. The probability that a perturbation appears is higher on the day of the earthquake and then gradually decreases when the time before the earthquake increases. The spatial distribution of perturbations shows that the probability of perturbations appearing southeast of the epicenter before an earthquake is a little bit higher and that there is an obvious trend because perturbations appear west of the conjugated point of an earthquake.     

Thermal and TEC anomalies detection using an intelligent hybrid system around the time of the Saravan,Iran, (Mw = 7.7) earthquake of 16 April 2013

A powerful earthquake of M_w = 7.7 struck the Saravan region (28.107° N, 62.053° E) in Iran on 16 April 2013. Up to now nomination of an automated anomaly detection method in a non linear time series of earthquake precursor has been an attractive and challenging task. Artificial Neural Network (ANN) and Particle Swarm Optimization (PSO) have revealed strong potentials in accurate time series prediction. This paper presents the first study of an integration of ANN and PSO method in the research of earthquake precursors to detect the unusual variations of the thermal and total electron content (TEC) seismo-ionospheric anomalies induced by the strong earthquake of Saravan. In this study, to overcome the stagnation in local minimum during the ANN training, PSO as an optimization method is used instead of traditional algorithms for training the ANN method. The proposed hybrid method detected a considerable number of anomalies 4 and 8 days preceding the earthquake. Since, in this case study, ionospheric TEC anomalies induced by seismic activity is confused with background fluctuations due to solar activity, a multi-resolution time series processing technique based on wavelet transform has been applied on TEC signal variations. In view of the fact that the accordance in the final results deduced from some robust methods is a convincing indication for the efficiency of the method, therefore the detected thermal and TEC anomalies using the ANN + PSO method were compared to the results with regard to the observed anomalies by implementing the mean, median, Wavelet, Kalman filter, Auto-Regressive Integrated Moving Average (ARIMA), Support Vector Machine (SVM) and Genetic Algorithm (GA) methods. The results indicate that the ANN + PSO method is quite promising and deserves serious attention as a new tool for thermal and TEC seismo anomalies detection.     

Genetic algorithm for TEC seismo-ionospheric anomalies detection around the time of the Solomon (Mw = 8.0) earthquake of 06 February 2013

On 6 February 2013, at 12:12:27 local time (01:12:27 UTC) a seismic event registering M_w 8.0 struck the Solomon Islands, located at the boundaries of the Australian and Pacific tectonic plates. Time series prediction is an important and widely interesting topic in the research of earthquake precursors. This paper describes a new computational intelligence approach to detect the unusual variations of the total electron content (TEC) seismo-ionospheric anomalies induced by the powerful Solomon earthquake using genetic algorithm (GA). The GA detected a considerable number of anomalous occurrences on earthquake day and also 7 and 8 days prior to the earthquake in a period of high geomagnetic activities. In this study, also the detected TEC anomalies using the proposed method are compared to the results dealing with the observed TEC anomalies by applying the mean, median, wavelet, Kalman filter, ARIMA, neural network and support vector machine methods. The accordance in the final results of all eight methods is a convincing indication for the efficiency of the GA method. It indicates that GA can be an appropriate non-parametric tool for anomaly detection in a non linear time series showing the seismo-ionospheric precursors variations.     

Removing Solar Radiation Based on the Empirical Mode Decomposition Method for Seismo-Ionospheric Anomaly Before the M9.0 Tohoku Earthquake

The ionospheric effect from solar activity can be seen as the background in the process of detecting the ionospheric precursor prior to strong earthquakes. The ionospheric variation induced by the forthcoming earthquake can be covered by the strong solar background during the period of high solar activity. The issue of how to remove the ionospheric effect from solar radiation is of outstanding significance. In this paper, a method of Empirical Mode Decomposition (EMD) is used to eliminate the solar background. As a case study, the global ionospheric map TEC before the M9.0 Tohoku earthquake on 11 March 2011 is analyzed. After the effect of solar radiation is removed using the EMD method, the precursor of the imminent earthquake is more obvious. The ionospheric anomaly had a local character and only appeared close to the earthquake epicenter while the useful signals were covered by the solar radiation background with traditional method, which implies that the EMD method is effective in eliminating solar radiation background.      

An Adaptive Network-based Fuzzy Inference System for the detection of thermal and TEC anomalies around the time of the Varzeghan,Iran, (Mw = 6.4) earthquake of 11 August 2012

Anomaly detection is extremely important for forecasting the date, location and magnitude of an impending earthquake. In this paper, an Adaptive Network-based Fuzzy Inference System (ANFIS) has been proposed to detect the thermal and Total Electron Content (TEC) anomalies around the time of the Varzeghan, Iran, (M_w = 6.4) earthquake jolted in 11 August 2012 NW Iran. ANFIS is the famous hybrid neuro-fuzzy network for modeling the non-linear complex systems. In this study, also the detected thermal and TEC anomalies using the proposed method are compared to the results dealing with the observed anomalies by applying the classical and intelligent methods including Interquartile, Auto-Regressive Integrated Moving Average (ARIMA), Artificial Neural Network (ANN) and Support Vector Machine (SVM) methods. The duration of the dataset which is comprised from Aqua-MODIS Land Surface Temperature (LST) night-time snapshot images and also Global Ionospheric Maps (GIM), is 62 days. It can be shown that, if the difference between the predicted value using the ANFIS method and the observed value, exceeds the pre-defined threshold value, then the observed precursor value in the absence of non seismic effective parameters could be regarded as precursory anomaly. For two precursors of LST and TEC, the ANFIS method shows very good agreement with the other implemented classical and intelligent methods and this indicates that ANFIS is capable of detecting earthquake anomalies. The applied methods detected anomalous occurrences 1 and 2 days before the earthquake. This paper indicates that the detection of the thermal and TEC anomalies derive their credibility from the overall efficiencies and potentialities of the five integrated methods.     

应用空间技术 GPS 监测地震

空间定位技术GPS 可以达到很高的相对定位精度(10~(-6)～10~(-8)),是监测地壳运动、进行地震预报的一种新手段。中国华北是一个地震活跃地区。应用GPS 技术在华北布设了一个地震监测网,经过精心观测、优化设计和严密平差,获得了高精度的第一期成果。华北GPS 网的建立将为大震预测和地球科学各领域的研究提供重要的科学依据。     

Is there an abnormal enhancement of atmospheric aerosol before the 2008 Wenchuan earthquake?

There are extensive reports of ionospheric disturbances before the great 2008 Wenchuan earthquake, which are possibly explained by seismogenic electric field hypotheses linked with the aerosols injected in atmosphere. This paper attempts to investigate the possible change of atmospheric aerosol optical depth (AOD) associated with this earthquake by using MODIS data from both Terra and Aqua satellites. The result shows a clear enhancement of AOD along the Longmenshan faults 7 days before the quake, which is 1 day and 4 days earlier than the reported negative and positive ionospheric disturbances, respectively, and is 1 day earlier than or quasi-synchronism with other reported atmospheric anomalies including air temperature, outgoing longwave radiation and relative humidity. Particularly, the spatial distribution of AOD enhancement is very local and it is correlated well with the active faults and surface ruptures. We suggest that this unique enhancement could be associated with the Lithosphere–Atmosphere–Ionosphere coupling process during the preparation of the Wenchuan earthquake.     

Change in refractivity of the atmosphere and large variation in TEC associated with some earthquakes,observed from GPS receiver

The present study reports the analysis of GPS based TEC for our station Surat (21.16°N, 72.78°E) located at the northern crest of equatorial anomaly region in India at times close to some earthquake events (M ? 5) during the year 2009 in India and its neighbouring regions. The TEC data used in the study are obtained from GPS Ionospheric Scintillation and TEC Monitoring (GISTM) system. The TEC data has been analysed corresponding to 11 earthquakes in low solar activity period and quiet geomagnetic condition. We found that, out of 11 cases of earthquakes (M > 5) there were seven cases in which enhancement in TEC occurred on earthquake day and in other four cases there was depletion in TEC on earthquake day. The variation in refractivity prior to earthquake was significant for the cases in which the epicentre lied within a distance of 600 km from the receiving station. By looking into the features on temporal enhancement and depletion of TEC a prediction was made 3–2 days prior to an earthquake (on 28 October 2009 in Bhuj – India). The paper includes a brief discussion on the method of potentially identifying an impending earthquake from ionospheric data.     

Fault rupture model of the 2008 Dangxiong (Tibet,China) Mw 6.3 earthquake from Envisat and ALOS data

On October 6, 2008, an Mw 6.3 earthquake occurred in Dangxiong county, southern Tibetan Plateau. In this study, Synthetic Aperture Radar (SAR) images from Envisat ASAR C-band descending Track 176 and ALOS PALSAR L-band ascending Track 500 are processed to generate the coseismic deformation caused by the earthquake. To estimate the source model, a downhill simplex non-linear inversion method is used to determine the fault rupture geometry, and an automatic fault discretization technique is employed to divide the fault plane to construct the optimal slip model, in which the uncertainties of the fault parameters are assessed by a Monte Carlo method. The inversion results show that the earthquake strikes almost south–north and has a normal faulting focal mechanism with rake angle and slip of −111.7° and 1.33 m, respectively. Peak slip of 2.15 m is located at a depth of 7.5 km. The estimated geodetic moment is 4.06 × 10¹⁸ N m (Mw 6.37), 71.2% of which is released in the depth range 4.5–11 km. The slip model suggests that coseismic slip also takes place at some fault patches near the earth’s surface and postseismic afterslip occurs below the coseismic rupture area after the earthquake.     

Ionospheric precursors observed at low latitudes around the time of koyna earthquake

On December 11, 1967 at 05:21 LT, an immense earthquake of magnitude 6.7 struck Koyna, the Indian province of Maharashtra. Its epicenter was located at geographic latitude 17.37°N and longitude 73.75°E with depth of about 3 km. Ground based measurements show variation in the critical frequency of ionospheric F2 layer (foF2) before and after the shock. In the present study the behavior of F2-region of ionosphere has been examined over the equatorial and low latitudinal region ionosphere during the month of December 1967 around the time of Koyna earthquake. For this purpose, the ionospheric data collected with the help of ground-based ionosondes installed at Hyderabad (located close to the earthquake epicenter) Ahmedabad, Trichirapulli, Kodaikanal and Trivendrum have been utilized. The upper and lower bound of Interquartile range (IRQ) are constructed to monitor the variations in foF2 other than day-to-day and diurnal pattern for finding the seismo-ionospheric precursors. Some anomalous electron density variations are observed between post midnight hours to local pre-noon hours at each station. These anomalies are strongly time dependent and appeared a couple of days before the main shock. The period considered in this study comes under the quiet geomagnetic conditions. Hence, the observed anomalies (which are more than the usual day-to-day variability) over all stations are likely to be associated with this imminent earthquake. The possible mechanism to explain these anomalies is the effect of seismogenic electric field generated just above the surface of earth within the earthquake preparation zone well before the earthquake due to emission of radioactive particles and then propagated upward, which perturbs the F-region ionosphere.     

GPS total electron content and surface latent heat flux variations before the 11 March 2011 M9.0 Sendai earthquake

The present paper describes the variations of the GPS total electron content (TEC) from the International GNSS service network and surface latent heat flux (SLHF) from the Scientific Computing Division of the National Center for Atmospheric Research (NCAR) before the 11 March 2011 M9.0 Sendai earthquake, respectively. The analysis shows pronounced enhancements in the GPS TEC and SLHF a few days prior to the earthquake event. The maximum increase in the GPS TEC was about 30 TECu with an extended spatial distribution on a geomagnetically quiet day (Dst ? −20 nT, between two moderate geomagnetic storms), 8 March, 3 days prior to the earthquake. This giant positive disturbance was possibly associated with the impending disastrous earthquake and contributed from the enhanced solar radiation. Moreover, there were several anomalous regions of SLHF on the global map, but an area of enhanced SLHF very close to the epicenter. The purpose of this paper is to report the existence of the changes in surface and ionosphere, and show the potential application of multi-source data to identify seismic precursors.     

Space technologies for short-term earthquake warning

Recent theoretical and experimental studies explicitly demonstrated the ability of space technologies to identify and monitor the specific variations at near-earth space plasma, atmosphere and ground surface associated with approaching severe earthquakes (named as earthquake precursors) which appear several days (from 1 to 5) before the seismic shock over the seismically active areas. Several countries and private companies are in the stage of preparation (or already launched) the dedicated spacecrafts for monitoring of the earthquake precursors from space and for short-term earthquake prediction. The present paper intends to outline the optimal algorithm for creation of the space-borne system for the earthquake precursors monitoring and for short-term earthquake prediction. It takes into account the following:

• 1.Selection of the precursors in the terms of priority, considering their statistical and physical parameters.
• 2.Configuration of the spacecraft payload.
• 3.Configuration of the satellite constellation (orbit selection, satellite distribution, operation schedule).
• 4.Different options of the satellite systems (cheap microsatellite or comprehensive multisatellite constellation).

Taking into account that the most promising are the ionospheric precursors of earthquakes, the special attention is devoted to the radiophysical techniques of the ionosphere monitoring. The advantages and disadvantages of such technologies as vertical sounding, in-situ probes, ionosphere tomography, GPS TEC and GPS MET technologies are considered.     

Seismo-ionospheric anomalies in total electron content of the GIM and electron density of DEMETER before the 27 February 2010 M8.8 Chile earthquake

In this paper we examine the pre-earthquake ionospheric anomalies by the total electron content (TEC) extracted from GIM (global ionospheric map) and the electron density (Ne) observed by the DEMETER (Detection of Electro-Magnetic Emissions Transmitted from Earthquake Regions) satellite during the 2010 M8.8 Chile earthquake. Temporal variations show the nighttime TEC and Ne simultaneously increase 9–19 days before the earthquake. A cross-comparison of data recorded during the period of 1 February to 3 March in 2006–2010 confirms the above temporal anomalies specifically appear in 2010. The spatial analyses show that the anomalies tend to appear over the epicenter.     

Multi precursors analysis associated with the powerful Ecuador (MW = 7.8) earthquake of 16 April 2016 using Swarm satellites data in conjunction with other multi-platform satellite and ground data

After DEMETER satellite mission (2004–2010), the launch of the Swarm satellites (Alpha (A), Bravo (B) and Charlie (C)) has created a new opportunity in the study of earthquake ionospheric precursors. Nowadays, there is no doubt that multi precursors analysis is a necessary phase to better understand the LAIC (Lithosphere Atmosphere Ionosphere Coupling) mechanism before large earthquakes. In this study, using absolute scalar magnetometer, vector field magnetometer and electric field instrument on board Swarm satellites, GPS (Global Positioning System) measurements, MODIS-Aqua satellite and ECMWF (European Centre for Medium-Range Weather Forecasts) data, the variations of the electron density and temperature, magnetic field, TEC (Total Electron Content), LST (Land Surface Temperature), AOD (Aerosol Optical Depth) and SKT (SKin Temperature) have been surveyed to find the potential seismic anomalies around the strong Ecuador (M_w = 7.8) earthquake of 16 April 2016. The four solar and geomagnetic indices: F10.7, D_st, K_p and a_p were investigated to distinguish whether the preliminary detected anomalies might be associated with the solar-geomagnetic activities instead of the seismo-ionospheric anomalies. The Swarm satellites (A, B and C) data analysis indicate the anomalies in time series of electron density variations on 7, 11 and 12 days before the event; the unusual variations in time series of electron temperature on 8 days preceding the earthquake; the analysis of the magnetic field scalar and vectors data show the considerable anomalies 52, 48, 23, 16, 11, 9 and 7 days before the main shock. A striking anomaly is detected in TEC variations on 1 day before earthquake at 9:00 UTC. The analysis of MODIS-Aqua night-time images shows that LST increase unusually on 11 days prior to main shock. In addition, the AOD variations obtained from MODIS measurements reach the maximum value on 10 days before the earthquake. The SKT around epicentral region presents anomalous higher value about 40 days before the earthquake. It should be noted that the different lead times of the observed anomalies could be acknowledged based on a reasonable LAIC earthquake mechanism. Our results emphasize that the Swarm satellites measurements play an undeniable role in progress the studies of the ionospheric precursors.     

Surface latent heat flux anomalies quasi-synchronous with ionospheric disturbances before the 2007 Pu’er earthquake in China

One of various mechanisms of pre-earthquake lithosphere–atmosphere–ionosphere coupling as possible explanation of the seismo-ionospheric effects is connected with the release of latent heat. Abnormal variations of ionospheric electromagnetic parameters possibly related to the 2007 Ms 6.4 Pu’er earthquake in China were reported. This paper attempts to examine whether there were abnormal changes of surface latent heat flux (SLHF) linked with this pre-earthquake ionospheric disturbances. The spatio-temporal statistical analyzes of multi-years SLHF data from USA NCEP/NCAR Reanalysis Project reveal that local SLHF enhancements appeared 11, 10 and 7 days before the Pu’er earthquake, respectively. As contrasted to the formerly reported local ionospheric Ne enhancement 9 days before the shock observed by DEMETER satellite, it is discovered that the SLHF and Ne anomalies are quasi-synchronous and have good spatial correspondence with the epicentre and the local active faults. This is valuable for understanding the seismogenic coupling processes and for recognizing earthquake anomaly with multiple parameters from integrated Earth observation system.     

Investigation of GPS-TEC measurements using ANN method indicating seismo-ionospheric anomalies around the time of the Chile (Mw = 8.2) earthquake of 01 April 2014

In 94 km NW of Iquique in Chile (19.610°S, 70.776°W) a powerful earthquake of M_w = 8.2 took place at 23:46:47 UTC (20:46:47 LT) on April 01, 2014. Using GPS-TEC (Total Electron Content) measurements, potential unusual variations around the time and location of the Chile earthquake have been detected based on the median and Artificial Neural Network (ANN) methods. The indices D_st, K_p, A_p and F10.7 were used to distinguish pre-earthquake anomalies from the other anomalies related to the solar-geomagnetic activities. Using the median method, striking anomalies in time series of TEC data are observed 4 days before the earthquake at 14:00 and 16:00 UTC. The ANN method detected a number of anomalies, 4 (02:00 and 16:00 UTC) and 13 (24:00 UTC) days preceding the earthquake. The results indicate that the ANN method due to its capability of non linear learning is quite promising and deserves serious attention as a robust predictor tool for seismo-ionospheric anomalies detection.     

Study of the 2013 Lushan M7.0 earthquake coseismic ionospheric disturbances

On April 20, 2013, an earthquake of M7.0 occurred in Lushan, Sichuan province, China. This paper investigates the coseismic ionospheric anomalies using GPS (Global Positioning System) data from 23 reference stations in Sichuan province that are a part of the Crustal Movement Observation Network of China (CMONOC). The recorded results show that a clear ionospheric anomaly occurred within 15 min after the earthquake near the epicenter, and the occurrence time of the anomalies recorded by various stations is related to the distance from the epicenter. The maximum anomaly is 0.25 TECu, with a 2 min duration and the distance of the recording station to the epicenter is 83 km. Acoustic waves generated by the crustal vertical movement during the earthquake propagate up to the height of the ionosphere lead to the ionospheric anomaly, and the propagation speed of the acoustic wave is calculated as 0.72 ± 0.04 km/s based on the propagation time and propagation distance, consistent with the average speed of sound waves within a 0–450 km atmospheric height.     

Investigation of precursory phenomena in the ionosphere,atmosphere and groundwater before large earthquakes of M > 6.5

Precursory phenomena in the ionosphere, atmosphere and groundwater before large earthquakes (M > 6.5) are extensively investigated toward the earthquake prediction. Upward tornado type seismic clouds occurred near the epicenter associated with strong LF-VLF radio noises from lightning discharges in the evening of January 9, 1995 [Yamada, T., Oike, K. On the increase of electromagnetic noises before and after the 1995 Hyogo-Ken Nanbu earthquake. In: Hayakawa M. (Ed.), Atmospheric and Ionospheric Electromagnetic Phenomena Associated with Earthquakes. TERRAPUB, Tokyo, pp. 417–427, 1999] and anomalous foEs increases up to 10 MHz were detected at Shigaraki, 90 km of the epicenter and at Kokubunji, 500 km east of the epicenter [Ondoh, T. Anomalous sporadic-E layers observed before M7.2 Hyogo-ken Nanbu earthquake; Terrestrial gas emanation model. Adv. Polar Upper Atmos. Res. 17, 96–108, 2003; Ondoh, T. Anomalous sporadic-E ionization before a great earthquake, Adv. Space Research 34, 1830–1835, 2004] associated with strong ELF noises from lightning discharges in the daytime on January 15, 1995 [Hata, M., Fujii, T., Takumi, I. EM precursor of large-scale earthquakes in Japan, in: Abstracts of International Workshop on Seismo Electromagnetics (IWSE 2005), Univ. Electro-Communications, Chofu, Tokyo, Japan, March 15–17, pp. 182–186, 2005] before the M7.2 Hyogoken–Nanbu earthquake of January 17, 1995. The anomalous foEs increases occurred at epicentral distances within 500 km that are the same as those of the terrestrial gas emanations along active faults before large earthquakes [King, C.-Y. Gas geochemistry applied to earthquake prediction: An overview. J. Geophys. Res. 91(B12), 12269–12281, 1986]. The anomalous foEs increases seem to be a seismic precursor because geomagnetic and solar conditions were very quiet all day on January 15,1995 and the normal foEs in Japanese winter is below 6 MHz. No significant pre-seismic geomagnetic field variation was detected at epicentral distance of 100 km before this earthquake [Ondoh, T., Hayakawa, M. Anomalous occurrence of sporadic-E layers before the Hyogoken–Nanbu earthquake, M7.2 of January 17, 1995. In: Hayakawa, M. (Ed.), Atmospheric and Ionospheric Electromagnetic Phenomena Associated with Earthquakes, TERRAPUB, Tokyo, pp. 629–639, 1999; Ondoh, T., Hayakawa, M. Seismo discharge model of anomalous sporadic E ionization before great earthquakes. In: Hayakawa, M., O.A. Molchanov, (Eds.), Seismo Electromagnetics: Lithosphere–Atmosphere–Ionosphere Couplings, TERRAPUB, Tokyo, pp. 385–390, 2002; Ondoh. T., Hayakawa, M. Synthetic study of precursory phenomena of the M7.2 Hyogo-ken Nanbu earthquake. Phys. Chem. Earth 31, 378–388, 2006]. The foF2 decrease and h’F increase occurred before the M7.8 Hokkaido Nansei-Oki earthquake of July 12,1993 in a geomagnetic quiet period [Ondoh, T. Ionospheric disturbances associated with great earthquake of Hokkaido southwest coast, Japan of July 12, 1993. Phys. Earth Planet. Interiors. 105, 261–269, 1998; Ondoh, T. Seismo ionospheric phenomena. Adv. Space Res. 26, 8, 1267–1272, 2000]. Characteristic phase changes at terminator times of Omega 10.2 kHz waves passing 70 km of the epicenter extended toward darker local times by 1 h for 3 days before this earthquake due to lowering of the wave reflection height or ion density increases in the D region [Hayakawa, M., Molchanov, O. A., Ondoh, T., Kawai, E. The precursory signature effect of the Kobe earthquake on VLF subionospheric signals. J. Commun. Res. La., 43, 00. 169–180, 1996]. The radon concentration in the atmosphere over Ashiya fault, Kobe [Yasuoka, Y., Shinogi, M. Anomaly in atmospheric radon concentration: a possible precursor of the 1995 Kobe, Japan, earthquake. Health Phys. 72(5), 759–761, 1997] and in the groundwater at 17 m well in Nishinomiya, Japan [Igarashi, G., Saeki, S., Takahata, N., Sumikawa, K., Tasaki, S., Sasaki, Y., Takahashi, M., Sano, Y. Ground-water radon anomaly before the Kobe earthquake in Japan. Science 269, 60–61, 1995] had gradually increased since 2 months before the M7.2 earthquake, increased suddenly in December 1994, and rapidly returned to the normal low level of October, 1994 [Yasuoka, Y., Shinogi, M. 1997. Anomaly in atmospheric radon concentration: a possible precursor of the 1995 Kobe. Japan, earthquake. Health Phys. 72(5), 759–761.]. Radon concentration changes in the groundwater before the M 7.0 Izu-Oshima-kinkai earthquake, Japan on January 14, 1978 [Wakita, H., Nakamura, Y., Notsu, K., Noguchi, M., Asada, T. 1980. Radon anomaly: a possible precursor of the 1978 Izu-Oshima-kinkai earthquake. Science 207, 882–883] and the M6.8 Chengkung earthquake, Taiwan on December 10, 2003 [Kuo, T., Fan, K., Chen, W., Kuochen, H., Han, Y., Wang, C., Chang, T., Lee, Y. Radon anomaly at the Antung Hot Spring before the Taiwan M6.8 Chengkung earthquake. Proceedings, Thirty-First Workshop on Geothermal Reservoir Engineering, Stanford University, Stanford, California, January 30–February 1, 2006, SGP-TR-179, 2006] are also investigated to find common features of the groundwater radon concentration changes before large earthquakes (M > 6.5) in comparison with those before the M7.2 Hyogoken–Nanbu earthquake. Groundwater radon concentrations before the 3 large earthquakes had shown common characteristic changes of gradually initial ones from the normal level since about 2 months before the earthquake onsets, rapid decreases down to the minimum, and quick increases up to the maximum at 7–20 days before the earthquake onsets, respectively. These are very useful characteristics of pre-seismic radon anomaly for the earthquake prediction or warning. Promising observations toward the earthquake prediction are also discussed.