陶 敏 謝碧文 齊澤民 楊品紅 李 斌①

(1. 內(nèi)江師范學(xué)院生命科學(xué)學(xué)院 長(zhǎng)江上游魚(yú)類資源保護(hù)與利用四川省重點(diǎn)實(shí)驗(yàn)室 內(nèi)江 641112; 2. 湖南文理學(xué)院 水產(chǎn)高效健康生產(chǎn)湖南省協(xié)同創(chuàng)新中心 環(huán)洞庭湖水產(chǎn)健康養(yǎng)殖及加工湖南省重點(diǎn)實(shí)驗(yàn)室 動(dòng)物學(xué)湖南省高校重點(diǎn)實(shí)驗(yàn)室 常德415000)

水質(zhì)是水生態(tài)系統(tǒng)健康的重要指標(biāo), 進(jìn)行水質(zhì)評(píng)價(jià)是水體開(kāi)發(fā)利用、保護(hù)與生態(tài)恢復(fù)的前提與基礎(chǔ)(焦珂?zhèn)サ? 2015)。早期的水質(zhì)評(píng)價(jià)多針對(duì)理化指標(biāo)(Bozzetti et al, 2004)。然而由于環(huán)境與生物之間的復(fù)雜相互作用, 僅僅用理化指標(biāo)無(wú)法完整體現(xiàn)水體的健康狀況。作為水域的初級(jí)生產(chǎn)者, 浮游植物是整個(gè)水生態(tài)系統(tǒng)物質(zhì)循環(huán)和能量流動(dòng)的基礎(chǔ)(王愛(ài)愛(ài)等,2014), 在維持水生態(tài)系統(tǒng)平衡中起著重要作用。浮游植物群落結(jié)構(gòu)對(duì)水環(huán)境條件變化敏感, 其種類組成和分布格局對(duì)環(huán)境變化具有重要指示作用, 已被廣泛應(yīng)用于水質(zhì)評(píng)價(jià)中(況琪軍等, 2005; 王愛(ài)愛(ài)等,2014; 茍婷等, 2015; 楊敏等, 2015)。

沱江為長(zhǎng)江上游重要支流, 流域內(nèi)人口密度高,工農(nóng)業(yè)發(fā)達(dá), 污染嚴(yán)重, 富營(yíng)養(yǎng)化程度在四川省主要河流中最高(劉佳等, 2012)。2003年的特大水污染事件(孫莉等, 2005)引起了學(xué)者們對(duì)沱江水質(zhì)的密切關(guān)注, 從各種理化指標(biāo)的角度對(duì)沱江水質(zhì)進(jìn)行評(píng)價(jià)(謝賢健等, 2009; 肖玖金等, 2011; 周實(shí)際等, 2011; 劉建平, 2013; 陳雨艷等, 2015; 周開(kāi)錫等, 2015), 而生物學(xué)監(jiān)測(cè)方面的研究仍停留在20世紀(jì)80、90年代(李進(jìn), 1986; 楊昌述, 1987)。本研究對(duì)沱江中下游10個(gè)斷面開(kāi)展采樣調(diào)查, 從理化指標(biāo)和浮游植物兩個(gè)角度出發(fā)進(jìn)行水質(zhì)評(píng)價(jià), 為沱江水體健康狀況提供重要的基礎(chǔ)數(shù)據(jù), 也為水體的利用及河流生態(tài)修復(fù)提供參考資料。

1 材料與方法

1.1 樣品采集與分析

1.1.1 采樣點(diǎn)布設(shè) 根據(jù)研究江段周邊環(huán)境特點(diǎn),在城市區(qū)(蓮花山、資陽(yáng)一橋、西林渡口、沱橋)、工業(yè)區(qū)或建筑工地附近(麻柳壩工業(yè)園、資州大橋)、城郊上游或鄉(xiāng)鎮(zhèn)(順河場(chǎng)鎮(zhèn)、五里店水電站、銀山鎮(zhèn)、二水廠)布設(shè)采樣斷面 10個(gè)(圖 1)。各采樣點(diǎn)環(huán)境描述參考陶敏等(2016)。

圖1 采樣點(diǎn)分布圖Fig.1 Map of sampling sites

1.1.2 樣品采集 于2013年8月(豐水期)與11月(枯水期)在研究江段進(jìn)行采樣和調(diào)查。表層水樣以辛德勒采水器(Schindler-Patalas trap)采集。水化學(xué)及葉綠素樣品直接裝入干凈的樣品瓶密封, 冷藏帶回實(shí)驗(yàn)室。同時(shí)將1L浮游植物定量樣品裝入樣品瓶并立即加入魯哥氏試劑(Lugol’s reagent) 現(xiàn)場(chǎng)固定。浮游植物定性樣品用 25號(hào)浮游生物網(wǎng)在水面下作“∞”字型撈取并過(guò)濾收集, 加入4%甲醛現(xiàn)場(chǎng)固定。

1.1.3 浮游植物樣品的處理與分析 將浮游植物定量樣品在浮游生物沉淀器中避光靜置沉降 48h, 以虹吸法去除上層清液后, 濃縮定容至 50mL, 加入數(shù)滴甲醛溶液保存。計(jì)數(shù)方法參照章宗涉等(1991), 種類鑒定參考胡鴻鈞等(2006)。

1.1.4 理化指標(biāo)的分析測(cè)定 使用塞氏盤(pán)(Secchi disc)、便攜式溶解氧測(cè)量?jī)x(YSI-550A)、電導(dǎo)率儀(Mettler Toledo SevenGo SG3)等設(shè)備對(duì)研究區(qū)域透明度(Secchi Depth, SD)、水溫、pH、溶解氧(DO)和電導(dǎo)率(Conductivity)等進(jìn)行現(xiàn)場(chǎng)測(cè)定。總氮(TN)、氨氮(NH4-N)、硝酸鹽氮(NO3-N)、總磷(TP)、正磷酸鹽(PO4-P)、高錳酸鹽指數(shù)(CODMn)、葉綠素a(Chl a)含量在實(shí)驗(yàn)室中按照《水和廢水監(jiān)測(cè)分析方法》(國(guó)家環(huán)境保護(hù)總局, 2002)進(jìn)行分析。

1.2 評(píng)價(jià)方法

1.2.1 優(yōu)勢(shì)種與優(yōu)勢(shì)度 浮游植物優(yōu)勢(shì)種依據(jù)各種的優(yōu)勢(shì)度指數(shù)(Y)值進(jìn)行鑒別(劉成士等, 2015)：

式中, Pi為第i種浮游植物豐度與樣品中浮游植物總豐度的比值(ni/N), fi為第i種浮游植物在各采樣點(diǎn)出現(xiàn)的頻度。Y≥0.02的物種為優(yōu)勢(shì)種。

1.2.2 綜合營(yíng)養(yǎng)狀態(tài)指數(shù) 以總氮(TN)、總磷(TP)、高錳酸鹽指數(shù)(CODMn)、葉綠素a(Chl a)和透明度(SD)為評(píng)價(jià)指標(biāo), 采用湖泊(水庫(kù))富營(yíng)養(yǎng)化評(píng)價(jià)方法及分級(jí)技術(shù)規(guī)定(中國(guó)環(huán)境監(jiān)測(cè)總站, 2001)中的綜合營(yíng)養(yǎng)狀態(tài)指數(shù) TLI(∑)評(píng)價(jià)方法與標(biāo)準(zhǔn)對(duì)沱江營(yíng)養(yǎng)狀態(tài)進(jìn)行分析。

式中, Wj為第 j種參數(shù)的營(yíng)養(yǎng)狀態(tài)指數(shù)的相關(guān)權(quán)重;TLI(j)為第 j種參數(shù)的營(yíng)養(yǎng)狀態(tài)指數(shù)。TLI(∑)≤30時(shí)為貧營(yíng)養(yǎng)狀態(tài); 30＜TLI(∑)≤50為中營(yíng)養(yǎng)狀態(tài);TLI(∑)＞50 則為富營(yíng)養(yǎng)化狀態(tài), 其中 50＜TLI(∑)≤60為輕度富營(yíng)養(yǎng)化, 60＜TLI(∑)≤70為中度富營(yíng)養(yǎng)化,TLI(∑)＞70 為重度富營(yíng)養(yǎng)化。

1.2.3 多樣性指數(shù) 同時(shí)運(yùn)用 Shannon-Wiener多樣性指數(shù)(H′)、Margalef豐富度指數(shù)(d)、Pielou均勻度指數(shù)(J)進(jìn)行評(píng)價(jià)。計(jì)算公式分別為

式中, S為樣品中物種種類總數(shù), n為樣品中第i種浮游植物的豐度, N為樣品中浮游植物總豐度。H’＞3時(shí)為輕或無(wú)污染, 1—3為中污染, 其中2—3為β-中污,1—2為 α-中污, 在 0—1之間為重污染。d＞5時(shí)為清潔, 在4—5之間為寡污型, 3—4之間為β-中污, 1—3為α-中污, 0—1為重污染。J值在0.8—1.0之間時(shí)為清潔, 0.5—0.8為寡污型, 0.4—0.5為β-中污, 0.3—0.4為α-中污, 0—0.3為重污染(況琪軍等, 2005; 周廣杰等, 2006; 郭春燕等, 2010; Spatharis et al, 2011; 劉成士等, 2015; 楊敏等, 2015; 茍婷等, 2015)。

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

本文統(tǒng)計(jì)分析采用 SPSS 19.0與 Salford Predictive Modeler (SPM) 7.0軟件實(shí)現(xiàn); 基本圖形繪制在Origin 9.0軟件中完成。

1.3.1 分類回歸樹(shù) 為分析營(yíng)養(yǎng)狀態(tài)與生物多樣性與浮游植物密度之間的關(guān)系, 以H′、J、d、TLI(∑)以及采樣點(diǎn)類型(城市、農(nóng)村、工業(yè)區(qū))為預(yù)測(cè)指標(biāo), 采用分類回歸樹(shù)(classification and regression tree, CART)評(píng)價(jià)這5個(gè)預(yù)測(cè)指標(biāo)的對(duì)浮游植物密度的重要性。對(duì)浮游植物密度進(jìn)行常用對(duì)數(shù)(lg)轉(zhuǎn)換, 采用R2評(píng)價(jià)CART預(yù)測(cè)效果。CART分析采用Salford Predictive Modeler (SPM) 7.0。

1.3.2 聚類分析 為探尋浮游植物的空間分布差異, 根據(jù)浮游植物種類和密度, 利用SPSS 19.0軟件對(duì)不同采樣斷面的浮游植物進(jìn)行聚類分析。浮游植物密度經(jīng)過(guò)平方根轉(zhuǎn)換。

2 結(jié)果與討論

2.1 沱江環(huán)境參數(shù)與營(yíng)養(yǎng)狀態(tài)

2.1.1 主要環(huán)境參數(shù) 從季節(jié)來(lái)看, 沱江水體中的浮游植物密度、葉綠素a、TN與TP數(shù)值總體呈現(xiàn)出豐水期高于枯水期的趨勢(shì); 與之相反, 透明度則表現(xiàn)為枯水期較高(圖2、圖3)。這一現(xiàn)象可能與豐水期較大的地表徑流量帶來(lái)的豐富氮磷(夏品華等, 2011)及該時(shí)期適宜的水溫、光照有關(guān)。其余參數(shù)無(wú)明顯的季節(jié)差異。

圖2 葉綠素a、透明度、高錳酸鹽指數(shù)、總氮和總磷的空間分布Fig.2 Spatial distribution of chlorophyll a, Secchi depth (SD), CODMn, TN and TP

圖3 浮游植物種類數(shù)與密度分布Fig.3 Spatial distribution of phytoplankton species and density

從空間分布來(lái)看, 豐水期研究區(qū)域上游河段(資陽(yáng)市區(qū)) 的葉綠素 a、浮游植物密度與 TP等參數(shù)值均高于中、下游河段(資中、內(nèi)江市區(qū)), 見(jiàn)圖2、圖3。這可能與沱江上游磷礦開(kāi)采有關(guān)(施澤明等, 2012),在豐水期高氮含量的情況下, 浮游植物的生長(zhǎng)主要受到磷的限制(田時(shí)彌等, 2015), 故而上游水域因地表徑流帶來(lái)較多的磷, 浮游植物與葉綠素含量均較高?？菟竟?jié)葉綠素a、浮游植物密度與TN在中、下河段較高(圖2、圖3), 說(shuō)明枯水期由于地表徑流沖刷力度降低, 帶入沱江中的氮素逐減少并成為藻類生長(zhǎng)限制因子。在兩個(gè)季節(jié)中, 透明度變化與浮游植物密度變化規(guī)律不相符(圖2、圖3), 可能與沱江水體中泥沙含量有關(guān)。

2.1.2 綜合營(yíng)養(yǎng)狀態(tài)指數(shù) 根據(jù)綜合營(yíng)養(yǎng)狀態(tài)指數(shù) TLI(∑)的大小對(duì)研究江段水體進(jìn)行富營(yíng)養(yǎng)化狀態(tài)評(píng)價(jià)(圖4a), 沱江水質(zhì)總體處于中-富營(yíng)養(yǎng)狀態(tài)。與四川省其它主要河流如岷江、嘉陵江等相比, 沱江的TLI(∑)較高, 富營(yíng)養(yǎng)化較為嚴(yán)重(劉佳等, 2012)?？菟竟?jié)水體整體為中營(yíng)養(yǎng)狀態(tài)(＜50), 水質(zhì)良好; 豐水期多存在各種程度的污染, 尤其是位于資陽(yáng)市區(qū)的 1—3號(hào)采樣點(diǎn), 富營(yíng)養(yǎng)化程度較高(＞60)。豐水期較大的降水沖刷使得更多的氮磷營(yíng)養(yǎng)進(jìn)入水體, 尤其是較接近上游磷礦的江段, 富營(yíng)養(yǎng)化水平較高, 而中、下游水體的自凈作用在一定程度上緩解了富營(yíng)養(yǎng)化進(jìn)程。

圖4 各種參數(shù)的空間分布情況Fig.4 Spatial distribution of various indicators

2.2 浮游植物群落結(jié)構(gòu)分析與水質(zhì)評(píng)價(jià)

2.2.1 浮游植物群落組成與優(yōu)勢(shì)種 調(diào)查期間共檢出浮游植物7門(mén)88種(屬), 其中綠藻種類數(shù)最多為37種(屬), 硅藻23種(屬), 藍(lán)藻15種(屬), 其它門(mén)類13種(屬)。浮游植物種類數(shù)在空間分布上變化不大,但在時(shí)間分布上有較大差異(圖 3), 該差異主要由藍(lán)藻和綠藻引起。豐水期各斷面種類數(shù)變化范圍為30—52種, 以6#最低, 4#最高; 枯水期為12—28種,3#最低, 9#最高。根據(jù)圖3, 豐水期浮游植物密度主要由藍(lán)藻(35%)、硅藻(35%)和綠藻(24%)組成; 枯水期硅藻成為優(yōu)勢(shì)類群(44%), 隱藻數(shù)量增多(18%), 藍(lán)藻和綠藻水平相當(dāng)(15%)。浮游植物群落結(jié)構(gòu)是評(píng)價(jià)水體富營(yíng)養(yǎng)化的重要參數(shù)。一般認(rèn)為甲藻、隱藻和硅藻易在中營(yíng)養(yǎng)型水體中占優(yōu)勢(shì), 綠藻和藍(lán)藻則易在富營(yíng)養(yǎng)型水體中占優(yōu)勢(shì)(金相燦等, 1990)。沱江豐、枯水季節(jié)浮游植物群落結(jié)構(gòu)分別為綠藻+硅藻+藍(lán)藻和硅藻+綠藻+隱藻, 由此推斷該江段目前應(yīng)處于中到富營(yíng)養(yǎng)狀態(tài), 尤其是豐水期, 富營(yíng)養(yǎng)化程度更高, 這與前文營(yíng)養(yǎng)狀態(tài)分析結(jié)果一致。

研究江段浮游植物優(yōu)勢(shì)種較多且優(yōu)勢(shì)度均不高(表 1), 表明其群落結(jié)構(gòu)較復(fù)雜也較為完整。硅藻在兩個(gè)季節(jié)均有優(yōu)勢(shì), 其中梅尼小環(huán)藻(Cyclotellameneghiniana)和花環(huán)小環(huán)藻(C. operculata)存在豐水期高峰, 脆桿藻(Fragilariasp.)、線形舟形藻(Navicula graciloides)和針桿藻(Synedrasp.)則在枯水期占優(yōu)勢(shì)。藍(lán)藻和隱藻有明顯的季節(jié)差異： 藍(lán)藻為豐水期優(yōu)勢(shì)類群, 優(yōu)勢(shì)種為湖泊假魚(yú)腥藻(Pseudanabaena limnetica)、拉氏擬柱孢藻(Cylindrospermopsis raciborskii)、微小平裂藻(Merismopedia tenuissima)和細(xì)小平裂藻(M. minima); 隱藻則主要在枯水季節(jié)占優(yōu)勢(shì), 包括尖尾藍(lán)隱藻(Chroomonas acuta)、卵形隱藻(Cryptomonas ovata)和嚙蝕隱藻(C. erosa)。其它類群如綠藻和金藻優(yōu)勢(shì)種類少且優(yōu)勢(shì)度均較低。豐水期優(yōu)勢(shì)度最高的梅尼小環(huán)藻和微小平裂藻均為 α-中污指示種(沈蘊(yùn)芬等, 1990), 且富營(yíng)養(yǎng)化水體中常見(jiàn)的絲狀藍(lán)藻湖泊假魚(yú)腥藻(黃成等, 2011)等出現(xiàn)頻度亦較高, 說(shuō)明該季節(jié)沱江已處于富營(yíng)養(yǎng)化狀態(tài)?？菟竟?jié)優(yōu)勢(shì)種則反映出寡污-中污狀態(tài)(沈蘊(yùn)芬等, 1990)。

表1 沱江浮游植物優(yōu)勢(shì)種及優(yōu)勢(shì)度Tab.1 Dominant species and dominance degree of phytoplankton in Tuojiang River

2.2.2 浮游植物群落聚類分析 聚類分析(圖 5)結(jié)果顯示, 研究江段中游3個(gè)位于五里店水電站壩下的斷面(5#、6#、7#)聚集為集合 1; 下游天宮堂水電站壩上 3個(gè)斷面(8#、9#、10#)與中上游的 4#斷面共同組成集合2; 上游南津驛水電站壩上3個(gè)斷面(1#、2#、3#)則聚為集合3。3個(gè)集合的藍(lán)藻、綠藻、硅藻、隱藻比例分別為15.6%、19.4%、50%、8.7%和33.7%、20.3%、30.5%、8.7%以及21.5%、17.7%、41.4、16.9%,硅藻在壩下集合 1中所占比例最大, 藍(lán)藻+綠藻在壩上集合2中最多, 隱藻則在壩上集合3中最多。由于水電壩的攔截作用, 大壩上下游之間生態(tài)條件會(huì)發(fā)生較大變化, 壩上湖泊區(qū)由于流速和透明度的變化其浮游植物群落結(jié)構(gòu)也顯著區(qū)別與壩下區(qū)域(韓博平,2010; 楊敏等, 2015), 因此與壩下集合 1相比, 壩上集合2與集合3中湖泊型藍(lán)藻、綠藻類及隱藻增多,河流型硅藻類有所下降(邱光勝等, 2011), 這一變化趨勢(shì)與嘉陵江等鄰近水體接近(楊敏等, 2015)。

單因素方差分析(One-way ANOVA)顯示, 浮游植物密度與NH+4數(shù)值均表現(xiàn)出從集合1到集合3明顯升高的趨勢(shì)(P＜0.05, 表 2), 其余指標(biāo)差異不顯著(P＞0.05, 表 2), 說(shuō)明水電壩的建設(shè)與空間分布格局對(duì)水體多種營(yíng)養(yǎng)鹽濃度作用較小, 但對(duì)與浮游植物群落結(jié)構(gòu)已有影響。由于水電壩的攔截作用, 大壩上游附近區(qū)域流速減緩, NH+4等營(yíng)養(yǎng)積累(趙晨等, 2014), 適宜浮游植物增殖, 故壩上集合2與3中浮游植物密度與NH+4均高于壩下集合1; 同為壩上采樣點(diǎn), 較上游的集合3中NH+4濃度最高可能與沱江上游的企業(yè)污水排放有關(guān)(孫莉等, 2005),喜好NH+4的浮游植物(de Aldaet al, 1996)密度也隨之更高。

圖5 采樣點(diǎn)聚類分析樹(shù)狀圖Fig.5 Cluster analysis for the sampling sites

表2 沱江不同采樣斷面集群水體特征比較Tab.2 Comparison of biological and physiochemical parameters among sample sites in 3 groups

2.2.3 多樣性指數(shù) 沱江浮游植物多樣性指數(shù)變化趨勢(shì)如圖 4b—d。豐水期中各指數(shù)變化較為一致,物種多樣性指數(shù)H′平均3.66, 大致處于清潔狀態(tài), 最高值出現(xiàn)在 7#斷面(4.14), 最低為 6#斷面(1.98); 物種豐富度指數(shù)d平均6.21, 整體為清潔, 1#斷面最高(7.47), 6#斷面最低(4.39); 均勻度指數(shù)J平均0.69, 多為清潔-寡污狀態(tài), 與多樣性指數(shù)H′一樣, 最高與最低值分別出現(xiàn)在6#(0.84)和7#(0.40)。6#與10#等斷面物種多樣性與均勻度指數(shù)較低呈中度污染狀態(tài), 是由于細(xì)胞很小的群體藍(lán)藻, 如平裂藻屬(Merismopediaspp.)的優(yōu)勢(shì)造成。盡管該類藻細(xì)胞數(shù)量較大, 但對(duì)浮游植物生物量和水質(zhì)的影響卻較小, 形成水華的可能性也較低, 根據(jù)此結(jié)果豐水期沱江處于清潔-寡污狀態(tài)。枯水期各指數(shù)一致性較差： 多樣性指數(shù)H′(3.13—4.21)與均勻度指數(shù)J(0.80—0.93)各點(diǎn)差異不明顯; 物種豐富度指數(shù)d最低值出現(xiàn)在 3#(3.12), 最高值為 8#斷面(5.40), 多數(shù)點(diǎn)處于輕度到中度污染狀態(tài)。三個(gè)指數(shù)中, 物種豐富度指數(shù)d(圖4)與浮游植物種類數(shù)(圖 3)變化趨勢(shì)一致性較好, 說(shuō)明該指數(shù)對(duì)沱江中浮游植物群落物種數(shù)變化較為敏感, 因此按物種豐富度指數(shù)d評(píng)價(jià)標(biāo)準(zhǔn), 枯水季節(jié)沱江各斷面多處于輕度到中度污染狀態(tài)。但是應(yīng)當(dāng)看到, 盡管從多樣性來(lái)說(shuō)沱江污染程度較低, 其實(shí)營(yíng)養(yǎng)水平和浮游植物群落結(jié)構(gòu)卻呈現(xiàn)出明顯的富營(yíng)養(yǎng)化狀態(tài)。由于環(huán)境的復(fù)雜性, 從不同角度評(píng)價(jià)水體狀態(tài)均有其局限性(易斌等, 2007), 綜合各種指標(biāo)進(jìn)行判定, 沱江污染程度仍不容樂(lè)觀。

2.3 浮游植物密度與營(yíng)養(yǎng)和物種多樣性的關(guān)系

浮游植物群落結(jié)構(gòu)受外界環(huán)境(營(yíng)養(yǎng)水平)與內(nèi)部因素(種內(nèi)與種間關(guān)系, 即群落結(jié)構(gòu))的共同作用,從數(shù)據(jù)上來(lái)說(shuō), 浮游植物密度與營(yíng)養(yǎng)狀態(tài)指數(shù)以及多樣性指數(shù)有某種關(guān)聯(lián)。以浮游植物密度為因變量,綜合營(yíng)養(yǎng)狀態(tài)指數(shù)和多樣性指數(shù)為自變量所作的分類回歸樹(shù)模型R2=0.93, 說(shuō)明這兩類指數(shù)能較好地解釋浮游植物密度的變化。該模型將浮游植物密度數(shù)據(jù)作了3次細(xì)分(圖6)。第一個(gè)分支將4個(gè)物種豐富度較高(d＞6.56)并且浮游植物密度較高(未進(jìn)行 lg轉(zhuǎn)換的平均值為 4.6×106cell/L)的樣品與其它樣品分離開(kāi)來(lái); 第二個(gè)分支又將 5個(gè)具有較低營(yíng)養(yǎng)水平(TLI≤43.9)且浮游植物密度較低(平均值為0.4 ×106cell/L)的樣品從數(shù)據(jù)集中分離出來(lái); 第三個(gè)分支將剩余的 11個(gè)樣品以J=0.85為標(biāo)準(zhǔn)分為兩部分, 均勻度低(≤0.85)的 8個(gè)樣品浮游植物密度相對(duì)較高(平均值為2.6×106cell/L), 均勻度較高(＞0.85)的 3 個(gè)樣品浮游植物密度較低(平均 1.0 ×106cell/L)。Shannon-Wiener多樣性指數(shù)(H′)與采樣點(diǎn)位置類型未被納入模型中。該模型顯示, 營(yíng)養(yǎng)水平與物種豐富度較高, 但均勻度偏低的水體中, 浮游植物密度易出現(xiàn)高值。在營(yíng)養(yǎng)物質(zhì)富足的情況下, 浮游植物種類較多, 但有一部分為偶見(jiàn)種(俞秋佳, 2014), 少數(shù)優(yōu)勢(shì)種的大量繁殖降低了群落均勻度, 這些優(yōu)勢(shì)種若繼續(xù)發(fā)展, 水體就有發(fā)生水華的潛在風(fēng)險(xiǎn)(馮波等, 2008)。對(duì)漢江硅藻水華的研究表明, 較豐富的物種多樣性有可能是水華初期的特點(diǎn)(潘曉潔等, 2014)。因此在營(yíng)養(yǎng)水平較高的水體中,為預(yù)警水華的發(fā)生, 豐富的浮游植物種類不僅不能成為放松警惕的原因, 更應(yīng)成為密切監(jiān)測(cè)的理由。

圖6 浮游植物密度與綜合營(yíng)養(yǎng)狀態(tài)指數(shù)、多樣性指數(shù)間的關(guān)系Fig.6 Classification and regression tree of phytoplankton abundance versus trophic level index and diversity index x為未進(jìn)行l(wèi)g轉(zhuǎn)換的浮游植物密度平均值; n為樣本數(shù)量

3 結(jié)論

根據(jù)綜合營(yíng)養(yǎng)狀態(tài)指數(shù)判斷沱江水體處于中-富營(yíng)養(yǎng)狀態(tài), 富營(yíng)養(yǎng)化程度較高。由于雨水淋溶沖刷作用和上游磷礦的影響, 豐水期污染程度高于枯水期,上游斷面高于較下游斷面。

從浮游植物群落組成來(lái)看, 沱江為中-富營(yíng)養(yǎng)型水體, 其中豐水期為綠藻+硅藻+藍(lán)藻的形式, 枯水期為硅藻+綠藻+隱藻的形式, 仍表現(xiàn)出豐水期富營(yíng)養(yǎng)化程度高于枯水期的狀態(tài)。優(yōu)勢(shì)種類的組成顯示, 沱江中浮游植物群落結(jié)構(gòu)較為復(fù)雜和完整, 豐水期為α-中污染, 枯水期則為寡污-中污狀態(tài)。由于沿江梯級(jí)大壩的影響, 沱江水文條件與 NH4+分布具有明顯空間異質(zhì)性, 導(dǎo)致了大壩上下游之間、上游斷面與下游斷面之間浮游植物群落結(jié)構(gòu)與密度的差異, 整體表現(xiàn)出大壩上游湖泊型藻類數(shù)量增加、上游斷面浮游植物密度高于下游斷面的趨勢(shì)。盡管多樣性指數(shù)顯示沱江污染狀態(tài)較輕, 綜合考慮各類指標(biāo)的局限性, 沱江污染程度仍不容樂(lè)觀。

沱江中浮游植物密度受到營(yíng)養(yǎng)和群落結(jié)構(gòu)的顯著影響： 在營(yíng)養(yǎng)水平高, 物種豐富度較大但均勻度較低的情況下, 容易暴發(fā)水華。因此對(duì)于營(yíng)養(yǎng)物質(zhì)含量豐富并且浮游植物種類數(shù)較多的水體, 不能放松警惕, 應(yīng)加強(qiáng)監(jiān)控。

王愛(ài)愛(ài), 馮 佳, 謝樹(shù)蓮, 2014. 汾河中下游浮游藻類群落特征及水質(zhì)分析. 環(huán)境科學(xué), 5(3)： 915—923

中國(guó)環(huán)境監(jiān)測(cè)總站, 2001. 湖泊(水庫(kù))富營(yíng)養(yǎng)化評(píng)價(jià)方法及分級(jí)技術(shù)規(guī)定. 北京： 中國(guó)環(huán)境科學(xué)出版社

田時(shí)彌, 楊 揚(yáng), 喬永民等, 2015. 珠江流域東江干流浮游植物葉綠素 a時(shí)空分布及與環(huán)境因子的關(guān)系. 湖泊科學(xué),27(1)： 31—37

馮 波, 陳俊輝, 盧伙勝, 2008. 茂名市大放雞島海域浮游植物群落的季節(jié)變化. 海洋漁業(yè), 30(1)： 67—73

劉 佳, 徐 亮, 張秋勁等, 2012. 四川省水體富營(yíng)養(yǎng)狀況時(shí)空分布調(diào)查. 中國(guó)環(huán)境監(jiān)測(cè), 28(5)： 6—8

劉成士, 黃 燕, 劉 瑜等, 2015. 青山水庫(kù)輪蟲(chóng)群落結(jié)構(gòu)特征的時(shí)空差異對(duì)水質(zhì)變化的影響. 海洋與湖沼, 46(3)：694—702

劉建平, 2013. 沱江富順段水質(zhì)分析及污染狀況研究. 四川環(huán)境, 32(2)： 23—26

孫 莉, 朱鴻斌, 張成云等, 2005. 一起工業(yè)廢水污染沱江水源水事故的調(diào)查. 環(huán)境與健康雜志, 22(3)： 192—193

李 進(jìn), 1986. 沱江水質(zhì)污染的微生物監(jiān)測(cè)評(píng)價(jià). 四川環(huán)境,5(3)： 26—32

楊 敏, 張 晟, 劉朔孺, 2015. 草街水庫(kù)蓄水后嘉陵江浮游植物群落特征及水質(zhì)評(píng)價(jià). 環(huán)境科學(xué), 36(7)： 2480—2486

楊昌述, 1987. 沱江干流底棲動(dòng)物調(diào)查及水質(zhì)生物學(xué)評(píng)價(jià). 環(huán)境科學(xué), 8(5)： 50—56

肖玖金, 譚周亮, 李旭東等, 2011. 基于物元分析法的沱江地表水水質(zhì)評(píng)價(jià). 安全與環(huán)境學(xué)報(bào), 11(1)： 104—108

邱光勝, 胡 圣, 葉 丹等, 2011. 三峽庫(kù)區(qū)支流富營(yíng)養(yǎng)化及水華現(xiàn)狀研究. 長(zhǎng)江流域資源與環(huán)境, 20(3)： 311—316

況琪軍, 馬沛明, 胡征宇等, 2005. 湖泊富營(yíng)養(yǎng)化的藻類生物學(xué)評(píng)價(jià)與治理研究進(jìn)展. 安全與環(huán)境學(xué)報(bào), 5(2)： 87—91

沈韞芬, 章宗涉, 龔循矩等, 1990. 微型生物監(jiān)測(cè)新技術(shù). 北京： 中國(guó)建筑工業(yè)出版社, 1—524

陳雨艷, 余 恒, 向秋實(shí)等, 2015. 沱江流域水環(huán)境質(zhì)量分析.四川環(huán)境, 34(2)： 85—89

茍 婷, 馬千里, 許振成等, 2015. 貝江浮游藻類群落特征及富營(yíng)養(yǎng)化風(fēng)險(xiǎn)分析. 環(huán)境科學(xué), 36(3)： 946—954

國(guó)家環(huán)境保護(hù)總局, 2002. 水和廢水監(jiān)測(cè)分析方法(第四版).北京： 中國(guó)環(huán)境科學(xué)出版社

易 斌, 吳玲玲, 林 端等, 2007. 營(yíng)養(yǎng)指數(shù)法和浮游植物多樣性指數(shù)法在評(píng)價(jià)富營(yíng)養(yǎng)化上的差異. 見(jiàn)： 中國(guó)海洋學(xué)會(huì)赤潮研究與防治專業(yè)委員會(huì)第二屆學(xué)術(shù)研討會(huì). 上海： 中國(guó)海洋學(xué)會(huì)

金相燦, 屠清瑛, 1990. 湖泊富營(yíng)養(yǎng)化調(diào)查規(guī)范(第二版). 北京： 中國(guó)環(huán)境科學(xué)出版社

周廣杰, 況琪軍, 胡征宇等, 2006. 三峽庫(kù)區(qū)四條支流藻類多樣性評(píng)價(jià)及“水華”防治. 中國(guó)環(huán)境科學(xué), 26(3)： 337—341

周開(kāi)錫, 李 敏, 吳 建, 2015. 沱江河內(nèi)江段主要監(jiān)控?cái)嗝嫠|(zhì)變化趨勢(shì)分析. 見(jiàn)： 2015年中國(guó)環(huán)境科學(xué)學(xué)會(huì)學(xué)術(shù)年會(huì)論文集(第一卷). 深圳： 中國(guó)環(huán)境科學(xué)學(xué)會(huì)

周實(shí)際, 朱 俊, 張 科等, 2011. 沱江內(nèi)江城區(qū)段水質(zhì)分析.環(huán)境科學(xué)與管理, 36(10)： 114—116

趙 晨, 董世魁, 劉世梁等, 2014. 漫灣大壩上下游沉積物重金屬與營(yíng)養(yǎng)元素分布特征及環(huán)境風(fēng)險(xiǎn)評(píng)價(jià). 環(huán)境科學(xué)學(xué)報(bào), 34(9)： 2417—2425

胡鴻鈞, 魏印心, 2006. 中國(guó)淡水藻類—系統(tǒng)、分類及生態(tài). 北京： 科學(xué)出版社, 23—915

俞秋佳, 2014. 蘇州河浮游植物群落變化特征及其影響因子研究. 上海： 華東師范大學(xué)碩士學(xué)位論文

施澤明, 倪師軍, 張成江等, 2012. 沱江流域磷礦開(kāi)采和加工過(guò)程中放射性環(huán)境問(wèn)題探討. 地球科學(xué)進(jìn)展, 27(10)：1134—1139

夏品華, 李秋華, 林 陶等, 2011. 貴州高原百花湖水庫(kù)湖沼學(xué)變量特征及環(huán)境效應(yīng). 環(huán)境科學(xué)學(xué)報(bào), 31(8)：1660—1669

郭春燕, 馮 佳, 謝樹(shù)蓮, 2010. 山西晉陽(yáng)湖浮游藻類分布的時(shí)空格局及水質(zhì)分析. 湖泊科學(xué), 22(2)： 251—255

陶 敏, 王永明, 謝碧文等, 2016. 沱江浮游生物群落時(shí)空分布及相關(guān)環(huán)境因子分析. 水生生物學(xué)報(bào), 40(2)： 301—312

黃 成, 侯 偉, 顧繼光等, 2011. 珠江三角洲城市周邊典型中小型水庫(kù)富營(yíng)養(yǎng)化與藍(lán)藻種群動(dòng)態(tài). 應(yīng)用與環(huán)境生物學(xué)報(bào), 17(3)： 295—302

章宗涉, 黃祥飛, 1991. 淡水浮游生物研究方法. 北京： 科學(xué)出版社

韓博平, 2010. 中國(guó)水庫(kù)生態(tài)學(xué)研究的回顧與展望. 湖泊科學(xué),22(2)： 151—160

焦珂?zhèn)? 周啟星, 2015. 基于水質(zhì)與生物指標(biāo)的松花江流域水生態(tài)健康評(píng)價(jià). 生態(tài)學(xué)雜志, 34(6)： 1731—1737

謝賢健, 蘭代萍, 2009. 基于因子分析法的沱江流域地表水水質(zhì)的綜合評(píng)價(jià). 安徽農(nóng)業(yè)科學(xué), 37(3)： 1304—1306

潘曉潔, 朱愛(ài)民, 鄭志偉等, 2014. 漢江中下游春季浮游植物群落結(jié)構(gòu)特征及其影響因素. 生態(tài)學(xué)雜志, 33(1)： 33—40

Bozzetti M, Schulz U H, 2004. An index of biotic integrity based on fish assemblages for subtropical streams in southern Brazil. Hydrobiologia, 529(1)： 133—144

de Alda J A G O, Tapia M I, Franck F et al, 1996. Changes in nitrogen source modify distribution of excitation energy in the cyanobacterium Phormidium laminosum. Physiologia Plantarum, 97(1)： 69—78

Spatharis S, Roelke D L, Dimitrakopoulos P G et al, 2011.Analyzing the (mis) behavior of Shannon index in eutrophication studies using field and simulated phytoplankton assemblages. Ecological Indicators, 11(2)： 697—703