羅丹竹，金 鑫

(1.桂林理工大學(xué)廣西環(huán)境污染控制理論與技術(shù)重點(diǎn)實(shí)驗(yàn)室，廣西 桂林 541004；2.桂林理工大學(xué)巖溶地區(qū)水污染控制與用水安全保障協(xié)同創(chuàng)新中心，廣西 桂林 541004)

土壤顆粒組成是指粒徑大小不同的顆粒在土壤中所占的比例[1]，常用黏粒、粉粒和砂粒在土壤中所占的比例表達(dá)[2]。降雨過程中，雨滴擊打土壤表面，土壤表層結(jié)構(gòu)在團(tuán)聚體破碎、分散等作用力下發(fā)生破壞，土壤顆粒沿坡面向各個(gè)方向發(fā)生遷移[3]，同時(shí)也隨著下滲水流通過孔隙向土壤下層移動(dòng)，造成原位置的土壤顆粒組成發(fā)生變化[4]。土壤顆粒遷移是一種普遍發(fā)生在土壤中的物理過程[5]，在土壤發(fā)育和土壤物質(zhì)遷移過程中起著重要的作用[6]。土壤機(jī)械組成是影響水分入滲的主要因素[7]，而不同的土壤顆粒組成，其產(chǎn)流產(chǎn)沙規(guī)律及土壤的抗蝕性能存在一定的差異[8]。一般而言，土壤顆粒組成中細(xì)顆粒含量越多，其黏結(jié)能力越強(qiáng)[9]，內(nèi)聚力增大[10]，在一定程度上能夠促使土壤形成團(tuán)粒結(jié)構(gòu)，從而增大土壤的抗打擊能力。逯海葉等[11]的研究指出，土壤顆粒級(jí)配是影響抗剪強(qiáng)度最重要的內(nèi)在因素，細(xì)粒含量越多，土壤的抗剪強(qiáng)度越大。

陳洪松等[12]的研究表明，在相同條件的降雨下，不同大小的坡度形成的徑流量不同，導(dǎo)致土壤侵蝕量及侵蝕泥沙顆粒的分布特征存在差別。一些研究結(jié)果表明，徑流在土壤表面的傳輸特性決定了先搬運(yùn)表層細(xì)顆粒[13]，細(xì)顆粒比粗顆粒易被侵蝕，侵蝕泥沙中黏粒含量顯著增加。李朝霞等[14]通過人工模擬降雨試驗(yàn)觀察了3種紅壤的土壤顆粒結(jié)構(gòu)及其土壤團(tuán)聚體在降雨過程中的變化，發(fā)現(xiàn)隨著降雨時(shí)間的增長(zhǎng)，土壤中粗顆粒和細(xì)顆粒的相互之間的分離現(xiàn)象非常嚴(yán)重，降雨時(shí)間的長(zhǎng)短會(huì)對(duì)土壤顆粒組成產(chǎn)生影響。陳曉燕等[15]通過人工模擬降雨試驗(yàn)研究紫色土降雨前后土壤顆粒組成的變化情況，發(fā)現(xiàn)降雨強(qiáng)度的大小對(duì)土壤顆粒組成產(chǎn)生影響。張騰飛等[16]在野外觀察天然降雨對(duì)坡面土壤顆粒組成的影響，發(fā)現(xiàn)坡面土壤顆粒組成中粗顆粒的含量增高。黃麗等[17]對(duì)紫色土坡地土壤顆粒的遷移進(jìn)行了研究，發(fā)現(xiàn)降雨后的表層土壤顆粒組成與降雨前相比，細(xì)顆粒含量減少，粗顆粒含量增加。李裕元[18]對(duì)土壤顆粒結(jié)構(gòu)系數(shù)與粒徑及坡長(zhǎng)的位置變化關(guān)系進(jìn)行了分析，認(rèn)為黏粒和粉粒在坡面的移動(dòng)最明顯，且粉粒隨徑流沿坡面向下坡移動(dòng)并在坡下位置發(fā)生沉積，黏粒主要是隨徑流直接從坡面流失，其流失量與坡度和坡長(zhǎng)有關(guān)。陳明華[19]對(duì)紅壤區(qū)土壤進(jìn)行可蝕性分析發(fā)現(xiàn)，最易流失的是粗粉粒和細(xì)粉粒。還有一些研究發(fā)現(xiàn)，降雨過程中，侵蝕量與泥沙顆粒組成密切相關(guān)，主要表現(xiàn)為侵蝕量較大時(shí)，泥沙顆粒分布與降雨前表層土壤的粒徑分布更接近[20]。

從目前的研究來看，在眾多降雨徑流對(duì)土壤顆粒組成影響的研究中，系統(tǒng)開展雨強(qiáng)和坡度這兩個(gè)主要影響因素對(duì)坡面土壤顆粒組成影響的研究尚不多見。因此，本研究利用人工模擬降雨裝置和變坡土槽開展實(shí)驗(yàn)，通過觀測(cè)不同雨強(qiáng)和坡度組合條件下，降雨前后黃土坡面土壤顆粒組成的變化情況，探索坡面土壤顆粒組成在雨強(qiáng)和坡度綜合影響下的變化規(guī)律。

1 材料與方法

1.1 實(shí)驗(yàn)材料

實(shí)驗(yàn)采用人工模擬降雨裝置以及變坡土槽開展。降雨裝置采用旋轉(zhuǎn)下噴式噴頭，以揚(yáng)程為50m的潛水泵提供動(dòng)力，降雨高度為6m，雨強(qiáng)范圍為30～200mm/h，雨滴直徑為0.1～6mm；實(shí)驗(yàn)土槽尺寸為4m×1.2m×0.8m(長(zhǎng)×寬×深)，坡度變化范圍0～30°，底部具備滲漏功能，末端設(shè)出水口用于收集坡面徑流。

實(shí)驗(yàn)土壤為陜西省榆林市子洲縣岔巴溝流域的黃土。岔巴溝流域地貌為黃土丘陵溝壑區(qū)，土壤以黃土母質(zhì)發(fā)育的黃綿土為主，土壤顆粒組成中以粉粒為主要部分，土壤結(jié)構(gòu)疏松，植被稀少，土壤易發(fā)生侵蝕。

1.2 實(shí)驗(yàn)方法

實(shí)驗(yàn)雨強(qiáng)為60、90、120mm/h，坡度為10°、20°，采用不同的雨強(qiáng)和坡度進(jìn)行組合，降雨歷時(shí)60min，共進(jìn)行6場(chǎng)降雨實(shí)驗(yàn)。

1.3 實(shí)驗(yàn)步驟

(1)每場(chǎng)降雨實(shí)驗(yàn)進(jìn)行前，在上、中、下三個(gè)坡位隨機(jī)各選擇1個(gè)取樣點(diǎn)，用直徑為3cm的取土鉆分別在3個(gè)取樣點(diǎn)采集表層土壤樣品，將土樣自然風(fēng)干，標(biāo)號(hào)并裝袋。

(2)調(diào)整實(shí)驗(yàn)土槽的坡度到設(shè)計(jì)坡度后進(jìn)行降雨實(shí)驗(yàn)。

(3)降雨結(jié)束后在上坡、中坡、下坡三個(gè)坡位隨機(jī)各選擇1個(gè)取樣點(diǎn)，采集表層土壤樣品，待土樣自然風(fēng)干，標(biāo)號(hào)并裝袋。

(4)稱量2g經(jīng)自然風(fēng)干的土樣，加入4mL濃度為0.5%的六偏磷酸鈉進(jìn)行分散，靜置過夜，用激光粒度儀進(jìn)行土壤顆粒組成測(cè)定。每個(gè)土壤樣品進(jìn)行3次測(cè)定，取算術(shù)平均值以消除誤差。

2 實(shí)驗(yàn)結(jié)果

根據(jù)激光粒度儀的測(cè)定結(jié)果，得到每份樣品的黏粒、粉粒、砂粒的體積百分比，將同場(chǎng)次降雨前后上、中、下三個(gè)坡位的測(cè)定結(jié)果進(jìn)行算術(shù)平均，得到相應(yīng)場(chǎng)次降雨前后的坡面土壤顆粒組成的體積百分比。表1為60、90、120mm/h雨強(qiáng)下，降雨前后10°、20°坡面土壤顆粒組成的體積百分比變化結(jié)果。

表1 降雨前后坡面表層土壤機(jī)械組成變化

從實(shí)驗(yàn)結(jié)果可知，降雨后的坡面表層土壤顆粒組成變化的普遍特點(diǎn)是：黏粒含量均為下降，但變換幅度小于0.5%；粉粒含量下降相對(duì)較為明顯，下降幅度在1.9%～6.39%之間；砂粒含量相對(duì)有較為明顯的上升，上升幅度在1.44%～7.48%之間；坡面表層土壤在雨后發(fā)生粗化。

2.1 坡度對(duì)坡面表層土壤顆粒組成變化的影響

同雨強(qiáng)不同坡度條件下降雨前后的土壤顆粒組成變化如圖1所示。

圖1 同雨強(qiáng)不同坡度條件下降雨前后的土壤顆粒組成變化

從圖1可以看出，黏粒在土壤顆粒組成中所占的比例在所有場(chǎng)次降雨前后的變化均較小，表明降雨徑流所搬運(yùn)攜帶的黏粒比例與原狀土非常接近，侵蝕過程對(duì)黏粒的分選作用很小，坡度對(duì)黏粒含量變化的影響不顯著；粉粒和砂粒含量的變化表明，降雨后流失的粉粒多于砂粒，降雨徑流對(duì)粉粒和砂粒有較為明顯的分選作用，直接導(dǎo)致坡面表層土壤粉粒含量下降，砂粒含量上升，坡面表層土壤粗化，但坡度對(duì)粉粒和砂粒含量變化的影響不一致。

在60mm/h和90mm/h雨強(qiáng)時(shí)，隨著坡度增加，坡面土壤中粉粒占比下降的幅度和砂粒占比上升的幅度增大；120mm/h雨強(qiáng)下，隨著坡度的增加，粉粒和砂粒的變化幅度縮小，與60mm/h和90mm/h雨強(qiáng)條件下的趨勢(shì)相反，表明在坡面徑流量足夠大的情況下，隨著坡度增加使流速增大，徑流的侵蝕搬運(yùn)能力增大，對(duì)土壤顆粒的分選作用減弱，大雨強(qiáng)大坡度情況下出現(xiàn)泥沙俱下現(xiàn)象。

2.2 雨強(qiáng)對(duì)坡面表層土壤顆粒組成變化的影響

同坡度不同雨強(qiáng)條件下降雨前后的土壤顆粒組成變化如圖2所示。

圖2 同坡度不同雨強(qiáng)條件下降雨前后的土壤顆粒組成變化

總體來看，同坡度下，雨強(qiáng)對(duì)土壤顆粒組成的變化影響表現(xiàn)不一致。在10°時(shí)，隨著雨強(qiáng)增大，坡面表層土壤粉粒和砂粒含量變化幅度增大，即在10°坡度時(shí)，隨著雨強(qiáng)增大，降雨徑流的分選作用增大，有更多的砂粒留在了原地，而更多的粉粒被搬離原位，說明在10°坡面時(shí)，雨強(qiáng)增大導(dǎo)致徑流量增加，增大的侵蝕搬運(yùn)能力還不足以對(duì)砂粒的輸移產(chǎn)生明顯影響；20°坡時(shí)，隨著雨強(qiáng)增大，坡面表層土壤粉粒和砂粒含量變化幅度逐步減小，與10°坡面的趨勢(shì)相反，說明在大坡度情況下，隨著雨強(qiáng)增大，侵蝕搬運(yùn)能力進(jìn)一步增加，降雨徑流對(duì)土壤顆粒的分選作用降低。10°和20°坡面隨著雨強(qiáng)增大，對(duì)土壤顆粒分選作用相反的現(xiàn)象，說明坡度和雨強(qiáng)綜合作用下的坡面土壤顆粒組成變化較為復(fù)雜。

2.3 雨強(qiáng)和坡度對(duì)降雨前后坡面土壤顆粒組成分形維數(shù)的影響

土壤顆粒組成直接用不同粒徑體積百分比表達(dá)，難以進(jìn)行準(zhǔn)確的定量化描述，而引入分形理論可以在一定程度上對(duì)土壤顆粒組成定量化[21]。

因此，為定量化描述坡面土壤降雨前后的土壤機(jī)械組成變化，本文使用土壤分形維數(shù)量化土壤質(zhì)地，根據(jù)激光粒度分析儀所獲的粒徑體積的數(shù)據(jù)和土壤分形維數(shù)模型來確定土壤分形維數(shù)[22]。根據(jù)研究[23]，分形維數(shù)越小表示土壤質(zhì)地越粗，而土壤質(zhì)地越細(xì)分形維數(shù)越大，土壤粒徑組成的分形維數(shù)可以作為表征土壤結(jié)構(gòu)的重要指標(biāo)。

土壤分形維數(shù)模型為[23]：

(1)

式中，V(r

分別以lg[V(r

不同條件下降雨前后坡面表層土壤分形維數(shù)的變化結(jié)果見表2。

表2 不同條件下降雨前后坡面表層土壤分形維數(shù)

從表2的分形維數(shù)變化結(jié)果來看，所有場(chǎng)次降雨后坡面表層土壤的分形維數(shù)均低于降雨前，降雨前的平均分形維數(shù)為2.421，降雨后的平均值為2.329，也即降雨徑流使得坡面土壤粗化。分形維數(shù)變化最大的是60mm/h雨強(qiáng)20°坡面，分形維數(shù)下降了0.150，其余場(chǎng)次的分形維數(shù)下降幅度在0.043～0.113之間。不同雨強(qiáng)不同坡度條件下，降雨后坡面土壤粗化程度從大到小的排列順序?yàn)椋孩?0mm/h雨強(qiáng)20°坡面，②90mm/h雨強(qiáng)20°坡面，③120mm/h雨強(qiáng)10°坡面，④90mm/h雨強(qiáng)10°坡面，⑤120mm/h雨強(qiáng)20°坡面，⑥60mm/h雨強(qiáng)10°坡面。

2.3.1坡度對(duì)降雨前后坡面土壤分形維數(shù)的影響

同雨強(qiáng)不同坡度條件下降雨前后的土壤分形維數(shù)變化如圖3所示。

圖3 同雨強(qiáng)不同坡度條件下降雨前后的土壤分形維數(shù)變化

同雨強(qiáng)下，隨著坡度的增加，分形維數(shù)下降幅度的變化趨勢(shì)不一致。在60mm/h雨強(qiáng)下，分形維數(shù)下降幅度從10°坡面的0.043增加到20°坡面的0.150，土壤粗化程度增加；90mm/h雨強(qiáng)條件下，坡度對(duì)土壤分形維數(shù)變化的影響與60mm/h雨強(qiáng)類似，同樣表現(xiàn)為坡度增加分形維數(shù)下降幅度增大，但下降的增幅小于60mm/h雨強(qiáng)，從10°坡面的0.077增加到20°坡面的0.113；120mm/h雨強(qiáng)下，隨著坡度的增加，坡面土壤分形維數(shù)的下降幅度的變化趨勢(shì)與60mm/h和90mm/h雨強(qiáng)相反，呈現(xiàn)出隨著坡度的增加，分形維數(shù)下降幅度減小，下降幅度從10°坡面的0.110減小到20°坡面的0.053，土壤粗化程度降低。

2.3.2降雨強(qiáng)度對(duì)降雨前后坡面土壤分形維數(shù)的影響

同坡度不同雨強(qiáng)條件下降雨前后的土壤分形維數(shù)變化如圖4所示。

從圖4中可以看出，隨著雨強(qiáng)的增大，10°和20°坡度條件下的坡面土壤分形維數(shù)下降幅度的變化趨勢(shì)相反，說明雨強(qiáng)對(duì)土壤分形維數(shù)變化的影響同時(shí)受到坡度變化的影響。10°坡面分形維數(shù)下降幅度從60mm/h雨強(qiáng)的0.043增加到120mm/h雨強(qiáng)的0.110，而20°坡面分形維數(shù)下降幅度從60mm/h雨強(qiáng)的0.150下降到120mm/h雨強(qiáng)的0.053。坡面土壤粗化是雨強(qiáng)和坡度綜合作用的結(jié)果，大坡度大雨強(qiáng)對(duì)土壤粗化作用小于大坡度小雨強(qiáng)，而小坡度大雨強(qiáng)對(duì)土壤粗化作用大于小坡度小雨強(qiáng)。

圖4 同坡度不同雨強(qiáng)條件下降雨前后的土壤分形維數(shù)變化

3 分析與討論

3.1 同雨強(qiáng)條件下，隨著坡度的增加，土壤顆粒組成變化的趨勢(shì)不一致現(xiàn)象分析

小雨強(qiáng)(60、90mm/h)條件下，隨著坡度增加，坡面土壤粗化程度(分形維數(shù)下降說明土壤變粗，砂粒含量上升也說明土壤變粗)上升，而大雨強(qiáng)(120mm/h)條件下，隨著坡度增加，坡面土壤粗化程度下降。這種土壤粗化趨勢(shì)不一致是徑流流速和流量均發(fā)生變化后，坡面徑流的水流功率變化造成的。

根據(jù)研究[25]，水流功率是計(jì)算水流輸沙強(qiáng)度的最好指標(biāo)。水流功率是決定水流對(duì)土壤顆粒沖刷、搬運(yùn)能力和分選作用大小的關(guān)鍵因素，其與流量成正比，也與水流速度的平方成正比。坡面薄層水流流速與流量的1/4次冪和坡度的3/8次冪成正比[26]。在小雨強(qiáng)條件下，隨著坡度增加，水流流速也相應(yīng)增加的，但流量小于大雨強(qiáng)條件下的流量。因此，隨著坡度的增加，小雨強(qiáng)下水流功率增加的幅度明顯低于大雨強(qiáng)下水流功率的增加幅度。對(duì)于坡面土壤顆粒來說，顆粒越粗越容易被分離，越不容易被搬運(yùn)；顆粒越細(xì)越不容易被分離，但越容易被搬運(yùn)[27]。60mm/h和90mm/h雨強(qiáng)條件下，隨著坡度增加，水流功率有所增大，分離的粉粒含量增加，但水流功率還沒有增大到能夠搬運(yùn)更多砂粒的程度，搬運(yùn)的砂粒比例并沒有明顯增加，而坡度增加后提供的可供搬運(yùn)粉粒增加，造成了在小雨強(qiáng)條件下，隨著坡度增加，對(duì)于粉粒和砂粒的分選作用增加的現(xiàn)象。在大雨強(qiáng)條件下，隨著坡度增加，增加的水流功率足夠大，對(duì)于砂粒的搬運(yùn)能力明顯提高，分選作用減小。黏粒的分離主要是伴隨著砂粒和粉粒的分離而分離，水流并不存在單獨(dú)增加對(duì)黏粒分離作用的情況，使得黏粒含量在降雨后并未發(fā)生明顯變化。

3.2 同坡度下，隨著雨強(qiáng)的增加，土壤顆粒組成變化的趨勢(shì)不一致現(xiàn)象分析

10°坡面和20°坡面隨著雨強(qiáng)的增加，土壤粗化的變化趨勢(shì)相反，此現(xiàn)象仍可用水流功率的變化進(jìn)行解釋。在雨強(qiáng)增加幅度相同的情況下，坡度成為水流功率的決定因素，隨著雨強(qiáng)的增加，10°坡面的水流功率總體小于20°坡面，即10°坡面隨著雨強(qiáng)的增加流量增加，導(dǎo)致水流功率增大，增加了對(duì)粉粒的分離，但水流功率對(duì)砂粒的搬運(yùn)能力還未發(fā)生質(zhì)變，使得在坡度較小的情況下，隨著雨強(qiáng)的增加土壤粗化程度上升；20°坡面隨著雨強(qiáng)的增加，水流功率能夠粉粒更多的粉粒，同時(shí)水流功率也增大到能夠搬運(yùn)更多比例的砂粒程度，造成在20°坡度時(shí)，隨著雨強(qiáng)的增加土壤粗化程度下降的現(xiàn)象。

3.3 土壤分形維數(shù)變化結(jié)果分析

根據(jù)按照分形維數(shù)變化進(jìn)行的土壤粗化程度排序結(jié)果，中小雨強(qiáng)大坡度的坡面土壤粗化程度大于中大雨強(qiáng)小坡度的坡面，大雨強(qiáng)大坡度和小雨強(qiáng)小坡度的坡面粗化程度最低。小雨強(qiáng)小坡度條件下，水流功率以及挾沙力均最小，對(duì)細(xì)顆粒的剝離能力最弱，水流并沒有更多的細(xì)顆粒土壤可搬運(yùn)，也沒有搬運(yùn)更多粗顆粒的能力，使得小雨強(qiáng)小坡度條件下，坡面土壤的粗化程度最低；大雨強(qiáng)大坡度條件下，水流功率和總的挾沙能力都是最大的，可以剝離的細(xì)顆粒最多，而能夠搬運(yùn)的粗顆粒也最多，分選作用減弱，坡面土壤粗化程度相對(duì)較低；其他坡度和雨強(qiáng)組合條件下，水流功率增加則分離粉粒和搬運(yùn)砂粒的能力均有所增加，粉粒被分離后大部分可以被水流搬運(yùn)走，而水流功率與砂粒的啟動(dòng)輸移條件對(duì)比關(guān)系決定了是否能夠搬運(yùn)更多砂粒，這種在坡度和雨強(qiáng)綜合作用下的水流功率變化造成了不同條件下坡面土壤粗化程度的不同。

3.4 土壤顆粒組成變化的臨界條件

實(shí)驗(yàn)的結(jié)果顯示，雨強(qiáng)和坡度條件改變，土壤顆粒組成變化的趨勢(shì)不一致。由于趨勢(shì)出現(xiàn)反轉(zhuǎn)，則必然存在反轉(zhuǎn)的臨界點(diǎn)。從現(xiàn)象來看，可能在10°和20°之間存在臨界坡度，或者在90mm/h和120mm/h之間存在臨界雨強(qiáng)。臨界現(xiàn)象實(shí)質(zhì)上是坡度和雨強(qiáng)的改變引起了流量、水深、流速、切應(yīng)力、挾沙力、水流功率、水流流態(tài)等一系列的復(fù)雜變化，也影響了坡面的產(chǎn)流和入滲。因此，有必要進(jìn)一步開展有關(guān)臨界條件及其作用機(jī)理的研究。

4 結(jié)語

(1)所有場(chǎng)次降雨實(shí)驗(yàn)雨后的坡面土壤均為砂粒含量上升，粉粒含量下降，而黏粒含量雖有下降，但變化很小。

(2)分形維數(shù)變化排序的結(jié)果反映出坡度變化對(duì)土壤粗化程度的影響大于雨強(qiáng)變化產(chǎn)生的影響。

(3)隨著雨強(qiáng)的增加，不同坡度坡面的土壤顆粒組成變化的趨勢(shì)不一致，而隨著坡度的增加，不同雨強(qiáng)條件下的坡面土壤顆粒組成變化的趨勢(shì)也不一致，提示可能存在影響坡面土壤顆粒組成變化的臨界坡度或者臨界雨強(qiáng)。