旧版本|1980年1月1日~2024年12月31日中国各省市区县地表PM2.5质量浓度

新版本在这里:https://rstata.duanshu.com/#/brief/course/67f02f4cf0cd496f8b643f8308bb6284

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

1980年1月~2024年12月各省市区县地表 PM2.5 质量浓度:https://rstata.duanshu.com/#/course/76698bc838ef49c5ab57633aa64b4e06

今天再为大家分享一份日度数据:1980年1月1日~2024年12月31日中国各省市区县地表 PM2.5 质量浓度。

数据源为 NASA 里面以小时为度量单位的卫星反演数据:https://disc.gsfc.nasa.gov/datasets/M2T1NXAER_5.12.4/summary。里面包含这些常用指标:

  1. 地表PM2.5质量浓度
  2. 地表SO4质量浓度
  3. 地表SO2质量浓度
  4. 地表黑碳质量浓度
  5. 地表有机碳质量浓度
  6. 地表粉尘质量浓度
  7. 地表二甲硫醚质量浓度

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

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

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

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

根据 NASA 网站上 FAQ 里面的指引,PM2.5 的浓度应该按照如下方法计算:

PM2.5 = DUSMASS25 + OCSMASS+ BCSMASS + SSSMASS25 + SO4SMASS * (132.14/96.06)

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

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

数据概览

为了方便大家使用,我把数据汇总成了中国各省份、各城市、各区县的地表 PM2.5 质量浓度面板数据。数据预览如下:

  • 1980年1月1日~2024年12月31日各省份地表 PM2.5 质量浓度日度数据(微克每立方米)

  • 2024年各城市地表 PM2.5 质量浓度日度数据(微克每立方米)

  • 2024年各区县地表 PM2.5 质量浓度日度数据(微克每立方米)_第一部分

绘图展示

中国各省 PM2.5 浓度日度面板数据

对中国范围的栅格数据数据分省份(使用 2021 年行政区划数据)裁剪汇总即可得到各省的面板数据(1980年1月1日~2024年12月31日),下图展示了 2024 年中国各省的年均 PM2.5 浓度分布(从日度数据平均得到,单位:微克/立方米):

中国各市 PM2.5 浓度日度面板数据

对中国范围的栅格数据数据分城市(使用 2021 年行政区划数据)裁剪汇总即可得到各市的面板数据(1980年1月1日~2024年12月31日),下图展示了 2024 年中国各市的年均 PM2.5 浓度分布(从日度数据平均得到,单位:微克/立方米):

中国各区县 PM2.5 浓度日度面板数据

对中国范围的栅格数据数据分区县(使用 2021 年行政区划数据)裁剪汇总即可得到各区县的面板数据(1980年1月1日~2024年12月31日),下图展示了 2024 年中国各区县的年均 PM2.5 浓度分布(从日度数据平均得到,单位:微克/立方米):

中国历年的年均 PM2.5 浓度变化

我们还可以直接计算全国历年的年均 PM2.5 浓度变化:

处理方法

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

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

点击这里跳转到 RStata 短书平台获取附件:旧版本|1980年1月1日~2024年12月31日中国各省市区县地表PM2.5质量浓度

评论