鄱阳湖区乐安河流域营养盐负荷影响因素分析*

夏 雨，鄢帮有，方 豫

(1：南昌大学理学院，南昌 330031) (2：江西省山江湖开发治理委员会办公室，南昌 330046) (3：江西省科学技术厅，南昌 330046)

鄱阳湖区乐安河流域营养盐负荷影响因素分析*

夏 雨1，2，鄢帮有3**，方 豫2

(1：南昌大学理学院，南昌 330031) (2：江西省山江湖开发治理委员会办公室，南昌 330046) (3：江西省科学技术厅，南昌 330046)

鄱阳湖流域；乐安河流域；农业面源污染；富营养化

湖泊水资源作为一种与人类生存和发展密切相关的独特资源，在国民经济和社会可持续发展中占有重要位置[1].河流是连接流域与湖泊的重要纽带，与流域、湖泊之间有着重要的物质与能量交换，河流的水质特点影响着下游湖泊的水质状况[2].

鄱阳湖是中国第一大淡水湖，是江西的母亲湖，承载着保护全球生物多样性、调蓄洪水及供给淡水的生态服务功能，为江西经济社会发展提供了重要的物质基础和生态资源.鄱阳湖被誉为“一湖清水”，但其水环境仍面临严峻挑战.据《1999-2009年江西省环境状况公报》等资料显示，鄱阳湖水系和湖体水质总体都呈下降趋势.1980s鄱阳湖水质以Ⅰ、Ⅱ类为主，平均占85%；1990s仍以Ⅰ、Ⅱ类为主，平均占70%；2003年以后，Ⅰ、Ⅱ类水只占50%；2006年后，鄱阳湖水全年优于Ⅲ类的不到六成，属于Ⅲ类的有两成多，劣于Ⅲ类的则逼近两成，整体上呈现出中度营养化的状态，且枯水期富营养化程度比丰水期更为严重.

湖泊的污染主要包括流域内的点源(工业污水排放、城市污水等)和面源(农业面源、流域水土流失、大气降尘等).随着对水体点源污染的有效控制，水污染负荷结构发生了明显变化，面源污染已成为水环境的一大污染源或首要污染源[3].抚仙湖[4]、滇池[5]、丹江口水库[6]、三峡水库[7]、太湖[8]和鄱阳湖[9]等都不同程度遭受面源污染.

面源污染中的氮、磷营养物累积是引起下游水体富营养化的主要原因[10].农业面源污染影响因素复杂[11]，包括土地利用方式、农田耕作、农事活动和田间水肥管理、地形地貌、土壤植被、气候、水文特征[12].研究表明我国已成为世界上藻类水华暴发最严重的国家之一[13].近年来，很多湖泊都开展了营养盐调查研究[14-16]，但由于面源污染起源于分散、多样的地区，地理边界和发生位置难以识别和确定，随机性强、成因复杂、潜伏周期长[17]，在富营养化评估方面，数据获取困难，方法的集成与模型研发还较为落后[18-19].

鄱阳湖的污染物主要由五河(赣、抚、信、饶、修)输入，五河输入的总磷(TP)和总氮(TN)分别占入湖总量的80.5%和66.4%[20]，因此五河等主要入湖河道的水质对鄱阳湖水环境有重要影响.乐安河是鄱阳湖重要的入湖水系，受沿岸工业废水、生活污水等的污染，该河的许多污染指标超标，其污染甚至波及到鄱阳湖的水质[21-22].本文选择主要营养盐指标为研究对象，通过野外实测数据和数学分析对该流域水系中营养物时空分布及影响因素等特征进行定量分析，从而评估乐安河流域内水质富营养化水平的现状，为制定针对性污染控制措施提供科学依据，从而有效消减入湖营养物负荷，保护鄱阳湖“一湖清水”.

1 研究区概况

乐安河，位于江西东北部(图1)，是鄱阳湖第4大入湖水系饶河的最大支流，发源于江西省婺源县鄣公山南麓，自东北向西南流，石镇街以上集水面积8367km2，主河长279km，多年年平均径流量为99.67×108m3.其流域是鄱阳湖区降水最多的地区之一，年平均降水量1687～1870mm，但季节分配不均.流域土壤属赣东北山地丘陵，农业土壤主要是水稻土.近年来流域内面源污染产出量迅速增加，河道水质下降趋势加剧，已呈现中营养水平[23].

在乐安河流域，土地利用覆盖类型多样，从图2可以看出，区域高程相差约1450m，呈东北部高和西南部低的地势，所经地区以丘陵、山地为主，林地和耕地面积之和占研究流域总面积接近95%，无大城镇分布，仅有零散分布的居民点.

图1 乐安河流域在江西省的位置Fig.1 Location of Le’an River watershed in Jiangxi Province

图2 乐安河流域DEM图Fig.2 Digital Elevation Model of Le’an River watershed

图3 乐安河流域监测点分布Fig.3 Distribution of monitoring stations in Le’an River watershed

2 数据与方法

2.1 实验设计

表1 相关矩阵Tab.1 Correlation Matrix

TNNO-3-NNH+4-NTPTN1.00000.90640.77690.1956NO-3-N0.90641.00000.58780.0172NH+4-N0.77690.58781.0000-0.0045TP0.19560.0172-0.00451.0000

表2 相关矩阵的特征值Tab.2 Eigenvalues of the Correlation Matrix

指标特征值差值贡献率累计贡献率TN2.532254351.518966280.63310.6331NO-3-N1.013288060.590773210.25330.8864NH+4-N0.422514850.390572110.10560.9920TP0.031942740.00801.0000

2.2 各营养盐主成分分析

首先采用主成分分析法对主要营养物质进行分析，通过相关矩阵得到各物质形态之间的相关性.

2.3 各营养盐物质形态聚类分析

图4 17个监测点水质的平均连接树状结构图Fig.4 Tree structure chart for water quality of 17 monitoring stations

聚类结果将采样点根据富营养化程度划分为4类：低富营养区、较低富营养区、中富营养区和重富营养区.从图5可以明显看出从上游到下游富营养化程度总体呈不断升高的趋势.其中LZ5121位于流域的最上游，该点上游地区林地与草地面积总和占95%以上，营养盐产出少且被较好地截留和消减，平均水质可达到Ⅱ类.在黄色类型中，LG511和LG516上游集水域林地和草地面积总和占80%，LG513、LG514、LZ509、LZ510的集水域中林地和草地面积总和超过90%，因此仍维持较好水质，TN浓度在0.6mg/L左右，略高于Ⅱ类水.橙色类型区域位于流域中部，TN浓度上升到0.7～0.8mg/L，为Ⅲ类水，LZ5041和LZ505两个子流域中耕地面积分别为27%和30%，林地为64%和65%，比上游地区明显下降，而耕地又是面源污染物的主要来源.LZ507和LZ515子流域的林地面积均为85%，但在采样时发现，LZ507监测点附近有化工企业，因此该地的高污染物输出强度应该来源于点源排放；而LZ515位于赋春镇下游，该镇有茶园3.33km2，大面积的茶叶种植必然对下游水质造成影响.LG508是干流上监测点，接受中上游地区的污染，污染浓度相应上升.同理，LG501和LG506是干流上最后两个点，流域内所有污染都汇集于此处，氮浓度明显提高，平均高于0.9mg/L，接近Ⅳ类水.

图5 水质空间分布Fig.5 Spatial distribution of water quality

3 主要研究结果

影响面源污染的因素众多，而降水状况和人类活动可在短期内对污染物浓度产生显著影响.首先根据每个监测点TN浓度的年平均值分析乐安河流域氮浓度的年际变化，可反映降雨量及降雨强度对面源氮浓度的影响.

图6 各监测点的TN年平均浓度Fig.6 Annual average concentration of TN in 17 monitoring stations

鄱阳湖流域在枯水期和丰水期的降雨量和径流量变化巨大，且对应不同的农事季节，对面源污染物产出及输移影响极大.本研究将每年的4-9月定为丰水期，10月至次年3月定为枯水期.并根据流域内地形地貌将17个子流域分为两类：上游东北部低中山区(LG508、LZ509、LZ5092、LZ510、LG511、LZ5121、LG513、LG514、LZ515、LZ516)，这10个子流域以林地为主，位于流域的上游，林地面积都超过80%，而耕地面积只有10%左右；下游西南部丘陵区(LG501、LZ502、LZ504、LZ5041、LZ505、LG506、LZ507)，这些子流域耕地面积都超过30%.将不同类型的监测数据与径流量和时间进行对比.

图7 上、下游监测点氮平均浓度Fig.7 Average nitrogen concentration in upstream and downstream

从图6和图7可见，在取样监测的3年中，TN浓度不断升高.2009年江西省平均降水量1392.0mm，比多年平均少15.0%，属偏枯水年份.流域内径流量相应减少，营养物流失较少，这一年大部分监测点的平均水质都在Ⅱ类以内，有6个监测点的平均水质为低水平Ⅲ类，Ⅲ类水主要出现在径流量最小的12月.

2010年江西省极端天气过程频繁，暴雨强度大，主汛期全省平均降雨量比1998的大洪水年同期还多14%，全年降水量多年平均增加27.3%，属丰水年份.降雨形成地表径流携带大量营养物进入河道，导致水体中氮浓度增加，所有监测点只在8月水量最大时平均值为Ⅱ类水质，全年平均为Ⅲ类水质.

2011年鄱阳湖流域水文条件特殊.1-5月，江西省降水量比常年同期均值减少49%，为1950年以来同期最少；全省年平均降水量比年均值少20.4%，属枯水年份.而6月份，北部出现典型的“旱涝急转”现象，平均降水量较常年同期平均偏多183%，瞬间增大的地表径流量，使流域中的营养物质在短时间内带入河道水体，所有监测点氮平均浓度在这个月急剧升高，并在7月达到全年的高峰值，导致这一年所有监测点的氮平均浓度显著升高，仅有两个点为Ⅲ类，8个点为Ⅳ类，6个点都为Ⅴ类，还出现一个劣Ⅴ类水质点.

根据上、下游地区不同时期的采样值可以看出，下游的数值明显高于上游，说明林地对面源污染有消减和抑制作用，而耕地是主要的面源污染产出源.在2009和2010年两个水文相对正常的年份，上、下游氮浓度变化趋势基本一致，都是从4月开始显著上升，在6月后又迅速下降，并在8-9月达到一年中的最低值，之后又急剧上升，在12月份达到最大值.

TP浓度的时空分布规律不如氮的明显.总体上来说，也是从上游到下游逐渐升高，枯水期高于丰水期.由于磷主要以颗粒态存在，其浓度受降雨影响程度低于氮，而河底沉积物的释放也使其浓度分布规律性较差.

可见降雨量及降雨强度是乐安河流域营养盐浓度年际变化的主要影响因素，而流域内的农业活动是流域营养盐浓度年内变化的主要影响因素.

4 结论

2) 降雨对营养物浓度有重大影响，降雨量和降雨强度的增加都会明显促进营养物的流失，因此偏枯水的2009年，年平均水质维持在Ⅱ类左右；而在丰水的2010年，年平均水质则下降至Ⅲ类.但同时，长时间的降雨也会使湖泊水位升高，水体容量增大，对污染物能起到一定的稀释作用.因此在极枯的年份，由于水量大量减少，即使少量污染物也能呈现较高浓度；而短时间内的强降雨则会使地表径流瞬间增大，携带大量污染物进入水体，使污染物浓度急剧升高，这就导致2011年全年营养物浓度都偏高，并在6-7月异常增大，出现劣Ⅴ类水质.说明降雨是流域水质年际变化的主要原因.

3) 通过一年中不同月份降雨量、人类活动及氮浓度的对比分析可以看出，在主汛期，湖泊河流水量达到最大时对污染物浓度稀释作用影响较大；而其他时期，农业活动强度与氮浓度基本呈正比，是引起流域内营养盐浓度年内变化的主要原因.这也表明该流域内社会经济发展使污染物排放大量增加，在水文正常年份，充足的水量极大地稀释了污染物浓度，使流域呈现较好水质，而在枯水年甚至枯水期水质则显著下降.

综上所述，流域内产生并随水流进入鄱阳湖的污染物是影响其水质的主要原因，而随着社会经济发展，污染物排放还在大量增加；且随着全球气候波动引起的极端气候事件频发和长江流域重大水利工程的建设运行，鄱阳湖及其周边湖泊缺水将呈现常态化，也将导致水环境容量下降.因此仅靠湖泊的自净能力已远远不能消耗流域内产生的污染物，因此，在五河入湖口实施污染物削减措施和工程，控制入湖污染物总量，是保护鄱阳湖“一湖清水”的重要途径.

[1] 李世杰，窦鸿身，舒金华等.我国湖泊水环境问题与水生态系统修复的探讨.中国水利，2006，13(6)：14-17.

[2] 马 倩，刘俊杰，高明远.江苏省入太湖污染量分析(1998-2007年).湖泊科学，2010，22(1)：29-34.

[3] 胡 斌，树 庄，陈春瑜等.滇池流域土壤氮磷分布特征及关键影响因素研究.土壤学报，2012，49(6)：1178-1184.

[4] 夏天翔，潘继征，刘雪华等.抚仙湖水体NP变化及其非点源污染特征.农业环境科学学报，2008，27(4)：1340-1345.

[5] 孙治旭，王红华，严婷婷.滇池沿岸农业面源污染治理模式研究.农业环境与发展，2008，25(5)：87-89.

[6] 黄文敏，毕永红，胡征宇等.丹江口水库农业面源污染特征研究.环境科学与管理，2012，37(1)：33-38.

[7] 王丽靖，郑丙辉，李子成.三峡库区及上游流域面源污染特征与防治策略.长江流域资源与环境，2009，18(8)：783-788.

[8] 张红举，陈 方.太湖流域面源污染现状及控制途径.水资源保护，2010，26(3)：87-90.

[9] 吴罗发.鄱阳湖区农业面源污染时空分布研究.江西农业学报，2012，24(9)：159-163.

[10] Wu HM， Zhang J， Li PZetal. Nutrient removal in constructed microcosm wetlands for treating polluted river water in northern China.EcologicalEngineering， 2011， 37(4)：560-568.

[11] Sun C， Wu HJ. Pollution from animal husbandry in China： a case study of the Han River Basin.WaterSciencesandTechnology， 2012， 66(4)：872-878.

[12] 陶 春，高 明，徐 畅等.农业面源污染影响因子及控制技术的研究现状与展望.土壤，2010，42(3)：336-343.

[13] 王 娟.藻类水华的发生及控制技术研究现状及展望.安徽农业科学，2011，39(4)：2212-2214.

[14] 赵林林，朱广伟，顾 钊等.太湖水体氮、磷赋存量的逐月变化规律研究.水文，2013，33(5)：28-33，45.

[15] 谢琳萍，蒲新明，孙 霞等.荣成湾营养盐的时空分布特征及其影响因素分析.海洋通报，2013，32(1)：19-27.

[16] 黄海宁，林 彩，暨卫东等.九龙江河口区营养盐分布特征及其影响因素分析.海洋通报，2012，31(3)：290-296.

[17] 王建兵，程 磊.农业面源污染现状分析.江西农业大学学报：社会科学版，2008，7(3)：35-39.

[18] 胡春华，楼 倩，丁文军等.鄱阳湖氮、磷营养盐的滞留效应研究.环境污染与防治，2012，34(9)：1-4.

[19] 张天锡，倪跃晓.九峰水库富营养化评价及入库营养盐总量控制研究.中国科技信息，2012，(7)：52.

[20] 闵 骞.鄱阳湖水环境现状与思考.鄱湖人家网，2008.9.http：//www.cnplph.com/sf_901ADCB56A864DE0947AE6F51 6758C8D_187_cnplph.html.

[21] 黄学平，刘 峰，万卫红.水环境容量与总量控制及其在乐安河的研究应用.西部资源，2011，(6)：76-77.

[22] 曾 艳，万金保.乐安河流域管理的研究.大江周刊：论坛，2010，(2)：11-12.

[23] 高海鹰，庄 霞，张 奇.鄱阳湖乐安河流域非点源氮污染时空变化特征分析.长江流域资源与环境，2011，20(5)：597-602.

[24] 李正兆，高海鹰，张 奇等.抚仙湖流域典型农田区地下水硝态氮污染及其影响因素.农业环境科学学报，2008，27(1)：286-290.

[25] Withers PJA， Sharpley AN. Preface： Characterization and apportionment of nutrient and sediment sources in catchments.JournalofHydrology， 2008， 350：127-130.

Nutrient loading and its controlling factors in Le’an River watershed， Lake Poyang basin

XIA Yu1，2， YAN Bangyou3& FANG Yu2

(1：SchoolofScience，NanchangUniversity，Nanchang330031，P.R.China)(2：OfficeforMountain-River-LakeDevelopmentandManagementCommitteeofJiangxiProvince，Nanchang330046，P.R.China)(3：DepartmentofScienceandTechnology，JiangxiProvince，Nanchang330046，P.R.China)

Lake Poyang watershed； Le’an River watershed； diffuse pollution； eutrophication

*江西省科学合作计划项目(20142BDH80016)和江西省主要学科学术和技术带头人培养计划专项(20123BCB22008)联合资助.2014-02-21收稿；2014-08-20收修改稿.夏雨(1980～)，女，博士，副研究员；E-mail：416522685@qq.com.

**通信作者；E-mail：ybangyou@sina.com.