湿地具有强大的生态净化作用,被誉为地球之肾。近年来,随着城市人口的集聚和建设用地的扩张,许多地区的自然湿地不断减少,为保护湿地生态系统和促进城市可持续发展,开展城市湿地公园建设,并将其纳入城市公园系统中,已经成为城市湿地保护管理的重要途径[1]。目前,虽然对湿地公园的动物物种多样性进行了较多研究[2],对城市湿地公园生态系统也进行了一些报道[3-4],但关于湿地公园浮游植物多样性的报道较为鲜见。
浮游植物在水体中扮演着多种重要角色,既是水域生态系统的初级生产者,又是鱼类等经济动物的重要食物来源,同时浮游植物的光合作用还是水体中溶解氧的主要生产方式之一,对水体中的生物环境具有重要影响[5]。此外,浮游植物也经常被作为评估水生生态系统健康状况及化学污染物对水生生态系统影响的重要指标之一。
目前,常用于浮游植物分类的方法主要为传统的种属分类法和功能群分类法。传统分类法是基于亲缘关系的系统分类法,而功能群分类法是由Reynolds等[6]基于适应性和生境特征所提出,后经Padisk等[7]补充完善,现已划分出39个功能群。两种方法在水环境质量评价中均被广泛应用,但不同生境类型对评价方法的适用性有明显的差异。赵耿楠[8]研究发现,传统分类方式反映的生境特征与实际状况存在偏差,而功能群法能够较好地反映实际状况;而郭芳等[9]对南亚热带地区8座水库的浮游植物研究则表明,传统分类方式要优于功能群分类法。因此,对不同水环境特征进行分析时要考虑方法的适用性和合理性。东莞生态园原是东江水系的一部分,随着东江的冲刷淤积,逐渐与东江隔离,成为一块低洼地。20世纪80年代以来,随着工业化和城市化的推进,该区域逐渐沦为污水汇集区。2006年,东莞市对该区域统筹整合建设了东莞生态园产业园区,实施一系列湿地修复工程,经过七年多的生态建设和修复,该园区在2013年12月获批成为国家城市湿地公园。东莞生态园是东莞国家城市湿地公园的重要组成部分,兼具蓄洪滞洪、生物保育及休闲观光等多重功能。但东莞生态园湿地2012—2013年水质整体为劣Ⅴ类[10],近期对该生态园的鱼类调查显示,耐污染外来种罗非鱼占优势[4],说明该湿地生态环境现状依然不容乐观。目前,相关研究中对东莞生态园湿地的报道较为鲜见,故有必要对其包含浮游植物在内的水生生物进行深入研究。本研究中,通过传统分类法和功能群分类法对东莞生态园浮游植物群落结构进行研究,对不同水质评价方法进行比较,以期为水生态系统保护、修复和管理提供科学参考。
东莞生态园位于广东省东莞市中心城区的东部,2011年建成并投入使用,由自然水体、人工河道、湖泊及湿地组成,总面积为2.83 km2,其中水域面积为1.78 km2。自然河道主要有东江干流、东江支流寒溪河,人工河道主要有东引河、南畲朗排渠及其支渠、大圳埔排渠及其支渠和新建排渠,人工湖包括拓塘而形成的月湖及中央岛群。月湖是东莞生态园核心景区之一,也是东莞国家城市湿地公园面向公众开放的重要水上游览区,西临中央岛群,南临广深铁路,北面为南畲朗河;中央岛群内水系发达,池塘湿地密布,水网纵横,呈现出塘、渠和河3种景观类型。该湿地地处亚热带季风气候区,年平均气温约为23.6 ℃,平均水深为4 m,气温温和,雨量充沛,光照充足。近年来,随着周边污水的输入和水体自净能力的降低,湖泊水质在成湖后不断恶化。
1.2.1 样品采集与处理 分别于2019年4、7、9、12月在东莞生态园进行浮游植物调查,因受特殊情况影响,未能在2020年更具冬季特征的1月采样,故采用孙莹蓓等[11]提出的方法,将4、7、9、12月分别归为春、夏、秋、冬四季。依据进出湿地的排渠与河流、湿地形状、水流等特征,在中央岛群水域和月湖共设置7个采样位点,其中S1、S2、S3和S4位于中央岛群,S5、S6和S7位于月湖(图1)。
图1 东莞生态园采样点分布
Fig.1 Distribution of sampling sites in Dongguan Ecological Park
按以上月份在每个采样位点进行浮游植物样品采集。采集表层和底层混合水样1 L,加入15 mL鲁哥氏液固定并带回实验室,利用光学显微镜对浮游植物物种进行鉴定,并对藻细胞计数。浮游植物物种鉴定依照《中国淡水藻类:系统、分类及生态》[12],并计算其丰度及生物量[13]。
采用便携式多参数水质检测分析仪(YSI)测定水温(WT)、溶解氧(DO)、电导率(SPC)和总溶解性固体(TDS)水质参数;采用便携式超声波探测仪(型号XHX-SM-5)测量水深(depth),采用赛事盘(Secchi disk)测定水体透明度(SD),采用便携式pH-902计测定水体pH值。采集表层和底层混合水样1 L,根据《水和废水监测分析方法》分析水体总氮(TN)、氨氮(NH+4-N)、亚硝酸盐氮硝酸盐氮总磷(TP)、磷酸盐高锰酸盐指数(CODMn)和叶绿素 a(Chl-a)含量。
1.2.2 功能群划分 将东莞生态园7个站点所鉴定的浮游植物,参照Reynolds等[6]和Padisk等[7]提出的功能群分类法进行划分。为进一步了解浮游植物功能群与环境因子间的关系,将功能群总生物量占比>5%的功能群定义为优势功能群。1.2.3 多样性指数及优势度指数 采用Shannon-Wiener多样性指数(H′)、Pielou均匀度指数(J)、Margalef丰富度指数(D)及优势度指数(Y)计算浮游植物群落特征[14],计算公式为
H′=-∑(ni/N)ln(ni/N),
(1)
J=H′/lnS,
(2)
D=(S-1)/lnN,
(3)
Y=ni/N×fi。
(4)
式中:ni为第i种藻类的个数;N为所有藻类总个数;S为样品中藻类种类数;fi为第i种藻类在各站点出现的频率。将Y≥0.02的藻类定为优势种。
1.2.4 水质评价
1)采用综合营养状态指数(TLI)对东莞生态园水体营养状态进行评价。以Chl-a(μg/L)、TP(mg/L)、TN(mg/L)、 SD(m)和CODMn(mg/L)为参数,各参数的营养状态指数计算公式为
ITL(Chl-a)=10[2.5+1.086ln(Chl-a)],
(5)
ITL(TP)=10[9.436+1.624ln(TP)],
(6)
ITL(TN)=10[5.453+1.694ln(TN)],
(7)
ITL(SD)=10[5.118-1.94ln(SD)],
(8)
ITL(CODMn)=10[0.109+2.661ln(CODMn)],
(9)
ITL(∑)=∑Wj×ITL。
(10)
式中:ITL(∑)为综合营养状态指数(TLI);ITL和Wj分别为第j种参数的营养状态指数及其相关权重。营养状态分级标准见表1[14]。
表1 水质评价标准
Tab.1 Water quality evaluation standard
注:VI—指数值;LWQ—水质等级,下同。
Note:VI—value of index;LWQ—level of water quality,et sequentia.
生态状态指数 Q多样性指数diversity index综合营养状态指数trophic level index丰度/(106 ind.·L-1)abundance生物量/(mg·L-1)biomass叶绿素a/(μg·L-1)chlorophyll-aVILWQH′JDLWQVILWQVILWQVILWQVILWQ4^5极好>30.8^1.0>3无污染<30贫营养<10贫营养<3贫营养≤1贫营养3^4好2^30.5^0.8>3轻污染30^50中营养10^40中营养3^5中营养1^10中营养2^3中等1^20.3^0.51^3中污染50^60轻度富营养41^80中度富营养5^7中度富营养10^26轻度富营养1^2耐受0^10^0.30^1重污染60^70中度富营养81^100富营养7^10富营养26^160中度富营养0^1差>70重度富营养>100重度富营养>10重度富营养>160重度富营养
2)采用生态状况指数(Q)对东莞生态园水质进行评价[15],其计算公式为
(11)
式中:n为浮游植物功能群的数量;N为浮游植物总生物量;ni为第i个功能群的生物量;Fi因子为第i个功能群的赋值。Q指数分级标准见表1[15]。
3)采用浮游植物多样性指数对东莞生态园水质进行评价,分级标准见表1[14]。
4)采用单因子对东莞生态园水质进行评价,根据《地表水环境质量标准(GB 3838—2002)》,最终水质评价类别为最差的水质指标所能达到的类别。
5)采用浮游植物现存量[16]和叶绿素a含量[17]对东莞生态园水质进行评价,分级标准见表1。
1.2.5 聚类分析及冗余分析 利用Primer 5.0软件,采用浮游植物优势种丰度数据构建物种矩阵,进行聚类分析(CA),通过非度量尺度分析(NMDS)结果判断分类可信度。采用Canoco 4.5软件进行去趋势对应分析(DCA),得到优势种和优势功能群4个排序轴中梯度长度最大值分别为1.5和1.55,因此,选择冗余分析(RDA),并采用蒙特卡洛置换检验(P<0.05,n=999)对影响浮游植物优势种的环境变量进行识别和筛选,并用Canodraw 4.5软件作图。
采用ArcGIS 10.0软件绘制采样点图,利用Origin 2022软件绘制浮游植物现存量和功能群时空分布图。利用SPSS 26.0软件进行方差分析,并用Duncan法进行组间多重比较。
各站点CODMn含量较高,水质有机污染较为严重,单因子评价结果显示,各个站点水质均为劣Ⅴ类;其次,各站点总氮含量也较高,单因子评价结果显示,S2样点水质为Ⅴ类,其余6个样点水质均为劣Ⅴ类(表2)。说明东莞生态园整体水质为劣Ⅴ类,有机污染严重,主要污染物质为氮。
表2 2019年东莞生态园各站点主要理化因子及自然环境概况
Tab.2 Main physical and chemical factors and overview of natural environment of each station in Dongguan Ecological Park in 2019
湖区lake area站点site位置situation水深/mdepth水温/℃WT透明度/mSDpH值pH value溶解氧/(mg·L-1)DO总氮/(mg·L-1)TN总磷/(mg·L-1)TP高锰酸盐/(mg·L-1)CODMn氨氮/(mg·L-1)NH+ 4-N磷酸盐/(mg·L-1)PO-3-P中央岛群centralislandsS1出水口2.20±0.0727.88±3.730.32±0.068.16±0.187.86±0.942.48±0.570.04±0.0117.26±2.120.27±0.040.03±0.01S2出水口2.28±0.1727.78±3.700.31±0.088.14±0.188.25±0.931.64±0.290.06±0.0216.57±2.100.28±0.080.04±0.01S3敞水区2.48±0.1827.65±3.660.37±0.088.18±0.188.83±0.872.14±0.440.06±0.0116.66±2.310.28±0.060.03±0.01S4敞水区2.65±0.5327.80±3.550.34±0.078.15±0.188.95±1.132.06±0.300.05±0.0116.05±1.240.25±0.070.03±0.01S5敞水区4.85±0.0527.20±3.130.34±0.108.13±0.228.02±1.363.15±0.790.04±0.0117.79±2.600.29±0.050.03±0.01月湖moonlakeS6敞水区3.73±0.0927.85±3.220.39±0.048.12±0.218.23±1.272.18±0.290.05±0.0116.42±1.960.29±0.060.03±0.01S7入水口4.80±0.1427.50±3.160.36±0.098.17±0.258.64±1.362.76±0.490.05±0.0116.02±1.590.33±0.070.03±0.01
2.2.1 物种组成及优势种 本次调查共鉴定出浮游植物7门137种(属),其中,绿藻门(Chlorophyta)物种数最多,共70种,占总种类数的51.09%,蓝藻门(Cyanophyta)和硅藻门(Bacillariophyta)次之,分别为28种(20.44%)和24种(17.52%),其他5个门的物种数较少,共占10.95%。
根据Y≥0.02判定优势种,东莞生态园各季节中共有20种优势种,隶属于蓝藻门(17种)和绿藻门(3种)。全年总体优势种共有12种,均为蓝藻;春季优势种共有10种,其中假鱼腥藻(Pseudoanabaena sp.,0.327)和拉氏拟柱胞藻(Cylindrospermopsis raciborskii,0.122)优势度较高;夏季优势种共有13种,其中假鱼腥藻(0.199)的优势度最高;秋季优势种共有14种,其中隐球藻(Aphanocapsa sp.,0.160)和假鱼腥藻(0.143)的优势度较高;冬季优势种共有7种,其中隐杆藻(Aphanothece sp.,0.248)、隐球藻(0.122)和微小平裂藻(Merismopedia tenuissima,0.113)优势度较高;假鱼腥藻在全年占主要优势地位(表3)。
表3 东莞生态园浮游植物优势种优势度
Tab.3 Dominant species and dominance of phytoplankton in Dongguan Ecological Park
注:—表示该物种的优势度低于0.02。
Note:—indicates that the dominance of this species is less than 0.02.
门 phylum优势种 dominant species代号 code春季 spring夏季 summer秋季 autumn冬季 winter总体 total蓝藻门(Cyanophyta)阿氏颤藻(Oseillatoria agardhii)S10.0550.0650.026—0.030颤藻(Oscillatoria sp.)S20.052—0.032——平裂藻(Merismopedia sp.)S3——0.028——点形平裂藻(Merismopedia punctata)S4—0.0200.0420.0430.025微小平裂藻(Merismopedia tenuissima)S5———0.113—尖头藻(Raphidiopsis sp.)S60.0830.0460.0530.0200.049节旋藻(Arthrospira sp.)S7—0.034———拉氏拟柱胞藻(Cylindrospermopsis raciborskii)S80.1220.0720.034—0.056束球藻(Gomphosphaeria sp.)S9—0.0680.029—0.023水华束丝藻(Aphanizomenon flosaquae)S100.0430.0530.037—0.027依沙束丝藻(Aphanizomenon issatschenkoi)S110.0650.062——0.033假鱼腥藻(Pseudoanabaena sp.)S120.3270.1990.1430.0910.191鞘丝藻(Lyngbya sp.)S130.020————细鞘丝藻(Planktolyngbya subtilis)S140.0230.069——0.028小席藻(Phormidium tenue)S15—0.0480.0250.0310.021隐杆藻(Aphanothece sp.)S16———0.2480.034隐球藻(Aphanocapsa sp.)S170.060—0.1600.1220.070绿藻门(Chlarophyta)小球藻(Chlorella sp.)S18—0.0210.034——韦氏藻(Westella sp.)S19——0.034——狭形纤维藻(Ankistrodesmus angustus)S20—0.0310.032——
2.2.2 物种组成的时空差异 聚类分析显示,东莞生态园浮游植物群落可划分为3类,即春季(CⅠ)、冬季(CⅡ)和夏秋季(CⅢ)各聚为一类(图2)。ANOSIM分析表明,3个群落的物种组成有极显著性差异(R=0.763,P<0.01),且不同群落间均存在显著性差异(P<0.05)。NMDS分析显示,压力系数为0.04,表明数据拟合良好,上述群落划分可行。
SP—春季;SU—夏季;AU—秋季;WI—冬季。ML—月湖;CI—中央岛群,下同。
SP—spring;SU—summer;AU—autumn;WI—winter.ML—moon lake;CI—central islands,et sequentia.
图2 东莞生态园浮游植物物种聚类
Fig.2 Clustering of phytoplankton species in Dongguan Ecological Park
SIMPER百分率相似性分析显示,CⅠ、CⅡ和CⅢ群落内相似性分别为86.07%、87.02%和85.86%。CⅠ组中主要特征种为假鱼腥藻(贡献率为10.11%,丰度占32.69%)、拉氏拟柱胞藻(贡献率为8.90%,丰度占12.17%)和尖头藻(贡献率为8.40%,丰度占8.31%)等,特征种对组内平均相似性贡献率为90.40%;CⅡ组中主要特征种为隐杆藻(贡献率为7.87%,丰度占24.77%)、假鱼腥藻(贡献率为7.72%,丰度占9.82%)和隐球藻(贡献率为7.31%,丰度占12.25%)等,特征种对组内平均相似性贡献率为82.24%;CⅢ组中主要特征种为假鱼腥藻(贡献率为7.82%,丰度占18.02%)、尖头藻(贡献率为6.95%,丰度占4.85%)和拉氏拟柱胞藻(贡献率为6.73%,丰度占6.16%)等,特征种对组内平均相似性贡献率为84.33%。
2.2.3 浮游植物现存量的时空差异 2019年东莞生态园浮游植物年平均丰度为(16.83±1.97)×107 ind./L,年平均生物量为(23.00±3.00)mg/L。双因素方差分析表明,丰度(F=8.046,P<0.05)和生物量(F=11.708,P<0.01)均存在显著的季节性差异,但不存在显著的空间差异(丰度F=0.695,P>0.05;生物量F=0.002,P>0.05)。季节与湖区对丰度(F=3.533,P<0.05)和生物量(F=3.522,P<0.05)均存在显著交互作用(图3)。
标有不同字母者表示同一湖区不同季节间有显著性差异(P<0.05),标有相同字母者表示组间无显著性差异(P>0.05)。
The means with different letters in the same lake area are significant differences in different seasons at the 0.05 probability level,and the means with the same letter are not significant differences.
图3 东莞生态园浮游植物现存量的时空分布
Fig.3 Temporal and spatial distribution of phytoplankton standing crop in Dongguan Ecological Park
多重比较结果表明:月湖丰度夏季显著高于秋、冬季(P<0.05),春季丰度显著高于秋季(P<0.05),其他季节间均无显著性差异(P>0.05);中央岛群丰度夏、冬季显著高于春、秋季(P<0.05),其他季节间均无显著性差异(P>0.05)。月湖生物量夏季显著高于春、秋、冬季(P<0.05),其他季节间均无显著性差异(P>0.05);中央岛群生物量夏季显著高于春季(P<0.05),其他季节间均无显著性差异(P>0.05)。
2.3.1 功能群划分和优势功能群 东莞生态园浮游植物共划分成30个功能群,分别为A、B、C、D、E、F、G、H1、J、K、LM、Lo、M、MP、N、NA、P、S1、S2、SN、TB、TC、W0、W1、W2、X1、X2、X3、XPh和Y。功能群A、E、M和TB出现频率小于20%,为东莞生态园偶见或罕见功能群;功能群N、NA、X2、XPh和Y出现频率为20%~50%,为少见功能群;其余功能群出现率均大于50%,为常见功能群(表4)。D、F、H1、J、Lo、MP、P、S1、SN、W1、X1和Y为优势功能群,其中S1、D、P、SN和X1为全年的优势功能群(图4)。
表4 东莞生态园浮游植物功能群组成
Tab.4 Functional groups composition of phytoplankton in Dongguan Ecological Park
注:*表示优势功能群。 Note:* represents dominant functional groups.
功能群functional groups代表属(种)representative genus (species)生境特征habitat characteristics出现频率/%frequency of occurrenceF值F valueA根管藻属(Rhizosolenia)贫营养、洁净深水水体10.715B小环藻属(Cyclotella)中营养、中小型或大型浅水水体96.434C梅尼小环藻(Cyclotella meneghiniana)富营养、中小型水体78.575D∗针杆藻属(Synedra)、菱形藻属(Nitzschia)浑浊、营养盐丰富水体100.004E鱼鳞藻属(Mallomonas)贫营养、小型浅水水体14.293F∗卵囊藻属(Oocystis)、蹄形藻属(Krichneriella)中富营养、洁净、混合强水体89.293G空球藻属(Eudorina)、实球藻属(Pandorina)富营养、停滞水体57.140H1∗束丝藻属(Aphanizomenon)、鱼腥藻属(Anabaena)富营养、分层、含氮低浅水水体96.431J∗栅藻属(Scenedesmus)、四角藻属(Tetraedron)富营养、混合强浅水水体96.432K隐球藻属(Aphanocapsa)、隐杆藻属(Aphanothece)富营养浅水水体89.292LM束球藻属(Gomphosphaeria)富到超富营养、中小型水体71.430Lo∗多甲藻属(Peridinium)、拟多甲藻属(Peridiniopsis)贫到富营养、中到大型水体100.005M微囊藻属(Microcystis)富到超富营养、透明度高、小到中型水体3.570MP∗颤藻属(Oscillatoria)、舟形藻属(Navicula)经常性搅动、浑浊浅水水体100.004N鼓藻属(Cosmarium)持续或半持续混合水体46.435NA角星鼓藻属(Staurastrum)贫到中营养、低纬度净水水体28.574P∗脆杆藻属(Fragilaria)、直链藻属(Melosira)持续或半持续混合、富营养浅水水体96.435S1∗假鱼腥藻属(Pseudanabaena)混合浑浊、透明度低水体100.000S2螺旋藻属(Spirulina)温暖、高碱性浅水水体64.292SN∗尖头藻属(Raphidiopsis)、拟柱孢藻属(Cylindrospermopsis)温暖、混合水体100.000TB异极藻属(Gomphonema)强急流水体14.293TC鞘丝藻属(Lyngya)富营养、静水或流动缓慢、具有挺水植物水体92.862W0衣藻属(Chlamydomonas)富含有机质或水生生物腐败物水体100.002W1∗扁裸藻属(Phacus)、裸藻属(Euglena)有机污染浅水水体82.142W2囊裸藻属(Trachelomonas)、陀螺藻属(Strombomonas)中营养浅水水体78.573X1∗小球藻属(Chlorella)、纤维藻属(Ankistrodesmus)超富营养浅水水体100.003X2蓝隐藻属(Chroomonas)、金杯藻属(Kephyrion)中到富营养浅水水体32.143.5X3布纹藻属(Gyrosigma)、弓形藻属(Schroederia)贫营养、混合浅水水体64.294XPh网球藻属(Dictyosphaerium)含钙高、光照好、碱性小型水体50.003Y∗隐藻属(Cryptomonas)静水、牧食强度低水体42.863
图4 浮游植物优势功能群生物量占比的时空分布
Fig.4 Spatiotemporal distribution of relative biomass of dominant functional groups of phytoplankton
2.3.2 功能群的时空分布 春季,S1和SN功能群占比最大,分别为50.38%和12.42%;夏季,S1、 D和P功能群占比最大,分别为26.32%、21.07%和15.92%;秋季,S1、P和D功能群占比最大,分别为17.59%、15.82%和12.51%;冬季,S1、D和P功能群占比最大,分别为15.48%、12.29%和12.08%。春、夏、秋3个季节,月湖和中央岛群水系浮游植物优势功能群占比差别不大;冬季,月湖X1和F功能群占比较中央岛群高,分别为9.18%和6.50%,中央岛群Y和J功能群占比较月湖高,分别为13.00%和7.12%。因此,东莞生态园浮游植物功能群总体演替趋势为S1+SN+D+MP+H1(春季)→S1+D+P+X1+SN(夏季)→S1+P+D+X1+W1+H1(秋季)→S1+D+P+Y+Lo+J+X1+F(冬季)(图4)。
Q指数评价显示,S7站点为耐受水质,其余站点为中等水质,东莞生态园整体水质处于中等状态;营养状态指数评价显示,所有站点均处于中度富营养状态;H′和J指数评价显示,所有站点均处于轻污染状态;D指数评价显示,所有站点均处于中污染状态;浮游植物丰度和生物量评价显示,所有站点均处于重度富营养状态;叶绿素a评价显示,所有站点均处于中度富营养状态(表5)。
表5 各指数对东莞生态园的水质评价
Tab.5 Water quality evaluation in Dongguan Ecological Park according to different indices
站点site生态状况指数Q综合营养状态指数TLIShannon-Wiener指数H′Pielou指数JMargalef指数D丰度/(106 ind.·L-1)abundance生物量/(mg·L-1) biomass叶绿素a/(μg·L-1) chlorophyll-aVILWQVILWQVILWQVILWQVILWQVILWQVILWQVILWQS12.41中等64.38中度富营养2.74轻污染0.68轻污染2.89中污染142.2重富营养19.48重富营养36.83中度富营养S22.21中等63.77中度富营养2.64轻污染0.68轻污染2.56中污染168.6重富营养20.60重富营养30.97中度富营养S32.42中等64.51中度富营养2.73轻污染0.69轻污染2.77中污染203.9重富营养27.54重富营养37.44中度富营养S42.13中等64.23中度富营养2.79轻污染0.71轻污染2.75中污染176.1重富营养21.68重富营养36.04中度富营养S52.28中等65.67中度富营养2.84轻污染0.73轻污染2.53中污染201.4重富营养29.07重富营养39.19中度富营养S62.57中等62.82中度富营养2.76轻污染0.74轻污染2.10中污染118.2重富营养19.35重富营养31.26中度富营养S71.89耐受64.42中度富营养2.71轻污染0.71轻污染2.29中污染172.9重富营养22.24重富营养31.25中度富营养总体2.27中等64.26中度富营养2.74轻污染0.71轻污染2.55中污染169.1重富营养22.85重富营养34.71中度富营养
2.5.1 浮游植物优势种与环境因子的关系 RDA分析显示,轴1和轴2的特征值分别为0.266和0.157,物种与环境因子的相关性分别达到了0.951和0.924,共解释了东莞生态园浮游植物累积变异的72.95%。其中,和TP为显著影响浮游植物群落时空分布的关键环境因子(P<0.05);尖头藻(Raphidiopsis sp.)(S6)、水华束丝藻(Aphanizomenon flosaquae)(S10)、拉氏拟柱胞藻(S8)和阿氏颤藻(Oseillatoria agardhii)(S1)等大部分浮游植物均与水温、NH+4-N和呈正相关,微小平裂藻(S5)、隐球藻(S17)和隐杆藻(S16)与水温、NH+4-N和呈负相关,与TP呈正相关,颤藻(Oscillatoria sp.)(S2)和鞘丝藻(Lyngbya sp.)(S13)与pH呈正相关,平裂藻(Merismopedia sp.)(S3)、点形平裂藻(Merismopedia punctata )(S4)和微小平裂藻 (S5)与pH呈负相关,与TP呈正相关(图5(a))。
S1~S20为优势种代码,见表3。
The dominant species coded from S1 to S20 in the figure see Tab.3.
图5 浮游植物优势种和优势功能群与环境因子的RDA分析
Fig.5 RDA analysis between dominant species and dominant functional groups of phytoplankton and environmental factors
2.5.2 浮游植物优势功能群与环境因子的关系 RDA分析显示,轴1和轴2的特征值分别为0.270和0.163,优势功能群与环境因子的相关性分别达到0.89和0.79,共解释了东莞生态园优势功能群累积变异的87.72%。其中,和TP为显著影响优势功能群时空分布的关键环境因子(P<0.05);MP、SN、S1和H1功能群与WT、和pH呈正相关,与TP呈负相关,而Y、Lo、F和W1功能群则与之相反(图5(b))。
浮游植物的数量和种类组成是评价水体水质的重要指标之一。在本次对东莞生态园的调查中,共鉴定出137种浮游植物,隶属于7个门类,其中绿藻门物种数最多,不同季节共出现20种优势种,隶属于蓝藻门(17种)和绿藻门(3种),全年优势种共12种,均为蓝藻门,这一结果与广东省其他水库和湿地公园的情况类似[14-15,18-20]。这主要是因为东莞生态园位于亚热带地区,光照充足且水温较高,调查期间平均水温为27.66 ℃,有利于蓝藻生长[19],致使其全年占优势。蓝藻最适生长温度为25~35 ℃,对高温的耐受能力强于其他藻类[21]。此外,东莞生态园是拓塘建湖而成,湖泊的水源主要为地下水和大气降水,水质在成湖初期较好,但随着周边污水排放量增加和水体自净能力下降,湖水逐渐呈富营养化[10]。富营养化水体中,水文气象条件对浮游植物的影响更大,若氮磷含量较高,蓝藻则可成为全年的优势类群[11]。除蓝藻外,绿藻最适生长温度为18~30 ℃[22],绿藻对强太阳辐射和高温的适应能力也较强[18]。因此,在夏、秋两季,绿藻门的韦氏藻、小球藻和狭形纤维藻也能成为优势种。综合来看,东莞生态园属于蓝藻-绿藻型群落结构,是较典型的中富营养型水体[23]。
热带和亚热带水体浮游植物现存量丰水期通常高于枯水期[14]。东莞生态园全年处于富营养化状态,夏季平均水温为33.09 ℃,显著高于其他3个季节,导致夏季浮游植物丰度和生物量也显著高于其他3个季节。但处于丰水期的秋季浮游植物现存量低于夏、冬两季,可能是秋季有机污染严重,抑制了浮游植物的生长。东莞生态园中央岛群和月湖浮游植物丰度和生物量无显著性差异,可能主要与两个湖区有沟渠相连且四周均有排渠有关。东莞生态园中央岛群周边有新开渠、坑尾一渠排渠、南畲朗排渠,月湖周边有文庙排渠、下沙二渠、埔心排渠和大圳埔排渠,其中南畲朗排渠贯穿整个生态园,从而导致东莞生态园整体营养水平一致。值得注意的是,东莞生态园浮游植物丰度和生物量在季节和湖区间存在显著的交互作用,这可能是因为不同季节水体流速和营养盐的输入不同引起的各个湖区间营养盐的差异导致。
浮游植物功能群与水体环境相互响应,不同浮游植物功能群能够反映不同水体环境。本研究表明,东莞生态园共划分出30个功能群,其中功能群D、F、H1、J、Lo、MP、P、S1、SN、W1、X1和Y为优势功能群,东莞生态园处于富营养化状态,营养盐丰富,平均水深为3.31 m,平均透明度为0.34 m,导致喜生长于富营养水体、浅水水体和透明度较低水体的功能群成为优势功能群。其中,S1、D、P、SN和X1为全年的优势功能群,S1(27.34%)、D(14.89%)和P(12.53%)功能群全年生物量占比均较高,是东莞生态园最具代表性的功能群。以假鱼腥藻属为主的S1功能群代表混合、浑浊、透明度低水体,假鱼腥藻对温度的耐受性非常高[24],假鱼腥藻等是不具备或很少具有伪空泡的丝状藻类,不仅对低营养盐浓度的环境有较好的适应能力[25],在高氮磷比水体中也能成为优势种群[26],S1功能群在4个季节均为优势功能群且生物量占比均最高,与假鱼腥藻为全年最主要的优势种一致。D功能群代表浑浊、营养盐丰富水体,且以针杆藻属(Synedra)和菱形藻属(Nitzschia)为主,对水温适应性较广,东莞生态园夏季高温、多雨,大量降雨会使水体发生扰动,更适宜以耐冲刷、耐受低光照的D功能群藻类生长,故其在夏季的生物量占比较高。P功能群适宜生活在富营养浅水水体,具有广温性,且更适应于流动性水体,东莞生态园夏季多雨,导致底泥中营养物质上涌至水体,且流速较快,其中脆杆藻属(Fragilaria)和直链藻属(Melosira)由于单个生物量较大,往往能在水体中沉降下来,故其在夏季的生物量占比较高。春季总磷含量显著低于其他3个季节,拉氏拟柱胞藻与总磷呈负相关[27],故拉氏拟柱胞藻在春季成为主要优势种。夏季高温、多雨,水体发生扰动,营养盐丰富,导致喜富营养水体的浮游植物优势种和功能群迅速生长,X1功能群逐渐代替SN功能群成为优势功能群。秋季有机污染严重,导致对有机污染耐受程度较高的W1功能群成为优势功能群。冬季平均水温为17.76 ℃,水位低、流动性差,属于低温、静水环境,此外,冬季滤食性鱼类摄食率降低,对藻类的滤食作用减少,导致Y、F、J和Lo功能群形成优势[28]。
浮游植物群落结构特征及时空分布还受其他环境因子如透明度、降雨和多种营养盐含量的影响。本研究中,RDA分析显示,环境因子对功能群时空分布的解释率高于优势种,表明功能群对环境因子的响应比优势种更敏感,这与严广寒等[29]对洞庭湖浮游植物的研究结果一致。东莞生态园主要浮游植物优势种和优势功能群均与水温和TP呈正相关。
大量研究表明,水温是影响浮游植物群落结构的重要环境因子,水温能控制水体中各营养盐的溶解度,同时较高水温也能促进藻类的呼吸和光合作用[22]。王亚坤等[14]对松木山水库的研究表明,水温是影响浮游植物群落结构的主要环境因子。孙平宇等[15]对广州海珠国家湿地公园河湖连通水系的研究表明,影响浮游植物功能群的主要因子为水温。有研究表明,隐杆藻和隐球藻在20 ℃左右水温时能够成为优势种[14],与本研究中冬季优势种为隐杆藻和隐球藻的结果一致。本研究中,微小平裂藻与水温呈负相关,这也是导致冬季微小平裂藻为优势种的原因。本研究表明,Y、Lo、F和J功能群与水温呈负相关,与冬季优势功能群一致。由此可见,春、夏、秋季优势种和优势功能群向冬季演替的主要因子为水温。
本研究中,东莞生态园主要浮游植物种类和优势功能群与pH呈正相关,与TP呈负相关。东莞生态园水体pH平均值为8.15,水体偏碱性,蓝藻在偏碱性水体中适宜生长[30],主要是因为在弱碱性条件下,浮游植物更易获得碳源,有利于藻类进行光合作用。Reynolds等[6]指出,氮磷比是浮游植物生长的重要影响因子,当氮磷比为16时,浮游植物生长最适宜,当氮磷比小于16时,浮游植物生长受氮限制,反之受磷限制。东莞生态园的水系四季氮磷比均大于16,浮游植物生长受磷限制。颤藻、拉氏拟柱胞藻、假鱼腥藻和尖头藻等与pH呈正相关,与TP呈负相关;MP和SN功能群与pH呈正相关,与TP呈负相关。春季pH高于其他3个季节,而TP则低于其他3个季节,故拉氏拟柱胞藻在春季成为优势种,以尖头藻属和拟柱孢藻属为主的SN和MP功能群成为主要优势功能群。浮游植物对3种氮源和中NH+4的吸收速率更高[25]。本研究中,NH+4-N也是影响浮游植物的主要因子,隐球藻与NH+4-N呈显著负相关,秋季和冬季NH+4-N含量较低,故秋季和冬季隐球藻成为优势种。由此可见,春、夏、秋季优势种和功能群演替的主要因子为pH、TP和NH+4-N。
单一水质评价方法往往不能客观反映水体的真实营养状态,本研究中综合采用单因子评价方法、营养状态指数、功能群Q指数、浮游植物现存量、叶绿素a和3种多样性指数对东莞生态园水质进行了水质评价。浮游植物现存量评价显示,东莞生态园水质处于重富营养状态;单因子水质评价显示,该生态园水质为劣Ⅴ类;H′和J指数评价显示,该生态园水质处于轻污染状态;D指数评价显示,该生态园水质为中污染。有学者指出,利用浮游植物进行水质评价时,多样性指数比生物量评价更为准确[31],且H′和J两种多样性指数在水质状况较好的情况下更为适用,而D指数用来表征污染较为严重水体的水质状况较为理想[32]。本研究表明,东莞生态园浮游植物现存量、CODMn和TN含量均较高,因此,H′和J两种指数的轻污染评价结果明显与实际不符,不适合用于该园水质评价。单因子水质评价法虽然简单直观,但运用于水环境质量评价时常采用最差的单项指标来判定水体的优劣,评价结果可能会导致污染程度比实际情况严重[33],因此,这可能说明东莞生态园的水质并未达到最严重污染的级别,即单因子水质评价的劣Ⅴ类和现存量评价的最差级别水质也不太符合该园水质实际。综上,H′、J指数、浮游植物现存量及单因子水质评价方法均不适合东莞生态园的水质评价,而D丰富度指数方法可能较适合。本研究中更多结果表明,营养状态指数、功能群Q指数、叶绿素a和D指数的评价结果较为一致,这些指数更能反映东莞生态园水质现状,表明该园的水质处于中度富营养状态或受到了中度污染。
本研究中,单因子水质评价发现,东莞生态园整体水质为劣Ⅴ类,有机污染严重,主要污染物质为氮,这与刘黾等[10]在2012—2013年同样采用水化学方法对东莞生态园水质评价的结果一致,说明该园水质无明显变化。
本研究表明,东莞生态园有机污染严重,总氮为主要污染负荷,藻类异常增殖。在仍存在入园排渠的情况下,作者建议严格控制外来污染,加强水体自身净水能力,并提出一些保护、修复建议如下:1)控制外来污染,如加强生活废水和工厂污水的脱氮处理,减少营养盐排放,同时加强排渠湿地建设以起到对污染物的拦截、吸收和净化作用;2)采用生物操纵技术,如通过放养鲢、鳙和贝类,加强异常增殖藻类的生物控制作用;3)加强对公园湖泊岸带挺水植物及浅水区沉水植物等的培植,强化水草对污染物质的吸收净化作用;4)东莞生态园建设时间尚短,其水生态系统可能仍处在不稳定状态,应加强对其水生态环境的长期监测,并对生态环境问题进行系统研究,为东莞生态园水生态环境保护和修复提供更多技术支持。
1)东莞生态园浮游植物群落物种数以绿藻门为主,全年以蓝藻为优势种;浮游植物群落结构呈明显的季节性变化;浮游植物现存量在季节间存在显著性差异,但季节差异性在2个湖区的表现不尽一致。
2)东莞生态园浮游植物共划分为30个功能群,其中优势功能群12个,分别为D、F、H1、J、Lo、MP、P、S1、SN、W1、X1和Y。东莞生态园浮游植物优势功能群存在显著的季节差异,总体演替趋势为S1+SN+D+MP+H1(春季)→S1+D+P+X1+SN(夏季)→S1+P+D+X1+W1+H1(秋季)→S1+D+P+Y+Lo+J+X1+F(冬季)。
3)东莞生态园有机污染严重,主要污染物质为氮,多种水质评价结果显示,该园水质已呈中度富营养化或处于中污染状态。
4)浮游植物功能群比浮游植物优势种对环境因子更敏感,为显著影响浮游植物优势种和优势功能群的关键环境因子。
[1] 郭子良,张曼胤,崔丽娟,等.中国国家城市湿地公园的建设现状及其趋势分析[J].湿地科学与管理,2018,14(1):42-46.
GUO Z L,ZHANG M Y,CUI L J,et al.Current status and trends of development of national urban wetland parks in China[J].Wetland Science &Management,2018,14(1):42-46.(in Chinese)
[2] 张庆辉,赵捷,朱晋,等.中国城市湿地公园研究现状[J].湿地科学,2013,11(1):129-135.
ZHANG Q H,ZHAO J,ZHU J,et al.Current status of studies on urban wetland park in China[J].Wetland Science,2013,11(1):129-135.(in Chinese)
[3] 王勇,钟国,宋波,等.花溪十里河滩国家城市湿地公园鸟类多样性研究[J].贵州师范大学学报(自然科学版),2019,37(1):59-64.
WANG Y,ZHONG G,SONG B,et al.Preliminary study on Huaxi national urban wetland park bird diversity[J].Journal of Guizhou Normal University (Natural Sciences Edition),2019,37(1):59-64.(in Chinese)
[4] 李一帆,安建国,裴瑞华,等.东莞生态园湿地鱼类群落结构[J].水产学报,2013,47(5):131-141.
LI Y F,AN J G,PEI R H,et al.Structural attribute estimation of fish community in Dongguan Ecological Park wetland[J].Journal of Fisheries of China,2013,47(5):131-141.(in Chinese)
[5] 田琪,陈政.洞庭湖浮游植物群落结构调查与分析[J].内陆水产,2007,32(8):30-32.
TIAN Q,CHEN Z.Investigation and analysis of phytoplankton co- mmunity structure in Dongting Lake[J].Inland Fisheries,2007,32(8):30-32.(in Chinese)
[6] REYNOLDS C S,HUSZAR V,KRUK C,et al.Towards a functional classification of the freshwater phytoplankton[J].Journal of Plankton Research,2002,24(5):417-428.
[7] PADISK J,CROSSETTI L O,NASELLI-FLORES L.Use and misuse in the application of the phytoplankton functional classification:a critical review with updates[J].Hydrobiologia,2009,621(1):1-19.
[8] 赵耿楠.渭河干流及其秦岭北麓典型支流浮游藻类和固着藻类群落结构特征与环境分析[D].西安:西安理工大学,2020.
ZHAO G N.Community structure characteristics and environmental analysis of phytoplankton and periphyton in the sediment-laden Weihe River and its tributaries from the northern foot of the Qinling Mountains[D].Xi’an:Xi’an University of Technology,2020.(in Chinese)
[9] 郭芳,顾继光,赵剑,等.类群划分方法对南亚热带水库夏季浮游植物群落与环境响应关系的影响[J].环境科学,2020,41(11):5050-5059.
GUO F,GU J G,ZHAO J,et al.Effects of phytoplankton classifications on the relationship between phytoplankton community and environment in summer subtropical reservoirs,southern China[J].Environmental Science,2020,41(11):5050-5059.(in Chinese)
[10] 刘黾,陈金峰,刘燕,等.东莞生态园湿地水质模糊识别评价与对应分析[J].广东农业科学,2015,42(16):110-116.
LIU M,CHEN J F,LIU Y,et al.Fuzzy identification evaluation and correspondence analysis of water quality in Dongguan Eco-Park wetland[J].Guangdong Agricultural Sciences,2015,42(16):110-116.(in Chinese)
[11] 孙莹蓓,罗潋葱,韦桂峰,等.大沙河水库浮游植物种类组成及季节性变化[J].资源科学,2014,36(8):1732-1738.
SUN Y B,LUO L C,WEI G F,et al.Species composition and seasonal variation in phytoplankton at Dashahe Reservoir[J].Resources Science,2014,36(8):1732-1738.(in Chinese)
[12] 胡鸿钧,魏印心.中国淡水藻类:系统、生态及分类[M].北京:科学出版社,2006.
HU H J,WEI Y X.Freshwater algae in China: system,ecology and classification[M].Beijing: Science Press,2006.(in Chinese)
[13] HILLEBRAND H,DÜRSELEN C D,KIRSCHTEL D,et al.Biovolume calculation for pelagic and benthic microalgae[J].Journal of Phycology,1999,35(2):403-424.
[14] 王亚坤,汪洋,裴瑞华,等.广东东莞松木山水库浮游植物群落结构特征及其与环境因子的关系[J].大连海洋大学学报,2022,37(1):129-139.
WANG Y K,WANG Y,PEI R H,et al.Community structure and relationship with environmental factors of phytoplankton in Songmushan Reservoir in Dongguan City of Guangdong Province[J].Journal of Dalian Ocean University,2022,37(1):129-139.(in Chinese)
[15] 孙平宇,胡韧,欧林坚,等.广州海珠国家湿地公园河湖连通水系浮游植物功能群特征及其对环境因子的响应[J].生态学杂志,2023,42(11):2655-2664.
SUN P Y,HU R,OU L J,et al.Characteristics of phytoplankton functional groups and their responses to environmental factors in the connected river-lake system of Guangzhou Haizhu National Wetland Park [J].Chinese Journal of Ecology,2023,42(11):2655-2664.(in Chinese)
[16] 况琪军,马沛明,胡征宇,等.湖泊富营养化的藻类生物学评价与治理研究进展[J].安全与环境学报,2005,5(2):87-91.
KUANG Q J,MA P M,HU Z Y,et al.Study on the evaluation and treatment of lake eutrophication by means of algae biology[J].Journal of Safety and Environment,2005,5(2):87-91.(in Chinese)
[17] 娄保锋,周正,苏海,等.鄱阳湖营养水平关键指标时空分布特征及适宜控制标准[J].湖泊科学,2023,35(3):897-909.
LOU B F,ZHOU Z,SU H,et al.Temporal and spatial characteristics of key indicators of nutritional level and control standards in Lake Poyang[J].Journal of Lake Sciences,2023,35(3):897-909.(in Chinese)
[18] 王超,高越超,王沛芳,等.广东长潭水库富营养化与浮游植物分布特征[J].湖泊科学,2013,25(5):749-755.
WANG C,GAO Y C,WANG P F,et al.Eutrophication and characteristics of phytoplankton distributions in Changtan Reservoir,Guangdong Province[J].Journal of Lake Sciences,2013,25(5):749-755.(in Chinese)
[19] 何森华,丘泓欣,邓婉璐,等.广东孔江国家湿地公园水体中的浮游植物群落结构特征[J].湿地科学,2021,19(2):218-225.
HE S H,QIU H X,DENG W L,et al.Structure characteristics of phytoplankton community in the water of Guangdong Kongjiang National Wetland Park[J].Wetland Science,2021,19(2):218-225.(in Chinese)
[20] 朱明明,范存祥,吴中奎,等.海珠国家湿地公园浮游植物群落结构时空变化[J].生态学杂志,2020,39(5):1501-1508.
ZHU M M,FAN C X,WU Z K,et al.Spatial and temporal dynamics of phytoplankton community structure in Haizhu National Wetland Park[J].Chinese Journal of Ecology,2020,39(5):1501-1508.(in Chinese)
[21] NALEWAJKO C,MURPHY T P.Effects of temperature,and availability of nitrogen and phosphorus on the abundance of Anabaena and Microcystis in Lake Biwa,Japan:an experimental approach[J].Limnology,2001,2(1):45-48.
[22] 方丽娟,刘德富,杨正健,等.水温对浮游植物群落结构的影响实验研究[J].环境科学与技术,2014,37(sup 2):45-50.
FANG L J,LIU D F,YANG Z J,et al.Effects of water temperature on the phytoplankton community structure[J].Environmental Science &Technology,2014,37(sup 2):45-50.(in Chinese)
[23] ZOU H J,REN H U,HAN B P.Structure and dynamics of phytoplankton ccommunity in Hedi Reservoir,South China[J].Journal of Tropical and Subtropical Botany,2010,18(2):196-202.
[24] 韦洁琳,崔玉静,李亦真,等.南水北调东线山东境内调蓄湖库蓝藻种群分布特征及环境驱动因子[J].环境科学,2023,44(11):6137-6148.
WEI J L,CUI Y J,LI Y Z,et al.Distribution characteristics and environmental driving factors of cyanobacteria community in impounded lakes and reservoirs in Shandong on the east route of south-to-north water diversion project [J].Environmental Science,2023,44(11):6137-6148.(in Chinese)
[25] 袁一文,肖利娟,韩博平.热带大型水库蓝藻群落季节动态特征:以大沙河和高州水库为例[J].生态环境学报,2015,24(12):2027-2034.
YUAN Y W,XIAO L J,HAN B P.Seasonal dynamics of cyanobacteria assemblage in tropical large reservoirs,south China—using Dashahe and Gaozhou Reservoirs as examples[J].Ecology and Environmental Sciences,2015,24(12):2027-2034.(in Chinese)
[26] 王飞虎,杨越,张洋,等.黄河流域河南段浮游藻类功能群分布特征及影响因子[J/OL].水生生物学报,2023.https://kns.cnki.net/kcms/detail/42.1230.Q.20230327.1014.005.html.
WANG F H,YANG Y,ZHANG Y,et al.Distribution of phytoplankton functional groups and their relationship with environmental factors in the Yellow River basin within Henan Province [J/OL].Acta Hydrobiologica Sinica,2023.https://kns.cnki.net/kcms/detail/42.1230.Q.20230327.1014.005.html.(in Chinese)
[27] 杨祎明.拉氏拟柱孢藻对氮磷营养的响应及其生态优势研究[D].北京:中国科学院大学,2018.
YANG Y M.Study on the response of Paecilomyces laplacis to nitrogen and phosphorus nutrition and its ecological advantage[D].Beijing:University of Chinese Academy of Sciences,2018.(in Chinese)
[28] 李晓莉,陶玲,毛梦哲,等.温度和氮磷浓度对平裂藻和栅藻生长及竞争的影响[J].水生生物学报,2015,39(6):1217-1223.
LI X L,TAO L,MAO M Z,et al.The effect of the temperature and the N/P concentration on the growth and competition of Merismopedia and Scenedesmus[J].Acta Hydrobiologica Sinica,2015,39(6):1217-1223.(in Chinese)
[29] 严广寒,殷雪妍,汪星,等.基于种群、功能群对比分析洞庭湖浮游植物群落驱动因素及水质评价[J].环境科学,2023,44(11):6125-6136.
YAN G H,YIN X Y,WANG X,et al.Driving factors analyze of phytoplankton community by comparison of population and function groups and water quality evaluation in Dongting Lake[J].Environmental Science,2023,44(11):6125-6136.(in Chinese)
[30] 杨波,储昭升,金相灿,等.CO2/pH对三种藻生长及光合作用的影响[J].中国环境科学,2007,27(1):54-57.
YANG B,CHU Z S,JIN X C,et al.Influence of CO2/pH on growth and photosynthesis of three kinds of algae[J].China Environmental Science,2007,27(1):54-57.(in Chinese)
[31] 潘鸿,唐宇宏.草型浅水湖泊富营养化评价的浮游植物指标比较[J].安徽农学通报,2016,22(sup 1):67-69,110.
PAN H,TANG Y H.Comparison of eutrophication evaluation by phytoplankton indexes in a shallow aquatic vegetation lake[J].Anhui Agricultural Science Bulletin,2016,22(sup 1):67-69,110.(in Chinese)
[32] 余海霞,来勇,赵佳佳,等.城市河道浮游植物群落多样性指数对水质的指示作用[J].化学与生物工程,2016,33(9):32-36.
YU H X,LAI Y,ZHAO J J,et al.Water quality indication of phytoplankton community diversity index in urban rivers[J].Chemistry &Bioengineering,2016,33(9):32-36.(in Chinese)
[33] 李喆,姜作发,霍堂斌,等.黑龙江中游浮游植物多样性动态变化及水质评价[J].中国水产科学,2012,19(4):671-678.
LI Z,JIANG Z F,HUO T B,et al.Dynamics of phytoplankton diversity and water quality evaluation in the middle reaches of Heilongjiang River[J].Journal of Fishery Sciences of China,2012,19(4):671-678.(in Chinese)