鄒生根，帥曉鵬，劉年生，葉金山，劉曉華*

（1.永豐縣官山林場(chǎng)，江西 吉安 331506；2.江西省林業(yè)科學(xué)院）

江西國(guó)家級(jí)自然保護(hù)區(qū)種子植物相似性的聚類研究

鄒生根1，帥曉鵬1，劉年生1，葉金山2，劉曉華2*

（1.永豐縣官山林場(chǎng)，江西 吉安 331506；2.江西省林業(yè)科學(xué)院）

建立自然保護(hù)區(qū)是生物資源就地保護(hù)的有效手段。江西省國(guó)家級(jí)自然保護(hù)區(qū)數(shù)量較多，但缺乏它們之間的系統(tǒng)研究。本研究在校對(duì)及統(tǒng)計(jì)了分布在江西省12個(gè)國(guó)家級(jí)自然保護(hù)區(qū)野生種子植物狀況的基礎(chǔ)上，選用Rogers 和Tonimoto指數(shù)，通過(guò)聚類分析，得出了基于種子植物種級(jí)水平各保護(hù)區(qū)之間的關(guān)系，較屬級(jí)水平的聚類支序圖的聚類效果更佳，該省12個(gè)保護(hù)區(qū)聚在三大分支（贛北平原區(qū)、贛南山地區(qū)及贛中周邊山地）。另外，本研究還建立了江西10個(gè)森林類型國(guó)家級(jí)自然保護(hù)區(qū)野生種子植物物種—面積島嶼效應(yīng)模型，S=11.885×0.542A，較高的z值（0.542）反映了保護(hù)區(qū)種子植物生境異質(zhì)性高的特點(diǎn)。

江西??；國(guó)家級(jí)自然保護(hù)區(qū)；種子植物；聚類分析；物種—面積關(guān)系

江西省12個(gè)國(guó)家級(jí)自然保護(hù)區(qū)（官山、九嶺山、武夷山、陽(yáng)際峰、馬頭山、井岡山、九連山、齊云山、桃紅嶺、廬山、鄱陽(yáng)湖、南磯山），各保護(hù)區(qū)有關(guān)野生種子植物均做了詳實(shí)的記載［1］，為開展保護(hù)區(qū)野生種子植物聚類分析提供了數(shù)據(jù)基礎(chǔ)［2-13］。涂飛云等［1］利用江西省12個(gè)國(guó)家級(jí)自然保護(hù)區(qū)種子植物資料，從屬級(jí)水平，基于Sorenson指數(shù)進(jìn)行了研究。結(jié)果表明，某些國(guó)家級(jí)自然保護(hù)區(qū)（官山、九嶺山、井岡山、武夷山、馬頭山、陽(yáng)際峰）之間的關(guān)系沒有得到解決，并提出基于種子植物種級(jí)水平可能解決它們之間的關(guān)系。許多學(xué)者［14-16］以Rogers和Tonimoto指數(shù)對(duì)動(dòng)物地理分布進(jìn)行聚類分析，并取得較佳的地理區(qū)劃。張顯理等［16］研究了寧夏哺乳動(dòng)物地理分布，認(rèn)為Rogers和Tonimoto指數(shù)擬合最佳，優(yōu)于Sorenson指數(shù)2a/（2a+b+c）、Jaccard指數(shù)a/（a+b+c）的擬合效果。

涂飛云等［1］統(tǒng)計(jì)了江西省12個(gè)國(guó)家級(jí)自然保護(hù)區(qū)的面積及保護(hù)區(qū)種子植物物種數(shù)量狀況，然而，卻未對(duì)物種數(shù)—面積關(guān)系做進(jìn)一步研究。許多研究［17-20］表明，物種數(shù)與面積關(guān)系遵循島嶼物種—面積經(jīng)驗(yàn)?zāi)Ｐ蚐=c×zA，其中z值在一定程度上反映了各保護(hù)區(qū)島嶼化程度，也能更客觀反映區(qū)域植物地理狀況。

本研究目的∶一是基于種子植物屬、種水平，選用Rogers和Tonimoto指數(shù)，期望解決它們之間的關(guān)系，為區(qū)域植物地理區(qū)劃提供科學(xué)參考；二是通過(guò)對(duì)種子植物物種—面積關(guān)系的研究，進(jìn)一步認(rèn)識(shí)區(qū)域植物地理特征。

1 方法

參照涂飛云等［1］相關(guān)資料，對(duì)江西省12個(gè)國(guó)家級(jí)自然保護(hù)區(qū)野生種子植物屬、種進(jìn)行了核對(duì)與統(tǒng)計(jì)，詳見表1。

表1 12個(gè)國(guó)家級(jí)自然保護(hù)區(qū)種子植物屬、種和面積表

在統(tǒng)計(jì)的基礎(chǔ)上，按照張顯理等［16］的方法進(jìn)行二元屬性對(duì)屬、種編碼，即將有分布編碼為“1”，無(wú)分布的編碼為“0”，選用Rogers和Tonimoto指數(shù)，S=（a+d）/（a+2b+2c+d），其中 a，b，c，d分別是“1，1”，“1，0”，“0，1”，“0，0”匹配的物種數(shù)，建立二元數(shù)據(jù)指數(shù)矩陣。聚類分析均采用SPSS17.0軟件，聚類方法選用組內(nèi)聯(lián)接，區(qū)間度量標(biāo)準(zhǔn)使用Euc1idean距離。

假定各保護(hù)區(qū)是獨(dú)立的島嶼，遵循島嶼效應(yīng)物種-面積模型S=c×zA，式中S為物種數(shù)，A為面積，c、z為常數(shù)。研究認(rèn)為，物種—面積關(guān)系以對(duì)數(shù)函數(shù)擬合的線性方程1gS=1gc+z×1gA的效果好且應(yīng)用廣泛［21-23］。因此，本研究采用對(duì)數(shù)函數(shù)進(jìn)行擬合，并進(jìn)行顯著性檢驗(yàn)。本研究分兩個(gè)層面：一是江西省12個(gè)國(guó)家自然保護(hù)區(qū)種子植物物種與面積關(guān)系研究；二是因鄱陽(yáng)湖和南磯山是濕地類型自然保護(hù)區(qū)，生境類型差異大，故開展其余10個(gè)森林類型的國(guó)家級(jí)自然保護(hù)區(qū)物種—面積關(guān)系研究。相關(guān)性分析的顯著性檢驗(yàn)采用GraphPad Prism 5軟件。

2 結(jié)果

2.1 保護(hù)區(qū)種子植物屬的相似性聚類

從圖1看，鄱陽(yáng)湖與南磯山聚在一支，再同桃紅嶺聚在一支；九連山與齊云山聚在一起；陽(yáng)際峰、馬頭山聚在一支；官山、九嶺山聚一支，同這兩支關(guān)系最近的是武夷山，其次是廬山、井岡山；陽(yáng)際峰（或馬頭山）同官山（或九嶺山）之間的關(guān)系并不清楚。

圖1 基于種子植物屬的相似性自然保護(hù)區(qū)聚類支序圖

2.2 保護(hù)區(qū)種子植物種的相似性聚類

由圖2可見，保護(hù)區(qū)聚為三大支，分別是贛北平原區(qū)（鄱陽(yáng)湖、南磯山、桃紅嶺、廬山）、贛南山地區(qū)（九連山、齊云山）及贛中周邊山地（陽(yáng)際峰、馬頭山、官山、九嶺山、武夷山、井岡山）。保護(hù)區(qū)之間，鄱陽(yáng)湖和南磯山聚在一支，再和桃紅嶺聚在一起，其次是廬山；九連山和齊云山聚在一支；陽(yáng)際峰與馬頭山聚在一支，再同官山聚在一起，其次是九嶺山、武夷山、井岡山。

2.3 物種—面積關(guān)系研究

由圖3可知，江西省12個(gè)國(guó)家級(jí)自然保護(hù)區(qū)的種子植物物種與保護(hù)區(qū)面積之間不存在顯著的相關(guān)性（F=2.479，P=0.146）；而排除鄱陽(yáng)湖和南磯山后剩余的10個(gè)國(guó)家級(jí)自然保護(hù)區(qū)種子植物物種與保護(hù)區(qū)面積之間存在顯著的相關(guān)性（F=57.47，P＜0.000 1），這10個(gè)國(guó)家級(jí)自然保護(hù)區(qū)種子植物物種（對(duì)數(shù)）—面積（對(duì)數(shù)）擬合方程：1gS=0.542× 1gA+1.075，R2=0.878。由此推導(dǎo)出物種—面積方程：S=11.885×0.542A。

圖2 基于種子植物種的相似性自然保護(hù)區(qū)聚類支序圖

圖3 保護(hù)區(qū)種子植物物種—面積關(guān)系圖

3 討論

3.1 保護(hù)區(qū)種子植物屬、種相似性聚類

鄱陽(yáng)湖、南磯山聚在一支，齊云山和九連山聚在另一支，該結(jié)果與涂飛云等［1］研究結(jié)果一致。各保護(hù)區(qū)基于屬級(jí)水平相似性研究，其結(jié)果顯示陽(yáng)際峰（或馬頭山）和官山（或九嶺山）之間的關(guān)系并不清楚；地處贛北平原的廬山與地處贛西南的井岡山相似性高，不符合它們的地理分布和氣候特點(diǎn)。物種分布反映物種對(duì)特定地理區(qū)域的適應(yīng)性。作者認(rèn)為相比屬級(jí)水平，種級(jí)水平更能反映區(qū)域地理特征及關(guān)系?；诜N子植物相似性聚類關(guān)系得到全面的呈現(xiàn)，可劃為3大地理區(qū)域，贛北平原區(qū)、贛中周邊山地、贛南山地區(qū)。其中分布在贛中部及周邊山地區(qū)的保護(hù)區(qū)種子植物相似性關(guān)系較為復(fù)雜，如分布在武夷山脈的武夷山，與九嶺山、官山的相似性高，而非同處武夷山脈且距離較近的馬頭山、陽(yáng)際峰［1］。原因可能是武夷山保護(hù)區(qū)面積大，海拔跨度大，隨著海拔的升高，表現(xiàn)出高緯度種子植物特點(diǎn)［1］，與緯度偏北的九嶺山、官山的較為相似。李言闊等［24］對(duì)江西13個(gè)保護(hù)地的兩棲類聚類分析表明，地處武夷山脈的馬頭山與九嶺山脈的九嶺山的相似性高，而非距離較近的陽(yáng)際峰，這可能也是個(gè)印證，但需要進(jìn)一步的研究。

3.2 物種—面積相關(guān)性

研究認(rèn)為z值能反映研究類群的所在島嶼的生境異質(zhì)性［18，25］。有學(xué)者［26-27］針對(duì)z值確定了3個(gè)不同的等級(jí)∶生物地理區(qū)之間（0.50～1.00）、生物地理區(qū)內(nèi)（0.20～0.49）和島嶼（或隔離生境斑塊）（0.10～0.19）。本研究中的z值（0.542）處于生物地理區(qū)之間的水平，該值高于Rosenzweig［27］記載的許多研究有關(guān)z值的范圍（0.10～0.40），并且顯著大于研究區(qū)域?yàn)閸u嶼時(shí)的z值（0.35）［27］。影響z值大小的因素有生物類群［28］、擴(kuò)散限度［28-29］、生境異質(zhì)性［25］及其相互之間的共同效應(yīng)等因素［28-29］。Aranda等［28］研究表明，種子植物取得的物種—面積模型中的z值高于蕨類植物、苔蘚植物。Shen 等［29］認(rèn)為，研究類群的擴(kuò)散限度及其生境異質(zhì)性共同影響著物種—面積關(guān)系。高的z值（0.542）可能既反映了種子植物的擴(kuò)散特點(diǎn)，又反映了保護(hù)區(qū)的生境異質(zhì)性高的特點(diǎn)。

物種—面積關(guān)系可以估算區(qū)域的物種數(shù)量、評(píng)價(jià)區(qū)域物種多樣性的喪失［30-31］。本研究利用10個(gè)森林類型國(guó)家級(jí)自然保護(hù)區(qū)種子植物物種數(shù)—面積數(shù)據(jù)，得到關(guān)系方程，并由此推算出江西省森林類型自然保護(hù)區(qū)的種子植物多樣性狀況，可為保護(hù)區(qū)種子植物多樣性研究提供參考。

［1］涂飛云，田徑，章蜜，等.江西國(guó)家自然保護(hù)區(qū)種子植物區(qū)系及相似性比較研究［J］. 江西林業(yè)科技，2014，24（2）：21-26.

［2］葛剛，陳少鳳，萬(wàn)文豪，等.裸子植物GYMNOSPERMAE、被子植物ANGIOSPERMAE［M］ // 劉信中，吳和平. 江西官山自然保護(hù)區(qū)科學(xué)考察與研究. 北京： 中國(guó)林業(yè)出版社，2005.

［3］候?qū)W良，陳少鳳，楊清培，等.植物資源［M］//李振基，吳小平，陳小麟，等. 江西九嶺山自然保護(hù)區(qū)綜合科學(xué)考察報(bào)告. 北京： 科學(xué)出版社，2009.

［4］劉信中，姚振生，陳鳳彬，等.江西武夷山自然保護(hù)區(qū)種子植物名錄［M］ // 劉信中，方福生. 江西武夷山自然保護(hù)區(qū)科學(xué)考察集. 北京： 中國(guó)林業(yè)出版社，2001.

［5］郭英榮，魯順保，張艷杰，等.江西陽(yáng)際峰自然保護(hù)區(qū)生物多樣性編目［M］ // 郭英榮，江波，王英永. 江西陽(yáng)際峰自然保護(hù)區(qū)綜合科學(xué)考察報(bào)告. 北京： 科學(xué)出版社，2010.

［6］劉仁林，譚策銘，裘利洪.種子植物名錄［M］//劉信中，傅清. 江西馬頭山自然保護(hù)區(qū)科學(xué)考察與稀有植物群落研究. 北京： 中國(guó)林業(yè)出版社，2006.

［7］劉仁林.井岡山自然保護(hù)區(qū)種子植物名錄［M］.北京：中國(guó)科學(xué)技術(shù)出版社，2013.

［8］姚振生，劉信中，賴學(xué)文，等.九連山自然保護(hù)區(qū)種子植物名錄［M］ // 劉信中，肖忠優(yōu)，馬建華. 江西九連山自然保護(hù)區(qū)科學(xué)考察與森林生態(tài)系統(tǒng)研究. 北京： 中國(guó)出版社，2002.

［9］裘利洪，劉仁林，季春峰，等.種子植物名錄［M］//劉小明，郭英榮，劉仁林. 江西齊云山自然保護(hù)區(qū)綜合科學(xué)考察集. 北京： 中國(guó)林業(yè)出版社，2010.

［10］覃海寧，于勝祥，武建勇，等.桃紅嶺自然保護(hù)區(qū)種子植物名錄［M］ // 蔣志剛. 江西桃紅嶺梅花鹿國(guó)家級(jí)自然保護(hù)區(qū)生物多樣性研究. 北京： 清華大學(xué)出版社，2009.

［11］葛剛，劉以珍，萬(wàn)文豪，等.植物區(qū)系（種子植物區(qū)系）［M］// 劉信中，王瑯. 江西省廬山自然保護(hù)區(qū)生物多樣性考察與研究. 北京： 科學(xué)出版社，2010： 539-567.

［12］萬(wàn)文豪，紀(jì)偉濤，吳建東，等.植物名錄［M］//吳英豪，紀(jì)偉濤. 江西鄱陽(yáng)湖國(guó)家級(jí)自然保護(hù)區(qū)研究. 北京： 中國(guó)林業(yè)出版社，2002.

［13］葛剛，陳少鳳，郁文彬，等.植物資源［M］//劉信中，樊三寶，胡斌華. 江西南磯山濕地自然保護(hù)區(qū)綜合科學(xué)考察集. 北京： 中國(guó)林業(yè)出版社，2006.

［14］徐寧，魏剛，李德俊.貴州各動(dòng)物地理省兩棲動(dòng)物分布聚類的初步探討［J］. 動(dòng)物學(xué)雜志，1992，27（3）： 16.

［15］楊再學(xué).貴州省嚙齒動(dòng)物地理分布的聚類分析［J］.四川動(dòng)物，1996，15（3）： 119-120.

［16］張顯理，邵月芬，張大治.寧夏哺乳動(dòng)物地理分布的聚類分析 ［J］.寧夏大學(xué)學(xué)報(bào)（自然科學(xué)版），1999，20（3）：285-288.

［17］Preston FW.The canonical distribution of commonness and rarity：Part I ［J］. Ecology，1962，43： 185-215.

［18］Macarthur RH，Wilson EO.An equilibrium theory of insular zoogeography［J］. Evolution，1963，17（4）： 373-387.

［19］May RM.Patterns of species abundance and distribution［M］//Cody MC，Diamond JM（eds.）. Ecology and evolution of communities. Cambridge： Belknap Press，1975.

［20］Harte J，Green J.Self-similarity in the distribution and abundance of species［J］. Science，1999，284（5412）： 334-336.

［21］Tj rve，E.Shapes and functions of species-area curves（II）：a review of new models parameterizations ［J］. Journal of biogeography，2009，36（8）： 1435-1445.

［22］孫雀，盧劍波，張鳳鳳，等.物種多樣性與島嶼面積的關(guān)系［J］. 生態(tài)學(xué)報(bào)，2009，29（5）： 2195-2202.

［23］要元媛，閆明，畢潤(rùn)成.山西霍山植物群落種-面積曲線與物種多樣性的關(guān)系［J］. 生態(tài)學(xué)雜志，2013，32（1）： 39-44.

［24］李言闊，單繼紅，龔瑜.江西省兩棲類動(dòng)物多樣性與地理區(qū)劃［J］. 動(dòng)物學(xué)雜志，2013，48（6）： 919-925.

［25］Desmet P，Cowling R.Using the species-area relationship to set baseline targets for conservation［J］. Ecology and Society，2004，9（2）： 301-303.

［26］Brown JH，Lomolino MV.Biogeography［M］.Sunderland，MA： Sinauer Associates，1998.

［27］Rosenzweig M.Species diversity in space and time［M］. Cambridge，UK： Cambridge University Press，1995.

［28］Aranda SC，Gabriel R，Borges PAV，et al.How do different dispersal modes shape the species-area relationship？ Evidence for between-group coherence in the Macaronesian flora［J］. Global Ecology and Biogeography，2013，22（4）： 483-493.

［29］Shen GC，Yu MJ，Hu XS，et al.Species-area relationships explained by the joint effects of dispersal limitation and habitat heterogeneity［J］. Ecology，2009，90（11）： 3033-3041.

［30］Malcolm JR，Liu CR，Neilson RP，et al.Global warming and extinctions of endemic species from biodiversity hotspots ［J］. Conservation Biology，2006，20（2）： 538-548.

［31］唐志堯，喬秀娟，方精云.生物群落的種—面積關(guān)系［J］.生物多樣性，2009，17（6）： 549-559.

（責(zé)任編輯：楊婷婷）

Clustering Study on the Similarity of Seed Plants in National Nature Reserves in Jiangxi Province

ZOU Shenggen1，SHUAI Xiaopeng1，Liu Niansheng1，YE Jinshan2，LIU Xiaohua2*
（1.Guanshan Forest Farm of Yongfeng County，Ji’an 221506，Jiangxi，China；2.Jiangxi Academy of Forestry）

Establishing nature reserves is an effective means of in-situ conservation of biological resources.There are a dozen of national nature reserves in Jiangxi Province which were lack of systematic study.In this paper，based on check and statistics of the wild seed plants in the 12 national nature reserves distributed in Jiangxi Province，a clustering analysis was made with Rogers&Tonimoto indexes.The results showed the clustering relationship among the reserves based on the species-level of seed plants and the clustering effect is better than that of the genus-level clustering cladogram.The 12 reserves cluster in 3 classes（North Jiangxi Plain Region，South Jiangxi Mountainous Region&Mountainous Region Surrounding Central Jiangxi）.Besides，the species-area island effect models of the wild seed plants were established in the 10 forest-type national nature reserves in the province，S=11.885×0.542A，with a relatively high z value（0.542），indicating characteristics of high heterogeneity of seed plant habitats in the reserves.

Jiangxi Province；National nature reserve；Seed plants；Clustering analysis；Species-area relationship

S718.3

A

2095-0152（2016）03-0010-04

2016-03-30

2016-04-14

鄒生根，男，林業(yè)工程師，主要從事營(yíng)造林工作。

劉曉華（1980-），女，從事野生動(dòng)植物研究。E-mail：410623173@qq.com