Wave-like perturbations in the ionospheric F2-layer observed after the Ms8.1 Samoa earthquake of September 29, 2009

A study of the critical frequency foF2 variations after the large earthquake (M_s = 8.1) which occurred on 29 September, 2009 in the region of Samoa Islands in the Pacific Ocean is carried out using data of the ionospheric station of Kwajalein. The epicenter of the earthquake was located at about 184 km southwest from Apia (the capital of West Samoa). It was found that wave-like perturbations of foF2 were observed for ∼3 h above the station (located approximately 3560 km northwest from the epicenter). The amplitude of the disturbance was as large as ∼20% of the average magnetic quiet day foF2 values. A comparison of the observed perturbations of foF2 with the ones detected at Stanford ionospheric station after the Alaska earthquake of 28 March 1964 (M_s = 8.4) showed a close similarity of the wave-like perturbations of foF2 in both cases.     

Initial accuracy and reliability of current BDS-3 precise positioning,velocity estimation,and time transfer (PVT)

The primary system of Chinese global BeiDou satellite system (BDS-3) was completed to provide global services on December 27, 2018; this was a key milestone in the development process for BDS in terms of its provision of global services. Therefore, this study analyzed the current performance of BDS-3, including its precise positioning, velocity estimation, and time transfer (PVT). The datasets were derived from international GNSS monitoring and assessment system (iGMAS) tracking networks and the two international time laboratories in collaboration with the International Bureau of Weights and Measures (BIPM). With respect to the positioning, the focus is on the real-time kinematic (RTK) positioning and precise point positioning (PPP) modes with static and kinematic scenarios. The results show that the mean available satellite number is 4.8 for current BDS-3 system at short baseline XIA1–XIA3. The RTK accuracy for three components is generally within cm level; the 3D mean accuracy is 8.9 mm for BDS-3 solutions. For the PPP scenarios, the convergence time is about 4 h for TP01 and BRCH stations in two scenarios. After the convergence, the horizontal positioning accuracy is better than cm level and the vertical accuracy nearly reaches the 1 dm level. With respect to kinematic scenarios, the accuracy stays at the cm level for horizontal components and dm level for the vertical component at two stations. In terms of velocity estimation, the horizontal accuracy stays at a sub-mm level, and the vertical accuracy is better than 2 mm/s in the BDS-3 scenario, even in the Arctic. In terms of time and frequency transfer, the noise level of BDS-3 time links can reach 0.096 ns for long-distances link NT01–TP02 and 0.016 ns for short-distance links TP01–TP02. Frequency stability reaches 5E–14 accuracy when the averaging time is within 10,000 s for NT01–TP02 and 1E–15 for TP01–TP02.     

Retrieval of vegetation equivalent water thickness from reflectance using genetic algorithm (GA)-partial least squares (PLS) regression

This study aimed to investigate the performance of genetic algorithms coupled with partial least squares (GA-PLS) modeling of spectral reflectance in retrieving equivalent water thickness (EWT) at leaf and canopy level. A genetic algorithm was used to identify a subset of spectral bands sensitive to the variation in EWT, and PLS was then applied to relate the identified bands to EWT values. GA-PLS was applied to leaf level reflectance available from LOPEX dataset, and canopy data includes reflectance simulated by a leaf radiative transfer model PROSPECT and a canopy radiative transfer model SAILH and acquired by airborne visible/infrared imaging spectrometer (AVIRIS). The results indicate that GA-PLS has the capability of retrieving EWT from leaf and canopy reflectance, and achieved good estimation accuracy, i.e. low root mean square errors (RMSE) and high squared correlation coefficients (R²). For the retrieval at leaf level, the estimation accuracy can be as good as RMSE = 0.0019 g/cm² and R² = 0.939 or better. For the retrieval at canopy level, the model accuracy is RMSE = 0.0061 g/cm² and R² = 0.966 or better when PROSPECT-SAILH simulated reflectance was used; when AVIRIS image spectra were used, the model accuracy is RMSE = 0.0094 g/cm² and R² = 0.8734 for the calibration, and RMSE = 0.0132 g/cm² and R² = 0.7756 for the validation. These results from GA-PLS modeling support the conclusion that GA-PLS has the potential to be applied to AVIRIS, Hyperion and HyMap imagery for retrieving EWT. The selected bands for the AVIRIS datasets differ from those for the LOPEX and PROSPECT-SAILH simulated datasets, and this inconsistence of the selected bands for different datasets indicates that the GA-PLS method has the advantage of tuning the optimum bands for PLS regression and accommodating the effects of confounding factors.     

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.     

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.     

Study of aerosol anomaly associated with large earthquakes (M>6)

We present the variation of unusual atmospheric phenomena, aerosols, to understand the preseismic irregularities for two major earthquakes in Japan. We consider aerosol optical depth and Angstrom exponent data retrieved from the Moderate Resolution Imaging Spectroradiometer (MODIS) instrument onboard the Terra satellite to establish possible connections between earthquakes and the generation of aerosols. Variation of the aerosol parameters shows significant changes before the April 15, 2016, Kumamoto earthquake ($M=7.0,h=10$ km) and the November 21, 2016, Fukushima earthquake ($M=6.9$ and $h=9$ km), where M indicates the Richter magnitude and h indicates the focal depth. To identify the source of the aerosol particles, we use the Hybrid Single-Particle Lagrangian Integrated Trajectory model (HYSPLIT-4). This model uses both Lagrangian and Eulerian approaches to compute trajectories and establish a source-receptor relationship. We compute backward trajectories to check whether the aerosol generated near the epicenter is due to the preseismic processes or is transported from other areas. From our results, we conclude the fine-mode aerosols are generated in the vicinity of the epicenter, 3–7 days before the earthquakes.     

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.     

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.     

Ionospheric foF2 variations prior to strong earthquakes in Taiwan area

Many studies of the seismo-ionospheric coupling effects have been reported. On 17 July 1998(M=6.2), 20 September 1999 (M=7.3) and 22 October 1999 (M=6.4) three large earthquakes respectively struck Rei-Li, Chi-Chi and Chia-Yi in central Taiwan. The three earthquakes severely damaged structures, heavily changed landforms and disturbed geophysical environments. This paper examines variations of the ionospheric penetration frequency, foF2, observed by Chung-Li ionosonde station (25.0° N, 121.1° E) several days before the three earthquakes. The mean- and median-based statistical techniques are introduced to investigate the ionospheric electron density prior to the three earthquakes. Results show that the foF2 decrease significantly before the three earthquakes.     

强震前异常电场激发声重波对电离层的影响

震前地震孕育期地表异常增强的电场,通过大气电导率传输到电离层高度.该异常电场通过非稳态局部加热,可以在电离层高度激发声重波.基于该理论,利用一维时变中纬电离层物理模型,模拟了该扰动源对电离层电子密度的影响.结果表明,重力波引起的中性风速度扰动对电离层电子密度分布影响甚微,该机理无法解释震前电离层异常扰动现象.      

The manifestation of seismic activity in 557.7 nm emission variations of the earth's upper atmosphere

For the Baikal seismic zone (East Siberia region, 52°N, 103°E) the Earth's upper atmosphere radiation parameters were analysed for 557.7 nm emission line of atomic oxygen [OI] (de-excitation heights 85–115 km) during earthquakes. The analysis revealed an increase in mean night-time 557.7 nm emission intensities prior to earthquakes, with a decrease in intensities during post-earthquake days. Preliminary results of a spectral analysis suggest possible changes in spectra (from a few to tens of minutes) of 557.7 nm emission prior to and after earthquakes. The results are compared with those obtained for the earthquakes in the Central Asian region.     

Evaluation of SARAL/AltiKa performance using GNSS/IMU equipped buoy in Sajafi,Imam Hassan and Kangan Ports

Performance of SARAL/AltiKa mission has been evaluated within 2016 altimeter calibration/validation framework in Persian Gulf through three campaigns conducted in the offshore waters of Sajafi, Imam Hassan and Kangan Ports, while the altimeter overflew the passes 470, 111 and 25 on 13 Feb, 7 March and 17 June 2016, respectively. As the preparation, a lightweight buoy was equipped with a GNSS receiver/choke-ring antenna and a MEMS-based IMU to measure independent datasets in the field operations. To obtain accurate sea surface height (SSH) time series, the offset of the onboard antenna from the equilibrium sea level was predetermined through surveying operations as the buoy was deploying in the onshore waters of Kangan Port. Accordingly, the double-difference carrier phase observations have been processed via the Bernese GPS Software v. 5.0 so as to provide the GNSS-derived time series at the comparison points of the calibration campaigns, once the disturbing effects due to the platform tilt and heave have been eliminated. Owing to comparing of the SSH time series and the associating altimetry 1?Hz GDR-T datasets, the calibration/validation of the SARAL/AltiKa has been performed in the both cases of radiometer and ECMWF wet troposphere corrections so as to identify potential land contamination. An agreement of the present findings in comparison with those attained in other international calibrations sites confirms the promising feasibility of Persian Gulf as a new dedicated site for calibration/validation of ongoing and future altimetry missions.     

Investigation of the relationship between ionospheric foF2 and earthquakes

Variations of the ionospheric F2 region critical frequency (foF2) have been investigated statistically before earthquakes during 1980–2008 periods in Japan area. Ionosonde data was taken from Kokubunji station which is in the earthquake preparation zone for all earthquakes. Standard Deviations and Inter-Quartile Range methods are applied to the foF2 data. It is observed that there are anomalous variations in foF2 before earthquakes. These variations can be regarded as ionospheric precursors and may be used for earthquake prediction.     

Precursory signature of several major earthquakes studied using 40 kHz low frequency signal

To investigate the precursory signature of earthquakes on low frequency (LF) signal propagation, six earthquakes, having magnitude greater than equal to 6.5 and depth less than equal to 30 km, are being studied. The base line level of 40 kHz signal, transmitted from JJY station, Japan, is analysed with respect to V_d statistical parameter. Results show that the V_d parameter values starts fluctuating from its ambient levels before and during the days of the earthquakes, with significant variation starting 1–3 days prior to the earthquake concerned. This present study is an approach for identifying the precursory signatures of earthquakes on LF signal propagation using a new methodology with V_d parameter.     

Identification of seismic activity sources on the subsatellite track by ionospheric plasma disturbances detected with the Sich-2 onboard probes

Using a cylindrical Langmuir probe and the authors’ proprietary two-channel pressure transducer, ionospheric plasma parameter distributions along the orbit of the Sich-2 satellite (Ukraine, 2011–2012) were measured. This paper is concerned with identifying the space–time location of ionospheric plasma disturbance sources, including the epicenters of actual earthquakes (before or during the satellite flyover) and incipient earthquakes on the subsatellite track, from the measured distributions of the electron density and temperature and the neutral particle temperature along the satellite orbit. To do this, the measured ionospheric plasma parameter distributions are connected to the coordinates on the subsatellite track.It is shown that local disturbances in the electron density and temperature and neutral particle temperature distributions in the satellite orbit in the ionosphere may serve as indicators of seismic activity on the subsatellite track. The epicenters of incipient earthquakes may be set off from other plasma parameter disturbance sources associated with seismic activity using information provided by special monitoring and survey centers that monitor the current seismic situation.     

Variations of equatorial electrojet as possible seismo-ionospheric precursor at the occurrence of TEC anomalies before strong earthquake

The measurements of GPS signal delays show that the local areas of increased/decreased Total Electron Content (TEC) of the ionosphere can be observed before strong earthquakes. The main possible cause of these TEC disturbances is the vertical plasma drift under the action of zonal electric field. The spatial pattern of electric potentials for such electric field was proposed. The model calculations were done to investigate the efficiency of the proposed mechanism. The calculation results revealed the agreement with TEC variations observed before strong earthquakes and showed that the equatorial electrojet variations can be considered as precursors of earthquakes.     

Sunspot numbers,interplanetary magnetic field,and cosmic ray intensity at earth: Nexus for the twentieth century

The pivotal role played by the interplanetary magnetic field (B) in modulating galactic cosmic ray (GCR) intensity in the heliosphere is described. We show that the inverse correlation observed by Forbush (1958) between GCRs and sunspot numbers (SSNs) is reflected in high correlation between SSNs and B (cc = 0.94). The SSN data are available since 1700 and the derived B data since 1835. The paleo-cosmic ray data are available for several millennia in the form of ¹⁰Be radionuclide sequestered in polar ice. The data of the ion chambers (ICs) at the Cheltenham–Fredericksburg–Yakutsk (CFY) sites are combined to create a data string for 1937–1988. In turn, these data are used to extend the measurements of the low energy GCR ions (>0.1 GeV) at balloon altitudes at high latitudes in Russia to 1937. These data are then correlated to B and the fit parameters are used to extend the low energy ion data to 1900, creating the instrumental era GCR time series for the twentieth century. The derived GCR time series is compared to ¹⁰Be measured at two sites in Greenland, namely Dye 3 and NGRIP for 1900–2000 to check the internal consistency of datasets for the long-term trend. We find that the annual mean rate (%) for 1965 at NGRIP is an outlier. We replace it with the mean of 1964 and 1965 rates and construct a new re-normalized time series at NGIP, improving the agreement with the derived instrumental era GCR time series for the twentieth century as well. This should encourage its use by heliophysics community for varied applications.     

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).     

Ionospheric perturbations due to earthquakes as determined from VLF and GPS-TEC data analysis at Agra,India

Employing SoftPAL receiver, amplitude variations of VLF transmitter signals NWC (19.8?kHz) and NPM (21.4?kHz) are analyzed at Agra station in India (Geograph. lat. 27.2°N, long. 78°E) ±15?days from five major earthquakes of magnitude M?=?6.9–8.5 occurred in Indian subcontinent during the years 2011–2013. We apply nighttime fluctuation (NF) method and show that in almost all cases the trend decreases and dispersion and NF increase on the same days corresponding to each earthquake about 11–15?days prior to the main shock. Assuming that the ionospheric perturbations are caused by atmospheric gravity waves (AGW), we also calculate AGW modulation index for each case and find its values increased on the days amplitude fluctuations take place. Its value is decreased in one case only where the perturbations may be attributed to penetration of seismogenic electric field.In order to support the above results we also present GPS-TEC data analyzed by us corresponding to three of the above earthquakes. We study the TEC anomalies (unusual enhancements) and find that in one case the precursory period is almost the same as that found in NF method.     

An early warning system for coastal earthquakes

Earthquakes are very common natural hazards, and every year a few damaging events occur throughout the globe. Satellite remote sensing data are found to be useful in providing information about changes in land, ocean, atmosphere and ionosphere. Recent studies have shown that some of the parameters observed from the remote sensing data are associated with impending coastal earthquakes. An automatic system “CQuake” has been developed to carry out spatial and temporal data mining analysis in real time. CQuake performs retrospective analysis of earthquakes, and performs forecast for predefined regions of the world based on the analysis of Surface latent heat flux (SLHF) or any other geophysical parameters. Details of CQuake, and an example of the Colima earthquake of January 22, 2003, are discussed.