劉瑛娜,劉寶元,2?,張科利,2,鄭袁志

(1.北京師范大學(xué)地表過程與資源生態(tài)國家重點(diǎn)實(shí)驗(yàn)室,100875,北京；2.北京師范大學(xué)地理學(xué)與遙感科學(xué)學(xué)院,100875,北京；3.榮豐控股集團(tuán)股份有限公司,100070,北京)

微型小區(qū)在土壤可蝕性估算中的應(yīng)用

劉瑛娜1,劉寶元1,2?,張科利1,2,鄭袁志3

(1.北京師范大學(xué)地表過程與資源生態(tài)國家重點(diǎn)實(shí)驗(yàn)室,100875,北京；2.北京師范大學(xué)地理學(xué)與遙感科學(xué)學(xué)院,100875,北京；3.榮豐控股集團(tuán)股份有限公司,100070,北京)

土壤可蝕性反映土壤對(duì)侵蝕的敏感程度,是評(píng)價(jià)土壤侵蝕程度和預(yù)報(bào)土壤流失量的重要指標(biāo)。方便可行且可信的土壤可蝕性測(cè)定與估算方法在土壤可蝕性的研究中具有重要意義。采用微型小區(qū)對(duì)黑土、褐土、黃土、紫色土和紅壤5種土壤可蝕性K值進(jìn)行估算,并基于標(biāo)準(zhǔn)小區(qū)天然降雨的徑流泥沙觀測(cè)資料,判斷微型小區(qū)能否代替標(biāo)準(zhǔn)小區(qū)進(jìn)行土壤可蝕性的估算。結(jié)果表明:微型小區(qū)估算的土壤可蝕性K值比標(biāo)準(zhǔn)小區(qū)估算K值小,但具有很好的規(guī)律性,可以進(jìn)行校正。黑土、紫色土、褐土和紅壤在天然降雨條件下微型小區(qū)計(jì)算的K值與標(biāo)準(zhǔn)小區(qū)間存在良好的線性關(guān)系,可對(duì)微型小區(qū)條件下的K值進(jìn)行修訂以達(dá)到標(biāo)準(zhǔn)小區(qū)的要求；容易產(chǎn)生細(xì)溝的黃土不宜采用微型小區(qū)進(jìn)行K值的估算。

土壤可蝕性；人工模擬降雨；微型小區(qū)；標(biāo)準(zhǔn)小區(qū)

作為反映土壤對(duì)侵蝕敏感程度的重要屬性[1-2],土壤可蝕性因子是通用土壤流失方程(USLE,Universal Soil Loss Equation)中的重要參數(shù)。根據(jù)T.C. Olson等[3]的定義,土壤可蝕性K值需要通過野外標(biāo)準(zhǔn)徑流小區(qū)觀測(cè)資料來計(jì)算。要求觀測(cè)年限長,存在數(shù)據(jù)累積慢,建設(shè)及維護(hù)費(fèi)用昂貴等問題。為此,國內(nèi)外學(xué)者一直努力通過建立土壤性質(zhì)和K值之間的經(jīng)驗(yàn)關(guān)系來解決大范圍、多類型的土壤可蝕性K值估算問題[4-5]。由于土壤特性的差異,現(xiàn)有的國外數(shù)學(xué)模型和圖解法不能直接用于我國土壤可蝕性值計(jì)算[6-7]。張科利等[4]收集全國5個(gè)主要水蝕區(qū)野外小區(qū)的實(shí)測(cè)資料對(duì)諾謨圖方法[8]、EPIC模型中的計(jì)算方法和M.A.Shirazi等[9]公式估算的土壤可蝕性K值進(jìn)行驗(yàn)證,結(jié)果表明:3種方法估算的K明顯大于野外小區(qū)實(shí)測(cè)值,不能直接用于估算我國的土壤可蝕性K值。翟偉峰等[10]通過對(duì)比實(shí)驗(yàn)小區(qū)野外實(shí)測(cè)與諾謨圖估算的土壤可蝕性K值,認(rèn)為估算值明顯小于實(shí)測(cè)值,諾謨圖法不能直接應(yīng)用于東北典型黑土區(qū)。高德武[11]采用實(shí)測(cè)法、經(jīng)驗(yàn)公式法和查圖法(諾謨圖法)對(duì)黑龍江省白漿土、黑土和暗棕壤K值進(jìn)行對(duì)比分析,結(jié)果表明查圖法和經(jīng)驗(yàn)公式法求得的K值很接近,但比實(shí)測(cè)法求得的K值偏小30%。修訂K值估算公式需要大量的標(biāo)準(zhǔn)徑流小區(qū)資料,野外觀測(cè)周期長、成本高,如采用人工模擬降雨可以在短時(shí)間內(nèi)獲得大量土壤侵蝕數(shù)據(jù)；但目前通用的人工模擬降雨機(jī)有效降雨面積有限[12],在20 m長的標(biāo)準(zhǔn)小區(qū)上進(jìn)行人工模擬降雨實(shí)驗(yàn)難度很大,成本也高,因此,很有必要探索利用微型小區(qū)來估算土壤可蝕性K值的方法。

筆者對(duì)微型小區(qū)5種土壤可蝕性K值進(jìn)行估算,以天然降雨條件下的我國5種典型土壤標(biāo)準(zhǔn)徑流小區(qū)實(shí)測(cè)數(shù)據(jù)為依據(jù),對(duì)比天然降雨條件下標(biāo)準(zhǔn)小區(qū)與微型小區(qū)K值,討論能否采用微型小區(qū)作為測(cè)定土壤可蝕性的簡易方法,為人工模擬降雨在微型小區(qū)上測(cè)定土壤可蝕性提供依據(jù)。

1 材料與方法

1.1 標(biāo)準(zhǔn)小區(qū)與微型小區(qū)天然降雨觀測(cè)

表1 5種土壤的理化性質(zhì)Tab.1 Physical and chemical properties of five types of soils

5個(gè)室外標(biāo)準(zhǔn)徑流小區(qū)和5個(gè)微型小區(qū)布設(shè)于北京師范大學(xué)房山綜合實(shí)驗(yàn)基地,均為人工填充小區(qū),模擬耕作表層土壤狀況填充我國5大水蝕區(qū)的典型土壤類型,分別為黑龍江省嫩江縣的黑土,北京市密云縣的褐土,陜西省榆林市的黃土,四川省南充市的紫色土,福建省安溪縣的紅壤。表1示出5種土壤的理化性質(zhì)。標(biāo)準(zhǔn)小區(qū)寬2.1 m,坡長20 m,坡度5°,裸露翻耕,按照USLE標(biāo)準(zhǔn)小區(qū)建設(shè)[13]。微型小區(qū)寬1 m,坡長1 m,坡度14.1°,裸露翻耕,是人工降雨試驗(yàn)常用的微型小區(qū)尺寸[14]。每個(gè)小區(qū)都配有分流桶和集流桶收集徑流和泥沙,侵蝕性降雨發(fā)生后,對(duì)各徑流小區(qū)進(jìn)行徑流量和侵蝕產(chǎn)沙量的觀測(cè)。標(biāo)準(zhǔn)小區(qū)觀測(cè)年限為2006—2011年,微型小區(qū)為2007—2011年。實(shí)驗(yàn)基地氣象園內(nèi)裝有虹吸式雨量計(jì)、自動(dòng)雨量計(jì)和普通雨量筒,用以觀測(cè)降雨并記錄降雨過程,每年的降雨特征見表2。

表2 實(shí)驗(yàn)基地2006—2011年侵蝕性降雨特征Tab.2 Erosive rainfall characteristics of the experimental base during 2006—2011

1.2 土壤可蝕性K值計(jì)算

標(biāo)準(zhǔn)小區(qū)與微型小區(qū)估算土壤可蝕性K值,都采用中國土壤侵蝕模型(CSLE)[4],CSLE的表達(dá)式為

式中:A為土壤水蝕模數(shù),t/hm2；R為降雨侵蝕力因子,MJ·mm/(hm2·h)；K為土壤可蝕性因子,t·hm2· h/(hm2·MJ·mm)；L為坡長因子；S為坡度因子；B為生物措施因子；E為工程措施因子；T為耕作措施因子。

降雨侵蝕力因子R采用經(jīng)典算法計(jì)算。

式中:E為一次降雨總動(dòng)能,MJ/hm2；I30為該次降雨的最大30 min降雨強(qiáng)度,mm/h；虹吸式雨量計(jì)詳細(xì)記錄的每次降雨過程數(shù)據(jù)用以計(jì)算R,每次侵蝕性天然降雨的R累加后作為每年的R。

坡長因子L和坡度因子S的計(jì)算公式為:

式中:θ為坡度值,(°)；λ為坡長,m；m為坡長指數(shù),根據(jù)坡度不同取值不同。

室外標(biāo)準(zhǔn)小區(qū)和微型小區(qū)都是裸地并經(jīng)過翻耕處理,因此其B、E、T取值都為1。觀測(cè)每次侵蝕性天然降雨的侵蝕產(chǎn)沙量,用以計(jì)算小區(qū)每年的土壤侵蝕模數(shù)A。

2 結(jié)果與分析

2.1 天然降雨條件下標(biāo)準(zhǔn)小區(qū)估算K值

根據(jù)房山實(shí)驗(yàn)基地標(biāo)準(zhǔn)小區(qū)2006—2011年共計(jì)61次的侵蝕泥沙與降雨特性觀測(cè)數(shù)據(jù),并采用式(2)～(5)計(jì)算標(biāo)準(zhǔn)小區(qū)5種土壤的可蝕性K值。標(biāo)準(zhǔn)小區(qū)6年實(shí)測(cè)資料估算5種土壤的K值由大到小為紫色土、黑土、紅壤、黃土和褐土(值分別為0.066、0.050、0.039、0.027、0.025 t·hm2·h/(hm2· MJ·mm))(圖1)。土壤性質(zhì)的差異造成5種土壤的可蝕性K值差異較大[15],這5種土壤中黑土與紫色土的黏粒(＜0.002 mm)質(zhì)量分?jǐn)?shù)高(表1),質(zhì)地最細(xì),更容易被侵蝕。

圖1 標(biāo)準(zhǔn)小區(qū)估算5種土壤K值年際變化Fig.1 Changes in K value estimated by unit plot of five types of soils

標(biāo)準(zhǔn)小區(qū)估算的5種土壤類型的土壤可蝕性K值大小在每年都有所波動(dòng)(圖1),其中褐土的年際變化最小(0.017 7～0.032 4 t·hm2·h/(hm2·MJ· mm)),估算的K值較集中,其他4種土壤的K值年際變化較大。由于降雨特性對(duì)土壤侵蝕量大小的影響很大,因而每年的土壤可蝕性K值因降雨不同而變化[16]。由圖2可見降雨侵蝕力越大標(biāo)準(zhǔn)小區(qū)估算的土壤可蝕性K值越大,5種土壤的K值均與降雨侵蝕力具有很好的線性相關(guān)關(guān)系:紫色土、黑土和黃土的土壤可蝕性隨降雨侵蝕力的增加的變化程度大；褐土的土壤可蝕性K值與降雨侵蝕力擬合曲線的斜率明顯要低于其他4種土壤,說明褐土的土壤可蝕性隨降雨侵蝕力增加的變化程度小,褐土的機(jī)械組成與其他4種土壤相比,黏粒、粉粒和砂粒的質(zhì)量分?jǐn)?shù)相差不大,不同降雨對(duì)褐土顆粒的分選侵蝕差異不大。土壤可蝕性是體現(xiàn)土壤抗侵蝕能力的一種特性,反映的是某一特定土壤的長期平均狀況,受降雨特性等因素年季變化的影響,某一年的實(shí)測(cè)數(shù)據(jù)無法反映土壤可蝕性K值的長期平均狀況,應(yīng)采用長期觀測(cè)小區(qū)多年總侵蝕模數(shù)與總降雨侵蝕力來計(jì)算每種土壤可蝕性K值,以消除土壤前期含水量和土壤性質(zhì)等隨季節(jié)變化的影響。

圖2 標(biāo)準(zhǔn)小區(qū)估算5種土壤K與降雨侵蝕力間的關(guān)系Fig.2 Relationship between K value estimated by unit plot and rainfall erosivity

2.2 天然降雨條件下微型小區(qū)估算K值

基于微型小區(qū)2007—2011年共計(jì)47次的侵蝕泥沙與降雨特性觀測(cè)數(shù)據(jù),采用式(2)～(5)計(jì)算微型小區(qū)5種土壤的可蝕性K值。微型小區(qū)觀測(cè)資料估算的5種土壤可蝕性K值在0.005～0.034 t· hm2·h/(hm2·MJ·mm)范圍內(nèi)變化:黃土K值最小,為0.005 t·hm2·h/(hm2·MJ·mm),與其他4種土壤K值差異很大(圖3)；紫色土、黑土K值較大,紅壤次之,褐土較小,黃土最小。

微小區(qū)天然降雨估算的K每年亦有波動(dòng)(圖3):其中黃土的 K值變化范圍較小(0.002 3～0.012 t·hm2·h/(hm2·MJ·mm))。由圖4所示,估算的黃土的K值與降雨侵蝕力沒有相關(guān)關(guān)系；其他4種土壤的K值均與降雨侵蝕力具有線性關(guān)系,降雨侵蝕力越大土壤可蝕性K值越大；紫色土、黑土隨降雨侵蝕力的增加,變化程度較大；紅壤和黃土次之。

2.3 標(biāo)準(zhǔn)小區(qū)與微型小區(qū)估算K值關(guān)系

將5種土壤在天然降雨條件下,標(biāo)準(zhǔn)小區(qū)與微型小區(qū)實(shí)測(cè)計(jì)算得到的土壤可蝕性K值進(jìn)行比較分析(圖5),結(jié)果表明5種土壤微型小區(qū)估算的土壤可蝕性K值均比標(biāo)準(zhǔn)小區(qū)估算的K值小。微型小區(qū)的坡長為1 m,標(biāo)準(zhǔn)小區(qū)的坡長為20 m,坡長變短,徑流對(duì)土壤的分離和搬運(yùn)能力降低,對(duì)土壤的侵蝕作用隨之降低。將微型小區(qū)侵蝕量采用坡長公式(3)和(4)訂正到標(biāo)準(zhǔn)小區(qū)上,計(jì)算的K值仍然小很多,因此在采用微型小區(qū)進(jìn)行土壤可蝕性K值的估算時(shí)必須進(jìn)行修訂。另一個(gè)原因可能是微小區(qū)觀測(cè)的徑流和侵蝕量有損失,因?yàn)槲⑿^(qū)尺寸小,降雨擊濺起的水滴可能會(huì)飛到小區(qū)外面,造成收集到的徑流量和泥沙量小于實(shí)際產(chǎn)流產(chǎn)沙量,在以細(xì)溝間侵蝕為主的微小區(qū)上,侵蝕量較小,擊濺損失的比例可能相當(dāng)可觀[17]。

圖3 微型小區(qū)估算5種土壤K值年際變化Fig.3 Changes in K value estimated by micro-plot of five types of soils

圖4 微型小區(qū)估算5種土壤K與降雨侵蝕力間的關(guān)系Fig.4 Relationship between K value estimated by micro-plot and rainfall erosivity

采用回歸方法分析5種土壤標(biāo)準(zhǔn)小區(qū)估算的K值(Ku)與微型小區(qū)估算的K值(Km)的關(guān)系,結(jié)果見表3。可以看出,黑土、褐土、紫色土和紅壤標(biāo)準(zhǔn)小區(qū)估算的K值與微型小區(qū)估算的K值之間存在良好的線性關(guān)系,而黃土沒有相關(guān)關(guān)系。微型小區(qū)估算的黃土K值變化范圍較小,受降雨影響很小,而標(biāo)準(zhǔn)小區(qū)估算黃土的K值隨降雨侵蝕力增加的變化程度大。標(biāo)準(zhǔn)小區(qū)黃土在侵蝕力較大的強(qiáng)降雨作用下容易產(chǎn)生細(xì)溝,使侵蝕模數(shù)明顯增大[18],而微型小區(qū)黃土由于較短坡長的影響則不容易產(chǎn)生細(xì)溝[19]。土壤的黏粒質(zhì)量分?jǐn)?shù)對(duì)土壤入滲能力有很大影響,隨著黏粒質(zhì)量分?jǐn)?shù)增大,入滲能力會(huì)遞減[20]。由于黃土顆粒組成具有較低的黏粒質(zhì)量分?jǐn)?shù)和較高的粉粒(0.002～0.05 mm)質(zhì)量分?jǐn)?shù),使得其入滲能力高,黃土產(chǎn)流量明顯小于其他4種土壤類型,再加上短坡長徑流能量降低的影響,導(dǎo)致微型小區(qū)黃土的侵蝕泥沙極少,從而計(jì)算得到的K值很低。黑土、紫色土和紅壤3種土壤的擬合直線相近,與褐土的差異較大。土壤性質(zhì)的不同,使土壤侵蝕過程和機(jī)制的變化有很大不同。

圖5 微型小區(qū)估算的K值與標(biāo)準(zhǔn)小區(qū)估算的K值間的關(guān)系Fig.5 Relationship between K value estimated by micro-plot and unit plot

表3 標(biāo)準(zhǔn)小區(qū)估算的值Ku與微型小區(qū)估算的Km值的相關(guān)性分析Tab.3 Correlation analysis between K value estimated by unit plot and micro-plot

3 結(jié)論與討論

1)采用微型小區(qū)方法比標(biāo)準(zhǔn)小區(qū)估算的土壤可蝕性K值小,必須修訂。土壤性質(zhì)不同,修訂公式有所不同?？讈喥降萚21]采用短坡長小區(qū)進(jìn)行土壤可蝕性模擬時(shí),計(jì)算結(jié)果明顯偏小,必須進(jìn)行修訂,此結(jié)果與筆者研究結(jié)果一致。

2)微型小區(qū)不適用于易于產(chǎn)生細(xì)溝的黃土等土壤的可蝕性的估算。

3)土壤可蝕性為標(biāo)準(zhǔn)小區(qū)上單位降雨侵蝕力所引起的土壤流失量,反映的是某一特定土壤的長期平均狀況,多年平均值需要有足夠大且觀測(cè)歷時(shí)較長(一般不少于20年)數(shù)據(jù)庫資料,以消除降雨特性等因素年際變化的影響。筆者采用的天然降雨條件下標(biāo)準(zhǔn)小區(qū)6年和微型小區(qū)5年的實(shí)測(cè)K值數(shù)據(jù)每年都存在一定差異,說明用較長歷時(shí)資料來研究K值的必要性,筆者研究的觀測(cè)數(shù)據(jù)時(shí)段明顯不夠,還需要更長時(shí)間的觀測(cè)資料來進(jìn)行補(bǔ)充修正。

4)筆者研究標(biāo)準(zhǔn)小區(qū)估算5種土壤的K值明顯大于全國各侵蝕區(qū)野外布設(shè)的徑流小區(qū)觀測(cè)數(shù)據(jù)計(jì)算的5種土壤的K值。主要原因是筆者為了對(duì)比研究試驗(yàn)方法,不研究每種土壤的K值,將不同地區(qū)的耕層土壤放在一個(gè)地方進(jìn)行實(shí)驗(yàn),忽略不同地區(qū)降雨及成土等各種因素的影響。

[1] Hudson N.Soil Conservation[M].AMES:Iowa State U-niversity Press,1995:93

[2] 史學(xué)正,鄧西海.土壤可蝕性研究現(xiàn)狀及展望[J].中國水土保持,1993(5):25- 29

[3] Olson T C,Wischmeier W H.Soil erodibility evaluation for soils on the runoff and erosion stations[J].Soil Science,Society of American Proceedings,1963,27(5): 590- 592

[4] 張科利,彭文英,楊紅麗.中國土壤可蝕性值及其估算[J].土壤學(xué)報(bào),2007,44(1):7- 13

[5] 張文太,于東升,史學(xué)正,等.中國亞熱帶土壤可蝕性K值預(yù)測(cè)的不確定性研究[J].土壤學(xué)報(bào),2009,46(2): 185- 191

[6] United States Department of Agriculture.EP IC-erosion/ productivity impact calculator 1.Model Documentation, Technical Bulletin Number 1768[R].Washington,D C: USDA-ARS,1990

[7] 史學(xué)正,于東升,邢廷炎.用田間實(shí)測(cè)法研究我國亞熱帶土壤的可蝕性K值[J].土壤學(xué)報(bào),1997,34(4): 399- 405

[8] Wischmeier W H,Johnson C B,Cross B V.A soil erodibility nomograph for farmland and construction sites[J]. Soil and Water conservation,1971,26(5):189- 193

[9] Shirazi M A,Boersma L.A unifying quantitative analysis of soil texture[J].Soil Science Society of America Journal,1984,48:142- 147

[10]翟偉峰,許林書.東北典型黑土區(qū)土壤可蝕性K值研究[J].土壤通報(bào),2011,42(5):1209- 1213

[11]高德武.黑龍江省土壤流失方程中土壤可蝕性因子(K)的研究[J].國土與自然資源研究,1993(3):40- 43

[12]謝云,林小鵑,劉瑛娜,等.槽式擺噴頭下噴式人工模擬降雨機(jī)的雨強(qiáng)及其空間分布率定[J].水土保持通報(bào),2008,28(4):1- 6

[13]Wischmeier W H,Smith D D.Predicting rainfall erosion losses:a guide to conservation planning[R].Washington,D C:USDA,1978

[14]Liu B Y,Zhang K L,Xie Y.An empirical soil loss equation[C]∥Proceedings of 12th ISCO Conference.Vol.Ⅱ Process of soil Erosion and Its Environment Effect.Beijing:Tsinghua University Press,2002:21- 25

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

[16]張科利,蔡永明,劉寶元,等.土壤可蝕性動(dòng)態(tài)變化規(guī)律研究[J].地理學(xué)報(bào),2001,56(6):673- 681

[17]Fu S H,Liu B Y,Liu H P.The effect of slope on interrill erosion at short slopes[J].Catena,2011,84:29- 34

[18]Meyer L D,Harmon W C.Multiple-intensity rainfall simulator for erosion research on row sideslopes[J].Transaction of the ASAE,1979,22:100- 103

[19]Lal R.坡長與坡度對(duì)水土流失的影響[J].水土保持科技情報(bào),1991(4):26- 29

[20]李卓,吳普特,馮浩,等.不同粘粒含量土壤水分入滲能力模擬試驗(yàn)研究[J].干旱地區(qū)農(nóng)業(yè)研究,2009,27 (3):71- 77

[21]孔亞平,張科利,楊紅麗.土壤可蝕性模擬研究中的坡長選定問題[J].地理科學(xué),2005,25(3):374- 378

(責(zé)任編輯:郭雪芳)

Soil erodibility estimation using micro-plots

Liu Yingna1,Liu Baoyuan1,2,Zhang Keli1,2,Zheng Yuanzhi3
(1.State Key Laboratory of Earth Surface Processes and Resource Ecology,Beijing Normal University,100875,Beijing,China；2.School of Geography,Beijing Normal University,100875,Beijing,China；3.Rongfeng Holding Grop Co.,Ltd,100070,Beijing,China)

Soil erodibility reflects the vulnerability or susceptibility of soil to erosion and is frequently used in the evaluation of soil erosion and the prediction of soil loss.A feasible and reliable method for soil erodibility measurement and estimation is critically required in the research of soil erosion.The objectives of the current study were to estimate the erodibilities of five typical soils in China using microplots and to compare the results with those obtained with the traditional unit plot method.Each of the five soils,i.e.,purple soil,black soil,red soil,cinnamon soil and loess soil,was packed into a 20 m×2.1 m unit plot and a 1 m×1 m micro-plot,with continuously in a clean-tilled fallow condition and tillage performed upslope and downslope.Runoff and sediment yield were measured for the unit plots in 61 natural rainfall events during the years 2006—2011 and for the micro-plots in 47 events during 2007—2011.Annual soil erodibilities of each soil were calculated as the K values in the Chinese Soil Loss Equation(CSLE)for both unit plot and micro-plot.And,the calculated formula of K value was presented as K=A/(R·LS·B·E·T).The results showed that the mean annual K values generated withthe micro-plots were all lower than those obtained with the unit plots.However,both plot methods yielded generally the highest K values for the purple and black soils,followed by the red soil,and the lowest values for the cinnamon and loess soils.In addition,significant correlations of annual K value were detected between the two methods for all the soils except for the loess soil,because rills were prone to develop in the unit plot of loess soil but not in the micro-plot due to the short slope length.Significant linear regression equations between two K values,i.e.,obtained with the unit plot and the micro-plot methods,were generated for purple soil,black soil,red soil and cinnamon soil.These findings suggest that micro-plots are applicable in soil erobibility estimation for five typical soils in China except for the loess soil.

soil erodibility；simulated rain；micro-plot；unit plot

S157.1

A

1672-3007(2015)04-0103-06

2014- 11- 02

2015- 06- 25

項(xiàng)目名稱:中央高?；究蒲袠I(yè)務(wù)費(fèi)專項(xiàng)資金“東北典型小流域黑土埋藏年齡與成土年齡研究”(2013NT36)；中國博士后科學(xué)基金資助項(xiàng)目“黑龍江省典型小流域黑土埋藏年齡及成土速率研究”(2014M550646)；國家自然科學(xué)基金“地表過程模型與模擬創(chuàng)新研究群體科學(xué)基金”(41321001)

劉瑛娜(1982—),女,博士,工程師。主要研究方向:土壤侵蝕及其環(huán)境影響。E-mail:liuyingna@bnu.edu.cn

?通信作者簡介:劉寶元(1958—),男,博士,教授。主要研究方向:土壤侵蝕與水土保持。E-mail:baoyuan@bnu.edu.cn