肖 海，夏振堯，楊?lèi)偸妫河勒?，詹檳?/p>

（1.三峽大學(xué) 土木與建筑學(xué)院，湖北 宜昌443002；2.三峽地區(qū)地質(zhì)災(zāi)害與生態(tài)環(huán)境湖北省協(xié)同創(chuàng)新中心，湖北 宜昌443002；3.三峽大學(xué) 水利與環(huán)境學(xué)院，湖北 宜昌443002）

由于長(zhǎng)期不合理的土地利用方式，加之較頻繁的侵蝕性降雨，導(dǎo)致三峽庫(kù)區(qū)水土流失嚴(yán)重，不僅致使庫(kù)區(qū)土層日趨貧瘠，當(dāng)?shù)剞r(nóng)業(yè)生態(tài)環(huán)境逐漸惡化，而且水庫(kù)泥沙淤積給三峽工程和長(zhǎng)江中下游的安危造成極大威脅［1］。小流域是水土保持研究的基本單元，對(duì)其侵蝕產(chǎn)沙空間分布的了解是預(yù)測(cè)和防治流域水土流失的重點(diǎn)［2］。

稀土元素（REE）由于其性質(zhì)穩(wěn)定，能夠被土壤顆粒強(qiáng)烈吸附，植物富集有限，淋溶遷移不明顯，對(duì)環(huán)境無(wú)放射性危害，有較低的土壤背景值，并可人工同時(shí)施放多種稀土元素示蹤，使研究的目的性和精確度提高等特點(diǎn)［3］，因此已被應(yīng)用于小流域泥沙來(lái)源示蹤方面的研究。吳普特和劉普靈［4］最先對(duì)坡面—溝坡連續(xù)體模型REE施放方法進(jìn)行了試驗(yàn)，以期為黃土高原小流域泥沙來(lái)源觀(guān)測(cè)提供理論基礎(chǔ)。石輝等［5］通過(guò)室內(nèi)模擬試驗(yàn)，將REE示蹤法引入黃土高原小流域泥沙來(lái)源研究。其研究結(jié)果證明REE示蹤法可比較滿(mǎn)意地說(shuō)明小流域侵蝕產(chǎn)沙空間分布情況。Polyakov等［6－7］將稀土元素均勻播撒在小流域不同地貌部位，從而定量分析了小流域不同地貌部位土壤侵蝕及沉積的空間分布。綜上所述，利用REE示蹤技術(shù)研究小流域侵蝕產(chǎn)沙空間分布規(guī)律具有一定的優(yōu)勢(shì)，但是其研究成果多集中在黏粒和粉粒含量較高的黃土區(qū)及粉砂壤土地區(qū)。由于不同土壤顆粒粒徑對(duì)于REE元素吸附能力存在較大差異［3，8］，因此 REE示蹤技術(shù)在土壤砂粒含量較高且泥沙分選較強(qiáng)烈的三峽庫(kù)區(qū)小流域是否能夠應(yīng)用仍需進(jìn)一步研究。

本研究以三峽庫(kù)區(qū)主要土壤類(lèi)型紫色土為研究對(duì)象，通過(guò)建立小流域微縮模型，結(jié)合人工模擬降雨試驗(yàn)及REE示蹤技術(shù)，計(jì)算該小流域不同地貌部位侵蝕量及相對(duì)貢獻(xiàn)比例，并分析溝道系統(tǒng)與坡面系統(tǒng)相對(duì)貢獻(xiàn)比例。開(kāi)展本項(xiàng)研究，能夠?yàn)镽EE示蹤技術(shù)在三峽庫(kù)區(qū)野外小流域的應(yīng)用提供技術(shù)支撐，并為該地區(qū)水土流失預(yù)測(cè)模型的建立提供理論依據(jù)。

1 材料與方法

1.1 建立小流域模型

研究區(qū)位于三峽庫(kù)區(qū)湖北宜昌段，以秭歸縣水田壩鄉(xiāng)的王家橋小流域（31°12′—31°15′N(xiāo)，110°40′—110°43′E）為原型，通過(guò)實(shí)地勘察，按1∶100的比例尺將其微縮模型建于三峽大學(xué)人工模擬降雨試驗(yàn)場(chǎng)。小流域模型外圍邊界用紅磚及水泥砌成，模型中下層鋪設(shè)砂土以便透水。上層鋪設(shè)10cm厚的紫色土層，其鋪設(shè)過(guò)程如下：首先，使用逐步稀釋法將不同稀土元素氧化物（La2O3，CeO2，Nd2O3，Sm2O3，Eu2O3，Tb4O7，Ho2O3，Yb2O3共8種）分別均勻混合于其中；其次，根據(jù)流域地形圖，以坡面和坡溝分界處每50cm取一點(diǎn)作為控制點(diǎn)，并根據(jù)控制點(diǎn)在模型中描繪溝緣線(xiàn)；然后，確定小流域中不同地貌部位的分區(qū)；接著，以鋁皮隔開(kāi)各分區(qū)，內(nèi)部插上鐵絲控制釋放深度，將混合有稀土元素的紫色土按照1.3g／cm3的容重逐層地均勻鋪設(shè)于各分區(qū)之中，鋪設(shè)過(guò)程中適當(dāng)灑水；最后，鋪設(shè)完成后用小鐵鏟等工具修整小流域微地貌，灑水養(yǎng)護(hù)，用防雨布遮擋靜置三個(gè)月，以備試驗(yàn)。

建成后小流域模型面積約為20m2，其最高處約1.3m，最低處高約0.4m，位于小流域出口處。由一條主溝、兩條支溝組成。主溝長(zhǎng)約5.5m，上半主溝溝底坡度為10°～15°，溝壁坡度為20°～35°；下半主溝溝底坡度為0°～9°，溝壁坡度為15°～30°。支溝位于主溝兩側(cè)，左支溝長(zhǎng)約2.2m，溝底坡度小于5°，溝壁坡度在25°～35°；右支溝長(zhǎng)約1.9m，溝底坡度小于5°，溝壁坡度在15°～30°。各坡面坡度在10°以?xún)?nèi)，下半坡面坡度整體比上半坡面坡度要大。綜合考慮了坡、溝系統(tǒng)的分布以及稀土元素的種類(lèi)和用量后，將小流域劃分成14個(gè)小區(qū)域。分別用8個(gè)符號(hào)代表8種泥沙來(lái)源地：Ⅰ為L(zhǎng)a2O3區(qū)；Ⅱ?yàn)閅b2O3區(qū)；Ⅲ為CeO2區(qū)；Ⅳ為Sm2O3區(qū)；Ⅴ為Eu2O3區(qū)；Ⅵ為T(mén)b4O7區(qū)；Ⅶ為Nd2O3區(qū)；Ⅷ為Ho2O3區(qū)。

1.2 稀土元素釋放濃度

按以下公式計(jì)算試驗(yàn)小流域各區(qū)域稀土元素釋放濃度［9］，結(jié)果見(jiàn)表1。

式中：Cj——j區(qū)第j 種元素的施放濃度（mg／kg）；Bj——土壤中第j種元素的背景值濃度（mg／kg）；λj——第j區(qū)相對(duì)侵蝕量的最小期望值；K——考慮到其它因素的綜合保證系數(shù)。

表1 稀土元素釋放量計(jì)算結(jié)果

1.3 人工降雨過(guò)程及處理

人工降雨試驗(yàn)在三峽大學(xué)人工模擬降雨試驗(yàn)場(chǎng)開(kāi)展。60，90，120mm／h三次降雨時(shí)間依次為2012年10月27日、2012年11月1日以及2012年11月10日。降雨設(shè)施采用南京南林電子科技生產(chǎn)的NLJY-10-01型便攜式人工模擬降雨系統(tǒng)，降雨高度4m，由9組大中小噴頭均勻噴出雨滴，降雨均勻度大于0.86，降雨強(qiáng)度為60，90，120mm／h。利用均勻分布于小流域中的雨量筒，通過(guò)多次測(cè)定平均降雨強(qiáng)度對(duì)降雨強(qiáng)度率定。每次降雨開(kāi)始前，先對(duì)小流域進(jìn)行霧化降雨使土壤表層濕潤(rùn)，直至接近產(chǎn)生地表徑流。然后，將降雨強(qiáng)度調(diào)至設(shè)定值，在降雨開(kāi)始產(chǎn)流后開(kāi)始計(jì)時(shí)，以2min為時(shí)間間隔接取全部泥沙樣品以待測(cè)試，同時(shí)照相，并對(duì)試驗(yàn)過(guò)程進(jìn)行記錄。降雨30 min結(jié)束后，覆蓋防雨布以備下次試驗(yàn)使用。泥沙樣品接收后靜置至泥沙沉淀，倒掉上部清液，并將下部泥沙烘干稱(chēng)重保存于自封袋中。

1.4 稀土元素含量分析

本研究采用ICP-MS方法測(cè)定土壤及泥沙樣品中的REE含量，其樣品制備方法如下：取0.025g樣品放入聚四氟乙烯（PTFE）材質(zhì)的消解罐之中，然后緩緩加入0.5ml優(yōu)級(jí)純濃硝酸（HNO3）和1ml優(yōu)級(jí)純濃氫氟酸（HF），蓋好罐子后放入鋼套之中置于烘箱中用185℃加熱24小時(shí)。待冷卻后開(kāi)蓋，取出PTFE內(nèi)罐放入電熱板上用130℃的溫度蒸發(fā)至近干，然后加入1ml優(yōu)級(jí)純HNO3繼續(xù)蒸發(fā)至近干，再重復(fù)加酸、蒸干一次。之后往PTFE內(nèi)罐中加入5ml的1∶1的HNO3，封蓋后在放入鋼套中置于130℃的烘箱中復(fù)溶。三個(gè)小時(shí)后取出，待PTFE內(nèi)罐冷卻后倒出液體至50ml的容量瓶中定容，搖晃均勻后取10 ml溶液放置于10ml離心管（EP管）待測(cè)。將EP管接入電感耦合等離子體質(zhì)譜儀（IPC-MS，賽默飛世爾科技公司X Series 2）測(cè)量8種REE元素的含量。每批樣品測(cè)定之前，進(jìn)行標(biāo)準(zhǔn)曲線(xiàn)測(cè)定。為了保證分析精度，測(cè)量每個(gè)樣品時(shí)，加入濃度為10μg／L的Rh和Re元素作為內(nèi)標(biāo)一起測(cè)量。

1.5 數(shù)據(jù)處理

小流域中各分區(qū)土壤的侵蝕量、貢獻(xiàn)率及計(jì)算誤差用以下公式計(jì)算［9］：

式中：wj——j示蹤區(qū)的侵蝕量（kg）；Rj——侵蝕泥沙中第j種稀土元素的實(shí)測(cè)濃度（mg／kg）；Bj——土壤中第j種稀土元素的背景值（mg／kg）；W——所接泥沙實(shí)際重量（kg）；γj——j區(qū)的相對(duì)侵蝕量（kg）；σ——侵蝕量的計(jì)算誤差。

2 結(jié)果與分析

2.1 徑流量及土壤侵蝕量變化

表2表示三次降雨徑流量及侵蝕量對(duì)比，可以看出，徑流量隨著降雨強(qiáng)度的增大而增加，但是增加幅度相差較大。雖然三次降雨強(qiáng)度均相差30mm／h，但與60mm／h相比，90mm／h徑流量增加325.106kg，120mm／h相對(duì)于90mm／h只增加261.616kg，增幅減少了19.53%。分析認(rèn)為，第一場(chǎng)60mm／h降雨前雖然進(jìn)行了預(yù)濕潤(rùn)，但其下層土壤并未飽和，降雨過(guò)程中入滲較大，從而使得徑流量減小。而隨著第一次降雨的進(jìn)行，土壤基本趨于飽和，后兩次入滲相對(duì)較小，徑流量增加，從而導(dǎo)致在三場(chǎng)降雨過(guò)程中徑流增幅減小。

與徑流量變化規(guī)律相同，三場(chǎng)降雨侵蝕量隨著降雨強(qiáng)度的增大而增加，且增加幅度相差較大：降雨強(qiáng)度為90mm／h侵蝕量比60mm／h增加20.852kg，120mm／h侵蝕量比90mm／h增加1.79kg，增幅較前兩場(chǎng)減少91.42%，這個(gè)增幅遠(yuǎn)小于同期徑流量的增幅。60mm／h降雨結(jié)束后坡面僅出現(xiàn)了細(xì)溝雛形，細(xì)溝并未完全發(fā)育，因此侵蝕量較少。而在隨后的90mm／h降雨強(qiáng)度下，侵蝕量陡增，主要是由于在前場(chǎng)降雨的基礎(chǔ)上，一方面坡面細(xì)溝繼續(xù)發(fā)育，而且較大的降雨強(qiáng)度導(dǎo)致了較快的發(fā)育速度［10］；另一方面較大的徑流量使得溝道侵蝕迅速增加［11］。在120 mm／h降雨強(qiáng)度下，由于坡面細(xì)溝以及溝道發(fā)育已經(jīng)較為穩(wěn)定，降雨強(qiáng)度及徑流量的增加并不能導(dǎo)致侵蝕量的顯著增加，因此在整個(gè)降雨過(guò)程中侵蝕量增幅不大。

2.2 侵蝕量和相對(duì)貢獻(xiàn)率分析

由于本研究中小流域侵蝕是在連續(xù)三場(chǎng)降雨的情況下發(fā)生，每次降雨之后并不恢復(fù)原始地貌，下次降雨是在上次降雨的基礎(chǔ)上進(jìn)行。因此，隨著流域內(nèi)溝道的發(fā)育，各部位的侵蝕量及相對(duì)貢獻(xiàn)率均發(fā)生了變化（圖1）。由圖1可知：整體上三場(chǎng)降雨中侵蝕量的空間分布規(guī)律一致，均是主溝溝道Ⅴ、Ⅵ區(qū)最大，靠近主溝的Ⅱ、Ⅳ及Ⅶ次之，位于小流域邊緣的Ⅰ、Ⅲ及Ⅶ三區(qū)侵蝕量相對(duì)較小。分析認(rèn)為，這主要與各部位的坡度及位置有關(guān)，Ⅰ、Ⅲ及Ⅶ三區(qū)位于小流域周邊，無(wú)其他區(qū)域匯集徑流流于這三個(gè)區(qū)域，而且坡度較小，因此發(fā)生的侵蝕不多；而靠近主溝的Ⅱ、Ⅳ及Ⅷ三區(qū)雖然坡度也不大，但是由Ⅰ、Ⅲ及Ⅶ三區(qū)產(chǎn)生的徑流匯集沖刷，因此其侵蝕量相對(duì)于Ⅰ、Ⅲ及Ⅶ三區(qū)均略有增大；而主溝溝道Ⅴ、Ⅵ區(qū)由于坡度大而且有各區(qū)徑流匯集于其中對(duì)溝底及溝壁進(jìn)行沖刷產(chǎn)生侵蝕，因而這兩區(qū)侵蝕量遠(yuǎn)大于其他各區(qū)。這個(gè)規(guī)律表明小流域內(nèi)侵蝕以主溝為為中心，向四周逐漸減弱。然而各區(qū)侵蝕量的變化均隨雨強(qiáng)增大呈現(xiàn)不一樣的規(guī)律：Ⅰ、Ⅱ、Ⅲ、Ⅶ四區(qū)侵蝕量隨雨強(qiáng)增加而增大，而Ⅳ、Ⅴ、Ⅵ、Ⅷ四區(qū)侵蝕量卻隨雨強(qiáng)增加呈現(xiàn)先增加后減小的趨勢(shì)，同時(shí)在降雨強(qiáng)度從60mm／h增加到90 mm／h時(shí)，各區(qū)侵蝕量增幅較大，這是因?yàn)橛陱?qiáng)增大，雨滴能量及徑流量增加導(dǎo)致坡面和溝道侵蝕量均大幅度增加［10－11］。而雨強(qiáng)從90mm／h到120mm／h時(shí)變化幅度較小，分析認(rèn)為，經(jīng)過(guò)前兩場(chǎng)降雨的侵蝕，小流域內(nèi)坡面細(xì)溝及溝道發(fā)育較為成熟，從而使得坡面各區(qū)侵蝕量增幅不大，且溝道系統(tǒng)各區(qū)侵蝕量還略有減少。

從侵蝕量的相對(duì)貢獻(xiàn)率來(lái)看，各區(qū)的貢獻(xiàn)率規(guī)律大致相同，均是主溝溝道Ⅴ、Ⅵ區(qū)最大，靠近主溝的Ⅱ、Ⅳ及Ⅶ次之，位于小流域邊緣的Ⅰ、Ⅲ及Ⅶ三區(qū)侵蝕量相對(duì)較小，這與各區(qū)侵蝕量規(guī)律變化一致。三次降雨強(qiáng)度下，各區(qū)相對(duì)侵蝕量變化較為復(fù)雜，坡面系統(tǒng)Ⅰ、Ⅱ及Ⅲ三區(qū)以及支溝Ⅶ區(qū)隨雨強(qiáng)增加而先大幅度增大而后小幅度增大，坡面系統(tǒng)Ⅳ區(qū)隨雨強(qiáng)增加呈現(xiàn)先大幅度增加后小幅度減少的趨勢(shì)，而主溝系統(tǒng)Ⅴ、Ⅵ兩區(qū)以及支溝Ⅷ區(qū)則呈現(xiàn)隨雨強(qiáng)增加而減小的規(guī)律。分析認(rèn)為降雨強(qiáng)度從60mm／h增加到90 mm／h時(shí)，可以觀(guān)察到60mm／h降雨強(qiáng)度下坡面所產(chǎn)生的細(xì)溝雛形快速發(fā)育，因此坡面及溝道侵蝕量大幅度增加，但坡面系統(tǒng)侵蝕量增加幅度略大于溝道系統(tǒng)，而在120mm／h降雨強(qiáng)度下，由于坡面細(xì)溝以及溝道發(fā)育已經(jīng)較為穩(wěn)定，雨強(qiáng)的增大并不能夠大幅度增加坡面侵蝕量，相反溝道系統(tǒng)中侵蝕量還略有減少，因而坡面系統(tǒng)各區(qū)相對(duì)貢獻(xiàn)率整體略有增大，而溝道系統(tǒng)整體略有減小。

圖1 不同降雨強(qiáng)度下小流域各分區(qū)土壤侵蝕量及貢獻(xiàn)率變化

2.3 溝道與坡面侵蝕分析

表3表明，溝道系統(tǒng)侵蝕泥沙貢獻(xiàn)率遠(yuǎn)大于坡面系統(tǒng)，尤其是60mm／h降雨強(qiáng)度中溝道系統(tǒng)侵蝕貢獻(xiàn)率達(dá)到79.241%，遠(yuǎn)大于坡面系統(tǒng)侵蝕貢獻(xiàn)率的20.759%；90mm／h降雨強(qiáng)度下，溝道系統(tǒng)侵蝕貢獻(xiàn)率降低到65.382%，坡面系統(tǒng)則上升到34.678%；120mm／h降雨強(qiáng)度下，溝道系統(tǒng)侵蝕減小至60.881%，而坡面系統(tǒng)則增至39.119%。整體上來(lái)看，一方面溝道侵蝕貢獻(xiàn)率始終占主導(dǎo)地位，這可能是由于，坡面系統(tǒng)坡度均在10°以下，而溝道系統(tǒng)中溝壁及溝底坡度較大，尤其是溝壁坡度均在20°以上。對(duì)于其他條件相同的坡面來(lái)說(shuō)，坡度小于某一臨界坡度時(shí)，坡度越大，土壤侵蝕量越大［12－13］。另外，由于坡面徑流匯集到溝道中，其徑流流量和流速遠(yuǎn)大于坡面，因而溝道系統(tǒng)的侵蝕量大于坡面系統(tǒng)［14］。另一方面，溝道系統(tǒng)的泥沙貢獻(xiàn)率逐漸減小，坡面系統(tǒng)泥沙貢獻(xiàn)率則逐漸增大。這主要是由于，第一場(chǎng)降雨結(jié)束后，小流域表面土壤變得疏松，并初步形成了匯水區(qū)，留下坡面細(xì)溝的雛形。第二場(chǎng)及第三場(chǎng)降雨在此基礎(chǔ)上，細(xì)溝持續(xù)發(fā)育，已經(jīng)在坡面上匯集成徑流對(duì)坡面進(jìn)行沖刷。同時(shí)隨著雨強(qiáng)的增大，雨滴能量增大，坡面徑流增加，也會(huì)加快細(xì)溝的發(fā)育，從而導(dǎo)致坡面侵蝕量急劇增加［15］。因而溝道系統(tǒng)侵蝕量和總侵蝕量雖然均有所增加，但由于坡面侵蝕增加更快，從而導(dǎo)致坡面系統(tǒng)泥沙貢獻(xiàn)率逐漸增大。

從表3中還可以看出，在溝道系統(tǒng)中，主溝與支溝的侵蝕泥沙貢獻(xiàn)率同樣存在差異。支溝侵蝕泥沙占總侵蝕量的相對(duì)貢獻(xiàn)率變化略有減小但較為穩(wěn)定，基本維持在17%左右。而主溝相對(duì)總侵蝕泥沙的貢獻(xiàn)率降幅較大，三次降雨強(qiáng)度下相對(duì)貢獻(xiàn)率依此從61.229% 降為 48.520%，44.523%，總 降幅達(dá)到16.706%，這幾乎與溝道系統(tǒng)相對(duì)貢獻(xiàn)率降幅18.333%相當(dāng)，且主溝與溝道系統(tǒng)侵蝕泥沙貢獻(xiàn)率的每次變化幅度幾乎一致，由此表明主溝侵蝕是溝道系統(tǒng)中最主要的部分，能夠反映溝道系統(tǒng)的侵蝕規(guī)律。

在坡面系統(tǒng)之中，坡面上部與下部的侵蝕泥沙貢獻(xiàn)率亦有所差別。下坡面侵蝕泥沙貢獻(xiàn)率比上坡面要大，兩者貢獻(xiàn)率均隨降雨強(qiáng)度增加而增加，當(dāng)兩者增幅不一致，上坡面從第一次降雨到第三次相對(duì)貢獻(xiàn)率增加9.347%，尤其是第二、三場(chǎng)之間僅增加1.491%，而小坡面從第一次降雨到第三次相對(duì)貢獻(xiàn)率增加9.012%，這與下坡面增幅差不多，但上坡面第二、三場(chǎng)之間相對(duì)貢獻(xiàn)率增幅為2.956%，比下坡面相對(duì)貢獻(xiàn)率增幅近2倍。同時(shí)二者相對(duì)貢獻(xiàn)率之間的差值隨著降雨強(qiáng)度增加呈現(xiàn)先增加后減小的規(guī)律。分析認(rèn)為：這主要是由于第一次降雨時(shí)，上下坡面均以濺蝕為主，下坡面總體坡度大于上坡面，且在后期在下坡面形成細(xì)溝的雛形，因此呈現(xiàn)下坡面侵蝕貢獻(xiàn)率大于上坡面，而第二場(chǎng)降雨時(shí)，由于降雨強(qiáng)度的增加，上下坡面侵蝕均會(huì)增加，而下坡面除了降雨強(qiáng)度增加引起的侵蝕增加外，下坡面原有的細(xì)溝雛形在大雨強(qiáng)作用下快速發(fā)育成細(xì)溝，因此呈現(xiàn)下坡面侵蝕貢獻(xiàn)率快速增加上下坡面侵蝕貢獻(xiàn)率差值增大的現(xiàn)象，而第三場(chǎng)降雨時(shí)，隨著降雨強(qiáng)度及降雨歷時(shí)的增加，細(xì)溝逐漸向上半坡面發(fā)育，發(fā)生溯源侵蝕，從而導(dǎo)致坡面上部侵蝕量迅速增大［6］，從而出現(xiàn)上下坡面侵蝕貢獻(xiàn)率差值減小的情況。

表3 溝道系統(tǒng)與坡面系統(tǒng)相對(duì)貢獻(xiàn)率對(duì)比

2.4 試驗(yàn)精確度分析

利用公式（4）對(duì)試驗(yàn)誤差進(jìn)行計(jì)算，計(jì)算結(jié)果見(jiàn)表4。由表4可知，三場(chǎng)降雨的實(shí)驗(yàn)誤差分別為33.172%，25.136%和18.014%，說(shuō)明利用 REE 示蹤技術(shù)計(jì)算土壤侵蝕量具有一定的精度和準(zhǔn)確性，并且隨著降雨的進(jìn)行，試驗(yàn)誤差逐漸減小，這與黑土地區(qū)試驗(yàn)結(jié)果相似［16］。但這一結(jié)果與黃土地區(qū)的研究成果相比［17］，其誤差值相對(duì)較大。這個(gè)可能是由于紫色土土壤細(xì)顆粒含量低于黃土［18］，而細(xì)顆粒土壤與粗顆粒土壤對(duì)于稀土元素吸附能力存在較大差異造成［2］。

表4 實(shí)驗(yàn)誤差分析

3 結(jié)論與建議

（1）三場(chǎng)降雨徑流量及侵蝕量均隨雨強(qiáng)增加，但徑流量的增加幅度遠(yuǎn)大于侵蝕量。

（2）本研究通過(guò)模擬人工降雨試驗(yàn)，結(jié)合REE示蹤技術(shù)，對(duì)三峽庫(kù)區(qū)模擬小流域侵蝕產(chǎn)沙空間分布進(jìn)行了研究。結(jié)果表明，三場(chǎng)降雨中各區(qū)侵蝕量規(guī)律一致，均是隨雨強(qiáng)增大而增加，且均是主溝溝道Ⅴ、Ⅵ區(qū)最大，靠近主溝的Ⅱ、Ⅳ及Ⅶ次之，位于小流域邊緣的Ⅰ、Ⅲ及Ⅶ三區(qū)侵蝕量相對(duì)較小。

（3）溝道系統(tǒng)土壤侵蝕貢獻(xiàn)率大于坡面系統(tǒng)，并呈現(xiàn)溝道系統(tǒng)的泥沙貢獻(xiàn)率逐漸減小，坡面系統(tǒng)泥沙貢獻(xiàn)率則逐漸增大的趨勢(shì)。

（4）三場(chǎng)降雨的侵蝕示蹤結(jié)果具有一定準(zhǔn)確性，但相比黃土區(qū)研究結(jié)果其誤差值略高。主要是侵蝕土壤顆粒大小與原始土壤不同以及不同顆粒的稀土元素吸附濃度不一樣引起的，今后可以通過(guò)考慮這兩個(gè)因素對(duì)實(shí)驗(yàn)數(shù)據(jù)進(jìn)行修正，以減小實(shí)驗(yàn)誤差。

（5）利用稀土元素示蹤三峽庫(kù)區(qū)小流域模型土壤侵蝕，可以比較準(zhǔn)確地描述不同地形部位的侵蝕量、相對(duì)貢獻(xiàn)率及各自隨時(shí)間變化的過(guò)程，從而揭示小流域土壤侵蝕變化規(guī)律。由此為進(jìn)一步開(kāi)展三峽庫(kù)區(qū)小流域野外土壤侵蝕規(guī)律的研究提供一種可靠的方法。

［1］ 熊平生，謝世友，謝金寧.初探三峽水庫(kù)濕地面臨的問(wèn)題及其對(duì)策［J］.國(guó)土與自然資源研究，2004（4）：60-62.

［2］ Kimoto A，Nearing M A，Shipitalo M J，et al.Multiyear tracking of sediment sources in a small agricultural watershed using rare earth elements［J］.Earth Surface Processes and Landforms，2006，31（14）：1763-1774.

［3］ Zhang X C，F(xiàn)riedrich J M，Nearing M A，et al.Potential use of rare earth oxides as tracers for soil erosion and aggregation studies［J］.Soil Science Society of America Journal，2001，65（5）：1508-1515.

［4］ 吳普特，劉普靈.溝坡侵蝕REE示蹤法試驗(yàn)研究初探［J］.水土保持研究，1997，4（2）：69-74.

［5］ Shi H，Tian J，Liu P，et al.A study on sediment sources in a small watershed by using REE tracer method［J］.Science in China Series E：Technological Sciences，1997，40（1）：12-20.

［6］ Polyakov V O，Kimoto A，Nearing M A，et al.Tracing sediment movement on a semiarid watershed using rare earth elements［J］.Soil Science Society of America Journal，2009，73（5）：1559-1565.

［7］ Polyakov V O，Nearing M A，Shipitalo M J.Tracking sediment redistribution in a small watershed：implications for agro-landscape evolution［J］.Earth surface processes and landforms，2004，29（10）：1275-1291.

［8］ Kimoto A，Nearing M A，Zhang X C，et al.Applicabil

ity of rare earth element oxides as a sediment tracer for

coarse-textured soils［J］.Catena，2006，65（3）：214-221.［9］ 李雅琦，吳普特，劉普靈，等.REE示蹤法研究土壤侵蝕

的室內(nèi)模擬試驗(yàn)［J］.水土保持研究，1997，4（2）：26-33.［10］ 劉剛，楊明義，劉普靈，等.7Be示蹤坡耕地次降雨細(xì)溝

與細(xì)溝間侵蝕［J］.農(nóng)業(yè)工程學(xué)報(bào)，2009，25（5）：47-53.［11］ 丁文峰，李勉，張平倉(cāng)，等.坡溝系統(tǒng)侵蝕產(chǎn)沙特征模擬試驗(yàn)研究［J］.農(nóng)業(yè)工程學(xué)報(bào)，2006，22（3）：10-14.

［12］ 曹銀真.黃土地區(qū)梁峁坡的坡地特征與土壤侵蝕［J］.地理研究，1983，2（3）：19-29.

［13］ 張會(huì)茹，鄭粉莉，耿曉東.地面坡度對(duì)紅壤坡面土壤侵蝕過(guò)程的影響研究［J］.水土保持研究，2009，16（4）：52-54，59.

［14］ Morgan R P C.Soil erosion and conservation（3rdedition）［M］.Wiley-Blackwell，2005.

［15］ 耿曉東，鄭粉莉，劉力.降雨強(qiáng)度和坡度雙因子對(duì)紫色土坡面侵蝕產(chǎn)沙的影響［J］.泥沙研究，2010（6）：48-53.

［16］ 王寧，楊春雨，張剛，等.黑土區(qū)土壤侵蝕的REE示蹤法研究［J］.地理科學(xué)，2008，28（4）：565-570.

［17］ 劉普靈，田均良，周佩華，等.土壤侵蝕稀土元素示蹤法操作技術(shù)研究［J］.水土保持研究，1997，4（2）：10-16.

［18］ 耿曉東.主要水蝕區(qū)坡面土壤侵蝕過(guò)程與機(jī)理對(duì)比研究［D］.北京：中國(guó)科學(xué)院研究生院，2011.