华北平原典型区域农牧系统氮素流动及其环境效应——以河北省为例

佟丙辛，张华芳，高肖贤，侯勇，马文奇



华北平原典型区域农牧系统氮素流动及其环境效应——以河北省为例

佟丙辛1，张华芳1，高肖贤1，侯勇2，马文奇1

（1河北农业大学资源与环境科学学院，河北保定 071001；2中国农业大学资源与环境学院/植物-土壤相互作用教育部重点实验室，北京 100193）

【目的】以华北平原典型地区河北省为例，明确农牧系统氮素养分流动特征和环境效应，分析农牧系统氮素循环利用率和农牧业结合的程度，探讨农牧系统氮素的优化管理途径，为农牧业养分循环和绿色可持续发展提供科学依据。【方法】以河北省“农田-畜牧”生产系统为边界，在整理统计资料、文献数据和调研数据的基础上，利用物质流分析方法，分别定量1980—2015年河北省农田体系、畜牧体系和农牧系统的各个氮素输入和输出项，同时利用氮素利用率、氮素盈余量和氮素回田率等指标分析氮素流动特征与环境效应。农牧系统的氮素输入项主要包括化肥、生物固氮、干湿沉降、灌溉、人粪尿还田、外源饲料；氮素输出项主要包括农田体系主产品的本地消费、外销，畜牧体系主产品的本地消费、外销，农牧系统大气排放、水体排放；内部的氮素循环项主要包括农田系统副产品还田、农田系统主、副产品作为本地饲料、畜牧系统副产品还田。【结果】1980—2015年，河北省农田生产体系氮素年均输入量增加1.9倍，而作物收获氮量仅增长1.5倍，导致农田氮盈余量和损失量分别增加1.7和1.9倍，氮素利用率由47.2%降至41.4%。与有机肥氮投入相比，化肥氮投入占农田总氮投入60%以上，近年来接近70%。区域间农田养分平衡差异大，氮素输入方面，邯郸市和唐山市较高（＞600 kg·hm-2）,承德市最低（＜200 kg·hm-2）;氮素盈余方面，唐山市最高，为267.8 kg·hm-2，衡水市最低，为51.6 kg·hm-2。畜牧体系氮素输入量也明显增加，在2005年达到最大值，为1980年的7.7倍；畜禽产品和粪尿氮素产生量同时呈现增加趋势，尤其是粪尿氮素产生量由21.8×104t增加到115.3×104t；各区域间动物粪尿氮素产生量存在明显差异，其中氮素产生量最高为邯郸市（377.3 kg·hm-2），最低为衡水市（122.6 kg·hm-2）。外源饲料氮素依赖率由60.5%增至72.7%，畜牧粪尿氮素还田率由70.4%降至30.2%，但畜牧体系氮素利用率由6.4%增至16.3%。从农牧系统整体来看，1980—2015年氮素输入总累计量高达9 038.9×104t，化肥氮素投入量约占总氮素投入量的55.7%，外源饲料氮素投入量占总氮素投入量的33.1%，农牧产品累计总输出氮为2 537.4×104t，占总累计输入氮量的28.1%，向大气、水体累计排放的总氮量高达4061.2×104t，约占总累计输入氮量的44.9%。【结论】1980—2015年河北省农牧系统氮素投入量大幅度增加，氮素富集和环境排放严重，氮素利用率偏低，不同区域单位面积氮素平衡存在较大差异，农田生产与畜禽生产之间养分循环严重脱节。因此，应该充分利用本地饲料资源，提高有机肥的还田率，走农牧结合的道路，从而降低因“农牧分离”造成的“高投入-低效率”代价，促进农牧业可持续发展。

氮素；农牧系统；养分循环；环境效应；养分利用率；河北省

0 引言

【研究意义】随着人类生活水平的不断提高，氮素作为重要的养分资源，在农牧业中投入数量不断增加，在满足人们对粮食、畜产品需求的同时，也引发了全球性和区域性的氮环境污染问题[1-3]。因此，在明确区域内农牧业氮素流动特征基础上，实现氮素的优化管理，是农牧业可持续发展亟需解决的问题。华北平原是中国重要的农牧业生产基地，其农牧业氮素利用状况直接影响国内农牧业可持续发展。河北省作为华北平原典型代表，农林牧渔业生产规模位居全国第3位，是中国农牧产品主要产区之一[4]。研究该区域农牧系统的氮素流动特征及其环境效应，对于优化区域养分资源配置、协调农牧业发展与生态环境之间的关系，促进农牧业可持续发展等具有重要指导意义。【前人研究进展】目前，国外对农牧系统氮素养分的平衡及管理问题已有大量的研究[5-6]，这些研究通过对“农田-畜牧”系统的养分流动状况进行综合评价，明确养分流动特征，对于解决氮素养分损失和环境污染问题具有重要意义。国内在国家尺度、典型区域、典型的农田或畜牧系统等方面也开展了一些研究，王激清等[7]借助物质流法建立中国农田生态系统氮素模型，估算中国不同地区的氮养分输入输出以及养分盈余并分析养分产生的环境效应；Ma等[8-9]建立和利用NUFER模型，在国家和区域尺度上从食物链的角度对氮素平衡、环境排放、流动规律及利用效率进行了分析；侯勇等[10]用物质流动的方法，对京郊典型集约化“农田-畜牧”生产系统氮素流动特征进行了研究；张华芳等[11]用类似的方法对河北省农牧系统磷素流动特征及其环境效应进行了分析。这些研究引进物质流分析方法，构建养分流动模型，对于揭示农田、区域和整个食物链养分流动特征和环境效应打下了方法基础，也为农牧系统养分优化管理提供了依据。【本研究切入点】目前，华北地区关于农牧系统氮养分流动特征的研究比较少，对于农牧结合体系下的氮素输入输出、流动特征及其环境效应尚不明确。【拟解决的关键问题】借助NUFER模型，应用物质流分析方法，分析河北省1980—2015年间农牧系统氮养分流动特征，明确不同区域间养分流动在时间和空间上的差异，为华北乃至全国农牧系统的氮素养分优化管理提供借鉴。

1 材料与方法

1.1 研究区域概况

河北省地处北纬36°05′—42°37′、东经113°11′—119°45′，是中国农牧业生产大省。20世纪80年代至今，农作物种植面积已由833.5×104hm2上升到874×104hm2，粮食作物的产量从1.5×107t升至3.4×107t，肉蛋奶总产量由5.5×105t升至1.3×107t，在中国农产品主产区中占有重要地位，以2015年为例，当年粮食产量居全国第8位，农牧业生产总产值居全国第5位[4,12]。

1.2 研究体系界定

以河北省“农田-畜牧”生产体系为研究边界，氮素为研究对象，核算1980—2015年间“农田-畜牧”生产体系每年氮素输入、输出和内部循环项（图1）。体系输入项主要包括：化肥、外源饲料（系统外输入饲料）、农田灌溉、生物固氮和干湿沉降等；体系输出项主要包括：动植物产品的产出（划分为本地居民消费和外销）和氮素环境（大气和水体）排放；而作物主副产品的本地饲用、作物副产品（秸秆残茬）还田和畜禽粪尿还田作为体系内部循环项。

1.3 养分流动项的计算和数据获取方法

1.3.1 农田子系统 （1）输入项计算方法和数据来源 农田输入氮素量包括：化肥氮素量、有机肥氮素量、灌溉氮素量、干湿沉降氮量和生物固氮量。各项算法和数据来源如下。

NIfer（化肥氮素投入量）：来源于中国统计年鉴（2008—2015年）和中国60年农业统计资料[4,12]；复合肥中氮按照文献数据[8,13]给出的比例折算。

图1 农牧生产体系氮素流动模式示意图

NIManure（有机肥输入氮量）：NIManure= NIAM+NISt+NIHM（1）

式中，NIAM代表动物粪尿氮还田量，NISt代表作物秸秆氮还田量，NIHM代表人粪尿氮还田量。

NIAM=QTYAM×NODischarge×NMTR（2）

式中，QTYAM代表动物饲养数量，NODischarge代表动物粪尿排氮量，NMTR代表动物粪尿还田比例[11,14]。

动物粪尿排氮量（NODischarge）参见1.3.2公式（17），还田比例来自笔者课题组文献积累数据和调研数据[15-16]。

NISt=NOSY×NSTR（3）

式中，NOSY代表作物副产品收获氮量，NSTR代表副产品还田比例[11,14]。作物副产品收获氮量，参见1.3.1公式（7），还田比例来自笔者课题组积累数据和文献数据[8,13]。

人粪尿氮还田量（NIHM）由本地家庭消费氮量和还田比例相乘获得，本地家庭消费氮量见1.3.3部分，还田比例采用NUFER模型数据[11,14]。

NIIW（灌溉水输入氮量）：NIIW=CCA×NIIWA（4）

式中，CCA代表作物种植面积，NIIWA代表年单位面积灌溉氮量[11,14]。年单位面积灌溉氮量来自笔者课题组积累数据和文献数据[8,13]。

NIBNF（生物固氮输入氮量）及NIdep（干湿沉降氮量）采用NUFER模型的算法[11,14]。

NIBNF=NFCCA×NINFCA（5）

式中，NIBNF代表生物固氮输入氮量，NFCCA代表固氮作物（花生、大豆等）种植面积，NINFCA代表单位面积生物固定氮量[11,14]。单位面积生物固氮量来自笔者课题组积累数据和文献数据[8,13]。NIdep与NIBNF算法相似，由耕地面积与年单位面积干湿沉降氮量相乘得出。

（2）输出项计算方法和数据来源 农田输出项包括作物收获氮素量、水体排放氮素量和气体排放氮素量。作物收获氮量包括作物主产品收获氮素量和作物副产品收获氮素量。

NOGY（作物主产品收获氮量）=NOYG×NYNC（6）

式中，NOYG代表主产品产量，NYNC代表主产品氮含量[11,14]。

NOYS（作物副产品收获氮量）=NOYG×NRCS×NSC（7）

式中，NOYG代表主产品产量，NRCS代表谷草比，NSC代表副产品含氮量[11,14]。

主要作物主产品的产量参考中国统计年鉴（2008—2015年）、中国60年农业统计资料[4,12]获得，谷草比、作物籽粒含氮量、作物秸秆含氮量、蔬菜氮素含量以及水果氮素含量的相关系数来源于NUFER模型[8,13,17]。

采用NUFER模型方法计算大气和水体氮素排放，其中地表径流、淋洗、侵蚀、氨挥发、氧化亚氮排放产生的氮素损失量，为各个损失项的损失系数与农田表观氮盈余量的乘积[17]。地表径流、淋洗、侵蚀、氨挥发、氧化亚氮排放损失系数参考NUFER模型[8,13,17]。农田表观氮盈余量为农田氮总输入量与作物产品收获氮量之差。

1.3.2 畜牧子系统 （1）输入项计算方法和数据来源 畜牧系统输入项主要包括本地饲料氮投入量和外源饲料氮投入量。

NILC（本地饲料氮投入量）=（NOYG+NOSY）×NRF（8）

式中，NOYG代表作物主产品收获氮量，NOSY代表作物副产品收获氮量，NRF代表饲用比例[11,14]。

NIEC（外源饲料氮投入量）=NIAM-NILC（9）

NIAM=NOAMMP+NOAMSP+NOManure+NOAMB（10）

式中，NIAM代表动物体系氮需求量（等于动物体系氮素输出量），NOAMMP代表肉蛋奶产品氮量，NOAMSP代表动物副产品氮量，NOManure代表动物粪尿氮量，NOAMB代表动物活体含储存氮量，这些项目的计算方法见1.3.2（2）输出项算法，外源饲料主要包括牧草和饲料添加剂等的外地进口。

（2）输出项计算方法和数据来源 动物体系氮素输出包括动物主产品氮量、副产品氮量、粪尿氮量和动物活体氮素的变化量。

NOAMMP（动物主产品氮量）=NOAMMeat+ NOAMMilk+NOAMEgg+NOAMAP（11）

NOAMMeat（肉氮量）=QTYAM×G.WAMA×NORAMMeat×NOAMMeatNC（12）

NOAMMilk（奶氮量）=NOAMMilkP×NOAMMilkNC（13）

NOAMEgg（蛋氮量）=NOAMEggP×NOAMEggNC（14）

NOAMAP（水产品氮量）=NOAMAPP×NOAMAPNC（15）

式中，QTYAM代表动物数量，G.WAMA代表单个动物活体质量，NORAMMeat代表肉所占比例，NOAMMeatNC代表肉含氮量，NOAMMilkP代表奶产量，NOAMMilkNC代表奶含氮量，NOAMEggP代表蛋产量，NOAMEggNC代表蛋含氮量，NOAMAPP代表水产品产量，NOAMAPNC代表水产品含氮量[11,14]。

NOAMSP（动物副产品氮量）=QTYAM×G.WAMA×（NOAMRB×NOAMNCB+NOAMRO×NOAMNCO） （16）

式中，G.WAMA代表单个动物活体质量，NOAMRB代表骨所占比例，NOAMNCB代表骨含氮量，NOAMRO代表其他副产品所占比例，NOAMNCO代表其他副产品含氮量[11,14]。

NOManure（粪尿氮量）=QTYAM×QTYAM×NODischargeA（17）

式中，NODischargeA代表单个动物排泄氮量[11,14]。

NOAMB（活体氮量），其为年际之间动物体系存栏数量的不同而引起的活体氮量变动。

NOAMB=QTYAM×NOAMBNC（18）

式中，NOAMBNC代表动物活体（可食用部分）氮含量[11,14]。

畜禽饲养数量和畜禽产品（肉蛋奶等）产量来源于中国60年农业统计资料及中国统计年鉴[4,12]；畜禽粪尿含氮量、饲料养分的饲用比例来自笔者课题组积累数据和文献数据[15,18-20]；动物活体可食部分氮含量、畜禽出栏酮体重以及肉、骨、副产物比例和含氮量均来源于笔者课题组积累数据[4,13,15-16]。

动物粪尿向水体及向大气排放氮量采用NUFER模型方法计算，由粪尿产生氮数量及相应的向水体、大气排放系数的乘积得到[17]。

1.3.3 农牧系统 整个系统中输入项如化肥、外源饲料、干湿沉降、生物固氮、人粪尿还田量、灌溉水的氮量和输出项如水体、大气排放氮量的计算方法和数据来源见农田和畜牧子系统部分。而动植物产品输出项的计算如下：

NCLC=NCLCCC+NCLCCU+NCLCAC+NCLCAU（19）

NCLCCC=PNLC×PCLCA×CPCN（20）

式中，NCLC代表本地农村及城镇动植物食品氮总消费量，NCLCCC与NCLCCU分别代表农村动植物食品氮本地消费量，NCLCAC以及NCLCAU分别代表城镇动植物食品氮本地消费量[11,14]；PNLC代表农村人口数，PCLCA人均食品植物性氮产品消费量，CPCN代表植物食品氮含量[11,14]；本地农村（城镇）动植物性氮产品与本地农村植物性氮产品消费量算法相同。

NCEC= NOLCM-NCLC-NFLCM（21）

式中，NCEC代表动植物产品外销氮数量，NOLCM代表本地主产品产生氮量，NCLC代表本地食品氮消费量，NFLCM代表动植物产品用作饲料氮量[11,14]。

1.4 评价指标

农田氮盈余量=农田输入氮素量-农田输出氮素量；农田氮养分利用率（%）=作物收获氮素量（主产品+副产品）÷农田输入氮素量×100；畜牧氮养分利用率（%）=（畜牧主产品氮+畜牧副产品氮）÷畜牧体系投入氮×100；本地饲料氮素利用率（%）=本地饲料投入氮素÷本地饲料资源氮素×100[11,14]；外部饲料氮素依赖率（%）=外部饲料氮素÷饲料总投入氮素×100；动物粪尿氮素回田率（%）=粪尿还田氮素÷粪尿产生氮素×100；农牧系统氮素利用率（%）=（植物性食品氮+动物性食品氮）÷农牧系统总投入氮量×100。

2 结果

2.1 农田子系统不同年代氮素输入和输出核算

1980—2015年间河北省农田体系的氮素总输入量和总输出量均呈现明显的增加趋势。总输入量变化可以分为2个阶段，1980—1998年为快速增长阶段，由95.1×104t增至277.8×104t，增长了1.9倍；之后为基本稳定阶段。化肥一直是农田体系氮素输入的主要来源，其投入量一直占总投入量的60%以上，近年来接近70%；而有机肥氮素投入1980—2005年整体呈现增加趋势，2006年以后呈现减少趋势。从氮素总输出量来看，总体趋势与氮素总输入一致，呈现先增加后平稳的趋势。作物收获量作为最重要的氮素输出项，36年间由46.8×104t增至112.2×104t，增长约1.5倍；氨挥发是主要的氮素损失项，约占氮素总输出项的30%左右，36年其年排放总量由21.6×104t增至57.8×104t，增长了1.7倍；此外，径流、淋洗和反硝化造成的氮素损失也在不断增加；氮素水体和大气损失量36年增加了1.9倍。氮素可计算的总输出量占总输入量的比例由1980年的75.4%降至2015年的69.9%，导致农田氮素大量盈余，年均总盈余量由73.2×104t增至190.9×104t，增加了1.6倍（图2）。

2.2 农田子系统不同区域氮素输入和输出核算

2015年农田体系不同区域单位面积氮素输入量不同，邯郸市和唐山市较高（＞600 kg·hm-2），张家口市和承德市较低（＜200 kg·hm-2）。其中化肥氮素投入量邯郸市最高，为463.1 kg·hm-2，张家口市最低，为65.3 kg·hm-2。作物收获氮素量邯郸市最多，为365.1 kg·hm-2，张家口市最少，仅为48.5 kg·hm-2；氮素盈余量唐山市最高，为267.8 kg·hm-2，衡水市最低，为51.6 kg·hm-2（图3）。上述结果表明，河北省不同区域间农田氮素的分配和利用存在着严重的不均衡。

图2 1980—2015年河北省农田系统氮素平衡

图3 2015年河北省不同区域单位面积农田系统氮素平衡

2.3 畜牧子系统不同年代氮素输入和输出核算

河北省畜牧系统氮素输入可分为2个阶段，1980—2005年间，河北省畜牧体系氮素总输入量明显增加，在2005年达到最高，约为201.8×104t，与1980年的氮素输入量相比增加了7.7倍，随后几年的输入量保持在140×104t左右，趋于稳定。相比于本地饲料氮素输入量，系统外源饲料氮素输入量大幅度增加，且明显高于本地饲料氮素的输入量。其中，外源饲料氮素输入量在畜牧体系氮素总输入量中所占的比例呈递增趋势，由60.5%（1980年）增至71.7%（2015年）。36年间，畜牧体系动物产品氮素输出量和动物粪尿氮素资源量，同样表现出明显的增加趋势，尤其是粪尿氮素资源量由21.8×104t增至115.3×104t（图4）。

图4 1980—2015年河北省畜牧体系氮素平衡

2.4 畜牧子系统不同区域氮素输入和输出核算

2015年河北省不同区域畜牧体系单位面积的氮素养分流动数量表现为：本地饲料的氮素供应量最高的是邯郸市，约为105.0 kg·hm-2，最低的是张家口市，约为18.4 kg·hm-2；对于外源饲料的氮素需求量最多的是邯郸市，约为406.2 kg·hm-2，最少的是沧州市，约为105.7 kg·hm-2；从畜牧系统的氮素输出项来说，动物产品的氮素输出量邯郸市最高，约为133.9 kg·hm-2，最低的是沧州市，约为32.8 kg·hm-2；动物粪尿产生的氮素量邯郸市最高，约为377.3 kg·hm-2，衡水市最低，为122.6 kg·hm-2（图5）。畜牧体系氮素的来源和利用也存在严重的区域不均衡。

2.5 农牧系统不同年代氮素利用和去向的变化

河北省农牧系统1980年和2015年总输入氮量分别为88.9×104t和325.7×104t，产品总输出氮量分别为35.1×104t和109.7×104t，产品总输出占总输入的比例分别为39.5%和33.7%。与1980年相比，2015年农牧系统氮素输入量和输出量分别增加2.7和2.1倍，产品输出量占输入量比例降低了近6%。从输入项构成看，1980年到2015年化肥氮素投入比例由66.9%降至57.0%，而外源饲料氮素投入比例则由15.9%增至31.0%。在氮素总输出中，1980年到2015年农田产品氮素输出量比例由95.4%降至82.1%，动物产品占比由4.6%升至17.9%。同时，农田产品和动物产品外销氮素也有大幅度提高。氮素向大气和水体的排放是主要的氮素损失项，1980年分别为27.2×104t和6.2×104t，2015年分别为77.5×104t和80.4×104t，分别增加了1.8和12.0倍。向水体排放氮素占氮素总损失的比例由1980年18.5%升至2015年50.9%；同时，畜牧业的环境损失贡献在加大，其贡献率由1980年的19.5%上升到2015年的51.0%（图6-a、6-b）。农牧结合体系36年养分积累量结果表明，河北省农牧系统累计输入总氮素量高达9 038.9×104t，化肥氮素投入量约占总氮素投入量的55.7%，外源饲料氮素投入量占总氮素投入量的33.1%，农牧产品累计输出总氮素量为2 537.4×104t，占累计输入总氮素量的28.1%，向大气、水体排放的总累计氮素量高达4 061.2×104t，约占累计输入总氮量的44.9%（图6-c）。

2.6 河北省农牧系统氮素利用和循环状况

1980—2015年间农田体系所产生的本地饲料氮利用率在缓慢增加，由45.6%增加到60.5%。同时，随着外源饲料投入比例的不断增加，畜牧生产体系对外源饲料的依赖程度在不断增加，由1980年的60.5%增至2005年的85.6%，之后有所下降并趋于平稳，维持在75.0%左右。1980—2015年间畜牧体系动物粪尿氮素还田比率逐渐降低，由70.4%降至30.2%（图7）。36年间农田系统的氮素利用率由47.2%降至41.4%。畜牧体系的氮素利用率从1980年的6.4%增至2015年的16.3%。农牧结合体系的氮素利用率总体呈现逐渐降低的趋势，由1980年的36.0%降至2015年的31.3%（图8）。

图5 2015年不同区域单位面积畜牧体系氮素平衡

3 讨论

3.1 河北省农田系统氮养分流动特征

农田养分流动特征反映了系统内的养分输入输出量变化、养分利用率、养分管理效果以及对环境的影响，一直都是研究者所关注的重要问题[7,18]。在区域尺度上，国内外很多研究对农田系统氮素平衡进行了探索，例如Krishna等[19]在对印度土表养分平衡（2001—2002年）进行核算后，对国家不同区域的养分盈余状况进行了系统的阐述，并明确了其优化管理方向；国内一些学者则研究了全国及其各个省域农田养分平衡状况[7,17,20]。与这些研究相比，河北省农田氮素流动具有典型的特征：（1）河北省农田体系氮素输入量增速快，36年增长了2倍；区域间极不平衡，单位面积农田氮素输入量最高的邯郸市与最低的张家口市相差近500 kg·hm-2，多数区域单位面积氮肥投入量高于作物推荐量150—180 kg·hm-2水平[21]；（2）2015年河北省农田化肥氮素投入量约占总投入氮量的70%，而国内一些研究结果显示，化肥氮素投入占农田氮素投入的比例分别为47.4%（全国尺度）[22]、56.1%（全国尺度）[7]、61.5%（华北地区）[23]，与这些研究结果相比，河北省化肥氮素投入所占比例偏高；（3）河北省农田氮素利用率逐年降低，至2015年为41.4%，接近李书田等[23]在华北地区研究结果（41.3%），低于Jarvis等[24]关于欧洲西北地区的典型种植体系的研究结果（53.0%）；（4）荷兰MINAS规定[18]的农田氮素盈余量为60 kg·hm-2（沙质土壤）、100 kg·hm-2（黏质土壤），在河北省的11个行政市中，有6个市域的农田氮素盈余量＞100 kg·hm-2，平均农田氮素盈余量为133.1 kg·hm-2，与王激清等[7]研究结果一致；与Sun等[20]建议的100 kg·hm-2（农田氮素年盈余量潜在环境污染指标）和180 kg·hm-2（高风险环境污染的指标）相比，河北省农田系统因氮素大量盈余所带来的环境风险不容忽视。由此可见，河北省农田体系氮素流动具有投入量大、化肥氮素投入占比高、氮素利用率低、区域之间差异大、环境风险高等特征。

图7 1980—2015年河北省本地饲料氮利用率、外源饲料氮依赖率和动物粪尿氮回田率的变化

图8 1980—2015年河北省农田系统、畜牧系统和农牧系统氮素利用率的变化

3.2 河北省畜牧系统氮养分流动特征

近些年来，中国畜牧业养殖方式由农户散养型向小区养殖型、集约化养殖型进行转变，由此带来一系列的大气、水体等环境问题，分析其原因是因为动物粪尿还田率低、管理措施不配套，所导致的氮素损失及其他环境污染[14,17,25-27]。本研究表明，1980—2015年间由于河北省农牧系统生产规模的不断扩大，2015年畜牧体系动物粪尿的氮素产生量大量增加，是1980年的5.3倍。与1980年相比，畜牧系统2015年的外源饲料氮素投入量增加了86.9×104t，而动物产品的氮素输出量却只增加了22.3×104t，外源饲料氮素依赖率从1980年的60.5%增至2005年的85.6%，之后有所下降并趋于平稳，维持在75.0%左右。畜牧体系的氮素利用率，从1980年的6.4%增加到2015年的16.3%，与Jonker等[28-29]的研究结果（25%—35%）相比偏低。河北省畜牧系统不同区域单位面积粪尿氮素产出量平均值为208.7×104t，其中邯郸市、秦皇岛市、唐山市、石家庄市、保定市的单位面积氮素粪尿氮素产出量均高于欧盟单位面积耕地氮素限量标准170 kg·hm-2[30-31]，说明河北省存在因粪尿氮素过量而导致的环境污染风险。综上所述，河北省畜牧体系氮素流动具有外源氮素依赖率高、养分利用率低、粪尿氮素还田率低、区域间不均衡、环境风险高等特征。

3.3 河北省农牧系统氮素流动特征

如何走好种养结合的道路，在充分保障粮食安全和畜牧业高效发展的同时，降低生产过程中的资源环境代价，实现环境友好、清洁生产，已成为目前中国农牧业可持续发展的首要任务[30-34]。欧洲、美国等发达国家通过出台一系列的环境法律法规限制农牧系统的氮素输入量，提高体系内氮素的利用效率，从而减缓氮素所带来的污染问题[32,35]。本研究表明，河北省农牧系统36年累积氮素投入量为9 038.9×104t，化肥氮素投入量约占总氮素投入量的55.8%，外源饲料氮素投入量占总氮素投入量的31.5%，且外源饲料氮的依赖程度在不断增大。1980—2015年间畜禽粪尿氮还田率不断下降，由70.4%大幅度下降到30.2%，明显低于欧盟国家的68%[36]。导致畜禽粪尿氮素还田率低的主要原因包括：（1）畜禽粪尿在“饲舍-储藏-处理”环节管理不合理，导致大量氮素向大气排放，特别是氨挥发损失严重，占总排泄氮量的26.8%；（2）在中大型规模养殖场/小区，缺乏粪污处理设施，未经处理的粪污直接向水体排放（占总排泄氮量的19.0%）；（3）农田生产大量依靠化肥投入，有机肥施用比例较低。由于大量的化肥氮素投入、高外源饲料氮素依赖率和较低的动物粪尿氮素还田率，河北省农牧结合体系中的氮素利用率不断下降，降至2015年的31.2%，明显低于中国农牧系统平均氮素利用率（41%）[37]。由此可见，河北省农牧生产体系在氮素利用和循环方面严重脱节。

3.4 河北省农牧系统生产过程中的环境影响

农牧结合体系生产过程中，对环境的影响主要体现在氮素向大气和水体的排放。本研究结果表明，河北省农牧系统向大气和水体排放的氮素总量分别为33.4×104t（1980年）、157.9×104t（2015年），气体氮素损失比例由1980年的81.4%降至2015年的49.1%；水体排放损失占氮素总损失的比例呈现上升趋势，由18.6%增加到50.9%，与前人研究结果较为一致。刘晓利[37]研究表明，2001年中国农牧系统氮素损失以气态损失为主，占总损失量的66.9%；马林[17]对2005年中国农牧系统氮素损失量进行了估算，其中气态损失占比为52.6%，进入水体的氮素损失为47.4%。在总损失的贡献方面，Zhang等[38]2010年研究结果表明大气氮素损失有47.5%来自畜牧业；Zhou等[39]对京津冀地区氨挥发损失贡献比率进行了研究，其中畜牧业贡献率为56.7%，本研究取得类似研究结果，畜牧体系对农牧系统氮素损失的贡献率由1980年的19.5%上升到2015年的51.0%，说明畜牧业在氮素损失方面的贡献率不断提高。36年间，河北省农牧系统向大气、水体排放的总积累氮素量高达4 061.2×104t，约占总积累输入氮量的44.9%，远远超出欧盟排放量（30%，其中气体26%，水体4%）[36]。由此可见，河北省农牧系统生产过程中对大气和水体的环境影响问题不容忽视。

3.5 优化途径

针对河北省农牧结合体系存在的高氮素投入、低植物产品和动物产品生产氮素利用率、高环境氮素排放，农牧结合体系氮素利用严重脱节的问题，建议从以下途径进行优化：（1）农田系统应减少化学肥料氮素的投入[21]，采取合理的施肥技术[40]（测土配方施肥技术），根据农作物需求确定施肥量；在控制施肥量的同时采用合理的施肥方法，例如沟施、深施，结合前氮后移技术，以及合适的作物种植管理技术，提高氮素利用率[41]；（2）畜牧系统应该从饲料营养和饲养管理方面来入手[28-29,32]，提高动物产品氮素的利用效率，从氮素投入源头和动物养分利用过程来降低后期动物粪尿中的氮素含量；对动物粪尿进行处理，例如建立沼气池生产沼气、进行堆肥、对粪尿进行覆盖等处理，减少粪尿管理不合理造成的氮素损失[42]；（3）农牧结合体系中，应该减少化学肥料的投入[43]，将动物粪尿进行处理后回田，有机肥和化肥配合施用，提高本地农作物产品的数量和品质，提高本地饲料的饲用比例，实现真正的农牧结合；（4）国家应加快制定相关法律法规，借鉴欧盟环境保护法规（硝酸盐法规、欧盟水框架法规等）对氮素肥料施用和动物粪尿储藏与处理进行明确规定，实现清洁生产[30,44]。

4 结论

（1）河北省农田体系氮素流动具有投入量大、化肥氮素投入占比高、氮素利用率低、区域之间差异大、环境风险高等特征。从1980年到2015年，河北省农田体系氮素总输入量、总输出量均明显增加，其中1980年到1996年为快速增长阶段并达到最高峰，之后趋于稳定；农田产品氮素输出量的增长幅度远小于氮素输入量的增长幅度，导致氮素盈余增加、氮素利用率持续降低。

（2）河北省畜牧体系氮素流动具有外源氮素依赖率高、养分利用率低、粪尿氮素还田率低、区域间不均衡、环境风险高等特征。河北省畜牧体系本省饲料（氮素）投入量难以满足畜牧业发展的需求，饲料主要来源于外省饲料的购入。此外，虽然畜禽生产环节的氮素利用率在逐年增加，但是较欧美发达国家，其利用率仍较低。大量粪尿氮素资源的产生及低回田率，导致向大气、水体的大量氮素损失，加剧环境风险。

（3）河北省农牧系统在氮素利用和循环方面严重脱节。较高的化肥氮素投入、高的外源饲料氮素依赖率、低的动物粪尿氮还田率和畜牧业在农牧系统氮素总损失中贡献率的增加都表明农田系统和畜牧系统未能很好地结合，畜牧业氮素损失的贡献越来越大，形成了现在的高环境代价、低产出、低利用率的农牧系统。河北省应该充分利用本地饲料资源，提高有机肥的还田率，走农牧结合的道路，从而降低因“农牧分离”造成的“高投入-低效率”代价，促进农牧业可持续发展。

[1] Nosengo N. Fertilized to death., 2003, 425(6961): 894-895.

[2] Cai G X, Chen D L, Ding H, Pacholski A, Fan X H, Zhu Z L. Nitrogen losses from fertilizers applied to maize, wheat and rice in the North China Plain., 2002, 63: 187-195.

[3] 朱兆良, 孙波, 杨林章, 张林秀.我国农业面源污染的控制政策和措施. 科技导报, 2005, 23(4): 47-51.

Zhu Z L, Sun B, Yang L z, Zhang L X. Policy and countermeasures to control non-point pollution of agriculture in China., 2005, 23(4): 47-51. (in Chinese)

[4] 中华人民共和国农业部. 新中国60年农业统计资料. 北京: 中国农业出版社, 2009.

Ministry of Agriculture of the People’s Republic of China.. Beijing: China Agriculture Press, 2009. (in Chinese)

[5] Oenema O, Witzke H P, Klimont Z, Lesschen J P, Velthof G L. Integrated assessment of promising measures to decrease nitrogen losses fromagriculture in EU-27., 2009, 133(3/4): 280-288.

[6] Velthof G L, Oudendag D, Witzke H P, Asman W A H, Klimont Z, Oenema O. Integrated assessment of nitrogen losses from agriculture in EU-27 using MITERRA-EUROPE., 2009, 38(2): 402-417.

[7] 王激清, 马文奇, 江荣风, 张福锁. 中国农田生态系统氮素平衡模型的建立及其应用. 农业工程学报, 2007, 23(8): 210-215.

Wang J Q, Ma W Q, Jiang R F, Zhang F S. Development and application of nitrogen balance model of agro-ecosystem in China., 2007, 23(8): 210-215. (in Chinese)

[8] Ma L, Ma W Q, Velthof G L, Wang F H, Qin W, Zhang F S, Oenema O. Modeling nutrient flows in the food chain of China., 2010, 39(4): 1279-1289.

[9] 马林, 魏静, 王方浩, 马文奇, 张福锁. 中国食物链氮素资源流动特征分析. 自然资源学报, 2009, 24(12): 2104-2114.

Ma L, Wei J, Wang F H, Ma W Q, Zhang F S. Analysis on the feature of nitrogen flow from food chain perspective in China., 2009, 24(12): 2104-2114. (in Chinese)

[10] 侯勇, 高志岭, 马文奇, Heimann L, Roelcke M, Nieder R. 京郊典型集约化“农田-畜牧”生产系统氮素流动特征. 生态学报, 2012, 32(4): 1028-1036.

Hou Y, Gao Z L, Ma W Q, Heimann L, Roelcke M, Nieder R. Nitrogen flows in intensive “crop-livestock” production systems typically for the peri-urban area of Beijing., 2012, 32(4): 1028-1036. (in Chinese)

[11] 张华芳, 高肖贤, 侯勇, 任月同, 马文奇. 河北省农牧系统磷素流动及其环境效应. 河北农业大学学报, 2013, 36(2): 17-22.

Zhang H F, Gao X X, Hou Y, Ren Y T, Ma W Q. Phosphorus flow in agro-livestock system and its environmental effect of Hebei Province., 2013, 36(2): 17-22. (in Chinese)

[12] 中华人民共和国国家统计局. 中国统计年鉴. 北京: 中国统计出版社, 2008-2015.

National Bureau of statistics of the People’s Republic of China.. Beijing: China Statistical Press, 2008-2015. (in Chinese)

[13] Ma L, Velthof G L, Wang F H, Qin W, Zhang W F, Liu Z, Zhang Y, Wei J, Lesschen J P, Ma W Q, Oenema O, Zhang F S. Nitrogen and phosphorus use efficiencies and losses in the food chain in China at regional scales in 1980 and 2005., 2012, 434: 51-61.

[14] 张华芳. 河北省农牧生产体系氮磷养分流动特征及调控途径[D]. 保定: 河北农业大学, 2013.

Zhang H F. Nitrogen and phosphorus flow in agro-livestock system and strategies of optimization in Hebei province[D]. Baoding: Hebei Agricultural university, 2013. (in Chinese)

[15] 许俊香, 刘晓利, 王方浩, 张福锁, 马文奇, 马林. 我国畜禽生产体系中磷素平衡及其环境效应. 生态学报, 2005, 25(11): 2911-2918.

Xu J X, Liu X L, Wang F H, Zhang F S, Ma W Q, Ma L. Phosphorus balance and environmental effect of animal production in China., 2005, 25(11): 2911-2918. (in Chinese)

[16] 赵路, 魏静, 马林, 王方浩, 马文奇. 河北省不同养殖模式的畜禽粪尿资源及污染风险分析. 农业环境科学学报, 2009, 28(3): 544-548.

Zhao L, Wei J, Ma L, Wang F H, Ma W Q. The assessment and comparison on animal manure pollution risk of different breeding patterns in Hebei., 2009, 28(3): 544-548. (in Chinese)

[17] 马林. 中国食物链氮素流动规律及调控策略[D]. 保定: 河北农业大学, 2009.

Ma L. Mechanism and regulatory strategies of nitrogen flow in food chain of China[D]. Baoding: Hebei Agricultural university, 2009. (in Chinese)

[18] Oenema O, Kros H, Vries W. Approaches and uncertainties in nutrient budgets: implications for nutrient management and environmental policies., 2003, 20: 3-16.

[19] Krishna Prasad V, Badarinath K V S, Yonemura S, Tsuruta H. Regional inventory of soil surface nitrogen balances in Indian agriculture (2000-2001)., 2004, 73: 209-218.

[20] Sun B, Shen R P, Bouwman A F. Surface N balances in agricultural crop production systems in China for the period 1980-2015., 2008, 18(3): 304-315.

[21] 朱兆良. 农田中氮肥的损失与对策. 土壤与环境, 2000, 9(1): 1-6.

Zhu Z L. Loss of fertilizer N from plants-soil system and the strategies and techniques for its reduction., 2000, 9(1): 1-6. (in Chinese)

[22] 陈敏鹏, 陈吉宁. 中国种养系统的氮流动及其环境影响. 环境科学, 2007, 28(10): 2342-2349.

Chen M P, Chen J N. Nitrogen flow in farming-feeding system and its environmental impact in China., 2007, 28(10): 2342-2349. (in Chinese)

[23] 李书田, 金继运.中国不同区域农田养分输入、输出与平衡. 中国农业科学, 2011, 44(20): 4207-4229.

Li S T, Jin J Y. Characteristics of nutrient input/output and nutrient balance in different regions of China., 2011, 44(20): 4207-4229. (in Chinese)

[24] Jarvis S, Hutchings N, Brentrup F, Olesen J E, van der Hoek K W. Nitrogen flows in farming systems across Europe//Sutton M A.. Cambridge, UK: Cambridge University Press, 2011: 211-227.

[25] 王方浩, 马文奇, 窦争霞, 马林, 刘晓利, 许俊香, 张福锁. 中国畜禽粪便产生量估算及环境效应. 中国环境科学, 2006, 26(5): 614-617.

Wang F H, Ma W Q, Dou Z X, Ma L, Liu X L, Xu J X, Zhang F S. The estimation of the production amount of animal manure and its environmental effect in China., 2006, 26(5): 614-617. (in Chinese)

[26] 黄鸿翔, 李书田, 李向林, 姚杰, 曹卫东, 王敏, 刘荣乐. 我国有机肥的现状与发展前景分析. 土壤肥料, 2006(1): 3-8.

Huang H X, Li S T, Li X L, Yao J, Cao W D, Wang M, Liu R L. Analysis on the status of organic fertilizer and its development strategies in China., 2006(1): 3-8. (in Chinese)

[27] 全国农业技术推广服务中心. 中国有机肥料养分志. 北京: 中国农业出版社, 1999.

National Agricultural Technology Extension Service Center.. Beijing: China Agriculture Press, 1999. (in Chinese)

[28] Jonker J S, Kohn R A, High J. Dairy herd management practices that impact nitrogen utilization efficiency., 2002, 85: 1218-1226.

[29] Huhtanen P, Nousiainen J I, Rinne M, Kytölä K, Khalili H. Utilization and partition of dietary nitrogen in dairy cows fed grass silage-based diets., 2008, 91: 3589-3599.

[30] Oenema O. Governmental policies and measures regulating nitrogen and phosphorus from animal manure in European agriculture., 2004, 82(E. Suppl.): E196-E206.

[31] Wang F H, Dou Z X, Ma L, Ma W Q, Sims J T, Zhang F S. Nitrogen mass flow in China’s animal production system and environmental implications., 2010, 39: 1537-1544.

[32] Lenis N P, Jongbloed A W. New technologies in low pollution swine diets: Diet manipulation and use of synthetic amino acids, phytase and phase feeding for reduction of nitrogen and phosphorus excretion and ammonia emissions., 1999, 12(2): 305-327.

[33] Osei E, Hauck L, Jones R, ogg c, Keplinger k. Livestock and the environment: Lessons from a national pilot project. 2008, http://tiaer.tarleton.edu.

[34] Steinshamn H, Thuen E, Bleken M A, Brenøe U T, Ekerholt G, Yri C. Utilization of nitrogen (N) and phosphorus (P) in an organic dairy farming system in Norway., 2004, 104(3): 509-522.

[35] Neeteson J J. Nitrogen and phosphorus management on Dutch dairy farms: Legislation and strategies employed to meet the regulations., 2000, 30: 566-572.

[36] Oenema O, Oudendag D, Velthof G L. Nutrient losses from manure management in the European Union., 2007, 112: 261-272.

[37] 刘晓利.我国“农田-畜牧-营养-环境”体系氮素养分循环与平衡[D]. 保定: 河北农业大学, 2002.

Liu X L. Nitrogen cycling and balance in “agriculture-livestock- nutrition-environment” system of China[D]. Baoding: Hebei Agricultural university, 2002. (in Chinese)

[38] Zhang Y, Dore A J, Ma L, Liu X J, Ma W Q, Cape J N, Zhang F S. Agricultural ammonia emissions inventory and spatial distribution in the North China Plain., 2010, 158(7): 490-501.

[39] Zhou Y, Cheng S Y, Lang J L, chen d S, zhao b b, liu c, xu r, li t t. A comprehensive ammonia emission inventory with high-resolution and its evaluation in the Beijing-Tianjin-Hebei (BTH) region, China., 2015, 106: 305-317.

[40] 张福锁. 测土配方施肥技术. 北京: 中国农业大学出版社, 2011: 3-22.

Zhang F S.. Beijing: China Agricultural University Press, 2011: 3-22. (in Chinese)

[41] 张福锁, 马文奇, 陈新平. 养分资源综合管理理论与技术概论. 北京: 中国农业大学出版社, 2006.

Zhang F S, Ma W Q, Chen X P.. Beijing: China Agricultural University Press, 2006. (in Chinese)

[42] Hou Y, Velthof G L, Oenema O. Mitigation of ammonia, nitrous oxide and methane emissions from manure management chains: a meta-analysis and integrated assessment., 2015, 21: 1293-1312.

[43] Zhu Z L, Chen D L. Nitrogen fertilizer use in China–contributions to food production, impacts on the environment and best management strategies., 2002, 63: 117-127.

[44] 王敬国, 林杉, 李保国. 氮循环与中国农业氮管理. 中国农业科学, 2016, 49(3): 503-517.

Wang J G, LIN S, LI B G. Nitrogen cycling and management strategies in Chinese agriculture., 2016, 49(3): 503-517. (in Chinese)

（责任编辑 岳梅）

Nitrogen flow and environmental effects of crop-livestock system in typical area of North China Plain-A Case Study in Hebei province

Tong Bingxin1, Zhang Huafang1, Gao Xiaoxian1, Hou Yong2, Ma Wenqi1

(1College of Resources and Environmental Sciences, Hebei Agricultural University, Baoding 071001;2College of Resources and Environmental Sciences, China Agricultural University/Key Laboratory of Plant-Soil Interactions, Ministry of Education, Beijing 100193)

【Objective】The objective of this study is to evaluate the characteristics of nitrogen (N) flow and environmental effects in crop-livestock system in Hebei Province, North China Plain at regional level for the period of 1980-2015, analyze the nitrogen cycle use and recycling efficiency and the degree of combination of agriculture and animal husbandry, discuss the improved management strategies of nitrogen, and to provide a basis for the sustainable development of crop-livestock system in China. 【Method】Taking the crop-livestock system in Hebei province as the research boundary, based on the material flow analysis (MFA) method, the input-output balance, recycling and loss of N in the crop system, livestock system and crop-livestock system of Hebei Province were analyzed from the year 1980 to 2015, supplemented with the data from national and regional statistics, literature and farmers’ surveys. N use efficiency, surplus and recycling ratio were considered as the indexes to analyze N flow characteristics and environmental effects. The N input items of crop-livestock system mainly included chemical fertilizer, biological nitrogen fixation, deposition, irrigation, human feces and urine reused as fertilizer and imported feed; the N output items mainly included the locally consumed and exported main product and by-product of crop and livestock systems, N lost to water body and air; the internal N recycling items mainly included by-product of crop system returned to field, main product and by-product of crop system as local animal feed, and livestock manure returned to field. 【Result】In the crop system, the total N input and the amount of N in harvested crops increased by 1.9 and 1.5 times from 1980 to 2015. The N surplus and loss from crop system increased by 1.7 and 1.9 times, respectively, leading to the decrease in N use efficiency (from 47.2% to 41.4%). Compared with the organic fertilizer N input, the chemical N fertilizer input accounted for more than 60% of the total N input in farmland, and in recent years it was close to 70%. The N balance of farmland varied greatly among regions. In terms of N input, Handan and Tangshang were higher (＞600 kg·hm-2), and Chengde was the lowest (＜200 kg·hm-2). In terms of N surplus, Tangshang was the highest (267.8 kg·hm-2), while Hengshui was the lowest (51.6 kg·hm-2). The amount of N input in livestock system also increased significantly, reaching a maximum on 2005, which was 7.7 times higher than that in 1980. The amount of N in animal products and animal excretion increased rapidly. Especially the N production of feces and urine increased from 21.8×104t to 115.3×104t. The N produced by manure was the highest in Handan (377.3 kg·hm-2), and Hengshui was the lowest production (122.6 kg·hm-2). The N dependence rate of exogenous feed increased from 60.5% to 72.7%, the N return rate of livestock feces and urine decreased from 70.4% to 30.2%, but the N utilization efficiency of livestock system increased from 6.4% to 16.3%. In the crop-livestock system, the 36-year accumulated N input was as high as 9 038.9×104t. The chemical N fertilizer input accounted for about 55.7% of the total N input, and N input via imported feed occupied 33.1% of total N input. For the total accumulated N output, the N output via crop and animal products (2 537.4×104t), accounting for 28.1% of the total input N. The 36-year accumulated N loss into air and water was 4 061.2×104t, accounting for 44.9% of the total N input. 【Conclusion】 From 1980 to 2015, the amount of N input in crop-livestock system in Hebei Province increased significantly, N enrichment and environmental emissions were serious, N use efficiency was low, and there were great differences in N balance per unit area in different regions. The N cycle between crop production and livestock production is seriously disjointed. Therefore, we should make full use of local feed resources, improve the return rate of organic fertilizer, and take the road of combination of agriculture and animal husbandry, so as to reduce the cost of “high input and low efficiency” caused by “separation of agriculture and animal husbandry” and promote sustainable development of agriculture and animal husbandry.

nitrogen; crop-livestock system; nutrient cycling; environmental effect; nutrient use efficiency; Hebei Province

2017-05-31；

2017-07-29

国家自然科学基金面上项目（31272247）、国家国际科技合作专项项目（2015DFG91990）

佟丙辛，E-mail：tbx19880806@126.com。

马文奇，E-mail：mawq@hebau.edu.cn