基于气候系统指数的月尺度霾日定量预测方法研究

付伟 司红君 何冬燕 卢尧 刘蕾 杨琼琼 邹莹瑾

摘要 利用1955—2018年芜湖市国家气象观测站资料，1980—2018年国家气候中心气候系统指数资料，对芜湖市近64 a的霾日气候序列进行重建，在此基础上使用线性趋势和Mann-Kendall方法，系统分析了芜湖市霾日的气候特征。以芜湖市为例，借助多元逐步回归方法，尝试研究了一种以气候系统指数为自变量，霾日为因变量，建立霾日预测方程并对月尺度霾日进行定量预测的方法。结果表明：1）重建的霾日时间序列能够更客观的反映芜湖市霾日实际的长期变化特征。近64 a霾的气候特征：年日数显著增加，并在1980年左右出现突变;季日数在冬、秋季较多，夏季最少，四季均呈显著增多趋势，增速秋季最快，夏季最慢;月日数在1、12、11月较多，7、8月较少。中度以上霾的气候特征：年日数显著增加，表现出渐变特征;季日数冬季最多，夏季最少，除夏季外均显著增多趋势，增速冬季最快，秋、春季次之;月日数在1、12、11月较多，5、6月较少，8月则从未出现过。2）各月霾日、中度以上霾日与多项月气候系统指数表现出显著的相关性，使用这些指数，计算出各月霾日、中度以上霾日的预测方程，最后对方程效果进行检验。

关键词 霾日;重建;气候系统指数;多元逐步回归;月尺度定量预测

霾是大量极细微的干尘粒等颗粒均匀浮游在空中，使水平能见度小于10 km的空气普遍浑浊现象（中国气象局，2010）。过去，形成霾的干气溶胶粒子主要来自自然界，但随着人类活动加剧，经济快速发展，城镇化率提高，近二三十年中国中东部地区人为排放的大气气溶胶显著增加（张小曳等，2013），这其中包括大量的细粒子（PM2.5），造成霾天气日益增多。近年来，长三角地区霾日在时间和空间上均呈增多趋势（史军和崔林丽，2013）。安徽霾日在20世纪80年代初之后发生激增，东部霾日数明显高于西部（邓学良等，2015），这些可观测到的霾污染增多，严重危及人类健康和交通安全等（Chameides et al.，1999）。

近年来科技工作者对霾的识别、成分特征、污染影响机制、气象影响条件等进行了大量的研究（吴兑，2008;吴丹等，2016;Yin and Wang，2016，2017）。霾天气时，空气中的硫酸盐气溶胶是造成能见度下降的主要原因（Malm，1992），有机物、水溶性无机离子和元素碳是PM2.5中影响大气能见度，影响大气消光的主要化学组分（Cao et al.，2012）。天气也是影响霾污染发生发展的重要因素，有利的天气形势可以造成稳定的大气层结，混合层降低，形成近地面逆温，地面风速小，大气垂直方向湍流交换减弱的气象条件，造成霾天气的维持和加剧（周莉等，2019），而冷空气南下造成污染物伴随输入至长三角地区也是造成重污染发生的有利天气条件，这种污染往往发生在重污染过程的初期（韩博威和马晓燕，2020）。为了解决霾造成的一系列社会经济问题，对霾进行了解和控制，珠三角、长三角、京津冀等地建立了霾的地面和遥感监测体系（向嘉敏等，2019），并开展霾预报预测研究和业务，主要预报方法包括逐步回归（尹志聪等，2019），动态统计（张自银等，2018），基于模式的MOS法（尤佳红等，2014）等，但是中长期预测尤其是定量预测的方法相对较少。大气环流、海温等气候系统指数常被用于气象要素的气候预测（陈丽娟等，2019），以及其他与气候密切相关事件的预测（王纯枝等，2019）。目前基于气候系统指数的霾日中长期定量预测方法并不多见，而霾和天气气候密切相关（陶丽等，2016），因此本文拟研究使用气候系统指数作为因子定量预测月尺度霾日的方法，并以蕪湖市为例开展。

芜湖市位于安徽省东部，霾的气候特征及预测对其经济发展、人民生活、长江生态环境保护有积极意义，但与其相关的研究成果大都时间相对较早，针对性不强。本研究首先对芜湖市霾的气候特征进行分析，再建立芜湖市月尺度霾日的预测方程，为其中长期空气质量预测业务提供客观定量的参考。但长期以来中国气象台站的观测业务中，对于霾的记录没有统一标准，尤其是作为霾辅助判据的相对湿度取值范围过宽，标准不一，大量霾被记成了轻雾（吴兑，2006），同时也受观测员主观因素影响，造成霾的观测记录可对比性较差，因此用长时间序列气候资料进行霾统计时，不能直接使用天气现象观测数据，要建立统一的定量标准（吴兑，2012），需要先对霾日的时间序列进行重建，再进行其他分析。因此本文对芜湖市霾日的时间序列进行重建（吴兑，2008），然后再进行气候特征分析和预测方程的研究。

1 资料与方法

1.1 资料来源

本文使用1955—2018年芜湖市国家气象观测站的观测资料重建霾日时间序列及气候特征分析，包括逐日能见度、相对湿度、天气现象数据。观测数据1955—1985年为每日02、08、14、20时4次观测，1986—2014年8月为每日08、14、20时3次观测，2014年9月—2018年12月为每日24次观测。其中能见度数据1955—1979年为等级观测，其他为距离观测。

用于霾日预测方程计算的130项气候系统指数为1980—2018年逐月的88项大气环流指数、26项海温指数、16项其他指数。数据来自国家气候中心网站百项气候系统指数集（https：//cmdp.ncc-cma.net/Monitoring/cn_index_130.php）。

1.2 霾日重建方法

1.2.1 能见度、相对湿度数据的处理

日平均值计算时，每日3次或4次观测的，取其平均值;每日24次观测的，取02、08、14、20时观测平均值。能见度等级观测数据，需将其按照表1换算为对应的距离中值（史军和吴蔚，2010）。另外，2014年9月之后能见度观测改为自动观测，需将之前人工观测数据乘以0.75的系数进行订正，以确保数据具有可对比性（中国气象局气象探测中心，2016）。

1.2.2 霾日的统计

采用日均法对霾日进行统计（史军和吴蔚，2010;吴兑等，2010），定义日平均能见度<10 km，日平均相对湿度<90 %，并排除降水、吹雪等其他天气现象导致低能见度事件的情况，记为1个霾日。参照《霾的观测和预报等级》（QX/T 113-2010），根据能见度距离（V，单位：km）将霾划分为轻微、轻度、中度和重度4个等级，划分标准依次为：5≤V<10、3≤V<5、2≤V<3、V<2。

1.3 霾日的气候特征分析方法

使用线性趋势分析和Mann-Kendall气候突变检验（简称M-K检验，魏凤英，2007）对芜湖市霾日的气候特征进行分析。季节分析时，划分如下：3—5月春季，6—8月夏季，9—11月秋季，12—次年2月冬季。

1.4 霾日的月尺度定量预测方程计算方法

使用多元逐步回归方法计算月尺度霾日的预测方程。计算方程前，首先使用月霾日时间序列与对应的前1—12月130项气候系统指数时间序列分别进行相关分析，并进行显著性检验，遴选出显著相关的气候系统指数。计算方程时，应变量为月霾日时间序列，自变量为遴选出的所有与其显著相关的月气候系统指数时间序列。计算出方程后，方程和系数进行相关性和显著性检验。各月霾日依次计算。由于20世纪70年代海气关系出现转变（江志红和屠其璞，2001），为保证建立的方程有一定气候意义，因此建立方程时仅使用1980年之后的数据。

2 结果分析

2.1 霾日的重建及气候特征

2.1.1 霾日重建值和实测值对比

对比1955—2018年芜湖市霾日的重建值与芜湖市国家观测站人工记录的霾天气现象实测值（图1a），1978年之前霾日较少，重建值和实测值差别不大。1978年之后重建值明显增加，而实测值仍变化不大，两者差距明显。对比芜湖市所在区域霾日变化特征的相关研究结果（史军和崔林丽，2013;邓学良等，2015），认为重建的霾日时间序列更能够客观反映霾日的实际变化特征，许多学者也通过研究证明了重建的可行性和必要性（史军和吴蔚，2010;杨志彪等，2014;王凯等，2016）。因此下文使用重建值进一步进行研究。

2.1.2 霾日的年际变化特征

1955—2018年，芜湖市年霾日总体呈显著增加趋势（图1），年均增加霾日2.08 d，1978年之前霾日较少，随后开始快速增长，其中2015年霾日最多，达191 d。2013年中国出台《大气污染防治行动计划》（简称“大气十条”）后，京津冀、长三角等重点区域不断加大减排措施力度，取得了实质性的进展，2016年起冬季污染物的排放量得到了有效消減，这是2015年之后空气质量得到改善的主要原因（张小曳等，2020）。另外，“大气十条”实施后的2014、2015年长三角等重点区域的气象条件较2013年变差，2016、2017年转好，也在一定程度上促成2015年霾日数明显高于相邻年份的特征。M-K检验表明芜湖市1980年之前，霾日变化不显著，增多减少交替出现，但1980年之后显著增多，存在突变特征（图2a）。

1955—2018年芜湖市轻微、轻度、中度和重度霾日数分别占总霾日数的88.3 %、9.1 %、2.1 %和0.5 %，各等级霾日均呈显著增加趋势（图1b），分别年均增加1.704、0.281、0.078、0.017 d。中度以上的霾日虽占比小，但影响大，因此除了总霾日，本文还对中度和重度霾日的气候特征和预测方程进行分析和计算，又因重度霾次数少，故将二者合称为中度以上霾日统一进行分析。

1980年之前，芜湖市无中度以上霾发生，之后零星出现，1991年起明显增多，年均4 d，其中2015年最多，达27 d，包括重度霾5 d。M-K检验显示其没有表现出突变特征（图2b）。

2.1.3 霾日的季、月变化特征

1955—2018年，霾日（图3a）和中度以上霾日（图3b）均表现出冬季最多，秋、春、夏季依次减少的特征，分别为23.890 d、17.328 d、14.359 d、9.266 d和1.093 d、0.422 d、0.141 d、0.047 d。各季节霾日均表现出显著增多的趋势（表2），增速秋季最快，夏季最慢。中度以上霾日在冬、秋、春季表现出显著增多的趋势，增速依次减小（表2），夏季则没有显著变化趋势。

一年之中霾日最多的是1月（图略），年均9.047 d，其次是12月（8.922 d）和11月（7.297 d）。4—9月霾日相对较少，其中7月最少（2.359 d）。中度以上霾日（图略）1月和12月最多，年均0.453 d，其次是11月（0.266 d）和2月（0.188 d）。4—9月很少出现中度以上霾，其中8月出现天数为0，5、6月也仅出现过1 d。

2.2 霾日的月尺度定量预测

2.2.1 预测方程的计算

月尺度霾日的预测，实际业务应用时常需提前1～2个月制作，如预测2020年11月的霾日数，9或10月就需要计算出预测结果，而当月气候系统指数的及时获取具有一定不确定性，为了业务应用时霾日计算数据完整，本文计算某月（T）霾日时，使用（T-3至T-14）这12个月的气候系统指数建立预测方程。含有缺测数据的气候系统指数将被剔除。

首先，使用月（T）霾日的37 a时间序列（1982年1月—2018年12月）分别与其对应的T-3至T-14月逐月130项气候系统指数（共1 560个）的37 a时间序列（1980年10月—2018年9月）求相关系数（R），并进行0.05的显著性检验，遴选显著相关的气候系统指数。

其次，月（T）霾日37 a时间序列作为应变量，遴选出的所有与其显著相关气候系统指数37 a时间序列作为自变量带入多元逐步回归模型，计算出该月（T）霾日的预测方程。计算出的方程和系数进行相关性和显著性检验。各月霾日数的预测方程依次计算，其中遴选出的作为自变量带入模型的显著相关气候系统指数数量不等。中度以上霾日的预测方程计算方法相同。

逐步回归是一种线性回归模型自变量选择方法，其基本思想是将变量一个一个引入，引入的条件是其偏回归平方和经检验是显著的。同时，每引入一个新变量后，对已入选回归模型的老变量逐个进行检验，将经检验认为不显著的变量删除，以保证所得自变量子集中每一个变量都是显著的。此过程经过若干步直到不能再引入新变量为止。这时回归模型中所有变量对因变量都是显著的（耿修林和张琳，2006），因此计算出的各最优方程中变量的个数是不同的。计算出的最优各月霾日、中度以上霾日的预测方程见表3、4，由于部分最优方程的自变量非常多，方程较长，受篇幅所限，这些方程使用10个变量的较优方程代替展示。所有方程均通过了0.001信度的显著性水平检验，其常数项、自变量系数均通过了0.05信度的显著性水平检验。由于8月没有中度以上霾日出现，因此没有预测方程。方程中气候系统指数使用字母加数字的方式表示，其中Atm代表大气环流指数（IAtm），Oce代表海温值数（IOce），Ext代表其他指数（IExt），其后数字代表指数的编号，括号中的数字代表指数的月份，各指数具体名称及定义可查阅国家气候中心网站百项气候系统指数集，如IAtm65（5）表示大气环流指数中的第65项，5月的西藏高原-1指数。

2.2.2 预测方程的检验

首先对各月霾日、中度以上霾日预测方程在1982—2018年的模拟效果进行检验，各月检验样本37个。月霾日误差范围在-7.2～7.4 d，68%的误差在2 d以内，其中2、6、9月误差较小，10、11、12月误差较大。月中度以上霾日误差范围在-1.6～1.8 d，明显优于霾日的模拟效果，其中4—7月的误差均在1 d以内，模拟效果很好，而3月相对较差。

接下来将2017年11月—2019年10月的逐月气候系统指数代入预测方程，对2019年1—12月芜湖市各月霾日、中度以上霾日的数量进行预测实验，并与实况对比检验（表5），多数月霾日的绝对误差在1～3 d，2、6、12月的预测误差相对较大。11月—次年3月霾日相对较多的月份，预测值普遍偏少。2019年实况仅1、2月出现中度以上霾，预测误差分别为0.7 d和-1.5 d，其余月份中有7个月预测值为负值，记为0，误差为0，预报基本准确。总体来说，预测方程可以较好地对月尺度霾日进行定量预测，尤其是月中度以上霾日，具有较好的参考价值。

3 结论与讨论

1）使用日均法对1955—2018年芜湖市霾日进行重建，对比霾日天气现象实测值，并参考其他相关研究结果，表明重建的霾日时间序列能够更客观地反映芜湖市霾日实际的长期变化特征。

2）1955—2018年蕪湖市年霾日呈显著增多趋势，年均增加2.080 d，在1980年左右出现突变特征。季霾日，年均数量从多到少依次为冬、秋、春、夏季，各季霾日均呈显著增多趋势，增速从快到慢依次为秋、冬、春、夏季。月霾日，1、12、11月较多，7、8月较少。

3）1955—2018年中度以上霾日也呈显著增多趋势，年均增加0.095 d，但没有出现突变特征。季中度以上霾日，数量从多到少依次为冬、秋、春、夏季，冬、秋、春季的中度以上霾日数表现出显著增多的趋势，增速依次减小，夏季则没有显著的变化趋势。月中度以上霾日数，1、12、11月较多，5、6月较少，8月从未出现过。

4）各月霾日、中度以上霾日与多项月气候系统指数表现出显著的相关性，使用这些显著相关的月气候系统指数作为自变量，月霾日、中度以上霾日作为因变量，用多元逐步回归方法计算出其月尺度预测方程。借助方程，可以提前查阅气候系统指数，对月霾日和中度以上霾日数进行计算、预测。预测方程进行了模拟效果检验、预测实验，结果显示预测方程具有较好的参考价值。

本研究表明，借助多元逐步回归方法，使用气候系统指数可以作为因子直接定量预测月尺度霾日，方法在气候相似的江淮等地区具有推广价值。这种统计方法通常历史拟合较好，但预报准确率不稳定，因此预测方程还需要进一步检验和完善，同时统计方法缺乏动力学意义（王冀等，2011），后期可以根据相关分析等结果，进一步分析具有动力学意义的气候系统指数，更深入地开展这些指数对霾天气影响机理的研究。

参考文献（References）

Cao J J，Wang Q Y，Chow J C，et al.，2012.Impacts of aerosol compositions on visibility impairment in Xi′an，China[J].Atmos Environ，59：559-566.doi：10.1016/j.atmosenv.2012.05.036.

Chameides W L，Yu H，Liu S C，et al.，1999.Case study of the effects of atmospheric aerosols and regional haze on agriculture：an opportunity to enhance crop yields in China through emission controls？[J].Proc Natl Acad Sci USA，96（24）：13626-13633.doi：10.1073/pnas.96.24.13626.

陈丽娟，顾薇，龚振淞，等，2019.影响2018年汛期气候的先兆信号及预测效果评估[J].气象，45（4）：553-564. Chen L J，Gu W，Gong Z S，et al.，2019.Precursory signals of the 2018 summer climate in China and evaluation of real-time prediction[J].Meteor Mon，45（4）：553-564.doi：10.7519/j.issn.1000-0526.2019.04.010.（in Chinese）.

邓学良，石春娥，姚晨，等，2015.安徽霾日重建和时空特征分析[J].高原气象，34（4）：1158-1166. Deng X L，Shi C E，Yao C，et al.，2015.Research of reconstruction and characteristics of hazes in Anhui[J].Plateau Meteor，34（4）：1158-1166.doi：10.7522/j.issn.1000-0534.2014.00007.（in Chinese）.

耿修林，张琳，2006.计量经济学[M].北京：科学出版社. Geng X L，Zhang L，2006.Econometrics[M].Beijing：Science Press.（in Chinese）.

韩博威，马晓燕，2020.2014—2018年冬季长三角强霾事件及天气形势影响分析[J].环境科学学报，40（7）：2333-2345. Han B W，Ma X Y，2020.Analysis of the severe haze events in the Yangtze River Delta during the winter of 2014—2018 and the impact of the weather situation on severe haze[J].Acta Sci Circumstantiae，40（7）：2333-2345.doi：10.13671/j.hjkxxb.2020.0075.（in Chinese）.

江志红，屠其璞，2001.国外有关海气系统年代际变率的机制研究[J].地球科学进展，16（4）：569-573. Jiang Z H，Tu Q P，2001.Review of overseas study on mechanisms of interdecadal variability of air-sea system[J].Adv Earth Sci，16（4）：569-573.doi：10.3321/j.issn：1001-8166.2001.04.020.（in Chinese）.

Malm W C，1992.Characteristics and origins of haze in the continental United States[J].Earth-Sci Rev，33（1）：1-36.doi：10.1016/0012-8252（92）90064-Z.

史军，吴蔚，2010.上海霾气候数据序列重建及其时空特征[J].长江流域资源与环境，19（9）：1029-1036. Shi J，Wu W，2010.Data reconstruction and spatial-temporal characteristics of haze days in Shanghai[J].Resour Environ Yangtze Basin，19（9）：1029-1036.（in Chinese）.

史军，崔林丽，2013.长江三角洲城市群霾的演变特征及影响因素研究[J].中国环境科学，33（12）：2113-2122. Shi J，Cui L L，2013.Characteristics and influencing factors of haze in the Yangtze River Delta region[J].China Environ Sci，33（12）：2113-2122.（in Chinese）.

陶丽，周宇桐，李瑞芬，2016.我国霾日和API分布特征及典型大城市中它们与气象条件关系[J].大气科学学报，39（1）：110-125. Tao L，Zhou Y T，Li R F，2016.Characteristics of haze days and air pollution index （API） and their relationship with weather conditions at four large cities/district of China[J].Trans Atmos Sci，39（1）：110-125.doi：10.13878/j.cnki.dqkxxb.20150721001.（in Chinese）.

王纯枝，张蕾，郭安红，等，2019.基于大气环流的稻纵卷叶螟气象预测模型[J].应用气象学报，30（5）：565-576. Wang C Z，Zhang L，Guo A H，et al.，2019.Long-term meteorological prediction model on the occurrence and development of rice leaf roller based on atmospheric circulation[J].J Appl Meteorol Sci，30（5）：565-576.doi：10.11898/1001-7313.20190505.（in Chinese）.

王冀，宋瑞艷，郭文利，2011.统计降尺度方法在北京月尺度预测中的应用[J].气象，37（6）：693-700. Wang J，Song R Y，Guo W L，2011.Application of statistical downscaling method to forecasts in monthly scale[J].Meteor Mon，37（6）：693-700.（in Chinese）.

王凯，方思达，李兰，2016.基于序列重建的湖北省大范围持续性霾过程特征分析[J].气象与环境学报，32（1）：40-45. Wang K，Fang S D，Li L，2016.Characteristics of persistent haze periods using large range sequence reconstruction in Hubei Province[J].J Meteorol Environ，32（1）：40-45.doi：10.11927/j.issn.1673-503X.2016.01.006.（in Chinese）.

魏凤英，2017.现代气候统计诊断与预测技术[M].北京：气象出版社. Wei F Y，2017. Modern climate statistical diagnosis and Prediction Technology[M].Beijing：China Meteorological Press.（in Chinese）.

吴丹，于亚鑫，夏俊荣，等，2016.南京市灰霾天气的长时间变化特征及其气候原因探讨[J].大气科学学报，39（2）：232-242. Wu D，Yu Y X，Xia J R，et al.，2016.Long-term variation in haze days and related climatic factors in Nanjing[J].Trans Atmos Sci，39（2）：232-242.doi：10.13878/j.cnki.dqkxxb.20150610001.（in Chinese）.

吴兑，2008.霾与雾的识别和资料分析处理[J].环境化学，27（3）：327-330. Wu D，2008.Discussion on the distinction between haze and fog and analysis and processing of data[J].Environ Chem，27（3）：327-330.doi：10.3321/j.issn：0254-6108.2008.03.011.（in Chinese）.

吴兑，2012.近十年中国灰霾天气研究综述[J].环境科学学报，32（2）：257-269. Wu D，2012.Hazy weather research in China in the last decade：a review[J].Acta Sci Circumstantiae，32（2）：257-269.doi：10.13671/j.hjkxxb.2012.02.011.（in Chinese）.

吴兑，吴晓京，李菲，等，2010.1951—2005年中国大陆霾的时空变化[J].气象学报，68（5）：680-688. Wu D，Wu X J，Li F，et al.，2010.Temporal and spatial variation of haze during 1951-2005 in Chinese mainland[J].Acta Meteorol Sin，68（5）：680-688.（in Chinese）.

向嘉敏，祝善友，张桂欣，等，2019.灰霾遥感监测研究进展[J].遥感技术与应用，34（1）：12-20. Xiang J M，Zhu S Y，Zhang G X，et al.，2019.Progress in haze monitoring by remote sensing technology[J].Remote Sens Technol Appl，34（1）：12-20.（in Chinese）.

杨志彪，陈冬冬，徐向明，2014.霾观测判识标准定量化对雾霾观测记录的影响[J].气象科技，42（2）：198-203. Yang Z B，Chen D D，Xu X M，2014.Influence of identifying index quantification on fog and haze observation records[J].Meteorol Sci Technol，42（2）：198-203.doi：10.19517/j.1671-6345.2014.02.003.（in Chinese）.

尤佳红，束炯，陈亦君，等，2014.基于MOS的杭州秋冬季空气污染预报和霾诊断[J].中国环境科学，34（7）：1660-1666. You J H，Shu J，Chen Y J，et al.，2014.Autumn and winter air pollution forecasting and haze events diagnosis in Hangzhou：a MOS-based approach[J].China Environ Sci，34（7）：1660-1666.（in Chinese）.

Yin Z C，Wang H J，2016.Seasonal prediction of winter haze days in the north central North China Plain[J].Atmos Chem Phys，16（23）：14843-14852.doi：10.5194/acp-16-14843-2016.

Yin Z C，Wang H J，2017.Statistical prediction of winter haze days in the North China plain using the generalized additive model[J].J Appl Meteorol Climatol，56（9）：2411-2419.doi：10.1175/jamc-d-17-0013.1.

尹志聰，王会军，段明铿，2019.近几年我国霾污染实时季节预测概要[J].大气科学学报，42（1）：2-13. Yin Z C，Wang H J，Duan M K，2019.Outline of the real-time seasonal haze pollution prediction in China in recent years[J].Trans Atmos Sci，42（1）：2-13.doi：10.13878/j.cnki.dqkxxb.20181224001.（in Chinese）.

张小曳，孙俊英，王亚强，等，2013.我国雾-霾成因及其治理的思考[J].科学通报，58（13）：1178-1187. Zhang X Y，Sun J Y，Wang Y Q，et al.，2013.Factors contributing to haze and fog in China[J].Chin Sci Bull，58（13）：1178-1187.doi：10.1360/972013-150.（in Chinese）.

张小曳，徐祥德，丁一汇，等，2020.2013～2017年气象条件变化对中国重点地区PM2.5质量浓度下降的影响[J].中国科学：地球科学，50（4）：483-500. Zhang X Y，Xun X D，Ding Y H，et al.，2020.The impact of meteorological changes from 2013 to 2017 on PM2.5 mass reduction in key regions in China [J].Sci China：Earth Sci，50（4）：483-500.（in Chinese）.

张自银，赵秀娟，熊亚军，等，2018.基于动态统计预报方法的京津冀雾霾中期预报试验[J].应用气象学报，29（1）：57-69. Zhang Z Y，Zhao X J，Xiong Y J，et al.，2018.The fog/haze medium-range forecast experiments based on dynamic statistic method[J].J Appl Meteorol Sci，29（1）：57-69.doi：10.11898/1001-7313.20180106.（in Chinese）.

中國气象局气象探测中心，2016.地面气象观测业务技术规定实用手册[M].北京：气象出版社. Meteorological Observation Center，China Meteorological Administration，2016.Standard business procedure for landserface meteorology observation[M].Beijing：China Meteorological Press.（in Chinese）.

中国气象局，2010.霾的观测和预报等级：QX/T 113-2010[S].北京：气象出版社. Observation and forecasting levels of haze：QX/T 113-2010[S].Beijing：China Meteorological Press.（in Chinese）

周莉，兰明才，张舰齐，等，2019.2017年冬季长沙一次重度霾污染天气过程气象成因分析[J].气象与环境学报，35（3）：45-51. Zhou L，Lan M C，Zhang J Q，et al.，2019.Characteristics and causes of a heavy haze event in Changsha during winter of 2017[J].J Meteorol Environ，35（3）：45-51.doi：10.3969/j.issn.1673-503X.2019.03.006.（in Chinese）.

In this article，the time series of haze days were constructed on the basis of past 64 years data which was collected from the national meteorological observatory of Wuhan and 130 climate system indexes of the National Climate Center during the time spread from 1955 to 2018.The linear trend analysis and Mann-Kendall were validated to analyze the climatic characteristics of haze days in Wuhan during this specific period.In the next stage of this research method，the stepwise multiple regression analysis was adopted in order to develop the prediction equations and quantitative monthly forecast scale of haze days.While in this piece of research，climate system index and haze days were considered as independent and dependent variables accordingly.The consequences of this research indicated that the reconstructed time series of haze days contemplated the actual long-term variations in haze days characteristics of Wuhan.During the selected span of time，the climatic characteristics in terms of annual haze days expanded significantly while the unanticipated variations were observed around the 1980s.It can also be observed that the span of seasonal haze days in winter and autumn comparatively more expanded than the summer.The number of seasonal haze days significantly increased in winter and autumn and in contrast decreased spring and summer seasons.The aggregate of monthly haze days in November，December，and January significantly increased while on other side declined in July and August.In the consideration of serious climatic conditions of haze days significantly increased and owned the gradient feature in winter and less effective in summer.The aggregate of moderate or more serious monthly haze days in November，December，and January significantly increased while on the other side declined in May，June，and August.Based on the outcome it can be concluded that discrete monthly haze days and harsh haze days owned a significant correlation with many monthly climate system indexes.The results showed that this prediction method had a great effect in Wuhu and wealth simplifying broadly in similar areas.

haze days;reconstruction;climate system index;stepwise multiple regression;monthly scale quantitative prediction

doi：10.13878/j.cnki.dqkxxb.20200725001

（责任编辑：袁东敏）