-
断层在形成之前随着应力作用聚集大量的能量,在形成过程中释放大量的能量,使得地壳发生形变与错断,直接或间接造成地震[1]。强烈的破坏性地震往往都形成了明显的地表破裂(地震形变带或地震断层)。实际监测与研究均表明,地壳变化主要是以断层或断裂系统错动变形为核心的形变活动[2-16],绝大多数大震就发生在这样的空间部位上[17],例如2008年汶川大地震[18-21]。如今,地震预测的诸多问题都是通过研究断层活动与形变规律展开的,断层形变活动的监测与分析对地震研究有着重要意义。
基于GPS、InSAR等观测技术的广泛应用,我们能够较容易地了解到断层及临近区域地壳运动与形变情况,随着GPS观测资料不断积累,人们能够不断深入地研究断层或断层系统的活动习性、彼此关系、时空差异及相应的物理含义[22-25]。
1973年,Savage等为了描述板块间相对运动及其断层的形变活动给出了断层形变活动分析模型[26],以其简便性、有效性使用至今。然而,在分析板内断层错动变形活动时发现其变形活动比板缘断层活动更加复杂,以至于存在模型描述失败或计算结果不佳的情况。为此,有学者对模型进行一定程度的系统构建与完善,使其适合震间与震时、单一断层与多断层系统、整体形变活动与细部形变特征等的描述,同时明确了有关参数的含义[27]。本文利用完善后的模型对鲜水河断裂带及周围区域的形变活动和特点进行了描述和分析。
-
Savage等给出了描述断层活动的数学模型[26]:
$$\mu (x) = \frac{{{V_{\rm b}}}}{\pi }\arctan \left(\frac{x}{D}\right)$$ (1) 其中:
${V_{\rm b}}$ 和$D$ 分别为断层两侧板块的相对运动速率和断层活动的锁定深度。 -
使用中我们发现,原始模型有时不能获得很好的描述,这与地质构造环境、断层的活动状态等密切相关。结合板内形变的复杂情况及应用实践,完善后的模型则为:
$$\mu (x) = {\mu _0} + \frac{{{V_{\rm b}}}}{\pi }\arctan \left(\frac{{x - x{}_{\rm c}}}{D}\right)$$ (2) 式中:
${\mu _0}$ 、${V_{\rm b}}$ 、${x_{\rm c}}$ 和$D$ 分别为研究区的整体运动速率、断层两侧块体的相对运动速率、形变中心位置和广义锁定深度。这些参数及其误差可基于最小二乘法通过反演计算获得。 -
鲜水河断裂带是中国大陆地震活动最强的断裂带之一,位于中国西南地区,起始于甘孜北面的英达,经朱楼、旦都、瓦各及炉霍南侧延伸至道孚,止于康定以南,由多条断层组成;断裂带走向NW40°~50°,总体略呈NE向微凸的弧状,长约300 km;西北段相对比较简单,大致为单线延伸几何结构[28-29](图1)。古生代时期,鲜水河断裂带开始形成,经过缓慢发展,到了第四纪,该断裂带发生了复杂的构造演化,现今成为大型的左旋走滑断裂带[30]。近年来,该地区发生过多次5级以上地震,表明该带及附近区域构造活动较为强烈。
-
为了监测与研究鲜水河断裂带的形变特征,“中国地壳运动观测网络”和“中国大陆构造环境监测网络”等国家重大科学工程在研究区进行了多期GPS重复测量。本文以2008年汶川8.0级地震为节点,分别选取1999—2007年和2009—2013年2个时段的地壳运动场计算结果进行分析。选择研究区域东南的四川盆地作为参考块体,取得其他区域相对于四川盆地的运动场结果,再进行恰当的滤波计算[12-13, 31],其结果见图2。考虑到汶川大地震断裂带破裂的长度与空间分布,选择具有针对性的区域结果进行投影等有关计算[14],具体区域与计算结果如图2~3所示(其中:红色线表示鲜水河断裂带,鲜水河以东的灰线表示龙日坝断裂带,鲜水河以西的灰线疑似形变中心,该形变中心存在与否及位置有待模型计算验证)。关于从场角度的描述有关研究已给出[12],不再赘述,在这主要从上述模型的角度探讨断裂带的活动与动态变化特征等。
从图2a和3a可以看出,汶川大震前后的1999—2007年和2009—2013年鲜水河断裂带整体左旋走滑活动清晰可见,但量值上并不完全等同(后面的分析中将给出具体结果)。这可能跟汶川大震的孕育与发生密切相关。另外,垂直于断裂带的运动性质也具有显著的分段特征,它们大致以龙日坝断层及其向西的几何延伸线(灰线)为界,断裂带北段在2个时期均呈压性,南端均呈张性(图2b,3b)。但从相对于参考基准的角度观察,2个时段的结果是有明显差异的,从运动的大小与方向上都表明了这一点,至少体现了大震前后在较大空间范围内存在着运动的差异。此外,由于西部地区的运动差异使我们怀疑可能有类似断层的变形中心存在,当然这仅是一个初步的假定,我们沿着龙日坝断裂带进一步向西延伸,希望能够进一步去了解它。
通过模型计算,从图4(粗黑线为断裂带的中心位置,粗灰线为形变中心位置,带有误差棒的黑点为观测值,蓝色曲线为模型值)可以看出,汶川大震前(1999—2007年)左旋走滑量约为9.2 mm/a(图4a),震后(2009—2013年)约为6.9 mm/a;此外,形变中心震前偏离了鲜水河断裂带的位置,震后则回归于鲜水河断裂带。这说明大震前后鲜水河断裂带的左旋走滑活动存在着一定的差别。这与参考文献[25]给出的汶川地震前横跨龙门山断裂中段和西南段的GPS站速度剖面(图2)非常一致。参考文献[32]给出的鲜水河不同段滑动速率结果也说明,鲜水河断裂滑动速率从北到南逐渐减小[32],也与本文得到的结果一致。
如前所述,鲜水河断裂带的张压活动具有分段性,北段为压性,南段为张性。因此,利用模型时也按北、南两段进行计算。在计算北段震前的滑动时发现,此时该段不具有断层滑动活动的模式,也就是说失去了断层的活动特征,这显然是一种不正常的活动模式(图5)。由此可知,震后该段恢复了断层属性的活动模式。断裂带北段的压性活动约为2.1 mm/a,但形变中心却偏离了断裂带约40 km。
图6为鲜水河断裂带南段2个时段的结果。汶川大震前(图6a)该段的张性运动约为3.7 mm/a,而震后(图6b)却约为8.3 mm/a,两者之间出现了显著的差别,震后的张性活动比震前大了约4.6 mm/a。这应归属于汶川大震的孕育与发生,否则难以给出合理的解释。
图7给出了鲜水河断裂带以西,平行于龙日坝断裂带的GPS速度剖面观测与模型计算结果。可以看出,汶川地震前后,区域形变并不是块体型的准均匀形变模式,而是具有类似于断裂带形变中心的形变模式。然而,目前我们在龙日坝断裂带延伸线或平行该线方向上并没有发现断裂带的存在。上述形变特征反映的是类似断层大体与鲜水河断裂带平行,并广泛地分布于这一区域,未观测到可辨的断层性走滑活动。由此推测,这种现象可能与地质构造的特殊环境和周边的动力环境密切相关。总体而言,汶川地震前后该区域平行于延伸线的运动分量存在较大的差异(图7),震前的右旋性旋剪形变大于震后。模型计算结果显示,理论上差异达6.5 mm/a左右,这是汶川地震前后包括鲜水河在内的该区域模型化差异活动最突出的变化结果。由此说明,一个大震的孕育直到临震前在较大的空间范围上都可能存在一些形变异常。
-
有分析表明,汶川大震临震之前不仅龙门山断裂带的活动已处于完全闭锁状态,甚至其北部的龙日坝断裂带的活动也近闭锁[27]。这就是说鲜水河断裂带以东地区的形变与构造活动临震前是一种显著的异常形变状态。那么,是否仅局限于此而在更大的空间范围上就不存在异常现象呢?本文的分析结果表明,龙门山和龙日坝断裂带附近的鲜水河断裂带及其以西更大的空间范围上也存在着可识别形变异常现象。这一现象的揭示让我们深入认识到,大震孕育至临震状态时形变异常可能不仅仅在发震断裂带上,其周边地区也可能会伴随产生与其相协调的形变异常现象。因此,异常的识别与判定应是系统性的、多角度的、甚至具有一定空间尺度的。
关于鲜水河断裂带以西的区域存在着疑似断层型的形变活动模式比较特殊,该模式可能与特殊的地质构造环境和周边的动力环境有密切的关系,需要今后能够通过观测与研究进一步认识类似的异常变化。就大震前兆异常的分析而言,跟踪监测与研究它的动态变化对于分析某种构造是否孕育着大震是很有意义的。既可以从较大空间尺度资料中提取针对性的动态形变变化信息,同时也有助于对周边动力环境变化的分析。然而,是什么机理使该区域存在类似断层的形变活动模式却未形成与之相应的地壳破裂等疑问,还有待于深入探索。
The Sliding of Xianshui River Fault Before and After Wenchuan MS8.0 Earthquake and Its Dynamic Changes of Western Crustal Deformation
-
摘要: 基于鲜水河断裂带及邻近区域GPS资料,以发生汶川MS8.0大震的2008年为时间节点,利用视错动变形系统分析模型研究了鲜水河断裂带及邻区不同时段的滑动与区域形变活动特征,结果表明:汶川大震前后鲜水河断裂带构造活动的动态变化是显著的,其西部地区地壳形变的动态变化也较清晰,这些不仅说明了汶川大震的孕育在较大的空间范围存在不同特征的形变与构造活动,同时其西部地区也存在着疑似断裂带形变中心的地壳动态变化,这可能与特定的构造和动力环境密切相关,为我们识别大震前兆提供了更多的可能。Abstract: Based on the GPS repetition measurement data of the Xianshui river fault and adjacent area, taking the 2008 Wenchuan Ms
8.0 earthquake as the time node, the characteristics of the sliding and regional deformation activities of Xianshui river fault and the adjacent area are studied by using the dynamic deformation system analysis model. The following results and views are obtained: The dynamic change of the tectonic activity of the Xianshui river fault before and after the Wenchuan earthquake is remarkable, and the dynamic change of crustal deformation in the west of the fault is also clear. These not only indicate that there are different characteristics of deformation and tectonic activities in the larger space range of Wenchuan earthquake. At the same time, there is similar deformation center in the west of the fault. This may be closely related to the specific structure and dynamic environment, and providing more possibilities for us to identify the precursors of large earthquakes. -
[1] 丁国瑜, 田勤俭, 孔凡臣, 等. 活断层分段: 原则、方法及应用[M]. 北京: 地质出版社, 1993. [2] 江在森, 马宗晋, 张希, 等. GPS初步结果揭示的中国大陆水平应变场与构造变形[J]. 地球物理学报, 2003(3): 352-358. doi: 10.3321/j.issn:0001-5733.2003.03.012 [3] 杨国华, 谢觉民, 韩月萍. 华北主要构造单元及边界带现今水平形变与运动机制[J]. 地球物理学报, 2001(5): 645-653. doi: 10.3321/j.issn:0001-5733.2001.05.008 [4] 王敏, 沈正康, 牛之俊, 等. 现今中国大陆地壳运动与活动块体模型[J]. 中国科学(D辑: 地球科学), 2003(S1): 21-32, 209. [5] Shen Z K, Wang M, Li Y X, et al. Crustal deformation along the Altyn Tagh fault system, western China, from GPS[J]. Journal of Geophysical Research: Solid Earth, 2001, 106(B12): 30607-30621. doi: 10.1029/2001JB000349 [6] Shen Z K, Lü J N, Wang M, et al. Contemporary crustal deformation around the southeast borderland of the Tibetan Plateau[J]. Journal of Geophysical Research: Solid Earth, 2005, 110(B11): B11409. [7] Zhang P Z. A review on active tectonics and deep crustal processes of the Western Sichuan region, eastern margin of the Tibetan Plateau[J]. Tectonophysics, 2013, 584: 7-22. doi: 10.1016/j.tecto.2012.02.021 [8] 杨博, 占伟, 刘志广. 晋冀蒙交界区水平向构造活动的基本特征及动态变化[J]. 地震研究, 2013, 36(4): 472-477, 553. doi: 10.3969/j.issn.1000-0666.2013.04.010 [9] 刘琦. 汶川8.0级和芦山7.0级地震对区域地壳形变场及应力场的影响[D]. 中国地震局地质研究所, 2018. [10] Tong X, Sandwell D T, Smith-Konter B. High-resolution interseismic velocity data along the San Andreas Fault From GPS and InSAR[J]. Journal of Geophysical Research: Solid Earth, 2013, 118(1): 368-389. [11] Lindsey E O, Fialko Y. Geodetic slip rates in the southern San Andreas Fault system: Effects of elastic heterogeneity and fault geometry[J]. Journal of Geophysical Research: Solid Earth, 2013, 118(2): 689-697. doi: 10.1029/2012JB009358 [12] Kaneko Y, Fialko Y, Sandwell D T, et al. Interseismic deformation and creep along the central section of the North Anatolian Fault (Turkey): InSAR observations and implications for rate-and-state friction properties[J]. Journal of Geophysical Research: Solid Earth, 2013, 118(1): 316-331. doi: 10.1029/2012JB009661 [13] Savage J C, Simpson R W. Clustering of GPS velocities in the Mojave Block, southeastern California[J]. Journal of Geophysical Research: Solid Earth, 2013, 118(4): 1747-1759. doi: 10.1029/2012JB009699 [14] Garthwaite M C, Wang H, Wright T J. Broadscale interseismic deformation and fault slip rates in the central Tibetan Plateau observed using InSAR[J]. Journal of Geophysical Research: Solid Earth, 2013, 118(9): 5071-5083. doi: 10.1002/jgrb.50348 [15] Müller M D, Geiger A, Kahle H G, et al. Velocity and deformation fields in the North Aegean domain, Greece, and implications for fault kinematics, derived from GPS data 1993-2009[J]. Tectonophysics, 2013, 597-598: 34-49. doi: 10.1016/j.tecto.2012.08.003 [16] Wang Y Z, Wang M, Shen Z K, et al. Inter-seismic deformation field of the Ganzi-Yushu fault before the 2010 Mw6.9 Yushu earthquake[J]. Tectonophysics, 2013, 584: 138-143. doi: 10.1016/j.tecto.2012.03.026 [17] 马杏垣. 中国岩石圈动力学图集[M]. 北京: 地质出版社, 1989: 16-58. [18] 刘静, 张智慧, 文力, 等. 汶川8级大地震同震破裂的特殊性及构造意义——多条平行断裂同时活动的反序型逆冲地震事件[J]. 地质学报, 2008, 82(12): 1707-1722. doi: 10.3321/j.issn:0001-5717.2008.12.007 [19] 王敏. 基于GPS同震位移场约束反演2008年5.12汶川大地震破裂空间分布[J]. 地球物理学报, 2009, 52(10): 2519-2526. doi: 10.3969/j.issn.0001-5733.2009.10.010 [20] 刘杰, 易桂喜, 张致伟, 等. 2013年4月20日四川芦山M7.0级地震介绍[J]. 地球物理学报, 2013, 56(4): 1404-1407. [21] 赵静, 江在森, 武艳强, 等. 汶川地震前龙门山断裂带闭锁程度和滑动亏损分布研究[J]. 地球物理学报, 2012, 55(9): 2963-2972. [22] 鲁小飞, 谭凯, 赵斌, 等. 汶川地震前后红河断裂闭锁程度和滑动亏损研究[J]. 大地测量与地球动力学, 2019, 39(5): 464-468, 475. [23] 占伟, 武艳强, 梁洪宝, 等. GPS观测结果反映的尼泊尔Mw7.8地震孕震特征[J]. 地球物理学报, 2015, 58(5): 1818-1826. doi: 10.6038/cjg20150532 [24] 赵静, 刘杰, 任金卫, 等. 汶川8.0级地震前GPS与跨断层资料反映的运动与变形演化特征[J]. 中国地震, 2018, 34(3): 414-427. doi: 10.3969/j.issn.1001-4683.2018.03.003 [25] 赵静, 刘杰, 任金卫, 等. 汶川、芦山地震前龙门山断裂带地壳形变特征对比分析[J]. 地球学报, 2019, 40(1): 186-198. doi: 10.3975/cagsb.2018.110201 [26] Savage J C, Burford R O. Geodetic determination of relative plate motion in central California[J]. Journal of Geophysical Research, 1973, 78(5): 832-845. doi: 10.1029/JB078i005p00832 [27] 杨博, 万文妮. 断层及邻区视错动变形系统分析模型[J]. 大地测量与地球动力学, 2018, 38(8): 796-800. [28] 康四林. 鲜水河断裂GPS监测研究[J]. 科技创新与应用, 2014(26): 62-63. [29] 杨珍. 基于InSAR数据的川西地区断层活动研究[D]. 西安: 长安大学, 2017. [30] 郭长宝, 杜宇本, 张永双, 等. 川西鲜水河断裂带地质灾害发育特征与典型滑坡形成机理[J]. 地质通报, 2015, 34(1): 121-134. doi: 10.3969/j.issn.1671-2552.2015.01.010 [31] 杨国华, 杨博, 陈欣, 等. 华北现今三维形变场空间变化的基本特征[J]. 武汉大学学报(信息科学版), 2013, 38(1): 31-35. [32] 李铁明, 祝意青, 杨永林, 等. 综合利用多种地壳形变观测资料计算鲜水河断裂带现今滑动速率[J]. 地球物理学报, 2019, 62(4): 1323-1335. doi: 10.6038/cjg2019L0723 -