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.

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.

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.

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.

Redox evolution of western Tianshan subduction zone and its effect on deep carbon cycle

Knowing the phase relations of carbon-bearing phases at high-pressure(HP) and high-temperature(HT) condition is essential for understanding the deep carbon cycle in the subduction zones.In particular,the phase relation of carbon-bearing phases is also strongly influenced by redox condition of subduction zones,which is poorly explored.Here we summarized the phase relations of carbon-bearing phases(calcite,aragonite,dolomite,magnesite,graphite,hydrocarbon) in HP metamorphic rocks(marble,metapelite,eclogite) from the Western Tianshan subduction zone and high-pressure experiments.During prograde progress of subduction,carbonates in altered oceanic crust change from Ca-carbonate(calcite) to Ca,Mg-carbonate(dolomite),then finally to Mgcarbonate(magnesite) via Mg-Ca cation exchange reaction between silicate and carbonate,while calcite in sedimentary calcareous ooze on oceanic crust directly transfers to high-pressure aragonite in marble or amorphous CaCO3 in subduction zones.Redox evolution also plays a significant effect on the carbon speciation in the Western Tianshan subduction zone.The prograde oxygen fugacity of the Western Tianshan subduction zone was constrained by mineral assemblage of garnet-omphacite from FMQ-1.9 to FMQ-2.5 at its metamorphic peak(maximum P-T) conditions.In comparison with redox conditions of other subduction zones,Western Tianshan has the lowest oxygen fugacity.Graphite and light hydrocarbon inclusions were ubiqutously identified in Western Tianshan HP metamorphic rocks and speculated to be formed from reduction of Fe-carbonate at low redox condition,which is also confirmed by high-pressure experimental simulation.Based on petrological observation and high-pressure simulation,a polarized redox model of reducing slab but oxidizing mantle wedge in subduction zone is proposed,and its effect on deep carbon cycle in subduction zones is further discussed.

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.

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.

Serpentinites of the Ural: Mineralogical Features, Petrophysical Properties and Subduction Processes

The article presents the results of study of composition, structure and properties of three genetic types of serpentinites from Ural: developed by chromite-bearing ultrabasic rocks, by metamorphic zones of carbonate rocks and within the zone of weathering of ultrabasic rocks. The samples were selected from deposits, located along the Main Ural Fault - the Paleozoic subduction zone (named GUR). Peculiarities of microstructure, chemical composition and properties of serpentinites formed in different geological conditions were investigated and their comparative study was held. Three groups of serpentinites with oriented and non-oriented structure, which formed by different protholites, were devided: mantle, lithosphere and crust.

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.

Seismic detection of the X-discontinuity beneath the Ryukyu subduction zone from the SdP conversion phase

The X-discontinuity,which appears at the depth of approximately 300 km,is an important seismic interface with positive velocity contrasts in the upper mantle.Detecting its presence and topography can be useful to understand phase transformations of relevant mantle minerals under the high-temperature and high-pressure circumstance of the Earth's interior.In this study,we detect the X-discontinuity beneath the Ryukyu subduction zone using five intermediate-depth events recorded by the dense Alaska Regional Network(AK).The X-discontinuity is successfully revealed from the robust slant stacking of the secondary down-going and converting Sd P phases.From the depth distribution of conversion points,we find that the X-discontinuity's depth ranges between 269 km and 313 km,with an average depth of 295 km.All the conversion points are located beneath the down-dipping side of the Philippine Sea slab.From energy comparisons in vespagrams for observed and synthetic seismograms,the strong converted energy is more likely from a thin high-velocity layer,and the S-wave velocity jumps across the X-discontinuity are up to 5% to 8% with an average of 6.0%.According to previous petrological and seismological studies,the X-discontinuity we detected can be interpreted as the phase transformation of coesite to stishovite in eclogitic materials within the oceanic crust.

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.

Thermal evolution of an ancient subduction interface revealed by Lue Hf garnet geochronology, Halilba■? Complex(Anatolia)

The thermal structure of subduction zones exerts a major influence on deep-seated mechanical and chemical processes controlling arc magmatism, seismicity, and global element cycles. Accretionary complexes exposed inland may comprise tectonic blocks with contrasting pressureetemperature(Pe T)histories, making it possible to investigate the dynamics and thermal evolution of former subduction interfaces. With this aim, we present new Lue Hf geochronological results for mafic rocks of the Halilbag?Complex(Anatolia) that evolved along different thermal gradients. Samples include a lawsoniteeepidote blueschist, a lawsoniteeepidote eclogite, and an epidote eclogite(all with counter-clockwise Pe T paths),a prograde lawsonite blueschist with a "hairpin"-type Pe T path, and a garnet amphibolite from the overlying sub-ophiolitic metamorphic sole. Equilibrium phase diagrams suggest that the garnet amphibolite formed at w0.6 -0.7 GPa and 800 -850℃, whereas the prograde lawsonite blueschist records burial from 2.1 GPa and 420℃ to 2.6 GPa and 520℃. Well-defined Lue Hf isochrons were obtained for the epidote eclogite(92.38 ± 0.22 Ma) and the lawsoniteeepidote blueschist(90.19 ± 0.54 Ma),suggesting rapid garnet growth. The lawsoniteeepidote eclogite(87.30 ± 0.39 Ma) and the prograde lawsonite blueschist(ca. 86 Ma) are younger, whereas the garnet amphibolite(104.5 ± 3.5 Ma) is older.Our data reveal a consistent trend of progressively decreasing geothermal gradient from granulite-facies conditions at ~104 Ma to the epidote-eclogite facies around 92 Ma, and the lawsonite blueschist-facies between 90 Ma and 86 Ma. Three Lue Hf garnet dates(between 92 Ma and 87 Ma) weighted toward the growth of post-peak rims(as indicated by Lu distribution in garnet) suggest that the HP/LT rocks were exhumed continuously and not episodically. We infer that HP/LT metamorphic rocks within the Halilbag?Complex were subjected to continuous return flow, with "warm" rocks being exhumed during the tectonic burial of "cold" ones. Our results, combined with regional geological constraints, allow us to speculate that subduction started at a transform fault near a mid-oceanic spreading centre. Following its formation, this ancient subduction interface evolved thermally over more than 15 Myr, most likely as a result of heat dissipation rather than crustal underplating.

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.

How many subductions in the Variscan orogeny?Insights from numerical models

We developed a 2 D numerical model to simulate the evolution of two superposed ocean-continent-ocean subduction cycles with opposite vergence,both followed by continental collision,aiming to better understand the evolution of the Variscan belt.Three models with different velocities of the first oceanic subduction have been implemented.Striking differences in the thermo-mechanical evolution between the first subduction,which activates in an unperturbed system,and the second subduction,characterised by an opposite vergence,have been enlighten,in particular regarding the temperature in the mantle wedge and in the interior of the slab.Pressure and temperature(P-T) conditions predicted by one cycle and two cycles models have been compared with natural P-T estimates of the Variscan metamorphism from the Alps and from the French Massif Central(FMC).The comparative analysis supports that a slow and hot subduction well reproduces the P-T conditions compatible with data from the FMC,while P-T conditions compatible with data of Variscan metamorphism from the Alps can be reproduced by either a cold or hot oceanic subduction models.Analysing the agreement of both double and single subduction models with natural P-T estimates,we observed that polycyclic models better describe the evolution of the Variscan orogeny.

Major Miocene geological events in southern Tibet and eastern Asia induced by the subduction of the Ninetyeast Ridge

Cenozoic adakitic rocks in the Gangdese changed from barren continental melts to ore-forming slab melts at * 23 Ma. The distribution and chemical characteristics of the ore-forming adakites point to an association with the Ninetyeast Ridge. The subduction of the thick,rigid Ninetyeast Ridge changed the geometry and rheology of the eastern Tibetan Plateau lithosphere and asthenosphere, restrained the eastward escape of asthenospheric mantle as well as continental fragments, and promoted the uplift and building of the Tibetan Plateau, which consequently changed the tectonic and climatic regimes in eastern Asia.

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.

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.

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.

New Zircon U-Pb Age of the Babu Ophiolites in Southeast Yunnan,China and Constrains of Plate Subduction Time

Objective The Babu ophiolite in Malipo County of southeastern Yunnan is interpreted as remanant ocean crust and represents a possible branch of Paleo-Tethyan Ocean in South China. It consists mainly of mafic and ultramafic rocks. These rocks are very important to understand the evolution of the Paleo-Tethyan Ocean. However, the Babu ophiolite is still disputed and the mafic and ultramafic rocks have been inferred to be part of the Emeishan large igneous province （LIP） by some researchers. In this paper, we present zircon U-Pb data on the metabasalts in Malipo to reveal the formation time of mafic and ultramafic rocks and their tectonic nature.