我国农业氨排放估算方法研究进展

发布时间：2019-10-25

　　摘要: 农业大量施用氮肥以及持续扩大的禽畜养殖业是我国氨污染的最大来源。近年来大气环境问题备受关注，氨排放研究的重要性日益凸显，如何客观、科学定量的评估我国区域氨排放的问题尤为重要。通过检索已报道的国内外氨排放估算的研究进展，对我国的氨排放研究进行梳理，比对了氨排放主要估算方法的特点，对其所使用数据类型、获取途径，参数的定量取值方法及不确定性产生等方面进行了分析。针对国内氨排放估算存在计算方法单一，排放因子本地化不足等问题，提出了进一步改善的意见和建议。研究结果以期为我国做好氨排放控制基础研究，开展控制技术试验，制定相关政策文件，加强政府引导和扶持等提供科学依据。

　　关键词: 农业氨排放; 估算方法; 排放因子; 不确定性分析

　　农牧业源是农业活动直接排放氨( NH3 ) 的排放源的统称[1]，包括畜禽养殖、化肥施用、生物质燃烧、秸秆堆肥等方面。在实际研究农业氨排放过程中，往往也包括了人体粪便这一排放源[2]。已有研究表明，我国大气氨排放主要来自化肥施用与畜禽养殖，两者排放量之和占人为源氨排放总量比值达 80%[3]。农业生态系统的氨排放是全球氮素循环的重要组成部分，在使农作物营养物质大量流失的同时，也对环境产生了重要的影响[4]。作为大气中唯一的常见气态碱，NH3 易溶于水，能与大气中硫酸气溶胶能够形成( NH4 ) 2 SO4 或 NH4HSO4 [5]，这些二次颗粒物的产生对大气 PM2.5污染和霾的形成有着重要影响。高层大气中氨参与了一系列自由基反应，氨的排放量也存在加剧温室效应的可能性[6]。此外，氨在土壤酸化及水体富营养化上也有着直接或间接的影响[5]。

　　从 20 世纪末开始，国内外学者对 NH3排放清单及其对大气污染影响的研究越来越多，也出现较多关于农业氨排放的研究工作[4，7-8]。在欧洲及美国，畜禽养殖和氮肥施用的 NH3排放量占总排放量的 80%—90%，在大部分亚洲国家二者则占到总量的 77%左右[5]。同时，相关行政部门建立了包括 NH3在内的排放清单[9]。具有代表性的氨排放研究中，Paina 等[10]采用排放因子法估算了英国农业氨排放量，该地区全年农业氨排放量为 197 Gg，畜禽养殖与化肥施用氨排放分别占排放总量的 31%和 16%。在实际过程中，面对较为复杂的生态系统氮循环，模型法在综合考虑氨的排放、迁移、转化过程时，具有一定优势。目前，国外使用广泛的模型估算法由早期的排放因子法发展而来。英国开发的国家氨减排措施评价体系( Nation Ammonia Reduction Strategy Evaluation System，NARSES) ，是一个用于估算农业氨排放规模、时空分布规律以及检测相关政策方案实行可能性的模型。1991 年农业氨排放被加入其中的区域空气污染信息和模拟模型( Regional Air Pollution Information and Simulation Model，RAINS Model) ，由国际应用系统分析学会( IIASA) 开发[9]。Klimont [11]运用 RAINS 模型对 1990 年和 1995 年中国氨排放总量进行估算，结果显示，1990 与 1995 年中国氨排放总量为 970 万 t 和 1170 万 t，预计到 2030 年 NH3的排放量将增加到近 2000 万 t。氨排放主要贡献来自氮肥施用和牲畜，分别占 90 年代总排放量的 52%和 41%。化肥施用氨排放的比例预计在 2030 年将增加到约 61%，而牲畜的份额则下降到 33%。Streets 等[12]基于 RAINS 模型，参考了 Klimont 等[13]的计算方法，估算 2000 年的中国 NH3 排放量为 13. 6Tg，其中 50% 来自化肥施用，占农业氨排放的 88%。丹麦的氨排放清单体模型( Danmark Ammonia Emission Inventory Model，DanAm) ，在建立各排放源排放因子时，考虑了季节因素对排放因子的影响，同时化肥氮施用的氨排放因子为综合试验结果与经验值得到。除此之外，由政府间气候专门委员会 ( Intergovernmental Panel on Climate Change，IPCC) 提供的 IPCC 方法[14]( 活动水平数据模型) 是目前国际上应用最广泛的模型，其划分了详细的氨排放源，给出了各种氨排放源和氨排放估算的指导方法，同时提供了大量全球各地可以选用的默认参数及排放因子[7]。目前，欧洲国家定期报告氨排放量估算，并承诺按照规定的路径实现国家排放限额[1，15]。在美国，一些报告中有排放要求[16-17]，美国环境保护局最近提出将氨纳入空气质量标准[18]。对于发展中国家来说，近五年来发表的关于氨排放的研究越来越多，尤其是在中国[8，19-20]。尽管如此，对于南美洲而言，关于氨排放量的唯一信息是全球数据库中报告的信息[21]。最近关于拉丁美洲和加勒比地区短期气候污染物的研究报告( 只有摘要可供利用) ，采用 GAINS 模型以及国家一级的信息估算了氨排放量[22]。

　　氨减排问题随着我国环境问题的凸显，显得愈来愈重要，如何客观、科学定量的评估我国区域氨排放量的问题尤为重要。因此，本研究在查阅国内外相关研究的基础上，从国内氨排放的估算方法发展历程、估算必需的数据类型与来源、参数取值以及不确定性产生等方面，比较了具有代表性的氨排放研究方法与结果，并与国外研究进行对比，提出了使区域氨排放估算进一步量化、精准的改进建议，以期为我国做好氨排放控制基础研究，制定相关管理政策等提供科学依据。

　　1、我国氨排放估算研究进展

　　我国氨排放估算领域起步于 20 世纪 90 年代。表 1 对比了我国各个地区的氨排放研究结果。早期相关的氨排放研究，大多从氨排放角度建立排放清单，排放清单以氮肥施用、畜禽养殖等农业氨排放源为主[3.7]。早期由于国内各类源排放因子数据的缺乏，研究者一方面参考欧洲地区相关研究，选用其中具有代表性的排放因子计算氨排放。例如，王文兴等[3]计算得到 1991 年全国氨的排放总量为 8918 Gg，其中畜禽、氨肥施用、人粪便与氮肥生产的排氨量占比分别为 64%、18%、17%和 1%，全国平均氨排放强度为 9 kg /hm2 。徐新华[23]采用类似方法，使用国外排放因子进行计算了江浙沪地区人为氨排放量。由于这些研究使用的排放因子来自于国外，研究结果的准确性可能会有偏差。另一方面，通过实地试验得到观测结果对氨的排放进行研究。朱兆良等[26]和蔡贵信等[27]在江苏丹阳、河南封丘等地运用了15 N 示踪技术和微气象学的方法对氨的挥发结果进行测定，研究发现石灰性稻田和酸性稻田土壤中氨挥发的情况差异明显，在酸性稻田区域，碳氨和尿素的氨挥发率分别为 19.5%和 8.8%; 而在石灰性稻田土壤中，碳氨和尿素的氨挥发率分别达到 39%和 30%。孙庆瑞等[24]在估算氮肥施用氨排放时，选取朱兆良等[26]和蔡贵信等[27]的氨排放观测数据作为我国氮肥施用氨排放因子，并将计算结果与欧洲地区氨排放量进行比较，结果显示，同期中国氨排放量大于全欧洲的排放量。以上为使用我国实地试验结果数据作为相关源氨排放因子的早期研究实例。Xing 和 Zhu [28]基于运用微气象学方法得到的研究结果，计算了 1990 年我国的农田氨排放量，并根据统计数据计算了不同氮肥( 尿素和碳酸氢铵) 分别在不同耕作方式下的氨挥发率。由于不同地区的气候条件，地理环境等影响农业发展的因素差异较大，加上已报道的排放因子数据非常有限，所以依靠单一的排放因子得到的计算结果存在较大的不确定性。因此，在后来的研究中，越来越多的研究者在考虑区域差异性和排放因子本地化等方面进行了探索。

　　Zheng 等[29]首次将模型法运用于我国的氨排放估算中。在对亚洲地区氮循环研究中，Zheng 等通过建立区域氮循环模型 IAP-N-1.0 模型，分析了 1961—2030 年亚洲各国家、地区的氮收支情况。在计算国内氨排放时，施肥农田部分采用已报道的国内旱地与水田施肥农田的氨排放因子，畜禽养殖排放因子选用 IPCC 推荐值，人体粪便部分则根据文献值计算。李富春等[7]在 IAP-N 模型基础上，综合考虑农田氮的输入量与氨的排放量。从粪便管理、施肥农田、秸秆燃烧等方面计算了川渝地区氨的排放量，并将计算结果分成 3 个时间段，分析该地区时间与空间的氨排放分布规律。张美双等[30]采用 NARSES 模型，对 2001 年我国种植业氮肥施用氨排放量进行估算，得到我国氨排放强度时空分布。房效凤等[5]在排放因子法的基础上，引入模型法对排放因子进行修正，在计算畜禽养殖氨排放时，通过 RAINS 模型计算出畜禽的 NH3实际排放因子，结合 NARSES 模型对氮肥施用氨排放排放因子进行修正，计算出 2011 年上海市农业源氨排放清单。类似地，Huang 等[2]在结合本地实验结果和修正排放因子的基础上，编制得到 2006 年我国氨排放清单。排放因子通过考虑环境温度、土壤酸度等参数得到，这能够使排放因子更加符合区域的地理环境。Kang 等[31]则在 Huang 等[2]的基础上，参考其估算方法，计算了 1980—2012 年我国的氨排放清单。此外，Wang 等[4]在 2012 年至 2013 年期间建立了全国范围内稻田氨排放的监测网络，并使用标准化的测量方法连续 2 年测量氨排放。该网络包括东北，东南和长江流域等中国主要水稻种植区域。结果表明，排放的氨占施用氮素的比例达 17.7%，2013 年中国稻田的氨总排放量估计为 1.7 Tg N/a。总体来讲，我国农业氨排放估算方法可以分为 3 个发展阶段( 图 1) ， 2011 年至今，国内氨排放研究领域取得一定的进步，排放因子的本地化与模型法估算得到了推广，研究结果的不确定性由早期( 80 年代末—90 年代末) 仅有单一的定性评估发展为定量评估，现有的利用多因素校正排放因子的方法和模型估算法具有较高的准确性。与国外的氨排放估算研究对比，国内在估算过程中仍是更多地依赖原有的单一排放因子法以及由国外开发的模型估算法，而对于适用于我国实地农业发展状况的模型以及根据不同地区影响因素修正的排放因子还需要进一步完善。

　　2 主要的估算方法与特征分析

　　排放因子法即根据排放源的活动水平与排放因子相乘，估算出该排放源的氨排放量。单一的排放源活动水平数据可获得性较高与可使用的本地排放因子数据较少，是排放因子法使用普遍的主要原因。早期研究所使用的排放因子法，多采用单一的排放源活动水平，如根据统计年鉴中牛、羊等畜禽的年末存栏量乘以对应的排放因子，排放因子则直接使用国外数据或多个国外氨排放因子的平均值[3，24]，由此计算得到我国或部分地区的氨排放水平。这种计算方法产生结果的不确定性，主要来源于国外排放因子与我国各地区实际排放因子之间的差异。房效凤等[5]与张美双等[30]在各自研究氨排放过程，引入模型法对排放因子进行修正，使排放因子更加完善，一定程度上提高了计算结果的准确性。Zheng 等[32]与李富春等[7]的研究中，根据 IPA-N 模型，从氮素循环的整体角度，计算了由氮肥总消耗、生物固氮、大气氮沉降回田等多个部分组成的农田生态系统氨排放，进一步提高了氮素输入量的准确性，但由于其各自研究中使用的排放因子都比较单一，主要排放源的排放因子缺乏进一步校正，依然存在一定的不确定性。Wu 等[20]运用 NARSES 模型计算福建省农田生态系统氨排放，通过考虑土壤 pH、耕地方式以及温度等影响因素进行排放因子的修正，而畜禽养殖、人体氨排放估算用到的活动水平数据则来源于政府提供的统计年鉴数据。与其他模型相比，NARSES 模型能够结合不同地域的特点对排放因子进行多因素修正，并且由此得到的氨排放时空分布精度较高。表 2 汇总了我国多个地区农业氨排放量估算的研究实例，对比分析了其方法特征。同时，参考通过不同途径获取的排放因子数据的不确定性评估方法，其不确定度数值参考了 TRACE-P 清单的经验数值的上限[12，33]，我们确定了各研究中排放因子的不确定度( 表 2) 。

　　单一的排放因子法仍是目前我国应用较为普遍的计算方法。在前人的研究基础上，许多研究者考虑到不同区域存在的差异性，对排放因子进行参数化，使其与该地区的实际环境更加符合。考虑到氨的排放是氮素循环的一个环节，建立了适用于我国农业状况的模型进行综合计算。同时，在氨排放的实地观测方面也进行了探索，为我国不同地区的氨排放研究提供了数据支持。但由于我国国土面积广大，各地区农业发展水平差异明显，仍存在( 1) 来自于实地测得的数据仍不够充足; ( 2) 已有的研究结果缺乏进一步验证与评估; ( 3) 广泛的模型估算法，在应用于不同地区时往往没有做出进一步调整; ( 4) 排放源的活动水平数据受限于统计资料的缺乏和不统一等问题，这导致重点排放源的识别不够且在同程度上降低了计算结果的准确性。

　　3 数据收集及参数取值

　　3.1 活动水平数据

　　目前国内氨排放研究中的活动水平数据包括: 主要农作物的播种面积和产量，主要畜禽( 牛、猪、羊、鸡等) 的饲养量，行政区划面积，耕地面积，乡村人口数，氮肥消费量，化石燃料消耗量等。大多来源于全国、省级统计年鉴，各部门统计资料及相关数据中心等，如中国统计年鉴、中国农村/农业统计年鉴、省级统计年鉴、省级农牧业统计资料、统计资料汇编、中国科学院资源环境科学数据中心、中国农业科学研究院数据中心等。数据收集过程中存在行政区划变异，指标不统一，部分数据缺失等问题。相关国际组织的数据库也可作为补充，Zheng 等[29]在运用模型法估算我国氨排放量时，使用的农作物数据来源于联合国粮农组织( FAO) 数据库。

　　3.2 参数取值方法

　　排放因子的选取对模型法及单一排放因子法估算氨排放量的准确性有着重要影响。氮肥施用的平均氨排放因子与各类氮肥的排放因子及多种氮肥施用比例有着较大关联。蔡贵信等[27]探索了我国碳氨及尿素等氮肥在水稻田中的主要损失途径，发现酸性粘质水稻土上这两种氮肥的氨挥发率为 20%和 9%。受限于各类氮肥施用比例数据的缺失，部分研究者直接从文献中获取氮肥施用的平均氨排放因子并加入计算之中。孙庆瑞等[24]根据调查资料计算得到我国氮肥使用比例; 王文兴等[3]从文献中获取到我国氮肥的使用比例。后续的氨排放研究较多借鉴这两份文献中的氮肥施用比例，但两者数据年代较早，是否与我国现阶段使用情况相符还需要进一步证实。董艳强等[35]则将多种氮肥的生产比例作为长三角地区氮肥的使用份额。周静等[36]根据调研资料及统计数据得到苏州地区多种氮肥的使用情况。Zhang 等[37]在估算我国氮肥施用产生的氨排放时，基于 2005 年我国县级调查数据得到各类氮肥的使用比例，精确程度较高。Huang 等[2]选用已有的测量结果作为基准的排放因子，在此基础上进行修正，这提高了计算结果的准确性。杨志鹏[6]基于物质流方法，参考 RAINS 模型，基于清单建立必需的基础数据，计算了畜禽养殖的氨排放因子。房效凤等[5]也根据 RAINS 模型计算了畜禽养殖的排放因子。从排放因子的不确定度来看，通过合理的实地测试所得到的结果具有较高的准确性，同时为相关研究提供重要的参照和数据支持，而使用相关模型法，通过参数的校正确定排放因子也是一种较好的提高排放因子准确性的方法。Wu 等[20]利用 NARSES 模型，计算了福建省的农田氮肥施用氨排放因子。这一模型考虑了土壤 pH、耕地使用方式、施肥率、降雨以及温度等因素。与直接使用文献中的排放因子数据计算相比，采用结合本地特点的参数得到的结果更加合理。但目前各个地区的可获得的参数不够完善，如施肥率和施肥方式等校正因子在一些地区出现数据缺失的情况，而使用省级或者全国性的平均数据则缺乏代表性。已有文献报道相关氮肥的使用比例及排放因子见表 3。

　　畜禽养殖氨排放是农业氨排放的重要组成部分，本研究选取牛、山羊、绵羊、猪、家禽、兔及马等几项指标，对其进行排放因子数据收集( 表 4) 。在计算畜禽养殖氨排放时，不同省份的活动水平存在差异。而同一排放源排放量也表现出明显的差别。国内在估算畜禽养殖氨排放时，选用的排放因子较多为文献中全国范围性的数据，但不同省份的畜禽养殖种类以及养殖条件均会随当地环境而改变。另外，李富春等[7]对放牧与非放牧部分牛、山羊以及绵羊氨排放因子的分配系数做了探索。

　　我国农业氨排放估算方法研究进展相关论文期刊推荐：《生态学报》(半月刊)创刊于1981年，是由中国科学技术协会主管、中国生态学学会中国科学院生态环境研究中心主办的我国生态学及生态学各分支学科研究领域的综合性学术期刊。主要报道生态学领域最新的基础理论和原始创新性研究成果，促进国内外学术交流和学术争鸣，推动我国的生态学研究发展，培养造就生态科学研究人才，为知识创新服务，为实施国家可持续发展和科教兴国服务。