基于GNSS的青藏高原东北缘地壳运动场及强震趋势研究

利用“中国大陆构造环境监测网络”GNSS数据研究1998—2018年青藏高原东北缘排除同震影响等干扰后的速度场、主应变率场、最大剪切应变率场、面应变场等的变化,活动断裂滑动速率变化、跨活动断裂基线变化等。将研究区域内的二级块体再分区,获得各次级块体内部的应变率变化;获取研究区域地壳运动场的趋势性、动态特征。研究结果显示,阿尔金断裂带中东段、祁连块体和柴达木块体交界、巴颜喀拉块体与羌塘块体交界、祁连块体南边界中段、海原—六盘山断裂带和西秦岭北缘断裂带西段的逆冲运动,祁连块体北边界西段、庄浪河断裂的左旋走滑运动,祁连块体北边界东段、西秦岭北缘断裂带东段的左旋逆走滑运动,都属于造成一定程度地壳变形的持续性局部应变增强活动。阿尔金断裂带东段、东昆仑断裂带中西段、祁连块体北边界、庄浪河断裂北段、海原断裂南段、六盘山断裂北段、西秦岭北缘断裂带东段可能存在闭锁,未来十年可能发生M_S6.0以上地震。     

青藏高原东北缘活动断裂剪切模量及应力状态数值模拟

作为控制断层两盘相对运动的重要因素,断裂带介质力学性能与断层面上的滑动速率及应力状态、区域地壳运动速度场等密切相关.受印度板块北东向推挤以及阿拉善地块和鄂尔多斯地块的阻挡作用,青藏高原东北缘构造变形复杂.本文在综合区域动力学环境、活动断裂空间展布以及下地壳黏滞性结构的基础上构建了青藏高原东北缘三维有限元动力学模型;以GPS速度场为约束模拟研究了断层剪切力学性能对区域地壳运动速度场图像的控制作用,进而在最优模型基础上分析了当前青藏高原东北缘不同断裂的应力状态.结果显示：阿尔金断裂东段和广义海原断裂对区域地壳运动速度场控制作用强烈,但二者剪切力学性能相反,阿尔金断裂东段断层剪切模量与周边地壳介质相当,而广义海原断裂断层剪切模量可低至周边地壳介质剪切模量的1/10000;六盘山断裂和西秦岭北缘断裂对区域地壳运动速度场的控制作用较弱,模拟结果显示二者均具有较强的剪切力学性能.基于最佳模型的应力状态分析指出：阿尔金断裂东段,广义海原断裂西段的木里—江仓断裂、中段的金强河—毛毛山—老虎山断裂、东段的六盘山断裂,以及西秦岭北缘断裂中西段当前应力率水平较高,且与前人给出的青藏高原东北缘高闭锁区域吻合.动力学上的高应力率与运动学上的强闭锁良好吻合,预示着这些断裂是地震危险分析值得关注的区域.     

青藏高原东北部现今地壳形变

<正>地壳空间变形模式是探究青藏高原隆升的关键,青藏高原东北部作为高原隆升和扩展的前缘,是由不同构造背景且被活动断裂分割的地块组成。因此,有必要研究不同地块在应变从高原内部向外迁移、分配中的各自行为,及地壳垂直变形对侧向扩展的响应。对GPS水平速度场和应变率场的分析表明:研究区应变向东的耗散主要沿东昆仑和海原断裂。而其应     

GPS观测研究现今青藏高原地壳形变机制——来自阿尔金断裂三维运动场及高原地壳减薄的证据

自50~55 Ma以来,印度次大陆向北与欧亚大陆碰撞后形成喜马拉雅—青藏高原造山带,碰撞导致地壳增厚致使高原大幅隆升,改变了亚洲大陆岩石圈的构造格局,也对东亚地区的气候和环境产生了巨大影响。阿尔金断裂作为青藏高原北缘的主控边界断裂,其运动学性质在20世纪70年代备受关注,不同量级的滑动速率引出了块体运动与东向逃逸和连续变形与地壳增厚两种端元模型。约10~15 Ma以来,在青藏高原南部与北部出现地堑与裂谷,为高原东西向拉张运动提供了证据,表明青藏高原开始经历地壳减薄过程。青藏高原形成以来形变场经历怎样变化,长时间尺度的地质学构造过程与现今GPS观测是否能够统一?10~15 Ma以来青藏高原地壳减薄过程造成高原高程怎样的变化?青藏高原北缘,尤其是跨阿尔金断裂具有怎样的现今三维地壳变形场,地壳应变是如何在北阿尔金断裂、祁漫塔格断裂和阿尔金断裂之间分配的?青藏高原北缘与塔里木盆地具有怎样的力学性质,对跨阿尔金断裂构造形变场造成怎样的影响?最后,GPS观测得到的现今地表形变场能够对青藏高原形变模式的争论作出何种解答?上述科学问题的解答,对于研究青藏高原隆升与变形过程具有十分重要的意义。本研究分为两部分。第一部分是青藏高原北缘三维震间运动场的观测与研究。在青藏高原北缘跨阿尔金断裂中段自建9个GPS连续台站并开展观测,根据区域研究特点设计无人值守的观测台站,具有低成本投入、高质量观测的特点。上述连续GPS台站的建立填补了青藏高原北缘,尤其是在阿尔金无人区地壳形变观测研究的空白,积累了宝贵的连续GPS数据;截止2015年7月,共有4年的连续GPS观测。数据分析结果证明,设计建站方法行之有效,GPS台站稳定、观测数据质量稳定、数据连续性稳定。结合使用中国大陆构造环境监测网络在研究区及邻域GPS连续台站数据作位置时间序列与速度场解算,获得青藏高原北缘地区跨阿尔金断裂中段现今三维形变场。使用三维线弹性后向滑移(backslip)块体运动模型,反演塔里木块体、北阿尔金块体、柴达木块体和祁漫塔格块体的三维块体运动。结果表明,北阿尔金山相对于塔里木盆地有(1.32±0.2)mm/a的抬升速率,相对于柴达木盆地具有(0.73±0.3)mm/a的抬升速率,可解释为北阿尔金块体存在显著的造山过程;阿尔金断裂有(8.21±0.60)mm/a的左旋走滑速率、(0.66±0.60)mm/a的缩短速率;祁漫塔格断裂有(0.53±0.60)mm/a的左旋走滑速率、(1.53±0.60)mm/a的缩短速率;北阿尔金断裂有(0.87±0.60)mm/a的左旋速率、(0.69±0.60)mm/a的缩短速率。同时,阿尔金断裂中、西两段滑动速率基本一致,约为8.0~10.0mm/a。定量研究结果支持连续形变与地壳增厚模型,表明相对塔里木块体,青藏高原北缘地区正在抬升、增厚,以北阿尔金山地区最为明显,抬升速率约达1.3mm/a。跨青藏高原北缘的阿尔金断裂、北阿尔金断裂和祁漫塔格断裂近200km的宽泛变形带内,南北向地壳缩短并不明显,缩短量仅约为2.9mm,且近一半缩短量发生在祁漫塔格山南侧。GPS观测阿尔金断裂车尔臣河段(~86°E)剖面表明,断裂两侧存在非对称变形特征。本文采用非对称变形模型反演GPS速度剖面数据,获得断裂两侧塔里木盆地和青藏高原北部的地壳介质剪切模量差异。结果显示,塔里木盆地地壳介质剪切模量约为青藏高原北部剪切模量1.53倍,相应S波波速比值为1.24,与Yang等人得到的地壳和上地幔三维VSV模型结果一致。地震学研究结果认为,青藏高原北部与东部地区在中地壳存在低速层,局部区域可能发生部分熔融;Hacker等进一步确认羌塘地块中地壳到深部地壳存在熔融现象。本文的研究运用了与地震学完全不同的资料,通过大地测量方法推导青藏高原北部与塔里木盆地的地壳介质力学性质差异,得到与地震学研究得到的S波波速比及其构造物理学解释相当一致的结果。成果为青藏高原力学演化模型提供新的约束。本论文第二部分内容是使用覆盖青藏高原及周边的GPS速度场,计算青藏高原内部应变率场。GPS观测速度场不仅显示了南东东-北西西向的地壳拉张过程,也揭示了青藏高原内部更加重要的地壳减薄过程。结果显示,青藏高原北部和南部的垂向应变率(减薄应变率)分别为(8.9±0.8)nanostrain/a和(7.4±1.2)nanostrain/a,青藏高原西南部的垂向应变率为(12.0±3.2)nanostrain/a,表明青藏高原内部大尺度范围应变率测量结果的一致性。并且青藏高原内部的拉张应变率观测也相当一致,青藏高原北部,沿着N114±1°E主应变方向的拉张应变率为(21.9±0.4)nanostrain/a;高原南部沿着N93±1°E主应变方向的拉张应变率为(16.9±0.2)nanostrain/a;高原西南部沿着N74±3°E主应变方向的拉张应变率为(22.2±1.8)nanostrain/a。如果地壳减薄开始于10~15 Ma,并且现今观测得到应变率适用于整个时间跨度,那么地壳累积减薄5.5~8.5km。应用Airy地壳均衡理论,青藏高原的平均高程将下降~1km。青藏高原北部、南部和西南部相似的垂向应变速率也表明,在3个区域的地壳拉张、正断裂运动和地壳减薄过程由相同的物理机制所支配。综合上述两部分研究成果,发现青藏高原现今垂向运动在高原内部和边缘地区存在很大差别。高原内部地区正在经历地壳减薄,而高原边缘地区正在经历不同程度的增厚与隆升。青藏高原北缘地区的垂向应变率约5~20nanostrain/a,如果考虑重力均衡作用,对应的垂向隆升速率在0.04~0.14mm/a左右。但是,对于局部地区如北阿尔金块体,其底部受到塔里木盆地南缘下插挠曲板块的支持,在没有重力均衡情况下,垂向隆升速率可能达到1mm/a。喜马拉雅地区呈现不同水平的垂向形变,垂向应变强烈(约10~80nanostrain/a),山脉底部受到印度下插板片的支持,无法通过重力均衡假定由垂向应变率估计隆升速率。但由GPS与水准数据约束的俯冲板片模型推测山脉隆升速率达到约7mm/a。而对于祁连山地区,GPS应变率推测得到垂向应变率约20~40nanostrain/a,应用地壳均衡理论,平均隆升速率为0.15~0.3mm/a;而由于逆冲推覆构造与褶皱变形带的存在,中下地壳有可能仍存在弹性变形,不能实现完全重力均衡,实际隆升速率有可能高于这一估计。本文研究给出青藏高原不同地区三维形变场与形变速率的定量估计,是对连续形变与地壳增厚形变模型的重要修正。结果并不支持块体运动与东向逃逸模型,并认为高原南北双向俯冲模型中的塔里木块体南向俯冲几乎不存在。     

GPS观测结果反映的尼泊尔M_w7.8地震孕震特征

针对2015年4月25日尼泊尔Mw7.8地震的孕震特征,本文首先对覆盖尼泊尔及周边地区的5套GPS水平速度场结果进行了融合,得到了近似统一参考框架下的速度场结果;在此基础上通过对此次地震震源区及周边地区的速度场、应变率场、基线时间序列分析,识别了震前变形特征.GPS应变率场结果显示,喜马拉雅主边界断裂存在大范围挤压应变积累,震源区处于近南北向应变积累高值过渡区.跨喜马拉雅构造带的GPS基线时间序列结果表现为持续缩短现象,表明印度板块与欧亚板块之间的持续挤压变形特征,2012年以来的缩短增强现象反映了印度板块对青藏块体的推挤增强作用明显.距离震中较近的西藏南部GPS同震位移结果以南向运动为主且指向震中,反映了青藏高原存在逆冲应变释放现象.综合此次尼泊尔地震前变形和同震应变释放特征,认为此次地震的孕震区域和同震应变释放区域均较大,将会对青藏高原的地壳变形与强震孕育产生深远影响.     

青藏高原东北缘与主要断裂带相关的构造应力率和应变率场的数值模拟

新生代以来,青藏高原东北缘发育了一系列大型走滑和逆冲断裂,构造活动强烈,大震发生频繁.这些主要断裂如何影响区域构造应力、应变场的分配及地震危险性一直没得到充分探讨.本文建立了三维黏弹性有限元数值模型,计算了该地区在现今构造加载作用下的地壳应力率场和应变率场.结果显示,受阿拉善块体和鄂尔多斯块体的阻挡,青藏高原东北缘水平主压应力率及水平主压应变率的方向整体呈顺时针旋转,这种旋转在东南方向更为显著.区域整体以北东-南西向的水平主压应力率为主,并伴随有北西-南东向水平引张应力率.模型0km深度处水平主压应变率的最大值在东昆仑断裂附近,约为4×10-8a-1;在稳定的阿拉善块体和鄂尔多斯块体内部较小,约为1×10-8a-1.计算得到的区域应力状态和实际震源机制解具有较好的一致性,表明青藏高原东北缘以走滑和逆冲型应力场为主.阿尔金断裂带、东昆仑断裂以及海原断裂均表现为最大剪应变率向东逐渐减小,并且最大剪应变率值的减小被其端部的造山隆起和地壳缩短所吸收.在阿尔金断裂的西段、东昆仑断裂的西段到中段,以及海原断裂的西段到中段,未来有较高的地震危险性.     

巴颜喀拉块体北东地区现今水平运动与变形

本文利用GPS数据研究了巴颜喀拉块体北东地区现今水平运动与变形特征。 在球坐标系中解算了各应变分量, 分析了应变率场的空间分布特征, 并与地球物理学和地震地质学研究结果进行了综合对比分析。 最新的GPS速度场结果表明, 巴颜喀拉块体北东地区与高原整体运动性质一样具有顺时针向南东方向旋转的特征, 自西向东和北东方向测站水平运动速度呈现明显的衰减特征。 应变场结果显示, 研究区以北东向的主压应变为主, 伴随着近北西向的张性应变。 应变较强的区域主要分布在活动块体的边界断裂东昆仑断裂带的东段塔藏段和龙门山断裂带上。 东昆仑断裂带东段塔藏段的主压应变明显, 结合地震地质和活动构造资料, 认为东昆仑断裂带东段塔藏段的运动性质自西向东发生了改变, 水平滑动速率逐渐减小, 垂向运动逐渐增强。 研究区GPS速度场和应变场的这一变形特征表明, 青藏高原内部的块体运动特征较为明显, 变形主要集中在作为活动块体边界的活动断裂带上, 边界断裂带的运动特征在调节活动块体间的相互运动中起着重要作用。     

联合GPS和InSAR研究海原-六盘山断裂现今的地壳变形特征

海原—六盘山构造区为青藏高原东北部构造变形最为显著的区域之一,历史强震活动频繁,是研究青藏高原NE向扩展的重要窗口和地震孕育过程的理想场所。文中处理了跨海原-六盘山断裂2014—2020年期间2个轨道的时序Sentinel-1A/B SAR数据,获得了该区域InSAR视线向现今的地壳形变场。融合公开发表的近十多年时间尺度的水平GPS地壳运动速度场,获得了研究区高密度地壳水平形变场。对比GPS、水准和InSAR观测结果,以及GPS-InSAR融合的高密度水平形变场,分析讨论了该区域的地壳形变、应变场特征及其与构造之间的对应关系。主要结论如下：1)GPS和InSAR观测表明,1920年海原8.5级大地震的震后黏弹性松弛效应在海原断裂南侧至今仍较为明显;2)GPS-InSAR高分辨率水平形变场表明,狭义海原断裂左旋滑动速率的递减主要发生在中东段,而中西段递减并不显著,可能与海原断裂向六盘山断裂之间由左旋走滑向逆冲推覆构造转换有关;3)六盘山断裂中—南段的地壳垂直形变和水平形变场特征均显示,该段断裂可能处于强震孕育的中晚期,根据反演得到的断层运动参数和地震地质资料,估算六盘山断裂中—南段发生强...     

青藏高原东北缘主要断裂滑动速率及其动力学意义

本文利用青藏高原东北缘地区1991—2015年的GPS速度场资料, 基于弹性球面块体模型获得了区域活动断裂的滑动速率, 并讨论了断裂滑动速率分配的动力学意义。 反演结果表明, 青藏高原东北缘地区主要块体以北东向并兼顺时针旋转运动为主; 区域断裂平均闭锁深度为17 km; 另外, 各主要断裂滑动速率也不尽相同。 其中, 阿尔金断裂、 东昆仑断裂左旋走滑速率为10～12 mm/a, 祁连—海原断裂左旋走滑速率为3～5 mm/a, 鄂拉山断裂、 拉脊山断裂右旋走滑速率为1～3 mm/a。 阿尔金断裂、 祁连—海原断裂、 东昆仑断裂的走滑速率被其端部的山脉隆起和逆冲断裂所吸收和转换, 鄂拉山断裂和拉脊山断裂则起到了调节块体间运动平衡的作用。     

中国大陆西部GNSS变形特征及其与M≥7.0强震孕育的关系

基于中国大陆地壳运动观测网络GNSS速度场结果,通过最小二乘配置建模、速度残差检验、应变率场分析等,研究了西部地区地壳变形特征及其与M≥7.0强震孕育的关系。主要认识包括：①GNSS应变率场结果显示青藏高原西部地区（92.5°E以西）呈现明显的EW向拉张变形特征,青藏高原东部（92.5°E～100°E）则表现为显著的EW向挤压应变积累。②GNSS旋转率场显示中国大陆西部呈现由南向北逐渐衰减的交替旋转变形现象,藏南地区为大范围顺时针旋转变形特征,旋转率极值达4.5×10^-8rad/a;柴达木地块表现为逆时针旋转变形特征,极值达-1.0×10^-8rad/a;塔里木地块表现为顺时针旋转变形特征,极值达1.0×10^-8rad/a。③2001～2018年中国大陆西部的7次M≥7.0强震均发生在与其构造背景一致的应变积累高值区边缘,呈现一定的孕震晚期特征。因此,最大剪应变率高值区边缘和大型走滑断裂交界部位、张应变率高值区边缘与大型正断层的交界部位、压应变率高值区边缘与大型逆断层交界部位是未来强震需要关注的地点。④2001年昆仑山口西地震导致了青藏高原东部地区较大尺度EW向挤压应变增强现象,在一定程度上反映了巴颜喀拉地块东向运动增强引起的变形调整过程,有利于汶川地震、芦山地震的孕育发生。     

青藏高原东南缘南段现今变形特征研究

本文以青藏高原东南缘南段1999—2017年的GPS速度场为主,结合小震分布、历史地震和活断层探测等资料,首先,基于Okada断层位错模型反演了研究区域主要活断层的滑动速率;其次,以断层滑动速率和GPS速度场观测资料作为约束,利用DEFNODE负位错方法反演了研究区域的块体内部变形及主要活断层的闭锁程度和滑动亏损;最后,计算研究区域现今应变率场,并结合Pms和XKS剪切波分裂结果,探讨分析了青藏高原东南缘的动力学特征.研究结果表明:(1)红河断裂带现今滑动速率明显低于南华—楚雄—建水断裂和无量山断裂;(2)红河断裂带的元江—元阳段、鹤庆—洱源段和小江断裂带北段处于强闭锁状态,南华—楚雄—建水断裂带和无量山断裂带中—北段的闭锁程度强于南段;(3)青藏高原东南缘南段现今地壳变形表现为近E-W向的拉张和近N-S向的挤压,最大剪切方向与Pms和XKS剪切波分裂的快波方向呈一定角度,表明地壳与地幔处于完全解耦状态,而中-下地壳低速层可能是壳幔解耦的主要原因之一;(4)青藏高原东南缘的整体变形受控于印度板块的推挤、印缅俯冲带的深源俯冲以及缅甸微板块与巽他板块的后撤/回退的共同作用.     

巴颜喀拉块体东部活动块体的划分、形变特征及构造意义

基于活动块体的基本概念,综合对研究区内活动断裂带空间展布、地震活动性等资料的分析将巴颜喀拉块体东部及邻区划分为巴颜喀拉块体（I）、华南块体（Ⅱ）、川滇块体（Ⅲ）和西秦岭块体（IV）等4个一级块体.利用GPS形变场、地球物理场等资料结合F检验法,将巴颜喀拉块体划分为阿坝（I₁）、马尔康（I₂）和龙门山（I₃）3个次级块体,将西秦岭块体划分为岷县（IV₁）和礼县（IV₂） 2个次级块体.利用分布在各个块体内部的GPS测站,计算各活动块体及块体边界断裂带的运动变形特征.结果表明：各活动块体的整体运动包括平移和旋转运动;东昆仑断裂带、甘孜—玉树断裂带和鲜水河断裂带的滑动速率明显高于龙门山断裂带的滑动速率;巴颜喀拉块体东部走向北西或北西西的边界断裂表现出左旋拉张的特性;走向北东的边界断裂带,除成县—太白断裂带外,均表现出右旋走滑兼挤压的活动特征.巴颜喀拉块体的东向运动存在自西向东的速度衰减,衰减主要被龙日坝断裂带和岷江断裂带分解吸收,其中龙日坝断裂带的水平右旋分解非常明显,约为~4.8±1.6 mm/a,岷江断裂带的水平分解较弱.龙门山断裂带被马尔康、龙门山和岷县等次级块体分成南、中、北三段,龙门山断裂带中段上的主压应变率要明显小于龙门山断裂带南段上的应变率,其北西侧变形幅度从远离断裂带较大到靠近断裂带逐渐减小,表明其在震前已经积累了较高的应变能,有利于发生破裂滑动.汶川地震后,地表破裂带和余震分布揭示的断裂带运动性质自南西向北东由以逆冲运动为主,逐渐转为逆冲兼走滑的特征可能与龙门山断裂带中段所受主压应力方向自南西向北东的变化有关.马尔康、龙门山和岷县3个次级块体与华南块体之间较低的相对运动速度以及龙门山断裂带低应变率、强闭锁的特征都决定了汶川地震前龙门山断裂带低滑动速率的运动特征.     

三维数值模拟在华北地区现今构造变形分析中的应用

基于并行版ANSYS数值模拟软件平台,根据华北地区活动地块划分及活动断裂分布,结合GPS资料等,确定模型的几何边界范围是99.8°~121.4°E,27.9°~42.3°N,模型的边界范围包括华北活动地块区的绝大部分以及周边青藏、西域、南华和东北亚地块的部分地区。模型单元平均边长为25km,模型共划分为单元416 582个,节点582 392个。主要研究成果为:1)地壳运动速度模拟与分析对比:华北地区地壳运动速率的整体分布呈现由东向西逐渐减小、自北向南逐渐增大的特点,数值模拟结果与GPS观测运动速度大部分较吻合;2)断层运动性质模拟与分析对比:华北地区有限元模型考虑了华北地区已知所有晚更新世以来的活动断裂,模拟结果所得的断层运动性质与活动断裂地质调查的运动性质符合得较好;3)应变场模拟与分析对比:通过数值模拟,计算了1999—2004年和2004—2007年华北地区地表第1主应变和第3主应变(最大压应变)的大小和方向,与前人利用震源机制结果、GPS测量结果等资料反演得到的华北地区主应变场较一致。     

PROGRESS IN APPLICATION OF GNSS TO DIVISION OF ACTIVE TECTONIC BLOCKS IN CONTINENTAL CHINA

Chinese scientists proposed that large earthquakes that occurred in mainland China are controlled by the movement and deformation of active tectonic blocks. This scientific hypothesis explains zoned phenomenon of seismicity in space. The active tectonic blocks are intense active terranes formed in late Cenozoic and late Quaternary, and the tectonic activity of block boundaries is the intensest. Global Navigation Satellite System(GNSS)has advantages of high spatio-temporal resolution, broad coverage, and high accuracy, and is utilized to monitor contemporary crustal deformation. High accuracy and resolution of GNSS velocity field within mainland China and vicinities provided by previous studies clearly demonstrate that different active tectonic blocks behave as different patterns of movement and deformation, and block interaction boundaries have intense tectonic deformation. The paper firstly introduces the GPS networks operated by the Crustal Movement Observation Network of China(CMONOC)since 1999, and GNSS data processing methods, including GAMIT, BERNESE and GIPSY/OASIS, and discusses the advantages of using South China block as a regional reference frame for GNSS velocity field, then proposes three strategies of block division, F-test, quasi-accurate detection(QUAD), and clustering analysis. Furthermore, we introduce rigid and non-rigid block motions. Rigid block motion can be denoted by translation and rotation, while non-rigid block motion can be described by rigid motion and internal strain deformation. Internal strain deformation can be divided into uniform and linear strains. We also review the usage of F-test to distinguish whether the block acts as rigid deformation or not. In addition, combining with recent GNSS velocity results, we elaborate the characteristics of present movement of rigid block, such as the South China, Tarim, Ordos, Alashan, and Northeast China, and that of non-rigid block, such as the Tibetan plateau, Tian Shan, and North China plain. Especially, the Tibetan plateau and Tian Shan seem to deform continuously with significant internal deformation. In order to enrich and perfect the active tectonic block hypothesis, we should carefully design dense GNSS networks in inner blocks and block boundaries, optimize utilizing other space geodesy technologies such as InSAR, and strengthen combining study of geodesy, seismogeology and geophysics. Through systematic summary, this paper is very useful to employing GNSS to investigate characteristics of block movement and dynamics of large earthquakes happening in block interaction boundaries.     

GPS观测及断裂晚第四纪滑动速率所反映的青藏高原北部变形

断裂晚第四纪滑动速率及现今GPS观测揭示了青藏高原向北扩展与高原边缘隆升的运动特征.主要断裂晚第四纪滑动速率及跨断裂GPS应变速率的结果表明,青藏高原北部边缘的断裂以低滑动速率（＜10 mm/a）为主,特别是两条边界断裂：阿尔金断裂和海原—祁连山断裂.两条主要边界断裂上的滑动速率分布显示了断裂间滑动速率转换及调整特征.阿尔金断裂自95°E以西的8～12 mm/a稳定滑动速率,向东逐渐降低到最东端的约1～2 mm/a,而海原断裂自哈拉湖一带开始发育后滑动速率为1～2 mm/a,到祁连一带（101°E以东）增大到相对稳定的4～5 mm/a,直到过海原后转向六盘山一带,滑动速率降低到1～3 mm/a,甚至更低.滑动速率的变化及分布特征显示,阿尔金断裂滑动主要是通过祁连山内部隆起及两侧新生代盆地变形引起的缩短来吸收的,海原—祁连山断裂的低滑动速率及沿断裂运动学特征表明断裂尾端的陇西盆地变形及六盘山的隆起是断裂左旋走滑速率的主要吸收方式.这一变形特征表明,青藏高原北部边缘的变形模式是一种分布式的连续变形,变形发生自高原内部,边界断裂的走滑被高原内部变形所吸收.     

Movement of Ordos block and alternation of activity along its boundaries

The GPS data in and around the Ordos block area indicate that the left-lateral slip rate along the northern or southern margin of the Ordos block is about twice or three times as fast as the right-lateral slip rate along the eastern or western margin of the block. However, many researchers stressed the dextral-slip of the eastern or western boundaries of the Ordos block, and suggested that the block as a whole rotated counterclockwise based on the available geological data. Focusing on the inconsistency, we reexamine the late Cenozoic deformation pattern in the Ordos region based on seismicity data and geodesy data (GPS and leveling) around it. The results indicate that the rigid block-like motion appears to be the basic characteristic of the kinematics of the Ordos region, and this motion is absorbed by the displacement of the faults around the block. When the faults along the northern and southern boundaries of the Ordos block are active, its eastern boundary is inactive. However, if the faults along the eastern boundary are active, the northern and southern are inactive. In recent years, the northern and southern boundaries of the Ordos block are in active. But in the long term, the Ordos block is moving southeastward relative to the Alxa and Yinshan blocks because of the strong pushing of the Tibetan Plateau on its southwestern side, and this deformation is accommodated by the counterclockwise rotation of the block itself.     

Movement of Ordos block and alternation of activity along its boundaries

The GPS data in and around the Ordos block area indicate that the left-lateral slip rate along the northern or southern margin of the Ordos block is about twice or three times as fast as the right-lateral slip rate along the eastern or western margin of the block. However, many researchers stressed the dextral-slip of the eastern or western boundaries of the Ordos block, and suggested that the block as a whole rotated counterclockwise based on the available geological data. Focusing on the inconsistency, we reexamine the late Cenozoic deformation pattern in the Ordos region based on seismicity data and geodesy data (GPS and leveling) around it. The results indicate that the rigid block-like motion appears to be the basic characteristic of the kinematics of the Ordos region, and this motion is absorbed by the displacement of the faults around the block. When the faults along the northern and southern boundaries of the Ordos block are active, its eastern boundary is inactive. However, if the faults along the eastern boundary are active, the northern and southern are inactive. In recent years, the northern and southern boundaries of the Ordos block are in active. But in the long term, the Ordos block is moving southeastward relative to the Alxa and Yinshan blocks because of the strong pushing of the Tibetan Plateau on its southwestern side, and this deformation is accommodated by the counterclockwise rotation of the block itself.     

SPATIAL AND TEMPORAL DISTRIBUTION OF SLIP RATE DEFICIT ACROSS HAIYUAN-LIUPAN SHAN FAULT ZONE CONSTRAINED BY GPS DATA

As the northeast boundary of the Tibetan plateau, the Haiyuan-Liupan Shan fault zone has separated the intensely tectonic deformed Tibetan plateau from the stable blocks of Ordos and Alxa since Cenozoic era. It is an active fault with high seismic risk in the west of mainland China. Using geology and geodetic techniques, previous studies have obtained the long-term slip rate across the Haiyuan-Liupan Shan fault zone. However, the detailed locking result and slip rate deficit across this fault zone are scarce. After the 2008 Wenchuan M_S8.0 earthquake, the tectonic stress field of Longmen Shan Fault and its vicinity was changed, which suggests that the crustal movement and potential seismic risk of Haiyuan-Liupan Shan fault zone should be investigated necessarily. Utilizing GPS horizontal velocities observed before and after Wenchuan earthquake(1999~2007 and 2009~2014), the spatial and temporal distributions of locking and slip rate deficit across the Haiyuan-Liupan Shan fault zone are inferred. In our model, we assume that the crustal deformation is caused by block rotation, horizontal strain rate within block and locking on block-bounding faults. The inversion results suggest that the Haiyuan fault zone has a left-lateral strike-slip rate deficit, the northern section of Liupan Shan has a thrust dip-slip rate deficit, while the southern section has a normal dip-slip rate deficit. The locking depths of Maomao Shan and west section of Laohu Shan are 25km during two periods, and the maximum left-lateral slip rate deficit is 6mm/a. The locking depths of east section of Laohu Shan and Haiyuan segment are shallow, and creep slip dominates them presently, which indicates that these sections are in the postseismic relaxation process of the 1920 Haiyuan earthquake. The Liupan Shan Fault has a locking depth of 35km with a maximum dip-slip rate deficit of 2mm/a. After the Wenchuan earthquake, the high slip rate deficit across Liupan Shan Fault migrated from its middle to northern section, and the range decreased, while its southern section had a normal-slip rate deficit. Our results show that the Maomao Shan Fault and west section of Laohu Shan Fault could accumulate strain rapidly and these sections are within the Tianzhu seismic gap. Although the Liupan Shan Fault accumulates strain slowly, a long time has been passed since last large earthquake, and it has accumulated high strain energy possibly. Therefore, the potential seismic risks of these segments are significantly high compared to other segments along the Haiyuan-Liupan Shan fault zone.     

Pattern of latest tectonic motion and its dynamics for active blocks in Sichuan-Yunnan region, China

Based on the concept of "active blocks" and spatial distribution of historical earthquakes with surface ruptures as well as major and subordinate active faults. The Sichuan-Yunnan region can be divided into four first-order blocks. They are the Markam block (I), the Sichuan-Yunnan rhombic block (II), Baoshan-Pu'er block (III), and Mizhina-Ximeng block (IV). Cut by sub-ordinate NE-trending active faults, the Sichuan-Yunnan rhombic block (II) can be further divided into two sub-blocks: the northwestern Sichuan sub-block (II1) and the middle Yunnan sub-block (II2), while the Baoshan- Pu'er block (III) can be further divided into three sub-blocks: Baoshan sub-block (III1), Jinggu sub-block (III2), and Mengla sub-block (III3). A quantitative study of offset landforms is carried out and the basic types of active faults and their long-term slip rates along the major boundaries of active blocks of different orders in the Sichuan-Yunnan region are determined, through slip vector analysis, the motion states of the active blocks are clarified and the deformation coordination on the block margins is discussed. It is suggested that the tectonic motion of the blocks in this region is a complex or superimposition of three basic types of motions: southeastward sliding, rotating on vertical axis, and uplifting. The Markam block (I), the northwestern Sichuan sub-block (II1), and middle Yunnan sub-block (II2) have a southeastward horizontal sliding rate of 1-5 mm/a, clockwise rotating angular rate of 1.4-4(/Ma, and uplifting rate of about 1 mm/a. The Baoshan-Pu'er (III) and Mizhina-Ximeng (IV) blocks have also been extensively clockwise rotated. This pattern of motion is a strain response to the collision between the Indian and Eurasian plates and the localized deformation and differential slip on the block margins associated with the northward motion of the Indian Plate. Because a set of transverse thrusts between the blocks absorbs and transforms some components of eastward or southeastward sliding motion, the eastward escape or extrusion of the Tibetan Plateau is limited as "imbricated thrusting transformation-limited extrusion model".     

青藏高原西北部区域地壳形变、构造地貌与孕震构造模型研究——以2008年与2014年新疆于田7.3级地震为例

为了清晰认识发生于青藏高原西北部2008年与2014年的两次于田MS7.3地震发震构造环境与构造地貌特征,本文利用DEM(数字高程模型)数据分析"喀喇昆仑—西昆仑—康西瓦地区"的地形地貌特征,结合区域活动断裂研究资料、相对于塔里木盆地的两期GPS速度场资料和区域运动学特征等讨论两次MS7.3地震所处的青藏高原西北部区域构造环境和地壳运动学特征,分析喀喇昆仑断裂、阿尔金断裂康西瓦段、龙木错-邦达错断裂及贡嘎错断裂所围限的西昆仑地块的地质构造背景、阿尔金断裂西南端发震断裂活动性及孕震环境等发震构造基本条件;进而利用"地形剖面"方法及断裂分布特征分析震源区的地形地貌特征,给出晚第四纪以来的地貌形态与发震构造的关系,从区域构造地貌学和GPS地壳运动学的角度探讨中上地壳变形特征及孕震过程;最后讨论区域孕震构造、克尔牙张性裂谷演化过程和地球动力学背景等。通过地形剖面及区域地貌综合分析新疆于田2008年MS7.3拉张型发震构造和2014年MS7.3走滑拉张型地震的发震构造特点的区别,认为2014年发生的地震可能与2008年MS7.3地震同震库伦应力变化、触发过程及震后变形过程密切相关,并且青藏高原西北部地区存在明显的东西向拉张性构造单元,可能与青藏高原10~15 Ma以来的地壳减薄过程有关。