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.

We present perturbations in very low frequency (VLF) signals received at Ionospheric & Earthquake Research Centre (IERC) (Lat. 22.50 ${}^{°}$N, Long. 87.48 ${}^{°}$E) during and prior to two earthquakes, one on 11 March 2011 at 11:16:24 LT (M = 9) in Honshu, Japan and another on 12 May 2015 at 12:35:19 LT (M = 7.3) in Kodari, Nepal. The VLF signal transmitted from JJI (22.2 kHz) in Japan (Lat. 32.05 ${}^{°}$N, Long. 131.51 ${}^{°}$E) showed strong shift in VLF-sunrise terminator times towards night-time few days prior to both the earthquakes. These two earthquakes took place near the transmitter JJI and receiver IERC respectively. We simulated the VLF sunrise terminator time shifts using the long wavelength propagation capability (LWPC) code. To effectively represent the D-region ionospheric variabilities, we assumed a mean dynamic perturbation over the path and presented them with a set of effective Wait’s parameters (β_eff, $h_{\mathit{eff}}^{\prime }$). We have successfully reproduced the temporal trend of the normalized VLF signal amplitude at VLF sunrise terminators for a few days around both the earthquakes. We used Wait’s exponential model for estimating the altitude profile of D-region electron density ( $N_e(h)$) at VLF sunrise terminator times around both the earthquakes and quantified the changes of those $N_e(h)$ profiles.

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.

This study examines the response of the thermal channel within the Lithosphere–Atmosphere–Ionosphere Coupling (LAIC) mechanism during the notable earthquake in Crete, Greece, on 27 September 2021. We analyze spatio-temporal profiles of Surface Latent Heat Flux (SLHF), Outgoing Longwave Radiation (OLR), and Atmospheric Chemical Potential (ACP) using reanalysis data from the National Oceanic and Atmospheric Administration (NOAA) satellite. Anomalies in these parameters are computed by removing the background profile for a non-seismic condition. Our findings reveal a substantial anomalous increase in these parameters near the earthquake’s epicenter 3 to 7 days before the main shock. The implications of these observations contribute to a deeper understanding of the LAIC mechanism’s thermal channel in seismic events.