É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) : Soc Geol France
Résumé : The Lesser Antilles subduction zone has produced no recent strong thrust earthquakes, making it difficult to quantify the seismic hazard from such events. The Lesser Antilles arc has a low subduction rate and an accretionary wedge that is very wide at its southern end. To investigate the effect of the wedge on seismogenesis, numerical models of forearc thermal structure were constructed along six transects perpendicular to the arc in order to determine the thermally predicted width of the seismogenic zone. The geometry of each section is constrained by published seismic profiles and crustal models derived from gravity and seismic data and by earthquake hypocenters at depth. A major constraint on the deep part of the model is that mantle temperature beneath the volcanic arc should achieve a temperature of 1,100 degrees C to generate partial melts. Predicted surface heat flow is compared to the available heat flow observations. Thermal modeling results indicate a systematic southward increase in the width of the seismogenic zone, more than doubling in width from north to south and corresponding to a dramatic southward increase in forearc width (distance from the arc to the deformation front of the accretionary wedge). The minimum width of the seismogenic zone (distance between the intersections of the subduction interface with the 150 degrees C and 350 degrees C isotherms) increases from about 80 km, north of 16 degrees N, to 230 km, at 13 degrees N. The maximum width (between the 100 degrees C and 450 degrees C isotherms) ranges from about 150 km in the north to up to 320 km in the south. This large variation in the width of the seismogenic zone is a consequence of the increasing width of the accretionary wedge to the south, caused by the increased thickness of sediment on the subducting plate. There is good agreement between the thermally predicted seismogenic limits and the sparse distribution of recorded thrust earthquakes, which are observed only in the northern portion of the arc. Possible scenarios for mega-thrust earthquakes are discussed. Depending on the segment length (along-strike) of the rupture plane, the occurrence of an event of magnitude 8-9 cannot be excluded.
L’absence de grands séismes récents à mécanismes chevauchants dans la zone de subduction des Petites Antilles rend difficile l’évaluation de l’aléa sismique lié à de tels événements. L’arc des Petites Antilles est caractérisé par une faible vitesse de subduction et par la présence d’un prisme d’accrétion très développé à son extrémité méridionale. Afin d’évaluer les effets de la largeur de ce prisme sur la genèse des séismes, nous avons étudié six sections perpendiculaires à l’arc, du nord au sud de celui-ci, pour déterminer la largeur de la zone sismogène prédite par les modèles thermiques appliqués à chacune de ces coupes. La géométrie de ces dernières est contrainte par les profils sismiques publiés, par les modèles de structure crustale déduits des données gravitaires et sismiques, et enfin par la distribution des hypocentres des séismes. Un contrôle important permettant de tester la validité des modèles thermiques en profondeur est qu’une température minimale de 1 100oC, compatible avec la fusion partielle du manteau hydraté, doit être atteinte sous l’arc volcanique actif. Par ailleurs, le flux thermique en surface prédit par ces modèles doit être compatible avec les mesures de flux de chaleur. Les modèles thermiques retenus d’après ces critères montrent une augmentation du simple au double vers le sud de la largeur de la zone sismogène, qui correspond à un élargissement considérable de la taille du domaine avant-arc. En effet, la largeur minimale de la zone sismogène (définie comme la distance entre les intersections de l’interface des plaques avec les isothermes 150o et 350oC) augmente d’environ 80 km au nord de 16oN jusqu’à 230 km à 13oN. La largeur maximale de cette zone (définie par les intersections de l’interface avec les isothermes 100o et 450oC) augmente, quant à elle, d’environ 150 km au nord jusqu’à 320 km au sud de l’arc. Cette variation considérable est la conséquence de l’augmentation de la largeur du prisme d’accrétion, elle-même causée par l’accumulation croissante des sédiments déposés sur la plaque plongeante. Les largeurs de la zone sismogène prédites à l’aide des modèles thermiques sont en bon accord avec les rares données disponibles sur les séismes à mécanismes chevauchants dans la partie nord de l’arc. Les scénarios possibles relatifs à des méga-séismes de ce type n’excluent pas de futurs événements atteignant des magnitudes de 8 à 9.
Bulletin De La Societe Geologique De France (0037-9409) (Soc Geol France), 2013 , Vol. 184 , N. 1-2 , P. 47-59

Droits : 2013 Societe Geologique de France
http://archimer.ifremer.fr/doc/00140/25140/29432.pdf
DOI:10.2113/gssgfbull.184.1-2.47
http://archimer.ifremer.fr/doc/00140/25140/

Éditeur(s) : HAL CCSD American Geophysical Union
Résumé : International audience
The brittle to ductile transition (BDT) in rocks may strongly influence their transport properties (i.e., permeability, porosity topology ...) and the maximum depth and temperature where hydrothermal fluids may circulate. To examine this transition in the context of Icelandic crust, we conducted deformation experiments on a glassy basalt (GB) and a glass-free basalt (GFB) under oceanic crust conditions. Mechanical and micro-structural observations at a constant strain rate of 10(-5) s(-1) and at confining pressure of 100-300 MPa indicate that the rocks are brittle and dilatant up to 700-800 degrees C. At higher temperatures and effective pressures the deformation mode becomes macroscopically ductile, i.e., deformation is distributed throughout the sample and no localized shear rupture plane develops. The presence of glass is a key component reducing the sample strength and lowering the pressure of the BDT. In the brittle field, strength is consistent with a Mohr-Coulomb failure criterion with an internal coefficient of friction of 0.42 for both samples. In the ductile field, strength is strain rate-and temperature-dependent and both samples were characterized by the same stress exponent in the range 3 < n < 4.2 but by very different activation energy Q(GB) = 59 +/- 15 KJ/mol and Q(GFB) = 456 +/- 4 KJ/mol. Extrapolation of these results to the Iceland oceanic crust conditions predicts a BDT at similar to 100 degrees C for a glassy basalt, whereas the BDT might occur in non-glassy basalts at deeper conditions, i.e., temperatures higher than 550 +/- 100 degrees C, in agreement with the Icelandic seismogenic zone.
ISSN: 0148-0227

hal-00760905
https://hal.archives-ouvertes.fr/hal-00760905
DOI : 10.1029/2011JB008884

Éditeur(s) : HAL CCSD Elsevier
Résumé : International audience
At the North-Ecuadorian convergent margin (1°S–1.5°N), the subduction of the rough Nazca oceanic plate leads to tectonic erosion of the upper plate and complex seismogenic behavior of the megathrust. We used three selected pre-stack depth migrated, multi-channel seismic reflection lines collected during the SISTEUR cruise to investigate the margin structure and decipher the impact of the subducted Atacames seamounts on tectonic erosion, interseismic coupling, and seismogenesis in the region of the 1942 Mw7.8 earthquake.This dataset highlights a subducted ∼30×40 km∼30×40 km, double-peak seamount that belongs to the Atacames seamount chain and that is associated with a deep morphologic re-entrant containing mass transport deposits.The seamount subduction uplifted the margin basement by ∼1.6 km and pervasively broke the margin by deep and intense reverse faulting ahead of the seamount, a process that is likely to weaken considerably the margin. In the seamount wake, the basement reverse fault system rotated counter-clockwise. This faulted basement is overlain with slope sediment sliding along listric normal faults that sole out onto the BSR. This superposition of deep tectonic contraction within the basement and shallow gravitational extension deformation within the sediment highlights the key role of gas hydrate on outer slope erosion.In addition to long-term regional basal erosion, the margin basement has thinned locally by an extra 0.8–1 km in response to the subduction of the Atacames seamount chain and hydrofracturing by overpressured fluids at the margin toe. This pervasively and deeply fractured margin segment is associated with a seismically quiet and GPS-modeled low interseismic coupling corridor that terminates downdip near the 1942 epicenter and locked zone. We suggest that the deeply buried double-peak Atacames seamount triggered the 1942 earthquake ahead of its leading flank. This result supports previous studies proposing that subducted seamounts provide unfavorable conditions for locking the updip segment of the plate boundary limiting the updip extent of seismogenic zones, but may favor large subduction earthquakes at greater depths.
ISSN: 0012-821X

hal-01256399
https://hal.archives-ouvertes.fr/hal-01256399
DOI : 10.1016/j.epsl.2015.10.043

Éditeur(s) : HAL CCSD
Résumé : National audience
The Lesser Antilles subduction zone has produced no recent strong thrust earthquakes, making it difficult to quantify the seismic hazard from such events. The Lesser Antilles arc has a low subduction rate and an accretionary wedge that is very wide at its southern end. To investigate the effect of the wedge on seismogenesis, numerical models of forearc thermal structure were constructed along six transects perpendicular to the arc in order to determine the thermally predicted width of the seismogenic zone. The geometry of each section is constrained by published seismic profiles and crustal models derived from gravity and seismic data and by earthquake hypocenters at depth. A major constraint on the deep part of the model is that mantle temperature beneath the volcanic arc should achieve a temperature of 1,100°C to generate partial melts. Predicted surface heat flow is compared to the available heat flow observations. Thermal modeling results indicate a systematic southward increase in the width of the seismogenic zone, more than doubling in width from north to south and corresponding to a dramatic southward increase in forearc width (distance from the arc to the deformation front of the accretionary wedge). The minimum width of the seismogenic zone (distance between the intersections of the subduction interface with the 150°C and 350°C isotherms) increases from about 80 km, north of 16°N, to 230 km, at 13°N. The maximum width (between the 100°C and 450°C isotherms) ranges from about 150 km in the north to up to 320 km in the south. This large variation in the width of the seismogenic zone is a consequence of the increasing width of the accretionary wedge to the south, caused by the increased thickness of sediment on the subducting plate. There is good agreement between the thermally predicted seismogenic limits and the sparse distribution of recorded thrust earthquakes, which are observed only in the northern portion of the arc. Possible scenarios for mega-thrust earthquakes are discussed. Depending on the segment length (along-strike) of the rupture plane, the occurrence of an event of magnitude 8-9 cannot be excluded.
Bulletin Société Géologique de France

insu-00814056
https://hal-insu.archives-ouvertes.fr/insu-00814056
DOI : 10.2113/gssgfbull.184.1-2.47

Éditeur(s) : HAL CCSD Springer
Résumé : Subduction is a major process that plays a first-order role in the dynamics of the Earth. The sinking of cold lithosphere into the mantle is thought by many authors to be the most important source of energy for plates driving forces. It also deeply modifies the thermal and chemical structure of the mantle, producing arc volcanism and is responsible for the release of most of the seismic energy on Earth. There has been considerable achievements done during the past decades regarding the complex interactions between the various processes acting in subduction zones. This volume contains a collection of contributions that were presented in June 2007 in Montpellier (France) during a conference that gave a state of the art panorama and discussed the perspectives about "Subduction Zone Geodynamics". The papers included in this special volume offer a unique multidisciplinary picture of the recent research on subduction zones geodynamics. They are organized into five main topics: Subduction zone geodynamics, Seismic tomography and anisotropy, Great subduction zone earthquakes, Seismogenic zone characterization, Continental and ridge subduction processes. Each of the 13 papers collected in the present volume is primarily concerned with one of these topics. However, it is important to highlight that papers always treat more than one topic so that all are related lighting on different aspects of the complex and fascinating subduction zones geodynamics.
https://hal.archives-ouvertes.fr/hal-00808864

ISBN : 978-3-540-87974
hal-00808864
https://hal.archives-ouvertes.fr/hal-00808864

Éditeur(s) : Elsevier Science Bv
Résumé : We present the results from a new grid of deep penetration multichannel seismic (MCS) profiles over the 280-km-long north-central segment of the Lesser Antilles subduction zone. The 14 dip-lines and 7 strike-lines image the topographical variations of (i) the subduction interplate décollement, (ii) the top of the arcward subducting Atlantic oceanic crust (TOC) under the huge accretionary wedge up to 7 km thick, and (iii) the trenchward dipping basement of the deeply buried forearc backstop of the Caribbean upper plate. The four northernmost long dip-lines of this new MCS grid reveal several-kilometres-high topographic variations of the TOC beneath the accretionary wedge offshore Guadeloupe and Antigua islands. They are located in the prolongation of those mapped on the Atlantic seafloor entering subduction, such as the Barracuda Ridge. This MCS grid also provides unexpected evidences on huge along-strike topographical variation of the backstop basement and of the deformation style affecting the outer forearc crust and sediments. Their mapping clearly indicates two principal areas of active deformation in the prolongation of the major Barracuda and Tiburon ridges and also other forearc basement highs that correspond to the prolongation of smaller oceanic basement highs recently mapped on the Atlantic seafloor. Although different in detail, the two main deforming forearc domains share similarities in style. The imaged deformation of the sedimentary stratification reveals a time- and space-dependent faulting by successive warping and unwarping, which deformation can be readily attributed to the forearc backstop sweeping over the two obliquely-oriented elongated and localized topographical ridges. The induced faulting producing vertical scarps in this transport does not require a regional arc-parallel extensional regime as proposed for the inner forearc domain, and may support a partitioned tectonic deformation such as in the case of an outer forearc sliver. A contrasted reflectivity of the sedimentary layering at the transition between the outer forearc and accretionary domains was resolved and used to define the seaward edge of the outer forearc basement interpreted as being possibly a proxy to the updip limit of the interplate seismogenic zone. Its mapping documents along-arc variations of some tens of kilometres the subduction backstop with respect to the negative gravity anomaly commonly taken as marking the subduction trench. With the exception of the southernmost part, the newly mapped updip limit reaches 25 km closer to the trench, thus indicating a possible wider seismogenic zone over almost the whole length of the study area.
Tectonophysics (0040-1951) (Elsevier Science Bv), 2013-09 , Vol. 603 , P. 32-54

Droits : 2013 Elsevier B.V. All rights reserved
http://archimer.ifremer.fr/doc/00139/25003/23109.pdf
DOI:10.1016/j.tecto.2013.05.028
http://archimer.ifremer.fr/doc/00139/25003/

Éditeur(s) : HAL CCSD
Résumé : International audience
Mylonite-bearing crustal shear zones are exhumed extensions at depth of (potentially seismogenic) brittle faultsat shallow structural level. By concentrating deformation, shear zones play a major role in the rheology and mechanicalbehavior of the continental crust. Analysis of shear zone geometry and microstructures is a useful toolto decipher deformation kinematics, though finite strain estimation is still relatively qualitative. Timing of thedeformation has also been the focus of numerous studies using the complete spectrum of geochronometers. However,only rare examples provide geochronological data that can be unambiguously linked with deformation. Thus,time-scales over which major mylonite zones develop and remain active under ductile conditions as well as theirstrain rates often remain poorly documented. Several sections across the South Armorican Shear Zone (SASZ),a crustal-scale, several km-thick dextral shear zone were investigated by structural, petrological and Ar/Ar radiochronologicalmethods. Finite strain profiles on these sections showed a strong deformation partitioning towardsthe ultramylonite core of the SASZ, with a decrease, of both the C’/S mean angle (from ca. 45 to less than5) and the Quartz grain-size (from ca. 150 to 5 m). In parallel, conventional and in-situ Ar/Ar datings on thefabric-forming white micas were performed on compositionally complex white-micas. Inherited magmatic muscovitescarried by the foliation and newly-formed, syn-tectonic substituted phengites located along the shear bandsyielded 2 age groups for the less deformed samples. Along the shear zone cross-section, the difference betweenthese two ages evolves from 10-15 Ma in weakly deformed domains to zero toward the SASZ core. In addition, thevery fine-grained syn-tectonic phengites yielded the same ca. 300-298 Ma age irrespectively of their distance to thecore of the shear zone, even in slightly deformed rocks distant by several km. These results support a new modelof crustal-scale shear zone activity, where the whole body of the shear zone is active throughout its history, i.e.the actively deforming zone does not get progressively narrower. The strain distribution nevertheless evolves withtime, with an increasingly larger fraction partitioned in the ultramylonitic core. Besides, combination of the finitestrain estimations and the Ar/Ar ages enables us to assess the range of strain rates associated with mylonitization.
EGU General Assembly 2014
Vienne, Austria

insu-01565057
https://hal-insu.archives-ouvertes.fr/insu-01565057

Éditeur(s) : HAL CCSD AGU and the Geochemical Society
Résumé : International audience
Based on global earthquake catalogs, the hypocenters, nodal planes, and seismic moments of worldwide subduction plate interface earthquakes were extracted for the period between 1900 and 2007. Assuming that the seismogenic zone coincides with the distribution of 5.5 <= M < 7 earthquakes, the subduction interface seismogenic zones were mapped for 80% of the trench systems and characterized with geometrical and mechanical parameters. Using this database, correlations were isolated between significant parameters to identify cause-effect relationships. Empirical laws obtained in previous studies were revisited in light of this more complete, accurate, and uniform description of the subduction interface seismogenic zone. The seismogenic zone was usually found to end in a fore-arc mantle, rather than at a Moho depth. The subduction velocity was the first-order controlling parameter for variations in the physical characteristics of plate interfaces, determining both the geometry and mechanical behavior. As such, the fast subduction zones and cold slabs were associated with large and steep plate interfaces, which, in turn, had large seismic rates. The subduction velocity could not account for the potential earthquake magnitude diversity that was observed along the trenches. Events with M-w >= 8.5 preferentially occurred in the vicinity of slab edges, where the upper plate was continental and the back-arc strain was neutral. This observation was interpreted in terms of compressive normal stresses along the plate interface. Large lateral ruptures should be promoted in neutral subduction zones due to moderate compressive stresses along the plate interface that allow the rupture to propagate laterally.
ISSN: 1525-2027

hal-00617681
https://hal.archives-ouvertes.fr/hal-00617681
https://hal.archives-ouvertes.fr/hal-00617681/document
https://hal.archives-ouvertes.fr/hal-00617681/file/heuretG32011.pdf
DOI : 10.1029/2010GC003230

Éditeur(s) : HAL CCSD Universidad Complutense de Madrid
Résumé : Based on new seismological and recent faulting data, we analyze five main sectors in the western part of the Betic Cordillera from an active tectonic perspective: i. The northwestern frontal sector of the cordillera, including the Morón de la Frontera area, with shallow (< 15 km depth) seismogenic active NE-SW reverse and strike-slip faults that do not reach the surface. ii. The Teba area, in the External Zones, which is deformed by sparse N-S to NW-SE normal faults with associated wedges of Quaternary sediments and active strike-slip faults at intermediate crustal levels. iii. The Ronda Basin, representing the major Late Miocene intramontane basin of the western Betic Cordillera, is mainly deformed by NNE-SSW and WNW-ESE folds and locally by NW-SE normal faults. iv. The internal part of the orogen, which is characterized by shallow (<40 km) and intermediate (40 to 120 km depth) seismicity beneath Málaga and the Alborán Sea, ENE-WSW folding, relief uplift, river incision and scarce active faults with surface expression. v. The Spanish Atlantic coast, which coincides with the western end of the Betic orogen, where Pliocene to Quaternary shallowmarine sedimentary rocks are deformed by meso-scale NW-SE oriented normal faults and few ENE-WSW reverse faults. These data support the recent tectonic activity of the western Betic Cordillera determined by the NW-SE Eurasian-African convergence with variable features from the mountain front, characterized by shallow deformations, up to the Internal Zones that mainly undergo uplift and intermediate seismicity.
ISSN: 1698-6180

hal-00750359
https://hal.archives-ouvertes.fr/hal-00750359
DOI : 10.5209/rev_JIGE.2012.v38.n1.39211

Éditeur(s) : HAL CCSD
Résumé : Fault zones are major discontinuities in sedimentary basins. Understanding their role on fluid migrations is an essential issue to (i) characterize the mechanisms and kinematics of deformation and (ii) to determine the parameters which control the distribution of energetic or mineral resources. This work applies to faulting under a temperature range of 200-250°C, corresponding to that of deeply buried reservoirs as well as potentially seismogenic fault zones. The studied faults are normal faults affecting the “Grès d’Annot” Formation, a Priabonian-Rupelian siliclastic turbidite succession of the Alpine foreland basin. They are located in two distinct areas: the Moutière-Restefond area in the eastern part of the basin and the Estrop area in the western part. Vitrinite reflectance outside fault zones indicates maximal temperatures of 240-260°C at Moutière-Restefond and 170-200°C at Estrop, i.e. burial depths around 8 km and 6 km, respectively, assuming a mean geothermal gradient of 30°C/km. Burial was due to underthrusting below the Embrunais-Ubaye Alpine nappes, the difference in burial between the two areas traducing the westward taper of the nappe front.The studied faults affect alternating arkosic sandstone beds and pelite layers with offsets from centimeters to decameters. Two types of core zones are recognized: (1) the foliated core zones of the Moutière-Restefond area, where deformation dominated by pressure solution and synkinematic phyllosilicate neoformation resulted in a foliated fabric, and (2) the core zones of the Estrop area, characterized by dilatant quartz-mineralized veins and breccia. Fault and host rock petrophysical properties were measured on drill plugs, using the water porosity technique and nitrogen permeability technique under 2 Mpa of confining pressure. In the Moutière-Restefond fault zones, plug axes were oriented in three orthogonal directions corresponding to the main deformation axes (X, Y, Z) of the foliated arkose to characterize the effect of structural anisotropy. In the Estrop fault zones, where deformation does not show equivalent anisotropy, the plug axes were oriented according to the main fault plane orientation (i.e., parallel to azimuth, parallel to dip and normal). In each fault zone, the porosity of fault rocks is equivalent to that of host rocks. In the Moutière-Restefond fault zones, the plugs parallel to the X and Z axes of deformation have a permeability similar to that of the host rocks, whereas permeability parallel to the Y axis (i.e. parallel to the foliation and the main fault plane) is one order of magnitude higher. By contrast, the permeability values in the Estrop fault zones are isotropic and equivalent to the host rock values. This study highlights the importance of various types of deformation mechanisms and related fluid-rock interactions in determining different fault rock types within a similar protolith, ultimately controlling distinct petrophysical properties of fault zones.
Les zones de failles constituent des discontinuités majeures dans la croûte supérieure. Ce mémoire s’intéresse aux failles actives à des profondeurs de 5-8 km, i.e. dans des conditions pouvant jouer un rôle important (i) sur le comportement mécanique de de la croûte, en particulier sur son potentiel sismogénique et (ii) sur la compartimentation des réservoirs géologiques profondément enfouis.Les failles étudiées sont des failles normales (rejet jusqu’à 50 mètres) affectant des alternances de bancs de grès arkosiques et de pélites dans les Grès d’Annot, une formation turbiditique d’âge priabonien-rupélien du bassin d’avant-chaine alpin. Les Grès d’Annot ont été enfouis sous les nappes de l’Embrunais-Ubaye peu de temps après leur dépôt et exhumés au Miocène moyen-supérieur. L’étude a été réalisée dans deux secteurs des Alpes de Haute-Provence, (1) le secteur de l’Estrop (chaine des Trois Evéchés, dans la partie occidentale du bassin) et (2) le secteur de Point Vert-Restefond (alentours du col de la Bonette, dans la partie orientale du bassin). La réflectance de la vitrinite dans les Grès d’Annot indique des températures maximales de 160-190°C dans le secteur (1) et de 240-260 °C dans le secteur (2) , soit des profondeurs d’enfouissement de 5-6 km et 8-10 km, respectivement, dans l’hypothèse d’un gradient géothermique moyen de 25-30°C/km. A partir des observations microstructurales, de données d’inclusions de fluides, de modélisations thermodynamiques sur chlorites syncinématiques et des indices de cristallinité de l’illite, on a établi que les failles du secteur (1) ont fonctionné à une température d’environ 200°C pour une pression de 0.6Kbar et celles du secteur (2) dans une gamme de température de 200-220°C et de pression comprise entre 0.6 et 1.1 Kba. L’activité de la faille de Point-Vert a été datée par la méthode 40Ar/39Ar au Miocène inférieur (21 Ma). Ces failles sont interprétées comme associées aux premiers stades du soulèvement résultant de l’activité des chevauchements dans le socle sous-jacent.Deux types principaux de roches de failles caractérisent les failles étudiées. Les brèches dilatantes sont majoritaires dans les zones de failles du secteur (1) alors que les arkoses foliées sont majoritaires dans les zones de failles du secteur (2). Les brèches dilatantes témoignent de mécanismes de déformation dominés par la fracturation transgranulaire et la pression-solution des silicates. Ce dernier mécanisme se manifeste principalement par d’abondantes précipitations cycliques de quartz, les stylolites restant peu développés. Les géométries structurales sont associées à l’activation de plans T de Riedel déterminant l’ouverture de relais extensifs et un régime localement dilatant. Nous associons ce type de roche de faille à une faible pression effective au cours de la déformation.La déformation cisaillante exprimée dans les arkoses foliées témoigne de mécanismes de déformation dominés par (i) la fracturation intragranulaire et transgranulaire, (ii) la pressionsolution/ précipitation des silicates, principalement le quartz, (iii) la néoformation de micas blancs liée à l’altération des feldspaths et (iv) la précipitation de chlorite. Ce travail montre également que les transformations minéralogiques liées aux interactions fluide-roche dans la zone de faille (forte augmentation de la proportion des phyllosilicates au dépend des feldspaths) changent radicalementla composition minéralogique de la roche. Nous associons ce type de roche de faille à une forte pression effective au cours de la déformation.La pétrophysique montre que les brèches de failles (secteur 1) ont une perméabilité isotrope et équivalente à celle du protolithe (10-3mD sous 2 MPa de confinement) tandis que les roches foliées (secteur 2) montrent une augmentation de perméabilité (10-2mD) selon l’axe Y de la déformation, d’un ordre de grandeur supérieure aux perméabilités suivant les axes X et Z et dans la roche hôte (10-3mD). Ce drain préférentiel potentiel est associé à l’orientation des phyllosilicates dans la foliation.
https://hal.archives-ouvertes.fr/tel-01310382

tel-01310382
https://hal.archives-ouvertes.fr/tel-01310382
https://hal.archives-ouvertes.fr/tel-01310382/document
https://hal.archives-ouvertes.fr/tel-01310382/file/Th%C3%A8seCAVAILHES2012.pdf

Éditeur(s) : HAL CCSD Oxford University Press (OUP)
Résumé : International audience
The Ryukyu subduction is known to generate very few earthquakes in its central segment contrarily to its two extremities. We focus in this study on the southernmost part of the Ryukyu subduction zone offshore east Taiwan. Our first task was to build a homogeneous earthquake catalogue for the period 1900–2007. The new catalogue provides homogenized M′W magnitudes and shows that several M′W≥ 7.0 earthquakes occurred offshore Hualien and Suao cities. We then focused on the 1920 June 5 earthquake (reported surface wave magnitude 8.1) previously located beneath the accretionary prism. The revised moment magnitude has been estimated in our catalogue at 7.7 ± 0.2. It is the biggest earthquake ever recorded in the Taiwan area but the fault that has produced this earthquake has not yet been identified with confidence. We relocated this event using regional phases (seismological bulletins archived at the Central Weather Bureau of Taiwan) about 50 km NNE and shallower of its former location, that is, within the Ryukyu Arc basement. According to earthquake bulletin information, revised magnitude, new hypocentral determination and known regional faults, we propose four potential active faults as candidates for the slip associated to this event: (1) the interplate seismogenic zone (ISZ), (2) an out-of-sequence thrust cutting through the forearc and branching on the ISZ, (3) a NS strike-slip fault cutting through the Ryukyu arc and (4) a N–S, westward dipping thrust fault, affecting the Philippine Sea Plate east of the Luzon Arc. The best compromise is to consider a rupture along the ISZ with a shallow nucleation possibly along a splay-fault followed by a downward and lateral propagation of the rupture that would explain the lack of significant seafloor motion and subsequent tsunami. We also estimate the maximum seismic coupling of the ISZ in the region east of Taiwan to about 0.4. In parallel, the evidences of aseismic slip occurring along the ISZ allow us to conclude that this region should only be affected by M < 8 earthquakes.
ISSN: 0956-540X

hal-00545459
https://hal.archives-ouvertes.fr/hal-00545459
https://hal.archives-ouvertes.fr/hal-00545459/document
https://hal.archives-ouvertes.fr/hal-00545459/file/theunisseGJI2010.pdf
DOI : 10.1111/j.1365-246X.2010.04813.x

É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 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
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
I use GPS interseismic velocities and classic 2D elastic half-space models with a screw dislocation to estimate the long-term fault slip rate, locking depth, and the offset between the surface fault trace and the location of the dislocation below the seismogenic zone for 13 segments along 8 major strike-slip faults. Using deduced strike-slip rates and the position of the dislocation to normalize the interseismic velocities to facilitate comparison of spatial patterns of deformation, I show that no substantial differences can be detected, ruling out a large asymmetry in interseismic velocities across the 8 faults used in this study. Only the Carrizo Plain segment of the San Andreas Fault shows a significant asymmetry that cannot be explained by shifting the position of the dislocation at depth relative to the fault trace. However, the resulting perturbation is less than 10% of total strike-slip rate. Fault traces are usually curved, defining a concave side. When the dislocation at depth is significantly offset from the fault trace, the shift is always toward the block on the concave side of the fault trace. This suggests that the fault zone in the lower crust may develop a simpler geometry more consistent with relative motion across the fault than its upper seismogenic part constrained by the structural complexity of the brittle crust. Since the faults used in this study are at different times in their interseismic period, comparing the interseismic velocity fields across them allows identification of possible variations of the interseismic velocities with time. When normalized by slip rate and dislocation location, all the faults show the same interseismic strain with no significant differences between deduced locking depths. These comparisons suggest that if temporal variations occur as suggested by some dynamic earthquake cycle models, they are small and below the accuracy of the available geodetic measurements.
ISSN: 0040-1951

hal-01171753
https://hal.archives-ouvertes.fr/hal-01171753
DOI : 10.1016/j.tecto.2015.01.013

Éditeur(s) : HAL CCSD Nature Publishing Group
Résumé : International audience
Temperature and fluid pressure conditions control rock deformation and mineralization on geological faults, and hence the distribution of earthquakes1. Typical intraplate continental crust has hydrostatic fluid pressure and a near-surface thermal gradient of 31 ± 15 degrees Celsius per kilometre2, 3. At temperatures above 300–450 degrees Celsius, usually found at depths greater than 10–15 kilometres, the intra-crystalline plasticity of quartz and feldspar relieves stress by aseismic creep and earthquakes are infrequent. Hydrothermal conditions control the stability of mineral phases and hence frictional–mechanical processes associated with earthquake rupture cycles, but there are few temperature and fluid pressure data from active plate-bounding faults. Here we report results from a borehole drilled into the upper part of the Alpine Fault, which is late in its cycle of stress accumulation and expected to rupture in a magnitude 8 earthquake in the coming decades4, 5. The borehole (depth 893 metres) revealed a pore fluid pressure gradient exceeding 9 ± 1 per cent above hydrostatic levels and an average geothermal gradient of 125 ± 55 degrees Celsius per kilometre within the hanging wall of the fault. These extreme hydrothermal conditions result from rapid fault movement, which transports rock and heat from depth, and topographically driven fluid movement that concentrates heat into valleys. Shear heating may occur within the fault but is not required to explain our observations. Our data and models show that highly anomalous fluid pressure and temperature gradients in the upper part of the seismogenic zone can be created by positive feedbacks between processes of fault slip, rock fracturing and alteration, and landscape development at plate-bounding faults.
ISSN: 0028-0836

hal-01685757
https://hal.archives-ouvertes.fr/hal-01685757
DOI : 10.1038/nature22355

Éditeur(s) : HAL CCSD American Geophysical Union
Résumé : International audience
Phyllosilicate content and related permeability of fault zones form primary controls on hydraulic and mechanical behavior of the brittle crust. Hence, understanding and predicting the localization of these ubiquitous minerals is a major issue for fundamental and applied geosciences. We describe normal fault zones cutting a foreland arkosic turbiditic formation suffering high-T diagenesis and formed under conditions (~200°C) typical of deeply buried reservoirs and common within the seismogenic interval. Microstructural analyses show a large proportion of phyllosilicates (up to 34%) in the fault rock, derived from near-complete feldspar alteration and disaggregation during deformation. This study shows that even faults with offsets (~20 cm) much lower than bed thickness can have such large feldspar-to-phyllosilicate transformation ratios, implying that the origin of the phyllosilicates is purely transformation related. These results imply that the potential sealing capacity and strength of faults could be predicted from the host rock feldspar content. Where sealing capacity and fault strength can be related to phyllosilicate content, these properties can then also be inferred from the predicted phyllosilicate content: this opens up new horizons concerning the hydraulic and the mechanical behavior of the upper crust.
ISSN: 0094-8276

hal-00903574
https://hal.archives-ouvertes.fr/hal-00903574
DOI : 10.1002/grl.50829

Éditeur(s) : HAL CCSD Elsevier
Résumé : International audience
Paleostress orientations were calculated from fault populations at 24 sites along the SW-NE segment and five sites along the E-W. Yasin segment of the Karakoram-Kohistan Suture Zone in NW Pakistan They demonstrate the importance of combined thrusting and strike-slip faulting. However, several paleostress tensor directions are distinguished, a dominant NW-SE compression and a minor E-W compression are compatible with the recent evolution of this part of the Hindu Kush. From the lack of both systematic overprinting-relationships and spatial trend (the two tensors were obtained at different locations) we conclude that in each location any of these two shortening directions can dominate Heterogeneously distributed extension is found in some places and is likely due to local conditions These paleostress tensors substantiate a transpressional regime due to far-field Himalayan compression and document the long-term background of the seismogenic deformation in this region.
ISSN: 0743-9547

hal-00496459
https://hal.archives-ouvertes.fr/hal-00496459
DOI : 10.1016/j.jseaes.2010.01.004

Éditeur(s) : HAL CCSD Elsevier
Résumé : International audience
We use 3-D mechanical modeling representing faults as planar surfaces with frictional properties that obey Coulomb-failure process to explore the long-term slip rates of the Altyn Tagh fault and Kunlun faults in the north Tibetan plateau. Crustal theology is simplified as an elastoplastic upper crust and a viscoelastic lower crust. Far-field GPS velocities and late Quaternary fault slip rates are used to constrain the model results. Rheological tests show that effective fault friction lower than 0.1-0.08 leads to high slip rates that fit with geologically and geodetically determined slip rates of the Kunlun fault (10-11.7 +/- 1.5 mm/yr). Meanwhile, the modeled Altyn Tagh fault reaches slip rates similar to 13.7 mm/yr to similar to 7.8 mm/yr in its central portion, between ranges of the geological slip rates. Associated with high slip rates, our model predicts that central Tibet (similar to 84 degrees E-95 degrees E) from the Altyn Tagh fault to the north of the Himalayan arc accommodates north-south shortening and east-west extension rates of similar to 10-12 mm/yr and similar to 8-10 mm/yr, respectively. We also question the widely accepted idea that interseismic strain is driven at the base of the seismogenic zone by a screw dislocation. If this assumption fails, the presented model implies that interseismic strain around large strike-slip faults could be distributed in a much broader way if the lithosphere deforms as a thin elastic plate rather than an elastic half-space with an embedded dislocation. If this distributed deformation is ignored, and the instantaneous surface deformation field modeled as that resulting from slip on a dislocation below a specified depth embedded in an elastic half-space, the estimated slip rate will inevitably be lower than the true long-term slip rate. This appears to explain why geodetic slip rates proposed for the Altyn Tagh fault (510 mm/yr) are lower than some of the geological slip rates. (C) 2008 Elsevier B.V. All rights reserved.
ISSN: 0012-821X

hal-00411926
https://hal.archives-ouvertes.fr/hal-00411926
DOI : 10.1016/j.epsl.2008.06.049