李佳靜，劉威，鄧培雁，

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

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

0 前言

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

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

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

1 材料和方法

1.1 采樣時(shí)間與位置

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

1.2 采樣與分析方法

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

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

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

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

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

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

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

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

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

2 結(jié)果與分析

2.1 硅藻群落種類組成

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

2.2 硅藻優(yōu)勢(shì)種

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

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

2.3 硅藻群落聚類分析

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

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

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

枯水期 DO 隨著類群 A1、A2、A3 在空間分布上的變化呈現(xiàn)依次降低的趨勢(shì)，其他水化學(xué)指標(biāo)相反。豐水期 DO 隨著類群 B1、B2 的變化呈上升趨勢(shì)，其他指標(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ū)動(dòng)因子，枯水期為 DO、CODMn、As、Hg、Cd 和 Pb；豐水期為 DO、CODMn、NH4+-N、As 和 Cd。

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

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

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

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

表3 豐水期各類群硅藻優(yōu)勢(shì)種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%，說(shuō)明枯、豐水期的 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)。同時(shí)結(jié)合水質(zhì)因子與前 2 軸的相關(guān)性分析結(jié)果(表5)，可以看出枯水期主要水質(zhì)驅(qū)動(dòng)因子對(duì)硅藻群落分布的影響程度依次表現(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)勢(shì)種與水質(zhì)驅(qū)動(dòng)因子的對(duì)應(yīng)關(guān)系趨勢(shì)在不同水文季節(jié)表現(xiàn)基本一致。其中，硅藻優(yōu)勢(shì)種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)養(yǎng)負(fù)荷高或重金屬含量高的水環(huán)境。

3 討論

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

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

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

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

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

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

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

4 結(jié)論

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

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

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

附表1 硅藻種類及縮寫(xiě)Appendix 1 Diatom species and abbreviation