詹小豪，王旭航，葉諾楠，吳初平，袁位高，伊力塔

（1. 浙江農(nóng)林大學(xué) 林業(yè)與生物技術(shù)學(xué)院，浙江 杭州 311300；2. 中國電建集團(tuán)華東勘測設(shè)計(jì)研究院有限公司，浙江 杭州 311122；3. 浙江省林業(yè)科學(xué)研究院，浙江 杭州 310023）

植物種群的空間分布格局是指種群在水平空間上的配置和分布狀況[1]，除了種群自身特性和環(huán)境條件因素與其形成密切相關(guān)外[2]，種間關(guān)系也是促使其形成的主要?jiǎng)恿χ籟3]，在判斷植物種間關(guān)系時(shí)，生態(tài)學(xué)家們通常會運(yùn)用種間聯(lián)結(jié)性分析來探究物種之間的內(nèi)在聯(lián)系[4]。此外，生態(tài)位的重疊程度在一定程度上能夠反映物種聯(lián)結(jié)關(guān)系以及空間配置關(guān)系[5?6]。對植物種間關(guān)系及空間分布格局的研究，有助于認(rèn)識種群與生境的相互關(guān)系、空間資源獲取能力與生態(tài)適應(yīng)對策，預(yù)測群落消長動(dòng)態(tài)，是深入了解維持樹木物種共存機(jī)制，并了解產(chǎn)生空間格局的過程的重要手段[7?9]，對于正確認(rèn)識群落的組成、功能及演替規(guī)律具有重要意義。次生林是中國森林資源的主體，它既保持著原始森林的物種組成成分與生境，又與原始森林在結(jié)構(gòu)組成、林木生長、生產(chǎn)力、林分環(huán)境和生態(tài)功能等諸多方面有著顯著的不同，原始林退化導(dǎo)致的次生林面積擴(kuò)大將引起森林生態(tài)系統(tǒng)中生物多樣性的下降[10]。如何恢復(fù)和保護(hù)次生林群落的物種多樣性成為生態(tài)學(xué)家面臨的重要問題。而對森林木本植物的空間分布格局及種間關(guān)系的研究有助于揭示群落結(jié)構(gòu)的形成機(jī)制與潛在的生態(tài)學(xué)過程，對次生林的經(jīng)營撫育具有一定指導(dǎo)意義[11]。以松林、松闊混交林和常綠闊葉林中幼齡林為主體的天然次生林是浙江省建德市的森林資源主體，主要由被過度干擾破壞的天然林地逐漸恢復(fù)、演替而來[12]。次生林中仍保有原始森林的部分物種，但不同類型次生林在群落結(jié)構(gòu)、林木年齡組成、植物生產(chǎn)力及生境等方面都存在顯著差異[13]。本研究對浙江建德的次生林群落進(jìn)行每木調(diào)查，分析主要樹種間的相互關(guān)系及空間分布格局特征，以增強(qiáng)對該區(qū)域次生林群落結(jié)構(gòu)特征和空間分布格局的認(rèn)識，促進(jìn)區(qū)域生物多樣性保護(hù)與木本植物資源的可持續(xù)利用，預(yù)測群落的演替方向并探索其驅(qū)動(dòng)力來源，進(jìn)而為該區(qū)域次生林群落的恢復(fù)、改造及撫育經(jīng)營提供依據(jù)，更好發(fā)揮其生態(tài)社會效益。

1 研究區(qū)概況

研究區(qū)位于浙江省建德市新安江林場(29°29′N，119°16′E)，屬亞熱帶北緣季風(fēng)氣候，溫暖濕潤，四季分明，年平均氣溫為16.9 ℃，最冷月平均氣溫為4.7 ℃，最熱月平均氣溫為29.2 ℃，年均降水量為1 504.0 mm。據(jù)《中國植被》區(qū)劃，該地區(qū)森林植被屬亞熱帶常綠闊葉林北部亞地帶，地帶性植被為常綠闊葉林。主要森林植被類型有常綠闊葉林、落葉闊葉林、常綠落葉闊葉混交林、針闊混交林等[14]。

2 研究方法

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

在實(shí)地踏查的基礎(chǔ)上，在浙江建德典型常綠闊葉林、松闊混交林和松林群落中分別設(shè)置面積為100 m×100 m的樣地(表1)，用木樁、塑料帶進(jìn)行圍封標(biāo)記。將每塊樣地劃分為25小塊 (20 m×20 m)，并以此為基本單位對喬木層進(jìn)行每木調(diào)查。記錄樣地群落類型、海拔、坡度、坡向、土壤等環(huán)境因子；記錄樣方中胸徑≥5 cm的木本植物的基礎(chǔ)數(shù)據(jù)，包括種名、樹高、胸徑、冠幅及其相對應(yīng)的坐標(biāo)[15?16]。

2.2 重要值

重要值(IV)計(jì)算[17]：IV=(相對多度+相對顯著度+相對頻度)/3×100%。

2.3 生態(tài)位重疊

采用Pianka重疊指數(shù)計(jì)算物種間的生態(tài)位重疊系數(shù)，公式[18]如下：

其中：Oik為種i和種k的生態(tài)位重疊系數(shù)[1]；w為劃分的資源位總數(shù)；Pij和Pkj分別是種i和種k在資源位j中的重要值占該物種在整個(gè)資源中總重要值的比例。Oik的值域?yàn)閇0，1]，當(dāng)2個(gè)種對群落中所有的資源都不存在共享狀態(tài)時(shí)，該種對之間生態(tài)位完全不重疊，其值為0；當(dāng)2個(gè)種對群落中所有資源利用完全重疊時(shí)，該種對的生態(tài)位重疊程度為1，達(dá)到最大狀態(tài)。

2.4 空間分布格局

2.4.1 種內(nèi)空間分布格局 Ripley’sK(r)函數(shù)是進(jìn)行種內(nèi)空間分布格局分析的基礎(chǔ)函數(shù)，公式如下[19]：

其中：r為分析的空間尺度(m)，A為樣方面積(m2)，n為植物種個(gè)體數(shù)，upq表示點(diǎn)p和點(diǎn)q之間的距離(m)；Ir(upq)為指示函數(shù)，當(dāng)upq≤r時(shí)，Ir(upq)=1，當(dāng)upq＞r時(shí)，Ir(upq)=0；epq為權(quán)重值，用于邊緣校正。為了更直觀地解釋實(shí)際的空間格局，通常用Ripley’sL(r)函數(shù)表示：

當(dāng)L(r)=0，種群分布類型為隨機(jī)分布，對種群的聚集分布研究基于隨機(jī)分布的基礎(chǔ)，以樣地中的任意一點(diǎn)為圓心，r為半徑畫圓，如果在圓中出現(xiàn)的個(gè)體數(shù)多于隨機(jī)狀態(tài)下的個(gè)體數(shù)，那么表示該種群呈現(xiàn)聚集分布；當(dāng)圓中出現(xiàn)的個(gè)體數(shù)少于隨機(jī)分布狀態(tài)下的個(gè)體數(shù)時(shí)，該種群呈現(xiàn)均勻分布[20?22]。采用Monte-Carlo模擬99%置信區(qū)間，進(jìn)行結(jié)果偏離隨機(jī)狀態(tài)的顯著性檢驗(yàn)。若L(r)值位于置信區(qū)間之上，種群呈聚集分布，L(r)值位于置信區(qū)間之下，種群呈均勻分布，L(r)值位于置信區(qū)間之內(nèi)，種群呈隨機(jī)分布[23]。

2.4.2 種間空間分布格局 種間空間分布格局基礎(chǔ)函數(shù)計(jì)算公式[24]如下：

其中：n1和n2分別為種1和種2的個(gè)體數(shù)，f和g分別代表種1和種2的個(gè)體。同樣，用L12(r)取代K12(r)，公式為：

當(dāng)L12(r)=0，表明2個(gè)種在r尺度下無關(guān)聯(lián)；當(dāng)L12(r)＞0，表明兩者為空間正關(guān)聯(lián)；當(dāng)L12(r)＜0，表明兩者為空間負(fù)相關(guān)。采用Monte-Carlo模擬99%置信區(qū)間，進(jìn)行結(jié)果偏離隨機(jī)狀態(tài)的顯著性檢驗(yàn)。當(dāng)L12(r)值位于置信區(qū)間之上，2個(gè)變量顯著正相關(guān)；L12(r)值位于置信區(qū)間之下，2個(gè)變量顯著負(fù)相關(guān)；L12(r)值位于置信區(qū)間之內(nèi)，2個(gè)變量相互獨(dú)立[23]。本研究根據(jù)實(shí)際樣地面積，在參考同類研究及毗鄰地區(qū)森林群落空間格局研究的基礎(chǔ)上[25?26]，將格局分析的尺度限定為0～25.0 m。

2.5 種間聯(lián)結(jié)

2.5.1 總體聯(lián)結(jié)性檢驗(yàn) 采用方差比率法來測定各類型次生林群落中主要樹種間的總體聯(lián)結(jié)性，并利用統(tǒng)計(jì)量W來檢驗(yàn)總體聯(lián)結(jié)是否顯著，計(jì)算公式[27]如下：

其中：VR為方差比率，N為樣方總數(shù)，為樣方內(nèi)出現(xiàn)的目標(biāo)物種總數(shù)；t為樣方中物種的平均數(shù)，S為總物種數(shù)，Pi為物種i的頻度。以VR作為不同類型次生林中主要樹種總體聯(lián)結(jié)性指數(shù)，在獨(dú)立性零假設(shè)條件下，VR期望值為1，即當(dāng)VR=1時(shí)認(rèn)為種間無聯(lián)結(jié)；若VR＜1則表示物種間存在正聯(lián)結(jié)；若VR＜1則表示物種間存在負(fù)聯(lián)結(jié)。對于VR偏離1的顯著程度采用統(tǒng)計(jì)量W來驗(yàn)證，W=VRN，若W＜則說明種間總體聯(lián)結(jié)性性顯著，若則說明種間總體聯(lián)結(jié)性不顯著[28]。

2.5.2 種間聯(lián)結(jié)性分析 采用經(jīng)Yates連續(xù)校正系數(shù)糾正的χ2統(tǒng)計(jì)量對種間聯(lián)結(jié)性進(jìn)行定性研究。公式[18]如下：

2.5.3 種間相關(guān)性測定 本研究以多度作為Spearman秩相關(guān)系數(shù)的數(shù)量指標(biāo)，對種對間的線性關(guān)系做定量分析，計(jì)算公式[31]如下：

3 結(jié)果與分析

3.1 重要值

根據(jù)重要值大于5%選出各植被類型主要樹種。常綠闊葉林中主要樹種是青岡Cyclobalanopsis glauca、石櫟Lithocarpus glaber、苦櫧Castanopsis sclerophylla及木荷Schima superba等4種；松闊混交林的主要樹種分別是馬尾松Pinus massoniana、苦櫧、杉木Cunninghamia lanceolata、檵木Loropetalum chinense、石櫟等5種；松林的主要樹種分別是馬尾松、杉木、苦櫧、化香Platycarya strobilacea、楓香Liquidambar formosana等5種。

3.2 生態(tài)位重疊

從表2可以看出：在常綠闊葉林中，主要樹種之間生態(tài)位重疊程度最高的是木荷-石櫟(0.861 4)、其次是苦櫧-青岡(0.615 4)，最低的是苦櫧-木荷(0.377 2)；在松闊混交林中，重疊程度最高的是檵木-苦櫧(0.801 8)，其次是馬尾松-石櫟(0.753 1)，最低的是石櫟-杉木(0.340 9)；在松林中，重疊程度最高的是楓香-化香 (0.757 5)，其次是楓香-苦櫧 (0.720 1)，最低的是苦櫧-馬尾松 (0.442 7)。

表2 不同森林類型主要樹種生態(tài)位重疊數(shù)值矩陣Table 2 Niche overlap of dominant tree species of different forest types

3.3 空間分布格局

3.3.1 種內(nèi)空間分布格局 在所有研究樣地中，主要樹種在0～25.0 m的尺度內(nèi)均表現(xiàn)為顯著的聚集分布 (P＜0.05)(圖 1～3)。

圖1 常綠闊葉林主要樹種空間分布格局Figure 1 Spatial distribution pattern of dominant tree species in evergreen broad-leaved forest

圖2 松闊混交林主要樹種空間分布格局Figure 2 Spatial distribution pattern of dominant tree species in mason pine and broad-leaved mixed forest

圖3 松林主要樹種空間分布格局Figure 3 Spatial distribution pattern of dominant tree species in mason pine forest

3.3.2 種間空間分布格局 常綠闊葉林主要樹種空間分布格局特征如下 (圖 4)：青岡-石櫟在 0～25.0 m 尺度上均呈現(xiàn)負(fù)相關(guān)，其中6.0～10.0 m、15.0～17.0 m尺度上負(fù)相關(guān)達(dá)到顯著(P＜0.05)；青岡-木荷在0～25.0 m尺度上均為負(fù)相關(guān)，其中8.0～14.0 m尺度范圍內(nèi)達(dá)到顯著(P＜0.05)；石櫟-木荷在0～25.0 m尺度上均呈現(xiàn)顯著的正相關(guān)(P＜0.05)。其他種對的空間關(guān)聯(lián)性則以不顯著相關(guān)為主。

圖4 常綠闊葉林主要樹種間空間分布格局圖Figure 4 Spatial distribution pattern of dominant tree species in evergreen broad-leaved forest

松闊混交林主要樹種空間分布格局特征如下(圖5)：馬尾松-石櫟在0～25.0 m尺度上呈現(xiàn)正相關(guān)，10.0～25.0 m尺度內(nèi)達(dá)到顯著(P＜0.05)；苦櫧-石櫟在0～20.0 m尺度范圍內(nèi)為正相關(guān)，其中7.0～8.0 m尺度上達(dá)到顯著(P＜0.05)，而在20.0～25.0 m尺度范圍內(nèi)呈現(xiàn)不顯著負(fù)相關(guān)。其他種對的空間關(guān)聯(lián)性以不顯著相關(guān)或無關(guān)聯(lián)為主。

圖5 松闊混交林主要樹種間空間分布格局圖Figure 5 Spatial distribution pattern of dominant tree species in mason pine and broad-leaved mixed forest

松林主要樹種間空間分布格局特征如下(圖6)：馬尾松-化香在0～25.0 m尺度上為負(fù)相關(guān)，其中4.0 m時(shí)達(dá)到顯著(P＜0.05)；杉木-苦櫧在0～25.0 m尺度范圍內(nèi)為正相關(guān)，其中15.0～25.0 m尺度上達(dá)到顯著(P＜0.05)；杉木-化香在0～25.0 m尺度范圍內(nèi)為負(fù)相關(guān)，其中7.5～17.0 m尺度上達(dá)到顯著(P＜0.05)；化香-楓香在0～25.0 m尺度范圍內(nèi)表現(xiàn)為正相關(guān)，其中12.0～13.0 m尺度上達(dá)到顯著(P＜0.05)。除此之外其他種對的空間關(guān)聯(lián)性表現(xiàn)為不顯著相關(guān)或無關(guān)聯(lián)。

圖6 松林主要樹種間空間分布格局圖Figure 6 Spatial distribution pattern of dominant tree species in mason pine forest

3.4 種間聯(lián)結(jié)

3.4.1 總體聯(lián)結(jié)性 各類型次生林的方差比率(表3)結(jié)果表明：3種類型次生林主要樹種的總體關(guān)聯(lián)性均為不顯著正關(guān)聯(lián)。

表3 不同森林類型主要樹種的總體關(guān)聯(lián)性Table 3 Overall interspecific associations among dominant tree species of different forest types

3.4.2 種間聯(lián)結(jié)性檢驗(yàn) χ2檢驗(yàn)結(jié)果表明(表4)：樣地中所有類型次生林的主要樹種種對均呈不顯著的正聯(lián)結(jié) (χ2＜3.841， a d?bc＞0)。

表4 不同森林類型的種間聯(lián)結(jié) χ2統(tǒng)計(jì)量矩陣Table 4 Value of χ2 of different forest types

3.4.3 種間相關(guān)性分析 據(jù)圖7可知：在常綠闊葉林中，呈負(fù)相關(guān)的種對有苦櫧-木荷、木荷-青岡和青岡-石櫟等3對，其中木荷-青岡和青岡-石櫟種對的負(fù)相關(guān)達(dá)極顯著(P＜0.01)；呈正相關(guān)的有苦櫧-青岡、苦櫧-石櫟和木荷-石櫟等3對，其中木荷-石櫟種對的正相關(guān)達(dá)到極顯著(P＜0.01)。在松闊混交林中，呈現(xiàn)負(fù)相關(guān)的種對有檵木-杉木、檵木-石櫟2對；其余8對均呈正相關(guān)，其中馬尾松-石櫟的正相關(guān)達(dá)極顯著(P＜0.01)。在松林中，呈現(xiàn)負(fù)相關(guān)的種對有楓香-杉木、化香-馬尾松、化香-杉木、化香-苦櫧和馬尾松-苦櫧等5對，其中化香-杉木的負(fù)相關(guān)達(dá)顯著(P＜0.05)；其余5對呈正相關(guān)，其中楓香-苦櫧的正相關(guān)達(dá)顯著(P＜0.05)。

圖7 不同森林類型種間 Spearman 秩相關(guān)半矩陣圖Figure 7 Semi-matrix of Spearman’s rank correlation coefficients of dominant tree species of different forest types

4 討論

森林在向頂級群落演替的過程中，植物種內(nèi)空間分布格局一般會由聚集分布漸漸過渡為隨機(jī)分布[12]，而在本研究中，0～25.0 m尺度上不同類型次生林的主要樹種均呈顯著(P＜0.05)聚集分布，這說明樣地群落還未進(jìn)入演替成熟期。物種自身的生物學(xué)特性是造成群落演替過程中物種聚集分布的主要原因，即相同物種對環(huán)境等條件有著相似的需求[32]。樣地群落種間生態(tài)位重疊程度也能佐證其所處的演替階段，頂級群落中樹種間的生態(tài)位重疊程度一般處在較低的水平，即樹種間的競爭并不活躍，這是因?yàn)榻?jīng)過長時(shí)間演替后，群落趨于穩(wěn)定，內(nèi)部樹種間達(dá)到了一種相對平衡的狀態(tài)[33?34]，主要樹種間的生態(tài)位重疊程度仍普遍較高，群落整體呈正關(guān)聯(lián)，正相關(guān)種對仍有較高占比。這反映了群落中的主要樹種具有相似的環(huán)境需求和生態(tài)適應(yīng)性，由此造成了樣地中主要樹種的聚集分布。此外，種子的擴(kuò)散限制和生境異質(zhì)性也被認(rèn)為是影響物種分布的主要因素[20]，樣地中主要樹種的種子傳播方式以重力傳播為主，種子傳播距離有限，多聚集在母樹周圍，離母樹越遠(yuǎn)種子越少，這也導(dǎo)致了樹種的聚集分布。在松林中，馬尾松和杉木作為建群早期先鋒種最先入侵林地，在光競爭中占據(jù)優(yōu)勢，對光的有效利用是形成林分空間格局的主要決定因素[35]，因而馬尾松和杉木與群落中其他優(yōu)勢樹種多呈負(fù)空間關(guān)聯(lián)。楓香幼樹稍耐蔭，所以與群落中其他優(yōu)勢樹種多呈正空間關(guān)聯(lián)。在松闊混交林中，主要樹種間空間關(guān)聯(lián)顯著程度較低，這或許與群落中的針葉樹種同群落中其他較晚發(fā)育起來的常綠闊葉樹種產(chǎn)生了垂直分層現(xiàn)象有關(guān)，這種垂直分層現(xiàn)象能夠減輕群落內(nèi)的光競爭，進(jìn)而影響群落內(nèi)的樹種空間分布格局。馬尾松和杉木作為早期先鋒樹種入侵林地后，逐漸改善了立地條件，為其他樹種進(jìn)入群落創(chuàng)造了條件，使得松林能夠向松闊混交林方向演替。隨后由于馬尾松和杉木生物學(xué)特性在垂直結(jié)構(gòu)上與其他樹種產(chǎn)生了分層現(xiàn)象，加之密度制約和擴(kuò)散限制等因素[12]共同作用，群落中的常綠闊葉樹種逐漸獲得更多的環(huán)境資源并對針葉樹種和落葉闊葉樹種的幼苗更新產(chǎn)生負(fù)作用，松闊混交林逐漸向常綠闊葉林演替。這種演替機(jī)制與XIANG等[9]的研究結(jié)果一致，即在不發(fā)生干擾的情況下，次生林群落中的不耐蔭落葉樹種將逐漸被耐蔭的常綠闊葉樹種所取代。

通過比較不同森林類型主要樹種組成可以發(fā)現(xiàn)：3類森林群落的主要樹種組成存在差異，落葉闊葉樹種楓香和化香僅出現(xiàn)在松林主要樹種中，而松闊混交林中已經(jīng)沒有落葉樹種，但出現(xiàn)了常綠闊葉樹種石櫟和檵木。在常綠闊葉林中，主要樹種則由木荷及青岡、石櫟和苦櫧等3種殼斗科Fagaceae樹種組成。樹種組成的變化是群落演替的結(jié)果，森林經(jīng)營管理可根據(jù)此對松林和松闊混交林進(jìn)行適當(dāng)撫育，伐除部分針葉樹種及清理林內(nèi)枯立木，改善林分內(nèi)的光照條件，促進(jìn)常綠闊葉樹種的生長。在常綠闊葉林中，則可通過修枝來改善林內(nèi)光照條件，也可通過為森林土壤施加養(yǎng)分來改善目標(biāo)樹種的營養(yǎng)狀況，通過人為栽種苗木等手段來加快木荷、石櫟等目標(biāo)樹種的更新、生長和郁閉，促進(jìn)森林群落向頂級群落發(fā)展并促進(jìn)生物多樣性的恢復(fù)。

本研究通過方差比率法、χ2檢驗(yàn)和Spearman秩相關(guān)系數(shù)檢驗(yàn)，對不同類型次生林群落種間的聯(lián)結(jié)性和相關(guān)性進(jìn)行了分析。其中，χ2檢驗(yàn)結(jié)果是由種對的二元數(shù)據(jù)轉(zhuǎn)換計(jì)算而來，反映的是物種能否共存和共存的概率，不能表達(dá)低顯著度種對的內(nèi)在相關(guān)性情況和強(qiáng)度。Spearman秩相關(guān)系數(shù)屬定量檢驗(yàn)方法，能夠在一定程度上檢驗(yàn)種對數(shù)量關(guān)系上的變化，對種對間的相關(guān)性及其顯著性水平更為敏感。Spearman秩相關(guān)系數(shù)屬定量檢驗(yàn)與定性的方差比率法、χ2檢驗(yàn)結(jié)合使用能夠更全面地反映物種的種間關(guān)系[36]。

5 結(jié)論

通過測度常綠闊葉林、松闊混交林和松林群落中主要樹種的生態(tài)位重疊程度、空間分布格局和種間聯(lián)結(jié)分析可以發(fā)現(xiàn)：本研究各群落喬木層主要樹種在空間分布格局中獨(dú)立性相對較強(qiáng)，物種間雖存在比較相似的環(huán)境資源需求，但種間聯(lián)結(jié)關(guān)系比較松散，群落演替尚未進(jìn)入成熟期。此外也可以看出：樹種間的空間分布格局、聯(lián)結(jié)性與種群的生態(tài)位重疊之間存在密切關(guān)聯(lián)。一般情況下，群落內(nèi)優(yōu)勢樹種種間正聯(lián)結(jié)性越強(qiáng)，其生態(tài)位重疊程度越高，種間負(fù)聯(lián)結(jié)性越強(qiáng)，生態(tài)位重疊程度越低[37]。樣地群落的總體關(guān)聯(lián)性檢驗(yàn)結(jié)果顯示：各群落種間總體關(guān)聯(lián)性呈不顯著的正相關(guān)。χ2檢驗(yàn)發(fā)現(xiàn)：多數(shù)種對都呈不顯著的正聯(lián)結(jié)，這與各群落物種生態(tài)位重疊程度的分析結(jié)果基本吻合。各群落主要樹種間的Spearman秩相關(guān)分析結(jié)果與種間空間關(guān)聯(lián)性結(jié)果也有比較一致的表現(xiàn)，這也說明植物種間關(guān)系對植物種群空間分布格局的形成有重要意義。本研究可為進(jìn)一步揭示物種自身生物特性、環(huán)境條件及種間關(guān)系等綜合作用下的種群空間分布格局形成機(jī)制提供依據(jù)，但群落中植物種間關(guān)系及空間分布格局是處在動(dòng)態(tài)變化過程中的，本研究結(jié)果僅反映某一特定時(shí)間節(jié)點(diǎn)的群落狀態(tài)，有待長期的群落動(dòng)態(tài)監(jiān)測研究。