許洪嬌， 胡海波， 初 磊， 燕 超， 王超騏

(南京林業(yè)大學(xué)南方現(xiàn)代林業(yè)協(xié)同創(chuàng)新中心，江蘇省水土保持與生態(tài)修復(fù)重點(diǎn)實(shí)驗(yàn)室，210037，南京)

近年來，南方紅壤區(qū)水土流失狀況雖整體好轉(zhuǎn)，但局部惡化，減少的幅度與治理面積不相稱，每年人為新增水土流失面積達(dá)到2 256 km2，水土流失治理任務(wù)仍然十分艱巨[1]。降雨是產(chǎn)生土壤侵蝕的主要自然因素[2-4]，主要包括降雨強(qiáng)度、降雨量、降雨歷時等因子?？萋湮飳邮巧稚鷳B(tài)系統(tǒng)的重要組成部分，在攔蓄降雨、保護(hù)地表土壤等方面都具有重要意義[5-6]。國內(nèi)外學(xué)者對不同森林類型枯落物的水文特征和生態(tài)效應(yīng)進(jìn)行了大量研究，取得了許多重要成果，但有關(guān)林下枯落物覆蓋程度對于土壤侵蝕的影響研究較少。張楨堯等[7]對栓皮櫟林枯落物減流減沙效應(yīng)研究發(fā)現(xiàn)，枯落物層能延長坡面徑流產(chǎn)生的時間，降低坡面流流速，同時還能有效地減少坡面的產(chǎn)流量(6.60%～8.17%)和產(chǎn)沙量(11.08%～65.08%)。張芝萍等[8]對北亞熱帶毛竹林枯落物水土保持功能進(jìn)行研究，發(fā)現(xiàn)室內(nèi)人工模擬降雨條件下，降雨強(qiáng)度為10～50 mm/h時，有枯落物覆蓋的土壤侵蝕量比裸露坡面減少35.58～256.16 g/m2。

目前，國內(nèi)外學(xué)者對于枯落物層能夠調(diào)節(jié)水文功能的研究，大多是采用浸泡法或通過室內(nèi)人工模擬降雨試驗(yàn)測定枯落物的最大持水量、持水率和持水過程[9-10]，而涉及枯落物減蝕能力以及不同枯落物覆蓋對土壤水蝕過程影響的研究較少[11-12]。筆者主要研究在不同降雨強(qiáng)度時香樟(Cinnamomumcamphora)枯落物量的減流減沙功能, 能夠揭示枯落物覆蓋的坡面侵蝕機(jī)制，加深人們對水土保持和土壤侵蝕的理解。

1 材料與方法

1.1 試驗(yàn)裝置與材料

實(shí)驗(yàn)于2018年6月在南京林業(yè)大學(xué)溧水白馬基地進(jìn)行，白馬基地位于南方紅壤區(qū)，多年平均降雨量(1967—2017年)為1 072.90 mm(南京站)，汛期為5—9月。試驗(yàn)地開闊空曠，便于安裝器材及開展試驗(yàn)。香樟為南方地區(qū)常見樹種，因此選取它作為研究對象。試驗(yàn)所用模擬降水裝置為南京南林電子科技有限公司研制的NLJY- 10型人工模擬降雨控制系統(tǒng)，本次實(shí)驗(yàn)每次設(shè)置3個重復(fù)，以減少實(shí)驗(yàn)誤差。準(zhǔn)備4個1 m×2 m×0.5 m的長方形鋼槽作為降水試驗(yàn)小區(qū)。本試驗(yàn)用土均來自于南京林業(yè)大學(xué)溧水白馬基地，以南京林業(yè)大學(xué)校園內(nèi)的香樟林枯落物(枝葉比例2：8)作為研究材料。為減少葉片表面攜帶泥沙、水分，導(dǎo)致試驗(yàn)誤差，將枯落物用清水洗凈,在65 ℃下烘干至恒質(zhì)量，分成50、150和250 g分別編號、裝袋，備用。

1.2 實(shí)驗(yàn)設(shè)計(jì)與測定

本試驗(yàn)在4個徑流小區(qū)分層填土，每10 cm一層，共分3層，每層填裝完畢充分壓實(shí)，保持土壤密度為1.34 g/cm左右。填土前，在土槽的底部平鋪17 cm細(xì)沙，保證良好的透水性。填土?xí)r，注意坡面保持平整，每層填裝完壓實(shí)后刮毛，制造一定的地表糙度，防止土層間滑動。通過不同量枯落物覆蓋設(shè)置各處理組，共設(shè)計(jì)4種枯落物覆蓋量，分別為0、50、100和250 g/m2；設(shè)計(jì)3個降雨強(qiáng)度，分別為30、60、90 mm/h，坡度均為5°，降雨均勻度超過90%。每次降雨結(jié)束后，土槽需要被靜置至土壤含水率與上一次實(shí)驗(yàn)相一致時再重新開始實(shí)驗(yàn)。為保證試驗(yàn)精度，每個處理均進(jìn)行重復(fù)試驗(yàn)取平均值分析(表1)。

表1 枯落物覆蓋量及降雨強(qiáng)度設(shè)計(jì)Tab.1 Litter coverage and rainfall intensity design

降雨前，將枯落物均勻鋪在試驗(yàn)土體上，保證厚度均勻。實(shí)驗(yàn)時間為產(chǎn)流開始后的1 h，由于試驗(yàn)初期產(chǎn)流率變化較大，因此前10 min每2 min取1次徑流泥沙樣，產(chǎn)流10 min以后每 5 min取1次徑流泥沙樣，用聚乙烯塑料桶取徑流水樣，用量筒來測量徑流的體積，然后將徑流靜置12 h，倒掉上層清液，先用烘干法測定泥沙質(zhì)量，然后稱量測得產(chǎn)沙量。試驗(yàn)小區(qū)設(shè)置4個雨量筒，以便測定實(shí)際降雨量。測定指標(biāo)包括降水量、降水強(qiáng)度、降水歷時、徑流量和產(chǎn)沙量等。

1.3 數(shù)據(jù)處理與分析

水動力學(xué)參數(shù)采用Microsoft Excel軟件進(jìn)行數(shù)據(jù)統(tǒng)計(jì)及繪圖，采用SPSS進(jìn)行顯著性分析。徑流減少率為枯落物覆蓋條件下的產(chǎn)流量相比裸坡的減少量與裸坡產(chǎn)流量的比值，計(jì)算公式[12]為：

Lp=((Sf-Sp)/Sf)×100%。

(1)

式中：Lp為徑流減少率，%；Sf為裸地產(chǎn)流量，L/(min·m2)；Sp為覆蓋坡面產(chǎn)流量，L/(min·m2)。

產(chǎn)沙減少率為枯落物覆蓋條件下坡面產(chǎn)沙率相比裸坡的減少值與裸坡產(chǎn)沙率的比值,計(jì)算公式[13]為：

Np=((Rf-Rp)/Rf)×100%。

(2)

式中：Np為產(chǎn)沙減少率，%；Rf為裸地產(chǎn)沙量，g/(min·m2)；Rp為覆蓋坡面產(chǎn)沙量，g/(min·m2)。

2 結(jié)果與分析

2.1 枯落物覆蓋下坡面產(chǎn)流過程特征

2.1.1 枯落物覆蓋的產(chǎn)流過程 從圖1～3可以看出，產(chǎn)流率均是隨著降雨強(qiáng)度的增大而增大，隨著枯落物的生物量增大而減小，產(chǎn)流率在前5 min迅速提高，在第20～25 min達(dá)到穩(wěn)定。當(dāng)降雨強(qiáng)度為30 mm/h時，達(dá)到穩(wěn)定時的裸坡徑流量為0.605 L/(min·m2)?？萋湮锔采w的生物量為50、100和250 g/m2時，達(dá)到穩(wěn)定時的產(chǎn)流率分別為0.469、0.387和0.326 L/(min·m2)。相比于裸坡，枯落物覆蓋量為50、100和250 g/m2時，達(dá)到穩(wěn)定時的產(chǎn)流率分別減少22.48%、36.03%和46.12%。當(dāng)降雨強(qiáng)度增大到60 mm/L時，徑流達(dá)到穩(wěn)定時，枯落物覆蓋量增多，減少的產(chǎn)流率依次7.70%、7.65%和15.83%。當(dāng)降雨強(qiáng)度達(dá)到90 mm/h時，產(chǎn)流率則分別減少0.70%、2.30%和12.45%。可以看出隨著降雨強(qiáng)度增大，枯落物對產(chǎn)流率的影響程度逐漸降低，降雨強(qiáng)度在30和60 mm/h時枯落物對于徑流量的影響大于90 mm/h。

圖1 降雨強(qiáng)度為30 mm/h時不同枯落物覆蓋量 產(chǎn)流率隨降雨歷時的變化Fig.1 Changes in runoff yield rate with different litter coverage under rainfall intensity of 30 mm/h

圖2 降雨強(qiáng)度為60 mm/h時不同枯落物覆蓋程度 下產(chǎn)流率隨降雨歷時的變化Fig.2 Changes in runoff yield rate with different litter coverage under rainfall intensity of 60 mm/h

圖3 降雨強(qiáng)度為90 mm/h時不同枯落物覆蓋程度下 產(chǎn)流率隨降雨歷時的變化Fig.3 Changes in runoff yield rate with different litter coverage under rainfall intensity of 90 mm/h

2.1.2 枯落物覆蓋的產(chǎn)流特征 如圖4所示，隨著降雨強(qiáng)度增大，徑流量減少率逐漸降低。在同等的降雨強(qiáng)度下，隨著枯落物量的增大，對坡面的減流作用逐漸增強(qiáng)。當(dāng)降雨強(qiáng)度為30 mm/h時，枯落物量在250 g/m2時，平均產(chǎn)流率僅為0.29 L/(min·m2)；而降雨強(qiáng)度達(dá)到90 mm/h時，產(chǎn)流率可以達(dá)到1.13 L/(min·m2)。徑流量隨著降雨強(qiáng)度的增大不斷增大，枯落物的減流作用逐漸下降。而在相同枯落物量的覆蓋下，當(dāng)降雨強(qiáng)度為30 mm/h時，枯落物對于徑流量的減少最為明顯，徑流量的減少率可以達(dá)到34.21%，當(dāng)降雨強(qiáng)度增加到90 mm/h、枯落物量為250 g/m2時，徑流量的減少率僅為18.23%。

圖4 不同降雨強(qiáng)度時枯落物量對徑流量的影響Fig.4 Effects of litter coverage on runoff reduction rate under different rainfall intensity

以降雨強(qiáng)度I和枯落物量G為自變量，覆蓋枯落物的坡面徑流量Q為因變量，運(yùn)用SPSS 18.0軟件進(jìn)行回歸分析，得到如下多元回歸方程式：

Q=0.980I-0.186G-0.977(I>30，0

(3)

式中：I為降雨強(qiáng)度，mm/h；G為枯落物量，g/m2；Q為坡面徑流量，L/(min·m2)。從式(3)可以看出，坡面徑流量與降雨強(qiáng)度和枯落物量呈線性相關(guān)關(guān)系，通過對比系數(shù)可以得知，降雨強(qiáng)度對坡面徑流量的影響高于枯落物對徑流量的影響。

2.2 枯落物覆蓋下坡面產(chǎn)沙過程特征

2.2.1 枯落物覆蓋的產(chǎn)沙過程 從圖5～7可以看出，產(chǎn)沙率-降雨歷時的變化曲線在不同量的枯落物覆蓋下相似，從坡面產(chǎn)沙開始，短時間內(nèi)產(chǎn)沙率迅速提高，在保持平穩(wěn)后下降，最終趨于穩(wěn)定。產(chǎn)沙率同時受降雨強(qiáng)度和枯落物量的制約，當(dāng)降雨強(qiáng)度在30 mm/h時，枯落物覆蓋量達(dá)到50 g/m2就可以起到良好的防蝕效果，最終產(chǎn)沙率穩(wěn)定在0.078 4 g/(min·m2)。當(dāng)枯落物覆蓋量達(dá)到250 g/m2時，產(chǎn)沙率穩(wěn)定時僅為0.019 g/(min·m2)。隨著降雨強(qiáng)度增大，枯落物對產(chǎn)沙率的影響逐漸減小。當(dāng)降雨強(qiáng)度為30 mm/h，枯落物量分別在50、100和250 g/m2時，產(chǎn)沙率分別降低59.57%、86.35%和97.58%。當(dāng)降雨強(qiáng)度上升到60 mm/h，枯落物量分別在50、100和250 g/m2時，產(chǎn)沙率分別降低51.32%、56.56%和68.10%。當(dāng)降雨強(qiáng)度達(dá)到90 mm/h，達(dá)到穩(wěn)定后的產(chǎn)沙率分別降低21.59%、42.38%和66.45%。

圖5 降雨強(qiáng)度為30 mm/h時枯落物量對產(chǎn)沙率的影響Fig.5 Effects of litter coverage on sediment yield rate under rainfall intensity of 30 mm/h

圖6 降雨強(qiáng)度為60 mm/h時枯落物量對產(chǎn)沙率的影響Fig.6 Effects of litter coverage on sediment yield rate under rainfall intensity of 60 mm/h

圖7 降雨強(qiáng)度為90 mm/h時枯落物量對產(chǎn)沙率的影響Fig.7 Effects of litter coverage on sediment yield rate under rainfall intensity of 90 mm/h

2.2.2 枯落物覆蓋的產(chǎn)沙特征 從實(shí)驗(yàn)結(jié)果可以看出，產(chǎn)沙量隨著降雨強(qiáng)度的增大而增大，枯落物對于降低產(chǎn)沙量有明顯作用。當(dāng)降雨強(qiáng)度為30 mm/h時，裸坡的總產(chǎn)沙量為69.55 g，而當(dāng)枯落物的生物量增加到250 g/m2時，總產(chǎn)沙量僅為2.273 g。在同等降雨強(qiáng)度下，產(chǎn)沙率隨著枯落物生物量的增加而降低；在枯落物的生物量相同時，產(chǎn)沙率隨著降雨強(qiáng)度的增加而不斷增加?？萋湮飳τ诋a(chǎn)沙量的減小作用隨著降雨強(qiáng)度的不斷增加而減弱，當(dāng)降雨強(qiáng)度為30 mm/h，枯落物的生物量達(dá)到250 g/m2時，相比較于裸坡，產(chǎn)沙量可降低96.73%；當(dāng)降雨強(qiáng)度為90 mm/h時，枯落物的生物量為250 g/m2，產(chǎn)沙量僅降低49.51%。

以降雨強(qiáng)度I和枯落物量G為自變量，覆蓋枯落物坡面的總產(chǎn)沙量S為因變量，運(yùn)用SPSS18.0軟件進(jìn)行回歸分析，得到如下多元回歸方程式：

S=3.969I-0.572G-50.137 (I>30，0

(4)

從式(4)可以看出，隨著枯落物量的增大，產(chǎn)沙量逐漸減小；隨著降雨強(qiáng)度的增大，產(chǎn)沙量增大。通過斜率比較發(fā)現(xiàn)，降雨強(qiáng)度對產(chǎn)沙量的影響明顯高于枯落量對產(chǎn)沙量的影響。

2.3 坡面產(chǎn)沙和產(chǎn)流的關(guān)系

根據(jù)統(tǒng)計(jì)分析可以得到，無論裸坡還是有枯落物覆蓋的坡面，產(chǎn)沙量和產(chǎn)流量呈顯著相關(guān)關(guān)系。使用Matlab軟件進(jìn)行擬合分析，從圖8可以看出，在裸坡和有枯落物覆蓋的坡面上，產(chǎn)沙量和產(chǎn)流量均呈現(xiàn)二次函數(shù)關(guān)系。以產(chǎn)流量x為自變量,以產(chǎn)沙量Y為因變量，在裸坡的條件下R2=1，在有枯落物覆蓋的情況下R2=0.935 8。

圖8 產(chǎn)沙量與產(chǎn)流量的關(guān)系Fig.8 Relationship between sediment yield and runoff yield

在裸坡的情況下，擬合出的曲線方程為：

Y=78.45-2.29x+0.069 83x2。

(5)

在有枯落物覆蓋的情況下，擬合出的曲線方程為：

Y=36.34-2.531x+0.062 3x2。

(6)

從上列方程可以看出，隨著產(chǎn)流量的增大，產(chǎn)沙量逐漸增多。且隨著產(chǎn)流量的增大，產(chǎn)沙量增加的速度也隨之增加。

3 結(jié)論與討論

1)在一定的坡度和降雨強(qiáng)度下，產(chǎn)流率在短時間內(nèi)迅速升高，隨后保持穩(wěn)定；產(chǎn)沙則在短時間內(nèi)迅速升高，隨后降低到一定的速率后保持穩(wěn)定。在一定的坡度和降雨強(qiáng)度條件下，坡面徑流量隨枯落物量的增大而減小，30～90 mm/h降雨強(qiáng)度下徑流量減少率為4.06%～34.21%，降雨強(qiáng)度越小，徑流量減少率越高。在一定的坡度和降雨強(qiáng)度條件下，坡面產(chǎn)沙量隨枯落物量的增大而減小，30～90 mm/h降雨強(qiáng)度下產(chǎn)沙減少率為25.78%～96.73%，降雨強(qiáng)度越小，產(chǎn)沙減少率越高。因此，枯落物對降低產(chǎn)沙率可以起到良好的效果，這是因?yàn)榭萋湮飳酉魅趿私涤陝幽?，防止雨滴直接打擊土壤，且大大減輕了徑流沖刷作用。隨著降雨強(qiáng)度的增大，枯落物對產(chǎn)沙率的影響逐漸減小，這主要是因?yàn)榻涤陱?qiáng)度是影響侵蝕過程最主要的因素，當(dāng)降雨強(qiáng)度增大時，徑流增多、侵蝕能力增強(qiáng)，相對而言枯落物的影響減小。

2)枯落物的減沙效應(yīng)高于枯落物的減少徑流作用。產(chǎn)沙量與產(chǎn)流量呈二次函數(shù)關(guān)系；在枯落物量較大的情況下，降雨強(qiáng)度較強(qiáng)時依然可以起到良好的減沙作用。在裸坡條件與有覆蓋的條件下，產(chǎn)流量均與產(chǎn)沙量呈二次函數(shù)關(guān)系，這對于水土流失量具有很好的預(yù)測作用。在何玉廣等[13]對于模擬降雨條件下側(cè)柏林地枯落物對坡面產(chǎn)流產(chǎn)沙的影響中得出，在裸坡條件下，產(chǎn)沙量和產(chǎn)流量呈指函數(shù)關(guān)系；在有枯落物覆蓋的條件下，產(chǎn)沙率和產(chǎn)流量呈線性相關(guān)關(guān)系，與本研究結(jié)果不同。這可能是因?yàn)樵诤斡駨V等實(shí)驗(yàn)中，采用的是野外側(cè)柏人工林，與本試驗(yàn)人造徑流小區(qū)不同所致。