环境科学  2023, Vol. 44 Issue (5): 3003-3016   PDF    
1950~2019年中国季节平均最高气温时空演变特征及其大气环流影响定量化分析
苏越1, 路春燕1,2,3, 黄雨菲1,2, 苏艳琳1,3, 王自立1,3, 雷依凡1     
1. 福建农林大学计算机与信息学院, 福州 350002;
2. 福建省资源环境空间信息统计研究中心, 福州 350002;
3. 生态与资源统计福建省高校重点实验室, 福州 350002
摘要: 全球变暖和人类活动加剧引起区域气候不稳定,高温事件发生频次和持续性增强,生态环境保护和经济社会发展面临严峻威胁.采用Mann-Kendall检验法和一元线性回归法探讨1950~2019年中国季节平均最高气温(AMT)时空演变特征,并结合线性相关、偏线性相关以及小波分析揭示中国各季节AMT与大气环流间的联系.结果表明:①时间尺度上,各季节AMT均呈显著上升趋势,春、夏、秋和冬季AMT上升幅度分别为1.21、0.08、1.81和0.25℃,其突变时间主要集中于20世纪90年代至21世纪初.②空间分布上,除夏季以外,其它季节AMT平均变化率由南向北逐渐增加,但其增加程度各不相同,其中春冬季节东北和西北地区变化速率最快.③各季节AMT与大气环流因子间存在复杂关系,其能量强度在不同时频域上具有显著差异.其中太平洋年代际涛动与夏季AMT呈显著负相关,北大西洋涛动对夏秋冬季AMT变化产生较大程度正向驱动,北极振荡对各季节AMT均具有明显正向驱动作用,厄尔尼诺-南方涛动因年份不同对春夏季AMT的短期变化存在显著正向或负向影响.研究结果可为我国制定科学有效的气候变化应对方案提供理论依据和技术参考.
关键词: 季节平均最高气温(AMT)      大气环流      Mann-Kendall检验      交叉小波      小波相干      遥相关分析     
Quantitative Analysis of Spatio-temporal Evolution Characteristics of Seasonal Average Maximum Temperature and Its Influence by Atmospheric Circulation in China from 1950 to 2019
Su Yue1 , Lu Chun-yan1,2,3 , Huang Yu-fei1,2 , Su Yan-lin1,3 , WANG Zi-li1,3 , LEI Yi-fan1     
1. College of Computer and Information Sciences, Fujian Agriculture and Forestry University, Fuzhou 350002, China;
2. Research Centre of Resource and Environment Spatial Information Statistics of Fujian Province, Fuzhou 350002, China;
3. Key Laboratory of Ecology and Resources Statistics of Fujian Province Universities, Fuzhou 350002, China
Abstract: Global warming and intensified human activities have led to regional climate instability with increasing frequency and the persistence of high-temperature climate events. Eco-environmental protection and socio-economic development have been faced with rigorous threats. Taking the monthly maximum temperatures from 1950 to 2019 as the basic data source, the spatial-temporal evolution characteristics of seasonal average maximum temperature (AMT) were discerned using the Mann-Kendall test and unary linear regression method in China from 1950 to 2019. Combined with linear correlation, partial linear correlation, and wavelet analysis, the correlation between seasonal AMT characteristics and atmospheric circulations was analyzed quantitatively. The results showed that: ① the AMT in all seasons had a significant upward trend, with an increase of 1.21, 0.08, 1.81, and 0.25℃ in spring, summer, autumn, and winter, respectively. The abrupt change times of the AMT were concentrated in the 1990s to the early 21st century. ② In terms of spatial distribution, except for in summer, the average trend rates of AMT in other seasons increased gradually from south to north, although the increasing degrees were different. Among them, the AMT change rate in spring-winter was the fastest in northeast and northwest China. ③ There were complex correlations between the AMT of every season and atmospheric circulation factors, and the distribution of the interrelation energy varied significantly in different frequency domains. Specifically, the Pacific Decadal Oscillation had a significant negative correlation with AMT in summer. The North Atlantic Oscillation had an active effect on AMT changes in summer, autumn, and winter. The Arctic Oscillation had a significant positive driving effect on AMT in all seasons, and there were significant positive or negative influences on the short-or long-term changes of AMT in spring and summer due to the different EI Niño-Southern Oscillation years. These results could provide a theoretical basis and technical reference for China to formulate scientific and effective response plans of climate change.
Key words: seasonal average maximum temperature(AMT)      atmospheric circulation      Mann-Kendall test      cross wavelet      wavelet coherence      tele-correlation analysis     

在全球气候变化和人类活动强度增加的双重驱动下, 高温天气频发且持续性增强, 生态系统平衡和社会经济发展均面临前所未有的危机和挑战, 高温气象灾害已成为世界各国共同关注的问题之一[1].联合国政府间气候变化专门委员会(Intergovernmental Panel on Climate Change, IPCC)在第五次气象评估报告指出, 近百年来全球气温呈现显著变暖趋势, 持续性高温事件所导致的干旱、森林火灾等气象灾害发生频次显著增加[2, 3].2006年中国川渝地区[4]和2018年欧洲西、北部[5]均遭遇严重的高温伏旱, 致使当地农作物严重受损.同时, 持续性高温可通过降低植被的含水量增大森林大火的可燃性和传播率等, 导致火灾快速升级至极端水平, 如2010年我国东北地区[6]和2019年澳大利亚[7], 因高温引起的持续性干旱造成特大森林火灾.此外, 高温天气频发对人类身体健康与物质财产安全亦产生多种不利影响, 一方面, 高温不仅影响人体舒适度, 导致工作效率降低, 同时易诱发多种生理疾病, 如中暑、心脑血管疾病等, 严重可致人死亡[8, 9].另一方面, 高温事件易对社会经济造成严重损失, 如2003年和2010年欧洲发生高温热灾事件, 导致数十亿美金的经济损失[10]; 1988~2004年高温引起的旱灾致使我国经济损失近756.88亿元[11].因此, 明确高温变化规律对于维持生态系统平衡与人类社会经济可持续发展均具有重要意义.

针对高温变化问题, 国内外学者已展开了研究, 并取得众多成果.早期研究主要采用气象站点数据, 多关注于最高气温变化趋势的分析, 注重在时间尺度上量化高温演变过程.任国玉等[12]利用中国600多个气象站资料分析1951~2004年全国气候变化趋势, 结果指出我国最高气温不断上升, 夏季炎热期延长; 陈建宇等[13]利用3个气象站点资料研究陕西榆林地区1954~2013年最高气温年际变化趋势, 结果表明该区最高气温显著上升; Gocic等[14]对1980~2010年塞尔维亚地区各气象变量的年际变化趋势进行分析, 发现该区最高气温呈显著上升趋势.随着研究不断深入, 在“时”与“空”多尺度分析的基础上, 高温天气变化驱动因素的研究日益细致, 逐渐探究其与大气环流间的联系, 如北大西洋涛动[15](North Atlantic Oscillation, NAO)、太平洋年代际涛动[16](Pacific Decadal Oscillation, PDO)、北极振荡[17](Arctic Cscillation, AO)和厄尔尼诺-南方涛动[18](EI Nin~o-Southern Oscillation, ENSO)等大尺度大气环流因子已相继被研究.此外, 随着空间技术不断发展, 已有学者采用空间降尺度方法对低分辨率的网格数据集进行处理, 以获得高分辨率的网格气候数据集研究区域尺度的气候特征[19].尽管以往对于高温变化驱动因素的研究取得了众多成果, 但研究多采用气象站点数据[12, 13, 17], 主要关注一个或两个大气环流因子与高温间的相关关系[15~17], 且在时间尺度上多以年计量[13, 17, 20], 而针对较小时间尺度的研究较少, 由此使得研究结果往往存在片面性, 不易发现细微差异; 同时, 截至目前针对我国全国范围最高气温变化趋势及其影响因素定量化分析相对较少, 且大气环流对我国最高气温的影响尚不完全清楚.

就研究方法而言, 目前专家学者多采用M-K检验[17, 21]、累积距平[22, 23]、线性回归[13, 24]、滑动T检验[25, 26]、R/S方法[20, 26]和小波分析法[27, 28]等对区域气候变化进行研究, 其中M-K检验具有研究样本无需服从正态分布, 不易受少数异常值干扰, 且其计算过程简便, 检验范围广的优势, 更适用于研究气候变化的趋势及突变情况[22, 29]; 小波分析法能够分析时频域中两序列多时间尺度遥相关关系, 并揭示两序列在高/低能量区的位相关系及其周期性特征, 更适用于分析气候变化与大气环流间的联系[17, 28], 但综合利用M-K检验和小波分析法进行气候变化的研究却少有报道.

鉴于以上, 本研究利用空间降尺度方法生成的高分辨率气温数据集, 在分析1950~2019年我国季节平均最高气温(average maximum temperature, AMT)时空变化趋势及突变特征的基础上, 综合利用交叉小波(cross wavelet transform, XWT)和小波相干(wavelet coherence, WTC)等多种方法探究各季节AMT与典型大气环流因子间的多尺度相关关系, 以期为中国高温气象灾害评估和预测提供科学参考和理论依据.

1 材料与方法 1.1 研究区概况

我国位于东半球北部, 介于73°33′E~135°05′E, 3°51′N~53°33′N间, 国土面积约9.60×106 km2, 地貌类型多样, 地势西高东低, 呈三级阶梯状分布.境内山脉众多, 按其走向可分为东西走向、东北-西南走向、西北-东南走向、南北走向和弧形山系等5种类型.整体上, 我国经度与纬度跨度广, 加之地势高低不同, 地形及山脉走向多样, 形成了复杂的气候类型[30].统计数据显示, 2019年我国平均气温为10.34℃, 较常年同期升高0.79℃, 且有64个站点日最高气温突破历史极值, 并发生不同程度的高温气象灾害(http://www.cma.gov.cn/).

1.2 数据来源及预处理

1950~2019年中国逐月最高气温数据来自国家地球系统科学数据中心(http://www.geodata.cn), 该数据集空间分辨率为0.008 333 3°(约1 km), 数据的绝对误差为1.28℃.该数据具有时间跨度长、数据完整性好和可靠性强等特点, 地理范围涵盖除南海岛礁与岛屿外的中国陆地范围.通过对月度最高气温进行均值处理得到各季节AMT数据.需说明的是, 本研究根据国家气象划分方法对季节进行划分, 即:3~5月为春季, 6~8月为夏季, 9~11月为秋季, 12月至翌年2月为冬季.

为探究各季节AMT与大气环流间的耦合关系, 本研究选取研究时段内4种最具代表性的大尺度大气环流因子进行分析, 即北极振荡(AO)、北大西洋涛动(NAO)、太平洋年代际涛动(PDO)和厄尔尼诺-南方涛动(ENSO), 各大气环流因子均以指数形式表示, 逐月AO、NAO、PDO和ENSO数据均来源于美国国家气象局气候预报中心(https://www.cpc.ncep.noaa.gov/products/precip/CWlink/MJO/climwx.shtml), 并通过均值处理得到各季节4个大气环流因子数据, 大气环流因子数据的详细信息见表 1.

表 1 各大气环流因子数据详细信息 Table 1 Detailed data information of atmospheric circulation factors

1.3 研究方法 1.3.1 Mann-Kendall检验

为明确中国季节AMT的时间变化趋势及突变特征, 本研究采用M-K检验分析各季节AMT的变化率与突变点.M-K检验包括M-K趋势检验和M-K突变检验.假设某要素的时间序列为x1, x2, …, xn, 则M-K趋势检验的基本原理可表示为[29]

(1)
(2)
(3)
(4)
(5)

式中, S服从正态分布, sgn为符号函数, β为斜率.在给定的α置信水平上, 如果|Zc| ≥Z(1-α/2), 则气象因素的变化趋势显著.如果β>0, 则表示气候要素呈上升趋势, 反之亦然.

M-K突变检验[31]的具体计算过程为:

(6)
(7)

在此基础上, 若原序列随机且独立, 则定义统计量UFk

(8)
(9)
(10)

引用公式(6)~(10)到反序列xn, xn-1, …, x1中, 同时使UBk=-UFk(k=n, n-1, …, 1), 若UFk曲线和UBk曲线出现交点, 且交点位于置信区间内, 则该点为突变点.

1.3.2 年际变化趋势分析

一元回归分析模型[32]可通过回归方程的斜率(Slope)反映空间要素在时间尺度的变化趋势和强度.在本研究中利用一元回归分析模型计算1950~2019年研究区各季节AMT的年际平均变化趋势, 其计算公式如下:

(11)

式中, n为所研究时间序列长度; Pi为第i年某季节AMT.

1.3.3 相关分析

为客观地量化各大气环流因子与各季节AMT的相关程度, 本研究综合利用线性相关系数(即Pearson相关系数)[33]与偏线性相关系数[34]进行分析, 其计算公式如下:

(12)

式中, n为所研究时间序列长度; h为第h年; xi为变量xi的样本均值; xj为变量xj的样本均值.

(13)

式中, rij.hml为排除变量xhxmxl后, 变量xixj的偏相关系数.

1.3.4 小波分析

大气环流是影响区域气候变化的重要因素, 其变化与异常波动在空间上可同时或先/后对距离遥远的区域产生作用, 由此产生的相关性称为遥相关[35].基于小波变换与交叉谱构建的交叉小波(XWT)方法能够分析时频域中两序列多时间尺度遥相关关系, 揭示两序列在高能量区的位相关系[27].具体原理为:假设两时间序列分别为X={xi|i=1, 2, …, n}与Y={yi|i=1, 2, …, n}, 连续小波变换分别为WnX(s)和WnY(s), 则交叉小波谱可定义为:

(14)

式中, WnY*(s)为WnY(s)的复共轭; 交叉小波功率谱为WnXY(s)的绝对值.

尽管交叉小波(XWT)可揭示两序列在高能量区的位相关系, 但对低能量区位相关系解析不足[27].为此, 本研究采用小波相干(WTC)对各季节AMT与大气环流因子间的遥相关关系加以分析.具体定义如下:

(15)

式中, |S[s-1WXY(s)]|2为两时间序列在某一频率下波动振幅的交叉积; S[s-1 |WnX(s)| 2]和S[s-1 |WnY(s)| 2] 为两时间序列振动波的振幅.

2 结果与分析 2.1 季节平均最高气温变化趋势及突变特征

基于M-K检验得到1950~2019年中国地区各季节AMT变化趋势及突变情况(表 2图 1).由表 2可知:在α=0.05的显著性水平下, 近70年中国各季节AMT均呈显著上升趋势.其中, 近70年各季节AMT趋势倾向率按降序排列分别为:春季[0.19℃ ·(10 a)-1]、秋季[0.10℃ ·(10 a)-1]、夏季[0.08℃ ·(10 a)-1]和冬季[0.02℃ ·(10 a)-1].由图 1可知:尽管各季节AMT总体均呈升温态势, 但波动振荡现象亦与之共存, 且发生突变年份各异.春季、夏季、秋季和冬季AMT分别在2002及2003、2005及2008、1994和1985年产生突变, 且均为增温突变.总体而言, 中国各季节AMT发生增温突变的时间略有差异, 但突变时间均集中于20世纪末至21世纪初, 且在突变后均呈显著增温趋势.

表 2 1950~2019年各季节平均最高气温变化趋势及其显著性检验 Table 2 Trend analysis and significance test of seasonal average maximum temperature from 1950 to 2019

图 1 1950~2019年各季节平均最高气温突变M-K曲线 Fig. 1 Mutation M-K curves of seasonal average maximum temperature from 1950 to 2019

2.2 季节平均最高气温空间演变特征

通过年际平均变化趋势分析得到基于像元的各季节AMT变化率.由图 2可知, 1950~2019年各季节AMT变化趋势存在较大程度的空间异质性.整体上, 除夏季以外, 其它季节AMT的年际平均倾向率由南向北均逐渐增加, 各季节AMT的变化趋势各不相同, 其中西北和东北地区增温幅度较大, 降温区域主要集中在华东、华中、华南和西南地区.就年际平均倾向率变化范围而言, 近70年春季AMT的平均变化范围最大[-0.84~1.12℃ ·(10 a)-1], 其次为秋季[-0.45~0.81℃ ·(10 a)-1], 再次为冬季[-0.30~0.90℃ ·(10 a)-1], 夏季变化范围最小[-0.14~0.82℃ ·(10 a)-1].

基于自然资源部标准地图服务网站GS(2020)4617号标准地图制作, 底图无修改; 饼图表示各季节AMT年际平均倾向率占研究区总面积的百分数 图 2 1950~2019年各季节平均最高气温变化趋势空间分布 Fig. 2 Spatial distribution of the tendency rate for seasonal average maximum temperature from 1950 to 2019

空间上, 春季AMT的升温速率由西南向东北逐渐增加, 且增温区域约占我国陆地总面积的91.1%, 主要分布在西北、东北、华北和华东地区, 而降温区域主要集中在青藏高原和云贵高原区域.夏季AMT的平均倾向率由中部向四周逐渐增加, 且绝大多数地区的平均变化速率处于-0.15~0.15℃ ·(10 a)-1之间, 变化幅度相对较小.秋季AMT升温幅度大于0.30℃ ·(10 a)-1的区域主要集中在塔里木盆地(约占我国陆地总面积的11.1%), 降温区域主要集中在青藏高原.冬季AMT的年际平均倾向率空间分布与春季相似, 由南向北逐渐增加, 其中平均变化速率大于0.30℃ ·(10 a)-1地区约占我国陆地总面积的48.5%, 主要集中在西北、华北和东北地区, 增幅相对较大.

2.3 季节平均最高气温与大气环流因子相关关系

利用线性相关分析和偏线性相关分析计算可知(表 3), 在α=0.01或α=0.05的显著性水平下, 各季节AMT均至少与一个大气环流因子呈显著线性或偏线性相关, 即春、秋和冬季AMT均与AO呈显著正相关, 夏季AMT与PDO呈显著负相关.

表 3 各季节平均最高气温与大气环流因子间的相关关系1) Table 3 Correlativity between seasonal average maximum temperature and atmospheric circulation factors

2.4 季节平均最高气温与大气环流因子遥相关关系

基于交叉小波和小波相干分析得到各季节AMT与大气环流因子间的遥相关关系结果, 如图 3所示.由图 3(a)可知, 春季AMT与AO、ENSO和PDO在高/低能量区均存在显著共振阶段(主要集中在2.07~4.38 a周期), 其中, 春季AMT与AO在高能量区存在4个显著共振[图 3(a1)], 即1954~1961、1966~1973、1978~1990和2006~2014年, 且在各共振阶段二者对应频段位相关系均呈正相关; 与ENSO在高能量区存在1968~1973、1978~1988和2009~2012年3个显著共振, 第一阶段春季AMT对应频段位相超前ENSO变化约1/4周期(约3个月), 后两个阶段二者对应频段位相关系呈负相关[图 3(a2)]; 与PDO在低能量区存在1979~2000和2004~2014年显著共振阶段, 且两阶段二者对应频段位相关系呈正/负相关[图 3(a8)].

U型黑色细实线为小波边界效应影响锥, 在该曲线外的功率由于受到边界效应的影响而不予考虑; 粗黑实线表示通过置信水平为95%的红噪声检验, 两者相关性显著; ←表示AMT与大气环流因子间变化位相相反, 即两者呈负相关关系, →表示两者变化位相一致, 即两者呈正相关关系, ↑表示AMT落后大气环流因子变化1/4周期, ↓表示AMT超前大气环流因子变化1/4周期; 色柱表示能量密度的相对变化; XWT:AMT(春季)-AO表示春季AMT与AO的交叉小波谱, WTC:AMT(春季)-AO表示春季AMT与AO的小波相干谱, 以此类推 图 3 各季节平均最高气温与各大气环流因子的交叉小波谱和小波相干谱 Fig. 3 Cross wavelet spectrum and wavelet coherence spectrum of seasonal average maximum temperature and atmospheric circulation factors

夏季AMT与各大气环流因子具有不同的共振周期性特征.夏季AMT与ENSO在高/低能量区分别存在1988~1999年(以2.32~3.68 a为周期)和2008~2016年(以2.32~3.90 a为周期)显著共振[图 3(b2)图 3(b6)].夏季AMT与PDO在高能量区存在3个间歇性共振阶段, 即1953~1957、1957~1978和1989~1998年, 各阶段夏季AMT与PDO对应频段位相关系均呈负相关[图 3(b4)].在低能量区1961~1984年(以7.80~11.69 a为周期)夏季AMT与AO存在显著共振, 位相关系表现为夏季AMT超前AO变化约3个月[图 3(b5)].在低能量区2008~2012年(以2.07~2.92 a为周期)夏季AMT与NAO存在显著共振, 且二者对应频段位相关系呈正相关[图 3(b7)].

秋季AMT与各大气环流因子在不同阶段存在显著的年际和年代际尺度共振周期.在高能量区1980~1996年(以2.76~6.19 a为周期)秋季AMT与AO存在显著共振, 其能量谱呈反“L”状分布特征; 在低能量区存在1966~2000年(以12.38~22.38 a为周期)显著共振, 且二者在高/低能量区对应频段位相关系均呈正相关[图 3(c1)图 3(c5)].秋季AMT与ENSO、PDO在高能量区分别存在1979~1997年(以2.32~6.19 a为周期)和1984~1998年(以3.90~6.19 a为周期)显著共振[图 3(c2)图 3(c4)].秋季AMT与NAO在高/低能量区1980~1996年(以2.07~8.26 a为周期)和1966~2003年(以12.38~14.73 a为周期)分别存在显著共振, 且二者在高/低能量区的对应频段位相关系均呈正相关[图 3(c3)图 3(c7)].

冬季AMT与各大气环流因子在不同时频域二者间的相关性与周期性存在差异.冬季AMT与AO在高/低能量区1970~1989年(以6.95~9.28 a为周期)和1963~1991年(以4.91~11.03 a为周期)分别存在显著共振, 两个阶段冬季AMT对应频段位相关系均超前AO变化约3个月[图 3(d1)图 3(d5)].冬季AMT与PDO在高能量区存在2个显著共振阶段, 即1975~1979年(以2.46~3.68 a为周期)和1997~2003年(以2.19~2.76 a为周期), 且各阶段二者对应频段位相关系均呈负相关; 在低能量区二者存在3个显著共振阶段, 即1976~1987、1977~1980和1995~2000年[图 3(d4)图 3(d8)].

总体上, 各季节AMT与各大气环流因子均存在显著共振, 在不同时频域上相关程度与共振周期存在差异, 其中春季AMT与各大气环流因子的联系均体现在较短年代际周期尺度上, 且与AO主要呈正相关; 与ENSO呈负相关或超前大气环流因子变化; 与PDO在短/长期年际周期尺度呈负/正相关; 而NAO与春季AMT关系较弱.夏季AMT与AO、ENSO、NAO呈正相关或超前大气环流因子变化, 而与PDO表现为负相关或滞后其变化.秋季AMT受到AO和NAO较长年代际周期变化的正相关影响, 而ENSO和PDO主要影响秋季AMT短期年际周期特征, 且相关关系相对较弱.冬季AMT与AO、NAO均存在显著共振阶段, 且对应频段位相关系呈现冬季AMT超前AO和NAO变化约3个月, 而ENSO和PDO与冬季AMT年际变化关系较弱.

3 讨论 3.1 中国季节平均最高气温影响因素分析

在人类活动强度增强和全球气候变化加剧的背景下, 高温天气已成为各国政府和专家学者共同担忧的问题之一.为此, 结合我国季节AMT时空变化趋势以及与大气环流因子间的耦合关系(图 1~3表 2表 3), 本研究对影响我国季节AMT变化的驱动因素进行探讨分析.总体上, 我国季节AMT受到长期或短期、直接或间接的人为和自然因素的影响(图 4), 明确驱动力因素可有助于确定未来研究工作的重心.

图 4 大气环流系统与人类活动影响中国平均最高气温变化的理论框架 Fig. 4 Theoretical framework of atmospheric circulation system and human activities influencing average maximum temperature change in China

3.1.1 气候系统内部变率作用

大气环流因子对区域气候具有重要影响, 在空间上不仅可作用于大气环流周边地区, 亦可以遥相关的形式对距离遥远的地区产生影响[36].相关研究表明, AO作为气候系统中一种重要的内部变率和动力过程, 主要影响北极地区与北半球中高纬度地区间气压的反向波动[37].冬季AO通过影响西伯利亚高压和东亚大槽改变冬季风, 并与北半球中高纬度地面气温呈显著正相关, 当AO为正相位增强时, 我国北方和西北地区易出现暖冬[38].本研究结果显示(图 2表 3), 冬季AMT与AO呈显著正相关, 且年际平均倾向率高于0.30℃ ·(10 a)-1主要出现在北方和西北地区.

NAO被认为是AO在北大西洋区域的一种表现形式, 具有显著的年代际共振周期[39].范可等[40]研究发现, 当NAO处于正相位时, 我国东北、华北和华南等地显著增温, 西南和青藏高原地区显著降温.龚道溢等[41]研究指出, 冬季NAO位相在20世纪80年代中期出现转折, 并造成气候系统出现年代际变化.在本研究中, 冬季AMT在1985年产生增温突变, 且与NAO存在年代际尺度显著共振周期(图 1图 3), 这与以往研究结果相吻合.

ENSO所引起的赤道东太平洋地区的风场和海面温度振荡是全球海洋和大气相互作用最强的信号[42].当发生ENSO事件时, 我国华北地区易产生高温, 东北地区易出现冷夏; 同时, 江淮、江南和华南等地气候变化也因受ENSO事件影响而发生气温突变[29].相关资料显示[43, 44], 在1986~1988、1994~1995、1997~1998、2002~2003、2004~2005和2006~2007年均有ENSO现象发生.本研究中各季节AMT发生突变的时间均与以上厄尔尼诺现象发生时间相一致(图 1).

PDO是一种类似于ENSO的太平洋气候变动过程, 具有显著的年际和年代际振荡周期, 并对东亚大气环流产生重要影响[45].沈皓俊等[46]研究发现PDO在2000年左右发生显著突变; 梁苏洁等[47]研究发现在1978~2005年冬季太平洋北部海温升高, 表现为PDO正相位.徐忆菲等[48]研究指出, 在年际尺度上, PDO超前于我国冬季气温变化, 且当PDO为正相位时, 我国冬季气温普遍偏高.在本研究中, 我国各季节AMT增温突变时间主要集中于20世纪90年代末至21世纪初, 并在之后呈显著增温态势(图 1); 冬季AMT与PDO在1995~2000年(以2.07~3.10 a为周期)存在显著共振, 对应频段位相表现为冬季AMT滞后PDO变化约3个月(图 3).

综上所述, 尽管各大气环流因子对高温天气作用方式和产生效果各不相同, 但气候系统内部变率对我国各季节AMT的变化趋势与突变特征均具有重要影响.

3.1.2 人类活动的影响

随着人类社会的发展, 人口急剧增加与经济社会快速发展使得我国在城市化过程、下垫面改变和能源消耗等方面发生结构性变化[49, 50], 这些因素均直接或间接地影响我国季节AMT变化.对我国1949~2019年人口数量和能源消耗总量进行分析可知(图 5), 二者在近70年均呈逐年上升趋势.人口增长与经济发展, 一方面使得我国下垫面不断改变, 从而影响区域高温天气.自改革开放以来, 我国沿海湿地因城市化进程加快、人类活动剧增而遭受巨大损失, 填海造地使得沿海湿地面积累计损失约占我国沿海湿地总面积的51.2%, 其中我国北部沿海地区填海造地超过填海造地总面积的70%[51].于竹筱等[52]研究1990~2018年我国林地的时空变化特征, 结果发现我国林地减少主要发生在东北地区、西南地区和新疆北部.本研究结果显示在1950~2019年我国各季节AMT均呈持续上升趋势且存在较大程度的空间异质性(表 2图 2), 其中东北地区、新疆北部和北部沿海地区升温速率较快, 可见我国自然植被破坏区域与AMT升高区域空间位置相一致.另一方面, 人口数量递增, 加之科学技术发展和生产规模迅速扩增, 在一定程度促进能源(如石油、天然气等)消耗[53], 由此使得大气中二氧化碳浓度相应增高; 加之人为活动对自然植被的破坏降低了它们吸收二氧化碳的能力[54].过量的二氧化碳能够强烈吸收地面辐射中的红外线, 致使区域气温不断上升, 产生“温室效应”.Skytt等[55]研究指出, 冬季供暖会增加以甲烷为主的天然气消耗, 产生更多温室气体, 从而使得冬季显著变暖.同时, 在社会生产活动中过度使用氟氯化碳等化学用品会产生大量氮氧化物, 导致空气污染恶化, 严重破坏臭氧层, 从而使得低层大气变暖, 高层大气变冷, 加剧区域气温升高[56].本研究中, 春冬季AMT在北方和西北地区增温速率较快(图 2), 这与我国北方和西北地区集体供暖关系密切.

数据来源于中国统计年鉴, 人口数量为1949~2019年数据, 标准煤消费总量为1949~2019年数据 图 5 中国地区总人口数量和标准煤消费总量变化 Fig. 5 Changes in total population and energy consumption in China

此外, 随着自然环境保护重要性逐步被政府和民众认知, 我国先后启动了多项重要生态系统保护和修复工程, 如“三北森林保护工程”、“退耕还林、封山绿化”和“草原生态保护计划”等[57].Zheng等[58]研究指出, 由于“三北”防护林和退耕还林等沙地荒漠化防治工程的实施, 使得我国毛乌素沙地植被总体得到改善, 该地区近地表温度平均降低0.12~0.32℃; Yu等[59]利用陆地-大气耦合区域气候模式评估2001~2017年中国植被绿化对气候的影响, 结果发现我国大部分地区, 尤其是华东地区, 夏季绿化趋势显著.本研究结果显示(图 2), 华东、华南和西南地区各季节AMT年际平均倾向率相较其他地区更为平缓, 夏季东北、西北和华北东部地区与其他季节相比, 其AMT年际平均倾向率呈减缓趋势.由此可见, 植被修复和绿化工程的实施有效提高了我国植被覆盖率, 生态系统碳循环得以良性促进, 有助于减缓气温升高.

由以上可知, 我国各季节AMT时空演变受到长期或短期以及直接或间接的人类活动影响.人类活动具有双重效应, 一方面, 人口与能源消耗增加、城市快速扩张等因素加剧我国季节AMT增温和突变; 另一方面, 生态保护工程实施对气温升高起到一定的缓解作用.

3.2 与以往研究结果的对比分析

相较于以往研究多以年为计量尺度, 本研究在季节尺度上探讨我国AMT时空演变特征, 并结合多个大气环流因子定量化分析其与我国季节AMT变化的联系, 较先前仅从气象要素本身或结合单个大气环流因子分析区域高温变化特征更为系统和准确.同时, 本研究从全国尺度刻画长时间序列的季节AMT时空演变特征, 拓展了以往区域性空间研究的尺度, 可更全面地了解我国高温整体变化情况.此外, 本研究结果显示各季节AMT均呈显著升温态势[0.02~0.19℃ ·(10 a)-1], 其突变时间主要集中在20世纪90年代至21世纪初, 且与AO和PDO的相关性较强、与ENSO和NAO的相关性较弱, 这与已有研究的结果相一致[22, 26, 60, 61].对比以往针对区域的研究, 由于其区域范围、时间尺度和数据来源与本研究存在不同, 其研究结果与本研究略有差异, 例如:史雯雨等[20]基于金沙江流域39个气象站点1960~2016年逐日平均最高气温数据分析发现该区季节平均最高气温均呈波动上升趋势[0.08~0.25℃ ·(10 a)-1], 贾玲等[62]基于疏勒河流域1951~2018年月最高气温数据得出研究区季节最高气温均呈上升趋势[0.01~0.36℃ ·(10 a)-1], 这与本研究在具体数值上存在略微差异; 李虹雨等[63]基于观测站点数据分析1951~2014年内蒙古地区季节最高气温变化特征, 结果显示研究区整体呈升温趋势, 但其变化速率的空间分布在局部地区与本研究略有不同.

3.3 不确定性分析

由于高温变化是一个长期而复杂的过程, 探究高温天气在较长时间序列下的变化特征具有重要意义, 而本研究的分析时段为近70年, 相对于地球宏观周期性气候变化而言, 时长较短, 研究结果存在不确定性.此外, 影响高温天气变化的因素众多, 例如, 平流层水汽增加、太阳活动增强和火山活动等亦可使得区域辐射持续上升, 进而影响区域气温变化.Solomon等[64]研究发现平流层水汽增加导致全球地表变暖加速; Huang等[65]和郭纪君等[66]认为太阳活动对全球变暖起着决定性作用, 当太阳黑子活动增强时, 地球表面气温呈升温态势; 温飞等[67]研究指出火山活动和太阳辐射量变化等均能显著影响地表气候变化.但由于本研究缺少相关实测数据, 无法对以上因素作进一步探讨.此外, 本研究仅从统计分析来探讨我国季节AMT与大气环流因子间的关系, 而各大气环流因子位相间的变化机制未作深入分析, 具有一定的局限性.鉴于以上, 今后研究需着重以下工作的开展:重建和预测历史高温天气变化, 模拟和分析大气环流因子间相互作用机制, 以此对高温天气的变化趋势以及影响因素进行全面的分析和研究, 以实现更精确的判断.

4 结论

(1) 1950~2019年我国各季节AMT均呈显著上升趋势, 春季升温速率最快, 秋季、夏季和冬季次之.各季节AMT产生增温突变的时间各不相同, 主要集中于20世纪90年代至21世纪初.

(2) 我国各季节AMT年际平均倾向率存在明显的空间异质性.除夏季以外, 其它季节AMT的年际平均倾向率均由南向北逐渐增加, 西北和东北地区增温幅度最为明显.

(3) 各季节AMT与大气环流因子存在复杂联系, 且在不同时频域作用强度各异.其中, 春季AMT与AO呈正相关, 夏季AMT与PDO呈负相关, 秋季AMT主要受各大气环流因子较长年际周期影响, 冬季AMT与AO、NAO存在显著共振周期.

参考文献
[1] 孙颖. 全球变暖趋势日益明晰[J]. 中国科学基金, 2021, 35(2): 226-228. DOI:10.16262/j.cnki.1000-8217.2021.02.011
[2] IPCC. Climate Change 2013: The Physical Science Basis. Contribution of Working Group I to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change[M]. Cambridge & New York: Cambridge University Press, 2014.
[3] Press F. Earth science and society[J]. Nature, 2008, 451(7176): 301-303. DOI:10.1038/nature06595
[4] 彭京备, 张庆云, 布和朝鲁. 2006年川渝地区高温干旱特征及其成因分析[J]. 气候与环境研究, 2007, 12(3): 464-474.
Peng J B, Zhang Q Y, Cholaw B. On the characteristics and possible causes of a severe drought and heat wave in the Sichuan-Chongqing region in 2006[J]. Climatic and Environmental Research, 2007, 12(3): 464-474. DOI:10.3969/j.issn.1006-9585.2007.03.026
[5] 孙劭, 李多, 王遵娅, 等. 2018年全球重大天气气候事件及其成因[J]. 气象, 2019, 45(4): 533-542.
Sun S, Li D, Wang Z Y, et al. Global major weather and climate events in 2018 and the possible causes[J]. Meteorological Monthly, 2019, 45(4): 533-542.
[6] 贾斌英. 2010—2019年中国东北地区森林火灾时空特征研究[J]. 防护林科技, 2021(4): 18-21.
Jia B Y. Temporal and spatial characteristics of forest fires in northeast China from 2010 to 2019[J]. Protection Forest Science and Technology, 2021(4): 18-21. DOI:10.13601/j.issn.1005-5215.2021.04.007
[7] Schweinsberg S, Darcy S, Beirman D. 'Climate crisis' and 'bushfire disaster': Implications for tourism from the involvement of social media in the 2019-2020 Australian bushfires[J]. Journal of Hospitality and Tourism Management, 2020, 43: 294-297. DOI:10.1016/j.jhtm.2020.03.006
[8] Song J, Pan R B, Yi W Z, et al. Ambient high temperature exposure and global disease burden during 1990-2019: An analysis of the global burden of disease study 2019[J]. Science of the Total Environment, 2021, 787. DOI:10.1016/j.scitotenv.2021.147540
[9] Schär C, Jendritzky G. Hot news from summer 2003[J]. Nature, 2004, 432(7017): 559-560. DOI:10.1038/432559a
[10] Li X, Ren G Y, Wang S Y, et al. Change in the heatwave statistical characteristics over China during the climate warming slowdown[J]. Atmospheric Research, 2021, 247. DOI:10.1016/j.atmosres.2020.105152
[11] 冯相昭, 邹骥, 马珊, 等. 极端气候事件对中国农村经济影响的评价[J]. 农业技术经济, 2007, 160(2): 19-25.
Feng X Z, Zou J, Ma S, et al. The impact of extreme climate on rural economy in China: An evaluation[J]. Journal of Agrotechnical Economics, 2007, 160(2): 19-25. DOI:10.3969/j.issn.1000-6370.2007.02.003
[12] 任国玉, 初子莹, 周雅清, 等. 中国气温变化研究最新进展[J]. 气候与环境研究, 2005, 10(4): 701-716.
Ren G Y, Chu Z Y, Zhou Y Q, et al. Recent progresses in studies of regional temperature changes in China[J]. Climatic and Environmental Research, 2005, 10(4): 701-716.
[13] 陈建宇, 邵天杰, 赵景波. 陕西榆林地区1954—2013年极端气温变化特征[J]. 中国沙漠, 2016, 36(5): 1426-1434.
Chen J Y, Shao T J, Zhao J B. Extreme temperature changes from 1954 to 2013 in Yulin, Shaanxi, China[J]. Journal of Desert Research, 2016, 36(5): 1426-1434.
[14] Gocic M, Trajkovic S. Analysis of changes in meteorological variables using Mann-Kendall and Sen's slope estimator statistical tests in Serbia[J]. Global and Planetary Change, 2013, 100: 172-182. DOI:10.1016/j.gloplacha.2012.10.014
[15] Zuo J Q, Ren H L, Li W J, et al. Interdecadal variations in the relationship between the winter North Atlantic oscillation and temperature in south-Central China[J]. Journal of Climate, 2016, 29(20): 7477-7493. DOI:10.1175/JCLI-D-15-0873.1
[16] 翟盘茂, 余荣, 郭艳君, 等. 2015/2016年强厄尔尼诺过程及其对全球和中国气候的主要影响[J]. 气象学报, 2016, 74(3): 309-321.
Zhai P M, Yu R, Guo Y J, et al. The strong EI Niño in 2015/2016 and its dominant impacts on global and China's climate[J]. Acta Meteorologica Sinica, 2016, 74(3): 309-321.
[17] 田立鑫, 韩美, 徐泽华, 等. 近50年淮河流域气温时空变化及其与PDO的关系[J]. 水土保持研究, 2019, 26(6): 240-248.
Tian L X, Han M, Xu Z H, et al. Temporal and spatial change of temperature in the Huaihe River Basin during the past 50 years and its relationship with PDO[J]. Research of Soil and Water Conservation, 2019, 26(6): 240-248.
[18] 马双梅, 祝从文, 刘伯奇. 2019年4~6月云南持续性高温天气的大气环流异常成因[J]. 大气科学, 2021, 45(1): 165-180.
Ma S M, Zhu C W, Liu B Q. Possible causes of persistently Extreme-Hot-Days-Related circulation anomalies in Yunnan from April to June 2019[J]. Chinese Journal of Atmospheric Sciences, 2021, 45(1): 165-180.
[19] Peng S Z, Ding Y X, Liu W Z, et al. 1 km monthly temperature and precipitation dataset for China from 1901 to 2017[J]. Earth System Science Data, 2019, 11(4): 1931-1946. DOI:10.5194/essd-11-1931-2019
[20] 史雯雨, 杨胜勇, 李增永, 等. 近57年金沙江流域气温变化特征及未来趋势预估[J]. 水土保持研究, 2021, 28(1): 211-217.
Shi W Y, Yang S Y, Li Z Y, et al. Variation characteristics and the future trend estimation of temperature in Chinsha River Basin over the past 57 years[J]. Research of Soil and Water Conservation, 2021, 28(1): 211-217.
[21] 刘丹, 赵景波. 黑龙江哈尔滨地区近60年极端气温变化研究[J]. 地球与环境, 2017, 45(6): 587-599.
Liu D, Zhao J B. Extreme temperature changes in Harbin region, Heilongjiang Province, China over the past 60 years[J]. Earth and Environment, 2017, 45(6): 587-599.
[22] Huang Y F, Lu C Y, Lei Y F, et al. Spatio-temporal variations of temperature and precipitation during 1951-2019 in arid and semiarid region, China[J]. Chinese Geographical Science, 2022, 32(2): 285-301. DOI:10.1007/s11769-022-1262-8
[23] 江铭诺, 李冠毅, 何婉文, 等. 1951—2020年广州地区极端气温指数年月尺度变化研究[J]. 热带气象学报, 2022, 38(2): 239-252.
Jiang M N, Li G Y, He W W, et al. A study on the annual and monthly variations of the extreme temperature index in Guangzhou from 1951 to 2020[J]. Journal of Tropical Meteorology, 2022, 38(2): 239-252.
[24] 李卓群, 刘星才. 1961—2019年辽宁省高温天气变化特征[J]. 应用生态学报, 2021, 32(11): 4059-4067.
Li Z Q, Liu Q C. Variations of high temperature from 1961 to 2019 in Liaoning Province, China[J]. Chinese Journal of Applied Ecology, 2021, 32(11): 4059-4067.
[25] 李春兰, 包玉海, 包刚, 等. 近57年锡林浩特市平均最高、最低气温变化特征分析[J]. 水土保持研究, 2013, 20(5): 257-260, 265.
Li C L, Bao Y H, Bao G, et al. Changes characteristic of mean maximum and minimum temperature in Xilinhot City in recent 57 years[J]. Research of Soil and Water Conservation, 2013, 20(5): 257-260, 265.
[26] 者萌, 张雪芹, 沈鹏珂, 等. 京津冀地区1957-2017年气温变化时空格局[J]. 水土保持研究, 2020, 27(5): 194-201.
Zhe M, Zhang X Q, Shen P K, et al. Spatial-temporal pattern of temperature variation in Beijing-Tianjin-Hebei region over the period 1957-2017[J]. Research of Soil and Water Conservation, 2020, 27(5): 194-201.
[27] 徐乔婷, 陈涟, 范月华, 等. 基于SPEI指数的兰州干旱特征与气候指数的关系[J]. 水文, 2021, 41(2): 56-62.
Xu Q T, Chen L, Fan Y H, et al. Relationship between Lanzhou drought and climate indices based on SPEI index[J]. Journal of China Hydrology, 2021, 41(2): 56-62.
[28] Li Y L, Wen Y L, Lai H X, et al. Drought response analysis based on cross wavelet transform and mutual entropy[J]. Alexandria Engineering Journal, 2020, 59(3): 1223-1231.
[29] 高弋斌, 路春燕, 钟连秀, 等. 1951—2016年中国沿海地区气温与降水量的时空特征[J]. 森林与环境学报, 2019, 39(5): 530-539.
Gao Y B, Lu C Y, Zhong L X, et al. Temporal and spatial characteristics of temperature and precipitation in China's coastal areas from 1951 to 2016[J]. Journal of Forest and Environment, 2019, 39(5): 530-539.
[30] 郑景云, 卞娟娟, 葛全胜, 等. 1981~2010年中国气候区划[J]. 科学通报, 2013, 58(30): 3088-3099.
Zheng J Y, Bian J J, Ge Q S, et al. The climate regionalization in China for 1981-2010[J]. Chinese Science Bulletin, 2013, 58(30): 3088-3099.
[31] 庞家泰, 段金亮, 张瑞, 等. 2000-2019年渭河流域植被覆盖度时空演变特征及气候响应[J]. 水土保持研究, 2021, 28(5): 230-237.
Pang J T, Duan J L, Zhang R, et al. Characteristics of spatiotemporal evolution and climate response of vegetation cover in the wei river basin from 2000 to 2019[J]. Research of Soil and Water Conservation, 2021, 28(5): 230-237.
[32] 李明龙, 贾梦丹, 孙天成, 等. 三峡库区非点源污染氮磷负荷时空变化及其来源解析[J]. 环境科学, 2021, 42(4): 1839-1846.
Li M L, Jia M D, Sun T C, et al. Spatiotemporal change and source apportionment of non-point source nitrogen and phosphorus pollution loads in the Three Gorges reservoir area[J]. Environmental Science, 2021, 42(4): 1839-1846.
[33] 王景才, 郭佳香, 徐蛟, 等. 近55年淮河上中游流域气候要素多时间尺度演变特征及关联性分析[J]. 地理科学, 2017, 37(4): 611-619.
Wang J C, Guo J X, Xu J, et al. Multi-time scales change characteristics and relationship of meteorological variables in the upper and middle regions of the Huaihe River Basin in recent 55 years[J]. Scientia Geographica Sinica, 2017, 37(4): 611-619.
[34] 段家乐, 巨天珍, 黄蕊蕊, 等. 粤港澳大湾区吸收性气溶胶的解析[J]. 环境科学, 2021, 42(8): 3652-3662.
Duan J L, Ju T Z, Huang R R, et al. Analysis of ultraviolet aerosol index in Guangdong-Hong Kong-Macao Greater Bay Area[J]. Environmental Science, 2021, 42(8): 3652-3662.
[35] Wallace J M, Gutzler D S. Teleconnections in the geopotential height field during the northern hemisphere winter[J]. Monthly Weather Review, 1981, 109(4): 784-812.
[36] 谭本馗, 陈文. 中高纬度大气遥相关动力学及其对东亚冬季气候影响的研究进展[J]. 气象学报, 2014, 72(5): 908-925.
Tan B K, Chen W. Progress in the study of the dynamics of extratropical atmospheric teleconnection patterns and their impacts on East Asian climate[J]. Acta Meteorologica Sinica, 2014, 72(5): 908-925.
[37] Wu B Y, Wang J. Possible impacts of winter arctic oscillation on Siberian high, the East Asian winter monsoon and sea-ice extent[J]. Advances in Atmospheric Sciences, 2002, 19(2): 297-320.
[38] 龚道溢, 王绍武. 近百年北极涛动对中国冬季气候的影响[J]. 地理学报, 2003, 58(4): 559-568.
Gong D Y, Wang S W. Influence of arctic oscillation on winter climate over China[J]. Acta Geographica Sinica, 2003, 58(4): 559-568.
[39] 李忠贤, 于怡秋, 邓伟涛, 等. 春季北大西洋三极型海温异常变化及其与NAO和ENSO的联系[J]. 气象科学, 2019, 39(6): 721-730.
Li Z X, Yu Y Q, Deng W T, et al. Variation characteristics of North Atlantic Tri-Polar Sea surface temperature in spring and its relationship with NAO and ENSO[J]. Journal of the Meteorological Sciences, 2019, 39(6): 721-730.
[40] 范可, 刘辉. 近百年东亚冬季气温及其大气环流变化型态[J]. 大气科学, 2013, 37(2): 383-394.
Fan K, Liu H. Winter temperature over east Asia and atmospheric circulation patterns in the last 100 years[J]. Chinese Journal of Atmospheric Sciences, 2013, 37(2): 383-394.
[41] 龚道溢, 王绍武. 北大西洋涛动指数的比较及其年代际变率[J]. 大气科学, 2000, 24(2): 187-192.
Gong D Y, Wang S W. The North Atlantic oscillation index and its interdecadal variability[J]. Chinese Journal of Atmospheric Sciences, 2000, 24(2): 187-192.
[42] 叶许春, 许崇育, 张丹, 等. 长江中下游夏季降水变化与亚洲夏季风系统的关系[J]. 地理科学, 2018, 38(7): 1174-1182.
Ye X C, Xu C Y, Zhang D, et al. Variation of summer precipitation and its connection with Asian monsoon system in the middle-lower Yangtze River basin[J]. Scientia Geographica Sinica, 2018, 38(7): 1174-1182.
[43] 吴国丽, 李思航, 刘子洲, 等. 2015/2016年超强厄尔尼诺期间热带西太平洋海洋环流变化特征[J]. 海洋科学, 2020, 44(1): 8-18.
Wu G L, Li S H, Liu Z Z, et al. Variations in the tropical western pacific gyre during the 2015/2016 super El Niño event[J]. Marine Sciences, 2020, 44(1): 8-18.
[44] 严韬, 徐明洁, 葛非凡, 等. 厄尔尼诺/拉尼娜事件对区域气温的影响与预测——以沈阳地区为例[J]. 安徽农业大学学报, 2019, 46(1): 57-64.
Yan T, Xu M J, Ge F F, et al. The effect and forecaster of the El Niño/La Niña events on climate in Shenyang[J]. Journal of Anhui Agricultural University, 2019, 46(1): 57-64.
[45] Fang C F, Wu L X, Zhang X. The impact of global warming on the pacific decadal oscillation and the possible mechanism[J]. Advances in Atmospheric Sciences, 2014, 31(1): 118-130.
[46] 沈皓俊, 游庆龙, 王朋岭, 等. 1961—2014年中国高温热浪变化特征分析[J]. 气象科学, 2018, 38(1): 28-36.
Shen H J, You Q L, Wang P L, et al. Analysis on heat waves variation features in China during 1961—2014[J]. Journal of the Meteorological Sciences, 2018, 38(1): 28-36.
[47] 梁苏洁, 丁一汇, 赵南, 等. 近50年中国大陆冬季气温和区域环流的年代际变化研究[J]. 大气科学, 2014, 38(5): 974-992.
Liang S J, Ding Y H, Zhao N, et al. Analysis of the interdecadal changes of the wintertime surface air temperature over mainland China's mainland and regional atmospheric circulation characteristics during 1960-2013[J]. Chinese Journal of Atmospheric Sciences, 2014, 38(5): 974-992.
[48] 徐忆菲, 彭丽霞, 李季, 等. 1951—2013年我国冬季气温年代际变化与PDO的关系[J]. 气象科技, 2017, 45(4): 669-677.
Xu Y F, Peng L X, Li J, et al. Relationship between PDO and interdecadal change of winter temperature in China from 1951 to 2013[J]. Meteorological Science and Technology, 2017, 45(4): 669-677.
[49] 李宇, 周德成, 闫章美, 等. 中国主要城市的城市化对局地增温的贡献[J]. 环境科学, 2022, 43(5): 2822-2830.
Li Y, Zhou D C, Yan Z M, et al. Contribution of urbanization to local warming in major cities of China[J]. Environmental Science, 2022, 43(5): 2822-2830.
[50] Sun Y X, Liu S L, Dong Y H, et al. Spatio-temporal evolution scenarios and the coupling analysis of ecosystem services with land use change in China[J]. Science of the Total Environment, 2019, 681: 211-225.
[51] Tian B, Wu W T, Yang Z Q, et al. Drivers, trends, and potential impacts of long-term coastal reclamation in China from 1985 to 2010[J]. Estuarine, Coastal and Shelf Science, 2016, 170: 83-90.
[52] 于竹筱, 张红旗, 许尔琪. 基于"三生"视角的中国林地分类提取[J]. 自然资源学报, 2021, 36(5): 1136-1148.
Yu Z X, Zhang H Q, Xu E Q. Classification and extraction of forest land in China based on the perspective of "production-living-ecology"[J]. Journal of Natural Resources, 2021, 36(5): 1136-1148.
[53] 曹孜, 陈洪波. 城市化和能源消费的门槛效应分析与预测[J]. 中国人口·资源与环境, 2015, 25(11): 59-68.
Cao Z, Chen H B. Threshold effects analysis and prediction for China's urbanization and energy consumption[J]. China Population, Resources and Environment, Resources and Environment, 2015, 25(11): 59-68.
[54] 陈朝军, 袁道先, 程海, 等. 人类活动和气候变化触发了中国西南石漠化的扩张[J]. 中国科学: 地球科学, 2021, 64(11): 1950-1963.
Chen C J, Yuan D, Cheng H, et al. Human activity and climate change triggered the expansion of rocky desertification in the karst areas of Southwestern China[J]. Science China Earth Sciences, 2021, 64(10): 1761-1773.
[55] Skytt T, Nielsen S N, Jonsson B G. Global warming potential and absolute global temperature change potential from carbon dioxide and methane fluxes as indicators of regional sustainability-A case study of Jämtland, Sweden[J]. Ecological Indicators, 2020, 110. DOI:10.1016/j.ecolind.2019.105831
[56] Luo Z B, Lam S K, Fu H, et al. Temporal and spatial evolution of nitrous oxide emissions in China: assessment, strategy and recommendation[J]. Journal of Cleaner Production, 2019, 223: 360-367.
[57] Cai D W, Ge Q S, Wang X M, et al. Contributions of ecological programs to vegetation restoration in arid and semiarid China[J]. Environmental Research Letters, 2020, 15(11). DOI:10.1088/1748-9326/abbde9
[58] Zheng Y, Dong L, Xia Q, et al. Effects of revegetation on climate in the Mu Us Sandy Land of China[J]. Science of the Total Environment, 2020, 739. DOI:10.1016/j.scitotenv.2020.139958
[59] Yu L X, Liu Y, Liu T X, et al. Impact of recent vegetation greening on temperature and precipitation over China[J]. Agricultural and Forest Meteorology, 2020, 295. DOI:10.1016/j.agrformet.2020.108197
[60] 周莹, 马龙, 梁珑腾. 内蒙古气温突变及其影响因子间的响应关系[J]. 水土保持研究, 2021, 28(3): 184-193.
Zhou Y, Ma L, Liang L T. Response of temperature mutation in inner Mongolia to its influencing factors[J]. Research of Soil and Water Conservation, 2021, 28(3): 184-193.
[61] 梁珑腾, 马龙, 刘廷玺, 等. 1951—2014年中国北方地区季节气温突变与变暖停滞年份的时空变异性[J]. 自然资源学报, 2018, 33(12): 2149-2166.
Liang L T, Ma L, Liu Y X, et al. The seasonal spatiotemporal variation of the temperature mutation and warming hiatus over northern China during 1951-2014[J]. Journal of Natural Resources, 2018, 33(12): 2149-2166.
[62] 贾玲, 孙栋元, 牛最荣, 等. 疏勒河流域最高、最低气温变化规律[J]. 水土保持研究, 2022, 29(1): 281-287.
Jia L, Sun D Y, Niu Z R, et al. Pattern of maximum and minimum temperature variation in Shule River basin[J]. Research of Soil and Water Conservation, 2022, 29(1): 281-287.
[63] 李虹雨, 马龙, 刘廷玺, 等. 1951-2014年内蒙古地区气温、降水变化及其关系[J]. 冰川冻土, 2017, 39(5): 1098-1112.
Li H Y, Ma L, Liu Y X, et al. Change and relationship of temperature and precipitation in Inner Mongolia during 1951-2014[J]. Journal of Glaciology and Geocryology, 2017, 39(5): 1098-1112.
[64] Solomon S, Rosenlof K H, Portmann R W, et al. Contributions of stratospheric water vapor to decadal changes in the rate of global warming[J]. Science, 2010, 327(5970): 1219-1223.
[65] Huang C, Rao Z G, Li Y X, et al. Holocene summer temperature in arid central Asia linked to millennial-scale north Atlantic climate events and driven by centennial-scale solar activity[J]. Palaeogeography, Palaeoclimatology, Palaeoecology, 2020, 556. DOI:10.1016/j.palaeo.2020.109880
[66] 郭纪君, 郭治龙, 王宁练. 1880年以来地球陆地-海洋温度指数中的太阳黑子活动周期研究[J]. 高原气象, 2020, 39(4): 851-858.
Guo J J, Guo Z L, Wang N L. Study on the cycles of sunspots in the Earth's land-ocean temperature index since 1880[J]. Plateau Meteorology, 2020, 39(4): 851-858.
[67] 温飞, 游爱华, 薛积彬. 过去2000年中国东南地区气温与降水的集成重建与分析[J]. 热带地理, 2021, 41(4): 778-789.
Wen F, You A H, Xue J B. Integrated reconstruction and analysis of temperature and precipitation changes in southeastern China during the last 2000 years[J]. Tropical Geography, 2021, 41(4): 778-789.