近20年来广东省农业面源污染负荷时空变化与来源分析

　　摘要：广东省农业面源污染负荷产生量大，对区域生态环境造成严重影响。采用清单分析法分析了近20年(1999—2019年)广东省农业面源污染负荷时空变化特征，探讨了农业面源污染的来源情况，并分析了农业生产投入强度、农业面源污染负荷和农业面源污染指标的关系。结果表明，近20年广东省农业面源污染总负荷下降6.08%，其中化学需氧量(COD)、总氮(TN)和总磷(TP)的污染负荷增幅分别为−11.88%、4.99%和 26.17%，耕地化肥和农药投入强度分别上升112.19%和60.38%。珠三角地区是广东省农业面源污染负荷最高的地区，其次分别是粤北、粤西和粤东地区。畜禽养殖是COD的主要来源，化肥和畜禽养殖是TN的主要来源，畜禽养殖和水产养殖是TP的主要来源，且水产养殖污染物排放占比呈现出明显上升趋势。不同区域的污染物来源存在一定差异，粤西、粤北和粤东地区COD和TP主要来源是畜禽养殖，TN的主要来源是化肥;珠三角地区水产养殖业成为TN和TP污染负荷的主要来源。广东省面源污染负荷总量下降主要源于城镇化水平的提高和农村人口比例减少。总体而言，广东省面源污染存在时间阶段性变化与空间差异，应当采取全面治理与针对性污染治理相结合的治理措施，全面加强化肥减量措施和养殖业污染治理，重点加强珠三角地区水产养殖污染治理和粤北地区农村生活污染治理。

　　关键词：面源污染;污染负荷;清单分析法;时空变化;污染来源分析

　　葛小君; 黄斌; 袁再健; 王栋栋; 王泉泉; 陈佳村; 谢真越环境科学2021-12-17

　　农业面源污染指的是农业活动过程中畜禽和水产养殖废水排放、化肥与农药流失和农村人类活动等因素导致污染物进入环境造成的污染[1,2]，具有排放随机、过程复杂和不易被监测与控制等特征[3]。中国是世界最大的集约化农业发展中国家之一[4]，从 1978—2017 年，我国化肥年施用量从 884 万 t 增加到 5859 万 t，农药年施用量由 0.73 万 t 增加到 166 万 t，但化肥和农药利用率不足 1/3，畜禽养殖粪污有效处理率不足 1/2[5]。农用化肥和农药的大量投入造成了严重的农业面源污染[6,7]，我国农业的化学需氧量(COD)、总氮(TN)和总磷(TP)排放量分别占全部污染物排放总量的 49.77%、51.97%和 67.22%[8]。

　　广东省是我国农业生产较发达的省份，同时也是农业面源污染问题突出的地区之一[5,9]，全省农业污染排放和农村生活排污占水体 COD、TN 和 TP 排放量的 60.47%、56.74%和 77.21%[10]。近年来，广东省城市化进程加快，城市化率从 1999 年的 31.18%上升到 2019 年的 71.40%，农业产业结构发生较大改变[4]，土地利用格局的改变会使得农业面源污染排放量与排放特征发生变化[11]。但以往广东省农业面源污染研究集中在种植业、畜禽养殖业和农村人口这些方面[12,13]，对水产养殖、农作物固废和农村生活所产生的污染负荷有所忽视 [4,14]，且大多数农业面源污染负荷估算的研究集中在某个流域或者某个典型区域，因此，有必要对广东省的农业面源污染负荷的时空变化与来源进行全面分析与评估，以期为广东省农业生态环境保护及治理提供理论依据。

　　1 材料与方法 1.1 研究方法与数据来源

　　本研究采用清单分析法估算广东省农业面源污染负荷，该方法不用考虑复杂的污染过程，能够在较长的时间和较大的空间尺度上对面源污染负荷进行较准确的评估[13]。空间上将广东省划分为粤东(汕头市、汕尾市、潮州市和揭阳市)、粤西(湛江市、茂名市和阳江市)、粤北(韶关市、河源市、梅州市、清远市和云浮市)和珠三角(广州市、佛山市、肇庆市、深圳市、东莞市、惠州市、珠海市、中山市和江门市)这 4 个区域进行对比研究。将农业面源污染源分为 3 类：第一类是耕作过程中施用的化肥、农药和塑料等化学品以及农作物产品收获后产生的作物废物;第二类是畜禽养殖和水产养殖产生的养殖业污染物;第三类是农村地区居民的生活废弃物[5]。农田固废很少被纳入评价面源污染的指标中，但是由于在还田过程中微生物分解养分和径流冲刷联合作用[14]，农田固废仍具有潜在污染量[15,16]，因此本研究将农田固废作为评价农业面源污染负荷指标之一。为了在统一尺度上比较各类污染物，根据文献[17]确定化学需氧量(COD)、总氮(TP)和总磷(TP)为主要的污染负荷指标，并以化肥、畜禽养殖、水产养殖、农作物固废和农村生活作为农业面源污染的污染源。化肥、农药和地膜的过度使用会增加温室气体排放和加剧耕地退化[16]，但是由于其排污系数难以获得[5]，因此引入排放强度来判断农业耕地投入要素污染现状。

　　本研究数据来源于文献[18,19]。农作物、化肥、农村生活和畜禽养殖的产污和排污系数参考文献[17,23,24];水产养殖产污和排污系数则依据文献中广东渔业各养殖类型占比，选取池塘养殖作为主要养殖方式[25]，并根据文献进行系数修正[17,20]，其中是否引入沉积物对污染物吸附对水产养殖污染负荷排污系数的取值范围影响较大，有研究表明，鱼塘清淤和内源性养分释放同样会产生环境污染[19,21]，因此为了避免高估或低估污染负荷本文中取化学分析法和物料平衡法所给出的两种排污系数的平均值以确定本文水产养殖产污排污系数。由于化肥施用后并不产生 COD 污染，在计算 COD 负荷时，化肥不作为计算指标。具体的计算指标、排污系数及其来源见表 1 和表 2，其中给出主要的污染负荷贡献单位的排污系数。

　　1.1.1 污染负荷计算

　　TPLi= ∑ EUactivity= ∑ ∑ EUclass= ∑ ∑ ∑ EUunit×EUA (1) TPL=TPLCOD+TPLTN+TPLTP (2) 式中，TPLi 为农业面源污染指标 i 的污染负荷量，i 包括化学需氧量(COD)、总氮(TN)和总磷(TP);TPL(Total Pollution Loads)为农业面源污染总污染负荷; EUactivity为农业活动的污染源的污染负荷;EUclass 为各类农业活动污染源中各类别的污染负荷;EUunit 为农业各类农业活动污染源中各类别基本单元的污染负荷;EUA 为 EUunit的排污量，EUA 由以下公式计算出。 EUA= ∑ EUi n i=1 ρi(1-ηi)Ci(EUi , S) (3) EUA= ∑ PEi n i=1 Qi (4) 式中，EUi 为基本单位;ρi 为污染物产生系数;ηi 为相关资源利用效率(%);PEi 为农业面源污染物的产量;Ci 为污染物排放系数，由(EUi, S)所决定的;(EUi, S)为基本单位在不同的区域环境特点、水文条件、降雨、土壤、植被覆盖和各种社会政策的综合作用。该公式可以简化为式(4)，Qi为 i 的排污系数，即农业面源污染物的产量乘以排污系数。

　　1.1.2 投入强度计算

　　EIi= Ei AL (5) 式中，EIi 为 i 的单位投入强度(kg·hm-2)，Ei 为投入要素 i 的施用量，AL 为 i 单位的作用面积(hm2)。根据文献[12,29]和本文数据计算结果，确定化肥和农药面源污染程度分级。针对化肥污染，将污染等级分为 4 个等级，分别是化肥年投入强度小于 225kg·hm-2 的区域为无污染区，225~500 kg·hm-2 的区域为轻度污染区，500~1000 kg·hm-2 的区域为中度污染区，大于 1000 kg·hm-2 的区域为重度污染区;对于农药投入强度，将污染等级分为 3 个等级，施用量小于 10 kg·hm-2 的区域为轻度污染区，10~30 kg·hm-2 的区域为中度污染区，大于 30 kg·hm-2 的区域为重度污染区。

　　1.2 分析方法

　　农业生产要素投入量和农业污染负荷指标间的相关性采用 Pearson 相关分析计算，显著性和极显著性水平分布设定为 P<0.05 和 P<0.01。采用 GraphPad Prism 9.0 进行数据处理和作图。使用 ArcGIS10.7 进行空间分布计算。

　　2 结果与分析 2.1 农业面源污染负荷年际变化特征

　　从图 1 可以看出，广东省农业面源污染总负荷由 1999 年的 236.63104 t 下降到 2019 年的 222.40104 t，降幅为 6.08%。COD 污染负荷由 169.05104 t 下降到 148.96104 t，降幅为 11.88%;TN 污染负荷从 59.92104 t 增加到 62.91104 t，增幅为 4.99%，TP 污染负荷从 7.65104 t 增加到 9.66104 t，增幅为 26.17%。

　　2.2 农业面源污染负荷空间变化特征

　　从广东省不同区域的污染负荷贡献率情况可以看出(图 2)，2019 年总污染负荷贡献率表现为：珠三角>粤西>粤北>粤东。2019 年珠三角 COD、TN 和 TP 污染负荷贡献率相较于 1999 年分别下降了 1.82%、4.19%和 5.23%，但仍为 4 个区域中最高;粤东污染负荷贡献率为 4 个区域中最低，COD、TN 和 TP 污染负荷贡献率下降了 3.06%、3.68%、3.22%;2016 年之前粤西和粤北地区污染负荷占基本保持不变，2016 年后污染负荷贡献率上升趋势较为明显，相比于 1999 年，2019 年污染负荷贡献率分别增加了 4.48%和 1.31%。

　　2.3 农业面源污染无来源特征分析 2.3.1 广东省农业面源污染负荷来源

　　1999~2019 年广东省 COD 的主要来源是畜禽养殖，TN 的主要来源是化肥和畜禽养殖， TP 的主要来源是畜禽养殖和水产养殖，而农作物废弃物排污占比相对较低(图 3)。值得注意的是水产养殖对不同污染物负荷的贡献均呈现出逐年上升的趋势，COD 污染负荷贡献率由 7.25%上升至 15.15%，TN 污染负荷贡献率从 10.07%上升至 16.60%，TP 污染负荷贡献率从 16.89%上升至 23.79%。

　　2.3.2 区域农业面源污染负荷来源

　　粤西和粤北 COD 和 TN 污染负荷贡献率结构较为类似，且 20 年间各污染物污染负荷贡献率变化幅度不明显，畜禽养殖是 COD 和 TP 的主要来源，化肥是 TN 的主要来源(图 4)。珠三角地区畜禽养殖 COD、TN 和 TP 的污染负荷贡献率明显下降，水产养殖占污染负荷贡献率明显上升，相比于 1999 年，2019 年水产养殖 COD 污染负荷贡献率上升 11.71%，TN 污染负荷贡献率达到 25.04%，并且成为了 TP 污染负荷的主要来源。粤东地区污染负荷结构变化不大，畜禽养殖是 COD 和 TP 的主要来源，化肥产污是 TN 主要来源。

　　2.4 农业生产投入强度与面源污染负荷关系 2.4.1 农业生产要素投入时空变化特征

　　近 20 年广东省耕地面积减少 38.44%(表 3)，化肥和地膜的施用量分别增加 30.63%和 96.90%;农药施用量在 2012 年达到最高值(113.88103 t)， 2016 年后开始下降，2019 年农药施用相较于 2016 年量减少了 28.24%。广东省近 20 年化肥平均投入强度属轻度污染状态(333.44 kg·hm-2)，投入强度较高的城市主要有佛山市、深圳市、中山市、广州市、东莞市、茂名市、湛江市、潮州市和汕头市(图 5)，所有城市化肥投入强度均有提高，全省总体提高了 67.61%，达到了中度污染水平(558.87 kg·hm-2)。广东省近 20 年农药平均投入强度属于中度污染(14.81 kg·hm-2)，投入强度涨幅较高的城市主要有佛山市、东莞市、中山市、珠海市、深圳市和潮州市。相比于 1999~2008 年，2009~2019 年农药投入强度提高了 64.23%，农药投入涨幅较高的城市有佛山市、东莞市和潮州市，广州市农药投入强度有所下降。

　　2.4.2 农业生产投入强度与农业面源污染负荷关系

　　相关性分析结果表明(表 4)，COD 污染负荷与 TN 和 TP 污染负荷之间存在极显著相关性(P<0.01);COD、TN 和 TP 污染负荷与化肥、农膜和农药投入强度之间均呈负相关，其中 COD 与化肥投入强度呈显著负相关(P<0.05)。 染负荷计算单位以及城镇化率进行相关性分析(表 5)，其中农业面源污染负荷与面源污染负荷计算单位呈极显著相关(P<0.01)，与城镇化率呈极显著负相关(P<0.01);投入强度与城镇化率呈极显著正相关(P<0.01)，与农村人口数和耕地面积呈极显著负相关(P<0.01)，投入强度与投入总量并不存在显著相关关系，并且投入强度与畜禽养殖产量和粮食作物产量呈负相关关系，其中农药和农膜投入强度与畜禽养殖产量和粮食作物产量呈极显著负相关 (P<0.01)。

　　3 讨论

　　近 20 年来，广东省面源污染负荷变化呈现一定的时间阶段性特征，国内经济体制改革和国际经济形势的变化以及环境保护政策的改变对农业生产和农产品市场的影响是导致污染负荷变化的主要原因[30]。1999—2006 年经历了市场调节增长的阶段，在自由市场机制下，农产品的商品属性得到充分的展现[5]，从而刺激了农业产品的生产，污染负荷随之增加。而在 2007 年污染负荷出现大幅度的降低，这是由于金融危机出现所导致畜禽养殖数量的急剧减少造成的[31]。2007—2015 年期间，随着经济的复苏以及一系列有利的农业政策出台，使得污染负荷缓慢增长。2015—2019 年期间，环境保护政策的完善与实施力度的加大，又使得污染负荷逐年降低。化肥、农药和地膜的使用量受经济形势影响较小，这可能与耕地逐渐转变为建设用地或者撂荒导致的耕地面积减小有关[32,33]，为保证农作物产量投入了更多的化肥和农药。相关行动方案的颁布使得化肥与农药的使用量近几年有所下降[34,35]，但是地膜投入量并未出现降低的趋势，说明种植业生产要素的投入很大程度上受到农业政策的影响。

　　珠三角地区污染负荷贡献率最高，说明珠三角地区存在着较大的农业面源污染问题[36]。珠三角地区作为华南地区经济中心和人口最稠密地区，高速的发展导致农业用地极大地被压缩[4]，而珠三角地区人均肉类产品消费远高于全国平均水平[37]，为了保证供应，化肥的施用量大幅增加并且畜禽养殖的规模不断扩大;其次，考虑到食物的运输距离成本，珠三角周边地区农业生产的积极性较高[4];此外，居民饮食结构的改变也会对蔬菜这类施肥量远高于其他农作物的农产品需求量增加[12,13]。

　　从污染来源方面看，化肥是广东省 TN 污染负荷最主要的来源[5]，畜禽养殖业是 COD 和 TP 污染负荷的主要来源，渔业养殖对于污染负荷的占比呈现明显上升趋势。广东是全国主要的水产养殖区域[38.39]，并且目前水产养殖业规模仍在持续扩大[40~42]，2019 年珠三角地区水产养殖占广东省 41.15%，相比于 1999 年，水产养殖总产量上升 112.98%，然而水产养殖过程中饲料的利用效率普遍较低[43]，导致大量的养分扩散至水域或者滞留在养殖区域内，水产养殖的污染物排放逐渐成为珠三角地区 TN 和 TP 污染负荷的主要来源[42,44]。粤北和粤西地区农业生产方式主要以畜禽养殖和种植业为主，所以污染来源主要是畜禽粪污以及化肥施用。粤东地区由于农村人口比重较高，使得农村生活的污染排放在 COD 污染负荷来源中的占比逐渐上升，此外粤东地区位于沿海地带，水产养殖规模的扩大使得水产养殖排放的总磷占比也逐年上升。

　　城镇化的提高会导致耕地面积的缩减和耕地使用率的降低以及农村劳动力的转移[45~47]，为了增加粮食产量，耕地投入要素强度相应增加，因此城镇化率与耕地投入要素强度呈正相关;但是高城镇化率会导致当地农业生产活动减弱[47]，降低当地的农业面源污染负荷(表 5)。因此城镇化率的增加在减少当地的农业面源污染负荷的同时增加了耕地投入要素强度。目前广东省化肥利用率较低[6]，较高的化肥投入强度会增加 N2O 等温室气体的排放[16]。较高强度的化肥和地膜施用，会导致农田持久性有机污染物和微塑料的大量富集[48,49]。因此在处理高城镇化率的城市农业面源污染问题时，化肥、农药和地膜是重点关注对象[50]，并且需要考虑提高耕地生产投入要素利用效率的问题[51]。

　　基于广东省面源污染负荷时空特征及变化趋势，对于广东省面源污染问题，需针对不同的污染物类型来采取相应的措施。在 COD 污染的管控方面，应着重于畜禽养殖以及农村生活污水处理上，而对于总氮和总磷污染负荷的管控，要从化肥的施用管理入手，其次应减少散户畜禽养殖并逐步完善规模化畜禽养殖业粪污处理系统。最后，需要增加水产养殖面源污染排放的监管和治理力度[52,53]，完善水产养殖业中污染物扩散机理，加强水产养殖污染控制技术的研究，制定有效的预警防治措施。此外，农业生产过程中产生的新型污染物(例如抗生素、微塑料和持久性有机污染物等)也需要引起重视[54,55]。

　　4 结论

　　(1)近 20 年来，广东省农业面源污染总污染负荷和 COD 负荷有所下降，TN 和 TP 的污染负荷上升。农业面源污染负荷排放存在阶段性变化特征，经济和农业政策对农业面源污染负荷影响较大，尤其是对种植业生产要素的影响较为明显。

　　(2)从空间上来看，珠三角地区是广东省农业面源污染负荷最高的地区。城镇化率与农村人口和农业面源污染指标呈现极显著关系。

　　(3)畜禽养殖是 COD 的主要来源，化肥和畜禽养殖是 TN 的主要来源，畜禽养殖和水产养殖是 TP 的主要来源，水产养殖污染物排放占比呈现出明显上升趋势，且不同区域的农业生产模式造成的污染物来源存在差异。

　　(4)广东省面源污染负荷总量虽有下降，但农业生产投入强度仍然呈上升趋势，面源污染问题依旧严峻。