Complex slab structure and arc magmatism beneath the Japanese Islands

A dense nationwide seismic network recently constructed in Japan has resulted in the production of a large amount of high-quality data that have enabled the high-resolution imaging of deep seismic structures in the Japanese subduction zone. Seismic tomography, precise locations of earthquakes, and focal mechanism research have allowed the identification of the complex structure of subducting slabs beneath Japan, revealing that the subducting Philippine Sea slab underneath southwestern Japan has an undulatory configuration down to a depth of 60–200 km, and is continuous from Kanto to Kyushu without disruption or splitting, even within areas north of the Izu Peninsula. Analysis of the geometry of the Pacific and Philippine Sea slabs identified a broad contact zone beneath the Kanto Plain that causes anomalously deep interplate and intraslab earthquake activity. Seismic tomographic inversions using both teleseismic and local events provide a clear image of the deep aseismic portion of the Philippine Sea slab beneath the Japan Sea north of Chugoku and Kyushu, and beneath the East China Sea west of Kyushu down to a depth of ∼450 km. Seismic tomography also allowed the identification of an inclined sheet-like seismic low-velocity zone in the mantle wedge beneath Tohoku. A recent seismic tomography work further revealed clear images of similar inclined low-velocity zones in the mantle wedge for almost all other areas of Japan. The presence of the inclined low-velocity zones in the mantle wedge across the entirety of Japan suggests that it is a common feature to all subduction zones. These low-velocity zones may correspond to the upwelling flow portion of subduction-induced convection systems. These upwelling flows reach the Moho directly beneath active volcanic areas, suggesting a link between volcanism and upwelling.     

Aseismic deep subduction of the Philippine Sea plate and slab window

We have made great efforts to collect and combine a large number of high-quality data from local earthquakes and teleseismic events recorded by the dense seismic networks in both South Korea and West Japan. This is the first time that a large number of Korean and Japanese seismic data sets are analyzed jointly. As a result, a high-resolution 3-D P-wave velocity model down to 700-km depth is determined, which clearly shows that the Philippine Sea (PHS) plate has subducted aseismically down to ∼460 km depth under the Japan Sea, Tsushima Strait and East China Sea. The aseismic PHS slab is visible in two areas: one is under the Japan Sea off western Honshu, and the other is under East China Sea off western Kyushu. However, the aseismic PHS slab is not visible between the two areas, where a slab window has formed. The slab window is located beneath the center of the present study region where many teleseismic rays crisscross. Detailed synthetic tests were conducted, which indicate that both the aseismic PHS slab and the slab window are robust features. Using the teleseismic data recorded by the Japanese stations alone, the aseismic PHS slab and the slab window were also revealed (Zhao et al., 2012), though the ray paths in the Japanese data set crisscross less well offshore. The slab window may be caused by the subducted Kyushu-Palau Ridge and Kinan Seamount Chain where the PHS slab may be segmented. Hot mantle upwelling is revealed in the big mantle wedge above the Pacific slab under the present study region, which may have facilitated the formation of the PHS slab window. These novel findings may shed new light on the subduction history of the PHS plate and the dynamic evolution of the Japan subduction zone.     

Amount of Asian lithospheric mantle subducted during the India/Asia collision

Body wave seismic tomography is a successful technique for mapping lithospheric material sinking into the mantle. Focusing on the India/Asia collision zone, we postulate the existence of several Asian continental slabs, based on seismic global tomography. We observe a lower mantle positive anomaly between 1100 and 900 km depths, that we interpret as the signature of a past subduction process of Asian lithosphere, based on the anomaly position relative to positive anomalies related to Indian continental slab. We propose that this anomaly provides evidence for south dipping subduction of North Tibet lithospheric mantle, occurring along 3000 km parallel to the Southern Asian margin, and beginning soon after the 45 Ma break-off that detached the Tethys oceanic slab from the Indian continent. We estimate the maximum length of the slab related to the anomaly to be 400 km. Adding 200 km of presently Asian subducting slab beneath Central Tibet, the amount of Asian lithospheric mantle absorbed by continental subduction during the collision is at most 600 km. Using global seismic tomography to resolve the geometry of Asian continent at the onset of collision, we estimate that the convergence absorbed by Asia during the indentation process is ~ 1300 km. We conclude that Asian continental subduction could accommodate at most 45% of the Asian convergence. The rest of the convergence could have been accommodated by a combination of extrusion and shallow subduction/underthrusting processes. Continental subduction is therefore a major lithospheric process involved in intraplate tectonics of a supercontinent like Eurasia.     

Origin of Late Cenozoic Abaga–Dalinuoer basalts,eastern China: Implications for a mixed pyroxenite–peridotite source related with deep subduction of the Pacific slab

Continental intraplate basalts (15.42–0.16 Ma) from Abaga–Dalinuoer volcanic field (ADVF) in central Inner Mongolia of eastern China, as a part of Cenozoic volcanic province along eastern margin of the Eurasian continent, provide a good opportunity to explore potential links between deep subduction of the Pacific slab and continental intraplate volcanism. In this study, we report an integrated dataset of whole-rock K–Ar ages, major and trace elements and Sr–Nd–Pb isotopes, and olivine major and minor elements for the Abaga–Dalinuoer basalts (ADBs), and propose that mantle source lithology of the ADB magmas may consist of both pyroxenite and peridotite. The ADBs display low SiO₂ (42.3–50.2 wt.%), high MgO (7.3–11.4 wt.%) and moderate K₂O + Na₂O (3.8–6.4 wt.%), and can be subdivided into basanites, alkali basalts and tholeiitic basalts that are all characterized by ocean island basalt (OIB)-like rare earth elements (REE) and enrichment in both large ion lithosphile elements (LILE) and high field strength elements (HFSE). Olivine phenocrysts have higher Ni and Fe/Mn and lower Mn, Ca and Ca/Fe relative to those from peridotite melts, but exhibit clearly lower Ni contents (< 2500 ppm) compared with expected Ni range (> 3000 ppm) for olivines crystallized from olivine-free pyroxenite melts. Estimated compositions of the ADB primary magmas, together with olivine compositions, suggest an iron-rich mantle source related with silica-deficient pyroxenite that is most likely derived from deeply subducted Pacific oceanic crust. Additionally, peridotite and recent subducted sediments are also required to account for high Ni and MgO in primary magmas together with their trace elements and Sr–Nd–Pb isotope systematics. We suggest that a mixed pyroxenite–peridotite source lithology can better match observed whole-rock and olivine signatures in the ADB, compared with either peridotite only or olivine-free pyroxenite only source lithology. In our model, pyroxenite melts would either react with or mechanically mix with peridotite, and the proportion of pyroxenite melts may range from 30% to 45% for mechanical mixing scenario. A continuous spectrum from tholeiitic to alkali melts revealed by melt-peridotite reaction experiment can explain formation of primary magmas of basanites, alkali basalts and tholeiitic basalts by increasing melting degree of a similar mantle source. Relatively higher ²⁰⁶Pb/²⁰⁴Pb of the ADB may suggest more significant role of recent (< 0.5 Ga) subducted Pacific oceanic materials, in contrast to other Cenozoic basalts in eastern China (e.g., Changbai basalts) that exhibit varying contributions from ancient (> 1.5 Ga) subducted continental sediments. We emphasize that coupled geochemical and geodynamic links (i.e., subduction polarity) between deeply subducted Pacific slab and continental intraplate volcanism in eastern China may exist, which directly support the involvement of deeply subducted Pacific materials in petrogenesis of the ADB. From the perspective of plate motion kinetics, decompression partial melting of upwelling fragmented Pacific slab (silica-deficient pyroxenite + recent subducted sediments) may be triggered by rollback of deeply subducted Pacific slab during Late Cenozoic times. Continental intraplate volcanism in the ADVF generally started with termination of opening of the Japan Sea, suggesting that deep subduction of the Pacific slab may have been an important geodynamic mechanism responsible for tectono-magmatic evolution of northeastern Asia. We suggest that the ADBs have the potential to shed light on genetic links between continental intraplate volcanism and deep subduction of the Pacific slab in geochemical and geodynamic processes.     

Geochemistry of Neoproterozoic mafic intrusions in the Panzhihua district (Sichuan Province,SW China): Implications for subduction-related metasomatism in the upper mantle

Metasomatism above subduction zones is an important process that produces heterogeneous mantle and thus a diversity of igneous rocks. In the Panzhihua district, on the western margin of the Yangtze Block (SW China), two Neoproterozoic mafic intrusions, one olivine gabbro and one hornblende gabbro, have identical ages of 746 ± 10 and 738 ± 23 Ma. Both of the gabbros are tholeiitic in composition and have arc-like geochemical compositions. The hornblende gabbros have K₂O concentrations ranging from 0.70 to 1.69 wt.% and show enrichment of Rb, Ba, U, Th and Pb and depletion of Nb,Ta and Ti. They have variable ⁸⁷Sr/⁸⁶Sr ratios (0.7045–0.7070) with constant ɛNd(t) values (−0.12 to −0.93). The olivine gabbros have relatively low K₂O (0.19–0.43 wt.%), are depleted in Rb and Th relative to Ba and U, and have obvious negative Nb–Ta and Zr–Hf anomalies on primitive mantle-normalized trace element diagrams. Their ɛNd(t) values range from −0.64 to −1.73 and initial ⁸⁷Sr/⁸⁶Sr ratios from 0.7070 to 0.7075. Both types of gabbro experienced fractional crystallization of clinopyroxene, plagioclase, amphibole and minor Fe–Ti oxide. The parental magmas of the olivine and hornblende gabbros were formed by about 20% partial melting of garnet–spinel lherzolite and spinel lherzolite, respectively. According to trace elemental ratios, the hornblende gabbros were probably derived from a source strongly modified by subducted slab fluids, whereas the olivine gabbros came from a mantle source modified by subducted slab melts. The close association of the olivine gabbros and hornblende gabbros suggests that a steep subduction zone existed along the western margin of the Yangtze Block during Neoproterozoic time. Thus, the giant Neoproterozoic magmatic event in South China was subduction-related.     

Thermochemical structure of the North China Craton from multi-observable probabilistic inversion: Extent and causes of cratonic lithosphere modification

We present a 3D thermochemical model of the North China Craton (NCC) from the surface down to 350 km by jointly inverting surface wave phase velocity data, geoid height, surface heat flow and absolute elevation with a multi-observable probabilistic inversion method. Our model reveals a thin (~ 65–100 km) and chemically fertile lithosphere (87 < Mg# < 90) beneath the Eastern NCC, consistent with independent results from mantle xenoliths, and supports the idea that the Eastern NCC experienced significant lithospheric destruction and refertilization during the Phanerozoic. In contrast, beneath the Trans-North China Orogen, Inner Mongolia Suture Zone and Yinshan belt, we observe a more heterogeneous (chemically and thermally) lithosphere, indicating that these areas have been partly involved in lithospheric modification and mechanical erosion at multiple scales. A cold and chemically refractory (Mg# > 90) lithospheric mantle is imaged beneath the central TNCO and Ordos Block, reaching depths > 260 km. This lithospheric “keel” is surrounded to the east by a high-temperature sublithospheric anomaly that originates at depths > 280 km. The spatial distribution of this anomaly and its correlation with the location of recent volcanism in the region suggest that the anomaly represents a deep mantle upwelling being diverted by the cratonic keel and spreading onto regions of shallow lithosphere. Our results indicate that the present-day thermochemical structure beneath the NCC is the result of a complex interaction between a large-scale return flow associated with the subduction of the Pacific slab and the shallow lithospheric structure.     

Subduction-related metasomatic mantle source in the eastern Central Asian Orogenic Belt: Evidence from amphibolites in the Xilingol Complex,Inner Mongolia,China

The Central Asian Orogenic Belt (CAOB) formed mainly in the Paleozoic due to the closure of the Paleo-Asian oceanic basins and accompanying prolonged accretion of pelagic sediments, oceanic crust, magmatic arcs, and Precambrian terranes. The timing of subduction–accretion processes and closure of the Paleo-Asian Ocean has long been controversial and is addressed in a geochemical and isotopic investigation of mafic rocks, which can yield important insight into the geodynamics of subduction zone environments. The Xilingol Complex, located on the northern subduction–accretion zone of the CAOB, mainly comprises strongly deformed quartzo-feldspathic gneisses with intercalated lenticular or quasi-lamellar amphibolite bodies. An integrated study of the petrology, geochemistry, and geochronology of a suite of amphibolites from the complex constrains the nature of the mantle source and the tectono-metamorphic events in the belt. The protoliths of these amphibolites are gabbros and gabbroic diorites that intruded at ca. 340–321 Ma with positive εHf_(t) values ranging from + 2.89 to + 12.98. Their T_DM1 model ages range from 455 to 855 Ma and peak at 617 Ma, suggesting that these mafic rocks are derived from a depleted continental lithospheric mantle. The primitive magma was generated by variable degrees of partial melting of spinel-bearing peridotites. Fractionation of olivine, clinopyroxene and hornblende has played a dominant role during magma differentiation with little or no crustal contamination. The mafic rocks are derived from a Late Neoproterozoic depleted mantle source that was subsequently enriched by melts affected by slab-derived fluids and sediments, or melts with a sedimentary source rock. The Carboniferous mafic rocks in the northern accretionary zone of the CAOB record a regional extensional event after the Early Paleozoic subduction of the Paleo-Asian Ocean. Both addition of mantle-derived magmas and recycling of oceanic crust played key roles in significant Late Carboniferous (ca. 340–309 Ma) vertical crustal growth in the CAOB. Amphibolite–facies metamorphism (P = 0.34–0.52 GPa, T = 675–708 °C) affected these mafic rocks in the Xilingol Complex at ca. 306–296 Ma, which may be related to the crustal thickening by northward subduction of a forearc oceanic crust beneath the southern margin of the South Mongolian microcontinent. The final formation of the Solonker zone may have lasted until ca. 228 Ma.     

A Late Cretaceous and Cenozoic reconstruction of the Southwest Pacific region: Tectonics controlled by subduction and slab rollback processes

A Cenozoic tectonic reconstruction is presented for the Southwest Pacific region located east of Australia. The reconstruction is constrained by large geological and geophysical datasets and recalculated rotation parameters for Pacific–Australia and Lord Howe Rise–Pacific relative plate motion. The reconstruction is based on a conceptual tectonic model in which the large-scale structures of the region are manifestations of slab rollback and backarc extension processes. The current paradigm proclaims that the southwestern Pacific plate boundary was a west-dipping subduction boundary only since the Middle Eocene. The new reconstruction provides kinematic evidence that this configuration was already established in the Late Cretaceous and Early Paleogene. From ∼ 82 to ∼ 52 Ma, subduction was primarily accomplished by east and northeast-directed rollback of the Pacific slab, accommodating opening of the New Caledonia, South Loyalty, Coral Sea and Pocklington backarc basins and partly accommodating spreading in the Tasman Sea. The total amount of east-directed rollback of the Pacific slab that took place from ∼ 82 Ma to ∼ 52 Ma is estimated to be at least 1200 km. A large percentage of this rollback accommodated opening of the South Loyalty Basin, a north–south trending backarc basin. It is estimated from kinematic and geological constraints that the east–west width of the basin was at least ∼ 750 km. The South Loyalty and Pocklington backarc basins were subducted in the Eocene to earliest Miocene along the newly formed New Caledonia and Pocklington subduction zones. This culminated in southwestward and southward obduction of ophiolites in New Caledonia, Northland and New Guinea in the latest Eocene to earliest Miocene. It is suggested that the formation of these new subduction zones was triggered by a change in Pacific–Australia relative motion at ∼ 50 Ma. Two additional phases of eastward rollback of the Pacific slab followed, one during opening of the South Fiji Basin and Norfolk Basin in the Oligocene to Early Miocene (up to ∼ 650 km of rollback), and one during opening of the Lau Basin in the latest Miocene to Present (up to ∼ 400 km of rollback). Two new subduction zones formed in the Miocene, the south-dipping Trobriand subduction zone along which the Solomon Sea backarc Basin subducted and the north-dipping New Britain–San Cristobal–New Hebrides subduction zone, along which the Solomon Sea backarc Basin subducted in the west and the North Loyalty–South Fiji backarc Basin and remnants of the South Loyalty–Santa Cruz backarc Basin subducted in the east. Clockwise rollback of the New Hebrides section resulted in formation of the North Fiji Basin. The reconstruction provides explanations for the formation of new subduction zones and for the initiation and termination of opening of the marginal basins by either initiation of subduction of buoyant lithosphere, a change in plate kinematics or slab–mantle interaction.     

Petrology and geochemistry of the Karaj Dam basement sill: Implications for geodynamic evolution of the Alborz magmatic belt

The northeastward subduction of the Neo-Tethyan oceanic lithosphere beneath the Iranian block produced vast volcanic and plutonic rocks that now outcrop in central (Urumieh–Dokhtar magmatic assemblage) and north–northeastern Iran (Alborz Magmatic Belt), with peak magmatism occurring during the Eocene. The Karaj Dam basement sill (KDBS), situated in the Alborz Magmatic Belt, comprises gabbro, monzogabbro, monzodiorite, and monzonite with a shoshonitic affinity. These plutonic rocks are intruded into the Karaj Formation, which comprise pyroclastic rocks dating to the lower–upper Eocene. The geochemical and isotopic signatures of the KDBS rocks indicate that they are cogenetic and evolved through fractional crystallization. They are characterized by an enrichment in LREEs relative to HREEs, with negative Nb–Ta anomalies. Geochemical modeling using Sm/Yb versus La/Yb and La/Sm ratios suggests a low-degree of partial melting of a phlogopite–spinel peridotite source to generate the KDBS rocks. Their low I_Sr = 0.70453–0.70535, ɛNd (37.2 Ma) = 1.54–1.9, and T_DM ages ranging from 0.65 to 0.86 Ga are consistent with the melting of a Cadomian enriched lithospheric mantle source, metasomatized by fluids derived from the subducted slab or sediments during magma generation. These interpretations are consistent with high ratios of ²⁰⁶Pb/²⁰⁴Pb = 18.43–18.67, ²⁰⁷Pb/²⁰⁴Pb = 15.59, and ²⁰⁸Pb/²⁰⁴Pb = 38.42–38.71, indicating the involvement of subducted sediments or continental crust. The sill is considered to have been emplaced in an environment of lithospheric extension due to the slab rollback in the lower Eocene. This extension led to localized upwelling of the asthenosphere, providing the heat required for partial melting of the subduction-contaminated subcontinental lithospheric mantle beneath the Alborz magmatic belt. Then, the shoshonitic melt generates the entire spectrum of KDBS rocks through assimilation and fractional crystallization during the ascent of the magma.     

Lithospheric structure beneath the East China Sea revealed by Rayleigh-wave phase velocities

We explore the variations of Rayleigh-wave phase-velocity beneath the East China Sea in a broad period range (5–200 s). Rayleigh-wave dispersion curves are measured by the two-station technique for a total of 373 interstation paths using vertical-component broad-band waveforms at 32 seismic stations around the East China Sea from 6891 global earthquakes.The resulting maps of Rayleigh-wave phase velocity and azimuthal anisotropy provide a high resolution model of the lithospheric mantle beneath the East China Sea. The model exhibits four regions with different isotropic and anisotropic patterns: the Bohai Sea, belonging to the North China Craton, displays a continental signature with fast velocities at short periods; the Yellow Sea, very stable unit associated with low deformation, exhibits fast velocities and limited anisotropy; the southern part of the East China Sea, with high deformation and many fractures and faults, is related to slow velocities and high anisotropic signature; and the Ryukyu Trench shows high-velocity perturbations and slab parallel anisotropy.     

Adakitic rocks derived from the partial melting of subducted continental crust: Evidence from the Eocene volcanic rocks in the northern Qiangtang block

Eocene is a critical time for the elevation of Tibetan Plateau and global climate change, and previous studies suggested that the Eocene elevation was caused by intra-continental subduction of the Songpan–Garze block beneath the Qiangtang block. This paper reports zircon U–Pb age and geochemistry of the Eocene volcanic rocks from the Zuerkenwula mountain area in the northern part of Qiangtang block, and proposes that both slab break-off of the Neo-Tethys oceanic slab along the Bangong–Nujiang suture and intra-continental subduction of the Songpan–Garze block beneath the Qiangtang block caused the extensive partial melting of lithospheric mantle and subducted Songpan–Garze continental crust, which resulted in the significant elevation of the Tibetan Plateau. The volcanic rocks have LA-ICP MS U–Pb zircon age of 40.25 ± 0.15 Ma (MSWD = 2.1, 2σ), which is contemporaneous with the Eocene eclogites in the Great Himalayan and K-rich lavas in the southeastern Tibet. They display some adakitic characteristics with SiO₂ = 57.44 to 68.72%, TiO₂ = 0.38 to 0.81%, Na₂O = 2.89 to 4.35%, K₂O = 2.77 to 4.48%, Al₂O₃ = 13.92 to 18.22%, A/CNK = 0.69 to 1.03, MgO = 0.27 to 5.86% with Mg# ranging from 13.2 to 72.0, strongly depleted in heavy rare earth elements (HREEs) (Yb = 0.92 to 1.51 ppm and Y = 10.1 to 24.1 ppm), in combination with their positive Sr anomalies, high Sr/Y ratios and no significant Eu anomalies, which suggest a garnet-in and plagioclase-free source residue. These volcanic rocks can be divided into high-Mg# (> 45) and low-Mg# (< 45) groups. Both of the two groups share evolved Sr–Nd–Pb isotopic compositions with ⁸⁷Sr/⁸⁶Sr = 0.707412–0.708284; ε_Nd(t) = ? 4.0 to ? 5.7; ²⁰⁶Pb/²⁰⁴Pb = 18.7499–18.8189, ²⁰⁷Pb/²⁰⁴Pb = 15.7189–15.7384; ²⁰⁸Pb/²⁰⁴Pb = 39.166–39.262. The geophysical data and regional geological setting suggest that the low-Mg# adakitic rocks were derived from the decompression melting of a subducted lower continental crust, when low-Mg# adakitic melts in the overlying peridotite mantle wedge captured some olivine crystals, resulting in their elevated Mg# and MgO values.     

Tectonic controls on post-subduction granite genesis and emplacement: The late Caledonian suite of Britain and Ireland

Rates of magma emplacement commonly vary as a function of tectonic setting. The late Caledonian granites of Britain and Ireland are associated with closure of the Iapetus Ocean and were emplaced into a varying regime of transpression and transtension throughout the Silurian and into the early Devonian. Here we evaluate a new approach for examining how magma volumes vary as a function of tectonic setting. Available radiometric ages from the late Caledonian granites are used to calculate probability density functions (age spectra), with each pluton weighted by outcrop area as a proxy for its volume. These spectra confirm an absence of magmatic activity during Iapetus subduction between c. 455 Ma and 425 Ma and a dominance of post-subduction magmas between c. 425 Ma and 380 Ma. We review possible reasons why, despite the widespread outcrop of the late Caledonian granites, magmatism appears absent during Iapetus subduction. These include shallow angle subduction or extensive erosion and tectonic removal of the arc.In contrast to previous work, we find no strong difference in the age or major element chemistry of post-subduction granites across all terranes. We propose a common causal mechanism in which the down-going Iapetus oceanic slab peeled back and detached beneath the suture following final Iapetus closure. The lithospheric mantle was delaminated beneath the suture and for about 100 km back beneath the Avalonian margin. While magma generation is largely a function of gravitationally driven lithosphere delamination, strike slip dominated kinematics in the overlying continental crust is what modulated granitic magma emplacement. Early Devonian (419–404 Ma) transtension permitted large volumes of granite emplacement, whereas the subsequent Acadian (late Early Devonian, 404–394 Ma) transpression reduced and eventually suppressed magma emplacement.     

Geophysical constraints on geodynamic processes at convergent margins: A global perspective

Convergent margins, being the boundaries between colliding lithospheric plates, form the most disastrous areas in the world due to intensive, strong seismicity and volcanism. We review global geophysical data in order to illustrate the effects of the plate tectonic processes at convergent margins on the crustal and upper mantle structure, seismicity, and geometry of subducting slab. We present global maps of free-air and Bouguer gravity anomalies, heat flow, seismicity, seismic Vs anomalies in the upper mantle, and plate convergence rate, as well as 20 profiles across different convergent margins. A global analysis of these data for three types of convergent margins, formed by ocean–ocean, ocean–continent, and continent–continent collisions, allows us to recognize the following patterns. (1) Plate convergence rate depends on the type of convergent margins and it is significantly larger when, at least, one of the plates is oceanic. However, the oldest oceanic plate in the Pacific ocean has the smallest convergence rate. (2) The presence of an oceanic plate is, in general, required for generation of high-magnitude (M > 8.0) earthquakes and for generating intermediate and deep seismicity along the convergent margins. When oceanic slabs subduct beneath a continent, a gap in the seismogenic zone exists at depths between ca. 250 km and 500 km. Given that the seismogenic zone terminates at ca. 200 km depth in case of continent–continent collision, we propose oceanic origin of subducting slabs beneath the Zagros, the Pamir, and the Vrancea zone. (3) Dip angle of the subducting slab in continent–ocean collision does not correlate neither with the age of subducting oceanic slab, nor with the convergence rate. For ocean–ocean subduction, clear trends are recognized: steeply dipping slabs are characteristic of young subducting plates and of oceanic plates with high convergence rate, with slab rotation towards a near-vertical dip angle at depths below ca. 500 km at very high convergence rate. (4) Local isostasy is not satisfied at the convergent margins as evidenced by strong free air gravity anomalies of positive and negative signs. However, near-isostatic equilibrium may exist in broad zones of distributed deformation such as Tibet. (5) No systematic patterns are recognized in heat flow data due to strong heterogeneity of measured values which are strongly affected by hydrothermal circulation, magmatic activity, crustal faulting, horizontal heat transfer, and also due to low number of heat flow measurements across many margins. (6) Low upper mantle Vs seismic velocities beneath the convergent margins are restricted to the upper 150 km and may be related to mantle wedge melting which is confined to shallow mantle levels.     

Subduction of a spreading ridge within the Bangong Co–Nujiang Tethys Ocean: Evidence from Early Cretaceous mafic dykes in the Duolong porphyry Cu–Au deposit,western Tibet

The intrusion of mafic dykes into a near-trench accretionary prism, and continental margin magmatism with characteristics that differ from those of adjacent arc magmatism, are direct manifestations of the subduction of a spreading ocean ridge and the formation of a slab window. In this paper, we investigated mafic dykes intruded into the accretionary prism that hosts the Duolong porphyry Cu–Au deposit (DCAD) of western Tibet. LA-ICP-MS analysis of U–Pb in zircon indicates that the dykes formed during the Early Cretaceous (126–127 Ma). The dykes are characterized by ε_Hf(t) values from + 2.44 to + 11.8. Twenty-nine mafic dyke samples were divided into three groups based on their locations and geochemical compositions: group I has Nb = 8.31–10.2 ppm, Nb/La = 0.71–1.20, and Nb/U = 21.4–37.9; group II has Nb = 40.5–52.6 ppm, Nb/La = 0.84–1.58, and Nb/U = 18.8–47.8; and group III has Nb = 65.7–105 ppm, Nb/La = 1.35–2.08, and Nb/U = 36.5–73.8. Group I is classified as Nb-enriched basalts (ENBs), whereas groups II and III are classified as high-Nb basalts (HNBs). Both the ENBs and HNBs were derived from an adakite-metasomatized mantle wedge that subsequently underwent crystallization of olivine and clinopyroxene. The compositional variations of the studied dykes resulted mainly from mantle source heterogeneity. The volume of the slab melts gradually increases from group I (ENBs) to group II (HNBs) and group III (HNBs), leading to gradually increasing incompatible element concentrations. Considering their geochemical characteristics and field relationships, as well as the unique characteristics of continental margin magmatism in the DCAD, we propose that the dykes emplaced in an extensional accretionary prism were derived from the northward subduction of a spreading ridge in the Bangong Co–Nujiang Tethys Ocean during the Early Cretaceous. The ridge subduction event was also responsible for the generation of coeval adakites, intermediate–felsic intrusions, the Maierze bimodal volcanic rocks as well as the metallogenesis of the DCAD.     

Regimes of subduction and lithospheric dynamics in the Precambrian: 3D thermomechanical modelling

Comparing the early Earth to the present day, geological–geochemical evidence points towards higher mantle potential temperature and a different type of tectonics. In order to investigate possible changes in Precambrian tectonic styles, we conduct 3D high-resolution petrological–thermomechanical numerical modelling experiments for oceanic plate subduction under an active continental margin at a wide range of mantle potential temperature T_P (∆ T_P = 0 − 250 K, compared to present day conditions). At present day mantle temperatures (∆ T_P = 0 K), results of numerical experiments correspond to modern-style subduction, whereas at higher temperature conditions important systematic changes in the styles of both lithospheric deformation and mantle convection occur. For ∆ T_P = 50 − 100 K a regime of dripping subduction emerges which is still very similar to present day subduction but is characterised by frequent dripping from the slab tip and a loss of coherence of the slab, which suggests a close relationship between dripping subduction and episodic subduction. At further increasing ∆ T_P = 150 − 200 K dripping subduction is observed together with unstable dripping lithosphere, which corresponds to a transitional regime. For ∆ T_P = 250 K, presumably equivalent to early Archean, the dominating tectonic style is characterised by small-scale mantle convection, unstable dripping lithosphere, thick basaltic crust and small plates. Even though the initial setup is still defined by present day subduction, this final regime shows many characteristics of plume-lid tectonics. Transition between the two end-members, plume-lid tectonics and plate tectonics, happens gradually and at intermediate temperatures elements of both tectonic regimes are present. We conclude, therefore, that most likely no abrupt geodynamic regime transition point can be specified in the Earth's history and its global geodynamic regime gradually evolved over time from plume-lid tectonics into modern style plate tectonics.     

Tectono-magmatic evolution of Late Jurassic to Early Cretaceous granitoids in the west central Lhasa subterrane,Tibet

There is ongoing debate as to the subduction direction of the Bangong–Nujiang Ocean during the Mesozoic (northward, southward or bidirectional subduction). Arc-related intermediate to felsic intrusions could mark the location of the subduction zone and, more importantly, elucidate the dominant geodynamic processes. We report whole rock geochemical and zircon U–Pb and Hf isotopic data for granitoids from the west central Lhasa subterrane (E80° to E86°). All rocks show metaluminous to peraluminous, calc-alkaline signatures, with strong depletion of Nb, Ta and Ti, enrichment of large ion lithophile elements (e.g., Cs, Rb, K), a negative correlation between SiO₂ and P₂O₅, and a positive correlation between Rb and Th. All these features are indicative of I-type arc magmatism. New zircon U–Pb results, together with data from the literature, indicate continuous magmatism from the Late Jurassic to the Early Cretaceous (160 to 130 Ma). Zircon U–Pb ages for samples from the northern part of the west central Lhasa subterrane (E80° to E82°30′) yielded formation ages of 165 to 150 Ma, whereas ages of 142 to 130 Ma were obtained on samples from the south. This suggests flat or low-angle subduction of the Bangong–Nujiang Ocean, consistent with a slight southward decrease in zircon εHf(t) values for Late Jurassic rocks. Considering the crustal shortening, the distance from the Bangong–Nujiang suture zone, and a typical subduction zone melting depth of ~ 100 km, the subduction angle was less than 14° for Late Jurassic magmatism in the central Lhasa interior, consistent with flat or low-angle subduction. Compared with Late Jurassic rocks (main εHf(t) values of − 16 to − 7), Early Cretaceous rocks (145 to 130 Ma) show markedly higher εHf(t) values (mainly − 8 to 0), possibly indicating slab roll-back, likely caused by slab foundering or break-off. Combined with previously published works on arc magmatism in the central Lhasa and west part of the southern Qiangtang subterranes, our results support the bidirectional subduction of the Bangong–Nujiang Ocean along the Bangong–Nujiang Suture Zone, and indicates flat or low-angle southward subduction (165 to 145 Ma) followed by slab roll-back (145 to 130 Ma).     

Hot lherzolite exhumation,UHT migmatite formation,and acid volcanism driven by Miocene rollback of the Banda Arc,eastern Indonesia

The northern Banda Arc, eastern Indonesia, exposes upper mantle/lower crustal complexes comprising lherzolites and granulite facies migmatites of the ‘Kobipoto Complex’. Residual garnet–sillimanite granulites, which contain spinel + quartz inclusions within garnet, experienced ultrahigh-temperature (UHT; > 900 °C) conditions at 16 Ma due to heat supplied by lherzolites exhumed during slab rollback in the Banda Arc. Here, we present U–Pb zircon ages and new whole-rock geochemical analyses that document a protracted history of high-T metamorphism, melting, and acid magmatism of a common sedimentary protolith. Detrital zircons from the Kobipoto Complex migmatites, with ages between 3.4 Ga and 216 Ma, show that their protolith was derived from both West Papua and the Archean of Western Australia, and that metamorphism of these rocks on Seram could not have occurred until the Late Triassic. Zircons within the granulites then experienced three subsequent episodes of growth – at 215–173 Ma, 25–20 Ma, and at c. 16 Ma. The population of zircon rims with ages between 215 and 173 Ma document significant metamorphic (± partial melting) events that we attribute to subduction beneath the Bird's Head peninsula and Sula Spur, which occurred until the Banda and Argo continental blocks were rifted from the NW Australian margin of Gondwana in the Late Jurassic (from c. 160 Ma). Late Oligocene-Early Miocene collision between Australia (the Sula Spur) and SE Asia (northern Sulawesi) was then recorded by crystallisation of several 25–20 Ma zircon rims. Thereafter, a large population of c. 16 Ma zircon rims grew during subsequent and extensive Middle Miocene metamorphism and melting of the Kobipoto complex rocks beneath Seram under high- to ultrahigh-temperature (HT–UHT) conditions. Lherzolites located adjacent to the granulite-facies migmatites in central Seram equilibrated at 1280–1300 °C upon their exhumation to 1 GPa (~ 37 km) depth, whereupon they supplied sufficient heat to have metamorphosed adjacent Kobipoto Complex migmatites under UHT conditions at 16 Ma. Calculations suggesting slight (~ 10 vol%) mantle melting are consistent with observations of minor gabbroic intrusions and scarce harzburgites. Subsequent extension during continued slab rollback exhumed both the lherzolites and adjacent granulite-facies migmatites beneath extensional detachment faults in western Seram at 6.0–5.5 Ma, and on Ambon at 3.5 Ma, as recorded by subsequent zircon growth and ⁴⁰Ar/³⁹Ar ages in these regions. Ambonites, cordierite- and garnet-bearing dacites sourced predominantly from melts generated in the Kobipoto Complex migmatites, were later erupted on Ambon from 3.0 to 1.9 Ma.     

P–T history of garnet-websterites in the Sharyzhalgai complex,southwestern margin of Siberian craton: evidence for Paleoproterozoic high-pressure metamorphism

The Saramta massif in the Paleoproterozoic Sharyzhalgai complex, the southwestern margin of the Siberian craton, is mainly composed of spinel-peridotites with garnet-websterites; it is enclosed within granitic gneisses and migmatites with mafic intercalations of granulite-facies grade. The garnet-websterites occur as lenses or layers intercalated within spinel-harzburgite and spinel-lherzolite. They consist mainly of clinopyroxene (Cpx), garnet (Grt), and orthopyroxene (Opx): Grt often includes Cpx, Opx, and pargasite (Prg). Opx also occurs as kelyphite with plagioclase (Pl), spinel, olivine, Prg, and biotite. Relationships between textures and chemical compositions of these minerals suggest the following P–T stages: stage 1 (pre-peak), 0.9–1.5 GPa at 640–780 °C; stage 2 (peak), 2.3–3.0 GPa at 920–1030 °C as the minimum estimate; and stage 3 (post-peak), 750–830 °C at 0.5–0.9 GPa. Finally, the garnet-websterites are veined with lower amphibolite- to greenschist-facies minerals (stage 4).These results suggests that the Saramta massif was carried to depths of c. 100 km by subduction, and metamorphosed under eclogite-facies conditions in the Paleoproterozoic, despite the commonly held view that high geothermal gradients in those times would have prevented such deep subduction. Paleoproterozoic plate subduction at the southwestern margin of the Siberian craton might have caused subduction-zone magmatism and mantle metasomatism similar to those in the Phanerozoic.     

Evidence for Ordovician subduction-related magmatism in the Truong Son terrane,SE Laos: Implications for Gondwana evolution and porphyry Cu exploration potential in SE Asia

The Truong Son Fold Belt (TSFB) is characterised by Late Carboniferous-Late Triassic metamorphic, volcanic and plutonic rocks, the product of accretion of the Indochina Terrane onto the South China Terrane and a range of composite subduction, collision and extensional events. This study discusses geochronological and geochemical data obtained from a dioritic intrusion and rhyolitic tuff mapped in the Donken area of SE Laos, and previously assigned to the Permian Antoum Granodiorite rock suite within the TSFB. Magmatic zircon U-Pb Q-ICP-MS dating undertaken in this study suggests ages of ca 470 ± 2 Ma for the diorite and ca 476 ± 1.5 Ma for a proximal rhyolitic tuff.Whole-rock geochemistry of both units suggests a subduction-related island arc environment, with calc-alkaline and tholeiitic affinities for the diorite and tuff respectively. The intrusion also exhibits an adakitic signature (high Sr, low Y and HREE contents) suggesting that Ordovician magmatism also occurred within the Indochina Terrane, associated with an enigmatic, early Gondwana subduction event. This intrusion appears part of a broader, bilateral Early Ordovician magmatism, newly linked to the south-east subduction of the Tamky-Phuoc Son Ocean underneath the Kontum terrane, and a north-west subduction beneath the Truong Son terrane. Significantly, sub-economic hydrothermal Cu mineralisation observed within the dioritic intrusion, hints at the presence of local Ordovician, porphyry-style base metal enrichment.     

Ultra-depleted peridotite xenoliths in the Northern Taihang Mountains: Implications for the nature of the lithospheric mantle beneath the North China Craton

We report the finding of peridotite xenoliths in the Early Cretaceous Longmengou olivine-bearing diabase (138 Ma) in the Northern Taihang Mountains in the central North China Craton. Based on the modal proportions of olivine, clinopyroxene, amphibole and anorthite, these peridotite xenoliths can be divided into three zones: clinopyroxene-bearing olivine zone (COZ), olivine-clinopyroxene zone (OCZ), and amphibole-bearing anorthite-clinopyroxene zone (AACZ). The core of olivine grains in clinopyroxene-bearing olivine zone have higher Mg^# (> 95), SiO₂ (41.80–42.53 wt%) and lower CaO (< 0.07 wt%), FeO (3.91–4.54 wt%) than the rim (Mg^# = 92.5–93.4, SiO₂ = 41.27–41.98 wt%, CaO = 0.20–0.34 wt%, and FeO = 7.02–8.87 wt%), suggesting that rim is reaction product. The core of olivine grains with higher Mg^# (> 95) and lower NiO content (< 0.04 wt%) in the clinopyroxene-bearing olivine zone was derived from ultra-depleted mantle subsequently altered by high Mg^# melts/magma with low Ni. Two generations of olivine grains occur in the OCZ where the first generation shows exsolution of ilmenite and magnetite rods containing up to 0.35 wt% TiO₂, and was likely derived from garnet peridotite hydrated by water. The second generation shows high Mg^# (96.2–97.1) and cataclastic texture, and was possibly formed by decomposition of the COZ. The occurrence of aluminous spinel suggests the role of melts with extremely high Al and Mg. Clinopyroxene in the AACZ shows systematic core-rim compositional variation with CaO and SiO₂ contents increasing towards the rim, and MgO and Fe₂O₃ concentrations decreasing from the core to the rim, indicating that the amphibole-bearing anorthite-clinopyroxene zone is a product of the reaction between mantle xenoliths and mafic magma. Plagioclase with high An value (92.0–99.95, average 97.79) indicates that the metasomatic melts have high Ca/Na and Al/Si ratios, possibly produced by the partial melting of ultra-depleted mantle under “wet” conditions. Combined with the data on other mantle xenoliths discovered in the NCC, our results suggest that the Mesozoic lithospheric mantle beneath the North Taihang Mountains within the central NCC is composed of ultra-depleted Archean and Paleoproterozoic peridotites and dunites modified by complex melts. We also propose that the destruction of eastern part of the NCC mainly occurred during Early Cretaceous, and that the boundary of the lithospheric destruction coincides with the Taihang Mountains.