Éditeur(s) : HAL CCSD Oxford University Press (OUP)
Résumé : International audience
We analyse P-wave receiver functions across the western Gulf of Aden and southern Red Sea continental margins in Western Yemen to constrain crustal thickness, internal crustal structure and the bulk seismic velocity characteristics in order to address the role of magmatism, faulting and mechanical crustal thinning during continental breakup. We analyse teleseismic data from 21 stations forming the temporary Young Conjugate Margins Laboratory (YOCMAL) network together with GFZ and Yemeni permanent stations. Analysis of computed receiver functions shows that (1) the thickness of unextended crust on the Yemen plateau is ∼35km; (2) this thins to ∼22km in coastal areas and reaches less than 14km on the Red Sea coast, where presence of a high-velocity lower crust is evident. The average Vp/Vs ratio for the western Yemen Plateau is 1.79, increasing to ∼1.92 near the Red Sea coast and decreasing to 1.68 for those stations located on or near the granitic rocks. Thinning of the crust, and by inference extension, occurs over a ∼130-km-wide transition zone from the Red Sea and Gulf of Aden coasts to the edges of the Yemen plateau. Thinning of continental crust is particularly localized in a <30-km-wide zone near the coastline, spatially co-incident with addition of magmatic underplate to the lower crust, above which on the surface we observe the presence of seaward dipping reflectors (SDRs) and thickened Oligo-Miocene syn-rift basaltic flows. Our results strongly suggest the presence of high-velocity mafic intrusions in the lower crust, which are likely either synrift magmatic intrusion into continental lower crust or alternatively depleted upper mantle underplated to the base of the crust during the eruption of the SDRs. Our results also point towards a regional breakup history in which the onset of rifting was synchronous along the western Gulf of Aden and southern Red Sea volcanic margins followed by a second phase of extension along the Red Sea margin.
ISSN: 0956-540X

hal-00825194
https://hal.archives-ouvertes.fr/hal-00825194
https://hal.archives-ouvertes.fr/hal-00825194/document
https://hal.archives-ouvertes.fr/hal-00825194/file/RF_Hakim2012_GJI.pdf
DOI : 10.1093/gji/ggt072

Éditeur(s) : HAL CCSD AGU and the Geochemical Society
Résumé : International audience
New magnetotelluric and receiver transfer function studies provide insights from the upper to the lower crust of the eastern Betic Cordillera, which is deformed by large folds, normal faults, and a major transcurrent left-lateral fault, the Carboneras Fault Zone (CFZ). Receiver function analysis determines a NNW dipping Moho reaching 20 degrees that increases in depth, from 20 km south of the CFZ up to 34 km in the Sierra de Los Filabres. In addition, seismic discontinuities determined in the upper crust are interpreted as major contacts between metamorphic complexes that are detached and folded. The MT inversion model reveals a conductive zone, also representing a crustal seismic discontinuity, associated with the Alpujarride/Nevado-Filabride contact and fitting the N vergent geometry of the Sierra Alhamilla antiform. A small flexure at Moho coincides with the CFZ, as revealed by the Bouguer anomaly trend, in agreement with the receiver function results. Moreover, the Bahr strike and tipper angle at the stations placed closest to the CFZ clearly reveal the continuity of the CFZ at least down to approximately 15 km in depth, crossing all the detected crustal discontinuities up to the Moho. The lack of a clear Moho offset associated with the Carboneras Fault supports the idea that some large strike-slip faults tend to accommodate the deformation by a broadening fault zone at lower crustal levels. Its nucleation could occur at the base of a thin crust, where melting processes critically reduced the lithospheric strength during the late Miocene, to then propagate upward, reaching the topographic surface. Northward, the lithosphere comprised moderately larger strength, and the crustal discontinuities favored the development of larger folds with kilometric amplitude instead of strike-slip faults since the late Miocene.
ISSN: 1525-2027

hal-00671024
https://hal.archives-ouvertes.fr/hal-00671024
DOI : 10.1029/2010GC003274

Éditeur(s) : HAL CCSD Elsevier
Résumé : International audience
The 12 May 2008 the Sichuan earthquake was one of the worst natural disasters in China. Here, we present a detailed picture of the epicentral area inferred geophysical data. Based on the analysis of teleseismic data acquired by a dense seismic network, we highlight an abrupt 20 km Moho offset between the Sichuan Basin and the Tibetan plateau and a horizontal discontinuity at similar to 15 km depth, which may connect with the ruptured zone of the 2008 event. We obtain low mean crustal velocity ratios, which suggest the absence of a thick and extensive zone of partial melt within the crust beneath the eastern part of the Songpan-Garze Terrane. All our data support the idea that the Longmen Shan range mark a zone of active strain localization due to the rigid Yangtze craton resisting eastward displacement of the Tibetan Plateau.
ISSN: 0040-1951

hal-00534471
https://hal.archives-ouvertes.fr/hal-00534471
DOI : 10.1016/j.tecto.2009.11.010

Éditeur(s) : HAL CCSD Oxford University Press (OUP)
Résumé : International audience
Rifting in a cratonic lithosphere is strongly controlled by several interacting processes including crust/mantle rheology, magmatism, inherited structure and stress regime. In order to better understand how these physical parameters interact, a 2 yr long seismological experiment has been carried out in the North Tanzanian Divergence (NTD), at the southern tip of the eastern magmatic branch of the East African rift, where the southward-propagating continental rift is at its earliest stage. We analyse teleseismic data from 38 broad-band stations ca. 25 km spaced and present here results from their receiver function (RF) analysis. The crustal thickness and Vp/Vs ratio are retrieved over a ca. 200 × 200 km2 area encompassing the South Kenya magmatic rift, the NTD and the Ngorongoro-Kilimanjaro transverse volcanic chain. Cratonic nature of the lithosphere is clearly evinced through thick (up to ca. 40 km) homogeneous crust beneath the rift shoulders. Where rifting is present, Moho rises up to 27 km depth and the crust is strongly layered with clear velocity contrasts in the RF signal. The Vp/Vs ratio reaches its highest values (ca. 1.9) beneath volcanic edifices location and thinner crust, advocating for melting within the crust. We also clearly identify two major low-velocity zones (LVZs) within the NTD, one in the lower crust and the second in the upper part of the mantle. The first one starts at 15–18 km depth and correlates well with recent tomographic models. This LVZ does not always coexist with high Vp/Vs ratio, pleading for a supplementary source of velocity decrease, such as temperature or composition. At a greater depth of ca. 60 km, a mid-lithospheric discontinuity roughly mimics the step-like and symmetrically outward-dipping geometry of the Moho but with a more slanting direction (NE–SW) compared to the NS rift. By comparison with synthetic RF, we estimate the associated velocity reduction to be 8–9 per cent. We relate this interface to melt ponding, possibly favouring here deformation process such as grain-boundary sliding (EAGBS) due to lithospheric strain. Its geometry might have been controlled by inherited lithospheric fabrics and heterogeneous upper mantle structure. We evidence that crustal and mantle magmatic processes represent first order mechanisms to ease and locate the deformation during the first stage of a cratonic lithospheric breakup.
ISSN: 0956-540X

hal-01622078
https://hal.archives-ouvertes.fr/hal-01622078
https://hal.archives-ouvertes.fr/hal-01622078/document
https://hal.archives-ouvertes.fr/hal-01622078/file/plasmanGJI2017.pdf
DOI : 10.1093/gji/ggx177

Éditeur(s) : HAL CCSD
Résumé : International audience
The Corinth Rift is the most seismically active zone in Europe. Since the year 2000, the Corinth Rift Laboratory (CRL, http://crlab.eu) consisting in a multidisciplinary natural observatory, aims at understanding the mechanics of faulting and earthquake nucleation in the Rift. Recent studies have improved our view about fault geometry and mechanics within CRL but there is still a critical need for a better knowledge of the structure at depth.In this project, we aim to analyze the complete seismological database (13 years of recordings) of CRL by using recently developed methodologies of structural imaging in order to determine, with high resolution, the local 3D structure and the earthquake locations. We perform an iterative joint determination of 3D velocity model and earthquake coordinates. P and S velocity models are estimated through a first arrival traveltime tomography by using the adjoint state method proposed by Taillandier et al. (2009). It consists in the minimization of the cost function between observed and theoretical arrival times,which is achieved by estimating the gradient of the cost function. Earthquake locations are also estimated with a linearized inversion approach by computing the gradient of the misfit function. Iterations are repeated until the cost function no longer decreases. We present preliminary results consisting in: (1) the adjustment of a 1D velocity model that is used as initial model of the 3D tomography and (2) a first attempt of the joint determination of 3D velocity model and earthquake locations.We also perform a preliminary receiver function analysis of teleseismic data recorded by the broadband stations of the CRL network. The multiscale approach (Gesret et. al., 2010) analysis of teleseismic P to S converted waves allows to analyze the signal toward the higher frequencies and thus to characterize the structure at depth with a high resolution. Our multiscale approach should give us reliable depth of interfaces beneath seismometers in the CRL area. In this first attempt, we adjust the 1D velocity model that produces a synthetic RF as similar as possible to the stack of observed RF for a subset of data. The identified interfaces could then be compared with heterogeneities imaged by the tomography.
AGU Fall Meeting
San Fransico, United States

hal-01417184
https://hal-mines-paristech.archives-ouvertes.fr/hal-01417184

É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 AGU and the Geochemical Society
Résumé : International audience
We image the lithospheric and upper asthenospheric structure of western continental Yemen with 24 broadband stations to evaluate the role of the Afar plume on the evolution of the continental margin and its extent eastwards along the Gulf of Aden. We use teleseismic tomography to compute relative P-wave velocity variations in southwestern Yemen down to 300 km depth. Published receiver function analysis suggest a dramatic and localized thinning of the crust in the vicinity of the Red Sea and the Gulf of Aden, consistent with the velocity structure that we retrieve in our model. The mantle part of the model is dominated by the presence of a low velocity anomaly in which we infer partial melting just below thick Oligocene flood basalts and recent off-axis volcanic events (from 15 Ma to present). This low velocity anomaly could correspond to an abnormally hot mantle and could be responsible for dynamic topography and recent magmatism in western Yemen. Our new P-wave velocity model beneath western Yemen suggests the young rift flank volcanoes beneath margins and on the flanks of the Red Sea rift are caused by focused small-scale diapiric upwelling from a broad region of hot mantle beneath the area. Our work shows that relatively hot mantle, along with partial melting of the mantle, can persist beneath rifted margins after breakup has occurred.
ISSN: 1525-2027

hal-01015496
https://hal.archives-ouvertes.fr/hal-01015496
https://hal.archives-ouvertes.fr/hal-01015496/document
https://hal.archives-ouvertes.fr/hal-01015496/file/Korostelev_etal_GGG_2014.pdf
DOI : 10.1002/2014GC005316

É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 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 American Geophysical Union
Résumé : International audience
The northern Canadian Cordillera (NCC) is an active orogenic belt in northwestern Canada characterized by deformed autochtonous and allochtonous structures that were emplaced in successive episodes of convergence since the Late Cretaceous. Seismicity and crustal deformation are concentrated along corridors located far (>200 to ~800 km) from the convergent plate margin. Proposed geodynamic models require information on crust and mantle structure and strain history, which are poorly constrained. We calculate receiver functions using 66 broadband seismic stations within and around the NCC and process them to estimate Moho depth and P-to-S velocity ratio (Vp/Vs) of the Cordilleran crust. We also perform a harmonic decomposition to determine the anisotropy of the subsurface layers. From these results, we construct simple seismic velocity models at selected stations and simulate receiver function data to constrain crust and uppermost mantle structure and anisotropy. Our results indicate a relatively flat and sharp Moho at 32 ± 2 km depth and crustal Vp/Vs of 1.75 ± 0.05. Seismic anisotropy is pervasive in the upper crust and within a thin (~10–15 km thick) sub-Moho layer. The modeled plunging slow axis of hexagonal symmetry of the upper crustal anisotropic layer may reflect the presence of fractures or mica-rich mylonites. The subhorizontal fast axis of hexagonal anisotropy within the sub-Moho layer is generally consistent with the SE-NW orientation of large-scale tectonic structures. These results allow us to revise the geodynamic models proposed to explain active deformation within the NCC.
ISSN: 2169-9313

Droits : info:eu-repo/semantics/OpenAccess
hal-01685550
https://hal.archives-ouvertes.fr/hal-01685550
https://hal.archives-ouvertes.fr/hal-01685550/document
https://hal.archives-ouvertes.fr/hal-01685550/file/Tarayoun_et_al-2017-Journalsur Geophysical_Research__Solid_Earth.pdf
DOI : 10.1002/2017JB014284

Éditeur(s) : HAL CCSD Elsevier
Résumé : International audience
Seismic anisotropy originating within the continental crust is used to determine kinematic flow lines within active mountain belts and is widely attributed to regionally aligned mica. However, naturally deformed micaceous rocks commonly show composite (e.g. S-C) fabrics. It is necessary therefore to understand how both varying mica content and differing intensities of multiple foliations impact on seismic interpretations in terms of deformation fields. An outcrop analogue for granitic mid-crustal deformed zones is used here to calibrate the seismic response against both parameters. Seismic responses are modelled using crystallographic preferred orientations for polymineralic, micaceous granitic gneisses, measured using Electron Back-Scatter Diffraction. The sample results are generalised by modelling the effects of variations in modal composition and the relative importance of deformation fabrics of variable orientation, so-called rock and fabric recipes. The maximum P- and S-anisotropy are calculated at 16.6% and 23.9% for single-foliation gneisses but for mixed (i.e. S-C) foliation gneisses these values reduce to 5.8% and 7.5% respectively. Furthermore, mixtures of multiple foliations generate significant variations in the geometry of the seismic anisotropy. This effect, coupled with the geographical orientation of fabrics in nature. can generate substantial variations in the orientation and magnitude of seismic anisotropy (especially for shear waves) as measured for the continental crust using existing receiver function and teleseismic near-vertical incidence methods. Thus, maps of seismic anisotropy varying with depth in deforming continents need not imply necessarily depth-varying deformation kinematics and tectonic decoupling.
ISSN: 0012-821X

hal-00445247
https://hal.archives-ouvertes.fr/hal-00445247
DOI : 10.1016/j.epsl.2009.09.035

Éditeur(s) : HAL CCSD Elsevier
Résumé : International audience
Nitrogen (N) transfer from legume trees to associated crops is a key factor for the N economy of low-input agroforestry systems. In this work, we presented a new approach to estimate N transfer based on the 15N content of root exudates and N released by root turnover of the donor plant (Gliricidia sepium) and the temporal change of the 15N content of the receiver plant (Dichanthium aristatum). The study was carried out in greenhouse using two isotopic methods: 15N leaf feeding (LF) and the natural 15N abundance (NA). Measurements of exudate 15N were made at several dates before and after tree pruning. A time-dependent box model was devised to quantify N transfer in time and to make comparisons between the isotopic methods. In NA, although tree roots and exudates presented a similar 15N signature before tree pruning, exudates were strongly depleted in 15N after pruning. In LF, exudates were always depleted in 15N in relation to tree roots. Hence, the current assumption used in N transfer studies concerning the equal 15N/14N distribution in tissues of the donor plant and in its excreted N was not confirmed in our study. Before pruning, N transfer functioned as a two-N-source system (soil N and exudates N) and both isotopic methods provided similar estimates: 11–12% for LF and 10–15% for NA. Calculations performed with the model indicated that N transfer occurred with small or nil fractionation of 15N in exudates. After pruning, there was a third N source associated with N released from tree root turnover. During this period, the isotopic signature of the receiver plant showed a transient state due to the progressive decrease of 15N content of that N source. The amount of N derived from the tree represented 65% of the total N content of the grass at the end of the experiments.
ISSN: 1161-0301

hal-01082183
https://hal-agroparistech.archives-ouvertes.fr/hal-01082183
DOI : 10.1016/j.eja.2007.03.003

Éditeur(s) : HAL CCSD Elsevier
Résumé : International audience
The lithospheric mantle is intrinsically heterogeneous and anisotropic. These two properties govern the repartition of deformation, controlling intraplate strain localization and development of new plate boundaries. Geophysical and geological observations provide clues on the types, ranges, and characteristic length scales of heterogeneity and anisotropy in the lithospheric mantle. Seismic tomography points to variations in geothermal gradient and hence in rheological behavior at scales of hundreds of km. Seismic anisotropy data substantiate anisotropic physical properties consistent at scales of tens to hundreds of km. Receiver functions imply lateral and vertical heterogeneity at scales < 10 km, which might record gradients in composition or anisotropy. Observations on naturally deformed peridotites establish that compositional heterogeneity and Crystal Preferred Orientations (CPOs) are ubiquitous from the mm to the km scales. These data allow discussing the processes that produce/destroy heterogeneity and anisotropy and constraining the time scales over which they are active. This analysis highlights: (i) the role of deformation and reactive percolation of melts and fluids in producing compositional and structural heterogeneity and the feedbacks between these processes, (ii) the weak mechanical effect of mineralogical variations, and (iii) the low volumes of fine-grained microstructures and difficulty to preserve them. In contrast, olivine CPO and the resulting anisotropy of mechanical and thermal properties are only modified by deformation. Based on this analysis, we propose that strain localization at the plate scale is, at first order, controlled by large-scale variations in thermal structure and in CPO-induced anisotropy. In cold parts of the lithospheric mantle , grain size reduction may contribute to strain localization, but the low volume of fine-grained domains limits this effect.
ISSN: 0040-1951

hal-01277174
https://hal.archives-ouvertes.fr/hal-01277174
DOI : 10.1016/j.tecto.2015.07.026

Éditeur(s) : HAL CCSD Oxford University Press (OUP)
Résumé : International audience
Our knowledge and understanding of the 3-D lithospheric structure of the Himalayas and the Tibetan Plateau is still challenging although numerous geophysical studies have been performed in the region. The GOCE satellite mission has the ambitious goal of mapping Earth's gravity field with unprecedented precision (i.e. an accuracy of 1-2 mGal for a spatial resolution of 100 km) to observe the lithosphere and upper-mantle structure. Consequently, it gives new insights in the lithospheric structure beneath the Himalayas and the Tibetan Plateau. Indeed, the GOCE gravity data now allow us to develop a new strategy for joint gravimetry-seismology inversion. Combined with teleseismic data over a large region in a joint inversion scheme, they will lead to lithospheric velocity-density models constrained in two complementary ways. We apply this joint inversion scheme to the Hi-CLIMB (Himalayan-Tibetan Continental Lithosphere during Mountain Building) seismological network which was deployed in South Tibet and the Himalayas for a 3-yr period. The large size of the network, the high quality of the seismological data and the new GOCE gravity data set allow us to image the entire lithosphere of this active area in an innovative way. We image 3-D low velocity and density structures in the middle crust that fit the location of discontinuous low S-velocity zones revealed by receiver functions in previous geophysical studies. In the deeper parts of our velocity model we image a positive anomaly interpreted to be the heterogenous Indian lithosphere vertically descending beneath the centre of the Tibetan Plateau.
ISSN: 0956-540X

hal-00858024
https://hal.archives-ouvertes.fr/hal-00858024
DOI : 10.1093/gji/ggs110

Éditeur(s) : HAL CCSD
Résumé : International audience
The Corinth Rift is the most seismically active zone in Europe. The area is characterized by very localized NSextension at a rate of1.5cm/year, the occurrence of frequent and intensive microseismic crises and occasionalmoderate to large earthquakes like in 1995 (Mw=6.1).Since the year 2000, the Corinth Rift Laboratory (CRL, http://crlab.eu) consisting in a multidisciplinary naturalobservatory, aims at understanding the mechanics of faulting and earthquake nucleation in the Rift. Recent studieshave improved our view about fault geometry and mechanics within CRL, but there is still a critical need for a betterknowledge of the structure at depth both for the accuracy of earthquake locations and for mechanical interpretationof the seismicity.In this project, we aim to analyze the complete seismological database (13 years of recordings) of CRL by usingrecently developed methodologies of structural imaging, in order to determine at the same time and with highresolution, the local 3D structure and the earthquake locations.We perform an iterative joint determination of 3D velocity model and earthquake coordinates. In a first step, Pand S velocity models are determined using first arrival time tomography method proposed by Taillandier et al.(2009). It consists in the minimization of the cost function between observed and theoretical arrival times which isachieved by the steepest descent method (e.g. Tarantola 1987). This latter requires computing the gradient of thecost function by using the adjoint state method (Chavent 1974). In a second step, earthquakes are located in thenew velocity model with a non-linear inversion method based on a Bayesian formulation (Gesret et al. 2015). Step1 and 2 are repeated until the cost function no longer decreases.We present preliminary results consisting in: (1) the adjustement of a 1D velocity model that is used as initialmodel of the 3D tomography and (2) a first attempt of the joint determination of 3D velocity model and earthquakelocation.In addition to the tomographic imaging, we perform a preliminary receiver function analysis of teleseismic datarecorded by the broadband stations of the CRL network. The RF analysis should provide the interface depthsbeneath seismometers and increase the imaging resolution of the upper crustal structures provided by the 3Dtomography. In this first attempt, we adjust the 1D velocity model that produces a synthetic RF as similar aspossible to the observed RF for a subset of data. We compare the identified interfaces with structures imaged bythe tomography.
EGU general assembly 2016
Vienne, Austria

hal-01306224
https://hal-mines-paristech.archives-ouvertes.fr/hal-01306224

Éditeur(s) : HAL CCSD Oxford University Press (OUP)
Résumé : International audience
Energy dissipation is observed on seismic data when a wave propagates through a porous medium, involving different frequency regimes depending on the nature of rock and fluid types. We focus here on the role of partial fluid saturation in unconsolidated porous media, looking in particular at P-wave phase velocity and attenuation. The study consists in running an experiment in a sand-filled tank partially saturated with water. Seismic propagation in the tank is generated in the kHz range by hitting a steel ball on a granite plate. Seismic data are recorded by buried accelerometers and injecting or extracting water controls the partial saturation. Several imbibition/drainage cycles were performed between the water and gas residual saturations. A Continuous Wavelet Transform applied on seismic records allowed us to extract the direct P-wave at each receiver.We observe an hysteresis in phase velocities and inverse quality factors between imbibition and drainage. Phase velocities and inverse quality factors are then jointly inverted to get a final poro-visco-elastic model of the partially saturated sand that satisfactorily reproduces the data. The model formulation consists in generalizing the Biot theory to effective properties of the fluid and medium (permeability and bulk modulus) in order to properly explain the phase velocity variation as a function of the saturation. The strong level of attenuation measured experimentally is further explained by an anelastic effect due to grain to grain sliding, adding to Biot's losses. This study shows that fluid distribution at microscopic scale has strong influence on the attenuation of direct P-waves at macroscopic scale and confirms that seismic prospection may be a powerful tool for the characterization of transport phenomena in porous media.
ISSN: 0956-540X

hal-00779920
https://hal.archives-ouvertes.fr/hal-00779920
https://hal.archives-ouvertes.fr/hal-00779920/document
https://hal.archives-ouvertes.fr/hal-00779920/file/Barriere_et_al_HAL.pdf
DOI : 10.1111/j.1365-246X.2012.05691.x

Éditeur(s) : HAL CCSD Geological Society of London
Résumé : International audience
Low-velocity zones ('bright spots') imaged by the INDEPTH seismic experiment in southern Tibet are extensively interpreted as widespread partial melt within the crust, which has given a strong support for the channel flow model. These suggest that a continuous seismic low-velocity zone underlies Tibet on the large scale. Here we take advantage of the Hi-CLIMB seismic experiment which includes a dense south-north profile and a lateral 2D seismic network to assess the vertical and the horizontal extension of low-velocity zones in southern Tibet. Several approaches including migration, amplitude analysis and waveform inversion of receiver functions are performed to detect crustal low-velocity zones using this new seismological dataset. Our results reveal localized and discontinuous low-velocity zones in Tibet. They indicate that the vertical extension of the low-velocity zones is about 10 km, and their maximum horizontal length appears to be c. 50 km. Our study suggests a partial correlation between the location of these low-velocity zones and the spatial distribution of Tibetan grabens. These results, especially the non-continuity of low-velocity zones, together with the observed regular value of mean crustal VP/VS ratio, question the existence of widespread partial melt of the southern Tibetan crust and, therefore, the viability of the channel flow model.
Growth and Collapse of the Tibetan Plateau

hal-00853971
https://hal.archives-ouvertes.fr/hal-00853971
DOI : 10.1144/SP353.6

Éditeur(s) : HAL CCSD Seismological Society of America
Résumé : International audience
On 9 October 2007, an unusual sequence of earthquakes began in central British Columbia about 20 km west of the Nazko cone, the most recent (circa 7200 yr) volcanic center in the Anahim volcanic belt. Within 25 hr, eight earthquakes of magnitude 2.3-2.9 occurred in a region where no earthquakes had previously been recorded. During the next three weeks, more than 800 microearthquakes were located (and many more detected), most at a depth of 25-31 km and within a radius of about 5 km. After about two months, almost all activity ceased. The clear P- and S-wave arrivals indicated that these were high-frequency (volcanic-tectonic) earthquakes and the b value of 1.9 that we calculated is anomalous for crustal earthquakes but consistent with volcanic-related events. Analysis of receiver functions at a station immediately above the seismicity indicated a Moho near 30 km depth. Precise relocation of the seismicity using a double-difference method suggested a horizontal migration at the rate of about Graphic, with almost all events within the lowermost crust. Neither harmonic tremor nor long-period events were observed; however, some spasmodic bursts were recorded and determined to be colocated with the earthquake hypocenters. These observations are all very similar to a deep earthquake sequence recorded beneath Lake Tahoe, California, in 2003-2004. Based on these remarkable similarities, we interpret the Nazko sequence as an indication of an injection of magma into the lower crust beneath the Anahim volcanic belt. This magma injection fractures rock, producing high-frequency, volcanic-tectonic earthquakes and spasmodic bursts.
ISSN: 0037-1106

hal-00854504
https://hal.archives-ouvertes.fr/hal-00854504
DOI : 10.1785/0120100013

Éditeur(s) : HAL CCSD Elsevier
Résumé : International audience
This paper provides a synthesis of current data and interpretations on the crustal structure of the Pyrenean-Cantabrian orogenic belt, and presents new tectonic models for representative transects. The Pyrenean orogeny lasted from Santonian (~84 Ma) to early Miocene times (~20 Ma), and consisted of a spatial and temporal succession of oceanic crust/exhumed mantle subduction, rift inversion and continental collision processes at the Iberia-Eurasia plate boundary. A good coverage by active-source (vertical-incidence and wide-angle reflection) and passive-source (receiver functions) seismic studies, coupled with surface data have led to a reasonable knowledge of the present-day crustal architecture of the Pyrenean-Cantabrian belt, although questions remain. Seismic imaging reveals a persistent structure, from the central Pyrenees to the central Cantabrian Mountains, consisting of a wedge of Eurasian lithosphere indented into the thicker Iberian plate, whose lower crust is detached and plunges northwards into the mantle. For the Pyrenees, a new scheme of relationships between the southern upper crustal thrust sheets and the Axial Zone is here proposed. For the Cantabrian belt, the depth reached by the N-dipping Iberian crust and the structure of the margin are also revised.The common occurrence of lherzolite bodies in the northern Pyrenees and the seismic velocity and potential field record of the Bay of Biscay indicate that the precursor of the Pyrenees was a hyperextended and strongly segmented rift system, where narrow domains of exhumed mantle separated the thinned Iberian and Eurasian continental margins since the Albian-Cenomanian. The exhumed mantle in the Pyrenean rift was largely covered by a Mesozoic sedimentary lid that had locally glided along detachments in Triassic evaporites. Continental margin collision in the Pyrenees was preceded by subduction of the exhumed mantle, accompanied by the pop-up thrust expulsion of the off-scraped sedimentary lid above. To the west, oceanic subduction of the Bay of Biscay under the North Iberian margin is supported by an upper plate thrust wedge, gravity and magnetic anomalies, and 3D inclined sub-crustal reflections. However, discrepancies remain for the location of continent-ocean transitions in the Bay of Biscay and for the extent of oceanic subduction. The plate-kinematic evolution during the Mesozoic, which involves issues as the timing and total amount of opening, as well as the role of strike-slip drift, is also under debate, discrepancies arising from first-order interpretations of the adjacent oceanic magnetic anomaly record.
ISSN: 0040-1951

insu-01683989
https://hal-insu.archives-ouvertes.fr/insu-01683989
https://hal-insu.archives-ouvertes.fr/insu-01683989/document
https://hal-insu.archives-ouvertes.fr/insu-01683989/file/teixell-tectono-2018.pdf
DOI : 10.1016/j.tecto.2018.01.009

É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