农业生产中，氮肥的投入极大提高了粮食产量，有资料表明因氮肥投入而增产的作物养育的人口数量超过了35亿^[1-2]。同时，氮肥的大量施用也导致了严重的环境问题，如水体富营养化、土壤酸化、N₂O和NH₃排放增加等^[3-6]。山东省是我国小麦和玉米的生产大省，根据2017年山东省统计年鉴，化肥单位面积平均用量达到466.54 kg N·hm⁻²·a^-1，而德州市农户的氮肥平均施用量已经高达591 kg N·hm⁻²·a^-1，远高于发达国家为防止化肥对水体造成污染而设置的225 kg N·hm⁻²·a^-1的安全上限^[7]。第一次全国污染源普查结果显示，山东省种植业源污染物排放占全国农业源比例接近10%。因此，在保证粮食作物产量的前提下，科学施用氮肥已经成为山东省农业相关部门亟待解决的问题。

合理的施氮量是科学施肥最为重要的内容^[8]，而根据最佳氮素盈余量可以推导出合理的施氮量范围^[9-10]。农田氮素盈余是氮素管理中的重要指标，是评价氮素投入生产力、环境影响和土壤肥力变化的最佳指标^[11]。某一生产单位合理的氮素盈余量实质上就是氮素损失后该生产单位可以自身建立起循环且不导致大气、土体和水体污染的氮素含量^[9]。近年来，已有的氮素盈余参考值往往是华北地区、长江中下游地区等大尺度区域下得出的笼统结果^[10]，如在巨晓棠等^[9]文章中，综述了华北地区氮素损失与氮素盈余情况，提出华北地区冬小麦-夏玉米轮作体系下氮素盈余量推荐值为80 kg N·hm⁻²·a^-1。氮素盈余值由于尺度不同、区域不同，其值往往也有很大的差异。因此，针对某一具体的小尺度区域，系统性的研究农田氮素盈余，根据盈余量参考指标体系进行优化施肥，对当地农田氮素优化管理十分重要。另外，氮素损失量是评估氮素扩散到环境中的重要指标^[11]，并结合氮素盈余量，根据二者的关系建立氮素污染指标体系可以有效地评价农田污染情况。但因土壤存在较大的空间变异性，氮素分布和流失受诸多因素影响，找到可以把复杂的因素简单化的评价方法非常关键，欧氏距离法在评价农业可持续发展的许多研究中被应用^[12]，并且评价结果符合实际情况，本文在其基础上进行改进，构建了氮素管理指标评价方法，可以直观地表现德州市氮素管理水平。

结合文献荟萃和长期定位试验，本文研究了山东省德州市冬小麦-夏玉米轮作体系土壤氮素盈余指标体系的构建与评价方法，旨在建立该地区合理的氮盈余指标体系和评价方法，指出德州市氮素管理存在的问题并给出相应的建议，为该地技术的推广应用和政策的制定提供科学依据，也为其他小尺度区域研究提供参考。

1 材料与方法 1.1 研究地概况

本研究以华北地区德州市为研究区域，研究该地区冬小麦-夏玉米轮作体系土壤氮素平衡与循环。德州市位于山东省西北部的黄河下游冲积平原，粮食作物以小麦、玉米为主。本文数据多来自于课题组在山东省德州市平原县黄河涯镇德州现代农业科技示范园的试验数据，进而更好地解析冬小麦-夏玉米轮作体系典型农田氮素盈余状况，更新氮素盈余指标相关数据，为以后研究提供参考。本试验平台所采取的种植方式、氮肥管理方式等与德州市农户所采取的措施一致，且德州市土壤肥力主要为中产土壤，整个德州市全氮含量标准误差和碱解氮标准误差均较低，分别为0.08 g·kg^-1和27.4 mg·kg^-1[13]，因此具有一定代表性，对德州市的施氮现状和土壤氮素平衡研究有重要的参考意义。试验地土壤类型为砂壤土，pH值为8.54，碱解氮为26.56 mg·kg^-1，年平均降水量547.50 mm，主要集中在6—8月。

1.2 试验设计

本试验共设5个施氮水平，分别以CK（常规施肥：570 kg N·hm^-2·a^-1）、N0（未施肥：0 kg N·hm^-2·a^-1）、N1（减量10%：514 kg N·hm^-2·a^-1）、N2（减量20%：456 kg N·hm^-2·a^-1）和N3（减量30%：400 kg N·hm^-2·a^-1）来表示，其中小麦季施氮量分别为315、0、284、252、221 kg N·hm-（2分基肥和追肥两次施入，比例为1:1），玉米季施氮量分别为255、0、230、204、179 kg N·hm^-2。各处理磷钾肥用量为：玉米季（P₂O₅ 45 kg·hm^-2，K₂O 60 kg·hm^-2）和小麦季（P₂O₅ 270 kg·hm^-2，K₂O 0 kg· hm^-2）。为防止试验小区之间侧渗和串流，各个试验小区用预制板隔离，地下埋深60 cm，田埂包高20 cm。试验所选用的玉米品种为鲁宁184，小麦品种为鲁宁22。试验地的其他农作管理措施如除草、病虫害防治等与当地常规管理措施一致。

1.3 研究方法

本研究以农田土壤为研究界面，研究华北平原冬小麦-夏玉米轮作体系上的氮素平衡。氮素盈余量公式为

(1)

式中，A表示氮素每年在土壤中的氮素盈余量；I_Total表示输入总量；O_Total表示输出总量。三者的单位均为kg N·hm^-2·a^-1。

I_Total为各输入项之和，公式为

(2)

式中：I_f为氮肥输入氮量；I_st为秸秆还田输入氮量；I_se为种子输入氮量；I_i为灌溉水输入氮量；I_a为大气沉降输入氮量。上述5项参数的单位均为kg N·hm^-2·a^-1。

O_Total为各损失项之和，公式为

(3)

式中：O_g为小麦玉米籽粒吸氮量；O_st为小麦玉米秸秆吸氮量；O_r为损失氮量。上述3项参数的单位均为kg N·hm^-2·a^-1。

氮输入项和输出项数据具体来源和算法：

（1）氮输入项

氮肥施用量的数据为课题组在试验田实际施用量；

秸秆还田带入氮=秸秆产量×秸秆还田率×秸秆氮含量

式中秸秆产量为田间实测产量；

灌溉带入氮=整个轮作周期灌溉量×灌溉水氮含量；

种子带入氮=播种量×作物籽粒氮含量；

干湿沉降来自文献汇总结果。

（2）氮输出项

籽粒移出氮=作物籽粒产量×籽粒携出氮量；

秸秆移出氮=秸秆产量×秸秆氮含量；

损失氮包括地表径流、淋洗、氨挥发和N₂O排放，其数据为田间实测数据。

在本试验条件下，输入项中，干湿沉降数据是根据2010—2017年已发表的关于华北地区冬小麦-夏玉米轮作体系下沉降的文章汇总分析的结果；非共生固氮数据参考赵荣芳等^[14]；种子含氮量为田间实际播种量乘以种子含氮量；灌溉水中氮含量为灌溉水量乘以其含氮量（采集水样分析测定其氮含量）。输出项中，作物吸收氮量和氨挥发、N₂O排放以及淋洗、径流等损失氮量均为田间实测数据。作物收获时每个监测点收获2 m×3 m=6 m²的样方进行籽粒和秸秆测产，采集样品利用H₂SO₄-H₂O₂消煮、蒸馏定氮方法测定籽粒和秸秆含氮量；氨挥发是采用通气法进行捕获，具体采集操作和计算方法见文献[15-16]；N₂O排放采用静态箱法进行采集，采集的气样利用气相色谱（Agilent7890，USA）测定，方法详见陈哲等^[17]；收集淋洗和径流的水样并测定水量和水样含氮量进而计算淋洗和径流氮素损失量。

1.4 氮素管理水平的评价方法

为了更加直观表现德州市氮素管理存在的问题及今后改进方向，我们采用基于欧氏距离的评价方法对其氮素管理指标进行评价。首先是建立目标系统（The goal system），在本研究中，我们参考巨晓棠等^[9-10]的各指标参考值作为目标系统，这一系统就是未来氮素盈余指标的方向。目标系统就是该系统评价的标尺，一个系统现阶段的状态与目标系统之间的差距就可以衡量该系统的发展状态，用D表示，如公式（4）所示：

(4)

式中：D代表系统与目标系统的距离，f_i代表任一系统，F代表目标系统。此外，引入最差状态系统，是目标系统的完全对立面，所有指标均为最差，基于此，将这些取值称之为零系统（The zero system）。通常，零系统的取值为0，但对于氮素损失量（氨挥发、N₂O排放和淋洗）的目标系统是0，这样对应的零系统取值就不能为0，要根据情况设定合理的零系统。这样一来，被评价系统的所有指标均处于零系统与目标系统之间，进而方便进行欧氏距离运算。

引入欧氏距离（Euclidean distance）的概念，用dIG来代表被评价系统与目标系统的加权欧氏距离，如公式（5）所示：

(5)

式中：X_ik代表被评价系统的第k个评价指标，X_Gk代表目标系统的第k个评价指标值。S_k代表第k个指标值的标准差。由于上述所说的X_ik和X_Gk数据数值量纲不同，因此需要将其标准化，如公式（6）所示：

(6)

式中：x代表原始值，m代表均值，S代表标准差。如果施肥量等指标超出区间范围的情况，将按公式（7）处理：

(7)

式中：Y代表变换后的值；A代表目标系统的指标参考值；X代表被评价系统指标值。如果转换后结果为负值，那么Y直接取值0（认为评价指标值过高，导致对被评价系统影响效果为负）。由此，可以建立氮素管理指标评价计算公式，公式如（8）所示。

(8)

式中：A代表氮素管理水平；d_IG代表被评价系统与目标系统的距离；d_ZG代表零系统与目标系统的距离，即最大距离。A值处在[0, 1]，这样可以符合通常的认知习惯，A值越大，氮素管理水平越接近理论目标状态，即氮素各管理指标越好。0 < A < 0.60表示低级氮素管理水平，0.60≤A < 0.80表示中级氮素管理水平，0.80≤ A≤1表示高级氮素管理水平。

1.5 数据分析

利用Microsoft Office Excel 2003进行数据整理和图形绘制。

2 结果与分析 2.1 德州市全年和各季氮素盈余量

全年、冬小麦季和夏玉米季氮素盈余表见表 1和表 2。在整个冬小麦-夏玉米轮作周期中，由于秸秆还田，秸秆氮素输入和输出项相等，因此表 1中省略了输入项中的秸秆输入氮量和输出项中的秸秆吸收氮量。从结果中可以看出，冬小麦-夏玉米轮作农田氮素输入主要来源于化肥投入；而氨挥发和氮素淋洗是重要的氮素损失去向。通过降低化肥施用量，均不同程度地降低了氮素损失量。从表 2可以发现，各处理的化肥氮投入量基本满足作物生长籽粒和秸秆所需。除化肥氮的投入，还有干湿沉降、非共生固氮、秸秆还田等其他氮素投入方式，这就造成了土壤氮素盈余，加大了环境污染风险。过量的氮肥施用是对资源的一种浪费，未来要合理引导农户适当施肥，以此实现在保证不减产的情况下优化施肥的目标。而在输出项中，氨挥发和淋洗是氮素损失的主要途径，尤其CK处理，两项输出项占其总氮素输出的31.60%。因此，如何有效降低氨挥发和减少氮素淋洗是该地氮素管理的重要难题。

2.2 德州市氮素管理水平的评价

德州市是华北地区冬小麦-夏玉米轮作体系下典型农田区域，对德州市氮素管理水平进行评价，有利于发现德州市氮素管理存在的问题以及今后改进的方向，使之实现优化施氮的目标。同时，本评价结果也可以为山东省或者其他区域的冬小麦-夏玉米轮作体系氮素管理提供借鉴和参考。参考杨世琦等^[12]文章建立目标系统的指标参考值，具体见表 3。由于各处理输入项中干湿沉降、种子、灌溉、非共生固氮和秸秆等带入的氮素一致，且为了清晰易懂，本评价体系将氮素输入项中的干湿沉降、非共生固氮、种子和灌溉等投入的氮统一为一个整体，即其他氮素总输入。通过1.4的评价方法得出德州市氮素管理水平，结果见图 1。N0处理虽然A值较高，但其原因主要是氮素损失量较低。虽然N0处理可以减少氮素损失、增强土壤保护能力，但由于其作物吸氮量也非常低，在实际中不能推广应用。因此，本研究仅将N0处理作为参考，不与其他处理共同分析。由图 1可知，德州市农户常规施肥（CK）处理的A值最低，氮肥减量10%~30%（N1、N2和N3）处理的A值略有增加，但与CK处理一样，均属于中级氮素管理水平。这意味着德州市常规氮素管理水平较低，存在问题较多，不能仅依靠减少化肥施用量来实现优化施肥的目的。为了进一步分析德州市当前氮素管理存在的问题，采用雷达图的方法直观地分析CK处理与零系统和目标系统各指标之间的距离，进而发现其存在的问题。结果（图 2）显示，氨挥发量、N₂O排放量和氮素淋洗量极高对CK处理的A值呈负贡献，也是导致其A值较低的主要原因，即氮素损失是导致A值较低的主要原因。而氮肥施用量、其他氮素总输入和籽粒吸收等处于中间水平，可以有计划地改善这些指标。德州市未来氮素施用主要关注点为减少氮素损失，目前可以通过施用新型缓释肥料、配施有机肥和推广测土平衡施肥技术等方法来解决。

2.3 德州市土壤氮素盈余指标体系的构建

基于上述研究结果，分别探讨全年、冬小麦季和夏玉米季的施氮量、氮素盈余量与氮素损失量（包括氨挥发、N₂O排放、淋洗和径流）的关系，以便建立合理的氮素盈余指标体系。图 3表示的是整个轮作周期施氮量与氮素盈余量及氮素损失量之间的关系，通过关系式可知：当盈余量为50 kg N·hm^-2·a^-1时，施氮量为352.87 kg N · hm^-2 · a^-1，且此时氮素损失量为112.74 kg N·hm^-2·a^-1，对环境危害较小，因此当氮素盈余量为30~70 kg N·hm^-2·a^-1时为德州市氮素盈余推荐值。根据三者之间的关系，建立德州市冬小麦-夏玉米轮作体系下氮素盈余指标体系（见表 4）。

计算出整个轮作周期氮素盈余体系，进一步可以计算各单季对应的结果。图 4表示冬小麦季施氮量与氮素盈余量和氮素损失量之间的关系。当氮素盈余量为30 kg N·hm^-2时，对应的施氮量为195.83 kg N·hm^-2，此时损失量为76.25 kg N·hm^-2，对环境危害较小。因此将30 kg N·hm^-2设置为冬小麦季氮素盈余量推荐值，氮素盈余量在10~50 kg N·hm^-2范围内即为合理。同样的，依据不同盈余量的施氮量和氮素损失量所造成的污染程度设置5个等级，分别为极低、低、合理、高和极高，提出相应的管理建议（见表 4）。

与冬小麦季相同，图 5表示夏玉米季施氮量与氮素盈余量和氮素损失量之间的关系。当氮素盈余量为40 kg N ·hm^-2时，对应的施氮量为182.22 kg N · hm^-2，此时氮素损失量为57.08 kg N·hm^-2，对环境危害小。因此将40 kg N·hm^-2设置为夏玉米季的氮素盈余量推荐值，氮素盈余量在20~60 kg N·hm^-2范围内即为合理。同样的，依据不同氮素盈余量所对应的施氮量和氮素损失量所造成的污染程度设置5个等级，分别为极低、低、合理、高和极高，提出相应的管理建议（表 4）。

3 讨论

土壤氮素平衡包括氮素进入土壤、从土壤中的损失及土壤氮库的变化，可以反映各种氮素损失途径及土壤氮素盈亏。过去的研究中，只是考虑了生育期有机氮的矿化量，而未考虑到无机氮的固持量，这样得出的结果是不完全的^[18]。因此，本研究补充了无机氮固持量的内容，进而得出更加准确的结果。当氮素盈余量为50 kg N·hm^-2·a^-1时，所对应的氮素损失量也较小，施氮量为352.87 kg N·hm^-2·a^-1，此时与巨晓棠等^[9]推算的华北地区全年合理施氮量（350 kg N·hm^-2· a^-1）一致，因此当氮素盈余量为30~70 kg N·hm^-2·a^-1时，为德州市氮素盈余量推荐值，这一推荐值略低于郭婧妤^[19]推荐的60~100 kg N·hm^-2·a^-1，主要原因是德州市长期高氮量施用导致土壤氮素含量本身较高且德州市氨挥发强度显著高于华北地区平均值^[20]，氮素损失量较大。因此，德州市当前的氮素盈余量推荐值应略低于华北地区平均值。在小麦季时，当盈余量为30 kg N·hm^-2时，对应的施氮量为195.83 kg N·hm^-2。这一结果与巨晓棠^[21]推算的华北冬小麦季氮肥推荐量（196 kg N·hm^-2）一致，符合现状条件下农田耕作措施的合理施肥量。在夏玉米季时，当盈余量为40 kg N·hm^-2时，对应的施氮量为182.22 kg N·hm^-2。这一结果与吴良泉等^[22]推荐华北地区夏玉米季氮肥施用量178 kg N·hm^-2一致。为了实现上述目标，德州市冬小麦-夏玉米轮作农田未来施肥方案应按季所需，适当降低氮肥施用量，在保证补充土壤氮库的同时，减少氮素资源浪费，以实现土壤氮素达到相对理想的平衡状态。

利用欧氏距离法对德州市氮素管理水平进行评价，发现常规处理（CK处理）未达到高级氮素管理水平，这主要是由于其极高的氮素损失尤其是氨挥发导致的，而且氮肥减量10%~30%处理的A值依然处于中级氮素管理水平，这说明仅仅依靠减氮的方法不能够根本解决氮素施用存在的问题。CK处理全年氨挥发达到96.71 kg N·hm^-2·a^-1，这与许多研究一致^[23-25]，氨挥发主要集中于作物生长的前期，这主要是由于较大的氮源输入，使得氨挥发速率增加^[26]。另外，7月初至9月初土壤温度和湿度较高，也会促进土壤的氨挥发^[27]。另一方面，CK处理的氮肥投入已经达到了作物所需氮素含量，同时补充的其他氮素输入，极大地增加了氮素损失的风险。未来，德州市农户应从减少氮肥施用量和应用缓释肥等新型肥料减少氮素损失以及提高氮素利用率等多方面入手，即从源头减少氮素损失，提高氮素利用率，缓解氮素流失风险。此外，农作管理措施中要注意降低灌溉量，特别是避免一次性大量灌溉，以减少氮素淋失。运用本评价方法，可以直观地看出当前氮素管理水平与目标管理存在的差距，便于农户理解以及地方政府推广一些利农政策；且本评价方法与其他评价方法相比，最大的优势在于其避免了指标权重分配和专家打分主观因素的影响^[28-29]，降低了人为主观因素的影响。

综合来说，针对德州市氮素管理面临的问题，给出以下建议：第一，根据氮素盈余指标体系，减少氮肥投入量，建议在冬小麦-夏玉米轮作体系下全年施氮量在320~380 kg N·hm^-2·a^-1，其中小麦季和玉米季均为160~210 kg N·hm^-2，小麦季和玉米季可以根据情况适当调整，但全年总量不宜超过380 kg N·hm^-2·a^-1的上限。第二，可以使用缓释肥料等新型肥料，以此减少氮素短时间内的大量损失。第三，合理的施肥方法与施肥时期，注重深施、基肥与追肥的配合等以达到减少氮素损失的目的。第四，依靠科技手段减少氮素损失，提高氮肥利用率，如选择合适的肥料品种、有机肥与无机肥配施等方式，减少氮肥在施用前期的大量损失。第五，注意改善农作管理措施，如避免一次性大量灌溉等等。

4 结论

（1）德州市冬小麦-夏玉米轮作体系整个周期氮素盈余量推荐值为30~70 kg N·hm^-2·a^-1；德州市冬小麦季氮素盈余量推荐值为10~50 kg N·hm^-2；德州市夏玉米季氮素盈余量推荐值为20~60 kg N·hm^-2。建议德州市冬小麦-夏玉米轮作体系下全年施氮量在320~380 kg N·hm^-2·a^-1。

（2）在本研究中，氮肥是氮素输入的主要来源，而氨挥发以及淋洗是氮素输出的主要去向。冬小麦季的氮素污染程度高，夏玉米季氮素污染程度相对较低，综合整个轮作周期氮素损失量较大，应特别注意冬小麦季的氮素管理。

（3）运用欧氏距离法分析出德州市氮素管理水平与理论目标值存在的差距，德州市氮素管理水平（A值）为0.63，属于中级偏低氮素管理水平，造成各处理A值较低的原因主要是高氨挥发量、高施氮量和高淋洗量等。德州市未来应采取控制氮肥施用量、科学施氮、施用新型肥料和改善粗放农作管理措施等方法。