李硕，李晓欣，赵善丽，3，郝曼，柏兆海，马林，3*

（1.中国科学院遗传与发育生物学研究所农业资源研究中心/河北省土壤生态学重点实验室/中国科学院农业水资源重点实验室，石家庄 050022；2.河北大学生命科学学院/生命科学与绿色发展研究院，河北 保定 071002；3.中国科学院大学，北京 100049）

华北平原作为我国粮食和畜禽产品的主产区，其粮食和肉蛋奶产量分别占全国总产量的34.4%和37.0%[1]。集约化玉米-小麦轮作体系产量不断提高的同时，耕地利用强度高、水肥资源投入量大、耕地质量下降快、盈余养分和温室气体大量排放等严重危害生态环境安全[2-3]。规模化畜禽生产快速发展导致饲料严重依赖外地引进、农牧分离问题凸显、畜禽粪尿循环利用率下降[4-6]。近年来，在玉米-小麦轮作体系中采用的农田养分增产提效技术提高了作物产量和养分利用效率[7-8]，但未充分考虑养殖产生的固液粪肥的合理施用。因此，构建农牧（种养）结合的“粮-饲”轮作体系，筛选固液粪肥科学施用方法已经成为华北平原农牧业绿色发展的关键。

根据畜禽饲料养分需求设计轮作体系能缓解我国畜禽生产对进口饲料的过度依赖，亦能大幅降低农牧生产中的氮素损失及其对环境的危害[9-10]。在华北平原进行集约化“粮-饲”轮作种植能够在有限水、热等资源条件下获取更多的优质饲草，有利于规模化畜禽的可持续生产[10-11]。随着我国“粮改饲”政策的实施，区域种植结构将大幅改变[11]。黑麦草、甜高粱和高丹草具有较高的饲料营养价值，多用于规模化的“粮-饲”轮作高效生产[12-13]。目前，我国在饲草方面的研究，主要侧重于饲草光能利用效率的改善及其干物质和粗蛋白产量的提高[14-15]，但是缺少不同“粮-饲”轮作系统氮素平衡的研究。

加强畜禽生产中的粪尿管理，促进粪尿养分高效利用，能够降低环境污染风险。研究发现合理施用固液粪肥可提高土壤碳、氮储量，并减少养分环境流失[16]；陶晓婷等[17]发现，在小麦越冬期补灌60～120 m3·hm-2的猪场肥水即可满足穗期30～60 kg N·hm-2化肥的氮素营养需求；也有研究表明，在常规施氮水平（225 kg N·hm-2）下，固态猪粪、鸡粪单施均可获得与化肥处理相当的小麦、玉米产量，且小麦产量随牛粪与化肥配施比例的提高而有所提升[18]。Zhang 等[19]经Meta 分析发现旱作生产中粪肥替代部分化肥可增加6.6%的作物产量，同时减少农田NH3和N2O 排放，有利于农业的可持续生产。有机、无机肥料配施亦可提高水稻、大麦及其他经济作物产量[20-22]。目前针对固态/液态粪肥施用的研究主要侧重于单一作物的农业生产[17-22]，对“粮-饲”轮作方式下的施用增产及其环境效应的研究很少。

本文依托中国科学院栾城农业生态试验站的粮-饲轮作试验平台，探讨了不同“粮-饲”轮作方式（籽粒玉米-小麦、青贮玉米/豇豆间作-填闲黑麦草、甜高粱-黑麦草、青贮玉米-黑麦草、高丹草-黑麦草）和不同施肥方式（化肥表施、固态粪肥替代50%化肥表施、液态粪肥替代50%化肥注射施用）对植物地上部生物量及其氮携出量、土壤硝态氮残留量和土壤氮素表观平衡的影响，为丰富该地区“粮-饲”轮作方式、加强畜禽粪肥农田合理利用、促进农牧可持续生产提供理论依据。

1 材料与方法

1.1 试验地概况

试验地位于太行山山前平原的中国科学院栾城农业生态试验站农牧轮作长期定位试验基地（37°53′N，114°40′E，海拔50.1 m）。该试验区属于典型暖温带半干旱季风气候，年均气温12.51 ℃，年均降水量536 mm，其中70%集中在7—9 月，年均蒸散量850 mm，年均日照时长2 340 h，年均无霜期206 d，年均有效积温4 665 ℃。土壤类型为褐土类灰黄土种。耕层（0～20 cm）土壤质地为粉砂壤土，其中砂粒、粉粒和黏粒含量分别为13.8%、66.3%和19.9%。耕层土壤有机碳和速效氮、速效磷、有效钾含量分别为7.2 g·kg-1和92.0、9.5、295.5 mg·kg-1。

1.2 试验设计

试验在2017—2018 年期间开展，采用完全随机区组设计，包括2个因素：（1）轮作方式，即籽粒玉米-小麦（T1）、青贮玉米/豇豆间作-填闲黑麦草（T2）、甜高粱-黑麦草（T3）、青贮玉米-黑麦草（T4）和高丹草-黑麦草（T5）；（2）施肥方式，即化肥表施（C）、固态粪肥表施替代50%化肥氮（S，简称：固态粪肥表施）和液态粪肥注射施用替代50%化肥氮（L，简称：液态粪肥注射）。试验共15个处理，每个处理设4个重复，各处理小区面积为90 m2（6 m×15 m），小区之间保护带宽1 m。籽粒玉米、青贮玉米、豇豆于6 月初播种，10月初收获，品种分别为郑单958、远缘1号和天优一号，其中T1和T4处理玉米行间距均为60 cm，T2处理玉米行间距为75 cm；甜高粱和高丹草于6月初播种，7 月底和10 月初各刈割一次，品种分别为F438 和冀草2号，其中甜高粱行间距为40 cm，高丹草行间距为20 cm；小麦于10 月初播种，次年6 月初收获，品种为石新633，行间距为20 cm；黑麦草于10 月初播种，次年4 月中旬和6 月初各刈割一次，品种为冬牧70，行间距为20 cm。籽粒玉米-小麦轮作方式在玉米和小麦收获后，秸秆全部粉碎还田。试验用粪肥取自河北省石家庄市栾城区典型集约化牛场，经固液分离后分别用于田间，固态和液态粪肥理化性质如表1所示。

夏秋季作物生长期施肥量为300 kg N·hm-2、180 kg P2O5·hm-2和180 kg K2O·hm-2，冬春季作物生长期施肥量为 240 kg N·hm-2、120 kg P2O5·hm-2和 120 kg K2O·hm-2，其中磷、钾肥以基肥形式一次性施用于作物播种前，氮肥分基肥和追肥两次施用，基追比为1∶1。籽粒玉米和青贮玉米追肥期为大喇叭口期，小麦追肥期为拔节期，其他作物追肥期为第一次刈割期，人工表施追肥后立即灌溉。基肥按照等氮原则施用，固态和液态粪肥的磷、钾不足时以化肥补足，化肥和固态粪肥采用人工撒施，随后旋耕入土；液态粪肥则由课题组研制的专用机具（肥水注射机具）注射到5～10 cm 土层；化肥选用尿素（N，46%）、过磷酸钙（P2O5，12%）和硫酸钾（K2O，53%）。玉米苗期、拔节期、大喇叭口期，以及小麦苗期、分蘖期、拔节期各灌水1次，每次灌水量约500 m3·hm2，采用畦灌方式。由于甜高粱和高丹草第一次刈割正处于玉米大喇叭口期，因此甜高粱、高丹草和玉米灌溉时期、方式和灌溉量保持一致；黑麦草前两次灌溉与小麦保持一致，而第3 次安排在刈割期；在4 月中旬填闲黑麦草小区不进行刈割和灌溉操作。根据田间病虫害状况及时施用农药，进行防治。其他田间管理措施保持一致。

表1 新鲜固态和液态粪肥理化特性Table 1 Physical and chemical characteristics of fresh solid and liquid farm manure

1.3 样品采集与氮素测定

1.3.1 土壤样品采集与测定

作物收获后，每个小区随机选取一个样点，采集0～200 cm 土层土样，层间距20 cm，测定土壤含水量和-N、-N 含量。土壤含水量采用烘干法测定；土壤-N 和-N 测定是将新鲜土样混匀后过2 mm筛，然后用1 mol·L-1KCl溶液浸提（水土比5∶1）过滤，流动分析仪（OI，美国）测定。

1.3.2 植株样品采集与测定

在作物收获或刈割期，每个小区随机采集2 m2（2 m×1 m）植株样品，风干，脱粒。小麦、甜高粱、高丹草和黑麦草植株样品分籽粒和秸秆称质量，而玉米样品分籽粒、玉米穗轴和秸秆3 部分称质量，然后以采样面积折算生物量。取部分植株样品烘干、粉碎，用浓H2SO4·H2O2消解，采用凯氏定氮法测定植株样品的氮含量。

1.4 数据计算与统计分析

本研究探讨一个完整轮作周期内农田土壤氮素表观平衡，其中农田土壤氮输入指化肥和畜禽粪肥所带入的氮素，未考虑种子、沉降和灌溉等形式的氮输入；氮输出指作物通过收获或刈割移出农田的氮总量（地上部氮携出量-秸秆还田氮量），未考虑淋溶流失、氨挥发、反硝化等途径氮损失。参考司东霞等[23]的方法计算周年作物氮肥偏生产力和0～200 cm 土层氮素表观平衡参数：

试验数据绘图和统计分析分别采用Excel 2010和DPS 7.05统计软件进行，Duncan多重比较和T检验分析不同处理的影响差异，多因素方差分析因素间的交互效应，差异显著性水平为0.05。

2 结果与分析

2.1 不同处理地上部生物量

在夏秋季（图1a），T5×C 处理地上部生物量最高（45.4 Mg·hm-2），显著高于其他处理。针对3 种施肥方式整体而言，轮作方式T3、T4 和T5 地上部生物量较T1 分别增加44.8%、10.9%和71.0%，T2 地上部生物量较T1 降低13.3%；针对5 种轮作体系整体而言，固态粪肥表施和液态粪肥注射处理地上部生物量较化肥表施分别显著降低7.4%和11.0%。在冬春季（图1b），T1×C 处理地上部生物量最高（16.1 Mg·hm-2），显著高于除T4×C、T1×L 和T1×S 外的其他处理；轮作方式 T2、T3、T4 和 T5 地上部生物量较 T1 分别显著降低34.0%、23.4%、17.4%和22.5%，而施肥方式及其与轮作方式的交互作用对地上部生物量却无明显影响。如图1c 所示，周年地上部生物量T5×C 处理最高（58.2 Mg·hm-2）；轮作方式T3 和T5 较 T1 分别显著增加17.7%和33.9%，而T2 较T1 却显著降低21.5%；固态粪肥表施和液态粪肥注射处理地上部生物量与化肥表施相比显著降低7.3%和9.2%。

2.2 地上部氮携出量与偏生产力

在夏秋季（图2a），地上部氮携出量与其地上部生物量趋势一致，T5×C 处理最高（476.6 kg·hm-2）；与T1相比，轮作方式T5和T3地上部氮携出量分别显著增加141.2%和125.6%，而T2 却降低19.6%；与化肥表施相比，固态粪肥表施和液态粪肥注射处理地上部氮携出量分别显著降低12.6%和15.0%。在冬春季（图2b），T1×C 处理地上部氮携出量最高（225.1 kg·hm-2）；轮作方式T2 地上部氮携出量较T1 显著降低20.2%，而施肥方式及其与轮作方式的交互作用对地上部氮携出量均无明显影响。从图2c可知，周年地上部氮携出量T5×C处理最高（691.8 kg·hm-2）；与T1相比，轮作方式T3和T5显著增加52.8%和61.0%，而T2却显著降低20.0%；与化肥表施相比，固态粪肥表施和液态粪肥注射处理显著降低10.0%和11.4%。

在夏秋季（图3a），T5×C 处理作物氮肥偏生产力最高（151.3 kg·kg-1），显著高于其他处理；与T1 相比，轮作方式T2、T4、T3 和T5 的作物氮肥偏生产力分别显著提高118.1%、264.4%、179.3%、330.4%；与化肥表施相比，固态粪肥表施和液态粪肥注射处理作物氮肥偏生产力分别显著降低7.6%和11.1%。在冬春季（图3b），T4×C 处理作物氮肥偏生产力最高（59.5 kg·kg-1），与T5×C、T4×L、T5×L、T3×S 和T4×S 均无显著差异，而显著高于其他处理；与T1 相比，轮作方式T3、T4 和T5 作物氮肥偏生产力分别显著提高61.5%、74.1%和63.4%，而施肥方式及其与轮作方式的交互作用对作物氮肥偏生产力却无明显影响。从图3c可知，周年作物氮肥偏生产力，T5×C 处理最高（107.9 kg·kg-1），而T1×L 处理最低（30.9 kg·kg-1），且后者显著低于除T1×C 和T1×S 外的其他各处理；轮作方式T2、T4、T3 和T5 较T1 分别显著提高119.8%、132.9%、175.1%和212.9%；固态粪肥表施和液态粪肥注射处理较化肥表施分别显著降低7.4%和9.6%。

2.3 土壤-N残留量

在夏秋季作物收获期（图4a～图4c），随土壤深度增加，各处理土壤-N 含量在 0～40 cm 土层均呈先降低后稳定的趋势，而在40～200 cm 各土层均保持相对稳定的趋势；在冬春季作物收获期（图4d～图4f），随土壤深度增加，各处理土壤-N 含量在0～40 cm土层均呈先降低后稳定的趋势，而在40～200 cm 各土层均保持相对稳定的趋势。在夏秋季作物收获期（图5a），与 T1 相比，轮作方式 T4、T3 和 T5 的 0～200 cm 土层NO-3-N 残留量分别显著降低16.8%、33.9% 和43.5%，而T2 和T1 之间无显著差异；与化肥表施相比，固态粪肥表施和液态粪肥注射处理0～200 cm 土层NO-3-N 残留量分别显著降低20.9%和16.1%。在冬春季作物收获期（图5b），轮作方式T2、T3 和T5 的0～200 cm 土 层 NO-3-N 残 留 量 较 T1 分 别 显 著 降 低22.3%、28.4%和40.1%，而T4 与T1 之间无显著差异；与化肥表施相比，固态粪肥表施和液态粪肥注射处理0～200 cm 土层-N 残留量分别显著降低34.4%和22.2%。

2.4 氮素表观平衡

从表2 可以看出，在夏秋季，与T1 相比，轮作方式T2、T3、T4和T5作物氮输出量分别显著增加21.9%（25.8 kg·hm-2）、241.1%（258.0 kg·hm-2）、62.0%（72.9 kg·hm-2）和265.8%（312.8 kg·hm-2），而T2却显著降低49.0%（137.9 kg·hm-2）；与化肥表施相比，固态粪肥表施和液态粪肥注射处理作物氮输出量显著降低12.9%（36.4 kg·hm-2）和14.4%（440.8 kg·hm-2）。在冬春季，T2 将地上部生物量全部还田，无生物量移除农田，其作物氮输出量视为0；与T1 相比，轮作方式T3、T4 和T5 冬春季作物氮输出量分别显著增加20.7%、27.0%和23.4%，且T3、T4 和T5 之间无显著差异。与T1 相比，轮作方式 T3、T4 和 T5 周年作物氮总输出量显著增加 113.3%、41.6% 和 124.8%，而 T2 却降低49.0%；与化肥表施相比，固态粪肥表施和液态粪肥注射处理显著降低9.9%和12.0%。

由表2 可知，与 T1 相比，轮作方式 T2、T3 和 T5 土壤残留无机氮显著降低20.6%、24.2%和36.8%，而T4与T1之间无显著差异；与化肥表施相比，固态粪肥表施和液态粪肥注射处理的土壤残留无机氮显著降低31.5%和20.6%。在化肥表施中，轮作方式T3和T5土壤氮表观损失小于0，而其他各处理均大于0。与T1相比，轮作方式T3、T4 和T5 土壤氮表观损失分别显著降低 81.5%（218.8 kg·hm-2）、30.8%（82.7 kg·hm-2）和 74.1%（199.1 kg·hm-2），而 T2 却增加 10.0%（26.9 kg·hm-2）；与化肥表施相比，固态粪肥表施和液态粪肥注射处理的氮表观损失显著增加281.0%（173.1 kg·hm-2）和265.9%（163.8 kg·hm-2）。

3 讨论

3.1 不同轮作方式的农田土壤氮素表观平衡

高丹草-黑麦草轮作方式（T5）周年地上部生物量最高，其次为甜高粱-黑麦草（T3）和籽粒玉米-小麦轮作方式（T1），而0～200 cm 土层无机氮残留量和系统氮表观损失却具有相反趋势。本研究中夏秋季轮作方式T2 较大的行间距造成玉米株数降低，进而导致其地上部生物量显著低于T1；冬春季轮作方式T2 实行不施肥的填闲黑麦草操作，进而导致其生物量亦小于T1。尽管冬春季T3、T4和T5轮作方式黑麦草刈割两次，但其地上部生物量仍显著低于T1，表明黑麦草干物质生产潜力低于小麦，这与杨羚等[11]和石鹏飞等[24]研究结果一致。然而，夏秋季T3、T4和T5轮作方式地上部生物量显著高于T1，表明青贮玉米、甜高粱和高丹草的干物质生产潜力明显高于籽粒玉米，这与Lyons等[12]和Samarappuli等[13]研究结果一致。在夏秋季高丹草地上部生物量最高，其次为甜高粱和青贮玉米，这一现象可归因于作物根系结构的不同[13]。研究表明，甜高粱根系可向下延伸至1.6～2.0 m，而玉米根系在土壤分布较浅，主要集中于0～30 cm 土层[25-26]。因此，与玉米相比，甜高粱能够吸收利用深层土壤中的水分和养分，进而生产较多生物量[27]。本研究中无论是夏秋季还是冬春季，轮作方式T1 在0～200 cm 土层的-N 残留量均显著高于其他轮作方式，而轮作方式T5和T3在0～200 cm 土层的-N 残留量却显著低于其他轮作方式，表明轮作方式T5 和T3 更利于作物吸收土壤无机氮，特别是-N。轮作方式T5 和T3 地上部氮携出量较T1 显著增加141.2%和125.6%，表明甜高粱和高丹草能够显著提高粗蛋白产量，利于农牧可持续生产。研究表明，作为青贮玉米的替代作物，甜高粱生长周期短，对养分和农药需求低，蛋白和糖类携出量高，能够在较为干旱和贫瘠的土壤上正常生长[28]，已成为一种具有巨大发展潜力的饲用作物[15，28-29]。本研究结果表明，高丹草的地上部生物量和氮携出量均高于甜高粱，是由于高丹草作为甜高粱和苏丹草的杂交品种具有更强的耐瘠和抗旱能力，而强大的再生能力使得其生物产量高于甜高粱[30-31]。冬春季轮作方式T3、T4 和T5 地上部氮携出量与T1 之间无显著差异，表明冬春季种植黑麦草可提供较多粗饲料。就周年氮肥偏生产力和作物氮总输出量而言，轮作方式T5 和T3 显著高于其他轮作体系，表明甜高粱-黑麦草和高丹草-黑麦草两种轮作方式，特别是后者，可提供的粗饲料最多，更利于华北地区农牧业的可持续发展。就氮表观损失而言，在化肥表施中，轮作方式T3 和T5 土壤氮表观损失小于0，其原因在于：夏秋季T3 和T5 处理甜高粱和高丹草具有较高生物量，氮素携出量高于化肥施用量，作物需要吸收土壤累积氮，以满足自身生长所需；未考虑土壤氮的矿化、大气干湿沉降以及灌溉输入氮；未考虑作物根茬等还田部分的植株残体氮归还，一定程度上低估了氮的投入。该现象表明，在甜高粱-黑麦草和高丹草-黑麦草两种轮作体系中，化肥表施不能满足夏秋季植株快速生长的需要，需适当增加化肥投入量。

表2 不同轮作和施肥方式对作物-土壤氮素表观平衡的影响Table 2 Effects of different rotation system and fertilization management practices on apparent N balance in crop-soil systems

3.2 不同有机肥施用方式的农田土壤氮素表观平衡

与化肥表施相比，固态粪肥表施和液态粪肥注射处理周年地上部生物量分别降低7.3%和9.2%，0～200 cm 土层无机氮残留量显著降低31.5%和20.6%，系统氮表观损失较大。本研究中，与化肥表施相比，固态粪肥表施和液态粪肥注射显著降低夏秋季地上部生物量，而对冬春季地上部生物量无明显影响。这一结果与李燕青等[18]和林治安等[32]在华北平原籽粒玉米-小麦轮作农田施用化肥和粪肥后作物产量趋势大致相同，即试验前期化肥处理的作物产量高于粪肥处理，而随着试验期的延长，粪肥处理的生物量达到甚至超过化肥处理的生物量。Zhang等[19]发现旱作生产中粪肥替代部分化肥可使作物产量提高6.6%，而粪肥完全替代化肥时作物产量却降低9.6%，这主要是由于所施肥料的养分时效性，化肥所提供的矿质氮可供作物直接吸收利用，而固态粪肥所供给的有机氮需要经过矿化缓慢释放，因此具有较强的后效性[17，33]。液态粪肥用于基肥，与作物生长对养分的需求不匹配，经固液分离后的液态粪肥-N 高达3.5 g N·kg-1，在碱性土壤中易以NH3形式排放到大气中，不利于作物后期生长[34]。李硕等[35]在华北平原的研究表明猪场液态粪肥农田注射施用亦可显著提高作物产量，其原因在于液态粪肥注射施用有助于-N被土壤吸附固定，便于后期矿化和作物吸收利用。尽管固液粪肥在短期可能降低作物产量，但其长期与化肥配施可提高作物产量[32]，主要是由于：固液粪肥具有一定的培肥作用，且含有作物生长所需的氮磷钾及中微量元素，能够更好地促进作物生长；土壤微生物能够对其有机养分起到一定的调节作用，使其更能匹配作物生长对养分的需求规律[18]。本研究中，固液粪肥农田施用之所以未能提高作物产量，可能是由于试验期限较短，固液粪肥培肥效果尚未显现；本研究各处理施用肥料量均高，且在作物生长期均实行2 次施肥，氮素供应基本符合作物生长需求规律。因此，固液粪肥对作物氮素吸收利用及其生长状况的效应尚需在长期定位试验中进行验证和深入研究。

本研究中夏秋季固态粪肥表施和液态粪肥注射处理地上部氮携出量较化肥表施处理显著降低12.6%和15.0%，表明固态粪肥表施和液态粪肥注射下粗蛋白产量有所降低，而冬春季固液粪肥各轮作方式下地上部氮携出量均不受施肥方式的影响。综合夏秋季和冬春季，固态粪肥表施和液态粪肥注射处理周年作物氮肥偏生产力较化肥表施仅降低7.4%和9.6%，周年作物氮总输出量亦仅降低12.0%和9.9%，表明固液粪肥可在一定程度上替代部分化肥。土壤氮素表观平衡反映农田土壤氮输入和输出之间的动态平衡，可表征氨挥发、-N淋失等的潜在危害，亦可衡量农田氮养分状况[23-24，39]。本研究氮输入主要来自化肥和固液粪肥，未考虑作物种子、前茬作物残体、灌溉、大气沉降以及试验期间土壤有机氮的矿化。研究中，固态粪肥表施和液态粪肥注射施用下周年作物氮总输出量和土壤残留无机氮较化肥表施均显著降低，表明固态粪肥表施和液态粪肥注射施用能降低土壤残留氮量和淋失风险。在华北平原地下水位持续降低的状况下，向农田施用液态粪肥能够为种子发芽及其前期生长提供急需的水分供给，可在一定程度上节省农田水消耗，降低生产成本，提高种植经济收益，亦可在一定程度上缓解规模化养殖场养殖废水处置难题，利于农牧业可持续发展。因此，在农业生产中，应加大固液粪肥农田施用力度，特别是液态粪肥注射施用。

3.3 华北平原不同轮作系统和施肥方式的养分管理策略

华北平原，特别是河北省农田生产体系在1980—2015 年期间氮素年均输入量远超作物输出氮量，畜牧体系畜禽粪尿氮素年均产生量由2.18×105Mg 增加到1.15×106Mg，而畜牧粪尿氮素还田率由70.4%降至30.2%[40]。可见华北平原农田生产与畜禽生产之间养分循环严重脱节，需要充分利用本地饲料资源，提高畜禽固液粪肥农田施用比例，才能有效解决因“农牧分离”而引发的“高投入、低效率、高污染”的生产问题。在本研究中，轮作方式T5 和T3 土壤残留无机氮和系统氮表观损失较T1 显著降低，且0～200 cm 土层的-N 残留量明显低于起始无机氮累积量（0～200 cm），表明轮作方式T3 和T5 在540 kg N·hm-2的施氮量下仍消耗了大量土壤无机氮，进而在一定程度上降低土壤-N 残留量。Li 等[41]研究表明，在籽粒玉米-小麦轮作方式下，通过采用固态粪肥替代50%化肥表施和液态粪肥替代50%化肥注射施用两种固液粪肥施用方式可显著降低农业生产潜在环境危害，并能提高经济效益和生态效益。综合考虑轮作和施肥方式，本研究周年作物氮总输出量在CT3、LT3、LT5、ST3 和LT5 各处理间无显著差异，且与CT5 差距均小于100 kg·hm-2，表明采用甜高粱-黑麦草和高丹草-黑麦草两种轮作方式生产的同时，选用固态粪肥表施替代50%化肥氮和液态粪肥注射施用替代50%化肥氮两种有机肥施用方式是系统生产力最高且环境效应最佳的绿色种养一体化生产模式。

4 结论

（1）与籽粒玉米-小麦轮作相比，高丹草-黑麦草轮作显著提高了周年地上部生物量，且增幅最高，0～200 cm 土层无机氮残留量和系统氮表观损失降幅最大，其次为甜高粱-黑麦草轮作方式。

（2）与化肥表施相比，固态粪肥表施和液态粪肥注射处理周年地上部生物量分别降低了7.3%和9.2%，0～200 cm 土层无机氮残留量降低了31.5%和20.6%，氮表观损失增加。

（3）在固态粪肥替代50%化肥和液态粪肥注射施用替代50%化肥的施肥方式下，发展甜高粱-黑麦草和高丹草-黑麦草两种轮作方式，有望成为华北平原生产力较高且环境效应较佳的绿色种养一体化生产模式。