将15N纳入多分辨率排放清单以模拟中国珠江三角洲地区排放NOx中δ15N的时空变化

《Atmosphere》：Incorporating 15N into the Multi-Resolution Emission Inventory to Simulate the Spatiotemporal Variations of δ15N in Emitted NOx over the Pearl River Delta Region, China

　　

氮氧化物（NO_x）由一氧化氮（NO）和二氧化氮（NO₂）组成，是大气硝酸盐（NO₃^?）的重要前体物；后者是细颗粒物（PM_2.5）的主要组分之一，并对空气质量、人类健康和生态系统产生关键影响。排放清单能够提供NO_x来源的精细时空信息，而稳定同位素技术则可为源解析提供额外约束。研究人员采用同位素质量平衡模型，将稳定氮同位素（¹⁴N、¹⁵N）纳入华南地区广泛使用的中国多分辨率排放清单（MEIC）中，重点关注中国城市化和工业化程度最高的地区之一——珠江三角洲（PRD）区域。2008年MEIC清单表明，华南地区NO_x排放具有显著空间异质性，以交通源为主，并主要集中在PRD及其他城市群。随后，研究人员将模拟得到的排放NO_x同位素组成与大气观测结果进行比较，以评估排放源在控制大气硝酸盐（NO₃^?）中的作用。结果发现，模拟的δ¹⁵N(NO_x)值总体上低估了观测到的δ¹⁵N(NO₃^?)值。该差异表明，未来需开展支持¹⁵N的空气质量模拟，更好表征源贡献及大气过程，从而改进源解析、排放清单评估以及对活性氮循环的认识。

该研究发表于《Atmosphere》，围绕华南尤其是珠江三角洲（PRD）地区NO_x排放的稳定氮同位素特征展开，目标是将排放清单与同位素示踪方法结合，以提升对大气硝酸盐来源及活性氮循环过程的认识。研究背景在于，NO_x是大气硝酸盐（NO₃^?）形成的关键前体，而NO₃^?又是细颗粒物（PM_2.5）的重要组成部分，直接关系空气质量与人体健康；同时，NO₃^?及其他活性氮沉降还会引发酸雨、土壤酸化、水体富营养化和生物多样性下降等生态环境问题。珠江三角洲作为全球高度城市化和工业化区域，受强烈人类活动驱动，NO_x排放强度高，是中国无机氮沉降热点区域之一，因此准确识别其NO_x来源，对制定有效减排策略具有重要意义。现有排放清单虽能提供较高分辨率的源时空分布信息，并广泛用于空气质量模拟、污染归因和治理决策，但不同源部门和区域层面仍存在明显不确定性。与此同时，稳定氮同位素δ¹⁵N因不同NO_x源具有差异化“指纹”而成为重要源解析工具，但此前尚缺乏将δ¹⁵N系统性嵌入排放清单、进而开展大尺度分源模拟的研究。基于这一缺口，研究人员将¹⁵N引入中国多分辨率排放清单（MEIC），在华南区域及PRD内模拟排放NO_x的δ¹⁵N时空变化，并与大气NO₃^?同位素观测对比，评估排放源对区域大气硝酸盐的控制作用。

研究人员采用的关键技术方法主要包括以下几个方面。首先，构建双层空间模拟域，即覆盖华南的外层区域和聚焦PRD的内层区域，两者水平分辨率均为36 km × 36 km。其次，选用2008年MEIC人为源排放清单作为基础数据，并依据分部门周变化和小时变化谱，将月排放量处理为逐小时NO_x排放，同时按行业特征进行空间分配。再次，系统整理文献中不同源部门NO_x的δ¹⁵N端元值，建立源特异性同位素参数库。最后，利用同位素质量平衡模型，将各网格中不同源的NO_x排放通量与对应δ¹⁵N端元耦合，模拟区域排放NO_x的δ¹⁵N分布，并将结果与PRD内广州城市站和鼎湖山农村站的大气NO₃^? δ¹⁵N观测值进行比较。

在结果与讨论部分，论文首先给出“3.1. Spatiotemporal Variability of NO_x Emission Rates and Source Fractions”的结果。该部分表明，2008年华南NO_x排放呈现显著且季节上相对稳定的空间异质性。高排放热点主要集中在PRD及其他城市群，而近海区域、偏远内陆和广大农村地区排放较低；多数网格排放量低于5 t N d^?1，仅少数热点网格超过20 t N d^?1甚至30 t N d^?1。季节差异方面，10—12月部分热点略强，7—9月相对较弱，但季节变化幅度明显小于空间差异，说明排放格局主要受大型人为源空间分布控制。分源结果显示，交通源是华南范围内最重要且空间影响最广的NO_x来源，在PRD、粤东沿海、北部湾、海南岛及部分内陆区域均表现出较高贡献；工业源也具有明显区域重要性，且分布广泛，主要见于工业带和城郊制造业区；电厂源总体占比有限，但在包含大型电厂的少数网格中具有很高贡献；居民源在大多数区域贡献较小，仅在缺乏交通或大型工业竞争源的少数内陆地区表现出中等占比。综合而言，交通源主导了区域NO_x排放格局，工业源次之，电厂源高度局地化，居民源整体作用有限。

随后，“3.2. Simulated Spatiotemporal Variability of δ¹⁵N(NO_x) Values)”部分展示了排放NO_x同位素组成的模拟结果。基于四类主要源部门的网格化排放贡献，研究人员计算得到2008年华南δ¹⁵N(NO_x)约在?8‰至+2‰之间，具有显著空间异质性，且四季总体空间格局较为一致。区域内广泛分布的δ¹⁵N(NO_x)值主要落在?4‰至?2‰之间，仅少数局地网格接近正值。这一分布特征与交通源的显著贡献相一致，因为交通源δ¹⁵N端元范围较宽且中位值偏负。相比之下，?8‰至?4‰的更低δ¹⁵N(NO_x)值主要出现在内陆和西北部区域，指示居民源贡献相对更高。季节上，4—6月和7—9月在部分热点区域略显富集，1—3月则相对更亏损，但这些变化仍明显弱于稳定存在的空间对比。

接着，“3.3. Model–Observation Comparisons”部分对模拟值与观测值进行了比较。由于大气中NO₃^?是更易获取的同位素分析对象，研究人员选取PRD区域两个站点——广州城市站和鼎湖山农村站——的δ¹⁵N(NO₃^?)观测数据，与对应网格中的模拟δ¹⁵N(NO_x)进行对照。广州观测值范围为?3.9‰至+7.9‰，鼎湖山为+5.1‰至+11.1‰。比较结果表明，两地模拟δ¹⁵N(NO_x)普遍低于观测δ¹⁵N(NO₃^?)，但广州站低估程度相对较轻，鼎湖山站低估更为明显。论文指出，鼎湖山偏差较大可能部分与其西侧邻近电厂有关，该电厂显著影响局地NO_x排放，而当前“仅排放”框架尚未纳入大气传输造成的混合效应。研究进一步归纳了模型—观测偏差的三方面来源：其一，源特异性δ¹⁵N端元本身存在不确定性，排放清单也可能存在偏差，且将复杂排放简化为宽泛部门会带来误差；其二，大气传输、混合和沉降过程会弱化区域排放δ¹⁵N信号；其三，NO_x氧化生成NO₃^?过程中可能发生光解同位素效应、动力学同位素效应和平衡同位素效应。由此，该研究明确指出，当前模型并非用于直接预测大气δ¹⁵N(NO₃^?)，而是作为概念验证，证明将¹⁵N纳入NO_x排放清单并模拟排放来源驱动的δ¹⁵N(NO_x)格局是可行的，并可用于评估排放源对PRD大气NO₃^?的控制作用。剩余偏差则进一步说明，大气硝酸盐同位素组成并非仅由局地NO_x排放决定，还受后续传输、混合、沉降和光化学转化共同影响。

讨论部分的核心在于，本研究成功建立了支持¹⁵N的排放清单框架，使排放清单从传统质量通量描述扩展到同位素维度，为区域尺度活性氮排放、迁移与转化研究提供了新的数据基础。研究结果强调，若仅依赖排放强度与源占比，仍不足以完整解释大气硝酸盐的同位素特征；只有进一步耦合大气传输、沉降以及对流层光化学过程，才能更准确地重建NO₃^?的δ¹⁵N组成。因此，该研究的重要意义在于为后续开展过程分辨的空气质量模拟奠定基础，也为排放清单评估和源解析方法改进提供了方向。

研究结论部分可译为：2008年MEIC清单显示，华南地区NO_x排放具有显著但季节上持续存在的空间异质性，热点主要集中于PRD及其他城市群。交通源是主导且空间分布最广的来源，工业排放在区域上也较为重要，电厂排放则呈局地化特征，而居民源贡献比例较小。研究人员进一步结合MEIC分部门NO_x排放与既有研究整理的源特异性δ¹⁵N值，模拟了δ¹⁵N(NO_x)。华南地区尤其是PRD广泛出现的?4‰至?2‰ δ¹⁵N(NO_x)值，与交通源的显著贡献相一致。与PRD内广州城市站和鼎湖山农村站实测δ¹⁵N(NO₃^?)相比，模拟δ¹⁵N(NO_x)明显偏低，表明仅基于排放的框架不足以再现大气NO₃^?同位素组成。这种差异很可能源于排放清单不确定性、大气传输与沉降，以及NO_x氧化过程中同位素分馏的综合影响。未来研究将把¹⁵N引入CMAQ模型，以解析大气活性氮排放、传输、化学转化和沉降过程中的影响。