之前给大家分享过根据卫星反演数据计算的 PM2.5、SO2、黑碳、粉尘、SO4、二甲硫醚等指标，这些指标都是 MERRA-2 数据可以直接计算出来的。那么 PM1 和 PM10 质量浓度如何计算呢？ MERRA-2 的常见问题中给出了 PM1 和 PM10 质量浓度的计算公式：
- PM1 = (1.375 * SO4 + BCphobic + BCphilic + OCphobic + OCphilic + 0.7 * DU001 + SS01 + SS002) * AIRDENS
- PM10 = (1.375 * SO4 + BCphobic + BCphilic + OCphobic + OCphilic + DU001 + DU002 + DU003 + 0.74 * DU004 + SS01 + SS002 + SS003 + SS004) * AIRDENS
MERRA-2 的常见问题：https://gmao.gsfc.nasa.gov/reanalysis/MERRA-2/FAQ/
按照该指引，我们首先找到 aer_Nv 数据的下载地址：
https://disc.gsfc.nasa.gov/datasets/M2I3NVAER_5.12.4/summary?keywords=aer_Nv
然后下载了中国范围的数据等待处理。由于数据量很大，所以今天我们先处理下 PM1 数据。今天分享给大家的就是处理得到的 1980～2025 年各省市区县 PM1 质量浓度日度面板数据。
数据概览
直接下载到的数据为 NetCDF 格式的数据，分辨率为 0.5 ° x 0.625 °，通过计算可以得到 PM1 的栅格数据，然后分省市区县进行裁剪汇总即可得到中国各省市区县的 PM1 质量浓度面板日度数据。其中，城市和区县的数据由于观测值超过了 100 万，所以我按年进行了拆分。附件中包含了如下文件：
- 1980~2025年全国平均PM1质量浓度日度数据.xlsx
- 1980~2025年各省份PM1质量浓度日度数据.xlsx
- 城市数据（文件夹）
- 区县数据（文件夹）
各文件预览如下：
绘图展示
为了大家更加直观的感受这份数据，我使用 Stata 绘制了下面的图表。
各省份 PM1 质量浓度
2025 年 10 月各省份月平均 PM1 质量浓度：
各城市 PM1 质量浓度
2025 年各城市年平均 PM1 质量浓度：
各区县 PM1 质量浓度
2025 年 12 月 31 日各区县日平均 PM1 质量浓度：
连续性检查
最后为了观察数据的连续性，我绘制了各省份历年 PM1 的均值趋势：
连续性还不错！
处理方法
该数据的处理方法类似之前讲解的逆温数据处理：
逆温数据是如何处理的？使用 R 语言完成整个过程：https://rstata.duanshu.com/#/course/55b47a7a83124138b918e7475f4141d9
感兴趣的小伙伴可以参考课程里面的方法学习。
数据引用格式
由于该数据包含较多 RStata 处理的内容，在研究中使用该数据请使用清晰的方式注明数据来源于 RStata 或者 RStata 数据中心，并使用如下格式引用：
RStata 数据中心: 1980～2025 年各省市区县 PM1 质量浓度日度面板数据. 2026. https://tidyfriday.cn/rsdb2/
英文文献可以使用下面的格式引用：
RStata Data Center: Panel Data on Daily PM₁ Mass Concentration by Province, City, District, and County in China, 1980–2025. 2026. https://tidyfriday.cn/rsdb2/
点击这里跳转到 RStata 短书平台获取附件：1980～2025 年各省市区县 PM1 质量浓度日度面板数据
评论