The study of the preparation phase of large earthquakes is essential to understand the physical processes involved, and potentially useful also to develop a future reliable short-term warning system. Here we analyse electron density and magnetic field data measured by Swarm three-satellite constellation for 4.7 years, to look for possible in-situ ionospheric precursors of large earthquakes to study the interactions between the lithosphere and the above atmosphere and ionosphere, in what is called the Lithosphere-Atmosphere-Ionosphere Coupling (LAIC). We define these anomalies statistically in the whole space-time interval of interest and use a Worldwide Statistical Correlation (WSC) analysis through a superposed epoch approach to study the possible relation with the earthquakes. We find some clear concentrations of electron density and magnetic anomalies from more than two months to some days before the earthquake occurrences. Such anomaly clustering is, in general, statistically significant with respect to homogeneous random simulations, supporting a LAIC during the preparation phase of earthquakes. By investigating different earthquake magnitude ranges, not only do we confirm the well-known Rikitake empirical law between ionospheric anomaly precursor time and earthquake magnitude, but we also give more reliability to the seismic source origin for many of the identified anomalies.

We analyse Swarm satellite magnetic field and electron density data one month before and one month after 12 strong earthquakes that have occurred in the first 2.5 years of Swarm satellite mission lifetime in the Mediterranean region (magnitude M6.1+) or in the rest of the world (M6.7+). The search for anomalies was limited to the area centred at each earthquake epicentre and bounded by a circle that scales with magnitude according to the Dobrovolsky’s radius. We define the magnetic and electron density anomalies statistically in terms of specific thresholds with respect to the same statistical quantity along the whole residual satellite track (|geomagnetic latitude| ≤ 50°, quiet geomagnetic conditions). Once normalized by the analysed satellite tracks, the anomalies associated to all earthquakes resemble a linear dependence with earthquake magnitude, so supporting the statistical correlation with earthquakes and excluding a relationship by chance.

The geomagnetic storm occurred on 25 August 2018, i.e. during the minimum of solar cycle 24, is currently the strongest ever probed by the first China Seismo‐Electromagnetic Satellite (CSES‐01). By integrating the in situ measurements provided by CSES‐01 (orbiting at altitude of 507 km) and by Swarm A satellite (orbiting at ca. 470 km) with ground‐based observations (ionosondes, magnetometers and Global Navigation Satellite System receivers), we investigate the ionospheric response at low‐ and mid‐ latitudes over Brazil. Specifically, we investigate the electrodynamic disturbances driven by solar wind changes, by focusing on the disturbances driving modifications of the Equatorial Electrojet (EEJ). Our proposed multi‐sensor technique analysis mainly highlights the variations in the topside and bottomside ionosphere, and the interplay between Prompt Penetrating Electric Fields and Disturbance Dynamo Electric Fields resulting in EEJ variations. Thanks to this approach and leveraging on the newly available CSES‐01 data, we complement and extend what recently investigated in the Western South American sector, by highlighting the significant longitudinal differences, which mainly come from the occurrence of a daytime counter‐EEJ during both 25 and 26 August at Brazilian longitudes and during part of 26 August only in the Peruvian sector. In addition, the increased thermospheric circulation driven by the storm has an impact on the EEJ during the recovery phase of the storm. The observations at the CSES‐01/Swarm altitudes integrated with the ground‐based observation recorded signatures of Equatorial Ionospheric Anomaly (EIA) crests formation and modification during daytime coupled with the positive ionospheric storm effects at mid‐latitude.

Potential earthquake precursors include, among others, electromagnetic fields, gas emissions, Land Surface Temperature (LST), Sea Surface Temperature (SST), and Surface Air Temperature (SAT) anomalies. These observables have been individually studied, before earthquakes, by many researchers. The ionospheric studies concerning earthquakes (EQs) using magnetic data from Low Earth Orbit (LEO) satellites are increasingly being used to detect ionospheric anomalies before large EQs. Also, LST, SST, and SAT values retrieved from Moderate Resolution Imaging Spectroradiometer (MODIS) Terra and Aqua satellites and Modern-Era Retrospective analysis for Research and Applications Version 2 (MERRA-2) are considered as physical precursors before EQs. In this work, we jointly analyze magnetic, MODIS, and MERRA-2 data in space and time around the epicenters before the selected EQs in Mexico, Japan, Chile, and Indonesia. Our analyses present interesting findings where anomalies in temperature and magnetic field, preceding the considered EQs, are confirmed through different methods. Particularly, we utilize the Fast Fourier Transform (FFT) and the Discrete Cosine Transform (DCT) for analyzing magnetic data over the designated EQs regions. We use the magnetic data acquired by Swarm satellites in the top side ionosphere along with MODIS and MERRA-2. Five case studies are described to prove the effectiveness of our analyses. Precursory anomalies were observed using these methods in different anomalous days from the considered four regions of interest around the epicenter. It is concluded that these methods could be effective and reliable in detecting anomalies preceding the upcoming EQs.