段佳鑫, 董有浦, 吳 可, 汝 鑫, 張東越, 王 丹

滇中河網(wǎng)特征及其對斷層活動的響應(yīng)

段佳鑫, 董有浦*, 吳 可, 汝 鑫, 張東越, 王 丹

(昆明理工大學(xué) 國土資源工程學(xué)院, 云南 昆明 650093)

滇中地區(qū)作為青藏高原東南緣構(gòu)造活動最強(qiáng)烈的地區(qū)之一, 發(fā)育了許多活動走滑斷裂, 并且在該區(qū)域上密布著不同級別的河網(wǎng), 為利用河網(wǎng)分析斷裂活動提供了有利條件。為揭示滇中河網(wǎng)發(fā)育與斷裂活動內(nèi)在聯(lián)系, 本文基于30 m分辨率的SRTM-3數(shù)字高程模型(DEM)在ArcGIS 環(huán)境下提取了滇中地區(qū)金沙江、珠江、紅河內(nèi)的20條支流的河流網(wǎng)絡(luò), 分析其典型的河網(wǎng)參數(shù)(分支比和長度比), 通過這些比值與Horton定律描述的一般性河流網(wǎng)絡(luò)進(jìn)行對比。同時, 基于ArcGIS空間分析技術(shù)及Matlab程序腳本系統(tǒng), 對這20條河流流域縱剖面進(jìn)行定量化研究, 分析其河道基巖侵蝕力模型, 得到-擬合圖解。結(jié)果表明: 滇中地區(qū)金沙江、珠江、紅河中大多數(shù)與斷裂帶重合的高級別河流的支流, 河網(wǎng)發(fā)育明顯地受到了斷裂活動的影響; 在河道基巖侵蝕力模型-擬合圖中呈現(xiàn)為上凸和下凹相間排列, 這一現(xiàn)象在金沙江流域的各支流最為顯著, 造成這一現(xiàn)象的原因可能是研究區(qū)內(nèi)的河道主要受走滑斷裂控制, 當(dāng)斷裂以壓扭性為主時, 河道的隆升量大于剝蝕量,-擬合曲線表現(xiàn)為上凸; 當(dāng)斷裂以張扭性為主時, 河道的隆升量小于剝蝕量,-擬合曲線表現(xiàn)為下凹。通過對河網(wǎng)構(gòu)造地貌特征的研究可以使我們更好地理解該地區(qū)的構(gòu)造活動特征。

滇中; 走滑斷層; 河流地貌; Horton定律; 河道基巖侵蝕力模型; 構(gòu)造活動

0 引 言

新生代印度板塊和歐亞板塊的碰撞, 導(dǎo)致青藏高原地殼發(fā)生大規(guī)?？s短, 青藏高原東南緣的塊體沿大型走滑斷裂系向南東方向逃逸擠出(Tapponnier et al., 1982; Leloupet al., 2007)。川滇地塊沿著鮮水河?小江斷裂向東擠出(徐錫偉等, 2003; 張培震, 2008; Heet al, 2009), 在滇中地區(qū)發(fā)育了大量的南北向左行走滑斷裂(圖1b)。前人針對這些走滑斷裂開展了廣泛而又細(xì)致的構(gòu)造活動研究, 揭示鮮水河?小江斷裂、普渡河斷裂、易門斷裂、元謀?綠汁江斷裂和紅河斷裂等在晚更新世?全新世仍在活動(Replumaz et al., 2001; 向宏發(fā)等, 2006; Schoenbohm et al., 2006; 朱航和聞學(xué)澤, 2012; Wang et al., 2016; Zhang et al., 2016; 周春景等, 2016; 黃小巾等, 2018)。但前人對滇中斷裂帶的研究更多關(guān)注在幾何特征、運(yùn)動學(xué)和古地震活動性等方面。

在滇中地區(qū)主要分布金沙江、珠江、紅河三大流域, 其多條支流與斷裂帶相重合(圖2)。河流作為陸地地形地貌演化過程中最活躍的影響因子之一, 通常能對構(gòu)造變形做出很好的響應(yīng)(Fan et al., 2018)。在大型河流系統(tǒng)的發(fā)展過程中, 構(gòu)造活動可以控制河流的發(fā)育和形態(tài)(Clarket al., 2004; 文力等, 2018), 同時河流地貌也可以反映構(gòu)造活動(Kirby et al., 2003; 張會平等, 2006; 王一舟等, 2013; Valente et al., 2019)。前人利用滇中部分地區(qū)的流域地貌特征揭示了云貴高原抬升, 對斷錯水系進(jìn)行了研究(虢順民等, 1984; Schoenbohm et al., 2004; 韓竹軍等, 2016;樊云龍等, 2018), 但對于滇中地區(qū)斷裂帶的構(gòu)造活動對河網(wǎng)發(fā)育形態(tài)的影響還缺乏深入研究。滇中地區(qū)各斷裂帶上強(qiáng)烈的構(gòu)造活動對河網(wǎng)發(fā)育是否也有影響？這些強(qiáng)烈的構(gòu)造活動又使滇中地區(qū)的構(gòu)造地貌呈現(xiàn)怎樣的狀態(tài)呢？本文基于30 m分辨率的數(shù)字高程模型(digital elevation model, 簡稱DEM), 運(yùn)用Horton定律和基巖河道侵蝕力模型來梳理滇中地區(qū)(婭都紫云斷裂以西, 紅河斷裂以北及永勝斷裂以東)的河流地貌特征, 對分析滇中地區(qū)構(gòu)造活動的特征具有重要意義。

F1. 永勝斷裂; F2. 綠汁江斷裂; F3. 武定?易門斷裂; F4. 普渡河斷裂; F5. 小江斷裂; F6. 尋甸斷裂; F7. 彌勒?宗師斷裂; F8. 文山斷裂; F9. 建水?dāng)嗔? F10. 紅河斷裂; F11. 個舊斷裂; F12. 則木河斷裂。

C1. 漁鮑江; C2. 龍川江; C3. 勐果河; C4. 普渡河; C5. 小江; C6. 牛欄江; C7. 以禮河; C8. 則木河; C9. 枯木河; C10. 萬馬河; Z1. 黃泥河; Z2. 南盤江; Z3. 北盤江; Z4. 清水河; H1. 南渡河; H2. 小河底河; H3. 齋江; H4. 南滲河; H5. 勐拉河; H6. 石羊河。

1 研究區(qū)地質(zhì)背景

滇中地區(qū)平均海拔2000 m, 大地構(gòu)造位置屬于青藏高原東南緣, 其南鄰思茅地塊, 東接滇東南加里東褶皺帶, 是典型的歷史強(qiáng)震多發(fā)區(qū)(Huangfu and Qin, 2006; He et al., 2009; Wang et al., 2015)。北部的鮮水河?小江斷裂、普渡河斷裂、易門斷裂和元謀?綠汁江斷裂在新構(gòu)造以來均以左旋走滑為主(圖1、2), 且至今仍在活動(何宏林等, 1993; 徐錫偉等, 2003; Wang et al., 2008; He et al., 2009; 吳中海等, 2015; 韓竹軍等, 2016); 南部以紅河斷裂為界與印支地塊相接, 其在35 Ma左右發(fā)生過大規(guī)模左行剪切運(yùn)動(Tapponnier et al., 1982), 自5 Ma以來為右旋走滑(虢順民等, 1984; 向宏發(fā)等, 2006; Wang et al., 2012; 李西等, 2016; Zhang et al., 2019)。滇中地區(qū)基底為元古界變質(zhì)巖, 上覆有晚古生代和中生代陸相沉積(Wang et al., 2008; Wu et al., 2015), 新近紀(jì)以來, 發(fā)生顯著的斷塊差異升降運(yùn)動, 形成了典型的盆山地貌格局(Wang et al., 1998; 吳中海等, 2015)。

2 研究方法

2.1 Horton定律

Horton (1945)對美國數(shù)十個流域的河網(wǎng)發(fā)育及演變模式進(jìn)行系統(tǒng)研究后發(fā)現(xiàn), 在任何一個流域內(nèi), 各級水道的河流數(shù)量, 平均長度與河流級別之間存在一定的幾何級數(shù)變化的關(guān)系(Horton, 1945; Strahlar, 1957)。近年來, 由于數(shù)字高程模型的發(fā)展, 很多學(xué)者利用Horton定律來描述河網(wǎng)發(fā)育特征, 研究表明天然河流在不受外界因素影響時符合該定律, 且呈良好的線性關(guān)系(劉懷湘和王兆印, 2008; 劉樂等, 2015)。如果河網(wǎng)發(fā)育受到外界因素如構(gòu)造活動影響, 則會有部分河道偏離Horton定律(Kalyanjit et al., 2012)。Horton定律如下:

式中:為河流級別;N為河流數(shù)目;L為平均長度; e、1、2、1、2為常數(shù)。公式(1)為河數(shù)定律, 公式(2)為長度定律。從公式(1)中可推導(dǎo)出:

公式(3)說明, 對于Horton定律描述的河流, 不同級別的河流數(shù)目之比為常數(shù), 這個常數(shù)B稱為分支比; 不同級別的河流的平均長度之比為長度比(L), 這些比值稱為Horton比。大量研究發(fā)現(xiàn)對于不受外界因素干擾的天然河流各級別的Horton比值相近, 如B一般在3~5之間, 平均值為4;L在1.5~3之間, 平均值為2(Abrahams et al., 1984; 劉懷湘和王兆印, 2008)。

2.2 河道基巖侵蝕力模型

河道的發(fā)育和演化能夠記錄活動造山帶地區(qū)基巖的隆升以及河流的下切過程, 基巖河道侵蝕模型就是通過對河道縱剖面上的坡度變化與河流集水面積之間的相互關(guān)系來推導(dǎo)河流與基巖河床之間的應(yīng)力作用(Kirby et al., 2003; 陳彥傑等, 2006)。本文基于30 m分辨率的DEM數(shù)據(jù), 在ArcGIS 和Matlab環(huán)境下, 對河流侵蝕的結(jié)果進(jìn)行分析。

當(dāng)基巖河道的抬升速率等于河流的下切速率時, 河流發(fā)育達(dá)到均衡狀態(tài), 河道高程保持不變, 這時定義基巖河床侵蝕力模型為:

式中: d/d代表河床高程隨時間變化;為構(gòu)造抬升作用引起的地殼抬升速率;為河流下切侵蝕速率;為河道上游的流域面積;為河道坡度;為侵蝕系數(shù)(包括巖層強(qiáng)度、河床物質(zhì)、河道寬度、徑流量、土石流失率等);和為常數(shù)(代表水文、侵蝕作用等)。

則:

則公式(6)可變?yōu)?

公式(7)兩側(cè)取對數(shù)可得:

公式(8)中, 參數(shù)為均衡河道縱剖面的下凹程度指數(shù), 而K則為均衡河道縱剖面的陡度指數(shù), 可直接在-(Slope-Area)的雙對數(shù)坐標(biāo)關(guān)系圖中以直線回歸的方法得出回歸直線的斜率和截距即為與K。本文利用-雙對數(shù)曲線圖作為分析判斷地形是否達(dá)到均衡狀態(tài)的主要依據(jù)。當(dāng)?shù)匦芜_(dá)到均衡狀態(tài)時, 其集水盆地內(nèi)河流的-在雙對數(shù)圖上會呈現(xiàn)一條直線形態(tài), 當(dāng)造山帶內(nèi)部基巖隆升速率大于河流的下切侵蝕速率時, 也即河床高度隨著時間的變化而逐漸升高, 即d/d>0,-在雙對數(shù)圖上會呈現(xiàn)上凸形態(tài), 此時造山帶地形發(fā)育處于均衡前狀態(tài); d/d<0,-在雙對數(shù)圖上會呈現(xiàn)下凹形態(tài), 造山帶地形發(fā)育處于過均衡狀態(tài)(圖3)(Kirby et al., 2003; 陳彥傑等, 2006)。

3 研究結(jié)果

3.1 Horton定律結(jié)果

對DEM數(shù)據(jù)進(jìn)行預(yù)處理后, 用上述方法分別提取滇中地區(qū)金沙江、珠江和紅河流域中的20條支流(圖2)的河網(wǎng)屬性及各級河網(wǎng)屬性比。河網(wǎng)屬性在對數(shù)坐標(biāo)下與河流級別呈線性關(guān)系, 擬合公式為一條直線, 由公式(1)與公式(4)可計(jì)算Horton比中的分支比, 而河流級別在ArcGIS選擇閾值為500的匯流面積自動篩選后分級, 由于流域面積的大小不同, 金沙江的部分支流漁鮑河(C1)、小江(C5)、則木河(C8)、枯木河(C9)被自動劃分為5級, 龍川江(C2)、普渡河(C4)、以禮河(C7)自動劃分為6級, 其余的劃分為4級; 珠江支流黃泥河(Z1)、清水河(Z4)被自動劃分為5級, 南盤江(Z2)、北盤江(Z3)被自動劃分為7級和6級; 紅河的支流南渡河(H1)、勐拉河(H5)、石羊河(H6)自動劃分為5級, 小河底河(H2)、齋江(H3)、南滲河(H4)被自動分為6級(圖2)。圖4和圖5分別給出了三大流域20條河流河網(wǎng)不同級別河流數(shù)目和平均長度的關(guān)系。

3.1.1 分支比

在這20條支流中, 金沙江流域各支流的低級別河流的分支比符合Horton定律。而漁鮑河(C1)、龍川江(C2)、普渡河(C4)、小江(C5)、牛欄江(C6)、則木河(C8)、枯木河(C9)高級別河流的分支比偏離了線性關(guān)系, 且分支比并不在3~5之間, 多為2。勐果河(C3)、以禮河(C7)高級別河流的分支比分別為7和8。而萬馬河(C10)低級別河流與高級別河流的分支比都符合Horton定律。其中, C1~C9高級別河流與斷裂帶相互重合(圖2), 如: 支流漁鮑河(C1)、枯木河(C9)高級別河流與永勝斷裂南段相互重合, 龍川江(C2)與綠汁江斷裂相互重合, 勐果河(C3)與武定?易門斷裂相互重合, 普渡河(C4)與普渡河斷裂相互重合, 小江(C5)與小江斷裂相互重合, 牛欄江(C6)與魯?shù)閿嗔严嗷ブ睾? 以禮河(C7)與會澤斷裂相互重合, 則木河(C8)與則木河斷裂相互重合。

珠江流域各支流的低級別河流的分支比符合Horton定律, 但南盤江(Z2)和清水河(Z4)2條高級別河流的分支比偏離了線性關(guān)系, 且不在3~5之間, 其多為2或8。南盤江(Z2)最高級別河流與小江斷裂南段、彌勒斷裂相互重合, 清水河(Z4)與丘北斷裂相互重合。

紅河流域各支流的低級別河流分支比均符合Horton定律, 南渡河(H1)和石羊河(H6)2條高級別河流的分支比偏離線性關(guān)系, 且它們的分支比并不在3~5以內(nèi), 也多為2, 其中南渡河(H1)與綠汁江斷裂重合, 石羊河(H6)與紅河斷裂北段相重合(圖4, 表1)。

圖3 河流水力侵蝕模型S-A雙對數(shù)圖解(據(jù)陳彥傑等, 2006)

圖4 滇中各支流不同級別河流數(shù)量N與級別w的關(guān)系(河流名稱見圖2)

3.1.2 長度比

金沙江流域各支流中低級河流的長度比都符合Horton定律, 但C1~C9高級別河流偏離線性關(guān)系, 其長度比多為1.5以下。珠江流域的北盤江(Z3)高級別河流長度比偏離了Horton定律, 在1.5以下。紅河流域支流中H1~H6幾條低級別河流符合Horton定律, 而勐拉河(H5)和石羊河(H6)高級別河流偏離線性關(guān)系, 且長度比并不在1.5~3之間, 高級別河流的長度比多為3以上或在1.5以下(圖5, 表1)。

3.2 河道基巖侵蝕力模型

對滇中地區(qū)三大流域內(nèi)的20條支流的基巖河道進(jìn)行統(tǒng)計(jì), 建立log和log的關(guān)系圖(圖6~8)。結(jié)果發(fā)現(xiàn):-關(guān)系圖呈現(xiàn)明顯分段性, 表現(xiàn)為既有上凸也有下凹, 部分呈直線形態(tài)。各支流的基巖河道侵蝕力模型結(jié)果如下:

3.2.1 金沙江流域

金沙江支流中漁鮑河(C1)、勐果河(C3)和以禮河(C7)的-擬合曲線呈現(xiàn)上游為直線形態(tài), 到中下游后呈上凸、下凹相間排列, 其河道中既有均衡地形, 也存在隆升和沉降地形; 龍川江(C2)、普渡河(C4)、小江(C5)、牛欄江(C6)、以禮河(C7)、則木河(C8)和枯木河(C9)幾條支流從上游開始就呈明顯的上凸、下凹相間排列的形態(tài), 幾條支流既有隆升地形也有沉降地形; 萬馬河(C10)呈現(xiàn)直線形態(tài), 為均衡地形。

3.2.2 珠江流域

珠江支流中江底大河(Z1)和清水河(Z4)的-擬合曲線上游為直線形態(tài), 到中下游后呈上凸、下凹相間排列, 河道中既有均衡狀態(tài), 也存在隆升地形和沉降地形; 南盤江(Z2)和北盤江(Z3)支流從上游開始就呈明顯的上凸、下凹相間排列的形態(tài), 既有隆升地形也有沉降地形。

3.2.3 紅河流域

紅河支流中石羊河(H6)的-擬合曲線上游為直線形態(tài), 到中下游后呈上凸、下凹相間排列, 河道中既有均衡狀態(tài)也存在隆升地形和沉降地形; 南渡河(H1)、小河底河(H2)、齋河(H3)和石羊河(H6)支流從上游開始就呈明顯的上凸、下凹相間排列的形態(tài), 既有隆升地形也有沉降地形; 南滲河(H4)、勐拉河(H5)呈現(xiàn)直線形態(tài), 為均衡地形。

4 討 論

4.1 河流發(fā)育是否受構(gòu)造活動的影響

滇中大多數(shù)支流的河網(wǎng)中低級別河流的分支比和長度比都呈良好的線性關(guān)系, 符合Horton定律, 但隨著級別的增加, 分支比和長度比都逐漸發(fā)生了偏離, 且河流級別越高, 其偏離得越多。Kalyanjit et al. (2012)發(fā)現(xiàn)水系分支比對構(gòu)造活動的響應(yīng)較為敏感, 可以反映出區(qū)域的構(gòu)造活動強(qiáng)弱: 分支比偏離程度越高, 構(gòu)造活動性越強(qiáng); 分支比偏離程度低, 構(gòu)造活動性則越弱。

金沙江支流中的漁鮑河(C1)、龍川江(C2)、勐果河(C3)、普渡河(C4)、小江(C5)、牛欄江(C6)、以禮河(C7)、則木河(C8)、枯木河(C9)高級別河流分別與程海?賓川斷裂東側(cè)、綠汁江斷裂、武定?易門斷裂、普渡河斷裂、小江、魯?shù)閿嗔选蓴嗔?、則木河斷裂、程海?賓川斷裂南段相互重合(圖2), 這些斷裂構(gòu)造活動強(qiáng)烈(徐錫瑋等, 2003; 王閻昭等, 2008; Wu and Hu, 2019), 可能導(dǎo)致河網(wǎng)數(shù)目與長度發(fā)生改變。而萬馬河(C10)不與斷裂帶相互重合, 受構(gòu)造活動影響小, 分支比、長度比則均不偏離Horton定律。

珠江支流中南盤江(Z2)和北盤江(Z3)高級別河流分別與彌勒斷裂、楊家大坡斷裂相互重合。彌勒斷裂在第四紀(jì)以來屬于弱構(gòu)造活動斷層(Wu et al., 2015; Wu and Hu, 2019), 但南盤江(Z2)的第4級河流與小江斷裂南段、曲江斷裂相互重合, 小江斷裂與曲江斷裂在第四紀(jì)來構(gòu)造活動強(qiáng)烈(聞學(xué)澤等, 2011; Wang et al., 2014), 因此南盤江(Z2)的河網(wǎng)發(fā)育受構(gòu)造活動影響。

圖5 滇中各支流水系不同級別河流平均長度L(m)與級別w的關(guān)系(河流名稱見圖2)

表1 滇中各支流不同級別河流的分支比RB與長度比RL

紅河支流中南渡河(H1)、小河底河(H2)、南滲河(H4)、石羊河(H6)與綠汁江斷裂、建水?dāng)嗔?、個舊斷裂、紅河斷裂北段相互重合。由于這些斷裂帶構(gòu)造活動強(qiáng)烈, 與其相互重合的支流的分支比、長度比偏離Horton定律。其中小河底河(H2)最高級別長度比較大, 而分支比卻較小, 可能是由于流域受到斷裂帶擠壓變得狹長。綜上, 不與斷裂帶重合或與構(gòu)造活動弱重合的河網(wǎng)基本上遵循Horton的指數(shù)變化規(guī)律, 而與強(qiáng)活動斷裂帶相互重合的支流中高級別河流的分支比和長度比大多偏離Horton定律。

4.2 走滑斷層對地貌發(fā)育的影響

滇中地區(qū)金沙江、珠江、紅河三大流域的-雙對數(shù)曲線呈現(xiàn)上凸、下凹相間排列的特征, 反映滇中地區(qū)既有地形隆升又有沉降。而前人在大型擠壓性斷裂控制下的天山北麓河、白龍江、龍門山等地區(qū)的研究結(jié)果(趙洪壯等, 2009; 蘇琦等, 2016; 閆亮等, 2018)與滇中地區(qū)的河道基巖侵蝕力模型存在偏差, 可能是由于滇中地區(qū)主要受大型走滑斷裂控制(Wang et al., 1998), 大型走滑斷層系中常伴隨著斷塊差異性升降運(yùn)動, 出現(xiàn)拉張性斷陷盆地和擠壓性逆沖?褶皺拱起共同存在的情況(Burchfiel et al., 1987)。

金沙江流域主要受綠汁江斷裂、武定?易門斷裂、普渡河斷裂、小江斷裂等強(qiáng)地震活動性的SN向走滑斷裂控制(Wang et al., 1998; He et al., 2009)。前人研究發(fā)現(xiàn)這些走滑斷裂發(fā)生張扭性運(yùn)動, 形成了滇池、玉溪、元謀等第四紀(jì)斷陷盆地(Wang et al., 1998; He et al., 2009; 裴向軍等, 2019), 同時也發(fā)育了擠壓或逆沖斷層(Ren, 2013; Wang et al., 2014; 樊云龍等, 2018)。

珠江流域主要受小江斷裂、南盤江斷裂、北盤江斷裂和曲靖斷裂等走滑斷裂活動的影響。形成了彌勒盆地(Li et al., 2015)和北盤江流域的差異隆升(Wang et al., 1998; 樊云龍等, 2018)。

紅河流域受建水?dāng)嗔选⑶瓟嗔?、紅河斷裂北段、綠汁江斷裂和個舊斷裂等影響, 形成了許多第四紀(jì)斷陷盆地(Wang et al., 2012), 而在曲江?石屏斷裂帶剖面上又表現(xiàn)為自北向南推覆的逆沖斷裂帶結(jié)構(gòu)(圖9a)(韓新民等, 1982; 韓新民和毛玉平, 1993; 樊云龍等, 2018)。

以小江為例, 小江流域主要是受小江斷裂活動控制。小江斷裂由多條左行走滑斷裂組成, 斷裂以西的湯丹沖斷塊及周邊地區(qū)在第四紀(jì)發(fā)生隆升, 小江中段地區(qū)地殼以整體抬升為主(Wang et al., 1998; 王剛和王二七, 2005)。另外, 小江斷裂還控制了許多伸展盆地(宋方敏等, 1997)。小江流域的上凸區(qū)與小江中段的逆沖斷裂以及湯丹沖斷塊相對應(yīng); 下凹區(qū)則與嵩明盆地、東川盆地相對應(yīng)(圖9b、c)。當(dāng)走滑斷裂表現(xiàn)為壓扭為主時, 可能造成隆升量大于剝蝕量,-雙對數(shù)曲線呈現(xiàn)為上凸; 當(dāng)走滑斷裂表現(xiàn)為伸展為主時, 可能造成隆升量小于侵蝕量,-雙對數(shù)曲線呈現(xiàn)為下凹。

圖6 金沙江流域內(nèi)各支流水力侵蝕模型S-A雙對數(shù)圖

圖7 珠江流域內(nèi)各支流水力侵蝕模型S-A雙對數(shù)圖

圖8 紅河流域內(nèi)各支流水力侵蝕模型S-A雙對數(shù)圖

(a) 滇中地區(qū)構(gòu)造活動簡圖; (b) 小江流域構(gòu)造活動簡圖; (c) 小江流域S-A雙對數(shù)曲線。

5 結(jié) 論

通過對滇中地區(qū)的20條干流流域縱剖面進(jìn)行定量化研究, 獲得以下認(rèn)識:

(1) 滇中地區(qū)低級別河流屬性較好地符合Horton定律, 而高級別河流屬性與Horton定律有較大偏離, 且偏離Horton定律的高級別河流與走滑斷裂帶相互重合, 滇中地區(qū)河網(wǎng)發(fā)育明顯受到斷裂活動的影響。

(2) 滇中地區(qū)河流網(wǎng)的-擬合圖呈現(xiàn)出明顯分段性, 在同一條河流內(nèi)表現(xiàn)為上凸區(qū)和下凹區(qū)相間排列。其原因可能是研究區(qū)內(nèi)的河道主要受走滑斷裂控制, 當(dāng)斷裂以壓扭性為主時, 河道的隆升量大于剝蝕量, 在-擬合曲線表現(xiàn)為上凸; 當(dāng)斷裂為張扭性為主時, 河道的隆升量小于剝蝕量, 在-擬合曲線表現(xiàn)為下凹。通過對河網(wǎng)地貌特征的研究可以使我們更好地理解該地區(qū)的構(gòu)造活動特征。

致謝:衷心感謝兩位匿名審稿專家對本文提出的建設(shè)性意見！

陳彥傑, 宋國城, 陳昭男. 2006. 非均衡山脈的河流水力侵蝕模型. 科學(xué)通報, 51(7): 865–870.

樊云龍, 潘保田, 胡振波, 任大銀, 陳起偉, 劉芬良, 李宗盟. 2018. 云貴高原北盤江流域構(gòu)造地貌特征分析. 地球科學(xué)進(jìn)展, 33(7): 751–761.

虢順民, 張靖, 李祥根, 向宏發(fā), 陳鐵牛, 張國偉. 1984. 云南紅河斷裂帶北段斷裂位錯與地震重復(fù)發(fā)生的時間間隔. 地震地質(zhì), 6(1): 1–12.

韓新民, 柴天俊, 肖九安, 李盛德. 1982. 石屏?建水?dāng)嗔阎卸蔚男禄顒优c地震. 地震研究, 5(2): 220–225.

韓新民, 毛玉平. 1993. 石屏–建水?dāng)嗔褞磥砣陜?nèi)七級以上大地震危險性分析. 地震研究, 16(1): 52–59.

韓竹軍, 董紹鵬, 毛澤斌, 呼楠, 譚錫斌, 袁仁茂, 郭鵬. 2016. 小江斷裂帶南段全新世活動的地質(zhì)地貌證據(jù)與滑動速率. 地震地質(zhì), 39(1): 2–18.

何宏林, 方仲景, 李坪. 1993. 小江斷裂帶西支斷裂南段新活動初探. 地震研究, 16(3): 291–298.

黃小巾, 吳中海, 黃小龍, 羅睿潔. 2018. 滇西北程海?賓川斷裂帶第四紀(jì)分段活動性的構(gòu)造地貌表現(xiàn)與限定. 地球科學(xué), 43(12): 4651–4670.

李西, 冉勇康, 陳立春, 王虎, 于江, 張彥琪, 謝英情. 2016. 紅河斷裂帶南段全新世地震活動證據(jù). 地震地質(zhì), 38(3): 598–604.

劉懷湘, 王兆印. 2008. 河網(wǎng)形態(tài)與環(huán)境條件的關(guān)系. 清華大學(xué)學(xué)報(自然科學(xué)版), 48(9): 1408–1412.

劉樂, 王兆印, 余國安, 栗騰. 2015. 青藏高原河網(wǎng)統(tǒng)計(jì)規(guī)律及高原抬升的影響. 清華大學(xué)學(xué)報(自然科學(xué)版), 55(9): 964–970.

裴向軍, 李天濤, 黃潤秋, 王雙. 2019. 巧家拉分盆地結(jié)構(gòu)特征及其形成演化過程分析. 西南交通大學(xué)學(xué)報, 54(2): 278–286.

蘇琦, 梁明劍, 袁道陽, 謝虹, 吳趙. 2016. 白龍江流域構(gòu)造地貌特征及其對滑坡泥石流災(zāi)害的控制作用. 地球科學(xué), 41(10): 1758–1770.

宋方敏, 汪一鵬, 俞維賢. 1998. 中國活斷層研究專輯//小江活動斷裂帶. 北京: 地震出版社.

王剛, 王二七. 2005. 擠壓造山帶中的伸展構(gòu)造及其成因——以滇中地區(qū)晚新生代構(gòu)造為例. 地震地質(zhì), 27(2): 188–199.

王閻昭, 王恩寧, 沈正康, 王敏, 甘衛(wèi)軍, 喬學(xué)軍, 孟國杰, 李鐵明, 陶瑋, 楊永林, 程佳, 李鵬. 2008. 基于GPS資料約束反演川滇地區(qū)主要斷裂現(xiàn)今活動速率. 中國科學(xué)(D輯), 38(5): 582–597.

王一舟, 張會平, 俞晶星, 張竹琪, 鄭文俊, 鄭德文. 2013. 祁連山洪水壩河流域地貌特征及其構(gòu)造指示意義. 第四紀(jì)研究, 33(4): 737–745.

文力, 魏鵬飛, 常華進(jìn), 孫小舟, 李權(quán)國, 杜紅旺, 劉琛, 李學(xué)敏, 秦玉莉. 2018. 青藏高原周邊地區(qū)河流分形特征與地貌、構(gòu)造活動耦合關(guān)系. 地質(zhì)通報, 37(6): 965–974.

聞學(xué)澤, 杜方, 龍鋒, 范軍, 朱航. 2011. 小江和曲江?石屏兩斷裂系統(tǒng)的構(gòu)造動力學(xué)與強(qiáng)震序列的關(guān)聯(lián)性. 中國科學(xué): 地球科學(xué), 41(5): 713–724.

吳中海, 龍長興, 范桃園, 周春景, 馮卉, 楊振宇. 2015. 青藏高原東南緣弧形旋扭活動構(gòu)造體系及其動力學(xué)特征與機(jī)制. 地質(zhì)通報, 34(1): 1–31.

向宏發(fā), 萬景林, 韓竹軍, 虢順民, 張晚霞, 陳立春, 董興權(quán). 2006. 紅河斷裂帶大型右旋走滑運(yùn)動發(fā)生時代的地質(zhì)分析與FT測年. 中國科學(xué)(D輯), 36(11): 977– 987.

徐錫偉, 聞學(xué)澤, 鄭榮章, 馬文濤, 宋方敏, 于貴華. 2003. 川滇地區(qū)活動塊體最新構(gòu)造變動樣式及其動力來源. 中國科學(xué)(D輯), 33(S1): 151–162.

閆亮, 李勇, 趙國華, 梁明劍, 顏照坤, 邵崇建, 李奮生, 云錕, 馬超, 張威. 2018. 青藏高原東緣龍門山構(gòu)造帶晚第四紀(jì)構(gòu)造隆升作用的河流地貌響應(yīng). 第四紀(jì)研究, 38(1): 232–246.

張會平, 楊農(nóng), 張?jiān)罉? 孟暉. 2006. 岷江水系流域地貌特征及其構(gòu)造指示意義. 第四紀(jì)研究, 26(1): 126–135.

張培震. 2008. 青藏高原東緣川西地區(qū)的現(xiàn)今構(gòu)造變形、應(yīng)變分配與深部動力過程. 中國科學(xué)(D輯), 38(9): 1041–1056.

趙洪壯, 李有利, 楊景春, 呂紅華. 2009. 天山北麓河流縱剖面與基巖侵蝕模型特征分析. 地理學(xué)報, 65(5): 563–570.

周春景, 吳中海, 馬曉雪, 李家存, 王繼龍. 2016. 滇中城市群重要活動斷裂與區(qū)域地殼穩(wěn)定性評價. 地質(zhì)力學(xué)學(xué)報, 22(3): 454–461, 463–477.

朱航, 聞學(xué)澤. 2012. 小江?則木河斷裂帶大地震序列的靜應(yīng)力觸發(fā)作用. 地球科學(xué), 37(1): 199–206.

Abrahams A D. 1984. Channel networks: A geomorphologicalperspective., 20(2): 161–188.

Burchfiel B C, Hodges K V and Royden L H. 1987. Geology of panamint valley-saline valley pull-apart system, California; Palinspastic evidence for low-angle geometry of a Neogene range-bounding fault., 92: 10422–10426.

Clark M K, Schoenbohm L M and Royden L H, Whipple K X, Burchfiel B C, Zhang X, Tang W, Wang E and Chen L. 2004. Surface uplift, tectonics, and erosion of eastern Tibet from large-scale drainage patterns., 23(1), TC1006.

Fan N N, Chu Z X, Jiang L G, Hassan M A, Lamb M P and Liu X N. 2018. Abrupt drainage basin reorganization following a Pleistocene river capture., 9(1), 3756.

He J K, Lu S J and Wang X G. 2009. Mechanical relation between crustal rheology, effective fault friction, and strike-slip partitioning among the Xiaojiang fault system, southeastern Tibet., 34(3): 363–375.

Horton R E. 1945. Erosional development of streams and their drainage basins: Hydrophysical applications of quantitative morphology., 56(3): 275–370.

Huangfu G and Qin J Z. 2006. Study of the seismicity of strong earthquakes in Yunnan area., 20(4): 449–459.

Kalyanjit S, Jha L K and Tiwari B K. 2012. Morphometric analysis of a highland microwatershed in eastkhasi hills district of Meghalays, India: Using remote sensing and geographic informationsystem (GIS) techniques., 5(5): 142–150.

Kirby E, Whipple K X, Tang W Q and Chen Z L. 2003. Distribution of active rock uplift along the eastern margin of the Tibetan Plateau: Inferences from Bedrock Channel Longitudinal profiles., 108(B4): 148–227.

Leloup P H, Tapponnier P, Lacassin R and Searle M P. 2007. Discussion on the role of the Red River shear zone, Yunnan and Vietnam, in the continental extrusion of SE Asia., 164(6): 1253– 1260.

Li S H, Deng C L, Dong W, Sun L, Liu S Z, Qin H F, Yin J Y, Ji X P and Zhu R X. 2015. Magnetostratigraphy of the Xiaolongtan Formation Bearing Lufengpithecus Keiyuanensis in Yunnan, Southwestern China: Constraint on the initiation Time of the Southern Segment of the Xianshuihe-Xiaojiang Fault., 655(1): 213–226.

Ren Z K. 2013. Geometry and deformation features of the most recent co-seismic surface ruptures along the Xiaojiang fault and its tectonic implications for the Tibetan Plateau., 77: 21–30.

Replumaz A, Lacassin R, Tapponnier P and Leloup P H. 2001. Large river offsets and Plio-Quaternary dextral slip rate on the Red River fault (Yunnan, China)., 106(B1): 819–836.

Schoenbohm L M, Burchfiel B C and Chen L. 2006. Propagation of surface uplift, lower crustal flow, and Cenozoic tectonics of the southeastern margin of the Tibetan Plateau., 34(10): 813–816.

Schoenbohm L M, Whipple K X, Burchfiel B C and Chen L. 2004. Geomorphic constraints on surface uplift, exhumation, and plateau growth in the Red River region, Yunnan Province, China., 116(7–8): 895–909.

Strahlar A N. 1957. Quantitative analysis of watershed gommphology., 38(6): 913–920.

Tapponnier P, Peltzer G, Le Dain A Y and Armijo R. 1982. Propagating extrusion tectonics in Asia: New insights from simple experiments with plasticine., 10(12): 611–616.

Valente E, Buscher J T, Jourdan F, Petrosino P, Reddy S M, Tavani S, Corradetti A and Ascione A. 2019. Constraining mountain front tectonic activity in extensional setting from geomorphology and Quaternary Stratigraphy: A case study from the Matese ridge, Southern Apennines., 219: 47–67.

Wang E C, Burchfiel B C, Royden L H, Chen L Z, Chen J S, Li W X and Chen Z L. 1998. Late Cenozoic Xianshuihe- Xiaojiang, Red River, and Dali Fault systems of southwestern Sichuan and central Yunnan, China., 327: 1–108.

Wang F, Wang M, Wang Y, Wang Y Z and Shen Z K. 2015. Earthquake potential of the Sichuan-Yunnan region, western China., 107: 232–243.

Wang S F, Jiang G G, Xu T D and Zheng D W. 2012. The Jinhe-Qinghe Fault an inactive branch of the Xianshuihe- Xiaojiang Fault Zone, eastern Tibet., 544–545: 93–102.

Wang Y, Zhang B, Hou J, Hou J J and Xu X W. 2014. Structure and tectonic geomorphology of the Qujiang Fault at the intersection of the Ailao Shan-Red River Fault and the Xianshuihe-Xiaojiang Fault system, China., 634: 156–170.

Wang Y, Zhang B, Schoenbohm L M, Zhang J J, Zhou R J, Hou J J and Ai S. 2016. Late Cenozoic tectonic evolution of the Ailao Shan-Red River fault (SE Tibet): Implications for kinematic change during plateau growth., 35(8): 1969–1988.

Wang Y Z, Wang E N, Shen Z K, Shen Z K, Wang M, Gan W J, Qiao X J, Meng G J, Li T M, Tao W, Yang Y L, Cheng J and Li P. 2008. GPS-constrained inversion of present-day slip rates along major faults of the Sichuan- Yunnan region, China.(), 51(9): 1267–1283.

Wu G J, Tan H B, Yang G L and Shen C Y. 2015. Research on the relationship between geophysical structural features and earthquakes in Mid-Yunnan and the surrounding area., 6(5): 384–391.

Wu Z H and Hu M M. 2019. Neotectonics, active tectonics and earthquake geology: Terminology, applications and advances., 127: 1–15.

Zhang Y, He Z T, Ma B Q, Wang J Y, Zhang H and Wang J P. 2019. Geological and geomorphic evidence for activity of the Mengzi fault along the southeastern margin of the Tibetan plateau., 171: 233–245.

Zhang Y S, Yao X, Yu K, Du G L and Guo C B. 2016. Late-Quaternary slip rate and seismic activity of the Xianshuihe Fault zone in Southwest China., 90(2): 525–536.

Characteristics of River Networks in the Central Yunnan Province and Their Responses to Tectonic Activities

DUAN Jiaxin, DONG Youpu*, WU Ke, RU Xin, ZHANG Dongyue and WANG Dan

(Faculty of Land and Resource Engineering,, Kunming 650093, Yunnan, China)

The central Yunnan province is located in the Southeast of the Tibetan Plateau. Being one of the most tectonically active regions, it is cut by several active strike-slip faults. It is an ideal place to investigate the relationship between the active tectonics and river systems. Here, the SRTM-DEM with spatial resolution of 30 m was used as the basic dataset, and the ArcGIS and MATLAB software were used to the extract profile characteristics of 20 tributaries from the Jinsha River, the Pearl River, and the Red River in the central Yunnan province. First, the networks of the 20 tributaries were mapped and analyzed. Secondly, the relationship between tectonic uplift and fluvial incision was analyzed by steep index (). The results show that (1) the networks are represented by the Horton Laws for the streams of the lower orders but deviate dramatically at the highest order. (2)-plots all exhibit alternate convex and concave forms. It can be inferred that the developments of high-order rivers were significantly affected by the strike-slips faults. The analysis of fluvial geomorphology can help us to further investigate the tectonic activities and characteristics of the central Yunnan province.

the central Yunnan province; strike-slip fault; fluvial geomorphology; Horton law; bedrock channel incision model; tectonic activity

2019-08-29;

2020-01-17

國家自然科學(xué)基金項(xiàng)目(41802215、41762017、41672206)與云南省教育廳科學(xué)研究基金項(xiàng)目(KKJB201721011)聯(lián)合資助。

段佳鑫(1994–), 女, 碩士研究生, 構(gòu)造地質(zhì)學(xué)專業(yè)。Email: 729234235@qq.com

董有浦(1983–), 男, 副教授, 主要從事構(gòu)造地質(zhì)學(xué)等相關(guān)領(lǐng)域的研究。Email: dongypsd@126.com

P542

A

1001-1552(2021)02-0296-012

10.16539/j.ddgzyckx.2021.02.003