2015—2019年北京市農(nóng)業(yè)氧化亞氮排放變化與分析*

王安吉 ，馬文林 ※，吳建繁 ，賈小紅

（1.北京建筑大學(xué)環(huán)境與能源工程學(xué)院，北京 100044；2.北京應(yīng)對(duì)氣候變化研究和人才培養(yǎng)基地，北京 100044；3.北京低碳農(nóng)業(yè)協(xié)會(huì)，北京 100107；4.北京市耕地建設(shè)保護(hù)中心，北京 100101）

0 引言

《巴黎協(xié)定》提出到2100年末將全球平均氣溫升幅控制在不超過工業(yè)化前水平2°C，并努力控制在1.5°C之內(nèi)[1]。CO2、CH4和N2O是導(dǎo)致全球升溫起主要作用的前三位溫室氣體。其中，N2O能長(zhǎng)期存在于大氣中，需要150年才能被完全分解[2]。依據(jù)IPCC發(fā)布的第四次氣候變化評(píng)估報(bào)告，N2O的全球增溫潛勢(shì)為CO2的298倍[3]，這個(gè)數(shù)值為同年公布的甲烷增溫潛勢(shì)的11.92倍。在目前的情況下，大氣中N2O含量每增加1倍，全球氣溫將提高0.44℃[4]。2019年全球農(nóng)業(yè)總計(jì)排放溫室氣體59億tCO2e[5]，N2O的貢獻(xiàn)率為34.4%。同年，中國(guó)農(nóng)業(yè)的溫室氣體排放量占到全球農(nóng)業(yè)總排放的11.2%，其中N2O排放占農(nóng)業(yè)自身總排放的45.9%[5]，而中國(guó)農(nóng)田土壤N2O直接和間接排放分別占全國(guó)N2O排放總量的55.8%和15.1%[6]。由此可見，為了實(shí)現(xiàn)《巴黎協(xié)定》的承諾，正確對(duì)待農(nóng)業(yè)N2O排放問題具有十分重要的意義。

目前，有許多研究者針對(duì)我國(guó)農(nóng)業(yè)溫室氣體排放問題進(jìn)行了研究。韋良煥等[7]依據(jù)《省級(jí)溫室氣體清單編制指南（試行）》（以下簡(jiǎn)稱“省級(jí)清單”）研究了中國(guó)省級(jí)區(qū)域農(nóng)業(yè)N2O的排放情況；李陽(yáng)等[8]運(yùn)用STIPART模型預(yù)測(cè)了中國(guó)省域農(nóng)業(yè)源N2O排放量；李庚欣等[9]研究了湖北省農(nóng)業(yè)碳排放空間分布特征；徐聰[10]研究了華北平原長(zhǎng)期施用氮肥與秸稈還田的N2O排放特征；李長(zhǎng)卓[11]利用Meta分析研究了耕作方式與麥田N2O排放的關(guān)系；陳煒等[12]依據(jù)種植業(yè)物料投入情況分析了1997—2015年中國(guó)種植業(yè)碳排放時(shí)空分布規(guī)律的特征；白義鑫等[13]分析了貴陽(yáng)市農(nóng)業(yè)碳排放與農(nóng)業(yè)GDP的脫鉤關(guān)系；Htun等[13]研究了中國(guó)西北地區(qū)秸稈與氮素管理對(duì)雨養(yǎng)農(nóng)業(yè)N2O與CH4排放的影響[14]；Su等[15]研究了中國(guó)東北地區(qū)典型雨養(yǎng)玉米田N2O通量的時(shí)態(tài)模式；Yang等[16]研究了中國(guó)西北地區(qū)雨養(yǎng)麥田不同施肥方式下的N2O排放情況；趙曉強(qiáng)等[17]分析了2006—2015年山西省農(nóng)業(yè)溫室氣體排放情況；王心宇等[18]分析了財(cái)政投入與農(nóng)業(yè)減排的關(guān)系；馬文林等[19]研究了長(zhǎng)期施糞肥對(duì)青貯玉米田N2O排放的影響；萬(wàn)小南等[20]研究了秸稈還田對(duì)楊凌區(qū)冬小麥—夏玉米田N2O排放的關(guān)系，上述研究人員分別采用不同研究方法、從不同角度對(duì)農(nóng)業(yè)溫室氣體排放做了相關(guān)研究。對(duì)于北京區(qū)域而言，張哲瑜等[21]對(duì)北京市規(guī)模化生豬養(yǎng)殖企業(yè)溫室氣體排放進(jìn)行了計(jì)算；劉月仙等[22]分析了北京市畜禽養(yǎng)殖溫室氣體的時(shí)空變化。目前尚缺少針對(duì)北京地區(qū)農(nóng)業(yè)N2O排放的全面研究。文章以IPCC《國(guó)家溫室氣體清單指南》（2006）（以下簡(jiǎn)稱“國(guó)家清單”）、“省級(jí)清單”和北京市相關(guān)地方標(biāo)準(zhǔn)《種植農(nóng)產(chǎn)品溫室氣體排放核算指南》（DB11/T 1564-2018）與《溫室氣體排放核算指南 畜牧養(yǎng)殖企業(yè)》（DB11/T 1422-2017）為核算方法，利用2015—2019年北京市種植業(yè)與畜牧業(yè)相關(guān)數(shù)據(jù)，對(duì)2015—2019年北京市農(nóng)業(yè)N2O排放量進(jìn)行核算并對(duì)核算結(jié)果進(jìn)行相應(yīng)分析。

1 研究方法

1.1 核算范圍與核算對(duì)象選取

該研究以2015—2019年為評(píng)估期，核算目標(biāo)是評(píng)估期內(nèi)北京全市與各區(qū)種植業(yè)和畜牧業(yè)的N2O排放量。

1.1.1 核算范圍與核算對(duì)象

種植業(yè)N2O排放核算范圍包括農(nóng)田土壤施肥和秸稈處置利用導(dǎo)致的N2O排放直接排放、氮揮發(fā)間接排放和氮淋溶/徑流間接排放，畜牧業(yè)N2O排放放核算范圍為糞便管理導(dǎo)致的N2O排放。

1.1.2 核算對(duì)象選擇

施肥N2O排放核算對(duì)象包括化肥與有機(jī)肥。根據(jù)調(diào)查，北京市主要在小麥、玉米與蔬菜種植中使用有機(jī)肥。依據(jù)《北京統(tǒng)計(jì)年鑒2020》中列出的主要農(nóng)作物，將秸稈還田N2O排放核算對(duì)象確定為水稻、小麥、玉米、蔬菜、大豆、薯類、花生和棉花。

依據(jù)《北京統(tǒng)計(jì)年鑒2020》與各區(qū)統(tǒng)計(jì)年鑒中列出的主要畜禽種類，結(jié)合北京市畜牧養(yǎng)殖特點(diǎn)，將糞便管理N2O排放對(duì)象確定為豬、肉牛、奶牛、羊和家禽。

1.2 核算方法與數(shù)據(jù)來(lái)源

1.2.1 核算方法

該文參考“國(guó)家清單”“省級(jí)清單”、DB11/T 1564-2018和DB11/T 1422-2017進(jìn)行種植業(yè)和畜牧業(yè)N2O排放核算，秸稈還田N2O排放按式（1）計(jì)算，施肥N2O排放按式（2）至（5）計(jì)算，糞便管理N2O排放按式（6）計(jì)算。

式(1)中，E秸稈為作物秸稈還田產(chǎn)生的N2O排放，單位為萬(wàn)tCO2e；F秸稈為作物秸稈還田氮量，單位為kgN；EF1為農(nóng)田N2O直接排放因子，單位為kgN2O-N/kgN；Frac淋溶/徑流為氮淋溶/徑流損失系數(shù)；EF3為農(nóng)田氮投入淋溶/徑流N2O間接排放因子，單位為kgN2O-N/kgN；R為秸稈還田率；為N2O-N轉(zhuǎn)化為N2O的系數(shù)；298為N2O增溫潛勢(shì)；10-3為kg轉(zhuǎn)化為t的系數(shù)。

式（2）中，EFN為農(nóng)田施肥產(chǎn)生的N2O排放，單位為萬(wàn)tCO2e；EFND為農(nóng)田施肥N2O直接排放，單位為萬(wàn)tCO2e；EFNI-ATD為肥料因氮揮發(fā)產(chǎn)生的N2O排放，單位為萬(wàn)tCO2e；EFNI-L為肥料因氮淋溶/徑流產(chǎn)生的N2O排放，單位為萬(wàn)tCO2e。

式（3）中，EFND為農(nóng)田施肥產(chǎn)生的N2O直接排放，單位為萬(wàn)tCO2e；F化肥為農(nóng)田施入化肥含氮量，單位為kgN；F有機(jī)肥為農(nóng)田施入有機(jī)肥含氮量，單位為kgN；EF1為農(nóng)田N2O直接排放因子，單位為kgN2O-N/kgN；為N2O-N轉(zhuǎn)化為N2O的系數(shù)；298為N2O增溫潛勢(shì)；10-3為kg轉(zhuǎn)化為t的系數(shù)。

式（4）中，EFN-ATD為肥料因氮揮發(fā)產(chǎn)生的N2O排放，單位為萬(wàn)tCO2e；F化肥為農(nóng)田施入化肥含氮量，單位為kgN；FracGASF為化肥氮揮發(fā)系數(shù)；F有機(jī)肥為農(nóng)田施入有機(jī)肥含氮量，單位為kgN/年；FracGASM為有機(jī)物料氮揮發(fā)系數(shù)；EF2為農(nóng)田施肥氮揮發(fā)N2O間接排放因子，單位為kgN2O-N/kgN；為N2O-N轉(zhuǎn)化為N2O的系數(shù)；298為N2O增溫潛勢(shì)；10-3為kg轉(zhuǎn)化為t的系數(shù)。

式（5）中，EFN-L為肥料因氮淋溶/徑流產(chǎn)生的N2O排放，單位為萬(wàn)tCO2e；F化肥為農(nóng)田施入化肥含氮量，單位kgN；F有機(jī)肥為農(nóng)田施入有機(jī)肥含氮量，單位kgN；FracLEACH-(H)為氮淋溶徑流損失系數(shù)；EF3為農(nóng)田氮投入淋溶/徑流N2O間接排放因子，單位為kgN2O-N/kgN；為N2O-N轉(zhuǎn)化為N2O的系數(shù)；298為N2O增溫潛勢(shì)；10-3為kg轉(zhuǎn)化為t的系數(shù)。

式（6）中，EMN為動(dòng)物糞便管理N2O排放，單位為萬(wàn)tCO2e；AP為動(dòng)物年末存欄量，單位為頭；EFMN為動(dòng)物的糞便管理N2O排放因子，單位為kgN2O/(頭·年)；298為N2O增溫潛勢(shì)；10-3為kg轉(zhuǎn)化為t的系數(shù)。

1.2.2 數(shù)據(jù)來(lái)源

（1）活動(dòng)數(shù)據(jù)：活動(dòng)數(shù)據(jù)來(lái)源包括《北京統(tǒng)計(jì)年鑒》（2016—2020）、《北京區(qū)域統(tǒng)計(jì)年鑒》（2016—2020）、北京市各區(qū)的地方統(tǒng)計(jì)年鑒以及北京市農(nóng)業(yè)相關(guān)管理部門的調(diào)研數(shù)據(jù)，如表1所示。將活動(dòng)數(shù)據(jù)根據(jù)來(lái)源分為3個(gè)優(yōu)先級(jí)，第一級(jí)為各區(qū)統(tǒng)計(jì)年鑒，第二級(jí)為《北京統(tǒng)計(jì)年鑒》和《北京區(qū)域統(tǒng)計(jì)年鑒》，第三級(jí)為北京市農(nóng)業(yè)相關(guān)管理部門的調(diào)研數(shù)據(jù)。由于東城、西城和石景山區(qū)不進(jìn)行農(nóng)業(yè)生產(chǎn)，因此核算范圍不包括此3區(qū)。

表1 數(shù)據(jù)來(lái)源

（2）排放因子：該文所用各項(xiàng)排放因子及相關(guān)參數(shù)如表2至表4所示。

表2 主要農(nóng)作物參數(shù)[23]

表4 動(dòng)物糞便管理N2O排放因子[25]

表3 土壤氮投入排放因子[24]

2 結(jié)果與分析

2.1 種植業(yè)N2O排放

2.1.1 秸稈還田N2O排放

秸稈還田率采用秸稈肥料化利用情況的調(diào)研數(shù)據(jù)。評(píng)估期北京市農(nóng)作物秸稈還田N2O排放情況如圖1所示。從圖1可以看出，2015—2019年全市秸稈還田N2O排放量從1.30萬(wàn)tCO2e下降至0.77萬(wàn)tCO2e，降幅高達(dá)40.65%，年均減排率為12.27%。2019年秸稈還田N2O直接排放量為0.55萬(wàn)tCO2e，約為其間接排放量（0.22萬(wàn)tCO2e/年）的2.5倍。

圖1 2015—2019年北京市秸稈還田N2O排放

2.1.2 施肥N2O排放

（1）化肥N2O排放：評(píng)估期北京市農(nóng)作物種植施化肥的N2O排放情況如圖2所示。其中，N2O間接排放包含氨揮發(fā)和氮淋溶滲濾兩種途徑。從圖2可以看出，2015—2019年北京市化肥N2O排放量從20.34萬(wàn)tCO2e下降至9.58萬(wàn)tCO2e，呈穩(wěn)定下降趨勢(shì)，降幅高達(dá)52.91%，年均減排率為17.16%。2019年化肥N2O直接排放量為6.10萬(wàn)tCO2e，約為其間接排放量（3.48萬(wàn)tCO2e）的1.75倍。

圖2 2015—2019年北京市化肥N2O排放

（2）有機(jī)肥N2O排放：有機(jī)肥含氮量統(tǒng)一以1%計(jì)。評(píng)估期北京市農(nóng)作物種植施有機(jī)肥的N2O排放情況如圖3所示。從圖3可以看出，2015—2019年北京市有機(jī)肥N2O排放量從19.52萬(wàn)tCO2e下降至5.28萬(wàn)tCO2e，總體呈下降趨勢(shì)，降幅高達(dá)73%，年均減排率為27.88%。2019年有機(jī)肥N2O直接排放量為3.02萬(wàn)tCO2e，約為其間接排放量（2.26萬(wàn)tCO2e）的1.34倍。

圖3 2015—2019年北京市有機(jī)肥N2O排放

2.1.3 種植業(yè)N2O總排放

（1）全市種植業(yè)N2O排放：評(píng)估期北京市種植業(yè)N2O總排放情況如圖4所示。由圖4可知，2015—2019年北京種植業(yè)N2O排放量從41.16萬(wàn)tCO2e下降至15.63萬(wàn)tCO2e，總體呈下降趨勢(shì)，排放降幅高達(dá)62.0%。施肥是種植業(yè)N2O排放的重要排放源（年均排放占比95.6%）。在各類種植業(yè)N2O排放源中，化肥在種植業(yè)N2O排放中占比最大（年均53.8%），其次為有機(jī)肥（年均42.7%），秸稈還田（年均3.5%）最小。2019年這三者產(chǎn)生的N2O排放量較2015年分別減少了52.9%、73.0%和40.7%。

圖4 2015—2019年北京市種植業(yè)N2O排放

（2）各區(qū)種植業(yè)N2O排放：2019年北京市各區(qū)種植業(yè)N2O排放情況如圖5所示，相對(duì)于2015年的減排率見表5。從圖5可以看出，2019年大興區(qū)種植業(yè)N2O排放量（3.65萬(wàn)tCO2e）居全市首位，占全市種植業(yè)N2O排放量的23.4%。根據(jù)表5，通州區(qū)2019年較2015年種植業(yè)N2O減排率最高（78.1%），年均減排率為31.6%。

表5 2019年北京市各區(qū)種植業(yè)N2O減排率（以2015年為基準(zhǔn)年）

圖5 2019年北京市各區(qū)種植業(yè)N2O排放

2.2 畜牧業(yè)N2O排放

2.2.1 全市畜牧業(yè)N2O排放情況

評(píng)估期北京市畜牧業(yè)N2O總排放情況和各畜種N2O排放比例分別如圖6和圖7所示。

從圖6可以看出，2015—2019年北京市畜牧業(yè)N2O排放量從32.04萬(wàn)tCO2e下降至10.59萬(wàn)tCO2e，總體呈遞減趨勢(shì)，總減排率為66.9%，年均減排率24.2%；在各畜種N2O排放中，豬N2O排放的減排量最高，減排8.17萬(wàn)tCO2e，降幅92.0%。

圖6 2015—2019年北京市畜牧業(yè)N2O排放

圖7展示了評(píng)估期各類畜禽種N2O排放占比情況，可以得出評(píng)估期家禽N2O排放都居全市畜牧業(yè)各畜種N2O排放首位，年均占比42.8%。

圖7 2015—2019年北京市畜牧業(yè)各畜種N2O排放比例

2.2.2 各區(qū)畜牧業(yè)N2O排放情況

2019年各區(qū)畜牧業(yè)N2O排放情況及其較2015年減排情況分別見圖8和表6。

從圖8看，豐臺(tái)區(qū)的排放量降為0，因其將畜禽養(yǎng)殖全部取締；平谷區(qū)畜牧業(yè)N2O排放量為3.29萬(wàn)tCO2e，居全市首位，對(duì)全市畜牧業(yè)N2O貢獻(xiàn)率達(dá)31.1%。其余各區(qū)對(duì)全市畜牧業(yè)N2O排放貢獻(xiàn)率由大到小排序依次為順義區(qū)、密云區(qū)、延慶區(qū)、房山區(qū)、通州區(qū)、大興區(qū)、懷柔區(qū)、昌平區(qū)、海淀區(qū)、朝陽(yáng)區(qū)和門頭溝區(qū)。

圖8 2019年北京市各區(qū)畜牧業(yè)N2O排放

由表6可知，豐臺(tái)區(qū)畜牧業(yè)N2O排放減排率最高，為100%。其余各區(qū)按畜牧業(yè)N2O減排率由高到低排序依次為懷柔區(qū)、門頭溝區(qū)、海淀區(qū)、大興區(qū)、通州區(qū)、昌平區(qū)、朝陽(yáng)區(qū)、延慶區(qū)、順義區(qū)、延慶區(qū)、密云區(qū)和平谷區(qū)。

表6 2019年北京市各區(qū)畜牧業(yè)N2O減排率（以2015年為基準(zhǔn)年）

2.3 農(nóng)業(yè)N2O排放

2.3.1 全市農(nóng)業(yè)N2O排放情況

評(píng)估期北京市農(nóng)業(yè)N2O排放情況如圖9所示。由圖9可看出，2015—2019年北京市農(nóng)業(yè)N2O排放總體呈下降趨勢(shì)，從73.2萬(wàn)tCO2e下降至26.22萬(wàn)tCO2e，減少了46.98萬(wàn)tCO2e，減排率高達(dá)64.2%，年均減排率為22.6%。

圖9 2015—2019年北京市農(nóng)業(yè)N2O排放

2.3.2 各區(qū)農(nóng)業(yè)N2O排放情況

2019年各區(qū)農(nóng)業(yè)N2O排放如圖10所示，較2015年減排情況如表7所示。

表7 2019年較2015年北京市各區(qū)農(nóng)業(yè)N2O減排率

從圖10可知，順義區(qū)農(nóng)業(yè)N2O排放量為4.85萬(wàn)tCO2e，居各區(qū)農(nóng)業(yè)N2O排放之首，對(duì)全市農(nóng)業(yè)N2O貢獻(xiàn)率為18.5%；豐臺(tái)區(qū)的排放最低，為0.03萬(wàn)tCO2e；懷柔區(qū)的排放量較2015年降幅最為顯著，減排率達(dá)82.7%，年均減排率為35.5%。

圖10 2019年北京市各區(qū)農(nóng)業(yè)N2O排放

3 討論與結(jié)論

3.1 討論

（1）不同區(qū)存在較大排放差異：北京市各區(qū)農(nóng)業(yè)N2O排放差異較為顯著，其主要原因?yàn)楦鲄^(qū)肥料施用量與畜禽養(yǎng)殖量存在顯著差異。以2019年豐臺(tái)區(qū)與順義區(qū)為例，豐臺(tái)區(qū)化肥施用折純量為74.8t、無(wú)畜禽養(yǎng)殖，順義區(qū)化肥施用折純量為13 390.5t、各類畜禽76.8萬(wàn)只。這兩個(gè)區(qū)由于化肥施用與畜禽養(yǎng)殖規(guī)模的不同使2019年農(nóng)業(yè)N2O產(chǎn)生了3.9萬(wàn)tCO2e的排放差異。

（2）種植規(guī)模和氮投入水平不斷降低：依據(jù)《北京統(tǒng)計(jì)年鑒》數(shù)據(jù)，2015年北京市農(nóng)作物播種面積為17.7萬(wàn)hm2，2019年減少到9.2萬(wàn)hm2，減少率為48.02%；同期全市農(nóng)田化肥施用折純氮量從4.852 4萬(wàn)t減少到2.248 13萬(wàn)t，下降率為53.67%。對(duì)二者進(jìn)行比較，氮投入量下降率大于播種面積的減少率，說(shuō)明北京種植業(yè)N2O排放減少的主要原因有兩方面，最主要的是播種面積的下降，同時(shí)單位面積氮投入水平的下降，也具有一定的影響。2016年北京市發(fā)布了《北京市“十三五”時(shí)期都市現(xiàn)代農(nóng)業(yè)發(fā)展規(guī)劃》，規(guī)劃中明確要減少化肥使用量、逐步退出高耗水糧食作物種植，這對(duì)北京市化肥氮投入量的減少產(chǎn)生了積極影響。

在北京市有種植業(yè)生產(chǎn)的13個(gè)區(qū)中，只有門頭溝區(qū)2019年的種植業(yè)N2O排放比2015年增加。依據(jù)《北京區(qū)域統(tǒng)計(jì)年鑒》，2015年門頭溝區(qū)施用化肥折純量為148.4t，2019年增加到931.7t，使得該區(qū)2019年化肥N2O排放較2015年增加了0.11萬(wàn)tCO2e。而調(diào)查數(shù)據(jù)表明該區(qū)同期有機(jī)肥投入量和秸稈還田量變化都不大。因此，2019年門頭溝區(qū)化肥施用量增加是該區(qū)種植業(yè)N2O排放增加的根本原因。對(duì)此，可以通過有機(jī)肥替代化肥、施用緩釋肥與長(zhǎng)效肥等方法減少N2O的排放[19,26]，此外有研究表明施用生物炭可提高土壤健康水平并降低因細(xì)菌途徑產(chǎn)生的N2O[27]。

（3）畜牧業(yè)N2O排放的快速減少：評(píng)估期北京市畜牧業(yè)N2O排放年均減排率高于種植業(yè)，這與禁養(yǎng)政策和畜禽疫病捕殺有著緊密的關(guān)聯(lián)。2016年環(huán)境保護(hù)部與原農(nóng)業(yè)部共同發(fā)布了《畜禽養(yǎng)殖禁養(yǎng)區(qū)劃定技術(shù)指南》，2017年北京市出臺(tái)了《北京市畜禽養(yǎng)殖禁養(yǎng)區(qū)劃定工作方案》《關(guān)于盡快開展禁養(yǎng)區(qū)內(nèi)規(guī)模化畜禽養(yǎng)殖場(chǎng)（小區(qū)）關(guān)閉或搬遷工作的函》等系列文件，明確了禁養(yǎng)區(qū)劃定原則和范圍以及關(guān)閉或搬遷工作的具體措施。隨后各區(qū)先后劃定了畜禽禁養(yǎng)區(qū)，北京市畜禽快速縮小。

畜禽養(yǎng)殖存欄量的減小是畜牧業(yè)N2O減排的直接因素。依據(jù)《北京統(tǒng)計(jì)年鑒》數(shù)據(jù)，2016—2019年全市豬存欄量分別為上一年的99.8%、67.9%、40.5%和29%。2018年受非洲豬瘟影響，大量生豬遭到捕殺，豬年底存欄量為45.4萬(wàn)頭，較2017年減少了69萬(wàn)頭，直到2019年生豬養(yǎng)殖規(guī)模依然在減小。由此得出，在不減少畜禽產(chǎn)品生產(chǎn)量前提下，如果可以提高畜禽生產(chǎn)效率，就可以減少養(yǎng)殖量，起到減排的效果。因此，培育高產(chǎn)品種、提高飼養(yǎng)管理水平、降低畜禽養(yǎng)殖病死率、增加繁殖率等措施，對(duì)減少畜禽養(yǎng)殖N2O排放具有顯著作用。

（4）不確定性分析：受核算對(duì)象、方法與數(shù)據(jù)來(lái)源的影響，該文核算結(jié)果存在一定的不確定性。此次核算采用的所有排放因子均為現(xiàn)行標(biāo)準(zhǔn)中所列的缺省值，因而核算結(jié)果較真實(shí)排放量有一定差異。文中所用活動(dòng)數(shù)據(jù)為公開發(fā)表的年鑒數(shù)據(jù)或從北京農(nóng)業(yè)管理部門調(diào)研所得，可信度較高。但因部分活動(dòng)數(shù)據(jù)和排放因子數(shù)值與其他研究者所用參數(shù)存在差異，使得核算結(jié)果不相同，甚至于有較大差距。

例如，2019年的一項(xiàng)研究得出2016年北京市畜禽糞便N2O排放量為646 t N2O[7]，折合19.25萬(wàn)t CO2e，該文核算得出同年北京畜禽N2O排放量為29.27萬(wàn)t CO2e，造成這種差異的最主要原因是由于選用了不同的排放因子。以奶牛糞便管理N2O排放因子為例，該研究選取DB11/T 1422-2017提供的缺省值1.94 kg/(頭·年)，而該研究取“省級(jí)清單”中奶牛、非奶牛和水牛糞便管理N2O排放因子之平均值，經(jīng)計(jì)算，該平均值為0.88 kg/(頭·年)。

此外，數(shù)據(jù)來(lái)源對(duì)于核算結(jié)果的影響同樣不可忽視。該文所用數(shù)據(jù)為北京市及各市轄區(qū)統(tǒng)計(jì)年鑒數(shù)據(jù)與北京市農(nóng)業(yè)管理部門提供的調(diào)研數(shù)據(jù)，并且在數(shù)據(jù)收集方面考慮了有機(jī)肥使用情況與秸稈利用情況。而2019年的研究[7]不僅針對(duì)北京市，還包括中國(guó)其他省份，因而該研究以《中國(guó)統(tǒng)計(jì)年鑒》作為數(shù)據(jù)來(lái)源，且氮肥折純比例采用了寧夏的數(shù)據(jù)[28]。這些數(shù)據(jù)來(lái)源的不同，使該文對(duì)農(nóng)用地N2O排放的核算結(jié)果是已有的研究[7]結(jié)果的5倍。

3.2 結(jié)論

該研究對(duì)2015—2019年北京市農(nóng)業(yè)源N2O排放進(jìn)行了核算，并對(duì)核算結(jié)果進(jìn)行了分析，主要結(jié)論如下。

（1）2015—2019年北京市農(nóng)業(yè)N2O排放呈下降趨勢(shì)。種植業(yè)為北京市農(nóng)業(yè)N2O最大排放源，其對(duì)北京市2015—2019年農(nóng)業(yè)N2O年均貢獻(xiàn)率為58.3%。

（2）2015—2019年北京市種植業(yè)和畜牧業(yè)的N2O排放均呈下降趨勢(shì)，它們的最大排放源分別是化肥與家禽，二者對(duì)種植業(yè)和畜牧業(yè)的N2O排放年均貢獻(xiàn)率分別為54.1%和42.8%。

（3）2019年全市畜牧業(yè)N2O與種植業(yè)N2O排放量均較2015年有所下降。門頭溝區(qū)是2019年較2015年全市唯一種植業(yè)N2O排放有所增加的地區(qū)，原因是化肥施用量的增加。