中国西北及周边地区上地幔密度异常驱动小尺度对流

假设地震层析成像提供的地震波速异常对应于上地幔物质的密度异常分布,而该密度异常直接源于上地幔热对流相应的温度扰动. 在给定边界条件下,利用三维傅里叶变换,在波数域内求解控制流体行为的运动方程和连续性方程,得到上地幔小尺度对流流场. 利用密度异常驱动上地幔小尺度对流的数学 物理模型,采用胥颐、刘福田等提供的地震层析成像数据计算得到了我国西北及周边地区上地幔对流模式. 结果表明,对流流场的顶部在岩石圈较薄的盆地区域呈现上升发散流动特征,如塔里木盆地、柴达木盆地、哈萨克斯坦块体及准噶尔盆地;岩石圈较厚的山脉则对应了会聚下降的流动特征,如天山山脉、昆仑山山脉和祁连山山脉. 同时,塔里木盆地处于拉张状态,驱动其上地幔物质南下向青藏高原北部西昆仑运动,以及北上向天山下部流动,这可能是天山隆升的原因之一.      

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…     

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

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

南海西南海盆构造演化的热模拟研究

南海西南海盆的张裂和海底扩张是白垩纪末至中始新世南海形成过程中最重要的构造事件.本文采用三维有限单元法对该区的热演化过程进行了模拟计算.通过对变形、温度结构的计算,研究了西南海盆张裂变形、海底扩张持续时间、地幔物质上升、地壳岩墙沿扩张中心的挤入扩张活力、岩浆活动等.计算结果表明：由于其深部动力学条件不足,海盆一次扩张持续时间在10～15Ma之间,其后地幔物质的上升活动逐渐停止,地壳失去扩张动力,使得扩张中心成为残留扩张中心的死亡裂谷,而未构成中脊或中隆带.虽然该处地幔物质上升的潜力不足,但伴随局部的断裂,尤其是盆、缘边界的拆离拉张,仍能产生相当强烈的岩浆喷溢活动,导致此区海盆成型之后的海山崛起.     

上地幔变黏度小尺度对流的数值研究

基于二维模型，利用有限元方法，研究上地幔-岩石圈系统的变黏度小尺度对流. 考虑该系统的黏度随温度以指数形式变化，数值结果表明，当黏度随温度变化较剧烈时，由于低温高黏度，系统的最上部物质不参与对流，系统发育形成一个类似于岩石圈的静止盖层. 计算表面热流、地形起伏及重力异常与对流格局有较好的相关性，高热流、上升地形对应于对流的上升区，反之低热流、下降地形与对流的下降区对应.     

裂谷盆地构造-热演化模拟中几个问题的讨论

裂谷盆地的构造-热演化模拟是在岩石圈尺度计算裂谷盆地形成演化过程中的热历史和沉降史.拉张模型实现了构造和热的完美结合,在描述裂谷盆地沉降和热流演化方面取得了很大的成功.本文使用二维运动学模型,通过有限元方法,在拉格朗日坐标系下进行拉张背景下的构造热演化模拟,探讨了拉张模型中初始地壳、岩石圈厚度、软流圈对流、模型上边界对构造热演化的影响,以及载水和载沉积物两种情况下盆地侧翼抬升的差异.     

Rheology of the mantle and tectonics of the oceanic lithospheric plates

The modern concepts of the rheology of viscous mantle and brittle lithosphere, as well as the results of the numerical experiments on the processes in a heated layer with a viscosity dependent on pressure, temperature, and shear stress, are reviewed. These dependences are inferred from the laboratory studies of olivine and measurements of postglacial rebound (glacial isostatic adjustment) and geoid anomalies. The numerical solution of classical conservation equations for mass, heat, and momentum shows that thermal convection with a highly viscous rigid lithosphere develops in the layer with the parameters of the mantle with the considered rheology under a temperature difference of 3500 K, without any special additional conditions due to the self-organization of the material. If the viscosity parameters of the lithosphere correspond to dry olivine, the lithosphere remains monolithic (unbroken). At a lower strength (probably due to the effects of water), the lithosphere splits into a set of separate rigid plates divided by the ridges and subduction zones. The plates submerge into the mantle, and their material is involved in the convective circulation. The results of the numerical experiment may serve as direct empirical evidence to validate the basic concepts of the theory of plate tectonics; these experiments also reveal some new features of the mantle convection. The probable structure of the flows in the upper and lower mantle (including the asthenosphere), which shows the primary role of the lithospheric plates, is demonstrated for the first time.     

Influence of low viscosity asthenosphere on mantle flows

Numerical simulation in recent years has revealed that the cold lithosphere, whose viscosity is three to four orders of magnitude higher than that of the underlying mantle, behaves during mantle convection as a stagnant lid. On the basis of model calculations, this paper shows how convection changes over to this regime with increasing viscosity. Spatially fixed high viscosity inclusions and those moving with the convective flow have fundamentally different effects on the structure of convective flows. The model calculations indicate that anomalously low viscosity asthenospheric regions also lead to a specific regime of convection. With a decrease in the viscosity by more than three orders of magnitude, a further reduction in the viscosity of the region ceases to influence the structure of convection in the outer region. The boundary with this region behaves as a freely permeable boundary. In the low viscosity asthenospheric region itself, autonomous convection can arise under certain conditions.     

The structure of thermal convection and geotherms beneath a continent and ocean

The formation of the thermal cross section of the lithosphere and mantle upon the interaction between the mantle convection and the immobile continent surrounded by the oceanic lithosphere is studied by numerical modeling. The convective temperature and velocity fields and then the averaged geotherms for subcontinental and suboceanic regions up to the boundary with the core are calculated from the solution of convection equations with a jump in viscosity in the continental zone. Using the experimental data on the solidus temperature in the rocks of the upper mantle, the average thickness of the continental and oceanic lithosphere is estimated at 190 and 30 km, respectively. The effect of a hot spot formed in the subcontinental upper mantle at a depth of 250–500 km, which has not been previously noted, is revealed. Although the temperature in this zone is typically assumed to be close to adiabatic, the calculations show that it is actually higher than adiabatic by up to 200°C. The physical mechanism responsible for this effect is associated with the accumulation of convective heat beneath the thermally insulating layer of the continental lithosphere. The revealed anomalies can be important in studying the phase and mineral transformations at the base of the lithosphere and in the regional geodynamical reconstructions.     

中国大陆及邻区岩石圈三维流变结构

依据地震波速得到的上地幔温度和气象台站记录的地表温度为约束,结合地表热流和热导率观测数据,利用有限元方法计算了中国大陆及邻区岩石圈三维热结构.基于此温度结果和GPS观测得到的应变率数据,以滑动摩擦、脆性破裂和蠕变三种强度机制为约束,计算得到了中国大陆及邻区岩石圈三维流变结构.结果显示:弱强度和低等效黏滞性系数的下地壳在中国大陆及邻区普遍存在,并且下地壳的流变强度和等效黏滞性系数比上地壳和岩石圈地幔一般要低1~2个数量级;中国大陆范围内青藏高原存在着厚度最大、强度最低的下地壳;青藏高原的岩石圈强度和等效黏滞性系数比华北、华南和印度板块的都要低;岩石圈流变结构的横向分布特征与重力梯度带和地形过渡带比较一致.     

Neotectonic Crustal Uplift on the Continents and its Possible Mechanisms. The Case of Southern Africa

According to a large volume of data an intensive crustal uplift began in the Oligocene over most of the continental areas after a long period of relative tectonic stability. This Neotectonic uplift formed most of the present positive topographic features on the continents, and its strong acceleration took place during the last several million years. In many regions the uplift was associated with magmatism. The main methods of studying the Neotectonic uplift are considered together with the data on the uplift of Southern Africa. In this area the uplift took place in the Early Miocene (up to 300 m) and in the Late Pliocene and Pleistocene (up to 900 m). It occurred without stretching or shortening of the crust. Rapid erosion of the lower part of mantle lithosphere by a plume material is proposed as a mechanism of the uplift. This material ascended from below and rapidly spread along the base of the lithosphere. Its spreading for 1000 km during a few million years is possible only under a low viscosity of normal asthenosphere (1019 Pa s) and a much lower viscosity of a plume material (2 × 1016 Pa s). As in Southern Africa, in most of the regions the Neotectonic uplift was associated with insignificant shortening or stretching of the crust. This indicates that in some regions a plume material ascended from below and rapidly spread along the base of the lithosphere and eroded the mantle lithosphere in vast areas beneath the continents. In regions with a hot asthenosphere a strong weakening of the mantle lithosphere which allows its erosion can be associated with a high temperature of the plume material. In regions where the asthenosphere is at moderate temperature weakening of the mantle lithosphere can result from infiltration of volatiles from the plume material.     

青藏高原-天山地区岩石层构造运动的地幔动力学机制

利用全球重力大地水准面异常、板块绝对运动及全球地震层析成像数据,计算了青藏高原-天山地区岩石层下部地幔大尺度对流格局以及此种尺度对流驱动下岩石层内应力场分布;同时,利用区域均衡重力异常数据反演青藏高原中、北部到天山地区上地幔小尺度对流模型．结果表明,大尺度的地幔物质运移过程可能驱动着中国大陆岩石层整体从西部以南北方向为主的运动转向东部地区以北东和南东方向的运动;而该区域上地幔小尺度上升流动支持了现代青藏高原和天山地区的抬升运动．提出和讨论了青藏高原隆升的“断离隆升-挤压隆升-对流隆升”三阶段模式,并探讨了大陆岩石层构造运动的地幔深部动力学背景．     

Geophysical observations pertaining to solid-state convection in the terrestrial planets

We present a broad-based review of the observational evidence that pertains to or otherwise implies solid-state convection to be occurring (or have occurred) in the interiors of the terrestrial planets.For the Earth, the motion of the plates is prima facie evidence of large-scale mantle convection. Provided we understand upper-mantle thermal conductivity correctly, heat flow beneath the old ocean basins may be too high to be transported conductively from the upper mantle through the base of the lithosphere and therefore convection on a second smaller scale might be operative. The horizontal scale of plate dimensions implies, due to typical cell aspect ratios observed in convection, that the motion extends to the core-mantle boundary. Improved global data coverage and viscoelastic modeling of isostatic rebound due to Pleistocene deglaciation imply a uniform mantle viscosity, and thus indicate that whole-mantle convection could exist. Additionally, there is some seismic evidence of lithospheric penetration to depths deeper than 700 km. We discuss some salient features and assumption boundedness of arguments for convection confined to the upper mantle and for convection which acts throughout the mantle since the vertical length scale has a profound effect upon the relevance of geophysical observations. The horizontal form of mantle convection may be fully three-dimensional with complex planform and, therefore, searching for correlative gravity patterns in the ocean basins may not be useful without additional geophysical constraints. Many long-wavelength gravity anomalies may arise from beneath the lithosphere and must be supported dynamically, although thermal convection is not a unique explanation. Topography is an additional geophysical constraint, but for wavelengths greater than a few hundred kilometers, a general lack of correlation exists between oceanic residual gravity and topography, except at specific locations such as Hawaii. Theoretical calculations predict a complex relationship between these two observational types. Oceanic gravity data alone shows no regular planform and there is no correlation with any small-scale convective pattern predicted by laboratory experiments.All of the observational evidence argues against Martian plate tectonics occurring now or over much of the history of this planet, but lack of plate tectonics is not an argument against interior convection. The Tharsis uplift on Mars may have resulted from convective processes in the mantle, and the present-day gravity anomaly associated with Tharsis must be supported by the finite strength of the lithosphere or by mantle convection. Stresses imparted by the present topographic load would be greater than a kilobar, in excess of long-term finite strength. Observed fracture patterns are probably a direct result of this load, and the key question concerns the level of resultant strain relief. The global topographic and geomorphic dichotomy between the northern and southern hemisphere required a solid-state flow process to create the accompanying center-of-figure to center-of-mass offset.Lunar heat flow values, in analogy with oceanic heat flow on the Earth, strongly imply a convective mechanism of heat transport in the interior which, based on seismic Q values, is limited to the lower mantle. The presence of moonquakes in this region does not preclude solid-state convective processes. Lunar conductivity profiles provide no information on convection because of the difficulty in conductivity modeling, uniqueness of models, and the uncertainty in the conductivity-temperature relationship. The excess oblateness of the lunar figure over the hydrostatic value does not require convective support; in fact, such a mechanism is unlikely.The presence of a dipole magnetic field on Mercury does not provide a constraint on mantle convection unless its existence can be inextricably linked to a molten core. The non-hydrostatic shape of the equatorial figure, required for the observed $\text{3}\text{2}$ resonance between Mercury's rotational and orbital periods, is most likely related to surface processes, as opposed to convection. The $\text{3n}\text{2}$ resonance implies escape from a 2n resonance and, therefore, is related to the question of a molten core. Further dynamical data is needed to constrain interior models.Interpretation of limited radar imagery for the surface of Venus is enigmatic in terms of plate tectonics and therefore interior convection. Linear tensional and possibly compressional features are observed, but there are also crustal regions which appear to show large impact structures and are thus geologically old and may not have been recycled.     

Lithosphere types in North China: Evidence from geology and geophysics

Deep-seated materials from lithosphere are the ba- sic parameters and the foundation for geodynamic and continental dynamic studies. Division of lithosphere types and their deep-seated materials and structure can provide important evidence in interpreting the com- plex phenomena derived from the processes of forma- tion and evolution of continents, in evaluating the mineral resource potential, in predicting geological disasters and in the research of the continental dy- namic process. Huge lit…     

Convection model of post-glacial uplift

A new approach to analytical and numerical study of the process of the post-glacial uplifting of the Earth’s surface was proposed within the framework of a viscous model. Displacement of the Earth’s surface is considered as the motion of the density boundary due to chemico-density convection. It is shown that the incorporation of the non-Newtonian rheology at observed velocities of post-glacial uplifts requires an obligatory presence of faults in the lithosphere and gives rise to quasi-uniform motion of the mantle material, whose viscosity under the lithosphere is, on the average, sufficiently small and amounts to ～10¹⁹ Pa. The study of the stability of the constructed model of the post-glacial uplift considered as the chemico-density convection relative to the thermal convection shows that the velocity of thermal convection developing in the presence of a quasiuniform mantle flow related to the post-glacial recovery is ～1 m/yr.     

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

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

辽河裂谷的深部构造与裂谷活动的侧向迁移

对横穿辽河裂谷的闾阳－海城－东沟剖面的深部地震探测、重力测量、大地电磁测深、大地热流测量资料做了进一步分析,发现辽河裂谷地带壳幔间具有多个薄层过渡的特征,与东西两侧的辽东台隆和燕山台褶带明显不同。壳幔过渡带的薄层结构、速度结构、温度结构以及壳内低速层和软流圈顶面构造特点表明辽河裂谷的构造演化过程自西向东侧向迁移,当前上地幔物质运动的前锋不在裂谷之下,而是向东移至金州地震带附近的析木一带,深部物质运动是辽河裂谷形成及其活动性侧向迁移的主要动力学因素     

Influence of mantle convection to the crustal movement pattern in the northeastern margin of the Tibetan Plateau based on numerical simulation

Based on the recent observations about the movement and rheological structure of the lithosphere and deformation pattern of the crust, we developed a three-dimensional finite element model for the northeastern margin of the Tibetan Plateau. The model considered the impacts of both external and internal conditions, including mantle convection, gravitational potential energy and block interactions. We compared the simulated surface movement rates to the observed GPS velocities, and the results revealed that crustal movement gradually decreased toward the edge of the plateau. The factors controlling this pattern are the interactions of adjacent blocks, gravitational potential energy of the plateau, and also mantle convection as well. Additionally, according to the observation that there was an apparent difference between the horizontal movement rate of the lithosphere and convective velocity of the underlying mantle, and also based on the results of seismic anisotropy studies that suggest different strengths and deformation regimes of the lithosphere in different tectonic blocks, we proposed that the impact of mantle convection on the lithosphere may have varied in space, and introduced a parameter named mantle convection intensity factor in numerical simulations. Our simulation results show consistent surface movement rates with GPS observations, which further supports the viewpoint of seismic anisotropy studies, i.e., the degree of coupling between the crust and mantle varies significantly among different blocks.     

渭河盆地岩石圈热结构与地热田热源机理

岩石圈热结构是盆地现今地温场研究的重要延伸和扩展,是了解大陆岩石圈构造变形及演化等大陆动力学问题的重要窗口,更是地热田热源机理研究的核心问题.本次工作,在系统分析渭河盆地现今地温场和水动力系统基础上,编制了渭河盆地大地热流分布等值线图;通过实测生热率等热物性参数,利用一维稳态热传导方程计算了研究区岩石圈热结构,并分析了渭河盆地岩石圈热结构特征和地热田热源机理.结果表明,渭河盆地现今大地热流值分布范围为62.5~80.2mW·m-2,平均为70.8±4.8mW·m-2,西部明显高于东部,西安坳陷最高,咸礼凸起次之;渭河断裂并不是控热断裂,其沟通作用引起的水热循环一定程度上影响了浅部热量再分配,对渭河盆地地温场并没有起到明显的控制作用.西安坳陷—咸礼凸起地壳热流介于32.2~37.5mW·m-2之间,平均为34.6mW·m-2;地幔热流分布范围为33.8~38.9mW·m-2,平均为36.0mW·m-2;地壳热流和地幔热流的总体变化趋势一致,西安坳陷高于咸礼凸起,分析认为西安坳陷沉积层厚度大于后者,且沉积层放射性生热率更大,是造成西安坳陷地壳热流高于咸礼凸起的原因,而西安坳陷相比咸礼凸起更高的地幔热流,表明西安坳陷深部活动性强于咸礼凸起.西安坳陷和咸礼凸起地壳/地幔热流比值相近,介于0.93~1.01之间,平均为0.96,"热"岩石圈厚度约为95~101km.渭河盆地岩石圈热结构特征与鄂尔多斯盆地在很大程度上具有相似性,暗示着二者具备相似的深部稳定性,这与渤海湾盆地为代表的中国东部中—新生代主动裂谷盆地岩石圈热结构特征截然不同,表明渭河盆地为被动伸展裂陷.从鄂尔多斯盆地、渭河盆地、山西裂谷到华北盆地,"热"岩石圈厚度的有序变化表明太平洋板块俯冲引起的地幔对流对华北地块深部动力学行为的影响主要发生在太行山以东,而太行山以西的鄂尔多斯盆地和渭河盆地则影响甚微,这种空间差异影响从侧面暗示着华北克拉通破坏过程的有序性.综合分析渭河盆地地质—地球物理资料认为,岩石圈表层伸展破裂、深部重力均衡调整进而引起软流圈被动上涌,其产生的相对高地幔热流的热传导和深大断裂沟通的水体热对流相互叠加作用,共同构成了渭河盆地中—低温地热田的热源机理.