Lithosphere Ionosphere Coupling Associated with Seismic Swarm in the Balkan Peninsula from ROB-TEC and GPS

Abstract

With the rapid development of global navigation satellite systems (GNSS) and their increasingly wide range of applications in atmospheric science, total electron content (TEC) data are widely used in the theoretical study of layer coupling related to seismicity. This study detected and analyzed pre-earthquake ionospheric anomalies (PEIA) by using TEC data from the Royal Observatory of Belgium (ROB), and analyzed coseismic ionospheric disturbance (CID) with vertical TEC (VTEC) from the GPS stations in earthquake preparation areas. The results show that PEIA appear to increase continuously from 08:00–12:00 UT in the 3 days before a seismic swarm of Mw > 5.0. The ionosphere over the seismogenic zones exhibited large-scale anomalies when multiple seismogenic zones of the Balkan Peninsula spatially and temporally overlapped. Moreover, the TEC around the earthquake centers showed a positive anomaly lasting for 7 h. In a single seismogenic zone in Greece, the TEC around the earthquake center reached over +3.42 TECu. In addition, the CID observed from GPS stations shows that with the increase in the number of earthquakes, the ionosphere over the seismogenic area is more obviously disturbed, and after three strong earthquakes, TEC suddenly decreased over the seismogenic area and formed a phenomenon similar to an ionospheric hole. We conclude that a lithosphere–atmosphere–ionosphere coupling mechanism existed before the seismic swarm appeared in the Balkan Peninsula. Earthquake-induced VTEC anomalies occurred more frequently within a 3–10 day window before the earthquake. This phenomenon is particularly evident when multiple seismogenic zones overlap spatiotemporally.