REVIEW 4 major objections 5 minor 1 cited by
On extracting coherent seismic wavefield using variational symmetric autoencoders
T0 review · 4 major / 5 minor · reviewed 2026-08-12 · deepseek-v4-flash
Pith's one-line read A variational symmetric autoencoder disentangles shared earthquake source information from station-specific scattering, and reconstructs clean virtual seismograms without deconvolution.
desk verdict A genuinely useful new formulation for extracting source time functions without deconvolution, but the paper undersells its own contribution by leaning on qualitative validation and skipping baselines. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The central object is the Symmetric Autoencoder (SymAE), a variational autoencoder whose latent space is split into a coherent code $s$ shared by a whole group of waveforms and nuisance codes $p_k$ for each waveform. The load-bearing mechanism is the conditional independence assumption $P(x_j^k, x_j^l | s, p_k, p_l) = P(x_j^k | s, p_k) P(x_j^l | s, p_l)$ (eq. 9), which lets the encoder factorize the approximate posterior and accumulate source information across receivers as a product of Gaussians; the resulting precision-weighted mean and variance (eqs. 17-18) are a probabilistic generalization of nonlinear stacking. A time-shift transformer module, built on spatial-transformer-style localization networks, aligns waveforms before encoding so the model can handle timing variations without explicit cross-correlation.
What would settle it
Take the synthetic setup of Fig. 1 and add a common scattered coda, identical across all receivers in a group, on top of the per-receiver path effects; if the recovered source time function changes or contains that common coda, or if the eq. (19) ranking starts preferring contaminated traces, the conditional-independence assumption is broken and the extraction is not purely source information.
Extended reading notes
Core claim
The central claim is that SymAE learns a disentangled latent representation in which one component, $s$, accumulates all information about the earthquake source shared by the seismograms in a group, while the remaining components $p_k$ carry waveform-specific nuisance. The accumulation is done by multiplying the per-seismogram approximate posteriors, which for Gaussians reduces to precision-weighted averages (eqs. 17-18). The paper defines an informativeness measure based on the KL divergence between $Q(s|X_j)$ and $Q(s|x_j^k)$ (eq. 19), shows that low-divergence seismograms carry reliable source content, and solves an optimization (eq. 21) for the nuisance code that minimizes this divergence while enforcing sparsity and normalization, thereby generating a virtual seismogram for each earthquake. On synthetic data the extracted source signatures match the ground-truth sources, and on real deep earthquakes the SymAE-derived source time functions are clearer and capture complex multi-event ruptures that traditional deconvolution obscures, all without an explicit deconvolution step.
Load-bearing premise
The method assumes that once the shared earthquake source is known, the recordings at different stations carry no additional shared information, so any leftover common path, noise, or site effect would be folded into the learned source and silently treated as source signal.
Editorial extensions
If this is right
- Analysis of complex earthquakes with multiple rupture episodes becomes possible without deconvolution, because the source code $s$ is learned from the waveform population rather than inverted from a known path response.
- The KL informativeness score (eq. 19) gives a principled way to rank, weight, or discard seismograms according to source content, potentially improving receiver-function and ambient-noise stacking workflows.
- Virtual seismograms formed by swapping path codes between similar earthquakes (eq. 20) provide a data-driven redatuming tool, so a clear recording from one event can help visualize a smaller or more obscured event with a comparable source.
- Source-similarity structure recovered from the information-loss heatmap (Fig. 6) aligns with known rupture complexity, suggesting the model can serve as an unsupervised grouping tool for earthquakes.
Reading between the lines
- Editorial inference: If eq. (9) is violated, for example when receivers share near-source structure or common mantle heterogeneity, that correlated path energy will leak into $s$, and because the same encoder scores informativeness, the leakage may not be visible in the KL metric; an external comparison against back-projection or moment-tensor inversions would expose it.
- Editorial inference: The same accumulation-by-conjunction recipe applies directly to receiver-based groups: treating a station's waveforms across many earthquakes as the group would let SymAE extract a coherent receiver function or site response without stacking individual teleseisms.
- Editorial inference: A natural stress test is to add a synthetic common scattered coda to all receivers in a group; if the recovered source time function changes or the informativeness ranking shifts toward contaminated traces, the method's core independence assumption is falsified.
- Editorial inference: The sparsity weight in eq. (21) controls a trade-off between temporal resolution and noise suppression; tuning it per target frequency band may let seismologists tailor virtual source functions to specific rupture processes.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper presents a variational symmetric autoencoder (SymAE) that disentangles a set of seismograms into a shared source code and per-seismogram path codes. It introduces a KL-divergence-based informativeness metric (Eq. 19) and an optimization procedure (Eq. 21) to generate virtual seismograms enriched in coherent source information. The method is applied to synthetic data with known source time functions and to teleseismic P-wave records from deep-focus earthquakes, with comparisons to traditional deconvolution. The central claim is that SymAE extracts coherent source information without deconvolution, enabling analysis of complex earthquakes with multiple rupture episodes.
Significance. If the claims are substantiated, the paper offers a novel nonlinear stacking/representation framework for seismic source extraction, building on prior SymAE work and extending it with a practical informativeness criterion and a latent-space optimization for virtual seismograms. The variational formulation in Sections 2-3 is standard and internally consistent, and the synthetic experiment with known sources is a well-posed falsifiable test. The paper also states that all codes and data are open access, which is a strength for reproducibility, though no repository is provided. However, the current experimental validation is insufficient to establish the central claims, particularly because the synthetic results are only visually compared and the real-data comparison is qualitative.
major comments (4)
- [§3, Eq. (9)] The conditional independence assumption P(x_k, x_l | s, p_k, p_l) = P(x_k | s, p_k) P(x_l | s, p_l) is load-bearing for the precision-weighted accumulation in Eqs. (17)-(18). For teleseismic P waves, receiver waveforms share path-related structure beyond the source (e.g., near-source heterogeneity, common mantle attenuation, source-side scattering) that has no dedicated latent variable in the model. When Eq. (9) is violated, shared path energy has nowhere to go except the shared code s, so the learned 'source' code is actually source plus common path. The synthetic experiment uses independent receiver-specific path effects and therefore cannot detect this failure. I recommend adding a synthetic experiment with correlated or partially shared path components across receivers and measuring whether the recovered s matches the true source under those conditions.
- [§3.3-3.5, Eqs. (19)-(21)] Both the informativeness metric H and the optimal virtual seismogram objective in Eq. (21) are evaluated using the same trained encoder Q(s|x) that was fit under Eq. (9). If s contains shared path energy, these self-evaluations will not reveal the failure; H will be small even when the virtual seismogram is not source-dominated. The synthetic validation in Fig. 9 is only visual: no quantitative recovery metric (e.g., correlation or L2 error between the optimized virtual seismogram and the true band-limited source), no error bars, and no random-seed variability. Please add quantitative metrics with multiple seeds and, if possible, an out-of-sample decoder or a separate validation encoder to break the circularity.
- [§5, Fig. 11] The real-data comparison against SCARDEC deconvolution is purely qualitative. The claim that SymAE 'captures complex features due to multiple events' is supported only by visual inspection. There is no external source time function, no quantitative similarity metric, and no uncertainty analysis. The statement that 'the half-durations of the source functions extracted using SymAE closely align with those determined from the raw seismograms' is not quantified. Please provide a quantitative comparison (e.g., cross-correlation, misfit, or at least a table of half-durations with uncertainties) and, if available, compare with independent source models for these events.
- [§3.3-3.4 and Fig. 5] The text contradicts itself on the sign of H. Section 3.3 states that seismograms with lower H values are more informative (Eq. 19), while Section 3.4 states that 'seismograms that exhibit higher H values (refer to eq. 19) are considered informative' and Fig. 5's caption repeats this. This reversal is not cosmetic: it governs the interpretation of Figs. 4-6 and the redatuming results. Please reconcile the definition and apply it consistently throughout.
minor comments (5)
- [§3.5, Eq. (21)] The display of the constraints in Eq. (21) is garbled; please write the zero-mean and unit-variance constraints clearly (e.g., sum_n (y_j[m])_n = 0 and (1/D) sum_n ((y_j[m])_n - mean)^2 = 1).
- [§3.4, Fig. 6] The heatmap description is ambiguous: the text says 'higher values after normalization indicate greater similarity' while the caption says lighter colors indicate lower information loss; these are opposite unless the normalization inverts the KL values. Please clarify which quantity is plotted.
- [§5] The paper states that SymAE is trained on 'roughly 5000 displacement seismograms (all components)' but does not specify the number of receivers, the time window, preprocessing details, network architecture (layer counts, latent dimensions), or the values of hyperparameters α and γ. Please add these details for reproducibility.
- [§8] The data availability section says all data and codes are open access, but no repository or DOI is provided. Please include a link or accession code.
- [References] There are several reference formatting errors, e.g., the Weaver 2001 entry is garbled ('RL Weaver. Lobkis 0 i. Ultrasonic without a sources Thermal fluctuation correlations at MHz frequencies...'). Please correct the reference list.
Circularity Check
Minor self-referential informativeness metric, but the main source-extraction claim is independently grounded by the synthetic test and by comparison with deconvolution.
-
self definitional
[Sec. 3.3, Eq. (19); used again in Sec. 3.5, Eq. (21)]
"After training SymAE, we can identify which specific seismograms contribute most significantly to the accumulated source information. This is accomplished by calculating the Kullback-Leibler (KL) divergence between Q(s | Xj), which represents the accumulated distribution, and Q(s | xk_j ), the distribution for an individual seismogram xk_j : H[xk_j ] = DKL(Q(s | Xj) ∥ Q(s | xk_j )). (19)"
The informativeness metric is defined entirely from the trained encoder's own outputs: Q(s|Xj) is the precision-weighted product of the individual posteriors (Eqs. 16-18), so H[xk_j] measures, by construction, how close one waveform's learned posterior is to the product of all learned posteriors. It is not an independent measure of source content. The same self-defined distance is then reused in Eq. (20) and as the primary objective in Eq. (21) to generate and evaluate virtual seismograms, so the model is its own judge at the metric level. The synthetic experiment (Fig. 9), which compares extracted waveforms with known sources, is an external check and prevents this self-referentiality from becoming load-bearing for the main source-extraction claim.
full rationale
The derivation of the SymAE objective (Eq. 14) follows from standard variational inference, and no step reduces a prediction to a fitted parameter. The conditional independence in Eq. 9 is an explicit modeling assumption rather than a hidden input, and the paper tests this framework in a synthetic setting with known sources and against a standard deconvolution method (Fig. 11). Citations to Bouchacourt et al. (2018) and Bharadwaj et al. (2022) are used for background and for the standard Gaussian-conjugate accumulation result, so they are not load-bearing self-citations. The one notable self-referential aspect is that the informativeness metric H in Eq. 19 and the Eq. 21 objective are evaluated with the trained encoder's own posteriors; this could conceal failures if correlated path energy leaks into the source code s. However, because the synthetic test provides an independent target and the real-data comparison is made against an external deconvolution method, the central claim retains independent content. Overall this is a minor self-referential evaluation step rather than a circular derivation, giving a score of 2.
Assumptions & free parameters
free parameters (3)
- alpha (sparsity weight, eq. 21) =
not reported
- gamma (time-shift regularization, eq. 22) =
not reported
- Latent dimensions and network architecture =
not reported
assumptions (5)
- domain assumption Conditional independence of waveforms given the shared source code (eq. 9)
- domain assumption Prior independence of source and path codes (eq. 8)
- standard math Factored approximate posterior (eq. 13)
- standard math Gaussian likelihood and diagonal-covariance Gaussian posteriors
- domain assumption The learned Q(s|x_k) is a faithful estimator of source evidence
invented entities (2)
-
Latent source code s
-
Per-seismogram latent path code p_k
Cite this review
Pith. "Pith review of On extracting coherent seismic wavefield using variational symmetric autoencoders." pith.science (2026). https://pith.science/paper/GN4UZUMC
@misc{pith2026241115613,
author = {Pith},
title = {Pith review of: On extracting coherent seismic wavefield using variational symmetric autoencoders},
year = {2026},
howpublished = {\url{https://pith.science/paper/GN4UZUMC}},
note = {Machine review of arXiv:2411.15613}
}
read the original abstract
We discuss the variational formulation of the Symmetric Autoencoder (SymAE) and its role in achieving disentanglement within the latent space to extract coherent information from a collection of seismic waveforms. Disentanglement involves separating the latent space into components for coherent information shared by all waveforms and components for waveform-specific nuisance information. SymAE employs a generative model that independently generates waveforms based on coherent and nuisance components, and an inference model that estimates these components from observed wavefield. By assuming the independence of waveforms conditioned on coherent information, the model effectively accumulates this information across multiple waveforms. After training, a metric based on Kullback-Leibler divergence is used to evaluate the informativeness of individual waveforms, enabling latent-space optimization and the generation of synthetic seismograms with enhanced signal-to-noise ratios. To demonstrate the efficacy of our proposed method, we applied it to a data set of teleseismic displacement waveforms of the P wave from deep-focus earthquakes. By training the SymAE model on high-magnitude events, we successfully identified seismograms that contained robust source information. Furthermore, we generated high-resolution virtual seismograms enriched with relevant coherent source information and less influenced by scattering noise, allowing a deeper understanding of the characteristics of the earthquake source. Importantly, our method extracts coherent source information without relying on deconvolution, which is often used in traditional source imaging. This enables the analysis of complex earthquakes with multiple rupture episodes, a capability that is not easily achievable with conventional approaches.
Figures
Figures from the paper (8 more)
Forward citations
Cited by 1 Pith paper
-
Enhanced receiver function imaging of crustal structures using symmetric autoencoders
A symmetric variational autoencoder separates coherent crustal signal from per-earthquake noise in receiver functions, producing cleaner virtual receiver functions than standard averaging.
Reference graph
Works this paper leans on
-
[8]
doi:10.1016/j.jseaes.2015.06.007
ISSN 1367-9120. doi:10.1016/j.jseaes.2015.06.007. RL Weaver. Lobkis 0 i. Ultrasonic without a sources Thermal fluctuation correlations at MHz frequencies Physical Review Letters, 87(13):134301,
-
[1995]
doi:10.1038/374152a0. T. Zheng, Z. Ding, J. Ning, L. Chang, X. Wang, F. Kong, and et al. Crustal azimuthal anisotropy beneath the southeastern tibetan plateau and its geodynamic implications. Journal of Geophysical Research: Solid Earth, 123:9733–9749,
-
[2004]
Auto-encoding variational bayes
Diederik P Kingma and Max Welling. Auto-encoding variational bayes. arXiv preprint arXiv:1312.6114,
-
[2005]
The 2018 fiji mw 8.2 and 7.9 deep earthquakes: One doublet in two slabs
Zhe Jia, Zhichao Shen, Zhongwen Zhan, Chenyu Li, Zhigang Peng, and Michael Gurnis. The 2018 fiji mw 8.2 and 7.9 deep earthquakes: One doublet in two slabs. Earth and Planetary Science Letters, 531:115997,
work page 2018
-
[2015]
Comparison of stacking methods regarding processing and computing of geoscientific depth data
Claus-Peter Rückemann. Comparison of stacking methods regarding processing and computing of geoscientific depth data. In Proceedings of the Fourth International Conference on Advanced Geographic Information Systems, Applications, and Services (GEOProcessing 2012), volume 30, pages 35–40,
work page 2012
-
[2017]
Understanding disentangling in beta-vae
Christopher P Burgess, Irina Higgins, Arka Pal, Loic Matthey, Nick Watters, Guillaume Desjardins, and Alexander Lerchner. Understanding disentangling in beta-vae. arXiv preprint arXiv:1804.03599,
-
[2018]
doi:10.1029/2018JB015995. B. Dalai, P. Kumar, U. Srinu, and M. K. Sen. De-noising receiver function data using the unsupervised deep learning approach. Geophysical Journal International, 229(2):737–749, December
-
[2020]
doi:https://doi.org/10.1016/j.epsl.2019.115997
ISSN 0012-821X. doi:https://doi.org/10.1016/j.epsl.2019.115997. URL https://www.sciencedirect.com/science/ article/pii/S0012821X19306892. Diederik P Kingma and Jimmy Ba. Adam: A method for stochastic optimization. arXiv preprint arXiv:1412.6980,
Show all 12 references
-
[2021]
doi:10.1093/gji/ggab494. W. Bloch, M. G. Bostock, and P. Audet. A cascadia slab model from receiver functions. Geochemistry, Geophysics, Geosystems, 24(10):e2023GC011088,
-
[2022]
Spatial transformer networks
Max Jaderberg, Karen Simonyan, Andrew Zisserman, and Koray Kavukcuoglu. Spatial transformer networks. In Advances in neural information processing systems, pages 2017–2025,
2017
-
[2023]
Jiafu Hu, Haiyan Yang, Guangquan Li, and Hengchu Peng
doi:10.1029/2023GC011088. Jiafu Hu, Haiyan Yang, Guangquan Li, and Hengchu Peng. A review on the analysis of the crustal and upper mantle structure using receiver functions. Journal of Asian Earth Sciences, 111:589–603,
-
[2024]
Learning physics for unveiling hidden earthquake ground motions via conditional generative modeling
Pu Ren, Rie Nakata, Maxime Lacour, Ilan Naiman, Nori Nakata, Jialin Song, Zhengfa Bi, Osman Asif Malik, Dmitriy Morozov, Omri Azencot, et al. Learning physics for unveiling hidden earthquake ground motions via conditional generative modeling. arXiv preprint arXiv:2407.15089,
Reviewed August 12, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.