On the role of heat flow, lithosphere thickness and lithosphere density on gravitational potential stresses

Gravitational potential stresses (GPSt) are known to play a first-order role in the state of stress of the Earth's lithosphere. Previous studies focussed mainly on crust elevation and structure and little attention has been paid to modelling GPSt using realistic lithospheric structures. The aim of the present contribution is to quantify gravitational potential energies and stresses associated with stable lithospheric domains. In order to model realistic lithosphere structures, a wide variety of data are considered: surface heat flow, chemical depletion of mantle lithosphere, crustal thickness and elevation. A numerical method is presented which involves classical steady-state heat equations to derive lithosphere thickness, geotherm and density distribution, but additionally requires the studied lithosphere to be isostatically compensated at its base. The impact of varying surface and crustal heat flow, topography, Moho depth and crust density on the signs and magnitudes of predicted GPSt is systematically explored. In clear contrast with what is assumed in most previous studies, modelling results show that the density structure of the mantle lithosphere has a significant impact on the value of the predicted GPSt, in particular in the case of thick lithospheres. Using independent information from the literature, the method was applied to get insights in the state of stress of continental domains with contrasting tectono-thermal ages. The modelling results suggest that in the absence of tectonic stresses Phanerozoic and Proterozoic lithospheres are spontaneously submitted to compression whereas Archean lithospheres are in a neutral to slightly tensile stress state. These findings are in general in good agreement with global stress measurements and observed geoid undulations.     

The thermal structure and thickness of continental roots

We compare heat flow data from the Precambrian shields in North America and in South Africa. We also review data available in other less well-sampled Shield regions. Variations in crustal heat production account for most of the variability of the heat flow. Because of this variability, it is difficult to define a single average crustal model representative of a whole tectonic province. The average heat flow values of different Archean provinces in Canada, South Africa, Australia and India differ by significant amounts. This is also true for Proterozoic provinces. For example, the heat flow is significantly higher in the Proterozoic Namaqua–Natal Belt of South Africa than in the Grenville Province of the Canadian Shield (61 vs. 41 mW m⁻² on average). These observations indicate that it is not possible to define single value of the average heat flow for all provinces of the same crustal age. Large amplitude short wavelength variations of the heat flow suggest that most of the difference between Proterozoic and Archean heat flow is of crustal origin. In eastern Canada, there is no good correlation between the local values of heat flow and heat production. In the Archean, Proterozoic and Paleozoic provinces of eastern Canada, heat flow values through rocks with the same heat production are not significantly different. There is therefore no evidence for variations of the mantle heat flow beneath these different provinces. After removing the local crustal heat production from the surface heat flow, the mantle (Moho) heat flow was estimated to be between 10–15 mW m⁻² in the Archean, Proterozoic and Paleozoic provinces of eastern Canada. Estimates of the mantle heat flow in the Kaapvaal craton of South Africa may be slightly higher (≈17 mW m⁻²). Large-scale variations of bulk crustal heat production are well-documented in Canada and imply significant differences of deep lithospheric thermal structure. In thick lithosphere, surficial heat flow measurements record a time average of heat production in the lithospheric mantle and are not in equilibrium with the instantaneous heat production. The low mantle heat flow and current estimates of heat production in the lithospheric mantle do not support a mechanical (conductive) lithosphere thinner than 200 km and thicker than 330 km. Temperature anomalies with surrounding oceanic mantle extend to the convective boundary layer below the conductive layer, and hence to depths greater than these estimates. Mechanical and thermal stability of the lithosphere require the mantle part of the lithosphere to be chemically buoyant and depleted in radiogenic elements. Both characteristics are achieved simultaneously by partial melting and melt extraction.     

Mesozoic lithospheric deformation in the North China block: Numerical simulation of evolution from orogenic belt to extensional basin system

Horizontal extension of a previously thickened crust could be the principal mechanism that caused the development of widespread extensional basins throughout the North China block (Hua-Bei region) during the Mesozoic. We develop here a regional tectonic model for the evolution of the lithosphere in the North China block, based on thin sheet models of lithospheric deformation, with numerical solutions obtained using the finite element method. The tectonic evolution of this region is defined conceptually by two stages in our simplified tectonic model: the first stage is dominated by N–S shortening, and the second by E–W extension. We associate the N–S shortening with the Triassic continental collision between the North and South China blocks, assuming that the Tan-Lu Fault system defines the eastern boundary of the North China block. The late Mesozoic E–W extension that created the Mesozoic basin systems requires a change in the regional stress state that could have been triggered by either or both of the following factors: First, gravitational instability of the lithosphere triggered by crustal convergence might have removed the lower layers of the thickened mantle lithosphere and thus caused a rapid increase in the local gravitational potential energy of the lithosphere. Secondly, a change to the constraining stress on the eastern boundary of the North China block, that might have been caused by roll-back of the subducting Pacific slab, could have reduced the E–W horizontal stress enough to activate extension. Our simulations show that widespread thickening of the North China block by as much as 50% can be explained by the collision with South China in the Triassic and Jurassic. If convergence then ceases, E–W extension can occur in the model if the eastern boundary of the region can move outwards. We find that such extension may occur, restoring crustal thickness of order 30 km within a period of 50 Myr or less, if the depth-averaged constitutive relation of the lithosphere is Newtonian, and if the Argand number (the ratio of buoyancy-derived stress to viscous stress) is greater than about 4. Widespread convective thinning of the lithosphere is not required in order to drive the extension with these parameters. If, however, the lithospheric viscosity is non-Newtonian (with strain-rate proportional to the third power of stress) the extensional phase would not occur in a geologically plausible time unless the Argand number were significantly increased by a lithospheric thinning event that was triggered by crustal thickening ratios as low as 1.5.     

Crustal and mantle strengths in continental lithosphere: is the jelly sandwich model obsolete?

The relative importance of the contribution of the lower crust and of the lithospheric mantle to the total strength of the continental lithosphere is assessed systematically for realistic ranges of layer thickness, composition, and temperature. Results are presented as relative strength maps, giving the ratio of the lower crust to upper mantle contribution in terms of crustal thickness and surface heat flow. The lithosphere shows a “jelly sandwich” rheological layering for low surface heat flow, thin to average crustal thickness, and felsic or wet mafic lower crustal compositions. On the other hand, most of the total strength resides in the seismogenic crust in regions of high surface heat flow, crust of any thickness, and dry mafic lower crustal composition.     

The evolution of the Arabian lower crust and lithospheric mantle — Geochemical constraints from southern Syrian mafic and ultramafic xenoliths

Geochemical compositions of lower crustal and lithospheric mantle xenoliths found in alkali basaltic lavas from the Harrat Ash Shamah volcanic field in southern Syria place constraints on the formation of the Arabian–Nubian Shield in northern Arabia. Compositions of lower crustal granulites are compatible with a cumulate formation from mafic melts and indicate that they are not genetically related to their host rocks. Instead, their depletion in Nb relative to other incompatible elements points to an origin in a Neoproterozoic subduction zone as recorded by an average depleted mantle Sm–Nd model age of 630 Ma.Lithospheric spinel peridotites typically represent relatively low degree (< 10%) partial melting residues of spinel lherzolite with primitive mantle compositions as indicated by major and trace element modelling of clinopyroxene and spinel. The primary compositions of the xenoliths were subsequently altered by metasomatic reactions with low degree silicate melts and possibly carbonatites. Because host lavas lack these signatures any recent reaction of the lherzolites with their host magma can be ruled out. Sm–Nd data of clinopyroxene from Arabian lithospheric mantle lherzolites yield an average age of 640 Ma suggesting that the lithosphere was not replaced since its formation and supporting a common origin of the Arabian lower crustal and lithospheric mantle sections.The new data along with published Arabian mantle xenolith compositions are consistent with a model in which the lithospheric precursor was depleted oceanic lithosphere that was overprinted by metasomatic processes related to subduction and arc accretion during the generation of the Arabian–Nubian Shield. The less refractory nature of the northern Arabian lithosphere as indicated by higher Al, Na and lower Si and Mg contents of clinopyroxenes compared to the more depleted nature of the south Arabian lithospheric mantle, and the comparable low extent of melt extraction suggest that the northern Arabian lithosphere formed in a continental arc system, whereas the lithosphere in the southern part of Arabia appears to be of oceanic arc origin.     

岩石圈伸展的壳/幔拆离模型(Parallel Extension Tectonics):华北克拉通东部早白垩世岩石圈减薄与破坏机理

华北克拉通岩石圈减薄和破坏机理长期以来存在争议,基于岩石学、岩石地球的化学分析研究突出强调深部过程的重要性。前人提出了两种重要模式:包括以拆沉作用为代表的top-down tectonics模型和以热-机械侵蚀与化学侵蚀,或地幔置换、交代作用的bottom-up tectonics模型。然而,对于这两种模式而言尚存在许多无法合理解释的问题,比如在此深部过程中,区域性岩石圈伸展有多大的贡献?地壳伸展构造是作为深部过程的响应,还是同为岩石圈伸展的产物?本文基于早白垩世东亚地区(尤其是华北克拉通东部地区)伸展构造与岩浆活动的综合分析,揭示出华北克拉通东部不同地区伸展构造变形与岩浆活动之间的时、空和成因关系有一定的差异。但整体上看,岩石圈伸展起着主导作用,控制着岩浆上侵和就位,在拆离断层下盘侵入形成各种规模的花岗质为主的侵入体,或于上盘喷发形成火山-沉积岩盆地。在伸展构造发育的不同阶段,可以有伸展早期、伸展期及伸展期后的岩浆活动。岩浆活动的强度及岩浆源区特点有显著的时空变化。一方面,在同一地区不同演化阶段其源区有很大的差异。表现为主体上是早期以古老下地壳源为主,随着壳/幔伸展作用演化,逐渐向混合源或独立幔源的演化。同时,不同地区岩浆源区的变化规律也显著不同。以胶辽地区为例,胶东整体上是壳幔混合源区对于岩浆演化有重要贡献;而辽东地区具有显著的源区演化特点:从剪切早期古老下地壳源区为主,并伴有幔源物质加入,剪切期古老下地壳为主,到剪切晚期和剪切期后以新生下地壳为主。本文认为岩石圈伸展的壳/幔拆离模型(Parallel Extension Tectonics),可以合理地解释华北克拉通及邻区早白垩世构造-岩浆活动性。在该模型中,遭受伸展的华北克拉通岩石圈发生壳-幔拆离作用。在岩石圈伸展作用期间,地壳层次的拆离作用与岩石圈地幔层次上的拆离作用可以是耦合的或者是解耦的,从而导致华北克拉通岩石圈减薄过程中在地壳尺度上的拆离作用与变质核杂岩的剥露有三种不同的类型:同岩浆活动型伸展(C型:Co-magmatism mode extension)、无岩浆活动型伸展(A型:Amagmatism mode extension)和多阶段混合型(M型:Multi-mode extension)。     

Effects of mantle upwelling in a compressional setting: the Atlas Mountains of Morocco

We discuss the implications of a lithospheric model of the Moroccan Atlas Mountains based on topography, heat flow, gravity and geoid anomalies, taking into account the regional geology. The NW African cratonic lithosphere, some 160–180 km thick, thins to c. 80 km beneath the Atlas fold-thrust belts, in contrast with the shortening regime prevailing there since the early Cenozoic. This fact explains several geological and geophysical features as high topography with modest tectonic shortening, the occurrence of alkaline magmatism contemporaneous to compression, the absence of large crustal roots to support elevation, the scarce development of foreland basins, and a marked geoid high. The modelled lithosphere thinning is related to a thermal upwelling constrained between the Iberia–Africa convergent plate boundary and the Saharan craton.     

Crustal‐mantle lithosphere decoupling as a control on Variscan metamorphism in the Eastern Alps

Mica schists of the Variscan Austroalpine Ötztal basement (Eastern Alps) contain synkinematically grown garnet porphyroblasts showing three zones with complex inclusion trails. Staurolite, kyanite and sillimanite grew over the main schistosity that envelops the garnet porphyroblasts. The mineral‐forming reaction sequences modelled in the MnKFMASH and KFMASH systems show that garnet of zone 1 grew at approximately 1300 MPa, whereas zone 2 and zone 3 grew during substantial decompression and partial thermal relaxation. For the growth of staurolite, kyanite and sillimanite, this modelling shows that thermal relaxation continued after decompression under only slightly changing pressure conditions, until the peak temperature conditions of Variscan metamorphism were reached. The result of this modelling is a pressure‐temperature‐time (P–T–t) heating path that is consistent with a tectonic interpretation implying thickening of the lithosphere during continental collision and subsequent mantle ‐ crustal lithosphere decoupling plus extension. A direct implication of this model is that the heat budget of Variscan metamorphism was controlled by the convective replacement of the thickened lithospheric root with asthenospheric material.     

东海陆架盆地西湖凹陷岩石圈热流变性质

为了探讨东海陆架盆地西湖凹陷岩石圈热流变性质,本文以实测地温数据为依据,模拟西湖凹陷岩石圈热结构,在此基础上,应用流变学原理模拟确定西湖凹陷岩石圈流变性质。结果表明,西湖凹陷岩石圈为一个冷地壳-热地幔、强地壳-弱地幔的"奶油蛋糕"型岩石圈。西湖凹陷平均地表热流密度为71 m W/m~2,地幔热流密度为40~65 m W/m~2,对地表热流密度的贡献度达73%~79%,地表热流受地幔热流控制,莫霍面温度在700℃左右,热岩石圈平均厚度为66 km。西湖凹陷岩石圈流变分层明显,上、中地壳基本为脆性层,下地壳和岩石圈上地幔为韧性层,岩石圈总流变强度平均约为2.65′10~(12) N/m,其中地壳流变强度为2.12′10~(12) N/m,地幔流变强度为5.29′10~(11) N/m,有效弹性厚度为11.7~14.5 km,地壳的流变性质控制了岩石圈的流变行为。此外,西湖凹陷岩石圈总强度较低,在构造应力作用下易于变形,且存在壳幔解耦现象。西湖凹陷岩石圈热状态及流变性质决定了西湖凹陷东部地区主要以浅部地壳的断层滑动和地层破裂来调节深部的构造应力。     

西藏及其周围地区地壳、地幔地震层析成像——印度板块大规模俯冲于西藏高原之下的证据

文中介绍有关西藏—喜马拉雅碰撞带的一项地震层析成像研究。根据一个用天然地震数据产生的全球波速模型 ,印度板块有可能以近水平状俯冲于整个西藏高原之下至 16 5～ 2 6 0km深度。西藏岩石圈具有低波速地壳和高波速下岩石圈 (75～ 12 0km深 )。在 12 0～ 16 5km深度范围 ,西藏岩石圈与俯冲的印度板块之间有一层低速软流圈物质。高原中部从地表到 310km深处有一低速体 ,说明地幔物质有可能穿过俯冲板块的脆弱部位上隆。这些结果以及野外实测的地壳缩短值说明高原的抬升得助于印度板块的近水平俯冲。我们推论俯冲印度板块的升温上浮以及上覆软流层的存在是造成西藏高原高海拔抬升以及内部地表仍相对平坦的主要原因。2 0 0 1年 1月 2 6日在印度西部发生的毁灭性大地震有可能是俯冲应力在印度板块后缘薄弱处引发的岩石圈大断裂。     

Geophysical constraints on mesozoic disruption of North China Craton by underplating‐triggered lower‐crust flow of the Archaean lithosphere

The mechanism of the disruption, both lithospheric thinning and oceanization of the commonly accepted long‐term‐stable Archaean craton, is still an open question. The available models, all imply a bottom to top process. With the construction of a 1660‐km‐long transect across the eastern North China Craton (NCC), we demonstrate that both the P‐wave velocity and density in the lowermost crust beneath the central section are significantly higher than in the corresponding parts of the south and north sections on the transect. These features are interpreted as geophysical signature of lower crustal underplating, which supplies sufficiently high gravitational potential energy to trigger lateral flow of the lower crust. This magma underplating‐triggered bilateral lower crust flow may facilitate the lithospheric thinning by means of asthenosphere upwelling and decompression melting, which infill the gap produced by the lower crust flow. The underplating‐triggered lower crustal flow can provide an alternative mechanism to explain the NCC lithosphere disruption, which highlights the crustal feedback to Archaean lithosphere disruption, from top to bottom.     

Paleolithospheric structure revealed by continental geoid anomalies

Lithospheric geoid anomalies record changes in elevation and potential energy experienced by continental lithosphere. Estimates of local isostatic equilibration and potential energy, in tandem with lithosphere-related geoid anomalies, can be used to estimate paleolithospheric thickness, providing a clearer understanding of how and why continental topography is developed. We employ several simplifying assumptions about the crustal and mantle lithosphere density and structure (and readily acknowledge that our results are therefore first-order approximations) to predict the pre-orogenic structure of the lithosphere. At the outset we emphasize that while this approach does not provide an exhaustive evaluation of the deformation mechanism, it does serve to quantify the relative role played by the variations in the crustal and upper mantle components of the lithosphere. In this way we are able to use independent measurement of lithospheric geoid anomalies, current (post-orogenic) elevation and lithospheric structure, and paleoelevation information to estimate topographic development and structural support over time. Application of this technique to the southwestern United States indicates that the uplift of the Colorado Plateau is the result of processes in both the crust and mantle lithosphere and that the lithosphere of the pre-orogenic Southern Basin and Range was thinned relative to the Northern Basin and Range and Colorado Plateau. Although we use the southwestern U.S. as an example, this method can help constrain uplift mechanisms for any region for which the structure and geoid anomaly of the modern lithosphere is well understood.     

Gravity models of two-level collision of lithospheric plates in northeastern Asia

Structural forms of emplacement of crustal and mantle rigid sheets in collision zones of lithospheric plates in northeastern Asia are analyzed using formalized gravity models reflecting the rheological properties of geological media. Splitting of the lithosphere of moving plates into crustal and mantle constituents is the main feature of collision zones, which is repeated in the structural units irrespective of their location, rank, and age. Formal signs of crustal sheet thrusting over convergent plate boundaries and subduction of the lithospheric mantle beneath these boundaries have been revealed. The deep boundaries and thickness of lithospheric plates and asthenospheric lenses have been traced. A similarity in the deep structure of collision zones of second-order marginal-sea buffer plates differing in age is displayed at the boundaries with the Eurasian, North American, and Pacific plates of the first order. Collision of oceanic crustal segments with the Mesozoic continental margin in the Sikhote-Alin is characterized, as well as collision of the oceanic lithosphere with the Kamchatka composite island arc. A spatiotemporal series of deep-seated Middle Mesozoic, Late Mesosoic, and Cenozoic collision tectonic units having similar structure is displayed in the transitional zone from the Asian continent to the Pacific plate.     

Geochemistry of the Cenozoic Potassic Volcanic Rocks in the West Kunlun Mountains and Constraints on Their Sources

The geochemical characteristics of the Cenozoic volcanic rocks from the north Pulu, east Pulu and Dahongliutan regions in the west Kunlun Mountains are somewhat similar as a whole. However, the volcanic rocks from the Dahongliutan region in the south belt are geochemically distinguished from those in the Pulu region (including the north and east Pulu) of the north belt. The volcanic rocks of the Dahongliutan region are characterized by relatively low TiO2 abundance, but more enrichment in alkali, much more enrichment in light rare earth elements and large ion lithosphile elements than those from the Pulu region. Compared with the Pulu region, volcanic rocks from the Dahongliutan region have relatively low 87Sr/86Sr ratios, and high εNd, 207Pb/204Pb and 208Pb/204Pb. Their trace elements and isotopic data suggest that they were derived from lithospheric mantle, consisting of biotite- and hornblende-bearing garnet lherzolite, which had undertaken metasomatism and enrichment. On the primitive mantle-normali     

由震源机制和地震波各向异性探讨青藏高原岩石圈变形

本文据青藏高原天然地震震源参数和地震波各向异性资料，讨论了高原岩圈不同圈层的变形特性。     

A rapid method to map the crustal and lithospheric thickness using elevation, geoid anomaly and thermal analysis. Application to the Gibraltar Arc System, Atlas Mountains and adjacent zones

We present a method based on the combination of elevation and geoid anomaly data together with thermal field to map crustal and lithospheric thickness. The main assumptions are local isostasy and a four-layered model composed of crust, lithospheric mantle, sea water and the asthenosphere. We consider a linear density gradient for the crust and a temperature dependent density for the lithospheric mantle. We perform sensitivity tests to evaluate the effect of the variation of the model parameters and the influence of RMS error of elevation and geoid anomaly databases. The application of this method to the Gibraltar Arc System, Atlas Mountains and adjacent zones reveals the presence of a lithospheric thinning zone, SW–NE oriented. This zone affects the High and Middle Atlas and extends from the Canary Islands to the eastern Alboran Basin and is probably linked with a similarly trending zone of thick lithosphere constituting the western Betics, eastern Rif, Rharb Basin, and Gulf of Cadiz. A number of different, even mutually opposite, geodynamic models have been proposed to explain the origin and evolution of the study area. Our results suggest that a plausible slab-retreating model should incorporate tear and asymmetric roll-back of the subducting slab to fit the present-day observed lithosphere geometry. In this context, the lithospheric thinning would be caused by lateral asthenospheric flow. An alternative mechanism responsible for lithospheric thinning is the presence of a hot magmatic reservoir derived from a deep ancient plume centred in the Canary Island, and extending as far as Central Europe.     

Stress and strain during fault-controlled lithospheric extension—insights from numerical experiments

Two-dimensional finite element techniques are used to study the temporal evolution and spatial distribution of stress and strain during lithospheric extension. The thermomechanical model includes a pre-existing fault in the upper crust to account for the reactivation of older tectonic elements. The fault is described using contact elements which allow for independent meshing of hanging wall and foot wall as well as simulation of large differential displacements between the fault blocks. Numerical models are run for three different initial temperature distributions representing extension of weak, moderately strong and strong lithosphere and three different extension velocities. In spite of the simple geodynamic boundary conditions selected, i.e., wholesale extension at a constant rate, stress and strain vary substantially throughout the lithosphere. In particular, in case of the weak lithosphere model, lower crustal flow towards the locus of maximum upper crustal extension results in the formation of a lower crustal dome while maintaining a subhorizontal Moho relief. The core of the dome experiences hardly any internal deformation, although it is the part of the lower crust which is exhumed the most. Stress fields in the lower crustal dome vary significantly from the regional trend underlining mechanical decoupling of the lower crust from the rest of the lithosphere. These differences diminish if cooler temperatures and, hence, stronger rheologies are considered. Lithospheric strength also exerts a profound control on the basin architecture and the surface expressions of extension, i.e., rift flank uplift and basin subsidence. If the lower crust is sufficiently weak, its flow towards the region of extended upper crust can provide a threshold value for the maximum subsidence which can be achieved during the syn-rift stage. In spite of continuous regional extension, corresponding burial history plots show exponentially decreasing subsidence rates which would traditionally be interpreted in terms of lithospheric cooling during the post-rift stage. The models provide templates to genetically link the surface and sub-surface expressions of lithospheric extension, for which usually no contemporaneous observations are possible. In particular, they help to decipher the information on the physical state of the lithosphere at the time of extension which is stored in the architecture and subsidence record of sedimentary basins.     

中国东部壳-幔、岩石圈-软流圈之间的相互作用带：特征及转换时限

中国东部中—新生代,下部岩石圈存在壳与幔、岩石圈与软流圈两个相互作用带,它们是重要的岩浆源区,在层圈相互作用中,热和物质的交换及其动力学过程是引起中、新生代岩石圈内部层圈间的厚度调整、岩石圈不均匀减薄以及区域构造-岩浆-成矿作用的重要机理。大陆内部的壳-幔作用有3种类型:地幔来源的底侵熔体与下地壳的作用;下地壳拆沉进入弱化(weakening)了的岩石圈地幔二者发生的作用以及陆-陆碰撞深俯冲带的壳-幔相互作用。它们形成的火山岩组合有一定的差别,但源区都含有地壳组分。岩石圈-软流圈的作用带也是重要的岩浆源区,源区是以软流圈地幔为主,基本不含地壳组分。东部岩石圈的减薄时间大体与出现大规模软流圈来源的玄武岩喷发的时间一致,也与上述两类层圈作用转换的时间一致,大约在100Ma以后。     

Stagnant lids and mantle overturns:Implications for Archaean tectonics,magmagenesis,crustal growth,mantle evolution,and the start of plate tectonics

The lower plate is the dominant agent in modern convergent margins characterized by active subduction,as negatively buoyant oceanic lithosphere sinks into the asthenosphere under its own weight.This is a strong plate-driving force because the slab-pull force is transmitted through the stiff sub-oceanic lithospheric mantle.As geological and geochemical data seem inconsistent with the existence of modernstyle ridges and arcs in the Archaean,a periodically-destabilized stagnant-lid crust system is proposed instead.Stagnant-lid intervals may correspond to periods of layered mantle convection where efficient cooling was restricted to the upper mantle,perturbing Earth's heat generation/loss balance,eventually triggering mantle overturns.Archaean basalts were derived from fertile mantle in overturn upwelling zones(OUZOs),which were larger and longer-lived than post-Archaean plumes.Early cratons/continents probably formed above OUZOs as large volumes of basalt and komatiite were delivered for protracted periods,allowing basal crustal cannibalism,garnetiferous crustal restite delamination,and coupled development of continental crust and sub-continental lithospheric mantle.Periodic mixing and rehomogenization during overturns retarded development of isotopically depleted MORB(mid-ocean ridge basalt)mantle.Only after the start of true subduction did sequestration of subducted slabs at the coremantle boundary lead to the development of the depleted MORB mantle source.During Archaean mantle overturns,pre-existing continents located above OUZOs would be strongly reworked;whereas OUZOdistal continents would drift in response to mantle currents.The leading edge of drifting Archaean continents would be convergent margins characterized by terrane accretion,imbrication,subcretion and anatexis of unsubductable oceanic lithosphere.As Earth cooled and the background oceanic lithosphere became denser and stiffer,there would be an increasing probability that oceanic crustal segments could founder in an organized way,producing a gradual evolution of pre-subduction convergent margins into modern-style active subduction systems around 2.5 Ga.Plate tectonics today is constituted of:(1)a continental drift system that started in the Early Archaean,driven by deep mantle currents pressing against the Archaean-age sub-continental lithospheric mantle keels that underlie Archaean cratons;(2)a subduction-driven system that started near the end of the Archaean.