《1. 引言》

1. 引言

天气条件是重污染事件发生的重要条件，大气对污染物的稀释扩散能力由于气象条件的不同会出现较大的差异[1]。目前，对于霾天气增多的气象/气候条件的研究大多基于霾日数变化与地面气象要素（温度、风速、相对湿度、降水、海平面气压等）的统计关系[2‒13]。近些年来，影响中国天气、气候的东亚冬季风系统、ENSO（厄尔尼诺-南方涛动）事件[14‒15]以及北极海冰[16‒19]的变化对我国霾天气增多的影响逐渐得到关注。需要指出的是，气象-气溶胶反馈也发挥着重要作用。气溶胶可通过减少到达地表的太阳辐射降低近地面风速和边界层高度，从而进一步加剧近地面霾污染程度[20]。

中国东部地区冬季霾污染的形成和累积也与气候变暖密切相关。吴萍等[21]的研究结果指出，由于气候变暖，冬季气温升高，近地面相对湿度减小，大气层结稳定度的增强以及东亚冬季风强度减弱是造成中国中东部地区冬季霾增加的主要原因。Zhang等[22]指出气候变暖导致的北京风速下降与北京局地霾污染气象条件的转差显著相关。此外，全球气候变暖背景下的天气条件使北京的霾更为频繁[23]。以往研究[24‒26]表明，加速发展的工业化、现代化以及人为排放的温室气体和气溶胶等人类活动可以通过调节气象条件引发各种气象灾害。但人类活动究竟对霾天气及产生这类天气的气象条件是否产生了一定的影响？如果有，对这些影响能否进行量化，人类活动的贡献究竟有多大？关于这些问题，都是目前气候变化与污染关系研究的科学问题，也是备受政府决策部门和公众关注的重要问题。

随着气候模式水平和计算机能力的提高，利用模式进行人类活动对全球和区域尺度平均温度、极端温度和降水影响等方面的研究逐渐增多[27‒33]。此外，气溶胶污染对气候异常的重要影响也引起了广泛关注。例如，Garrett和Zhao [34‒35]研究了中纬度地区的人为污染如何导致北极冬季和春季强烈增暖。Zhao等[36]指出气溶胶污染可以扩大热带气旋的降水面积。此外，Yang等[37]发现珠江三角洲的人为污染对香港冬季地表气温的年代际变化有显著影响。Li等[25]利用观测和模式资料总结了人类活动对天气和气候的影响。利用气候模式开展与大气污染相关的研究，以往工作主要集中在评估模式对大气污染扩散条件的模拟能力[38]及对未来情景的预估研究[39]。

近两年由世界气候研究计划（WCRP）组织的20世纪气候检测归因计划（C20C + Detection and Attribution Project）[40]开展了大量考虑人类活动影响和自然强迫影响的检测与归因敏感性试验，与CMIP5† 气候模式比较计划（http://www.cesm.ucar.edu/experiments/cmip5.html）相比，这个检测与归因计划的试验设计更适合研究人类活动对全球气候的影响。利用这些试验结果，目前已开展了人类活动对极端天气和气候事件影响方面的研究[41‒48]。

然而，京津冀作为霾多发区域[26]，人类活动（包括交通运输、工业生产、土地利用以及能源消耗等）对该地区霾污染气象条件影响的检测与归因研究还比较少。因此，本文将利用综合考虑人类活动和自然强迫的影响以及单独考虑自然强迫的影响两种模拟试验结果，来分析和讨论近年来人类活动对与京津冀冬季霾气象条件的影响，C20C计划的模拟试验结果为开展此项研究提供了坚实的基础。本文其他部分安排如下：第2节介绍本研究所采用的数据与方法；第3节分析人类活动对霾气象条件的影响；第4节研究动力和热力气象因子对霾气象条件的贡献程度；第5节给出本文的主要结论及有关讨论。

《2. 数据与方法》

2. 数据与方法

《2.1. 模式试验和数据》

2.1. 模式试验和数据

本文选用了参加20世纪气候检测归因计划中的CAM5.1-1degree模式的输出结果进行分析研究。CAM5.1-1degree模式是CESM1.0.3地球系统模式的大气部分，水平分辨率是1.25°×0.9375°（经度和纬度），垂直分辨率为30层。动力学方程基于CAM5.1模式有限元内核所构建，陆面模式为CLM4.0，是CESM1.0.3中的陆面模型。

利用CAM5.1-1degree模式在C20C + D＆A计划下完成的两个基准情景下的模拟试验，其中一个试验使用观测到的边界条件[包括观测到的温室气体浓度、硫酸盐气溶胶（或有害气溶胶）、黑碳气溶胶、沙尘气溶胶、海盐气溶胶、臭氧、地表覆盖/使用、海表温度（SST）和海冰浓度（SIC）]来驱动CAM5.1模式，得到包括所有观测的外界强迫下的真实历史情景（后文简称为All-Hist）；另一个是非现实情景，试验中除了太阳辐射量和火山强迫与All-Hist相同外，温室气体浓度、气溶胶和臭氧浓度的数值使用的是1855年估计数值，这些数值在模拟过程中保持不变，SST则是根据多个CMIP5模式模拟的历史变化（仅受自然强迫的历史变化驱动）而来，海冰浓度是根据观测到的温度与冰的关系进行调整，使其与较冷的温度保持一致，此试验中的土地利用/覆盖变化与All-Hist试验保持一致，代表没有人类活动干扰下的历史情景（后文简称为Nat-Hist）。每种情景下的试验都进行了400个样本的集合，开始日期分别为1959年、1996年和2010年的1月1日，所有模拟均从相同的初始状态开始，只考虑了三维温度场的微小的均匀干扰[40]。值得注意的是，从1959年1月1日开始的一些模拟实际上是在1958年开始的，具体见表1。

两种情景下的模拟试验输出数据包括：逐日和逐月的温度、降水、气压、风速、湿度和高度场等多个要素‡。

《2.2. 方法》

2.2. 方法

根据已有研究结果，与京津冀地区冬季霾日频数气候变率显著关联的气象条件主要有动力和热力两类，其中，动力因子包括东亚冬季风和近地面风速等，热力因子包括地面气温、湿度以及逆温等。因此本文试图从这几个因子入手，对比现实情景下（All-Hist）和非现实情景下（Nat-Hist）上述这些因子的变化，来考察人类活动是否对造成京津冀地区冬季霾气象条件产生了影响。关于模式在全球范围内模拟能力的评估可参考Risser等[47]的研究工作。分析时，本文将两种情景下模拟出的气象条件在时空演变上表现出的差异认为是主要由人类活动影响所导致的。

《》

† http://www.cesm.ucar.edu/experiments/cmip5.html.

本文研究的区域为京津冀地区（36°~42.5° N, 113°~120° E），站点分布及1961‒2016年京津冀地区56个冬季（当年12月至次年2月）霾日总数分布见图1。从图1可以看出，京津冀西部和南部及环首都地区都是霾日偏多的地区，近56年冬季部分地区霾日总数超过500 d以上，平均每个冬季约有10 d的霾天气。这里霾日的判定主要利用国家气象信息中心提供的逐日天气现象观测资料，并满足以下判定条件：当天定时值（北京时02、08、14、20时）能见度出现小于10 km的记录，当天相应时次相对湿度小于90%，当天有霾天气现象记录且未出现降水、吹雪、雪暴、扬沙、沙尘暴、浮尘和烟幕等其他能导致低能见度的事件。

《图1》

图1 京津冀地区站点及1961-2016年冬季霾日总数分布（单位：d）。

此外，考虑到2013年以后气象观测仪器的更换，根据Pei等[49]的方法将2013年1月1日之后的自动观测能见度数值做了订正处理，即原自动观测值除以0.766。图2给出了1961‒2016年京津冀地区冬季霾日数的逐年演变（平均值为4.6 d）。从中选出冬季霾日数 ≥ 5 d的16个典型年（1980年、1981年、1982年、1998年、1999年、2001年、2005年、2006年、2007年、2008年、2010年、2011年、2013年、2014年、2015年和2016年）进行具体分析，以便更好地分析京津冀地区冬季霾气象条件与人类活动的关系。从图2可以明显看出，京津冀地区冬季霾日数 ≥ 5 d的年份都出现在1980年后，特别是在2013年以后，冬季霾日数呈现出突然增加的趋势，冬季霾日数可达20 d以上。造成这种现象的原因除与我国经济快速发展导致的排放增加有关外，与天气气候条件的变化也密切相关，主要表现为东亚冬季风减弱、近地面风速减弱、地面气温升高、大气相对湿度降低以及大气稳定度和逆温增加等[21,50‒52]。

《图2》

图2 1961-2016年京津冀地区冬季霾日发生次数时序图（单位：d），虚线表示霾日数5日平均线。

《3. 人类活动对京津冀地区冬季霾气象条件的影响》

3. 人类活动对京津冀地区冬季霾气象条件的影响

《3.1. 动力条件》

3.1. 动力条件

东亚冬季风是北半球冬季最为活跃的环流系统之一，对我国的天气和气候有着重要的影响。东亚冬季风主要成员包括低层风场上亚洲东岸引导极地冷空气南下的西北气流、海平面气压场上的西伯利亚高压和500 hPa位势高度上中国大陆东岸至日本上空的东亚大槽等。强东亚冬季风常带来较强的冷空气和偏北大风，有利于污染物的稀释和扩散。反之，弱东亚冬季风与冬季霾污染增多密切联系[25]。因此，地面风速与霾天气的发生密切联系。这里选取500 hPa高度场在（25°~40° N, 110°~130° E）区域的平均值来代表东亚季风的强度，文中将该区域平均的500 hPa高度场经过标准化处理后将其结果乘以-1，以此来定义东亚冬季风指数。东亚冬季风指数越大表示冬季风越强，反之亦然[21,53]。

《3.1.1. 东亚冬季风》

3.1.1. 东亚冬季风

东亚大槽是北半球冬季对流层中层重要的环流系统。图3给出了考虑了人类活动强迫的现实情景（All-Hist）与非现实情景（Nat-Hist）下京津冀地区冬季500 hPa高度场变化时序图。考虑人类活动强迫的现实情景（All-His）下，京津冀地区冬季高度场明显偏高，特别是在20世纪90年代以后这种差别更为显著，说明人类活动对对流层中层高度场的影响十分明显。此外，两种情景下京津冀地区冬季高度场有着较为一致的阶段性变化，即20世纪60~70年代为高度场降低期，其后则经历了近20年的高度场波动性增高期。从21世纪开始，高度场的变化进入一个较为平稳的时期，变化幅度较小。对比京津冀冬季霾日偏多年冬季的高度场变化，可以清楚看到，人类活动影响下，霾日偏多年京津冀基本处在高度场正距平控制下。高度场的变化，一方面会伴随有风场的调整，另一方面通过气流的上升冷却和下沉增温作用影响地面气温的变化。关于风场和气温的变化将在以下章节详细分析。

《图3》

图3 现实情景（All-Hist，红色线）与非现实情景（Nat-Hist，蓝色线）下京津冀地区冬季500 hPa高度场变化时序图（滑动平均）（相对于Nat-Hist情景下1981‒2010年平均值的距平）。阴影区表示100个样本的25%~75%的范围，实线表示100个样本的中值，灰色点线为霾日偏多年份。1 gpm = 9.8 J·kg-1

进一步分析东亚冬季风指数变化时序图（图4）可以发现，20世纪80年代之前，两种情景下，东亚冬季风的强度都有增强的趋势且Nat-Hist情景下东亚冬季风强度总体要强于All-Hist情景下的结果。此后，All-Hist情景下东亚冬季风强度减弱明显；Nat-Hist情景下东亚冬季风强度则表现出逐渐加强的趋势。从20世纪80年代中期开始，All-Hist情景下东亚冬季风年代际减弱趋势非常明显，这与观测结果基本一致[20,54]。此外，霾日偏多的年份基本上都处在东亚冬季风减弱期。这表明在人类活动影响下，东亚冬季风的减弱更为明显，冷空气不活跃且势力偏弱，大气停滞-静稳日偏多，大气扩散颗粒物的能力偏弱。与此同时，这一时期中国经济进入快速发展阶段，城市人口、工业生产及机动车辆排放到大气中的气溶胶粒子大量增加，因而也进一步加剧了京津冀地区霾天气的发生与发展。

《图4》

图4 现实情景（All-Hist，红色线）与非现实情景（Nat-Hist，蓝色线）下东亚冬季风指数变化时序图（滑动平均）（相对于Nat-Hist情景下1981‒2010年平均值的距平）。阴影区表示100个样本的25%~75%的范围，实线表示100个样本的中值，灰色点线为霾日偏多年份。

《3.1.2. 近地面风速的变化》

3.1.2. 近地面风速的变化

东亚冬季风的强弱直接影响近地面风速的大小，近地面风力条件是污染物颗粒稀释扩散的主要动力条件。图5的结果表明，20世纪80年代中期以来影响京津冀地区的近地面纬向风（925 hPa）在人类活动强迫情景（All-Hist）和非现实情景（Nat-Hist）下均表现出明显减小趋势，且其演变趋势较为一致。但考虑人类活动强迫情景（All-Hist）下，近地面纬向风减小趋势更加明显[图5（a）]。对近地面经向风的分析表明，南风在人类活动强迫情景（All-Hist）下的模拟结果明显较非现实情景（Nat-Hist）下偏强[图5（b）]。这说明在人类活动强迫下，南风明显增强，从而更有利于京津冀南部及以南地区的污染物向北输送，不利于近地面附近的污染物向区域外扩散。上述分析表明在人类活动影响下，与东亚冬季风减弱相伴随的低层西风减弱、南风加强是造成霾天气增加的原因之一。

《图5》

图5 京津冀地区All-Hist（红色线）与Nat-Hist（蓝色线）情景下冬季地面纬向风（U）（a）和经向风（V）（b）时序变化图（滑动平均）（相对于Nat-Hist情景下1981‒2010年平均值的距平）。阴影区表示100个样本的25%~75%的范围，实线表示100个样本的中值，灰色点线为霾日偏多年份。

《3.2. 热力条件》

3.2. 热力条件

3.1节分析表明，人类活动确实对与霾相关的动力气象条件产生了一定的影响。本小节将主要分析人类活动强迫对与霾相关的热力气象因子的影响。

《3.2.1. 气温的变化》

3.2.1. 气温的变化

图6为京津冀地区在现实情景（All-Hist，红色线）与非现实情景（Nat-Hist，蓝色线）下冬季地面气温变化时间序列图。可以看出，与东亚冬季风的年代际变化一致，在All-Hist情景下，京津冀地区冬季气温在20世纪80年代中期也发生了一次年代际变化，80年代中期之前东亚冬季风偏强，京津冀地区气温相对偏低，20世纪90年代之后，京津冀地区冬季气温逐渐进入一个相对偏高阶段，之后东亚冬季风偏弱10年左右。从图6中还可以看出，在20世纪90年代初期，霾天气发生较少的年份，京津冀地区冬季气温明显偏低，随后温度突然上升，21世纪初期一直保持温度较高的状态，相应霾日偏多的年份也逐渐增加。对比Nat-Hist情景下的京津冀地区冬季温度的变化可看出，从20世纪70年代开始，京津冀地区冬季的温度基本没有发生明显的变化，基本维持在平均水平附近。

《图6》

图6 现实情景（All-Hist，红色线）与非现实情景（Nat-Hist，蓝色线）下京津冀地区冬季地面气温时序变化图（相对于Nat-Hist情景下1981‒2010年平均值的距平）。阴影区表示100个样本的25%~75%的范围，实线表示100个样本的中值，灰色点线为霾日偏多年份。

《3.2.2. 相对湿度的变化》

3.2.2. 相对湿度的变化

相对湿度是霾形成的重要条件之一，与其他季节相比，在冬季，京津冀地区降水日数明显偏少[13]，相对湿度的作用主要表现为对雾霾形成过程的增湿作用[55]。在现实情景（All-Hist）和非现实情景（Nat-Hist）下，京津冀地区的相对湿度有着较为相似的演变特征（图7）。从20世纪80年代中后期至20世纪90年代中期，两种情景下相对湿度均表现出明显的年代际增湿特征。此后，相对湿度突然减小，进入相对偏低的阶段（也即雾少霾多的阶段）[51]。同时也看出，在All-Hist情景下，相对湿度增加明显大于Nat-Hist情景，表明人类活动可能使得京津冀地区的相对湿度增加更明显，这种相对湿度的增加在一定程度上会使得北方干燥大气中颗粒物的吸水性增加，导致霾天气增多。

《图7》

图7 现实情景（All-Hist，红色线）与非现实情景（Nat-Hist，蓝色线）下京津冀地区冬季地面相对湿度时序变化图（相对于Nat-Hist情景1981‒2010年平均值的距平）。阴影区表示100个样本的25%~75%的范围，实线表示100个样本的中值，灰色点线为霾日偏多年份。

《3.2.3. 逆温的变化》

3.2.3. 逆温的变化

霾天气的强弱除了与对流层中低层大气层结不稳定性有关外，也与近地面的逆温状况有关。为了分析大气近地面的逆温状况，图8给出了现实情景（All-Hist，红色线）与非现实情景（Nat-Hist，蓝色线）下京津冀冬季温度垂直变化剖面图。在All-Hist情景下，对流层中层气温表现出较为一致的逐渐增强趋势，特别是从20世纪90年代中后期开始，对流层气温的增加有逐渐向高层发展的趋势。21世纪以来，中层气温增加更加明显，增温幅度达1.0 ℃以上。而在非现实情景下（Nat-His）情景下，除20世纪60~70年代外，对流层中层大气并没有表现出持续的、明显的增暖现象，这说明在全球变暖背景下，低层大气的热力状况不利于污染物的扩散，逆温特征不明显。

《图8》

图8 现实情景（All-Hist）（a）与非现实情景（Nat-Hist）（b）下京津冀冬季温度垂直变化剖面图（相对于Nat-Hist情景下1981‒2010年平均值的距平）。引自Wang等[9]。

《4. 动力和热力气象条件对霾日数的贡献》

4. 动力和热力气象条件对霾日数的贡献

上述分析表明，大气的动力和热力作用都会对雾霾天气的发生有重要影响。虽然与人类活动有关的气溶胶与霾天气有密切联系[50]，但人类活动影响下的气象条件对雾霾天气的贡献有多大是一个需要进一步研究的问题。这一节分别对现实情景（All-Hist）和非现实情景（Nat-Hist）下的动力和热力气象条件对霾产生的贡献进行分析。

如前所述，动力条件包括：500 hPa高度场（H500）、东亚冬季风指数（EAWMI）、近地面纬向风（U）和经向风（V）；热力条件包括：温度（Tas）、相对湿度和逆温。针对上述因子，文中利用多元线性回归方法建立京津冀地区1961‒2016年京津冀冬季霾日数的线性回归方程。为了检验各个气象因子的独立性，文中计算了这些因子之间的相关系数。计算后发现，H500与EAWMI之间存在显著负相关，相关系数为-0.91。此外，H500与Tas也呈显著正相关，相关系数为0.74。因此，在建立多元线性回归方程时不考虑动力因子H500。其他因素之间也存在一定的相关性，但相关系数较小。

现实情景下的线性回归方程为：

Haze days = -2.735EAWMI - 2.527U - 12.179V + 2.244Tas-3.266

非现实情景下的线性回归方程为：

Haze days = 2.1777EAWMI + 4.1085U - 0.5391V - 4.2665Tas + 0.6966

其中，霾日数单位为d；EAWMI单位为位势米（gpm, 1 gpm = 9.8 J·kg-1）；纬向风（U）和经向风（V）单位为m·s-1；温度单位为℃。

图9给出了由上述两个回归方程计算得到的现实情景（All-Hist）和非现实情景（Nat-Hist）下动力和热力气象因子共同以及两个因子单独的京津冀冬季霾日多元线性回归的逐年变化曲线以及观测的京津冀冬季霾日数变化曲线。可看出两种情景（表2）下，利用所有因子回归的京津冀冬季霾日数与观测值之间具有较好的一致性，二者的相关系数分别为0.48和0.37（均通过了置信水平为99%的检验），基本上能够反映观测到的京津冀冬季霾日数随着时间增加的变化趋势，回归值对观测值方差的解释可达22.5%和14.2%，这一结果说明人类活动对京津冀冬季霾日的变化有明显的正贡献。当只考虑动力气象因子时，二者的相关系数分别为0.45和0.35，回归值对观测值方差的解释可达20.1%和12.4%。需要注意的是，All-Hist情景下动力因子的方差贡献明显较Nat-Hist情景下偏大。当只考虑热力因子时，二者的相关系数均为0.29，两种情景下的方差贡献基本相当。

《图9》

图9 1961‒2016年霾日数变化曲线（虚线；单位：d）以及利用气象动力和热力因子构建的霾日数变化多元线性回归曲线（实线；单位：d）。（a）所有因子；（b）动力因子；（c）热力因子。

《表2》

表2 所有因子及动力和热力因子统计分析结果

Statistical analysisAll factorsDynamical factorsThermal factors
Nat-HistAll-HistNat-HistAll-HistNat-HistAll-Hist
F test**—***—
Analytic variance14.2%22.5%12.4%20.1%8.5%8.4%
TCC0.370.480.350.450.290.29

从上述的统计分析结果可看出，人类活动（All-Hist情景下）对京津冀冬季霾日增多的贡献大于只考虑自然强迫，其中动力气象因子对京津冀冬季霾日产生的贡献更大；热力气象因子对霾日的贡献在两种情景下相差不大。

《5. 结论与讨论》

5. 结论与讨论

本文利用利用20世纪气候检测归因计划的CAM5.1-1degree模式的模拟结果，即综合考虑人类活动和自然强迫（All-Hist）以及单独考虑自然强迫（Nat-Hist）的两种气候数值模式的试验结果，对近年来人类活动对与京津冀地区冬季霾事件相关的气象条件的影响进行了分析和讨论，主要结论为：

（1）考虑人类活动影响后，与京津冀冬季霾事件相关的动力条件明显不利于霾的扩散，其中主要包括东亚冬季风明显减弱、对流层低层异常南风增强、近地面风速减弱。

（2）考虑人类活动影响后，与京津冀冬季霾事件相关的热力条件包括地面气温明显增加、水汽输送增强和湿度增加、逆温增强等均明显利于霾的形成和累积。

（3）利用多元回归对影响京津冀霾日的动力和热力气象因子的相对贡献的分析表明，人类活动对京津冀霾日增多的贡献大于自然强迫，其中人类活动影响下动力气象因子对京津冀霾日产生的贡献大于自然强迫；热力气象因子对霾日的贡献在两种情景下相差不大。

上述分析主要是定性诊断人类活动对与霾相关气象条件的影响，未来还将进一步利用最优指纹法等多种统计工具对人类活动对霾气象条件的贡献进行进一步的定量检测和归因研究。同时，还可对其他影响霾发生的气象因子进行分析；本研究只选取了一个模式的多个样本，后续研究可进一步利用多模式多样本的结果进行分析。

除了气象条件之外，气溶胶排放源[56]、臭氧[57]以及植树造林[58]等其他因素也与我国的霾变率有关。因此，进一步探索这些因子的贡献将有助于更好地了解我国霾天气的形成。