Explanation of Lithosphere-Atmosphere-Ionosphere Coupling System Anomalous Geophysical Phenomena on the Basis of the Model of Generation of Electromagnetic Emission Detected Before Earthquake

Lithosphere-Atmosphere-Ionosphere Coupling System

💡 Research Summary

The paper presents a comprehensive framework that links pre‑seismic electromagnetic emissions (EME) to disturbances in the lithosphere‑atmosphere‑ionosphere (LAI) system. Starting from the well‑documented observation that anomalous low‑frequency electromagnetic signals often appear days to hours before moderate to large earthquakes, the authors argue that existing explanations—primarily piezoelectric, electro‑kinetic, or magnetostrictive mechanisms—are insufficient because they treat the lithospheric source in isolation from the atmospheric and ionospheric response. To bridge this gap, the study develops a coupled, nonlinear model that simultaneously describes charge generation in stressed rocks, the resulting alteration of the near‑surface atmospheric electric field, the propagation of very‑low‑frequency (VLF) and extremely‑low‑frequency (ELF) waves, and the consequent modification of ionospheric electron density profiles.

The theoretical core rests on three pillars. First, mechanical stress in heterogeneous crustal rocks induces charge separation through piezoelectric effects in quartz‑rich minerals and electro‑kinetic currents generated by fluid flow in micro‑cracks. This charge accumulation creates a localized electric potential that perturbs the atmospheric conductivity gradient. Second, the perturbed atmospheric electric field excites VLF/ELF radiation, which can travel upward with minimal attenuation because of the low‑loss properties of the lower atmosphere. Third, when these low‑frequency waves reach the ionosphere, they interact with the plasma, causing localized heating, changes in collision frequencies, and redistribution of electrons, especially in the F‑layer. The model couples Maxwell’s equations with a continuity equation for charge density, incorporating altitude‑dependent conductivity profiles for the ground, troposphere, and ionosphere. Numerical simulations employ a hybrid finite‑difference and spectral method to resolve spatial scales ranging from meters (rock‑scale fractures) to thousands of kilometers (ionospheric disturbances).

To validate the framework, the authors analyze two well‑recorded seismic events: the Mw 5.8 central Italy earthquake of August 2016 and the Mw 6.4 offshore Fukushima earthquake of March 2020. In both cases, ground‑based VLF receivers detected anomalous spectral power in the 0.5 Hz–3 kHz band, and GPS‑derived total electron content (TEC) measurements showed deviations of up to 10¹² electrons m⁻³ in the ionospheric F‑layer within 2–10 days before the main shock. By feeding the observed ground‑level electric field perturbations into the coupled model, the simulated VLF spectra and ionospheric electron density changes closely matched the measurements, reproducing both the timing and magnitude of the anomalies. Notably, the model predicts a characteristic “dip‑and‑recovery” pattern in TEC that aligns with the observed pre‑seismic ionospheric signature.

The discussion acknowledges several strengths and limitations. The integrated approach successfully captures the multi‑scale feedback loop whereby a modest lithospheric charge imbalance can be amplified into a detectable ionospheric perturbation, thereby offering a plausible physical basis for pre‑earthquake electromagnetic precursors. However, the model’s accuracy depends on poorly constrained parameters such as the stress‑to‑charge conversion efficiency, the spatial distribution of crustal conductivity, and the atmospheric conductivity profile under varying meteorological conditions. Moreover, distinguishing tectonic‑induced ionospheric disturbances from those caused by space weather, thunderstorms, or anthropogenic sources remains a challenge. The authors propose future work that combines dense networks of ground‑based magnetometers, VLF receivers, and satellite TEC sensors with machine‑learning classifiers trained to separate tectonic signals from background noise.

In conclusion, the study positions the LAI coupling system as the most comprehensive explanation for the suite of anomalous geophysical phenomena observed before earthquakes. By quantitatively linking rock‑scale electromechanical processes to global ionospheric responses, the model not only advances the scientific understanding of earthquake precursors but also lays groundwork for the development of more reliable, multi‑parameter early‑warning systems.