中国栽培草地氮磷流动空间特征

魏志标，柏兆海，马林，张福锁



中国栽培草地氮磷流动空间特征

魏志标1,2，柏兆海2，马林2，张福锁1

（1中国农业大学资源与环境学院/植物-土壤相互作用教育部重点实验室，北京 100193；2中国科学院遗传与发育生物学研究所农业资源研究中心/ 河北省节水农业重点实验室/中国科学院农业水资源重点实验室，石家庄 050021）

【目的】中国畜牧业快速发展，对牧草需求日益增加，种植业结构调改已成为必然趋势。定量研究栽培草地的氮磷流动特征，可为优化牧草施肥和提高牧草产量提供科学依据。【方法】以苜蓿、黑麦草和燕麦草为研究对象，建立中国主要栽培牧草氮磷养分输入(输出)数据库，利用NUFER模型定量中国主要栽培牧草氮磷平衡账户、利用率和环境排放特征。【结果】（1）2014年，全国苜蓿、黑麦草和燕麦草草地氮输入(输出)总量分别为1 547、236和67 Gg，磷输入(输出)总量分别为323、44和16 Gg。单位面积苜蓿、黑麦草和燕麦草草地的氮输入(输出)量分别为326、427和217 kg N·hm-2，磷输入(输出)量分别为49、18和35 kg P·hm-2。生物固氮是苜蓿草地氮的最主要输入项，占总输入量的51%，氮肥占黑麦草草地氮总输入量的93%，占燕麦草地总输入量的84%。磷肥是3种牧草磷的最主要输入项。（2）2014年，全国苜蓿、黑麦草和燕麦草的氮利用率分别为64%、93%和69%；磷利用率分别为28%、77%和34%。不同区域氮磷利用率差异很大。(3)2014年，单位面积苜蓿、黑麦草和燕麦草草地的氮环境排放量分别为23、4.0和9.9 kg N·hm-2，磷的环境排放量分别为2.6、3.8和2.6 kg P·hm-2。西南地区苜蓿和黑麦草的氮磷环境排放量较高，西部各地区燕麦草草地氮磷环境排放量高于其他区域。(4)2014年，黑麦草草地氮呈亏缺态，其余栽培草地氮磷均有不同程度的盈余。中国西部地区苜蓿的土壤氮磷累积量高于东部地区；黑麦草的土壤氮磷累积量没有明显规律；青藏高原区燕麦草的土壤氮磷累积量大于其他区域。【结论】2014年，全国栽培苜蓿、黑麦草和燕麦草草地的氮磷输入(输出)量差异较大。生物固氮是苜蓿获取氮素的重要途径，田间管理需减少氮肥投入，燕麦和黑麦草主要通过化肥获取所需养分，应注重肥料施用。3种栽培牧草氮的利用率均高于60%，磷的利用率相对较低。黑麦草草地氮呈亏缺态，需增施氮肥，其余栽培草地氮磷均有不同程度的盈余，应控制肥料过多施用。

栽培草地；氮；磷；养分流动；养分利用率；土壤养分累积

0 引言

【研究意义】栽培草地的种植规模和水平是衡量一个国家畜牧业发展程度的标准[1]。过去20年间，中国栽培草地面积增加了两倍以上[2-3]。中国主要栽培牧草超过70种，几乎在全国各地均有分布，面积约2.2×107hm2，分为多年生和一年生两种类型[4]。苜蓿是中国最主要栽培牧草，占多年生牧草总量的28%，燕麦和黑麦草是最主要的一年生牧草，分别占一年生牧草总量的50%和12%。2016年中央一号文件提出要扩大粮改饲试点，加快建设现代饲草产业体系[5]。集约化管理的栽培草地生产力是天然草地生产力的10—20倍[6-7]。通过建设栽培草地，可以将大部分的天然草地从目前的放牧压力中解放出来，充分发挥其生态功能[2,8]。【前人研究进展】同美国等畜牧业发达的国家相比，不合理的田间管理措施导致中国的牧草产量偏低[4,9]。播种、施肥和灌溉作为重要的栽培管理措施，会显著影响牧草的产量[10-12]。不同地区需要因地制宜的选择适宜的节水灌溉方式。有研究表明，相比畦灌和不灌溉，西北地区采用滴灌的方式牧草可以增产12%和108%，水分利用率比畦灌提高19%[13]。关于养分管理方面的文章大多是围绕着某一特定的地区或者牧草种类做氮磷肥效试验，或者是研究牧区土壤的基础肥力。陈萍等[14]分析了碱性条件下苜蓿的优化施肥方案，指出苜蓿的最佳产量施肥量为氮肥60 kg·hm-2，磷肥135 kg·hm-2；肖相芬等分别报道了黑麦草和燕麦草的优化施肥量[15-16]；顾梦鹤等[17]在青藏高原的研究表明，在人工建植的混播草地上，施肥可以有效抑制杂草的入侵和生长，有利于青藏高原高寒地区的经济发展和生态环境建设。关于中国主要草地土壤养分累积量的文章也有部分报道[18]，Gu等[19]认为中国草地氮素有部分盈余，氮的利用率为59%；Zhou等[20]综述了全球草地磷流动的结果，认为栽培草地磷累积量分布在-46—126 kg·hm-2的范围内。总体而言，中国栽培草地的氮磷养分管理技术规范尚不健全，氮磷肥施用不合理，不同区域间肥料的使用量差异很大[9]，迫切需要对草地养分流动情况有更深入的了解。【本研究切入点】苜蓿、黑麦草和燕麦草是中国栽培牧草的典型种类，分别属于多年生牧草、一年生草本和一年生饲用作物[4]，种植面积较广，三者的养分流动情况可以反映中国牧草生产的总体状况。本研究通过建立中国不同区域苜蓿、黑麦草和燕麦草氮磷输入和输出的数据库，使用NUFER模型，估算牧草生产过程中的环境排放。【拟解决的关键问题】定量中国不同区域天然草地的养分平衡账户；分析不同区域天然草地氮磷利用率；评估不同区域天然草地的氮磷环境排放。

1 材料与方法

1.1 数据来源

本研究利用NUFER模型计算，模型包括输入模块、输出模块和计算模块[21-22]。输入模块包括变量和参数。化肥用量、牧草播种面积、作物产量等属于模型变量，来源于统计数据和专家咨询等。本研究先后咨询了中国农业大学、河北省农林科学院、西藏农牧科学研究院、山西省农业科学院、西北农林科技大学、西南林业大学等单位的相关专家，获取了内蒙、新疆、西藏、云南、四川等地区牧草的氮磷肥料施用量。参数例如牧草氮磷含量、养分环境排放因子来源于文献数据。模型输出模块包括氮磷平衡帐户、氮磷利用率和各种形态氮磷环境排放，包括挥发（NH3）、氧化亚氮排放（N2O）、反硝化（N2）、氮磷径流、侵蚀、淋洗等。

1.2 模型计算

养分账户的输入项包括粪尿输入、大气沉降、生物固氮（N）、灌溉带入、播种带入；输出项包括牧草收获、挥发（NH3）、氧化亚氮排放（N2O）、反硝化（N2）、径流、侵蚀、淋洗、累积。具体的养分输入和输出公式如下：

Iftotal=Iffertilizer+Ifdeposition+Iffixation+Ifirrigation+Ifseeding（1）

Ifseeding=SR×SeedN(P)C/100 （2）

式中，Iftotal代表栽培草地氮（磷）总输入量（kg·hm-2），Iffertilizer代表栽培草地中氮（磷）肥施用量（kg·hm-2），其数据基于专家咨询和文献搜集获得，各区域的具体量见表1；Ifdeposition代表氮（磷）沉降（kg·hm-2），Iffixation代表栽培草地生物固氮（kg·hm-2），Ifirrigation代表灌溉带入栽培草地中的氮（磷）（kg·hm-2），以上数据通过文献总结获得[23-27]，具体数据见表2；Ifseeding代表播种带入栽培草地中的氮（磷）（kg·hm-2），其数据通过栽培牧草播种量和种子氮（磷）含量估算；栽培牧草播种量（kg·hm-2）和种子氮（磷）含量（%）分别由SR和SeedN(P)C表示，通过文献总结获得，二者的具体数据见表3。

Oftotal=Offorage+OfNH3+OfN2O+Ofdenitrification+Ofrunoff+Oferosion+Ofleaching+Ofaccumulation（3）

Offorage=yield×forageN(P)C/100 （4）

式中，Oftotal代表栽培草地氮（磷）总输出量（kg·hm-2），Offorage表示栽培牧草氮（磷）输出量（kg·hm-2），其数值基于牧草产量和牧草氮（磷）含量估算；牧草产量（kg·hm-2）和栽培牧草氮（磷）含量（%）分别由yield和forageN(P)C表示，二者的数值分别源于《中国草业统计2014》和已发表的文献，各地区的具体量见表3和表4；OfNH3表示氨挥发，OfN2O表示氧化亚氮排放，Ofdenitrification表示反硝化，Ofrunoff表示氮（磷）径流，Oferosion表示氮（磷）侵蚀，Ofleaching表示氮（磷）淋洗，Ofaccumulation表示栽培草地土壤氮（磷）累积（kg·hm-2），以上数据由NUFER模型模拟得到[21-22]，该模型可用于模拟食物链养分流动情况。

N(P)UEf=（Offorage/Iftotal）×100% （5）

式中，N(P)UEf表示栽培草地氮（磷）利用率，N(P)UEf由牧草氮（磷）输出量和栽培草地氮（磷）输入量的比值得到。

2 结果

2.1 2014年全国主要栽培草地氮磷流动特征

2014年全国苜蓿草地氮总输入（输出）量为1 547 Gg，生物固氮和氮肥输入到苜蓿草地中的氮占总输入量的51%和37%，另有12%通过氮沉降、播种和灌溉3种方式输入到苜蓿草地。牧草携出占苜蓿草地氮输出量的63%，苜蓿草地土壤氮累积量为294 Gg。中国苜蓿草地磷的总输入量为323 Gg，磷肥输入占总输入量的99%，苜蓿草地磷总输入量的68%累积到土壤中。2014年单位面积苜蓿草地氮输入（输出）量为326 kg N·hm-2，氮的总损失量为55 kg N·hm-2，氨挥发、氧化亚氮排放、径流、侵蚀、淋洗和反硝化6种途径的损失量分别为32、1.3、2.0、0.29、7.7和12 kg N·hm-2，苜蓿草地氮利用率为64%。全国苜蓿草地土壤氮累积量为62 kg N·hm-2。全国单位面积苜蓿草地磷输入量为68 kg P·hm-2，沉降、灌溉和播种3种途径输入到土壤中磷的总和不足0.50 kg P·hm-2，相比磷肥的施入量，这3项值可以忽略不计。苜蓿带走的磷为19 kg P·hm-2，苜蓿草地磷的利用率为28%。通过径流、淋洗和侵蚀损失到环境中的磷量分别为0.26、2.2和0.10 kg P·hm-2，全国苜蓿草地土壤磷的累积量很大，为47 kg P·hm-2（表5）。

2014年全国黑麦草草地氮总输入（输出）量为236 Gg，氮肥是最主要输入途径，93%的氮通过该途径输入，另外7%的氮通过氮沉降、灌溉和种子3种途径输入到土壤，牧草携出是氮输出的一种主要途径，占总输出量的93%。黑麦草草地磷的总输入（输出）量为44 Gg，磷肥输入到黑麦草草地中的磷占总输入量的99%，牧草携出是最重要的磷输出项，占总输出量的77%。2014年单位面积黑麦草草地的氮输入量为427 kg N·hm-2。牧草携出的氮为396 kg N·hm-2，通过氨挥发和氧化亚氮损失的量分别为99和4.0 kg N·hm-2。黑麦草生产中的氮利用率为93%，草地土壤氮呈亏缺态，亏损量为73 kg N·hm-2。单位面积黑麦草草地磷的输入（输出）量为79 kg P·hm-2。黑麦草生产的磷利用率为77%。通过径流、侵蚀和淋洗3种途径损失的磷分别为1.2、2.5和0.065 kg P·hm-2，黑麦草土壤磷累积量为15 kg P·hm-2（表5）。

表1 2014年中国不同区域苜蓿、黑麦草和燕麦草的施肥量

>a表示专家咨询获得的数据，其他的部分源于已发表的文章（附录1）

arepresents data acquired by expert consultation, other data derived from published data (appendix 1)

表2 2014年中国不同区域栽培草地氮磷沉降、灌溉及生物固氮量

表中的沉降量指湿沉降，苜蓿的生物固氮量占总吸氮量的80%，数据源于已发表的文章

The deposition rate refers to wet deposition. Alfalfa acquires around 80% of their nitrogen from N fixation. All data used derived from published articles

2014年全国燕麦草地氮的总输入（输出）量为67 Gg，84%的氮通过氮肥输入到燕麦草草地，其余16%的氮通过氮沉降、灌溉和种子3种途径输入，牧草携出是氮输出的一种主要途径，占总输出量的69%。燕麦草地磷的总输入和总输出量为16 Gg，几乎所有的磷都是以磷肥的形式输入到燕麦草地，占总输入量的98%，牧草携出是最重要的磷输出项。2014年单位面积燕麦草草地氮输入（输出）量为217 kg N·hm-2，单位面积燕麦草每年的氮带走量为149 kg N·hm-2，通过氨挥发、氧化亚氮排放、径流、侵蚀、淋洗和反硝化6种途径损失到环境中的量分别为46、1.8、1.9、0.069、3.4和2.7 kg N·hm-2。燕麦草草地土壤氮累积量为12 kg N·hm-2。燕麦草生产中的氮利用率为69%。磷账户方面，单位面积燕麦草地磷的输入（输出）量为53 kg P·hm-2。单位面积燕麦草磷带走量为18 kg P·hm-2，通过径流、侵蚀和淋洗3种途径损失的磷分别为0.81、1.7和0.14 kg P·hm-2，土壤磷累积量为33 kg P·hm-2。燕麦草生产中的磷利用率为34%（表5）。

表3 2014年中国苜蓿、黑麦草和燕麦草的种子带入养分量和牧草养分含量

种子带入土壤中的养分含量=播种量×种子养分含量，其中播种量源于已发表的文章[28]，种子养分含量数据源于EPIC模型；牧草养分含量通过文献总结获得（附录2）

nutrient input by seed= Seeding rate × nutrient content. Seeding rate used derived from published articles, nutrient content of the seed derived from EPIC model; Forage nutrient content derived from published articles (appendix 2)

表4 2014年中国不同区域苜蓿、黑麦草和燕麦草的养分输出量

牧草养分携出=牧草产量×牧草养分含量，牧草产量源于中国草业统计2014[4]

Nutrient output by forage=forage yield × forage nutrient content, forage content derived from China Grass Industry Statistics 2014[4]

表5 2014年中国苜蓿、黑麦草和燕麦草氮磷流动特征

研究区苜蓿、黑麦草和燕麦草的面积分别为4.7×106、0.55×106和0.31×106hm2

the total research area for alfalfa, ryegrass and oat grass is 4.7×106, 0.55×106and 0.31×106hm2, respectively

2.2 2014年中国各地区主要栽培草地氮磷平衡账户

2.2.1 苜蓿氮磷平衡账户 各地区氮投入量的变化范围为194—668 kg N·hm-2，生物固氮是苜蓿氮输入的一个主要方式。氮输入量最多的5个地区分别为贵州、青海、四川、云南和广西，输入量分别为668、563、540、529和519 kg N·hm-2。山东苜蓿草地的氮输入量最小，为194 kg N·hm-2。总体而言，西部和西南地区的总投入量大于华北和东北地区。牧草的输出量方面，基本上氮输入量越大的区域，苜蓿带走的养分越多，各地区牧草携出的量分布在86—568 kg N·hm-2，江苏和西藏分别是牧草携出最多和最少的区域（图1-a）。

磷输入量和氮有相似的规律，各地区差异很大，西南和东北地区输入量较高，磷输入量最高的省份为贵州，其输入量为122 kg P·hm-2；华北和东北地区输入量低，输入量最低的省份为吉林，其磷输入量不足20 kg P·hm-2。磷输入量较多的地区一般牧草携出量也偏高，牧草携出的磷分布在8—50 kg·hm-2（图1-b）。

a：苜蓿氮账户the N budget of alfalfa；b：苜蓿磷账户the P budget of alfalfa

2.2.2 黑麦草氮磷平衡账户 2014年各地区单位面积黑麦草氮投入量差异较大，变化范围为149—623 kg N·hm-2，这与苜蓿草地的氮输入量规律相似，但是黑麦草和苜蓿氮输入项的组成有很大的区别。黑麦草作为禾本科牧草，其氮输入中不包含生物固氮，氮肥是黑麦草氮输入的主要来源。黑麦草草地氮输入最高和最低的地区分别为四川和重庆。牧草氮输出量方面，各地区牧草携出氮的量分布在165—654 kg N·hm-2，湖南和山东分别是牧草氮带走最多和最少的地区。与苜蓿不同，部分氮输入量不多的区域，牧草携出的养分却很多。以湖北为例，湖北的氮输入量在所有的省份中排在第11位，但是该省的牧草产量却位居全国第3位（低于湖南和贵州）（图2-a）。

2014年各地区单位面积黑麦草草地磷输入量和氮有相似的规律，各区域差异很大。磷输入量最高的为贵州，其输入量为146 kg P·hm-2；输入量最低的为重庆，其输入量仅为12 kg P·hm-2。牧草携出的磷分布在25—100 kg P·hm-2（图2-b）。

2.2.3 燕麦草氮磷平衡账户 2014年各地区单位面积燕麦草氮投入量差异较大，变化范围为76—334 kg N·hm-2，氮肥是燕麦草氮输入的主要来源，这与黑麦草相同。燕麦草地氮输入最高和最低的地区分别为青海和宁夏。牧草氮的输出量方面，燕麦草的牧草氮携出量远低于黑麦草，各地区燕麦草氮携出量分布在71—220 kg N·hm-2，云南和新疆分别是牧草携出氮最多和最少的区域（图3-a）。

2014年各地区单位面积燕麦草地磷输入量和氮有相似的规律，青藏高原区输入量较高。磷输入量最高的地区为西藏，其输入量为90 kg P·hm-2；输入量最低的地区为宁夏，其输入量为14 kg P·hm-2。牧草磷携出量分布在8—26 kg P·hm-2（图3-b）。

2.3 2014年中国各地区主要栽培草地氮磷利用率

2014年中国各地区苜蓿、黑麦草和燕麦草的氮磷利用率如表6。苜蓿氮利用率分布在22%—110%，江苏的氮利用率最高，达到了110%，安徽、四川、湖北、重庆和天津的氮利用率也超过了90%。内蒙古、吉林和西藏的氮利用率不足50%，其中西藏的氮利用率仅有22%，主要原因是这些地区的牧草产量较低。苜蓿的磷利用率低于氮，且不同地区间差异很大，利用率最高的为江苏（128%），四川的磷利用率也超过了100%。磷利用率最低的地区是西藏，仅为8.0%。新疆、青海等西部地区的磷利用率不足30%。

a：黑麦草氮账户the N budget of ryegrass；b：黑麦草磷账户The P budget of ryegrass

a：燕麦草氮平衡账户the N budget of oat grass；b：燕麦草磷平衡账户the P budget of oat grass

2014年各地区黑麦草氮利用率远高于苜蓿，分布在60%—227%，重庆和四川分别是利用率最高和最低的地区。江西的黑麦草氮利用率超过了200%，江苏、安徽、湖北、河南、湖南、贵州、广西和西藏的氮利用率在100%—200%。黑麦草磷利用率同样高于苜蓿，但区域间差异很大，分布在55%—410%，其中重庆的利用率达到了410%，江苏和湖南的磷利用率也较高（＞200%），湖北、山东、广西和福建的利用率分布在100%—200%。黑麦草磷利用率最低的为贵州（55%）。

2014年各地区燕麦草的氮利用率在苜蓿和黑麦草之间，分布在22%—180%，四川、湖北、重庆、云南和内蒙古的氮利用率高于100%。西藏和新疆两地的利用率较低，分别为22%和23%。燕麦草的磷利用率低于氮，分布在10%—119%。四川和内蒙古两地的利用率高于100%，西藏、新疆、陕西、山西、青海等西部地区的磷利用率不足30%。

表6 2014年中国各地区苜蓿、黑麦草和燕麦草氮磷利用率

2.4 2014年中国各地区主要栽培草地氮磷环境损失

2.4.1 苜蓿氮磷环境损失 西南地区和青藏高原区氮的环境损失量较高，其中，西藏的氮环境损失量最高（126 kg N·hm-2）。华北和长江中下游地区的损失量较少，以江苏为例，该省氮的环境损失量仅为7.4 kg N·hm-2。在各损失项中，氨挥发是最主要部分，而通过侵蚀损失的氮很少，基本上不超过1.0 kg N·hm-2（图4-a）。磷排放的区域分布规律与氮相似，同样是西南地区高于其他地区。贵州的磷环境排放量远高于其他区域（18 kg P·hm-2），主要原因有两方面，一是当地种植苜蓿时，磷肥的投入量较大，二是当地的喀斯特地貌容易发生侵蚀。各项磷的环境损失中，侵蚀是最主要的部分，其次是径流，通过淋洗损失的养分最少（图4-b）。

2.4.2 黑麦草氮磷环境损失 四川省氮环境损失量最高，为198 kg N·hm-2，重庆市的氮环境损失量最少，为31 kg N·hm-2。在各损失项中，氨挥发是最主要部分（图5-a）。各地区磷的损失量变化范围为0.72—24 kg P·hm-2，西南地区黑麦草磷的环境损失远高于其他地区，其原因与贵州地区苜蓿磷的环境损失项相同。在磷的各损失项中，侵蚀是最主要的部分（图5-b）。

2.4.3 燕麦草氮磷环境损失 西部地区的氮环境损失量明显高于其他区域，西藏、新疆和青海3省（自治区）的氮环境排放量超过100 kg N·hm-2，相比之下，四川和重庆等南部省（直辖市）的氮损失量较少，不足30 kg N·hm-2。此外，山西的氮环境排放量明显高于周边的其他省份，主要原因是山西的氮肥施用量明显高于其他省份（图6-a）。单位面积燕麦草磷损失规律与氮损失相似，均为西部地区高于其他区域。云南的磷环境排放量最高，超为5.0 kg P·hm-2。侵蚀是最主要的磷损失项（图6-b）。

2.5 2014年中国地区主要栽培草地单位面积土壤氮磷累积量

2.5.1 苜蓿草地土壤氮磷累积量 2014年，除江苏、安徽和四川外，其余地区苜蓿氮土壤累积量均为正值。西藏、青海、广西、云南、贵州和辽宁的土壤累积量超过90 kg N·hm-2，陕西、内蒙古、新疆、云南、宁夏、山西、甘肃、河北、黑龙江、辽宁、北京、湖南、河南、山东、江西、天津、重庆和湖北的土壤累积量分布在0—90 kg N·hm-2（图7）。

图4 2014年中国不同地区苜蓿氮磷环境排放

图5 2014年中国不同地区黑麦草氮磷环境排放

全国26个有苜蓿分布地区中，24个地区的磷累积量为正值。西北地区苜蓿草地土壤磷累积量超过40 kg P·hm-2，此外，广西和贵州两地的磷累积量同样超过40 kg P·hm-2。西北地区苜蓿草地土壤磷累积量过高的原因是通过施肥输入土壤中的磷较多，而西藏、陕西、广西和山西的磷累积量高于贵州是因为贵州通过收获牧草输出的磷较多。河北、黑龙江、山东、湖南、云南、江西、北京、安徽、河南、吉林、重庆、天津、辽宁和湖北等地的苜蓿土壤磷累积量分布在0—40 kg P·hm-2。江苏和四川的苜蓿土壤磷累积量为负值，亏缺量分别为14和15 kg P·hm-2，主要原因是这两个地区苜蓿肥料投入量不高，但收获苜蓿携出的磷较多（图7）。

2.5.2 黑麦草草地土壤氮磷累积量 全国共有16个地区有黑麦草分布，仅四川、山东和福建3省的氮累积量为正值。陕西、云南、广东、安徽、河南、广西、西藏、江苏、贵州、湖北、江西、重庆和湖南的黑麦草土壤氮累积量为负值，主要原因是通过收获黑麦草携出的氮较多，其中湖南、重庆和江西的土壤亏缺量超过200 kg N·hm-2（图8）。

贵州、西藏、广东、四川、陕西、云南、河南和江西的黑麦草土壤磷累积量为正值。湖南、重庆、江苏、福建、广西、湖北、山东和安徽的黑麦草土壤磷累积量为负值。其中，湖南、重庆和江苏的亏缺量超过30 kg P·hm-2，主要原因是收获牧草携出了大量的磷（图8）。

2.5.3 燕麦草地土壤氮磷累积量 燕麦草在西部和南部的12个地区有分布，其中西藏和新疆两地燕麦草地土壤氮累积量超过70 kg N·hm-2，主要原因是当地通过施肥带入土壤中的氮很高，远大于其他地区，山西、青海和陕西燕麦草地土壤的氮同样呈盈余态，累积量分布在0—70 kg N·hm-2。宁夏、甘肃、内蒙古、湖北、重庆、四川和云南的燕麦草土壤氮呈亏缺态（图9）。

图6 2014年中国不同地区燕麦草氮磷环境排放

图7 2014年中国不同地区苜蓿草地土壤氮磷累积量

图8 2014年中国不同地区黑麦草草地土壤氮磷累积量

同土壤氮累积相似，西藏、新疆、山西和青海的燕麦草地磷有较高的累积，累积量＞40 kg P·hm-2，造成土壤过多磷累积的原因为磷肥输入过多。陕西、甘肃、重庆、宁夏、湖北和云南的土壤磷也呈盈余态，累积量在0—40 kg P·hm-2。内蒙古和四川燕麦草地土壤磷略有亏缺（图9）。

图9 2014年中国不同地区燕麦草地土壤氮磷累积量

3 讨论

3.1 中国栽培草地氮磷流动特征

本研究发现，除黑麦草草地土壤氮外，中国主要栽培草地土壤氮磷养分均有不同程度的累积。而黑麦草土壤氮亏缺的主要原因是黑麦草产量较高，通过收获牧草从土壤中带走了大量的养分。本文中牧草产量数据为2014年数据，当年中国黑麦草平均产量为16 t·hm-2，湖南的黑麦草产量更是达到了27 t·hm-2（表6），这就导致在计算中2014年黑麦草携出了大量的养分，使土壤的氮磷呈亏缺状态。总结已发表的文献，发现通过肥效和播种实验获得的黑麦草干草的平均产量为10 t·hm-2，远低于2014年中国草业统计中的产量（16 t·hm-2），按此产量标准来计算，中国的黑麦草草地土壤氮则与2.5.2中计算的结果相反，多数地区呈养分盈余态。

此外，本研究还发现不同区域氮磷流动特征差异很大。关于主要栽培牧草生长规律和营养需求规律的研究相对滞后，营养调控技术主要参考大田作物的施肥技术及国外的营养调控技术，缺乏合理的技术规范[9]。加之中国以粮为纲的思想，人们对于牧草科学施肥的重要性认识不足，农民认为牧草也是草，不需要管理，生产上采用盲目性的经验施肥比较普遍[9,29]。国家牧草产业体系的人员于2012年进行的调研进一步证实了上述的观点，他们对全国主要苜蓿产区进行了调研，其中有18个样点不施氮肥，而另外7个样点的氮肥施用量超过150 kg N·hm-2，不同地区磷肥的使用比例变化也很大，此外，施肥的样点大多根据经验沿用大田作物的管理技术[9]。

由表6可知，部分地区的氮磷利用率超过100%，这意味着当地栽培草地的土壤氮磷处于亏缺态。但是当考虑整个轮作系统时，土壤养分亏缺量或许比较小。以黑麦草为例，黑麦草主要分布在中国南方地区，与水稻轮作[18]。过高的施肥量会使水稻季土壤养分出现相当的累积，当施氮量超过190 kg·hm-2时，稻田表观盈余量超过130 kg·hm-2[30]。因而，部分地区栽培草地养分亏缺是否会造成生产的不可持续有待进一步分析。本文中黑麦草施肥量超过300 kg N·hm-2，但黑麦草氮利用率依旧很高，原因可能是本研究未考虑前茬作物（水稻）生育期的养分盈余量。

3.2 中国栽培草地养分管理优化策略

基于本研究结果，对于中国主要栽培牧草养分管理的优化策略包括以下几个方面：

首先，中国苜蓿和燕麦草生产中可以适当的减少氮磷肥料的使用，但是黑麦草生产中应补施氮肥。全国苜蓿和燕麦草的土壤氮累积量分别为50和12 kg N·hm-2，磷累积量分别为46和33 kg P·hm-2。两种草地土壤的氮磷略有盈余，但相比大田作物，盈余量不大。以中国华北平原小麦-玉米轮作体系为例，该种植模式下，土壤每年的氮和磷累积量分别为227 kg N·hm-2和53 kg P·hm-2[31]。黑麦草草地土壤氮呈亏缺态，亏缺量为73 kg N·hm-2，因而应补施尿素150 kg·hm-2，施肥时期可以选择每次刈割后。黑麦草草地土壤略有盈余，无需补充磷肥，但是在实际的管理中，应考虑牧草的轮作方式。

其次，牧草养分流动的空间特征差异较大，需因地制宜的调控养分管理策略。苜蓿的主产区土壤的氮磷有部分盈余，当地应适当减少氮磷肥料的施用量，而四川、江苏等省份苜蓿草地土壤氮和磷有不同程度的亏缺，应补施肥料。全国大部分黑麦草种植区的土壤氮亏缺严重，应增施氮肥。黑麦草主产区（四川、贵州、云南和江西）的土壤磷基本平衡，无需进行养分的调控。各地区燕麦草养分管理方面，很多氮投入不足的区域，土壤磷却呈盈余态，燕麦草的肥料施用结构不平衡，在实际生产中应注意减施磷肥，或者将氮肥和家畜粪尿配合施用。

最后，氨挥发是氮肥损失的最主要途径，可选择有机肥做基肥，延长肥效，减少牧草生育后期氮肥的表施。部分降雨量高且土壤保水性较差的地区在实际生产中应注意防止磷损失。很多区域的牧草种植体系，氨挥发损失量超过100 kg N·hm-2，有很大的减排空间，因而在氮肥的施用方式上，氨态氮肥要避免表施。侵蚀是磷损失的最主要途径，由于磷在土壤中的移动性较弱，因此总体的损失量不太高，但是部分降雨量高且土壤保水性较差的地区仍有较大的损失量。

3.3 模型的不确定性

本研究利用NUFER模型，模拟了中国苜蓿、黑麦草和燕麦草地氮磷流动特征，但是由于模型的算法和参数获取的缺陷，也存在一些不确定性。

牧草施肥量方面，由于农民对牧草施肥管理的盲目性，获取不同区域农户施肥数据的难度较大，本章所采用的施肥量数据，是通过专家咨询和文献搜集两种方式协同获取的。其中青藏高原区的情况特殊，由于历史原因，重牧轻草的观念在当地比较普遍[9]。当地栽培草地的肥料来源多为畜禽粪尿，草地的管理方式多为底肥施用牛羊粪38 t·hm-2，配施一定量的二铵，约150 kg·hm-2，后期追施少量尿素，约75 kg·hm-2（专家咨询），本研究中青海、西藏苜蓿、黑麦草和燕麦草的施肥量均按上述的标准计算。

模型参数获取方面，不同种类栽培草地养分排放因子由综述文献获得。但是关于燕麦和黑麦草环境排放参数的报道很少，考虑到上述两种作物与小麦等大田作物特性相似，且田间管理也与大田作物相差不大，故本研究中燕麦和黑麦草的排放参数沿用大田作物的数据。

4 结论

（1）2014年，全国栽培苜蓿、黑麦草和燕麦草生产体系的氮输入（输出）量差异较大。生物固氮是苜蓿获取氮素的重要途径，田间管理需减少氮肥投入，燕麦和黑麦草主要通过化肥获取所需养分，应注重肥料施用；（2）2014年，3种栽培牧草氮的利用率均高于60%，磷的利用率相对较低；（3）2014年，不同区域苜蓿、黑麦草和燕麦草的氮磷利用率和环境排放量差异较大；（4）2014年，黑麦草生产氮呈亏缺态，需增施氮肥，其余栽培草地氮磷均有不同程度的盈余，应控制肥料过多施用。

[1] 李凌浩, 路鹏, 顾雪莹, 毛小涛, 高树琴, 张英俊. 人工草地建设原理与生产范式. 科学通报, 2016, 61(2): 193-200.

LI L H, LU P, GU X Y, MAO X T, GAO S Q, ZHANG Y J. Principles and paradigms for developing artificial pastures., 2016, 61(2): 193-200. (in Chinese)

[2] 方精云, 白永飞, 李凌浩, 蒋高明, 黄建辉, 黄振英, 张文浩, 高树琴. 我国草原牧区可持续发展的科学基础与实践. 科学通报, 2016, 61(2): 155-164.

FANG J Y, BAI Y F, LI L H, JIANG G M, HUANG J H, HUANG Z Y, ZHANG W H, GAO S Q. Scientific basis and practical ways for sustainable development of China’s pasture regions., 2016, 61(2): 155-164. (in Chinese)

[3] 沈海花, 朱言坤, 赵霞, 耿晓庆, 高树琴, 方精云. 中国草地资源的现状分析. 科学通报, 2016, 61(2): 139-154.

SHEN H H, ZHU Y K, ZHAO X, GENG X Q, GAO S Q, FANG J Y. Analysis of current grassland resources in China., 2016, 61(2): 139-154. (in Chinese)

[4] 全国畜牧总站. 中国草业统计. http://www.digitalgrss.cn, 2015.

National Animal Husbandry Station. China Grass Industry Statistics. http://www.digitalgrss.cn, 2015. (in Chinese)

[5] 中共中央国务院. 中央一号文件. http://www.farmer.com.cn/uzt/ ywj/gea/201601/t20160128_1176612_2.htm, 2016.

The State Council of China. The central first document. http://www. farmer.com.cn/uzt/ywj/gea/201601/t20160128_1176612_2.htm, 2016. (in Chinese)

[6] 洪绂曾, 王元素. 中国南方人工草地畜牧业回顾与思考. 中国草地学报, 2006, 28(2): 71-75, 78.

HONG F Z, WANG Y S. The strategical thinking of grassland agriculture in Southern China., 2006, 28(2): 71-75, 78. (in Chinese)

[7] 张新时, 唐海萍, 董孝斌, 李波, 黄永梅, 龚吉蕊. 中国草原的困境及其转型. 科学通报, 2016, 61(2): 165-177.

ZHANG X S, TANG H P, DONG X B, LI B, HUANG Y M, GONG J R. The dilemma of steppe and it’s transformation in China., 2016, 61(2): 165-177. (in Chinese)

[8] O’MARA F P. The role of grasslands in food security and climate change., 2012, 110(6): 1263-1270.

[9] 国家牧草产业技术体系. 中国栽培草地. 北京: 科学出版社, 2015.

National Technical System of Grass Industry.. Beijing: Science Press, 2015. (in Chinese)

[10] RIVERO D S, VILLAMIL A M S, MAHECHA O M, TERAN A M, RIVERO M S, NAVARRO A M S. Evaluation of the effect of two types of fertilizer on the growth, development and productivity of hydroponic green forage oat (L.) and ryegrass (Lam.) as a biomass source., 2016, 50: 385-390.

[11] 李小坤, 鲁剑巍, 刘晓伟, 刘威, 潘福霞, 魏云霞. 配方施肥对一年生黑麦草产草量及品质的影响. 草业科学, 2011, 28(9): 1666-1670.

LI X K, LU J W, LIU X W, LIU W, PAN F X, WEI Y X. Effect of fertilization on the yield and forage quality of annual ryegrass., 2011, 28(9): 1666-1670. (in Chinese)

[12] 张菁, 张于卉, 马力, 周鹏, 安渊. 播期和施肥对冬闲田紫花苜蓿生长和品质的影响. 中国草地学报, 2015, 37(6): 35-41.

ZHANG J, ZHANG Y H, MA L, ZHOU P, AN Y. Effect of sowing data and applying fertilizer on growth and quality of alfalfa in winter fallow land., 2015, 37(6): 35-41. (in Chinese)

[13] 陶雪, 苏德荣, 寇丹, 乔阳. 西北旱区灌溉方式对苜蓿生长及水分利用率的影响. 草地学报, 2016, 24(1): 114-120.

TAO X, SU D R, KOU D, QIAO Y. Effect of irrigation methods on growth and water use efficiency of alfalfa in arid Northwest China., 2016, 24(1): 114-120. (in Chinese)

[14] 陈萍, 沈振荣, 迟海峰, 杨东海, 谷传申, 杜婉君, 鲍举文. 不同施肥处理对紫花苜蓿产量和株高的影响. 作物杂志, 2013(1): 91-94.

CHEN P, SHEN Z R, CHI H F, YANG D H, GU C S, DU W J, BAO J W. Effects of different fertilization on output and plant height of alfalfa., 2013(1): 91-94. (in Chinese)

[15] 肖相芬, 周川姣, 周顺利, 胡跃高, 任长忠, 郭来春, 王春龙. 燕麦氮吸收利用特性与适宜施氮量的定位研究. 中国农业科学, 2011, 44(22): 4618-4626.

XIAO X F, ZHOU C J, ZHOU S L, HU Y G, REN C Z, GUO L C, WANG C L. Characteristics of N absorption, utilization and optimum N-fertilizer rate based on a long-term fertilization experiment of oats., 2011, 44(22): 4618-4626. (in Chinese)

[16] LI W X, LU J W, CHEN F, LI X K. Effect of NPK application on yield, nutrients and water utilization under sudangrass and ryegrass rotation regime., 2010, 9(7): 1026-1034.

[17] 顾梦鹤, 王涛, 杜国祯. 施肥对高寒地区多年生人工草地生产力及稳定性的影响. 兰州大学学报(自然科学版), 2010, 46(6): 59-63.

GU M H, WANG T, DU G Z. Effects of fertilization on productivity and stability of cultivated grassland in an alpine region., 2010, 46(6): 59-63. (in Chinese)

[18] 杨曾平, 聂军, 廖育林, 周兴, 谢坚, 鲁艳红, 纪雄辉, 吴家梅, 谢运河. 不同施肥量对稻田一年生黑麦草产量及氮磷钾吸收的影响. 中国农学通报, 2015, 31(30): 173-180.

YANG Z P, NIE J, LIAO Y L, ZHOU X, XIE J, LU Y H, JI X H, WU J M, XIE Y H. Effect of different fertilizations on biomass yield and N, P, K absorption of annual ryegrass in rice field., 2015, 31(30): 173-180. (in Chinese)

[19] GU B J, JU X T, CHANG S X, GE Y, CHANG J. Nitrogen use efficiencies in Chinese agricultural systems and implications for food security and environmental protection., 2017, 17(4): 1217-1227.

[20] ZHOU Q K, DARYANTO S, XIN Z M, LIIU Z M, LIU M H, CUI X, WANG L X. Soil phosphorus budget in global grasslands and implications for management., 2017, 144: 224-235.

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

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

[22] 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.

[23] LIU X J, ZHANG Y, HAN W X, TANG A H, SHEN J L, CUI Z L, VITOUSEK P, ERISMAN J W, GOULDING K, CHRISTIE P, FANGMEIER A, ZHANG F S. Enhanced nitrogen deposition over China., 2013, 494(7438): 459-462.

[24] ZHU J, WANG Q, HE N, SMITH M D, ELSER J J, DU J, YUAN G, YU G, YU Q. Imbalanced atmospheric nitrogen and phosphorus depositions in China: Implications for nutrient limitation., 2016, 121(6): 1605-1616.

[25] 陈敏鹏, 陈吉宁. 中国种养系统的氮流动及其环境影响. 环境科学, 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)

[26] 王激清, 马文奇, 江荣风, 张福锁. 中国农田生态系统氮素平衡模型的建立及其应用. 农业工程学报, 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)

[27] CARLSSON G, HUSS-DANELL K. Nitrogen fixation in perennial forage legumes in the field., 2003, 253(2): 353-372.

[28] 国家牧草产业技术体系. 牧草标准化生产管理技术规范. 北京: 科学出版社, 2014.

National Technical System of Grass Industry.Beijing: Science Press, 2014. (in Chinese)

[29] 任继周, 胥刚, 李向林, 林慧龙, 唐增. 中国草业科学的发展轨迹与展望. 科学通报, 2016, 61(2): 178-192.

REN J Z, XU G, LI X L, LIN H L, TANG Z. Trajectory and prospect of China’s prataculture., 2016, 61(2): 178-192. (in Chinese)

[30] 罗付香, 林超文, 朱永群, 姚莉, 张建华, 王谦. 氮肥用量对四川紫色丘陵区水稻氮素吸收及土壤氮素表观平衡的影响. 西南农业学报, 2016, 29(8): 1093-1097.

LUO F X, LIN C W, ZHU Y Q, YAO L, ZHANG J H, WANG Q. Effect of nitrogen fertilizer amount on rice nitrogen uptake and soil nitrogen apparent balance of Sichuan purple hilly region., 2016, 29(8): 1093-1097. (in Chinese)

[31] VITOUSEK P M, NAYLOR R, CREWS T, DAVID M B, DRINKWATER L E, HOLLAND E, JOHNES P J, KATZENBERGER J, MARTINELLI L A, MATSON P A, NZIGUHEBA G, OJIMA D, PALM C A, ROBERTSON G P, SANCHEZ P A, TOWNSEND A R, ZHANG F S. Nutrient imbalances in agricultural development., 2009, 324(5934): 1519-1520.

（责任编辑 岳梅）

附录一：栽培草地肥料施用量数据来源

[S1] 侯湃, 刘自学, 刘艺杉, 苏爱莲, 王跃栋. 北京平原区紫花苜蓿施肥组合试验. 草业科学, 2014, 31(1): 144-149.

[S2] 李源, 王赞, 高洪文, 刘贵波. 磷钾配合施肥对紫花苜蓿产量的影响. 华北农学报, 2008, 23(S2): 315-318.

[S3] 张丽娟. 根瘤菌、施肥等技术对紫花苜蓿生产性能的影响[D]. 太谷: 山西农业大学, 2005.

[S4] 韩思训, 王森, 高志岭, 王殿武. 不同施肥条件下苜蓿产量、氮素累积量及肥料氮素利用率研究. 华北农学报, 2014, 29(6): 220-225.

[S5] 王秀领, 徐玉鹏, 闫旭东, 赵忠祥. 磷钾肥对紫花苜蓿鲜草产量的影响. 河北农业科学, 2007(2): 64-66.

[S6] 谢勇. 坝上地区紫花苜蓿施肥研究及其土壤养分丰缺指标体系的初步建立[D]. 雅安: 四川农业大学, 2010.

[S7] 谢勇, 孙洪仁, 张新全, 刘爱红, 孙雅源, 刘治波, 丁宁. 坝上地区紫花苜蓿氮磷、钾肥料效应与推荐施肥量. 中国草地学报, 2012, 34(2): 52-57.

[S8] 万井尉. 紫花苜蓿人工草地施肥试验. 现代畜牧兽医, 2015(1): 10-14.

[S9] 王喜春. 牧草不同施肥量及其经济效益的观察. 吉林畜牧兽医, 1981(2): 5-6.

[S10] 柴凤久, 许金玲, 李红, 刘崇生. 紫花苜蓿施肥试验研究. 中国草地, 2004, 26(2): 81-82.

[S11] 柴凤久, 高海娟, 刘泽东, 尤海洋. 施肥和灌水对紫花苜蓿草产量的影响. 草业与畜牧, 2012(7): 21-23.

[S12] 曲善民, 孙旭东, 杜广明, 杜吉到, 郑殿峰, 冯乃杰. 不同时期定量施肥对苜蓿生产性能的影响. 当代畜牧, 2009(10): 50-52.

[S13] 张东铁, 唐凤兰, 朱瑞芬, 刘凤歧, 刘杰淋, 陈积山, 韩微波. 盐碱地苜蓿播量、施肥期和施肥量对其产量的影响. 黑龙江农业科学, 2011(12): 116-119.

[S14] 储国良, 丁剑英, 付反生, 王敬根, 王全洪, 刘亚柏. 南方丘陵坡地紫花苜蓿施肥效应初探. 中国乳业, 2002(5): 17-19.

[S15] 姜慧新, 沈益新, 翟桂玉, 刘展生. 施磷对紫花苜蓿分枝生长及产草量的影响. 草地学报, 2009(5): 588-592.

[S16] 李圣昌, 朱树森. 紫花苜蓿施用磷肥试验研究. 草业科学, 1993, 10(1): 41-43.

[S17] 温洋, 金继运, 黄绍文, 左余宝, 李志杰. 不同磷水平对紫花苜蓿产量和品质的影响. 土壤肥料, 2005(2): 21-24.

[S18] SHEN Y X, JIANG H X, ZHAI G Y, CAI Q S. Effects of cutting height on shoot regrowth and forage yield of alfalfa (L.) in a short-term cultivation system., 2013, 59(2): 73-79.

[S19] 朱华敏, 李海梅, 杨国锋, 孙娟, 刘志英. 施肥对紫花苜蓿生产性能的影响. 北方园艺, 2013(14): 174-177.

[S20] 胡华锋, 肖金帅, 郭孝, 介晓磊, 刘世亮, 化党领, 鲁剑巍, 刘芳. 氮磷钾肥配施对黄河滩区紫花苜蓿产量和品质的影响. 湖南农业大学学报(自然科学版), 2009(2): 178-180, 188.

[S21] 刘继先, 李伟. 平衡施肥对紫花苜蓿产量和品质的影响. 南方农业, 2012, 6(5): 6-8.

[S22] 张健. 三峡库区牧草种植区划及适生牧草栽培利用技术研究[D]. 重庆: 西南大学, 2007.

[S23] 张健, 沈贵平. 不同施肥处理对紫花苜蓿鲜草产量和营养品质的影响研究. 草业与畜牧, 2008(9): 1-4.

[S24] 邓蓉, 张定红, 陈武, 向清华. 施肥对黔中地区混播草地牧草生长性能的影响. 畜牧与兽医, 2004, 36(3): 12-15.

[S25] 黄梅. 不同施肥处理对紫花苜蓿产量影响及效益分析. 当代畜牧, 2015(11): 92-93.

[S26] LI W C, LU J W, CHEN F, LI X K. Effect of NPK application on yield, nutrients and water utilization under sudangrass and ryegrass rotation regime., 2010, 9(7): 1026-1034.

[S27] 杨曾平, 聂军, 廖育林, 周兴, 谢坚, 鲁艳红, 纪雄辉, 吴家梅, 谢运河. 不同施肥量对稻田一年生黑麦草产量及氮磷钾吸收的影响. 中国农学通报, 2015, 31(30): 173-180.

[S28] 曹景勤. 第四纪红色粘土区氮磷肥交互作用对黑麦草产量和养分累积的影响. 土壤通报, 1995, 26(3): 120-121.

[S29] 刘经荣, 张美良, 郑群英, 石庆华, 王少先. 不同施N水平对黑麦草产量和品质的效应. 江西农业大学学报, 2003, 25(6): 844-847.

[S30] 方华舟. 施肥对黑麦草产量和质量的影响. 中国农村小康科技, 2007(8): 41, 79-80.

[S31] 李文西, 鲁剑巍, 鲁君明, 陈防, 李小坤. 苏丹草-黑麦草轮作中氮磷钾配施对饲草产量与养分积累的影响. 草原与草坪, 2009(6): 5-10.

[S32] 李文西. 苏丹草-黑麦草轮作制中连续施肥对饲草生长与土壤肥力的影响[D]. 武汉: 华中农业大学, 2010.

[S33] 李小坤, 李文西, 鲁剑巍, 黄元仿, 刘光文. 施肥对黑麦草产量和氮磷钾养分吸收的影响. 中国农学通报, 2006, 22(7): 331-334.

[S34] 李小坤, 鲁剑巍, 鲁君明, 陈防, 黄元仿. 磷肥用量对黑麦草产量及经济效益的影响. 草业科学, 2006, 23(10): 18-22.

[S35] 李小坤, 鲁剑巍, 刘晓伟, 刘威, 潘福霞, 魏云霞. 配方施肥对一年生黑麦草产草量及品质的影响. 草业科学, 2011, 28(9): 1666-1670.

[S36] 刘晓伟, 杨娟, 李文西, 李小坤, 鲁剑巍. 氮磷钾肥施用对多花黑麦草越冬期生理指标和产草量的影响. 草地学报, 2010, 18(4): 584-588.

[S37] 杨娟. 施肥对两种轮作牧草的产量和品质效果及其机理研究[D]. 武汉: 华中农业大学, 2009.

[S38] 鲁剑巍, 李小坤, 梁友光, 鲁君明, 陈防. 平衡施肥对黑麦草生长及产量的影响. 水利渔业, 2004, 24(2): 20-22.

[S39] 占丽平, 鲁剑巍, 杨娟, 李文西, 李小坤, 刘晓伟. 施肥对黑麦草生长和产量的影响. 草业科学, 2011, 28(2): 260-265.

[S40] 罗婷, 李林. 多花黑麦草不同品种草产量研究. 中国农学通报, 2010, 26(8): 259-263.

[S41] 申晓萍, 李仕坚, 朱梅芳, 樊露晓. 刈割次数对南方冬闲田种植黑麦草产量及品质的影响. 安徽农业科学, 2009, 37(35): 17445-17446, 17461.

[S42] 孙学映, 刘春英, 朱体超, 陈光蓉. 黑麦草高产栽培技术研究. 南方农业学报, 2013, 44(7): 1105-1109.

[S43] 陈国南, 尚以顺. 几个栽培因子对多花黑麦草贵草2号鲜草产量的影响. 种子, 2009, 28(5): 84-85.

[S44] 陈培燕, 左相兵, 尚以顺, 赵熙贵, 陈国南. 特高黑麦草不同用种量及施肥处理的产量效果比较. 贵州畜牧兽医, 2006, 30(2): 6-7.

[S45] 陈志祯, 黄河斋, 周泽英, 何胜礼, 薛万珍, 陈克庆, 罗国琼. 不同施肥水平对稻田免耕黑麦草生长及产量的影响. 贵州农业科学, 2011, 39(4): 68-70.

[S46] 何静, 李辰琼, 陈瑞祥. 贵草1号多花黑麦草密度和施肥研究. 草业与畜牧, 2010(4): 9-10.

[S47] 皇甫江云, 刘贵林, 赵楠, 赵熙贵. 杰威多花黑麦草及其在贵州喀斯特山区的栽培管理及利用技术. 四川草原, 2005(11): 61-62.

[S48] 尚以顺, 舒健虹, 陈燕萍, 熊先勤, 赵熙贵. 不同施肥水平及播种量对一年生黑麦草A2003产草量的影响. 贵州农业科学, 2008, 36(5): 135-136, 146.

[S49] 田松军, 何胜江, 雷荷仙, 付强, 何俊, 赵相勇. 不同施肥水平对冬闲田牧草产量的影响. 贵州畜牧兽医, 2006, 30(4): 3-5.

[S50] 王柏渊. 施肥和刈割对黔南区轮作多花黑麦草产草量及体外消化的研究[D]. 雅安: 四川农业大学, 2009.

[S51] 张明均, 张子琴. 多花黑麦草在贵州地区的栽培技术. 农技服务, 2008, 25(4): 97-98.

[S52] 池银花, 张爱华. 氮肥对禾本科牧草——黑麦草产量和品质的影响. 福建畜牧兽医, 2003, 25(6): 29-30.

[S53] 池银花, 张爱华. 追施氮肥对多花黑麦草产量和品质的影响. 福建畜牧兽医, 2006, 28(2): 10-11.

[S54] 辛国荣, 杨中艺, 徐亚幸, 陈三有, 郭仁东. “黑麦草-水稻”草田轮作系统的研究Ⅴ.稻田冬种黑麦草的优质高产栽培技术. 草业学报, 2000, 9(2): 17-23.

[S55] WANG S P, WANG Y F, SCHNUG E, HANEKLAUS S, FLECKENSTEIN J. Effects of nitrogen and sulphur fertilization on oats yield, quality and digestibility and nitrogen and sulphur metabolism of sheep in the Inner Mongolia Steppes of China., 2002, 62(2): 195-202.

[S56] 蒲珉锴. “红试2号”燕麦新品系品比试验及生产性能评价[D]. 雅安: 四川农业大学, 2012.

[S57] 才让吉, 王巧玲, 王贵珍, 刘丽, 花立民. 不同播量、氮磷肥互作对高寒牧区燕麦产量和品质的影响. 中国草食动物科学, 2015, 35(5): 30-33, 61.

[S58] 刘刚. 青藏高原饲用燕麦种质资源评价与筛选[D]. 兰州: 甘肃农业大学, 2006.

[S59] 蒲小鹏. 氮肥对高寒地区燕麦草地生产性能及土壤酶活性的影响[D]. 兰州: 甘肃农业大学, 2005.

[S60] WU G L, WANG M R, GAO T, HU T M, DAVISON G. Effects of mowing utilization on forage yield and quality in five oat varieties in alpine area of the eastern Qinghai-Tibetan Plateau., 2010, 9(4): 461-466.

附录二：栽培牧草养分含量数据来源

[SS1] 张丽娟. 根瘤菌、施肥等技术对紫花苜蓿生产性能的影响[D]. 太谷: 山西农业大学, 2005.

[SS2] 曲善民, 孙旭东, 杜广明, 杜吉到, 郑殿峰, 冯乃杰. 不同时期定量施肥对苜蓿生产性能的影响. 当代畜牧, 2009(10): 50-52.

[SS3] 温洋, 金继运, 黄绍文, 左余宝, 李志杰. 不同磷水平对紫花苜蓿产量和品质的影响. 土壤肥料, 2005(2): 21-24.

[SS4] 胡华锋, 肖金帅, 郭孝, 介晓磊, 刘世亮, 化党领, 鲁剑巍, 刘芳. 氮磷钾肥配施对黄河滩区紫花苜蓿产量和品质的影响. 湖南农业大学学报(自然科学版), 2009(2): 178-180, 188.

[SS5] 刘继先, 李伟. 平衡施肥对紫花苜蓿产量和品质的影响. 南方农业, 2012, 6(5): 6-8.

[SS6] 张健. 三峡库区牧草种植区划及适生牧草栽培利用技术研究[D]. 重庆: 西南大学, 2007.

[SS7] 范志东. 毛乌素沙地紫花苜蓿水肥耦合效应的分析研究[D]. 呼和浩特: 内蒙古农业大学, 2012.

[SS8] 贾珺, 韩清芳, 周芳, 贾志宽, 王俊鹏, 杨保平. 氮磷配比对旱地紫花苜蓿产量构成因子及营养成分的影响. 中国草地学报, 2009, 31(3): 77-82.

[SS9] 刘晓静, 张进霞, 李文卿, 范俊俊. 施肥及刈割对干旱地区紫花苜蓿产量和品质的影响. 中国沙漠, 2014, 34(6): 1516-1526.

[SS10] 马孝慧, 阿不来提·阿不都热依木, 孙宗玖, 赵清, 卡米力.氮磷、钾、硫肥对苜蓿产量和品质影响. 新疆农业大学学报, 2005, 28(1): 18-21.

[SS11] 汪茜. 施肥对苜蓿低产田增产效应研究[D]. 兰州: 甘肃农业大学, 2015.

[SS12] 王聪明. 高产优质紫花苜蓿施肥效应研究[D]. 呼和浩特: 内蒙古农业大学, 2011.

[SS13] 王永宁, 李跃进, 张曙光, 刘慧. 河套地区盐碱地苜蓿施肥试验初探. 北方农业学报, 2016, 44(3): 18-22.

[SS14] 杨恒山, 曹敏建, 李春龙, 李凤山. 苜蓿施用磷、钾肥效应的研究. 草业科学, 2003, 20(11): 19-22.

[SS15]张进霞. 氮磷调控对紫花苜蓿生产性能及其土壤养分特性的影响[D]. 兰州: 甘肃农业大学, 2014.

[SS16] 张学洲, 兰吉勇, 张荟荟, 张学沣, 张一弓, 王玉, 李学森. 不同施肥配比对多叶型紫花苜蓿产量、品质和效益的影响. 现代农业科技, 2016(4): 270-273.

[SS17] 赵云, 谢开云, 杨秀芳, 李向林. 氮磷钾配比施肥对敖汉苜蓿产量和品质的影响. 草业科学, 2013, 30(5): 723-727.

[SS18] 李荣霞. 不同施肥水平对紫花苜蓿产量、营养吸收及土壤肥力的影响[D]. 乌鲁木齐: 新疆农业大学, 2007.

[SS19] 蔺蕊. 北疆盐碱地苜蓿施肥参数与高产施肥初步研究[D]. 乌鲁木齐: 新疆农业大学, 2004.

[SS20] LI W X, LU J W, CHEN F, LI X K. Effect of NPK application on yield, nutrients and water utilization under sudangrass and ryegrass rotation regime., 2010, 9(7): 1026-1034.

[SS21] 杨曾平, 聂军, 廖育林, 周兴, 谢坚, 鲁艳红, 纪雄辉, 吴家梅, 谢运河. 不同施肥量对稻田一年生黑麦草产量及氮磷钾吸收的影响. 中国农学通报, 2015, 31(30): 173-180.

[SS22] 曹景勤. 第四纪红色粘土区氮磷肥交互作用对黑麦草产量和养分累积的影响. 土壤通报, 1995, 26(3): 120-121.

[SS23] 刘经荣, 张美良, 郑群英, 石庆华, 王少先. 不同施N水平对黑麦草产量和品质的效应. 江西农业大学学报, 2003, 25(6): 844-847.

[SS24] 方华舟. 施肥对黑麦草产量和质量的影响. 中国农村小康科技, 2007(8): 79-80, 41.

[SS25] 李文西. 苏丹草—黑麦草轮作制中连续施肥对饲草生长与土壤肥力的影响[D]. 武汉: 华中农业大学, 2010.

[SS26] 李小坤, 李文西, 鲁剑巍, 黄元仿, 刘光文. 施肥对黑麦草产量和氮磷钾养分吸收的影响. 中国农学通报, 2006, 22(7): 331-334.

[SS27] 李小坤, 鲁剑巍, 鲁君明, 陈防, 黄元仿. 磷肥用量对黑麦草产量及经济效益的影响. 草业科学, 2006, 23(10): 18-22.

[SS28] 李小坤, 鲁剑巍, 刘晓伟, 刘威, 潘福霞, 魏云霞. 配方施肥对一年生黑麦草产草量及品质的影响. 草业科学, 2011, 28(9): 1666-1670.

[SS29] 杨娟. 施肥对两种轮作牧草的产量和品质效果及其机理研究[D]. 武汉: 华中农业大学, 2009.

[SS30] 申晓萍, 李仕坚, 朱梅芳, 樊露晓. 刈割次数对南方冬闲田种植黑麦草产量及品质的影响. 安徽农业科学, 2009, 37(35): 17445-17446, 17461.

[SS31] 田松军, 何胜江, 雷荷仙, 付强, 何俊, 赵相勇. 不同施肥水平对冬闲田牧草产量的影响. 贵州畜牧兽医, 2006, 30(4): 3-5.

[SS32] 王柏渊. 施肥和刈割对黔南区轮作多花黑麦草产草量及体外消化的研究[D]. 雅安: 四川农业大学, 2009.

[SS33] 张明均, 张子琴. 多花黑麦草在贵州地区的栽培技术. 农技服务, 2008, 25(4): 97-98.

[SS34] 池银花, 张爱华. 氮肥对禾本科牧草—黑麦草产量和品质的影响. 福建畜牧兽医, 2003, 25(6): 29-30.

[SS35] 池银花, 张爱华. 追施氮肥对多花黑麦草产量和品质的影响. 福建畜牧兽医, 2006, 28(2): 10-11.

[SS36] 辛国荣, 杨中艺, 徐亚幸, 陈三有, 郭仁东. “黑麦草—水稻”草田轮作系统的研究Ⅴ.稻田冬种黑麦草的优质高产栽培技术. 草业学报, 2000, 9(2): 17-23.

[SS37] 陈勇, 罗富成, 甄莉, 毛华明, 陈勇, 左应梅, 朱琳梅, 呼红平. 春播多花黑麦草产量和营养价值变化规律的研究. 黑龙江畜牧兽医, 2007(2): 56-58.

[SS38] 丁成龙, 顾洪如, 许能祥, 程云辉. 不同刈割期对多花黑麦草饲草产量及品质的影响. 草业学报, 2011, 20(6): 186-194.

[SS39] 韩东辰, 田雨, 娄玉杰, 管博. 两种有机肥对黑麦草生长和营养品质的影响. 黑龙江畜牧兽医, 2016(3): 127-129, 132.

[SS40] 李静. 施肥对牧草产量影响的试验报告. 当代畜牧, 2015(5): 92-93.

[SS41] 李小梅, 陈天峰, 孙娟娟, 罗燕, 张新全, 闫艳红. 氮肥水平对多花黑麦草产量和品质的影响. 中国草地学报, 2016, 38(1): 54-60.

[SS42] 潘永年. 一年生黑麦草的氮肥增产效益. 草业科学, 1989, 6(1): 55-57.

[SS43] 王文石, 田雨, 娄玉杰. 施用鸡粪对特高多花黑麦草SPAD值、干草产量和蛋白质含量的影响. 草业学报, 2016, 25(3): 240-244.

[SS44] 夏劲伟. 氮肥对一年生黑麦草产量及品质影响的研究[D]. 南京: 南京农业大学, 2000.

[SS45] 张艾青, 黄尧瑶, 初晓辉, 毕玉芬, 姜华, 段新慧, 任健, 单贵莲. 施氮水平对多花黑麦草生产性能和营养品质的影响. 草业与畜牧, 2014(4): 1-5.

[SS46] 张瑞珍, 张新跃, 何光武, 何丕阳, 王元清. 不同刈割高度对多花黑麦草产量和品质的影响. 草业科学, 2008, 25(8): 68-72.

[SS47] 张文明. 甘肃半干旱地区裸燕麦吸肥特性及施肥效应的研究[D]. 兰州: 甘肃农业大学, 2007.

[SS48] WANG S P, WANG Y F, SCHNUG E, HANEKLAUS S, FLECKENSTEIN J. Effects of nitrogen and sulphur fertilization on oats yield, quality and digestibility and nitrogen and sulphur metabolism of sheep in the Inner Mongolia Steppes of China., 2002, 62(2): 195-202.

[SS49] 蒲珉锴. “红试2号”燕麦新品系品比试验及生产性能评价[D]. 雅安: 四川农业大学, 2012.

[SS50] 才让吉, 王巧玲, 王贵珍, 刘丽, 花立民. 不同播量、氮磷肥互作对高寒牧区燕麦产量和品质的影响. 中国草食动物科学, 2015, 35(5): 30-33, 61.

[SS51] 蒲小鹏. 氮肥对高寒地区燕麦草地生产性能及土壤酶活性的影响[D]. 兰州: 甘肃农业大学, 2005.

[SS52] 拉巴. 施有机肥对青引1号燕麦增效试验研究. 草业与畜牧, 2015(05): 20-21+43.

[SS53] 刘文辉, 周青平, 贾志锋, 梁国玲. 施钾对青引1号燕麦草产量及根系的影响. 植物营养与肥料学报, 2010, 16(2): 419-424.

[SS54] 魏巍, 拉巴, 杨文才, 白玛嘎翁, 周娟娟. 氮磷肥配施对青引1号燕麦产量和品质的影响. 作物杂志, 2016(1): 120-124.

[SS55] 肖小平, 王丽宏, 叶桃林, 张帆, 曾昭海, 胡跃高, 杨光立. 施N量对燕麦“保罗”鲜草产量和品质的影响. 作物研究, 2007(1): 19-21.

Spatial characteristics of Nitrogen and phosphorus Flow in cultivated grassland of China

WEI Zhibiao1,2, BAI Zhaohai2, MA Lin2, ZHANG Fusuo1

(1College of Resources and Environmental Sciences, China Agricultural University/Key Laboratory of Plant-Soil Interactions, Ministry of Education, Beijing 100193;2Center for Agricultural Resources Research, Institute of Genetics and Developmental Biology, Chinese Academy of Sciences/Hebei Key Laboratory of Water-Saving Agriculture/Key Laboratory of Agricultural Water Resources, Chinese Academy of Sciences, Shijiazhuang 050021)

【Objective】With the rapid development of animal husbandry in China and the increasing demand for forage, the adjustment of planting structure has become an inevitable trend. The objective of this study is to quantify the nitrogen (N) and phosphorus (P) spatial flow of cultivated forage, and to provide a scientific basis for increasing grassland productivity. 【Method】Taking alfalfa, ryegrass and oat grass as target forage species, the database of N and P input and output was established, the nutrient balance budgets, nutrient use efficiency and environmental emissions were quantified in cultivated grassland of China by using NUFER model. 【Result】 (1) In 2014, the total input (output) of N in alfalfa, ryegrass and oat grassland of China was 1 547, 236 and 67 Gg, respectively. The N input (output) of per unit area in alfalfa, ryegrass and oat grass grassland of China was 326, 427 and 217 kg N·hm-2, respectively. The total input (output) of P in alfalfa, ryegrass and oat grass grassland of China was 323, 44 and 16 Gg and the P input (output) of per unit area was 49, 18 and 35 kg P·hm-2, respectively. N fixation was the largest input item for alfalfa which accounted for 51% of the total N input into alfalfa grassland. N fertilizer was the largest input item for oat grass and ryegrass which accounted for 93% and 84% of the total N input into ryegrass and oat grass grassland, respectively. P fertilizer was the most important input of P in these kinds of forage. (2) In 2014, the nitrogen use efficiency (NUE) of alfalfa, ryegrass and oat grass grassland was 64%, 93% and 69%, respectively. The phosphorus use efficiency (PUE) was 28%, 77% and 34%, respectively. The NUE and pue varied greatly in different regions. (3) In 2014, the N loss of alfalfa, ryegrass and oat grass grassland per unit area was 23, 4.0 and 9.9 kg N·hm-2and the P loss of them was 2.6, 3.8 and 2.6 kg P·hm-2, respectively. The N and P loss of alfalfa and ryegrass grassland in Southwest China was higher than that in other regions. The N and P loss of oat grass grassland in Western China was higher than that in other regions. (4) In 2014, the N and P accumulation of cultivated grassland were positive except for N accumulation in ryegrass grassland. The accumulation rate of alfalfa grassland in Western China was higher than that in Eastern China. There was no regular regional change of N and P accumulation for ryegrass grassland. The N and P accumulation rate of oat grass grassland in Tibet Plateau was higher than that in other regions. 【Conclusion】In 2014, there was a large gap of N and P input (output) for alfalfa, ryegrass and oat grass grassland in different regions of China. N fixation was the major way for alfalfa to acquire N nutrition, so alfalfa grassland management should reduce the input of N fertilizer. Oat grass and ryegrass mainly get the required nutrients by chemical fertilizer, so it should pay attention to fertilizer application for oat grass and ryegrass grassland. The NUE of target forage species was higher than 60%, but PUE was relatively low. Ryegrass grassland was deficient in N, and it is necessary to apply more nitrogen fertilizer. other cultivated grasslands had surplus of N and P in varying degrees, so the excessive application of fertilizer should be controlled.

cultivated grassland; nitrogen; phosphorus; nutrient flow; nutrient use efficiency; soil nutrient accumulation

2017-07-31；

2017-09-29

国家重点研发计划（2016YFD0800106，2016YFD0200105）、国家自然科学基金面上项目（31572210）、河北省杰出青年基金（D2017503023）、中国科学院百人计划项目

魏志标，E-mail：weizb1993@126.com。

马林，E-mail：malin1979@sjziam.ac.cn