Full-waveform inversion (FWI) is a critical technique for deriving high-resolution velocity models in geophysical studies. However, the lack of low-frequency components in seismic data can result in significant cycle-skipping issues, leading to inaccurate inversion outcomes. Fortunately, the envelope of seismic data contains abundant low-frequency information that can be utilized to improve the inversion of large-scale background velocity models. Within the FWI framework, current envelope inversion (EI) methods can be classified into two primary categories: (1) conventional envelope inversion, which employs the waveform Fréchet derivative and is optimally suited for weak scattering scenarios, and (2) direct envelope inversion, which utilizes the envelope Fréchet derivative and demonstrates effectiveness in inverting models characterized by strong contrasts. This review paper addresses the following key aspects: (1) an overview of the principles underlying conventional envelope inversion and direct envelope inversion methodologies, alongside a detailed comparative analysis; (2) an evaluation of the limitations and effectiveness of both conventional and direct envelope inversion techniques when applied to weak and strong scattering media, supported by comprehensive numerical test cases; and (3) a thorough review of the historical development of envelope inversion, exploring its critical challenges, current solutions, and future perspectives in the field.
Seismic data often face challenges in their utilization due to noise contamination, incomplete acquisition, and limited low-frequency information, which hinder accurate subsurface imaging and interpretation. Traditional processing methods rely heavily on task-specific designs to address these challenges and fail to account for the variability of data. To address these limitations, we present a generative seismic foundation model (GSFM), a unified framework based on generative diffusion models (GDMs), designed to tackle multi-task seismic processing challenges, including denoising, backscattered noise attenuation, interpolation, and low-frequency extrapolation. GSFM leverages a pre-training stage on synthetic data to capture the features of clean, complete, and broadband seismic data distributions and applies an iterative fine-tuning strategy to adapt the model to field data. By adopting a target-oriented diffusion process prediction, GSFM improves computational efficiency without compromising accuracy. Synthetic data tests demonstrate that GSFM surpasses benchmarks with equivalent architectures in all tasks and achieves performance comparable to traditional pre-training strategies, even after their fine-tuning. Also, field data tests suggest that our iterative fine-tuning approach addresses the generalization limitations of conventional pre-training and fine-tuning paradigms, delivering significantly enhanced performance across diverse tasks. Furthermore, GSFM’s inherent probabilistic nature enables effective uncertainty quantification, offering valuable insights into the reliability of processing results.
Recent studies of high-resolution seismic tomography of source zones of large crustal earthquakes, megathrust earthquakes, and intraslab earthquakes are reviewed, which shed new light on seismogenic structures and fluids in subduction zones. Large crustal earthquakes generally occurred in high-velocity (high-V) zones in the brittle upper crust, whereas low-velocity and high Poisson’s ratio anomalies exist in the lower crust and upper (or uppermost) mantle, which may reflect fluids released from dehydration of the subducting slab. The fluids may trigger large crustal earthquakes. The interplate megathrust zone exhibits prominent structural heterogeneities. Large megathrust earthquakes generally occurred in high-V areas, reflecting strongly coupled patches (or asperities) in the megathrust zone due to the subduction of seamounts or topographic plateaus in the incoming oceanic plate. The megathrust seismogenesis may be affected or controlled by structural anomalies in both the upper and lower plates, as well as hot upwelling flows in the subslab mantle. Lower-velocity anomalies are revealed in source zones of large intraslab earthquakes, which are attributed to the process of dehydration embrittlement resulting from dehydration of hydrous minerals in the slab, which may trigger the mainshock and aftershock sequences by enhancing pore pressures along preexisting faults and fractures in the slab. All these results indicate that fluids play an important role in the generation of most earthquakes in subduction zones.
Three-dimensional (3-D) magnetotelluric (MT) inversion outcomes are influenced by various user-defined configurations, yet comprehensive analyses and rigorous testing of them are rarely conducted, particularly for 3-D large-scale surveys. Additionally, different inversion frameworks will inevitably yield varying “preferred” models, highlighting the need to investigate the reasons behind these discrepancies, which can uncover model uncertainties linked to algorithm-specific choices. This study focuses on analyzing and comparing factors from the above two categories that contribute to uncertainties in inversions and interpretations. First, to explore the user-defined configurations that lead to differing results, we performed multiple tests using the USArray MT data. Specifically, we assessed the impact of regularization parameters, prior information, projection methods, regularization types, data types and frequencies, initial models, mesh settings, and distortion correction on inversion outcomes by the two most widely used inversion frameworks: ModEM and FEMTIC. After thoroughly analyzing these factors, we present a new electrical model for the Northwestern U.S. using FEMTIC code. The primary conductivity variation of this model aligns well with a previous one obtained by utilizing ModEM. However, large-scale differences between the two “preferred” solutions are still observed. We conclude that the discrepancies in geologically active high-conductivity zones are mainly due to different selection method for regularization parameter and mesh settings, and difference in average resistivities of the lower upper mantle primarily stem from differing degrees of dependence on the initial model, a distinction that fundamentally stems from differences in regularization assumptions and optimization algorithms.
Tsunami is the oceanic gravity waves produced by mass displacements in large shallow offshore thrust earthquakes. The January 1, 2024 (M_textrm{W})7.5 Noto Peninsula earthquake excited a tsunami to spread across the East Sea (Sea of Japan), arriving at the east coast of the Korean Peninsula. The influence of oceanic gravity waves on the coastal medium is investigated. The mass loading by tsunami induces ground tilting, producing transient long-period ground motions to be polarized in coastline-perpendicular directions. The tsunami-induced ground motions are well recorded in inland seismometers nearby the coast. The wavetrain durations and spectral contents of the tsunami-induced seismic signals in seismometers share with those of the tsunami waves in tide gauges, suggesting the same source of energy. The amplitudes of tsunami-induced ground motions are proportional to the tsunami heights, being modulated by the distance from the coast and medium properties. The discriminative tsunami-induced ground motions produce dynamic stress changes that are effective at shallow depths, reaching 0.81 kPa on the coast. A large runup height may induce dynamic stress changes effective to depths.

