Petrologic and geochemical role of serpentinite in subduction zones and plate interface domains

2015 ◽  
Vol 37 ◽  
pp. 61-64
Author(s):  
Marco Scambelluri ◽  
Enrico Cannaò ◽  
Mattia Gilio ◽  
Marguerite Godard
Author(s):  
Anne-Aziliz Pelleter ◽  
Gaëlle Prouteau ◽  
Bruno Scaillet

Abstract We performed phase equilibrium experiments on a natural Ca-poor pelite at 3 GPa, 750-1000 °C, under moderately oxidizing conditions, simulating the partial melting of such lithologies in subduction zones. Experiments investigated the effect of sulphur addition on phase equilibria and compositions, with S contents of up to ∼ 2.2 wt. %. Run products were characterized for their major and trace element contents, in order to shed light on the role of sulphur on the trace element patterns of melts produced by partial melting of oceanic Ca-poor sediments. Results show that sulphur addition leads to the replacement of phengite by biotite along with the progressive consumption of garnet, which is replaced by an orthopyroxene-kyanite assemblage at the highest sulphur content investigated. All Fe-Mg silicate phases produced with sulphur, including melt, have higher MgO/(MgO+FeO) ratios (relative to S-free/poor conditions), owing to Fe being primarily locked up by sulphide in the investigated redox range. Secular infiltration of the mantle wedge by such MgO and K2O-rich melts may have contributed to the Mg and K-rich character of the modern continental crust. Addition of sulphur does not affect significantly the stability of the main accessory phases controlling the behaviour of trace elements (monazite, rutile and zircon), although our results suggest that monazite solubility is sensitive to S content at the conditions investigated. The low temperature (∼ 800 °C) S-bearing and Ca-poor sediment sourced slab melts show Th and La abundances, Th/La systematics and HFSE signatures in agreement with the characteristics of sediment-rich arc magmas. Because high S contents diminish phengite and garnet stabilities, S-rich and Ca-poor sediment sourced slab melts have higher contents of Rb, B, Li (to a lesser extent), and HREE. The highest ratios of La/Yb are observed in sulphur-poor runs (with a high proportion of garnet, which retains HREE) and beyond the monazite out curve (which retains LREE). Sulphides appear to be relatively Pb-poor and impart high Pb/Ce ratio to coexisting melts, even at high S content. Overall, our results show that Phanerozoic arc magmas from high sediment flux margins owe their geochemical signature to the subduction of terrigenous, sometimes S-rich, sediments. In contrast, subduction of such lithologies during Archean appears unlikely or unrecorded.


2021 ◽  
Author(s):  
Nadaya Cubas ◽  
Philippe Agard ◽  
Roxane Tissandier

<p>Predicting the spatial extent of mega-earthquakes is an essential ingredient of earthquake hazard assessment. In subduction zones, this prediction mostly relies on geodetic observations of interseismic coupling. However, such models face spatial resolution issues and are of little help to predict full or partial ruptures of highly locked patches. Coupling models are interpreted in the framework of the rate-and-state friction laws. However, these models are too idealized to take into account the effects of a geometrically or rheologically complex plate interface. In this study, we show, from the critical taper theory and a mechanical analysis of the topography, that all recent mega-earthquakes of the Chilean subduction zone are surrounded by distributed interplate deformation emanating from either underplating or basal erosion. This long-lived plate interface deformation builds up stresses ultimately leading to earthquake nucleation. Earthquakes then propagate along a relatively smooth surface and are stopped by segments of heterogeneously distributed deformation. Our results are consistent with long-term features of the subduction margin, with observed short-term deformation as well as physical parameters of recovered subducted fragments. They also provide an explanation for the apparent mechanical segmentation of the megathrust, reconciling many seemingly contradictory observations on the short- and long-term deformation. Consequently, we propose that earthquake segmentation relates to the distribution of deformation along the plate interface and that slip deficit patterns reflect the along-dip and along-strike distribution of the plate interface deformation. Topography would therefore mirror plate interface deformation and could serve to improve earthquake rupture prediction.</p>


Lithos ◽  
2019 ◽  
Vol 334-335 ◽  
pp. 69-87 ◽  
Author(s):  
Hugues Raimbourg ◽  
Vincent Famin ◽  
Giulia Palazzin ◽  
Asuka Yamaguchi ◽  
Romain Augier ◽  
...  
Keyword(s):  

2020 ◽  
Author(s):  
Mario D'Acquisto ◽  
Matthew Herman ◽  
Rob Govers

<div> <p>During and after a large megathrust earthquake, the overriding plate above the rupture zone moves oceanward. Enigmatically, the post-seismic motion of the overriding plate after several recent large earthquakes, further along strike from the rupture zone, was faster in the landward direction than before the event. Previous studies interpreted these changes as the result of increased mechanical coupling along the megathrust interface, transient slab acceleration, or bulk postseismic deformation with elastic bending mentioned as a possible underlying mechanism. Before invoking additional mechanisms, it is important to understand the contribution of postseismic deformation processes that are inherent features of megathrust earthquakes. We thus aim to quantify and analyse the deformation that produces landward motion during afterslip and viscous relaxation. </p> </div><div> <p>We use velocity-driven 3D mechanical finite element models, in which large megathrust earthquakes occur periodically on the finite plate interface. The model geometry is similar to most present-day subduction zones, but does not exactly match any specific subduction zone. </p> </div><div> <p>The results show increased post-seismic landward motion at (trench-parallel) distances greater than 450 km from the middle of the ruptured asperity. Similar patterns of landward motion are generated by viscous relaxation in the mantle wedge and by deep afterslip on the shear zone downdip of the brittle megathrust interface. Landward displacement due to postseismic relaxation largely accumulates at exponentially decaying rates until ~6 Maxwell relaxation times after the earthquake. The spatial distribution and magnitude of the velocity changes is broadly consistent with observations related to both the 2010 Maule and the 2011 Tohoku-oki earthquakes.  </p> </div><div> <p>Further model experiments show that patterns of landward motion due to afterslip and to viscous relaxation are insensitive to the locking pattern of the megathrust. However, the locking distribution does affect the magnitudes of the displacements and velocities. Results show that the increased landward displacement due to postseismic deformation scales directly proportionally to seismic moment. </p> </div><div> <p>We conclude that the landward motion results from in-plane horizontal bending of the overriding plate and mantle. This bending is an elastic response to oceanward tractions near the base of the plate around the ruptured asperity, causing extension locally and compression further away along-trench. This elastic in-plate bending consistently contributes to earthquake-associated changes in surface velocities for the biggest megathrust earthquakes, producing landward motion along strike from the rupture zone.</p> </div>


2020 ◽  
Author(s):  
Whitney Behr ◽  
Adam Holt ◽  
Thorsten Becker ◽  
Claudio Faccenna

<p>Tectonic plate velocities predominantly result from a balance between the potential energy change of the subducting slab and viscous dissipation in the mantle, bending lithosphere, and slab–upper plate interface. A range of observations suggest that slabs may be weak, implying a more prominent role for plate interface dissipation than previously thought. Behr & Becker (2018) suggested that the deep interface viscosity in subduction zones should be strongly affected by the relative proportions of sedimentary to mafic rocks that are subducted to depth, and that sediment subduction should thus facilitate faster subduction plate speeds. Here we use fully dynamic 2D subduction models built with the code ASPECT to quantitatively explore how subduction interface viscosity influences: a) subducting plate sinking velocities, b) trench migration rates, c) convergence velocities, d) upper plate strain regimes, e) dynamic topography, and f) interactions with the 660 km mantle transition zone.  We implement two main types of models, including 1) uniform interface models where interface viscosity and slab strength are systematically varied, and 2) varying interface models where a low viscosity sediment strip of finite width is embedded within a higher viscosity interface. Uniform interface models indicate that low viscosity (sediment-lubricated) slabs have substantially faster sinking velocities prior to reaching the 660, especially for weak slabs, and also that they achieve faster ‘steady state’ velocities after 660 penetration. Even models where sediments are limited to a strip on the seafloor show accelerations in convergence rates of up to ~5 mm/y per my, with convergence initially accommodated by trench rollback and later by slab sinking. We discuss these results in the context of well-documented plate accelerations in Earth’s history such as India-Asia convergence and convergence rate oscillations along the Andean margin.</p><p>References: Behr, W. M., & Becker, T. W. (2018). Sediment control on subduction plate speeds. <em>Earth and Planetary Science Letters</em>, <em>502</em>, 166-173.</p>


2020 ◽  
Author(s):  
Frederique Rolandone ◽  
Jean-Mathieu nocquet ◽  
Patricia Mothes ◽  
Paul Jarrin ◽  
Mathilde Vergnolle

<p>In subduction zones, slip along the plate interface occurs in various modes including earthquakes, steady slip, and transient accelerated aseismic slip during either Slow Slip Events (SSE) or afterslip. We analyze continuous GPS measurements along the central Ecuador subduction segment to illuminate how the different slip modes are organized in space and time in the zone of the 2016 Mw 7.8 Pedernales earthquake. The early post-seismic period (1 month after the earthquake) shows large and rapid afterslip developing at discrete areas of the megathrust and a slow slip event remotely triggered (∼100 km) south of the rupture of the Pedernales earthquake. We find that areas of large and rapid early afterslip correlate with areas of the subduction interface that had hosted SSEs in years prior to the 2016 earthquake. Areas along the Ecuadorian margin hosting regular SSEs and large afterslip had a dominant aseismic slip mode that persisted throughout the earthquake cycle during several years and decades: they regularly experienced SSEs during the interseismic phase, they did not rupture during the 2016 Pedernales earthquake, they had large aseismic slip after it. Four years after the Pedernales earthquake, postseismic deformation is still on-going. Afterslip and SSEs are both involved in the postseimsic deformation. Two large aftershocks (Mw 6.7 & 6.8) occurred after the first month of postseismic deformation in May 18, and later in July 7 2016 two other large aftershocks (Mw 5.9 & 6.3) occurred, all were located north east of the rupture. They may have triggered their own postseismic deformation. Several seismic swarms were identified south and north of the rupture area by a dense network of seismic stations installed during one year after the Pedernales earthquakes, suggesting the occurrence of SSEs. Geodetically, several SSEs were detected during the postseismic deformation either in areas where no SSEs were detected previously, or in areas where regular seismic swarms and repeating earthquakes were identified. The SSEs may have been triggered by the stress increment due to aftershocks or due to afterslip.</p>


2014 ◽  
Vol 376 ◽  
pp. 20-30 ◽  
Author(s):  
Jason Harvey ◽  
Carlos J. Garrido ◽  
Ivan Savov ◽  
Samuele Agostini ◽  
José Alberto Padrón-Navarta ◽  
...  

2019 ◽  
Vol 10 (1) ◽  
Author(s):  
W. P. Schellart ◽  
Z. Chen ◽  
V. Strak ◽  
J. C. Duarte ◽  
F. M. Rosas

Abstract The India-Asia collision has formed the highest mountains on Earth and is thought to account for extensive intraplate deformation in Asia. The prevailing explanation considers the role of the Pacific and Sunda subduction zones as passive during deformation. Here we test the hypothesis that subduction played an active role and present geodynamic experiments of continental deformation that model Indian indentation and active subduction rollback. We show that the synchronous activity and interaction of the collision zone and subduction zones explain Asian deformation, and demonstrate that east-west extension in Tibet, eastward continental extrusion and Asian backarc basin formation are controlled by large-scale Pacific and Sunda slab rollback. The models require 1740 ± 300 km of Indian indentation such that backarc basins form and central East Asian extension conforms estimates. Indentation and rollback produce ~260–360 km of eastward extrusion and large-scale clockwise upper mantle circulation from Tibet towards East Asia and back to India.


2021 ◽  
Author(s):  
Cailey Condit ◽  
Victor Guevara ◽  
Melodie French ◽  
Adam Holt ◽  
Jonathan Delph

<p>Feedbacks amongst petrologic and mechanical processes along the subduction plate boundary play a central role influencing slip behaviors and deformation styles. Metamorphic reactions, resultant fluid production, deformation mechanisms, and strength are strongly temperature dependent, making the thermal structure of these zones a key control on slip behaviors.</p><p> </p><p>Firstly, we investigate the role of metamorphic devolatilization reactions in the production of Episodic Tremor and Slip (ETS) in warm subduction zones. Geophysical and geologic observations of ETS hosting subduction zones suggest the plate interface is fluid-rich and critically stressed, which together, suggests that this area is a zone of near lithostatic pore fluid pressure.  Fluids and high pore fluid pressures have been invoked in many models for ETS. However, whether these fluids are sourced from local dehydration reactions in particular lithologies, or via up-dip transport from greater depths remains an open question. We present thermodynamic models of the petrologic evolution of four lithologies typical of the plate interface along predicted pressure–temperature (P-T) paths for the plate boundary along Cascadia, Nankai, and Mexico which all exhibit ETS at depths between 25-65 km. Our models suggest that 1-2 wt% H<sub>2</sub>O is released at the depths of ETS along these subduction segments due to punctuated dehydration reactions within MORB, primarily through chlorite and/or lawsonite breakdown. These reactions produce sufficient in-situ fluid across this narrow P-T range to cause high pore fluid pressures. Punctuated dehydration of oceanic crust provides the dominant source of fluids at the base of the seismogenic zone in these warm subduction margins, and up-dip migration of fluids from deeper in the subduction zone is not required to produce ETS-facilitating high pore fluid pressures. These dehydration reactions not only produce metamorphic fluids at these depths, but also result in an increased strength of viscous deformation through the breakdown of weak hydrous phases (e.g., chlorite, glaucophane) and the growth of stronger minerals (e.g., garnet, omphacite, Ca-amphibole). Lastly, we present preliminary data on viscosity along warm subduction paths showing the locations of these dehydration pulses correlate with viscosity increases in mafic lithologies along the shallow forarc.</p>


Sign in / Sign up

Export Citation Format

Share Document