相華，郭偉，竇冰，殷旭旺，王博涵

（1.濟(jì)南市水文局，山東 濟(jì)南 250014；2.大連海洋大學(xué)水產(chǎn)與生命學(xué)院，遼寧省水生生物學(xué)重點(diǎn)實(shí)驗(yàn)室，遼寧 大連 116023）

在水生生物群落的研究中，浮游植物群落一直是研究的主要對(duì)象之一[1,2]。浮游植物群落的物種數(shù)量、物種豐度、優(yōu)勢(shì)度和多樣性指數(shù)等結(jié)構(gòu)特征經(jīng)常作為評(píng)價(jià)水生態(tài)系統(tǒng)的指標(biāo)之一[3]，這些年來(lái)，早已經(jīng)被廣泛的應(yīng)用于不用國(guó)家、不同城市、不同環(huán)境的水生態(tài)系統(tǒng)狀況的評(píng)價(jià)之中[4-6]。

浮游植物功能群的研究方法考慮到了浮游植物群落差異性的不同，根據(jù)不同種類的浮游植物群落對(duì)不同環(huán)境因子的限制、個(gè)體大小、空間分布的不同及對(duì)不同環(huán)境的適應(yīng)性的不同等多個(gè)因素進(jìn)行劃分[7,8]。浮游植物功能群的劃分是進(jìn)一步分析了浮游植物的群落結(jié)構(gòu)特點(diǎn)與其適宜生活的水環(huán)境的特征，包括：水層、水深、營(yíng)養(yǎng)化等多個(gè)因素。功能群的劃分跟以往研究的群落形態(tài)的劃分方法不同。功能群的劃分可以更加準(zhǔn)確、明了地顯現(xiàn)出不用浮游植物類群與不同水環(huán)境因子的關(guān)聯(lián)性。目前為止，浮游植物功能群對(duì)水生態(tài)學(xué)的研究及環(huán)境評(píng)價(jià)已經(jīng)有了十分廣泛的應(yīng)用和研究[9-13]。

小清河流域是貫穿山東省多個(gè)城市的主要流域之一，而其中的黃河流域是山東的渤海水系河流，他起源于濟(jì)南的泉水群。小清河?xùn)|流經(jīng)的地域包括濟(jì)南槐蔭區(qū)、天橋區(qū)、歷城區(qū)、章丘區(qū)，濱州市的鄒平、高青、桓臺(tái)、博興等縣區(qū)，濰坊市的青州市至壽光市的羊角溝匯入渤海，全長(zhǎng)共233 km，是一條防洪除澇、灌溉的綜合利用河道。

近年來(lái)，小清河流域水的生態(tài)環(huán)境受到了很大的影響及破壞。本次研究，以小清河流域流經(jīng)濟(jì)南的部分水體作為主要研究對(duì)象，分析了2018 年和2019 年小清河流域的浮游植物的空間特征，以期找出影響功能群的主要的水環(huán)境因子，并為小清河流域的環(huán)境保護(hù)和治理提供基本數(shù)據(jù)支持。

1 材料和方法

1.1 研究區(qū)域劃定及采樣點(diǎn)的設(shè)置

小清河流域的地理、地貌特征比較復(fù)雜，為了更好研究小清河流域中浮游植物群落結(jié)構(gòu)的空間差異，根據(jù)兩個(gè)季節(jié)中共同采集的點(diǎn)位、濟(jì)南市水環(huán)境監(jiān)測(cè)的點(diǎn)位及各縣、區(qū)的主要水資源等多個(gè)因素，共設(shè)置了9 個(gè)具有代表性的采樣點(diǎn)（圖1）。本次采樣點(diǎn)位多流經(jīng)市區(qū)，河流的支流較短，水土流失比較嚴(yán)重，水中泥沙含量較大，底質(zhì)主要為基石。

圖1 小清河流域采樣點(diǎn)分布Fig.1 The sampling sites distribution in Xiaoqing River Basin

1.2 浮游植物樣品采集、鑒定及水體環(huán)境因子測(cè)定

參照其他人員的采樣方式，在每一個(gè)采樣點(diǎn)采集2 L 的水樣，用魯哥氏液固定。水樣送回實(shí)驗(yàn)室后，將樣品靜置48 h 后，用虹吸法將水樣濃縮至100 mL，在顯微鏡下鑒定物種，記錄和整理數(shù)據(jù)[14-17]。

按照殷旭旺、王博涵等[17,18]的采樣方法測(cè)定透明度、渾濁度、水溫等水環(huán)境因子。采集的2 L 水樣，用低溫保存，在48 h 內(nèi)送回實(shí)驗(yàn)室。按照標(biāo)準(zhǔn)方法[17,18]測(cè)定總氮、硝酸態(tài)氮、氨氮含量等多個(gè)水環(huán)境因子。

1.3 浮游植物功能群的劃分

結(jié)合Reynolds、閔文武和王博涵等[9-12,19,20]對(duì)浮游植物功能群劃分的標(biāo)準(zhǔn)方法，對(duì)小清河流域浮游植物群落進(jìn)行功能群的劃分。根據(jù)浮游植物群落的特點(diǎn)、自身特性、生活環(huán)境等不同因素，共劃分出16類浮游植物功能群（表1）。

表1 清河流域浮游植物功能群組成Tab.1 Functional group composition of phytoplankton in Xiaoqing River Basin

1.4 數(shù)據(jù)統(tǒng)計(jì)及分析方法

計(jì)算浮游植物多樣性指數(shù)及優(yōu)勢(shì)種。

香農(nóng)威納指數(shù)（H'）：H'=－Σ（ni/N）×log2（ni/N）

均勻度指數(shù)（J）：J=H'/lns

浮游植物群落優(yōu)勢(shì)種的確定[13,18-20]：

Y=（ni/N）×fi。

群落及環(huán)境因子數(shù)據(jù)的相關(guān)性分析采用SPSS 19.0 進(jìn)行分析，群落特征值分析及功能群特征值分析采用Excel 2016 分析、Biodiversity pro 2.0 計(jì)算浮游植物群落的香農(nóng)威納指數(shù)，Canoco 4.5 進(jìn)行環(huán)境因子的主成成分分析（PCA）和選出的水環(huán)境因子和浮游植物功能群的典范對(duì)應(yīng)分析（CCA）的制圖，在Arc Map 10.2 上制作小清河流域圖[18-20]。

2 結(jié)果與分析

2.1 浮游植物群落結(jié)構(gòu)及其功能群的劃分

2.1.1 浮游植物物種組成及其群落結(jié)構(gòu)

本次調(diào)查中，小清河流域在2018 年和2019 年中分別鑒定出32 和54 種浮游植物。2018 年小清河流域共鑒定出硅藻門13 種，綠藻門8 種，其他物種門類相對(duì)較少，分別是：藍(lán)藻門5 種、裸藻門3 種、隱藻門2 種和甲藻門1 種。2019 年小清河流域共鑒定出硅藻門24 種，綠藻門19 種，其他物種門類相對(duì)較少，分別是：裸、隱藻各3 種、藍(lán)藻門2 種、甲、黃、金藻門各1 種。2018 年小清河流域藍(lán)藻門密度占比較大，占總密度的89.30%，其他門類密度較少，綠藻門占總密度的7.22%，硅藻門2.69%，裸藻門0.35%，隱藻門0.31%，甲藻門0.12%。2019 年小清河流域浮游植物門類密度與2018 年有所不同，硅藻門密度占比較大，占總密度的60.22%，其他門類密度較少，綠藻門占21.85%，隱藻門16.41%，裸藻門1.00%，裸藻門0.41%，甲藻門、黃藻門和金藻門占比最少（圖2）。

圖2 不同年份小清河流域浮游植物群落的種類組成及其所占比例Fig.2 Species number and percentage of phytoplankton in Xiaoqing River Basin in different years

2018 年春季浮游植物的優(yōu)勢(shì)種有3 種，以綠藻門和藍(lán)藻門為主，分別為四尾柵藻（Scenedesmus quadricauda）（Y=0.026）、優(yōu)美平列藻（Merismopedia elegans）（Y=0.206）和銀灰平裂藻（M.glauca）（Y=0.587）；2019 年的優(yōu)勢(shì)種共有4 種，以硅藻門和綠藻門為主，分別為尖針桿藻（Synedra acus）（Y=0.076）、具星小環(huán)藻（Cyclotella stelligera Cl.et Grun.）（Y=0.047）、弧形短縫藻（Eunotia arcus）（Y=0.028）和小球藻（Chlorella vulgaris Beijerinck）（Y=0.047）。兩年中，浮游植物群落優(yōu)勢(shì)種無(wú)顯著相關(guān)性（P>0.05）。

2018 年，全流域各采樣點(diǎn)浮游植物物種數(shù)平均值為7.11 種；全流域各采樣點(diǎn)浮游植物物種密度平均值為7.06×107cell/L；香農(nóng)威納指數(shù)的變化范圍是2.32～3.32，平均值為2.78；均勻度指數(shù)變化范圍為0.63～0.90，平均值為0.75。2019 年，全流域各采樣點(diǎn)浮游植物物種數(shù)平均值為15.33 種；全流域浮游植物物種密度平均值為1.42×107cell/L；香農(nóng)威納指數(shù)的變化范圍是0.86～3.18，平均值為2.27；均勻度指數(shù)變化范圍為0.27～0.93，平均值為0.59（表3）。2018 和2019 年浮游植物物種數(shù)、物種密度、均勻度指數(shù)、香農(nóng)威納指數(shù)差異均不顯著（P>0.05）。

表3 不同年份小清河流域浮游植物群落的群落特征數(shù)值Tab.3 Community characteristic value of phytoplankton in Xiaoqing River Basin in different years

2.1.2 浮游植物功能群的劃分及研究

本次研究中，共劃分為16 類功能群（表1），其中，2018 年共劃分出9 類功能群，分別為：P、MP、Lo、D、C、W1、Y、X1 和J。2019 年共劃分出16 類功能群，分別為：B、P、MP、Lo、D、A、C、W1、Y、X1、X2、J、T、F、N 和E。P、MP、Lo、D、C、W1、Y、X1、J，這9 類是兩個(gè)年份中共同出現(xiàn)的功能群，2019 年特有的功能群有：B、X2、T、F、N、E 和A。

就功能群各類的數(shù)目而言，2018 和2019 年均以MP 和J 數(shù)目較多，但功能群密度卻不相同。單看功能群物種數(shù)：2018 年中，MP 數(shù)目最高，為9 種，J為6 種，其他功能群數(shù)量相對(duì)較少，Lo 為4 種，P、W1 和Y 為3 種，X1 為2 種，C 和D 數(shù)目最少，僅為1 種。2019 年J 數(shù)目最高，為13 種，其次為MP，為12 種，D 為8 種，Y 為4 種，W1 和P 為3 種，X1 為2 種，其他為1 種。單看功能群物種密度，兩年差異性較大，2018 年Lo 密度占有絕大優(yōu)勢(shì)，占總密度的88.72%，其次為J 類，占總密度的6.01%，其他功能群密度均相對(duì)較少；2019 年中，MP 和D 類密度相對(duì)較多，分別占總密度的24.57%和23.31%，其次為Y和X1，占總密度的16.46%和15.67%，其他功能群密度相對(duì)較少（圖3）。浮游植物功能群的數(shù)目和密度在2018 年和2019 年均無(wú)顯著相關(guān)性（P>0.05）。

圖3 不同年份小清河流域浮游植物功能群各類數(shù)目和密度Fig.3 The number and density of phytoplankton functional groups in Xiaoqing River Basin between 2018 and 2019

2.2 浮游植物功能群與水環(huán)境因子的相關(guān)性

2.2.1 小清河流域內(nèi)水環(huán)境因子特征

2018 年和2019 年小清河流域測(cè)定出的水環(huán)境因子中，氯化物表現(xiàn)出較顯著的時(shí)間差異（P<0.01），其他16 個(gè)水環(huán)境因子均無(wú)顯著差異性（P>0.05）。這兩年春季整體水溫均在25℃左右，水體溫度較高，pH 平均分別為8.28 和8.47，水體呈弱堿性，2018年中，水中氯化物、氨氮、亞硝酸鹽氮、高錳酸鹽指數(shù)和氟化物的含量高于2019 年，其他水環(huán)境因子水中含量均低于2019 年（表4）。

表4 2018 和2019 年春季小清河流域水環(huán)境因子比較（平均值±標(biāo)準(zhǔn)差）Tab.4 Comparison of environment factors（Mean± SD）in Xiaoqing River Basin in spring of 2018 and 2019

2.2.2 浮游植物功能群與水環(huán)境因子的相關(guān)性

水環(huán)境因子主成分分析結(jié)果表明（圖4）：2018年，通過(guò)數(shù)據(jù)分析，找出兩主軸變化分?jǐn)?shù)較大的5種水環(huán)境因子，分別為pH、氯化物、溶解氧、總氮和硝酸鹽氮；2019 年，通過(guò)數(shù)據(jù)分析，找出兩主軸變化分?jǐn)?shù)較大的5 種水環(huán)境因子分別為電導(dǎo)率、氯化物、硫酸鹽、總氮和硝酸鹽氮。氯化物、總氮和硝酸鹽氮是2018 年和2019 年選出相同的水環(huán)境因子。

圖4 小清河流域不同年份水環(huán)境因子主成分分析（PCA）Fig.4 Principal component analysis of water environmental factors in Xiaoqing River Basin

2018 年和2019 年浮游植物功能群和選出的水體環(huán)境因子的CCA 分析結(jié)果（圖5）顯示：2018 年，溶解氧在軸1 上對(duì)浮游植物功能群的影響最大，對(duì)軸1 呈負(fù)相關(guān)性，并且主要影響LO功能群；pH 主要對(duì)軸1 浮游植物功能群的影響最大，呈負(fù)相關(guān)性?？偟拖跛猁}氮對(duì)軸1 的影響較大，呈正相關(guān)性，主要影響D、C、J 和P 類功能群。2019 年，總氮對(duì)浮游植物功能群的影響最大，主要影響B(tài)、D、C、J 和MP 類功能群。氯化物和電導(dǎo)率對(duì)軸1 呈負(fù)相關(guān)性，主要影響X1 和B 類功能群。

圖5 小清河流域不同年份浮游植物與環(huán)境因子CCA 分析Fig.5 Canonical correspondence analysis of species-environmental relationship in Xiaoqing River Basin between 2018 and 2019

2.3 浮游植物群落水生態(tài)健康評(píng)價(jià)

近些年來(lái)都用浮游植物的香農(nóng)威納指數(shù)來(lái)評(píng)價(jià)水生態(tài)環(huán)境的健康，這也是最常用的評(píng)價(jià)水生態(tài)系統(tǒng)的方法之一。此類方法與其他水生生物群落評(píng)價(jià)方法一致（表5）。從浮游植物評(píng)價(jià)結(jié)果看，本次研究中，2018 年和2019 年采集到的浮游植物群落的香農(nóng)威納指數(shù)平均值分別為2.78 和2.27，均勻度指數(shù)平均值分別為0.75 和0.59。單從浮游植物群落多樣性指數(shù)的角度來(lái)看，2018 年和2019 年小清河流域水體質(zhì)量一般[1,2,16,17]。

表5 基于浮游植物多樣性指數(shù)的小清河流域水生態(tài)健康標(biāo)準(zhǔn)Tab.5 Ecological health of water in Xiaoqing River Basin according to diversity index of phytoplankton

3 討論

3.1 小清河流域浮游植物功能群及其群落結(jié)構(gòu)特征

調(diào)查結(jié)果表明，小清河流域浮游植物群落結(jié)構(gòu)有一定的空間和時(shí)間差異，2018 年物種數(shù)為32 種，物種數(shù)以硅藻門為主，密度以藍(lán)藻門為主，優(yōu)勢(shì)種為優(yōu)美平列藻、銀灰平裂藻和四尾柵藻，藍(lán)藻密度占絕對(duì)優(yōu)勢(shì)。這主要是由于本次鑒定出的平裂藻細(xì)胞數(shù)量較多。濟(jì)南春季水溫已達(dá)到25℃左右，十分適合平裂藻的生長(zhǎng)和繁殖。據(jù)監(jiān)測(cè)，2018 年水中氮、磷等多項(xiàng)水環(huán)境因子含量同樣相對(duì)較高，特別適合藍(lán)藻門的物種快速繁殖，其密度占主導(dǎo)優(yōu)勢(shì)，其他物種密度相對(duì)較少。相比2018 年，2019 年物種數(shù)要高很多，共鑒定出54 種，物種數(shù)和密度均以硅藻門為主，優(yōu)勢(shì)種為小球藻、尖針桿藻、具星小環(huán)藻和弧形短縫藻，優(yōu)勢(shì)種中硅藻門也占絕對(duì)優(yōu)勢(shì)。2019 年小清河流域的水量比2018 年少，而硅藻門類物種大多喜歡生活在活水體的淺水環(huán)境中，這種環(huán)境適合硅藻門的生長(zhǎng)和繁殖[21]。許多研究已經(jīng)證明了浮游植物群落的多樣性指數(shù)能在某種程度上評(píng)價(jià)水生態(tài)系統(tǒng)的健康情況[22,23]。本次調(diào)查中，兩年中小清河流域浮游植物香農(nóng)威納指數(shù)平均值均在2～3 之間，均勻度指數(shù)平均值0.5～0.8 之間。兩年中這兩個(gè)指數(shù)平均值均較低，根據(jù)評(píng)價(jià)水質(zhì)標(biāo)準(zhǔn)來(lái)看[24]，小清河流域水體質(zhì)量一般，水體呈中度污染。

2018 和2019 年小清河流域共劃分出16 類浮游植物功能群，2018 年MP 和J 類物種數(shù)最高，Lo類密度最高，2019 年J 和MP 類物種數(shù)較高，MP 和D 類密度較高。鑒定出的MP 類的物種主要在淺水的環(huán)境中，這類功能群喜歡生活在不穩(wěn)定的水體中。小清河流域的水流量歷年來(lái)都比較少，這一特點(diǎn)則剛好適合這類功能群藻種的生長(zhǎng)、繁殖和生活，很多硅藻門的物種也都屬于MP 這一類功能群。硅藻門在這兩年的物種數(shù)均比較高，這也是MP 類的物種在這兩年中物種數(shù)均較高的原因。J 類功能群適宜生活在富營(yíng)養(yǎng)的水體中，這兩年中J 類的物種數(shù)都很高，測(cè)定的水環(huán)境因子中氮、磷、高錳酸鉀指數(shù)等多個(gè)水環(huán)境因子含量均過(guò)高，水體富營(yíng)養(yǎng)化比較嚴(yán)重，特別適合J 類的生長(zhǎng)。J 類主要以綠藻門為主，水中綠藻門的常見(jiàn)物種也比較多，故兩年中J類的物種數(shù)都比較高。本次野外調(diào)查為春季，測(cè)定水體水溫已達(dá)到25℃，陽(yáng)光充足，這樣良好的環(huán)境會(huì)導(dǎo)致喜歡陽(yáng)光的D 類功能群大量出現(xiàn)。尖針桿藻、針形菱形藻等藻類屬于D 類功能群，尖針桿藻又是2019 年樣品中的優(yōu)勢(shì)種，這也是2019 年D 類功能群密度占比高于2018 年的原因之一。2018 年中Lo 類功能群密度占絕對(duì)優(yōu)勢(shì)，Lo 對(duì)各類水體的適應(yīng)性十分廣泛，各類水體均能很好生活。Lo 代表的物種為色球藻和平裂藻，2018 年優(yōu)美平列藻、銀灰平裂藻為優(yōu)勢(shì)種，其細(xì)胞數(shù)量遠(yuǎn)遠(yuǎn)高于其他物種，并在2018 年中占絕對(duì)優(yōu)勢(shì)。鑒定出的B、C、E、F、J、P、N、X1、X2 和W1 等功能群都喜歡生活在富營(yíng)養(yǎng)化的水體中，這一點(diǎn)也說(shuō)明了小清河流域部分水體呈現(xiàn)出富營(yíng)養(yǎng)化的狀態(tài)[9-12,19,20]。

3.2 浮游植物功能群與水環(huán)境因子的關(guān)系

綜合考慮后，運(yùn)用浮游植物功能群與水環(huán)境因子進(jìn)行典范對(duì)應(yīng)分析更具有說(shuō)服力。浮游植物功能群與環(huán)境因子CCA 結(jié)果表明，2018 年，溶解氧含量是影響小清河流域浮游植物功能群的主要環(huán)境因子。溶解氧含量與浮游植物功能群呈負(fù)相關(guān)性，并顯著影響LO。pH 與浮游植物功能群呈負(fù)相關(guān)性。溶解氧一直以來(lái)都是影響浮游植物群落生長(zhǎng)繁殖的必不可少的重要指標(biāo)。溶解氧含量直接影響浮游植物群落的自身狀況，水中的含量如果較低，很可能導(dǎo)致浮游植物群落大量死亡[25,26]。小清河流域水體透明度較低，渾濁度較高，水體較為渾濁，嚴(yán)重影響了水中的光合作用[27]。總氮和硝酸鹽氮對(duì)軸1 的影響較大，主要影響D、C、J 和P 類功能群。氮含量一直都是衡量水生態(tài)環(huán)境受到營(yíng)養(yǎng)物污染程度的十分重要指標(biāo)[19,20,28]，水中氮含量較高，C、D、J 和P 等功能類群大量迅速繁殖，群落密度增加。與2018 年相同，2019 年，氮含量對(duì)浮游植物功能群的影響最大，主要影響B(tài)、D、C、J 和MP 類功能群。硝酸態(tài)氮含量對(duì)D、Lo、T 和A 功能群也有一定的影響。水中氯化物的含量在一定程度上也會(huì)影響浮游植物群落的生活，含量較高也會(huì)導(dǎo)致浮游植物群落的死亡。電導(dǎo)率主要影響水中陰陽(yáng)離子和懸浮物的含量，對(duì)浮游植物群落有一定影響。2019 年中，氯化物和電導(dǎo)率主要影響X1 和B 類功能群。小清河流域受到的人為影響，農(nóng)業(yè)和工業(yè)開(kāi)采的影響都較為嚴(yán)重，對(duì)整個(gè)流域的水生態(tài)系統(tǒng)影響較為嚴(yán)重。氮及其他元素含量的變化，對(duì)浮游植物的生長(zhǎng)繁殖有很直接的影響，2018 年和2019 年中，水中檢測(cè)出的氮含量相對(duì)較高，且對(duì)很多浮游植物類群有顯著影響[29-32]。整體看來(lái)，不同季節(jié)中，影響浮游植物群落功能群的驅(qū)動(dòng)因子不同，影響浮游植物功能群的驅(qū)動(dòng)因子也各不相同。

根據(jù)本次調(diào)查結(jié)果，建議加強(qiáng)防治小清河流域點(diǎn)源和面源的污染。在農(nóng)業(yè)生產(chǎn)中，建立合理的氮、磷施肥體系。加強(qiáng)山區(qū)水土保持的管理，減少水土流失。在工業(yè)生產(chǎn)中，加強(qiáng)工廠污染的治理。針對(duì)化工等企業(yè)廢水、廢氣的排放導(dǎo)致的水體自凈能力減弱，應(yīng)加強(qiáng)的環(huán)境保護(hù)意識(shí)，同時(shí)應(yīng)增加對(duì)污水治理的投入性建設(shè)，改善污水處理設(shè)備，實(shí)現(xiàn)排污達(dá)標(biāo)。逐步減輕廢水、污水、廢氣等排放對(duì)小清河流域的危害。