京津冀地区重污染天气过程的污染扩散气象条件探析

　　1前言
　　近年来，我国东部地区出现大气污染的天数迅速增加，主要污染物是人类排放的气溶胶和气态污染物通过光化学过程转化形成的二次气溶胶。污染形成的主要原因是污染物的大量排放，但受客观气象条件的影响，当出现不利于污染物输送扩散的天气时，易造成污染物的集聚、堆积。出现严重污染时空气质量极差，大气能见度降低，给社会和人类带来严重影响。京津冀地区人口密集，以汽车工业、电子工业、机械工业、冶金工业为主，是全国主要的高新技术和重工业基地，造成了该地区大气污染物排放量大，一旦出现不利于污染物扩散的气象条件，极易发生严重的大气污染现象。
　　稳定的大气环流形势为大气污染的发生和发展提供了有利的大气环境场。有研究表明，很多大气污染的发生与反气旋的停滞密切相关。反气旋在某地区长期滞留期间，地面风速可以下降到很小，近地面环流减弱，因此不足以吹散积累起来的污染物，造成空气质量下降。WangandAngell（1999）定义了空气滞留日，从天气学角度描述这种有利于大气污染发生的气象环流背景条件。当某地满足海平面地转风小于8ms，对流层中层（500hPa）水平风速小于13ms，并且没有降水，则定义为一个空气滞留日。其中，如果对流层中低层（850hPa至地面）出现逆温，将海平面地转风适当放宽10（小于8。8ms）。用海平面气压梯度表示近地面的平均流场特征比用近地面风速更具有代表性，因为后者受到地表等复杂因素的影响，其本质是由模式的边界层参数化分析得来的。从天气学意义来讲，当某地区发生滞留现象时，天气系统长期维持或是缓慢地移动，地面风速和500hPa风速小，说明对流层中下层对流较弱，使得大气水平扩散能力小，污染物不易被输送到其他地区而造成长期积累发生污染。另一方面，若850hPa以下有异常逆温存在，更加稳定的低层大气使得空气的垂直混合减弱，污染物的垂直输送得到抑制。加上没有降水的湿沉降作用，这些异常大气环流背景条件都有利于污染天气的维持和发展。吴兑等使用区域近地层风矢量和算法定义气流停滞区，分析了2013年1月典型雾霾个例气流停滞区的形成过程，认为气流停滞区的形成反映了区域平流输送条件，气流停滞区造成污染物的停滞、积累，结果证明霾过程中，华北平原均出现明显的气流停滞区，区域矢量和很小，不利于空气中污染物的水平扩散。
　　2研究方法
　　2。1污染气象条件分析
　　根据WangandAngell（1999）的研究，海平面地转风和500hPa风速是表征空气滞留的两个重要天气条件。当海平面地转风和500hPa风速都分别小于一定值时，说明天气系统在某地区持续或缓慢地移动，对流层中低层水平风速小，大气水平输送弱，就不利于污染物的扩散，因此本文将通过分析污染过程中海平面地转风和500hPa风速的特征了解污染过程中的天气环流背景。
　　2。2中尺度数值模拟
　　WRF（WeatherResearchandForecastModel）是由美国国家大气研究中心（NCAR）和国家环境预测中心（NCEP）联合开发的新一代中尺度数值天气预报模式，在目前的天气预报业务和大气科学研究中都得到了广泛应用。随着WRF中边界层参数化的不断改进，WRF对大气边界层模拟能越来越准确地描述真实大气的边界层特征。本文WRF模拟使用NCEP11分辨率的FNL资料提供初始场和边界条件，设置3层嵌套的模拟区域，网格分辨率分别为45、15、5公里。如图5。1所示，第一层区域覆盖东亚大陆和周围海域；第二层区域覆盖了整个中国；最内层区域包括了京津冀、山西省以及山东北部地区。垂直方向分36层，模式顶高在50hPa，采用上疏下密分层以更精确地模拟大气边界层的结构，其中前6层所在高度大约分别为：8，29，53，82，115和152米。第一、二层嵌套选用KainFritsch积云对流参数化方案（Kain，2004），最里层5公里水平分辨率能够解析积云对流方程，因此不需要使用参数化。第一层选用Thompson微物理过程（Thompson，Field，RasmussenandHall，1989），第二、三重采用WSM6微物理过程（HongandLin，2006）。近地面层方案需要和边界层方案配合使用，第一层区域近地面和边界层方案分别采用Eta方案和MYJ方案（Janjic，2002），其他两层都分别采用MM5方案（Beljaars，1995）和ACM2方案（Pleim，2007）。其他物理过程参数化方案设置如下：长波辐射采用RRTM方案（Mlaweretal。，1997）；短波辐射采用Dudhia方案（Dudhia，1989），每30分钟调用一次辐射方案；陆面过程使用Noah方案（ChenandDudhia，2001）。通过在45km的粗网格中运用张弛逼近（Nudging）技术，对2013年12月19日和2014年2月1926日两个重污染过程分别进行了连续的数值模拟。WRF分别从201311300800LST和201402180800LST开始积分，其中前16小时作为模式启动的调整时间，第一层网格6h输出一次结果，第二、三层每1h输出一次结果。
　　本文所做的研究中用到的资料有：（1）NCEPNCAR日平均再分析资料（Kalnay等，1996），时间范围是20032012年每年的2月和12月，以及2013年的12月和2014年的2月，水平分辨率为2。52。5，垂直分辨率为17层，本文分析中用到了500hPa的风分量资料和海平面气压资料；（2）NCEPFNL分析资料（finalanalysesdata），时间范围为北京时间2013年11月30日08时到2013年12月10日08时和2014年2月18日08时到2014年2月27日08时，水平分辨率为11，时间分辨率为每日4次（02、08、14、20时）。（3）中华人民共和国环境保护部发布的2013年12月和2014年2月各城市空气质量指数（AQI）数据。
　　3结果分析
　　3。1通风量
　　WRF数值模式是高分辨率的动力降尺度模式。通过WRF的第三层嵌套5km分辨率模拟的气象要素场结果，运用（4）式计算这次污染过程平均大气边界层通风量分布（图4。1。2）。由图可见，太行山以东的京津冀地区，通风量为明显低值分布，石家庄、保定、衡水等地的通风量均不足1000m2s。通风量体现了边界层高度和边界层内风速对污染扩散的综合作用，混合层高度的日变化特征显示，混合层高度在白天14时左右达到日最大值。根据数值模拟结果，对2013年12月19日重污染天气过程中石家庄、邢台、北京和天津等城市14时平均通风量气象条件的分析结果显示，京津冀地区各城市日最大混合层高度平均值在448m649m之间，区域平均值约539m；各城市日最大通风量平均值在1211m2s3360m2s之间，区域平均值约2543m2s。
　　3。22014年2月1926日重污染天气过程分析
　　3。2。1空气滞留气象条件
　　京津冀地区2014年2月1925日严重污染过程中500hPa风速和海平面地转风距平分布（图4。2。1）也表现出明显低于10年气候平均值的异常特征。500hPa风速的低值区分布在京津冀地区以及辽东半岛，比10年气候平均值低12ms以上，与此次污染过程中石家庄、邢台、保定、北京等城市的重污染较为一致。通过对这次污染过程中京津冀地区各城市空气滞留气象条件统计分析得知，2014年2月1926日京津冀地区严重污染天气发生期间，500hpa平均风速较近10年同期下降了约4959，平均下降55。1，平均海平面地转风下降约642，平均下降22。9。六级污染天数达到5天及以上的石家庄、邢台和北京的500hPa平均风速和海平面地转风分别降低51。356。6和8。736。2。
　　3。3地形对河北南部城市大气污染的作用
　　石家庄、邢台、衡水等京津冀南部城市位于太行山脉的东侧，为了分析山脉阻挡对污染的作用，利用模式垂直分层加密对流层低层的优点，设置模式第三层嵌套的垂直分层为36层，其中第5、8、11和17层所在高度大约为100m、250m、500m和1000m，得到2013年12月19日污染过程中每日08时（北京时间）的各高度层水平风速分布。图4。3为2013年12月2日、5和8日的各层水平风速图，其他图略。通过水平风速和空气质量指数对比发现，二者之间有着直接的关系。2日京津冀南部各城市均出现五级污染，08时的水平风速图在太行山以东为明显的风速低值区，平均风速约3ms，称之为风影区。太行山脉大部分海拔在1200m以上，因此在约1000m高度上仍然可看到有风影区的存在。持续到5日08时，之前的风影区变为明显相反的风速大值区，累积的污染物得到部分扩散、稀释，5日京津冀的污染明显得到缓解。6日到8日，风影区又逐渐显现出，从8日08时的风速图可看出明显的风影区，污染程度也随之加剧，8日达到此次污染过程的AQI峰值，以后随着气象条件逐渐变好，堆积的污染物不断被输送和沉降，空气质量变好，能见度提高。由此可见，太行山对西风气流的阻挡是造成河北南部城市大气污染的重要原因。
　　4结论
　　本文通过再分析全球大气环流模式格点资料和WRF中尺度数值模拟，从天气学和大气边界层两个角度分析了2013年12月和2014年2月我国两次严重大气污染过程中空气停滞气象条件和通风量的特征及其作用，得到研究结论如下：
　　a。两次严重大气污染天气发生期间，500hPa等压面上的平均风速和平均海平面地转风表现为明显的气候异常特征，500hpa平均风速较近10年同期分别下降了约30。6和55。1，平均海平面地转风均下降约23。
　　b。京津冀地区在两次严重污染天气过程中，区域平均日最大混合层高度分别为539m和618m，区域平均日最大通风量分别为2543m2s和2212m2s。
　　c。太行山对西风气流的阻挡是造成河北南部城市大气污染的重要原因。
　　本文仅仅分析了城市大气污染发生时的污染扩散天气条件，没有考虑降水带来的污染物湿清楚过程，以及温度、湿度、辐射等对二次光化学污染形成的影响，这些重要的污染气象条件还有待进一步研究。