危诗敏,冯鑫媛,王式功,张峥子(.成都信息工程大学大气科学学院,高原大气与环境四川省重点实验室,四川 成都60225;2.安徽省黄山市气象局,安徽 黄山 24502)
逆温是气温随高度的升高而增加的一种现象[1].逆温层的存在抑制了污染物在垂直方向上的扩散,导致近地面污染物浓度升高,空气质量比预期变差[2-5].在欧洲城市,有研究指出贴地逆温和脱地逆温都与颗粒物浓度超标的时期相吻合[6].Rendón等[7]指出河谷内逆温层的存在、演化及其与城市热岛的相互作用强烈地影响着城市河谷外污染物浓度.姚青等[5]在一次天津重污染过程中发现,PM2.5质量浓度出现峰值时伴随有强逆温出现,逆温强度甚至高达7.69 ℃/1 00m.吴蒙等[8]指出珠三角地区污染日情况下逆温层的高度主要集中在 0~1000m,其中贴地逆温出现的频率多在 20%~35%之间.张敏等[9]研究表明,下半夜有逆温出现时,全天的 PM2.5浓度要明显高于无逆温时,尤其是在逆温消散比较晚的情况下,PM2.5日均浓度将快速增长.对地中海东岸Haifa湾的加密观测试验表明,地形陡峭地区有时会出现多层逆温现象,并且逆温强度与其内部精细结构有关[10].Wolf对卑尔根河谷地区的研究指出存在两种类型的逆温,即近地面逆温和高空逆温,它们由不同的物理机制产生[11].国内城市有时也会出现多个逆温层.乌鲁木齐 5000m以下的大气存在多个逆温层的特征[12];天津夜间温度层结复杂,也常出现多个逆温层的结构,并对大气污染物的扩散有重要影响[9];广州地区在多个逆温层出现时,容易形成灰霾天气[13];京津冀一次污染过程中有双重逆温层的现象出现,强度逐渐增加,最高达到 3.6 ℃/1 00m,是造成PM2.5爆发性增长的原因之一[14].
四川盆地 SO2、NOx和烟(粉)尘排放总量占全国总排放量的 7.3%、5.3%和 3.9%[15].高污染物排放量以及盆地地形造成的不利于大气扩散的条件,使四川盆地被列为中国第四大污染最严重的地区[16].四川盆地空气质量季节变化明显,冬季最差[17],且具有明显的二次污染的污染特征[18].首要污染物以 PM2.5、O3、PM10为主[19].其中,PM2.5污染最为严重,尤其是盆地西部和南部地区[20].同时,盆地西北部城市群受到高原东移低值系统影响,也常常发生空气污染事件[21].
逆温作为重要的污染气象条件之一,是四川盆地大气污染研究中重点关注的方面之一[22-23].四川盆地边界层到对流层中低层的逆温有 91.4%出现在5500m以下[24].Feng等[24]的研究表明,四川盆地逆温可按其底高分为3类:贴地逆温(SI-底高为0m的逆温)、脱地逆温(EI-底高大于0m且小于等于1000m的逆温)、对流层低层逆温(LTI-底高大于1000m且小于等于5500m的逆温).其中SI和EI合称边界层逆温.四川盆地因其特殊的地形和气候条件,以及紧邻青藏高原东缘的特殊位置(尤其是盆沿—盆底海拔高度差最大处超过2500m——远高于盆底区域边界层的顶部),使得这一地区的逆温特征具有独特性.表现为边界层到对流层中低层的逆温最容易出现在两个高度,即距地面600m高度以下和2200~3500m;LTI年平均频率(63.0%)和厚度(264.7m)最大,SI的强度(1.8℃/100m)最大;而对流层低层逆温频发是四川盆地逆温特征有别于国内其他三大重污染区(即京津冀、长三角、珠三角地区)的独特之处[24].值得注意的是,若将同一统计时段内(例如1月份),四川盆地SI、EI和LTI的发生频率相加,三者之和有可能明显超过百分之百,暗示了多个逆温层现象可能会阶段性呈现.关于重污染个例的零星研究也证明了多个逆温层现象的存在[25].然而,对于四川盆地这一特殊的逆温结构尚未见详细报道,缺少多层逆温的结构特征及其季节性差异的研究,尤其缺乏多个逆温层的配置关系与协同作用对大气污染综合影响的研究.因此,本文将深入研究四川盆地多层逆温的结构特征,搞清不同类型和结构的逆温对盆地大气污染的影响.以期为四川盆地环境空气质量预报和大气污染防治提供重要的科学依据.
四川盆地共有4个高空气象观测站——成都市温江站、达州市达川站、宜宾市宜宾站和重庆市沙坪坝站(图1).本文利用4站2015年1月1日~2018年12月31日每日两次(即北京时08:00和20:00)的规定等压面和特性层探空观测数据.从国家生态环境部的全国城市空气质量实时发布平台获取成都、达州、宜宾、重庆4个城市同期逐日空气质量指数(AQI)以及大气污染物质量浓度的逐时监测资料,其中包括6种大气污染物:PM2.5、PM10、SO2、NO2、O3和CO.
图1 四川盆地探空站点地理位置分布及地形(m)Fig.1 Geographical distribution of radiosonde sites and topography of Sichuan Basin
逆温层定义为温度随高度的升高而增加的气层[1].逆温出现频率定义为出现逆温的时次数占总有效观测时次数的百分比.逆温层厚度定义为逆温层顶高度与逆温层底高度之差,即D=H2-H1.其中H1逆温层底的高度,H2逆温层顶的高度.逆温强度则定义为高度每上升100m时温度的增量,即I=(t2-t1)/D×100.其中 t1逆温层底的温度(℃),t2逆温层顶的温度(℃).鉴于SI、EI和LTI层数均以 1层为主(见下文),本文将多层逆温定义为同时出现两类或三类逆温,而一般逆温则指仅出现一种类型逆温.根据多层逆温中不同类型逆温的配置情况,将其垂直结构划分为4种亚型,分别为第Ⅰ类多层逆温(SI+ EI),第Ⅱ类多层逆温(SI+LTI),第Ⅲ类多层逆温(EI+ LTI)以及第Ⅳ类多层逆温(SI+EI+LTI).
四川盆地逆温层数的频率分布如表1所示.绝大多数为 1~3层,占总数的 96.7%(08:00)和 98.2%(20:00).其中,同一时刻仅出现一个逆温层的情况最常见,占 47.4%(08:00)和 58.7%(20:00);出现 2层逆温的分别占36.6%和30.9%;同时出现4层及以上逆温层的概率很低,仅占 3.3%(08:00)和 1.7%(20:00).08:00最多为8层,出现在夏季;20是最多为6层,出现在春季.从季节分布来看,单层逆温夏季最多,冬季最少.而 2~4层逆温均为冬季最多,夏季最少.逆温层数的日变化也较为明显,20:00单层逆温明显多于08:00,而2~4层逆温则表现为08:00多于20:00,这与南京[1]的逆温层数变化特征较为一致.
表1 2015~2018年四川盆地逆温层数的频率分布(%)Table 1 Frequency distribution of the number of inversion layers in the Sichuan Basin during 2015~2018(%)
进一步统计了不同类型逆温的层数.根据定义,SI仅有1层.EI若出现几乎只有1层(表2),极少出现2层及以上(仅占1.2%~3.0%).
表2 2015~2018年四川盆地EI层数的频率分布(%)Table 2 Frequency distribution of the number of EI layers in the Sichuan Basin during 2015~2018(%)
LTI则以 1~2层占绝大多数(表 3),08:00和20:00只出现1层LTI分别占65.0%和65.6%,只出现2层的分别占28.3%和28.7%.08:00LTI最多可达6层,20:00可达4层.
表3 2015~2018年四川盆地LTI层数的频率分布(%)Table 3 Frequency distribution of the number of LTI layers in the Sichuan Basin during 2015~2018(%)
上述分析表明,四川盆地存在同时出现多个逆温层的现象.其中,SI(底高为0m的逆温)若出现则仅有1层,而EI和LTI均以1层为主.基于此,本文在前述“资料与方法”部分中根据是否同时出现两类或以上的逆温作为标准,将逆温划分为多层逆温和一般逆温,此划分依据是合理的.
2015~2018年四川盆地多层逆温和一般逆温的出现频率分别为20.1%和54.2%,表明四川盆地的多层逆温是一种并不少见的现象(图2).其中,多层逆温出现频率成都最高(36.1%),宜宾最低.盆地多层逆温 08:00出现频率(28.9%)高于 20:00(11.4%),而一般逆温则刚好相反,20:00出现频率(58.6%)高于08:00(49.9%).
图2 2015~2018年四川盆地多层逆温和一般逆温出现频率Fig.2 Frequency of multi-layer inversion and general inversion in the Sichuan Basin during 2015~2018
四川盆地多层逆温和一般逆温逐月出现频率的变化曲线均呈“U”型分布(图3),与四川盆地AQI的月变化较为一致0.其中,多层逆温冬季最频发,尤其是12月,频率达 30.6%;夏季频率最低,最小值出现在 6月,为 12.3%(图 3(a)).其中,成都冬季多层逆温频率为51.6%,甚至高于一般逆温出现频率.四川盆地一般逆温出现频率的季节变化也表现为冬季最频发,其中 1月出现频率最高,夏季频率最低(42.4%),但最小值出现在9月(图3(b)).此外,成都一般逆温出现频率夏季高、秋季低,季节变化与其他城市存在差异.
图3 2015~2018年四川盆地多层逆温和一般逆温出现频率的月变化Fig.3 Monthly frequency of multi-layer inversion and general inversion in the Sichuan Basin during 2015~2018
按逆温垂直结构将多层逆温分为4类亚型,分别统计它们的特征并与一般逆温的3类亚型进行比较.
2015~2018年四川盆地仅LTI出现频率(43.2%)明显高于其他逆温,四城市中宜宾的出现频率最高、成都最低.其次是第Ⅱ类多层逆温,出现频率为12.2%,其中成都的出现频率(31.1%)明显高于其他城市.第Ⅲ类多层逆温再次之,频率为 7.1%,达州出现频率最高(15.1%).仅SI和仅EI出现频率相近,分别为5.7%和5.4%,其中仅SI的出现频率成都最高(12.5%),而仅EI的出现频率达州最高(9.5%).第Ⅰ、Ⅳ类多层逆温出现频率最低,仅为0.3%和0.5%(图4).
图4 2015~2018年四川盆地不同结构逆温的出现频率Fig.4 Frequency of the inversions with different structures in the Sichuan Basin during 2015~2018
四川盆地的多层逆温之中,以第Ⅱ类多层逆温为主,占多层逆温总数的 60.5%;其次是第Ⅲ类多层逆温,占 35.2%;第Ⅳ类、第Ⅰ类多层逆温很少,分别占2.7%和1.6%.
不同多层逆温结构出现频率的季节变化(图 5)存在明显差异.第Ⅰ类多层逆温出现频率夏季最高、秋冬两季最低(图 5(a)).而第Ⅱ类多层逆温出现频率的季节变化则相反,表现为冬季最高、夏季最低,月变化曲线呈“U”型分布(图5(b)),其中,出现频率最大值在2月(22.3%),最小值在7月(4.7%).相比而言,第Ⅲ、IV类多层逆温的季节变化不明显,除 12月出现频率最大(12.6%和1.3%)外,其余月份变化不大(图5(c)和(d)).
图5 2015~2018年四川盆地不同结构多层逆温出现频率的月变化Fig.5 Monthly variation of frequency of multi-layer inversion with different structures in Sichuan Basin during 2015~2018
不同多层逆温的平均厚度和强度也存在差异,如图6所示.第II、Ⅲ类多层逆温的较厚(图6(a)),平均厚度分别为232.4m和236.7m,其次是第Ⅳ类多层逆温(193.6m),第Ⅰ类多层逆温最薄(174.3m).4类多层逆温均为冬季最厚,夏/春季最薄.就平均强度而言,第Ⅰ类多层逆温最强,为1.7 ℃/1 00m(图6(b)),但出现频率低(0.3%);第Ⅱ、Ⅳ类多层逆温次之,分别为1.5℃/100m 和 1.4 ℃/1 00m;第Ⅲ类多层逆温强度最弱,为 1.0 ℃/1 00m.从季节变化来看,除第I类多层逆温外,其它 3类多层逆温均为冬季最强、夏季最弱.其中,第Ⅱ、III、IV 类多层逆温冬季平均强度分别为1.8℃/100m、1.1℃/100m和1.6℃/1 00m.而第Ⅰ类多层逆温夏季最强(1.8℃/1 00m),春季最弱(1.4 ℃/1 00m),这与其逆温出现频率的季节变化较为一致.
图6 2015~2018年四川盆地不同结构多层逆温平均厚度(a)及平均强度(b)的季节变化Fig.6 Seasonal variation of(a)the average thickness and(b)the average intensity of multi-layer inversion with different structures in the Sichuan Basin during 2015~2018
逆温对四川盆地大气污染有重要影响.这里以各空气质量等级所占的百分比来衡量某种逆温条件下大气污染发生的可能性.有逆温时,四川盆地大气污染(即空气质量III级及以上)发生率远高于无逆温时(图 7).无逆温时,空气质量优良率高达 95.4%,污染率仅 4.6%.而多层逆温时大气污染发生率(39.1%)明显高于一般逆温(25.0%),表明逆温的叠加作用对大气污染的影响更大.各类逆温结构对大气污染水平的影响存在一定的差异.第Ⅰ、Ⅱ类多层逆温时大气污染发生的可能性最高(分别为 59.6%和46.8%),其次是第Ⅲ、Ⅳ类多层逆温以及仅 SI,污染发生率分别为30.1%、30.6%和33.5%,而仅EI、仅LTI的大气污染发生率最低(分别为14.5%和25.3%).此外,就中度污染及重度污染而言,多层逆温时中度和重度污染发生率均高于一般逆温,其中,第 II、IV类多层逆温的中度和重度污染发生率最高,分别为15.9%和15.6%;而仅LTI和仅EI出现时,中度和重度污染发生的可能性最小,分别占6.7%和2.3%.四川盆地极少出现严重污染,只在仅 LTI和第Ⅱ类多层逆温条件下出现过(发生率仅 0.12%和 0.04%),并且均发生在1月.
图7 2015~2018年四川盆地不同逆温条件下各空气质量等级所占比例Fig.7 Proportion of air quality grades under different temperature inversion conditions in Sichuan Basin during 2015~2018
总而言之,四川盆地的多层逆温是一种并不少见的现象,且逆温叠加作用对四川盆地大气污染的影响更大.研究表明,国内其它污染较重的地区或城市(例如京津冀、长三角、珠三角地区)也存在多个逆温层同时出现的现象.然而,这种现象大多发生在边界层内,即本文中的第Ⅰ类多层逆温,边界层逆温与对流层低层逆温的多层逆温配置很少出现,与四川盆地的多层逆温特征具有明显的差异.这些研究通常将多层逆温定义为出现两层或两层以上的逆温.北京地区多个逆温层大多出现在500m高度内0,600~800m 高度以及 1200~1400m 高度也有出现,明显抑制大气垂直方向上的扩散,使大气污染物不断地累积0.南京08:00郊区近地层的湍流发展稳定,会出现多个逆温层0.广州地区同时存在的多个逆温层主要出现在 1000m内,当贴地逆温与脱地逆温同时存在时,容易出现灰霾天气0.兰州边界层内也经常出现多个逆温层的情况,出现高度分别为 300m 和800m0.
相比而言,四川盆地多层逆温以边界层逆温叠加对流层低层逆温的配置为主要特征.这与四川盆地LTI频发这一有别于其他三大重污染区(京津冀、长三角、珠三角)逆温特征的独特之处0有关.边界层逆温由于贴近地面,对于近地面污染物扩散的抑制作用十分明显.这不仅对于四川盆地,对于其它很多城市和地区也是如此[0-0].然而,四川盆地多层逆温的独特之处在于边界层逆温之上常叠加 LTI.当 SI之上叠加 LTI时,重度污染发生率分别是仅 SI的 3.3倍(第II型多层逆温)和2.8倍(第IV型多层逆温)(图7).边界层内逆温层(尤其是SI)的存在抑制了近地面污染物向上扩散和运移,而边界层之上常出现在2200~3500m高度的 LTI则好似覆盖在四川盆地盆周山地上空的盖子,抑制其下方次级环流的向上发展,进一步加剧了对污染物扩散的抑制作用.未来很有必要对四川盆地多层逆温的形成机制及对空气污染的影响机制开展深入研究.
大气污染物浓度存在明显的季节变化[31],引入去除季节差异的污染物浓度[32],深入分析不同逆温条件对各种污染物的影响程度及其差异.四川盆地多层逆温以及仅SI对当地大气污染物的扩散有明显地抑制作用;而仅EI和仅LTI对污染物的抑制作用不明显(图8),但发生中度和重度污染的可能仍有 2.3%和6.7%(图 7).这可能与冷锋过境期间出现的锋面逆温有关.通常而言,冷锋有利于当地大气污染物的清除[33].此外,降水对大气污染物具有明显的清除作用[34].剔除前期明显降水后,其结果与上述结果基本一致.
图8 2015~2018年四川盆地不同逆温条件下6种大气污染物剔除季节差异的平均浓度Fig.8 The average concentration of six air pollutants under different temperature inversion conditions in the Sichuan Basin during 2015~2018(in the removal seasonal differences)
逆温对各种污染物的作用及影响程度还存在差异.其中,对于PM2.5和PM10,抑制作用最明显的是第Ⅳ类多层逆温,剔除季节性差异后的浓度分别为16.6μg/m3和26.8μg/m3,其次是第Ⅰ、Ⅱ类多层逆温以及仅 SI,而第Ⅲ类多层逆温的抑制作用相对较弱.对于 O3,抑制作用最明显的也是第Ⅳ类多层逆温,剔除季节性差异后的浓度为8.5μg/m3,其次是第Ⅰ、Ⅱ类多层逆温.对 SO2抑制作用最明显的是第Ⅰ类多层逆温,剔除季节性差异后的浓度为 5.8μg/m3,其余逆温条件的影响程度均较小.对于 NO2和 CO,抑制作用最明显的分别为第IV类多层逆温和第I类多层逆温.逆温对不同污染物浓度影响差异的具体原因尚不清楚,有可能与污染物本身的特性以及逆温所在的高度有关,还有待研究.
3.1 四川盆地存在同时出现多个逆温层的特点,以 1~3层为主.其中,EI以 1层为主,而 LTI以 1~.2层居多.
3.2 四川盆地多层逆温是一种并不少见的现象,2015~2018年多层逆温年平均出现频率达到20.1%.以边界层逆温叠加对流层低层逆温的配置为主要特征,这是四川盆地相比与国内其他三大重污染区(京津冀、长三角、珠三角)逆温特征的独特之处.季节变化明显,冬季多、夏季少.第Ⅱ、Ⅲ类多层逆温占绝大多数,分别占多层逆温总数的60.5%和35.2%.
3.3 逆温的叠加作用对四川盆地大气污染的影响更大.多层逆温出现时,大气污染发生的概率高于一般逆温,多层逆温尤其加大了中度及重度污染出现的可能性.当SI之上叠加LTI时,重度污染发生率分别是仅SI的3.3倍(第II型多层逆温)和2.8倍(第IV型多层逆温).
3.4 不同逆温条件对 6种大气污染物浓度的作用及其影响程度存在明显差异.第I、IV类多层逆温对污染物扩散的抑制作用最为明显.第 II类多层逆温(即 SI+LTI的配置)在冬季(污染最重的季节)出现频率高、强且厚,叠加抑制作用明显,应是逆温对四川盆地大气污染影响预报和研究的重点关注对象.