鳳陽山常綠闊葉林優勢種木荷的種群結構與動態

彭 輝，周紅敏，王宏翔，張弓喬

(1. 鳳陽山管理處，浙江 麗水 323000；2. 龍泉市林業科學研究院，浙江 龍泉 323700；3. 廣西大學林學學院，廣西 南寧 530004；4. 中國林業科學研究院林業研究所，北京 100091)

種群結構是植物種群的最基本特征，種群的結構與動態是種群生態學研究的核心之一，是揭示闡明種群更新與維持機制的關鍵[1]。分析種群現實特征并對其未來發展進行預測，能為該物種及其群落的科學保護及可持續經營提供理論依據。生命表、存活曲線及生存分析等是研究種群結構和動態的重要方法[2]，其中，編制靜態生命表可以展現種群的動態特征，利用生存分析函數可以展示種群生存規律[3]，結合時間序列分析等可以進一步探討種群生態學過程，預測其發展的趨勢[4]。

木荷（Schima superbaGardn. et Champ.）是山茶科（Theaceae）木荷屬（SchimaReinw.）的常綠喬木[5]。木荷習性偏陽，常與杉木（Cuninghamia lanceolata（lamb.） Hook.）等混交[6]，是我國亞熱帶山地最重要的鄉土闊葉造林樹種和重要的工業用材樹種之一。木荷是中亞熱帶常綠闊葉林的優勢樹種，也是常綠闊葉林恢復過程中的先鋒樹種[7]。在常綠闊葉林森林群落穩定和演替過程中扮演重要角色，但目前關于木荷天然種群在常綠闊葉林演替過程中的生存現狀和數量動態研究較少。鑒于此，本研究以鳳陽山常綠闊葉林木荷種群為研究對象，通過建立靜態生命表，繪制存活曲線，并結合種群數量動態量化研究及時間序列分析等方法深度剖析木荷種群的更新特征及生態對策，以期對木荷天然種群進行科學保護，進而為鳳陽山常綠闊葉林植被恢復與重建提供科學依據。

1 研究區概況

研究區位于武夷山系洞宮山山脈中段的浙江省龍泉市南的鳳陽山-百山祖國家級自然保護區內（119°06′～119°15′ E、27°46′～27°58′ N） ，面積151 km2。土壤以黃壤土為主，地帶性植被為中亞熱帶常綠闊葉林，植被垂直分布明顯。區內年平均氣溫12.3 ℃，年日照1 515.5 h，平均相對濕度80%，年蒸發量1 171.0 mm，年降水量2 438.2 mm。保護區內有種子植物164 科、659 屬、1 488 種，是南方紅豆杉（Taxus wallichianavar. mairei（Lemée H. Léveillé.） L. K. Fu et Nan Li）、 伯 樂 樹（Bretschneidera sinensisHemsl.）等6 種國家一級重點保護植物及多種珍稀瀕危植物重要分布區[8]。

在全面考察的基礎上，依照美國CTFS （the Center for Tropical Forest Science） 的森林生物多樣性監測規范，在鳳陽山典型常綠闊葉林內設立80 m × 100 m 的固定樣地（海拔1 300 m，坡度20°）。樣地內劃分為5 m × 5 m 的小樣方進行每木檢尺，調查胸徑（DBH）≥1 cm 的所有喬木樹種，掛牌并記錄每棵樹的樹種、位置坐標（x，y）、胸徑、樹高、冠幅等因子；灌草層記錄各物種的高度、蓋度、多度等指標。

2 研究方法

2.1 劃分木荷種群齡級

在相同的環境條件下，同一樹種的徑級和齡級對環境的反應規律具有一致性[9]。因此，在木荷種群年齡難以準確獲取的情況下，參照亞熱帶常綠闊葉林相關研究，結合鳳陽山0.8 hm2樣地內實際情況[10]，以徑級結構代替年齡結構來分析木荷種群的數量動態變化。按照DBH 的大小等級，把徑級I（DBH < 5 cm）作為幼樹，之后按DBH 大小以5 cm為一級劃分（上限排外），把徑級II （5 cm≤DBH <10 cm）作為小樹、徑級III -徑級Ⅴ（10 cm≤DBH < 25 cm）作為中樹、徑級Ⅵ-徑級Ⅸ （25 cm≤DBH < 45 cm）作為大樹，共分為4 個生長階段[11]。

2.2 編制靜態生命表

因現實中無法對某一種群的生活史進行跟蹤分析，特定時間生命表只能根據同一時間的種群個體徑級進行編制，進而反映世代重疊的年齡動態歷程[12]。在依據自然種群的不同齡級進行分析時，多因采取空間代替時間、徑級代替齡級等會導致死亡率出現負值的情形。大多數研究會對數據進行均滑處理，但因此反而忽略了種群真實的生態過程。為此，本研究采用各齡級實際個體數編制靜態生命表，并繪制種群結構圖。計算公式如下：

式中：ax是x齡級內的個體數，a0是ax的初始值，lx是在x齡級標準化后的存活個體，一般為1 000，dx和qx是從x齡級到x+ 1 齡級間標準化的死亡數和死亡率，Lx是x到x+ 1 齡級的存活個體，Tx是從x齡級起的個體總數，ex是x齡級的期望壽命值，kx是損失度（消失率），Sx是存活率。

2.3 擬合存活曲線

種群各齡級個體的存活狀況及過程可以通過存活曲線直觀反映[13]，一般以齡級為橫坐標，以實際存活個體數ax和存活個體對數值l nax為縱坐標，繪制種群存活曲線圖，采用Hett 等[14]提出的數學模型，依據Deevey 對存活曲線的分類，若指數方程Nx=N0e?bx擬合效果好，存活曲線符合DeeveyⅡ型描述，種群中各齡級個體的死亡率相同，為對角線型。若符合冪函數Nx=N0x?b方程，則存活曲線為DeeveyⅢ型[15]，種群中的幼齡個體具有較高的死亡率，隨著齡級增加死亡率降低且穩定。

2.4 量化種群動態

根據陳曉德[16]的方法，以種群兩個相鄰齡級個體的數量變化指數（Vn）、種群年齡結構數量動態變化指數（Vpi）及考慮外部干擾時種群年齡結構動態指數（V′pi）等客觀準確的量化分析木荷種群的動態變化。公式為：

式中：Vn是種群的n齡級到n+ 1 齡級個體數量變化動態，Vpi代表整個種群的數量動態變化指數，Sn和Sn+1分別是相鄰齡級的各組的個體數，k表示種群的齡級數量，max/min（）表示取其極大/極小值，Vn、Vpi取值在[?1,1]之間，取負值、零、正值時分別表示種群相鄰兩個齡級間（整個種群）是衰退、穩定、增長的動態變化。當有外部干擾時，用i表示整個種群年齡結構動態，其還與齡級數量（k）及各齡級個體數（Sn）相關，則：也是衡量隨機干擾對種群結構動態影響的指標，

2.5 時間序列模型

本研究中運用時間序列分析中的一次平均推移法，對木荷種群的年齡結構在未來2、4、6、8 個齡級發展趨向進行預測[17]，模型為：

式中：n為需預測的未來時間年限，t為齡級，為未來n年時t齡級的種群存活數量，是近期的n個觀測值在t時刻的平均值，稱為第n周期的移動平均，Xk為當前k齡級種群的存活數。

3 結果與分析

3.1 木荷種群年齡結構

本研究樣地中共有木荷692 株，種群密度為0.086 5 株?m?2，幼樹和小樹的比例較高（圖1），占整個種群的84.97%。小樹階段（徑級II）個體量急劇下降，占整體的20.8%，表明樣地中木荷種群存在激烈的種內和種間競爭，幼齡個體被大量淘汰，只有少數個體能進入下一生長階段。總體看，木荷年齡結構呈典型“倒J 形”種群，屬于增長型種群。

圖1 鳳陽山木荷種群的齡級結構Fig. 1 Size class structure ofSchima superbapopulation

3.2 種群生命表

對木荷種群靜態生命表的分析（表1）可知：木荷種群不同齡級個體存活數差別較大，特別是Ⅰ-Ⅱ齡級，存活個體下降迅速，在環境及競爭因素作用下淘汰了近68%的個體。期望壽命（）反映了齡級內種群個體的平均生存能力，木荷種群在Ⅲ-Ⅳ齡級期望壽命較高，死亡率較低，在Ⅶ齡級后個體生命期望值較低。整體上，種群期望壽命雖呈現出小幅波動，但期望值總體偏低，呈下降趨勢。木荷種群的死亡率（）二者有相同的變化趨勢，均在Ⅵ徑級出現高峰，死亡率為71.4% 消失率為1.253，而非呈單調遞減趨勢。

表1 木荷種群靜態生命表Table 1 Static life table ofSchima superbapopulation

3.3 種群存活曲線

存活曲線有3 種類型，Ⅰ 型是凸曲線，該類型早期死亡率低，但達到平均生理年齡時，會在短期內個體全部死亡，Ⅱ型、Ⅲ型分別為直線型和凹線型，Ⅱ型特征為種群各年齡段死亡率基本相同，Ⅲ型是前期種群個體死亡率高，一旦進入某一齡級死亡率低而穩定。因此，分別以實際存活個體數ax和其對數值 lnax做存活曲線圖（圖2）。由圖2可知：ax做的存活曲線圖在凹型和直線型之間，lnax做圖近似于直線型，因此，采用指數方程Nx=N0e?bx和冪函數Nx=N0x?b方程，進行擬合效果檢驗，結果見表2。對比4 個擬合方程可知：以ax和l nax為因變量擬合的指數模型中F值和R2均明顯高于相應的冪函數模型，因此，認為木荷種群的存活曲線更符合Deevey Ⅱ型描述，這表明木荷種群各齡級具有近乎相同的死亡率，與在靜態生命表中的分析一致。

圖2 鳳陽山木荷種群的存活曲線Fig. 2 survival curve ofSchima superbapopulation

表2 木荷種群存活曲線方程Table 2 The equation of survival curve forSchima superbapopulation

3.4 種群動態與時間序列

木荷種群相鄰齡級數量變化動態指數隨著齡級的增大表現不同（表3），種群動態指數大多為正值，說明種群多數時期屬于增長的結構動態，只是在V7時處于穩定階段，均大于0，說明即使受到外界隨機干擾時，種群仍能處于增長狀態，但對隨機干擾有較高的敏感度。

表3 木荷種群存動態變化指數Table 3 Dynamic index ofSchima superbapopulation

以各齡級的實際存活數作為基數，以時間序列的一次平均推移法對木荷種群未來2、4、6、8 個齡級后的存活個體數量進行預測，從表4 可知：在未來4 個齡級后，隨著時間推移各齡級個體數量均逐漸增加，且增加的趨勢增大，其中，在Ⅵ齡級，個體數量經過2、4、6 個齡級后分別增加了57.1%、257.1%、1 542.9%，而Ⅷ齡級經過4、6、8 個齡級后個體增加了2、7.5、42 倍。說明今后一段時期的演替過程中，木荷種群仍會在常綠闊葉林中占據優勢地位。

表4 木荷種群存動態變化指數Table 4 Time sequence ofSchima superbapopulation

4 討論

優勢物種是塑造群落結構的骨架[18]，對群落的外貌和類型具有重要調控作用。種群結構動態變化規律，反映了其對環境因子的適應[19]。木荷喜光，種子有種翅，在樣地中廣泛分布。在鳳陽山常綠闊葉林當前的演替階段中，木荷種群為典型的倒“J”形，屬于增長型種群，其中，齡級Ⅰ和齡級Ⅱ個體占種群數量的比例較大，說明木荷擁有大量的幼苗儲備，這與李立等[11]的研究結果一致。小樹階段（齡級Ⅱ）競爭最激烈，個體出現死亡高峰，該階段植物對其對環境資源，特別是光照、水分、土壤養分及生長空間的需求激增，是更新和演替過程中最敏感的階段[20]，受Jazen-Connell 效應制約，同種個體的競爭激烈，自疏明顯，結果為其他物種入侵和定居提供空間，從而提高、維持群落的多樣性。

靜態生命表及死亡率和消失率的分析結果顯示，種群的期望壽命隨著齡級增大呈下降趨勢但略有波動，死亡率和消失率有著相同的變化趨勢，在Ⅵ徑級出現高峰，可能是因為木荷種群已經開始進入生理死亡年齡階段。存活曲線更符合Deevey-Ⅱ型描述，說明外界的干擾可能一直存在，導致種群各年齡段具有相似的死亡率。相鄰齡級的動態指數均大于零（V7=0），為正值，說明木荷種群雖有輕微波動但一直處在增長動態變化中，種群動態指數和受外界干擾時的動態指數（Vpi和）均為正值，說明木荷種群為增長型種群，但對外界干擾敏感指數較高，暗示種群仍面臨較大的外界干擾。

時間序列是分析預測種群未來年齡結構的重要手段[21]，研究發現，在經歷2、4、6、8 個齡級后，木荷種群個體數量均呈現上升趨勢，特別是幼齡個體充足是，增長的潛力更大，說明充足的幼苗儲備是各齡級個體增加的前提保障，幼苗多則種群增長潛力大。

5 結論

木荷作為中亞熱帶常綠闊葉林的優勢種，對維持群落穩定具有重要意義，與清涼峰保護區木荷穩定型種群相比[22]，鳳陽山保護區的木荷種群與古田山木荷種群更相似，均為增長型種群，更新良好。鳳陽山木荷種群的年齡結構合理，但是對外界環境變化敏感。一旦環境遭到破壞，必將嚴重影響到木荷的正常生長發育。因此，必須鞏固加強現有的群落結構和生境保護，避免人為干擾造成的破壞。