陳 瑜， 景 明， 高 昕， 李繼影

(江蘇省蘇州環(huán)境監(jiān)測中心, 江蘇 蘇州 215004)

0 引言

太湖（北緯30°55′40" ～31°32′58"和東經(jīng)119°52′32" ～120°36′10"）橫跨江、浙2 省，位于長江三角洲的南緣，是中國第三大淡水湖。太湖的主要進(jìn)出河流有50 余條[1]，環(huán)太湖出入湖水量是對太湖水量變化影響最大的一個(gè)因子[2]。 所處的太湖流域?yàn)榈湫偷钠皆泳W(wǎng)地區(qū)，河道總長約12 萬km， 密度達(dá)3.3 km/km2，0.5 km2以上的大小湖泊189 個(gè)，多年平均水資源量176 億m3，人均本地水資源量僅為全國平均水平的1/5[3]。 太湖在給人們帶來物質(zhì)文化財(cái)富的同時(shí)，也承受著人類活動帶來的污染，其中太湖的主要外源污染是由河流輸入[4-5]。 水利部太湖流域管理局自2002年起實(shí)施引江濟(jì)太調(diào)水工程，通過望虞河將長江水引入太湖，“引江濟(jì)太” 工程使得太湖換水周期縮短，水質(zhì)改善，但也為太湖環(huán)境治理帶來了新挑戰(zhàn)[6]。 “引江濟(jì)太”一定程度上抬高了太湖TP 濃度；入湖河流TP 與考核目標(biāo)有較大差距，入湖負(fù)荷仍處于高位[7]。 毛新偉[8]的研究表明引水期長江及入湖的TN 和TP 濃度基本高于貢湖灣。

硅藻( Diatoms) 是一種真核單細(xì)胞微藻類。作為水體中重要的初級生產(chǎn)者[9]，硅藻在水生生態(tài)系統(tǒng)和全球碳循環(huán)中具有舉足輕重的地位[10]。 研究發(fā)現(xiàn)，硅藻對環(huán)境變化尤其是有機(jī)污染和富營養(yǎng)化反應(yīng)很快，對水環(huán)境因子變化非常敏感[11]，溫度、營養(yǎng)鹽、浮游動物捕食等環(huán)境的變化都會使硅藻的種類及數(shù)量發(fā)生相應(yīng)的改變[12-13]。硅藻廣泛存在于不同營養(yǎng)類型水體中，硅藻與環(huán)境因子相互影響，多項(xiàng)研究表明，二者之間的相互作用在水生生態(tài)系統(tǒng)中產(chǎn)生了重要影響。 2000年歐盟水框架指導(dǎo)委員會(Water Framework Directive of the European Union) 將硅藻推薦為當(dāng)前水環(huán)境整治決策中可以用來確定營養(yǎng)水平的一個(gè)有效的生物指標(biāo)[14]。 因其對環(huán)境的指示作用比藍(lán)藻、綠藻更靈敏[15-17]，作為良好的水生態(tài)環(huán)境指示生物[18]，硅藻已被國內(nèi)外廣泛應(yīng)用于監(jiān)測湖泊和河流水質(zhì)的變化[19-21]。

本研究擬通過對太湖主要出入湖河道及太湖湖體（蘇州轄區(qū)）15 個(gè)點(diǎn)位的著生硅藻和水環(huán)境因子進(jìn)行調(diào)查研究，探究硅藻與環(huán)境的關(guān)系，期望能夠?yàn)樘h(huán)境監(jiān)測和治理以及水源地保護(hù)提供技術(shù)支持。

1 材料與方法

1.1 采樣點(diǎn)位

太湖的主要進(jìn)出河流眾多， 其中望虞河是太湖最大入湖河道之一， 望虞河是引江濟(jì)太調(diào)水的重要工程。 太湖在蘇州轄區(qū)的出湖河道有木光河、胥江、太浦河和蘇東河等。 本研究在出入湖河道及太湖水域共設(shè)置15 個(gè)采樣點(diǎn)，見圖1。

圖1 太湖東部水域采樣點(diǎn)分布

1.2 樣品采集

1.2.1 水質(zhì)檢測

現(xiàn)場使用便攜式水質(zhì)分析儀 （YSI Professional Plus,Yellow Springs,OH,USA）測定pH 值、水溫（T）、溶解氧（DO）、濁度（NTU）和電導(dǎo)率（Cond）等；透明度（SD）用塞式透明度盤測定。 實(shí)驗(yàn)室測定水樣用采水器采集水下50 cm 及湖底以上50 cm 處水樣，混勻取1 L，密封立即運(yùn)回實(shí)驗(yàn)室分析測定總氮（TN）、總磷（TP）、氨氮（NH3-N）、五日生化需氧量（BOD5）、高錳酸鹽指數(shù)（CODMn）、化學(xué)需氧量（CODCr）以及葉綠素a（Chla）等指標(biāo)，樣品保存及測定方法均參考《水和廢水監(jiān)測分析方法（第4 版）》(國家環(huán)境保護(hù)總局，2002)[22]。

1.2.2 著生硅藻采樣

將硅藻計(jì)（菲迪環(huán)境GZJ1-10）放置于采樣點(diǎn)1個(gè)月后取樣， 用毛刷或硬膠皮將基質(zhì)上所著生的藻類及其他生物，全部刮到盛有蒸餾水的玻璃瓶中，并用蒸餾水將基質(zhì)沖洗多次，用魯哥氏液固定，定容至30 mL 備用。 采樣方法及樣品保存方法均參考《水和廢水監(jiān)測分析方法（第4 版）》。

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

根據(jù)著浮游硅藻多度數(shù)據(jù)計(jì)算各采樣點(diǎn)Shannon-Wiener 多樣性指數(shù)H′， 采用Mcnaughton優(yōu)勢度指數(shù)來判定優(yōu)勢種的組成，選取Mcnaughton優(yōu)勢度指數(shù)（Y ＞0.02）的藻類為優(yōu)勢種[23]。

Shannon-Wiener 生物多樣性指數(shù)按下式計(jì)算：

式中：s 為樣品中的種類數(shù)；ni為樣品中第i 種生物的個(gè)體數(shù)；n 為樣品中生物總個(gè)體數(shù)；ni/n 為第i 種藻類的細(xì)胞數(shù)占所有藻類總細(xì)胞數(shù)的比值；fi為第i種藻類在樣點(diǎn)中出現(xiàn)的頻率。

根據(jù)PETR 等[24-26]對Canoco5 軟件使用方法的說明，將藻類數(shù)據(jù)和水質(zhì)數(shù)據(jù)導(dǎo)入Canoco 5 進(jìn)行初步分析， 保留可以進(jìn)一步分析的水質(zhì)數(shù)據(jù)（forward selection results）進(jìn)行約束（constrained）分析。根據(jù)項(xiàng)目效應(yīng)（term effects）結(jié)果，綜合簡單效應(yīng)結(jié)果（simple effects）和附條件的效應(yīng)結(jié)果（conditional effects）選取水質(zhì)因子進(jìn)入RDA 分析。 進(jìn)行冗余分析（Redundancy Analysis，RDA)分析時(shí)，除pH 值外的所有水體理化數(shù)據(jù)和著生藻類相對多度數(shù)據(jù)均進(jìn)行數(shù)據(jù)對數(shù)轉(zhuǎn)換[lg(x+1)]。以上分析在Canoco5 上進(jìn)行，浮游藻類群落結(jié)構(gòu)圖ArcMap 10.2.2 上完成.

2 著生硅藻群落結(jié)構(gòu)分析

太湖東部水體出入湖河道及太湖水體中共鑒定出著生藻類4 門63 屬， 其中硅藻門32 屬80 種，藍(lán)藻門9 屬，裸藻門1 屬，綠藻門21 屬。硅藻群落中分類單元最多的種屬為舟形藻屬（Navicula）10 種、異極藻屬（Gomphonema）12 種、菱形藻屬（Nitzschia）13種。 優(yōu)勢種為兩棲菱形藻（Nitzschia amphibia）、具球異極藻（Gomphonema sphaerophorum）、線形舟形藻（Navicula graciloides）、谷皮菱形藻（Nitzschia palea）、奇異菱形藻（Nitzschia paradoxa）、 胡斯特橋彎藻（Cymbella hustedtii）、窄異極藻（Gomphonema angustatum）、顆粒直鏈藻（Melosira granulata）、絲葉菱形藻（Nitzschia filiformis）、鈍脆桿藻（Fragilaria brevistriata）、梅尼小環(huán)藻（Cyclotella meneghiniana）。

河道水體中的物種豐富度比湖體高， 物種數(shù)最多的點(diǎn)位為T13，T15 和T4， 分別為35 種、27 種和26 種；豐富度最少的點(diǎn)位為T6 和T11，均為11 種。著生硅藻密度沒有明顯的規(guī)律可循， 藻類密度平均值為6.35×105個(gè)/cm2, 大于107個(gè)/cm2有3 個(gè)，T5，T13 和T15； 小 于104個(gè)/cm2有2 個(gè)，T9 和T11。Shannon-Wiener 指數(shù)均值為2.87，入湖河道及湖體的Shannon-Wiener 指數(shù)值普遍大于3，而出湖河道的Shannon-Wiener 指數(shù)普遍在2 左右。 見表1。

表1 各采樣點(diǎn)位的Shannon-Wiener 指數(shù)、藻類密度和物種豐富度

3 著生藻類聚類分析

將處理后的藻密度數(shù)據(jù)采用Ward 方法，以Euclidean 距離對采樣點(diǎn)位進(jìn)行聚類分析[27]，分析結(jié)果見圖2。 組1 包含6 個(gè)點(diǎn)位，組2 包含9 個(gè)點(diǎn)位。組1 的點(diǎn)位多為湖體點(diǎn)位和距離湖體較近的出湖河道點(diǎn)位，優(yōu)勢種為胡斯特橋彎藻（Cymbella hustedtii）、顆粒直鏈藻（Melosira granulata）、肘狀針桿藻（Synedra ulna）、具球異極藻（Gomphonema sphaerophorum）、鈍脆桿藻（Fragilaria brevistriata）、尖針桿藻（Synedra acus）、新月形橋彎藻（Cymbella cymbiformis）； 組2 的點(diǎn)位多為入湖河道點(diǎn)位和距離湖體較遠(yuǎn)的出湖河道點(diǎn)位，優(yōu)勢種為兩棲菱形藻（Nitzschia amphibia）、具球 異極藻（Gomphonema sphaerophorum）、線形舟形藻（Navicula graciloides）、谷皮菱形藻（Nitzschia palea）、奇異菱形藻（Nitzschia paradoxa）、胡斯特橋彎藻 （Cymbella hustedtii）、 窄異極藻（Gomphonema angustatum）、 顆粒直鏈藻（Melosira granulata）、絲葉菱形藻（Nitzschia filifor-mis）、梅尼小環(huán)藻（Cyclotella meneghiniana）。 組1 的豐富度均值為13， 比組2 小9；2 組的Shannon-Wiener 指數(shù)均值差別不大分別為2.88 和2.86；組2 的藻密度均值約為組1 的5 倍。

圖2 根據(jù)藻類種類和密度進(jìn)行采樣點(diǎn)聚類分析

4 著生硅藻生態(tài)指示型分析

VAN D 等[28]通過分析硅藻生存環(huán)境中的pH 值、氧含量、鹽度、水生環(huán)境腐殖度和營養(yǎng)狀況建立了淡水硅藻生態(tài)指示值清單。按照指示清單，將2 組硅藻所對應(yīng)的需氧量、 水生環(huán)境腐殖度和營養(yǎng)狀況這3個(gè)生態(tài)指示類型進(jìn)行統(tǒng)計(jì)， 計(jì)算了2 個(gè)分組不同生態(tài)類型物種的占比情況。 見表2（硅藻生態(tài)型分類及數(shù)據(jù)計(jì)算由OMNIDIA 6 完成）。 由表2 可知，從組1的硅藻類型主要為極度需氧型和需氧型， 組2 中的硅藻主要為中度需氧型和需氧型； 水生環(huán)境腐殖度指示結(jié)果表明，組1 硅藻主要是寡污型，其次是α-中度腐生型和β-中度腐生型， 組2 主要是α-中度腐生型和β-中度腐生型；組1 硅藻所指示的營養(yǎng)狀況從寡營養(yǎng)型至富營養(yǎng)型均有分布， 富營養(yǎng)型相對較多，組2 所指示的營養(yǎng)狀況則主要為富營養(yǎng)型。

表2 各分組藻類生態(tài)型

5 環(huán)境因子分析

RDA 分析結(jié)果表明，影響藻類群落的主要環(huán)境因子為水溫，TN，DO，CODCr，TP 和pH 值。 以上因子在個(gè)點(diǎn)位的分布見圖3。 橫坐標(biāo)為點(diǎn)位信息，縱坐標(biāo)為水質(zhì)數(shù)據(jù)值，其中TP 值為實(shí)驗(yàn)測得結(jié)果乘以20。 將水質(zhì)數(shù)據(jù)和著生藻類聚類結(jié)果進(jìn)行綜合分析，組1 和組2 的水溫，pH 值和CODCr均值差別不大；組1 的DO 均值比組2 高1.56；組2 的TN 均值是組1 的2 倍，TP 是組1 的3 倍。

圖3 各采樣點(diǎn)位重要水質(zhì)因子狀況

6 環(huán)境因子與硅藻群落RDA 分析

根據(jù)Term Effects 篩選出影響15 個(gè)采樣點(diǎn)位著生藻類分布的環(huán)境因子為DO，TP，TN，CDOCr， 見圖4。 它們的解釋度和P 值依次為DO（18.0%，0.004），TP（16.7%，0.016），TN（16.2%，0.01），CODCr（9.0%，0.184）；影響組1 著生藻類分布的環(huán)境因子為水溫，pH 值 和DO， 見 圖5。 解釋度和P值依次為T（38.5%，0.012），pH 值（19.2%，0.256），DO （17.9%，0.32）；影響組2 著生藻類分布的環(huán)境因子為CODCr，TN，TP，見圖6。解釋度和P 值依次為CODCr（20.2%，0.042），TN（19.8%，0.027），TP（11.0%，0.348）。 結(jié)合環(huán)境因子分析和圖4 可知全流域的TP，TN 及CODCr處于較高水平，組1 中的點(diǎn)位T5，T6，T7 和T11 與DO 相關(guān)性較大，組2 中的點(diǎn)位和TP，TN 及CODCr的相關(guān)性較大（本文RDA 分析根據(jù)Term Effects 篩選，首先將所有的環(huán)境因子和每個(gè)點(diǎn)位的藻類進(jìn)行分析，再從中選取解釋度和P 值高的即為影響藻類分布的主要環(huán)境因子）。

圖4 15 個(gè)采樣點(diǎn)位藻類與環(huán)境因子RDA 分析

圖5 組1（遠(yuǎn)離湖體點(diǎn)位）藻類與環(huán)境因子RDA 分析

圖6 組2（近湖體點(diǎn)位）藻類與環(huán)境因子RDA 分析

7 討論

太湖的出入河道眾多，其中望虞河作為“引江濟(jì)太”的主要輸水通道,對太湖主要區(qū)域水體交換的影響起著至關(guān)重要的作用[29]。 入湖河道增加了太湖的流動性、減少了太湖的換水周期，也給太湖環(huán)境治理帶來了新挑戰(zhàn)。 朱偉等[30]認(rèn)為河水中大量營養(yǎng)鹽進(jìn)入太湖，磷大多數(shù)以顆粒態(tài)存在并沉積在湖泊中，在適宜條件下釋放入水體。

本次的研究表明， 影響太湖主要出入湖河道著生硅藻群落的環(huán)境因子為DO，TP，TN，CODCr。 基于著生硅藻的群落結(jié)構(gòu)特征， 將太湖河道著生藻類群落劃分為2 個(gè)類群。 組1 的點(diǎn)位多為湖體點(diǎn)位和距離湖體較近的出湖河道點(diǎn)位， 組2 的點(diǎn)位多為入湖河道點(diǎn)位和距離湖體較遠(yuǎn)的出湖河道點(diǎn)位。

對2 組著生硅藻的生態(tài)指示型進(jìn)行分析， 組1硅藻所指示的營養(yǎng)狀況從寡營養(yǎng)型至富營養(yǎng)型均有分布， 組2 所指示的營養(yǎng)狀況則主要為富營養(yǎng)型。RDA 結(jié)果表明影響組1 和組2 著生藻類分布的環(huán)境因子分別為水溫，DO，pH 值，CODCr，TN 和TP。 水質(zhì)指標(biāo)也印證了硅藻的分組結(jié)果， 組1 的DO 均值優(yōu)于組2；組2 的TN 和TP 均值是組1 的2 倍以上。黃漪平[2]的研究表明太湖主要的外源污染來自河流，人們生產(chǎn)生活產(chǎn)生的污染物質(zhì)經(jīng)由河道輸入到湖泊中。 本研究的結(jié)果顯示入湖河道的CODCr，TN，TP 等水質(zhì)因子比湖體中高， 近湖體的河道水質(zhì)與湖體相差不大， 而離湖體較遠(yuǎn)的河道水質(zhì)因子比湖體低也表明太湖具有一定的自我凈化能力，。

由于受限于采樣次數(shù)和采樣點(diǎn)位量， 本研究難免有不足之處，在今后的工作中不斷收集、匯總并分析，以期能夠更真實(shí)、更準(zhǔn)確的反映硅藻和水質(zhì)的關(guān)系。