薛永偉+王茹

摘要：利用侵蝕-生產(chǎn)力評價(jià)模型（erosion-productivity impact calculator，簡稱EPIC）模型對紫色土小流域土壤可蝕性K值進(jìn)行估算，研究紫色土丘陵區(qū)小流域的土壤可蝕性K值空間變異特征。結(jié)果表明：（1）土壤顆粒組成的變化是影響紫色土小流域土壤可蝕性大小的主要因素，沙粒含量越高，土壤可蝕性K值愈大，粉粒、黏粒、有機(jī)質(zhì)的含量越高，土壤可蝕性K值愈?。唬?）流域中6種不同土地利用條件下可蝕性指標(biāo)K均值排序?yàn)榛牟莸?柑橘林地>馬尾松林地>槐樹林地>甘薯地>玉米地；（3）研究區(qū)紫色土小流域土壤可蝕性K值總體分布趨勢是自東向西隨著流域海拔的降低而降低，呈條帶狀分布，南北兩側(cè)高，中間低。隨著海拔的降低，除槐樹林外，5種地類的可蝕性指標(biāo)K值大致呈減小趨勢，表現(xiàn)為K上游>K中游>K下游；在槐樹林地則表現(xiàn)出相反的趨勢，表現(xiàn)為K上游

關(guān)鍵詞：土壤可侵性；紫色土丘陵區(qū)；小流域；EPIC模型；土地利用條件；土壤侵蝕預(yù)報(bào)；可蝕性K值

中圖分類號： S157 文獻(xiàn)標(biāo)志碼： A

文章編號：1002-1302（2016）11-0488-04

土壤可蝕性K值的大小表示土壤被沖蝕的難易程度，是影響土壤流失量的內(nèi)在因素，它也是定量研究土壤侵蝕的基礎(chǔ)[1]。土壤可蝕性是定量計(jì)算土壤流失的重要指標(biāo)，是土壤侵蝕預(yù)報(bào)模型中的必要參數(shù)[2]。1990年Wischmeier建立了侵蝕-生產(chǎn)力評價(jià)模型（erosion-productivity impact calculator，簡稱EPIC），簡化了可蝕性K值的計(jì)算[3]。我國學(xué)者研究表明，土壤可蝕性K值與土壤的理化性質(zhì)有著密切的聯(lián)系[4-7]。

我國學(xué)者通過對東北黑土地[8-10]、黃土高原地區(qū)[11-13]、南方丘陵山區(qū)[14-16]土壤可蝕性的研究，總結(jié)上述地區(qū)流域可蝕性因子同可蝕性K值的關(guān)系。雖然在紫色土區(qū)也有相關(guān)的土壤可蝕性研究，但是主要集中在人工模擬降水條件下小區(qū)尺度的土壤可蝕性變化規(guī)律及其影響因素和大空間尺度的可蝕性空間變異[17-19]。本研究以戴家溝小流域?yàn)樵囼?yàn)地點(diǎn)，通過采用可蝕性計(jì)算模型（EPIC模型），探討紫色土小流域中不同土地利用類型方式和空間因素對土壤可蝕性K值的影響規(guī)律，其目的是為了更好地了解小流域尺度土壤可蝕性空間分布特征，為小流域綜合治理中的土壤侵蝕預(yù)報(bào)和水土保持規(guī)劃提供科學(xué)依據(jù)。

1 材料與方法

1.1 研究區(qū)概況

戴家溝小流域位于重慶市萬州區(qū)長嶺鎮(zhèn)五橋河流域，閉合面積為0.67 km2，溝口地理位置為108°30′19″E、30°44′21″N，距匯入長江的出口15.45 km；地貌以淺、中丘為主，海拔在427.9～811.3 m之間，年均氣溫17 ℃，年均降水量 1 100 mm；出露侏羅系沙溪廟組紫色砂泥巖地層，巖層傾角15°～30°，土壤為灰棕紫泥，保水抗旱抗侵蝕能力較弱。流域內(nèi)分布著以馬尾松（Pinus massoniana Lamb）、洋槐（Robinia pseudoacacia L.）為主的林地，以柑橘（Citrus reticulata Banco）為主的經(jīng)果林，以玉米（Zea mays L.）、甘薯（Lpomoea batatas）為主的農(nóng)地、荒草地。流域內(nèi)多陡坡，25°以上的陡坡地面積占流域總面積的55.7%，可耕作土地比重低，豐富的降水和易流失的土壤為水土流失提供了動力。

1.2 土壤樣品采集及處理方法

結(jié)合1 ∶10 000地形圖，在實(shí)地調(diào)查繪制戴家溝小流域土地利用類型的基礎(chǔ)上，選擇林地（馬尾松林地、洋槐林地和柑橘林地）、農(nóng)地（玉米地、甘薯地）、荒草地3種土地利用類型為研究對象，主要采集小流域內(nèi)上、中、下游不同坡位、不同土地利用中0～15 cm表層土壤樣品。根據(jù)戴家溝小流域內(nèi)不同土地利用類型面積的大小，并參考張金池等在流域中土樣的采集方法[16]，首先選定1個(gè)樣點(diǎn)，以所選樣點(diǎn)為中心、10～15 m 為半徑范圍內(nèi)隨機(jī)采集3個(gè)表層土樣，共采集54個(gè)土壤樣品。土壤顆粒組成采用英國馬爾文公司生產(chǎn)的Mastersize2000激光粒度儀進(jìn)行測試，該型號儀器所測定的顆粒組成是按照美國土壤粒級制的分類標(biāo)準(zhǔn)，計(jì)算出沙粒（0.05 mm<粒徑≤2.00 mm）、粉粒（0.002 mm<粒徑≤0.05 mm）、黏粒（≤0.002 mm）3個(gè)粒徑的體積分?jǐn)?shù)；土壤有機(jī)質(zhì)含量采用重鉻酸鉀-硫酸法測定[20]，且土壤有機(jī)質(zhì)含量為土壤有機(jī)碳含量的1.724倍。

1.3 土壤可蝕性因子K值的計(jì)算

從圖1、表1可以看出，土壤可蝕性K值與沙粒含量呈極顯著正相關(guān)，隨著沙粒含量的增加，土壤可蝕性K值大致呈線性上升趨勢；而土壤可蝕性K值與粉粒、黏粒含量呈極顯著負(fù)相關(guān)，隨著粉粒、黏粒含量的增加，土壤可蝕性K值呈線性下降趨勢，與黏粒含量相比，粉粒含量與土壤可蝕性關(guān)系更密切。由圖1可見，隨著有機(jī)質(zhì)含量的增加，土壤可蝕性K值相應(yīng)下降，但相關(guān)性較低。因此可以認(rèn)為，土壤顆粒組成的變化是影響紫色土小流域土壤可蝕性差異的主要因素，且土壤可蝕性的大小主要與土壤中沙粒、粉粒含量密切相關(guān)，沙粒含量越高，越易發(fā)生侵蝕，土壤可蝕性K值愈大。

沙粒含量越高，土壤可蝕性K值越大，可能原因是沙粒含量越大，則土壤顆粒中的粉粒、黏粒含量越低。研究表明，土壤中的有機(jī)質(zhì)是土壤水穩(wěn)性結(jié)構(gòu)的膠結(jié)劑，能增加土壤疏松度、通透性，可增強(qiáng)團(tuán)聚性，從而增強(qiáng)土壤的抗侵蝕能力。而有機(jī)質(zhì)主要吸附于土壤顆粒中的黏粒上，因此黏粒含量越低，土壤中的有機(jī)質(zhì)含量也越低，表現(xiàn)為沙粒含量增高，土壤容易被侵蝕。

2.2 土地利用對土壤可蝕性K值的影響特征

研究區(qū)內(nèi)荒草地和馬尾松林主要分布在流域南北兩側(cè)的山頂，柑橘林地、槐樹林地主要分布在中坡位置，下坡位置主要分布農(nóng)耕地。根據(jù)EPIC模型，按經(jīng)典統(tǒng)計(jì)方法分析，分別計(jì)算出6種不同土地利用類型土壤的可蝕性指標(biāo)K值，其計(jì)算差異性分析結(jié)果見圖2。

通過采用EPIC模型計(jì)算小流域中6種不同土地利用類型表層土壤的可蝕性K值，得出流域中可蝕性指標(biāo)K均值從大到小排序?yàn)榛牟莸?柑橘林地>馬尾松林地>槐樹林地>甘薯地>玉米地。從圖2可以看出，在上游位置，6種土地利用類型中槐樹林地的可蝕性K值最低，同時(shí)玉米地的可蝕性K值要大于甘薯地。進(jìn)一步比較相同流域不同土地利用類型的土壤可蝕性差異的顯著性，由圖2可以看出，玉米地、甘薯地的土壤可蝕性差異不顯著，它們與其他4種土地利用類型的土壤可蝕性表現(xiàn)出顯著性差異（P<0.05），而荒草地、柑橘林地、槐樹林地、馬尾松林地基本上與其他5種土地利用類型的土壤可蝕性差異顯著（P<0.05），僅在流域的中游位置，槐樹林地、馬尾松林地的土壤可蝕性表現(xiàn)出差異不顯著現(xiàn)象，這有待進(jìn)一步研究。

通過分析不同地類土壤表層的理化性質(zhì)發(fā)現(xiàn)：（1）荒草地的沙粒含量在6種地類中最高，而有機(jī)質(zhì)含量則最低。這是由于荒草地在流域中處于母質(zhì)出露較多、土壤發(fā)育未完全成熟、土層薄、養(yǎng)分比較貧瘠的坡頂位置，同時(shí)荒草地的覆蓋度低、根系淺，在雨滴的動能打擊下，土壤顆粒被破壞，在重力作用下向下坡遷移，因此荒草地表層土沙粒含量增多，而粉粒、黏粒含量減少，可見荒草地對于土壤侵蝕的敏感程度最高。

（2）流域中柑橘林地的土壤可蝕性K值較高，高于馬尾松林地、槐樹林地。主要原因是土壤可蝕性K值還與土壤中的有機(jī)質(zhì)積累量有關(guān)，研究表明，植被能改善土壤性質(zhì)、增加有機(jī)質(zhì)含量，使土壤可蝕性K值降低[19]。橘園主要由人為管理，林下多數(shù)為無覆蓋狀態(tài)，或者種有甘薯，因此橘園林地有機(jī)質(zhì)積累較馬尾松林地、槐樹林地少。同時(shí)，由于橘園的郁閉度相對較低，發(fā)生侵蝕性降水時(shí)，林下土壤細(xì)顆粒容易被沖走，使得土壤的可蝕性要高于馬尾松林地、槐樹林地。通過對比機(jī)械組成變化和有機(jī)質(zhì)含量也可以發(fā)現(xiàn)：與荒草地相比，柑橘林地的沙粒含量較低，而粉粒、黏粒、有機(jī)質(zhì)含量則較高，但是與馬尾松林地、槐樹林地相比，柑橘林的沙粒含量明顯高于馬尾松林、槐樹林地，而柑橘林的有機(jī)質(zhì)含量要明顯低于另外2種林地。

（3）馬尾松林地、槐樹林地的可蝕性K值總體呈高于玉米地和甘薯地的趨勢。這是因?yàn)殡m然2種林地的有機(jī)質(zhì)含量高于農(nóng)耕地，但是由于玉米地、甘薯地主要分布在農(nóng)戶附近，主要為梯平地，與其他地類相比，受到人為因素的影響相對較大，人為活動對地表的擾動使土壤大顆粒物質(zhì)含量減少，粉粒、黏粒含量相對增加。因此，侵蝕性因子對玉米地和甘薯地的影響要小于馬尾松林地、槐樹林地。

（4）上游、中游馬尾松林地的可蝕性K值大于槐樹林地。主要是因?yàn)轳R尾松林地的郁閉度較小。遇到降水時(shí)，馬尾松林的枝葉降低雨滴動能的效果明顯要弱于槐樹林的枝葉，而馬尾松林地覆蓋度較低，雨滴的擊濺作用破壞了馬尾松林地土壤表層的結(jié)構(gòu)；同時(shí)，馬尾松林地主要生長于流域的山頂，土層較薄，沙粒含量要高于槐樹林地，且枯枝落葉難以分解，有機(jī)質(zhì)含量較槐樹林地低，因此馬尾松林地表現(xiàn)出比槐樹林地更容易受到侵蝕的特征。

（5）在流域的中游位置馬尾松林地和槐樹林地的可蝕性差異不顯著。發(fā)生異變主要是由于這一流域位置的馬尾松林地、槐樹林地的沙粒、黏粒含量差異不明顯。雖然馬尾松林地土壤可蝕性的平均值大于槐樹林地，但是由于這一流域位置的馬尾松林地、槐樹林地所有樣本測定的土壤沙粒、黏粒含量差異不明顯，因而可蝕性K值差異也不顯著。因此馬尾松林地和槐樹林地在流域上游、下游的可蝕性差異顯著，流域中游的可蝕性的差異反而不顯著。

（6）整個(gè)流域中玉米地和甘薯地的可蝕性差異性不顯著。這可能是由于2種地類都受到人為因素的影響，2種作物的管理措施差異性較小，且土壤的沙粒、黏粒含量差異不大，因而其可蝕性差異性不顯著。

2.3 土壤可蝕性K值的空間響應(yīng)特征

研究區(qū)海拔自東向西逐漸降低、南北兩側(cè)高，中間低，研究區(qū)內(nèi)土壤可蝕性K值總體分布趨勢是自東向西隨著流域海拔的降低而降低，呈條帶狀分布，南北兩側(cè)高、中間低。通過流域的分析可知，隨著海拔的降低，除槐樹林地外5種地類的可蝕性指標(biāo)K值大致呈減小趨勢，排序?yàn)镵上游>K中游>K下游，槐樹林地的土壤可蝕性則表現(xiàn)出相反的趨勢，排序?yàn)镵上游

從上游到下游，隨著海拔的降低，可蝕性指標(biāo)K值大致呈減小趨勢，根據(jù)陳利頂?shù)鹊难芯縖21]，上游位置為整個(gè)流域的“源”景觀，屬于能促進(jìn)土壤侵蝕過程發(fā)展的景觀類型，在降水侵蝕、徑流沖刷的雙重作用下，上游位置的地表結(jié)構(gòu)被破壞后，土壤表層中細(xì)顆粒的粉粒、黏粒被沖刷至流域的中、下游。此外，在流域的下游處于“匯”景觀位置，接納來自流域上游和中游遷移而來的細(xì)顆粒、養(yǎng)分，能阻止或延緩?fù)寥狼治g過程的發(fā)展。通過機(jī)械組成分析也可知，6種地類的沙粒含量隨著海拔的降低而減少，粉粒、黏粒、有機(jī)質(zhì)含量隨著海拔的降低而增加，土壤可蝕性隨著海拔降低表現(xiàn)為K上游>K中游>K下游，由此得出以下結(jié)論。

（1）柑橘林地、槐樹林地都表現(xiàn)出在不同流域位置土壤可蝕性差異顯著（P<0.05）。因?yàn)橥寥揽晌g性K值與土壤中的有機(jī)質(zhì)積累量有關(guān)，植被能改善土壤性質(zhì)，增加有機(jī)質(zhì)含量，使土壤可蝕性K值降低。而橘園主要由人為管理，林下多為無覆蓋狀態(tài)，郁閉度相對較低，發(fā)生侵蝕性降水時(shí)，林下土壤細(xì)顆粒容易被沖走，集中至下游，肥力也跟著降低，從而使土壤可蝕性從流域的上游到下游迅速減?。换睒淞值耐寥揽晌g性隨著海拔的降低反而增大是因?yàn)樵谥?、下游位置由于人為活動的影響，致使槐樹林地郁閉度低，為疏林地，土壤表層結(jié)構(gòu)較差，隨著海拔的降低，沙粒含量反而增加。所以槐樹林在上游位置的土壤表層有機(jī)質(zhì)含量要遠(yuǎn)遠(yuǎn)高于流域中、下游位置。（2）玉米地、甘薯地的土壤可蝕性從流域的上游到下游逐漸減小，但是由于它們的耕地類型主要為梯平地，且受到人工施肥、翻耕等管理因素的影響較大，因此受到流域位置“源”“匯”景觀生態(tài)學(xué)的影響較小，2種農(nóng)耕地的土壤可蝕性在相同流域中表現(xiàn)為差異不顯著。（3）雖然荒草地、馬尾松林都處于流域的上坡位置，但是由于地面覆蓋度以及根系的不同，且在上游位置由于荒草地的覆蓋度低、荒草的淺根性特征，使固土效果較差，土壤母質(zhì)還在發(fā)育，因此水土流失現(xiàn)象嚴(yán)重，在中游、下游位置，都有細(xì)粒的泥沙沉積、養(yǎng)分的積累，因此荒草地上游位置與中、下游位置的土壤可蝕性差異顯著（P<0.05），而中、下游位置的土壤可蝕性差異不顯著。（4）馬尾松林在中游位置，土層相對于下游位置較薄，地表覆蓋情況較上游變化不大，表現(xiàn)為土壤表層的沙粒含量變化差異不明顯，因此上游、中游的馬尾松林地土壤可蝕性差異不明顯，而在下游位置由于“匯”景觀的影響，土壤中的細(xì)顆粒物質(zhì)大量沉積，同時(shí)馬尾松林地的枯枝落葉層能夠通過腐殖質(zhì)的大量分解，使流域下游的馬尾松林地有機(jī)質(zhì)含量要大大高于上游、中游。因此，通過不同流域位置的土壤可蝕性顯著性分析可知，上游、中游的馬尾松林地土壤表層的可蝕性差異不顯著，而與下游的馬尾松林土壤表層的可蝕性差異顯著（P<0.05）。

3 結(jié)論與討論

3.1 結(jié)論

通過基于EPIC模型對6種土地利用類型的表層土壤可蝕性的計(jì)算，同時(shí)分析土壤機(jī)械組成、有機(jī)質(zhì)含量與土壤可蝕性K值的相關(guān)性得出以下結(jié)論。

（1）通過分析土壤可蝕性與土壤機(jī)械組成、有機(jī)質(zhì)含量的相關(guān)性可知，土壤顆粒組成的變化是影響紫色土小流域土壤可蝕性大小差異的主要因素，且可蝕性大小主要與土壤中沙粒、粉粒含量密切相關(guān)。沙粒含量越高，越易發(fā)生侵蝕，土壤可蝕性K值愈大；粉粒、黏粒、有機(jī)質(zhì)含量越高，越不容易發(fā)生侵蝕，土壤可蝕性K值愈小。

（2）流域中可蝕性指標(biāo)K均值從大到小依次為荒草地>柑橘林地>馬尾松林地>槐樹林地>甘薯地>玉米地。下坡位置的玉米地和甘薯地由于土地利用方式為梯平地，且受到人為因素的影響，抵抗土壤侵蝕能力較強(qiáng)，而流域中荒草地的土壤抗侵蝕能力最弱，這與張金池等得出的結(jié)論[16]剛好相反。

（3）研究區(qū)紫色土小流域土壤可蝕性K值總體分布趨勢是自東向西隨著流域海拔的降低而降低，呈條帶狀分布，南北兩側(cè)高、中間低。隨著海拔的降低（除槐樹林外），5種地類的可蝕性指標(biāo)K值大致呈減小趨勢，表現(xiàn)為K上游>K中游>K下游。在槐樹林地則表現(xiàn)出相反的趨勢，表現(xiàn)為K上游

3.2 討論

通過對比發(fā)現(xiàn)，對于不同土壤類型，其土壤可蝕性受到的土壤機(jī)械組成變化的影響也不同。黃土可蝕性的大小主要與土壤中粉粒、黏粒含量密切相關(guān)。黃土中黏粒含量愈高，土壤抵抗侵蝕的能力越強(qiáng)，可蝕性指標(biāo)值越小。粉粒含量愈高，越易發(fā)生侵蝕，土壤可蝕性指標(biāo)值愈大。張金池等認(rèn)為，土壤可蝕性與土壤顆粒組成和有機(jī)碳含量有關(guān)，植被能改善土壤性質(zhì)，增加有機(jī)碳含量，使土壤可蝕性K值降低[16]。紫色土作為一種成土作用迅速、礦物組成復(fù)雜、礦質(zhì)養(yǎng)分含量豐富，質(zhì)地松軟，易于風(fēng)化，抗蝕性較差的一種土壤類型，其可蝕性的大小主要與土壤中沙粒、粉粒含量密切相關(guān)，紫色土丘陵區(qū)不同土地利用方式土壤的可蝕性差異也比較大。通過研究可知，紫色土小流域的土壤可蝕性表現(xiàn)出不同于其他土壤類型的規(guī)律。例如，沙粒含量越高，越易發(fā)生侵蝕，土壤可蝕性K值愈大。由于人為因素的影響，梯平地中的農(nóng)耕地可蝕性K值反而要低于森林覆蓋區(qū)的，而馬尾松林地作為針葉林地的可蝕性K值要比作為闊葉林地的槐樹林地土壤可蝕性K值大。這和張金池研究的黃棕壤、紅壤小流域的可蝕性分布規(guī)律不同。

流域中的土壤可蝕性受到氣象因素、地形因素、植被類型和人為因素等多種因素的影響。本研究通過采用EPIC模型計(jì)算紫色土小流域不同土地利用方式的土壤可蝕性K值，探討植被類型以及空間分布對流域中可蝕性變化規(guī)律的影響。由于自然情況下的土壤可蝕性的影響因子復(fù)雜，目前國內(nèi)外研究學(xué)者對不同土質(zhì)的土壤可蝕性實(shí)測值主要是結(jié)合野外徑流小區(qū)觀測，采用通用土壤流失方程進(jìn)行推算，而這種土壤可蝕性估算方法需要有降水觀測資料，具有一定的限制性，不能適應(yīng)國內(nèi)復(fù)雜地形的小流域土壤可蝕性計(jì)算。目前土壤可蝕性估算方法主要有土壤可蝕性諾莫方程估算、EPIC估算、Shirazi公式估算，這些直接應(yīng)用于紫色土丘陵區(qū)土壤可蝕性的估算，需要進(jìn)行一定的修正，而侵蝕生產(chǎn)力模型EPIC需求的資料相對較少，且估算精度較高，應(yīng)用廣泛。土壤可蝕性作為評價(jià)土壤對侵蝕敏感程度的一個(gè)重要指標(biāo)，其在不同坡度和坡長條件下，對紫色土小流域土壤可蝕性大小的影響以及在紫色土丘陵區(qū)對比上述3種估算的差異性有待于更深入研究和論證。

參考文獻(xiàn)：

[1]Wischmeier W H. Predicting rainfall erosion losses：a guide to conservation planning[M]. Washington DC：USDA，Agricultural Handbook，1978：537.

[2]劉寶元，張科利，焦菊英. 土壤可蝕性及其在侵蝕預(yù)報(bào)中的應(yīng)用[J]. 自然資源學(xué)報(bào)，1999，14（4）：345-350.

[3]United States Department of Agriculture. EPIC-erosion/productivity impact calculator 1. Model documentation[Z].

[4]方學(xué)敏，萬兆惠，徐永年. 土壤抗蝕性研究現(xiàn)狀綜述[J]. 泥沙研究，1997，17（2）：87-91.

[5]宋 陽，劉連友，嚴(yán) 平，等. 土壤可蝕性研究評述[J]. 干旱區(qū)地理，2006，29（1）：124-131.

[6]田積瑩，黃義端. 子午嶺連家砭地區(qū)土壤物理性質(zhì)與土壤抗侵蝕性指標(biāo)的初步研究[J]. 土壤學(xué)報(bào)，1964，12（3）：286-296.

[7]楊玉盛. 不同利用方式下紫色土可蝕性的研究[J]. 水土保持學(xué)報(bào)，1992，6（3）：52-58.

[8]高德武. 黑龍江土壤流失方程中土壤可蝕性因子（K）的研究[J]. 國土與自然資源研究，1993（3）：40-43.

[9]范昊明，蔡強(qiáng)國，崔 明. 東北黑土漫崗區(qū)土壤侵蝕垂直分帶性研究[J]. 農(nóng)業(yè)工程學(xué)報(bào)，2005，21（6）：8-11.

[10]繆馳遠(yuǎn)，劉寶元，劉 剛，等. 東北典型黑土區(qū)剖面粒徑分布特征及其可蝕性研究[J]. 水土保持學(xué)報(bào)，2008，22（3）：20-23.

[11]朱顯謨. 黃土地區(qū)植被因素對于水土流失的影響[J]. 土壤學(xué)報(bào)，1960，8（2）：110-121.

[12]蔣定生. 黃土區(qū)不同利用類型土壤抗沖刷能力的研究[J]. 土壤通報(bào)，1979（4）：20-23.

[13]周佩華，武春龍. 黃土高原土壤抗沖性的試驗(yàn)研究方法探討[J]. 水土保持學(xué)報(bào)，1993，7（1）：29-34.

[14]邢廷炎，史學(xué)正，于東升. 我國亞熱帶土壤可蝕性的對比研究[J]. 土壤學(xué)報(bào)，1998，35（3）：296-302.

[15]張黎明. 我國南方不同類型土壤可蝕性K值及相關(guān)因子研究[D]. 海口：華南熱帶農(nóng)業(yè)大學(xué)，2005.

[16]張金池，李海東，林 杰，等. 基于小流域尺度的土壤可蝕性K值空間變異[J]. 生態(tài)學(xué)報(bào)，2008，28（5）：2199-2206.

[17]張 兵，蔣光毅，陳正發(fā)，等. 紫色丘陵區(qū)土壤可蝕性因子研究[J]. 土壤學(xué)報(bào)，2010，47（2）：354-358.

[18]侯大斌. 川渝地區(qū)土壤可蝕性評價(jià)[D]. 雅安：四川農(nóng)業(yè)大學(xué)，2001.

[19]李月臣，劉春霞，趙純勇，等. 三峽庫區(qū)（重慶段）土壤侵蝕敏感性評價(jià)及其空間分異特征[J]. 生態(tài)學(xué)報(bào)，2009，29（2）：788-796.

[20]中國科學(xué)院南京土壤研究所土壤物理研究室.土壤理化分析[M]. 上海：上?？茖W(xué)技術(shù)出版社，1980.

[21]陳利頂，傅伯杰，趙文武. “源”“匯”景觀理論及其生態(tài)學(xué)意義[J]. 生態(tài)學(xué)報(bào)，2006，26（15）：1444-1449.