應(yīng)用浮游植物生物完整性指數(shù)評價長江上游河流健康

譚 巧，馬芊芊，李斌斌，呂紅健，付 梅，姚維志

(西南大學(xué)動物科技學(xué)院，西南大學(xué)漁業(yè)資源環(huán)境研究中心，水產(chǎn)科學(xué)重慶市重點實驗室，重慶 400716)

應(yīng)用浮游植物生物完整性指數(shù)評價長江上游河流健康

譚 巧，馬芊芊，李斌斌，呂紅健，付 梅，姚維志

(西南大學(xué)動物科技學(xué)院，西南大學(xué)漁業(yè)資源環(huán)境研究中心，水產(chǎn)科學(xué)重慶市重點實驗室，重慶 400716)

2013年11月-2014年6月對長江上游宜賓至江津段五個斷面的浮游植物進行了四次調(diào)查，共鑒定出浮游植物6門38屬95種。對13個備選指標進行分布范圍、判別能力及Pearson相關(guān)性分析，構(gòu)建了適合長江上游的浮游植物生物完整性指數(shù)(P-IBI)指標體系，即香農(nóng)多樣性指數(shù)、Margalef豐富度指數(shù)、Pielou均勻度指數(shù)、浮游植物密度、硅藻密度百分比等參數(shù)指標。采用3分制法、4分制法和比值法分別對生物指標計分，評價結(jié)果顯示：高莊橋、羊石處于“健康”狀態(tài)；白沙處于“亞健康”狀態(tài)；而江安、德感處于“亞健康”向“一般”過渡狀態(tài)。綜合來看，三種評價方法所反映各樣點的健康狀況基本一致，只是4分制法和比值法在劃分評價等級上更細致，評價結(jié)果更精確。Pearson相關(guān)性分析顯示：TN、pH與P-IBI值呈顯著負相關(guān)(P<0.05)。

長江上游；浮游植物；生物完整性指數(shù)；生態(tài)系統(tǒng)健康評價

河流健康的概念由美國的《清潔水法》于1972年首次提出，該法案認為物理、化學(xué)和生物的完整性是河流健康的標準，其中完整性指維持生態(tài)系統(tǒng)的自然結(jié)構(gòu)和功能的狀態(tài)[1]。近年來，隨著人們對河流健康重要性認識的加深，國內(nèi)外學(xué)者在其定義、內(nèi)涵及研究方法等方面進行了大量的探索[2-9]。Karr[10]認為一個良好的生態(tài)系統(tǒng)必定存在一個完整的生物群落，提出用生物完整性指數(shù)(Index of Biotic Integrity，IBI)來評價河流健康。IBI最初以魚類為研究對象，后又應(yīng)用到其它生物[1-2,11]。目前，在河流健康評價中主要應(yīng)用著生藻類、底棲無脊椎動物以及魚類生物完整性指數(shù)，浮游生物完整性指數(shù)的應(yīng)用相對較少。蔡琨等[12]應(yīng)用浮游植物生物完整性指數(shù)(P-IBI)評價了冬季太湖水生態(tài)健康質(zhì)量，認為浮游植物生物完整性指數(shù)盡管存在局限性，但用其來評價水生態(tài)健康是可行的，其研究結(jié)果也表明P-IBI是相對科學(xué)準確的。沈強等[13]應(yīng)用浮游生物完整性指數(shù)對浙江水源地水質(zhì)進行了評價。徐玉萍[14]、馬亮[15]分別對上海市河道、潭江的浮游植物生物完整性做過研究。雖然針對長江上游浮游生物方面的研究較多[16-19]，但局限于群落結(jié)構(gòu)及多樣性評價，未見應(yīng)用浮游生物完整性指數(shù)進行健康評價的報道。本研究將在此方面做出探討，希望能為今后的研究提供幫助。

1 材料與方法

1.1 采樣點設(shè)置與采樣時間

以長江上游為研究對象，在宜賓至江津段共設(shè)置了5個采樣斷面，采樣點從上游到下游依次為高莊橋左岸(金沙江段，距宜賓航道里程5.1 km)、江安右岸(距宜賓航道里程64.4 km)、羊石右岸(距宜賓航道里程213.8 km)、白沙右岸(距宜賓航道里程261.4 km)、德感左岸(距宜賓航道里程299.7 km)。其中高莊橋、羊石位于長江上游珍稀特有魚

類國家級自然保護區(qū)的核心區(qū)，江安位于緩沖區(qū)，白沙、德感位于實驗區(qū)。采樣時間為2013年11月、2014年1、4、6月，每月采樣一次。

1.2 樣品采集與處理

采集定量標本，使用有機玻璃采水器在水面下0.5 m處取水樣1 L，現(xiàn)場加入15 mL魯哥氏液固定。另外，用25#浮游生物網(wǎng)在表層至0.5 m深處作“∞”形緩慢拖動2 min左右以撈取定性標本。先對除硅藻門外其他門類的種類進行鑒定，再用酸液消解除去硅藻細胞內(nèi)含物，以準確鑒定硅藻種類。種類鑒定參照《中國淡水藻類》[21]和《淡水微型生物圖譜》[22]。

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

1.3.1 參考點與污染點的區(qū)分

根據(jù)香農(nóng)多樣性指數(shù)(H)水質(zhì)劃分標準[23]，將 H≥3的采樣點定為參考點，H<3的采樣點定為污染點。

1.3.2 備選指標的確定與篩選

選取了13個應(yīng)用廣泛且入選率較高的指標作為備選指標(表1)。對這13個備選指標進行分布范圍、判別能力以及相關(guān)性分析，以篩選出合適的參數(shù)指標。

表1 P-IBI備選指標描述及對干擾的反映Tab.1 Description of candidate metrics for P-IBI and their expected direction of response to disturbance

續(xù)表1

1.3.3 P-IBI指標的計分及評價標準

(1) 3分制法

對隨干擾減小的指標，高于參考點的25%分位數(shù)值(對隨干擾增大的指標，低于參考點的75%分位數(shù)值)，計6分；低于25%分位數(shù)值(高于75%分位數(shù)值)的分布范圍進行二等分，分別計3分和0分。

(2) 4分制法

對隨干擾減小的指標，以所有點的95%分位數(shù)值為最佳期望值(對隨干擾增大的指標，以所有點的5%分位數(shù)值為最佳期望值)，低于(高于)此值的分布范圍進行三等分，依次計6、4、2、0分。

(3) 比值法

對隨干擾減小的指標，各指標的分值等于指標值除以所有點的95%分位數(shù)值；對隨干擾增大的指標，以所有點的5%分位數(shù)值為最佳期望值進行計算，方法為：(最大值-指數(shù)值)/(最大值-5%分位數(shù)值)。該法規(guī)定，經(jīng)計算后的分值的分布范圍為0～1，若大于1，則記為1。

(4) 評價標準

以參考點IBI值的25%分位數(shù)作為健康評價的標準，若采樣點的IBI值大于25%分位數(shù)值，則表示該樣點受到的干擾很小，是健康的；對小于25%分位數(shù)值的分布范圍進行四等分，就可以確定出健康、亞健康、一般、較差、極差五個健康等級。

2 結(jié)果

2.1 浮游植物的群落結(jié)構(gòu)特征

本研究共鑒定出浮游植物6門38屬95種(變種)，其中硅藻門22屬67種，占藻類總種數(shù)的71%；綠藻門10屬18種，占19%；隱藻門1屬4種，占4%；藍藻門2屬2種，占2%；裸藻門1屬2種，占2%；甲藻門2屬2種，占2%。

2.2 備選指標的篩選

分析各指標值的分布范圍，刪除分布范圍太窄或可預(yù)測環(huán)境變化值的范圍太小的指標，結(jié)果備選指標全部適合。用箱線圖法進行判別能力分析(圖1)，香農(nóng)多樣性指數(shù)(M3)、Margalef豐富度指數(shù)(M4)、Pielou均勻度指數(shù)(M5)、浮游植物密度(M9)、硅藻密度百分比(M11)、前三優(yōu)勢種百分比(M13)等6個參數(shù)指標在參考點與污染點間有顯著差異，可以進入下一步分析。Pearson相關(guān)性檢驗(SPSS 19.0)(表2)顯示，M4、M9、M11與其它指標間的相關(guān)性不顯著，直接保留；M3、M5分別與M13顯著相關(guān)，舍棄M13，保留M3、M5。

對所有備選指標進行篩選后，確定長江上游P-IBI指標體系由香農(nóng)多樣性指數(shù)(M3)、Margalef豐富度指數(shù)(M4)、Pielou均勻度指數(shù)(M5)、浮游植物密度(M9)、硅藻密度百分比(M11)構(gòu)成。

圖1 各入選P-IBI指標在參考點和污染點的箱線圖Fig.1 Discriminatory power of selected metrics for reference and impaired sites using box & whisker plots表2 備選P-IBI指標的Pearson相關(guān)性分析Tab.2 Pearson’s correlating coefficient between six candidate biological metrics

備選指標M3M4M5M9M11M13M31M40.65531M5-0.36680.19901M90.6799-0.0490-0.67991M11-0.3639-0.3199-0.21310.24971M13-0.8298*-0.4124-0.9266*0.58010.30201

注：*表示P<0.05。

2.3 P-IBI指標分值計算及計分標準

將篩選出的指標分別應(yīng)用3分制法、4分制法、比值法計算了各指標分值(表3、4、5)。

表3 5個P-IBI指標在參照點中的分布及3分制法計分標準Tab.3 Distribution of the five metrics in reference sites and score standards according to three score system

表4 5個P-IBI指標在所有樣點中的分布及4分制法計分標準Tab.4 Distribution of the five metrics in all the sampling sites and score standards according to four score system

表5 比值法計算5個P-IBI指標分值的公式Tab.5 Formulas for calculating five metrics scores according to ratio score system

2.4 P-IBI指標體系評價標準及評價結(jié)果

2.4.1 3分制法

根據(jù)3分制法的計分標準(表3)，將篩選所得指標的分值進行加和，即獲得P-IBI值。以參考點的25%分位數(shù)為最佳期望值，確定P-IBI指標體系評價標準(表6)，進而得出各采樣點浮游植物生物完整性評分及健康評價(表7)。四次采樣中高莊橋評價為三個“健康”、一個“亞健康”；江安評價為兩個“亞健康”、兩個“一般”；羊石評價為三個“健康”、一個“亞健康”；白沙評價為四個“健康”；德感評價為兩個“健康”，兩個“一般”。

表6 長江上游P-IBI 3分制法評價標準Tab.6 Assessment criteria for biological integrity of phytoplankton in the upper reaches of the Yangtze River according to three score system

表7 上江上游P-IBI 3分制法評價結(jié)果Tab.7 P-IBI results for each sampling sites in the upper reaches of the Yangtze River according to three score system

續(xù)表7

2.4.2 4分制法

4分制法的評價標準與3分制法相同，據(jù)此標準得出各采樣點浮游植物生物完整性評分及健康評價(表8)。四次采樣中高莊橋評價為四個“健康”；江安評價為兩個“健康”、一個“亞健康”、一個“一般”；羊石評價為四個“健康”；白沙評價為三個“健康”、一個“亞健康”；德感評價為兩個“健康”，一個“亞健康”、一個“一般”。

表8 長江上游P-IBI 4分制法評價結(jié)果Tab.8 P-IBI results for each sampling sites in the upper reaches of the Yangtze River according to four score system

續(xù)表8

2.4.3 比值法

確定比值法的P-IBI指標體系評價標準(表9)，進而得出本次試驗所有采樣點浮游植物生物完整性評分及健康評價(表10)。四次采樣中高莊橋評價為三個“健康”、一個“亞健康”；江安評價為一個“健康”、兩個“亞健康”、一個“一般”；羊石評價為四個“健康”；白沙評價為三個“健康”、一個“亞健康”；德感評價為兩個“健康”，一個“亞健康”、一個“一般”。

表9 長江上游P-IBI比值法評價標準Tab.9 Assessment criteria for biological integrity of phytoplankton in the upper reaches of the Yangtze River according to ratio score system

表10 長江上游P-IBI比值法評價結(jié)果Tab.10 P-IBI results for each sampling sites in the upper reaches of the Yangtze River according to ratio score system

續(xù)表10

2.5 P-IBI與水質(zhì)參數(shù)的一致性

2.5.1 主要水質(zhì)參數(shù)分析

表11 各采樣點主要水質(zhì)參數(shù)Tab.11 Main chemical and physical information of sampling sites

2.5.2 P-IBI與環(huán)境因子之間的相關(guān)性

表12 P-IBI與環(huán)境因子之間的相關(guān)性Tab.12 Pearson’s correlating coefficient between P-IBI and environmental factors

注：*表示P<0.05。

3 討論

3.1 不同評分方法的比較

本研究采用3分制法、4分制法和比值法三種評分方法來計算P-IBI值。比較三種方法評價結(jié)果，4分制法與比值法高度相似，僅高莊橋、江安采樣點存在一個“健康”與“亞健康”的區(qū)別；而3分制法與4分制法、比值法相比，其評價結(jié)果主要為“健康”和“一般”，較少出現(xiàn)“亞健康”狀態(tài)，顯得跨度較大。3分制法是根據(jù)參考點的指標值分布范圍，將其劃分為三部分，由高到低分別計6、3、0分；而4分制法是根據(jù)所有點的指標值分布范圍，將其劃分為四部分，分別計6、4、2、0分，因此，評價結(jié)果的差別可能體現(xiàn)在對各指標計分時劃分的細致程度上。裴雪姣等[24]在應(yīng)用魚類完整性評價遼河流域健康中也認為，1、3、5 賦值法 IBI 分值為各參數(shù) 1、3、5 記分的累加，分值連續(xù)性較小，間斷性較大，對最高和最低等級間的區(qū)分性比較明顯，而比值法的分值則連續(xù)性較強，間斷性較小。王備新[23]的研究表明：用比值法統(tǒng)一各生物參數(shù)量綱后，依據(jù)所有樣點數(shù)據(jù)建立的評價標準與依據(jù)參照點數(shù)據(jù)所得標準相比，提高了對參照點的準確率，結(jié)果優(yōu)于3分制法和4分制法。從本研究的結(jié)果來看，3分制法所得P-IBI值分別與4分制法、比值法所得P-IBI值呈顯著正相關(guān)(P<0.05)，而4分制法所得P-IBI值與比值法所得P-IBI值呈極顯著正相關(guān)(P<0.01)，因此，它們所反映各樣點的健康狀況基本一致，只是4分制法和比值法在劃分評價等級上更細致，評價結(jié)果更精確。

3.2 應(yīng)用P-IBI評價長江上游河流健康狀況

長江上游珍稀特有魚類自然保護區(qū)是我國最長的河流型自然保護區(qū)，該江段自然生境相對完整，但近年來，受采砂、航運、涉水工程、工業(yè)廢水、生活污水的影響，長江上游河流健康度逐年下降[25]。高莊橋在宜賓主城上游、羊石鎮(zhèn)遠離大城市，且周圍人口密度不大，工農(nóng)業(yè)活動相對較少，在本研究中評價為“健康”；江安江段位于宜賓市下游，前方有金沙江、岷江匯入，受到城鎮(zhèn)生活污水、工業(yè)廢水、航運等影響，評價結(jié)果為“亞健康”向“一般”過渡階段；白沙江段評價為“亞健康”；德感江段處于城鎮(zhèn)，受到人類干擾相對較大，評價為“亞健康”向“一般”過渡狀態(tài)。表明保護區(qū)內(nèi)高莊橋、羊石江段水質(zhì)比較清潔，而江安、白沙、德感江段受到一定干擾，這與各采樣點主要水質(zhì)參數(shù)分析的結(jié)果吻合。因此，應(yīng)用浮游植物生物完整性指數(shù)評價長江上游河流健康狀況是可行的。

目前，浮游植物生物完整性指數(shù)主要被應(yīng)用于評價一個相對靜止的水體如水庫[13, 26]、湖泊[12, 27]，評價江河健康狀況的研究比較少。分析原因認為：浮游植物漂流在江河中，樣本對采樣點的代表性和采樣的均一性受到質(zhì)疑；生活周期短，不能反映采樣點長時間的生態(tài)健康狀況[12]。但是，浮游植物對環(huán)境變化反應(yīng)靈敏，能在短時間內(nèi)做出響應(yīng)，可以通過設(shè)置合理的采樣時間、采樣頻次來構(gòu)建浮游植物生物完整性。Sommer等[28]提出的PEG(plankton ecology group)模型認為，浮游植物群落存在冬春季時隱藻和硅藻占優(yōu)勢，夏季綠藻占優(yōu)勢，夏末秋初是藍藻占優(yōu)勢，而到秋季時硅藻數(shù)量再次上升的自然演替規(guī)律。針對浮游植物群落的季節(jié)性演替規(guī)律，本研究設(shè)置采樣時間為11月、次年1、4、6月，研究結(jié)果較好地反映了長江上游宜賓至江津段健康狀況。建議以后的研究根據(jù)浮游植物群落的季節(jié)性演替規(guī)律增加采樣頻次，保證樣品的代表性和采樣的均一性，構(gòu)建更加科學(xué)合理的浮游植物生物完整性指數(shù)。

[1] Karr J R.Defining and measuring river health [J].Freshwater Biol, 1999, 41(2):221-234.

[2] An K, Park S S, Shin J Y.An evaluation of a river health using the index of biological integrity along with relations to chemical and habitat conditions [J].Environ Int, 2002, 28 (5):411-420.

[3]Scrimgeour G J, Wicklum D.Aquatic ecosystem health and integrity: Problems and potential solutions [J].J North Am Benthol Soc, 1996, 15(2):254-261.

[4]Meyer J L.Stream health: Incorporating the human dimension to advance stream ecology [J].J North Am Benthol Soc, 1997, 16(2):439-447.

[5]Costanza R, Mageau M.What is a healthy ecosystem? [J].Aquat Ecol, 1999, 33(1):105-115.

[6]Norris R H, Thoms M C.What is river health? [J].Freshwater Biol, 1999, 41(2):197-209.

[7]文伏波, 韓其為, 許炯心, 等.河流健康的定義與內(nèi)涵[J].水科學(xué)進展, 2007,18(1):140.

[8]董哲仁.河流健康的內(nèi)涵[J].中國水利, 2005,(4):15-18.

[9]劉昌明, 劉曉燕.河流健康理論初探[J].地理學(xué)報, 2008,63(7):683-692.

[10]Karr J R.Assessment of biotic integrity using fish communities [J].Fisheries, 1981, 6(6):21-27.

[11]Kerans B L, Karr J R.A benthic index of biotic integrity (B-IBI) for rivers of the Tennessee Valley[J].Ecol Appl, 1994, 4(4):768-785.

[12]蔡 琨, 秦春燕, 李繼影, 等.基于浮游植物生物完整性指數(shù)的湖泊生態(tài)系統(tǒng)評價——以2012年冬季太湖為例[J].生態(tài)學(xué)報, 2016, 36(5):1431-1441.

[13]沈 強, 俞建軍, 陳 暉, 等.浮游生物完整性指數(shù)在浙江水源地水質(zhì)評價中的應(yīng)用[J].水生態(tài)學(xué)雜志, 2012, 33(2):26-31.

[14] 徐玉萍.上海市河道生物完整性指數(shù)評價研究及生態(tài)治理效果評價[D].上海：上海海洋大學(xué), 2014.

[15] 馬 亮.潭江浮游植物完整性指數(shù)及水質(zhì)評價[D].河北保定：河北大學(xué), 2012.

[16]陳燕琴, 申志新, 劉玉婷, 等.長江上游曲麻萊至玉樹段春秋季浮游植物群落結(jié)構(gòu)及多樣性評價[J].長江流域資源與環(huán)境, 2013, 22(10):1325-1332.

[17]余海英.長江上游珍稀、特有魚類國家級自然保護區(qū)浮游植物和浮游動物種類分布和數(shù)量研究[D].重慶：西南大學(xué), 2008.

[18]劉曉娟, 李云峰, 茹輝軍, 等.長江上游珍稀、特有魚類國家級自然保護區(qū)干流段浮游植物群落結(jié)構(gòu)的年度變化[J].淡水漁業(yè), 2015, 45(1):52-57.

[19]何 滔, 劉建虎, 張春霖, 等.長江上游支流抱龍河浮游生物現(xiàn)狀及多樣性評價[J].淡水漁業(yè), 2014, 44(3):51-55.

[20]胡鴻鈞，李堯英，魏印心，等.中國淡水藻類[M].上海：上?？茖W(xué)技術(shù)出版社, 1980.

[21]周鳳霞, 陳劍虹.淡水微型生物圖譜[M].北京：化學(xué)工業(yè)出版社, 2005.

[22]魏復(fù)盛.水和廢水監(jiān)測分析方法[M].第4版.北京：中國環(huán)境科學(xué)出版社, 2002.

[23]王備新.大型底棲無脊椎動物水質(zhì)生物評價研究[D].南京：南京農(nóng)業(yè)大學(xué), 2003.

[24]裴雪姣, 牛翠娟, 高 欣, 等.應(yīng)用魚類完整性評價體系評價遼河流域健康[J].生態(tài)學(xué)報, 2010, 30(21):5736-5746.

[25]危起偉，李堯英，魏印心，等.長江上游珍稀特有魚類國家級自然保護區(qū)科學(xué)考察報告[M].北京：科學(xué)出版社, 2012.

[26]婁方瑞, 劉明濤, 程光平, 等.巖灘水庫水質(zhì)與浮游生物完整性灰色關(guān)聯(lián)分析[J].廣東農(nóng)業(yè)科學(xué), 2015, (17):147-152.

[27]賀方兵.東部淺水湖泊水生態(tài)系統(tǒng)健康狀態(tài)評估研究——以湖北省為例[D].重慶：重慶交通大學(xué), 2015.

[28]Sommer U, Gliwicz Z M, Lampert W, et al.The PEG-model of seasonal succession of planktonic events in fresh waters [J].Arch Hydrobiol, 1986, 106(4):433-471.

(責(zé)任編輯：鄧 薇)

Ecological health assessment of the upper reaches of the Yangtze River, based on biotic integrity index of phytoplankton

TAN Qiao, MA Qian-qian, LI Bin-bin, Lü Hong-jian,FU Mei, YAO Wei-zhi

(SchoolofAnimalScienceandTechnology；ResearchCentralofFisheryResourcesandEnvironment；KeyLaboratoryofAquaticScienceofChongqing；SouthwestUniversity,Chongqing, 400716,China)

In the present study, phytoplankton samples were collected from the upper reaches of the Yangtze River during four consecutive seasons, from November 2013 to November 2014.The results showed that there were 95 phytoplankton species belonging to 6 phyla and 38 genera found in the investigating region.And based on the distributing range，discriminatory power and Pearson's correlation analysis of the 13 candidate indexes, five metrics, including Shannon-Wiener index, Margalef index, Pielou index, density of phytoplankton, percentage of bacillariophyta were selected to establish biotic integrity index of phytoplankton (P-IBI) evaluation system.Three score system, four score system and ratio score system were used to get metrics into a uniform score for all the sampling points, and the results showed that Gao Zhuangqiao and Yang Shi were in the state of health, Bai Sha in sub-health, and Jiang An and De Gan in the sub-health to good-fair transitional phase.On the whole, the health state of each sampling sites reflected by the three scoring methods almost the same.However, four score system and ratio score system were more meticulous on the grade classification and evaluation, thus the evaluation results of four score system and ratio score system were more accurate.The Pearson′s correlation analysis showed that the P-IBI was significantly negative correlated with TN and pH (P<0.05).

the upper reaches of the Yangtze River; phytoplankton ; index of biotic integrity；ecosystem health assessment

2016-09-10；

2017-01-05

美國大自然保護協(xié)會(TNC)“長江上游珍稀特有魚類國家級自然保護區(qū)保護行動實施情況評估研究”；農(nóng)業(yè)部長江流域漁政監(jiān)督管理辦公室“長江上游保護區(qū)干流江段棲息地評估與修復(fù)研究”；中央高?；究蒲袠I(yè)務(wù)費專項資金資助(XDJK2015A001)

譚 巧(1991-)，男，碩士研究生，專業(yè)方向為漁業(yè)資源與環(huán)境。E-mail:752043576@qq.com

姚維志。E-mail: yaowz@swu.edu.cn

X832

A

1000-6907-(2017)03-0097-08