É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 Elsevier
Résumé : Primary arc melts may form through fluxed or adiabatic decompression melting in the mantle wedge, or via a combination of both processes. Major limitations to our understanding of the formation of primary arc melts stem from the fact that most arc lavas are aggregated blends of individual magma batches, further modified by differentiation processes in the sub-arc mantle lithosphere and overlying crust. Primary melt generation is thus masked by these types of second-stage processes. Magma-hosted peridotites sampled as xenoliths in subduction zone magmas are possible remnants of sub-arc mantle and magma generation processes, but are rarely sampled in active arcs. Published studies have emphasised the predominantly harzburgitic lithologies with particularly high modal orthopyroxene in these xenoliths; the former characteristic reflects the refractory nature of these materials consequent to extensive melt depletion of a lherzolitic protolith whereas the latter feature requires additional explanation.Here we present major and minor element data for pristine, mantle-derived, lava-hosted spinel-bearing harzburgite and dunite xenoliths and associated primitive melts from the active Kamchatka and Bismarck arcs. We show that these peridotite suites, and other mantle xenoliths sampled in circum-Pacific arcs, are a distinctive peridotite type not found in other tectonic settings, and are melting residues from hydrous melting of silica-enriched mantle sources. We explore the ability of experimental studies allied with mantle melting parameterisations (pMELTS, Petrolog3) to reproduce the compositions of these arc peridotites, and present a protolith (‘hybrid mantle wedge’) composition that satisfies the available constraints. The composition of peridotite xenoliths recovered from erupted arc magmas plausibly requires their formation initially via interaction of slab-derived components with refractory mantle prior to or during the formation of primary arc melts. The liquid compositions extracted from these hybrid sources are higher in normative quartz and hypersthene (i.e., they have a more silica-saturated character) in comparison with basalts derived from prior melt-depleted asthenospheric mantle beneath ridges. These primary arc melts range from silica-rich picrite to boninite and high-Mg basaltic andesite along a residual spinel harzburgite cotectic. Silica enrichment in the mantle sources of arc-related, subalkaline picrite-boninite-andesite suites coupled with the amount of water and depth of melting, are important for the formation of medium-Fe (‘calc-alkaline’) andesite-dacite-rhyolite suites, key lithologies forming the continental crust.
ISSN: 0016-7037

hal-01553520
https://hal.archives-ouvertes.fr/hal-01553520
DOI : 10.1016/j.gca.2016.09.030

É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 Elsevier
Résumé : International audience
To better understand the microstructural evolution of a "serpentinized" mantle rock and the influence of various hydrous phases on the seismic properties of the mantle wedge, we have conducted the detailed microstructural analyses of a sample of tremolite-chlorite-antigorite schist collected from the Moses Rock dike (central part of the Colorado Plateau). We performed differential effective media (DEM) modelling to study the effect of three hydrous phases forming two-phase aggregates with olivine, considering the crystallographic preferred orientation (CPO) of each phase and the shape ratio of the hydrous phases. We have demonstrated that in a partially serpentinized peridotite, the olivine CPO characteristic of [100](010) dislocation glide is still preserved, and the high-temperature asthenospheric flow is preserved with a foliation normal to that of antigorite schist. The transformation of olivine into antigorite occurs predominantly (~ 75%) by the relationship (100)ol || (001)atg with [001]ol || [010]atg, with the (010)ol || (001)atg and [001]ol || [010]atg relationship observed in areas of weak antigorite CPO. Chlorite results from the phase transformation of olivine in a relatively static environment, as shown by the correlation between the olivine-chlorite CPOs with (100)ol || (100)ch, (010)ol || (001)ch and (001)ol || (010)ch. The fluid percolation that caused the localized metasomatism and partial hydration of the mantle occurred possibly along trans-lithospheric shear zones. The presence of chlorite induces the most important drop on the P-wave velocities and may help to explain some local low velocities in the fore-arc mantle wedges, but is unlikely to be of global importance due to its very high Vp/Vs ratio ~ 1.9. On the other hand, antigorite is the only phase that causes important modification on the propagation directions of P and S-waves, and the only phase to explain the polarization of the fastest shear waves parallel to the subduction trench.
ISSN: 0040-1951

hal-00858061
https://hal.archives-ouvertes.fr/hal-00858061
DOI : 10.1016/j.tecto.2013.03.022

Éditeur(s) : HAL CCSD Oxford University Press (OUP)
Résumé : International audience
The mantle above a subducting slab is the site of complex interactions between deformation, partial melting, fluid migration and magma transport. To constrain these interactions and their effects on olivine deformation, we analyze microstructures, crystal preferred orientations, and water contents of peridotite xenoliths entrained by andesites of the Avacha volcano, southern Kamchatka arc. These xenoliths are refractory spinel harzburgites that have coarse-grained microstructures with widely spaced subgrain boundaries and sinuous grain boundaries in olivine, consistent with deformation by dislocation creep under low deviatoric stress (< 13 MPa) and with a significant contribution from diffusional processes. Analysis of crystal preferred orientations (CPO) indicates dominant activation of high-temperature, low-stress {0kl}[100] slip systems in olivine and of(100)[001] in orthopyroxene. In most samples, coarse opx crystals, elongated parallel to the lineation, enclose small olivine grains in crystallographic continuity with neighbouring crystals, indicating crystallization of orthopyroxene at the expense of olivine as a result of reactive percolation of Si-rich fluids coeval with the high-temperature deformation. Secondary crystallization of interstitial orthopyroxene led locally to development of opx-rich lenses parallel to the foliation, characterized by a decrease in olivine grain size and dispersion of the olivine CPO without changing the dominant slip systems. Half of the samples also show acicular orthopyroxene aggregates with which is associated a fine-grained matrix composed of rounded strain-free olivine, orthopyroxene, spinel, and rare amphibole crystals. This matrix occurs pervasively along grain boundaries or forms millimeter-scale irregular lenses and anastomosing veinlets that crosscut the coarse crystals and their ductile deformation structures. Both acicular orthopyroxene and the fine-grained matrix are interpreted as resulting from reactive transport of H2O-rich fluids under static conditions, probably in the lithospheric mantle. Infrared analyses show that olivine contains 1-8 center dot 6 ppm by weight of water. These low water contents are similar to those observed in spinel peridotites from other subduction zones and probably record both the low solubility of water in olivine at low pressure and dehydration during exhumation of the xenoliths. Water contents in orthopyroxene are highly variable (25-506 ppm H2O), probably recording spatially heterogeneous interaction with fluids or melts and compositional disequilibrium in the studied samples. Change in the dominant percolation mechanism from porous flow to fracturing suggests cooling, consistent with the low temperatures estimated from pyroxene thermometry (< 800-900 degrees C). The Avacha xenoliths therefore record pervasive deformation of a region of the mantle under asthenospheric conditions, followed by its accretion to the base of the lithosphere, probably as a result of cooling of the mantle wedge. Percolation of Si-rich fluids or hydrous melts is recorded at all stages; this probably enhanced diffusion and lowered deviatoric stresses during ductile deformation, but did not change the dominant slip direction in olivine from [100] to [001].
ISSN: 0022-3530

hal-00456053
https://hal.archives-ouvertes.fr/hal-00456053
DOI : 10.1093/petrology/egp085

É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
Résumé : This work provides new constraints on the interactions between deformation and percolation of melt or fluids, as well as on the implications of these processes for the seismic properties of the mantle wedge. It is based on the analysis of spinel peridotites from the massif of Ronda (Spain) and two xenolith suites from active subduction zones (Kamchatka, Papua New Guinea). A coupled structural, petrological and geochemical study of these samples shows that they underwent a reactive percolation of melts or Si-rich fluids, which was synchronous to a deformation event occurring under high temperature and low stress conditions, consistent with the PT conditions at the base of the lithosphere or in the asthenosphere. This reactive percolation is responsible for enrichment in pyroxenes, mainly orthopyroxene, which is often localized in bands parallel to the high-temperature foliation. This enrichment is associated with grain size reduction and dispersion of the crystallographic orientation of olivine. The dominant slip system in olivine is {0kl}[100], which results in fast S-wave polarization parallel to the flow direction in the mantle. The decrease in the intensity of the olivine crystal preferred orientations associated with the enrichment in pyroxene results in significant decrease of the anisotropy that may induce an underestimation of the thickness of the anisotropic layer by up to 33%. The observed orthopyroxene enrichment also lowers the Vp/Vs ratio, but cannot explain Vp/Vs < 1.7 mapped locally in the fore-arc mantle in Japan and the Andes. Such low Vp/Vs ratios may however be explained by considering the intrinsic anisotropy of the peridotites, which is generally ignored in large-scale Vp/Vs ratio mapping of the mantle wedge. Infrared analyses show that olivine from the both xenolith collections contains less water than the theoretical saturation calculated for their estimated equilibrium temperature in the spinel stability field. These low water content are similar to those observed in spinel peridotites from other subduction zones and probably record both the low solubility of water in olivine at relatively low pressure and dehydration during exhumation of the xenoliths. These measured water contents as well as theoretical saturation estimations are not sufficient to change the dominant slip direction of in olivine from [100] to [001].
Ce travail apporte de nouvelles contraintes sur les interactions entre déformation et processus d'hydratation et de percolation de magma ou de fluides dans le manteau à l'aplomb d'une zone de subduction et leurs implications sur les propriétés sismiques dans le coin mantellique. Il se base sur l'analyse de péridotites à spinelle provenant du massif de Ronda (Espagne) et deux séries de xénolites issues de zones de subduction actives (Kamchatka, Papouasie-Nouvelle-Guinée). L'étude structurale, pétrologique et géochimique de ces échantillons montrent qu'ils ont subit une percolation réactive de magma ou de fluide synchrone d'une déformation de haute température basse contrainte cohérente avec les condition PT de la base de la lithosphère ou de l'asthénosphère. Cette percolation réactive est responsable d'un enrichissement en pyroxènes, qui est souvent localisé dans des bandes parallèles à la foliation. Cet enrichissement est associé à une dispersion de l'orientation cristallographique de l'olivine. Les systèmes de glissement dominant dans l'olivine sont dans tous les cas {0kl}[100], ce qui implique que la direction de polarisation rapide des ondes S dans la partie supérieure du coin mantellique est parallèle à la direction d'écoulement dans le manteau. La décroissance de l'intensité des OPR de l'olivine associée à l'enrichissement en pyroxènes a pour conséquence une réduction non négligeable de l'anisotropie sismique qui peut induire jusqu'à 33% d'erreur sur l'interprétation de l'épaisseur de la couche anisotrope. Un enrichissement en orthopyroxène peut aussi entraîner une baisse du rapport Vp/Vs, mais ne peut expliquer les rapports Vp/Vs <1,7 cartés dans certains avant-arcs. Cependant de telles valeurs peuvent être expliquées si l'anisotropie des péridotites du coin mantellique est prise en compte. Les analyses par spectroscopie infra-rouge montrent que les olivines des deux séries de xénolites étudiées contiennent moins d'eau que la saturation théorique calculée pour leurs températures d'équilibre dans le domaine du spinelle. Ces faibles concentrations en eau sont similaires à celles observées dans les olivines des péridotites à spinelle des autres zones de subduction. Elles enregistrent probablement la faible solubilité de l'eau dans l'olivine à des pressions relativement basses et la déshydratation au cours de l'exhumation des xénolites. Ces teneurs en eau mesurées ainsi que les estimations de la saturation théorique ne sont pas suffisantes pour changer la direction de glissement dominante dans l'olivine de [100] à [001].
https://tel.archives-ouvertes.fr/tel-00689832

tel-00689832
https://tel.archives-ouvertes.fr/tel-00689832
https://tel.archives-ouvertes.fr/tel-00689832/document
https://tel.archives-ouvertes.fr/tel-00689832/file/These-Soustelle2010.pdf

Éditeur(s) : HAL CCSD Elsevier
Résumé : International audience
The Tso Morari serpentinites in the Ladakh area, northwest Himalaya, originated from the forearc mantle overlying the northward subducting Neo-Tethys lithosphere and the margin of the Indian continent. The serpentinites are characterized by high concentration of fluid-mobile elements (FME: As, Sb, B, Li, and U) compared to ophiolitic or abyssal serpentinites. The Pb isotopic compositions of serpentinites show influence of the subducted Indian continental lithosphere. Trace element concentrations of antigorite determined in situ with Laser Ablation High Resolution Inductively Coupled Mass Spectrometer (LA-HR-ICP-MS) show high contents of FME including Pb, in contrast to the spatially associated iron oxides. Rare earth elements (REE) and compatible elements, such as Sc and Co, remained immobile during the hydration, allowing the identification of the primary minerals (olivine or orthopyroxene) from which serpentine formed. Serpentinized olivine displays higher Sb and As concentrations (up to 1000 x PM) than serpentinized orthopyroxenes that are enriched in Pb, Cs and Li (2 to up to 10 x PM).;We propose that the observed FME distribution in two types of serpentine reflect the differential incorporation of FME during the downward movement of the serpentinite along the subduction plane. At temperature lower than 400 degrees C, at shallow depths, olivine is preferentially serpentinized and incorporates elements that are fluid soluble at low temperatures, such as Sb and As. Above 400 degrees C, orthopyroxene is hydrated and incorporates Pb, Cs, Li and possibly Ba. Boron and U are incorporated in both types of serpentine suggesting that they are released from slabs at temperatures around 300-400 degrees C. The serpentine acts as a sink for water, but also for FME and transports them to deeper and hotter levels in the mantle, down to the isotherm 600-650 degrees C where dehydration occurs.
ISSN: 0009-2541

hal-00475565
https://hal.archives-ouvertes.fr/hal-00475565
DOI : 10.1016/j.chemgeo.2009.10.002

Éditeur(s) : HAL CCSD European Geosciences Union
Résumé : The properties of the subduction interplate domain are likely to affect not only the seismogenic potential of the subduction area but also the overall subduction process, as it influences its viability. Numerical simulations are performed to model the long-term equilibrium state of the subduction interplate when the diving lithosphere interacts with both the overriding plate and the surrounding convective mantle. The thermomechanical model combines a non-Newtonian viscous rheology and a pseudo-brittle rheology. Rock strength here depends on depth, temperature and stress, for both oceanic crust and mantle rocks. I study the evolution through time of, on one hand, the brittle-ductile transition (BDT) depth, zBDT, and, on the other hand, of the kinematic decoupling depth, zdec, simulated along the subduction interplate. The results show that both a high friction and a low ductile strength at the asthenospheric wedge tip shallow zBDT. The influence of the weak material activation energy is of second order but not negligible. zBDT becomes dependent on the ductile strength increase with depth (activation volume) if the BDT occurs at the interplate decoupling depth. Regarding the interplate decoupling depth, it is shallowed (1) significantly if mantle viscosity at asthenospheric wedge tip is low, (2) if the difference in mantle and interplate activation energy is weak, and (3) if the activation volume is increased. Very low friction coefficients and/or low asthenospheric viscosities promote zBDT = zdec. I then present how the subducting lithosphere age affects the brittle-ductile transition depth and the kinematic decoupling depth in this model. Simulations show that a rheological model in which the respective activation energies of mantle and interplate material are too close hinders the mechanical decoupling at the down-dip extent of the interplate, and eventually jams the subduction process during incipient subduction of a young (20-Myr-old) and soft lithosphere under a thick upper plate. Finally, both the BDT depth and the decoupling depth are a function of the subducting plate age, but are not influenced in the same fashion: cool and old subducting plates deepen the BDT but shallow the interplate decoupling depth. Even if BDT and kinematic decoupling are intrinsically related to different mechanisms of deformation, this work shows that they are able to interact closely. Comparison between modelling results and observations suggests a minimum friction coefficient of 0.045 for the interplate plane, even 0.069 in some cases, to model realistic BDT depths. The modelled zdec is a bit deeper than suggested by geophysical observations. Eventually, the better way to improve the adjustment to observations may rely on a moderate to strong asthenosphere viscosity reduction in the metasomatised mantle wedge.
ISSN: 1869-9510

hal-00795559
https://hal.archives-ouvertes.fr/hal-00795559
DOI : 10.5194/se-3-467-2012, 2012

Éditeur(s) : HAL CCSD Elsevier
Résumé : The uppermost mantle in back arc regions is the site of complex interactions between partial melting, melt percolation, and fluid migration. To constrain these interactions and evaluate their consequences on geochemical cycles, we carried out an in situ trace element and water study of a suite of spinel peridotite xenoliths from two regions of the Japan back arc system, Ichinomegata (NE Japan) and Oki-Dogo (SW Japan), using LA-ICPMS and FTIR spectrometry, respectively. This study provides the first full dataset of trace element and hydrogen compositions in peridotites including analyses of all their main constitutive silicate minerals: olivine, orthopyroxene and clinopyroxene. The Ichinomegata peridotites sample a LREE-depleted refractory mantle (Mg# olivine = 0.90; Cr# spinel = 0.07–0.23; Yb clinopyroxene = 7.8–13.3 × C1-chondrite, and La/Yb clinopyroxene = 0.003–0.086 × C1-chondrite), characterized by Th-U positive anomalies and constant values of Nb/Ta. The composition of the studied Ichinomegata samples is consistent with that of an oceanic mantle lithosphere affected by cryptic metasomatic interactions with hydrous/aqueous fluids (crypto-hydrous metasomatism). In contrast, the Oki-Dogo peridotites have low Mg# olivine (0.86–0.93) and a broad range of compositions with clinopyroxene showing “spoon-shaped” to flat, and LREE-enriched patterns. They are also characterized by their homogeneous compositions in the most incompatible LILE (e.g., Rb clinopyroxene = 0.01–0.05 × primitive mantle) and HFSE (e.g., Nb clinopyroxene = 0.01–2.16 × primitive mantle). This characteristic is interpreted as resulting from various degrees of melting and extensive melt-rock interactions. FTIR spectroscopy shows that olivine in both Ichinomegata and Oki-Dogo samples has low water contents ranging from 2 to 7 ppm wt. H2O. In contrast, the water contents of pyroxenes from Ichinomegata peridotites (113–271 ppm wt. H2O for orthopyroxene, and 292–347 ppm wt. H2O for clinopyroxene) are significantly higher than in Oki-Dogo peridotites (9–35 ppm wt. H2O for orthopyroxene, and 15–98 ppm wt. H2O for clinopyroxene). This indicates a relationship between melt-rock interaction and water concentrations in pyroxenes. Our study suggests that the water content of the Japan mantle wedge is controlled by the late melt/fluid/rock interactions evidenced by trace element geochemistry: a mechanism triggered by magma-rock interactions may have acted as an efficient dehydrating process in the Oki-Dogo region while the Ichinomegata mantle water content is controlled by slab-derived crypto-hydrous metasomatism.
ISSN: 0016-7037

hal-01553105
https://hal.archives-ouvertes.fr/hal-01553105
DOI : 10.1016/j.gca.2017.03.042

Éditeur(s) : HAL CCSD Elsevier
Résumé : International audience
Before attaining the mantle wedge, where they trigger partial melting, volatiles released from dehydration reactions in the slab have to migrate across a relatively cold (<750 degrees C), peridotite-layer above the incoming slab. In order to unravel the mechanisms allowing for this initial stage of fluid transport, we performed a detailed field and microstructural study of metamorphic prograde peridotites in the Cerro del Almirez ultramafic massif (Betic Cordillera, Spain), where evidences of one of the most important dehydration reactions in subduction zones, the high-pressure antigorite breakdown (P = 1.6-1.9 GPa and T approximate to 680 degrees C), can be mapped in the field. This reaction led to arborescent growth of centimeter-size olivine and orthopyroxene, producing a chlorite-harzburgite with a spinifex-like texture. Microstructural observations and crystal preferred (CPO) mapping show no evidences of solid-state deformation during the prograde growth of olivine and orthopyroxene at the expenses of antigorite. However, a few tens to a hundred meters away from the reaction front, the metamorphic texture is partially obliterated by grain-size reduction in roughly planar conjugate zones, a few mm to meters wide. Grain size reduction zones (GSRZ) are characterized by (1) sharp contacts with undeformed spinifex-like texture domains, (2) important reduction of the olivine grain size (60-250 mu m), (3) olivine color change from brownish to colorless, (4) decrease in the modal amount of orthopyroxene, and (5) at the mm- to cm-scale, irregular shapes and abrupt terminations. Field and microstructural observations exclude that relative displacement took place across these GSRZ. Changes in modal composition imply reactions with fluids undersaturated in silica. Analysis of olivine crystal-preferred orientations (CPO) in GSRZ shows patterns similar, but more dispersed, than those in neighboring spinifex-like domains. It also reveals mm- to cm-scale discrete domains with rather homogeneous crystallographic orientations suggesting inheritance from the preexisting spinifex-like olivines in the host peridotite. Misorientation angles between neighboring grains in the GSRZ show peaks at similar to 5-10 degrees and similar to 20 degrees, but rotations are not crystallographically controlled. Based on these observations, we rule out the formation of the GSRZ by dynamic recrystallization during dislocation creep and propose that they record brittle deformation (microcraking) of the spinifex-like chlorite-harzburgite, probably induced by hydrofracturing at high pressure and relative low temperature conditions (680-710 degrees C). High-pressure hydrofracturing can, thus, be invoked as an efficient mechanism for fluid flow across the cold top-slab mantle layer, hence allowing the slab-derived fluids to ingress in the wedge.
ISSN: 0012-821X

hal-00534035
https://hal.archives-ouvertes.fr/hal-00534035
DOI : 10.1016/j.epsl.2010.06.029

Éditeur(s) : HAL CCSD American Geophysical Union
Résumé : International audience
This work presents seismic properties derived from a detailed microstructural study of supra-subduction mantle xenoliths extracted by the active Avacha volcano in the southern Kamchatka arc. These peridotites underwent low-stress deformation and syn- to late-kinematic reactive percolation of hydrous melts or fluids leading to local orthopyroxene enrichment and dispersion of the olivine CPO in the shallow mantle wedge. All samples show an axial [ 100] olivine CPO implying fast S-wave polarization parallel to the flow direction in the mantle. Orthopyroxene-enrichment by reactive fluids/melts percolation lowers the Vp/Vs ratio, but cannot explain the Vp/Vs ratios <1.7 mapped in the Japan and Andes forearc mantle. These low Vp/Vs ratios may however be explained by considering the intrinsic anisotropy of the peridotites.
ISSN: 0094-8276

hal-00512636
https://hal.archives-ouvertes.fr/hal-00512636
DOI : 10.1029/2010GL043450

Éditeur(s) : HAL CCSD Elsevier
Résumé : International audience
We report chlorine (Cl) and fluorine (F) abundances in minerals, interstitial glasses, and melt inclusions in 12 andesite-hosted, spinel harzburgite xenoliths and crosscutting pyroxenite veins exhumed from the sub-arc lithospheric mantle beneath Avacha volcano in the Kamchatka Arc (NE Russia). The data are used to calculate equilibrium mineral-melt partition coefficients (Dmineral/meltDmineral/melt) for Cl and F relevant to subduction-zone processes and unravel the history of volatile depletion and enrichment mechanisms in an arc setting.Chlorine is ∼100 times more incompatible in pyroxenes (View the MathML sourceDClmineral/melt = 0.005–0.008 [±0.002–0.003]) than F (View the MathML sourceDFmineral/melt = 0.50–0.57 [±0.21–0.24]), which indicates that partial melting of mantle sources leads to strong depletions in Cl relative to F in the residues. The data set in this study suggests a strong control of melt composition on View the MathML sourceDCl,Fpyroxene/melt, in particular H2O contents and Al/(Al + Si), which is in line with recent experiments. Fluorine is compatible in Ca-amphibole in the ‘wet’ sub-arc mantle (View the MathML sourceDFamphibole/melt = 3.5–3.7 [±1.5]) but not Cl (View the MathML sourceDClamphibole/melt = 0.03–0.05 [±0.01–0.03]), indicating that amphibole may fractionate F from Cl in the mantle wedge. The inter-mineral partition coefficients for Cl and F in this study are consistent amongst different harzburgite samples, whether they contain glass or not. In particular, disseminated amphibole hosts much of the Cl and F bulk rock budgets of spinel harzburgites (View the MathML sourceDClamphibole/pyroxene up to 14 and View the MathML sourceDFamphibole/pyroxene up to 40). Chlorine and fluorine are variably enriched (up to 1500 ppm Cl and 750 ppm F) in the parental arc picrite and boninite melts of primitive pyroxenite veins (and related melt inclusions) crosscutting spinel harzburgites.Based on the data in this study, the main inferences on the behaviour of Cl and F during melting and metasomatic processes in the sub-arc mantle are as follow:(i) Melting models show that most depleted mantle protoliths of intra-oceanic arc sources can have extremely low Cl/F (0.002–0.007) before being overprinted by subduction-derived components.(ii) Chlorine has a higher percolation distance in the mantle than F. Even for small fluid or melt volumes, Cl and F signatures of partial melting are overprinted by those of pervasive percolation, which increases Cl/F in percolating agents and bulk peridotites during chromatographic interaction and/or amphibole-forming metasomatic reactions. These processes ultimately control the bulk Cl and F compositions of the residual mantle lithosphere beneath arcs, and likely in other tectonic settings.(iii) Fluxed melting models suggest that Cl enrichment in arc picrite and boninite melts in this study, and in many arc melt inclusions reported in the literature, could be related to the infiltration of high Cl/F fluids derived from subducted serpentinite or altered crust in mantle wedge sources. However, these high Cl/F signatures should be re-evaluated with new models in light of the possible overprint of pervasive percolation effects in the mantle. The breakdown of amphibole (and/or mica) in the deep metasomatised mantle at higher pressure and temperature conditions than in the slab may explain, at least in part, the positive correlations between F abundances and Cl/F in primitive arc melt inclusions and slab depth.
ISSN: 0016-7037

hal-01553610
https://hal.archives-ouvertes.fr/hal-01553610
DOI : 10.1016/j.gca.2016.10.035

Éditeur(s) : HAL CCSD Elsevier
Résumé : International audience
Major and trace element signatures and Sr-Nd-Pb isotope data for muscovite (Ms)-bearing amphibolite blocks and associated muscovite-bearing trondhjemite and quartz-muscovite rocks from the Sierra del Convento melange (eastern Cuba) indicate that Proto-Caribbean oceanic crust underwent wet partial melting processes during Mesozoic subduction and after accretion to the upper plate. Trace element normalized patterns of Ms-bearing amphibolites are enriched in light rare earth elements (LREE) and large-ion lithophile elements (LILE) and evidence variable trace element transfer from the Proto-Caribbean subducting slab to the mantle wedge. Ms-bearing trondhjemites show LREE enrichment and HREE depletion and have geochemical features similar to adakites, including SiO2>56 wt.%, high Na2O contents (5.5-9.0 wt.%) and high Sr/Y (16-644). We consider that the trondhjemites represent primary and natural melts formed by deep partial melting of the subducting slab which did not significantly react with the mantle wedge before intrusion in the subduction channel. The Ms-bearing trondhjemites show different geochemical and petrological signatures compared with the Ms-free tonalites-trondhjemites from Sierra del Convent, that are interpreted as primary slab melts. These differences support the idea that partial melting processes in the Sierra del Convento subduction channel were triggered by the infiltration of fluids derived from three distinct subducted sources: sediments, altered mid-ocean ridge basalts (MORB), and serpentinites. In this scenario, the pegmatitic quartz-muscovite rocks, which are highly enriched in LILE, probably represent the crystallization products of fluids derived by differentiation of the trondhjemitic melts.
ISSN: 0024-4937

hal-00669881
https://hal.archives-ouvertes.fr/hal-00669881
DOI : 10.1016/j.lithos.2011.07.011

Éditeur(s) : HAL CCSD Elsevier
Résumé : In this study, I examine the strain modes of the forearc, arc and back-arc domains in arc-continent collisional settings leading to arc material subduction, delamination and/or accretion. The study focusses on two well-documented colliding island arcs: the Izu–Bonin–Mariana (IBM) arc in Japan and the Luzon arc in Taiwan, both carried by the Philippine Sea plate. Firstly, there is a body of evidence that both the IBM and the Luzon arcs were built on the same Late Jurassic to Early Cretaceous “proto-Philippine Sea Plate” crust. Their internal structure is thus more heterogeneous than expected from Paleogene or Neogene supposedly “intra-oceanic” island arcs. Secondly, those arc systems and proximal “back-arcs” have similar seismic characteristics attesting either for the presence of a middle crust with continental velocities and/or serpentinized uppermost mantle that facilitate crustal shortening/slivering and subsequent decoupling from the rest of the subducting plate. It is shown that the proximal back-arc domain (called “rear-arc” in case of paleoarc activity), overlying the mantle wedge and the subducting slab, may lose its strength if slab-derived hydration occur. Decoupling then occurs below the Moho. Arc delamination likely occurs in mid-crustal levels because middle-crust, heated by nearby magmatism, becomes weak. Accretion of arc material onto the upper plate depends on the characteristics of the arc itself and the geodynamic configuration. Most of the accreted material is probably underplated rather than frontally accreted.
ISSN: 0743-9547

hal-01243080
https://hal.archives-ouvertes.fr/hal-01243080
DOI : 10.1016/j.jseaes.2013.07.043

Éditeur(s) : HAL CCSD Oxford University Press (OUP)
Résumé : International audience
Island arc picrites and boninites are magnesian magmatic rocks believed to be generated by high degrees of melting of depleted mantle sources fluxed by subduction-derived, volatile-rich components. These magmas can be probes of both the mantle wedge protoliths and subduction components, but are rare among other, usually more evolved, types of arc lavas. Furthermore, many arc picrites and boninites show evidence for late-stage differentiation prior to or during eruption, masking their primary, mantle-derived geochemical signatures. We report textural and chemical data on spinel-hosted melt inclusions of mantle origin in amphibole-bearing websterite veins cross-cutting spinel harzburgite xenoliths from the active andesitic Avacha volcano (south Kamchatka, Russia). The data are used to constrain the composition and origin of melts that formed the websterite veins in the sub-arc lithospheric mantle. The melt inclusions typically contain euhedral orthopyroxene and clinopyroxene and occasionally minor amphibole in silicate glass. The melt inclusions were homogenized using heating stages and gas-mixing furnaces. The homogenized glasses range from subalkaline primitive silica-rich picrite and high-Ca boninite (>15 wt % MgO, 48–54 wt % SiO2) to rhyolite. High-Ca boninite glasses have moderate volatile and low heavy rare earth element contents and elevated Cs, Rb, Ba, U, Sr, and Li abundances, with extremely high U/Th. In turn, the glasses display no negative spikes in the high field strength elements Nb, Ta, Zr, Hf, and Ti. We show that the silica-rich picrite and high-Ca boninite liquids in this study formed by high degrees of melting (>25%), at volatile under-saturation, of hybrid melt-depleted but silica-rich mantle sources at ≥1·5 GPa. The hybrid sources formed in two stages: first, by extraction of ∼15% melt from the convecting mantle to form a refractory protolith, which was subsequently enriched in silica via interaction with subduction-derived components prior to or during remelting in the mantle wedge. The subduction-derived components were enriched in fluid-mobile elements and probably oxidized. Overall, our results suggest that silica-rich picrites and high-Ca boninites can be primary melts in mature subduction zones and differentiate within the mantle wedge and the deep arc crust to form more evolved andesite magmas.
ISSN: 0022-3530

hal-01553676
https://hal.archives-ouvertes.fr/hal-01553676
DOI : 10.1093/petrology/egw066

Éditeur(s) : HAL CCSD Elsevier
Résumé : The ultramafic–mafic complexes located at the base of obducted paleo-island arcs are commonly interpreted as evidence for intra-crustal fractionation of primitive (high-Mg#) mantle melts. The present study, realized on the Jijal mafic–ultramafic basal section of the Cretaceous Kohistan oceanic arc (N Pakistan), discards a crystal fractionation model from a single parental magma to take into account the geochemical and isotope variations observed between the plutonic crust and the underlying ultramafic section. The basal ultramafic rocks, i.e. dunites, wehrlites and clinopyroxenites, show high Mg#, extremely depleted REE patterns, no marked HFSE anomalies and high LILE enrichment. Cpx and leached whole rocks from this unit yield restricted εNd values (+ 5.5 to + 6.8), a large εSr range (− 5.9 to + 8.8) and mainly radiogenic 207Pb/204Pb ratios. The overlying mafic section, i.e. gabbroic rocks, displays low Mg# (46 < Mg# < 54), enriched REE patterns and marked HFSE anomalies. These plutonic rocks have also very homogeneous Sr and Pb isotopic compositions defining a restricted domain, which does not overlap that of the ultramafic samples. Moreover, Cpx from the basal ultramafics yield an Sm–Nd isochron at 117 ± 7 Ma interpreted as dating crystallization of the ultramafic section. Together, these results indicate that the mafic rocks on one hand and the ultramafic rocks on the other hand, originate from separate sources with specific emplacement underneath the subduction zone. A lherzolite lens, collected from the top of the mafic section, plots on the early Cretaceous Sm–Nd isochron and therefore yields the same 143Nd/144Nd initial ratio. This suggests that the ultramafic part of the lower Kohistan arc complex sampled the initiation stages of the subduction with production of magmas with boninite-like features. At this stage, the mantle wedge has been already modified by percolation of the first slab-derived fluids. On such a scenario the lherzolite lens would correspond to the lithospheric mantle underneath the arc before its initiation. A multi-stage history is proposed for the evolution for the Kohistan arc through not, vert, similar 30 Ma. Stage #1 reflects the spontaneous initiation of subduction associated with extensive boninitic affinity magmatism, stages #2 and #3 are related to arc building by tholeiitic melts and, a fourth step (stage #4) corresponds to the intra-crustal differentiation.
ISSN: 0012-821X

hal-00407206
https://hal.archives-ouvertes.fr/hal-00407206
DOI : 10.1016/j.epsl.2007.06.026

Éditeur(s) : HAL CCSD Elsevier
Résumé : International audience
We investigated the response of the crystal structure, lattice parameters, and unit-cell volume of hydrous layered silicate phlogopite at conditions relevant to subduction zone settings. We have used first principles simulation based on density functional theory to calculate the equation of state and full elastic constant tensor. Based on the generalized gradient approximation, the full single crystal elastic constant tensor with monoclinic symmetry shows significant anisotropy with the compressional elastic constants: c11 = 181 GPa, c22 = 185 GPa, c33 = 62 GPa, the shear elastic constants c44 = 14 GPa, c55 = 20 GPa, c66 = 68 Ga, and c46 = −6 GPa; the off diagonal elastic constants c12 = 48 GPa, c13 = 12 GPa, c23 = 12 GPa, c15 = −16 GPa, c25 = −5 GPa and c35 = −1 GPa at zero pressure. The elastic anisotropy of phlogopite is larger than most of the layered hydrous phases relevant in the subduction zone conditions. The shear anisotropy, AVS for phlogopite is ∼77% at zero pressure condition and although it decreases upon compression it remains relatively high compared to other hydrous phases relevant in the subduction zone settings. We also note that the shear elastic constants for phlogopite are relatively low. Phlogopite also has a high isotropic bulk VP/VS ratio ∼2.0. However, the VP/VS ratio also exhibits significant anisotropy with values as low as 1.49. Thus, phlogopite bearing metasomatized mantle could readily explain unusual VP/VS ratio as observed from seismological studies from the mantle wedge regions of the subduction zone.
ISSN: 0031-9201

hal-01115670
https://hal.archives-ouvertes.fr/hal-01115670
DOI : 10.1016/j.pepi.2014.05.004

Éditeur(s) : HAL CCSD Elsevier
Résumé : International audience
Despite decades of petrological and geochemical studies, the nature and setting of obducted ophiolites remain controversial: the influence of supra-subduction zone environments on pre-existing oceanic lithosphere is yet to assess, and the processes leading to subduction/obduction initiation are still poorly constrained. Our study documents successive influx of slab-derived fluids and progressive strain localization within the upper mantle in a supra-subduction environment during the first few My of the subduction history. We focus on strongly sheared mafic amphibolites intruding peridotites near the mantle–crust transition of the New Caledonia obducted ophiolite and ~ 50 to 100 m above the basal thrust contact of the ophiolite. These m- to hm-long and several m-thick shear bands are interpreted as inherited small-scale intrusions of mafic melts, probably dikes or sills, which were derived from a moderately refractory mantle source refertilized by supra-subduction zone fluids. 40Ar/39Ar age constraints on pargasite at ca. 90 Ma suggest that they could be inherited from the former Pacific west-dipping subduction.Secondary deformation of these mafic intrusions is intimately associated to three major stages of fluid infiltration: (1) the first stage of deformation and metasomatism is marked by syn-kinematic growth of Ca-amphibole (at 700–800 °C and 3–5 kbar) with a distinctive supra-subduction zone signature, and controlled later channelization of aqueous fluids. 40Ar/39Ar dating on magnesio-hornblende indicates that this deformation episode occurred at ca. 55 Ma, coincident with east-dipping subduction initiation; (2) the main metasomatic stage, characterized by the development of a phlogopite-rich matrix wrapping peridotites and amphibolite boudins, points to the percolation of alkali-rich aqueous fluids at still high temperature (650–750 °C); (3) the last, low temperature (< 600 °C) metasomatic stage results in the formation of deformed veinlets containing talc, chlorite and serpentine.
ISSN: 0024-4937

hal-01306664
https://hal.archives-ouvertes.fr/hal-01306664
https://hal.archives-ouvertes.fr/hal-01306664/document
https://hal.archives-ouvertes.fr/hal-01306664/file/Soret_2016_Strain_localization.pdf
DOI : 10.1016/j.lithos.2015.11.022

Éditeur(s) : HAL CCSD Elsevier
Résumé : International audience
The present study aims at better deciphering the different mechanisms involved in the functioning of the subduction interplate. A 2D thermo-mechanical model is used to simulate a subduction channel, made of oceanic crust, free to evolve. Convergence at constant rate is imposed under a 100 km thick upper plate. Pseudo-brittle and non-Newtonian behaviours are modelled. The influence of the subduction channel strength, parameterized by the difference in activation energy between crust and mantle () is investigated to examine in detail the variations in depth of the subduction plane down-dip extent, . First, simulations show that numerical resolution may be responsible for an artificial and significant shallowing of if the weak crustal layer is not correctly resolved. Second, if the age of the subducting plate is 100 Myr, subduction occurs for any . The stiffer the crust is, that is, the lower is, the shallower is (60 km depth if kJ/mol) and the hotter the fore-arc base is. Conversely, imposing a very weak subduction channel ( J/mol) leads there to an extreme mantle wedge cooling and inhibits mantle melting in wet conditions. Partial kinematic coupling at the fore-arc base occurs if kJ/mol. If the incoming plate is 20 Myr old, subduction can occur under the conditions that the crust is either stiff and denser than the mantle, or weak and buoyant. In the latter condition, cold crust plumes rise from the subduction channel and ascend through the upper lithosphere, triggering (1) partial kinematic coupling under the fore-arc, (2) fore-arc lithosphere cooling, and (3) partial or complete hindrance of wet mantle melting. then ranges from 50 to more than 250 km depth and is time-dependent if crust plumes form. Finally, subduction plane dynamics is intimately linked to the regime of subduction-induced corner flow. Two different intervals of are underlined: 80–120 kJ/mol to reproduce the range of slab surface temperature inferred from geothermometry, and 10–40 kJ/mol to reproduce the shallow hot mantle wedge core inferred from conditions of last equilibration of near-primary arc magmas and seismic tomographies. Therefore, an extra process controlling mantle wedge dynamics is needed to satisfy simultaneously the aforementioned observations. A mantle viscosity reduction, by a factor 4–20, caused by metasomatism in the mantle wedge is proposed. From these results, I conclude that the subduction channel down-dip extent, , should depend on the subduction setting, to be consistent with the observed variability of sub-arc depths of the subducting plate surface.
ISSN: 0031-9201

hal-01623155
https://hal.archives-ouvertes.fr/hal-01623155
DOI : 10.1016/j.pepi.2017.05.008