Generation of arc basalt magmas and thermal structure of the mantle wedge in subduction zones

1983 ◽  
Vol 88 (B7) ◽  
pp. 5815-5825 ◽  
Author(s):  
Yoshiyuki Tatsumi ◽  
Masanori Sakuyama ◽  
Hiroyuki Fukuyama ◽  
Ikuo Kushiro
Keyword(s):  
Subduction Zones ◽  
Mantle Wedge ◽  
Thermal Structure

Exhumation of subducted mafic rocks in a dynamically evolving thermal structure: constraints from phase equilibria modelling

2021 ◽  
Author(s):  
Rilla C. McKeegan ◽  
Victor E. Guevara ◽  
Adam F. Holt ◽  
Cailey B. Condit
Keyword(s):  
Phase Equilibria ◽  
Subduction Zone ◽  
Subduction Zones ◽  
Mantle Wedge ◽  
Dynamic Models ◽  
Thermal Structure ◽  
Temporal Dependence ◽  
Subduction Channel ◽  
Mafic Rocks ◽  
Point Of No Return

<p>The dominant mechanisms that control the exhumation of subducted rocks and how these mechanisms evolve through time in a subduction zone remain unclear. Dynamic models of subduction zones suggest that their thermal structures evolve from subduction initiation to maturity. The series of metamorphic reactions that occur within the slab, resultant density, and buoyancy with respect to the mantle wedge will co-evolve with the thermal structure. We combine dynamic models of subduction zone thermal structure with phase equilibria modeling to place constraints on the dominant controls on the depth limits of exhumation. This is done across the temporal evolution of a subduction zone for various endmember lithologic associations observed in exhumed high-pressure terranes: sedimentary and serpentinite mélanges, and oceanic tectonic slices.</p><p>Initial modeling suggests that both serpentinite and sedimentary mélanges remain positively buoyant with respect to the mantle wedge throughout all stages of subduction (up to 65 Myr), and for the spectrum of naturally constrained ratios of mafic blocks to serpentinite/sedimentary matrix. In these settings, exhumation depth limits and the “point of no return” (c. 2.3 GPa) are not directly limited by buoyancy, but potentially rheological changes in the slab at the blueschist-eclogite transition stemming from: the switch from amphibole-dominated to pyroxene-dominated rheology and/or dehydration embrittlement. These mechanisms may increase the possibility of brittle failure and hence promote detachment of the slab top into the subduction channel. For the range of temperatures recorded by exhumed serpentinite mélanges, the locus of dehydration for altered MORB at the slab top coincides with the point of no return (2.3 GPa) between 35 and 40 Myr, suggesting a strong temporal dependence on deep exhumation in the subduction channel. </p><p>Tectonic slices composed of 50% mafic rocks and 50% serpentinized slab mantle show a temporal dependence on the depth limits of positive buoyancy. For the range of temperatures recorded by exhumed tectonic slices, the upper pressure limit of positive buoyancy is ~2 GPa, and is only crossed between ~30 and 40 Myr after subduction initiation. Some exhumed tectonic slices record much higher pressures (2.5 GPa); thus, other mechanisms or lithologic combinations may also play a significant role in determining the exhumation limits of tectonic slices. </p><p>Future work includes constraining how the loci of dehydration vary through time for different degrees of oceanic crust alteration, how exhumation limits and mechanisms may change with different subducting plate ages, and calculating how initial exhumation velocities may vary through time. Further comparison with the rock record will constrain the parameters that control the timing and limits of exhumation in subduction zones.</p>

scholarly journals Three-dimensional thermal structure of subduction zones: effects of obliquity and curvature

2012 ◽  
Vol 4 (2) ◽  
pp. 919-941 ◽  
Author(s):  
A. K. Bengtson ◽  
P. E. van Keken
Keyword(s):  
Cross Section ◽  
Cross Sections ◽  
Subduction Zones ◽  
Mantle Wedge ◽  
Thermal Structure ◽  
Critical Role ◽  
Three Dimensional ◽  
Arc Volcanism ◽  
Oblique Subduction ◽  
Dehydration Reactions

Abstract. Quantifying the precise thermal structure of subduction zones is essential for understanding the nature of metamorphic dehydration reactions, arc volcanism, and intermediate depth seismicity. High resolution two-dimensional (2-D) models have shown that the rheology of the mantle wedge plays a critical role and establishes strong temperature gradients in the slab. The influence of three-dimensional (3-D) subduction zone geometry on thermal structure is however not yet well characterized. A common assumption for 2-D models is that the cross-section is taken normal to the strike of the trench with a corresponding velocity reduction in the case of oblique subduction, rather than taken parallel to velocity. A comparison between a full 3-D Cartesian model with oblique subduction and selected 2-D cross-sections demonstrates that the trench-normal cross-section provides a better reproduction of the slab thermal structure than the velocity-parallel cross-section. An exception is found in the case of strongly curved subduction, such as in the Marianas, where strong 3-D flow in the mantle wedge is generated. In this case it is shown that the full 3-D model should be evaluated for an accurate prediction of the slab thermal structure.

Numerical simulation of subduction zone pressure-temperature-time paths: constraints on fluid production and arc magmatism

1991 ◽  
Vol 335 (1638) ◽  
pp. 341-353 ◽  
Keyword(s):  
Heat Transfer ◽  
Partial Melting ◽  
Subduction Zone ◽  
Subduction Zones ◽  
Mantle Wedge ◽  
Thermal Structure ◽  
Oceanic Lithosphere ◽  
Transfer Model ◽  
Temperature Time ◽  
Subducting Slab

The location and sequence of metamorphic devolatilization and partial melting reactions in subduction zones may be constrained by integrating fluid and rock pressure-temperature-time ( P-T-t ) paths predicted by numerical heat-transfer models with phase diagrams constructed for metasedimentary, metabasaltic, and ultramafic bulk compositions. Numerical experiments conducted using a two-dimensional heat transfer model demonstrate that the primary controls on subduction zone P-T-t paths are: (1) the initial thermal structure; (2) the amount of previously subducted lithosphere; (3) the location of the rock in the subduction zone; and (4) the vigour of mantle wedge convection induced by the subducting slab. Typical vertical fluid fluxes out of the subducting slab range from less than 0.1 to 1 (kg fluid) m -2 a -1 for a convergence rate of 3 cm a -1 . Partial melting of the subducting, amphibole-bearing oceanic crust is predicted to only occur during the early stages of subduction initiated in young (less than 50 Ma) oceanic lithosphere. In contrast, partial melting of the overlying mantle wedge occurs in many subduction zone experiments as a result of the infiltration of fluids derived from slab devolatilization reactions. Partial melting in the mantle wedge may occur by a twostage process in which amphibole is first formed by H 2 O infiltration and subsequently destroyed as the rock is dragged downward across the fluid-absent ‘hornblende-out’ partial melting reaction.

scholarly journals Role of warm subduction in the seismological properties of the forearc mantle: An example from southwest Japan

2021 ◽  
Vol 7 (28) ◽  
pp. eabf8934
Author(s):  
Changyeol Lee ◽  
YoungHee Kim
Keyword(s):  
Subduction Zones ◽  
Mantle Wedge ◽  
Numerical Models ◽  
Thermal Structure ◽  
Free Fluid ◽  
Southwest Japan ◽  
Seismic Properties ◽  
Delay Times ◽  
Forearc Mantle ◽  
Nonvolcanic Tremors

A warm slab thermal structure plays an important role in controlling seismic properties of the slab and mantle wedge. Among warm subduction zones, most notably in southwest Japan, the spatial distribution of large S-wave delay times and deep nonvolcanic tremors in the forearc mantle indicate the presence of a serpentinite layer along the slab interface. However, the conditions under which such a layer is generated remains unclear. Using numerical models, we here show that a serpentinite layer begins to develop by the slab-derived fluids below the deeper end of the slab-mantle decoupling interface and grows toward the corner of the mantle wedge along the interface under warm subduction conditions only, explaining the large S-wave delay times in the forearc mantle. The serpentinite layer then allows continuous free-fluid flow toward the corner of the mantle wedge, presenting possible mechanisms for the deep nonvolcanic tremors in the forearc mantle.

scholarly journals Advances in the thermal and petrologic modeling of subduction zones

Geosphere ◽  
2020 ◽  
Vol 16 (4) ◽  
pp. 936-952 ◽  
Author(s):  
Simon M. Peacock
Keyword(s):  
Heat Flux ◽  
Subduction Zones ◽  
Mantle Wedge ◽  
Dynamic Models ◽  
Seismic Velocity ◽  
Thermal Structure ◽  
Metamorphic Rocks ◽  
Ridge Subduction ◽  
Dehydration Reactions ◽  
Flux Measurements

Abstract In the two decades since Subduction: Top to Bottom was published in 1996, improved analytical and numerical thermal-petrologic models of subduction zones have been constructed and evaluated against new seismological and geological observations. Advances in thermal modeling include a range of new approaches to incorporating shear (frictional, viscous) heating along the subduction interface and to simulating induced flow in the mantle wedge. Forearc heat-flux measurements constrain the apparent coefficient of friction (μ′) along the plate interface to <∼0.1, but the extent to which μ′ may vary between subduction zones remains challenging to discern owing to scatter in the heat-flux measurements and uncertainties in the magnitude and distribution of radiogenic heat production in the overriding crust. Flow in the mantle wedge and the resulting thermal structure depend on the rheology of variably hydrated mantle rocks and the depth at which the subducting slab becomes coupled to the overlying mantle wedge. Advances in petrologic modeling include the incorporation of sophisticated thermodynamic software packages into thermal models and the prediction of seismic velocities from mineralogic and petrologic models. Current thermal-petrologic models show very good agreement between the predicted location of metamorphic dehydration reactions and observed intermediate-depth earthquakes, and between the predicted location of the basalt-to-eclogite transition in subducting oceanic crust and observed landward-dipping, low-seismic-velocity layers. Exhumed high-pressure, low-temperature metamorphic rocks provide insight into subduction-zone temperatures, but important thermal parameters (e.g., convergence rate) are not well constrained, and metamorphic rocks exposed at the surface today may reflect relatively warm conditions in the past associated with subduction initiation or ridge subduction. We can anticipate additional advances in our understanding of subduction zones as a result of further testing of model predictions against geologic and geophysical observations, and of evaluating the importance of advective processes, such as diapirism and subduction-channel flow, that are not captured in hybrid kinematic-dynamic models of subduction zones but are observed in fully dynamical models under certain conditions.

scholarly journals Reconciling mantle wedge thermal structure with arc lava thermobarometric determinations in oceanic subduction zones

2016 ◽  
Vol 17 (10) ◽  
pp. 4105-4127 ◽  
Author(s):  
Alexander Perrin ◽  
Saskia Goes ◽  
Julie Prytulak ◽  
D. Rhodri Davies ◽  
Cian Wilson ◽  
...  
Keyword(s):  
Subduction Zones ◽  
Mantle Wedge ◽  
Thermal Structure

scholarly journals Three-dimensional thermal structure of subduction zones: effects of obliquity and curvature

Solid Earth ◽  
2012 ◽  
Vol 3 (2) ◽  
pp. 365-373 ◽  
Author(s):  
A. K. Bengtson ◽  
P. E. van Keken
Keyword(s):  
Cross Section ◽  
Cross Sections ◽  
Subduction Zones ◽  
Mantle Wedge ◽  
Thermal Structure ◽  
Critical Role ◽  
Three Dimensional ◽  
Arc Volcanism ◽  
Oblique Subduction ◽  
Dehydration Reactions

Abstract. Quantifying the precise thermal structure of subduction zones is essential for understanding the nature of metamorphic dehydration reactions, arc volcanism, and intermediate depth seismicity. High resolution two-dimensional (2-D) models have shown that the rheology of the mantle wedge plays a critical role and establishes strong temperature gradients in the slab. The influence of three-dimensional (3-D) subduction zone geometry on thermal structure is however not yet well characterized. A common assumption for 2-D models is that the cross-section is taken normal to the strike of the trench with a corresponding velocity reduction in the case of oblique subduction, rather than taken parallel to velocity. A comparison between a full 3-D Cartesian model with oblique subduction and selected 2-D cross-sections demonstrates that the trench-normal cross-section provides a better reproduction of the slab thermal structure than the velocity-parallel cross-section. An exception is found in the case of a strongly curved trench, such as in the Marianas, where strong 3-D flow in the mantle wedge is generated. In this case it is shown that the full 3-D model should be evaluated for an accurate prediction of the slab thermal structure. The models demonstrate that the use of a dynamic slab and wedge, separated by a kinematic boundary, yields good results for describing slab velocities in 3-D.

scholarly journals Ophiolitic Pyroxenites Record Boninite Percolation in Subduction Zone Mantle

Minerals ◽  
2019 ◽  
Vol 9 (9) ◽  
pp. 565 ◽  
Author(s):  
Véronique Le Roux ◽  
Yan Liang
Keyword(s):  
Subduction Zone ◽  
Subduction Zones ◽  
Mantle Wedge ◽  
Major Element ◽  
Early Stage ◽  
Diffusion Modeling ◽  
Melt Infiltration ◽  
Melt Percolation ◽  
Back Arc ◽  
Parental Melts

The peridotite section of supra-subduction zone ophiolites is often crosscut by pyroxenite veins, reflecting the variety of melts that percolate through the mantle wedge, react, and eventually crystallize in the shallow lithospheric mantle. Understanding the nature of parental melts and the timing of formation of these pyroxenites provides unique constraints on melt infiltration processes that may occur in active subduction zones. This study deciphers the processes of orthopyroxenite and clinopyroxenite formation in the Josephine ophiolite (USA), using new trace and major element analyses of pyroxenite minerals, closure temperatures, elemental profiles, diffusion modeling, and equilibrium melt calculations. We show that multiple melt percolation events are required to explain the variable chemistry of peridotite-hosted pyroxenite veins, consistent with previous observations in the xenolith record. We argue that the Josephine ophiolite evolved in conditions intermediate between back-arc and sub-arc. Clinopyroxenites formed at an early stage of ophiolite formation from percolation of high-Ca boninites. Several million years later, and shortly before exhumation, orthopyroxenites formed through remelting of the Josephine harzburgites through percolation of ultra-depleted low-Ca boninites. Thus, we support the hypothesis that multiple types of boninites can be created at different stages of arc formation and that ophiolitic pyroxenites uniquely record the timing of boninite percolation in subduction zone mantle.

Sign in / Sign up

Export Citation Format

Share Document