Atmosphere Ionosphere Coupling Research Articles

The Generation of Seismogenic Anomalous Electric Fields in the Lower Atmosphere, and Its Application to Very-High-Frequency and Very-Low-Frequency/Low-Frequency Emissions: A Review

The purpose of this paper is, first of all, to review the previous works on the seismic (or earthquake (EQ)-related) direct current (DC) (or quasi-stationary) electric fields in the lower atmosphere, which is likely to be generated by the conductivity current flowing in the closed atmosphere–ionosphere electric circuit during the preparation phase of an EQ. The current source is electromotive force (EMF) caused by upward convective transport and the gravitational sedimentation of radon and charged aerosols injected into the atmosphere by soil gasses during the course of the intensification of seismic processes. The theoretical calculations predict that pre-EQ DC electric field enhancement in the atmosphere can reach the breakdown value at the altitudes 2–6 km, suggesting the generation of a peculiar seismic-related thundercloud. Then, we propose to apply this theoretical inference to the observational results of seismogenic VHF (very high frequency) and VLF/LF (very low frequency/low frequency) natural radio emissions. The formation of such a peculiar layer initiates numerous chaotic electrical discharges within this region, leading to the generation of VHF electromagnetic radiation. Earlier works on VHF seismogenic radiation performed in Greece have been compared with the theoretical estimates, and showed a good agreement in the frequency range and intensity. The same idea can also be applied, for the first time, to seismogenic VLF/LF lightning discharges, which is completely the same mechanism with conventional cloud-to-ground lightning discharges. In fact, such seismogenic VLF/LF lightning discharges have been observed to appear before an EQ. So, we conclude in this review that both seismogenic VHF radiation and VLF/LF lightning discharges are regarded as indirect evidence of the generation of anomalous electric fields in the lowest atmosphere due to the emanation of radioactive radon and charged aerosols during the preparation phase of EQs. Finally, we have addressed the most fundamental issue of whether VHF and VLF/LF radiation reported in earlier works is either of atmospheric origin (as proposed in this paper) or of lithospheric origin as the result of microfracturing in the EQ fault region, which has long been hypothesized. This paper will raise a question regarding this hypothesis of lithospheric origin by proposing an alternative atmospheric origin outlined in this review. Also, the data on seismogenic electromagnetic radiation and its inference on perturbations in the lower atmosphere will be suggested to be extensively integrated in future lithosphere–atmosphere–ionosphere coupling (LAIC) studies.

Integrated Analysis of Multi-Parameter Precursors to the Fukushima Offshore Earthquake (Mj = 7.3) on 13 February 2021 and Lithosphere–Atmosphere–Ionosphere Coupling Channels

Integrated Analysis of Multi-Parameter Precursors to the Fukushima Offshore Earthquake (Mj = 7.3) on 13 February 2021 and Lithosphere–Atmosphere–Ionosphere Coupling Channels

The Preparation Phase of the 2023 Kahramanmaraş (Turkey) Major Earthquakes from a Multidisciplinary and Comparative Perspective

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.

A PLL-Based Doppler Method Using an SDR-Receiver for Investigation of Seismogenic and Man-Made Disturbances in the Ionosphere

The article describes in detail the equipment and method for measuring the Doppler frequency shift (DFS) on an inclined radio path, based on the principle of the phase-locked loop using an SDR receiver for the investigation of seismogenic and man-made disturbances in the ionosphere. During the two M7.8 earthquakes in Nepal (25 April 2015) and Turkey (6 February 2023), a Doppler ionosonde detected co-seismic and pre-seismic effects in the ionosphere, the appearances of which are connected with the various propagation mechanisms of seismogenic disturbance from the lithosphere up to the ionosphere. One day before the earthquake in Nepal and 90 min prior to the main shock, an increase in the intensity of Doppler bursts was detected, which reflected the disturbance of the ionosphere. A channel of geophysical interaction in the system of lithosphere–atmosphere–ionosphere coupling was traced based on the comprehensive monitoring of the DFS of the ionospheric signal, as well as of the flux of gamma rays in subsoil layers of rocks and in the ground-level atmosphere. The concept of lithosphere–atmosphere–ionosphere coupling, where the key role is assigned to ionization of the atmospheric boundary layer, was confirmed by a retrospective analysis of the DFS records of an ionospheric signal made during underground nuclear explosions at the Semipalatinsk test site. A simple formula for reconstructing the velocity profile of the acoustic pulse from a Dopplerogram was obtained, which depends on only two parameters, one of which is the dimension of length and the other the dimension of time. The reconstructed profiles of the acoustic pulses from the two underground nuclear explosions, which reached the height of the reflection point of the sounding radio wave, are presented.

The Preparation Phase of the 2022 ML 5.7 Offshore Fano (Italy) Earthquake: A Multiparametric–Multilayer Approach

This paper presents an analysis of anomalies detected during the preparatory phase of the 9 November 2022 ML = 5.7 earthquake, occurring approximately 30 km off the coast of the Marche region in the Adriatic Sea (Italy). It was the largest earthquake in Italy in the last 5 years. According to lithosphere–atmosphere–ionosphere coupling (LAIC) models, such earthquake could induce anomalies in various observable variables, from the Earth’s surface to the ionosphere. Therefore, a multiparametric and multilayer approach based on ground and satellite data collected in each geolayer was adopted. This included the revised accelerated moment release method, the identification of anomalies in atmospheric parameters, such as Skin Temperature and Outgoing Longwave Radiation, and ionospheric signals, such as Es and F2 layer parameters from ionosonde measurements, magnetic field from Swarm satellites, and energetic electron precipitations from NOAA satellites. Several anomalies were detected in the days preceding the earthquake, revealing that their cumulative occurrence follows an exponential trend from the ground, progressing towards the upper atmosphere and the ionosphere. This progression of anomalies through different geolayers cannot simply be attributed to chance and is likely associated with the preparation phase of this earthquake, supporting the LAIC approach.

Variability of Ionospheric Total Electron Content Over Morocco During the Godzilla Sand and Dust Storm of June 2020

During sand and dust storm (SDS) events, atmospheric suspension and transport of sand and dust brings a reasonable amount of electrification in the atmosphere which plays a very important role in the atmosphere-ionosphere coupling. The Godzilla SDS began on 5<sup>th</sup> June 2020 in Algeria following a decrease in pressure and spread to other areas across the Sahara between 6<sup>th</sup> and 28<sup>th</sup> June 2020. Using SDS data from Copernicus Sentinel-5P satellite mission and Vertical Total Electron Content (VTEC) data from four GNSS receiver stations: IFR1 (Ifrane Seismic), MELI (Melilla), TETN (Tetouan) and OUCA (Ouca) over Morocco, we investigate the possible ionospheric TEC variability over the four GNSS receiver stations during the Godzilla SDS event which was tracked using the Sentinel-5P Satellite mission. Solar wind parameters: Horizontal component of Interplanetary Magnetic Field (IMF-Bz), interplanetary Electric Field (IEF-Ey) and solar wind speed (V) and geomagnetic indices: Disturbance Storm Time (Dst) and Planetary K (Kp) indices were examined and showed very minimal geomagnetic influence during the period. We observed major ionospheric disturbances over the four Global Navigation Satellite System (GNSS) receiver stations on 16<sup>th</sup>, 17<sup>th</sup>, 18<sup>th</sup>, 21<sup>st</sup>, 22<sup>nd</sup>, 23<sup>rd</sup> 25<sup>th</sup> and 26<sup>th</sup> June 2020: the period with the Sentinel-5P Aerosol Index (SAI) of more than 4 as recorded by the Sentinel-5P Satellite engine. The daily VTEC values over the four GNSS receiver stations recorded continuous electron density perturbations during these days. Apart from the ionospheric TEC perturbations, significant enhancements and decreases in daily maximum VTEC values over the four GNSS receiver stations were also noted. These were attributed to the changes in the atmospheric electric fields generated by the SDS event. The VTEC plots for each day exhibited similar trends, hence exhibited the same ionospheric dynamics. VTEC depletions of depths 3 to 6 TECU over all the four GNSS receiver stations were noted on 12<sup>th</sup>, 14<sup>th</sup>, 17<sup>th</sup>, 20<sup>th</sup> and 25<sup>th</sup> June 2020. Nighttime VTEC enhancements were also noted and majorly occurred between 20:00 and 21:00 UT on 9<sup>th</sup>, 13<sup>th</sup>, 15<sup>th</sup>, 17<sup>th</sup>, 19<sup>th</sup>, 20<sup>th</sup> and 21<sup>st</sup> June 2020. This was attributed to the development of the electron avalanche processes including dust and electron absorption or losses and the active conversion to electron dissociative attachment leading to electron excitation. In conclusion, the Godzilla SDS of June 2020 led to the electron density perturbations over Morocco.

Pre-Earthquake Oscillating and Accelerating Patterns in the Lithosphere–Atmosphere–Ionosphere Coupling (LAIC) before the 2022 Luding (China) Ms6.8 Earthquake

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.

Editorial of Special Issue “Remote Sensing Observations to Improve Knowledge of Lithosphere–Atmosphere–Ionosphere Coupling during the Preparatory Phase of Earthquakes”

We launched this Special Issue with the aim of collecting papers that use satellite data and new methodologies to understand the preparatory phase of medium–large earthquakes in the world [...]

Thermal Anomalies Observed during the Crete Earthquake on 27 September 2021

This study examines the response of the thermal channel within the Lithosphere–Atmosphere–Ionosphere Coupling (LAIC) mechanism during the notable earthquake in Crete, Greece, on 27 September 2021. We analyze spatio-temporal profiles of Surface Latent Heat Flux (SLHF), Outgoing Longwave Radiation (OLR), and Atmospheric Chemical Potential (ACP) using reanalysis data from the National Oceanic and Atmospheric Administration (NOAA) satellite. Anomalies in these parameters are computed by removing the background profile for a non-seismic condition. Our findings reveal a substantial anomalous increase in these parameters near the earthquake’s epicenter 3 to 7 days before the main shock. The implications of these observations contribute to a deeper understanding of the LAIC mechanism’s thermal channel in seismic events.

The Analysis of Ionospheric TEC Anomalies Prior to the Jiuzhaigou Ms7.0 Earthquake Based on BeiDou GEO Satellite Data

The position between BeiDou geostationary Earth orbit (GEO) satellites and ground-based receiving stations can roughly be considered to be constant with negligible fluctuations; thus, the total electron content (TEC) data over a fixed ionospheric piercing point (IPP) can be continuously acquired, which is advantageous for monitoring ionospheric disturbances. Focused on the Jiuzhaigou Ms7.0 earthquake that occurred on 8 August 2017, the TEC data inverted by the BeiDou GEO satellite were analyzed to extract ionospheric disturbances potentially associated with the earthquake. It was found that significant anomalies in ionospheric TEC occurred 10–11 days, 6–7 days, and 1–9 h prior to the earthquake, which was mainly located in the southeast and southwest directions within about 2500 km distance from the epicenter. Comparing the spatial and temporal characteristics between the ionospheric disturbance and the radon gas near the surface, the atmospheric electric field, and the spectrum of TEC data, it was considered that the chemical and acoustic–gravity wave pathway may play an important role in the lithosphere–atmosphere–ionosphere coupling (LAIC) mechanism.

Synchronized and Co-Located Ionospheric and Atmospheric Anomalies Associated with the 2023 Mw 7.8 Turkey Earthquake

Earth observations from remotely sensed data have a substantial impact on natural hazard surveillance, specifically for earthquakes. The rapid emergence of diverse earthquake precursors has led to the exploration of different methodologies and datasets from various satellites to understand and address the complex nature of earthquake precursors. This study presents a novel technique to detect the ionospheric and atmospheric precursors using machine learning (ML). We examine the multiple precursors of different spatiotemporal nature from satellites in the ionosphere and atmosphere related to the Turkey earthquake on 6 February 2023 (Mw 7.8), in the form of total electron content (TEC), land surface temperature (LST), sea surface temperature (SST), air pressure (AP), relative humidity (RH), outgoing longwave radiation (OLR), and air temperature (AT). As a confutation analysis, we also statistically observe datasets of atmospheric parameters for the years 2021 and 2022 in the same epicentral region and time period as the 2023 Turkey earthquake. Moreover, the aim of this study is to find a synchronized and co-located window of possible earthquake anomalies by providing more evidence with standard deviation (STDEV) and nonlinear autoregressive network with exogenous inputs (NARX) models. It is noteworthy that both the statistical and ML methods demonstrate abnormal fluctuations as precursors within 6 to 7 days before the impending earthquake over the epicenter. Furthermore, the geomagnetic anomalies in the ionosphere are detected on the ninth day after the earthquake (Kp > 4; Dst < −70 nT; ap > 50 nT). This study indicates the relevance of using multiple earthquake precursors in a synchronized window from ML methods to support the lithosphere–atmosphere–ionosphere coupling (LAIC) phenomenon.

Lower Ionospheric variations during the intense tectonic activity in the broader area of Greece on October of 2020

In this paper, we investigate the Lower ionospheric variations from TEC observations during the intense seismic activity of October 2020 in the area of Greece. The data were analysed using both, statistical analysis of TEC variations in order to detect uneven gross variations and Discrete Fourier analysis in order to investigate the TEC turbulence. The results of this investigation indicate that the High-Frequency limit fo of the ionospheric turbulence content, increases as aproaching the occurrence time of the earthquake, pointing to the earthquake epicenter, in accordance with our previous investigations. We conclude that the Lithosphere Atmosphere Ionosphere Coupling, LAIC, mechanism through acoustic or gravity waves could explain this phenomenology. In addition, the statistical analysis shows an excess greater than 3𝜎 from the mean TEC values, one and seven days before the earthquake. Since no major disturbance of the geomagnetic field occured during these days, we conclude that we probably observed precursory ionospheric variations in accordance to analogous findings from the variation of VH/VHF electromagnetic wave propagrations over strong earthquake areas.

Automatic Identification and Statistical Analysis of Data Steps in Electric Field Measurements from CSES-01 Satellite

The spaceborne Electric Field Detector (EFD) is one of the payloads of the China Seismo-Electromagnetic Satellite (CSES-01), which can measure electric field data at near-Earth orbit for investigating fundamental scientific topics such as the dynamics of the top-side ionosphere, lithosphere–atmosphere–ionosphere coupling, and electromagnetic field emissions possibly associated with earthquake occurrence. The Extremely Low-Frequency (ELF) waveform shows anomalous step variations, and this work proposes an automatic detection algorithm to identify steps and analyze their characteristics using a convolutional neural network. The experimental results show that the developed detection method is effective, and the identification performance reaches over 90% in terms of both accuracy and area under the curve index. We also analyze the rate of the occurrence of steps in the three components of the electric field. Finally, we discuss the stability of the statistical results on steps and their relevance to the probe’s function. The research results provide a guideline for improving the quality of EFD data, and further applications in monitoring the low-Earth electromagnetic environment.

The Correlation between Ionospheric Electron Density Variations Derived from Swarm Satellite Observations and Seismic Activity at the Australian–Pacific Tectonic Plate Boundary

Swarm electron density (Ne) observations from the Langmuir probe (LP) can detect ionospheric disturbances at the altitude of a satellite. Along-track satellite observations provide a large number of very short observations of different places in the ionosphere, where Ne is disturbed. Moreover, different perturbations occupy various Ne signal frequencies. Therefore, such short signals are more recognizable in two dimensions, where aside from their change in time, we can observe their diversity in the frequency domain. Spectral analysis is an essential tool applied here, as it enables signal decomposition and the recognition of composite patterns of Ne disturbances that occupy different frequencies. This study shows a high-resolution application of short-term Fourier transform (STFT) to Swarm Ne observations in the Papua New Guinea region in the vicinity of earthquakes, tsunamis, and related general seismic activity. The system of tectonic plate junctions, including the Pacific–Australian boundary, is located orthogonally to Swarm track footprints. The selected wavelengths of seismically induced ionospheric disturbances detected via Swarm are compared with the three sets of three-month records of seismic activity: in the winter solstice of 2016/2017, when seismic activity was highest, and in the summer solstice and vernal equinox of 2016, which were calmer. Moreover, more Swarm data records are analyzed at the same latitudes for validation purposes, in a place where there are no tectonic plate boundaries that are orthogonal to the Swarm orbital footprint. Additional validation is supplied through Swarm Ne observations from completely different latitudes, where the Swarm orbital footprint orthogonally crosses a different subducting plate boundary. Aside from the seismic energy, the solar radio flux (F10.7), equatorial plasma bubbles (EPBs), and geomagnetic ap and Dst indices are also reviewed here. Their influence on the ionospheric Ne is also found in Swarm observations. Finally, the Pearson correlation coefficient (PCC), applied to the pairs of 3-month time series created from Swarm Ne variations, seismic energy, ap, Dst, and F10.7, summarizes the graphical inspection of mutual correlations. It points to the predominant correlation of Swarm Ne disturbances with seismicity, especially during nighttime. We show that most of the Ne disturbances at a selected wavelength of 300 km correlate more with seismicity than with geomagnetic and solar indices. Therefore, Swarm LP can be assessed as being capable of observing the lithosphere–atmosphere–ionosphere coupling (LAIC) from the orbit.

Analysis of Ocean–Lithosphere–Atmosphere–Ionosphere Coupling Related to Two Strong Earthquakes Occurring in June–September 2022 on the Sea Coast of Philippines and Papua New Guinea

Scientific progress in the context of seismic precursors reveals a systematic mechanism, namely lithosphere–atmosphere–ionosphere coupling (LAIC), to elaborate the underlying physical processes related to earthquake preparation phases. In this study, a comprehensive analysis was conducted for two earthquakes that occurred on the sea coast through tidal force fluctuation to investigate ocean–lithosphere–atmosphere–ionosphere coupling (OLAIC), based on oceanic parameters (i.e., sea potential temperature and seawater salinity), air temperature and electron density profiles. The interrupted enhancement and diffusion process of thermal anomalies indicate that the intensity of seismic anomalies in the atmosphere is affected by the extent of land near the epicenter. By observing the evolution of the ocean interior, we found that the deep water was lifted and formed upwelling, which then diffused along the direction of plate boundaries with an “intensification-peak-weakening” trend under the action of the accelerated subduction of tectonic plates. Furthermore, the analysis shows that the seismic anomalies have two propagation paths: (i) along active faults, with the surface temperature rising as the initial performance, then the air pressure gradient being generated, and finally the ionosphere being disturbed; (ii) along plate boundaries, upwelling, which is the initial manifestation, leading to changes in the parameters of the upper ocean. The results presented in this study can contribute to understanding the intrinsic characteristics of OLAIC.

The Analysis of Lithosphere–Atmosphere–Ionosphere Coupling Associated with the 2022 Luding Ms6.8 Earthquake

Taking the Luding Ms6.8 earthquake (EQ) on 5 September 2022 as a case study, we investigated the potential seismic anomalies of the ionosphere, infrared radiation, atmospheric electrostatic field (AEF), and hot spring ions in the seismogenic region. Firstly, we analyzed the multi-parameter anomalies in the ionosphere around the epicenter and found synchronous anomalous disturbances in the ground parameters, namely the global ionospheric map (GIM), GPS, TEC, and satellite parameters, such as the He+ and O+ densities on 26 August under relatively quiet solar–geomagnetic conditions (F10.7 < 120 SFU; Kp < 3; Dst > −30 nT; |AE| < 500 nT). Next, both the anomaly analysis of the infrared radiation and AEF, and the survey results of the Luding EQ scientific expedition on the hot spring ions showed pre-seismic anomalous variations at different time periods in the seismogenic region. The characteristics of Earth’s multi-sphere coupling anomalies in temporal evolution and spatial distribution were obvious, which validated the Lithosphere–Atmosphere–Ionosphere Coupling (LAIC) mechanism. Finally, combining the analysis results and the LAIC mechanism, we suggested that the multi-sphere coupling anomalies were more likely associated with the Luding Ms6.8 EQ, and that the differential motion and the regional crustal stress accumulation between the Chuandian block and the Bayan Har block might have led to this EQ. Furthermore, remote sensing and ground-based monitoring technologies can play an important role in corroborating and compensating each other, while further study of the multi-sphere coupling mechanism will provide a clearer understanding of the seismogenic process for major EQs.

Observation of Ultra-Low-Frequency Wave Effects in Possible Association with the Fukushima Earthquake on 21 November 2016, and Lithosphere–Atmosphere–Ionosphere Coupling

The study presents seismogenic ULF (ultra-low-frequency) wave effects, as observed at our own new magnetic observatory at Asahi (geographic coordinates: 35.770° N, 140.695° E) in Chiba Prefecture. Our target earthquake (EQ) is a huge one offshore of Fukushima prefecture (37.353° N, 141.603° E) with a magnitude (M) of 7.4, which occurred at 20.59 h on November 21 UT, 2016. As a sampling frequency of 1 Hz was chosen for our induction magnetometer, we could detect both ULF wave effects: ULF radiation from the lithosphere, and the ULF depression effect, indicative of lower ionospheric perturbations. Observing the results of polarization analyses, we detected clear enhancements in ULF (frequency = 0.01–0.03 Hz) lithospheric radiation 14 days, 5 days, and 1 day before the EQ, and also observed a very obvious phenomenon of ULF (0.01–0.03 Hz) depression just 1 day prior to the EQ, which is regarded as the signature of lower ionospheric perturbations. These findings suggest that pre-EQ seismic activity must be present in the lithosphere, and also that the lower ionosphere was very much perturbed by the precursory effects of the Fukushima EQ. These new observational effects from our station have been compared with our previous investigations on different seismogenic topics for the same EQ, including the ULF observations at another magnetic observatory at Kakioka, belonging to the Japan Meteorological Agency (JMA), about 50 km north of our Asahi station, subionospheric VLF/LF propagation data (Japanese and Russian data), AGW (Atmospheric gravity wave) activity in the stratosphere, and satellite observation of particle precipitations. We have found that seismogenic anomalies of different parameters tend to happen just around the EQ day, but mainly before the EQ, and have found the chain-like tendency of the effects of the lithosphere, which seem to propagate upwards the lower ionosphere. Finally, we will try to gain a better understanding of the physical phenomena or mechanisms of the lithosphere–atmosphere–ionosphere coupling (LAIC) process during the EQ preparation phase.

Energetic Particle Precipitation Influence on Tidal Variations of Thermosphere Parameters in September 2017

The role of energetic particle precipitation in the formation of thermospheric tides is investigated. Using the Global Self-consistent Model of the Thermosphere, Ionosphere, and Protonosphere (GSM TIP) and two models of energetic particle precipitation, APM (Auroral Precipitation Model) and AIMOS 2.0.-AISstorm (Atmospheric Ionization Module Osnabrück 2.0—Atmospheric Ionization during Substorms), we performed simulations for the period 3–12 September 2017. This period covers both geomagnetically quiet days and the storm of 7–8 September. The analysis shows that migrating tides diurnal wave (DW) and semidiurnal wave (SW) prevail on quiet days for both versions of the simulations. On the day of maximum storm development on 8 September and the day after—9 September, the nonmigrating components of tidal waves, mainly DW0 and SW0, are intensified. There are also increasing differences in the spatial structure of tides between the two versions of simulations, especially between diurnal and semidiurnal tides at 154 km. On the disturbance days, the amplitudes of diurnal and semidiurnal tides at 154 km in the APM version are markedly greater than the corresponding values in the AIMOS version. Thus, it should be noted that the integral characteristics of the ionization function from precipitations are more important for the formation of DW and SW than its spatial structure or the features of temporal variations. A comparison of the total electron content (TEC) maps for the two versions of the simulation showed that the AIMOS version nicely reproduces the experimental data for a quiet time. The APM version is less accurate for quiet time but simulates the quantitative increase to disturbed conditions at high and middle latitudes better. For model reproduction of observed TEC variations at low and equatorial latitudes, it is not enough to consider the thermospheric source of thermospheric tides. In this case, the role of atmosphere–ionosphere coupling is very important.

Are There One or More Geophysical Coupling Mechanisms before Earthquakes? The Case Study of Lushan (China) 2013

Several possible lithosphere–atmosphere–ionosphere coupling mechanisms before earthquake occurrence are presented in the literature. They are described by several models with different interaction channels (e.g., electromagnetic, mechanics, chemical, thermal), sometimes in conflict with each other. In this paper, we search for anomalies six months before the Lushan (China) 2013 earthquake in the three geo-layers looking for a possible view of the couplings and testing if one or another is more reliable to describe the observations. The Lushan earthquake occurred in China’s Sichuan province on 20 April 2013, with a magnitude of Mw = 6.7. Despite the moderate magnitude of the event, it caused concern because its source was localized on the southwest side of the same fault that produced the catastrophic Wenchuan event in 2008. This paper applies a geophysical multi-layer approach to search for possible pre-earthquake anomalies in the lithosphere, atmosphere, and ionosphere. In detail, six main increases in the accumulated seismic stress were depicted. Anomalous geomagnetic pulsations were recorded in the Chengdu observatory, sometimes following the increased stress. Atmosphere status and composition were found to be anomalous in several periods before the earthquake, and, spatially, the anomalies seem to appear firstly far from the upcoming earthquakes and later approaching the Longmenshan fault where the Lushan earthquakes nucleated. The Formosat-3 data identified interesting anomalies in the altitude or electron content of the ionospheric F2 peak in correspondence with seismic and atmospheric anomalies 130 days before the earthquake. In addition, the total electron content showed high anomalous values from 12 to 6 days before the earthquake. We compared the anomalies and tried to explain their correspondences in different geo-layers by the lithosphere–atmosphere–ionosphere coupling models. In particular, we identified three possible couplings with different mechanisms: a first, about 130 days before the earthquake, with a fast (order of one day) propagation delay; a second, about 40 days before the earthquake occurrence, with a propagation delay of few days and a third from 2.5 weeks until one week before the event. Such evidence suggests that the geo-layers could interact with different channels (pure electromagnetic or a chain of physical-chemical processes) with specific propagation delays. Such results support the understanding of the preparation for medium and large earthquakes globally, which is necessary (although not sufficient) knowledge in order to mitigate their impact on human life.

Scrutinizing and rooting the multiple anomalies of Nepal earthquake sequence in 2015 with the deviation–time–space criterion and homologous lithosphere–coversphere–atmosphere–ionosphere coupling physics

Abstract. The continuous increasing of Earth observations benefits geosciences and seismicity study but increases greatly the difficulties in understanding and discriminating multiple source data. Although the lithosphere–coversphere–atmosphere-ionosphere (LCAI) coupling paradigm and the deviation–time–space (DTS) criterion were presented for better searching for and understanding the potential seismic anomalies from multiple observations, the strict consistency of spatiotemporal characteristics and homologous physics of multiple-parameter seismic anomalies has not been investigated sufficiently. With the 2015 Nepal earthquake sequence being a typical case, the reported multi-parameter anomalies were systematically reviewed, and their space–time characteristics were summarized thoroughly in this study. Numerical simulation with refined geological structures in three-dimensional space revealed the inhomogeneous crustal stress field alteration (CSFA) along the faults and around the hypocenters of the 2015 Nepal earthquake sequence, which is expected to be the root of the seismic anomalies. The stress-activated positive charge carriers would have given rise to different responses near the ground surface (coversphere), including the microwave dielectric reduction, the additional infrared radiation, and the atmospheric ionization, which subsequently affected the physical properties of the atmosphere and the ionosphere and resulted in abnormal phenomena therein. Based on the DTS criterion and LCAI coupling paradigm, the seismic anomalies of the 2015 Nepal earthquakes were scrutinized strictly, and the retained anomalies were rooted carefully to the regional CSFA as well as its local blocking. Therefore, an integrated LCAI coupling framework with strict space–time correspondence and homologous physics in CSFA was proposed for the 2015 Nepal earthquake sequence. This research provides a definite philosophy as well as a practical solution for scrutinizing the rootable seismic anomalies from multi-parameter observations of earthquakes, which is of scientific meanings for searching earthquake precursors and reaching earthquake prediction.