We present a comprehensive analysis of pre-seismic anomalies as computed from the ground and space-based techniques during the recent Samos earthquake in Greece on 30 October 2020, with a magnitude M = 6.9. We proceed with a multi-parametric approach where pre-seismic irregularities are investigated in the stratosphere, ionosphere, and magnetosphere. We use the convenient methods of acoustics and electromagnetic channels of the Lithosphere–Atmosphere–Ionosphere-Coupling (LAIC) mechanism by investigating the Atmospheric Gravity Wave (AGW), magnetic field, electron density, Total Electron Content (TEC), and the energetic particle precipitation in the inner radiation belt. We incorporate two ground-based IGS GPS stations DYNG (Greece) and IZMI (Turkey) for computing the TEC and observed a significant enhancement in daily TEC variation around one week before the earthquake. For the space-based observation, we use multiple parameters as recorded from Low Earth Orbit (LEO) satellites. For the AGW, we use the SABER/TIMED satellite data and compute the potential energy of stratospheric AGW by using the atmospheric temperature profile. It is found that the maximum potential energy of such AGW is observed around six days before the earthquake. Similar AGW is also observed by the method of wavelet analysis in the fluctuation in TEC values. We observe significant energetic particle precipitation in the inner radiation belt over the earthquake epicenter due to the conventional concept of an ionospheric-magnetospheric coupling mechanism by using an NOAA satellite. We first eliminate the particle count rate (CR) due to possible geomagnetic storms and South Atlantic Anomaly (SAA) by the proper choice of magnetic field B values. After the removal of the statistical background CRs, we observe a significant enhancement of CR four and ten days before the mainshock. We use Swarm satellite outcomes to check the magnetic field and electron density profile over a region of earthquake preparation. We observe a significant enhancement in electron density one day before the earthquake. The parameters studied here show an overall pre-seismic anomaly from a duration of ten days to one day before the earthquake.

We focus on the possible thermal channel of the well-known Lithosphere–Atmosphere–Ionosphere Coupling (LAIC) mechanism to identify the behavior of thermal anomalies during and prior to strong seismic events. For this, we investigate the variation of Surface Latent Heat Flux (SLHF) as resulting from satellite observables. We demonstrate a spatio-temporal variation in the SLHF before and after a set of strong seismic events occurred in Kathmandu, Nepal, and Kumamoto, Japan, having magnitudes of 7.8, 7.3, and 7.0, respectively. Before the studied earthquake cases, significant enhancements in the SLHF were identified near the epicenters. Additionally, in order to check whether critical dynamics, as the signature of a complex phenomenon such as earthquake preparation, are reflected in the SLHF data, we performed a criticality analysis using the natural time analysis method. The approach to criticality was detected within one week before each mainshock.

Atmospheric disturbances caused by seismic activity are a complex phenomenon. The Lithosphere–Atmosphere–Ionosphere Coupling (LAIC) mechanism gives a detailed idea to understand these processes to study the possible impacts of a forthcoming earthquake. The atmospheric gravity wave (AGW) is one of the most accurate parameters for explaining such LAIC process, where seismogenic disturbances can be explained in terms of atmospheric waves caused by temperature changes. The key goal of this work is to study the perturbation in the potential energy associated with stratospheric AGW prior to many large earthquakes. We select seven large earthquakes having Richter scale magnitudes greater than seven ( ) in Japan (Tohoku and Kumamoto), Mexico (Chiapas), Nepal, and the Indian Ocean region, to study the intensification of AGW using the atmospheric temperature profile as recorded from the Sounding of the Atmosphere using Broadband Emission Radiometry (SABER) satellite. We observe a significant enhancement in the potential energy of the AGW ranging from 2 to 22 days prior to different earthquakes. We examine the conditions of geomagnetic disturbances, typhoons, and thunderstorms during our study and eliminate the possible contamination due to these events.

In this work we present the statistical and criticality analysis of the very low frequency (VLF) sub-ionospheric propagation data recorded by a VLF/LF radio receiver which has recently been established at the University of West Attica in Athens (Greece). We investigate a very recent, strong (M6.9), and shallow earthquake (EQ) that occurred on 30 October 2020, very close to the northern coast of the island of Samos (Greece). We focus on the reception data from two VLF transmitters, located in Turkey and Israel, on the basis that the EQ’s epicenter was located within or very close to the 5th Fresnel zone, respectively, of the corresponding sub-ionospheric propagation path. Firstly, we employed in our study the conventional analyses known as the nighttime fluctuation method (NFM) and the terminator time method (TTM), aiming to reveal any statistical anomalies prior to the EQ’s occurrence. These analyses revealed statistical anomalies in the studied sub-ionospheric propagation paths within ~2 weeks and a few days before the EQ’s occurrence. Secondly, we performed criticality analysis using two well-established complex systems’ time series analysis methods—the natural time (NT) analysis method, and the method of critical fluctuations (MCF). The NT analysis method was applied to the VLF propagation quantities of the NFM, revealing criticality indications over a period of ~2 weeks prior to the Samos EQ, whereas MCF was applied to the raw receiver amplitude data, uncovering the time excerpts of the analyzed time series that present criticality which were closest before the Samos EQ. Interestingly, power-law indications were also found shortly after the EQ’s occurrence. However, it is shown that these do not correspond to criticality related to EQ preparation processes. Finally, it is noted that no other complex space-sourced or geophysical phenomenon that could disturb the lower ionosphere did occur during the studied time period or close after, corroborating the view that our results prior to the Samos EQ are likely related to this mainshock.

Lithosphere–atmosphere–ionosphere coupling (LAIC) is studied through various physical or chemical quantities, obtained from different sources, which are observables of the involved complex processes. LAIC has been proposed to be achieved through three major channels: the chemical, the acoustic, and the electromagnetic. Accumulated evidence supporting the acoustic channel hypothesis has been published, while atmospheric gravity waves (AGWs) play a key role in LAIC as the leading mechanism for the transmission of energy from the lower atmosphere to the stratosphere and mesosphere, associated with atmospheric disturbances observed prior to strong earthquakes (EQs). The seismogenic AGW is the result of temperature disturbances, usually studied through stratospheric potential energy (E_P). In this work, we examined 11 cases of significant EQs (M > 6.7) that occurred during the last 10 years at different geographic areas by analyzing the temperature profile at the wider location of each one of the examined EQs. The “Sounding of the Atmosphere using Broadband Emission Radiometry” (SABER) instrument, part of the “Thermosphere Ionosphere Mesosphere Energetics Dynamics” (TIMED) satellite, data were employed to compute the potential energy (E_P) of the AGW. Using the temperature profile, we first calculated E_P and determined the altitudes’ range for which prominent pre-seismic disturbances were observed. Subsequently, the E_P time series at specific altitudes, within the determined “disturbed” range, were for the first time analyzed using the criticality analysis method termed the “natural time” (NT) method in order to find any evidence of an approach to a critical state (during a phase transition from a symmetric phase to a low symmetry phase) prior to the EQ occurrence. Our results show criticality indications in the fluctuation of E_P a few days (1 to 15 days) prior to the examined EQs, except from one case. In our study, we also examined all of the temperature-related extreme phenomena that have occurred near the examined geographic areas, in order to take into account any possible non-seismic influence on the obtained results.

In this study, we present a six-year (2014–2020) statistical analysis of VLF subionospheric propagation data at 19 VLF receivers from the VLF transmitter with call name JJI and frequency 22.2 kHz, all located in Japan. Moderate and strong earthquakes (EQs) ($M_L\ge 4.5$ and depth ≤ 50 km) that occurred in the wider area around Japan during the same time period with the available VLF data are investigated. The terminator times’ (TT) shift in VLF amplitude data as a possible precursor of an EQ are statistically examined, focusing on the correlation with seismic activity. The concept of the effective EQ magnitude ($M_{eff})$ is used in order to define the total EQ energy possibly affecting the midpoint of each path for each day. It is important to note that dates when geomagnetic storms or solar flares occurred as well as dates corresponding to the already known winter effect on TT statistics in the north–south direction were excluded. The cross-correlation between TT statistical anomalies and seismic activity, represented by $M_{eff}$, was extracted. Maximum cross-correlation values were found for most of the cases prior to the subsequent seismic activity, indicating a link between the ionospheric anomalies and the subsequent seismicity. Finally, the wide temporal range of the cross-correlation maxima temporal locations is justified by the inhomogeneity of the lower ionosphere, coupled with the anisotropy of the preseismic effect of the impending seismicity, highlighting the complexity of the EQ preparation processes.