Taking the Luding Ms6.8 earthquake (EQ) on 5 September 2022 as a case study, we investigated the potential seismic anomalies of the ionosphere, infrared radiation, atmospheric electrostatic field (AEF), and hot spring ions in the seismogenic region. Firstly, we analyzed the multi-parameter anomalies in the ionosphere around the epicenter and found synchronous anomalous disturbances in the ground parameters, namely the global ionospheric map (GIM), GPS, TEC, and satellite parameters, such as the He⁺ and O⁺ densities on 26 August under relatively quiet solar–geomagnetic conditions (F10.7 < 120 SFU; Kp < 3; Dst > −30 nT; |AE| < 500 nT). Next, both the anomaly analysis of the infrared radiation and AEF, and the survey results of the Luding EQ scientific expedition on the hot spring ions showed pre-seismic anomalous variations at different time periods in the seismogenic region. The characteristics of Earth’s multi-sphere coupling anomalies in temporal evolution and spatial distribution were obvious, which validated the Lithosphere–Atmosphere–Ionosphere Coupling (LAIC) mechanism. Finally, combining the analysis results and the LAIC mechanism, we suggested that the multi-sphere coupling anomalies were more likely associated with the Luding Ms6.8 EQ, and that the differential motion and the regional crustal stress accumulation between the Chuandian block and the Bayan Har block might have led to this EQ. Furthermore, remote sensing and ground-based monitoring technologies can play an important role in corroborating and compensating each other, while further study of the multi-sphere coupling mechanism will provide a clearer understanding of the seismogenic process for major EQs.

In this paper, we focused on the characteristics of the seismo-ionospheric effects related to two successive earthquakes, namely, the earthquakes in 2022 in Taitung Sea, Taiwan, China, with magnitudes (M) of 6.7 and 6.3, at 23.45° N, 121.55° E and 23.39° N, 121.52° E and with the same focal depth of 20 km, which were detected by the China Seismo-Electromagnetic Satellite (CSES). By applying the sliding interquartile range method to electron density (Ne) data acquired by the Langmuir probe (LAP) onboard the CSES and the grid total electron content (TEC) data obtained from the Center for Orbit Determination in Europe (CODE), positive anomalies were found under quiet geomagnetic conditions on 2–3 March and 8–9 March 2022—that is, 19–20 and 13–14 d before the earthquakes, respectively, and the global ionospheric mapping (GIM) TEC data suggested that anomalies may also have been triggered in the magnetic conjugate area 13–14 d prior to the earthquakes occurrences. In addition, the CSES Ne data showed enhancements 3 and 5 d before the earthquakes occurred. Furthermore, 138 earthquakes with M ≥ 5.0 that occurred in Taiwan and the surrounding region during the period February 2019 to March 2022 were statistically analyzed using the CSES Ne data. The results show that most of the Ne anomalies were positive. Moreover, the greater the earthquake magnitude, the greater the frequency of the anomalies; however, the amplitude of the anomalies did not increase with the earthquake magnitude. The anomalies were concentrated during the period of 10 d before to 5 d after the earthquakes. No increase in the amplitude of anomalies was observed as the time of the earthquakes approached. Finally, based on evidence relating to earthquake precursor anomalies, we conclude that it is possible that earthquakes in Taiwan and the surrounding region affect the ionosphere through the geochemical, acoustic, and electromagnetic channels, as described by the lithosphere–atmosphere–ionosphere coupling (LAIC) model, and that the two studied earthquakes in Taiwan may have induced ionospheric effects through the geochemical channel.

Based on the multi-data of the global ionospheric map (GIM), ionospheric total electron content (TEC) inversed from GPS observations, the critical frequency of the F2 layer (f_OF₂) from the ionosonde, electron density (Ne), electron temperature (Te), and He⁺ and O⁺ densities detected by the China Seismo-Electromagnetic Satellite (CSES), the temporal and spatial characteristics of ionospheric multi-parameter perturbations were analyzed around the Maerkang Ms6.0 earthquake swarm on 9 June 2022. The results showed that the seismo-ionospheric disturbances were observed during 2–4 June around the epicenter under quiet solar-geomagnetic conditions. All parameters we studied were characterized by synchronous changes and negative anomalies, with a better consistency between ionospheric ground-based and satellite observations. The negative ionospheric anomalies for all parameters appeared 5–7 days before the Maerkang Ms6.0 earthquake swarm can be considered as significant signals of upcoming main shock. The seismo-ionospheric coupling mechanism may be a combination of two coupling channels: an overlapped DC electric field and an acoustic gravity wave, as described by the lithosphere–atmosphere–ionosphere coupling (LAIC). In addition, in order to make the investigations still more convincing, we completed a statistical analysis for the ionospheric anomalies of earthquakes over Ms6.0 in the study area (20°~40° N, 92°~112° E) from 1 January 2019 to 1 July 2022. The nine seismic events reveal that most strong earthquakes are preceded by obvious synchronous anomalies from ground-based and satellite ionospheric observations. The anomalous disturbances generally appear 1–15 days before the earthquakes, and the continuity and reliability of ground-based ionospheric anomaly detection are relatively high. Based on the integrated ionospheric satellite–ground observations, a cross-validation analysis can effectively improve the confidence level of anomaly identification and reduce the frequency of false anomalies.

This study presents the spatio–temporal evolution of the electric and magnetic fields recorded by the China Seismo–Electromagnetic Satellite (CSES) and the thermal infrared remote sensing data observed by the Chinese stationary meteorological satellites Feng Yun–2G (FY–2G) associated with the 2021 Mw7.3 Maduo earthquake. Specifically, we analyzed the power spectrum density (PSD) data of the electric field in the extremely low frequency (ELF) band, the geomagnetic east–west vector data, and the temperature of brightness blackbody (TBB) data to investigate the spatio–temporal evolution characteristics under quiet space weather conditions (Dst > −30 nT and Kp < 3). Results showed that (1) the TBB radiation began to increase notably along the northern fault of the epicenter ~1.5 months prior to the occurrence of the earthquake. It achieved its maximum intensity on 17 May, and the earthquake occurred as the anomalies decreased. (2) The PSD in the 371 Hz–500 Hz and 700 Hz–871 Hz bands exhibited anomaly perturbations near the epicenter and its magnetic conjugate area on May 17, with particularly notable perturbations observed in the latter. The anomaly perturbations began to occur ~1 month before the earthquake, and the earthquake occurred as the anomalies decreased. (3) Both the magnetic –east–west component vector data and the ion velocity Vx data exhibited anomaly perturbations near the epicenter and the magnetic conjugate area on 11 May and 16 May. (4) The anomaly perturbations in the thermal infrared TBB data, CSES electric field, and magnetic field data all occurred within a consistent perturbation time period and spatial proximity. We also conducted an investigation into the timing, location, and potential causes of the anomaly perturbations using the Vx ion velocity data with magnetic field –east–west component vector data, as well as the horizontal –north–south and vertical component PSD data of the electric field with the magnetic field –east–west component vector data. There may be both chemical and electromagnetic wave propagation models for the “lithosphere—atmosphere—ionosphere” coupling (LAIC) mechanism of the Maduo earthquake.

Previous earthquake polarization (as the ratio of vertical and horizontal components) studies using geomagnetic data were all performed with ground data. The advantage of satellite data is that it is not limited by geography. Therefore, in this work, we tried to select 12 typical earthquakes in Northeast Asia with Ms > 5.0 and an epicenter depth ≤ 40 km within the longitude 105° E–145° E and latitude 38° N–58° N ranges from December 2018 to January 2023 for analysis by using the satellite data of the high-precision magnetometer (HPM) payload onboard the China Seismo-Electromagnetic Satellite (CSES) for the first time in a quiet magnetic environment. The geomagnetic three-component vector data were investigated, and the minimum study period was divided into 10 s intervals. Fourier transform was performed to obtain 0.01–0.2 Hz geomagnetic three-component dynamic spectra, and the time series of the polarization (as the ratio of vertical and horizontal components) data was then obtained. The average value of the polarization data over four years was used to obtain the time series of the polarization perturbation amplitude, after which joint research was conducted. The results showed that (1) earthquakes with larger magnitudes are more likely to exhibit anomaly perturbations recorded by satellites; (2) among all earthquakes with anomalies, the horizontal east–west component perturbation is the largest, the vertical component perturbation is the smallest, and the east–west component may be the dominant component in seismic anomaly observations; (3) the applicability of the polarization method to space-based earthquake-related data is limited; (4) the perturbation amplitude of polarization data can be used as a reference for extracting seismic anomalies; and (5) ion velocity Vx data from the plasma analyzer package (PAP) can be considered to approximately verify the physical mechanism of the anomaly perturbation of the horizontal component in the ionospheric magnetic field, and the two kinds of data (PAP and HPM) can be combined in seismic prediction research.