周文英 梁明劍 魏婭玲 宮 悅 廖 程 左 洪 楊 力 薛 蓮 毛澤斌

1）四川省地震局, 成都 610041

2）中國地震局成都青藏高原地震研究所, 成都 610041

3）成都市永利聯(lián)科技有限公司, 成都 610000

引言

青藏高原東緣地質(zhì)構(gòu)造復(fù)雜，地形地貌反差顯著，具有高海拔、高原面和深切河谷共存的特點(diǎn)，是青藏高原向東擴(kuò)展生長活動構(gòu)造與地貌特征最復(fù)雜的地方，在橫穿龍門山50～60 km 寬的地勢斜坡帶構(gòu)成了東緣地區(qū)高差達(dá)4 500～5 000 m 的陡峻地貌邊界帶（張?jiān)罉虻龋?008）。巴顏喀拉塊體自西向東或南東向運(yùn)動過程中受到華南地塊的強(qiáng)烈阻擋，導(dǎo)致青藏高原東緣龍門山斷裂帶到龍日壩斷裂帶之間的整體縮短和隆升（徐錫偉等，2008）。東緣地區(qū)相繼發(fā)生了2008 年汶川8.0 級地震、2013 年蘆山7.0 級地震和2017 年九寨溝7.0 級地震，地震造成的隆升和強(qiáng)烈破壞對青藏高原東緣龍門山地貌具有重塑作用（Parker 等，2011；李勇等，2015；鄭立龍，2015；高青等，2018）。

青藏高原東緣主要發(fā)育有岷江、大渡河、安寧河、涪江等水系，這些水系可分為貫通型水系和山前型水系，均發(fā)源于高原內(nèi)部，在東緣地區(qū)形成深切河谷，并最終匯入長江。河流的演化受構(gòu)造運(yùn)動和氣候的共同影響和制約，同時(shí)地表過程也對構(gòu)造運(yùn)動和氣候具有反饋與記錄作用（賈營營等，2010），尤其是在構(gòu)造活躍地區(qū)，流域盆地的發(fā)育對構(gòu)造隆升、氣候和地形地貌等控制因素的變化最敏感（Clark 等，2000；Kirby 等，2003；張會平，2006；劉靜等，2009）。流域內(nèi)不同級別水系的發(fā)育與演化不僅記錄了區(qū)域內(nèi)不同程度的構(gòu)造活動，也記載了氣候變化、海平面升降等外部變化信息。至今，國內(nèi)外學(xué)者利用典型流域地貌參數(shù)對青藏高原東緣，尤其是跨龍門山邊界地帶的流域地貌開展了許多研究（賈營營等，2010；Kirby 等，2011；李?yuàn)^生等，2015；閆亮等，2018；梁歐博等，2018；李明等，2019；李偉等，2021），研究結(jié)果表明，流域盆地的發(fā)育受區(qū)域構(gòu)造活動控制，而面積-高程積分、河道陡峭指數(shù)等地貌指數(shù)也對區(qū)內(nèi)構(gòu)造活動具有明顯的指示作用。然而，這些研究主要集中在青藏高原東緣邊界地帶-龍門山構(gòu)造帶中，對向高原延伸的過渡地帶和腹地構(gòu)造地貌的研究仍較少。

2022 年6 月10 日，四川省馬爾康市發(fā)生6.0 級震群地震，地震發(fā)生在強(qiáng)烈隆升的龍門山及其西緣高原腹地的過渡地帶。區(qū)內(nèi)發(fā)育的主要河流為大渡河支流-腳木足河，河流發(fā)源于4 000～4 500 m 的川西北高原，走向NW 向，展布方向NW-NE 向，流經(jīng)寬谷地貌、落差超1 000 m 的高山峽谷后匯入大金川河。馬爾康震群位于NE 和NW 向的斷裂交匯區(qū)域，斷裂構(gòu)造相對復(fù)雜，中小地震活動較活躍。本文采用SRTM 數(shù)據(jù)，利用數(shù)字地貌分析方法提取腳木足河流域坡度、地形起伏度、面積-高程積分和河道陡峭指數(shù)等地貌參數(shù)，同時(shí)結(jié)合區(qū)域降水量和巖性數(shù)據(jù)，分析流域盆地構(gòu)造地貌特征及其控制因素，進(jìn)一步探討研究區(qū)地震活動-構(gòu)造變形-地貌變化之間的關(guān)系。

1 研究區(qū)概況

1.1 地震構(gòu)造環(huán)境

青藏高原東緣是地殼活動與地震活動強(qiáng)烈的地區(qū)之一，其中，巴顏喀拉塊體向東運(yùn)動，受到華南塊體的阻擋，在龍門山地區(qū)形成了宏大的地貌與構(gòu)造邊界帶，也是近年來我國大陸強(qiáng)震主體活動地區(qū)之一。研究區(qū)位于巴顏喀拉塊體內(nèi)部（圖1（a）），塊體以龍日壩斷裂為界，東側(cè)為龍門山次級塊體，構(gòu)造隆升強(qiáng)烈，河谷深切，地形落差大，如圖1（b）中的地形剖面A-A’所示；西側(cè)阿壩次級塊體構(gòu)造變形以整體抬升為主，差異性隆升不強(qiáng)，僅表現(xiàn)在局部，如圖1（c）中的地形剖面B-B’所示（徐錫偉等，2008；陳長云等，2013；詹艷等，2021）。玉樹-甘孜-鮮水河斷裂帶、東昆侖斷裂帶和龍門山斷裂帶分別構(gòu)成了巴顏喀拉塊體的南邊界、北邊界和東邊界，邊界斷裂也是強(qiáng)震、大地震孕育與發(fā)生最活躍的地方，已有學(xué)者圍繞這些邊界斷裂開展了詳細(xì)研究（Royden 等，1997；Meng 等，2006；張培震等，2009；Hubbard 等，2009；Fu 等，2011），關(guān)于塊體東緣的變形機(jī)制衍生出了眾多的構(gòu)造變形和演化模式。除邊界斷裂外，巴顏喀拉塊體內(nèi)部還發(fā)育有眾多的活動構(gòu)造，位于塊體東部的馬爾康6.0 級震群所在區(qū)域主要以NE 向與NW 向活動斷裂為主，包括NW 向的松崗斷裂、達(dá)日斷裂、阿壩斷裂及NE 向的龍日壩斷裂。其中，達(dá)日斷裂為全新世活動斷裂，斷裂中段仍保留有1947 年達(dá)日7?級地震的地表破裂帶，斷裂晚第四紀(jì)以來的平均滑動速率為（3.0±0.6）mm/a（梁明劍等，2014，2020）。阿壩斷裂為瑪多-甘德斷裂的東段，也是全新世活動斷裂（熊仁偉等，2010），筆者在阿壩盆地附近發(fā)現(xiàn)斷裂上仍保留有地震地表破裂帶。根據(jù)最新地震地質(zhì)調(diào)查結(jié)果，松崗斷裂為晚更新世活動斷裂（孫東等，2010），尤以馬爾康震群以南段的構(gòu)造地貌最為顯著。龍日壩斷裂為全新世活動斷裂（徐錫偉等，2008；Ren 等，2013；何建軍等，2016），可分為北段、中段和南段，均表現(xiàn)為右旋走滑兼逆沖的運(yùn)動性質(zhì)（徐錫偉等，2008），地震地質(zhì)獲得斷裂晚更新世以來的平均右旋走滑速率為（5.4±2.0）mm/a，垂直滑動速率約為0.7 mm/a（徐錫偉等，2008），南段全新世平均滑動速率為0.6 mm/a（何建軍等，2016）。馬爾康6.0 級震群發(fā)生在龍日壩斷裂南段與松崗斷裂交匯處附近，斷裂構(gòu)造相對復(fù)雜，次級斷裂較發(fā)育，也是中小地震活躍的地段。

圖1 研究區(qū)位置及周邊主要活動斷裂分布Fig.1 The distribution map of the main active faults, earthquakes and the location of the study area

1.2 地形地貌概況

研究區(qū)跨青藏高原東部腹地的丘狀高原地貌區(qū)與東緣的高山峽谷深切割地貌區(qū)，地貌單元較復(fù)雜，地勢總體西北高、東南低，東緣甚至在龍門山與四川盆地之間形成落差達(dá)4 500～5 000 m 的地貌邊界帶（圖1（a））。其中，高原東部腹地主要為研究區(qū)水系發(fā)源地，為水系上游流經(jīng)地段，以丘狀山體、嶺平谷寬、“U”形河谷等地貌為特征，地形起伏度??；東緣地區(qū)為水系出山河段或下游河段，地形地貌以高山峽谷、“V”形河谷、高聳山體為特征，地形起伏大，落差顯著，也是大地震導(dǎo)致地貌重塑的主要地帶。

馬爾康震群位于高原腹地與龍門山過渡地帶，區(qū)域內(nèi)發(fā)育的主干河流為腳木足河，該河流是大渡河的主要支流之一，發(fā)源于青海省巴顏喀拉山脈，先后流經(jīng)青海省的久治縣、班瑪縣，四川省的壤塘縣、阿壩縣、馬爾康市，在馬爾康市與梭磨河匯合后流入大金川河（圖2）。腳木足河干流的總體走向?yàn)镹W 向，支流眾多，呈羽毛狀分布；上游河流淺切于高原面上，河谷寬淺，水流緩慢而曲折，具有海拔高、起伏度低、河流下切有限的特征；下游河谷坡度陡增，河流下切較深，形成高差超過1 000 m 的陡峭峽谷。

圖2 研究區(qū)地形地貌Fig.2 The topography map showing the rives, active faults and earthquake in the study area

2 研究數(shù)據(jù)與方法

2.1 數(shù)據(jù)來源與處理

SRTM 90 m 分辨率高程數(shù)據(jù)來源于地理空間數(shù)據(jù)云（http://www.gscloud.cn/）。利用下載的utm_srtm_57_06 數(shù)據(jù)提取河流和流域分布，從中提取出研究區(qū)范圍-馬爾康6.0 級震群所在的腳木足河流域。巖性數(shù)據(jù)主要來源于1∶50 萬區(qū)域地質(zhì)圖，經(jīng)過矢量化獲取了研究區(qū)的巖性數(shù)據(jù)。降雨量數(shù)據(jù)來源于國家科技基礎(chǔ)條件平臺-國家地球系統(tǒng)科學(xué)數(shù)據(jù)中心（http://www.geodata.cn），數(shù)據(jù)集名稱為中國0.25°分辨率年降水量數(shù)據(jù)集（1948-2016 年），通過求解這68 年的平均值獲得年平均降雨量，最后裁剪出研究區(qū)范圍。為確保多源空間數(shù)據(jù)能夠疊加展示，所有數(shù)據(jù)統(tǒng)一采用WGS_1984_UTM_Zone_47N 坐標(biāo)系。

2.2 原理與方法

2.2.1 坡度與地形起伏度

坡度與地形起伏度是地貌分析中最常見的參數(shù)，其對地貌的發(fā)育程度和地貌成因具有重要指示作用。其中，地表任一點(diǎn)的坡度是指過該點(diǎn)切平面與水平面之間的夾角。地形起伏度（RA）又稱相對高度，其能夠定量判斷流域盆地的切割剝蝕程度（梁歐博等，2018），可深刻表征區(qū)域構(gòu)造活動強(qiáng)度差異，常用于造山帶、高原山脈等發(fā)育演化特征研究（張會平，2006）。一般而言，在構(gòu)造活動發(fā)育、剝蝕強(qiáng)烈的地區(qū)，其區(qū)域地形起伏度較大（蘇琦，2015）。地形起伏度表達(dá)式為：

其中，Hmax和Hmin分別為一定統(tǒng)計(jì)單元內(nèi)的最高和最低高程。

本文利用ArcGIS 柵格表面分析中的坡度功能提取研究區(qū)坡度，利用ArcGIS 鄰域分析中的焦點(diǎn)統(tǒng)計(jì)（Focal Statistics）功能計(jì)算地形起伏度，其中窗口大小設(shè)置為適宜分析青藏高原的12 網(wǎng)格×12 網(wǎng)格（韓海輝等，2012），統(tǒng)計(jì)類型選擇Range。

2.2.2 面積-高程積分

流域盆地面積-高程積分最先是由Strahler（1952）提出來的，其能夠指示地貌演化的發(fā)育階段，間接反映構(gòu)造活動抬升的差異情況（Pike 等，1971）。經(jīng)過不斷地發(fā)展，現(xiàn)已成為指示地貌發(fā)育階段及構(gòu)造活動的重要地貌參數(shù)，流域面積-高程積分HI值對數(shù)字高程模型（Digital Elevation Model, DEM）的分辨率不具有依賴性。求取HI值的方法主要有積分曲線法、體積比例法及起伏比法（Pike 等，1971），各方法計(jì)算結(jié)果一致，其中起伏比法是最高效簡捷的方法（常直楊等，2015），在水系發(fā)育、流域數(shù)量較多的地區(qū)得到廣泛應(yīng)用（梁歐博等，2018；洪艷等，2019；周文英等，2022），本文采用該方法計(jì)算各亞流域盆地HI值：

其中，hmean、hmin和hmax分別為流域內(nèi)高程的平均值、最小值和最大值。

本文研究流程如下：（1）以ArcGIS 水文分析模塊（Hydrology）為基礎(chǔ)，利用空間建模工具，實(shí)現(xiàn)填洼、流向計(jì)算、匯流量計(jì)算、河網(wǎng)提取、河流Strahler 分級、流域盆地提取的半自動化提取，經(jīng)過多次試驗(yàn)，選擇與2017 年版天地圖（1∶100 萬）水系數(shù)據(jù)吻合度較好的河流水系作為分析基礎(chǔ)，進(jìn)而確定研究區(qū)提取水系的閾值（約10 km2），獲得相應(yīng)的亞流盆地。（2）由于HI具有明顯的空間尺度效應(yīng)，大量研究表明，較大的流域盆地HI值能夠較好地反映區(qū)域構(gòu)造活動性特征（陳彥杰等，2005），較小流域盆地的HI值更有可能反映局部構(gòu)造作用、地層巖性和地表過程的差異（梁明劍等，2014），甚至可能反映近期活動的褶皺構(gòu)造（Chen 等，2003）?？紤]到流域盆地HI值應(yīng)盡可能反映區(qū)域的構(gòu)造活動及亞流域盆地范圍的覆蓋度，最終選擇3 級以上河流對應(yīng)的246 個(gè)亞流域盆地（圖3，21 個(gè)3 級亞流域盆地，124 個(gè)4 級亞流域盆地，64 個(gè)5級亞流域盆地，37 個(gè)6 級亞流域盆地）作為分析對象，保證了亞流域盆地盡可能多地覆蓋全流域范圍。（3）利用ArcGIS 區(qū)域統(tǒng)計(jì)（Zonal Statistics）分析功能計(jì)算每個(gè)亞流域盆地的最大高程值、最小高程值和平均高程值，再利用式（2）求取HI值，同時(shí)將HI值賦值到相對應(yīng)的亞流域盆地幾何中心點(diǎn)。（4）采用反距離加權(quán)空間插值法得到區(qū)域HI等值分區(qū)圖，再利用研究區(qū)范圍圖裁剪出研究區(qū)HI等值分區(qū)圖。

圖3 亞流域分布Fig.3 Distribution map of sub-basins

2.2.3 河道陡峭指數(shù)

河道陡峭指數(shù)是可以反映區(qū)域構(gòu)造隆升的地貌參數(shù)，同時(shí)也是示蹤抬升速率較高地區(qū)或構(gòu)造較活躍地區(qū)的構(gòu)造活動參數(shù)（Kirby 等，2011）。在構(gòu)造活動強(qiáng)烈的地區(qū)，普遍發(fā)育了基巖河道或沉積物較薄的基巖水系，在這樣的條件下，河道中某點(diǎn)高程是由基巖抬升速率U和侵蝕速率E決定的。當(dāng)基巖河道處于均衡狀態(tài)時(shí)，河道中的高程將不隨時(shí)間而改變，此時(shí)認(rèn)為基巖抬升速率和侵蝕速率相等，此時(shí)河流比降S由下式求得（Howard 等，1983；Wang 等，2017）：

其中，K為侵蝕系數(shù)，A為流域面積，m為面積指數(shù)，n為坡度指數(shù)。

令ks= （U/K）1/n，?=m/n，則可得：

其中，ks為陡峭指數(shù)，?為河道凹曲度或曲率。

為使不同流域及流域面積指數(shù)下的陡峭指數(shù)具有可比性，一般需要先確定研究區(qū)域內(nèi)的參考凹度指數(shù)?ref，從而得到歸一化河道陡峭指數(shù)ksn。本文利用Schwanghart 等（2014）開發(fā)的TopoToolbox 2 工具，將選定的匯流面積在10 km2以上的所有河道坡度-面積（S-A）進(jìn)行線性回歸（圖4），得到?平均值為0.46，該值在0.35～0.65 范圍（Snyder 等，2000）內(nèi)，因此將0.46 作為研究區(qū)?ref。將平滑窗口設(shè)置為1 000 m，對河道數(shù)據(jù)進(jìn)行平滑處理進(jìn)而消除異常點(diǎn)（Whipple 等，2007），提取研究區(qū)內(nèi)所有河道的歸一化河道陡峭指數(shù)ksn。

圖4 流域河道坡度-面積線性回歸關(guān)系Fig.4 Slope versus area in a log-log plot

3 研究結(jié)果

3.1 坡度與地形起伏度

研究區(qū)坡度為0°～67.5°，為便于分析，使用自然裂點(diǎn)法將其劃分為5 個(gè)等級。其中，坡度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%。研究區(qū)坡度（圖5）與地形起伏度（圖6）呈現(xiàn)的分布特征相似，在龍日壩斷裂左側(cè)區(qū)域大致呈中間高、邊緣低的特征，在龍日壩斷裂右側(cè)區(qū)域大部分是高值區(qū)，僅北部梭磨河河源地區(qū)坡度和地形起伏度較低?？傮w來說，成片的高坡度、高起伏度區(qū)域主要集中在東南區(qū)域（圖7（a）、7（b）），低坡度、低起伏度區(qū)域主要集中在壤塘縣以東區(qū)域（圖7（c））和甘德南緣斷裂-阿壩盆中斷裂連線附近區(qū)域（圖7（d））。

圖5 研究區(qū)坡度分布Fig.5 The slope distribution map of study area

圖6 研究區(qū)地形起伏度分布Fig.6 The relief amplitude distribution map of study area

3.2 面積-高程積分

整個(gè)流域HI值為0.28～0.68（圖8），大部分亞流域HI值為0.35～0.6。從HI值分布特征來看，低值區(qū)集中在阿壩盆地等局部地區(qū)，次低值區(qū)集中在班瑪縣以北區(qū)域、阿壩縣附近及壤塘縣以東區(qū)域，大致沿班瑪縣-阿壩縣呈帶狀分布，高值區(qū)域主要沿松崗斷裂兩側(cè)分布，在達(dá)日斷裂和松崗斷裂夾持區(qū)域內(nèi)存在一片低HI值區(qū)域。依據(jù)河流走向，上游區(qū)域HI值普遍小于下游區(qū)域HI值。

圖8 研究區(qū)HI 分布Fig.8 HI values distribution map of study area

3.3 河道陡峭指數(shù)

為更清晰地展示ksn空間分布情況，將研究區(qū)ksn分為7 個(gè)等級（圖9）。研究區(qū)ksn呈現(xiàn)明顯的區(qū)域特征，即從NW 向到SE 向逐漸增大，其中ksn高值區(qū)主要集中在主河道下游河段，ksn次高值區(qū)主要集中在主河道中游河段，ksn低值區(qū)分布在主河道和支流河道上游河段。

圖9 研究區(qū)ksn 分布Fig.9 Spatial distribution map of channel steepness index

4 討論

4.1 降水對流域地貌的影響

分析面積-高程積分和河道陡峭指數(shù)空間分布差異時(shí)，需綜合考慮構(gòu)造活動、巖性和降水的影響。一般而言，降水量越大，侵蝕系數(shù)越大，如果基巖抬升速率一定，對應(yīng)的陡峭系數(shù)越?。簹W博等，2018）。研究區(qū)年均降水量為428～656 mm/a（圖10），且呈現(xiàn)由西北到東南增加的特征，龍日壩斷裂兩側(cè)的降水量呈現(xiàn)東多西少的差異。區(qū)域內(nèi)面積-高程積分和河道陡峭指數(shù)并不存在由西北到東南增加的性質(zhì)，且在龍日壩斷裂兩側(cè)均呈現(xiàn)高值分布。Liu 等（2011）的研究成果也表明，研究區(qū)內(nèi)的現(xiàn)代侵蝕速率相對較?。ǎ?.1 mm/a）。這些均表明降水并不是導(dǎo)致區(qū)域地貌指數(shù)呈現(xiàn)空間差異性的主要因素。

圖10 研究區(qū)年均降水量Fig.10 Spatial distribution of annual precipitation of study area

4.2 巖性對流域地貌的影響

研究區(qū)主要分布有三疊系，中統(tǒng)為板巖夾長石石英砂巖、頁巖，上統(tǒng)為淺變質(zhì)的長石石英砂巖；第四系集中分布在阿壩盆地，區(qū)內(nèi)少量夾雜分布有花崗巖（圖11）。依據(jù)梁歐博等（2018）對巖石堅(jiān)硬程度的劃分可知，區(qū)內(nèi)巖石主要為板巖和砂巖，抗侵蝕能力較強(qiáng)，河流侵蝕作用較弱。綜合分析認(rèn)為，研究區(qū)面積-高程積分在阿壩一帶呈現(xiàn)低值的原因可能與之發(fā)育的第四紀(jì)盆地有關(guān)，其余區(qū)域面積-高程積分和河道陡峭指數(shù)的空間差異與巖性的相關(guān)性較小。

圖11 研究區(qū)巖性分布Fig.11 Geological map of study area

4.3 斷裂構(gòu)造活動對流域地貌的影響

由于巖性和降雨不是導(dǎo)致地貌參數(shù)區(qū)域變化的主要因素，為此可用面積-高程積分和河道陡峭指數(shù)分析區(qū)域內(nèi)的構(gòu)造抬升速率，由此可知腳木足河下游較上游區(qū)域具有更高的基巖河道隆升速率。已有研究表明，達(dá)日、松崗和龍日壩斷裂具有晚第四紀(jì)活動性（徐錫偉等，2008；Ren 等，2013；梁明劍等，2014，2020；何建軍等，2016）。龍日壩斷裂帶兩側(cè)面積-高程積分和河道陡峭指數(shù)均呈現(xiàn)一定的高值，這表明這些區(qū)域存在構(gòu)造隆升作用；龍日壩斷裂帶東側(cè)呈現(xiàn)面積-高程積分和河道陡峭指數(shù)高值的原因可能是靠近龍門山斷裂帶，構(gòu)造活動變形強(qiáng)烈，差異性隆升明顯（徐錫偉等，2008；陳長云等，2013）；龍日壩斷裂帶西側(cè)呈現(xiàn)面積-高程積分和河道陡峭指數(shù)高值的原因可能與達(dá)日、松崗和龍日壩斷裂的共同作用有關(guān)。在構(gòu)造活躍地區(qū)，構(gòu)造活動是區(qū)域作用的主要力量，面積-高程積分＞0.43 的區(qū)域以隆升為主（洪艷等，2019），松崗斷裂與甘德南緣斷裂階區(qū)處存在的面積-高程積分高值區(qū)可能是2 條左旋走滑斷裂呈右階排列導(dǎo)致的擠壓隆升區(qū)，也可能存在局部復(fù)雜構(gòu)造。區(qū)域內(nèi)中強(qiáng)地震主要集中在龍日壩斷裂西側(cè)、面積-高程積分高值區(qū)中的低值帶內(nèi)，這可能與區(qū)域差異性構(gòu)造活動有關(guān)。

5 結(jié)論

馬爾康6.0 級震群發(fā)生在青藏高原東部腹地和龍門山過渡地帶，震區(qū)位于NE 向的龍日壩斷裂和NW 向的松崗斷裂、達(dá)日斷裂交匯處附近，斷裂構(gòu)造相對復(fù)雜，次級斷裂較發(fā)育，也是現(xiàn)今地震活躍的地段。震群所在區(qū)域主要發(fā)育的主干河流為腳木足河，河流整體呈NW 走向，支流較多，先后流經(jīng)川西北高原、高山峽谷地區(qū)，水系走向和展布方向受構(gòu)造控制。本文利用SRTM3 數(shù)據(jù)提取了馬爾康6.0 級震群所在腳木足河流域的坡度、起伏度、面積-高程積分和河道陡峭指數(shù)等地貌參數(shù)，發(fā)現(xiàn)腳木足河流域地貌參數(shù)分布特征與構(gòu)造變形、地震活動具有密切關(guān)系，主要獲得以下認(rèn)識：

（1）腳木足河流域坡度、起伏度、面積-高程積分和河道陡峭指數(shù)整體表現(xiàn)出西北低、東南高的相似變化規(guī)律，西北區(qū)域是河流的上游區(qū)域，屬于青藏高原東緣腹地，地貌參數(shù)呈現(xiàn)低值；東南區(qū)域是河流下游流入龍門山地區(qū)，屬于高山峽谷地貌，地貌參數(shù)呈現(xiàn)高值，這些結(jié)果表明下游區(qū)域較上游區(qū)域具有更大的地形起伏度、坡度，這與龍門山地區(qū)構(gòu)造活動與差異性隆升強(qiáng)烈有關(guān)。

（2）腳木足河流域坡度、起伏度、面積-高程積分和河道陡峭指數(shù)的高值區(qū)域主要分布在龍日壩斷裂兩側(cè)，沿松崗斷裂兩側(cè)展布，龍日壩斷裂東側(cè)為高原腹地向高山峽谷區(qū)轉(zhuǎn)換的區(qū)域，構(gòu)造活動變形強(qiáng)烈，差異性隆升明顯（徐錫偉等，2008；陳長云等，2013），應(yīng)該是各地貌參數(shù)呈現(xiàn)高值的主要原因；龍日壩斷裂西側(cè)是達(dá)日斷裂、松崗斷裂和龍日壩斷裂交匯的區(qū)域，構(gòu)造變形復(fù)雜，中強(qiáng)地震活躍，地貌破壞與重塑作用相對強(qiáng)烈，這可能是地貌參數(shù)在此處呈現(xiàn)高值的原因。值得注意的是面積-高程積分還在松崗斷裂與甘德南緣斷裂階區(qū)處存在一個(gè)高值區(qū)，推測其原因?yàn)? 條左旋走滑斷裂呈右階排列導(dǎo)致的擠壓隆升區(qū)或存在局部復(fù)雜構(gòu)造。

（3）研究區(qū)坡度、起伏度、面積-高程積分和河道陡峭指數(shù)等構(gòu)造地貌參數(shù)的空間分布差異與降水和巖性的關(guān)系較小，構(gòu)造活動是控制區(qū)域流域地貌空間差異的主要因素，這反映了巴顏喀拉塊體東緣向高原腹地過渡地帶的地貌對構(gòu)造變形仍存在顯著的反饋?zhàn)饔?，而青藏高原腹地的地貌與構(gòu)造變形的影響關(guān)系仍需進(jìn)一步研究。

致謝 審稿專家對本文提出了建設(shè)性的修改意見，在此表示衷心感謝。感謝國家科技基礎(chǔ)條件平臺-國家地球系統(tǒng)科學(xué)數(shù)據(jù)中心（http://www.geodata.cn）提供數(shù)據(jù)支撐。