Experimental and Finite Element Analysis of Corroded High-Pressure Pipeline Repaired by Laminated Composite

Repairs of corroded high-pressure pipelines are essential for fluids transportation under high pressure.One of the methods used in their repairs is the use of layered composites.The composite used must have the necessary strength.Therefore,the experiments and analytical solutions presented in this paper are performed according to the relevant standards and codes,including ASME PCC-2,ASME B31.8S,ASME B31.4,ISO 24817 and ASME B31.G.In addition,the experimental tests are replicated numerically using the finite element method.Setting the strain gauges at different distances from the defect location,can reduce the nonlinear effects,deformation,and fluctuations due to the high pressure.The direct relationship between the depth of an axial defect and the stress concentration is observed at the inner side edges of the defect.Composite reparation reduces the non-linearities related to the sharp variation of the geometry and a more reliable numerical simulation could be performed.

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.

In situ characterization of forearc serpentinized peridotite from the Sulu ultrahigh-pressure terrane: Behavior of fluid-mobile elements in continental subduction zone

Serpentinized peridotites in the Yangkou(YK),Suoluoshu(SLS) and Hujialin(HJL) areas in the Sulu ultrahighpressure terrane represent the relic of ancient subcontinental lithospheric mantle below the North China Craton.Their protoliths,harzburgite and dunite,were variably hydrated by aqueous fluids released from subducting Yangtze continent.The rocks are enriched in fluid-mobile elements(FME) including Sb(42–333 times the depleted mantle value) and Pb(30–476 times).The degrees of the FME enrichment are comparable to that of the Himalayan forearc serpentinites,and greater than forearc mantle serpentinites from Marianas,suggesting that the degrees of FME enrichment in the forearc serpentinites are greater in continental subduction zones than those in the oceanic subduction zones.Lizardite after olivine in the SLS serpentinite shows higher degrees of enrichment in Sb and As than those for antigorite after both olivine and orthopyroxene in the YK area.The antigorite has highly enriched in Pb,U,Cs,and LREE,but not for the lizardite.The abundance of FME in two different species of serpentine reflects the different temperature of hydration.At temperature lower than 300 ℃,formed lizardite at shallow depths of the mantle wedge incorporates elements that are fluid mobile at low temperatures,such as Sb and As.When the temperature greater than 300 ℃,formed antigorite at a relatively deep mantle wedge incorporate more FME from the subducting continental slab(or fragments),including Pb,U,Cs,LREE as well as Sb and As.The eventual breakdown of antigorite(600–700 ℃) in prograde metamorphism would discharge water as well as FME into the subducting channel and/or the overlying mantle.

Metamorphic gradient modification in the Early Cretaceous Northern Andes subduction zone:A record from thermally overprinted high-pressure rocks

New field observations and petrological data from Early Cretaceous metamorphic rocks in the Central Cordillera of the Colombian Andes allowed the recognition of thermally overprinted high-pressure rocks derived from oceanic crust protoliths.The obtained metamorphic path suggests that the rocks evolved from blueschist to eclogite facies towards upper amphibolite to high-pressure granulite facies transitional conditions.Eclogite facies conditions,better recorded in mafic protoliths,are revealed by relic lawsonite and phengite,bleb-to worm-like diopside-albite symplectites,as well as garnet core composition.Upper amphibolite to high pressure granulite facies overprinting is supported by coarse-grained brown-colored Ti-rich amphibole,augite,and oligoclase recrystallization,as well as the record of partial melting leucosomes.Phase equilibria and pressure-temperature(P-T)path modeling suggest initial high-pressure metamorphic conditions M1 yielding 18.2-24.5 kbar and 465-580℃,followed by upper amphibolite to high pressure granulite facies overprinting stage M2 yielding 6.5-14.2 kbar and 580-720℃.Retrograde conditions M3 obtained through chlorite thermometry yield temperatures ranging around 286-400℃at pressures below 6.5-11 kbar.The obtained clockwise P-T path,the garnet zonation pattern revealing a decrease in X_(grs)/X_(prp)related to Mg#increment from core to rim,the presence of partial melting veins,as well as regional constraints,document the modification of the thermal structure of the active subduction zone in Northern Andes during the Early Cretaceous.Such increment of the metamorphic gradient within the subduction interface is associated with slab roll-back geodynamics where hot mantle inflow was triggered.This scenario is also argued by the reported trench-ward magmatic arc migration and multiple extensional basin formation during this period.The presented example constitutes the first report of Cretaceous roll-back-related metamorphism in the Caribbean and Andean realms,representing an additional piece of evidence for a margin-scale extensional event that modified the northwestern border of South America during the Early Cretaceous.

Fluid flow on centimeter-scale in deep paleosubduction zones in western Tianshan,China:Evidence from high-pressure veins

High-pressure(HP)veins were extensively developed in western Tianshan high-pressure(HP)metamorphic belt.The HP vein and host-rocks were analyzed by electronic microprobe to trace the origin of vein-forming fluids.Analytical data show that the immediately adjacent host-rocks of the studied HP vein are eclogites and gradually turned into blueschist as the distance from the veins increases,which indicates that the vein-forming fluid was derived from adjacent host-rocks;the boundaries between the vein and the host-rocks are sharp,which indicates that the fracture of the host-rocks is brittle during the vein-forming process.It is suggested that this type of HP veins is precipitated from the liquid formed by the dehydration of the host-rocks during the prograde metamorphism from blueschist to eclogite facies,which results in hydrofracturing of the rocks and provides the space for the vein to precipitate.The width of the eclogite-facies host-rocks is usually 1-2 cm,which provides the direct evidence that the fluid flow is on centimeter-scale.

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.

An assessment of the subduction rate in the CMIP6 historical experiment

Subduction process is a dynamical bridge for the exchanges of heat between the atmosphere and subsurface ocean water,which is regarded as a central proxy for the ocean climate studies.Given its key indicator in climate signals,it is of importance to examine the ability of a model to simulate the global subduction rate before investigating the climate dynamics.In this paper,we evaluated the ability of 21 climate models from Coupled Model Intercomparison Project Phase 6(CMIP6)in simulating the subduction rate.In general,the simulation ability of the models to the subduction climatology is better than that to the long-term variation trend.Based on the comprehensive analysis of climatology distribution and long-term trend of the subduction rate,GISS-E2-1-G performs better in reproducing the subduction rate climatology and IPSL-CM6A-LR can simulate positive long-term trend for both the global mean subduction rate and the lateral induction term in the Antarctic Circumpolar Current(ACC)region.However,it is still challenging to capture both the distribution characteristics of the subduction climatology and the long-term temporal trend for the 21 CMIP6 models.In addition,the model results demonstrate that,the ACC area is the major region contributing to the long-term trend of the global mean subduction rate.The analysis in this paper indicates that the poor simulation ability of reproducing the long-term trend of global mean subduction rate might be attributed to the ocean dynamics,for example,the zonal velocity at the bottom mixed layer and zonal gradient of mixed layer depth.

Aseismic ridge subduction and flat subduction:Insights from three-dimensional numerical models

Flat subduction can significantly influence the distribution of volcanism,stress state,and surface topography of the overriding plate.However,the mechanisms for inducing flat subduction remain controversial.Previous two-dimensional(2-D)numerical models and laboratory analogue models suggested that a buoyant impactor(aseismic ridge,oceanic plateau,or the like)may induce flat subduction.However,three-dimensional(3-D)systematic studies on the relationship between flat subduction and buoyant blocks are still lacking.Here,we use a 3-D numerical model to investigate the influence of the aseismic ridge,especially its width(which is difficult to consider in 2-D numerical models),on the formation of flat subduction.Our model results suggest that the aseismic ridge needs to be wide and thick enough to induce flat subduction,a condition that is difficult to satisfy on the Earth.We also find that the subduction of an aseismic ridge parallel to the trench or a double aseismic ridge normal to the trench has a similar effect on super-wide aseismic ridge subduction in terms of causing flat subduction,which can explain the flat subduction observed beneath regions such as Chile and Peru.

A new transition metal diphosphideα-MoP_(2) synthesized by a high-temperature and high-pressure technique

Monoclinicα-MoP_(2),with the OsGe2-type structure(space group C2/m,Z=4)and lattice parameters a=8.7248(11)Å,b=3.2322(4)Å,c=7.4724(9)Å,andβ=119.263°,was synthesized under a pressure of 4~GPa at a temperature between 1100℃and 1200℃.The structure ofα-MoP_(2) and its relationship to other transition metal diphosphides are discussed.Surprisingly,the ambient pressure phase orthorhombicβ-MoP_(2)(space group Cmc21)is denser in structure thanα-MoP_(2).Room-temperature high-pressure x-ray diffraction studies exclude the possibility of phase transition fromβ-MoP_(2) to α-MoP_(2),suggesting thatα-MoP_(2) is a stable phase at ambient conditions;this is also supported by the total energy and phonon calculations.

Characterization of subduction initiation

Compression is required for all kinds of subduction initiations,which may cause either subsidence or uplift,depending on the ages of the oceanic plates.Subduction initiations associated with the old oceanic crust tend to amplify preexisting subsidence by compression,whereas those associated with young oceanic plates may result in uplift.

Density Structure of the Papua New Guinea-Solomon Arc Subduction System

The Papua New Guinea-Solomon(PN-SL)arc is one of the regions with active crustal motions and strong geological actions.Thus,its complex subduction system makes it an ideal laboratory for studying the initiation mechanism of plate subduction.However,the PN-SL subduction system has not yet been sufficiently studied,and its density structure has yet to be revealed.In this paper,we used the free-air gravity data,Parker-Oldenburg density surface inversion method,and the genetic algorithm density inversion method to obtain the density structure of an approximately 1000-km-long northwest-southeast line crossing the PN-SL subduction system under the constraints of the CRUST1.0 global crustal model,onshore seismic data,and the LLNL-G3Dv3 global P-wave velocity model.The density structure shows that density differences between the plates on the two sides of the trench could play a significant role in plate subduction.

Numerical Simulation Analysis of the Transformer Fire Extinguishing Process with a High-Pressure Water Mist System under Different Conditions

To thoroughly study the extinguishing effect of a high-pressure water mist fire extinguishing system when a transformer fire occurs,a 3D experimental model of a transformer is established in this work by employing Fire Dynamics Simulator(FDS)software.More specifically,by setting different parameters,the process of the highpressure water mist fire extinguishing system with the presence of both diverse ambient temperatures and water mist sprinkler laying conditions is simulated.In addition,the fire extinguishing effect of the employed high-pressure water mist system with the implementation of different strategies is systematically analyzed.The extracted results show that a fire source farther away fromthe centerline leads to a lower local temperature distribution.In addition,as the ambient temperature increases,the temperature above the fire source decreases,while the temperature and the concentrationof theupperflue gas layer bothdecrease.Interestingly,after thehigh-pressurewatermist sprinkler begins to operate,both the temperature distribution above the fire source and the concentration of the flue gas decrease,which indicates that the high-pressure water mist system plays the role of cooling and dust removal.By comparing various sprinkler laying methods,it is found that the lower sprinkler height has a better effect on the temperature above the fire source,the temperature of the upper flue gas layer,and the concentration of the flue gas.Moreover,when the sprinkler is spread over thewhole transformer,the cooling effect on both the temperature above the fire source and the temperature of the upper flue gas layer is good,whereas the change in the concentration of the flue gas above the fire source is not obvious compared to the case where the sprinkler is not fully spread.

High-pressure new phases of V–N compounds

The high-pressure diagram of V–N compounds is enriched by proposed seven new stable high-pressure phases.The P-1-VN_4with the armchair N-rich structure may be quenched to ambient conditions.The formed N–N covalent bond plays an important role for the structural stability of N-chain.The charge transfer results in a V–N ionic bond interaction,which further improves the stability of N-chain structure.The P-1-VN_4,P4mnc-VN_8,and Immm-VN_(10)with the outstanding detonation properties have potential application in explosive field.

The discovery of TiO_(2)-Ⅱ,the α-PbO_(2)-structured high-pressure polymorph of rutile,in the Suizhou L6 chondrite

We report the discovery of TiO_(2)-Ⅱ in the unmelted rock of the shocked Suizhou L6 chondrite.Natural TiO_(2)-Ⅱ was previously found in ultrahigh-pressure metamorphic and mantle-derived rocks,terrestrial impact structures,and tektite.Our microscopic,Raman spectroscopic,electron microprobe and transmission electron microscopic investigations have revealed:(1) All observed TiO_(2)-Ⅱ grains are related with ilmenite and pyrophanite;(2) TiO_(2)-Ⅱ occurs as needle-and leaf-shaped inclusions in llmenite and patch-,tape-shaped body in pyrophanite;(3)The composition of TiO_(2)-Ⅱ is identical with that of its precursor rutile;(4) The Raman spectrum of TiO_(2)-Ⅱ is in good agreement with that of natural and synthesized α-PbO_(2)-type TiO_(2);(5) TiO_(2)-Ⅱ occurs mainly in the form of well-ordered nano-domains and small mis-orientation among the domains can be observed.(6) All electron diffraction reflections from TiO_(2)-Ⅱ can be indexed to α-PbO_(2)structure in space group Pbcn with lattice parameters of a=4.481 ?,b=5.578 A and c=4.921 A;(7) The exsolution inclusions of rutile from host ilmenite are mostly connected with an alternation process along the lamellar twinning plane of ilmenite induced by shockinduced high pressure and high temperature;(8) The P-T regime of 20-25 GPa and 1000 ℃ estimated for the Suizhou unmelted rock is suitable for phase transition of rutile into TiO_(2)-Ⅱ phase.

Magmatism-Related Thermal Simulation of Volcanic Arcs in the Molucca Sea Bidirectional Subduction System

The bidirectional subduction system,island arc magmatic activities,and thermal structure of the forearc basin in the Molucca Sea are taken into consideration in this study.The active volcanic arcs on both sides of the bidirectional subduction zone in the Molucca Sea are undergoing arc-arc collisions.We applied a finite element thermal simulation method to reconstruct the thermal evolution history of the Molucca Sea Plate based on geophysical data.Then,we analyzed the thermodynamic characteristics of island arc volcanism on both sides of the bidirectional subduction zone.The results showed that at 10Myr,the oceanic ridge of the Molucca Sea Plate was asymmetrically biased to the west,causing this bidirectional subduction to be deeper in the west than in the east.Furthermore,the oceanic ridge subducted under the Sangihe arc at 5.5Myr,causing intermittent cessation of volcanic activities.Due to the convergence of bidirectional subduction,the geothermal gradient in the top 3km depth of the forearc area between the Sangihe and Halmahera arcs decreased from about 60℃km^(−1) at 4Myr to about 38℃km^(−1) today.Finally,within the 45–100 km depth range of the sliding surface of the subduction,anomalously high-temperature zones formed due to shear friction during the bidirectional subduction.

Carbonatite Occurrence in Ambaji-Sendra Belt of NW Indian Shield: Evidence of Carbonatitic Magmatism in the Subduction Setting

The paper represents a new discovery of a late Mesoproterozoic lenticular and discontinuous, carbonatitic body exposed at Basantgarh, Sendra and near the Abu-road area of the Ambaji-Sendra belt of the South Delhi Fold Belt. It is medium to coarse-grained and light to dark coloured compact rock. The common associates of the carbonatitic rock are felsic rocks, rich in alkalies. Carbonatite contains more than 50% carbonate minerals, the majority of which are calcite, dolomite, ankerite, augite-aegirine augite and plagioclase. It is classified as calcite carbonatite of the sövite variety due to its coarse-grained character, chemically as calico-carbonatite and magnesio-carbonatite and even as silico-carbonatites for having more than 20% SiO2. The ∑REE contents of calico-carbonatite samples are nearly 100 times greater than magnesio-carbonatite. Chondrite normalised REE profiles of calcio-carbonatites are LREE enriched with nearly flat HREEs whereas the magnesio-carbonatite is characterised by flat REE patterns. The mantle-normalized incompatible trace element spidergram of Ambaji-Sendra belt carbonatites shows distinct negative anomalies of Ba, Nb, Ta, P, Sm, Eu, Ti and Y and positive at U and Pb by calcio-carbonatite whereas the magnesio-carbonatite displays negative kinks at K, Zr, Nb, Ta and Ti and positive at Th, Pb and Sr. The variable and/or contrasting enrichment/depletion in various elements in the two types of Ambaji-Sendra belt carbonatite is attributed either to significant differences in the type and modal proportion of different accessory mineral species or selective incorporation of metasomatic fluids during the subduction process. The chemical attributes of Ambaji-Sendra belt carbonatite suggest its emplacement in subduction settings.

Coupling between the Cenozoic west Pacific subduction initiation and decreases of atmospheric carbon dioxides

At the beginning of the Cenozoic,the atmospheric CO_(2)concentration increased rapidly from~2000 ppmv at 60 Ma to~4600 ppmv at 51 Ma,which is 5–10 times higher than the present value,and then continuous declined from~51 to 34 Ma.The cause of this phenomenon is still not well understood.In this study,we demonstrate that the initiation of Cenozoic west Pacific plate subduction,triggered by the hard collision in the Tibetan Plateau,occurred at approximately 51 Ma,coinciding with the tipping point.The water depths of the Pacific subduction zones are mostly below the carbonate compensation depths,while those of the Neo-Tethys were much shallower before the collision and caused far more carbonate subducting.Additionally,more volcanic ashes erupted from the west Pacific subduction zones,which consume CO_(2).The average annual west Pacific volvano eruption is 1.11 km~3,which is higher than previous estimations.The amount of annual CO_(2)absorbed by chemical weathering of additional west Pacific volcanic ashes could be comparable to the silicate weathering by the global river.We propose that the initiation of the western Pacific subduction controlled the long-term reduction of atmospheric CO_(2)concentration.

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.

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.