朱曉潔, 褚忠信*, 于廣科,毛士博, 姜瑾斐, 朱龍海

(1. 中國海洋大學(xué) 海洋地球科學(xué)學(xué)院, 山東 青島 266100; 2. 海底科學(xué)與探測(cè)技術(shù)教育部重點(diǎn)實(shí)驗(yàn)室, 山東 青島 266100)

山區(qū)河流的襲奪演化和分水嶺的遷移對(duì)流域地貌演變發(fā)育過程起著重要的作用,同時(shí),河流水系演化也是地球表層循環(huán)系統(tǒng)中的重要組成部分[1].分水嶺遷移的結(jié)果是在一定的時(shí)間內(nèi)使河道發(fā)生改道[2],進(jìn)而引起河流發(fā)生襲奪.分水嶺遷移通常多發(fā)生在不對(duì)稱分水嶺區(qū)域,隨著時(shí)間的推移,下切侵蝕速率快、位置較低的一側(cè)河流率先切入并侵蝕破壞分水嶺,使分水線不斷地向侵蝕能力較弱、位置較高的一側(cè)河流移動(dòng),最終導(dǎo)致原先完整的分水嶺斷開并發(fā)生遷移[3-4],因此對(duì)斷開的分水線進(jìn)行連接恢復(fù)是重建襲奪前水系的關(guān)鍵環(huán)節(jié).對(duì)于河流襲奪前水系的重建,已有學(xué)者根據(jù)最陡路徑法[5]、流域面積插值法[6]等對(duì)其進(jìn)行重建研究.

距今約8萬年前,在柴汶河上游流域發(fā)育一起較為年輕的河流襲奪事件[6-9],其襲奪地貌特征保存完整，兼具充足的地貌學(xué)及沉積學(xué)證據(jù),在全球范圍內(nèi)具有稀缺性與典型性.由于數(shù)字高程模型(DEM)囊括了豐富的地形、地貌及水文等空間分布信息[7],因此基于DEM數(shù)據(jù)對(duì)流域地形地貌特征進(jìn)行提取已成為獲取流域信息的主要手段[8].本文對(duì)柴汶河河流襲奪前的水系進(jìn)行重建分析,以期有助于進(jìn)一步了解河流襲奪理論、認(rèn)識(shí)水系變遷規(guī)律、豐富流域地貌認(rèn)知及規(guī)劃河流的開發(fā)和治理.

1 研究區(qū)域概況

此研究區(qū)域位于山東省淄博市沂源縣大張莊鎮(zhèn)西南邊緣,117°56′E～118°2′E,35°56′N～36°3′N之間,地處魯中腹地,西北毗鄰濟(jì)南市萊蕪區(qū),西南與新泰市、蒙陰縣相連,東北地處沂源縣境內(nèi),屬4縣(市)交界位置[9].晚更新世時(shí)期,海拔較低的黃河支流——古柴汶河向上溯源侵蝕襲奪了海拔較高的淮河支流——古沂河[6],使古沂河流向由東北急劇向西南作了近90°的大拐彎,在形態(tài)上呈倒勾狀彎曲,因此古沂河上游匯入古柴汶河,形成現(xiàn)如今的柴汶河流域,如圖1所示.

圖1 柴汶河流域位置

柴汶河流域[10]源于沂源縣西南部牛欄峪一帶,東流至周科峪后西折向西南流,穿新泰市東周水庫,經(jīng)新泰、寧陽縣境至大汶口入主流牟汶河,其干支流都是源短流急的山洪河流,洪水漲落迅猛,平時(shí)只有涓涓細(xì)流.流域內(nèi)地勢(shì)東高西低,北高南低,北有泰山,東靠魯山和蒙山,西、南為丘陵和平原.柴汶河全長(zhǎng)116 km,流域面積1 944 km2,是汶河的最大支流.

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

2.1 數(shù)據(jù)來源本文采用的實(shí)驗(yàn)數(shù)據(jù)由地理空間數(shù)據(jù)云提供(http://www.gscloud.cn/),為ASTER GDEM-DEM數(shù)據(jù),空間分辨率為30 m×30 m,空間參考為WGS_1984_UTM_Zone_50N.對(duì)所獲得的DEM數(shù)據(jù)利用ArcGIS 10.2軟件進(jìn)行拼接裁剪處理,得到河流襲奪研究區(qū)的原始DEM數(shù)據(jù),面積約為80 km2.同時(shí),在原始DEM數(shù)據(jù)上疊加部分水系流域，得到柴汶河上游流域水系分布,如圖2所示.

圖2 柴汶河上游流域附近水系分布

2.2 最佳集水面積閾值確定將河流襲奪研究區(qū)的原始DEM數(shù)據(jù)加載到ArcGIS軟件中,對(duì)原始的DEM進(jìn)行填洼[11]、流向提取、匯流累積量計(jì)算[12]以及河網(wǎng)水系的提取[13-15],再通過均值變點(diǎn)分析法[16]確定與實(shí)際水系接近的河網(wǎng)水系[17].均值變點(diǎn)分析法目前主要應(yīng)用于提取地形起伏度的最佳單元[18-21],同時(shí)也逐漸應(yīng)用于確定河流的最佳集水面積閾值[22].

2.3 分水線提取與恢復(fù)分水線通常為分水嶺最高點(diǎn)之間的連線,同時(shí)也為水流的起源點(diǎn)[23].因此，對(duì)DEM數(shù)據(jù)進(jìn)行地表徑流模擬計(jì)算時(shí),通過提取匯流累積量的零值可得到分水線[24].

斷開的分水線在DEM數(shù)據(jù)上主要表現(xiàn)為高程值由高變低,因此將分水線的高程值作為恢復(fù)襲奪前分水嶺的單一因子,即通過抬高斷開的分水線的高程值恢復(fù)襲奪前的分水嶺.已知斷開分水線兩端點(diǎn)的坐標(biāo)及高程值(x1,y1,z1)、(x2,y2,z2),利用C++編寫程序?qū)嚅_的分水線之間的柵格依次重新賦予高程值,其原理如圖3所示.

圖3 分水線高程值重新賦值原理

首先,基于均值替換法[25],計(jì)算被連接的分水線兩端點(diǎn)(A、B)高程的平均值(即中點(diǎn)值)代替襲奪前的原中點(diǎn)(C)高程值;再基于二次曲線內(nèi)插法相似原理,設(shè)Zi為被連接的分水線對(duì)應(yīng)的各柵格單元高程值,可得到某一柵格單元(D、E)高程值的計(jì)算公式:

Z

(1)

Z

(2)

3 結(jié)果與分析

3.1 河網(wǎng)水系提取基于ArcGIS水文分析模塊,分別設(shè)置匯流累積量閾值自小到大為100、200、300、…、2 600共26組數(shù)據(jù)進(jìn)行水系的提取,以100、600和1 600匯流累積量閾值為例,僅提取柴汶河上游流域河網(wǎng)水系,如圖4(a)～(c)所示.根據(jù)結(jié)果可知,隨著匯流累積量閾值的增大,提取出來的河網(wǎng)水系逐漸稀疏,河流總長(zhǎng)度逐漸變短,水系分級(jí)逐漸不詳細(xì),在后續(xù)研究分析時(shí)會(huì)對(duì)結(jié)果的準(zhǔn)確度產(chǎn)生影響[26].因此,為提高結(jié)果準(zhǔn)確度,采用均值變點(diǎn)分析法確定當(dāng)集水面積閾值為0.54 km2時(shí)(即匯流累積量閾值為600時(shí)),根據(jù)提取出來的水系與1∶5 000的實(shí)際河網(wǎng)水系(如圖4(d)所示)進(jìn)行對(duì)比發(fā)現(xiàn)較為吻合,最接近于實(shí)際河網(wǎng)水系.

3.2 分水線恢復(fù)河流襲奪與分水嶺遷移是相輔相成的2個(gè)過程,河流襲奪必然引起分水嶺的遷移;同樣,分水嶺遷移會(huì)引起水系重組,進(jìn)而發(fā)生河流襲奪.分水嶺是分隔2個(gè)流域的界線,是地球內(nèi)外營力共同作用的產(chǎn)物.當(dāng)分水嶺兩側(cè)的水系侵蝕能力不同時(shí),分水嶺會(huì)向侵蝕能力弱的一側(cè)遷移,當(dāng)?shù)孛惭莼聊骋浑A段時(shí),侵蝕能力強(qiáng)的水系會(huì)在某處率先侵蝕切穿分水嶺,并匯入侵蝕能力弱的水系,進(jìn)而發(fā)生流域重組.

因此,斷開的分水嶺在DEM數(shù)據(jù)上最直觀的表現(xiàn)是:某處被切穿的分水嶺的高程明顯低于兩端的高程,在地勢(shì)上呈“凹”字形態(tài),如圖5(a)所示.根據(jù)均值變點(diǎn)分析法確定的最佳集水面積閾值提取河網(wǎng)水系，并與分水線圖層進(jìn)行疊加,同時(shí)結(jié)合野外地形實(shí)地考察分析,判斷確定需要重新被連接的分水線,如圖5(b)和(c)所示.

(a) 斷開的分水嶺示意圖 (b) 無人機(jī)航拍襲奪灣位置 (c) 需連接的分水線位置
圖 5 柴汶河上游襲奪區(qū)域分水線恢復(fù)連接位置

Fig. 5 The connection position of the dividing line of the captured area in the upper reaches of the Chaiwen River

根據(jù)分水線恢復(fù)原理,對(duì)被連接分水線途經(jīng)柵格單元的高程值重新賦值.圖6中方格為柵格單元,數(shù)值為高程值.圖6(a)中填充的深色方格表示斷開的分水線兩端端點(diǎn),斜直線表示需要被重新連接的分水線,斜直線途經(jīng)的方格表示需要被重新賦值的柵格單元.由圖6(a)可知,在未恢復(fù)斷開的分水線之前,斷開的分水線兩端點(diǎn)高程值均高于需要被重新連接的分水線所途徑的高程值,呈“兩端高、中間低”的形態(tài),證明本文判斷的需要被重新連接的分水線是正確的.圖6(b)中填充的深色方格表示已被重新賦值的柵格單元,方格內(nèi)的數(shù)值為賦值后的新高程值.將分水線恢復(fù)前后的DEM高程值對(duì)比可知,原先斷開的分水線的高程值已由最初的低值重新賦值為高值,在恢復(fù)分水線前的高程值上抬高了約0.3～15.5 m,從斷開的分水線一端的高程值378 m逐漸抬高,最終逐漸接近另一端點(diǎn)高程值383 m.至此,柴汶河流域河流襲奪前的分水線已被成功恢復(fù).同時(shí),基于上述原理重新生成新的DEM,即柴汶河河流襲奪前的DEM數(shù)據(jù).

(a) 分水線恢復(fù)前的DEM值 (b) 分水線恢復(fù)后的DEM值圖 6 柴汶河上游襲奪位置分水線恢復(fù)前后的DEM值

3.3 襲奪前水系流域恢復(fù)河流襲奪是引起水系結(jié)構(gòu)改變,進(jìn)而誘發(fā)流域盆地重組的重要原因之一.恢復(fù)襲奪前的水系流域有助于進(jìn)一步了解河流襲奪產(chǎn)生的機(jī)理,以及河流襲奪對(duì)地貌產(chǎn)生的影響.流域的邊界線通常是相鄰兩個(gè)流域之間的分水線,因此根據(jù)恢復(fù)分水線成功后得到的河流襲奪前DEM數(shù)據(jù),對(duì)其進(jìn)行水系提取和流域劃分.將襲奪前的流域劃分為被奪河流域與襲奪河流域,并疊加新提取的水系圖層,從而得到柴汶河流域河流襲奪前的水系分布與流域劃分(如圖7(a)所示).為了便于對(duì)比與分析,同時(shí)提取并劃分了襲奪后的水系及流域(如圖7(b)所示).

(a) 襲奪前水系流域劃分 (b) 襲奪后水系流域劃分圖 7 柴汶河上游襲奪前后水系流域劃分

通過對(duì)襲奪前后的水系流域進(jìn)行對(duì)比分析,并根據(jù)河流襲奪前水系流域分布圖中的水系形態(tài)可以推斷得出:在河流襲奪發(fā)生前,被奪河古沂河自西南向東北方向流經(jīng)大張莊、田莊水庫后最終匯入沂河干流;河流襲奪發(fā)生后,古沂河在襲奪灣處發(fā)生了近90°的大轉(zhuǎn)彎后匯入了古柴汶河,襲奪的同時(shí)在襲奪灣生成具有高勢(shì)差的裂點(diǎn),裂點(diǎn)一分為二,分別沿被奪河、倒流河向上遷移.倒淌河因失去源頭供給,且裂點(diǎn)遷移對(duì)其河道的調(diào)整,使流域內(nèi)地勢(shì)發(fā)生倒轉(zhuǎn),呈西南低東北高,遂則改向匯入古柴汶河,成為柴汶河流域的一部分.同時(shí),河流襲奪也產(chǎn)生了斷頭河(如圖2所示),并與倒淌河之間形成了新的分水高地(即風(fēng)口),將流域分為柴汶河流域和沂河流域.

4 討論

4.1 河流襲奪與流域地貌形態(tài)演變關(guān)系流域地貌形態(tài)是由于地表受到侵蝕和堆積形成的,其形態(tài)的差異反映了流域被侵蝕的程度以及流域地貌演化的階段[27].河谷的形成與發(fā)育問題是地質(zhì)地貌學(xué)長(zhǎng)期關(guān)注的經(jīng)典課題,有關(guān)河谷發(fā)育的學(xué)說已成為地貌學(xué)發(fā)展的支持基石和理論源泉[28].在山地和丘陵地區(qū),由于地質(zhì)構(gòu)造復(fù)雜、巖層堅(jiān)硬,河谷的形態(tài)主要由地質(zhì)因素決定,表現(xiàn)為寬谷與峽谷沿河交替分布,河谷發(fā)展過程在基巖山地河谷的橫剖面發(fā)展過程中最為清楚.

以襲奪灣為界,分別選取古沂河與古柴汶河的上游部分河谷作橫剖面,剖面位置如圖7(b)所示.從橫剖面(如圖8所示)A-A’、B-B’和C-C’中可以看出,古沂河上游的河谷剖面為不對(duì)稱“V”字形.其河谷的凹岸發(fā)生了侵蝕,坡面較陡;凸岸逐漸堆積,坡面較緩,且出現(xiàn)了河灣和水小谷寬的現(xiàn)象.這是由于河流發(fā)生襲奪后,水流量減少,河流的下蝕能力減弱,側(cè)蝕能力加強(qiáng),呈現(xiàn)出原河谷不相適應(yīng)的現(xiàn)象.

圖8 柴汶河上游部分河谷橫剖面

在古柴汶河的上游,以下切侵蝕作用為主,河流深切入基巖,形成河身直,河床坡度陡,斷面狹窄的“V”型谷.剖面D-D’為“V”型谷的初級(jí)發(fā)育狀態(tài)——嶂谷,兩壁陡峭,常有基巖或礫石灘露出水面以上.嶂谷進(jìn)一步發(fā)展形成峽谷(剖面E-E’),橫剖面呈明顯的“V”字形,出現(xiàn)谷中谷現(xiàn)象,這是因?yàn)樵诤恿饕u奪未發(fā)生前,該段河谷河流流量較小,下蝕較弱,以側(cè)蝕為主,從而形成較寬的河谷地貌;在襲奪發(fā)生后,被奪河的河水匯入,使河流流量加大,流速變快,河流溯源侵蝕和下蝕較為強(qiáng)烈,從而使河床加深變窄,形成落差較大的谷中谷地貌.F-F’橫剖面處的河谷形態(tài)“V”字形逐漸展寬,河谷曲線光滑,有從“V”型谷逐漸向“U”型谷演變的趨勢(shì).

4.2 河流襲奪與河流縱剖面關(guān)系河流縱剖面是河流作用形成的,河流為了適應(yīng)外在條件的變化,每一時(shí)刻任一河段不僅進(jìn)行侵蝕,同時(shí)也發(fā)生堆積,因此河流縱剖面存在自動(dòng)調(diào)整作用[29].同樣,以襲奪灣為界,分別提取古沂河和古柴汶河上游的河流縱剖面(如圖9所示).結(jié)果可知，山區(qū)相較于平原地區(qū)的河流縱剖面較陡,起伏較大,多跌水和瀑布,兩岸常有許多山嘴突出,使河床岸線犬牙交錯(cuò),呈凹凸不規(guī)則型.河流縱比降不僅反映出流經(jīng)該河床的水流流速,還反映巖性的變化或構(gòu)造運(yùn)動(dòng),同時(shí)也能間接反映河流對(duì)地面侵蝕的能力.通過計(jì)算分別得出被奪河的河流縱比降為10.7‰,襲奪河的河流縱比降為15.5‰,古柴汶河上游的河流縱比降大于古沂河上游.當(dāng)流域面積不大且其他條件相當(dāng)時(shí),河流縱比降越大,水流速度就越大,代表水動(dòng)能就越大,河流對(duì)地面的下切侵蝕能力就越強(qiáng)[30].古柴汶河上游對(duì)地面的侵蝕能力要大于古沂河上游,因此能率先切穿分水嶺,襲奪古沂河上游.結(jié)合上述對(duì)其襲奪前水系流域的恢復(fù)研究可知,該研究結(jié)果是對(duì)海拔較低、侵蝕能力強(qiáng)的古柴汶河向上溯源侵蝕襲奪了海拔較高、侵蝕能力弱的古沂河的又一佐證.

圖9 柴汶河上游部分河流縱剖面

5 結(jié)論

本文主要基于ArcGIS軟件為技術(shù)支撐,結(jié)合野外調(diào)查,利用DEM基礎(chǔ)數(shù)據(jù)對(duì)柴汶河襲奪河流域進(jìn)行河網(wǎng)水系的提取以及河流襲奪前水系流域的恢復(fù)重建,取得了較好的結(jié)果.主要認(rèn)識(shí)如下:

1) 在利用ArcGIS軟件提取河網(wǎng)水系的過程中,基于不同集水面積閾值提取河網(wǎng)水系時(shí)會(huì)得到不同的結(jié)果,因此會(huì)對(duì)后續(xù)的水文分析產(chǎn)生影響.本文采用均值變點(diǎn)分析法確定了柴汶河流域的最佳集水面積閾值為0.54 km2,即設(shè)置匯流累積量為600時(shí),提取的河網(wǎng)水系最接近于實(shí)際水系,進(jìn)一步印證了均值變點(diǎn)分析法在確定最佳集水面積閾值時(shí)具有一定的可靠性.

2) 在對(duì)襲奪前流域的重建過程中,主要基于二次曲線內(nèi)插法來恢復(fù)其襲奪前的分水嶺,得到新的流域邊界,完成襲奪前水系的重建;同時(shí)尚未考慮其區(qū)域巖性、地質(zhì)構(gòu)造、構(gòu)造運(yùn)動(dòng)、氣候變化等因素的影響.最終得到的結(jié)果說明這一想法應(yīng)用在襲奪前河流的恢復(fù)重建方面是科學(xué)的,并且為襲奪前流域的恢復(fù)重建提供了重要進(jìn)展.

3) 通過對(duì)柴汶河襲奪前水系流域的恢復(fù)重建以及流域地貌演化分析,河流襲奪對(duì)于被奪河與襲奪河的河谷形態(tài)發(fā)育產(chǎn)生了重要的影響,同時(shí)古柴汶河上游的河流縱比降大于古沂河上游,對(duì)地面侵蝕能力強(qiáng),故能首先切穿分水嶺,襲奪古沂河,該研究結(jié)果對(duì)位于較低的河流通過向上溯源侵蝕襲奪較高的鄰近河流提供了又一例證.