劉 譞，林 舟，丁 超

浙江大學地球科學學院，杭州310012

1 引言

均衡狀態(tài)時的河道縱剖面表現(xiàn)為平滑的形態(tài)，然而自然界的河流常因區(qū)域構(gòu)造或環(huán)境的擾動而呈現(xiàn)出不均衡狀態(tài)（Duvall et al.,2004；Whipple,2004；Cyr et al.,2014）。裂點（knickpoint）則位于不均衡河道縱剖面上緩坡段與陡坡段的轉(zhuǎn)折處（Howard et al., 1994; Crosby and Whipple, 2006；Hayakawa and Oguchi, 2006；Wobus et al., 2006），其發(fā)育常受控于流域盆地內(nèi)巖性差異、構(gòu)造抬升差異、侵蝕基準面變化或氣候變化等控制因素（Duvall et al.,2004；Cyr et al.,2014）。裂點的發(fā)育在活動造山帶流域內(nèi)表現(xiàn)得尤為明顯（Wesnousky et al., 1999；Kirby et al., 2010）。在活動造山帶區(qū)域，構(gòu)造、氣候及侵蝕過程相互作用，影響活動造山帶區(qū)域河流地貌的演化（Peizhen et al.,2001；Snyder et al., 2003；Duvall et al., 2004；Egholm et al., 2009； Whipple, 2009； DiBiase et al., 2010，2011；Ouimet et al., 2010；Herman et al., 2013），為裂點的發(fā)育提供了有利條件。因此活動造山帶流域內(nèi)河道縱剖面的演化及裂點的發(fā)育和溯源遷移，是理解該區(qū)域地貌演化過程的重要依據(jù)（Duvall et al.,2004；Kirby and Whipple,2012；Cyr et al.,2014）。

岷江上游流域位于青藏高原東緣的活動造山帶區(qū)域，河流普遍處于不均衡狀態(tài)且普遍發(fā)育有裂點，記錄了區(qū)域的河流地貌發(fā)育過程和特征。前人對岷江上游流域的構(gòu)造活動和地貌演化做了大量的調(diào)查研究工作（Kirby et al.,2000,2011；張岳橋等，2005；Kong et al.,2011；Liu et al.,2011；Zhang et al., 2011；Li et al., 2012；Thompson et al.,2015；Airaghi et al., 2017），這些工作主要從河流下切、隆升階段、地貌特征、斷層分布、冰川作用等方面探討了岷江上游流域地貌特征及構(gòu)造活動的特征和規(guī)律。前人的研究表明了河流下切作用與區(qū)域階段性隆升的相關(guān)性，地貌參數(shù)與區(qū)域差異性隆升的關(guān)聯(lián)，但對于岷江上游流域區(qū)域構(gòu)造活動、巖性差異、冰川作用等因素對河流縱剖面形態(tài)特征和演化具體的影響研究還比較薄弱。

本文擬對岷江上游流域支流的裂點進行提取和分析。為客觀地提取裂點，先將河道縱剖面置于同比例下，再結(jié)合基于河流水力侵蝕模型的裂點提取方法，獲得了較為客觀統(tǒng)一的裂點分布。在此基礎(chǔ)上，根據(jù)裂點的空間分布特征，分析區(qū)域內(nèi)裂點的類型和形成原因，探討冰川作用、巖性、構(gòu)造等因素對研究區(qū)河流縱剖面形態(tài)演化的具體影響。通過對流域內(nèi)裂點的特征進行分析，有助于尋找裂點形成的控制因素，認識區(qū)域的地質(zhì)演化。

2 研究區(qū)概況

2.1 地質(zhì)概況

岷江上游流域位于青藏高原東緣造山帶，處于岷山隆起帶、龍門山構(gòu)造帶、松潘-甘孜褶皺帶的結(jié)合部位（圖1），是青藏高原強烈隆升的區(qū)域（Kirby et al., 2003；Zhang et al., 2004；Xu et al.,2009；Gray and Pysklywec, 2012）。流域內(nèi)有岷江斷裂帶，龍門山斷裂帶等主要的斷裂系統(tǒng)。岷江斷裂帶大致呈南北向沿岷江上游的西側(cè)分布，是傾向北西的高角度逆沖斷裂帶，晚新生代以來構(gòu)造活動強烈（張岳橋等，2012；李峰等，2018）。而龍門山斷裂帶位于岷江流域的南側(cè)，大致呈南西—北東向沿岷江延伸，為傾向北西并帶有一定右旋走滑分量的高角度逆沖斷裂帶，分為茂縣—汶川斷裂、映秀—北川斷裂和灌縣—安縣斷裂三大分支斷裂，第四紀以來構(gòu)造活動強烈，是地震頻發(fā)的區(qū)域（張培震等，2008；Zhang et al.,2009；Wang et al.,2014；譚錫斌等，2015）。研究區(qū)內(nèi)各地層單元的界線較明顯，巖性組成較復(fù)雜，北部和中部廣泛分布有中生界變質(zhì)巖，龍門山構(gòu)造帶地區(qū)主要出露前中生界深變質(zhì)巖，局部地區(qū)發(fā)育花崗巖侵入體。該區(qū)域特殊的地質(zhì)狀況和構(gòu)造背景，為具體分析裂點的成因提供了理想的場所。

2.2 地貌概況

岷江上游流域地形表現(xiàn)為高坡度、高地形起伏的特點，流域內(nèi)平均海拔3000 m 以上，最高海拔超過5000 m，岷山、龍門山與四川盆地之間構(gòu)成顯著的地形陡變帶，記錄了印度和歐亞板塊撞擊后高原東向擴張的巨大形變（England and Molnar,1990；Clark and Royden，2000；Holt et al.,2000），沿地形陡變帶的許多區(qū)域相對地形高差大于3000 m。流域內(nèi)高海拔區(qū)域發(fā)育有冰川地貌，以U 型谷為主要特征；低海拔區(qū)域則由于岷江水系長期的深切割作用，受到了較為強烈的河流侵蝕作用（Kirby et al.,2003;Kirby and Ouimet,2011；張岳橋等，2005；張會平等，2006），形成溝谷縱橫的河流侵蝕地貌，以V型谷為主要特點，這一特征明顯受控于青藏高原東緣的隆升所造成的岷江水系主河道及其支流侵蝕作用的增強。

圖1 研究區(qū)分辨率30 m的DEM及河流、斷裂分布圖（疊溪為1933年疊溪古地震震中，映秀為2008年汶川地震震中，斷裂分布據(jù)Kirby et al.,2011）Fig.1 Map showing the topography,rivers,and faults in the study area（Diexi was the epicenter of the 1933 Diexi ancient earthquake,and Yingxiu was the epicenter of the 2008 Wenchuan earthquake,faults from Kirby et al.,2011）

2.3 流域概況

岷江發(fā)源于岷山隆起帶西側(cè)的高原腹地，主干道由北向南流經(jīng)龍門山構(gòu)造帶，與黑水河、雜谷腦河、魚溪三大主要支流交匯后，最終向東流入成都平原。根據(jù)岷江主干道與斷層的相對位置，將岷江上游主干道分為三段：岷江上游上段位于茂縣—汶川斷裂以北；岷江上游中段位于茂縣—汶川斷裂上，與斷裂帶重合；岷江上游下段位于茂縣—汶川斷裂以南（圖1）。研究區(qū)南部的岷江上游中、下段流域、魚溪流域、雜谷腦河流域受2008年汶川地震的影響廣泛發(fā)育有地震滑坡，而研究區(qū)北部黑水河及岷江上游上段流域則較少，地震滑坡對流域內(nèi)的河流地貌演化產(chǎn)生了較大的影響（Gorum et al.,2011）。因而岷江上游流域的河流地貌保留了許多構(gòu)造活動、河流演化的信息，對研究岷江上游流域地貌演化提供了豐富的資料。

3 研究方法

3.1 理論基礎(chǔ)

河流水力侵蝕模型被廣泛應(yīng)用于河流地貌研究中（Kirby et al., 2003; Whipple, 2004; Wobus et al., 2006；Kirby and Whipple, 2012），盡管基巖河流的發(fā)育環(huán)境不同，但局部河段仍可以表達為河道坡度S和上游流域面積A的冪函數(shù)關(guān)系：

式中：ks表示河道縱剖面的陡峭指數(shù)，θ表示河道縱剖面的凹度指數(shù)。通過設(shè)定θ的值（通常為0.45），可以得到標準化后的不同河流的陡峭指數(shù)ksn，進而進行不同河流的對比（Wobus et al.,2006）。這種分析方法已在不同的研究區(qū)域得到了廣泛的應(yīng)用（Korup,2006;Korup et al.,2010；張會平等，2011；Kirby and Whipple, 2012；王乃瑞等，2015；陳苗等，2018；Xue et al.,2018）。

當河道均衡狀態(tài)被破壞而處于瞬時狀態(tài)時，河道縱剖面會出現(xiàn)裂點，將河道分隔為陡峭指數(shù)不同的河段。 裂點可以分為坡斷型裂點（slope-break knickpoint）和垂階型裂點（vertical-step knickpoint），通過坡度—流域面積雙對數(shù)圖可辨別出來，破斷型裂點位于雙對數(shù)圖中數(shù)據(jù)點呈錯斷狀分布的位置，垂階型裂點則位于雙對數(shù)圖中數(shù)據(jù)點呈長釘狀分布的位置（圖2a,b）（Haviv et al.,2010； Kirby and Whipple, 2012； 王 乃 瑞 等 ，2015）。因構(gòu)造隆升或侵蝕基準面下降形成的裂點為坡斷型裂點，垂直遷移速率是恒定的，因此同一期侵蝕基準面變化所產(chǎn)生的裂點將處于同樣的海拔高度（Niemann et al.,2001;Berlin and Anderson,2007;Wei et al.,2015）。因滑坡、巖性差異等因素形成的裂點則常為垂階型裂點，在空間位置上通常是固定的，基本不會向上游遷移（Kirby et al.,2003； Crosby and Whipple, 2006； Korup, 2006;Korup et al.,2010；Berlin and Anderson,2007；王乃瑞等，2015）。通過解譯坡斷型裂點和垂階型裂點，可以幫助解譯區(qū)域的構(gòu)造活動信息和地貌演化信息且較為可靠（王乃瑞等，2015；Xue et al.,2018）。

3.2 數(shù)據(jù)采集和處理

地貌特征研究使用分辨率30 m 的ALOS DEM數(shù)據(jù)，雖然精度可能受到DEM 數(shù)據(jù)分辨率的限制，但一些研究表明通過30 m DEM數(shù)據(jù)可以獲取到河道形態(tài)及地形地貌的重要特征（Snyder et al.,2000；Duvall et al., 2004；陳苗等, 2018；Xue etal.,2018）。為研究地質(zhì)因素對河流縱剖面形態(tài)的影響，搜集了1:20 萬地質(zhì)圖，確定了主要的地層界線及斷裂的位置（圖1；圖3b）。岷江上游流域河流分支較多，支流流域面積大小不一，因而選取了流域內(nèi)5級以下共213條支流，包含岷江主干道兩側(cè)支流，及黑水河、雜谷腦河、魚溪三大支流的兩側(cè)支流，基本能覆蓋岷江上游流域內(nèi)的各個區(qū)域，便于全面了解和分析研究區(qū)河流縱剖面的形態(tài)特征及裂點的分布特征。河流的終點統(tǒng)一為岷江上游流域的出水口，即岷江主干道最終流入四川盆地的位置。

圖2 （a）坡斷型裂點的河道縱剖面與雙對數(shù)圖；（b）垂階型裂點的河道縱剖面與雙對數(shù)圖（150號、213號支流位置參考圖2）Fig.2 （a）Longitudinal profile and the double logarithmic plot of slope-drainage area of slope-break knickpoints;（b）Longitudinal profile and the double logarithmic plot of slope-drainage area of vertical-step knickpoints（Refer to Figure 2 for the position of tributary 150 and tributary 213）

根據(jù)坡度—流域面積雙對數(shù)圖，結(jié)合直觀的河流縱剖面形態(tài)，就可識別出河道中裂點的類型和空間信息。但這種提取方法在操作過程中易因河道剖面比例的不一致提取出規(guī)模不一的裂點。因此在提取過程中先將河道縱剖面置于同比例下，再結(jié)合相應(yīng)河道坡度—流域面積雙對數(shù)圖以獲得規(guī)模一致的裂點。

利用ArcGIS 軟件，首先計算研究區(qū)地形起伏度分布情況，之后提取5級以下的支流流域及河道縱剖面并利用Origin軟件將其以同比例顯示。再利用MATLAB 腳本程序提取河道縱剖面及流域面積，為進一步降低由數(shù)據(jù)產(chǎn)生的精度問題，用80 m 的移動窗口對河道進行平滑，每隔12.192 m 垂直距離計算河道坡度，繪制出坡度—流域面積雙對數(shù)圖。依據(jù)坡度—流域面積雙對數(shù)圖中數(shù)據(jù)點分布特征，配合同比例顯示的河流縱剖面，識別出規(guī)模一致的裂點及其類型， 再根據(jù)公式（logS=-θlogA+logks），設(shè)置θ=0.45 分段進行擬合得到標準化陡峭指數(shù)ksn。

4 結(jié)果與解析

4.1 裂點的平面分布特征

通過對研究區(qū)裂點分布與地形起伏度分布進行對比可以發(fā)現(xiàn)（圖2a）：研究區(qū)坡斷型裂點和垂階型裂點都有分布；起伏度較大的區(qū)域裂點分布較為密集，多為垂階型裂點，主要位于研究區(qū)南部；部分5級以下支流上游靠近分水嶺的區(qū)域則出現(xiàn)了裂點密集分布的特征。

研究區(qū)巖性分布較為復(fù)雜，但基于本文研究的需要，主要考慮巖石抗侵蝕效能的差別，根據(jù)地質(zhì)圖對研究區(qū)內(nèi)的巖性進行統(tǒng)計后我們發(fā)現(xiàn)：震旦系主要出露抗侵蝕能力較強的白云巖、硅質(zhì)巖；寒武系、奧陶系、志留系、泥盆系、三疊系主要出露抗侵蝕能力較差的千枚巖；石炭系、二疊系主要出露抗風化能力較強的灰?guī)r；抗侵蝕能力較強的花崗巖在研究區(qū)也有出露。因巖性對裂點發(fā)育的影響主要體現(xiàn)在巖體抗侵蝕效能的差異上，所以將抗侵蝕效能差異較大的地層間的裂點視為受巖性影響較大的裂點，這一類裂點通常為垂階型裂點，不會往上游遷移；因斷層形成的垂階型裂點同樣應(yīng)該出現(xiàn)在斷層線附近，不會向上游遷移（Whipple, 2004； Kirby and Whipple,2012；王乃瑞等，2015）。經(jīng)過比對（圖2b），發(fā)現(xiàn)處于地層分界線上的裂點主要是位于花崗巖與其他地層的分界線上的垂階型裂點，符合巖性形成裂點的裂點類型，如213號支流上海拔較低的垂階型裂點位于花崗巖與其他地層的分界線上（圖3b），花崗巖體相對其他地層較大的抗侵蝕能力可能是這部分垂階型裂點形成的主要原因；極個別垂階型裂點位于斷層線上，這些垂階型裂點則可能受到斷層的影響。

圖3 （a）地形起伏度與岷江上游流域裂點分布；（b）巖性分布與岷江上游流域裂點分布Fig.3 （a）Topographic relief and distribution of knickpoints in the upper reaches of Minjiang River;（b）Lithology and distribution of knickpoints in the upper reaches of Minjiang River

基于冰川的研究表明，冰川的流域面積越大，末端向下游延伸越低（Anderson et al.,2006）。高原東緣現(xiàn)代冰川的高度已退縮到較高的海拔（～4100 m），但古冰川作用形成的U 形谷則可向下延伸到3800 m 或以下的高度，在3300 m 左右仍可見到殘留的冰磧地貌（蘇珍等，2002；周尚哲等，2017）。將裂點分布與地形起伏度分布對比（圖2a；圖4a），排除巖性影響的裂點后，發(fā)現(xiàn)位于地形起伏度較低區(qū)域的裂點大部分處于高海拔區(qū)域（＞3400 m），這些高海拔裂點所在區(qū)域低地形起伏度的特點符合冰川谷地U形谷的谷地特征。而關(guān)于冰川谷地形態(tài)和冰川谷地的構(gòu)造模擬研究表明，在冰川谷上常發(fā)育冰蝕地貌和冰磧地貌形成水洼、陡坎且分布密集（MacGregor et al., 2000；Korup et al.,2010；周尚哲等，2017）。以125 號支流為例，在對靠近分水嶺的高海拔裂點進行觀察時，在河道縱剖面上也發(fā)現(xiàn)了相似的特征，表現(xiàn)為階梯狀的小陡坎，在雙對數(shù)圖上表現(xiàn)為垂階型裂點（圖4c）。而從遙感影像上也可以看到（圖4b），125 號支流上的高海拔裂點都位于冰川發(fā)育區(qū)域附近地形起伏度較低的位置（圖4a），海拔都在3400 m 以上（圖3c）。綜合以上分析，本文認為這些高海拔裂點發(fā)育的主要原因是冰川作用，而冰川谷內(nèi)形成的多個裂點，則可能受到冰蝕作用和冰磧作用的共同影響。

圖4 （a）125號支流高海拔裂點與地形起伏度分布；（b）125號支流高海拔裂點遙感影像；（c）125號支流河道縱剖面及流域面積-坡度雙對數(shù)圖像（125號支流位置參考圖2a）Fig.4 (a)High-altitude knickpoints and topographic relief of Tributary 125;(b)Remote sensing image of high altitude knickpoints of Tributary 125;(c)Longitudinal profile and the double logarithmic plot of slope-drainage area of Tributary 125（Refer to Figure 2a for the position of tributary 125）

4.2 裂點在河道縱剖面上的分布特征

岷江上游流域主要有黑水河、雜谷腦河、魚溪三大支流，發(fā)源于流域西側(cè)。在排除巖性影響的裂點和高海拔區(qū)域冰川作用影響的裂點后，本文將從河道縱剖面上裂點的分布特征來分析裂點的成因。

4.2.1 呈層狀分布的裂點特征

排除受巖性影響的裂點和高海拔（＞3400 m）區(qū)域受冰川作用影響的裂點之后，從岷江主干道及主要支流黑水河、雜谷腦河、魚溪的裂點分布特征（圖5a）可以看到，規(guī)模較大的破斷型裂點在海拔1300 m、2100 m 左右明顯成層分布。在黑水河流域我們也發(fā)現(xiàn)了相似的特征，黑水河兩側(cè)支流規(guī)模較小的坡斷型裂點在海拔1700 m、2100 m、2400 m處成層分布（圖5b）。經(jīng)過比對，發(fā)現(xiàn)這類坡斷型裂點與巖性界線、構(gòu)造線基本無對應(yīng)關(guān)系，可以視為因侵蝕基準面下降導(dǎo)致河流溯源侵蝕而形成的裂點。這種裂點的形成可能受到構(gòu)造或氣候的影響（Whipple，2001；Schoenbohm et al.,2004；Kirby and Whipple,2012）。

值得注意的是，隆升導(dǎo)致的侵蝕基準面相對下降所產(chǎn)生的裂點為坡斷型裂點，會由侵蝕基準面處形成并往上游遷移，垂直遷移速率是相似的，因此多期隆升將形成多期處于相近海拔高度的坡斷型裂點（Niemann et al., 2001; Berlin and Anderson,2007）。岷江上游流域處于青藏高原東緣強烈隆升的區(qū)域，且流域內(nèi)低海拔區(qū)域成層分布的坡斷型裂點符合這種特征。雖然氣候變化導(dǎo)致的降水變化也可能對河道縱剖面形態(tài)造成影響，但岷江上游流域相對于整個青藏高原東緣來說范圍較小，降水在研究區(qū)內(nèi)基本一致，對河流縱剖面的影響較小，因此降水變化不是這類裂點形成的主要原因。由此可見，岷江上游流域內(nèi)低海拔區(qū)域成層分布的坡斷型裂點可能主要因區(qū)域構(gòu)造隆升而形成。

4.2.2 呈階梯狀分布的裂點特征

“階梯—深潭（step-pool system）”體系一直是水利的研究中的熱點（徐江和王兆印，2004；王兆印等，2006），但在構(gòu)造地貌的研究中的相關(guān)內(nèi)容較少。這種地貌通常出現(xiàn)在山區(qū)陡坡河流中，由一段陡坡和一段緩坡加上深潭相間連接而成，呈現(xiàn)出階梯狀的特征（徐江和王兆印，2004）?；潞捅浪斐傻暮拥雷枞仔纬傻湫偷碾A梯—深潭系統(tǒng)（張康等，2011），這種特征在河道縱剖面的形態(tài)上表現(xiàn)為階梯狀的折點，在構(gòu)造地貌學中稱為裂點，其發(fā)育的間隔隨河道坡度的增大而減小，因此在坡度越陡的河道中出現(xiàn)越密集，相距距離越短（徐江和王兆印，2004）。這類因滑坡和崩塌形成的裂點通常為垂階型裂點（Korup,2006； Korup et al., 2010； Kirby and Whipple,2012）。

1933 年岷江上游疊溪一帶發(fā)生的7.5級地震使該區(qū)域發(fā)生了大規(guī)模的崩塌和滑坡，使該區(qū)域發(fā)育有大量的階梯—深潭系統(tǒng)并保留至今（四川省地震局，1983；李海龍等，2015）。堰塞湖是大型的階梯—深潭系統(tǒng)（張康等，2011），從疊溪區(qū)域的遙感影像上可以看到（圖6b），除岷江主干道上的大型堰塞湖外，松坪溝支流（即圖6中156號支流）內(nèi)還存在3個串珠狀的堰塞湖，形成了典型的階梯—深潭系統(tǒng)，與河道上的階梯狀垂階型裂點有較好的對應(yīng)關(guān)系。而岷江上游上段疊溪區(qū)域（圖6a）的水系流域發(fā)育的較多階梯狀垂階型裂點與階梯—深潭系統(tǒng)的河道縱剖面特征相符，垂階型裂點也符合滑坡和崩塌形成的裂點類型。經(jīng)過比對，這些裂點與巖性界線和構(gòu)造線基本不相關(guān)（圖2b），同時也處于低海拔區(qū)域，不受冰川作用的影響。這些特征表明地震滑坡形成的疊溪一帶規(guī)模較大的階梯—深潭系統(tǒng)在河流縱剖面上表現(xiàn)為大量的階梯狀垂階型裂點，這部分垂階型裂點的發(fā)育主要受到地震滑坡的影響。

圖5 （a）岷江主干道及主要支流河道縱剖面及裂點分布；（b）黑水河支流河道縱剖面及裂點分布；（c）96號支流河道縱剖面及雙對數(shù)圖（96號支流位置參考圖1)Fig.5 （a）Longitudinal profile and knickpoints distribution of the main channel and main tributaries of Minjiang River;（b）Longitudinal profile and distribution of knickpoints in tributaries of Heishui River;（c）Longitudinal profile and the double logarithmic plot of slope-drainage area of Tributary 96（Refer to Figure 1 for the position of tributary 96）

基于青藏高原東緣滑坡的研究表明岷江上游流域在2008 年汶川地震后廣泛發(fā)育有滑坡，滑坡主要分布在研究區(qū)南部岷江上游中段和下段流域、雜谷腦河流域和魚溪流域（Gorum et al.,2011）。研究區(qū)南部廣泛分布的垂階型裂點符合滑坡形成裂點的裂點類型（圖2a）。以汶川地震震中附近的岷江下段流域裂點分布特征（圖7a,d）為例，該區(qū)域裂點海拔較低且都為垂階型裂點，不處于冰川作用的影響區(qū)域，符合滑坡形成的裂點類型。這些裂點也基本與巖性交界和構(gòu)造線無關(guān)（圖2b）。大多數(shù)裂點都位于坡度較陡的河段，呈現(xiàn)出陡坡和緩坡相間分布的特征，在河道縱剖面上表現(xiàn)為階梯狀，符合滑坡和崩塌形成的階梯—深潭系統(tǒng)的特征。而通過區(qū)域內(nèi)典型河道8號支流2008 年汶川地震前后河道縱剖面的對比也可以看出（圖7b, c），地震后8 號支流河道縱剖面（圖7c）上出現(xiàn)了較多的階梯狀裂點?？梢娽航嫌瘟饔蚰喜窟@些階梯狀垂階型裂點的形成和發(fā)育主要受到地震滑坡的影響。

圖6 （a）疊溪區(qū)域河流縱剖面及裂點分布；（b）疊溪區(qū)域遙感影像；（c）156號支流河道縱剖面及雙對數(shù)圖（156號支流位置參考圖1)Fig.6 （a） Remote sensing image of the dammed lake by Diexi ancient earthquake;（b） Longitudinal profile and knickpoints distribution of rivers in Diexi area;（c） Longitudinal profile and the double logarithmic plot of slope-drainage area of Tributary 156（Refer to Figure 1 for the position of tributary 156）

圖7 （a）岷江上游下段河道縱剖面及裂點分布；（b）汶川地震前8號支流河道縱剖面；（c）汶川地震后8號支流河道縱剖面及雙對數(shù)圖；（d）12號支流河道縱剖面及雙對數(shù)圖（8號、12號支流位置參考圖1）Fig.7 （a） Longitudinal profile and distribution of knickpoints in the lower section of Minjiang River;（b） Longitudinal profile of Tributary 8 before Wenchuan earthquake;（c） Longitudinal profile of Tributary 8 after Wenchuan earthquake;（d）Longitudinal profile and the double logarithmic plot of slope-drainage area of Tributary 12（Refer to Figure 1 for the position of tributary 8,tributary 12）

5 結(jié)論

選擇青藏高原東緣的岷江上游流域為研究區(qū)，本文以30 m DEM數(shù)據(jù)為基礎(chǔ)，結(jié)合MATLAB腳本程序確定了岷江上游流域裂點的分布，并對裂點發(fā)育的影響因素進行了詳細分析。得出以下結(jié)論：

（1）岷江上游流域位于活動造山帶區(qū)域，流域內(nèi)水系對構(gòu)造活動和氣候變化等因素的響應(yīng)較為敏感。裂點作為河流非均衡狀態(tài)下的瞬時地貌特征，在揭示上述作用機制時有較為明顯的指示意義，可以作為研究區(qū)域地貌演化的切入點。

（2）對岷江上游流域裂點的特征分析表明，研究區(qū)水系普遍處于不均衡的狀態(tài)并廣泛發(fā)育有裂點：因巖性差異形成的裂點主要分布在花崗巖與其他地層的交界線上；因冰川作用形成的裂點位于地形起伏度較小的高海拔區(qū)域，河道縱剖面上表現(xiàn)為連續(xù)的小陡坎，反映了冰川谷地的冰蝕地貌和冰磧地貌特征；在岷江上游流域內(nèi)低海拔區(qū)域成層分布的坡斷型裂點可能主要因區(qū)域多期次的構(gòu)造隆升而形成，同期形成的坡斷型裂點最終位于相近的海拔高度；在河道縱剖面上表現(xiàn)為階梯狀的垂階型裂點則與階梯—深潭系統(tǒng)相關(guān)，主要因地震事件導(dǎo)致的大面積地震滑坡和崩塌阻塞河道而形成，在震中附近區(qū)域分布密集。研究結(jié)果也說明在對河道縱剖面分析時有必要更多考慮到滑坡及冰川等因素對現(xiàn)代河道縱剖面產(chǎn)生的重要影響。

致謝：感謝兩位審稿人對本文提出的寶貴修改意見。此外，在稿件撰寫過程中，浙江大學饒剛老師、楊蓉老師給予了幫助與指導(dǎo)；浙江大學地球科學學院提供了儀器和設(shè)備的支持。在此一并表示衷心的感謝。