Evolution and diversity of subduction zones controlled by slab width

Subducting slabs provide the main driving force for plate motion and flow in the Earth's mantle, and geodynamic, seismic and geochemical studies offer insight into slab dynamics and subduction-induced flow. Most previous geodynamic studies treat subduction zones as either infinite in trench-parallel extent (that is, two-dimensional) or finite in width but fixed in space. Subduction zones and their associated slabs are, however, limited in lateral extent (250-7,400 km) and their three-dimensional geometry evolves over time. Here we show that slab width controls two first-order features of plate tectonics-the curvature of subduction zones and their tendency to retreat backwards with time. Using three-dimensional numerical simulations of free subduction, we show that trench migration rate is inversely related to slab width and depends on proximity to a lateral slab edge. These results are consistent with retreat velocities observed globally, with maximum velocities (6-16 cm yr(-1)) only observed close to slab edges (<1,200 km), whereas far from edges (>2,000 km) retreat velocities are always slow (<2.0 cm yr(-1)). Models with narrow slabs (< or =1,500 km) retreat fast and develop a curved geometry, concave towards the mantle wedge side. Models with slabs intermediate in width ( approximately 2,000-3,000 km) are sublinear and retreat more slowly. Models with wide slabs (> or =4,000 km) are nearly stationary in the centre and develop a convex geometry, whereas trench retreat increases towards concave-shaped edges. Additionally, we identify periods (5-10 Myr) of slow trench advance at the centre of wide slabs. Such wide-slab behaviour may explain mountain building in the central Andes, as being a consequence of its tectonic setting, far from slab edges.     

Laboratory models of the thermal evolution of the mantle during rollback subduction

The subduction of oceanic lithosphere plays a key role in plate tectonics, the thermal evolution of the mantle and recycling processes between Earth's interior and surface. Information on mantle flow, thermal conditions and chemical transport in subduction zones come from the geochemistry of arc volcanoes, seismic images and geodynamic models. The majority of this work considers subduction as a two-dimensional process, assuming limited variability in the direction parallel to the trench. In contrast, observationally based models increasingly appeal to three-dimensional flow associated with trench migration and the sinking of oceanic plates with a translational component of motion (rollback). Here we report results from laboratory experiments that reveal fundamental differences in three-dimensional mantle circulation and temperature structure in response to subduction with and without a rollback component. Without rollback motion, flow in the mantle wedge is sluggish, there is no mass flux around the plate and plate edges heat up faster than plate centres. In contrast, during rollback subduction flow is driven around and beneath the sinking plate, velocities increase within the mantle wedge and are focused towards the centre of the plate, and the surface of the plate heats more along the centreline.     

西太平洋俯冲带岩石圈变形研究

全球海洋岩石圈的最大弯曲与大地震发生在俯冲带。当弯曲应力超过岩石承受范围，就会产生正断层和地震，海水沿着正断层进入上地幔并发生蛇纹石化，引发浅源地震并可能造成灾难性海啸。选取西太平洋最具代表性的日本、伊豆-小笠原、马里亚纳和雅浦俯冲带以及汤加-克马德克俯冲带，归纳近些年的地球物理观测及地球动力学模拟的结果，对比分析了不同俯冲带挠曲正断层的分布特征，并探讨了俯冲板块变形与地震之间的相关性,以揭示俯冲板片弯曲变形及相应的正断层与潜在板块水化特征。     

Subduction dynamics and the origin of Andean orogeny and the Bolivian orocline

The building of the Andes results from the subduction of the oceanic Nazca plate underneath the South American continent. However, how and why the Andes and their curvature, the Bolivian orocline, formed in the Cenozoic era (65.5?million years (Myr) ago to present), despite subduction continuing since the Mesozoic era (251.0-65.5?Myr ago), is still unknown. Three-dimensional numerical subduction models demonstrate that variations in slab thickness, arising from the Nazca plate's age at the trench, produce a cordilleran morphology consistent with that observed. The age-dependent sinking of the slab in the mantle drives traction towards the trench at the base of the upper plate, causing it to thicken. Thus, subducting older Nazca plate below the Central Andes can explain the locally thickened crust and higher elevations. Here we demonstrate that resultant thickening of the South American plate modifies both shear force gradients and migration rates along the trench to produce a concave margin that matches the Bolivian orocline. Additionally, the varying forcing along the margin allows stress belts to form in the upper-plate interior, explaining the widening of the Central Andes and the different tectonic styles found on their margins, the Eastern and Western Cordilleras. The rise of the Central Andes and orocline formation are directly related to the local increase of Nazca plate age and an age distribution along the margin similar to that found today; the onset of these conditions only occurred in the Eocene epoch. This may explain the enigmatic delay of the Andean orogeny, that is, the formation of the modern Andes.     

Bending-related faulting and mantle serpentinization at the Middle America trench

The dehydration of subducting oceanic crust and upper mantle has been inferred both to promote the partial melting leading to arc magmatism and to induce intraslab intermediate-depth earthquakes, at depths of 50-300 km. Yet there is still no consensus about how slab hydration occurs or where and how much chemically bound water is stored within the crust and mantle of the incoming plate. Here we document that bending-related faulting of the incoming plate at the Middle America trench creates a pervasive tectonic fabric that cuts across the crust, penetrating deep into the mantle. Faulting is active across the entire ocean trench slope, promoting hydration of the cold crust and upper mantle surrounding these deep active faults. The along-strike length and depth of penetration of these faults are also similar to the dimensions of the rupture area of intermediate-depth earthquakes.     

Seismic evidence of negligible water carried below 400-km depth in subducting lithosphere

Strong evidence exists that water is carried from the surface into the upper mantle by hydrous minerals in the uppermost 10-12?km of subducting lithosphere, and more water may be added as the lithosphere bends and goes downwards. Significant amounts of that water are released as the lithosphere heats up, triggering earthquakes and fluxing arc volcanism. In addition, there is experimental evidence for high solubility of water in olivine, the most abundant mineral in the upper mantle, for even higher solubility in olivine's high-pressure polymorphs, wadsleyite and ringwoodite, and for the existence of dense hydrous magnesium silicates that potentially could carry water well into the lower mantle (deeper than 1,000?km). Here we compare experimental and seismic evidence to test whether patterns of seismicity and the stabilities of these potentially relevant hydrous phases are consistent with a wet lithosphere. We show that there is nearly a one-to-one correlation between dehydration of minerals and seismicity at depths less than about 250?km, and conclude that the dehydration of minerals is the trigger of instability that leads to seismicity. At greater depths, however, we find no correlation between occurrences of earthquakes and depths where breakdown of hydrous phases is expected. Lastly, we note that there is compelling evidence for the existence of metastable olivine (which, if present, can explain the distribution of deep-focus earthquakes) west of and within the subducting Tonga slab and also in three other subduction zones, despite metastable olivine being incompatible with even extremely small amounts of water (of the order of 100?p.p.m. by weight). We conclude that subducting slabs are essentially dry at depths below 400?km and thus do not provide a pathway for significant amounts of water to enter the mantle transition zone or the lower mantle.     

Thermochemical structures beneath Africa and the Pacific Ocean

Large low-velocity seismic anomalies have been detected in the Earth's lower mantle beneath Africa and the Pacific Ocean that are not easily explained by temperature variations alone. The African anomaly has been interpreted to be a northwest-southeast-trending structure with a sharp-edged linear, ridge-like morphology. The Pacific anomaly, on the other hand, appears to be more rounded in shape. Mantle models with heterogeneous composition have related these structures to dense thermochemical piles or superplumes. It has not been shown, however, that such models can lead to thermochemical structures that satisfy the geometrical constraints, as inferred from seismological observations. Here we present numerical models of thermochemical convection in a three-dimensional spherical geometry using plate velocities inferred for the past 119 million years. We show that Earth's subduction history can lead to thermochemical structures similar in shape to the observed large, lower-mantle velocity anomalies. We find that subduction history tends to focus dense material into a ridge-like pile beneath Africa and a relatively more-rounded pile under the Pacific Ocean, consistent with seismic observations.     

The depth distribution of azimuthal anisotropy in the continental upper mantle

The most likely cause of seismic anisotropy in the Earth's upper mantle is the lattice preferred orientation of anisotropic minerals such as olivine. Its presence reflects dynamic processes related to formation of the lithosphere as well as to present-day tectonic motions. A powerful tool for detecting and characterizing upper-mantle anisotropy is the analysis of shear-wave splitting measurements. Because of the poor vertical resolution afforded by this type of data, however, it has remained controversial whether the splitting has a lithospheric origin that is 'frozen-in' at the time of formation of the craton, or whether the anisotropy originates primarily in the asthenosphere, and is induced by shear owing to present-day absolute plate motions. In addition, predictions from surface-wave-derived models are largely incompatible with shear-wave splitting observations. Here we show that this disagreement can be resolved by simultaneously inverting surface waveforms and shear-wave splitting data. We present evidence for the presence of two layers of anisotropy with different fast-axis orientations in the cratonic part of the North American upper mantle. At asthenospheric depths (200-400 km) the fast axis is sub-parallel to the absolute plate motion, confirming the presence of shear related to current tectonic processes, whereas in the lithosphere (80-200 km), the orientation is significantly more northerly. In the western, tectonically active, part of North America, the fast-axis direction is consistent with the absolute plate motion throughout the depth range considered, in agreement with a much thinner lithosphere.     

Mid-mantle deformation inferred from seismic anisotropy

With time, convective processes in the Earth's mantle will tend to align crystals, grains and inclusions. This mantle fabric is detectable seismologically, as it produces an anisotropy in material properties--in particular, a directional dependence in seismic-wave velocity. This alignment is enhanced at the boundaries of the mantle where there are rapid changes in the direction and magnitude of mantle flow, and therefore most observations of anisotropy are confined to the uppermost mantle or lithosphere and the lowermost-mantle analogue of the lithosphere, the D" region. Here we present evidence from shear-wave splitting measurements for mid-mantle anisotropy in the vicinity of the 660-km discontinuity, the boundary between the upper and lower mantle. Deep-focus earthquakes in the Tonga-Kermadec and New Hebrides subduction zones recorded at Australian seismograph stations record some of the largest values of shear-wave splitting hitherto reported. The results suggest that, at least locally, there may exist a mid-mantle boundary layer, which could indicate the impediment of flow between the upper and lower mantle in this region.     

A perovskitic lower mantle inferred from high-pressure, high-temperature sound velocity data

The determination of the chemical composition of Earth's lower mantle is a long-standing challenge in earth science. Accurate knowledge of sound velocities in the lower-mantle minerals under relevant high-pressure, high-temperature conditions is essential in constraining the mineralogy and chemical composition using seismological observations, but previous acoustic measurements were limited to a range of low pressures and temperatures. Here we determine the shear-wave velocities for silicate perovskite and ferropericlase under the pressure and temperature conditions of the deep lower mantle using Brillouin scattering spectroscopy. The mineralogical model that provides the best fit to a global seismic velocity profile indicates that perovskite constitutes more than 93 per cent by volume of the lower mantle, which is a much higher proportion than that predicted by the conventional peridotitic mantle model. It suggests that the lower mantle is enriched in silicon relative to the upper mantle, which is consistent with the chondritic Earth model. Such chemical stratification implies layered-mantle convection with limited mass transport between the upper and the lower mantle.     

Arc-parallel flow in the mantle wedge beneath Costa Rica and Nicaragua

Resolving flow geometry in the mantle wedge is central to understanding the thermal and chemical structure of subduction zones, subducting plate dehydration, and melting that leads to arc volcanism, which can threaten large populations and alter climate through gas and particle emission. Here we show that isotope geochemistry and seismic velocity anisotropy provide strong evidence for trench-parallel flow in the mantle wedge beneath Costa Rica and Nicaragua. This finding contradicts classical models, which predict trench-normal flow owing to the overlying wedge mantle being dragged downwards by the subducting plate. The isotopic signature of central Costa Rican volcanic rocks is not consistent with its derivation from the mantle wedge or eroded fore-arc complexes but instead from seamounts of the Galapagos hotspot track on the subducting Cocos plate. This isotopic signature decreases continuously from central Costa Rica to northwestern Nicaragua. As the age of the isotopic signature beneath Costa Rica can be constrained and its transport distance is known, minimum northwestward flow rates can be estimated (63-190 mm yr(-1)) and are comparable to the magnitude of subducting Cocos plate motion (approximately 85 mm yr(-1)). Trench-parallel flow needs to be taken into account in models evaluating thermal and chemical structure and melt generation in subduction zones.     

Seismic reflection imaging of two megathrust shear zones in the northern Cascadia subduction zone

At convergent continental margins, the relative motion between the subducting oceanic plate and the overriding continent is usually accommodated by movement along a single, thin interface known as a megathrust. Great thrust earthquakes occur on the shallow part of this interface where the two plates are locked together. Earthquakes of lower magnitude occur within the underlying oceanic plate, and have been linked to geochemical dehydration reactions caused by the plate's descent. Here I present deep seismic reflection data from the northern Cascadia subduction zone that show that the inter-plate boundary is up to 16 km thick and comprises two megathrust shear zones that bound a >5-km-thick, approximately 110-km-wide region of imbricated crustal rocks. Earthquakes within the subducting plate occur predominantly in two geographic bands where the dip of the plate is inferred to increase as it is forced around the edges of the imbricated inter-plate boundary zone. This implies that seismicity in the subducting slab is controlled primarily by deformation in the upper part of the plate. Slip on the shallower megathrust shear zone, which may occur by aseismic slow slip, will transport crustal rocks into the upper mantle above the subducting oceanic plate and may, in part, provide an explanation for the unusually low seismic wave speeds that are observed there.     

Recycled dehydrated lithosphere observed in plume-influenced mid-ocean-ridge basalt

A substantial uncertainty in the Earth's global geochemical water cycle is the amount of water that enters the deep mantle through the subduction and recycling of hydrated oceanic lithosphere. Here we address the question of recycling of water into the deep mantle by characterizing the volatile contents of different mantle components as sampled by ocean island basalts and mid-ocean-ridge basalts. Although all mantle plume (ocean island) basalts seem to contain more water than mid-ocean-ridge basalts, we demonstrate that basalts associated with mantle plume components containing subducted lithosphere--'enriched-mantle' or 'EM-type' basalts--contain less water than those associated with a common mantle source. We interpret this depletion as indicating that water is extracted from the lithosphere during the subduction process, with greater than 92 per cent efficiency.     

Experimental evidence for the existence of iron-rich metal in the Earth's lower mantle

The oxidation state recorded by rocks from the Earth's upper mantle can be calculated from measurements of the distribution of Fe3+ and Fe2+ between the constituent minerals. The capacity for minerals to incorporate Fe3+ may also be a significant factor controlling the oxidation state of the mantle, and high-pressure experimental measurements of this property might provide important insights into the redox state of the more inaccessible deeper mantle. Here we show experimentally that the Fe3+ content of aluminous silicate perovskite, the dominant lower-mantle mineral, is independent of oxygen fugacity. High levels of Fe3+ are present in perovskite even when it is in chemical equilibrium with metallic iron. Silicate perovskite in the lower mantle will, therefore, have an Fe3+/total Fe ratio of at least 0.6, resulting in a whole-rock ratio of over ten times that of the upper mantle. Consequently, the lower mantle must either be enriched in Fe3+ or Fe3+ must form by the disproportionation of Fe2+ to produce Fe3+ plus iron metal. We argue that the lower mantle contains approximately 1 wt% of a metallic iron-rich alloy. The mantle's oxidation state and siderophile element budget have probably been influenced by the presence of this alloy.     

Redox freezing and melting in the Earth's deep mantle resulting from carbon-iron redox coupling

Very low seismic velocity anomalies in the Earth's mantle may reflect small amounts of melt present in the peridotite matrix, and the onset of melting in the Earth's upper mantle is likely to be triggered by the presence of small amounts of carbonate. Such carbonates stem from subducted oceanic lithosphere in part buried to depths below the 660-kilometre discontinuity and remixed into the mantle. Here we demonstrate that carbonate-induced melting may occur in deeply subducted lithosphere at near-adiabatic temperatures in the Earth's transition zone and lower mantle. We show experimentally that these carbonatite melts are unstable when infiltrating ambient mantle and are reduced to immobile diamond when recycled at depths greater than ～250?kilometres, where mantle redox conditions are determined by the presence of an (Fe,Ni) metal phase. This 'redox freezing' process leads to diamond-enriched mantle domains in which the Fe(0), resulting from Fe(2+) disproportionation in perovskites and garnet, is consumed but the Fe(3+) preserved. When such carbon-enriched mantle heterogeneities become part of the upwelling mantle, diamond will inevitably react with the Fe(3+) leading to true carbonatite redox melting at ～660 and ～250 kilometres depth to form deep-seated melts in the Earth's mantle.     

Stability of magnesite and its high-pressure form in the lowermost mantle

Carbonates are important constituents of marine sediments and play a fundamental role in the recycling of carbon into the Earth's deep interior via subduction of oceanic crust and sediments. Study of the stability of carbonates under high pressure and temperature is thus important for modelling the carbon budget in the entire Earth system. Such studies, however, have rarely been performed under appropriate lower-mantle conditions and no experimental data exist at pressures greater than 80 GPa (refs 3-6). Here we report an in situ X-ray diffraction study of the stability of magnesite (MgCO(3)), which is the major component of subducted carbonates, at pressure and temperature conditions approaching those of the core-mantle boundary. We found that magnesite transforms to an unknown form at pressures above approximately 115 GPa and temperatures of 2,100-2,200 K (depths of approximately 2,600 km) without any dissociation, suggesting that magnesite and its high-pressure form may be the major hosts for carbon throughout most parts of the Earth's lower mantle.     

Evidence for recycled Archaean oceanic mantle lithosphere in the Azores plume

The compositional differences between mid-ocean-ridge and ocean-island basalts place important constraints on the form of mantle convection. Also, it is thought that the scale and nature of heterogeneities within plumes and the degree to which heterogeneous material endures within the mantle might be reflected in spatial variations of basalt composition observed at the Earth's surface. Here we report osmium isotope data on lavas from a transect across the Azores archipelago which vary in a symmetrical pattern across what is thought to be a mantle plume. Many of the lavas from the centre of the plume have lower 187Os/188Os ratios than most ocean-island basalts and some extend to subchondritic 187Os/188Os ratios-lower than any yet reported from ocean-island basalts. These low ratios require derivation from a depleted, harzburgitic mantle, consistent with the low-iron signature of the Azores plume. Rhenium-depletion model ages extend to 2.5 Gyr, and we infer that the osmium isotope signature is unlikely to be derived from Iberian subcontinental lithospheric mantle. Instead, we interpret the osmium isotope signature as having a deep origin and infer that it may be recycled, Archaean oceanic mantle lithosphere that has delaminated from its overlying oceanic crust. If correct, our data provide evidence for deep mantle subduction and storage of oceanic mantle lithosphere during the Archaean era.     

Seismic evidence for catastrophic slab loss beneath Kamchatka

In the northwest Pacific Ocean, a sharp corner in the boundary between the Pacific plate and the North American plate joins a subduction zone running along the southern half of the Kamchatka peninsula with a region of transcurrent motion along the western Aleutian arc. Here we present images of the seismic structure beneath the Aleutian-Kamchatka junction and the surrounding region, indicating that: the subducting Pacific lithosphere terminates at the Aleutian-Kamchatka junction; no relict slab underlies the extinct northern Kamchatka volcanic arc; and the upper mantle beneath northern Kamchatka has unusually slow shear wavespeeds. From the tectonic and volcanic evolution of Kamchatka over the past 10 Myr (refs 3-5) we infer that at least two episodes of catastrophic slab loss have occurred. About 5 to 10 Myr ago, catastrophic slab loss shut down island-arc volcanic activity north of the Aleutian-Kamchatka junction. A later episode of slab loss, since about 2 Myr ago, seems to be related to the activity of the world's most productive island-arc volcano, Klyuchevskoy. Removal of lithospheric mantle is commonly discussed in the context of a continental collision, but our findings imply that episodes of slab detachment and loss are also important agents in the evolution of oceanic convergent margins.     

Continental collision slowing due to viscous mantle lithosphere rather than topography

Because the inertia of tectonic plates is negligible, plate velocities result from the balance of forces acting at plate margins and along their base. Observations of past plate motion derived from marine magnetic anomalies provide evidence of how continental deformation may contribute to plate driving forces. A decrease in convergence rate at the inception of continental collision is expected because of the greater buoyancy of continental than oceanic lithosphere, but post-collisional rates are less well understood. Slowing of convergence has generally been attributed to the development of high topography that further resists convergent motion; however, the role of deforming continental mantle lithosphere on plate motions has not previously been considered. Here I show that the rate of India's penetration into Eurasia has decreased exponentially since their collision. The exponential decrease in convergence rate suggests that contractional strain across Tibet has been constant throughout the collision at a rate of 7.03?×?10(-16)?s(-1), which matches the current rate. A constant bulk strain rate of the orogen suggests that convergent motion is resisted by constant average stress (constant force) applied to a relatively uniform layer or interface at depth. This finding follows new evidence that the mantle lithosphere beneath Tibet is intact, which supports the interpretation that the long-term strain history of Tibet reflects deformation of the mantle lithosphere. Under conditions of constant stress and strength, the deforming continental lithosphere creates a type of viscous resistance that affects plate motion irrespective of how topography evolved.     

The Earth's mantle

Seismological images of the Earth's mantle reveal three distinct changes in velocity structure, at depths of 410, 660 and 2,700 km. The first two are best explained by mineral phase transformations, whereas the third-the D" layer-probably reflects a change in chemical composition and thermal structure. Tomographic images of cold slabs in the lower mantle, the displacements of the 410-km and 660-km discontinuities around subduction zones, and the occurrence of small-scale heterogeneities in the lower mantle all indicate that subducted material penetrates the deep mantle, implying whole-mantle convection. In contrast, geochemical analyses of the basaltic products of mantle melting are frequently used to infer that mantle convection is layered, with the deeper mantle largely isolated from the upper mantle. We show that geochemical, seismological and heat-flow data are all consistent with whole-mantle convection provided that the observed heterogeneities are remnants of recycled oceanic and continental crust that make up about 16 and 0.3 per cent, respectively, of mantle volume.