白茹茹， 張加瓊,2， 鄧鑫欣， 李志鳳， 南 瓊

(1.西北農林科技大學水土保持研究所，黃土高原土壤侵蝕與旱地農業(yè)國家重點實驗室，陜西 楊凌 712100；2.中國科學院水利部水土保持研究所，陜西 楊凌 712100；3.西北農林科技大學資源與環(huán)境學院,陜西 楊凌 712100)

土壤交錯/復合侵蝕是在水力、風力、重力、凍融等多種營力共同或交錯作用下，土壤被破壞、剝離和搬運的復雜過程；交錯/復合侵蝕過程有著顯著的疊加放大作用，以時間交錯、空間疊加為典型特征的風水交錯侵蝕是最為典型的類型之一。風水交錯侵蝕主要發(fā)生在生態(tài)環(huán)境較脆弱的干旱與半干旱地區(qū)，該區(qū)域的土壤侵蝕程度較單一風蝕或水蝕影響區(qū)域更嚴重。因此，加強風水交錯侵蝕的特征、過程、機理等研究，有益于干旱與半干旱區(qū)域水土保持的科學開展。

前人通過野外觀測(小區(qū)法)、室內模擬試驗(風洞和降雨模擬試驗)、元素示蹤等方法對不同過程的風水交錯侵蝕(風蝕后水蝕，水蝕后風蝕)的速率、風蝕和水蝕對總侵蝕的貢獻、風力與水力的相互關系等開展了大量研究，并取得了重要成果。對先水蝕后風蝕的交錯侵蝕(水—風交錯侵蝕)研究發(fā)現(xiàn)，前期水蝕能抑制后繼風蝕，徑流沖刷通過選擇性搬運較細顆粒、地表形成物理結皮、改變微地貌等方式，減小后繼風蝕。對風蝕后水蝕的交錯侵蝕(風—水交錯侵蝕)研究發(fā)現(xiàn)，前期風蝕對后繼水蝕主要為促進作用，風蝕一方面帶走坡面土壤細顆粒，造成土壤粗化和土壤抗侵蝕能力減??；另一方面，土壤顆粒的剝離和搬運造成微地貌變化，尤其是較大風速下，形成風蝕凹痕，從而促進水蝕。前人的研究基本在徑流方向與風向相同的條件下開展，而缺乏對風向與徑流方向相反(風水反向)條件下風水交錯侵蝕的研究。尤其是風—水交錯侵蝕研究，尚未見風水反向條件下的研究報道。事實上，發(fā)生在迎風坡面徑流方向與風向相反的風水交錯侵蝕與徑流方向與風向相同的交錯侵蝕過程一樣具有代表性。因此，在徑流方向與風向相反的條件下研究風水交錯侵蝕過程與機理對系統(tǒng)的理解風水交錯侵蝕機制具有重要意義。

黃土高原風蝕水蝕交錯帶是受風水交錯侵蝕影響的典型區(qū)域，本研究選擇該區(qū)域的典型土壤，采用人工模擬降雨與風洞試驗相結合的方法，研究風向與徑流方向相反條件下的風—水交錯侵蝕特征，通過與僅水蝕下的侵蝕特征對比，分析前期風蝕對后繼水蝕侵蝕的影響，量化風蝕和水蝕對交錯侵蝕的貢獻，為明確風水交錯侵蝕中風蝕與水蝕的相互影響研究奠定基礎。

1 材料與方法

本研究采用人工模擬降雨和風洞試驗相結合的方法，以黃土高原水蝕風蝕交錯帶典型的砂質壤土(砂黃土)為例，開展風—水交錯侵蝕研究。砂質壤土采自神木六道溝流域，質地均一，砂粒、粉粒和黏粒含量分別為53.0%，42.1%和4.9%，平均粒徑34.8 μm，分選系數(shù)1.84Φ(34.8 μm)。砂質壤土在去除草根、礫石等雜物并過1 cm篩后，分層填裝于鋼制土槽(1.2 m×1.0 m×0.1 m)。土槽填裝時，首先在底部鋪設1層孔徑較小的透水紗布，再均勻平鋪1層細沙(1 cm厚)，風干的砂質壤土(含水量約為1.5%)按照容重1.35 g/cm按3 cm分層(每層3 cm)填裝，每層填裝完成后打毛表面以消除填裝土層間的分層。

人工模擬降雨和風洞試驗在黃土高原土壤侵蝕與旱地農業(yè)國家重點實驗室人工降雨大廳與風洞實驗室開展。降雨試驗使用側噴式降雨裝置，降雨高度為16 m，降雨強度在40～260 mm/h連續(xù)可調，降雨均勻度>80%，最大持續(xù)降雨時間12 h。風洞設備為室內直流吹氣式風洞(圖1)，橫截面積為1.0 m×1.2 m，由動力段、調節(jié)段、整流段、試驗段和集沙段組成，風速在2～15 m/s范圍內連續(xù)可調。吹風前進行風速率定，達到設計風速(±0.2 m/s)。

圖1 風洞結構示意

試驗時，先進行風洞試驗后再進行人工模擬降雨試驗。風洞試驗在9，12，15 m/s風速下分別持續(xù)吹蝕15 min，吹蝕前后均使用地秤(精度為1.0 g，量程200 kg)稱量土槽的總重，用于土壤風蝕量計算。對不同風速下完成吹蝕試驗的土槽，在不同雨強(60，90 mm/h)和坡度(5°，10°和15°)條件下開展人工模擬降雨試驗。所有試驗從產(chǎn)流開始計時，試驗時長為60 min。產(chǎn)流過程中，按2 min間隔收集徑流泥沙過程樣品。降雨結束后稱量徑流泥沙樣品，并在105 ℃恒重并稱量(精度0.01 g)泥沙樣品，計算降雨過程的侵蝕產(chǎn)沙量。同時，以只進行降雨而無前期吹蝕的試驗作為對照。每組試驗重復2次，當2次重復結果存在明顯差異的試驗，進行第3次重復。

依據(jù)徑流量，第min的徑流強度(g/(min·m))和整場降雨平均徑流強度(g/m)計算公式分別為：

(1)

(2)

式中：為第min降雨收集到的徑流量(g)；為徑流收集時間(min)，本研究中取值為2，4，6，…，30；為降雨坡面(土槽)面積(m)。

依據(jù)泥沙量，第min的水蝕速率(g/(min·m))和整場降雨平均水蝕速率(g/m)計算公式分別為：

(3)

(4)

式中：為第min降雨收集到的泥沙量(g)。

為量化前期風蝕對后繼水蝕速率的影響，將受前期風蝕影響的水蝕速率或徑流強度相對僅水蝕試驗條件下的變化率定義為風蝕對水蝕速率/徑流強度的影響率()，即：

(5)

式中：、分別為受前期風蝕影響的水蝕速率和徑流強度(g/(min·m))；、分別為僅水蝕條件下的水蝕速率和徑流強度(g/(min·m))。

2 結果與分析

2.1 前期風蝕對后繼水蝕徑流過程的影響

前期風蝕改變了后繼水蝕的徑流強度，但未明顯改變徑流過程的總體變化趨勢。無論是否受前期風蝕的影響，后繼徑流強度均呈現(xiàn)隨降雨歷時增加逐漸增大而后逐漸趨于穩(wěn)定的變化特征。經(jīng)前期不同風速吹蝕后，后繼水蝕的產(chǎn)流時間較僅水蝕提前了7～15 min(表1)。前期風蝕導致后繼水蝕徑流強度增大，尤其是12，15 m/s風速吹蝕后，后繼水蝕徑流強度在60 mm/h雨強下較僅水蝕分別增大191.1～340.1，280.0～695.9 g/(min·m)，在90 mm/h雨強下較僅水蝕分別增大161.5～672.6，323.4～942.3 g/(min·m)(圖2)。此外，前期風蝕對后繼水蝕徑流強度的影響在小雨強下表現(xiàn)更明顯。60 mm/h雨強下風蝕對水蝕徑流強度的影響率較90 mm/h雨強下均略增大。在相同風速吹蝕后，坡度從5°增大到15°時，60 mm/h雨強下風蝕對水蝕徑流強度的影響率較90 mm/h雨強分別增大了0～0.2，0.1～0.2，0.2～0.5倍。

表1 受前期風蝕影響的后繼降雨產(chǎn)流時間

2.2 前期風蝕對后繼水蝕速率變化的影響

總體上，前期風蝕導致后繼水蝕速率大多較僅水蝕顯著增大(<0.05)，且其增加程度隨風速的增大而增大。在9，12，15 m/s風速下，前期風蝕導致后繼水蝕速率較相同條件的僅水蝕試驗分別增大了0～1.8，0.6～3.0，1.0～4.6倍。此外，較小雨強(60 mm/h)下，前期風蝕對后繼水蝕速率的影響表現(xiàn)得更明顯(圖3)。60 mm/h雨強下風蝕對水蝕速率的影響率較90 mm/h雨強下均增大。當坡度從5°增大到15°時，在9，12，15 m/s風速吹蝕后，60 mm/h雨強下風蝕對水蝕速率的影響率較90 mm/h雨強分別增大了0.2～1.8，1.1～1.9，0.1～0.2倍?？梢?，在不同條件的風—水交錯侵蝕過程中，水蝕和風蝕的交互效應受試驗條件的影響，可能存在交互效應發(fā)生改變的臨界條件。

若以產(chǎn)沙時刻起，將水蝕過程分為前期(20 min)、中期(20～40 min)和后期(40～60 min)3個階段，隨著前期風蝕風速增大，后繼水蝕速率從前期到中期都呈現(xiàn)增大趨勢，中期到后期的變化相對復雜。從中期到后期，60 mm/h雨強下水蝕速率均呈增長趨勢；90 mm/h雨強下，水蝕速率在小坡度(5°)時均呈增長趨勢，而在較大坡度，前期風速較大時，后期水蝕速率較中期減小。與大雨強下，尤其是坡度較大時，坡面在降雨過程中易形成跌坎，導致產(chǎn)沙前期水蝕速率凸增有關(圖4)。對比前期到中期與中期到后期水蝕速率的變化速度發(fā)現(xiàn)，除12 m/s風速吹蝕后在60 mm/h雨強，5°條件下降雨的試驗外，無論是否受前期風蝕的影響，當坡度從5°增加到10°和15°時，水蝕速率的變化速度在不同雨強下均表現(xiàn)出先增長后減小的趨勢。

2.3 水蝕對風-水交錯侵蝕的貢獻

整體上，水蝕對總侵蝕的貢獻隨試驗條件的變化而明顯改變。60 mm/h雨強下風蝕貢獻占主導。90 mm/h雨強下大坡度(10°和15°)以水蝕為主，小坡度(5°)以風蝕為主。此外，交錯侵蝕中，隨前期風蝕風速增大，水蝕對交錯侵蝕的貢獻明顯減小。90 mm/h雨強下，當前期風蝕風速從9 m/s增大到15 m/s時，水蝕對交錯侵蝕的貢獻減小，比例最高達到21%，27%和28%；在60 mm/h雨強下，水蝕對交錯侵蝕的貢獻減小，比例最高達到4%，12%和43%(圖5)。證明了水蝕和風蝕的交互效應與試驗條件的相關性。

圖2 前期風蝕對后繼水蝕徑流過程的影響

注：圖柱上方不同小寫字母表示相同降雨強度下不同風速處理間差異顯著(P<0.05)。

水蝕在交錯侵蝕的比例均高于僅水蝕在僅風蝕與僅水蝕的侵蝕量之和中的比例(即：水蝕/風水交錯總侵蝕>僅水蝕/僅水蝕+僅風蝕)。與僅水蝕在僅水蝕與僅風蝕侵蝕量之和的比例相比，在9，12，15 m/s吹蝕后，60 mm/h雨強下，水蝕貢獻率較單一風蝕和水蝕過程總和中的水蝕比例分別增大了3%～8%，1%～11%，0%～11%，90 mm/h雨強下，水蝕貢獻率較單一風蝕和水蝕過程總和中的水蝕比例分別增大了1%～2%，6%～12%，9%～30%。可見，風速較大時(12，15 m/s)，前期風蝕對地表的改變對后期水蝕速率增大有明顯影響。

3 討 論

風—水交錯侵蝕中，受前期風蝕影響，后繼水蝕的徑流強度增大、水蝕加劇，這與前人在風水同向條件下的研究結果一致。已有研究發(fā)現(xiàn)，受風蝕影響的水蝕過程中產(chǎn)流時間卻相對僅水蝕試驗不同程度地滯后，本研究發(fā)現(xiàn)，前期風蝕影響導致產(chǎn)流時間相對僅水蝕試驗提前，可能與風蝕對微地貌改變程度、地表粗糙度和受試土壤性質的差異有關。首先，前期風蝕改變了地表微形態(tài)，在較大風速(12，15 m/s)下，前期風蝕后地表產(chǎn)生有固定走向的明顯風蝕凹痕，為徑流匯集提供了有利的地形條件，這是導致后繼徑流強度增大、水蝕加劇的重要原因(圖6)。前人研究中，風蝕后地表并未出現(xiàn)有固定走向的明顯風蝕凹痕，僅有少量零星的凹痕和條紋，未連通的微起伏對徑流的形成有一定的阻礙作用，延緩了徑流的形成。其次，風蝕導致地表粗化，增大土壤可蝕性，導致徑流強度增大、水蝕加劇。地表粗化對徑流的影響有兩面性，一方面可能因地表粗糙度增大而延滯徑流形成的時間；另一方面，風蝕過程中部分較細顆粒阻塞粗顆粒之間的空隙，影響水分入滲，導致產(chǎn)流時間提前。本研究的砂質壤土土質疏松、持水性差、入滲快，在降雨過程中初始產(chǎn)流時間較長，而經(jīng)吹蝕后，表層(0—1 cm)土壤細顆粒(黏粒和粉粒)含量減少，砂粒含量增多，風蝕導致表層土壤明顯粗化(表2)。在水蝕過程中產(chǎn)流時間提前，表明較細顆粒阻塞粗顆粒之間的空隙對水分入滲有明顯影響，經(jīng)前期風蝕作用后，后繼水蝕的入滲率減小(表3)?？紤]到不同的土壤在顆粒組成、團聚體結構、水分入滲等方面均有差異，風蝕過程中坡面微地形、土壤顆粒、水分入滲等對吹蝕的響應不同，因此，要量化上述因素的影響程度，還需進一步研究。

圖4 前期風蝕對后繼水蝕過程的影響

Zhang等基于Cs-137示蹤技術對結果顯示，神木六道溝流域粉質砂壤土迎風坡和背風坡面的風水交錯侵蝕(包含風—水和水—風侵蝕)以水蝕為主，徑流方向與風向相反的迎風坡(N坡和NW坡)，平均水蝕速率約為0.03 g/(min·m)，占迎風坡總侵蝕量的81.4%；而本試驗迎風坡的平均水蝕速率達到189.4 g/(min·m)，遠大于自然坡面多年平均侵蝕速率。此外，在徑流方向與風向相同的背風坡面(SE坡和S坡)，Zhang等觀測到風積(平均沉積速率<0.01 g/(min·m)，而本研究結果均為侵蝕，這與研究開展的環(huán)境、手段、尺度等直接相關。室內交錯侵蝕模擬的試驗土槽一般尺寸較小(如本研究中的1.2 m×1.0 m×0.1 m)，侵蝕難以達到沉積發(fā)生的坡長條件，在野外環(huán)境下，坡面較長且上風向坡面會對下風向坡面的侵蝕或沉積造成影響。

此外，野外條件下坡面地形、植被覆蓋、土地利用等因素也會影響風水交錯侵蝕結果。在未來的研究中可以通過考慮坡形、植被等因素，研究風水交錯侵蝕中的侵蝕和沉積特征，從而充分認識風水交錯過程。

圖5 交錯侵蝕中水蝕對總侵蝕的貢獻

圖6 風蝕試驗后的地表形態(tài)相對于風蝕前的變化

表2 風蝕前后表層(0-1 cm)砂質壤土機械組成

表3 風蝕前后降雨入滲率變化

鄧鑫欣等在風向與徑流方向相反條件下，采用相同的土壤類型在相同試驗條件(風速、降雨強度和坡度)開展的水—風交錯侵蝕得出，前期水蝕對后繼風蝕有明顯的抑制作用，本研究的風—水交錯侵蝕中，前期風蝕對后繼水蝕有明顯的促進作用。在相同風速下，水蝕在風—水交錯總侵蝕中的貢獻隨雨強的增大逐漸增大，而水—風交錯總侵蝕一直以水蝕為主(約78%的試驗組中水蝕的貢獻>50%)。相同試驗下，風—水侵蝕的總侵蝕量均較水—風侵蝕總侵蝕量大，在9，12和15 m/s風速下，風—水交錯侵蝕速率比水—風交錯侵蝕速率分別增大了129.4～146.4，640.4～834.4，1 884.3～2 239.6 g/(min·m)(圖7)?？梢?，風力與水力疊加過程不同的交錯侵蝕中，風蝕和水蝕的相互作用(促進或抑制)、風蝕和水蝕在總侵蝕中的貢獻、總侵蝕量均存在差異。

在未來的研究中，要明確風蝕交錯侵蝕過程，揭示風蝕和水蝕的相互影響，首先，必須系統(tǒng)開展不同過程(風—水、水—風)的風水交錯侵蝕，甚至應該考慮多次疊加過程的交錯侵蝕過程(如風—水—風—水、水—風—水—風)，為揭示風蝕和水蝕的互作奠定基礎；其次，本研究發(fā)現(xiàn)，不同試驗條件下，風蝕和水蝕的交互效應不同，加之不同的試驗中變量類別及各變量的梯度較多，且不同變量的數(shù)值、單位等均存在差異，直接采用變量進行交互效應分析存在巨大的難度。因此，未來的研究中應考慮采用統(tǒng)一的變量表達方式(如侵蝕力、侵蝕力的能量)量化復合侵蝕過程中的交互效應。

圖7 風水和水風交錯侵蝕的總侵蝕速率

4 結 論

通過室內風洞和模擬降雨試驗，研究了砂質壤土坡面在風向和徑流方向相反條件下的風—水交錯侵蝕特征，結果顯示，與僅水蝕試驗相比，前期風蝕導致后繼水蝕的產(chǎn)流時間提前、徑流強度增大、產(chǎn)沙速率增大，水蝕在交錯侵蝕中的貢獻較水蝕在僅風蝕與僅水蝕侵蝕量之和的比例增大。前期風蝕對后繼水蝕存在明顯的促進作用，前期風蝕導致后繼水蝕增大0～4.6倍。不同試驗條件下，前期風蝕對后繼水蝕速率變化的影響程度隨試驗條件的改變發(fā)生差異變化，表明水蝕和風蝕的交互效應受試驗條件的影響。風向與徑流方向相反條件下，前期風蝕明顯促進后繼水蝕，而前期水蝕明顯抑制后繼風蝕。與水—風交錯侵蝕相比，風—水交錯侵蝕的總侵蝕量更大，水蝕在總侵蝕中的比例變化更復雜。未來研究中應系統(tǒng)研究不同交錯侵蝕過程、多次交錯過程疊加、地形影響的交錯侵蝕，在此基礎上，考慮采用統(tǒng)一的變量表達方式量化風蝕和水蝕的交互效應，為揭示風蝕和水蝕的互作機理奠定基礎。