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.

In this work, we present 16 directional and 27 intensity high‐quality values from Iberia. Moreover, we have updated the Iberian archeomagnetic catalogue published more than 10 years ago with a considerable increase in the database. This has led to a notable improvement of both temporal and spatial data distribution. A full vector paleosecular variation curve from 1000 BC to 1900 AD has been developed using high‐quality data within a radius of 900 km from Madrid. A hierarchical bootstrap method has been followed for the computation of the curves. The most remarkable feature of the new curves is a notable intensity maximum of about 80 μT around 600 BC, which has not been previously reported for the Iberian Peninsula. We have also analyzed the evolution of the paleofield in Europe for the last three thousand years and conclude that the high maximum intensity values observed around 600 BC in the Iberian Peninsula could respond to the same feature as the Levantine Iron Age Anomaly, after travelling westward through Europe.

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.