地球物理学报  2016, Vol. 59 Issue (11): 4022-4034   PDF    
GPS揭示的郯庐断裂带中南段闭锁及滑动亏损
李彦川1,2 , 单新建1 , 宋小刚1 , 姜宇1,2 , 甘卫军1 , 屈春燕1 , 王振杰2     
1. 地震动力学国家重点实验室, 中国地震局地质研究所, 北京 100029;
2. 地球科学与技术学院, 中国石油大学(华东), 青岛 266580
摘要: 利用华北地区2009-2014年GPS水平运动速度场数据,采用块体负位错模型反演了郯庐断裂带中南段断层深部滑动速率、断层闭锁程度分布、断层滑动亏损速率分布及地震矩积累率,结合地表应变率分布,对郯庐断裂带中南段深、浅部形变、应变特征以及华北地区的地壳形变模式进行了分析.结果表明:郯庐断裂中南段的北端主要为右旋走滑特性,南端则表现为右旋走滑兼拉张性运动,断层滑动速率在0.9 mm·a-1至1.2 mm·a-1,且沿断层走向由北至南逐次增大.断层闭锁程度分布沿走向分布不均一,断层闭锁深度由最北端的27 km增加到中段的32 km,至最南端变为5 km,断层闭锁最深处与1668年郯城MS8.5震中位置相对应.断层滑动亏损速率沿走向由0.9 mm·a-1增加到1.2 mm·a-1,沿倾向由地表至深部逐渐减小为0 mm·a-1.地震矩积累率在郯庐断裂带中南段郯城附近较大,而地表对应区域为第二应不变分量的低值区.华北地区地壳变形以块体运动为主,块体内部应变及断层闭锁产生的负位错效应次之;郯庐断裂带中南段断层形变沿走向呈条带状分布,形变宽度单侧小于50 km,形变量不超过1 mm·a-1,且上盘形变略大于下盘.
关键词: 郯庐断裂      GPS      断层闭锁      滑动亏损     
Fault locking and slip rate deficit on the middle and southern segment of the Tancheng-Lujiang fault inverted from GPS data
LI Yan-Chuan1,2, SHAN Xin-Jian1, SONG Xiao-Gang1, JIANG Yu1,2, GAN Wei-Jun1, QU Chun-Yan1, WANG Zhen-Jie2     
1. State Key Laboratory of Earthquake Dynamics, Institute of Geology, CEA, Beijing 100029, China;
2. School of Geosciences, China University of Petroleum(East China), Qingdao 266580, China
Abstract: By using GPS-derived horizontal velocities of 2009-2014 and a dislocation model, we inverted the slip rate, fault coupling and slip rate deficit on the middle and southern segments of the Tanlu fault. We also calculated the moment accumulation rate based on the fault slip rate and the fault locking depth and analyzed crustal deformation in East China. Along the middle and southern segments of the Tanlu fault, the results predict 0.9 mm·a-1 of right-lateral strike-slip at the northern end, gradually changing to 1.2 dextrorotary extension at the southern end. The overall pattern of the interseismic locking depth is unevenly distributed along the middle and southern segments of the Tanlu fault. We estimated fault coupling down to 27 km depth or more for the northern part, and 25~32 km for the central part which is coincident with the epicenter of the 1668 Tancheng earthquake. The southern part, however, is only locked from 5 km to 25 km depth. Similar to spatial distribution of the fault coupling, the slip rate deficit changes from a lower value in the northern part to a higher level in the southern part, ranging from 0.9 mm·a-1 to 1.2 mm·a-1 along the middle and southern segments of the Tanlu fault. And it gradually decreases to 0 mm·a-1 along the dip. The seismic moment accumulation rate per unit length of the fault is high on the fault plane beneath Tancheng city; however, it is negatively correlated with the strain rate on the surface. Our results also show that tectonic deformation occurring in East China is dominated by block rotation. While most sites in East China are affected by elastic strain within blocks to some degree. However, only sites in a narrow strip (less than 50 km with velocities less than 1 mm·a-1) along the middle and southern segments of the Tanlu fault are affected by elastic strain loading due to fault coupling. Compared with sites on the footwall, sites on the hanging wall are affected more by strain loading due to fault coupling..
Key words: Tanlu fault      GPS      Fault locking      Slip rate deficit     
1 引言

郯庐断裂带横跨中国东部,是亚洲东部著名的深大断裂带之一,其新活动构造诱发了一系列近代地震,对现代地震活动起着重要的控制作用(方仲景等,1980王小凤等,2000刘保金等,2015).郯庐断裂带是由多组呈斜列分布的多条断裂组成,按构造习性、演化历史和地震活动性可分为北、中、南三段,其中渤海湾以南称为中南段(国家地震局地质研究所,1987Xu,1993万桂梅等,2009),该段作为郯庐断裂的主体,总体走向北18—25°东,是一条中生代的大陆裂谷构造,是一条从地表直达上地幔的贯通断裂(高维明等,1980),第四纪以来活动断裂最为发育,是我国主要的地震活动带之一(吴大铭等,1981),先后发生了多次M6以上地震(图 1).2008年5月12日汶川8.0地震和2011年3月11日日本宫城地震分别对郯庐断裂带中南段造成了不同程度的影响(刘东旺等,2009王敏等,2011杨少敏等,2011尹海涛等,2013;尹京苑等;2015),加强对该断裂现今地壳形变的研究和认识,对评价区域地震危险性具有重要的价值.更好地量化郯庐断裂的地壳形变,一定意义上要依赖于断层的滑动速率;同时,断层滑动速率为我们辨别地壳形变模式、了解应力加载过程、探索大地震复发周期及评估该区域地震危险性提供了不可或缺的约束(Cavalié et al., 2008).此外,了解断层的闭锁、滑动亏损及应力积累的空间分布,也为我们更好地理解断层的地震活动性奠定了基础(Jolivet et al., 2013).

图 1 研究区构造及历史地震分布黑色实线表示郯庐断裂,黑色空心圈为M≤5历史地震分布,白色实心圈为M>5历史地震分布,黄色线框为图 6微震剖面位置,左上角插图表示研究区域所在位置. Fig. 1 Tectonics and historical earthquakes in the study area Black line represents the Tanlu fault. Black circles are earthquakes with magnitude equal or smaller than 5 and white circles correspond to earthquakes with magnitude greater than 5. Yellow rectangle marks location of profiles of micro-earthquakes in Figure 6. Inset map on the upper left shows the study area.

利用现代大地测量测量手段(GPS)对郯庐断裂震间地壳形变及应变的研究已有很多(郭良迁和应绍奋,1998黄立人和郭良迁,1998许才军等,2002Wang et al., 2009郭良迁等,2011Wang et al., 2011刘晓霞等,2012张希等,2013朱泽等,2014),这些研究结果不仅提供了郯庐断裂中南段现今滑动速率,也勾勒出了该断裂的活动特征,但仍有一些问题有待解决,如不同研究结果得到的郯庐断裂活动习性不一致;未有对断层闭锁程度和滑动亏损空间连续变化情况的研究;依据块体-应变模型对郯庐断裂形变的研究,未考虑断层闭锁产生的负位错效应;地表应变、断层闭锁程度、断层滑动亏损、地震矩积累率等断层深、浅部过程的对应关系有待于进一步揭示.本文利用2009—2014年密集的GPS水平速度场结果,采用块体负位错模型(TDEFNODE程序)(McCaffrey,2009)反演了郯庐断裂现今滑动速率、断层闭锁、断层滑动亏损速率及地震矩积累率,结合该区域的应变结果,分析了郯庐断裂中南段深浅部地壳形变特征,并将实测GPS速度分解,量化了该区域的地壳形变,最后讨论了反演结果的不确定性.

2 数据及模型 2.1 GPS数据

作为现代地壳运动获取的主要手段,GPS被广泛应用于形变监测及地球物理反演研究(Wang et al., 2001).考虑到2008年5月12日汶川地震同震、震后效应可能对周围流动测站造成持续影响(Zhao et al., 2015),收集了研究区中国地壳运动观测网络、中国大陆构造环境监测网络2009—2014年流动站数据(分别在2009、2011和2013年进行了复测)和1998—2014年连续站的GPS数据,采用GAMIT10.4软件(Herring et al., 2010)联合中国大陆周围10个IGS站进行宽松约束处理,其中高度角设为10°,加入EGM08重力模型并考虑了固体潮、大气延迟效应、电离层改正等;将解算得到的单天解与SOPAC全球单天解进行合并,剔除时间序列中跳值点;为了得到“干净”的时间序列,分别在时间序列中对2001年可可西里地震、2004年苏门答腊地震、2008年汶川地震、2010年玉树地震、2011年日本311地震、2013年芦山地震和2014年于田地震的同震及震后效应进行了剔除;选取全球均匀分布的50个稳定IGS站,采用七参数转换,得到ITRF2008下速度;最终利用欧亚块体的欧拉矢量(Altamimi et al., 2011),将ITRF2008下速度场转换为稳定欧亚板块下的速度场(Zhao et al., 2015)(图 2).在模型反演前,对方向、大小明显偏离区域运动背景的GPS测站进行了剔除,最终选定了653个GPS测站数据,其中流动GPS测站水平方向平均误差为0.40 mm·a-1,连续运行GPS测站水平向平均误差为0.26 mm·a-1.

图 2 稳定欧亚板块下的GPS速度场 黑色箭头表示处理得到的GPS速度,蓝色箭头为反演得到的郯庐断裂滑动速率,误差椭圆为95%置信区间,红色线框为图 3 GPS速度剖面,灰色线框表示断裂分段;右上角插图中红色框线为研究区所在位置,黑色实线为划分的块体边界. Fig. 2 GPS velocities with respect to the stable Eurasian plate Black arrows are processed GPS velocities. Blue arrows are slip rates of the fault from inversion. Error ellipses show with 95% confidence level. Red rectangle marks the profile of velocities in Figure 3 and the grey is flags for fault segmentation. The red rectangle in the inset map on the upper left shows the study area and black bold lines mark boundaries of blocks.
2.2 反演模型

文中用到的块体负位错模型(McCaffrey, 2002,2005)假设地壳形变是由块体旋转、块体内部应变、块体边界断层闭锁产生的负位错效应共同作用下的地表弹性变形之和(公式1),该模型结合了块体旋转、块体均匀应变和Okada位错,可在大地测量(GPS、InSAR、水准)、地质(如断层滑动速率)及地震(如断层同震滑动方位角)等数据的约束下,采用网格搜索和模拟退火法同时反演块体旋转欧拉极、断层滑动速率、块体边界断层闭锁系数(PHI)等.断层闭锁系数取值在0~1,PHI=0表示断层自由滑动,断层两盘的相对运动不产生应变能积累;PHI=1表示断层完全闭锁,断层两盘相对运动量将完全转化为应变能;PHI处于0和1之间表示断层处于部分滑动状态,两盘相对运动量部分转化为应变能(McCaffrey et al., 2007).模型由开源程序TDEFNODE实现(http://web.pdx.edu/~mccaf/www/defnode/; McCaffrey,2009).模型原理、反演过程及参数控制,相关文献已有详细的描述(McCaffrey,2002赵静等,2012).公式(1)为

(1)

其中Vsf表示实测地表速度,Vbr为块体旋转引起的速度,Vis为块体内部应变引起的速度,Vfs为断层闭锁负位错效应引起的速度.

前人根据地质、地震及GPS速率资料将我国东部划分为华南和华北两个一级活动地块,其中华北地块划分为两个亚板块:华北平原和鲁东—黄海块体(以下简称鲁东块体)(张培震等,2003Wang et al., 2011),两者以郯庐断裂中南段为界(图 2).参考前人研究成果(Zhang et al., 2005),在反演过程中,设置华北平原块体和鲁东块体内部存在均匀应变,华南块体为刚性;同时,简化郯庐带为单一深大断裂,断层走向为SSW,倾角为65°,倾向为NW(邓启东,2007).模型中断层是由一系列节点组成的面,我们的断层模型沿断层走向共有10排节点,节点之间的平均距离约为70 km,节点沿断层面分别设在0.1 km、5 km、10 km、15 km、20 km、25 km、35 km和45 km深度.尚未有研究发现郯庐断裂带中南段存在地表浅层蠕滑现象,因此在反演时对断层闭锁施加强约束(McCaffrey,2002; Wang et al., 2003),即在地表完全闭锁,在45 km深度以下自由滑动,处于两者之间的节点在反演时闭锁系数沿断层倾向随深度递减(图 3).在讨论部分,我们将对模型反演结果的不确定性做进一步分析.

图 3 郯庐断裂中南段简化示意图 黑色圆点为节点所在位置,红色线表示PHI在地表完全闭锁,蓝色表示PHI自由滑动,绿色表示PHI反演时为自由参数. Fig. 3 Sketch of simplified geometry and boundary conditions for the model Block dots are fault nodes. The red line on the top indicates the Tanlu fault fully locked at surface and the blue line on the bottom indicates the Tanlu fault is freely slipping below 45 km. The PHI in the green area between the surface and 4 5 km at depth would be inverted from GPS data.
3 郯庐断裂中南段断层闭锁及滑动亏损特征 3.1 断层滑动速率

断层滑动速率是根据最优化模型反演的块体旋转欧拉矢量计算而来(McCaffrey,2002),反映的是大地测量时间尺度的断层在深部的运动情况.我们得到了郯庐断裂中南段断层现今滑动(图 2,表 1),为方便分析,将郯庐断裂带中南段分为A(35.5°N—37.5°N)、B(33°N—35.5°N)、C(31°N—33°N)三段(图 1).

表 1 郯庐断裂中南段断层滑动速率(右旋拉张为正)、闭锁深度、滑动亏损及地震矩积累率 Table 1 Slip rates(right-lateral with tension is positive),locking depth,slip rate deficit and seismic moment accumulation rates on the middle and southern segments of the Tanlu fault

表 1图 2中反演结果可以看出,郯庐断裂带中南段呈现右旋走滑兼拉张运动的特性,断裂运动速率为0.9~1.2 mm·a-1,且运动速率从南向北逐渐递减.断裂在南部以右旋走滑为主,张性形变为辅,拉张分量向北逐步递减,到最北端基本为0,断层变为纯右旋走滑.为验证模型的可靠性,横跨断裂做了三个速度剖面a、b和c,剖面估值(表 1图 4)与模型反演值具有较好的一致性.平行断裂方向c剖面值与模型结果相差较大,可能与该区域处于块体交界的地方,形变相对较为复杂有关(Thatcher,2007).刘晓霞等(2012)利用2007—2009年两期GPS速度场资料,基于块体的整体旋转与均匀应变模型反演了郯庐断裂中南段的运动,结果显示断裂在北段呈右旋拉张、南段呈左旋挤压,但走滑分量占主导;此外郭良迁等(2011)利用1999—2009年华北地区GPS速度场为约束,同样采用块体-应变模型反演了郯庐断裂的滑动速率,得到郯庐断裂带运动速度从北向南逐渐减小(1.24 ~1.06 mm·a-1),呈现右旋挤压特征.上述结论与我们反演的结果相互矛盾,原因前者用到的数据均为2009年以前的数据,而我们反演时GPS数据跨度为2009—2014年,期间发生了2011年3月11日日本宫城MW9.0巨震,位于我国的GPS连续站显示同震在华东地区造成了3~10 mm的水平位移(主要为东西向)(王敏等,2011殷海涛等,2013),导致东北和华北地区一系列北北东走向的断裂产生了不同程度的张性应变,对郯庐断裂起到了拉张作用.这些远场同震位移虽然不大,却对郯庐断裂(中南段)运动性质产生了影响.

3.2 断层闭锁程度分布

图 5给出了郯庐断裂带中南段闭锁程度沿断层走向的三维分布,表 1给出了统计的A、B和C段的闭锁深度.可以看出,断层闭锁深度在沿走向从北向南呈现先深后浅的特征,断层闭锁最深可达32 km,且正好位于1668年7月25日郯城MS8.5震中附近.地震实际上是在区域构造应力作用下,应变在活动断裂带上不断积累并达到极限状态后而突发失稳破裂的结果(Scholz, 1990,1998张培震等,2013);而断层只有在闭锁(或部分闭锁)的状态下,断层两盘相对运动才会产生亏损,转化为能量并逐渐累积,直至某次地震(或无震蠕滑)进行释放,并进入下一个地震周期(Wallace et al., 2004).由此可以推论地震发生在断层闭锁的区域.周翠英等(2013)采用郯城地震震源区重新定位的现今中小震数据,反演得到1668年郯城地震震源断层下界面深度约为32 km,上界埋深约为4 km.这与我们反演断层闭锁深度较为吻合.值得思考的是,现今距1668年郯城地震离逝仅347年,如果地震发生时断层发生了完全破裂,那么断层怎样在短时间内重新形成闭锁的.依旧精定位结果,周翠英等(2013)给出了郯城大地震断层面现代小震的深度分布,结果显示震源以南(对应图 5中B、C段交接处)的深度较浅(25 km左右),向北渐次加深,最深达33 km.上述结果与我们模型反演的结果基本一致.

图 4 GPS速度剖面图 (a)(b)(c)剖面位置在图 2中,左图为平行断裂速率,右图为垂直断裂速率;图中黄色方框用于估计断层滑动速率,红色线为中值. Fig. 4 GPS slip rates along and across the middle and southern segments of the Tanlu fault Locations of each profile are shown in Figure 2. Figures on the left show fault-parallel velocities and fault-normal velocities on the right. The stations in light yellow rectangle would be used for estimating slip rates.
图 5 郯庐断裂中南段断层闭锁程度分布 Fig. 5 Locking degrees of the middle and southern segments of the Tanlu fault

收集了1970年以来沿郯庐断裂带中南段分布的现代微震(图 1),并将震源深度进行垂直展布(图 6),结果显示几乎所有微震均发生于30 km以内,这表明在断层以下30 km内应变能在不断地释放,也反映了断层处在闭锁状态.文中反演的断层闭锁深度深于当代地震的震源深度,这证明了我们的反演结果具有合理性.张希等(2013)利用2009—2011年全国GPS水平运动场资料,借助负位错反演得到郯庐断裂带平均闭锁深度仅为16 km,该结果与我们反演结果不一致,这可能是模型及数据集不一致导致的;朱泽等(2014)虽然用同样的方法对郯庐带中南段进行了反演,但其模型设置仅为两排节点,且断层倾角(SE50°)与现有的认识(邓启东,2007)存在不一致,但结果与我们的反演结果基本一致.

图 6 历史微震剖面 Fig. 6 Profile of the historical micro-earthquakes
3.3 断层滑动亏损分布

量化断层面上应变能的分布及积累速率,对该断层地震危险性的判断具有重要意义(Jolivet et al., 2013).断层滑动亏损速率由断层滑动速率乘以断层面闭锁系数得到(图 7).由于断层闭锁产生的滑动亏损,将以应变能的方式在断层面及附近累积(McCaffrey et al., 2007).从图 7可以看出,滑动亏损速率整体分布于断层闭锁程度基本一致,且沿断层走向由南向北量值逐次减小(1.2~0.9 mm·a-1);断层倾向上,滑动亏损速率从地面到深部逐渐变小,这与实际地球物理现象是吻合的(Wang et al., 2003).

图 7 郯庐断裂带中南段断层滑动亏损速率分布 Fig. 7 Slip rate deficits of the middle and southern segments of the Tanlu fault

滑动亏损速率仅表示断层两盘滑动量转化为应变能快慢的量,但应变能积累的快并不表示应变能总量就一定大,2008年5月12日汶川地震前龙门山断裂带两盘相对低速的滑动速率就是一个很好的例证(Zhang,2013).为对郯庐断裂带中南段地震危险性有较为客观的分析,我们对断裂滑动亏损和闭锁深度进行积分(公式2),计算了各分段(A、B和C)的地震矩积累率(表 1).公式为

(2)

其中M0为沿断层走向单位长度的地震矩积累率,j为某段断层地震矩积累率,l为断层长度,PHI为断层面闭锁系数,η为岩石剪切模量(取30 GPa),dj为断层闭锁深度,vj为断层面滑动速率.

表 1中,地震矩积累速率与断层闭锁深度及滑动亏损速率成相关,反映的是地震矩积累的快慢程度.郯庐断裂中南段中B段地震矩积累速率较大,这可能反映了该段能量积累较快.刘晓霞等(2012)根据块体模型反演GPS数据的结果推断郯城附近具有发生强震的孕育背景,秦四清等(2014)运用孕震断层多锁固段脆性破裂理论及相关预测方法,在此段得到了相似的结论.

4 华北地区应变特征

文中用到的块体负位错模型利用地表形变速率,反演的是断层深部构造性变特征.为将断层深部特征与地表形变进行对比分析,我们首先采用L曲线方法,将华北地区的GPS速度场进行0.4°×0.4°插值,并计算了既包含正应变又包含剪应变的应变第二不变分量(Wang and Wright, 2012)(图 8).图 8中,第二应变不变分量高值主要分布在环首都圈区域,该地区存在一系列的次级断裂;同时,第二应变不变分量的高值区也对应着主压应变的高值区.沿郯庐带走向由北向南,A段两侧主应变发生了方向改变,且该区域第二应变不变分量相对较高;郯庐带B段主应变几乎为零,第二应变不变分量也处于低值区;C段主应变率是连续变化的,且以主张应变为主,这与王敏等(2011)结论相同,推测该段可能仍在受日本311宫城MW9.0大地震的影响.

将地表应变分布与断层在深部的闭锁深度、滑动亏损速率及断层地震矩积累速率进行比较(图 8,表 1),发现前者与后三者并不是正相关,尤其是在郯庐带B段,即断层深部高地震矩积累率在地表的映射关系却是低应变积累率,一种可能的解释是该段在浅地壳已经积累了一定的应变能,应变累积速度反而减小了,由此推测,郯庐断裂带中南段B段(即郯城附近)可提供地震孕育的环境.需要强调的是,对某条断层地震危险性的分析,不仅需要科学地分析其现今形变特征,完整的断层的结构及古地震研究也是必不可少的(Zhang,2013).

图 8 内插GPS速度场及应变场分布 (a)图中黑色箭头为实测GPS速度,红色箭头为内插得到的GPS速度,误差椭圆为95%置信区间; (b)图中背景颜色表示第二应变不变分量,交叉箭头为主应变. Fig. 8 Velocity(a)and strain(b)fields in East China Arrows in black are GPS velocities derived from data processing,and arrows in red are derived from interpolation using the L-curve method. Background color in right panel denotes the magnitude of the second invariant of the strain rate tensor,and black arrows show the compressional and extensional principal strain rates.
5 讨论 5.1 GPS速度分解

以GPS为约束,前人运用的块体模型(黄立人和郭良迁,1998许才军等,2002Zhang et al., 2005; Wang et al., 2009; 郭良迁等,2011刘晓霞等,2012)和考虑断层闭锁产生弹性负位错效应的块体模型(Wang et al., 2011张希等,2013朱泽等,2014)对华北地区的地壳形变进行了研究,均对数据进行了很好的拟合,但未对GPS数据进行各分量的分解,即块体旋转、块体应变和断层闭锁的GPS分量.我们在反演过程中,将华北地区实测GPS速度场转化到稳定华南块体框架下,并将该框架下理论的GPS速度场包含的各分量进行了分解(图 9).

图 9中,相对于稳定的华南块体,华北地区的GPS速度(黑色箭头)呈逆时针旋转,这与江在森等(2009)给出的华南基准速度场类似.郯庐断裂中南段南端(C段)两侧,GPS速度呈现右旋拉张运动,向北则两盘差异运动不明显.由速度场的分布可以得出,华北地区主要以块体运动为主导(绿色箭头),且华北块体和华东块体之间的相对运动,在郯庐断裂中南段南端表现为右旋走滑兼拉张特性,往北则相对运动不突出,这与我们模型计算的郯庐断裂中南段运动性质是一致的.块体内部应变部分的GPS速度(红色箭头)主要分布于块体边界,且量值小于1.0 mm·a-1,呈背离块体中心运动,这与反演模型假设块体为均匀应变有关,即远离块体中心,应变分量变大(李延兴等,2001).断层闭锁影响的GPS速度分量,沿断层走向两侧展布,量值不足1.0 mm·a-1,且在断层两盘呈反方向运动,这是由于断层闭锁产生弹性负位错阻碍两盘的相对运动,与断层右旋走滑的运动性质是相吻合的;图 9中断层闭锁影响的GPS分量,主要分布于郯庐断裂带中南段的上盘(西盘),这符合华北地区整体运动的动力学背景,即上盘是主动盘(张培震等,2013);远离断裂两侧,GPS速度呈现快速下降,可推断郯庐断裂带中南段变形宽度很窄(<50 km),从图 4中GPS速度剖面可推测出同样结论.

图 9 相对于稳定华南块体的GPS速度场分解 图中黑色箭头表示理论上相对于华南块体的速度,绿色箭头为块体旋转部分,红色箭头为块体内部应变部分,蓝色箭头为断层闭锁引起的负位错效应部分. Fig. 9 Decomposition of velocity field with respect to stable South China Block arrows are velocities relative to the stable South China block. Green arrows are velocities induced by block rotation. Red arrows are velocities which lead to block internal strain. Blue arrows are velocities due to fault locking.
5.2 反演结果的不确定性分析

统计了GPS速度残差及其分布(图 10),绝大部分速度残差(<2 mm·a-1)在其误差椭圆内,残差方向具有随机性,残差分布符合高斯正态分布;华北和鲁东块体内部应变残差均小于1.0nstrain/a,证明了模型反演结果内符合精度较高.值得注意的是,由于模型假设块体内部形变为均匀形变,未考虑形变、地壳介质、波速结构等空间分布的不均一性,这可能是反演结果不确定性的来源之一(Zhang et al., 2013).此外,沿块体边界以及块体内部次级断层周围(图 10),GPS速度残差相对较大(约2 mm·a-1),这是由于块体边界形变的复杂性引起的(Thatcher,2007; Wang et al., 2011),同时也说明块体负位错模型并非可完全模拟该地区地壳形变,这也给反演结果带来了不确定性.断层近场GPS观测站点的分布(特别是断裂带附近站点的数量和位置)也对反演结果有一定影响(赵静等,2013).

图 10 块体应变及应变残差,GPS速度残差分布 误差椭圆为1.5倍中误差,95%的置信区间,蓝色箭头为块体内部主应变,绿色箭头为应变残差. Fig. 10 Residuals(observed minus calculated velocities)of GPS velocities,strain in blocks and its residuals GPS velocity residuals are within 95% confidence ellipses in 1.5 sigma. Blue cross represents the strain and yellow cross represents its residuals.

模型参数的设置(如是否设置块体内部应变、参考框架的选择等)对反演结果也会有影响.统计了不同参数设置下的反演GPS的结果精度(表 2,表 3),可以看出,不同参考框架对反演结果的影响相对较小,不同稳定块体作为参考框架时对反演结果几乎不产生影响;块体是否为刚性块体(即是否存在内部应变)对反演结果有较大的影响;反演结果对断层倾角不敏感,断层初始深度对结果有较大影响.总之,我们利用块体负位错模型对华北地区的地壳形变取得了较好的描述,但模型结果亦受不确定因素的影响.

表 2 不同模型参数设置 Table 2 Different models for inversion
表 3 断层倾角及初始深度结果对比 Table 3 Different fault dips and initial depths for inversion
6 结论

利用华北地区2009—2014年GPS水平向速度场数据,采用块体负位错模型对郯庐断裂带中南段断层深部滑动速率、断层深部闭锁系数、滑动亏损速率及地震矩积累率进行了反演,结合地表应变分布,对郯庐断裂带中南段形变、应变及华北地区地壳形变模式进行了分析,取得了如下认识:

(1) 模型估值及横跨断层GPS速度剖面均显示,受2011年3月11日日本宫城MW9.0地震影响,郯庐断裂带中南段表现为右旋走滑兼拉张的运动性质,滑动速率为0.9~1.2 mm·a-1,且断层沿走向逐次递增,在北端以右旋走滑运动为主,南端兼拉张特性.

(2) 郯庐断裂带中南段断层闭锁深度沿走向不均一分布,呈先加深后变浅(27~32~5 km),断层最深闭锁段落与1668年郯城MS8.0地震震中位置接近;断层面上滑动亏损速率沿断层走向呈递增(0.9~1.2 mm·a-1),沿倾向从地表至深部逐次减小;地震矩积累率计算结果显示,在郯庐断裂带中南段郯城附近地震矩积累较快;地表第二不变应变分量在郯城附近为低值区,推测该段可能已经积聚了较多的应变能.

(3) 华北地区地壳形变以块体旋转为主,块体内部应变和断层闭锁效应分别次之.郯庐断裂带中南段因断层闭锁产生的形变(<1 mm·a-1)沿断层两侧分布,单侧形变宽度不足50 km,且上盘形变量略大于下盘.

参考文献
Altamimi Z, Collilieux X, Métivier L. 2011. ITRF2008:an improved solution of the international terrestrial reference frame. J. Geod. , 85 (8) : 457-473. DOI:10.1007/s00190-011-0444-4
Cavalié O, Lasserre C, Doin M P, et al. 2008. Measurement of interseismic strain across the Haiyuan fault (Gansu, China), by InSAR. Earth Planet. Sci. Lett. , 275 (3-4) : 246-257. DOI:10.1016/j.epsl.2008.07.057
Deng Q D. China Active Tectonic Map(1:4000000) (in Chinese). (in Chinese) Beijing: Seismological Press, 2007 .
Fang Z J, Ding M L, Ji F J, et al. 1980. Geological analysis of the seismicity in the Tancheng-Lujiang fault zone, East China. Seismology and Geology (in Chinese) , 2 (4) : 39-45.
Gao W M, Li J L, Sun Z Y. 1980. Formation and evolution of the Yihe-Suhe continental rift. Seismology and Geology (in Chinese) , 2 (3) : 11-18.
Guo L Q, Bo W J, Yang G H, et al. 2011. Characteristics of horizontal deformation-strain field in north China from 1999 to 2009. Journal of Geodesy and Geodynamics (in Chinese) , 31 (3) : 15-19.
Guo L Q, Ying S F. 1998. Primary study on horizontal deformation and tectonic activity using GPS results in northern part of North China. Earthquake Research in China (in Chinese) , 14 (4) : 11-19.
Herring T A, King R W, McClusky S C. 2010. Documentation of the GAMIT and GLOBK software release 10.4. Cambridge:Mass. Inst. of Technol.
Huang L R, Guo L Q. 1998. GPS measurement (1995-1996) in the northern area of north China and its explanation. Seismology and Geology (in Chinese) , 20 (4) : 423-430.
Institute of Geology, C EA. Tan-Lu Fault Zone (in Chinese). (in Chinese) Beijing: Seismological Press, 1987 .
Jiang Z S, Fang Y, Wu Y Q, et al. 2009. The dynamic process of regional crustal movement and deformation before Wenchuan MS8.0 earthquake. Chinese J. Geophys. (in Chinese) , 52 (2) : 505-518.
Jolivet R, Lasserre C, Doin M P, et al. 2013. Spatio-temporal evolution of aseismic slip along the Haiyuan fault, China:implications for fault frictional properties. Earth Planet. Sci. Lett. : 377-378-23-33. DOI:10.1016/j.epsl.2013.07.020
Li Y X, Huang C, Hu X K, et al. 2001. The rigid and elastic-plastic model of the blocks in intro-plate and strain status of principal blocks in the continent of China. Acta Seismologica Sinica (in Chinese) , 24 (6) : 565-572.
Liu B J, Feng S Y, Ji J F, et al. 2015. Fine lithosphere structure beneath the middle-southern segment of the Tan-Lu fault zone. Chinese J. Geophys. (in Chinese) , 58 (5) : 1610-1621. DOI:10.6038/cjg20150513
Liu D W, Liu Z M, Shen X Q, et al. 2009. Study on the Seismic Risk in Centric-South Sector of the Tan-Lu Fault Zone after the MS8.0 Wenchuan Earthquake. Recent Developments in World Seismology , 4 (4) : 85-85.
Liu X X, Jiang Z S, Wu Y Q. 2012. Motion and deformation state of the Tancheng-Lujiang fault zone derived from GPS data. Earthquake (in Chinese) , 32 (4) : 1-10.
McCaffrey R. 2002. Crustal block rotations and plate coupling. //Plate Boundary Zones. Washington, DC:AGU : 101-122. DOI:10.1029/GD030p0101
McCaffrey R. 2005. Block kinematics of the Pacific-North America plate boundary in the southwestern United States from inversion of GPS, seismological, and geologic data. J. Geophys. Res. , 110 (B7) : B07401. DOI:10.1029/2004JB003307
McCaffrey R. 2009. Time-dependent inversion of three-component continuous GPS for steady and transient sources in northern Cascadia. Geophys. Res. Lett. , 36 : L07304. DOI:10.1029/2008GL036784
McCaffrey R, Qamar A I, King R W, et al. 2007. Fault locking, block rotation and crustal deformation in the Pacific Northwest. Geophys. J. Int. , 169 (3) : 1315-1340. DOI:10.1111/j.1365-246X.2007.03371.x
Qin S Q, Li G L, Xue L, et al. 2014. Analysis of the future earthquake situation for some seismic zones in the northeast China, north China and Taiwan. Progress in Geophysics (in Chinese) , 29 (2) : 535-554. DOI:10.6038/pg20140210
Scholz C H. The Mechanics of Earthquakes and Faulting. New York: Cambridge University Press, 1990 .
Scholz C H. 1998. Earthquakes and friction laws. Nature , 391 (6662) : 37-42. DOI:10.1038/34097
Thatcher W. 2007. Microplate model for the present-day deformation of Tibet. J. Geophys. Res. , 112 (B1) : B01401. DOI:10.1029/2005JB004244
Wallace L M, Stevens C, Silver E, et al. 2004. 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. Geophys. Res. , 109 (B5) : B05404. DOI:10.1029/2003JB002481
Wan G M, Tang L J, Jin Z W, et al. 2009. Progresses and problems in the study of Tancheng-Lujiang fault zone. Geological Review (in Chinese) , 55 (2) : 251-259.
Wang H, Liu M, Cao J L, et al. 2011. Slip rates and seismic moment deficits on major active faults in mainland China. J. Geophys. Res. , 116 (B2) : B02405. DOI:10.1029/2010JB007821
Wang H, Wright T J. 2012. Satellite geodetic imaging reveals internal deformation of western Tibet. Geophys. Res. Lett. , 39 : L07303. DOI:10.1029/2012GL051222
Wang K L, Wells R, Mazzotti S, et al. 2003. A revised dislocation model of interseismic deformation of the Cascadia subduction zone. J. Geophys. Res. , 108 (B1) . DOI:10.1029/2001JB001227
Wang M, Li Q, Wang F, et al. 2011. Far-field coseismic displacements associated with the 2011 Tohoku-oki earthquake in Japan observed by Global Positioning System. Chinese Sci. Bull. , 56 (23) : 2419-2424. DOI:10.1007/s11434-011-4588-7
Wang Q, Zhang P Z, Freymueller J T, et al. 2001. Present-day crustal deformation in China constrained by global positioning system measurements. Science , 294 (5542) : 574-575. DOI:10.1126/science.1063647
Wang W, Yang S M, Wang Q. 2009. Crustal block rotations in Chinese mainland revealed by GPS measurements. Earthquake Science , 22 (6) : 639-649. DOI:10.1007/s11589-009-0639-4
Wang X F, Li Z J, Chen B L, et al. The Tanlu Fault System (in Chinese). (in Chinese) Beijing: Geological Publishing House, 2000 .
Wu D M, Zhang Y M, Fang Z J., et al. 1981. On the activity of the Tancheng-Lujiang fault zone in China. Seismology and Geology (in Chinese) , 3 (4) : 15-26.
Xu C J, Li Z H, Wang H. 2002. The temporal and spatial variation characteristics of crustal deformation of active tectonic blocks in north China. Journal of Geodesy and Geodynamics (in Chinese) , 22 (2) : 33-40.
Xu J W. 1993. The Tancheng-Lujiang Wrench Fault System. Chichester:Wiley .
Yang S M, Nie Z S, Jia Z G, et al. 2011. Far-field coseismic surface displacement caused by the MW9.0 Tohoku earthquake. Geomatics and Information Science of Wuhan University (in Chinese) , 36 (11) : 1336-1339.
Yin H T, Gan W J, Huang B, et al. 2013. Study on the effects of Japan M9.0 huge earthquake on the crustal movement of Shandong area. Chinese J. Geophys. (in Chinese) , 56 (5) : 1497-1505. DOI:10.6038/cjg20130508
Yin J Y, Liu L, Wang H L, et al. 2015. The research on the impact of displacement field of Japan 3·11 earthquake on the south central section of Tanlu fault zone. Seismology and Geology (in Chinese) , 37 (3) : 928-938.
Zhang P Z. 2013. Beware of slowly slipping faults. Nat. Geosci. , 6 (5) : 323-324. DOI:10.1038/ngeo1811
Zhang P Z, Deng Q D, Zhang G M, et al. 2003. Active tectonic blocks and strong earthquakes in the continent of China. Science in China Series D:Earth Sciences, 46(S2):13-24.
Zhang P Z, Deng Q D, Zhang Z Q, et al. 2013. Active faults, earthquake hazards and associated geodynamic processes in continental China. Scientia Sinica Terrae (in Chinese) , 43 (10) : 1607-1620.
Zhang P Z, Gan W J, Shen Z K, et al. 2005. A coupling model of rigid-block movement and continuous deformation:patterns of the present-day deformation of China's continent and its vicinity. Acta Geologica Sinica , 79 (6) : 748-756.
Zhang X, Hao M, Jia P, et al. 2013. Inversion of negative dislocation for GPS horizontal movement and features of strain accumulation in main tectonic areas of China. Journal of Seismological Research (in Chinese) , 36 (1) : 1-8.
Zhang Z Q, McCaffrey R, Zhang P Z. 2013. Relative motion across the eastern Tibetan Plateau:contributions from faulting, internal strain and rotation rates. Tectonophysics , 584 : 240-256. DOI:10.1016/j.tecto.2012.08.006
Zhao B, Huang Y, Zhang C H, et al. 2015. Crustal deformation on the Chinese mainland during 1998-2014 based on GPS data. Geodesy and Geodynamics , 6 (1) : 7-15. DOI:10.1016/j.geog.2014.12.006
Zhao J, Jiang Z S, Wu Y Q, et al. 2012. Study on fault locking and fault slip deficit of the Longmenshan fault zone before the Wenchuan earthquake. Chinese J. Geophys. (in Chinese) , 55 (9) : 2963-2972. DOI:10.6038/j.issn.0001-5733.2012.09.015
Zhao J, Jiang Z S, Wu Y Q, et al. 2013. Analysis of reliability and stability of inversion result with negative dislocation model of Defnode. Journal of Geodesy and Geodynamics (in Chinese) , 33 (1) : 21-24.
Zhou C Y, Diao G L, Geng J, et al. 2013. 3-D characteristics inversion of hypocenter fault-plane of the 1668 Tancheng great earthquake. Progress in Geophysics (in Chinese) , 28 (6) : 2814-2824. DOI:10.6038/pg20130602
Zhu Z, Meng G J, Su X N, et al. 2014. Study on characteristics of tectonic blocks and faults in the north China using block-dislocation model. Journal of Geodesy and Geodynamics (in Chinese) , 34 (4) : 39-42.
邓启东. 中国活动构造图(1:4000000). 北京: 地震出版社, 2007 .
方仲景, 丁梦林, 计凤桔, 等. 1980. 郯城-庐江断裂带地震活动的地质分析. 地震地质 , 2 (4) : 39–45.
高维明, 李家灵, 孙竹友. 1980. 沂沭大陆裂谷的生成与演化. 地震地质 , 2 (3) : 11–18.
国家地震局地质研究所. 郯庐断裂. 北京: 地震出版社, 1987 .
郭良迁, 薄万举, 杨国华, 等. 2011. 华北1999-2009年水平形变应变场特征. 大地测量与地球动力学 , 31 (3) : 15–19.
郭良迁, 应绍奋. 1998. 从GPS复测结果初步研究华北北部的水平形变及构造活动. 中国地震 , 14 (4) : 11–19.
黄立人, 郭良迁. 1998. 华北北部GPS观测(1995-1996)及其结果解释. 地震地质 , 20 (4) : 423–430.
江在森, 方颖, 武艳强, 等. 2009. 汶川8.0级地震前区域地壳运动与变形动态过程. 地球物理学报 , 52 (2) : 505–518.
李延兴, 黄珹, 胡新康, 等. 2001. 板内块体的刚性弹塑性运动模型与中国大陆主要块体的应变状态. 地震学报 , 24 (6) : 565–572.
刘保金, 酆少英, 姬计法, 等. 2015. 郯庐断裂带中南段的岩石圈精 细结构. 地球物理学报 , 58 (5) : 1610–1621. DOI:10.6038/cjg20150513
刘东旺, 刘泽民, 沈小七, 等. 2009. 汶川8.0级地震后郯庐断裂带中南段地震危险性研究. 国际地震动态 , 4 (4) : 85.
刘晓霞, 江在森, 武艳强. 2012. 利用GPS资料研究郯庐带现今运动及变形状态. 地震 , 32 (4) : 1–10.
秦四清, 李国梁, 薛雷, 等. 2014. 东北、华北与台湾地区某些地震区地震趋势分析. 地球物理学进展 , 29 (2) : 535–554. DOI:10.6038/pg20140210
万桂梅, 汤良杰, 金文正, 等. 2009. 郯庐断裂带研究进展及存在问题探讨. 地质论评 , 55 (2) : 251–259.
王敏, 李强, 王凡, 等. 2011. 全球定位系统测定的2011年日本宫城MW9. 0级地震远场同震位移. 科学通报 , 56 (20) : 1593–1596.
王小凤, 李中坚, 陈柏林, 等. 郯庐断裂带. 北京: 地质出版社, 2000 .
吴大铭, 张裕明, 方仲景, 等. 1981. 论中国郯庐断裂带的活动. 地震地质 , 3 (4) : 15–26.
许才军, 李志才, 王华. 2002. 华北地区活动地块运动时空变化特征. 大地测量与地球动力学 , 22 (2) : 33–40.
杨少敏, 聂兆生, 贾志革, 等. 2011. GPS解算的日本MW9.0级地震的远场同震地表位移. 武汉大学学报:信息科学版 , 36 (11) : 1336–1339.
殷海涛, 甘卫军, 黄蓓, 等. 2013. 日本M9.0级巨震对山东地区地壳活动的影响研究. 地球物理学报 , 56 (5) : 1497–1505. DOI:10.6038/cjg20130508
尹京苑, 刘岚, 王华林, 等. 2015. 日本3·11大地震对郯庐断裂带(中南段)位移场影响的数值模拟研究. 地震地质 , 37 (3) : 928–938.
张培震, 邓起东, 张国民, 等. 2003. 中国大陆的强震活动与活动地块. 中国科学 , 33 (B04) : 12–20.
张培震, 邓起东, 张竹琪, 等. 2013. 中国大陆的活动断裂、地震灾害及其动力过程. 中国科学:地球科学 , 43 (10) : 1607–1620.
张希, 郝明, 贾鹏, 等. 2013. 全国主要构造区GPS水平运动负位错反演与应变积累特性. 地震研究 , 36 (1) : 1–8.
赵静, 江在森, 武艳强, 等. 2012. 汶川地震前龙门山断裂带闭锁程度和滑动亏损分布研究. 地球物理学报 , 55 (9) : 2963–2972. DOI:10.6038/j.issn.0001-5733.2012.09.015
赵静, 江在森, 武艳强, 等. 2013. Defnode负位错模型反演结果的可靠性和稳定性分析. 大地测量与地球动力学 , 33 (1) : 21–24.
周翠英, 刁桂苓, 耿杰, 等. 2013. 1668年郯城大地震震源断层三维特征反演. 地球物理学进展 , 28 (6) : 2814–2824. DOI:10.6038/pg20130602
朱泽, 孟国杰, 苏小宁, 等. 2014. 利用块体-位错模型反演华北地区块体运动和断裂活动. 大地测量与地球动力学 , 34 (4) : 39–42.