我国典型湖泊磷浓度对水体指标及蓝藻水华爆发趋势的潜在影响——基于Meta分析

王金航 李京霖 江德霜 武会会 王曙光 林爱军*

(1.北京化工大学 北京市环境污染控制与资源化工程研究中心， 北京 100029；2.北京市可持续发展科技促进中心， 北京 100035)

引 言

湖泊水体富营养化是我国面临的主要水环境问题[1]。我国第二次湖泊现状调查显示，85.4%的湖泊达到富营养化标准，40.1%的湖泊达到严重富营养化标准[2]。随着水体富营养化日益严重，湖泊水华现象频繁发生，尤其在我国的滇池、太湖和巢湖最为严重。水华是由藻类、浮游动植物和细菌等爆发性繁殖或高度聚集而引起水体变色的现象[2]。在各类水华中，蓝藻水华在淡水中占据优势，被公认是全球最严重的湖泊水环境问题之一。蓝藻水华的爆发会导致水体缺氧、水质变坏，同时蓝藻还会产生蓝藻毒素，直接威胁湖泊水环境和人类健康[3]。

蓝藻水华爆发的主要原因是水体富营养化，同时易受气象、水文、人类活动等因素的影响。Dokulil等[4]认为营养物质浓度、水温、光照等条件是影响蓝藻优势形成的关键。孔繁翔等[5]提出了蓝藻水华成因的“四阶段理论”，探索了符合我国湖泊特征的蓝藻水华爆发的规律和机制。徐潇等[6]从水- 土- 气- 生多角度研究了各环境因素对蓝藻水华爆发机制的影响。目前，我国对蓝藻水华爆发的研究多集中在单一研究区域，缺乏多维度的蓝藻水华爆发规律的研究。

近年来，Meta分析作为一种统计方法，可通过汇总初步研究中的数据，得出比单一研究更可靠的综合分析结果，避免单一研究的局限性[7]。该方法已广泛应用于医学领域，并逐步应用于土壤、水、大气等环境科学领域[8-10]。徐晓晔等[11]通过Meta分析函数构建了我国水资源价值移位模型，揭示了我国地理空间格局的水资源价值的差异性。Khan等[12]通过Meta分析研究了印度河流域的气候变化与流域水文的关系，为印度水供应不足制定策略。Erni-Cassola等[13]利用Meta分析揭示了海洋不同采样区域中常见塑料聚合物相对丰度的总体趋势，为海洋微塑料治理提供了依据。这些研究表明Meta分析可全面揭示水环境流域特征，为水环境保护和管理策略的制定提供支撑，但目前没有使用Meta分析对流域环境水华的爆发规律进行研究的报道。Conley等[14]提出氮磷双控是防治湖泊富营养化和蓝藻水华爆发的关键，但磷是最主要的限制性养分。因此，本文从1989—2019年公开发表的相关研究中提取583对有效数据，利用Meta分析，以我国典型湖泊(太湖、巢湖、滇池、白洋淀和松花湖)为对象，研究了磷浓度对水质指标以及蓝藻水华爆发的影响，为蓝藻水华爆发的防治提供参考，对水环境保护具有一定的现实意义。

1 研究方法

1.1 文献收集

利用中国知网(CNKI)、万方和维普等国内期刊数据库，以“蓝藻爆发”、“富营养化”、“磷”和“湖泊”为关键词，检索1989—2019年发表的有关“我国湖泊水体磷浓度与蓝藻水华爆发”的文献，对其进行Meta分析。涉及湖泊包括太湖、巢湖、滇池、白洋淀和松花湖。为达到本研究的目的和减少筛选文献带来的偏差，制定以下6条文献筛选标准：1)在我国湖泊进行的实验性或监测性研究须设有实验组和对照组，并以最低磷浓度作为对照组；2)在研究过程中，同一项研究须控制水体环境、气象条件、实验时长一致；3)所测点位至少观测一个月以上，并取其均值，避免短期及气候条件等导致实验误差；4)至少包含一项以下参数：总磷(TP)含量、溶解氧(DO)含量、水体透明度(SD)、高锰酸盐指数(CODMn)和叶绿素a(Chla)含量；5)对于报道同一个研究的重复文献，只选择其中一篇；6)相关指标参数须以均值表示，并包含标准差(标准误或置信区间)和样本容量，必要时可从图表中提取相关数据。基于上述筛选标准，由最初的5 000余篇文献，最终获得15篇有效文献(见表1)，并利用拉依达准则(3σ准则)剔除异常值后获得583对有效数据。同时，数据的预处理结果显示，失安全数为4.57×1011(P<0.1)，大于Hedges等[30]提出的参考限值5×1010，即文献的发表偏倚很小，基于此数据的Meta分析结果稳定。

表1 研究的4项指标及其筛选文献Table 1 Four research indexes and their screened literatures

1.2 数据整理

对每篇文献中湖泊的地理位置、气候条件、TP含量、DO含量、SD、CODMn、Chla含量的观测数据进行记录，同时也包括对照组和实验组的均值、标准差、样本容量。对表格形式的数据直接进行提取，对图片形式的数据使用Origin 2018软件(OriginLab，Northampton，MA，USA)将其数字化。

为研究不同磷浓度下湖泊水质指标的变化及蓝藻水华爆发的趋势，将磷的质量浓度分为4类：0～0.1、0.1～0.2、0.2～0.3 mg/L及>0.3 mg/L。此外，有研究表明气候因素对蓝藻水华爆发有较大影响[4-5]。因此，根据我国蓝藻水华爆发区域的主要气候特征，将亚热带季风气候地区(年均气温13～20 ℃)和温带季风气候地区(年均气温11.7 ℃)设为亚组进行Meta分析。

1.3 数据分析

采用MetaWin 2.0软件(Sinauer Associates Inc，Sunderland，MA，USA)中的随机效应模型计算合并效应值及其95%置信区间。Meta分析的效应值可分析某一指标的实验组相对于对照组的变化幅度，以反应比(r)的自然对数(lnr)为效应值(E)，其计算公式如下。

E=lnr=ln (XE/XC)=lnXE-lnXC

(1)

式中，XE和XC分别表示实验组和对照组中研究指标的均值。若E=0，则说明磷浓度未引起实验组与对照组之间的差异；若E>0，则说明磷浓度对所研究的指标参数产生正效应：若E<0，则说明磷浓度对所研究的指标参数产生负效应。

效应值方差(v)的计算公式如下。

(2)

式中，SE和SC分别表示实验组和对照组中研究指标的标准差；NE和NC分别表示实验组和对照组的样本容量。为了得出实验组相对于对照组的效应，对效应值均值进行加权处理，分析效应值及其95%置信区间(95% CI)。若效应值的95% CI与效应值的“0”点相交，则认为磷浓度对该指标的影响不显著(P>0.05)；反之，影响效果显著(P<0.05)。

湖泊水体中磷浓度变化后，水体各项指标的变化率通过以下公式[31]获得。

D=(eE-1)×100%=(eln r-1)×100%

(3)

式中，D为实验组相对于对照组的指标变化率。

采用Origin 2018软件进行相关图件的制作，使用SPSS 25.0软件(IBM，Armonk，NY，USA)进行Pearson相关性分析，P<0.05为相关性显著，P<0.01为相关性极显著。

2 结果与讨论

2.1 湖泊中水质指标对磷浓度的总体响应

湖泊中DO含量、SD、CODMn和Chla含量对磷浓度的总体响应如图1所示。Meta分析结果表明，湖泊中磷浓度对SD产生负效应，其变化率为-38.57%；对CODMn和Chla含量产生正效应，其变化率分别为17.35%和97.76%。DO含量效应值的95%置信区间与效应值的“0”点相交，说明水体中磷浓度对其影响不显著(P>0.05)。但是，水体中磷浓度对SD、CODMn和Chla含量的影响显著(P<0.05)。

括号内的第一个数字表示样本容量，第二个数字表示指标变化率。图1 湖泊中DO含量、SD、CODMn和Chla含量对磷浓度的总体响应Fig.1 The overall response of DO content, SD, CODMn and Chla content to phosphorus concentration in lakes

太湖、巢湖、滇池、白洋淀和松花湖5个不同地区的湖泊中，水体中DO含量、SD、CODMn和Chla含量对磷浓度的响应如表2所示，从表2中可以看出5个湖泊中DO含量对磷浓度的响应均不显著(P>0.05)。在上述湖泊(不包含滇池)中磷浓度对SD均产生负效应，且效果显著(P<0.05)，说明磷浓度升高会显著降低水体透明度。5个湖泊中仅太湖中磷浓度对CODMn的效应显著，产生正效应，而对其他湖泊的效果均不显著。这可能是由于对太湖的相关研究主要集中在2007—2011年，正值2007年太湖蓝藻污染事件，湖泊整体水质较差，因此CODMn效应显著高于其他湖泊[32]。在Chla含量对磷浓度的响应方面，磷浓度对湖泊中Chla含量产生显著的正效应，说明磷浓度升高是导致蓝藻水华爆发的关键因素[2]。

表2 我国典型湖泊中DO含量、SD、CODMn和Chla含量对磷浓度的响应Table 2 Response of DO content, SD, CODMn and Chla content to phosphorus concentration of typical lakes in China

2.2 湖泊中磷浓度对水体指标的影响

湖泊中磷浓度是影响水体富营养化程度及蓝藻水华爆发的主要因素，此外，温度、光照时间、水利条件、风速和风向等自然因素也是影响湖泊蓝藻水华爆发的关键因素[2]。太湖、巢湖和滇池属于亚热带季风气候地区，白洋淀和松花湖属于温带季风气候地区，故在分析湖泊中磷浓度变化的同时，也根据湖泊的气候特点对各指标进行对比分析，以明确不同气候地区的湖泊中蓝藻水华爆发的规律。

2.2.1DO含量

DO含量反映水体的污染程度，是影响水环境健康的重要因素[33]。湖泊中磷浓度对DO含量的影响如图2所示。可以看出，在2种气候条件下当磷的质量浓度为0～0.1 mg/L时磷浓度对DO含量产生正效应，这可能是由于适量的磷可为水生植物提供养分，促进光合作用，从而使DO含量增加。随着磷浓度的升高，磷浓度对DO含量产生负效应。在磷的质量浓度为0.2～0.3 mg/L(亚热带季风气候)和>0.3 mg/L(温带季风气候)时负效应最大，DO含量变化率分别为-12.74%和-27.32%。浮游植物会通过光合作用产生氧气，从而增加水体中DO含量，但其呼吸作用又会消耗氧气[33]。随着磷含量逐渐升高，浮游动物和藻类迅速繁殖，占据大量湖泊表层面积，下层浮游动物和藻类进行光合作用的强度减弱。此外藻类死亡分解也会消耗氧气，当光合作用的产氧量小于呼吸作用和藻类分解的耗氧量时，水体DO含量下降，进而危害水生生物健康和水生态系统安全[2]。

括号内的第一个数字表示样本容量，第二个数字表示DO含量变化率。图2 磷浓度对DO含量的影响Fig.2 Effect of phosphorus concentration on DO content

在自然条件下，湖泊水体中的DO含量会处于一个相对稳定的范围，水温是影响其变化的主要因素[34]。温度的改变可以在很大程度上影响有机物的产生和积累，同时也会影响水体的氧动态和氮磷循环过程[35]。从图2可以看出亚热带和温带季风气候地区的湖泊中DO含量随磷浓度变化的趋势相似。在磷的质量浓度为0～0.1 mg/L时，磷浓度对亚热带季风气候地区DO含量的正效应显著(P<0.05)，而对温带季风气候地区的正效应不显著(P>0.05)。这可能是由于亚热带季风气候地区的光照强度和水温较高，有利于水生植物的光合作用，从而导致水体DO含量升高[33]。

2.2.2SD

SD是反映水体浑浊程度的重要参数，在研究水环境变化和水体富营养化等方面都具有重要意义[36]。磷浓度对SD的影响如图3所示，磷在不同浓度下均会对SD产生显著负效应(P<0.05)，且随磷浓度的升高负效应先增大后减小。在2种气候条件下，当磷的质量浓度为0.1～0.2 mg/L时磷浓度对SD的负效应最明显，SD变化率分别为-42.59%(亚热带季风气候)和-57.33%(温带季风气候)。Chang等[36]构建了水体透明度评价体系，指出SD受悬浮颗粒物和Chla的影响最大。磷浓度的增加会导致藻类加速生长，产生大量Chla，降低水体透明度。但藻类植物的生长在一定程度上可以起到净化水体的作用，Brune等[37]提出菌藻共生体系可通过藻类光合作用和好氧菌代谢，对水体进行净化，同时藻类和菌株还会吸附水中的悬浮颗粒物，改善水体透明度。因此，在磷浓度增加时，负效应不会持续降低，而是呈现轻微的反弹趋势。

括号内的第一个数字表示样本容量，第二个数字表示SD变化率。图3 磷浓度对SD的影响Fig.3 Effect of phosphorus concentration on SD

亚热带和温带季风气候地区的湖泊所呈现的规律相似。在磷的质量浓度为0～0.1 mg/L和0.1～0.2 mg/L时，温带季风气候地区的SD负效应明显高于亚热带季风气候地区。这样的差异主要是由温度、光照、降水量不同导致的，亚热带季风气候地区的年均温度较高，更利于浮游动物和藻类的繁殖，对水体的净化效果更好，所以负效应较低。

2.2.3CODMn

CODMn是我国评价地表水质的一项重要指标，能在一定程度上反映地表水质量[38]。磷浓度对CODMn的影响如图4所示，磷在不同浓度下对CODMn均产生正效应，在磷的质量浓度为0.2～0.3 mg/L(亚热带季风气候)时CODMn变化率最大，为74.46%。CODMn增加是水体中有机物含量增加导致的，藻类和浮游动物的大量繁殖会使CODMn升高[39]。韩志萍等[40]发现湖泊中蓝藻生物量与CODMn呈正相关，且高峰期一致。其他相关研究[40-41]均发现CODMn对蓝藻水华爆发有显著影响，因此CODMn的管控是控制蓝藻水华爆发的重要措施之一。

括号内的第一个数字表示样本容量，第二个数字表示CODMn变化率。图4 磷浓度对CODMn的影响Fig.4 Effect of phosphorus concentration on CODMn

2.2.4Chla含量

Chla含量作为浮游植物存量的表征指标[42]，可有效反映水体富营养化情况和蓝藻水华爆发趋势。例如，Aubriot等[43]利用遥感技术监测Chla浓度，建立了蓝藻水华爆发预警系统，为普拉塔河的蓝藻水华爆发提供了阻控依据。磷浓度对Chla含量的影响如图5所示，Chla含量对磷浓度的响应趋势与CODMn基本相似。在亚热带季风气候条件下，随磷浓度的增加正效应逐渐增大，在磷的质量浓度>0.3 mg/L时效应值最大，Chla含量变化率为326.99%。在温带季风气候条件下，在磷的质量浓度为0.1～0.2 mg/L时正效应最大，Chla含量变化率为154.77%。

括号内的第一个数字表示样本容量，第二个数字表示Chla含量变化率。图5 磷浓度对Chla含量的影响Fig.5 Effect of phosphorus concentration on Chla content

在不同气候地区，磷浓度对Chla含量均产生正效应，在磷的质量浓度为0～0.1 mg/L时，亚热带季风气候地区的Chla含量变化率(86.60%)高于温带季风气候地区(67.21%)，说明此时其发生水体富营养化的风险相对较高。不同气候地区的水温差异较大，水温可直接影响浮游植物体内酶的活性，调控浮游植物的光合作用和呼吸代谢速率，进而影响湖泊中Chla的浓度[43]。

2.3 相关性分析

为了论证Meta分析的结果以及磷浓度与水体各指标之间的关系，对各指标间的相关性进行了分析，结果如表3所示。DO含量、SD、CODMn与TP含量呈现显著或极显著相关性，其中DO含量和SD与TP含量呈负相关，CODMn与TP含量呈正相关。水体中的营养物质是蓝藻生长繁殖的物质基础，营养物质含量决定着蓝藻的繁殖速度，对水华的发生有至关重要的影响[2]。结果表明TP含量是影响湖泊水质指标的重要因素，限制湖泊中磷浓度是防治水体富营养化和维护湖泊水质的重要手段。

表3 我国典型湖泊的水质指标相关系数Table 3 Correlation coefficients of water quality indexes of typical lakes in China

但是相关性分析结果显示TP含量与Chla含量的相关性不显著，这与前文Meta分析中TP含量会显著影响Chla含量的结论存在出入，其原因分析如下：1)本文选择5个不同地区的湖泊进行相关性分析，不同湖泊间TP含量和Chla含量的关系可能受自然因素(水利条件、气象条件等)和人为因素(实验差异等)的影响，其间的相关性被破坏；在Meta分析中由于选用同一篇文献的实验组与对照组进行分析，可以排除一定的人为因素和自然因素的影响。2)研究变量的差异。Meta分析为实验组与对照组间的对比分析，而相关性分析属于TP含量和Chla含量间关系的分析，故会产生一定的差异。

3 结论

(1)湖泊中磷浓度对SD产生负效应，使SD降低38.57%；对DO含量、CODMn和Chla含量产生正效应，使其分别增加了1.39%、17.35%和97.76%，其中对DO含量的影响不显著(P>0.05)。

(2)随着湖泊中磷浓度的增加，各指标呈现不同的变化趋势。磷浓度对DO含量的正效应逐渐转变为负效应，在磷的质量浓度为0.2～0.3 mg/L(亚热带季风气候)和>0.3 mg/L(温带季风气候)时对DO含量的负效应最大；对SD的负效应随磷浓度的升高先增大后略微减小，在磷的质量浓度为0.1～0.2 mg/L时负效应最大；CODMn和Chla含量对磷浓度的响应趋势基本相似，在亚热带季风气候条件下，磷的质量浓度分别在0.2～0.3 mg/L和>0.3 mg/L时正效应最大。

(3)在不同气候条件下，在磷浓度变化时湖泊的水质和蓝藻爆发情况不尽相同。在磷的质量浓度为0～0.1 mg/L时亚热带季风气候地区湖泊的DO含量、CODMn和Chla含量对磷浓度的响应程度高于温带季风气候地区，更易受磷浓度的影响，蓝藻水华爆发的风险相对较高。

在不同湖泊中，蓝藻水华的影响因素间交叉和组合的作用是不同的，因此，不同水体蓝藻水华爆发的原因和机制可能有所不同。在下一步的研究中建议采用网状Meta分析，进行多因素的混合效应研究。同时，在讨论特定水体的蓝藻水华状况时，有必要结合当地实际情况对其进行详细分析。