Abstract:Taking into account the motional induction effect, in which the Earth crust that has finite electrical conductivity vibrates in the ambient geomagnetic field resulting in motionally induced electric current, we derive semi-analytical solutions of seismo-electromagnetic signals generated by an earthquake source in 3-D multi-layered media, which consists of an air half-space and multiple solid layers. First, both the elastic and electromagnetic (EM) wave-fields involved in the governing equations, which have the form of Maxwell’s equations coupled with elastodynamic equations, are expanded by a set of vector basis functions in cylindrical coordinate system. Then, we reorganize the transformed governing equations expressed by expansion coefficients and obtain corresponding first-order ordinary differential equations for the wave-fields in air and solid media. The expansion of the motionally induced electric current and the reorganization of Maxwell’s equations are the most important part, and also the most complicated and tedious part of this work. Thereafter, we solve the first-order ordinary differential equations through the Luco–Apsel–Chen generalized reflection and transmission method gaining solutions of the expansion coefficients. Finally, we obtain the frequency–space-domain semi-analytical solutions written as integrations of corresponding expansion coefficients over wavenumber domain, which can be numerically calculated by the discrete wavenumber method. The time-domain solutions can be achieved by further applying the discrete inverse Fourier transform. To have a numerical stability at any high frequency, we adopt the analytical regularization approach in the derivation process by introducing two artificial interfaces with infinitely small distance from the source. On the basis of the semi-analytical solutions, we can tell that only EM fields of TM mode (in which magnetic fields are transversely polarized) will be induced by SH waves, whereas EM fields of both TE mode (in which magnetic fields are transversely polarized) and TM mode will be induced by P and SV waves. The derived semi-analytical solutions can be used to calculate seismo-electromagnetic signals either below or above the free surface.

Seismic Wave Fields in a Spherically Symmetric Earth: An Analytical Solution

Modelling Seismic Wave Propagation in a Two-Dimensional Cylindrical Whole-Earth Model Using the Pseudospectral Method

Complete Synthetic Seismograms Based on A Spherical Self-Gravitating Earth Model with an Atmosphere-Ocean-Mantle-Core Structure

Analytic and numerical solutions to the seismic wave equation in continuous media

An Efficient Approach of the Pseudospectral Method for Modelling of Geometrically Symmetric Seismic Wavefield

Seismic stress path induced by obliquely incident waves and its simulation

SH-wavefield Simulation for a Laterally Heterogeneous Whole-Earth Model Using the Pseudospectral Method

Three-dimensional Anisotropic Seismic Wave Modelling in Spherical Coordinates by a Collocated-Grid Finite-Difference Method

An Analytical Solution to the Scattering of Cylindrical SH Waves by a Partially Filled Semi-Circular Alluvial Valley: Near-Source Site Effects

A Two‐dimensional Hybrid Method for Modeling Seismic Wave Propagation in Anisotropic Media

SIMULATING SEISMIC WAVE PROPAGATION IN TRANSVERSELY ISOTROPIC MEDIA BY USING PSEUDO-SPECTRAL METHOD

Symplectic Stereomodelling Method for Solving Elastic Wave Equations in Porous Media

A 3D Optimized Frequency–Wavenumber (FK), Time–Space Optimized Symplectic (TSOS) Hybrid Method for Teleseismic Wave Modeling

Soil seismic analysis for 2D oblique incident waves using exact free-field responses by frequency-based finite/infinite element method

Mathematical Modeling of Seismic Wave Kinematics in Complex Media

Quasi‐spherical Approach for Seismic Wave Modeling in a 2‐D Slice of a Global Earth Model with Lateral Heterogeneity

A Low-Dispersive Symplectic Partitioned Runge-Kutta Method for Solving Seismic-Wave Equations: I. Scheme and Theoretical Analysis

Semi-analytical Solutions of Seismo-Electromagnetic Signals Arising from the Motional Induction in 3-D Multi-Layered Media: Part II—numerical Investigations

Seismic wave-field modeling of 3D irregularities in layered half-space under P-, S-, and Rayleigh waves with arbitrary incident azimuth

Near-Field Seismic Motion: Waves, Deformations and Seismic Moment

3D Sensitivity Kernels with Full Attenuation Computed by a Combination of the Strong Stability Preserving Runge-Kutta Method and the Scattering-Integral Method