袁 方，張振師，張 朋，卜崇峰，2

（1.西北農(nóng)林科技大學(xué) 水土保持研究所，陜西 楊凌712100；2.中國科學(xué)院水利部水土保持研究所，陜西 楊凌712100；3.中國電建集團(tuán)西北勘測設(shè)計研究院有限公司，西安710065）

生物結(jié)皮（Biological soil crusts）是土壤顆粒與微生物、藻類、地衣以及苔蘚等組成的復(fù)雜聚合體，是干旱半干旱廣泛存在的活性地被物［1］。它的存在可以富集土壤養(yǎng)分，提高土壤抗蝕性［2］，降低土壤的風(fēng)蝕、水蝕量，影響水分再分配［3－4］，具有重要的生態(tài)功能［5］和良好的水土保持作用。通徑分析是簡單相關(guān)分析的繼續(xù)，是把自變量與因變量之間的相關(guān)性分解為該自變量對因變量的直接影響和通過其他相關(guān)自變量對因變量間接影響的分析方法［6］。其不僅具有回歸與相關(guān)分析的作用，且進(jìn)一步地將自變量與因變量的相關(guān)統(tǒng)計量剖分為直接通徑系數(shù)和間接通徑系數(shù)，揭示影響現(xiàn)象變異的主要因素所產(chǎn)生的直接和間接作用，從而區(qū)分自變量的相對重要性及其與因變量的關(guān)系，為研究和分析的問題提供更全面更完善的決策依據(jù)［7］。前人對生物結(jié)皮空間分布及其影響因子的研究大多通過野外調(diào)查，對所得各類指標(biāo)進(jìn)行直觀的統(tǒng)計分析［8－9］，但鮮有學(xué)者利用通徑分析研究各因子對生物結(jié)皮的直接或間接影響及其作用路徑。

黃土高原水蝕風(fēng)蝕交錯區(qū)位于黃土高原北部，屬于干旱半干旱地區(qū)草原帶［10］，是黃土高原水土流失最為嚴(yán)重的地區(qū)［11］，其環(huán)境條件類型復(fù)雜多變、景觀破碎，水分條件比較惡劣。自退耕還林（草）工程實施以來，黃土高原地區(qū)生物結(jié)皮大面積發(fā)育［8］，其在防治水土流失和促進(jìn)生態(tài)演替方面作用顯著。而生物結(jié)皮在自然條件下的發(fā)育十分緩慢，容易受到環(huán)境因素（土壤、水分條件、植被、海拔、坡向等［12－13］）和人類干擾活動［14］的影響，且各影響因素對生物結(jié)皮的影響程度及作用存在差異，故通過生物結(jié)皮生長狀況的野外調(diào)查，了解生物結(jié)皮的空間分布規(guī)律及其影響因子，明確各影響因子的具體作用路徑，以期為更好地促進(jìn)生物結(jié)皮的形成和發(fā)育奠定基礎(chǔ)。

為此，本研究選取黃土高原北部水蝕風(fēng)蝕交錯區(qū)六道溝小流域的一處完整梁峁坡面，通過全面調(diào)查和測算分析，應(yīng)用SPSS 12.0及 Excel 2007軟件，采用通徑分析的方法研究生物結(jié)皮空間分布與環(huán)境因子之間的關(guān)系，旨在為野外生物結(jié)皮發(fā)育的環(huán)境評價和快速培育提供理論支持。

1 材料與方法

1.1 研究區(qū)概況

研究區(qū)位于陜西省神木縣以西14km處的六道溝小流域（110°21′—110°23′E，38°46′—38°51′N），海拔1 094～1 274m，流域面積6.89km2，主溝道南北走向，長4.2 km，屬窟野河二級支流。地處黃土高原向毛烏素沙地、森林草原向干旱草原的過渡地帶，屬水蝕作用向風(fēng)蝕作用過渡的典型水蝕風(fēng)蝕交錯帶，年侵蝕模數(shù)達(dá)10 000t/km2［15］。氣候?qū)僦袦貛О敫珊禋夂颍杭靖珊瞪儆甓囡L(fēng)沙，風(fēng)蝕嚴(yán)重，夏秋多雨且暴雨頻繁，水蝕強(qiáng)烈。年均降雨量400mm，其中6—9月份的降雨量約占全年的70%～80%，年均氣溫7～9℃，盛行西北風(fēng)，年平均風(fēng)速為2.2m/s［16］。流域內(nèi)的主要植被有檸條（Caragana korshinskiiKom.）、長芒草（StipabungeanaTrin.）、苜蓿（MedicagosativaLinn）、達(dá) 烏 里 胡 枝 子 （Lespedeza davurica（Laxm.）Schindl.）、阿爾泰狗娃花（Heteropappusaltaicus（Willd）Novopokr）和茭蒿（ArtemisiagiraldiiPamp）等。流域東側(cè)主要為黃土覆蓋，占流域面積的86.5%；西側(cè)多為固定沙丘，占流域面積的13.5%［17］。

1.2 研究方法

本研究在流域內(nèi)選定一處撂荒30a以上的典型的梁峁坡面，于2014年6—10月開展生物結(jié)皮空間分布的系統(tǒng)調(diào)查工作。其中于6月，7月份對生物結(jié)皮的蓋度、厚度和抗剪強(qiáng)度進(jìn)行調(diào)查。利用etrex HD型手持GPS對野外每個調(diào)查點進(jìn)行定位，利用方格網(wǎng)法測算該點周圍半徑10m左右范圍內(nèi)生物結(jié)皮的蓋度（%），厚度（mm）用游標(biāo)卡尺進(jìn)行原位測定，重復(fù)3～5次；抗剪強(qiáng)度（kPa）用袖珍剪力儀（BWT2XZJL）原位測定，重復(fù)3～5次；于8月中旬到9月上旬對植被類型和植被蓋度進(jìn)行調(diào)查；于10月上旬（雨季末）在調(diào)查區(qū)內(nèi)選擇幾處典型樣地，對樣地內(nèi)的土壤進(jìn)行采樣分析，取表層0—5cm土層，帶回室內(nèi)105℃下烘24h，測算土壤含水率（%）。利用GIS插值功能將所有野外調(diào)查點及調(diào)查數(shù)據(jù)進(jìn)行制圖，共得到1 342個樣本的結(jié)皮蓋度和厚度、地形、土壤、植被等指標(biāo)。

數(shù)據(jù)的統(tǒng)計分析利用SPSS 12.0及 Excel 2007軟件進(jìn)行通徑分析。通徑分析是研究多個自變量與因變量間多重線性關(guān)系的一種重要的統(tǒng)計方法，是相關(guān)分析和回歸分析的結(jié)合，將自變量與因變量的簡單相關(guān)系數(shù)分解為直接通徑系數(shù)和間接通徑系數(shù)，分別表示某一自變量對因變量的直接影響及通過其他自變量對因變量的間接影響，決定系數(shù)表示自變量對因變量的決定程度，通過比較各變量的決定系數(shù)大小，對各變量的作用進(jìn)行排序，為統(tǒng)計決策提供更全面更完善的依據(jù)［18］；袁志發(fā)等［19］提出決策系數(shù)的概念，其反映了自變量通過其他自變量的相關(guān)網(wǎng)對因變量的綜合決定作用，指出在復(fù)雜的信息路徑中，選擇什么樣的路徑對因變量最好。

2 結(jié)果與分析

2.1 生物結(jié)皮的空間分布概況

生物結(jié)皮蓋度的空間分布具有明顯的分異性和連續(xù)性。結(jié)皮的平均蓋度在30%以上，主要分布在沙土區(qū)，呈成片的連續(xù)分布；而在黃土區(qū)生物結(jié)皮則呈離散的零星分布，蓋度大多在20%以下，且主要分布在坡地的邊緣和末端。沙土區(qū)生物結(jié)皮的空間分布具有明顯的發(fā)源地，大都處在沙土和黃土交界處的低洼處。相對來說，生物結(jié)皮的厚度和抗剪強(qiáng)度的空間變異性卻并不大。沙土區(qū)和黃土區(qū)生物結(jié)皮的平均厚度分別為12.8mm和13.1mm，平均抗剪強(qiáng)度分別為4.76kg/cm2，4.96kg/cm2。t檢驗結(jié)果表明：黃土區(qū)生物結(jié)皮的厚度及平均抗剪強(qiáng)度和沙土區(qū)的差異并不顯著（p＞0.05）。

2.2 生物結(jié)皮蓋度與環(huán)境因子間的關(guān)系

黃土高原水蝕風(fēng)蝕交錯區(qū)復(fù)雜的環(huán)境特征影響著生物結(jié)皮的生長發(fā)育，結(jié)皮蓋度是描述生物結(jié)皮發(fā)育狀況的一個重要指標(biāo)。為研究結(jié)皮蓋度與環(huán)境因子之間的關(guān)系，根據(jù)逐步回歸結(jié)果，選取沙土、沉積物運移指數(shù)、坡向、沙蒿蓋度、海拔及土壤含水率作為影響因變量生物結(jié)皮蓋度的自變量，進(jìn)行相關(guān)及通徑分析。

2.2.1 生物結(jié)皮蓋度及各環(huán)境因子的相關(guān)分析對結(jié)皮蓋度及沙土、沉積物運移指數(shù)、坡向、沙蒿蓋度、海拔及土壤含水率采用SPSS 12.0軟件進(jìn)行相關(guān)分析，得出上述環(huán)境因子與結(jié)皮蓋度的相關(guān)及回歸系數(shù)（表1）。

表1 環(huán)境因子與結(jié)皮蓋度的相關(guān)與回歸分析結(jié)果

由表1可以看出：6個環(huán)境因子與結(jié)皮蓋度相關(guān)系數(shù)的絕對值大小依次為：沙土＞沙蒿蓋度＞海拔＞坡向＞沉積物運移指數(shù)＞土壤含水率，且沙土和沙蒿蓋度與結(jié)皮蓋度呈正相關(guān)關(guān)系，其余均與結(jié)皮蓋度呈負(fù)相關(guān)關(guān)系，除土壤含水率以外各環(huán)境因子均與結(jié)皮蓋度的相關(guān)性達(dá)到極顯著水平（p＜0.01）。此外，各環(huán)境因子之間也存在顯著的相關(guān)性，即共線性；各環(huán)境因子與結(jié)皮蓋度之間標(biāo)準(zhǔn)化回歸方程達(dá)到極顯著水平（p＜0.01），R2＝0.706，這說明結(jié)皮蓋度關(guān)于上述6個環(huán)境因子的通徑分析是必要的。

2.2.2 生物結(jié)皮蓋度及各環(huán)境因子的通徑分析 表2和圖1a展示了不同環(huán)境因子對結(jié)皮蓋度的直接和間接影響。從表2可以看出：（1）沙土對結(jié)皮蓋度的直接影響最大（0.720），土壤含水率（0.526）及沙蒿蓋度（0.235）次之，即從直接影響看沙土、土壤含水率及沙蒿蓋度是影響結(jié)皮蓋度的主要因子。（2）沙蒿蓋度通過沙土對結(jié)皮蓋度的間接影響最大（0.422），土壤含水率通過沙土對結(jié)皮蓋度的間接影響（－0.418）次之，且影響為負(fù)。（3）沙土通過其他環(huán)境因子對結(jié)皮蓋度的間接影響為0.067（小于直接影響0.720），說明沙土是影響結(jié)皮蓋度的最直接因子；沙蒿蓋度通過其他環(huán)境因子對結(jié)皮蓋度的間接影響為0.403，直接影響為0.235，且其僅通過沙土對結(jié)皮蓋度的間接影響就已經(jīng)達(dá)到0.422，即沙蒿蓋度對結(jié)皮蓋度的間接影響大于直接影響且以通過沙土的間接影響最為主要；土壤含水率對結(jié)皮蓋度的直接影響為0.526，通過其他環(huán)境因子的間接影響為－0.612（小于0），導(dǎo)致土壤含水率與結(jié)皮蓋度呈負(fù)相關(guān)關(guān)系（相關(guān)系數(shù)為－0.086），故實際上低估了土壤含水率對結(jié)皮蓋度的直接影響。

2.2.3 各環(huán)境因子對結(jié)皮蓋度的決定系數(shù)和決策系數(shù) 分析各環(huán)境因子對結(jié)皮蓋度的決定系數(shù)和決策系數(shù)（表3）可知：（1）6個環(huán)境因子的直接決定系數(shù)的大小依次為沙土＞土壤含水率＞海拔＞坡向＞沉積物運移指數(shù)＞沙蒿蓋度，即沙土及土壤含水率對結(jié)皮蓋度的直接決定作用最大；（2）共同決定系數(shù)中，沙土和土壤含水率的共同決定系數(shù)最大（－0.44），其次是沙土和沙蒿蓋度（0.198）以及沙土和海拔（0.176）。（3）從決策系數(shù)來看，6個環(huán)境因子決策系數(shù)絕對值的大小依次為：沙土＞土壤含水率＞海拔＞沙蒿蓋度＞坡向＞沉積物運移指數(shù)，即沙土、土壤含水率、海拔和沙蒿蓋度是主要的決策變量，其中除土壤含水率的決策系數(shù)為負(fù)值之外，其余均為正值。

表2 各環(huán)境因子對結(jié)皮蓋度的通徑系數(shù)

圖1 生物結(jié)皮蓋度、厚度與環(huán)境因子間的通徑關(guān)系

表3 各環(huán)境因子對結(jié)皮蓋度的決定系數(shù)和決策系數(shù)

因此，綜上所述，沙土、土壤含水率、海拔及沙蒿蓋度為影響結(jié)皮蓋度的主要因素，但各自對生物結(jié)皮蓋度的主要作用路徑不同；沙土的直接作用（0.72），沙蒿蓋度通過沙土的間接作用（0.422）是其主要的作用路徑，兩者都對結(jié)皮蓋度有促進(jìn)作用；土壤含水率的直接正作用（0.526）和間接負(fù)作用（－0.612）相當(dāng)，因而總體貢獻(xiàn)并不明顯；海拔的直接和間接作用均為負(fù)（－0.358和－0.244），抑制作用明顯。

2.3 生物結(jié)皮厚度與環(huán)境因子間的關(guān)系

結(jié)皮厚度是描述和衡量生物結(jié)皮發(fā)育狀況的另外一個重要指標(biāo)。為研究結(jié)皮厚度與環(huán)境因子之間的關(guān)系，根據(jù)逐步回歸結(jié)果，選取坡向、沉積物運移指數(shù)、海拔、沙蒿蓋度、小葉楊蓋度及硬質(zhì)早熟禾蓋度作為影響因變量生物結(jié)皮厚度的自變量，進(jìn)行相關(guān)及通徑分析。

2.3.1 結(jié)皮厚度及各環(huán)境因子的相關(guān)分析 對結(jié)皮厚度及坡向、沉積物運移指數(shù)、海拔、沙蒿蓋度、小葉楊蓋度及硬質(zhì)早熟禾蓋度采用SPSS 12.0軟件進(jìn)行相關(guān)和回歸分析，得出上述6個環(huán)境因子與結(jié)皮厚度的相關(guān)及回歸系數(shù)（表4）。

由表4可以看出，6個環(huán)境因子與結(jié)皮厚度相關(guān)系數(shù)的絕對值大小依次為：坡向＞沙蒿蓋度＞海拔＞小葉楊蓋度＞硬質(zhì)早熟禾蓋度＞沉淀物運移指數(shù)，且沙蒿蓋度、小葉楊蓋度及硬質(zhì)早熟禾蓋度與結(jié)皮厚度呈正相關(guān)關(guān)系，坡向、海拔及沉積物運移指數(shù)與結(jié)皮厚度均呈負(fù)相關(guān)關(guān)系，各環(huán)境因子均與結(jié)皮厚度的相關(guān)性達(dá)到極顯著水平（p＜0.01）。此外，各環(huán)境因子之間也存在著不同程度的顯著相關(guān)性，即共線性；6個環(huán)境因子與結(jié)皮厚度之間標(biāo)準(zhǔn)化回歸方程達(dá)到極顯著水平（p＜0.01），R2＝0.552，說明結(jié)皮厚度關(guān)于上述6個環(huán)境因子的通徑分析仍然是有必要的。

表4 環(huán)境因子與結(jié)皮厚度的相關(guān)與回歸分析

2.3.2 生物結(jié)皮厚度及各環(huán)境因子的通徑分析 表5和圖1b展示了不同環(huán)境因子對結(jié)皮厚度的直接和間接影響。從表5可以看出：（1）直接通徑系數(shù)的大小依次為沉積物運移指數(shù)（－0.237）＜坡向（－0.215）＜海拔（－0.196）＜0＜硬質(zhì)早熟禾蓋度（0.193）＜小葉楊蓋度（0.233）＜沙蒿蓋度（0.329），即沙蒿蓋度對結(jié)皮厚度的直接影響最大，沙蒿蓋度越大，結(jié)皮厚度越大。（2）硬質(zhì)早熟禾蓋度通過海拔對結(jié)皮厚度的間接影響（0.109）最大，海拔通過硬質(zhì)早熟禾蓋度對結(jié)皮厚度的間接負(fù)影響（－0.108）最大。（3）坡向通過其他環(huán)境因子對結(jié)皮厚度的間接負(fù)影響（－0.311）最大，其中以通過沙蒿蓋度的間接影響（－0.103）最為突出，其導(dǎo)致坡向與結(jié)皮厚度之間呈負(fù)相關(guān)關(guān)系（相關(guān)系數(shù)為－0.526）。

表5 各環(huán)境因子對結(jié)皮厚度的通徑系數(shù)

2.3.3 各環(huán)境因子對結(jié)皮厚度的決定系數(shù)和決策系數(shù) 分析各個環(huán)境因子對結(jié)皮厚度的決定系數(shù)和決策系數(shù)（表6）可以看出：（1）6個環(huán)境因子直接決定系數(shù)的大小依次為沙蒿蓋度＞沉積物運移指數(shù)＞小葉楊蓋度＞坡向＞海拔＞硬質(zhì)早熟禾蓋度，即沙蒿蓋度對結(jié)皮厚度的直接決定作用最大（0.108），沉積物運移指數(shù)和小葉楊蓋度次之，分別為0.056，0.054；（2）共同決定系數(shù)中，沙蒿蓋度和坡向的共同決定系數(shù)最大（0.044），硬質(zhì)早熟禾蓋度與坡向和海拔的共同決定作用次之（均為0.042）；（3）從決策系數(shù)來看，6個環(huán)境因子決策系數(shù)的大小依次為：沙蒿蓋度＞坡向＞海拔＞小葉楊蓋度＞硬質(zhì)早熟禾蓋度＞沉積物運移指數(shù)，且各因子的決策系數(shù)均為正，對結(jié)皮厚度均起促進(jìn)作用，說明沙蒿蓋度、坡向及海拔是比較重要的環(huán)境決策變量，沉積物運移指數(shù)的作用最小；（4）6個環(huán)境因子對結(jié)皮厚度的總決定作用達(dá)55.2%（表6），說明這些環(huán)境因子是影響結(jié)皮厚度的主要因素，對結(jié)皮蓋度的影響占相對重要的地位，因此在培育生物結(jié)皮中應(yīng)予注意。

綜上所述，沙蒿蓋度、坡向及海拔是影響生物結(jié)皮空間分布和發(fā)育的關(guān)鍵因素，但三者對生物結(jié)皮蓋度的主要作用路徑不同，沙蒿蓋度的直接作用是主要的作用路徑，而坡向通過沙蒿蓋度、海拔通過硬質(zhì)早熟禾蓋度的間接作用是主要作用路徑。

表6 各環(huán)境因子對結(jié)皮厚度的決定系數(shù)和決策系數(shù)

3 結(jié)論與討論

本研究的結(jié)果表明，生物結(jié)皮在黃土區(qū)和沙土區(qū)的空間分布具有明顯差異。在沙土區(qū)呈連續(xù)分布，多發(fā)源于沙土和黃土交界處的低洼帶，平均蓋度達(dá)30%以上；而在黃土區(qū)呈離散的零星分布，蓋度大多在20%以下，這與前人的結(jié)論一致［8］，其原因可能是，由于低洼帶具有較好的水分條件，其有利于沙土區(qū)生物結(jié)皮的連片發(fā)育；而黃土區(qū)的植被條件較好，植被及其枯落物的覆蓋使生物結(jié)皮的發(fā)育空間受到限制［20］。

各環(huán)境因子對結(jié)皮蓋度都有不同程度的影響。各環(huán)境因子的決策系數(shù)除土壤含水率為負(fù)值之外，其余均為正值，而造成其決策系數(shù)為負(fù)的原因可能是：土壤含水率高的地方被更多的維管植物占據(jù)，若維管植物大量繁衍，將大大降低生物結(jié)皮的生存空間，使得結(jié)皮蓋度有所減少［12］。故土壤質(zhì)地（沙土）、土壤含水率、海拔及沙蒿蓋度可作為影響結(jié)皮蓋度的主要因子，四者是影響生物結(jié)皮空間分布和發(fā)育的關(guān)鍵因素，這與前人的結(jié)論一致［13］。土壤質(zhì)地（沙土）是影響結(jié)皮蓋度的直接因子，對其直接決定作用最大，沙蒿蓋度是影響結(jié)皮蓋度的間接因子，對其間接決定作用最大，兩者均對結(jié)皮蓋度起促進(jìn)作用，土壤含水率的直接正作用和間接負(fù)作用相當(dāng)，因而總體貢獻(xiàn)并不明顯；海拔的直接和間接作用均為負(fù)。主要因子對生物結(jié)皮蓋度的作用路徑不同，土壤質(zhì)地（沙土）的直接作用、沙蒿蓋度通過土壤質(zhì)地（沙土）的間接作用是其主要的作用路徑。

各環(huán)境因子對結(jié)皮厚度的影響程度不同。沙蒿蓋度、坡向及海拔可作為影響結(jié)皮厚度的主要因子。沙蒿蓋度是直接影響因子，對結(jié)皮厚度的直接決定作用最大，且沙蒿蓋度越大，結(jié)皮厚度越大，原因可能是沙蒿獨特的灌叢群落結(jié)構(gòu)為生物結(jié)皮的定居和增殖提供了相對穩(wěn)定的環(huán)境［20］，有利于生物結(jié)皮的發(fā)育。坡向及海拔是影響結(jié)皮厚度的間接因子，其間接決定作用大于直接決定作用，且生物結(jié)皮厚度陰坡大于陽坡，這與前人的結(jié)論一致［9］；海拔對結(jié)皮厚度起促進(jìn)作用，而前人就海拔對結(jié)皮厚度影響的研究還鮮有報道。但三者對生物結(jié)皮厚度的主要作用路徑不同，沙蒿蓋度的直接作用、坡向通過沙蒿蓋度及海拔通過硬質(zhì)早熟禾蓋度的間接作用是其主要的作用路徑。

此外本研究中，坡向、沉積物運移指數(shù)、海拔、沙蒿蓋度、小葉楊蓋度及硬質(zhì)早熟禾蓋度對生物結(jié)皮厚度的總的決定作用僅為55.2%，其可能是因為結(jié)皮厚度本身的變化就不大，影響生物結(jié)皮厚度的因素較為復(fù)雜，且主要是發(fā)育年限［20］，而本研究區(qū)域內(nèi)均為撂荒30a以上的樣地，故并沒有考慮這個因素。除此之外，本研究僅從地形因子（坡向、沉淀物運移指數(shù)及海拔）和植被因子（沙蒿蓋度、小葉楊蓋度及硬質(zhì)早熟禾蓋度）兩方面對結(jié)皮厚度進(jìn)行研究，沒有更多的考慮土壤因子。本研究僅從環(huán)境因子這一較為宏觀的方面出發(fā)，分析其對生物結(jié)皮的影響，以獲得其主要的影響因子，而對影響生物結(jié)皮的微觀因子及其影響機(jī)理尚待深入研究。

［1］ Greene R S B，Eldridge D J.Microbiotic soil crusts：A view of their roles in soil and ecological processes in the rangelands of Australia［J］.Australian Journal of Soil Research，1994，32（3）：389－415.

［2］ 李聰會，朱首軍，陳云明，等.黃土丘陵區(qū)生物結(jié)皮對土壤抗蝕性的影響［J］.水土保持研究，2013，20（3）：6－10.

［3］ 李新榮，張元明，趙允格.生物土壤結(jié)皮研究：進(jìn)展、前沿與展望［J］.地球科學(xué)進(jìn)展，2009，24（1）：11－24.

［4］ 熊好琴，段金躍，王妍，等.毛烏素沙地生物結(jié)皮對水分入滲和再分配的影響［J］.水土保持研究，2011，18（4）：82－87.

［5］ Fernando T M，Cristina E，Bowker M A.Biological crusts as a model system for examining the biodiversity ecosystem function relationship in soils［J］.Soil Biol.Biochem.，2009，42：1－13.

［6］ 吳明山，胥輝.通徑分析在林業(yè)多元線性模型共線性分析中的應(yīng)用研究［J］.西部林業(yè)科學(xué)，2008，37（3）：83－86.

［7］ 趙益新，陳巨東.通徑分析模型及其在生態(tài)因子決定程度研究中的應(yīng)用［J］.四川師范大學(xué)學(xué)報：自然科學(xué)版，2007，30（1）：120－123.

［8］ 卜崇峰，張朋，葉菁，等.陜北水蝕風(fēng)蝕交錯區(qū)小流域苔蘚結(jié)皮的空間特征及其影響因子［J］.自然資源學(xué)報，2014，29（3）：490－499.

［9］ 焦雯珺，朱清科，張宇清，等.陜北黃土區(qū)退耕還林地生物結(jié)皮分布及其影響因子研究［J］.北京林業(yè)大學(xué)學(xué)報，2007，29（1）：102－107.

［10］ 桑廣書.黃土高原歷史時期植被變化［J］.干旱區(qū)資源與環(huán)境，2005，19（4）：54－58.

［11］ 李國英.統(tǒng)一思想全面推進(jìn)黃河流域水土保持工作［J］.中國水土保持，2003（12）：9－10.

［12］ 盧曉杰，張克斌，李瑞.北方農(nóng)牧交錯帶生物結(jié)皮的主要影響因子探討［J］.水土保持研究，2007，14（6）：1－4.

［13］ 吳玉環(huán)，高謙，于興華.生物土壤結(jié)皮的分布影響因子及其監(jiān)測［J］.生態(tài)學(xué)雜志，2003，22（3）：38－42.

［14］ 呂建亮，廖超英，孫長忠，等.黃土地表藻類結(jié)皮分布影響因素研究［J］.西北林學(xué)院學(xué)報，2010，25（1）：11－14.

［15］ 成向榮，黃明斌，邵明安.神木水蝕風(fēng)蝕交錯帶主要人工植物細(xì)根垂直分布研究［J］.西北植物學(xué)報，2007，27（2）：321－327.

［16］ 李勉，李占斌，劉普靈，等.黃土高原水蝕風(fēng)蝕交錯帶土壤侵蝕坡向分異特征［J］.水土保持學(xué)報，2004，18（1）：63－65.

［17］ 賈恒義，雍紹萍，王富乾.神木地區(qū)的土壤資源［C］∥中國科學(xué)院西北水土保持研究所集刊.西安：陜西科學(xué)技術(shù)出社，1993.

［18］ 張琪，叢鵬，彭勵.通徑分析在Excel和SPSS中的實現(xiàn)［J］.農(nóng)業(yè)網(wǎng)絡(luò)信息，2007（3）：109－110.

［19］ 袁志發(fā)，周靜芋，郭滿才，等.決策系數(shù)—通徑分析中的決策指標(biāo)［J］.西北農(nóng)林科技大學(xué)學(xué)報：自然科學(xué)版，2001，29（5）：131－133.

［20］ 孟杰.黃土高原水蝕交錯區(qū)生物結(jié)皮的時空發(fā)育特征研究［D］.陜西楊凌：西北農(nóng)林科技大學(xué)，2011.