張 莉，張青峰，2，徐 健，張 翔

（1.西北農(nóng)林科技大學(xué) 資源環(huán)境學(xué)院，陜西 楊凌 712100；2.中國(guó)科學(xué)院 水利部 水土保持研究所黃土高原土壤侵蝕與旱地農(nóng)業(yè)國(guó)家重點(diǎn)實(shí)驗(yàn)室，陜西 楊凌 712100；3.中國(guó)石油長(zhǎng)慶油田分公司，西安 710018）

黃土坡耕地地表糙度是由于人為土地管理和土壤侵蝕共同作用形成高低起伏和凹凸不平的地表形態(tài)，嚴(yán)重影響地表的下滲、地表產(chǎn)流、產(chǎn)沙及地表演化過(guò)程［1－5］。與此相反，地表洼地填水和泥沙攔截、輸移等過(guò)程也會(huì)對(duì)地表糙度產(chǎn)生影響，而這一相互影響過(guò)程是隨地表糙度空間分布特征的變化而變化的［1］，因此，開(kāi)展坡面侵蝕過(guò)程中地表糙度空間分布特征的研究，對(duì)于深入了解坡面水蝕過(guò)程具有重要意義。空間異質(zhì)性是生態(tài)學(xué)過(guò)程和格局在空間分布上的差異性和復(fù)雜性［6］。近年來(lái)已逐漸成為生態(tài)學(xué)、土壤學(xué)等領(lǐng)域研究的熱點(diǎn)問(wèn)題［7－8］。而針對(duì)坡耕地地表糙度空間異質(zhì)性的研究較少［9－10］，尤其是針對(duì)不同水蝕過(guò)程中地表糙度空間異質(zhì)性的研究［10］。研究不同水蝕過(guò)程中地表糙度的空間異質(zhì)性，能夠深入了解水蝕過(guò)程中地表糙度的空間結(jié)構(gòu)特征，空間變異性，尺度及其復(fù)雜程度，對(duì)進(jìn)一步揭示水蝕過(guò)程中地表糙度與土壤侵蝕的相互耦合關(guān)系具有重要意義。

地統(tǒng)計(jì)學(xué)是以傳統(tǒng)統(tǒng)計(jì)學(xué)為基礎(chǔ)發(fā)展起來(lái)的一門(mén)新的空間分析方法，適用于區(qū)域化變量空間分布特征的描述。將地統(tǒng)計(jì)學(xué)與GIS技術(shù)相結(jié)合，一方面能有效地揭示變量的空間分布、變異及相關(guān)性特征，另一方面還能將空間格局與生態(tài)過(guò)程聯(lián)系起來(lái)，更好地解釋空間格局對(duì)生態(tài)過(guò)程與功能的影響，被認(rèn)為是研究變量的空間異質(zhì)性最為有效的方法之一［11］。在區(qū)域土壤屬性的空間變異性研究中得到了廣泛地應(yīng)用［6，12－13］。本文試圖采用地統(tǒng)計(jì)學(xué)與 GIS技術(shù)相結(jié)合的方法，分析比較黃土人工鋤耕坡面地表糙度在不同侵蝕（雨前—濺蝕—片蝕—細(xì)溝侵蝕）演化過(guò)程中的空間異質(zhì)性特征，揭示水蝕過(guò)程中地表糙度的空間異質(zhì)性，以期為進(jìn)一步明確水蝕過(guò)程中地表糙度與土壤侵蝕的相互耦合關(guān)系奠定基礎(chǔ)。

1 研究方法

1.1 實(shí)驗(yàn)設(shè)計(jì)

本試驗(yàn)在黃土高原土壤侵蝕與旱地農(nóng)業(yè)國(guó)家重點(diǎn)實(shí)驗(yàn)室土壤侵蝕人工模擬降雨大廳進(jìn)行，所用降雨設(shè)備為側(cè)噴式人工降雨器，降雨均勻度大于80%，降雨高度為16m，可以滿足所有的雨滴達(dá)到終點(diǎn)速度。試驗(yàn)用土槽規(guī)格為2.0m×1.0m×0.5m，以陜西省楊凌區(qū)坡耕地原狀表層土（0—20cm）為實(shí)驗(yàn)用土，土壤容重為1.30g／cm3，含水率約為10%。為了使布設(shè)的耕作措施更接近自然狀況，試驗(yàn)中聘請(qǐng)長(zhǎng)期從事相同耕作措施農(nóng)業(yè)生產(chǎn)的農(nóng)民進(jìn)行操作。試驗(yàn)雨強(qiáng)為60mm／h，坡度為20°，設(shè)置3個(gè)重復(fù)，進(jìn)行連續(xù)分段（雨前BR、濺蝕階段SpE、面蝕階段ShE、細(xì)溝侵蝕階段RE）人工模擬降雨實(shí)驗(yàn)［14］。同時(shí)，采用高分辨率、非接觸式激光掃描儀［15］（徠佧，垂直誤差小于3 mm）分別測(cè)量分段降雨前后的耕作地表糙度（采樣間隔2cm×2cm），進(jìn)而構(gòu)建微地表數(shù)字高程模型（M－DEM）［15］。

1.2 數(shù)據(jù)分析

運(yùn)用SPSS 17.0分別對(duì)不同侵蝕階段坡面相對(duì)高程數(shù)據(jù)進(jìn)行描述性統(tǒng)計(jì)分析，并采用單樣本K－S（柯?tīng)柲缏宸颉鼓β宸颍┓ㄔ讦粒?.05顯著水平上對(duì)統(tǒng)計(jì)量進(jìn)行正態(tài)分布檢驗(yàn)。將服從正態(tài)分布的原始數(shù)據(jù)在GS＋軟件中進(jìn)行地統(tǒng)計(jì)分析，擬合半方差函數(shù)模型、繪制半方差函數(shù)圖，進(jìn)行空間異質(zhì)性分析。半方差函數(shù)計(jì)算公式如下:

式中:γ（h）——半方差函數(shù)；h——樣點(diǎn)間距即步長(zhǎng)（lag）；N （h）——間距為h的樣點(diǎn)數(shù)；Z（xi），Z（xi＋h）——高程Z（x）在空間位置xi和xi＋h處的觀測(cè)值。

半變異函數(shù)曲線通?？捎玫臄M合模型有球狀模型（Spherical）、高斯模型（Gaussian）、指數(shù)模型（Exponential）和線性模型（Linear）等。模型的最優(yōu)化選擇可用決定系數(shù)R2、殘差平方和RSS及變程來(lái)決定。R2越大，擬合曲線模型越精確；RSS越小，說(shuō)明擬合曲線模型的效果越好［16］，越能反映地表糙度的空間變異性。通過(guò)半方差函數(shù)擬合參數(shù)可以分析地表糙度的空間異質(zhì)性特征。

分形維數(shù)D的大小是變量復(fù)雜程度的一種度量［17］，其計(jì)算公式為:

式中:m——雙對(duì)數(shù)半方差圖的斜率，可通過(guò)對(duì)數(shù)據(jù)點(diǎn)的直線擬合來(lái)得到。D值的大小反映了變量的小尺度變異和大尺度變異的相對(duì)重要性。D值越大，說(shuō)明小尺度上的變異顯著，由隨機(jī)性因素引起的異質(zhì)性占有較大的比重，它所表現(xiàn)的空間分布越復(fù)雜；D值越小，說(shuō)明大尺度上的變異控制著地表糙度的配置，空間分布相對(duì)簡(jiǎn)單，空間自相關(guān)性較強(qiáng)，結(jié)構(gòu)性好。

Krigring插值方法是地統(tǒng)計(jì)學(xué)中最常用的插值方法之一，它以空間自相關(guān)為基礎(chǔ)，利用原始數(shù)據(jù)和半方差函數(shù)的結(jié)構(gòu)性，對(duì)坡面內(nèi)未知高程點(diǎn)的相對(duì)高程進(jìn)行最優(yōu)無(wú)偏估值，得到相關(guān)坡面地表糙度的空間變異分布狀況［18］。

2 結(jié)果與分析

2.1 地表糙度的描述性統(tǒng)計(jì)特征

運(yùn)用經(jīng)典統(tǒng)計(jì)方法對(duì)各個(gè)坡面相對(duì)高程進(jìn)行描述性統(tǒng)計(jì)特征分析，獲取變量的平均值、變異系數(shù)等，以此來(lái)衡量不同侵蝕階段地表糙度的平均狀況和變異程度，當(dāng)CV≤10%時(shí)，為弱變異性；10%＜CV＜100%時(shí)，為中等變異性；當(dāng)CV≥100%時(shí)，為強(qiáng)變異性［19］。結(jié)果如表1所示，其中，平均值反映變量的絕對(duì)變異，變異系數(shù)反映變量的相對(duì)變異。

表1 黃土坡耕地地表糙度的基本特征統(tǒng)計(jì)值

從平均值來(lái)看，表現(xiàn)為ShE＞SpE＞BR＞RE。從變異系數(shù)CV來(lái)看，各侵蝕階段的變異系數(shù)CV均在10%以下，屬于弱變異性，其中，RE階段最小，為6.4%，BR階段最大，為6.5%，SpE和ShE階段相等，介于BR和RE之間。經(jīng)K－S檢驗(yàn)，各侵蝕階段相對(duì)高程數(shù)據(jù)均符合正態(tài)分布。

2.2 地表糙度的空間結(jié)構(gòu)特征

由上述描述性統(tǒng)計(jì)特征分析可知，各侵蝕階段坡面相對(duì)高程數(shù)據(jù)均符合正態(tài)分布，故可對(duì)其進(jìn)行半方差函數(shù)分析［20］（圖1）。

圖1 不同侵蝕階段地表糙度的半方差函數(shù)

由圖1可知，在整體上，各侵蝕階段的半方差函數(shù)值均隨采樣間距的增加而增加，在一定的尺度范圍內(nèi)又保持相對(duì)穩(wěn)定的趨勢(shì)。其中，BR、SpE、ShE階段在一定的尺度范圍內(nèi)表現(xiàn)出相似的空間變異特征。

在0—10cm的較小范圍內(nèi)，各侵蝕階段的半方差值均隨采樣間距的增加而增加，表現(xiàn)為BR＞SpE＞ShE＞RE。在20—50cm的范圍內(nèi)，各侵蝕階段的半方差值保持相對(duì)穩(wěn)定的趨勢(shì)。BR、SpE、ShE在60—170cm的范圍內(nèi)又逐漸增加，達(dá)到最大值，在170—180cm的范圍內(nèi)逐漸減?。籖E階段在60—160 cm的范圍內(nèi)逐漸增加，達(dá)到最大，在160—180cm的范圍內(nèi)半方差值又隨采樣間距的增加而減小。半方差函數(shù)圖的這種變化說(shuō)明在不同侵蝕階段，黃土坡耕地地表糙度的空間異質(zhì)性存在多尺度變異特征，且不同尺度上的變異疊加在一起，形成“巢狀”的方差結(jié)構(gòu)［21］。

由表2半方差函數(shù)擬合結(jié)果可以看出，BR和SpE階段的空間變異以指數(shù)模型擬合效果較好，ShE階段以球狀模型擬合效果較好，RE階段以高斯模型擬合效果較好，其中，除了BR階段的R2較小，RSS較大之外，其他階段的理論模型擬合程度均較高，能夠很好地反映地表糙度的空間結(jié)構(gòu)特征。

塊金值C0、基臺(tái)值C＋C0、結(jié)構(gòu)比C／（C＋C0）以及變程是半方差函數(shù)中的幾個(gè)重要參數(shù)，通過(guò)對(duì)這些參數(shù)的分析可以比較地表糙度的空間異質(zhì)性特征，其中，塊金值C0反映地表糙度由隨機(jī)因素引起的空間異質(zhì)性，即受不確定性因素影響的程度。不同侵蝕階段地表糙度隨機(jī)部分的空間異質(zhì)性相同?；_(tái)值（C＋C0）反映地表糙度在研究范圍內(nèi)的總的空間變異強(qiáng)度，其值越高，表示系統(tǒng)總的空間異質(zhì)性越高。由表2得出，在RE階段地表糙度的空間異質(zhì)性最強(qiáng)。結(jié)構(gòu)比C／（C＋C0）用來(lái)量化地表糙度由自相關(guān)部分引起的空間異質(zhì)性程度的高低［22］。當(dāng)C／（C＋C0）小于25%時(shí)，表明變量具有較弱的空間自相關(guān)性，當(dāng)C／（C＋C0）在25%～75%之間時(shí)，變量具有中等的空間自相關(guān)性，當(dāng)C／（C＋C0）大于75%時(shí)，變量具有較強(qiáng)的空間自相關(guān)性［23］。通過(guò)對(duì)比分析不同侵蝕階段地表糙度的結(jié)構(gòu)比C／（C＋C0）發(fā)現(xiàn)，各侵蝕階段的C／（C＋C0）均在25%～75%之間，說(shuō)明黃土人工鋤耕坡面地表糙度在各侵蝕階段均具有中等程度的空間自相關(guān)性。變程反映地表糙度空間自相關(guān)變異的尺度范圍，在變程內(nèi)地表糙度存在空間相關(guān)性，反之則不存在。隨降雨侵蝕的不斷進(jìn)行，各侵蝕階段的空間自相關(guān)范圍也不同，ShE階段的變程最大，達(dá)到1.29m，表明在面蝕階段，坡面具有較好的空間連續(xù)性，在2cm×2cm的研究尺度下其空間變異性并不明顯，有待更小尺度的研究；RE階段的變程較小，表明在細(xì)溝侵蝕階段，影響坡面空間變異過(guò)程的作用尺度較?。?4］。SpE階段的變程最小，為0.28m，說(shuō)明濺蝕階段坡面的空間連續(xù)性較差，在以后的研究中應(yīng)取更小的采樣間距。

分形維數(shù)D表征地表糙度空間分布的復(fù)雜程度。由表2可以看出，在不同的侵蝕階段，地表糙度的分形維數(shù)較大（1.80±0.02），接近于2，這個(gè)較大的D值說(shuō)明黃土人工鋤耕坡面地表糙度在不同水蝕過(guò)程中的空間變異主要發(fā)生在較小的尺度上［21］。其中，不同侵蝕階段地表糙度的復(fù)雜性表現(xiàn)為BR≈ShE＞SpE＞RE，表明在BR和ShE階段，地表糙度的空間格局較復(fù)雜，而在SpE和RE階段地表糙度的空間格局相對(duì)簡(jiǎn)單，結(jié)構(gòu)性較好。

表2 黃土坡耕地地表糙度的半方差函數(shù)擬合曲線模型及參數(shù)

Christien等［25］認(rèn)為一個(gè)具有異質(zhì)性的系統(tǒng)屬性其在空間上的表現(xiàn)形式為斑塊狀態(tài)（Patchiness），而非均勻（Uniform）或隨機(jī)（Random）分布。由表1、表2可得出，黃土坡耕地地表糙度在不同水蝕過(guò)程中均存在空間異質(zhì)性現(xiàn)象，表明將地統(tǒng)計(jì)學(xué)與GIS技術(shù)相結(jié)合來(lái)研究黃土坡耕地地表糙度的空間異質(zhì)性具有可行性，這為今后在黃土區(qū)進(jìn)行地表糙度空間異質(zhì)性的研究提供一定的方法支持。

2.3 地表糙度的空間分布特征

Krigring插值圖能夠直觀反映地表糙度的空間異質(zhì)性及空間分布格局。利用Krigring插值法對(duì)所有坡面相對(duì)高程進(jìn)行最優(yōu)化內(nèi)插，繪制黃土人工鋤耕坡面地表糙度在各侵蝕階段的空間分布圖（圖2）并進(jìn)行分析，進(jìn)一步了解不同水蝕過(guò)程中地表糙度的空間分布特征。Krigring插值各參數(shù)的設(shè)置參照表2中的半方差函數(shù)模型參數(shù)。

圖2 不同侵蝕階段地表糙度的空間分布

由圖2可以看出，不同水蝕階段地表糙度空間分布的復(fù)雜程度和斑塊等級(jí)不同，具有明顯的斑塊狀或連續(xù)分布的特點(diǎn)，通過(guò)分形維數(shù)可以衡量斑塊的復(fù)雜程度，其中，BR、ShE階段的分形維數(shù)較大，其空間分布較復(fù)雜，在圖中呈現(xiàn)明顯的斑塊狀分布，且差異明顯；RE階段的分形維數(shù)較小，在圖中的空間分布較簡(jiǎn)單，呈現(xiàn)出較明顯的連續(xù)分布狀況。這正好和分形維數(shù)分析所得的結(jié)果相吻合。表明SpE、ShE和RE階段對(duì)地表糙度變化的影響不同，ShE階段對(duì)地表糙度變化的影響最大。

3 結(jié)論

（1）由描述性統(tǒng)計(jì)特征分析得出，黃土人工鋤耕坡面地表糙度在各侵蝕階段均具有較弱的空間變異性，這與趙龍山等［9］的研究結(jié)果相似。

（2）從空間結(jié)構(gòu)特征來(lái)看，黃土人工鋤耕坡面地表糙度在各侵蝕階段均具有中等程度的空間自相關(guān)性，空間異質(zhì)性受隨機(jī)性因素和結(jié)構(gòu)性因素的共同作用，其自相關(guān)距分別為0.46m，0.28m，1.29m，1.04m；地表糙度在RE階段具有較強(qiáng)的空間異質(zhì)性。分形維數(shù)分析表明地表糙度在不同水蝕過(guò)程中均具有較大的分形維數(shù)（1.80±0.02），表明不同侵蝕階段地表糙度的空間變異主要發(fā)生在較小尺度上，其中，不同侵蝕階段地表糙度的復(fù)雜性表現(xiàn)為BR≈ShE＞SpE＞RE，表明在BR和ShE階段，地表糙度的空間格局較復(fù)雜，而在SpE和RE階段，地表糙度的空間格局相對(duì)簡(jiǎn)單，結(jié)構(gòu)性較好。

（3）從空間分布特征來(lái)看，BR、ShE階段地表糙度的空間分布較復(fù)雜，在圖中呈現(xiàn)明顯的斑塊狀分布，且差異明顯；SpE階段地表糙度的空間分布較簡(jiǎn)單，RE階段糙度的空間分布最簡(jiǎn)單，呈現(xiàn)出較明顯的連續(xù)分布狀況。這正好和分形維數(shù)分析所得的結(jié)果相吻合。表明SpE、ShE和RE對(duì)地表糙度變化的影響不同，ShE階段對(duì)地表糙度變化的影響最大。

地表糙度受人為耕作管理措施［10］、降雨［26］、雨強(qiáng)和坡度［27］等多種因素的影響，本文僅對(duì)固定坡度（20°）、雨強(qiáng)（60mm／h）及耕作管理措施（人工鋤耕）坡面地表糙度的空間異質(zhì)性進(jìn)行了研究，對(duì)于不同坡度、雨強(qiáng)及耕作措施條件下地表糙度的空間異質(zhì)性還有待深入系統(tǒng)地研究，另外，本文對(duì)于分段降雨中各侵蝕階段的界定主要依據(jù)個(gè)人的主觀判斷，存在一定的不科學(xué)性，在以后的研究中應(yīng)尋求一種較為客觀的判定方法。另外，從地表糙度的空間結(jié)構(gòu)特征和空間分布特征得出，SpE、ShE和RE對(duì)地表糙度變化的影響不同，ShE階段對(duì)地表糙度變化的影響最大，但不同水蝕過(guò)程對(duì)地表糙度變化的影響程度如何，如何度量，本文對(duì)這一問(wèn)題的研究還不夠深入，在以后的研究中應(yīng)努力解決這一問(wèn)題。

［1］鄭子成，何淑勤，吳發(fā)啟.坡面水蝕過(guò)程中地表糙度的研究進(jìn)展［J］.節(jié)水灌溉，2008（8）:8－12.

［2］Burwell R E，Larson W E.Infiltration as influenced by tillage－induced random roughness and pore space［J］.Soil Science Society of America Journal，1969，33（3）:449－452.

［3］吳發(fā)啟，趙曉光，劉秉正.緩坡耕地侵蝕環(huán)境及動(dòng)力機(jī)制分析［M］.西安:陜西科學(xué)技術(shù)出版社，2000.

［4］Takken I，Govers G，Jetten V，et al.Effects of tillage on runoff and erosion patterns［J］.Soil and Tillage Research，2001，61（1）:55－60.

［5］Martin Y，Valeo C，Tait M.Centimetre－scale digital representations of terrain and impacts on depression storage and runoff［J］.Catena，2008，75（2）:223－233.

［6］蘇松錦，劉金福，何中聲，等.格氏栲天然林土壤養(yǎng)分空間異質(zhì)性［J］.生態(tài)學(xué)報(bào)，2012，32（18）:5673－5682.

［7］趙安玖，胡庭興，陳小紅.森林類型對(duì)土壤表層有機(jī)碳空間異質(zhì)性的影響［J］.自然資源學(xué)報(bào)，2009，24（10）:1748－1756.

［8］郭旭東，傅伯杰，馬克明，等.基于GIS和地統(tǒng)計(jì)學(xué)的土壤養(yǎng)分空間變異特征研究:以河北省遵化市為例［J］.應(yīng)用生態(tài)學(xué)報(bào)，2000，11（4）:557－563.

［9］趙龍山，宋向陽(yáng)，張青峰，等.黃土坡耕地地表微地形空間變異性研究［J］.干旱地區(qū)農(nóng)業(yè)研究，2011，29（1）:146－150.

［10］高樹(shù)靜，趙龍山，梁心藍(lán)，等.人工降雨條件下坡耕地地表糙度的時(shí)空變異分布研究［J］.水土保持研究，2010，17（3）:12－16.

［11］王軍，傅伯杰，邱揚(yáng)，等.黃土高原小流域土壤養(yǎng)分的空間異質(zhì)性［J］.生態(tài)學(xué)報(bào)，2002，22（8）:1173－1178.

［12］趙軍，張久明，孟凱，等.地統(tǒng)計(jì)學(xué)及在黑土區(qū)域土壤養(yǎng)分空間異質(zhì)性分析中的應(yīng)用:以海倫市為例［J］.水土保持學(xué)報(bào)，2004，18（6）:53－57.

［13］郭旭東，傅伯杰，馬克明，等.基于GIS和地統(tǒng)計(jì)學(xué)的土壤養(yǎng)分空間變異特征研究:以河北省遵化市為例［J］.應(yīng)用生態(tài)學(xué)報(bào)，2000，11（4）:557－563.

［14］張莉，張青峰，鄭子成，等.基于 M－DEM的黃土人工鋤耕坡面水系分維特征研究［J］.水土保持研究，2012，19（5）:7－11.

［15］趙龍山，張青峰，梁心藍(lán)，等.基于GIS的坡耕地?cái)?shù)字高程模型的建立與應(yīng)用［J］.農(nóng)業(yè)工程學(xué)報(bào)，2010，26（11）:317－322.

［16］Si J H，F(xiàn)eng Q，Yu T，et al.Spatial heterogeneity of soil nutrients in Ejina oasis［J］.Chinese Journal of E－cology，2009，28（12）:2600－2606.

［17］李小昱，雷延武，王為.農(nóng)田土壤特性的空間變異性及分形特征［J］.干旱地區(qū)農(nóng)業(yè)研究，2000，18（4）:61－65.

［18］薛冬，姚槐應(yīng)，何振立，等.紅壤酶活性與肥力的關(guān)系［J］.應(yīng)用生態(tài)學(xué)報(bào)，2005，16（8）:1455－1458.

［19］鄭紀(jì)勇，邵明安，張興昌.黃土區(qū)坡面表層土壤容重和飽和導(dǎo)水率空間變異特征［J］.水土保持學(xué)報(bào)，2004，18（4）:53－56.

［20］Vieira，S R Hatfield J L，Nielsen D R.Geostatistical theory and application to variability of some agronomical properties［J］.Hilgardia，1982，51（3）:74－75.

［21］陳玉福，于飛海，董鳴.毛烏素沙地沙生半灌木群落的空間異質(zhì)性［J］.生態(tài)學(xué)報(bào)，2000，20（4）:568－572.

［22］王存國(guó)，韓士杰，張軍輝，等.長(zhǎng)白山闊葉紅松林表層土壤水分空間異質(zhì)性的地統(tǒng)計(jì)學(xué)分析［J］.應(yīng)用生態(tài)學(xué)報(bào)，2010，21（4）:849－855.

［23］劉春利，胡偉，賈宏福，等.黃土高原水蝕風(fēng)蝕交錯(cuò)區(qū)坡地土壤剖面飽和導(dǎo)水率空間異質(zhì)性［J］.生態(tài)學(xué)報(bào)，2012，32（4）:1211－1219.

［24］張偉，陳洪松，王克林，等.典型喀斯特峰叢洼地坡面土壤養(yǎng)分空間變異性研究［J］.農(nóng)業(yè)工程學(xué)報(bào)，2008，24（1）:68－73.

［25］Ettema C H，Wardle D A.Spatial soil ecology［J］.Trends in Ecology ＆ Evolution，2002，17（4）:177－183.

［26］鄭子成，何淑勤，吳發(fā)啟.降雨對(duì)地表糙度影響的研究［J］.水土保持研究，2003，10（2）:151－154.

［27］張霞，賈蓮蓮，李占斌，等.模擬降雨條件下地表糙度動(dòng)態(tài)變化特征研究［J］.水土保持研究，2012，19（4）:16－20.