中国天然草地氮磷流动空间特征

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



中国天然草地氮磷流动空间特征

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

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

【目的】定量研究天然草地的氮磷流动空间特征，为优化牧草施肥和提高牧草产量提供科学依据。【方法】建立中国天然草地氮磷养分输入(输出)数据库，利用NUFER模型定量中国天然草地氮磷平衡账户、利用率和环境排放特征。【结果】(1)2013年，全国天然草地氮和磷的输入(输出)总量分别为5 034 Gg N和318 Gg P，单位面积的输入(输出)量分别为19 kg N·hm-2和1.2 kg P·hm-2。氮沉降和畜禽粪尿磷分别占氮和磷输入总量的49%和89%。各区域天然草地氮和磷输入(输出)量变化范围分别为7.0—70 kg N·hm-2和0.12—8.0 kg P·hm-2；(2)2013年，天然草地氮和磷养分利用率分别为105%和191%，各区域间差异很大。中国各地区天然草地的氮利用率变化范围为67%—141%，磷利用率的变化范围为75%—538%；(3)2013年，天然草地氮和磷的环境损失量分别为1.7 kg N·hm-2和0.059 kg P·hm-2，氨挥发和侵蚀分别是氮和磷的主要损失途径。西南和东北地区天然草地氮损失量较多，部分区域的损失量超过8.0 kg N·hm-2；西北地区氮损失量较少，平均不足3.0 kg N·hm-2；青藏高原区氮损失量最少，不足1.0 kg N·hm-2。磷的环境排放空间规律与氮排放相似；(4)2013年，全国天然草地土壤氮和磷的亏缺总量分别为706 Gg N和315 Gg P，单位面积亏缺量分别为2.7 kg N·hm-2和1.2 kg P·hm-2。北方和西南部分地区天然草地的氮土壤累积量为负值，重庆、吉林和辽宁的土壤氮亏损量超过20 kg N·hm-2；西部和西南部分省份天然草地的氮土壤累积量为正值，其中广西和云南的土壤氮累积量超过5.0 kg N·hm-2。除广西和贵州外，其他区域天然草地磷养分均有不同程度的亏缺，重庆天然草地磷的亏缺量最大，为3.7 kg P·hm-2。【结论】2013年，全国天然草地的氮和磷输入量较小，约50%的氮素通过氮沉降输入系统，约90%的磷素通过畜禽粪尿磷输入系统；全国天然草地土壤的氮和磷呈亏缺状态，养分利用率高于100%，当前草地系统不可持续，应注意补施氮磷养分；全国天然草地单位面积氮和磷的环境损失量较小，西南地区天然草地的氮和磷环境损失量大于其他区域。各区域天然草地氮磷流动空间特征差异较大。

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

0 引言

【研究意义】中国天然草地的面积为4亿公顷，是中国国土面积的42%，耕地面积的3.2倍[1-2]。草地不仅具有生产功能，而且还有诸多的生态功能，据估计，通过优化放牧和管理，草地系统可以吸收1.5亿吨CO2当量的温室气体[3-4]。中国天然草地牧草的产量基本不足1 000 kg·hm-2，除宁夏外其他地区的载畜量均有不同程度的超载[5]。天然草地养分流动特征的研究可以为优化草地施肥提供理论依据。【前人研究进展】氮和磷是牧草生产的主要限制因素，其循环对草地系统有重要影响[6-7]。张杰琦等[8]在青藏高原的研究表明，添加氮可以提高高寒草甸的植被高度和盖度，添加磷同样对天然草地有重要影响；闫钟清等[9]研究了人类活动对中国草地系统氮循环关键过程的影响，认为人类活动，如放牧、开垦和火烧，是草原生态系统维持稳定的重要干扰因子，可能会增加硝态氮的流失和N2O的排放。草地系统磷流动方面的研究较少，Chen等[10]分析了中国农业系统的磷流动及其对水体和人类生活的影响，研究指出，尽管中国农业系统每年磷的盈余量为15 kg P·hm-2，但草地生产系统呈现严重的磷亏缺，每年的亏损量接近10 kg P·hm-2，严重影响草地的生产和生态功能。当前相关方面的研究尚有很多不足[11-12]。天然草地基本不施肥且载畜量较高，管理和利用方式不合理，约90%的草地出现不同程度的退化，草原区沙化土地占国土面积的3.9%，局部治理，整体恶化的趋势仍在继续[13-15]。【本研究切入点】天然草地氮磷流动特征研究对于国家政策制定和农户生产均具有重要意义，而当前在这方面的研究尚需完善，本研究通过建立中国各区域的天然草地氮磷输入和输出的数据库，借助NUFER模型，估算各区域天然草地氮磷流动特征、利用率和环境排放特征。【拟解决的关键问题】定量中国不同区域天然草地的养分平衡账户；分析不同区域天然草地氮磷利用率；评估不同区域天然草地的氮磷环境排放。

1 材料与方法

1.1 数据来源

本研究利用NUFER模型计算，模型包括输入模块、输出模块和计算模块[16-17]。输入模块包括变量和参数。化肥用量、牧草播种面积、作物产量等属于模型变量，来源于统计数据。参数例如牧草氮磷含量、养分环境排放因子来源于文献数据。模型输出模块包括氮磷平衡帐户、氮磷利用率和各种形态氮磷环境排放，包括挥发（NH3）、氧化亚氮排放（N2O）、反硝化（N2）、氮磷径流、侵蚀、淋洗等。

1.2 模型计算

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

Igtotal=Igmanure+Igdeposition+Igfixation（1）

式中，Igtotal表示天然草地氮（磷）总输入量（kg·hm-2），Igmanure代表天然草地中粪尿氮（磷）输入（kg·hm-2），Igdeposition代表氮（磷）沉降（kg·hm-2）；Igfixation代表草地生物固氮（kg·hm-2）。具体数值见表1。

Ogtotal=Oggrass+OgNH3+OgN2O+Ogdenitrification+Ogrunoff+Ogerosion+Ogleaching+Ogaccumulation（2）

Oggrass=yield×grassN(P)C/100 （3）

式中，Ogtotal表示天然草地氮（磷）总输出量（kg·hm-2），Oggrass表示天然草地牧草氮（磷）输出量（kg·hm-2），其数值基于牧草产量和牧草氮含量估算；牧草产量（kg·hm-2）和牧草氮（磷）含量（%）分别由yield和grassN（P）C表示，各地区的具体量见表2；OgNH3表示氨挥发，OgN2O表示氧化亚氮排放，Ogdenitrification表示反硝化，Ogrunoff表示氮（磷）径流，Ogerosion表示氮（磷）侵蚀，Ogleaching表示氮（磷）淋洗，Ogaccumulation表示天然草地氮（磷）土壤累积（kg·hm-2），以上数据由NUFER模型模拟得到[16-17]，该模型可模拟食物链养分流动情况。

N(P)UEg=（Oggrass/Igtotal）×100% （4）

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

表1 2013年中国不同区域天然草地氮磷养分输入量

表中数据为2013年的值，基于发表的文章。草地生物固氮量基于草产量估算，青藏高原区天然草地生物固氮量占总吸氮量的9.0%，其他区域占20%

Date used derived from published articles and based on the results of 2013. N fixation of natural grassland was calculated based on grass yield. N fixation of natural grassland in Tibet Plateau accounts for 9.0% of total N uptake, while other regions account for 20%

表2 2013年中国不同区域天然草地牧草产量、氮磷含量和牧草氮磷携出量

天然草地产量数据源于全国草原监测报告；牧草氮和磷含量通过文献总结获得（附录）

Natural grassland yield data used derived from The National Grassland Report. The content of N and P of forage acquired by summarizing published articles (appendix)

2 结果

2.1 2013年中国天然草地氮磷流动特征

2013年中国天然草地氮总输入（输出）量为5 034 Gg N，氮沉降是氮输入天然草地的一种主要途径，占总输入量的49%，牧草携出是天然草地氮输出的最主要方式，占总输出量的105%。全国天然草地磷总输入（输出）量为318 Gg P，动物粪尿是天然草地磷输入的一种主要途径，占总输入量的89%，另外有11%的磷以沉降的形式输入到天然草地，牧草携出是天然草地磷输出的最主要方式，占总输出量的191%。

2013年中国单位面积天然草地氮总输入量为19 kg N·hm-2，牧草携出氮为20 kg N·hm-2，天然草地氮利用率为105%，以氨挥发等环境损失的方式进入到大气或水体中的量为1.7 kg N·hm-2，天然草地土壤氮累积量为-2.7 kg N·hm-2。单位面积天然草地磷输入（输出）量较少，为1.2 kg P·hm-2，动物粪尿是天然草地磷输入的一种主要途径（1.1 kg P·hm-2），以沉降的形式输入到单位面积天然草地的量为0.14 kg P·hm-2。在单位面积天然草地磷的各输出项中，牧草携出量为2.3 kg P·hm-2，以径流和侵蚀等环境损失的方式进入到地表和地下水体中的量分别为0.029和0.061 kg P·hm-2，磷利用率为191%，草地土壤磷累积量为-1.2 kg P·hm-2（表3）。

2.2 2013年中国各地区天然草地氮磷平衡账户

2013年中国各地区单位面积天然草地氮输入量变化范围为7.0—70 kg N·hm-2，南方和东北地区天然草地氮输入量较大，内蒙、新疆、甘肃、西藏和青海5大牧区的氮输入量均较小，其中青藏高原区的输入量最少，西藏和青海的氮输入量为7.4和7.0 kg N·hm-2。各地区天然草地的氮输出量变化范围为6.3—91 kg N·hm-2，氮输入量越大的地区牧草携出的养分越多，重庆通过牧草携出的氮养分最多，为91 kg N·hm-2（图1-a）。

2013年中国各地区单位面积天然草地磷输入量变化范围为0.12—8.0 kg P·hm-2，西南和东北地区磷输入量较高，输入量最高的省份为贵州（8.0 kg P·hm-2）；5大牧区单位面积的输入量较低，最低的省份为青海，输入量仅为0.21 kg P·hm-2。磷输入量较多的区域一般牧草携出量也偏高，牧草携出量分布在0.61—11 kg P·hm-2（图1-b）。

表3 2013年全国天然草地氮磷流动特征

图1 2013年中国不同地区天然草地氮磷平衡账户

2.3 2013年中国各地区天然草地氮磷利用率

2013年中国各地区天然草地的氮利用率变化范围为67%—141%，由于青海、重庆、吉林、辽宁、内蒙古和甘肃的牧草携出量大于当地天然草地的氮总输入量，故上述区域的氮利用率超过100%，其中青海天然草地的氮利用率最高，为141%；其余9个地区的牧草携出量小于氮的总输入量，氮利用率小于100%，云南的氮利用率最低（67%）（表4）。

2013年中国地区天然草地的磷利用率高于氮利用率，变化范围为75%—538%。青海的磷利用率最高，达到538%，主要原因是当地天然草地的磷输入量非常低。内蒙和甘肃的磷利用率也较高，超过200%；仅广西、贵州和云南的磷利用率小于100%（表4）。

2.4 2013年中国各地区天然草地氮磷环境损失

西南和东北地区天然草地氮损失较多，贵州、云南、广西和吉林的氮损失量分别为10、9.5、9.4和7.5 kg N·hm-2；西北地区的氮损失量较少，平均不足3.0 kg N·hm-2；青藏高原区氮损失量最少，不足1.0 kg N·hm-2。氨挥发是氮损失的主要途径，辽宁、吉林、贵州、重庆和广西等地区通过氨挥发损失的氮大于5.0 kg N·hm-2（图2-a）。磷的环境排放与氮排放规律相似，西南和东北地区的排放量较多。通过淋洗损失的磷很少，通过侵蚀损失的量较多，但基本不超过1.0 kg P·hm-2（图2-b）。

2.5 2013年中国各地区单位面积天然草地氮磷土壤累积

北方和部分西南地区的氮土壤累积量为负值，重庆、吉林和辽宁的土壤氮亏损量超过20 kg N·hm-2。西部和西南部分地区的氮土壤累积量为正值，其中广西和云南两地的土壤氮累积量超过5.0 kg N·hm-2。西藏、宁夏和贵州的氮土壤累积量介于0—2.0 kg N·hm-2，土壤氮基本处于平衡状态（图3）。

除广西和贵州外，其余地区均有不同程度的磷亏缺。其中新疆和云南草地土壤磷亏缺量不足0.1 kg P·hm-2，土壤磷基本处于平衡状态，重庆天然草地磷的亏缺量最大，为3.7 kg P·hm-2。广西 和贵州天然草地的土壤磷累积量分别为1.3和0.26 kg P·hm-2（图3）。

图2 2013年中国不同地区天然草地氮磷损失

图3 2013年中国不同地区天然草地氮磷土壤累积

表4 2013年中国不同地区天然牧草氮磷利用率

3 讨论

3.1 中国天然草地氮磷流动特征

本研究发现，中国天然草地氮磷养分账户均为负值，草地土壤养分略有亏缺，草地严重退化与此相关。草地氮磷的利用率超过100%，青海等部分地区的磷利用率甚至超过200%。Gu等[21]研究表明，中国草地系统氮利用率为59%—95%，与本研究结果存在差异。主要有两方面原因，首先Gu等的研究对象包括天然和栽培草地两种类型，由于化肥的投入，栽培草地的养分利用率低于天然草地。另外，限于数据的问题，本研究认为天然草地的产量等于草地的氮磷移除量，故本研究的草地养分利用率较高。下一步的研究中，可以区分放牧区和没有采食的地区或禁牧区，将这部分草地的牧草输出设为0，以减小研究结果的误差。

中国天然草地土壤氮亏缺量为2.7 kg N·hm-2，且各区域间差异很大，这与前人研究结果较为一致。Giese等[22]在内蒙古的计算结果表明，放牧区氮累积量为负值，禁牧区的氮累积量为正值；Chen等[23]认为中国天然草地的土壤氮亏缺量为-31.7 kg N·hm-2，养分输出方面，牧草携出氮量为80 kg N·hm-2，粪尿氮的输入量仅为13 kg N·hm-2，两方面的原因导致其计算的氮账户为负值。在本研究所涉及的地区中，吉林省牧草携出量与Chen等计算的牧草氮带走量相当，而吉林天然草地的氮累积量为-24 kg N·hm-2，略高于Chen等的估算，原因是本研究中吉林的粪尿输入量为29 kg N·hm-2，高于Chen等的估算。

动物粪尿归还草地量占天然草地氮总输入量的33%（表3），Bouwman等[24]研究表明，粪尿施用量占亚洲农用地氮总输入量的31%。Bouwman等的研究对象包括耕地和草地，而本文的研究对象是天然草地。耕地和人工草地的粪尿施用比例远小于天然草地[11,17]，本研究的粪尿归还天然草地的比例实际低于Bouwman等的研究结果。其一是因为中国的氮沉降量较大（约占天然草地氮总输入量的50%）；其二是因为中国牲畜粪尿环境损失量较高[7]。

本研究表明，中国天然草地的土壤磷亏缺量为1.2 kg P·hm-2。Sattari等[12]对亚洲草原系统磷流动特征分析的结果表明，草原磷亏缺量为3.7 kg P·hm-2，低于本研究计算的结果。Sattari等认为每年牧草携出的磷为6.6 kg P·hm-2, 按照磷含量0.22%计算，每年牧草的产量为3 000 kg·hm-2，而中国主要牧区牧草实际产量不足1 000 kg·hm-2，因此，Sattari等对亚洲草原区牧草的产量有所高估，进而低估了草原磷累积量。Chen等[10]使用模型估算了中国草地磷流动特征，表明中国草地的磷累积量为-4.7 kg P·hm-2，其估算的累积量结果更低。原因除高估了中国天然草地的产量外（3 441 kg·hm-2），他们还认为仅有16%的动物粪尿归还到草地。事实上中国牧区多为天然放牧，除了环境损失和牧民取暖消耗掉部分动物排泄物外，多数粪尿重新归还草地。

3.2 中国天然草地养分管理优化策略

基于本研究结果，对于中国天然草地养分管理的优化策略应该包括以下两个方面：

其一，中国天然草地应适当的补施氮肥和磷肥。当前天然草地退化严重，不仅制约了畜牧业的发展，还严重影响草地的生态功能[5,25-26]。以往的研究认为气候和放牧强度是影响草地生产力的最主要因素[27-28]，本研究发现中国天然草地土壤氮磷整体上略有亏缺。长期以来，中国天然草地基本不施肥，且载畜量较高，草地土壤养分的管理不合理可能是天然草地退化的一个更直接原因。

其二，不同牧草种植区应采用不同的氮磷养分管理策略，西南和西部地区天然草地土壤氮磷养分基本平衡，草地退化并非由土壤养分管理不合理所导致。辽宁、吉林和重庆天然草地土壤氮亏缺严重，每年应补施氮20—30 kg N·hm-2。草原磷账户方面，除广西和贵州外，全国草原磷账户均为负值。为了增加草地的生产力，应适当补施磷肥。内蒙和重庆应在控制载畜量的基础上，补施部分磷肥，防止草地退化。

3.3 模型的不确定性

本研究利用NUFER模型模拟了中国天然草地氮磷流动特征，但是由于模型的算法和参数的获取的缺陷，也存在一些不确定性，主要表现在以下几方面：本研究中氮沉降的数据源于文献[18]。大气按沉降的形式主要分为干沉降和湿沉降，其中报道最多的是湿沉降。由于下垫面土地利用方式和测定方法的不同，现有研究关于氮沉降量的报道差异很大。其中，对氮湿沉降的研究比较系统，使用该数值作为本研究的沉降量虽然会对总沉降量造成低估，但各区域的相对大小可以保证。

磷沉降方面，磷沉降同样分为干沉降和湿沉降，沉降形式主要为无机磷。对湿沉降的研究相对比较系统，因此本研究所用的数据来自Zhu等[19]2016年对全国各地区磷湿沉降的研究结果。当前关于磷沉降的研究较少，且报道的数值变化较大，本研究中，全国的平均值为0.21 kg P·hm-2。近年来关于中国北方地区特别是内蒙古沙尘的报道逐渐增多，大量的沙尘会导致磷以干沉降的形式进入草原系统[29-31]。但不能忽略的是，由于利用不合理，放牧强度过大等原因，内蒙古草地的沙化问题不断加剧[32]。本研究中，内蒙地区的草原磷沉降为0.11 kg P·hm-2（湿沉降量），虽然该数值未包含干沉降的部分，但是笔者认为干沉降进入系统的磷和通过风蚀等途径向系统外输出的磷在总量上可以抵消，因此不予考虑。

天然草地粪尿输入量根据Bai等[7]的算法获得，利用中国各区域牛、羊、马等草食动物的养殖量计算粪尿产生量，再基于粪尿归还草地的比例获取天然草地粪尿输入量。养分的环境排放方面，当前关于天然草地环境排放参数的系统报道很少，本研究根据使用NUFER模型中的算法，根据植被覆盖等因子，计算相应的环境排放量[16-17]。

4 结论

（1）2013年，全国天然草地的氮和磷输入量较小，氮沉降和畜禽粪尿磷分别是最主要的氮和磷输入项；（2）2013年，全国天然草地氮和磷的养分利用率高于100%，各区域氮和磷利用率差异较大；（3）2013年，全国天然草地单位面积氮和磷的环境损失量较小，各区域的环境损失量差异较大，西南地区天然草地的氮和磷环境损失量大于其他区域；（4）2013年，全国天然草地土壤氮和磷呈亏缺状态，当前草地系统不可持续，应注意补施氮磷养分。

[1] 方精云, 白永飞, 李凌浩, 蒋高明, 黄建辉, 黄振英, 张文浩, 高树琴. 我国草原牧区可持续发展的科学基础与实践. 科学通报, 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)

[2] 沈海花, 朱言坤, 赵霞, 耿晓庆, 高树琴, 方精云. 中国草地资源的现状分析. 科学通报, 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)

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

[4] 何念鹏, 韩兴国, 于贵瑞. 内蒙古放牧草地土壤碳固持速率和潜力. 生态学报, 2012, 32(3): 844-851.

HE N P, HAN X G, YU G R. Soil carbon sequestration rates and potential in the grazing grasslands of Inner Mongolia., 2012, 32(3): 844-851. (in Chinese)

[5] 张新时, 唐海萍, 董孝斌, 李波, 黄永梅, 龚吉蕊. 中国草原的困境及其转型. 科学通报, 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)

[6] 王敬国, 林杉, 李保国. 氮循环与中国农业氮管理. 中国农业科学, 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)

[7] BAI Z H, MA L, JIN S Q, MA W Q, VELTHOF G L, OENEMA O, LIU L, CHADWICK D, ZHANG F S. Nitrogen, phosphorus, and potassium flows through the manure management chain in China., 2016, 50(24): 13409-13418

[8] 张杰琦, 李奇, 任正炜, 杨雪, 王刚. 氮素添加对青藏高原高寒草甸植物群落物种丰富度及其与地上生产力关系的影响. 植物生态学报, 2010, 34(10): 1125-1131.

ZHANG J Q, LI Q, REN Z W, YANG X, WANG G. Effect of nitrogen addition on species richness and relationship between species richness and aboveground productivity of alpine meadow of the Qinghai- Tibetan Plateau, China., 2010, 34(10): 1125-1131. (in Chinese)

[9] 闫钟清, 齐玉春, 董云社, 彭琴, 孙良杰, 贾军强, 曹丛丛, 郭树芳, 贺云龙. 草地生态系统氮循环关键过程对全球变化及人类活动的响应与机制. 草业学报, 2014, 23(6): 279-292.

YAN Z Q, QI Y C, DONG Y S, PENG Q, SUN L J, JIA J Q, CAO C C, GUO S F, HE Y L. Nitrogen cycling in grassland ecosystems in response to climate change and human activities., 2014, 23(6): 279-292. (in Chinese)

[10] CHEN M, CHEN J, SUN F. Agricultural phosphorus flow and its environmental impacts in China., 2008, 405(1/3): 140-152.

[11] GU B J, JU X T, CHANG J, GE Y, VITOUSEK P M. Integrated reactive nitrogen budgets and future trends in China., 2015,112(28): 8792-8797.

[12] SATTARI S Z, BOUWMAN A F, MARTINEZ RODRIGUEZ R, BEUSEN A H, VAN ITTERSUM M K. Negative global phosphorus budgets challenge sustainable intensification of grasslands., 2016, 7: 10696.

[13] 周道玮, 姜世成, 王平. 中国北方草地生态系统管理问题与对策. 中国草地, 2004, 26(1): 58-65.

ZHOU D W, JIANG S C, WANG P. The problems and strategies for grassland management in Northern China., 2004, 26(1): 58-65. (in Chinese)

[14] NIU K, HE J S, LECHOWICZ M J. Grazing-induced shifts in community functional composition and soil nutrient availability in Tibetan alpine meadows., 2016, 53(5): 1554-1564.

[15] LIU M, LIU G H, WU X, WANG H, CHEN L. Vegetation traits and soil properties in response to utilization patterns of grassland in Hulun Buir City, Inner Mongolia, China., 2014, 24(4): 471-478.

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

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

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

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

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

[20] YANG B, QIAO N, XU X, OUYANG H. Symbiotic nitrogen fixation by legumes in two Chinese grasslands estimated with the15N dilution technique., 2011, 91(1): 91-98.

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

[22] GIESE M, BRUECK H, GAO Y Z, LIN S, STEFFENS M, KOGEL-KNABNER I, GLINDEMANN T, SUSENBETH A, TAUBE F, BUTTERBACH-BAHL K, ZHENG X H, HOFFMANN C, BAI Y F, HAN X G. N balance and cycling of Inner Mongolia typical steppe: a comprehensive case study of grazing effects., 2013, 83(2): 195-219.

[23] CHEN M, SUN F, SHINDO J. China’s agricultural nitrogen flows in 2011: Environmental assessment and management scenarios., 2016, 111: 10-27.

[24] BOUWMAN L, GOLDEWIJK K K, VAN DER HOEK K W, BEUSEN A H W, VAN VUUREN D P, WILLEMSS J, RUFINO M C, STEHFEST E. Exploring global changes in nitrogen and phosphorus cycles in agriculture induced by livestock production over the 1900-2050 period., 2013, 110(52): 20882-20887.

[25] 李小坤, 鲁剑巍, 陈防. 牧草施肥研究进展. 草业学报, 2008, 17(2): 136-142.

LI X K, LU J W, CHEN F. Primary study on fertilizer application of forage., 2008, 17(2): 136-142. (in Chinese)

[26] 倪珍. 松嫩平原退化草地土壤跳虫多样性恢复潜力及管理对策研究[D]. 长春: 吉林大学, 2015.

NI Z. Recovery of soil Collembola biodiversity and its service to management policy-making in degraded grassland in Songnen Plain, China[D]. Changchun: Jilin University, 2015. (in Chinese)

[27] 张鲜花, 朱进忠, 孙宗玖, 靳瑰丽, 郑伟, 古伟容. 放牧强度对草地牧草现存量及养分动态的影响. 草业科学, 2014, 31(1): 116-124.

ZHANG X H, ZHU J Z, SUN Z J, JIN G L, ZHENG W, GU W R. Influence of grazing intensity on the aboveground biomass and nutrient dynamics of community., 2014, 31(1): 116-124. (in Chinese)

[28] 杨凯. 气候变化对藏北地区草地生产力的影响模拟[D]. 北京: 中国农业科学院, 2007.

YANG K. Stimulation study of climate change impacts on grassland productivity in Northern Tibet[D]. Beijing: Chinese Academy of Agricultural Sciences, 2007. (in Chinese)

[29] DONG W, ZHENG J, GUO Q. Particle-size speciation of Pu isotopes in surface soils from Inner Mongolia (China) and its implications for Asian dust monitoring., 2017, 120: 133-136.

[30] PARK S U, JU J W, LEE I H, JOO S J. A parameterization of dust emission (PM10) fluxes of dust events observed at Naiman in Inner Mongolia using the monitored tower data., 2016, 140: 420-431.

[31] SONG H Q, ZHANG K S, PIAO S L, WAN S Q. Spatial and temporal variations of spring dust emissions in northern China over the last 30 years., 2016, 126: 117-127.

[32] TANG Z, AN H, DENG L, WANG Y, ZHU G, SHANGGUAN Z. Effect of desertification on productivity in a desert steppe., 2016, 6: 27839.

（责任编辑 岳梅）

附录：天然草地氮磷养分含量数据来源

[S1] CHENG X, AN S, CHEN J, LI B, LIU Y, LIU S. Spatial relationships among species, above-ground biomass, N, and P in degraded grasslands in Ordos Plateau, northwestern China., 2007, 68(4): 652-667.

[S2] GONG X Y, CHEN Q, DITTERT K, TAUBE F, LIN S. Nitrogen, phosphorus and potassium nutritional status of semiarid steppe grassland in Inner Mongolia., 2011, 340(1/2): 265-278.

[S3] GONG Z, KAWAMURA K, ISHIKAWA N, INABA M, ALATENG D. Estimation of herbage biomass and nutritive status using band depth features with partial least squares regression in Inner Mongolia grassland, China., 2016, 62(1): 45-54.

[S4] QI J, NIE Z N, JIAO T, ZHANG D G. Phosphorus and defoliation interact and improve the growth and composition of the plant community and soil properties in an alpine pasture of Qinghai-Tibet Plateau., 2015, 10(10): e0141701.

[S5] WU G, WEI J, DENG H B, ZHAO J Z. Nutrient cycling in an alpine tundra ecosystem on Changbai mountain, northeast China., 2006, 32(2): 199-209.

[S6] XU S J, FAN X Y, WANG L L, ZHANG X F, AN L Z. The patterns of nitrogen and phosphorus stoichiometry across communities along altitudinal gradients in Qilian Mountains, China., 2015, 62: 58-65.

[S7] 曹丽霞. 丛枝菌根对退化羊草典型草原修复效应的研究[D]. 呼和浩特: 内蒙古农业大学, 2015.

[S8] 陈军强, 张蕊, 侯尧宸, 马丽娜, 丁路明, 龙瑞军, 尚占环. 亚高山草甸植物群落物种多样性与群落C、N、P生态化学计量的关系. 植物生态学报, 2013, 37(11): 979-987.

[S9] 樊江文, 张良侠, 张文彦, 钟华平. 中国草地样带植物氮磷元素空间格局及其与气候因子的关系. 草地学报, 2014(1): 1-6.

[S10] 高金龙, 侯尧宸, 白彦福, 孟宝平, 杨淑霞, 胡远宁, 冯琦胜, 崔霞, 梁天刚. 基于高光谱数据的高寒草甸氮磷钾含量估测方法研究——以青海省贵南县及玛沁县高寒草甸为例. 草业学报, 2016, 25(3): 9-21.

[S11] 侯金文, 谢边疆, 乔志强, 刘贵山, 闫晓婷, 张玉萍. 2007年包头市达茂旗天然草场牧草营养成分浅析. 农村牧区机械化, 2008(2): 15-16.

[S12] 靳发兰. 不同利用方式下山地高草草甸草地群落变化趋势[D]. 乌鲁木齐: 新疆农业大学, 2007.

[S13] 母悦, 耿元波. 内蒙古羊草草原植物营养元素的含量特征. 生态环境学报, 2015, 24(7): 1118-1124.

[S14] 秦建蓉, 马红彬, 王丽, 虎巧能, 沈艳, 许冬梅. 宁夏荒漠草原植物群落特征对不同轮牧开始时间的响应. 草业科学, 2016, 33(5): 963-971.

[S15] 王艳萍. 沙图仕特性及可可西里牧草-土壤中主要养分含量的研究[D]. 兰州: 甘肃农业大学, 2005.

[S16] 吴发莉. 不同区域和季节对高寒天然牧草品质的影响[D]. 雅安: 四川农业大学, 2014.

[S17] 吴海艳, 马玉寿, 董全民, 孙小弟, 施建军, 王彦龙, 盛丽. 黄河源区藏嵩草沼泽化草甸地上生物量及营养季节动态研究. 草业科学, 2009, 26(1): 8-12.

[S18] 乌鲁木山·布仁巴依尔. 不同放牧强度对克氏针茅典型草原土壤、植被的影响[D]. 呼和浩特: 内蒙古农业大学, 2013.

[S19] 徐长林. 坡向对青藏高原东北缘高寒草甸植被构成和养分特征的影响. 草业学报, 2016, 25(4): 26-35.

[S20] 薛艳锋, 郝力壮, 刘书杰. 玉树州藏嵩草草地牧草营养价值评定与营养载畜量. 草业科学, 2015, 32(10): 1660-1667.

[S21] 杨阔, 黄建辉, 董丹, 马文红, 贺金生. 青藏高原草地植物群落冠层叶片氮磷化学计量学分析. 植物生态学报, 2010, 34(1): 17-22.

[S22] 杨霞, 王珍, 运向军, 卫智军. 不同降雨年份和放牧方式对荒漠草原初级生产力及营养动态的影响. 草业学报, 2015, 24(11): 1-9.

[S23] 张均. 西藏那曲地区不同月份草地牧草营养价值评定及绒山羊营养补饲研究[D]. 雅安: 四川农业大学, 2005.

[S24] 宗子贺, 李俊有, 赵红. 赤峰市2011年天然草场牧草营养成分监测评估. 内蒙古农业科技, 2012(3): 110.

[S25] LU X Y, YAN Y, FAN J H, CAO Y Z, WANG X D. Dynamics of above- and below-ground biomass and C, N, P accumulation in the alpine steppe of Northern Tibet., 2011, 8(6): 838-844.

[S26] YAO Z Y, ZHAO C Y, YANG K S, LIU W C, LI Y, YOU J D, XIAO J H. Alpine grassland degradation in the Qilian Mountains, China - A case study in Damaying Grassland., 2016, 137: 494-500.

[S27] ZHANG J, WANG X J, WANG J P, WANG W X. Carbon and nitrogen contents in typical plants and soil profiles in Yanqi Basin of Northwest China., 2014, 13(3): 648-656.

[S28] 阿拉坛格日乐. 浅谈2013年呼伦贝尔市天然牧草营养成分. 现代农业, 2015(4): 97.

[S29] 董全民, 赵新全, 李青云, 马玉寿, 李有福, 李发吉. 小嵩草高寒草甸的土壤养分因子及水分含量对牦牛放牧率的响应Ⅱ冬季草场土壤营养因子及水分含量的变化. 土壤通报, 2005, 36(4): 493-500.

[S30] 郭春华, 张均, 王康宁, 意西多吉, 吴玉江, 索朗达. 高寒草地生物量及牧草养分含量年度动态研究. 中国草地学报, 2007, 29(1): 1-5.

[S31] 黄晓霞, 韩京萨, 刘全儒, 江源, 和克俭. 小五台亚高山草甸植物化学成分研究. 草业科学, 2009, 26(3): 1-8.

[S32] 金花. 基于3S技术支持的草地营养与载畜量评价研究[D]. 呼和浩特: 内蒙古农业大学, 2008.

[S33] 李超, 安沙舟, 周小丽, 范天文, 李海. 昭苏马场季节草地牧草经济类群营养成分初步分析. 新疆农业科学, 2012, 49(9): 1681-1687.

[S34] 李赟. 长期围封对亚高山草地土壤和植被的影响[D]. 乌鲁木齐: 新疆农业大学, 2009.

[S35] 刘兴波. 天然牧草养分对草地利用强度与加工方式的响应[D]. 呼和浩特: 内蒙古农业大学, 2015.

[S36] 陆阿飞. 三江源区不同地区天然草地牧草营养成分分析. 黑龙江畜牧兽医, 2015(15): 139-141.

[S37] 那守海, 孙国荣, 阎秀峰, 李景信, 李凤芹, 李艳波. 松嫩草地6种牧草地上部分养分含量与土壤因子的关系. 黑龙江畜牧兽医, 1996(3): 1-4.

[S38] 孙鹏飞, 崔占鸿, 刘书杰, 柴沙驼, 郝力壮, 王迅. 三江源区不同季节放牧草场天然牧草营养价值评定及载畜量研究. 草业学报, 2015, 24(12): 92-101.

[S39] 塔娜, 那日苏, 王海, 拾涛. 荒漠草原不同经营规模牧户草场地上生物量和牧草养分含量比较. 中国草地学报, 2013, 35(5): 156-162.

[S40] 王建林, 钟志明, 王忠红, 余成群, 张宪洲, 胡兴祥, 沈振西, 大次卓嘎. 青藏高原高寒草原生态系统植被磷含量分布特征及其影响因素. 草地学报, 2014, 22(1): 27-38.

[S41] 王鑫, 胡玉昆, 热合木都拉.阿迪拉, 李凯辉, 范永刚, 柳妍妍. 高寒草地主要类型土壤因子特征及对地上生物量的影响. 干旱区资源与环境, 2008, 22(3): 196-200.

[S42] 王岩春. 阿坝县国家退牧还草工程项目区围栏草地恢复效果的研究[D]. 雅安: 四川农业大学, 2007.

[S43] 许涛. 玛曲县高寒人工草地植被群落和土壤特性研究[D]. 兰州: 甘肃农业大学, 2013.

[S44] 许志信, 巴图朝鲁, 段淳清. 草甸草原六种牧草贮藏养分含量变化规律的研究. 中国草地, 1993(6): 22-25.

[S45] 阎萍, 潘和平. 不同季节牧草营养成分与牦牛血液激素含量变化的研究. 甘肃农业大学学报, 2004, 39(1): 50-52.

[S46] 张伟娜. 不同年限禁牧对藏北高寒草甸植被及土壤特征的影响[D]. 北京: 中国农业科学院, 2015.

spatial Characteristics of Nitrogen and phosphorus Flow in natural 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】The objective of this study is to quantify the spatial characteristics of nitrogen (N) and phosphorus (P) flow in natural grassland of China, and to provide a scientific basis for optimizing forage fertilization and increasing forage yield. 【Method】The input and output database of N and P in natural grassland of China was established, the nutrient balance budget, nutrient use efficiency and environmental emission characteristics in natural grassland of China were quantified by using NUFER model. 【Result】(1) In 2013, the total input of N and P in natural grassland of China was 5 034 Gg N and 318 Gg P, respectively. The N and P input per unit area in natural grassland was 19 kg N·hm-2and 1.2 kg P·hm-2, respectively. The N deposition and the manure P applied accounted for 49% and 89% of total N and P input. The range of input (output) of N and P in natural grassland in different regions was 7.0-70 kg N·hm-2and 0.12-8.0 kg P·hm-2; (2) In 2013, the use efficiency of N and P in natural grassland of China was 105% and 191%, respectively. There were great differences among different regions. The range of use efficiency of N and P was 67%-141%, 75%-538%, respectively; (3) In 2013, the environmental loss of N and P in natural grassland of China was 1.7 kg N·hm-2and 0.059 kg P·hm-2, respectively. Ammonia emission and erosion were the main loss pathways of N and P in natural grassland, respectively. The N loss in natural grassland of Southwest and Northeast China was over 8.0 kg·hm-2. The N loss of Northwest China was relatively low, which was less than 3.0 kg·hm-2on average. The N loss in Tibet plateau was the lowest, which was less than 1.0 kg·hm-2. Spatial regularity of environmental P emission was similar to that of N emission; (4) In 2013, the total N and P deficit in natural grassland of China was 706 Gg N and 315 Gg P, respectively. The deficit per unit area of N and P was 2.7 kg N·hm-2and 1.2 kg P·hm-2, respectively. The value of soil N accumulation in North and Southwest China was negative, and the soil N deficit in Chongqing, Jilin and Liaoning was more than 20 kg N·hm-2. The value of soil nitrogen accumulation in West and Southwest China was positive, and the soil N accumulation in Guangxi and Yunnan was more than 5.0 kg N·hm-2. Soil P accumulation in natural grassland of China was different from N accumulation. With the exception of Guangxi and Guizhou, all of other regions had a negative soil P accumulation. The P deficit of natural grassland in Chongqing was the largest, which was -3.7 kg P·hm-2.【Conclusion】In 2013, the total input of N and P in natural grassland of China was relatively low, about 50% of the N was input through the N deposition and 90% of the P was input by manure P applied. The soil accumulation of N and P in natural grassland of China was negative in 2013, and the N and P use efficiency was higher than 100%. The current grassland system was unsustainable and applying moderate chemical fertilizers was sensible. In 2013, the environmental loss of N and P in natural grassland of China was relatively low. The nutrient loss in Southwest China was larger than that in other regions. The spatial flow characteristics of N and P in natural grassland varied greatly among different regions.

natural 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