-
氮素作物生长和产量形成的必需营养元素[1]。近二十年来我国化学氮肥用量连续升高,当前我国氮肥消耗量占全世界氮肥用量的30%左右[2],高量的氮肥投入,不仅造成了资源的浪费,同时带来了一系列的环境问题,如温室气体排放、氨挥发、地下水硝酸盐污染、水体富营养化等[3-5]。与此同时,我国畜牧业也蓬勃发展,但由于我国农牧分离的现状,畜禽粪便等有机资源未得到充分利用,造成了更多的环境污染[6]。众多学者认为,有效利用有机资源,农牧结合是我国氮肥减施增效的重要措施之一[7-9]。2015年以来,我国先后发布《到2020 年化肥使用量零增长行动方案》,《国务院办公厅关于加快推进畜禽养殖废弃物资源化利用的意见》,《畜禽粪污资源化利用行动方案 (2017—2020 年)》等行动方案,倡导构建种养一体化循环发展机制,提高畜禽粪便利用率,走农业绿色发展道路。因此,将农牧系统作为一个整体,定量分析农牧系统氮素流动途径,是氮素资源高效利用和综合管理的重要内容。
目前,国外对农牧系统氮素循环和环境排放的研究较多,其中以物质流分析[10](Material flow analysis,MFA) 应用最广,它根据质量守恒定律,在一定时间和空间范围内对特定系统的物质及养分流动和贮存进行系统性分析,定量其在各个环节的转化和移动过程,如Isermann等[11]利用物质流分析计算了德国“作物-畜禽-家庭-污水垃圾处理”系统氮素流动帐户,分析了不同营养需求下德国氮素平衡的前景。Ma等[12]利用物质流方法建立了Nufer模型 (Nutrient flows in food chain,environment and resources use),在全国尺度上定量了食物链和农牧系统氮素平衡、环境排放、利用效率以及流动特征。此后,佟丙辛等[9]、张建杰等[13]、张晓萌等[14]利用Nufer模型定量了河北、山西和东北三省农牧系统氮素流动特征,证明Nufer模型在我国区域尺度具有良好的适用性,也提出了提高有机肥还田率,农牧结合等提高区域氮素利用效率的措施。
氮素输入输出平衡 (平衡施氮) 是研究化肥减施潜力的主要方法之一,如米国华等[15]通过此方法计算出东北地区氮肥节约潜力为9~25 kg/hm2,赵亚南等[16]的研究也表明通过平衡施氮原则河南地区小麦节氮潜力为16.4%~36.7%,玉米节氮潜力为30.0%。平衡施氮可降低土壤氮素盈余,降低环境损失。此外,有机替代也是重要的化学氮肥减施措施,宋大利等[17]研究表明我国畜禽粪便全量还田理论上可减少氮肥37.3%;耿维等[18]对安徽畜牧养殖产生的畜禽粪便进行分析,表明其替代化肥潜力为21.8%。农牧系统中可通过平衡施氮和有机替代综合应该更大限度地降低氮肥施用量,Jin等[19]研究表明,通过对综合应用这两种措施,我国氮肥可减施50%以上。当前对京津冀地区氮素养分流动的研究较少,且节肥潜力的研究也均未考虑对“农田–畜牧”整体系统的影响。
京津冀地区是除长三角和珠三角外的第三大经济、政治、交通、文化、科技一体化发展区域,该区域人口密集,养殖业密度大,环境质量要求高,生态压力大[20]。该区域农牧系统整体氮素流动情况、氮肥减施潜力、活性氮减排潜力还不明确,对其研究具有重要的现实意义。本文以京津冀地区农牧系统为研究对象,采用Nufer模型建立氮素账户,定量其流动途径和利用效率,并设置情景分析,探讨平衡施氮 (减少土壤氮盈余和累积) 和有机替代 (增加部分有机肥还田量以替代化学氮肥) 措施下化学氮肥的减施潜力,为该地区实现化肥零增长和种养结合提供科学依据。
-
京津冀地区 (113°27′~119°50′E,36°05′~42°40′N),是我国重要集约化种植区和养殖区。该区域包括北京市、天津市以及河北省三个行政区,土地面积为21.6万 km2,约占全国的2.7%,人口达1.1亿;耕地面积6.54万km2,约占区域总面积的30%,该地区是全国范围内较为集中的肉类消费区和畜禽养殖区,畜禽存栏总量约占全国6.9%,养殖率明显高于全国平均水平[21-22]。2015年4月30日中共中央政治局通过了《京津冀协同发展规划纲要》,提出京津冀一体化协同发展,要求三地在经济、文化、环境等方面协同治理、共同进步。该区域作为未来中国区域发展“三大战略”之一,在中国社会经济发展中具有重要的战略地位。
-
Nufer(Nutrient flows in food chain,environment and resources use) 模型是国家和区域尺度的食物链系统氮磷流动模型[12],该模型基于“土壤-作物-畜禽-食物消费”4个系统,计算各系统间及系统内氮磷流动量、利用率及环境损失量。本研究以京津冀地区“农田-畜牧”生产系统为研究边界,研究系统氮素输入、输出状况及氮素养分在系统内和系统间流动情况 (图1)。
-
系统氮素输入项包括化肥、沉降、生物固氮和饲料等的输入,输出项中有流出系统的供人食用的植物产品、动物产品中所含氮素,环境排放损失的氮素 (包括氨挥发,反硝化脱氮过程的 N2O、N2排放等损失,通过淋溶、地表径流与土壤侵蚀等途径进入到水体,也包括畜禽粪便直接排入水体或闲置的损失),及其他部分 (包括植物生产的籽粒和秸秆在存放和运输过程的损失,和动物生产中的动物骨头、皮毛等部分)。在系统内部,涉及到氮素在农田和畜牧两个子系统的内部循环,主要包括秸秆还田、植物籽粒和秸秆作为饲料进入畜牧系统、畜禽粪便还田等。
-
农田系统氮素输入量 = 化肥 + 畜禽粪便还田 + 秸秆还田 + 沉降 + 生物固氮。其中,畜禽粪便还田量 = 畜禽排泄量 − 畜禽损失量 (包括氨挥发、水体排放与堆置,见畜牧系统);秸秆还田量 = 秸秆氮素携出量 × 还田比例;沉降量 = 耕地面积 × 单位面积年沉降量;生物固氮 = 粮食作物与豆科作物面积 × 单位面积年固氮量。
农田系统氮素输出量 = 籽粒氮携出量 + 秸秆氮携出量 + 氨挥发 + 氧化亚氮 + 淋溶 + 径流和侵蚀 + 反硝化氮气 + 土壤累积。其中,籽粒氮量 = 年产量 × 氮含量;秸秆氮量 = 年籽粒产量 × 谷草比 × 氮含量;环境排放量 = 化肥 × 排放系数 + 畜禽粪便还田量 × 排放系数 (以不同的排放形态分别计算,如氨挥发、氧化亚氮、淋溶等);土壤累积量 = 输入量− 籽粒 − 秸秆 − 环境排放。
-
畜牧系统氮素输入量 = 籽粒氮素投入量 + 秸秆氮素投入量 + 外地输入饲料氮素量。其中,籽粒氮素投入量 = 农田籽粒氮素量 × 饲用比例;秸秆 = 农田秸秆氮量 × 饲用比例;外地输入饲料氮素量 = 本地畜禽氮素输入量 (需求量) − 籽粒氮素投入量 − 秸秆氮素投入量。
畜牧系统氮素输出量 = 动物产品氮素携出量 + 畜禽粪便氮素携出量 + 动物其他部分氮素携出量。其中,动物产品氮素携出量 = 动物产品产量 × 氮含量;畜禽粪便氮素携出量 = 畜禽量 × 排泄量 × 排泄物氮含量;动物其他部分氮素携出量 = 畜禽量 × 畜禽骨头和皮毛氮含量。
-
本研究所用京津冀地区数据 (化肥投入量、农作物播种面积、畜牧业各类动物存栏和出栏量、农作物产量、肉蛋奶产量、人口数量、食物人均消费量等) 均来源于“河北省2018年经济统计年鉴”、“天津市2018年经济统计年鉴”、“北京市2018年经济统计年鉴”和国家统计局数据网站 (http://data.stats.gov.cn/index.htm)。化肥氮素投入中,复合肥按N∶P2O5∶K2O=1∶1∶1的比例计算氮量[23];地区氮沉降量参考尹兴等[24]的研究;谷物和豆科作物固氮量参考李贞宇等[25]的研究;作物籽粒和秸秆含氮量参考马林[26]的研究,作物秸秆谷草比参考牛新胜等[6]的研究,粮食作物籽粒用途比例根据FAO(Food and Agriculture Organization of the United Nations,http://www.fao.org/home/en/) 统计数据确定,秸秆利用比例根据高利伟等[27]的研究确定,化肥不同氮素损失系数根据胡春胜等[5]的研究确定,各畜禽种类的产物系数和氮素含量根据耿维等[18]的研究确定,畜禽粪便还田、生产、存储氮素损失系数根据侯勇等[28]的研究确定,畜禽粪便还田比例根据宋大利等[17]的研究确定。其他所需参数来源于Nufer模型[12,26,29]。
-
运行Nufer模型,查看农田系统、畜牧系统和农牧系统养分流动方向与流动量。将土壤累积、畜禽粪便水体排放与堆置这两种最容易调控的途径作为主要调控对象,减少化学氮肥投入,计算化学氮肥减施潜力。
从降低化肥氮投入20%(S1-N20) 开始,以10%为一个阶段,逐级降低化肥氮投入 (S2-N30、S3-N40),至土壤累积氮出现亏缺 (S4-N50),根据平衡施氮原则[30-31],即氮肥投入量计算中不能长期耗竭土壤养分,则上一个减氮量为当前减排潜力。再将畜禽粪便水体排放和推置量的1/2增加至还田畜禽粪便中 (S5-N40+M),以增加有机替代比例,后续减少化肥氮投入直至至出现亏缺 (S6-N50+M,S7-N60+M),从而得出减氮潜力。
-
农田系统氮素利用率 (NUEc,%) = 作物籽粒氮素量/农田输入氮素量 × 100;
畜牧系统氮素利用率 (NUEl,%) = 动物产品氮素量/畜牧体系输入氮素量 × 100;
农牧系统氮素利用率 (NUEc+l,%) = (植物产品氮素量+动物产品氮素量)/农牧系统总输入氮素量 × 100;
氮素养分循环再利用率 (NRR,%) = (作物秸秆还田氮素量 + 作物秸秆作饲料量氮素量 + 畜禽粪便还田氮素量)/(作物秸秆收获氮素量+畜禽粪便产生氮素量) × 100;
食物氮代价 (NCf,kg/kg) = 农牧系统总输入氮素量/(植物产品氮素量 + 动物产品氮素量)。
-
京津冀地区氮素资源量如表1所示,资源总量为324.3万t,其中化肥、秸秆和畜禽粪便所占比例分别为59.1%、12.5%和28.5%。河北地区资源最多占总资源92.3%、天津占5.1%、北京占2.7%。并且,三种资源分别占当地资源总量百分比有所差异,其中化肥资源由高到低为河北 (59.3%) > 北京 (58.5%) > 天津 (55.7%),秸秆资源为河北 (12.8%) > 天津 (10.6%) > 北京 (3.8%),畜禽粪便为北京 (37.6%) > 天津 (33.7%) > 河北 (27.9%)。京津冀地区氮素需求量分别为植物生产160.5万t,动物生产41.3万t,人口消费42.2万t,三类需求量均以河北最高,但三个地区资源量和需求量占京津冀总资源量和需求量的比例并不匹配,北京和天津需求比例高于资源比例,河北需求比例低于资源比例,河北动植物产品存在外输北京和天津的现象。在农牧生产体系中三种需求量存在交互关系,以最大的植物生产来看,京津冀地区氮素资源量高于需求量,该地区氮素资源可以满足地区需求量,并且有很大节约潜力。
表 1 京津冀地区氮素资源量与需求量
Table 1. Nitrogen resources and demand in Beijing-Tianjin-Hebei region
项目
Item河北
Hebei天津
Tianjin北京
Beijing京津冀
Whole region资源量 Resources (× 104 t) 化肥 Fertilizer 177.3 9.2 5.1 191.6 秸秆 Straw 38.3 1.7 0.3 40.4 畜禽粪便 Manure 83.6 5.5 3.3 92.3 总计 Total 299.2 16.4 8.7 324.3 需求量 Demand (× 104 t) 植物生产 Plant production 151.3 7.5 1.8 160.5 动物生产 Animal production 36.2 3.1 2.0 41.3 人口消费 Human consumption 28.6 6.5 7.2 42.2 比例 Proportion (%) 资源 Resources 92.3 5.1 2.7 100.0 人口消费 Human consumption 67.8 15.4 17.1 100.0 注(Note):植物生产需求量是籽粒氮携出量与秸秆氮携出量之和,动物生产是动物产品氮携出量和动物其他部分氮携出量之和,人口消费需求量是人均消费量 (包括动物产品消费与植物产品消费) 乘以人口数量 Plant production demand refer to the sum of grain N and straw N. Animal production refer to the sum of animal product N and animal other parts N. Population consumption demand is the per capita consumption (including consumption of animal products and plant products) multiplied by the population; 资源比例与人口消费比例是指三地区各自占京津冀地区总资源量和人口消费量的比例 Proportions of resource and human consumption refer to the proportions of each region accounting for the total resources and human consumption of the whole region. -
京津冀地区农田系统氮素输入总量为296.1万t,河北、天津、北京分别为274.4、14.7、7.1万t,占总量的92.6%、5.0%、2.4%(表2)。三个地区总输入量中大部分氮素来自化肥,化肥氮素投入占总输入量的64.7%。其次为畜禽粪便和沉降,分别为12.3%和11.5%。输出项中除籽粒和秸秆外,土壤累积和氨挥发是最大的输出项,分别占输出总量的17.3%和12.8%。各地区中土壤累积、氨挥发、氧化亚氮排放、淋溶、径流与侵蚀、反硝化占比最高的地区均为北京,分别占该地区总输出量的38.6%、15.1%、1.4%、10.1%、2.7%、6.8%,且土壤累积量所占比例高于籽粒。农田系统氮素利用效率 (NUEc) 由高到底为河北 (41.2%) > 天津 (38.9%) > 北京 (20.7%),京津冀地区总体NUEc为40.6%。
表 2 京津冀农田系统氮素流动与平衡
Table 2. Nitrogen flow and balance of cropping system in Beijing-Tianjin-Hebei region
项目
Item河北
Hebei天津
Tianjin北京
Beijing京津冀
Whole region输入 Input (× 104t) 化肥 Fertilizer 177.3 9.2 5.1 191.6 畜禽粪便 Manure 33.0 2.2 1.2 36.4 还田秸秆 Returned straw 18.0 0.8 0.2 18.9 沉降 Deposition 31.9 1.8 0.5 34.1 生物固氮 Biological N fixation 14.2 0.7 0.2 15.1 合计 Total 274.4 14.7 7.1 296.1 输出 Output (× 104t) 籽粒 Grain 112.9 5.7 1.5 120.1 秸秆 Straw 38.3 1.7 0.3 40.4 氨挥发 NH3 34.9 1.9 1.1 37.8 氧化亚氮 N2O 3.4 0.2 0.1 3.7 淋溶 Leaching 25.0 1.3 0.7 27.0 径流与侵蚀 Runoff and erosion 6.3 0.3 0.2 6.8 反硝化 Denitrification 8.0 0.5 0.5 9.0 土壤累积 Accumulation 45.5 3.0 2.7 51.2 合计 Total 274.4 14.7 7.1 296.1 农田系统氮素
利用效率 NUEc (%)41.2 38.9 20.7 40.6 -
京津冀地区畜牧系统氮素投入总量为133.6万t,河北、天津、北京分别为119.7、8.6、5.3万t(表3);氮素输入项中饲料外地输入占比最大 (49.3%),其在三个地区中分别占总量的46.5%、64.0%和88.9%;籽粒饲料和秸秆饲料的氮占总输入量的37.3%和13.4%。京津冀地区畜牧系统输出中,畜禽粪便是最大输出项,为92.3万t,占总量的69.1%,且畜禽粪便中氨挥发、进入水体和堆置、还田三种途径分别占总输出量的19.9%、21.9%和27.3%;动物产品 (肉蛋奶) 是第二大输出项,以动物产品输出氮素占总输入量定义的畜牧系统氮素利用效率 (NUEl) 最高为北京 (33.4%),其次为天津 (29.7%) 和河北 (24.3%),京津冀地区整体为25.0%。
表 3 京津冀畜牧系统氮素养分流动与平衡
Table 3. Nitrogen flow and balance of livestock system in Beijing-Tianjin-Hebei region
项目
Item河北
Hebei天津
Tianjin北京
Beijing京津冀
Whole region输入 Input (× 104 t) 籽粒 Grain 47.1 2.3 0.4 49.9 秸秆 Straw 17.0 0.8 0.1 17.9 饲料外地输入 Feed import 55.6 5.5 4.7 65.8 合计 Total 119.7 8.6 5.3 133.6 输出 Output (× 104 t) 动物产品 Animal products 29.1 2.5 1.8 33.4 其他部分 Other parts 7.1 0.5 0.3 7.9 畜禽粪便总量 Total manure 83.6 5.5 3.3 92.3 氨挥发 NH3 volatilization 24.0 1.7 1.0 26.7 水体排放与堆置
via water and manure stacking26.5 1.7 1.0 29.3 还田 Returning to field 33.0 2.2 1.2 36.4 合计 Total 119.7 8.6 5.3 133.6 畜牧系统氮素
利用效率 NUEl (%)24.3 29.7 33.4 25.0 注(Note):动物产品是指肉、蛋、奶等主产品 Animal products refer to main products such as meat, eggs and milk; 其他部分是指骨头、皮毛等不可食用部分 other parts refer to inedible parts such as bones and fur. -
京津冀地区农牧系统氮素投入总量为306.6万t,河北、天津、北京分别为279.0、17.2、10.4万t,占总量的91.0%、5.6%、3.4% (表4)。三个地区农牧系统氮素输入中以化肥占比最高,为62.5%,其次为饲料外地输入,占21.5%,沉降和生物固氮占总输入11.1%和4.9%。输出项中植物产品、土壤累积、农田氨挥发和动物产品是主要的输出项,分别为70.2、51.2、37.8和33.4万t,占总输出的22.9%、16.7%、12.3%和10.9%。各地区各输出项占当地总输入量的比例中,动物产品、土壤累积、农田反硝化、畜牧氨挥发、畜牧水体排放与堆置以北京地区最高,其余输出项以河北地区最高,天津处于中间位置。农牧系统氮素利用效率 (NUEf+l) 由高到低为天津 (34.6%) > 河北 (34.0%) > 北京 (26.8%),京津冀地区为33.8%。秸秆与畜禽粪便的循环再利用率 (NRR) 以河北地区最高,为55.8%,其次天津,为51.7%,北京为42.9%,京津冀地区整体为55.2%。食物氮代价 (NCf) 由高到低为北京、河北、天津,分别为3.7、2.9、2.9 kg/kg,京津冀地区整体为3.0 kg/kg。
表 4 京津冀农牧系统氮素养分流动与平衡
Table 4. Nitrogen flow and balance of crop-livestock system in Beijing-Tianjin-Hebei region
项目 Item 河北 Hebei 天津 Tianjin 北京 Beijing 京津冀 Whole region 输入 Input (× 104t) 化肥 Fertilizer 177.3 9.2 5.1 191.6 沉降 Deposition 31.9 1.8 0.5 34.1 生物固氮 Biological N fixation 14.2 0.7 0.2 15.1 饲料外地输入 Feed import 55.6 5.5 4.7 65.8 合计 Total 279.0 17.2 10.4 306.6 输出 Output (× 104t) 植物产品 Plant products 65.8 3.4 1.0 70.2 动物产品 Animal products 29.1 2.5 1.8 33.4 土壤累积 Accumulation 45.5 3.0 2.7 51.2 农田氨挥发 NH3-farmland 34.9 1.9 1.1 37.8 农田氧化亚氮 N2O-farmland 3.4 0.2 0.1 3.7 农田淋溶 Leaching-farmland 25.0 1.3 0.7 27.0 农田径流和侵蚀 Runoff and erosion-farmland 6.3 0.3 0.2 6.8 农田反硝化 Denitrification-farmland 8.0 0.5 0.5 9.0 畜禽粪便氨挥发 NH3-livestock manure 24.0 1.7 1.0 26.7 畜禽粪便水体排放与堆置 Into water and stack- livestock manure 26.5 1.7 1.0 29.3 其他部分 Other parts 10.5 0.7 0.3 11.4 合计 Total 279.0 17.2 10.4 306.6 农牧系统氮素利用效率 NUEc+l (%) 34.0 34.6 26.8 33.8 农牧系统氮素循环再利用率 NRR (%) 55.8 51.7 42.9 55.2 食物氮代价 NCf (kg/kg) 2.9 2.9 3.7 3.0 注(Note):其余部分包括植物生产的籽粒和秸秆在存放和运输过程的损失,和动物生产中的动物骨头、皮毛等部分 The other parts include the loss of seeds and straw produced by plants during storage and transportation, and animal bones and fur in animal production process. 京津冀地区农牧系统环境排放量和百分比如表5所示,环境排放总量为140.4万t,占系统总输出量的45.8%,河北、天津、北京环境排放量分别为128.2、7.6、4.6万t,分别占各地区系统总输出量的45.9%、44.3%、44.0%。其中损失量最大为氨挥发 (包括农田与畜牧氨挥发),京津冀地区合计64.5万t,占总输出量的21.0%,其次为畜牧水体排放与堆置,排放量29.3万t,排放率9.5%。各地区中也均以氨挥发与畜牧水体排放与堆置排放百分比最高。
表 5 京津冀农牧系统氮素环境排放量和百分比
Table 5. Nitrogen environmental emissions of crop-livestock system in Beijing-Tianjin-Hebei region
项目 Item 河北 Hebei 天津 Tianjin 北京 Beijing 京津冀 Whole region (× 104 t) (%) (× 104 t) (%) (× 104 t) (%) (× 104 t) (%) 农田 Farmland 氨挥发 NH3 volatilization 34.9 12.5 1.9 11.1 1.1 10.3 37.8 12.3 氧化亚氮 N2O emission 3.4 1.2 0.2 1.1 0.1 1.0 3.7 1.2 淋溶 Leaching 25.0 9.0 1.3 7.5 0.7 6.9 27.0 8.8 径流和侵蚀 Runoff and erosion 6.3 2.3 0.3 2.0 0.2 1.8 6.8 2.2 反硝化 Denitrification 8.0 2.9 0.5 3.1 0.5 4.6 9.0 2.9 畜牧业 Livestock 氨挥发 NH3 volatilization 24.0 8.6 1.7 9.6 1.0 9.6 26.7 8.7 畜牧水体排放和堆置 Into water and stack-livestock 26.5 9.5 1.7 9.9 1.0 9.9 29.3 9.5 合计 Total 128.2 45.9 7.6 44.3 4.6 44.0 140.4 45.8 注(Note):百分比是指相应损失途径氮素损失量占系统总输入量的百分数 Percentage refers to the percentage of corresponding nitrogen loss to the total N input of the system. -
不同情景下京津冀地区氮素流动指标结果如表6所示,将化肥氮投入降低20%~40% (S1-N20~S3-N40),农牧系统总投入量也在降低,土壤累积量由当前情况的51.2万t逐渐降至8.1万t,降低幅度为84.2%,化肥投入继续降低至50%时 (S4-N50),土壤氮素出现亏缺现象 (–2.7%),从氮素投入平衡和不损耗土壤养分的角度看,化学氮肥可减施40% (S3-N40)。随氮素投入降低,氮素效率NUEc与NUEc+l逐渐增加,S3-N40情景下这两个效率分别达到54.7%和45.1%,较S0-cur分别增加34.9%和33.3%;食物氮代价和环境排放量逐渐降低,S3-N40情景下分别为2.2 kg/kg和106.8万t,分别较S0降低25.0%和23.9%。在S3-N40基础上,将水体排放与堆置畜禽粪便的1/2还田 (S5-N40+M) 以替代化学氮肥,土壤重新出现累积并且效率降低,即重新出现化肥减施潜力,可继续减施氮肥 (S6-N50+M和S7-N60+M)。其中S7-N60+M情景 (氮肥减施60%) 土壤累积也出现亏缺,因此在增加有机肥还田量替代氮肥情况下氮肥可减施50% (S6-N50+M),在此情景下,NUEc、NUEc+l继续增加至55.9%和49.2%,较S0增加37.7%和45.4%;食物氮代价和环境排放量继续降低至2.0 kg/kg和89.6万t,分别较S0-current降低31.2%和36.2%;氮素循环再利用率 (NRR) 也增加20%。因此,以平衡施氮原则减化学氮肥投入量时,减施潜力为40%,继续增加部分畜禽粪便还田量以替代化学氮肥 (有机替代) 时,减施潜力为50%。
表 6 京津冀农牧系统氮素流动情景分析
Table 6. Scenario analysis of nitrogen flow of crop-livestock system in Beijing-Tianjin-Hebei region
项目 Item S0-current S1-N20 S2-N30 S3-N40 S4-N50 S5-N40+M S6-N50+M S7-N60+M 化肥 Fertilizer (× 104 t) 191.6 153.3 134.1 115.0 95.8 115.0 95.8 76.6 畜禽粪便 Manure (× 104 t) 36.4 36.4 36.4 36.4 36.4 51.1 51.1 51.1 饲料外地输入 Feed import (× 104 t) 133.6 133.6 133.6 133.6 133.6 133.6 133.6 133.6 农牧系统总投入 Total input (× 104 t) 306.6 268.3 249.2 230.0 210.8 230.0 210.8 191.7 土壤累积 Soil accumulation (× 104 t) 51.2 29.6 18.9 8.1 –2.7 16.9 6.2 –4.6 NUEc (%) 40.6 46.6 50.3 54.7 60.0 51.3 55.9 61.3 NUEl (%) 25.0 25.0 25.0 25.0 25.0 25.0 25.0 25.0 NUEc+l (%) 33.8 38.6 41.6 45.1 49.2 45.1 49.2 54.1 NRR (%) 55.2 55.2 55.2 55.2 55.2 66.2 66.2 66.2 NCf (kg/kg) 3.0 2.6 2.4 2.2 2.0 2.2 2.0 1.8 环境排放 N emission (× 104 t) 140.4 123.6 115.2 106.8 98.4 98.0 89.6 81.2 环境排放占投入总量比例 (%)
Percentage of N emission45.8 46.1 46.2 46.4 46.7 42.6 42.5 42.4 注(Note):S0-cur—当前状况 Current situation; S1-N20—化肥氮减施 20% Reducing chemical N input by 20%; S2-N30—化肥氮减施 30% Reducing chemical N input by 30%; S3-N40—化肥氮减施 40% Reducing chemical N input by 40%; S4-N50—化肥氮减施 50% Reducing chemical N input by 50%; S5-N40+M—化肥氮减施 40% 增施有机肥 N 50%Reducing chemical N input by 40% and increasing manure N by 50%; S6-N50+M—化肥氮减施 50% 增施有机肥 N 50% Reducing chemical N input by 50% and increasing manure N by 50%; S7-N60+M—化肥氮减施 60% 增施有机肥 N50% Reducing chemical N input by 50% and increasing manure N by 50%. 此外各环境排放途径排放量和占比如图2所示,随化肥减施量的增加,各途径氮素排放量随之降低,S3-N40与S0-current相比,降低量最高的为农田氨挥发 (11.5万t),其次为农田淋溶 (10.8万t),降低量分别占其各自途径排放量的30.4%和40%,占环境排放总量的8.2%和7.7%;S6-N50+M与S0-current相比,降低量最高的途径为水体排放或堆置 (14.6万t),其次为农田淋溶 (13.5万t),降低量分别占其各自途径排放量的50.0%和50.0%,占环境排放总量的10.4%和9.6%。
-
京津冀地区农田系统氮素输入量为296.1万t,利用率为40.6%;畜牧系统氮素输入量为 133.6万t,利用率为25.0%;农牧系统整体氮素输入量为306.6万t,利用率为33.8%,其低于我国整体农牧系统氮素利用率 (41%)[32]。大量的化学氮肥投入,高外源饲料依赖率和低畜禽粪便还田量是造成利用低的主要原因[9],同时,这与京津冀地区人口氮需求总量大及养分投入量大也有很大关系。高量的氮肥投入及低效率,造成了氮素的环境排放量达140.4万t,占总氮素投入的45.8%,超出欧盟的30%的排放标准[33]。在环境排放中,氨挥发 (包括农田氨挥发与畜禽饲养阶段的氨挥发) 占较大比例,为21.0%,对总排放量贡献46%,这与马林[26]在全国尺度的估算接近,气态损失占总损失的52.6%。本研究中氮素输入项中未包含灌溉水及种子带入氮素,京津冀地区灌溉氮素和种子氮素约分别为4.73 万t和1.86 万t[34],若计入其中,氮素利用率将略降低 (33.1%)。
在三个地区中,北京地区农田系统氮素利用率低,畜牧系统利用率高,与该地区耕地面积小,农田化肥用量大,畜牧饲料外地输入量高,人口食物资源需求多有很大关系,其食物氮代价也最高 (表4)。河北地区资源量、需求量均最多,但其资源需求占京津冀地区比例与资源提供比例不匹配,承担了向京津地区输出食物的任务,起到了疏解首都生产压力的作用,在食物生产过程中产生的环境排放均留在了河北地区,造成河北地区环境排放率高。
-
在众多氮素损失途径中,农田氮素累积及畜禽粪便水体排放与堆置是最易调控的两个损失途径。并且,农田氮素累积 (盈余) 与氮素损失之间存在显著相关关系,有研究指出,氮素盈余与活性氮排放呈线性或指数关系[35],降低土壤累积是养分高效利用,降低环境损失的必要措施。本研究中土壤累积占总养分投入量的17.3%,占化肥投入量的26.7%,这引起了该地区很多环境问题,如平均土壤硝酸盐累积量超过150 kg/hm2[36],水体质量低于4级 (5级为最低等级)[37]等。在本研究中,应用输入输出平衡施氮原则,将土壤累积量降至较低情景时 (情景S3-N40),土壤累积量降低84.2%,农田氮素利用率和农牧系统氮素利用率分别增加14.1和11.3个百分点,环境排放降低23.9%,此时节氮潜力为40%。大量研究均指出,采取合理的养分管理措施,在不损失产量的情况下,节约氮肥30%~40%是完全可以实现的,如采用土壤-作物综合管理[34]、养分专家系统[38]、测土配方施肥技术[39]等措施。由于我国农牧分离的现状及畜牧养殖系统的粗犷管理,大量的畜禽粪便在清洗饲舍过程中直接排入水体中,也有部分堆置在养殖场内所含养分也均进入环境中[40]。在氮肥平衡施用的基础上增加1/2水体排放与堆置的畜禽粪便还田进入农田 (情景S6-N50+M) 时,农田系统和农牧系统氮素利用率分别增加15.3和15.4个百分点,环境排放降低36.2%,此时,节氮潜力为50%。畜禽粪便水体排放不仅造成养分资源浪费也严重威胁水体安全,有研究指出,畜禽粪便水体排放对流域水体氮和磷富营养化的贡献率分别为10%和20%[40]。农牧分离是农牧生产系统氮素利用率低下的重要原因[8,13],系统内部循环率低,造成资源的大量浪费与向环境排放[41]。在S6情景下,增加畜禽粪便还田量,增强了系统内部养分循环,氮素循环再利用率较当前状态增加20.0%,食物氮代价也降低了33.3%,既增加了氮肥利用效率,又减少资源消耗和环境排放。在此情景时,化肥投入与有机肥投入比例约为2∶1,相当于有机肥占比为33.3%,田昌玉等[42]也指出的有机肥比例在34%时可达到最大经济收益。因此,此情景是完全可以实现的。
平衡施氮,有机替代,农牧结合是解决当前我国氮肥过量施用问题的重要措施。欧美等发达国家已经提供了很多先进的经验[8,18],如美国实施的基于种养结合的畜禽粪便综合养分管理计划,通过玉米种植户和养殖场签订合同的方式实现农场间的农牧结合,从而提高氮素在农场间的循环利用;丹麦严格规定中小型畜禽养殖场将种植业和养殖业有机结合实现养分资源的重新利用,这些国家畜禽粪便还田率均高于50%,而我国还不足30%[6]。大量的有机资源排入环境中,造成大气污染及水体污染。合理利用有机资源,农牧结合,可提高氮素利用效率,降低环境损耗。
-
本研究是建立在统计数据基础上,以模型模拟方式进行的农牧系统氮素流动分析和化肥减施潜力探究,其还存在一些局限性,如本研究所用数据为统计数据,无法进行统计检验,但可以反映区域尺度的整体情况;本研究所用参数,大多数为华北地区通用参数,针对京津冀地区可能存在微小偏差;对于农牧结合的方式,未进行细致研究,在具体应用时还应与土壤消纳能力及资源的空间匹配方面结合。并且研究中所列其他氮素损失途径,也具有很大的节肥潜力,如畜禽粪便存贮阶段的氨挥发等[43-44],应进一步深入研究。
-
京津冀地区农牧系统总体氮素输入量为306.6万t,农田、畜牧、农牧系统氮素利用率分别为40.6%、25.0%和33.8%,环境排放量为140.4万t,占总输入量的45.8%,其中氨挥发为占比最大的排放途径。优化农田系统肥料氮素投入,使农田系统土壤氮素累积将为较低水平时,化学氮肥有40%的减施潜力,增加部分畜禽粪便还田量以替代化肥时,化学氮肥有50%的减施潜力。
基于Nufer模型的京津冀农牧系统氮素平衡状况与化学氮肥减施潜力分析
Nitrogen balance and potential reduction of nitrogen fertilizer input in Beijing-Tianjin-Hebei crop-livestock system based on NUFER model
-
摘要:
【目的】 为明确京津冀地区农牧系统氮素流动及平衡特征,解析化学氮肥减施潜力,从而为该地区化肥零增长和农业绿色发展提供依据。 【方法】 通过收集文献数据及2017年统计数据,采用Nufer (nutrient flows in food chain, environment and resources use) 模型,以“农田-畜牧”系统为研究对象,量化京津冀地区农牧系统氮素流动特征及利用效率,并针对输入输出平衡施氮和增加畜禽粪便还田量以替代化肥两种措施,通过情景分析评估该地区化学氮肥减施潜力。 【结果】 京津冀地区农田系统和畜牧系统氮素输入量分别为296.1和133.6万t,农牧系统氮素总体输入量为306.6万t,其中化学氮肥是最大的输入项,占总输入量的62.5%。农田系统、畜牧系统氮素利用率分别为40.6%和25.0%,农牧系统整体氮素利用率为33.8%,农牧系统氮素循环再利用率为55.2%,生产单位氮含量食物产品所需氮素投入量为3.0 kg/kg。该地区中,土壤氮素累积量为51.2万吨 (占总输入量的16.7%),环境排放总量为140.4万t (占总输入量的45.8%)。环境排放中,氨挥发 (包括农田和畜禽粪便氨挥发)、畜禽粪便水体排放与堆置为占比最大的两种损失途径,分别占总氮投入量的21.0%和9.5%。氮素土壤累积和畜禽粪便水体排放与堆置为两种最易调控的损失途径。通过输入输出平衡降低氮素土壤累积量,此地区有40%的氮肥减施潜力,此时,农牧系统氮素总投入量为230.0万t,农田和农牧系统氮素利用率分别较当前增加14.1和11.3个百分点;继续优化氮素管理,增加部分畜禽粪便还田量 (减少水体排放与堆置部分) 以替代化肥,此地区有50%的氮肥减施潜力,且在此状态下,农牧系统氮素总投入量为210.8万t,较当前降低31.2%,农田系统和农牧系统氮素利用率分别增加15.3和15.4个百分点,氮素环境排放量降低36.2%,氮素循环再利用率增加20.0%,食物氮代价降低33.3%。 【结论】 化学氮肥施用量大,农牧分离是京津冀地区农牧系统氮素利用效率低的主要原因。在农牧系统氮素管理中,平衡作物生产中的氮素投入,系统氮素投入的减少潜力有40%。如果能进一步合理利用有机资源,增加畜禽粪便还田率,化学氮肥减施的总潜力可在50%以上,且可提高氮素利用效率,有效降低氮素的环境排放总量。 Abstract:【Objectives】 This research aims to examine the nitrogen (N) flow and balance characteristics of crop-livestock system in Beijing-Tianjin-Hebei (BTH) region, and analyze the potential reduction of chemical N fertilizer application. This study could provide a theoretical basis for zero growth of chemical fertilizer inputs and the sustainable agricultural development in the region. 【Methods】 Using literature data and statistical data of 2017, the N flow characteristics and N-use efficiency (NUE) of crop-livestock system in the BTH region was quantified through NUFER (nutrient flows in food chain, environment and resources use) model. The reduction potential of chemical N fertilizer use in this region was evaluated by scenario analysis via two measures of N balance and manure substitution. 【Results】 The N input in crop and livestock systems of BTH region were 296.1 and 133.6 million tons, respectively. The total N input of the crop-livestock system in BTH region was 306.6 million tons and chemical N fertilizer was the highest contributor (62.5%). The NUE in crop and livestock systems were 40.6% and 25.0%, respectively, while the overall NUE of the crop-livestock system was 33.8%. The N recycling rate (NRR) of crop-livestock system was 55.2%, while cost of N for producing food containing a unit N (NCf) was 3.0 g/kg. The soil N accumulation was 51.2 × 104 t (16.7% of the total N input), and the total N emission was 140.4 × 104 t (45.8% of the total N input). The two main pathways of the total N emission include ammonia volatilization in field and livestock manure, and the N loss via water and livestock manure stacking, and they account for 21.0% and 9.5% of the total N input, respectively. These two pathways were mostly managed for N emission reduction. By balancing N input and output to decrease N accumulation in farmland, there would be a 40% potential reduction, total N input in crop-livestock system would decrease to 230.0 × 104 t, and NUEs of crop and crop-livestock systems would increase by 14.1 and 11.3 percentage points, respectively. By further optimizing N management via increase in the N input from livestock manure, which reduces N loss via water and manure stacking at the same time, the potential reduction could be as high as 50%. Thus, the total N input of crop-livestock system would be 31.2% (210.8 × 104 t) lower than the current level; the NUEs of crop and crop-livestock systems would increase by 15.3% and 15.4%, respectively; the total N emission would decrease by 36.2%; the N recycling rate would increase by 20.0%, and the NCf would decrease by 33.3%. 【Conclusions】 The heavy application of chemical N fertilizers and the separation of crop and livestock systems are the main causes of the relatively low NUE and high N emission in Beijing-Tianjin-Hebei region. By balancing the N input and output in crop system to reduce soil N accumulation, the potential reduction of N input in crop-livestock system would be 40%. Further optimization of N management by increasing the rate of returning livestock manure to the field would raise the potential reduction of N input to 50%, increase NUE, and decrease total N emission in the whole region. -
表 1 京津冀地区氮素资源量与需求量
Table 1. Nitrogen resources and demand in Beijing-Tianjin-Hebei region
项目
Item河北
Hebei天津
Tianjin北京
Beijing京津冀
Whole region资源量 Resources (× 104 t) 化肥 Fertilizer 177.3 9.2 5.1 191.6 秸秆 Straw 38.3 1.7 0.3 40.4 畜禽粪便 Manure 83.6 5.5 3.3 92.3 总计 Total 299.2 16.4 8.7 324.3 需求量 Demand (× 104 t) 植物生产 Plant production 151.3 7.5 1.8 160.5 动物生产 Animal production 36.2 3.1 2.0 41.3 人口消费 Human consumption 28.6 6.5 7.2 42.2 比例 Proportion (%) 资源 Resources 92.3 5.1 2.7 100.0 人口消费 Human consumption 67.8 15.4 17.1 100.0 注(Note):植物生产需求量是籽粒氮携出量与秸秆氮携出量之和,动物生产是动物产品氮携出量和动物其他部分氮携出量之和,人口消费需求量是人均消费量 (包括动物产品消费与植物产品消费) 乘以人口数量 Plant production demand refer to the sum of grain N and straw N. Animal production refer to the sum of animal product N and animal other parts N. Population consumption demand is the per capita consumption (including consumption of animal products and plant products) multiplied by the population; 资源比例与人口消费比例是指三地区各自占京津冀地区总资源量和人口消费量的比例 Proportions of resource and human consumption refer to the proportions of each region accounting for the total resources and human consumption of the whole region. 表 2 京津冀农田系统氮素流动与平衡
Table 2. Nitrogen flow and balance of cropping system in Beijing-Tianjin-Hebei region
项目
Item河北
Hebei天津
Tianjin北京
Beijing京津冀
Whole region输入 Input (× 104t) 化肥 Fertilizer 177.3 9.2 5.1 191.6 畜禽粪便 Manure 33.0 2.2 1.2 36.4 还田秸秆 Returned straw 18.0 0.8 0.2 18.9 沉降 Deposition 31.9 1.8 0.5 34.1 生物固氮 Biological N fixation 14.2 0.7 0.2 15.1 合计 Total 274.4 14.7 7.1 296.1 输出 Output (× 104t) 籽粒 Grain 112.9 5.7 1.5 120.1 秸秆 Straw 38.3 1.7 0.3 40.4 氨挥发 NH3 34.9 1.9 1.1 37.8 氧化亚氮 N2O 3.4 0.2 0.1 3.7 淋溶 Leaching 25.0 1.3 0.7 27.0 径流与侵蚀 Runoff and erosion 6.3 0.3 0.2 6.8 反硝化 Denitrification 8.0 0.5 0.5 9.0 土壤累积 Accumulation 45.5 3.0 2.7 51.2 合计 Total 274.4 14.7 7.1 296.1 农田系统氮素
利用效率 NUEc (%)41.2 38.9 20.7 40.6 表 3 京津冀畜牧系统氮素养分流动与平衡
Table 3. Nitrogen flow and balance of livestock system in Beijing-Tianjin-Hebei region
项目
Item河北
Hebei天津
Tianjin北京
Beijing京津冀
Whole region输入 Input (× 104 t) 籽粒 Grain 47.1 2.3 0.4 49.9 秸秆 Straw 17.0 0.8 0.1 17.9 饲料外地输入 Feed import 55.6 5.5 4.7 65.8 合计 Total 119.7 8.6 5.3 133.6 输出 Output (× 104 t) 动物产品 Animal products 29.1 2.5 1.8 33.4 其他部分 Other parts 7.1 0.5 0.3 7.9 畜禽粪便总量 Total manure 83.6 5.5 3.3 92.3 氨挥发 NH3 volatilization 24.0 1.7 1.0 26.7 水体排放与堆置
via water and manure stacking26.5 1.7 1.0 29.3 还田 Returning to field 33.0 2.2 1.2 36.4 合计 Total 119.7 8.6 5.3 133.6 畜牧系统氮素
利用效率 NUEl (%)24.3 29.7 33.4 25.0 注(Note):动物产品是指肉、蛋、奶等主产品 Animal products refer to main products such as meat, eggs and milk; 其他部分是指骨头、皮毛等不可食用部分 other parts refer to inedible parts such as bones and fur. 表 4 京津冀农牧系统氮素养分流动与平衡
Table 4. Nitrogen flow and balance of crop-livestock system in Beijing-Tianjin-Hebei region
项目 Item 河北 Hebei 天津 Tianjin 北京 Beijing 京津冀 Whole region 输入 Input (× 104t) 化肥 Fertilizer 177.3 9.2 5.1 191.6 沉降 Deposition 31.9 1.8 0.5 34.1 生物固氮 Biological N fixation 14.2 0.7 0.2 15.1 饲料外地输入 Feed import 55.6 5.5 4.7 65.8 合计 Total 279.0 17.2 10.4 306.6 输出 Output (× 104t) 植物产品 Plant products 65.8 3.4 1.0 70.2 动物产品 Animal products 29.1 2.5 1.8 33.4 土壤累积 Accumulation 45.5 3.0 2.7 51.2 农田氨挥发 NH3-farmland 34.9 1.9 1.1 37.8 农田氧化亚氮 N2O-farmland 3.4 0.2 0.1 3.7 农田淋溶 Leaching-farmland 25.0 1.3 0.7 27.0 农田径流和侵蚀 Runoff and erosion-farmland 6.3 0.3 0.2 6.8 农田反硝化 Denitrification-farmland 8.0 0.5 0.5 9.0 畜禽粪便氨挥发 NH3-livestock manure 24.0 1.7 1.0 26.7 畜禽粪便水体排放与堆置 Into water and stack- livestock manure 26.5 1.7 1.0 29.3 其他部分 Other parts 10.5 0.7 0.3 11.4 合计 Total 279.0 17.2 10.4 306.6 农牧系统氮素利用效率 NUEc+l (%) 34.0 34.6 26.8 33.8 农牧系统氮素循环再利用率 NRR (%) 55.8 51.7 42.9 55.2 食物氮代价 NCf (kg/kg) 2.9 2.9 3.7 3.0 注(Note):其余部分包括植物生产的籽粒和秸秆在存放和运输过程的损失,和动物生产中的动物骨头、皮毛等部分 The other parts include the loss of seeds and straw produced by plants during storage and transportation, and animal bones and fur in animal production process. 表 5 京津冀农牧系统氮素环境排放量和百分比
Table 5. Nitrogen environmental emissions of crop-livestock system in Beijing-Tianjin-Hebei region
项目 Item 河北 Hebei 天津 Tianjin 北京 Beijing 京津冀 Whole region (× 104 t) (%) (× 104 t) (%) (× 104 t) (%) (× 104 t) (%) 农田 Farmland 氨挥发 NH3 volatilization 34.9 12.5 1.9 11.1 1.1 10.3 37.8 12.3 氧化亚氮 N2O emission 3.4 1.2 0.2 1.1 0.1 1.0 3.7 1.2 淋溶 Leaching 25.0 9.0 1.3 7.5 0.7 6.9 27.0 8.8 径流和侵蚀 Runoff and erosion 6.3 2.3 0.3 2.0 0.2 1.8 6.8 2.2 反硝化 Denitrification 8.0 2.9 0.5 3.1 0.5 4.6 9.0 2.9 畜牧业 Livestock 氨挥发 NH3 volatilization 24.0 8.6 1.7 9.6 1.0 9.6 26.7 8.7 畜牧水体排放和堆置 Into water and stack-livestock 26.5 9.5 1.7 9.9 1.0 9.9 29.3 9.5 合计 Total 128.2 45.9 7.6 44.3 4.6 44.0 140.4 45.8 注(Note):百分比是指相应损失途径氮素损失量占系统总输入量的百分数 Percentage refers to the percentage of corresponding nitrogen loss to the total N input of the system. 表 6 京津冀农牧系统氮素流动情景分析
Table 6. Scenario analysis of nitrogen flow of crop-livestock system in Beijing-Tianjin-Hebei region
项目 Item S0-current S1-N20 S2-N30 S3-N40 S4-N50 S5-N40+M S6-N50+M S7-N60+M 化肥 Fertilizer (× 104 t) 191.6 153.3 134.1 115.0 95.8 115.0 95.8 76.6 畜禽粪便 Manure (× 104 t) 36.4 36.4 36.4 36.4 36.4 51.1 51.1 51.1 饲料外地输入 Feed import (× 104 t) 133.6 133.6 133.6 133.6 133.6 133.6 133.6 133.6 农牧系统总投入 Total input (× 104 t) 306.6 268.3 249.2 230.0 210.8 230.0 210.8 191.7 土壤累积 Soil accumulation (× 104 t) 51.2 29.6 18.9 8.1 –2.7 16.9 6.2 –4.6 NUEc (%) 40.6 46.6 50.3 54.7 60.0 51.3 55.9 61.3 NUEl (%) 25.0 25.0 25.0 25.0 25.0 25.0 25.0 25.0 NUEc+l (%) 33.8 38.6 41.6 45.1 49.2 45.1 49.2 54.1 NRR (%) 55.2 55.2 55.2 55.2 55.2 66.2 66.2 66.2 NCf (kg/kg) 3.0 2.6 2.4 2.2 2.0 2.2 2.0 1.8 环境排放 N emission (× 104 t) 140.4 123.6 115.2 106.8 98.4 98.0 89.6 81.2 环境排放占投入总量比例 (%)
Percentage of N emission45.8 46.1 46.2 46.4 46.7 42.6 42.5 42.4 注(Note):S0-cur—当前状况 Current situation; S1-N20—化肥氮减施 20% Reducing chemical N input by 20%; S2-N30—化肥氮减施 30% Reducing chemical N input by 30%; S3-N40—化肥氮减施 40% Reducing chemical N input by 40%; S4-N50—化肥氮减施 50% Reducing chemical N input by 50%; S5-N40+M—化肥氮减施 40% 增施有机肥 N 50%Reducing chemical N input by 40% and increasing manure N by 50%; S6-N50+M—化肥氮减施 50% 增施有机肥 N 50% Reducing chemical N input by 50% and increasing manure N by 50%; S7-N60+M—化肥氮减施 60% 增施有机肥 N50% Reducing chemical N input by 50% and increasing manure N by 50%. -
[1] 张福锁, 王激清, 张卫峰, 等. 中国主要粮食作物肥料利用率现状与提高途径[J]. 土壤学报, 2008, 45(5): 148–157. Zhang F S, Wang J Q, Zhang W F, et al. Nutrient use efficiencies of major cereal crops in China and measures for improvement[J]. Acta Pedologica Sinica, 2008, 45(5): 148–157. [2] 李欠欠, 李雨繁, 高强, 等. 传统和优化施氮对春玉米产量、氨挥发及氮平衡的影响[J]. 植物营养与肥料学报, 2015, 21(3): 571–579. Li Q Q, Li Y F, Gao Q, et al. Effect of conventional and optimized nitrogen fertilization on spring maize yield, ammonia volatilization and nitrogen balance in soil-maize system[J]. Journal of Plant Nutrition and Fertilizers, 2015, 21(3): 571–579. doi: 10.11674/zwyf.2015.0303 [3] 薛利红, 何世颖, 段婧婧, 等. 基于养分回用-化肥替代的农业面源污染氮负荷削减策略及技术[J]. 农业环境科学学报, 2017, 36(7): 1226–1231. Xue L H, He S Y, Duan J J, et al. Agricultural non-point source pollution source nitrogen load reduction strategy and technology of nutrient reusing in agricultural fields to replace fertilizer[J]. Journal of Agro-Environment Science, 2017, 36(7): 1226–1231. doi: 10.11654/jaes.2017-0949 [4] 熊正琴, 张晓旭. 氮肥高效施用在低碳农业中的关键作用[J]. 植物营养与肥料学报, 2017, 23(6): 1433–1440. Xiong Z Q, Zhang X X. Key role of efficient nitrogen application in low carbon agriculture[J]. Journal of Plant Nutrition and Fertilizers, 2017, 23(6): 1433–1440. doi: 10.11674/zwyf.17282 [5] 胡春胜, 张玉铭, 秦树平, 等. 华北平原农田生态系统氮素过程及其环境效应研究[J]. 中国生态农业学报, 2018, 26(10): 1501–1514. Hu C S, Zhang Y M, Qin S P, et al. Nitrogen processes and related environmental effects on agro-ecosystem in the North China Plain[J]. Chinese Journal of Eco-Agriculture, 2018, 26(10): 1501–1514. [6] 牛新胜, 巨晓棠. 我国有机肥料资源及利用[J]. 植物营养与肥料学报, 2017, 23(6): 1462–1479. Niu X S, Ju X T. Organic fertilizer resources and utilization in China[J]. Journal of Plant Nutrition and Fertilizers, 2017, 23(6): 1462–1479. doi: 10.11674/zwyf.17430 [7] 李书田, 刘晓永, 何萍. 当前我国农业生产中的养分需求分析[J]. 植物营养与肥料学报, 2017, 23(6): 1416–1432. Li S T, Liu X Y, He P. Analyses on nutrient requirements in current agriculture production in China[J]. Journal of Plant Nutrition and Fertilizers, 2017, 23(6): 1416–1432. doi: 10.11674/zwyf.17393 [8] Ma L, Bai Z, Ma W, et al. Exploring future food provision scenarios for China[J]. Environmental Science & Technology, 2019, 53(3): 1385–139. [9] 佟丙辛, 张华芳, 高肖贤, 等. 华北平原典型区域农牧系统氮素流动及其环境效应—以河北省为例[J]. 中国农业科学, 2018, 51(3): 442–455. Tong B X, Zhang H F, Gao X X, et al. Nitrogen flow and environmental effects of crop–livestock system in typical area of North China Plain—A case study in Hebei Province[J]. Scientia Agricultura Sinica, 2018, 51(3): 442–455. doi: 10.3864/j.issn.0578-1752.2018.03.005 [10] Leip A, Marchi G, Koeble R, et al. Linking an economic model for European agriculture with a mechanistic model to estimate nitrogen and carbon losses from arable soils in Europe[J]. Biogeosciences, 2008, 5: 73–94. doi: 10.5194/bg-5-73-2008 [11] Isermann K, Isermann R. Food production and consumption in Germany: N flows and N emissions[J]. Nutrient Cycling in Agroecosystems, 1998, 52(2–3): 289–301. [12] Ma L, Ma W, Velthof G, et al. Modeling nutrient flows in the food chain of China[J]. Journal of Environmental Quality, 2010, 39(4): 1279–1289. doi: 10.2134/jeq2009.0403 [13] 张建杰, 郭彩霞, 李莲芬, 等. 农牧交错带农牧系统氮素流动与环境效应—以山西省为例[J]. 中国农业科学, 2018, 51(3): 456–467. Zhang J J, Guo C X, Li L F, et al. Nutrient flow and environmental effects on crop-livestock system in farming–pastoral transition zone —A case study in Shanxi Province[J]. Scientia Agricultura Sinica, 2018, 51(3): 456–467. doi: 10.3864/j.issn.0578-1752.2018.03.006 [14] 张晓萌, 王寅, 焉莉, 等. 东北地区农牧系统氮、磷养分流动特征[J]. 中国农业科学, 2018, 51(3): 417–429. Zhang X M, Wang Y, Yan L, et al. Characteristics of nitrogen and phosphorus nutrient flow in farming and animal husbandry system in northeast China[J]. Scientia Agricultura Sinica, 2018, 51(3): 417–429. doi: 10.3864/j.issn.0578-1752.2018.03.003 [15] 米国华, 伍大利, 陈延玲, 等. 东北玉米化肥减施增效技术途径探讨[J]. 中国农业科学, 2018, 51(14): 2758–2770. Mi G H, Wu D L, Chen Y L, et al. The ways to reduce chemical fertilizer input and increase fertilizer use efficiency in maize in Northeast China[J]. Scientia Agricultura Sinica, 2018, 51(14): 2758–2770. doi: 10.3864/j.issn.0578-1752.2018.14.013 [16] 赵亚南, 徐霞, 黄玉芳, 等. 河南省小麦、玉米氮肥需求及节氮潜力[J]. 中国农业科学, 2018, 51(14): 2747–2757. Zhao Y N, Xu X, Huang Y F, et al. Nitrogen requirement and saving potential for wheat and maize in Henan Province[J]. Scientia Agricultura Sinica, 2018, 51(14): 2747–2757. doi: 10.3864/j.issn.0578-1752.2018.14.012 [17] 宋大利, 侯胜鹏, 王秀斌, 等. 中国畜禽粪尿中养分资源数量及利用潜力[J]. 植物营养与肥料学报, 2018, 24(5): 1131–1148. Song D L, Hou S P, Wang X B, et al. Nutrient resource quantity of animal manure and its utilization potential in China[J]. Journal of Plant Nutrition and Fertilizers, 2018, 24(5): 1131–1148. doi: 10.11674/zwyf.17415 [18] 耿维, 孙义祥, 袁嫚嫚, 等. 安徽省畜牧业环境承载力及粪便替代化肥潜力评估[J]. 农业工程学报, 2018, 34(18): 252–260, 315. Geng W, Sun Y X, Yuan M M, et al. Evaluation of environmental carrying capacity and manure organic fertilizer instead chemical fertilizer potential of animal husbandry in Anhui[J]. Transactions of the Chinese Society of Agricultural Engineering, 2018, 34(18): 252–260, 315. doi: 10.11975/j.issn.1002-6819.2018.18.031 [19] Jin X, Bai Z, Oenema O, et al. Spatial planning needed to drastically reduce nitrogen and phosphorus surpluses in China’s agriculture[J]. Environmental Science & Technology, 2020, 54(19): 11894–11904. [20] 冼超凡, 王莉雁, 逯非, 等. 基于“氮补偿”机制及数据包络分析的活性氮减排评价——以京津冀区域为例[J]. 应用生态学报, 2017, 28(8): 2545–2553. Xian C F, Wang L Y, Lu F, et al. Evaluation of reactive nitrogen emissions reduction in Beijing-Tianjin-Hebei agglomeration based on "N offset" mechanism and data envelopment analysis[J]. Chinese Journal of Applied Ecology, 2017, 28(8): 2545–2553. [21] 凌文翠, 范玉梅, 方瑶瑶, 等. 京津冀地区畜禽养殖业抗生素污染现状分析[J]. 环境工程技术学报, 2018, 8(4): 390–397. Ling W C, Fan Y M, Fang Y Y, et al. Antibiotics pollution of livestock and poultry breeding in Beijing-Tianjin-Hebei region[J]. Journal of Environmental Engineering Technology, 2018, 8(4): 390–397. [22] 李金祥. 科技助推京津冀地区畜禽养殖废弃物处理与利用[J]. 北方牧业, 2018, (3): 10. Li J X. Science and technology promote the treatment and utilization of livestock and poultry breeding waste in Beijing-Tianjin-Hebei region[J]. North of China livestock husbandry, 2018, (3): 10. [23] 张卫峰, 李亮科, 陈新平, 等. 我国复合肥发展现状及存在的问题[J]. 磷肥与复肥, 2009, 24(2): 14–16. Zhang W F, Li L K, Chen X P, et al. The present status and existing problems in China's compound fertilizer development[J]. Phosphate & Compound Fertilizer, 2009, 24(2): 14–16. doi: 10.3969/j.issn.1007-6220.2009.02.005 [24] 尹兴, 张丽娟, 刘学军, 等. 河北平原城市近郊农田大气氮沉降特征[J]. 中国农业科学, 2017, 50(4): 698–710. Yin X, Zhang L J, Liu X J, et al. Nitrogen deposition in suburban croplands of Hebei Plain[J]. Scientia Agricultura Sinica, 2017, 50(4): 698–710. doi: 10.3864/j.issn.0578-1752.2017.04.010 [25] 李贞宇, 王旭, 高志岭, 等. 我国不同区域小麦施肥资源环境影响的生命周期评价[J]. 农业环境科学学报, 2010, 29(7): 1417–1422. Li Z Y, Wang X, Gao Z L, et al. Life cycle assessment of wheat fertilization in different regions of China[J]. Journal of Agro-Environment Science, 2010, 29(7): 1417–1422. [26] 马林. 中国食物链氮素流动规律及调控策略[D]. 河北保定: 河北农业大学博士学位论文, 2010. Ma L. Mechanism and regulatory strategies of nitrogen flow in food chain of China[D]. Baoding, Hebei: PhD Dissertation of Agricultural University of Hebei, 2010. [27] 高利伟, 马林, 张卫峰, 等. 中国作物秸秆养分资源数量估算及其利用状况[J]. 农业工程学报, 2009, 25(7): 173–179. Gao L W, Ma L, Zhang W F, et al. Estimation of nutrient resource quantity of crop straw and its utilization situation in China[J]. Transactions of the Chinese Society of Agricultural Engineering, 2009, 25(7): 173–179. doi: 10.3969/j.issn.1002-6819.2009.07.032 [28] 侯勇, 高志岭, 马文奇, 等. 京郊典型集约化“农田-畜牧”生产系统氮素流动特征[J]. 生态学报, 2012, 32(4): 24–32. Hou Y, Gao Z L, Ma W Q, et al. Nitrogen flows in intensive “crop–livestock” production systems typically for the peri–urban area of Beijing[J]. Acta Ecologica Sinica, 2012, 32(4): 24–32. [29] 马林, 马文奇, 张福锁, 等. 中国食物链养分流动与管理研究[J]. 中国生态农业学报, 2018, 26(10): 1494–1500. Ma L, Ma W Q, Zhang F S, et al. Nutrient flow and management in the food chain in China[J]. Chinese Journal of Eco-Agriculture, 2018, 26(10): 1494–1500. [30] 巨晓棠. 理论施氮量的改进及验证—兼论确定作物氮肥推荐量的方法[J]. 土壤学报, 2015, 52(2): 249–261. Ju X T. Improvement and validation of theoretical N rate (TNR) —Discussing the methods for N fertilizer recommendation[J]. Acta Pedologica Sinica, 2015, 52(2): 249–261. [31] Haas G, Caspari B, Kopke U. Nutrient cycling in organic farms: stall balance of a suckler cow herd and beef bulls[J]. Nutrient Cycling in Agroecosystems, 2002, 64(3): 225–230. doi: 10.1023/A:1021442806311 [32] 刘晓利. 我国“农田–畜牧–营养–环境”体系氮素养分循环与平衡[D]. 河北保定: 河北农业大学硕士学位论文, 2005. Liu X L. Cycle and balance of nitrogen nutrients in my country's "farmland–livestock–nutrition–environment" system[D]. Baoding, Hebei: MS thesis of Agricultural University of Hebei, 2005. [33] Oenema O, Oudendag O, Velthof G. Nutrient losses from manure management in the European Union[J]. Livestock Science, 2007, 112: 261–272. doi: 10.1016/j.livsci.2007.09.007 [34] He W, Jiang R, He P, et al. Estimating soil nitrogen balance at regional scale in China’s croplands from 1984 to 2014[J]. Agricultural Systems, 2018, 167: 125–135. doi: 10.1016/j.agsy.2018.09.002 [35] Chen X, Cui Z, Fan M, et al. Producing more grain with lower environmental costs[J]. Nature, 2014, 514: 486–491. doi: 10.1038/nature13609 [36] Lu J, Bai Z, Velthof G, et al. Accumulation and leaching of nitrate in soils in wheat–maize production in China[J]. Agricultural Water Management, 2019, 212: 407–415. doi: 10.1016/j.agwat.2018.08.039 [37] Bai Z, Lu J, Zhao H, et al. Designing vulnerable zones of nitrogen and phosphorus transfers to control water pollution in China[J]. Environmental Science & Technology, 2018, 52: 8987–8988. [38] Zhang J, He P, Xu X, et al. Nutrient Expert improves nitrogen efficiency and environmental benefits for winter wheat in China[J]. Agronomy Journal, 2018, 110: 696–706. doi: 10.2134/agronj2017.05.0291 [39] 张卫红, 李玉娥, 秦晓波, 等. 应用生命周期法评价我国测土配方施肥项目减排效果[J]. 农业环境科学学报, 2015, 34(7): 1422–1428. Zhang W H, Li Y E, Qin X B, et al. Evaluation of greenhouse gas emission reduction by balanced fertilization in China using life cycle assessment[J]. Journal of Agro-Environment Science, 2015, 34(7): 1422–1428. doi: 10.11654/jaes.2015.07.027 [40] 王磊, 香宝, 苏本营, 等. 京津冀地区农业面源污染风险时空差异研究[J]. 农业环境科学学报, 2017, 36(7): 1254–1265. Wang L, Xiang B, Su B Y, et al. Spatial-temporal variation of agricultural non-point source pollution risk in Beijing-Tianjin-Hebei region, China[J]. Journal of Agro-Environment Science, 2017, 36(7): 1254–1265. doi: 10.11654/jaes.2016-1539 [41] 陈轩敬, 宫雅慧, 谢军, 等. 近20年重庆市农牧生产体系氮素流动特征及驱动力分析[J]. 中国农业科学, 2018, 51(3): 468–480. Chen X J, Gong Y H, Xie J, et al. Nitrogen flow of crop-livestock production system and its driving forces in Chongqing over the past 20 years[J]. Scientia Agricultura Sinica, 2018, 51(3): 468–480. doi: 10.3864/j.issn.0578-1752.2018.03.007 [42] 田昌玉, 林治安, 唐继伟, 等. 基于最佳经济效益的冬小麦–夏玉米轮作体系有机肥氮替代率的长期演变[J]. 植物营养与肥料学报, 2019, 25(10): 1623–1632. Tian C Y, Lin Z A, Tang J W, et al. Evolution of optimum substitution rate of chemical nitrogen by organic nitrogen for the optimum nitrogen for the optimum economic benefit in long-term winter wheat-summer maize system[J]. Journal of Plant Nutrition and Fertilizers, 2019, 25(10): 1623–1632. doi: 10.11674/zwyf.18473 [43] 曹玉博, 邢晓旭, 柏兆海, 等. 农牧系统氨挥发减排技术研究进展[J]. 中国农业科学, 2018, 51(3): 566–580. Cao Y B, Xing X X, Bai X H, et al. Review on ammonia emission mitigation techniques of crop-livestock production system[J]. Scientia Agricultura Sinica, 2018, 51(3): 566–580. doi: 10.3864/j.issn.0578-1752.2018.03.014 [44] Wei S, Bai Z, Qin W, et al. Environmental, economic and social analysis of peri-urban pig production[J]. Journal of Cleaner Production, 2016, 129: 596–607. doi: 10.1016/j.jclepro.2016.03.133 -