Signatures of downgoing plate-buoyancy driven subduction in Cenozoic plate motions

The dynamics of plate tectonics are strongly related to those of subduction. To obtain a better understanding of the driving forces of subduction, we compare relations between Cenozoic subduction motions at major trenches with the trends expected for the simplest form of subduction. i.e., free subduction, driven solely by the buoyancy of the downgoing plate. In models with an Earth-like plate stiffness (corresponding to a plate–mantle viscosity contrast of 2–3 orders of magnitude), free plates subduct by a combination of downgoing plate motion and trench retreat, while the slab is draped and folded on top of the upper-lower mantle viscosity transition. In these models, the slabs sink according to their Stokes’ velocities. Observed downgoing-plate motion–plate-age trends are compatible with >80% of the Cenozoic slabs sinking according to their upper-mantle Stokes’ velocity, i.e., subducting-plate motion is largely driven by upper-mantle slab pull. Only in a few cases, do young plates move at velocities that require a higher driving force (possibly supplied by lower-mantle–slab induced flow). About 80% of the Cenozoic trenches retreat, with retreat accounting for about 10% of the total convergence. The few advancing trench sections are likely affected by regional factors. The low trench motions are likely encouraged by low asthenospheric drag (equivalent to that for effective asthenospheric viscosity 2–3 orders below the upper-mantle average), and low lithospheric strength (effective bending viscosity ～2 orders of magnitude above the upper-mantle average). Total Cenozoic trench motions are often very oblique to the direction of downgoing-plate motion (mean angle of 73°). This indicates that other forces than slab buoyancy exert the main control on upper-plate/trench motion. However, the component of trench retreat in the direction of downgoing plate motion (≈ slab pull) correlates with downgoing-plate motion, and this component of retreat generally does not exceed the amount expected for free buoyancy-driven subduction. High present-day slab dips (on average about 70°) are compatible with largely upper-mantle slab-pull driven subduction of relatively weak plates, where motion partitioning and slab geometry adjust to external constraints/forces on trench motion.     

上地幔俯冲板块的动力学过程:数值模拟

大洋板块俯冲到地幔转换带,进而可形成不同的形态:板块可以停滞在660km不连续面,抑或穿过地幔转换带进入下地幔.这些不同的俯冲模式可进一步影响到海沟的运动.为更好地理解上地幔中俯冲板片的变形行为以及俯冲过程与海沟运动之间的关系,本文通过建立一系列高精度二维热-力学自由俯冲的数值模型,揭示了俯冲板块在上地幔中的变形方式及其与地幔转换带之间的相互作用过程.模拟结果显示,在俯冲板块与地幔转换带的相互作用过程中,其动力学过程可以分为以海沟后撤主导、海沟前进主导以及稳定型海沟等三种主要动力学类型.对于年龄较老,厚度较大的俯冲板块容易形成海沟后撤型俯冲,俯冲板块停滞在660km不连续面.相反,年龄较小,塑性强度较小的板块容易形成海沟前进型俯冲,俯冲板块穿越660km不连续面.     

The influence of trench migration on slab penetration into the lower mantle

A two-dimensional numerical convection model in cartesian geometry is used to study the influence of trench migration on the ability of subducted slabs to penetrate an endothermic phase boundary at 660 km depth. The transient subduction history of an oceanic plate is modelled by imposing plate and trench motion at the surface. The viscosity depends on temperature and depth. A variety of styles of slab behaviour is found, depending predominantly on the trench velocity. When trench retreat is faster than 2–4 cm/a, the descending slab flattens above the phase boundary. At slower rates it penetrates straight into the lower mantle, although flattening in the transition zone may occur later, leading to a complex slab morphology. The slab can buckle, independent of whether it penetrates or not, especially when there is a localised increase in viscosity at the phase boundary. Flattened slabs are only temporarily arrested in the transition zone and sink ultimately into the lower mantle. The results offer a framework for understanding the variety in slab geometry revealed by seismic tomography.     

Deep geodynamics of far field interconti- nental back-arc extension: Formation of Cenozoic volcanoes in northeastern China

Introduction Northeastem China has the most strong Cenozoic volcanism in China (Liu, 1999), where dis-tributes more than 500 Cenozoic volcanoes, including sleeping volcanoes of Tianchi Lake (Celes-tial Pond) of Changbai Mountain, and Wudalianchi (Five linked Lakes) (LIU, 1999). Vo lcano ofTianchi Lake of Changbai Mountain consists of basaltic rocks of shield-forming stage andtrachytes and pantellerites in cone-forming stage. It is suggested by study of REE, incompatibleelements a…     

Assembly and dispersal of supercontinents: The view from Avalon

Recent models for the post-750 Ma Rodinian supercontinent dispersal (e.g. Hoffman, 1991) envision that cratons margined by Grenvillian belts, were reorganized before ca 540 Ma to form the Gondwanan supercontinent. Laurentia and Baltica distanced themselves from Gondwana by moving out of the Rodinian cratonal cluster. West Gondwana, of which Avalon was a part during the late Proterozoic to Cambrian cratonal assembly, consisted mainly of Africa and South America.The main geological evidence is presented for: (1) a transition from continental platform conditions to those of a subduction-related volcanic arc regime in Late Proterozoic time during the dispersal of the Rodinian supercontinent, and the resulting assembly of the Gondwanan supercontinent; and (2) a second transition that marked a reversal from the volcanic arc regime to marine platformal environments by early Cambrian time.Evidence for progressive instability of the continental shelf margining the Rodinian supercontinent is contained in late Proterozoic olistostromes, mylonite zones, calc-alkaline magmatism, and arc-derived clastic rocks, some being glacigenic, during three phases of the Avalonian orogeny.By early Cambrian time the reversal from a tectonically unstable volcanic arc regime to more stable platformal conditions took place as Avalon, Armorica and related microcontinental blocks rifted from Gondwana. These Gondwanan fragments sequentially come into collision, first with each other and Baltica, and then with Laurentia in Mid to Late Paleozoic time as Pangaea was being assembled.     

Numerical models of subduction of the oceanic crust with basaltic plateaus

Light continents and islands characterized by a crustal thickness of more than 30 km float over a convective mantle, while the thin basaltic oceanic crust sinks completely in subduction zones. The normal oceanic crust is 7 km thick. However, anomalously thick basaltic plateaus forming as a result of emplacement of mantle plumes into moving oceanic lithospheric plates are also pulled into the mantle. One of the largest basaltic plateaus is the Ontong Java plateau on the Pacific plate, which arose during the intrusion of a giant superplume into the plate ～100 Myr ago. Notwithstanding its large thickness (averaging ～30 km), the Ontong Java plateau is still experiencing slow subduction. On the basis of numerical modeling, the paper analyzes the oceanic crust subduction process as a function of the mantle convection vigorousness and the density, thickness, viscosity, and shape of the crust. Even a simplified model of thermocompositional convection in the upper mantle is capable of explaining the observed facts indicating that the oceanic crust and sediments are pulled into the mantle and the continental crust is floating on the mantle.     

Deep geodynamics of far field intercontinental back-arc extension: Formation of Cenozoic volcanoes in northeastern China

There are three cases of variation of trench location possible to occur during subduction: trench fixed, trench advancing, and trench retreating. Retreat of trench may lead to back-arc extension. The Pacific plate subducts at low angle beneath the Eurasia plate, tomographic results indicate that the subducted Pacific slab does not penetrate the 670 km discontinuity, instead, it is lying flat above the interface. The flattening occurred about 28 Ma ago. Geodynamic computation suggests: when the frontier of the subducted slab reaches the phase boundary of lower and upper mantle, it may be hindered and turn flat lying above the boundary, facilitates the retreat of trench and back-arc extension. Volcanism in northeastern China is likely a product of such retreat of subduction, far field back-arc extension, and melting due to reduce of pressure while mantle upwelling. Foundation item: National Natural Science Foundation of China (40234042 and 40174027).     

Cyclical one-way continental rupture-drift in the Tethyan evolution: Subduction-driven plate tectonics

Numerous continents have rifted and drifted away from Gondwana to repeatedly open ocean basins over the past-500 millionyears.These Gondwana-derived continents drifted towards and collided with components of the Eurasian continent to successively close the preexisting oceans between the two.Plate tectonics satisfactorily describes the continental drift from Gondwana to Eurasia but does not define the geodynamic mechanism of continuously rifting to collisions of continents in the Tethy an Realm.After reappraisal of geological records of the rift,collision and subduction initiation from the surface and various geophysical observations from depth,we propose that Eurasia-directed subducting oceanic slabs would have driven Tethyan system in the Phanerozoic.The Eurasia-directed subduction would have dragged the passive Gondwana margin to rift and drift northwards,giving birth to new oceans since the Paleozoic.The closure of preexisting oceans between the Gondwana-derived continents and Eurasia led to continental collisions,which would have induced the initiation of oceanic subduction in the Tethyan Realm.Multiple episodic switches between collision-subduction-rift repeatedly led to the separation of continental fragments from Gondwana and dragged them to drift towards Eurasia.The final disappearance of Neo-Tethy s would have induced collision of the Gondwana-derived continents with the Eurasian continent,giving rise to the Cenozoic Alpine-Zagros-Himalayan collisional system.Therefore,the Eurasia-directed oceanic subduction would have acted as a 'one-way train' that successively transferred the ruptured Gondwana continental fragments in the south,into the terminal in the north.In this regard,the engine of this "Tethyan one-way train" is the negative buoyancy of subducting oceanic slabs.     

西太平洋板块向我国东北地区深部俯冲的数值模拟

本文采用依赖温度的黏度结构以及考虑海洋板块和大陆板块厚度差异等特征,以太平洋板块向欧亚板块会聚速率作为板块速度的主要约束,通过变化海沟后撤速度模型,数值模拟西太平洋板块向中国东北的俯冲过程.结果表明,要产生类似于中国东北之下低角度的板片俯冲,海沟后撤是重要条件;而上下地幔黏度的较大差异是决定俯冲板片不穿透660 km相变面的决定因素;西太平洋板块向欧亚板块的俯冲应早于70 Ma B.P.,海沟后撤速度可能小于一些地质学家估计的45 mm/a, 而且可能是分阶段变化的;速度场表明运动学模型的反过程：大陆岩石圈之下物质的不断水平向东的流动和推挤可能成为海沟后撤的力源之一,地幔物质的这种东向流动可能与印度板块挤压碰撞欧亚板块有关,沿欧亚板块东缘的扩张构造可能是太平洋-欧亚板块运动和印度-欧亚板块运动的综合效应.     

From incipient island arc to doubly‐vergent orogen: A review of geodynamic models and sedimentary basin‐fills of southern Central America

Southern Central America is a Late Mesozoic/Cenozoic island arc that evolved in response to the subduction of the Farallón Plate beneath the Caribbean Plate in the Late Cretaceous and, from the Oligocene, the Cocos and Nazca Plates. Southern Central America is one of the best studied convergent margins in the world. The aim of this paper is to review the sedimentary and structural evolution of arc‐related sedimentary basins in southern Central America, and to show how the arc developed from a pre‐extensional intra‐oceanic island arc into a doubly‐vergent, subduction orogen. The Cenozoic sedimentary history of southern Central America is placed into the plate tectonic context of existing Caribbean Plate models. From regional basin analysis, the evolution of the southern Central American island arc is subdivided into three phases: (i) non‐extensional stage during the Campanian; (ii) extensional phase during the Maastrichtian‐Oligocene with rapid basin subsidence and deposition of arc‐related, clastic sediments; and (iii) doubly‐vergent, compressional arc phase along the 280 km long southern Costa Rican arc segment related to either oblique subduction of the Nazca plate, west‐to‐east passage of the Nazca–Cocos–Caribbean triple junction, or the subduction of rough oceanic crust of the Cocos Plate. The Pleistocene subduction of the Cocos Ridge contributed to the contraction but was not the primary driver. The architecture of the arc‐related sedimentary basin‐fills has been controlled by four factors: (i) subsidence caused by tectonic mechanisms, linked to the angle and morphology of the incoming plate, as shown by the fact that subduction of aseismic ridges and slab segments with rough crust were important drivers for subduction erosion, controlling the shape of forearc and trench‐slope basins, the lifespan of sedimentary basins, and the subsidence and uplift patterns; (ii) subsidence caused by slab rollback and resulting trench retreat; (iii) eustatic sea‐level changes; and (iv) sediment dispersal systems.     

SEASAT geoid anomalies and the Macquarie Ridge complex

The seismically active Macquarie Ridge complex forms the Pacific-India plate boundary between New Zealand and the Pacific-Antarctic spreading center. The Late Cenozoic deformation of New Zealand and focal mechanisms of recent large earthquakes in the Macquarie Ridge complex appear consistent with the current plate tectonic models. These models predict a combination of strike-slip and convergent motion in the northern Macquarie Ridge, and strike-slip motion in the southern part. The Hjort trench is the southernmost expression of the Macquarie Ridge complex. Regional considerations of the magnetic lineations imply that some oceanic crust may have been consumed at the Hjort trench. Although this arcuate trench seems inconsistent with the predicted strike-slip setting, a deep trough also occurs in the Romanche fracture zone.Geoid anomalies observed over spreading ridges, subduction zones, and fracture zones are different. Therefore, geoid anomalies may be diagnostic of plate boundary type. We use SEASAT data to examine the Macquarie Ridge complex and find that the geoid anomalies for the northern Hjort trench region are different from the geoid anomalies for the Romanche trough. The Hjort trench region is characterized by an oblique subduction zone geoid anomaly, e.g., the Aleutian-Komandorski region. Also, limited first-motion data for the large 1924 earthquake that occurred in the northern Hjort trench suggest a thrust focal mechanism. We conclude that subduction is occurring at the Hjort trench. The existence of active subduction in this area implies that young oceanic lithosphere can subduct beneath older oceanic lithosphere.     

The Role of History-Dependent Rheology in Plate Boundary Lubrication for Generating One-Sided Subduction

We have developed a two-dimensional dynamical model of asymmetric subduction integrated into the mantle convection without imposed plate velocities. In this model we consider that weak oceanic crust behaves as a lubricator on the thrust fault at the plate boundary. We introduce a rheological layer that depends on the history of the past fracture to simulate the effect of the oceanic crust. The thickness of this layer is set to be as thin as the Earth's oceanic crust. To treat 1-kilometer scale structure at the plate boundary in the 1000-kilometer scale mantle convection calculation, we introduce a new numerical method to solve the hydrodynamic equations using a couple of uniform and nonuniform grids of control volumes. Using our developed models, we have systematically investigated effects of basic rheological parameters that determine the deformation strength of the lithosphere and the oceanic crust on the development of the subducted slab, with a focus on the plate motion controlling mechanism. In our model the plate subduction is produced when the friction coefficient (0.004–0.008) of the modeled oceanic crust and the maximum strength (400 MPa) of the lithosphere are in plausible range inferred from the observations on the plate driving forces and the plate deformation, and the rheology experiments. In this range of the plate strength, yielding induces the plate bending. In this case the speed of plate motion is controlled more by viscosity layering of the underlying mantle than by the plate strength. To examine the setting of the overriding plate, we also consider the two end-member cases in which the overriding plate is fixed or freely-movable. In the case of the freely-movable overriding plate, the trench motion considerably changes the dip angle of the deep slab. Especially in the case with a shallow-angle plate boundary, retrograde slab motion occurs to generate a shallow-angle deep slab.     

Global scale patterns of continental fragmentation: Wilson's cycles as a constraint for long-term sea-level changes

We propose to characterize land–ocean distributions over Late Proterozoic to Phanerozoic times from measurement of perimeters and areas of continental fragments, based on paleomagnetic reconstructions. These measurements serve to calculate geophysically constrained breakup and scatter indexes of continental land masses from 0 to 1100 Ma. We then provide quantitative investigation and modelling of relationships between scatter of continental landmasses and mean age of the oceanic lithosphere during Mesozoic times, which appears to range from 56 to 62 Ma over the last 170 My. We then inverse the scatter of continental landmasses in terms of global oceanic crust mean age over the last 600 My, i.e. back in times where no measurement of seafloor accretion history is possible because of subduction. We finally show that the inferred evolution of oceanic lithosphere mean age over the Phanerozoic remarkably correlates in time with long-term sea-level changes since the Cambrian.     

Mantle flow with existence of plates and generation of the toroidal field

The observed plate velocities contain two types of motions. The poloidal component is related to the formation of ridges and subduction zones and the toroidal field expresses the shearing of surface plates. One very important consideration in modeling flow in the earth's mantle is the existence and motion of the lithospheric plates. The motion of plates represents a large-scale circulation with strong viscous coupling to the mantle underneath. The mantle flow probably is neither a purely free convection driven by buoyancy forces due to nonadiabatic temperature gradients in the mantle nor a forced convection generated by boundary forces, but a mixed convection that combines the effects of boundary and buoyancy forces. We present, in this paper, the mixed convection model resulting in a surface velocity field that contains both the observed poloidal and toroidal components.     

Effects of relative plate motion on the deep structure and penetration depth of slabs below the Izu-Bonin and Mariana island arcs

An increasing number of seismological studies indicate that slabs of subducted lithosphere penetrate the Earth's lower mantle below some island arcs but are deflected, or, rather, laid down, in the transition zone below others. Recent numerical simulations of mantle flow also advocate a hybrid form of mantle convection, with intermittent layering. We present a multi-disciplinary analysis of slab morphology and mantle dynamics in which we account explicitly for the history of subduction below specific island arcs in an attempt to understand what controls lateral variations in slab morphology and penetration depth. Central in our discussion are the Izu-Bonin and Mariana subduction zones. We argue that the differences in the tectonic evolution of these subduction zones—in particular the amount and rate of trench migration—can explain why the slab of subducted oceanic lithosphere seems to be (at least temporarily) stagnant in the Earth's transition zone below the Izu-Bonin arc but penetrates into the lower mantle below the Mariana arc. We briefly speculate on the applicability of our model of the temporal and spatial evolution of slab morphology to other subduction zones. Although further investigation is necessary, our tentative model shows the potential for interpreting seismic images of slab structure by accounting for the plate-tectonic history of the subduction zones involved. We therefore hope that the ideas outlined here will stimulate and direct new research initiatives.     

Lagging mantle convection,the geoid and mantle structure

In numerical models of convection incorporating migration of a simulated subduction zone, the main descending flow lags far behind the migrating trench, and a geoid low is associated with the main descending flow. This provides physical plausibility for the suggestion by Chase and Sprowl that present very long-wavelength (degree 2–4) geoid lows are associated with Mesozoic trench locations, and suggests further that the present long-wavelength geoid, deep mantle structure and hotspot distribution may be straightforward consequences of plate evolution since the Upper Paleozoic.     

The generation of plate tectonics from mantle convection

In the last decade, significant progress has been made toward understanding how plate tectonics is generated from mantle dynamics. A primary goal of plate-generation studies has been the development of models that allow the top cold thermal boundary layer of mantle convection, i.e. the lithosphere, to develop broad and strong plate-like segments separated by narrow, weak and rapidly deforming boundaries; ideally, such models also permit significant strike-slip (toroidal) motion, passive ridges (i.e. pulled rather than pried apart), and self-consistent initiation of subduction. A major outcome of work so far is that nearly all aspects of plate generation require lithospheric rheologies and shear-localizing feedback mechanisms that are considerably more exotic than rheologies typically used in simple fluid-dynamical models of mantle flow. The search for plate-generating behavior has taken us through investigations of the effects of shear weakening (‘stick-slip’) and viscoplastic rheologies, of melting at ridges and low-viscosity asthenospheres, and of grain-size dependent rheologies and damage mechanics. Many such mechanisms, either by themselves or in combination, have led to self-consistent fluid-mechanical models of mantle flow that are remarkably plate-like, which is in itself a major accomplishment. However, many other important problems remain unsolved, such as subduction intiation and asymmetry, temporal evolution of plate geometry, rapid changes in plate motion, and the Archaean initiation of the plate-tectonic mode of convection. This paper presents a brief review of progress made in the plate-generation problem over the last decade, and discusses unresolved issues and future directions of research in this important area.     

Linking the Alps and Apennines subduction systems: New constraints revealed by high-resolution teleseismic tomography

We report a new model of the upper mantle structure beneath Italy obtained by means of P-wave teleseismic tomography. Besides the recent and remarkable development of the Italian Seismic Network, a high model resolution has been achieved improving the inversion method upon the ACH method used in previous investigations and picking high quality arrival times with the Multi-Channel Cross-Correlation technique. The finer details of our Vp model yield new insights into the heterogeneous structure of the Adria continental lithosphere involved in the collision between the Africa and Europe plates. A wide low Vp anomaly located in the northern Adria mantle, facing the Alpine high Vp slab, supports the idea that the Adria lithosphere has been hydrated and thinned during the Alpine subduction. We argue that this mantle softening may have played a key role in favoring the subsequent delamination of the Adria lithosphere in the northern Apennines. We hypothesize that delamination of continental lithosphere previously thinned in a back-arc setting may be considered a key process to favor subduction polarity reversal and recycling of continental material into the mantle circulation. Conversely, in the central-southern Apennines, the velocity structure is consistent with the existence of a deeper oceanic slab that flattens at the base of the upper mantle, in agreement with the widely accepted geodynamic evolution of the central Mediterranean by slab retreat and back-arc spreading. The oceanic slab is discontinuously detached from the surface plate, suggesting a different structure of the Adria lithosphere, which resists subduction instead of favoring delamination.     

Mantle heterogeneity beneath the southern Mid-Atlantic Ridge: trace element evidence for contamination of ambient asthenospheric mantle

We report new trace element data for an extensive suite of quench basalt glasses dredged from the southern Mid-Atlantic Ridge (MAR) between 40°S and 52.5°S. Ratios between highly incompatible trace elements are strongly correlated and indicate a systematic distribution of incompatible element enriched mid-ocean ridge basalt (MORB) (E-type: Zr/Nb=5.9-19, Y/Nb=0.9-8.4, (La/Sm)_n=1.0-2.9) and incompatible element depleted MORB (N-type: Zr/Nb=30-69, Y/Nb=11-29, (La/Sm)_n=0.48-0.79) along this section of the southern MAR. A notable feature of N-type MORB from the region is the higher than usual Ba/Nb (4-9), La/Nb (1.2-2.4) and primitive mantle normalised K/Nb ratios (>1). Ba/Nb ratios in E-type MORB samples from 47.5 to 49°S are especially elevated (>10). The occurrence and geographic distribution of E-type MORB along this section of the southern MAR can be correlated with the ridge-centred Shona and off-axis Discovery mantle plumes. In conjunction with published isotope data for a subset of the same sample suite [Douglass et al., J. Geophys. Res. 104 (1999) 2941], a model is developed whereby prior to the breakup of Gondwana and the opening of the South Atlantic Ocean, the underlying asthenospheric mantle was locally contaminated by fluids/melts rising from the major Mesozoic subduction zone along the south-southwest boundary of Gondwana, leaving a subduction zone geochemical imprint (elevated (K/Nb)_n and ⁸⁷Sr/⁸⁶Sr ratios, decreased ¹⁴³Nd/¹⁴⁴Nd ratios). Subsequent impingement of three major mantle plume heads (Tristan/Gough, Discovery, Shona) resulted in heating and thermal erosion of the lowermost subcontinental lithosphere and dispersal into the convecting asthenospheric mantle. With the opening of the ocean basin, continued plume upwelling led to plume-ridge interactions and mixing between geochemically enriched mantle derived from the Shona and Discovery mantle plumes, material derived from delamination of the subcontinental lithosphere, and mildly subduction zone contaminated depleted asthenospheric mantle.     

Subduction tectonics vs.plume tectonics——Discussion on driving forces for plate motion

Plate tectonics describes the horizontal motions of lithospheric plates,the Earths outer shell,and interactions among them across the Earths surface.Since the establishment of the theory of plate tectonics about half a century ago,considerable debates have remained regarding the driving forces for plate motion.The early"Bottom up"view,i.e.,the convecting mantledriven mechanism,states that mantle plumes originating from the core-mantle boundary act at the base of plates,accelerating continental breakup and driving plate motion.Toward the present,however,the"Top down"idea is more widely accepted,according to which the negative buoyancy of oceanic plates is the dominant driving force for plate motion,and the subducting slabs control surface tectonics and mantle convection.In this regard,plate tectonics is also known as subduction tectonics."Top down"tectonics has received wide supports from numerous geological and geophysical observations.On the other hand,recent studies indicate that the acceleration/deceleration of individual plates over the million-year timescale may reflect the effects of mantle plumes.It is also suggested that surface uplift and subsidence within stable cratonic areas are correlated with plumerelated magmatic activities over the hundred-million-year timescale.On the global scale,the cyclical supercontinent assembly and breakup seem to be coupled with superplume activities during the past two billion years.These correlations over various spatial and temporal scales indicate the close relationship and intensive interactions between plate tectonics and plume tectonics throughout the history of the Earth and the considerable influence of plumes on plate motion.Indeed,we can acquire a comprehensive understanding of the driving forces for plate motion and operation mechanism of the Earth's dynamic system only through joint analyses and integrated studies on plate tectonics and plume tectonics.