Computation of Large Anisotropic Seismic Heterogeneities (CLASH)

• Published in: Geophysical Journal International
• Main Authors: Beucler, Éric, Montagner, Jean-Paul
• Other Authors: Laboratoire de Planétologie et Géodynamique UMR 6112 (LPG), Université d'Angers (UA)-Université de Nantes - UFR des Sciences et des Techniques (UN UFR ST), Université de Nantes (UN)-Université de Nantes (UN)-Institut national des sciences de l'Univers (INSU - CNRS)-Centre National de la Recherche Scientifique (CNRS), Institut de Physique du Globe de Paris (IPGP), Institut national des sciences de l'Univers (INSU - CNRS)-Université Paris Diderot - Paris 7 (UPD7)-Université de La Réunion (UR)-Institut de Physique du Globe de Paris (IPG Paris)-Centre National de la Recherche Scientifique (CNRS)
• Format: Article in Journal/Newspaper
• Language: English
• Published: HAL CCSD 2005
• Subjects: azimuthal anisotropy; global phase velocity maps; Rayleigh surface waves; seismic tomography; [SDU]Sciences of the Universe [physics]; Antarc*; Antarctica
• Online Access: https://insu.hal.science/insu-03603142; https://insu.hal.science/insu-03603142/document; https://insu.hal.science/insu-03603142/file/165-2-447.pdf; https://doi.org/10.1111/j.1365-246X.2005.02813.x

International audience A general tomographic technique is designed in order (i) to operate in anisotropic media; (ii) to account for the uneven seismic sampling and (iii) to handle massive data sets in a reasonable computing time. One modus operandi to compute a 3-D body wave velocity model relies on surface wave phase velocity measurements. An intermediate step, shared by other approaches, consists in translating, for each period of a given mode branch, the phase velocities integrated along ray paths into local velocity perturbations. To this end, we develop a method, which accounts for the azimuthal anisotropy in its comprehensive form. The weakly non-linear forward problem allows to use a conjugate gradient optimization. The Earth's surface is regularly discretized and the partial derivatives are assigned to the individual grid points. Possible lack of lateral resolution, due to the inescapable uneven ray path coverage, is taken into account through the a priori covariances on parameters with laterally variable correlation lengths. This method allows to efficiently separate the 2ψ and the 4ψ anisotropic effects from the isotropic perturbations. Fundamental mode and overtone phase velocity maps, derived with real Rayleigh wave data sets, are presented and compared with previous maps. The isotropic models concur well with the results of Trampert & Woodhouse. Large 4ψ heterogeneities are located in the tectonically active regions and over the continental lithospheres such as North America, Antarctica or Australia. At various periods, a significant 4ψ signature is correlated with the Hawaii hotspot track. Finally, concurring with the conclusions of Trampert & Woodhouse, our phase velocity maps show that Rayleigh wave data sets do need both 2ψ and 4ψ anisotropic terms.