欧亚板块均衡大地水准面异常及其动力学特征

采用弹性挠曲模型,利用球谐函数方法计算了欧亚板块均衡大地水准面异常.均衡大地水准面异常消除了地形起伏、地壳物质补偿的影响,主要反映了地幔物质密度异常以及核幔边界起伏和热状态的不均匀分布特征,而这种物质密度和热状态的不均匀分布导致了地幔物质的流动.针对滤波后的均衡大地水准面异常分别计算了欧亚板块大尺度和中尺度的地幔对流,并进行了构造和地球动力学分析.结果表明,现今欧亚板块的地壳运动与地幔对流有着密切的联系.     

中国及邻区大地水准面异常的场源深度探讨

利用中国及邻区地形、地震层析成像、沉积层底面、Moho面及岩石层底面资料,正演计算出中国及邻区岩石圈大地水准面异常;再从全阶大地水准面异常中扣除正演模拟得到的岩石圈大地水准面异常与不同阶次波段的大地水准面进行比较,寻求表示中国及邻区地幔物质不均匀的大地水准面异常的最佳阶次为2-60阶. 结果表明,对应于岩石圈的大地水准面异常的重力位球谐函数阶数为61-20阶;下地幔重力位球谐函数阶数为2-6阶;而-60阶重力位球谐函数则表示中国及邻近区域上地幔大地水准面异常.     

黏滞分层地幔中密度异常驱动对流模型的研究

在地震层析成像计算的地幔密度异常直接驱动地幔对流的新方法的基础上,发展了在上、下地幔不同黏性结构框架下,密度异常驱动地幔对流的物理模型.利用 Grands和S12 WM13等地震层析成像模型推得的地幔密度异常分布,设置板块绝对运动极型场为运动上边界,考虑深度660km地震波不连续面为界的上、下地幔之间存在黏滞性的差异,直接反演了不同黏滞系数的双层地幔结构下地幔对流的模式.研究中选取地幔平均密度为ρ=5500kg/m3, 上层地幔平均黏滞系数为μ=1021Pa·s,计算了上、下地幔黏滞系数之比为1∶1, 1∶10, 1∶100和1∶1000时地幔大圆剖面、以及区域剖面上的流场.结果表明,两种模型在球谐展开1～13阶的范围内其对流的基本格局相似.当下地幔黏滞性超过上地幔的100倍时,下地幔流场速度与上地幔的流场速度相比显著减小,但是对流仍然表现出单层对流环的基本格局.论文还用 240km深度球面上的对流格局讨论了对流和全球构造之间的关系.     

利用多种地球物理观测资料直接反演地幔对流模型

假定地幔为一个均匀的、粘滞系数为常数、同时均匀分布放射性热源的流体球层,其内部存在的对流则由流体力学3个基本方程:运动方程、能量方程和连续性方程确定.如果假定地幔处于低瑞利数的状态（临界瑞利数1.5倍左右）,那么上述方程中的非线性项可以忽略不计.作为一类可能的模型,本文计算一组用6个边界条件确定6个未知数的线性方程组.这些条件包括板块绝对运动极型场、地球大地水准面异常和地震层析结果提供的地幔密度分布横向不均匀相应的“刚性地球”水准面异常等.模型计算表明:1.地幔中流体运动格局不仅受地幔热动力学参数（瑞利数）控制,而且强烈地受边界条件的影响.2.若不限定下边界为等温边界,则上、下地幔之间并不呈现出活动性明显差异;但是在模型瑞利数加大到一定值时,核-幔边界附近将出现一些局部的小尺度对流环.3.当模型瑞利数从很小增加时,对流格局将发生变化,这些格局可能反应由地幔热动力学参数决定的地幔固有特性.4.当瑞利数为50000和80000时,核-幔边界形变与PcP波得到的结果吻合较好.     

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

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

青藏地区大地水准面形态及其与构造动力学的关系

利用密合法求解青藏地区大地水准面模型，研究青藏地区大地水准面主要特征及其与该区内部构造和活动的相关性．表明青藏地区大地水准而异常与青藏高原隆起的形状极为相似，大地水准面整体负异常与下地幔相关，形状起伏与上地幔相对应，岩石层反映的大地水准而异常与青藏高原活动构造有较好的对应关系．     

地幔对流的数值模拟及其与表面观测的关系

本文从基本的热对流方程出发,并结合地幔对流特点,特别考虑到自重及非线性影响,探讨地幔对流及其与表面观测的关系,发展了相应的数值方法.结果表明,计算得到的长波大地水准面、地表地形、板块速度场水平散度与观测值符合程度较好.上、下地幔的非绝热温度异常与由地震层析得到的地震波速异常显示一定的相关性.地幔内部的流动呈现复杂形态,反映了高瑞利数对流的特征.     

弹性岩石层,区域重力异常和上地幔小尺度对流

在小尺度地幔对流模型的基础上,建立了顾及岩石层与地幔耦合的地幔对流模型,考虑到岩石层对上地幔小尺度对流的弹性响应,将上地幔视为一均匀的等粘滞系数的牛顿粘滞流体,其对注能量来源于下地幔,用热流体动力学基本方程得到地上地幔对流在壳幔边界处的垂向应力作为弹性板弯曲方程的垂直加载,来耦合弹性岩石层和可流动地幔,推导了两者耦合下区域重力异常和上地幔小尺度对流的相关方程,用以反演上地幔小尺度对流场模式和岩石层底部拖力场格局,对比了有,无弹性岩石层影响的模型之间的差异,并对模型应用范围,特别要板块边界水平均造力可能产生的影响,作了进一步探讨。     

上地幔密度异常驱动小尺度对流及实验模型

建立了由密度异常驱动上地幔小尺度对流的数学 物理模型, 发展了利用地震层析成像数据反演上地幔小尺度对流的基本理论和方法. 该模型建立在三维直角坐标系框架上, 假设地震层析成像所显示的地震波速度异常对应于上地幔物质密度异常, 而该密度异常反映了上地幔小尺度热对流系统的温度异常场. 模型首先将地震层析成像确定的地震波速度异常转换为密度异常, 并视其为对流的驱动力; 进而利用三维傅立叶变换, 在波数域内, 在给定的边界条件下, 求解控制流体行为的运动方程和连续性方程, 最后求得对流的流场. 为检验本研究提出的理论和方法的有效性, 本文使用了两个简单的实验模型: 热体和冷体模型; 俯冲断离( break off)板片模型, 计算了其驱动的地幔流场. 结果表明, 本文提供的理论和方法, 可以直接应用于与区域岩石层构造动力学相关的上地幔小尺度对流的研究.      

蒙古-贝加尔地区上地幔小尺度对流及 地球动力学意义

上地幔小尺度对流是控制区域地球动力学过程的主要机制之一,蒙古-贝加尔地区的一些区域动力学过程被认为与上地幔小尺度对流相关.本文目的在于利用重力资料研究蒙古-贝加尔地区的上地幔小尺度对流,并探讨其与构造动力学的关系.基于区域均衡重力异常与上地幔小尺度对流的相关方程,本文利用区域均衡重力异常资料反演了蒙古-贝加尔地区上地幔小尺度对流流场及作用于岩石层底部的应力场.结果显示,蒙古-贝加尔地区地幔流场及对流应力场呈现非常复杂的图像,流场及应力场分布与地表构造具有很好的相关性.西伯利亚地台和蒙古褶皱带下地幔流场和对流应力场均较弱,这与这些地区现今较弱的构造活动性是一致的.贝加尔裂谷区下存在地幔上升流,对流应力场呈拉张状态,但应力场的幅值较小(约8 MPa),表明地幔对流不是贝加尔裂谷开裂的主要控制因素.Hangay高原、阿尔泰和戈壁-阿尔泰下存在地幔上升流,对流应力场为拉张状态,这一方面可能构成Hangay高原隆升的深部动力机制,另一方面,也为Amurian板块西边界划分提供了动力背景.     

Hotspot distribution, gravity, mantle tomography: evidence for plumes

Thermal convection is the motor of Earth dynamics and therefore is the link between plate motions, hotspots, seismic velocity variations in the mantle, and anomalies of the gravity field. Small scale mantle anomalies, such as plumes, do, however, generally escape detection by tomographic methods. It is attempted to approach the problem of detection in a somewhat statistical manner. Correlations are sought between spherical harmonic expansions of the fields under study: the hotspot distribution, mantle velocity variations, gravity, heat flow. Using spherical harmonic representations of global fields implies integration and averaging over the whole globe. Thus, although relationships may remain masked in the space domain by a multitude of effects, tendencies may become visible in the spectra or in appropriate averages.The main results are the following: There is a significant long wavelength (n=2,3) negative correlation between the hotspot density and the P-wave velocity variation in the lower mantle. Positive hotspot density of degree 2 to 9 generally correlates with low seismic velocity in all depths of the upper mantle and with positive gravity. This fits well with plume-type convection. These results are also confirmed regionally for a number of individual mid-ocean ridges and hotspots. The hotspot density and the free air anomalies are distinctly positive above regions of low velocity extending to great depth. The effect is not distinct at ridges with shallow velocity anomalies. In a general way, we suggest that the antipodal upwellings (Pacific, Africa) are divided by downwelling currents around the shrinking Pacific. Plate boundaries can easily move away from their past connections with the deeper mantle. Small scale plume currents seem to be depicted in the hotspot expansion. © 1999 Elsevier Science Ltd. All rights reserved.     

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.     

Geoid anomalies over Gorringe Ridge,North Atlantic Ocean

Gorringe Ridge is a strong uplifted block of oceanic crust and upper mantle lying at the eastern end of the Azores-Gibraltar plate boundary. The geoid over this structure derived from Seasat altimeter data exhibits a 9-m height anomaly with a north-south lateral extension smaller than 200 km. An attempt is made to interpret this geoid together with the gravity anomalies and with the seismicity, which has been compiled as a function of depth.It is first shown that the flexure of the oceanic lithosphere due to the ridge loading does not provide a good fit of the geoid anomalies and probably should be discarded, as it assumes a continuous unfractured elastic plate.Models involving local heterogeneities are then tested. The comparison of the observed geoid anomalies with the anomalies due to the uncompensated relief indicates that the topographic high has no shallow compensation.Uncompensated models, previously proposed to explain the gravity anomalies, are tested using the geoid. One model (Purdy and Bonnin, in Bonnin [11]), which involves an uplift of upper mantle material at depth, generates too strong geoid anomalies and must be discarded. Another model, which represents a nascent subduction zone (Le Pichon et al. [25]), fits both the gravity and geoid anomalies, but leads to difficulties in explaining the deep seismicity north of Gorringe Ridge.A model in isostatic equilibrium is also able to fit both gravity and geoid anomalies. This model involves a deep root of density 3.0 g cm^?3, as has been previously proposed for many oceanic ridges and plateaus. This model is compatible with the deep seismicity, but the origin of this low-density material at great depth is up to now an unresolved question.More likely, dynamical models taking into account the forces induced by the convection related to the slow plate convergence in this area will have to be considered.     

欧亚地区均衡残差大地水准面和上地幔强度

首先计算了欧亚地区均衡残差大地水准面．基于地幔热对流的内负荷理论和最新全球层析成像结果，探讨了欧亚地区中波长均衡残差大地水准面的地球动力学意义．研究结果表明，中波长均衡残差大地水准面主要受上地幔粘滞度和岩石层强度的影响，进而得出欧亚地区一些古老地盾和构造稳定地区的上地幔与年轻山脉及构造活动地区的上地幔结构存在着差异．这个差异主要是占老地盾和构造稳定地区，如波罗的海地盾、中西伯利亚地台、东欧等区域，冷却的上地幔已穿透地幔较深，上地幔与岩石层之间耦合较好;而年轻山脉和构造活动区，如帕米尔、天山、贝加尔活动带、青藏高原、日本海周围地区，在上地幔可能存在着热物质即粘滞度很低的软流层，上地幔与岩石层耦合程度较差，甚至有可能解耦．从欧亚地区上地幔属性的差异，可以解释该地区的一些地球动力学问题．     

Tomographic structure of East Asia:Ⅱ.Stagnant slab above 660 km discontinuity and its geodynamic implications

P-wave arrival times of both regional and teleseismic earthquakes were inverted to obtain mantle structures of East Asia.No fast(slab) velocity anomalies was not find beneath the 660-km discontinuity through tomography besides a stagnant slab within the transition zone.Slow P-wave velocity anomalies are present at depths of 100-250 km below the active volcanic arc and East Asia.The western end of the flat stagnant slab is about 1 500 km west to active trench and may also be correlated with prominent surface topographic break in eastern China.We suggested that active mantle convection might be operating within this horizontally expanded "mantle wedge" above both the active subducting slabs and the stagnant flat slabs beneath much of the North China plain.Both the widespread Cenozoic volcanism and associated extensional basins in East Asia could be the manifestation of this vigorous upper mantle convection.Cold or thermal anomalies associated with the stagnant slabs above the 660-km discontinuity have not only caused a broad depression of the boundary due to its negative Clapeyron slope but also effectively shielded the asthenosphere and continental lithosphere above from any possible influence of mantle plumes in the lower mantle.     

The modern geoid and ancient plate boundaries

The non-hydrostatic geoid is dominated by three large anomalies: an area of high gravity potential in the equatorial Pacific; another stretching from Greenland through Africa to the southwest Indian Ocean; and a semi-continuous low region passing from Hudson's Bay through Siberia to India and on to Antarctica. None of these three high-amplitude (greater than 60 m) and long-wavelength anomalies corresponds to present-day plate boundaries. However, if the modern geoid is plotted over the positions of continents and plate boundaries at 125 Ma B.P. (reconstructed relative to hotspots) a strong correlation emerges. The modern geoidal low corresponds in position to the areas of subduction surrounding the Pacific 125 Ma ago. The geoidal high now centered on Africa is entirely contained within ancient Pangaea, and the equatorial Pacific high overlies the location of the spreading centers preserved in the magnetic anomalies of the central Pacific. The most plausible cause of the large geoidal undulations is lower mantle convection only weakly coupled to plate motions. The correspondence between modern geoid and ancient plate boundaries implies either that the coupling was much more intimate in the past, or that there is a lag of at least 100 Ma in response of the lower mantle to upper mantle conditions.     

华北克拉通及东邻西太平洋活动大陆边缘地区的P波速度结构：对岩石圈减薄动力学过程的探讨

本文通过地震层析成像研究获得了华北克拉通及其东邻地区（30°N-50°N,95°E -145°E）1°×1°的P波速度扰动图像.结果显示,在西太平洋俯冲带地区,上地幔中西倾的板片状高速异常体与其上方的低速异常区构成俯冲带与上覆地幔楔的典型速度结构式样.俯冲板片高速体在约300~400 km深度范围内被低速物质充填,暗示俯冲板片可能发生了断离.在华北克拉通地区的上地幔中发现三个东倾排列的高速异常带.在此基础上,本文构建了华北克拉通及其东邻西太平洋活动大陆边缘地区的上地幔速度结构模式图,并据此探讨克拉通岩石圈减薄与西太平洋活动大陆边缘的深部动力学联系.本文认为,太平洋板片的俯冲（断离）,触发热地幔物质上涌并在上覆地幔楔中形成对流,使克拉通岩石圈受到改造（底侵与弱化）.随着俯冲板片后撤,地幔楔中的对流场以及对岩石圈改造的影响范围均随之东移,最终导致华北克拉通岩石圈自下而上、从西向东分三个阶段依次拆沉减薄.这一模式能很好地解释现今克拉通岩石圈自西向东呈台阶状减薄的深部现象.