王 濤 伊麗努爾·阿力甫江 李思穎 阿依謝姆古麗·孜比不拉 張心蕊 高騰飛

1 新疆大學(xué)地理與遙感科學(xué)學(xué)院，烏魯木齊 830017 2 新疆克州氣象局，阿圖什 845350

提 要：基于東帕米爾高原3個(gè)氣象站近58年(1961—2018年)逐月降水量資料，借助創(chuàng)新性趨勢(shì)分析(innovative trend analysis,ITA)法、Mann-Kendall(M-K)非參數(shù)檢驗(yàn)和Sen’s斜率估計(jì)法分析了東帕米爾高原降水量的變化趨勢(shì)。結(jié)果表明：三種分析方法得出的結(jié)論具有一致性，近年來東帕米爾高原降水量呈增加趨勢(shì)。Sen’s斜率估計(jì)法分析得出有69.44%的月份降水呈增加趨勢(shì)，11.11%的月份呈減少趨勢(shì)；M-K法分析得出有75%的月份降水呈增加趨勢(shì)，22.22%的月份呈減少趨勢(shì)；ITA法分析得出有94.44%的月份降水呈增加趨勢(shì)，5.56%的月份呈減少趨勢(shì)。三種方法在分析過程中也存在著差異性，如利用Sen’s斜率估計(jì)法分析塔什庫爾干氣象站點(diǎn)時(shí)有5個(gè)月無變化趨勢(shì)，與其他兩種方法分析結(jié)果差異較大。而ITA和M-K 法較為準(zhǔn)確地辨識(shí)時(shí)間序列的變化趨勢(shì)，且不受異常值影響；三種分析方法中，顯著性趨勢(shì)最強(qiáng)的是ITA法，其次是M-K 法，Sen’s斜率估計(jì)法表現(xiàn)較為一般。其中，ITA分析法有63.89% 的月份通過了0.10的顯著性水平檢驗(yàn)，M-K有25% 的月份通過了0.10的顯著性水平檢驗(yàn)。相較于傳統(tǒng)分析方法，ITA的優(yōu)勢(shì)在于，它能夠給出不同月份不同級(jí)別降水量增加或減少的范圍和趨勢(shì)，所以可以有效地識(shí)別傳統(tǒng)檢驗(yàn)無法確定的隱藏的微趨勢(shì)，更有助于分析水文氣象方面數(shù)據(jù)系列的變化趨勢(shì)。

引 言

近百年來全球氣候變化對(duì)各地區(qū)降水產(chǎn)生了顯著的影響，受影響最為顯著的是干旱半干旱區(qū)(馬柱國(guó)和符淙斌，2005)。IPCC(2014)的第五次評(píng)估報(bào)告表明，陸地上越來越多的地區(qū)出現(xiàn)降水量的增長(zhǎng)。受全球氣候變化影響，區(qū)域間氣候差異較大，整體氣溫和降水變化與局部地區(qū)氣溫和降水變化呈現(xiàn)出的趨勢(shì)可能不同。而我國(guó)氣候復(fù)雜多變，氣溫和降水變化可能與其他地區(qū)呈現(xiàn)出差異性，尤其中國(guó)西北干旱半干旱地區(qū)近年來降水變化顯著。目前關(guān)于氣候變化的研究，丁一匯(2008)、姚世博等(2017)、嚴(yán)中偉等(2020)主要是從氣候要素，如氣溫、降水變化特征來分析氣候變化的趨勢(shì)和原因。左洪超等(2004)、王遵婭等(2004)、張先恭和李小泉(1982)、丁一匯和戴曉蘇(1994)研究表明中國(guó)氣候變化與全球氣候變化趨勢(shì)基本一致。西北地區(qū)的降水變化也引起了眾多學(xué)者的關(guān)注(楊金虎等，2006；陳冬冬和戴永久，2009；王澄海等，2021)。楊陽等(2022)得出我國(guó)西北干旱半干旱地區(qū)降水量呈增加趨勢(shì)。王寶靈(1997)得出西北大部分地區(qū)大氣可降水量呈增加趨勢(shì)，新疆等地增加趨勢(shì)明顯。黃小燕等(2018)也分析出了西北地區(qū)6月降水量1961—1990年間明顯遞增。

新疆是西部面積最大的省份，研究新疆及其周邊地區(qū)近年來降水量變化具有重大意義。左洪超等(2004)研究得出中國(guó)新疆地區(qū)降水量變化最為明顯。降水量的增加在范圍和強(qiáng)度上與國(guó)內(nèi)其他地區(qū)相比最大(熱孜宛古麗·麥麥提依明等，2016；楊維西，2008)，增濕趨勢(shì)十分明顯(崔彩霞，2001；姜大膀等，2009)。新疆腹地的和田地區(qū)1961—2000年降水量總體呈增多趨勢(shì)(玉蘇甫·阿布都拉等，2006)。位于新疆省西南邊緣的東帕米爾高原氣候變化也明顯受到影響，降水量的不斷變化可能會(huì)形成洪水和極端天氣事件，并且對(duì)生態(tài)環(huán)境產(chǎn)生了較大影響。毛煒嶧等(2006)分析得出1961—2015年帕米爾高原的年降水量呈顯著線性增加趨勢(shì)。努爾帕提曼·買買提熱依木和帕爾哈提·阿不都拉(2015)分析得出1960—2014年塔什庫爾干縣年降水量呈緩慢增加的趨勢(shì)。關(guān)于東帕米爾高原降水量變化研究較少，研究方法集中于M-K檢驗(yàn)法、一元線性趨勢(shì)分析法和累積距平法等，基于創(chuàng)新性趨勢(shì)分析(innovative trend analysis，ITA)法的研究還未實(shí)踐。ITA法首先是對(duì)數(shù)據(jù)處理做出直觀的圖像變化趨勢(shì)，圖像將數(shù)據(jù)點(diǎn)分為上半三角增加區(qū)、下半三角減少區(qū)和無明顯變化區(qū)，且將圖像縱向分為“低”組、“中”組和“高”組，能夠非常直觀地觀察出降水量在不同階段不同范圍的變化趨勢(shì)(黃丹等，2018)。它彌補(bǔ)了M-K只能檢驗(yàn)單調(diào)變化趨勢(shì)的不足；它更易于理解、計(jì)算，分析結(jié)果更精確和更清楚，表現(xiàn)在于：ITA不需要假設(shè)，保證了對(duì)趨勢(shì)存在性的直觀檢查可能性；它可以用于給定時(shí)間序列中的詳細(xì)趨勢(shì)的確定和分類解釋，用數(shù)字觀察可視化趨勢(shì)之外的變化，給出的D值(表征趨勢(shì)增減性)有助于確定最小值隱藏的趨勢(shì)，還可以防止其他誤導(dǎo)性的判斷(如無限變化率等)，而圖形檢驗(yàn)可以直觀地檢驗(yàn)出單調(diào)和非單調(diào)的變化趨勢(shì)，允許對(duì)三部分做圖形化的評(píng)估。因此，本研究利用ITA法研究分析東帕米爾高原塔什庫爾干、吐爾尕特、烏恰3個(gè)氣象站58年(1961—2018年)的月降水變化趨勢(shì)。對(duì)了解東帕米爾高原近年來降水變化特征具有一定意義。

1 研究區(qū)概況

東帕米爾高原位于中亞東南部、中國(guó)新疆維吾爾自治區(qū)的西南邊緣，38°～41°N、73°～76°30′E。北起天山、東以木吉谷地和塔什庫爾干谷地為界，西與昆侖山相接、南至喀拉昆侖山。高原海拔為4000～7700 m，擁有許多高峰，最高峰可達(dá)7719 m，平均海拔在4500 m以上，主要山峰均在6000 m以上。

東帕米爾高原屬嚴(yán)寒的強(qiáng)烈大陸性高山氣候，冬季漫長(zhǎng)，絕對(duì)最低氣溫達(dá)-50℃，歷史上布侖戈?duì)?953年最低氣溫曾低達(dá)-53.3℃，歷年最高氣溫不超過20℃，年降水量?jī)H為75～100 mm，在喀拉湖盆地更是少至30 mm。本區(qū)處在西風(fēng)和西南季風(fēng)冷濕氣流進(jìn)入青藏高原西北邊緣的通道山區(qū)，因而較為寒冷。由于本區(qū)地處亞歐大陸腹地，遠(yuǎn)離海洋，降水較少，故氣候寒冷且干燥(曾磊等，2013)。多年來極端升溫過程造成該區(qū)域夏季融雪(冰)洪水發(fā)生風(fēng)險(xiǎn)提高，并且極端降水形成的混合型洪水及其衍生地質(zhì)災(zāi)害風(fēng)險(xiǎn)也在加劇(毛煒嶧等，2019)，地質(zhì)災(zāi)害風(fēng)險(xiǎn)的加劇越來越影響人類正常的生產(chǎn)生活，所以研究此區(qū)域近年來氣候變化有一定的實(shí)際意義。

2 研究方法

2.1 Mann-Kendall檢驗(yàn)法

Mann-Kendall(M-K)趨勢(shì)檢驗(yàn)法是一種非參數(shù)檢驗(yàn)，用于檢測(cè)氣象水文時(shí)間序列的趨勢(shì)和突變檢驗(yàn)(Kendall，1938；Mann，1945)。該方法已被世界氣象組織應(yīng)用于研究處理氣候數(shù)據(jù)的時(shí)間變化。其優(yōu)點(diǎn)是不需要樣本遵從一定的分布，也不受少數(shù)異常值的干擾，檢驗(yàn)?zāi)芰?yōu)于參數(shù)檢驗(yàn)，并且無須事先假定樣本的統(tǒng)計(jì)分布，人為性少，定量程度高，計(jì)算較簡(jiǎn)便。α為顯著性水平，本文取α=0.10。

(1)

(2)

式中：xi和xj是在i和j時(shí)刻的數(shù)據(jù)值，n表示數(shù)據(jù)集的長(zhǎng)度(數(shù)據(jù)集中數(shù)據(jù)的個(gè)數(shù))。而S的正值表示增加趨勢(shì)，負(fù)值表示減少趨勢(shì)(Mann，1945；黃丹等，2018；魏鳳英，2007)。

(3)

式中：m是并列組數(shù)(時(shí)間序列中的相等數(shù)據(jù))，摘要符號(hào)(∑)表示所有并列組的總和。tk是第m組中觀察值的數(shù)目(李紅梅等，2008；艾雅雯等，2020)。用式(3)計(jì)算時(shí)間序列方差后，標(biāo)準(zhǔn)Zc值根據(jù)式(4)計(jì)算，Zc為正表示該序列具有增加趨勢(shì)，而Zc為負(fù)表示該序列呈減少趨勢(shì)(魏鳳英，2007)。

(4)

2.2 Sen’s斜率估計(jì)法

Sen’s斜率估計(jì)法計(jì)算序列的斜率β，β表示序列的平均變化率和序列的趨勢(shì)，當(dāng)β>0時(shí)，序列呈上升趨勢(shì)；β=0時(shí)，序列趨勢(shì)不明顯；β<0時(shí)，序列呈下降趨勢(shì)。估計(jì)n個(gè)樣本中N對(duì)數(shù)據(jù)的趨勢(shì)斜率：

(5)

式中：xj和xi分別是第j個(gè)和i個(gè)時(shí)間的樣本序列值(j>i)，Median為取中值函數(shù)。

2.3 創(chuàng)新性趨勢(shì)分析(ITA)法

圖1 以創(chuàng)新性趨勢(shì)分析(ITA)法為例：(a)非單調(diào)趨勢(shì)，(b)3組數(shù)據(jù)點(diǎn)分布趨勢(shì)

(5)

3 分析結(jié)果

3.1 M-K法分析結(jié)果

選取東帕米爾高原3個(gè)代表性的氣象站(塔什庫爾干、吐爾尕特和烏恰)58年(1961—2018年)的逐月降水資料，采用M-K趨勢(shì)檢驗(yàn)3個(gè)氣象站的降水量進(jìn)行趨勢(shì)分析，其計(jì)算結(jié)果見表1。

從表1可知：

表1 1961—2018年東帕米爾高原3個(gè)氣象站點(diǎn)月降水量變化(單位：mm)

(1)塔什庫爾干站6月和8月降水量的Zc值>0，分別為2.119和2.361，均通過了0.05的顯著性水平檢驗(yàn)，呈顯著增加趨勢(shì)。而其他月份均未通過0.10顯著性水平檢驗(yàn)，呈增加或減少趨勢(shì)。其中3、5、10和11月的Zc值均>0，呈增加趨勢(shì)。而其他月份均未通過0.10顯著性水平檢驗(yàn)，占33.333%；而1、2、4、7、9和12月的Zc值均<0，呈減少趨勢(shì)，約占50%。

(2)吐爾尕特站降水量除了1月(Zc=-0.085)呈減少趨勢(shì)外，2—12月的Zc值都>0，呈增加趨勢(shì)，約占91.667%。其中，4月(Zc=1.733)通過了0.10的顯著性水平檢驗(yàn)，3月(Zc=2.119)、6月(Zc=2.197)和10月(Zc=2.348)均通過了0.05的顯著性水平檢驗(yàn)，占25%；而12月(Zc=2.871)通過了0.01的顯著性水平檢驗(yàn)，呈顯著增加趨勢(shì)。

(3)烏恰站的降水量在6月和7月(Zc=2.295和Zc=2.439)，通過了0.05的顯著性水平檢驗(yàn)，呈顯著增加趨勢(shì)，占16.667%。除了3月(Zc=-0.268)為減少趨勢(shì)，該站1、2、4、5月和8—12月這9個(gè)月均沒通過0.10的顯著性水平檢驗(yàn)，但Zc值都>0，因此說明該站在這10個(gè)月降水量呈增加趨勢(shì)，占83.333%。

基于M-K檢驗(yàn)法3個(gè)站點(diǎn)的多年降水量中夏季(6—8月)增加趨勢(shì)均顯著，對(duì)全年降水貢獻(xiàn)最大；秋季(9—11月)和冬季(12月和次年1—2月)次之；春季變化趨勢(shì)不明顯，甚至塔什庫爾干站和吐爾尕特站1月和2月與烏恰站3月呈減少趨勢(shì)。與曾磊(2012)分析的1957—2010年東帕米爾高原不同季節(jié)降水量變化趨勢(shì)一致：“研究區(qū)不同季節(jié)降水量除春季外，夏季、秋季和冬季表現(xiàn)出一定的升高趨勢(shì)并且夏季降水量增幅最為明顯”。

3.2 Sen’s斜率估計(jì)法趨勢(shì)分析結(jié)果

由表1可知，塔什庫爾干站在1—2月、9月、11—12月β為0，表明沒有變化趨勢(shì)；3月、5—6月、8月和10月β為正值，表明降水量呈增加趨勢(shì)；4月和7月β為負(fù)值，表明降水量呈減少趨勢(shì)。

吐爾尕特站在1月β為0，表明沒有變化趨勢(shì)；2—12月β為正值，表明降水量呈增加趨勢(shì)。而烏恰站在2月和11月β為0，表明沒有變化趨勢(shì)；3月β為負(fù)值，呈減少趨勢(shì)；1月、4—12月β為正值，表明呈增加趨勢(shì)?？傊?，58年來東帕米爾高原降水量變化趨勢(shì)不明顯。這與胡凡盛(2018)研究結(jié)果相似，東帕米爾高原1976—2016年降水量線性擬合系數(shù)低，整體呈緩慢增加趨勢(shì)。

3.3 ITA法趨勢(shì)分析結(jié)果

圖2～圖4顯示1∶1(45°)線與數(shù)據(jù)點(diǎn)的位置關(guān)系，由圖可知逐月降水量的變化趨勢(shì)。由圖2可知，塔什庫爾干站3月的數(shù)據(jù)點(diǎn)均落在1∶1線的上方(在增加區(qū)域)，整體呈增加趨勢(shì)；且數(shù)據(jù)點(diǎn)落在“低”組和“中”組的增加區(qū)域，就說明這3個(gè)月降水量的增加是在0～20、0～40和0～60 mm范圍內(nèi)增加。在1月大部分?jǐn)?shù)據(jù)點(diǎn)落在1∶1線的下方(在減少區(qū)域)，且數(shù)據(jù)點(diǎn)落在“低”組和“高”組的減少區(qū)域；說明1月降水量的減少是在0～30 mm范圍內(nèi)減少。

從圖2可知，塔什庫爾干站降水量在2、3、5、6、8、11和12月呈增加趨勢(shì)，這些月的數(shù)據(jù)點(diǎn)落在1∶1線上和1∶1線的上方，落在“低”組和“中”組內(nèi)，增加范圍在0～60 mm內(nèi)。1月降水量呈顯著減少趨勢(shì)，通過了0.10的顯著性水平檢驗(yàn)，1月在“低”組和“高”組減少趨勢(shì)顯著，減少范圍在0～30 mm。4月大部分?jǐn)?shù)據(jù)落在1∶1線上，說明無變化趨勢(shì)，但有些數(shù)據(jù)落在1∶1線下的“低”組的減少區(qū)域，也有一些數(shù)據(jù)落在“中”組的增加區(qū)域。7、9和10月部分?jǐn)?shù)據(jù)落在1∶1線的上方，部分?jǐn)?shù)據(jù)落在1∶1線的下方，上方是在“中”組的增加區(qū)域，下方是在“低”組的減少區(qū)域，這種情況就無法正確閱讀出該月的降水量變化趨勢(shì)，必須結(jié)合ITA指數(shù)進(jìn)行分析。從表3 可知，2、3、5、6、8、11和12月，7個(gè)月顯著水平檢驗(yàn)均通過了0.10，降水量呈顯著增加趨勢(shì)。

圖2 塔什庫爾干1961—2018年1—12月ITA趨勢(shì)分析結(jié)果

如圖3所示，吐爾尕特站點(diǎn)，1月、3—6月和10—12月這8個(gè)月均通過了0.10顯著性水平檢驗(yàn)，呈顯著增加趨勢(shì)。1月一部分?jǐn)?shù)據(jù)點(diǎn)落在“低”組的減少區(qū)域，但十分靠近1∶1線，還有一部分?jǐn)?shù)據(jù)點(diǎn)落在“低”組的增加區(qū)域，但距離1∶1線較遠(yuǎn)，且通過了0.10顯著性水平檢驗(yàn)，所以1月降水量是呈顯著增加的。3—6月和10—12月這7個(gè)月的數(shù)據(jù)點(diǎn)大部分都落在“低”組和“高”組的增加區(qū)域，所以這些月降水量是呈顯著增加的。7月的“低”組的增加區(qū)域和“中”組的減少區(qū)域，分別在0～50和 0～100 mm內(nèi)減少和增加。

圖3 同圖2，但為吐爾尕特

如圖4所示，烏恰站點(diǎn)1月、5—10月和12月全部都通過了0.10顯著性水平檢驗(yàn)，呈顯著增加趨勢(shì)。5月和7—8月數(shù)據(jù)點(diǎn)幾乎全部都落在在“低”組、“中”組和“高”組的增加區(qū)域內(nèi)，分別在0～90、0～120和0～150 mm范圍內(nèi)增加。6、10和12月的數(shù)據(jù)點(diǎn)大部分都落在“低”組的增加區(qū)域內(nèi)，少部分?jǐn)?shù)據(jù)點(diǎn)落在“中”組的增加區(qū)域內(nèi)。9月則是大部分?jǐn)?shù)據(jù)點(diǎn)落在“低”組的增加區(qū)域內(nèi)，少部分落在“中”組的減少區(qū)域內(nèi)，但是通過了0.10顯著性水平檢驗(yàn)，所以也呈顯著增加趨勢(shì)。11月大部分?jǐn)?shù)據(jù)落在1∶1線的下方，落在“低”組和“中”組的減少區(qū)域內(nèi)，所以說該月呈顯著減少趨勢(shì)。

圖4 同圖2，但為烏恰

表2 東帕米爾高原3個(gè)氣象站的ITA指數(shù)

3.4 三種趨勢(shì)分析法的比較

通過ITA和M-K法進(jìn)行對(duì)比得出，3個(gè)氣象站在36個(gè)月(每個(gè)站點(diǎn)有12個(gè)月尺度降水量)的降水量中變化顯著的月份分別是23個(gè)月和9個(gè)月，其中ITA顯著趨勢(shì)為64%，M-K顯著趨勢(shì)為25%。Sen’s斜率估計(jì)檢驗(yàn)36個(gè)月的降水量，其中有8個(gè)月斜率為0，表明無變化趨勢(shì)，而剩余月份的斜率估計(jì)與Zc值的正(負(fù))變化情況一致，表明Sen’s斜率估計(jì)檢驗(yàn)和M-K法有較高的擬合度。通過M-K和Sen’s斜率估計(jì)法與ITA法對(duì)比分析，塔什庫爾干站除6月和8月之外所有月份的降水量沒有顯著變化趨勢(shì)；而通過ITA方法，降水量在2、3、5、6、8、11和12月呈顯著增加趨勢(shì)。通過M-K和Sen’s斜率估計(jì)檢驗(yàn)法分析吐爾尕特氣象站點(diǎn)3、4、6、10和12月的降水量呈顯著增加趨勢(shì)，而通過ITA分析法得知1月、3—6月和10—12月的降水量均呈顯著增加趨勢(shì)；烏恰氣象站點(diǎn)，通過M-K和Sen’s斜率估計(jì)檢驗(yàn)法得知6月和7月呈顯著增加趨勢(shì)，而其他10個(gè)月無顯著趨勢(shì)，通過ITA分析法得知1月、5—10月和12月呈顯著增加趨勢(shì)，而11月呈顯著減少趨勢(shì)。

4 結(jié) 論

本文分別利用M-K檢驗(yàn)法、Sen’s斜率估計(jì)和ITA法，針對(duì)東帕米爾高原3個(gè)代表氣象站點(diǎn)，1961—2018年逐月降水資料，分析得到東帕米爾高原的降水變化整體趨勢(shì)。結(jié)果如下：

(1)1961—2018年，三種方法均得出東帕米爾高原年降水量呈現(xiàn)出顯著的增加趨勢(shì)，與尚華明等(2021)分析東帕米爾高原近幾十年來降水增幅大結(jié)論一致。但顯著性趨勢(shì)最強(qiáng)的是ITA法，其次是M-K法，Sen’s斜率估計(jì)表現(xiàn)較為一般。ITA法分析了不同月份不同級(jí)別降水量趨勢(shì)，能夠?qū)σ恍┰路莸碾[藏趨勢(shì)做出回答，比較系統(tǒng)全面分析出了各個(gè)階段的變化趨勢(shì)；M-K與ITA趨勢(shì)變化大體一致，3個(gè)站點(diǎn)的多年降水量均呈現(xiàn)出增加趨勢(shì)，尤其是顯著月份與ITA高度擬合；Sen’s斜率估計(jì)與兩種方法得出結(jié)論有一定差異性，線性擬合系數(shù)低，原因可能是Sen’s斜率估計(jì)法是在眾多斜率中尋找中值，而中值斜率的正負(fù)不能很好地評(píng)估時(shí)間序列的增減趨勢(shì)，而且Sen’s斜率估計(jì)法不能很好地排除異常值的影響，但分析得出整體呈緩慢增加趨勢(shì)。

(2)季節(jié)特征上，夏季(6—8月)降水增強(qiáng)幅度最大，6月和8月最為顯著,與王寶靈(1997)分析1961—1990中國(guó)西北地區(qū)30年間6月降水量遞增的結(jié)果吻合；冬季(12月、次年1月和2月)為其次，1月增強(qiáng)幅度最小；春季(3—5月)增強(qiáng)幅度最小，一年中夏季降水量貢獻(xiàn)最大，春季最小。原因可能是冰川加速消融(IPCC，2007)導(dǎo)致的徑流增加，潛在植被變化(Jiang，2008)導(dǎo)致的土壤水分平衡過程變化等。

(3)ITA的優(yōu)勢(shì)在于，它能夠根據(jù)圖形非常直觀地觀察出降水量在不同階段的變化趨勢(shì)，彌補(bǔ)了M-K和Sen’s斜率估計(jì)法只能檢驗(yàn)單調(diào)變化趨勢(shì)的不足，不需要假設(shè)，保證了對(duì)趨勢(shì)存在性的直觀檢查可能性，能夠確定最小值隱藏的趨勢(shì)，彌補(bǔ)了M-K 和Sen’s斜率估計(jì)法對(duì)短時(shí)間序列分析的低精確性。