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.

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.