Physical conditions producing slab stagnation: Constraints of the Clapeyron slope, mantle viscosity, trench retreat, and dip angles

Recent seismic tomography has revealed various morphologies in the subducted lithosphere. In particular, significant flattening and stagnation of slabs around the 660-km boundary are seen in some areas beneath the northwestern Pacific subduction zones. We examined the cause of slab stagnation in terms of the Clapeyron slope of the phase transformation from ringwoodite to perovskite + magnesiowüstite, trench retreat velocity, dip angles, and high viscosity of the lower mantle based on two-dimensional (2-D) numerical simulations of thermal convection. In particular, we examined the conditions necessary for slab stagnation assuming a very small absolute value of the Clapeyron slope, which were proposed based on recent high-pressure, high-temperature (high P–T) experiments. Our calculations show that slabs tend to stagnate above the 660-km boundary with an increasing absolute value of the Clapeyron slope, viscosity jump at the boundary, and trench retreat velocity and a decreasing initial dip angle. Stagnant slabs could be obtained numerically for a realistic range of parameters obtained from high P–T experiments and other geophysical observations combining buoyancy, high lower-mantle viscosity, and trench retreat. We found that a low dip angle of a descending slab at the bottom of the upper mantle plays an important role in slab stagnation. Two main regimes underlie slab stagnation: buoyancy-dominated and viscosity-dominated regimes. In the viscosity-dominated regime, it is possible for slabs to stagnate above the 660-km boundary, even when the value of the Clapeyron slope is 0 MPa/K.     

Model calculations of the thermal fields of subducting lithospheric slabs and partial melting

A new numerical model that simulates a downgoing slab is used to study the conditions required to produce melting on its upper surface. Models with dip angles of 26.6° and 45°, rates of subduction of 0.7 and 5.6 cm y⁻¹, varying heat sources and rising material from the top of the slab are included. The results indicate that melting will not be greatly affected by dip angle, though the rate of subduction and the amount of shearstrain heating are important. When melting occurs, material rising from the top of the slab may produce high heat flow values at the surface of the earth on the continental side of the ocean trench, if the process continues sufficiently long. The sinking slab produces a positive gravity anomaly on the continental side of subduction, which is reduced in amplitude when rising material is present.     

Global trench migration velocities and slab migration induced upper mantle volume fluxes: Constraints to find an Earth reference frame based on minimizing viscous dissipation

Since the advent of plate tectonics different global reference frames have been used to describe the motion of plates and trenches. The difference in plate motion and trench migration between different reference frames can be substantial (up to 4 cm/yr). This study presents an overview of trench migration velocities for all the mature and incipient subduction zones on Earth as calculated in eight different global reference frames. Calculations show that, irrespective of the reference frame: (1) trench retreat always dominates over trench advance, with 62–78% of the 244 trench segments retreating; (2) the mean and median trench velocity are always positive (retreating) and within the range 1.3–1.5 cm/yr and 0.9–1.3 cm/yr, respectively; (3) rapid trench retreat is only observed close to lateral slab edges (< 1500 km); and (4) trench retreat is always slow far from slab edges (> 2000 km). These calculations are predicted by geodynamic models with a varying slab width, in which plate motion, trench motion and mantle flow result from subduction of dense slabs, suggesting that trench motion is indeed primarily driven by slab buoyancy forces and that proximity to a lateral slab edge exerts a dominant control on the trench migration velocity. Despite these four general conclusions, significant differences in velocities between such reference frames remain. It is therefore important to determine which reference frame most likely describes the true absolute velocities to get an understanding of the forces driving plate tectonics and mantle convection. It is here proposed that, based on fluid dynamic considerations and predictions from geodynamic modelling, the best candidate is the one, which optimises the number of trench segments that retreat, minimizes the trench–perpendicular trench migration velocity (v_T) in the centre of wide (> 4000 km) subduction zones, maximizes the number of retreating trench segments located within 2000 km of the closest lateral slab edge, minimizes the average of the absolute of the trench–perpendicular trench migration velocity (|v_T|) for all subduction zones on Earth, and minimizes the global upper mantle toroidal volume flux (_To) that results from trench migration and associated lateral slab migration (i.e. slab rollback or slab roll-forward). Calculations show that these conditions are best met in one particular Indo-Atlantic hotspot reference frame, where 75% of the subduction zones retreat, v_T in the centre of wide subduction zones ranges between − 3.5 and 1.8 cm/yr, 83% of the trench segments located within 2000 km of the closest lateral slab edge retreat, the average of |v_T| is 2.1 cm/yr, and _To = 456 km³/yr (lower limit) and 539 km³/yr (upper limit). Inclusion of all the incipient subduction zones on Earth results in slightly greater fluxes of 465 km³/yr (lower limit) and 569 km³/yr (upper limit). It is also found that this reference frame is close to minimizing the total sub-lithospheric upper mantle volume flux (_K) associated with motion of continental keels located below the major cratons. It is stressed, however, that _K is an order of magnitude smaller than _To, and thus of subordinate importance. In conclusion, the Indo-Atlantic hotspot reference frame appears preferable for calculating plate velocities and plate boundary velocities.     

Structural characteristics off Miyagi forearc region, the Japan Trench seismogenic zone, deduced from a wide-angle reflection and refraction study

The Japan Trench subduction zone, located east of NE Japan, has regional variation in seismicity. Many large earthquakes occurred in the northern part of Japan Trench, but few in the southern part. Off Miyagi region is in the middle of the Japan Trench, where the large earthquakes (M > 7) with thrust mechanisms have occurred at an interval of about 40 years in two parts: inner trench slope and near land. A seismic experiment using 36 ocean bottom seismographs (OBS) and a 12,000 cu. in. airgun array was conducted to determine a detailed, 2D velocity structure in the forearc region off Miyagi. The depth to the Moho is 21 km, at 115 km from the trench axis, and becomes progressively deeper landward. The P-wave velocity of the mantle wedge is 7.9–8.1 km/s, which is typical velocity for uppermost mantle without large serpentinization. The dip angle of oceanic crust is increased from 5–6° near the trench axis to 23° 150 km landward from the trench axis. The P-wave velocity of the oceanic uppermost mantle is as small as 7.7 km/s. This low-velocity oceanic mantle seems to be caused by not a lateral anisotropy but some subduction process. By comparison with the seismicity off Miyagi, the subduction zone can be divided into four parts: 1) Seaward of the trench axis, the seismicity is low and normal fault-type earthquakes occur associated with the destruction of oceanic lithosphere. 2) Beneath the deformed zone landward of the trench axis, the plate boundary is characterized as a stable sliding fault plain. In case of earthquakes, this zone may be tsunamigenic. 3) Below forearc crust where P-wave velocity is almost 6 km/s and larger: this zone is the seismogenic zone below inner trench slope, which is a plate boundary between the forearc and oceanic crusts. 4) Below mantle wedge: the rupture zones of thrust large earthquakes near land (e.g. 1978 off Miyagi earthquake) are located beneath the mantle wedge. The depth of the rupture zones is 30–50 km below sea level. From the comparison, the rupture zones of large earthquakes off Miyagi are limited in two parts: plate boundary between the forearc and oceanic crusts and below mantle wedge. This limitation is a rare case for subduction zone. Although the seismogenic process beneath the mantle wedge is not fully clarified, our observation suggests the two possibilities: earthquake generation at the plate boundary overridden by the mantle wedge without serpentinization or that in the subducting slab.     

Comment on “The potential influence of subduction zone polarity on overriding plate deformation,trench migration and slab dip angle” by W.P. Schellart

The Schellart's [Schellart, W.P., 2007, The potential influence of subduction zone polarity on overriding plate deformation, trench migration and slab dip angle. Tectonophysics, 445, 363–372.] paper uses slab dip and upper plate extension for testing the westward drift. His analysis and discussion are misleading for the study of the net rotation of the lithosphere since the first 125 km of subduction zones are sensitive also to other parameters such upper plate thickness, geometry and obliquity of the subduction zone with respect to the convergence direction. The deeper (> 125 km) part cannot easily be compared as well because E- or NE-directed subduction zones have seismic gaps between 270–630 km. Moreover the velocity of subduction hinge cannot be precisely estimated and it does not equal to backarc spreading due to accretionary prism growth and asthenospheric intrusion at the subduction hinge. It is shown here that hinge migration in the upper plate or lower plate reference frames supports a general global polarization of the lithosphere in agreement with the westward drift of the lithosphere. The W-directed subduction zones appear controlled by the slab–mantle interaction with slab retreat imposed by the eastward mantle flow. The opposite E-NE-directed subduction zones seem rather mainly controlled by the convergence rate, plus density, thickness and viscosity of the upper and lower plates. Finally, the geological and geophysical asymmetries recorded along subduction and rift zones as a function of their polarity with respect to the tectonic mainstream are not questioned in the Schellart's paper, but they rather represent the basic evidence for the westward drift of the lithosphere.     

Influence of the precollisional stage on subduction dynamics and the buried crust thermal state: Insights from numerical simulations

At continental subduction initiation, the continental crust buoyancy may induce, first, a convergence slowdown, and second, a compressive stress increase that could lead to the forearc lithosphere rupture. Both processes could influence the slab surface P–T conditions, favoring on one side crust partial melting or on the opposite the formation of ultra-high pressure/low temperature (UHP-LT) mineral. We quantify these two effects by performing numerical simulations of subduction. Water transfers are computed as a function of slab dehydration/overlying mantle hydration reactions, and a strength decrease is imposed for hydrated mantle rocks. The model starts with an old oceanic plate ( 100 Ma) subducting for 145.5 Myr with a 5 cm/yr convergence rate. The arc lithosphere is thermally thinned between 100 km and 310 km away from the trench, due to small-scale convection occuring in the water-saturated mantle wedge. We test the influence of convergence slowdown by carrying on subduction with a decreased convergence rate (≤ 2 cm/yr). Surprisingly, the subduction slowdown yields not only a strong slab warming at great depth (> 80 km), but also a significant cooling of the forearc lithosphere at shallower depth. The convergence slowdown increases the subducted crust temperature at 90 km depth to 705 ± 62 °C, depending on the convergence rate reduction, and might thus favor the oceanic crust partial melting in presence of water. For subduction velocities ≤ 1 cm/yr, slab breakoff is triggered 20–32 Myr after slowdown onset, due to a drastic slab thermal weakening in the vicinity of the interplate plane base. At last, the rupture of the weakened forearc is simulated by imposing in the thinnest part of the overlying lithosphere a dipping weakness plane. For convergence with rates ≥ 1 cm/yr, the thinned forearc first shortens, then starts subducting along the slab surface. The forearc lithosphere subduction stops the slab surface warming by hot asthenosphere corner flow, and decreases in a first stage the slab surface temperature to 630 ± 20 °C at 80 km depth, in agreement with P–T range inferred from natural records of UHP-LT metamorphism. The subducted crust temperature is further reduced to 405 ± 10 °C for the crust directly buried below the subducting forearc. Such a cold thermal state at great depth has never been sampled in collision zones, suggesting that forearc subduction might not be always required to explain UHP-LT metamorphsim.     

Comparison of gravity anomaly between mature and immature intra-oceanic subduction zones in the western Pacific

The western Pacific hosts major subduction systems such as Izu–Bonin–Mariana and Tonga–Kermadec, but also less conspicuous systems such as Yap, Mussau and Hjort trenches which constitute the young, incomplete, or ultraslow-member in the evolutionary spectrum of subduction zones. We used satellite-derived gravity data to compare well-developed and immature subduction systems. It is shown that at spatial resolution > 10–20 km or so, the satellite data have accuracy comparable to ship-board gravity measurements over intra-oceanic subduction zones. In the isostatic residual gravity anomaly map, the width of non-isostatically-compensated region of the mature subduction zones is much wider than that of immature ones. More importantly, when the gravitational attraction due to seafloor is removed, a large difference exists between the mature and immature subduction zones in the overriding plate side. Mature subduction zones exhibit broad low gravity anomalies of ~ 200–250 mGal centered at distances of 150–200 km from the trench which are not found over immature subduction zones. The cause of the broad low gravity anomalies over mature subduction zones is debatable due to lack of information on the deep crust and upper mantle structure and property. We discuss the following four causes: (1) serpentinization of the upper mantle beneath the forearc; (2) presence of partial melt in the mantle wedge caused by release of volatiles from the slab, frictional heating and distributed by mantle circulation; (3) difference in density structure between the overriding and subducting plates caused by difference in age and thermal structures with and without compositional stratification between crust and mantle; and (4) anomalous thickness of the arc not explained by isostasy. Our analysis suggests that serpentinization cannot explain the observed gravity anomaly which appears ~ 150–200 km from the trench. Although the extent and distribution of partial melt within the mantle wedge remain in question, to our best estimate, partial melting contributes little (< 50 mGal) to the total negative gravity anomaly. The difference in density structure reflecting temperature difference can only explain less than half of the low gravity anomaly. The sinking of lighter crustal material produces a large negative anomaly in the forearc but its location does not match the observed gravity anomaly. It appears that one cannot explain the total difference in gravity anomaly without invoking anomalous thickness of the arc. Although we could not identify the sole or combination of factors that give rise to the low gravity anomaly in mature subduction zones, the comparison of gravity anomalies between mature and immature subduction zones is likely to provide an important constraint for understanding the evolution and structure of subduction zones as more complementary evidences become available.     

Variations of slab dip and overriding plate tectonics during subduction: Insights from analogue modelling

We present small-scale laboratory models of oceanic subduction in which plates motion is imposed by lateral boundary conditions. The oceanic plate moves trenchward at constant speed and subducts below a fixed overriding plate. In this configuration, the long-term process of subduction is not steady-state. Slab interaction with the upper mantle-lower mantle boundary results in periods of slab flattening during which the dip of the slab diminishes, followed by periods of slab steepening. The overriding plate tectonic regime is influenced by the dynamics of subduction, slab anchoring favouring trench perpendicular shortening. When the slab is anchored, slab flattening further favours shortening, while slab steepening favours extension or smaller shortening rates. Non-steady-state long-term subduction may explain part of the variability of slab geometries evidenced by statistical analyses of present-day subduction zones. Experiments suggest that, despite boundary conditions applied on the converging plates do not change, tectonics pulses within the overriding plate may be caused by this non-steady-state behaviour.     

Subduction of the Woodlark Basin at New Britain Trench, Solomon Islands region

The Woodlark Basin, located south of the Solomon Islands arc region, is a young (5 Ma) oceanic basin that subducts beneath the New Britain Trench. This region is one of only a few subduction zones in the world where it is possible to study a young plate subduction of several Ma. To obtain the image of the subducting slab at the western side of the Woodlark Basin, a 40-day Ocean Bottom Seismometer (OBS) survey was conducted in 1998 to detect the micro-seismic activity. It was the first time such a survey had been performed in this location and over 600 hypocenters were located. The seismic activity is concentrated at the 10–60 km depth range along the plate boundary. The upper limit just about coincides with the leading edge of the accretionary wedge. The upper limit boundary was identified as the up-dip limit of the seismogenic zone, whereas the down-dip limit of the seismogenic zone was difficult to define. The dip angle of the plate at the high seismicity zone was found to average about 30°. Using the Cascadia subduction zone for comparison, which is a typical example of a young plate subduction, suggests that the subduction of the Woodlark Basin was differentiated by a high dip angle and rather landward location of the seismic front from the trench axis (30 km landward from the trench axis). Furthermore, as pointed out by previous researchers, the convergent margin of the Solomon Islands region is imposed with a high stress state, probably due to the collision of the Ontong Java Plateau and a rather rapid convergence rate (10 cm/year). The results of the high angle plate subduction and inner crust earthquakes beneath the Shortland Basin strongly support the high stress state. The collision of the Ontong Java Plateau, the relatively rapid convergence rate, and moderately cold slab as evidenced by low heat flow, rather than the plate age, may be dominantly responsible for the geometry of the seismogenic zone in the western part of the Woodlark Basin subduction zone.     

Seismicity, gravity anomalies and lithospheric structure of the Andaman arc, NE Indian Ocean

The Andaman arc in the northeastern Indian Ocean defines nearly 1100 km long active plate margin between the India and Burma plates where an oblique Benioff zone develops down to 200 km depth. Several east-trending seismologic sections taken across the Andaman Benioff Zone (ABZ) are presented here to detail the subduction zone geometry in a 3-D perspective. The slab gravity anomaly, computed from the 3-D ABZ configuration, is a smooth, long-wavelength and symmetric gravity high of 85 mGal amplitude centering to the immediate east of the Nicobar Island, where, a prominent gravity “high” follows the Nicobar Deep. The Slab-Residual Gravity Anomaly (SRGA) and Mantle Bouguer Anomaly (MBA) maps prepared for the Andaman plate margin bring out a double-peaked SRGA “low” in the range of − 150 to − 240 mGal and a wider-cum-larger MBA “low” having the amplitude of − 280 to − 315 mGal demarcating the Andaman arc–trench system. The gravity models provide evidences for structural control in propagating the rupture within the lithosphere. The plate margin configuration below the Andaman arc is sliced by the West Andaman Fault (WAF) as well as by a set of sympathetic faults of various proportions, often cutting across the fore-arc sediment package. Some of these fore-arc thrust faults clearly give rise to considerably high post-seismic activity, but the seismic incidence along the WAF further east is comparatively much less particularly in the north, although, the lack of depth resolution for many of the events prohibits tracing the downward continuity of these faults. Tectonic correlation of the gravity-derived models presented here tends to favour the presence of oceanic crust below the Andaman–Nicobar Outer Arc Ridge.     

Numerical modeling of thermal structure, circulation of H2O, and magmatism–metamorphism in subduction zones: Implications for evolution of arcs

Numerical models on thermal structure, convective flow of solid, generation and transportation of H₂O-rich fluid in subduction zones are consolidated to have a comprehensive view of the subduction zone processes: heat balance, circulation of H₂O magmatism–metamorphism, growth of arcs and continental margins. A large scale convection model with steady subduction of a cold old slab (130 Myr old) predicts rapid ( 100 Myr) cooling of subduction zones, resulting in cessation of magmatism. The model also predicts that the mantle temperature beneath arcs and continental margins is greatly affected by the effective temperature of the subducting slab, i.e., the age of the subducting slab. If subduction of a young hot slab, including ridge subduction, occurs every 60 to 120 Myr as is suggested for eastern Asia, the average temperature beneath arcs is increased by about 300 °C, which may explain the long-lasting magmatism in eastern Asia. Associated with subduction of young slabs and ridges, thermal structure and circulation of H₂O are greatly modified to cause a transition from (1) normal arc magmatism, (2) forearc mantle melting, to (3) slab melting to produce a significant amount (100 km³) of granitic melts, associated with both high-P/T and low-P/T type metamorphism. The last stage of (3) can result in formation of a granitic batholith belt and a paired metamorphic belts. Synthesis of the numerical models and observations suggest that episodic subduction of young slabs and ridges can explain heat source for generating a large amount of granitic magmas of batholiths, synchronous formation of batholith and regional metamorphic belts, and P–T conditions of the paired metamorphism. Even the high-P/T metamorphism requires an elevated geothermal structure in the forearc region, associated with ridge subduction. Although the emplacement of the batholiths and the regional metamorphic belts, and the mass balance in subduction zones are not well constrained at present, the episodic event associated with ridge subduction is thought to be essential for net growth of arcs and continental margins, as well as for the long-term heat balance in subduction zones.     

On the role of subducting oceanic plateaus in the development of shallow flat subduction

Oceanic plateaus, aseismic ridges or seamount chains all have a thickened crust and their subduction has been proposed as a possible mechanism to explain the occurrence of flat subduction and related absence of arc magmatism below Peru, Central Chile and at the Nankai Trough (Japan). Their extra compositional buoyancy could prohibit the slab from sinking into the mantle. With a numerical thermochemical convection model, we simulated the subduction of an oceanic lithosphere that contains an oceanic crustal plateau of 18-km thickness. With a systematic variation, we examined the required physical parameters to obtain shallow flat subduction. Metastability of the basaltic crust in the eclogite stability field is of crucial importance for the slab to remain buoyant throughout the subduction process. In a 44-Ma-old subducting plate, basalt must be able to survive a temperature of 600–700 °C to keep the plate buoyant sufficiently long to cause a flat-slab segment. We found that the maximum yield stress in the slab must be limited to about 600 MPa to allow for the necessary bending to the horizontal. Young slabs show flat subduction for larger parameter ranges than old slabs, since they are less gravitationally unstable and show less resistance against bending. Hydrous weakening of the mantle wedge area and lowermost continent are required to allow for the necessary deformation of a change in subduction style from steep to flat. The maximum flat slab extent is about 300 km, which is sufficient to explain the observed shallow flat subduction near the Nankai Trough (Japan). However, additional mechanisms, such as active overthrusting by an overriding continental plate, need to be invoked to explain the flat-slab segments up to 500 km long below Peru and Central Chile.     

Neocene and Quaternary extent and geometry of the subducted Pacific Plate beneath North Island, New Zealand: Implications for Kaikoura tectonics

The extent and geometry of the obliquely subduced oceanic Pacific Plate beneath North Island, New Zealand, for five million year intervals through the mid-Miocene to Quaternary, are presented in a series of maps and cross-sections. These show that the subducted plate progressively increased its extent from NE to SW beneath the North Island, and in the more northern regions where it was first emplaced, concomitantly increased its dip from 10° to 50°.The changing extent and geometry of the subducted slab has been established from the age pattern of orogenic andesites and from the geochemical K₂O-h parameter of depth of magma generation. The radiometric dates show a migration of the volcanic front back towards the trench at an average rate of 20 km/My. The trenchward migration is explained by a model of increasing slab dip which is corroborated by the K₂O data calibrated against the presently active arc (Taupo Volcanic Zone). With the exception of northern Coromandel Peninsula, the andesitic magmas were generated at 85–100 km depth. The interpretation of the dates adopted here indicates that the subducted slab originated at the NE-SW trending Kermadec-Hikurangi Trench, and implies a different and much simpler evolution of the Australia-Pacific plate boundary in the vicinity of North Island than other recent models.Subduction geometry has been found elsewhere to be a principal influence upon the state of stress and deformational style in an over-riding plate. The possibility is explored that the timing, nature and pattern of the Neogene to Quaternary Kaikoura Orogeny in North Island is due to this influence. Apart from the effect of oblique subduction in eastern North Island, there is an accord between the onset of deformation and the emplacement sequence of the shallow slab beneath North Island, and between the change in subduction geometry and a progressive north to south change in northern North Island from compression to extension.     

Evidence from the Mesta half-graben, SW Bulgaria, for the Late Eocene beginning of Aegean extension in the Central Balkan Peninsula

The ENE-tilted Mesta half-graben contains a 3-km-thick section of Priabonian (Late Eocene) to Oligocene sedimentary and volcanic rocks that rest unconformably on basement metamorphic rocks along its west side. Basal strata dip 50–60° E and dip at progressively lower angles upward, indicating synrotational deposition. The southern part of the half-graben contains nested volcanic caldera complexes, formed during the deposition of the middle part of the sedimentary sequence, which have been rotated by about half the total rotation of the sedimentary succession. The half-graben is bounded on the east by a fault that steepens from more deeply exposed structural levels in the south (8–18° W) to shallower exposed structural levels in the north (70° W) and together with the rotation of Paleogene strata during deposition indicate the Mesta half-graben is underlain by a listric detachment fault, the Mesta detachment. Subhorizontal Middle Miocene strata that unconformably overlie tilted Paleogene strata yield an upper age limit to the extension. West and northwest of the Mesta half-graben are many other NNW-trending NE-tilted Paleogene half-grabens which we suggest are part of an important extended area in SW Bulgaria and eastern Macedonia that lies above one or more west-dipping detachment faults and date the beginning of Aegean extension in the southern Balkan region as at least as old as Priabonian. The Mesta detachment is oblique to the trend of a contemporaneous Paleogene magmatic arc in the southern Balkans and the origin of the detachment is probably related to gravitationally induced spreading of thickened hot arc crust and Hellenic trench roll back.     

Mantle dynamics of Western Pacific and East Asia: Insight from seismic tomography and mineral physics

Recent results of high-resolution seismic tomography and mineral physics experiments are used to study mantle dynamics of Western Pacific and East Asia. The most important processes in subduction zones are the shallow and deep slab dehydration and the convective circulation (corner flow) processes in the mantle wedge. The combination of the two processes may have caused the back-arc spreading in the Lau basin, affected the morphology of the subducting Philippine Sea slab and its seismicity under southwest Japan, and contributed to the formation of the continental rift system and intraplate volcanism in Northeast Asia, which are clearly visible in our tomographic images. Slow anomalies are also found in the mantle under the subducting Pacific slab, which may represent (a) small mantle plumes, (b) upwellings associated with the slab collapsing down to the lower mantle, or (c) sub-slab dehydration associated with deep earthquakes caused by the reactivation of large faults preserved in the slab. Combining tomographic images and earthquake hypocenters with phase diagrams in the systems of peridotite + water, we proposed a petrologic model for arc volcanism. Arc magmas are caused by the dehydration reactions of hydrated slab peridotite that supply water-rich fluids to the mantle wedge and cause partial melting of the convecting mantle wedge. A large amount of fluids can be released from hydrated MORB at depths shallower than 55 km, which move upwards to hydrate the wedge corner under the fore-arc, and never drag down to the deeper mantle along the slab surface. Slab dehydration reactions at 120 km depth are the antigorite-related 5 reactions which supply water-rich fluids for forming the volcanic front. Phase A and Mg-surssasite breakdown reactions at 200 and 300 km depths below 700 °C cause the second and third arcs, respectively. Moreover, the dehydration reactions of super-hydrous phase B, phases D and E at 500–660 km depths cause the fluid transportation to the mantle boundary layer (MBL) (410–660 km depth). The stagnant slabs extend from Japan to Beijing, China for over 1000 km long, indicating that the arc–trench system covers the entire region from the Japan trench to East Asia. We propose a big mantle wedge (BMW) model herein, where hydrous plumes originating from 410 km depth cause a series of intra-continental hot regions. Fluids derived from MBL accumulated by the double-sided subduction zones, rather than the India–Asia collision and the subsequent indentation into Asia, are the major cause for the active tectonics and mantle dynamics in this broad region.     

Slab–mantle flow interaction: influence on subduction dynamics and duration

We investigate the influence of mantle flow relative to the lithosphere on subduction dynamics. We use 2D thermo‐mechanical models assuming incompressible non‐Newtonian fluid rheology. Different mantle flow velocities consistent with absolute plate motion models are tested, as well as both directions of flow, either sustaining or opposing slab dip. The effects of different inflow/outflow velocity profiles, slab strengths and upper–lower mantle viscosity contrasts are also evaluated. Slab dip deviations between models with opposite mantle flow directions range from 37° for relatively strong slabs (η_max = 10²⁵ Pa s) to 50° for weaker slabs (η_max = 10²⁴ Pa s), accounting for a significant amount of natural slab dip variability. For imposed mantle flow supporting the slab, the initial stage of slab steepening is followed by a stage of continuous slab dip decrease. This slab shallowing eventually leads to mantle wedge closure, subduction cessation and slab break‐off, possibly driving subduction flips.     

Convergence structures of the Peru trench between 10°S and 14°S

High resolution seafloor studies of the Peru Trench between 10°S and 14°S with the GLORIA long-range side-scan sonar system show that the Nazca plate is broken by numerous normal faults as it bends into the trench. These bending-induced faults strike subparallel to the trench axis and overprint and cut across spreading fabric structures of the plate. They commonly form grabens having widths and spacings of 3–5 km and extend for as much as 100 km along strike. Vertical displacements are generally 200 m or more by the time they reach the trench axis. Turbidite deposits are found in the trench north of 11.5°S. Both turbidite and pelagic sediments are folded and temporarily accreted to the base of the overriding plate along the length of the trench axis. They are apparently subsequently implaced in the grabens by slumping and subducted with the Nazca plate. The Mendaña Fracture Zone, which intersects the trench between 9°40′S and 10°35′S, appears to be the locus of a seaward propagating rift that is forming in response to subduction-induced extensional stresses in the Nazca plate.     

The dynamics of big mantle wedge, magma factory, and metamorphic–metasomatic factory in subduction zones

The East Asian continental margin is underlain by stagnant slabs resulting from subduction of the Pacific plate from the east and the Philippine Sea plate from the south. We classify the upper mantle in this region into three major domains: (a) metasomatic–metamorphic factory (MMF), subduction zone magma factory (SZMF), and the ‘big mantle wedge’ (BMW). Whereas the convection pattern is anticlockwise in the MMF domain, it is predominantly clockwise in the SZMF and BMW, along a cross section from the south. Here we define the MMF as a small wedge corner which is driven by the subducting Pacific plate and dominated by H₂O-rich fluids derived by dehydration reactions, and enriched in large ion lithophile elements (LILE) which cause the metasomatism. The SZMF is a zone intermediate between MMF and BMW domains and constitutes the main region of continental crust production by partial melting through wedge counter-corner flow. Large hydrous plume generated at about 200 km depth causes extensive reduction in viscosity and the smaller scale hydrous plumes between 60 km and 200 km also bring about an overall reduction in the viscosity of SZMF. More fertile and high temperature peridotites are supplied from the entrance to this domain. The domain extends obliquely to the volcanic front and then swings back to the deep mantle together with the subducting slab. The BMW occupies the major portion of upper mantle in the western Pacific and convects largely with a clockwise sense removing the eastern trench oceanward. Sporadic formation of hydrous plume at the depth of around 410 km and the curtain flow adjacent to the trench cause back arc spreading. We envisage that the heat source in BMW could be the accumulated TTG (tonalite–trondhjemite–granodiorite) crust on the bottom of the mantle transition zone. The ongoing process of transportation of granitic crust into the mantle transition zone is evident from the deep subduction of five intra-oceanic arcs on the subducting Philippine Sea plate from the south, in addition to the sediment trapped subduction by the Pacific plate and Philippine Sea plate. The dynamics of MMF, SZMF and BMW domains are controlled by the angle of subduction; a wide zone of MMF in SW Japan is caused by shallow angle subduction of the Philippine Sea plate and the markedly small MMF domain in the Mariana trench is due to the high angle subduction of Pacific plate. The domains in NE Japan and Kyushu region are intermediate between these two. During the Tertiary, a series of marginal basins were formed because of the nearly 2000 km northward shift of the subduction zone along the southern margin of Tethyan Asia, which may be related to the collision of India with Asia and the indentation. The volume of upper mantle under Asia was reduced extensively on the southern margin with a resultant oceanward trench retreat along the eastern margin of Asia, leading to the formation of a series of marginal basins. The western Pacific domain in general is characterized by double-sided subduction; from the east by the oldest Pacific plate and from the south by the oldest Indo-Australian plate. The old plates are hence hydrated extensively even in their central domains and therefore of low temperature. The cracks have allowed the transport of water into the deeper portions of the slab and these domains supply hydrous fluids even to the bottom of the upper mantle. Thus, a fluid dominated upper mantle in the western Pacific drives a number of microplates and promote the plate boundary processes.     

An insight into the Bucaramanga nest

We used the local seismicity for the period of 1993 to 2001, in the northeast of Colombia to show the existence of two slabs in the north and south of the Bucaramanga nest. The northern slab has a dip angle of about 25° and the southern slab has a 50° dip angle, while the dip in the Bucaramanga nest is about 29°. In order to explain the nature of the Bucaramanga nest, we proposed the scenario of collision between these two slabs. Using a 3D Finite Element Model (FEM) we show that collision can concentrate, modify and perturb the stress field. The active process of dehydration embrittlement at intermediate depths and the concentrated stress field in the collision zone may explain the high rate of seismic activity inside the Bucaramanga nest. The perturbed and modified stress field resulting from the simultaneous effect of collision between two subducted slabs and subduction of the lithosphere under its own weight can explain the variation in the focal mechanism of micro-earthquakes and the complexity in the source of the moderate size earthquakes in the Bucaramanga nest.     

Pressure–temperature evolution of lawsonite eclogite in Sivrihisar; Tavşanlı Zone–Turkey

The Cretaceous blueschist belt, Tavşanlı Zone, representing the subducted and exhumed northern continental margin of the Anatolide–Tauride platform is exposed in Western Anatolia. The Sivrihisar area east of Tavşanlı is made up of tectonic units consisting of i) metaclastics and conformably overlying massive marbles (coherent blueschist unit), ii) blueschist-eclogite unit, iii) marble–calcschist intercalation and iv) metaperidotite slab. The metaclastics are composed of jadeite–lawsonite–glaucophane and jadeite–glaucophane–chloritoid schists, phengite phyllites, and calcschists with glaucophane–lawsonite metabasite layers. The blueschist-eclogite unit representing strongly sheared, deeply buried and imbricated tectonic slices of accreted uppermost levels of the oceanic crust with minor metamorphosed serpentinite bodies consists of lawsonite-bearing eclogitic metabasites (approximately 90% of the field), lawsonite eclogites, metagabbros, serpentinites, pelagic marbles, omphacite–glaucophane–lawsonite metapelites and metacherts. The mineral assemblage of the lawsonite eclogite (garnet + omphacite > 70%) is omphacite, garnet, lawsonite, glaucophane, phengite and rutile. Lawsonite eclogite lenses are enclosed by garnet–lawsonite blueschist envelopes.Textural evidence from lawsonite eclogites and country rocks reveals that they did not leave the stability field of lawsonite during subduction and exhumation. The widespread preservation of lawsonite in eclogitic metabasites and eclogites can be attributed to rapid subduction and subsequent exhumation in a low geothermal gradient of the oceanic crust material without experiencing a thermal relaxation. Peak P–T conditions of lawsonite eclogites are estimated at 24 ± 1 kbar and 460 ± 25 °C. These P–T conditions indicate a remarkably low geotherm of 6.2 °C/km corresponding to a burial depth of 74 km.