胡 芳,曾馥平,杜 虎,彭晚霞,張 芳，譚衛(wèi)寧，宋同清*

1 中國(guó)科學(xué)院亞熱帶農(nóng)業(yè)生態(tài)研究所亞熱帶農(nóng)業(yè)生態(tài)過(guò)程重點(diǎn)實(shí)驗(yàn)室, 長(zhǎng)沙 410125 2 中國(guó)科學(xué)院環(huán)江喀斯特生態(tài)系統(tǒng)觀測(cè)研究站, 環(huán)江 547100 3 廣西壯族自治區(qū)木論國(guó)家級(jí)自然保護(hù)區(qū)管理局, 環(huán)江 547200

物種多樣性一般是指物種豐富度、多度或者是兩者在特定區(qū)域下的組合[1]。物種多樣性是當(dāng)前生態(tài)學(xué)研究中十分重要的內(nèi)容和熱點(diǎn)之一[2],不同的取樣尺度下物種多樣性可能不同[3],因?yàn)榭臻g尺度對(duì)物種多樣性通常具有非線性影響,此外,異質(zhì)性隨著空間尺度的增加而增大,在空間尺度增加的過(guò)程中常會(huì)導(dǎo)致較小尺度的信息丟失[4]。不同的多樣性指數(shù)自身性質(zhì)不同,如物種豐富度和多度在不同空間尺度反應(yīng)不同。一個(gè)尺度適合這個(gè)指數(shù),但是不一定適合另外的指數(shù)。盡管生態(tài)學(xué)家現(xiàn)在敏銳地意識(shí)到,將信息從一個(gè)尺度轉(zhuǎn)移到另一個(gè)尺度是必不可少的,并在過(guò)去的20多年已經(jīng)取得了很大進(jìn)步[5- 6],尺度效應(yīng)仍是一個(gè)具有挑戰(zhàn)性的問(wèn)題,因?yàn)槠鋵?duì)生物多樣性研究以及保護(hù)具有重大意義[7]。因此,分析不同空間尺度內(nèi)物種多樣性將有助于自然生態(tài)系統(tǒng)的保護(hù)與管理[1]。

物種多樣性在不同氣候帶的地理位置上表現(xiàn)出明顯的變化,例如,西雙版納20 hm2的熱帶森林動(dòng)態(tài)監(jiān)測(cè)樣地有468種[8],而長(zhǎng)白山闊葉紅松林25 hm2森林大樣地僅51種[9]。即使在相似的氣候帶下,物種多樣性也有差異,如茂蘭喀斯特原始森林的1 hm2的樣地中共有木本植物63科140屬199種[10],而木論喀斯特原始森林的0.8 hm2樣地僅有43科91屬123種[11]。在相同的研究區(qū)域隨不同的取樣尺度的變化,物種多樣性也不同[12- 16]。因此,研究物種多樣性時(shí)非常有必要考慮到尺度[7]。

在熱帶雨林[12- 13]、南亞熱帶常綠闊葉林[14]、溫帶針闊混交林[3]以及溫帶次生林[15]中開(kāi)展的研究表明物種多樣性具有尺度依賴性[12- 14]?？λ固氐貐^(qū)普遍具有景觀異質(zhì)性高、土壤總量少、土層薄、巖石裸露率高、土層貧瘠等特點(diǎn),生境和植被具有高度的異質(zhì)性,植被隨海拔變化垂直現(xiàn)象明顯[11,16- 17],已有張忠華等[10]對(duì)喀斯特地區(qū)兩個(gè)1 hm2的樣地的多樣性進(jìn)行了研究,在更大尺度上的研究顯得尤為重要。本文以木論喀斯特常綠落葉闊葉混交林的25 hm2動(dòng)態(tài)監(jiān)測(cè)樣地為研究對(duì)象,探討了不同空間尺度下各物種多樣性指標(biāo)的空間分布特征,以期為喀斯特森林在空間尺度和空間位置上的多樣性模式提供重要的見(jiàn)解,同時(shí)為物種的多樣性的管理和保護(hù)提供一個(gè)合理的科學(xué)基礎(chǔ)。

1 材料與方法

1.1 研究區(qū)概況

木論國(guó)家自然保護(hù)區(qū)(25°07′01″—25°12′22″ N,107°54′01″—108°05′51″ E) 位于廣西環(huán)江毛南族自治縣的西北部的喀斯特峰叢洼地集中連片區(qū),南北長(zhǎng)10.75 km,東西寬19.80 km,總面積約89.69 km2,屬中亞熱帶季風(fēng)氣候,年均日照為1451 h,無(wú)霜期為310天,年平均氣溫為19.3 ℃,年降雨量豐沛,海拔420 m處年均降雨量為1529.2 mm(海拔420 m),平均蒸發(fā)量為1571.1 mm,平均相對(duì)濕度為79%。研究區(qū)屬于亞熱帶石灰?guī)r區(qū)常綠落葉闊葉混交林生態(tài)系統(tǒng),是我國(guó)生物區(qū)系相匯交錯(cuò)區(qū)的中心,生境異質(zhì)性極高,區(qū)系成分復(fù)雜,物種種類豐富,生態(tài)環(huán)境脆弱,是目前世界上喀斯特地區(qū)保存最完好、面積最大的原生林。林區(qū)巖石裸露率高達(dá)80%,甚至部分區(qū)域達(dá)到90%以上,土層淺且分布不均,區(qū)內(nèi)主要以碳酸鹽巖發(fā)育的深色或棕褐色石灰土為主,典型景觀單元為峰叢洼地。

1.2 樣地設(shè)置及調(diào)查

圖1 木論樣地等高線圖 Fig.1 Contour map of Mulun forest plot

自2013年初經(jīng)過(guò)多次實(shí)地勘察,反復(fù)測(cè)量并調(diào)整樣地邊框,按照CTFS(Center for Tropical Forest Science)全球森林生物多樣性監(jiān)測(cè)規(guī)范[18],于2014年在木論自然保護(hù)區(qū)內(nèi)建設(shè)完成25 hm2森林動(dòng)態(tài)監(jiān)測(cè)樣地。樣地為正方形,東西、南北邊長(zhǎng)均為500 m(圖1)[19]。整個(gè)樣地被劃分成625個(gè)20 m×20 m樣方,樣方4角用水泥樁作永久標(biāo)記。每個(gè)20 m ×20 m樣方又細(xì)分為16個(gè)5 m×5 m小樣方。野外調(diào)查以20 m×20 m樣方為單元,記錄了其中所有胸徑≥1 cm的木本植物樹(shù)種名稱、胸徑、坐標(biāo),計(jì)劃今后每隔5年復(fù)查一次。首次調(diào)查結(jié)果表明,樣地共有胸徑(DBH)≥1 cm的獨(dú)立木本植物個(gè)體108,667株,隸屬61科147屬227種。

1.3 數(shù)據(jù)處理與分析

本研究將整個(gè)25 hm2樣地劃分為以下6個(gè)空間尺度：0.0025 hm2(5 m×5 m),0.01 hm2(10m×10 m),0.04 hm2(20 m×20 m),0.25 hm2(50 m×50 m),1 hm2(100 m×100 m)和6.25 hm2(250 m×250 m)。6個(gè)空間尺度的樣方數(shù)見(jiàn)表1。分別計(jì)算6個(gè)空間尺度中樣方內(nèi)的物種多樣性指數(shù),包括：多度、豐富度、Shannon-Wiener指數(shù)、Simpson指數(shù)以及Pielou均勻度指數(shù)等。

各物種多樣性指數(shù)[20]、方差、變異系數(shù)計(jì)算公式如下：

1)豐富度指數(shù)S(Species richness)：S=樣方內(nèi)物種數(shù)

4)Pielou均勻度指數(shù)計(jì)算公式為：E=H/ln(S)

6)變異系數(shù)：CV=(sd/μ)×100%

式中,Pi為單個(gè)樣方內(nèi)第i個(gè)物種的多度占所有多度之和的比例,S為單個(gè)樣方內(nèi)的物種數(shù),n為樣方總數(shù),Xj為計(jì)算第j個(gè)樣方物種多樣性指標(biāo)的數(shù)值,μ為平均值,sd為標(biāo)準(zhǔn)差。

根據(jù)計(jì)算結(jié)果,繪制各多樣性指數(shù)的空間分布圖,以對(duì)整個(gè)樣地不同取樣尺度下的物種多樣性空間分布特點(diǎn)及空間變化進(jìn)行分析。計(jì)算物種豐富度和多度在每一取樣尺度下的Spearman秩相關(guān)系數(shù),探討物種豐富度和多度的相關(guān)性隨尺度變化的分異。同時(shí)利用多樣性指標(biāo)的方差和變異系數(shù)對(duì)不同取樣面積的自然對(duì)數(shù)進(jìn)行作圖,以探討取樣尺度的變化對(duì)物種多樣性格局的影響。

數(shù)據(jù)統(tǒng)計(jì)分析和物種多樣性空間分布格局制圖均通過(guò)R 3.3.1軟件完成,物種豐富度與多度Spearman秩相關(guān)檢驗(yàn)通過(guò)SPSS 21.0軟件完成,物種多樣性指數(shù)空間變異制圖通過(guò)ORINGIN 8.0軟件完成。

2 結(jié)果

2.1 多度和物種豐富度的空間分布格局

圖2 多度在6種取樣尺度下的空間分布特征Fig.2 Spatial distribution of abundance index at six scales

圖3 物種豐富度在6種取樣尺度下的空間分布特征Fig 3 Spatial distribution of species richness index at six scales

多度和物種豐富度隨著取樣空間尺度的增加而增加(圖2,圖3)。在小尺度上(5 m×5 m至20 m×20 m)上,多度和物種豐富度的空間分布格局表現(xiàn)出相同的趨勢(shì),即樣地中部洼地區(qū)域的多度和物種豐富度明顯低于東部山脊以及西北角；在中等尺度(50 m×50 m),這種分布差異不明顯；在較大尺度(100 m×100 m至250 m×250 m)上,多度和物種豐富度的空間分布格局明顯不同,在250 m×250 m的尺度上豐富度的最大值出現(xiàn)在東南角,而多度的最大值出現(xiàn)在東北角。此外,多度和物種豐富度的秩相關(guān)分析結(jié)果表明,在5 m×5 m至20 m×20 m尺度上,多度和物種豐富度呈現(xiàn)出極顯著(P<0.001)的正相關(guān)關(guān)系,當(dāng)尺度超過(guò)100 m×100 m時(shí),兩者的相關(guān)性則不再顯著(表1)。

表1 多度和物種豐富度的Spearman秩相關(guān)分析結(jié)果

2.2 物種多樣性指數(shù)和均勻度指數(shù)的空間分布格局

從圖4、圖5、圖6中可以看出,Shannon-Weiner指數(shù)、Simpson指數(shù)和Pielou均勻度指數(shù)等在不同的尺度上表現(xiàn)出相同的趨勢(shì)。在最大尺度250 m×250 m(6.25 hm2)上,3種多樣性指數(shù)均是樣地東南部的最高,并且西半部多樣性明顯地于東半部；在中等尺度(100 m×100 m至20 m×20 m)上,3個(gè)指數(shù)均在在西北角地區(qū)現(xiàn)了高值,整體表現(xiàn)為山坡部分高于洼地、山頂區(qū)域,符合“中間高度膨脹(mid-altitude bulge)”理論。在小尺度上(10 m×10 m至5 m×5 m),每個(gè)樣方下3個(gè)指數(shù)的數(shù)值都很低,而且十分接近,甚至很難從灰度圖中區(qū)分出數(shù)值較高的樣方。

圖4 Shannon-Weiner指數(shù)在6種取樣尺度下的空間分布特征Fig.4 Spatial distribution of Shannon-Weiner index at six scales

圖5 Simpson指數(shù)在6種取樣尺度下的空間分布特征Fig.5 Spatial distribution of Simpson index at six scales

圖6 Simpson指數(shù)在6種取樣尺度下的空間分布特征Fig.6 Spatial distribution of Simpson index at six scales

2.3 多度和物種豐富度空間變異隨尺度的變化規(guī)律

多度和物種豐富度的方差和變異系數(shù)隨著取樣尺度的增加有著類似的變化趨勢(shì)(圖7)。多度和物種豐富度都隨著尺度的增加表現(xiàn)出單峰型的變化特征,且在100 m×100 m尺度上達(dá)到最大值。多度和物種豐富度的變異系數(shù)均隨著尺度的增加呈現(xiàn)出線性減小的趨勢(shì),而多度在小于100 m×100 m尺度上相對(duì)于物種豐富度較線性平緩減小。

圖7 多度和物種豐富度在6種尺度下的空間變異Fig.7 Spatial variation of species richness and abundance at six scales

2.4 物種多樣性指數(shù)和均勻度指數(shù)空間變異隨尺度的變化規(guī)律

隨著取樣尺度的增加,Shannon-Weiner指數(shù)、Simpson指數(shù)和Pielou均勻度指數(shù)的方差都有著明顯不同的變化趨勢(shì)(圖8)。Shannon-Weiner指數(shù)、Simpson指數(shù)和Pielou均勻度指數(shù)的方差隨著尺度的增加均表現(xiàn)出單峰型的變化特征,且都是在100 m×100 m尺度上達(dá)到最大值。Shannon-Weiner指數(shù)、Simpson指數(shù)和Pielou均勻度指數(shù)的變異系數(shù)都隨著尺度的增加呈現(xiàn)出線性減小的趨勢(shì),且在5 m×5 m至20 m×20 m尺度上急劇下降,之后則比較平緩。

圖8 Shannon-Weiner指數(shù)、Simpson指數(shù)和Pielou均勻度指數(shù)在6種尺度下的空間變異Fig.8 Spatial variation of Shannon-weiner,Simpson and Pielou evenness index at six scales

3 討論

物種豐富度指數(shù)、Shannon-Weiner指數(shù)、Simpson指數(shù)、 Pielou均勻度指數(shù)等均沒(méi)有疊加的特性,特定的空間尺度上觀察到的高的值可能在另一個(gè)空間尺度上表現(xiàn)為低,這是因?yàn)楫?dāng)尺度增大時(shí),小尺度的物種豐富度不能簡(jiǎn)單的疊加。物種多樣性指數(shù)在中尺度(100 m×100 m至20 m×20 m)上,整體表現(xiàn)為山坡高于洼地區(qū)域和山頂,這與前人研究結(jié)果相同,木論喀斯特森林木本植物物種豐富度沿海拔梯度的曲線符合單峰格局,峰值出現(xiàn)在坡腰位置[21]。這可能是由于樣地處于亞熱帶區(qū)域,海拔最高為651.4 m,溫度隨海拔的升高下降不明顯,仍然適合植物的生長(zhǎng),低海拔洼地地區(qū),土層較深厚、濕度條件較好,但太陽(yáng)直射光照時(shí)間不長(zhǎng),僅限于一些喜蔭和耐蔭性強(qiáng)的物種,導(dǎo)致物種較少；高海拔的山頂,土層較薄,山體陡峭,日照時(shí)間長(zhǎng)導(dǎo)致土壤中的水分蒸發(fā)旺盛,不利于植物生長(zhǎng),因而物種較少；而中等海拔的山坡地帶,土層厚度、溫濕條件都相對(duì)而言適合植物生長(zhǎng),因而物種相對(duì)豐富[21]。物種豐富度和多度指數(shù)在小尺度、中尺度以及的尺度上表現(xiàn)出明顯的差異性,這是因?yàn)樾〕叨壬衔锓N多樣性主要是與樣方內(nèi)木本植物個(gè)體數(shù)量有關(guān),隨著尺度的增加,新生境的出現(xiàn)導(dǎo)致物種多樣性隨著空間尺度的變化表現(xiàn)出明顯的差異性[22]。本研究結(jié)果表明物種豐富度和多度之間的相關(guān)性隨空間取樣尺度的增加逐漸變得不再顯著,類似的趨勢(shì)同樣出現(xiàn)在其他地區(qū)的熱帶山地雨林[13]、溫帶森林[3]。在本研究表明喀斯特常綠落葉闊葉混交林物種多樣性指數(shù)均隨著尺度的變化均表現(xiàn)出不連續(xù)性。這一結(jié)果意味著在生物多樣性保護(hù)和管理時(shí)所謂的生物多樣性“熱點(diǎn)”地區(qū)的評(píng)估需要在多個(gè)尺度上進(jìn)行評(píng)估[7]。另外,本研究結(jié)果顯示,不同的物種多樣性指數(shù)得出的結(jié)果不完全相同,這與多樣性指數(shù)的本身性質(zhì)有關(guān),其中物種豐富度和多度主要受生境異質(zhì)性(坡向和巖石裸露率)的影響[23],Shannon-Weiner多樣性指數(shù)在群落類型、物種個(gè)體密度和生境差異等方面的反映結(jié)果較好,Simpson優(yōu)勢(shì)度指數(shù)是反映群落優(yōu)勢(shì)度的較好指標(biāo),而Pielou均勻度指數(shù)受物種豐富度的影響,即對(duì)樣方面積的大小較敏感[24]。但就目前的各種物種多樣性指數(shù)都不能完全客觀地反映群落物種多樣性的現(xiàn)狀[25],Iudin等認(rèn)為可以利用多維分形的方法來(lái)選擇合適的物種多樣性指數(shù)[26]。在木論喀斯特常綠落葉闊葉混交林木本植物多樣性的測(cè)定中,本研究所采用的物種多樣性指數(shù),盡管表現(xiàn)形式有所不同,但都可以明確反映物種數(shù)量的變化。

本文中用方差和變異系數(shù)來(lái)描述物種多樣性指數(shù)的空間變異,較好地描述了木論喀斯特常綠落葉闊葉混交林的物種多樣性空間分布特征。本研究中,物種多樣性指數(shù)的方差隨尺度的增加都表現(xiàn)出單峰分布的特征,在100 m×100 m尺度上達(dá)到最大值,與物種豐富度和多度格局分離尺度相一致,這也印證了生態(tài)學(xué)家通過(guò)方差分析中產(chǎn)生的拐點(diǎn)來(lái)確定相應(yīng)尺度下的生態(tài)過(guò)程的合理性。以往有研究認(rèn)為這種單峰特征可能與樣地內(nèi)不同的稀有種組成有很大的關(guān)系[27],本樣地中有稀有種有100種,占總樹(shù)種數(shù)的39.37%[22],而海南尖峰嶺樣地內(nèi)的稀有物種數(shù)量為60種,占總物種數(shù)的20.7%[27],該區(qū)域多樣性研究中沒(méi)有出現(xiàn)類似的現(xiàn)象[12]；類似的單峰曲線在其他研究中也出現(xiàn)過(guò)[14],但是不同的森林類型其峰值的尺度不同。例如：海南尖峰嶺熱帶山地雨林中物種豐富度的方差在20 m×20 m尺度上達(dá)到峰值[13],而在長(zhǎng)白山溫帶森林中物種豐富度的峰值出現(xiàn)在125 m×125 m尺度上[3]。在類似的森林類型中,張忠華等研究表明物種豐富度的方差在25 m×25 m尺度上出現(xiàn)峰值[10],這可能是由于不同的生境異質(zhì)性導(dǎo)致的。同時(shí),Shannon-Weiner指數(shù)、Simpson指數(shù)、Pielou均勻度指數(shù)的變異系數(shù)在5 m×5 m至20 m×20 m尺度上線性急劇減小,之后則變化平緩,表明小尺度比大尺度有更高的可變性。以上表明該區(qū)域的物種多樣性具有區(qū)域依賴特征,和鼎湖山亞熱帶常綠闊葉林[14]、長(zhǎng)白山溫帶森林[9]以及海南尖峰嶺熱帶森林[13]樣地的研究結(jié)果一致。

張忠華等曾研究表明喀斯特常綠落葉闊葉混交林大多數(shù)物種是聚集分布的,且聚集度隨空間尺度的變化而變化[28]；本樣地前人研究結(jié)果表明也有類似的結(jié)果[29],這在其他研究區(qū)域也是常見(jiàn)的[30-31]。一些研究者認(rèn)為隨機(jī)分布和規(guī)則分布取決于種內(nèi)或種間對(duì)于資源的競(jìng)爭(zhēng),而聚集分布則主要是生境異質(zhì)性和種子傳播方式導(dǎo)致的[32-33]。木論25 hm2喀斯特常綠落葉闊葉混交林樣地的樣地海拔在442.6—651.4 m之間,最大高差達(dá)208.8 m。坡度范圍0.12°—66.97°,地形起伏較大,生境的異質(zhì)程度相當(dāng)高,導(dǎo)致物種多樣性在不同的尺度上分布格局不同。對(duì)于喀斯特常綠落葉闊葉混交林來(lái)說(shuō),植被類型、海拔高度、人為干擾、演替階段、生境條件差異均對(duì)物種多樣性指數(shù)產(chǎn)生一定的影響[34-35]。因此,想要了解物種多樣性隨空間尺度的變化模式,我們需要在類似的森林生態(tài)系統(tǒng)做更多的研究。

致謝：感謝廣西壯族自治區(qū)木論國(guó)家自然保護(hù)區(qū)管理工作人員和中南林業(yè)科技大學(xué)2012級(jí)生態(tài)學(xué)專業(yè)同學(xué)對(duì)樣地調(diào)查給予的幫助。