1. 引言
全球变暖将导致区域水循环的变化,引起水资源时空分布改变,进而影响区域生态环境和社会经济的发展。昆仑山北坡和田地区的水资源系统主要以冰雪融水为基础,在全球气候变化背景下,表现出的脆弱性非常明显 [1]。刘时银 [2] 分析和田地区冰川数为564.0条,面积6812.67 km2,冰储量632.66 km3,分别占新疆冰川总量的27.25%、30.11%和29.36%;比第一次冰川编目,冰川减少了272条,面积减少了601.67 km2。沈永平 [3] 分析,玉龙喀什河冰川融水量为59.3%,喀拉喀什河,冰川融水量为45.7%,雪线高度在5760~6120 m。和田河两条支流出山径流的年内分配极为不均,但是夏季径流比例分别高达80.7%和72.9%,和田河地表水资源主要来源于冰川融化,高度集中在夏季。近年来许多学者在新疆昆仑山北坡地区径流变化趋势研究方面开展了比较全面的研究工作,尤其针对气温、降水和冰川融化与径流变化之间的关系作了系统分析 [4] - [9],针对昆仑山北坡年径流量变化的原因讨论相对较少。应用多年连续监测记录昆仑山北坡年径流量,出山口降水量的变化及昆仑山北坡平原区年平均气温,700~400 hPa气温、年降水量、极端降水日数、高温日数资料,深入开展昆仑山北坡年径流量对气候变化的响应研究,对实现区域水资源可持续利用,维护流域生态环境,促进和田流域绿洲经济发展具有重要的意义。
2. 资料与方法
以和田地区皮山县,和田市、于田、民丰气象站1961~2018年58年气温,降水、蒸发、高温日数资料,平均值为WMO规定的1981~2010年的标准气候平均值。
水文资料为玉龙喀什河,喀拉喀什河、克里雅河、皮山河,1961~2017年57年径流量,平均值为WMO规定的1981~2010年的标准气候平均值。1) 利用极端降水采用百分位阈值法,把和田地区1961~2018年逐年日降水量序列的第90个百分位阈值的30年平均值定义为极端降水事件的阈值。2) 线性趋势倾向系数由最小二乘法求的。3) 5年滑动,累积距平、M-K突变检验法。
3. 结果与讨论
3.1. 和田地区气温变化特征
和田地区年平均气温为12.1℃。气温呈线性上升趋势(图1),上升率0.45℃/10a,高于全国(0.23℃/10a)和全球(0.13℃/10a)同期变化水平。四季平均气温均呈上升趋势,秋季升温速率最大,0.41℃/10a,其次是冬季,为0.39℃/a年。1997年开始和田地区变暖尤其明显,年平均气温连续22年持续偏高。
从各区域情况看,皮山县、和田市、于田县、民丰县各分区年平均气温变化趋势均呈现显著上升趋势,分别上升了0.28℃/10a、0.45℃/10a、0.16℃/10a和0.43℃/10a。
3.2. 高空700~400 hPa气温
张鲜鹤 [10] 分析中国冰川高度结构特征中指出,昆仑山冰川平衡线高度在5582 m,500 hPa的平均高度为5441 gpm左右,沈永平 [3] 分析,玉龙喀什河冰川融水量为59.3%,喀拉喀什河,冰川融水量为45.7%,雪线高度在5760~6120 m。700~400 hPa (3000~7000米)平均气温变化能够反映出高山区的气候变化特征。从和田高空700~400 hPa年平均气温标准化序列曲线(图2)看出和田市高空700~400 hPa年平均气温总体呈上升趋势,上升率0.168℃/10a,气温变化趋势与地面气温变化趋势一致。5年滑动曲线反映出1961~1993年为偏冷阶段,1994~2018年为偏暖阶段。其峰点是2016年(平均气温为−11℃,距平1.2℃),谷点是1989年(平均气温为−13.3℃,距平−1.1℃)。
![](//html.hanspub.org/file/23-2950928x8_hanspub.png?20220916085257271)
Figure 1. Variation of the average discharge annual temper ature over the Hotan region
图1. 和田地区年平均气温
![](//html.hanspub.org/file/23-2950928x9_hanspub.png?20220916085257271)
Figure 2. Variation of the average discharge 700~400 hPa temperature over the Hotan region
图2. 和田高空700~400 hPa年平均气温标准化序列曲线
3.3. 平原区年降水量的变化特征
和田地区平均年降水量为48.2 mm。1961年以来和田地区年降水量呈线性上升趋势,上升率5.1 mm/10a。1961年以来降水量化大致可以分为4个阶段:1961~1986年和1997~2001为偏少阶段,1987~1996年和2002~2018为偏多阶段。降水量最大出现在2010年,为147 mm,最小出现在1985年,为10.6 mm。从各区域的情况看,皮山县、和田市、于田县、民丰县各分区年降水量分别增多6.8 mm/10a、4.5 mm/10a、3.8 mm/10a和5.3 mm/10a,皮山县增多速率最大。
3.4. 极端降水日数变化特征
和田地区1961~2018年和田地区年降水量标准化序列曲线(图3),极端降水日数标准化序列曲线(图4)可以看出,极端降水日数呈增多的趋势,上升率0.5 d/10a。极端降水日数变化大致可以分为4个阶段,1961~1986年和1997~2000年为偏少阶段,1987~1996年和2001~2018年为偏多阶段,极端降水日数最多的年份出现(2010年,为19 d)。
![](//html.hanspub.org/file/23-2950928x10_hanspub.png?20220916085257271)
Figure 3. Variation of the average discharge annual precipitation over the Hotan region
图3. 和田地区年降水量标准化序列曲线
![](//html.hanspub.org/file/23-2950928x11_hanspub.png?20220916085257271)
Figure 4. Variation of the average discharge of annual extreme precipitation over the Hotan region
图4. 和田地区年极端降水日标准化序列曲线
3.5. 高温天气日数变化特征
和田地区1961~2018年高温天气日数标准化序列曲线(图5),可以看出,近58年来和田地区高温天气日数呈线性增加趋势,增加率0.9 d/10a。从各区域情况看,皮山县、和田市、民丰县各分区年平均高温天气日数变化趋势均呈增加趋势,分别增加了0.7 d/10a,0.8 d/10a、2.1 d/10a、于田县−0.2 d/10a。高温天气日数最多的年份出现(1994年,为38 d);最少数的年份出现(1972年,为6 d)。
![](//html.hanspub.org/file/23-2950928x12_hanspub.png?20220916085257271)
Figure 5. Variation of the average discharge of annual mean high temperature days in over the Hotan region
图5. 和田地区年平均高温日数标准化序列曲线
4. 气温,降水突变和趋势检验分析
4.1. 年平均气温突变和趋势检验分析
对和田地区年平均气温进行Mann-Kendall突变检验(见图6)。由图6可知,曲线UF从1961~1968年都是处于0刻度线上下波动,从1969~1996年大多数年份是处于0刻度线以上,年平均气温变化是呈上升的趋势,曲线UF与UB相交点在1997年,且交点在临界线之间,说明1997年发生了突变。1999年增加趋势通过了0.05显著性检验(Z0.05 = 1.96),特别是2016年通过了0.001的显著性检验(Z0.001 = 2.56)。这表明20世纪90年代后期以来,和田地区年平均气温上升趋势是十分显著的。年平均气温Z值为5.63,通过了置信度99%的显著检验。
![](//html.hanspub.org/file/23-2950928x13_hanspub.png?20220916085257271)
Figure 6. M-K nonparametric test of abrupt change for the average temperature over the Hotan region
图6. 和田年平均气温M-K突变检验
4.2. 高空700~400 hPa气温突变和趋势检验分析
对和田市高空700~400 hPa年平均气温进行Mann-Kendall突变检验(见图7)。由图8可知,曲线UF从1961~1984年处于波动式上升阶段,700~400 hPa年平均气温变化呈上升的趋势,曲线UF与UB相交点在1984年,且交点在临界线之间;1984年是突变开始的时间。1988年增加趋势通过了0.05显著性检验(Z0.05 = 1.96),特别是1998年通过了0.001的显著性检验(Z0.001 = 2.56)。700~400 hPa年平均气温从1990年开始大多数年份是处于0刻度线以上,和田市高空700~400 hPa气温Z值为4.23,通过了置信度99%的显著检验。表明和田市高空700~400 hPa年平均气温上升趋势十分显著。
![](//html.hanspub.org/file/23-2950928x14_hanspub.png?20220916085257271)
Figure 7. M-K nonparametric test of 700~400 hPa temperature over the Hotan region
图7. 和田700~400 hPa气温M-K突变检验
4.3. 平原区年降水量突变和趋势检验分析
和田地区平原区年平均降水量进行Mann-Kendall突变检验,如图8可知,1961~1968年UF曲线0刻度线以上,1969~1991年UF曲线0刻度线以下,1992~2002年UF曲线0刻度线上下波动,2003年开始处于0刻度线以上,和田地区年平均降水量变化是呈上升的趋势,曲线UF与UB相交点在2010年,平原区年平均降水量在2010年发生了突变。但是和田站1961~2018年降水Z值为1.69没有通过0.05显著性检验(Z0.05 = 1.96),增加趋势不显著。
![](//html.hanspub.org/file/23-2950928x15_hanspub.png?20220916085257271)
Figure 8. M-K nonparametric test of annualprecipitation over the Hotan region
图8. 和田地区平原区降水量M-K突变检验
4.4. 和田地区昆仑山北坡河流年径流量突变和趋势检验分析
和田地区昆仑山北坡河流年径流量进行Mann-Kendall突变检验,如(图9)可知,1961~1967年和1975~2004年UF线在零线以下振动,和田地区昆仑山北坡河流年径流量偏少阶段,1967~1974年和2005~2017年UF曲线零线以上,和田地区昆仑山北坡河流年径流量偏多阶段。1994年以后上升趋势十分明显,线UF与UB相交点在2009年,昆仑山北坡河流年径流量2009年发生了由少到多的突变。2017年增加趋势通过了0.05显著性检验(Z0.05 = 1.96)。玉龙喀什河2010年发生了突变,Z值为1.79没有通过0.05显著性检验(Z0.05 = 1.96),增加趋势不显著;喀拉喀什河2009年发生了突变,Z值为1.79没有通过0.05显著性检验(Z0.05 = 1.96),增加趋势不显著;克里雅河2007年发生了突变,Z值为3.39通过了0.001的显著性检验(Z0.001 = 2.56),增加趋势显著;皮山河2003年发生了突变,Z值为3.26通过了0.001的显著性检验(Z0.001 = 2.56),增加趋势显著。这表明20世纪90年代中期以来,和田地区昆仑山北坡河流年径流量上升趋势。克里雅河和皮山河Z值置信度99%的显著检验,增加趋势十分显著。
4.5. 昆仑山北坡出山口降水量突变和趋势检验分析
和田地区昆仑山北坡出山口山区年平均降水量进行Mann-Kendall突变检验,如(图9、图10)可知,曲线UF从1961~1986年处于0刻度线上下波动,从1987~2018年处于0刻度线以上,昆仑山北坡出山口降水量呈上升的趋势,1986年为转折点,曲线UF与UB相交点在2002年,昆仑山北坡出山口降水量在2002年发生了突变。但是Z值为1.86没有通过0.05显著性检验(Z0.05 = 1.96),增加趋势不显著。
![](//html.hanspub.org/file/23-2950928x16_hanspub.png?20220916085257271)
Figure 9. M-K nonparametric test of annual runoff on the north slope of Kunlun Mountain over the Hotan region
图9. 和田地区昆仑山北坡河流径流量
![](//html.hanspub.org/file/23-2950928x17_hanspub.png?20220916085257271)
Figure 10. M-K nonparametric test the change of precipitation at the mountain pass on the north slope of Kunlun Mountain over the Hotan region
图10. 出山口降水量 M-K非参数突变检验
5. 和田地区昆仑山北坡年径流量和出山口降水量变化特征
5.1. 和田地区昆仑山北坡河流年径流量变化特征
由图11可见,由于气候变暖,以积雪融水补给为主的和田地区昆仑山北坡河流年径流量变化是呈波动上升的趋势,线性倾向率为19.06 m3∙s−1/10a;年径流年际之间差异波动显著,年最大径流量和年最小径流量分别为781 m3∙s−1 (2010年)、320.8 m3∙s−1 (1965年)。1961~2017年的57a中,有35a为负距平,22a为正距平。自1961年以来和田地区昆仑山北坡年平均流量变化大致可以分为4个阶段:1961~1973年为枯水阶段,1974~1981年为丰水阶段,1982~2000年为枯水阶段,2001~2017为丰水阶段,2010年以来昆仑山北坡河流年径流量增加趋势非常明显。其中,玉龙喀什河增加趋势显著,线性倾向率为26.7 m3∙s−1/10a、喀拉喀什河线性倾向率为23.5 m3∙s−1/10a、克里亚河线性倾向率为18.9 m3∙s−1/10a、皮山河线性倾斜率为7.1 m3∙s−1/10a。
![](//html.hanspub.org/file/23-2950928x18_hanspub.png?20220916085257271)
Figure 11. Variation of the average discharge of annual runoff on the north slope of Kunlun Mountain
图11. 昆仑山北坡河流年径流量
5.2. 和田地区昆仑山北坡出山口降水量变化趋势分析
由图12可见近57年和田地区昆仑山北坡出山口降水量变化的主要特征是年际波动比较剧烈,但总体上呈增加趋势,以9.823 mm∙(10a)−1的倾向率增加。降水量在1961~2001年41a中,有27a为负距平,14a为正距平,2002~2017年16a中有6a为负距平,10a为正距平,但到2010年开始降水量不但增加趋势明显,而且降水变率增大显著。降水量最大出现在2010年,为303 mm,偏多166.5 mm,最小出现在1985年,为45.8 mm,偏少91.6 mm。其中,玉龙喀什河出山口为2.9 mm/10a,克里亚河出山口为13.8 mm/10a、皮山河出山口为12.7 mm/10a。
![](//html.hanspub.org/file/23-2950928x19_hanspub.png?20220916085257271)
Figure 12. Variation of the average discharge of annual precipitation the north slope of Kunlun Mountain
图12. 昆仑山北坡河流出山口年降水量标准化序列曲线
6. 和田地区昆仑山北坡气候变化对径流量的影响
6.1. 温度和降水变化对径流量的影响
和田地区昆仑山北坡河流发源于昆仑山脉,径流均以冰雪融水为主。图12显示,昆仑山北坡径流量,高空700~400 hPa气温和山区降水量标准化变量累积距平曲线变化具有明显的同位相特征。昆仑山北坡的径流量累积距平曲线1961~2000年波动式下降,属于枯水期,2000年为拐点,2001~2017年累积距平曲线持续上升,属于丰水期;高空700~400 hPa气温累积距平在1961~1988年期间呈波动式下降趋势,1989年为拐点,1990~2017年期间呈上升趋势;上升趋势显著。昆仑山北坡出山口降水量累积距平在1961~1970年期间波动式下降,1971~1974年上升,1975~1986年下降趋势十分明显;1987~1993年累积距平曲线波动式上升,1994~2001年期间波动式下降,2001为转折点,2002年开始上升趋势十分显著。
6.2. 昆仑山北坡径流量,高空700~400 hPa气温和山区降水突变前后变化分析
采用M-K非参数检验法,累积距平对昆仑山北坡径流量,高空700~400 hPa气温和山区降水进行突变分析。昆仑山北坡的径流量2009年发生突变,高空700~400 hPa气温在1984年发生突变、山区降水量在2002年发生突变。径流量突变前(1961~2009年)平均为414.9 m3∙s−1,突变后(2010~2017年)平均为547.7 m3∙s−1,突变后比突变前增多了132.8 m3∙s−1;山区降水量突变前(1961~2002年)平均为123.8 mm,突变后(2003~2017年)平均为172.4 mm,突变后比突变前增多了48.6 mm;高空700~400 hPa气温突变前(1961~1984年)平均为−12.5℃,突变后(1985~2017年)平均为−12.0℃,突变后比突变前上升了0.5℃。
6.3. 昆仑山北坡径流量,高空700~400 hPa气温、山区降水异常年分析
采用了标准化距平大于1.0的年份为偏多(偏高),距平小于−1.0的年份为偏少(偏低);距平大于2.0的年份为异常偏多(偏高),距平小于−2.0的年份为异常偏少(偏低)依此选定昆仑山北坡径流量,高空700~400 hPa气温、山区降水典型的异常偏多(偏高),异常偏少(偏低)年。陈亚宁 [1] 最新研究结果显示,西北干旱区河川径流量的长时期序列检测发现,一些冰川融水补给较大的河流在1994年前后普遍出现了一次径流“突变型”增加的过程与气候异常年份较一致。昆仑山北坡径流量5个异常偏多年中4个出现在1994年以后,9个偏多年中5个出现在1994年以后,出现在1994年以前;1个异常偏少年份(1965年)和9个偏少年中7个出现在1994年以前,1994年以后只出现了2个偏少年。高空700~400 hPa气温1个异常偏高年份(2016年),2个异常偏低年份(1962,1989);7个偏高年份1994年以后;13个偏低年中11个出现在1994年以前,2个偏低年出现1994年以后。山区降水量2个异常偏多年份(2010、2017年),没有出现异常偏少年份;8个偏多年中6个出现在1987年以后,2个出现在1987年以前,6个偏少年中4个出现在1987年以前,2个出现在1987年以后。
6.4. 突变前后昆仑山北坡径流量,高空700~400 hPa气温、平原区降水量、山区降水量、平原区极端降水日数相关分析
采用M-K非参数检验法,累积距平对昆仑山北坡径流量,高空700~400 hPa气温和山区降水进行突变分析。昆仑山北坡的径流量2009年发生突变,高空700~400 hPa气温在1984年发生突变、山区降水量在2002年发生突变、平原区降水量在2010年发生突变、平原区极端降水日数1986年发生突变。700~400 hPa气温与昆仑山北坡径流量进行Pearson相关性分析;突变前(1961~1983年)相关系数为0.233,相关不显著;突变后(1984~2017年)相关系数为0.377,通过了0.05的显著性水平检验。尤其是400 hPa气温与昆仑山北坡径流量的相关最明显,突变前相关系数0.339,相关不显著,突变后相关系数0.439,通过了0.01的显著性水平检验。昆仑山北坡径流量与平原区降水量突变前(1961~2009年)负相关,相关系数为−0.005,相关不显著;突变后(2010~2017年)正相关,相关系数0.806,通过了0.05的显著性水平检验。昆仑山北坡径流量与山区降水量突变前(1961~2001年)负相关,相关系数为−0.01,相关不显著;突变后(2002~2017年)正相关,相关系数0.687,通过了0.05的显著性水平检验。昆仑山北坡径流量与平原区极端降水日数突变前(1961~1986年)正相关,相关系数为0.105,相关不显著;突变后相关系数0.719,通过了0.05的显著性水平检验。昆仑山北坡径流量与高空700~400 hPa气温通过了0.01的显著性检验,昆仑山北坡径流量与山区,平原降水量通过了0.05信度水平。这说明,对昆仑山北坡径流量来说,高空700~400 hPa气温是近57年来径流量增长的最主要因素。
7. 结果与讨论
7.1. 和田地区昆仑山北坡流域暖干向暖湿的转型明显
施雅风等 [11] 对西北近年来的降水、径流、洪水、湖泊等变化分析,发现1987年以来西北西部和中部地区气候发生了由暖干向暖湿的转型。和田地区年平均气温上升率0.45℃/10a,高于全国(0.23℃/10a)和全球(0.13℃/10a)同期变化水平。1997年以来和田地区变暖尤其明显,年平均气温连续22年持续偏高;700~400 hPa年平均气温上升率0.168℃/10a,平原区降水量上升率5.1 mm/10a、山区降水量9.8 mm/10a、极端降水日数上升率0.5 d/10a、高温天气日数上升率0.9 d/10a、昆仑山北坡河流年径流量上升率为19.06 m3∙s−1/10a。和田地区昆仑山北坡流域气候正在暖干向暖湿的转型明显。
7.2. 冰川和径流对气候异常的响应
胡凡盛 [12] 分析1976~2016年西昆仑地区冰川变化遥感监测指出西昆仑地区冰川面积退缩了892.38 km2,退缩速率0.28%∙a−1。气温和降水分别影响冰川的消融和积累过程,是冰川发育气候条件最主要的因素。过去57年来,和田地区年平均气温,700~400 hPa年平均气温和山区,平原区降水都有所增长,气温升高有利于冰川消融可以促使冰川退缩,降水的增加有利于冰川的积累。高空700~400 hPa气温的异常变化与西昆仑地区冰川面积退缩有高度的相应,700~400 hPa气温偏高年份(2016年),7个偏高年份(1998、1999、2006、2007、2010、2013,2015)。400 hPa气温与昆仑山北坡径流量的相关最明显,突变前相关系数0.339,相关不显著,突变后相关系数0.439,通过了0.01的显著性水平检验。颜伟 [13] 分析MODIS数据的2000~2013年西昆仑山玉龙喀什河流域积雪面积变化研究指出低,高山区平均、最大积雪面积呈微弱增加趋势,高山区最小积雪面积显著增加,倾向率为65.877 km2/a;秋季,2009、2010年是低山区积雪面积的峰值年份。昆仑山北坡径流量4个异常偏多年出现在(2006、2010、2012、2013),4个偏多年(2005、2015、2016、2017)。径流量异常偏多与高山区最小积雪面积显著增加高度相应,2009年秋季低山区积雪面积的峰值年与2010年昆仑山北坡诸河夏季特大洪水相对应,高空400 hPa气温是昆仑山北坡冰川消融,冰川面积退缩和径流量增多的最主要因素。
8. 结论
1) 1961~2017年,和田昆仑山北坡河属冰雪融水补给型河流,昆仑山北坡径流量以19.06 m3∙s−1∙(10a)−1的倾向率递增,高空700~400 hPa气温以0.168℃∙(10a)−1倾向率增加,山区降水夏季流量以0.186 m3∙s−1∙(10a)−1的倾向率递增、山区降水量以9.823 mm∙(10a)−1的倾向率增加。
2) 昆仑山北坡的径流量2009年发生突变,玉龙喀什河2010年发生了突变,喀拉喀什河2009年发生了突变,克里雅河2007年发生了突变,皮山河2003年发生了突变。高空700~400 hPa气温在1984年发生突变、山区降水量在2001年发生突变。突变前后昆仑山北坡径流量,高空700~400 hPa气温、山区降水量差异明显。
3) 昆仑山北坡径流量与700~400 hPa气温,山区和平原降水量、平原区极端降水日数进行Pearson相关性分析。昆仑山北坡径流量与各要素突变前后差异较大。昆仑山北坡径流量与高空700~400 hPa气温,山区、平原降水量、平原区极端降水日通过了0.05的显著性检验;尤其是400 hPa气温与昆仑山北坡径流量的相关最明显,突变前相关系数0.339,相关不显著,突变后相关系数0.439,通过了0.01的显著性水平检验。这说明,对昆仑山北坡径流量来说,高空400 hPa气温是近57年来昆仑山北坡径流量增长的最主要因素。
4) 和田地区年平均气温上升率0.45℃/10a,高于全国(0.23℃/10a)和全球(0.13℃/10a)同期变化水平。1997年以来和田地区变暖尤其明显,年平均气温连续22年持续偏高。700~400 hPa年平均气温上升,平原区降水量、山区降水量、极端降水日数、高温天气日数增加、径流量增多。说明和田地区昆仑山北坡流域气候正在暖干向暖湿的转型明显。
5) 高空700~400 hPa气温的异常变化与西昆仑地区冰川面积退缩有高度的相应,昆仑山北坡径流量对高空700~400 hPa气温和山区降水量的变化都有很好的响应。径流量异常偏多与高山区最小积雪面积显著增加高度相应,2009年秋季低山区积雪面积的峰值年与2010年昆仑山北坡诸河夏季特大洪水相对应,高空400 hPa气温是昆仑山北坡冰川消融,冰川面积退缩和径流量增多的最主要因素。
基金项目
新疆少数民族科技人才特殊培养计划基金项目(2017D03026)唯一资助。
NOTES
*通讯作者。