地幔对流与岩石层板块的相互耦合及影响──（Ⅱ）地幔混合对流理论及其应用

观测到的板块运动包含两种能量分布几乎相等的运动形态：极型场和环型场．纯粹的由热驱动的地幔自由对流不能预期和解释环型运动的产生．本文提出地幔混合对流理论，既考虑了热驱动的自由对流，也考虑了由板块自身激发的强迫对流．根据板块处于动力学平衡状态的观测事实，建立了相应的模型．数值结果表明，根据混合对流模型所预期的板块速度场，既能产生极型场，也能产生环型场，而且在空间分布特征及功率谱分布上与观测资料符合相当好．地幔物质的上升流动基本和洋脊对应，而下降流动和俯冲带对应．     

板块绝对运动及地幔热对流

本文以板块绝对运动AM1-2模型为边界条件探讨了不同的瑞利数下地幔热对流模型.结果表明,瑞利数小于10000（529.41）时,地幔对流呈现以板块驱动图式,运动的极型场和环型场由板块运动激发,两种场占有差不多相同的功率.当瑞利数增加到接近或略超过最低临界值时（约1.5倍）,对流呈现出复杂状态:1.板块运动速率小于下伏地幔对流速率;2.区域性的双层对流环出现;3.对流谱成分发生变化;4.环型场仅在地幔很浅的区域中起作用,而在地幔深部对流图式中影响很小.     

地幔对流环型场的激发(Ⅰ)

讨论了地幔内部的粘滞度及施加在地表和CMB 的边界条件对地幔对流环型场的激发分析表明,当粘滞度侧向均匀时,环型场与极型场自然解耦,且环型场不影响重力位,当粘滞度侧向不均时,环型场与极型场耦合在一起.两者共同影响重力位.当引入板块运动速度时,边界条件非零,也能激发环型场;对侧向均匀粘滞度地幔,零边界条件不能激发环型场     

地幔对流与岩石层板块的相互耦合及影响──（Ⅰ）球腔中的自由热对流

板块运动是地幔对流的主要证据之一．同时，作为地球动力系统中一个相对独立部分，板块自身的存在和运动对地幔内部物质的流动形态有巨大影响．地幔内部的流动由两部分组成：一是由内部非绝热温度差异造成的自由对流解；另一部分是由在地表运动的板块所激发．作为系列工作的第一部分，本文研究球腔中的自由热对流问题．得到了对地幔对流研究有实际意义的下边界为自由、上边界为刚性情况下的临界瑞利数值，不同的瑞利数时球腔内流场和温度场的分布形态等．     

地震波速度约束的大横向黏度变化的地幔对流研究

本文将横向黏度变化提高到3个量级,获得了地震波速度结构约束下的大横向黏度变化的地幔浅部的极型场和环型场对流图像.与小横向黏度变化下的结果相比,本文的结果具有显著的改善.对极型场对流图像,主要体现在本文结果能更清楚地解释太平洋板块、大洋洲和南美洲以及东太平洋洋中脊处的现今运动状态;对环型场对流图像,能更合理地解释北太平洋...     

三维变粘度地幔对流的理论研究

在朱涛等的"球层中的非线性自由热对流-变粘度模型"研究中(2004),我们建立了变粘度地幔对流模型,不仅获得了描述地幔物质垂直运移特征的极型场,而且获得了描述地幔物质水平运移的环型场.在计算中,假定在常粘度背景下的小粘度横向扰动形式非常简单,仅仅随着纬度而变化,其它参数取值非常接近真实地球,依此从理论上获得了变粘度情形下环型场的某些认识.但是,它还不能被直接用来解释或联系实际的地球构造运动及其分布,因为在模型中还没有引入任何实际地球物理观测资料来约束模型.建立模型的目的是希望它能为合理认识和理解现今的地球构造现象及其动力学演化过程提供帮助,为此,必须使得地幔对流模型更加合理、更加接近真实情况.一般情况下,主要采用两种方式将地球物理资料引入地幔对流模型:     

三维地震波速结构约束下的变黏度地幔对流及其动力学意义

建立了三维黏度扰动下的变黏度地幔对流模型,并提供了在引入地幔的三维地震波速度结构下相应的求解方法. 依此反演了瑞利数Ra = 106时,两种不同边界条件下的极、环型场对流图像,这有助于深化对地幔物质流动和大地构造运动的深部动力学过程的认识和理解. 研究结果表明,不但地幔浅部的极型场对流图像显示出了与大地构造运动的相关性并揭示了其深部动力学过程,更重要的是,地幔浅部的环型场对流图像首次为我们认识和理解板块构造的水平与旋转运动提供了重要的信息：环型场速度剖面中在赤道附近存在一条大致南东东—北西西向的强对流条带,可能与环赤道附近大型剪切带的形成相关,进而表明可能是该带强震发生的深部动力学背景;在南北半球存在的旋转方向相反的对流环表明它们整体上可能存在差异旋转.     

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

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

地幔对流

地幔对流是地幔中．特别是地幔软流层中发生的热对流。地幔对流是一种自然对流．既是发生在地幔中的一种传热方式（通过物质运动传递热量）．又是一种地幔物质的运动过程（由物质内部密度差或温度差所驱使的）．是地球内部向地球表面输送能量、动量和质量的一种有效途径。地幔对流是一个复杂的系统．是在缓慢的进行的,对流活动的时间可达几千万年,甚至几亿年。地幔对流的流动形态可以不同。     

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

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

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.     

The patterns of high-degree thermal free convection and its features in a spherical shell

Introduction Mantle convection is thought to be the most effective way of heat transportation in the earth and the source of driving lithospheric plates (Elasser, 1971). The velocity field of plate motion can be split into poloidal and toroidal parts, which are corresponding to the vertical (i.e. radial) and horizontal motions, respectively, in global model. The toroidal component is manifested in the existences of transform faults and the poloidal one in the presences of convergence and diver…     

Small-scale magnetic buoyancy and magnetic pumping effects in a turbulent convection

We determine the nonlinear drift velocities of the mean magnetic field and nonlinear turbulent magnetic diffusion in a turbulent convection. We show that the nonlinear drift velocities are caused by three kinds of the inhomogeneities; i.e., inhomogeneous turbulence, the nonuniform fluid density and the nonuniform turbulent heat flux. The inhomogeneous turbulence results in the well-known turbulent diamagnetic and paramagnetic velocities. The nonlinear drift velocities of the mean magnetic field cause the small-scale magnetic buoyancy and magnetic pumping effects in the turbulent convection. These phenomena are different from the large-scale magnetic buoyancy and magnetic pumping effects which are due to the effect of the mean magnetic field on the large-scale density stratified fluid flow. The small-scale magnetic buoyancy and magnetic pumping can be stronger than these large-scale effects when the mean magnetic field is smaller than the equipartition field. We discuss the small-scale magnetic buoyancy and magnetic pumping effects in the context of the solar and stellar turbulent convection. We demonstrate also that the nonlinear turbulent magnetic diffusion in the turbulent convection is anisotropic even for a weak mean magnetic field. In particular, it is enhanced in the radial direction. The magnetic fluctuations due to the small-scale dynamo increase the turbulent magnetic diffusion of the toroidal component of the mean magnetic field, while they do not affect the turbulent magnetic diffusion of the poloidal field.     

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.     

Hotspots,mantle convection and plate tectonics: A synthetic calculation

A model is developed that unifies vigorous hotspots with global-scale mantle convection and plate tectonics. The convection dynamics are assumed to generate flow patterns that emerge as closely packed polygonal cells in approaching the asthenosphere, and whose geometry is completely determined by a defining set of vigorous hotspots. Overlying viscously coupled rigid plates are driven with unique velocities (Euler vectors) at which the area integral of the shear forces is zero; these velocities are dynamically stable. The computed plate velocities, resulting from convection based on 15 hotspots, are compared with the velocities of plate motion models AM1-2 (Minster andJordan, 1978) and HS-NUVEL1 (Gripp andGordon, 1990), which combine transform fault geometries, magnetic anomalies and seismic data. The comparison shows a striking agreement for a majority of the plates. Geophysical implications of this numerical exercise are discussed.