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?

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.

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.

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.

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.

An Orthogonal Collision Dynamic Mechanism of Wave-Like Uplift Plateaus in Southern Asia

In southern Asia, there are three large-scale wave-like mountains ranging from the Tibetan Plateau westward to the Iranian Plateau and the Armenian Plateau. On the southern side between plateaus, there are the Indian Peninsula and the Arabian Peninsula. What dynamic mechanisms form the directional alignment of the three plateaus with the two peninsulas remains a mystery. In the early stages of the Earth’s geological evolution, the internal structure of the Earth was that the center was a solid core, and the outmost layer was a thin equatorial crust zone separated by two thick pristine continents in polar areas, while the middle part was a deep magma fluid layer. Within the magma fluid layer, thermal and dynamic differences triggered planetary-scale vertical magma cells and led to the core-magma angular momentum exchange. When the core loses angular momentum and the magma layer gains angular momentum, the movement of upper magma fluids to the east and the tropical convergence zone (TCZ) drives the split and drift of two thick pristine continents, eventually forming the current combination of these plateaus and peninsulas and their wave-like arrangement along the east-west direction. Among them, the horizontal orthogonal convergence (collision) of upper magma fluids from the two hemispheres excited the vertical shear stress along the magma TCZ, which is the dynamic mechanism of mountain uplifts on the north side and plate subductions on the south side. To confirm this mechanism, two examples of low-level winds are used to calculate the correspondence between cyclone/anticyclonic systems generated by the orthogonal collision of airflows along the atmospheric TCZ and satellite-observed cloud systems. Such comparison can help us revisit the geological history of continental drift and orogeny.

洋中脊动力学与俯冲带地震-岩浆-成矿事件远程效应

板块俯冲带和大陆岩浆弧的深部过程与极端地质事件之间存在密切的关系。板块俯冲与造山过程会导致地震、岩浆活动和成矿等事件的发生。这些极端事件的发生与俯冲过程中的壳幔相互作用、地幔楔形成、岩石圈部分熔融、构造-岩浆活动等因素密切相关。然而,人们对洋中脊新生地壳的不均匀性或先天“缺陷”对以上的极端事件的长远影响和远程效应了解甚少。在洋中脊新的地壳形成过程中,由于受到板块扩张、压力降低、软流圈物质上涌等因素的作用,导致新生地壳温度升高、孔隙度和裂缝发育、密度降低、结构复杂的正反馈过程。因此,新生地壳在密度、强度、温度、厚度等方面存在非均质性。这些地壳的差异性将影响和决定板块在扩张和俯冲过程中的行为,并对板块俯冲作用形成的地震、岩浆和成矿等事件产生远程影响。以太平洋俯冲和安第斯造山带为例研究发现,板块运动速度、板块俯冲角度、板片撕裂、岩石圈厚度、Moho面深度等的突变与地震、火山和斑岩矿床的时空分布存在远程关联效应,这些认知对预测板块俯冲-碰撞带发生的极端地质事件的时空分布具有重要意义。

上覆板块对大陆深俯冲构造响应:以胶辽地区为例

在大陆深俯冲形成的超高压造山带中,由于超高压岩石折返构造的强烈改造,板块汇聚过程中的收缩构造难以在造山带内部识别,进而影响了对造山过程的完整认识。上覆板块,由于其构造部位的特殊性,较少受折返构造的影响,较大可能地记录了板块汇聚过程的收缩构造,成为理解板块汇聚过程和演化的关键地区。本文以华南板块大陆深俯冲过程中苏鲁超高压造山带上覆华北大陆岩石圈的构造响应为切入点,以胶辽地区为研究靶区,运用多尺度构造解析的方法,综合前人研究成果开展三叠纪收缩构造的流变学特点研究,结合造山带内部变形特征,探讨苏鲁超高压造山带及其上覆板块构造演化过程和动力学机制。研究结果表明,上覆板块在华南-华北汇聚过程中经历了早期上部向NE剪切的逆冲构造和晚期上部向SW剪切的反冲构造。早期收缩构造表现为胶北地体莱西单元早-中三叠世上部指向NE的韧性变形,辽东地体中-晚三叠世上部指向NE的韧性变形和极性向NE的褶皱-逆冲构造;从南向北,变形层次渐浅,时间渐新。晚期收缩构造主要表现为胶北地体林寺山单元晚三叠世早期上部指向SW的韧性变形和辽东地体层次较浅的上部向SW剪切的反冲构造。综合苏鲁超高压造山带晚三叠世晚期造成超高压变质岩石折返的伸展构造,我们建立了华南-华北板块汇聚的三阶段演化模型。我们认为华南板块的低角度-平板俯冲可能是造成上覆华北板块强烈收缩构造发育的重要原因。此外,华南-华北NE-SW向汇聚,但俯冲板块造山后折返方向为NW-SE向,造就了现今苏鲁超高压造山带的构造格局。

华南地区地幔过渡带结构及其动力学意义

华南地区同时受特提斯构造域和太平洋构造域影响,是研究板块相互作用的最佳场所之一.为研究华南地区深部动力学过程,本文基于国家固定台网30°N以南的宽频带台站数据,利用接收函数方法开展了华南地区地幔过渡带结构研究,并结合已有地质、地球物理资料,讨论了华南地区深部构造单元划分及其浅部构造响应.研究结果表明:(1)扬子克拉通下方地幔过渡带结构接近全球平均水平,相对稳定;(2)海南岛地区地幔过渡带厚度偏薄,且410-km和660-km界面分别显著下沉和上升,可能与海南地幔柱密切相关;(3)太平洋构造域和特提斯构造域深部构造边界显著,青藏高原东南缘主要受特提斯构造域影响,其地幔过渡带中可能残留有拆沉的岩石圈板片;(4)雪峰山以东地区的地幔过渡带主要受太平洋构造域影响,但南北也存在局部差异,南部华夏地块下方可能不存在滞留板片,太平洋板块的俯冲、后撤与东南沿海地区地壳减薄和大规模出露的中-新生代岩浆岩密切相关.