Remnant Paleoproterozoic Subduction or Lithospheric Drip Initiation at the Yilgarn Craton Margin:Constraints from P-wave Tomography

The Paleoproterozoic was a critical time in whether modern-style plate tectonics had become globally dominant(e.g.,Wan et al.,2020).The Capricorn Orogen witnessed the assembly of the Pilbara and Yilgarn Cratons and an exotic microcontinent,the Glenburgh Terrane,to form the West Australia Craton(WAC)through two collisional orogenic events,the 2215–2145 Ma Ophthalmian and 2005–1950 Ma Glenburgh Orogenies(Johnson et al.,2013;Fig.1).Compared to other Proterozoic orogenic belts in Australia,the Capricorn Orogen preserves‘complete'opposing continental margin successions,together with intervening arc fragments associated with oceanic closure and foreland basins associated with collisional loading(Cawood et al.,2009).

Deep structure of the Southeast Asian curved subduction system and its dynamic process

The deep structure,material circulation,and dynamic processes in the Southeast Asia have long been an elusive scientific puzzle due to the lack of systematic scientific observations and recognized theoretical models.Based on the deep seismic tomography using long-period natural earthquake data,in this study,the deep structure and material circulation of the curved subduction system in Southeast Asia was studied,and the dynamic processes since 100 million years ago was reconstructed.It is pointed out that challenges still exist in the precise reconstruction of deep mantle structures of the study area,the influence of multi-stage subduction on deep material exchange and shallow magma activity,as well as the spatiotemporal evolution and coupling mechanism of multi-plate convergence.Future work should focus on high-resolution land-sea joint 3-D seismic tomography imaging of the curved subduction system in the Southeast Asia,combined with geochemical analysis and geodynamic modelling works.

Coseismic Coulomb stress changes induced by a 2020-2021 M_(W)>7.0 Alaska earthquake sequence in and around the Shumagin gap and its influence on the Alaska-Aleutian subduction interface

Three M_(W)>7.0 earthquakes in 2020-2021 occurred in the Shumagin seismic gap and its adjacent area of the Alaska-Aleutian subduction zone,including the Mw7.8 Simeonof thrust earthquake on July 22,2020,the M_(W)7.6 Sand Point strike-slip earthquake on October 19,2020,and the M_(W)8.2 Chignik thrust earthquake on July 29,2021.The spatial and temporal proximity of these three earthquakes prompts us to probe stress-triggering effects among them.Here we examine the coseismic Coulomb stress change imparted by the three earthquakes and their influence on the subduction interface.Our results show that:(1)The Simeonof earthquake has strong loading effects on the subsequent Sand Point and Chignik earthquakes,with the Coulomb stress changes of 3.95 bars and 2.89 bars,respectively.The Coulomb stress change caused by the Sand Point earthquake at the hypocenter of the Chignik earthquake is merely around 0.01 bars,suggesting the negligible triggering effect on the latter earthquake;(2)The triggering effects of the Simeonof,Sand Point,and Chignik earthquakes on aftershocks within three months are not well pronounced because of the triggering rates of 38%,14%,and 43%respectively.Other factors may have played an important role in promoting the occurrence of these aftershocks,such as the roughness of the subduction interface,the complicated velocity structure of the lithosphere,and the heterogeneous prestress therein;(3)The three earthquakes caused remarkable coseismic Coulomb stress changes at the subduction interface nearby these mainshocks,with an average Coulomb stress change of 3.2 bars in the shallow region directly inwards the trench.

The Zoulang Nanshan Caledonian Subduct ion Complex in the Northern Qilian Mountains and Its Dynamics

The front of the Zoulang Nanshan Caledonian volcanic island arc zone in the northern Qilian Mountains is a forearc accretionary terrane, composed of multiple accretionary volcanic island arcs, flysch accretionary wedges,high-pressure metamorphosed detachment zones and remnants of ophiolites. It resulted from the northeastward subduction of the Early Palaeozoic Qilan oceanic crust beneath the Alxa block. High-pressure metamorphism, which occurred during the subduction, progressed through three stages: the initial stage of medium T-high P,the main stage of temperature decrease and pressure increase, and the lag stage of pressure decrease and temperature increase. Finally the paper presents a retrotrench subduction dynamic model indicative of northward subduction of the central Qilian block and southward accretion of the Alxa block during the period of 450-500 Ma.

Petrology and metamorphism of glaucophane eclogites in Changning-Menglian suture zone, Bangbing area, southeast Tibetan Plateau: An evidence for Paleo-Tethyan subduction

High/ultrahigh-pressure(HP/UHP)metamorphic complexes,such as eclogite and blueschist,are generally regarded as significant signature of paleo-subduction zones and paleo-suture zones.Glaucophane eclogites have been recently identified within the Lancang Group characterized by accretionary mélange in the Changning-Menglian suture zone,at Bangbing in the Shuangjiang area of southeastern Tibetan Plateau.The authors report the result of petrological,mineralogical and metamorphism investigations of these rocks,and discuss their tectonic implications.The eclogites are located within the Suyi blueschist belt and occur as tectonic lenses in coarse-grained garnet muscovite schists.The major mineral assemblage of the eclogites includes garnet,omphacite,glaucophane,phengite,clinozoisite and rutile.Eclogitic garnet contains numerous inclusions,such as omphacite,glaucophane,rutile,and quartz with radial cracks around.Glaucophane and clinozoisite in the matrix have apparent optical and compositional zonation.Four stages of metamorphic evolution can be determined:The prograde blueschist facies(M_(1)),the peak eclogite facies(M_(2)),the decompression blueschist facies(M_(3))and retrograde greenschist facies(M_(4)).Using the Grt-Omp-Phn geothermobarometer,a peak eclogite facies metamorphic P-T condition of 3000–3270 MPa and 617–658℃ was determined,which is typical of low-temperature ultrahigh-pressure metamorphism.The comparison of the geological characteristics of the Bangbing glaucophane eclogites and the Mengku lawsonite-bearing retrograde eclogites indicates that two suites of eclogites may have formed from significantly different depths or localities to create the tectonic mélange in a subduction channel during subduction of the Triassic Changning-Menglian Ocean.The discovery of the Bangbing glaucophane eclogites may represent a new oceanic HP/UHP metamorphic belt in the Changning-Menglian suture zone.

Lithium Isotopic Geochemistry in Subduction Zones:Retrospects and Prospects

Subduction zones involve many complex geological processes, including the release of slabderived fluids, fluid/rock interactions, partial melting, isotopic fractionations, elemental transporting, and crust/mantle interactions. Lithium（Li） isotopes（~6Li and ~7Li） have relative mass difference up to 16%, being the largest among metal elements. Thus, Li isotopes have advantage to interprete trace various geological processes. Most importantly, during crust/mantle interactions in deep subduction zones, surface materials and mantle rocks usually have distinct Li isotopic compositions. Li isotopes can be potential tracer for subduction processes, from the onset of subduction to the release of Li from subducted slabs and interaction with mantle wedge, as well as the fate of Li in slab-derived fluids and residual slabs. Moreover, the Li isotopic composition of subducting output materials can provide useful information for understanding global Li circulation. With developments in measurement and expansion of Li isotopic database, Li isotopic geochemistry will provide more inference and be a powerful tracer for understanding subduction-related processes. This work retrospected the application of Li isotopes in tracing successive subduction processes, and made some prospects for further studies of Li isotopes.

Enlightenment of the Mariana Fore-arc Sedimentary Basin Evolution to the Subduction Process

The Mariana subduction structure is a hot topic in ocean-ocean subduction zone research,and its subduction mechanism has attracted wide attention from experts and scholars in China and abroad.Based on the multi-channel seismic data of survey line MGL1204 in the Mariana fore-arc and DSDP ocean drilling data,this paper studies the development and evolution characteristics of the structure and strata in the Cenozoic Mariana fore-arc sedimentary basin.The Cenozoic strata are divided into six seismic sequences,with the possible era of each seismic sequence discerned,and the relationship between fault development and earthquakes analyzed.The episodic activity of the volcanic chain of the Mariana island arc is thought to control the tectonic and stratigraphic development pattern of the Cenozoic sedimentary basin in the fore-arc.Between 16°N-19°N and 146°E-151°E,the maximum thickness of the sedimentary center of the Cenozoic fore-arc sedimentary basin in Mariana is about 2360 m.Normal faults are developed in the area and some broke to the seabed,indicating that the Mariana island arc is still in the post-arc expansion stage.The application of multi-channel seismic sections in structural and stratigraphic evolution study is an important means to elucidating the Mariana subduction mechanism.

Topography of the 660-km discontinuity within the Izu-Bonin subduction zone and evidence of slab penetration near the Bonin Super Deep Earthquake(~680 km)

The Izu-Bonin subduction zone in the Northwest Pacific is an ideal location for understanding mantle dynamics such as cold lithosphere subduction. The slab produces a lateral thermal anomaly, inducing local topographic changes at the boundary of a post-spinel phase transformation, considered to be the origin of the ‘660-km discontinuity.’ In this study, the short-period(1–2 Hz) S-to-P conversion phase S660P was used to obtain the fine-scale structure of the discontinuity. More than 100 earthquakes that occurred from the 1980s to the 2020s and were recorded by high-quality seismic arrays in the United States and Europe were analyzed. A discontinuity in the ambient mantle with an average depth of ~670 km was found beneath the 300–400-km event zone in the northern Bonin region near 33°N. Meanwhile, the ‘660-km discontinuity’ has been pushed upward, away from the slab, possibly because of a hot upwelling mantle plume. In the central part of the subduction zone, the 660-km discontinuity is depressed to an average depth of(690 ± 5) km within the slab at approximately 150 km below the coldest slab core, indicating a(300 ± 100) ℃ cold anomaly estimated using a post-spinel transformation Clapeyron slope of(-2.0 ± 1.0) MPa/K. In southern Bonin near 28°N, the discontinuity was found to be further depressed at an average depth of(695 ± 5) km below the deepest event and with a focal depth of ~550 km. The discontinuity is located where the slab bends abruptly to become sub-horizontal toward the west-southwest. Near the zone of the isolated Bonin Super Deep Earthquake, which occurred at ~680 km on May 30,2015, the discontinuity is depressed to ~700 km, suggesting a near-vertical penetrating slab and an S-to-P conversion in the coldest slab core, where a large low-temperature anomaly should exist.

Tectono-magmatic evolution of Tethyan oceanic lithosphere in supra subduction zone fore arc regime:Geochemical fingerprints from crust-mantle sections of Naga Hills Ophiolite

The Naga Hills Ophiolite(NHO)belt in the Indo-Myanmar range(IMR)represents a segment of Tethyan oceanic crust and upper mantle that was involved in an eastward convergence and collision of the Indian Plate with the Burmese Plate during the Late Cretaceous-Eocene.Here,we present a detailed petrological and geochemical account for the mantle and crustal sections of NHO,northeastern India to address(i)the mantle processes and tectonic regimes involved in their genesis and(ii)their coherence in terms of the thermo-tectonic evolution of Tethyan oceanic crust and upper mantle.The NHO suite comprises well preserved crustal and mantle sections discretely exposed at Moki,Ziphu,Molen,Washelo and Lacham areas.The ultramafic-mafic lithologies of NHO are mineralogically composed of variable proportions of olivine,orthopyroxene,clinopyroxene and plagioclase.The primary igneous textures for the mantle peridotites have been overprinted by extensive serpentinisation whereas the crustal section rocks reflect crystal cumulation in a magma chamber.Chondrite normalised REE profiles for the cumulate peridotite-olivine gabbro-gabbro assemblage constituting the crustal section of NHO show flat to depleted LREE patterns consistent with their generation from depleted MORB-type precursor melt in an extensional tectonic setting,while the mantle peridotites depict U-shaped REE patterns marked by relative enrichment of LREE and HREE over MREE.These features collectively imply a dual role of depleted MORB-type and enriched arc-type mantle components for their genesis with imprints of melt-rock and fluid-rock interactions.Tectonically,studied lithologies from NHO correspond to a boninitic to slab-proximal Island Arc Tholeiite affinity thereby conforming to an intraoceanic supra subduction zone(SSZ)fore-arc regime coherent with the subduction initiation process.The geochemical attributes for the crustal and mantle sections of NHO as mirrored by Zr/Hf,Zr/Sm,Nb/Ta,Zr/Nb,Nb/U,Ba/Nb,Ba/Th,Ba/La and Nd/Hf ratios propound a two-stage petrogenetic process:(i)a depleted fore arc basalt(FAB)type tholeiitic melt parental to the crustal lithologies was extracted from the upwelling asthenospheric mantle at SSZ fore-arc extensional regime thereby rendering a refractory residual upper mantle;(ii)the crust and upper mantle of the SSZ fore arc were progressively refertilised by boninitic melts generated in response to subduction initiation and slab-dehydration.The vestiges of Tethyan oceanic lithosphere preserved in NHO represent an accreted intra-oceanic fore arc crust and upper mantle section which records a transitional geodynamic evolution in a SSZ regime marked by subduction initiation,fore arc extension and arc-continent accretion.

Ocean-Continent Subduction within the Paleotethyan Archiopeligic Ocean from Mineralogical Evidence of Muztag Ophiolite, Eastern Kunlun Mountain, China

Former studies show that the Muztag ophiolite, outcropped in the East Kunlun area of Xinjiang, formed in a supra-subduction zone environment. This study is to gain more information about the type of subduetion zone. Through field work, thin section observation and microprobe analysis, petrological and mineralogical characteristics of the metamorphic peridotites of this ophiolite are obtained. Although the olivines of metamorphic peridotites appear in three textures of metamorphic relict, metamorphic recrystallizations and orthopyroxene-melting crystallizations by thin-section observations, they have stable and low Fo range of 87.8- 89.5 by microprobe analysis. The orthopyroxenes show metamorphic relict and melting relict textures, with a low En of 88-90 and a wide range of Al2O3 content, from 2.90 wt% to 5. 13 wt%. The spinels develop anhedral-subhedral textures, with Cr^# （=Cr/（Cr＋AI）） focusing on two ranges of 0. 508-0. 723 and 0. 100-0. 118, respectively. Based on these petrological and mineralogical observations, and combined with the era and tectonic setting for the Muztag ophiolite, it can be concluded that the ophiolite formed in a supra-subduction zone where the oceanic crust subducted down to the continental are with a thick continental crust, and resulted from ocean-continent subducion within the Paleotethyan arehiopelagic ocean in the East Kunlun area of Xinjiang.

Phase Stability and Reactions of Subducting CaCO_(3)under Upper Mantle Conditions

CaCO_(3)is an important component of marine sediments and one of the major deep-carbon carriers at subduction zones.Some subducted CaCO_(3)can be dissolved in subduction fluids and recycled back to the surface via arc volcanoes degassing.At the same time,there still remain large amounts of CaCO_(3)and its reaction products,which could be further transported into Earth's deep interior.These internal processes link atmosphere,hydrosphere and biosphere with the deep solid Earth,modifying the environments of our planet.In this review,we summarize current understanding from high pressure-temperature experiments and field petrological observations on the physical and chemical properties of CaCO_(3).In particular,the phase stability and reactions of CaCO_(3)largely control the migration and reservation of oxidized carbon in subducting slabs.Finally,we present several critical but unsolved questions on CaCO_(3)subducting in the deep mantle.

Chemical composition of sediments from the subducting Cocos Ridge segment at the Southern Central American subduction zone

Subducted sediments play an important role in the magmatism at subduction zones and the formation of mantle heterogeneity,making them an important tracer for shallow crustal processes and deep mantle processes.Therefore,ascertaining the chemical compositions of different subduction end-members is a prerequisite for using subducted sediments to trace key geological processes.We reports here the comprehensive major and trace element analyses of 52 samples from two holes(U1414 A and U1381 C)drilled on the subducting Cocos Ridge segment at the Southern Central American(SCA)subduction zone during Integrated Ocean Drilling Program(IODP)Expedition 344.The results show that the SCA subducting sediments contain 51%(wt%)Ca CO_(3),27%(wt%)terrigenous material,16%(wt%)opal,and 6%(wt%)mineral-bound H2 O+.Compared to the global trenches subducting sediment,the SCA subducting sediments are enriched in biogenic elements(Ba,Sr,and Ca),and depleted in high field strength elements(Nb,Ta,Zr,Hf,and Ti)and alkali elements(K,Rb,and Cs).Meanwhile,the sediments in this area were affected by the carbonate crash event,which could have been caused by a~800 m rise in the carbonate compensation depth at 11 Ma in the Guatemala Basin.The reason for the sedimentary hiatus at Hole U1381 C may be the closure of the Panama Isthmus and the collision between the Cocos Ridge and the Middle America Trench.In addition,the sediments from the subducting Cocos Ridge segment have influenced the petrogenesis of volcanic lavas erupted in the SCA.

Subduction zone geochemistry

Crustal recycling at convergent plate boundaries is essential to mantle heterogeneity.However,crustal signatures in the mantle source of basaltic rocks above subduction zones were primarily incorporated in the form of liquid rather than solid phases.The physicochemical property of liquid phases is determined by the dehydration behavior of crustal rocks at the slab-mantle interface in subduction channels.Because of the significant fractionation in incompatible trace elements but the full inheritance in radiogenic isotopes relative to their crustal sources,the production of liquid phases is crucial to the geochemical transfer from the subducting crust into the mantle.In this process,the stability of specific minerals in subducting crustal rocks exerts a primary control on the enrichment of given trace elements in the liquid phases.For this reason,geochemically enriched oceanic basalts can be categorized into two types in terms of their trace element distribution patterns in the primitive mantle-normalized diagram.One is island arc basalts(IAB),showing enrichment in LILE,Pb and LREE but depletion in HFSE such as Nb and Ta relative to HREE,The other is ocean island basalts(OIB),exhibiting enrichment in LILE and LREE,enrichment or non-depletion in HFSE but depletion in Pb relative to HREE.In either types,these basalts show the enhanced enrichment of LILE and LREE with increasing their incompatibility relative to normal mid-ocean ridge basalts(MORB).The thermal regime of subduction zones can be categorized into two stages in both time and space,The first stage is characterized by compressional tectonism at low thermal gradients.As a consequence,metamorphic dehydration of the subducting crust prevails at forearc to subarc depths due to the breakdown of hydrous minerals such as mica and amphibole in the stability field of garnet and rutile,resulting in the liberation of aqueous solutions with the trace element composition that is considerably enriched in LILE,Pb and LREE but depleted in HFSE and HREE relative to normal MORB.This provides the crustal signature for the mantle sources of IAB.The second stage is indicated by extensional tectonism at high thermal gradients,leading to the partial melting of metamorphically dehydrated crustal rocks at subarc to postarc depths.This involves not only the breakdown of hydrous minerals such as amphibole,phengite and allanite in the stability field of garnet but also the dissolution of rutile into hydrous melts.As such,the hydrous melts can acquire the trace element composition that is significantly enriched in LILE,HFSE and LREE but depleted in Pb and HREE relative to normal MORB,providing the crustal signature for the mantle sources of OIB.In either case,these liquid phases would metasomatize the overlying mantle wedge peridotite at different depths,generating ultramafic metasomatites such as serpentinized and chloritized peridotites,and olivine-poor pyroxenites and hornblendites.As a consequence,the crustal signatures are transferred by the liquid phases from the subducting slab into the mantle.

Plate subduction, oxygen fugacity, and mineralization

Plate subduction is the largest natural factory that processes elements,which controls recycling and mineralization of a variety of elements.There are three major ore deposit belts in the world:the circumPacific,the centralAsian,and the Tethys belts.All the three belts are closely associated with plate subductions,the mechanism remains obscure.We approached this problem from systematic studies on the behaviours of elements during geologic processes.This contribution summaries the recent progress of our research group.Our results suggest that porphyry Cu deposits form through partial melting of subducted young oceanic crust under oxygen fugacities higher than AFMQ^+1.5,which is promoted after the elevation of atmospheric oxygen at ca.550 Ma.Tin deposits are associated with reducing magmatic rocks formed as a consequence of slab rollback.The Neo-Tethys tectonic regime hosts more than 60%of the world's total Sn reserves.This is due to the reducing environment formed during the subduction of organic rich sediments.For the same reason,porphyry Cu deposits formed in the late stages during the closure of the Neo-Tethys Ocean.Tungsten deposits are also controlled by slab rollback,but is not so sensitive to oxygen fugacity.Subduction related W/Sn deposits are mostly accompanied by abundant accessory fluorites due to the breakdown of phengite and apatite.Decomposition of phengite is also significant for hard rock lithium deposits,whereas orogenic belt resulted from plate subduction promote the formation of Li brine deposits.Cretaceous red bed basins near the Nanling region are favorable for Li brines.Both Mo and Re are enriched in the oxidationreduction cycle during surface processes,and may get further enriched once Mo-,Re-enriched sediments are subducted and involved in magmatism.During plate subduction,Mo and Re fractionate from each other.Molybdenum is mainly hosted in porphyry Mo deposits and to a less extent,porphyry Cu-Mo deposits,whereas Re is predominantly hosted in porphyry Cu-Mo deposits and sedimentary sulfide deposits.

Global kinematics of tectonic plates and subduction zones since the late Paleozoic Era

Detailed global plate motion models that provide a continuous description of plate boundaries through time are an effective tool for exploring processes both at and below the Earth's surface. A new generation of numerical models of mantle dynamics pre-and post-Pangea timeframes requires global kinematic descriptions with full plate reconstructions extending into the Paleozoic(410 Ma). Current plate models that cover Paleozoic times are characterised by large plate speeds and trench migration rates because they assume that lowermost mantle structures are rigid and fixed through time. When used as a surface boundary constraint in geodynamic models, these plate reconstructions do not accurately reproduce the present-day structure of the lowermost mantle. Building upon previous work, we present a global plate motion model with continuously closing plate boundaries ranging from the early Devonian at 410 Ma to present day.We analyse the model in terms of surface kinematics and predicted lower mantle structure. The magnitude of global plate speeds has been greatly reduced in our reconstruction by modifying the evolution of the synthetic Panthalassa oceanic plates, implementing a Paleozoic reference frame independent of any geodynamic assumptions, and implementing revised models for the Paleozoic evolution of North and South China and the closure of the Rheic Ocean. Paleozoic(410-250 Ma) RMS plate speeds are on average ～8 cm/yr, which is comparable to Mesozoic-Cenozoic rates of ～6 cm/yr on average.Paleozoic global median values of trench migration trend from higher speeds(～2.5 cm/yr) in the late Devonian to rates closer to 0 cm/yr at the end of the Permian(～250 Ma), and during the Mesozoic-Cenozoic(250-0 Ma) generally cluster tightly around ～1.1 cm/yr. Plate motions are best constrained over the past 130 Myr and calculations of global trench convergence rates over this period indicate median rates range between 3.2 cm/yr and 12.4 cm/yr with a present day median rate estimated at～5 cm/yr. For Paleozoic times(410-251 Ma) our model results in median convergence rates largely～5 cm/yr. Globally,～90% of subduction zones modelled in our reconstruction are determined to be in a convergent regime for the period of 120-0 Ma. Over the full span of the model, from 410 Ma to 0 Ma,～93% of subduction zones are calculated to be convergent, and at least 85% of subduction zones are converging for 97% of modelled times. Our changes improve global plate and trench kinematics since the late Paleozoic and our reconstructions of the lowermost mantle structure challenge the proposed fixity of lower mantle structures, suggesting that the eastern margin of the African LLSVP margin has moved by as much as ～1450 km since late Permian times(260 Ma). The model of the plate-mantle system we present suggests that during the Permian Period, South China was proximal to the eastern margin of the African LLSVP and not the western margin of the Pacific LLSVP as previous thought.

Zoned Zircon from Eclogite Lenses in Marbles from the Dabie-Sulu UHP Terrane,China: A Clear Record of Ultra-deep Subduction and Fast Exhumation

Eclogite lenses in marbles from the Dabie-Sulu ultrahigh-pressure （UHP） terrane are deeply subducted meta-sedimentary rocks. Zircons in these rocks have been used to constrain the ages of prograde and UHP metamorphism during subduction, and later retrograde metamorphism during exhumation. Inherited （detrital） and metamorphic zircons were distinguished on the basis of transmitted light microscopy, cathodoluminescence （CL） imaging, trace element contents and mineral inclusions. The distribution of mineral inclusions combined with CL imaging of the metamorphic zircon make it possible to relate zircon zones （domains） to different metamorphic stages. Domain 1 consists of rounded, oblong and spindly cores with dark-luminescent images, and contains quartz eclogite facies mineral inclusion assemblages, indicating formation under high-pressure （HP） metamorphic conditions of T = 571-668℃ and P =1.7-2.02 GPa. Domain 2 always surrounds domain 1 or occurs as rounded and spindly cores with white-luminescent images. It contains coesite eclogite facies mineral inclusion assemblages, indicating formation under UHP metamorphic conditions of T = 782-849℃ and P 〉 5.5 GPa. Domain 3, with gray-luminescent images, always surrounds domain 2 and occurs as the outermost zircon rim. It is characterized by low-pressure mineral inclusion assemblages, which are related to regional amphibolite facies retrograde metamorphism of T = 600- 710℃ and P = 0.7-1.2 GPa. The three metamorphic zircon domains have distinct ages; sample H1 from the Dabie terrane yielded SHRIMP ages of 245 ± 4 Ma for domain 1, 235 ± 3 Ma for domain 2 and 215± 6 Ma for domain 3, whereas sample H2 from the Sulu terrane yielded similar ages of 244 ± 4 Ma, 233 ± 4 Ma and 214 ± 5 Ma for Domains 1, 2 and 3, respectively. The mean ages of these zones suggest that subduction to UHP depths took place over 10-11 Ma and exhumation of the rocks occurred over a period of 19-20 Ma. Thus, subduction from - 55 km to 〉 160 km deep mantle depth took place at rates of approximately 9.5-10.5 km/Ma and exhumation from depths 〉160 km to the base of the crust at -30 km occurred at approximately 6.5 km/Ma. We propose a model for these rocks involving deep subduction of continental margin lithosphere followed by ultrafast exhumation driven by buoyancy forces after break-off of the UHP slab deep within the mantle.

The dominant driving force for supercontinent breakup: Plume push or subduction retreat?

Understanding the dominant force responsible for supercontinent breakup is crucial for establishing Earth＇s geodynamic evolution that includes supercontinent cycles and plate tectonics. Conventionally,two forces have been considered： the push by mantle plumes from the sub-continental mantle which is called the active force for breakup, and the dragging force from oceanic subduction retreat which is called the passive force for breakup. However, the relative importance of these two forces is unclear. Here we model the supercontinent breakup coupled with global mantle convection in order to address this question. Our global model features a spherical harmonic degree-2 structure, which includes a major subduction girdle and two large upwelling（superplume） systems. Based on this global mantle structure,we examine the distribution of extensional stress applied to the supercontinent by both subsupercontinent mantle upwellings and subduction retreat at the supercontinent peripheral. Our results show that：（1） at the center half of the supercontinent, plume push stress is ～3 times larger than the stress induced by subduction retreat;（2） an average hot anomaly of no higher than 50 K beneath the supercontinent can produce a push force strong enough to cause the initialization of supercontinent breakup;（3） the extensional stress induced by subduction retreat concentrates on a ～600 km wide zone on the boundary of the supercontinent, but has far less impact to the interior of the supercontinent. We therefore conclude that although circum-supercontinent subduction retreat assists supercontinent breakup, sub-supercontinent mantle upwelling is the essential force.

Geochemical cycling during subduction initiation:Evidence from serpentinized mantle wedge peridotite in the south Andaman ophiolite suite

The ophiolite suite from south Andaman Islands forms part of the Tethyan Ophiolite Belt and preserves the remnants of an ideal ophiolite sequence comprising a basal serpentinized and tectonised mantle peridotite followed by ultramafic and mafic cumulate units, basaltic dykes and spilitic pillow basalts interlayered with arkosic wacke. Here, we present new major, trace, rare earth(REE) and platinum group(PGE) element data for serpentinized and metasomatized peridotites(dunites) exposed in south Andaman representing the tectonized mantle section of the ophiolite suite. Geochemical features of the studied rocks, marked by Al_2 O_3/TiO_2 > 23, LILE-LREE enrichment, HFSE depletion, and U-shaped chondrite-normalized REE patterns with(La/Sm)N > 1 and(Gd/Yb)N <1, suggest contributions from boninitic mantle melts. These observations substantiate a subduction initiation process ensued by rapid slab roll-back with extension and seafloor spreading in an intraoceanic fore-arc regime. The boninitic composition of the serpentinized peridotites corroborate fluid and melt interaction with mantle manifested in terms of(i) hydration, metasomatism and serpentinization of depleted, MORB-type, sub-arc wedge mantle residual after repeated melt extraction; and(ii) refertilization of refractory mantle peridotite by boninitic melts derived at the initial stage of intraoceanic subduction. Serpentinized and metasomatized mantle dunites in this study record both MOR and intraoceanic arc signatures collectively suggesting suprasubduction zone affinity. The elevated abundances of Pd(4.4-12.2 ppb) with highΣPPGE/∑IPGE(2-3) and Pd/Ir(2-5.5) ratios are in accordance with extensive melt-rock interaction through percolation of boninitic melts enriched in fluid-fluxed LILE-LREE into the depleted mantle after multiple episodes of melt extraction. The high Pd contents with relatively lower Ir concentrations of the samples are analogous to characteristic PGE signatures of boninitic magmas and might have resulted by the infiltration of boninitic melts into the depleted and residual mantle wedge peridotite during fore-arc extension at the initial stage of intraoceanic subduction. The PGE patterns with high Os + Ir(2-8.6 ppb)and Ru(2.8-8.4 ppb) also suggest mantle rejuvenation by infiltration of melts derived by high degree of mantle melting. The trace, REE and PGE data presented in our study collectively reflect heterogeneous mantle compositions and provide insights into ocean-crust-mantle interaction and associated geochemical cycling within a suprasubduction zone regime.

Reproducing past subduction and mantle flow using high-resolution global convection models

Plate subduction drives both the internal convection and the surface geology of the solid Earth.Despite the rapid increase of computational power,it remains challenging for geodynamic models to reproduce the history of Earth-like subduction and associated mantle flow.Here,based on an adaptive approach of sequential data assimilation,we present a high-resolution global model since the mid-Mesozoic.This model incorporates the thermal structure and surface kinematics of tectonic plates based on a recent plate reconstruction to reproduce the observed subduction configuration and Earth-like convection.Introduction of temperature-and composition-dependent rheology allows for incorporation of many natural complexities,such as initiation of subduction zones,reversal of subduction polarity,and detailed plate-boundary dynamics.The resultant present-day slab geometry well matches Benioff zones and seismic tomography at depths < 1500 km,making it possible to hindcast past subduction dynamics and mantle flow.For example,the model produces a flat Farallon slab beneath North America during the Late Cretaceous to Early Cenozoic,a feature that has been geodynamically challenging to reproduce.This high-resolution model can also capture details of the 4-D evolution of slabs and the ambient mantle,such as temporally and spatially varying mantle flow associated with evolving slab geometry and buoyancy flux,as well as the formation of shallow slab tears due to subduction of young seafloors and the resulting complex mantle deformation.Such a geodynamic framework serves to further constrain uncertain plate reconstruction in the geological past,and to better understand the origin of enigmatic mantle seismic features.

Arc magmatism as a window to plate kinematics and subduction polarity:Example from the eastern Pontides belt,NE Turkey

The Eastern Pontides orogenic belt in the Black Sea region of Turkey offers a critical window to plate kinematics and subduction polarity during the closure of the Paleotethys. Here we provide a brief synthesis on recent information from this belt. We infer a southward subduction for the origin of the Eastern Pontides orogenic belt and its associated late Mesozoic--Cenozoic magmatism based on clear spatial and temporal variations in Late Cretaceous and Cenozoic arc magmatism, together with the exis- tence of a prominent south-dipping reverse fault system along the entire southern coast of the Black Sea. Our model is at variance with some recent proposals favoring a northward subduction polarity, and illus- trates the importance of arc magmatism in evaluating the geodynamic milieu associated with convergent margin orocesses.