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.     

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.     

Variations of ionospheric total electron content before three strong earthquakes in the Qinghai-Tibet region

We investigate the ionospheric total electron content (TEC) anomalies occurred in the Qinghai-Tibet region before three large earthquakes (M > 7.0). The temporal and spatial TEC variations were used to detect the ionospheric possible precursors of these earthquakes. We identified two TEC enhancements in the afternoon local time 9 days and 2–3 days before each earthquake, between which a TEC decrement occurred 3–6 days before earthquakes. These anomalies happened in the area of about 30° in latitude and the maximum is localized equatorward from the epicenters. These TEC anomalies can be found in all three earthquakes regardless the geomagnetic conditions. The features of these anomalies have something in common and may have differences from those caused by geomagnetic storms. Our results suggest that these ionospheric TEC perturbations may be precursors of the large earthquakes.     

A nonlinear background removal method for seismo-ionospheric anomaly analysis under a complex solar activity scenario: A case study of the M9.0 Tohoku earthquake

A precise determination of ionospheric total electron content (TEC) anomaly variations that are likely associated with large earthquakes as observed by global positioning system (GPS) requires the elimination of the ionospheric effect from irregular solar electromagnetic radiation. In particular, revealing the seismo-ionospheric anomalies when earthquakes occurred during periods of high solar activity is of utmost importance. To overcome this constraint, a multiresolution time series processing technique based on wavelet transform applicable to global ionosphere map (GIM) TEC data was used to remove the nonlinear effect from solar radiation for the earthquake that struck Tohoku, Japan, on 11 March, 2011. As a result, it was found that the extracted TEC have a good correlation with the measured solar extreme ultraviolet flux in 26–34 nm (EUV_26–34) and the 10.7 cm solar radio flux (F_10.7). After removing the influence of solar radiation origin in GIM TEC, the analysis results show that the TEC around the forthcoming epicenter and its conjugate were significantly enhanced in the afternoon period of 8 March 2011, 3 days before the earthquake. The spatial distributions of the TEC anomalous and extreme enhancements indicate that the earthquake preparation process had brought with a TEC anomaly area of size approximately 1650 and 5700 km in the latitudinal and longitudinal directions, respectively.     

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.     

Ionospheric anomalies of local earthquakes detected by GPS TEC measurements using data from Tashkent and Kitab stations

This paper reports the ionospheric anomalies observed during strong local earthquakes (M?5.0) which occurred mostly in and around Uzbekistan in seismically active zones, during years 2006 to 2009 within approximately 1000 km distance from the observing GPS stations located in Tashkent and Kitab, Uzbekistan. The solar and geomagnetic conditions were quiet during occurrence of the selected strong earthquakes. We produce Total Electron Content (TEC) time series over both sites and apply them to detect anomalous TEC signals preceding or accompanying the local earthquakes. The results show anomalous increase or decrease of TEC before or during the earthquakes. In general the anomalies occurred 1–7 days before the earthquakes as ionospheric electromagnetic precursors. To identify the anomalous values of TEC we calculated differential TEC (dTEC). dTEC is obtained by subtracting monthly averaged diurnal vTEC from the values of observed vTEC at each epoch. This procedure removes normal diurnal variations of vTEC. The present results are in good agreement with the previous observations on ionospheric earthquake precursors reported by various researchers.     

Use of principal component analysis in the identification of the spatial pattern of an ionospheric total electron content anomalies after China’s May 12, 2008, M = 7.9 Wenchuan earthquake

This paper uses principal component analysis (PCA) to determine the spatial pattern of total electron content (TEC) anomalies in the ionosphere post the China’s Wenchuan Earthquake of 12 May, 2008 (UTC) (M_w = 7.9). PCA and image processing are applied to the global ionospheric map (GIM) with transforms conducted for the time period 08:00–10:00 UT on 12 May 2008. Results show that at a height of approximately 200 km the anomaly is widespread and less intense; however, it becomes more localized with height reaching maximum intensity and localization at an altitude of 300 km. The spatial distribution is remarkably similar to that reported for a TEC anomaly previously identified as a precursor anomaly on May 9, 2008 for the same time period. Potential causes of the results are discussed with emphasis given to vertical acoustic gravity waves based on the spatial pattern identified.     

Ionospheric disturbances on 8 September, 2010: Was it connected with the incoming moderate Chongqing earthquake?

Earthquake prediction stimulates the searches for a correlation between seismic activity and ionospheric anomalies. Contrary to common focuses on strong earthquakes, we report the ionospheric disturbances, 2 days before a moderate Ms = 4.7 Chongqing earthquake (29.4°N, 105.5°E, depth = 7.0 km, occurred at 21:21 LT, 10 September, 2010) with the data of ground-based ionosondes and IGS receivers. The data covering the period under the quiet geomagnetic conditions and a geomagnetic storm was analyzed with upper and lower bounds. It is found that there were significant enhancements of foF2 and total electron content (TEC) on the afternoon of 8 September, 2010, with a limited area close to the epicentre, which was different from the feature of ionospheric perturbations triggered by the geomagnetic storm on 15 September. Taking into account the heliogeomagnetical condition, we conclude that the observed ionospheric enhancements were very likely associated with the forthcoming moderate Chongqing earthquake, which implies that the relationship between the amplitudes of ionospheric disturbances and earthquakes is very complicated.     

A morphological study of GPS-TEC data at Agra and their comparison with the IRI model

The temporal and seasonal variations of Total Electron Content (TEC) are studied at Agra (Geographic Lat. 27.17°N, Long. 78.89°E, Dip: 41.4°), India, which is in the equatorial anomaly region, for a period of 12 months from 01 January to 31 December, 2007 using a Global Positioning System (GPS) receiver. The mean TEC values show a minimum at 0500 h LT (LT = UT + 5.5 h) and a peak value at about 1400 h LT. The lowest TEC values are observed in winter whereas largest values are observed in equinox and summer. Anomalous variations are found during the period of magnetic disturbances. These results are compared with the TEC derived from IRI-2007 using three different options of topside electron density, NeQuick, IRI01-corr, and IRI-2001. A good agreement is found between the TEC obtained at Agra and those derived from IRI models.     

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.     

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.     

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.     

Ionospheric anomalies observed over South Korea preceding the Great Tohoku earthquake of 2011

We investigated the ionospheric anomalies observed before the Tohoku earthquake, which occurred near the northeast coast of Honshu, Japan on 11 March, 2011. Based on data from a ground-based Global Positioning System (GPS) network on the Korean Peninsula, ionospheric anomalies were detected in the total electron content (TEC) during the daytime a few days before earthquake. Ionospheric TEC anomalies appeared on 5, 8 and 11 March. In particular, the ionospheric disturbances on 8 March evidenced a remarkable increase in TEC. The GPS TEC variation associated with the Tohoku earthquake was an increase of approximately 20 total electron content units (TECU), observed simultaneously in local and global TEC measurements. To investigate these pre-earthquake ionospheric anomalies, space weather conditions such as the solar activity index (F10.7) and geomagnetic activity indices (the Kp and Dst indices) were examined. We also created two-dimensional TEC maps to visual the spatial variations in the ionospheric anomalies preceding the earthquake.     

Anomalous TEC variations associated with the strong Pakistan-Iran border region earthquake of 16 April 2013 at a low latitude station Agra,India

Employing a dual frequency GPS-receiver, ionospheric total electron (TEC) measurements have been in progress at Agra station in India (Geograph. lat. 27.2°N, long. 78°E) since 1 April 2006. In this paper, the TEC data have been analyzed for a period of one month from 1 April-1 May 2013 to examine the effect of multiple earthquakes, some of which occurred on the same day of 16 April 2013, and others occurred in the same month of April, 2013 in India and neighboring countries. We process the data using quartile and epoch analysis based statistical techniques and show that out of all the earthquakes, the one of the largest magnitude (M = 7.8) that occurred on Pakistan-Iran border caused anomalous enhancements and depletions in TEC 1–9 days before the occurrence of main sock. The E × B drift mechanism is suggested for the anomalies to occur in which the seismogenic electric field E is generated in a process suggested by Pulinets (2004).     

Surface temperature estimation in Singhbhum Shear Zone of India using Landsat-7 ETM+ thermal infrared data

Land surface temperature (LST) is an important factor in global change studies, heat balance and as control for climate change. A comparative study of LST over parts of the Singhbhum Shear Zone in India was undertaken using various emissivity and temperature retrieval algorithms applied on visible and near infrared (VNIR), and thermal infrared (TIR) bands of high resolution Landsat-7 ETM+ imagery. LST results obtained from satellite data of October 26, 2001 and November 2, 2001 through various algorithms were validated with ground measurements collected during satellite overpass. In addition, LST products of MODIS and ASTER were compared with Landsat-7 ETM+ and ground truth data to explore the possibility of using multi-sensor approach in LST monitoring. An image-based dark object subtraction (DOS3) algorithm, which is yet to be tested for LST retrieval, was applied on VNIR bands to obtain atmospheric corrected surface reflectance images. Normalized difference vegetation index (NDVI) was estimated from VNIR reflectance image. Various surface emissivity retrieval algorithms based on NDVI and vegetation proportion were applied to ascertain emissivities of the various land cover categories in the study area in the spectral range of 10.4–12.5 μm. A minimum emissivity value of about 0.95 was observed over the reflective rock body with a maximum of about 0.99 over dense forest. A strong correlation was established between Landsat ETM+ reflectance band 3 and emissivity. Single channel based algorithms were adopted for surface radiance and brightness temperature. Finally, emissivity correction was applied on ‘brightness temperature’ to obtain LST. Estimated LST values obtained from various algorithms were compared with field ground measurements for different land cover categories. LST values obtained after using Valor’s emissivity and single channel equations were best correlated with ground truth temperature. Minimum LST is observed over dense forest as about 26 °C and maximum LST is observed over rock body of about 38 °C. The estimated LST showed that rock bodies, bare soils and built-up areas exhibit higher surface temperatures, while water bodies, agricultural croplands and dense vegetations have lower surface temperatures during the daytime. The accuracy of the estimated LST was within ±2 °C. LST comparison of ASTER and MODIS with Landsat has a maximum difference of 2 °C. Strong correlation was found between LST and spectral radiance of band 6 of Landsat-7 ETM+. Result corroborates the fact that surface temperatures over land use/land cover types are greatly influenced by the amount of vegetation present.     

Regional modeling of the ionosphere using adaptive neuro-fuzzy inference system in Iran

In recent years, new techniques and algorithms such as Artificial Neural Networks (ANNs), Fuzzy Inference Systems (FIS) and Genetic Algorithm (GA) have been used as alternative statistical tools in modeling and forecasting issues. These methods have been extensively used in the field of geosciences and atmospheric physics. The main purpose of this paper is to combine FIS and ANNs for local modeling of the ionosphere Total Electron Content (TEC) in Iran. An Adaptive Neuro-Fuzzy Inference System (ANFIS) is developed for TEC modeling. Also, Multi-Layer Perceptron ANN (MLP-ANN) and ANN based on Radial Base Functions (RBF) have been designed for analyzing ANFIS results. Observations of 29 Global Positioning System (GPS) stations from the Iranian Permanent GPS Network (IPGN) have been used in 3 different seasons in 2015 and 2016. These stations are located at geomagnetic low latitudes region. Out of these 29 stations, 24 stations for training and 5 stations for testing and validating were selected. The relative and absolute errors have been used to evaluate the accuracy of the proposed model. Also, the results of this paper are compared with the International Reference Ionosphere model (IRI2016). The maximum values of the average relative error for RBF, MLP-ANN, ANFIS and IRI2016 methods are 13.88%, 11.79%, 10.06%, and 18.34%, respectively. Also, the maximum values of the average absolute error for these methods are 2.38, 2.21, 1.5 and 3.36 TECU, respectively. Comparison of diurnal predicted TEC from the ANFIS, RBF, MLP-ANN and IRI2016 models with GPS-TEC revealed that the ANFIS provides more accurate predictions than the other methods in the test area.     

The effect of total solar eclipse of October 3, 2005, on the total electron content over Europe

GPS observations from EUREF permanent GPS network were used to observe the response of TEC (Total Electron Content) to the total solar eclipse on October 3, 2005, under quiet geomagnetic conditions of the daytime ionosphere. The effect of the eclipse was detected in diurnal variations and more distinctly in the variations of TEC along individual satellite passes. The trough-like variations with a gradual decrease and followed by an increase of TEC at the time of the eclipse were observed over a large region. The depression of TEC amounted to 3–4 TECU. The maximum depression was observed over all stations located at the maximum path of the solar eclipse. The delay of a minimum level of TEC with respect to the maximum phase of the eclipse was about 20–30 min.     

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.     

Analysis of electron content variations over Japan during solar minimum: Observations and modeling

Comparative analysis of GPS TEC data and FORMOSAT-3/COSMIC radio occultation measurements was carried out for Japan region during period of the extremely prolonged solar minimum of cycle 23/24. COSMIC data for different seasons corresponded to equinox and solstices of the years 2007–2009 were analyzed. All selected electron density profiles were integrated up to the height of 700 km (altitude of COSMIC satellites), the monthly median estimates of Ionospheric Electron Content (IEC) were retrieved with use of spherical harmonics expansion. Monthly medians of TEC values were calculated from diurnal variations of GPS TEC estimates during considered month. Joint analysis of GPS TEC and COSMIC data allows us to extract and estimate electron content corresponded to the ionosphere (its bottom and topside parts) and the plasmasphere (h > 700 km) for different seasons of 2007–2009. Percentage contribution of ECpl to GPS TEC indicates the clear dependence from the time and varies from a minimum of about 25–50% during day-time to the value of 50–75% at night-time. Contribution of both bottom-side and topside IEC has minimal values during winter season in compare with summer season (for both day- and night-time). On average bottom-side IEC contributes about 5–10% of GPS TEC during night and about 20–27% during day-time. Topside IEC contributes about 15–20% of GPS TEC during night and about 35–40% during day-time. The obtained results were compared with TEC, IEC and ECpl estimates retrieved by Standard Plasmasphere–Ionosphere Model that has the plasmasphere extension up to 20,000 km (GPS orbit).