配方施肥对水稻生产和经济、环境效益的影响

2022-01-24韩雪梅杨林章俞映倞徐红艳朱秀红杨梖程浩

农业环境科学学报 2021年12期

韩雪梅，杨林章，俞映倞*，徐红艳，朱秀红，杨梖，程浩

（1.苏州市吴中区农业环境与土壤肥料站，江苏苏州 215104；2.江苏省农业科学院农业资源与环境研究所，农业部长江下游平原农业环境重点实验室，南京 210014；3.苏州市吴中区越溪街道经济服务中心，江苏苏州 215104；4.苏州市吴中区农技植保推广站，江苏苏州 215104；5.南京万德斯环保科技股份有限公司，南京 211122）

测土配方施肥是联合国开发计划署推行的一种作物科学施肥管理技术，其在国际领域的推广应用可追溯至上世纪中叶。我国对测土配方施肥的初探始于1992年，2005年我国明确了“提高农业综合生产能力”的目标，构建“测、配、产、供、施”多环节链性联动体系，由政府提供技术和经费上的双向支持[1−2]，推动测土配方施肥技术在广大种植区的应用。受社会发展的影响，如今单纯性增产、保产已不再是种植业唯一的需求，兼顾生产、经济、环保等多方向的综合性需求应运而生。但是，测土配方施肥方案的制定依然沿用初期的老方式、老办法，对综合性需求中除生产以外其他需求的响应存在不足。

具体来说，现阶段综合性需求以保产、稳产为前提，经济上需优化提升投产比，环保上需降低种植区农业面源污染排放，生产可持续性上需维持、提升土壤地力。首先，现有配方施肥对作物生产的需求回应较好，其对水稻生物量、籽粒产量的正向影响已得到广泛证实[3−5]。其次，配方施肥对经济需求的回应主要体现为养分增效。肥料利用率的提升是诸多研究展现施肥方案有效性的重要指标[4，6−7]，然而对投产比（单位质量大米产出的肥料投入）的变化缺少关注。此外，配方施肥对环保和生产可持续需求未作出明确回应。作为肥料优化运筹的技术手段之一，配方施肥理论上可大幅度避免土壤中速效态氮磷的大量盈余，在种植区域内部减少可流失的物质源，以此降低氮磷流失风险[8−11]，然而专门针对测土配方施肥方案实施对氮磷流失影响的研究并未见报道。少数研究关注到测土配方施肥有提升土壤有机质及速效态氮磷钾含量的作用[12]，但是多数测土配方施肥相关研究对土壤前后地力变化的关注度不够，而这恰恰是支撑种植业生产的关键，也是对已有配方进行调整的依据。综上所述，因地制宜地将生产以外的经济、环保等因子纳入考虑，及时对当前配方施肥模式进行调整，成为当前种植业生产对测土配方施肥提出的新要求。

江苏省苏州市吴中区紧邻太湖，境内河道湖泊众多，水网密布。该区域农耕历史悠久，自古就是“鱼米之乡”。种植系统的肥料施用是该区域农业面源污染的主要来源之一。水稻作为该区典型粮食作物，种植面积达2 000 hm2以上，稻田系统的肥料施用、养分利用效率等直接关系到周边流域水体的水质情况[13]。在维持稻米产出的同时，控制稻田氮磷流失，是寻找农业生产和环境保护两者平衡的关键。该区域积极推进“两减六治三提升”专项行动，而测土配方施肥就是落地“治理太湖水环境”工作的主要途径之一。2015—2016 年吴中区农业生产部门开展稻田田间试验和土壤地力情况摸底，明确土壤现有养分含量、水稻产量目标、水稻生长肥料需求等参数，以“大配方、小调整”[14]为原则，应用县域测土配方施肥专家系统，制定肥料配方方案，并经专家组论证微调，最终确定主推氮、磷、钾比例为16%∶12%∶17%的配方肥（由华昌化工股份有限公司生产）作为基肥，搭配尿素为蘖穗肥，作为该区域水稻种植的配方施肥模式。该模式自2017 年稻季起开展全区推广，覆盖吴中区主要农村合作社和种植类工商资本企业，已连续执行4 个种植季。与常规施肥模式相比，当前配方施肥模式的养分利用效率如何？是否降低稻田系统的面源污染负荷？在经济投产上是否具有优势？连续使用后土壤地力现况如何？这些问题都有待验证。

对此，本研究选取连续实施配方施肥模式4 个种植季的稻田区域，设置配方施肥模式和常规施肥模式条件下的缺素处理，探讨配方施肥模式对水稻养分利用、农田氮磷排放、农户经济收益的影响，分析配方施肥模式的经济效益、环境效益和生产可持续性，基于当前地力现况提出配方施肥模式调整建议，为农业绿色发展工作的推进提供支持。

1 材料与方法

1.1 供试土壤及作物

1.1.1 供试土壤

本研究开展于吴中区临湖镇湖桥村（31°8′12″N，120°30′28″ E），位于长江三角洲太湖流域，属亚热带湿润季风气候，年均温15.7 ℃，平均降水1 100 mm。该区域水稻土多由湖底沉积物发育而成，耕作层深度18 cm，经过百年水稻耕作，近5 a以水稻−绿肥为主要轮作模式。试验开展前耕作层土壤基本理化性质：有机质14.9 g·kg−1，全氮1.79 g·kg−1，全磷0.44 g·kg−1，碱解氮839 mg·kg−1，有效磷14.2 mg·kg−1，速效钾54.7 mg·kg−1，阳离子交换量15.7 cmol·kg−1，容重1.23 g·cm−3，pH 7.88。

1.1.2 供试作物

本研究供试水稻品种为南粳46 号。该品种为江苏省农业科学院粮食作物研究所选育而成的晚粳，具有株型紧凑、分蘖力中等偏强等特点，是供试区域乃至太湖流域的主要水稻栽种品种之一。

1.2 试验设计及田间管理

试验设置常规施肥（CF）和配方施肥（SF）2 种模式，每种模式分别包含缺氮（PK）、缺磷（NK）、缺钾（NP）和不缺素（NPK）处理，另设1个不施用任何肥料的空白处理，共计9个处理。每种处理设3个重复，小区面积64 m2（8 m×8 m）。其中，缺素处理氮磷钾养分投入分别施用尿素、过磷酸钙和氯化钾；常规施肥的不缺素处理使用氮、磷、钾比例为15%∶15%∶15%，常规复合肥作为基肥，配方施肥的不缺素处理使用氮、磷、钾比例为16%∶12%∶17%的配方复合肥为基肥；蘖肥和穗肥期施用尿素完成氮补充。2种施肥模式的养分投入量以当地实际种植的氮磷钾投入情况为指导，即常规施肥模式施用常规复合肥525 kg·hm−2和尿素450 kg·hm−2，配方施肥模式施用配方复合肥350 kg·hm−2和尿素375 kg·hm−2。具体肥料养分投入量如表1所示。

表1 试验处理及其氮磷钾肥投入量（kg·hm−2）Table 1 Treatments and their nitrogen，phosphorus and potassium input（kg·hm−2）

稻季试验开展于2020 年，基肥、蘖肥和穗肥施用时间分别为6 月8 日、7 月1 日和7 月29 日；6 月9 日移栽水稻幼苗至试验区内，间距为18 cm×30 cm，11 月10日收获。

1.3 样品采集及指标测定

1.3.1 植株样品采集及指标测定

收获时，统计各小区平均分蘖数，以此为标准，采集每个小区典型植株样品（仅地上部分）4 穴。植株样品穗和茎叶部分分开，于105 ℃烘箱内杀青30 min，在80 ℃温度下分别烘至恒质量，称量、记数，用于生物量核算。此后，植株样品穗和茎叶分别磨碎，过20目筛，测定其中氮（凯式法，Foss Scino KT260，中国）、磷（钼锑抗比色法，SHIMADZU UV−1800，日本）、钾（火焰光度计法，上海光谱，中国）含量。各处理理论产量通过籽粒考种计算而来。

1.3.2 土壤样品采集及指标测定

2020年稻季种植前采集5个试验区域耕作层（0～18 cm）土壤样品，作为初始土壤，用于土壤基本理化性质分析（详见1.1.1），水稻收获后采集常规施肥和配方施肥2 种模式下的不缺素处理土壤，检测土壤全氮（同植株1.3.1）、全磷（同植株1.3.1）和有机质（重铬酸钾氧化−比色法）含量。

1.3.3 径流液的采集及氮磷损失的测定

水稻生长期发生较大规模降雨且引起排水，则视为1 次径流。产流2 h 后，于各小区排水口采集对应处理的径流液，由降雨量、田面水高度变化的差值确定径流深，用于单次径流量的估算。试验期间共采集到5次降水径流液（6月13日、6月29日、8月4日、9月16日和9月17日），并未发生烤田排水。

1.4 数据处理及统计分析

1.4.1 作物氮磷钾养分吸收量的计算

各小区每穴植株茎叶和穗分别称质量，记录干物质量和对应氮磷钾含量，用于作物氮磷钾吸收量的计算：

1.4.2 作物氮磷钾肥利用效率的计算

（1）氮磷钾表观利用效率

（2）氮磷钾农学利用效率

（3）偏生产力

记录单次降水前、降水后田面水深和降水量数据，计算径流深。使用流动分析仪（Skalar，荷兰）检测径流液中氮、磷含量。通过单次径流氮磷损失量的累计获得整个稻季径流氮/磷损失量。

1.4.4 生产经济收支的情况

本研究的经济测算以吴中区临湖镇苏州市湖嘉生态农业有限公司规模化生产实际投入为依据。生产支出包括物料投入和人工水电投入，物料包括肥料和种苗，人工水电包含农机使用、灌溉用水电和田间管理相关人工。由于该公司以大米为产出进行售卖，因此，支出上还包含稻米烘干和出米包装两个后期加工环节，收入上则以成品米售卖获益为计。此处计算不包括土地租赁费用。

1.4.5 数据处理

研究采用SPSS 19.0 中Duncan 法对作物产量、产量结构、径流氮磷流失量等数据进行差异显著性分析（P＜0.05），采用T 检验对作物氮吸收量、氮肥利用效率等数据进行差异显著性分析（P＜0.05）；采用Excel软件对数据进行制图。

2 结果与分析

2.1 水稻生产及养分利用情况

2.1.1 产量及产量结构

两种肥料模式下，NPK 处理的产量差异并不显著，说明配方施肥模式未对水稻产量产生不利影响，而养分缺素投入处理的水稻产量均出现一定程度的下降（5%～44%，图1）。这其中，氮和磷的缺失对产量的影响达到显著水平。此外，PK 处理的产量已接近CK处理，可见当前地力状况下，氮肥投入对水稻种植产出起决定性作用。

水稻单位面积的穗数及穗粒数是受投入缺素影响最为明显的指标，也是直接影响水稻产量的关键。配方施肥模式较常规施肥模式在养分投入缺素时，相对穗粒数更多、结实率也更高（表2），因此，配方施肥模式NK 和NP 处理的理论产量高于常规施肥模式下的该处理。在缺氮条件下，配方施肥模式无显著优势，暗示着氮素是保障产出的关键。

表2 水稻产量结构Table 2 The yield structure of rice

2.1.2 养分吸收情况及利用效率

常规施肥模式下，磷、钾元素的缺失降低了8%～13%的水稻氮吸收，而配方施肥模式条件下水稻氮吸收对钾的缺失无明显响应（图2a）。氮的缺失降低了水稻磷吸收16%～52%（图2b），降低了钾吸收4%～63%（图2c）。此外，两种施肥模式下水稻对养分的利用仅在钾元素上存在显著差异，即配方施肥显著提高了水稻对钾的吸收量（图2c）。

表观利用效率数据显示，配方施肥模式较常规施肥模式分别提高了水稻对磷、钾肥利用效率6.35个和28.59个百分点（表3）。这说明，配方施肥模式下的磷钾肥转化利用效率更高。此外，农学利用效率数据显示，配方施肥模式较常规施肥模式显著提高了氮肥对水稻产量贡献，即每公斤氮肥投入可多生产2.81 kg水稻。从总体养分投入（偏生产力）角度上看，单位质量的肥料投入在配方施肥模式下能够输出更多的水稻产量。对肥料利用率的综合提升，体现了配方施肥模式中氮磷钾养分的投入量和配比更符合水稻生长的需求，能够支撑水稻生长的养分吸收，这也一定程度解释了配方施肥模式在减少养分投入条件下保障水稻产出的原因。

表3 肥料利用效率Table 3 Fertilizer use efficiency

2.2 氮磷径流流失情况

在径流深一致的情况下，配方施肥模式显著降低了稻季36%径流途径的氮流失量，对径流途径的磷流失影响不显著（图3）。从氮肥投入占比上看，配方施肥模式径流途径氮损失约为4%，而常规施肥模式为6%。将CK 处理的氮流失量作为土壤本底进行扣除（减去CK 处理对应数值），得出配方施肥模式的因肥料引发的氮流失量仅为常规施肥模式的一半。

2.3 水稻生产经济收支情况

配方施肥模式较常规施肥模式显著降低了35%的肥料经济投入，且增加了4%左右的产出（表4）。因此，在其他支出都相同的情况下，可为农户增加2 558元·hm−2的净利润。

3 讨论

3.1 配方施肥模式的经济、环境效益分析

测土配方施肥以养分补偿学说为原则，通过调整投入的肥料物质及投入量，协调土壤、肥料、作物三者的物质关系。社会的高速发展潜移默化地影响着生产者和管理者心中生产、经济和环境的权重，测土配方施肥的工作侧重也随之发生转移。推广之初，提高作物单产、保障粮食供给是测土配方施肥的主要目标。随着种植规模化、产业化的推进，促进节本增效、获取最优经济收益逐渐取代了对高产的一味追求。“十三五”期间，农业绿色发展对农业生产提出了环境友好和发展可持续的新要求，削减肥料流失量、培肥土壤也纳入了种植生产的目标之中。因此，对配方施肥模式的评判离不开对其经济和环境效益的综合考量。

本研究中配方施肥模式较常规施肥模式降低每吨大米生产所需的肥料投入38%（图4a），减少每吨大米生产产生的氮磷流失量35%（图4b）。由此可见，在现地力条件下，配方施肥模式较常规施肥模式能够更好兼顾生产、经济和环境的多方效益。

本研究所在的吴中区通过自主招标配方肥生产、实行配方肥购买和售卖补贴，将农户获取配方肥的单价降低至常规复合肥的65%。同时，配方施肥模式在追肥期的肥料施入量较常规施肥模式减少了17%，进一步降低了用于肥料购买的投入成本，极大提高了农户对配方施肥模式的选择意愿。对经济不发达地区而言，政府能够用于补贴配方肥购买销售的资金并不宽裕，便捷的配方肥购置途径和基层技术培训对于让农户选择推荐的施肥模式显得尤为重要[15−16]。而配方施肥模式本身能否通过提升作物产出提高农户净收益，才是推广、落实测土配方工作的关键[3，5]。值得注意的是，本研究中整个水稻生产的资金投入在13 000元·hm−2左右，而肥料投入为1 000～1 500 元·hm−2（表4），仅占物料投入的35%～45%，总投入的7%～11%，和前人研究中肥料化肥在农业生产中的成本占比数据相差较远[17]。如此低的占比，一定程度削弱了农户对肥料投入量的关注，使农户比起精准控量减少投入，更愿意选择过量投施确保作物生长养分需求。事实上，过量或不恰当的肥料施用不仅不能提高作物产量，反而会带来一定程度的减产[18]，增加土壤中未被作物吸收利用的肥料源氮磷存量，提升氮磷流失风险[19]。本研究结果显示，氮肥是影响水稻产量最主要的元素，同样的结果在油菜系统的测土配方施肥中也得到了验证[20]。因此，氮相较于磷对于产量和污染风险的调控更为重要。本研究中两种施肥模式的径流氮损失比例在4%～6%范围内（图3），与种植业氮肥平均径流损失系数5%左右的前期研究结论相符，说明配方模式能够更好呼应水稻对氮素的需求。径流发生时土壤环境中的矿质态氮留存量较低，实现稻田氮排放的有效削减[10−11]。随着各地政府对农业面源污染治理的重视，减少肥料施用量成为实现种植业源头削减的主要抓手，测土配方“按需、按量、按时”施肥的原则[21]，为“减什么”“减多少”“怎么减”给出了答案[8，10]，而测土配方工作的落实也是推动农业面源污染治理的重要工作组成[9]。

表4 两种肥料施用模式的水稻生产收支（元·hm−2）Table 4 Income and expenditure of rice planting in two fertilizer patterns（yuan·hm−2）

3.2 配方施肥模式的生产可持续性分析及改进建议

现有配方施肥模式已沿用4 a，有关施肥配方调整周期的研究此前并未见报道。就本研究而言，对比2017 年使用配方模式前的土壤养分指标数据，土壤全氮、全磷未出现显著变化，但土壤有机质含量下降显著（表5）。这说明，现有配方施肥模式可满足水稻生长需求，但是，即使实行秸秆还田，也未能弥补土壤中有机物质耗损缺口。此外，本研究配方施肥模式下水稻地上部分氮、磷吸收量分别为169 kg·hm−2和45 kg·hm−2（图2），仅占氮磷投入量的66%和72%，可见在养分投入配比不变的情况下，配方施肥模式的养分投入量依然存在削减的余地。因此，在肥料品种上，建议适当增加一定比例的有机肥，以30%～70%占比为宜[22]；在肥料养分投入上，建议削减蘖肥、穗肥的尿素投入量，将一定比例的配方肥由基肥期移至穗肥期，并在基肥期使用尿素进行氮补充，在保证磷钾投入量的同时，减少氮投入总量。考虑到多环节的推行需要一定操作周期[1，15]，太过频繁地调整缺乏可行性，因此，建议以4～5 a为周期，进行配方施肥模式的调整。

表5 实施配方施肥模式前后土壤氮、磷、有机质含量的变化（g·kg−1）Table 5 Changes of soil nitrogen，phosphorus and organic matter after using special formulated fertilizer patterns（g·kg−1）