Thermal anomalies detection before strong earthquakes (<i>M</i> > 6.0) using interquartile, wavelet and Kalman filter methods

Thermal anomaly is known as a significant precursor of strong earthquakes, therefore Land Surface Temperature (LST) time series have been analyzed in this study to locate relevant anomalous variations prior to the Bam (26 December 2003), Zarand (22 February 2005) and Borujerd (31 March 2006) earthquakes. The duration of the three datasets which are comprised of MODIS LST images is 44, 28 and 46 days for the Bam, Zarand and Borujerd earthquakes, respectively. In order to exclude variations of LST from temperature seasonal effects, Air Temperature (AT) data derived from the meteorological stations close to the earthquakes epicenters have been taken into account. The detection of thermal anomalies has been assessed using interquartile, wavelet transform and Kalman filter methods, each presenting its own independent property in anomaly detection. The interquartile method has been used to construct the higher and lower bounds in LST data to detect disturbed states outside the bounds which might be associated with impending earthquakes. The wavelet transform method has been used to locate local maxima within each time series of LST data for identifying earthquake anomalies by a predefined threshold. Also, the prediction property of the Kalman filter has been used in the detection process of prominent LST anomalies. The results concerning the methodology indicate that the interquartile method is capable of detecting the highest intensity anomaly values, the wavelet transform is sensitive to sudden changes, and the Kalman filter method significantly detects the highest unpredictable variations of LST. The three methods detected anomalous occurrences during 1 to 20 days prior to the earthquakes showing close agreement in results found between the different applied methods on LST data in the detection of pre-seismic anomalies. The proposed method for anomaly detection was also applied on regions irrelevant to earthquakes for which no anomaly was detected, indicating that the anomalous behaviors can be related to impending earthquakes. The proposed method receives its credibility from the overall capabilities of the three integrated methods

First results on the systematic search of land surface temperature anomalies as earthquakes precursors

Every year, earthquakes cause thousands of casualties and high economic losses. For example, in the time frame from 1998 to 2018, the total number of casualties due to earthquakes was larger than 846 thousand people, and the recorded economic losses were about USD 661 billion. At present, there are no earthquake precursors that can be used to trigger a warning. However, some studies have analyzed land surface temperature (LST) anomalies as a potential earthquake precursor. In this study, a large database of global LST data from the Geostationary Operational Environmental Satellite (GOES) and AQUA satellites during the whole year 2020 has been used to study the LST anomalies in the areas affected by earthquakes. A total of 1350 earthquakes with a magnitude larger than M4 were analyzed. Two methods widely used in the literature have been used to detect LST anomalies in the detrended LST time series: the interquartile (IQT) method and the standard deviation (STD). To the authors’ knowledge, it is the first time that the confusion matrix (CM), the receiver operating characteristic curve (ROC), and some other figures of merit (FoM) are used to assess and optimize the performance of the methods, and to select the optimum combination that could be used as a proxy for their occurrence. A positive anomaly was found a few days before the studied earthquakes, followed by the LST decrease after the event. Further studies over larger regions and more extended periods will be needed to consolidate these encouraging results.This work was sponsored by project “GENESIS: GNSS Environmental and Societal Missions—Subproject UPC”, Grant PID2021-126436OB-C21, sponsored by MCIN/AEI/10.13039/501100011033/ and EU ERDF “A way to do Europe”. Badr-Eddine Boudriki Semlali received support in the form of an FI grant: 2021 FI_B 00471 from Generalitat de Catalunya—FI AGAUR 2021.Peer ReviewedPostprint (published version

From visual comparison to robust satellite techniques: 30 years of thermal infrared satellite data analyses for the study of earthquake preparation phases

This review paper reports the main contributions and results achieved after more than 30 years of studies on the possible relationships among space-time variation of Earth’s thermally emitted radiation, measured by satellite sensors operating in the Thermal InfraRed (TIR) spectral range (8-14 m), and earthquake occurrence. Focus will be given on the different existing methods/models to: 1) discriminate a possible pre-seismic TIR anomaly from all the other TIR signal fluctuations; 2) correlate such anomalies with space, time and magnitude of earthquakes; 3) physically justify such a correlation

Study of outgoing longwave radiation anomalies associated with Haiti earthquake

The paper presents an analysis by using the methods of Eddy field calculation mean and wavelet maxima to detect seismic anomalies within the outgoing longwave radiation (OLR) data based on time and space. The distinguishing feature of the method of Eddy field calculation mean is that we can calculate "the total sum of the difference value" of "the measured value" between adjacent points, which could highlight the singularity within data. The identified singularities are further validated by wavelet maxima, which using wavelet transformations as data mining tools by computing the maxima that can be used to identify obvious anomalies within OLR data. The two methods has been applied to carry out a comparative analysis of OLR data associated with the earthquake recently occurred in Haiti on 12 January 2010. Combining with the tectonic explanation of spatial and temporal continuity of the abnormal phenomena, the analyzed results have indicated a number of singularities associated with the possible seismic anomalies of the earthquake and from the comparative experiments and analyses by using the two methods, which follow the same time and space, we conclude that the singularities observed from 19 to 24 December 2009 could be the earthquake precursor of Haiti earthquake

Coarse-graining research of the thermal infrared anomalies before earthquakes in the Sichuan area on Google Earth engine

Seismo-induced Thermal infrared (TIR) anomalies has been proposed as a significant precursor of earthquakes. Several methods have been proposed to detect Thermal infrared anomalies that may be associated with earthquakes. However, there is no comparison of the influence for Thermal infrared extraction methods with a long time statistical analysis. To quantify the effects of various techniques used in Thermal infrared anomaly extraction, in this paper, we offer a complete workflow of their comparative impacts. This study was divided into three parts: anomaly detection, statistical analysis, and tectonic factor research. For anomaly detection, daily continuous nighttime surface temperature (ConLST) data was obtained from the Google Earth Engine (GEE) platform, and each different anomaly detection method was used to detect Thermal infrared outliers in the Sichuan region (27°-37°N, 97°-107°E). During statistical analysis, The heated core model was applied to explore Thermal infrared anomalies which is to filter anomalies unrelated to earthquakes by setting time-space-intensity conditions. The 3D error diagram offers scores to assume the best parameter set using training-test-validation steps. In the final part, we considered information on stresses, active faults, and seismic zones to determine the optimal parameters for extracting the Thermal infrared anomalies. The Kalman filter method detected the highest seismic anomaly frequency without considerating the heating core condition. The Autoencoder and Isolation Forest methods obtain the optimal alert type and parameter set to determine if the anomaly is likely earthquake-related. The RST method performs optimally in the final part of the workflow when it considers physical factors such as active faults, seismic zones, and stresses. However, The six methods we have chosen are not sufficient to contain the entire Thermal infrared anomaly extraction. The consideration of tectonic factors in the research remains poorly developed, as statistical methods were not employed to explore the role of constructive factors. Nevertheless, it is a significant factor in comparing anomaly extraction methods and precursor studies

Searching for possible seismogenic signatures in ionosphere by an entropy-based analysis of magnetic satellite data: A case study.

The importance of detecting possible electromagnetic signatures due to large earthquakes is self-evident, signatures which can be either anticipating, simultaneous or subsequent with respect to the main shock. Taking advantage of the present low Earth’s orbiting CHAMP satellite, we apply an “ad hoc” technique based on the Information Theory, to the satellite magnetic data with the aim at extracting eventual time anomalies. This technique has small time-space resolution using a preliminary wavelet analysis in order to detect shorter-wavelength anomalies. Some examples are given for magnetic satellite data taken over periods including the times of two large earthquakes, one being the Sumatra region event on 26 December 2004 (M=9.1)

Thermal Radiation Anomalies Associated with Major Earthquakes

Recent developments of remote sensing methods for Earth satellite data analysis contribute to our understanding of earthquake related thermal anomalies. It was realized that the thermal heat fluxes over areas of earthquake preparation is a result of air ionization by radon (and other gases) and consequent water vapor condensation on newly formed ions. Latent heat (LH) is released as a result of this process and leads to the formation of local thermal radiation anomalies (TRA) known as OLR (outgoing Longwave radiation, Ouzounov et al, 2007). We compare the LH energy, obtained by integrating surface latent heat flux (SLHF) over the area and time with released energies associated with these events. Extended studies of the TRA using the data from the most recent major earthquakes allowed establishing the main morphological features. It was also established that the TRA are the part of more complex chain of the short-term pre-earthquake generation, which is explained within the framework of a lithosphere-atmosphere coupling processes