1、第 卷 第 期 年 月人 民 长 江 ，收稿日期：基金项目：国家自然科学基金项目“输水明渠寡营养水体沉积物 水界面营养盐活化机制与藻类协同演化”（）；长江生态环境保护修复联合研究二期项目（）作者简介：张 洪，男，研究员，博士，主要从事流域物质循环与管理调控研究工作。：文章编号：（）引用本文：张洪，薛雪，郁达伟，等 鄱阳湖水位对沉积物磷释放的影响及总磷考核建议 人民长江，（）：“水生态文明建设”专题 特邀作者简介张 洪，男，年生，研究员，博士，现任中国科学院生态环境研究中心流域物质循环与调控课题组组长，中国科学院大学资环学院水质学教研室副主任，国家生态环境保护专业技术青年拔尖人才，主要从事流域水
2、系磷迁移转化过程研究与修复工程实践。近年来，主持承担了国家自然科学基金、国家重点研发课题、中国科学院重点部署课题、国家水体污染控制与治理科技重大专项（子课题）余项，发表论文 余篇，申请 获授权专利 余项。担任环境工程技术学报青年编委，环境工程学报栏目特邀学术编辑，东南大学、济南大学校外兼职导师等。鄱阳湖水位对沉积物磷释放的影响及总磷考核建议张洪，薛雪，郁 达 伟，魏 源 送，桂 双 林，魏 志 华（中国科学院生态环境研究中心 环境水质学国家重点实验室，北京；中国科学院大学，北京；江西省科学院 能源研究所，江西 南昌；南昌市新建区环境保护监测站，江西 南昌）摘要：为探讨干湿交替对湖泊沉积物内源磷
3、释放的影响，结合基础资料和模拟试验，研究了鄱阳湖水位波动及其驱动下沉积物中磷的释放过程，并对湖区总磷（）考核提出建议。结果表明：年鄱阳湖星子站水位高幅波动，最高水位差分别为 ，和 ，导致水面面积也呈现高度动态性，最低水位时的裸露面积分别达到 ，和 ，大面积沉积物会经历干湿交替过程；沉积物磷释放模拟试验结果显示，南矶山和青岚湖断面沉积物经历裸露 暴晒 淹水后，淹水初期沉积物 的平均释放通量分别高达 （）和 （）；综合考虑鄱阳湖水体中 组成特征和水位波动的影响，建议湖区国控断面以活性磷酸盐（）或者溶解性总磷（）替代 作为磷考核指标，并动态调整磷考核标准的阈值，当星子站水位低于 时以河流标准考核。关
4、 键 词：沉积物磷释放；总磷考核；水位变化；鄱阳湖中图法分类号：文献标志码：第 期 张 洪，等：鄱阳湖水位对沉积物磷释放的影响及总磷考核建议 引 言长江流域磷来源复杂、“三磷”问题突出，已对流域水体中磷浓度产生了重要影响，而磷浓度又是评价河湖水生态系统健康与否的重要指标。近年来，长江流域水质总体呈现好转趋势，但总磷已成为流域水体首要污染因子。其中，年，长江流域总磷作为首要超标因子的断面占比为，高于氨氮占比，也高于高锰酸盐指数占比 。作为长江流域最大通江湖泊的鄱阳湖，其水体磷污染问题也十分凸显，总磷普遍超标且浓度逐年增加，成为制约鄱阳湖流域（江西省）经济社会可持续发展的重要因素。随着长江大保护工
5、作的推进，外源污染得到有效控制后，内源磷释放将成为影响鄱阳湖水体总磷浓度的关键因素。磷很难通过气态形式从水体中脱除，因此，水体沉积物是流域磷的主要蓄积库，同时其内源磷释放也会对水环境产生重大影响。水体沉积物中磷可通过两种途径释放到上覆水中：自由扩散过程，将附着在沉积物颗粒物上的磷先释放到间隙水中，再通过浓度梯度或外力扰动，将磷扩散到上覆水中；再悬浮过程，在外界自然作用下，沉积物中含或吸附磷的颗粒物质再悬浮到上覆水中，可增大上覆水中磷含量。沉积物中磷释放过程受到环境因素的显著影响，包括沉积物 上覆水性质、水体水量变化、扰动作用等，。目前，鄱阳湖沉积物中磷释放影响因素研究已涵盖、溶解氧、温度等，但
6、有关水位波动引起沉积物干燥和再淹没过程对磷释放的研究未见报道。长江干流上游水库群运行对鄱阳湖与长江之间水沙交换过程产生了重要影响，推动江湖关系演变，加剧了鄱阳湖区秋旱现象频繁发生，以及鄱阳湖的水位波动。加之 年鄱阳湖流域降雨量较同期偏少，到 月 日：，鄱阳湖星子站水位跌破了 ，标志着鄱阳湖正式进入极枯水位。湖区通江水体面积仅，相较丰水期面积（），萎缩了 。到 月 日，星子站水位退至 ，刷新 年有记录以来历史最低水位。水位下降导致鄱阳湖湖盆形成大面积裸露，沉积物在高温暴晒下迅速干燥。随着后期水位回升，沉积物干湿交替过程将会对内源磷循环产生显著影响。因此，本文选择鄱阳湖南昌湖区典型断面，采集不同水
7、位沉积柱样品，开展干湿交替对沉积物中磷释放过程影响的模拟研究，为探明鄱阳湖沉积物中磷的生物地球化学循环规律提供依据，同时为实现鄱阳湖上覆水总磷考核达标提供支撑。研究区域和方法 研究区域鄱阳湖位于江西省境内，近 最大水面面积 （湖口水文站水位 ），属吞吐型湖泊，是中国第一大淡水湖，和洞庭湖有着“长江之肾”的美称。鄱阳湖湖区北部与长江相接，位于江西省北部，长江中下游南岸。以鄱阳湖星子站为界，北湖湖道狭长，宽度仅有 。鄱阳湖是一个过水性、吞吐型、季节性湖泊，它承接容纳赣江、抚河、信江、饶河、修河五大河流入湖，影响着长江径流流量和流域生态环境质量。鄱阳湖湖区容水量受季节性降水影响很大，具有“洪期一片是
8、湖、枯期一线似河”的独特景观，洪水期湖泊面积有时可达枯水期湖泊面积的几十倍。本研究选取星子水文站、南矶山断面和青岚湖断面为实验点，如图 所示。南矶山断面（，）位于南昌市新建区南矶乡，断面地处赣江中支和南支汇入鄱阳湖开放水域冲积形成的赣江三角洲前缘，是赣江两大支流的河口与鄱阳湖大水体之间的水陆过渡地带。青岚湖断面（，）位于南昌市进贤县英山圩以北的外青岚湖水域，与城区紧密相连，是鄱阳湖众多卫星湖之一。星子站是鄱阳湖具有标志性的水文监测站，一直以来记录着鄱阳湖不同时期的水位高程，因其位于五大河流汇合必经卡口而被设为水文监测站。图 研究区域采样点 沉积物磷释放潜力模拟和分析方法采用室内静态模拟方法估算
9、沉积物磷释放潜力。年 月，采集鄱阳湖南昌湖区柱状沉积物，连同上覆水一起运回实验室。利用高 ，内径为 的亚克力圆柱形实验管，在鄱阳湖流域的青岚湖和南 人 民 长 江 年矶山 个国控断面各采集 根沉积柱样品和 湖水。根沉积柱中 根作为实验对照柱，其余 根为平行实验柱。实验对照柱保证了实验点沉积物原位状态，未做任何处理，一直处于水淹状态（上覆水）。而 个平行样为模拟沉积物干湿交替过程对沉积物磷释放的影响，运回实验室静置 后，用注射器取走沉积柱全部上覆水（见图）。为保证实验的准确性，所有沉积柱底泥控制为 厚度。图 模拟装置 将抽干上覆水的沉积柱（个平行样）放置于室外曝晒直至底泥晒干，使沉积柱表层沉积物
10、处于干旱状态。平行实验组曝晒完成后加入对应采样点湖水至，和对照实验柱一起放置实验室进行常温静态培养，分别在第，天后，用注射器采集 水样装入离心管后，放入冰箱保存待测。水样取完后沿管壁缓缓注入相同体积的超纯水作为补充。所取水样用 滤膜过滤，用于活性磷酸盐（）和总溶解性磷（）测定，总磷（）测定水样无需过滤。水样中、测定参考 水质 总磷的测定 钼酸铵分光光度法。其中，和 测定时需对水样进行消解，参考 水质 金属总量的消解 微波消解法。为准确表达干湿交替对沉积物磷释放的影响，个平行样测样结果需扣除对照样测样浓度。数据处理方法沉积物静态释放通量 的计算方法见公式（）：（）（）（）（）单位面积释放量 计算
11、方法见公式（）：（）式中：为释放通量，（）；为单位面积沉积物的磷释放量，；为柱中上覆水体积，；，为初次、第 次和 次采样时某物质浓度，；为添加水样中的物质浓度，；为第 次采样时采样体积，；为沉积柱样中沉积物 水界面接触面积，；为释放时间，。水位 水质数据决策树模型分析方法采用 年鄱阳湖星子站水位数据和南矶山国控断面水质数据，基于决策树（）模型对该阈值进行了建模分析。建模主要采用 语言 包和 包，将数据集划分为训练集（）和验证集（），在条件推断的框架下建立二叉树回归推断，检验水位和水质类别之间统计无关的零假设直到假设不能拒绝，基于 值选择水位对水质最强的影响值作为一个节点。递归上述步骤，直至所有
12、无关假设被拒绝，则决策树完成建立。从而建立的每个二叉树根节点，即为水质类别突变点对应的水位阈值。所建立的模型对水质类别的预测正确率平均为 。统计和作图方法相关统计计算使用 软件，相关性分析采用 软件，图形处理使用 软件。结果与讨论 鄱阳湖多年水位波动过程分析根据 星子站逐月水位、年逐日水位资料，以及不同水位对应的水面面积，分析鄱阳湖水位 裸露面积的定量关系，结果如图 所示。由图（）可知：近 ，鄱阳湖水位呈现下降的趋势，由 年平均 降至 年 ，对应水面面积也呈现明显的减少趋势。鄱阳湖水位随季节和年际时间尺度变化均非常大，尤其 年多次跌破历史最低水位，与长江流域枯水特征基本保持一致。近年来，鄱阳湖
13、洪旱频繁，对人们生活及生态环境造成很大影响。洪旱事件与水位密切相关，鄱阳湖水位年内变化特征十分明显，高水位基本出现在洪水期（月）。年，鄱阳湖星子站最高水位分别为 ，和 ；最低水位分别为，和 ，且低水位大幅提前（见图（）。水位的高幅波动，导致鄱阳湖水面面积也呈现出高度的动态性，进而呈现不同尺度的裸露范围。由图（）可知：鄱阳湖水位与其裸露面积呈现很好的线性相关。计算可得，近 最低水位时鄱阳湖裸露面积分别达到 ，和 ，而本研究模拟期 月份，裸露面积分别为 ，和 。巨大的水面面积差，使 第 期 张 洪，等：鄱阳湖水位对沉积物磷释放的影响及总磷考核建议图 鄱阳湖多年水位波动（年）及其与湖盆裸露面积关系
14、得鄱阳湖湖区大量沉积物经历干湿交替过程，会对湖区磷的生物地球化学循环过程产生重要影响。水位驱动下鄱阳湖沉积物磷释放分析图 展现了干湿交替驱动下鄱阳湖沉积物中磷释放通量的变化过程。由图 可知：南矶山和青岚湖断面沉积物在经历裸露 暴晒 淹水后，沉积物中磷均出现不同程度释放过程。整体上，个断面的沉积物在模拟第 天出现大幅释放过程，随后迅速降至微释放状态，而对照样沉积物磷通量变化甚微。淹水初期南矶山断面沉积物、和 的平均释放通量分别为 ，（）和 （）；青岚湖断面分别为 ，（）和 （）。整个模拟期间，两个断面沉积物单位面积 释放量均值分别为 和 ，结合 年 月鄱阳湖平均裸露面积，估算可知整个湖盆裸露区水
15、位上升后，由于干湿交替导致的 释放量预计在 之间。已有研究估算，鄱阳湖最大支流赣江 年输入约 ，但其输入伴随大量河流泥沙颗粒，且流动水体水环境稀释容量较大。而湖体干湿交替诱导的沉积物磷释放，属于零容量的直接输入，会对局部区域上覆水 浓度产生高强度冲击，尤其是对国控断面 考核的影响不容忽视。水位变化是驱动湖泊沉积物中氮磷等营养盐迁移转化的重要因子，进而对湖泊上覆水水质产生重要影响。本研究通过模拟鄱阳湖沉积物干湿交替过程磷释放效益，发现 个断面退水后裸露的沉积物再次淹水，会诱导短时间集中释放过程。其中，磷释放组成以颗粒态为主，分别占比 和，占比较低，仅为 和。已有研究表明，在干湿交替过程中，沉积物
16、中磷赋存形态及其含量均会发生变化，增加沉积物中各形态磷的释放风险。而这一过程将直接影响到鄱阳湖国控断面 的考核达标。注：图中 条线代表 个平行样。图 鄱阳湖沉积物磷释放通量模拟结果 鄱阳湖水位与水体磷浓度关联性分析水位 水质关系的决策树模型分析结果表明（见图），水位在极大程度上决定了鄱阳湖南昌湖区断面（南矶山）的水质优劣。随着星子站水位从 经 人 民 长 江 年图 鄱阳湖（南矶山断面）水质类别天数占比与水位（星子站）的关系 （）（）上升到 以上，对应的鄱阳湖断面水质类别总体上从地表劣类经地表类、地表类上升到地表类以上。在水位 以上，断面水质达到地表类及以上的概率稳定提高。近（年）在星子站水位高于 的 中，有 的天数（）断面水质为地表类，个别天数（十余天）会波动到地表类或地表类。而最近（年）在星子站水位低于 的 里，有 的天数（）断面水质为劣类，只有个别天数能达到地表类或地表类，尚未达到过地表类以上。上述水位 水质关系分析的结果再次表明，水资源是鄱阳湖断面水环境达标的基本条件。已有较多研究表明，鄱阳湖周期性水位变化是水体磷浓度变化的重要原因，枯水期 浓度较丰水期和平水期高，与本研究结果一