A proper methodology aimed at surface wave tomography

When applying a methodology for obtaining the 3D shear-wave velocity structure of a medium from surface wave dispersion data, the problem must be considered with caution since one inverts path-averaged velocities and the use of any inversion method entails some drawbacks such as lack of uniqueness, unwarranted stability and constraints affecting the data. In order to avoid the application of consecutive inversions and to overcome these drawbacks, we propose alternative mapping methods, for example spatial prediction methods, or else the use of an algorithm that, from a mathematical viewpoint, can be understood through the application of the orthogonal projection theorem onto convex sets (POCS). Among the first ones, we try inverse weighted distance interpolation. The POCS algorithm we have used discretises a second order differential equation for the velocity field with boundary conditions. All these imaging techniques aimed at volumetric modelling and the visualisation of data are discussed, and finally we show some results based on ray path velocities obtained previously by inversion of phase and group velocities of Rayleigh waves propagating across the Iberian peninsula.