Coupling between different layers of the lithosphere under horizontal drag underneath

Based on a layered Theological model of the lithosphere, the velocity and stress distributions in the lithosphere under horizontal drag underneath were calculated using viscoelastic finite element method of plain strain with finite deformation. In the simulation, different conditions of drag and blocking were assumed to study their influences on the stress distribution and the coupling between different layers. Blocking depth has little influence on the stress level in the whole area and the coupling between different layers, but influences the stress state in the area around the blocking. The area covered by the high stress anomaly becomes larger when the blocking depth becomes deeper, but the magnitude of the value of the maximum shear stress decreases. The greater the viscosity differences between different layers of the lithosphere, the greater the possibility of decoupling between them. Under the drag of normal mantle convection (the convection velocity is about 20 cm · a^?1), a lithosphere with a Theological structure similar to that of North China could not have decoupling between different layers, while could have stress distribution with magnitude of several MPa to tens of MPa and could have anomalous areas with stress accumulation if the geological structure is complicated.     

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.     

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

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

Three-dimensional numerical thermal and rheological modelling in the central Fennoscandian Shield

A three-dimensional model for the central Fennoscandian Shield was constructed for analysing the thermal, the rheological and the structural conditions in the lithosphere. The mesh covers a rectangular area in the southern Finland with horizontal dimensions of 500 km × 400 km and a depth extent of 100 km. Structural boundaries are derived from the several deep seismic soundings carried out in the area. Constructed model is first used in the calculation of the thermal and the rheological models and secondly in analysing the stress and the deformational conditions with the obtained rheology. Thermal and structural models are solved with the finite element method. The calculated surface HFD is between 40 and 48 mW m⁻² in the Proterozoic southern part and below 40 mW m⁻² in the older and northern Archaean part of the model. The calculated rheological strength shows a layered structure with two individual rheologically weak layers in the crust and strong layer in the upper part of the lower crust. The minimum brittle–ductile transition (BDT) depth is around 10 km in the southern part of the model while in the north and north-eastern parts the BDT depth is around 45–50 km. Comparison with the focal depth data shows that as most of the earthquakes occur no deeper than the depth of 10 km are they located in the brittle regime. Resulting stress conditions and possible regions of deformation after the model is subjected to pressure of 50 MPa reveals that the stress field is quite uniformly distributed in different crustal layers and that the elastic parameters control more the state of the stress than the applied rheological structure. In the upper crust, the stress intensity has values between 42 and 45 MPa whereas in the middle crust the values are around 50 MPa. Comparison of the 3-D model with earlier 2-D models shows that some differences in the results are to be expected.     

区域重力异常和上地幔小尺度对流相关方程及对流拖曳力场

建立了区域重力均衡异常的物理模型，假定区域重力均衡异常源于上地幔密度分布横向不均匀以及由上地幔小尺度对流产生的边界形变．导出了区域重力均衡异常和上地幔小尺度对流之间的相关方程，讨论了该方程在地球动力学研究中的应用．应用几种模型计算上地幔小尺度对流产生的作用于岩石层底部的拖曳力场．结果表明，将这一相关方程用于利用区域重力均衡异常研究地壳及上地幔应力状态，将有广阔的前景．     

华北地区中东部岩石圈挠曲与均衡特性以及地震活动性分析

华北地区中东部涵盖北京、天津以及即将建设的雄安新区等大型城市,区内发育了张渤地震带等多条大型活动断裂,地震活动性较强,历史上发生过多次6级以上地震.本文利用Fan小波的布格重力异常一致性方法研究该区的岩石圈有效弹性厚度和均衡调整初始加载比分布,同时基于均衡调整方法计算该区垂向构造应力分布,并将以上结果与历史地震活动进行统计分析.岩石圈挠曲分析表明,华北地区中东部的岩石圈有效弹性厚度为10~65 km,分布特征为自东南向西北逐渐减小.均衡调整初始加载比为0.5~0.8,表明现今的岩石圈挠曲状态主要由莫霍面加载形成.该区地壳承载的垂向构造应力约为-20~20 MPa,中西部地区垂向构造应力向上,东北和西南地区向下.统计分析结果显示,华北地区中东部的地震活动性随着岩石圈有效弹性厚度和均衡调整初始加载比的增加而减弱,垂向构造应力零值区域地震活动性较弱.雄安新区的岩石圈有效弹性厚度大约为15 km,均衡调整初始加载比为0.5~0.6,垂向构造应力为15~20 MPa,岩石圈参数对应的地震活动性较强,相关结果对于新区建设具有一定参考价值.     

The global stress field in the lithosphere obtained from the satellite gravitational harmonics

In this paper, the lithosphere is considered to be a homogeneous elastic spherical shell for the sake of simplicity and the stress equations for the base of the lithosphere are taken as boundary conditions. Then the stress equations are obtained for use in the computation of the stress field in the lithosphere with the satellite gravitational harmonic coefficients. The 5 × 5° global stress field in the lithosphere is computed from harmonics of 2–30°. The directions of principal stresses of this stress field agree favourably with the directions of principal stresses indicated by mid-plate earthquake mechanisms, in situ stress measurements and sensitive geological features. This result indicates that the drag forces exerted on the base of the lithosphere, due to gravitational mantle convection, may be among the driving forces of plate motion and a major source for the stress field in the lithosphere.     

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

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

岩石圈流变强度与中国大陆构造运动关系的探讨

以GPS观测资料和地震学研究成果为约束，针对不同流变参数的中国大陆岩石圈模型，数值模拟了岩石粘度与中国大陆板块边界作用强度的关系，探讨了陆-陆碰撞对中国大陆分层岩石圈运动的驱动机制．给出了陆-陆碰撞驱动力、附加地形与山根浮力及热浮力对中国大陆构造运动的驱动特点．印度板块、太平洋板块和菲律宾板块对中国大陆驱动的边界作用强度之比约是4：1．25：1，所引起的水平主压应力主要集中在坚硬岩石层；而附加地形等垂直方向作用力在水平方向产生的最大主压应力则主要集中在软弱岩石层．这种垂直方向上的作用力在高原南部地区阻碍陆-陆碰撞向北的推挤运动，在高原东北部增加对其它块体的推挤作用。     

中国大陆与各活动地块、南北地震带实测应力特征分析

本文以"中国大陆地壳应力环境基础数据库"为基础,补充了迄今为止查阅到的中国大陆水压致裂法与应力解除法的实测地应力数据,在1474个测点上得到3586条数据,研究区经度范围75°E－130°E,纬度范围18°N－47°N,深度范围0~4000 m,基本覆盖了中国大陆的各活动地块与南北地震带各段等研究区.本文采用等深度段分组归纳的方法解决了实测地应力数据样本数量沿深度分布的不均匀问题,给出了中国大陆与各研究区地壳浅层测量深度范围内应力量值、方位特征.结果显示:(1)中国大陆地壳浅层最大水平应力、最小水平应力、垂直应力随深度呈线性增加;(2)中国大陆地区侧压系数随深度的变化特征为:浅部离散,随着深度增加而集中,并趋向0.68,D=465 m是水平作用为主导向垂直作用为主导的转换深度,K_av=1;(3)中国大陆水平差应力在地表为3 MPa左右,随深度增加以5.8 MPa/km的梯度增大;(4)在测量深度范围内,中国大陆各研究区最大水平应力中间值(深度为2000 m时的统计回归值)从大到小的顺序是:青藏地块63.6 MPa、南北带北段57.3 MPa、华南地块51.4 MPa、华北地块50.5 MPa、南北带中段47.9 MPa、西域地块47.5 MPa、南北带南段45.4 MPa、东北地块44.8 MPa,总体表现为"西强东弱"的基本特征,反映了印度板块与欧亚板块的强烈碰撞是中国大陆构造应力场强度总体特征的主要来源;(5)与其他研究区相比较,青藏地块地壳在从南向北的挤压作用下呈现出明显的"浅弱深强"特点;(6)最大水平应力方向的总体特征,基本以青藏高原为中心,呈辐射状展布,由西向东,从近N-S方向逐步顺时针旋转至NNE-SSW、NE-SW、NEE-SWW、NW-SE方向,与深部的震源机制解研究结果有一致性.     

Implications of melting for stabilisation of the lithosphere and heat loss in the Archaean

The increased depth and volume of melting induced in a higher temperature Archaean mantle controls the stability of the lithosphere, heat loss rates and the thickness of the oceanic crust. The relationship between density distributions in oceanic lithosphere and the depth of melting at spreading centres is investigated by calculating the mineral proportions and densities of residual mantle depleted by extraction of melt fractions. The density changes related to compositional gradients are comparable to those produced by thermal effects for lithosphere formed from a mantle which is 200°C or more hotter than modern upper mantle. If Archaean continental crust formed initially above oceanic lithosphere, the compositional density gradients may be sufficient to preserve a thick Archaean continental lithosphere within which the Archaean age diamonds are preserved. The amount of heat advected by melts at mid-ocean ridges today is small but heat advected by melting becomes proportionally more important as higher mantle temperatures lead to a greater volume of melt and as the rate of production of oceanic plates increases. Archaean tectonics could have been dominated by spreading rates 2–3 times greater than now and with mantle temperatures between ca. 1600°C and 1800°C at the depth of the solidus. Mid-ocean ridge melting would produce a relatively thick but light refractory lithosphere on which continents could form, protected from copious volcanism and high mantle temperatures.     

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

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

汶川8.0级地震环境剪应力特征研究

基于甘肃强震固定台和流动观测记录的汶川8.0级地震主震及余震加速度资料,选用三分向记录均完整的87次地震(震级范围为3.0~8.0级),根据用位错理论二维断裂模式推导的震源峰值加速度与环境剪应力关系式计算汶川地震序列的环境剪应力值,探讨其环境剪应力场特征。结果表明,环境剪应力和矩震级有较好的相关性,3~5级地震对应的应力值多在3~6 MPa,5~6级多在6~9 MPa,6~7级多在9MPa以上,且环境剪应力对震源深度也有较强的依赖性。     

长白山火山区壳幔S波速度结构研究

利用面波层析成像和远震接收函数方法对长白山地区的地壳上地幔速度结构进行了研究。结果表明:长白山火山区附近存在岩石圈减薄、上地幔软流圈增厚以及上地幔S波速度降低等与上地幔高温物质有关的现象,它表明长白山的岩浆系统一直延伸到上地幔软流圈范围。天池火山区地壳内部存在明显的S波低速层,在离天池火山口较近的WQD台附近,低速层顶部埋深约8km,厚度近20km,S波最小速度约2.2km/s。在距离天池火山北部50km的EDO台地壳中没有明显的低速层。火山区S波速度结构总体表现出距离天池越近,地壳的V_P/V_S越大,低速层的厚度和幅度增加的特征,表明天池火山口附近地壳内部存在高温物质或岩浆囊。CBS台站不同方位的接收函数及反演结果表明,地表低速层厚度以及莫霍面深度存在随方位的变化。地表低速层在南部方向明显较厚,莫霍面深度在南部天池火山口方向存在小幅度抬升。CBS台站附近特殊的近地表速度结构可能是该台站记录的火山地震波形主频较低的主要因素。天池火山口附近莫霍面的小幅度抬升意味着存在与火山作用有关的壳幔物质交换通道     

峰值速度和加速度对环境剪应力的依赖性

从导出的地震定标律和地震破裂过程的断裂力学模式出发,得到了震源平均位移、速度和加速度谱的表达式,进而又推导出震源的峰值位移d_m,速度v_m和加速度a_m的表达式:d_m=k_dM_0~2/~3τ_0~(2/3),v=k_vM_0~(1/3)τ_0~(4/3),a_m=k_aτ_0~2式中M_0是地震矩,τ_0是环境剪应力值,k_d,k_v,k_a为适当的常数.我们选用了66个地震的观测资料,这些地震的矩震级范围包括了从1级左右的极微震,3-5级的小震,直到6-7级的大震;地震矩从10~9-10~(20)Nm,跨越了10个数量级,并用这些地震检验了上述公式. 令所有地震的平均应力为5MPa,定出常数k_d,k_v,k_a,进而由速度和加速度观测资料求得66个地震的环境剪应力τ_0值,这些数值相当稳定.多数极微震的τ_0值在2-4MPa之间;小震的τ_0值多数在4-8MPa左右;大震的τ_0值为10MPa左右。τ_0值对震源深度和断层类型有明显的依赖性.一般深度很浅的小震和极微震,τ_0值很低;正断层地震的τ_0值相对较低;逆断层地震的τ_0值较高;走滑断层地震的τ_0值则居中.     

A laboratory study of teleseismic Pn waves

In this paper, the influences of mechanical properties of low velocity layers and thickness of high velocity layers in the lower lithosphere on the horizontal propagaties of seismic waves are discussed by means of altrasonic seismic model experiments; It is shown that the structure of the lower lithosphere consisting of soft low velocity layers and thin high velocity layers is favourable for the horizontal propagation of seismic energy, therefore, it is a very possible structure pattern of the lower lithosphere responsible for Pn and later high velocity arrivals observed in long range distances; Such pattern is also consistent with studies of rheological stratification of the lower lithosphere beneath mid-ocean-ridges and continental rifts. The Chinese version of this paper appeared in the Chinese edition ofActa Seismologica Sinica,15, 97–102, 1993. The project is supported by the Chinese joint Seismological Science Foundation.     

Three-dimensional numerical thermal and rheological modelling in the central Fennoscandian Shield

A three-dimensional model for the central Fennoscandian Shield was constructed for analysing the thermal, the rheological and the structural conditions in the lithosphere. The mesh covers a rectangular area in the southern Finland with horizontal dimensions of 500 km × 400 km and a depth extent of 100 km. Structural boundaries are derived from the several deep seismic soundings carried out in the area. Constructed model is first used in the calculation of the thermal and the rheological models and secondly in analysing the stress and the deformational conditions with the obtained rheology. Thermal and structural models are solved with the finite element method. The calculated surface HFD is between 40 and 48 mW m⁻² in the Proterozoic southern part and below 40 mW m⁻² in the older and northern Archaean part of the model. The calculated rheological strength shows a layered structure with two individual rheologically weak layers in the crust and strong layer in the upper part of the lower crust. The minimum brittle–ductile transition (BDT) depth is around 10 km in the southern part of the model while in the north and north-eastern parts the BDT depth is around 45–50 km. Comparison with the focal depth data shows that as most of the earthquakes occur no deeper than the depth of 10 km are they located in the brittle regime. Resulting stress conditions and possible regions of deformation after the model is subjected to pressure of 50 MPa reveals that the stress field is quite uniformly distributed in different crustal layers and that the elastic parameters control more the state of the stress than the applied rheological structure. In the upper crust, the stress intensity has values between 42 and 45 MPa whereas in the middle crust the values are around 50 MPa. Comparison of the 3-D model with earlier 2-D models shows that some differences in the results are to be expected.     

中国及其邻区地球三维结构初始模型的建立

对人工地震测深及天然地震面波体波三维层折反演数据进行统一处理，建立了中国及其邻区地球三维结构初始模型．此模型图像表明，中国及其邻区地球各圈层横向变化明显．岩石圈及软流圈内速度分布主要反映这一区域自古生代以来板块及地块拼合模式．各主要板块或地块（塔里木、扬子、中朝、青藏、哈萨克斯坦、印度、印度支那）岩石圈增厚或有很深的地慢根，板块或地块间的造山带岩石圈减薄，软流圈速度降低。下地幔底部及核幔边界D″层出现高速异常，表明古太平洋及古特提斯洋俯冲板块因重力坍塌已进入地球深层，形成亚洲超级下降地幔柱。这一下降地幔柱引起地球表层物质向中亚、东亚地区集中，印度半岛、青藏高原、新疆、蒙古至贝加尔一带，成为全球岩石圈最大的汇聚场所．     

华北地区地下介质波速比值(VP/VS)研究

文中利用分别独立反演得到的华北地区 16个台的S波和P波速度结构 ,计算了各台站下方分层速度结构各层的波速比值 ,发现在强震区的波速比总体上高于弱震区 ,且地壳中的一些低速层位会出现 >1 8的波速比值 ,而在弱震区或稳定块体内部 ,各层波速比一般都低于正常波速比 1 732 ,平均值甚至 <1 7;在过渡区域 ,如位于盆地边缘或稳定块体边缘的台站 ,其波速比值有高有低 ,一般情况下 >1 732。在此基础上 ,最后我们对位于不同地质构造体的台站下方的介质性质进行了讨论     

俯冲带耦合作用对苏门答腊地区应变场影响的三维数值模拟

采用有限元方法模拟了俯冲带耦合作用对巽他弧及其邻区的影响.根据模拟结果,对比GPS、地震和地质学观测数据,定量分析了苏门答腊及其周边地区的应变强度和主应变方向的分布特征,据此探讨了该区构造特征、地震发生模式与耦合面积之间的关系.模型由具有黏弹性性质的岩石圈和软流圈上地幔组成,其中岩石圈包括了大陆岩石圈和大洋岩石圈以及俯冲至上地幔中的俯冲板片.研究结果如下:(1) 通过对不同俯冲带耦合面积模拟,发现苏门答腊前弧伴随耦合面积的增加应变强度逐渐增大,而增大的应变强度又影响了其周边地区的应变分布特征,因此整个苏门答腊前弧呈现出明显的分段性,这与该区地震破裂模式有较好的对应.(2)苏门答腊北部地区主应变方向与南部相比存在一定的差异,该差异是俯冲带的俯冲方向、俯冲速度、俯冲形态以及不同区域间耦合面积共同作用的结果.(3)虽然苏门答腊2004年地震主震区处于弱耦合状态,但从本文模拟的结果中可以看到,在俯冲作用下该区依然存在垂直向下的位移,这为地震激发海啸提供了有利的构造环境.