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.

Seismogenic Quasi-Stationary Electric Fields and Currents from Large-Scale Sources on the Earth's Surface: Comparison of Model Representations

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

Source parameters inversion of the global large earthquakes using 3D SEM Green's functions: Strain Green's function calculation and validation

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

Numerical Modeling on the Propagation of Seismic Electromagnetic Signals

Electromagnetic ULF Fields from an Underground Seismic Source on the Earth Surface and in the Ionosphere

Tomography of the Source Zone of the Great 2011 Tohoku Earthquake

Seismic Electromagnetic Signals (sems) Explained by A Simulation Experiment Using Electromagnetic Waves

The Earth's oscillating electric field (T = 1 day) in relation to the occurrence time of large EQs (Ms>5.0R). A postulated theoretical physical working model and its statistical validation

Electromagnetic Signals Associated With Earthquakes: A Review of Observations, Data Processing, and Mechanisms in China

Physics of Electromagnetic Phenomena Associated with the Rupture of A Finite Fault Model

A Multidimension Source Model for Generating Broadband Ground Motions with Deterministic 3D Numerical Simulations

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

Rethinking earthquake-related DC-ULF electromagnetic phenomena: Towards a physics-based approach

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

Electric current generation associated with the deformation rate of a solid: Preseismic and coseismic signals

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

Can seismogenic atmospheric current influence the ionosphere?

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

On a preseismic electric field "strange attractor" like precursor observed short before large earthquakes

Controlled Analogue Experiments on Propagation of Seismic Electromagnetic Signals