今天给大家分享一份 1980 年 1 月到 2025 年 12 月中国各省市区县的 SO2 质量浓度月度数据。由 NASA 的 M2TMNXAER_5.12.4 数据处理得到（卫星反演数据）。处理方法类似之前介绍的逆温数据处理：https://rstata.duanshu.com/#/course/55b47a7a83124138b918e7475f4141d9。
数据概览
为了方便大家使用，我把数据汇总成了全国、各省份、各城市、各区县的地表 SO2 质量浓度面板数据。
- 1980年1月~2025年12月各城市地表SO2质量浓度(微克每立方米).xlsx
- 1980年1月~2025年12月各区县地表SO2质量浓度(微克每立方米)_第二部分.xlsx
- 1980年1月~2025年12月各区县地表SO2质量浓度(微克每立方米)_第一部分.xlsx
- 1980年1月~2025年12月各省份地表SO2质量浓度(微克每立方米).xlsx
- 1980年1月~2025年12月全国地表SO2质量浓度(微克每立方米).xlsx
数据预览如下：
绘图展示
中国各省 SO2 质量浓度
对中国范围的栅格数据数据分省份（使用 2021 年行政区划数据）裁剪汇总即可得到各省的面板数据（1980年1月～2025年12月），下图展示了 2025 年中国各省的年均 SO2 浓度分布（从月度数据平均得到，单位：微克/立方米）：
中国各市 SO2 质量浓度
对中国范围的栅格数据数据分城市（使用 2021 年行政区划数据）裁剪汇总即可得到各市的面板数据（1980年1月～2025年12月），下图展示了 2025 年中国各市的年均 SO2 浓度分布（从月度数据平均得到，单位：微克/立方米）：
中国各区县 SO2 质量浓度
对中国范围的栅格数据数据分区县（使用 2021 年行政区划数据）裁剪汇总即可得到各区县的面板数据（1980年1月～2025年12月），下图展示了 2021 年中国各区县的年均 SO2 浓度分布（从月度数据平均得到，单位：微克/立方米）：
连续性检查
为了对数据质量进行检查，我们可以观察下数据的连续性：
连续性还不错！
处理方法
这份数据的处理类似之前介绍的逆温数据，所以可以参考那个课程学习：
逆温数据是如何处理的？使用 R 语言完成整个过程：https://rstata.duanshu.com/#/course/55b47a7a83124138b918e7475f4141d9
代码
为方便大家学习，附件中还提供了绘图代码以供参考。
数据引用格式
由于该数据包含较多 RStata 处理的内容，在研究中使用该数据请使用清晰的方式注明数据来源于 RStata 或者 RStata 数据中心，并使用如下格式引用：
RStata 数据中心: 1980 年 1 月～2025 年 12 月各省市区县地表 SO2 质量浓度. 2026. https://tidyfriday.cn/rsdb2/
英文文献可以使用下面的格式引用：
RStata Data Center: Panel Data on Surface SO₂ Mass Concentration by Province, City, District, and County in China, January 1980–December 2025. 2026. https://tidyfriday.cn/rsdb2/
点击这里跳转到 RStata 短书平台获取附件：1980 年 1 月～2025 年 12 月各省市区县地表 SO2 质量浓度
评论