中国地震  2018, Vol. 34 Issue (3): 414-427
汶川8.0级地震前GPS与跨断层资料反映的运动与变形演化特征
赵静1,2, 刘杰1, 任金卫3, 江在森3, 闫伟1, 岳冲1, 苑争一1     
1. 中国地震台网中心, 北京市西城区三里河南横街5号 100045;
2. 中国地震局地质研究所(地震动力学国家重点实验室), 北京市朝阳区华严里甲1号 100029;
3. 中国地震局地震预测研究所(中国地震局地震预测重点实验室), 北京 100036
摘要:为了研究与总结2008年5月12日汶川8.0级地震前GPS与跨断层资料反映的龙门山断裂带及其周边地区的运动、构造变形、应变积累演化过程,以及汶川地震临震阶段可能的物理机制,本文综合1999~2007期GPS速度场、1999~2008年大尺度GPS基线时间序列、1985~2008年跨断层短水准等资料进行了相关分析与讨论。结果表明:①GPS速度剖面结果显示,宽达500km的川西高原在震前有明显的连续变形,而四川盆地一侧和跨龙门山断裂带基本没有变形趋势,表明震前川西高原在持续不断地为已经处于闭锁状态的龙门山断裂带提供能量积累。②GPS应变率结果显示,震前龙门山断裂带中北段的NW侧EW向挤压变形明显,变形幅度从远离断裂带较大到靠近断裂带逐渐减小,而断裂带变形微弱;龙门山断裂带西南段周边形成了显著的EW向挤压应变集中区,应变积累速率明显大于中北段。③断层闭锁程度反演结果显示,除了汶川地震的震源位置闭锁相对较弱,且西南段有大概20km宽度断层在12~22.5km深度为蠕滑状态以外,震前整条龙门山断裂基本处于强闭锁状态。④大尺度GPS基线结果显示,跨南北地震带区域的NE向基线从2005年开始普遍出现压缩转折,反映NE向地壳缩短的相对运动增强。⑤跨断层短水准场地结果显示,震前年均垂直变化速率和形变累积率很低,表明断层近场垂向活动很弱、闭锁较强。通过以上分析认为,在相对小尺度的地壳变形中,震前龙门山断裂带深浅部均处于强闭锁状态,断裂带水平与垂直变形都很微弱,这可能经历了一个缓慢的过程,而且越是临近地震的发生,微弱变形的范围可能越大;在相对大尺度的地壳变形中,震前龙门山断裂带西侧的巴颜喀拉块体东部地区经历了地壳缓慢且持续的缩短挤压变形,为龙门山断裂带应变积累持续提供了动力支持。
关键词汶川地震    GPS资料    跨断层资料    断层闭锁与变形    块体运动与变形    
Movement and Deformation Evolution Characteristics Reflected by GPS and Cross-fault Data before the MS8.0 Wenchuan Earthquake
Zhao Jing1,2, Liu Jie1, Ren Jinwei3, Jiang Zaisen3, Yan Wei1, Yue Chong1, Yuan Zhengyi1     
1. China Earthquake Networks Center, Beijing 100045, China;
2. State Key Laboratory of Earthquake Dynamics, Institute of Geology, CEA, Beijing 100029, China;
3. CEA Key Laboratory of Earthquake Prediction, Institute of Earthquake Forecasting, CEA, Beijing 100036, China
Abstract: In order to study and summarize movement, tectonic deformation and strain accumulation evolution characteristics of the Longmenshan fault and its surrounding area reflected by GPS and cross-fault data before the MS8.0 Wenchuan earthquake, and study the possible physical mechanism late in the seismic cycle of the Wenchuan earthquake, we use GPS velocity during 1999~2007, GPS baseline time series on large scale during 1999~2008 and cross-fault leveling data during 1985~2008 to make analysis and discussion. (1)GPS velocity profiles show that obvious continuous deformation across the eastern Tibet before the earthquake distributed across a zone at least 500km wide, while there is little deformation in the Sichuan Basin and Longmenshan fault zone, which means that the eastern Tibet provides energy accumulation for the locked Longmenshan fault zone continuously.(2)GPS strain rates show that the east-west compression deformation is large northwest of the mid-northern segment of the Longmenshan fault zone, and deformation amplitude decreases gradually from far field to near fault zones, and there is little deformation in fault zones. The east-west compression deformation is significant surrounding the southwestern segment of the Longmenshan fault zone, and strain accumulation rate is larger than that of the mid-northern segment. (3)Fault locking indicates nearly whole Longmenshan fault is locked before the earthquake except that the hypocenter of the earthquake is weaker locked, and a 20km width patch in the southwestern segment between 12km to 22.5km depth is in creeping state. (4)GPS baseline time series in the northeast direction on large scale become compressive generally from 2005 in the north-south seismic zone, which reflects that relative compression deformation enhances. (5)The cross-fault leveling data show that annual vertical change rate and deformation trend accumulation rate in the Longmenshan fault zone are little, which indicates that vertical activity near the fault is very weak and the fault is tightly locked. By analyzing GPS and cross-fault data before the Wenchuan earthquake, we consider that the Longmenshan fault is tightly locked from the surface to the deep, and the horizontal and vertical deformation are weak surrounding the fault in relative small scale crustal deformation. The process of weaken deformation may be slow, and when close to large earthquake, the weaken deformation area may be larger. Continuous and slow compression deformation across the eastern Tibet before the earthquake provides dynamic support for strain accumulation in the Longmenshan fault zone in relative large scale crustal deformation.
Key words: The Wenchuan earthquake     GPS data     Cross-fault data     Fault locking and deformation     Block movement and deformation    
0 引言

2008年5月12日汶川8.0级地震发生在巴颜喀拉块体与华南块体交界的NE向龙门山断裂带上,地震基本为沿着断裂带的NE向单侧破裂。汶川地震是有历史记载以来首例发生在大陆内部低滑动速率断裂上的高角度逆冲型8级强震,其孕育机制和发生机理与其他特大逆冲型地震可能是不同的(张培震等,2009),而要认识汶川地震的孕育过程,仅研究其发震断裂——龙门山断裂带是不够的,巨大地震孕育过程涉及更大的时空尺度,因此需要对与发震断裂相关、相邻的地质构造单元在震前不同时段、不同尺度的运动与变形状态进行观测和研究,才能获得完整、可靠的信息(中国地震局监测预报司,2009),这包括垂直断裂带的GPS速度场剖面、地表大范围的GPS应变率场、断层面上闭锁程度分布、大尺度GPS基线时间序列、近场跨断层短水准等。至2018年5月12日,汶川地震已发生10年,在这期间许多科学工作者利用震前GPS资料与跨断层资料对龙门山断裂带的闭锁与变形特征、巴颜喀拉块体与华南块体的运动与变形演化过程等进行了研究,Chen等(2000)张培震等(20082009)Hubbard等(2009)江在森等(2009)杜方等(2009)Zhang等(2010)赵静等(2012)Wu等(2015)利用GPS速度剖面方法对龙门山断裂带远、近场变形特征进行了相关分析;江在森等(2009)武艳强等(2011)赵静等(2012)Wei等(2015)Wu等(2015)对龙门山断裂带及其周边地区的应变率场进行了相关分析;Shen等(2005)Meade(2007)赵静等(2012)建立块体与断层模型,对龙门山断裂带运动特征、闭锁程度和滑动亏损空间分布进行了相关研究;江在森等(2009)利用大尺度GPS基线对南北地震带大范围地壳变形进行了相关分析;焦青等(2008)苏琴等(2009)程万正等(2011)Bo等(2011)利用跨断层短水准资料对龙门山断裂带垂向运动进行了相关分析;杨国华等(2006)陈祖安等(2011)陈运泰等(2013)徐晶等(2017)分析了远大震和周边强震对龙门山断裂带及附近地区的影响。除此之外,前人还在其他方面开展了较丰富的研究工作。上述研究为我们进一步认识汶川地震前龙门山断裂带近场与远场的变形特征、研究汶川地震孕育过程等提供了帮助,但是这些研究大部分只涉及了一种或几种方法,没有利用多种资料、方法将地表和断层深部运动与变形、断层水平方向和垂直方向运动与闭锁等进行综合研究,因此缺乏针对龙门山断裂带远近场、深浅部、三维立体的综合性分析。本文的目的是对这方面的10年研究成果进行汇总,使读者了解汶川地震前GPS与跨断层资料的变化特征和目前的认识。

GPS和跨断层观测数据可以有效监测断裂带及其周边区域不同尺度地壳三维动态变形特征,而断层深浅部运动与闭锁作用很大程度上影响着地壳变形,因此本文将主要依据1999~2007期GPS速度场,综合前人跨断层速度剖面结果分析断裂带远、近场的运动变形特征;综合前人最小二乘配置方法计算结果分析周边区域的应变率场,采用DEFNODE负位错反演程序估算汶川地震前龙门山断裂带的闭锁程度和滑动亏损分布;综合前人大尺度GPS基线时间序列结果分析南北地震带地壳运动特征;利用1985~2008年跨断层短水准资料分析龙门山断裂带不同场地年均垂直变化速率、断层形变累积率空间分布。然后利用上述相关研究结果综合分析汶川地震前龙门山断裂带近场和远场的水平与垂直变形演化特征、地表与断层深部运动与变形特征,探讨汶川地震的孕育过程和发震机理,以期对大陆内部低滑动速率断裂上的高角度逆冲型8级强震有更全面、深入的认识。

1 跨断层GPS速度剖面

利用横跨断裂带的GPS速度剖面可以分析断裂带的远、近场变形特征(Savage et al,1973Meade et al,2005),汶川地震后多位学者(张培震等,20082009江在森等,2009杜方等,2009Zhang et al,2010赵静等,2012Wu et al,2015)利用该方法对龙门山断裂带进行了相关分析。图 1为垂直于龙门山断裂带的GPS速度剖面,剖面NW侧数据结果显示,平行于龙门山断裂带的水平右旋走滑分量和垂直于龙门山断裂带的水平缩短分量都揭示了宽达500km的川西高原震前遭受着明显的连续变形,其右旋走滑速率为9~11mm/a,水平缩短速率为4~5mm/a;川西高原中速度剖面大体上的线性梯度表明其变形和应变是准均匀的(张培震等,2009)。剖面SE侧数据结果显示,四川盆地一侧GPS站速度都在0值附近,表明四川盆地内部无明显变形趋势,且跨整个龙门山断裂带的GPS站速度与四川盆地内部一致,基本看不到变形趋势(江在森等,2009)。以上结果表明川西高原震前经历着强烈的变形,其在变形幅度和空间范围上远大于龙门山断裂带和四川盆地(张培震等,2009)。

图 1 龙门山地区活动构造与汶川地震前GPS速度剖面(据Zhang等(2010)) (a)龙门山地区活动构造与GPS速度场(相对华南块体);(b)汶川震前横跨龙门山断裂带GPS速度剖面
2 GPS应变率场

由GPS速度场计算得到的应变率场可以反映地壳上部的张压和剪切等变形特征,因此汶川地震前龙门山断裂带及其附近地区的应变率场结果能够反映地震孕育的部分特征。选取25°~34°N、99°~107°E范围内的GPS水平速度场数据,采用最小二乘配置方法(武艳强等,2009江在森等,2010Wu et al,2011)计算得到了该区域的主应变率、EW向应变率和最大剪应变率的分布(图 2)。主应变率结果(图 2(a))表明龙门山断裂带西南段的应变积累速率远大于中北段,具有明显的分区特征,整条断裂带表现为挤压兼右旋剪切变形特征。断裂带中北段的NW侧变形幅度从远离断裂带到靠近断裂带逐渐减小,表明川西高原的长期变形在持续不断地为闭锁的龙门山断裂带提供应变积累,并且在震前较长时期已经积累了很高的应变能,而龙门山断裂带由于本身强闭锁状态和四川盆地的强烈阻挡,只发生微弱的震间变形,使得有GPS观测以来无法观测到汶川地震前的断裂带显著变形(赵静等,2012)。断裂带中北段的SE侧变形很微弱,表明作为支撑单元的四川盆地在汶川地震前对川西高原的挤压起到了稳定的阻挡作用,这也是造成龙门山断裂带闭锁和高应变积累的必要条件(张培震等,2009)。EW向应变率结果(图 2(b))表明受巴颜喀拉块体东向运动影响,龙门山断裂带西南段周边区域形成了显著的EW向挤压应变集中区;龙门山断裂带中北段的NW侧挤压变形明显,而断裂带的SE侧则属于EW向拉张变形区,该结果表明龙门山断裂带作为EW向变形的分界带,具备长期积累挤压应变的背景(赵静等,2012),而汶川地震就发生在EW向强烈挤压变形的边缘区域。最大剪应变率结果(图 2(c))表明鲜水河-安宁河-则木河-小江断裂带上最大剪应变率较高(计算时断层滑动速率没有扣除,因此断层蠕滑可能会影响应变率值计算),而在龙门山断裂中北段剪应变并不明显(Wu et al,2015),表明汶川震前龙门山断裂中北段的右旋剪切运动也较弱,断层处于闭锁状态。此外,从整个中国大陆西部的EW向变形特征(Wei et al,2015Wu et al,2015)来看,青藏块体在92.5°以西地区主要表现为EW向拉张,反映了青藏高原地壳物质的东向流动;在92.5°E以东地区主要表现为EW向挤压,一直持续到龙门山断裂带附近,反映了青藏高原地壳物质东向流动受到了华北和华南块体的阻挡作用,形成了宽达900多千米的挤压区域,而龙门山断裂带作为挤压区域的边界,一直接受着来自西部青藏高原大范围挤压应变能的积累。

图 2 1999~2007期GPS应变率场结果(赵静等,2012Wu et al,2015) (a)主应变率分布;(b)EW向应变率分布;(c)最大剪应变率分布

汶川地震前川滇地区EW向应变率场的三期结果(图 3)反映了龙门山断裂带周边地区EW向挤压变形的演化特征。1999~2001年(图 3(a)),EW向挤压变形高值区的面积较小、最高值区位于鲜水河断裂带;龙门山断裂带西南段挤压变形较强,而中北段挤压变形微弱。2001~2004年(图 3(b)),EW向挤压变形高值区的量值和面积均有所增加,最大挤压变形仍然位于鲜水河断裂带并向东移动;龙门山断裂带挤压变形范围有所增加,西南段挤压变形进一步增强,中北段附近挤压变形较强。2004~2007年(图 3(c)),EW向挤压变形的量值明显增大,最大挤压变形继续往东迁移至龙门山断裂带与鲜水河断裂带交界处;在EW向应变率场逐渐东向迁移的过程中,龙门山断裂带的EW向挤压变形逐渐增强,到2004~2007年达到最高(武艳强等,2011)。

图 3 川滇地区EW向应变率场结果 (a)1999~2001期;(b)2001~2004期;(c)2004~2007期
3 断层闭锁

在前面分析跨断层GPS速度剖面、GPS应变率场的基础上,我们建立了块体与断层物理模型对龙门山断裂带在汶川地震前的闭锁程度和滑动亏损空间三维分布等进行了相关反演计算,试图从断层深部运动与变形角度对汶川地震有更深一步的认识。所用DEFNODE负位错反演程序的基本原理和使用方法见参考文献(McCaffrey,20022005Savage et al,2001赵静等,20122015),该程序假定块体内部点的运动为块体旋转、块体内部均匀应变和由于断层闭锁而引起的地表弹性变形之和,因此该程序能够利用GPS数据,采用网格搜索和模拟退火方式反演计算块体旋转欧拉极、块体内部均匀应变、块体边界断层闭锁程度和滑动亏损速率三维分布等,其中,断层闭锁程度为断层滑动亏损速率与断层长期滑动速率的比值,比值为0表示断层完全蠕滑,比值为1表示断层完全闭锁,比值介于0~1之间,表示断层部分闭锁,并存在部分蠕滑运动。我们首先对该程序的反演可靠性和稳定性进行了验证,结果表明在不同断层倾角条件下反演结果均具有高度可靠性;权值f对反演结果影响微弱,反演结果也具有较好的稳定性;沿经度、纬度方向均为1°间隔分布的GPS速度场数据基本能够约束发震断层的震间变形形态(赵静等,2013)。

所用数据为汶川地震前的1999~2007期GPS速度场(王敏研究员提供),共有186个测站数据参与反演,因此GPS的空间分辨率应该可以约束汶川地震发震断层的震间变形形态。考虑到龙门山断裂带与鲜水河断裂带、安宁河-则木河-小江断裂带在青藏高原东缘构成了交叉口,三条断裂带之间的运动与变形相互影响,因此在建立块体与断层物理模型时,将整个研究区域以上述三条断裂带为边界,划分为巴颜喀拉地块、川滇地块和华南地块。反演所得龙门山断裂带闭锁程度和滑动亏损空间分布情况见图 4图 4(a)结果显示:整条龙门山断裂基本处于强闭锁状态,只有汶川地震的震源位置处闭锁相对较弱,闭锁程度为0.88左右;西南段有大概20km宽度断层在12~22.5km深度为蠕滑状态,为整条断裂带上唯一不闭锁的位置,而芦山地震则发生在强闭锁至不闭锁的过渡部位(赵静等,2018)。图 4(b)结果显示:平行龙门山断裂的滑动亏损速率均为右旋亏损速率,北东段亏损速率最大,为5.8~6.0mm/a;中南段完全闭锁段落的亏损速率为5.7~5.8mm/a,较北东段稍有减弱。图 4(c)结果显示:垂直龙门山断裂的滑动亏损速率均为挤压亏损速率,自北东端至南西端滑动亏损速率逐渐增大,其中,北东段亏损速率最小,为0.6~1.1mm/a;中南段完全闭锁段落的亏损速率为2.0~3.6mm/a。

图 4 汶川地震前龙门山断裂带闭锁程度与滑动亏损空间分布(据赵静等(2018)) (a)断层闭锁程度;(b)平行断层滑动亏损速率;(c)垂直断层滑动亏损速率
4 大尺度GPS基线

单个GPS站的时间序列会受到参考框架、季节性天气等影响,而这些影响在一定程度上是相关的,因此需要找出一些具有物理意义、受共有误差影响较小的参数来描述地壳变形信息,提取两个GPS站间基线的时间序列与参考框架相关程度较小且可削弱一些同源误差影响,能够凸显地壳变形信息,而采用大地线弧长作为基线时间序列又避免了GPS法向观测误差(其非构造影响更突出)的影响。通过滤波可削弱多种年周期及更短频段的非构造因素的影响,保留相对低频的具有稳定性的趋势变化,所反应的地壳运动信息可能更为可靠(江在森等,2009),因此提取两个GPS站间大地线弧长低频信息的方法在地壳变形分析中较为常用。汶川地震后,江在森等(2009)对中国大陆GPS基准站连续观测数据进行了处理,结果表明跨南北地震带区域的NE向GPS基线时间序列自2005年以来出现长趋势普遍转折现象(图 5 (a)),反映南北地震带NE向地壳缩短的相对运动增强、青藏地块相对华南地块的NE向运动增强,其中与西南地区的KMIN(昆明)、XIAG(下关)、LUZH(泸州)等站相关的北东或北北东基线均显示,2004年12月26日苏门答腊MW9.0地震后,在2004~2005年先相对伸长,而后在2005~2006年转为缩短挤压变化,如KMIN(昆明)-YANC(盐池)基线(图 5 (b)),表明西南地区受苏门答腊地震同震影响明显(杨国华等,2006)。青藏块体和华南块体内部的NE向基线并没有出现缩短增强的趋势。同时,LUZH(泸州)-LHAS(拉萨)-YANC(盐池)GPS基准站构成的三角形变形单元应变参数时间序列结果(江在森等,2009)表明,南北地震带大尺度地壳变形的变化为近EW向地壳缩短和NE向右旋剪切活动的增强,这与龙门山断裂带发震所需要的构造动力一致,因此苏门答腊地震可能促进了汶川地震的发生。

图 5 南北地震带及其附近区域部分NE向GPS基线时间序列去线性后曲线图(江在森等,2009) (a)LHAS(拉萨)-YANC(盐池);(b)KMIN(昆明)-YANC(盐池)
纵坐标为基线长度变化量;细线为最小二乘配置结果;粗线为小波(DB6-a9)趋势部分;图右下标出了基线长期线性速率
5 跨断层短水准结果

图 6为汶川地震前龙门山断裂带跨断层短水准场地分布与高差变化,短水准场地包括灌县、七盘沟、耿达、双河、蒲江等5个场地,其水准数据时间跨度为1985~2008年。综合分析认为,汶川地震前5个场地跨断层形变资料基本没有出现明显的短期异常(苏琴等,2009Bo et al,2011)。

图 6 汶川地震前龙门山断裂带跨断层短水准场地分布(a)与高差变化(b) 北边的灌县、七盘沟、耿达场地以圆点表示;中间双河场地以方块表示;南边的蒲江以菱形表示

图 6可以看出:耿达A-B测线在震前出现了大幅拉张异常,但经异常核实认为该变化为A点受环境因素影响产生下沉,且异常出现时段与当地修建房屋时段非常吻合,因此不作为地震前兆异常(苏琴等,2009);灌县3-4测线在1992~1994年出现显著拉张变化,最大幅度达到2.47mm,1994年后断层处于弱活动状态,2005~2008年有一个先挤压后拉张的微弱变化;蒲江1-2测线自1987年以后呈现持续拉张状态,1993年9月~1995年9月出现拉张速率增强又减弱的异常过程(焦青等,2008),1995~2002年曲线恢复走势,表明该观测断层段张性活动恢复,2002~2004年曲线拉平,水准变化明显减弱,断层活动消失,而从2004年底开始,水准拉张变化再次增强,2006年下半年由拉张转挤压,之后汶川地震发生(程万正等,2011);七盘沟和双河场地水准年变周期较为规律,但同时呈现一定的压性变化。

为进一步研究汶川地震前跨断层短水准变化,我们计算了水准年均形变速率(图 7),从几个场地的年均变化速率结果可以看出,整个龙门山断裂带上断层的年均变化速率都在0附近波动变化,说明整个断层垂向运动很弱,处于强闭锁的状态,这与1999~2007期GPS反演所得断层闭锁程度结果是一致的。蒲江场地和灌县场地水准年均形变速率分别在2006年和2007年出现了小幅度的变化,可能表明断层的活动性在汶川震前稍有增强。但是纵观龙门山断裂带上所有断层形变资料不难看出,汶川地震前基本没有出现明显的短期突跳或速率改变异常。

图 7 汶川地震前龙门山断裂带跨断层短水准年均垂直变化速率

为研究汶川地震前龙门山断裂带周边地区的跨断层短水准垂直变化空间分布结果,我们计算了断层垂直形变累积率Dc(图 8),Dc为断层形变长期平均年变速率与平均年变幅度之比(无量纲),其绝对值范围在0~1之间,越趋于1,表示曲线线性趋势越好且波动变化越小;越趋于0,表示曲线往复变化或曲线平稳且无形变和应变累积(岳冲等,2017)。图 8结果显示:鲜水河断裂带上的侏倭、沟普、安顺场场地Dc绝对值一直趋于1,表明其长期处于形变和应变累积过程,而其他场地的形变和应变累积速率较低,综合鲜水河断裂带8处场地结果可以看出,断裂带北西段闭锁较弱,而南东段闭锁较强。龙门山断裂带上耿达、七盘沟、双河场地具有很小的形变累积率,结合图 7可以看出这三处场地的年均垂直变化速率曲线的波动性即扰动性并不强,表明这三处场地活动很弱;蒲江和灌县场地的形变累积率稍大,结合图 67结果可以看出,灌县场地形变累积率稍大主要是由1992~1996年观测数据引起的。综合龙门山断裂带5处场地结果,可以看出整个断裂带的形变累积率很低,表明断层没有明显的速率改变、处于强闭锁状态。

图 8 汶川地震前断层形变累积率空间分布等值线(1990~2008年结果)
6 讨论与结论

(1) 利用汶川地震同震位移除以GPS方法获得的震间断层滑动速率,得到汶川地震的复发间隔为2500~5000年(赵静等,2012)、3190~5952年(张培震等,2008)、2000~10000年(Burchfiel et al,2008),但是十年尺度GPS结果获得的断层滑动速率是否能够代表过去几千年甚至更长时间断层的滑动速率,是计算地震复发间隔的关键所在(Yin,2010)。像喜马拉雅断层这样的大断裂,远离断层500km左右还能够积累弹性应变(Feldl et al,2006),而断层近场GPS观测不到明显的变形,当积累的能量释放时却能够促使断层上大地震的发生。龙门山断裂带与喜马拉雅断裂带有一定的相似性,均为大陆内部逆冲断裂带,只不过龙门山断裂带为仰冲断裂带(主动盘为上盘),喜马拉雅断裂带为俯冲断裂带(主动盘为下盘),因此除了青藏高原东缘的龙门山断裂带以外,巴颜喀拉块体东部应该有一个大范围应变积累区域,而跨断层GPS速度剖面结果就显示宽达500km的整个川西高原震前经历着空间上连续分布的强烈变形。地震地质研究结果给出的汶川地震前龙门山断裂带滑动速率比GPS结果给出的速率要高,因此利用汶川地震同震位移除以地震地质方法获得的断层滑动速率,得到汶川地震的复发间隔为2233~4167年,比GPS结果获得的复发间隔要短(张培震等,2008)。

(2) 大地震释放的弹性应变来源于震前地壳运动造成的长期应变积累,汶川地震可能经历了几千年的弹性应变能积累过程,而跨断层资料的观测时间只有不到30年、GPS资料的观测时间只有10年,都处于汶川地震周期的最后阶段,因此可能在该阶段龙门山断裂带附近地壳弹性应变积累已经趋于极限,正处于一个十分缓慢积累甚至已经无法积累的状态,这可能是导致前兆异常特别是近断层异常不易突出显现的原因,也是导致我们在中长期尺度上低估了其地震危险性的重要原因之一。汶川地震后,我们认识到断层在低运动和变形速率下并不一定代表低的潜在地震危险性,特别是逆冲型断层可能反映了其正处在强闭锁和应变能高度积累的状态。由于汶川地震前龙门山断裂带北西盘没有布设GPS连续观测站,因此我们对震前断层附近的地壳连续变形并没有获得较好的观测,也无法用连续的GPS资料分析断层附近连续变形的演化特征。汶川地震后中国地震局地震预测研究所江在森研究员的研究团队认识到,未发生破裂的龙门山断裂带南段依然有发生强震的危险性,因此他们在南段周边布设了10个GPS连续观测站,并成功捕捉到了芦山7.0级地震前龙门山断裂带南段的连续变形演化特征(刘晓霞等,2015),使我们对大陆内部逆冲型断层的震前变形和发震机理等有了更进一步的认识。

(3) 在地震学研究中,震间变形因涉及时间窗口长和各种干扰因素的影响,已成为具有挑战性的问题之一。在汶川地震前大致半个世纪内,龙门山断裂带周边区域发生了多次6~7级强震,主要包括1976年8月16日松潘7.2级、8月22日松潘6.7级、8月23日松潘7.2级地震序列,1941年6月12日泸定-天全一带6.0级地震和1970年2月24日大邑6.2级地震。有关大地震同震效应和震后效应引起的库仑应力变化研究结果显示,松潘地震序列对2008年汶川8.0级地震、2013年芦山7.0级地震和整个龙门山断裂带影响微弱(徐晶等,2017)。1941年泸定-天全一带6.0级地震和1970年大邑6.2级地震所释放的地震波能量太小,仅为MW6.8(汶川地震和芦山地震之间未破裂段落的震级估值)地震所释放地震波能量的6%~13%,不足以改变龙门山断裂带西南段危险断层的地震估值(陈运泰等,2013)。地震的影响既有同震相变,也有震后蠕变调整,虽然上述研究指出这些强震对汶川地震的影响可能较弱,主要因素仍是汶川地震本身能量的积累,但本文给出的龙门山断裂带及周边地区的震间变形仍然是包括这些地震影响的,这也是不同地震之间变形差异的一个原因。

(4) 横跨断裂带的GPS速度剖面结果显示,不论是垂直于龙门山断裂的水平缩短分量还是平行于龙门山的水平右旋走滑分量,都揭示出宽达500km的川西高原在震前遭受着明显的、连续的变形;而四川盆地一侧GPS站速度和跨龙门山断裂处的GPS站速度都在0值附近,基本看不到变形趋势,表明汶川地震前川西高原在持续不断地为已经处于闭锁状态的龙门山断裂带提供能量积累。

(5) GPS主应变率结果显示,龙门山断裂带西南段的应变积累速率明显大于中北段,整条断裂带表现为挤压兼右旋剪切变形特征。龙门山断裂带中北段的NW侧变形幅度从远离断裂带较大到靠近断裂带逐渐减小,而断裂带由于本身强闭锁状态,不发生或只发生微弱的震间形变。GPS面应变率结果显示,龙门山断裂带西南段形成了显著的EW向挤压应变集中区;以龙门山断裂带中北段为界,断裂带的NW侧挤压变形明显,而断裂带的SE侧则属于EW向拉张变形区。GPS剪应变率结果显示,鲜水河-安宁河-则木河-小江断裂带上最大剪应变率很高,而在龙门山断裂中北段剪应变并不明显。GPS应变率多种结果表明龙门山断裂带作为巴颜喀拉块体和华南块体运动与变形的分界带,具备长期积累挤压和右旋剪切应变的背景。

(6) 闭锁程度和滑动亏损三维分布反演结果显示,汶川地震前整条龙门山断裂带基本处于强闭锁状态,只有汶川地震的震源位置处闭锁相对较弱(闭锁程度为0.88左右),西南段有大概20km宽度断层在12~22.5km深度为蠕滑状态,而芦山地震则发生在强闭锁至不闭锁的过渡部位。平行龙门山断裂带的滑动亏损速率均为右旋亏损速率,北东段亏损速率为5.8~6.0mm/a,中南段完全闭锁段落的亏损速率为5.7~5.8mm/a;垂直龙门山断裂带的滑动亏损速率均为挤压亏损速率,北东段亏损速率为0.6~1.1mm/a,中南段完全闭锁段落的亏损速率为2.0~3.6mm/a。此外,随机选取的2组只包含80%GPS数据反演计算所得结果与原始结果基本一致,仅在西南段蠕滑位置处稍有差异;权值f由最优值1.74调整为1.64和1.84后,反演计算所得结果与原始结果基本一致,表明反演结果具有较好的可靠性和稳定性。

(7) 大尺度GPS基线结果显示,跨南北地震带区域的NE向GPS基线时间序列从2005年以来长趋势出现普遍转折,反映NE向地壳缩短的相对运动增强,其中与西南地区站点相关的NE或NNE基线显示,2004年苏门答腊MW9.0地震后,在2004~2005年先相对伸长,而后在2005~2006年转为缩短挤压变化。GPS基准站构成的三角形变形单元应变参数时间序列结果显示,南北地震带大尺度地壳变形的变化为近EW向地壳缩短和NE向右旋剪切活动的增强,这与龙门山断裂带发震所需要的构造动力一致。

(8) 沿龙门山断裂带跨断层短水准场地(灌县、七盘沟、耿达、双河、蒲江场地)相关解算结果表明,汶川地震前观测原始曲线和年均垂直变化速率基本没有出现明显的短期突跳或速率改变异常;整条龙门山断裂带的形变累积率很低,表明断层近场活动很弱,处于垂向活动强闭锁状态。GPS水平速度场反演的断层面强闭锁状态和跨断层短水准资料反映的断层垂向强闭锁状态相一致,表明汶川8.0级逆冲型强震发生前,无论断层深浅部还是断层水平向与垂直向运动,都处于强闭锁状态,不能够自由运动。

(9) 通过以上汶川8.0级地震前GPS与跨断层资料的相关分析,可以看出在相对小尺度的地壳变形中,汶川地震前龙门山断裂带深浅部均处于强闭锁状态,断裂带水平与垂直变形都很微弱,从超过2000年的复发间隔来看,这种现象可能已经经历了一个长时间的缓慢过程,而且越是临近地震的发生,微弱变形的范围可能越大。同时在相对大尺度的地壳变形中,汶川地震前龙门山断裂带西侧的巴颜喀拉块体东部地区经历了地壳缓慢且持续的缩短挤压变形,为龙门山断裂带应变积累持续提供了动力支持。2001年昆仑山口西8.1级地震的发震断层发生左旋走滑错动,引起东昆仑断裂带南侧巴颜喀拉块体进一步东扩和一定规模的变形,并在东侧受到四川盆地的阻挡,使得龙门山断裂带库仑破裂应力增大,促进了其挤压应变能积累增强(江在森等,2009陈祖安等,2011)。2004年苏门答腊MW9.0地震有利于印度板块向北运动,而川滇和华南地块南部向西南方向运动加快(杨国华等,2006),同时造成南北地震带大尺度NE向地壳缩短的相对运动明显增强,这种大尺度的影响与龙门山断裂带发生右旋剪切破裂所需要的动力作用一致,因此苏门答腊地震可能会促进龙门山断裂带右旋剪切应变能积累增强(江在森等,2009)。

致谢: 感谢审稿专家提出的中肯意见和建议。
参考文献
陈运泰、杨智娴、张勇等, 2013, 从汶川地震到芦山地震, 中国科学:地球科学, 43(6): 1064-1072.
陈祖安、林邦慧、白武明等, 2011, 2001年8.1级昆仑山大震破裂过程及对2008年汶川8.0级大震孕育发生影响的研究, 地球物理学报, 54(1): 108-120.
程万正、官致君、苏琴等, 2011, 汶川MS8.0地震前四川地区前兆异常及其统计分析, 地震学报, 33(3): 304-318. DOI:10.3969/j.issn.0253-3782.2011.03.004
杜方、闻学泽、张培震等, 2009, 2008年汶川8.0级地震前横跨龙门山断裂带的震间形变, 地球物理学报, 52(11): 2729-2738.
江在森、方颖、武艳强等, 2009, 汶川8.0级地震前区域地壳运动与变形动态过程, 地球物理学报, 52(2): 505-518.
江在森、刘经南, 2010, 应用最小二乘配置建立地壳运动速度场与应变场的方法, 地球物理学报, 53(5): 1109-1117. DOI:10.3969/j.issn.0001-5733.2010.05.011
焦青、杨选辉、许丽卿等, 2008, 汶川8.0级地震前后龙门山断裂活动特征浅析, 大地测量与地球动力学, 28(4): 7-11, 37.
刘晓霞、武艳强、江在森等, 2015, GPS观测揭示的芦山MS7.0地震前龙门山断裂带南段变形演化特征, 中国科学:地球科学, 45(8): 1198-1207.
苏琴、朱航、杨永林, 2009, 耿达短水准异常与汶川8.0级地震, 大地测量与地球动力学, 29(增刊Ⅰ): 103-105.
武艳强、江在森、杨国华等, 2009, 利用最小二乘配置在球面上整体解算GPS应变场的方法及应用, 地球物理学报, 52(7): 1707-1714.
武艳强、江在森、杨国华等, 2011, 汶川地震前GPS资料反映的应变率场演化特征, 大地测量与地球动力学, 31(5): 20-25, 29.
徐晶、邵志刚、刘静等, 2017, 基于库仑应力变化分析巴颜喀拉地块东端的强震相互关系, 地球物理学报, 60(10): 4056-4068.
杨国华、江在森、王敏等, 2006, 印尼地震对我国川滇地区地壳水平活动的影响, 大地测量与地球动力学, 26(1): 9-14, 20.
岳冲、闫伟、李晓帆等, 2017, 四川断层活动性分析及汶川、芦山地震活动关联性研究, 大地测量与地球动力学, 37(9): 888-892.
张培震、闻学泽、徐锡伟等, 2009, 2008年汶川8.0级特大地震孕育和发生的多单元组合模式, 科学通报, 54(7): 944-953.
张培震、徐锡伟、闻学泽等, 2008, 2008年汶川8.0级地震发震断裂的滑动速率、复发周期和构造成因, 地球物理学报, 51(4): 1066-1073.
赵静、江在森、牛安福等, 2015, 川滇菱形块体东边界断层闭锁程度与滑动亏损动态特征研究, 地球物理学报, 58(3): 872-885.
赵静、江在森、武艳强等, 2012, 汶川地震前龙门山断裂带闭锁程度和滑动亏损分布研究, 地球物理学报, 55(9): 2963-2972.
赵静、江在森、武艳强等, 2013, Defnode负位错模型反演结果的可靠性和稳定性分析, 大地测量与地球动力学, 33(1): 21-24.
赵静、任金卫、江在森等, 2018, 龙门山断裂带西南段闭锁与变形特征, 地震研究, 41(2): 216-225.
中国地震局监测预报司, 2009, 汶川8.0级地震科学研究报告, 205. 北京: 地震出版社.
Bo W J, Yang G H, Zhan W, et al, 2011, Preparatory mechanism of MS8.0 Wenchuan earthquake evidenced by crust-deformation data, Geod Geodyn, 2(2): 23-28. DOI:10.3724/SP.J.1246.2011.00023
Burchfiel B C, Royden L H, van der Hilst R D, et al, 2008, A geological and geophysical context for the Wenchuan earthquake of 12 May 2008, Sichuan, People's Republic of China, GSA Today, 18(7): 4-11. DOI:10.1130/GSATG18A.1
Chen Z, Burchfiel B C, Liu Y, et al, 2000, Global Positioning System measurements from eastern Tibet and their implications for India/Eurasia intercontinental deformation, J Geophys Res, 105(B7): 16215-16227. DOI:10.1029/2000JB900092
Feldl N, Bilham R, 2006, Great Himalayan earthquakes and the Tibetan Plateau, Nature, 444(7116): 165-170. DOI:10.1038/nature05199
Hubbard J, Shaw J H, 2009, Uplift of the Longmen Shan and Tibetan Plateau.and the 2008 Wenchuan(M=7.9)earthquake, Nature, 458(7235): 194-197. DOI:10.1038/nature07837
McCaffrey R, 2002, Crustal block rotations and plate coupling, In: Stein S, Freymueller J T, Plate Boundary Zones, 101~122, Washington: American Geophysical Union. http://onlinelibrary.wiley.com/doi/10.1029/GD030p0101/summary
McCaffrey R, 2005, Block kinematics of the Pacific-North America plate boundary in the southwestern US from inversion of GPS. seismological, and geologic data, J Geophys Res, 110(B7).
Meade B J, Hager B H, 2005, Block models of crustal motion in southern California constrained by GPS measurements, J Geophys Res, 110(B3). DOI:10.1029/2004JB003209
Meade B J, 2007, Present-day kinematics at the India-Asia collision zone, Geology, 35(1): 81-84. DOI:10.1130/G22924A.1
Savage J C, Burford R O, 1973, Geodetic determination of relative plate motion in central California, J Geophys Res, 78(5): 832-845. DOI:10.1029/JB078i005p00832
Savage J C, Gan W J, Svarc J L, 2001, Strain accumulation and rotation in the eastern California Shear Zone, J Geophys Res, 106(B10): 21995-22007. DOI:10.1029/2000JB000127
Shen Z K, Lü J N, Wang M, et al, 2005, Contemporary crustal deformation around the southeast borderland of the Tibetan Plateau, J Geophys Res, 110(B11).
Wei W X, Jiang Z S, Wu Y Q, 2015, Cognitions and questions regarding crustal deformation and location forecasts of strong earthquakes, Geod Geodyn, 6(3): 210-219. DOI:10.1016/j.geog.2015.03.006
Wu Y Q, Jiang Z S, Yang G H, et al, 2011, Comparison of GPS strain rate computing methods and their reliability, Geophys J Int, 185(2): 703-717. DOI:10.1111/j.1365-246X.2011.04976.x
Wu Y Q, Jiang Z S, Zhao J, et al, 2015, Crustal deformation before the 2008 Wenchuan MS8.0 earthquake studied using GPS data, J Geodyn, 85: 11-23. DOI:10.1016/j.jog.2014.12.002
Yin A, 2010, A special issue on the great 12 May 2008 Wenchuan earthquake(MW7.9):observations and unanswered questions, Tectonophysics, 491(1~4): 1-9.
Zhang P Z, Wen X Z, Shen Z K, et al, 2010, Oblique.high-angle, Listric-reverse faulting and associated development of strain:the Wenchuan earthquake of May 12, 2008, Sichuan, China, Annu Rev Earth Planet Sci, 38: 353-382. DOI:10.1146/annurev-earth-040809-152602
汶川8.0级地震前GPS与跨断层资料反映的运动与变形演化特征
赵静, 刘杰, 任金卫, 江在森, 闫伟, 岳冲, 苑争一