Analyse des intercorrélations du champ d'onde sismique ambiant - Application à la tomographie de l'Ouest de la France

The continuous excitation of the Earth’s surface due to oceanic wave-seafloor coupling, is widely used for both analyzing the seismic ambient wavefield properties and for imaging the deep structures. A new statistical approach based on the phase redundancy of cross-correlations is developed to provide quantitative information on the ubiquitous seismic signal. Two dense arrays of broadband seismometers have been temporarily deployed over the southern and western France, and over the northern Spain, providing a good opportunity to reveal the crustal and uppermost mantle features beneath this region. A Monte-Carlo Markov chain inversion algorithm is used to translate, for each station pair, the frequency dependent energy diagram of the ambient noise cross-correlation, into probability density distribution of the shear wave velocity as a function of depth. The 1D velocity profiles are then combined together using an anisotropic regionalization scheme, to build a 3D velocity model with an unprecedented lateral resolution of 75 km. The isotropic velocity maps correlate well with the geology up to 10 km depth. While negative velocity anomalies are associated with the sedimentary basins, the crystalline Variscan basement displays positive anomalies. At 20 km depth, the Bay of Biscay exhibits positive velocity anomalies, whose eastern boundaries can beinterpreted as the ocean-continent transition. The overall crustal structure below the Armorican Massif appears to be heterogenous at the subregional scale.L'énergie sismique générée constamment dans les océans peut être utilisée pour l'analyse du champ d'onde ambiant et pour l'imagerie des structures crustales. Une nouvelle approche statistique basée sur la redondance en phase des intercorrélations est développée pour caractériser le champ d'onde sismique. Le déploiement inédit d'un réseaude sismomètres à large-bande passante dans le Nord de l'Espagne, le Sud et l'Ouest de la France, ainsi que l'utilisation de stations permanentes européennes, permet une couverture sismique inédite de ces régions parintercorrélation des enregistrements continus. La méthode d'imagerie utilise un algorithme de type Monte-Carlo par chaînes de Markov. Pour chaque paire de capteurs, le diagramme de dispersion de la fonction de Green empirique est transcrit en densité de probabilité de vitesse sismique cisaillante en fonction de la profondeur. Une régionalisationanisotrope permet d'obtenir un modèle de vitesse 3D sous le Massif Armoricain et le Golfe de Gascogne avec une résolution latérale d'environ 75 km. Depuis la surface jusqu'à 10 km de profondeur, les cartes de vitesse isotropes sont très bien corrélées avec la géologie : des anomalies négatives sont observées au niveau des bassins sédimentaires, et des anomalies positives au niveau des socles cristallins affleurants, fragments de la chaîne hercynienne. À 20 km de profondeur, la structure à l'aplomb du Golfe de Gascogne se singularise par une anomalie de vitesse positive dont les limites peuvent être interprétées comme la transition entre les domaines océaniques et continentaux. La structure crustale du Massif Armoricain apparaît également hétérogène à l'échelle sub-régionale