A Spatially Self-Adaptive Multiparametric Anomaly Identification Scheme Based on Global Strong Earthquakes

Earthquake forecasting aims to determine the likelihood of a damaging earthquake occurring in a particular area within a period of days to months. This provides ample preparation time for potential seismic hazards, resulting in significant socioeconomic benefits. Surface and atmospheric parameters derived from satellite thermal infrared observations have been utilized to identify pre-earthquake anomalies that may serve as potential precursors for earthquake forecasting. However, the correlation between these anomalies and impending earthquakes remains a significant challenge due to high false alarm and missed detection rates. To address this issue, we propose a spatially self-adaptive multiparametric anomaly identification scheme based on global strong earthquakes to establish the optimal recognition criteria. Each optimal parameter exhibits significant spatial variability within the seismically active region and indicates transient and subtle anomaly signals with a limited frequency of occurrences (<10 for most regions). In comparison to the fixed criterion for identifying anomalies, this new scheme significantly improves the positive Matthew’s correlation coefficient (MCC) values from ~0.03 to 0.122–0.152. Additionally, we have developed a multi-parameter anomaly synthesis method based on the best MCC value of each parameter anomaly. On average, the MCC increased from 0.143 to 0.186, and there are now more earthquake-prone regions with MCC values > 0.5. Our research emphasizes the critical importance of a multiparametric system in earthquake forecasting, where each geophysical parameter can be assigned a distinct weight, and the findings specifically identify OLR, including all-sky and clear-sky ones, as the most influential parameter on a global scale, highlighting the potential significance of OLR anomalies for seismic forecasting. Encouraging results imply the effectiveness of utilizing multiparametric anomalies and provide some confidence in advancing our knowledge of operational earthquake forecasting with a more quantitative approach.