Quasi-Newton inversion of seismic first arrivals using source finite bandwidth assumption: Application to subsurface characterization of landslides

Characterizing the internal structure of landslides is of first importance to assess the hazard. Many geophysical techniques have been used in the recent years to image these structures, and among them is seismic tomography. The objective of this work is to present a high resolution seismic inversion algorithm of first arrival times that minimizes the use of subjective regularization operators. A Quasi-Newton P-wave tomography inversion algorithm has been developed. It is based on a finite frequency assumption for highly heterogeneous media which considers an objective inversion regularization (based on the wave propagation principle) and uses the entire source frequency spectrum to improve the tomography resolution.The Fresnel wavepaths calculated for different source frequencies are used to retropropagate the traveltime residuals, assuming that in highly heterogeneous media, the first arrivals are only affected by velocity anomalies present in the first Fresnel zone. The performance of the algorithm is first evaluated on a synthetic dataset, and further applied on a real dataset acquired at the Super-Sauze landslide which is characterized by a complex bedrock geometry, a layering of different materials and important changes in soil porosity (e.g. surface fissures). The seismic P-wave velocity and the wave attenuation are calculated, and the two tomographies are compared to previous studies on the site.

► We implement a new Q-N inversion of first arrivals using Fresnel wavepaths.
► We consider the source finite bandwidth to improve the tomography's resolution.
► We invert the seismic P-wave attenuation.
► This algorithm permits to recover accurately the paleotopography.
► We created a geological model of the profile containing geotechnical structures.