聶曉剛，梁 博，楊東升，喻 武，萬(wàn) 丹，方江平

(西藏農(nóng)牧學(xué)院 a資源與環(huán)境學(xué)院,b高寒水土保持研究中心,c高原生態(tài)研究所，西藏 林芝 860000)

土壤抗蝕性是結(jié)合土壤物理、化學(xué)性質(zhì)和生物酶等多重因素于一體，用以體現(xiàn)土壤抵抗水的分散和懸浮能力的綜合性因素指標(biāo)體系，是評(píng)定土壤抗蝕力的重要參數(shù)之一[1]。土壤抗蝕性的大小除了與土壤理化性質(zhì)等內(nèi)在因素有關(guān)外，還受植被類型等外部因素的影響[2]。氣候、地形、植被等的不同，形成了各種類型的土壤及環(huán)境條件，導(dǎo)致土壤性質(zhì)存在明顯的差異，稱之為土壤空間差異，而研究不同尺度上土壤的空間差異，不但對(duì)了解土壤的形成過(guò)程、結(jié)構(gòu)和功能具有重要的理論意義[3]，而且對(duì)了解土壤抗蝕性的空間格局等也具有重要的參考價(jià)值。海拔作為高山地區(qū)主要的自然環(huán)境特征指標(biāo)，是影響區(qū)域光、熱、水、氣的環(huán)境因子之一，并直接影響著生境的氣候生態(tài)學(xué)特征，通過(guò)氣候環(huán)境的改變使土壤物理、化學(xué)和生物等性質(zhì)發(fā)生變化，表現(xiàn)為不同的垂直地帶土壤的有機(jī)質(zhì)礦化程度、土壤機(jī)械組成以及微生物數(shù)量[4]等的不同，從而進(jìn)一步影響植物群落結(jié)構(gòu)和類型的演變[5-6]，最終引起生態(tài)系統(tǒng)功能的改變[7]。

青藏高原是世界上中低緯度面積最大的高海拔凍土區(qū)[8-9]，其能量和水分循環(huán)對(duì)亞洲季風(fēng)系統(tǒng)的形成和演化具有十分重要的作用[10]。近年來(lái)，隨全球氣候變化和人類活動(dòng)的影響，青藏高原土壤物理過(guò)程、水文循環(huán)、熱量循環(huán)等過(guò)程變化劇烈[11]，生態(tài)環(huán)境破壞嚴(yán)重，如何加強(qiáng)青藏高原國(guó)家生態(tài)安全屏障建設(shè)，減少水土流失，已成為該地區(qū)重要的生態(tài)環(huán)境熱點(diǎn)問(wèn)題之一[12]。色季拉山作為青藏高原東南部典型的山脈體系，其植被、土壤類型沿海拔高度的升高有明顯的分布規(guī)律，是研究高寒區(qū)土壤物質(zhì)分布特征與植被、土壤、氣候響應(yīng)的的天然實(shí)驗(yàn)室[13]。目前關(guān)于土壤抗蝕性的研究主要集中在中低海拔地區(qū)，而針對(duì)高海拔寒區(qū)土壤抗蝕性的研究相對(duì)較少。為此，本試驗(yàn)以色季拉山西坡3 200～4 600 m海拔土壤為研究對(duì)象，借鑒國(guó)內(nèi)外常用的土壤抗蝕性評(píng)價(jià)指標(biāo)，分析不同海拔土壤抗蝕性的空間差異特征，以期為以色季拉山為代表的藏東南高原山地水土流失的預(yù)防治理提供理論依據(jù)。

1 材料與方法

1.1 研究區(qū)概況

色季拉山位于青藏高原東南部的林芝地區(qū)(94°28′～94°51′E，29°21′～29°50′N)，海拔2 200～5 300 m[13]，年平均日照時(shí)數(shù)1 150 h[14]，年均氣溫-0.7 ℃，極端最低溫-31.6 ℃(1991年1月)，極端最高溫24 ℃(2000年7月)；整個(gè)山脈呈東北-西南走向，地勢(shì)西高東低。年降水總量600～1 000 mm，且約75%的降水量集中在5-10月份[15]。隨著海拔升高，色季拉山垂直氣候帶明顯，從山腳到山頂分別為山地溫帶-亞高山寒溫帶-高山寒帶。在氣候帶影響下，山體植被生長(zhǎng)茂盛，類型豐富多樣且垂直譜帶明顯，從山頂?shù)缴侥_分別為山頂凍原稀疏植被-亞高山寒帶灌叢、草甸-亞高山寒溫帶暗針葉林-山地溫帶暗針葉林-山地暖溫帶針闊混交林、松林。山體土壤類型多種多樣，垂直分布譜帶明顯而完整[16]。

1.2 土壤樣品的采集與處理

2016年9月中旬，沿318國(guó)道采集色季拉山西坡海拔3 200～4 600 m土壤土樣。按照典型選樣方法，每2 00 m一個(gè)梯度，沿不同海拔梯度分別選取樣點(diǎn)，每個(gè)樣點(diǎn)按照五點(diǎn)法設(shè)置，分別采集0～10，10～20，20～30 cm土層環(huán)刀土和原狀土，帶回實(shí)驗(yàn)室用于土壤理化性質(zhì)測(cè)定。各采樣點(diǎn)周邊基本特征詳見(jiàn)表1 。

1.3 土壤抗蝕性指標(biāo)體系

為揭示不同空間土壤抗蝕性的差異及土壤各因子對(duì)抗蝕性的影響，本研究采用指標(biāo)體系法，從土壤的內(nèi)涵出發(fā)，堅(jiān)持指標(biāo)靈敏性、獨(dú)立性和協(xié)調(diào)性準(zhǔn)則[19-20]，對(duì)色季拉山土壤抗蝕性進(jìn)行定量評(píng)價(jià)。本研究共選擇3類11個(gè)土壤理化指標(biāo)作為土壤抗蝕性評(píng)價(jià)指標(biāo)。

1)土壤團(tuán)粒類指標(biāo)。包括粒徑>0.25 mm水穩(wěn)性團(tuán)聚體含量(X1)、粒徑>0.5 mm水穩(wěn)性團(tuán)聚體含量(X2)、粒徑>0.25 mm非水穩(wěn)性團(tuán)聚體含量(X3)、土壤水穩(wěn)性指數(shù)(X4)、團(tuán)聚體平均質(zhì)量直徑(X5)、團(tuán)聚體破壞率(X6)。

2)土壤基本物理類指標(biāo)。包括非毛管孔隙度(X7)、毛管孔隙度(X8)、土壤體積質(zhì)量(X9)、吸濕水含量(X10)。

3)有機(jī)膠體類指標(biāo)。有機(jī)質(zhì)含量(X11)。

表1 西藏色季拉山西坡土壤采樣點(diǎn)的基本特征Table 1 Basic characteristics of sampling points in west slope of the Sejila Mountain

注：研究區(qū)根據(jù)海拔可分為山地溫帶(3 200～3 600 m)、亞高山寒溫帶(3 800～4 200 m)、高山寒帶(4 400～4 600 m)[16-18]。

Note:According to elevation, the Sejila Mountain can be divided into mountain temperate zone (3 200-3 600 m), subalpine cold temperate zone (3 800-4 200 m) and alpine frigid zone (4 400-4 600 m)[16-18].

1.4 指標(biāo)測(cè)定及計(jì)算方法

采用文獻(xiàn)[21-22]的方法測(cè)定粒徑>0.25 mm水穩(wěn)性團(tuán)聚體含量(X1)、粒徑>0.5 mm水穩(wěn)性團(tuán)聚體含量(X2)、粒徑>0.25 mm非水穩(wěn)性團(tuán)聚體含量(X3)。采用環(huán)刀法[23]測(cè)定非毛管孔隙度(X7)、毛管孔隙度(X8)、土壤體積質(zhì)量(X9)、吸濕水含量(X10)。采用文獻(xiàn)[24]的方法測(cè)定有機(jī)質(zhì)含量。

土壤水穩(wěn)性指數(shù)(X4)采用靜水崩解法測(cè)定，其計(jì)算公式[25]為：

X4=(∑PiKi+Pj)/A。

式中：Pi為第i分鐘內(nèi)分散的土壤顆粒數(shù);Ki為第i分鐘的校正系數(shù);Pj為10 min內(nèi)未分散的土壤顆粒數(shù)；A為供試土壤顆粒總數(shù)；i為時(shí)間，i=1,2,…,10 min。

團(tuán)聚體平均質(zhì)量直徑(X5)采用邱麗萍等[26]方法的測(cè)定，其計(jì)算公式為：

團(tuán)聚體破壞率(X6)采用土壤團(tuán)粒分析法測(cè)定，其計(jì)算公式[27]為：

X6=(粒徑>0.25 mm團(tuán)聚體干篩分析值-粒徑>0.25 mm團(tuán)聚體濕篩分析值)/粒徑>0.25 mm團(tuán)聚體干篩分析值×100%。

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

以不同土層各抗蝕性指標(biāo)的平均值和變異系數(shù)描述色季拉山不同海拔區(qū)域土壤抗蝕性的空間變異性。利用SPSS 17.0、Origin 9.0軟件和單因素方差(ANOVA)分析進(jìn)行圖形制作和數(shù)據(jù)處理，利用主成分分析法選取抗蝕性最優(yōu)指標(biāo)，最后通過(guò)構(gòu)建評(píng)價(jià)模型分析比較色季拉山不同海拔區(qū)域土壤抗蝕性的變化規(guī)律。

2 結(jié)果與分析

2.1 色季拉山不同海拔區(qū)域土壤抗蝕性指標(biāo)的空間分異特征

2.1.1 海拔3 200～3 600 m山地溫帶土壤 由表2可以看出，在6個(gè)土壤團(tuán)粒類指標(biāo)中，除團(tuán)聚體破壞率(X6)外，其余指標(biāo)均隨著海拔上升總體表現(xiàn)出逐漸增大的趨勢(shì)，其中在海拔為3 200～34 00 m時(shí)增幅最大；團(tuán)聚體破壞率(X6)表現(xiàn)為3 200 m海拔最大，3 600 m海拔次之，3 400 m海拔最小，這可能與研究區(qū)人為干擾強(qiáng)度相關(guān)。土壤基本物理指標(biāo)隨海拔變化表現(xiàn)復(fù)雜，其中非毛管孔隙度(X7)和毛管孔隙度(X8)隨海拔的增加呈先增后減的趨勢(shì)；土壤體積質(zhì)量(X9)則隨海拔增加而逐漸減小，吸濕水含量(X10)隨海拔升高逐漸增加。有機(jī)膠體類指標(biāo)有機(jī)質(zhì)含量(X11)則隨海拔的升高而逐漸增大。

表2 西藏色季拉山海拔3 200～3 600 m山地溫帶土壤的抗蝕性Table 2 Soil-erodibility indexes at 3 200-3 600 m altitude of the Sejila Mountain

2.1.2 海拔3 800～4 200 m亞高山寒溫帶土壤 由表3可知，當(dāng)海拔為3 800～4 200 m時(shí)，在研究區(qū)土壤團(tuán)粒類指標(biāo)中，除團(tuán)聚體破壞率(X6)隨海拔上升呈增加趨勢(shì)之外，其余指標(biāo)總體均呈降低趨勢(shì)。在海拔為3 800～4 200 m時(shí)，土壤基本物理類指標(biāo)變化較為復(fù)雜，隨海拔高度增加，非毛管孔隙度(X7)先降后升，毛管孔隙度(X8)先增后降，二者變化規(guī)律相反，這主要與土壤表層植被根系的穿插、切割作用有關(guān)；土壤體積質(zhì)量(X9)逐漸增加， 吸濕水含量(X10)逐漸減小。對(duì)有機(jī)膠體類指標(biāo)而言，有機(jī)質(zhì)含量則隨海拔上升總體呈減小趨勢(shì)。

表3 西藏色季拉山海拔3 800～4 200 m亞高山寒溫帶土壤的抗蝕性Table 3 Soil-erodibility index at 3 800-4 200 m altitude of the Sejila Mountain

2.1.3 海拔4 400～4 600 m高山寒帶土壤 由表4可知，當(dāng)海拔為4 400～4 600 m時(shí)，在土壤團(tuán)粒類指標(biāo)中，隨海拔上升，除粒徑>0.25 mm非水穩(wěn)性土壤團(tuán)聚體含量(X3)、土壤水穩(wěn)性指數(shù)(X4)和土壤團(tuán)聚體平均質(zhì)量直徑(X5)呈降低趨勢(shì)外，其余指標(biāo)均呈增加趨勢(shì)，這主要與凍融作用有關(guān)。土壤基本物理類指標(biāo)中除土壤體積質(zhì)量(X9)隨海拔上升呈增加趨勢(shì)外，其余指標(biāo)均呈減小趨勢(shì)。有機(jī)質(zhì)含量隨海拔的升高呈遞減趨勢(shì)。

表4 西藏色季拉山海拔4 400～4 600 m高山寒帶土壤的抗蝕性Table 4 Soil-erodibility index at 4 400-4 600 m altitude of the Sejila Mountain

2.2 色季拉山不同土層土壤抗蝕性指標(biāo)的空間分異特征

不同土層由于生物、淋溶作用及土壤本身性質(zhì)的影響，土壤抗蝕性能也存在明顯的空間差異。由表2～4可以看出， 在0～30 cm土層，土壤抗蝕性指標(biāo)隨土層深度增加變化規(guī)律不明顯，且變異系數(shù)差異較大，其中土壤團(tuán)粒類指標(biāo)粒徑>0.25 mm水穩(wěn)性團(tuán)聚體含量、粒徑>0.5 mm水穩(wěn)性團(tuán)聚體含量、粒徑>0.25 mm非水穩(wěn)性團(tuán)聚體含量、水穩(wěn)性指數(shù)、團(tuán)聚體平均質(zhì)量直徑、團(tuán)聚體破壞率的變異系數(shù)分別為0.85%～67.46%,1.53%～98.63%,0.48%～23.02%,19.58%～82.09%,9.53%～120.2%和11.65%～123.54%；土壤物理類指標(biāo)非毛管孔隙度、毛管孔隙度、土壤體積質(zhì)量、吸濕水含量的變異系數(shù)分別為7.50%～26.06%，5.08%～50.79%，4.15%～45.69%和2.16%～44.48%；土壤有機(jī)膠體類指標(biāo)有機(jī)質(zhì)含量的變異系數(shù)為17.86%～93.63%。

2.3 色季拉山土壤抗蝕性最優(yōu)指標(biāo)選取及評(píng)價(jià)模型的建立

2.3.1 最優(yōu)指標(biāo)的選取 由于各抗蝕性指標(biāo)間信息重疊，相互間又有關(guān)聯(lián)，因此采用主成分分析法，在不損失較多信息的基礎(chǔ)上以較少主成分來(lái)代替原有指標(biāo)，通過(guò)構(gòu)建土壤抗蝕性評(píng)價(jià)模型，并計(jì)算綜合主成分值[28]，以評(píng)價(jià)色季拉山西坡不同海拔高度土壤的抗蝕性能。本研究選取11個(gè)土壤抗蝕性指標(biāo)，分別對(duì)0～10,10～20,20～30 cm土層進(jìn)行主成分分析,結(jié)果見(jiàn)表5。

表5 西藏色季拉山不同土層土壤抗蝕性指標(biāo)因子旋轉(zhuǎn)后的荷載和累計(jì)方差貢獻(xiàn)率Table 5 Load and cumulative variance contribution rate after factor rotating of soil anti-erodibility index at different soil layers

由表5可知，不同土層各提煉出3個(gè)主成分，其累積方差貢獻(xiàn)率分別為88.41%，90.61%和86.47%，信息損失量為9.39%～13.53%，滿足主成分分析的要求。但由于不同土層各指標(biāo)對(duì)主成分貢獻(xiàn)大小不同，造成指標(biāo)評(píng)價(jià)體系的垂直空間差異。由表5可以看出，在0～10 cm土層，第1主成分中粒徑>0.5 mm水穩(wěn)性團(tuán)聚體含量、土壤團(tuán)聚體平均質(zhì)量直徑、團(tuán)聚體破壞率的荷載值較大，說(shuō)明研究區(qū)0～10 cm土層粒徑>0.5 mm水穩(wěn)性團(tuán)聚體含量、土壤團(tuán)聚體平均質(zhì)量直徑指標(biāo)值越大，團(tuán)聚體破壞率越小，抗蝕性越強(qiáng)；第2主成分主要解釋了粒徑>0.25 mm土壤非水穩(wěn)性團(tuán)聚體含量、土壤體積質(zhì)量、吸濕水含量等信息;第3主成分主要解釋了非毛管孔隙度、毛管孔隙度信息。在10～20 cm土層，第1主成分中粒徑>0.25 mm土壤非水穩(wěn)性團(tuán)聚體含量、團(tuán)聚體破壞率、土壤體積質(zhì)量、吸濕水含量的荷載絕對(duì)值較高，說(shuō)明10～20 cm土層粒徑>0.25 mm非水穩(wěn)性團(tuán)聚體含量、吸濕水含量越高，團(tuán)聚體破壞率、土壤體積質(zhì)量越小，土壤抗蝕性越強(qiáng)；第2主成分中，粒徑>0.5mm水穩(wěn)性團(tuán)聚體含量、土壤團(tuán)聚體平均質(zhì)量直徑等指標(biāo)的荷載值較大；第3主成分主要解釋了非毛管孔隙度、毛管孔隙度等信息。在20～30 cm土層，第1主成分主要解釋了粒徑>0.25 mm水穩(wěn)性團(tuán)聚體含量、粒徑>0.5 mm水穩(wěn)性團(tuán)聚體含量、土壤團(tuán)聚體平均質(zhì)量直徑等信息，表明20～30 cm土層中粒徑>0.25 mm水穩(wěn)性團(tuán)聚體含量、粒徑>0.5 mm水穩(wěn)性團(tuán)聚體含量、土壤團(tuán)聚體平均質(zhì)量直徑越大，土壤抗蝕性越強(qiáng)；第2主成分主要解釋了土壤體積質(zhì)量、吸濕水含量、有機(jī)質(zhì)含量等信息;第3主成分主要解釋了土壤水穩(wěn)性指數(shù)的信息。

2.3.2 土壤抗蝕性綜合評(píng)價(jià)模型的構(gòu)建 以每個(gè)主成分對(duì)應(yīng)特征值占所提取主成分總特征值的比重作為權(quán)重，得到綜合主成分模型為：

式中：Wi、Fi分別為第i個(gè)因子的權(quán)重和隸屬度，n為參評(píng)因子數(shù)[29]。

由上述模型最終得到0～10 cm、10～20 cm和20～30 cm土層土壤抗蝕性能的評(píng)價(jià)模型分別為：

Y0～10=0.589Y1+0.280Y2+0.131Y3,Y10～20=0.631Y1+0.265Y2+0.104Y3,Y20～30=0.664Y1+0.186Y2+0.150Y3。

式中：Y1、Y2和Y3分別表示不同土層第1，2，3主成分值。

據(jù)此計(jì)算色季拉山不同海拔土壤抗蝕性能的綜合主成分值為：

Y=Y0～10+Y10～20+Y20～30。

通過(guò)綜合主成分評(píng)價(jià)模型,可以得出色季拉山不同海拔土壤綜合抗蝕性強(qiáng)弱順序?yàn)?3 600 m(土壤抗蝕性能的綜合主成分值為1.94，下同)>3 800 m(1.79)>3 400 m(1.59)>4 400 m(-0.24)>4 000 m(-0.50)>4 600 m(-1.01)>3 200 m(-1.33)>4 200 m(-2.23)。

3 討論與結(jié)論

不同海拔高度由于氣候條件、植被及土壤類型的差異，造成土壤理化性質(zhì)的不同，從而使其土壤抗蝕性能也具有明顯的空間差異。本研究結(jié)果表明，在色季拉山西坡海拔3 200～4 600 m區(qū)域，用于土壤抗蝕性能評(píng)價(jià)的土壤理化指標(biāo)中，除土壤體積質(zhì)量、團(tuán)聚體破壞率及毛管孔隙度、非毛管孔隙度外，其余土壤抗蝕性指標(biāo)值總體隨海拔上升表現(xiàn)出先增大后減小再緩慢增大的空間分異特征。具體可分為3個(gè)變化區(qū)間，即當(dāng)海拔為3 200～3 600 m時(shí)，隨海拔升高而增加；當(dāng)海拔為3 600～4 200 m時(shí)，隨海拔升高而減??；當(dāng)海拔為4 200～4 600 m時(shí)，隨海拔升高而緩慢增加。

本研究結(jié)果表明，在色季拉山西坡海拔3 200～3 600 m區(qū)域，土壤抗蝕性隨海拔的升高而增加，這與廖超林等[30]的研究結(jié)果相一致，主要由于該海拔區(qū)間水熱條件充足，植被類型豐富，隨海拔上升人為擾動(dòng)減少，植被根系的纏繞、固結(jié)對(duì)土壤的改良作用明顯。而在海拔3 600～4 200 m區(qū)域，隨海拔的升高土壤抗蝕性能逐漸減弱，這主要與不同研究區(qū)間氣候、土壤及植被類型等的差異較大有關(guān)。研究區(qū)地處青藏高原東南部，區(qū)內(nèi)隨著海拔增高，氣候從亞高山寒溫帶向高山寒帶過(guò)渡，急尖長(zhǎng)苞冷杉、林芝云杉等針葉物種所占比例增多，在低溫潮濕條件下，隨著表層枯落物分解，有機(jī)酸增多，腐殖質(zhì)、鐵鋁氧化物淋溶作用增強(qiáng)，形成氧化硅相對(duì)富集的灰化層，導(dǎo)致土壤結(jié)構(gòu)穩(wěn)定性降低，抗蝕性能減弱。當(dāng)海拔為4 200～4 600 m時(shí)，隨海拔增加土壤抗蝕性總體增大，而這與干冷氣候條件下土壤微生物活動(dòng)受到限制及表層土壤有機(jī)質(zhì)分解緩慢而積累較易有關(guān)。茍小林等[31]的研究發(fā)現(xiàn)，季節(jié)性凍融會(huì)促進(jìn)土壤可溶性有機(jī)碳的淋溶，在非生長(zhǎng)季節(jié)最高淋溶輸入量可達(dá)7.30 kg/hm2，從而對(duì)深層土壤理化性質(zhì)有明顯改良，土壤結(jié)構(gòu)穩(wěn)定性、抗蝕性能增強(qiáng)。

本研究中，比較色季拉山土壤抗蝕性指標(biāo)的空間變異系數(shù)可以發(fā)現(xiàn)，在海拔3 400～3 800 m區(qū)域，土壤團(tuán)粒類指標(biāo)的變異系數(shù)相對(duì)較低，而在海拔3 600 m區(qū)域土壤基本物理類指標(biāo)的變異系數(shù)相對(duì)較高；在海拔4 200～4 600 m區(qū)域，土層間土壤團(tuán)粒類指標(biāo)變異系數(shù)較高，而土壤基本物理類指標(biāo)的變異系數(shù)卻相對(duì)較??；有機(jī)膠體類指標(biāo)有機(jī)質(zhì)含量的變異系數(shù)則分別在海拔4 600和4 200 m處出現(xiàn)最小值和最大值。

主成分分析發(fā)現(xiàn)，色季拉山不同土層間描述土壤抗蝕性的最優(yōu)指標(biāo)存在較大差異，在0～10 cm土層，以土壤團(tuán)聚體平均質(zhì)量直徑、團(tuán)聚體破壞率和粒徑>0.5 mm水穩(wěn)性團(tuán)聚體含量3個(gè)指標(biāo)表現(xiàn)最好；在10～20 cm土層，土壤吸濕水含量、體積質(zhì)量和粒徑>0.25 mm非水穩(wěn)性土壤團(tuán)聚體含量3個(gè)指標(biāo)具有較強(qiáng)的代表性；在20～30 cm土層，粒徑>0.25 mm水穩(wěn)性團(tuán)聚體含量、粒徑>0.5 mm水穩(wěn)性土壤團(tuán)聚體含量、團(tuán)聚體平均質(zhì)量直徑3個(gè)指標(biāo)表現(xiàn)出較好的評(píng)價(jià)性能。

本研究構(gòu)建了色季拉山不同海拔區(qū)域土壤抗蝕性綜合評(píng)價(jià)模型，由其可得不同海拔高度下土壤抗蝕性由強(qiáng)到弱依次表現(xiàn)為3 600 m>3 800 m>3 400 m>4 400 m>4 000 m>4 600 m>3 200 m>4 200 m。

[參考文獻(xiàn)]

[1] Chen J,He B,Wang X,et al.The effects of Herba Andrographitis hedgerows on soil ero-dibility and fractal features on sloping cropland in the Three Gorges Reservoir Area [J].Environmental Science and Pollution Research,2013,20(10):7063-7070.

[2] Cotler H,Ortega-Larrocea M P.Effects of land use on soil erosion in a tropical dry forest ecosystem,Chamela watershed,Mexico [J].Catena,2006,65(2):107-117.

[3] 王政權(quán),王慶成.森林土壤物理性質(zhì)的空間異質(zhì)性研究 [J].生態(tài)學(xué)報(bào),2000,20(6):945-950.

Wang Z Q,Wang Q C.The spatial heterogeneity of soil physical properties in forests [J].Acta Ecologica Sinica,2000,20(6):945-950.

[4] Biederbeck V O,Curtin D,Bouman O T,et al.Soil microbial and biochemical properties after ten years of fertilization with urea and anhydrous ammonia [J].Canadian Journal of Soil Science,1996,76(1):7-14.

[5] 黃建輝.物種多樣性的空間格局及其形成機(jī)制初探 [J].生物多樣性,1994,48(2):103-107.

Huang J H.The spatial pattern of species diversity and its forming mechanism [J].Biodiversity Science,1994,48(2):103-107.

[6] Begon M,Harper J L,Townsend C R.Ecology:individuals populations and communities [M].2nd ed.Boston:Blackwell Scientific Publications,2000.

[7] 黃建輝,白永飛,韓興國(guó).物種多樣性與生態(tài)系統(tǒng)功能:影響機(jī)制及有關(guān)假說(shuō) [J].生物多樣性,2001,9(1):1-7.

Huang J H,Bai Y F,Han X G.Effects of species diversity on ecosystem functioning:mechanisms and hypotheses [J].Biodiversity Science,2001,9(1):1-7.

[8] 楊梅學(xué),姚檀棟,何元慶.青藏高原土壤水熱分布特征及凍融過(guò)程在季節(jié)轉(zhuǎn)換中的作用 [J].山地學(xué)報(bào),2002,20(5):553-558.

Yang M X,Yao C D,He Q Y.The role of soil moisture-energy distribution and melting-freezing processes on seasonal shift in Tibetan plateau [J].Mountain Research,2002,20(5):553-558.

[9] 周幼吾,郭東信,邱國(guó)慶,等.中國(guó)凍土 [M].北京:科學(xué)出版社,2000.

Zhou Y W,Guo D X,Qiu G Q,et al.China’s permafrost [M].Beijing:Science Press,2000.

[10] 趙 勇,錢永甫.青藏高原地表熱力異常與我國(guó)江淮地區(qū)夏季降水的關(guān)系 [J].大氣科學(xué),2007, 31(1):145-154.

Zhao Y,Qian Y F.Relationships between the surface thermal anomalies in the Tibetan plateau and the rainfall in the Jianghuai area in summer [J].Chinese Journal of Atmospheric Sciences,2007,31(1):145-154.

[11] 范繼輝,魯旭陽(yáng),王小丹.藏北高寒草地土壤凍融循環(huán)過(guò)程及水熱分布特征 [J].山地學(xué)報(bào),2014,32(4):385-392.

Fan J H,Lu X Y,Wang X D.The freezing-thawing processes and soil moisture-energy distribution in permafrost active layer,northern Tibet [J].Mountain Research,2014,32(4):385-392.

[12] 孫鴻烈,鄭 度,姚檀棟,等.青藏高原國(guó)家生態(tài)安全屏障保護(hù)與建設(shè) [J].地理學(xué)報(bào),2012,67(1):3-12.

Sun H L,Zhen D,Yao C D,et al.Protection and construction of the national ecological security [J].Acta Geographica Sinica,2012,67(1):3-12.

[13] 劉合滿,曹麗花,張 華,等.色季拉山山地酸性棕壤土壤氮素的分布特征 [J].中南林業(yè)科技大學(xué)學(xué)報(bào),2013,33(10):126-129.

Liu H M,Cao L H,Zhang H,et al.Distribution characteristics of soil nitrogen in mountain acid brownsoil in Sejila Mountain,Tibet [J].Journal of Central South University of Forestry & Technology,2013,33(10):126-129.

[14] 李為虎,楊小林,馬和平.色季拉山植物群落物種多樣性研究 [J].西南林學(xué)院學(xué)報(bào),2008,28(5):3-6.

Li W H,Yang X L,Ma H P.Species diversity of sejila mountain forest plant community [J].Journal of Southwest Forestry College,2008,28(5):3-6.

[15] 何建清,張格杰,張新軍,等.西藏色季拉山土壤放線菌的區(qū)系 [J].山地學(xué)報(bào),2008,26(4):503-508.

He J Q,Zhang G J,Zhang X J,et al.Study on actinomycetes distribution of Serkyimla Mountain in Tibet [J].Journal of Mountain Science,2008,26(4):503-508.

[16] 方江平.西藏色季拉山土壤的性狀與垂直分布 [J].山地學(xué)報(bào),1997(4):228-233.

Fang J P.Properties and vertical distribution of soil on Shergyl Mountain in Xizang [J].Mountain Research,1997(4):228-233.

[17] 張啟翔.中國(guó)觀賞植物種質(zhì)資源 [M].北京：中國(guó)林業(yè)出版社,2011.

Zhang Q X.China ornamental plant germplasm resources [M].Beijing:China Forestry Publishing House,2011.

[18] 柴 勇,樊國(guó)盛,李鄉(xiāng)旺,等.西藏色季拉山種子植物垂直帶譜的劃分與分布特點(diǎn)研究 [J].廣西植物,2004,24(2):107-112.

Chai Y,Fan S G,Li X W,et al.Study on vertical distributional belts and their floristic characters of seed plants from Shegyla Mountains of Xizang(Tibet),China [J].Guihaia,2004,24(2):107-112.

[19] 中國(guó)林業(yè)科學(xué)研究院林業(yè)研究所.中國(guó)森林土壤 [M].北京：科學(xué)出版社,1986.

Forestry Science Research Institute of Forestry Research Institute of China.China’s forest soil [M].Beijing:Science Press,1986.

[20] 水利水電部農(nóng)村水利水土保持司.水土保持實(shí)驗(yàn)規(guī)范 [S].北京:水利水電部,1988.

Department of Water Resources and Hydropower Rural Water Conservancy Department of Soil and Water Conservation.Soil and water conservation experiment specification [S].Beijing:The Ministry of Water Resources and Hydropower,1988.

[21] 劉國(guó)彬.黃土高原土壤抗沖性研究及有關(guān)問(wèn)題 [J].水土保持研究,1997(S1):91-101.

Liu G B.Soil anti-scourability research and its perspectives in Loess Plateau [J].Research of Soil and Water Conservation,1997(S1):91-101.

[22] 劉光崧.土壤理化分析與剖面描述 [M].北京:中國(guó)標(biāo)準(zhǔn)出版社,1996.

Liu G S.Soil physical and chemical analysis and section describes [M].Beijing:China Standard Press,1996.

[23] 中國(guó)林業(yè)科學(xué)研究院林業(yè)研究所森林土壤研究實(shí)驗(yàn)室.LY/T 1215-1999 森林土壤水分-物理性質(zhì)的測(cè)定 [S].北京：中國(guó)標(biāo)準(zhǔn)出版社,1999.

China institute of forestry science research institute of forestry forest soil research laboratory.LY/T 1215-1999 Determination of forest soil water Physical properties [S].Beijing：China Standard Publishing House,1999.

[24] 中國(guó)農(nóng)學(xué)會(huì).土壤農(nóng)業(yè)化學(xué)分析方法 [M].北京:中國(guó)農(nóng)業(yè)科技出版社,2000.

Chinese Association of Agricultural Science Societies.Soil agricultural chemical analysis methods [M].Beijing:China Agricultural Science and Technology Publishing House,2000.

[25] 史曉梅,史東梅,文卓立.紫色土丘陵區(qū)不同土地利用類型土壤抗蝕性特征研究 [J].水土保持學(xué)報(bào),2007,21(4):63-66.

Shi X M,Shi D M,Wen Z L.Study on soil anti-erodibility of different land utilization types in purple soil hilly region [J].Journal of Soil and Water Conservation,2007,21(4):63-66.

[26] 邱莉萍,張興昌,張晉愛(ài).黃土高原長(zhǎng)期培肥土壤團(tuán)聚體中養(yǎng)分和酶的分布 [J].生態(tài)學(xué)報(bào),2006,26(2):364-372.

Qiu L P,Zhang X C,Zhang J A.Distribution of nutrients and enzymes in Loess Plateau soil aggregates after longterm fertilization [J].Acta Ecologica Sinica,2006,26(2):364-372.

[27] 李立青,尹澄清,何慶慈,等.城市降水徑流的污染來(lái)源與排放特征研究進(jìn)展 [J].水科學(xué)進(jìn)展,2006,17(2):288-294.

Li L Q,Yin C Q,He Q C,et al.Progress in research on the sources and characterization of urban rainfall runoff [J].Advances in Water Science,2006,17(2):288-294.

[28] 白秀梅,韓有志,郭漢清.關(guān)帝山不同植被恢復(fù)類型土壤抗蝕性研究 [J].水土保持學(xué)報(bào),2014,28(2):79-84.

Bai X M,Han Y Z,Guo H Q.Study on soil anti-erodibility of different vegetation restoration types in Guandi Mountain [J].Journal of Soil and Water Conservation,2014,28(2):79-84.

[29] 黃 婷,岳西杰,葛璽祖,等.基于主成分分析的黃土溝壑區(qū)土壤肥力質(zhì)量評(píng)價(jià):以長(zhǎng)武縣耕地土壤為例 [J].干旱地區(qū)農(nóng)業(yè)研究,2010,28(3):141-147.

Huang T,Yue X J,Ge X Z,et al.Effects of exogenous spermidine on photosynthesis of tomato seedlings under drought stress [J].Agricultural Research in the Arid Areas,2010,28(3):141-147.

[30] 廖超林,何毓蓉,徐 佩.泥石流源地土壤團(tuán)聚體抗蝕特征研究:以蔣家溝為例 [J].地球與環(huán)境,2005,33(4):65-69.

Liao C L,He Y R,Xu P.A Study on antierodibility of soil aggregates in debris flow source:a case study in Jiangjiagou [J].Earth and Environment,2005,33(4):65-69.

[31] 茍小林,吳福忠,楊萬(wàn)勤,等.季節(jié)性凍融格局變化對(duì)高山森林土壤DOC淋洗的影響 [J].水土保持學(xué)報(bào),2013,27(6):201-210.

Gou X L,Wu F Z,Yang W Q,et al.Effect of changes in seasonal freeze-thaw pattern on DOC loss from leaching in the alpine forest soil [J].Journal of Soil and Water Conservation,2013,27(6):201-210.