基于精细模式气溶胶与WRF模式估算PM2.5质量浓度
南京师范大学地理科学学院,南京 210023
The estimation of PM2.5 mass concentration based on fine-mode aerosol and WRF model
School of Geography Science, Nanjing Normal University, Nanjing 210023, China
通讯作者: 查 勇(1963-),男,教授,主要研究方向为大气遥感。Email:yzha@njnu.edu.cn。
责任编辑: 李瑜
收稿日期: 2020-07-10 修回日期: 2021-08-31 网络出版日期: 2021-06-15
基金资助: |
|
Received: 2020-07-10 Revised: 2021-08-31 Online: 2021-06-15
作者简介 About authors
韦 耿(1996-),男,硕士研究生,主要研究方向为大气颗粒物质量浓度估算。Email:
将气象要素加入到基于气溶胶光学厚度(aerosol optical depth,AOD)的近地面大气颗粒物浓度估算是目前热门的技术手段之一。获取了江苏省南京市2014年3月—2019年2月期间的AOD,精细模式分数 (fine-mode fraction,FMF)和PM2.5质量浓度数据,并结合天气研究和预报 (weather research and forecast,WRF)模式得到的气象模拟数据,对南京市PM2.5的质量浓度进行反演。结果表明,相比于AOD与PM2.5进行相关性分析,通过FMF校正得到的精细气溶胶光学厚度AODf与PM2.5的相关性分析能够取得更高的拟合系数,R2最高达到0.40。利用随机森林模型,引入含不同高度的气象因子对PM2.5质量浓度建立反演模型,得到的拟合系数与各误差指标均优于仅含近地面气象因子的模型,表明PM2.5质量浓度受到多因子共同作用的影响,能较好地为利用多源数据反演PM2.5质量浓度提供依据和参考。
关键词:
The addition of meteorological factors to the estimation of near-ground atmospheric particulate concentration based on AOD is one of the most popular techniques nowadays. In this paper, AOD (Aerosol Optical Depth), FMF (Fine-Mode Fraction) and PM2.5mass concentration data from March 2014 to February 2019 in Nanjing were obtained, and the mass concentration of PM2.5 in Nanjing was retrieved in combination with the meteorological simulation data from WRF (Weather Research and Forecast) model. The results show that, compared with correlation between AOD and PM2.5, the correlation analysis of fine aerosol optical depth AODf and PM2.5 obtained by FMF correction can obtain a higher fitting coefficient, and the maximum R2 reaches 0.40. By adding meteorological factors on different heights into random forest model to establish an inversion model for PM2.5 mass concentration, the obtained fitting coefficients and various error indicators are better than those from models with only near-surface meteorological factors, which indicates that the PM2.5 mass concentration is affected by the combined effect of multiple factors, thus the result can provide a basis and reference for inversion of PM2.5 mass concentration by using multi-source data.
Keywords:
本文引用格式
韦耿, 侯钰俏, 韩佳媚, 查勇.
WEI Geng, HOU Yuqiao, HAN Jiamei, ZHA Yong.
0 引言
大气颗粒物污染的遥感监测主要是通过对气溶胶的遥感监测进行研究,气溶胶光学厚度(aerosol optical depth,AOD)是气溶胶最重要的物理参数之一,是确定气溶胶气候效应的重要因素。利用AOD与大气颗粒物建立线性回归模型复杂程度较低、应用性较强,是早期估算地面大气颗粒物浓度的最常规方法[11,12]。但部分研究证明气溶胶光学特性的多样性以及粒径分布、组成和形状的变化会影响AOD和大气颗粒物之间的相关性[13,14],大气中细粒子的消光作用对总消光作用中所占比例往往容易被忽视,因此通过引入精细模式分数(fine mode fraction,FMF)可将大气中的细气溶胶与总气溶胶分离,利用FMF校正得到的精细模式气溶胶与大气细颗粒物PM2.5进行相关性分析,能够有效提高气溶胶反演PM2.5的精度[15,16]。气象因子能够影响到气溶胶从形成到消散的完整过程[17],是影响大气颗粒物时空变化较为显著的因素之一。当前诸多研究通过引入多种数据综合分析,证明能够取得更好的相关性,建立更精准的PM2.5反演数学模型。纪晓璐等[18]通过将MODIS气溶胶数据与风速数据相结合,得出环渤海区域AOD年内变化与风速整体呈反比关系,且AOD与空气质量指数相关性较高; Ni等[19]通过MODIS数据结合温度、气压、湿度、能见度、风速等数据,通过一般验证和交叉验证对京津冀地区的PM2.5进行反演,反演精度较好,R2分别达到0.68和0.54; 侯俊雄等[20]通过将天气研究和预报 (weather research and forecast,WRF)模式地表气象因子反演了京津冀地区PM2.5的质量浓度,同样取得较好结果。
尽管引入FMF能够提高通过AOD反演PM2.5质量浓度的精度,但该类型研究尤其是国内相关研究依旧较少,并且大多数多元数据综合分析研究中,仅考虑将近地面气象因子加入分析中,而未考虑到高空气象因子对近地面大气颗粒物质量浓度的影响。本文以南京市为研究区域,通过利用FMF与AOD结合,讨论FMF对PM2.5质量浓度估算的影响,并将通过WRF模式得到的气象因子数据(包括近地面及高空的温度、风速与相对湿度,和边界层高度)作为辅助因子,分为4个季节加入到AOD对PM2.5的反演中,分析和讨论引入近地面和不同高度的气象因子对AOD反演大气颗粒物的可行性。
1 研究区概况与数据源
1.1 研究区概况
选择江苏省南京市作为研究区,南京市是中国江苏省政治、文化中心,是中国东部地区重要的经济中心和工业化城市之一,该地区属亚热带季风气候,四季分明,属于典型的低山丘陵地貌。随着城市化和工业化的加速发展,导致汽车尾气、工业污染等人为污染物排放不断增加,并且外来沙尘、化石燃料燃烧产生的烟尘又会进一步提高本地大气污染程度,空气质量问题较为严峻。
1.2 MODIS卫星数据
与SPOT,TM等卫星影像相比,MODIS传感器有更多的传感器通道,更好的动态性,且空间覆盖广、可靠性高[21],数据由美国国家航空航天局每日向全球免费发布(
1.3 大气颗粒物数据
南京市目前共有9个国控级别空气质量检测站点,包含多种污染气体的质量浓度数据,监测频率为1 h,监测时间为每日0—23时,并对全社会发布每小时的整点数据。PM2.5质量浓度数据来自中国生态环境部全国城市空气质量实时发布平台。对于出现部分数据缺测的现象,若缺失值较少,则采用缺失值前后最相近的数据进行补充。若缺失值前后6 h内均没有有效数据,则不使用该时刻的数据集。
1.4 WRF气象数据
在某特定时间区域内大气环境的空气污染程度主要受到包括气象因素、地形地貌等自然因素影响。对于某固定区域而言,地形地貌处于稳定状态,因此大气污染程度主要受垂直自净能力(边界层高度)、水平自净能力(水平风速)、沉降(降水、湍流作用、重力沉降等)、湿度等气象因素影响。
WRF模式是由美国国家环境预报中心及美国国家大气研究中心等机构开发的新一代中尺度数值模式,采用标准F90程序进行编写,具有更加丰富的参数化方案及精细的物理过程,在中尺度天气系统中的风场、高度场、垂直速度场等方面的模拟值精确度较高[24],对我国的区域模拟效果同样较好[25],因此本文选择WRF模式模拟气象因子。WRF模式采用的是FNL(final operational global analysis)全球分析资料数据,FNL数据由二进制格点方式生成,其同化面较广,在气象领域研究较为广泛,能够为气象、气候模拟提供初始场,对极端天气也能够取得较好的研究结果[26]。其空间分辨率为1°×1°,时间分辨率为6 h,无法满足与高时空分辨率的气溶胶数据的匹配,因此通过WRF模型对FNL数据进行降尺度处理(具体参数如表1所示),空间上输出3 km×3 km分辨率,时间上以每小时整点输出,提取的气象要素为1 000 hPa(视为近地面),850 hPa,700 hPa与500 hPa气压面的温度(T)、平均风速(WS)和相对湿度(RH),以及边界层高度(PBLH),经裁剪、提取、计算均值等处理后,获取南京市范围内气象因子数据集。
表1 WRF模式部分参数方案
Tab.1
参数 | 采用方案 |
---|---|
地图投影 | Lambert投影 |
垂直分层嵌套层数 | 14层 |
嵌套层数 | 3层 |
微物理过程 | Lin 方案 |
积云参数化 | Kain-Fritsch 方案 |
边界层参数化 | Mellor-Yamada-Janjic湍流动能方案 |
陆面过程 | Noah 方案 |
短波辐射 | Dudhia 方案 |
长波辐射 | RRTM方案 |
2 目标方法与流程
以卫星过境时刻为准,生成同时刻(或时间差控制在30 min内)的AOD,FMF,PM2.5和模拟气象因子数据集,并按季节分组。首先分季节对PM2.5和AODm,PM2.5和AODf进行线性回归分析,计算仅采用AOD作为输入参数时大气颗粒物的模型拟合效果,并讨论精细模式分数对PM2.5质量浓度反演的影响; 其次利用随机森林模型,分别通过AODf与近地面气象因子(1 000 hPa气压层处的温度、平均风速和相对湿度,以及边界层高度)、AODf与近地面及不同高度气象因子(1 000 hPa,850 hPa,700 hPa,500 hPa气压面的平均风速和相对湿度,以及近地面温度和边界层高度),构建PM2.5质量浓度的反演模型,探讨含高空气象因子对PM2.5质量浓度反演的影响。其中随机选取2/3的数据进行建模,剩余1/3的数据则用于验证。采用决定系数R2,均方根误差(root mean square error,RMSE)、平均绝对误差(mean absolute error,MAE)与平均绝对百分比误差(mean absolute percentage error,MAPE),共4 个指标为依据,对两个模型进行评估,并利用随机森林模型获取不同因子的评分,分析其对PM2.5质量浓度的影响。
3 结果
3.1 FMF对AOD估算PM2.5的影响分析
图1
图1
南京市AOD和PM2.5质量浓度的年均、季均值
Fig.1
Annual and quarterly average concentrations of AOD and PM2.5. in Nanjing
图2与图3分别显示了4个季节里通过垂直与湿度订正AODm、经过垂直、湿度订正与精细模式分数校正后得到的AODf与PM2.5质量浓度的拟合关系。经对比可发现经过FMF校正后的AODf与PM2.5拟合系数均有提高,RMSE也更小(均通过显著性检验)。春季的拟合结果提升比较明显,因春季中国北方沙尘天气活跃,在偏北风影响下向东南部地区扩散,大气中粗颗粒占比较高。通过FMF进行校正能够有效减少粗颗粒对分析结果的影响,R2由0.23提升至0.33。夏季雨水充沛,对大气中粗颗粒物的冲刷能力较强,且高温高湿度环境有利于污染气体向大气颗粒物的转化作用[27],夏季大气中主要为细颗粒物,因此夏季R2达到0.37和0.40,效果好于其他季节,但FMF校正前后提升却并不明显。秋季和冬季由于大气环境趋于平稳,偶有沙尘输入现象,且由于入冬后受到北方化石燃料燃烧增加的影响,总体受外源污染物输入增多,相关性表现一般。此外由于受到PM2.5质量浓度极高值的影响,RMSE明显大于春季和夏季,冬季出现重度污染天气多于其他季节,即使经过精细模式分数校正,RMSE依然高达37.66 μg/m 3。估算PM2.5的RMSE随PM2.5均值的提高而上升,该趋势与Tian等[28]利用MODIS估算PM2.5得出的结论一致。经过精细模式分数校正,AOD与PM2.5的相关性分析更优,但由于卫星获取的气溶胶数据受到云的干扰较强,大部分情况下并不能完全覆盖研究区,虽然已将只有零散分布的图像进行剔除,但是研究区范围内有效数据的缺失在一定程度上增加了精细模式分数对PM2.5质量浓度反演的不确定性。
图2
图3
3.2 基于WRF模式的多元回归分析
表2 不同季节的多元回归模型
Tab.2
仅含地面气象因子 | 包含不同高度气象因子 | |||||||
---|---|---|---|---|---|---|---|---|
季节 | R2 | RMSE | MAE | MAPE/% | R2 | RMSE | MAE | MAPE/% |
春季 | 0.57 | 14.97 | 15.56 | 34.03 | 0.71 | 10.03 | 10.98 | 28.35 |
夏季 | 0.71 | 11.45 | 13.50 | 34.85 | 0.74 | 6.83 | 7.68 | 23.19 |
秋季 | 0.62 | 18.77 | 18.52 | 39.07 | 0.71 | 13.27 | 14.15 | 38.92 |
冬季 | 0.65 | 20.59 | 25.51 | 41.35 | 0.73 | 19.92 | 20.12 | 37.22 |
图4
图4
仅含近地面气象因子模型的验证结果
Fig.4
Verification results of model with only near-ground meteorological factors
图5-1
图5-1
包含不同高度气象因子的验证结果
Fig.5-1
Verification results of model including meteorological factors on different altitudes
图5-2
图5-2
包含不同高度气象因子的验证结果
Fig.5-2
Verification results of model including meteorological factors on different altitudes
将WRF模拟气象因子加入回归分析中,各季节都取得较好的结果,说明气象因子对PM2.5的质量浓度有一定的影响。图6中,RH1000,RH850,RH700,RH500分别代表1 000 hPa,850 hPa,700 hPa和500 hPa气压面处的相对湿度(%),WS1000,WS850,WS700,WS500分别代表1 000 hPa,850 hPa,700 hPa,500 hPa气压面处的风速(m/s),T2代表距地面2 m处的温度(℃)。从季节上来看,春季提升效果较大,R2提高到0.71; 夏季提升效果虽然较小,但R2依旧最高,达到0.74,且各项误差同样为各季节最小。冬季的R2为0.73,虽然高于春季与秋季,但由于受到重污染天气期间PM2.5质量浓度极高值的影响,导致各项误差明显大于其他季节,RMSE与MAE分别达到19.92 μg/m 3与20.12 μg/m 3。从模型上来看,包含不同高度气象因子的模型的各项指标均优于仅含地面气象因子的模型,说明高空气象因子对PM2.5质量浓度拥有一定程度的影响,且对于高值估算方面,估算值与实测值更加接近。
图6
图6
四季不同影响因子对PM2.5质量浓度的影响
Fig.6
Effect of different influence factors on PM2.5 mass concentration in four seasons
利用随机森林模型拟合PM2.5质量浓度与各影响因子的关系,获得各影响因子的评分(图6),各因子在不同季节上的重要性存在一定差异。AOD对PM2.5质量浓度的影响程度较高,其中在夏季表现最为明显。其次是边界层高度,代表了大气在垂直方向上的扩散能力,有利于大气颗粒物的扩散,各季节得分均比较稳定。除夏季外,近地面和高空风速的影响程度总体高于湿度,水平扩散能力对PM2.5质量浓度的影响总体大于在湿度影响下气体向颗粒物转化作用的影响。而夏季处于一年中温度最高的时期,在高空中的高温高湿度环境会促进气体向颗粒物的转化作用,且夏季盛行东南季风,受北方沙尘输入影响较小,近地面风速的影响力明显高于高空。夏季近地面气象因素影响力高于高空气象因子影响力,与将高空气象因子加入随机森林模型进行回归,R2提升较小的趋势相对应。
4 结论与展望
1)2014—2019年间南京市年均PM2.5质量浓度总体呈逐年下降的趋势,其中夏季PM2.5质量浓度最低,为40.53 μg/m 3,冬季则最高,达到71.16 μg/m 3。南京市年均AOD值则稳定在0.7左右,季均值则呈现由春季向冬季逐渐较小的趋势,且春季夏季均值明显高于秋季与冬季。
2)通过利用MODIS获取的FMF对气溶胶光学厚度AOD进行校正,能够有效提高AODf与PM2.5质量浓度的相关性,且均方根误差也均有降低。春季由于受到北方粗颗粒沙尘输入影响,因此经过校正后能够有效减少大气粗颗粒对PM2.5质量浓度反演的影响,因此效果提升最为明显。夏季AODf与PM2.5质量浓度的相关性最高,RMSE最小。而秋季冬季因受到重度污染天气的影响,均方根误差相对较高。但由于卫星受云的干扰性较大,尤其无法获取发生雨水天气时的有效数据,且有效值并不能完全覆盖整个研究区,局部范围的有效数据作为整个研究区的平均数据,可能存在一定程度的不确定性。
3)通过引入含不同高度的气象因子能够有效提高利用AOD反演PM2.5质量浓度的相关性,且包含不同高度气象因子的模型的各项指标均优于仅含地面气象因子的模型。在各气象因子中,各季节风速和边界层高度,对PM2.5质量浓度影响程度较高,风速和湿度在不同季节对PM2.5质量浓度的影响有所不同,除夏季外,风速影响程度高于湿度; 而在夏季,近地面气象因子影响程度则明显高于高空。在引入高空气象因子分析后,秋季和冬季R2虽有提升,但均方根误差依然较高。因此在未来的研究中,针对重度污染天气的大气颗粒物估算,则需要进一步考虑更多大气环境因素、污染源排放等影响。
参考文献
2013 Southeast Asian smoke haze:Fractionation of particulate-bound elements and associated health risk
[J].DOI:10.1021/es405533d URL [本文引用: 1]
Emergency admissions for cardiovascular and respiratory diseases and the chemical composition of fine particle air pollution
[J].DOI:10.1289/ehp.0800185 URL [本文引用: 1]
Long and short term exposure to PM2.5 and mortality
[J].
DOI:10.1097/EDE.0b013e318294beaa
PMID:23676266
[本文引用: 1]
Many studies have reported associations between ambient particulate matter (PM) and adverse health effects, focused on either short-term (acute) or long-term (chronic) PM exposures. For chronic effects, the studied cohorts have rarely been representative of the population. We present a novel exposure model combining satellite aerosol optical depth and land-use data to investigate both the long- and short-term effects of PM2.5 exposures on population mortality in Massachusetts, United States, for the years 2000-2008. All deaths were geocoded. We performed two separate analyses: a time-series analysis (for short-term exposure) where counts in each geographic grid cell were regressed against cell-specific short-term PM2.5 exposure, temperature, socioeconomic data, lung cancer rates (as a surrogate for smoking), and a spline of time (to control for season and trends). In addition, for long-term exposure, we performed a relative incidence analysis using two long-term exposure metrics: regional 10 × 10 km PM2.5 predictions and local deviations from the cell average based on land use within 50 m of the residence. We tested whether these predicted the proportion of deaths from PM-related causes (cardiovascular and respiratory diseases). For short-term exposure, we found that for every 10-µg/m increase in PM 2.5 exposure there was a 2.8% increase in PM-related mortality (95% confidence interval [CI] = 2.0-3.5). For the long-term exposure at the grid cell level, we found an odds ratio (OR) for every 10-µg/m increase in long-term PM2.5 exposure of 1.6 (CI = 1.5-1.8) for particle-related diseases. Local PM2.5 had an OR of 1.4 (CI = 1.3-1.5), which was independent of and additive to the grid cell effect. We have developed a novel PM2.5 exposure model based on remote sensing data to assess both short- and long-term human exposures. Our approach allows us to gain spatial resolution in acute effects and an assessment of long-term effects in the entire population rather than a selective sample from urban locations.
Remote sensing of tropospheric aerosols from space:Past,present and future
[J].DOI:10.1175/1520-0477(1999)080<2229:RSOTAF>2.0.CO;2 URL [本文引用: 1]
Dynamic issues of urban atmospheric pollution models
[J].
Influence domain and climate effect related to aerosol of urban community around Beijing
[J].
Community-level spatial heterogeneity of chemical constituent levels of fine particulates and implications for epidemiological research
[J].
Integrated aerosol optical thickness,gaseous pollutants and meteorological parameters to estimate ground PM2.5 concentration
[J].
Characteristics of aerosol observed during two severe haze events over Korea in June and October 2004
[J].DOI:10.1016/j.atmosenv.2006.03.050 URL [本文引用: 1]
Sptio-temporal variability of satellite derived aerosol optical thickness over Northeast Asia in 2004
[J].DOI:10.1016/j.atmosenv.2007.01.048 URL [本文引用: 1]
Intercomparison between satellite-derived aerosol optical thickness and PM2.5 mass:Implications for air quality studies
[J].DOI:10.1029/2003GL018174 URL [本文引用: 1]
Qualitative and quantitative evaluation of MODIS satellite sensor data for regional and urban scale air quality
[J].DOI:10.1016/j.atmosenv.2004.01.039 URL [本文引用: 1]
Correlating MODIS aerosol optical thickness data with ground-based PM2.5 observations across Texas for use in a real-time air quality prediction system
[J].DOI:10.1016/j.atmosenv.2005.08.036 URL [本文引用: 1]
Estimating ground-level PM2.5 in the eastern United States using satellite remote sensing
[J].DOI:10.1021/es049352m URL [本文引用: 1]
Montreux,Switzerland,
Satellite-based PM2.5 estimation using fine-mode aerosol optical thickness over China
[J].DOI:10.1016/j.atmosenv.2017.09.023 URL [本文引用: 1]
大气气溶胶粒径分布特征与气象条件的相关性分析
[J].
Relationship between meteorological conditions and particle size distribution of atmospheric aerosols
[J].
基于MODIS数据的环渤海地区气溶胶时空变化特征分析
[J].
Temporal and spatial variation characteristics analysis of aerosol around Bohai Sea region based on MODIS data
[J].
Spatio-temporal pattern estimation of PM 2.5 in Beijing-Tianjin-Hebei region based on MODIS AOD and meteorological data using the back propagation neural network
[J].
融机器学习与WRF大气模式的PM2.5预报方法
[J].
PM2.5 forecasting method based on machine learning and WRF hybrid model
[J].
基于MODIS的城市大气颗粒物污染指数研究
[J].
Particulate matter indices derived from MODIS data for indicating urban air pollution
[J].
MODIS气溶胶产品在北京监测PM2.5质量浓度中的应用
[J].
Application of MODIS aerosol products in monitoring PM2.5 concentration in Beijing
[J].
Remote sensing of atmospheric fine particulate matter (PM2.5) mass concentration near the ground from satellite observation
[J].DOI:10.1016/j.rse.2015.02.005 URL [本文引用: 2]
Precipitation simulations using WRF as a nested regional climate model
[J].DOI:10.1175/2009JAMC2186.1 URL [本文引用: 1]
WRF模式对2008年1月我国南方冻雨极端天气过程的数值模拟研究
[D].
WRF numerical simulation of freezing rain extreme weather process in southern China in January 2008
[D].
WRF模式气象数据在遥感反演PM2.5中的应用研究
[J].
Application of WRF model meteorological data in PM2.5 remote sensing inversion
[J].
相对湿度对大气气溶胶消光系数的影响
[J].
The effect of relative humidity on the extinction coefficient of aerosols
[J].
A semi-empirical model for predicting hourly ground-level fine particulate matter (PM2.5) concentration in southern Ontario from satellite remote sensing and ground-based meteorological measurements
[J].DOI:10.1016/j.rse.2009.09.011 URL [本文引用: 1]
/
〈 |
|
〉 |
