Éditeur(s) : HAL CCSD Elsevier
Résumé : International audience
We explore the use of seismic noise produced by rivers to monitor the bed load transport in the case of a low-discharge braided river in the French Alps: the "torrent de St Pierre". For this purpose, we deployed two dedicated seismic networks during summers 2007 and 2008, for which the characteristics of the recorded continuous signal are similar despite changes in the sensor locations. For dry weather conditions, only melting of nearby glaciers controls the supply of water to the stream. In these conditions, the river hydrology and the seismic energy in the 2-80 Hz frequency band both follow a diurnal fluctuation similar to the thermal amplitude. In contrast during rainfall episodes, the temperature variation fails to explain the hydrodynamic changes. Dense cloud covers reduce glacier melting and the recorded seismic energy denotes bursts of high-frequency seismic noise well correlated with water level data. Comparisons between the recorded seismic signals and the collected hydrological and sediment load data indicate that a frequency band of 3-9 Hz best explains the water level changes and thus the seismic waves coming from the flow turbulence. These analyses also reveal the presence of a seismic noise threshold that might be linked to the water shear stress exerted by the flowing water. Using the seismic energy in this frequency band as a proxy of the fluvial shear stress, the seismic-hydrologic relationship may be sensitive to variations in bed load transport. The spectral content of the seismic energy shows patterns consistent with the mobilization of sediment particles. From the interpretations of the seismic wave attenuation of river sources, we finally propose that stations at a distance from the stream less than 50 m are able to record most sediment particles. Farther stations are still useful during extreme events when largest grain sizes are mobilized. More generally this study demonstrates the feasibility of using the river seismic signal to survey bed load transport in various river types from small braided mountain rivers like the "torrent de St Pierre" to the large entrenched Himalayan rivers.
ISSN: 0022-1694

hal-00644789
https://hal.archives-ouvertes.fr/hal-00644789
DOI : 10.1016/j.jhydrol.2011.07.014

É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
B-type olivine fabrics are pervasive within highly depleted dunites of the small-sized Imono peridotite body located within the subduction-type Sanbagawa metamorphic belt of the southwest Japan arc. The dunites contain various microstructures, ranging from porphyroclastic to fine-grained intensely sheared textures. The Mg/(Mg+Fe) atomic ratios (Fo number) of olivine within these dunites are consistently around 0.9, as are the Cr/(Cr+Al) atomic ratios (Cr number) of chromian spinel, suggesting their evolution from a highly depleted magma (boninite). These data provide strong thermal constraints on the formation of the highly depleted dunites,as their formation requires hot, hydrous, shallow mantle (> 1250 degrees C at < 30 km depth) in the mantle wedge. Because the Sanbagawa metamorphic belt finally entrained these peridotites during progressive retrogression, B-type olivine fabrics probably developed in the fore-arc side of the subduction zone, above or along the subducting slab, possibly in association with dehydration fluids derived from the slab. The previously documented small magnitude of S-wave splitting can be explained by the seismic properties of B-type peridotites within an anisotropic layer of approximately several kilometers in thickness, oriented by flow parallel to the subducting slab, under maximum temperatures of 880-1030 degrees C depending on the flow stress. These findings indicate that such a B-type layer could constitute a dominant source of seismic S- and P-wave anisotropy in mantle wedge regions. (c) 2008 Elsevier B.V. All rights reserved.
ISSN: 0012-821X

hal-00412594
https://hal.archives-ouvertes.fr/hal-00412594
DOI : 10.1016/j.epsl.2008.06.014

É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
Résumé : East African Rift (EAR) is the divergent plate boundary. EAR exposes different stages of extension, from early stage rifting in Tanzania to oceanic accretion in Afar (Ethiopia). Manyara basin is the southernmost rift system of the east branch of EAR with recent volcanism (< 1.5 Ma) and a seismic swarm in the lower crust (20 – 40 km). Due to its location and tectonic setting, the Manyara basin offers the opportunity to study the earliest stage of rift initiation. Manyara volcanism is composed of several types of hyper-alkaline lavas as Mg-nephelinites (Mg# > 55) (Labait, Kwaraha), calciocarbonatite (Kwaraha) and evolved nephelinites (Mg# < 35) (Hanang).Mg-nephelinites (Labait and Kwaraha) are primary lavas mainly composed of olivine and clinopyroxene (cpx). Geochemical modelling from trace elements suggests that these primary magmas result from a degree of partial melting < 1 % from a CO2-garnet-phlogopite-bearing peridotite. These magmas have an asthenospheric source at depth > 120 km (lava carries xenoliths with equilibrium conditions > 4 GPa). The minerals were crystallized from a magma with a low H2O content (0.1 and 0.5 wt% H2O). The calciocarbonatite and evolved nephelinites are derived from Mg-nephelinites by fractional crystallization and immiscibility processes. Hanang nephelinites are silica- and alkaline-rich lavas (44.2 – 46.7 wt % SiO2, 9.5 –12.1 wt % Na2O+K2O, respectively) composed by cpx, Ti-garnet, nepheline, apatite and titanite. Complex zonation of cpx (e.g. abrupt change of Mg#, Nb/Ta, and H2O) and trace element patterns of nephelinites record magmatic differentiation involving open system with carbonate-silicate immiscibility and primary melt replenishment. The low H2O content of cpx (3 – 25 ppm wt. H2O) indicates that at least 0.3 wt % H2O was present at depth during carbonate-rich nephelinite crystallization at 340 – 640 MPa and 1050 – 1100 °C. The study of hosted-nepheline melt inclusions from Hanang allows constraining the late magmatic evolution of nephelinites during storage and magma ascent. Melt inclusions are composed by a silicate trachytic glass, a carbonate phase and a shrinkage bubble. Trachytic glass contains high content in CO2 (0.43 wt %, SIMS analyses), sulfur (0.21 – 0.92 wt % S), chlorine (0.28 –0.84 wt % Cl) and H2O low content (< 0.1 wt %, Raman analyses). Immiscibility process leading to the formation of carbonate occurs in a closed system during rapid magma ascent between 200 – 500 MPa. The carbonate phase is a Ca-Na-K-S-rich and anhydrous carbonate (33 wt % CaO, 20 wt % Na2O, 3 wt % K2O, and 3 wt % S). The pre-immiscible liquid has a phonolitic composition with 6 ± 1.5 wt % CO2 at 700 MPa. A preliminary study of melt inclusions by XANES spectroscopy and whole rocks by Mössbauer spectroscopy was used to determine these Manyara lavas were formed at oxidizing conditions (~ ΔFMQ +1.5).The early stage rifting volcanism (Manyara Basin) is characterized by CO2-rich and H2O-poor magmas from at least 120 km below the rift escarpment. The presence of CO2-rich magmas and the small amount of volcanic rocks erupted at the surface may indicate that the storage and percolation of these magmas at depth is a potential trigger for deep seismic swarms.
Le rift Est africain (REA) est une frontière de plaque en extension. Ce rift présente plusieurs stades d’extension, de l’initiation du rift en Tanzanie jusqu’à l’accrétion océanique en Afar. Le bassin de Manyara se situe le plus au sud de branche Est du REA. Il est caractérisé par la présence de volcanisme récent (< 1,5 Ma) et d’un essaim sismique dans la croûte inférieure (20 – 40 km). De par sa localisation et son contexte tectonique, le bassin de Manyara offre l’opportunité d’étudier le stade le plus précoce de l’initiation du rift. Le bassin de Manyara est composé de plusieurs types de laves hyperalcalines, les néphélinites magnésiennes (Mg# > 55) (Labait, Kwaraha), de calciocarbonatite (Kwaraha) et des néphélinites différenciées (Mg# < 35) (Hanang).Les néphélinites magnésiennes (Labait et Kwaraha) sont des laves primaires composées d’olivines et de clinopyroxènes (cpx). La modélisation géochimique des éléments en trace suggère que ces magmas primaires résultent d'un degré de fusion partielle ≤ 1 % à partir d'une péridotite à grenat et phlogopite. Ces magmas proviennent d’une profondeur > 120 km (présence de xénolites avec des conditions d’équilibre > 4 GPa). Les minéraux ont cristallisés à partir d’un magma pauvre en eau (0,1 et 0,5 pds % H2O). La calciocarbonatite et les néphélinites différenciés sont issues des néphélinites magnésiennes par cristallisation fractionnée et processus d’immiscibilité. Les néphélinites du Hanang sont riches en éléments alcalins (9,5 – 12,1 pds % Na2O+K2O) et en silice (44,2 – 46,7 pds% SiO2) et sont composés de cpx, grenat, néphéline, titanite et apatite. La zonation complexe dans les cpx (par exemple, changement brusque de Mg#, Nb/Ta, et H2O) implique une différenciation magmatique en système ouvert avec immiscibilité de liquide carbonaté et silicaté ainsi qu’un remplissage de la chambre magmatique avec des liquides primaires. La faible teneur en eau des cpx (3 – 25 ppm H2O) indique la présence d’un magma pauvre en eau (0,3 pds % H2O) lors de la cristallisation des cpx à des conditions crustales (340 – 640 MPa et 1050 – 1100 °C). L’étude des inclusions vitreuses dans les néphélines de Hanang permet de contraindre l'évolution magmatique tardive des néphélinites et le comportement des éléments volatils (CO2, H2O, S, F, Cl) lors du stockage et de la remontée du magma. Les inclusions vitreuses sont composées d’un verre trachytique, d’une phase carbonatée et d’une bulle de rétraction. Le verre trachytique contient du CO2 (0,43 pds % CO2, analyses SIMS), du soufre (0,21 à 0,92 pds% S), du chlore (0,28 – 0,84 pds % Cl) et très peu d’H2O (< 0,1 pds % H2O, analyses Raman). Le processus d’immiscibilité conduisant à la formation du carbonate se produit dans un système fermé pendant l'ascension rapide du magma, entre 200 – 500 MPa. La phase carbonatée est un carbonate anhydre et riche en Ca-Na-K-S (33 pds % CaO, 20 pds % Na2O, 3 pds % K2O, et 3 pds % S). Le liquide pré-immiscible a une composition phonolitique avec 6 ± 1,5 pds % CO2 à une pression de 700 MPa. Une étude préliminaire des inclusions par spectroscopie XANES et des roches par spectroscopie Mössbauer a permis de déterminer que les laves de Manyara se sont formées à conditions oxydantes (~ ∆FMQ +1,5).À l’initiation du rift, le volcanisme dans le bassin de Manyara est caractérisé par des magmas riches en CO2 et pauvres en H2O issus d’au moins 120 km de profondeur sous l'escarpement du rift. La présence de ces magmas riches en CO2 et la faible quantité de roches volcaniques émises à la surface peuvent indiquer que le piégeage et la percolation de ces magmas en profondeur est un déclencheur potentiel des essaims sismiques profonds.
https://tel.archives-ouvertes.fr/tel-01563231

NNT : 2016MONTT114
tel-01563231
https://tel.archives-ouvertes.fr/tel-01563231
https://tel.archives-ouvertes.fr/tel-01563231v2/document
https://tel.archives-ouvertes.fr/tel-01563231/file/2016_BAUDOUIN_archivage.pdf

Éditeur(s) : Gauthier-villars
Résumé : Seismic, Seabeam, magnetic and gravimetrie surveys in the southern part of the N-Fiji basin were conducted during the SEAPSO LEG III Cruise. The active spreading axis was recognized between 17°30'S and 21°S. It consists of a continuous N-S Ridge cross-cut by transverse directions oriented N25 and N45. North of 17°30'S, the spreading System is more complex due to the existence of a possible triple junction located around 15° S. The eastern part of the N-Fiji basin, west of Viti Levu is mostly characterized by alternating highs and deeps (<4,000 m). This area must be regarded no more as a spreading center, but as complex boundary suffering transverse tectonics.
La mission SEAPSO III a permis de localiser une zone d'accrétion active, d'orientation méridienne et pratiquement continue entre 17°30'S et 21° S. Des accidents transverses, qui encadrent ou recoupent cette dorsale, s'orientent autour de deux directions principales N25 et N45. Au Nord de 17 30'S le système d'accrétion, beaucoup plus complexe, annonce la triple jonction reconnue aux alentours de 15" S. La partie orientale du bassin Nord-Fidjien (à l'Ouest de l'île de Viti Levu) est caractérisée par l'alternance de zones hautes et de dépressions dont certaines dépassent 4000 m de profondeur. Il s'agit d'une zone complexe de déformations liées à des accidents cisaillants transverses.
Comptes Rendus de l'Academie des Sciences Serie II (1251-8069) (Gauthier-villars), 1986-06 , Vol. 303 , N. 1 , P. 93-98

Droits : Académie des Sciences
http://archimer.ifremer.fr/doc/00184/29492/27835.pdf
http://archimer.ifremer.fr/doc/00184/29492/

Éditeur(s) : HAL CCSD Bureau de recherches géologiques et minières
Résumé : National audience
Les tremblements de Terre sont sans doute l’un des phénomènes les plus catastrophiques associés à la dynamique de notre planète. Leur occurrence peut être associée d’une part au chargement lent et progressif des failles dû au mouvement des plaques tectoniques, mais également à la proximité d’autres séismes qui peuvent localement provoquer une déstabilisation du milieu. Les séismes sont classiquement interprétés dans le cadre de la théorie du cycle sismique (Fig. 1). La période asismique, correspondant au chargement séculaire des contraintes sur le plan des failles, est appelée période inter-sismique. Sa durée, entre 50 et 1000 ans, dépend du contexte géodynamique et notamment de la vitesse de chargement. La phase co-sismique correspond au séisme, i.e. à la phase de nucléation et de propagation de la rupture sismique. La période de temps associée est extrêmement courte et varie suivant les séismes entre quelques secondes à plusieurs minutes pour les méga-tremblements de terre comme celui de Sumatra en 2004 ou plus récemment celui du Japon en 2011. Enfin, suite à la rupture sismique, on considère une phase post-sismique qui correspond à un réajustement mécanique rapide du milieu qui peut s’effectuer sur plusieurs semaines, mois ou même années. La courte échelle de temps des observations sismologiques ou géodésiques (<100 ans) par rapport à la durée d’un cycle sismique (plusieurs centaines d’années) est un facteur limitant. Ainsi, aucun cycle sismique complet n’a pour l’instant pu être documenté.
ISSN: 0292-8477

insu-00913121
https://hal-insu.archives-ouvertes.fr/insu-00913121
https://hal-insu.archives-ouvertes.fr/insu-00913121/document
https://hal-insu.archives-ouvertes.fr/insu-00913121/file/geochronique_cattinetal.pdf

Éditeur(s) : HAL CCSD Seismological Society of America
Résumé : International audience
Abstract The continuous record of large surface-rupturing earthquakes along the Dead Sea fault brings unprecedented insights for paleoseismic and archaeoseismic research. In most recent studies, paleoseismic trenching documents the late Holocene faulting activity, while tectonic geomorphology addresses the long-term behavior (>10 ka), with a tendency to smooth the effect of individual earthquake rupture events (Mw >7). Here, we combine historical, archaeological, and paleoseismic investigations to build a consolidated catalog of destructive surface-rupturing earthquakes for the last 14 ka along the left-lateral Jordan Valley fault segment. The 120- km-long fault segment limited to the north and the south by major pull-apart basins (the Hula and the Dead Sea, respectively) is mapped in detail and shows five subsegments with narrow stepovers (width < 3 km). We conducted quantitative geomorphology along the fault, measured more than 20 offset drainages, excavated four trenches at two sites, and investigated archaeological sites with seismic damage in the Jordan Valley. Our results in paleoseismic trenching with 28 radiocarbon datings and the archaeoseismology at Tell Saydiyeh, supplemented with a rich historical seismic record, document 12 surface-rupturing events along the fault segment with a mean interval of ∼1160 yr and an average 5 mm=yr slip rate for the last 25 ka. The most complete part of the catalog indicates recurrence intervals that vary from 280 yr to 1500 yr, with a median value of 790 yr, and suggests an episodic behavior for the Jordan Valley fault. Our study allows a better constraint of the seismic cycle and related short-term variations (late Holocene) versus long-term behavior (Holocene and late Pleistocene) of a major continental transform fault.
ISSN: 0037-1106

hal-00747401
https://hal.archives-ouvertes.fr/hal-00747401
DOI : 10.1785/0120100097

É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 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 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) : HAL CCSD Elsevier
Résumé : International audience
This work focuses on the analysis of a unique set of seismological data recorded by two temporary networks of seismometers deployed onshore and offshore in the Central Lesser Antilles Island Arc from Martinique to Guadeloupe islands. During the whole recording period, extending from January to the end of August 2007, more than 1300 local seismic events were detected in this area. A subset of 769 earthquakes was located precisely by using HypoEllipse. We also computed focal mechanisms using P-wave polarities of the best azimuthally constrained earthquakes. We detected earthquakes beneath the Caribbean forearc and in the Atlantic oceanic plate as well. At depth seismicity delineates the Wadati-Benioff Zone down to 170 km depth. The main seismic activity is concentrated in the lower crust and in the mantle wedge, close to the island arc beneath an inner forearc domain in comparison to an outer forearc domain where little seismicity is observed. We propose that the difference of the seismicity beneath the inner and the outer forearc is related to a difference of crustal structure between the inner forearc interpreted as a dense, thick and rigid crustal block and the lighter and more flexible outer forearc. Seismicity is enhanced beneath the inner forearc because it likely increases the vertical stress applied to the subducting plate.
ISSN: 0040-1951

hal-01052679
https://hal.archives-ouvertes.fr/hal-01052679
DOI : 10.1016/j.tecto.2011.08.006

Éditeur(s) : HAL CCSD
Résumé : The aim of this thesis is to analyze the variations in slip rates along strike-slip faults, in intracontinental domain, producing large earthquakes (M > 7.5). To illustrate these variations, this work has focused in two area located in the most tectonically active continental domains in the world: in the western part of Mongolia (Bogd and Bulnay faults) and in the northern part of Iran (Astaneh and Tabriz faults). Using morphotectonic and paleoseismological analyses, the kinematics, the slip rates and the recurrence times have been estimated and allow us to describe the characteristics of seismic cycle along these faults. In Mongolia, the slip rates are estimated at ~ 1,2 and ~ 2,6 mm/yr along the Bogd and the Bulnay faults, respectively, with no variations of geological slip rates over the late Pleistocene-Holocene period. These two faults present characteristic slips and similar recurrence times between large earthquakes. The paleoseismological investigations suggest that a cluster occurred 3000 years ago, similar to the seismic cluster recorded in Mongolia during the XX century. In Iran, the slip rate was estimated to 2 mm/yr along the Astaneh fault and the recurrence times are ranging from 1600 to 2200 years, associated with offsets comprised between 3 and 4,5 m. We have also estimated a slip rate of ~ 7 mm/yr along the Tabriz fault, in agreement with the present day rate estimated by GPS, suggesting no variations in the slip rate over the past 45 ka.
Ce travail de thèse a pour but d'analyser les variations de vitesses sur des grandes failles décrochantes en contexte intracontinental, capables de produire des séismes de très fortes magnitudes (M > 7.5). Afin d'illustrer ces variations d'activités, cette analyse a été effectuée sur deux zones d'études situées en domaine continental et sismiquement actives: la région ouest de la Mongolie (failles de Bogd et Bolnay) et le nord de l'Iran (failles d'Astaneh et de Tabriz). À partir d'une approche morphotectonique et paléosismologique, les cinématiques, les vitesses de failles et les intervalles de récurrence entre les séismes majeurs ont été estimés, permettant d'analyser les caractéristiques du cycle sismique sur chacune des failles. En Mongolie, les failles de Bogd et Bolnay présentent respectivement des vitesses de ~ 1,2 et 2,6 mm/an, qui semblent être constantes sur la période Pleistocène supérieur-Holocène. Ces deux failles présentent également des glissements caractéristiques et des intervalles de temps similaires entre les séismes majeurs. Les analyses paléosismologiques suggèrent qu'un essaim sismique comparable à celui enregistré au XXème siècle a eu lieu il y a environ 3000 ans. En Iran, une vitesse géologique de 2 mm/an a été estimée sur la faille d'Astaneh et les données paléosismologiques suggèrent des intervalles de récurrence de 1800 ans, associés à des déplacements en surface compris entre 3 et 4,5 m. Nous avons également estimé une vitesse de 7 mm/an sur la faille de Tabriz, en accord avec les données GPS, suggérant que la vitesse sur cette faille est constante depuis 45 ka.
https://tel.archives-ouvertes.fr/tel-00591794

tel-00591794
https://tel.archives-ouvertes.fr/tel-00591794
https://tel.archives-ouvertes.fr/tel-00591794/document
https://tel.archives-ouvertes.fr/tel-00591794/file/these_rizza.pdf

Éditeur(s) : Amer Geophysical Union
Résumé : Seismic noise is an indirect source of information on ocean waves. Using a model of noise generation and propagation, seismic stations can be separated into those that are mostly sensitive to local sea states, and those that integrate sources from a large oceanic area. The model also provides a classification of noise-generating sea states into three classes. The analysis of Central California seismic noise data, well correlated with local waves, reveals that class I events dominate in summer, caused by a single wind-sea system, and for which ocean wave spectral levels are proportional to seismic spectral levels to an exponent b similar or equal to 0.9. In winter, noise is dominated by class II generation, for which coastal reflection is important, with a wave spectral density roughly proportional to the seismic spectral density to an exponent b similar or equal to 0.7. Sporadic events of class III probably produce some of the strongest noise events in Central California and need to be properly screened. These events are caused by opposed wave systems that are usually the wind-sea and a swell. This noise classification can be used to improve on the correlation between measured and estimated wave heights(up to r = 0.93 for daily averages). For other locations, where remote oceanic sources are recorded, a significant wave height estimated from the seismic noise compares well with area-averaged satellite data or wave model results(r > 0.85 for daily averages). These analyses pave the way for quantitative uses of seismic records, including the reconstruction of past wave climates, and the calibration of wave hindcasts.
Journal Of Geophysical Research-oceans (0148-0227) (Amer Geophysical Union), 2012-05 , Vol. 117 , N. C05002 , P. 19 pp.

Droits : 2012 AGU
http://archimer.ifremer.fr/doc/00083/19443/17051.pdf
DOI:10.1029/2011JC007449
http://archimer.ifremer.fr/doc/00083/19443/

É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
This study analyses the interactions between deformation and reactive fluid and melt percolation, and their effects on sub-arc mantle seismic properties based on microstructural observations on mantle xenoliths extracted by the Tubaf and Edison seamounts close to the Lihir Island, in the Papua New Guinea archipelago. These xenoliths sample an oceanic lithosphere, which has experienced high-temperature deformation in the presence of fluids or melts. This was followed by metasomatism under static conditions. Syn-kinematic percolation of reactive Si-rich melts or fluids in peridotites has produced pyroxene-enrichment, grain size reduction, and dispersion of olivine crystal preferred orientation (CPO). Fourier transform infrared spectroscopy analyses show that olivine has very low water contents (1-4 wt. ppm H2O), similar to spinel peridotites from other subduction zones. These low values may record both low water solubility in olivine at low pressure and dehydration during transport and exhumation. Water contents in pyroxenes are highly variable and likely result from spatially heterogeneous melt or fluid percolation. Analysis of olivine CPO indicates dominant activation of both (010)[100] and (001)[100] slip systems, which are characteristic of deformation under high temperature, low stress, low pressure and low to moderate hydrous conditions. Fast S-wave polarization and P- and Rayleigh propagation directions are thus parallel to the mantle flow direction. The pyroxene enrichment by melt-rock reactions is accompanied by dispersion of olivine CPO and induces a significant decrease of the maximum S-wave and P-wave anisotropy in the peridotites. The calculated seismic properties also show that the lowest Vp/Vs ratios (< 1.7) mapped in fore-arc mantle may only be explained by taking in consideration the CPO-induced elastic anisotropy of the peridotites.
ISSN: 0040-1951

hal-00950058
https://hal.archives-ouvertes.fr/hal-00950058
DOI : 10.1016/j.tecto.2013.09.024

Éditeur(s) : HAL CCSD Elsevier
Résumé : International audience
We explore the relations between deformation, annealing, and melt percolation during rifting and the effect of these processes on seismic anisotropy by analyzing the microstructures and crystal preferred orientations (CPO) in a suite of mantle xenoliths from Mega, in the southern end of the Ethiopian rift. Previous geochemical studies on these xenoliths showed evidence for interactions with variable melt types and volumes during the rifting process. The peridotites have dominantly coarse-porphyroclastic microstructures, but coarse granular or partially recrystallized microstructures also occur. The olivine CPO, characterized by orthorhombic to fiber-[100] patterns and moderate intensities, the common occurrence of (100) tilt walls, and the predominance of <0vw> rotation axes accommodating low angle misorientations in olivine support deformation by dislocation creep with dominant activation of the [100](010) system. Annealing (static recrystallization) of variable intensity followed this deformation. Modal enrichment in pyroxenes in > 60% of the studied peridotites corroborates extensive, but spatially heterogeneous reactive melt percolation leading to refertilization of the lithospheric mantle beneath the southern Ethiopian rift. The common interstitial shapes of the pyroxenes and lack of correlation between the pyroxenes and the olivine CPOs in many samples suggest that part of the refertilization is post-kinematic. However, there is no simple relation between reactive melt percolation and annealing of the olivine deformation microstructure. Comparison with data from other xenolith localities points to changes in the metasomatic imprint in the lithospheric mantle along the East African rift system correlated with the evolution in the rift maturity. Seismic properties averaged over all samples show typical lithospheric mantle patterns with fast propagation of P- and polarization of the fast S-waves parallel to the lineation. The anisotropy is moderate (< 6% for P-waves and < 4% for S-waves), implying that aligned melt must contribute to the strong seismic anisotropy measured in the Main Ethiopian rift, north of Mega.
ISSN: 0040-1951

hal-01355918
https://hal.archives-ouvertes.fr/hal-01355918
DOI : 10.1016/j.tecto.2016.05.027

É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
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) : 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/