中国省域化肥面源污染时空格局演变与分组预测

减少化肥面源污染的同时保持农业产值持续增长是实现农业产业生态化和农业高质量发展的必然要求。中国各省均制定并实施了化肥零增长行动计划,但进展和成效并不一致,并可能相互影响。论文运用化肥流失系数法对中国1997—2018年31个省(市、自治区)化肥面源污染排放强度进行核算,再运用空间自相关和热点分析对其进行时空格局演变分析,揭示化肥面源污染的时空演变规律,探讨区域间的相互影响。根据时空格局特征将全国分为热点区、冷点区和非热(冷)点区,在考虑相邻省份间空间异质性和相关性的条件下,分组模拟和预测化肥面源污染排放强度与人均农业产值间的环境库兹涅茨曲线(EKC)时间路径。结果表明:1)化肥面源污染排放强度省际差异较大,表现为空间正自相关,呈集聚模式。热点分析显示,化肥面源污染时空格局相对稳定,热点区主要集中在中南部,长江中下游地区尤其显著,黄淮海地区近年热点程度下降较明显;冷点区主要集中在西部地区和黑龙江。2)基于时空格局分组的环境库兹涅茨曲线(EKC)趋势模拟表明,各组均存在显著的非线性EKC关系但趋势和拐点差异明显,热点区为"倒U型",冷点区和非热(冷)点区为"倒N型",多数省份正处于曲线上升阶段且距拐点较远。3)产业结构的调整和转移促使区域间存在化肥面源污染空间溢出效应,要从整体上把握区域间的协同治理。根据研究结果,提出热点区应研发推广适用施肥设备,提高化肥利用率;冷点区应保护性耕作,增施有机肥;非热(冷)点区应合理调整农业产业结构,注重种养循环。区域间则应当通过生态补偿、排污权交易等方式实现协同治理。     

洞庭湖平原耕地面源污染时空特征及脱钩效应

[目的] 揭示2009—2019年洞庭湖平原耕地面源污染时空演变规律,探究耕地面源污染与农业生产的脱钩关系,为洞庭湖平原耕地面源污染治理提供科学参考。[方法] 采用化肥农药施用强度和粮食单产分别描述洞庭湖平原耕地面源污染和农业生产状况,利用空间分析法和脱钩模型揭示2009—2019年洞庭湖平原耕地面源污染和农业生产的时空演变规律及其脱钩关系。[结果] ①2009—2019年洞庭湖平原化肥、农药施用强度均呈现“先升后降”的趋势,且具有明显的空间差异性,化肥、农药生产率总体上呈现出上升趋势; ②2009—2019年洞庭湖平原粮食单产呈现出倒U形趋势,且具有明显的空间差异性; ③洞庭湖平原粮食单产与化肥农药施用强度的脱钩关系呈现出明显的阶段性特征：第一阶段为波动期(2009—2013年),第二阶段为强脱钩期(2014—2015年),第三阶段为弱负脱钩和衰退脱钩期(2016—2019年); ④洞庭湖平原农业生产与耕地面源污染的脱钩关系在县级空间尺度上逐渐呈现出较高的趋同性和聚集特征。[结论] 应进一步加强耕地面源污染源头防治,并采取多种举措提高粮食单产,最终实现耕地面源污染与农业生产的强脱钩。     

中国粮食主产区化肥施用时空特征及生态经济合理性分析

研究粮食主产区化肥施用状况能够为化肥减量增效、防治化肥损失引起的面源污染、保障国家粮食安全提供决策参考。该研究在分析中国粮食主产区化肥施用强度的时空变化特征的基础上,根据大田粮食作物对氮磷钾的养分要求,提出化肥施用生态经济合理性评价方法,并对中国粮食主产区化肥施用生态经济合理性进行评价。研究表明:从1993到2017年,中国粮食主产区化肥施用强度从205.5 kg/hm~2增加到319.9 kg/hm~2,其中氮、磷、钾化肥的施用强度分别从133.8、48.0、23.7 kg/hm~2增加到162.6、84.2、73.2 kg/hm~2;总化肥及氮、磷、钾化肥施用生态经济适宜量分别从217.4、108.7、54.4、54.4 kg/hm~2增加到300.6、150.3、75.2、75.2 kg/hm~2;氮肥施用强度在1993年已经超过了其生态经济适宜量,磷肥从2001年开始超过了其生态经济适宜量,钾肥在2009年后逼近其生态经济适宜量;1993-2007年间,河南、山东、江苏一直是高强度施肥省份,其他次高强度施肥的省份由高强度施肥省份逐步向周围扩散。2017年,河南为化肥施用生态经济不合理区,化肥施用明显过量;河北、江苏、安徽、山东和湖北则属于经济合理生态不合理区,化肥施用有所过量;湖南、吉林、辽宁和内蒙古施肥适当,属于施肥生态经济合理区;黑龙江、四川和江西为施肥生态合理经济不合理区,化肥施用不足。因此,中国粮食主产区化肥施用强度及生态经济合理性差异较大,总体施肥过量,部分地区存在施肥不足的情况。节肥增效、防治农业面源污染要坚持分类指导、因地减量的原则,完善配方测土施肥政策和技术手段,大力推进信息技术和有机肥利用,走有机无机肥相结合道路。     

河南省农业面源污染的影响因素分析

河南省是典型的农业大省,在粮食高产的背后,农业面源污染十分严重,威胁着地表水环境的安全。通过构建实证分析模型发现,农业经济规模、农户环保支付意愿、种植业结构均对农业面源污染排放量和排放强度有正向影响,农业结构对农业面源污染排放量有负向影响;农业技术进步对农业面源污染排放量有正向影响,对农业面源污染排放强度有负向影响。基于此,提出如下对策建议:构建农业面源污染与工业点源污染的"双控"模式;制定农业经济发展与面源污染治理的"双赢"政策;制定农业经济发展与面源污染治理的"双赢"政策;建立技术创新和制度创新的"双创"途径;打造政府引导与社会公众广泛参与的多管齐下治污措施。     

化肥施用技术与政策的国际对比

主要内容为化肥施用技术与化肥施用管控政策的国际比较。改革开放以来,化肥的广泛利用支撑我国粮食产量快速增长,化肥在保障我国粮食安全中扮演了重要角色。但化肥也是水体和土壤污染源的主要起因。21世纪以来,化肥不合理施用造成的农地和食品污染,已经成为我国亟待解决的重大环境问题。科学合理的利用化肥,事关我国生态、粮食和食品安全。2015年,中央的"一号文件"中明确指出:加快转变农业发展方式,从主要追求产量和依赖资源消耗的粗放经营到数量质量效益并重;粮食要"量"增,更要"质"优;加强农业面源污染治理,深入开展测土配方施肥。     

洞庭湖平原化肥面源污染生态环境风险及脱钩效应

[目的]推动实现农业经济增长与化肥面源污染生态环境风险的脱钩,为有效治理耕地面源污染和耕地可持续利用提供理论支持。[方法]基于生态环境风险评价模型,从县域尺度评价洞庭湖平原各地区化肥面源污染生态环境风险程度,采用非参数核密度估计洞庭湖平原化肥面源污染生态环境风险区域差异,进而采用脱钩模型进一步分析了洞庭湖平原农业经济增长与耕地化肥面源污染生态环境风险的脱钩关系。[结果](1) 2009—2019年洞庭湖平原化肥面源污染生态环境风险指数呈现先升后降的趋势,化肥面源污染生态环境风险指数先由2009年的0.749升至2014年的0.762,然后降至2019年的0.732,经历了由中度风险转为高度风险,后又降为中度风险的过程;中部地区污染更严重;(2)化肥面源污染生态环境风险区域差异明显,呈现两极分散化集聚特征,内部差距先扩大后缩小;(3)洞庭湖平原农业经济增长与化肥面源污染生态环境风险的脱钩关系呈现出明显的阶段性特征：第一阶段为波动期(2009—2013年),第二阶段为强脱钩期(2014—2019年)。[结论]在面源污染治理工作的推进下,洞庭湖平原化肥面源污染生态环境风险经历了先上升后下降的...     

2005-2019年洞庭湖平原耕地面源污染迁移轨迹及空间格局

洞庭湖平原因不合理的农业生产活动导致严重的耕地面源污染问题，揭示洞庭湖平原近10余年耕地面源污染迁移轨迹及空间格局，可为洞庭湖平原下一步科学防治耕地面源污染提供参考。该研究以洞庭湖平原湖南省部分的21个区（县、市）为实例，利用清单分析法构建农田化肥、人畜排泄物和农田固废3类污染单元，测算化学需氧量（Chemical Oxygen Demand，COD）、总氮（Total Nitrogen，TN）和总磷（Total Phosphorus，TP）3类污染物排放情况，并运用空间分析、重心模型和冷热点分析研究2005－2019年洞庭湖平原耕地面源污染轨迹迁移变化及空间格局。结果表明：1）从时序上看，2005－2019年间3类污染物排放总量和排放强度总体上均呈现下降趋势；从空间分布上看，COD排放强度等级分布更为集中，洞庭湖平原耕地面源污染程度呈现出"东部>西部>中部"的空间格局；2）从重心模型分析来看，3类污染物排放强度重心在东西方向上的移动更显著，并且重心移动速度呈现阶段变化；3）从冷热点分析来看，3类污染物空间分布格局差异性显著，COD最为明显，TP次之，TN最小，热点区域大多分布于洞庭湖平原东部，冷点区域大多分布于洞庭湖平原中部。研究说明了在生态文明建设过程中，洞庭湖平原耕地面源污染情况有所好转，并且各地区耕地面源污染程度和治理都具有差异性，为今后因地制宜治理耕地面源污染提供理论依据和思路。     

黄石市农业面源污染的解析及其空间异质性研究

为弄清黄石农业面源污染的源分配以及在空间上的分布差异,进而为削减黄石农业面源污染提供理论依据,该文采用清单分析方法,研究了黄石市化肥施用、有机肥施用、秸秆遗弃、畜禽养殖、水产养殖、生活污水、生活垃圾、地表径流中COD、BOD5、TN、TP的污染负荷及其在空间上的分布差异.结果表明,黄石市农业面源污染排放COD、BOD5、TN、TP的实物总量分别为6.12万、3.55万、1.65万、0.43万t/a,相应的等标排放量分别为3061、8874、16542、21726 f/a.农业面源污染负荷已占全市水环境的半壁江山,接近水环境容量,必须予以削减.主要污染源是畜禽养殖、化肥与生活污水,主要污染物是TP、TN.因农业面源污染引起COD、BOD5、TN、TP的排放浓度为26.04、15.10、7.04、1.85 mg/L,达到严重污染.按压力确定的污染程度是大冶市>市区>阳新县,按响应确定的污染程度是市区≈大冶市>阳新县.     

安徽省巢湖流域众兴水库小流域农业面源污染调查与评价

[目的]对安徽省巢湖流域众兴水库小流域农业面源污染情况进行调查分析和评价,了解该地区农业面源污染的现状和特征,以期为制定针对性防控措施提供参考依据。[方法]应用等标污染负荷法,通过问卷调查形式对巢湖流域众兴水库小流域内进行实地的走访调研,对水库周边3个村域的种植、养殖及人居面源污染进行分类调查和核算分析。[结果](1)众兴水库项目区周边3个村落面源污染排放源中,养殖业污染物的排放量和污负荷率最大,污染物排放量达到109.00t,污负荷率最大为57.52%;种植业化肥施用带来的面源污染物总量达25.52t,负荷率为13.47%。生活污染源污染负荷为29.01%,污染物总量达到54.97t。(2)在化学需氧量(COD),TN,TP这3个评价因子中,以COD的污染负荷率最高,达74.26%;TN污染负荷率为21.13%;TP的污染负荷率为4.61%。[结论]小流域内养殖业的污染物排放是威胁众兴水库水源地饮用水安全的主要构成因子,也是农业面源污染防控的重点。     

面源污染下三峡库区农业生态环境效率及影子价格测算

农业的现代化进程产生了较为严重的农业面源污染,将农业的生态环境效应纳入农业生产效率评价框架是研究的新方向。采用非径向、非角度的SBM-DEA模型,结合方向性距离函数,考察了2003-2015年期间重庆三峡库区农业生态环境效率,引入影子价格测算了各区(县)农业面源污染的边际减排成本。研究结果表明:2003-2015年期间,重庆三峡库区农业生态环境效率趋于稳步提升态势,但整体水平较低,平均值仅为0.508,全样本期投入无效率和非期望产出无效率平均值分别为27.3%和24.7%。农业生态环境效率较高的地区集中在都市核心区和库尾都市外围区,江北和大渡口的农业生态环境效率值为0.7;而库尾低山丘陵区、库中平行岭谷区、腹地山地丘陵区和腹地外围山地区的农业生态环境效率相对较低,尤其是巫溪、万州与开县农业生态环境效率值仅为0.4左右。从效率损失的原因看,化肥、农药的过度投入及其相应的逆生态环境效应表现突出,在投入无效率的来源分解中贡献最大,达35.82%;大渡口、南岸、江北、北碚与渝北农业面源污染排放的影子价格居前,分别为3.503、3.478、3.370、3.317、3.302亿元/万t,污染减排空间较小;万州、丰都与开县农业面源污染排放影子价格较低,分别为2.167、2.321和2.492亿元/万t。在制定区域排放目标时,应考虑区域农业生态环境效率的空间异质性,针对不同区域的农业生态环境效率和影子价格制定具有针对性的减排策略,尤其应当把生态环境效率改进余地较大和边际减排成本相对较低的地区作为减排重点地区。     

中国耕地利用净碳汇时空演进及综合分区

耕地利用净碳汇是影响耕地利用进程的重要因素，分区优化耕地利用是实现中国农业绿色、可持续发展的有效手段。基于耕地利用碳源-碳汇视角测算了2005-2020年中国31省（市、区）耕地利用净碳汇量，利用核密度估计分析其动态演进趋势，并探究耕地利用净碳排放与农业经济的脱钩效应；结合脱钩结果和粮食产区进行综合分区优化，借助灰色模型预测了2021-2035年各综合分区的净碳汇量及脱钩状态。结果表明：1）研究期内耕地利用净碳汇量呈波动式上升，年均增速3.740%，粮食主产区碳汇能力优于其他粮食产区；2）核密度估计显示，中国耕地利用净碳汇地区差距逐渐缩小，但仍存在区域发展不均衡，粮食主产区、产销平衡区和主销区均呈向碳汇高值区的扩散态势；3）研究期内耕地利用净碳排放与农业经济脱钩状态主要为强脱钩，2020年除北京、上海、广东和海南外，其他省（市、区）均处于实现经济效益和生态效益协调的综合分区；4）预测2021-2035年中国耕地利用碳汇能力保持平稳增长，各综合分区脱钩状态主要为强脱钩。并对全国及各综合分区提出针对性优化策略，推进实现耕地利用低碳、绿色、高质量发展。     

基于缺口率模型的区域粮食安全风险评估

深入分析各区域的粮食供需趋势对于保障中国未来的粮食安全至关重要。该文首先分析了影响中国粮食供需安全的风险源,提出了粮食主产区、粮食主销区和产销平衡区的粮食缺口率风险评估模型并划分了相应的风险等级,然后以洞庭湖区为例进行了实证分析。结果表明:自然灾害是影响洞庭湖区粮食供给安全的最主要因素,人口增长则是影响该区域粮食安全的普遍因素;对于大多数的县/市,资源和投入要素对未来的粮食供给将产生积极的影响。当自然灾害为一般情景或乐观情景时,洞庭湖区的粮食安全状况较为乐观;当自然灾害为悲观情景时,湖区部分县/市的粮食安全形势非常严峻。但就湖区整体而言,通过湖区内部调剂即可以满足高风险县/市的粮食需求。研究结果可为相关部门的政策制定提供客观依据,为粮食安全的风险预警和防范提供参考。     

基于 SWAT模型的小流域面源污染负荷时空分异研究--以重庆市万州区陈家沟小流域为例

以三峡库区陈家沟小流域为研究区,通过实地监测调查和资料收集建立空间和属性数据库,采用2006、2007年监测数据建立SWAT模型,并进行模型参数率定,利用2008、2009年的监测数据进行模型验证。应用验证后的SWAT模型模拟分析研究区面源污染负荷的时空分异规律以及不同土地利用类型污染负荷的差异性,结果表明：空间上,小流域内不同子流域的污染负荷存在差异;时间上,污染负荷分布与降雨周期、农药化肥施用时间密切相关;不同土地利用类型面源污染负荷存在差异,耕地负荷最大,林地最小。针对面源污染负荷的时空分异性,综合考虑土地利用方式和小流域生态、经济要求,提出了小流域面源污染防治措施配置建议。     

内蒙古河套灌区粮食生产灰水足迹评价

随着面源污染和土壤盐渍化成为重要的环境污染源和粮食安全的制约要素,如何评价和量化农业生产的负面效应,成为亟需解决的问题。水足迹理论的出现使该负面效应的量化成为可能,水足迹包括蓝水、绿水和灰水足迹,灰水足迹可以表征不同类型负面效应的大小,粮食生产灰水足迹反映了单位粮食产量的负面效应。该文应用水足迹理论,以内蒙古河套灌区为研究区,给出粮食生产灰水足迹的计算方法,选取环境最大允许浓度Cmax和本底浓度Cnat,逐项计算各项灰水足迹,并根据短板原理得出总灰水足迹,计算分析河套灌区粮食生产灰水足迹。结果表明:2005-2008年面源污染的灰水足迹为0.55~0.58亿m3;积盐的灰水足迹从2005年的4.570亿m3减少到2008年的1.825亿m3。总灰水足迹从2006年的5.872亿m3,减少到2008年的1.825亿m3,总灰水足迹在总水足迹中比例小于10%,有逐年降低的趋势。2005-2008年粮食生产灰水足迹分别为0.129、0.159、0.062和0.043 m3/kg,粮食生产灰水足迹有逐年降低的趋势,2008年的粮食生产灰水足迹仅相当于2006年的27.04%。节水灌溉等新技术的推广是其主要原因,节水灌溉可以减少无效灌溉水量和水分的无效蒸发量,从而减小粮食生产灰水足迹。在此基础上,给出减少粮食生产灰水足迹的措施,即节水灌溉、种植业结构调整、合理确定地下水位和合理使用化肥、农药。研究成果较好地量化了大型灌区粮食生产的负面效应和粮食生产灰水足迹,可为其他粮食主产区农业可持续发展及制定农业产业政策提供参考。     

水稻保护性耕作模式“水田茬地免耕少耕”的效益分析

根据广东省农业面源污染治理项目关于种植业中的化肥、农药减量防控示范工程和农业保护性耕作的有关要求,于2018年在广东省清远市连州市西岸镇冲口蔬菜生产农民专业合作社进行早(中)稻保护性耕作模式水田茬地免耕、少耕同步深施基肥机插秧技术示范试验。试验结果表明,引进先进的机械配套设备用于水田茬地免耕、少耕同步深施基肥机插秧技术,能够有效提高化肥、农药利用率,减少化肥、农药使用量,减少农业面源污染,减少劳动用工和生产成本,促进水稻生产绿色高效。     

基于LMDI方法的恩施州化肥施用驱动因素分析

近年来由化肥过量施用造成的农业面源污染问题日益突出,加大农业污染治理力度,探究化肥施用强度的变化特征及其驱动因素,可以为减肥增效、推进农业绿色发展提供科学依据。研究以湖北省恩施州为实证对象,在研究其时空变化特征的基础上,通过综合考虑化肥施用效率、区域施肥结构、农业技术发展、复种指数、耕地面积等因素,构建对数平均迪氏指数法(Logarithmic Mean Divisia Index,LMDI)指数分解模型与面板回归模型。结果表明:1)近10 a来,恩施州化肥施用量呈现先增长后平稳再下降趋势,氮、磷、钾比例由2007年的1∶0.30∶0.15增长到2018年的1∶0.39∶0.26,逐渐向世界发达国家水平靠拢。截至2018年,恩施州化肥施用共有1个低强度区、2个中强度区、4个高强度区、1个超高强度区。2)2016年恩施州化肥施用量比上年减少7 008.32 t,2017年减少3 341.65 t,2018年减少15 059.81 t。其中化肥施用效率提高使化肥施用量在近3 a分别减少11 159.44、14 051.95和8 644.21 t,复种指数使化肥施用量在2017年减少5 229.76 t,2018年减少30 889.98 t,2016年效应数值不大。农业技术发展使化肥施用量增加了4 089.71、11 136.97和12 215.74t,区域结构的驱动效应不大。3)经济发展、农业结构调整、金融规模扩大对化肥施用量起正向作用,人口数量增加对化肥施用量起负向作用。建议转变农业技术发展方向由增产至高效增产,合理调整有效耕地面积,进一步提升化肥施用效率,培育生态安全绿色发展模式,促进化肥减量增效。     

四川省化肥面源污染环境风险评估及趋势模拟

系统评估农业化肥过量施用对生态环境的潜在风险危害,可为制定下一阶段农业面源污染防治工作提供决策参考。本文基于2000—2015年化肥施用量和耕地面积等数据,采用化肥施用环境风险评价模型,探讨了四川省化肥施用及环境风险的时空变化特征;通过设置延续现状和政策干预两种情景,模拟了四川省2016—2018年执行化肥施用零增长行动期间的化肥施用环境风险变化趋势。结果表明:四川省化肥施用总量在2000—2015年间整体呈增加态势,2015年四川省化肥施用强度为376.63 kg?hm~(-2),超过了中国生态市建设的化肥施用强度250 kg?hm~(-2)的生态标准。四川省2015年化肥施用的环境总风险指数为0.69,处于中等风险程度;氮、磷、钾肥施用的环境风险指数分别为0.69、0.72和0.46。钾肥施用的环境风险处于安全状态;氮肥和磷肥的环境风险区域集中在成都平原经济区、川南经济区和川东北经济区,呈现中等风险程度并有成片聚集特征。延续现状情景下,四川各地区化肥施用环境总风险指数处于0.48~0.69,存在不同程度增加。与延续现状情景相比,四川在完成化肥施用零增长行动的任务规划前提下,2016—2018年政策干预情景的化肥施用环境总风险呈现下降趋势。以关键时间点2018年来看,川东北经济区化肥施用环境总风险下降最明显,较2015年执行化肥零增长时下降4.24%。但四川省多数地区仍存在不同程度的化肥施用环境污染风险,当地生态环境具有明显压力。建议省级部门加强对成都平原经济区和川东北经济区的化肥面源污染监测,根据各地区差异特征制定不同施肥方案,协调粮食增产和生态保护间的关系。     

乐清湾养殖区表层沉积物重金属含量分布及污染评价

以乐清湾养殖区表层沉积物中Cu等7种重金属为研究对象,在野外实地调查取样及重金属含量测定基础上,以中国海洋沉积物质量(GB18668-2002)一类标准值为评价标准,利用单因子指数评价法、内梅罗综合污染指数评价法和沉积物质量基准(sediment quality guideline,SQG),对单项重金属污染程度、多种重金属综合污染效应及重金属潜在生物毒性风险进行了评价。结果表明,乐清湾养殖区表层沉积物重金属含量分布差别明显,所有采样点表层沉积物均受到Cu污染,局部受到严重Hg污染;大部分养殖区表层沉积物重金属污染达到轻度等级;所有养殖区均可发生由Cu偶尔引起的不利生物毒性效应,个别养殖区可发生由Hg频繁引发的不利生物毒性效应。乐清湾养殖区表层沉积物重金属主要来源于沿岸电镀企业废水排放及农药化肥造成的农业面源污染。开展Cu、Hg污染治理修复是实现乐清湾养殖环境安全的迫切需求。     

赤水河流域昭通段化肥减量增效工作现状及建议

赤水河流域(昭通段)是赤水河流域的源头,开展化肥减量增效工作,推广科学施肥可提高肥料利用率,协调作物养分平衡,减轻化肥面源污染,对赤水河流域保护及农业可持续发展具有很好的促进作用.近年来,赤水河流域(昭通段)化肥减量化工作成效明显.2020年,化肥用量(折纯)18413 t,比上年减少3％;化肥利用率达40％以上;完成...     

重庆市渝北区农业面源污染河长制综合防控对策

据调查测算,2016年渝北区规模以上畜禽养殖产生粪污24.94万吨,其中有部分畜禽粪污未经处理直接排放,或露天堆放,经雨水冲淋四处流失,排入沟渠、农田、河(库)等,造成水质恶化,污染水环境。在几种水体养殖中以池塘养殖污染最重。2014—2016年全区化肥平均单位面积施用量为591.3 kg·hm~(-2)(折纯量),是国际公认的化肥施用安全上限的2.63倍。2014—2016年全区平均农药施用总量为93 t·a~(-1),单位施用药量为0.26 kg·667 m~(-2)。分析各类污染的成因。提出农业面源污染防控建议对策:层层落实防控工作,抓好宣传发动,加大政策资金扶持力度,制订各种制度措施整治和防控畜禽养殖污染、水产养殖污染、化肥污染、农药污染。