中國典型森林生態系統喬木層群落物種多樣性的空間分布格局及其影響因素

吳安馳,鄧湘雯,,任小麗,項文化,,張 黎, ,葛 蓉, ,牛忠恩, ,何洪林,,*,何立杰

1 中南林業科技大學生命科學與技術學院,長沙 410004 2 南方林業生態應用技術國家工程實驗室,長沙 410004 3 中國科學院地理科學與資源研究所,生態系統網絡觀測與模擬重點實驗室, 北京 100101 4 中國科學院大學,北京 100049 5 沈陽農業大學,沈陽 110866

生物多樣性是指一個地區內生命形成的豐富程度,是生物以及與環境有規律地結合所構成穩定的生態綜合體以及與此相關的各種生態過程的綜合。它包括物種多樣性、遺傳多樣性及生態系統多樣性3個層次。其中,物種多樣性是生物多樣性的核心[1]。森林是陸地生態系統生物多樣性最重要的資源庫,但森林砍伐、森林退化、氣候變化等因素正威脅全球大約一半的樹種[2]。植物群落物種多樣性的空間分布格局是各種生態梯度變化的綜合反映,研究物種多樣性沿環境梯度的變化格局以及影響因子是生物多樣性研究的重要議題之一[3- 4]。

國內外學者圍繞森林物種多樣性的分布格局進行了廣泛、深入的研究[5- 9],全球和區域范圍內低海拔地區不同生態系統的研究結果表明,植物群落物種多樣性具有明顯的緯度和海拔梯度變化規律[10]。其中,針對我國山地植物物種多樣性垂直分布格局研究結果表明,在人為干擾較少的山地,木本植物物種豐富度一般隨著海拔的升高而逐漸減少[11- 13];但王志恒等[14]在高黎貢山的研究發現,全部種子植物和不同生活型植物的物種豐富度隨海拔升高呈現先增后減的趨勢,即在中海拔地段的豐富度最大。在大尺度上,森林群落物種豐富度一般隨緯度增加而減少[15]。馮建孟等[16]對中國種子植物物種多樣性大尺度分布格局進行研究,發現從隨緯度增加,物種密度呈顯著遞減趨勢,而物種豐富度的遞減趨勢不夠明顯;陳圣賓等[17]對中國種子植物特有屬地理分布格局進行了研究,發現種子植物屬豐富度及特有屬比例具有很強的空間變異性,且表現出顯著的緯度梯度性。目前,關于影響物種多樣性分布格局的因素爭議較大,在較小的空間尺度上,決定物種多樣性分布格局的主要因素是生物、生態過程(如競爭、共生和遷徙等)、地形和土壤因子等[18- 20]。在大尺度上,很多學者提出了若干理論或假說解釋物種多樣性的空間格局,如水熱動態假說[21- 22]、環境能量假說[23- 24]、生理忍耐假說[25- 26]、氣候穩定性假說[27- 28]和生產力假說[29- 31]等。雖然以往的研究對各種假說進行了大量的驗證,但不同區域物種多樣性的主導因素及各種假說的適用性存在很大差異。因此,物種多樣性分布格局及其影響因素的研究有待深入。

中國疆域遼闊,地形復雜,有豐富的木本植物以及多樣的環境梯度,從溫暖潮濕的熱帶到極度干旱的沙漠和寒冷的北方地區[32],適合開展大尺度物種多樣性格局的研究。衡量不同區域物種多樣性分布格局影響因素的相對貢獻率成為現階段物種多樣性格局研究的重點[33]。同時,使用不同區域、不同生態系統的實測數據是研究物種多樣性地理格局的形成機制以及假說驗證重要手段。本文基于國家生態系統觀測研究網絡(CNERN)的13個典型森林生態系統綜合觀測場樣地喬木層調查數據,利用物種豐富度指數、多樣性指數和均勻度指數3類指標來測度森林群落喬木層的物種多樣性,并結合氣候因子和植物特性數據進行統計分析,旨在揭示:(1)中國典型森林生態系統喬木層的物種多樣性隨經緯度變化的空間分布格局;(2)確定植物特性和氣候因子對喬木層的物種多樣性分布格局的相對貢獻率,解釋中國典型森林生態系統喬木層的物種多樣性格局形成的主導因子。

1 研究區域和研究方法

1.1 研究區概況

本研究選擇國家生態系統研究網絡(http://www. cnern. org. cn)13個典型森林生態系統綜合觀測場樣地,涵蓋熱帶(尖峰嶺和西雙版納)、亞熱帶(鼎湖山、哀牢山、會同、神農架和貢嘎山)和溫帶(秦嶺、茂縣、北京、長白山、帽兒山和大興安嶺)(圖1),基本覆蓋我國從南到北各種氣候帶和典型森林植被類型。各個綜合觀測場樣地經過長期的封山育林,人為干擾少,植被保護良好,對每個地區植被類型具有典型代表性。

圖1 中國13處典型森林生態系統綜合性樣地分布地點及物種多度直方圖Fig.1 Locations and species abundance histograms of the comprehensive sites of 13 typical forest ecosystem in China插圖中y坐標為物種數量,x坐標以2為底的物種多度進行對數轉換,即物種個體數分別為1,2,3—4,5—8,9—16等進行分組

1.2 野外觀測場設置和布局

13個森林生態站的綜合觀測場樣地均按照統一規范和方法設置,根據物種面積曲線確定最小面積樣方的大小,且位于所在地區內最具有代表性的生態系統的典型地段[34]。對每各站綜合觀測場樣地內喬木樹種的所有個體進行每木檢尺,記錄喬木層的物種、株數、胸徑、樹高、生長狀況和空間相對位置等指標,用卷尺記錄物種個體空間坐標精確到±15 cm。各個森林生態系統的生物量數據通過相對生長方程計算法(allometric method)計算獲得,即通過胸徑、樹高與各器官地上、地下生物量構建經驗方程估算森林生態系統的生物量。同時,記錄樣地的經緯度、海拔和郁閉度等,樣地的基本信息[35]見表1。由于在較大空間格局上,喬木層的物種多樣性與海拔沒有顯著的相關性(P>0.05),本文的研究未涉及海拔對喬木層物種多樣性的影響。

表1 中國13處典型森林生態系統綜合觀測場樣地的基本信息

九節Psychotriarubra; 大葉白顏Gironnierasubaequalis; 番龍眼Pometiapinnata; 千果欖仁Terminaliamyriocarpa; 錐CastanopsischinensisHance; 木荷Schimasuperba; 厚殼桂Cryptocaryachinensis; 硬殼柯Lithocarpushancei; 木果柯Lithocarpusxylocarpus; 變色錐Castanopsisrufescens; 紅拷Castanopsishystrix; 青岡Cyclobalanopsisglauca; 刨花楠Machiluspauhoi; 米心青岡Fagusengleriana; 多脈青岡Cyclobalanopsismultinervis; 曼青岡Cyclobalanopsisoxyodon; 冷杉Abiesfabri; 五尖槭AcermaximowicziiPax; 糙皮樺Betulautilis; 華山松Pinusarmandii; 油松Pinustabuliformis; 銳齒槲櫟Quercusaliena; 銳齒礫Quercusalienavar.acuteserrata; 遼東櫟Quercuswutaishanica; 白樺Betulaplatyphylla; 紅松Pinuskoraiensis; 色木槭Acermono; 假色槭Acerpseudo-sieboldianum; 水曲柳Fraxinusmandschurica; 興安落葉松Larixgmelinii

1.3 物種多樣性測度方法

根據中國典型森林生態系統的綜合觀測場樣地喬木層調查數據,本研究選用物種豐富度(S)、Shannon-Wiener指數(H′)、Simpson指數(D)和 Pielou指數(E)4個常用指數測度群落喬木層的物種多樣性[36],具體計算公式如下:

物種豐富度(S):

S=出現在樣地內的物種數

(1)

Simpson指數(D):

(2)

Shannon-Wiener指數(H′):

H′=-∑PiInPi

(3)

Pielou指數(E):

E=H′/InS

(4)

式中,S為樣方中物種數,Pi為群落中第i個物種的重要值。物種重要值(important value,IV)用來反映物種在群落中的相對重要性,計算公式為:

重要值(IV)=(相對多度+相對優勢度+相對高度)/3

(5)

1.4 氣候數據和NDVI指數

根據各森林生態站實際觀測的月均溫和月降水數據(2000—2010年平均值),計算年均溫(annual mean temperature, MAT)、最暖月均溫(mean temperature of the warmest month, MTWM)、最冷月均溫(mean temperature of the coldest month, MTCM)、溫度年較差(annual variation of mean annual temperature, AVMAT)、年降水量(annual mean precipitation, MAP)、最干燥季降水量(precipitation of driest quarter, PDQ)、最濕潤季降水量(precipitation of wettest quarter, PWQ)、實際蒸散量(actual evapotranspiration, AET)和潛在蒸散量(potential evapotranspiration, PET),其中潛在蒸散量和實際蒸散量根據魏煥奇等[37]方法計算。這些指標不僅與物種多樣性有著密切關系,也是驗證各假說的關鍵指標[38]。

歸一化植被指數(normalized difference vegetation index, NDVI)可以很好地反映植被生長狀況,與植被覆蓋度和生產力顯著相關[39]。本研究使用的NDVI數據來自美國國家航空航天局NASA網站(http://glovis. usgs. gov)下載的MOD13Q1數據產品,其空間分辨率為250 m,根據綜合觀測場樣地的地理坐標在Arcgis 10.0軟件中提取2000—2010年均值用于數據分析。

1.5 數據處理與分析

運用一元回歸模型檢驗物種多樣性指數與經緯度的關系,并對物種多樣性指數與各自變量進行單因素方差分析(one-way ANOVA),分別比較物種多樣性指數和各自變量之間Pearson相關性系數,所用的統計顯著水平值設為P<0.05。

將所有因子分為植物特性、能量和水分三大類。其中,植物特性包括NDVI指數、郁閉度和生物量;能量因子包括年均溫、最暖月均溫、最冷月均溫、溫度年較差、實際蒸散量和潛在蒸散量;水分因子包括年降水量、最干燥季降水量和最濕潤季降水量。對三類因子分別進行主成分分析(principal component analysis, PCA),提取主成分作為各類因子的替代因子。

數據處理和分析采用Excel 2010、SPSS 18.0和R 3.3.3(R-Development Core Team 2017)軟件完成。物種豐富度在分析之前經過自然對數轉換。方差分解(variation partitioning)分析采用R語言的“Vegan”數據包中“Varpart”功能完成。運用韋恩圖顯示各類因子的純效應(pure effect)和共同效應(share effect)[40- 41]。

物種的個體數轉換成以2為底的對數,并進行合并,即按照物種的個體數分別為1, 2, 3—4, 5—8, 9—16等進行分組,按照這種方法作物種多度直方圖[42]。表明中國典型森林生態系統綜合觀測場樣地喬木層的常見物種和稀有種差異(圖1)。統計分析在R軟件包“sads”中完成。

2 結果與分析

2.1 物種多樣性指數隨經度和緯度的分布格局

中國典型森林生態系統的喬木層物種多樣性指數存在較大的空間變異(圖2)。其中,物種豐富度與緯度(R2=0.765,P<0.001)和經度(R2=0.410,P<0.05)均呈顯著負相關,即隨經緯度的升高物種豐富度呈逐漸下降趨勢。而Shannon-Wiener指數(R2= 0.694,P<0.001)、Simpson指數(R2=0.592,P<0.01)和Pielou指數(R2=0.315,P<0.05)僅與緯度呈顯著負相關,與經度無顯著的負相關關系(P>0.05)。物種豐富度、Shannon-Wiener指數和Simpson指數最高值位于熱帶的西雙版納,最低值位于寒溫帶的大興安嶺。Pielou指數波動較大,例如溫帶長白山的Pielou指數值(0.86)高于亞熱帶哀牢山的Pielou指數值(0.81)。雖然Shannon-Wiener指數、Simpson指數和Pielou指數沒有物種豐富度隨經緯度變化趨勢明顯,但結果表現的趨勢一致。

圖2 物種多樣性指數隨緯度和經度梯度的變化,物種豐富度數據均經過自然對數轉換Fig.2 Species diversity indices along the longitude and latitude gradient change. Species richness data were ln-transformed

2.2 物種多樣性指數與各因子的單因素相關分析

物種多樣性指數與各因子的單因素相關性分析結果顯示(表2):物種豐富度與年均溫、最冷月均溫、溫度年較差和潛在蒸散量之間的相關性最強(MAT:R2=0.803,P<0.001; TEMP:R2=0.827,P<0.001; AVMAT:R2=-0.725,P<0.001; PET:R2=0.734,P<0.001),與實際蒸散量、NDVI指數和郁閉度的相關性次之(AET:R2=0.607,P<0.01; NDVI:R2=0.675,P<0.01; CB:R2=0.457,P<0.01),與最暖月均溫、年降水量、最濕潤季降水量和生物量的相關性較低(P<0.05)。Shannon-Wiener指數與年均溫、最冷月均溫和潛在蒸散量的相關性最高(MAT:R2=0.767,P<0.001; TEMP:R2=0.726,P<0.001; PET:R2=0.718,P<0.001),與溫度年較差、實際蒸散量、NDVI指數和郁閉度相關性次之(AVMAT:R2=-0.628,P<0.01; AET:R2= 0.651,P<0.01; NDVI:R2=0.597,P<0.01; CB:R2=0.537,P<0.01),與最暖月均溫、年降水量、最濕潤季降水量和生物量的相關性較弱(P<0.05)。Simpson指數與年均溫、最冷月均溫、溫度年較差、實際蒸散量和潛在蒸散量最強(MAT:R2=0.615,P<0.01; TEMP:R2=0.632,P<0.01; AVMAT:R2=-0.510,P<0.01; AET:R2=0.617,P<0.01; PET:R2=0.466,P<0.01),與最暖月均溫和NDVI指數呈較低顯著正相關(P<0.05)。Pielou指數僅與年均溫、最冷月均溫、實際蒸散量、潛在蒸散量和郁閉度之間具有顯著相關性(P<0.05)。

表2 中國典型森林生態系統物種多樣性指數與各因子之間的的單因素相關關系

*P<0.05; **P<0.01; ***P<0.001; ns,P>0.05. MAT:年均溫annual mean temperature; MTWM: 最暖月均溫mean temperature of the warmest month; MTCM: 最冷月均溫mean temperature of the coldest month; AVMAT: 溫度年較差annual variation of mean annual temperature; AET: 實際蒸散量actual evapotranspiration; PET: 潛在蒸散量potential evapotranspiration; MAP: 年降水量annual mean precipitation; PDQ: 最干燥季降水量precipitation of driest quarter; PWQ: 最濕潤季降水量precipitation of wettest quarter; DNVI: 植被歸一化指數normalized difference vegetation index; CB: 郁閉度crown density; B: 生物量biomass

物種多樣性指數與植物特性或氣候因子的單因素相關性相差較大。其中,物種豐富度和Shannon-Wiener指數與多個因子具有顯著相關性;Simpson指數與能量因子的相關性較高,與水分因子不具有顯著相關性(P>0.05);Pielou指數與各因子之間的相關性相對較低。溫度年較差與物種豐富度、Shannon-Wiener指數和Simpson指數均呈強負相關(P<0.01),但與Pielou指數呈不顯著負相關(R2=-0.239,P>0.05)。其中最干燥季降水量與物種多樣性指數均沒有顯著相關性(P> 0.05)。

2.3 植物特性和氣候因子對物種多樣性分布格局的影響

主成分分析表明,能量因子第一個主成分的特征值為5.101,貢獻率為85.020%,以最冷月均溫(0.990)的載荷量最大;水分因子第一主成分的特征值為2.443,貢獻率為81.432%,以年降水量的載荷量最大,為0.990;植物特性第一主成分的特征值為2.157,貢獻率為71.889%,生物量的載荷量最大,為0.881(表3),其成分矩陣見表4。三類因子的第一主成分指標載荷量較大(載荷量絕對值均大于0.7)(表4),所以分別選擇三類因子第一主成分代表各類因子主要信息的綜合反映。

表3 植物特性、能量與水分因子主成分的特征值和貢獻率

表4 植物特性、能量和水分因子的主成分矩陣

縮寫代表含義同表2

進一步通過方差分解分析,結果表明能量因子和水分因子的共同作用對物種豐富度、Shannon-Wiener指數、Simpson指數和Pielou指數的解釋率最大,分別為42%、38%、26%和15%;植物特性、能量因子和水分因子三者共同作用對物種豐富度、Shannon-Wiener指數、Simpson指數和Pielou指數的解釋率次之,分別為25%、27%、24%和14%。能量因子對物種豐富度、Shannon-Wiener指數和Simpson指數單獨解釋率分別為17%、13%和17%,對Pielou指數沒有產生任何影響。植物特性與能量因子或者水分因子兩者之間的共同作用以及植物特性和水分因子獨立作用對物種多樣性指數空間分布格局的解釋率較小,甚至無影響(圖3)。

圖3 方差分解分析植物特性、能量和水分因子對物種多樣性指數的純效應和共同效應Fig.3 Individual and shared effects of plant traits, energy and moisture factors to species diversity indices as derived from variation partitioning analysis

3 討論

3.1 物種多樣性的經緯度格局

在中國典型森林生態系統中喬木層物種豐富度與緯度相關性最強,即隨緯度升高物種豐富度呈顯著下降,與已有學者[16,43]對我國木本植物和種子植物的物種豐富度與緯度分布規律的研究結果一致。其中,喬木層的Shannon-Wiener指數、Simpson指數和Peilou指數空間分布格局隨緯度呈顯著負相關,與物種豐富度的變化規律一致。在熱帶的西雙版納和尖峰嶺地區,稀有物種和中間物種個體數量較多,物種多樣性指數整體較高。亞熱帶森林受復雜地形因子和水熱條件變化的影響,森林結構比較復雜。鼎湖山、哀牢山和神農架的森林群落郁閉度約在85%—92%,林齡超過200年,都是典型的地帶性頂極森林群落,物種多樣性指數僅次于熱帶森林群落。從暖溫帶到寒溫帶地區,稀有物種逐漸減少,常見物比例增大,群落結構變得簡單,喬木層物種多樣性逐漸減少。Peilou指數隨緯度變化趨勢與中國東部32個森林群落Pielou均勻度指數的研究結果一致[44]。在長白山群落中優勢物種比例較大,稀有物種較少,Peilou指數最高;在海拔較高的哀牢山Peilou指數較高,2—4個偶見種和常見物種個體數較多,稀有物種個體數較少,哀牢山喬木層的群落物種分布比較均勻。

喬木層的物種多樣性指數在經度梯度上波動較大,除了物種豐富度與經度呈顯著負相關(P<0.05)以外,Shannon-Wiener指數、Simpson指數和Pielou指數在經度上均呈不顯著負相關關系。在經度上,喬木層物種多樣性指數在西雙版納和109—113°E之間明顯較高,較低值位于喜馬拉雅-橫斷山區和北方地區。結果表明在中國典型森林生態系統中喬木層物種多樣性南北方向的聯系強度大于東西方向的聯系,與中國高山植物區系的地理空間格局研究一致[45]。

3.2 氣候因子對物種多樣性格局的影響

水熱動態假說認為物種多樣性的大尺度格局由水分和能量共同決定。這一假說通常用潛在蒸散量和溫度表示熱能,而水分一般指液態水,通常用一個地區的年降水量表示,也有用實際蒸散量反映一個地區的水熱動態平衡[21- 22]。能量和水分的空間分異越大,意味著在能量和水分梯度上存在越多樣的生境滿足來不同物種水熱的需要,為植物造就更加多樣的生境,從而支持更多的物種得以共存和分布,故物種多樣性較高。方差分解結果表明,能量和水分因子的共同作用對物種豐富度、Shannon-Wiener指數、Simpson指數和Pielou指數的解釋率分別達到42%、38%、26%和15%,這在一定程度上解釋能量和水分因子共同作用對物種多樣性指數的空間分布格局起決定作用,進而影響物種多樣性分布格局,支持了水熱動態假說。

能量因子對喬木層的物種豐富度、Shannon-Wiener指數和Simpson指數單獨解釋率分別為17%、13%和17%,在單因素相關分析中年均溫、最冷月均溫和潛在蒸散量差異均顯著影響物種多樣性指數的分布格局(P<0.05),而且相對于其他因子的空間分異而言,其相關性更高。所以環境能量假說適合解釋中國典型森林生態系統喬木層物種多樣性的分布格局。寒冷忍耐假說認為,在寒冷地區,很多物種由于不能忍受冬季的寒冷而無法生存。因此,隨著冬季溫度的降低,其物種多樣性逐漸減少[38,46]。本研究結果表明最冷月均溫與物種多樣性指數均呈強正相關,說明冬季低溫對物種多樣性的分布格局構成了很強的限制作用,比如在大興安嶺最低溫度為-29℃,主要喬木樹種只有興安落葉松(Larixgmelinii)和白樺(Betulaplatyphylla);由于絕對大部分物種的祖先都是在濕熱的環境中進化,極端環境條件導致植物個體水平的死亡,通過影響物種的種群動態,最終影響到物種的分布范圍和物種多樣性的分布格局[47- 48]。本研究結果同樣支持寒冷忍耐假說。

根據氣候穩定性假說,在氣候穩定地區物種的生態位相對較窄,所以就有更多的空間容納更多的物種;反之,在氣候波動較少的地區,可能資源的供應相對而言較不穩定,物種的生態位相對較寬,能容納的物種相對較少,所以穩定的自然資源更有利于物種的進化與適應[27- 28]。這一假說通常用溫度年較差來表示,說明該地區溫度的季節波動性,用于表征氣候的穩定性。物種豐富度、Shannon-Wiener指數和Simpson指數均與溫度年較差有很強的負相關(P<0.01),Pielou指數與溫度年較差呈不顯著負相關。這說明溫度年較差在很大程度上決定了物種多樣性的空間分布格局。隨著緯度的增加,溫度年較差呈遞減趨勢,這也解釋了本研究中喬木層的物種多樣性由南向北的遞減趨勢。

3.3 植物特性對物種多樣性格局的影響

生產力假說常采用歸一化植被指數、實際蒸散量和年降水量作為生產力的替代指標[29- 31]。單因素相關分析表明,物種豐富度與NDVI指數、實際蒸散量和年降水量均有顯著正相關(P<0.05),其中物種豐富度與NDVI指數的相關性最顯著(R2=0.675,P<0.01),表明NDVI指數與實際蒸散量或年降水量相比能更好的解釋物種豐富度的大尺度格局,支持物種豐富度格局的生產力假說。NDVI指數與Shannon-Wiener指數和Simpson指數也具有顯著相關性。由于Pielou指數反映群落內的個體密度和均勻情況,NDVI指數反映植被覆蓋程度,Pielou指數和NDVI指數呈不顯著相關性。此外,郁閉度通過改變林內光照條件、植物物種組成和群落結構,從而形成森林群落物種共存以及維持多樣性的重要場所,但郁閉度和物種多樣性不存在線性關系[49]。南方森林群落結構復雜,森林群落結構動態變化大,而北方喬木層居優勢地位,群落結構穩定,郁閉度相對南方較低。在13個典型森林喬木層的群落郁閉度與物種豐富度、Shannon-Wiener指數和Pielou指數均呈顯著正相關,與Simpson指數呈不顯著的正相關。生物量與物種豐富度和Shannon-Wiener指數均呈顯著正相關(P<0.05),生物量高的植物群落其物種多樣性也相對較高,這說明群落中的物種多樣性和生物量是相互影響。在一定程度上促進森林群落生物量的發展,同時也有利于提高群落物種多樣性增加,植物多樣性的保護應優先保護植物生物量較高的區域。進一步通過方差分解發現,植物特性、能量因子和水分因子三者共同作用對喬木層的物種多樣性指數能解釋到14%—27%的變異,植物特性對喬木層的物種豐富度和Shannon-Wiener指數單獨解釋率分別為7%和4%。所以,在解釋物種多樣性大尺度分布格局時,也必須充分考慮植物因素的影響。

綜上所述,本研究主要考慮植物特性和氣候因子的影響,仍有一些影響物種多樣性空間分布格局的因素未被涉及,如微環境(坡度、坡向、海拔等)、土壤因子(有機碳、全氮、全磷等)、邊緣效應、區域歷史、生境異質性、人為干擾以及種子傳播等。未來應該綜合考慮各種生態因子,從而更加全面的解釋中國典型森林生態系統物種多樣性空間分布格局及其影響因素。

4 結論

基于中國典型森林生態系統喬木層群落調查數據,結合植物特性和氣候數據,研究大尺度物種多樣性的分布格局及其影響因素。主要結論如下:

(1)喬木層的物種多樣性指數在地理格局上存在較大差異,均隨經緯度上升而下降。其中,物種多樣性指數在緯度梯度上均呈顯著相關性,在經度上波動較大,表明喬木層物種多樣性在緯度梯度上變化規律更為明顯。

(2)物種多樣性指數與各因子的單因素相關性相差較大。物種豐富度、Shannon-Wiener指數和Simpson指數與能量因子相關性較強,Pielou指數僅與年均溫、最冷月均溫、實際蒸散量、潛在蒸散量和郁閉度存在顯著相關性。其中,Simpson指數、Pielou指數與水分因子無顯著相關性,物種豐富度和Shannon-Wiener指數均與DNVI指數、郁閉度和生物量具有顯著的相關性。

(3)植物特性和氣候因子對喬木層的物種多樣性指數分布格局的解釋率存在較大差異。能量和水分因子的共同作用對喬木層物種多樣性指數的解釋率最高,植物特性、能量和水分因子三者的共同作用解釋率次之。其中,能量因子單獨解釋率高于植物特性或水分因子,植物特性與能量因子或水分因子之間的共同作用以及植物特性和水分因子單獨作用對喬木層的物種多樣性指數解釋率較小,甚至無影響。

致謝：感謝國家生態系統觀測研究網絡(CNERN)、中國生態系統研究網絡(CERN)以及野外觀測站的支持。