桫欏種群生態學研究概況

王新靚，盧 杰

(西藏農牧學院 高原生態研究所，西藏 林芝 860000； 西藏高原森林生態教育部重點實驗室，西藏 林芝 860000； 西藏林芝高山森林生態系統國家野外科學觀測研究站，西藏 林芝 860000； 西藏自治區高寒植被生態安全重點實驗室，西藏 林芝 860000)

桫欏是唯一擁有樹狀莖的古老的蕨類植物，與恐龍同一時期存活過，在第四紀冰川襲擊中存活至今。由于其生存環境受到破壞且繁殖條件較為苛刻，數量愈發稀少，被評定為國家二級保護植物。桫欏樹形優美，具有良好的觀賞性，是有名的觀賞植物。其葉與莖干含有豐富的黃酮類化合物，對心血管與內分泌疾病及癌癥的治療有良好的功效。我國對桫欏的研究可以追溯到20世紀80年代，主要針對桫欏的繁育與人工栽培[1-2]，后有學者研究了桫欏葉與莖干的化學成分[3-7]、形態結構[8-10]、種群與群落結構[11-15]等。

種群生態學是研究種群數量動態與環境相互作用關系的科學，對一種植物的種群生態學研究主要是從其種群結構與動態特征、空間分布格局、生態位與種間聯結性等多個方面進行。

1 桫欏種群的生境

大多數桫欏種群生長在海拔較低、降水豐沛、水熱條件良好的河邊溝谷及峭壁上。貴州赤水[13]，福建瓜溪[16]，四川峨眉山[15]、畫稿溪[12]與榮縣，云南威信[17]等長有桫欏的地方，都是中亞熱帶或亞熱帶濕潤氣候，降雨豐富，有些地方年降水量最高可達2 000 mm，年平均氣溫較高，且無霜期長，能夠為桫欏提供良好的水熱條件。桫欏的伴生種在各個地區的差別不大，毛竹、杉木與芭蕉是桫欏群落中常見的伴生種，在福建、貴州、四川與云南等桫欏分布區都有。有學者發現，毛竹常生長在桫欏幼苗附近，野芭蕉生長在桫欏成熟個體附近[13]。

2 種群結構及動態特征

種群在生態學中占據重要的地位[18]。種群結構中，垂直結構是種群的高度結構，而水平結構范圍比較廣泛，包括種群的徑級結構、年齡結構和冠幅結構等[19]。其中徑級結構是重要的劃分植株等級結構的指標，通常適用于原木的等級劃分，一般通過測量胸徑或基徑(DBH)進行劃分，不同種類的植物劃分等級標準也不同。如巴山榧樹(Torreyafargesii)的徑級劃分，便將其胸徑大小劃分為11個等級[20]，而高山松(Pinusdensata)的徑級結構劃分依據是參照其親本的徑級劃分情況并結合生活史特征，將其徑級結構劃分為12個等級[19]。種群年齡結構有增長型、穩定型、成熟型和衰退型，且種群的年齡結構與徑級、高度級關系密切。有些情況下，徑級或高度級可代替齡級來分析植物種群的年齡結構，但是對于桫欏來說，由于其沒有徑向生長特征，因此無法使用徑級替代齡級，故而通常選取莖干高度來對桫欏的齡級進行劃分。一般是在習水自然保護區樣地中，依據個體高度，將桫欏種群劃分為10個等級[21]。

對桫欏種群結構與動態分析的研究主要集中在四川習水與峨眉山、貴州赤水及海南省等地的自然保護區中，因為保護區內有大片的桫欏種群，能夠滿足種群結構的調查條件。種群結構的研究方法有很多，通常采用空間序列代替時間變化方法來分析桫欏的種群結構與動態特征。徐德靜等通過野外實地調查，采用繪制種群靜態生命表、存活曲線等方法來研究貴州習水自然保護區的桫欏種群生存狀況，計算種群齡級間動態指數(Vn)等指數，運用時間序列預測模型來預測該地區桫欏種群未來的發展趨勢。研究表明，桫欏種群的存活曲線呈Deevey-II型，年齡結構上中齡級占比最多(58.5%)，幼齡級次之(12.7%)，老齡級最少(1%)，整體呈不規則的金字塔形，動態為先衰退后增長；種群年齡結構數量動態變化指數V’pi(考慮外部因素干擾)的值趨近于零，說明該桫欏種群不穩定，極易受到外界環境因子的影響，計算該種群對完全隨機干擾的最大風險概率Pmax=12.5%，表明其對外界干擾敏感，抗干擾能力弱，種群發展不穩定[21]。謝春平等分析了海南島霸王嶺黑桫欏自然保護區中的黑桫欏種群特征，設置樣方，測量并記錄黑桫欏的株高、莖粗、冠幅、長勢等指標，對該群落內出現的其他種群進行登記，分析了優勢種、建群種。結果發現，黑桫欏種群存活曲線擬合是Deevey-II型[22]，種群年齡結構是中齡級多、老幼齡級少的金字塔型。宗秀虹等對赤水桫欏國家級自然保護區的桫欏種群進行研究，采取相似的方法對該樣地的桫欏進行種群結構及動態特征分析，比較發現，兩個不同地區的不同桫欏種群的齡級組成相同，都是中齡級>幼齡>老齡級，而生存分析曲線走勢也基本相同[23]。另外針對四川峨眉山自然保護區中的桫欏也有種群結構相關研究[15]。以上這幾個地區的桫欏種群存活曲線都是Deevey-II型(直線型)，與靜態生命表的存活率曲線相近。此外還有人使用生存分析函數，如生存函數、積累死亡率函數、死亡密度函數、危險率函數來輔助生命表的分析，為生命表的應用提供理論依據[24]。

3 空間分布格局及動態

植物種群的空間分布格局受植物自身的生物學特性、種內及種間關系、生境等的影響[25]。空間分布格局的研究方法有很多，如核密度、空間自相關、熱點分析等GIS空間分析法[26]，基于一些函數如Ripley’s K(r)函數的點格局分析法[27]，運用各種數學模型的指數分析法[28]等。

桫欏的空間分布格局常用分散度或聚集度指標來表示。貴州赤水桫欏自然保護區桫欏的空間分布格局研究使用計算分散度S2的方法，通過分散度的大小將種群空間格局分為均勻型、集群型和隨機型3類[13]。桫欏種群的空間分布格局受水熱條件的影響極大，當水熱因子較為均一時，桫欏的空間分布格局為隨機型，當水熱因子影響存在差異時，桫欏的空間分布格局為集群型。生長在不同地段的桫欏，其空間分布格局也不同，因為不同地段桫欏的生長環境與種群結構不同[16]。方差均值比與Morisita指數相結合，同樣可以判定桫欏種群的聚集程度[29]。龍文興等通過比較這兩個指數的大小來判定哪個樣地中桫欏的集群程度高。還有人將分形理論引入植物種群空間分布格局的研究中[29]。部分與整體以某種方式相似的形體稱為分形，早在2004年，有人以福建省瓜溪桫欏自然保護區的桫欏為研究對象，分析其種群生態學特性，在空間分布格局的研究中引入了分形模型，對種群格局分形維數進行計算，主要包括計盒維數、信息維數和關聯維數。分形模型對桫欏這種個體數量小、設置樣地困難的瀕危植物的種群生態學研究具有極大的幫助。周文娟以四川榮縣為研究地，通過公式計算得出幾個樣地桫欏種群的聚集強度，種群計盒維數擬合曲線得出種群對生態空間的占據能力，信息維數擬合直線得出種群對空間占據程度的強弱[30]。

4 生態位與種群聯接性

生態位這一概念于1917年首次被Grinnell[31]提出，隨后經Elton[32]進一步完善，給出了較為全面的定義。生態位是生態學研究中的一大熱點，生態位理論主要包括生態位寬度和生態位重疊，能夠解決具體的生態學問題[33]。一般來說，生態位寬度越大，表示該物種環境適應能力越強，具有較強的競爭力，往往在群落中處于優勢地位[23]。生態位重疊程度通常在一定程度上反映了物種對同等級資源的利用程度及空間配置關系。

自然界中，特定生態環境中存在的任何一個植物種群都是與其他植物種群相適應且共同存在的[20]，在植物群落的生長發育過程中，各個植物種群之間相互爭奪生存環境和資源，久而久之不同植物種群之間產生了聯系[21]，這就是種間聯結。種間聯結性的測定，對于珍稀瀕危植物的種質資源保護及其恢復具有重要的指導意義[22]。測定種間聯結性需要以2×2聯列表為基礎[23]。

對貴州赤水桫欏國家自然保護區桫欏群落的生態位與種間聯結性進行分析發現，生態位寬度大小為毛竹(Phyllostachysheterocycle)>桫欏>毛桐(Mallotusbarbatus)，說明毛竹數量多且范圍廣，在該群落中占據主導地位，是喬木層的建群種，對群落環境有著重要影響，且對環境資源的利用能力最強。桫欏的生態位寬度較毛竹小，但與除毛竹以外的其他種群相比，生態位是比較寬的，因此也能夠充分利用資源[33]。分析優勢種的總體聯結性計算發現，該區域群落物種間總體聯結性為不顯著負聯結，即群落優勢物種間總體聯結性弱。

貴州赤水與習水市具有典型的丹霞地貌，對該地貌的桫欏群落進行調查，分析了生態位寬度并計算了該群落的生態位相似比例。其中，生態位寬度使用兩種不同的指數來計算，得到各種群的生態位寬度大小不同，排列順序也不同，此時需結合兩種計算方式得到的結果進行分析，得出生態位相近的種群對環境利用能力相似，有些種群的生態位寬度為0，這是由于它們只出現在一個群落類型中[34]。

5 其他植物干擾對桫欏的影響

在其他物種對桫欏的影響研究中，最常見的是竹類，如毛竹、楠竹(Phyllostachyspubescens)、慈竹(Neosinocalamusaffinis)等。毛竹是禾本科剛竹屬，靠地下莖無性繁殖，因而繁殖速度很快，能夠迅速擴張，給群落其他種群的生存帶來了極大的威脅[35]。研究發現，毛竹在該群落中生態位寬度大于桫欏，二者生態位重疊程度高，存在很大的資源競爭[34]。毛竹對桫欏的感染程度分為輕度、中度和重度[36]。瞿歡歡等研究毛竹擴張對桫欏根系可塑性的影響，結果表明，重度干擾樣地的毛竹根系趨向于分布在肥沃的表層土壤[37]。為了更好地保護桫欏生長的原生境，必須采取一定的措施阻止毛竹的繼續擴張。楠竹和慈竹是我國西南地區發展經濟時引種的，具有生長快、材質好、用途多的特征，經濟價值很高。黃茹等研究了生長在四川樂山的西壩桫欏峽谷的楠竹與慈竹，發現這些竹類對桫欏的空間分布影響極大，方差均值比 (ID)、聚集度大小指數 (ICS)、Green指數(GI)等7項指數的測定計算表明，受到竹類入侵的桫欏種群空間分布呈聚集分布，未受到入侵的桫欏種群分布則比較規則[38]。

植物干擾會造成桫欏種群幼苗分布格局的改變，且不同植物引發的分布格局改變的類型也不同。

6 生態勢與生態場

除以上幾種桫欏種群生態學研究外，還有一項特別的研究——生態場(ecological field)。生態場理論是基于物理學的“場”理論，于1985年首次提出的[39]，是數學與物理等學科與生態學相互交叉、相互滲透的一種體現。天然狀態下植物的3種空間分布格局，即集群分布、隨機分布與均勻分布是生態場潛在作用規律性的表面特征[40]。

由于生態場這一概念提出的時間較短，相關研究成果較少[40]。有人運用生態場理論研究了日本落葉松林冠層光生態場[41]；還有人將生態場應用于東靈山遼東櫟林的研究中[39]；研究發現，生態場理論能夠更好地描述城市不同空間結構綠地的生態效應及其差異[42]。宋萍等運用桫欏種群生態場模型及種群生態場等勢線等工具，分析了桫欏種群生態場的作用強度與作用范圍，發現在個體數量差異不大的情況下，桫欏種群內生態場作用的強弱與植株個體的大小呈正相關關系[43]。隨后又對桫欏個體生態場梯度與生態勢進行了研究，發現桫欏個體生態場梯度的大小由多種因素決定，且植株個體小，生態場作用范圍就小，反之亦然[44]。

7 結論與展望

多數桫欏生活在海拔較低、降水豐富、氣候濕潤、水熱條件良好的河邊溝谷與峭壁，說明桫欏的生長對溫度、濕度與水分要求較高。桫欏種群的存活曲線大多是Deevey-II型，即直線型，與靜態生命表的存活率曲線相近。年齡結構都是中齡級>老齡級>幼齡級，整個種群結構呈現一種極不穩定的狀態，抗干擾能力弱。桫欏空間分布多為隨機型和集群型，生長在不同地段的桫欏空間分布格局也不同，這取決于桫欏種群的生境與種群結構。常與桫欏有生態位競爭的植物是竹類，如毛竹、楠竹等，這些對資源的掠奪對桫欏幼苗的影響很大，是造成桫欏幼苗存活率低的重要因素之一。干擾植物會改變桫欏幼苗的空間分布格局，且不同干擾植物對幼苗分布格局的影響不同。

種群是生態學研究、組成物種與物種進化的基本單元。桫欏處于瀕危的生存狀況，保護工作刻不容緩，由此提出以下幾點建議：一是加大宣傳力度，增加人們對桫欏的了解，加強保護意識。有些桫欏分布區的居民不認識、不了解桫欏，需要對其展開瀕危野生動植物科普宣傳與普法活動，使其正確認識桫欏的價值，自覺主動地參與到保護桫欏種群及整個生態環境的行動中。二是加強保護區內桫欏的保護。建立分布區內桫欏種群檔案庫，將桫欏的生長、發育、繁殖及死亡等數據記錄下來，對桫欏種群實行統一化、系統化的檔案管理，確切掌握桫欏種群的動態，制定相應且合理的保護措施。三是加大人工繁育桫欏的研究力度。桫欏繁育一直以來是個難題。程治英等對桫欏的種質保存技術與快速繁殖進行了研究[1]；王俊浩等為了恢復鼎湖山桫欏華南植物種子庫的作用，對桫欏進行了引種與孢子繁殖培育，成功培育出100余株桫欏幼苗，為植物的遷地保護與研究打下了良好的基礎[2]。四是適時介入桫欏幼苗的生長。桫欏瀕危的原因之一是桫欏幼苗成活率不高，這是因為其他植物如毛竹等對其生長干擾太大，競爭過于激烈，以致于幼苗沒有足夠的資源存活下來，因此需對桫欏生長區域進行定期巡查，針對桫欏面對的生存威脅及時進行干預，提高桫欏幼苗的存活率。