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.

The Tonga-Kermadec subduction zone represents one of the most active areas from both seismic and volcanic points of view. Recently, two planetary-scale geophysical events took place there: the 2019 M7.2 earthquake (EQ) with the epicentre in Kermadec Islands (New Zealand) and the astonishing 2022 eruption of Hunga Tonga-Hunga Ha’apai (HTHH) volcano. Based on the Lithosphere-Atmosphere-Ionosphere Coupling (LAIC) models, we analysed the three geolayers with a multi-parametric approach to detect any effect on the occasion of the two events, through a comparison aimed at identifying the physics processes that interested phenomena of different nature but in the same tectonic context. For the lithosphere, we conducted a seismic analysis of the sequence culminating with the main shock in Kermadec Islands and the sequence of EQs preceding the HTHH volcanic eruption, in both cases considering the magnitude attributed to the released energy in the lithosphere within the respective Dobrovolsky area. Moving to the above atmosphere, the attention was focused on the parameters—gases, temperature, pressure—possibly influenced by the preparation or the occurrence of the events. Finally, the ionosphere was examined by means of ground and satellite observations, including also magnetic and electric field, finding some interesting anomalous signals in both case studies, in a wide range of temporal and spatial scales. The joint study of the effects seen before, during and after the two events enabled us to clarify the LAIC in this complex context. The observed similarities in the effects of the two geophysical events can be explained by their slightly different manifestations of releasing substantial energy resulting from a shared geodynamic origin. This origin arises from the thermodynamic interplay between a rigid lithosphere and a softer asthenosphere within the Kermadec-Tonga subduction zone, which forms the underlying tectonic context.

On 4 March 2021, a devastating M8.1 earthquake struck the Kermadec Islands of New Zealand. Given the tremendous energy released during the event, we sought to investigate the event’s potential impact on the ionosphere and the inner Van Allen Belt using data from the high-energy electron detectors on board the NOAA-18 satellite. The survey was also extended to the strongest shallow M6.5+ earthquakes occurring between 150° and 190° in longitude, and between −5° and −35° in latitude over the previous ten years. In nearly all cases, evident electron fluxes entering the loss cone were observed. To explore the possibility of a connection between ionospheric signals and tectonic events in this intensely active region, we analyzed electron losses from the inner Van Allen Belt, taking into account latitude, longitude, day/night times, and proximity to the South Atlantic Anomaly. Compared to previous studies, here only the most significant loss phenomena persistent in the ionosphere were considered. Particular interest was reserved for the intense electron loss events that had a duration spanning from a few to several minutes and occurred several hours before and after strong seismic events. Thereafter, time series of electron counting rates and strong Southern Pacific earthquakes were transformed into binary series, and the series multiplication was investigated. The results suggest four peaks of association, including a first couple between electron perturbations detected for ascending semi-orbits and seismic events and a second one between electron perturbations detected in the southern ionosphere and seismic events. They both anticipated the occurrence of earthquakes, occurring around 4 h before them. Other couples were observed between electron perturbations detected for descending semi-orbits and seismic events and between electron perturbations detected in the northern ionosphere and seismic events. They both occurred around 3 h after the occurrence of earthquakes. The case of perturbations anticipating seismic events has the intriguing properties of sustaining the hypothesis that a physical interaction occurred around 6 h before seismic events as in the West Pacific case. A physical model of electrons detected far several thousands of km from the earthquake epicenters was also presented. However, a simulation of random seismic events suggested that the null hypothesis cannot be fully rejected for these associations, prompting many more analyses and case studies.