冰川冻土, 2022, 44(1): 34-45 doi: 10.7522/j.issn.1000-0240.2022.0018

第二次青藏高原综合科学考察研究

1979—2100年青藏高原夏季大气0 ℃层高度变化分析

张怡,, 段克勤,, 石培宏

陕西师范大学 地理科学与旅游学院,陕西 西安 710119

Analysis of the 0 ℃ level height variation over the Qinghai-Tibet Plateau in summer from 1979 to 2100

ZHANG Yi,, DUAN Keqin,, SHI Peihong

School of Geography and Tourism,Shaanxi Normal University,Xi’an 710119,China

通讯作者: 段克勤,教授,主要从事全球变化研究. E-mail: kqduan@snnu.edu.cn

收稿日期: 2021-03-30   修回日期: 2021-12-10  

基金资助: 第二次青藏高原综合科学考察研究项目.  2019QZKK0201
国家自然科学基金项目.  41571062
中国博士后基金项目.  2017M610622

Received: 2021-03-30   Revised: 2021-12-10  

作者简介 About authors

张怡,硕士研究生,主要从事气候变化与自然灾害研究.E-mail:zhangyi19@snnu.edu.cn , E-mail:zhangyi19@snnu.edu.cn

摘要

大气0 ℃层高度是决定青藏高原冰冻圈消融状态的重要指标。基于ERA5再分析资料,分析了1979—2019年青藏高原夏季大气0 ℃层高度时空变化,发现青藏高原夏季大气0 ℃层高度介于4 423~5 972 m之间,以高原中南部(30°~32° N,83.5°~88.5° E)为高值中心,呈纬向分带状向四周逐渐降低。过去41 a青藏高原夏季大气0 ℃层高度总体呈持续上升趋势,高原北部上升趋势大于南部,祁连山地区上升趋势最为明显,为60 m∙(10a)-1,而在高原西南部略呈下降趋势。平均而言,青藏高原夏季地面温度每升高1 ℃,大气0 ℃层高度升高122 m。利用CMIP6模式数据,预估在SSP1-2.6、SSP2-4.5、SSP3-7.0和SSP5-8.5四种社会共享路径情景下,2020—2100年期间青藏高原夏季大气0 ℃层高度都呈现升高趋势,但不同情景下升高趋势在空间上差别较大。相对于1979—2014年参考时段,在四种情景下,到2081—2100年青藏高原夏季平均大气0 ℃层高度将分别升高265 m、394 m、576 m 和729 m;相对应的是到2081—2100年,在高原上处于夏季大气0 ℃层高度以下的冰川面积分别为第二次冰川编目数据的79%、86%、94%和98%。仅从夏季大气0 ℃层高度变化角度看,在SSP5-8.5情景下,到本世纪末期,预估除帕米尔高原和昆仑山西北部地区外,青藏高原其他地区的冰川在夏季将不存在积累区。

关键词: 青藏高原 ; 大气0 ℃层高度 ; ERA5 ; CMIP6

Abstract

As an Asian water tower, the Qinghai-Tibet Plateau is extremely sensitive to global warming due to its unique cryosphere landform. At present, due to the warming of the plateau, a series of ecological and environmental problems have been caused, such as glacier retreat, frozen soil melting and so on. In recent years, the near-surface air temperature of the Qinghai-Tibet Plateau has been increasing at the rate of 0.3~0.4 ℃ (10a)-1, which is twice the global average. The near-surface air temperature increase will affect the upper atmospheric temperature in the form of long wave radiation. An important indicator of the change of upper air temperature is the 0 ℃ level height. As the lowest height of free air temperature of 0 ℃, the 0 ℃ level height is also an important indicator of global change. The height of atmospheric 0 ℃ layer determines the height of water vapor phase transition in the air. If the 0 ℃ level height is higher than the surface, the cryosphere on the surface is in a state of melting, otherwise it is in a state of freezing. Therefore, the 0 ℃ level height directly determines the state of melting of the cryosphere on the surface of the plateau. Therefore, quantifying the change of the 0 ℃ level height over the Qinghai-Tibet Plateau is very important for understanding the climate change of the plateau and the change state of the cryosphere. Based on ERA5 reanalysis data, the temporal and spatial variation of the 0 ℃ level height in summer over the Qinghai-Tibet Plateau from 1979 to 2019 is analyzed. It is found that the 0 ℃ level height in summer over the Qinghai-Tibet Plateau ranges from 4 423 to 5 972 m, with the central and southern part of the Plateau (30~32° N, 83.5~88.5° E) as the high value center, and gradually decreases in zonal direction. In the past 41 years, the 0 ℃ level height in summer over the Qinghai-Tibet Plateau has been increasing continuously. The increasing trend in the north of the plateau is greater than that in the south. The increasing trend in the Qilian Mountains is the most obvious, which is 60 m∙(10a)-1, while it is slightly decreasing in the southwest of the plateau. On average, the 0 ℃ level height increases by 122 m for every 1 ℃ rise in summer near-surface air temperature over the Qinghai-Tibet Plateau. Using CMIP6 model data, it is estimated that the 0 ℃ level height in summer over the Qinghai-Tibet Plateau will increase during 2020—2100 under SSP1-2.6, SSP2-4.5, SSP3-7.0 and SSP5-8.5 social shared path scenarios, but the increasing trend under different scenarios is quite different in space. Compared with the reference period of 1979—2014, under the four scenarios, the summer mean the 0 ℃ level height of the Qinghai-Tibet Plateau will increase by 265 m, 394 m, 576 m and 729 m respectively in 2081—2100; Correspondingly, from 2081 to 2100, the area of glaciers below 0 ℃ in summer on the plateau is 79%, 86%, 94% and 98% of the second glacial catalogue data, respectively. Only from the perspective of the 0 ℃ level height change in summer, it is estimated that the accumulation area of glaciers in the Qinghai-Tibet Plateau will disappear in summer under SSP5-8.5 scenario until the end of this century, except for the Pamir Plateau and the northern Kunlun Mountains.

Keywords: Qinghai-Tibet Plateau ; the 0 ℃ level height ; ERA5 ; CMIP6

PDF (9131KB) 元数据 多维度评价 相关文章 导出 EndNote| Ris| Bibtex  收藏本文

本文引用格式

张怡, 段克勤, 石培宏. 1979—2100年青藏高原夏季大气0 ℃层高度变化分析[J]. 冰川冻土, 2022, 44(1): 34-45 doi:10.7522/j.issn.1000-0240.2022.0018

ZHANG Yi, DUAN Keqin, SHI Peihong. Analysis of the 0 ℃ level height variation over the Qinghai-Tibet Plateau in summer from 1979 to 2100[J]. Journal of Glaciology and Geocryology, 2022, 44(1): 34-45 doi:10.7522/j.issn.1000-0240.2022.0018

0 引言

青藏高原作为亚洲水塔,其独特的冰冻圈地貌对全球变暖极其敏感。目前因高原变暖,已造成冰川退缩、冻土消融等一系列的生态环境问题1-5。近年来青藏高原地面温度以0.3~0.4 ℃∙(10a)-1的升温速率上升,是全球平均增温速率的2倍6-7,地面增温会以长波辐射形式影响其上的高空大气温度。表征高空大气温度变化的一个重要指标就是大气0 ℃层高度,作为自由大气温度为0 ℃的最低高度,大气0 ℃层高度决定了空中水汽相变发生的高度。若大气0 ℃层高度高于地表,则地表的冰冻圈处于消融状态,故大气0 ℃层高度可反映高原地表冰冻圈的消融状态。因此,量化青藏高原地区大气0 ℃层高度变化,对认识高原的气候变化,以及冰冻圈的变化状态十分关键。

研究指出在热带地区大气0 ℃层高度变化与海表温度密切相关,且是影响热带地区冰川变化的重要因素8-10。近年来在我国西北、西南、黄河流域和高亚洲等地区,也进行了大气0 ℃层高度变化的研究,发现大气0 ℃层高度呈升高趋势,并与地面温度的变化呈显著正相关关系11-20。在青藏高原地区,夏季大气0 ℃层高度的变化与冰川消融量21、零平衡线高度22以及冰川物质平衡23之间存在很好的线性关系,并直接影响冰川、积雪与冻土区的消融面积和强度,进而对区域径流产生重要影响12,可诱发洪水24-26及断流1227-32等灾害。通过建立夏季径流量对大气0 ℃层高度变化的响应关系,也可预测评估夏季河流径流量的变化33-37。因此,夏季大气0 ℃层高度的变化不仅是大气热状况的反映,也是影响冰冻圈各要素变化的重要指标38

采用探空站数据,以及ERA-Interim1739和NCEP/NCAR92330-3139-40等再分析资料,对青藏高原不同区域的夏季大气0 ℃层高度变化进行了研究,但对青藏高原全局的研究相对薄弱。探空数据获取的大气0 ℃层高度可信度高,但由于青藏高原探空站稀少,无法获得整个高原面上大气0 ℃层高度的变化。再分析资料时空分布连续,由其可获得较高精度的大气0 ℃层高度变化3941-45。ERA5作为最新一代的再分析资料,与上一代的ERA-Interim再分析数据相比,其空间分辨率更高,数据质量更好46-48。青藏高原夏季大气0 ℃层高度未来如何变化,目前尚无研究,也缺乏对夏季大气0 ℃层高度变化对青藏高原冰川变化影响方面的研究。基于此,本文的目的是:(1)综合探空站数据、ERA5再分析资料以及CMIP6模式数据,分析1979—2100年高原夏季大气0 ℃层高度变化的特征及规律。(2)从高原夏季大气0 ℃层高度变化角度,预估青藏高原冰川未来变化的可能态势。

1 数据与方法

青藏高原平均海拔在4 000 m以上,为体现高原的高海拔特性,选取25°~40° N,70°~105° E范围内,海拔2 500 m以上的区域作为研究区(图1)。

图1

图1   研究区范围和探空站分布

Fig. 1   Location of the study area and adiosonde station


1.1 数据来源

本研究采用欧洲中期天气预报中心发布的第五代全球气候大气再分析产品ERA5(https://www.ecmwf.int),空间分辨率为0.25°×0.25°、垂直方向分为37个层,选择1 000 hPa、600 hPa、500 hPa和400 hPa四个高度标准气压层的逐月平均位势高度和温度数据。为了探究地面温度与大气0 ℃层高度的相关性,选取逐月平均距离地面2 m处的空气温度(简称地面温度)数据,时间跨度为1979—2019年。

为检验ERA5再分析资料在青藏高原地区的适用性,以区域内有探空资料的8个台站(图1)为参考站(https://data.cma.cn/),选取地表、700 hPa、500 hPa和400 hPa逐月平均位势高度和温度计算大气0 ℃层高度。为了确保数据质量,使用最大惩罚F法49-50对原始探空数据进行均一性检验和订正2338,以消除或减少观测的系统误差。

为探讨未来大气0 ℃层高度的变化,采用国际耦合模式比较计划第六阶段(CMIP6)的MIROC6模式模拟结果(https://esgf-node.llnl.gov/projects/cmip6/),其空间分辨率为1.4°×1.4°,模拟时段为1979—2100年。选取SSP1-2.6、SSP2-4.5、SSP3-7.0和SSP5-8.5四种共享社会经济路径情景下51-52,三个标准气压层(600 hPa、500 hPa和400 hPa)的逐月平均位势高度和温度数据。为与探空站和ERA5再分析资料比较,采用双线性插值方法将MIROC6模式数据插值到0.25°×0.25°的经纬网格上。

1.2 研究方法

大气0 ℃层高度可通过已知气压层的高度和空气温度进行线性插值得到。逐一判断各月、各时次0 ℃层所在位置上下的两个标准气压层,设气压层之间温度在垂直方向上呈线性变化,由线性插值计算得到各月大气0 ℃层高度。青藏高原冬半年近地表温度往往低于0 ℃,理论上的大气0 ℃层高度在地表以下,故本文只研究夏季(6—8月)大气0 ℃层高度的变化。

计算大气0 ℃层高度的公式为:

H=Tdown×Hdown-HupTup-Tdown+Hdown

式中:H为大气0 ℃层的高度(m);HdownHup分别为大气0 ℃层所在位置的上下两个标准气压层高度(m);TdownTup分别为HdownHup对应的温度(℃)。

气候倾向率反映了各气候因子长期趋势变化的方向和程度。采用最小二乘法计算气候要素与时间的线性回归系数a,使气候要素的变化可用一元线性回归方程表示53,其计算公式为:

Xt=at+bt=1979,1980,…,2100

式中:Xt 为气候要素拟合值,ab为回归系数。

a>0时,说明随t的增加X呈上升趋势;当a<0时,说明随t的增加X呈下降趋势。a值的大小反映上升或下降的倾向程度。将10称为气候倾向率。大气0 ℃层高度和地面温度倾向率的单位分别是m∙(10a)-1和℃∙(10a)-1。本文运用t检验法对气候倾向率进行显著性检验,选取信度为0.05和0.01的显著性水平,分别对应95%和99%的置信度。

2 结果分析

2.1 ERA5再分析资料和MIROC6模式数据在青藏高原地区的适用性评价

采用双线性插值方法,将探空站邻近的4个格点数据插值到站点位置。图2对比了由ERA5和MIROC6数据计算得到的夏季大气0 ℃层高度与探空站观测值随时间变化。表1进一步从平均值、相关系数、绝对偏差和相对偏差评估了计算值和观测值之间的差异。可见由ERA5获得的大气0 ℃层高度,无论从平均值、变化趋势以及年际波动上,与由探空站观测值得到的大气0 ℃层高度变化极为一致,均通过了信度为0.01的显著性检验,两者之间的绝对偏差都小于50 m(在那曲站偏大,为97 m)。因此,由ERA5再分析资料计算青藏高原夏季大气0 ℃层高度是可信的,完全可以弥补在高原观测数据的不足,可作为研究青藏高原夏季大气0 ℃层高度变化的代用资料。

图2

图2   ERA5再分析资料和MIROC6模式数据与探空站数据对比

Fig. 2   Comparison of ERA5 reanalysis data and MIROC6 model data with radiosonde data


表1   ERA5再分析资料和MIROC6模式数据与探空站数据夏季大气0 ℃层高度的误差指标 (m)

Table1  1 Error index of the summer 0 ℃ level height between ERA5 reanalysis data, MIROC6 model data and radiosonde data

合作格尔木都兰西宁玉树甘孜那曲昌都
平均值探空站5 0375 1265 1364 8855 4975 5175 8335 638
ERA55 0235 1025 1244 9305 4605 5005 7375 593

MIROC6

校正前

4 9704 9464 9964 8755 3205 3895 6275 428

MIROC6

校正后

5 0255 1125 1284 9155 4625 5025 7455 596
相关系数ERA50.960.950.970.940.940.970.930.90
CMIP60.190.230.280.260.250.150.020.09
绝对偏差ERA5-14-24-1245-37-17-96-45

MIROC6

校正前

-67-180-140-10-177-128-206-210

MIROC6

校正后

-12-14-830-35-15-88-42
相对偏差ERA5-0.28%-0.47%-0.23%0.92%-0.67%-0.31%-1.65%-0.80%

MIROC6

校正前

-1.33%-3.51%-2.73%-0.20%-3.22%-2.32%-3.53%-3.72%
MIROC6校正后-0.24%-0.27%-0.16%0.61%-0.64%-0.27%-1.51%-0.74%

新窗口打开| 下载CSV


MIROC6全球模式数据由于模式本身的不完备性,以及分辨率较低,由其计算并插值得到夏季大气0 ℃层高度,虽然也能够较好地表现出与探空数据一致的变化趋势(图2),但仍存在偏差(表1)。为保证未来模拟精度,以ERA5再分析资料为观测事实,采用双线性插值方法将MIROC6模式数据插值到0.25°×0.25°的经纬网格上,利用经验性系数校正方法,计算MIROC6每个格点基于ERA5对应格点的多年相对偏差平均值,得到每个格点的夏季大气0 ℃层高度校正系数,进而对MIROC6模式数据进行偏差校正,并以探空站数据为参考对校正后的MIROC6模式数据进行点尺度的评估,结果如表1所示。校正后的MIROC6模式数据能够有效减少偏差,对比校正前有着较大提升。因此,校正后的MIROC6模式数据可用于预估未来青藏高原夏季大气0 ℃层高度的变化趋势。

2.2 1979—2019年青藏高原夏季大气0 ℃层高度的时空变化

图3是由ERA5再分析资料计算得到的高原1979—2019年夏季大气0 ℃层高度平均值的空间分布,可见大气0 ℃层高度以高原中南部(30°~32° N,83.5°~88.5° E)为高值中心,呈纬向带逐渐向四周降低。在高原中南部,大气0 ℃层高度最高值可达到5 972 m,而在帕米尔高原西北部最低仅为4 423 m,祁连山地区均低于5 000 m。

图3

图3   青藏高原夏季大气0 ℃层高度空间分布

Fig. 3   Spatial distribution of the 0 ℃ level height over the Qinghai-Tibet Plateau in summer


在1979—2019年期间,青藏高原夏季大气0 ℃层高度气候倾向率在夏季以及6、7、8月表现出相似的空间分布(图4)。空间上夏季、6月、7月和8月大气0 ℃层高度升高趋势分别在-6~80 m∙(10a)-1、 -10~56 m∙(10a)-1、-7~92 m∙(10a)-1和-16~103 m∙(10a)-1之间。整体上在高原北部大气0 ℃层高度变化存在不规则的纬度地带性特征,北部升高幅度大于南部,在祁连山地区夏季大气0 ℃层高度升高最快,为60 m∙(10a)-1,通过了信度为0.01的显著性检验。以(31° N,85° E)为中心的高原中南部地区呈微弱的下降趋势,但未通过95%信度检验。过去41 a在祁连山地区(36.5°~39.5° N,93.5°~103° E)夏季大气0 ℃层高度按线性趋势升高了约246 m,而在帕米尔高原(35 °~40° N,71°~76° E)仅仅升高了约55 m。

图4

图4   青藏高原夏季(a)和6月(b)、7月(c)、8月(d)大气0 ℃层高度气候倾向率空间分布,红点表示通过了信度为0.01的显著性检验,黑点表示通过了信度为0.05的显著性检验

Fig.4   Spatial distribution of climate tendency rate of the 0 ℃ level height in summer (a), June (b), July (c) and August (d) over the Qinghai-Tibet Plateau, the red dot indicates that it has passed the significance test with the reliability of 0.01, and the black dot indicates that it has passed the significance test with the reliability of 0.05


图1的研究区为轮廓,其中共有4 876个EAR5再分析格点,把这些格点上大气0 ℃层高度进行平均,可反映青藏高原整体夏季大气0 ℃层高度的变化(图5)。在1979—2019 年期间,在夏季和6、7、8月,大气0 ℃层高度的变化趋势分别为27 m∙(10a)-1和17 m∙(10a)-1、35 m∙(10a)-1、30 m∙(10a)-1,按照线性趋势,在过去41 a平均的大气0 ℃层高度分别升高了110 m和68 m、143 m、123 m,其中夏季和7月、8月都通过了信度为0.01的显著性检验。

图5

图5   1979—2019年青藏高原夏季(a)和6月(b)、7月(c)、8月(d)大气0 ℃层高度年际变化趋势

Fig. 5   Interannual variation trend of the 0 ℃ level height in summer (a), June (b), July (c) and August (d) over the Qinghai-Tibet Plateau from 1979 to 2019


2.3 青藏高原夏季大气0 ℃层高度与地面温度的关系

在1979—2019年期间青藏高原夏季地面温度总体呈上升趋势[图6(a)],其空间变化与大气0 ℃层高度变化[图4(a)]相似,都呈现出北快南慢的态势。其中祁连山地区升温速率最高,可达0.67 ℃∙(10a)-1,通过了信度为0.01的显著性检验,相应[图4(a)]中夏季大气0 ℃层高度在祁连山升高也最快。

图6

图6   青藏高原夏季地面温度气候倾向率空间分布(a)以及大气0 ℃层高度和地面温度年际变化曲线(b),红点表示通过了信度为0.01的显著性检验,黑点表示通过了信度为0.05的显著性检验

Fig.6   Spatial distribution of climate tendency rate of summer near-surface air temperature (a) and interannual variation curves of the 0 ℃ level height and near-surface air temperature (b) over the Qinghai-Tibet Plateau, the red dot indicates that it has passed the significance test with the reliability of 0.01, and the black dot indicates that it has passed the significance test with the reliability of 0.05


在青藏高原夏季大气0 ℃层高度比地表高出一般在2 km以内,即高原夏季大气0 ℃层位于大气边界层内。对流层大气能量主要靠吸收来自地表的长波辐射,并通过对流活动使地表能量同上层大气进行能量交换,同时对流层中下层温室气体(主要为水汽)增加也可吸收大量的长波辐射,因此造成对流层中下层大气温度的改变,进而影响到大气0 ℃层高度变化。[图6(b)]显示青藏高原夏季平均地面温度与大气0 ℃层高度随时间的变化特征高度一致,两者之间的相关系数为0.95(P<0.01),通过了信度为0.01的显著性检验,表明青藏高原地区夏季大气0 ℃层高度上升与地面温度持续升高存在显著相关性。通过最小二乘法线性拟合得到以下关系:大气0 ℃层高度=122×地面温度+4380,即夏季地面温度每升高1 ℃,夏季大气0 ℃层高度升高122 m。在1979—2019年期间,高原夏季平均地面温度按线性趋势升高了1.23 ℃,对应的平均夏季大气0 ℃层高度升高了150 m。

2.4 青藏高原夏季大气0 ℃层高度未来变化趋势预估

为预估未来不同共享社会经济路径下青藏高原夏季大气0 ℃层高度变化趋势,采用校正后的MIROC6模式的未来模拟结果。在SSP1-2.6、SSP2-4.5、SSP3-7.0和SSP5-8.5情景下,2020—2100年期间高原夏季大气0 ℃层高度都表现为升高趋势(图7),但不同情景下升高趋势在空间上差别较大。在SSP1-2.6情景下高值中心位于青藏高原中部,而在SSP5-8.5情景下,高值中心位于喜马拉雅山一带,最高值可达99 m∙(10a)-1。4种情景下夏季大气0 ℃层高度升高趋势分别在2~21 m∙(10a)-1、28~40 m∙(10a)-1、61~76 m∙(10a)-1和82~99 m∙(10a)-1之间,均通过了信度为0.01的显著性检验。从低排放情景到高排放情景,升高速率不断增加,如帕米尔高原的气候倾向率从SSP1-2.6情景下的6 m∙(10a)-1上升到SSP5-8.5情景下的86 m∙(10a)-1

图7

图7   在SSP1-2.6(a)、SSP2-4.5(b)、SSP3-7.0(c)、SSP5-8.5(d)情景下大气0 ℃层高度气候倾向率空间分布,所有区域均通过了信度为0.01的显著性检验

Fig. 7   Spatial distribution of climate tendency rate of the 0 ℃ level height under SSP1-2.6 (a), SSP2-4.5 (b), SSP3-7.0 (c) and SSP5-8.5 (d) scenarios, all regions passed the significance test with reliability of 0.01


图8是1979—2100年期间在4种SSP情景下,青藏高原夏季平均大气0 ℃层高度变化。在21世纪前期受辐射强迫差异的影响较小,4种情景下大气0 ℃层高度变化幅度不明显。从2050年开始,4种情景下升高幅度开始出现明显的差异。SSP1-2.6情景下高原夏季大气0 ℃层高度在2060年后基本保持稳定,2081—2100年夏季大气0 ℃层高度比参考时段(1979—2014年)要高265 m;而在SSP2-4.5、SSP3-7.0和SSP5-8.5三种情景下,高原平均夏季大气0 ℃层高度均呈逐年上升趋势,相对于1979—2014年参考时段,到2081—2100年分别升高了394 m、576 m和729 m。

图8

图8   MIROC6模式模拟的青藏高原在不同SSP情景下,1979—2100年夏季大气0 ℃层高度变化

Fig. 8   The variation of the 0 ℃ level height over the Qinghai-Tibet Plateau in summer from 1979 to 2100 under different SSP scenarios simulated by MIROC6 model


3 讨论

在全球变暖背景下,随着地面温度的上升,一方面增加了向上的长波辐射,另一方面对流活动更加活跃,将底层的热量传向高层,致使高层大气温度升高,大气0 ℃层高度也随之上升。近年全球各地区的大气0 ℃层高度整体上都呈升高趋势11-1316-19,相比较,青藏高原地区升高速率更显著,这与高原增温速率是全球平均的2倍6-7相一致。

夏季大气0 ℃层高度的快速升高,对以冰冻圈地貌为主的青藏高原而言,必然会造成消融面积的扩大,引起积雪、冰川和冻土的加速融化。当夏季大气0 ℃层高度超过冰川顶部时,意味着夏季冰川完全处于消融状态,冰川积累区消失并将加速消融。在不同情景下,选取21世纪中叶(2041—2060年)和21世纪末期(2081—2100年)两个时间段,计算青藏高原夏季冰川表面高度低于大气0 ℃高度的面积占冰川总面积的比例(图9)。在SSP1-2.6、SSP2-4.5、SSP3-7.0和SSP5-8.5情景下,2041—2060年处于大气0 ℃层高度以下的冰川面积分别是2014年第二次冰川编目54时的78%、80%、82%和87%,而到2081—2100年,面积比则分别增加到了79%、86%、94%和98%。仅从夏季大气0 ℃层高度变化角度看,在SSP5-8.5情景下,到本世纪末期,预估除帕米尔高原和昆仑山西北部地区外,青藏高原其他地区冰川在夏季都将不存在积累区。

图9

图9   MIROC6模拟的青藏高原在不同SSP情景下2041—2060年(a)和2081—2100年(b)夏季冰川积累区变化

Fig. 9   The changes of glacier accumulation area in summer of 2041—2060 (a) and 2081—2100 (b) were simulated by MIROC6 model under different SSP scenarios


以祁连山冰川为例,更能说明夏季大气0 ℃层高度的变化对青藏高原冰冻圈的影响。祁连山地区(36.5°~39.5° N,93.5°~103° E)在1980s夏季大气0 ℃层高度约为4 839 m,到2010s则升高到5 026 m,仅仅30 a夏季大气0 ℃层高度升高了约187 m,意味着祁连山地区夏季冰川积累区急速缩小,而消融区急速扩大。基于第一次和第二次冰川编目数据,发现近50 a祁连山海拔4 000 m以下的冰川已完全消失,海拔4 350~5 100 m区间冰川面积减少量占冰川面积总损失的84.24%54。祁连山地区冰川顶部海拔一般不超过5 200 m。在SSP1-2.6、SSP2-4.5、SSP3-7.0和SSP5-8.5情景下,夏季大气0 ℃层高度在祁连山地区将分别在2083年、2066年、2052年和2038年达到5 200 m,仅从夏季大气0 ℃层高度变化角度看,夏季整个冰川都处于消融状态,这与从物质能量平衡角度预估的祁连山冰川变化结论一致55。这对受冰川融水影响较大的西北干旱区河流径流而言,短期冰川融水的增加会调节河流径流,而当冰川得不到足够的降水补充时,冰川的退化必然影响到下游的河流径流,需尽快制定规划以应对水资源的变化。

4 结论

本文利用探空数据、ERA5再分析资料和CMIP6的MIROC6模式数据,分析了1979—2100年青藏高原夏季大气0 ℃层高度的时空变化。发现在过去41 a青藏高原夏季大气0 ℃层高度整体呈上升趋势,但空间上呈现北快南慢的态势,在祁连山地区夏季大气0 ℃层高度升高最快。过去41 a在祁连山地区夏季大气0 ℃层高度按线性趋势升高了约246 m,而在帕米尔高原仅仅升高了约55 m。平均而言,夏季地面温度每升高1 ℃,夏季大气0 ℃层高度升高122 m,在1979—2019年期间,高原夏季平均地面温度按线性趋势升高了1.23 ℃,对应的平均夏季大气0 ℃层高度升高了150 m。

预估在SSP1-2.6、SSP2-4.5、SSP3-7.0和SSP5-8.5情景下,2020—2100年期间青藏高原夏季大气0 ℃层高度都表现为升高趋势,这一升高趋势对青藏高原冰川变化造成重大影响。在本世纪中叶(2041—2060年)和末期(2081—2100年)两个时间段,因夏季大气0 ℃层高度的升高,在SSP1-2.6、SSP2-4.5、SSP3-7.0和SSP5-8.5情景下,2041—2060年处于大气0 ℃层高度以下的冰川面积比分别是2014年第二次冰川编目时的78%、80%、82%和87%;而到2081—2100年,面积比则分别增加到了79%、86%、94%和98%,意味着青藏高原地区冰川的积累区在夏季将大幅度缩小。

参考文献

Yao TandongXue YongkangChen Delianget al.

Recent third pole’s rapid warming accompanies cryospheric melt and water cycle intensification and interactions between monsoon and environment: multidisciplinary approach with observations, modeling, and analysis

[J]. Bulletin of the American Meteorological Society, 20191003): 423-444.

[本文引用: 1]

Yao TandongWu GuangjianXu Baiqinget al.

Asian Water Tower change and its impacts

[J]. Bulletin of Chinese Academy of Sciences, 20193411): 1203-1209.

姚檀栋邬光剑徐柏青.

“亚洲水塔”变化与影响

[J]. 中国科学院院刊, 20193411): 1203-1209.

Di YangpingZhang YangjianZeng Huiet al.

Effects of changed Asian water tower on Tibetan Plateau ecosystem: a review

[J]. Bulletin of Chinese Academy of Sciences, 20193411): 1322-1331.

底阳平张扬建曾辉.

“亚洲水塔”变化对青藏高原生态系统的影响

[J]. 中国科学院院刊, 20193411): 1322-1331.

Gao YongpengYao XiaojunLiu Shiyinet al.

Spatial-temporal variation of glacier resources in the Hexi interior from 1956 to 2017

[J]. Journal of Glaciology and Geocryology, 2019416): 1313-1325.

高永鹏姚晓军刘时银.

1956—2017年河西内流区冰川资源时空变化特征

[J]. 冰川冻土, 2019416): 1313-1325.

Zhu RongChen JizuSun Weijunet al.

Ice temperature changing with elevations of a polythermal valley glacier in the northern Tibetan Plateau: a case study on the Laohugou Glacier No.12, Qilian Mountains

[J]. Journal of Glaciology and Geocryology, 2019416): 1292-1301.

[本文引用: 1]

朱荣陈记祖孙维君.

青藏高原北部多温型山谷冰川不同海拔处冰温变化研究——以祁连山老虎沟12号冰川为例

[J]. 冰川冻土, 2019416): 1292-1301.

[本文引用: 1]

Chen DeliangXu BaiqingYao Tandonget al.

Assessment of past, present and future environmental changes on the Tibetan Plateau

[J]. Chinese Science Bulletin, 20156032): 3025-3035.

[本文引用: 2]

陈德亮徐柏青姚檀栋.

青藏高原环境变化科学评估: 过去、现在与未来

[J]. 科学通报, 20156032): 3025-3035.

[本文引用: 2]

Wang BinBao QingHoskins Bet al.

Tibetan Plateau warming and precipitation changes in East Asia

[J]. Geophysical Research Letters, 20083514): L14702.

[本文引用: 2]

Diaz H FGraham N E.

Recent changes in tropical freezing heights and the role of sea surface temperature

[J]. Nature, 19963836596): 152-155.

[本文引用: 1]

Bradley R SKeimig F TDiaz H Fet al.

Recent changes in freezing level heights in the Tropics with implications for the deglacierization of high mountain regions

[J]. Geophysical Research Letters, 20093617): L17701.

[本文引用: 1]

Diaz H FEischeid J KDuncan Cet al.

Variability of freezing levels, melting season indicators, and snow cover for selected high-elevation and continental regions in the last 50 years

[J]. Climatic Change, 20035933-52.

[本文引用: 1]

Huang XiaoyanZhang MingjunWang Shengjieet al.

Variation of 0 ℃ isotherm height and ground temperature in summer in northwest China during the past 50 years

[J]. Acta Geographica Sinica, 2011669): 1191-1199.

[本文引用: 2]

黄小燕张明军王圣杰.

中国西北地区近50年夏季0 ℃层高度及气温时空变化特征

[J]. 地理学报, 2011669): 1191-1199.

[本文引用: 2]

Qiang FangZhang MingjunWang Shengjieet al.

Variation of free-air 0 ℃Isotherm height over tomur peak in the Tianshan Mountain, Xinjiang Uygur Autonomous Region

[J]. Research of Soil and Water Conservation, 2016231): 325-331.

[本文引用: 2]

强芳张明军王圣杰.

新疆天山托木尔峰地区夏季大气0 ℃层高度变化

[J]. 水土保持研究, 2016231): 325-331.

[本文引用: 2]

Zhang GuangxingYang LianmeiYang Qing.

Changing trend and abrupt change of the 0 ℃ level height in summer in Xinjiang from 1960 to 2002

[J]. Journal of Glaciology and Geocryology, 2005273): 376-380.

[本文引用: 1]

张广兴杨莲梅杨青.

新疆43 a来夏季0 ℃层高度变化和突变分析

[J]. 冰川冻土, 2005273): 376-380.

[本文引用: 1]

Zhao AifangZhang MingjunSun Meipinget al.

Changes in 0 ℃ isotherm height of southwest China during 1960-2010

[J]. Acta Geographica Sinica, 2013687): 994-1006.

赵爱芳张明军孙美平.

1960—2010年中国西南地区0 ℃层高度变化特征

[J]. 地理学报, 2013687): 994-1006.

Dong YuanchangWu YaoWang Yuancheng.

The spatial and temporal distribution characteristics of the height of the melting layer over Sichuan Basin and its surroundings in summer

[J]. Plateau and Mountain Meteorology Research, 2018382): 9-14.

董元昌吴遥王源程.

四川盆地及周边地区夏季0 ℃层高度时空分布特征

[J]. 高原山地气象研究, 2018382): 9-14.

Ma XueningZhang MingjunWang Shengjieet al.

Variation in summer 0 ℃ level height and its relationships with temperature and precipitation over the Yellow River Basin

[J]. Resources Science, 20113312): 2302-2307.

[本文引用: 1]

马雪宁张明军王圣杰.

黄河流域夏季0 ℃层高度变化及与地面气温和降水量的关系

[J]. 资源科学, 20113312): 2302-2307.

[本文引用: 1]

Zhou PanpanZhang MingjunWang Shengjieet al.

Variation and its influences of 0 ℃ isotherm height in summer over high Asia

[J]. Plateau Meteorology, 2017362): 371-383.

[本文引用: 1]

周盼盼张明军王圣杰.

高亚洲地区夏季0 ℃层高度变化及其影响特征研究

[J]. 高原气象, 2017362): 371-383.

[本文引用: 1]

Wang LiweiZhang MingjunGao Feng.

Variation of

0 ℃ atmospheric height in the headwaters of Changjiang River in summer during 1977-2010[J]. Plateau Meteorology, 2014333): 769-774.

王立伟张明军高峰.

1977—2010年长江源区夏季大气0 ℃层高度变化

[J]. 高原气象, 2014333): 769-774.

Li GuochangLiu ShixiangZhang Cunjieet al.

Analysis on 0 ℃ level height change in summer over northeast side of Qilian Mountain

[J]. Arid Meteorology, 2006243): 31-34 41.

[本文引用: 1]

李国昌刘世祥张存杰.

祁连山东北侧夏季零度气温层高度变化研究

[J]. 干旱气象, 2006243): 31-34 41.

[本文引用: 1]

Huang XiaoyanWang XiaopingWang Jinsonget al.

Spatio-temporal changes of 0 ℃ isotherm height in China during summer half year of 1970—2012

[J]. Meteorological Monthly, 2017433): 286-293.

[本文引用: 1]

黄小燕王小平王劲松.

1970—2012年夏半年中国大气0 ℃层高度时空变化特征

[J]. 气象, 2017433): 286-293.

[本文引用: 1]

Cheng YingLi DongliangHu Wenchaoet al.

Relationship between glacial thaw of Qilian Mountain and upper temperature

[J]. Plateau Meteorology, 2002212): 217-221.

[本文引用: 1]

程瑛李栋梁胡文超.

祁连山冰川消融与高空气温变化的关系

[J]. 高原气象, 2002212): 217-221.

[本文引用: 1]

Zhang GuangxingSun ShufangZhao Linget al.

The response of the Glacier No.1 to the height change of the 0 ℃ level in summer at the riverhead of the Urumqi River, Tianshan mountains

[J]. Journal of Glaciology and Geocryology, 2009316): 1057-1062.

[本文引用: 1]

张广兴孙淑芳赵玲.

天山乌鲁木齐河源1号冰川对夏季0 ℃层高度变化的响应

[J]. 冰川冻土, 2009316): 1057-1062.

[本文引用: 1]

Wang ShengjieZhang MingjunPepin N Cet al.

Recent changes in freezing level heights in High Asia and their impact on glacier changes

[J]. Journal of Geophysical Research: Atmospheres, 20141194): 1753-1765.

[本文引用: 3]

Zhang JunlanLuo JiWang Rongmei.

Combined analysis of the spatiotemporal variations in snowmelt (ice) flood frequency in Xinjiang over 20 years and atmospheric circulation patterns

[J]. Arid Zone Research, 2021382): 339-350.

[本文引用: 1]

张俊兰罗继王荣梅.

近20 a新疆升温融雪(冰)型洪水频次时空变化及大气环流型分析

[J]. 干旱区研究, 2021382): 339-350.

[本文引用: 1]

Mao WeiyiChen PengxiangShen Yongping.

Characteristics and effects of the extreme maximum air temperature in the summer of 2015 in Xinjiang under global warming

[J]. Journal of Glaciology and Geocryology, 20162): 291-304.

毛炜峄陈鹏翔沈永平.

气候变暖背景下2015年夏季新疆极端高温过程及其影响

[J]. 冰川冻土, 20162): 291-304.

Fu HuaJia LihongXiao Jidonget al.

Classification of snowmelt flood and analysis on its formation causes in the kumalak river basin

[J]. Arid Zone Research, 2011283): 433-437.

[本文引用: 1]

傅华贾丽红肖继东.

阿克苏地区库玛拉克河流域融雪洪水分型及成因

[J]. 干旱区研究, 2011283): 433-437.

[本文引用: 1]

Zhang Guangxing.

The response of annual runoff to the height change at the zero temperature level in summer over Xinjiang

[J]. Acta Geographica Sinica, 2007623): 279-290.

[本文引用: 1]

张广兴.

新疆夏季0 ℃层高度变化对河流年径流量的影响

[J]. 地理学报, 2007623): 279-290.

[本文引用: 1]

Yang PengChen YaningLi Weihonget al.

Analysis of changes in runoff and drying in the Tarim River from 2003 to 2012

[J]. Resources Science, 2015373): 485-493.

杨鹏陈亚宁李卫红.

2003—2012年新疆塔里木河径流量变化与断流分析

[J]. 资源科学, 2015373): 485-493.

Sun GuiliChen YaningLi Weihonget al.

The response of glacial lake outburst floods to climate change in the Yarkant River, Xinjiang

[J]. Journal of Glaciology and Geocryology, 2010323): 580-586.

孙桂丽陈亚宁李卫红.

新疆叶尔羌河冰川湖突发洪水对气候变化的响应

[J]. 冰川冻土, 2010323): 580-586.

Wang YongliYusup·AbdullaMa Hongwuet al.

Response of summer average discharge in the Hotan River to changes in regional 0 ℃ level height

[J]. Advances in Climate Change Research, 200843): 151-155.

[本文引用: 1]

王永莉玉苏甫·阿布都拉马宏武.

和田河夏季流量对区域0 ℃层高度变化的响应

[J]. 气候变化研究进展, 200843): 151-155.

[本文引用: 1]

Zhang JunlanDuan Jianjun.

Analyses on springtime runoff variance of Aksu River and its climate causes

[J]. Plateau Meteorology, 2009282): 465-473.

[本文引用: 1]

张俊岚段建军.

阿克苏河流域春季径流变化及气候成因分析

[J]. 高原气象, 2009282): 465-473.

[本文引用: 1]

Gong WeihuaWang YanguoZhao Chengyiet al.

The hydrological characteristics in Aksu River in lowflow year of 2009 and its effects on the ecological environment of the Tarim River mainstream

[J]. Journal of Glaciology and Geocryology, 2010323): 602-608.

[本文引用: 1]

龚伟华王彦国赵成义.

2009年阿克苏河枯水年水文特征及其对塔里木河干流生态环境的影响

[J]. 冰川冻土, 2010323): 602-608.

[本文引用: 1]

Chen ZhongshengChen YaningLi Weihong.

Response of runoff to change of atmospheric 0 ℃ level height in summer in arid region of Northwest China

[J]. Scientia Sinica Terrae, 20124211): 1770-1780.

[本文引用: 1]

陈忠升陈亚宁李卫红.

中国西北干旱区夏季径流量对大气0 ℃层高度变化的响应

[J]. 中国科学: 地球科学, 20124211): 1770-1780.

[本文引用: 1]

Shang LiHuang YuyingMao Weiyi.

Features of the snow and ice meltwater flood caused by high temperature in the Southern Xinjiang Region during the summer of 2015

[J]. Journal of Glaciology and Geocryology, 2016382): 480-487.

商莉黄玉英毛炜峄.

2015年夏季南疆地区高温冰雪洪水特征

[J]. 冰川冻土, 2016382): 480-487.

Pan ShukunZhang MingjunWang Shengjieet al.

Relationship between streamflow in summer at the headwaters of Urumqi River in the Tianshan mountains and the 0 ℃ isotherm height

[J]. Resources Science, 2012348): 1565-1573.

潘淑坤张明军王圣杰.

天山乌鲁木齐河源区夏季径流量与0 ℃层高度的关系

[J]. 资源科学, 2012348): 1565-1573.

Gong HengruiShi YuFeng Zhimin.

Relationship between the 0 ℃ layer height and the streamflow of the Urumqi River in the period of spring snowmelt

[J]. Arid Zone Research, 2010271): 69-74.

宫恒瑞石玉冯志敏.

春季融雪期0 ℃层高度与乌鲁木齐河径流量的关系

[J]. 干旱区研究, 2010271): 69-74.

Mao WeiyiWu JunChen Chunyan.

Relationship of 0 ℃ level height and summer flood of Aksu River, Xinjiang

[J]. Journal of Glaciology and Geocryology, 2004266): 697-704.

[本文引用: 1]

毛炜峄吴钧陈春艳.

0 ℃层高度与夏季阿克苏河洪水的关系

[J]. 冰川冻土, 2004266): 697-704.

[本文引用: 1]

Zhang YinshengGuo Y.

Variability of atmospheric freezing-level height and its impact on the cryosphere in China

[J]. Annals of Glaciology, 20115258): 81-88.

[本文引用: 2]

Dong LeiZhang MingjunWang Shengjieet al.

The freezing level height in the Qilian Mountains, northeast Tibetan Plateau based on reanalysis data and observations, 1979—2012

[J]. Quaternary International, 201538060-67.

[本文引用: 3]

Huang XiaoyanWang ShengjieWang Jinsonget al.

Spatio-temporal changes in free-air freezing level heights in Northwest China, 1960-2012

[J]. Quaternary International, 2013313/314130-136.

[本文引用: 1]

Zhu ZhiShi ChunxiangGu Junxiaet al.

Ground temperature spatial and temporal variation in Qinghai-Xizang Plateau in recent 10 years

[J]. Science Technology and Engineering, 20202010): 3828-3837.

[本文引用: 1]

朱智师春香谷军霞.

近10 a来青藏高原地表温度时空变化特征分析

[J]. 科学技术与工程, 20202010): 3828-3837.

[本文引用: 1]

Qin YanhuiWu TonghuaLi Renet al.

The applicability of ERA-Interim land surface temperature dataset to map the permafrost distribution over the Tibetan Plateau

[J]. Journal of Glaciology and Geocryology, 2015376): 1534-1543.

秦艳慧吴通华李韧.

ERA-Interim地表温度数据集在青藏高原冻土分布制图应用的适用性评估

[J]. 冰川冻土, 2015376): 1534-1543.

He DongyanTian HongDeng Weitao.

Applicability analysis of three reanalysis surface temperature data over the Tibetan Plateau

[J]. Transactions of Atmospheric Sciences, 2013364): 458-465.

何冬燕田红邓伟涛.

三种再分析地表温度资料在青藏高原区域的适用性分析

[J]. 大气科学学报, 2013364): 458-465.

Li RuiqingShihua Han Boet al.

Preliminary comparison and analyses of air temperature at 2 m height between three reanalysis data-sets and observation in the east of Qinghai-Xiang Plateau

[J]. Plateau Meteorology, 2012316): 1488-1502.

李瑞青吕世华韩博.

青藏高原东部三种再分析资料与地面气温观测资料的对比分析

[J]. 高原气象, 2012316): 1488-1502.

Cao YangChen HongbinLi Junet al.

Characteristic and correlation between surface temperature and 0 ℃ isotherm height derived from ERA-interim reanalysis and radiosonde data

[J]. Plateau Meteorology, 2017366): 1608-1618.

[本文引用: 1]

曹杨陈洪滨李军.

利用再分析与探空资料对0 ℃层高度和地面气温变化特征及其相关性的分析

[J]. 高原气象, 2017366): 1608-1618.

[本文引用: 1]

Zhu JingYuan Huizhen.

Applicability of ERA reanalysis data of land surface temperature in Zhejiang Province

[J]. Meteorological Science and Technology, 2019472): 289-298.

[本文引用: 1]

朱景袁慧珍.

ERA再分析陆面温度资料在浙江省的适用性

[J]. 气象科技, 2019472): 289-298.

[本文引用: 1]

Meng XianguiGuo JunjianHan Yongqing.

Preliminarily assessment of ERA5 reanalysis data

[J]. Journal of Marine Meteorology, 2018381): 91-99.

孟宪贵郭俊建韩永清.

ERA5再分析数据适用性初步评估

[J]. 海洋气象学报, 2018381): 91-99.

Jing HuiZhao LinMeng Xianhonget al.

Effects of shallow residual layer on the development of the convective boundary layer in Naqu area

[J]. Plateau Meteorology, 2020396): 1318-1328.

[本文引用: 1]

景慧赵林孟宪红.

那曲地区浅薄残余层对对流边界层发展的影响

[J]. 高原气象, 2020396): 1318-1328.

[本文引用: 1]

Wang X L.

Accounting for autocorrelation in detecting mean shifts in climate data series using the penalized maximal t or F test

[J]. Journal of Applied Meteorology and Climatology, 2008479): 2423-2444.

[本文引用: 1]

Wang X L.

Penalized maximal F test for detecting undocumented mean shift without trend change

[J]. Journal of Atmospheric and Oceanic Technology, 2008253): 368-384.

[本文引用: 1]

O’Neill B CTebaldi Cvan Vuuren D Pet al.

The scenario model intercomparison project (ScenarioMIP) for CMIP6

[J]. Geoscientific Model Development, 201699): 3461-3482.

[本文引用: 1]

Zhang LixiaChen XiaolongXin Xiaoge.

Short commentary on CMIP6 scenario model intercomparison project(ScenarioMIP)

[J]. Climate Change Research, 2019155): 519-525.

[本文引用: 1]

张丽霞陈晓龙辛晓歌.

CMIP6情景模式比较计划(ScenarioMIP)概况与评述

[J]. 气候变化研究进展, 2019155): 519-525.

[本文引用: 1]

Wei Fengying. Modern climate statistical diagnosis and forecasting technique[M]. 2nd. BeijingChina Meteorological Press2007.

[本文引用: 1]

魏凤英. 现代气候统计诊断与预测技术[M]. 2版. 北京气象出版社2007.

[本文引用: 1]

Sun MeipingLiu ShiyinYao Xiaojunet al.

Glacier changes in the Qilian Mountains in the past half century: based on the revised First and Second Chinese Glacier Inventory

[J]. Acta Geographica Sinica, 2015709): 1402-1414.

[本文引用: 2]

孙美平刘时银姚晓军.

近50年来祁连山冰川变化——基于中国第一、二次冰川编目数据

[J]. 地理学报, 2015709): 1402-1414.

[本文引用: 2]

Duan KeqinYao TandongShi Peihonget al.

Simulation and prediction of equilibrium line altitude of glaciers in the eastern Tibetan Plateau

[J]. Scientia Sinica Terrae, 2017471): 104-113.

[本文引用: 1]

段克勤姚檀栋石培宏.

青藏高原东部冰川平衡线高度的模拟及预测

[J]. 中国科学: 地球科学, 2017471): 104-113.

[本文引用: 1]

/