Abstract:Summary Full waveform inversion has emerged as the state-of-the art high resolution seismic imaging technique, both in seismology for global and regional scale imaging and in the industry for exploration purposes. While gaining in popularity, full waveform inversion, at an operational level, remains a heavy computational process involving the repeated solution of large-scale 3D wave propagation problems. For this reason it is a common practice to focus the interpretation of the results on the final estimated model. This is forgetting full waveform inversion is an ill-posed inverse problem in a high dimensional space for which the solution is intrinsically non-unique. This is the reason why being able to qualify and quantify the uncertainty attached to a model estimated by full waveform inversion is key. To this end, we propose to extend at an operational level the concepts introduced in a previous study related to the coupling between ensemble Kalman filters and full waveform inversion. These concepts had been developed for 2D frequency-domain full waveform inversion. We extend it here to the case of 3D time-domain full waveform inversion, relying on a source subsampling strategy to assimilate progressively the data within the Kalman filter. We apply our strategy to an ocean bottom cable field dataset from the North Sea to illustrate its feasibility. We explore the convergence of the filter in terms of number of elements, and extract variance and covariance information showing which part of the model are well constrained and which are not. Analyzing the variance helps to gain insight on how well the final estimated model is constrained by the whole full waveform inversion workflow. The variance maps appears as the superposition of a smooth trend related to the geometrical spreading and a high resolution trend related to reflectors. Mapping lines of the covariance (or correlation matrix) to the model space helps to gain insight on the local resolution. Through a wave propagation analysis, we are also able to relate variance peaks in the model space to variance peaks in the data space. Compared to other posterior-covariance approximation scheme, our combination between Ensemble Kalman filter and full waveform inversion is intrinsically scalable, making it a good candidate for exploiting the recent exascale high performance computing machines.

TSUFLIND-EnKF inversion model applied to tsunami deposits for estimation of transient flow depth and speed with quantified uncertainties

Inversion of Tsunamis Characteristics from Sediment Deposits Based on Ensemble Kalman Filtering

A model for tsunami flow inversion from deposits (TSUFLIND)

A Numerical Study of Tsunami Wave Impact and Run-Up on Coastal Cliffs Using a CIP-based Model

Tracking of time-evolving sound speed profiles in shallow water using an ensemble Kalman-particle filter.

Predicting soil settlement with quantified uncertainties by using ensemble Kalman filtering

Quantifying uncertainties in fault slip distribution during the Tōhoku tsunami using polynomial chaos

Simulation of tsunami induced by a submarine landslide in a glaciomarine margin: the case of Storfjorden LS-1 (southwestern Svalbard Islands)

Numerical simulation of submarine landslide tsunamis based on the smoothed particle hydrodynamics model

GeoClaw-STRICHE: A coupled model for Sediment TRansport In Coastal Hazard Events

Quantifying Uncertainties in the 2004 Sumatra-Andaman Earthquake Source Parameters by Stochastic Inversion

Validation of Linearity Assumptions for Using Tsunami Waveforms in Joint Inversion of Kinematic Rupture Models: Application to the 2010 Mentawai M-W 7.8 Tsunami Earthquake

Resolving Shallow Moment Releases of Megathrust Earthquakes by Empirical Green’s Function‐Based Inversions: A New Approach for Tsunami Warning

Landslide-induced tsunami simulation based on progressive landslide-shallow water equation coupling model: 1946 Aleutian tsunami case

Local uncertainty quantification for 3D time-domain full waveform inversion with ensemble Kalman filters: application to a North sea OBC dataset

Understanding the tsunami with a simple model

Dynamics of tsunami wave propagation in uncertain environment

Numerical Implementation of Wave Friction Factor into the 1D Tsunami Shallow Water Equation Model

X-ray tomography applied to tsunami deposits: Optimized image processing and quantitative analysis of particle size, particle shape, and sedimentary fabric in 3D

Effects of Dispersion in Tsunami Green'S Functions and Implications for Joint Inversion with Seismic and Geodetic Data: A Case Study of the 2010 Mentawai M-W 7.8 Earthquake

Adjoint Synthesis for Trans‐Oceanic Tsunami Waveforms and Simultaneous Inversion of Fault Geometry and Slip Distribution