·，

(新疆伊犁水文勘測局，新疆 伊犁 835000)

20世紀(jì)以來的全球變暖趨勢已成為公認(rèn)的事實[1]，在全球氣溫趨于變暖的大背景下，不同地區(qū)的氣候也發(fā)生了不同程度的變化，引起了不同國家、學(xué)者的關(guān)注，為更好的應(yīng)對氣候變化所引起的一系列問題，投入了大量的精力與財力，尋求較好的解決途徑與應(yīng)對措施。由于水資源對全球變暖趨勢有著敏感的響應(yīng)，并且水資源的合理配置制約著社會經(jīng)濟發(fā)展，因此目前不同學(xué)者關(guān)于全球氣候變暖對水資源的響應(yīng)方面展開了大量的研究，取得一定成果。伴隨著全球氣候的變暖、人類的強烈活動，極端水文事件的頻發(fā)，給人類的發(fā)展帶來了巨大的損失，給當(dāng)?shù)鼐用竦纳a(chǎn)和生活帶來的一定的危害，這就急需去重新認(rèn)識人為及外界因素干擾下的水文循環(huán)規(guī)律，尤其是像新疆這種大多數(shù)據(jù)河流發(fā)源于高山冰川積雪消融區(qū)的區(qū)域。因此，開展該區(qū)域的水文規(guī)律分析有助于所在區(qū)域加強冰雪融水的開發(fā)利用、水文過程的再現(xiàn)以及水資源合理開發(fā)利用等。

水文規(guī)律分析主要分為序列的趨勢性分析、序列的突變性分析以及水文序列的周期性變化規(guī)律分析的三個方面，其中對序列規(guī)律分析中被大量采用的方法有M-K法、小波分析方法、最小二乘法等。總而言之，由于地理位置和氣候的差異，我國各地氣候長期變化具有不同的特征。目前，新疆地區(qū)與內(nèi)地有明顯的氣候差異，新疆地區(qū)主要河流易受氣候變化及人類活動的影響，但有關(guān)學(xué)者對新疆南北疆地區(qū)的氣候變化規(guī)律也有不少研究[8-10]，主要涉及氣候變化特征研究等內(nèi)容。如趙澤珍等[11]于2012年對新疆伊犁河流域的水文序列變化規(guī)律進行研究，其中包括流域的氣溫序列、降水序列以及徑流序列的變化趨勢。葉佰生[12]于1997年利用近三十年氣象資料對天山山區(qū)伊犁河流域中低山區(qū)的五個典型水文站點詳細(xì)分析了降水和氣溫的變化規(guī)律，總結(jié)出該地區(qū)中低山帶的降水和氣溫與其對應(yīng)高度的關(guān)系夏季較好,冬季較差的結(jié)論。馬建勇[13]于2012年對北疆地區(qū)近半個世紀(jì)以來的20個氣象站點的氣溫、降水要素進行了全面詳細(xì)的分析。

由于近年來受人類活動、氣候變化以及水資源開發(fā)利用等因素的影響，使不同區(qū)域所在河流發(fā)生的極端水文災(zāi)害事件的頻次增加、強度增強，給所在區(qū)域帶來了極大自然災(zāi)害損失，如研究區(qū)所在的鞏乃斯河流域。因此，為了更好的認(rèn)識所在區(qū)域的產(chǎn)匯流規(guī)律、指導(dǎo)所在流域雨、洪資源的合理開發(fā)利用、自然災(zāi)害的預(yù)估、預(yù)警與應(yīng)對等，對所在研究區(qū)的水文氣象序列的變化特征以及趨勢分析具有重要研究意義。

1 研究區(qū)概況及資料來源

鞏乃斯河為伊犁河上的第三大支流，河流起源于艾肯達坂區(qū)，從東向西橫穿鞏乃斯河谷，最后匯入伊犁河。流域全長280 km，集水總面積為7 707 km2；流域平均高程2 940 m，自西向東由800 m遞增至4 000 m以上,使西來的水汽逐漸抬升促使降水的形成，自西向東多年平均年降水量為3 600～850 mm。流域控制斷面處為則克臺水文站，水文站處于東經(jīng)83°16′、北緯43°31′處，海拔為840 m。

本文選取新疆鞏乃斯河流域則克臺水文站1972—2015年的年平均氣溫、降水資料系列進行趨勢分析，數(shù)據(jù)來源于新疆水文年鑒。

2 研究方法

本文借助線性趨勢分析法[2]、滑動平均法[3]對流域的氣溫、降水時間序列進行趨勢性變化分析；借助M-K法[4]對流域的氣溫、降水時間序列的突變情況進行分析；借助小波變換分析法[5]揭示了流域的氣溫、降水時間序列的周期性演變規(guī)律。

3 分析結(jié)果

3.1 年際變化趨勢分析

根據(jù)對研究區(qū)年氣溫序列的統(tǒng)計分析可知，鞏乃斯河流域多年氣溫的平均值為8.4℃，其中2006年出現(xiàn)最高氣溫10.1℃，1972年出現(xiàn)最低氣溫5.6℃。從氣溫序列的變化趨勢圖(圖1)中可知，研究區(qū)氣溫序列在1972～2015年持續(xù)增長，年平均氣溫變化傾向率為0.65℃/10a。由氣溫的5年滑動平均曲線可知，其變化可明顯分為上升-下降-上升-下降四個階段，1983年之前序列呈上升趨勢，1984～1988年之間序列呈下降趨勢；1989～2010年之間氣溫變化呈持續(xù)上升趨勢，且從1994年開始?xì)鉁孛黠@高于多年均值，但從2010年以后氣溫變化呈下降趨勢。

根據(jù)研究區(qū)年降水序列統(tǒng)計分析可知，流域多年平均降水量為514.9 mm，主要集中于3～8月份，占全年的63.5%；1998年為特豐水年，雨量為913.0 mm，1995年為特枯水年，雨量為316.4 mm。對降水年際變化進行線性分析如圖1所示。由圖1可以看出，鞏乃斯河流域年降水量波動幅度較大，但總體上降水量逐年增加，1972～2012年總的降水傾向率值為17.2 mm/10a。由5年滑動平均降水曲線分析得出，流域降水變化趨勢可分為兩個時段， 1972～2002年降水量明顯波動增加，其降水傾向率的值為55.8 mm/10a，這一階段降水量較為豐富。2002年后降水量下降，并以2003年為臨界時，2003年以后降水量明顯減少，降水傾向率的值為-75.8 mm/10a。

圖1 鞏乃斯河流域 1972～2015 年

3.2 Mann-Kendall突變分析

M-K法是一種非參數(shù)統(tǒng)計檢驗方法，判定其突變的參數(shù)為α=0.05、U0.05=±1.96。如果|UFm|>U0.05，則代表序列有顯著的變化。當(dāng)UF或UB大于零，代表有上升的變化，反之則有下降的變化。若兩值相交于參數(shù)的區(qū)域內(nèi)，相交處則表示為序列的突變點。本研究主要借助M-K法對鞏乃斯河流域的年氣溫、年降水序列進行突變特性分析，確定序列相應(yīng)的突變時間。

從M-K檢驗曲線(圖2)中可以看出：近44年來，年均氣溫具有突變現(xiàn)象，氣溫從較冷的階段向較暖階段進行越變；從UF曲線可見，1970～1980年氣溫處于持續(xù)上升時期，1984至1989年呈氣溫逐年降低，1991年又開始出現(xiàn)較顯著的持續(xù)上升趨勢；同時在置信區(qū)間內(nèi)UF曲線與UB線于1991年左右相交，故認(rèn)為1991年為突變年份。年平均氣溫在1998年以后上升超過α=0.05的置信水平線，表明鞏乃斯河流域平均氣溫從1998年開始上升趨勢是較為明顯的，這與全球性的氣候變暖是相一致的。

圖2 鞏乃斯河流域年均氣溫及

年降水量變化情況由圖2(b)可知，UF線與UB線都在α=0.05的信度下，UF線在信度的上下限之間，自1979年以前UF值<0，說明從1972～1979年之間降水序列呈減少-增加兩個變化趨勢；1980～1997年之間降水序列呈波動性減少趨勢，自1998～2010年之間UF值>0，降水呈波動增長趨勢，但不顯著；2011年以后UF值>0，但降水呈減少趨勢。UF和UB兩條曲線的交點出現(xiàn)在1979年、1982年和2012年處，說明降水序列在這三年有突變情況。

3.3 研究區(qū)氣溫及降水的主周期特性分析

3.3.1 氣溫序列的小波變換分析及主周期分析

將則克臺水文站年氣溫序列通過小波變換后，同時基于suffer8.0軟件做系數(shù)的等值線圖(圖3、圖4)。從圖3中分析得到氣溫序列在17～25年和26～32年的時間尺度上有周期性演變規(guī)律， 26～32年的尺度上周期震蕩較為明顯，從1972年至今表現(xiàn)為高低交替3個中心。在17～25年的尺度上也發(fā)生了3次低、高的周期變換?；跉鉁啬７降戎稻€圖(圖4)，得到最大的模方值出現(xiàn)在24～32年的時間尺度上，表示該尺度的規(guī)律最顯著。

圖3 氣溫序列小波系數(shù)實部等值線圖

圖4 氣溫小波系數(shù)模方等值線圖

從氣溫小波方差圖(圖5)中發(fā)現(xiàn)，鞏乃斯河流域年均氣溫方差的第一峰值是22年，代表該時間尺度下氣溫的震動周期最強，是第一主周期；小波方差值的第二峰值是29年，為年氣溫序列變化的第二主周期。

基于則克臺氣溫序列小波系數(shù)方差的檢驗(圖5)得知第一主周期(最大峰值)出現(xiàn)在22年時間尺度上，說明在22年前后的周期變換規(guī)律最為突出，可以代表氣溫序列的周期變換趨勢?；跉鉁匦蛄兄髦芷趫D(圖6)，在22年尺度上，有2次偏冷期到偏暖期的變換，平均變換周期在15年左右；第二主周期(第二峰值)出現(xiàn)在29年時間尺度上，也有2次偏冷期到偏暖期的變換，平均變換周期在19年左右。

圖5 氣溫變化的小波方差圖

圖6 溫度序列22年和29年時間尺度周期趨勢圖

3.3.2 年降水的小波變換分析及主周期特性分析

基于則克臺站44年年降水序列通過小波變換系數(shù)后得到相應(yīng)的等值線圖(圖7)。可看出，該序列在于5～13年，15～22年和23～32年的時間尺度上，存在演變規(guī)律。自20世紀(jì)80年代中期開始，降水在5～13年的時間尺度上完成了6次枯水-豐水交替過程；在15～22年的尺度上完成了了枯水-豐水交替的3次循環(huán)，在23～32年的尺度上完成2次降水枯水-豐水交替過程。

圖7 降水小波系數(shù)實部等值線圖

從年降水小波系數(shù)模方等值線圖(圖8)，可得到6～11年時間尺度有最大的模方值，表示該時間尺度下降水周期變化規(guī)律最顯著。從研究區(qū)年降水小波方差圖9中可以看出，鞏乃斯河流域年降水的方差的第一峰值是27年，第二峰值是8年，第三峰值是17年，說明年降水量在這三年的時間尺度上周期振蕩最強，分別為年降水量序列的第一、二和三主周期。

基于則克臺降水序列小波變換系數(shù)方差的檢驗(圖9)得知第一主周期(最大峰值)出現(xiàn)在27年時間尺度上，說明27年左右的周期變化規(guī)律最為顯著，第二至第三峰值分別出現(xiàn)在8年和17年時間尺度上，依次為第二、三主周期。這3個周期的變化規(guī)律可以表示降雨序列的周期變化趨勢?；诮邓蛄兄髦芷谮厔輬D(圖10)，可以分析出第一主周期和第二主周期分別有2次和6次降水偏少期到偏多期的變換過程，平均周期在18年和7年左右。

圖8 降水小波系數(shù)模方等值線圖

圖9 降水變化的小波方差圖

圖10 降水序列27年和8年時間尺度周期趨勢圖

4 結(jié)語

為科學(xué)合理的指導(dǎo)研究區(qū)的水資源合理開發(fā)利用，本文借助了氣候傾向率、滑動平均法、M-K法以及小波分析法對鞏乃斯河流域近44年的平均氣溫及降水變化特征進行了分析，可以看出：

(1)鞏乃斯河流域近44年氣溫波動中有明顯上升，平均每10年增加0.65℃，高于全球氣溫增長率。氣溫冷暖變化交替出現(xiàn)，在1989～2010年之間氣溫變化呈持續(xù)上升趨勢，且從1994年開始?xì)鉁孛黠@高于多年均值，但從2010年以后氣溫變化明顯下降。研究區(qū)年降水量波動變化比較大，但總體表現(xiàn)為增加，其降水傾向率是17.2 mm/10a。降水量在1972～2002年之間波動呈明顯增加，其降水傾向率是55.8 mm/10a，這一階段降水量較為豐富。進入21世紀(jì)后，降水減少，并以2003年為臨界時，2003年以后降水量呈明顯減少，降水傾向率為-75.8 mm/10a。

(2)鞏乃斯河流域氣溫序列的突變時間在1991年左右，于1998年以后氣溫增暖趨勢超過了α=0.05的置信水平線，上升趨勢十分顯著；研究區(qū)年降水量UF線在置信水平上下限之間，說明降水減少-偏多變化趨勢交替出現(xiàn)，但不顯著；年降水量在1979、1982、2012年存在著突變。

(3)研究區(qū)氣溫和降水序列在近44年的變化趨勢中包含了多個不同時間尺度的周期變化特性。年平均氣溫存在著明顯的17～25年和26～32年的時間尺度上的周期震蕩。從其周期變換的情況中可以推斷，鞏乃斯河流域的年氣溫在未來幾年內(nèi)仍會偏高，屬于暖期。年平均降水量在5～13年，15～22年及23～32年有較為明顯的周期性變化規(guī)律。15～22年及23～32年較5～13年振動周期表現(xiàn)較弱，在5～13年時間尺度上的周期振蕩特別明顯，其間年降水經(jīng)歷了6個少-多-少的循環(huán)過程。小波系數(shù)方差結(jié)果揭示了27年在1972～2015年整個時間域內(nèi)，周期振動最強，為降水序列的第一主周期。

以上表明：全球氣候的變化影響著研究區(qū)的自然水文循環(huán)規(guī)律，但由于部分流域水文監(jiān)測站點分布較稀疏，設(shè)備較落后，為了更好的指導(dǎo)所在流域雨、洪資源的合理開發(fā)利用及應(yīng)對，因此要合理布置氣象數(shù)據(jù)監(jiān)測站點，提高分析手段，以便為所在流域水文氣象要素變化特征及趨勢預(yù)測進行全面分析。