1980年1月1日～2025年12月31日中国各省市区县地表 SO2 质量浓度

之前给大家分享过一份月度的中国各省市区县地表 SO2 质量浓度的卫星反演数据：

1980年1月～2025年12月各省市区县地表 SO2 质量浓度：https://rstata.duanshu.com/#/brief/course/72351a45d60c4c63af3de00db732221c

有些小伙伴表示自己还需要日度的，不过因为那个数据源就是月度，所以没办法处理成日度的。最近有一个小伙伴在 NASA 里面找到了一个小时的卫星反演的数据：https://disc.gsfc.nasa.gov/datasets/M2T1NXAER_5.12.4/summary。

同样，我从中选取了以下几个比较常用的指标进行处理：

地表 PM2.5 质量浓度
地表SO4质量浓度
地表SO2质量浓度
地表黑碳质量浓度
地表有机碳质量浓度
地表粉尘质量浓度
地表二甲硫醚质量浓度

处理方法类似之前介绍的逆温数据处理：

逆温数据是如何处理的？使用 R 语言完成整个过程：https://rstata.duanshu.com/#/course/55b47a7a83124138b918e7475f4141d9

就是先从 NASA 下载 nc 格式的数据，使用 R 语言读取为栅格数据再分区域平均得到各省市区县的平均质量浓度。感兴趣的小伙伴可以学习逆温数据处理的课程，另外最近如果有时间的话也会推出专题课程讲解该数据的处理。

由于每个指标的处理都非常耗费时间（日度数据需要处理上万个文件），因此这次先处理地表 SO2 质量浓度指标。

由于日度数据的栅格数据文件太多了，所以这一次没有准备栅格数据。

另外就是为了让大家了解推文中的图片是如何绘制的，附件中为大家保留了绘制这些图表的 R 语言代码和相关数据。详细绘制方法可以学习课程：

使用 R 语言绘制中国夜间灯光地图：https://rstata.duanshu.com/#/course/a430bddc4100426fb4dfeeec38c994d9

数据概览

为了方便大家使用，我把数据汇总成了中国各省份、各城市、各区县的地表 SO2 质量浓度面板数据。

1980年1月1日~2025年12月31日各省份地表 SO2 质量浓度日度数据(微克每立方米).xlsx
城市数据（文件夹）
区县数据（文件夹）

数据预览如下：

绘图展示

中国各省 SO2 质量浓度

对中国范围的栅格数据数据分省份（使用 2021 年行政区划数据）裁剪汇总即可得到各省的面板数据（1980年1月1日～2025年12月31日），下图展示了 2025 年中国各省的年均 SO2 浓度分布（从日度数据平均得到，单位：微克/立方米）：

中国各市 SO2 质量浓度

对中国范围的栅格数据数据分城市（使用 2021 年行政区划数据）裁剪汇总即可得到各市的面板数据（1980年1月1日～2025年12月31日），下图展示了 2025 年中国各市的年均 SO2 浓度分布（从日度数据平均得到，单位：微克/立方米）：

中国各区县 SO2 质量浓度

对中国范围的栅格数据分区县（使用 2021 年行政区划数据）裁剪汇总即可得到各区县的面板数据（1980年1月1日～2025年12月31日），下图展示了 2025 年中国各区县的年均 SO2 浓度分布（从日度数据平均得到，单位：微克/立方米）：

连续性检查

我们还可以直接计算全国历年的年均 SO2 浓度变化：

连续性还不错！

处理方法

这份数据的处理类似之前介绍的逆温数据，所以可以参考那个课程学习，另外最近如果有时间的话也会推出专题课程讲解该数据的处理。

逆温数据是如何处理的？使用 R 语言完成整个过程：https://rstata.duanshu.com/#/course/55b47a7a83124138b918e7475f4141d9

数据引用格式

由于该数据包含较多 RStata 处理的内容，在研究中使用该数据请使用清晰的方式注明数据来源于 RStata 或者 RStata 数据中心，并使用如下格式引用：

RStata 数据中心: 1980年1月1日～2025年12月31日中国各省市区县地表 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 1, 1980 – December 31, 2025. 2026. https://tidyfriday.cn/rsdb2/

点击这里跳转到 RStata 短书平台获取附件：1980年1月1日～2025年12月31日中国各省市区县地表 SO2 质量浓度