第10届地幔对流模拟和岩石圈动力学国际讨论会概述

第10届地幔对流模拟和岩石圈动力学国际讨论会于2007年9月2日至9月7日在法国Cayy-le-Rouet市举行,来自13个国家的近百名科学家参加了这次讨论会.     

大陆岩石圈、地幔底部异常体与地幔对流相互作用的数值模拟

在地球表层存在着占地表面积约30%的具有低固有密度、高黏度的大陆岩石圈.由于其特殊的物理化学性质,大陆岩石圈通常不直接参与下方的地幔对流,但其与地幔对流格局有着重要的相互影响.大量研究显示,在中太平洋和非洲的下地幔底部,存在着两块占核幔边界（CMB）面积约20%的高密度热化学异常体（由于其剪切波速度较低,常称作低剪切波速度省（LSVPs））.LSVPs的演化既受地幔对流的影响,同时也影响地幔物质运动的格局和动力学过程.本文系统研究了存在大陆岩石圈,下地幔LSVPs的地幔对流模型.模拟结果显示：（1）当大陆体积较小时,其边缘常伴随着俯冲,大陆区域地幔常处于下涌状态,其上地幔温度较低,大陆岩石圈在水平方向处于压应力状态.随着大陆体积的增大,大陆边缘的俯冲逐渐减弱,大陆区域地幔由下涌转为上涌,其上地幔温度较高,大陆岩石圈水平方向处于拉应力状态.（2） 岩石圈与软流圈边界（LAB）在大陆下方较深,温度较低;在海洋区域较浅,温度较高.随着大陆体积的增大,陆洋之间LAB深度、温度的差异逐渐减小.（3）大陆区域地幔底部LSVPs物质的丰度与大陆的体积呈正相关.当大陆体积较小时,大陆下方的LSVPs丰度比海洋区域少.随着大陆体积的增大,大陆下方LSVPs的丰度逐渐增大.（4）海洋地区地表热流高,且随时间波动大,大陆地区地表热流低,随时间波动较小;LSVPs区域的核幔边界热流低.     

地幔对流对全球岩石圈应力产生与分布的作用

利用动力学模拟方法研究地幔对流对于大尺度岩石圈内部应力场形成的作用. 地幔物质内部的密度横向非均匀及表面板块运动引起地幔流动,并在岩石圈底部产生一个应力场. 该应力场作为面力将造成岩石圈本身变形,从而产生岩石圈内部的应力分布. 模拟计算结果表明,大部分俯冲带及大陆碰撞带区域应力均呈现挤压特征,如环太平洋俯冲带及印度－欧亚碰撞带等;而东太平洋洋脊、大西洋洋脊及东非裂谷处应力状态均表现为拉张;并且绝大多数热点位置处于应力拉张区域,这与目前对全球构造应力状态的理解是一致的. 计算的岩石圈内部最大水平主压应力的方向与观测表现出相当的一致,其结果总体上吻合得较好,然而在局部区域（例如西北太平洋的俯冲带、青藏高原等地区）存在着较大的差异. 研究表明,地幔对流是造成岩石圈内部大尺度应力状态及分布的一个重要因素.     

中国东部地幔矿物中富硅、碱熔融包裹体：对岩石圈演化的启示

中国东部华北、东北和东南沿海几个地区的新生代玄武岩的尖晶石相和石榴石相橄榄岩捕虏体矿物中熔融包裹体研究表明, 玻璃化学成分为富硅(SiO₂=60%~68%)、碱(K₂O+Na₂O= 5%~11%), 特别是富钾(K₂O>Na₂O), 以及H₂O和CO₂等挥发分(2%~7%)的中酸性硅酸盐, 属于高钾钙碱性系列的英安岩和安山岩, 少量属于橄榄安粗岩系列. 熔融包裹体中高Al和Ca透辉石子矿物是熔体在高温高压下结晶的产物, 而非脱玻化重结晶的产物. 研究指出, 这种富钾(一般K₂O>3%)中酸性硅酸盐更具大陆特点, 与主矿物和寄主玄武岩均无成因关系, 是被捕获的中生代岩石圈壳-幔相互作用之熔体, 同时暗示中新生代时期中国东部大陆岩石圈地幔可能经历了破碎和更新的过程. 这无疑给中国东部大陆岩石圈演化过程提供重要启示, 同时指出富硅、碱熔体应对中国东部地幔化学非均一性负责.     

地幔对流的实验研究:非立柱状幔柱和地幔涡旋

地幔对流的物理模拟实验结果表明 ,在地幔介质和温度非均匀分布的复杂条件下 ,热卷流 (地幔柱 )往往由立柱状转变为非立柱状 (含斜柱状、涡旋状等 )。在忽略科里奥利力的情况下 ,板块的下插和滞积下沉、岩石圈根的存在以及地幔介质粘度的非均匀分布等都可能构成不同形状的障碍 -导流体 ,导致地幔的涡旋运动。软流圈中的水平涡旋环带属于对数螺线型 ,环带旋转半径及线速度逐渐减小 ,最终在旋转中心处下沉 ,而旋转角速度大致保持恒定     

地幔对流的数值实验研究

数值实验是地幔对流研究的基本方法。本文概述了近二十年地幔对流研究的主要进展,介绍并评价了每个研究阶段所提出的新概念和新思想。最后,着重讨论了上地幔与全地幔对流的论争。     

金星非单调冷却热演化历史分析

金星表层年龄和构造活动特点表明其岩石层在最近的地质历史时期经历过广泛的更新.这种全球性的表层改造与其内部热演化历史进程密切相关.如果金星存在相变形成的上、下地幔,依据现今所了解的金星物理性质和参量化的热对流理论,并且考虑金星地幔相变边界层状态对对流的控制作用,我们计算了金星热演化历史.结果表明,金星的热演化历史是一种非单调的冷却过程,在这种非单调的热演化历史进程中,金星地幔会出现大体等周期的翻转.由于参数选取的不同,翻转时金星上、下地幔的温差随时间可能出现稳定变化、逐渐加强、逐渐减弱三种不同演化模式,目前尚不能确定实际金星热演化历史究竟是哪一种模式.金星地幔相变边界层的穿透对流可能是推动其表层岩石层全球性更新的关键,导致其表层火山活动和地表构造以大致500 Ma时间间隔更新和重造.     

地幔对流的数值实验研究

数值实验是地幔对流研究的基本方法。本文概述了近二十年地幔对流研究的主要进展,介绍并评价了每个研究阶段所提出的新概念和新思想。最后,着重讨论了上地幔与全地幔对流的论争。     

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

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

大陆岩石圈地幔形成与演化研究的新进展

研究大陆岩石圈地幔形成与演化的主要途径是利用岩石－地球化学探针研究幔源橄榄岩的年龄、组成及其演化，目前多将岩石圈地幔划分为大洋型和克拉通型，并认识到，岩石圈地幔的矿和相、成分、年龄在垂向上具有分层性；根据研究资料，目前已提出了以熔体抽取为岩石圈地幔形成的主要方式，并通过交代、拆沉等过程使岩石圈地幔遭受改造，中国东部岩石圈地幔的成分主体上是饱满的，与其上的太古代地壳并不协调，一些研究者提出了岩石圈拆沉和地幔置换两种具有代有性的模型来解释岩石圈地幔成分的这种特点，人们对中国东部岩石圈地幔作了大量的研究工作，然而，若要形成较为合理的解释，还需应用新技术（如，Re/Os同位素技术）做进一步的工作。     

The gravity field and crustal thickness of Venus

The gravity and topography of Venus obtained from observations of the Magellan mission, as well as the gravity and topography from our numerical mantle convection model, are discussed in this paper. We used the hypothesis that the geoid of degrees 2–40 is produced by sublithospheric mantle density anomalies that are associated with dynamical process within the mantle. We obtained the model dynamical admittance(the geoid topography ratio based on a convection model) by a numerical simulation of the Venusian mantle convection, and used it to correct the dynamical effect in the calculation of crustal thickness. After deducting the dynamical effect, the thickness of the Venusian crust is presented. The results show that the gravity and topography are strongly correlated with the Venusian mantle convection and the Venusian crust has a significant influence on the topography. The Venusian crustal thickness varies from 28 to 70 km. Ishtar Terra, and Ovda Regio and Thetis Regio in western Aphrodite Terra have the highest crustal thickness(larger than 50 km). The high topography of these areas is thought to be supported by crustal compensation and our results are consistent with the hypothesis that these areas are remnants of ancient continents. The crustal thickness in the Beta, Themis, Dione, Eistla, Bell, and Lada regiones is thin and shows less correlation with the topography, especially in the Atla and Imdr regiones in the eastern part of Aphrodite Terra. This is consistent with the hypothesis that these highlands are mainly supported by mantle plumes. Compared with the crustal thickness calculated with the dynamical effect, our results are more consistent with the crust evolution and internal dynamical process of Venus.     

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.     

Thinning and destruction of the cratonic lithosphere:A global perspective

It has been proposed that the North China Craton(NCC)was thinned up to a thickness of100 km during the Phanerozoic,and underwent an associated craton destruction.Evidently,it is an important topic worthy of future study to understanding the mechanism of cratonic destruction and its role played in the continental evolution.After synthesized the global cratons of India,Brazil,South Africa,Siberia,East Europe(Baltic)and North America,we found that lithospheric thinning is common in the cratonic evolution,but it is not always associated with craton destruction.Most cratons was thinned by thermal erosion of mantle plume or mantle upwelling,which,however,may not cause craton destruction.Based on the studies of the North American and North China Cratons,we suggest that oceanic subduction plays an important role in caton destruction.Fluids or melts released by dehydration of the subducted slabs metasomatize the mantle wedge above and trigger extensive partial melting.More importantly,the metasomatized mantle lost its original rigidity and make craton easier to be deformed and then to be destoyed.Therefore,we suggest that the widespread crust-derived granite and large-scale ductile deformation within the continental crust can be regarded as the petrological and structural indicators of craton destruction,respectively.     

大陆岩石圈流变学研究进展

大陆岩石圈流变学是固体地球科学领域的基本问题,也是理解大陆变形的动力学过程和构造演化史的关键.本文对近年来大陆岩石圈流变学的研究进展进行了综述.重点讨论大陆岩石圈流变学的纵向分层和横向分块特性及流变学与地震活动性、壳－幔解耦及下地壳流动变形的地球动力学意义和中国大陆地区岩石圈流变学研究状况.最后,对大陆岩石圈流变学研究存在的问题和未来的研究方向做了一定的探讨和展望.     

Upper mantle flow beneath the Northwest of China and its lithospheric dynamics

ＵｐｐｅｒｍａｎｔｌｅｆｌｏｗｂｅｎｅａｔｈｔｈｅＮｏｒｔｈｗｅｓｔｏｆＣｈｉｎａａｎｄｉｔｓｌｉｔｈｏｓｐｈｅｒｉｃｄｙｎａｍｉｃｓＪＩＡＮ－ＨＵＡＨＵＡＮＧＩ（黄建华）；ＸＩＡ－ＨＵＡＣＨＡＮＧＩ（常筱华）ａｎｄＲＯＮＧ－ＳＨＡＮＦＵＩ傅容珊）（...     

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.