淡水是一种有限且分布不均的资源,地表水是最重要的淡水来源[1]。在降水不规律或呈现季节性特征的干旱和半干旱地区,淡水分布不均且资源稀缺,建造水库以蓄积水资源是一种至关重要的解决方案。水库在为人类提供饮用水、农田灌溉和水力发电的同时,还承担着其他多重社会和生态功能。在社会功能方面,水库通过雨季蓄洪、旱季调水,有效缓解极端气候对地区安全的威胁;在生态功能方面,稳定的水域环境为动植物提供了良好的栖息场所,维护了区域生物多样性。由于气候变化和人类活动的干预,很多淡水水域的水质正在显著恶化[2-3],水体的物理、化学和生物特性发生了改变,影响了水生生物群落的结构和功能。浮游植物是水域生态系统初级生产者之一,具有种类多、繁殖快、生命周期短、生态转换效率高等特点,并通过食物网的能量流动影响水域生物的群落结构[4-5]。同时,浮游植物的种类组成、丰度等群落结构特征也受到水域环境条件的影响,如温度、光照和营养盐等[6-8]。已有研究表明,浮游植物的群落结构特征能够反映水体生态的健康状况[9]。
石佛寺水库位于辽宁省沈阳市沈北新区,始建于2003年5月,水库面积50 km2,其中生态面积22 km2,水域面积16 km2,平均水深6 m,是辽河干流上唯一一座控制性水利工程,也是大型河段型平原湿地水库。石佛寺水库主要以防洪和生态补水为核心功能,同时兼具休闲、游憩的功能,对辽河中下游起到了重要的生态调节作用。本研究分别于2024年3、7、9、12月对石佛寺水库及邻近辽河河段开展了浮游植物群落结构调查、水质监测,并对水质状况进行评价,以期为辽河水域有效管护,以及提高水域生态服务功能提供重要参考。
2024年3月29日、7月18日、9月24日和12月10日在石佛寺水库及邻近的辽河河段布设8个调查站位,其中水库设置6个站位(S1~S6),河段2个站位(S7、S8)。8个调查站位信息见表1、图1。
表1 石佛寺水库及邻近辽河河段调查站位坐标
Tab.1 Geographical coordinates of survey station in the Shifosi Reservoir and the adjacent sections of the Liaohe River
调查站位survey station东经east longitude北纬north latitudeS1123°29'28″42°09'32″S2123°28'46″42°09'08″S3123°27'52″42°09'36″S4123°27'11″42°09'40″S5123°27'34″42°10'32″S6123°26'20″42°10'25″S7123°25'58″42°09'57″S8123°29'20″42°11'09″
图1 石佛寺水库及邻近辽河河段调查站位
Fig.1 Survey station in the Shifosi Reservoir and the adjacent sections of the Liaohe River
1.2.1 浮游植物采集与分析 使用水生80型采水器在表层(水面下1 m)、底层(距水底1 m内)处采集水样,充分混合后,取1 L水样,立即加入鲁哥氏剂10 mL、混匀后,带回实验室,静置24 h,移去上清液,定容至50 mL,在电子显微镜(ECLIPSE 50i,Nikon)下,按照文献[10]鉴定浮游植物种类并统计数量。
1.2.2 水质样品采集与测定 现场使用水质分析仪(Quanta,美国哈希公司)测定溶解氧(DO)、水温(WT)、浊度(TUR)。使用水生80型采水器采集表层、底层水样,充分混合后取1 L水样带回实验室内,分别采用碱性过硫酸钾消解-紫外分光光度法、钼酸铵分光光度法、水杨酸分光光度法、重铬酸盐法、高锰酸钾氧化法测定水样总氮(TN)、总磷(TP)、氨氮(NH3-N)含量,以及化学需氧量(COD)和高锰酸盐指数(CODMn);采用分光光度法测定叶绿素a含量。
1.2.3 浮游植物优势度计算 优势度计算公式如下
Y=ni/N×fi。
(1)
式中:Y为优势度;ni为第i种浮游植物的个体数;fi为第i种在各站位出现的频率;N为样品中全部浮游植物的个体数。判别标准为:Y>0.02的种类定为优势种[11]。
1.2.4 群落特征指数计算 表2为基于浮游植物群落特征指数的水质评价标准[12]。浮游植物的Margalef丰富度指数(D)、Shannon-Wiener多样性指数(H′)、Pielou均匀度指数(J)计算公式如下
D=(S-1)/lnN,
(2)
(3)
J=H′/lnS。
(4)
表2 基于浮游植物群落特征指数的水质评价标准
Tab.2 Water quality evaluation criteria based on phytoplankton community characterization indices
类别categoryDH'J0~11~22~33~4>40~11~22~3>30~0.30.3~0.50.5~0.8判断标准evaluation criteria严重污染型α-中污染型β-中污染型寡污染型清洁型严重污染型α-中污染型β-中污染型清洁或轻污型重污染型中污染型清洁或轻污型
式中:S为样品中浮游植物种类数;N为样品中全部浮游植物的个体数;pi为第i种的个体数与总个体数的比值。
采用Excel 2016和Origin 2018软件进行数据整理、统计分析和图表绘制。采用SPSS 27.0软件进行月份和水体类型对浮游植物丰度影响的双因素方差分析,以P<0.05为显著水平。将数据进行对数转换后,采用SPSS 27.0软件对浮游植物群落特征指数与环境因子进行Pearson相关性分析,评估各变量间的相关性;将数据进行对数转换后,采用Canoco 5.0软件对浮游植物优势种丰度、环境因子进行去趋势对应分析(detrended correspondence analysis,DCA),如果第一轴排序梯度长度小于4,则采用冗余分析(RDA)方法,如果长度大于4,则采用典范对应分析(CCA)方法,以此探究浮游植物优势种丰度与环境因子的相关性[13]。
2.1.1 浮游植物物种组成及丰度 在石佛寺水库6个站位中共检测出浮游植物168种(属),分属于7个门。其中以硅藻门(Bacillariophyta)种类最多,有58种(属),绿藻门(Chlorophyta)有55种(属),蓝藻门(Cyanophyta)有25种(属),裸藻门(Euglenophyta)有15种(属),甲藻门(Pyrrophyta)有6种(属),隐藻门(Cryptophyta)有4种(属),金藻门(Chrysophyta)有5种(属)。在4次调查中,水库内发现的浮游植物种类数依次为70、115、121、71种(属)(表3)。
表3 石佛寺水库浮游植物优势种
Tab.3 Dominant phytoplankton species in the Shifosi Reservoir
门类 phylum种类 species3月 March7月 July9月 September12月 December蓝藻门(Cyanophyta)针晶蓝纤维藻镰刀型 (Dactylococcopsis rhaphidioides f.falciformis)+针晶蓝纤维藻 (Dactylococcopsis rhaphidioides)+隐藻门(Cryptophyta)尖尾蓝隐藻(Chroomonas acuta)+卵形隐藻(Cryptomonas ovata)+裸藻门(Euglenophyta)囊裸藻(Trachelomonas sp.)+旋转囊裸藻(Trachelomonas volvocina)+梅尼小环藻(Cyclotella meneghiniana)++硅藻门(Bacillariophyta)新星形冠盘藻小型变种(Stephanodiscus astraea)++舟形藻(Navicula sp.)++针形纤维藻(Ankistrodesmus acicularis)+球衣藻(Chlamydomonas globosa)+绿藻门(Chlorophyta)小球藻(Chlorella vulgaris)++++蛋白核小球藻(Chlorella protothecoides)++++绿球藻(Chlorella pyrenoidosa)+四尾栅藻(Scenedesmus quadricauda)+
注:+表示为优势种,下同。
Note:+means the dominant species,et sequentia.
在辽河河段2个站位中共检测出浮游植物114种(属),分属于7个门。其中硅藻门种类最多,有50种(属),绿藻门有36种(属),蓝藻门有14种(属),裸藻门有9种(属),甲藻门有1种(属),隐藻门有3种(属),金藻门有1种(属)。在4次调查中,河段内发现的浮游植物种类数依次为45、40、70、27种(属)。水库、河段浮游植物种类组成均以硅藻门最多,绿藻门其次,蓝藻门为第三位(表4)。
表4 辽河河段浮游植物优势种
Tab.4 Dominant phytoplankton species in the sections of Liaohe River
门类 phylum种类 species3月 March7月 July9月 September12月 December蓝藻门(Cyanophyta)小颤藻(Oscillatoria tenuis)+隐藻门(Cryptophyta)卵形隐藻(Cryptomonas ovata)+甲藻门(Pyrrophyta)薄甲藻(Glenodinium pulvisculus)+梅尼小环藻(Cyclotella meneghiniana)+++扭曲小环藻(Cyclotella comta)+颗粒直链藻(Melosira granulata)+硅藻门(Bacillariophyta)新星形冠盘藻小型变种(Stephanodiscus astraea)+舟形藻(Navicula sp.)++短小舟形藻(Navicula nivalis)+尖布纹藻(Cyrosigma acuminatum)+针形纤维藻(Ankistrodesmus acicularis)+球衣藻(Chlamydomonas globosa)++绿藻门(Chlorophyta)小球藻(Chlorella vulgaris)++++蛋白核小球藻(Chlorella protothecoides)++++四尾栅藻(Scenedesmus quadricauda)+
石佛寺水库浮游植物丰度在不同季节存在差异,总体表现为秋季>夏季>春季>冬季。年平均丰度为(610.54±115.26)×104 ind./L。辽河河段的浮游植物丰度变化趋势与水库相似,即在秋季达到峰值,冬季最低,年平均丰度为(604.59±173.64)×104 ind./L(图2)。
图2 浮游植物丰度季节变化
Fig.2 Seasonal variation of phytoplankton abundance
双因素方差分析结果显示,季节对浮游植物丰度具有极显著性影响(P<0.01),而地理位置(水库和河段)及两者的交互作用影响均不显著(P>0.05)。说明浮游植物丰度的变化主要受季节因素驱动,未受地理位置影响。
从群落组成来看,水库的浮游植物群落以硅藻门和绿藻门为主,两者丰度合计占浮游植物总丰度的70%以上。冬季硅藻门相对丰度相比其他季节为最高,达到41.96%,主要是新星形冠盘藻小型变种(Stephanodiscus astraea)数量显著增加,其相对丰度为52.73%。春季绿藻门相对丰度相比其他季节为最高,达到了60.49%,主要优势种为小球藻(Chlorella vulgaris)和蛋白核小球藻(Chlorella protothecoides),其中小球藻的相对丰度为25.59%,蛋白核小球藻的相对丰度为35.69%。春季蓝藻门相对丰度相比其他季节为最高,但总体占比较低,为14.05%。辽河河段与水库的群落组成总体一致,均以硅藻门和绿藻门为主,但河段中硅藻门的相对丰度占比相比水库要高。冬季绿藻门相比其他季节为最高,达到62.23%,优势种以小球藻和蛋白核小球藻为主,其中小球藻的相对丰度为38.20%,蛋白核小球藻的相对丰度为43.82%。夏季硅藻门的相对丰度相比其他季节为最高,达到45.91%,优势种以梅尼小环藻(Cyclotella meneghiniana)为主,其相对丰度为23.29%。春季蓝藻门的相对丰度相比其他季节为最高,达到5.73%。
2.1.2 浮游植物优势种 石佛寺水库和辽河河段浮游植物优势种结果见表3和表4。4次调查中,水库浮游植物优势种类数分别为6、7、7、4种(属);河段优势种类数分别为5、8、7、5种(属)。水库浮游植物优势种数量为7月和9月最多,河段浮游植物优势种数量7月最多。
浮游植物群落特征指数的季节变化如图3~图5所示。总体来看,石佛寺水库及辽河河段的Margalef丰富度指数(D)、Shannon-Wiener多样性指数(H′)和Pielou均匀度指数(J)均呈现出秋季最高、冬季最低的变化规律。水库年均D、H′和J分别为(4.58±1.45)、(2.76±0.51)和(0.69±0.04);辽河下游河段分别为(4.30±1.68)、(2.65±0.53)和(0.68±0.04);上游河段分别为(4.33±1.96)、(2.63±0.63)和(0.68±0.05),上下游河段数据较为接近。
图3 浮游植物丰富度指数(D)的季节变化
Fig.3 Seasonal variation of phytoplankton Margalef abundance index
图4 浮游植物多样性指数(H′)的季节变化
Fig.4 Seasonal variation of phytoplankton Shannon diversity index
图5 浮游植物均匀度指数(J)的季节变化
Fig.5 Seasonal variation of phytoplankton Pielou evenness index
基于浮游植物群落特征指数的水质评价标准,以D和H′指数为判定指标,水库和河段的水质在春、夏、秋季主要为清洁或轻污型,冬季则多为β-中污染型,上游河段冬季出现α-中污染型。以J指数判定,水库和河段全年均处于清洁或轻污水平。总体上,水库和河段在秋季水质状况最佳,冬季相对较差。
2.3.1 叶绿素a单因子评价 石佛寺水库叶绿素a含量季节变化如图6所示。总体表现为冬季最高,其次为秋季,春季最低。邻近水库的上、下游河段的季节变化趋势与水库不同,呈现为秋季最高、冬季次之、夏季最低。根据《地表水资源质量评价技术规程》(SL 395—2007)[14],叶绿素a值0.5~1.0 μg/L为贫营养,1.1~10.0 μg/L为中营养,10.1~26.0 μg/L为轻度富营养,26.1~160.0 μg/L为中度富营养,160.1~1 000.0 μg/L为重度富营养的标准,水库和河段均介于中营养到轻度富营养水平。
图6 叶绿素a含量的季节变化
Fig.6 Seasonal variation of chlorophyll a content
2.3.2 水质等级分析 环境因子季节变化见。石佛寺水库和辽河河段的高锰酸盐指数、总磷、总氮、化学需氧量、溶解氧、氨氮含量、浊度和水温存在时空差异。水库中靠近河段的S5、S6站位的总磷、总氮和氨氮浓度较高,并与河段站位S7、S8接近,显著高于水库站位S1~S4(P<0.05)。夏秋季水库总磷升高,冬季总氮、氨氮升高,表明河段输入可能是导致水库氮、磷含量高的主要因素。同时,S5、S6站位的浊度水平亦显著高于水库其他站位,说明水库对流域水质具有净化作用。总体上,四季S7、S8站位的浊度和营养盐水平显著高于水库站位,反映了河流携带陆源营养盐进入水库的特征。
依据《地表水环境质量标准》(GB 3838—2002)[15],以高锰酸盐指数为判定指标,水库水质处于Ⅱ~Ⅳ类水平;以化学需氧量为判定指标,处于Ⅰ~Ⅳ类水平;以溶解氧为判定指标,处于Ⅰ~Ⅲ类水平;以氨氮为判定指标,处于Ⅰ~Ⅳ类水平;以总磷为判定指标,处于Ⅱ~Ⅴ类水平;以总氮为判定指标,处于Ⅱ~劣Ⅴ类水平,局部达到劣Ⅴ类水平,表明水库水体存在一定氮污染。辽河上下游河段水质总体较差,以总氮为判定指标,均为劣Ⅴ类水平;其余指标多处于Ⅰ~Ⅲ类水平。总体而言,水库水质优于河段,主要污染因子为总氮和总磷。
2.4.1 Pearson相关性分析 石佛寺水库及辽河河段浮游植物与环境因子的Pearson相关性分析结果见表6和表7。石佛寺水库浮游植物丰富度、多样性和均匀度指数均与水温呈极显著正相关,与氨氮呈显著负相关,与浊度呈极显著正相关;多样性和丰富度指数与高锰酸盐指数呈显著正相关;多样性和丰富度指数与总氮呈极显著负相关;均匀度指数与总氮呈显著负相关。辽河河段浮游植物多样性、丰富度指数与总磷呈极显著正相关,丰度与总磷呈显著正相关;丰富度和多样性指数与水温呈显著正相关,均匀度指数与水温呈极显著正相关。
表5 水质参数季节变化
Tab.5 Seasonal variation of environmental factors
调查时间surveytime站位site高锰酸盐指数/(mg·L-1)CODMn化学需氧量/(mg·L-1)COD溶解氧/(mg·L-1)DO氨氮/(mg·L-1)NH3-N总磷/(mg·L-1)TP总氮/(mg·L-1)TN浊度/NTUTUR水温/℃WTS13.28±0.11BCd10.51±0.45Cd7.34±0.10Db0.133±0.002Dc0.061±0.002Da2.21±0.04Ea9.01±0.78Dc7.66±0.02FcS23.13±0.11CDd9.21±0.15Dc7.80±0.08Ba0.129+0.003Dc0.053±0.001Ec2.02±0.04Fb10.88±1.00Db7.78±0.02EcS33.11±0.05CDd9.48±0.16Dd7.20±0.05DEc0.129±0.002Dc0.071±0.001Cc2.11±0.03EFc26.01±2.55Cb8.01±0.02Cc3月S43.02±0.07Dd9.41±0.11Dc8.37±0.06Aa0.119±0.002Dc0.042±0.001Fc3.01±0.08Dc4.77±0.69Dc7.20±0.02HcMarchS53.44±0.09Bd10.89±0.16BCc5.62±0.05Gd1.049±0.023Aa0.100±0.003Ab5.11±0.06Ab37.85±2.20Bc8.69±0.01BcS63.35±0.08Bd10.51±0.22Cc6.50±0.08Fc0.710±0.023Ca0.090±0.002Bb4.31±0.18Bb11.35±1.24Db7.90±0.02DcS73.28±0.13BCc11.11±0.12Bd7.07±0.04Ec0.737±0.008Ba0.089±0.004Bb4.12±0.16Cb23.12±1.87Cc9.52±0.02AcS83.64±0.09Ab12.10±0.49Ad7.57±0.17Cab0.713±0.011Ca0.100±0.002Ac4.40±0.05Bc85.35±10.37Ac7.35±0.02GcS17.97±0.13Ac24.41±0.91Ab5.94±0.04Cc0.060±0.001Cd0.022±0.001Fb1.51±0.04Db25.29±2.48Fb27.78±0.05DaS26.09±0.12Ba19.12±0.48BCa6.81±0.05Bc0.071±0.001Ad0.021±0.001Fc0.82±0.03Ec2.18±0.24Hc27.79±0.06DaS34.98±0.11DEb16.42±0.36Eb6.90±0.06Bd0.068±0.002Bd0.053±0.002Eb2.41±0.07Cb9.59±1.42Gc27.49±0.06Ea7月S44.87±0.13Eb15.52±0.56EFa6.06±0.11Cc0.049±0.001Dd0.092±0.002Dc3.12±0.05ABc37.81±4.92Eb27.30±0.06FaJulyS55.51±0.10Cb19.80±0.54Ba6.92±0.10Bc0.071±0.001Ad0.173±0.002Ad3.03±0.07Bd322.42±1.63Ca27.74±0.10DaS64.62±0.10Fb14.79±0.25Fb7.51±0.13Ab0.069±0.001ABd0.109±0.004Cd3.02±0.04Bd74.89±0.45Da28.85±0.01AaS75.00±0.09DEa17.91±0.82Da7.61±0.04Ab0.069±0.001ABd0.150±0.004Bc3.21±0.08Ac1 036.28±5.34Aa27.99±0.05CaS85.10±0.08Da18.61±0.28CDa6.96±0.12Bc0.070±0.002Ad0.108±0.003Cd3.01±0.12Bd952.31±4.94Ba28.12±0.03BaS15.15±0.10BCb17.09±0.31Bc5.05±0.04Fd0.222±0.006Cb0.103±0.001Dc1.11±0.04Ec44.96±4.41Da17.27±0.07CbS25.01±0.07Cb15.01±0.34Db5.70±0.11Ed0.269±0.004Bb0.041±0.001Fb0.92±0.01Fc67.90±2.32Ca15.81±0.05FbS37.64±0.08Aa18.71±0.90Aa7.98±0.06Ba0.320±0.004Ab0.102±0.004Da0.51±0.02Gd70.22±5.96Ca17.39±0.05Bb9月S45.03±0.06Ca16.06±0.18Ca7.04±0.05Db0.209±0.005EFb0.088±0.002Eb3.43±0.04Db32.93±2.08Db16.12±0.06EbSeptemberS55.04±0.10Ca14.63±0.40Db8.12±0.07Ba0.209±0.003Ec0.161±0.005Bc4.49±0.15Bc112.14±10.23Ba16.14±0.02EbS65.26±0.06Bd16.20±0.71Ca8.01±0.09Ba0.214±0.003DEc0.141±0.002Cc4.11±0.05Cc70.59±8.20Ca16.45±0.06DbS75.01±0.01Ca14.91±0.52Db8.83±0.14Aa0.202±0.001Fc0.174±0.005Ab4.21±0.06Cb107.68±3.47Bb16.44±0.03DbS85.03±0.07Ca16.51±0.21BCb7.68±0.09Ca0.218±0.006CDc0.138±0.005Cb4.79±0.11Ab176.6±22.75Ab17.90±0.04AbS18.49±0.07Aa27.39±0.40Aa9.00±0.12Aa0.411±0.011Da0.069±0.001Bd0.90±0.03Dd4.71±0.54Cc2.55±0.01AdS24.71±0.06Bc19.21±0.66Ba7.02±0.15Eb0.290±0.005Ea0.022±0.001Ea3.61±0.17Ca14.10±0.22Cc2.55±0.01AdS33.79±0.01Ec14.49±0.31Dc7.81±0.05Bb0.663±0.017Aa0.081±0.003Aa5.81±0.07Aa7.63±0.62Cc0.45±0.01Dd12月S43.92±0.04Dc13.10±0.20Fb7.18±0.08DEb0.584±0.010Ca0.041±0.001Da5.52±0.06Ba58.02±3.02Ba1.39±0.01BdDecemberS53.79±0.08Ec13.91±0.41DEb7.38±0.18CDb0.657±0.012Ab0.070±0.002Ba5.61±0.24ABa6.11±0.73Cd0.43±0.01EdS64.26±0.07Cc15.73±0.42Ca7.79±0.20Ba0.273±0.007Eb0.069±0.003Ba5.70±0.07ABa7.22±0.88Cb0.45±0.01DdS73.65±0.07Fb12.21±0.25Gc7.47±0.19Cb0.668±0.019Ab0.081±0.002Aa5.71±0.07ABa6.87±0.35Cd0.58±0.01CdS83.72±0.07EFb13.40±0.21EFc7.39±0.08CDb0.617±0.011Bb0.061±0.001Ca5.58±0.13ABa94.99±5.50Ac0.41±0.01Fd
注:数据中标有相同大写字母代表同季节不同站位间无显著性差异(P>0.05),不同大写字母代表有显著性差异(P<0.05)。数据中标有相同小写字母代表同站位不同季节无显著性差异(P>0.05),不同小写字母代表同站位不同季节有显著性差异(P<0.05),下同。
Note:identical uppercase letters above the data indicate no significant difference (P>0.05) between different sites in a same season,while different uppercase letters indicate significant differences (P<0.05).Identical lowercase letters indicate no significant difference (P>0.05) at the same site across different seasons,while different lowercase letters indicate significant differences (P<0.05),et sequentia.
表6 石佛寺水库浮游植物群落特征指数与环境因子的Pearson相关系数
Tab.6 Pearson correlations of phytoplankton community indices with environmental factors in the Shifosi Reservoir
环境因子 environmental factor丰度 abundance丰富度指数 D均匀度指数 J多样性指数 H'氨氮NH3-N-0.140-0.435*-0.411*-0.462*总磷TP0.0090.2660.3500.312总氮TN-0.274-0.540**-0.424*-0.508**水温WT0.3510.773**0.774**0.809**浊度TUR0.2300.549**0.625**0.587**高锰酸盐指数CODMn0.2600.471*0.3770.436*
注:*表示在0.05水平上显著相关;**表示在0.01水平上显著相关,下同。
Note:* means significant correlation at the 0.05 level;** indicates significant correlation at the 0.01 level,et sequentia.
表7 辽河河段浮游植物群落特征指数与环境因子的Pearson相关系数
Tab.7 Pearson correlations of phytoplankton community indices with environmental factors in the sections of Liaohe River
环境因子 environmental factor丰度 abundance丰富度指数 D均匀度指数 J多样性指数 H'氨氮NH3-N-0.413-0.514-0.647-0.568总磷TP0.771*0.852**0.6770.853**总氮TN0.040-0.336-0.677-0.429水温WT0.4160.744*0.836**0.798*浊度TUR0.1980.2860.4540.338高锰酸盐指数CODMn0.5870.6190.6330.652
2.4.2 冗余分析(RDA) 根据DCA分析结果(表8),第一轴排序梯度长度小于4,采用冗余分析(RDA)探讨浮游植物与环境因子间的关系。石佛寺水库浮游植物优势种与环境因子的RDA分析结果如图7所示,第1、2排序轴的特征值分别为0.181、0.112,共解释了29.26%的物种信息,环境变量解释度为41.6%。水库最重要的环境影响因子为水温(p=0.002),其次为高锰酸盐指数(p=0.004)。球衣藻、囊裸藻和旋转囊裸藻的丰度与水温呈正相关关系,小球藻和蛋白核小球藻的丰度与水温呈负相关关系;针晶蓝纤维藻镰刀型和针晶蓝纤维藻的丰度与高锰酸盐指数呈负相关关系。
表8 浮游植物优势种与环境因子的DCA分析
Tab.8 DCA analysis results of dominant phytoplankton species and environmental factors
水体地点water location轴1梯度长度axis 1 gradient length水库 reservoir1.57邻近河段 river sections1.63
Dac f,针晶蓝纤维藻镰刀型;Dac,针晶蓝纤维藻;Chr,尖尾蓝隐藻;Cry,卵形隐藻;Tra,囊裸藻;Tra v,旋转囊裸藻;Cyc m,梅尼小环藻;Ste,新星形冠盘藻小型变种;Nav,舟形藻;Ank,针形纤维藻;Chl g,球衣藻;Chl v,小球藻;Chl p,蛋白核小球藻;Chl py,绿球藻;Sce,四尾栅藻。WT,水温;TUR,浊度;TP,总磷;TN,总氮;CODMn,高锰酸盐指数;NH3-N,氨氮。
Dac f,Dactylococcopsis rhaphidioides f.falciformis;Dac,Dactylococcopsis rhaphidioides;Chr,Chroomonas acuta;Cry,Cryptomonas acuta;Tra,Trachelomonas sp.;Tra v,Trachelomonas volvocina;Cyc m,Cyclotella meneghiniana;Ste,Stephanodiscus astraea;Nav,Navicula sp.;Ank,Ankistrodesmus acicularis;Chl g,Chlamydomonas globosa;Chl v,Chlorella vulgaris;Chl p,Chlorella protothecoides;Chl py,Chlorella pyrenoidosa;Sce,Scenedesmus quadricauda;WT,water temperature;TUR,turbidity;TP,total phosphorus;TN,total nitrogen;CODMn,permanganate index;NH3-N,ammonia nitrogen.
图7 石佛寺水库浮游植物优势种与主要环境因子的冗余分析
Fig.7 Redundancy analysis of phytoplankton dominant species and major environmental factors in Shifosi Reservoir
辽河河段浮游植物优势种与环境因子的RDA分析结果如图8所示,第1、2排序轴的特征值分别为0.315、0.275,共解释了58.95%的物种信息,环境变量解释度为90.6%。河段最重要的环境影响因子为高锰酸盐指数(p=0.012),其次为水温(p=0.04)和氨氮(p=0.05)。扭曲小环藻和尖布纹藻的丰度与高锰酸盐指数呈正相关关系,蛋白核小球藻的丰度与高锰酸盐指数呈负相关关系;短小舟形藻和小颤藻的丰度与水温呈正相关关系,小球藻的丰度与水温呈负相关关系;蛋白核小球藻的丰度与氨氮呈正相关关系,扭曲小环藻和尖布纹藻的丰度与氨氮呈负相关关系。
Osc,小颤藻;Cry,卵形隐藻;Gle,薄甲藻;Cyc m,梅尼小环藻;Cyc c,扭曲小环藻;Mel g,颗粒直链藻;Ste,新星形冠盘藻小型变种;Nav,舟形藻;Nav n,短小舟形藻;Cyr,尖布纹藻;Ank,针形纤维藻;Chl g,球衣藻;Chl v,小球藻;Chl p,蛋白核小球藻;Sce,四尾栅藻。WT,水温;TUR,浊度;TP,总磷;TN,总氮;CODMn,高锰酸盐指数;NH3-N,氨氮。
Osc,Oscillatoria tenuis;Cry,Cryptomonas acuta;Gle,Glenodinium pulvisculus;Cyc m,Cyclotella meneghiniana;Cyc c,Cyclo tella comta;Mel g,Melosira granulata;Ste,Stephanodiscus astraea;Nav,Navicula sp.;Nav n,Navicula nivalis;Cyr,Cyrosigma acuminatum;Ank,Ankistrodesmus acicularis;Chl g,Chlamydomonas globosa;Chl v,Chlorella vulgaris;Chl p,Chlorella protothecoides;Sce,Scenedesmus quadricauda;WT,water temperature;TUR,turbidity;TP,total phosphorus;TN,total nitrogen;CODMn,permanganate index;NH3-N,ammonia nitrogen.
图8 辽河河段浮游植物优势种与主要环境因子的冗余分析
Fig.8 Redundancy analysis of dominant phytoplankton species and major environmental factors in the section of Liaohe River
石佛寺水库浮游植物中,以硅藻门(58种,占34.52%)、绿藻门(55种,占32.74%)为优势门类,这与2016—2020年水库的调查结果一致[16],但优势种群的组成存在一定差异。分析原因应该与水域环境的年际变化有关。辽河河段同样以硅藻门(50种,占43.86%)、绿藻门(36种,占31.58%)为优势门类。水库和河段在种类组成及丰度上接近,这应该与水体连通有关。然而,从优势种及其相对丰度分布来看,二者存在一定差异,表明局部水文、水质条件(如营养盐浓度、流速)对群落结构具有调控作用。
石佛寺水库的总氮、总磷含量高于相关参考阈值(总氮为0.2 mg/L,总磷为0.02 mg/L)[17],表明水库水体营养盐含量较高,需注意发生水华的潜在风险。此外,浮游植物群落中的优势种群常作为评估水体营养状态的重要指标[18]。本次调查发现,水库中浮游植物优势种有球衣藻(Chlamydomonas globosa)、梅尼小环藻(Cyclotella meneghiniana)、四尾栅藻(Scenedesmus quadricauda)、尖尾蓝隐藻(Chroomonas acuta)、旋转囊裸藻(Trachelomonas volvocina)等,这些物种多常见于营养程度相对较高的水体,表明水库水质需持续关注[19]。
在生态学研究中,群落特征指数是评价浮游植物群落稳定性和水质状况的重要指标[20]。D、H′和J值越高,水域生态系统越稳定、水质越好[21]。生物群落特征指数高,表明水体生态系统的物种组成丰富且稳定性良好;如果指数低,则反映浮游生物群落结构趋于简单化,生态系统稳定性降低。综合D、H′指数可见,春、夏季石佛寺水库、辽河上下游河段水质整体表现为轻度污染状态,冬季石佛寺水库、辽河上下游河段均达到中度污染状态,秋季则全部为清洁或轻度污染状态。分析原因为农业活动在春夏季较为频繁,且降雨丰富,陆源氮、磷等营养盐随降雨汇入水体,加上适宜的温度,浮游植物快速生长,一直持续到秋季。冬季为辽河的枯水期,同时低温也抑制了浮游植物的繁殖、生长,进而降低了浮游植物的多样性。此外,冬季H′结果显示,上游河段为α-中污染型,而下游河段为β-中污染型,上游污染程度高于下游。这种差异表明水库作为缓冲区起到了净化水质的作用,阻隔了污染物向下游的扩散。然而,水库的存在虽有助于改善水质,但在水库中积累了大量营养盐,增加了汛期发生内源释放且引发二次污染的风险。
叶绿素a是反映水域环境中浮游植物生物量的重要指标。夏季水库叶绿素a含量显著高于河段,可能与夏季汛期水动力条件差异导致的浊度变化有关。调查结果显示,河段水体浊度约为水库水体的8.1倍,导致光照受限,影响藻类光合作用的高效进行;而水库相对稳定的水文条件则更有利于浮游植物的生长。
浮游植物群落结构受水体生态功能和地理特征共同制约,其时空动态则受到水文水质、营养盐供给、气候变化及人类活动等多种因素的综合影响[22]。石佛寺水库位于中国北方地区,水温季节性差异显著,因此水温成为浮游植物生长的重要环境影响因子[23-24]。该水库浮游植物与环境因子的RDA结果显示,球衣藻、囊裸藻和旋转囊裸藻等多种藻类的丰度与水温呈正相关关系,表明较高的水温有利于浮游植物的生长与代谢活动。小球藻和蛋白核小球藻在四季均为优势种,说明其具有较强的环境适应能力和稳定的生态位。RDA结果显示,这两种藻类的丰度与水温呈负相关,表明其丰度在低温季节相对更高。这种现象并非其最适生长温度较低,而是因为其具有较强的温度适应性和环境耐受力,使其在气温较低、其他藻类生长受限时仍能保持较高丰度,从而呈现出相对优势。作为浮游植物生长和代谢的关键限制因子,氮元素在藻类物质合成与能量传递过程中发挥着不可替代的作用。研究表明[25],总氮水平与浮游植物群落的演替过程及空间分布密切相关。Pearson相关性分析显示,石佛寺水库浮游植物群落特征指数与总氮呈负相关关系,这与已有研究的结论一致,说明水库营养盐含量高,氮元素已不是浮游植物生长的限制因子[26-28]。已有研究表明,总磷也是影响浮游植物群落结构的关键因子[29-33],在对南京玄武湖[34]和白洋淀水域生态系统[35]的研究中,总磷浓度被证实与浮游藻类生物量存在正相关关系。此外,也有研究认为氨氮是影响浮游植物的重要因子[36-38],本研究亦得出相同结论。浊度反映了水体中悬浮物对光线的阻碍程度,与水体中泥沙量、浮游生物量有关[39],已有研究表明,浊度为影响浮游植物群落的关键因子[40-41]。高锰酸盐指数作为综合反映水体中有机物和还原性无机物污染程度的指标,通常与浮游植物生长呈正相关关系[42]。相比水库,河段水流速度较快,陆源输入的营养盐和有机物易被稀释,局部影响相对较小,因此河段浮游植物丰度与高锰酸盐指数未表现出显著相关性。在一些湖泊、水库生态系统中,通常浮游植物的丰度与水体中氮、磷等营养盐的含量呈显著正相关关系,反映了浮游植物生长受环境因子制约的基本规律[43],然而,在石佛寺水库的相关性分析中,并未发现浮游植物丰度与营养盐之间存在正相关关系,这可能与水库每年增殖放流数十万尾鲢鳙苗种有关。鲢鳙对浮游植物的摄食压力,影响了浮游植物的群落结构,扰动了浮游植物丰度与营养盐之间的相关性。
1)石佛寺水库和辽河河段浮游植物均以硅藻门和绿藻门种类为主,多样性较高,群落结构表现出季节差异。
2)石佛寺水库和辽河河段水质总体为清洁或轻污染水平,仍有发生水华的风险,需加强外源污染管控并提升水体自净能力。
3)通过Pearson分析、RDA分析得出,水温、总氮、浊度、高锰酸盐指数和氨氮是影响石佛寺水库浮游植物群落结构的主要环境因子;水温、氨氮、高锰酸盐指数和总磷是影响辽河河段浮游植物群落结构的主要环境因子。
4)水库水质优于邻近的辽河河段,水库生态环境对水质调控发挥了重要作用。
[1] SHIKLOMANOV AND JOHN C RODDA E B I A.World Water Resources at the Beginning of the Twenty-First Century[M].Cambridge:Cambridge University Press,2003.
[2] LOUCIF K,NEFFAR S,MENASRIA T,et al.Physico-chemical and bacteriological quality assessment of surface water at Lake Tonga in Algeria[J].Environmental Nanotechnology,Monitoring &Management,2020,13:100284.
[3] BETIKU O C,SARJEANT K C,NGATIA L W,et al.Evaluation of microbial diversity of three recreational water bodies using 16S rRNA metagenomic approach[J].Science ofthe Total Environment,2021,771:144773.
[4] LI Y,FANG L C,CAO G R,et al.Reservoir regulation-induced variations in water level impacts cyanobacterial bloom by the changing physiochemical conditions[J].Water Research,2024,259:121836.
[5] LIU B Q,CHEN M R,WANG C S.The dynamics of trophic cascades on phytoplankton induced by mesozooplankton in coastal water,Daya Bay,northern South China Sea[J].Microorganisms,2024,12(10):1982.
[6] CHEN T T,ZHANG Y X,SONG S Q,et al.Diversity and seasonal variation of marine phytoplankton in Jiaozhou Bay,China revealed by morphological observation and metabarcoding[J].Journal of Oceanology and Limnology,2022,40(2):577-591.
[7] RAJAPAKSHA R P,WU M L,WANG Y T,et al.Long-term alterations of nutrient dynamics and phytoplankton communities in Daya Bay,South China Sea[J].Marine Pollution Bulletin,2024,208:116955.
[8] SUN R,WEI J L,ZHANG S S,et al.The dynamic changes in phytoplankton and environmental factors within Dongping Lake (China) before and after the South-to-North Water Diversion Project[J].Environmental Research,2024,246:118138.
[9] HU X Y,HU M,ZHU Y,et al.Phytoplankton community variation and ecological health assessment for impounded lakes along the eastern route of China’s South-to-North Water Diversion Project[J].Journal of Environmental Management,2022,318:115561.
[10] 胡鸿钧,魏印心.中国淡水藻类:系统、分类及生态[M].北京:科学出版社,2006.HU H J,WEI Y X.The freshwater algae of China[M].Beijing:Science Press,2006.(in Chinese)
[11] 罗民波,陆健健,王云龙,等.东海浮游植物数量分布与优势种[J].生态学报,2007,27(12):5076-5085.LUO M B,LU J J,WANG Y L,et al.Horizontal distribution and dominant species of phytoplankton in the East China Sea[J].Acta Ecologica Sinica,2007,27(12):5076-5085.(in Chinese)
[12] 王国涛,陈斌斌,王敏,等.浙南大罗山天河水库浮游植物群落结构的季节变化及对水质的指示作用[J].应用生态学报,2021,32(6):2227-2240.WANG G T,CHEN B B,WANG M,et al.Seasonal variation of phytoplankton community structure in the Tianhe Reservoir of the Daluoshan Mountain in southern Zhejiang Province,China and its indication to water quality[J].Chinese Journal of Applied Ecology,2021,32(6):2227-2240.(in Chinese)
[13] 刘川,闫金霞,刘家乐,等.青龙湖浮游植物群落与环境因子相关性分析[J].人民黄河,2021,43(7):104-109.LIU C,YAN J X,LIU J L,et al.Correlation analysis of phytoplankton community and environmental factors in Qinglong Lake[J].Yellow River,2021,43(7):104-109.(in Chinese)
[14] 中华人民共和国水利部.地表水资源质量评价技术规程:SL 395—2007[S].北京:中国水利水电出版社, 2007.Ministry of Water Resources of the People’s Republic of China.Technical regulations for surface water resource quality assessment:SL 395-2007[S].Beijing:China Water Resources and Hydropower Press, 2007.(in Chinese)
[15] 国家环境保护总局,国家质量监督检验检疫总局.地表水环境质量标准:GB 3838—2002[S].北京:中国环境科学出版社,2002.State Environmental Protection Administration of the People’s Republic of China,General Administration of Quality Supervision,Inspection and Quarantine of the People’s Republic of China.Environmental quality standards for surface water:GB 3838—2002[S].Beijing:China Environmental Science Press,2002.(in Chinese)
[16] 孙婧妍.石佛寺水库浮游植物时间分布相关研究[J].地下水,2022,44(2):89-91.SUN J Y.Study on temporal distribution of phytoplankton in shifosi reservoir[J].Ground Water,2022,44(2):89-91.(in Chinese)
[17] 李如忠,刘科峰,钱靖,等.合肥市区典型景观水体氮磷污染特征及富营养化评价[J].环境科学,2014,35(5):1718-1726.LI R Z,LIU K F,QIAN J,et al.Nitrogen and phosphate pollution characteristics and eutrophication evaluation for typical urban landscape waters in Hefei City[J].Environmental Science,2014,35(5):1718-1726.(in Chinese)
[18]
V.Rotifers as indicators of water quality[J].Hydrobiologia,1983,100(1):169-201.
[19] 袁晓泉,周胜福,杨秋菊,等.抚河浮游生物多样性和功能群组成分析[J].江西水产科技,2024(1):34-40.YUAN X Q,ZHOU S F,YANG Q J,et al.Analysis of plankton diversity and functional group composition in Fuhe River[J].Jiangxi Fishery Science and Technology,2024(1):34-40.(in Chinese)
[20] 刘钢,孟云飞,吴丹,等.青藏高原可鲁克湖浮游动物群落结构特征及水质评价[J].大连海洋大学学报,2018,33(3):379-386.LIU G,MENG Y F,WU D,et al.Characteristics of zooplankton community structure and water quality assessment in Keluke Lake in the Qinghai-Tibetan Plateau[J].Journal of Dalian Ocean University,2018,33(3):379-386.(in Chinese)
[21] 白海锋,孔飞鹤,王怡睿,等.北洛河流域浮游动物群落结构时空特征及其与环境因子相关性[J].大连海洋大学学报,2021,36(5):785-795.BAI H F,KONG F H,WANG Y R,et al.Spatio-temporal characteristics and relationship with environmental factors of zooplankton community structure in the Beiluo River basin[J].Journal of Dalian Ocean University,2021,36(5):785-795.(in Chinese)
[22] MA C X,YU H X.Phytoplankton community structure in reservoirs of different trophic status,Northeast China[J].Chinese Journal of Oceanology and Limnology,2013,31(3):471-481.
[23] KASPRZAK P,LATHROP R C,CARPENTER S R.Influence of different sized Daphnia species on chlorophyll concentration and summer phytoplankton community structure in eutrophic Wisconsin lakes[J].Journal of Plankton Research,1999,21(11):2161-2174.
[24] 邵海燕,王卿,高春霞,等.长江口浮游植物群落特征及影响因素分析[J].大连海洋大学学报,2024,39(1):124-133.SHAO H Y,WANG Q,GAO C X,et al.Analysis of phytoplankton community characteristics and influencing factors in the Yangtze River Estuary[J].Journal of Dalian Fisheries University,2024,39(1):124-133.(in Chinese)
[25] 叶琳琳,吴晓东,刘波,等.太湖西北湖区浮游植物和无机、有机氮的时空分布特征[J].湖泊科学,2017,29(4):859-869.YE L L,WU X D,LIU B,et al.Temporal and spatial distributions of phytoplankton and inorganic and organic nitrogen in northwest region of Lake Taihu[J].Journal of Lake Sciences,2017,29(4):859-869.(in Chinese)
[26] 蒋玉玲,李青,陈晓宏,等.广州市大金钟水库浮游植物与环境因子的关系[J].水资源保护,2011,27(1):46-50.JIANG Y L,LI Q,CHEN X H,et al.Relationship between phytoplankton and environmental factors in Dajinzhong reservoir of Guangzhou[J].Water Resources Protection,2011,27(1):46-50.(in Chinese)
[27] 孙莹蓓,罗潋葱,韦桂峰,等.大沙河水库浮游植物种类组成及季节性变化[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)
[28] 王亚坤,汪洋,裴瑞华,等.广东东莞松木山水库浮游植物群落结构特征及其与环境因子的关系[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)
[29] 沈会涛,刘存歧.白洋淀浮游植物群落及其与环境因子的典范对应分析[J].湖泊科学,2008,20(6):773-779.SHEN H T,LIU C Q.Canonical correspondence analysis of phytoplankton community and its environmental factors in the Lake Baiyangdian[J].Journal of Lake Sciences,2008,20(6):773-779.(in Chinese)
[30] 杨丽标,韩小勇,孙璞,等.巢湖藻类组成与环境因子典范对应分析[J].农业环境科学学报,2011,30(5):952-958.YANG L B,HAN X Y,SUN P,et al.Canonical correspondence analysisof algae community and its environmental factors in the Lake Chaohu,China[J].Journal of Agro-Environment Science,2011,30(5):952-958.(in Chinese)
[31] LALLY H,GORMALLY M,HIGGINS T,et al.Phytoplankton assemblages in four wetlands created on cutaway peatlands in Ireland[J].Biology and Environment:Proceedings of the Royal Irish Academy,2012,112B(2):207-216.
[32] 师艳丽,张依章,闫祯,等.天津市近岸海域秋季浮游动物群落结构与环境因子相关性[J].环境工程技术学报,2022,12(6):1995-2001.SHI Y L,ZHANG Y Z,YAN Z,et al.Zooplankton community structure and its relationship with environmental factors in autumn in Tianjin coastal marine areas[J].Journal of Environmental Engineering Technology,2022,12(6):1995-2001.(in Chinese)
[33] 郭坤,彭婷,罗静波,等.长湖浮游动物群落结构及其与环境因子的关系[J].海洋与湖沼,2017,48(1):40-49.GUO K,PENG T,LUO J B,et al.Community structure of zooplankton and the driving physicochemical factors in Changhu lake[J].Oceanologia et Limnologia Sinica,2017,48(1):40-49.(in Chinese)
[34] 付显婷,杨薇,赵彦伟,等.白洋淀浮游动物群落结构与水环境因子的关系[J].农业环境科学学报,2020,39(6):1271-1282.FU X T,YANG W,ZHAO Y W,et al.Community structure of zooplankton and its relationship with water environmental factors in Baiyangdian Lake[J].Journal of Agro-Environment Science,2020,39(6):1271-1282.(in Chinese)
[35] 梅彦昆,苏雅玲,董志国,等.富营养化浅水湖泊生态修复对水质与含碳温室气体通量的影响:以南京玄武湖为例[J].湖泊科学,2025,37(6):1978-1987.MEI Y K,SU Y L,DONG Z G,et al.Impacts of ecological restoration on water quality and carbon-containing greenhouse gas fluxes in a shallow,eutrophic lake:a case study of Lake Xuanwu,Nanjing[J].Journal of Lake Sciences,2025,37(6):1978-1987.(in Chinese)
[36] 韦梦琳,李法云,洪天宇,等.上海鹦鹉洲湿地与外围河道浮游植物群落时空差异及其影响因子[J].湖泊科学,2025(2):429-445.WEI M L,LI F Y,HONG T Y,et al.Spatio-temporal differences of phytoplankton communities and their driving factors in the Yingwuzhou Wetland and its surrounding canals,Shanghai[J].Journal of Lake Sciences,2025(2):429-445.(in Chinese)
[37] 于洪贤,曲翠,马成学.牡丹江浮游植物丰度与环境因子的相关性分析[J].湿地科学,2008,6(2):293-297.YU H X,QU C,MA C X.Correlation between phytoplankton abundance and its environmental factors in Mudan River[J].Wetland Science,2008,6(2):293-297.(in Chinese)
[38] 李利强,黄代中,熊剑,等.洞庭湖浮游植物增长的限制性营养元素研究[J].生态环境学报,2014,23(2):283-288.LI L Q,HUANG D Z,XIONG J,et al.Nutrient limiting phytoplankton growth in Dongting Lake[J].Ecology and Environment Sciences,2014,23(2):283-288.(in Chinese)
[39] SILVOSO J,IZAGUIRRE I,ALLENDE L.Picoplankton structure in clear and turbid eutrophic shallow lakes:a seasonal study[J].Limnologica,2011,41(3):181-190.
[40] ZHANG Y X,YU H X,LIU J M,et al.Investigating the differences in driving mechanisms for phytoplankton community composition under various human disturbances in cold regions[J].Journal of Cleaner Production,2024,461:142686.
[41] 周铭浩,李彤,邓滢,等.不同水期东江干流浮游植物群落结构特征及其与水环境因子的关系[J].中国环境监测,2025,41(1):121-133.ZHOU M H,LI T,DENG Y,et al.Phytoplankton community structural characteristics and its relationship with water environmental factors in the main stream of Dongjiang River during different water periods[J].Environmental Monitoring in China,2025,41(1):121-133.(in Chinese)
[42] 白海锋,王怡睿,宋进喜,等.渭河陕西段浮游植物群落结构时空变化与影响因子分析[J].环境科学学报,2021,41(8):3290-3301.BAI H F,WANG Y R,SONG J X,et al.Spatio-temporal characteristics and influencing factors of phytoplankton community structure in the Shaanxi Section of Weihe River,China[J].Acta Scientiae Circumstantiae,2021,41(8):3290-3301.(in Chinese)
[43] KASPRZAK P,PADIS
K J,KOSCHEL R,et al.Chlorophyll a concentration across a trophic gradient of lakes:an estimator of phytoplankton biomass [J].Limnologica,2008,38(3/4):327-338.