• ISSN 1673-5722
  • CN 11-5429/P

四川马尔康6.0级震群区流域地貌特征与控制因素研究

周文英 梁明剑 魏娅玲 宫悦 廖程 左洪 杨力 薛莲 毛泽斌

周文英,梁明剑,魏娅玲,宫悦,廖程,左洪,杨力,薛莲,毛泽斌,2023. 四川马尔康6.0级震群区流域地貌特征与控制因素研究. 震灾防御技术,18(4):745−756. doi:10.11899/zzfy20230409. doi: 10.11899/zzfy20230409
引用本文: 周文英,梁明剑,魏娅玲,宫悦,廖程,左洪,杨力,薛莲,毛泽斌,2023. 四川马尔康6.0级震群区流域地貌特征与控制因素研究. 震灾防御技术,18(4):745−756. doi:10.11899/zzfy20230409. doi: 10.11899/zzfy20230409
Zhou Wenying, Liang Mingjian, Wei Yaling, Gong Yue, Liao Cheng, Zuo Hong, Yang Li, Xue Lian, Mao Zebin. Regional Geomorphic Characteristics and its factors Controlling of the Maerkang M6.0 Earthquake Swarm in Sichuan Province, Eastern Tibetan Plateau[J]. Technology for Earthquake Disaster Prevention, 2023, 18(4): 745-756. doi: 10.11899/zzfy20230409
Citation: Zhou Wenying, Liang Mingjian, Wei Yaling, Gong Yue, Liao Cheng, Zuo Hong, Yang Li, Xue Lian, Mao Zebin. Regional Geomorphic Characteristics and its factors Controlling of the Maerkang M6.0 Earthquake Swarm in Sichuan Province, Eastern Tibetan Plateau[J]. Technology for Earthquake Disaster Prevention, 2023, 18(4): 745-756. doi: 10.11899/zzfy20230409

四川马尔康6.0级震群区流域地貌特征与控制因素研究

doi: 10.11899/zzfy20230409
基金项目: 四川省自然科学基金(2022NSFSC0211);中国地震局地震预测研究所基本科研业务费(CEAIEF2022050502);国家重点研发计划(2021YFC3000601);西藏拉萨地球物理国家野外科学观测站研究课题(NORSLS2104)
详细信息
    作者简介:

    周文英,女,生于1988年。硕士,工程师。主要从事构造地貌方面的研究。E-mail:zhouwenying10@163.com

    通讯作者:

    梁明剑,男,生于1979年。博士,高级工程师。主要从事活动构造方面的研究。E-mail:23800794@qq.com

Regional Geomorphic Characteristics and its factors Controlling of the Maerkang M6.0 Earthquake Swarm in Sichuan Province, Eastern Tibetan Plateau

  • 摘要: 2022年6月10日四川马尔康6.0级震群地震发生在青藏高原东缘,位于强烈隆升的龙门山及其西缘高原腹地过渡地带,断裂构造复杂,水系发育,小震与中强地震活动频发。基于90 m×90 m空间分辨率的SRTM3数据,利用数字地貌分析方法提取了马尔康6.0级震群所在的脚木足河流域坡度、起伏度、面积-高程积分和河道陡峭指数等地貌参数,结合降水和岩性数据进行综合分析。研究结果表明,脚木足河流域坡度、起伏度、面积-高程积分和河道陡峭指数整体表现出西北低、东南高的特征。构造活动是导致区域地貌参数呈现差异性的主要原因,地貌参数高值主要分布于龙日坝断裂两侧,沿松岗断裂两侧展布,龙日坝断裂东侧是构造变形与差异性隆升强烈的龙门山地区;龙日坝断裂西侧是多条断裂交汇的构造复杂区,中强地震活跃,地貌破坏与重塑作用相对强烈,这可能是地貌参数高值的主要原因。
  • 青藏高原东缘地质构造复杂,地形地貌反差显著,具有高海拔、高原面和深切河谷共存的特点,是青藏高原向东扩展生长活动构造与地貌特征最复杂的地方,在横穿龙门山50~60 km宽的地势斜坡带构成了东缘地区高差达4 500~5 000 m的陡峻地貌边界带(张岳桥等,2008)。巴颜喀拉块体自西向东或南东向运动过程中受到华南地块的强烈阻挡,导致青藏高原东缘龙门山断裂带到龙日坝断裂带之间的整体缩短和隆升(徐锡伟等,2008)。东缘地区相继发生了2008年汶川8.0级地震、2013年芦山7.0级地震和2017年九寨沟7.0级地震,地震造成的隆升和强烈破坏对青藏高原东缘龙门山地貌具有重塑作用(Parker等,2011李勇等,2015郑立龙,2015高青等,2018)。

    青藏高原东缘主要发育有岷江、大渡河、安宁河、涪江等水系,这些水系可分为贯通型水系和山前型水系,均发源于高原内部,在东缘地区形成深切河谷,并最终汇入长江。河流的演化受构造运动和气候的共同影响和制约,同时地表过程也对构造运动和气候具有反馈与记录作用(贾营营等,2010),尤其是在构造活跃地区,流域盆地的发育对构造隆升、气候和地形地貌等控制因素的变化最敏感(Clark等,2000Kirby等,2003张会平,2006刘静等,2009)。流域内不同级别水系的发育与演化不仅记录了区域内不同程度的构造活动,也记载了气候变化、海平面升降等外部变化信息。至今,国内外学者利用典型流域地貌参数对青藏高原东缘,尤其是跨龙门山边界地带的流域地貌开展了许多研究(贾营营等,2010Kirby等,2011李奋生等,2015闫亮等,2018梁欧博等,2018李明等,2019李伟等,2021),研究结果表明,流域盆地的发育受区域构造活动控制,而面积-高程积分、河道陡峭指数等地貌指数也对区内构造活动具有明显的指示作用。然而,这些研究主要集中在青藏高原东缘边界地带——龙门山构造带中,对向高原延伸的过渡地带和腹地构造地貌的研究仍较少。

    2022年6月10日,四川省马尔康市发生6.0级震群地震,地震发生在强烈隆升的龙门山及其西缘高原腹地的过渡地带。区内发育的主要河流为大渡河支流——脚木足河,河流发源于4 000~4 500 m的川西北高原,走向NW向,展布方向NW-NE向,流经宽谷地貌、落差超1 000 m的高山峡谷后汇入大金川河。马尔康震群位于NE和NW向的断裂交汇区域,断裂构造相对复杂,中小地震活动较活跃。本文采用SRTM数据,利用数字地貌分析方法提取脚木足河流域坡度、地形起伏度、面积-高程积分和河道陡峭指数等地貌参数,同时结合区域降水量和岩性数据,分析流域盆地构造地貌特征及其控制因素,进一步探讨研究区地震活动-构造变形-地貌变化之间的关系。

    青藏高原东缘是地壳活动与地震活动强烈的地区之一,其中,巴颜喀拉块体向东运动,受到华南块体的阻挡,在龙门山地区形成了宏大的地貌与构造边界带,也是近年来我国大陆强震主体活动地区之一。研究区位于巴颜喀拉块体内部(图1(a)),块体以龙日坝断裂为界,东侧为龙门山次级块体,构造隆升强烈,河谷深切,地形落差大,如图1(b)中的地形剖面A—A’所示;西侧阿坝次级块体构造变形以整体抬升为主,差异性隆升不强,仅表现在局部,如图1(c)中的地形剖面B—B’所示(徐锡伟等,2008陈长云等,2013詹艳等,2021)。玉树—甘孜—鲜水河断裂带、东昆仑断裂带和龙门山断裂带分别构成了巴颜喀拉块体的南边界、北边界和东边界,边界断裂也是强震、大地震孕育与发生最活跃的地方,已有学者围绕这些边界断裂开展了详细研究(Royden等,1997Meng等,2006张培震等,2009Hubbard等,2009Fu等,2011),关于块体东缘的变形机制衍生出了众多的构造变形和演化模式。除边界断裂外,巴颜喀拉块体内部还发育有众多的活动构造,位于块体东部的马尔康6.0级震群所在区域主要以NE向与NW向活动断裂为主,包括NW向的松岗断裂、达日断裂、阿坝断裂及NE向的龙日坝断裂。其中,达日断裂为全新世活动断裂,断裂中段仍保留有1947年达日7¾级地震的地表破裂带,断裂晚第四纪以来的平均滑动速率为(3.0±0.6)mm/a(梁明剑等,20142020)。阿坝断裂为玛多-甘德断裂的东段,也是全新世活动断裂(熊仁伟等,2010),笔者在阿坝盆地附近发现断裂上仍保留有地震地表破裂带。根据最新地震地质调查结果,松岗断裂为晚更新世活动断裂(孙东等,2010),尤以马尔康震群以南段的构造地貌最为显著。龙日坝断裂为全新世活动断裂(徐锡伟等,2008Ren等,2013何建军等,2016),可分为北段、中段和南段,均表现为右旋走滑兼逆冲的运动性质(徐锡伟等,2008),地震地质获得断裂晚更新世以来的平均右旋走滑速率为(5.4±2.0)mm/a,垂直滑动速率约为0.7 mm/a(徐锡伟等,2008),南段全新世平均滑动速率为0.6 mm/a(何建军等,2016)。马尔康6.0级震群发生在龙日坝断裂南段与松岗断裂交汇处附近,断裂构造相对复杂,次级断裂较发育,也是中小地震活跃的地段。

    图 1  研究区位置及周边主要活动断裂分布
    Figure 1.  The distribution map of the main active faults, earthquakes and the location of the study area

    研究区跨青藏高原东部腹地的丘状高原地貌区与东缘的高山峡谷深切割地貌区,地貌单元较复杂,地势总体西北高、东南低,东缘甚至在龙门山与四川盆地之间形成落差达4 500~5 000 m的地貌边界带(图1(a))。其中,高原东部腹地主要为研究区水系发源地,为水系上游流经地段,以丘状山体、岭平谷宽、“U”形河谷等地貌为特征,地形起伏度小;东缘地区为水系出山河段或下游河段,地形地貌以高山峡谷、“V”形河谷、高耸山体为特征,地形起伏大,落差显著,也是大地震导致地貌重塑的主要地带。

    马尔康震群位于高原腹地与龙门山过渡地带,区域内发育的主干河流为脚木足河,该河流是大渡河的主要支流之一,发源于青海省巴颜喀拉山脉,先后流经青海省的久治县、班玛县,四川省的壤塘县、阿坝县、马尔康市,在马尔康市与梭磨河汇合后流入大金川河(图2)。脚木足河干流的总体走向为NW向,支流众多,呈羽毛状分布;上游河流浅切于高原面上,河谷宽浅,水流缓慢而曲折,具有海拔高、起伏度低、河流下切有限的特征;下游河谷坡度陡增,河流下切较深,形成高差超过1 000 m的陡峭峡谷。

    图 2  研究区地形地貌
    Figure 2.  The topography map showing the rives, active faults and earthquake in the study area

    SRTM 90 m分辨率高程数据来源于地理空间数据云(http://www.gscloud.cn/)。利用下载的utm_srtm_57_06数据提取河流和流域分布,从中提取出研究区范围——马尔康6.0级震群所在的脚木足河流域。岩性数据主要来源于1∶50万区域地质图,经过矢量化获取了研究区的岩性数据。降雨量数据来源于国家科技基础条件平台——国家地球系统科学数据中心(http://www.geodata.cn),数据集名称为中国0.25°分辨率年降水量数据集(1948—2016年),通过求解这68年的平均值获得年平均降雨量,最后裁剪出研究区范围。为确保多源空间数据能够叠加展示,所有数据统一采用WGS_1984_UTM_Zone_47N坐标系。

    2.2.1   坡度与地形起伏度

    坡度与地形起伏度是地貌分析中最常见的参数,其对地貌的发育程度和地貌成因具有重要指示作用。其中,地表任一点的坡度是指过该点切平面与水平面之间的夹角。地形起伏度(RA)又称相对高度,其能够定量判断流域盆地的切割剥蚀程度(梁欧博等,2018),可深刻表征区域构造活动强度差异,常用于造山带、高原山脉等发育演化特征研究(张会平,2006)。一般而言,在构造活动发育、剥蚀强烈的地区,其区域地形起伏度较大(苏琦,2015)。地形起伏度表达式为:

    $$ RA=H_{{\rm{max}}}-H_{{\rm{min}}} $$ (1)

    其中,HmaxHmin分别为一定统计单元内的最高和最低高程。

    本文利用ArcGIS栅格表面分析中的坡度功能提取研究区坡度,利用ArcGIS邻域分析中的焦点统计(Focal Statistics)功能计算地形起伏度,其中窗口大小设置为适宜分析青藏高原的12网格×12网格(韩海辉等,2012),统计类型选择Range。

    2.2.2   面积-高程积分

    流域盆地面积-高程积分最先是由Strahler(1952)提出来的,其能够指示地貌演化的发育阶段,间接反映构造活动抬升的差异情况(Pike等,1971)。经过不断地发展,现已成为指示地貌发育阶段及构造活动的重要地貌参数,流域面积-高程积分HI值对数字高程模型(Digital Elevation Model, DEM)的分辨率不具有依赖性。求取HI值的方法主要有积分曲线法、体积比例法及起伏比法(Pike等,1971),各方法计算结果一致,其中起伏比法是最高效简捷的方法(常直杨等,2015),在水系发育、流域数量较多的地区得到广泛应用(梁欧博等,2018洪艳等,2019周文英等,2022),本文采用该方法计算各亚流域盆地HI值:

    $$ HI=(h_{{\rm{mean}}}-h_{{\rm{min}}})/(h_{{\rm{max}}}-h_{{\rm{min}}}) $$ (2)

    其中,hmeanhminhmax分别为流域内高程的平均值、最小值和最大值。

    本文研究流程如下:(1)以ArcGIS水文分析模块(Hydrology)为基础,利用空间建模工具,实现填洼、流向计算、汇流量计算、河网提取、河流Strahler分级、流域盆地提取的半自动化提取,经过多次试验,选择与2017年版天地图(1∶100万)水系数据吻合度较好的河流水系作为分析基础,进而确定研究区提取水系的阈值(约10 km2),获得相应的亚流盆地。(2)由于HI具有明显的空间尺度效应,大量研究表明,较大的流域盆地HI值能够较好地反映区域构造活动性特征(陈彦杰等,2005),较小流域盆地的HI值更有可能反映局部构造作用、地层岩性和地表过程的差异(梁明剑等,2014),甚至可能反映近期活动的褶皱构造(Chen等,2003)。考虑到流域盆地HI值应尽可能反映区域的构造活动及亚流域盆地范围的覆盖度,最终选择3级以上河流对应的246个亚流域盆地(图3,21个3级亚流域盆地,124个4级亚流域盆地,64个5级亚流域盆地,37个6级亚流域盆地)作为分析对象,保证了亚流域盆地尽可能多地覆盖全流域范围。(3)利用ArcGIS区域统计(Zonal Statistics)分析功能计算每个亚流域盆地的最大高程值、最小高程值和平均高程值,再利用式(2)求取HI值,同时将HI值赋值到相对应的亚流域盆地几何中心点。(4)采用反距离加权空间插值法得到区域HI等值分区图,再利用研究区范围图裁剪出研究区HI等值分区图。

    图 3  亚流域分布
    Figure 3.  Distribution map of sub-basins
    2.2.3   河道陡峭指数

    河道陡峭指数是可以反映区域构造隆升的地貌参数,同时也是示踪抬升速率较高地区或构造较活跃地区的构造活动参数(Kirby等,2011)。在构造活动强烈的地区,普遍发育了基岩河道或沉积物较薄的基岩水系,在这样的条件下,河道中某点高程是由基岩抬升速率U和侵蚀速率E决定的。当基岩河道处于均衡状态时,河道中的高程将不随时间而改变,此时认为基岩抬升速率和侵蚀速率相等,此时河流比降S由下式求得(Howard等,1983Wang等,2017):

    $$ S=(U/K)^{1/n}A^{-m/n } $$ (3)

    其中,K为侵蚀系数,A为流域面积,m为面积指数,n为坡度指数。

    ks= (U/K1/nƟ=m/n,则可得:

    $$ S=k_{{\rm{s}}} A^{-\theta } $$ (4)

    其中,ks为陡峭指数,Ɵ为河道凹曲度或曲率。

    为使不同流域及流域面积指数下的陡峭指数具有可比性,一般需要先确定研究区域内的参考凹度指数Ɵref,从而得到归一化河道陡峭指数ksn。本文利用Schwanghart等(2014)开发的TopoToolbox 2工具,将选定的汇流面积在10 km2以上的所有河道坡度-面积(S-A)进行线性回归(图4),得到Ɵ平均值为0.46,该值在0.35~0.65范围(Snyder等,2000)内,因此将0.46作为研究区Ɵref。将平滑窗口设置为1 000 m,对河道数据进行平滑处理进而消除异常点(Whipple等,2007),提取研究区内所有河道的归一化河道陡峭指数ksn

    Figure 4.  Slope versus area in a log-log plot
    图4 流域河道坡度-面积线性回归关系

    研究区坡度为0°~67.5°,为便于分析,使用自然裂点法将其划分为5个等级。其中,坡度0°~10.6°占总面积的16.06%,10.6°~18.5°占总面积的25.37%,18.5°~25.7°占总面积的26.80%,25.7°~33.4°占总面积的21.99%,33.4°~67.5°占总面积的9.78%。研究区坡度(图5)与地形起伏度(图6)呈现的分布特征相似,在龙日坝断裂左侧区域大致呈中间高、边缘低的特征,在龙日坝断裂右侧区域大部分是高值区,仅北部梭磨河河源地区坡度和地形起伏度较低。总体来说,成片的高坡度、高起伏度区域主要集中在东南区域(图7(a)、7(b)),低坡度、低起伏度区域主要集中在壤塘县以东区域(图7(c))和甘德南缘断裂——阿坝盆中断裂连线附近区域(图7(d))。

    图 5  研究区坡度分布
    Figure 5.  The slope distribution map of study area
    图 6  研究区地形起伏度分布
    Figure 6.  The relief amplitude distribution map of study area
    图 7  研究区典型流域地貌
    Figure 7.  Typical geomorphology in the study area

    整个流域HI值为0.28~0.68(图8),大部分亚流域HI值为0.35~0.6。从HI值分布特征来看,低值区集中在阿坝盆地等局部地区,次低值区集中在班玛县以北区域、阿坝县附近及壤塘县以东区域,大致沿班玛县—阿坝县呈带状分布,高值区域主要沿松岗断裂两侧分布,在达日断裂和松岗断裂夹持区域内存在一片低HI值区域。依据河流走向,上游区域HI值普遍小于下游区域HI值。

    图 8  研究区HI分布
    Figure 8.  HI values distribution map of study area

    为更清晰地展示ksn空间分布情况,将研究区ksn分为7个等级(图9)。研究区ksn呈现明显的区域特征,即从NW向到SE向逐渐增大,其中ksn高值区主要集中在主河道下游河段,ksn次高值区主要集中在主河道中游河段,ksn低值区分布在主河道和支流河道上游河段。

    图 9  研究区ksn分布
    Figure 9.  Spatial distribution map of channel steepness index

    分析面积-高程积分和河道陡峭指数空间分布差异时,需综合考虑构造活动、岩性和降水的影响。一般而言,降水量越大,侵蚀系数越大,如果基岩抬升速率一定,对应的陡峭系数越小(梁欧博等,2018)。研究区年均降水量为428~656 mm/a(图10),且呈现由西北到东南增加的特征,龙日坝断裂两侧的降水量呈现东多西少的差异。区域内面积-高程积分和河道陡峭指数并不存在由西北到东南增加的性质,且在龙日坝断裂两侧均呈现高值分布。Liu等(2011)的研究成果也表明,研究区内的现代侵蚀速率相对较小(<0.1 mm/a)。这些均表明降水并不是导致区域地貌指数呈现空间差异性的主要因素。

    图 10  研究区年均降水量
    Figure 10.  Spatial distribution of annual precipitation of study area

    研究区主要分布有三叠系,中统为板岩夹长石石英砂岩、页岩,上统为浅变质的长石石英砂岩;第四系集中分布在阿坝盆地,区内少量夹杂分布有花岗岩(图11)。依据梁欧博等(2018)对岩石坚硬程度的划分可知,区内岩石主要为板岩和砂岩,抗侵蚀能力较强,河流侵蚀作用较弱。综合分析认为,研究区面积-高程积分在阿坝一带呈现低值的原因可能与之发育的第四纪盆地有关,其余区域面积-高程积分和河道陡峭指数的空间差异与岩性的相关性较小。

    图 11  研究区岩性分布
    Figure 11.  Geological map of study area

    由于岩性和降雨不是导致地貌参数区域变化的主要因素,为此可用面积-高程积分和河道陡峭指数分析区域内的构造抬升速率,由此可知脚木足河下游较上游区域具有更高的基岩河道隆升速率。已有研究表明,达日、松岗和龙日坝断裂具有晚第四纪活动性(徐锡伟等,2008Ren等,2013梁明剑等,20142020何建军等,2016)。龙日坝断裂带两侧面积-高程积分和河道陡峭指数均呈现一定的高值,这表明这些区域存在构造隆升作用;龙日坝断裂带东侧呈现面积-高程积分和河道陡峭指数高值的原因可能是靠近龙门山断裂带,构造活动变形强烈,差异性隆升明显(徐锡伟等,2008陈长云等,2013);龙日坝断裂带西侧呈现面积-高程积分和河道陡峭指数高值的原因可能与达日、松岗和龙日坝断裂的共同作用有关。在构造活跃地区,构造活动是区域作用的主要力量,面积-高程积分>0.43的区域以隆升为主(洪艳等,2019),松岗断裂与甘德南缘断裂阶区处存在的面积-高程积分高值区可能是2条左旋走滑断裂呈右阶排列导致的挤压隆升区,也可能存在局部复杂构造。区域内中强地震主要集中在龙日坝断裂西侧、面积-高程积分高值区中的低值带内,这可能与区域差异性构造活动有关。

    马尔康6.0级震群发生在青藏高原东部腹地和龙门山过渡地带,震区位于NE向的龙日坝断裂和NW向的松岗断裂、达日断裂交汇处附近,断裂构造相对复杂,次级断裂较发育,也是现今地震活跃的地段。震群所在区域主要发育的主干河流为脚木足河,河流整体呈NW走向,支流较多,先后流经川西北高原、高山峡谷地区,水系走向和展布方向受构造控制。本文利用SRTM3数据提取了马尔康6.0级震群所在脚木足河流域的坡度、起伏度、面积-高程积分和河道陡峭指数等地貌参数,发现脚木足河流域地貌参数分布特征与构造变形、地震活动具有密切关系,主要获得以下认识:

    (1)脚木足河流域坡度、起伏度、面积-高程积分和河道陡峭指数整体表现出西北低、东南高的相似变化规律,西北区域是河流的上游区域,属于青藏高原东缘腹地,地貌参数呈现低值;东南区域是河流下游流入龙门山地区,属于高山峡谷地貌,地貌参数呈现高值,这些结果表明下游区域较上游区域具有更大的地形起伏度、坡度,这与龙门山地区构造活动与差异性隆升强烈有关。

    (2)脚木足河流域坡度、起伏度、面积-高程积分和河道陡峭指数的高值区域主要分布在龙日坝断裂两侧,沿松岗断裂两侧展布,龙日坝断裂东侧为高原腹地向高山峡谷区转换的区域,构造活动变形强烈,差异性隆升明显(徐锡伟等,2008陈长云等,2013),应该是各地貌参数呈现高值的主要原因;龙日坝断裂西侧是达日断裂、松岗断裂和龙日坝断裂交汇的区域,构造变形复杂,中强地震活跃,地貌破坏与重塑作用相对强烈,这可能是地貌参数在此处呈现高值的原因。值得注意的是面积-高程积分还在松岗断裂与甘德南缘断裂阶区处存在一个高值区,推测其原因为2条左旋走滑断裂呈右阶排列导致的挤压隆升区或存在局部复杂构造。

    (3)研究区坡度、起伏度、面积-高程积分和河道陡峭指数等构造地貌参数的空间分布差异与降水和岩性的关系较小,构造活动是控制区域流域地貌空间差异的主要因素,这反映了巴颜喀拉块体东缘向高原腹地过渡地带的地貌对构造变形仍存在显著的反馈作用,而青藏高原腹地的地貌与构造变形的影响关系仍需进一步研究。

    致谢 审稿专家对本文提出了建设性的修改意见,在此表示衷心感谢。感谢国家科技基础条件平台——国家地球系统科学数据中心(http://www.geodata.cn)提供数据支撑。

  • 图  1  研究区位置及周边主要活动断裂分布

    Figure  1.  The distribution map of the main active faults, earthquakes and the location of the study area

    图  2  研究区地形地貌

    Figure  2.  The topography map showing the rives, active faults and earthquake in the study area

    图  3  亚流域分布

    Figure  3.  Distribution map of sub-basins

    Figure  4.  Slope versus area in a log-log plot

    图  5  研究区坡度分布

    Figure  5.  The slope distribution map of study area

    图  6  研究区地形起伏度分布

    Figure  6.  The relief amplitude distribution map of study area

    图  7  研究区典型流域地貌

    Figure  7.  Typical geomorphology in the study area

    图  8  研究区HI分布

    Figure  8.  HI values distribution map of study area

    图  9  研究区ksn分布

    Figure  9.  Spatial distribution map of channel steepness index

    图  10  研究区年均降水量

    Figure  10.  Spatial distribution of annual precipitation of study area

    图  11  研究区岩性分布

    Figure  11.  Geological map of study area

  • 常直杨, 王建, 白世彪等, 2015. 面积高程积分值计算方法的比较. 干旱区资源与环境, 29(3): 171—175

    Chang Z. Y. , Wang J. , Bai S. B. , et al. , 2015. Comparison of hypsometric integral methods. Journal of Arid Land Resources and Environment, 29(3): 171—175. (in Chinese)
    陈长云, 任金卫, 孟国杰等, 2013. 巴颜喀拉块体东部活动块体的划分、形变特征及构造意义. 地球物理学报, 56(12): 4125—4141

    Chen C. Y. , Ren J. W. , Meng G. J. , et al. , 2013. Division, deformation and tectonic implication of active blocks in the eastern segment of Bayan Har block. Chinese Journal of Geophysics, 56(12): 4125—4141. (in Chinese)
    陈彦杰, 郑光佑, 宋国城, 2005. 面积尺度与空间分布对流域面积高度积分及其地质意义的影响. 地理学报, (39): 53—69

    Chen Y. J. , Zheng G. Y. , Song G. C. , 2005. Influence of area and space dependence for hypsometric integral and its geological implications. Journal of Geographical Science, (39): 53—69. (in Chinese)
    高青, 李奋生, 邵崇建, 2018. 汶川地震驱动的剥蚀作用及对龙门山地貌演化的影响. 成都理工大学学报(自然科学版), 45(3): 386—392

    Gao Q. , Li F. S. , Shao C. J. , 2018. The erosion driven by Wenchuan earthquake and its impact on the geomorphic evolution of Longmenshan, Sichuan, China. Journal of Chengdu University of Technology (Science & Technology Edition), 45(3): 386—392. (in Chinese)
    韩海辉, 高婷, 易欢等, 2012. 基于变点分析法提取地势起伏度——以青藏高原为例. 地理科学, 32(1): 101—104

    Han H. H. , Gao T. , Yi H. , et al. , 2012. Extraction of relief amplitude based on Change Point Method——A case study on the Tibetan Plateau. Scientia Geographica Sinica, 32(1): 101—104. (in Chinese)
    何建军, 任俊杰, 丁锐等, 2016. 青藏高原东缘龙日坝断裂带南段晚第四纪活动及其构造意义. 震灾防御技术, 11(4): 707—721

    He J. J. , Ren J. J. , Ding R. , et al. , 2016. Late quaternary activity of the Southern Segment of Longriba fault zone in Eastern Tibet and its tectonic implications. Technology for Earthquake Disaster Prevention, 11(4): 707—721. (in Chinese)
    洪艳, 赵银兵, 王运生等, 2019. 利用面积高程积分方法研究龙门山断裂带地貌隆升特征. 科学技术与工程, 19(30): 43—51

    Hong Y. , Zhao Y. B. , Wang Y. S. , et al. , 2019. Study on the geomorphic uplift characteristics of the Longmenshan fault zone using hypsometric integral method. Science Technology and Engineering, 19(30): 43—51. (in Chinese)
    贾营营, 付碧宏, 王岩等, 2010. 青藏高原东缘龙门山断裂带晚新生代构造地貌生长及水系响应. 第四纪研究, 30(4): 825—836

    Jia Y. Y. , Fu B. H. , Wang Y. , et al. , 2010. Late Cenozoic tectono-geomorphic growth and drainage response in the Longmen Shan fault zone, east margin of Tibet. Quaternary Sciences, 30(4): 825—836. (in Chinese)
    李奋生, 赵国华, 李勇等, 2015. 龙门山地区水系发育特征及其对青藏高原东缘隆升的指示. 地质论评, 61(2): 345—355

    Li F. S. , Zhao G. H. , Li Y. , et al. , 2015. The characteristics of drainage development in Longmen Mountains Area and its indication to the uplift of eastern margin of Qinghai-Xizang (Tibet) Plateau. Geological Review, 61(2): 345—355. (in Chinese)
    李明, 张世民, 姜大伟等, 2019. 龙门山中北段流域地貌特征及其构造意义. 震灾防御技术, 14(3): 640—651

    Li M. , Zhang S. M. , Jiang D. W. , et al. , 2019. Geomorphic features and tectonic significance of the middle and North Section of Longmenshan. Technology for Earthquake Disaster Prevention, 14(3): 640—651. (in Chinese)
    李伟, 谢超, 程宏宾等, 2021. 利用水系方位角和GPS数据研究龙门山后山断裂运动模式. 四川地震, (2): 7—10

    Li W. , Xie C. , Cheng H. B. , et al. , 2021. Study on the fault movement model of the backmountain of Longmen Mountain based on the azimuth angle and GPS data. Earthquake Research in Sichuan, (2): 7—10. (in Chinese)
    李勇, 周荣军, 赵国华等, 2015. 汶川地震驱动的隆升、剥蚀作用与龙门山地貌生长——以映秀红椿沟为例. 成都理工大学学报(自然科学版), 42(1): 5—17

    Li Y. , Zhou R. J. , Zhao G. H. , et al. , 2015. Uplift and erosion driven by Wenchuan earthquake and their effects on geomorphic growth of Longmen Mountains——A case study of Hongchun gully in Yingxiu, China. Journal of Chengdu University of Technology (Science & Technology Edition), 42(1): 5—17. (in Chinese).
    梁明剑, 周荣军, 闫亮等, 2014. 青海达日断裂中段构造活动与地貌发育的响应关系探讨. 地震地质, 36(1): 28—38

    Liang M. J. , Zhou R. J. , Yan L. , et al. , 2014. The relationships between neotectonic activity of the middle segment of Dari Fault and its geomorphological response, Qinghai Province, China. Seismology and Geology, 36(1): 28—38. (in Chinese)
    梁明剑, 杨耀, 杜方等, 2020. 青海达日断裂中段晚第四纪活动性与1947年M7¾地震地表破裂带再研究. 地震地质, 42(3): 703—714

    Liang M. J. , Yang Y. , Du F. , et al. , 2020. Late Quaternary activity of the central segment of the Dari Fault and restudy of the surface rupture zone of the 1974 M7¾ Dari earthquake, Qinghai Province. Seismology and Geology, 42(3): 703—714. (in Chinese)
    梁欧博, 任俊杰, 吕延武, 2018. 涪江流域河流地貌特征对虎牙断裂带活动性的响应. 地震地质, 40(1): 42—56

    Liang O. B. , Ren J. J. , Lü Y. W. , 2018. The response of fluvial geomorphologic characteristics of the Fujiang drainge basin to activity of the Huya fault zone. Seismology and Geology, 40(1): 42—56. (in Chinese)
    刘静, 曾令森, 丁林等, 2009. 青藏高原东南缘构造地貌、活动构造和下地壳流动假说. 地质科学, 44(4): 1227—1255

    Liu J. , Zeng L. S. , Ding L. , et al. , 2009. Tectonic geomorphology, active tectonics and lower crustal channel flow hypothesis of the southeastern Tibetan Plateau. Chinese Journal of Geology, 44(4): 1227—1255. (in Chinese)
    苏琦, 2015. 青藏高原东北缘典型流域地貌参数分析与构造变形探讨. 兰州: 中国地震局兰州地震研究所.

    Su Q. , 2015. Study on typical drainage basins along northeastern Tibetan Plateau and its tectonic deformation. Lanzhou: China Earthquake Administration Lanzhou Institute of Seismology. (in Chinese)
    孙东, 王道永, 吴德超等, 2010. 马尔康巴拉水电站近场区主要断裂活动性及对工程的影响. 工程地质学报, 18(6): 940—949

    Sun D. , Wang D. Y. , Wu D. C. , et al. , 2010. Activity and effect of main faults in near field of Bala Hydropower station in Maerkang. Journal of Engineering Geology, 18(6): 940—949. (in Chinese)
    熊仁伟, 任金卫, 张军龙等, 2010. 玛多—甘德断裂甘德段晚第四纪活动特征. 地震, 30(4): 65—73

    Xiong R. W. , Ren J. W. , Zhang J. L. , et al. , 2010. Late Quaternary active characteristics of the Gande Segment in the Maduo-Gande Fault zone. Earthquake, 30(4): 65—73. (in Chinese)
    徐锡伟, 闻学泽, 陈桂华等, 2008. 巴颜喀拉地块东部龙日坝断裂带的发现及其大地构造意义. 中国科学 D辑: 地球科学, 38(5): 529—542.

    Xu X. W. , Wen X. Z. , Chen G. H. , et al. , 2008. Discovery of the Longriba fault zone in eastern Bayan Har block, China and its tectonic implication. Science in China Series D: Earth Sciences, 51(9): 1209—1223.
    闫亮, 李勇, 赵国华等, 2018. 青藏高原东缘龙门山构造带晚第四纪构造隆升作用的河流地貌响应. 第四纪研究, 38(1): 232—246

    Yan L. , Li Y. , Zhao G. H. , et al. , 2018. The uplift process and the geomorphological features of the rivers system in the Longmen Shan since Late Quaternary. Quaternary Sciences, 38(1): 232—246. (in Chinese)
    詹艳, 梁明剑, 孙翔宇等, 2021.2021年5月22日青海玛多MS7.4地震深部环境及发震构造模式. 地球物理学报, 64(7): 2232—2252

    Zhan Y. , Liang M. J. , Sun X. Y. , et al. , 2021. Deep structure and seismogenic pattern of the 2021.5. 22 Madoi (Qinghai) MS7.4 earthquake.Chinese Journal of Geophysics,64(7):2232—2252.(in Chinese)
    张会平, 2006. 青藏高原东缘、东北缘典型地区晚新生代地貌过程研究[D]. 北京: 中国地质大学.

    Zhang H. P. , 2006. Study on late Cenozoic geomorphic processes of typical regions along the eastern and northeastern Tibetan margins. Beijing: China University of Geosciences. (in Chinese)
    张培震, 闻学泽, 徐锡伟等, 2009.2008年汶川8.0级特大地震孕育和发生的多单元组合模式. 科学通报, 54(7): 944—953 doi: 10.1360/csb2009-54-7-944

    Zhang P. Z. , Wen X. Z. , Xu X. W. , et al. , 2009. Tectonic model of the great Wenchuan earthquake of May 12, 2008, Sichuan, China. Chinese Science Bulletin, 54(7): 944—953. (in Chinese) doi: 10.1360/csb2009-54-7-944
    张岳桥, 杨农, 施炜等, 2008. 青藏高原东缘新构造及其对汶川地震的控制作用. 地质学报, 82(12): 1668—1678

    Zhang Y. Q. , Yang N. , Shi W. , et al. , 2008. Neotectonics of eastern Tibet and its control on the Wenchuan Earthquake. Acta Geologica sinica, 82(12): 1668—1678. (in Chinese)
    郑立龙, 2015. 汶川地震后绵远河流域地表隆升与侵蚀过程研究[D]. 成都: 成都理工大学.

    Zheng L. L. , 2015. The research of surface uplift and erosion process of Mianyuan river drainage basin after the Wenchuan Earthquake. Chengdu: Chengdu University of Technology. (in Chinese)
    周文英, 梁明剑, 左洪等, 2022. 濑溪河流域地貌对泸县MS6.0地震发震构造变形的响应及其指示作用. 地震学报, 44(2): 250—259

    Zhou W. Y. , Liang M. J. , Zuo H. , et al. , 2022. Response of the topography of Laixihe drainage to the structural deformation induced by MS6.0 Luxian earthquake and its indication. Acta Seismologica Sinica, 44(2): 250—259. (in Chinese)
    Chen Y. C. , Sung Q. C. , Cheng K. Y. , 2003. Along-strike variations of morphotectonic features in the Western Foothills of Taiwan: tectonic implications based on stream-gradient and hypsometric analysis. Geomorphology, 56(1—2): 109—137. doi: 10.1016/S0169-555X(03)00059-X
    Clark M. K. , Royden L. H. , 2000. Topographic ooze: building the eastern margin of Tibet by lower crustal flow. Geology, 28(8): 703—706. doi: 10.1130/0091-7613(2000)28<703:TOBTEM>2.0.CO;2
    Fu B. H. , Shi P. L. , Guo H. D. , et al. , 2011. Surface deformation related to the 2008 Wenchuan earthquake, and mountain building of the Longmen Shan, eastern Tibetan Plateau. Journal of Asian Earth Sciences, 40(4): 805—824. doi: 10.1016/j.jseaes.2010.11.011
    Howard A. D. , Kerby G. , 1983. Channel changes in Badlands. GSA Bulletin, 94(6): 739—752. doi: 10.1130/0016-7606(1983)94<739:CCIB>2.0.CO;2
    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
    Kirby E. , Whipple K. X. , Tang W. Q. , et al. , 2003. Distribution of active rock uplift along the eastern margin of the Tibetan Plateau: Inferences from bedrock channel longitudinal profiles. Journal of Geophysical Research: Solid Earth, 108(B4): 2217.
    Kirby E., Ouimet W., 2011. Tectonic geomorphology along the eastern margin of Tibet: Insights into the pattern and processes of active deformation adjacent to the Sichuan Basin. Geological Society, London, Special Publications, 353(1): 165—188.
    Liu-Zeng J., Wen L., Oskin M., et al., 2011. Focused modern denudation of the Longmen Shan margin, eastern Tibetan Plateau. Geochemistry, Geophysics, Geosystems, 12(11): Q11007.
    Meng Q. R. , Hu J. M. , Wang E. , et al. , 2006. Late Cenozoic denudation by large-magnitude landslides in the eastern edge of Tibetan Plateau. Earth and Planetary Science Letters, 243(1—2): 252—267. doi: 10.1016/j.jpgl.2005.12.008
    Parker R. N. , Densmore A. L. , Rosser N. J. , et al. , 2011. Mass wasting triggered by the 2008 Wenchuan earthquake is greater than orogenic growth. Nature Geoscience, 4(7): 449—452. doi: 10.1038/ngeo1154
    Pike R. J. , Wilson S. E. , 1971. Elevation-relief ratio, hypsometric integral, and geomorphic area-altitude analysis. GSA Bulletin, 82(4): 1079—1084. doi: 10.1130/0016-7606(1971)82[1079:ERHIAG]2.0.CO;2
    Ren J. J. , Xu X. W. , Yeats R. S. , et al. , 2013. Latest Quaternary paleoseismology and slip rates of the Longriba fault zone, eastern Tibet: Implications for fault behavior and strain partitioning. Tectonics, 32(2): 216—238. doi: 10.1002/tect.20029
    Royden L. H. , Burchfiel B. C. , King R. W. , et al. , 1997. Surface deformation and lower crustal flow in Eastern Tibet. Science, 276(5313): 788—790. doi: 10.1126/science.276.5313.788
    Schwanghart W. , Scherler D. , 2014. TopoToolbox 2 – MATLAB-based software for topographic analysis and modeling in Earth surface sciences. Earth Surface Dynamics, 2(1): 1—7. doi: 10.5194/esurf-2-1-2014
    Snyder N. P. , Whipple K. X. , Tucker G. E. , et al. , 2000. Landscape response to tectonic forcing: Digital elevation model analysis of stream profiles in the Mendocino triple junction region, northern California. Geological Society of America Bulletin, 112(8): 1250—1263. doi: 10.1130/0016-7606(2000)112<1250:LRTTFD>2.0.CO;2
    Strahler A. N. , 1952. Hypsometric (area-altitude) analysis of erosional topography. GSA Bulletin, 63(11): 1117—1142. doi: 10.1130/0016-7606(1952)63[1117:HAAOET]2.0.CO;2
    Wang Y. Z. , Zhang H. P. , Zheng D. W. , et al. , 2017. Coupling slope–area analysis, integral approach and statistic tests to steady-state bedrock river profile analysis. Earth Surface Dynamics, 5(1): 145—160 doi: 10.5194/esurf-5-145-2017
    Whipple K., Wobus C. W., Crosby B. T., et al., 2007. New tools for quantitative geomorphology: extraction and interpretation of stream profiles from digital topographic data [EB/OL]. (2007-01-01)[2023-11-10].https://www.researchgate.net/profile/Cameron-Wobus/publication/268056101_New_Tools_for_Quantitative_Geomorphology_Extraction_and_Interpretation_of_Stream_Profiles_from_Digital_Topographic_Data/links/54f9bf970cf29a9fbd7c51d4/New-Tools-for-Quantitative-Geomorphology-Extraction-and-Interpretation-of-Stream-Profiles-from-Digital-Topographic-Data.pdf.
  • 期刊类型引用(1)

    1. 李圣,徐雨帆,卿元华,高骏. 滇西澜沧江流域面积-高程积分及其地貌意义. 四川地震. 2025(01): 36-41 . 百度学术

    其他类型引用(0)

  • 加载中
图(11)
计量
  • 文章访问数:  133
  • HTML全文浏览量:  42
  • PDF下载量:  17
  • 被引次数: 1
出版历程
  • 收稿日期:  2022-08-03
  • 刊出日期:  2023-12-01

目录

/

返回文章
返回