袁志良，陳 云，韋博良，張斌強，汪東亞，葉永忠，*

(1.河南農(nóng)業(yè)大學生命科學學院，鄭州 450002;2.三門峽市陜州公園管理處，三門峽 472000)

物種如何在有限空間內(nèi)共存一直是群落生態(tài)學研究的核心內(nèi)容［1-4］，所有物種共存理論中，生態(tài)位理論和群落中性理論目前最受關注［5-9］。生態(tài)位理論認為，不同的物種適應于不同的生境，并受不同生境因子的限制，物種通過占有不同資源和時間、空間來實現(xiàn)共存，所以在生態(tài)位主導的植物群落中，不同的物種會進化至適應各自特定的環(huán)境，從而表現(xiàn)出一定的生境相關性［10-14］。然而中性理論主導的群落物種與各類生境間則沒有明顯的相關性存在［15-16］。因此，物種生境關聯(lián)研究有助于加深對物種共存機制的理解和認識。

近年來，國內(nèi)外森林樣地中陸續(xù)開展了一系列有關物種生境關聯(lián)的研究工作，但多集中在單一氣候帶內(nèi)，如熱帶［17-19］、亞熱帶［20-21］和溫帶［22］，相鄰氣候帶間生態(tài)過渡區(qū)內(nèi)的此類研究則少有報道。而生態(tài)過渡區(qū)群落不僅具有相鄰區(qū)域的一些特性，而且具有獨特性和較高的豐富度，加之其對全球氣候變化反應的敏感性，使生態(tài)過渡區(qū)物種的研究在全球物種共存機制研究中占有非常重要的地位。

寶天曼國家級自然保護區(qū)地處典型的暖溫帶-北亞熱帶過渡區(qū)，以及我國第二級地貌階梯向第三級地貌階梯過渡的邊緣，其植被屬暖溫帶落葉闊葉林向亞熱帶常綠闊葉林的過渡型［23］，復雜多樣的氣候、地形條件促進了不同植物區(qū)系的交匯和融合，使得本區(qū)成為我國同緯度生態(tài)結構保存最完整的地區(qū)，同時也是我國十大生物多樣性分布區(qū)之一。對其物種共存機制的研究將有助于進一步理解生態(tài)過渡區(qū)物種分布和共存特點，對實施物種多樣性保護具有重要意義。

許多研究表明地形對物種形成不同的格局有著重要的作用［24-26］，然而土壤是作用于植物的直接因素，地形是間接因素［21］，因此本研究選取地形與土壤兩類環(huán)境因子，研究了生態(tài)過渡區(qū)內(nèi)主要物種生境相關性，探討(1)物種生境特化在暖溫帶-北亞熱帶過渡區(qū)森林生態(tài)系統(tǒng)生物多樣性維持的作用;(2)對比地形因子和土壤因子對物種空間分布影響;(3)樣地中樹種與各類生境的相關性怎樣，哪類生境在樣地中發(fā)揮更重要的作用?

1 材料及研究方法

1.1 研究區(qū)概況

寶天曼國家級自然保護區(qū)位于河南省西南部，地處秦嶺東段，伏牛山南坡，地理坐標為111°46'55″—112°03'32″E，33°35'43″—33°20'12″N，南北長24.3 km，東西寬25.9 km，總面積23198 hm2。山體呈東南西北走向，最高峰寶天曼海拔1830 m。年均氣溫15.1℃，年均降水量885.6 mm，年均蒸發(fā)量99.6 mm。土壤可劃分為山地棕壤、山地黃棕壤和山地褐土3種，土壤pH值為6.5左右。

1.2 樣地調(diào)查

2008年至2009年的6—10月間，在保護區(qū)的核心區(qū)內(nèi)選取群落保護較好的區(qū)域，參照BCI 50 hm2樣地(CTFS)的技術規(guī)范［27］，建立了本樣地(100 m×100 m)(圖1)。用全站儀將樣地分成25個20 m×20 m的樣方，每個樣方又進一步分為16個5 m×5 m的小樣方。測量樣方內(nèi)所有胸徑≥1 cm的木本植物的胸徑和坐標，定位并掛牌。測量每個5 m×5 m的小樣方的4個頂點的海拔高度，以計算樣地微地形。

參照CTFS的土壤采樣方案進行土壤采樣［28］:以20 m×20 m樣方的4個頂點為基準點，在此基準點的8個方向(東、南、西、北、東南、東北、西南、西北)中隨機選取一個方向，并在此方向隨機選取一個坐標延伸5 m或10 m的樣點，共有72個土壤采集點，樣品帶回實驗室，參照《土壤農(nóng)化分析》［29］測定其全碳、全氮、有效鎂、有效鉀、有效鐵、有效錳、有效硫、有效鉬、有效磷、有效氮、有效鋅、有效銅、有效鈣、pH、有效鋁 15個土壤因子。

1.3 數(shù)據(jù)分析

選擇樣地內(nèi)胸徑≥1 cm、個體數(shù)≥5的43種木本植物作為研究對象。選擇4個地形因子和15個土壤因子作為兩類環(huán)境因子。用典范對應分析CCA分別研究兩類環(huán)境因子對樹種分布的影響差異，并用隨機化排列檢驗［30］來分析兩類環(huán)境因子對樹種分布影響的顯著性;用Torus-translation檢驗分別研究樹種對兩類環(huán)境下不同生境的選擇差異。數(shù)據(jù)分析采用R2.15.3軟件實現(xiàn)。

圖1 樣地等高線圖Fig.1 Contour map of the sample plot

1.3.1 環(huán)境因子對樹種分布的影響

CCA分析中的物種矩陣由每個樹種在10 m×10 m樣方中的多度組成。通過樣點的土壤因子對整個樣方的土壤因子用克里格［31］進行插值估計，得到每個10 m×10 m樣方土壤因子的估計值，作為土壤因子矩陣。用調(diào)查得到的海拔數(shù)據(jù)計算每個10 m×10 m樣方的平均海拔、凹凸度、坡度和坡向4個生境指標作為地形因子矩陣。

樣方平均海拔是樣方4個頂點海拔的平均值［17，32］;從樣方4 個頂點取 3 個頂點組成一個平面，4個頂點可以組合成4個不同的平面，這4個平面與樣方投影面夾角的平均值和與正北方向角度的平均值分別為樣方的坡度與坡向值［33］;樣方的凹凸度是樣方的平均海拔減去與該樣方相鄰的8個樣方平均海拔的平均值，處于樣地邊緣樣方的凹凸度為樣方中心的海拔減去4個頂點海拔的平均值［34］。

1.3.2 生境分類

在物種—地形生境相關性分析過程中，選取4個地形因子，結合地形圖的識別，利用多元回歸樹［35］作為地形生境劃分的輔助方法，把寶天曼1 hm2樣地分為溝谷，緩坡，陡坡和山脊4種生境類型(圖2)，其所占的樣方數(shù)分別為28、46、13和13個，總面積分別為2800、4600、1300 m2和 1300 m2。

圖2 寶天曼1 hm2樣地100個樣方4類地形生境分類后圖示Fig.2 The100 quadratsof the1-hm2 plot in Baotianman divided into four habitat types

15個土壤因子之間存在很大相關性，為了壓縮和提取出土壤因子的主要的信息，本研究對土壤因子進行主成分分析，從而獲得代表土壤因子大部分信息的主分量，并將其代表土壤生境狀況。本文提取的前兩個主分量代表99%土壤因子的信息，其中第一主分量(PC1)代表總體信息的92%，第二主分量(PC2)代表7%。因此，第一主分量為最重要的主分量，這一主分量與有效鋁、有效鐵具有正相關關系，與有效鈣、有效鎂有負相關關系。第二主分量與有效銅、有效錳、有效硫、全碳、全氮、pH具有正相關關系。兩個主分量的空間異質性如圖3所示。兩主分量空間變化趨勢明顯，結合此圖，以兩個主分量的中值劃分生境的方法［19，36］，把兩主分量分別劃為高濃度區(qū)，中濃度區(qū)和低濃度區(qū)3個生境，具體生境參數(shù)如表1所示。

圖3 寶天曼1 hm2樣地土壤因子2個主分量空間分布圖Fig.3 Maps of 2 PCs at Baotianman plot

表1 寶天曼1 hm2樣地中幾類土壤生境參數(shù)統(tǒng)計Table 1 The physical parameters of several soil habitat categories in the 1 hm2 Baotianman plot

1.3.3 樹種與兩類環(huán)境下不同生境的關聯(lián)分析

Torus-translation檢驗是目前最常用的檢驗物種與生境關聯(lián)的方法［12，19，36-38］。和常規(guī)的方法不同，這種方法的優(yōu)點在于能夠較大程度上排除生境和物種分布的空間自相關，從而使檢驗更加敏感［17］。其基本思想是計算出一個物種在每個生境中的真實分布在隨機分布條件下的概率，通過概率確定一個物種是否與某一類生境顯著相關。

2 結果分析

2.1 兩類環(huán)境因子對樹種分布的影響

以CCA的結果進行Monte Carlo檢驗，分析4個地形因子和15個土壤因子對物種分布的解釋量，顯著性都小于1，說明排序的結果可以接受環(huán)境因子對物種分布的解釋量。用隨機化排列檢驗對兩類環(huán)境因子前2個排序軸進行檢驗，排序軸較好的體現(xiàn)了生境梯度。

物種和兩類環(huán)境不同生境排序結果及特征值分別見圖4和表2，第1軸(CCA1)為水平軸，第2軸(CCA2)為垂直軸。地形因子對樹種分布的解釋量為7.3%，4個軸依次解釋了13.1%、4.3%、2.8%和2.3%的變異?？山忉尩臉浞N分布與土壤因子的關系為16.1%，第一軸解釋了14.1%，第二軸解釋了10.3%。物種在排序軸上的分布基本反映了群落空間分布特征隨環(huán)境梯度的變化趨勢。

2.2 樹種與兩類環(huán)境下不同生境的關聯(lián)

基于Torus-translation檢驗表明，在被檢驗的43個物種中，23.3%的物種與地形4個生境關聯(lián)，46.5%的物種與土壤6個生境關聯(lián)。與兩類環(huán)境因子下不同生境關聯(lián)的物種數(shù)占總物種數(shù)的比例如圖5所示。

正相關關系中，與高濃度區(qū)、中濃度區(qū)和低濃度區(qū)正相關物種數(shù)最多的是兩個主分量中的高濃度區(qū)，占正相關數(shù)的52.6%，明顯高于中濃度區(qū)和低濃度區(qū)正相關物種數(shù)，而在地形生境中，正相關物種數(shù)最多的是山脊，占正相關數(shù)的57%。負相關關系中，與土壤生境負相關物種數(shù)最多的仍是高濃度區(qū)，占負相關數(shù)的42.8%，其次為低濃度區(qū)，占26.3%，在地形4類生境中表現(xiàn)負相關的物種很少，只有4種。

表2 物種CCA排序結果及其特征值Table2 Results of CCA and eigenvalues of species

圖4 43種木本植物和地形因子與土壤因子CCA排序圖Fig.4 CCA biplot of the 43 species and topographic factors and soil factors

3 討論

3.1 土壤生境分化對暖溫帶-北亞熱帶過渡區(qū)物種空間分布有重要作用

通過對樣地土壤因子進行主成分分析，提取出來的前兩個主分量不能作為一個嚴格的土壤生境分類方法［21］，但是總體上反映了不同土壤因子特征。Torus-translation檢驗表明，有30.2%的物種空間分布與第一主分量顯著相關，說明第一主分量中有效鋁、有效鐵、有效鉬、全碳、全氮、有效磷、有效鉀、有效氮的空間異質性對樣地物種空間分布影響明顯。其他微量元素的空間異質性主要表現(xiàn)在第二主分量，雖然只代表了所有土壤因子信息的7.6%，卻對物種空間分布的作用不小，有30.2%的物種與其顯著相關。由此可以看出，微量元素雖然含量少，但其空間異質性對于物種空間分布具有重要的作用［39］。

圖5 43種木本植物與地形因子和土壤因子不同生境相關聯(lián)Fig.5 43 species of woody plants are associated in different habitats with topographic factors and soil factors

共有46.5%的物種空間分布與兩個主分量顯著正相關或顯著負相關，所以樣地土壤生境分化影響著物種的空間分布。John等分析了土壤因子空間異質性對熱帶森林3個大樣地物種空間分布的影響［28］(50 hm2BCI樣地、25 hm2La Planada樣地和25 hm2Yasuni樣地)，分別有40%、36%、51% 的物種空間分布和土壤因子存在顯著相關，證實了土壤養(yǎng)分在熱帶雨林中對物種空間分布起到重要的作用。張儷文分析了土壤因子空間異質性對亞熱帶古田山大樣地物種空間分布的影響，發(fā)現(xiàn)有88.2%的常見種的空間分布表現(xiàn)出與四個主分量有相關性，得出古田山樣地土壤因子的空間異質性影響著物種的空間分布［21］。寶天曼樣地內(nèi)近一半的物種空間分布與土壤生境表現(xiàn)出顯著正相關或顯著負相關，所以有理由推斷土壤生境分化對暖溫帶-北亞熱帶過渡區(qū)物種的空間分布也有著重要作用。

3.2 土壤因子與地形因子對物種分布影響比較

CCA分析顯示4種地形因子對物種空間分布的解釋量為7.3%，15種土壤因子解釋量為16.1%。通過Torus-translation方法與生境的相關性進行檢驗，23.3%的物種至少與地形中一類生境類型顯著相關，46.5%的物種至少與土壤中一類生境類型顯著相關。從CCA分析和Torus-translation檢驗結果可以發(fā)現(xiàn)，土壤因子對物種分布的影響要遠遠大于地形因子對物種分布的影響。這可能是因為樣地面積相對于其它樣地較小，地形變化不大，沒能顯現(xiàn)出地形對所有物種分布的影響。然而，很多物種與地形生境關聯(lián)性研究都表明地形作為環(huán)境因子只能解釋部分原因［20，40］。土壤的空間異質性和地形是密切相關的，地形造成了土壤中可獲取的水分和養(yǎng)分的梯度［41-43］，但土壤因子是作用于植物的直接因素，地形是間接因素。

通過對比4類地形生境劃分圖(圖2)與土壤因子2個主分量空間分布圖(圖3)，可發(fā)現(xiàn)第一主分量高濃度區(qū)主要分布在低海拔溝谷與緩坡生境，低濃度區(qū)主要分布在高海拔山脊和陡坡。也就是說溝谷、緩坡養(yǎng)分含量相對于山脊、陡坡更豐富，這與Clark等研究結果相似［44-45］;而第二主分量高濃度區(qū)主要分布在高海拔陡坡，低濃度區(qū)主要分布在低海拔溝谷和緩坡，這是因為第二主分量所包含的元素主要都屬于微量元素，含量都特別低。15個土壤因子不僅反應了樣地地形部分信息，還反應了很多土壤因素信息。本文從更直接的角度證實了寶天曼1 hm2樣地內(nèi)大部分物種存在生境分化，這不僅是對地形生境分化的利用，而且更多的是對土壤生境利用的分化，這與古田山張儷文研究結果一致［21］。

3.3 物種生境選擇差異與物種共存

由Torus-translation檢驗表明，地形和土壤生境分化是本區(qū)物種共存的主要機制之一。在4類地形生境表現(xiàn)顯著正相關的物種中，85.7%的物種與山脊和陡坡正相關，而對溝谷和緩坡也沒有明顯的排斥作用。因為樣地處于北坡，接受的陽光相對較少，相對于氮、磷、鉀等資源，光照成為樣地的限制因子，這可能是大部分物種會選擇在山脊和陡坡生境中生長的原因。

在土壤生境中，表現(xiàn)顯著正相關或負相關的物種主要集中在兩個主分量的高濃度區(qū)或低濃度區(qū)，與兩個高濃度區(qū)顯著正相關的物種占到所有與土壤生境相關物種數(shù)的50%，與低濃度區(qū)負相關的物種數(shù)占所有相關物種數(shù)的25%;而與高濃度區(qū)表現(xiàn)負相關的物種占所有相關物種數(shù)的30%，與低濃度區(qū)顯著正相關的的物種占所有相關物種數(shù)的25%。這說明物種對生境的選擇差異主要集中在土壤高濃度與低濃度區(qū)。其中與高濃度區(qū)表現(xiàn)負相關的物種有:鵝耳櫪(Carpinus turczaninowii)、銳齒櫟(Quercus aliena var.acutiserrata)、三椏烏藥(Lindera obtusiloba var.obtusiloba)、白檀(Symplocos paniculata)、樺葉莢蒾(Viburnum betulifolium)、四照花(Dendrobenthamia japonica var.chinensis)，與低濃度區(qū)表現(xiàn)正相關的物種有:柳葉衛(wèi)矛(Euonymus cornutus Hemsl)、白檀(Symplocos paniculata)、大葉樸(Celtis koraiersi)、樺葉莢蒾(Viburnum betulifolium)、水榆花楸(Sorbus alnifolia)。高濃度區(qū)有著豐富的植物生長所需的營養(yǎng)元素，如全碳、全氮、有效磷、有效鉀、有效氮、有效鉬等，這在一定程度上能夠解釋為什么樹種會更偏好在高濃度區(qū)，這與本文的假設相符;但有一部分物種傾向于選擇在低濃度區(qū)，并且主要是與第二主分量低濃度正相關，這部分物種在樣地內(nèi)大部分屬于優(yōu)勢種，耐受性強，生態(tài)位寬度大，并且可能與有效銅、有效錳、有效硫等微量元素以及pH有關。

4 結論

本文對地形因子和土壤因子在暖溫帶-北亞熱帶過渡區(qū)物種多樣性的維持機制進行研究，結果支持物種共存機制中的生態(tài)位理論，地形和土壤生境分化是寶天曼暖溫帶-北亞熱帶過渡區(qū)物種共存的重要機制之一，在一定程度上解釋了該區(qū)物種多樣性的維持。研究結果將有助于進一步理解生態(tài)過渡區(qū)物種分布和共存特點，對實施物種多樣性保護具有重要意義。

然而，CCA分析表明了土壤因子和地形因子只能解釋物種分布和共存的部分原因，其它的環(huán)境因素，如水分、光照等對物種分布和共存的影響都有待進一步研究。此外，本文主要從物種整體水平上研究了與各類生境關聯(lián)情況，如果考慮到物種本身的特性對生境選擇的差異影響，生態(tài)位理論可能解釋更多的物種共存機制。

［1］Bengtsson J，F(xiàn)agerstrom T，Rydin H.Competition and Coexistence in Plant-Communities.Trends in Ecology＆ Evolution，1994，9:246-250.

［2］Nakashizuka T.Species coexistence in temperate，mixed deciduous forests.Trends in Ecology＆ Evolution，2001，16:205-210.

［3］Wright SJ.Plant diversity in tropical forests:a review of mechanisms of species coexistence.Oecologia，2002，130:1-14.

［4］Wilson W G，Lundberg P，Vazquez DP，Shurin JB，Smith MD，Langford W，Gross K L，Mittelbach GG.Biodiversity and species interactions:extending Lotka-Volterra community theory.Ecology Letters，2003，6:944-952.

［5］Connell J H.Diversity in tropical rain forests and coral reefs.Science，1978，199:1302-1310.

［6］Leigh E G.Tropical Forest Ecology.Oxford:Oxford University Press，1999:1-260.

［7］Hubbell，Stephen P.The unified neutral theory of biodiversity and biogeography(MPB-32).Vol.32.Princeton University Press，2001:1-390.

［8］Whitfield J.Ecology:neutrality versus the niche.Nature，2002，417:480-481.

［9］Potts M D，Davies S J，Bossert W H，Tan S，Supardi M N N.Habitat heterogeneity and niche structure of trees in two tropical rain forests.Oecologia，2004，139:446-453.

［10］Hubbell SP，F(xiàn)oster R B.Diversity of canopy trees in a neotropical forest and implications for conservation.In:Tropical Rain Forest:Ecology and Management(eds Sutton SL，Whitmore TC，Chadwick AC).Oxford:Blackwell Scientific Publications，1983:25-41.

［11］Tilman D，Pacala S.The maintenance of species richness in plant communities//Ricklefs R E，Schluter D，eds.Species Diversity in Ecological Communities.Chicago:University of Chicago Press，1993:13-25.

［12］Debski I，Burslem D F R P，Palmiotto P A，Lafrankie J V，Lee H S，Manokaran N.Habitat preferences of Aporosa in two Malaysian forests:Implications for abundance and coexistence.Ecology，2002，83:2005-2018.

［13］Potts M D，Ashton P S，Kaufman L S，Plotkin J B.Habitat patterns in tropical rain forests:a comparison of 105 plots in northwest Borneo.Ecology，2002，83:2782-2797.

［14］Itoh A，Yamakura T，Ohkubo T，Kanzakim M，Palmiotto P A，Lafrankie JV，Ashton PS，Lee H S.Importance of topography and soil texture in the spatial distribution of two sympatric dipterocarp trees in a Bornean rain forest.Ecological Research，2003，18:307-320.

［15］Hubbell SP，F(xiàn)oster R B.Biology，chance and history and the structure of tropical rain forest tree communities.In:Community Ecology.Harper＆Row，New York，1986:314-319.

［16］Hubbell S P，F(xiàn)oster R B，O'Brien S T，Harms K E，Condit R，Wechsler B，Wright S J，de Lao S L.Light-gap disturbances，recruitment limitation，and tree diversity in a neotropical forest.Science，1999，283:554-557.

［17］Harms K E，Condit R，Hubbell SP，F(xiàn)oster R B.Habitat associations of trees and shrubs in a 50-hm2neotropical forest plot.Journal of Ecology，2001，89:947-959.

［18］Svenning JC.Microhabitat specialization in a species-rich palm community in Amazonian Ecuador.Journal of Ecology，1999，87:55-65.

［19］Gunatilleke C V S，Gunatilleke IA U N，Esufali S，Harms K E，Ashton P M S，Burslem D FR P，Ashton P S.Species-habitat associations in a Sri Lankan dipterocarp forest.Journal of Tropical Ecology，2006，22:371-384.

［20］Lai JS，Mi X C，Ren H B，Ma K P.Species-habitat associations change in a subtropical forest of China.Journal of Vegetation Science，2009，20:415-423.

［21］Zhang L W.The Effect of Spatial Heterogeneity of Environmental Factors on Species Distribution and Community Structure［D］.Beijing:Institute of Botany，Chinese Academy of Sciences，2010.

［22］Liu H F，Sang W G，Xue D Y.Topographical habitat variability of dominant species populations in a warm temperate forest.Chinese Journal of Ecology，2013，32(4):795-801.

［23］Shang F D，Wang P J，F(xiàn)eng G P，Zeng Y，Wang Z D.Study on the Characteristics and causes of formation of plant diversity in the Funiushan Transition Region.Journal of Henan University(Natural Science)，1998，28(1):54-60.

［24］Ohsawa M.Habitat differentiation and the ecological niche in vegetation.Physiological Ecology，Japan，1987，24:s15-s27.

［25］Kikuchi T，Miura O.Vegetation patterns in relation to micro-scale landforms in hilly land regions.Vegetatio，1993，106:147-154.

［26］Sakai A，Ohsawa M.Topographical pattern of the forest vegetation on a river basin in a warm-temperate hilly region，central Japan.Ecological Research，1994，9:269-280.

［27］Condit R.Tropical forest census plots:methods and results from Barro Colorado Island，Panama and a comparison with other plots.Springer，1998:1-224.

［28］John R，Dalling J W，Harms K E，Yavitt J B，Stallard R F，Mirabello M，Hubbell SP，Valencia R，Navarrete H，Vallejo M，F(xiàn)oster R B.Soil nutrients influence spatial distributions of tropical tree species.Proceedings of the National Academy of Sciences of the United States of America，2007，104:864-869.

［29］Bao SD.Soil and agricultural chemistry analysis.Agriculture Publication，Beijing，2000:355-356.

［30］Wan Y，Cohen J，Guerra R.A permutation test for the robust sib-pair linkage method.Annals of Human Genetics，2005，61:79-87.

［31］Su S，Li A W，Liu Q H.Ordinary Kriging Algorithm Applied for Interpolation.Journal of Southern Yangtze University(Natural Science Edition)，2004，3(1):19-21.

［32］Valencia R，F(xiàn)oster R，Villa G，Condit R，Svenning JC，Hernandez C，Romoleroux K，Losos E，Magard E，Balslev H.Tree species distributions and local habitat variation in the Amazon:a large plot in eastern Ecuador.Journal of Ecology，2004，92:214-229.

［33］Nee S.The neutral theory of biodiversity:do the numbers add up?.Functional Ecology，2005，19:173-176.

［34］Yamakura T，Kanzake M，Itoh A，Ohkubo T，Ogino K，Chai E OK，Lee H S，Ashton P S.Topography of a large-scale research plot established within a tropical rain forest at Lambir，Sarawak.Tropics，1995，5:41-56.

［35］De'ath G.Multivariate regression trees:a new technique for modeling species-environment relationships.Ecology，2002，83:1105-1117.

［36］Yamada T，Tomita A，Itoh A，Yamakura T，Ohkubo T，Kanzaki M，Tan S，Ashton P S.Habitat associations of Sterculiaceae trees in a Bornean rain forest plot.Journal of Vegetation Science，2006，17:559-566.

［37］DeWalt SJ，Ickes K，Nilus R，Harms K E，Burslem D F R P.Liana habitat associations and community structure in a Bornean lowland tropical forest.Plant Ecology，2006，186:203-216.

［38］Comita L S，Condit R，Hubbell SP.Developmental changes in habitat associations of tropical trees.Journal of Ecology，2007，95:482-492.

［39］Sollins.Factors inflluencing species composition in tropical lowland rain forest:does soil matter?Ecology，1999，79(1):23-30.

［40］Wang W，Luo Z R，Zhou R F，Xu D M，Ai J G，Ding B Y.Habitat associations of woody plant species in Baishanzu subtropical broad-leaved evergreen forest.Biodiversity Science，2011，19(2):134-142.

［41］Zak D R，Hairston A，Grigal D F.Topographic influences on nitrogen cycling within an upland pin oak ecosystem.Forest Science，1991，37:45-53.

［42］Enoki T，Kawaguchi H，Iwatsubo G.Nutrient-uptake and nutrient-use efficiency of Pinus thunbergii Parl，along a topographical gradient of soil nutrient availability.Ecological Research，1997，12:191-199.

［43］Hirobe M，Tokuchi N，Iwatsubo G.Spatial variability of soil nitrogen transformation patter ns along a forest slope in a Cryptomeriajaponica D.Don plantation.European Journal of Soil Biology，1998，34:123-131.

［44］Clark D B，Clark D A，Read JM.Edaphic variation and the mesoscale distribution of tree species in a neotropical rain forest.Journal of Ecology，1998，86:101-112.

［45］Clark D B，Palmer M W，Clark D A.Edaphic factors and the landscape scale distributions of tropical rain forest trees.Ecology，1999，80:2662-2675.

參考文獻:

［21］張儷文.環(huán)境空間異質性對物種空間分布和群落結構的影響［D］.北京:中國科學院植物研究所，2010.

［22］劉海豐，桑衛(wèi)國，薛達元.暖溫帶森林優(yōu)勢種群的地形生境變異性.生態(tài)學雜志，2013，32(4):795-801.

［23］尚富德，王磐基，馮廣平，曾穎，王正德.伏牛山南北自然過渡地帶植物多樣性的特征及其成因分析.河南大學學報，1998，28(1):54-60.

［29］鮑士旦.土壤農(nóng)化分析，第三版.北京:中國農(nóng)業(yè)出版社，2005:355-356.

［31］蘇姝.普通克里格法在空間內(nèi)插中的應用.江南大學學報，2004，3(1):19-21.

［40］王偉，駱爭榮，周榮飛，許大明，哀建國，丁炳揚.百山祖常綠闊葉林木本植物的生境相關性分析.生物多樣性，2011，19(2):134-142.