Episodic tectonic plate reorganizations driven by mantle convection

Periods of relatively uniform plate motion were interrupted several times throughout the Cenozoic and Mesozoic by rapid plate reorganization events [R. Hey, Geol. Soc. Am. Bull. 88 (1977) 1404-1420; P.A. Rona, E.S. Richardson, Earth Planet. Sci. Lett. 40 (1978) 1-11; D.C. Engebretson, A. Cox, R.G. Gordon, Geol. Soc. Am. Spec. Pap. 206 (1985); R.G. Gordon, D.M. Jurdy, J. Geophys. Res. 91 (1986) 12389-12406; D.A. Clague, G.B. Dalrymple, US Geol. Surv. Prof. Pap. 1350 (1987) 5-54; J.M. Stock, P. Molnar, Nature 325 (1987) 495-499; C. Lithgow-Bertelloni, M.A. Richards, Geophys. Res. Lett. 22 (1995) 1317-1320; M.A. Richards, C. Lithgow-Bertelloni, Earth Planet. Sci. Lett. 137 (1996) 19-27; C. Lithgow-Bertelloni, M.A. Richards, Rev. Geophys. 36 (1998) 27-78]. It has been proposed that changes in plate boundary forces are responsible for these events [M.A. Richards, C. Lithgow-Bertelloni, Earth Planet. Sci. Lett. 137 (1996) 19-27; C. Lithgow-Bertelloni, M.A. Richards, Rev. Geophys. 36 (1998) 27-78]. We present an alternative hypothesis: convection-driven plate motions are intrinsically unstable due to a buoyant instability that develops as a result of the influence of plates on an internally heated mantle. This instability, which has not been described before, is responsible for episodic reorganizations of plate motion. Numerical mantle convection experiments demonstrate that high-Rayleigh number convection with internal heating and surface plates is sufficient to induce plate reorganization events, changes in plate boundary forces, or plate geometry, are not required.     

Plate-kinematic Constraints from Mantle Bénard Convection

General kinematic implications for plate tectonics are determined for Rayleigh-Bénard convection of the mantle. The continuum of all possible configurations of Bénard polygons is probed by large random samples of global configurations (450,000 to 54,000,000), for each of which the Euler poles are determined on the basis of viscous coupling across the asthenosphere. Two computationally related methods lead first, to Euler pole restrictions for fourteen plates, and second, to restrictions on the Bénard cell configuration. Result No. 1: Euler poles occur in global preference-patterns, which are determined exclusively by the shape of the plate. The observational HS2-NUVEL1 model poles occur near regions preferred by Bénard convection (Eurasia excluded); the agreement is best for the most accurate observational poles. Result No. 2: Seven specific mantle Bénard cells are indicated by present-day plate motions. The upwelling centers correlate with hotspot domains; the major global subduction zones correlate with Bénard model downwelling. This result is independent of the Euler pole accuracy used in its determination, and is consistent with the distribution of low seismic p-wave propagation velocities determined by tomography, and with shear-wave splitting analysis within the asthenosphere. Conclusions: The results suggest that the bulk mantle is divided into less than ten Bénard convection cells globally (cf., Fohlmeister and Renka, 2002), each of which extends from the asthenosphere to the core-mantle boundary; turbulent flow, and other perturbations of the Bénard kinematics appear to be limited. These primally poloidal flow kinematics provide basal shear forces as a major component in driving plate tectonics, and are specifically configured for the directions of plate motions. The Bénard model is incomplete without a dynamic contribution from the lithosphere, which represents a separate convection layer of the distinct polar kinematics of rigid plates. The complete hybrid mechanism for driving plate tectonics includes lithospheric buoyancy dynamics, specifically from the subducting Pacific plate slabs to compensate for plate-slowing due to the back-flow sector of the Hawaiian convection cell, and collision-drag dynamics principally for smaller plates or continental margins.     

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.     

现今的地幔流动及其与板块构造的相互作用:从观测证据得到的推论

地震波速度结构层析成像和地震各向异性分析,是推测现今地幔流动的主要观测依据．从已有研究结果看,全球尺度的地幔流动的两个主要边界驱动力是顶部的冷却和底部的加热,地幔的密度和粘度控制流动速率．沿海沟由消减板带动的地幔下沉,和沿热点下面幔柱及洋中脊的地幔上升,是地幔垂直向流动的表现．GPS等测量显示的全球板块运动在一定条件下反映地幔顶部的水平流动．地幔柱可能有不同的根源深度,沿幔柱上升的地幔流在200—350km深度转变为水平流动．消减带附近有复杂的地幔流动格局,表明局部构造条件对地幔流动的影响．大陆下一般出现两个地幔各向异性层,较深的可能反映地幔流动,并与大陆根的状态有关．在不同构造环境下,地幔流动与板块构造之间有不同形态的相互作用关系,它可能驱动板块运动,也可能对板块运动产生阻力。     

An approximate method for determining the stability of two-scale flow in the mantle

Summary An approximate method of calculating the stability of two-dimensional convection rolls to cross-roll disturbances is evaluated by comparison with the results from exact analyses and good agreement is obtained. In particular, investigation of the second mode of the disturbance provides a qualitative estimate of the terms excluded in making the approximation. We conclude that this approximate method, sensibly used, gives a good indication of the stability of rolls to cross-roll disturbances.Application of this method is made to convection rolls heated partially or wholly from within, and to rolls in the presence of a longitudinal shear flow. The approximation indicates that, at high Rayleigh Numbers, growth rates are underestimated and so the amount of shear calculated to stabilise longitudinal rolls is a lower bound. Our results suggest that two-dimensional rolls are unlikely to be a stable flow pattern beneath even the fastest moving plates.     

地幔对流拖曳力对中国大陆岩石层变形的影响

采用较为符合实际岩石层变形的非线性幂指数本构关系,基于ANSYS有限元平台, 模拟了近20万年来中国大陆地区地表运动及演化过程,探讨了印度板块挤压作用和地幔对流拖曳力各自对于中国大陆地区地表形变运动格局的影响.模拟结果与观测数据的比较表明：在印度板块的挤压和地幔拖曳力联合作用下,中国及东亚大陆岩石层运动形变模式能够和现代GPS观测有较好的吻合; 印度大陆和欧亚大陆的碰撞以及印度大陆的持续向北推进、挤压所产生的应力环境,一直主导了以青藏高原为核心的我国西部地域岩石圈构造、运动和演化,但其影响随着远离青藏高原地区而逐渐变小;地幔对流产生的作用于岩石层底部的拖曳力是中国大陆（特别是远离碰撞带）岩石层运动构造变形的重要驱动力.然而在构造复杂和东部靠近太平洋板块的区域,模型预测结果和GPS观测还存在一定的差距,这说明在未来的中国大陆岩石层变形运动的数值模拟中,应当采用更为复杂的构造模型和驱动力因素.     

Seismic gaps and plate tectonics: Seismic potential for major boundaries

The theory of plate tectonics provides a basic framework for evaluating the potential for future great earthquakes to occur along major plate boundaries. Along most of the transform and convergent plate boundaries considered in this paper, the majority of seismic slip occurs during large earthquakes, i.e., those of magnitude 7 or greater. The concepts that rupture zones, as delineated by aftershocks, tend to abut rather than overlap, and large events occur in regions with histories of both long- and short-term seismic quiescence are used in this paper to delineate major seismic gaps.In detail, however, the distribution of large shallow earthquakes along convergent plate margins is not always consistent with a simple model derived from plate tectonics. Certain plate boundaries, for example, appear in the long term to be nearly aseismic with respect to large earthquakes. The identification of specific tectonic regimes, as defined by dip of the inclined seismic zone, the presence or absence of aseismic ridges and seamounts on the downgoing lithospheric plate, the age contrast between the overthrust and underthrust plates, and the presence or absence of back-arc spreading, have led to a refinement in the application of plate tectonic theory to the evaluation of seismic potential.The term seismic gap is taken to refer to any region along an active plate boundary that has not experienced a large thrust or strike-slip earthquake for more than 30 years. A region of high seismic potential is a seismic gap that, for historic or tectonic reasons, is considered likely to produce a large shock during the next few decades. The seismic gap technique provides estimates of the location, size of future events and origin time to within a few tens of years at best.The accompanying map summarizes six categories of seismic potential for major plate boundaries in and around the margins of the Pacific Ocean and the Caribbean, South Sandwich and Sunda (Indonesia) regions for the next few decades. These categories range from what we consider high to low potential for being the site of large earthquakes during that period of time. Categories 1, 2 and 6 define a time-dependent potential based on the amount of time elapsed since the last large earthquake. The remaining categories, 3, 4, and 5, are used for areas that have ambiguous histories for large earthquakes; their seismic potential is inferred from various tectonic criteria. These six categories are meant to be interpreted as forecasts of the location and size of future large shocks and should not be considered to be predictions in which a precise estimate of the time of occurrence is specified.Several of the segments of major plate boundaries that are assigned the highest potential, i.e., category 1, are located along continental margins, adjacent to centers of population. Some of them are hundreds of kilometers long. High priority should be given to instrumenting and studying several of these major seismic gaps since many are now poorly instrumented. The categories of potential assigned here provide a rationale for assigning prorities for instrumentation, for future studies aimed at predicting large earthquakes and for making estimates of tsunami potential.Lamont-Doherty Geological Observatory Contribution No. 2906.     

The influence of thermochemical convection on the fixity of mantle plumes

A general feature of both isochemical and thermochemical studies of mantle convection is that horizontal plume velocities tend to be smaller than typical convective velocities, however, it is not clear which system leads to a greater fixity of mantle plumes. We perform two- and three-dimensional numerical calculations and compare both thermochemical and isochemical cases with similar convective vigor to determine whether presence of a dense component in the mantle can lead to smaller ratios of horizontal plume velocity to surface velocity. We investigate different viscosity and density contrasts between chemical components in the thermochemical calculations, and we perform isochemical calculations with both free-slip and no-slip bottom boundary conditions. We then compare both visually and quantitatively the results of the thermochemical and isochemical calculations to determine which leads to greater plume fixity. We find that horizontal plume velocities for thermochemical calculations are similar to those from isochemical calculations with no-slip bottom boundary conditions. In addition, we find that plumes tend to be more fixed for isochemical cases with free-slip bottom boundary conditions for two-dimensional calculations, however, in three dimensions, we find that plume fixity is similar to that observed in thermochemical calculations.     

Upper mantle convection driving by density anomaly and a test model

Introduction Richter and Mckenzie (1978) supposed that there is a small-scale convection system in the mantle. For a long time lots of research provides observational data to infer the possibility of a small-scale convection in the upper mantle. For example, Haxby and Weissel (1986) discussed the relationship between SEASAT map and small-scale convection. Baudry and Kroenke (1991), Maia and Diament (1991) found that the geoid and bathymetry exhibit peaks in the 400~650 km range in the Pa…     

Possible effects of lateral viscosity variations induced by plate-tectonic mechanism on geoid inferred from numerical models of mantle convection

In a traditional analytical method, the convective features of Earth’s mantle have been inferred from surface signatures obtained by the geodynamic model only with depth-dependent viscosity structure. The moving and subducting plates, however, bring lateral viscosity variations in the mantle. To clarify the effects of lateral viscosity variations caused by the plate-tectonic mechanism, I have first studied systematically instantaneous dynamic flow calculations using new density-viscosity models only with vertical viscosity variations in a three-dimensional spherical shell. I find that the geoid high arises over subduction zones only when the vertical viscosity contrast between the upper mantle and the lower mantle is O(10³) to O(10⁴), which seems to be much larger than the viscosity contrast suggested by other studies. I next show that this discrepancy may be removed when I consider the lateral viscosity variation caused by the plate-tectonic mechanism using two-dimensional numerical models of mantle convection with self-consistently moving and subducting plates, and suggest that the observed geoid anomaly on the Earth’s surface is significantly affected by plate-tectonic mechanism as a first-order effect.     

Mass heterogeneities and convection in the earth's mantle inferred from gravity and core-mantle boundary irregularities

Mass heterogeneities in the earth's mantle are retrieved from the gravity data and the topography of the core-mantle boundary as well as the topography of the earth's surface. A mantle circulation induced by the heterogeneities is modelled by solving the Stokes problem for incompressible Newtonian fluid. The derived models of mantle motions correlate well with the plate tectonics and point at a close relation between the surface tectonic activity and the processes in the vicinity of the core-mantle boundary.     

Transitions in the style of mantle convection at high Rayleigh numbers

The pattern and style of mantle convection govern the thermal evolution, internal dynamics, and large-scale surface deformation of the terrestrial planets. In order to characterize the nature of heat transport and convective behaviour at Rayleigh numbers, Ra, appropriate for planetary mantles (between 10⁴ and 10⁸), we perform a set of laboratory experiments. Convection is driven by a temperature gradient imposed between two rigid surfaces, and there is no internal heating. As the Rayleigh number is increased, two transitions in convective behaviour occur. First we observe a change from steady to time-dependent convection at Ra≈10⁵. A second transition occurs at higher Rayleigh numbers, Ra≈5×10⁶, with large-scale time-dependent flow being replaced by isolated rising and sinking plumes. Corresponding to the latter transition, the exponent β in the power law relating the Nusselt number Nu to the Rayleigh number (NuRa^β) is reduced. Both rising and sinking plumes always consist of plume heads followed by tails. There is no characteristic frequency for the formation of plumes.     

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.     

Three-dimensional numerical simulations of crustal deformation and subcontinental mantle convection

3-D simulations of mantle convection allowing for continental crust are explored to study the effects of crustal thickening on lithosphere stability and of continents on large-scale mantle flow. Simulations begin with a crustal layer within the upper thermal boundary layer of a mantle convection roll in a 1 × 1 × 1 Cartesian domain. Convective stresses cause crust to thicken above a sheet-like mantle downwelling. For mild convective vigor an initial crustal thickness variation is required to induce 3-D lithospheric instability below the zone of crustal convergence. The amplitude of the required variation decreases with increasing convective vigor. Morphologically, instability is manifest in formation of drip-like thermals that exist within the large-scale roll associated with initial crustal thickening. A strong surface signature of the drips is their ability to cause deviations from local Airy compensation of topography. After the initial thickening phase, the crustal accumulation that forms serves as a model analog to a continent. Its presence leads to mantle flow patterns distinctly different from the steady-state roll that results in its absence. Large lateral thermal gradients are generated at its edge allowing this region to be the initiation site for continued small-scale thermal instabilities. Eventually these instabilities induce a restructuring of large-scale mantle flow, with the roll pattern being replaced by a square cell. Although preliminary and idealized, the simulations do show the fluid dynamical plausibility behind the idea that significant mantle variations can be generated along the strike of a largely 2-D mountain chain by the formation of the chain itself. The ability of a model continent to cause a change in fundamental convective planform also suggests that the effects of continental crust on mantle convection may be low-order despite the seemingly trivial volume of crust relative to mantle.     

地幔对流与深部物质运移研究的新进展

现代固体地球科学已经认识到,地幔对流不再是少数动力学家的假想,它是地幔热动力系统的主要构架.地幔对流和板块运动驱动机理关系的研究已经从简单的主动或被动驱动的讨论转向对统一热动力系统的探讨.包括地幔热柱在内的地幔对流的深入研究不仅成为研究地幔热动力系统演化的主线,也成为研究大陆形成和演化驱动机理的主线.与此同时,以地震层析成像为主体的地震、地球物理观测资料和以地幔岩石化学组份为主体的地球化学观测成为认识地幔对流的强有力的工具.然而,地球化学和地球物理观测之间存在明显的差异,一些依赖于地球化学数据构思的新的热动力学框架对地幔对流的研究构成了强烈的挑战.     

One-dimensional modelling of mantle flow

A one-dimensional model of flow between a fixed boundary at the bottom and a moving one on top with no net flow through vertical sections is tested for geophysically interesting mantle viscosity-depth functions. Such a model, although simplistic, may help in answering the question to what depth the return flow extends, at least in the case of moving plates measuring many thousand kilometers across, such as the Pacific plate.It the viscosity in the asthenosphere is less than three orders of magnitude smaller than that of the mantle below, the return flow extends to great depth and the asthenosphere is a zone of concentrated shear. If the viscosity contrast is greater, the return flow is concentrated in the asthenosphere. For a wide range of model parameters typical flow velocities below the asthenosphere are about one-tenth of the plate velocity. The pressure gradient required by the mantle flow may be manifest in gravity trends across moving plates, but no excessive gravity anomalies are required by the model if the absolute viscosity values conform to those inferred from post-glacial rebound data. A thinner and lower-viscosity layer is favored over a thicker and more viscous layer if both fit glacial rebound evidence. The present model may not be applicable if down to the core the viscosity is as low as about 10²¹ N s m^–2 with a free-slip bottom boundary.     

Small-scale upper mantle convection and orogeny of Tianshan Mountains in China

In this study, from the travel time data recorded in the Tianshan passive seismic array experiment, we present the P-wave velocity structure of the upper mantle down to 660 km along the Kuqa-Kuitun pro-file in terms of seismic tomography technique. Based on the P-wave velocity model, we derive the corresponding 2D upper mantle density model. The 2D small-scale convection of the upper mantle underneath the Tianshan Mountains in China driven by the density anomalies is simulated using the hybrid finite element method combining with the marker-in-cell technique. The main features of the upper mantle convection and the reciprocation between the convection and mountain building are in-vestigated. The results manifest that (1) in the upper mantle underneath the Junggar basin and North Tianshan exists a counterclockwise convection, which scale is ~ 500 km; (2) underneath the Tarim ba-sin and South Tianshan exists a clockwise northward convection, which is relatively weak; (3) the convective velocity at the top of the upper mantle underneath the Tianshan Mountains in China should not be less than 20 mm/a, while considering the dependent of convective velocity on the viscosity; (4) the northward extrusion of the Tarim block plays a key role in the Cenozoic Tianshan mountain building and the present-day tectonic deformation of the Tianshan range is related closely to the upper mantle convection; and (5) the northward subduction of the Tarim block does not influence obviously the up-per mantle convection.     

横向黏度变化的全地幔对流应力场初步研究

将地幔地震波速度异常转换为地幔横向黏度变化(达到3个数量级),在球坐标系下计算了瑞雷数为106、上边界为刚性、下边界为应力自由等温边界条件下的岩石层底部的地幔对流极型和环型应力场.结果表明,地幔对流极型应力场与地表大尺度构造具有良好的对应关系:俯冲带和碰撞带的应力呈现挤压状态,而洋中脊处的应力则呈现拉张状态.地幔对流环...     

Upper mantle convection beneath northwest China and its adjacent region driven by density anomaly

Introduction The northwest of China includes Tarim, Junggar and Qaidam basins, and Kunlun, Tianshan, Altun and Qilian mountains, as well as the north part of the Tibetan Plateau. For a long time, the study of lithosphere structures and dynamics in this area has been a popular topic in geoscience, and has yielded many results. For example, TANG (1994) and LI et al (1998) suppose that the tectonic structure of Tarim Basin is various with geological periods, which changes many times betwee…     

天山南北地块构造演化与地幔对流耦合动力机制

天山造山带南北分别于塔里木盆地和准噶尔盆地相接,经历古生代时期超级大陆裂解、南北天山洋裂开、洋盆持续扩张、洋壳俯冲消减、陆陆碰撞缝合过程及中新生代陆内再造山构造调整,是现今世界上较为活跃的陆内造山带,成为国内外大陆动力学研究的热点地带.在综合分析地质学、地球物理(地震剖面、重力异常、地震层析)、地球化学、岩石学及天文学...