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.

Dynamic subduction process of local plate revealed by Ibaraki earthquake sequence of 1982 in Japan

The kinematics and dynamics of plate tectonics are frontal subjects in geosciences and the strong earthquake occurred along the plate boundary result directly from plate movement. By analyzing Ibaraki earthquake sequence, it has been found that the focal fault plane shows a special image of grading expansion along the direction of strike and adjustment along the dip direction respectively. With the consideration of strike, dip and slip directions of focal mechanism, we have confirmed that Ibaraki earthquake belongs to a thrust fault earthquake occurred under the Japan Trench. The cause of the earthquake sequence is discussed in the paper. The study on the temporal-spatial distribution of the earthquake sequence with a time-scale between the year-scale spatial geodetic data and the second-scale moment tensor of the strong earthquake has indicated the dynamic process of Pacific Plate subduction under the Eurasia Plate. According to the average slip distance of earthquake and the velocity of plate movement, it is predicted that a strong earthquake might occur in recent years.

Reviewing subduction initiation and the origin of plate tectonics:What do we learn from present-day Earth?

The theory of plate tectonics came together in the 1960s,achieving wide acceptance after 1968.Since then it has been the most successful framework for investigations of Earth’s evolution.Subduction of the oceanic lithosphere,as the engine that drives plate tectonics,has played a key role in the theory.However,one of the biggest unanswered questions in Earth science is how the first subduction was initiated,and hence how plate tectonics began.The main challenge is how the strong lithosphere could break and bend if plate tectonics-related weakness and slab-pull force were both absent.In this work we review state-of-the-art subduction initiation(SI)models with a focus on their prerequisites and related driving mechanisms.We note that the plume-lithosphere-interaction and mantleconvection models do not rely on the operation of existing plate tectonics and thus may be capable of explaining the first SI.Reinvestigation of plate-driving mechanisms reveals that mantle drag may be the missing driving force for surface plates,capable of triggering initiation of the first subduction.We propose a composite driving mechanism,suggesting that plate tectonics may be driven by both subducting slabs and convection currents in the mantle.We also discuss and try to answer the following question:Why has plate tectonics been observed only on Earth?

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.

Formation of the adakite-like granitoid complex and porphyry copper-gold deposit in Shaxi from southern Tancheng-Lujiang fault belt: A clue to the West Pacific plate subduction

On the basis of the geological and geochemical studies, including chemical analysis of bulk rocks, rare-earth and trace element studies, fluid inclusion, and S and O isotopic analyses, the authors described the geological background of the deposit in detail and presented significant proofs for the conditions of formation of the Shaxi porphyry copper-gold deposit. Compared with other large and supper-large porphyry copper deposits in China and the adjacent Cu-Au mineralized areas, the ore-forming processes and conditions were analyzed; and the possibility of forming large porphyry copper deposits in the Shaxi area was discussed. The present study indicated that the ore-forming fluid and material were mainly of magmatic origin, while meteoric water played a certain role in the ore-forming processes. Interactions between subducting and overriding plates provided a major driving force for the formation of igneous rocks and the deposition of metal elements in East China since Jurassic. Based on the geo- chemical data of the Shaxi intrusive, it is found that the copper (gold) mineralization is closely related to the genesis of adakite-like intrusive in the Shaxi area. This adakite-like intrusive was formed in the subduction environment as a result of the subduction of the West Pacific plate toward the East China continent, where there is a great potentiality to form a large porphyry copper deposit.

Characteristics of Paleoproterozoic Subduction System in Western Margin of Yangtze Plate

Paleoproterozoic subduction strongly occurred in the western margin of Yangtze plate. The basalticandesite volcanics of Ailaoshan Group and Dibadu Formation had been formed during paleo QinghaiTibet oceanic plate subduction under the paleoYangtze plate. Their trace element geochemistry suggests that their forming environments are continentalmarginarc and back arcbasin respectively. Consequently, the Paleoproterozoic subduction system in the western margin of Yangtze plate was established. Ailaoshan Group and Dibadu Formation came from an enriched mantle source that was contaminated by crustal sediments carried by subducted slab, and formed the Paleoroterozoic metamorphic basement of western margin of Yangtze plate. Ailaoshan Group is actually western boundary of Yangtze plate.

Two-dimensional Numerical Modeling Research on Continent Subduction Dynamics

Continent subduction is one of the hot research problems in geoscience. New models presented here have been set up and two-dimensional numerical modeling research on the possibility of continental subduction has been made with the finite element software, ANSYS, based on documentary evidence and reasonable assumptions that the subduction of oceanic crust has occurred, the subduction of continental crust can take place and the process can be simplified to a discontinuous plane strain theory model. The modeling results show that it is completely possible for continental crust to be subducted to a depth of 120 km under certain circumstances and conditions. At the same time, the simulations of continental subduction under a single dynamical factor have also been made, including the pull force of the subducted oceanic lithosphere, the drag force connected with mantle convection and the push force of the mid-ocean ridge. These experiments show that the drag force connected with mantle convection is critical for continent subduction.

Why primordial continents were recycled to the deep:Role of subduction erosion

Geological observations indicate that there are only a few rocks of Archean Earth and no Hadean rocks on the surface of the present-day Earth.From these facts,many scientists believe that the primordial continents never existed during Hadean Earth,and the continental volume has kept increasing.On the other hand,recent studies reported the importance of the primordial continents on the origin of life,implying their existence.In this paper,we discussed the possible process that could explain the loss of the primordial continents with the assumption that they existed in the Hadean.Although depending on the timing of the initiation of plate tectonics and its convection style,subduction erosion,which is observed on the present-day Earth,might have carried the primordial continents into the deep mantle.

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.

Decadal Seismicity Prior to Great Earthquakes at Subduction Zones: Roles of Major Asperities and Low-Coupling Zones

Decadal forerunning seismic activity of magnitude Mw ≥ 5.0 is mapped for all 45 mainshocks of Mw 7.7 to 9.1 at subduction zones of the world from 1993 to mid 2020. The zones of high slip in nearly all great earthquakes were nearly quiescent beforehand and are identified as the sites of great asperities and zones of strong seismic coupling. Much forerunning activity occurred at smaller asperities along the peripheries of the rupture zones of many great and giant mainshocks. Those sizes of great asperities as ascertained from forerunning activity generally agree with the areas of high seismic slip as determined by others from geodetic and tide-gauge data and finite-source seismic modeling. Asperities are strong, well-coupled portions of plate interfaces. Different patterns of forerunning activity on time scales of about 5 to 45 years are attributed to either the sizes and spacing of asperities (or lack of). This permits many great asperities to be mapped decades before they rupture in great and giant shocks. Several poorly coupled subduction zones such as Java, Lesser Sunda, Marianas, Tonga and Kermadec are characterized by few great thrust earthquakes and little, in any forerunning activity. Rupture zones of many great and giant earthquakes are bordered either along strike, updip, or downdip by zones of lower plate coupling. Several bordering regions were sites of forerunning activity, aftershocks, and slow-slip events. The detection of forerunning and precursory activities of various kinds should be sought on the peripheries of great asperities as well as within zones of high co-seismic slip.

A-type granites induced by a breaking-off and delamination of the subducted Junggar oceanic plate,West Junggar,Northwest China

The A-type granites with highly positiveε_(Nd)(t)values in the West Junggar,Central Asian Orogenic Belt(CAOB),have long been perceived as a group formed under the same tectonic and geodynamic setting,magmatic sourceq and petrogenetic model.Geological evidence shows that these granites occurred at two different tectonic units related to the southeastern subduction of Junggar oceanic plate:the Hongshan and Karamay granites emplaced in the southeast of West Junggar in the Baogutu continental arc;whereas the Akebasitao and Miaoergou granites formed in the accretionary prism.Here the authors present new bulk-rock geochemistry and Sr-Nd isotopes,zircon U-Pb ages and Hf-O isotopes data on these granites.The granites in the Baogutu continental arc and accretionary prism contain similar zirconε_(Hf)(t)values(+10.9 to+16.2)and bulk-rock geochemical characteristics(high SiO_(2)and K_(2)O contents,enriched LILEs(except Sr),depleted Sr,Ta and Ti,and negative anomalies in Ce and Eu).The Hongshan and Karamay granites in the Baogutu continental arc have older zircon U-Pb ages(315-305 Ma)and moderate^(18)O enrichments(δ^(18)_(O_(zircon))=+6.41‰-+7.96‰);whereas the Akebasitao and Miaoergou granites in the accretionary prism have younger zircon U-Pb ages(305-301 Ma)with higher^(18)O enrichments(δ^(18)_(O_(zircon))=+8.72‰-+9.89‰).The authors deduce that the elevated^(18)O enrichments of the Akebasitao and Miaoergou granites were probably inherited from low-temperature altered oceanic crusts.The Akebasitao and Miaoergou granites were originated from partial melting of low-temperature altered oceanic crusts with juvenile oceanic sediments below the accretionary prism.The Hongshan and Karamay granites were mainly derived from partial melting of basaltic juvenile lower crust with mixtures of potentially chemical weathered ancient crustal residues and mantle basaltic melt(induced by hot intruding mantle basaltic magma at the bottom of the Baogutu continental arc).On the other hand,the Miaoergou charnockite might be sourced from a deeper partial melting reservoir under the accretionary prism,consisting of the low-temperature altered oceanic crust,juvenile oceanic sediments,and mantle basaltic melt.These granites could be related to the asthenosphere's counterflow and upwelling,caused by the break-off and delamination of the subducted oceanic plate beneath the accretionary prism Baogutu continental arc in a post-collisional tectonic setting.

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.

Effective Elastic Thickness of the Lithosphere in the Mariana Subduction Zone and Surrounding Regions and Its Implications for Their Tectonics

To understand the rheology,structure,and tectonics of the lithosphere in the Mariana subduction zone and surrounding regions,we calculated the effective elastic thickness of the lithosphere(Te)in these areas using the improved moving window admittance technique(MWAT)method.We find that smaller data grid spacing can better reflect Te variations in the subduction zone.The Te of the study region ranges from 0 to 47 km.The Te is reduced from 40 km on the seaward side of the outer-rise region to 1-2 km along the trench axis.The lithospheric breaking distance from the trench axis ranges from 0 to 250 km.We suggest that the intermediate Te values in seamounts and high Te values on the seaward side of the outer-rise region respectively reflect the‘fossil’rheological state and current lithospheric strength of the Pacific plate.The faulting induced by the downward bending of subducting plate not only ruptures the lithosphere but also contributes to the mantle serpentinization,significantly reducing the lithospheric strength.The largest breaking distance of the Ogasawara Plateau may be due to the increase in the mass load of the subducting plate in the Ogasawara Plateau and the significant horizontal bending force in the plate caused by the resistance of seamounts to subduction.Furthermore,a good positive correlation exists between the breaking distance and subduction dip angle along the trench axis.We suggest that the subducting plate with a larger breaking distance is likely to form a larger subduction angle.

Seismic-gravimetric analysis of the subducted Nazca plate 1 between 32°S and 36°S

The study region is seismically and tectonically characterized by the angle variations in the subduction of the Nazca plate. The results obtained from earthquakes location between 32° and 36°S latitude and67°-71°W longitude are presented in this work. The presence of a wedge of asthenospheric materials and the partial or total eclogitization of the subducted Nazca plate and its relation with isostatic cortex models published was analyzed. In addition, a gravimetric profile obtained from gravity forward modeling is presented at 33.5°S, proposing a new configuration at depths for the main tectonic components： Nazca plate, asthenospheric wedge and South American plate. Also, a new density scheme using recently published velocity models was obtained.

Ripple Tectonics—When Subduction Is Interrupted

Subduction plays a fundamental role in plate tectonics and is a significant factor in modifying the structure and topography of the Earth. It is driven by convection forces that change over a >100 Myr time scale. However, when an oceanic plateau approaches, it plugs the subduction, and causes slab necking and tearing. This abrupt change may trigger a series of geodynamic (tectonic, volcanic) and sedimentary responses recorded across the convergence boundary and its surrounding regions by synchronous structural modifications. We suggest that a large enough triggering event may lead to a ripple tectonic effect that propagates outwards while speeding up the yielding of localized stress states that otherwise would not reach their threshold. The ripple effect facilitates tectonic, volcanic, and structural events worldwide that are seemingly unrelated. When the world’s largest oceanic plateau, Ontong Java Plateau (OJP), choked the Pacific-Australian convergence zone at ~6 Myr ago, it induced kinematic modifications throughout the Pacific region and along its plate margins. Other, seemingly unrelated, short-lived modifications were recorded worldwide during that time window. These modifications changed the rotation of the entire Pacific plate, which occupies ~20% of the Earth’s surface. In addition, the Scotia Sea spreading stopped, global volcanism increased, the Strait of Gibraltar closed, and the Mediterranean Sea dried up and induced the Messinian salinity crisis. In this paper, we attribute these and many other synchronous events to a new “ripple tectonics” mechanism. We suggest that the OJPincipient collision triggered the Miocene-Pliocene transition. Similarly, we suggest that innovative GPS-based studies conducted today may seek the connectivity between tectonic, seismic, and volcanic events worldwide.

A new tectonic model for the Itmurundy Zone,central Kazakhstan:linking ocean plate stratigraphy,timing of accretion and subduction polarity

The Itmurundy Zone of Central Kazakhstan is a key structure in the core of the Kazakh Orocline representing a typical Pacific-type orogenic belt hosting accretionary complex,ophiolite massifs and serpentinite mélange.The main controversies in the existing tectonic models of the Itmurundy Zone are about the timing of subduction and accretion,the direction and kinematics of subduction and the number of oceanic plates.A new model for the early Paleozoic tectonic story of the Itmurundy Zone is postulated in this paper,based on new detailed geological and U-Pb detrital zircon age data,combined with previously documented geological,U-Pb age,microfossil,geochemical and isotope data from igneous rocks,deep-sea sediments and greywacke sandstones.The present study employs the Ocean Plate Stratigraphy(OPS)model to explain the tectonic processes involved in the evolution of the Itmurundy Zone and to pre-sent a holistic story of Ordovician oceanic plate(s),which accretion formed an accretionary complex.The detailed mapping allows distinguishing three types of OPS assemblages:(1)Chert-dominated,(2)OIB-hosting,and(3)MORB-hosting.The U-Pb ages of detrital zircons from sandstones of OIB and Chert types show unimodal distributions with similar main peaks of magmatism at 460-455 Ma in the provenance,and their maximum depositional ages(MDA)span 455-433 Ma.Two samples from OPS Type 3 show the peaks of magmatism both at ca.460 Ma and the MDA of 452 Ma and 459 Ma,respectively.The MDA of sandstones and microfossils data from chert show the younging of strata to the south and SE in Types 1 and 2 and to NEE for Type 3(in present coordinates)suggesting double-sided subduction to the NNW and SEE and,accordingly,the co-existence of pieces of two oceanic plates in Ordovician time.The U-Pb zircon data from both igneous and clastic rocks indicate a period of subduction erosion in early Ordovician time.As a whole,the accreted OPS units of the Itmurundy Zone record the timing of subduction and accretion from the early Ordovician to the early Silurian,i.e.,60 Ma at shortest.

Early Cenozoic Tectonics of the Tibetan Plateau

Geological mapping at a scale of 1：250000 coupled with related researches in recent years reveal well Early Cenozoic paleo-tectonic evolution of the Tibetan Plateau. Marine deposits and foraminifera assemblages indicate that the Tethys-Himalaya Ocean and the Southwest Tarim Sea existed in the south and north of the Tibetan Plateau, respectively, in Paleocene-Eocene. The paleo- oceanic plate between the Indian continental plate and the Lhasa block had been as wide as 900km at beginning of the Cenozoic Era. Late Paleocene transgressions of the paleo-sea led to the formation of paleo-bays in the southern Lhasa block. Northward subduction of the Tethys-Himalaya Oceanic Plate caused magma emplacement and volcanic eruptions of the Linzizong Group in 64.5-44.3 Ma, which formed the Paleocene-Eocene Gangdise Magmatic Arc in the north of Yalung-Zangbu Suture （YZS）, accompanied by intensive thrust in the Lhasa, Qiangtang, Hoh Xil and Kunlun blocks. The Paleocene- Eocene depression of basins reached to a depth of 3500-4800 m along major thrust faults and 680-850 m along the boundary normal faults in central Tibetan Plateau, and the Paleocene-Eocene depression of the Tarim and Qaidam basins without evident contractions were only as deep as 300-580 m and 600-830 m, respectively, far away from central Tibetan Plateau. Low elevation plains formed in the southern continental margin of the Tethy-Himalaya Ocean, the central Tibet and the Tarim basin in Paleocene-Early Eocene. The Tibetan Plateau and Himalaya Mts. mainly uplifted after the Indian- Eurasian continental collision in Early-Middle Eocene.

Plate convergence in the Indo-Pacific region

The Indo-Pacific convergence region is the best target to solve the teo remaining challenges of the plate tectonics theory,i.e.,subduction initiation and the driving force of plate tectonics.Recent studies proposed that the Izu-Bonin subduction initiation belongs to spontaneous initiation,which implies that it started from extension,followed by low angle subduction.Numerical geodynamic modeling suggests that the initiation of plate subduction likely occurred along a transform fault,which put the young spreading ridge in direct contact with old oceanic crust.This,however,does not explain the simultaneous subduction initiation in the west Pacific region in the Cenozoic.Namely,the subduction initiations in the Izu-BoninMariana,the Aleutian,and the Tonga-Kermadec trenches are associated with oceanic crusts of different ages,yet they occurred at roughly the same time,suggesting that they were all triggered by a maj or change in the Pacific plate.Moreover,low angle subduction induces compression rather than extension,which requires external compression forces.Given that the famous Hawaiian-Emperor bending occurred roughly at the same time with the onset of westward subductions in the west Pacific,we propose that these Cenozoic subductions were initiated by the steering of the Pacific plate,which are classified as induced initiation.Induced subduction initiation usually occurs in young ocean basins,forming single-track subduction.The closure s of Neo-Tethys Oceans were likely triggered by plume s in the south,forming northward subductions.Interestingly,the Indian plate kept on moving northward more than 50 Ma after the collision between the Indian and Eurasian continents and the break-off of the subducted oceanic slab attached to it.This strongly suggests that slab pull is not the main driving force of plate tectonics,whereas slab sliding is.

Early Cretaceous Tectonics and Evolution of the Tibetan Plateau

Selected geological data on Early Cretaceous strata, structures, magmatic plutons and volcanic rocks from the Kunlun to Himalaya Mountains reveal a new view of the Early Cretaceous paleo-tectonics and the related geodynamic movement of the Tibetan Plateau. Two major paleo- oceans, the Mid-Tethys Ocean between the Qiangtang and Lhasa blocks, and the Neo-Tethys Ocean between the Lhasa and Himalayan blocks, existed in the Tibetan region in the Early Cretaceous. The Himalayan Marginal and South Lhasa Seas formed in the southern and northern margins of the Neo- Tethys Ocean, the Central Tibet Sea and the Qiangtang Marginal Sea formed in the southern and northern margins of the Mid-Tethys Ocean, respectively. An arm of the sea extended into the southwestern Tarim basin in the Early Cretaceous. Early Cretaceous intensive thrusting, magmatic emplacement and volcanic eruptions occurred in the central and northern Lhasa Block, while strike- slip formed along the Hoh-Xil and South Kunlun Faults in the northern Tibetan region. Early Cretaceous tectonics together with magmatic K20 geochemistry indicate an Early Cretaceous southward subduction of the Mid-Tethys Oceanic Plate along the Bangoin-Nujiang Suture which was thrust ~87 km southward during the Late Cretaceous-Early Cenozoic. No intensive thrust and magmatic emplacement occurred in the Early Cretaceous in the Himalayan and southern Lhasa Blocks, indicating that the spreading Neo-Tethys Oceanic Plate had not been subducted in the Early Cretaceous. To the north, terrestrial basins of red-beds formed in the Hoh-Xil, Kunlun, Qilian and the northeastern Tarim blocks in Early Cretaceous, and the Qiangtang Marginal Sea disappeared after the Qiangtang Block uplifted in the late Early Cretaceous.