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.

Earthquakes are the most energetic phenomena in the lithosphere: their study and comprehension are greatly worth doing because of the obvious importance for society. Geosystemics intends to study the Earth system as a whole, looking at the possible couplings among the different geo-layers, i.e., from the earth’s interior to the above atmosphere. It uses specific universal tools to integrate different methods that can be applied to multi-parameter data, often taken on different platforms (e.g., ground, marine or satellite observations). Its main objective is to understand the particular phenomenon of interest from a holistic point of view. Central is the use of entropy, together with other physical quantities that will be introduced case by case. In this paper, we will deal with earthquakes, as final part of a long-term chain of processes involving, not only the interaction between different components of the Earth’s interior but also the coupling of the solid earth with the above neutral or ionized atmosphere, and finally culminating with the main rupture along the fault of concern. Particular emphasis will be given to some Italian seismic sequences.

Nowadays, the possibility that medium-large earthquakes could produce some electromagnetic ionospheric disturbances during their preparatory phase is controversial in the scientific community. Some previous works using satellite data from DEMETER, Swarm and, recently, CSES provided several pieces of evidence supporting the existence of such precursory phenomena in terms of single case studies and statical analyses. In this work, we applied a Worldwide Statistical Correlation approach to M5.5+ shallow earthquakes using the first 8 years of Swarm (i.e., from November 2013 to November 2021) magnetic field and electron density signals in order to improve the significance of previous statistical studies and provide some new results on how earthquake features could influence ionospheric electromagnetic disturbances. We implemented new methodologies based on the hypothesis that the anticipation time of anomalies of larger earthquakes is usually longer than that of anomalies of smaller magnitude. We also considered the signal’s frequency to introduce a new identification criterion for the anomalies. We find that taking into account the frequency can improve the statistical significance (up to 25% for magnetic data and up to 100% for electron density). Furthermore, we noted that the frequency of the Swarm magnetic field signal of possible precursor anomalies seems to slightly increase as the earthquake is approaching. Finally, we checked a possible relationship between the frequency of the detected anomalies and earthquake features. The earthquake focal mechanism seems to have a low or null influence on the frequency of the detected anomalies, while the epicenter location appears to play an important role. In fact, land earthquakes are more likely to be preceded by slower (lower frequency) magnetic field signals, whereas sea seismic events show a higher probability of being preceded by faster (higher frequency) magnetic field signals.

On 20 December 2021, after six quiet years, the Hunga Tonga–Hunga Ha’apai volcano erupted abruptly. Then, on 15 January 2022, the largest eruption produced a plume well registered from satellites and destroyed the volcanic cone previously formed in 2015, connecting the two islands. We applied a multi-parametric and multi-layer study to investigate all the possible pre-eruption signals and effects of this volcanic activity in the lithosphere, atmosphere, and ionosphere. We focused our attention on: (a) seismological features considering the eruption in terms of an earthquake with equivalent energy released in the lithosphere; (b) atmospheric parameters, such as skin and air temperature, outgoing longwave radiation (OLR), cloud cover, relative humidity from climatological datasets; (c) varying magnetic field and electron density observed by ground magnetometers and satellites, even if the event was in the recovery phase of an intense geomagnetic storm. We found different precursors of this unique event in the lithosphere, as well as the effects due to the propagation of acoustic gravity and pressure waves and magnetic and electromagnetic coupling in the form of signals detected by ground stations and satellite data. All these parameters and their detailed investigation confirm the lithosphere–atmosphere–ionosphere coupling (LAIC) models introduced for natural hazards such as volcano eruptions and earthquakes.

The coupling processes among the lithosphere, atmosphere, and ionosphere (LAI) during the earthquake preparation phase are still an open scientific debate. Comprehensive LAI coupling effects around the 2022 Ms6.8 Luding earthquake in China are investigated with a multi-parameter and multi-layer approach, including the b-value, revised accelerated moment release, Earth resistivity, ELF magnetic field emissions, atmospheric electric field, surface temperature, foF2 from ionosonde, GNSS TEC, Ne and magnetic field from CSES and Swarm satellites, and energetic electrons from CSES and NOAA satellites. It is found that the anomalies start from the lithospheric parameters as Earth resistivity and b-values 1–2 years before to reflect the local stress loading in the seismic region, then the ionospheric and atmospheric disturbances occur and accelerate −50 days before and −15 days before, and finally the electrons precipitate a few days before. The simultaneous perturbations in LAI illustrate the thermodynamic coupling channel, such as on 24 August, −12 days before. Meanwhile, the abundant developed ionospheric anomalies without atmospheric disturbances demonstrate the electromagnetic coupling way from the lithosphere to the ionosphere directly. Finally, the results demonstrate a two-way model of LAIC: one way is characterized by a slow chain of processes, of thermodynamic nature, starting from the ground and proceeding to the above atmosphere and ionosphere, showing an exponential trend in the cumulative number of anomalies; the second way is characterized by oscillating electromagnetic coupling between the lithosphere and ionosphere, showing intermittent fluctuations in the corresponding cumulative number of anomalies.

On 6 February 2023, Turkey experienced its most powerful earthquake in over 80 years, with a moment magnitude (Mw) of 7.7. This was then followed by a second earthquake of Mw 7.6 just nine hours later. According to the lithosphere–atmosphere–ionosphere coupling (LAIC) models, such a significant seismic activity is expected to cause anomalies across various observables, from the Earth’s surface to the ionosphere. This multidisciplinary study investigates the preparatory phase of these two major earthquakes by identifying potential precursors across the lithosphere, atmosphere, and ionosphere. Our comprehensive analysis successfully identified and collected various anomalies, revealing that their cumulative occurrence follows an accelerating trend, either exponential or power-law. Most anomalies appeared to progress from the lithosphere upward through the atmosphere to the ionosphere, suggesting a sequential chain of processes across these geospheres. Notably, some anomalies deviated from this overall trend, manifesting as oscillating variations. We propose that these anomalies support a two-way coupling model preceding major earthquakes, highlighting the potential role of fluid chemistry in facilitating these processes.