柳江流域不同水文季節(jié)水質(zhì)對附生硅藻群落的影響研究

李佳靜，劉威，鄧培雁，

1.華南師范大學(xué)化學(xué)與環(huán)境學(xué)院，廣州 510631

2.廣州大學(xué)環(huán)境科學(xué)與工程學(xué)院，廣州 510405

0 前言

硅藻是河流生態(tài)系統(tǒng)的重要初級生產(chǎn)者，在生態(tài)系統(tǒng)的物質(zhì)循環(huán)和能量流動中起著重要作用[1]。相對于底棲動物和魚類，硅藻具有種類繁多、生命周期短、對水環(huán)境變化敏感等特征，其群落結(jié)構(gòu)變化能較好地反映各種污染物對河流累積產(chǎn)生的綜合生態(tài)效應(yīng)[2]，常被廣泛用作水質(zhì)監(jiān)測和水生態(tài)健康評價的有效指示生物[3]。

水質(zhì)、地理因子、人為干擾等是影響硅藻群落形成時空分布特征的重要因素，而水質(zhì)被認(rèn)為是最主要的影響因素[4-6]，這是因為在流域尺度上，地理因子或者人為干擾都以水質(zhì)為介質(zhì)對附生硅藻群落起作用[7]。然而，河流作為一個動態(tài)系統(tǒng)，其水質(zhì)理化特征往往因人為活動和地形、降雨等自然因素的驅(qū)動而存在不同的時空差異性[8]。不同硅藻對水體中營養(yǎng)物、有機(jī)污染物和重金屬等的敏感程度和耐受能力不同，其群落組成在不同的理化條件下也會產(chǎn)生較大差異[9-10]。劉麟菲等[1]報道了渭河水系硅藻群落特征及其水質(zhì)影響因子在枯水期和豐水期具有顯著差異。Dalu 等[6]也強(qiáng)調(diào)了水文季節(jié)變化下的水質(zhì)因子對硅藻種類組成，群落演替的作用程度不同。因此，闡明水質(zhì)動態(tài)變化對硅藻群落時空變化的影響，才可更好地利用硅藻進(jìn)行河流的生態(tài)監(jiān)測與管理。

柳江流域是珠江流域西江水系的第二大支流，干流全長 751 km，面積 58520 km2，發(fā)源于貴州省獨山縣，跨越黔、桂、湘三省，是流域周圍居民生活和工農(nóng)業(yè)生產(chǎn)的重要水源[11]。同時，該流域?qū)儆谟暝葱秃恿鳎芗撅L(fēng)影響顯著，4—8 月為雨季，降水占全年總降水量的 72%，與其余月份具有明顯的豐、枯水期變化特征[12]。目前已有學(xué)者對柳江流域支流貝江[13]和龍江[14]的硅藻群落特征進(jìn)行了研究，也有學(xué)者基于附生硅藻對柳江河流構(gòu)建了水質(zhì)多指標(biāo)的評價體系[3]，但針對柳江流域硅藻群落結(jié)構(gòu)的時空分布特征及其與水質(zhì)因子關(guān)系的研究尚未見報道。為此，本文運(yùn)用聚類分析、主成分分析和冗余分析等方法，研究了水質(zhì)因子在不同水文季節(jié)變化下的硅藻優(yōu)勢種組成及群落結(jié)構(gòu)的時空分布特征，并闡明水質(zhì)動態(tài)變化對河流附生硅藻的影響，旨在為深入研究柳江流域硅藻群落的分布特征奠定基礎(chǔ)，以期為柳江流域水質(zhì)生物監(jiān)測評價及水質(zhì)的利用和管理提供理論依據(jù)。

1 材料和方法

1.1 采樣時間與位置

本研究于 2017 年 3月(枯水期)和 6 月(豐水期)分別對柳江流域的附生硅藻和水質(zhì)進(jìn)行采樣。根據(jù)柳江流域不同河段的地形地貌和環(huán)境狀況等特征，共設(shè)置 19 個采樣點(S1—S19)，覆蓋柳江流域的干流及支流，各采樣點地理位置分布如圖1所示。

1.2 采樣與分析方法

根據(jù) EN 13946 和 EN 14407 歐盟標(biāo)準(zhǔn)方法[13]，并結(jié)合采樣點實際情況進(jìn)行采樣，為減小誤差，確保采樣的代表性和準(zhǔn)確性，每個樣點采集 3 個平行樣本。附生硅藻采樣基質(zhì)選取能抵抗水流、地勢開闊處無樹蔭遮擋的石頭，在每個采樣點至少采集 5 塊石頭，用干凈的硬牙刷在石塊向陽面刷取硅藻，經(jīng)自來水沖洗至不銹鋼托盤，并收集刷洗下的硅藻樣品裝入廣口塑料瓶中，加入 5%甲醛溶液固定保存。

圖1 柳江流域采樣點分布Figure 1 Location of the sampling sites in Liujiang River Basin

硅藻樣品經(jīng) 30%的過氧化氫溶液加熱后消解，靜置去除上清液，用濃鹽酸進(jìn)行酸化，洗滌離心后制成封片。在光學(xué)顯微鏡下采用 1000 倍的油鏡鏡檢，視野內(nèi)所有完整及破損面積不超過 1/4 的硅藻細(xì)胞都要鑒定和計數(shù)，每個玻片鑒定的總數(shù)大于400 個，鑒定具體到硅藻的種等級。硅藻種類鑒定主要根據(jù) Krammer & Lange-Bertalot 鑒定體系(1986—1991)[15]開展。

水樣運(yùn)送回實驗室后，根據(jù)《水和廢水監(jiān)測分析方法(第 4 版)》[16]測定高錳酸鹽指數(shù)(CODMn)、氨氮(NH4+-N)、砷(As)、汞(Hg)、鎘(Cd)、鉛(Pb)等水質(zhì)指標(biāo)，水體溶解氧(DO)用多參數(shù)水質(zhì)分析儀(YSI Pro 2000 美國)現(xiàn)場測定。

1.3 數(shù)據(jù)處理與分析

Czekanowsiki 相似系數(shù)計算公式為:

優(yōu)勢種以優(yōu)勢度指數(shù)(Y＞0.02)判斷，其計算公式為:

式中，A、B分別為枯、豐水期水體中的硅藻種類數(shù)，C為 2 個水期水體中共有的種類數(shù)，fi為第i種的出現(xiàn)頻率，Pi為第i種的個體數(shù)與樣品中總個體數(shù)的比值。

選擇至少在2個采樣點中出現(xiàn)且相對豐度超過1%的硅藻種，其種類拉丁文名的縮寫見附錄。運(yùn)用主成分分析(PCA)篩選出枯、豐水期對硅藻群落影響較為顯著的水質(zhì)驅(qū)動因子。然后對硅藻相對多度進(jìn)行降趨勢對應(yīng)分析(Detrended correspondence analysis，DCA)，用物種單峰影響值梯度來判斷是線型分布還是單峰型分布，以判定影響硅藻群落結(jié)構(gòu)分布的水質(zhì)因子。單因素方差分析和聚類分析在IBM SPSS 19.0 軟件中處理，PCA 和 DCA 分析等在 Canoco 4.5 軟件中完成。

2 結(jié)果與分析

2.1 硅藻群落種類組成

共鑒定出附生硅藻137種，隸屬于6目9科34屬，舟形藻屬(Navicula)和曲殼藻屬(Achnanthes)在各采樣點中均有出現(xiàn)，為該流域的主要優(yōu)勢藻屬。枯、豐水期硅藻種數(shù)分別為 99 種和 102 種，種類相似系數(shù)為 63.68%。枯水期各采樣點硅藻豐富度的范圍為 13—35 種，平均值為 21 種，豐水期各采樣點硅藻豐富度的范圍為 16—42 種，平均值為31 種(圖2)。在水文季節(jié)變化上，表現(xiàn)為豐水期明顯高于枯水期。

2.2 硅藻優(yōu)勢種

圖2 各采樣點硅藻種類豐富度Figure 2 Species richness of diatoms at each sampling site

兩次調(diào)查共發(fā)現(xiàn)硅藻優(yōu)勢種為 11 種，其中枯水期有5種，豐水期有8種，優(yōu)勢度指數(shù)為0.021—0.363(表1)。硅藻優(yōu)勢種組成隨水文季節(jié)變化呈現(xiàn)明顯的更替狀態(tài)，極小曲殼藻(Achnanthes minutissima)和隱頭舟形藻(Navicula cryptocephala)在枯、豐水期都成為優(yōu)勢種，而Cyclotella fottii、具星小環(huán)藻(Cyclotella stelligera)和Cymbella laveis僅在枯水期為優(yōu)勢種，Amphora montana、Achnanthes petersenii、Cymbella minuta、變異直鏈藻(Melosira varians)、短喙舟形藻(Navicula rostellata)和谷皮菱形藻(Nitzschia palea)僅在豐水期為優(yōu)勢種。

2.3 硅藻群落聚類分析

在同等相似水平下，枯水期聚類分析將 19 個樣點劃分為 3 個硅藻類群(圖3)，類群 A1 優(yōu)勢種組成包括C.fottii、C.stelligera、A.minutissima和C.laveis等；類群 A2 優(yōu)勢種A.minutissima和Nitzschia valdecostata的豐度較高；類群 A3 優(yōu)勢種為N.cryptocephala、A.montana和A.minutissima等(表2)。豐水期 19 個樣點被劃分為 2 個硅藻類群(圖3)，類群 B1 優(yōu)勢種組成主要以N.palea為主；類群 B2 優(yōu)勢種為A.minutissima、A.montana和A.petersenii等(表3)。

表1 柳江流域硅藻優(yōu)勢種組成Table 1 Composition of dominant species of diatoms in Liujiang River Basin

2.4 水質(zhì)因子及其主成分分析

枯水期 DO 隨著類群 A1、A2、A3 在空間分布上的變化呈現(xiàn)依次降低的趨勢，其他水化學(xué)指標(biāo)相反。豐水期 DO 隨著類群 B1、B2 的變化呈上升趨勢，其他指標(biāo)則相反(表4)。單因素方差分析結(jié)果表明，水質(zhì)理化指標(biāo)在枯、豐水期存在顯著性差異(P＜ 0.05)，且 DO、CODMn、NH4+-N、Hg、Cd 和Pb 在豐水期的平均值均高于枯水期。主成分分析(PCA)結(jié)果表明，選定與兩主軸相關(guān)系數(shù)大于 0.3的因子為主要水質(zhì)驅(qū)動因子，枯水期為 DO、CODMn、As、Hg、Cd 和 Pb；豐水期為 DO、CODMn、NH4+-N、As 和 Cd。

2.5 硅藻群落與水質(zhì)因子的關(guān)系

DCA 顯示枯、豐水期的排序軸梯度最大值為2.545 和 2.414，均小于 3.0，表明柳江流域硅藻物種對水質(zhì)的響應(yīng)是線性的，故應(yīng)采用冗余分析(RDA)來分析硅藻群落與水質(zhì)因子之間的關(guān)系。

圖3 枯水期(a)和豐水期(b)各樣點硅藻群落的聚類分析Figure 3 The cluster analysis of diatom communities in dry season(a) and wet(b) season

表2 枯水期各類群硅藻優(yōu)勢種Table 2 Dominant species of diatoms in various groups in dry season

表3 豐水期各類群硅藻優(yōu)勢種Table 3 Dominant species of diatoms in various groups in wet season

表4 水質(zhì)理化特征(平均值±標(biāo)準(zhǔn)差)Table 4 Physical and chemical characteristics of water quality(Mean ± SD)

RDA 分析結(jié)果顯示，枯水期軸 1 與軸 2 的特征值為 0.116 和 0.099，豐水期前兩軸的特征值分別為 0.133 和 0.064，結(jié)合物種-環(huán)境關(guān)系的累積百分比，枯、豐水期分別解釋了硅藻群落變化的59.7%和 60.3%，說明枯、豐水期的 RDA 排序圖都能較好地反映硅藻種類分布與水質(zhì)因子之間的相關(guān)關(guān)系。在枯水期，重金屬因子 Pb、Cd、Hg、As 和CODMn位于第一軸的右側(cè)，與軸 1 呈正相關(guān)，DO位于第三象限，與軸 1 和軸 2 均呈負(fù)相關(guān)(圖4)；在豐水期，Cd、As 和 NH4+-N 位于第一軸的右側(cè)，與軸 1 呈正相關(guān)，DO 位于第一軸的左側(cè)，則呈負(fù)相關(guān)，CODMn位于軸 2 的下方，與軸 2 呈負(fù)相關(guān)(圖5)。同時結(jié)合水質(zhì)因子與前 2 軸的相關(guān)性分析結(jié)果(表5)，可以看出枯水期主要水質(zhì)驅(qū)動因子對硅藻群落分布的影響程度依次表現(xiàn)為 Pb＞Cd＞Hg＞ As＞CODMn＞DO，豐水期則表現(xiàn)為 As＞DO＞CODMn＞ NH4+-N＞Cd。

圖4 枯水期硅藻物種與水質(zhì)因子的RDA排序圖Figure 4 RDA ordination diagram of diatom species and water quality factors in dry season

圖5 豐水期硅藻物種與水質(zhì)因子的RDA排序圖Figure 5 RDA ordination diagram of diatom species and water quality factors in wet season

表5 RDA前2軸與水質(zhì)因子之間相關(guān)系數(shù)Table 5 Correlations between water quality factors and the first two axis of RDA

RDA 排序圖表明，硅藻優(yōu)勢種與水質(zhì)驅(qū)動因子的對應(yīng)關(guān)系趨勢在不同水文季節(jié)表現(xiàn)基本一致。其中，硅藻優(yōu)勢種A.minutissima、C.laveis、Aulacoseira ambigua和C.stelligera等都分布在軸2 的左側(cè)，主要與溶解氧呈正相關(guān)，表現(xiàn)出偏好溶解氧高，有機(jī)負(fù)荷程度低，重金屬污染小的水環(huán)境；而硅藻A.montana、A.petersenii、N.cryptocephala、M.varians、C.minuta、N.palea和N.rostellata等則分布在軸 2 的右側(cè)，主要與重金屬因子、CODMn和 NH4+-N 呈正相關(guān)，表明其更喜好有機(jī)營養(yǎng)負(fù)荷高或重金屬含量高的水環(huán)境。

3 討論

3.1 硅藻群落結(jié)構(gòu)特征

本研究中，硅藻群落在枯、豐水期具有明顯的空間分布特征，主要可能是受各采樣點的水質(zhì)因子所影響[17]，河水的水體理化特征往往因人為活動和降雨等自然因素的驅(qū)動而存在不同的時空差異性[8]。在柳江流域中，枯水期有 3 個不同的硅藻類群大致依次分布在該流域干流及支流段的上游、中游以及下游。其中，樣點 S1 和 S5 分別位于林地和景觀區(qū)中，S7、S8 和 S11 樣點周邊多為裸地，說明上游幾個采樣點受人為活動影響較小，水質(zhì)狀況良好，硅藻優(yōu)勢種組成的相似性較高；中游大多數(shù)采樣點(S3、S6、S9、S10 和 S15)主要分散在鄉(xiāng)村附近，人口較為密集，養(yǎng)殖、農(nóng)業(yè)和居民生活等活動明顯，水體主要受生活污水和農(nóng)業(yè)養(yǎng)殖用水排放的影響，使得流入水體的營養(yǎng)鹽等陸源物質(zhì)類型相似，營養(yǎng)水平上升，對硅藻群落影響較為一致；下游采樣點(S16、S18 和 S19)分布在城市附近，工業(yè)生產(chǎn)等人為活動頻繁，伴隨著中上游污染物的積累，導(dǎo)致水體中重金屬及有機(jī)物等污染物的含量增加，進(jìn)而使得硅藻群落組成的相似性也較高。 3 個硅藻類群的空間變化體現(xiàn)了采樣點周圍人為活動方式及強(qiáng)度的不同導(dǎo)致水質(zhì)的污染來源有所差異[18-19]，從而使得水質(zhì)逐漸變差，具有明顯的空間梯度變化。進(jìn)入豐水期，群落結(jié)構(gòu)演替為 2 個不同的硅藻類群，基本依次呈現(xiàn)在該流域的中下游和上游。豐水期雨量較大，地表徑流的增加促使有機(jī)物等陸源性物質(zhì)被沖刷匯入河流，導(dǎo)致各種污染物和離子等濃度的增加[20]，同時充沛的水量增強(qiáng)了河流水系之間的連通性，也促進(jìn)了中下游水體之間的流動，從而縮小了水質(zhì)在其空間分布上的差異，提高了群落結(jié)構(gòu)的相似性。因此，枯、豐水期硅藻群落結(jié)構(gòu)類群劃分的不同體現(xiàn)了地表徑流受雨季影響的變化，說明硅藻群落結(jié)構(gòu)及物種組成與水化學(xué)變量和水文季節(jié)的變化有關(guān)[21-22]，可能是人為活動和降雨等水文過程可間接改變水體中的各項理化指標(biāo)，進(jìn)而造成附生硅藻群落分布的差異性[23-24]。

根據(jù) Van 等[25]的硅藻生態(tài)指示劃分，C.stelligera、C.laveis和A.ambigua為貧污染性種類，主要出現(xiàn)在溶解氧高，有機(jī)污染負(fù)荷和營養(yǎng)水平較低的水體中，而A.montana、A.petersenii、N.cryptocephala、M.varians和G.parvulum為 α —中污染性種類，常生活在中營養(yǎng)狀況的水體中，能夠忍耐較高的有機(jī)污染物，與本文的RDA排序結(jié)果一致。研究發(fā)現(xiàn)，運(yùn)動性較強(qiáng)的N.palea是富營養(yǎng)化水體和重金屬污染水體的常見種[26]，在本文中，其主要分布在豐水期柳江流域的下游，偏好于有機(jī)營養(yǎng)負(fù)荷和重金屬含量高的水體。因此，硅藻群落組成的變化體現(xiàn)了這些種類對水質(zhì)環(huán)境變化的適應(yīng)性。本文硅藻優(yōu)勢種組成的季節(jié)演替伴隨著由更敏感的物種向更耐污染的物種組成轉(zhuǎn)變(表1)，與其他學(xué)者[27-28]研究的結(jié)論大體相同，反映了該流域豐水期水質(zhì)的污染程度高于枯水期，這與理化監(jiān)測結(jié)果一致(表4)。但值得注意的是，本文 RDA 排序分析顯示A.minutissima偏好于污染較小的水環(huán)境，但它在柳江流域各個采樣點的出現(xiàn)頻率卻為 100%，可能是由于該種的分布范圍較廣，其對水體環(huán)境并不敏感，無論在是否受到污染的水體中，它都有可能成為優(yōu)勢種[29]。

3.2 水質(zhì)因子對硅藻群落的影響

RDA 結(jié)果表明，水質(zhì)因子 Cd、As、DO 和CODMn在枯水期和豐水期都對柳江流域附生硅藻群落分布產(chǎn)生影響，但其影響程度不完全一致。由于不同硅藻物種對水體重金屬的敏感程度和耐受能力不同[30]，重金屬則是影響附生硅藻群落組成的重要因素[31]。在本研究中，枯水期附生硅藻分布受水體重金屬的影響程度較大。有研究報道，具星小環(huán)藻(C.stelligera)和C.laveis與重金屬濃度呈負(fù)相關(guān)[32]，與本文研究結(jié)果一致。這一現(xiàn)象說明枯水期硅藻優(yōu)勢種在空間分布上的變化可能是由于其對重金屬含量的變化較為敏感。研究發(fā)現(xiàn)，豐水期優(yōu)勢種谷皮菱形藻(N.palea)與重金屬 As、Cd 濃度之間呈正相關(guān)[33]，且在柳江流域中，Cd、As 是影響硅藻分布的主要重金屬因子，在一定程度上反映了該流域重金屬污染可能主要以 Cd、As 污染為主[12]。

溶解氧是衡量水體初級生產(chǎn)力的重要標(biāo)志，它在柳江流域也是影響附生硅藻的主要水質(zhì)因子之一，這與其他學(xué)者在其他區(qū)域的研究結(jié)果相似。陳向等[34]研究東江惠州河流段的硅藻群落時，認(rèn)為溶解氧對群落結(jié)構(gòu)起著重要作用。在本研究中，水體溶解氧含量在豐水期較高，其對硅藻影響較大。究其原因，一方面可能是豐水期水生植物和藻類的生長迅速，在光合作用下釋放大量的氧氣；另一方面可能是豐水期降雨頻繁，水體持續(xù)流動增加了空氣中氧的溶解效率，良好的氧化環(huán)境有利于有機(jī)物的礦化和營養(yǎng)鹽的再生，從而間接影響硅藻群落結(jié)構(gòu)的特征[35-36]。然而，在柳江流域下游，硅藻物種豐富度較高，但水體中溶解氧含量卻下降，可能是由于人為活動頻繁，大量污水排放和中上游污染物的積累致使水體中營養(yǎng)水平等上升，從而導(dǎo)致藻類及其他浮游生物的繁殖，消耗大量的氧氣，使得水體中溶解氧量下降[37]。

柳江流域硅藻群落在一定程度上也受水體中CODMn的影響，說明其在解釋硅藻群落結(jié)構(gòu)的時空變化中起著重要作用。CODMn是衡量水體中有機(jī)物含量的指標(biāo)，其值越大，有機(jī)物越多，微生物經(jīng)氧化分解產(chǎn)生的可被藻類利用的營養(yǎng)鹽越多，而相關(guān)研究認(rèn)為營養(yǎng)鹽也是影響硅藻生長的重要因子[38]。本文硅藻優(yōu)勢種與水體 CODMn存在一定的響應(yīng)關(guān)系，可能是由于人為活動和地表徑流導(dǎo)致水體中CODMn出現(xiàn)濃度梯度的變化，從而影響硅藻群落的分布。從研究結(jié)果看，本文所涉及的水質(zhì)因子在枯、豐水期分別解釋了 59.7%和 60.3%的硅藻群落變異，尚有部分變化還未能得到相應(yīng)的解釋。除此之外，其他如水溫[35]、透明度[18]及藍(lán)藻、浮游動物的生物作用[39]等也是影響硅藻種群密度和生物量的重要因子。因此，在探討硅藻群落結(jié)構(gòu)變化的機(jī)理時，還需要對其他影響硅藻種群特征的影響因素做進(jìn)一步的基礎(chǔ)生態(tài)學(xué)研究。

4 結(jié)論

(1)柳江流域共鑒定出附生硅藻 137種，隸屬于6目9科34屬，主要優(yōu)勢藻屬以舟形藻屬和曲殼藻屬為主。其中，硅藻種在枯水期和豐水期的變化不大，種類相似度為 63.68%，而各采樣點硅藻豐富度隨季節(jié)變化明顯，表現(xiàn)為豐水期高于枯水期。硅藻優(yōu)勢種組成在不同水期也存在差異，極小曲殼藻和隱頭舟形藻為枯、豐水期共有的優(yōu)勢種，而其余優(yōu)勢種發(fā)生了更替。

(2)聚類分析顯示，柳江流域19個采樣點在枯水期可劃分為 3個硅藻類群，在豐水期則可劃分為 2個硅藻類群，其群落結(jié)構(gòu)類群在空間和時間上劃分的不同，可能主要與采樣點周圍人為活動和地表徑流導(dǎo)致水質(zhì)的空間梯度變化有關(guān)。

(3)主成分分析(PCA)和冗余分析(RDA)結(jié)果表明，Cd、As、DO 和 CODMn在枯、豐水期都對柳江流域附生硅藻群落分布產(chǎn)生影響，但其影響程度發(fā)生變化。其中重金屬因子在枯水期對硅藻影響較大，而溶解氧在豐水期對硅藻影響較大。

附表1 硅藻種類及縮寫Appendix 1 Diatom species and abbreviation