2. 中国地质科学院岩溶地质研究所, 桂林 541004
2. Institute of Karst Geology, Chinese Academy of Geological Sciences, Guilin 541004, China
印度板块与欧亚板块的碰撞不仅造成了青藏高原的隆升与地壳缩短,还驱动着高原物质的东向挤出,使得青藏高原内部一系列大型走滑断裂开始启动(Molnar and Tapponnier, 1975; Tapponnier and Molnar, 1977; Harrison et al., 1992; Houseman and England, 1993).关于青藏高原的变形模式一直是一个争论的热点,现在较流行的观点是刚性微板块模型(Peltzer and Tapponnier, 1988; Thatcher, 2007)和连续变形模型(England and McKenzie, 1982; England and Houseman, 1986),这两个模型的区别在于高原内部大型走滑断裂滑动速率的大小.刚性微板块模型需要块体周边的断裂具有较高的走滑运动速率,而连续变形模型所对应的断裂滑动速率较小.甘孜—玉树断裂位于青藏高原中东部,东南与鲜水河断裂左行左阶斜列形成甘孜拉分盆地,西北与风火山断裂斜接;三者共同协调羌塘地体和巴颜喀拉地体之间的相对运动(图 1).
根据甘孜—玉树断裂带的几何形态和历史地震记录将断裂从东往西依次划分为甘孜段、玛尼干戈段、邓柯段、玉树段和当江段(周荣军等, 1996).前人基于地质地貌调查对甘孜—玉树断裂带的活动速率做过很多研究,李闽峰等(1995)根据阶地和洪积扇位错量及其位错起始年代获得该断裂带相古—竹庆段左旋滑动速率为3.3~5.0 mm·a-1,当江段为7.3 mm·a-1.周荣军等(1996)根据甘孜—当江段不同地貌体的位错及其相应的位错起始时间的限定获得该断裂带甘孜段、玛尼干戈段、邓柯段和当江段晚第四纪以来的左旋滑动速率分别为:3.4±0.3 mm·a-1、7.0±0.7 mm·a-1、7.2±1.2 mm·a-1和7.3±0.6 mm·a-1.闻学泽等(2003)根据不同的地貌位错量及其相关沉积物年龄获得该断裂带甘孜—相古段近50 ka以来平均左旋滑动速率为12±2 mm·a-1.徐锡伟等(2003a)根据阶地位错量及其位错起始时间获得该断裂带甘孜段和玛尼干戈段全新世左旋滑动速率为14±3 mm·a-1.Lin等(2011a, 2011b)根据2010年玉树地震同震地表位移和全新世古地震复发周期获得该甘孜—玉树断裂带玉树段平均左旋滑动速率为2~5 mm·a-1.石峰等(2013)沿断裂带甘孜—玛尼干戈段进行了断错地貌分析和晚第四纪滑动速率计算,获得该断裂段左旋走滑速率介于7.1~11.6 mm·a-1之间.综上所述,甘孜—玉树断裂带晚第四纪活动性总体研究程度较高,但是主要集中在断裂东南段,对于断裂西北当江段晚第四纪活动性及其几何结构的研究非常薄弱,不利于对甘孜—玉树断裂带活动性进行整体评价与分析.
当江断裂位于甘孜—玉树断裂带的最西段,据历史地震资料记载该断裂1738年发生过6.5级地震(国家地震局震害防御司, 1995),后来经过调查核实,修正为7.5级(周荣军等, 1997; 袁道阳等, 2011; 李彦宝等, 2016; Wu et al., 2017).本文基于遥感影像解译和数字高程地形模型(DEM)数据分析,通过对当江断裂进行详细的野外构造地貌填图,掌握了断裂的几何展布和运动性质;通过对研究区断错地貌的测量并结合相应的年代测试估算了断裂晚第四纪以来的滑动速率.研究成果可以为治多县地区的地震危险性评价提供基础资料,也有利于对甘孜—玉树断裂带的整体运动学特征进行探讨.
2 区域构造背景受印度板块和欧亚板块碰撞的远程效应影响,中-晚中新世以来在青藏高原东北部形成了一系列大型的走滑与逆冲断裂(来庆洲等, 2006; 杨农和张岳桥, 2010; 王伟涛等, 2014; 李安等, 2016)(图 1a).北缘的海原—六盘山断裂带全长约240 km,由西向东整体表现出由走滑兼挤压到走滑-挤压的构造转换样式(李强等, 2013).根据GPS水平速度场资料,海原断裂左旋滑动速率为5 mm·a-1,六盘山北段挤压缩短速率为2.4 mm·a-1(李强等, 2013).东昆仑断裂带全长约1600 km,地质调查结果与GPS数据均表明东昆仑断裂带晚第四纪以来的左旋滑动速率为10 mm·a-1(Zhang et al., 2004; Fu and Awata, 2007).东缘的龙门山断裂带作为巴颜喀拉地体和华南地体的边界断裂,经历了三叠纪的印支造山运动和新生代喜马拉雅运动,形成了复杂的龙门山褶皱逆冲断裂带和四川前陆盆地,以及在青藏高原东缘和成都平原这几十千米范围内造就了约4000 m的高差(Chen and Wilson, 1996).晚第四纪以来,龙门山逆冲断裂带主要活动的区域由西往东集中在汶川—茂县断裂段、北川—映秀断裂段和江油—灌县断裂段以及山前隐伏断裂段,呈现出由西向东传播变形的特征(徐锡伟等, 2008; 任俊杰等, 2012).龙门山断裂带以逆冲作用为主兼有右旋走滑分量,GPS数据显示晚第四纪以来的垂直运动速率为1~3 mm·a-1(Zhang et al., 2004; Zhang and Wei, 2011).在龙门山区域内的一些次级断裂,主要包括龙日坝断裂和南北向的岷山断裂,龙日坝断裂晚第四纪以来的右旋滑动速率约5.5 mm·a-1(Zhang and Wei, 2011);岷山断裂主要由虎牙断裂、岷山隆起、岷江断裂组成,前人研究结果表明虎牙断裂和岷江断裂都是左旋逆冲断裂;其中,虎牙断裂的走滑速率是1.4 mm·a-1,倾滑速率是0.3~0.5 mm·a-1;而岷江断裂的走滑速率约1 mm·a-1,倾滑速率是0.37~0.53 mm·a-1(Jones et al., 1984; 赵小麒等, 1994; 唐文清等, 2004; 周荣军等, 2006).分析青藏高原东南缘GPS速度场资料结果表明,青藏高原物质在向东挤出的过程中遇到稳定且刚性的华南块体的阻挡后,转向南东甚至向南运动.这种推动力同时带动川滇地体围绕喜马拉雅东构造结做顺时针转动,使得川滇地块东侧断裂以左旋走滑为主,西侧地块则以右旋走滑为主.其中甘孜—玉树、鲜水河、安宁河、则木河、大凉山和小江断裂构成了青藏高原东南部东向挤出的东边界,其左旋滑动速率分别为7.3~14.7 mm·a-1、8.9~17.1 mm·a-1、5.1±2.5 mm·a-1、2.8±2.3 mm·a-1、7.1±2.1 mm·a-1和9.4±1.2 mm·a-1(王阎昭等, 2008).
3 断错地貌特征当江断裂西起于多彩群地以西,往南东经当江荣村、当江、立新乡、哲达村、达约列一直到达彭董以东,走向NW-NWW,延伸约150km(图 2).此外,当江断裂东南段与玉树断裂右阶斜列,形成一挤压阶区,阶区内发育隆宝滩和哈秀次级断裂.在当江断裂中部,治多盆地以南发育一次级断裂(这里称为治多断裂),该次级断裂走向NW,延伸约30 km,从影像上测量该次级断裂与当江断裂之间的锐夹角为15°~20°.根据断裂沿线断错地貌特征表明,治多断裂是一条以左旋走滑为主兼有正断分量的晚第四纪活动断裂.
当江断裂北侧发育通天河一级水系,受通天河各级支流的侵蚀改造,研究区发育多级河流阶地.通过遥感影像解译,在当江荣村,通天河一级支流聂恰曲左拐约6 km;在当江,通天河一级支流登额曲左拐约16 km,反映了断裂的长期左旋走滑活动习性.当江断裂周围被一系列海拔约5000 m的山峰环绕,这些山峰常年冰雪覆盖,受季节影响在研究区发育多级晚第四纪地貌面.当江断裂沿线断错地貌发育,主要包括断层三角面、断层陡坎、挤压脊、闸门脊、马鞍状地貌、垭口以及系统左旋断错冲沟水系、阶地、洪积扇和山脊;断裂沿线还发育当江荣、当江和哲达等一系列串珠状第四纪断层谷地,表明该断裂是一条晚第四纪或全新世强烈活动的左旋走滑断裂.以下就几个典型的断错地貌点进行描述.
3.1 立新乡观察点(33°35′ 55.43″,96°03′53.15″)基于Google earth影像判读,在立新乡断裂线性特征非常清晰,并且横切自北往南流的大型洪积扇(图 3).野外实地调查发现该线性特征为一断续延伸的断层陡坎,断坎高度约0.5~1 m.被断裂横切的大型洪积扇其扇面上断头沟保存较好,通过Google earth影像可以清晰的识别出3条断头沟和一条完整的冲沟, 冲沟已被左旋断错,位错量为38±5 m;断头沟最大位错为330±20 m,代表该洪积扇的位错量.在该洪积扇最东侧登额曲河流两侧发育两级不对称河流阶地,并且都是堆积阶地,河流西岸阶地保存较好.通过影像解译和野外实地调查发现河流西岸T1、T2阶地均被左旋断错,野外用皮尺测量T2/T1阶地坎左旋位错29±1 m(图 4).
冲沟位移变形是走滑断裂运动的重要指标.遥感影像上显示断裂沿线断错冲沟非常普遍,其中当江谷地南冲洪积扇上冲沟左旋发育最好(图 5).野外调查结合Google earth影像解译对当江谷地南冲洪积扇上的典型冲沟位错量进行了实测(图 6),结果表明断裂错断规模不等的冲沟,位错量从几米至上百米不等.年轻冲沟左旋位移为3.6 m(图 7a)以及位移变形量为5m和9.8 m的小规模冲沟同步左旋位错(图 7b).规模较大的冲沟位错量有65 m(图 7c)、87m(图 7d),最大约204 m.在断裂经过处断塞塘呈串珠状展布,并且断裂沿线分布地震鼓包和断头沟等断错地貌(图 5和图 6).在当江南冲洪积扇中部,一小型年轻洪积扇被断裂左旋断错,形成洪积台地(图 6).洪积台地东西边界和中间均发育冲沟,并且台地遭受冲沟侵蚀切割,冲沟与台地同步左旋位错,洪积台地东侧冲沟位错量为65 m,代表台地东边界位错量,考虑到冲沟的侵蚀破坏作用,因此,洪积台地东边界位错量为65±5 m(图 7c).同样,台地中间的位错量为87±5 m(图 7d).
结合Google earth影像和野外调查发现在当江荣村东侧冲洪积扇扇面上,一系列自北往南流的小冲沟发生左旋位错,位错量有5 m(图 8a),3.7 m以及同一冲沟但已被废弃的老冲沟位错量为11 m(图 8b).在当江荣村东侧发育一洪积扇体,该扇体左旋位错在Google earth影像上非常清晰并且扇体规模与立新乡观察点处的洪积扇相当(图 8c和图 9a).该洪积扇扇体西边界受河流侵蚀改造,已无法恢复;扇体东边界经平移恢复,效果较好(图 9b中扇体东边界的白线),经测量左旋位错量约350 m.考虑到该洪积扇长期受自然侵蚀的影响,断裂两侧扇体的东边界也受到一定的侵蚀破坏,在这分配±20 m的误差范围来反映这种不确定性;因此,该洪积扇左旋位错量为350±20 m.
通天河一级支流聂恰曲河流从达生村经过,并在该处东西岸发育两级不对称河流阶地(图 10),其中T0和T1为堆积阶地,西岸受河流凹岸侵蚀破坏,只有局部保留,但仍可识别(图 11a);东岸阶地保留较好,发育宽平阶地面(图 11b).T2在河流东西两岸也有保存,东岸为堆积阶地,主要为宽平阶地面;西岸为基座阶地.Google earth影像解译和野外地质调查结果表明,河流两岸的阶地均被左旋断错,经皮尺测量T2/T1阶地坎左旋位错量为29±1 m(图 10和图 11).
本研究主要采用光释光(OSL)方法测年,由中国地震局地壳应力研究所光释光实验室专业技术人员完成.其中光释光样品主要采自阶地二元结构内粉细砂层或砾石层之上的细砂加积层.
本文测定的地貌面为T1、T2阶地和洪积扇,具体描述如下.
T1阶地样品采自当江荣村东侧东西向河流的北岸(图 9a),样品采自粉细砂层,取样深度为0.7 m(DJ-OSL-8)(图 12a,表 1),测得年龄为距今3.04±0.37 ka,为全新世晚期.
T2阶地样品采自当江荣村东侧东西向河流的北岸(图 9a)和当江荣村西侧聂恰曲河流的北岸.在当江荣村东侧东西向河流北岸,T2阶地面以下1.5 m处粉细砂层内取样(DJ-OSL-9)(图 12b,表 1),测得年龄为距今8.76±0.99 ka,为全新世早期.在当江荣村西侧聂恰曲河流北岸,T2阶地面之下约0.5 m处取了2个光释光(OSL)年龄样品(图 12c,表 1),采样位置是在砾石层之上的细砂加积层,其年代最接近T2阶地面的形成时代,获得的年代为距今10.39±1.21 ka(DJ-OSL-10) 和10.61±1.22 ka(DJ-OSL-11),为全新世早期.
洪积扇样品采自当江南通天河一级支流登额曲河流南岸冲洪积扇面前缘之下约1 m处(图 5a和图 12d,表 1),采样位置是在砾石层之上的细砂加积层,获得的年代为距今9.21±1.83 ka(DJ-OSL-1) 和9.88±1.88 ka(DJ-OSL-2),为全新世早期.由于该冲洪积扇前缘与河流T2阶地面共面,因此该测年结果与T2阶地面的形成时间非常接近(图 5).
5 晚第四纪滑动速率理论上,走滑断裂滑动速率可以用位移除以位移的起始时间获得(邓起东等, 2004).但要准确获得断裂的滑动速率并非易事,不同的研究者测定的同一条断裂的滑动速率可以相差3倍,主要原因在于很难确定位移的起始时间(张培震等, 2008).一种模式认为河流的漫滩在被废弃并形成新的阶地之前,一直对其上一级阶地陡坎进行侵蚀作用;因此,阶地陡坎的位移在河漫滩被废弃并形成新的阶地之前是不会被完整的保留下来的(Tapponnier et al., 2001).另一种模式认为当河漫滩废弃形成阶地后,断裂的走滑位移致使河流两侧的阶地陡坎都遭受错动,其中一侧的下游阶地陡坎被错入河道而遭受侵蚀;而另一侧的下游阶地陡坎被错离河道,受到上游河岸的保护而免遭侵蚀;因此,错离河道阶地陡坎的位移累积的起始时间就是该阶地陡坎上阶地面的形成时间(张培震等, 2008)(图 13a, b).在这种情况下,上阶地面的废弃时间不仅是该阶地陡坎位错起始时间的上限,同时也最接近该阶地陡坎位错的起始时间(张培震等, 2008);而错入河道的阶地陡坎一直遭受河流的侵蚀,假设河流的侵蚀作用使水平位移完全没有保留,那么只有当下一级阶地形成时,水平位移才开始累积,该阶地坎水平位移应为下阶地面形成以来累积的.在这种情况下,下阶地面的废弃时间应该最接近该阶地陡坎位错的起始时间(图 13c, d).此外,在阶地陡坎位移恢复的过程中,河流的宽度以及阶地面废弃时间的穿时性也具有非常重要的作用(Cowgill, 2007).较理想的情况是,阶地陡坎上阶地面的形成时间与下阶地面的形成时间间隔很小,这样就能很好的限定该阶地陡坎位错的起始时间.应用这种方法已经很好地获得了阿尔金断裂带中段的晚第四纪滑动速率为9~14 mm·a-1(Cowgill et al., 2009).
尽管前人已经建立了较完整的河流阶地演化模式,但是也要根据实际问题具体问题具体分析.根据上文描述,当江断裂沿线阶地位错保存较好的地方主要位于立新乡东侧登额曲河流西岸和达生村聂恰曲河流东岸,达生村聂恰曲河流西岸断裂南侧T1阶地被侵蚀殆尽没有保留,无法与断裂北侧的T2/T1阶地坎对应(图 10a, b和图 11a),只是作为一种判断标志.由于立新乡观察点和达生村观察点处的河道较宽,在此采用上下阶地面的形成时间共同限定阶地陡坎位移的起始时间.立新乡和达生村T2/T1阶地坎水平位移量均为29±1 m,研究区T1的形成时间为3.04±0.37 ka,T2的形成时间为8~10 ka,水平滑动速率约为7±3 mm·a-1.当江南中部冲洪积扇上洪积台地左旋位错量为65~87 m,年代为距今9.5±1.8 ka,据此估算的全新世早期以来的水平滑动速率为8±1 mm·a-1.
在研究区立新乡和当江荣村东侧观察点处,两个规模相当的大型洪积扇被断裂左旋断错,位错量分别为330±20 m和350±20 m,平均约340±20 m.目前并没有获得该大型洪积扇的年代学样品,但是考虑到高原内部区域气候变化的一致性,可以利用气候冷暖变化的时间对区域上冲洪积扇的形成时间进行估计.研究区海拔约5000 m,冰期内水不会自然流动,更不会发育冲洪积扇堆积,因此,研究区内冲洪积扇的形成主要是在暖期.已有的研究结果表明,青藏高原东部河流阶地与洪积扇的形成时间应该是在末次冰期间冰阶和全新世大暖期,这两个时期分别对应于深海氧同位素的MIS-3(60~30 ka BP)阶段和MIS-1(10 ka BP至现在)阶段(Thompson et al., 1997).姚檀栋等(1997)在昆仑山古里雅冰芯中得到的末次间冰期时间为58~32 ka BP,与深海氧同位素记录非常接近.根据气候变化所获得的年龄与区域内T4~T5阶地的形成年龄也非常接近(表 2).因此,在这里应用末次间冰期时间段(58~32 ka BP)作为研究区立新乡和当江荣村西侧大型洪积扇的形成时间的上下限年龄,计算得到当江断裂晚更新世中期以来的水平滑动速率约为8±2 mm·a-1.
综上所述,当江断裂全新世以来的滑动速率为7±3 mm·a-1,晚更新世中期以来的滑动速率为8±2 mm·a-1;计算得到的当江断裂晚更新世中期以来的水平滑动速率与断裂全新世以来的水平滑动速率值在误差范围内一致,说明当江断裂晚更新世以来具有长期的一致活动性,并且水平滑动速率为7±3 mm·a-1.该结果与李闽峰等(1995)和周荣军等(1996)根据当江南洪积扇位错量及其扇体沉积物年代获得的当江断裂约7.3 mm·a-1和7.3±0.6 mm·a-1的滑动速率值一致.由于青藏高原各地块的不均匀挤出,南侧羌塘地块较北侧巴颜喀拉地块逃逸快,使得其边界断裂甘孜—玉树断裂带晚第四纪以来表现为左旋走滑.基于大地测量学(GPS和InSAR)方法,最新GPS数据结果表明羌塘地体向东运动的速度为22~28 mm·a-1,而巴颜喀拉地体向东运动的速度是15~19 mm·a-1,羌塘地体相对于巴颜喀拉地体向东的运动速度为7~9 mm·a-1(Shen et al., 2005, 2009).根据跨甘孜—玉树断裂带的GPS数据结果表明该断裂现今左旋滑动速率为6.6±1.5 mm·a-1(王阎昭等, 2010).基于InSAR数据对跨玉树—当江段2003—2010年的地形变场研究结果表明该断裂段震间平均左旋走滑速率为6.4 mm·a-1(Liu et al., 2011).因此,本文所获得的当江断裂晚第四纪以来的左旋滑动速率也与大地测量结果较接近.
6 讨论分析 6.1 甘孜—玉树断裂带走滑运动量的分解与吸收各研究者根据不同地貌体的位错及其相应的位错起始时间的限定所报道的甘孜—玉树断裂晚第四纪以来的滑动速率分别为:甘孜段,3.4±0.3 mm·a-1(周荣军等, 1996),11.5±2.4 mm·a-1(闻学泽等, 2003),14.3±3 mm·a-1(徐锡伟等, 2003a);玛尼干戈段,3.3~5.0 mm·a-1(李闽峰等, 1995),7.0±0.7 mm·a-1(周荣军等, 1996),12.8±1.7 mm·a-1(闻学泽等, 2003),13.4±2 mm·a-1(徐锡伟等, 2003a);邓柯段,7.2±1.2 mm·a-1(周荣军等, 1996),11.3±1.8 mm·a-1(闻学泽等, 2003),玉树段,2~5 mm·a-1(Lin et al., 2011b);当江段,7.3 mm·a-1(李闽峰等, 1995),7.3±0.6 mm·a-1(周荣军等, 1996),7±3 mm·a-1(本文)(表 3).综合以上分析,发现甘孜—玉树断裂带在邓柯以东断裂平均滑动速率总体为12~14 mm·a-1,与鲜水河断裂的平均左旋滑动速率(10~15 mm·a-1)(Wen et al., 1988; Allen et al., 1991; Zhao et al., 1992; 徐锡伟等, 2003a)很接近;而在邓柯以西断裂平均滑动速率总体为2~7 mm·a-1,两者相差较大.对于甘孜—玉树断裂带西北段滑动速率递减的原因,本章给出以下两种解释:(1) 与不同研究者所选择的断错地貌标志以及采用的位错起始年代不同有关,对于邓柯以东12~14 mm·a-1的滑动速率值,仔细分析前人研究成果发现,他们应用的是下阶地面的废弃年代作为阶地坎位错的起始时间,所估算的滑动速率值相对偏大(张培震等, 2008).(2) 是由于在邓柯以西的玉树断裂和当江断裂段上发育很多分支断裂以及断裂之间形成大的阶区.分支断裂对主断裂的运动具有分解作用,这在走滑断裂上非常普遍,典型的有东昆仑断裂和阿尔金断裂.根据前人研究结果表明,东昆仑断裂从西段库赛湖段到东段的塔藏断裂,晚更新世以来的滑动速率从约11 mm·a-1逐渐的下降到≤1 mm·a-1(Van Der Woerd et al., 2000; Li et al., 2005; Kirby et al., 2007; Harkins and Kirby, 2008),滑动速率的下降主要是由于断裂东段尾端的构造转换,以及东昆仑断裂带东段帚状散开的几何形态所致(李陈侠等, 2009).对于阿尔金断裂,断裂从西往东滑动速率不断递减,其左旋滑动速率的下降主要发生在三联点上,走滑速率的减少量转化成了断裂南盘NW向逆冲垂直运动分量(徐锡伟等, 2003b).因此,走滑断裂带的帚状散开的几何形态和构造转换可以造成走滑断裂滑动速率的沿断裂走向上的递减.如北侧的东昆仑断裂,其左旋走滑运动在断裂东南部被近SN向的岷山断裂、虎牙断裂、龙日坝断裂、临江断裂和龙门山断裂所吸收(徐锡伟等, 2003a; 唐文清等, 2004; Zhang et al., 2004; 周荣军等, 2006; 杜义等, 2009),导致龙门山地区发生强震活动.
在前人研究的基础上,结合本次调查研究成果,玉树断裂和当江断裂分别发育巴塘乡次级断裂和治多次级断裂(图 2和图 14).Huang等(2015a)研究结果表明,玉树巴塘断裂在甘孜—玉树断裂带玉树段起着变形分解的作用,分解了玉树段内约1/3的走滑变形,致使玉树段相比于其他断裂段滑动速率偏小.本次野外调查发现在当江断裂中部和治多盆地以南发育治多次级断裂,推测对当江断裂晚第四纪以来的左旋滑动速率会有部分分解作用.此外,当江断裂与玉树断裂呈左旋右阶排列,形成挤压阶区(图 14).推测对两侧断裂滑动速率会有部分吸收,这有待于进一步研究.
对于板内走滑断裂终端样式可以归纳为四种:走滑、旋转、伸展和汇聚(Storti et al., 2003).对于走滑样式,主断裂上的位移主要沿内部次一级断裂分布;对于旋转样式,主要是通过反向断层来调节主断裂上的变形;而伸展和压缩样式则表明,在主断裂一侧发育正断层和挤压结构.本文研究结果表明,当江断裂作为甘孜—玉树断裂西北段终端,明显以走滑样式为主;而与甘孜—玉树断裂带呈左行左阶排列的鲜水河—则木河—小江断裂带,则是以旋转样式为主(吴中海等, 2015).或许甘孜拉分盆地以西更接近于连续变形模型,而甘孜拉分盆地以东更接近刚性微板块模型.事实上,就上地壳尺度来说,刚性微板块模型和连续变形模型之间并没有太大的差别,两者只是一种渐变过渡.随着微板块的继续分割,板块规模越来越小,断裂的数量越来越多,也就越来越接近连续变形模型(Thatcher, 2007).
7 结论(1) 当江断裂在遥感影像上线性特征非常清晰,断裂沿线系统断错冲沟、河流、阶地和洪积扇,是一条强烈活动的全新世左旋走滑断裂,最新活动时代为全新世晚期,距今约3.04 ka.
(2) 根据河流阶地演化与阶地坎位错量和阶地年代分析,以及断错洪积台地位错量和洪积台地面年代数据,估算了当江断裂全新世早期以来的左旋滑动速率为7±3 mm·a-1;根据断错洪积扇位错量和气候地层年代资料,估算了当江断裂晚更新世中期以来的左旋滑动速率为8±2 mm·a-1;两者在误差范围内一致.
(3) 由于青藏高原各地块的不均匀挤出,南侧羌塘地块较北侧巴颜喀拉地块逃逸快,其边界断裂甘孜—玉树断裂带晚第四纪以来表现为左旋走滑,滑动速率自南东向北西呈逐渐减小的趋势,这一方面与断裂带西北段帚状散开的几何结构相对应,另一方面与断裂带各断裂段之间所形成的阶区有关,阶区能够调节并吸收走滑运动量.
致谢本文是依托中央级公益性科研院所基本科研业务专项(DZJ2016-18) 和国家自然科学基金(41602222) 完成的,在论文成稿过程中,匿名审稿人提出了宝贵的修改意见,在此致以诚挚的谢意!
Allen C R, Luo Z L, Qian H, et al.
1991. Field study of a highly active fault zone:The Xianshuihe fault of southwestern China. Geological Society of America Bulletin, 103(9): 1178-1199.
DOI:10.1130/0016-7606(1991)103<1178:FSOAHA>2.3.CO;2 |
|
Chen S F, Wilson C J L.
1996. Emplacement of the Longmen Shan Thrust-Nappe Belt along the eastern margin of the Tibetan Plateau. Journal of Structural Geology, 18(4): 413-430.
|
|
Cowgill E.
2007. Impact of riser reconstructions on estimation of secular variation in rates of strike-slip faulting:revisiting the Cherchen River site along the Altyn Tagh Fault, NW China. Earth Planet. Sci. Lett., 254(3-4): 239-255.
DOI:10.1016/j.epsl.2006.09.015 |
|
Cowgill E, Gold R D, Chen X H, et al.
2009. Low Quaternary slip rate reconciles geodetic and geologic rates along the Altyn Tagh fault, northwestern Tibet. Geology, 37(7): 647-650.
DOI:10.1130/G25623A.1 |
|
Deng Q D, Chen L C, Ran Y K.
2004. Quantitative studies and applications of active tectonics. Earth Science Frontiers (in Chinese), 11(4): 383-392.
|
|
Department of Earthquake Disaster Prevention, State Seismological Bureau. 1995.
Historical Strong Earthquakes Catalogue in China (in Chinese). Beijing: Seismological Press: 277-278.
|
|
Du Y, Xie F R, Zhang X L, et al.
2009. The mechanics of fault slip of MS8.0 Wenchuan earthquake. Chinese J. Geophys. (in Chinese), 52(2): 464-473.
|
|
England P, McKenzie D.
1982. A thin viscous sheet model for continental deformation. Geophys. J. Int., 70(2): 295-321.
DOI:10.1111/j.1365-246X.1982.tb04969.x |
|
England P, Houseman G.
1986. Finite strain calculations of continental deformation 2. Comparison with the India-Asia collision zone. J. Geophys. Res., 91(B3): 3664-3676.
|
|
Fu B H, Awata Y.
2007. Displacement and timing of left-lateral faulting in the Kunlun Fault Zone, northern Tibet, inferred from geologic and geomorphic features. Journal of Asian Earth Sciences, 29(2-3): 253-265.
DOI:10.1016/j.jseaes.2006.03.004 |
|
Harkins N, Kirby E.
2008. Fluvial terrace riser degradation and determination of slip rates on strike-slip faults:An example from the Kunlun Fault, China. Geophysical Research Letters, 35(5): L05406.
|
|
Harrison T M, Copeland P, Kidd W S F, et al.
1992. Raising Tibet. Science, 255(5052): 1663-1670.
DOI:10.1126/science.255.5052.1663 |
|
Houseman G, England P.
1993. Crustal thickening versus lateral expulsion in the Indian-Asian continental collision. J. Geophys. Res.:Solid Earth, 98(B7): 12233-12249.
DOI:10.1029/93JB00443 |
|
Huang X M, Du Y, He Z T, et al.
2015a. Late Quaternary slip rate of the Batang Fault and its strain partitioning role in Yushu area, central Tibet. Tectonophysics, 653: 52-67.
DOI:10.1016/j.tecto.2015.03.026 |
|
Huang X M, Du Y, He Z T, et al.
2015b. Late Pleistocene-Holocene paleoseismology of the Batang fault (central Tibet plateau, China). Geomorphology, 239: 127-141.
DOI:10.1016/j.geomorph.2015.03.026 |
|
Jones L M, Han W B, Hauksson E, et al.
1984. Focal mechanisms and aftershock locations of the Songpan earthquakes of August 1976 in Sichuan, China. J. Geophys. Res.:Solid Earth (1978-2012), 89(B9): 7697-7707.
DOI:10.1029/JB089iB09p07697 |
|
Kirby E, Harkins N, Wang E Q, et al.
2007. Slip rate gradients along the eastern Kunlun fault. Tectonics, 26(2): TC2010.
|
|
Lai Q Z, Ding L, Wang H W, et al.
2006. Thermal history check and supervision of apatite fission-track in the expand process of eastern margin of Qinghai-Xizang Plateau. Sci. China Earth Sci. (in Chinese), 36(9): 785-796.
|
|
Li A, Wang X X, Zhang S M, et al. 2016. The Structural deformation of Baiyanghe thrust fault-fold belt in Jiuxi basin.//Bulletin of the Institute of Crustal Dynamics 28. Beijing:Seismological Press, 1-11.
|
|
Li C X, Xu X W, Wen X Z, et al.
2009. The segmental characteristics of geometrical structure of the East Kunlun Active Fault (Maqin-Maqu segment). Seismology & Geology (in Chinese), 31(3): 441-458.
|
|
Li H B, Van Der Woerd J, Tapponnier P, et al.
2005. Slip rate on the Kunlun Fault at Hongshui Gou, and recurrence time of great events comparable to the 14/11/2001, MW~7.9 Kokoxili earthquake. Earth and Planetary Science Letters, 237(1-2): 285-299.
DOI:10.1016/j.epsl.2005.05.041 |
|
Li M F, Xing C Q, Cai C X, et al.
1995. Research on activity of yushu fault. Seismology & Geology (in Chinese), 17(3): 218-224.
|
|
Li Q, Jiang Z S, Wu Y Q, et al.
2013. Present-day tectonic deformation characteristics of Haiyuan-Liupanshan fault zone. Journal of Geodesy & Geodynamics (in Chinese), 33(2): 18-22.
|
|
Li Y B, Wang H, Liu H G, et al.
2016. The surface rupture remains along the Dangjiang segment of Garzê-Yushu fault and the 1738 AD earthquake. Technology for Earthquake Disaster Prevention (in Chinese), 11(2): 207-217.
|
|
Lin A M, Rao G, Jia D, et al.
2011a. Co-seismic strike-slip surface rupture and displacement produced by the 2010 MW6.9 Yushu earthquake, China, and implications for Tibetan tectonics. J. Geodyn., 52(3-4): 249-259.
DOI:10.1016/j.jog.2011.01.001 |
|
Lin A M, Jia D, Rao G, et al.
2011b. Recurrent morphogenic earthquakes in the past millennium along the strike-slip Yushu Fault, Central Tibetan Plateau. Bulletin of the Seismological Society of America, 101(6): 2755-2764.
DOI:10.1785/0120100274 |
|
Liu Y, Xu C J, Li Z H, et al.
2011. Interseismic slip rate of the Garze-Yushu fault belt in the Tibetan Plateau from C-band In SAR observations between 2003 and 2010. Advances in Space Research, 48(12): 2005-2015.
DOI:10.1016/j.asr.2011.08.020 |
|
Molnar P, Tapponnier P.
1975. Cenozoic tectonics of Asia:effects of a continental collision. Science, 189(4201): 419-426.
DOI:10.1126/science.189.4201.419 |
|
Peltzer G, Tapponnier P.
1988. Formation and evolution of strike-slip faults, rifts, and basins during the India-Asia collision:an experimental approach. J. Geophys. Res.:Solid Earth, 93(B12): 15085-15117.
DOI:10.1029/JB093iB12p15085 |
|
Ren J J, Xu X W, Sun X Z, et al.
2012. Geological and geophysical evidences of late Quaternary activity of the range-front fault along the mid-segment of the Longmen Shan thrust belt. Chinese J. Geophys. (in Chinese), 55(6): 1929-1941.
DOI:10.6038/j.issn.0001-5733.2012.06.014 |
|
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: B11409.
|
|
Shen Z K, Sun J B, Zhang P Z, et al.
2009. Slip maxima at fault junctions and rupturing of barriers during the 2008 Wenchuan earthquake. Nat. Geosci., 2(10): 718-724.
DOI:10.1038/ngeo636 |
|
Shi F, Li A, Yang S P, et al.
2013. Research on late quaternary activity of the southeastern segment of Ganzi-Yushu fault zone. Seismology & Geology (in Chinese), 35(1): 50-63.
|
|
Storti F, Holdsworth R E, Salvini F.
2003. Intraplate strike-slip deformation belts. Geological Society, London, Special Publications, 210: 1-14.
DOI:10.1144/GSL.SP.2003.210 |
|
Tang W Q, Liu Y P, Chen Z L, et al.
2004. The preliminary study of the tectonic activities along the boundary faults around the Minshan uplift, western Sichuan. Sedimentary Geology & Tethyan Geology (in Chinese), 24(4): 31-34.
|
|
Tapponnier P, Molnar P.
1977. Active faulting and tectonics in China. Journal of Geophysical Research, 82(20): 2905-2930.
DOI:10.1029/JB082i020p02905 |
|
Tapponnier P, Ryerson F J, Van der Woerd J, et al.
2001. Long-term slip rates and characteristic slip:keys to active fault behaviour and earthquake hazard. C. R. Acad. Sci.-Ser. ⅡA-Earth Planet. Sci., 333(9): 483-494.
|
|
Thatcher W.
2007. Microplate model for the present-day deformation of Tibet. J. Geophys. Res.:Solid Earth, 112: B01401.
|
|
Thompson L G, Yao T D, Davis M E, et al.
1997. Tropical climate instability:The last glacial cycle from a Qinghai-Tibetan ice core. Science, 276(5320): 1821-1825.
DOI:10.1126/science.276.5320.1821 |
|
Van Der Woerd J, Ryerson F J, Tapponnier P, et al.
2000. Uniform slip-rate along the Kunlun Fault:Implications for seismic behaviour and large-scale tectonics. Geophysical Research Letters, 27(16): 2353-2356.
DOI:10.1029/1999GL011292 |
|
Wang W T, Zhang P Z, Zheng D W, et al.
2014. Late Cenozoic tectonic deformation of the Haiyuan fault zone in the northeastern margin of the Tibetan Plateau. Earth Science Frontiers (in Chinese), 21(4): 266-274.
|
|
Wang Y Z, Wang E N, Shen Z K, et al.
2008. GPS-constrained inversion of present-day slip rates along major faults of the Sichuan-Yunnan region, China. Science China Earth Sciences, 51(9): 1267-1283.
DOI:10.1007/s11430-008-0106-4 |
|
Wang Y Z, Wang M, Shen Z K, et al.
2010. Inter-seismic deformation field of the Ganzi-Yushu fault before the 2010 Yushu earthqauke. Seismology & Geology (in Chinese), 33(3): 525-532.
|
|
Wen X Z, Jia J K, Pan M, et al.
1988. Recent slip rates, earthquake recurrence intervals and strong seismic hazards on the northwestern segment of the Xianshuihe fault zone. Earthquake Research in China, 2(4): 432-451.
|
|
Wen X Z, Xu X W, Zheng R Z, et al.
2003. Average slip-rate and recent large earthquake ruptures along the Ganzi-Yushu fault. Sci. China (Ser.D) (in Chinese), 33(B04): 199-208.
|
|
Wu J W, Huang X M, Xie F R.
2017. Deformation of the most recent co-seismic surface ruptures along the Garzê-Yushu fault zone (Dangjiang segment) and tectonic implications for the Tibetan Plateau. Acta Geologica Sinica (English Edition), 91(2): 443-454.
DOI:10.1111/acgs.2017.91.issue-2 |
|
Wu Z H, Long C X, Fan T Y, et al.
2015. The arc rotational-shear active tectonic system on the southeastern margin of Tibetan Plateau and its dynamic characteristics and mechanism. Geological Bulletin of China (in Chinese), 34(1): 1-31.
|
|
Xu X W, Wen X Z, Zheng R Z, et al.
2003a. Pattern of latest tectonic motion and its dynamics for active blocks in Sichuan-Yunnan region, China. Sci. China (Ser. D), 46(2): 210-226.
|
|
Xu X W, Tapponnier P, Van D W J, et al.
2005. Late Quaternary sinistral slip rate along the Altyn Tagh fault and its structural transformation model. Sci. China (Ser. D), 48(3): 384-397.
DOI:10.1360/02yd0436 |
|
Xu X W, Wen X Z, Ye J Q, et al.
2008. The MS8.0 wenchuan earthquake surface ruptures and its seismogenic structure. Seismology & Geology (in Chinese), 30(3): 597-629.
|
|
Yang N, Zhang Y Q.
2010. Tission-track dating for activity of the Longmenshan fault zone and uplifting of the western Sichuan Plateau. Journal of Geomechanics (in Chinese), 16(4): 359-371.
|
|
Yao T D, Thompson L G, Shi Y F, et al.
1997. Climate variation since the Last Interglaciation recorded in the Guliya ice core. Science China Earth Sciences, 40(6): 662-668.
DOI:10.1007/BF02877697 |
|
Yuan D Y, Lei Z S, He W G, et al.
2011. Textual research of Yushu earthquake in 1738 A. D. in Qinghai Province and discussion their relationship with Yushu earthquake MS7.1 in 2010. Progress in Geophys. (in Chinese), 26(6): 1950-1958.
DOI:10.3969/j.issn.1004-2903.2011.06.008 |
|
Zhang K L, Wei D P.
2011. Implications from the kinematic pattern of the Longmenshan region. Tectonophysics, 504(1-4): 57-64.
DOI:10.1016/j.tecto.2011.01.019 |
|
Zhang P Z, Sheng Z K, Wang M, et al.
2004. Continuous deformation of the Tibetan Plateau from global positioning system data. Geology, 32(9): 809-812.
DOI:10.1130/G20554.1 |
|
Zhang P Z, Li C Y, Mao F Y.
2008. Strath terrace formation and strike-slip faulting. Seismology & Geology (in Chinese), 30(1): 44-57.
|
|
Zhao G, Li D, Wei W, et al. 1992. The late Quaternary slip rate and segmentation of the Xianshuihe active fault zone.//Proceedings of the PRC-USA Bilateral Symposium on the Xianshuihe fault zone. Beijing:Seismological Press, 41-57.
|
|
Zhao X L, Deng Q D, Chen S F.
1994. Tectonic geomorphology of the Minshan Uplift in western Sichuan, southwestern China. Seismology & Geology (in Chinese), 16(4): 429-439.
|
|
Zhou R J, Ma S H, Cai C X.
1996. Late Quaternary active features of the Ganzi-Yushu fault zone. Earthquake Research in China (in Chinese), 12(3): 250-260.
|
|
Zhou R J, Wen X Z, Cai C X, et al.
1997. Recent earthquakes and assessment of seismic tendency on the Garzê-Yushu fault zone. Seismology & Geology (in Chinese), 19(2): 115-124.
|
|
Zhou R J, Li Y, Densmore A L, et al.
2006. Active tectonics of the eastern margin of the Tibet Plateau. J. Mineral. Petrol. (in Chinese), 26(2): 40-51.
|
|
邓起东, 陈立春, 冉勇康.
2004. 活动构造定量研究与应用. 地学前缘, 11(4): 383–392.
|
|
杜义, 谢富仁, 张效亮, 等.
2009. 汶川MS8.0级地震断层滑动机制研究. 地球物理学报, 52(2): 464–473.
|
|
国家地震局震害防御司. 1995.
中国历史强震目录. 北京: 地震出版社: 277-278.
|
|
来庆洲, 丁林, 王宏伟, 等.
2006. 青藏高原东部边界扩展过程的磷灰石裂变径迹热历史制约. 中国科学:地球科学, 36(9): 785–796.
|
|
李安, 王晓先, 张世民等. 2016. 酒西盆地内部白杨河逆断裂-褶皱构造变形. //地壳构造与地壳应力文集28. 北京: 地震出版社, 1-11.
|
|
李陈侠, 徐锡伟, 闻学泽, 等.
2009. 东昆仑断裂东段玛沁-玛曲段几何结构特征. 地震地质, 31(3): 441–458.
|
|
李强, 江在森, 武艳强, 等.
2013. 海原-六盘山断裂带现今构造变形特征. 大地测量与地球动力学, 33(2): 18–22.
|
|
李闽峰, 邢成起, 蔡长星, 等.
1995. 玉树断裂活动性研究. 地震地质, 17(3): 218–224.
|
|
李彦宝, 王虎, 刘华国, 等.
2016. 甘孜-玉树断裂当江段地表破裂遗迹与1738年玉树西北地震的讨论. 震灾防御技术, 11(2): 207–217.
DOI:10.11899/zzfy20160204 |
|
任俊杰, 徐锡伟, 孙鑫喆, 等.
2012. 龙门山推覆构造带中段山前断裂晚第四纪活动的地质与地球物理证据. 地球物理学报, 55(6): 1929–1941.
DOI:10.6038/j.issn.0001-5733.2012.06.014 |
|
石峰, 李安, 杨晓平, 等.
2013. 甘孜-玉树断裂带东南段晚第四纪活动性研究. 地震地质, 35(1): 50–63.
|
|
唐文清, 刘宇平, 陈智梁, 等.
2004. 岷山隆起边界断裂构造活动初步研究. 沉积与特提斯地质, 24(4): 31–34.
|
|
王伟涛, 张培震, 郑德文, 等.
2014. 青藏高原东北缘海原断裂带晚新生代构造变形. 地学前缘, 21(4): 266–274.
|
|
王阎昭, 王恩宁, 沈正康, 等.
2008. 基于GPS资料约束反演川滇地区主要断裂现今活动速率. 中国科学:地球科学, 38(5): 582–597.
|
|
王阎昭, 王敏, 沈正康, 等.
2010. 2010年玉树地震震前甘孜-玉树断裂形变场分析. 地震地质, 33(3): 525–532.
|
|
闻学泽, 徐锡伟, 郑荣章, 等.
2003. 甘孜-玉树断裂的平均滑动速率与近代大地震破裂. 中国科学, 33(B04): 199–208.
|
|
吴中海, 龙长兴, 范桃园, 等.
2015. 青藏高原东南缘弧形旋扭活动构造体系及其动力学特征与机制. 地质通报, 34(1): 1–31.
|
|
徐锡伟, 闻学泽, 郑荣章, 等.
2003a. 川滇地区活动块体最新构造变动样式及其动力来源. 中国科学, 33(B04): 151–162.
|
|
徐锡伟, TapponnierP, Van Der Woerd J, 等.
2003b. 阿尔金断裂带晚第四纪左旋走滑速率及其构造运动转换模式讨论. 中国科学, 33(10): 967–974.
|
|
徐锡伟, 闻学泽, 叶建青, 等.
2008. 汶川MS8.0地震地表破裂带及其发震构造. 地震地质, 30(3): 597–629.
|
|
杨农, 张岳桥.
2010. 龙门山断裂活动和川西高原隆升历史的裂变径迹测年. 地质力学学报, 16(4): 359–371.
|
|
姚檀栋, ThompsonL G, 施雅风, 等.
1997. 古里雅冰芯中末次间冰期以来气候变化记录研究. 中国科学:地球科学, 27(5): 447–452.
|
|
袁道阳, 雷中生, 何文贵, 等.
2011. 1738年青海玉树地震考证及其与2010年玉树7. 1级地震关系探讨.地球物理学进展, 26(6): 1950–1958.
DOI:10.3969/j.issn.1004-2903.2011.06.008 |
|
张培震, 李传友, 毛凤英.
2008. 河流阶地演化与走滑断裂滑动速率. 地震地质, 30(1): 44–57.
|
|
赵小麒, 邓起东, 陈社发.
1994. 岷山隆起的构造地貌学研究. 地震地质, 16(4): 428–439.
|
|
周荣军, 马声浩, 蔡长星.
1996. 甘孜-玉树断裂带的晚第四纪活动特征. 中国地震, 12(3): 250–260.
|
|
周荣军, 闻学泽, 蔡长星, 等.
1997. 甘孜-玉树断裂带的近代地震与未来地震趋势估计. 地震地质, 19(2): 115–124.
|
|
周荣军, 李勇, DensmoreA L, 等.
2006. 青藏高原东缘活动构造. 矿物岩石, 26(2): 40–51.
|
|