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

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

全地幔三维异常密度的联合反演

由Hager(1981)提出的,叶正仁(1989)进一步发展的传播矩阵方法,提供了一种运用地球重力位、板块运动速度资料来探求地幔内部横向不均匀性的一条途径.但由于地球重力位、板块运动速度资料均是二维资料,故仅运用上述两种资料来反演三维异常密度会导致反演结果很强的不唯一性、不稳定性.作为80年代地球物理学最富有成效的成果之一就是利用地震P 波和S波的层析技术,提供了上、下地幔三维异常速度分布,从而为地球物理及地球动力学的研究提供了最新的先验知识.把传播矩阵方法、粘滞损耗最小原理运用于全地幔对流的地球动力学模型,并同时把长波地形起伏、板块运动速度、重力位资料以及傅容珊根据地震层析结果所合成的全地幔三维异常密度作为约束条件,联合反演全地幔的三维异常密度结构.计算结果显示:(1)从计算所及的密度横切剖面来看,至少在250—670km 深度范围内异常密度的高低与地表构造有一定的相关性:即全球的大部分发散边界和相当部分的大陆分别对应低异常密度区和高异常密度区,地幔的异常密度数值在1500km 深度附近达到最小.(2)从计算所得的密度纵切剖面来看,对应太平洋中脊、印度洋中脊、红海及亚丁湾地区,低密度的异常区从地表一直延伸到核一幔边界,不过太平洋中脊相对明显.除上述区域外全地幔的异常密度基本上(按数值大小及正负)分三层.     

利用低阶卫星重力位系数研究下地幔横向密度异常分布

本文试图采用卫星重力资料和一种新的反演方法来研究地幔的横向密度异常分布.先将密度异常△（r,,（?））在一个三维正交函数系下进行展开,其展开系数待定.然后,根据密度异常与重力扰动位之间的关系建立观测方程组,其中未知向量由密度异常展开系数组成,重力扰动观测向量由 GEM10B 重力模型中的位系数计算而得,并通过适当选取重力位系数的阶数,对观测向量进行滤波.最后,就下地幔（670km——CM 界面）作了实际计算.计算中,重力扰动位阶数取为2——11阶,密度异常展开式的截断阶数取为 K=4和 L=6,求解观测方程组时采用阻尼最小二乘法.结果表明:密度扰动值在670km 不连续面及核幔界面处达到极大值,且在环太平洋地区存在一高密度带,太平洋中部对应于一低密度区,这些特征与 Dziewonski 得到的下地幔三维波速异常分布特征相一致.但是,在南极地区、大西洋及印度洋部分地区,所得的密度异常分布与三维波速异常分布呈负相关,文章就其原因作了初步分析.      

地球重力模型中球谐函数阶数的异常源深度估计

利用大地水准面高度异常与重力异常的比值方法，通过选择合理的接近真实地球重力异常的异常源，为地球重力位球谐系数阶数与异常源埋深的关系进行了探讨，给出了异常源最大的深度范围。指出在核幔边界起伏的研究中，利用２－８阶重力位球谐系数的比较合理，这与Ｄｚｉｅｗｏｎｓｋｉ等人利用地球资料，在考虑地幔对流时推测的阶数具有较好的一致性。     

从克拉通破坏到板块动力模型的研究

华北克拉通破坏动力机制研究导致了全球动力系统及“板块动力模型”研究.板块运动最有可能的动力是地幔物质流动,但由于地幔物质流动的成因至今尚未查明,所以板块构造学说研究仍处于“运动模型”阶段,而没有进入“动力模型”阶段.如果地幔密度异常是驱动地幔物质流动的成因,那么就有可能基于重力学方法以“板块动力模型”的形式建立地幔密度异常驱动模式;软流圈中可能存在着动力特性不同的区块,地幔密度正异常代表物质盈余、区块内的物质要向区块外移动,地幔密度负异常代表物质亏损、区块外的物质要向区块内移动.本项目采用重力和地震资料相结合研究地球的整体分层,根据重力大地水准面联合地震波速度结构反演求解地幔密度异常,再根据地球正常密度假说和板块运动重力学机制的观点并与现有“板块运动模型”相结合,分析地幔密度异常动力区块,初步建立“基于重力学机制的板块动力模型”;为最终建立多学科机制的“全球板块动力模型”,迈出重要一步.     

全球地表热流的产生与分布

全球地表热流是反映地球内部热与动力学过程的一种主要能流.本文在三维球坐标框架下，就几个不同的粘度模型分别研究地幔内部密度异常(基于全球地震层析结果)以及板块运动激发的地幔流动的热效应及其对于观测地表热流产生和分布特征的贡献.由于地幔动力系统具有较高的Pe数，可以期望由板块运动激发的地幔流动将强烈地扰动地幔内部初始传导状态下的温度场以及地表热的热流分布.结果表明，与地幔内部密度异常产生的热效应相比，运动的板块及其激发的地幔流动在全球地表观测热流的产生和分布特征上起着更为重要的作用.观测到的大洋中脊处的高热流在很大程度上可以归因于板块激发的地幔流动的热效应.计算的平均温度剖面较好地揭示了岩石圈和D″层的温度特征，即温度随深度的剧烈变化，这与我们目前通过其他手段对岩石圈和D″层的温度结构了解是一致的.一个下地幔粘度比上地幔高出30倍的粘度结构(文中使用的粘度模型2)较之其余模型的拟合程度似乎更好.     

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

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

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

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

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

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

核－幔边界的动力学背景

根据传播矩阵方法，并把由联合反演得到的同时满足长波地形起伏、板块运动速度、重力位异常资料以及地震层析先验知识的全地幔三维异常密度作为载荷，以求取核－幔边界的动力学背景．计算结果显示：１．所求的核－幔边界起伏图像与Ｈａｇｅｒ等根据格林函数方法所求得的核－幔边界的起伏在全球范围内基本相符．２．核－幔边界处的环型场仅在数量上降为地表处的环型场的１／８左右，而极型场较地表处的极型场的流动图像有显著变化，数值也增为地表处极型场的３倍左右．     

Possible effects of lateral viscosity variations induced by plate-tectonic mechanism on geoid inferred from numerical models of mantle convection

In a traditional analytical method, the convective features of Earth’s mantle have been inferred from surface signatures obtained by the geodynamic model only with depth-dependent viscosity structure. The moving and subducting plates, however, bring lateral viscosity variations in the mantle. To clarify the effects of lateral viscosity variations caused by the plate-tectonic mechanism, I have first studied systematically instantaneous dynamic flow calculations using new density-viscosity models only with vertical viscosity variations in a three-dimensional spherical shell. I find that the geoid high arises over subduction zones only when the vertical viscosity contrast between the upper mantle and the lower mantle is O(10³) to O(10⁴), which seems to be much larger than the viscosity contrast suggested by other studies. I next show that this discrepancy may be removed when I consider the lateral viscosity variation caused by the plate-tectonic mechanism using two-dimensional numerical models of mantle convection with self-consistently moving and subducting plates, and suggest that the observed geoid anomaly on the Earth’s surface is significantly affected by plate-tectonic mechanism as a first-order effect.     

One-dimensional modelling of mantle flow

A one-dimensional model of flow between a fixed boundary at the bottom and a moving one on top with no net flow through vertical sections is tested for geophysically interesting mantle viscosity-depth functions. Such a model, although simplistic, may help in answering the question to what depth the return flow extends, at least in the case of moving plates measuring many thousand kilometers across, such as the Pacific plate.It the viscosity in the asthenosphere is less than three orders of magnitude smaller than that of the mantle below, the return flow extends to great depth and the asthenosphere is a zone of concentrated shear. If the viscosity contrast is greater, the return flow is concentrated in the asthenosphere. For a wide range of model parameters typical flow velocities below the asthenosphere are about one-tenth of the plate velocity. The pressure gradient required by the mantle flow may be manifest in gravity trends across moving plates, but no excessive gravity anomalies are required by the model if the absolute viscosity values conform to those inferred from post-glacial rebound data. A thinner and lower-viscosity layer is favored over a thicker and more viscous layer if both fit glacial rebound evidence. The present model may not be applicable if down to the core the viscosity is as low as about 10²¹ N s m^–2 with a free-slip bottom boundary.     

Lateral density anomaly distribution in the lower mantle obtained from satellite low order geopotential coefficients

In this paper, we attempt to use satellite gravity data and a new inversion method to study the lateral density anomaly distribution in the mantle. First, density difference Δρ(τ,θ,φ) is expanded in terms of a three—dimensional orthogonal function system, the coefficients of the expansion are to be determined. Then, a set of observation equations is established from the relationship between density anomaly and disturbing geopotential. In the equations the unknown vector contains the coefficients of density anomaly expansion, the observational vector is obtained by computing geopotential perturbations using the potential coefficients of GEM10B, and a filtering process is done for the observational values by properly selecting the harmonic degrees of geopotentical. Finally, the lateral density variations in the lower mantle (670 km toCM boundary) are investigated. In this case, the degrees of disturbing geopotential are selected as 2–11, the truncated degrees of density anomaly expansion are taken asL=6 andK=4, and the damping least squares method is used to solve the observation equations. The resulting model shows the high level of density perturbations at 670 km discontinuity and core — mantle boundary, a high — density zone circumscribing the Pacific and a lower—density region under the center of Pacific. These features are in agreement with the three—dimensional seismic velocity variation features by Dziewonski (1984). In the Antarctic region and some parts of Atlantic and Indian Ocean, however, the resulting density anomalies are negatively correlated with the seismic velocity anomalies, the cause resulting in these phenomena is preliminarily analysed in this paper. The Chinese version of this paper appeared in the Chinese edition ofActa Seismologica Sinica,13, 53–65, 1991. The principle and method represented in this paper can also be suitable to study the lateral density anomaly distribution in the earth’s crust and the upper mantle.     

黑水─泉州地学断面的重磁解释

讨论华南黑水─泉州地学断面的重磁解释。在解释中除应用常规的２．５维重磁异常反演外，还采用等密度线算法构制梯度层密度模型，用于分辨地壳密度细结构，以及计算自由地幔而深度用于分析上地幔密度的横向不均匀性。重磁模拟结果揭示了１０ｋｍ深度内的上地壳构造以及地壳与上地幔密度的分布。提出华南造山带以低密度上地幔为特征，它可能与上地幔的改造有关；四川盆地具有较高的上地幔密度，为未经改造的原始地幔。扬子克拉通与华南造山带的分界线与上地幔密度的分界线一致。根据布格异常、地表岩石密度和地形资料的综合分析，圈定出反映内生成矿作用深部标志的密度倒转区，可作为进一步找矿的远景区。     

黑水─泉州地学断面的重磁解释

讨论华南黑水─泉州地学断面的重磁解释。在解释中除应用常规的２．５维重磁异常反演外，还采用等密度线算法构制梯度层密度模型，用于分辨地壳密度细结构，以及计算自由地幔而深度用于分析上地幔密度的横向不均匀性。重磁模拟结果揭示了１０ｋｍ深度内的上地壳构造以及地壳与上地幔密度的分布。提出华南造山带以低密度上地幔为特征，它可能与上地幔的改造有关；四川盆地具有较高的上地幔密度，为未经改造的原始地幔。扬子克拉通与华南造山带的分界线与上地幔密度的分界线一致。根据布格异常、地表岩石密度和地形资料的综合分析，圈定出反映内生成矿作用深部标志的密度倒转区，可作为进一步找矿的远景区。     

Numerical model for the generation of the ensemble of lithospheric plates and their penetration through the 660-km boundary

In the kinematic theory of lithospheric plate tectonics, the position and parameters of the plates are predetermined in the initial and boundary conditions. However, in the self-consistent dynamical theory, the properties of the oceanic plates (just as the structure of the mantle convection) should automatically result from the solution of differential equations for energy, mass, and momentum transfer in viscous fluid. Here, the viscosity of the mantle material as a function of temperature, pressure, shear stress, and chemical composition should be taken from the data of laboratory experiments. The aim of this study is to reproduce the generation of the ensemble of the lithospheric plates and to trace their behavior inside the mantle by numerically solving the convection equations with minimum a priori data. The models demonstrate how the rigid lithosphere can break up into the separate plates that dive into the mantle, how the sizes and the number of the plates change during the evolution of the convection, and how the ridges and subduction zones may migrate in this case. The models also demonstrate how the plates may bend and break up when passing the depth boundary of 660 km and how the plates and plumes may affect the structure of the convection. In contrast to the models of convection without lithospheric plates or regional models, the structure of the mantle flows is for the first time calculated in the entire mantle with quite a few plates. This model shows that the mantle material is transported to the mid-oceanic ridges by asthenospheric flows induced by the subducting plates rather than by the main vertical ascending flows rising from the lower mantle.     

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.     

On the importance of advection in determining the local isotopic composition of the mantle

The equations for isotopic evolution in a deforming medium are derived and used to show that the local isotopic composition of the mantle depends primarily on the advection of heterogeneities by the flow field. Various examples of the role of advection are given using a two-scale model of mantle convection. The effect of small-scale flows, which have dimensions smaller than the plates themselves, is to convert any initially localized heterogeneity into thin spiral sheets. The isotopic properties of erupted lavas will depend on the relative size of the zone of partial melt to the spacing between such sheets. An idealized model for⁸⁷Sr/⁸⁶Sr variations from Pacific islands and the covariation of¹⁴³Nd/¹⁴⁴Nd-⁸⁷Sr/⁸⁶Sr is based on the combined effect of isotopic fluxes at subduction zones and advection by the return flow under the moving lithospheric plates.