Éditeur(s) : HAL CCSD Springer-Verlag
Résumé : The seismology is the most effective method to explore the Structure of subduction zones to great depth. The distinguishing feature of the mantle in the subduction regions is the presence of hydrated phases, which transport water into the Earth's interior and release it with dramatic local consequences, triggering earthquakes and melting. The seismological detection of these hydrous phases and geo-dynamic interpretation of flow in the hydrous mantle depend on knowledge of the anisotropic elastic properties and the characteristics of the wave propagation in anisotropic media. In this paper we briefly recall the distinguishing features of anisotropic wave propagation and the observable parameters. We suggest the Vp/Vs ratio is a physically sound parameter than can be observed by seismology in anisotropic regions of the Earth, whereas the Poisson's ratio, which is often quoted, is not directly observable and does not correspond to the characteristics of wave propagation in an anisotropic or isotropic medium. We report for the first time the ratios of Vp/Vs1 and Vp/Vs2, where Vs1 is the fastest and Vs2 slowest S-wave velocities of an anisotropic media. We present the current knowledge of the anisotropic seismic properties of hydrous minerals in the upper mantle, transition zone, and lower mantle that are stable along low temperature geotherms associated with subduction, and identify which minerals are likely to influence seismological observations because they have very high volume fractions, or very high anisotropies, or both of these.;In the upper mantle antigorite and talc are exceptionally anisotropic (Vp 71% Vs 68% and Vp 65%, Vs 68%, respectively) and chlorite is also very anisotropic for S-waves (Vp 35%, Vs 76%). Comparatively less anisotropic are hornblende (Vp 27%, Vs 31%), used as proxy for tremolite in subduction zones, and clinohumite (21.8%, Vs 15.9%). Brucite at 4 GPa (Vp 26.5%, Vs 30.9%) and the dense hydrogen magnesium silicate (DHMS) phase A at 9GPa (Vp 9.3%, Vs 17.6%) are the only hydrous minerals stable in the upper mantle that have had their elastic properties measured at in situ mantle pressures. Except for the phase A, all these minerals are more anisotropic than olivine for at least one parameter (Vp, Vs, Vp/Vs1, or Vp/Vs2). In the transition zone the major phases hydrous wadsleyite and ringwoodite C have moderate to weak anisotropies (Vp 16.3%, Vs 16.5%, and Vp 1.9%, Vs 4.4%, respectively). The DHMS Superhydrous B has a moderate anisotropy (Vp 6.9%, Vs 11.6%). At greater depth the DHMS phase D Lit 24 GPa (Vp 10.8%, Vs 18.0%) is the only hydrous phase that can transport hydrogen from the transition zone into the lower mantle. The phase D is more anisotropic for S than P waves like many of the hydrous phases. From these data it is clear that hydrous phases are in general very anisotropic. However, pressure can play a strong role in reducing anisotropy. It is the case for brucite and talc, in which increasing pressure from ambient to 4 GPa reduces the anisotropy by about 50% for both P and S waves. In contrast, the anisotropy of the DHMS phase A does not change significantly with pressure.;Our picture of the seismic anisotropy of hydrated minerals remains incomplete; the elastic properties of many have not been measured even at ambient conditions (e.g., 10 angstrom phase and phase E) or not measured in their true elastic symmetry (e.g., clinochlore). The majority of hydrous minerals have not been measured at high pressure. and none related to hydrated mantle have been measured at elevated temperature. We have shown in the few cases where hydrous minerals have been measured as a function of pressure, that this variable has an important effect oil the velocity distribution and in most cases reduces the degree of anisotropy, hence we Would expect seismic anisotropy to play a key role in the determination of the shallow structure of subduction zones in the upper mantle.
Subduction zone geodynamics

hal-00420893
https://hal.archives-ouvertes.fr/hal-00420893

Éditeur(s) : HAL CCSD Oxford University Press (OUP)
Résumé : Gravity data and P-wave teleseismic traveltime residuals from 29 temporary broad-band stations spread over the northern margin of the Gulf of Aden (Dhofar region, Oman) were used to image lithospheric structure. We apply a linear relationship between density and velocity to provide consistent density and velocity models from mid-crust down to about 250 km depth. The accuracy of the resulting models is investigated through a series of synthetic tests. The analysis of our resulting models shows: (1) crustal heterogeneities that match the main geological features at the surface; (2) the gravity edge effect and disparity in anomaly depth locations for layers at 20 and 50 km; (3) two low-velocity anomalies along the continuation of Socotra-Hadbeen and Alula-Fartak fracture zones between 60 and 200 km depth; and (4) evidence for partial melting (3-6 per cent) within these two negative anomalies. We discuss the presence of partial melting in terms of interaction between the Sheba ridge melts and its along-axis segmentation.
ISSN: 0956-540X

hal-00456048
https://hal.archives-ouvertes.fr/hal-00456048
DOI : 10.1111/j.1365-246X.2009.04438.x

Éditeur(s) : HAL CCSD European Geosciences Union
Résumé : This study presents a unique database of 172 plagioclase Crystallographic Preferred Orientations (CPO) of variously deformed gabbroic rocks. The CPO characteristics as a function of the deformation regime (magmatic or crystal-plastic) are outlined and discussed. The studied samples are dominantly from slow- and fast-spread present-day ocean crust, as well as from the Oman ophiolite. Plagioclase is the dominant mineral phase in the studied samples. Plagioclase CPOs are grouped into three main categories: Axial-B, a strong point alignment of (010) with a girdle distribution of [100]; Axial-A, a strong point maximum concentration of [100] with parallel girdle distributions of (010) and (001); and P-type, point maxima of [100], (010), and (001). A majority of CPO patterns are Axial-B and P-type, in samples showing either magmatic or crystal-plastic deformation textures. Axial-A CPOs are less common; they represent 21% of the samples deformed by crystal-plastic flow. Although fabric strength (ODF J index) does not show any consistent variation as a function of the CPO patterns, there is a significant difference in the relationship between the ODF and pole figures J indices; the magmatic type microstructures have high (010) pole figures J indices, which increase linearly with ODF J index, whereas the high [100] pole figures J indices of plastically deformed samples vary in a more scattered manner with ODF J index. The multistage nature of plastic deformation superposed on a magmatic structure compared with magmatic flow, and the large number of possible slip-systems in plagioclase probably account for these differences. Calculated seismic properties (P wave and S wave velocities and anisotropies) of plagioclase aggregates show that anisotropy (up to 12% for P wave and 14% for S wave) tends to increase as a function of ODF J index. In comparison with the olivine 1998 CPO database, the magnitude of P wave anisotropy for a given J index is much less than olivine, whereas it is similar for S wave anisotropy. Despite a large variation of fabric patterns and geodynamic setting, seismic properties of plagioclase-rich rocks have similar magnitudes of anisotropy. There is a small difference in the aggregate elastic symmetry, with magmatic microstructures having higher orthorhombic and hexagonal components, whereas plastic deformation microstructures have a slightly higher monoclinic component, possibly correlated with predominant monoclinic simple shear flow in plastically deformed samples. Overall, plots for CPO strength (ODF J index), pole figure strength, CPO symmetry and seismic anisotropy show significant scattering. This could be related to sampling statistics, although our database is a factor of ten higher than the olivine database of 1998, or it could be related to the low symmetry (triclinic) structure of plagioclase resulting in the addition of degrees of freedom in the processes creating the CPOs.
ISSN: 1869-9510

hal-00963439
https://hal.archives-ouvertes.fr/hal-00963439
DOI : 10.5194/se-4-511-2013

Éditeur(s) : HAL CCSD
Résumé : poster
International audience
Locations for earthquake recorded in the Lesser Antilles subduction zone are processed separately by regional observatories, NEIC and ISC. There is no earthquake location catalog available compiling all available phase arrival data. We propose a new best complete earthquake catalog by merging all available phase arrival data for better constrains on earthquake locations. ISC provides the phase arrival data of 29243 earthquakes (magnitude range from 1.4 to 6.4) recorded by PRSN (Porto Rico), SRC (British West Indies), and from FUNVISIS (Venezuela). We add phases data from IPGP observatories for 68718 earthquakes from magnitudes 0.1 to 7.5 (OVSG, Guadeloupe, recorded 53226 earthquakes since 1981, and OVSM, Martinique, recorded 29931 earthquakes since 1972). IPGP also provides the accelerometer waveform data of the GIS-RAP network. We achieved automatic picking on the GIS-RAP data using the Component Energy Correlation Method. The CECM provides high precision phase detection, a realistic estimation of picking error and realistic weights that can be used with manual pick weights. The CECM add an average of 3 P-waves and 2 S-waves arrivals to 3846 earthquakes recorded by the GIS-RAP network since 2002. The final catalog contains 84979 earthquakes between 1972 and 2013, 24528 of which we compiled additional data. We achieve earthquake location using NonLinLoc, regional P and S waves data and a set of one dimensional velocity models. We produce improved locations for 22974 earthquakes (better residuals, on equal or larger arrival dataset) and improved duration magnitudes for 6258 earthquakes (using duration data and improved locations). A subset of best constrained 15626 hypocenters (with more than 8 phases and an average RMS of 0.48_0.77s) could be used for structural analysis and earthquake local tomography. Relative locations are to be applied in order to image active faulting. We aim to understand coupling in the seismogenic zone as well as triggering mechanisms of intermediate depth seismicity like fluid migration beneath the volcanic arc.
Geophysical Research Abstracts
Vienne, Austria

hal-01059207
https://hal-brgm.archives-ouvertes.fr/hal-01059207

Éditeur(s) : HAL CCSD Elsevier
Résumé : International audience
Joint inversions are now commonly used in the earth sciences. They have been developed to better understand the structure of the earth, since they provide more constraints on the inverted parameters. We propose a new process to simultaneously invert several data sets in order to better image 3D crustal and upper mantle structures. Our inversion uses three kinds of data that present good complementarity: (1) P-wave receiver functions to provide Moho depth variations, (2) teleseismic delay times of P-waves to retrieve velocity anomalies in the crust and the upper mantle, and (3) gravity anomalies to image density variations at the lithospheric scale. We use a stochastic scheme, where receiver functions are first inverted. The interpolated resulting Moho depths are incorporated as a priori information into the joint inversion of teleseismic delay times and gravity anomalies process. Moreover, velocity and density can be linked by empirical relationships, which justifies the joint inversion of those parameters. In our stochastic approach, we perform a model space search for Moho variations, P-velocity, and density structure to find an acceptable fit to the three data sets. In order to preferentially sample the good data fit regions, we chose the neighborhood algorithm of Sambridge to optimistically survey the model space. We model the delay times with 3D raytracing using evenly spaced velocity-density nodes. We present here the first results given by this method on synthetic tests. (C) 2011 Elsevier Ltd. All rights reserved.
ISSN: 0098-3004

hal-00643158
https://hal.archives-ouvertes.fr/hal-00643158
DOI : 10.1016/j.cageo.2010.11.017

Éditeur(s) : HAL CCSD American Geophysical Union
Résumé : International audience
We report measurements of long-period P wave polarization (P pol) in Australia and Tahiti made by combining modeling of the polarization deviation and harmonic analysis. The analysis of the deviation of the horizontal polarization of the P wave as a function of event back azimuth may be used to obtain information about (1) sensor misorientation, (2) dipping discontinuities, (3) seismic anisotropy, and (4) velocity heterogeneities beneath a seismic station. The results from harmonic analysis and a grid search using Snell's law suggest the presence of a dipping seismic discontinuity beneath stations CTAO and CAN in Australia. These results are consistent with published receiver function studies for these stations. The P pol fast axis orientation is close to the N–S absolute plate motion direction at station TAU (Tasmania), which may be due to plate-motion-driven alignment of olivine crystals in the asthenosphere. Interestingly, measurements of SKS splitting at Tahiti (French Polynesia) show an apparent isotropy, whereas an inversion of P pol observations at PPTL seismic station located in Tahiti suggests the presence of two anisotropic layers. The fast axis azimuth is oriented E–W in the upper layer, and it is close to the NW–SE orientation in the lower layer. Since P pol orientations are used for real-time earthquake locations, especially in poorly instrumented areas such as the South Pacific, we show that the bias from anisotropy and sensor misorientation determined here can be corrected to improve the location accuracy, which yields fundamental data for rapid location necessary for effective tsunami warning.
ISSN: 2169-9313

hal-01389113
http://hal.univ-reunion.fr/hal-01389113
http://hal.univ-reunion.fr/hal-01389113/document
http://hal.univ-reunion.fr/hal-01389113/file/Fontaine_Ppol_JGR2009.pdf
DOI : 10.1029/2008JB005709

Éditeur(s) : HAL CCSD Oxford University Press (OUP)
Résumé : International audience
To quantify the seismic properties of lower crustal rocks and to better constrain the origin of the lower crustal seismic reflectivity, we determined the complete 3-D seismic properties of a lower crustal section. Eight representative samples of the main lithologic and structural units outcropping in the Val Sesia (Ivrea zone) were studied in detail. The seismic velocities were calculated using the single crystal stiffness coefficients and the lattice preferred orientation (LPO) of each mineral in all samples. The 21 stiffness coefficients characterizing the elastic behaviour of each rock are determined. Mafic and ultramafic rocks such as pyroxenite and pyroxene-bearing gabbros display complex shear wave properties. These rocks are weakly birefringent (maximum 0.1 kms-1) and it is difficult to find consistent relationships between the seismic properties and the rock structure. On the other hand, seismic properties of deformed felsic rocks are essentially controlled by mica. They display strong S-wave birefringence (0.3 km s-1) and relatively high Vp anisotropy (7.6 per cent). Amphibole also strongly influences the rock birefringence patterns. For both kind of rocks, the foliation is highly birefringent and the fast polarized shear wave is systematically oriented parallel to the foliation. We show that the number of mineral phases in the rock strongly controls the anisotropy. The seismic anisotropy has a complex role in the P-wave reflectivity. Compared to the isotropic case, anisotropy enhances the reflection coefficient for about 60 per cent of the possible lithological interfaces. For 40 per cent of the interfaces, the reflection coefficient is much lower when one considers the medium as anisotropic.
ISSN: 0956-540X

hal-01389122
http://hal.univ-reunion.fr/hal-01389122
http://hal.univ-reunion.fr/hal-01389122/document
http://hal.univ-reunion.fr/hal-01389122/file/Barruol_3D_velocities_GJI1993.pdf
DOI : 10.1111/j.1365-246X.1993.tb01519.x

Éditeur(s) : HAL CCSD
Résumé : Joint inversions are now commonly used in Earth Sciences. They have been developed to better understand Earth structure since they provide more constraints on the inverted parameters. We propose a new process to simultaneously invert several data sets in order to better image 3-D crustal and upper mantle structures. Our inversion uses 3 kind of data which present a good complementarity : (1) P receiver functions to provide Moho depth variations, (2) teleseismic delay times of P waves to retrieve velocity anomalies in the crust and upper mantle, (3) gravity anomalies to image density variations at lithospheric scale. In our scheme, receiver functions are first inverted. The resulting Moho depths are then interpolated and incorporated as a priori information into the joint inversion of teleseismic delay times and gravity anomalies process. In our new approach, we perform a model space search for Moho variations, P-velocity and density structure to find acceptable fit to the three data sets. In order to preferentially sample the good data fit region, we chose the Neighborhood algorithm of Sambridge to optimistically survey the model space. We model the delay times with a 3-D raytracing using velocity nodes evenly spaced and linked with density nodes via a linear relationship. This relationship can evolve through depth, simulating its pressure dependance. We test our scheme on synthetic examples, and apply the method in a continental rifting region where deformation processes are complex and badly known. This joint inversion thus provides more constraints on the inverted parameters, and allows us to better image crust - mantle interactions in the case of the Hangai dome (Central Mongolia). Indeed, in this area, we expect an asthenospheric upwelling to dynamically support the topography, but we image a density - velocity anomaly in the lower crust. The effects on topography and the origin (composition vs. thermal) of this anomaly are still under debate.
Les inversions conjointes sont maintenant communément utilisées en Sciences de la Terre. Elles ont été développées afin d'améliorer notre connaissance et compréhension de la structure interne de la Terre puisqu'elles apportent de plus en plus de contraintes sur les paramètres inversés. Dans cette thèse, nous proposons un nouveau processus d'inversion conjointe qui prend en compte trois paramètres différents et qui mène à l'obtention d'un modèle lithosphérique tridimensionnel. Notre méthode utilise trois types de données gravimétriques et sismologiques qui présentent une bonne complémentarité : (1) Les fonctions récepteur P afin d'obtenir les variations de profondeur du Moho, (2) les délais de temps d'arrivées P des téléséismes pour retrouver les anomalies de vitesse dans la croûte et le manteau supérieur, (3) les anomalies gravimétriques qui donnent accès aux variations de densité à l'échelle lithosphérique. Dans le schéma d'inversion proposé, nous inversons d'abord les fonctions récepteur. Les variations de profondeur du Moho qui en résultent sont alors interpolées puis incorporées comme information a priori dans le processus d'inversion conjointe gravimétrie - tomographie télésismique. Que ce soit pour les fonctions récepteur ou pour la partie conjointe, notre méthode d'inversion se base sur un algorithme de recherche afin de trouver la structure qui minimise l'écart aux trois types de données. Nous avons opté pour l'algorithme de voisinage qui présente l'avantage de concentrer sa recherche dans les régions qui minimisent l'écart aux données sans toutefois abandonner les autres régions afin de limiter les risques de convergence vers un minimum local. Nous utilisons un modèle 3-D constitué de noeuds de vitesse afin de modéliser les délais de temps. Ces noeuds de vitesse sont liés aux noeuds de densité grâce à une relation linéaire entre vitesse et densité. Cette relation peut évoluer avec la profondeur de manière à simuler les effets dus aux variations de pression. Nous avons réalisé plusieurs tests synthétiques afin d'évaluer le comportement de notre nouvelle méthode. Ils ont permis de montrer que les variations du Moho ainsi que les anomalies de vitesse et densité sont bien retrouvées. Ils ont également montré l'importance de la prise en compte de ces trois types de données dans un même schéma d'inversion, notamment en ce qui concerne la faculté de discerner les variations d'interfaces des anomalies volumiques. Nous avons ensuite appliqué la méthode dans une région continentale en phase de rifting où les processus de déformation sont complexes et mal connus. Notre nouvelle méthode a apporté davantage de contraintes sur les paramètres inversés et a permis d'obtenir un modèle correctement résolu de la croûte et du manteau supérieur sous le dôme de Hangai (Mongolie Centrale). Une anomalie asthénosphérique sous le dôme est souvent invoquée pour expliquer sa topographie élevée mais c'est une anomalie (en vitesse et densité) crustale (croûte inférieure) qui s'exprime dans nos modèles. Ses effets sur la topographie et son origine (compositionnelle et/ou thermique) sont encore débattus et nécessiteraient des études plus poussées (paramétriques, VP =VS, ...).
https://tel.archives-ouvertes.fr/tel-00783524

tel-00783524
https://tel.archives-ouvertes.fr/tel-00783524
https://tel.archives-ouvertes.fr/tel-00783524/document
https://tel.archives-ouvertes.fr/tel-00783524/file/Basuyau2011.pdf

Éditeur(s) : HAL CCSD
Résumé : The purpose of this work is to specify the deformation offshore East of Taiwan by characterizing major active faults in particular at the level of the Ryukyu forearc. The major problem associated with the study of this region is that earthquakes are located outside the permanent networks. Consequently, this led to some unprecise earthquake locations limiting tectonic and statistical analyses, and therefore the correct estimation of the seismic hazard. This work is included in the project ANR ACTS (Active Tectonics and Seismic Hazard in Taiwan). To achieve this goal, we worked on three different timescales for which a different localization method has been proposed to obtain a clear image of seismic deformation and highlight major faults offhore East of Taiwan. First, at the scale of the century (from the analysis of a seismicity catalogue with homogeneous magnitude), the instrumental historical earthquakes (1897-2007) show that the Ryukyu margin was affected by 6 earthquakes of magnitude greater than 7 with 4 of them for which the responsible fault is unknown. We propose a method of relative location to relocate the historic events from the time difference of P and S arrival times (Ts-Tp) at each station. The location is obtained by searching for "analogue" earthquakes in the recent instrumental catalogue (1991-2008) for which Ts-Tp at each station is close. This method has been applied to the largest (in terms of magnitude) earthquake ever recorded in Taiwan, the earthquake of June 5, 1920 (M7.7 +/- 0.2). A rupture along the interplate of the Ryukyu subduction with a possible nucleation at the downdip limit of a splay-fault is certainly responsible of this earthquake. Then, at the scale of the decade, the improvement of the method of absolute location based on MAXimum of Intersection of EDT (MAXI, EDT: Equal Differential Time) allows to better extract the erroneous arrival times and prevents the "trade-off" between depth and origin time, and between depth and epicentral position when the azimuthal coverage is low (gap Azimuthal > 180degre). Synthetic tests show the effectiveness of MAXI (using P-wave only) to determine the parameters x, y and z even when the azimuthal gap is important. In the latter case, the quality of the results is dependent on the velocity model to represent the 3D structure of the Earth. We propose an approach involving the use of a 3D a priori P-wave velocity model to locate earthquakes which are lateral and remote to the seismic network. An application to the case of Taiwan and Ecuador validates this approach. Finally, at the level of a few months, the marine geophysical campaign RATS (Ryukyu Arc: Tectonics and Seismology) was conducted in two stages. A passive seismology experiment from July to October, 2008 (RATS1) has been conducted over the forearc of the Ryukyus and then active seismic experiment (refraction and reflection) was conducted in May 2009 in a NNE - SSW line through the Ryukyu margin. These two combined experiments allow improving our knowledge about the crustal structure of the margin. At the level of the forearc, the basement of the forearc is characterized by vertical backstop and a very deformed basis certainly associated with a significant out of sequence deformation. In depth, the downgoing plate is certainly affected by a tear that controls the seismicity in the region of transition between the subduction and collision.
Le but de ce travail est de préciser la déformation en mer à l'est de Taïwan et de notamment caractériser les failles actives majeures en particulier au niveau de l'avant-arc des Ryukyus. Le problème majeur associé à l'étude de cette région est que les séismes sont localisés en dehors des réseaux permanents et que cela conduit à des localisations peu précises limitant les analyses statistiques, tectoniques et donc l'estimation correcte de l'aléa sismique. Ce travail est inclus dans le projet ANR ACTS (Active Tectonics and Seismic Hazard in Taiwan). Pour atteindre cet objectif, nous travaillons à trois échelles de temps différentes pour lesquelles une méthode de localisation différente a été proposée pour obtenir une image précise de la déformation sismique et mettre en évidence les failles majeures en mer à l'est de Taïwan. D'abord, à l'échelle du siècle (à partir de l'analyse d'un catalogue de sismicité homogène en magnitude), les séismes historiques instrumentaux (1897-2007) montrent que la marge des Ryukyus a été affectée par 6 séismes de magnitude supérieure à 7 dont 4 pour lesquels la faille responsable n'est pas connue. Nous proposons une méthode de localisation relative qui permet de relocaliser les évènements historiques à partir des différences de temps d'arrivées P et S (Ts-Tp) à chaque station. La localisation est obtenue en recherchant les séismes analogues dans le catalogue instrumental récent (1991-2008) pour lesquels Ts-Tp à chaque station est proche. Cette méthode a été appliquée au plus gros séisme (en terme de magnitude) jamais enregistré à Taïwan, le séisme du 5 juin 1920 (M7.7 $pm$ 0.2). Ce séisme a certainement pour origine une rupture sur l'interplaque sismogène de la subduction des Ryukyus avec une possible nucléation à la base d'une faille hors-séquence. Ensuite, à l'échelle de la dizaine d'années, l'amélioration de la méthode de localisation absolue du MAXimum d'Intersection EDT (MAXI, EDT: Equal Differential Time) permet de mieux extraire les temps d'arrivées erronés et empêche les "trade-off" entre profondeur et temps d'origine et, entre profondeur et position épicentrale lorsque la couverture azimutale est faible (gap azimutal > 180°). Des tests synthétiques montrent l'efficacité de MAXI (en utilisant les ondes P uniquement) à déterminer les paramètres x, y et z même lorsque le gap azimutal est important. Dans ce dernier cas, la qualité des résultats reste dépendant du modèle de vitesse à représenter la structure 3D de la Terre. Nous proposons une démarche associant l'utilisation d'un modèle de vitesse des ondes P 3D a priori avec la méthode MAXI pour localiser les séismes qui sont latéraux et éloignés du réseau. Une application au cas de Taïwan et en Équateur valide cette approche. Enfin, à l'échelle de quelques mois, la campagne de géophysique marine RATS (Ryukyu Arc: Tectonics and Seismology) a été menée en deux temps. Une expérience de sismologie passive de juillet à octobre 2008 (RATS1) a été menée au-dessus de l'avant-arc des Ryukyus puis une expérience de sismique active (réfraction et réflexion) a été menée en mai 2009 sur une ligne NNE-SSW à travers la marge des Ryukyus. Ces deux expériences combinées permettent d'améliorer notre connaissance de la structure crustale de la marge. Au niveau de l'avant-arc, le socle de l'avant-arc est caractérisé par une extrémité formant un butoir vertical et une base très déformée certainement associée à une importante déformation hors-séquence. En profondeur, la plaque plongeante est certainement affectée par une déchirure qui contrôle la sismicité dans cette région de transition entre la subduction et la collision.
https://tel.archives-ouvertes.fr/tel-00660737

tel-00660737
https://tel.archives-ouvertes.fr/tel-00660737
https://tel.archives-ouvertes.fr/tel-00660737/document
https://tel.archives-ouvertes.fr/tel-00660737/file/Thesis_small.pdf

Éditeur(s) : HAL CCSD Oxford University Press (OUP)
Résumé : International audience
To improve earthquake location, we create a 3-D a priori P-wave velocity model (3-DVM) that approximates the large velocity variations of the Ecuadorian subduction system. The 3-DVM is constructed from the integration of geophysical and geological data that depend on the structural geometry and velocity properties of the crust and the upper mantle. In addition, specific station selection is carried out to compensate for the high station density on the Andean Chain. 3-D synthetic experiments are then designed to evaluate the network capacity to recover the event position using only P arrivals and the MAXI technique. Three synthetic earthquake location experiments are proposed: (1) noise-free and (2) noisy arrivals used in the 3-DVM, and (3) noise-free arrivals used in a 1-DVM. Synthetic results indicate that, under the best conditions (exact arrival data set and 3-DVM), the spatiotemporal configuration of the Ecuadorian network can accurately locate 70 per cent of events in the frontal part of the subduction zone (average azimuthal gap is 289° ± 44°). Noisy P arrivals (up to ± 0.3 s) can accurately located 50 per cent of earthquakes. Processing earthquake location within a 1-DVM almost never allows accurate hypocentre position for offshore earthquakes (15 per cent), which highlights the role of using a 3-DVM in subduction zone. For the application to real data, the seismicity distribution from the 3-D-MAXI catalogue is also compared to the determinations obtained in a 1-D-layered VM. In addition to good-quality location uncertainties, the clustering and the depth distribution confirm the 3-D-MAXI catalogue reliability. The pattern of the seismicity distribution (a 13 yr record during the inter-seismic period of the seismic cycle) is compared to the pattern of rupture zone and asperity of the Mw = 7.9 1942 and the Mw = 7.7 1958 events (the Mw = 8.8 1906 asperity patch is not defined). We observe that the nucleation of 1942, 1958 and 1906 events coincides with areas of positive Simple Bouguer anomalies and areas where marine terraces are still preserved on the coastal morphology. From north to south: (1) the 1958 rupture zone is almost aseismic and is attributed to a zone of high coupling; (2) south of the Galera alignment (perpendicular to the trench), the 1942 rupture zone presents moderate seismicity, deeper on the seismogenic interplate zone, and abutting on the Jama cluster (to the south). This cluster is facing the Cabo Pasado cap and positive Bouguer anomalies on the overriding margin. We suspect that this cluster reflects a zone of local asperity (partial coupling). South of the Jama cluster, the spherical aseismic zone in the Bahia area is interpreted as having a low seismic coupling (steady creep motion or slow slip events). We suspect that the site that generated the three M > 7 events (1896, 1956 and 1998) correspond to a small patch of strong coupling. To the south, in the Manta-Puerto Lopez zone, the seismicity is mainly organized in earthquake swarms (1998, 2002, 2005). Although slow slip events have been observed in the area (Vallée et al. submitted), we infer from the coastline shape, the marine terraces and the high positive Bouguer anomalies that the seismicity here might reveal a significant amount of seismic coupling.
ISSN: 0956-540X

hal-00873478
https://hal.archives-ouvertes.fr/hal-00873478
DOI : 10.1093/gji/ggs083

Éditeur(s) : HAL CCSD Oxford University Press (OUP)
Résumé : Geophysical Journal International, v. 173, n. 2, p. 473-482, 2008. http://dx.doi.org/10.1111/j.1365-246X.2008.03718.x
International audience
The western Provence in southern France is an intraplate low deforming region, cut by large sinistral strike-slip faults as the Nîmes Fault. The deformation rate of this fault was estimated at 0.1 mm yr−1 from geological and morphological observations. Nevertheless, some large earthquakes occurred in this area in historical times inducing well-documented partial destructions on the roman monument 'Pont-du-Gard'. In order to investigate if this area is seismically active, we installed a temporary seismological network of 14 stations between the cities of Nîmes and Avignon for a 1-yr period (2002 July–2003 June). We recorded and located a total amount of 80 earthquakes of magnitude Ml smaller than 2.8, and of 153 quarry blasts. The seismic event locations are interpreted in terms of regional tectonics as well as of geometry at depth of the Nîmes Fault and of possible seismic activity of the Roquemaure Fault. At the beginning of 2002 September, a catastrophic storm occurred in the same area with cumulated precipitations at 600 mm in 28 hr, causing casualties and major inundations and damages. We identified a clear and sudden increase of the seismic activity as an immediate consequence of this exceptional meteorological phenomenon. We interpret this rainfall triggering of earthquakes as the response of the crust to an abrupt vertical loading. The geological context of this zone is characterized by the presence of sedimentary basins and of an important karstic network in a calcareous environment allowing the retention of a third of the total rainwater volume during few weeks after the catastrophic storm. Month by month inversion of the seismic wave velocities for only the first 1-km thick layer of the crust allows investigating the influence of the stored water. A clear decrease of both P- and S-velocities, and an increase of the Vp/Vs ratio are evidenced. The variations obtained are at 2, 6 and 4.5 per cent for P-, S-velocities and Vp/Vs ratio, respectively. These variations are discussed in confront of previous studies and can be directly related to the known water volume.
ISSN: 0956-540X

hal-00407884
https://hal.archives-ouvertes.fr/hal-00407884
https://hal.archives-ouvertes.fr/hal-00407884v2/document
https://hal.archives-ouvertes.fr/hal-00407884/file/rigoGJI2008.pdf
DOI : 10.1111/j.1365-246X.2008.03718.x

Éditeur(s) : HAL CCSD Elsevier
Résumé : International audience
The Gibraltar arc, spans a complex portion of the Africa-Eurasia plate boundary marked by slow oblique convergence and intermediate and deep focus seismicity. The seemingly contradictory observations of a young extensional marine basin surrounded by an arcuate fold-and-thrust belt, have led to competing geodynamic models (delamination and subduction). Geophysical data acquired in the past decade provide a test for these models and support a narrow east-dipping, subduction zone. Seismic refraction studies indicate oceanic crust below the western Gulf of Cadiz. Tomography of the upper mantle reveals a steep, east-dipping high P-wave velocity body, beneath Gibraltar. The anisotropic mantle fabric from SKS splitting shows arc-parallel "fast directions", consistent with toroidal flow around a narrow, westward retreating subducting slab. The accompanying WSW advance of the Rif-Betic mountain belt has constructed a thick pile of deformed sediments, an accretionary wedge, characterized by west-vergent thrust anticlines. Bathymetric swath‐mapping images an asymmetric embayment at the deformation front where a 2 km high basement ridge has collided. Subduction has slowed significantly since 5 Ma, but deformation of recent sediments and abundant mud volcanoes suggest ongoing activity in the accretionary wedge. Three possible origins for this deformation are discussed; gravitational spreading, overall NW-SE convergence between Africa and Iberia and finally a WSW tectonic push from slow, but ongoing roll-back subduction. In the absence of arc volcanism and shallow dipping thrust type earthquakes, evidence in favor of present-day subduction can only be indirect and remains the object of debate. Continued activity of the subduction offers a possible explanation for great (M > 8.5) earthquakes known to affect the area, like the famous 1755 Great Lisbon earthquake. Recent GPS studies show SW motion of stations in N Morocco at velocities of 3-6 mm/yr indicating the presence of an independent block, a "Rif-Betic-Alboran" microplate, situated between Iberia and Africa.
ISSN: 0040-1951

insu-00769762
https://hal-insu.archives-ouvertes.fr/insu-00769762
DOI : 10.1016/j.tecto.2012.08.038

Éditeur(s) : HAL CCSD Elsevier
Résumé : International audience
Imaging deep active volcanic areas remains a challenge in our understanding of their activity and evolution, especially in subduction zones. Study of magnetic anomalies is appropriate to access such dynamics in depth. The magnetic anomaly pattern of the Lesser Antilles Arc (LAA) subduction is studied through Curie Point Depth (CPD), interpreted as the depth of the 580 °C isotherm, and developed to better assess the deep thermal structure of the arc. The depth of the estimated CPD exhibits a complex topography. Keeping in mind the overall uncertainty associated with this method, a main doming is evidenced below the Guadeloupe archipelago. Its apex is shifted towards the ancient arc, suggesting a very hot state of the fore-arc/arc domain. To better understand the LAA thermal state, we perform 2D thermo-mechanical simulations of the subduction zone. Recalling that magnetite is a serpentinization by-product, we simulate water transfer triggered by slab dehydration to test the assumption of fore-arc serpentinization suggested by the positive magnetic anomaly in the vicinity of the Guadeloupe archipelago. In this area, the subduction-induced arc lithosphere hydration and related weakening trigger a fast heating of the upper plate by basal convective removal. This process of fast arc lithosphere thinning may apply where simultaneously the volcanic arc is split in two and normal convergence is high enough. As serpentinization strongly decreases P-wave velocity, we propose a new interpretation of a published seismic profile below Guadeloupe. The seismic layer previously interpreted as the arc lower crust may rather be a layer of serpentinized mantle, as supported by spatial correlations between gravimetric and magnetic anomalies. Consequently, at the scale of Guadeloupe Island, the fore-arc Moho would be shallower than initially assumed, with a dome shape more consistent with both the extensive deformation active since the Oligocene in the inner fore-arc and the CPD doming.
ISSN: 0040-1951

Droits : http://creativecommons.org/licenses/by-nc-nd/
insu-01505300
https://hal-insu.archives-ouvertes.fr/insu-01505300
https://hal-insu.archives-ouvertes.fr/insu-01505300/document
https://hal-insu.archives-ouvertes.fr/insu-01505300/file/Gailler-2017.pdf
DOI : 10.1016/j.tecto.2017.03.014

Éditeur(s) : HAL CCSD American Geophysical Union
Résumé : International audience
We report measurements of long-period P wave polarization (P-pol) in Australia and Tahiti made by combining modeling of the polarization deviation and harmonic analysis. The analysis of the deviation of the horizontal polarization of the P wave as a function of event back azimuth may be used to obtain information about (1) sensor misorientation, (2) dipping discontinuities, (3) seismic anisotropy, and (4) velocity heterogeneities beneath a seismic station. The results from harmonic analysis and a grid search using Snell's law suggest the presence of a dipping seismic discontinuity beneath stations CTAO and CAN in Australia. These results are consistent with published receiver function studies for these stations. The P-pol fast axis orientation is close to the N-S absolute plate motion direction at station TAU (Tasmania), which may be due to plate-motion-driven alignment of olivine crystals in the asthenosphere. Interestingly, measurements of SKS splitting at Tahiti (French Polynesia) show an apparent isotropy, whereas an inversion of P-pol observations at PPTL seismic station located in Tahiti suggests the presence of two anisotropic layers. The fast axis azimuth is oriented E-Win the upper layer, and it is close to the NW-SE orientation in the lower layer. Since P-pol orientations are used for real-time earthquake locations, especially in poorly instrumented areas such as the South Pacific, we show that the bias from anisotropy and sensor misorientation determined here can be corrected to improve the location accuracy, which yields fundamental data for rapid location necessary for effective tsunami warning.
ISSN: 2169-9313

hal-00413141
https://hal.archives-ouvertes.fr/hal-00413141
https://hal.archives-ouvertes.fr/hal-00413141/document
https://hal.archives-ouvertes.fr/hal-00413141/file/FontaineetalJGR2008.pdf
DOI : 10.1029/2008JB005709