1. 引言
逆温通常而言是指发生对流层内气温随高度变化而增加的现象 [1],业内对大气逆温层的研究主要集中于近地空间或大气边界层以内,逆温层阻碍了近地空间污染物的向上扩散,会极大地影响当地空气质量。胡义成等利用L波段雷达探空资料对乌鲁木齐近地逆温层特征进行研究发现逆温现象在乌鲁木齐季节变化特征明显 [1],于文立等对临江市2011~2017年低空逆温层特征进行了分析,研究表明临江市逆温主要以贴地逆温和低悬逆温为主 [2],张雪玲等对南阳地区边界层逆温特征研究发现逆温与PM10浓度存在必然联系,是影响南阳地区空气质量的重要因子 [3],业内有不少学者对各地近地空间的逆温层特征进行了大量研究 [4] - [10],重点强调了逆温层与污染物扩散之间的关系 [11] [12],张峻和王海军通过南方地区多个探空站点数据研究了1960~1999年冰冻天气过程中的低空逆温层特征,研究表明西南地区逆温强度明显强于其余地区 [13],孙建华和赵思雄对2008年南方地区持续雨雪冰冻灾害天气的层结结构研究发现,准静止锋和逆温层是此次冰冻天气的重要影响因素 [14]。
云贵高原独特的地理位置使得其在冬半年长时间受昆明准静止锋控制,贵阳市地处云贵高原东麓,域内海拔普遍在1000米以上,是来自偏北方向和偏东方向冷空气入侵云贵高原的前沿阵地,而来自印度洋的西南季风则经云贵高原西部到达贵州中部地区,冷暖空气层的特殊结构使得逆温层得以形成并维持,本文通过对近40年贵阳市探空数据进行气候分析和逆温层特征研究,总结了贵阳市近地空间逆温层厚度、强度以及逆温层底高和顶高的分布特征,并揭示了其周期性的变化规律,研究将有助于理解准静止锋的温廓结构,探究逆温层在贵阳地区不同灾害天气中的影响和作用。
2. 资料和方法
本文选取的资料包括:1) 贵阳市东山站(编号57816) 1980~2019年高空探测资料;2) 地面常规观测资料;3) 欧洲中心ERA5 0.25˚ × 0.25˚逐小时高空再分资料;4) 全球数字高程模型DEM地形数据。
探空观测每日进行两次,分别为北京时08时及20时,本文统计1980~2019年贵阳市东山探空站共计得28561时次探空数据,缺测率为2.2%,用Mann-kendall非参数检验方法 [15] [16] 对逆温时次随时间变化进行突变检验分析,用Morlet小波 [17] 分析了逆温发生的周期变化特征,本文在研究分析时已对时间序列进行统一标准化处理。
本文关注的逆温层为地面至地面以上5 km高度区间(即本文所指近地空间,下同)发生的任意逆温情况,逆温底高是指逆温层底的海拔高度,逆温顶高是指逆温层顶的海拔高度,逆温厚度是指逆温层的厚度(单位:m),逆温类型分为两类:贴地逆温是指逆温起始高度为地面海拔高度的逆温层;悬浮逆温是指逆温起始高度大于地面海拔高度逆温层,逆温强度:
(单位:℃/100m),其中
(T2为逆温层顶气温,T1为逆温层底气温);
(H2为逆温层顶海拔高度,H1为逆温层底海拔高度),逆温强度表征为在逆温层中每升高100米增加的温度。
本文所涉及的时间均为北京时,涉及的高度均为海拔高度。
3. 逆温时次变化特征
3.1. 年际变化
统计1980年至2019年贵阳市东山站探空数据发现,近40年贵阳地区近地空间逆温出现的总时次为18,790次,出现频率为66%,其中08时占比53%,20时占比47%,如图1所示,贵阳地区平均每年发生逆温时次474.25次,其中08时逐年平均251.425次,20时222.825次,08时在所有年份的逆温时次均要高于20时,从年代变化来看,1997年在所有年份中逆温时次出现最少,总时次仅324次,其中08时197次,20时127次,均大幅低于整体平均值,1992年在所有年份中逆温时次出现则最多,总时次为545次,其中08时284次,20时261次,由此可见从1992年到1997年总逆温时次有陡然下降的趋势。

Figure 1. Time distribution of low-level temperature inversion in Guiyang from 1980 to 2019
图1. 1980年至2019年贵阳近地空间逆温时次分布
Mann-kendall法是一种非参数的检验方法,样本不必遵从某一特定的分布,同时也不受个别异常值的干扰,能够客观地表征样本序列的整体变化趋势 [17]。如图2所示为贵阳地区近地空间逆温时次年变化序列的M-K突变检验,由图可见,UFk曲线和UBk曲线在近40年中总共有三次相交,意味着逆温时次有三个明显的突变时刻,分别是1980年代初、1990年代初以及2000年代末,具体而言,从1980~1990年代逆温时次呈现出上升趋势,而从1990年至2000年代则呈显著的下降趋势,2000年代以后则又呈缓慢上升趋势。

Figure 2. Mann-kendall mutation test of low-level temperature inversion in Guiyang from 1980 to 2019
图2. 1980年至2019年贵阳近地空间逆温时次Mann-kendall突变检验
根据逆温层的起始位置将逆温分为贴地逆温和悬浮逆温,在统计时若贴地逆温与悬浮逆温同时出现,则单独计数,但若同一时次有多层悬浮逆温,则悬浮逆温时次合计为一次,统计结果表明,近40年贵阳地区贴地逆温总发生时次为2961次,悬浮逆温总发生时次为17,596次,其中贴地逆温与悬浮逆温一起出现的情况共有1767次,由此可见对于贵阳地区而言,悬浮逆温是主要的逆温现象,占比超过了85%。如图3所示,
(a) 贴地逆温时次年际变化
(b) 贴地逆温时次M-K检验
(c) 悬浮逆温时次年际变化
(d) 悬浮逆温时次M-K检验
Figure 3. Annual variation and Mann-kendall mutation test in different temperature inversion types of low-level temperature inversion in Guiyang from 1980 to 2019
图3. 1980年至2019年贵阳近地空间不同逆温类型逆温时次年际变化及Mann-kendall突变检验
贴地逆温时次在近40年里整体呈下降趋势,08时占比70%,20时占比30%,二者相差十分显著,M-K检验显示在2000年至2010年期间有明显突变,突变以后下降趋势更加显著;悬浮逆温的情况与贴地逆温有所不同,其年际变化趋势与逆温总时次变化相一致,M-K检验则显示自2000年代中后期以后上升趋势显著增加,另外对于悬浮逆温而言,08时和20时逆温出现的频率相差不大,08时占比51%,20时占比49%。
3.2. 季节变化
表1显示了1980年至2019年贵阳地区不同季节逆温发生的总时次分布,冬春季发生频率高于夏秋季,图4是各季节逆温时次的年际变化及M-K突变检验,一个比较明显的特征是,冬季逆温时次在2000年代中期出现明显的突变,突变之后呈显著上升趋势,秋季逆温时次在2000年代中期同样有突变,突变

Table 1. Distribution of temperature inversion in different seasons from 1980 to 2019
表1. 1980年至2019年不同季节逆温时次分布
后呈缓慢上升趋势,春夏季则表现出了和秋冬季相反的特征,自1990年中后期以后,春夏季的逆温时次均呈现出波动的下降态势。
前文提到贴地逆温在近40年里整体呈下降趋势,悬浮逆温则经历了1990年代的下降后转为上升趋势,图5为贴地逆温及悬浮逆温不同季节的年际变化特征,可见对于贴地逆温而言,四季均呈现出相同的下降趋势,而对于悬浮逆温来说,如前文所叙与总体逆温时次一样,秋冬季呈上升态势,春夏季则呈下降态势。
(a) 贴地逆温时次各季节年际变化
(b) 悬浮逆温时次各季节年际变化
Figure 5. Annual variation of different types and different seasons of temperature inversion in Guiyang from 1980 to 2019
图5. 1980年至2019年贵阳近地空间不同类型逆温各季节时次年际变化
3.3. 月变化
图6显示了近40年贵阳地区近地空间不同类型逆温时次的月分布情况,为有效对比,将逆温时次数值进行了标准化处理,由图可见,悬浮逆温与总逆温时次呈现出高度一致的状态,整体呈U型分布,11月至2月是逆温发生的高频月份,6~8月为逆温发生的波谷月份,而贴地逆温则表现出了不同的分布特征,从1月至7月波动幅度较小,从8月至12月则表现出了显著的上升趋势,并同样在12月达到峰值。

Figure 6. Monthly distribution of diffrent types of low-level temperature inversion in Guiyang from 1980 to 2019 (normalized value)
图6. 1980年至2019年贵阳近地空间不同类型逆温时次月分布(标准化值)
3.4. 逆温时次的周期变化
小波分析是应用广泛的时频分析工具 [17] [18],其数学表达式见式(1):
(1)
式1中:
为小波系数;a为伸缩因子;b为平移因子;t是时间;
为任意平方可积函数,即研究对象的时间序列;
为一个基本小波函数,
为共轭函数,本文采用Morlet为母小波,对近40年贵阳地区逆温时次数据进行周期性分析。
图7和图8分别是近40年贵阳地区贴地逆温和悬浮逆温发生时次的小波分析综合图,两者表现出了不同的周期性特征。对贴地逆温层而言,能量图在周期2~4年,4~8年均出现了大值区,预示着贴地逆温层存在着2~4年的年际变化周期和4~8年的年际变化周期,特别是在2010年之前周期性特征较为明显,但2010年以后周期性特征不再显著;对于悬浮逆温层而言则存在着2年左右以及4~8年的年际变化周期和8~16年的年代际变化周期,其中2年左右及4~8年的年际变化周期在2010年之前较为显著,2010年之后不再显著,而8~16年的年代际变化周期则仍然较为显著。
4. 逆温指数统计
4.1. 逆温厚度
统计近40年贵阳地区近地空间逆温层厚度的分布发现,贴地逆温层平均厚度为345.7米,中位数为199米,悬浮逆温层平均厚度为471米,中位数为334米,从年代变化来看(图9),贴地逆温层的厚度在1990年至2000年代呈现出上升趋势,从2000年代以后则呈显著的下降趋势,对于悬浮逆温层而言,2000年代以前逆温层厚度波动振幅不大,但从2000年以后则呈显著下降趋势,M-K检验显示其在2000年代中后期存在着明显的突变特征。结合前文可知,贴地逆温层从2000年以后不仅发生时次开始减少,逆温厚度也显著降低,而悬浮逆温层从2000年代以后发生时次呈缓慢波动上升趋势,但逆温厚度却下降趋势明显。

Figure 7. Wavelet analysis of low-level temperature inversion from surface in Guiyang from 1980 to 2019
图7. 1980年至2019年贵阳近地空间贴地逆温时次小波分析

Figure 8. Wavelet analysis of low-level temperature inversion above surface in Guiyang from 1980 to 2019
图8. 1980年至2019年贵阳近地空间悬浮逆温时次小波分析
(a) 不同类型逆温厚度年际变化
(b) 贴地逆温厚度M-K检验
(c) 悬浮逆温厚度M-K检验
Figure 9. Annual variation of low-level temperature inversion thickness and Mann-kendall mutation test in Guiyang from 1980 to 2019
图9. 1980年至2019年贵阳近地空间不类型逆温厚度年际变化及Mann-kendall突变检验
4.2. 逆温强度
统计结果表明,贵阳地区近地空间近40年贴地逆温层的平均逆温强度为1.53℃/100m,中位数为0.62℃/100m,悬浮逆温层的平均逆温强度为0.97℃/100m,中位数0.65℃/100m,图10显示了逆温强度的年际变化及M-K检验,由此可以看到和逆温厚度特征一样,贴地逆温层整体逆温强度呈下降趋势,特别是在1990年代至2000年代下降趋势显著,2000年代以后则有缓慢上升的态势,对于悬浮逆温层而言,2010年代存在较为明显的突变,在此之前逆温强度波动幅度不大,在此之后则上升趋势较为明显,另外值得注意的是在1990年代以前,贴地逆温强度明显大于悬浮逆温,1990年代以后两者差距逐渐减小,并趋于一致。
(a) 不同类型逆温厚度年际变化
(b) 贴地逆温厚度M-K检验
(c) 悬浮逆温厚度M-K检验
Figure 10. Annual variation of different types of temperature inversion intensity and Mann-kendall mutation test in Guiyang from 1980 to 2019
图10. 1980年至2019年贵阳近地空间不同类型逆温强度年际变化及Mann-kendall突变检验
4.3. 逆温底高和逆温顶高
对于贴地逆温而言,逆温层底高即探空站海拔高度,统计近40年贵阳地区贴地逆温层顶高发现,其平均值为1557.7米,中位数为1421米,悬浮逆温层底高的平均值为2739.7米,中位数2569米,悬浮逆温顶高平均值3210.7米,中位数3075米,从图11的年际变化及M-K检验可以看到,贴地逆温层顶高在1990年代至2000年代呈上升趋势,2000年代以后下降趋势显著,对于悬浮逆温而言,逆温层底高整体呈上升趋势,而逆温层顶高却呈下降趋势,结合前文可知,这也是悬浮逆温层逆温厚度在近40年整体呈下降趋势的原因。
5. 逆温层成因初探
前文总结了贵阳机场近40年近地空间逆温时次及逆温指数的变化特征,形成逆温层有多种原因,包括辐射逆温、平流逆温、地形逆温、下沉逆温、锋面逆温和湍流逆温。辐射逆温是因地面的强烈辐射而失去热量,使紧贴其上的气层比上层空气有较大的降温,而形成温度随高度递增的现象,因此辐射逆温通常是造成贴地逆温的最主要原因。前文已经提到,贵阳地区贴地逆温层平均厚度为334米,中位数为199米,统计2010年至2019年贵阳地区辐射雾当日贵阳东山站08时探空数据发现,起雾日贴地层均出现了逆温现象,如图12(a)所示,辐射雾次数与总的逆温时数有着很好的相关性,两者相关系数计算所得为0.61,表现为显著性正相关,考虑到非辐射雾日近地面仍然可能存在辐射逆温现象,可以推论出贵阳地区贴地逆温现象主要以辐射逆温为主。图12(b)所示的2015年10月20日08时贵阳站探空图可作为贵阳地区贴地逆温层的温廓结构模型,表现出了浅薄逆温层的特点,此类逆温层通常在日出后随着太阳辐射增强而消散,具有明显的日变化特征,不易长时间维持。
(a) 2010~2019年
(b) 2015年10月20日08时贵阳站探空
Figure 12. (a) Comparison of radiation fog times and temperature inversion from surface times in Guiyang from 2010~2019; (b) Sounding map of Guiyang Station at 08:00:00 on October 20, 2015 (curve, black: temperature, green: dew point)
图12. (a) 2010年至2019年贵阳地区辐射雾次数与贴地逆温次数对比;(b) 2015年10月20日08时贵阳站探空(曲线,黑色:气温,绿色:露点)
由前文可知,对于贵阳地区而言,悬浮逆温是最主要的逆温现象,其发生频率超过了85%,张峻和王海军通过研究1960~1999年南方多个探空站数据发现,西南地区逆温强度较大,且多与低温雨雪冰冻天气相配合 [11],对于贵州地区而言,昆明准静止锋尤其活跃,特别是在冬季,2019年2月,贵州地区仅3日无准静止锋存在,选取准静止锋日ERA5高空要素场进行合成,如图13(a)所示,700 hPa贵州地区为一致的西南气流控制,风速达到了16 m/s及以上的急流强度,也是暖平流中心最大值集中区域所在,图13(b)显示了900 hPa的情况,贵州地区主要受来自偏东方向的冷平流所控制,这种底层冷上层暖的温廓结构在贵阳地区近地空间长时间维持,从图13(c)可以更加清晰看到冷暖平流在垂直方向上的配置,冷平流主要控制850 hPa及其以下区域,850 hPa至700 hPa则是暖平流集中的区域,图13(d)是2019年2月28日08时贵阳站探空图,可以看到逆温层从近地面的850 hPa附近一直延展到了800 hPa以上,结合前文所叙,这种特殊的冷暖结构冬季在贵州地区往往能长时间维持,一方面来自印度洋和西太平洋上的偏南暖湿气流持续从西部和南部进入贵州山区,另一方面冬季冷空气入侵频繁,但受地形阻挡主要以偏东北、偏东方向的回流形式进入贵州省,主要集中在近地层,且由于山区地形阻塞常常滞留不前,因此使得上暖下冷的底层结构稳定,有利于逆温层维持。
本文统计了贵阳地区近40年近地空间逆温层的年际变化和月变化,分析了其周期性特征和各项逆温指数的变化特征,就贵阳地区逆温层形成的原因进行了初步探讨,强调了辐射逆温在贴地逆温层和准静止锋锋面逆温在悬浮逆温中的贡献作用,但贵阳地区逆温层在气候尺度上的影响因子和变化原因则还有待于更进一步的研究。
6. 结论
本文利用1980~2019近40年探空数据对贵阳近地空间逆温层特征进行了统计分析,并对逆温层形成的原因进行了初步的探讨,得出以下结论:
1) 贵阳地区近地空间内逆温现象发生较为频繁,其全年发生频率为66%,08时出现频率高于20时,M-K检验显示近40年存在三个明显的突变时刻,1980至1990年代逆温发生时次处于上升趋势,1990年代至2000年代呈下降趋势,2000年以后又呈上升趋势。
2) 贵阳地区近地空间以悬浮逆温为主,占比超过了85%,贴地逆温多发生在08时,悬浮逆温在08时和20时发生时次相当,贴地逆温发生时次在近40年里整体呈下降态势,而悬浮逆温在1990年以前波动缓和,1990~2000年代处于下降趋势,2000年代以后有着明显的突变,突变以后呈上升趋势。
3) 贵阳地区逆温现象在冬春季的发生频率高于夏秋季,其中秋冬季发生时次在2000年代以后呈上升趋势,而春夏季在1990年中后期以后就呈下降趋势。
4) 贵阳地区贴地逆温时次在1~7月波动幅度不大,8~12月呈线性上升趋势;悬浮逆温全年则呈U型分布特征,11月至2月为发生逆温现象的高频月份,6~8月为发生逆温现象的波谷月份。
5) Morlet小波分析表明,贵阳地区贴地逆温时次存在着2~4年以及4~8年的年际变化周期,周期性特征在2010年之前较为显著,2010年之后不再显著;悬浮逆温时次则存在着2年左右和4~8年的年际变化周期以及8~16年的年代际变化周期,其中年际变化周期在2010年之前较为显著,2010年之后不再显著,年代际变化周期则始终较为显著。
6) 贵阳地区贴地逆温层的平均厚度为345.7米,平均逆温强度为1.53℃/100m,平均逆温顶高为1557.7米,悬浮逆温层的平均厚度为471米,平均逆温强度为0.97℃/100m,平均逆温底高为2739.7米,平均顶高为3210.7米;贴地逆温层在厚度、强度以及顶高均呈现出了整体下降的年变化趋势,悬浮逆温则厚度处于下降趋势,强度稳定并波动上升,底高整体处于上升趋势但顶高却处于下降趋势。
7) 贵阳地区近地空间的逆温现象以辐射逆温和准静止锋锋面逆温为主,辐射逆温多导致贴地逆温,而准静止锋锋面逆温多导致悬浮逆温。
NOTES
*第一作者。
#通讯作者。