ASSESSMENTS OF FORECAST PRODUCTS OF AIR QUALITY FORECASTING SYSTEM BASED ON GRAPES-CMAQ IN ZHONGSHAN CITY
-
摘要: 对基于GRAPES-CMAQ模式的广东省中山市空气质量预报系统2016年的AQI等级、首要污染物、PM2.5和O3浓度预报数据进行检验评估。预报系统对中山市全年的AQI等级24小时预报TS评分为0.66, TS评分在6—7月最高。当实况AQI等级为一级, TS为0.84, 当空气质量变差时, TS明显下降, 说明模式目前对空气重污染事件的预报能力不足。模式对中山市的主要污染物O3、PM和NO2的预报TS评分分别为0.21、0.24和0.24, 其中对O3和NO2的漏报率达到0.78和0.76, 而PM的空报率高达0.71。预报系统对PM2.5的浓度预报总体略偏高, 相关系数为0.59, 但在霾污染较严重的月份中预报效果并不稳定。对O3的预报全年相关系数达到0.7, 且在O3污染较严重的7—11月, 相关系数均能达到0.5以上, 但这几个月的O3预报浓度均比实况明显偏低。从气象预报入手分析预报系统对10月2—7日中山市一次连续O3污染过程漏报的原因, 发现GRAPES模式对这段时间内整个珠三角地区日间降水的预报范围偏大, 强度偏强, 且预报的地面风速也比实况要偏大, 导致模拟O3浓度比实况明显偏低。
-
关键词:
- 环境气象 /
- 空气质量 /
- GRAPES-CMAQ /
- 评估 /
- O3
Abstract: In this study, forecast products for 2016, including AQI level, primary pollutant, and PM2.5 and O3 concentrations, from an air quality forecasting system based on GRAPES-CMAQ model in Zhongshan city are assessed with the observations. The annual average TS score of forecasted AQI level was up to 0.66 and the best scores were found in June and July. When the observed AQI level was the first-class, the TS score was up to 0.84, but the TS score dropped evidently when the air pollution got worse, indicating that the model still could not give good performance in severe air pollution events. The annual average TS scores for the forecast of primary pollutants of O3, PM and NO2 were 0.21, 0.24 and 0.24, respectively. The PO for O3 and NO2 were up to 0.78 and 0.76, respectively, while the FAR for PM was 0.71. The annual average forecasted PM2.5 concentration was a little overestimated compared to the observations and the correlation coefficient was 0.59 for the whole year. However, the performances of the model were not stable among the months with heavy haze pollution. The correlation coefficient for O3 was 0.7 in 2016, and the coefficient of July to November with severe O3 pollution were all over 0.5. The forecasted O3 concentrations were all underestimated in these five months. An O3 pollution case during 2nd to 7th October 2016 was studied based on the forecasted data of GRAPES. It was found that the rainfall area and intensity in Pearl River Delta as well as the wind speed were all overestimated by GRPAES in the day time, leading to the underestimation of forecasted O3 concentrations in this event.-
Key words:
- environmental meteorology /
- air quality /
- GRAPES-CMAQ /
- assessment /
- O3
-
表 1 旱涝等级指标
GRAPES模式 CMAQ模式 嵌套区域 网格数 网格距/km 中心经纬度 嵌套区域 网格数 网格距/km 中心经纬度 1 283×184 27 121.86°E, 23.05°N 1 182×138 27 121.86°E, 23.05°N 2 233×163 9 117.55°E, 23.02°N 2 98×74 9 117.55°E, 23.02°N 3 172×130 3 113.38°E, 22.93°N 3 152×110 3 113.38°E, 22.93°N 表 2 空气质量预报系统对中山市2016年全年及各月份PM2.5浓度预报偏差
时间 平均误差 平均绝对误差 均方根误差 相关系数 2016全年 0.8 14 20 0.59 1月 -7 17 23 0.28 2月 10 20 26 0.28 3月 -4 20 26 0.33 4月 -12 16 20 0.50 5月 -4 14 18 0.54 6月 0 6 10 0.11 7月 -1 5 8 0.81 8月 1 11 16 0.57 9月 7 12 18 0.69 10月 2 10 16 0.50 11月 10 17 23 0.52 12月 9 19 27 0.54 表 3 空气质量预报系统对中山市2016年全年及各月份O3浓度预报偏差
时间 平均误差 平均绝对误差 均方根误差 相关系数 2016全年 -9 25 36 0.70 1月 -5 21 30 0.15 2月 -5 23 30 0.31 3月 4 26 36 0.64 4月 -13 21 31 0.66 5月 -15 24 32 0.77 6月 -7 19 25 0.57 7月 -4 20 29 0.83 8月 -15 33 46 0.69 9月 -14 29 39 0.83 10月 -18 35 49 0.53 11月 _6 25 35 0.72 12月 -11 28 37 0.62 表 4 实况和预报的中山市10月2—7日08—20时的平均温度、平均相对温度、平均降水量、平均风速以及主导风向
物理量 温度/℃ 相对湿度/% 降水量/mm 地面风速/(m/s) 地面主导风向 实况 28.7 72 0 1.9 N 预报 27.6 84 5.8 3.5 NE -
[1] 王自发, 谢付莹, 王喜全, 等.嵌套网格空气质量预报模式系统的发展与应用[J].大气科学, 2006, 30(5):778-790. [2] 邓涛, 吴兑, 邓雪娇, 等.珠江三角洲一次典型复合型污染过程的模拟研究[J].中国环境科学, 2012, 32(2):193-199. [3] 王自发, 李杰, 王哲, 等. 2013年1月我国中东部强霾污染的数值模拟和防控对策[J].中国科学:地球科学, 2014, 44(1):3-14. [4] CHE W, ZHENG J, ZHONG L, et al. Assessment of motor vehicle emission control policies using Model-3/CMAQ model for the Pearl River Delta region, China[J]. Atmos Env, 2011, 45(9):1 740-1 751. [5] ZHANG R, SARWAR G, FUNG J C H, et al. Role of photoexcited nitrogen dioxide chemistry on ozone formation and emission control strategy over the Pearl River Delta, China[J]. Atmos Res, 2013, 132-133(10):332-344. [6] 王茜, 吴剑斌, 林燕芬. CMAQ模式及其修正技术在上海市PM2.5预报中的应用检验[J].环境科学学报, 2015, 35(6):1 651-1 656. [7] 陈彬彬, 林长城, 杨凯, 等.基于CMAQ模式产品的福州市空气质量预报系统[J].中国环境科学, 2012, 32(10):1 744-1 752. [8] 程念亮, 李红霞, 孟凡, 等.山东省空气质量预报平台设计及其预报效果评估[J].环境污染与防治, 2015, 37(9):92-99. [9] 王茜, 伏晴艳, 王自发, 等.集合数值预报系统在上海市空气质量预测预报中的应用研究[J].气象, 2010, 2(4):1-6. [10] 邓涛, 邓雪娇, 吴兑, 等.珠三角灰霾数值预报模式与业务运行评估[J].气象科技进展, 2012, 2(6):38-44. [11] 谢敏, 钟流举, 陈焕盛, 等. CMAQ模式及其修正预报在珠三角区域的应用检验[J].环境科学与技术, 2012, 35(2):96-101. [12] MAI J H, DENG T, YU L L, et al. A modeling study of impact of emission control strategies on PM2. 5 reductions in Zhongshan, China, using WRF-CMAQ[J]. Advances in Meteorology, 2016, Article ID 5836070. [13] 麦健华, 邓涛, 于玲玲, 等.中山市旱季霾特征及数值模拟分析[J].环境科学学报, 2016, 36(6):2 170-2 179. [14] 李颖敏, 邓雪娇, 邓涛, 等.中山市2000-2014年霾天气特征及气象影响因子[J].中国环境科学, 2016, 36(6):1 638-1 644. [15] 陈德辉, 薛纪善, 沈学顺, 等.我国自主研制的全球/区域一体化数值天气预报系统GRAPES的应用与展望[J].中国工程科学, 2012, 14(9):46-54. [16] 钟水新, 陈子通, 黄燕燕, 等.地形重力波拖曳参数化方案在华南中尺度模式(GRAPES)中的应用试验[J].热带气象学报, 2014, 30(3):413-422. [17] 陈子通, 万齐林, 沈学顺, 等. GRAPES区域模式水汽平流方案的比较与改进[J].热带气象学报, 2010, 26(1):1-6. [18] 王占山, 李晓倩, 王宗爽, 等.空气质量模型CMAQ的国内外研究现状[J].环境科学与技术, 2013(S1):386-391. [19] HE K B. Multi-resolution Emission Inventory for China(MEIC): model framework and 1990-30 2010 anthropogenic emissions[R]. International Global Atmospheric Chemistry Conference. [20] 郑君瑜, 张礼俊, 钟流举, 等.珠江三角洲大气面源排放清单及空间分布特征[J].中国环境科学, 2009, 29(5):455-460. [21] 白爱民. AQI vs API——新老空气质量标准之对比[J].环境科学导刊, 2013, 32(6):95-97. [22] 中华人民共和国环境?げ?环境空气质量指数(AQI)技术规定(试行)[S]. HJ 633-2012. [23] 廖志恒, 孙家仁, 范绍佳, 等. 2006-2012年珠三角地区空气污染变化特征及影响因素[J].中国环境科学, 2015, 35(2):329-336. [24] LI Y C, YU J Z, HO S S H, et al. Chemical characteristics of PM2.5 and organic aerosol source analysis during cold front episodes in Hong Kong, China[J]. Atmos Res, 2012, 118(3):41-51. [25] LI J, LU K, LV W, et al. Fast increasing of surface ozone concentrations in Pearl River Delta characterized by a regional air quality monitoring network during 2006-2011[J]. J Env Sci, 2014, 14(1):23-26. [26] 刘健, 吴兑, 范绍佳, 等.前体物与气象因子对珠江三角洲臭氧污染的影响[J].中国环境科学, 2017, 37(3):813-820. [27] 张艳霞, 蒙伟光, 戴光丰, 等.城市冠层模式在GRAPES模式中的应用[J].热带气象学报, 2016, 32(3):311-321.