Éditeur(s) : HAL CCSD Elsevier
Résumé : International audience
Talc is a hydrous magnesium rich layered silicate that is widely disseminated in the Earth from the seafloor to over 100 km depth, in ultra-high pressure metamorphism of oceanic crust. In this paper we determine the single crystal elastic constants at pressures from 0 to 12 GPa of talc triclinic (C 1¯) and monoclinic (C2/c) polytypes using ab initio methods. We find that talc has an extraordinarily high elastic anisotropy at zero pressure that reduces with increasing pressure. The exceptional anisotropy is complemented by a negative Poisson's ratio for many directions in crystal space. Calculations show that talc is not only one of very few common minerals to exhibit auxetic behaviour, but the magnitude of this effect may be the largest reported so far for a mineral. The compression (Vp) and shear (Vs) wave velocity anisotropy is 80% and 85% for the triclinic polytype. At pressures where talc is known be stable in the Earth (up to 5 GPa) the Vp and Vs anisotropy is reduced to about 40% for both velocities, which is still a very high value. Vp is slow parallel to the c-axis and fast perpendicular to it. This remains unchanged with increasing pressure and is observed in both polytypes. The shear wave splitting (difference between fast and slow S-wave velocities) at low pressure has high values in the plane normal to the c-axis, with a maximum near the alow asterisk-axis in the triclinic and the b-axis in the monoclinic polytype. The c-axis is the direction of minimum splitting. The pattern of shear wave splitting does not change significantly with pressure. The volume fraction of talc varies between 11 and 41% for hydrated mantle rocks, but the lack of data on the crystallographic preferred orientation (CPO) precludes a detailed analysis of the impact of talc on seismic anisotropy in subduction zones. However, it is highly likely that CPO can easily develop in zones of deformation due to the platy habit of talc crystals. For random aggregates of talc, the isotropic Vp, Vs and Vp/Vs ratio have significantly lower values than those of antigorite and may explain low-velocity regions in the mantle wedge. Vp/Vs ratios are more complex in anisotropic media because there are fast and slow S-waves, resulting in Vp/Vs1 and Vp/Vs2 ratios for every propagation direction, making interpretation difficult in deformed polycrystalline talc with a CPO. Talc on the subduction plate boundary can strongly influence guided wave velocity as CPO would develop in this region of intense shearing. The very low coefficient of friction (< 0.1) of talc above 100 °C could also explain silent earthquakes at shallow depths (ca 30 km) along the subduction plate boundaries, frequently responsible for tsunami.
ISSN: 0012-821X

hal-00411481
https://hal.archives-ouvertes.fr/hal-00411481
DOI : 10.1016/j.epsl.2008.07.047

Éditeur(s) : HAL CCSD Springer Verlag (Germany)
Résumé : International audience
Surface displacements due to temporal changes in environmental mass redistributions are observable in the coordinate time series of many Global Navigation Satellite System (GNSS) sites. In this study, we investigated the effect of loading on estimates of tectonic velocity computed from campaign-style GNSS observations. The study region is in the Pyrenees mountain range between France and Spain (ResPyr campaigns). In this area, seismic activity is continuous and moderate and the expected amplitude of the horizontal tectonic velocity is less than 0.5 mm/yr. In order to determine the velocity, 4 sparse GNSS campaigns were carried out from 1995 to 2010. Considering this small rate of deformation, loading phenomena can contribute a non-negligible artifact to the velocity computation that could affect our geodynamical interpretation. In this investigation, we specifically considered the atmospheric, hydrological, and non-tidal ocean loading phenomena. The computed loading deformations for this region show the horizontal displacements are dominated by the non-tidal ocean loading (maximum 4 mm for the North and 3.1 mm for the East components); the main vertical contributions are due to the atmospheric and continental water storage loading (maximum 14.3 for the atmosphere and 8.1 mm for the hydrology, respectively). We have found that the dominant loading effect on the horizontal velocity is the non-tidal ocean loading (mean of 0.11 mm/yr), whereas the vertical component is dominated by the hydrological loading (mean of 0.21 mm/yr). Since the study area is in a mountainous region, we also analyzed the difference between the atmospheric and the topography dependent atmospheric loading models at our GNSS campaign sites. We did not find any significant difference between the two atmospheric loading models in terms of horizontal velocity. Finally, we performed simulations to identify the optimum timing and frequency of future GNSS campaigns in this area that would minimize the loading effects on tectonic velocity estimates.
ISSN: 0039-3169

hal-00963629
https://hal.archives-ouvertes.fr/hal-00963629
DOI : 10.1007/s11200-012-0458-2

Éditeur(s) : HAL CCSD American Geophysical Union (AGU) AGU and the Geochemical Society
Résumé : International audience
To study how an impacting plume modifies the mantle lithosphere, we analyzed the microstructures and crystal preferred orientations (CPO) of 29 peridotites and 37 pyroxenites that sample the mantle root of the Ontong Java Plateau (OJP) from 60 to 120 km depth. The peridotites show a strong compositional variability, but homogeneous coarse granular to tabular microstructures, except for those equilibrated at the shallowest and deepest depths, which are porphyroclastic. All peridotites have clear olivine CPO, with dominant fiber-[010] patterns. Low intragranular misorientations and straight grain boundaries in olivine suggest that, above 100 km depth, annealing often followed deformation. Calculated density and P wave velocities of the peridotites decrease weakly with depth. S wave velocities decrease faster, resulting in increasing Vp/Vs ratio with depth. Calculated densities and seismic velocity profiles are consistent with those estimated for normal mantle compositions under a cold oceanic geotherm. Enrichment in pyroxenites may further increase seismic velocities. The calculated seismic properties cannot therefore explain the low S waves velocities predicted by Rayleigh wave tomography and ScS data in the mantle beneath the OJP. Calculated P and S waves anisotropy is variable (2–12%). It is higher on average in the deeper section of the lithosphere. Because olivine has dominantly [010]-fiber CPO patterns, if foliations are horizontal, vertically propagating S waves and Rayleigh waves will sample very weak anisotropy in the OJP mantle lithosphere. Moreover, if the orientation of the lineation changes with depth, the anisotropy-induced contrast in seismic properties might produce an intralithospheric reflector marking the stratification of the OJP mantle root.
ISSN: 1525-2027

hal-01116428
https://hal.archives-ouvertes.fr/hal-01116428
DOI : 10.1002/2014GC005452

Éditeur(s) : HAL CCSD Elsevier
Résumé : International audience
Melt generation in the continental crust is thought to significantly influence seismic bulk velocities and anisotropy, although existing laboratory data provide limited constraints on such seismic attributes. In this study we measured ultrasonic compressional wave speeds (Vp) and anisotropy during sample compaction, mica-breakdown and melt-generation in synthetic, foliated quartz–muscovite aggregates. Measurements were performed at peak conditions of 300 MPa hydrostatic confining pressure and 750 °C, over a six hour time period, for three separate samples where wave propagation directions were aligned at 0°, 45° and 90° with respect to the foliation plane. The experiments are initially marked by sample compaction and rapid reduction in porosity (from 25% to <2 vol%), with a corresponding increase in Vp. During initial stages of the experiment Vp ranges from 4.65 to 5.45 km/s, depending on the direction of the P wave propagation, resulting in up to ∼16% anisotropy, which originates mainly from the preferred orientation of muscovite and possible incipient melt lenses oriented parallel to foliation. As the experiment progresses the bulk Vp decreases gradually and Vp anisotropy drops from ∼16% to ∼2% after five hours into the experiment. The drop in anisotropy is caused by muscovite-breakdown and increasing amount of partial melting of the aggregate. The bulk velocities are considerably lower than predicted for a mixture of muscovite and β-quartz (stable at 300 MPa and 750 °C). Thermodynamic modeling indicates that mica breakdown and production of melt with ∼2 wt% dissolved water in the aggregate can explain the low Vp, even though β-quartz is stable; the modeled Vp is ∼4.75 km/s at 300 MPa and 750 °C. These ultrasonic attributes signify movement of the melt to fill vacant pores and coating grain boundaries, effectively connecting melt films and reducing seismic anisotropy. The results may apply to regions where crustal melting is on-going, in particular regions of elevated temperature that are only weakly deforming, with coupled anomalously low bulk velocities and seismic anisotropy.
ISSN: 0012-821X

hal-01243569
https://hal.archives-ouvertes.fr/hal-01243569
DOI : 10.1016/j.epsl.2015.05.039

Éditeur(s) : HAL CCSD American Geophysical Union
Résumé : International audience
A set of synthetic aggregates, containing mixtures of calcite and muscovite, were compacted with varying uniaxial loads ranging from 20 MPa up to 400 MPa. Their elastic properties have been measured using compressional and shear waves in confining hydrostatic pressures up to 475 MPa. Measured seismic velocities are shown to depend on 1) the ratio of calcite to muscovite, 2) the uniaxial load used during sample manufacturing, and 3) the porosity amount prior to velocity measurements. The matrix framework may also affect the seismic velocities, but this effect is not easily quantifiable. In general, measured seismic velocities decrease with increasing muscovite content and porosity. Elastic properties have been calculated based on texture measurements, which were obtained by neutron diffraction. The calculated velocities are based on the calcite and muscovite single crystal elastic tensors, their orientation distribution functions, and their modal fractions. A large discrepancy is apparent between calculated and measured velocities, where the former always overestimate the actual velocities. Shear waves display less of a difference between calculated and measured values, than do compressional waves, indicating that waves propagating by shearing are less affected by sample porosity and matrix framework. A better agreement between measured and calculated seismic velocities arises when the initial porosity is taken into account, using a differential effective medium model. Seismic anisotropy is evident from both measured and calculated velocities, and is low for samples containing pure calcite, but becomes prominent as the muscovite concentration increases. The intensity of anisotropy further depends on the uniaxial load used during sample compaction and the initial porosity.
ISSN: 0148-0227

hal-00512666
https://hal.archives-ouvertes.fr/hal-00512666
DOI : 10.1029/2009JB006523

Éditeur(s) : HAL CCSD Oxford University Press (OUP): Policy P - Oxford Open Option A Oxford University Press (OUP)
Résumé : International audience
The Pyrenean mountain range is a slowly deforming belt with continuous and moderate seismic activity. To quantify its deformation field, we present the velocity field estimated from a GPS survey of the Pyrenees spanning 18 yr. The PotSis and ResPyr networks, including a total of 85 GPS sites, were installed and first measured in 1992 and 1995–1997, respectively, and remeasured in 2008 and 2010. We obtain a deformation field with velocities less than 1 mm yr−1 across the range. The estimated velocities for individual stations do not differ significantly from zero with 95 per cent confidence. Even so, we estimate a maximum extensional horizontal strain rate of 2.0 ± 1.7 nanostrain per year in a N–S direction in the western part of the range. We do not interpret the vertical displacements due to their large uncertainties. In order to compare the horizontal strain rates with the seismic activity, we analyse a set of 194 focal mechanisms using three methods: (i) the ‘r’ factor relating their P and T axes, (ii) the stress tensors obtained by fault slip inversion and (iii) the strain-rate tensors. Stress and strain-rate tensors are estimated for: (i) the whole data set, (ii) the eastern and western parts of the range separately, and (iii) eight zones, which are defined based on the seismicity and the tectonic patterns of the Pyrenees. Each of these analyses reveals a lateral variation of the deformation style from compression and extension in the east to extension and strike-slip in the west of the range. Although the horizontal components of the strain-rate tensors estimated from the seismic data are slightly smaller in magnitude than those computed from the GPS velocity field, they are consistent within the 2σ uncertainties. Furthermore, the orientations of their principal axes agree with the mapped active faults.
ISSN: 0956-540X

hal-01171752
https://hal.archives-ouvertes.fr/hal-01171752
DOI : 10.1093/gji/ggv052

Éditeur(s) : John Wiley & Sons
Résumé : Based on all available seismic refraction data, the seismic structure of the Western Approaches to the English Channel and the South Armorican continental shelf has been established. The correlation of seismic velocities showed a three-fold structure, related to geological layers as follows: 1. High-velocity, deep sequence (4.4-6.6km/sec): metamorphic basement and Lower to Middle Palaeozoic. 2. Intermediate sequence (3.6-4.3km/sec): Permo-Triassic (and Jurassic?) and possibly also Lower Cretaceous on the margin of the Western Approaches. 3. Upper sequence (1.9-3.6km/sec): Cretaceous and Cenozoic sediments. An acoustic unconformity appears to be present between the Upper sequence and the High velocity and Intermediate sequences. Different geomorphology characterizes these layers: the Palaeozoic is tectonized and generally follows the Hercynian structural trends; mainly Permo-Triassic sediments fill the depressions of the Palaeozoic floor, but only in the Western Approaches; eastwards-transgressive Cretaceous and Cenozoic sediments cover discordantly the preceding layers. [NOT CONTROLLED OCR]
Petroleum and the continental shelf of North-west Europe (John Wiley & Sons), 1975 , P. 29-43

Droits : John Wiley & Sons
http://archimer.ifremer.fr/doc/1975/publication-5465.pdf
http://archimer.ifremer.fr/doc/00000/5465/

Éditeur(s) : HAL CCSD Elsevier
Résumé : International audience
Seismic tomography and high surface heat flow models indicate an anomalously thin and hot lithospheric mantle beneath the Canadian Cordillera. To understand the thermal and deformational history of the lithospheric mantle under these conditions, we performed a study of the microstructures, crystal preferred orientations (CPO), and seismic anisotropy of mantle xenoliths collected in the Summit Lake volcano in the southern Canadian Cordillera. The xenoliths vary in composition between spinel lherzolites, cpx-rich lherzolites, and wehrlites. All samples are coarse-grained and show clear evidence for annealing, but deformation microstructures are well-preserved in the spinel lherzolites. Olivine CPO is characterized by a strong point-maximum of [100] and by a girdle distribution of [010[ and [001] axes normal to it; CPO strength is variable, spinel-lherzolites have the strongest CPO. Some samples show a weak orthopyroxene CPO consistent with the olivine one. However most pyroxenes have random CPO, interstitial shapes and undeformed appearance, which suggest a secondary origin by reactive melt percolation under static conditions, leading to crystallization of cpx + sp +/- opx at the expenses of olivine, followed by crystallization of cpx + sp +/- ol at the expenses of opx to produce the wehrlites. The studied xenoliths have typical upper mantle seismic properties patterns with fastest propagation of P-waves and polarization of S-waves parallel to the olivine [100] maxima. Maximum P- and S-waves anisotropy for lherzolites are 8 and 5% respectively. Dispersion of the olivine CPO due to neocrystallization of olivine results in weaker anisotropies (6.5 and 3.5%) in wehrlites. Seismic velocities calculated using a geotherm consistent with the equilibrium temperatures of the Summit Lake xenoliths are 3-4% lower than those expected in a normal, 100 km-thick lithosphere, consistent with seismic tomography models for the region. This high geothermal gradient may be partially explained by considering the anisotropy of thermal conductivity of olivine and subhorizontal foliations in the lithospheric mantle, which result in a decrease by up to 15% of the thermal conductivity in the vertical direction.
ISSN: 0040-1951

hal-00618332
https://hal.archives-ouvertes.fr/hal-00618332
DOI : 10.1016/j.tecto.2011.04.014

Éditeur(s) : HAL CCSD Elsevier
Résumé : International audience
To quantify the relationships between anisotropy. S-wave splitting and tectonics, we determined the seismic properties of lower-crustal and upper-mantle rocks outcropping in the lvrea Zone (Northem Italy). We obtained P-and S-wave seismic velocities by laboratory direct velocity measurements and/or by calculations based on the modal compositions of the rocks, the lattice preferred orientations (LPOs), and the single crystal stiffness coefficients. Measured P-and S-wave velocities (6.0-7 .5 km s-1 and 3.6-4.2 km s-1) are typical of the lower crust. The P-wave anisotropy is in the range 0-l 0%. Shear-wave birefringence is in the range 0.0-0.6 km s-1 , with typical values between 0.0 and 0.2 km s-1 • In many cases, the birefringence is clearly related to fabric elements (foliation, lineation). Mafic rocks such as anorthosite or pyroxene-bearing gabbros exhibit low P-wave anisotropies (< 5%) and low shear-wave birefringences (less than 0.1 km s-1). In contrast, the seismic properties of felsic rocks such as biotite-bearing gneisses and mafic rocks such as amphibolites display high V P anisotropy (10%) and strong birefringence (0.3 km s-1). Biotite and amphibole preferred orientations clearly control seismic anisotropy and particularly shear-wave splitting. In these rocks, maximum splitting is observed in directions parallel to the foliation with the fast split shear wave polarized parallel to the foliation plane. To have an overview of the seismic properties of this lower-crustal section at a broader scale, we calculated from our data the anisotropie seismic properties of several hypothetical samples that are perhaps more representative of the regional anisotropy than each sample individually. For instance, the average lower-crustal sample displays an anisotropy of 5.5% for P waves and a birefringence around 0.14 km s I for S waves propagating parallel to the foliation. We observe little splitting for waves propagating at high angle to the foliation.
ISSN: 0031-9201

hal-01389712
http://hal.univ-reunion.fr/hal-01389712
http://hal.univ-reunion.fr/hal-01389712/document
http://hal.univ-reunion.fr/hal-01389712/file/Barruol_Kern-PEPI96_hal.pdf
DOI : 10.1016/0031-9201(95)03124-3

Éditeur(s) : North-Holland Publishing Company, Amsterdam
Résumé : The Barents Sea is underlain by a thick (up to 5 km) sedimentary basin. Seismic refraction has outlined four main velocity discontinuities which have been correlated with geological units on the basis of the geological history of the region. The basin is underlain by a crystalline basement, the nature of which cannot be determined on the basis of seismic velocity alone. Metamorphosed Paleozoic units (velocities around 5-4 km/sec) lie over this basement. Their thickness is not well established but appears to reach 1 or 2 km in some cases. A very distinct and thick (up to 2.5 km) layer (4.1 km/sec) is found almost everywhere and is thought to correspond to a major discontinuity at the end ot the Paleozoic. This discontinuity is followed by a variable velocity layer (2.4-4 km/sec, up to 1.2 km thick) which includes Mesozoic and Cenozoic sediments and is limited to the south of the Barents Sea. The absence of thick Tertiary deposits support the idea that the Barents shelf was emergent at that time. Recent low-velocity sediments are found close to the shelf edge. [NOT CONTROLLED OCR]
Earth and Planetary Science Letters (North-Holland Publishing Company, Amsterdam), 1974 , Vol. 24 , P. 33-47

Droits : North-Holland Publishing Company, Amsterdam
http://archimer.ifremer.fr/doc/1974/publication-5447.pdf
http://archimer.ifremer.fr/doc/00000/5447/

Éditeur(s) : HAL CCSD American Geophysical Union
Résumé : International audience
The three-dimensional S-wave velocity structure beneath the South Pacific Superswell is obtained from joint broadband seismic experiments on the ocean floor and islands. We collected only approximately 800 relative times of long-period teleseismic SH-waves by using a waveform cross-correlation from 76 events occurring from January 2003 to May 2005. We conducted relative time tomography to obtain a 3D structure to depths of 1600 km. In the resultant image, we find a characteristic distribution of low- velocity regions. The most prominent features are a large doughnut-shaped low-velocity region at 800 km depth, and an elongated large low-velocity region beneath the Society to Pitcairn hotspots at 1200 km depth. Our model suggests that a large low shear velocity province rooted in the D00 extends upwards and culminates near the top of the lower mantle beneath the central part of the South Pacific Superswell although its perfect continuity is not still confirmed.
ISSN: 0094-8276

hal-01389114
http://hal.univ-reunion.fr/hal-01389114
http://hal.univ-reunion.fr/hal-01389114/document
http://hal.univ-reunion.fr/hal-01389114/file/Tanaka_PLUME_lower_mantle_GRL2009.pdf
DOI : 10.1029/2009GL037568

Éditeur(s) : HAL CCSD Geological Society of America
Résumé : International audience
Integrated Ocean Drilling Program Expedition 313 continuously cored uppermost Eocene to Miocene sequences on the New Jersey shallow shelf (Sites M27, M28, and M29). Previously, 15 Miocene (ca. 23-13 Ma) seismic sequence boundaries were recognized on several generations of multichannel seismic profi les using criteria of onlap, downlap, erosional truncation, and toplap. We independently recognize sequence boundaries in the cores and logs based on an integrated study of core surfaces, lithostratigraphy and process sedimentology (grain size, mineralogy, facies, and paleoenvironments), facies successions, stacking patterns, benthic foraminiferal water depths, downhole logs, core gamma logs, and chronostratigraphic ages. We use a velocitydepth function to predict the depths of seismic sequence boundaries that were tested by comparison with major core surfaces, downhole and core logs, and synthetic seismograms. Using sonic velocity (core and downhole), core density, and synthetic seismograms, we show that sequence boundaries correspond with acoustic impedance contrasts, although other stratal surfaces (e.g., maximum fl ooding and transgressive surfaces) also produce refl ections. Core data are suffi cient to link seismic sequence boundaries to impedance contrasts in 9 of 12 instances at Site M27, 6 of 11 instances at Site M28, and 8 of 14 instances at Site M29. Oligocene sequences have minimal lithologic and seismic expression due to deep-water locations on clinoform bottomsets. Miocene sequences (ca. 23-13 Ma) were sampled across several unconformity clinothems (prograding units) on topset, foreset, and bottomset locations. Excellent recovery allows core-seismic integration that confi rms the hypothesis that unconformities are a primary source of impedance contrasts. Our core-seismic-log correlations predict that key seismic surfaces observed in other sub surface investigations without core and/or well logs are stratal surfaces with sequence stratigraphic signifi cance.
EISSN: 1553-040X

insu-00946235
https://hal-insu.archives-ouvertes.fr/insu-00946235
DOI : 10.1130/GES00858.1

É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) : Université Nice Sophia Antipolis
Résumé : The Lesser Antilles is a case study of a very slow subduction (~2 cm/yr) of an old oceanic lithosphere (~84-100 Ma). The region presents a relatively low seismic activity, especially along the interplate contact, and the seismic hazard associated with a possible mega-thrust earthquake is still poorly known. This PhD thesis is a first step toward assessing the ability of the Lesser Antilles subduction zone to produce such a large subduction event. To do so, it aims at constraining the downdip width of the interplate's seismogenic zone. The lack of coverage of permanent seismological stations is a major limitation in the exploration of the Lesser Antilles subduction zone. It is due to the presence of only small aligned islands at far distances from the potentially seismogenic interplate area. Several oceanographic cruises were therefore planned that notably allowed the repeated deployment of ocean bottom seismometers; some of them being left for a few months of background seismicity recording. This thesis specifically focuses on two sets of wide-angle seismic data acquired offshore the Dominica and Martinique islands. From their analysis 3D and 2D tomographic models were produced respectively over the forearc region and across the whole subduction complex. These models constrain the plates' seismic structure as well as their geometry. They allow the discussion of how the imaged structures affect the subduction processes and give a first estimation of the downdip width of the seimogenic zone, defined as the segment of the interplate between the backstop and the upper plate's Moho. The joint interpretation of seismic models and earthquake localizations then refine this first assessment. Epicenter distribution from height months of seismic recording shows indeed that seismicity concentrates within the inner forearc region. The strong velocity gradient that characterize its basement suggests it is denser and more rigid than the more deformable outer forearc basement. The updip limit of the seismogenic zone could then lie arcward of the backstop at the contact of the interplate and the seaward end of the inner forearc crust. At depth, interplate earthquake mechanisms are observed between 35 and 45 km depth and interpreted to have occurred at the downdip limit of the seismogenic zone. The latter could reach a depth over 10,km deeper than the contact of the upper plate's Moho with the interplate, and therefore lie within the mantle wedge. All together, these results imply a large downdip width of the seismogenic zone (~70 km) offshore the Dominica and Martinique islands. Further work is, however, needed in order to fully comprehend the ability of the Lesser Antilles subduction zone to produce a possible mega-thrust earthquake. This would necessitate the evaluation of seismic coupling at the interplate contact and the possible segmentation of the seismogenic zone, for instance, due to the subduction of oceanic ridges.
Les Petites Antilles présentent un contexte géodynamique caractérisé par la subduction à très faible vitesse (~2 cm/an) d'une lithosphère océanique âgée (~84-100 Ma). L'activité sismique y est relativement faible, notamment à l'interplaque, où l'aléa sismique lié à un éventuel séisme de méga-chevauchement reste encore mal contraint. Cette thèse se veut être une première étape dans l'évaluation de la capacité de la zone de subduction des Petites Antilles à générer un tel évènement. Dans ce but, ces travaux tentent d'appréhender l'extension en profondeur du domaine sismogène de l'interplaque. Le manque de couverture des stations sismologiques permanentes est un inconvénient majeur dans l'exploration des Petites Antilles. Il s'explique en raison du peu de terres émergées et de leur éloignement de la zone potentiellement sismogène de l’interplaque. La région a donc fait l'objet de plusieurs campagnes océanographiques qui ont permis, notamment, le déploiement de sismomètres fond de mer (OBS); certains instruments étant restés immergés plusieurs mois afin de procéder à une écoute de la sismicité. Les travaux réalisés au cours de cette thèse se sont focalisés sur deux jeux de données de sismique grand-angle acquis au large des îles de la Dominique et de la Martinique. Leur analyse a permis la construction de modèles tomographiques 3D et 2D respectivement à l'échelle de l'avant-arc et de l'ensemble de la subduction. Ces modèles renseignent sur la structure sismique des plaques en convergence ainsi que sur leur géométrie. Ils permettent de discuter le rôle de la structure dans le fonctionnement de la subduction et d’obtenir une première estimation de l'extension en profondeur de la zone sismogène en considérant la portion de l'interplaque comprise entre le butoir et le Moho de la plaque supérieure. Dans un second temps, l'interprétation conjointe des modèles tomographiques et des localisations des séismes locaux a permis d'affiner cette estimation. Huit mois d'enregistrement de la sismicité montrent en effet une concentration des épicentres dans la région interne de l’avant-arc. Celle-ci présente un socle épais, à fort gradient de vitesse interprété comme plus dense et rigide que le socle de la région externe de l'avant-arc, plus déformable. La limite amont de la zone sismogène pourrait donc se situer en retrait de la position du butoir au contact de l'interplaque et de la limite entre ces deux zones de l'avant-arc. En profondeur, des mécanismes interplaques sont observés entre 35 et 45 km et interprétés comme des marqueurs de la limite aval de la zone sismogène. Cette dernière pourrait donc atteindre une profondeur jusqu'à 10 km supérieure à la limite précédemment évoquée et se trouver, par conséquent, au contact du manteau lithosphérique. L'ensemble de ces résultats suggèrent une extension en profondeur de la zone sismogène (i.e une largeur) de près de 70 km face aux îles de la Dominique et de la Martinique. Ces travaux doivent cependant être poursuivis afin d’évaluer pleinement la capacité de la zone de subduction des Petites Antilles à générer un éventuel séisme de méga-chevauchement. Le taux de couplage à l'interplaque doit être précisé ainsi que sa possible segmentation en lien, par exemple, avec à l'entrée en subduction des rides océaniques.

Droits : info:eu-repo/semantics/openAccess
http://archimer.ifremer.fr/doc/00172/28335/26628.pdf
http://archimer.ifremer.fr/doc/00172/28335/

Éditeur(s) : HAL CCSD Elsevier
Résumé : International audience
The southernmost part of the Ryukyu subduction, where the Philippine Sea Plate is subducting under the Eurasian Plate, is known to be a very seismically active region of transition from a north-dipping subduction along the Ryukyu subduction to an ~ SE-NW collision along the Taiwanese orogenic wedge. In this paper, we will focus on the Ryukyu forearc area close to Taiwan where the deformation is paroxysmal. In order to decipher the nature of the seismic deformation in this region, a three month passive experiment, combining 22 Ocean Bottom Seismometers and 51 onland stations, has been led. Starting from an a-priori heterogeneous model, we have obtained 801 well-located earthquake hypocenters, a precise P-wave tomography model and 14 focal mechanisms. The seismicity along the Ryukyu forearc is mainly located not only in the vicinity of the Interplate Seismogenic Zone (ISZ) but also within both the subducting PSP and the overriding plate. Seismicity within the upper-plate is essentially localized east of Nanao basin where NW-SE extension occurs, and northwest of the Hoping basin where strike-slip dominates. As revealed by both the P-wave velocity structure and the newly derived seismicity, we argue that a sub-vertical step offsetting the subducting PSP around 10 km may support the presence of a trench-parallel tear. The PSP also undergoes extension in its upper part that is probably caused by buckling and slab pull. The P-wave velocity structure reveals three other major features: (1) a continuity between the Central Range and the Ryukyu Arc with a shallower Moho (~ 30 km depth) between ~ 122.3°N and ~ 122.5°N along the Ryukyu Arc, (2) high P-wave velocities along the eastern side of the Central Range and, (3) two bodies with similar high crustal velocities (6.5-7.0 km/s) at 12-18 km depths, embedded within the Ryukyu arc basement, just north of Hoping Basin and north of the Nanao Basin.
ISSN: 0040-1951

hal-00795521
https://hal.archives-ouvertes.fr/hal-00795521
DOI : 10.1016/j.tecto.2012.04.011

Éditeur(s) : HAL CCSD Elsevier
Résumé : International audience
We investigated the deformation mechanisms and the seismic properties of 10 eclogite samples from different localities (Alps, Norway, Mali and eastern China) through the analysis of their microstructures and lattice-preferred orientations (LPO). These samples are representative of various types and intensity of deformation under eclogitic metamorphic conditions. Omphacite and garnet LPO were determined from electron backscatter diffraction (EBSD) technique. Garnet appears to be almost randomly oriented whereas omphacite develops strong LPO, characterized by the [001]-axes concentrated sub-parallel to the lineation, and the (010)-poles concentrated sub-perpendicular to the foliation. In order to analyze the deformation mechanisms that produced such omphacite LPO, we compare our observations to LPO simulated by viscoplastic self-consistent numerical models. A good fit to the measured LPO is obtained for models in which the dominant slip systems are 1/2h110i{11 ¯ 0}, [001] {110} and [001] (100). Dominant activation of these slip systems is in agreement with TEM studies of naturally deformed omphacite. Seismic properties of eclogite are calculated by combining the measured LPO and the single crystal elastic constants of omphacite and garnet. Although eclogite seismic anisotropies are very weak (less than 3% for both P-and S-wave), they are generally characterized by a maximum P-wave velocity sub-parallel to the lineation and by a minimum velocity approximately normal to foliation. The mean P-and S-wave velocities are high (respectively, 8.6 and 4.9 km/s). The S-wave anisotropy pattern displays complex relationships with the structural frame but the fast polarization plane generally tends to be parallel to the foliation. Calculated reflection coefficients show that an eclogite/crust interface is generally a good reflector (Rc > 0.1), whereas an eclogite body embedded in the upper mantle would be hardly detectable.
ISSN: 0040-1951

hal-01389715
http://hal.univ-reunion.fr/hal-01389715
http://hal.univ-reunion.fr/hal-01389715/document
http://hal.univ-reunion.fr/hal-01389715/file/Bascou_eclo-sismo_tphy2001_hal.pdf
DOI : 10.1016/S0040-1951(01)00156-1

Éditeur(s) : HAL CCSD American Geophysical Union
Résumé : International audience
Progress in seismic methodology and ambitious large-scale seismic projects are enabling high-resolution imaging of the continental crust. The ability to constrain interpretations of crustal seismic data is based on laboratory measurements on rock samples and calculations of seismic properties. Seismic velocity calculations and their directional dependence are based on the rock microfabric, which consists of mineral aggregate properties including crystallographic preferred orientation (CPO), grain shape and distribution, grain boundary distribution, and misorientation within grains. Single-mineral elastic constants and density are crucial for predicting seismic velocities, preferably at conditions that span the crust. However, high-temperature and high-pressure elastic constant data are not as common as those determined at standard temperature and pressure (STP; atmospheric conditions). Continental crust has a very diverse mineral composition; however, a select number of minerals appear to dominate seismic properties because of their high-volume fraction contribution. Calculations of microfabric-based seismic properties and anisotropy are performed with averaging methods that in their simplest form takes into account the CPO and modal mineral composition, and corresponding single crystal elastic constants. More complex methods can take into account other microstructural characteristics, including the grain shape and distribution of mineral grains and cracks and pores. Dynamic or active wave propagation schemes have recently been developed, which offer a complementary method to existing static averaging methods generally based on the use of the Christoffel equation. A challenge for the geophysics and rock physics communities is the separation of intrinsic factors affecting seismic anisotropy, due to properties of crystals within a rock and apparent sources due to extrinsic factors like cracks, fractures, and alteration. This is of particular importance when trying to deduce crustal composition and the state of deformation from seismic parameters.
ISSN: 8755-1209

Droits : info:eu-repo/semantics/OpenAccess
hal-01685568
https://hal.archives-ouvertes.fr/hal-01685568
https://hal.archives-ouvertes.fr/hal-01685568/document
https://hal.archives-ouvertes.fr/hal-01685568/file/Almqvist_et_al-2017-Reviewssur Geophysics.pdf
DOI : 10.1002/2016RG000552

Éditeur(s) : HAL CCSD AGU and the Geochemical Society
Résumé : International audience
The South Pacific region is characterized by a broadly elevated seafloor known as the South Pacific superswell. This region has a concentration of midplate volcanoes that experienced massive eruptions in the mid-Cretaceous period (90-120 Ma). These characteristics suggest the presence of a large-scale mantle plume beneath the South Pacific, called the South Pacific superplume. The geometry, origin depth, temperature, and composition of the superplume remain controversial, however, mainly due to the lack of seismological data that documents the mantle structure beneath the South Pacific. Seismic stations are sparse in the area due to its remote ocean environment. To obtain a better seismic image of the superplume, we deployed temporary broadband seismographs on oceanic islands and the seafloor in the South Pacific, which made possible the highest spatial resolution that has ever been achieved for the mantle structure beneath the region. The seismic image obtained from this new seismic data indicates that large-scale low-velocity anomalies (on the order of 1000 km in diameter), indicative of the superplume, are located from the bottom of the mantle to a depth of 1000 km, and small-scale low-velocity anomalies (on the order of 100 km in diameter) are present above it. A comparison of the seismic image with recent mantle convection studies based upon laboratory and numerical experiments suggests that the superplume may be a hot and chemically distinct mantle dome, and that the small-scale anomalies may be narrow plumes generated from the top of the dome. This model may explain various characteristics of hot spots in the South Pacific, such as the seafloor swell, short-lived hot spot chains, and the periodicity of massive eruptions.
ISSN: 1525-2027

hal-00443845
https://hal.archives-ouvertes.fr/hal-00443845
https://hal.archives-ouvertes.fr/hal-00443845/document
https://hal.archives-ouvertes.fr/hal-00443845/file/Suetsugu_PLUME_tomo_G32009_hal.pdf
DOI : 10.1029/2009GC002533

Éditeur(s) : HAL CCSD AGU and the Geochemical Society
Résumé : International audience
Although magmatism may occur during the earliest stages of continental rifting, its role in strain accommodation remains weakly constrained by largely 2-D studies. We analyze seismicity data from a 13 month, 39-station broadband seismic array to determine the role of magma intrusion on state-of-stress and strain localization, and their along-strike variations. Precise earthquake locations using cluster analyses and a new 3-D velocity model reveal lower crustal earthquakes beneath the central basins and along projections of steep border faults that degas CO2. Seismicity forms several disks interpreted as sills at 6-10 km below a monogenetic cone field. The sills overlie a lower crustal magma chamber that may feed eruptions at Oldoinyo Lengai volcano. After determining a new ML scaling relation, we determine a b-value of 0.87±0.03. Focal mechanisms for 65 earthquakes, and 13 from a catalogue prior to our array reveal an along-axis stress rotation of ∼60° in the magmatically active zone. New and prior mechanisms show predominantly normal slip along steep nodal planes, with extension directions ∼N90°E north and south of an active volcanic chain consistent with geodetic data, and ∼N150°E in the volcanic chain. The stress rotation facilitates strain transfer from border fault systems, the locus of early-stage deformation, to the zone of magma intrusion in the central rift. Our seismic, structural, and geochemistry results indicate that frequent lower crustal earthquakes are promoted by elevated pore pressures from volatile degassing along border faults, and hydraulic fracture around the margins of magma bodies. Results indicate that earthquakes are largely driven by stress state around inflating magma bodies.
ISSN: 1525-2027

hal-01667198
https://hal.archives-ouvertes.fr/hal-01667198
DOI : 10.1002/2017GC007027

Éditeur(s) : HAL CCSD Elsevier
Résumé : International audience
We use velocities from 65 continuous stations and 31 survey-mode GPS sites as well as kinematic modeling to investigate present day deformation along the Africa-Iberia plate boundary zone in the western Mediterranean region. The GPS velocity field shows southwestward motion of the central part of the Rif Mountains in northern Morocco with respect to Africa varying between 3.5 and 4.0 mm/yr, consistent with prior published results. Stations in the southwestern part of the Betic Mountains of southern Spain move west-southwest with respect to Eurasia (similar to 2-3 mm/yr). The western component of Betics motion is consistent with partial transfer of Nubia-Eurasia plate motion into the southern Betics. The southward component of Betics motion with respect to Iberia is kinematically consistent with south to southwest motion of the Rif Mountains with respect to Africa. We use block modeling, constrained by mapped surface faults and seismicity to estimate the geometry and rates of strain accumulation on plate boundary structures. Our preferred plate boundary geometry includes one block between Iberia and Africa including the SW Betics, Alboran Sea, and central Rif. This geometry provides a good fit to the observed motions, suggesting a wide transpressive boundary in the westernmost Mediterranean, with deformation mainly accommodated by the Gloria-Azores fault system to the West and the Rif-Tell lineament to the East. Block boundaries encompass aspects of earlier interpretations suggesting three main deformation styles: (i) extension along the NE-SW trending Trans-Alboran shear zone, (ii) dextral strike-slip in the Betics corresponding to a well defined E-W seismic lineament, and (iii) right lateral strike-slip motion extending West to the Azores and right-lateral motion with compression extending East along the Algerian Tell. We interpret differential motion in the Rif-Alboran-Betic system to be driven both by surface processes related the Africa-Eurasia oblique convergence and sub-crustal dynamic processes associated with the long history of subduction of the Neotethys ocean lithosphere. The dextral slip identified in the Betic Mountains in Southern Spain may be related to the offshore fault that produced the Great 1755 Lisbon Earthquake, and as such may represent a significant seismic hazard for the West Mediterranean region.
ISSN: 0012-821X

hal-00618337
https://hal.archives-ouvertes.fr/hal-00618337
DOI : 10.1016/j.epsl.2011.05.048