Subduction zone stresses: Constraints from mechanics and from topographic and geoid anomalies

There are, in principle, direct relations between several important phenomena associated with subduction zones: the depth of oceanic trenches, the magnitude of the net force from trenches acting on subducting plates, the distribution and fault plane orientations of earthquakes, the magnitude of stresses on subduction faults, the bathymetry of back-arc regions, and the magnitudes of gravity and geoid anomalies. These phenomena are related through the stresses transmitted through surface and subducted lithosphere, and are associated with the mass anomaly of the subducted lithosphere. Quantitative estimates suggest that observed trench depths imply a trench pull force on subducting plates which is comparable to the ridge push force but much less than the excess weight of the subducted lithospheric slab. It is further suggested that either the mass anomaly of subducted lithosphere is much less than would be expected on the basis of conventional thermal and compositional models or that (a) a large resistance acts on the upper part of slabs due to high-stress corner flow, and (b) the mass anomaly of the slab is 70–90% compensated either by a broad 1 km-deep back-arc depression or a low density mantle wedge above the slab or both.     

The origin of reversed geochemical zoning in the northern New Hebrides volcanic arc

Integration of petrographic and geochemical data on each island in the Banks Group, northern Vanuatu (New Hebrides) has revealed a decrease in K₂O (and related incompatible elements) across the islands, away from the New Hebrides trench. The correlation with depth to the Benioff zone is therefore the reverse of that typically found in island arcs. REE evidence and source modelling indicates that this variation represents a progressive depletion in LREE/HREE in the upper mantle, laterally away from the trench and a progressive increase in partial melting of the source, in the same direction. These variations in chemistry are attributed to an earlier west-dipping Miocene subduction system, and the variation in degree of partial melting to the location of the islands with respect to the active back-arc basin. An enrichment of the upper mantle in the Miocene is ascribed to the ascent of hydrous fluids enriched in incompatible elements, from the west-dipping Benioff zone, and subsequent reaction with the overlying upper mantle. The model thus envisages approximately contemporaneous development of the Central Chain volcanoes and the backarc basin in the late Pliocene, with partial melting of hydrous, laterally zoned upper mantle resulting from the convective thermal regime associated with the rifting apart of the back-arc basin. The role of the present subduction regime in magma generation is considered to be very limited.     

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.     

西北太平洋边缘海热流特征研究

西北太平洋各边缘海及其相应俯冲系统受深部构造活动等地质条件的控制,热流变化较大。在收集整理该区域最新的热流数据基础上,重点探讨西北太平洋俯冲带热结构相关理论、边缘海大洋岩石圈热演化理论模型和局部高异常热流的影响因素,总结了西北太平洋边缘海热流所反映的地质意义。研究结果表明,在西北太平洋“沟-弧-海”体系中,从“沟”到“弧”再到“边缘海”,热流密度呈“低-高-较高”的变化趋势,弧后地区整体表现为“均一高热”特征;千岛海沟、日本海沟和琉球海沟热流密度值在30.0 mW/m²左右,而对应的岛弧值为其2~3倍。弧后热流大小受到汇聚型俯冲带热结构的影响,俯冲带脱水作用导致的弧后上地幔黏度变化,地震速度降低,岩石圈弹性厚度减薄,引起小尺度地幔对流,形成弧后“均一高热”的热状态。热流的时空分布与岩石圈年龄也有关,随着岩石圈年龄增大,地表热流密度值会随之降低,热流密度值大小和离散性与其形成时间大致呈负相关。鄂霍次克海形成时代(30~65 Ma)较早,其热流密度值(86.8 mW/m²)和离散性(标准差3.727)相对较低;冲绳海槽目前还处于扩张阶段,其热流密度值(139.0 mW/m²)和离散性(标准差7.001)较高。浅层的地下水循环、断裂活动,深层的地幔部分熔融岩浆活动、弧后小尺度地幔对流、俯冲带拐角流等对局部异常热流起到一定程度的控制作用。     

Seismic anisotropy in the region of the Chile margin triple junction

We have measured shear wave splitting at three temporary three-component short period stations that were deployed in southern Chile above the subducted Chile Rise spreading centre (Taitao Peninsula and environs). Subduction of the Chile Rise has been occurring beneath South America for at least the past 14 m.y. Previously published models of the ridge subduction posit the existence of ‘slab windows’, asthenosphere-filled gaps between subducted lithosphere segments of the spreading ridge, through which mantle might flow. Our preliminary results include two consistent fast polarization directions of splitting in the study region. Delay times between fast and slow split shear waves average around 1.0 s for all phases (ScS, PcS, SKS, and SKKS) that we measured. Fast-axis azimuths vary systematically among the three stations: near the coast, fast axes are parallel to the spreading ridge segments of the Chile Rise (approximately N-trending). This splitting fast-axis direction probably reflects either along-axis asthenospheric flow or results from the preferential attenuation effects of aligned pockets of melt at the subducted ridge segment. At one inland station above the slab window, we find two splitting fast-axis directions, one parallel to the subducted Chile Rise ridge segments, and a second trending NW–SE. We infer that upper mantle deformation in the vicinity of a well developed slab window is complicated and probably involves two superposed directions of upper mantle deformation. One of these directions (NW–SE) may indicate anomalous flow of asthenospheric mantle in the vicinity of the slab window gap.     

Two-dimensional simulations of surface deformation caused by slab detachment

Detachment of the deeper part of subducted lithosphere causes changes in a subduction zone system which may be observed on the Earth's surface. Constraints on the expected magnitudes of these surface effects can aid in the interpretation of geological observations near convergent plate margins where detachment is expected. In this study, we quantify surface deformation caused by detachment of subducted lithosphere. We determine the range of displacement magnitudes which can be associated with slab detachment using numerical models. The lithospheric plates in our models have an effective elastic thickness, which provides an upper bound for rapid processes, like slab detachment, to the surface deformation of lithosphere with a more realistic rheology. The surface topography which develops during subduction is compared with the topography shortly after detachment is imposed. Subduction with a non-migrating trench system followed by detachment leads to a maximum surface uplift of 2–6 km, while this may be higher for the case of roll-back preceding detachment. In the latter situation, the back-arc basin may experience a phase of compression after detachment. Within the context of our elastic model, the surface uplift resulting from slab detachment is sensitive to the depth of detachment, a change in friction on the subduction fault during detachment and viscous stresses generated by sinking of the detached part of the slab. Overall, surface uplift of these magnitudes is not diagnostic of slab detachment since variations during ongoing subduction may result in similar vertical surface displacements.     

Architecture and composition of ocean floor subducted beneath northern Gondwana during Neoproterozoic to Cambrian: A palinspastic reconstruction based on Ocean Plate Stratigraphy (OPS)

The Blovice accretionary complex, Bohemian Massif, hosts well-preserved basaltic blocks derived from an oceanic plate subducted beneath the northern active margin of Gondwana during late Neoproterozoic to early Cambrian. The major and trace element and Hf–Nd isotope systematics revealed two different suites, tholeiitic and alkaline, whose composition reflects different sources of melts within a back-arc basin setting. The former suite has composition similar to mid-ocean ridge basalts (MORB), yet with striking enrichment in large-ion lithophile elements (LILE) and Pb paralleled by depletion in Nb, in agreement with its derivation from depleted mantle fluxed by subduction-related fluids. In contrast, the latter suite has composition similar to ocean island basalts (OIB) with variable contribution of ancient, recycled crustal material. We argue that both suites represent volcanic members of Ocean Plate Stratigraphy (OPS) and indicate that the oceanic realm consumed by the Cadomian subduction was a complex mosaic of intra-oceanic subduction zones, volcanic island arcs, and back-arc basins with mantle plume impinging the spreading centre. Hence, the basalt geochemistry implies that two distinct domains of oceanic lithosphere may have existed off the Gondwana’s continental edge: an outboard domain, made up of old and less buoyant oceanic lithosphere (remnants of the Mirovoi Ocean surrounding former Rodinia?) that was steeply subducted and generated the back-arcs, and young, hot, and more buoyant oceanic lithosphere generated in the back-arcs and later involved in accretionary complexes as dismembered OPS. Perhaps the best recent analogy of this setting is the Izu Bonin–Mariana arc–Philippine Sea in the western Pacific.     

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.     

Volcanic arcs as archives of plate tectonic change

Processing of the oceanic lithosphere in subduction zones gives rise to arc magmatism, and strong compositional links exist between trench input and arc output. Here we address the question whether these compositional links are sufficiently strong to allow for ‘tracing’ the composition of the sedimentary and igneous oceanic crust through the chemistry of arcs. The tracing approach hinges critically on whether key characteristics of the subducted slab are transmitted to arcs. Results from forward and inverse modeling, verified by observations from modern arc settings, demonstrate that elements Sr, Pb, Nd and Hf that are associated with radiogenic isotopes may preserve chemical characteristics of the subducted slab in arc magmas. The data indicate that the much thicker igneous subducted crust dominates the recycled flux to arcs. The flux from the highly enriched, but thin sediment layer is buffered, and may be even concealed, by the concomitant contributions from igneous crust, and/or subarc mantle, despite the much better visibility of sediment components in trace element and isotope space. Arc Pb and Pb isotopes are the most promising tracers that may capture the isotopic diversity of subducted MORB-type and OIB-type crust with sufficient temporal and spatial resolution. While arc Sr is also strongly controlled by the flux from the subducted crust, arc data may allow for distinguishing among radiogenic Sr recycled from altered oceanic crust or from subducted sediment in moderately radiogenic arcs (⁸⁷Sr/⁸⁶Sr < ~ 0.7045). Co-mingling of Nd and Hf from igneous subducted crust with mantle contributions mostly hinders the isotopic identification of subducted crust through arc chemistry. However, Nd and Hf may provide complementary information about the efficiency of recycling, and recycling via subduction erosion.The tracing approach appears feasible in Cenozoic arcs where much of the original subduction context is preserved. First results from the Izu Bonin and Central American arcs show that plate tectonic events like oceanic plate formation and destruction, subduction of hotspot tracks and the closure of oceanic gateways are recorded in the chemistry of arcs. A comparative evaluation of Cenozoic global arcs may hence significantly complement the information from the modern oceanic basins, help to obtain a more complete image of the oceanic crustal composition and implicate the geochemical processes by which it formed. Possibly, the tracing approach may also be useful in ancient, inactive arcs to obtain information on the composition of oceanic crust subducted in the geological past.     

Geochemical response of magmas to Neogene–Quaternary continental collision in the Carpathian–Pannonian region: A review

The Carpathian–Pannonian Region contains Neogene to Quaternary magmatic rocks of highly diverse composition (calc-alkaline, shoshonitic and mafic alkalic) that were generated in response to complex microplate tectonics including subduction followed by roll-back, collision, subducted slab break-off, rotations and extension. Major element, trace element and isotopic geochemical data of representative parental lavas and mantle xenoliths suggests that subduction components were preserved in the mantle following the cessation of subduction, and were reactivated by asthenosphere uprise via subduction roll-back, slab detachment, slab-break-off or slab-tearing. Changes in the composition of the mantle through time are evident in the geochemistry, supporting established geodynamic models.Magmatism occurred in a back-arc setting in the Western Carpathians and Pannonian Basin (Western Segment), producing felsic volcaniclastic rocks between 21 to 18 Ma ago, followed by younger felsic and intermediate calc-alkaline lavas (18–8 Ma) and finished with alkalic-mafic basaltic volcanism (10–0.1 Ma). Volcanic rocks become younger in this segment towards the north. Geochemical data for the felsic and calc-alkaline rocks suggest a decrease in the subduction component through time and a change in source from a crustal one, through a mixed crustal/mantle source to a mantle source. Block rotation, subducted roll-back and continental collision triggered partial melting by either delamination and/or asthenosphere upwelling that also generated the younger alkalic-mafic magmatism.In the westernmost East Carpathians (Central Segment) calc-alkaline volcanism was simultaneously spread across ca. 100 km in several lineaments, parallel or perpendicular to the plane of continental collision, from 15 to 9 Ma. Geochemical studies indicate a heterogeneous mantle toward the back-arc with a larger degree of fluid-induced metasomatism, source enrichment and assimilation on moving north-eastward toward the presumed trench. Subduction-related roll-back may have triggered melting, although there may have been a role for back-arc extension and asthenosphere uprise related to slab break-off.Calc-alkaline and adakite-like magmas were erupted in the Apuseni Mountains volcanic area (Interior Segment) from15–9 Ma, without any apparent relationship with the coeval roll-back processes in the front of the orogen. Magmatic activity ended with OIB-like alkali basaltic (2.5 Ma) and shoshonitic magmatism (1.6 Ma). Lithosphere breakup may have been an important process during extreme block rotations (60°) between 14 and 12 Ma, leading to decompressional melting of the lithospheric and asthenospheric sources. Eruption of alkali basalts suggests decompressional melting of an OIB-source asthenosphere. Mixing of asthenospheric melts with melts from the metasomatized lithosphere along an east–west reactivated fault-system could be responsible for the generation of shoshonitic magmas during transtension and attenuation of the lithosphere.Voluminous calc-alkaline magmatism occurred in the Cãlimani-Gurghiu-Harghita volcanic area (South-eastern Segment) between 10 and 3.5 Ma. Activity continued south-eastwards into the South Harghita area, in which activity started (ca. 3.0–0.03 Ma, with contemporaneous eruption of calc-alkaline (some with adakite-like characteristics), shoshonitic and alkali basaltic magmas from 2 to 0.3 Ma. Along arc magma generation was related to progressive break-off of the subducted slab and asthenosphere uprise. For South Harghita, decompressional melting of an OIB-like asthenospheric mantle (producing alkali basalt magmas) coupled with fluid-dominated melting close to the subducted slab (generating adakite-like magmas) and mixing between slab-derived melts and asthenospheric melts (generating shoshonites) is suggested. Break-off and tearing of the subducted slab at shallow levels required explaining this situation.     

Vorticity Analysis and Deformation History of the Mizil Gneiss Dome,Eastern Arabian Shield,Saudi Arabia

The Mizil gneiss dome is an elliptical structure consisting of an amphibolite-facies volcanosedimentary mantle and a gneissic granite core. This dome is located at the northern tip of the Ar Rayn terrane only a few kilometers from the eastern edge of the Arabian shield. Previous investigations have shown the intrusive core to be an adakitic diapir with a U–Pb zircon age of 689 ± 10 Ma; this age is 50–80 Ma years older than other granites in this terrane. Vorticity analysis was carried out on samples from the intrusive core and volcanosedimentary cover; the Passchier and Rigid Grain Net (RGN) methods were used to obtain the kinematic vorticity number (W_k) and the mean kinematic vorticity number (W_m). The W_k and W_m values show a marked increase towards the south; such a pattern indicates a N-S movement of the core pluton thus creating an inclined diapir tilted to the south. Analogue experiments simulating the flow of magma diapirs rising form a subducted slab through the mantle wedge have shown that supra-subduction zone oblique diapirs are produced close to the trench and are elongated normal to the convergence direction as is the case in the Mizil pluton. This effect was found to increase with increasing slab dip due to enhanced drag along the upper surface of the subducted lithospheric plate. Spontaneous subduction which is often associated with rollback resulting in back-arc extension and steep dipping slabs is thought to have occurred in the Mozambique Ocean by 700 Ma. The Mizil pluton is coeval with the back-arc Urd ophiolite from the adjacent Dawadimi terrane, and could therefore have been produced by incipient subduction of a relatively cold slab as observed in many Pacific margin adakites. The tectonic evolution of the eastern shield, as deduced from the Mizil dome and other data from Ar Rayn and neighboring terranes, begins with the subduction of >100 My-old lithosphere beneath the Afif terrane resulting in back-arc spreading and the splitting of the Ar Rayn arc from the Afif microplate, with the concomitant production of a small volume of adakite melt. Other arc terrane(s) docked east of Ar Rayn with the westward-directed subduction still going but a lower angles and greater depth due to trench jump; this phase produced the more prevalent non-adakitic group-1 granites. A major collisional orogeny affected the entire eastern shield between 620–600 Ma and sutured the eastern shield terranes with northern Gondwana.     

A proposal for the kinematic modelling of W-dipping subductions - possible applications to the Tyrrhenian-Apennines system

In W-dipping subduction zones there is a general eastward progression of the back-arc basin-accretionary wedge-foredeep complex. With the forward progression, early stages of the complex are revealed by slices of upper crust and sedimentary cover abandoned to the west left floating above a new section of mantle. A major shear zone should form at the new Moho separating upper crust slices of earlier accretionary stages and the eastward flowing mantle. The mantle wedging at the top of the subduction plane could be responsible for the uplift of the central parts of the belt. The retreating of the subduction hinge is interpreted as due to the push generated by the ‘eastward mantle flow detected in the hot spot reference frame. The foredeep depth is mainly a function of the radius of curvature of the subduction hinge. The frontal wedge is constructed from the stacking of the upper layers of the subducting plate and the syntectonic clastics that fill the foredeep which are progressively involved in thrusting and later by extension. In order to preserve volume balance, the lithosphere of the eastern plate before subduction has to be the same size as that which has been subducted: due to the longer length of the arc with respect to the original length of the linear margin between the two converging plates, laterally stretched subducted lithosphere is predicted at depth. W-dipping subductions usually have a short life probably due to their inherent capability to produce new lateral heterogeneities of the lithosphere (the thin back-arc) which are a key factor in controlling and generating new subductions (both E- and W-dipping). This model is applied to the Apennines-Tyrrhenian Sea system.     

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.     

Subduction of a fossil slow–ultraslow spreading ocean: a petrology-constrained geodynamic model based on the Voltri Massif,Ligurian Alps,Northwest Italy

Slow–ultraslow spreading oceans are mostly floored by mantle peridotites and are typified by rifted continental margins, where subcontinental lithospheric mantle is preserved. Structural and petrologic investigations of the high-pressure (HP) Alpine Voltri Massif ophiolites, which were derived from the Late Jurassic Ligurian Tethys fossil slow–ultraslow spreading ocean, reveal the fate of the oceanic peridotites/serpentinites during subduction to depths involving eclogite-facies conditions, followed by exhumation.

The Ligurian Tethys was formed by continental extension within the Europe–Adria lithosphere and consisted of sea-floor exposed mantle peridotites with an uppermost layer of oceanic serpentinites and of subcontinental lithospheric mantle at the rifted continental margins. Plate convergence caused eastward subduction of the oceanic lithosphere of the Europe plate and the uppermost serpentinite layer of the subducting slab formed an antigorite serpentinite-subduction channel. Sectors of the rather unaltered mantle lithosphere of the Adria extended margin underwent ablative subduction and were detached, embedded, and buried to eclogite-facies conditions within the serpentinite-subduction channel. At such P–T conditions, antigorite serpentinites from the oceanic slab underwent partial HP dehydration (antigorite dewatering and growth of new olivine). Water fluxing from partial dehydration of host serpentinites caused partial HP hydration (growth of Ti-clinohumite and antigorite) of the subducted Adria margin peridotites. The serpentinite-subduction channel (future Beigua serpentinites), acting as a low-viscosity carrier for high-density subducted rocks, allowed rapid exhumation of the almost unaltered Adria peridotites (future Erro–Tobbio peridotites) and their emplacement into the Voltri Massif orogenic edifice. Over in the past 35 years, this unique geologic architecture has allowed us to investigate the pristine structural and compositional mantle features of the subcontinental Erro–Tobbio peridotites and to clarify the main steps of the pre-oceanic extensional, tectonic–magmatic history of the Europe–Adria asthenosphere–lithosphere system, which led to the formation of the Ligurian Tethys.

Our present knowledge of the Voltri Massif provides fundamental information for enhanced understanding, from a mantle perspective, of formation, subduction, and exhumation of oceanic and marginal lithosphere of slow–ultraslow spreading oceans.     

Asymmetric dynamics at subduction zones derived from plate kinematic constraints

The lithospheric sinking along subduction zones is part of the mantle convection. Therefore, computing the volume of lithosphere recycled within the mantle by subducting slabs quantifies the equivalent amount of mantle that should be displaced, for the mass conservation criterion. The rate of subduction is constrained by the convergence rate between upper and lower plates and the motion of the subduction hinge H that may either converge or diverge relative to the upper plate. Here, starting from the analysis of the slab hinge kinematics, we evaluate the subduction rate at 31 subduction zones worldwide, useful to compute volumes of sinking lithosphere into the mantle. Our results show that ∼190 km³/yr and ∼88 km³/yr of lithospheric slabs are currently subducting below H-divergent and H-convergent subduction zones, respectively. We also propose supporting numerical models providing asymmetric volumes of the subducted lithosphere, using the subduction rate instead of plate convergence, as boundary condition. Furthermore, H-divergent subduction zones appear to be coincident with subductions having “westward”-directed slabs, whereas H-convergent subduction zones are mostly compatible with those that have “eastward-to-northeastward”-directed slabs. On the basis of this geographical polarity, our lithospheric volume estimation gives ∼214 km³/yr and ∼88 km³/yr of subducting lithosphere, respectively. This entails that W-directed subduction zones contribute more than twice in lithospheric sinking into the mantle with respect to E-to-NE-directed ones. In accordance with the conservation of mass principle, this volumetric asymmetry in the mantle suggests a displacement of ∼120 km³/yr of mantle material from west to east, providing a constraint for global asymmetric mantle convection.     

Depleted arc volcanism in the Alboran Sea and shoshonitic volcanism in Morocco: geochemical and isotopic constraints on Neogene tectonic processes

Samples of volcanic rocks from Alborán Island, the Alboran Sea floor and from the Gourougou volcanic centre in northern Morocco have been analyzed for major and trace elements and Sr–Nd isotopes to test current theories on the tectonic geodynamic evolution of the Alboran Sea. The Alborán Island samples are low-K tholeiitic basaltic andesites whose depleted contents of HFS elements (0.5×N-MORB), especially Nb (0.2×N-MORB), show marked geochemical parallels with volcanics from immature intra-oceanic arcs and back-arc basins. Several of the submarine samples have similar compositions, one showing low-Ca boninite affinity. ¹⁴³Nd/¹⁴⁴Nd ratios fall in the same range as many island-arc and back-arc basin samples, whereas ⁸⁷Sr/⁸⁶Sr ratios (on leached samples) are somewhat more radiogenic. Our data point to active subduction taking place beneath the Alboran region in Miocene times, and imply the presence of an associated back-arc spreading centre. Our sea floor suite includes a few more evolved dacite and rhyolite samples with (⁸⁷Sr/⁸⁶Sr)₀ up to 0.717 that probably represent varying degrees of crustal melting. The shoshonite and high-K basaltic andesite lavas from Gourougou have comparable normalized incompatible-element enrichment diagrams and Ce/Y ratios to shoshonitic volcanics from oceanic island arcs, though they have less pronounced Nb deficits. They are much less LIL- and LREE-enriched than continental arc analogues and post-collisional shoshonites from Tibet. The magmas probably originated by melting in subcontinental lithospheric mantle that had experienced negligible subduction input. Sr–Nd isotope compositions point to significant crustal contamination which appears to account for the small Nb anomalies.

The unmistakable supra-subduction zone (SSZ) signature shown by our Alboran basalts and basaltic andesite samples refutes geodynamic models that attribute all Neogene volcanism in the Alboran domain to decompression melting of upwelling asthenosphere arising from convective thinning of over-thickened lithosphere. Our data support recent models in which subsidence is caused by westward rollback of an eastward-dipping subduction zone beneath the westernmost Mediterranean. Moreover, severance of the lithosphere at the edges of the rolling-back slab provides opportunities for locally melting lithospheric mantle, providing a possible explanation for the shoshonitic volcanism seen in northern Morocco and more sporadically in SE Spain.     

Dynamic subsidence of Eastern Australia during the Cretaceous

During the Early Cretaceous Australia's eastward passage over sinking subducted slabs induced widespread dynamic subsidence and formation of a large epeiric sea in the eastern interior. Despite evidence for convergence between Australia and the paleo-Pacific, the subduction zone location has been poorly constrained. Using coupled plate tectonic–mantle convection models, we test two end-member scenarios, one with subduction directly east of Australia's reconstructed continental margin, and a second with subduction translated ~ 1000 km east, implying the existence of a back-arc basin. Our models incorporate a rheological model for the mantle and lithosphere, plate motions since 140 Ma and evolving plate boundaries. While mantle rheology affects the magnitude of surface vertical motions, timing of uplift and subsidence depends on plate boundary geometries and kinematics. Computations with a proximal subduction zone result in accelerated basin subsidence occurring 20 Myr too early compared with tectonic subsidence calculated from well data. This timing offset is reconciled when subduction is shifted eastward. Comparisons between seismic tomography and model temperature cross-sections, and an absence of subduction zone volcanism in eastern Australia in the Early Cretaceous provide support for the back-arc basin scenario.     

Ocean, continents and orogens

In this paper, the processes related to subduction and mountain building are discussed, and some new models and notions are proposed.At all known epochs, the Earth's surface comprised essentially migrating plates and large belts (of the order of 10,000 × 2000 km) where the lithosphere is mobilised, so that subductions and crustal resorption occur in complex structural patterns. Through time, these orogens start as island arcs and evolve into folded ranges.The formation and location of island-arc belts is shown to be related to the obliquity between rifts and continental margins i.e. to lateral relative movements of plates: torsion couples are developed bringing arcs to bulge through the inducted arc mechanism. This accounts for tensions prevailing in the back-arc basin while compressive stresses accumulate as vertical deformations in the arc-trench or uplift—subsidence couple. When the rupture limit is reached, a tangential tectonic phase occurs.It is suggested that the energy output (heat flow) varies with the pressure exerted by the top (elastic) lithosphere upon the underlying mantle. A compensation tends to be established between the inner flow and tectonic stress, as the latter brings about pressure variations through uplift-subsidence couples. These may therefore be related to the generation of paired metamorphic belts.The evolution from island arcs to completed folded ranges is briefly described, with special attention being paid to ophiolitic and crystalline basement nappes, as well as to strike-slip faults, which are the final expression of lateral relative movements.     

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.     

热与克拉通破坏

大陆克拉通是地球表面上相对稳定的构造单元。从地热学的角度考虑,克拉通岩石圈的稳定意味着地表热流等于对流地幔岩石圈底部提供的热流加上岩石圈内部由放射性衰变产生的热量。太古代稳定克拉通一般具有冷的地热特征,且处于热平衡状态。打破克拉通热平衡的因素有多种,如岩石圈地幔的放射性生热、来自深部的地幔柱、板块俯冲等。华北是全球古老克拉通遭受破坏最明显和最典型的地区,其破坏与中生代太平洋板块向东亚的深俯冲密切相关。古太平洋板块快速俯冲并停滞在地幔转换带脱水、形成宽约1 000 km的低粘大地幔楔,导致地幔对流增强。在活跃地幔对流的热侵蚀与橄榄岩-熔体相互作用共同作用下,华北克拉通在中生代期间迅速减薄。经过中生代加热减薄的华北克拉通岩石圈强度显著变弱,在俯冲板块后撤作用下,岩石圈拉张并进一步减薄,地表热流升高。华北克拉通破坏是一个漫长的历程,期间大地热流的演化特征呈现出由中生代以前的低值演化至新生代的高值再过渡到现今的中等状态（接近全球大陆平均值）。