达卡PM2.5负载水溶性无机离子的综合来源解析：PCA与PMF相结合的方法

《Atmospheric Environment: X》：Comprehensive Source Apportionment of PM 2.5-bound Water-Soluble Inorganic Ions in Dhaka: A Combined PCA and PMF Approach

　　

摘要
空气污染已成为达卡一个严重的环境问题，该市始终被评为世界上污染最严重的城市之一。本研究描述了2023年11月至2024年3月期间在达卡大学一处具有代表性的城市屋顶地点采集的细颗粒物（PM_2.5）中水溶性无机离子的化学成分和来源，研究发现主要阳离子，尤其是Na⁺、K⁺、Ca²⁺，以及阴离子，如F^-、Cl^-、NO₃^-、SO₄^2-，具有较高的污染水平，在12月和1月最为严重，超过世界卫生组织（WHO）限值的11-15倍。来源识别整合了Pearson相关筛选和仅作为探索性诊断步骤使用的主成分分析（PCA），以及使用EPA PMF 5.0受体模型进行的定量来源解析。分析成功解析出四个不同的来源，包括区域性二次气溶胶生成、当地富氟工业排放、地壳与建筑扬尘，以及海洋/生物质燃烧的混合影响。这些结果得到了HYSPLIT后向轨迹分析的佐证，该分析揭示冬季存在来自西北和西部的强跨境流入，为全球关于孟加拉国冬季PM_2.5动态和空气质量管理科学知识提供了新的科学数据。据研究人员所知，这项研究是迄今为止在达卡为期五个月期间对主要阳离子和阴离子进行的最全面的连续记录之一，并且是首批同时应用PCA、PMF和HYSPLIT分析于达卡水溶性无机离子（WSII）数据集的研究之一，识别出当地氟化物来源。

论文解读：达卡PM_2.5水溶性无机离子的综合来源解析

**一、研究背景、现存问题与研究必要性**

空气污染是人类文明以来最重要的环境问题之一，工业革命以来其影响日益加剧。达卡市的空气质量长期不佳，主要归因于人类活动的增加和工业的发展。根据IQAir发布的2024年世界空气质量报告，达卡在全球首都城市中PM_2.5年均浓度排名第三，达到78 μg/m³，是世界卫生组织（WHO）最新建议值（5 μg/m³）的15倍以上。达卡分区2023年人口约2320万，城市在极少的监管控制下快速扩张，导致严重的交通拥堵，机动车辆和砖窑分别约占城市细颗粒物（PM_2.5）排放的22%和36%。PM_2.5可深入肺部及血管，导致呼吸系统和心血管系统疾病死亡率上升，每增加10 μg/m³会显著增加缺血性心脏病和总体死亡风险。研究指出，空气污染使孟加拉国人预期寿命缩短6.7年，2005年至2018年间达卡有24000例过早死亡。尽管已有部分研究对达卡气溶胶特征进行了分析，但往往受限于时间窗口狭窄或化学物种覆盖不全。例如，近期一项研究虽提供了水溶性无机阴离子的见解，但仅限于2023年2月这一单一过渡月份。相比之下，本研究将采样时长扩展至五个月（2023年11月至2024年3月），涵盖冬季至季风前期，同时定量了主要阴离子和阳离子，并结合Pearson相关筛选、定量PMF受体模型及HYSPLIT后向轨迹案例研究，以提供该地区更完整的二次无机气溶胶（SIA）过程和跨境传输表征，为全球污染最严重的大城市之一的空气质量管理提供科学依据。

**二、研究内容与主要结论**

研究人员于2023年11月5日至2024年3月18日在孟加拉国达卡大学Mukarram Hussain Khundker Biggyan Bhaban大楼屋顶（距地面约45米，坐标23°43'42.70" N, 90°23'54.90" E）使用主动空气采样器（型号TH-150H）采集PM_2.5样品，采样持续24小时（上午8:00至次日上午8:00），流量为70 L/min，共获得22个有效样品。研究通过离子色谱（IC）定量了主要阳离子（Na⁺、K⁺、Ca²⁺）和阴离子（SO₄^2-、NO₃^-、Cl^-、F^-），并运用Pearson相关分析、PCA（仅作探索性步骤）和EPA PMF 5.0受体模型进行来源解析。研究发现PM_2.5浓度约为WHO 24小时指导值（15 μg/m³）的11-15倍，约为孟加拉国国家环境空气质量标准（NAAQS）24小时限值（65 μg/m³）的3-3.5倍。PMF模型成功解析出四个来源因子：因子1（贡献40.89%方差）为区域二次无机气溶胶（SIA），以硫酸盐和硝酸盐为主导；因子2（16.33%）为当地富氟工业燃烧排放；因子3（14.71%）为地壳和建筑扬尘，钙离子占主导；因子4（28.07%）为海洋/生物质燃烧混合来源，富含钠、钾和氯离子。HYSPLIT后向轨迹分析表明，冬季高污染事件中气团主要来自西北和西部方向，经印度次大陆（特别是印度河-恒河平原）长距离传输至达卡，证实了跨境传输对达卡冬季PM_2.5污染的重要贡献。该研究的科学意义在于提供了达卡冬季至季风前期五个月的连续WSII数据集，首次同时应用PCA、PMF和HYSPLIT三种方法于达卡WSII数据集，并识别出当地氟化物工业来源，为区域空气质量管理提供了重要的科学依据。

**三、主要关键技术方法**

本研究采用以下关键技术方法：（1）基于达卡大学屋顶（约45米高度）采样点，使用TH-150H高流量采样器配合PM_2.5切割头，以70 L/min流量在石英滤膜上连续采集24小时PM_2.5样品，共22个有效样品，并设置野外空白对照；（2）采用离子色谱（SYKAM S150系统）定量分析水溶性无机离子，阴离子（F^-、Cl^-、NO₃^-、SO₄^2-）检测限为0.033-0.098 μg/m³，阳离子（Na⁺、K⁺、Ca²⁺）检测限低于0.05 μg/m³；（3）使用NOAA HYSPLIT模型（版本5.3.2）计算24小时后向轨迹，起始高度为100、500和1000米，用于识别气团来源和传输路径；（4）采用Pearson相关分析（并用Spearman秩相关进行交叉验证）识别离子间线性关系；（5）应用主成分分析（PCA）作为探索性步骤，检验协方差结构；（6）采用EPA PMF 5.0正矩阵因子分解受体模型进行定量来源解析，通过100次Bootstrap和DISP分析验证模型稳定性（因子映射一致性>80%，0%因子交换）。

**四、研究结果与讨论**

**（一）达卡空气质量状况**

研究期间达卡AQI（空气质量指数）水平持续处于“不健康”至“非常不健康”等级。美国大使馆监测站数据显示，2023年11月至2024年3月各月平均PM_2.5浓度分别为175.77 ± 40.12、225.81 ± 44.15、259.65 ± 34.27、224.34 ± 45.83和186.48 ± 29.53 μg/m³，其中1月AQI平均值最高（252.94 ± 17.21）。本研究的PM_2.5实测值（约220-250 μg/m³）与美国大使馆数据（约176-260 μg/m³）高度吻合，两者均远超WHO限值11-17倍，验证了屋顶采样点的代表性。相对湿度（62%-73%）和低温等气象条件促进了冬季逆温层对PM_2.5的捕获和累积。

**（二）HYSPLIT轨迹影响**

HYSPLIT后向轨迹分析显示显著的季节差异：2023年11月（8日）和12月（19日）高污染事件期间，轨迹主要来自西北至西方向，穿过印度河-恒河平原和印度北部，表明长距离大陆传输。2024年1月和2月的轨迹主要经过印度中部和东部，传输距离较短，表明区域尺度和本地排放的共同影响。2024年3月18日，PM_2.5浓度显著降低，气团轨迹从东南偏南方向弯曲，经过孟加拉湾海域，较清洁的海洋气团对大陆污染羽流起到稀释作用。这些结果与同期AQI和PM_2.5变化相吻合，证实印度次大陆跨境传输是达卡冬季高PM_2.5污染事件的重要驱动因素。

**（三）PM_2.5中WSII浓度和分布**

研究表明，SO₄^2-是最丰富的离子，平均浓度为20,581.05 ± 7,180.25 ng/m³，显著高于2013年的5,300 ng/m³和2023年2月的11,039 ng/m³；NO₃^-平均浓度达14,166.65 ± 10,108.42 ng/m³，较2013年的7,750 ng/m³大幅上升。与全球城市相比，达卡的二次无机气溶胶（SIA）负担远超欧洲和沿海亚洲城市（如伯尔尼、横滨、新加坡），接近武汉、上海、广州等工业中心城市水平。Cl^-浓度从2013年的1,190 ng/m³升至当前的5,014.05 ± 3,371.81 ng/m³，K⁺浓度（3,610.73 ± 2,125.50 ng/m³）为所有对比站点最高值，Ca²⁺（3,697.66 ± 1,870.27 ng/m³）和F^-（106.13 ± 48.79 ng/m³）也处于较高水平。质量百分比分析显示，SO₄^2-占WSII总质量的40.42%，NO₃^-占29.18%，Cl^-占10.52%，Ca²⁺占7.48%，K⁺占7.31%，Na⁺占4.88%。

**（四）相关性分析**

Pearson相关矩阵（以Spearman秩相关交叉验证）揭示了三类重叠的来源类型：（1）盐性/海盐组分，以Na⁺和Cl^-为主，伴生K⁺和部分NO₃^-；（2）二次无机气溶胶组分，NO₃^-和SO₄^2-与Cl^-显著相关，SO₄^2-与Na⁺、K⁺无显著关联，表明其非海洋来源；（3）地壳/道路扬尘组分，以Ca²⁺为主要特征，与F^-、NO₃^-、Na⁺、K⁺呈弱但显著的相关性。

**（五）PCA结果**

PCA将七个标准化离子变量降至两个主成分（PC1 = 64.5%，PC2 = 18.4%，累计方差82.9%）。Ca²⁺、F^-、Na⁺和K⁺在PC1上载荷最强，SO₄^2-、NO₃^-和Cl^-在PC2上载荷最强。该双轴模式与Pearson相关矩阵确定的离子聚类总体一致。但PCA仅用作探索性工具，其载荷模式受到3个极端浓度样本（样本8、9和11）的显著影响，且PCA不施加非负约束、不考虑测量不确定度，因此不能用于识别和量化独立排放源。

**（六）PMF来源解析**

PMF模型通过系统测试3至6个因子，最终选定4因子模型，判定标准为：（i）Q_true/Q_theoretical比值接近1；（ii）模型无因子分裂；（iii）因子结构旋转稳定性高。100次Bootstrap运行证实因子高度可重现（因子1：93%，因子2：89%，因子3：84%，因子4：89%），DISP分析显示0%因子交换，最大ΔQ下降-0.665（<1.0）。解析出的四个来源为：因子1（40.89%）为区域二次无机气溶胶（SIA），由车辆、电厂等燃烧过程排放的NO₂和SO₂经大气氧化形成；因子2（16.33%）以高浓度氟化物为特征，归因于燃煤工业、砖窑、陶瓷和玻璃制造等高温工业过程；因子3（14.71%）以钙为主导，为地壳和建筑扬尘来源，反映道路和土壤扬尘再悬浮及密集建筑施工活动；因子4（28.07%）富含Na⁺、K⁺和Cl^-，代表海洋气团与生物质燃烧（包括农作物残渣露天焚烧和家庭生物燃料使用）的混合影响。

**（七）研究局限与未来展望**

研究的主要局限包括：样本量较小（n = 22），未来应扩展至完整年周期并采用高频率自动化监测；缺乏高分辨率本地工业排放清单，氟化物源解析主要依赖文献对比，建议对当地砖窑、陶瓷和玻璃制造设施进行直接采样和化学指纹分析；HYSPLIT轨迹分析仅为定性案例研究，未来应结合多年、更长持续时间（72-120小时）的集合轨迹以及潜在源贡献函数（PSCF）和浓度加权轨迹（CWT）等空间统计模型，以建立定量跨境传输预算。

**五、研究结论**

该研究对2023-2024年高污染月份达卡PM_2.5中水溶性无机离子（WSII）进行了初步评估。研究结果表明，四个关键因子对WSII有贡献：由NO_x和SO₂气体排放导致的区域二次硫酸盐和硝酸盐生成、高氟化物浓度的当地工业燃烧、再悬浮地壳和建筑扬尘，以及海洋气团和生物质燃烧的协同效应。Pearson相关筛选和PCA仅用作初步探索步骤，本文呈现的定量来源解析完全基于PMF模型，该模型不同于PCA，它纳入了非负约束和测量不确定度，因此能在更符合物理实际和质量平衡的条件下解析这些混合的自然和人为来源。此外，气团轨迹案例研究证实，印度次大陆产生的跨境空气污染对达卡的空气质量有显著影响，这与本地空气污染相结合，加剧了冬季的停滞状况。总之，这些结果凸显了对本地工业污染源采取特定监管干预措施以及开展区域合作的必要性，以减轻细颗粒物对该地区居民造成的严重健康威胁。