---
title: 'Seismic: Wave Analysis and Inversion'
url: https://www.emergentmind.com/topics/seismic
type: topic
---

# Seismic: Wave Analysis and Inversion

Seismic denotes the physics, observation, and analysis of wavefields generated by Earth’s vibrations and, in planetary contexts, by analogous vibrations of other bodies. In current research, the term spans earthquake-source characterization from far-field \(P\)- and \(S\)-waves, continuous monitoring with station networks and unconventional sensors, reconstruction and interpretation of incomplete or noisy seismic records, subsurface inversion and semantic annotation, engineered control of near-surface wave propagation, and extrapolation of seismological methods to icy ocean worlds and other planetary bodies [1808.03049] [1703.00561] [2206.01785].

## 1. Source physics, wave types, and inverse formulations

A central seismic problem is the reconstruction of source properties from observed wavefields. One explicit formulation treats the source as a localized tensorial force density,
\[
F_i = M_{ij}\,\partial_j \delta(\mathbf R-\mathbf R_0),
\]
with a symmetric seismic-moment tensor \(M_{ij}\). In a homogeneous isotropic body, the displacement field is decomposed into near-field and far-field parts, and the far field further separates into longitudinal and transverse components,
\[
\mathbf u=\mathbf u^n+\mathbf u^f, \qquad \mathbf u^f=\mathbf u_l^f+\mathbf u_t^f.
\]
The far-field amplitudes yield three observational quantities through the vectors \(\mathbf v_l\) and \(\mathbf v_t\), but a general symmetric moment tensor has six components, so additional physical structure is required for inversion [1808.03049].

One such structure is the Kostrov vectorial representation,
\[
M_{ij}=2\mu S u^0\left(s_i a_j+a_i s_j\right),
\]
where \(\mathbf s\) is the fault-normal direction, \(\mathbf a\) is the slip direction, \(u^0\) is slip magnitude, \(S\) is fault area, and \(\mu\) is shear modulus. This reduces the unknowns to four and makes the tensor traceless, \(M_{ii}=0\). The remaining closure is obtained from energy conservation together with a covariance condition that constrains \(\mathbf s\), \(\mathbf a\), and the source-to-receiver direction \(\mathbf n\) to lie in the same plane. The resulting framework yields an explicit resolved form of the moment tensor,
\[
M_{ij}=\frac{M}{1-m_4^2} \Big[m_i n_j+m_j n_i - m_4(m_i m_j+n_i n_j)\Big],
\]
and connects source observables to radiated energy, source duration, focal volume, and focal strain [1808.03049].

The same formulation introduces a geometric representation of the source mechanism through the quadratic form \(M_{ij}x_ix_j\). In coordinates aligned with \(\mathbf s\) and \(\mathbf a\), the level set is a rectangular hyperbola, termed the **seismic hyperbola**:
\[
M_{ij}x_ix_j = 2Muv.
\]
Its asymptotes align with the fault normal and slip direction, providing a compact geometric image of focal-region structure [1808.03049].

At smaller scales, source complexity can be resolved through explicit rupture and damage simulation. A hybrid finite-discrete element model of a rough laboratory fault shows that quasi-static shear loading still produces local dynamic seismic activities associated with stress concentration on interlocking asperities. In that study, **7,557 broken CCEs** were clustered into **1,561 seismic events**, with event magnitudes ranging from **\(-11.1\) to \(-4.4\)**. Larger events were mostly shear-mode, smaller events mostly tensile-mode, and damage was spatially heterogeneous, with gouge-forming regions adjacent to minimally damaged patches [2301.04033]. This suggests that even nominally simple fault slip can embed a strongly intermittent seismic source field.

## 2. Monitoring architectures and the ambient seismic environment

Modern seismic monitoring increasingly treats continuous waveform streams, rather than thresholded detections, as the primary observable. In SIGVISA, seismic monitoring is posed as Bayesian inference over a generative model,
\[
p(E,S) = p(S \mid E)\, p(E),
\]
where events \(E\) and station signals \(S\) are jointly modeled. The event prior uses a homogeneous Poisson process, while each station waveform is represented as a sum of phase-wise envelopes, wavelet-based modulation, and autoregressive background noise. Gaussian processes over arrival-time, amplitude, rise-time, decay parameters, and wavelet coefficients allow the model to exploit waveform similarity from historical seismicity while degrading gracefully to parametric envelopes in regions with no nearby training data. On western United States data, this system recovered **three times as many events as previous work** and reduced mean location errors by **a factor of four** [1703.00561].

Observation systems have also broadened beyond standard seismometer networks. A subsea telecom cable from Iceland to Ireland was converted into a distributed seismic sensor using per-span laser interferometry with **17 monitored spans** over about **1770 km**. The system detected clear \(P\)-waves, \(S\)-waves, and surface waves from multiple worldwide earthquakes, grouped arrivals into **13 phase families** within the same minute, and reported about **20 earthquakes detected on several spans**, with **15 detected across all spans**. The Japan example showed the expected sequence of arrivals, with the \(P\)-wave after about **12.3 minutes**, the \(S\)-wave at **22.5 minutes total**, and the surface wave after about **43.5 minutes**, consistent with global travel-time physics [2409.19827].

Seismic background and site characterization remain equally important. At the Sanford Underground Laboratory, an array of broadband instruments established the **4100-ft level** as a **world-class low seismic-noise environment**. Above about **1 Hz**, seismic noise at 4100 ft is about a **factor of 10 weaker in amplitude than at 300 ft**. The study linked the primary and secondary microseismic peaks in the **50 mHz–0.3 Hz** band to ocean-wave activity and showed that low-frequency wind correlation present at **300 ft** disappeared at **2000 ft** and **4100 ft**, indicating that surface wind-driven seismicity is not well described by simple homogeneous-Rayleigh-wave propagation [1006.0678].

Precision infrastructure can itself function as a seismic observer. In the European X-ray Free-Electron Laser, controller I/O data from phase-locked loops in link stabilization units were analyzed as a proxy for ground motion. The study separated earthquakes, ocean-generated microseism, and civilization noise using spectral fingerprints and external data. Distant earthquakes were detectable mainly below roughly **5 Hz**, ocean-generated microseism occupied **\(0.1\text{–}0.3\,\text{Hz}\)**, and civilization noise occupied **\(0.9\text{–}3.5\,\text{Hz}\)**. Even earthquakes about **5000 km away** produced measurable fluctuations, with controller-output jitter spans up to **\(45\,\text{fs}\)** and controller-input spans exceeding **\(4\,\text{as}\)**. The phase-locked loops eliminated **more than 99%** of interference, and the conclusion describes attenuation of **99.99%** of disturbances caused by seismic activity [2502.02453].

## 3. Computational seismic analysis, machine learning, and multimodal models

Seismic data processing has become a major application area for modern machine learning. For missing-trace reconstruction, DSPRecon formulates recovery of complete data \(X\) from masked observations \(M\) through an untrained U-net,
\[
X = f_{\theta}(Z), \qquad \text{s.t. } P_{\Omega}(X)=P_{\Omega}(M),
\]
with parameters optimized by
\[
\theta^{*} = \arg\min_{\theta}\left\|P_{\Omega}\big(f_{\theta}(Z)\big) - P_{\Omega}(M)\right\|_F^2.
\]
The method uses one undersampled seismic record, a fixed random noise tensor, and the architecture itself as a prior. Its reported signal-to-noise ratios were **32.68 dB** versus **19.11 dB** for SSA on synthetic single-shot pre-stack data with **50% randomly missing traces**, **33.09 dB** versus **24.27 dB** for post-stack data, and **35.91 dB** versus **15.32 dB** for regularly missing field traces compared with de-aliased Cadzow [1911.08784].

For seismic event detection, Seismic-Net casts continuous monitoring as sliding-window binary classification on **18,000-timestamp** windows with a **6,000-timestamp offset** at test time. The architecture is a deep densely connected 1D CNN with growth rate **\(k=12\)** and about **800K parameters**. Trained on Chimayó geyser data as a natural analog for \(\mathrm{CO}_2\) leakage, it achieved **precision 0.889** and **recall 0.923**, outperforming a Kernel SVM, a VGG-based CNN, and a ResNet-based CNN [1802.02241].

Synthetic waveform generation addresses label scarcity. SeismoGen uses a conditional GAN for **40 s, 3-channel** waveform windows at **40 Hz**, with three separate generator pipelines for BHE, BHN, and BHZ. The best synthetic-only classifier \(C_{S0}\) achieved **97.11%** accuracy on station V34A and **96.96%** on V35A, compared with **98.56%** and **98.48%** for a real-trained classifier. When synthetic augmentation was added to limited real training sets, classification accuracy improved in **23 out of 30 cases** on V35A and **27 out of 30 cases** on V34A, with largest gains **over 14%** and **over 17%**, respectively [1911.03966].

Representation learning has expanded from task-specific models to foundation models. SeisLM pretrains on a union of **eight seismic datasets**—ETHZ, INSTANCE, Iquique, STEAD, GEOFON, MLAAPDE, PNW, and OBST2024—using a masked contrastive objective with **\(K=100\)** negatives and **\(\kappa=0.1\)**. Two variants are reported: **SeisLM-base** with about **11.4M parameters** and **SeisLM-large** with about **90.7M parameters**. Fine-tuning improves event detection, phase identification, onset time regression, and foreshock–aftershock classification, with the largest gains in low-label regimes [2410.15765].

Seismic analysis is also becoming explicitly multimodal. MultiSeismo assembles **over 16K seismic events** spanning **2010 to 2023**, integrating waveform recordings, intensity maps, population exposure visualizations, and text in a standardized JSON format. The image archive contains **\(N = 91{,}313\)** images, and the instruction set MISCE/MICSE includes **17 instruction templates** and **more than 500K total instructions**. SeisModal, built by augmenting Unified-IO 2 with a pretrained Chronos-T5 time-series encoder, outperformed Phi 4 and Unified-IO on text, image, and especially time-series reasoning; in the reported instruction evaluation, its time-series scores were **BLEU 0.675**, **ROUGE 0.664**, and **BERT 0.824** [2605.26320].

## 4. Imaging, inversion, and interpretation of the subsurface

Seismic inversion addresses the recovery of quantitative subsurface properties from band-limited wavefields. In the post-stack setting considered by IntraSeismic, the forward model is
\[
d(t)= \frac{1}{2} w(t) * \frac{d}{d t} \log (I_p(t)),
\]
or
\[
\mathbf{d} = \mathbf{W D m},
\]
leading to the inverse problem
\[
\mathbf{m^*} = \argmin_{\mathbf{m}} \frac{1}{2}\left\|\mathbf{G m} - \mathbf{d}\right\|_2^2 + \mathcal{R}(\mathbf{m}),
\qquad \mathbf{G}=\mathbf{WD}.
\]
The method represents the impedance field implicitly as \(\boldsymbol{m}_{\Theta}=  \mathbf{F}_{\Theta}(\boldsymbol{x}) + \boldsymbol{m}_{\text{back}}\), using multiresolution hash encoding and a small MLP. On Marmousi, reported SNRs were **27.34 dB**, **25.00 dB**, and **22.93 dB** for noise levels \(\sigma=0\), \(0.1\), and \(0.2\); on a **\(180\times180\times180\)** SEAM Phase I subvolume, the method achieved **40.84 dB**, **36.92 dB**, and **34.5 dB**, converging in under **400 iterations** in each case. The coordinate-based parameterization also reduced storage from about **5.83 million voxel values** to about **1.16 million parameters**, giving an approximate compression ratio of **5** [2312.10568].

ContextSeisNet extends seismic processing with in-context learning. Instead of a fixed mapping \(f_{\theta}(X)=Y\), it conditions predictions on a support set of neighboring gather-label pairs,
\[
f^{\mathrm{ICL}}_{\theta}(X \mid V) = Y.
\]
Applied to demultiple processing, the method uses support examples from spatially related CDPs on the same line and performs inference in a single forward pass without retraining. The reported field-data result is comparable performance with **90% less training data** than the U-Net baseline, together with improved lateral consistency, improved near-offset performance, and more complete multiple removal [2512.11575].

Semantic interpretation has likewise moved toward dense annotation. SpiNet defines **12 commonly observed seismic patterns** grouped into **7 horizons**, **2 stratigraphic sequences**, and **3 structures**, and learns pixel-wise annotations with a deconvolutional encoder–decoder network. Training on **51×51** patches, augmented to **76,950 training images**, the model reached about **0.8** training accuracy and below **0.2** loss, and achieved **78%** pixel-wise annotation accuracy on held-out inline **#390**. Applied to the F3 cube of **651 inlines, 951 crosslines, and 463 samples per trace**, SpiNet annotated the full volume in about **10 minutes on one NVIDIA P4000 GPU** [1810.08517].

## 5. Structured media, forecasting, and seismic engineering

Seismic research also includes deliberate modification and anticipation of wave propagation. Seismic metamaterials translate ideas from electromagnetic metamaterials to **surface seismic wave control in sedimentary soils structured at the meter scale**. Large-scale in-situ experiments near Grenoble and Lyon demonstrated that borehole arrays can redistribute energy, reduce wave amplitude along some directions, change polarization, and act as a seismic filter or lens. The Lyon experiment used a grid of **23 boreholes**, about **2 m diameter** and **5 m depth**, with spacing about **5–7 m**, while a **17 ton** mass drop generated the source. A later reinterpretation emphasized **energy corridors**: some paths between boreholes acted as guiding channels rather than simple barriers [1912.12916].

The corresponding transfer function was defined as
\[
\mid T_{A/O} (\omega)\mid=\mid {\cal A}(\omega)/{\cal O}(\omega)\mid,
\]
and the homogenized effective-medium description used Willis-type constitutive relations,
\[
\Sigma=\C_{\rm eff}:\nabla{\bf A}_0 + {\bf S}^1_{\rm eff} \frac{\partial}{\partial t}{\bf A}_0,
\qquad
\Pi={\bf S}^2_{\rm eff}:\nabla{\bf A}_0 + \rho_{\rm eff} \frac{\partial}{\partial t}{\bf A}_0.
\]
In that framework, ambient seismic noise can be interpreted as a time modulation that makes structured soils behave like moving media. The reported application domain includes shielding, lensing, cloaking, Rayleigh-wave control, energy harvesting, and analog computation using ambient seismic noise [1912.12916].

Forecasting is a complementary engineering direction. SeismoGPT is a transformer-based autoregressive model for forecasting three-component seismic waveforms for future gravitational-wave detectors. It uses token length **16 samples**, context window **64 tokens**, **6 layers**, and **8 attention heads per layer**, and is trained with mean squared error on synthetic waveforms generated with Instaseis from the **ak135f\_2s** database. The paper reports good performance in the immediate prediction window, gradual degradation at longer horizons due to autoregressive error accumulation, and more robust late-forecast behavior for the array-based model than for the single-station version. The intended applications are Newtonian-noise mitigation, active seismic isolation, and real-time observatory control [2509.21446].

## 6. Planetary and icy-ocean-world seismology

Seismic methods have expanded from Earth to the Solar System. A broad review argues that seismic experiments constrain crust, mantle, and core structure, tectonic activity, thermal state, ocean depth, ice-shell thickness, and habitability-relevant layering. InSight on Mars demonstrated that even a **single seismometer** can constrain crustal thickness, mantle properties, and core radius while identifying active source regions such as Cerberus Fossae. The same review treats the Moon, Europa, Titan, Enceladus, Io, Mercury, Venus, Ceres, and the giant planets as potential seismic targets, each with distinct mission and environmental constraints [2206.01785].

Icy ocean worlds require a revised seismic lexicon because an ice shell overlies a global ocean. A dedicated framework introduces \(P\) and \(S\) for compressional and shear waves in ice or mantle, \(F\) for ocean legs, and \(K\) for compressional waves in the core, together with boundary markers such as \(e\) for the bottom of the ice shell and \(o\) for the bottom of the ocean. It also emphasizes guided phases such as flexural waves, longitudinal waves, Love waves, and the Crary phase. For a floating ice shell, the longitudinal-wave speed is
\[
v_L = 2v_S \sqrt{1 - \left(\frac{v_S}{v_P}\right)^2},
\]
and the Crary-phase frequency is
\[
f_{\mathrm{Cr}} = \frac{(n+1)\,v_S}{2d \sqrt{1 - \left(\frac{v_S}{v_P}\right)^2}}.
\]
These relations make guided-wave observations direct constraints on ice thickness and elastic structure [1705.03500].

Europa-specific seismicity models use a Gutenberg–Richter framework with cumulative seismic moment release between
\[
10^{16}\text{ and }10^{18}\,\mathrm{Nm/yr}.
\]
Across the explored models, ambient noise levels vary by about **80 dB**. Most mean noise estimates are below the self-noise floor of a typical high-frequency geophone, but larger events are expected to stand about **\(\sim 30\)–50 dB** above the mean background. The study also shows that autocorrelation of simulated noise can recover reflection-like arrivals, including an ocean-floor reflection at about **175 s** in one example [1705.03424].

Titan introduces a different microseismic regime because its methane–ethane seas can generate ambient seismic noise analogous to Earth’s ocean microseisms. Using wind-wave and loading models, one study concludes that storms of more than **2 m/s** wind speed would create a signal that is globally observable with a high-quality broadband sensor and observable to a **thousand kilometer** distance with an InSight-class space-ready seismometer. The modeled Titan microseisms lie mainly between about **0.05 and 0.2 Hz**, and are expected to be episodic rather than constant [1905.11251].

Taken together, these results show that seismic research is no longer confined to terrestrial earthquake seismology. It includes analytic source inversion, ambient-noise characterization, distributed sensing on cables and precision timing systems, data-driven reconstruction and interpretation, structured control of wave propagation, and planetary interior probing in environments where oceans, ice shells, and nonstandard guided phases fundamentally reshape the wavefield.

Source: https://www.emergentmind.com/topics/seismic