蔣萬祥, 陳　靜,*, 王紅妹, 賀詩水, 卓麗玲, 陳　青, 王洪凱, 蔡慶華

1 棗莊學院生命科學學院,棗莊　277160 2 中國科學院水生生物研究所淡水生態(tài)與生物技術(shù)國家重點實驗室,武漢　430072

生態(tài)系統(tǒng)功能是一個非常寬泛的概念,包括生態(tài)系統(tǒng)過程(如營養(yǎng)循環(huán)),生態(tài)系統(tǒng)服務(wù)(如提供漁獲物、棲息環(huán)境、過濾作用等),及對干擾的恢復力和抵抗力等[1- 2]。一直以來,河流生態(tài)系統(tǒng)功能研究絕大多數(shù)是基于傳統(tǒng)的物種分類開展的；然而,近年來的研究表明,生態(tài)系統(tǒng)功能主要受制于功能性狀的多樣性,即功能性狀的分布和豐度的時空格局[3]。功能性狀研究近20年的發(fā)展使得生態(tài)學家能夠從新的角度重新審視復雜的生態(tài)學過程。從個體、種群、群落到生態(tài)系統(tǒng),功能性狀已被證明是探索各類生態(tài)學前沿問題的有效手段[4]。

氣候變化引起整個流域內(nèi)環(huán)境特征包括水文特征、水體理化因子、河岸植被帶等隨季節(jié)發(fā)生相應的改變,并最終影響底棲動物的生長、繁殖和種群演替,使底棲動物群落結(jié)構(gòu)發(fā)生改變[5]；功能性狀研究對深入理解上述過程提供了一條新的路徑。相比大量關(guān)于功能性狀空間格局的研究,功能性狀時間動態(tài)方面研究遠遠不夠[4]。功能多樣性指數(shù)目前主要應用在陸生植被群落,水域生態(tài)系統(tǒng)遠比陸生植被系統(tǒng)復雜,且功能性狀的量化和獲取非常困難,因此,水生生物功能多樣性研究起步較晚[6]。溪流底棲動物功能性狀及多樣性相關(guān)研究在我國尚鮮見報道[7],而有關(guān)時間序列的研究更是未見報道。本文基于山東省新薛河底棲動物調(diào)查數(shù)據(jù),對功能性狀、功能多樣性的時空動態(tài)進行了統(tǒng)計分析,并就功能性狀和多樣性對環(huán)境壓力的響應進行了探討,以期為后續(xù)系統(tǒng)地開展相關(guān)研究(如土地利用、氣候變化等對生態(tài)系統(tǒng)功能的影響及生態(tài)系統(tǒng)服務(wù)功能等)提供背景數(shù)據(jù),為季節(jié)性河流的科學管理提供一定的理論支持。

1　材料與方法

1.1　區(qū)域概況與樣點設(shè)置

新薛河位于山東省西南部,屬于南四湖水系,全長89.6km,流域面積686km2。根據(jù)生境特征,沿水流方向共調(diào)查5種不同溪流生境的河段,每個河段沿水流方向設(shè)3個采樣點,共計15個。A河段,位于石嘴子水庫壩下,流態(tài)以緩流為主,沉積物較多。B河段位于棗莊市山亭區(qū)上游,河道頻繁出現(xiàn)干涸現(xiàn)象。C河段位于山亭區(qū)下游,城市污水在經(jīng)過沉積處理后直接排到該河段；D河段距離城區(qū)較遠,生境質(zhì)量相對較好。E河段間歇性地受家禽、衣物涮洗的影響。

1.2　底棲動物采集及鑒定

2012年4、10、12月,每月中旬用40目網(wǎng)徑、采樣面積0.09m2的索伯網(wǎng)(surber sampler)采集底棲動物,將網(wǎng)中的底棲動物揀出,裝標本瓶中,用10%的福爾馬林液固定,標本帶回實驗室后鑒定、計數(shù)[8- 10]。共采集標本42個(4月份B河段斷流)。

1.3　功能性狀和功能多樣性

功能性狀主要參考文獻[4,11-12],少數(shù)類群性狀根據(jù)標本和分類資料描述[8- 10]確定。選擇化性、漂移性、游泳能力、吸附能力、形狀、成熟個體大小、流態(tài)偏好、溫度偏好、生活型、營養(yǎng)習性等10個功能性狀納入分析。功能性狀選擇及其等級描述、代碼參考文獻[4,7],具體見表1。

表1　新薛河底棲動物功能性狀、性狀等級和代碼

功能多樣性采用Rao二次熵指數(shù)計算,該指數(shù)結(jié)合了物種相對豐度和生物性狀信息,將每個物種看作多維性狀空間中的點,主要計算物種距離的變異[13]。遵循Podani和Schemra的建議[14],計算功能多樣性時使用Gower距離來度量物種及生物性狀的距離。

式中,FDQ為Rao二次熵指數(shù),Pi和Pj分別代表物種i和物種j的相對豐度,S為物種數(shù),dij為物種i和物種j之間的距離,dij的取值范圍在0(物種i和物種j具有完全相同的生物性狀組成)和1(物種i和物種j未共有任何相似的生物性狀)之間。

1.4　數(shù)據(jù)處理

功能多樣性計算運用的軟件為R 3.1.3和FDiversity軟件包[15]。單因素方差分析(One-way ANOVA)、獨立樣本T檢驗(Independent samplesTtest)用于比較功能性狀、功能多樣性時空格局上的差異；使用的軟件為IBMSPSS 19.0。圖使用軟件Oringin8.0制作。

2　結(jié)果與分析

2.1　功能性狀組成

對功能性狀組成進行分析(表2),發(fā)現(xiàn)新薛河底棲動物10個功能性狀的組分(等級性狀)間均存在顯著差異(P<0.05)。具體為：化性表現(xiàn)為1世代/a占比最大,相對豐度為68.7%；漂移頻率較高的個體最多,相對豐度為60.8%；游泳能力表現(xiàn)為無游泳能力個體相對豐度高達97.1%；附著能力以無附著能力個體所占比重最大,為65.0%；成熟個體大小以小型個體為主,相對豐度為69.0%；流態(tài)偏好以沉積/侵蝕兼?zhèn)湫蛡€體所占比重最大,相對豐度為42.3%；溫度偏好以廣溫型個體所占比重最大,相對豐度為60.7%；形態(tài)方面以非流線型個體占絕大多數(shù),相對豐度為92.4%；生活型方面以穴居者占比最大,相對豐度為55.4%；營養(yǎng)習性方面以收集者占優(yōu),相對豐度為60.4%。

表2　新薛河底棲動物功能性狀組成和相對豐度/%

字母不同表示差異顯著(P<0.05)

2.2　等級性狀時空分布

2.2.1等級性狀空間格局

對等級性狀空間格局進行對比分析(圖1),結(jié)果表明33個等級性狀中有21個在空間序列上差異顯著(P<0.05)；涉及10個功能性狀中的9個,僅溫度偏好所屬等級性狀在各河段間差異不顯著(P>0.05)。具體為：身體形態(tài)表現(xiàn)為流線型在B、C河段相對豐度最??；非流線型個體在B、C河段豐度相對較高?；员憩F(xiàn)為1世代/a個體相對豐度在下游D、E河段高于上游A、B、C河段；2代或多代/a個體在D河段低于其他河段。漂移性表現(xiàn)為極少漂移個體相對豐度在A河段高于其他河段；一般漂移型個體相對豐度表現(xiàn)為A、B河段顯著低于C、D、E河段；高漂移性個體相對豐度在B、C河段高于其他河段。游泳能力表現(xiàn)為無游泳能力個體在C河段最高；強游泳能力個體在D、E河段高于其他河段。附著能力方面表現(xiàn)為無附著能力個體在B河段最高；有附著能力個體相對豐度在下游D、E河段高于其他河段。成熟個體大小方面,小型個體相對豐度B河段高于D河段；中等個體相對豐度D、E河段高于A、B、C河段。流態(tài)偏好方面,沉積型個體相對豐度在A、B、C河段高于D、E河段；侵蝕型個體在D、E河段相對豐度較高。生活型表現(xiàn)為穴居者相對豐度在A、B、C河段高于D、E河段；蔓生者在A河段高于其他河段；附著者在D、E河段相對豐度高于其他河段。營養(yǎng)習性方面,集食者在A、B、C河段豐度最大；濾食者在D、E河段豐度最大；捕食者在A河段高于其他河段。

圖1　各河段底棲動物等級性狀相對豐度Fig.1　Relative abundance of macroinvertebrates state traits in different river reaches圖中所示為差異顯著等級性狀,P <0.05；A代表A河段,B代表B河段,C代表C河段,D代表D河段,E代表E河段；各參數(shù)見表2

2.2.2等級性狀季節(jié)動態(tài)

對等級性狀時間格局進行對比分析(圖2),結(jié)果表明33個等級性狀中僅有12個在時間序列上差異顯著(P<0.05)；涉及10個功能性狀中的6個,化性、附著能力、形態(tài)、流態(tài)偏好4個功能性狀所屬等級性狀在各季節(jié)間差異不顯著(P>0.05)。具體為：強漂移性個體相對豐度10月≥4月≥12月。強游泳能力個體相對豐度在4月小于10和12月。個體大小方面,小型個體相對豐度在12月小于4月和10月；中等大小個體相對豐度12月大于10和4月；大型個體相對豐度10月高于其他月份。溫度偏好表現(xiàn)為喜熱型個體相對豐度在4月高于其他月份。生活型方面,穴居者相對豐度在4、10月高于12月；蔓生者相對豐度在4月份最大,10月份最小；附著者相對豐度12月份高于其他月份。營養(yǎng)習性方面,集食者相對豐度在4、10月高于12月；濾食者、捕食者相對豐度在12月份高于其他月份。

圖2　不同季節(jié)底棲動物等級性狀相對豐度Fig.2　Relative abundance of macroinvertebrates state traits among different seasons

2.3　功能多樣性

2.3.1功能多樣性空間格局

分析底棲動物功能多樣性空間格局(見圖3),發(fā)現(xiàn)D河段功能多樣性均值最高(0.119),E、A河段次之(分別為0.106、0.104),C、B河段最低(分別為0.094、0.067)；整體表現(xiàn)為D≥E、A≥C>B(P<0.05)。對各季節(jié)不同河段功能多樣性進行比較,結(jié)果表明,4月份各河段功能多樣性介于0.08—0.12之間,各河段功能多樣性指數(shù)差異不顯著(P>0.05)；10月份各河段功能多樣性介于0.05—0.15之間,表現(xiàn)為D≥C、E>A、B(P<0.05)；12月份各河段功能多樣性介于0.04—0.15之間,表現(xiàn)為A≥D、E≥C≥B(P<0.05)。

圖3　新薛河底棲動物功能多樣性空間格局Fig.3　The functional diversity among different river reaches in the New Xue River

2.3.2功能多樣性季節(jié)動態(tài)

在時間序列上(圖4),4月、10月、12月功能多樣性均值分別為0.107、0.095、0.100,差異不顯著(P>0.05)。具體到各河段： 4、10、12月功能多樣性均值在A河段分別為0.103、0.073、0.137；B河段分別為：0.073、0.060(4月份未采集)；C河段分別為0.107、0.087、0.090；D河段分別為：0.117、0.130、0.110；E河段分別為：0.103、0.113、0.103。方差分析和獨立樣本T檢驗表明,各季節(jié)功能多樣性僅在A河段差異顯著(P<0.05),表現(xiàn)為12月≥4月≥10月；其他河段功能多樣性于季節(jié)間差異均不顯著(P>0.05)。

圖4　新薛河底棲動物功能多樣性季節(jié)動態(tài)Fig.4　The functional diversity among different seasons in the New Xue River

3　討論

功能性狀能夠響應生存環(huán)境的變化,并對生態(tài)系統(tǒng)功能有一定影響[16]。物種在群落中的功能主要由其獲取和保存資源的能力、耐受競爭壓力和環(huán)境脅迫的性狀所決定[17]。性狀相對豐度的差異是生境過濾的結(jié)果,即相對豐度較高的等級性狀可認為是對區(qū)域環(huán)境適應較好的性狀[18]。通過方差分析確定了適宜新薛河生境的等級性狀,不僅增加了對底棲動物功能結(jié)構(gòu)的了解,同時,對研究區(qū)域河流狀況也有一定的反映。如漂移性能夠反映水生昆蟲規(guī)避不利環(huán)境的能力,高漂移性昆蟲豐度高,一定程度上反映了河流受干擾程度較高[19- 20]。附著能力、流態(tài)偏好、身體形態(tài)能夠反映河流水文條件,無吸附能力、侵蝕/沉積兼?zhèn)湫?、非流線型昆蟲豐度最高是對流域水文條件的綜合反映。一方面,就一個年度周期而言,流速在大多數(shù)時間應屬緩流,另一方面,季節(jié)性降水又使得水文條件出現(xiàn)波動。體型小的個體通常生活周期較短(r選擇),使群落在受到干擾后能后較快地恢復[21]；同時,個體小能夠增加底棲動物對水環(huán)境的利用能力(通常個體越小,生活需水量越少),對底棲動物規(guī)避旱季河道缺水意義重大[22]。營養(yǎng)習性是群落所在生態(tài)系統(tǒng)物質(zhì)循環(huán)和能量流動方式的反映,研究水域集食者所占比重最高,表明新薛河底棲動物在水域生態(tài)系統(tǒng)中主要扮演次級生產(chǎn)者的角色[23]。

生態(tài)系統(tǒng)功能本質(zhì)上依賴于物種所具有的功能性狀[24-25]。通過功能性狀來探討群落動態(tài)變化規(guī)律已經(jīng)成為有力可靠的方式,相比大量關(guān)于功能性狀格局的研究,功能性狀動態(tài)變化的研究遠遠不夠[26]。功能性狀在空間格局上存在較大差異是對環(huán)境變化的響應及在不同功能之間的權(quán)衡[27],分析結(jié)果同研究河段間生境高異質(zhì)性比較一致。河段A流速低、水深淺,由多個小型匯水區(qū)(pool)構(gòu)成,群落性狀組成主要表現(xiàn)為極少遷移型、蔓生型、捕食者個體所占比重最高,強游泳能力、有附著能力型、侵蝕型、濾食者個體相對較少。河段B最易出現(xiàn)斷流現(xiàn)象,水流以激流為主,其性狀特征主要表現(xiàn)為高遷移型、有附著能力、個體小型、侵蝕型、集食者個體少,流線型、濾食者個體占比較高。D、E河段流速相對較大,受人為干擾較小,群落性狀主要表現(xiàn)為強游泳型、吸附型、侵蝕型、附著型、濾食者個體相對豐度較高。

相對空間格局,等級性狀在時間序列上差異相對較小,化性、附著能力、形態(tài)、流態(tài)偏好等4個性狀所屬等級性狀在季節(jié)間差異均不顯著,表明底棲動物各功能性狀之間存在一定的權(quán)衡關(guān)系。權(quán)衡關(guān)系是功能性狀間各種關(guān)系中最普遍的一種[26],這種權(quán)衡關(guān)系是經(jīng)過自然篩選后形成的性狀組合,也稱“生態(tài)策略”,即物種沿一定的生態(tài)策略軸排列于最適應或最具競爭力的位置[28]。比如,Westoby分析了7個新熱帶森林中2134種木本植物的7種功能性狀間的關(guān)系,發(fā)現(xiàn)同一地點的物種共存其實是多種生活策略的穩(wěn)定組合[29]；Statzner等通過對歐洲地區(qū)底棲動物80年的監(jiān)測數(shù)據(jù)進行分析發(fā)現(xiàn),其生物形狀組成基本穩(wěn)定[30]。功能性狀能夠保持相對穩(wěn)定,主要受河流生境的唯一性決定[21]。在時間尺度上,本研究功能性狀組成雖然表現(xiàn)出一定的穩(wěn)定性,但相對上述研究也表現(xiàn)出一定的變異,主要是因為1)研究范圍較??；2)研究河流為季節(jié)性河流,B河段所在河段常有斷流、干涸現(xiàn)象；3)流域內(nèi)人口密集,干擾強度大。上述原因促使季節(jié)間生境質(zhì)量差異相對較大,導致底棲動物性狀組成發(fā)生變化；這和Snook等[31]研究結(jié)果較為一致。

功能多樣性可以用來分析生物對環(huán)境驅(qū)動因素的響應[27,32]。本研究底棲動物功能多樣性格局主要有2個特點：1)在空間格局上差異顯著；2)季節(jié)間差異不顯著(A河段除外)。由Rao二次熵指數(shù)[5,13,33]的計算過程可知,造成功能多樣性差異的主要原因是物種性狀的距離變異；一般而言,外部壓力越大,功能性狀趨同性越強,功能多樣性越小。功能多樣性在空間格局上差異顯著,是生境對適應該生境性狀的一類物種進行選擇的結(jié)果[4]。上述分析進一步表明,底棲動物功能多樣性對間歇性斷流(B河段)和污水的排放(C河段)等環(huán)境壓力具有良好的響應性。

功能多樣性季節(jié)動態(tài)比較符合生境模板理論,該理論認為,局地環(huán)境會對區(qū)域物種庫中的物種具有過濾篩選作用,使得群落生物性狀組成會逐漸同質(zhì)化,即在類似的生境中,即使物種組成不同,但其生物性狀組成是相似的[4]。此外,研究區(qū)域年降水量較少,降水頻率較低,河道底質(zhì)組成變化不大,加之性狀間權(quán)衡關(guān)系的影響,使得季節(jié)間功能性狀組成差異不大。A河段多樣性指數(shù)季節(jié)間變異顯著,主要是因為A河段靜水期、流水期交替出現(xiàn),水文節(jié)律變化相對較大所致。

致謝：馬相銘、李坤、李叢叢幫助野外調(diào)查，特此致謝。

參考文獻(References):

[1]Díaz S, Cabido M. Vive la différence: plant functional diversity matters to ecosystem processes. Trends in Ecology & Evolution, 2001, 16(11): 646- 655.

[2]Bremner J. Species′ traits and ecological functioning in marine conservation and management. Journal of Experimental Marine Biology and Ecology, 2008, 366(1): 37- 47.

[3]Elliott M, Quintino V. The estuarine quality paradox, environmental homeostasis and the difficulty of detecting anthropogenic stress in naturally stressed areas. Marine Pollution Bulletin, 2007, 54(6): 640- 645.

[4]Poff N L, Olden J D, Vieira N K M, Finn D S, Simmons M P, Kondratieff B C. Functional trait niches of North American lotic insect: traits-based ecological applications in light of phylogenetic relationships. Journal of the North American Benthological Society, 2006, 25(4): 730- 755.

[5]Rao C R. Diversity and dissimilarity coefficients: A unified approach. Theoretical Population Biology, 1982, 21(1):24- 43.

[6]張瀲波. 城鎮(zhèn)化對溪流底棲動物物種多樣性和功能多樣性影響的研究. 南京: 南京農(nóng)業(yè)大學, 2013:1- 18.

[7]蔣萬祥, 何逢志, 蔡慶華. 香溪河水生昆蟲功能性狀及功能多樣性空間格局. 生態(tài)學報, 2017, 37(6): 1861- 1870.

[8]Epler J H. Identification manual for the larval chironomidae (Diptera) of North and South Carolina. America: EPA Grant#X984170- 97. 2001: 1- 500.

[9]Morse J C, Yang L F, Tian L X. Aquatic Insects of China Useful for Monitoring Water Quality. Nanjing: Hohai University Press, 1994: 1- 568.

[10] Brinkhurst R O. Guide to the freshwater aquatic microdrile Oligochaetes of North America. Freshwater Science, 1987, 6(1):78- 79.

[11]Usseglio-Polatera P, Bournaud M, Richoux P, Tachet H. Biological and ecological traits of benthic freshwater macroinvertebrates: relationships and definition of groups with similar traits. Freshwater Biology, 2000, 43(2): 175- 205.

[12]Tachet H, Richoux P, Bournaud M, Usseglio-Polatera P. Invertebres d′Eau Douce: Systematique, Biologie, Ecologie. Paris: CNRS editions, 2002: 1- 552.

[13]Botta-Dukát Z. Rao′s quadratic entropy as a measure of functional diversity based on multiple traits. Journal of Vegetation Science, 2005, 16(5): 533- 540.

[14]Podani J, Schmera D. On Dendrogram-Based Measures of Functional Diversity. Oikos, 2006, 115(1): 179- 185.

[15]Casanoves F, Pla L, Di Rienzo J A, Díaz S. FDiversity: a software package for the integrated analysis of functional diversity. Methods in Ecology and Evolution, 2011, 2(3): 233- 237.

[16]賈鵬, 杜國禎. 生態(tài)學的多樣性指數(shù): 功能與系統(tǒng)發(fā)育. 生命科學, 2014, 26(2):153- 157.

[17]Wright I J, Reich P B, Westoby M, Ackerly D D, Baruch Z, Bongers F, Cavender-Bares J, Chapin T, Cornelissen J H C, Diemer M, Flexas J, Garnier E, Groom P K, Gulias J, Hikosaka K, Lamont B B, Lee T, Lee W, Lusk C, Midgley J J, Navas M L, Niinemets ü, Oleksyn J, Osada N, Poorter H, Poot P, Prior L, Pyankov V I, Roumet C, Thomas S C, Tjoelker M G, Veneklaas E J, Villar, R. The worldwide leaf economics spectrum. Nature, 2004, 428(6985): 821- 827.

[18]Menezes S, Baird D J, Amvm S. Beyond taxonomy: a review of macroinvertebrate trait-based community descriptors as tools for freshwater biomonitoring. Journal of Applied Ecology, 2010, 47(4): 711- 719.

[19]Wise E J. Seasonal distribution and life histories of Ephemeroptera in a Northumbrian river. Freshwater Biology, 1980, 10(2): 101- 111.

[20]Vannote R L, Sweeney B W. Geographic analysis of thermal equilibria: a conceptual model for evaluating the effect of natural and modified thermal regimes on aquatic insect communities. The American Naturalist, 1980, 115(5): 667- 695.

[21]Bêche L A, Mcelravy E P, Resh V H. Long-term seasonal variation in the biological traits of benthic-macroinvertebrates in two Mediterranean-climate streams in California, U.S.A. Freshwater Biology, 2006, 51(1): 56- 75.

[22]Townsend C R. The Patch Dynamics Concept of Stream Community Ecology. Freshwater Science, 1989, 8(1): 36- 50.

[23]蔣萬祥, 蔡慶華, 唐濤, 渠曉東. 香溪河水系大型底棲動物功能攝食類群生態(tài)學. 生態(tài)學報, 2009, 29(10): 5207- 5218.

[24]Diaz S, Hodgson J G, Thompson K, Cabido M, Cornelissen J H C, Jalili A, Montserrat-Martí G, Grime J P, Zarrinkamar F, Asri Y, Band S R, Basconcelo S, Castro-Díez P, Funes G, Hamzehee B, Khoshnevi M, Pérez-Harguindeguy N, Pérez-Rontomé M C, Shirvany F A, Vendramini F, Yazdani S, Abbas-Azimi R, Bogaard A, Boustani S, Charles M, Dehghan M, de Torres-Espuny L, Falczuk V, Guerrero-Campo J, Hynd A, Jones G, Kowsary E, Kazemi-Saeed F, Maestro-Martínez M, Romo-Díez A, Shaw S, Siavash B, Villar-Salvador P, Zak M R. The plant traits that drive ecosystems: evidence from three continents. Journal of vegetation science, 2004, 15(3): 295- 304.

[25]Hooper D U, Vitousek P M. The effects of plant composition and diversity on ecosystem processes. Science, 1997, 277(5330): 1302- 1305.

[26]劉曉娟, 馬克平. 植物功能性狀研究進展. 中國科學: 生命科學, 2015, 45(4): 325- 339.

[27]肖玉, 謝高地, 安凱, 魯春霞. 基于功能性狀的生態(tài)系統(tǒng)服務(wù)研究框架.植物生態(tài)學報, 2012, 36(4): 353- 362.

[28]Wright I J, Ackerly D D, Bongers F, Harms K E, Ibarra-Manriquez G, Martinez-Ramos M, Mazer S J, Muller-Landauhc H C, Paz H, Pitman N C A, Poorter L, Silman M R, Vriesendorp C F, Webb C O, Westobym, Wright S J. Relationships among ecologically important dimensions of plant trait variation in seven Neotropical forests. Annals of Botany, 2007, 99(5): 1003- 1015.

[29]Westoby M. A leaf-height-seed (LHS) plant ecology strategy scheme. Plant and soil, 1998, 199(2): 213- 227.

[30]Statzner B, Dolédec S, Hugueny B. Biological Trait Composition of European Stream Invertebrate Communities: Assessing the Effects of Various Trait Filter Types. Ecography, 2004, 27(4): 470- 488.

[31]Snook D L, Milner A M. Biological traits of macroinvertebrates and hydraulic conditions in a glacier-fed catchment (French Pyrenées). Archiv Für Hydrobiologie, 2002, 153(2): 245- 271.

[32]Mouchet M A, Villéger S, Mason N W H, Mouillot D. Functional diversity measures: an overview of their redundancy and their ability to discriminate community assembly rules. Functional Ecology, 2010, 24(4): 867- 876.

[33]Pavoine S, Ollier S, Pontier D. Measuring diversity from dissimilarities with Rao′s quadratic entropy: are any dissimilarities suitable? Theoretical Population Biology, 2005, 67(4): 231- 239.