保障粮食供应是我国经济社会可持续发展的重要基础，化肥施用是粮食增产的重要环节^[1]。但随着我国人口政策的改变，近年来我国人口增长速度有所加快，人均耕地面积随之减少，如何合理施用化肥以保证农作物产量成为了我国亟待解决的问题^[2]。根据联合国粮食及农业组织（Food and Agriculture Organization of the United Nations，FAO）和中国国家统计局提供的化肥施用量数据，2013年全球化肥施用总量为1.67×10⁸ t，中国农业化肥施用当量为5.91×10⁷ t，占全球化肥使用总量的35.48%。同时，单位面积化肥施用量超过世界平均水平三倍多，我国农田化肥过量施用面临的环境风险是其他国家无法相比的^[3-4]。因此，在我国进行化肥施用调研及时空尺度的分析有益于环境管理与污染控制。

近年来，由施肥不当或过量施肥带来的环境污染问题越来越突出。农田氮磷随地表径流进入受纳水体，或随地下淋溶进入土壤，引起水体富营养化和土壤污染等农业面源污染问题^[5]。而地下淋溶对于南方地区影响较小，本研究主要针对地表径流中的农田氮磷流失进行估算。随着工业污染的治理，农业氮磷流失成为水环境污染的主要来源^[6]。农田氮磷流失包括地表径流和地下淋溶两个主要方面，其中地表径流为氮磷流失的主要途径，是农田化肥进入周围水体最直接、最快速的方式^[7]，但由于其随机性强、周期长等特点，治理起来较为困难，因此控制农田氮磷流失是控制农业面源污染的重要途径。目前，我国氮肥利用率仅为30%~35%，磷肥利用率不足20%，远低于发达国家水平^[8]。为此，2016年中央一号文件提出了“加大农业面源污染防治力度，实施化肥农药零增长行动”的目标，对农田化肥调研以及对地表径流氮磷流失的估算与评估有助于实现对农田面源污染的风险防范。

本文以长江中下游所处的六省一市为研究区，以市为单位，开展农业氮磷流失（地表径流流失）估算方法和风险区域识别研究。该区域是我国重要的农业生产基地，但近年来亦产生较多的农业面源问题^[9]。由于化肥的不合理利用，导致富营养问题频发，水环境恶化，水质降低^[10]。长江中下游区域辽阔，因此本文在估算流失量的同时，综合冲刷、污染和河道三个过程进行合理地污染风险评估，以对研究区的地表氮磷污染进行重点区域识别。

1 材料与方法 1.1 试验区概况

研究区位于我国东南部，地跨湖南、湖北、江西、安徽、江苏、浙江、上海六省一市（图 1）。该区域大部分处于长江中下游平原，总面积约为9×10⁶ km², 海拔约为5~100 m，年均温度约为14~18 ℃，年降水量约为1000~1500 mm。研究区属于亚热带季风气候，以冲积平原为主，具有较为优越的农业气候条件^[11]。同时长江中下游地区是我国重要的粮、油、棉生产基地，亦为中国水资源最为丰富的地区，主要河流有长江、汉江、赣江等^[12]。大部分城市也是我国重要的工农业基地，是经济和科技文化发达地区，其中湖南、湖北、江苏省为我国的农业大省，农业生产活动较多^[13]。双季稻是该区域的主要种植模式，水田亦是造成化肥流失的主要区域。

1.2 数据资料

本研究的基础数据资料主要包括市级行政单元划分、农作物播种面积、市级化肥施用数据、省级化肥施用数据、流失系数以及确定流失系数所需的相关数据。主要数据及来源见表 1。

1.3 方法介绍 1.3.1 化肥施用量及施用强度的估算方法

本研究从中国统计年鉴中获取了2000—2015年六省一市农田化肥的施用数据（市级）。有些年份或地区数据存在缺失，为了提高精确性和可利用性，选取了数据相对完善、年份较近的2013年数据进行化肥施用的空间分布特征分析。同时，对缺失数据进行估算，估算方法分为以下两种：

一是只存在化肥施用总量，无氮肥、磷肥和复合肥施用数据的市级行政区，采取以下公式进行估算：

二是部分2013年统计数据缺失的市级行政区，取相邻年份（2012、2014）的化肥施用数据的平均值进行替代。

农田氮磷肥的总施用量包括单一元素肥料和复合肥中氮/磷有效成分之和。复合肥中氮磷钾有效成分采用通用型化肥比例15:15:15进行计算^[20]，即：

化肥施用强度表示为单位市级行政区面积的年化肥施用量（kg·hm^-2·a^-1），相对于化肥施用总量更能直观地反映一个地区的化肥施用状况。计算公式为：

1.3.2 地表流失量及流失率的计算

本研究以农业部出版的《全国农业面源污染流失系数手册》中的农田地表径流流失系数为基础，对研究区80个市级行政区的农田氮磷流失系数进行查算。根据5个要素对各田块的流失系数进行划分，即“农区-地形-梯田/非梯田-旱地/水田-种植类型”（例如：模式47为“南方山地丘陵区-陡坡地-梯田-旱地-露地蔬菜”），获得的氮磷流失系数分别为1.5%、0.85%。根据中国科学院资源环境科学数据中心获取到的数据，对研究区各种要素进行划分，从而查询《全国农业面源污染流失系数手册》获得各个模式的地表氮磷流失系数。并利用公式（4）将各个模式升尺度，获得各市级行政区的氮磷流失系数μ：

式中：β_i为某一田块下某农作物种植面积占市级行政区农作物总播种总面积的比例，%；μ_i为该田块的氮磷流失系数。

其中《全国农业面源污染流失系数手册》定义流失系数公式为：

式中：RL为常规处理下，即正常施肥状态下氮/磷肥流失量，×10⁵ t；CK为对照处理氮/磷肥的流失量，×10⁵ t；FA为氮/磷肥的施用量，×10⁵ t。因此，氮磷肥的流失量可由流失系数与氮磷肥施用量相乘获得。

1.3.3 面源污染风险评价方法的构建

本文按照氮磷指数法的思路，构建农田面源污染风险方法评价，该方法是综合考虑了影响农业面源污染的主要因子，评价流域内不同地区发生氮磷流失危险性高低的一种方法^[21]。该方法考虑三种影响氮磷流失过程的因子，即主要考虑冲刷过程（坡面迁移能力）、污染过程（污染物多少）和河道过程（造成水污染的可能性）。其中，冲刷过程主要考虑降雨因子，污染过程主要考虑氮磷流失量，河道过程主要考虑河网密度。降雨数据以研究区的年平均降雨量数据为准，氮磷流失量以本文计算为准，而河网密度则以水系数据为基础经简化计算获得（将河网长度除以行政区面积）。

面源农田污染风险等级构建如下：首先，市级行政区的农田氮磷流失量、河网密度、降雨等因素的数值进行归一化处理，见公式（6）和公式（7）；其次，将各因子归一化的值相乘，计算氮磷潜在危险性指数，为了简化计算过程，本研究将各影响因子赋予相同的权重；最后，由高到低进行排序，并将各行政区的风险值分为若干等级，在此基础上识别面源污染风险的重点区域。

式中：x_i为各个因子的实际值；x_min为80个市级行政区中的最小值；x_max为80个市级行政区中的最大值；x_i^*为各个因子归一化处理后所得的值。最终得到相应的污染分级方案如表 2所示。

2 结果与讨论 2.1 化肥施用的时空分布特征分析 2.1.1 化肥施用的时间变化特征

长江中下游地区六省一市施肥总量的时间变化趋势如图 2所示。自1996年至2015年20年内，施肥总量波动不大，总体呈现上升趋势，施肥总量范围在1.233×10⁸~1.513×10⁸ t之间，约占全国化肥施用总量的30%。由此说明，研究区所处的六省一市一直为我国的农业中心。自2011年起，施肥总量达到1.5×10⁸ t，2015年再次下降至1.4×10⁸ t，这与2014年我国提出的适时开征化肥税引导农民减少化肥的过度施用有关^[22]。其中，对于各类化肥施用来说，氮、磷、钾肥施用量基本保持不变，波动较小，氮、磷、钾、复合肥所占比例的均值分别为47.97%、16.07%、9.67%、26.29%。可以看出，氮肥施用量最高，钾肥施用量最低，而复合肥施用呈现明显上升趋势。在农田化肥所含有的三种元素中，氮、磷仍是农田化肥所含的主要元素，也是造成农业面源污染的潜在关键元素，因此氮磷流失是本研究所要关注的重点内容^[23]。同时，粮食产量自2003年至2012年呈现上升趋势，至2012年粮食产量已达1.508×10⁹ t，而当年我国粮食总产量为5.896×10⁹ t，实现粮食产量“九连增”，而研究区粮食总产量约为全国的25.58%^[24]。同样的，2012年至2015年研究区的粮食产量也呈现上升趋势，而1996年至2003年粮食产量呈现下降趋势。由此可以看出，化肥施用量的增加在一定程度上促进粮食增产，但化肥施用和粮食产量没有呈现完全的正相关，合理的化肥施用更能保障粮食产量和环境管理。

2.1.2 化肥施用的空间分布特征

长江中下游地区幅员辽阔，各个省市的社会经济和自然环境状况存在明显差异。因此，本研究利用地理信息系统进行空间分布特征分析，如图 3所示。为了更直观地显示市级行政区的氮磷施用状况，采用化肥施用强度作为评估指标。对于氮肥施用强度来说，湖南省是施用强度最大的省份，东部地区都为高氮肥施用强度区，西部地区处于较高氮肥施用强度地区；相比之下，江西省和上海市的大部分区域氮肥施用强度都处于较低水平；磷肥施用强度与氮肥施用强度呈现出明显差异。湖南省东部和湖北省的中部属于磷肥流失强度较高的区域，对这几个市级行政区依次排序为湘潭市>益阳市>常德市>宜昌市>岳阳市>襄樊市；而浙江省、江西省和上海市的大部分地区磷肥施用强度较低。综上所述，湖南省为六省一市中氮肥施用强度最大的省份，而磷肥施用强度较高的为湖南湖北两省的部分区域。因此，对于肥料施用控制管理，应具体位置具体分析，同一个省份不同市级单位的化肥施用量仍具有一定的差异。同时，上海市的氮磷施用强度一直处于较低水平，这与上海城市化水平较高、农田面积比例较小有关。

2.2 农田化肥氮磷流失特征

不同的流失率将直接导致氮磷流失量的不同，农田化肥氮磷流失量的空间分布如图 4所示。总体来说，湖南省是氮磷流失的主要区域。对于氮肥流失量，湖南省处于高流失量的市级行政区排序为：益阳市（7.66×10³ t）>岳阳市（7.31×10³ t）>衡阳市（7.16×10³ t）>永州市（6.05×10³ t）>长沙市（5.37×10³ t）>邵阳市（5.26×10³ t），其次湖南省的郴州市、株洲市、湘潭市以及湖北省的襄樊市和江苏省的徐州市、盐城市属于次高氮流失区；相比之下，江西省和浙江省的大部分区域和上海市都属于氮肥流失较少的区域。总体上来看，研究区氮肥流失量处于0.521×10³~1.214×10³ t范围的居多。对于磷肥流失量，湖南省具有较高流失量的市级行政区排序为：常德市（2.69×10³ t）>益阳市（1.97×10³ t）>岳阳市（1.81×10³ t）>衡阳市（1.73×10³ t），湖南省的其他市级行政区磷流失量也属于较高水平；同时，其他五省一市的磷流失量基本小于1000 t，尤其是浙江省的全部区域磷的流失量小于200 t。总体上来看，研究区磷肥流失量处于0.154×10³~0.381×10³ t范围的居多。

氮磷流失强度分布如图 5所示。对于氮肥施用强度来说，研究区流失强度的平均值为2.41 kg·hm^-2·a^-1；而湖南省大部分地区都属于氮肥流失强度较高地区，氮肥施用强度远高于平均水平，对这几个市级行政区依次排序为常德市>岳阳市>长沙市>株洲市>衡阳市>永州市>郴州市>娄底市>邵阳市，湖北省的襄樊市和宜昌市、江苏省的北部地区氮肥流失强度也处于较高水平；而江西省大部分地区氮肥流失强度都处于较低水平。对于磷肥施用强度，研究区流失强度的平均值为0.61 kg·hm^-2·a^-1。与氮肥流失强度相同，湖南省大部分地区都属于磷肥流失强度较高地区，其中湖南省的东北部属于流失强度最高的地区，对这几个市级行政区依次排序为湘潭市>益阳市>常德市>岳阳市>长沙市>株洲市；而安徽省、江西省和江苏省的大部分地区磷肥流失强度较低。综上所述，湖南省为氮肥、磷肥流失强度最大的一个省份，同时氮肥、磷肥的流失强度在空间分布上稍有不同。这与不同地区的种植模式有关，也与氮磷元素在自然界中的形态有关，氮以溶解态为主，磷以吸附态为主，氮肥与磷肥相比更易流失^[25]。

2.3 农田面源污染风险评估

利用1.3.3的方法对长江中下游地区六省一市开展农田面源污染重点区域识别，风险分级结果如图 6所示。总体上来看，大部分高流失强度区域亦是高污染风险区，但有些区域，例如江苏省的苏州市，磷肥流失强度为0.351 kg·hm^-2·a^-1，却处于较高流失风险区，这与苏州市临近长江，河网密度较大有关。同时，氮磷地表流失污染存在着一定的空间差异，但湖南省中东部地区同时处于氮磷流失的高风险区，上海市和安徽省的北部地区同时处于氮磷流失的低风险区。对于农田氮污染风险分区，湖南省的11个市级行政单元、安徽省的铜陵市和芜湖市、江苏省的淮安市和宿迁市、湖北省的鄂州市属于高污染风险区；湖南省的其他地区均属于较高污染风险区，其他五省也均有分布；相比之下，上海市、安徽省北部和江西省的大部分地区都处于氮的低污染风险区；同时，高污染风险区约占研究区面积的20.27%，湖南省高污染风险区约占总面积的75.19%。对于磷地表径流污染风险分区，湖南省的11个市级行政单元、安徽省的安庆市和铜陵市、湖北省的鄂州市、江西省的南昌市属于高污染风险区；湖南省的其他地区均属于较高污染风险区，除浙江省外其他四省也均有分布；相比之下，上海市、安徽省北部和南部、江西省东部部分区域处于氮的低污染风险区；同时，高污染风险区约占研究区面积的20.07%，湖南省高污染风险区与氮肥高污染风险区相同，约占总面积的75.19%。综上所述，湖南省为氮磷流失的高风险区域，上海市为氮磷流失的低风险区域，其他五省以中、低污染风险为主。因此，长江中下游地区所处的六省一市污染防治宜“以防为主，以治为辅”，但对于湖南省的高污染风险区域应加强控制措施^[12]。

3 结论

本研究以长江中下游所处的六省一市为研究区，调研了区域内市级行政区化肥数据及农作物播种面积，以农田氮磷流失系数为基础计算出对应的流失量，并利用地理信息系统进行了降雨量和河网密度的空间分析，确定了农田污染风险等级。该方法可直观地比较各市级行政区的氮磷流失量，识别农田污染高风险区，为农田面源防控和水环境治理提供理论基础。

案例研究结果表明，近20年来，研究区化肥施用量呈现上升趋势。湖南省是氮肥施用强度和流失强度最大的省份。湖南、湖北则为磷肥施用强度较高的区域，该区域氮磷流失强度的平均值分别为2.41、0.61 kg·hm^-2·a^-1。从风险分级的角度，农业化肥对周边水体的高污染风险区为湖南省，上海市为低风险区域，其他五省则以中、低污染风险为主。