跌水高度對元謀干熱河谷沖溝溝頭侵蝕產沙特征的影響初探*

張寶軍熊東紅楊 丹張 素校 亮蘇正安董一帆

（1 中國科學院山地災害與地表過程重點實驗室，成都 610041）

（2 中國科學院水利部成都山地災害與環境研究所，成都 610041）

（3 中國科學院大學，北京 100049）

跌水高度對元謀干熱河谷沖溝溝頭侵蝕產沙特征的影響初探*

張寶軍1，2，3熊東紅1，2?楊 丹1，2，3張 素1，2，3校 亮1，2，3蘇正安1，2董一帆1，2

（1 中國科學院山地災害與地表過程重點實驗室，成都 610041）

（2 中國科學院水利部成都山地災害與環境研究所，成都 610041）

（3 中國科學院大學，北京 100049）

元謀干熱河谷沖溝侵蝕強烈，以溝頭溯源侵蝕過程為主，陡立溝壁跌坎是活躍溝頭的重要形態特征之一。采用野外原位放水沖刷試驗，研究了25和50 cm兩種不同跌水高度對侵蝕產沙過程的影響。結果表明：2種跌水高度下，（1）徑流跌水沖刷力差異顯著。跌水勢能轉化量平均值分別為4.89和9.78 J s-1，跌水剪切力平均值分別為25.9和53.5 Pa，均呈現倍數關系。（2）跌穴發育形態特征及下游侵蝕量差異大。2個小區跌水最大下切深度分別為7.38和7.50 cm，平均下切深度分別為3.87 和5.16 cm；溝頭溝壁及下游溝床部位侵蝕量累積貢獻平均為41%和54%。（3）徑流含沙量差異顯著（p＜ 0.1）。整個試驗中2個小區徑流平均含沙量分別為7.51和18.76 g L-1，且徑流經過溝壁跌坎后，含沙量分別平均增大6.25 和25.49 g L-1。初步認為溝壁跌坎高度差異影響徑流沖刷力及其侵蝕產沙特征，但仍需進一步開展更多跌坎高度下的相關研究，為干熱河谷區沖溝溝頭溯源侵蝕動力學機制研究提供補充。

徑流沖刷力；侵蝕產沙特征；跌水高度；沖溝溝頭；元謀干熱河谷

土壤侵蝕是當今世界普遍關注的重大環境問題之一，由于土壤侵蝕導致的嚴重水土流失制約著區域農業經濟發展［1］。沖溝侵蝕作為一種重要的土壤侵蝕方式，受到國內外學者的廣泛關注［2］，已有學者針對溝頭溯源侵蝕過程展開了較多研究［3-4］，并認為這一過程對沖溝發生發展極為重要。在溝頭溯源侵蝕發育過程中，常由于土壤異質性和徑流不均一、不穩定等原因，導致溝頭形態變化不均衡，逐漸產生跌水落差［5］。當坡面集中流到達跌坎時，由于地形突變轉變為跌水（Jet flow），伴隨著徑流勢能向動能轉化，對跌坎底部土體直接沖蝕，導致溝頭形態變化劇烈［6］。Dey等［7］研究認為，在某種程度上跌水高度和徑流量可以決定溝頭溯源侵蝕過程。Bennett和Casalí［8］采用模擬水槽沖刷試驗，研究了不同跌水高度下（5～50 mm）的溝頭溯源侵蝕過程，發現隨著跌水高度增加，溯源侵蝕速率、最大下切深度和產沙量均逐漸增大。

金沙江干熱河谷區是我國西南特殊的脆弱生態環境類型區，該區域沖溝侵蝕問題突出，水土流失強烈，嚴重威脅著區域經濟發展和下游水利水電工程安全［9］。由于該區域巖層特性復雜、降雨季節性差異等原因，導致沖溝溝頭具有形態復雜多樣、溝壁跌坎陡立、溯源侵蝕速率快等特征［10-11］。目前，已有學者針對該區域沖溝侵蝕過程中的水動力條件及其侵蝕產沙效應展開了較多研究［12-14］。Su等［15-16］研究認為干熱河谷沖溝形態變化及侵蝕產沙，主要集中表現在一條溝上若干個活躍的溝頭部位，溝頭溯源侵蝕產沙約占整個沖溝活躍區產沙的62%～75%。溝壁跌坎作為溝頭重要形態特征參數，在溝頭不同發育階段表現出不同高度大?。?7］，而不同跌坎高度發育直接影響徑流沖刷特性，導致溝頭侵蝕產沙過程差異，但其具體如何影響溝頭徑流沖刷力及其侵蝕產沙效應尚不得而知，需開展深入研究。本研究采用野外原位放水沖刷試驗，結合近景攝影測量技術，初步研究了跌水高度差異對溝頭溯源侵蝕過程中的徑流水動力條件、形態變化及侵蝕產沙特征的影響，旨在為沖溝溯源侵蝕動力學機制研究提供補充，同時為該區域進行植被恢復實踐控制溝頭溯源侵蝕提供理論依據。

1 材料與方法

1.1 研究區概況

研究區位于云南元謀干熱河谷（25°23′～26°06′N，101°35′～102°06′E），該區屬南亞熱帶季風氣候，具有炎熱干燥、降水集中、干濕季分明的氣候特征：年均氣溫21.9℃，年均降水量615 mm，主要集中于6—10月的雨季，約占年降水量的90%，年均蒸發量高達3 911 mm，約為年均降水量的6.4倍，年均干燥度為2.8［18］。土壤類型以燥紅土和變性土為主，植被以稀樹灌木草叢為主，森林覆蓋率僅為3.4%～6.3%［19］。該區沖溝極為發育，水土流失嚴重，土壤侵蝕模數高達1.64×104t km-2a-1，溝壑密度3.0～5.0 km km-2，最大達7.4 km km-2［10］（圖1a）。沖溝發育的元謀組地層為第四紀河湖相沉積物，厚達695 m，具有結構松散、膠結度差、黏砂互層等特征［20］。特殊的巖土性質和氣候條件，決定該區域沖溝發育主要以溝頭溯源侵蝕過程為主，年均溯源侵蝕速率為50 cm a-1左右，最大達200 cm a-1。此外，由于溯源侵蝕過程中土體崩塌發生頻繁，溝頭常發育有陡立跌坎［21］（圖1b），其高度大小往往可以表征溝頭的不同活躍程度［17］。

圖1 元謀干熱河谷沖溝侵蝕強烈（a）及溝頭陡立跌坎發育（b）Fig. 1 Intensive gully erosion（a）and development of steep gully headcut（b）in Yuanmou dry-hot valley region

1.2 試驗小區

本研究選擇在中國科學院水利部成都山地災害與環境研究所與云南省農業科學院熱區生態農業研究所合建的元謀干熱河谷溝蝕崩塌觀測研究站（以下簡稱“元謀站”）內進行。2015年7月，在元謀站內選取了一個雨季常有崩塌發生、裸露的原位活躍溝頭，基于原始跌坎形態，修建了跌坎高度分別為25 cm和50 cm的兩個模擬試驗小區，其中溝壁人工一致修整為“平整、陡立、無內凹洞”的初始侵蝕形態。模擬試驗小區還包括穩流槽、上游集水區和下游溝床，其中集水區為裸露光板地，長度為5 m、寬度為1.5 m、坡度為11°；溝床長度為2 m、寬度為1.5 m、坡度為9°。穩流槽采用半開口設計，集水區和溝床部位采用淺V型設計，以模擬集中徑流沖刷過程。為了監測沖刷試驗過程中的水動力參數，集水區從穩流槽開口處，以1 m為間隔，設置0～1、1～2、2～3、3～4和4～5共計5個觀測斷面；溝床從溝壁跌坎底部開始，以1 m為間隔，設置0～1和1～2共計2個觀測斷面（如圖2所示）。野外試驗之前，在試驗小區集水區和溝床部位分別隨機采集表層土壤混合樣（0～5 cm），帶回實驗室采用常規方法分析土壤理化性質，包括土壤機械組成、容重和總孔隙度，具體見表1。

表1 試驗小區土壤基本理化性質Table 1 Basic physico-chemical properties of the experimental soil

圖2 試驗小區及放水沖刷試驗示意圖Fig. 2 Sketch of the simulated scouring experiment with two plots

1.3 試驗設計

本試驗采用野外模擬放水沖刷方法研究溝頭溯源侵蝕過程特征，屬連續觀測試驗，但需要說明的是，由于試驗小區構建困難及原位控制試驗難度大，未進行重復試驗。參照該研究區典型暴雨強度60 mm h-1（據元謀站氣象數據），及活躍溝頭發育的上游平均匯水面積120 m2計算［21］，設置放水流量為120 L min-1，流量通過試驗小區上方設計的穩流池和流量計控制（圖2）。徑流進入試驗小區前先通過穩流槽，可以保證徑流進入小區時的初始流速基本一致。整個試驗過程持續140 min，其中為監測溝頭地形變化將試驗分5次進行，第1次沖刷時間為20 min，之后4次分別為30 min。試驗開始前，在小區內預先均勻撒水，直至表面充分濕潤但又無地表徑流產生。

試驗中觀測徑流水動力參數和產沙特征，試驗之后監測小區地形變化。各斷面徑流水動力參數（徑流深、徑流寬和流速）的變化，按照一定時間間隔每次分別在斷面上、中、下部位測量3次取平均值［16］。考慮到試驗初期，小區表面地形變化明顯，可能導致徑流水動力特性變化大，因此在前20 min試驗中設置5 min的監測間隔，之后試驗設置10 min的監測間隔。水動力參數測定方法與參考文獻［16］一致，徑流寬采用精度1 mm薄鋼尺測定，徑流深采用精度1 mm直尺測定，流速采用KMnO4染色法測定。文中在分析徑流水動力特性時為單次試驗數據均值。試驗過程中的徑流產沙特征，采用500 ml集流瓶收集跌水和小區集流槽內徑流泥沙樣進行觀測，分別表征集水區和整個試驗小區的徑流含沙量。在每次沖刷試驗中，自產流開始至10 min內每2 min采集1次，之后每5 min采集1次，試驗結束后帶回實驗室采用烘干法測定，并換算為1 L的泥沙含量。

1.4 徑流水動力特性

坡面徑流到達跌坎后轉變為跌水，從能量角度出發，徑流跌水的沖刷能力可以采用勢能、動能和剪切力表示［22］。徑流轉變為跌水時，大量勢能向動能轉化，徑流勢能轉化量（ΔEp）可以采用公式計算：

式中，ρ為水的容重（k g m-3），q為徑流量（m3s-1），H為跌水落差（m），g為重力加速度（m s-2）。坡面徑流到達跌坎上沿轉變為射流時的水平出射動能可以表述為：

式中，Vbrink為射流水平出射流速［23］（m s-1），與坡面徑流流速和流態有關。坡面徑流流態采用弗羅德數（Fr）參數表征，它表示過水斷面上徑流的動能和勢能的對比關系，計算公式為：

式中，V為坡面徑流到達跌坎前的平均流速（m s-1），h為平均徑流深（m）；當Fr＞1時為急流，當Fr≤1時為緩流［24］。Vbrink可以采用下式計算［25］：

徑流跌水侵蝕能力也可以采用最大剪切力表示［26］，采用公式計算：

式中，Cf為摩擦系數：

式中，ν為徑流運動粘滯系數（m2s-1），主要與徑流溫度有關；q′為單寬徑流量（m2s-1）。

1.5 地形變化

沖刷試驗過程中試驗小區的地形及侵蝕量變化，采用高精度近景攝影測量技術獲取。近景攝影測量技術是一種基于運動結構（Structure from Motion，SfM）和多視角立體成像技術（MultiView-Stereo，MVS）發展形成的三維影像重建技術，它可以生成高質量、高分辨率的三維點云數據，最高精度可達毫米級［27］。沖刷試驗之前，在試驗小區邊界以1 m左右間隔設置約15個控制點，各控制點坐標（x，y，z）借助RTK-GPS測定（Trimble R8，美國，水平精度1 cm± 1 ppm，垂直精度2 cm± 1 ppm）［21］。在沖刷試驗之前和5次沖刷試驗后，分別采用Nikon D610相機從不同角度定焦（50 mm）采集小區地形照片。地形數據處理借助Agisoft PhotoScan 1.1.6 專業版軟件進行，該軟件可以實現三維模型重建過程的完全自動化。每次選取約30張采集照片導入軟件（足以覆蓋整個試驗小區），并添加控制點坐標進行對齊、優化，最后自動生成真實坐標的3D點云模型。在12次點云數據建立中，軟件模型依據控制點高程計算的均方根誤差（RMSE）介于0.002～0.009 m，精度較高（圖3）。最后將生成的點云數據導入ArcGIS 10.1軟件，利用3D分析工具中的TIN模型插值方法，生成數字高程模型（DEM），通過對比不同階段DEM分析小區地形變化量。

2 結 果

2.1 徑流水動力學特性

圖3 2號小區第2次試驗后的3D點云模型Fig. 3 3D point cloud model of Plot 2 after the second test

圖4 試驗小區坡面徑流流速（a）和弗羅德數（b）Fig. 4 Runoff flow velocity（a）and Froude number（b）of the two plots

圖4a、圖4b分別為試驗小區集水區和溝床徑流流速和弗羅德數隨時間的變化規律。可以看出，整個沖刷試驗過程，徑流流速和弗羅德數隨時間均呈現先減小后趨于穩定的趨勢，且該兩個參數均表現為集水區部位顯著高于溝床部位（p＜0.05）。其中，2個小區集水區平均流速分別為0.58和0.53 m s-1，溝床平均流速分別為0.40和0.43 m s-1。徑流Fr變化范圍在1.02～2.54之間，均大于1，屬于急流狀態。徑流流速和Fr隨時間的變化趨勢與侵蝕發育過程有關。在試驗初始階段，小區坡面平整，徑流流速大、流路寬、流深淺，弗羅德數大；之后隨著沖刷進行，表層松散土體被迅速侵蝕，并形成明顯溝道，徑流深度增加，加之徑流侵蝕不斷消耗能量，導致徑流流速變??；最后，小區坡面形態趨于穩定，徑流流速和弗羅德數也趨于穩定狀態。

由于試驗采用固定放水流量0.002 m3s-1，所以徑流在溝頭跌坎部位的勢能轉化量主要與跌水高度有關，動能主要與流速有關?？紤]到徑流流速和流態在上游集水區的沿程變化，采用集水區最下游斷面流速和弗羅德數，模擬計算徑流到達跌坎上沿的出射流速，并基于此分析徑流由坡面流轉換為跌水時的出射動能（圖5）。2個小區徑流到達跌坎上沿的水平方向出射流速和動能，與坡面流速和弗羅德數的變化規律一致，同樣表現為先減小后趨于穩定的趨勢，且2個小區之間無顯著差異（p＞0.1）。

圖5 徑流轉變為跌水的初始流速和動能Fig. 5 Initial flow velocity and kinetic energy while surface turning into jet

圖6a、圖6b分別表示坡面徑流轉變為跌水時的勢能轉化量及其剪切力變化。由于試驗小區集水區為裸露光板地，表層有鐵錳結膜覆蓋，質地堅硬，導致徑流下切侵蝕困難。沖刷試驗過程中，2個小區跌坎上沿徑流最大下切深度不足2 cm，對跌水高度變化影響較小，因此在計算跌水勢能轉化量時統一采用初始跌水高度。從圖6a可以看出，跌水高度50 cm的2號小區，其勢能轉化量顯著高于1號小區，兩者平均值分別為9.78和4.89 J s-1。基于能量守恒定律，并忽略徑流跌落過程中的能量耗散，徑流勢能將主要轉化為動能，增大跌水流速，經計算，2個小區跌水到達底部時的瞬時流速達到2.28 和3.18 m s-1，分別為出射流速的4.1倍和5.6倍。此外，分析2個小區跌水剪切力，平均值分別為25.9和53.5 Pa，2號小區同樣顯著高于1號小區（圖6b）。以上結果說明，坡面徑流到達溝壁跌坎這一地形陡變的特殊部位時，徑流勢能向動能轉化，導致徑流沖刷力顯著增大，并且隨著跌水高度的成倍增加，徑流勢能轉化量和跌水剪切力大小也將呈現倍數增大。

圖6 徑流跌水勢能轉化量（a）和剪切力（b）Fig. 6 Potential energy transformation（a）and shear stress（b）of jet flow

2.2 地形及侵蝕量變化特征

徑流跌水對溝壁跌坎底部土體的直接沖蝕，導致跌水坑（Plunge pool）發育。隨著沖刷進行，跌水坑寬度和深度逐漸增大，最終促進溯源侵蝕的發生。此次沖刷試驗過程中，由于2個小區集水區徑流流路偏向右側，導致跌水侵蝕主要發生在溝壁底部靠右1 m范圍內。根據試驗之前和每次沖刷試驗后的6次DEM數據，提取了2個小區每次試驗后跌穴形態發育特征（表2），并選取了沖刷試驗前、80 min和140 min試驗后地形數據，繪制了2個小區距離溝壁底部5cm處的溝床橫向剖面線，以直觀描述跌水坑的形態發育特征（圖7）。從圖7中看出，2個小區在試驗中發生跌水的明顯下切，導致跌穴發育，且1號小區更接近V型發育，2號小區更接近U型發育。而從表2中可以看出，隨著沖刷試驗進行，2個小區跌水最大和平均下切深度均逐漸增大，且2號小區平均下切深度顯著大于1號小區。其中1、2號小區在140 min試驗后跌穴最大下切深度分別為7.38 cm和7.50 cm，差異不顯著，但中間變化過程卻表現出極大差異。2號小區前20 min試驗后最大下切深度達5.81 cm，之后的120 min試驗僅增大了1.69 cm，而1號小區前20 min試驗后最大下切深度僅為1.65 cm，之后的120 min試驗逐漸增大至7.38 cm。2個小區平均下切深度隨時間也表現為相似的增長特征，在前20min試驗后平均下切深度分別為0.70 cm和3.69 cm，相差4倍左右，而在140 min試驗后分別為3.87 cm和5.16 cm，約為1.3倍。以上結果說明徑流經過不同的跌水高度后，產生不同的沖刷能力，導致溝壁跌坎底部土體沖刷侵蝕特征差異明顯。

表2 跌水下切深度隨沖刷時間的變化Table 2 Variation of incision depth of the jet flow with the scouring going on（cm）

圖7 跌穴形態發育特征Fig. 7 Morphological characteristics of plunge pools

此外，通過分析溝壁跌坎及下游溝床地形侵蝕量，發現2個小區之間差異顯著（圖8）。在初始20 min沖刷后，1號小區侵蝕量僅為0.003 m3，但2號小區高達0.067 m3；最后經過140 min沖刷，1號小區侵蝕量增加至0.072 m3，2號小區達到0.115 m3，約為1號小區的1.6倍。分析每次沖刷試驗后，溝壁跌坎及下游溝床侵蝕量對總侵蝕量的累積貢獻，發現1號小區在20 min時為34%，之后逐漸增長至45%左右，總體平均為41%，而2號小區在20 min時達到65%，之后降低至50%左右，平均為54%（表3）。據此可以推測，跌坎高度對溝頭溝壁及下游侵蝕產沙具有一定影響，但是否可以說明跌水高度越高，溝壁跌坎及下游溝床侵蝕量越大，對總侵蝕量的貢獻也越大，這一結論還需要在以后設置其他跌坎高度進行補充驗證。

圖8 跌坎上沿以下部位侵蝕量的變化Fig. 8 Soil erosion volume in sections below the border-top of headcut

表3 溝壁跌坎及溝床侵蝕量占小區總侵蝕量的比例Table 3 Contribution of the soil erosion at the gully headcut and in the bed to total soil erosion of the plot（%）

2.3 徑流輸沙特征

對于整個沖刷過程，2個小區徑流泥沙含量總體上隨時間呈現冪函數減小的趨勢，擬合曲線分別為：①y=53.39 x-0.48，N=49，R2=0.86；②y=93.09 x-0.36，N=49，R2=0.71。在整個沖刷試驗產流開始，2個小區徑流泥沙含量分別達到了最大值229.1 g L-1和259.8 g L-1，但經過10 min后分別急劇下降至8.5 g L-1和34.6 g L-1，并最終呈現趨于穩定的趨勢（圖9a）。但從圖9a又可以明顯看出，之后試驗過程中出現明顯的4次含沙量極大值點，但其大小隨沖刷次數呈現逐漸減小的趨勢。這是由于2次試驗間隔內小區表面受到干濕交替影響，導致表面土體松散，下次試驗產流開始時表土在集中徑流作用下極易被帶走。之后隨著表面浮土被沖刷，徑流泥沙含量又逐漸降低至穩定。

采用單因素方差分析進行均值比較，發現徑流到達溝頭跌坎時2個小區之間含沙量無顯著差異，但小區末端2個小區徑流含沙量卻在0.1水平上差異顯著，說明經過不同跌水高度后，徑流泥沙含量變化較大。從圖9a可以看出，在5次沖刷試驗產流開始時，2號小區徑流泥沙含量均顯著高于1號小區，特別是第2次和第3次試驗，1號小區泥沙含量分別為57.34 g L-1和31.80 g L-1，而2號小區高達126.02 g L-1和78.12 g L-1，均為1號小區的2倍多。除去首次試驗產流開始時的含沙量最大值及之后4次試驗的極大值，整個沖刷試驗過程中，2個小區徑流平均含沙量分別為7.51 g L-1和18.76 g L-1，相差約1.5倍。通過計算小區出口和集水區末端徑流含沙量差值，表示徑流經過溝壁跌坎后的泥沙變化量（圖9b）。可以看出，前80 min沖刷試驗中，2號小區徑流泥沙增大量明顯高于1號小區，最高增大了108.8 g L-1，平均增大25.49 g L-1，而1號小區最高僅增大了25.28 g L-1，平均增大6.25 g L-1，這說明徑流經過更高的跌水落差后，徑流沖刷力增大，導致下游更多泥沙被徑流帶走。

3 討 論

圖9 徑流總含沙量（a）和跌坎上沿以下徑流含沙量變化量（b）Fig. 9 Total runoff sediment concentration（a）and change of sediment concentration along the headcut（b）

陡立跌坎是干熱河谷沖溝溝頭的重要形態特征，其高度大小常被用作溝頭活躍程度的野外定性判定標準之一［17］。當坡面匯流到達跌坎時，會由于地形的突然轉折轉變為噴射流，伴隨著勢能向動能的快速轉化，對溝壁底部土體產生嚴重沖刷，促進溝頭溯源侵蝕發生發展。野外多年地形監測數據也表明，跌坎高度大的溝頭，年均溯源侵蝕速率越大，形態變化越明顯。本文研究了2種跌水高度下的溝頭侵蝕產沙過程特征，發現更高跌水高度下的徑流沖刷力、跌水下切深度和徑流含沙量均表現為更大，說明跌坎高度對溝頭溯源侵蝕特征會有影響，這一初步認識對干熱河谷溝頭溯源侵蝕過程研究具有重要的補充意義。同時本文試驗結果與Bennett和Casalí［8］研究基本一致，他們認為溯源侵蝕速率、最大下切深度和產沙量均隨跌水高度增加而逐漸增大，這一結論是基于室內水槽模擬沖刷試驗，研究了5～50 mm內10種跌水高度下溝頭溯源侵蝕過程得出的。與之相比，本研究不足之處在于跌水高度設置過少，但優勢在于采用的是野外原位控制實驗，更接近和反映自然真實情況，且跌水高度尺度更大，為其最大跌水高度的5倍和10倍。

盡管本研究初步查明了跌水高度會影響干熱河谷溝頭溯源侵蝕過程特征，但試驗設計上還存在一些不足，需要進一步改進和完善。首先，本研究選擇基于原位溝頭構建試驗小區，由于溝頭寬度和集水區面積限制，僅允許設置2個小區，如果增加小區設計就需要另外選擇溝頭。最終為保證試驗一致性，在溝頭原始跌坎基礎上暫時只設計了25 cm 和50 cm 兩種跌水高度。但兩種跌坎高度處理不足以得出非常確定的趨勢性結論，特別是更大跌水高度下，徑流沖刷力和侵蝕產沙特征呈現何種變化規律？是否存在影響的臨界跌坎高度？這些均需要在以后設置其他跌坎高度進行補充驗證。其次，本試驗是基于原位溝頭開展的野外控制試驗，其試驗難度遠大于室內水槽試驗，未能有效地進行重復試驗，以后試驗需考慮設計多組試驗模擬重復。最后，本試驗實際為溝頭溯源侵蝕過程的連續觀測試驗，但為了監測試驗過程中的地形變化，按照經驗將試驗分5次進行，其中會涉及干濕交替問題，對溝頭溯源侵蝕過程可能產生影響，本文對此暫未考慮。

基于以上分析，研究者將于后續研究中不斷完善本試驗中的不足，特別是增加其他跌坎高度設置，以期獲得更好、更完善的理論數據，旨在查明跌坎高度對溝頭溯源侵蝕過程的影響機理，為干熱河谷溝頭溯源侵蝕研究做出更有意義的補充。

4 結 論

本文采用野外放水沖刷試驗方法，初步研究了溝頭溝壁不同跌水高度對侵蝕產沙過程的影響。對比25 cm和50 cm跌坎高度，發現徑流跌水沖刷能力差異顯著，2種跌水高度下徑流平均勢能轉化量分別為4.89和9.78 J s-1，跌水平均剪切力分別為25.9 和53.5 Pa，均呈現倍數增大。徑流跌水沖刷力差異，導致溝壁底部跌穴發育形態特征及下游侵蝕量貢獻顯著不同。沖刷試驗中，2個小區跌水下切深度變化過程差異極大，試驗結束后最大下切深度分別為7.38和7.50 cm，平均為3.87和5.16 cm；2個小區溝壁跌坎及溝床對總侵蝕量的平均貢獻為41%和54%。此外，2個小區徑流含沙量在試驗中同樣差異顯著，除去試驗中含沙量最大值和極大值，平均含沙量分別為7.51和18.76 g L-1；特別地，在試驗前80 min，徑流經過溝壁跌坎后，2個小區徑流含沙量分別平均增大6.25 和 25.49 g L-1。以上研究結果說明，跌坎高度影響溝頭溯源侵蝕過程中的徑流沖刷力及其侵蝕產沙特征，但其中內在影響機理及更大跌坎高度下的溝頭溯源侵蝕特征還需進一步深入研究。

［1］Guo Q K，Hao Y F，Liu B Y. Rates of soil erosion in China：A study based on runoff plot data. Catena，2015，124：68—76

［2］Poesen J，Nachtergaele J，Verstraeten G，et al. Gully erosion and environmental change：Importance and research needs. Catena，2003，50：91—133

［3］Frankl A，Poesen J，Deckers J，et al. Gully head retreat rates in the semi-arid highlands of Northern Ethiopia. Geomorphology，2012，173/174：185—195

［4］Francis K R，Greg E T. The evolution of gully headcut morphology：A case study using terrestrial laser scanning and hydrological monitoring. Earth Surface Processes and Landforms，2015，40：1304—1317

［5］Prasad S N，R?mkens M J M. Energy formulations of head cut dynamics. Catena，2003，50（2/4）：469—487

［6］Flores-Cervantes J，Istanbulluoglu E，Bras R. Development of gullies on the landscape：A model of headcut retreat resulting from plunge pool erosion.Journal of Geophysical Research，2006，111（F1）：1—14

［7］Dey A K，Tsujimoto T，Kitamura T. Experimental investigations on different modes of headcut migration. Journal of Hydraulic Research，2007，45：333—346

［8］Bennett S J，Casalí J. Effect of initial step height on headcut development in upland concentrated flows. Water Resources Research，2001，37（5）：1475—1484

［9］Xiong D H，Yan D，Long Y，et al. Simulation of morphological development of soil cracks in Yuanmou Dry-hot Valley Region，Southwest China. Chinese Geographical Science，2010，20（2）：112—122

［10］柴宗新，范建容，劉淑珍. 金沙江下游元謀盆地沖溝發育特征和過程分析. 地理科學，2001，21（4）：339—343

Chai Z X，Fan J R，Liu S Z. Analysis on development characteristics and process of gully in Yuanmou Basin on lower reaches of Jinsha River（In Chinese）. Scientia Geographica Sinica，2001，21（4）：339—343

［11］王小丹，鐘祥浩，范建容，等. 金沙江干熱河谷元謀盆地沖溝溝頭形態學特征研究. 地理科學，2005，25 （1）：63—67 Wang X D，Zhong X H，Fan J R，et al. Study on the morphological characteristics of the gully heads in Yuanmou Basin，Arid River Valley of Jinsha River，China（In Chinese）. Scientia Geographica Sinica，2005，25（1）：63—67

［12］熊東紅，楊丹，翟娟，等. 元謀干熱河谷沖溝溝頭徑流水動力學特性及產沙效應初探. 水土保持學報，2012，26（6）：52—56

Xiong D H，Yang D，Zhai J，et al. Preliminary study on hydrodynamic characteristics of overland flow and sediment yields of gully heads in Yuanmou Hotdry Valley（In Chinese）. Journal of Soil and Water Conservation，2012，26（6）：52—56

［13］楊丹，熊東紅，張寶軍，等. 溝床草被對干熱河谷沖溝產沙特性影響的野外模擬試驗. 農業工程學報，2015，31（15）：124—132

Yang D，Xiong D H，Zhang B J，et al. Field experiment on impacts of grass belt length on characteristics of sediment yields and transport rates for gullies in Jinsha dry-hot valley region（In Chinese）. Transactions of the Chinese Society of Agricultural Engineering，2015，31（15）：124—132

［14］楊丹，熊東紅，張素，等. 金沙江干熱河谷沖溝溝床草被對徑流水動力特性的影響. 水土保持學報，2014，28（6）：77—81

Yang D，Xiong D H，Zhang S，et al. Effects of grass belt length to runoff hydrodynamic properties in gully beds in Jinsha River Dry-hot Valley Region（In Chinese）. Journal of Soil and Water Conservation，2014，28（6）：77—81

［15］Su Z A，Xiong D H，Dong Y F，et al. Simulated headward erosion of bank gullies in the dry-hot valley region of southwest China. Geomorphology，2014，204：532—541

［16］Su Z A，Xiong D H，Dong Y F，et al. Hydraulic properties of concentrated flow of a bank gully in the dry-hot valley region of Southwest China. Earth Surface Processes and Landforms，2015，40：1351—1363

［17］李佳佳，熊東紅，盧曉寧，等. 基于RTK－GPS技術的干熱河谷沖溝溝頭形態特征. 山地學報，2014，32 （6）：706—716

Li J J，Xiong D H，Lu X N，et al. Morphological characteristics of the Gully Head in Dry-hot Valley based on the RTK-GPS technology（In Chinese）. Journal of Mountain Science，2014，32（6）：706—716

［18］劉剛才，鄧偉，文安邦，等. 試論金沙江河谷建立溝蝕崩塌科學觀測研究站的重要性及學科方向. 山地學報，2010，28（3）：333—340

Liu G C，Deng W，Wen A B，et al. A review on the essentiality and subject orientation of establishing gully erosion and collapse experimental station in Jinsha River Valley（In Chinese）. Journal of Mountain Science，2010，28（3）：333—340

［19］楊艷鮮，紀中華，方海東，等. 元謀干熱河谷旱坡地復合生態農業模式效益研究初評. 水土保持研究，2005，12（4）：88—89，99

Yang Y X，Ji Z H，Fang H D，et al. Study and preliminary evaluation on the benefit of models of compound eco-agriculture on dry slope land in Yuanmou Dry Hot Valley（In Chinese）. Research of Soil and Water Conservation，2005，12（4）：88—89，99

［20］錢方，凌小惠. 元謀土林成因及類型的初步研究. 中國科學：B輯，1989（4）：412—418，449—450

Qian F，Ling X H. Preliminary study on the types and causes of Tulin in Yuanmou（In Chinese）. Science in China：Series B，1989（4）：412—418，449—450

［21］Dong Y F，Xiong D H，Su Z A，et al. Critical topographic threshold of gully erosion in Yuanmou dryhot valley in southwestern China. Physical Geography，2013，34（1）：50—59

［22］王民，崔靈周，李占斌，等. 模擬降雨條件下徑流侵蝕力與地貌特征的動態響應關系. 水利學報，2008，39 （9）：1105—1110

Wang M，Cui L Z，Li Z B，et al. Dynamic responserelationship between runoff erosivity and topographic feature under the condition of artificial simulated rainfall （In Chinese）. Journal of Hydraulic Engineering，2008，39（9）：1105—1110

［23］王海云，戴光清，張建民，等. 高水頭泄水建筑物摻氣設施研究綜述. 水利水電科技進展，2004，24（4）：46—48，56

Wang H Y，Dai G Q，Zhang J M，et al. Review of the study of aeration facilities for high water head discharge structures（In Chinese）. Advances in Science and Technology of Water Resources，2004，24（4）：46—48，56

［24］龔家國，王文龍，郭軍權. 黃土丘陵溝壑區淺溝水流水動力學參數實驗研究. 中國水土保持科學，2008，6（1）：93—100

Gong J G，Wang W L，Guo J Q. Ephemeral gully erosion experiment on hydrodynamics parameter of concentrated flow in hilly area of Loess Plateau（In Chinese）. Science of Soil and Water Conservation，2008，6（1）：93—100

［25］Hager W H. Hydraulics of plane free overfall. Journal of Hydraulic Engineering，1983，109（12）：1683—1697

［26］Stein O，Julien P，Alonso C. Mechanics of jet scour downstream of a headcut. Journal of Hydraulic Research，1993，31（6）：723—738

［27］Frankl A，Stala C，Abrahab A，et al. Detailed recording of gully morphology in 3D through imagebased modelling. Catena，2015，127：92—101

Preliminary Study on Effects of Headcut Height on Soil Erosion and Sediment Yield at Gully Heads in Yuanmou Dry-hot Valley Region

ZHANG Baojun1，2，3XIONG Donghong1，2?YANG Dan1，2，3ZHANG Su1，2，3XIAO Liang1，2，3
SU Zheng’an1，2DONG Yifan1，2
（1 Key Laboratory of Mountain Hazards and Surface Processes，Chinese Academy of Sciences，Chengdu 610041，China）
（2Institute of Mountain Hazards and Environment，Chinese Academy of Sciences，Chengdu 610041，China）
（3University of Chinese Academy of Sciences，Beijing 100049，China）

【Objective】Gully erosion is very severe in the Yuanmou Dry-hot Valley Region，dominated with headward erosion，which is characterized by steep gully headcuts. During the field investigations，headcut height is often regarded as one of the main indices of distinguishing the type of gully head activity. The abrupt drop of concentrated flow at headcut would cause a sudden change of the runoff in hydrodynamics，thus accelerating the headward erosion processes. It is，therefore，of great significance to carry out the preliminary study on effects of headcut height on characteristics of soil erosion and sediment yield in the gully headward processes，to the gully erosion research in this region as scientific supplement.【Method】To simulate different step heights，based on an in-situ active bank gully head，two experimental plots were constructed with the headcut height of 25 cm and 50 cm，separately. And in November 2015 a series of insitu scouring experiments were carried out with the flow set at 120 L min-1. The entire scouring test lasted 140 min，consisting of 5 test runs. During the scouring tests，cross-sections were set in the drainage area and along the gully bed at1 m intervals，for observation of changes of runoff in width，depth and velocity at regular intervals. Meanwhile，sediment samples were also collected at regular time intervals at the outlet of the drainage area and the plots，separately，for analysis of sediment content with the oven drying method. Besides，before and after the experiment，topographic data of the two plots were collected with the high-resolution close range digital photogrammetry measurements，for analysis of changes in landform，morphological parameters and net soil loss.【Result】Results showed that：（1）The scouring power of jet flow differed significantly between the two plots. When surface runoff turned into impinging jets，the mean potential energy transformation in the two plots was about 4.89 and 9.78 J s-1，separately with shearstress reaching 25.9 and 53.5 Pa，on average respectively，both in multiple relations. And the velocity at the instance when the flow touched the bottom of the headcut reached as high as 2.28 and 3.18 m s-1，separately，which was 4.1 and 5.6 times larger than that at the brink point respectively.（2）The two plots also differed significantly in morphological characteristics of plunge pools and soil erosion volume at lower end. At the end of the experiment，the incision depth of jet flow reached 7.38 and 7.50 cm in maximum，and was 3.87 and 5.16 cm on average，separately，in the two plots. Especially，the growth of incision depth with the experiment going on varied sharply between the two plots. Furthermore，during the scouring tests，the soil erosion cumulative contribution rate at the gully headwall and in the bed at the lower end reached as high as 45.1% and 64.9%，and was 41% and 54% on average.（3）The two plots also differed sharply in runoff sediment concentration（p＜0.1）. Excluding the maximum concentration at the beginning of each test，the mean sediment concentration was 7.51 and 18.76 g L-1，separately，about 1.5 times in difference. Especially，during the first previous 80 min of the scouring tests，when runoff flew over the headcut，the runoff sediment concentration increased by 25.28 and 108.78 g L-1in maximum and by 6.25 and 25.49 g L-1on average，separately，in the two plots. 【Conclusion】This study has tentatively concluded that the effect of gully headcut height affecting runoff scouring force and characteristics of soil erosion and sediment yield is probably the key factor influencing the headward erosion process. However，it is still not very clear about how the effect would vary in the case of higher headcuts. More experiments with a variety of headcut heights should be conducted，in an attempt to supplement more complete meaningful theories to the study on dynamics and mechanisms of gully headward erosion in dry-hot valley regions.

Runoff erosivity；Characteristics of soil erosion and sediment yield；Headcut height；Gully head；Yuanmou Dry-hot Valley Region

S157.1

A

10.11766/trxb201604240209

（責任編輯：檀滿枝）

* 國家重點基礎研究發展計劃（973）項目（2015CB452704）、國家自然科學基金項目（41571277）和中國科學院“西部之光”重點項目（Y4R2060060）共同資助 Supported by the National Basic Research Program of China（973 Program）（No. 2015CB452704），the National Natural Science Foundation of China（No. 41571277）and the Key Programme of the “Western Light” Talents Cultivation programme of the Chinese Academy of Sciences（No. Y4R2060060）

? 通訊作者 Corresponding author，E-mail：dhxiong@imde.ac.cn

張寶軍（1990—），男，山東臨沂人，博士研究生，主要從事沖溝侵蝕機理研究。E-mail：zhangbaojunaixiao@163. com

2016-04-24；

2016-07-18；優先數字出版日期（www.cnki.net）：2016-08-29