-
北京时间2015-04-25T 14:11:26，尼泊尔发生Mw7.8大地震，美国地质调查局(USGS)测定震中位置为(28.147°N，84.708°E)，震源深度为15 km，震源机制结果显示此次地震为逆冲型破裂，地震发生于喜马拉雅构造带中段，由印度板块与欧亚板块俯冲碰撞形成。全球地震台网宽频带波形数据快速反演的地震破裂模型显示，主震断层面最佳模型走向为295°，倾角为10°，释放的地震矩为8.1×1020 N·m(http://earthquake.usgs.gov/earthquakes/eventpage/us20002926)。该地震发生后，2015年4月25日、4月26日与5月12日，又分别发生了尼泊尔Mw6.6、Mw6.7和Mw7.3地震。
印度板块与欧亚板块的俯冲碰撞开始于距今约55~65 Ma，至今仍处于强烈推挤过程中，在长期的陆-陆俯冲碰撞过程中，喜马拉雅构造带上大地震不断发生，也形成了一个由北向南发展的前展式逆断裂系，分别为喜马拉雅主中央断裂(the main central thrust, MCT)、喜马拉雅主边界断裂(the main boundary thrust, MBT)和喜马拉雅主前缘断裂(the main frontal thrust, MFT)3个巨大逆断裂系，其中MBT和MFT晚第四纪以来仍在活动[1]，而且在地壳下部存在喜马拉雅主逆冲断层(the main Himalayan thrust, MHT)[2-4]。目前有关喜马拉雅构造带的相关研究较为丰富，各国学者利用GPS观测资料、水准数据、InSAR数据等对该地区不同区域之间的相对运动速率和MHT闭锁规律进行了研究[3, 5-7]。
虽然目前对断层闭锁和块体相对运动研究较多，但有关断层滑动亏损速率和断层远近场运动速率等的研究较少。地壳变形动态特征是断层深浅部动力作用的表现形式，因此可以利用地表GPS观测结果反演断层深浅部的闭锁程度、滑动亏损分布等特征，进而结合历史地震分布和大地震空区等判断断裂带未来一段时期内的地震危险性。目前，喜马拉雅构造带附近GPS观测资料已较为丰富。本文首先收集了覆盖印度板块、青藏高原及周边地区的6套GPS水平速度场并将其融合为同一参考框架的结果；然后以融合的速度场作为震间地壳变形的约束，通过DEFNODE负位错反演方法和剖面拟合方法分析了MHT的闭锁程度、闭锁断层面宽度、断层滑动亏损速率、远近场运动速率等，并讨论了该逆冲断层的地震危险性。
-
本文反演研究采用DEFNODE负位错反演程序[8]，该方法在国际上已得到了较为广泛的应用[9-17]，程序假定块体内部点的运动为块体旋转、块体内部整体均匀应变、块体边界由于断层闭锁产生的滑动亏损而引起的地表弹性变形之和。若块体内部不存在整体均匀应变，理论表达式为：
$$ \begin{array}{l} {V_i}\left( X \right) = \sum\limits_{b = 1}^B H \left( {X \in {\Delta _b}} \right)\left[ {_R{\mathit{\Omega} _b} \times X} \right] \cdot i - \\ \sum\limits_{k = 1}^F {\sum\limits_{n = 1}^{{N_k}} {\sum\limits_{j = 1}^2 {{\varphi _{nk}}{G_{ij}}} } } (X,{X_{nk}}){[_h}{\mathit{\Omega} _f} \times {X_{nk}}]\cdot j \end{array} $$ (1) 式中, 各参数的定义见文献[8, 18]，其中φnk为断层k上节点n处的闭锁程度，是断层滑动亏损速率与长期滑动速率的比值，若值为0，表明断层完全蠕滑; 若值为1，表明断层完全闭锁。
若块体内部存在整体均匀应变，应变率的计算采用Savage等[19]给出的式(2)。此时模型原理为式(1)加式(2)，其中式(2)表示块体内部均匀应变引起的速度量值大小。
$$ \left[ \begin{array}{l} {V_\lambda }\\ {V_\theta } \end{array} \right] = \left[ {\begin{array}{*{20}{c}} {R\sin {\theta _0}\Delta \lambda }&{R\Delta \theta }&0\\ 0&{R\sin {\theta _0}\Delta \lambda }&{R\Delta \theta } \end{array}} \right]\left[ \begin{array}{l} {\dot \varepsilon _\lambda }\\ {\dot \varepsilon _{\lambda \theta }}\\ {\dot \varepsilon _\theta } \end{array} \right] $$ (2) 式中, 各参数的定义见文献[9]。
利用GPS水平速度场数据反演块体旋转和断层闭锁程度等参数的过程中，参数拟合的不符值可由式(3)表征，其中各参数含义可参考文献[8, 10, 18]:
$$ {\chi _n}^2 = [\sum\limits_{i = 1}^n {{{({r_i}/f{\sigma _i})}^2}} ]/{\rm{dof}} $$ (3) 为得到使χn2达到最小的一组最佳参数值，本反演程序使用网格搜索和模拟退火算法求解，通过不断改变f值大小进行逐步择优，寻求最满意的f值，进而保证χn2≈1，此时模型能够最佳地解释观测数据。
-
根据前人研究结果，本文收集了覆盖印度板块、青藏高原及周边地区的6套GPS速度结果(75°E~110°E, 15°N~55°N)，详细的数据信息见表 1。在应用表 1结果进行地壳变形综合分析之前，应完成不同来源数据参考基准的统一。经过分析发现，表 1中第1套速度场[20]的分布范围较大且包含大部分国际GNSS服务(International GNSS Service, IGS)站点，便于与其他结果融合，选定其为参考速度场。速度场融合的具体方法见文献[25], 但时间尺度上未考虑观测时间点的差异性，因为目前求解速度场一般基于线性模型，认为速度在一段时间内都是线性的，尤其流动观测都是不定期观测，观测时间太短，获得的观测数据太少，只能用线性模型，因此观测时间点的差异性并不会对求解的速度场结果产生较大影响。表 1中公共站速度矢量差的均方差结果显示，融合前后的速度残差较小，均方差均小于2 mm/a，由于目前求解速度场自身的精度为1~2 mm/a左右，而本文融合的速度场精度与求解速度场自身的精度大致相当，表明融合结果满足研究震间变形的要求。
表 1 本文使用的GPS速度场资料
Table 1. GPS Velocity Field Used in This Work
序号 区域 点数 时间/年份 框架 融合时使用的公共站个数 融合后残差均方差/(mm·a-1) 资料来源 1 中国大陆 354 1991~2001 ITRF97 - - 文献[20] 2 青藏高原 45 1991~2001 ITRF97 40 0.49 文献[21] 3 尼泊尔喜马拉雅 97 1995~2006 ITRF2000 33 0.43 文献[6] 4 青藏高原及周边 726 1998~2004 ITRF2000-Eurasia 10 1.48 文献[22] 5 印度及邻区 229 1995~2007 ITRF2000 36 0.35 文献[23] 6 青藏高原 750 1998~2013 ITRF2008-Eurasia 11 1.72 文献[24] 注：ITRF2000-Eurasia为欧亚基准，是ITRF2000的结果去除了相对于欧亚板块的整体运动量。ITRF2008-Eurasia以此类推 图 1给出了融合后的ITRF97参考框架下研究区域的GPS速度场分布。结果显示，印度板块以北东向运动为主，青藏高原内部的地壳运动具有顺时针方向旋转的特点，自南向北逐渐由北东向偏转向东向，且量值有所减小，表明青藏高原内部存在近南北向挤压应变积累且东西向断层存在左旋剪切应变变形背景[25]。
图 1 研究区域GPS速度场结果
Figure 1. The GPS Velocity Field in the Study Region
-
本文研究区域以MHT为界，包含南部的印度板块和北部的青藏地块(如图 1所示，因本文主要研究MHT，青藏地块只包含鲜水河-玉树-玛尔盖查卡断裂带以南的区域)。数据约束方面, 包括320个GPS测点参与反演，其中印度板块46个, 青藏地块274个。根据相关文献[3, 6, 7]的研究结果，断层倾角设为10°。断裂总长度为1 616 km，断层面总宽度为241 km；沿着MHT走向共设8条等深线，深度依次为0.1、6、12、18、24、30、36和42 km，每条等深线上有16个节点。利用DEFNODE程序可反演计算得到每个节点处断层的闭锁程度，进而通过双线性插值方法计算相邻节点之间断层网格(沿走向方向长10 km，深度方向宽6 km)区域的闭锁程度。对断层节点φ(闭锁程度)值加约束为：节点均为独立节点，φ值从地表沿垂直等深线向下单调递减(0.1 km与6 km深度处的φ值相等)。
由于印度板块的内部变形很小[3, 6]，因此本文假定印度板块内部不存在整体均匀应变。Gan等[23]以鲜水河-玉树-玛尔盖查卡断裂带将青藏高原分为南北两大部分，南部区域即为本文的研究区域。由于研究区域距离MHT较近，印度板块的推挤作用、喜马拉雅主逆冲断层的闭锁作用是其运动与变形的动力源，因此将该区域作为一个大的块体来研究，暂不考虑区域内亚块体的内部变形差异。
首先，假定青藏地块内部的变形只是由于MHT的逆冲运动和闭锁作用引起的弹性变形(假定1)，其内部不存在整体均匀应变。对于假定1，通过不断试算得到了最优模型，其中GPS速度场误差权重因子f取为4.2，χn2≈0.991(观测值个数为640，自由度为522)。为了观察模型拟合效果，图 2(a)给出了残差分布情况。结果显示，假定1的最优模型中，印度板块测点拟合较好，除了块体东北角几个测点的残差稍大以外，其他测点的残差基本保持在5 mm/a以内(大部分测点残差处在误差范围之内)。青藏地块断层附近测点拟合很好，残差基本处在测点误差范围之内；而块体东部和西部远离断层、处于青藏地块的东西边缘处的测点拟合较差，这些测点由于距离断层较远，受断层闭锁的约束相对较弱，而且青藏地块南边受到喜马拉雅东、西构造节不同方向的挤压，造成其东、西边缘处分别具有东向、西向的运动特征，使得青藏块体的整体运动特征减弱。对假定1最优模型的残差进行计算，得到了青藏地块的残差应变，其中EEE=6.88×10-9·a-1，ENN=-5.08×10-9·a-1，ENE=-3.66×10-9·a-1，约为北东16°的主压和与之垂直方向的主张应变。这表明青藏地块内部除了由于MHT的逆冲运动和闭锁作用引起的弹性变形外，整体均匀应变的量值也不容忽视[7]。对于假定2，将假定1的最优模型中得到的整体均匀应变值作为约束数据加入到模型中，通过不断试算得到了最优模型，其中GPS速度场误差权重因子f取3.09，χn2≈1.001(自由度为528)。图 2(b)残差分布结果显示，假定2的最优模型中，印度板块的拟合情况与假定1基本一致；青藏地块东部和西部远离断层的测点拟合情况相对于假定1的拟合情况有了明显改善，其误差基本在10 mm/a以内，而且也不具有明显的规律性，断层附近测点拟合情况与假定1基本一致。综上分析，认为MHT附近的变形主要由断层逆冲运动和闭锁作用引起，同时青藏地块内部存在整体均匀应变，因此选择假定2的最优模型作为本研究的最优模型。
图 2 最优模型速度残差
Figure 2. Residual Velocities of the Optimal Model
-
表 2给出了假定2最优模型计算得到的印度板块与青藏地块的欧拉极，结果显示在ITRF97参考框架下，印度板块为逆时针旋转，青藏地块为顺时针旋转。表 2显示，Bettinelli等[6]计算的ITRF2000下印度板块的欧拉极和Ader等[3]计算的ITRF2005下印度板块的欧拉极与本文得到的欧拉极基本一致。
表 2 块体欧拉矢量
Table 2. Euler Vectors of Blocks
-
根据假定2的最优模型，得到MHT闭锁程度的空间分布(图 3)，图 3中震源机制解数据来源于http://www.globalcmt.org，其中黑色五角星为此次震中位置，历史地震破裂区域参考文献[3, 6, 26]的研究结果绘制。最优模型结果表明，整条MHT的闭锁深度基本达到18~24 km，断层面闭锁宽度达到102~136 km；而两次大地震之间的未破裂区域和未发生大地震的不丹地区断层基本为完全闭锁，闭锁深度为42 km、闭锁断层面宽度达到241 km。该结果与文献[3, 6]相比长度更长，除了尼泊尔地区的断层外，还包括了东部的不丹地区和西部尼泊尔与印度接壤地区的断层，且GPS速度场的数据更为丰富；但反演的断层闭锁深度和闭锁断层面宽度在两次大地震之间的未破裂区域结果有所差异。本文除了反演断层闭锁程度，还反演了断层滑动亏损速率的分布。
图 3 喜马拉雅主逆冲断层闭锁程度
Figure 3. Fault Locking of the Main Himalaya Thrust Fault
喜马拉雅构造带是一条强烈活动带，沿这条构造带发生过多次8级及8级以上大地震，图 3中标注了破裂范围较大的3次地震，其中1905年印度坎格拉8.0级地震的破裂区域参考文献[6], 1505年尼泊尔格尔纳利河8.2级地震的破裂区域参考文献[3], 1934年尼泊尔比哈尔邦8.1级地震的破裂区域参考文献[3, 6, 26]。这3次大地震的破裂段落附近，其闭锁深度和闭锁程度较研究区域的其他位置相对较弱。1934年尼泊尔比哈尔邦8.1级地震破裂段落中西段附近，其12 km深度以下的闭锁程度约为0.60；破裂段落中段附近基本为完全闭锁；破裂段落东段闭锁程度约为0.85。1505年尼泊尔格尔纳利河8.2级地震破裂段落附近，其24 km深度以下基本为完全蠕滑状态。1905年印度坎格拉8.0级地震破裂段落附近，其15 km深度以下基本为完全蠕滑状态。
两次大地震破裂区域之间的未破裂段落和未发生大地震的段落闭锁深度更深、闭锁断层面更宽。相关研究表明，此次尼泊尔地震的浅部破裂紧邻1934年8.1级地震的地表破裂[26, 27]，余震分布与1833年7.6级地震的宏观震中基本重合，其破裂填补了前两次地震破裂以西100 km左右的空区，可见此次地震是前两次地震向西继续延伸的结果[26, 28]。2015年尼泊尔Mw7.8大地震发生在强闭锁区域内，一定程度上证明了计算结果的可靠性。
震间期，MHT的闭锁部分吸收了块体俯冲挤压造成的弹性应变，直到有8级甚至更大地震发生时，弹性应变才得到释放[3, 29]，而大地震发生的频率也并不高[30]，Feldl和Bilham[31]曾估计喜马拉雅地区Mw8.6大地震的复发间隔约为500 a。根据USGS[26]的计算结果，此次地震的最大滑移量在4 m左右，结合§3.4中该位置处约17.5 mm/a的垂直断层滑动亏损速率，估算该位置处Mw7.8地震的复发周期约为230 a。目前来看，由于最近500多年只发生了1905年印度坎格拉8.0级、1934年尼泊尔比哈尔邦8.1级和2015年尼泊尔Mw7.8地震，因此整条段落尤其是尼泊尔西部与印度接壤处[30]和可能还没有破裂的不丹地区[26]依然有发生8级大地震的危险。
-
图 4(a)~4(c)分别为MHT总滑动亏损速率、垂直断层滑动亏损速率、平行断层滑动亏损速率。图 4(a)显示，MHT总滑动亏损速率呈现自东向西逐渐减小的特征，最东端滑动亏损速率约为28.2 mm/a，中间完全闭锁区域的滑动亏损速率逐渐过渡，到最西端滑动亏损速率降低到约为7.9 mm/a。
图 4 喜马拉雅主逆冲断层滑动亏损速率
Figure 4. Fault Slip Deficit of the Main Himalaya Thrust Fault
喜马拉雅构造带的地壳缩短吸收了印度板块和欧亚板块俯冲汇聚的很大一部分量值[6]。图 4(b)显示，MHT垂直断层挤压滑动亏损速率也呈现自东向西逐渐减小的特征，最东端挤压滑动亏损速率约为28.0 mm/a，中间完全闭锁区域的滑动亏损速率逐渐过渡，到最西端挤压滑动亏损速率降低到约为6.4 mm/a。其中，尼泊尔境内最东端挤压滑动亏损速率约为22.1 mm/a，最西端挤压滑动亏损速率约为11.0 mm/a，这与Ader[3]、Bettinelli等[6]、Lavé和Avouac[29]以及Ponraj等[4]的研究结果基本一致。MHT全新世时期的平均地质缩短速率为21±1.5 mm/a[29]，因此本文的计算结果与地质结果也较为吻合。
图 4(c)显示，MHT平行断层右旋滑动亏损速率基本上自东向西逐渐增加。其中，最东端右旋滑动亏损速率约为2.0 mm/a，在苏丹与尼泊尔之间的部分区域降低到约0.4 mm/a，进入尼泊尔后逐渐增加，到尼泊尔西部达到最大值7.0 mm/a，在印度约为3.3~5.2 mm/a。
另外，在获得融合的ITRF97参考框架下的GPS速度场后，将其统一扣除相对于欧亚板块的整体运动量，获得欧亚参考基准的GPS速度场; 然后将不同参考基准速度场都进行负位错反演，通过对比两种反演结果发现，其反演的闭锁程度、滑动亏损分布等结果几乎完全一样，这表明DEFNODE负位错反演结果与速度场的参考基准无关，与文献[32]一致，即参考基准的不同并不影响地壳的相对运动，因此也不影响块体之间断层的相对运动与闭锁。
-
图 5(a)~5(c)为MHT 3条剖面的GPS站速度场与最优模型计算结果的拟合情况(分别对应图 1(a)~1(c)剖面)，每条剖面长约1 600 km，宽390 km，基本跨越本文研究断层的西段、中段和东段。图 5(a)显示垂直断层的挤压速率约为14 mm/a，平行断层的右旋滑动速率约为7 mm/a。图 5(b)显示垂直断层的挤压速率约为19 mm/a、平行断层的右旋滑动速率约为5 mm/a。图 5(c)显示垂直断层的挤压速率约为26 mm/a，平行断层的右旋滑动速率约为3 mm/a。剖面拟合结果与断层面上的滑动亏损速率基本一致。
图 5 横跨喜马拉雅主逆冲断层剖面的GPS站速度和最优模型结果拟合
Figure 5. GPS Velocity Profiles Across the Main Himalaya Thrust Fault and the Model Fitting Results
通过分析发现，垂直和平行断裂带走向的挤压缩短变形和剪切变形主要发生在该断裂带的北侧，即青藏块体内部，而断裂带内部几乎无变形，即喜马拉雅主逆冲断裂带内部的GPS站速度与南侧印度板块中的几乎没有差异。另外，MHT的剖面结果反映的变形与2008年四川汶川MS8.0地震前，横跨龙门山断裂带中段的GPS变形图像[9, 18]非常相似，这也反映出喜马拉雅主逆冲断裂带存在很强的闭锁。
-
本文通过对印度板块、青藏高原及周边地区的6套GPS速度场数据进行融合，得到了近似统一参考框架下的速度场结果，然后以融合的速度场作为震间地壳变形的约束，采用DEFNODE负位错反演程序估算了MHT的闭锁程度和滑动亏损空间分布特征；结合剖面结果分析了断层远、近场的运动特征，并讨论了该逆冲断层的地震危险性。
1) 负位错模型反演结果表明，MHT的闭锁深度基本达到18~24 km，断层面闭锁宽度达到102~136 km；3次发生过大地震的破裂段落附近，其闭锁深度较浅、闭锁断层面较窄；两次大地震破裂区域之间的未破裂段落和未发生大地震的段落闭锁深度更深、闭锁断层面更宽。2015年尼泊尔MW7.8大地震就发生在两次大地震破裂段落之间闭锁更深、更强的区域。根据此次地震4 m左右[26]的最大滑移量的计算结果，结合该位置处约17.5 mm/a的垂直断层滑动亏损速率，我们估算该位置处MW7.8地震的复发周期约为230 a。由于最近500多年只发生了1905年印度坎格拉8.0级、1934年尼泊尔比哈尔邦8.1级和2015年尼泊尔MW7.8地震，因此整条段落尤其是尼泊尔西部与印度接壤处和可能还没有破裂的不丹地区依然有发生8级大地震的危险。此外，我们将不同参考基准速度场的反演结果进行对比，其反演的闭锁程度、滑动亏损分布等结果几乎完全一样，这表明反演结果与速度场的参考基准无关。
2) MHT总滑动亏损速率和垂直断层挤压滑动亏损速率呈现自东向西逐渐减小的特征，其中最东端垂直断层挤压滑动亏损速率约为28.0 mm/a，到最西端垂直断层挤压滑动亏损速率约为6.4 mm/a。平行断层右旋滑动亏损速率基本上自东向西逐渐增加，尤其在尼泊尔，右旋滑动亏损速率从约0.4 mm/a增大到约7.0 mm/a。
3) MHT的3条剖面拟合结果与相应位置处断层面上的滑动亏损速率基本一致。垂直和平行断裂带走向的挤压缩短变形和剪切变形主要发生在断裂带北侧的青藏块体内部，而断裂带内部和南侧印度板块几乎无变形，反映出MHT存在很强的闭锁。汶川地震与2015年尼泊尔地震都属于逆冲型地震，其中汶川地震属于仰冲型地震(主动盘为上盘)、尼泊尔地震属于俯冲型地震(主动盘为下盘)，这两次大地震的发生对加深对逆冲型地震的理解和认识起到了重要作用。
Characteristics of Fault Locking and Fault Slip Deficit in the Main Himalaya Thrust Fault
-
摘要: 基于融合的GPS速度场结果，使用DEFNODE负位错反演程序估算了喜马拉雅主逆冲断层（the main Himalayan thrust，MHT）的闭锁程度和滑动亏损空间分布，并结合剖面结果分析了断层远、近场的运动特征。结果表明，MHT的闭锁深度基本达到18~24 km，断层面闭锁宽度达到102~136 km，两次历史大地震破裂区域之间的未破裂段落和未发生大地震的段落闭锁深度更深，闭锁断层面更宽，2015年尼泊尔Mw7.8大地震就发生在两次大地震破裂区域之间的段落；MHT总滑动亏损速率和垂直断层挤压滑动亏损速率自东向西逐渐减小，平行断层右旋滑动亏损速率则基本上自东向西逐渐增加；MHT 3条剖面拟合结果也反映出其存在很强的闭锁。根据估算的此次Mw7.8地震的复发周期230年和最近500多年发生的大地震分布，认为MHT整条段落尤其是尼泊尔西部与印度接壤处和可能还没有破裂的不丹地区依然有发生8级大地震的危险。Abstract: An Mw7.8 earthquake occurred in Nepal on April 25, 2015. The epicenter is (28.147°N, 84.708°E). It occurred as a result of fault thrusting on the main thrust interface between the subducting India plate and the overriding Eurasia plate to the north. The movement and locking of the main Himalaya thrust (MHT)fault and its seismic potential need us to do some further study. We collected 6 sets of GPS horizontal velocity field covering Indian plate, Qinghai-Tibet plateau and adjacent areas, and transformed them to an approximately unified reference frame. Then we used the negative dislocation model of DEFNODE to invert the spatial fault locking and fault slip deficit in the MHT, and analyzed the far-near field movement of the fault combining with the GPS velocity profiles.The results show that the locking depth of the MHT is about 18-24 km, and the locking width of the fault plane reaches 102-136 km. The locking depth is deeper and the locking width of the plane is wider in the fault between two major historical earthquakes rupture area and the fault of no historical earthquakes. The Mw7.8 Nepal Earthquake is located in the area between two major historical earthquakes rupture area where the fault is relatively deeply and tightly locked. The total slip deficit and vertical compressional slip deficit of the MHT gradually reduce from east to west. The vertical compressional slip deficit of the easternmost segment of the fault is about 28.0 mm/a, which gradually reduces in the tightly locked area, and decreases to about 6.4 mm/a at the westernmost segment. The dextral strike-slip of parallel slip deficit changes from 0.4-7.0 mm/a from east to west in Nepal. The three profile model fitting results of the MHT are consistent with the slip deficit of corresponding fault plane, and it is very like the profile results across the middle segment of the Longmenshan fault before the 2008 Wenchuan earthquake, which reflects that the MHT is tightly locked. We estimated that the Mw7.8 earthquake recurrence is about 230 years according to the peak slip and the vertical compressional slip deficit at the epicenter. Because only the 1905 India Kangra 8.0 earthquake, the 1934 Nepal Bihar 8.1 earthquake and this Nepal Mw7.8 earthquake happened for recent 500 years, it is possible that M8 earthquake happens on this fault especially on the border of Nepal and India and possibly unruptured Bhutan.
-
图 5 横跨喜马拉雅主逆冲断层剖面的GPS站速度和最优模型结果拟合
Figure 5. GPS Velocity Profiles Across the Main Himalaya Thrust Fault and the Model Fitting Results
表 1 本文使用的GPS速度场资料
Table 1. GPS Velocity Field Used in This Work
序号 区域 点数 时间/年份 框架 融合时使用的公共站个数 融合后残差均方差/(mm·a-1) 资料来源 1 中国大陆 354 1991~2001 ITRF97 - - 文献[20] 2 青藏高原 45 1991~2001 ITRF97 40 0.49 文献[21] 3 尼泊尔喜马拉雅 97 1995~2006 ITRF2000 33 0.43 文献[6] 4 青藏高原及周边 726 1998~2004 ITRF2000-Eurasia 10 1.48 文献[22] 5 印度及邻区 229 1995~2007 ITRF2000 36 0.35 文献[23] 6 青藏高原 750 1998~2013 ITRF2008-Eurasia 11 1.72 文献[24] 注：ITRF2000-Eurasia为欧亚基准，是ITRF2000的结果去除了相对于欧亚板块的整体运动量。ITRF2008-Eurasia以此类推 
下载: 导出CSV
-
[1] 邓起东, 程绍平, 马冀, 等.青藏高原地震活动特征及当前地震活动形势[J].地球物理学报, 2014, 57(7):2025-2042 doi: 10.6038/cjg20140701 Deng Qidong, Cheng Shaoping, Ma Ji, et al. Seismic Activities and Earthquake Potential in the Tibetan Plateau[J]. Chinese J. Geophys., 2014, 57(7):2025-2042 doi: 10.6038/cjg20140701 [2] 许志琴, 王琴, 曾令森, 等.高喜马拉雅的三维挤出模式[J].中国地质, 2013, 40(3):671-680 http://kns.cnki.net/KCMS/detail/detail.aspx?filename=dizi201303003&dbname=CJFD&dbcode=CJFQ Xu Zhiqin, Wang Qin, Zeng Lingsen, et al. Three-Dimensional Extrusion Model of the Great Himalaya Slice[J]. Geology in China, 2013, 40(3):671-680 http://kns.cnki.net/KCMS/detail/detail.aspx?filename=dizi201303003&dbname=CJFD&dbcode=CJFQ [3] Ader T, Avouac J P, Jing L Z, et al. Convergence Rate Across the Nepal Himalaya and Interseismic Coupling on the Main Himalayan Thrust:Implications for Seismic Hazard[J]. J. Geophys. Res., 2012, 117(B4):389-399 doi: 10.1029/2011JB009071/abstract [4] Ponraj M, Miura S, Reddy C D, et al. Estimation of Strain Distribution Using GPS Measurements in the Kumaun Region of Lesser Himalaya[J]. Journal of Asian Earth Science, 2010, 39(6):658-667 doi: 10.1016/j.jseaes.2010.04.037 [5] 王琪, 游新兆, 王文颖, 等.跨喜马拉雅的GPS观测与地壳形变[J].地壳形变与地震, 1998, 18(3):43-50 http://www.cqvip.com/QK/95685X/1998003/3158852.html Wang Qi, You Xinzhao, Wang Wenying, et al. GPS Measurement and Current Crustal Movement Across the Himalaya[J]. Crustal Deformation and Earthquake, 1998, 18(3):43-50 http://www.cqvip.com/QK/95685X/1998003/3158852.html [6] Bettinelli P, Avouac J P, Flouzat M, et al. PlateMotion of India and Interseismic Strain in the Nepal Himalaya from GPS and DORIS Measurements[J]. Journal of Geodesy, 2006, 80(8-11):567-589 doi: 10.1007/s00190-006-0030-3 [7] Stevens V L, Avouac J P. Interseismic Coupling on the Main Himalayan Thrust[J]. Geophysical Research Letters, 2015, 42(14):5828-5837 doi: 10.1002/2015GL064845 [8] McCaffrey R. Crustal Block Rotations and Plate Coupling[J]. Plate Boundary Zones, AGU Geodynamics Series, 2002, 30:101-122 doi: 10.1029/GD030p0101/summary [9] 赵静, 武艳强, 江在森, 等.芦山地震前龙门山断裂带闭锁程度与变形动态特征研究[J].地震学报, 2013, 35(5):681-691 http://www.wenkuxiazai.com/doc/89d957974028915f814dc23c.html Zhao Jing, Wu Yanqiang, Jiang Zaisen, et al. Fault Locking and Dynamic Deformation of the Longmenshan Fault Zone Before the 2013 Lushan Ms7.0 Earthquake[J]. Acta Seismologica Sinica, 2013, 35(5):681-691 http://www.wenkuxiazai.com/doc/89d957974028915f814dc23c.html [10] 赵静, 江在森, 武艳强, 等. Defnode负位错模型反演结果的可靠性和稳定性分析[J].大地测量与地球动力学, 2013, 33(1):21-24 http://www.wenkuxiazai.com/doc/202faa4d227916888586d74e.html Zhao Jing, Jiang Zaisen, Wu Yanqiang, et al. Analysis of Reliability and Stability of Inversion Result with Negative Dislocation Model of Defnode[J]. Journal of Geodesy and Geodynamics, 2013, 33(1):21-24 http://www.wenkuxiazai.com/doc/202faa4d227916888586d74e.html [11] Manaker D M, Calais E, Freed A M, et al.Interseismic Plate Coupling and Strain Partitioning in the Northeastern Caribbean[J]. Geophys. J. Int., 2008, 174(3):889-903 doi: 10.1111/gji.2008.174.issue-3 [12] McCaffrey R. Block Kinematics of the Pacific-North America Plate Boundary in the Southwestern United States from Inversion of GPS, Seismological, and Geologic Data[J]. J. Geophys. Res., 2005, 110(B7):139-148 doi: 10.1029/2004JB003307/citedby [13] McCaffrey R, Qamar A I, King R W, et al. FaultLocking, Block Rotation and Crustal Deformation in the Pacific Northwest[J]. Geophys. J. Int., 2007, 169(3):1315-1340 doi: 10.1111/gji.2007.169.issue-3 [14] Wallace L M, Beavan J, McCaffrey R, et al. Subduction Zone Coupling and Tectonic Block Rotations in the North Island, New Zealand[J]. J. Geophys. Res., 2004, 109(B12):159-163 doi: 10.1029/2004JB003241/full [15] Wallace L M, Stevens C, Silver E, et al.GPS and Seismological Constraints on Active Tectonics and Arc-Continent Collision in Papua New Guinea:Implications for Mechanics of Microplate Rotations in a Plate Boundary Zone[J]. J. Geophys. Res., 2004, 109(B5):219-240 http://adsabs.harvard.edu/abs/2004JGRB..109.5404W [16] Wallace L M, Beavan J, McCaffrey R, et al.Balancing the Plate Motion Budget in the South Island, New Zealand Using GPS, Geological and Seismological Data[J]. Geophys. J. Int., 2007, 168(1):332-352 doi: 10.1111/gji.2007.168.issue-1 [17] 赵静, 江在森, 牛安福, 等.川滇菱形块体东边界断层闭锁程度与滑动亏损动态特征研究[J].地球物理学报, 2015, 58(3):872-885 doi: 10.6038/cjg20150316 Zhao Jing, Jiang Zaisen, Niu Anfu, et al. Study on Dynamic Characteristics of Fault Locking and Fault Slip Deficit in the Eastern Boundary of the Sichuan-Yunnan Rhombic Block[J]. Chinese J. Geophys., 2015, 58(3):872-885 doi: 10.6038/cjg20150316 [18] 赵静, 江在森, 武艳强, 等.汶川地震前龙门山断裂带闭锁程度和滑动亏损分布研究[J].地球物理学报, 2012, 55(9):2963-2972 doi: 10.6038/j.issn.0001-5733.2012.09.015 Zhao Jing, Jiang Zaisen, Wu Yanqiang, et al. Study on Fault Locking and Fault Slip Deficit of the Longmenshan Fault Zone Before the Wenchuan Earthquake[J]. Chinese J. Geophys., 2012, 55(9):2963-2972 doi: 10.6038/j.issn.0001-5733.2012.09.015 [19] Savage J C, Gan W J, Svarc J L. Strain Accumulation and Rotation in the Eastern California Shear Zone[J]. J. Geophys. Res., 2001, 106(B10):21995-22007 doi: 10.1029/2000JB000127 [20] Wang Q, Zhang P Z, Freymueller J T, et al. Present-Day Crustal Deformation in China Constrained by Global Positioning System Measurements[J]. Science, 2001, 294(5542):574-577 doi: 10.1126/science.1063647 [21] Chen Q Z, Freymueller J T, Wang Q, et al. A Deforming Block Model for the Present-Day Tectonic of Tibet[J]. J. Geophys. Res., 2004, 109, B01403, doi: 10.1029/2002JB002151 [22] Gan W J, Zhang P Z, Shen Z K, et al. Present-DayCrustal Motion with the Tibetan Plateau Inferred from GPS Measurements[J]. J. Geophys. Res., 2007, 112(B08):582-596 http://adsabs.harvard.edu/abs/2007JGRB..112.8416G [23] Banerjee P, Bürgmann R, Nagarajan B, et al. Intraplate Deformation of the Indian Subcontinent[J]. Geophys. Res. Lett., 2008, 35, L18301, doi: 10.1029/2008GL035468 [24] Liang S, Gan W, Shen C, et al. Three-Dimensional Velocity Field of Present-Day Crustal Motion of the Tibetan Plateau Derived from GPS Measurements[J]. J. Geophys. Res.:Solid Earth, 2013, 118, doi: 10.1002/2013JB010503 [25] 占伟, 武艳强, 梁洪宝, 等. GPS观测结果反映的尼泊尔Mw7.8地震孕震特征[J].地球物理学报, 2015, 58(5):1818-1826 doi: 10.6038/cjg20150532 Zhan Wei, Wu Yanqiang, Liang Hongbao, et al. Characteristics of the Seismogenic Model for the 2015 Nepal Mw7.8 Earthquake Derived from GPS Data[J]. Chinese J. Geophys., 2015, 58(5):1818-1826 doi: 10.6038/cjg20150532 [26] USGS. The April-May 2015 Nepal Earthquake Sequence[EB/OL]. http://earthquake.usgs.gov/learn/topics/Nepal_Slides.pdf, 2015 [27] Avouac J P, Meng L, Wei S, et al. Lower Edge of Locked Main Himalayan Thrust Unzipped by the 2015 Gorkha Earthquake[J]. Nature Geoscience. 2015, 8(9):1-6 https://www.nature.com/articles/ngeo2518 [28] 张勇, 许力生, 陈运泰. 2015年尼泊尔Mw7.9地震破裂过程:快速反演与初步联合反演[J].地球物理学报, 2015, 58(5):1804-1811 doi: 10.6038/cjg20150530 Zhang Yong, Xu Lisheng, Chen Yuntai. Rupture Process of the 2015 Nepal Mw7.9 Earthquake:Fast Inversion and Preliminary Joint Inversion[J]. Chinese J. Geophys., 2015, 58(5):1804-1811 doi: 10.6038/cjg20150530 [29] Lavé J, Avouac J P. Active Folding of Fluvial Terraces Across the Siwaliks Hills, Himalayas of Central Nepal[J]. J. Geophys. Res., 2000, 105(B3):5735-5770 doi: 10.1029/1999JB900292 [30] Bilham R, Larson K, Freymueller J, et al. GPS Measurements of Present-Day Convergence Across the Nepal Himalaya[J]. Nature, 1997, 386(6620):61-64 doi: 10.1038/386061a0 [31] Feldl N, Bilham R. Great Himalayan Earthquakes and the Tibetan Plateau[J]. Nature, 2006, 444(7116):165-170 doi: 10.1038/nature05199 [32] 江在森, 刘经南.应用最小二乘配置建立地壳运动速度场与应变场的方法[J].地球物理学报, 2010, 53(5):1109-1117 http://kns.cnki.net/KCMS/detail/detail.aspx?filename=dqwx201005012&dbname=CJFD&dbcode=CJFQ Jiang Zaisen, Liu Jingnan. The Method in Establishing Strain Field and Velocity Field of Crustal Movement Using Least Squares Collocation[J]. Chinese J. Geophys., 2010, 53(5):1109-1117 http://kns.cnki.net/KCMS/detail/detail.aspx?filename=dqwx201005012&dbname=CJFD&dbcode=CJFQ -
点击查看大图
图(5) / 表(2)
计量
- 文章访问数: 1645
- HTML全文浏览量: 62
- PDF下载量: 621
- 被引次数: 0
