江西省农田投入品输入时空演变特征及与种植业结构调整的关系.pdf

生态与农村环境学报，（）：收稿日期：基金项目：国家自然科学基金（）；江西省教育厅科技项目（）通信作者：温小珊，冯雅雯，杜超，等江西省农田投入品输入时空演变特征及与种植业结构调整的关系生态与农村环境学报，（）：，（）：江西省农田投入品输入时空演变特征及与种植业结构调整的关系温小珊，冯雅雯，杜 超，黄秀珍，王宵君，谢金招（赣南师范大学地理与环境工程学院，江西赣州；江西省城市固废低碳循环技术重点实验室，江西 赣州；赣州市果业发展中心，江西 赣州）摘要：厘清近 年间江西省农田投入品输入的时空演变特征，并探究其与种植业结构调整的关系，可为江西省农业生产提质增效和农业面源污染防治提供理论依据。采用重心模型和结构均衡度方法，基于 探究了江西省 年农田化肥、农药和农用薄膜输入的时空演变特征，并通过线性回归分析探讨了种植业结构调整与农田投入品输入量的关系。结果表明：（）年江西省化肥、农药和农用薄膜输入量整体呈先增后减的变化趋势，自 年实施了一系列农药化肥减量增效政策之后，江西省农田化肥和农药输入量明显降低，年化肥、农药输入量比 年分别下降 和 ，但是近年来农用薄膜输入量的减缓趋势并不明显；（）江西省化肥输入强度存在明显的空间聚集现象，总体呈现沿鄱阳湖流域、抚河流域和赣江流域聚集分布的趋势；农药输入程度在赣东北地区普遍较高；农用薄膜输入量在全省普遍增加，其中，赣东和赣南农用薄膜输入量增加尤其明显；（）江西省种植业结构调整表现为果园种植面积占比不断增大，从 年的 万 增长到 年的 万，占比由 增长至 。与化肥和农药输入量相比，果园种植面积扩张与农用薄膜输入量呈更强的相关性，表明江西省果园种植面积扩张与农用薄膜输入量增加关系密切。关键词：面源污染；演变特征；空间分布；种植业结构调整；江西省中图分类号：文献标志码：文章编号：（），（，；，；，）：，：（），第 期 温小珊等：江西省农田投入品输入时空演变特征及与种植业结构调整的关系 ；（），；（），：；广义的农业面源污染来源包括种植业（农田）、畜禽养殖业、水产养殖业和农村生活污水 个方面，其中，农田面源污染是重要的污染来源，主要是由化肥、农药和农用薄膜等农田投入品的过量输入产生。近年来，农业生产过程导致的农田面源污染问题依旧严峻，其结果是增加水体污染风险，影响耕地资源质量，并威胁农产品质量安全。江西省是我国重要的粮食主产区和全国首批 个生态文明试验区之一，长期的过量施肥和高度集约化的种植模式，导致该地区农田面源污染形式严峻，衍生出一系列生态环境问题。因此，系统探究江西省近 年来农田化肥、农药和农用薄膜输入的时空演变特征，精准刻画其输入轨迹，对于农业生产提质增效和农业面源污染防治具有重要指导意义。近年来，已有许多学者对全国及各省（区、市）农田投入品输入进行研究，得出一些结论。马兆嵘等分析了全国的农膜输入现状，发现我国农膜输入量巨大且在不断增加，由 年的 万 增长到 年的 万，增长 倍。葛小君等分析了广东省农业面源污染负荷时空变化特征，发 现 耕 地 化 肥 和 农 药 输 入 强 度 分 别 上 升 和 。刘存庆等基于 探究了菏泽市农业面源污染时空分布，发现 年之后菏泽市化肥、农药输入量都在稳步降低，但化肥输入强度仍高于国际化肥预警值（）。傅春等研究发现磷是江西省农田生产的主要污染物，年化肥输入量中占比最大的复合肥是主要污染源。熊昭昭等研究发现，年江西省农村地区 和 污染负荷分别为 和 ，其中，种植业带来的污染负荷贡献率达到 。类似地，针对某地区农业面源污染输入与输出负荷评价的研究虽时有报道，但是结合重心迁移模型和结构均衡度方法，系统探究污染物输入总量和输入强度在地理空间上移动规律和长时间线上演变特征的研究尚鲜见。此外，影响农田投入品输入的因素众多，其中，种植业结构调整的影响在以往的研究中没有受到足够重视。“耕地非粮化”现象在全国各地农业生产活动中有不同程度的发生，随着赣南脐橙和南丰蜜桔等江西特色柑橘产业的兴起，水果种植业在江西农业种植中的占比逐年加大，关于种植业结构调整对江西省农田投入品输入的影响，目前尚鲜见报道。基于此，为了全面把握江西省农田投入品输入特征及影响因素，以江西省 个地级市为研究单元，基于化肥、农药和农用薄膜输入量统计数据，采用重心模型和结构均衡度等模型方法，系统分析江西省近 年来农田投入品输入量和输入强度的时空演变规律，并着重探讨种植业结构调整对江西省农田投入品输入的影响，以期为江西省农业面源污染防治、种植业结构优化和各地区的协调发展提供科学依据。研究区域与数据来源 研究区概况江西省位于我国东南部，长江中下游南岸，处于北纬 、东经 之间，总面积达 万，设有 个地级市（图）。江西省东西南三面环山，中部为盆地丘陵，北部为鄱阳湖平原，属亚热带季风湿润气候区，年均气温为 ，年均降水量为 。江西省是我国重要的粮食生产基地之一，粮食作物面积总计达 万（年）。年，全省化肥输入总量从 万 增加到 万，农药输入总量从 万 增加到 万，农用薄膜输入总量从 万 增加到 万。生 态 与 农 村 环 境 学 报第 卷图 江西省行政区划及主要水系图 数据来源统计数据主要包括江西省及其 个地级市化肥输入折纯量，农药和农用薄膜输入量和主要种植结构面积数据主要来源于 年江西统计年鉴。研究方法 重心模型重心可反映全省各市化肥、农药和农用薄膜输入强度的总体分布情况。基于化肥、农药和农用薄膜输入的分散情况，在具体研究中几乎不可能测定研究区域内所有化肥、农药和农用薄膜输入的实际位置。因此，假设江西省各市区域为 个均质平面，各市级区域的输入量或输入强度重心集中于各市的几何重心。重心坐标计算公式为 ，（）。（）式（）（）中，为参与评价的样本数量；（，）为第 个市的重心坐标；为该单元属性。面源污染输入结构均衡度面源污染输入结构均衡度（）是基于信息熵函数构建的用于反映面源污染输入空间特征的量值，其值处于 之间。可反映数据的离散程度，越小表示指标数据的离散程度越大，计算公式为 ，（）。（）式（）（）中，为 省 年面源污染物输入强度；为参与评价的样本数量。自然断点法自然段点法是将数据集中不连续的地方作为分级依据对数据集合进行分级的方法。在使用 技术处理空间信息时，一般采用自然断点法对空间分布等级进行划分。分别基于化肥、农药和农用薄膜输入强度将面源污染划分为低、较低、中等、较高和高 个等级（表）。表 面源污染程度等级划分 等级化肥输入强度（）农药输入强度（）农用薄膜输入强度（）低，较低（，（，（，中等（，（，（，较高（，（，（，高 结果与分析 江西省农田投入品输入时序演变特征 江西省农田投入品输入总量的时序分布特征如图 所示，年间前期江西省化肥、农药和农用薄膜输入量均有波动，但整体上均呈现先增后减的变化趋势，增长率分别为 、和 。化肥输入总量变化率最大的为萍乡市，增长 ；宜春市化肥输入总量变化率最小，降低 。在 年之前，各市化肥输入总量呈现增长趋势，其中，赣州、吉安和抚州持续增长。各市化肥输入总量增长率总体呈下降趋势，景德镇、萍乡、九江、赣州、抚州和上饶近 年大多呈负增长；年前后，南昌、景德镇、萍乡、九江、宜春和上饶都呈明显负增长；赣州增长率在 年突降，降至近 年最低值（）。总体来看，江西省农田化肥输入量在 年之后下降明显，年比 年下降 。农药输入总量变化率最大的为鹰潭市，增长 第 期 温小珊等：江西省农田投入品输入时空演变特征及与种植业结构调整的关系；农药输入总量变化率最小的为南昌市，降低 。总体来看，各市农药输入总量呈现先增后减趋势，其中，南昌、景德镇、萍乡和新余农药输入量在 年前后大大增加，而近年来呈下降趋势；九江、抚州和上饶在 年之前呈背斜式增加，之后呈下降趋势；赣州、吉安农药输入量在 年出现明显的断崖式减量。各市农药输入总量增长率总体而言呈下降趋势。其中，南昌市增长率在 年迅速下降至负值，至次年又快速上升；景德镇、吉安、抚州增长率则在 年增长迅速；萍乡、九江、新余、鹰潭、抚州和上饶增长率则在 年前后上升显著。总体来看，江西省农药输入量与化肥输入量类似，年比 年明显下降 。图 江西省各市化肥、农药和农用薄膜输入总量 ，薄膜输入总量变化率最大的为鹰潭市，增长 ；薄膜输入总量变化率最小的为萍乡市，增加 。总体来看，各市薄膜输入总量呈现增长趋势，其中，景德镇薄膜输入量在 年出现大幅度减量，呈现一个凹槽；新余、鹰潭、赣州、吉安、宜春和上饶呈阶梯式增长，中间偶有高峰低谷：新余在 年出现一个增长高峰，年出现增长次高峰；鹰潭在 和 年前后出现两次增长高峰；年赣州薄膜输入量出现低谷，而吉安薄膜输入量出现高峰，各市薄膜输入总量增长率总体而言呈上升趋势。江西省农田投入品输入强度的时序分布特征 年江西省化肥、农药和农用薄膜输入强度（输入总量与农作物播种面积之比）见图。总 体 来 看，年 化 肥 输 入 强 度 从 增长到 ，农药输入强度从 增长到 ，农用薄膜输入强度从 增长到 。江西省化肥和农药输入强度均在 年之前大体上呈增长趋势，之后呈波动下降趋势；而江西省农用薄膜输入强度在 年达到小高峰后有所降低，但从 年至近几年整体上呈增长趋势。图 江西省各市化肥输入强度 图 江西省各市农药输入强度 生 态 与 农 村 环 境 学 报第 卷图 江西省各市农用薄膜输入强度 总体来看，各市化肥输入强度呈现先增后减趋势，其中，新余和鹰潭化肥输入强度波动较大，年之后输入强度均明显降低。各市农药输入强度总体上呈先增后减趋势，年之前，江西省各市农药输入强度呈波动增加趋势，年各市农药输入强度明显增加，年之后农药输入强度波动降低，而上饶、抚州降低程度不大，仍处于高输入强度。江西省农用薄膜输入强度在 年有段小高峰，除抚州、南昌外，各市薄膜输入强度总体而言呈上升趋势。江西省农田投入品输入强度的空间演变特征 江西省农田投入品输入强度的空间分布特征将 年设为基期，采用等距法每 年抽取相应年份，并将末年 年也纳入研究，最终选取、和 年农田投入品输入强度，运用 软件制图展示空间分布特征。如图 所示，近 年来江西省大部分城市化肥输入强度呈先增后减的变化趋势，年江西省只有萍乡属于较高污染等级，和 年江西省化肥输入强度的空间特征表现为沿鄱阳湖流域、抚河流域和赣江流域呈规律性聚集趋势，年江西省绝大部分城市农药输入强度降低或不变。具体来看，萍乡、南昌一直保持着较高的化肥输入强度。如图 所示，和 年江西省农药输入强度空间分异不明显，各市农药污染程度大多不高于中等水平；和 年均呈现赣东北地区农药输入强度高的态势，年江西省农药污染处在较高程度及以上的市高达 个，其中，上饶、南昌、鹰潭、抚州、新余和萍乡达到高农药污染程度，年全省除了景德镇和宜春外，其余地市农药污染均得到明显缓解。如图 所示，总体来看，年江西省各市农用薄膜输入强度总体呈逐渐增大的变化趋势，其中，赣东南地区薄膜污染尤其显著，其余地区农膜污染程度持续处于中等水平以下。具体来看，赣州市薄膜污染程度显著增大，从 年低污染等级上升至 年的高污染等级。江西省农田投入品输入强度的重心迁移 年江西省农田投入品输入强度重心迁移图见图。由图 可知，年化肥输入强度重心落在 、之间，均位于宜春境内偏东南部；农药输入强度重心在 、之间，大部分位于宜春市东南部，少部分位于吉安市东北部；农用薄膜输入强度重心在 、之间，大部分位于宜春市内，部分位于吉安市东北部。图 江西省化肥输入强度空间分布 第 期 温小珊等：江西省农田投入品输入时空演变特征及与种植业结构调整的关系图 江西省农药输入强度空间分布 图 江西省农用薄膜输入强度空间分布 图 面源污染输入强度重心随年份的迁移 具体来看，年化肥输入重心累计向西迁移 ，向南迁移 ，整体呈现由东北向西南迁移的态势，期间有 次向南移动，移动频次为，向西移动 次，移动频次为，整体表现为一段明显的波浪式曲折运动；农药输入强度重心累计向东迁移 ，向北迁移 ，整体呈现由西南向东北迁移的态势。分阶段来看，江西省农药输入强度的重心迁移轨迹表现出明显的南北走向波动，年之前总体向南迁移，之后重心基本向北迁移；农用薄膜输入强度重心累计向西迁移 ，向南迁移 ，整体呈现由东北向西南迁移的态势。农用薄膜输入强度重心东西迁移不断变化交替，表现出较明显的波动性特征，阶段性明显。年，农用薄膜输入强度重心迁移比较迅速，范围较大，但近年来，重心迁移范围和跨度明显减小。江西省农田投入品输入强度的区域差异性变化趋势如图 所示，年江西省化肥输入强度的结构均衡度一直处于较高状态，但在 生 态 与 农 村 环 境 学 报第 卷年，江西省化肥输入强度的结构均衡度达到最低值（），年之后，市域间施肥强度的差异呈现减小趋势。年江西省农药输入强度的结构均衡度波动比化肥输入强度结构均衡度大，整体呈现“”型变化趋势。年江西省农药输入强度的结构均衡度达到最高值（），表明这一年各市农药输入强度的差异性最小；年达到最低值（），表明市域间农药输入强度差异性大。年江西省农用薄膜输入强度结构均衡度呈先增后减的变化趋势，整体来说结构均衡度跨度较大。年农用薄膜输入强度的结构均衡度急剧上升，比上一年上升 ，之后薄膜输入强度的结构均衡度一直处于波动状态。图 面源污染输入强度的结构均衡度 江西省农田投入品输入与种植业结构调整的关系 随着社会经济的发展，人们的消费需求不断升级，种植业结构也在随之调整。种植业结构直接决定化肥、农药和农用薄膜输入量，是影响农田面源污染的重要因素之一。江西省是全国粮食主产区之一，农业种植主要分为粮食作物（包括谷物、豆类、薯类等）、经济作物（包括油料作物、棉花、麻类、甘蔗、烟叶、中药材、蔬菜等）和水果（脐橙、柑橘等）类。基于江西省种植特色，笔者将粮食作物、经济作物和水果 者的种植面积之和作为总量，再分别研究 者种植面积占比变化。如图 所示，江西省种植结构以粮食作物为主体，经济作物占比次之，水果占比最低。年江西省粮食作物种植面积持续处于高占比状态且波动不大，其占比均在 以上；经济作物种植面积明显缩减，从 年的 万 减少到 年的 万，占比由 减少到；水果种植面积明显扩张，从 年的 万 增长到 年的 万，增长率为 ，占比由 增长至 。图 江西省种植业结构面积占比 第 期 温小珊等：江西省农田投入品输入时空演变特征及与种植业结构调整的关系 为验证江西省种植业结构调整与农田投入品输入的关系，将不同种植类型面积分别与化肥、农药和农用薄膜输入量之间进行线性回归分析，结果见图。线性回归分析结果表明粮食作物种植面积与化肥、农药和农用薄膜输入量呈正相关，且相关系数 差异较小，依次为 、和 ；经济作物种植面积与化肥、农药和农用薄膜输入量呈负相关，依次为、和；果园面积与化肥、农药和农用薄膜输入量之间呈显著正相关，依次为 、和 ，其中，农用薄膜输入量与果园面积之间的线性相关程度最高。图 不同种植业面积与农田投入品输入总量的相关分析 讨论与结论 讨论通过对江西省农田投入品输入时空分布特征的分析可知，年江西省农药输入总量出现暴跌，主要是因为该年江西省开展了农药减量提升行动，各市积极响应该项工作。在农业经济发展早期，农户追求经济效益，但又缺乏科学指导，大量投入农业生产要素，从而忽视了农业生态环境的维护，农户在作物种植相对稳定时也会持续加大化肥输入量，以增加化肥输入量替代劳动力的投入。自 年到 年化肥施用量零增长行动方案实施后，化肥、农药和农用薄膜输入量明显得到控制，使得化学农药输入量减少。年江西省化肥、农药和农用薄膜的输入增长率均呈波动下降趋势，且 年之后 者增长率均呈减小趋势，近几年增长率甚至低至负值；这可能是由于生态文明建设不断推进，国务院研究部署了农业农村面源综合治理工作，且人们对绿色食品的需求增加，对食品安全更加关注。和 年江西省化肥输入强度的空间特征表现为沿鄱阳湖流域、抚河流域和赣江流域呈规律性聚集的趋势，究其原因，鄱阳湖流域和抚河流域是江西的水稻主产区，农民会过量施肥以保证水稻等粮食作物产量。因此，九江、南昌和抚州化肥输入强度普遍较高，而赣江流域所处的赣州市柑橘类水果种植业发达，果农为追求高产往往也会过量输入化肥。年江西省绝大部分地市农药输入强度降低或 生 态 与 农 村 环 境 学 报第 卷不变，由此可见江西省对化肥输入零增长行动的执行效果已然显现。具体来看，萍乡、南昌一直保持着较高的化肥输入强度；究其原因，萍乡占地面积小，农用地面积少造成农业生产压力大，农户往往会通过提高生产要素的投入来增产和提高收入，而南昌作为江西省水稻主产区之一，化肥输入量也普遍较高。赣州农用薄膜输入强度由 年的低污染等级升至 年的高污染等级，这是因为赣州市耕地面积大，且脐橙等果园种植业蓬勃发展，对农用薄膜的需求量逐渐增大。通过对江西省农田投入品输入强度空间分布特征分析可知，与 年相比，年江西省化肥输入强度重心向西南方迁移。结合图 发现，年吉安市、赣州市化肥输入强度等级增大，景德镇市化肥输入强度等级减小，导致化肥输入强度重心向西南方迁移，该重心迁移趋势反映出，在深入推进化肥减量增效的背景下，西南方向化肥输入强度呈现明显高于其他区域的增长趋势。江西省农药输入强度重心向东北方向迁移，结合图 发现，年上饶市、景德镇市和鹰潭市农药输入强度等级增大，导致农药输入强度重心向东北方向迁移。江西省农用薄膜输入强度重心向西南方向迁移，结合图 发现，年萍乡市、吉安市和赣州市农用薄膜输入强度等级增大，导致农用薄膜输入强度重心向西南方向迁移。年，江西省开始探索北粳南引，方克明等于 年在景德镇市开展了“籼改粳”晚稻产势与化肥减量增效施用方法试验，试验结果表明籼稻改粳稻很好地实现了化肥减量。柳开楼等在有关江西省不同地市化肥用量分析及化肥用量对农业产值贡献的研究中发现，赣州由于脐橙、油茶等作物种植面积较大，且为了提高种植效益，使得赣州市化肥输入量长期以来居高不下。因此，年景德镇市化肥输入强度等级减小而赣州市化肥输入强度等级增大进而导致江西省化肥输入强度重心向西南方向迁移。近 年来，江西省市域间化肥输入强度的差异性变化并不明显，但 年九江市化肥输入强度骤减，南昌市、新余市化肥输入强度加大，导致江西省化肥输入强度差异性增大。年后，化肥“零增长”行动使得江西省各市施肥强度往更加均衡的方向发展，因此，市域间施肥强度的差异呈现减小趋势。年南昌市、赣州市和吉安市农药输入强度大幅度下降，使得该年江西省市域间农药输入强度差异性明显，呈现个性化特征。年市域间农用薄膜输入强度的差异性急剧增大，可能是由于南昌市农用薄膜输入强度骤增，而萍乡市农用薄膜输入强度骤减，同时这可能也与各市种植结构存在差异有关，再加上市域间发展也存在不均衡问题，使得江西省各市农用薄膜输入强度差异明显。近 年来，江西省粮食作物种植持续处于高占比状态且波动不大，其占比均在 以上，相对比较稳定，而经济作物种植面积明显缩减，但是水果种植面积呈现明显扩张态势，占比由 年的 增长为 年的 。这可能是由于江西省赣南脐橙、南丰蜜桔等特色柑橘产业和其他特色水果享誉全国且经济效益高，刺激了种植户的种植意愿。由此可见，江西省种植业结构调整表现为在保证粮食安全生产的前提下对经济作物和果园种植结构进行优化。粮食作物种植面积与 种农田投入品输入量的相关系数均在 以上，说明粮食作物面积对江西省农田投入品输入量有积极影响，粮食作物面积越大，种农田投入品的输入量就越高。经济作物种植面积与化肥、农药、农用薄膜输入量的相关系数均为负值，分别为、和。相比于化肥和农药输入量，农用薄膜输入量与果园种植面积扩张呈现更强的正相关性，相关系数达到 ，表明江西省农用薄膜输入量增加与果业种植面积扩张关系密切。江西省粮食作物结构调整、高用药农作物对低用药农作物的种植替代，是化肥、农药输入量增长的重要影响因素。此外，随着江西省种植业结构的调整，农田面源污染防治工作要着重关注果园的农田投入品输入，尤其是农用薄膜的输入。不少地区果园采用地膜覆盖技术来减少水土流失，提高果园生产力。农用地膜等塑料制品的输入是土壤中微塑料的重要来源，微塑料在土壤中累积会破坏良好生态环境，降低土壤微生物群落多样性。因此，江西省要严格控制薄膜输入量并加大其回收利用力度，研发和推广可降解农膜，严格执行塑料农膜标准，健全农膜回收利用机制。结论基于、重心模型和结构均衡度等技术手段，对比分析了近 年江西省农田投入品输入的时空分布特征与区域差异，并着重探究了种植业结构调整对农田投入品输入的影响，得出以下结论：（）近 年来，江西省化肥、农药和农用薄膜输入量整体呈现先增后减的变化趋势。自 年实施了一系列农药化肥减量增效政策之后，江西省农田化肥和农药输入量明显降低。与 年相比，年化肥和农药输入量分 别 下 降 和 第 期 温小珊等：江西省农田投入品输入时空演变特征及与种植业结构调整的关系，但是近年来农用薄膜输入量的减缓趋势并不明显。（）江西省化肥输入强度在空间上的特征表现为沿鄱阳湖流域、抚河流域和赣江流域呈聚集性分布；农药输入强度的空间分布呈现赣东北地区普遍较高的态势；而农用薄膜的输入强度普遍增大，其中赣东南地区农用薄膜污染形势较为严峻。自 年以来，江西省化肥、农药输入强度的区域性差异减小，而农用薄膜输入强度呈现强烈的个性化特征。（）江西省种植业结构调整表现为在保障粮食作物播种面积的基础上，果园种植面积占比不断增大。粮食作物和经济作物种植面积对农田投入品输入量分别具有一定的正向和负向影响。与化肥和农药输入量相比，果园种植面积扩张与农用薄膜输入量呈现更强的正相关性，表明江西省果园种植面积扩张与农用薄膜输入量增加关系密切。参考文献：，：，（）：，：，：，：，（）：姜旭海，韩玲，李帆陕西省主要作物灰水足迹时空变化特征研究干旱地区农业研究，（）：，（）：刘亦晴，张建玲比较视角下江西生态文明试验区建设研究：基于福建、江西、贵州三个首批生态文明试验区的比较生态经济，（）：，：，（）：，（）：李淑英，黄秋萍，兰美华，等江西农业面源污染及其控制对策江西农业学报，（）：，（）：马兆嵘，刘有胜，张芊芊，等农用塑料薄膜使用现状与环境污染分析生态毒理学报，（）：，（）：葛小君，黄斌，袁再健，等近 年来广东省农业面源污染负荷时空变化与来源分析环境科学，（）：，（）：刘存庆，何苗，张翠静，等基于 的菏泽市农业面源污染时空分布能源环境保护，（）：，（）：傅春，陈毓迪，刘业忠，等江西省农田灰水足迹时空分析农业环境科学学报，（）：，（）：熊昭昭，王书月，童雨，等江西省农业面源污染时空特征及污染风险分析农业环境科学学报，（）：，（）：，（），（）：，：张世昙肥东县种植业面源污染现状与治理对策研究安徽农学通报，（）：，（）：，邱乐丰，张玲，徐保根，等种植结构非粮化对农田氮磷流失负荷的影响浙江农业学报，（）：，（）：卿娟，张起明，齐述华 江西省农业化肥施用量时 生 态 与 农 村 环 境 学 报第 卷空差异研究 江西科学，（）：，（）：张永树，杨振凯，訾璐，等中国艾滋病空间格局和时空演化分析地球信息科学学报，（）：，（）：朱磊，杨爱民，夏鑫鑫，等基于空间自相关的 年玛纳斯河流域耕地时空特征变化分析中国生态农业学报（中英文），（）：，（）：，：孙兆凯，冯欣，刘洋，等中国化肥施用强度的时空分异特征中国农业资源与区划，（）：，（）：徐新良，陈建洪，张雄一我国农田面源污染时空演变特征分析中国农业大学学报，（）：，（）：陈友权，陈禹涵深入推进农业种植结构调整民主与科学，（）：，（）：刘雪青，黄琦，揭玲，等 年江西省种植业内部结构分析现代农业科技，（）：，（）：王启，张辉，廖桂堂，等四川省主要农业投入品时空变化特征及影响因素生态与农村环境学报，（）：，（）：胡贤辉，张紫怡，王子诚，等洞庭湖平原耕地面源污染时空特征及脱钩效应水土保持通报，（）：，（）：王惠明，陈燕，刘晖，等江西省不同流域种植业面源污染现状分析中国农学通报，（）：，：，（）：黄斌，黄明，彭杨，等赣南脐橙园土壤状况及果实品质相关性试验基层农技推广，（）：，（）：方克明，肖欣，江麟，等“籼改粳”稻产势与化肥减量增效方法试验效果 农学学报，（）：，（）：柳开楼，万长艳，胡丹丹，等江西省不同地市化肥用量分析及对农业产值的贡献湖北农业科学，（）：，（）：，仇相玮，胡继连我国农药使用量增长的驱动因素分解：基于种植结构调整的视角生态与农村环境学报，（）：，：，（）：段四喜，杨泽，李艳兰，等洱海流域农业面源污染研究进展生态与农村环境学报，（）：，：，（）：，：？，（）：张秀玲，鄢紫薇，王峰，等微塑料添加对橘园土壤有机碳矿化的影响环境科学，（）：，（）：作者简介：温小珊（），女，江西赣州人，主要研究方向为南方红壤丘陵区农业面源污染防控。：（责任编辑：李祥敏）