青藏高原松潘—西秦岭—临夏盆地深地震反射剖面——采集、处理与初步解释

由于印度洋板块向亚欧板块俯冲使青藏高原不断隆起,其形成不仅导致了亚洲大陆内部强烈的晚新生代构造变形,还对其边缘地区的地貌格局产生重大影响.青藏高原东北缘是青藏高原向北东方向扩展的前缘部位,是印度与欧亚两大板块碰撞作用由近南北方向向北东、东方向转换的重要场所.本文利用2004年和2008年完成的深地震反射剖面资料,采用关键处理技术和参数开展唐克-合作剖面与合作-临夏剖面联线处理,获得总长约400 km的深地震反射剖面,完整揭示了西秦岭造山带及其两侧盆地的地壳结构和构造变形样式.结果显示西秦岭造山带下地壳向若尔盖逆冲推覆的深部构造特征;西秦岭下地壳北倾的强反射及其北侧南倾的强反射特征揭示出扬子与华北两个大陆板块在西秦岭造山带下的汇聚行为.Moho的埋深和起伏形态表明青藏高原东北缘地壳经历了高原隆升后强烈的伸展减薄作用.     

岩石圈根拖曳作用在青藏高原西部和东部造山中的比较研究

根据地球物理探测资料,基于二维模型,利用黏弹性有限元方法,研究青藏高原西、东剖面的地壳均衡和岩石圈根拖曳的构造应力机制.数值模拟结果表明,青藏高原西部（B B′剖面）的造山水平挤压力主要来源于岩石圈根的向下拖曳,印度板块向北挤压为次要因素,形成“山隆盆降”的地表形态;而青藏高原东部（A A’剖面）岩石圈根向下拖曳还不足以形成硬上地壳中挤压造山的主要力源.对比结果认为,青藏高原的深部层圈结构和应力体系在西、东部存在明显的差异,反映了高原内部造山演化的西、东分异特征.     

青藏高原东南缘地壳结构与变形机制研究进展

新生代青藏高原的隆升改变了整个亚洲的构造格局,对气候、环境均产生了重要的影响,但高原的隆升扩展机制众说纷纭.青藏高原东南缘作为扩展前缘,其构造演化对了解整个高原的扩展机制具有重要的意义.本文总结了近年来对青藏高原东南缘地壳结构研究的最新进展,特别是2011年中国地震科学探测台阵计划开展以来,利用密集地震台阵取得的新成果,探讨了青藏高原东南缘地壳的结构与变形机制.这些研究发现青藏高原的地壳由高原向外围减薄,但在高原边界断裂附近存在地壳厚度突变带;下地壳中存在两个独立的低速异常,一个位于松潘—甘孜块体下方,被高原的边界断裂所围限,另一个位于小江断裂带下方,呈NE-SW向展布.我们认为青藏高原东南缘下地壳物质被边界(丽江—小金河)断裂所围限,并没有继续向边缘流出,但是地壳挤出产生的应力作用继续向东南方向传递,造成了小江断裂带附近的地壳变形.     

甘肃东南地区构造活动与2013年岷县—漳县MS6.6级地震孕震机制

位于南北地震带中北段的甘东南地区,其构造变形和构造活动特征与青藏高原向北东方向的扩展密切相关,该地区复杂的构造几何形态主要受控于东昆仑断裂和西秦岭北缘断裂,区域新构造运动主要动力来源于青藏高原向北东的扩展.近年来,甘东南地区中强地震频发,本文主要通过对该地区构造活动特征、历史地震等资料的综合分析讨论,结合地球物理、地震学和野外调查等资料,认为青藏高原东北部东昆仑断裂的向北挤压和向东的运动是该地区构造应力集中的主要原因,也是该地区中强地震的主要孕震环境和机制,而西秦岭北缘断裂的走滑及向南北两侧逆冲“花状构造”是临潭—宕昌断裂带上中强地震频繁发生的一个重要动力因素.2013年7月22日发生在甘肃岷县—漳县的M_S6.6级地震正好位于临潭—宕昌断裂带中东段上,是该断裂分段不均匀活动的结果.     

龙门山南段前陆区晚第四纪构造变形样式

龙门山中南段前陆区是青藏高原东缘唯一发育新生代薄皮构造与沉积盆地的地段,研究其最新构造变形样式有助于深入理解青藏高原向东扩展的构造机理.论文通过青衣江河流阶地测量与古青衣江洪积扇形态重建,研究了龙门山南段前陆区晚第四纪活动构造格局及其活动性,取得了如下认识：（1）青衣江河流阶地纵剖面显示,龙门山南段前陆地区晚第四纪变形主要为褶皱作用,总体地壳缩短速率为2.5~3.9 mm·a^-1,远大于山区冲断带0.48~0.77 mm·a^-1的地壳缩短速率,地壳缩短主要由前陆地区吸收;（2）青衣江古洪积扇错断变形显示,龙门山南段前陆区活动构造表现为北西—南东向地壳缩短与近东西向的地壳缩短的叠加作用,两者分别受控于巴颜喀拉块体南东向推挤作用与川滇块体向东推挤作用;（3）自中新世初川滇块体向南东挤出,四川盆地西南角起到分流青藏高原物质的作用,其西南侧物质通过鲜水河—小江断裂带的左旋错动向南东方向分流,其西北侧物质通过龙门山断裂带的右旋错动向北东方向分流,迎面受到了最大的推挤作用,进而向前陆扩展形成了薄皮褶皱构造带.     

木里—盐源地区深部电性结构及构造意义

木里—盐源地区地处青藏高原东南缘,属于古特提斯洋构造域,是松潘—甘孜地块及扬子地块的交接地带,是研究青藏高原东南缘构造演化过程的重要区域.本文介绍的是横穿木里—盐源地区的大地电磁剖面,自北西向南东依次跨越锦屏山断裂、木里弧形构造区、丽江—小金河断裂、盐源盆地、金河—箐河断裂等构造.维性分析表明木里弧形构造区和金河—箐河...     

青藏高原东北缘和鄂尔多斯地区的均衡重力异常

青藏高原东北缘地区是青藏高原块体与华北板块和华南板块相互作用的交汇地区.在青藏高原隆升并向北东方向扩展的过程中,与华北板块和华南板块发生相互碰撞、挤压,在这里形成了强烈的构造变形和复杂的构造格局,这里也是中国大陆上强震多发的地区,己经多次发生8级和更大地震.     

维西—贵阳剖面重力异常与地壳密度结构特征

维西—贵阳剖面位于青藏高原东南缘,为青藏高原物质往东南逃逸、东构造结侧向挤压及华南地块北西西向推挤作用的重要地段.利用剖面观测的重力与GPS定位数据,结合区域背景重力场、地质构造及深部地球物理成果,反演研究剖面较为细化的地壳密度结构特征.观测研究表明:剖面布格重力异常总幅差变化达190×10-5 m·s-2,具"斜N"分段变化特征,从西往东呈上升(维西至攀枝花,水平梯变大)—下降(攀枝花至会泽,水平梯变较大)—上升(会泽至贵阳,水平梯变较小)态势;高程与布格重力异常比值的趋势性转折部位为康滇地轴核心和小江断裂带东侧,可能与先存构造或新生构造发育有关;剖面地壳密度结构可分上、中和下三层结构,各层底界面平均埋深分别约20km、35km和51km,金沙江—红河断裂带和鲜水河—小江断裂带为地壳结构相对简单与复杂的过渡带;地壳厚度西深东浅,可能是东构造结的侧向挤压所致;下地壳厚度变化相对较大,可能对地壳增厚起主要作用;华坪—攀枝花附近的Moho面隆起和上地壳高密度体的存在暗示上地幔往上底侵作用,对青藏高原物质向南东逃逸和东构造结的侧向挤压均起到一定阻挡作用;中地壳下伏有限低密度薄层有利于其上物质的南东逃逸和顺时针旋转,有利于其下物质受喜马拉雅东构造结作用下往东向运移.     

川滇地区壳幔介质地震各向异性特征研究

川滇地区位于青藏高原东缘,由于印度-欧亚大陆碰撞造成物质的侧向挤出和逃逸,造成了该地区的构造格局由近东西向转为近南北向.该地区的地震活动强烈,地质构造十分复杂,对该地区的各向异性研究,对于该区域主压应力场特征、断层性质、地震成因、板块运动和强震深部背景等问题的研究具有重要意义,尤其对青藏高原的隆升演化的动力学研究有重要意义.     

中国布格重力异常水平梯度图的判读和构造解释

分析和对比中国大陆布格重力异常水平梯度的彩色图像，得到一些新的构造特征信息. 大致以105°E为界，图像的东部以蓝色低梯度为背景，出现中、高梯度带状构造和大梯度差值边界的线状构造；西部以高、中梯度的带状构造为格架，镶嵌着低梯度区，区内有阴影线状构造. 总体上，重力异常水平梯度构造可概括为三大组. 第一组是NNE向的大兴安岭、太行山、武陵山高梯度带与斜切东部全区的郯庐大断裂带和NNW向的斜切西部全区的三条大断裂带；第二组为西部的NWW向与NEE向高、中梯度带的交叉构造格局，以及东部近于等间距的三条NWW向断续的阴影构造带和一条单一的构造线，而华南则仅有NE向的低梯度构造线；第三组是横贯全区的的几条近EW向和中部SN向宽而复杂的构造带，在后者的东、西两侧散臣着大间距、小长度、断续的SN向阴影构造线. 这三组重力梯度构造不仅反映了地貌起伏的信息，而且显示了大区域构造，特别是深部构造的密度差. 各组线状构造带的内部结构和各组的形成时间及演变历史可能不同，但它们相互复合构成现今总体的构造格局，对大尺度现代构造应力场作用都会有不同的响应，因此导致主应力方向的区域性变化.     

青藏高原东西向差异形变与隆升机制

高精度布格重力异常约束下的三维空间域挠曲形变模拟显示,大约以90°E为界,青藏高原东、西两部分的岩石圈强度存在明显的差异.在90°E以东,岩石圈有效弹性厚度为35~45 km,该岩层厚度可使刚性的上地壳与上地幔岩石通过中下地壳柔塑性地层的黏滞流动产生构造解耦;地壳处于区域均衡状态,下地壳热物质的流动膨胀是地壳隆升的主控要素.而在90°E以西,断裂带严重削弱了该区域的岩石圈机械强度,岩石圈有效弹性厚度小于15 km,向西逐渐减小,至喀喇昆仑断裂带变为零,断裂切穿莫霍面进入地幔,发生纯剪切构造形变;这里的地壳接近局部均衡,厚皮逆冲是地形隆升的主要因素.震源深度大于80 km的地幔地震大多发生在青藏高原西部,其岩石圈深部具有的脆裂特征很好地支持了岩石圈机械强度模拟的结果.     

PRELIMINARY APPLICATION OF FOCAL MECHANISM SOLUTIONS OF SMALL AND MEDIUM-SIZE EARTHQUAKES TO FAULT STABILITY ANALYSIS IN THE SOUTHEASTERN TIBETAN PLATEAU

Analysis of stress state of faults is helpful to understand crustal mechanical properties and seismicity. In the paper, we invert the horizontal crustal stress field in the southeastern Tibetan plateau using focal mechanism solutions of small and medium-size earthquakes, and apply them to estimate the stability of regional major faults. Firstly, we collect focal mechanism solutions of small and medium-sized earthquakes in the southeastern Tibetan plateau. The dataset includes more than 1 000 focal mechanism solutions in the past twenty years. Magnitudes of these earthquakes vary from M3.0 to M6.0. Most of the focal mechanism solutions were determined using waveform inversion technique. Although most of focal mechanism solutions in the southeastern Tibetan plateau are strike-slip faulting, their spatial pattern is different in sub-regions. Normal faulting earthquakes mainly occurred in the western Sichuan region, reverse faulting earthquakes mainly occurred in the boundary zone between the Tibetan plateau and the South China craton, and strike-slip faulting earthquakes mainly occurred in the central and southern Yunnan region. Next, we settle on a mesh with grid spacing of 0.5° in longitude and latitude in the region and invert the horizontal crustal stress field at each grid point. Spatial variation of the maximum principal stress axis in the southeastern Tibetan plateau shows a clockwise rotation around the eastern Himalaya syntax. The azimuth of maximum compressional stress axis is about 88.1° in the western Sichuan region, about 124.6° in the South China craton, and about 21.6° in the western and southern Yunnan region. The azimuth of regional maximum compressional stress is nearly parallel to the direction of terrain elevation gradient, and that of the minimum compressional stress is nearly parallel to the tangential direction of the topographic elevation contours. The spatial pattern reflects the control role of gravity spreading of the Tibetan plateau on the regional horizontal stress field. Finally, we analyzed regional fault stability based on these collected focal mechanism solutions. The fault instability parameter (I) is defined based on the Mohr-Coulomb criterion and indicates the degree of fault approximating to rupture. The instability parameters on fourteen major faults in the southeastern Tibetan plateau were calculated. Our results show that the stability of the Lianfeng-Zhaotong Fault is the lowest before 2014 in the region, which indicates the fault zone is close to rupture at that time. Our results provide a new useful tool to assess regional seismic potential using dense focal mechanism solutions.     

青藏高原东北缘及周边地区基于程函方程的面波层析成像

利用"中国地震科学台阵探测"在南北地震带北段布设的670套宽频带地震台站记录到的面波资料,使用新近发展的程函方程面波层析成像方法,获得了青藏高原东北缘及周边地区12~60 s周期范围比以往成像结果具有更高分辨率的瑞利面波相速度分布图像.青藏高原东北缘的祁连褶皱系西段、秦岭褶皱系西段和松潘一甘孜褶皱系,在16~60s周期范围内均显示出明显的低速异常分布,表明该地区的地壳力学强度较低,在强烈的构造应力作用下易发生形变.与西段不同,祁连褶皱系东段和秦岭褶皱系中段的相速度分布特征揭示,其中下地壳的速度明显高于高原内部区域.鄂尔多斯块体整体上表现为稳定块体具有的高速特征,但其西部边缘在中上地壳的速度比块体中部地区偏低,且存在一定的不均匀性.鄂尔多斯块体西北缘的临河断陷盆地和西缘的银川断陷盆地,在较短的周期范围内(12~20 s)表现为局部低速特征,但与银川断陷盆地不同,临河断陷盆地的低速特征可一直延续至60 s周期以上,表明该盆地下方地壳及上地幔速度明显偏低,可能与深部热作用有关.阿拉善块体与其北部地区的速度差异主要表现在中上地壳,这一现象值得今后进一步探讨.基于程函方程面波层析成像方法给出了青藏高原东北缘及周边地区高分辨率的成像结果,揭示了以往面波层析成像难以获得的深部细节特征,为该地区的深部构造研究提供了新的信息.     

喜马拉雅西构造结区深部地球物理场特征剖析

本文基于喜马拉雅西构造结地带区的重力场、电磁波场和地震波场特征,从不同学科的角度对喜马拉雅西构造结区深部地球物理场做了分析.不同学科的研究成果表明,西构造结区不仅存在着令人瞩目的中、深源地震,而且其重力场、电磁波场和壳、幔地震速度结构存在着高度的非均匀性.这不仅表明其在不同演化阶段和不同物质组成先后参与了岛弧与大陆、大陆与大陆碰撞的造山过程及其导致的横向结构特征的显著差异上,同时也反映出这种复杂过程在对深层演化进程所带来的影响,故导致了西构造结区的形成和长期演化的复杂过程.显然西构造区的介质属性、构造格局和深层过程均制约着、并强烈影响着高原西北缘及其西部的深、浅构造组构和运动学与动力学响应.故在青藏高原的动力学模型构建中占有着重要地位.由于西构造结区深部结构和动力学研究工作尚少,而又因自然条件和政治环境的制约,故该区深部物质与能量的交换和其物理-力学过程等核心问题尚有待于加强研究和探索.     

LIMITATION OF CURRENT TECTONIC DEFORMATION MODES IN THE WESTERN MARGIN OF ORDOS BASED ON SEISMIC ACTIVITY CHARACTERISTICS

Due to the interaction between the Tibetan plateau, the Alxa block and the Ordos block, the western margin of Ordos(33.5°~39°N, 104°~108°E)has complex tectonic features and deformation patterns with strong tectonic activities and active faults. Active faults with different strikes and characteristics have been developed, including the Haiyuan Fault, the Xiangshan-Tianjingshan Fault, the Liupanshan Fault, the Yunwushan Fault, the Yantongshan Fault, the eastern Luoshan Fault, the Sanguankou-Niushoushan Fault, the Yellow River Fault, the west Qinling Fault, and the Xiaoguanshan Fault. In this study, 7 845 earthquakes(M≥1.0)from January 1st, 1990 to June 30th, 2018 were relocated using the double-difference location algorithm, and finally, we got valid locations for 4 417 earthquakes. Meanwhile, we determined focal mechanism solutions for 54 earthquakes(M≥3.5)from February 28th, 2009 to September 2nd, 2017 by the Cut and Paste(CAP)method and collected 15 focal mechanism solutions from previous studies. The spatial distribution law of the earthquake, the main active fault geometry and the regional tectonic stress field characteristics are studied comprehensively. We found that the earthquakes are more spatially concentrated after the relocation, and the epicenters of larger earthquakes(M≥3.5) are located at the edge of main active faults. The average hypocenter depth is about 8km and the seismogenic layer ranges from 0 to 20km. The spatial distributions and geometry structures of the faults and the regional deformation feature are clearly mapped with the relocated earthquakes and vertical profiles. The complex focal mechanism solutions indicate that the arc-shaped tectonic belt consisting of Haiyuan Fault, Xiangshan-Tianjingshan Fault and Yantongshan Fault is dominated by compression and torsion; the Yellow River Fault is mainly by stretching; the west Qinling Fault is characterized by shear and compression. The structural properties of the fault structure are dominated by strike-slip and thrust, with a larger strike-slip component. The near-north-south Yellow River Fault is characterized by high angle NW dipping and normal fault motion. Based on small earthquake relocation and focal mechanism solution results, and in combination with published active structures and geophysical data in the study area, it is confirmed that the western margin of Ordos is affected by the three blocks of the Tibetan plateau, the Alax and the Ordos, presenting different tectonic deformation modes, and there are also obvious differences in motion among the secondary blocks between the active faults. The area south of the Xiangshan-Tianjingshan Fault has moved southeastward since the early Quaternary; the Yinchuan Basin and the block in the eastern margin of the Yellow River Fault move toward the SE direction.     

青藏高原上地幔速度结构及其动力学性质

利用地震层析成像结果分析了中国西部地区的上地幔速度结构,发现青藏高原北部至东南边缘上地幔顶部速度普遍偏低;随着深度的增加,低速区主要分布在羌塘、松潘—甘孜和云南西部地区,而印度大陆、塔里木、柴达木、鄂尔多斯和四川盆地均显示出较高的速度.上述速度分布与青藏高原及周边地区的岩石层结构和深部动力性质密切相关:其中羌塘地区的低速异常反映了青藏北部的地幔上涌和局部熔融,起因于印度大陆岩石层的向北俯冲;松潘—甘孜地区的低速异常与青藏东部的深层物质流动及四川盆地刚性岩石层的阻挡有关;而滇西地区的低速异常可能受到印缅块体向东俯冲作用的影响.以上三个区域构成青藏高原和周边地区的主要地幔异常区.相比之下,印度大陆、塔里木、柴达木、鄂尔多斯和四川盆地的高速异常反映了大陆构造稳定地区的岩石层地幔特点.根据速度变化推测,地幔上涌和韧性变形并非贯穿整个青藏高原,而是主要集中在羌塘、松潘—甘孜和滇西地区,上述构造效应不仅导致岩石层厚度减薄且引发了火山和岩浆活动.     

青藏高原东缘的地壳流及动力过程

黏滞性地壳流对地壳及上地幔变形作用及动力机制,是大陆新生代造山带的一个重要研究内容.青藏高原中下地壳存在部分熔融或含水物质的黏滞性流体,已为一系列地球物理及岩石学研究所证实.为研究青藏高原东缘地壳流的动力作用,本文用密集的被动源宽频带地震台的观测数据,反演了地壳上地幔精细速度结构和泊松比.研究表明,川西及滇西北高原的中地壳内普遍存在低速层,而高泊松比的地壳只分布在川西北地区.位于中地壳的黏滞性地壳流从青藏高原腹地羌塘高原流出,自北西向南东流入青藏高原东缘.这些黏滞性地壳流带动了上地壳块体水平移动,当它们受到刚强的四川盆地及华南地块阻挡时将发生分层作用,地壳流将分为二或更多分支不同方向的分流,向上的一支地壳流将对上地壳产生挤压,引起地面隆升,向下的一支地壳流将使莫霍面下沉加厚下地壳·黏滞性地壳流的运动在地壳中产生应变破裂发生强烈地震活动,地震的空间分布与震源机制也受到地壳流动力作用控制.     

青藏高原东北部多尺度重力场及其地球动力学意义

由于受到SN和EW双向挤压作用的综合影响,处于年轻高原与古老地块交接区的青藏高原东北部岩石圈强烈变形,构造活动十分活跃.为了整体性地了解青藏高原东北部重力场和深部动力学机制,本文基于EGM2008全球重力场模型数据,利用小波多尺度分析获得不同尺度的重力异常信息;同时反演了研究区的地壳厚度;通过构建穿越龙门山造山带和西秦岭造山带两条剖面的岩石圈密度结构模型,分析了地壳上地幔内不同介质的分布特征.研究结果表明,青藏高原东北部岩石圈显示十分复杂的塑性体的特征,其重力异常走向多以EW或SSE为主,反映了高原岩石圈物质向东运移的趋势;地壳厚度由西向东逐渐减薄,边缘造山带深部并没有发现"山根"痕迹,结合该地区低重力的特征,推测西秦岭—松潘构造结岩石圈发生过大规模地幔流底侵作用;地幔流上涌的动力可能来源于印度板块向欧亚大陆板块俯冲,激发了地幔流体侧向移动,在扬子地台和华北地台附近受到坚硬岩石圈的阻挡而被迫上移,并因此造成龙门山与西秦岭的隆升.     

滇西地区壳幔解耦与腾冲火山区岩浆活动的深部构造研究

根据青藏东部边缘的深部地球物理资料,分析了滇西地区壳幔耦合和腾冲火山区岩浆活动的深部构造特征,确认了地幔各向异性与上地幔速度结构(包括P波速度和S波速度)的内在联系,指出产生这一结果的原因与以腾冲火山区为中心的地幔热物质上涌有关:上地幔顶部平均温度升高导致介质强度降低,在印支块体的侧向挤压或印缅块体的向东俯冲作用下发生韧性变形,造成滇西地区地幔各向异性的快波方向与青藏东部地壳块体的旋转方向不一致.此外,鉴于中下地壳低速层的横向非均匀性,估计韧性流动并非贯通青藏高原的东部边缘,而是被不同的构造块体和边界断裂限定在局部地区.总体而言,滇西地区下地壳的地震波速度和电阻率偏低,具备发生韧性变形的构造条件.作为地壳和上地幔之间的解耦层,它使得青藏东部地壳块体旋转产生的构造应力未能传输至上地幔.腾冲火山区的地壳结构与不同时期的岩浆活动有关,火山区东侧的高速结构代表了上新世时期火山通道内冷凝固结的岩浆侵入体或难以挥发的高密度残留物质,火山区西侧的低速结构反映了更新世以来持续至今的岩浆活动,壳内岩浆源主要分布在10～20km的深度范围内,横向尺度约为15～20km,有可能通过地壳深部的断裂与上地幔岩浆源区相连,估计腾冲火山区下方的岩浆活动将持续进行.     

NUMERICAL SIMULATION OF LONG-TERM DEFORMATION OF TIBETAN PLATEAU AND SURROUNDING AREA

The subduction of the Indian plate underneath Eurasian plate results not only in deformation and movement of the elastic upper crust, but also flow of the ductile lower crust in the high temperature and high pressure which drags the brittle upper crust to move at the same time. These two actions work together producing the present movement and deformation field in Tibetan plateau. The dynamics progress has been verified by GPS observation data. Therefore, in a two-dimension plain model, only the elastic deformation with the boundary action at the upper crust cannot explain the deformation well, the drag force acted on the base of upper crust by the drag of ductile flow of the lower crust also need to be considered. However, it's hard to figure out the magnitude and direction of the drag force. Thus, we established a two-dimension plain elastic finite element model, with the equivalent-body force approach to simulate the drag force. With the internal GPS observation data of Tibetan plateau as constraint condition, we calculated inversely the drag force of key nodes in the model with trial method, and the other nodes in the model with bilinear interpolation method. Finally, we got the drag forces(nodal forces, unit:N) caused by the difference flow of ductile lower crust dragging the brittle upper crust, which are distributed mainly in the region of 86°~100°E and 26°~32°N, the direction is east and south, and the maximum reaches to 1e8N; in some areas in the western part of the study region at 31°~36°N and 76°~80°E, the direction is west, and the maximum reaches to 1e7N. All these work provides a new thought for further research on long-term dynamic mechanism of surface deformation in Tibetan plateau and its surrounding area.