耿 韌，張光輝，，郁耀闖，王 浩

（1.北京師范大學(xué) 地理學(xué)與遙感科學(xué)學(xué)院，北京100875；2.中國(guó)科學(xué)院 水利部水土保持研究所，陜西 楊凌712100）

土壤侵蝕是黃土高原面臨的主要環(huán)境問(wèn)題之一，它包括土壤分離、泥沙輸移和泥沙沉積三大過(guò)程［1］。在降雨擊濺和徑流沖刷作用下土壤顆粒脫離母質(zhì)的過(guò)程為土壤分離過(guò)程［2］，它為泥沙輸移和泥沙沉積提供了物質(zhì)準(zhǔn)備，是土壤侵蝕過(guò)程的起始階段。土壤分離過(guò)程是降雨特性、地形條件、坡面流水動(dòng)力條件、土壤性質(zhì)、植被覆蓋的函數(shù)［3］。土壤分離速率隨著坡面流流量、水深、流速的增大而增大，在特定降雨和下墊面條件下，坡面流流量、水深、流速等水動(dòng)力學(xué)特征又與坡面水文過(guò)程密切相關(guān)，無(wú)論是結(jié)皮的發(fā)育、地表積水、產(chǎn)流歷時(shí)、入滲及產(chǎn)流過(guò)程均受控于土壤物理性質(zhì)（土壤含水量、容重、粘結(jié)力、水穩(wěn)團(tuán)聚體）及其時(shí)空變化特征［4－6］。

土地利用方式顯著影響土壤分離過(guò)程。黃土高原5種典型土地利用方式對(duì)土壤分離速率的影響，發(fā)現(xiàn)農(nóng)地最容易分離，其分離速率分別是草地、灌木、荒地和林地分離速率的2.05，2.76，3.23，13.31倍，土壤細(xì)溝可蝕性和臨界剪切力均受到土地利用方式的顯著影響［7］。黃土高原典型土地利用條件下土壤分離速率的季節(jié)變化特征，結(jié)果發(fā)現(xiàn)無(wú)論是農(nóng)地還是非農(nóng)地，土壤分離速率均呈現(xiàn)顯著的季節(jié)波動(dòng)，主要驅(qū)動(dòng)力為農(nóng)事活動(dòng)、土壤硬化和植被根系生長(zhǎng)［8］。耕作、播種、鋤草、收獲等農(nóng)事活動(dòng)，擾動(dòng)地表從而增大土壤分離速率，影響的大小取決于農(nóng)事活動(dòng)對(duì)地表的擾動(dòng)程度；在降雨打擊及重力的作用下，耕作后疏松的土壤逐漸沉降，土壤容重和粘結(jié)力逐漸增大，導(dǎo)致土壤分離速率逐漸減?。恢脖桓稻哂泻軓?qiáng)的抗蝕能力，隨著根系密度的增大，土壤分離速率呈指數(shù)函數(shù)下降［9］。但目前關(guān)于土壤物理性質(zhì)季節(jié)變化的定量化研究較少，特別是對(duì)于土壤物理性質(zhì)變化顯著的農(nóng)地更是如此［10］，尚無(wú)法從機(jī)理上理解、模擬、預(yù)測(cè)土壤分離過(guò)程的季節(jié)變化，因此亟待展開(kāi)不同氣候、作物、田間管理?xiàng)l件下土壤物理性質(zhì)的季節(jié)變化特征的研究，為解釋、分析土壤分離過(guò)程的動(dòng)態(tài)變化，檢驗(yàn)或建立土壤分離季節(jié)變化預(yù)報(bào)方程或算法提供數(shù)據(jù)或理論支持。

1 試驗(yàn)區(qū)概況

試驗(yàn)在中國(guó)科學(xué)院水利部水土保持研究所安塞水土保持綜合試驗(yàn)站進(jìn)行。該站位于黃土高原丘陵溝壑區(qū)第二副區(qū)（109°19′24.4″E，36°51′29.9″N），海拔1 068～1 309m。屬暖溫帶半干旱氣候，年均氣溫8.8℃，年均降雨量500mm，60%的降水集中于6—9月，年際變化幅度較大。研究區(qū)土壤類型為黃綿土，黏粒含量為5.5%，粉粒含量為67.0%，沙粒含量為26.3%，有機(jī)質(zhì)含量為0.9% 。受強(qiáng)烈的水土流失影響，地形千溝萬(wàn)壑，支離破碎。研究區(qū)主要農(nóng)作物為玉米、谷子、大豆和土豆。

2 研究方法

在安塞水土保持綜合試驗(yàn)站選取玉米、谷子、大豆和土豆4種典型農(nóng)地，分別進(jìn)行土壤含水量、土壤容重、土壤水穩(wěn)性團(tuán)聚體和土壤粘結(jié)力的測(cè)定。4種農(nóng)地按照當(dāng)?shù)卮筇镒魑镞M(jìn)行施肥和管理，表1給出了4種農(nóng)地主要的農(nóng)事活動(dòng)及日期。試驗(yàn)從2012年4月8日開(kāi)始，持續(xù)到9月29日，時(shí)間間隔大約20d。

表1 4種農(nóng)地的農(nóng)事活動(dòng)

土壤含水量：用烘干稱重法測(cè)定，試驗(yàn)時(shí)用小鋁盒從上到下均勻采集0—5cm表層土壤，用烘干法測(cè)定土壤含水量，每種農(nóng)地每次重復(fù)測(cè)定5次，取其平均值作為本次測(cè)定的土壤含水量。

土壤容重：用體積為100cm3的不銹鋼環(huán)刀（高5.1cm，直徑5cm）取0—5cm表層土壤，使土充滿環(huán)刀，稱濕重后放置于烘箱中在105℃下烘干24h，每次每種農(nóng)地重復(fù)測(cè)定3次，取其平均值。

水穩(wěn)性團(tuán)聚體的測(cè)量步驟如下：

（1）采樣：用鐵鏟取耕層土壤。土壤濕度要求最好不沾鐵鏟。用鐵質(zhì)飯盒采集3個(gè)有代表性的原狀土樣，運(yùn)回實(shí)驗(yàn)室，輕輕的掰成10～12mm的土塊，攤平，使其自然風(fēng)干。

（2）干篩：取混合均勻的風(fēng)干土樣，過(guò)孔徑分別為10，7，5，3，2，1，0.5，0.25mm 篩子。分別稱重各孔徑對(duì)應(yīng)的團(tuán)聚體，計(jì)算各孔徑所占的質(zhì)量百分比。并依據(jù)各級(jí)孔徑所占百分比把各孔徑上團(tuán)聚體配成2份50g的土樣。

（3）濕篩：本文使用的濕篩分析儀為FT-3型電動(dòng)團(tuán)聚體分析儀，濕篩的孔徑為5，3，2，1，0.5，0.25 mm。先將濕篩組用鐵架夾住放入水桶，把濕篩組升到最高處，向桶內(nèi)緩慢注水，直到水面高于篩面1cm為止。將配好的50g土放入水中，震蕩30次，把濕篩組取出，用清水洗出各孔徑上的團(tuán)聚體，烘干后，稱重每個(gè)孔徑濕篩的團(tuán)聚體［11］，再把每個(gè)孔徑濕篩的團(tuán)聚體相加，然后取兩個(gè)重復(fù)的平均值。

粘結(jié)力：采用由荷蘭生產(chǎn)的微型粘結(jié)力儀（Durham Geo-enterprises，Inc.，UK），選取有代表性的平地，用噴壺噴濕地表，將粘結(jié)力儀的葉輪垂直插入土壤，順時(shí)針轉(zhuǎn)動(dòng)粘結(jié)力儀，當(dāng)扭力大于土壤最大抗剪強(qiáng)度時(shí)，粘結(jié)力儀葉輪開(kāi)始轉(zhuǎn)動(dòng)，此時(shí)指針?biāo)傅臄?shù)值即為此次測(cè)定的粘結(jié)力，每種農(nóng)地重復(fù)10次，取平均值。

本文用變異系數(shù)來(lái)描述農(nóng)地土壤物理性質(zhì)隨季節(jié)的波動(dòng)情況，它反映了數(shù)據(jù)的離散程度。用峰度和偏度描述數(shù)據(jù)的正態(tài)分布情況，峰度反映了數(shù)據(jù)在均值附近集中情況，偏度反映了數(shù)據(jù)偏離以平均值為中心的程度，它們由下式計(jì)算：比標(biāo)準(zhǔn)正態(tài)分布低，但γkurtosis＞0時(shí)是數(shù)據(jù)集中程度比標(biāo)準(zhǔn)正態(tài)分布高。當(dāng)γskewness＝0時(shí)是對(duì)稱分布，當(dāng)γskewness＜0時(shí)是數(shù)據(jù)的均值在峰值的左邊，當(dāng)γskewness＞0時(shí)是數(shù)據(jù)的均值在峰值的右邊。

數(shù)據(jù)處理使用IBM SPSS Statistics 19.0，并對(duì)數(shù)據(jù)進(jìn)行配對(duì)樣本T檢驗(yàn)分析不同農(nóng)地土壤物理性質(zhì)變化趨勢(shì)的差異性，Pearson相關(guān)分析檢驗(yàn)容重和粘結(jié)力的相關(guān)性。

3 結(jié)果與分析

式中：Cv——變異系數(shù)；γkurtosis——峰度；γskewness——偏度；SD——標(biāo)準(zhǔn)差；ˉX——樣本平均值；Xi——第i個(gè)樣本值；n——樣本個(gè)數(shù)。

根據(jù)Nielsen和Bouma提出的分類體系［12］，當(dāng)Cv≤0.1時(shí)是弱變異性，當(dāng)0.1＜Cv＜1時(shí)是中等變異性，當(dāng)Cv≥1是強(qiáng)變異性。當(dāng)γkurtosis＝0是標(biāo)準(zhǔn)正態(tài)分布的峰度系數(shù)，當(dāng)γkurtosis＜0時(shí)是數(shù)據(jù)集中程度

3.1 土壤含水量

圖1給出了4種農(nóng)地土壤含水量的季節(jié)變化，從圖中可以明顯的看出，4種農(nóng)地土壤含水量具有大致相同的變化趨勢(shì)，從4月初到9月下旬，土壤含水量出現(xiàn)了3個(gè)高峰和2個(gè)低谷，3個(gè)高峰均出現(xiàn)降雨之后，而土壤含水量的低谷則對(duì)應(yīng)著一段時(shí)間的干旱期。這表明該地區(qū)農(nóng)地表層土壤含水量的季節(jié)變化受降雨的影響強(qiáng)烈，土壤含水量變化趨勢(shì)基本上與降雨變化趨勢(shì)一致，這與蘇敏［13］和 Coronato［14］等的研究結(jié)果一致。土壤含水量的干濕交替，可能引起土壤表層出現(xiàn)較大縫隙，從而影響土壤抗蝕能力，導(dǎo)致土壤分離速率增大［15］。

表2給出了4種作物土壤物理性質(zhì)的統(tǒng)計(jì)特征值。就土壤含水量平均值而言，玉米地＞大豆地＞土豆地＞谷子地。4種農(nóng)地土壤含水量變異系數(shù)介于0.43～0.49之間，大于0.1小于1，屬于中等變異性。

表2 土壤物理性質(zhì)的統(tǒng)計(jì)特征值

從表2土壤含水量的峰度值、偏度值和單樣本的Kolmogorov-Smirnov檢驗(yàn)來(lái)看，各農(nóng)地土壤含水量呈正態(tài)分布的，可以進(jìn)行統(tǒng)計(jì)學(xué)分析。表3給出了4種農(nóng)地配對(duì)樣本T檢驗(yàn)的P值，從表中可以看出玉米和谷子之間的P值為0.04，表明在0.05的顯著性水平下，兩者之間變化趨勢(shì)存在顯著差異，而其他作物配對(duì)的P值遠(yuǎn)大于0.05，表明這些作物的土壤水分隨季節(jié)變化趨勢(shì)之間無(wú)顯著差異。各農(nóng)地表層土壤水分總體上隨季節(jié)變化趨勢(shì)大體相同，但無(wú)論在均值和變異系數(shù)還是配對(duì)樣本T檢驗(yàn)結(jié)果都表現(xiàn)出一定差異性，其中玉米和谷子季節(jié)變化趨勢(shì)有顯著差異，就土壤含水量平均值而言，玉米地是谷子地的1.12倍。這可能是因?yàn)樗鼈冎g有著不同的植被蓋度、不同的根系密度和不同的蒸散強(qiáng)度［16］所致。

圖1 土壤含水量季節(jié)變化

表3 不同農(nóng)地物理性質(zhì)的配對(duì)樣本T檢驗(yàn)

3.2 土壤容重

圖2給出了不同農(nóng)地土壤容重的季節(jié)變化，從圖中可以看出4種農(nóng)地土壤容重隨時(shí)間的變化趨勢(shì)大體相同。4月8日到4月10日，受播種的影響各農(nóng)地的容重迅速下降，達(dá)到最小值。隨后受降雨擊濺及土體自重的綜合影響，各農(nóng)地土壤均表現(xiàn)出一定程度的硬化過(guò)程，各農(nóng)地的容重均呈上升趨勢(shì)（5月23日之前），但上升的程度各不相同，玉米、谷子、大豆和土豆分別上升了6.0%，14.2%，3.6%，1.8%。到6月1日左右，受鋤草的影響各農(nóng)地土壤容重都呈下降趨勢(shì)，之后各農(nóng)地土壤容重出現(xiàn)再次上升趨勢(shì)，上升的幅度也隨作物類型的不同而有所差異。在7月22日之后，各農(nóng)地之間的變化趨勢(shì)較為混亂，但整體呈現(xiàn)上升趨勢(shì)，這一方面是由硬化過(guò)程引起，另一方面與土壤結(jié)皮的發(fā)育密切相關(guān)。結(jié)皮的發(fā)育直接受控于降雨對(duì)地表的打擊作用，而各種作物的蓋度存在一定差異，從而導(dǎo)致各農(nóng)地結(jié)皮發(fā)育的差異，從而造成土壤容重較為復(fù)雜的變化趨勢(shì)。9月22日之后，各農(nóng)地土壤容重迅速減小，其原因在于作物收獲對(duì)地表引起的擾動(dòng)，其擾動(dòng)程度取決于收獲方式，土豆的收獲過(guò)程對(duì)地表擾動(dòng)最大，其容重下降幅度也最大，為15.2%?？傮w而言，生長(zhǎng)季4種農(nóng)地的土壤容重整體呈增加趨勢(shì)，這與Alletto等人的研究結(jié)果一致［17］。土壤容重的增大，提高了土壤抗蝕能力，從而導(dǎo)致土壤分離速率下降［18］。

從表2土壤容重的平均值來(lái)看，4種農(nóng)地中谷子地的容重最大，土豆地的容重最小，玉米地和大豆地居中，這可能與作物種植的密度有關(guān)。從變異系數(shù)來(lái)看，玉米地、谷子地、大豆地和土豆地變異系數(shù)依次為0.07，0.05，0.05，0.07，都小于0.1，為弱變異性。表3給出了4種農(nóng)地的容重配對(duì)樣本T檢驗(yàn)結(jié)果。除了谷子和大豆外，其他配對(duì)組合的顯著性水平P值均大于0.05，表明這些農(nóng)作物在0.05顯著性水平下沒(méi)有顯著差異。谷子和大豆間的P值為0.03，小于0.05，表明谷子和大豆地的土壤容重隨時(shí)間變化的趨勢(shì)之間存在顯著差異。這可能與不同農(nóng)地的施肥水平、播種密度、農(nóng)事活動(dòng)的強(qiáng)度等因素有關(guān)［19］。

圖2 土壤容重季節(jié)變化

3.3 土壤水穩(wěn)性團(tuán)聚體

從圖3中可以看出，播種初期各農(nóng)地的土壤水穩(wěn)性團(tuán)聚體幾乎相同，隨時(shí)間的變化都是單調(diào)遞增趨勢(shì)，這一結(jié)果與 Mulla等人的研究成果一致［20］。團(tuán)聚體通過(guò)其物理、化學(xué)及生物的膠結(jié)作用提高土壤抗蝕能力，農(nóng)地團(tuán)聚體隨作物生長(zhǎng)呈上升趨勢(shì)表明，隨著作物的生長(zhǎng)，農(nóng)地的抗蝕能力呈逐漸增強(qiáng)的趨勢(shì)。然而，水穩(wěn)性團(tuán)聚體隨時(shí)間的增大幅度隨著作物類型的不同而存在差異，這可從表2中的極差得到驗(yàn)證，農(nóng)作物的水穩(wěn)性團(tuán)聚體最小值差別不大，極差值（最大值減去最小值）最小的土豆地為1.32，而極差值最大的玉米地為16.51，玉米地的極差值是土豆地的極差值得12.5倍。玉米地、谷子地和大豆地的變異系數(shù)分別為0.27，0.19，0.17，介于0.1和1之間，為中等變異性。而土豆地的變異系數(shù)只有0.03，為弱變異性。4種農(nóng)地水穩(wěn)性團(tuán)聚體的平均值從大到小依次為玉米、大豆、谷子和土豆，其值分別為21.9，19.8，19.8，15.4。

圖3 水穩(wěn)性團(tuán)聚體季節(jié)變化

由表3可知，玉米地和谷子地、玉米地和大豆地、谷子地和大豆地的顯著性水平P值都大于0.05，表明它們的水穩(wěn)性團(tuán)聚體季節(jié)變化趨勢(shì)無(wú)顯著性差異。而土豆和玉米、土豆和谷子、土豆和大豆的顯著性水平P值小于0.05，表明它們的水穩(wěn)性團(tuán)聚體隨季節(jié)的變化趨勢(shì)存在顯著差異。總體而言，玉米地、谷子地和大豆地水穩(wěn)性團(tuán)聚體季節(jié)變化差異不大，而土豆地和玉米地差異明顯。這可能是因?yàn)橥炼馆^為稀疏，在單位面積上積累的有機(jī)物較少，因此團(tuán)聚體含量少［21］。此外原本水穩(wěn)性團(tuán)聚體含量低的土豆地，由于植被蓋度低，裸露的地表受到更多的雨滴打擊作用，團(tuán)聚體容易遭到破壞［22］。

3.4 土壤粘結(jié)力

圖4給出了各農(nóng)地土壤粘結(jié)力隨時(shí)間變化趨勢(shì)。由圖中可以看出4種農(nóng)地的粘結(jié)力隨時(shí)間的變化具有相同的趨勢(shì)，土壤粘結(jié)力存在兩個(gè)高峰和三個(gè)低谷。三個(gè)低谷分別對(duì)應(yīng)著4月10日的農(nóng)地翻耕、6月4日的鋤草和9月28日的作物收獲，粘結(jié)力的減少顯然與農(nóng)事活動(dòng)的擾動(dòng)密切相關(guān)［23］。從表2可以看出無(wú)論是極差和平均值，還是變異系數(shù)，4種農(nóng)地粘結(jié)力之間差異不大。4種農(nóng)地的變異系數(shù)介于0.44和0.54之間，為中度變異性。

表3給出了4種農(nóng)地粘結(jié)力配對(duì)樣本T檢驗(yàn)結(jié)果，從表中可以看出P值介于0.16和0.53之間，都遠(yuǎn)大于0.05，表明各農(nóng)地粘結(jié)隨時(shí)間變化在0.05的顯著性水平下沒(méi)有顯著性差異，這在一定程度上說(shuō)明了各農(nóng)地抗蝕能力的相似性。比較圖4和圖2可以看出，四種農(nóng)地粘結(jié)力和容重的變化趨勢(shì)具有一定的相似性，土壤容重的減少，則對(duì)應(yīng)著粘結(jié)力的減少，反之亦然。對(duì)4種農(nóng)地的粘結(jié)力與土壤容重進(jìn)行Pearson相關(guān)分析，相關(guān)系數(shù)為0.68，顯著性水平P值小于0.01，兩者相關(guān)關(guān)系顯著。這可理解為土壤容重是土壤緊實(shí)度的指標(biāo)之一，土壤緊實(shí)度低則粘結(jié)力小，土壤緊實(shí)度高則粘結(jié)力大［23］。

圖4 粘結(jié)力季節(jié)分布

4 結(jié)論

本文研究了黃土高原4種典型農(nóng)地（玉米、谷子、大豆和土豆）的表層土壤物理性質(zhì)——土壤含水量、土壤容重、土壤水穩(wěn)性團(tuán)聚體和土壤粘結(jié)力的季節(jié)變化規(guī)律。得出以下結(jié)論：（1）土壤含水量季節(jié)變化受降水影響顯著，生長(zhǎng)季存在多個(gè)峰值和低谷，呈中等變異。（2）受土壤硬化的影響，土壤容重總體呈增大趨勢(shì)，但農(nóng)事活動(dòng)可引起土壤容重的顯著減小，影響程度取決于農(nóng)事活動(dòng)對(duì)地表的擾動(dòng)程度，4種農(nóng)地的土壤容重均屬弱變異。（3）土壤團(tuán)聚體隨著作物的生長(zhǎng)呈上升趨勢(shì)，土豆地團(tuán)聚體顯著低于其他3種農(nóng)地，除土豆地屬弱變異外，其他3種農(nóng)地的水穩(wěn)性團(tuán)聚體均屬于中等變異。（4）作物生長(zhǎng)季，土壤粘結(jié)力呈上升趨勢(shì)，但受農(nóng)事活動(dòng)的影響顯著，農(nóng)事活動(dòng)的擾動(dòng)會(huì)導(dǎo)致土壤粘結(jié)力顯著減小?？傮w而言，黃土高原典型農(nóng)地的土壤物理性質(zhì)具有較為相似的季節(jié)變化特征，主要影響因素包括農(nóng)事活動(dòng)、土壤硬化及根系生長(zhǎng)引起的有機(jī)質(zhì)積累，在生長(zhǎng)季農(nóng)地抗蝕能力呈增強(qiáng)趨勢(shì)。

［1］ 張光輝.坡面水蝕過(guò)程水動(dòng)力學(xué)研究進(jìn)展［J］.水科學(xué)進(jìn)展，2001，12（3）：39-402.

［2］ Owoputi L O，Stolte W J.Soil detachment in the physically based soil erosion process：a review［J］.Transactions of the ASAE，1995，38（4）：1099-1110.

［3］ Ciampalini R，Torri D.Detachment of soil particles by shallow flow：Sampling methodology and observations J .Catena199832137-53.

［4］ Hanks R J.Soil crusting and seedling emergence［J］.Transactions 7th int.Congr.Soil Sci.，1960，1：340-346.

［5］ De Ploey J，Poesen J.Aggregate stability，runoff generation and interrill erosion［C］∥Richards K S，Arnett R R，Ellis S.Geomorphology and Soils.London：University of Hull，1985：99-120.

［6］ Hu W，Shao M A，Si B C.Seasonal changes in surface bulk density and saturated hydraulic conductivity of natural landscapes［J］.European Journal of Soil Science，2012，63（6）：820-830.

［7］ Zhang G H，Liu G B，Tang K M，et al.Flow detachment of soils under different land uses in the Loess Plateau of China［J］.Transactions of the ASABE，2008，51（3）：883-890.

［8］ Zhang G，Tang M，Zhang X C.Temporal variation in soil detachment under different land uses in the Loess Plateau of China［J］.Earth Surface Processes and Landforms，2009，34（9）：1302-1309.

［9］ Gyssels G，Poesen J，Bochet E，et al.Impact of plant roots on the resistance of soils to erosion by water：a review［J］.Progress in Physical Geography，2005，29（2）：189-217.

［10］ 張光輝，劉國(guó)彬.黃土丘陵區(qū)小流域土壤表面特性變化規(guī)律研究［J］.地理科學(xué)，2001，21（2）：118-122.

［11］ Yoder R E.The significance of soil structure in relation to the tilth problem［J］.Soil Science Society of America Journal，1938，2（C）：21-33.

［12］ Nielsen D R，Bouma J.Soil spatial variability［M］.The Netherlands：Wageningen，1985：243.

［13］ 蘇敏，盧宗凡，李夠霞.陜北丘陵溝壑區(qū)主要農(nóng)作物水分利用與平衡［J］.水土保持研究，1996，3（2）：36-45.

［14］ Coronato F R，Bertiller M B.Precipitation and land-scape related effects on soil moisture in semi-arid rangelands of Patagonia［J］.Journal of Arid Environments，1996，34（1）：1-9.

［15］ Kemper W D，Trout T J，Brown M J，et al.Furrow erosion and water and soil management［J］.Transaction of the ASAE，1985，28（5）：1564-1572.

［16］ 傅伯杰，王軍，馬克明.黃土丘陵區(qū)土地利用對(duì)土壤水分的影響［J］.中國(guó)科學(xué)基金，1999（4）：22-227.

［17］ Alletto L，Coquet Y.Temporal and spatial variability of soil bulk density and near-saturated hydraulic conductivity under two contrasted tillage management systems［J］.Geoderma，2009，152（1）：85-94.

［18］ Bennett S J，Casali J，Robinson K M，et al.Characteristics of actively eroding ephemeral gullies in an experimental channel［J］.Transactions of the ASAE，2000，43（3）：641-649.

［19］ 于法展，李保杰，尤海梅，等.徐州泉山自然保護(hù)區(qū)人工林下土壤容重與孔隙度時(shí)空變化研究［J］.水土保持研究，2007，14（6）：162-164.

［20］ Mulla D J，Huyck L M，Reganold J P.Temporal variation in aggregate stability on conventional and alternative farms［J］.Soil Science Society of America Journal，1992，56（5）：1620-1624.

［21］ Bissonnais Y.Aggregate stability and assessment of soil crustability and erodibility：Ⅰ.Theory and methodology［J］.European Journal of Soil Science，1996，47（4）：425-437.

［22］ Al-Durrah M M，Bradford J M.The mechanism of raindrop splash on soil surfaces［J］.Soil Science Society of America Journal，1982，46（5）：1086-1090.

［23］ Bullock M S，Nelson S D，Kemper W D.Soil cohesion as affected by freezing，water content，time and tillage［J］.Soil Science Society of America Journal，1988，52（3）：770-776.