Abstract:For a shallow or ground surface source and receiver at the same level or close depth, it is very difficult or computationally inefficient to simulate seismoelectric or electroseismic wave-fields in stratified porous media by current reflectivity methods, such as the Luco-Apsel-Chen generalized reflection and transmission method (LAC GRTM). In this work, the peak-trough averaging method which has been proved effective and efficient in dealing with this kind of computational problem is adopted to update the seismoelectric and electroseismic modeling algorithm based on LAC GRTM. After thoroughly verifying the accuracy and computational efficiency of the updated algorithm, we utilize it to numerically investigate both the electroseismic and seismoelectric couplings. Snapshots of electroseismic wave-fields indicate evanescent electroseismic conversion, a reverse process of evanescent seismoelectric conversion, dominates at relatively larger ratios of seismic wavelength to interface depth, whereas the interfacial radiation electroseismic conversion is more prominent for the opposite situation. Our seismoelectric modeling results demonstrate that electric signals can arrive at the ground surface a few milliseconds earlier than their causative seismic signals due to evanescent seismoelectric conversion. This is the first modeling result considering source-receiver geometries on the surface capable of explaining similar phenomena reported in geophysical field observations of seismically induced electrokinetic effects over a long history. The updated algorithm offers an accurate and efficient tool for forward modeling and will benefit interpretations of field observations as well as future inversion studies. Plain Language Summary For the case of active seismoelectric or electroseismic measurements with source and receivers located on or near the ground surface, it is computationally difficult for current reflectivity methods, such as Luco-Apsel-Chen generalized reflection and transmission method (LAC GRTM), to simulate the wave-fields in stratified porous media. In this study, the peak-trough averaging method is adopted to update the LAC-GRTM-based seismoelectric and electroseismic modeling algorithms to solve this computational problem. We verify the accuracy and efficiency of the updated algorithm and apply it to investigate the electroseismic and seismoelectric couplings. Two kinds of electroseismic conversions were identified by our modellings. When the ratio of seismic wavelength to interface depth is large, evanescent electroseismic conversion will be dominant; otherwise, the interfacial radiation electroseismic conversion will become prominent. In our seismoelectric modellings, it is found electric signals recorded by ground surface receivers could be several milliseconds earlier than their causative seismic signals due to evanescent seismoelectric conversion. This is the first modeling result capable of explaining similar phenomena that have been reported in seismoelectric field measurements over a long history. Generally speaking, the updated algorithm is a useful forward modeling tool, which could benefit interpretations of field data and future inversion studies.

Fast Convergence Algorithm for Earthquake Prediction Using Electromagnetic Fields Excited by SLF/ELF Horizontal Magnetic Dipole and Schumann Resonance

Stochastic Finite‐Fault Simulation of the Ms 7.0 Lushan Earthquake Based on Frequency‐ and Distance‐Dependent Radiation Patterns

A Fast Time Domain Electric Field Integral Equations for Predicting Transient Responses of Dielectric Objects

Numerical Study of Global ELF Electromagnetic Wave Propagation with Respect to Lithosphere–Atmosphere–Ionosphere Coupling

Hybrid Time-Domain Adaptive Integration Method for Simulating Transient Electromagnetic Responses of an Arbitrary Structure with Large Scale

A Refinement and High-Degree Stable Marching-on-in-Degree Method for Wide-Band Scattering from Perfect Conductors

ELF Near-Field Propagation of a Vertical Electric Dipole Due to Lightning Discharges

Strength of the Electric Dipole Source Field in Multilayer Spherical Media

Pade approximant ELF VED fields in the Earth‐Ionosphere cavity

Seismoelectric and Electroseismic Modeling in Stratified Porous Media with a Shallow or Ground Surface Source

Far-field Prediction from Amplitude-Only Near-Field Measurements Using Equivalent Electric Currents

Forward Calculation of 3D Controlled-Source Electromagnetic Responses Based on Joint Application of Secondary Field and Coupled Potential Formulations

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

Numerical Simulation of Coseismic Electromagnetic Fields Associated with Seismic Waves Due to Finite Faulting in Porous Media

A New Numerical Technique for Simulating the Coupled Seismic and Electromagnetic Waves in Layered Porous Media

SLF/ELF Electromagnetic Field of a Horizontal Dipole in the Presence of an Anisotropic Earth-Ionosphere Cavity

Fast Calculation of the Electromagnetic Field by a Vertical Electric Dipole over a Lossy Ground and Its Application in Evaluating the Lightning Radiation Field in the Frequency Domain

Forward calculation of 3D CSEM responses based on joint application of secondary field and coupled potential formulations

Electrokinetic Effect Combined with Surface-Charge Assumption: a Possible Generation Mechanism of Coseismic EM Signals

Numerical Simulation of Seismo-Electromagnetic Fields Associated with a Fault in a Porous Medium

Fast Evaluating the Electromagnetic Fields Generated by a Vertical Electric Dipole over the Lossy Ground Using Sommerfeld Integral Without Truncation