亞喀斯特林地野外模擬降雨條件下地表徑流特征

羅雅雪， 張思琪， 顏 紅， 岳彩雯， 韋小茶， 周秋文

(貴州師范大學(xué) 地理與環(huán)境科學(xué)學(xué)院， 貴州 貴陽 550001)

1 研究背景

地表徑流是地表土壤侵蝕、面源污染以及養(yǎng)分流失的動力之一[1-4]，具有可以侵蝕地表土壤的能量，從而改變徑流區(qū)、流域的地貌特征[5-6]。亞喀斯特地貌是針對貴州等特殊的喀斯特生態(tài)環(huán)境，提出的一種介于典型喀斯特地貌以及非喀斯特地貌之間的地貌形態(tài)[7]。該地貌類型以丘陵谷地、淺中丘、淺切割低中山為主，地形起伏較緩。“亞喀斯特”地區(qū)喀斯特發(fā)育程度不高，土層較厚，由于其地質(zhì)環(huán)境的特殊性，生態(tài)系統(tǒng)較為脆弱，穩(wěn)定性差，同等條件下土壤侵蝕量較高，對外界干擾很敏感[8]。因此研究亞喀斯特地區(qū)地表徑流特征具有重要意義。

前人對不同地區(qū)的不同下墊面條件的地表徑流等方面做了大量的研究，并取得頗多的研究成果。李文華等[9]研究總結(jié)，地表徑流的調(diào)節(jié)受多種因素的影響，包括氣候條件、植被特征、土壤條件、地質(zhì)、地形特征以及人類活動影響等。李耀明等[10]在縉云山地區(qū)研究降雨因子對地表徑流的影響，研究表明，對地表徑流起著制約性的影響因素有降雨強度、降雨歷時和降雨量等綜合因子。張洪江等[11]在相同地區(qū)做不同植被類型對地表徑流系數(shù)的影響分析研究，得出植被類型對地表徑流系數(shù)有影響，從而影響地表徑流。唐小燕[12]在錢江源進行典型森林類型地表徑流研究。肖登攀等[13]對不同地表類型降雨入滲產(chǎn)流規(guī)律進行研究，結(jié)果顯示地表類型不同，地表產(chǎn)流特征也不同。

夏子翔等[14]在喀斯特地區(qū)進行流域巖性對地表徑流形態(tài)的影響研究，研究得出巖性對地表徑流形態(tài)有影響。涂成龍等[15]在典型喀斯特流域進行地表產(chǎn)流輸出特征研究。張喜等[16]選取黔中喀斯特地區(qū)的不同山地，對不同類型的森林做地表徑流相關(guān)研究，研究證明森林類型對地表徑流有影響。胡奕等[17]以我國西南喀斯特地區(qū)為研究區(qū)進行地表徑流相關(guān)研究，研究得出巖石裸露率、降雨強度、坡度和地下孔隙度均對地表徑流造成影響的結(jié)論。伏文兵等[18]對喀斯特耕地進行研究，結(jié)果表明地表徑流與降雨強度、坡度的變化呈極顯著正相關(guān)，與地下孔隙度的變化呈顯著負相關(guān)。綜上所述，已有的研究主要集中在典型的喀斯特地區(qū)或非喀斯特地區(qū)，針對亞喀斯特地區(qū)的研究很少。

鑒于亞喀斯特地貌是一種重要的地貌類型單元，在中國西南地區(qū)占有重要的地位[19]。因此，本研究采用人工模擬降雨實驗的方法，對亞喀斯特地貌地區(qū)的地表徑流特征進行分析研究，籍此揭示亞喀斯特地區(qū)在不同降雨條件下的地表徑流變化規(guī)律，為亞喀斯特地區(qū)的生態(tài)環(huán)境建設(shè)和可持續(xù)發(fā)展提供理論基礎(chǔ)和科學(xué)依據(jù)。

2 研究區(qū)概況

研究樣地位于貴州省貴陽市花溪區(qū)桐木嶺，地理位置為106°40′47″E，26°22′15.85″N。花溪區(qū)地處烏江與珠江分水嶺，地表景觀破碎，以山地和丘陵為主。屬于亞熱帶季風(fēng)濕潤氣候，全年平均氣溫為14.9℃，年極端最高氣溫為35.1℃，年極端最低溫度為-7.3℃。年降雨量1 178.3 mm，蒸發(fā)量738 mm，雨量充沛。研究樣地的巖性為中三疊統(tǒng)花溪組白云巖(含雜質(zhì))，主要土壤類型為黃壤，植被優(yōu)勢種為馬尾松。該地區(qū)成土母質(zhì)相對豐富，土層較厚，地貌和景觀均符合亞喀斯特地貌類型特征。

3 研究方法

3.1 徑流小區(qū)設(shè)置

在選取的亞喀斯特樣地內(nèi)，共設(shè)置9個面積為60 cm×100 cm的臨時徑流小區(qū)。樣地內(nèi)的9個徑流小區(qū)盡量連續(xù)分布，目的在于使徑流小區(qū)下墊面條件均一，以減少誤差的產(chǎn)生。徑流小區(qū)的組成材料包括兩塊帶有把手的長度為1 m的鋼板、一塊長為60 cm的鋼板，所有鋼板寬度均為20 cm。另有一塊60 cm長的V型不銹鋼制集流槽，集流槽一端深度為5 cm，另一端深度為10 cm，將10 cm深的一端作為出水口。設(shè)置徑流小區(qū)時，先將兩塊長為1 m的鋼板平行楔入土壤中(圖1(a))，兩塊鋼板的間隔為60 cm。再將60 cm長的鋼板垂直插入兩塊鋼板上頂端處，使其緊密接觸形成封閉的長方形。將鋼板插入土層內(nèi)10 cm左右的深度，然后在兩塊鋼板的下接口處挖一條溝槽，將集流槽放入該溝槽內(nèi)。為保證產(chǎn)生的地表徑流量都能被采集，集流槽上邊略低于地面，集流槽靠近徑流小區(qū)一側(cè)楔入土壤中。在集流槽的出水口處，挖一個較深的坑，用于放置容器，收集集流槽流出的水，完整的徑流小區(qū)如圖1(b)所示。為避免雨水淋，V型鋼板和水樣收集容器的上方用防水布蓋好，且不能覆蓋徑流小區(qū)的承雨范圍。此外，用防水布將試驗樣地兩側(cè)的其他樣地蓋好，避免其他樣地在未開展實驗前被雨水淋濕，影響地表徑流數(shù)據(jù)的收集。

3.2 人工模擬降雨實驗設(shè)置

本研究采用便攜式人工模擬降雨裝置開展實驗。該裝置在實驗前已經(jīng)過檢驗，雨滴的霧化效果良好，能較好地模擬自然降雨。人工降雨裝置安裝在徑流小區(qū)旁，地勢平坦且無其他植物影響噴灑過程。為使雨滴速度接近天然降雨，噴頭固定在距地面2.3 m處，模擬降雨過程如圖2所示。

圖1 徑流小區(qū)設(shè)置

圖2 人工模擬降雨實驗過程

本次實驗的唯一變量為雨強，共設(shè)置3個不同的雨強，分別為30 mm/h、75 mm/h、120 mm/h。在同一雨強下設(shè)置3個重復(fù)實驗，最終數(shù)據(jù)為3組實驗數(shù)據(jù)平均值，故樣地內(nèi)共需進行9場人工降雨。降雨歷時設(shè)置為30 min，噴頭開始噴灑雨霧時計時，待集流槽出口開始有連續(xù)水流時，記錄產(chǎn)流時間，并以3 min為一時間段，每隔3 min采集一次徑流樣品。每次收集水樣后即刻用量筒量取并記錄讀數(shù)，分別用瓶子裝好，并在瓶子上注明采樣時間及水樣的體積。實驗進行30 min時立即關(guān)閉噴灑系統(tǒng)，一塊樣地的實驗便完成。立即在與該樣地相鄰的另一樣地中做相同的實驗。

3.3 相關(guān)計算公式

產(chǎn)沙量與徑流系數(shù)的計算公式如下：

cs=(m-ρV)/V

(1)

c=w/(i·t·s)

(2)

式中：cs為產(chǎn)沙量，g/mL；m為含沙的水的質(zhì)量，g；ρ為水的密度g/cm3；V為含沙的水的體積，mL；c為徑流系數(shù)；w為徑流量，mL；i為降雨強度；mm/h；t為降雨時間，h；s為降雨面積，cm2。

4 結(jié)果與分析

4.1 不同雨強和降雨歷時的地表徑流量變化

降雨強度的不同，導(dǎo)致相同時間段內(nèi)產(chǎn)生的地表徑流量有所不同。以3min為一時間段進行分析，徑流量的變化如圖3所示。30mm/h雨強條件下，模擬降雨開始3min時的地表徑流量為最小值，只有28mL，一直到15min時產(chǎn)生的地表徑流量始終保持遞增的趨勢。15min時產(chǎn)生的地表徑流量為128mL，而在18min時，地表徑流量減少至112mL，在19min到24min內(nèi)，地表徑流量繼續(xù)呈現(xiàn)遞增狀態(tài)。地表徑流量在24min時達到最大值為137mL，24min后先減少后增加，地表徑流量最大值與最小值之間的差值為109mL。

圖3 不同雨強和降雨歷時條件下地表徑流量變化

75mm/h雨強條件下，模擬降雨開始3min時的地表徑流量為81mL，地表徑流量在24min時達到最大值312mL，整個產(chǎn)流過程中地表徑流量變化不穩(wěn)定，降雨結(jié)束時產(chǎn)生的地表徑流量為210mL，地表徑流量最大值與最小值之間的差值為231mL。當(dāng)雨強為120mm/h時，模擬降雨開始3min時的地表徑流量為157mL，地表徑流量在21min時達到最大值為497mL，在降雨結(jié)束時地表徑流量減至400mL，地表徑流量最大值與最小值之間的差值為340mL。

從圖3可看出，整個降雨過程的地表徑流量都表現(xiàn)為隨著雨強的增大而增加。在降雨開始3min時，雨強為120mm/h產(chǎn)生的地表徑流量是雨強為30mm/h時的5.61倍，是雨強為75mm/h時的1.94倍。當(dāng)降雨結(jié)束時，雨強為120mm/h產(chǎn)生的地表徑流量是雨強為30mm/h時的3.13倍，是雨強為75mm/h時的1.90倍。

4.2 不同雨強和降雨歷時條件下徑流系數(shù)變化

徑流系數(shù)是一定匯水面積內(nèi)總徑流量與降雨量的比值，能綜合反映流域內(nèi)自然地理要素對徑流的影響。由圖4可看出，降雨強度不同，導(dǎo)致不同降雨歷時條件下的徑流系數(shù)也不同。其中，雨強為30和75mm/h時的徑流系數(shù)變化趨勢基本相同，但是，當(dāng)降雨歷時為3和6min時，雨強75mm/h較雨強30mm/h的徑流系數(shù)大，而在其他降雨歷時中均表現(xiàn)為雨強75mm/h的徑流系數(shù)比雨強30mm/h的小。當(dāng)雨強為120mm/h時，徑流系數(shù)在降雨歷時為3和21min時均大于其他兩個雨強，除這兩個降雨歷時以外，雨強為120mm/h的徑流系數(shù)均小于雨強75mm/h和雨強30mm/h的徑流系數(shù)。不同雨強和降雨歷時下的徑流系數(shù)總體表現(xiàn)為雨強越大，徑流系數(shù)越小，但在降雨歷時為3和21min時出現(xiàn)不同的規(guī)律，產(chǎn)生此現(xiàn)象的原因是徑流收集時操作不當(dāng)引起的，最終導(dǎo)致徑流系數(shù)偏大。

圖4 不同雨強和降雨歷時條件下的徑流系數(shù)

降雨過程中，不同雨強的最大徑流系數(shù)不同。雨強為30mm/h時為0.0407，雨強為75mm/h時為0.0420，當(dāng)雨強增至120mm/h時，最大徑流系數(shù)為0.0435。以最大徑流系數(shù)來看，表現(xiàn)為雨強越大，徑流系數(shù)越大。但從圖4看出，徑流系數(shù)整體變化為雨強越大，徑流系數(shù)越小，分析原因為雨強越大，產(chǎn)生的地表徑流量越大，但不同雨強下不同降雨歷時條件的地表徑流增加量有所不同。

4.3 不同雨強和降雨歷時條件下產(chǎn)沙量變化

雨強大小的不同，對土壤的侵蝕能力可通過產(chǎn)沙量大小顯示。如圖5所示，雨強為30mm/h時的不同時段產(chǎn)沙量變化較大，降雨開始3min時產(chǎn)沙量為最大值0.0088g/mL，最小值在24min出現(xiàn)，產(chǎn)沙量為0.0004g/mL，最大值與最小值間的差值為0.0084g/mL，降雨30min時的產(chǎn)沙量為0.0017g/mL。雨強為75mm/h時，產(chǎn)沙量變化幅度較雨強為30mm/h時變化小，降雨開始3和9min時的產(chǎn)沙量相等且為最大值，數(shù)值為0.0056g/mL。當(dāng)降雨為18min時產(chǎn)沙量為最小值，降至0.0002g/mL，差值為0.0054，降雨30min時的產(chǎn)沙量為0.0007g/mL。當(dāng)雨強為120mm/h時，降雨6min時產(chǎn)沙量達最大值，產(chǎn)沙量為0.0075g/mL，27min時為最小值，產(chǎn)沙量為0.0003g/mL，差值為0.0072，降雨30min時的產(chǎn)沙量為0.004g/mL。雨強為30mm/h時，產(chǎn)沙量的變化趨勢、幅度均較大。雨強為75mm/h時的產(chǎn)沙量較雨強為30mm/h時變化幅度較小，折線起伏較小。雨強為120mm/h時，產(chǎn)沙量在降雨3～6min間的變化最大，從3min時的0.0006g/mL增至0.0075g/mL，又在9min時降至0.0012g/mL，9min以后的變化趨勢減緩。表明雨強為120mm/h時，不同時段的產(chǎn)沙量差值較小，產(chǎn)沙量相差不大。

總體而言，侵蝕產(chǎn)沙隨降雨歷時變化的幅度較地表徑流大，產(chǎn)沙量變化曲線呈現(xiàn)明顯的不規(guī)則波動。其原因是，土粒崩解具有隨機性和偶發(fā)性，加上泥沙的運移在短距離、短時間內(nèi)具有非勻速特征，因此在集水區(qū)面積較小情況下，產(chǎn)沙量呈現(xiàn)出明顯的不規(guī)則波動。

圖5 不同雨強和降雨歷時條件下的產(chǎn)沙量變化

4.4 不同雨強條件下產(chǎn)流時間變化

如表1所示，不同雨強條件下產(chǎn)流時間不同。雨強為30mm/h時，產(chǎn)流時間為2.37min，當(dāng)雨強為75mm/h時，產(chǎn)流時間降至1.29min，當(dāng)雨強增至120mm/h時，產(chǎn)流時間又增為1.51min。雨強為75mm/h時的產(chǎn)流時間較雨強為30mm/h時的產(chǎn)流時間少1.08min，較雨強為120mm/h時少0.22min，差距較雨強為30mm/h時小的多。由表1可知，雨強為75mm/h時產(chǎn)流所需時間最短，即產(chǎn)流最快，但是產(chǎn)流時間并未隨著降雨強度的增加而保持遞減的趨勢，雨強75mm/h與雨強120mm/h的產(chǎn)流時間差距不大，推測可能是由于實驗誤差導(dǎo)致的。一方面，設(shè)置徑流小區(qū)時集流槽與土壤切面接觸不緊密，導(dǎo)致降雨初期產(chǎn)生的少量地表徑流從集流槽與土壤切面間的縫隙漏掉；另一方面，可能由于徑流小區(qū)下墊面狀況導(dǎo)致，設(shè)置雨強120mm/h的3個徑流小區(qū)土壤表層均較疏松，或枯落物層較厚，有利于水分蓄滯和入滲，導(dǎo)致產(chǎn)流時間推遲。

4.5 不同雨強條件下累積徑流量與徑流系數(shù)變化

降雨強度的大小影響地表徑流量，導(dǎo)致降雨結(jié)束時不同雨強的累積地表徑流量不同。如表1所示，降雨結(jié)束時，雨強為30mm/h的累積地表徑流量為1 075mm，雨強為70mm/h時，累積地表徑流量為2 351mm，降雨強度增至120mm/h時，累積地表徑流量增加到3 526mm。總體而言，歷時30min的降雨結(jié)束時，雨強為120mm/h的累積地表徑流量是雨強為30mm/h時累積地表徑流量的3.28倍，是雨強為75mm/h時累積地表徑流量的1.50倍，顯示出雨強越大，產(chǎn)生的地表徑流量越多的規(guī)律。

表1 不同雨強的產(chǎn)流產(chǎn)沙情況

由表1可知，歷時30min的降雨結(jié)束時，雨強為30mm/h的徑流系數(shù)為0.119，雨強為75mm/h的徑流系數(shù)為0.104，雨強為120mm/h的徑流系數(shù)為0.098。降雨結(jié)束時徑流系數(shù)隨著雨強的增加呈現(xiàn)遞減趨勢，即為雨強越大，徑流系數(shù)越小。

5 結(jié) 論

本文以貴州省亞喀斯特地區(qū)為例，基于人工模擬降雨的方法，分析不同雨強下地表徑流的產(chǎn)流時間、累積地表徑流量、徑流系數(shù)、產(chǎn)流速率和產(chǎn)沙量，得出以下結(jié)論：

(1)在降雨實驗過程中，總體表現(xiàn)為地表徑流量隨雨強的增大而增加。

(2)歷時30min的降雨過程中，徑流系數(shù)整體表現(xiàn)為雨強越大，徑流系數(shù)越小。

(3)侵蝕產(chǎn)沙量隨降雨歷時變化的幅度較地表徑流大，產(chǎn)沙量變化曲線呈現(xiàn)明顯的不規(guī)則波動。

(4)產(chǎn)流時間表現(xiàn)為雨強30mm/h時最長，120mm/h時次之，雨強75mm/h時所需時間最短。

(5)累積地表徑流量隨著雨強的增加而增加，此時的徑流系數(shù)隨著雨強的增加而逐漸減小。

本項研究的對象為亞喀斯特林地，相對已有研究而言，對喀斯特、亞喀斯特地區(qū)地表徑流和侵蝕產(chǎn)沙的研究是一項有益的補充。此外，研究中以雨強作為唯一變量，未分析其他變量對地表徑流的影響，在今后的研究中可以增加坡度、土壤類型、森林類型等因素從多個角度分析亞喀斯特地貌的徑流產(chǎn)沙效應(yīng)，為喀斯特、亞喀斯特地區(qū)水土保持和生態(tài)保護建設(shè)等提供一定的理論參考。

[1] 魏 強,張秋良.地表徑流與土壤侵蝕研究進展[J].中國水土保持,2008(9):30-33.

[2] 李恒鵬,黃文鈺,楊桂山,等.太湖上游典型城鎮(zhèn)地表徑流面源污染特征[J].農(nóng)業(yè)環(huán)境科學(xué)學(xué)報,2006,25(6):l598-1602.

[3] 黃河仙,謝小立,王凱榮,等.不同覆被下紅壤坡地地表徑流及其養(yǎng)分流失特征[J].生態(tài)環(huán)境,2008,17(4):1645-1649.[4] 孫 閣.林地地表徑流的研究[J].水土保持學(xué)報,1989,3(2):52-56.

[5] 張慧芳.東山流域水土保持措施對地表徑流和土壤侵蝕的影響[J].山西水利,2012，28(11):22-41.

[6] 李秋芳,王克勤,王帥兵,等.不同治理措施在紅壤坡耕地的水土保持效益[J].水土保持通報,2012,32(6):196-200.[7] 胡 鋒,安裕倫,許 璟.“亞喀斯特”概念與景觀特征的初步探討——以貴州為例[J].地理研究,2015,34(8):1569-1580.

[8] 胡 鋒,安裕倫,趙海兵.基于土地利用轉(zhuǎn)型視角的“亞喀斯特”區(qū)域生態(tài)環(huán)境效應(yīng)研究——以黔中部分地區(qū)為例[J].地球與環(huán)境,2016,44(4):447-454.

[9] 李文華,何永濤,楊麗韞.森林對徑流影響研究的回顧與展望[J].自然資源學(xué)報,2001,16(5):398-406.

[10] 李耀明,王玉杰,儲小院,等.降雨因子對縉云山地區(qū)典型森林植被類型地表徑流的影響[J].水土保持研究,2009,16(4):244-249.

[11] 張洪江,孫艷紅,程 云,等.重慶縉云山不同植被類型對地表徑流系數(shù)的影響[J].水土保持學(xué)報,2006,20(6):11-13+45.

[12] 唐小燕.錢江源典型森林類型地表徑流和土壤侵蝕特征研究[D].杭州：浙江農(nóng)林大學(xué)，2012.

[13] 肖登攀,韓淑敏,楊艷敏,等.太行山低山丘陵區(qū)不同地表類型降雨入滲產(chǎn)流規(guī)律研究[J].水土保持研究,2009,16(5):35-39.

[14] 夏子翔,趙海兵,胡洪麗.喀斯特流域巖性對地表徑流形態(tài)的影響——以花溪河流域為例[J].安徽農(nóng)業(yè)科學(xué),2017,45(3):52-54.

[15] 涂成龍,陸曉輝,劉裕錄,等.典型喀斯特流域地表產(chǎn)流輸出特征[J].長江流域資源與環(huán)境,2016,25(12):1879-1885.

[16] 張 喜,薛建輝,許效天,等.黔中喀斯特山地不同森林類型的地表徑流及影響因素[J].熱帶亞熱帶植物學(xué)報,2007,15(6):527-537.

[17] 胡 奕,戴全厚,王佩將.喀斯特坡耕地產(chǎn)流特征及影響因素[J].水土保持學(xué)報,2012,26(6):46-51.

[18] 伏文兵,戴全厚,嚴友進.喀斯特坡耕地及其淺層孔(裂)隙土壤侵蝕響應(yīng)試驗研究[J].水土保持學(xué)報,2015,29(2):11-16+22.

[19] 許 璟,安裕倫,胡 鋒,等.基于植被覆蓋與生產(chǎn)力視角的亞喀斯特區(qū)域生態(tài)環(huán)境特征研究——以黔中部分地區(qū)為例[J].地理研究,2015,34(4):644-654.