魏 寧, 魏 霞, Huang Chihua

(1.西北農(nóng)林科技大學(xué) 理學(xué)院, 陜西 楊陵 712100; 2.蘭州大學(xué) 資源環(huán)境學(xué)院, 蘭州 730000; 3.Agronomy Dep., Purdue Univ., and USDA-ARS National Soil Erosion Research Lab., 275S Russell St, West Lafayette, IN 47907-2077,USA)

凍融是指土層由于溫度降到0℃以下和升至0℃以上而產(chǎn)生凍結(jié)和融化的一種物理地質(zhì)作用和現(xiàn)象[1-2]。土壤的凍結(jié)和融化實(shí)質(zhì)是土壤水的凍結(jié)和融化，即土壤中水分相態(tài)的變化過程。交替凍融指的是由于季節(jié)或晝夜熱量變化在凍土表層及以下一定深度形成的反復(fù)凍結(jié)——融化的土壤過程，它作為一種自然現(xiàn)象普遍存在于高緯度、高海拔地區(qū)以及大多中緯度地區(qū)。交替凍融對(duì)土壤物理性質(zhì)的影響，主要表現(xiàn)在水與冰轉(zhuǎn)化過程中體積脹縮對(duì)土壤團(tuán)聚結(jié)構(gòu)的破壞。交替凍融對(duì)土壤物理性質(zhì)的影響主要取決于交替凍融的速率、交替凍融溫度、土壤含水率、土壤容重和交替凍融循環(huán)的次數(shù)等[3]。研究表明，交替凍融作用下土壤理化性質(zhì)、結(jié)構(gòu)和質(zhì)地等均會(huì)發(fā)生改變，使土壤更容易成為侵蝕的物質(zhì)來源，進(jìn)而加劇土壤侵蝕程度[3-8]。既然交替凍融作用破壞土壤團(tuán)聚體結(jié)構(gòu)，降低土壤的團(tuán)聚體穩(wěn)定性[8-15]，那么必然會(huì)影響到土壤的臨界剪切力和土壤細(xì)溝可蝕性，進(jìn)而影響土壤侵蝕量。眾所周知，土壤臨界剪切力和土壤細(xì)溝可蝕性是表征土壤抗侵蝕能力的兩個(gè)重要參數(shù)，反映土壤對(duì)侵蝕外營力剝蝕和搬運(yùn)的易損性和敏感性，是影響土壤流失量的內(nèi)在因素和定量研究土壤侵蝕的基礎(chǔ)[16-23]。然而，目前有關(guān)凍融的研究主要集中在交替凍融作用對(duì)土壤團(tuán)聚體穩(wěn)定性[8,10-12]、土壤容重[13-14]、土壤抗沖性等[15]方面，對(duì)凍融作用影響下的土壤臨界剪切力和土壤細(xì)溝可蝕性這兩個(gè)參數(shù)的研究較少涉及。鑒于此，本文利用室內(nèi)凍融循環(huán)試驗(yàn)和變流量放水沖刷試驗(yàn)，研究凍融作用影響下土壤臨界剪切力和土壤細(xì)溝可蝕性的變化特征，旨在揭示交替凍融作用影響下土壤自身抗侵蝕能力的兩個(gè)重要參數(shù)——土壤臨界剪切力和土壤細(xì)溝可蝕性的變化趨勢(shì)與規(guī)律，為凍融作用影響下土壤侵蝕的預(yù)報(bào)提供重要參考。

1 材料與方法

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

本試驗(yàn)是在美國農(nóng)業(yè)部土壤侵蝕國家重點(diǎn)實(shí)驗(yàn)室進(jìn)行。供試土壤取自美國印第安納州西拉法葉普渡大學(xué)動(dòng)物科學(xué)研究和教育中心(87°00′54″W，40°29′24″N)的表層土壤(0—10 cm)，供試土壤>0.05,0.05～0.002，<0.002 mm的粒徑質(zhì)量百分比為14%，66%，20%。將采集的土壤剔除植物根系和雜物后自然風(fēng)干，然后過5 mm篩，配置成試驗(yàn)設(shè)計(jì)的土壤質(zhì)量含水率(10%和20%)，放置在凍融循環(huán)監(jiān)測(cè)系統(tǒng)中進(jìn)行交替凍融循環(huán)，將滿足凍融循環(huán)周期設(shè)計(jì)的土壤再次風(fēng)干，作為試驗(yàn)用土備用。凍融循環(huán)試驗(yàn)的最高溫度為4℃，最低溫度為-12℃，交替凍融的溫度變化速率分別是-12℃/24 h和4℃/24 h，凍融循環(huán)的周期為1，3,6。該凍融循環(huán)格局(溫度、溫度的變化速率、凍融循環(huán)周期等)的設(shè)計(jì)與該土壤所在地區(qū)的主要交替凍融方式相符合。

試驗(yàn)土槽[22-23]尺寸為3個(gè)0.5 m×0.045 m×0.13 m底部開孔的土槽，每個(gè)土槽的兩端有一個(gè)“V”型通道。裝土前，在土槽下部鋪填一層0.08 m厚的天然砂礫，以保持試驗(yàn)土的透水狀況接近天然坡面，再在天然沙層上鋪設(shè)一層粗纖維，然后裝填試驗(yàn)用土，試驗(yàn)用土的厚度為0.05 m，試驗(yàn)設(shè)計(jì)坡度3%。具體試驗(yàn)處理見表1，除此之外，為了保證降雨前期土壤初始含水率一致，以不經(jīng)受任何凍融作用的風(fēng)干土為參照試驗(yàn)(CK)進(jìn)行對(duì)照。

表1 試驗(yàn)處理

1.2 變流量人工放水沖刷試驗(yàn)

小流量的情況下，徑流剪切力小于土壤抗侵蝕的臨界剪切力，幾乎不會(huì)發(fā)生土壤侵蝕，當(dāng)土壤表面徑流的強(qiáng)度增加時(shí)，徑流剪切力增加，當(dāng)超過土壤臨界剪切力時(shí)就會(huì)產(chǎn)生侵蝕。因此依據(jù)土壤臨界剪切力出現(xiàn)的實(shí)際，本試驗(yàn)設(shè)計(jì)從小流量向大流量依次梯度變化。沖刷試驗(yàn)采用逐漸增大流量的方法進(jìn)行，在穩(wěn)流水箱與沖刷試驗(yàn)槽之間安裝小型流量計(jì)進(jìn)行流量控制。試驗(yàn)過程中，總共有100,150,200,250,300,350,400,450,500,550,600,700 ml/min 12個(gè)梯級(jí)流量。試驗(yàn)過程中每分鐘接取一個(gè)徑流樣品，每種流量下接取3個(gè)徑流樣品。沖刷方式為連續(xù)沖刷。

試驗(yàn)結(jié)束后，將試驗(yàn)徑流樣品稱重，加入飽和明礬溶液沉淀24 h，然后倒掉上清液，將倒掉清液的徑流泥沙樣放進(jìn)烘箱，105℃烘24 h，等徹底干燥后，取出稱重，計(jì)算徑流泥沙樣品。根據(jù)土壤侵蝕率和水流剪切力的關(guān)系確定土壤臨界剪切力和土壤細(xì)溝可蝕性。

1.3 土壤臨界剪切力、細(xì)溝可蝕性的確定

土壤臨界剪切力計(jì)算公式如下[24-26]：

Dr=Kr(τ-τc)

(1)

式中:Dr為土壤分離能力[kg/(m2·s)];Kr為土壤細(xì)溝可蝕性(s/m);τ為水流剪切力(Pa);τc為土壤臨界剪切力(Pa)。

水流剪切力計(jì)算公式如下[27-28]：

τ=ρgHS

(2)

式中:τ為水流剪切力(Pa);ρ為水的密度(kg/m3);g為重力加速度(m/s2);H為水深(m);S為水槽坡度(%)。

土壤分離能力計(jì)算公式如下[29-30]：

(3)

式中:Dr為土壤分離能力[kg/(m2·s)];Ms為時(shí)間t內(nèi)被分離的土壤干重(kg);t為坡面徑流沖刷時(shí)間(s);A為采樣器面積(m2)。

1.4 數(shù)據(jù)處理方法

采用Microsoft Excel 2010和IBM SPSS Statistics 20統(tǒng)計(jì)軟件對(duì)試驗(yàn)數(shù)據(jù)進(jìn)行統(tǒng)計(jì)分析。采用單因素方差分析(one-way ANOVA)和最小顯著差異法(least-significant difference，LSD)比較相同含水率不同凍融循環(huán)周期間的差異。

2 結(jié)果與分析

眾所周知，土壤侵蝕的發(fā)生必須具備兩個(gè)條件，即外部營力——侵蝕動(dòng)力和內(nèi)在物質(zhì)——侵蝕對(duì)象。對(duì)于水力侵蝕而言，侵蝕的外營力即是降雨及其形成的徑流，侵蝕對(duì)象即為地表土壤。徑流在沿坡面梯度方向運(yùn)動(dòng)過程中，在其運(yùn)動(dòng)方向上有徑流剪切力產(chǎn)生，徑流剪切力是徑流在坡面位置和時(shí)間的函數(shù)[24-25]。只有當(dāng)土壤表面徑流的強(qiáng)度增加，徑流剪切力增加，超過臨界剪切力τc時(shí)就會(huì)產(chǎn)生侵蝕。與之相伴隨的臨界剪力τc是表征土壤抗侵蝕能力的一個(gè)參數(shù)，其值相對(duì)比較固定，與土壤的比重、含水率、結(jié)構(gòu)等因素有關(guān)，是土壤侵蝕計(jì)算式中的一個(gè)非常重要的參數(shù)。土壤細(xì)溝可蝕性Kr也是表征土壤自身抗侵蝕能力大小的一個(gè)重要因子，是評(píng)價(jià)土壤對(duì)侵蝕敏感程度的重要指標(biāo)，表示了土壤被侵蝕的難易程度，反映土壤對(duì)侵蝕外營力剝蝕和搬運(yùn)的易損性和敏感性。但如何確定τc和Kr，至今仍沒有一個(gè)廣為接受的理論或試驗(yàn)的方法。目前，通常的做法是將侵蝕率與徑流剪應(yīng)力進(jìn)行回歸分析，取回歸直線在剪應(yīng)力軸上的截距為τc值，取回歸直線的斜率為Kr值。

因此，本文也采用目前比較通用的做法來確定土壤的臨界剪切力τc和土壤細(xì)溝可蝕性Kr。土壤侵蝕率采用“V”型水槽的放水沖刷試驗(yàn)，觀測(cè)得到土壤侵蝕量，進(jìn)而轉(zhuǎn)換得之。徑流剪切力由公式(2)確定。將計(jì)算的每種試驗(yàn)處理，不同交替凍融循環(huán)周期和不同土壤質(zhì)量含水率組合下的土壤的侵蝕率與徑流剪切力點(diǎn)繪在直接坐標(biāo)系中，見圖1。可以看出，土壤侵蝕率隨水流剪切力的增大而增大，土壤侵蝕率與水流剪切力呈線性關(guān)系。其函數(shù)的一般形式如下：

y=Ax-B

(4)

式中:y為土壤侵蝕率[g/(m2·s)];x為水流剪切力(Pa);A,B為系數(shù)。

為了進(jìn)一步分析不同交替凍融循環(huán)周期和不同的土壤質(zhì)量含水率對(duì)土壤臨界剪切力的影響，將圖1各種處理下的侵蝕率與徑流剪切力的關(guān)系式(1)進(jìn)行變形，變形為與公式Dr=Kr(τ-τc)[25]形式一致的關(guān)系式，如下所示：

y=A(x-B/A)

(5)

式中:y為土壤侵蝕率[g/(m2·s)];x為水流剪切力(Pa);A為土壤細(xì)溝可蝕性(s/m);B/A為水流臨界剪切力(Pa)。也就是說對(duì)于每場(chǎng)試驗(yàn),只有當(dāng)水流剪切力大于相應(yīng)的B/A時(shí),才會(huì)有侵蝕發(fā)生,該值也為每種處理下的土壤臨界剪切力值。

2.1 相同土壤含水率條件下凍融循環(huán)對(duì)土壤臨界剪切力和土壤細(xì)溝可蝕性的影響

表2列出了根據(jù)圖1中的各個(gè)擬合關(guān)系方程計(jì)算的土壤臨界剪切力和土壤細(xì)溝可蝕性值。由表2可知，參照試驗(yàn)CK的土壤臨界剪切力為0.300 Pa。當(dāng)土壤質(zhì)量含水率為10%，凍融循環(huán)周期為1,3,6時(shí)，對(duì)應(yīng)的土壤臨界剪切力的值分別為0.380,0.259,0.267 Pa。當(dāng)土壤質(zhì)量含水率為20%，凍融循環(huán)周期為1,3,6時(shí)，對(duì)應(yīng)的土壤臨界剪切力的值分別為0.338,0.238,0.242 Pa。分析發(fā)現(xiàn)，兩種土壤質(zhì)量含水率下，土壤經(jīng)過1個(gè)周期的凍融循環(huán)后，土壤臨界剪切力與參照試驗(yàn)的相應(yīng)值相比有所增大，當(dāng)土壤經(jīng)過3個(gè)周期的凍融循環(huán)后，土壤臨界剪切力顯著降低，隨著土壤的凍融循環(huán)周期增至6時(shí)，土壤臨界剪切力又有所增加，但其值仍小于參照試驗(yàn)CK的土壤臨界剪切力。

圖1 不同試驗(yàn)處理下土壤侵蝕率與徑流剪切力的關(guān)系

由表2可知，參照試驗(yàn)CK的土壤細(xì)溝可蝕性為0.003 0 s/m，當(dāng)土壤質(zhì)量含水率為10%，凍融循環(huán)周期為1,3,6時(shí)，對(duì)應(yīng)的土壤細(xì)溝可蝕性的值分別為0.007 9,0.002 7,0.001 5 s/m，即當(dāng)土壤經(jīng)歷1個(gè)周期的凍融循環(huán)后，土壤細(xì)溝可蝕性相對(duì)參照試驗(yàn)顯著增大，當(dāng)土壤經(jīng)歷3個(gè)和6個(gè)周期的凍融循環(huán)后，土壤臨界剪切力降低，且均小于參照試驗(yàn)CK的土壤細(xì)溝可蝕性值。同理，當(dāng)土壤質(zhì)量含水率為20%，凍融循環(huán)周期為1,3,6時(shí)，對(duì)應(yīng)的土壤細(xì)溝可蝕性的值分別為0.006 8,0.002 1,0.003 3 s/m，即當(dāng)土壤經(jīng)歷1個(gè)周期的凍融循環(huán)后，土壤細(xì)溝可蝕性增大，當(dāng)土壤繼續(xù)經(jīng)歷凍融循環(huán)作用，其周期增大至3時(shí)，土壤細(xì)溝可蝕性出現(xiàn)明顯的減小趨勢(shì)，當(dāng)凍融循環(huán)周期增大至6時(shí)，土壤細(xì)溝可蝕性又有所增加，但凍融循環(huán)周期為3,6時(shí)，土壤的細(xì)溝可蝕性均小于參照試驗(yàn)。

表2 不同試驗(yàn)處理下的土壤臨界剪切力

由以上分析可知，交替凍融循環(huán)周期對(duì)土壤臨界剪切力和土壤細(xì)溝可蝕性的影響顯著，且隨著交替凍融循環(huán)周期的不同，土壤臨界剪切力和土壤細(xì)溝可蝕性的影響趨勢(shì)與程度不同。兩種土壤含水率水平下，當(dāng)交替凍融循環(huán)周期為1時(shí)，交替凍融作用使土壤臨界剪切力和土壤細(xì)溝可蝕性顯著增大，當(dāng)交替凍融循環(huán)周期為3,6時(shí)，交替凍融作用使土壤臨界剪切力和土壤細(xì)溝可蝕性減小。

2.2 相同交替凍融循環(huán)周期下土壤含水率對(duì)土壤臨界剪切力和土壤細(xì)溝可蝕性的影響

由表2可知，當(dāng)交替凍融循環(huán)周期為1時(shí)，10%和20%兩種不同土壤含水率下的土壤臨界剪切力的值分別為0.380,0.338 Pa，均大于參照試驗(yàn)的0.300 Pa，即在兩種不同含水率下，經(jīng)過1個(gè)周期的凍融循環(huán)后，土壤臨界剪切力增大，且低含水率時(shí)土壤臨界剪切力增大的幅度較高含水率時(shí)土壤臨界剪切力增加的幅度大。當(dāng)交替凍融循環(huán)周期為3時(shí)，10%和20%兩種不同土壤含水率下的土壤臨界剪切力分別為0.259,0.238 Pa，均小于參照試驗(yàn)的0.300 Pa，且低含水率下的土壤臨界剪切力較高含水率下土壤臨界剪切力的值大。當(dāng)交替凍融循環(huán)周期為6時(shí)，10%和20%兩種不同土壤含水率下的土壤臨界剪切力分別為0.267,0.242 Pa，均小于參照試驗(yàn)的0.300 Pa，且10%含水率下的土壤臨界剪切力較20%含水率下的土壤臨界剪切力大。

同理，當(dāng)交替凍融循環(huán)周期為1時(shí)，10%和20%兩種不同土壤含水率下的土壤細(xì)溝可蝕性分別為0.007 9,0.006 8 s/m，均大于參照試驗(yàn)的0.003 0 s/m，即在兩種不同含水率處理下，土壤經(jīng)過1個(gè)周期的凍融循環(huán)后，土壤細(xì)溝可蝕性增大，且低含水率時(shí)土壤細(xì)溝可蝕性增大的幅度較高含水率時(shí)土壤細(xì)溝可蝕性增大的幅度大。當(dāng)交替凍融循環(huán)周期為3時(shí)，10%和20%兩種不同土壤含水率下的土壤細(xì)溝可蝕性分別為0.002 7,0.002 1 s/m，均小于參照試驗(yàn)的0.300 Pa，且低含水率下的土壤細(xì)溝可蝕性較高含水率下的土壤細(xì)溝可蝕性大。當(dāng)交替凍融循環(huán)周期為6時(shí)，10%和20%兩種不同土壤含水率下土壤細(xì)溝可蝕性分別為0.001 5,0.003 3 s/m，均小于參照試驗(yàn)的0.300 s/m，且10%含水率下的土壤細(xì)溝可蝕性較20%含水率下的土壤細(xì)溝可蝕性大。

綜上可知，不同土壤含水率對(duì)土壤臨界剪切力和土壤細(xì)溝可蝕性的影響程度不同，土壤質(zhì)量含水率水平為10%時(shí)的土壤臨界剪切力和土壤細(xì)溝可蝕性的變化幅度均大于土壤質(zhì)量含水率水平為20%時(shí)的土壤臨界剪切力和土壤細(xì)溝可蝕性的變化幅度。

3 討 論

本研究是探索交替凍融作用對(duì)土壤臨界剪切力和細(xì)溝可蝕性的影響，結(jié)果表明，交替凍融能夠改變土壤的臨界剪切力和細(xì)溝可蝕性。當(dāng)凍融循環(huán)周期為1時(shí)，10%和20%含水量處理下的土壤臨界剪切力和細(xì)溝可蝕性相比對(duì)照組都顯著增大。隨著凍融循環(huán)周期繼續(xù)增大至3時(shí)，兩種含水量處理下的土壤臨界剪切力和細(xì)溝可蝕性顯著降低，且明顯低于對(duì)照組。當(dāng)凍融循環(huán)周期繼續(xù)增大至6時(shí)，10%含水量的土壤臨界剪切力反而顯著增大，細(xì)溝可蝕性顯著減小，但其值小于控制組；20%含水量的土壤臨界剪切力和細(xì)溝可蝕性顯著增大。并且隨著凍融循環(huán)周期的增加，土壤臨界剪切力和細(xì)溝可蝕性的變化速率減小，即隨著凍融循環(huán)周期的增加，土壤臨界剪切力和細(xì)溝可蝕性趨于穩(wěn)定。已有研究結(jié)果也表明，歷經(jīng)一次凍融循環(huán)，土壤的團(tuán)聚體穩(wěn)定性顯著增大，隨著凍融循環(huán)周期的增加，土壤團(tuán)聚體的穩(wěn)定性降低，最終趨于穩(wěn)定。

此外，土壤臨界剪切力和土壤細(xì)溝可蝕性是表征土壤抗侵蝕能力的兩個(gè)重要參數(shù)，土壤細(xì)溝可蝕性越大，土壤越易侵蝕，土壤臨界剪切力越小，土壤越易侵蝕。從理論上來講，土壤臨界剪切力和細(xì)溝可蝕性對(duì)交替凍融的響應(yīng)應(yīng)該相反。但本研究的結(jié)果表明，在10%含水量時(shí)，隨著凍融循環(huán)周期的增加，土壤臨界剪切力和細(xì)溝可蝕性的變化趨勢(shì)相反，但在20%含水量時(shí)，隨著凍融循環(huán)周期的增加，土壤臨界剪切力和細(xì)溝可蝕性的變化趨勢(shì)相同。這可能是由于初始凍融含水量相比凍融循環(huán)周期對(duì)土壤臨界剪切力和細(xì)溝可蝕性影響較大緣故，已有研究也表明，對(duì)于同一種土壤，團(tuán)聚體穩(wěn)定性可能隨著F-T循環(huán)的增加而增加或降低，具體取決于凍融初始土壤含水量[10-12]。

4 結(jié) 論

(1) 交替凍融循環(huán)周期不同，其對(duì)土壤臨界剪切力和土壤細(xì)溝可蝕性的影響程度和趨勢(shì)也不同，當(dāng)凍融循環(huán)周期為1時(shí)，交替凍融作用使土壤臨界剪切力和土壤細(xì)溝可蝕性增大，當(dāng)凍融循環(huán)周期為3,6時(shí)，交替凍融作用使土壤臨界剪切力和土壤細(xì)溝可蝕性減小。

(2) 土壤含水率不同，其對(duì)土壤臨界剪切力和土壤細(xì)溝可蝕性的影響程度不同，本研究所涉及的兩種土壤質(zhì)量含水率下，當(dāng)交替凍融循環(huán)周期相同時(shí)，低含水率(10%)下的土壤臨界剪切力與土壤細(xì)溝可蝕性均大于高含水率(20%)下的土壤臨界剪切力與土壤細(xì)溝可蝕性。

(3) 各種試驗(yàn)處理下的土壤侵蝕率與徑流剪切力之間的關(guān)系仍然呈良好的線性關(guān)系。不同處理下，交替凍融作用對(duì)土壤臨界剪切力與土壤細(xì)溝可蝕性的影響趨勢(shì)相同。

致謝：感謝同在美國農(nóng)業(yè)部土壤侵蝕國家重點(diǎn)實(shí)驗(yàn)室吳秋菊博士和楊韶洋碩士在試驗(yàn)過程中給予的幫助!