王魏根

?

長(zhǎng)江中下游湖泊螺類beta多樣性分析

王魏根

阜陽師范學(xué)院生物與食品工程學(xué)院，阜陽 236037

采用兩種方法對(duì)長(zhǎng)江中下游19個(gè)湖泊螺類beta多樣性進(jìn)行了分析。首先采用Legendre方法分析單個(gè)螺類物種(SCBD)和單個(gè)湖泊(LCBD)在β多樣性中的貢獻(xiàn)。結(jié)果表明分布范圍居中的螺類, 如大沼螺(), SCBD值大, 分布范圍窄的和分布范圍廣的螺類SCBD值均會(huì)降低; 花馬湖LCBD值最大, 單個(gè)湖泊的LCBD值和其中分布的螺類物種數(shù)沒有顯著的相關(guān)性。二是采用Baselga方法分析螺類β多樣性中的周轉(zhuǎn)和嵌套成分。多位點(diǎn)計(jì)算結(jié)果表明長(zhǎng)江中下游湖泊中螺類β多樣性主要原因是物種在空間上的周轉(zhuǎn)。配對(duì)的分析結(jié)果表明不同湖泊對(duì)間的β多樣性周轉(zhuǎn)和嵌套組成比例不同, 并且存在著完全的周轉(zhuǎn)和完全的嵌套格局。

湖泊; 螺類; beta多樣性; 周轉(zhuǎn); 嵌套

1 前言

比較不同群落中物種組成及其變化是群落生態(tài)學(xué)研究的重要內(nèi)容, Whittaker(1960)首次將物種多樣性分成3個(gè)層次: α多樣性(局域物種的多樣性)、γ多樣性(區(qū)域尺度的物種多樣性)和β多樣性(群落組成沿樣帶或某一環(huán)境梯度、格局變化的測(cè)度)[1], 但反映不同采樣單元間群落組成非方向性、無梯度的變化也是描述β多樣性合適的方法[2]。由于Whittaker最初對(duì)beta多樣性定義內(nèi)容不很精確, 所以隨后的研究者根據(jù)自己研究目的不同定義了很多β多樣性內(nèi)容及其度量方法[3–4]。

Legendre提出一種通過計(jì)算由位點(diǎn)-物種構(gòu)成的群落組成資料矩陣的方差來估計(jì)β多樣性的方法[5]。總β多樣性(BD)等于總的離均差平方(SS)和除以其自由度(-1), 也就是方差。為了了解單一物種或單一位點(diǎn)在β多樣性中的貢獻(xiàn), 總β多樣性可以繼續(xù)分解成單一物種對(duì)總β多樣性的貢獻(xiàn)度SCBD (Species Contribution to Beta Diversity)和單個(gè)位點(diǎn)對(duì)總β多樣性的貢獻(xiàn)度LCBD(Local Contribu-tion to Beta Diversity)[2]。該方法既可以利用群落頻度數(shù)據(jù)和豐度數(shù)據(jù)直接計(jì)算, 也可以通過相異性矩陣來計(jì)算, 但相異性矩陣必須是Euclidean性質(zhì)距離, 只能計(jì)算LCBD[2]。某位點(diǎn)LCBD值大則表示該樣點(diǎn)群落組成和區(qū)域內(nèi)其它位點(diǎn)組成差異大, 也說明該位點(diǎn)具有較高的保護(hù)價(jià)值, 相反則說明該位點(diǎn)物種貧乏或者環(huán)境被破壞需要修復(fù)。某物種SCBD值大小則反映該物種對(duì)該區(qū)域β多樣性影響的相對(duì)重要性。

根據(jù)群落間β多樣性來源的相關(guān)理論, Baselga[6]提出β多樣性可分解成物種在空間上周轉(zhuǎn)和嵌套兩個(gè)部分, 分別代表在空間上的物種替代和物種得失兩個(gè)不同的過程, 并基于S?rensen相異性指數(shù), 提出了加性的分解方法: βsor=βsim+βnes, 即總β多樣性(βsor)可以分解成物種替換引起的β多樣性(βsim)和嵌套結(jié)構(gòu)導(dǎo)致的β多樣性(βnes)。該方法既可以利用群落組成配對(duì)的相異性矩陣資料計(jì)算, 也可以利用多位點(diǎn)相異性矩陣資料計(jì)算, 多位點(diǎn)方法總β多樣性、周轉(zhuǎn)和嵌套部分分別用大寫英文下標(biāo)βSOR、βSIM、和βNES表示[6]。在保護(hù)生物學(xué)中, 總β多樣性中周轉(zhuǎn)和嵌套組分相對(duì)比例大小分別暗示著不同的保護(hù)措施, 嵌套為主導(dǎo)區(qū)域需要優(yōu)先保護(hù)物種多樣性豐富度高的位點(diǎn), 周轉(zhuǎn)為主導(dǎo)區(qū)域則需要同時(shí)保護(hù)多個(gè)位點(diǎn)[16]。

長(zhǎng)江是中國(guó)最長(zhǎng)的河流, 其中下游地區(qū)分布著眾多湖泊, 湖泊總面積達(dá)15770 km2, 構(gòu)成了一個(gè)復(fù)雜的江湖復(fù)合生態(tài)系統(tǒng), 長(zhǎng)期演化中形成了獨(dú)特而多樣的生物類群, 被世界自然基金會(huì)列為全球200個(gè)需要優(yōu)先保護(hù)的生態(tài)區(qū)之一[7]。淡水螺類是底棲生物群落的重要組成部分, 是水生食物鏈的基本環(huán)節(jié)之一, 和人類生活關(guān)系也很密切。長(zhǎng)江中下游是我國(guó)淡水螺類資源分布集中的地區(qū), 特有種多。自20世紀(jì)50年代以來, 進(jìn)行過幾次零散的螺類資源調(diào)查, 但這些調(diào)查主要是了解各湖泊的區(qū)系組成[8]。近期研究表明, 江湖阻隔、過度利用和水體污染等人為因素已造成長(zhǎng)江中下游湖泊的螺類多樣性明顯下降, 特別是受阻隔湖泊的螺類物種較少。根據(jù)物種保護(hù)價(jià)值指數(shù),有5種處于受威脅等級(jí), 5種處于近危等級(jí)[9]。

β多樣性分解方法可以通過度量α多樣性變化的方法在一定尺度上來預(yù)測(cè)γ多樣性變化, 通過了解保持β多樣性背后的機(jī)制可以為較大區(qū)域內(nèi)生物多樣性保護(hù)提供重要啟示。生物多樣性保護(hù)實(shí)踐中, 對(duì)研究區(qū)域內(nèi)β多樣性進(jìn)行定量分析, 有利于量化生物多樣性損失并確定保護(hù)區(qū)的位置, 對(duì)生物入侵物種管理等都有重要意義[10]。因此我們采用2種β多樣性分解方法, 對(duì)長(zhǎng)江中下游地區(qū)湖泊淡水螺類群落組成資料進(jìn)行β多樣性分析, 探討該地區(qū)螺類β多樣變化的原因(周轉(zhuǎn)和嵌套)以及各個(gè)物種(SCBD)、各個(gè)湖泊(LCBD)在整個(gè)β多樣性中的貢獻(xiàn), 為該地區(qū)淡水螺類物種多樣性保護(hù)和資源合理利用提供科學(xué)依據(jù)。

2 材料和方法

2.1 研究區(qū)概況

長(zhǎng)江中下游地區(qū)屬于亞熱帶季風(fēng)氣候, 年平均氣溫13—20 ℃, 年平均降水量800—1600mm[7]。本研究共編輯整理了長(zhǎng)江中下游19個(gè)湖泊的螺類名錄資料進(jìn)行分析(見表1), 湖泊按面積由大到小順序編號(hào)為1—19, 依次是鄱陽湖(江西)[17–19]、洞庭湖(湖南)[17,20]、太湖(江蘇)[17,21–22]、巢湖(安徽)[23]、洪湖(湖北)[17,24]、梁子湖(湖北)[17]、菜子湖(安徽)[25]、南漪湖(安徽)[26]、大通湖(湖南)[27]、淀山湖(上海)[28]、牛山湖(湖北)[17]、保安湖(湖北)[17]、東湖(湖北)[17,29]、花馬湖(湖北)[30]、橋墩湖(湖北)、青菱湖(湖北)、大黃家湖(湖北)、后湖(湖北)和扁擔(dān)塘(湖北)[17]。這些湖泊處在長(zhǎng)江流域最上游的是洞庭湖(28°18′8—30°12′ N, 110°24′—113°6′E), 處在長(zhǎng)江流域最下游的是淀山湖(31°04′—31°12′ N, 120°54′—121°01′ E)(圖1)。

圖1 湖泊位置示意圖

表1 長(zhǎng)江中下游湖泊特征及螺類物種數(shù)

在這19個(gè)湖泊中, 共記錄了螺類名錄79種(表2), 隸屬于9個(gè)科, 種類最多的為田螺科25種, 最少的膀胱螺科1種[17-30]。這些螺類分布范圍差異較大,有的廣泛分布, 如田螺科的銅銹環(huán)棱螺(a)在19個(gè)湖泊中都有分布, 豆螺科的長(zhǎng)角涵螺()在18個(gè)湖泊中有分布。而分布稀少的僅在單個(gè)湖泊中出現(xiàn)記錄。在各個(gè)湖泊中, 以鄱陽湖記錄螺類種數(shù)最多為48種, 大黃家湖和后湖記錄螺類種數(shù)最少為5種。

表2 長(zhǎng)江中下游湖泊螺類組成概況

2.2 數(shù)據(jù)整理

為了盡可能全面了解長(zhǎng)江中下游湖泊中螺類物種名錄及其分布狀況, 我們依據(jù)多年已經(jīng)公開發(fā)表文獻(xiàn)的調(diào)查記錄進(jìn)行編輯整理成數(shù)據(jù)集。數(shù)據(jù)記錄成為有-無類型的二元數(shù)據(jù)表, 即某種螺類在該湖泊出現(xiàn)就記錄為“1”, 沒有出現(xiàn)則標(biāo)記為“0”。所有用于統(tǒng)計(jì)分析的螺類分類單元均為“種”。對(duì)于不同文獻(xiàn)中記錄的同一物種中文名有可能不同, 則按照其拉丁名統(tǒng)一為一個(gè)物種。有的文獻(xiàn)中記錄的物種沒有鑒定到種名, 這些物種則不包括在這次分析數(shù)據(jù)中。

2.3 β多樣性分解方法

設(shè)SS和SS分別為物種的方差和第個(gè)湖泊的方差, 則SCBD和LCBD計(jì)算公式分別為:

SCBD=SS/Total; LCBD=SS/Total。

采用Baselga[6]的加性分解方法, 利用S?rensen相異性指數(shù)度量β多樣性及其周轉(zhuǎn)和嵌套成分, 配對(duì)的方法計(jì)算sor、sim和nes公式如下:

式中是兩個(gè)湖泊共有的物種數(shù),是只在第一個(gè)湖泊存在的物種數(shù),是只在第二個(gè)湖泊存在的物種數(shù), 總β多樣性(βsor)通過S?rensen相異性指數(shù)計(jì)算而來, 周轉(zhuǎn)成分(βsim)由Simpson相異性指數(shù)計(jì)算得到, 嵌套部分(βnes)則可由前兩者之差計(jì)算得到。多位點(diǎn)方法計(jì)算總β多樣性(βSOR)、周轉(zhuǎn)(βSIM)和嵌套(βNES)組分計(jì)算公式請(qǐng)參閱文獻(xiàn)[6]。

2.4 數(shù)據(jù)分析

單一物種對(duì)總β多樣性的貢獻(xiàn)度(SCBD)和單個(gè)樣方對(duì)總β多樣性的貢獻(xiàn)度(LCBD)分析過程用R語言工具包adespatial中的beta.div函數(shù)完成[11]。即計(jì)算湖泊螺類組成矩陣總方差(BD)作為總β多樣性值, 然后再分解為單個(gè)物種和單個(gè)湖泊對(duì)總β的貢獻(xiàn)度。計(jì)算方差前, 螺類組成數(shù)據(jù)先進(jìn)行“hellinger”轉(zhuǎn)化, 使相異性矩陣具有Euclidean距離特征。LCBD值的顯著性檢驗(yàn)采用置換檢驗(yàn), 置換次數(shù)為999。為了探討SCBD值和螺類物種分布范圍的關(guān)系, 我們用SCBD值和物種分布湖泊數(shù)量進(jìn)行相關(guān)分析。為了探討LCBD值與湖泊面積、深度及其中螺類物種數(shù)之間的關(guān)系, 我們用LCBD值與各個(gè)湖泊面積、水深以及螺類物種數(shù)進(jìn)行相關(guān)分析。相關(guān)分析用R語言工具包betareg完成[12]。

β多樣性及其周轉(zhuǎn)和嵌套組分的分解過程采用R語言工具包betapart完成[13]。為了可視化不同湖泊對(duì)之間β多樣性關(guān)系, 我們利用配對(duì)計(jì)算的相異性矩陣進(jìn)行非度量多維尺度分析(NMDS)和聚類分析(單連接聚合聚類法)。

以上兩種β多樣性分析過程均在統(tǒng)計(jì)軟件R語言環(huán)境中進(jìn)行[14], 除了以上所述工具包外, 還利用了beanplot[15]、vegan工具包[16]和R基礎(chǔ)工具包。

3 結(jié)果與分析

3.1 單一螺類SCBD和單一湖泊LCBD

根據(jù)Legendre[5]方法計(jì)算結(jié)果表明, 長(zhǎng)江中下游19個(gè)湖泊螺類SS即總的離均差平方和為9.8039,BD即總β多樣性值為0.5447。

SCBD最小值為0.0020, 由瘦沼螺()、紅帶短溝蜷()、河環(huán)棱螺(貢獻(xiàn), 這三種螺類均只在鄱陽湖有分布。最大值為0.5133, 由大沼螺()貢獻(xiàn), 該螺在12個(gè)湖泊有分布。SCBD平均值為0.0127, 有26個(gè)物種超過平均值。Beta相關(guān)性分析表明, SCBD值和該物種分布的湖泊數(shù)量顯著相關(guān)(<0.001,2=0.6243)(圖2、3)。

LCBD值最小為0.0349, 由淀山湖貢獻(xiàn)。最大值為0.0711, 由花馬湖貢獻(xiàn)。各個(gè)湖泊平均貢獻(xiàn)度為0.0526, 有9個(gè)湖泊LCBD貢獻(xiàn)度超過平均(按順序依次是CZH、DHJH、CHH、DATH、HOUH、DTH、DH、NSHH、HMH)(圖4)。但置換檢驗(yàn)表明, 只有牛山湖(=0.047)和花馬湖LCBD(=0.013)統(tǒng)計(jì)概率達(dá)到顯著水平。相關(guān)分析結(jié)果表明, LCBD與湖泊面積、平均水深和螺類物種數(shù)三者關(guān)系均不顯著(>0.14)。

圖2 SCBD值分布圖

圖3 SCBD值與分布湖泊數(shù)關(guān)系

3.2 周轉(zhuǎn)(βsim)和嵌套(βnes)成分

根據(jù)Baselga[6]的分解方法, 首先采用多位點(diǎn)方法計(jì)算, βSOR、βSIM和βNES值分別為0.8779、0.7143和0.1636, 周轉(zhuǎn)和嵌套成分在總β多樣性中比例分別是81.4%和18.6%, 表明長(zhǎng)江中下游湖泊中螺類總β多樣性主要原因是物種在空間上的周轉(zhuǎn)。

其次, 采用配對(duì)方法計(jì)算βsim、βnes和βsor。βsor最小值為0.2414, 最大值為0.8824, 平均值為0.5810。βsim最小值為0, 最大值為0.8571, 平均值為0.3272, βsim在βsor中的比例平均為57.46%。βnes最小值為0, 最大值為0.8113, 平均值為0.2538, βnes在βsor中的比例平均為42.54%(圖5)。在βsim和βnes中均有0值出現(xiàn), 說明配對(duì)的湖泊之間螺類分布變化存在著完全的周轉(zhuǎn)(湖泊TH和NYH, LZH和NSHH, LZH和BAH, DATH和QDH, NSHH和BAH, DH和HMH, DHJH和HOUH之間, 嵌套部分βnes為0)以及完全的嵌套(湖泊PYH和CHH、LZH、CZH、DATH、QDH、DHJH、HOUH之間, DTH和DATH、HOUH之間, BAH和QLH之間, BDT和QLH、DHJH之間, 周轉(zhuǎn)部分βsim值為0)。

分別根據(jù)配對(duì)計(jì)算的βsim、βnes和βsor的值進(jìn)一步進(jìn)行非度量多維尺度分析(NMDS)和聚類分析(單連接聚合聚類法)。利用βsim聚類分析表明, 長(zhǎng)江中下游湖泊在0.10處分成3個(gè)類群, DH和HMH各為1群, 其余為一群。NMDS非約束排序分析的應(yīng)力函數(shù)值為0.2317, 表明二維點(diǎn)陣圖幾乎是隨機(jī)的, 不能很好解釋原始距離。利用βnes聚類分析表明, 長(zhǎng)江中下游湖泊在0.06處分成4個(gè)類群, DH和HMH為一群, PYH和DTH為一群, HH、BDT、TH、NYH為一群, 其余為第四群。其NMDS非約束排序分析的應(yīng)力函數(shù)值為0.0827, 表明二維點(diǎn)陣圖可以較好解釋長(zhǎng)江中下游湖泊螺類分布的關(guān)系。利用βsor聚類分析表明, 長(zhǎng)江中下游湖泊在0.45處分成4個(gè)類群, 第一群只有NSHH, 第二群包含DH和HMH, 第三群包含PYH和DTH, 其余14個(gè)湖泊為一群。NMDS非約束排序分析的應(yīng)力函數(shù)值為0.1594, 表明二維點(diǎn)陣圖對(duì)原始距離解釋能力一般(圖6)。

4 討論

長(zhǎng)江中下游湖泊是我國(guó)水生生物多樣性豐富地區(qū), 也是人口密度比較大, 農(nóng)業(yè)、工業(yè)和經(jīng)濟(jì)比較發(fā)達(dá)的地區(qū), 水生生物多樣性受到人類活動(dòng)很大的干擾[7]。螺類由于運(yùn)動(dòng)能力不強(qiáng), 更容易受到湖泊環(huán)境變化的影響, 圍湖造田、江河阻隔、水質(zhì)污染、水產(chǎn)養(yǎng)殖等都對(duì)螺類生存和分布造成很大影響, 加強(qiáng)螺類資源的保護(hù)是當(dāng)前很迫切的任務(wù)[9]。

圖4 LCBD值分布圖

圖5 βsim和βsne及其在βsor中所占比例分布

圖6 βsim, βsne and βsor NMDS分析和聚類分析

由于分布數(shù)量居中的物種帶來更多的變異(方差), 所以分布湖泊數(shù)量居中的螺類物種SCBD值應(yīng)最大, 而分布湖泊數(shù)很少和很多的物種SCBD值會(huì)降低, 即呈現(xiàn)單峰的關(guān)系[31]。本次研究中, 可能由于分布廣泛的螺類物種數(shù)很少而降低了其在相關(guān)分析中的作用, 因此檢驗(yàn)表明各螺類物種SCBD值和該物種分布湖泊數(shù)有顯著相關(guān)性(<0.001), 但從圖3中的散點(diǎn)分布和擬合線可以看出, 二者并不是呈簡(jiǎn)單的線性相關(guān), 如果分布廣泛的螺類物種數(shù)增多, 這種趨勢(shì)應(yīng)該更加明顯。螺類中的銅銹環(huán)棱螺在19個(gè)湖泊中都有分布, 直覺上該物種在19個(gè)湖泊中的β多樣性貢獻(xiàn)應(yīng)該為0, 但在Legendre方法中其SCBD值為0.0154, 原因是該方法中原始數(shù)據(jù)在計(jì)算前進(jìn)行了 “hellinger”轉(zhuǎn)化, 轉(zhuǎn)化后的數(shù)據(jù)和原始數(shù)據(jù)不是簡(jiǎn)單的線性關(guān)系, 因而轉(zhuǎn)化后的值在不同湖泊不一樣, 計(jì)算出的SCBD值不為0[32]。

在生態(tài)學(xué)上, LCBD的值反映了各湖泊螺類群落組成的獨(dú)特性程度。經(jīng)置換檢驗(yàn), 除牛山湖和花馬湖外, 其它湖泊LCBD值差異均不顯著, 說明螺類物種在這些湖泊中主要呈獨(dú)立、隨機(jī)的分布[5]?？赡茉蚴锹蓊愡\(yùn)動(dòng)能力弱, 活動(dòng)范圍小, 物種之間相互作用、相互競(jìng)爭(zhēng)影響不大, 因而成隨機(jī)分布。LCBD和湖泊中螺類的物種數(shù)沒有顯著的線性相關(guān)性, 在這些湖泊中, 鄱陽湖的螺類物種數(shù)最多, 但LCBD值最大的為花馬湖, 螺類種類數(shù)最少的是大黃家湖和后湖, 但LCBD值最小的是淀山湖。這可能和分布在湖泊中物種的稀有性有關(guān), 在淀山湖由于分布了較多的廣泛分布的物種反而降低了其LCBD值, 相反的, 分布較多的稀有物種則會(huì)增加該湖泊的LCBD值。

關(guān)于β多樣性起源的假說一般包括3個(gè)方面, 一是由于生物之間的相互作用, 在很大區(qū)域范圍內(nèi), 物種的組成是均勻分布的; 二是由于物種形成和擴(kuò)散歷史不同, 物種組成呈現(xiàn)為隨機(jī)的和自相關(guān)的形式; 三是由于環(huán)境的篩選、控制, 物種的分布和環(huán)境相關(guān)。如果物種組成變化是隨機(jī)的并且伴隨自相關(guān)的生物過程, 那么物種組成變化就主要體現(xiàn)在不同物種的替換, 即β多樣性可由周轉(zhuǎn), 即物種空間上的替換來解釋[5]。利用Baselga的加性分解方法, 對(duì)長(zhǎng)江中下游湖泊螺類β多樣性分解表明, 物種周轉(zhuǎn)成分在總β多樣性中占很大的比例(81.4%), 表明螺類β多樣性可能主要來源于物種隨機(jī)分布和自相關(guān)的過程?？赡茉蚴锹蓊惖倪\(yùn)動(dòng)緩慢, 擴(kuò)散能力現(xiàn)對(duì)較弱, 因而分布更多受物種形成歷史以及擴(kuò)散程度的影響。一般來說研究區(qū)域尺度不同, 影響β多樣性的主要因素不同。在宏觀尺度(不同進(jìn)化歷史的生物地理區(qū)域), 地球板塊運(yùn)動(dòng)造成的地理隔離使得物種形成和滅絕等是β多樣性主要推動(dòng)力; 在微觀尺度(環(huán)境異質(zhì)的棲息地), β多樣性主要受棲息地組成特性、物種生態(tài)位需求、干擾等因素影響; 本次研究是從中觀水平來分析長(zhǎng)江中下游湖泊螺類β多樣性, 因此影響因素主要是區(qū)域內(nèi)相對(duì)隔離的棲息地、氣候、物種擴(kuò)散能力等[33]。長(zhǎng)江中下游湖泊螺類在宏觀和微觀層次β多樣性變化規(guī)律有待進(jìn)一步研究。

圖5可以看出, Baselga配對(duì)的β多樣性分析方法中周轉(zhuǎn)和嵌套是兩個(gè)對(duì)立的過程, 結(jié)果表明不同湖泊對(duì)的β多樣性中周轉(zhuǎn)和嵌套組成比例不同。βnes在螺類物種數(shù)相等的湖泊間值為0, 例如湖泊LZH、NSHH、BAH分布螺類均為10種, DATH和QDH均為7種, DH和HMH均為29種, DHJH和HOUH均為5種。原因是由前文βsim和βsor計(jì)算公式可知b/(a+b)=2b/(2a+2b), 所以在配對(duì)的β多樣性分析中, βsim值和βsor值相等, 即這兩個(gè)湖泊beta多樣性完全是由于物種空間周轉(zhuǎn)[6]。在湖泊PYH和CHH、LZH、CZH、DATH、QDH、DHJH、HOUH之間, DTH和DATH、HOUH之間, BAH和QLH, BDT和QLH、DHJH之間則βsim值為0, 原因是在這些湖泊間, 螺類物種數(shù)少的湖泊的螺類也全部在螺類多的湖泊中, 即前者螺類是后者螺類的一個(gè)子集部分, 因而形成一個(gè)完全的嵌套格局。

根據(jù)不同湖泊對(duì)計(jì)算的βsim、βnes和βsor的值進(jìn)行聚類分析(單連接聚合聚類法), 結(jié)果顯示長(zhǎng)江中下游螺類可以分成3或4組, 主要反映了各湖泊中螺類分布種數(shù)的差異。利用βsim、βnes和βsor分別進(jìn)行非度量多維尺度分析(NMDS), 應(yīng)力函數(shù)結(jié)果顯示β多樣性不同組分二維點(diǎn)陣圖對(duì)原始距離的解釋能力不同, 利用βnes值NMDS分析的應(yīng)力函數(shù)值最小, 對(duì)原始數(shù)據(jù)的距離解釋能力最好。

理解一個(gè)區(qū)域內(nèi)生物多樣性如何分布才能更有效的保護(hù)該地區(qū)生物多樣性。我們利用多年長(zhǎng)江中下游湖泊螺類調(diào)查物種組成資料, 采用2種分析方法, 分析了長(zhǎng)江中下游湖泊中螺類分布的β多樣性, 表明不同湖泊、不同螺類物種對(duì)該區(qū)域螺類β多樣性貢獻(xiàn)度不同, 而且不同湖泊對(duì)之間螺類β多樣性形成主要原因不同, 因而整個(gè)區(qū)域螺類多樣性保護(hù)要采取不同的措施。例如, 在嵌套格局占優(yōu)勢(shì)的湖泊對(duì)中, 我們一般會(huì)保護(hù)物種數(shù)量多的湖泊, 但同時(shí)我們也要注意保護(hù)哪些面積不大但物種組成獨(dú)特即LCBD貢獻(xiàn)大的湖泊。因此, 將螺類分布的α多樣性和β多樣性結(jié)合起來, 可以為長(zhǎng)江中下游地區(qū)湖泊螺類資源利用和生物多樣性保護(hù)提供更合理的依據(jù)。

[1] WHITTAKER R H. Vegetation of the Siskiyou Mountains, Oregon and California[J]. Ecological Monographs, 1960(30): 279–338.

[2] LEGENDRE P, DE CACERES M. Beta diversity as the variance of community data: dissimilarity coefficients and partitioning [J]. Ecology Letters, 2013(16): 951–963.

[3] TUOMISTO H. A diversity of beta diversities: straightening up a concept gone awry, Part 1. Defining beta diversity as a function of alpha and gamma diversity [J]. Ecography, 2010(33): 2–22.

[4] ANDERSON M J, CRIST T O, CHASE J M, et al. Navigating the multiple meanings of beta diversity: a roadmap for the practicing ecologist[J] . Ecology Letters, 2011(14): 19–28.

[5] LEGENDRE P, BORCARD D, PERES-NETO P R. Analyzing beta diversity: partitioning the spatial variation of community composition data[J]. Ecological Monographs, 2005, 75(4): 435–450.

[6] BASELGA A. Partitioning the turnover and nestedness components of beta diversity[J]. Global Ecology and Biogeography, 2010(19): 134–143.

[7] 王洪鑄, 王海軍, 劉學(xué)勤, 等. 實(shí)施環(huán)境-水文-生態(tài)-經(jīng)濟(jì)協(xié)同管理戰(zhàn)略, 保護(hù)和修復(fù)長(zhǎng)江湖泊群生態(tài)環(huán)境[J]. 長(zhǎng)江流域資源與環(huán)境, 2015, 24(3): 353–357.

[8] 舒鳳月, 王海軍, 崔永德, 等. 長(zhǎng)江流域淡水軟體動(dòng)物物種多樣性及其分布格局[J]. 水生生物學(xué)報(bào), 2014, 38(1): 19–26.

[9] 舒鳳月, 王海軍, 潘保柱, 等. 長(zhǎng)江中下游湖泊貝類物種瀕危狀況評(píng)估[J]. 水生生物學(xué)報(bào), 2009, 33(6): 1051– 1058.

[10] JACOB B S, JAMES J G, WILLIAM E K, et al. How should beta-diversity inform biodiversity conservation? [J]. Trends in Ecology & Evolution, 2016, 31(1): 67–80.

[11] Dray S, Blanchet G, Borcard D, et al. adespatial: Multivariate Multiscale Spatial Analysis [EB/OL]., 2016.

[12] CRIBARI-NETO F, ZEILEIS A. Beta Regression in R[J]. Journal of Statistical Software, 2010(34): 1–24.

[13] ANDRéS B, DAVID L O C. betapart: an R package for the study of beta diversity[J]. Methods in Ecology and Evolution, 2012(3): 808–812.

[14] RDEVELOPMENT CORE TEAM. R: A Language and Environment for Statistical Computing [EB/OL]. R Foundation for Statistical Computing, Vienna, Austria. http://www. R-project. org/, 2016.

[15] KAMPSTRA P, BEANPLOT: a boxplot alternative for visual comparison of distributions[J]. Journal of Statistical Software, 2008(28): 1–9.

[16] OKSANEN J, F. BLANCHET G, KINDT R, et al. Vegan: community ecology package[EB/OL]. https://cran. r-project. org/web/packages/vegan/index. html, 2016.

[17] 吳小平, 梁彥齡, 王洪鑄, 等. 長(zhǎng)江中下游湖泊淡水貝類的分布及物種多樣性[J]. 湖泊科學(xué), 2000, 12(2): 111–118.

[18] 張銘華, 徐亮, 謝廣龍, 等. 鄱陽湖流域淡水貝類物種多樣性、分布與保護(hù)[J]. 海洋科學(xué), 2013, 37(8): 114–124.

[19] 陳曄光. 鄱陽湖及其周圍水域的淡水螺類[J]. 動(dòng)物學(xué)集刊, 1988(6): 69–75.

[20] 胡自強(qiáng). 洞庭湖及其周圍主要水域的螺類[J]. 湖南師范大學(xué)自然科學(xué)學(xué)報(bào), 1993, 16(1): 80–86.

[21] 蔡永久, 龔志軍, 秦伯強(qiáng). 太湖大型底棲動(dòng)物群落結(jié)構(gòu)及多樣性[J]. 生物多樣性, 2010, 18 (1): 50–59.

[22] 吳慶龍. 東太湖的貝類及其生物學(xué)[J]. 海洋湖沼通報(bào), 1993(4): 68–74.

[23] 胡菊英, 姚聞卿. 巢湖底棲動(dòng)物調(diào)查[J]. 安徽大學(xué)學(xué)報(bào)(自然科學(xué)版), 1981(2): 159–173.

[24] 李恒德, 呂代鈞. 洪湖底棲動(dòng)物調(diào)查報(bào)告[J]. 淡水漁業(yè), 1985(2): 25–29.

[25] 徐小雨, 周立志, 朱文中, 等. 安徽菜子湖大型底棲動(dòng)物的群落結(jié)構(gòu)特征[J]. 生態(tài)學(xué)報(bào), 2011, 31(4): 0943–0953.

[26] 陳立婧, 彭自然, 孫家平, 等. 安徽南漪湖大型底棲動(dòng)物群落結(jié)構(gòu)[J]. 動(dòng)物學(xué)雜志, 2008, 43(1): 63–68.

[27] 李德亮, 張婷, 余建波, 等. 典型養(yǎng)殖湖泊大通湖軟體動(dòng)物的時(shí)空分布格局[J]. 水生生物學(xué)報(bào), 2011, 35(1): 946– 954.

[28] 由文輝, 尤力群. 淀山湖軟體動(dòng)物群落的研究[J]. 華東師范大學(xué)學(xué)報(bào)(自然科學(xué)版), 1998(1): 103–109.

[29] 陳其羽, 梁彥齡, 宋貴保, 等. 武昌東湖軟體動(dòng)物的生態(tài)分布及種群密度[J]. 水生生物學(xué)集刊, 1975, 5(3): 371– 379.

[30] 陳其羽. 湖北省花馬湖軟體動(dòng)物的調(diào)查報(bào)告[J]. 海洋與湖沼, 1979, 10(1): 46–62.

[31] HEINO J, GR?NROOS M. Exploring species and site contributions to beta diversity in stream insect assemblages[J]. Oecologia, 2017(183): 151–160.

[32] RICOTTA C. Of beta diversity, variance, evenness, and dissimilarity [J]. Ecology and Evolution, 2017(7): 4835– 4843.

[33] PHILIP S B, SAUL A C, ADRIAN D M, et al. The spatial scaling of beta diversity[J]. Global Ecology and Biogeography. 2013(22): 639–647.

Analysis of the beta diversity of snails in the lakes of mid-lower reaches of the Yangtze River

WANG Weigen

School of Biology and Food Engineering of Fuyang Teachers College, Fuyang 236037, China

In this study, two kinds of partitioning measures were used to analyze the beta diversity of the snails in the lakes of mid-lower reaches of the Yangtze River. The first method proposed by Legendre was to explore the individual species contribution to total beta diversity (SCBD) and the single lake contribution to total beta diversity (LCBD). The results showed that the snail distributed in middle number of lakes, such as the, had the highest SCBD value, and the relationship between the number of lakes occupied by the snail and its SCBD value was hump-shaped. The Huama Lake had the highest LCBD value. There were no significant correlations between the LCBD values and the number of snail species in the same lake. The second approach proposed by Baselga was to partition turnover and nestedness components of beta diversity. Based on multiple-site analysis method, the results indicated that the turnover was main reason contribute to total beta diversity. In the case of pairwise patterns, the proportion of turnover and nestedness was different depended on the pair of lakes. There were some complete turnover and complete nestedness patterns between two lakes.

lake; snail; beta diversity; turnover; nestedness

10.14108/j.cnki.1008-8873.2018.06.016

Q958

A

1008-8873(2018)06-122-09

2017-11-24;

2018-01-18

安徽高校省級(jí)自然科學(xué)研究一般項(xiàng)目(2015KJ004); 阜陽師范學(xué)院自然科學(xué)研究項(xiàng)目(2015FSKJ01ZD)

王魏根(1975-), 男, 安徽宣城人, 碩士, 講師, 主要從事水生動(dòng)物學(xué)研究. E-mail: wweigen@163.com

王魏根. 長(zhǎng)江中下游湖泊螺類beta多樣性分析[J]. 生態(tài)科學(xué), 2018, 37(6): 122-130.

WANG Weigen. Analysis of the beta diversity of snails in the lakes of mid-lower reaches of the Yangtze River[J]. Ecological Science, 2018, 37(6): 122-130.