Role of Land Use Changes on Ammonia Emissions from Agricultural Ecosystems in the Yangtze River Delta Region from 2000 to 2018
NH3作为大气中最重要的碱性痕量气体, 是全球氮循环的关键成分[1], 可以中和二氧化硫和氮氧化物的氧化产物生成二次无机盐(SNA), 增加大气中细颗粒物(PM2.5)浓度.已有研究表明, 我国SNA的质量占总PM2.5质量可高达60%[2~5], 进而降低城市能见度[6], 在霾污染形成过程中具有十分重要的作用[7]; 在严重的霾天气中, PM2.5浓度受NH3排放的影响更大[8], 严重威胁人类健康[9, 10].NH3排放源包括农业、生物质燃烧、人体排泄物、化工生产、废弃物处置和机动车尾气等[11~13]; 其中, 农业源是主要的NH3排放源, 其中又以氮肥施用和畜禽养殖为主.近年来有研究表明, 氮肥施用排放的NH3占全球NH3排放总量的40%, 畜禽养殖排放的NH3占50%[14].
由于农业集约化生产的盛行, 20世纪NH3排放迅速增高[15~17].在过去的几十年间, 中国NH3的排放总量整体上处于增高趋势, 从1980年的4.7 Mt·a-1(以氮计)上升到2016年的11 Mt·a-1, 增长了2.4倍[18]; 从2000年的12.1 Mt·a-1增长到2015年的15.6 Mt·a-1, 年增长率达1.9%[19].国内学者对NH3排放清单的研究结果表明, 畜禽养殖和氮肥施用产生的NH3占我国总排放量的80%以上[12, 20~22].作为我国东部三大典型城市群之一, 长江三角洲(Yangtze River Delta, YRD)地区的PM2.5浓度受到农业源NH3排放的影响也较为显著. 2000~2017年江苏省农业源NH3排放平均78%来源于畜禽养殖, 22%来源于氮肥施用, 2012年NH3排放量出现最大值(837.64 kt·a-1)[7]. 2013年浙江省农业源的NH3排放主要来自氮肥施用(43%)[23], 2017年农田生态系统NH3排放量达到36.06 kt, 并以氮肥施用贡献最大(87%)[24]. 2006年安徽省NH3排放量为320.2 kt, 其中52.7%来自农田生态系统[12]. 2006~2014年, 氮肥施用是YRD NH3排放的最重要来源, 该来源排放量约占总NH3排放量的55%[25].城市化水平与NH3排放具有负相关关系, 因此, 土地利用变化对NH3排放的研究具有十分重要的意义[26].
过去20年间, 学者们对农业NH3排放的研究大多是对其排放变化趋势和时空分布特征的研究.随着快速城市化和机动化进程, 在过去十几年间, YRD城市群土地利用状况发生了巨大变化, 直接影响到农业NH3的排放.土地利用变化直接导致土壤本底NH3排放的变化, 间接影响氮肥施用总量进而导致氮肥利用率的变化.然而, 针对土地利用变化对NH3排放影响的研究还相对匮乏, 利用高分辨率卫星遥感资料反演我国典型城市群下垫面土地利用演变, 并进而研究其对农业源NH3排放的影响, 鲜见报道.为了弄清楚YRD地区历史土地利用变化对农田生态系统NH3排放的影响, 本研究利用Landsat卫星遥感手段解译高分辨率2000~2018年历史土地利用状况, 采用自下而上的方法, 结合情景分析, 探讨了YRD在过去20年间由于历史土地利用变化对农田生态系统NH3排放变化的影响, 以期为YRD区域PM2.5防控提供科学依据.
1 材料与方法
1.1 卫星遥感影像解译土地利用变化
本研究涵盖了上海及江苏、浙江和安徽三省的长江三角洲土地利用范围, 采用基于2000、2005、2010、2015和2018年Landsat系列卫星遥感影像解译的YRD土地利用数据, 该数据采用全国生态10a土地覆被Ⅰ、Ⅱ级分类系统, 将YRD土地覆被类型分为林地、草地、水域、耕地、建设用地及未利用地6大类.利用ArcGIS 10.6软件, 创建1 km×1 km格网, 提取各格网中的土地利用信息.
1.2 NH3排放计算方法
本文主要研究YRD农田生态系统NH3的排放, 主要涉及源类包括氮肥施用、土壤本底、固氮植物和秸秆堆肥, 本文采用排放因子法计算NH3排放量, 计算公式如下:
|
(1) |
式中, i和j分别为地区和源类别, A为活动水平数据, 主要包括氮肥施用量、耕地面积、固氮植物的种植面积和田间堆肥的秸秆量; EF(NH3)为NH3的排放因子.
1.2.1 活动水平
氮肥施用活动水平:为YRD广泛使用的氮肥种类, 包括尿素、碳铵、硝铵和硫铵等, 其中以尿素和碳铵施用最为广泛.氮肥施用总量来自2000~2018年各省市统计年鉴[27~31], 各种氮肥施用占比参考全国各种氮肥生产比例和董艳强等[32]的研究, 其中尿素占64%; 碳铵占24%; 硝酸铵、硫酸铵及其他氮肥各占4%.
土壤本底活动水平:为该地区的耕地面积, 数据来自2000、2005、2010、2015和2018年1∶50 000彩红外航空遥感图像解译的土地利用类型数据.
固氮植物活动水平:为该地区2000~2018年固氮植物(豆类、花生和绿肥)的种植面积, 数据来自2000~2018年各省市统计年鉴及农业农村局公布的数据.
秸秆堆肥活动水平:为该地区2000~2018年主要可用于堆肥的农作物(水稻、小麦、玉米、棉花、豆类、花生和油菜)的田间堆肥的秸秆量.利用该地区作物产量、谷草比、秸秆堆肥比例三者乘积来估算田间堆肥量.秸秆堆肥活动水平的计算公式如下:
|
(2) |
式中, j为作物类型, A为作物堆肥秸秆量; H为作物产量; K为秸秆与产品的比例; D为秸秆干物质比; C为秸秆堆肥率.其中, 作物产量来自文献[27~31], R和D值来自Zhou等[33]的研究.
1.2.2 排放因子
氮肥施用源EF(NH3)计算公式如下:
|
(3) |
式中, i为地区, EF为基准排放因子, P为施肥率校正系数, R为施肥方式校正系数.本文在计算YRD氮肥施用时依据YRD农时情况主要考虑了水稻、小麦、玉米、豆类和油菜种植过程所施基肥和追肥.上述5种作物的施肥时间主要发生在4~10月, 见表 1, 该时间段内YRD大部分区域的平均气温为20~30℃.
表 1
(Table 1)
表 1 长三角主要施肥作物及施肥时间
Table 1 Major crop types and fertilization time in the YRD
作物 |
水稻 |
小麦 |
玉米 |
豆类 |
油菜 |
基肥时间 |
移栽期(5~7月) |
播种期(10月) |
播种期(4月) |
播种期(3~4月) |
播种期(9~10月) |
追肥时间 |
分蘖期(5~8月) |
拔节期(4~5月) |
拔节期(5~6月) |
— |
抽薹期(4月) |
|
表 1 长三角主要施肥作物及施肥时间
Table 1 Major crop types and fertilization time in the YRD
|
根据文献[34], 由于土壤酸碱性的差异, 尿素、碳铵、硝酸铵、硫酸铵和其他氮肥的基准排放因子不同, YRD各亚区域土壤酸碱性见表 2; 每亩耕地施氮肥高于13 kg的地区, P取1.18, 其它地区取1.0; 在覆土深施时R取0.32, 表面撒施时R为1.0, 本文R取平均值0.66.
表 2
(Table 2)
表 2 长三角各地区土壤酸碱性
Table 2 Soil acidity-alkalinity in various regions of the YRD
地区 |
上海 |
苏北 |
苏南 |
浙北 |
浙南 |
皖北 |
皖南 |
酸碱性 |
中性 |
碱性 |
中性 |
碱性 |
碱性 |
中性 |
酸性 |
|
表 2 长三角各地区土壤酸碱性
Table 2 Soil acidity-alkalinity in various regions of the YRD
|
土壤本底源EF(NH3):根据文献[34], 土壤本底源EF(NH3)(以氨计)取0.18 g·(m2·a)-1.
固氮植物源EF(NH3):根据文献[34], 固氮植物(大豆、花生和绿肥)源EF(NH3)(以氨计)取值分别为0.10、0.12和0.13 g·(m2·a)-1.
秸秆堆肥源EF(NH3):根据文献[34], 秸秆堆肥源EF(NH3)(以氨计)取0.32 kg·t-1.
1.3 情景设置
农田生态系统NH3排放的历史变化, 一方面受到土地利用变化的影响, 另一方面还受到氮肥施用变化的影响.随着城市化进程的加快, 土地利用发生了变化; 同时随着科技进步、农业生产方式的改变, 氮肥施用量和施用方式也产生了变化.本研究依据2000~2018年YRD地区土地利用变化和氮肥施用变化设置了9种情景(见表 3), 比较不同情景下NH3的排放情况, 以探讨土地利用变化影响、氮肥施用变化影响以及两者共同的影响.情景1(S1)选用2000年土地利用和2000年氮肥施用率计算2000年农田生态系统NH3排放量; 情景2(S2)选用2005年土地利用和2000年氮肥施用率; 情景3(S3)选用2005年土地利用和2005年氮肥施用率计算2005年农田生态系统NH3排放量.分别比较情景1与情景2(S2-S1)、情景2与情景3(S3-S2)和情景1与情景3(S3-S1), 得出2000~2005年土地利用变化、氮肥施用变化和两者共同变化对农田生态系统NH3排放的影响.以相同的方式分别获得2000~2010、2000~2015和2000~2018年农田生态系统NH3排放的影响.
表 3
(Table 3)
表 3 氨排放计算情景设置
Table 3 Scenarios of ammonia emission calculation
情景 |
情景描述 |
年份 |
土地利用 |
氮肥施用 |
S1 |
2000年土地利用+2000年氮肥施用 |
2000 |
2000 |
S2 |
2005年土地利用+2000年氮肥施用 |
2005 |
2000 |
S3 |
2005年土地利用+2005年氮肥施用 |
2005 |
2005 |
S4 |
2010年土地利用+2000年氮肥施用 |
2010 |
2000 |
S5 |
2010年土地利用+2010年氮肥施用 |
2010 |
2010 |
S6 |
2015年土地利用+2000年氮肥施用 |
2015 |
2000 |
S7 |
2015年土地利用+2015年氮肥施用 |
2015 |
2015 |
S8 |
2018年土地利用+2000年氮肥施用 |
2018 |
2000 |
S9 |
2018年土地利用+2018年氮肥施用 |
2018 |
2018 |
|
表 3 氨排放计算情景设置
Table 3 Scenarios of ammonia emission calculation
|
2 结果与讨论
2.1 2000~2018年YRD土地利用变化
图 1所示为基于Landsat卫星遥感解译的2000~2018年1 km×1 km分辨率土地利用变化情况.研究区域面积共计35.8×104 km2, 从总体空间分布来看, 耕地主要分布在江苏、安徽平原及浙江北部; 林地主要分布在安徽西南部以及浙江丘陵一带; 建设用地主要呈现以上海为中心, 长江入海口和沿岸一带为轴带, 江苏的苏锡常和浙江的杭州湾城市连片发展的趋势; 水域主要分布在长江、太湖和洪泽湖等大面积湖泊区域.
由图 2可知, 2000~2018年YRD土地利用变化特征突出表现为, 城市建设用地迅速扩张和耕地面积锐减.城市建设用地自2000年到2010年增幅明显, 2000年建设用地面积为36 326.14 km2, 占比10.15%, 至2010年建设用地面积达53 044.62 km2, 占比14.81%, 2000~2018年建设用地共增长了77.86%.耕地类型是城镇建设用地主要的转出类型, 2000~2010年为建设用地快速增加的时期, 江苏省、上海市以及杭州湾等城市群连片发展, 迅速扩张, 占用了大量的耕地区域, 导致耕地面积大幅下降.耕地面积在2000年到2018年持续减少, 尤其在2000年至2010年降幅明显, 研究时段内耕地面积共减少了10%.林地面积从2005年开始保持小幅减少的趋势, 从2005年占比30.19%下降至2018年的30.06%.水域面积呈现先减少后增加的趋势, 2010~2018年增加了1 610.81 km2.草地和未利用地则变化较小.
由图 3、图 4和表 4可知, 随着YRD城市化进程的推进, 41个城市18年来耕地面积均呈减少的趋势.其中, 安徽省整体减幅最小, 上海市最为明显, 耕地面积减少了1 438 km2(34.06%), 这是由于2000~2018年上海市快速城市化所致.苏锡常地区经济发展迅速, 城市化比率大幅提高, 耕地面积也大幅减少, 分别减少了33%、28%和21%.安徽省亳州市、滁州市、六安市、池州市和淮北市耕地面积变化不大, 均减少不到10%, 以上城市农业活动较多.YRD耕地面积总量2000~2010年迅速减少, 2010年之后耕地面积减少速度放缓, 主要是因为2009年6月国土资源部提出了“保经济增长、保耕地红线”行动[35].
表 4
(Table 4)
表 4 2000~2018年长三角41个城市耕地面积/km2
Table 4 Cultivated land area of 41 cities in YRD from 2000 to 2018/km2
省份 |
城市 |
2000年 |
2005年 |
2010年 |
2015年 |
2018年 |
上海市 |
上海市 |
4 222 |
3 825 |
3 011 |
2 810 |
2 784 |
江苏省 |
南京市 |
4616 |
4081 |
3 626 |
3 651 |
3 592 |
无锡市 |
2332 |
2113 |
1 715 |
1 693 |
1 689 |
镇江市 |
2740 |
2518 |
7 511 |
7 541 |
7 503 |
常州市 |
2868 |
2634 |
2 324 |
2 296 |
2 254 |
苏州市 |
3599 |
3172 |
2 530 |
2 413 |
2 411 |
南通市 |
6369 |
5999 |
5 482 |
5 464 |
5 466 |
连云港市 |
5320 |
5312 |
5 091 |
5 113 |
5 092 |
淮安市 |
6534 |
6449 |
6 184 |
6 213 |
6 056 |
盐城市 |
11 049 |
11 001 |
10 233 |
10 086 |
10 162 |
扬州市 |
4 326 |
4 136 |
3 773 |
3 749 |
3 689 |
徐州市 |
7 846 |
7 749 |
2 219 |
2 218 |
2 215 |
泰州市 |
3 979 |
3 848 |
3 414 |
3 143 |
3 074 |
宿迁市 |
5 737 |
5 596 |
5 383 |
5 430 |
5 398 |
小计 |
67 314 |
64 608 |
59 486 |
59 008 |
58 601 |
浙江省 |
杭州市 |
2 644 |
2 310 |
2 098 |
2 057 |
2 076 |
宁波市 |
2 704 |
2 521 |
2 285 |
2 296 |
2 294 |
嘉兴市 |
2 830 |
2 733 |
2 494 |
2 424 |
2 411 |
湖州市 |
2 223 |
2 021 |
1 821 |
1 683 |
1 658 |
绍兴市 |
2 641 |
2 537 |
2 377 |
2 365 |
2 371 |
舟山市 |
272 |
252 |
230 |
220 |
237 |
温州市 |
2 102 |
2 073 |
1 934 |
1 916 |
1 932 |
金华市 |
3 119 |
3 035 |
2 829 |
2 797 |
2 792 |
衢州市 |
1 851 |
1 875 |
1 765 |
1 736 |
1 731 |
台州市 |
2 358 |
2 255 |
1 883 |
1 990 |
1 971 |
丽水市 |
2 390 |
2 270 |
2 184 |
2 145 |
2 157 |
小计 |
25 133 |
23 882 |
21 899 |
21 629 |
21 629 |
安徽省 |
合肥市 |
9 038 |
8 675 |
8 349 |
8 205 |
8 216 |
淮北市 |
2 303 |
2 254 |
2 219 |
2 172 |
2 150 |
亳州市 |
7 105 |
7 107 |
7 064 |
6 994 |
6 965 |
宿州市 |
7 354 |
7 246 |
7 126 |
7 062 |
7 071 |
蚌埠市 |
4 750 |
4 609 |
4 593 |
4 582 |
4 563 |
阜阳市 |
7 755 |
7 716 |
7 536 |
7 482 |
7 265 |
淮南市 |
2 027 |
1 918 |
1 803 |
1 795 |
1 787 |
滁州市 |
9 862 |
9 803 |
9 379 |
9 311 |
9 231 |
六安市 |
9 525 |
9 687 |
9 470 |
9 363 |
9 335 |
马鞍山市 |
3 064 |
2 957 |
2 716 |
2 657 |
2 534 |
芜湖市 |
4 474 |
4 427 |
4 290 |
4 197 |
4 150 |
宣城市 |
4 617 |
4 597 |
4 479 |
4 380 |
4 309 |
铜陵市 |
682 |
652 |
643 |
616 |
607 |
池州市 |
2 373 |
2 340 |
2 281 |
2 257 |
2 185 |
安庆市 |
6 794 |
6 737 |
6 654 |
6 588 |
6 521 |
黄山市 |
1 206 |
1 172 |
1 163 |
1 108 |
1 081 |
小计 |
82 930 |
81 896 |
79 765 |
78 768 |
77 971 |
合计 |
179 599 |
174 211 |
164 161 |
162 216 |
160 986 |
|
表 4 2000~2018年长三角41个城市耕地面积/km2
Table 4 Cultivated land area of 41 cities in YRD from 2000 to 2018/km2
|
2.2 2000~2018年农田生态系统NH3排放变化
图 5所示为基于土地利用变化和氮肥施用量变化计算的2000~2018年YRD地区农田生态系统NH3排放的变化.从中可见, 农田生态系统NH3排放较高的城市主要集中在苏北地区(盐城、淮安、宿迁、连云港和徐州市), 该区域地势平缓、水源丰沛, 利于农业发展, 耕地较多、农业活动相对较多.而浙江省城市整体排放较少, 因为该省城市化水平相对较高, 西南部多为山地丘陵, 难以形成集约化、规模化的农作物种植基地.安徽省则呈现出明显的区域差异, 皖北排放明显高于皖南, 淮河平原区地势平坦, 河网密集, 耕地占比大, NH3排放多; 南部地势较高, 多山地, 林地和草地较多, NH3排放明显少于北部.
由图 6可知, 2000~2018年YRD大部分城市的农田生态系统NH3排放量都在减少, 只有连云港市、淮安市、盐城市、嘉兴市、绍兴市、淮南市、马鞍山市、芜湖市和铜陵市NH3排放量有所增加, 氮肥施用量的增加是NH3排放量增加的主要原因.其中, 铜陵市和淮南市NH3排放量出现大幅增长的原因主要是氮肥施用总量明显提高所致.18年来, 部分城市农田生态系统NH3排放量呈现先增长后降低的趋势, 连云港、淮安、盐城、宿迁市和淮南等城市NH3排放量在2000~2010年均持续增高, 连云港市、盐城市、宿迁市、宣城市、池州市、六安市和安庆市在2010~2015年达到峰值随后降低, 淮安市和淮南市在2015~2018年达到峰值随后降低.有15个城市NH3排放量在2010~2015年出现峰值, 11个城市NH3排放量在2015~2018年出现峰值, 这种情况大多是由于氮肥施用引起.41个城市中上海、苏州、徐州、合肥和亳州年排放量均大幅减少, 减少量均在10 kt以上.此外, 温州市的年减排比最大, 减少了85%.但是, 18年来淮南市农田生态系统NH3排放量增高极大, 年排放量增加了23 305 t(299%), 主要是因为2010年以后该市在农业生产过程中大量增加氮肥用量, 18年来该地区氮肥施用量(折纯)增加了81 801 t·a-1.
江苏省的NH3排放量整体大幅降低, 在3个省中减排最大, 之所以出现这种趋势, 主要原因为江苏省农委自2015年制定了化肥使用量零增长行动计划[36], 近十年来全省氮肥施用量下降.随着技术的发展和农业生产结构的改善, 农业生产对肥料的需求逐渐减少[37], 部分城市氮肥施肥率校正系数有所降低, 2015年后氮肥施用量急剧减少, 2018年NH3排放量达到最低水平.浙江省全省减排最小, 但减排比例最高, 主要是因为该地区耕地减少最明显且氮肥施用减少力度最大.安徽省之所以减排较低, 是因为该省为农业大省, 农业仍作为重要的经济组成部分, 城市化进程较缓, 许多城市农业活动密集, 氮肥施用量及耕地面积减少不多.
由图 6和图 7可见, 尽管部分城市氮肥施肥率增多, 但由于耕地面积大量减少导致最终氮肥施用总量在减少, 故氮肥施用源NH3排放量占农田生态系统NH3排放的比例呈现减小的趋势, 而土壤本底源NH3排放量占比随之增长.近20年来, 秸秆堆肥源NH3排放占比持续增高, 主要是因为农业技术提高带来的农作物产量大幅提高及秸秆堆肥率的增高.
2.3 土地利用变化对农田生态系统NH3排放变化的影响
由表 5~7可知, 近20年来YRD地区土地利用变化对农田生态系统NH3排放的影响持续增大.相较于2000年, 2005年由于土地利用变化导致农田生态系统NH3排放量仅减少了23.2 kt(7.9%).到2010年土地利用变化产生的影响明显增大, 减少了66.8 kt(22.7%), 随后影响略有增大. 2000~2015年YRD大部分城市氮肥施肥率增加, 并且秸秆产率及堆肥率有所提高, 因此, 氮肥施用变化导致农田生态系统NH3排放量持续增长.但是2015年以后NH3排放量反而减小, 主要是因为随着氮肥品质和种植技术的提高, 以及氮肥施用量管控行动计划的实施, YRD氮肥施肥率降低.2015~2018年, 氮肥施用变化所推动的农田生态系统NH3减排最明显, 由6.9 kt(2.3%)增长到69.5 kt(23.6%), 主要因为该时间段内绝大部分城市氮肥施肥率明显减少所致.2000年以来, 土地利用变化和氮肥施用变化共同导致农田生态系统NH3排放量持续下降, 相较于2000年, 2015年农田生态系统NH3排放量下降27.2%, 其中土地利用变化和氮肥施用量的影响贡献分别为91.2%和8.8%; 2018年农田生态系统NH3排放量下降49.8%, 其中土地利用变化和氮肥施用量的影响贡献分别为52.6%和47.3%.表明伴随着快速城镇化, 耕地面积的减少带来NH3排放量明显下降; 同时, 随着农业作业方式的转变以及对氮肥施用量的控制, NH3的排放也得到了显著控制.
表 5
(Table 5)
表 5 9种排放情景模拟结果/kt
Table 5 Simulation results of nine emission scenarios/kt
情景 |
S1 |
S2 |
S3 |
S4 |
S5 |
S6 |
S7 |
S8 |
S9 |
参数 |
729 |
705 |
708 |
662 |
704 |
655 |
649 |
651 |
582 |
|
表 5 9种排放情景模拟结果/kt
Table 5 Simulation results of nine emission scenarios/kt
|
表 6
(Table 6)
表 6 情景模拟结果比较/kt
Table 6 Comparison of scenario simulation results/kt
项目 |
2005年 |
| 2010年 |
| 2015年 |
| 2018年 |
S2-S1 |
S3-S2 |
S3-S1 |
| S2-S1 |
S3-S2 |
S3-S1 |
| S6-S1 |
S7-S6 |
S7-S1 |
| S8-S1 |
S9-S8 |
S9-S1 |
参数 |
-23.2 |
2.1 |
-21.1 |
| -66.8 |
42.6 |
-24.2 |
| -73.2 |
-6.9 |
-80.1 |
| -77.4 |
-69.5 |
-146.9 |
|
表 6 情景模拟结果比较/kt
Table 6 Comparison of scenario simulation results/kt
|
表 7
(Table 7)
表 7 2005~2018年相较于2000年在不同影响因素下对农田生态系统氨排放变化的影响/%
Table 7 Changes in ammonia emissions from agriculture under different influencing factors during 2005-2018 compared to 2000/%
年份 |
土地利用 |
氮肥施用 |
土地利用+氮肥施用 |
贡献率(土地利用∶氮肥施用) |
2000 |
1 |
1 |
1 |
/ |
2005 |
-7.9 |
0.7 |
-7.2 |
/ |
2010 |
-22.7 |
14.5 |
-8.2 |
/ |
2015 |
-24.8 |
-2.3 |
-27.2 |
91.2∶8.8 |
2018 |
-26.2 |
-23.6 |
-49.8 |
52.6∶47.3 |
|
表 7 2005~2018年相较于2000年在不同影响因素下对农田生态系统氨排放变化的影响/%
Table 7 Changes in ammonia emissions from agriculture under different influencing factors during 2005-2018 compared to 2000/%
|
2.4 不确定性分析
土地利用类型数据的获取上, 利用红外航空遥感图像解译具有一定的不确定性.NH3排放的估算也具有一定的不确定性, 主要是因为精确活动水平获取的困难性和排放因子数据匹配的问题.本研究农田生态系统EF(NH3)来自文献[34], 根据该地区实际情况校正排放因子, 参数的选择还无法做到完全本地化, 因此最终采用的排放因子会在一定程度上影响NH3排放量的估算结果.NH3排放活动水平中的氮肥施用量、固氮植物种植面积、秸秆堆肥的作物产量来自各市统计年鉴(部分来自农业农村局公布数据或省统计年鉴), 由于个别城市氮肥施用量通过推算而来、固氮植物和秸秆堆肥数据不全, 各类氮肥使用比例选取相关文献中的数据, 以上因素对各来源的NH3排放量均具有一定影响; 耕地面积来自1∶50 000红外航空遥感影像解译, 在土地类型的分辨上存在一定的误差, 因此, 对土地利用变化和土壤本底源NH3排放量均有一定影响.另外, 本文在情景设置上虽然尽可能将土地利用的影响分离出来, 但是不能做到将其他间接潜在影响因素完全剔除.
3 结论
(1) 伴随着快速城市化进程, YRD耕地面积逐渐减少, 耕地面积从2000年的276 269 km2(49%)减少到2018年的244 001 km2(44%).土地利用变化对农田生态系统NH3排放的影响显著, 表现在直接带来土壤本底NH3排放变化, 以及间接导致氮肥施用量变化进一步带来NH3排放变化.2000~2018年, YRD氮肥施用产生的NH3排放量从690 kt·a-1减少到541 kt·a-1(减少22%); 土壤本底产生的NH3排放量从32 kt·a-1减少到29 kt·a-1(减少9%).
(2) 由于土地利用变化所驱动的NH3减排在2005~2010年十分明显, 2010年以后减排不大.相较于2000年, 由于氮肥施用变化导致2010年农田生态系统NH3排放增长了14.5%, 2015年随着管控力度的加强急剧变化, 减少了2.3%, 随后大幅减少, 2018年减排了23.6%.
(3) 近20年来YRD城镇化进程加快及农业技术提高, 土地利用和氮肥施用变化显著, 导致农田生态系统NH3排放量大幅减少.相较于2000年, 二者共同导致农田生态系统NH3排放减少比例在2005年和2010年仅为7.2%和8.2%, 到2018年达到49.8%, 其中土地利用变化和氮肥施用的影响贡献分别为52.6%和47.3%.研究结果表明, 伴随着快速城镇化, 耕地面积的减少带来NH3排放量明显下降; 同时, 随着农业作业方式的转变以及对氮肥施用量的控制, NH3的排放也得到了显著控制.
(4) 由于农业源NH3对于大气PM2.5贡献突出, 而农业源NH3排放的主要来源之一为氮肥施用, 建议控制氮肥施用量、优化氮肥施用方式, 进一步降低NH3排放水平.