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.     

Subduction-zone peridotites and their records of crust-mantle interaction

Subduction is the core process of plate tectonics. The mantle wedge in subduction-zone systems represents a key tectonic unit, playing a significant role in material cycling and energy exchange between Earth's layers. This study summarizes research progresses in terms of subduction-related peridotite massifs, including supra-subduction zone(SSZ) ophiolites and mantle-wedge-type(MWT) orogenic peridotites. We also provide the relevant key scientific questions that need be solved in the future. The mantle sections of SSZ ophiolites and MWT orogenic peridotites represent the mantle fragments from oceanic and continental lithosphere in subduction zones, respectively. They are essential targets to study the crust-mantle interaction in subduction zones. The nature of this interaction is the complex chemical exchanges between the subducting slab and the mantle wedge under the major control of physical processes. The SSZ ophiolites can record melt/fluid-rock interaction, metamorphism,deformation, concentration of metallogenic elements and material exchange between crust and mantle, during the stages from the generation of oceanic lithosphere at spreading centers to the initiation, development, maturation and ending of oceanic subduction at continental margins. The MWT orogenic peridotites reveal the history of strong metamorphism and deformation during subduction, the multiple melt/fluid metasomatism(including silicatic melts, carbonatitic melts and silicate-bearing C-HO fluids/supercritical fluids), and the complex cycling of crust-mantle materials, during the subduction/collision and exhumation of continental plates. In order to further reveal the crust-mantle interaction using subduction-zone peridotites, it is necessary to utilize high-spatial-resolution and high-precision techniques to constrain the complex chemical metasomatism, metamorphism,deformation at micro scales, and to reveal their connections with spatial-temporal evolution in macro-scale tectonics.     

Progress in numerical modeling of subducting plate dynamics

The core concerns of plate tectonics theory are the dynamics of subducting plates, which can be studied by integrating multidisciplinary fields such as seismology, mineral physics, rock geochemistry, geological formation studies, sedimentology, and numerical simulations. By establishing a theoretical model and solving it with numerical methods, one can replicate the dynamic effects of a subducting plate, quantifying its evolution and the surface response. Simulations can also explain the observations and experimental results of other disciplines. Therefore, numerical models are among the most important tools for studying the dynamics of subducting plates. This paper provides a review on recent advances in the numerical modeling of subducting plate dynamics. It covers various aspects, namely, the origin of plate tectonics, the initiation process and thermal structure of subducting slab, and the main subduction slab dynamics in the upper mantle, mantle transition zone, and lower mantle. The results of numerical models are based on the theoretical equations of mass, momentum, and energy conservation. To better understand the dynamic progress of subducting plates, the simulation results must be verified in comparisons with the results from natural observations by geology, geophysics and geochemistry. With the substantial increase in computing power and continuous improvement of simulation methods, numerical models will become a more accurate and efficient means of studying the frontier issues of Earth sciences, including subducting plate dynamics.     

南海瑞雷面波群速度层析成像及其地球动力学意义

南海处于欧亚板块、太平洋板块和印度—澳大利亚板块的交会区,是西北太平洋一系列边缘海中最大的边缘海。关于南海的打开以往研究提出了如印度板块与欧亚板块碰撞驱动挤出以及古南海俯冲拖拽等诸多模型。本文力图通过南海海盆及周边各向异性结构来约束南海演化机制。基于同济大学2012和2014年在南海中央海盆进行的两次被动源宽频带海底地震观测试验回收的10台OBS记录仪近1年的地震数据,本文采用三种不同的横波分裂方法,获取了中央海盆针对两次远震的XKS分裂结果以及南海周边20次区域地震提供的S震相分裂结果。SKS分裂结果显示,南海中央海盆下方存在快轴方向为NE-SW向的各向异性,其成因可能与海底扩张时期沿洋脊方向的地幔流以及南海海洋板块俯冲拖拽的地幔流有关。南海及其周边上地幔存在强各向异性,且不同方位观测到的各向异性不同,快轴方向与前人SKS横波分裂结果、GPS和板块运动一致,较好地对应了区域构造运动或者地幔对流模型。各向异性结果与印度—欧亚板块碰撞驱动挤出模型以及古南海俯冲板块拖拽模型预期结果一致,与理想的地幔柱上涌驱动模型不一致。由于海盆各向异性观测特别有限,各向异性结果不能证实亦不能证伪“大西洋型”海底扩张模型、弧后扩张模型和板缘破裂模型,后续还需要更多的观测结果来证实或证伪上述模型。     

Reconciling strong slab pull and weak plate bending: The plate motion constraint on the strength of mantle slabs

Although subducting slabs undergo a bending deformation that resists tectonic plate motions, the magnitude of this resistance is not known because of poor constraints on slab strength. However, because slab bending slows the relatively rapid motions of oceanic plates, observed plate motions constrain the importance of bending. We estimated the slab pull force and the bending resistance globally for 207 subduction zone transects using new measurements of the bending curvature determined from slab seismicity. Predicting plate motions using a global mantle flow model, we constrain the viscosity of the bending slab to be at most ~ 300 times more viscous than the upper mantle; stronger slabs are intolerably slowed by the bending deformation. Weaker slabs, however, cannot transmit a pull force sufficient to explain rapid trenchward plate motions unless slabs stretch faster than seismically observed rates of ~ 10^− 15 s^− 1. The constrained bending viscosity (~ 2 × 10²³ Pa s) is larger than previous estimates that yielded similar or larger bending resistance (here ~ 25% of forces). This apparent discrepancy occurs because slabs bend more gently than previously thought, with an average radius of curvature of 390 km that permits subduction of strong slabs. This gentle bending may ultimately permit plate tectonics on Earth.     

金星:认识早期地球的窗口

无论在行星大小、质量还是轨道速度等方面,金星都是太阳系中与地球最相似的行星.自1960年代初期开始,金星一直是人类深空探测的重要目标.本文简要地回顾了人类探索金星的历史,总结了对金星已有的认识,梳理了金星的主要科学问题,最后介绍了未来的国际探测计划,并建议了我国的金星探测目标.早期对金星的探测以苏联的金星计划(Венера)和美国的水手系列(Mariner)为代表,后期的探测器以欧盟、日本等国家的“金星快车(Venus Express)”、“拂晓号(Akatsuki)”为代表.这些探测结果为我们认识金星大气成分、地表地形和内部结构提供了重要的数据.金星的大气组成以CO₂为主,含少量N2,与现在地球的大气组成显著不同,类似早期地球的大气组成.虽然金星地表目前没有液态水,但部分理论模拟工作表明金星地表可能曾经有液态水.一系列探测器对金星地表成分的分析表明,金星地表主要由玄武岩组成.在地形地貌方面,由于金星特殊的地表环境,金星表面风化作用对地表地貌影响很小.金星的地表主要受控于比较年轻的火山作用,发育了许多不同于地球的地貌特征,主要包括区域平原、盾状火山平原、冕状地形以及瓦片状地形等,其动力学机制可能是地幔柱—岩石圈相互作用或地幔对流,至今未发现与板块构造相关的地貌.现阶段金星没有太多大型的、活跃的火山热点,虽然无法估测准确的火山活动速率,但相比地球来说火山活动速率小很多.在内部结构方面,金星具有与地球类似的核幔壳结构.金星的内部组成也与地球类似,例如金星地幔很可能是与地球相似的橄榄岩成分.不存在内部磁场和缺乏板块构造是金星区别于地球的两个重要特征.关于金星为什么没有自身磁场,主流观点是金星地核缺乏对流,无法演化出磁场.而针对金星为什么没有演化出板块构造,目前认为主要有三个可能的原因:地表温度过高,没有软流圈,金星缺乏液态水,其中液态水的缺乏接受度最广.从大气组成、地表岩石组合、构造作用等角度来看,金星都与早期地球非常相似,是我们理解类地行星演化的天然实验室.研究金星和地球为什么会朝不同方向演化,是深入理解包括系外行星在内的行星的宜居性形成与演变的重要途径。因此,金星一直是优先级别最高的深空探测目标之一.近几年,美国、俄罗斯以及欧洲等国家和地区分别针对金星目前主要的科学问题,例如金星是否存在早期海洋、金星的宜居性以及结构和重力场等,先后提出各自的金星探测计划.我国在新的国际竞争中应该、也必然有所作为.     

Assessment of the cooling capacity of plate tectonics and flood volcanism in the evolution of Earth, Mars and Venus

Geophysical arguments against plate tectonics in a hotter Earth, based on buoyancy considerations, require an alternative means of cooling the planet from its original hot state to the present situation. Such an alternative could be extensive flood volcanism in a more stagnant-lid like setting. Starting from the notion that all heat output of the Earth is through its surface, we have constructed two parametric models to evaluate the cooling characteristics of these two mechanisms: plate tectonics and basalt extrusion/flood volcanism. Our model results show that for a steadily (exponentially) cooling Earth, plate tectonics is capable of removing all the required heat at a rate of operation comparable to or even lower than its current rate of operation, contrary to earlier speculations. The extrusion mechanism may have been an important cooling agent in the early Earth, but requires global eruption rates two orders of magnitude greater than those of known Phanerozoic flood basalt provinces. This may not be a problem, since geological observations indicate that flood volcanism was both stronger and more ubiquitous in the early Earth. Because of its smaller size, Mars is capable of cooling conductively through its lithosphere at significant rates, and as a result may have cooled without an additional cooling mechanism. Venus, on the other hand, has required the operation of an additional cooling agent for probably every cooling phase of its possibly episodic history, with rates of activity comparable to those of the Earth.     

论青藏高原及邻区板片构造的一个新模式

本文首先论述了板块学说提出的过程和存在的一些不足与疑问,特别是该学说将Holmes（1948）的地幔热对流说作为驱使岩石圈板块运动的动力机制.而后又以青藏高原及邻区为例,根据区域地质、蛇绿岩和地质构造研究的成果,特别是地震测深研究的成果,详细地论证了本区不存在有大洋中脊扩张成为大洋盆地的新大洋和大洋板块简单的B型俯冲模式,但存在有海底扩张的陆间海和海洋地壳板片（蛇绿岩构造岩片）的仰冲以及大陆岩石圈板片复杂的A型俯冲新模式.新模式不是以地幔对流运动,而是以扩张分离A型俯冲的大陆岩石圈板片与软流圈之间的水平剪切相对运动机制作为它的躯动力.     

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.     

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

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

     

俯冲初始时板块分界面形状对俯冲过程的影响

自板块理论建立以来,俯冲一直是学者们关心的热点问题.前人结合地质、地球物理、实验室物理实验和数值模拟等多种手段对这一问题进行了大量的研究.以往的研究更为关注俯冲过程中板块的作用、地幔流动的规律和物质的迁移与相变等问题,却常忽视了俯冲是如何开始的这一基本问题.同时,由于相关数据资料较为有限,更限制了俯冲启动的相关研究.因此,本文选取俯冲启动问题中板块分界面形状对俯冲过程的影响这一问题,使用有限元的方法进行了数值模拟.我们选择针对倾斜型、垂直型和弯曲型三种不同形状的板块分界面建立对比模型,比较它们演化至10 Ma的过程我们发现：分界面几何形状的不同的确会对俯冲板块演化和海沟的深度产生影响.倾斜型模型的俯冲角度最大,海沟深度最深,俯冲深度最深;垂直型模型的俯冲角度和海沟深度仅次于倾斜型模型,俯冲深度最浅;弯曲型模型的俯冲角度最小,海沟深度最小,俯冲的深度介于倾斜型和垂直型之间.结合以上结论不难看出,俯冲角和海沟深度变化具有一定的相关性,俯冲角度越大,相应的海沟深度越大.

     

大尺度地幔动力学研究的现状和展望

这篇综述讨论大空间、大时间尺度的地幔动力学近几十年的发展和现状,着重讨论了相关的观测及其动力学意义.这些观测包括现在地球的板块运动的基本特性,中、长波重力异常及大地水准面异常,地震层析成像得到的地幔结构,以及过去10亿年超级大陆Pangea和Rodinia的形成、裂解和演化,及火山岩浆活动.关于地球动力学模型的讨论是围绕着这些相关的观测而进行的.涉及到的一些主要问题包括以下.第一,地幔动力学研究显示,地震层析成像得到的下地幔的二阶结构（比如核幔边界附近的LLSVP结构）,和俯冲带的快速异常体,可以解释为过去1亿年左右的板块运动和地幔对流的结果;第二,地幔三维结构作为地幔对流的驱动力,是导致中、长波重力及大地水准面异常的直接原因;结合地幔动力学模拟,观测的大地水准面异常对地幔黏性结构提供了强有力的约束,很可靠的结果之一是下地幔的黏性比上地幔要高至少一个量级,并且最近的研究确定软流圈的存在;第三,过去10亿年大陆块体经历过的Rodinia和Pangea两期超级大陆的形成和破裂是地幔动力学在地表的反映.地幔结构在Pangea形成过程中是一阶结构（即一个半球是冷的下降流,而另一个半球是热的上涌流）主导的,而现在的二阶为主导的地幔结构是Pangea形成后,破裂前或破裂过程中才形成的;地幔动力学和其他研究支持地幔结构在一阶和二阶间转换的1-2-1模型;第四,板块构造在地球上的起源和动力机制依然是充满争议和不确定的课题,但是这些问题同时也是重要的地球动力学基本问题.

     

Convection of a fluid with strongly temperature and pressure dependent viscosity

Plate tectonics on the Earth is a surface manifestation of convection within the Earth’s mantle, a subject which is as yet improperly understood, and it has motivated the study of various forms of buoyancy-driven thermal convection. The early success of the high Rayleigh number constant viscosity theory was later tempered by the absence of plate motion when the viscosity is more realistically strongly temperature dependent, and the process of subduction represents a continuing principal conundrum in the application of convection theory to the Earth. A similar problem appears to arise if the equally strong pressure dependence of viscosity is considered, since the classical isothermal core convection theory would then imply a strongly variable viscosity in the convective core, which is inconsistent with results from post-glacial rebound studies. In this paper we address the problem of determining the asymptotic structure of high Rayleigh number convection when the viscosity is strongly temperature and pressure dependent, i.e. thermobaroviscous. By a method akin to lid-stripping, we are able to extend numerical computations to extremely high viscosity contrasts, and we show that the convective cells take the form of narrow, vertically-oriented fingers. We are then able to determine the asymptotic structure of the solution, and it agrees well with the numerical results. Beneath a stagnant lid, there is a vigorous convection in the upper part of the cell, and a more sluggish, higher viscosity flow in the lower part of the cell. We then offer some comments on the possible meaning and interpretation of these results for planetary mantle convection.     

Regionalized models for the thermal evolution of the Earth

Traditional models for the heat loss in oceanic and continental regions are combined into a regionalized model for the thermal evolution of the Earth. The need for regionalization is obvious when one considers that the mantle loses 3 to 4 times as much heat per unit area in oceanic regions than in continental areas. The present-day rate of heat loss together with a geochemical estimate of the concentration of heat-producing elements in the Earth fixes the response time of the thermally convecting mantle. The response time in turn can be used to select the most reasonable representation for mantle convection in terms of the sensitivity of viscosity on temperature and layering versus mantle-wide circulation. Present geochemical estimates of the bulk composition of the Earth are most easily reconciled with the observed heat flow if the mantle is layered and its rheology is slightly less temperature dependent than generally assumed. The layered system can produce sufficiently high temperatures to explain the high-magnesian komatiites of the Archean. One difficulty with the models is that they predict widespread melting at shallow depth in the early stages of Earth history but do not address how such melting affects and alters the heat transfer mechanisms.     

汤加地区地幔结构及俯冲演化的地球动力学研究

汤加的俯冲板块结构复杂,板块形态的起因和形成机制不是很清晰.为了探索其演化,我们借助于板块运动学、海底年龄和关键的大地构造特征,用地球动力学模型模拟了本区域40 Ma之后的俯冲历史.通过研究我们发现,沿着海沟的分布位置,不同的黏度、流体的吸力与浮力比特征控制着板块的力学强度,从而影响着板块的俯冲角度和形态.我们的模型显示在俯冲板块形成前,研究区表面西侧的黏度剖面形成的弱化带为水平低黏度区,它对俯冲的发生有一定的引导作用,高黏度的上覆板块的存在和俯冲导致的地幔流体的吸力降低了俯冲板块的倾角,15°S附近的俯冲板块的分离发生在16-8 Ma之间.最后我们通过匹配实际的构造特征,包括贝尼奥夫带和波速结构剖面,用所得的相关性最佳的模型对本区的地幔结构随时间演化提供探索.

     

Numerical model for the generation of the ensemble of lithospheric plates and their penetration through the 660-km boundary

In the kinematic theory of lithospheric plate tectonics, the position and parameters of the plates are predetermined in the initial and boundary conditions. However, in the self-consistent dynamical theory, the properties of the oceanic plates (just as the structure of the mantle convection) should automatically result from the solution of differential equations for energy, mass, and momentum transfer in viscous fluid. Here, the viscosity of the mantle material as a function of temperature, pressure, shear stress, and chemical composition should be taken from the data of laboratory experiments. The aim of this study is to reproduce the generation of the ensemble of the lithospheric plates and to trace their behavior inside the mantle by numerically solving the convection equations with minimum a priori data. The models demonstrate how the rigid lithosphere can break up into the separate plates that dive into the mantle, how the sizes and the number of the plates change during the evolution of the convection, and how the ridges and subduction zones may migrate in this case. The models also demonstrate how the plates may bend and break up when passing the depth boundary of 660 km and how the plates and plumes may affect the structure of the convection. In contrast to the models of convection without lithospheric plates or regional models, the structure of the mantle flows is for the first time calculated in the entire mantle with quite a few plates. This model shows that the mantle material is transported to the mid-oceanic ridges by asthenospheric flows induced by the subducting plates rather than by the main vertical ascending flows rising from the lower mantle.     

Evolution of a section of the Africa-Europe plate boundary: Paleomagnetic and volcanological evidence from Sicily

New paleomagnetic data relative to Upper Cretaceous, Neogene and Quaternary volcanic rocks from eastern Sicily definitively indicate that Sicily is a part of the African plate, which collided with the European continental plate in Middle Miocene times. These data and the tectonic evolution of Sicily as inferred from the nature, age and distribution of volcanic products, are broadly consistent with the motions of Africa relative to Europe since the Upper Trias. During the Mesozoic, eastern Sicily was affected by extensional tectonics with associated alkali basaltic volcanism, and oceanic crust was produced in the meantime between the diverging African and European plates. Near the end of Mesozoic times the two plates started to converge with consequent consumption of oceanic crust. Different times of oceanic plate consumption along the Sicily-Calabria section of the plate boundary are suggested by the occurence of andesitic volcanism of different ages. The tectonic significance of late Tertiary to present basaltic activity in eastern Sicily is also discussed.     

Three-dimensional mantle lithosphere deformation at collisional plate boundaries: A subduction scissor across the South Island of New Zealand

The continental plate collision across the South Island of New Zealand is highly oblique (dextral) and bounded by oppositely verging ocean plate subduction zones. As such, the region can be considered as a type of ‘subduction scissor’. Within this tectonic context, we use three-dimensional computational geodynamic models to consider how convergent mantle lithosphere can be modified by scissor and strike–slip effects. Bounding subduction at both ends of the continental collision causes flow of the descending mantle lithosphere in the direction along strike of the model plate boundary, with thinning in the centre and thickening towards the subduction zones that bifurcates the continental mantle lithosphere root. With dipping bounding subduction, the mantle lithosphere root takes on a more complex morphology that folds over from one subduction polarity to the other, but remains as a continuous feature as it folds under the collision zone. In the absence of bounding subduction, the plate convergence causes a linear (along strike) mantle lithosphere root to develop. A rapid strike–slip motion between the converging plates transfers material in the plate boundary-parallel direction and tends to blur out features that develop in this direction—such as descending viscous instabilities. The along-strike variations in the morphology of the mantle lithosphere root that develop in the models—viz., thickening of the root towards the subduction edges, thinning in the centre—are consistent with recent, albeit poorly constrained, geophysical interpretations of the large-scale lithospheric structure of the South Island. We speculate that this reflects the nature of the evolution of the South Island collision as a limited continental segment of the plate boundary that it is dominated and guided by adjacent well-developed/developing ocean plate subduction.     

Recent geodynamic evolution of convergent plate margins and the reconstruction of fossil plate boundaries

Summary The discovery of paleoplates buried in the upper mantle leads to an interpretation of the subduction as a discontinuous process running in cycles and shifting the place of its operation in or against the direction of ocean floor spreading. This mechanism explains the distribution of calc-alkaline volcanism of different age in fossil convergent plate boundaries. The establishment of regular spatial correlation of the aseismic gap in the Wadati-Benioff zones with the distribution of calc-alkaline volcanism enables to reconstruct fossil plate boundaries and to define allochtonous terranes in apparently homogeneous continental plates. The hampering effect of the ocean floor morphology and of the fragments of continental plates approaching the trench, which substantially influences the rates of subduction and the geodynamic history of active continental margins in different domains along the trench, allows us to understand the complicated geological development of continental wedges in fossil convergent plate margins. The establishment of the segmented nature of active subduction zones and the dramatic morphology of the lower limit of the active subducted slab along the trench help us to interpret extensive lateral gaps in volcanic chains overlying active as well as fossil subduction zones.     

Nd-isotope systematics of ∼2.7 Ga adakites, magnesian andesites, and arc basalts, Superior Province: evidence for shallow crustal recycling at Archean subduction zones

An association of adakite, magnesian andesite (MA), and Nb-enriched basalt (NEB) volcanic flows, which erupted within ‘normal’ intra-oceanic arc tholeiitic to calc-alkaline basalts, has recently been documented in ∼2.7 Ga Wawa greenstone belts. Large, positive initial ?_Nd values (+1.95 to +2.45) of the adakites signify that their basaltic precursors, with a short crustal residence, were derived from a long-term depleted mantle source. It is likely that the adakites represent the melts of subducted late Archean oceanic crust. Initial ?_Nd values in the MA (+0.14 to +1.68), Nb-enriched basalts and andesites (NEBA) (+1.11 to +2.05), and ‘normal’ intra-oceanic arc tholeiitic to calc-alkaline basalts and andesites (+1.44 to +2.44) overlap with, but extend to lower values than, the adakites. Large, tightly clustered ?_Nd values of the adakites, together with Th/Ce and Ce/Yb systematics of the arc basalts that rule out sediment melting, place the enriched source in the sub-arc mantle. Accordingly, isotopic data for the MA, NEBA, and ‘normal’ arc basalts can be explained by melting of an isotopically heterogeneous sub-arc mantle that had been variably enriched by recycling of continental material into the shallow mantle in late Archean subduction zones up to 200 Ma prior to the 2.7 Ga arc. If the late Archean Wawa adakites, MA, and basalts were generated by similar geodynamic processes as their counterparts in Cenozoic arcs, involving subduction of young and/or hot ocean lithosphere, then it is likely that late Archean oceanic crust, and arc crust, were also created and destroyed by modern plate tectonic-like geodynamic processes. This study suggests that crustal recycling through subduction zone processes played an important role for the generation of heterogeneity in the Archean upper mantle. In addition, the results of this study indicate that the Nd-isotope compositions of Archean arc- and plume-derived volcanic rocks are not very distinct, whereas Phanerozoic plumes and intra-oceanic arcs tend to have different Nd-isotopic compositions.