1. 引言

全球变暖始终是受到世界关注的重要议题。近年来温室气体排放加剧，引起全球气温上升趋势不断增强，这一问题已经成为国际环境关注的焦点。2015年，全球众多国家齐聚法国巴黎，通过了划时代的《巴黎协定》。该协定旨在到本世纪末将地球升温控制在比工业化前水平高出不超过2℃，并争取进一步把增温限制在1.5℃以内，2℃和1.5℃温控目标迅速成为应对气候变化领域的热点话题，国际社会从各个领域和方向迅速开展行动 [1] 。自上个世纪以来，全球升温迅速。王绍武等学者探究发现1880~1991年间全球气候变暖的速度为0.5℃/100a [2] 。龚道溢和其他研究人员观察到，在过去的一百年间，世界范围内的平均气温呈上升状态，且自1980年代起温度上升速度明显加快 [3] 。唐国利等学者发现1910~2009年的变暖趋势为0.70℃~0.75℃/100a，截止2009年，已经变暖了0.8℃ [4] 。全球变暖是毋庸置疑的事实，但VV Klimenko、苏京志等学者认为海洋作为热量的主要吸收载体，使得全球变暖趋缓 [5] 。因此有必要厘清近年来全球气温变化的具体情况，明确全球升温趋势趋于变缓还是不断加速这一科学问题。

随着全球变暖加剧，极端气温事件越来越频繁。从1980年代末期起，东亚地区频遭暖冬侵袭，这一气候模式维持至2003年的冬天。此后东亚经历的寒冬事件逐渐增多，且寒冷程度逐渐加剧 [6] 。北半球中高纬度在21世纪初以来冷冬频繁 [7] ，对全球经济社会发展和民生造成严重影响。很多学者们认为由于全球气温上升，对大气环流发生正/负反馈作用，通过遥相关进一步影响其它地区的气候状况。如赵宗慈等学者发现若全球继续变暖，使得各个地区季风位置、范围和强度发生改变，从而导致极端降水时间频发 [8] 。单楷越等学者也发现人类温室气体排放是近年来台风频发的首要驱动因子 [9] 。因此本文想在此基础上，了解全球变暖对大气环流究竟产生何种程度的影响。综上，本文将使NCEP/NCAR再分析数据与部分中国科学院大气物理所发布的大气环流指数数据对近24a来全球气温的变化情况进行探究，并判断近年来全球变暖趋势是否趋于平缓，在此基础上分析全球气温变化与大气环流的相关程度，以期为学者们研究气候变化理论提供数据支撑。

2. 数据来源

2.1. NCEP/NCAR再分析资料

本文使用NCEP/NCAR再分析数据是发展最早且时序长的全球再分析资料。本文从Ncep网站注册下载了2000~2023年全球近地面2 m与500 hPa逐月平均气温数据，空间分辨率为1.875˚ × 1.875˚。

2.2. 逐月大气环流指数

气候系统监测指数集由中国科学院大气物理所公布。此数据集共包含130项气候系统指数。本文筛选主要大气环流指数20项，时间分辨率为月平均，起止时间为2000年1月~2023年12月。

3. 研究方法

3.1. 线性趋势法

采用线性趋势法对 [10] 分析2000~2023年全球气温的长期变化趋势，建立以时间(x)为自变量，气温(y)为因变量的一元线性回归方程：

(1)

式中：a为回归系数，正负代表气温随时间呈上升或下降趋势，用t检验检验方程是否具有显著性。

3.2. Mann-Kendall突变检验法

Mann-Kendall检验法是一种非参数检验方法，能够客观地表样本序列的整体变化趋势 [11] 。

对于具有n个样本量的时间序列x，构造秩序列：

$S_k={\displaystyle \underset{i=1}{\overset{k}{\sum }}{r}_i}\left(k=2,3,\cdots ,n\right)$ (2)

S_k是第i时刻气温数值大于第j时刻数值的累计数。在时间序列随机独立的假定下，定义统计量：

$U{F}_k=\frac{S_k-E\left(S_k\right)}{\sqrt{Var\left(S_k\right)}}\left(k=1,2,\cdots ,n\right)$ (3)

给定显著性水平a，根据正态分布表显示，若UF_k > U_a，则表明该序列存在明显的趋势变化特征。若UF_k < U_a，则表明该序列的趋势变化特征并不明显。按时间序列x逆序，再重复上述过程，同时满足UB_k = −UF_k，UB₁ = 0 $\left(k=1,2,\cdots ,n\right)$ 。

3.3. 相关分析

相关分析是一种研究现象间是否存在某种依存关系，并对具有依存关系的现象的相关方向及相关程度进行研究的理论和方法 [12] 。本文主要利用Pearson相关分析方法，对近24a来全球气温与主要大气环流进行相关分析。将r作为相关系数，用来表示最终的相关结果。计算公式如下：

$r=\frac{{\displaystyle {\sum }_{i=1}^n(x_i-\bar{x})(y_i-\bar{y})}}{\sqrt{{\displaystyle {\sum }_{i=1}^n{(x_i-\bar{x})}^2}\cdot {\displaystyle {\sum }_{i=1}^n{(y_i-\bar{y})}^2}}}$ (4)

r作为相关系数，其范围介于−1~1间。当全球气温与此大气环流的相关系数r > 0时，表明这二者呈正相关关系，在0~1范围内的r越大，代表研究区气温与主要大气环流的正相关关系越显著；当全球气温与此大气环流的相关系数r < 0时，表明这二者呈负相关关系，在−1~0范围内的r越小，代表研究区气温与主要大气环流的负相关关系越显著；当r = 0时，表明这全球气温与此大气环流呈无相关性的相互独立的关系。

4. 结果分析

4.1. 2000~2023年全球气温变化特征

4.1.1. 2000~2023年全球近地表气温变化特征

采用NCEP/NCAR再分析数据分析全球近地表年平均气温的变化特征(图1)。结果表明2000~2023年全球近地表年平均气温为7.589℃，呈显著升温趋势(P < 0.05)，变化速率为0.318℃/10a，2016年气温最高，为8.140℃，2000年气温最低，为6.997℃。R²为0.705，MK分析结果表明，全球近地表年均温近24a来未发生突变，表明总体升温趋势较为平缓。

4.1.2. 2000~2023年全球500 hPa气温变化特征

采用NCEP/NCAR再分析数据分析全球500 hPa年平均气温的变化特征(图1)。结果表明2000~2023年全球500 hPa年平均气温为−20.302℃，呈显著升温趋势(P < 0.05)，变化速率为0.265℃/10a，2016年气温最高，为−19.878℃，2000年气温最低，为−20.834℃。R²为0.634，MK分析结果表明，全球500 hPa年均温近24a来未发生突变，表明总体升温趋势较为平缓。

4.2. 2013~2023年全球气温变化特征

4.2.1. 2013~2023年全球近地表气温变化特征

采用NCEP/NCAR再分析数据分析全球近地表年平均气温的变化特征。结果表明2013~2023年全球近地表年平均气温为7.799℃，呈显著升温趋势(P < 0.05)，变化速率为0.343℃/10a，2016年气温最高，为8.140℃，2021年气温最低，为7.751℃。R²为0.343，MK分析结果表明，全球近地表年均温近10a来未发生突变，表明总体升温趋势较为平缓。相比于2000~2023年，近十年北极近地表平均气温高出0.210℃，变化速率超出0.025℃/10a，证明全球变暖趋势仍不断增强。

4.2.2. 2013~2023年全球500 hPa气温变化特征

采用NCEP/NCAR再分析数据分析全球500 hPa年平均气温的变化特征。结果表明2013~2023年全球500 hPa年平均气温为−20.142℃，呈显著升温趋势(P < 0.05)，变化速率为0.280℃/10a，2016年气温最高，为−19.878℃，2013年气温最低，为−20.375℃。R²为0.315，MK分析结果表明，全球500 hPa年均温近10a来未发生突变，表明总体升温趋势较为平缓。相比于2000~2023年，近十年北极500 hPa平均气温高出0.160℃，变化速率超出0.015℃/10a，证明全球500 hPa气温上升趋势仍不断加强，但相比于近地表，全球500 hPa升温趋势相对较缓。

4.3. 2000~2023年全球气温变化与大气环流的相关性

4.3.1. 全球近地表气温变化与大气环流的相关性

计算2000~2023年全球近地表年均温与20项主要大气环流指数相关系数，统计结果如表1所示。全球近地表年均温与编号为1~6的大气环流指数呈极显著正相关关系，与编号为8与18的大气环流指数呈显著正相关关系，相关系数介于0.442~0.884间；全球近地表年均温与编号为12、13与15的大气环流指数呈极显著负相关关系，与编号为11与14的大气环流指数呈显著负相关关系，相关系数介于−0.704~−0.457间。

注：^*、^**表示在0.05、0.01概率水平下检验。

4.3.2. 全球500 hPa气温变化与大气环流的相关性

计算2000~2023年全球500 hPa年均温与20项主要大气环流指数相关系数，统计结果如表2所示。全球500 hPa年均温与编号为1~6的大气环流指数呈极显著正相关关系，与编号为8、9与18的大气环流指数呈显著正相关关系，相关系数介于0.438~0.921间；全球500 hPa年均温与编号为12、13与15的大气环流指数呈极显著负相关关系，与编号为11的大气环流指数呈显著负相关关系，相关系数介于−0.667~−0.459间。

注：^*、^**表示在0.05、0.01概率水平下检验。

5. 结论

1) 2000~2023年全球近地表年平均气温为7.589℃，介于6.997℃~8.140℃间，变化速率为0.318℃/10a。MK分析结果表明，全球近地表年均温近24a来未发生突变。近10a全球近地表气温比2000~2023年全球近地表平均气温高出0.210℃，变化速率超出0.025℃/10a，证明近年来全球变暖趋势不断加强。

2) 2000~2023年全球500 hPa年平均气温为−20.302℃，介于−20.834℃~−19.878℃间，变化速率为0.265℃/10a。MK分析结果表明，全球500 hPa年均温近24a来未发生突变。近10a全球500 hPa气温比2000~2023年全球500 hPa平均气温高出0.160℃，变化速率超出0.015℃/10a，证明近年来全球500 hPa变暖程度逐渐加剧，但相比于近地表，升温幅度相对较小。

3) 近24年来，全球近地表与500 hPa最高气温均出现在2016年，最低气温均出现在2000年。年平均气温相差较大，为27.891℃；升温速率相差0.053℃/a，近地表气温的升温速率是500 hPa的1.2倍。MK分析结果表明，近24a来全球近地表与500 hPa年均温均未发生突变。

4) 全球近地表年均温与北半球、西太平洋和北太平洋的副高面积指数，北半球、西太平洋和太平洋的副高强度指数，西太平洋副高北界位置指数，欧亚经向环流指数呈正相关关系；与北半球和亚洲区的极涡面积指数，亚洲区、北半球和太平洋区极涡强度指数呈负相关关系。

5) 全球500 hPa年均温与北半球、西太平洋和北太平洋的副高面积指数，北半球、西太平洋和太平洋的副高强度指数，西太平洋和北太平洋的副高北界位置指数，欧亚经向环流指数呈正相关关系；与北半球和亚洲区的极涡强度指数、北半球和亚洲区的极涡面积指数呈负相关关系。

参考文献