Mechanisms of continental subduction and exhumation of HP and UHP rocks

We discuss possible scenarios of continental collision, and their relation to mechanisms of exhumation of HP and UHP rocks, inferred from thermo-mechanical numerical models accounting for thermo-rheological complexity of the continental lithosphere. Due to this complexity, mechanisms of continental convergence are versatile and different, in many aspects from those that control oceanic subduction. Elucidating these mechanisms from conventional observations is difficult, and requires additional constraints such as those derived from petrological data. Indeed, exhumation of HP/UHP rocks is an integral part of convergent processes, and burial/exhumation dynamics inferred from metamorphic P–T–t paths provides strong constraints on the collision scenarios. Metamorphic rocks also play an active role due to their contrasting physical properties (rheology, density, fluid transport capacity). Numerical thermo-mechanical experiments suggest that HP/UHP exhumation can only be produced in subduction contexts, as well as that long-lasting (> 10 Myr) continental subduction can only occur in case of cold strong lithospheres (T_Moho < 550 °C, the equivalent elastic thickness T_e > 50 km) and of relatively high convergence rates (> 3–5 cm yr^− 1 ). In this case, high density UHP material in the crustal part of subduction interface provides additional pull on the slab and is not always exhumed to the surface. In case of slower convergence and/or weaker lithosphere (T_e < 40 km), continental subduction is a transient process that takes a limited time span in the evolution of collision zone. Under these conditions, hot mechanically weak UHP rocks enhance decoupling between the upper and lower plate while their exhumation may be rapid (faster than convergence rate) and abundant. Therefore, the UHP exhumation paths can be regarded as sensitive indicators of subduction. Rheological changes and fluid exchanges associated with low-to-middle pressure phase transitions along the subduction interface, such as serpentinization during the oceanic phase and schisting, play a major role producing necessarily mechanical softening of the subduction interface and of the hydrated mantle wedge. The oceanic UHP rocks are exhumed thanks to mixing with low-density continental crustal units during transition from oceanic to continental subduction. At the continental phase, the UHP exhumation occurs as a result of a multi-stage process: at the deep stage (< 40 km depth) the exhumation is rapid and is driven by buoyancy of partly metamorphosed (or partly molten) UHP material often mixed with non-metamorphosed crustal volumes. At final stages, exhumation takes common slow path through the accretion prism mechanism and the erosional denudation. The experiments suggest that formation of UHP rocks requires that continental subduction starts at higher oceanic subduction rate. It then may progressively slow down until the lockup of the subduction interface and/or slab-break-off. A rate of ~ 1–2 cm yr^− 1 is generally sufficient to drive continental subduction during the first several Myr of convergence, but pertinent subduction requires faster convergence rates (> 3–5 cm yr^− 1). We suggest that most continental orogenic belts could have started their formation from continental subduction but this process has been generally limited in time.     

Contrasting modes of eclogite and blueschist exhumation in a retreating subduction system: The Tasmanides,Australia

Blueschists and eclogites located in the Tasmanides of eastern Australia preserve evidence of contrasting modes of exhumation. A review of structural, metamorphic, geochronological and geochemical data indicates that these HP metamorphic rocks can be sub-divided into three main groups: (i) eclogite–blueschists with calc-alkaline and tholeiitic affinities contained within thick sedimentary sequences (called continental HP rocks); (ii) moderate-pressure (< 9 kbar) blueschist of arc to MORB-type composition within sedimentary or serpentinite mélange zones (called accretionary HP rocks) and (iii) eclogites of MORB-type composition with or without a pervasive blueschist overprint contained within serpentinite (called exotic HP rocks). Three different modes of exhumation can be ascribed to the different rock types, namely: (i) exhumation influenced by the buoyancy of continental slabs; (ii) exhumation of accretionary HP rocks by corner flow and/or extensional collapse in the accretionary wedge or (iii) discontinuous exhumation of eclogites triggered by slab rollback and trench retreat. We suggest that a dominant west-dipping, eastward migrating subduction zone can explain the distribution and formation of HP metamorphic rocks in the Tasmanides.Thermobarometric and geochronological data from eclogites and blueschists in the Peel–Manning Fault System (New England Orogen) also provide evidence for discontinuous exhumation of subducted oceanic rocks. These data indicate that eclogites were exhumed from depths of ~ 70 km to ~ 30 km during the Ordovician (490–470 Ma), with terminal exhumation and exposure along the Peel–Manning Fault system probably occurring during the Permian. Based on these timing constraints, we suggest a model where HP rocks reside between depth-dependant exhumation circuits for considerable lengths of time.     

High-temperature,high-pressure granulites (retrogressed eclogites) in the central region of the Lewisian,NW Scotland: Crustal-scale subduction in the Neoarchaean

Eclogites and associated high-pressure (HP) rocks in collisional and accretionary orogenic belts preserve a record of subduction and exhumation, and provide a key constraint on the tectonic evolution of the continents. Most eclogites that formed at high pressures but low temperatures at > 10–11 kbar and 450–650 °C can be interpreted as a result of subduction of cold oceanic lithosphere. A new class of high-temperature (HT) eclogites that formed above 900 °C and at 14 to 30 kbar occurs in the deep continental crust, but their geodynamic significance and processes of formation are poorly understood. Here we show that Neoarchaean mafic–ultramafic complexes in the central granulite facies region of the Lewisian in NW Scotland contain HP/HT garnet-bearing granulites (retrogressed eclogites), gabbros, lherzolites, and websterites, and that the HP granulites have garnets that contain inclusions of omphacite. From thermodynamic modeling and compositional isopleths we calculate that peak eclogite-facies metamorphism took place at 24–22 kbar and 1060–1040 °C. The geochemical signature of one (G-21) of the samples shows a strong depletion of Eu indicating magma fractionation at a crustal level. The Sm–Nd isochron ages of HP phases record different cooling ages of ca. 2480 and 2330 Ma. We suggest that the layered mafic–ultramafic complexes, which may have formed in an oceanic environment, were subducted to eclogite depths, and exhumed as HP garnet-bearing orogenic peridotites. The layered complexes were engulfed by widespread orthogneisses of tonalite–trondhjemite–granodiorite (TTG) composition with granulite facies assemblages. We propose two possible tectonic models: (1) the fact that the relicts of eclogitic complexes are so widespread in the Scourian can be taken as evidence that a > 90 km × 40 km-size slab of continental crust containing mafic–ultramafic complexes was subducted to at least 70 km depth in the late Archaean. During exhumation the gneiss protoliths were retrogressed to granulite facies assemblages, but the mafic–ultramafic rocks resisted retrogression. (2) The layered complexes of mafic and ultramafic rocks were subducted to eclogite-facies depths and during exhumation under crustal conditions they were intruded by the orthogneiss protoliths (TTG) that were metamorphosed in the granulite facies. Apart from poorly defined UHP metamorphic rocks in Norway, the retrogressed eclogites in the central granulite/retrogressed eclogite facies Lewisian region, NW Scotland have the highest crustal pressures so far reported for Archaean rocks, and demonstrate that lithospheric subduction was transporting crustal rocks to HP depths in the Neoarchaean.     

From oceanic plateaus to allochthonous terranes: Numerical modelling

Large segments of the continental crust are known to have formed through the amalgamation of oceanic plateaus and continental fragments. However, mechanisms responsible for terrane accretion remain poorly understood. We have therefore analysed the interactions of oceanic plateaus with the leading edge of the continental margin using a thermomechanical–petrological model of an oceanic-continental subduction zone with spontaneously moving plates. This model includes partial melting of crustal and mantle lithologies and accounts for complex rheological behaviour including viscous creep and plastic yielding. Our results indicate that oceanic plateaus may either be lost by subduction or accreted onto continental margins. Complete subduction of oceanic plateaus is common in models with old (> 40 Ma) oceanic lithosphere whereas models with younger lithosphere often result in terrane accretion. Three distinct modes of terrane accretion were identified depending on the rheological structure of the lower crust and oceanic cooling age: frontal plateau accretion, basal plateau accretion and underplating plateaus.Complete plateau subduction is associated with a sharp uplift of the forearc region and the formation of a basin further landward, followed by topographic relaxation. All crustal material is lost by subduction and crustal growth is solely attributed to partial melting of the mantle.Frontal plateau accretion leads to crustal thickening and the formation of thrust and fold belts, since oceanic plateaus are docked onto the continental margin. Strong deformation leads to slab break off, which eventually terminates subduction, shortly after the collisional stage has been reached. Crustal parts that have been sheared off during detachment melt at depth and modify the composition of the overlying continental crust.Basal plateau accretion scrapes oceanic plateaus off the downgoing slab, enabling the outward migration of the subduction zone. New incoming oceanic crust underthrusts the fractured terrane and forms a new subduction zone behind the accreted terrane. Subsequently, hot asthenosphere rises into the newly formed subduction zone and allows for extensive partial melting of crustal rocks, located at the slab interface, and only minor parts of the former oceanic plateau remain unmodified.Oceanic plateaus may also underplate the continental crust after being subducted to mantle depth. (U)HP terranes are formed with peak metamorphic temperatures of 400–700 °C prior to slab break off and subsequent exhumation. Rapid and coherent exhumation through the mantle along the former subduction zone at rates comparable to plate tectonic velocities is followed by somewhat slower rates at crustal levels, accompanied by crustal flow, structural reworking and syndeformational partial melting. Exhumation of these large crustal volumes leads to a sharp surface uplift.     

Initial water budget: The key to detaching large volumes of eclogitized oceanic crust along the subduction channel?

We herein investigate the extent to which extensive hydration of the oceanic lithosphere influences the preservation and exhumation of large-scale ophiolite bodies from subduction zones. The Zermatt–Saas ophiolite (ZS, W. Alps), which was subducted during the late stages of oceanic subduction, preserves a complete section of Mesozoic Tethys oceanic lithosphere and particularly fresh eclogites, and represents, so far, the largest and deepest known portion of exhumed oceanic lithosphere. Pervasive hydrothermal processes and seafloor alteration led to the incorporation of large amounts of fluid bound in the hydrated upper layers of the oceanic crust (now as lawsonite eclogites, glaucophanites, and chloritoschists) and in associated ultramafic rocks.Internally, the ZS ophiolite is made up of a series of tectonic slices of oceanic crust (150–300 m thick) which are systematically separated by a 5 to 100 m thick layer of serpentinite. This stack of slices is separated from the underlying eclogitized continental crust (e.g., Monte Rosa) by a thick (~ 500 m) serpentinite sole. Field observations, textural relationships and pseudosection modelling reveal that lawsonite was abundant and widespread in mafic eclogites when the ophiolite detached from the slab at around 550 °C and 24 kbar.Comparison between fresh eclogitic samples and pseudosection modelling shows that (i) water remained in excess from burial to eclogitic peak conditions, (ii) the lightest eclogitized metabasalts correspond to the portions of oceanic crust where metasomatism was the strongest, (iii) crystallization of widespread hydrated parageneses (such as lawsonite, glaucophane and phengite) instead of garnet and omphacite decreased by 5 to 10% the rock density and subsequently enhanced its buoyancy.We propose that this density decrease acted as a ‘float’ which prevented the slices from an irreversible sinking in the mantle. These slices were subsequently detached from the downgoing slab and stacked in the serpentinized subduction channel at pressures between 15 and 20 kbar, in the epidote blueschist facies. Exhumation of the underlying, positively buoyant continental crust dragged this “frozen” nappe-stack from the subduction channel towards the surface.     

Tectonics of the North Qilian orogen,NW China

The Qilian Orogen at the northern margin of the Tibetan Plateau is a type suture zone that recorded a complete history from continental breakup to ocean basin evolution, and to the ultimate continental collision in the time period from the Neoproterozoic to the Paleozoic. The Qilian Ocean, often interpreted as representing the “Proto-Tethyan Ocean”, may actually be an eastern branch of the worldwide “Iapetus Ocean” between the two continents of Baltica and Laurentia, opened at ≥ 710 Ma as a consequence of breakup of supercontinent Rodinia.Initiation of the subduction in the Qilian Ocean probably occurred at ~ 520 Ma with the development of an Andean-type active continental margin represented by infant arc magmatism of ~ 517–490 Ma. In the beginning of Ordovician (~ 490 Ma), part of the active margin was split from the continental Alashan block and the Andean-type active margin had thus evolved to western Pacific-type trench–arc–back-arc system represented by the MORB-like crust (i.e., SSZ-type ophiolite belt) formed in a back-arc basin setting in the time period of ~ 490–445 Ma. During this time, the subducting oceanic lithosphere underwent LT-HP metamorphism along a cold geotherm of ~ 6–7 °C/km.The Qilian Ocean was closed at the end of the Ordovician (~ 445 Ma). Continental blocks started to collide and the northern edge of the Qilian–Qaidam block was underthrust/dragged beneath the Alashan block by the downgoing oceanic lithosphere to depths of ~ 100–200 km at about 435–420 Ma. Intensive orogenic activities occurred in the late Silurian and early Devonian in response to the exhumation of the subducted crustal materials.Briefly, the Qilian Orogen is conceptually a type example of the workings of plate tectonics from continental breakup to the development and evolution of an ocean basin, to the initiation of oceanic subduction and formation of arc and back-arc system, and to the final continental collision/subduction and exhumation.     

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.     

Boundaries of mantle–lithosphere domains in the Bohemian Massif as extinct exhumation channels for high-pressure rocks

Five domains (microplates) have been recognized by seismic anisotropy in the mantle lithosphere of the Bohemian Massif. The mantle domains correspond to major crustal units and each of the domains bears a consistent fossil olivine fabric formed before their Variscan assembly. The present-day mantle fabric indicates that this process consisted of at least three oceanic subductions, each followed by an underthrusting of the continental lithosphere. The seismic anisotropy does not detect remnants of the oceanic subductions, but it can trace boundaries of the preserved continental domains subsequently underthrust along the paths of previous oceanic subductions. The most robust continent–continent collision was followed by westward underthrusting of the Brunovistulian mantle lithosphere, still detectable by seismic anisotropy more than 100 km beneath the Moldanubian mantle lithosphere. Major occurrences of the high-pressure/ultra high-pressure (HP–UHP) rocks follow the ENE and NNE oriented sutures and boundaries of the mantle–lithosphere domains mapped from three-dimensional modeling of body-wave anisotropy. The HP–UHP rocks are products of oceanic subductions and the following underthrusting of the continental crust and mantle lithosphere exhumed along the mantle boundaries. The close relation of the mantle sutures and occurrences of the HP–UHP rocks near the paleosubductions testifies for models interpreting the granulite–garnet peridotite association by oceanic/continental subduction/underthrusting followed by the exhumation of deep-seated rocks. Our findings support the bivergent subduction model of tectonic development of the central part of the Bohemian Massif. The inferences from seismic anisotropy image the Bohemian Massif as a mosaic of microplates with a rigid mantle lithosphere preserving a fossil olivine fabric. The collisional mantle boundaries, blurred by tectonometamorphic processes in easily deformed overlying crust, served as major exhumation channels of the HP–UHP rocks.     

Effect of time-evolving age and convergence rate of the subducting plate on the Cenozoic adakites and boninites

Partial melting of subducting oceanic crust expressed as high-Mg volcanic rocks such as adakites and boninites has been actively studied for decades, and Lee and King (2010) reported that time-evolving subduction parameters such as the age and the subduction rate of the converging oceanic plate play important roles in transient partial melting of the subducting oceanic crust (e.g., Aleutians). However, few subduction model experiments have considered time-evolving subduction parameters, posing problems for studies of transient partial melting of subducting oceanic crust in many subduction zones. Therefore, we constructed two-dimensional kinematic–dynamic subduction models for the Izu–Bonin, Mariana, Northeast Japan, Kuril, Tonga, Java–Sunda, and Aleutian subduction zones that account for the last 50 Myr of their evolution. The models include the time-evolving age and convergence rate of the incoming oceanic plate, so the effect of time-evolving subduction parameters on transient partial melting of oceanic crust can be evaluated. Our model calculations revealed that adakites and boninites in the Izu–Bonin and Aleutian subduction zones resulted from transient partial melting of oceanic crust. However, the steady-state subduction model using current subduction parameters did not produce any partial melting of oceanic crust in the aforementioned subduction zones, indicating that time-evolving subduction parameters are crucial for modeling transient eruption of adakites and boninites. Our model calculations confirm that other geological processes such as forearc extension, back-arc opening, mantle plumes and ridge subduction are required for partial melting of the oceanic crust in the Mariana, Northeast Japan, Tonga, and southeastern Java–Sunda subduction zones.     

Early Mesozoic southward subduction history of the Mongol–Okhotsk oceanic plate: Evidence from geochronology and geochemistry of Early Mesozoic intrusive rocks in the Erguna Massif,NE China

In this paper we present new zircon U–Pb ages, Hf isotope data, and whole-rock major and trace element data for Early Mesozoic intrusive rocks in the Erguna Massif of NE China, and we use these data to constrain the history of southward subduction of the Mongol–Okhotsk oceanic plate, and its influence on NE China as a whole. The zircon U–Pb dating indicates that Early Mesozoic magmatic activity in the Erguna Massif can be subdivided into four stages at ~ 246 Ma, ~ 225 Ma, ~ 205 Ma, and ~ 185 Ma. The ~ 246 Ma intrusive rocks comprise a suite of high-K calc-alkaline diorites, quartz diorites, granodiorites, monzogranites, and syenogranites, with I-type affinities. The ~ 225 Ma intrusive rocks consist of gabbro–diorites and granitoids, and they constitute a bimodal igneous association. The ~ 205 Ma intrusive rocks are dominated by calc-alkaline I-type granitoids that are accompanied by subordinate intermediate–mafic rocks. The ~ 185 Ma intrusive rocks are dominated by I-type granitoids, accompanied by minor amounts of A-types. These Early Mesozoic granitoids mainly originated by partial melting of a depleted and heterogeneous lower crust, whereas the coeval mafic rocks were probably derived from partial melting of a depleted mantle modified by subduction-related fluids. The rock associations and their geochemical features indicate that the ~ 246 Ma, ~ 205 Ma, and ~ 185 Ma intrusive rocks formed in an active continental margin setting related to the southward subduction of the Mongol–Okhotsk oceanic plate. The ~ 225 Ma bimodal igneous rock association formed within an extensional environment in a pause during the subduction process of the Mongol–Okhotsk oceanic plate. Every magmatic stage has its own corresponding set of porphyry deposits in the southeast of the Mongol–Okhotsk suture belt. Taking all this into account, we conclude the following: (1) during the Early Mesozoic, the Mongol–Okhotsk oceanic plate was subducted towards the south beneath the Erguna Massif, but with a pause in subduction at ~ 225 Ma; and (2) the southward subduction of the Mongol–Okhotsk oceanic plate not only caused the intense magmatic activity, but was also favorable to the formation of porphyry deposits.     

The Western Sierras Pampeanas: Protracted Grenville-age history (1330–1030 Ma) of intra-oceanic arcs,subduction–accretion at continental-edge and AMCG intraplate magmatism

New U–Pb SHRIMP zircon ages combined with geochemical and isotope investigation in the Sierra de Maz and Sierra de Pie de Palo and a xenolith of the Precordillera basement (Ullún), provides insight into the identification of major Grenville-age tectonomagmatic events and their timing in the Western Sierras Pampeanas. The study reveals two contrasting scenarios that evolved separately during the 300 Ma long history: Sierra de Maz, which was always part of a continental crust, and the juvenile oceanic arc and back-arc sector of Sierra de Pie de Palo and Ullún. The oldest rocks are the Andino-type granitic orthogneisses of Sierra de Maz (1330–1260 Ma) and associated subalkaline basic rocks, that were part of an active continental margin developed in a Paleoproterozoic crust. Amphibolite facies metamorphism affected the orthogneisses at ca. 1175 Ma, while granulite facies was attained in neighbouring meta-sediments and basic granulites. Interruption of continental-edge magmatism and high-grade metamorphism is interpreted as related to an arc–continental collision dated by zircon overgrowths at 1170–1230 Ma. The next event consisted of massif-type anorthosites and related meta-jotunites, meta-mangerites (1092 ± 6 Ma) and meta-granites (1086 ± 10 Ma) that define an AMCG complex in Sierra de Maz. The emplacement of these mantle-derived magmas during an extensional episode produced a widespread thermal overprint at ca. 1095 Ma in neighbouring country rocks. In constrast, juvenile oceanic arc and back-arc complexes dominated the Sierra de Pie de Palo–Ullún sector, that was fully developed ca. 1200 Ma (1196 ± 8 Ma metagabbro). A new episode of oceanic arc magmatism at ～1165 Ma was roughly coeval with the amphibolite high-grade metamorphism of Sierra de Maz, indicating that these two sectors underwent independent geodynamic scenarios at this age. Two more episodes of arc subduction are registered in the Pie de Palo–Ullún sector: (i) 1110 ± 10 Ma orthogneisses and basic amphibolites with geochemical fingerprints of emplacement in a more mature crust, and (ii) a 1027 ± 17 Ma TTG juvenile suite, which is the youngest Grenville-age magmatic event registered in the Western Sierras Pampeanas. The geodynamic history in both study areas reveals a complex orogenic evolution, dominated by convergent tectonics and accretion of juvenile oceanic arcs to the continent.     

Metamorphism and tectonic evolution of the Lhasa terrane,Central Tibet

The Lhasa terrane in southern Tibet is composed of Precambrian crystalline basement, Paleozoic to Mesozoic sedimentary strata and Paleozoic to Cenozoic magmatic rocks. This terrane has long been accepted as the last crustal block to be accreted with Eurasia prior to its collision with the northward drifting Indian continent in the Cenozoic. Thus, the Lhasa terrane is the key for revealing the origin and evolutionary history of the Himalayan–Tibetan orogen. Although previous models on the tectonic development of the orogen have much evidence from the Lhasa terrane, the metamorphic history of this terrane was rarely considered. This paper provides an overview of the temporal and spatial characteristics of metamorphism in the Lhasa terrane based mostly on the recent results from our group, and evaluates the geodynamic settings and tectonic significance. The Lhasa terrane experienced multistage metamorphism, including the Neoproterozoic and Late Paleozoic HP metamorphism in the oceanic subduction realm, the Early Paleozoic and Early Mesozoic MP metamorphism in the continent–continent collisional zone, the Late Cretaceous HT/MP metamorphism in the mid-oceanic ridge subduction zone, and two stages of Cenozoic MP metamorphism in the thickened crust above the continental subduction zone. These metamorphic and associated magmatic events reveal that the Lhasa terrane experienced a complex tectonic evolution from the Neoproterozoic to Cenozoic. The main conclusions arising from our synthesis are as follows: (1) The Lhasa block consists of the North and South Lhasa terranes, separated by the Paleo-Tethys Ocean and the subsequent Late Paleozoic suture zone. (2) The crystalline basement of the North Lhasa terrane includes Neoproterozoic oceanic crustal rocks, representing probably the remnants of the Mozambique Ocean derived from the break-up of the Rodinia supercontinent. (3) The oceanic crustal basement of North Lhasa witnessed a Late Cryogenian (~ 650 Ma) HP metamorphism and an Early Paleozoic (~ 485 Ma) MP metamorphism in the subduction realm associated with the closure of the Mozambique Ocean and the final amalgamation of Eastern and Western Gondwana, suggesting that the North Lhasa terrane might have been partly derived from the northern segment of the East African Orogen. (4) The northern margin of Indian continent, including the North and South Lhasa, and Qiangtang terranes, experienced Early Paleozoic magmatism, indicating an Andean-type orogeny that resulted from the subduction of the Proto-Tethys Ocean after the final amalgamation of Gondwana. (5) The Lhasa and Qiangtang terranes witnessed Middle Paleozoic (~ 360 Ma) magmatism, suggesting an Andean-type orogeny derived from the subduction of the Paleo-Tethys Ocean. (6) The closure of Paleo-Tethys Ocean between the North and South Lhasa terranes and subsequent terrane collision resulted in the formation of Late Permian (~ 260 Ma) HP metamorphic belt and Triassic (220 Ma) MP metamorphic belt. (7) The South Lhasa terrane experienced Late Cretaceous (~ 90 Ma) Andean-type orogeny, characterized by the regional HT/MP metamorphism and coeval intrusion of the voluminous Gangdese batholith during the northward subduction of the Neo-Tethyan Ocean. (8) During the Early Cenozoic (55–45 Ma), the continent–continent collisional orogeny has led to the thickened crust of the South Lhasa terrane experiencing MP amphibolite-facies metamorphism and syn-collisional magmatism. (9) Following the continuous continent convergence, the South Lhasa terrane also experienced MP metamorphism during Late Eocene (40–30 Ma). (10) During Mesozoic and Cenozoic, two different stages of paired metamorphic belts were formed in the oceanic or continental subduction zones and the middle and lower crust of the hanging wall of the subduction zone. The tectonic imprints from the Lhasa terrane provide excellent examples for understanding metamorphic processes and geodynamics at convergent plate boundaries.     

Late Triassic alkaline complex in the Sulu UHP terrane: Implications for post-collisional magmatism and subsequent fractional crystallization

Continental subduction and its interaction with overlying mantle wedge are recognized as fundamental solid earth processes, yet the dynamics of this system remains ambiguous. In order to get an insight into crust–mantle interaction triggered by partial melting of subudcted continental crust during its exhumation, we carried out a combined study of the Shidao alkaline complex from the Sulu ultrahigh pressure (UHP) terrane. The alkaline complex is composed of shoshonitic to ultrapotassic (K₂O: 3.4–9.3 wt.%) gabbro, pyroxene syenite, amphibole syenite, quartz syenite, and granite. Field studies suggest that the mafic rocks are earlier than the felsic ones in sequence. LA-ICPMS zircon U–Pb dating on them gives Late Triassic ages of 214 ± 2 to 200 ± 3 Ma from mafic to felsic rocks. These ages are slightly younger than the Late Triassic ages (225–210 Ma) of the felsic melts from partial melting of the Sulu UHP terrane during exhumation. The alkaline rocks have wide ranges of SiO₂ (49.7–76.7 wt.%), MgO (8.25–0.03 wt.%), Ni (126.0–0.07 ppm), and Cr (182.0–0.45 ppm) contents. The contents of MgO, total Fe₂O₃, CaO, TiO₂ and P₂O₅ decrease with increasing SiO₂ contents. The contents of Na₂O, K₂O, and Al₂O₃ increase from gabbro to amphibole syenite, and decrease from amphibole syenite to granite, respectively. The alkaline rocks have characteristics of an arc-like pattern in trace element distribution, e.g., enrichment of LREE, LILE (Rb and Ba), Th and U, depletion of HFSE (Nb, Ta, P and Ti), and positive Pb anomalies. From the mafic rocks to the felsic rocks, the (La/Yb)_N ratios and the contents of the total REE, Sr and Ba decrease but the Rb contents increase. The alkaline rocks with high SiO₂ contents also display features of an A2-type granitoids, e.g., high contents of total alkalis, Zr and Nb and high ratios of Fe₂O₃^T/MgO, Ga/Al, Yb/Ta and Y/Nb, suggesting a post-collisional magmatism during exhumation of the Sulu UHP terrane. The alkaline rocks have homogeneous initial ⁸⁷Sr/⁸⁶Sr ratios (0.7058–0.7093) and negative ε_Nd(t) values (− 18.6 to − 15.0) for whole-rock. The Sr–Nd isotopic data remain almost unchanged with varying SiO₂ and MgO contents, suggesting a fractional crystallization (FC) process from the same parental magma. Our studies suggest a crust–mantle interaction in continental subduction interface as follows: (1) hydrous felsic melts from partial melting of subducted continental crust during its exhumation metasomatized the overlying mantle wedge to form a K-rich and amphibole-bearing mantle; (2) partial melting of the enriched lithospheric mantle generated the Late Triassic alkaline complex under a post-collisional setting; and (3) the alkaline magma experienced subsequent fractionational crystallization mainly dominated by olivine, clinopyroxene, plagioclase and alkali feldspar.     

Tectonic evolution of the North Qinling Orogen from subduction to collision and exhumation: Evidence from zircons in metamorphic rocks of the Qinling Group

The North Qinling Block (NQB) is an important segment of the Qinling Orogen in Central China. Here we report the results from SIMS geochronology and oxygen isotopes, as well as LA-MC-ICPMS Hf isotopic analyses on zircon grains from a suite of metamorphic rocks (felsic gneisses, garnet plagioclase amphibolites, and retrograde eclogite dikes) in the Qinling Group of the NQB. The age data show that these rocks underwent at least two episodes of metamorphism with the peak at 483–501 Ma, followed by 454–470 Ma retrograde metamorphism. These results are generally coeval with the periods of 500–480 Ma for peak metamorphism and 460–420 Ma for retrograde metamorphism previously obtained from the HP/UHP metamorphic rocks of the NQB. During the prograde and retrograde metamorphism, widespread fluid and melt circulation within the block has been identified from the geochemical features of the metamorphic zircons. The fluids that circulated in the felsic gneisses and retrograde eclogite dikes originated from the dehydration of altered oceanic basalts as inferred from the exceedingly low Th/U ratios, positive ε_Hf(t) (> 5) and extremely δ¹⁸O (10.01–13.91‰) values in metamorphic zircons. In contrast, the melt involved in the formation of garnet plagioclase amphibolites appears to have been derived from continental sediments interlayered with the oceanic basalts since zircons crystallized during the peak and retrograde metamorphism show typical magmatic features with high U and Th contents and Th/U ratios and enriched Hf (ε_Hf(t) = − 5.42 to − 0.18) and oxygen isotope composition (δ¹⁸O around 8‰). Geochronological and geochemical features of the magmatic cores of the clear core-rim textured zircons demonstrate that the protoliths of the gneisses were intermediate-acid volcanic rocks erupted before Neoproterozoic (800 Ma), which is further supported by the intrusion of basaltic magma of asthenospheric origin as represented by protoliths of retrograde eclogite dikes, with the oldest magmatic zircon formed at 789 Ma. The protoliths of garnet plagioclase amphibolites appear to be altered oceanic basalts but had been significantly affected by the melt during the metamorphism. Combined with the previous studies, the Qinling Group experienced overall subduction in the Early Paleozoic. The NQB as represented by the Qinling Group was most likely a discrete micro-block in the Neoproterozoic, and underwent deep subduction in the Cambrian (483–501 Ma) and exhumation in Ordovician (454–470 Ma). We propose that the NQB preserves a complete cycle of tectonic evolution of an orogen from an oceanic basin spreading, and micro-continent formation to deep subduction and exhumation.     

Formation and evolution of Precambrian continental lithosphere in South China

An overview is presented for the formation and evolution of Precambrian continental lithosphere in South China. This is primarily based on an integrated study of zircon U–Pb ages and Lu–Hf isotopes in crustal rocks, with additional constraints from Re–Os isotopes in mantle-derived rocks. Available Re–Os isotope data on xenolith peridotites suggest that the oldest subcontinental lithospheric mantle beneath South China is primarily of Paleoproterozoic age. The zircon U–Pb ages and Lu–Hf isotope studies reveal growth and reworking of the juvenile crust at different ages. Both the Yangtze and Cathaysia terranes contain crustal materials of Archean U–Pb ages. Nevertheless, zircon U–Pb ages exhibit two peaks at 2.9–3.0 Ga and ~ 2.5 Ga in Yangtze but only one peak at ~ 2.5 Ga in Cathaysia. Both massive rocks and crustal remnants (i.e., zircon) of Archean U–Pb ages occur in Yangtze, but only crustal remnants of Archean U–Pb ages occur in Cathaysia. Zircon U–Pb and Lu–Hf isotopes in the Kongling complex of Yangtze suggest the earliest episode of crustal growth in the Paleoarchean and two episodes of crustal reworking at 3.1–3.3 Ga and 2.8–3.0 Ga. Both negative and positive ε_Hf(t) values are associated with Archean U–Pb ages of zircon in South China, indicating both the growth of juvenile crust and the reworking of ancient crust in the Archean. Paleoproterozoic rocks in Yangtze exhibit four groups of U–Pb ages at 2.1 Ga, 1.9–2.0 Ga, ~ 1.85 Ga and ~ 1.7 Ga, respectively. They are associated not only with reworking of the ancient Archean crust in the interior of Yangtze, but also with the growth of the contemporaneous juvenile crust in the periphery of Yangtze. In contrast, Paleoproterozoic rocks in Cathaysia were primarily derived from reworking of Archean crust at 1.8–1.9 Ga. The exposure of Mesoproterozoic rocks are very limited in South China, but zircon Hf model ages suggest the growth of juvenile crust in this period due to island arc magmatism of the Grenvillian oceanic subduction. Magmatic rocks of middle Neoproterozoic U–Pb ages are widespread in South China, exhibiting two peaks at about 830–800 Ma and 780–740 Ma, respectively. Both negative and positive ε_Hf(t) values are associated with the middle Neoproterozoic U–Pb ages of zircon, suggesting not only growth and reworking of the juvenile Mesoproterozoic crust but also reworking of the ancient Archean and Paleoproterozoic crust in the middle Neoproterozoic. The tectonic setting for this period of magmatism would be transformed from arc–continent collision to continental rifting with reference to the plate tectonic regime in South China.     

Amphibole and phlogopite in “hybrid” metasomatic bands monitor trace element transfer at the interface between felsic and ultramafic rocks (Eastern Alps,Italy)

Ultramafic rocks in contact to gneisses in the Mt. Hochwart HP mélange (Eastern Italian Alps) preserve a series of metasomatic mineral zones. A phlogopitite with minor tremolite and accessory zircon and apatite forms close to the gneiss (Zone 1). Zone 2 consists of tremolite, phlogopite and anthophyllite followed by Zone 3 with anthophyllite plus minor chlorite and talc. Zone 3 grades into an amphibole–garnet peridotite lens. This reaction zone has been generated by infiltration of hydrous fluids at T of 660–700 °C and P < 1.2 GPa, which occurred during exhumation of coupled continental crust and mantle peridotites.The reaction zone between a trace element-rich (gneiss) and a trace element depleted reservoir (peridotite) allows assessment of local trace element mobility in aqueous fluids. We present the results of in situ LA-ICP-MS trace element analysis of minerals from the reaction zone. Phlogopite is the main host for Large Ion Lithophile Elements (LILE) and contributes significantly to the Li, Ti, Nb, Ta, Pb and Sc budget. Anthophyllite is the main host for Li whereas all other trace elements including Rare Earth Elements (REE) are preferentially incorporated into tremolite. Combined with the abundance of these minerals over the contact zone, the mineral trace element data suggests that the LILE and REE were mobile on a small scale of a few centimetres only. Limited mobility of Ta, which is generally regarded as barely mobile in fluids, is documented in elevated contents of Ta in anthophyllite coupled with low Nb/Ta. The high Li content in minerals throughout the reaction zone suggests that Li was the most mobile element.The studied metasomatic zones mirror geochemical processes occurring in subduction zones at the slab–mantle interface. Phlogopite crystallization at the slab–mantle interface is an efficient mechanism to filter LILE from the aqueous fluid. Thus, such reaction zones, forming at temperatures < 660–700 °C, likely prevents that the typical slab signature with enriched LILE is transported by aqueous fluids over long distances in the mantle wedge. However, if coupled to the downgoing slab, phlogopite- and tremolite-rich rocks from such reaction zones might be able to act as carriers of trace elements and water into deeper parts of the subduction zone.     

Late Paleozoic–Mesozoic tectonic evolution of the Trans-Altai and South Gobi Zones in southern Mongolia based on structural and geochronological data

The Trans-Altai Zone in southern Mongolia is characterized by thrusting of greenschist-facies Silurian oceanic rocks over Devonian and Lower Carboniferous volcano-sedimentary sequences, by E–W directed folding affecting the early Carboniferous volcanic rocks, and by the development of N–S trending magmatic fabrics in the Devonian–Carboniferous arc plutons. This structural pattern is interpreted as the result of early Carboniferous thick-skinned E–W directed nappe stacking of oceanic crust associated with syn-compressional emplacement of a magmatic arc. The southernmost South Gobi Zone represents a Proterozoic continental domain affected by shallow crustal greenschist-facies detachments of Ordovician and Devonian cover sequences from the Proterozoic substratum, whereas supracrustal Carboniferous volcanic rocks and Permian sediments were folded into N–S upright folds. This structural pattern implies E–W directed thin-skinned tectonics operating from the late Carboniferous to the Permian, as demonstrated by K–Ar ages ranging from ~ 320 Ma to 257 Ma for clay fractions separated from a variety of rock types. Moreover, the geographical distribution of granitoids combined with their geochemistry and SHRIMP U–Pb zircon ages form distinct groups of Carboniferous and Permian age that record typical processes of magma generation and increase in crustal thickness. The field observations combined with clay ages, the geochemical characteristics of the granitoids and their ages imply that the E–W trending zone affected by tectonism migrated southwards, leaving the Trans Altai Zone inactive during the late Carboniferous and Permian, suggesting that the two units were tectonically amalgamated along a major E–W trending strike slip fault zone. This event was related to late Carboniferous subduction that was responsible for the vast volume of granitoid magma emplaced at 300–305 Ma in the South Gobi and at 307–308 Ma in the Trans-Altai Zones. The formation and growth of the crust was initially due only to subduction and accretion processes. During the post-collisional period from 305 to 290 Ma the addition of heat to the crust led to the generation of (per-) alkaline melts. Once amalgamated, these two different crustal domains were affected by N–S compression during the Triassic and early Jurassic (185–173 Ma), resulting in E–W refolding of early thrusts and folds and major shortening of both tectonic zones.     

Zircon geochemistry of two contrasting types of eclogite: Implications for the tectonic evolution of the North Qaidam UHPM belt,northern Tibet

Compared to the extensively documented ultrahigh-pressure metamorphism at North Qaidam, the pre-metamorphic history for both continental crust and oceanic crust is poorly constrained. Trace element compositions, U–Pb ages, O and Lu–Hf isotopes obtained for distinct zircon domains from eclogites metamorphosed from both continental and oceanic mafic rocks are linked to unravel the origin and multi-stage magmatic/metamorphic evolution of eclogites from the North Qaidam ultrahigh-pressure metamorphic (UHPM) belt, northern Tibet.For continental crust-derived eclogite, magmatic zircon cores from two samples with U–Pb ages of 875–856 Ma have both very high δ¹⁸O (10.6 ± 0.5‰) and mantle-like δ¹⁸O (averaging at 5.2 ± 0.7‰), high Th/U and ¹⁷⁶Lu/¹⁷⁷Hf ratios, and steep MREE-HREE distribution patterns (chondrite-normalized) with negative Eu anomalies. Combined with positive ε_Hf (t) of 3.9–14.3 and T_DM (1.2–0.8 Ga and 1.3–1.0 Ga, respectively), they are interpreted as being crystallized from either subduction-related mantle wedge or recycled material in the mantle. While the metamorphic rims from the eclogites have U–Pb ages of 436–431 Ma, varying (inherited, lower, and elevated) oxygen isotopes compared with cores, low Th/U and ¹⁷⁶Lu/¹⁷⁷Hf ratios, and flat HREE distribution patterns with no Eu anomalies. These reflect both solid-state recrystallization from the inherited zircon and precipitation from external fluids at metamorphic temperatures of 595–622 °C (T_Ti-in-zircon).For oceanic crust-derived eclogite, the magmatic cores (510 ± 19 Ma) and metamorphic rims (442.0 ± 3.7 Ma) also show distinction for Th/U and ¹⁷⁶Lu/¹⁷⁷Hf ratios, and the REE patterns and Eu anomalies. Combined with the mantle-like δ¹⁸O signature of 5.1 ± 0.3 ‰ and two groups of model age (younger T_DM close to the apparent ages and older > 700 Ma), two possible pools, juvenile and inherited, were involved in mixing of mantle-derived magma with crustal components. The relatively high δ¹⁸O of 6.6 ± 0.3‰ for metamorphic zircon rims suggests either the protolith underwent hydrothermal alteration prior to the ~ 440 Ma oceanic crust subduction, or external higher δ¹⁸O fluid activities during UHP metamorphism at ~ 440 Ma.Therefore, the North Qaidam UHPM belt witnesses multiple tectonic evolution from Late Mesoproterozoic–Neoproterozoic assembly/breakup of the Rodinia supercontinent with related magmatic emplacement, then Paleozoic oceanic subduction, and finally transition of continental subduction/collision related to UHP metamorphism.     

The Meso-Neoarchaean Belomorian eclogite province: Tectonic position and geodynamic evolution

The aim of this paper is to review the main features of the Meso-Neoarchaean Belomorian eclogite province (BEP) in the northeastern Fennoscandian Shield, including regional and local geology, geochemistry, petrology and geochronology and to compare the Belomorian eclogites with Precambrian eclogites elsewhere. Two eclogite associations have been recognized within Belomorian TTG gneisses: (1) the subduction-type Salma association and (2) Gridino eclogitized mafic dykes. Protoliths of the Salma eclogites represent a sequence comprising gabbro, Fe–Ti gabbro and troctolites, formed at ~ 2.9 Ga in a slow-spreading ridge setting (like the Southwest Indian Ridge). The main subduction and eclogite-facies events occurred between ~ 2.87 and ~ 2.82 Ga. Injection of mafic magma into an active continental margin setting, recorded by the Gridino dyke swarm, is attributed to subduction of a mid-ocean ridge, commencing at 2.87 Ga. Crustal delamination of the active margin and subsequent involvement of the lower crust in subduction between 2.87 and 2.82 Ga ago caused high-pressure metamorphism of the Gridino dykes, culminating in eclogite-facies conditions between 2.82 and 2.78 Ga and accompanying amalgamation of the Karelia, Kola and Khetolamba blocks and formation of the Mesoarchaean Belomorian accretionary–collisional orogen. The clockwise P–T paths of the Salma and Gridino associations cross the granulite-facies P–T field. Detailed metamorphic studies indicate a complicated post-eclogite history with thermal events and fluid infiltration, related to plume activity at 2.72–2.70, ~ 2.4 and ~ 1.9 Ga. The eclogite assemblages were exhumed to mid-to-lower crustal depths at ~ 1.7 Ga, while erosion or younger tectonic events were responsible for final exhumation to the surface. Comparison of P–T–t paths and data for peak metamorphic parameters demonstrates the general similarity of the Archaean and Palaeoproterozoic eclogites worldwide and their association with anomalously “hot” environments. The occurrence of high-T conditions during eclogite-facies metamorphism can be attributed to either subduction of a mid-ocean ridge (Archaean, BEP) or to interaction with mantle plumes (Proterozoic).     

The multi-stage tectonic evolution of the Xitieshan terrane,North Qaidam orogen,western China: From Grenville-age orogeny to early-Paleozoic ultrahigh-pressure metamorphism

The geodynamic evolution of the early Paleozoic ultrahigh-pressure metamorphic belt in North Qaidam, western China, is controversial due to ambiguous interpretations concerning the nature and ages of the eclogitic protoliths. Within this framework, we present new LA-ICP-MS U–Pb zircon ages from eclogites and their country rock gneisses from the Xitieshan terrane, located in the central part of the North Qaidam UHP metamorphic belt. Xitieshan terrane contains clearly different protolith characteristics of eclogites and as such provides a natural laboratory to investigate the geodynamic evolution of the North Qaidam UHP metamorphic terrane. LA-ICP-MS U–Pb zircon dating of three phengite-bearing eclogites and two country rock gneiss samples from the Xitieshan terrane yielded 424–427 Ma and 917–920 Ma ages, respectively. The age of 424–427 Ma from eclogite probably reflects continental lithosphere subduction post-dating oceanic lithosphere subduction at ~ 440–460 Ma. The 0.91–0.92 Ga metamorphic ages from gneiss and associated metamorphic mineral assemblages are interpreted as evidence for the occurrence of a Grenville-age orogeny in the North Qaidam UHPM belt. Using internal microstructure, geochemistry and U–Pb ages of zircon in this study, combined with the petrological and geochemical investigations on the eclogites of previous literature’s data, three types of eclogitic protoliths are identified in the Xitieshan terrane i.e. 1) Subducted early Paleozoic oceanic crust (440–460 Ma), 2) Neoproterozoic oceanic crust material emplaced onto micro-continental fragments ahead of the main, early Paleozoic, collision event (440–420 Ma) and 3) Neoproterozoic mafic dikes intruded in continental fragments (rifted away from the former supercontinent Rodinia). These results demonstrate that the basement rocks of the North Qaidam terrane formed part of the former supercontinent Rodinia, attached to the Yangtze Craton and/or the Qinling microcontinent, and recorded a complex tectono-metamorphic evolution that involved Neoproterozoic and Early Paleozoic orogenies.