The short-term prediction of earthquakes is an essential issue connected with human life protection and related social and economic matters. Recent papers have provided some evidence of the link between the lithosphere, lower atmosphere, and ionosphere, even though with marginal statistical evidence. The basic coupling is hypothesized as being via the atmospheric gravity wave (AGW)/acoustic wave (AW) channel. In this paper we analyze a scenario of the low latitude earthquake (Mw = 6.9) which occurred in Indonesia on 5 August 2018, through a multi-instrumental approach, using ground and satellites high quality data. As a result, we derive a new analytical lithospheric–atmospheric–ionospheric–magnetospheric coupling model with the aim to provide quantitative indicators to interpret the observations around 6 h before and at the moment of the earthquake occurrence.

On 25 August 2018, a G3-class geomagnetic storm reached the Earth’s magnetosphere, causing a transient rearrangement of the charged particle environment around the planet, which was detected by the High-Energy Particle Detector (HEPD) on board the China Seismo-Electromagnetic Satellite (CSES-01). We found that the count rates of electrons in the MeV range were characterized by a depletion during the storm’s main phase and a clear enhancement during the recovery caused by large substorm activity, with the key role played by auroral processes mapped into the outer belt. A post-storm rate increase was localized at L-shells immediately above ∼3 and mostly driven by non-adiabatic local acceleration caused by possible resonant interaction with low-frequency magnetospheric waves.

In the polar ionosphere, the electric field is characterized by broadband and power law spectral densities at small/short spatio-temporal scales, which support a possible turbulent nature of the electric field fluctuations. Here, we investigate the multifractal character of the full three-dimensional electric field in the polar ionosphere as recorded on board the first Chinese Seismo-Electromagnetic Satellite (CSES-01). The results of our analysis prove a clear different degree of multifractality of the electric field fluctuations approaching either the polar cap trailing edge or the auroral region. The observed differences in the multifractal character are interpreted in terms of the different natures of the particle precipitation in the polar cap and in the auroral region. A possible link between the multifractal nature of electric field fluctuations, parallel to the geomagnetic field, and filamentation of field aligned currents (FACs) is established.

The storm onset on 7 September 2017, triggered several variations in the ionospheric electron density, causing severe phase fluctuations at polar latitudes in both hemispheres. In addition, although quite rare at high latitudes, clear amplitude scintillations were recorded by two Global Navigation Satellite System receivers during the main phase of the storm. This work attempted to investigate the physical mechanisms triggering the observed amplitude scintillations, with the aim of identifying the conditions favoring such events. We investigated the ionospheric background and other conditions that prevailed when the irregularities formed and moved, following a multi-observations approach. Specifically, we combined information from scintillation parameters and recorded by multi-constellation (GPS, GLONASS and Galileo) receivers located at Concordia station (75.10°S, 123.35°E) and SANAE IV base (71.67°S, 2.84°W), with measurements acquired by the Special Sensor Ultraviolet Spectrographic Imager on board the Defense Meteorological Satellite Program satellites, the Super Dual Auroral Radar Network, the Swarm constellation and ground-based magnetometers. Besides confirming the high degree of complexity of the ionospheric dynamics, our multi-instrument observation identified the physical conditions that likely favor the occurrence of amplitude scintillations at high latitudes. Results suggest that the necessary conditions for the observation of this type of scintillation in high-latitude regions are high levels of ionization and a strong variability of plasma dynamics. Both of these conditions are typically featured during high solar activity.

In the last few decades, the efforts of the scientific community to search earthquake signatures in the atmospheric, ionospheric and magnetospheric media have grown rapidly. The increasing amount of good quality data from both ground stations and satellites has allowed for the detections of anomalies with high statistical significance such as ionospheric plasma density perturbations and/or atmospheric temperature and pressure changes. However, the identification of a causal link between the observed anomalies and their possible seismic trigger has so far been prevented by difficulties in the identification of confounders (such as solar and atmospheric activity) and the lack of a global analytical lithospheric–atmospheric–magnetospheric model able to explain (and possibly forecast) any anomalous signal. In order to overcome these problems, we have performed a multi-instrument analysis of a low-latitude seismic event by using high-quality data from both ground bases and satellites and preserving their statistical significance. An earthquake (Mw = 7.2) occurred in the Caribbean region on 14 August 2021 under both solar quiet and fair weather conditions, thus proving an optimal case study to reconstruct the link between the lithosphere, atmosphere, ionosphere, and magnetosphere. The good match between the observations and novel magnetospheric–ionospheric–lithospheric coupling (M.I.L.C.) modeling of the event confirmed that the fault break generated an atmospheric gravity wave that was able to mechanically perturb the ionospheric plasma density, in turn triggering a variation in the magnetospheric field line resonance frequency.

In this work, we aim to characterize the effective scale height at the ionosphere F2-layer peak (H₀) by using in situ electron density (N_e) observations by Langmuir Probes (LPs) onboard the China Seismo-Electromagnetic Satellite (CSES—01). CSES—01 is a sun-synchronous satellite orbiting at an altitude of ~500 km, with descending and ascending nodes at ~14:00 local time (LT) and ~02:00 LT, respectively. Calibrated CSES—01 LPs N_e observations for the years 2019–2021 provide information in the topside ionosphere, whereas the International Reference Ionosphere model (IRI) provides N_e values at the F2-layer peak altitude for the same time and geographical coordinates as CSES—01. CSES—01 and IRI N_e datasets are used as anchor points to infer H₀ by assuming a linear scale height in the topside representation given by the NeQuick model. COSMIC/FORMOSAT—3 (COSMIC—1) radio occultation (RO) data are used to constrain the vertical gradient of the effective scale height in the topside ionosphere in the linear approximation. With the CSES—01 dataset, we studied the global behavior of H₀ for daytime (~14:00 LT) and nighttime (~02:00 LT) conditions, different seasons, and low solar activity. Results from CSES—01 observations are compared with those obtained through Swarm B satellite N_e-calibrated measurements and validated against those from COSMIC—1 RO for similar diurnal, seasonal, and solar activity conditions. H₀ values modeled by using CSES—01 and Swarm B-calibrated observations during daytime both agree with corresponding values obtained directly from COSMIC—1 RO profiles. Differently, H₀ modeling for nighttime conditions deserves further investigation because values obtained from both CSES—01 and Swarm B-calibrated observations show remarkable and spatially localized differences compared to those obtained through COSMIC—1. Most of the H₀ mismodeling for nighttime conditions can probably to be attributed to a sub-optimal spatial representation of the F2-layer peak density made by the underlying IRI model. For comparison, H₀ values obtained with non-calibrated CSES—01 and Swarm B N_e observations are also calculated and discussed. The methodology developed in this study for the topside effective scale height modeling turns out to be applicable not only to CSES—01 satellite data but to any in situ N_e observation by low-Earth-orbit satellites orbiting in the topside ionosphere.