1. 引言

牛口峪水库位于北京市房山区，原为燕山石化储存和处理炼化装置工业用水的水库，燕山石化公司对园内约8万平方米已停用的水库进行了改造修复，种植了芦苇(Phragmitesaustralias)、香蒲(Typhaangustafolia)、穗状狐尾藻(Myriophyllumspicatum)、篦齿眼子菜(Potamogetonpectinatus)等具有高净化能力和高耐受力的水生植物，来实现深度净化工业外排废水。此外，该水库还交替错落地建设了不同密度的湿草甸和植有漂浮水生植物的自由水面，形成了十分理想的鸟类栖息地 [1] 。2019年7月，牛口峪水库建成了“保育区”，面积达10万m²，用以保护栖息在水库内的各类野生动物。水体的理化因子和浮游生物群落是水生态系统的重要组成部分，共同在能量流动、物质循环和信息联系中发挥作用，可以反映水体富营养化状态和生物多样性状况 [2] [3] [4] [5] 。

该水库是北京市实现了污水达标排放领域内生态环境恢复与保护的有机结合的范例 [1] ，但截止目前，关于牛口峪水库的浮游生物的物种多样性及其环境评价报道较少，因此，对其开展水体理化因子和浮游生物群落研究，对于北京市其他地区污水处理和生态修复具有重要的指导意义。

2. 材料和方法

2.1. 采样时间与样站的设置

根据2017~2021年北京市水生野生动植物救护中心对牛口峪水库实施的监测规划，本团队分别于2017年6月(夏季)、9月(秋季)，2018年6月(夏季)、9月(秋季)，2021年5月(夏季)对牛口峪水库进行了水体理化因子和浮游生物采样调查，水体理化因子采集参照《养殖水环境化学实验》 [6] 进行，浮游生物采集方法参照《水生生物学》 [7] 进行。2017，2018年采样站点为B，2021年采样站位为A、B和C，详见图1。

图1. 牛口峪水库采样样站设置

2.2. 测定方法

2.2.1. 理化参数测定

对牛口峪水库水体物理特征(温度、透明度)、化学特征[pH、溶氧(DO)、叶绿素a (Chl-a)、化学需氧量(COD_Mn)、总氮(TN)、总磷(TP)、氨氮( ${\text{NH}}_4^{\text{+}}\text{-N}$)、亚硝酸氮(NO₂-N)、硝酸氮(NO₃-N)、活性磷(ADP, PO₄-P)、碱度(ALK)、硬度(TH)、氯离子( ${\text{Cl}}^{-}$)、硫酸根离子( ${\text{SO}}_4^{2-}$)、钾离子(K⁺)、钠离子(Na⁺)]等进行了测定，测定方法按照《水和废水监测分析方法》 [8] 和《养殖水环境化学实验》 [6] 进行。

2.2.2. 浮游植物样品采集

每个样站取水5 L (用水生-80型有机玻璃采水器的固定容积)，将采集后的水放入水桶中混合均匀。另取1 L塑料瓶，从混合水样中取1 L水作为浮游植物样品，立即加入15 mL鲁哥氏液固定并混合均匀。样品在实验室内放置48 h后用虹吸管进行浓缩，装入小塑料瓶中定容至50 mL后进行定性、定量分析。分析前需充分摇匀样品，后立即从中吸取0.1 mL混合液小心置于0.1 mL浮游植物计数框中，在15 × 40倍显微镜(麦克奥迪BA210)下进行观察计数。每瓶样品计数两片，取平均值。浮游植物参照《中国淡水藻类——系统、分类及生态》 [9] 、《中国淡水藻志》 [10] 等文献进行定性分析，定量分析参照《水生生物学》 [7] 提供的方法执行。

2.2.3. 浮游动物样品采集

小型浮游动物(原生动物和轮虫)的定量采用浮游植物定量的浓缩水样浓缩至50 ml，测定时充分摇匀，用定量吸管准确吸取1 ml置于浮游动物计数框内，并都在麦克奥迪BA210型显微镜下进行计数。大型浮游动物(枝角类、桡足类等)采水50 L，经过25#浮游生物网(300目)过滤，浓缩液用5%甲醛固定，测定时吸取1 ml在显微镜下计数，然后根据镜检的种类和个数来计算密度和生物量。浮游动物的种类鉴定参照《水生生物学》《中国动物志 节肢动物门 甲壳纲 淡水枝角类》《中国动物志 节肢动物门 甲壳纲 淡水桡足类》《中国淡水轮虫志》 [7] [11] [12] 著作进行鉴定，浮游动物的密度、生物量按照《湖泊富营养化调查规范》 [13] 进行计算和数据处理。浮游动物生物量计算参照《水生生物学》 [7] 进行。

2.2.4. 浮游生物数据处理

实验数据采用Excel软件进行处理。浮游生物多样性指数采用Shannon-Wiener指数和均匀度指数，优势种利用优势度进行确定。加权综合营养状态指数(TLI)参照《养殖水域生态学》 [14] 进行，评价标准为：70 < TLI ≤ 100为极富营养型，重度污染；60 < TLI ≤ 70为(中度)富营养型，中度污染；50 < TLI ≤ 60为(轻度)富营养型，轻度污染；0 < TLI ≤ 30为中营养型，水质良好；0 < TLI ≤ 30为贫营养型，水质为优。

3. 结果与分析

3.1. 理化指标

牛口峪水库理化指标测定和富营养化状态如表1和表2所示。该水库水型总体为氯化物类钠组I型水(ClNa I)，水型有所变化，2017年6月和2021年5月水型为CNa I，2017年9月水型CMg I，2018年6月和9月水型为ClNa I。牛口峪水库为碱性水体，透明度为0.4~1.5 m，水温为19.7℃~29.2℃。参照《地表水环境质量标准》(GB3823-2002)可知，牛口峪水库的化学需氧量和溶氧总体均值属于国家I类标准，但在2017年秋季溶氧为II类标准；总氮除2018年夏季为III类标准，其余皆为V类标准及以上；总磷平均为III类标准，2017年夏季和秋季、2018年夏季和秋季、2021夏季分别为III、I、IV、III、IV类标准。综合2018~2021年五次调查，加权综合营养状态指数(TLI)显示，牛口峪水库为轻度富营养型。

表1. 牛口峪水库水体的理化特征

表2. 牛口峪水库理化特征及健康评价

3.2. 浮游植物

3.2.1. 浮游植物种类组成

调查期间牛口峪水库浮游植物共发现112种，其中，硅藻最多为45种，占总种数的40.18%；绿藻次之为33种，占比29.46%；蓝藻门为22种，占比19.64%；甲藻门4种，占比3.57%，金藻门和隐藻门各3种，占比分别为2.68%，黄藻门和裸藻门各1种，占比分别为0.89%。发现高耐盐种类：钝顶节旋藻(Arthrospira platensis)、翼茧形藻(Amphiprora alata)、波罗的海布纹藻(Gyrosigma balticum)、洛伦菱形藻(Nitzschia lorenziana)、牟氏角毛藻(Chaetoceros muelleri)、小三毛金藻(Prymnesium parvum)。

3.2.2. 浮游植物密度

牛口峪水库浮游植物平均密度为5.35 × 10⁶ind./L，其中硅藻最多，占比为37.78%；绿藻次之，为31.08%；第三为蓝藻，占比22.62%；黄藻占比5.79%；隐藻占比0.95%；金藻占比0.93%；甲藻占比0.48%；裸藻占比0.37%。从图2可见，浮游植物密度显示出季节和年度变化，2017年夏季蓝藻密度最多，其余时间硅藻密度占优，同年秋季绿藻密度增大。

图2. 牛口峪水库浮游植物密度及其组成的时间变化

3.2.3. 浮游植物生物量

牛口峪水库浮游植物平均生物量为7.85 mg/L，其中硅藻最多，占比70.85%；绿藻次之，占比13.57%；金藻第三，占比4.74%；黄藻占比3.95%；蓝藻占比3.35%；甲藻占比1.55%；隐藻占比1.23%；裸藻占比0.76%。从图3可见，浮游植物生物量显示出季节和年度变化，浮游植物生物量秋季大于夏季。组成上除2018年秋季绿藻生物量较大外，均以硅藻占优。显示该水库浮游植物为硅藻–绿藻型。

图3. 牛口峪水库浮游植物生物量及其组成的时间变化

3.2.4. 浮游植物多样性指数与均匀度指数的变化

牛口峪水库浮游植物平均多样性指数为3.60；平均均匀度指数为0.53，整体来看，牛口峪水库浮游植物多样性介于较好与非常丰富之间。多样性指数和均匀度指数显示出季节和年度变化，夏季大于秋季(图4)。

图4. 牛口峪水库浮游植物多样性指数与均匀度指数的时间变化

3.2.5. 浮游植物优势种及优势度

牛口峪水库浮游植物优势种共有22种，优势种以绿藻和硅藻为主，主要优势种由尖脆杆藻(Fragilaria acus)、疏刺多芒藻(Golenkinia paucispina)、梅尼小环藻(Cyclotella meneghiniana)和浮球藻(Planktosphaeria gelatinosa)等。2021年夏季，优势种数量最多，其余年度和季节优势种较少。

3.3. 浮游动物

3.3.1. 浮游动物组成及其分布

浮游动物共发现111种，其中轮虫最多为48种，占总种数的43.24%；原生动物次之为41种，占比36.94%；桡足类为10种，占比9.01%，枝角类和兼性浮游动物各6种，分别占比5.41%。发现高耐盐种类：褶皱臂尾轮虫(Brachionus plicatilis)、变异臂尾轮虫(Brachionus variabilis)、哈瓦那臂尾轮虫(Brachionus havannaensis)，哈瓦那臂尾轮虫为北京首次发现，是中国轮虫的新记录种。

3.3.2. 浮游动物密度

牛口峪水库浮游动物平均密度为9430 ind./L，其中原生动物最多，占比为55.09%；轮虫次之，占比44.08%；桡足类占比0.4；枝角类占比0.35%；兼性浮游动物占比0.4%。从图5可见，牛口峪水库浮游动物密度有明显的季节和年度变化，2021年夏季密度最大，2019年秋季最小，密度组成随季节变化，原生动物和轮虫交替占优势。

图5. 牛口峪水库浮游动物密度及其组成的时间变化

3.3.3. 浮游动物生物量

牛口峪水库浮游动物生物量平均值为2.61 mg/L (图6)，其中轮虫最多，占比为53.44%；原生动物次之，占比为28.64%；桡足类占比为9.42%；枝角类占比3.68%；兼性浮游动物占比4.83%。浮游动物生物量及其组成有明显的季节和年间变化，2021年夏季生物量最大达4.48 mg/L，原生动物占优；2018年秋季次之，为3.87 mg/L，轮虫占优；2018年秋季最小，为1.01 mg/L，枝角类占优。总体看夏季生物量高于秋季。

图6. 牛口峪水库浮游动物生物量及其组成的时间变化

3.3.4. 浮游动物多样性指数与均匀度指数的水平变化

浮游动物多样性指数为3.42，均匀度指数为0.64，多样性指数和均匀度指数显示出季节和年度变化，2017年夏季两种指数优于秋季，而2018年秋季优于夏季，2017和2018年均优于2021年(图7)。

图7. 牛口峪水库浮游动物多样性指数与均匀度指数的时间变化

3.3.5. 浮游动物优势种及优势度

牛口峪水库浮游动物优势种共有30种，优势种由原生动物和轮虫组成，主要的优势种为小单环栉毛虫(Didinium balbianii nanum)、旋回侠盗虫(Strobilidium gyrans)、萼花臂尾轮虫(Brachionus calyciflorus)、褶皱臂尾轮虫、环顶巨腕轮虫(Hexathra fennica)、针簇多肢轮虫(Polyarthra trigla)。值得注意的是，本地调查中出现高耐盐种类如褶皱臂尾轮虫。

图8. 牛口峪水库生物指标与环境因子的CCA分析

3.4. 牛口峪水库生物与环境关系的CCA分析

牛口峪水库生物指标与环境因子的典范对应分析结果如图8所示。可见，浮游植物密度(PD)、生物量(PB)、多样性指数(PH)、均匀性指数(PJ)、浮游动物密度(ZD)、均匀性指数(ZJ)和多样性指数(ZH)与水温(WT)、碱度(TA,ALK)、硬度(TH)呈正相关；与pH、叶绿素a (Chl-a)、总氮(TN)、溶解氧(DO)、化学需氧(COD)量呈负相关。

4. 讨论

4.1. 合理性分析

本研究选择水体理化指标和水生生物群落特征作为研究对象。水体理化指标参照《地表水环境质量标准》(GB3823-2002)，指标越大则说明水体健康状况越差；综合五次采样的水体理化特征，TLI值越大说明水体营养程度越高，水体污染程度也越高；因此用水体理化指标和TLI来评价水体健康是较为合理的。水生生物密度反映了水体中浮游生物群落的多少，密度越大说明浮游生物群落在水体中越多，水生生物群落稳定性越高；生物量反映了水体中浮游生物的总个数或总干重，生物量越高，说明水体营养程度越高，水体食物网结构越稳定；多样性指数越大，说明水体中的物种越丰富，食物网结构越复杂，水体生态系统越稳定，水体抵御外部变化的能力越强；均匀度指数反映了各物种个体数目分配的均匀程度，均匀度指数越接近1，则各物种间的个体数越接近，生态系统就越稳定；因此水生生物的密度、生物量、多样性指数、均匀度指数可以反映水体的健康程度。

4.2. 牛口峪水库水质安全评价

牛口峪水库氮磷处在较高水平，其中总氮含量偏高，总磷升高，总氮和总磷增加原因可能与地表径流有关 [15] 。水体理化因子中，pH和碱度显示水体具有一定的盐碱化，这可能与牛口峪水库所属的燕山石化公司所处理的水的制剂中含有一定盐度有关，盐度升高可导致细菌、浮游生物等种类和丰度的减少，利用水中营养物质能力减少，可能会导致水体富营养化。TLI表明牛口峪为轻度富营养状态，污染程度总体处于轻度污染，但不容忽视的是从2017到2018年，再到2021年水质营养状态不断加深，其中水质中的总氮总磷也可以侧面反映一定的水质变化趋势，可能存在爆发水华的危险，针对蓝藻水华暴发的氮磷阈值，Gilbson等指出湖泊和水库中总氮、总磷浓度分别超过0.15和0.01 mg/L即有可能发生蓝藻水华，与本次研究结果对比，表明有爆发水华的可能 [16] 。

4.3. 牛口峪水库浮游生物群落变化特征

该水库浮游生物中有很多耐盐种类，如浮游植物钝顶节旋藻、翼茧形藻、波罗的海布纹藻、洛伦菱形藻、牟氏角毛藻、小三毛金藻，浮游动物褶皱臂尾轮虫、变异臂尾轮虫，浮游生物耐盐种类与水中高pH和碱度有关，这也进一步印证了水中含有一定盐度的猜想。本次调查的牛口峪水库浮游植物和浮游动物数量和丰度与北京市其他水库相一致 [3] [4] [5] 。

2021年浮游生物密度和生物量与2017年、2018年相比，浮游植物和浮游动物的密度和生物量均有所下降，浮游植物呈现密度小，生物量大趋势，而浮游动物则呈现密度大生物量小，这可能与季节和地理分布差异有关。浮游植物优势种以硅藻门和绿藻门为主，且大多数为优质藻类，这从侧面反映牛口峪水库具有一些优良饵料基础。浮游植物的群落组成也能反映水体所处的营养状态 [17] ，营养化加深水体中，绿藻门和硅藻门种类比例很大 [18] ，本次研究硅藻占比45%，绿藻占比33%，与上述得到的水质营养状态加深相呼应。浮游动物中原生动物和轮虫为优势种，这从侧面说明存在着一定的小型化趋势，小型化趋势可能是由于浮游动物面临较大的捕食压力所致，而浮游动物的小型化将会导致浮游动物对浮游植物的控制力减弱，同时潜在富营养化风险在增加 [19] 。牛口峪水库浮游动物平均多样性指数(H' = 3.42 > 3)和均匀度指数(J = 0.64, 0.5 < J < 0.8)结果均显示为轻污染 [20] ，这与水质综合营养状态为轻度富营养一致，但仍要重视潜在的浮游植物水华可能和浮游动物小型化。

4.4. 牛口峪水库水质和水生生物多样性可持续利用对策

首先，建议减少使用含盐量多的水处理制剂，这样可以减少水库的盐碱化程度，增加水生水生物种类，提高水体中氮磷的利用效率，减轻水体富营养化程度。其次，应该进一步控制氮磷，切断外界污染源，减少污染物进入水体的量，如周边的山体土壤中大量营养盐通过地表径流进入水体，造成水中营养物质的增加，容易使水体富营养化，所以应改善周边的生态环境，加大植树造林，减少土壤径流对水体的污染。

5. 结论

综上所述，牛口峪水库水质和生物多样性仍需改善，鉴于该水库是生态环境恢复与保护的有机结合的范例，应在控制污染源的基础上，大力进行优水生态工程和生物监测。

基金项目

国家自然科学基金项目(42077226)和北京市农业农村局项目(2017-2021)资助。

参考文献