A Monte Carlo multimodal inversion of surface waves

Abstract

S U M M A R YThe analysis of surface wave propagation is often used to estimate the S-wave velocity profile at a site. In this paper, we propose a stochastic approach for the inversion of surface waves, which allows apparent dispersion curves to be inverted. The inversion method is based on the integrated use of two-misfit functions. A misfit function based on the determinant of the Haskell-Thomson matrix and a classical Euclidean distance between the dispersion curves. The former allows all the modes of the dispersion curve to be taken into account with a very limited computational cost because it avoids the explicit calculation of the dispersion curve for each tentative model. It is used in a Monte Carlo inversion with a large population of profiles. In a subsequent step, the selection of representative models is obtained by applying a Fisher test based on the Euclidean distance between the experimental and the synthetic dispersion curves to the best models of the Monte Carlo inversion. This procedure allows the set of the selected models to be identified on the basis of the data quality. It also mitigates the influence of local minima that can affect the Monte Carlo results. The effectiveness of the procedure is shown for synthetic and real experimental data sets, where the advantages of the two-stage procedure are highlighted. In particular, the determinant misfit allows the computation of large populations in stochastic algorithms with a limited computational cost.

Time-lapse monitoring of climate effects on earthworks using surface waves

“…Three dispersion curves, each considered representative of parts of the seismic line (x < 35 m, 35 < x < 75 m, and x > 75 m), were selected and inverted. The MC inversion algorithm developed by Maraschini and Foti (2010) was used. This approach enables a computationally efficient search of the parameter space to be carried out.…”

Crustal and uppermost mantle S-wave velocity below the East European Craton in northern Poland from the inversion of ambient-noise records

“…It involves an iterative search technique to seek for the global minimum of the multidimensional error function (EF) in the parameter space without being influenced by the chosen starting model. The best model was identified by iterative inversions, comparing each SWD curve with a group of theoretical dispersion curves belonging to the same population of that SWD curve within a preset confidence limit (acceptable models, Maraschini and Foti 2010). The theoretical curves were generated from a set of representative Fig.…”

The geophysical characteristic of the lower lithosphere and asthenosphere in the marginal zone of the East European Craton

“…To find the solution of the inversion process, the MCNA involves a nonlinear iterative search technique to seek the global minimum of the error function without being influenced by the chosen starting model. Therefore, the final S-wave velocity models were evaluated by iterative inversions, comparing the dispersion curves with a group of theoretical ones belonging to the same population within a preset confidence limit (Maraschini and Foti 2010). The best model (Fig.…”