地幔内变粘度的滞弹性流体运动的三维有限元模型

建立了描述具有无限大Ｐｒａｎｄｌｔ数，变粘度的滞弹性流体运动的稳态和动态三维有限元公式，其中运动方程的有限元公式是由用应力分量表示的平衡方程，再利用牛顿流体的本构关系导出的。这种处理方式，不仅能够分析流体粘度变化所起的动力学作用，而且还能够直接计算温度分布和由流体运动所引起的扰动压力。最后，论证了导出的有限元公式和程序的正确性和可靠性，并给出了具有三维流动特征的算例     

地幔对流

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

地幔粘度结构的研究

地幔的流变性质已成为认识地球内部结构及动力学过程的核心问题之一。本文总结了近年来地幔粘度结构研究的方法,其中包括利用地球物理观测资料进行反演计算和实验室试验研究,重点讨论了利用冰期后回跳,板块运动速度和大地水准面异常资料反演地幔粘度结构２的方法和结果以及地幔矿物的实验结果,并对不同方法进行了比较总结。最后简单讨论了地幔粘度结构研究存在的问题和未来的研究方向。     

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

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

地幔动力学研究进展--地幔对流

分类而且系统地回顾了地幔对流研究的进展情况.发现虽然静态地幔对流模型(没有考虑地球自转或差异旋转效应)在拟合某些地表观测如重力异常、大地水准面异常以及地下应力场时取得了较满意的结果,但是却对于某些全球尺度的地质特征如0°、180°半球的非对称性显得难以适应,因此建议性地提出了今后地幔对流的可能发展方向.     

地幔对流的数值实验研究

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

地幔对流的数值实验研究

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

地幔对流（讲座）

地幔对流（讲座）李荫亭（中国北京１０００８１国家地震局地球物理研究所）引言板块构造学说的出现与发展一直同地幔对流的研究密切联系在一起，地幔对流研究已成为当代地球科学中的一个基础性课题。板块构造学说的驱动力问题是目前地球科学家最为关心的问题之一。国际地...     

Convection in the mantle with an endothermic phase transition

An endothermic phase transition at a depth of 660 km in the mantle partially slows down mantle flows. Many models considering the possibility of temporary layering of flows with separation of convection in the upper and lower mantle have been constructed over the past two decades. The slowing-down effect of the endothermic phase transition is very sensitive to the slope of the phase-equilibrium curve. However, laboratory measurements contain considerable uncertainties admitting both a partial convection layering and only an insignificant slowing down of a part of downgoing mantle flows. In this work, we present results of calculations of mantle flows within a wide range of phase-transition parameter values, determine ranges of one-and two-layer convection, and derive dependences of the amplitude and period of oscillations on phase-transition parameters.     

华北地区地壳上地幔S波三维速度结构

利用华北地区大型流动地震台阵的记录资料,采用近震和远震联合成像方法,得到了水平分辨率0.5°×0.5°、深至600km的S波速度结构.研究结果表明,上地壳S波速度结构与地表地质构造基本一致,燕山—太行山山脉均呈现高速异常,延庆—怀来盆地、大同盆地表现为低速异常,华北盆地内部的拗陷和隆起分别呈现低速和高速.唐山地区中地壳、山西裂陷盆地中下地壳存在明显的低速异常,可能分别与流体和热物质作用有关,有利于形成孕育强震的地质构造环境.90km的速度结构图像依然与地表的构造特征有较大的相关性,可能说明深部结构对地表构造有一定的控制作用.燕山隆起区岩石圈的厚度可达120～150km左右,华北盆地的岩石圈厚度可能在80km左右,太行山地区的岩石圈厚度介于两者之间.山西裂陷盆地上地幔低速层较厚,反映了该区不稳定的构造环境造成了地幔热物质的上涌.华北盆地下方220～320km出现的高速异常体,可能揭示了华北盆地上地幔仍然存在拆沉后残留的难熔、高密度的古老岩石圈地幔.研究区东部地幔转换带呈低速异常,推测可能与太平洋板块俯冲至该区下方地幔转换带前缘120°E左右的俯冲板块相变脱水有关.     

鄂尔多斯地区上地幔岩石圈三维速度结构面波反演研究

双平面波拟合法是一种新的面波成像方法,反演中考虑地震波场中的非平面波成分,提高反演的分辨率.本文利用双平面波拟合法,反演获得鄂尔多斯地区上地幔岩石圈的速度结构.所用资料为国家数字地震台网69个宽频带地震仪和北京大学34个流动数字地震台观测到的地震波面波资料.首先从面波记录中提取了研究区域20~125 s瑞利波相速度频散曲线,进而得到各个周期瑞利波相速度异常分布图.结果显示,短周期瑞利波相速度异常与地表的构造特征吻合较好,中长周期的瑞利波相速度可以反映出上地幔岩石圈的速度异常分布以及构造特征.由研究区20~125 s的瑞利波相速度分布图可以反演得到地表到地下200 km范围内的三维剪切波速度结构.结果显示,鄂尔多斯块体内部稳定均一,活化或改造的痕迹不明显;鄂尔多斯块体西南缘受到青藏高原的强烈作用,有大量地幔物质流动的痕迹存在;中央转换带下超过200 km深度存在地幔物质上涌,可能与太平洋板块的俯冲和青藏高原板块的挤压有关.     

3-D velocity structure beneath the crust and upper mantle of Aegean Sea region

The region of the Aegean Sea and the surrounding areas in the Eastern Mediterranean lies on the boundary zone between the Eurasian and the African plates. It is a zone of widespread extensive deformation and, therefore, reveals a high level of seismicity.Three-dimensional velocity structure, beneath the crust and upper mantle of the region between 33.0°N–43.0°N and 18.0°E–30.6°E, is determined.The data used are arrival times ofP-waves from 166 earthquakes, recorded at 62 seismological stations. In total, 3973 residual data are inverted.The resultant structure reveals a remarkable contrast of velocity. In the top crustal layer, low velocities are dominant in Western Turkey and on the Greek mainland, while a high velocity zone is dominant in the Ionian Sea and in the southern Aegean Sea.In the upper mantle, high velocity zones dominate along the Hellenic arc, corresponding to the subducting African plate and in the northern part of the region, corresponding to the subducting African plate and in the northern part of the region, corresponding to the margin of Eurasian plate.A low velocity zone is dominant in the Aegean Sea region, where large-scale extension and volcanic activity are predominant, associated with the subduction of the African plate.     

Simultaneous inversion for 3-D variations in shear and bulk velocity in the mantle

Until recently, most of the seismic tomographic modeling has been addressing the question of lateral heterogeneity either in P- or S-wave velocities. The S-wave velocity variations are larger and hence provide stronger signal on long-period waveforms. The direct P travel times, being the first arrivals, on the other hand, are most frequently reported in the International Seismological Centre (ISC) Bulletins. In mineral physics experiments, the variation in bulk velocity is more often measured. To better understand the differences between δv_P and δv_P patterns and better link the results from mineral physics to those of seismic tomography, we formulate the inverse problem in terms of relative perturbations in the shear velocity v_S = (μ/)^1/2 and bulk sound velocity v_Φ = (K/)^1/2. We use a large data set which consists of waveforms, waveform-derived travel times and travel times from the ISC Bulletins. The earthquakes are relocated using corrections for lateral heterogeneity. The events which cannot be reliably determined are discarded. The model is defined as spherical harmonics to Degree 12 horizontally and as Chebyshev polynomials to order 13 radially, for both shear and bulk sound velocity. The inversion is performed under smoothness constraint. The resolution tests and bootstrapping analysis indicate that the model is well recovered, particularly at long wavelength.

The results indicate a much larger variability of shear than bulk sound velocity. The model explains observations well. The most intriguing result obtained in this study is that the variations in shear velocity and bulk sound velocity are negatively correlated in the lowermost mantle. The explanation is not very clear. From the mineral physics point of view, it is not unlikely that this could be explained in terms of thermal variation, even though we are unwilling to rule out the possibility of large wavelength compositional variations.     

3-D spherical shell modeling of mantle flow and its implication for global tectogenesis

In order to study the relationship between mantle flow and global tectogenesis, we present a 3-D spherical shell model with incompressible Newtonian fluid medium to simulate mantle flow which fits the global tectogenesis quite well. The governing equations are derived in spherical coordinates. Both the thermal buoyancy force and the self-gravitation are taken into account. The velocity and pressure coupled with temperature are computed, using the finite-element method with a punitive factor. The results show that the lithosphere, as the boundary layer of the earth's thermodynamic system, moves with the entire mantle. Both its horizontal and vertical movements are the results of the earth's thermal motion. The orogenesis occurs not only in the collision zones at the plates' boundaries, but also occurs within the plates. If the core-mantle boundary is impermeable and the viscosity of the lower mantle is considerable, the vertical movement is mostly confined to the upper mantle. The directions of the asthenospheric movements are not fully consistent with those of the lithospheric movements. The depths of spreading movements beneath all ridges are less than 220 km. In some regions, the shear stresses, acting on the base of the lithosphere by the asthenosphere, are the main driving force; but in other regions, the shear stresses are the resisting force.     

3-D velocity images of the crust and upper mantle of the Tianshan area

（胥颐，朱介寿，刘志坚，张华卿，朱燕）３－ＤｖｅｌｏｃｉｔｙｉｍａｇｅｓｏｆｔｈｅｃｒｕｓｔａｎｄｕｐｐｅｒｍａｎｔｌｅｏｆｔｈｅＴｉａｎｓｈａｎａｒｅａ￥ＹｉＸＵ；Ｊｉｅ－ＳｈｏｕＺＨＵ；Ｚｈｉ－ＪｉａｎＬＩＵ；Ｈｕａ－ＱｉｎｇＺＨＡＮＧａｎｄＹａ...     

Seismic waves in 3-D: from mantle asymmetries to reliable seismic hazard assessment

A global cross-section of the Earth parallel to the tectonic equator (TE) path, the great circle representing the equator of net lithosphere rotation, shows a difference in shear wave velocities between the western and eastern flanks of the three major oceanic rift basins. The low-velocity layer in the upper asthenosphere, at a depth range of 120 to 200 km, is assumed to represent the decoupling between the lithosphere and the underlying mantle. Along the TE-perturbed (TE-pert) path, a ubiquitous LVZ, about 1,000-km-wide and 100-km-thick, occurs in the asthenosphere. The existence of the TE-pert is a necessary prerequisite for the existence of a continuous global flow within the Earth. Ground-shaking scenarios were constructed using a scenario-based method for seismic hazard analysis (NDSHA), using realistic and duly validated synthetic time series, and generating a data bank of several thousands of seismograms that account for source, propagation, and site effects. Accordingly, with basic self-organized criticality concepts, NDSHA permits the integration of available information provided by the most updated seismological, geological, geophysical, and geotechnical databases for the site of interest, as well as advanced physical modeling techniques, to provide a reliable and robust background for the development of a design basis for cultural heritage and civil infrastructures. Estimates of seismic hazard obtained using the NDSHA and standard probabilistic approaches are compared for the Italian territory, and a case-study is discussed. In order to enable a reliable estimation of the ground motion response to an earthquake, three-dimensional velocity models have to be considered, resulting in a new, very efficient, analytical procedure for computing the broadband seismic wave-field in a 3-D anelastic Earth model.